<div dir="ltr">Dear all,<div><br></div><div>I am trying to run an SCF calculation in cp.x (because I want to use the spin contamination analysis of cp.x), however reaching convergence with cp.x seems painfully slow - probably my fault as I am not familiar with cp.x.</div><div><br></div><div><div>System is a spin polarized calculation for a molecule in a vacuum box. The SCF calculation converges in 30 steps:</div><div><br></div><div>!    total energy              =    -813.86536682 Ry<br>     Harris-Foulkes estimate   =    -813.86536682 Ry<br>     estimated scf accuracy    <          1.3E-11 Ry<br></div><div><br></div><div>But the cp.x calculation is still running after 2000 steps </div><div><br></div><div>1997    0.000605472420187    0.0    0.00    -406.607110121158<br></div><div> 1998    0.000604666115183    0.0    0.00    -406.607351342414<br></div><div>1999    0.000603861731893    0.0    0.00    -406.607592243454<br></div><div>...</div><div><br></div><div>(NB: the energy difference between CP and PW at this point is ~0.650182333 Ry)</div><div><br></div><div>I am running the cp.x calculation with </div><div><br></div><div>calculation = 'scf'<br></div><div>and electron_dynamics = 'damp</div><div><br></div><div>Is that the right way to do this? If yes, am I missing some critical parameters for improving convergence?</div><div><br></div><div>I attached my input below - any help is much appreciated!</div><div><br></div><div>Best,</div><div>Chris </div><div><br></div><div>&CONTROL<br>  calculation = "scf" <br>  prefix = "fepc"<br>  pseudo_dir = "/home/chwolf/Work/DFT/00_pseudo/"<br>  outdir = "tmp"<br>  wf_collect = .true<br>  tprnfor = .true.<br>  restart_mode = "from_scratch"<br>  nstep      = 2000<br>/ <br>&SYSTEM<br>   ibrav = 0<br>   nat  = 57<br>   ntyp = 4<br>   nspin = 2<br>   occupations = "fixed"<br>   ecutwfc = 60 <br>   ecutrho = 240 <br>   tot_magnetization=1<br>   tot_charge=-1<br>   starting_magnetization(1) = 0.5<br>   starting_magnetization(2) = 0.0<br>   starting_magnetization(3) = 0.0<br>   starting_magnetization(4) = 0.0<br>   nr1b=24, nr2b=24, nr3b=12<br>/ <br>&ELECTRONS<br>  conv_thr=1e-10<br>  electron_maxstep = 350<br>  electron_dynamics = 'damp'<br>/ <br><br>CELL_PARAMETERS angstrom<br>29.47070000000000 0.00000000000000 0.00000000000000<br>0.00000000000000 29.47070000000000 0.00000000000000<br>0.00000000000000 0.00000000000000 15.000<br><br>ATOMIC_SPECIES<br>Fe    55.9330    Fe_ONCV_PBE_sr.upf<br>C     12.0107    C_ONCV_PBE_sr.upf<br>N     14.0067    N_ONCV_PBE-1.0.upf<br>H      1.0000    H_ONCV_PBE_sr.upf<br><br>ATOMIC_POSITIONS angstrom<br></div><div>...</div><div>..</div><div><br></div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">--<br>IBS Center for Quantum Nanoscience<br>Seoul, South Korea<blockquote type="cite" style="font-size:12.8px"><div dir="ltr"><div><div dir="ltr"></div></div></div></blockquote></div></div></div></div>