[CP2K-user] SCF does not converge for 1st step

Travis polla... at gmail.com
Fri Oct 11 16:08:38 UTC 2019


Hi,

&SCF
 ...
    &SMEAR  T
      METHOD  FERMI_DIRAC
      ELECTRONIC_TEMPERATURE  300
    &END SMEAR
 ...
&END SCF

-T

On Friday, October 11, 2019 at 12:58:02 PM UTC-3, Jiapeng Liu wrote:
>
> For your reference, the basis set I used is DZVP-MOLOPT-SR-GTH and 
> potential is GTH-PBE. Thanks for any help.
>
> On Friday, October 11, 2019 at 11:54:39 PM UTC+8, Jiapeng Liu wrote:
>>
>> Dear CP2K users,
>>
>> I am just learning CP2k now and I wanna use CP2K to optimize a Li solid 
>> state electrolyte. The problem is that the SCF run does not converge for 
>> the first step even after 180 scf steps. I am wondering if this is ok for 
>> the continue calculations or there is something wrong with my parameter 
>> setting. Can you kindly help me have a look. The input file is like 
>> following
>>
>> &GLOBAL
>>   PROJECT_NAME relax
>>   RUN_TYPE GEO_OPT
>>   PRINT_LEVEL LOW
>> &END GLOBAL
>>
>> &MOTION
>>   &GEO_OPT
>>     MAX_ITER  200
>>     OPTIMIZER CG
>>     MAX_DR    1.0E-3
>>     MAX_FORCE 4.5E-4
>>     RMS_DR    1.0E-3
>>     RMS_FORCE 3.0E-4
>>   &END
>> &END MOTION
>>
>> &FORCE_EVAL
>>   METHOD QS
>>   STRESS_TENSOR ANALYTICAL
>>   &DFT
>>     CHARGE 0
>>     MULTIPLICITY 1
>>     BASIS_SET_FILE_NAME  BASIS_MOLOPT
>>     POTENTIAL_FILE_NAME  GTH_POTENTIALS
>>
>>     &MGRID
>>       NGRIDS 5
>>       CUTOFF 600
>>       REL_CUTOFF 80
>>     &END MGRID
>>
>>     &QS
>>       METHOD GPW ! to optimize the geometry the GPW method will be used
>>       EPS_DEFAULT 1.0E-12
>>       EXTRAPOLATION USE_GUESS
>>     &END QS
>>
>>     &SCF
>>       MAX_SCF 200
>>       EPS_SCF 1.0E-6
>>       SCF_GUESS RESTART
>>       ADDED_MOS  200
>>       &MIXING T
>>         METHOD BROYDEN_MIXING
>>         ALPHA   0.4
>>         NMIXING  5
>>         NBUFFER  8
>>       &END MIXING
>>     &END SCF
>>
>>     &XC
>>       FUNCTIONAL_ROUTINE NEW
>>       DENSITY_CUTOFF  1.0E-12
>>       GRADIENT_CUTOFF 1.0E-12
>>       TAU_CUTOFF      1.0E-12
>>       &XC_FUNCTIONAL
>>         &PBE  ! PBE exchange-correlation functional
>>           PARAMETRIZATION ORIG
>>         &END PBE
>>       &END XC_FUNCTIONAL
>>    &END XC
>>    &POISSON
>>      POISSON_SOLVER  PERIODIC
>>      PERIODIC  XYZ
>>    &END POISSON
>>    &KPOINTS
>>      SCHEME  GAMMA  2  2  2
>>      FULL_GRID  .TRUE.
>>    &END KPOINTS
>>   &END DFT
>>
>>   &SUBSYS
>>     &CELL
>>       ABC .....
>>       ALPHA_BETA_GAMMA  ..........
>>       PERIODIC XYZ
>>       MULTIPLE_UNIT_CELL  .....
>>     &END CELL
>>     &COORD
>>       ...........
>>     &END COORD
>>
>>     &KIND Li
>>     ............
>>     &END KIND
>>   &END SUBSYS
>>   &PRINT
>>     &FORCES ON
>>     &END FORCES
>>   &END PRINT
>> &END FORCE_EVAL
>>
>> and the first scf run is shown like this
>>
>>  SCF WAVEFUNCTION OPTIMIZATION
>>
>>   Step     Update method      Time    Convergence         Total energy    
>> Change
>>   
>> ------------------------------------------------------------------------------
>>      1 NoMix/Diag. 0.40E+00    2.2    33.06169264     -1105.3440392037 
>> -1.11E+03
>>      2 Broy./Diag. 0.40E+00    2.5    18.58674968     -1325.2412363377 
>> -2.20E+02
>>      3 Broy./Diag. 0.40E+00    2.5    13.64590095     -1613.5521238999 
>> -2.88E+02
>>      4 Broy./Diag. 0.40E+00    2.5     6.33206879     -1767.5043485819 
>> -1.54E+02
>>      5 Broy./Diag. 0.40E+00    2.5     3.10815819     -1788.5141959477 
>> -2.10E+01
>>      6 Broy./Diag. 0.40E+00    2.5     3.46033998     -1840.7455766805 
>> -5.22E+01
>>      7 Broy./Diag. 0.40E+00    2.5     2.74195755     -1831.3589058361  
>> 9.39E+00
>>    ......................
>>    196 Broy./Diag. 0.40E+00    2.7     1.63271069     -1923.1236698860  
>> 6.48E+00
>>    197 Broy./Diag. 0.40E+00    2.7     1.72647036     -1924.2660229966 
>> -1.14E+00
>>    198 Broy./Diag. 0.40E+00    2.6     1.75474955     -1920.5269108168  
>> 3.74E+00
>>    199 Broy./Diag. 0.40E+00    2.6     1.57890425     -1926.1904852274 
>> -5.66E+00
>>    200 Broy./Diag. 0.40E+00    2.6     1.67473633     -1926.3874701171 
>> -1.97E-01
>>
>> Thanks very much.
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20191011/95308b95/attachment.htm>


More information about the CP2K-user mailing list