<div dir="ltr"><div>Hi Everyone,</div><div><br></div><div>I'm currently having issues achieving convergence for my SCF calculations with SOC on nickel intercalated TaSe2. Since it is a metal, I've tried increasing the total number of bands and that didn't really work. I've also tried using local Thomas Fermi mixing due to address its elongated unit cell, as well as increased the ecutwfc and k-point density. So far the best scf accuracy I can achieve it ~2e-6 eV. One thing I've noticed is when the scf accuracy reaches 2e-6 eV, it will begin to oscillate from between 3e-6 to 2e-6 eV. I'm not sure what else to try. Anyways, thank you for reading this post and if you have any suggestions or tips, I would greatly appreciate hearing them.</div><div><br></div><div>Here is my input</div><div>#################################################<br></div><div>&CONTROL<br>    calculation='scf'<br>    prefix='TaNiSe'<br>    restart_mode='from_scratch'<br>    outdir='./outdir/'<br>    pseudo_dir='/global/scratch/lolzen/qe-6.7/SSSP_precision_pseudos/'<br>/<br>&SYSTEM<br>    ibrav = 0, celldm(1)=12.991867, <br>    nat = 26, ntyp = 3,<br>    ecutwfc = 70,<br>    ecutrho = 700,<br>    occupations = 'smearing', degauss = 0.01, smearing = 'mv'</div><div>    nbnd = 660<br></div><div>    noncolin=.true.<br>    lspinorb=.true.<br>    starting_magnetization(1)=0.8<br>/<br>&ELECTRONS<br>    mixing_mode='local-TF'<br>    conv_thr = 1.0d-9,<br>/<br>ATOMIC_SPECIES<br>  Ni    58.6934    Ni.rel-pbe-spn-kjpaw_psl.1.0.0.UPF<br>  Ta    180.94788  Ta.rel-pbe-spfn-kjpaw_psl.1.0.0.UPF<br>  Se    78.96      Se.rel-pbe-dn-kjpaw_psl.1.0.0.UPF<br><br>CELL_PARAMETERS (alat)<br>   1.004956653   0.000000000   0.000000000<br>  -0.502478327   0.870317991   0.000000000<br>   0.000000000   0.000000000   1.808363590<br><br>ATOMIC_POSITIONS (crystal)<br>Ni            0.0000000000       -0.0000000000        0.0000000000<br>Ni           -0.0000000000       -0.0000000000        0.5000000000<br>Ta            0.0000000000        0.0000000000        0.7500000000<br>Ta            0.0000000000       -0.0000000000        0.2500000000<br>Ta            0.9795988645        0.4897994323        0.2500000000<br>Ta            0.5102005677        0.4897994323        0.2500000000<br>Ta            0.5102005677        0.0204011355        0.2500000000<br>Ta            0.0204011355        0.5102005677        0.7500000000<br>Ta            0.4897994323        0.5102005677        0.7500000000<br>Ta            0.4897994323        0.9795988645        0.7500000000<br>Se            0.6666666700        0.3333333300        0.3922175119<br>Se            0.3333333300        0.6666666700        0.6077824881<br>Se            0.3333333300        0.6666666700        0.8922175119<br>Se            0.6666666700        0.3333333300        0.1077824881<br>Se            0.6624039120        0.8312019560        0.3820242878<br>Se            0.1687980440        0.8312019560        0.3820242878<br>Se            0.1687980440        0.3375960880        0.3820242878<br>Se            0.3375960880        0.1687980440        0.6179757122<br>Se            0.8312019560        0.1687980440        0.6179757122</div><div>Se            0.8312019560        0.6624039120        0.6179757122<br>Se            0.3375960880        0.1687980440        0.8820242878<br>Se            0.8312019560        0.1687980440        0.8820242878<br>Se            0.8312019560        0.6624039120        0.8820242878<br>Se            0.6624039120        0.8312019560        0.1179757122<br>Se            0.1687980440        0.8312019560        0.1179757122<br>Se            0.1687980440        0.3375960880        0.1179757122<br><br>  <br>K_POINTS automatic<br> 13 13 7 0 0 0</div><div>#################################################</div><div>End of input</div><div><br></div><div>Last few iterations of output<br></div><div>#################################################</div><div> Davidson diagonalization with overlap<br>     ethr =  5.54E-10,  avg # of iterations =  4.2<br><br>     total cpu time spent up to now is    62608.6 secs<br><br>     total energy              =  -17171.63213074 Ry<br>     estimated scf accuracy    <       0.00000313 Ry<br><br>     total magnetization       =     0.00     0.00     2.38 Bohr mag/cell<br>     absolute magnetization    =     2.89 Bohr mag/cell<br><br>     iteration # 34     ecut=    70.00 Ry     beta= 0.70<br>     Davidson diagonalization with overlap<br>     ethr =  5.54E-10,  avg # of iterations =  4.0<br><br>     total cpu time spent up to now is    64690.6 secs<br><br>     total energy              =  -17171.63214316 Ry<br>     estimated scf accuracy    <       0.00000314 Ry<br><br>     total magnetization       =    -0.00     0.00     2.38 Bohr mag/cell<br>     absolute magnetization    =     2.89 Bohr mag/cell<br><br>     iteration # 35     ecut=    70.00 Ry     beta= 0.70<br>     Davidson diagonalization with overlap<br>     ethr =  5.54E-10,  avg # of iterations =  4.0<br><br>     total cpu time spent up to now is    66352.1 secs<br><br>     total energy              =  -17171.63194323 Ry<br>     estimated scf accuracy    <       0.00000339 Ry<br><br>     total magnetization       =     0.00    -0.00     2.38 Bohr mag/cell<br>     absolute magnetization    =     2.89 Bohr mag/cell<br><br>     iteration # 36     ecut=    70.00 Ry     beta= 0.70<br>     Davidson diagonalization with overlap<br>     ethr =  5.54E-10,  avg # of iterations =  5.3<br><br>     total cpu time spent up to now is    68566.2 secs<br><br>     total energy              =  -17171.63217226 Ry<br>     estimated scf accuracy    <       0.00000263 Ry<br><br>     total magnetization       =    -0.00    -0.00     2.38 Bohr mag/cell<br>     absolute magnetization    =     2.87 Bohr mag/cell<br><br>     iteration # 37     ecut=    70.00 Ry     beta= 0.70<br>     Davidson diagonalization with overlap<br>     ethr =  5.18E-10,  avg # of iterations =  5.2</div><div><br></div><div>     total cpu time spent up to now is    70693.8 secs<br><br>     total energy              =  -17171.63216140 Ry<br>     estimated scf accuracy    <       0.00000268 Ry<br><br>     total magnetization       =     0.00     0.00     2.38 Bohr mag/cell<br>     absolute magnetization    =     2.87 Bohr mag/cell<br><br>     iteration # 38     ecut=    70.00 Ry     beta= 0.70<br>     Davidson diagonalization with overlap<br>     ethr =  5.18E-10,  avg # of iterations =  6.2<br><br>     total cpu time spent up to now is    73119.5 secs<br><br>     total energy              =  -17171.63204582 Ry<br>     estimated scf accuracy    <       0.00000351 Ry<br><br>     total magnetization       =    -0.00     0.00     2.38 Bohr mag/cell<br>     absolute magnetization    =     2.92 Bohr mag/cell<br><br>     iteration # 39     ecut=    70.00 Ry     beta= 0.70<br>     Davidson diagonalization with overlap<br>     ethr =  5.18E-10,  avg # of iterations =  4.5<br><br>     total cpu time spent up to now is    75281.8 secs<br><br>     total energy              =  -17171.63218978 Ry<br>     estimated scf accuracy    <       0.00000285 Ry<br><br>     total magnetization       =     0.00    -0.00     2.38 Bohr mag/cell<br>     absolute magnetization    =     2.89 Bohr mag/cell<br><br>     iteration # 40     ecut=    70.00 Ry     beta= 0.70<br>     Davidson diagonalization with overlap<br>     ethr =  5.18E-10,  avg # of iterations =  6.0<br><br>     total cpu time spent up to now is    77599.0 secs<br><br>     total energy              =  -17171.63247747 Ry<br>     estimated scf accuracy    <       0.00000310 Ry<br><br>     total magnetization       =    -0.00     0.00     2.38 Bohr mag/cell<br>     absolute magnetization    =     2.89 Bohr mag/cell<br><br>     iteration # 41     ecut=    70.00 Ry     beta= 0.70<br>     Davidson diagonalization with overlap<br></div><div>#################################################</div><div>Thank you,</div><div>Stephen<br></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><i>University of California, Berkeley</i><div><i>Department of Letter and Sciences</i></div></div></div></div>