<div dir="ltr"><div>By the way: if I remember correctly,  CP and PW do not give the same energy for charged cells. The difference is given by the net charge times the quantity 'Delta V(G=0)' printed in the output of CP</div><div><br></div><div>Paolo<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Mar 12, 2021 at 10:23 AM Paolo Giannozzi <<a href="mailto:p.giannozzi@gmail.com">p.giannozzi@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Correct, but scf convergence in CP is not THAT slow. The default value of the time step is dt=1, if I remember correctly. A value around 10 should yield a much faster convergence.</div><div><br></div><div>Paolo<br></div><div><br></div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Mar 12, 2021 at 10:18 AM Giuseppe Mattioli <<a href="mailto:giuseppe.mattioli@ism.cnr.it" target="_blank">giuseppe.mattioli@ism.cnr.it</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br>
Dear Chris<br>
scf convergence with cp.x actually used to be painfully slow! For this  <br>
reason, pw.x and cp.x should be interoperable, that is, you should be  <br>
able to restart a converged pw.x calculation with cp.x. Have a look  <br>
into the CPV/Restart_example/ directory, where restart details are  <br>
explained.<br>
HTH<br>
Giuseppe<br>
<br>
Quoting Christoph Wolf <wolf.christoph@qns.science>:<br>
<br>
> Dear all,<br>
><br>
> I am trying to run an SCF calculation in cp.x (because I want to use the<br>
> spin contamination analysis of cp.x), however reaching convergence with<br>
> cp.x seems painfully slow - probably my fault as I am not familiar with<br>
> cp.x.<br>
><br>
> System is a spin polarized calculation for a molecule in a vacuum box. The<br>
> SCF calculation converges in 30 steps:<br>
><br>
> !    total energy              =    -813.86536682 Ry<br>
>      Harris-Foulkes estimate   =    -813.86536682 Ry<br>
>      estimated scf accuracy    <          1.3E-11 Ry<br>
><br>
> But the cp.x calculation is still running after 2000 steps<br>
><br>
> 1997    0.000605472420187    0.0    0.00    -406.607110121158<br>
>  1998    0.000604666115183    0.0    0.00    -406.607351342414<br>
> 1999    0.000603861731893    0.0    0.00    -406.607592243454<br>
> ...<br>
><br>
> (NB: the energy difference between CP and PW at this point is ~0.650182333<br>
> Ry)<br>
><br>
> I am running the cp.x calculation with<br>
><br>
> calculation = 'scf'<br>
> and electron_dynamics = 'damp<br>
><br>
> Is that the right way to do this? If yes, am I missing some critical<br>
> parameters for improving convergence?<br>
><br>
> I attached my input below - any help is much appreciated!<br>
><br>
> Best,<br>
> Chris<br>
><br>
> &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>
> ...<br>
> ..<br>
><br>
> --<br>
> IBS Center for Quantum Nanoscience<br>
> Seoul, South Korea<br>
<br>
<br>
<br>
GIUSEPPE MATTIOLI<br>
CNR - ISTITUTO DI STRUTTURA DELLA MATERIA<br>
Via Salaria Km 29,300 - C.P. 10<br>
I-00015 - Monterotondo Scalo (RM)<br>
Mob (*preferred*) +39 373 7305625<br>
Tel + 39 06 90672342 - Fax +39 06 90672316<br>
E-mail: <<a href="mailto:giuseppe.mattioli@ism.cnr.it" target="_blank">giuseppe.mattioli@ism.cnr.it</a>><br>
<br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">www.max-centre.eu</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
</blockquote></div><br clear="all"><br>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 206, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>
</blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 206, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>