[QE-users] Large estimated SCF accuracy for monolayer \alpha-RuCl3

Lingzhi ZHANG lingzhi at g.ecc.u-tokyo.ac.jp
Wed Jun 21 07:29:49 CEST 2023


Dear QE community, 

I meet a strange SCF convergence problem for my monolayer \alpha-RuCl3 system.

I’m conducting GGA+SOC+U+Constrained DFT in QE-7.0.
But the estimated SCF accuracy for my system is considerably large in the first and the following steps, more enormous than 10000 Ry. 

I have tried many methods such as modifying the mixing scheme from plain to TF, tuning the degauss, and trying different occupation and smearing schemes.
I also tried to start from the GGA+LSDA+U case, which works well and gets converged in 700 steps. 
But no scheme works for  GGA+SOC+U+Constrained DFT.
And even for just the GGA+SOC+U case, I can't obtain a small SCF accuracy.

I’d like to know what I can do for it and get a normal SCF process. 

My input is in the following. 
This structure was obtained from OpenMX, another DFT package based on the LCAO method.

&CONTROL
    calculation='scf', 
    prefix='mono', 
    pseudo_dir='/home/zhang/calculation/rucl3/pseudo', 
    restart_mode='from_scratch',
    outdir='./', 
    verbosity='high',
/
&SYSTEM
    ibrav= 0,
    nat= 16, ntyp= 3,
    occupations = 'smearing',
    smearing='mp', degauss=1d-9,
    ecutwfc= 100,
    noncolin=.true.,lspinorb=.true.,
    starting_magnetization(1)= 0.5,
    starting_magnetization(2)= 0.5,
    starting_magnetization(3)= 0.0,
    lda_plus_u = .true.,
    lda_plus_u_kind = 1,
    U_projection_type = 'ortho-atomic',
    Hubbard_U(1)=1.5,
    Hubbard_U(2)=1.5
/
&ELECTRONS
    electron_maxstep = 1600
    conv_thr = 1.0d-8
    mixing_mode = 'TF'
    mixing_beta = 4.0d-2
    mixing_ndim = 20
    startingpot = 'atomic'
    diagonalization = 'david'

/
ATOMIC_SPECIES
Ruu     101.07       Ru.rel-pbe-spn-kjpaw_psl.1.0.0.UPF
Rud     101.07       Ru.rel-pbe-spn-kjpaw_psl.1.0.0.UPF
Cl      35.45        Cl.pbe-n-kjpaw_psl.1.0.0.UPF
ATOMIC_POSITIONS {angstrom}
    Ruu     -2.4165100      0.0119000      2.4251100
    Ruu      0.0032200      0.0119000      0.0053700
    Cl       2.3833500     -0.0499700      0.0334100
    Cl      -2.3884800     -0.0499700      4.8052300
    Cl      -0.0381100      0.0580100      2.3837700
    Cl      -2.3751800     -0.0342100      0.0467100
    Cl      -2.4518700     -2.3681500      2.5014800
    Cl       0.0796000     -2.3681500     -0.0299800
    Rud      0.0105400      2.4538100     -2.4511300
    Rud      2.4302700      2.4538100     -4.8708700
    Cl       4.8104000      2.3919400     -4.8428300
    Cl       0.0385700      2.3919400     -0.0710100
    Cl       2.3889400      2.4999200     -2.4924700
    Cl       0.0518700      2.4077000     -4.8295300
    Cl      -0.0248200      0.0737600     -2.3747600
    Cl       2.5066500      0.0737600     -4.9062200
CELL_PARAMETERS {angstrom}
     4.8541000     4.8838200    -9.7524800
    -4.8762300     2.4419100     2.4270500
    10.0203300    10.0185400    10.0203300
K_POINTS {automatic}
 6 12 1 0 0 0



The following is part of the output
total energy = -2977.28363823 Ry estimated scf accuracy < 529768.11044287 Ry
total energy = -2977.28363823 Ry estimated scf accuracy < 529768.11044287 Ry
total energy = -3184.06809008 Ry estimated scf accuracy < 814811.34802588 Ry
total energy = -3274.94372181 Ry estimated scf accuracy < 746852.52429958 Ry
total energy = -3307.21972156 Ry estimated scf accuracy < 752237.62011255 Ry
total energy = -3291.26999775 Ry estimated scf accuracy < 757458.22549562 Ry
total energy = -3349.60629608 Ry estimated scf accuracy < 765501.04417510 Ry
total energy = -3337.49045478 Ry estimated scf accuracy < 768896.10678680 Ry
total energy = -3287.13640144 Ry estimated scf accuracy < 762797.72531670 Ry
total energy = -3333.32846680 Ry estimated scf accuracy < 751861.05720886 Ry
total energy = -3341.65785225 Ry estimated scf accuracy < 754604.69881806 Ry
total energy = -3350.72599779 Ry estimated scf accuracy < 749241.08182438 Ry
total energy = -3366.01560869 Ry estimated scf accuracy < 741018.08470080 Ry



With best regards

Lingzhi Zhang
Department of Applied Physics, 
University of Tokyo, Japan


More information about the users mailing list