[QE-users] Convergence in SCF calculation

mabola safar masrfa at googlemail.com
Wed Dec 8 22:12:43 CET 2021


Dear All,

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.

This is the input file

3sp
 &control
    calculation = 'scf'
    restart_mode = 'from_scratch'
    prefix='3sp_MnBi2Te4',
    pseudo_dir = '/home/qe',
    outdir='/home/qe',
    wf_collect = .false.,
    tprnfor=.true.,
    forc_conv_thr=1.0D-5
 /
 &system
    ibrav = 0, nat = 21, ntyp = 5,
    lspinorb = .true.,
    noncolin = .true.,
    ecutwfc = 100.0,
    ecutrho = 200.0,
    nspin = 1,
    occupations = 'smearing',
    smearing = 'mv',
    degauss = 0.02,
    nosym = .false.,
    starting_magnetization(1) = 0.9,
    starting_magnetization(2) = -0.9,
    starting_magnetization(3) = 0.9,
    starting_magnetization(4) = 0.0,
    starting_magnetization(5) = 0.0,
!    constrained_magnetization = 'atomic',
    angle1(1) = 0.0,
    angle1(2) = 0.0,
    angle1(3) = 0.0,
    angle1(4) = 0.0,
    angle1(5) = 0.0,
!    lambda = 5.0
 /
 &electrons
    mixing_beta=0.2,
    electron_maxstep = 400,
    conv_thr = 1.0D-9
 /
 &ions
 /
 &cell
    cell_dynamics='none'
 /
CELL_PARAMETERS {angstrom}
  4.3338000000  0.0000000000  0.0000000000
 -2.1669000000  3.75318087852 0.0000000000
  0.0000000000  0.0000000000 59.1826680000
ATOMIC_SPECIES
Mn1 54.938 Mn.rel-pbe-spn-kjpaw_psl.0.3.1.UPF
Mn2 54.938 Mn.rel-pbe-spn-kjpaw_psl.0.3.1.UPF
Mn3 54.938 Mn.rel-pbe-spn-kjpaw_psl.0.3.1.UPF
Bi 208.98 Bi.rel-pbe-dn-kjpaw_psl.1.0.0.UPF
Te 127.6 Te.rel-pbe-dn-kjpaw_psl.0.2.2.UPF
ATOMIC_POSITIONS {crystal}
Mn1  0.333333333  0.666666667  0.322583180 1 1 1
Mn2  0.000000000  0.000000000  0.552999740 1 1 1
Mn3  0.666666667  0.333333333  0.783416310 1 1 1
Bi  0.666666667  0.333333333  0.385864790 1 1 1
Bi  0.000000000  0.000000000  0.259301580 1 1 1
Bi  0.666666667  0.333333333  0.489718140 1 1 1
Bi  0.000000000  0.000000000  0.846697900 1 1 1
Bi  0.333333333  0.666666667  0.616281340 1 1 1
Bi  0.333333333  0.666666667  0.720134690 1 1 1
Te  0.666666667  0.333333333  0.295642880 1 1 1
Te  0.333333333  0.666666667  0.414749810 1 1 1
Te  0.000000000  0.000000000  0.349523470 1 1 1
Te  0.000000000  0.000000000  0.460833110 1 1 1
Te  0.333333333  0.666666667  0.230416550 1 1 1
Te  0.000000000  0.000000000  0.756476010 1 1 1
Te  0.000000000  0.000000000  0.645166370 1 1 1
Te  0.666666667  0.333333333  0.875582930 1 1 1
Te  0.333333333  0.666666667  0.526059450 1 1 1
Te  0.666666667  0.333333333  0.579940040 1 1 1
Te  0.666666667  0.333333333  0.691249680 1 1 1
Te  0.333333333  0.666666667  0.810356600 1 1 1
K_POINTS {AUTOMATIC}
11 11 1 0 0 0

This is the output file

............
     iteration #333     ecut=   100.00 Ry     beta= 0.20
     Davidson diagonalization with overlap
     ethr =  1.58E-10,  avg # of iterations =  4.0

     negative rho (up, down):  5.333E-03 1.542E-03

     total cpu time spent up to now is    86026.7 secs

     total energy              =  -12505.34455324 Ry
     estimated scf accuracy    <       0.00000122 Ry

     total magnetization       =     0.00     0.00     4.94 Bohr mag/cell
     absolute magnetization    =    14.89 Bohr mag/cell

     iteration #334     ecut=   100.00 Ry     beta= 0.20
     Davidson diagonalization with overlap
     ethr =  1.58E-10,  avg # of iterations =  4.2

     negative rho (up, down):  5.333E-03 1.523E-03

     total cpu time spent up to now is    86331.7 secs

     total energy              =  -12505.34451496 Ry
     estimated scf accuracy    <       0.00000066 Ry

     total magnetization       =     0.00    -0.00     4.94 Bohr mag/cell
     absolute magnetization    =    14.89 Bohr mag/cell

     iteration #335     ecut=   100.00 Ry     beta= 0.20
     Davidson diagonalization with overlap

Would anyone please provide any suggestions on how to make the SCF
calculation converge?

Thank you very much in advance.

Kind regards,

Kieran
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20211208/ff4fa592/attachment.html>


More information about the users mailing list