<div dir="ltr">Dear All,<div><br></div><div>I paste my input file for the SCF calculation below. The output file shows that the iterations reaches 335 steps but the calculator still did not converge.</div><div><br></div><div>This is the input file</div><div><br></div><div>3sp<br> &control<br>    calculation = 'scf'<br>    restart_mode = 'from_scratch'<br>    prefix='3sp_MnBi2Te4',<br>    pseudo_dir = '/home/qe',<br>    outdir='/home/qe',<br>    wf_collect = .false.,<br>    tprnfor=.true.,<br>    forc_conv_thr=1.0D-5<br> /<br> &system<br>    ibrav = 0, nat = 21, ntyp = 5,<br>    lspinorb = .true.,<br>    noncolin = .true.,<br>    ecutwfc = 100.0,<br>    ecutrho = 200.0,<br>    nspin = 1,<br>    occupations = 'smearing',<br>    smearing = 'mv',<br>    degauss = 0.02,<br>    nosym = .false.,<br>    starting_magnetization(1) = 0.9,<br>    starting_magnetization(2) = -0.9,<br>    starting_magnetization(3) = 0.9,<br>    starting_magnetization(4) = 0.0,<br>    starting_magnetization(5) = 0.0,<br>!    constrained_magnetization = 'atomic',<br>    angle1(1) = 0.0,<br>    angle1(2) = 0.0,<br>    angle1(3) = 0.0,<br>    angle1(4) = 0.0,<br>    angle1(5) = 0.0,<br>!    lambda = 5.0<br> /<br> &electrons<br>    mixing_beta=0.2,<br>    electron_maxstep = 400,<br>    conv_thr = 1.0D-9<br> /<br> &ions<br> /<br> &cell<br>    cell_dynamics='none'<br> /<br>CELL_PARAMETERS {angstrom}<br>  4.3338000000  0.0000000000  0.0000000000<br> -2.1669000000  3.75318087852 0.0000000000<br>  0.0000000000  0.0000000000 59.1826680000<br>ATOMIC_SPECIES<br>Mn1 54.938 Mn.rel-pbe-spn-kjpaw_psl.0.3.1.UPF<br>Mn2 54.938 Mn.rel-pbe-spn-kjpaw_psl.0.3.1.UPF<br>Mn3 54.938 Mn.rel-pbe-spn-kjpaw_psl.0.3.1.UPF<br>Bi 208.98 Bi.rel-pbe-dn-kjpaw_psl.1.0.0.UPF<br>Te 127.6 Te.rel-pbe-dn-kjpaw_psl.0.2.2.UPF<br>ATOMIC_POSITIONS {crystal}<br>Mn1  0.333333333  0.666666667  0.322583180 1 1 1<br>Mn2  0.000000000  0.000000000  0.552999740 1 1 1<br>Mn3  0.666666667  0.333333333  0.783416310 1 1 1<br>Bi  0.666666667  0.333333333  0.385864790 1 1 1<br>Bi  0.000000000  0.000000000  0.259301580 1 1 1<br>Bi  0.666666667  0.333333333  0.489718140 1 1 1<br>Bi  0.000000000  0.000000000  0.846697900 1 1 1<br>Bi  0.333333333  0.666666667  0.616281340 1 1 1<br>Bi  0.333333333  0.666666667  0.720134690 1 1 1<br>Te  0.666666667  0.333333333  0.295642880 1 1 1<br>Te  0.333333333  0.666666667  0.414749810 1 1 1<br>Te  0.000000000  0.000000000  0.349523470 1 1 1<br>Te  0.000000000  0.000000000  0.460833110 1 1 1<br>Te  0.333333333  0.666666667  0.230416550 1 1 1<br>Te  0.000000000  0.000000000  0.756476010 1 1 1<br>Te  0.000000000  0.000000000  0.645166370 1 1 1<br>Te  0.666666667  0.333333333  0.875582930 1 1 1<br>Te  0.333333333  0.666666667  0.526059450 1 1 1<br>Te  0.666666667  0.333333333  0.579940040 1 1 1<br>Te  0.666666667  0.333333333  0.691249680 1 1 1<br>Te  0.333333333  0.666666667  0.810356600 1 1 1<br>K_POINTS {AUTOMATIC}<br>11 11 1 0 0 0<br></div><div><br></div><div>This is the output file</div><div><br></div><div>............</div><div>     iteration #333     ecut=   100.00 Ry     beta= 0.20<br>     Davidson diagonalization with overlap<br>     ethr =  1.58E-10,  avg # of iterations =  4.0<br><br>     negative rho (up, down):  5.333E-03 1.542E-03<br><br>     total cpu time spent up to now is    86026.7 secs<br><br>     total energy              =  -12505.34455324 Ry<br>     estimated scf accuracy    <       0.00000122 Ry<br><br>     total magnetization       =     0.00     0.00     4.94 Bohr mag/cell<br>     absolute magnetization    =    14.89 Bohr mag/cell<br><br>     iteration #334     ecut=   100.00 Ry     beta= 0.20<br>     Davidson diagonalization with overlap<br>     ethr =  1.58E-10,  avg # of iterations =  4.2<br><br>     negative rho (up, down):  5.333E-03 1.523E-03<br><br>     total cpu time spent up to now is    86331.7 secs<br><br>     total energy              =  -12505.34451496 Ry<br>     estimated scf accuracy    <       0.00000066 Ry<br><br>     total magnetization       =     0.00    -0.00     4.94 Bohr mag/cell<br>     absolute magnetization    =    14.89 Bohr mag/cell<br><br>     iteration #335     ecut=   100.00 Ry     beta= 0.20<br>     Davidson diagonalization with overlap<br></div><div><br></div><div>Would anyone please provide any suggestions on how to make the SCF calculation converge?</div><div><br></div><div>Thank you very much in advance.</div><div><br></div><div>Kind regards,</div><div><br></div><div>Kieran</div></div>