Hi everyone,

I am try to run geometry relaxation for a CoS slab. However, forces
keep oscillating around 0.5 eV/Ang and never reach convergence.
I use values of MD.MaxCGDispl  between 0.1 and 0.001 Ang.
This is the setup of the calculations:

xc.functional         GGA
xc.authors            PBE
PAO.BasisSize       DZP
PAO.EnergyShift     0.005 Ry
SpinPolarized        .false.
MeshCutoff            750 Ry
MaxSCFIterations    250
DM.MixingWeight      0.01
DM.NumberPulay       8
DM.Tolerance         1.d-4
SolutionMethod       diagon
ElectronicTemperature  40 meV
MD.TypeOfRun         cg
MD.NumCGsteps        100
MD.MaxCGDispl         0.01  Ang
MD.MaxForceTol        0.04 eV/Ang
DM.NumberKick   100

Any suggestion is very appreciate.

Thanks,

Marco

Responder a