What smearing would be ideal for a bcc W with interstitial O atom. I am running 3x3x3 unit cells with one N interstitial.
I am having convergence issues in relax
and vc-relax
, but not in scf
. Could this be due to smearing? I have tried mv and cold. For an additional interstitial atom, should we change or consider charge differences? There aren't many examples of interstitial atom studies other than papers which don't mentioned anything in particular about these.
&control
calculation = 'vc-relax',
etot_conv_thr = 1e-05
forc_conv_thr = 1e-04
tprnfor = .true.
tstress = .true.
verbosity = 'high'
pseudo_dir = '.',
!prefix='w',
/
&system
ibrav = 1,
celldm(1) = 18.000,
nat = 55,
ntyp= 2,
ecutwfc = 50.0, ecutrho = 500.0,
occupations='smearing', smearing='mv', degauss=0.01,
!nbnd = 8,
!nspin=2,
starting_magnetization(1)= 0.6
/
&electrons
mixing_beta = 0.7
conv_thr=1e-6
/
&ions
ion_dynamics='bfgs'
/
&cell
cell_dynamics='bfgs'
cell_dofree='volume'
/
ATOMIC_SPECIES
W 183.84 W.pbe-spn-rrkjus_psl.1.0.0.UPF
O 15.999 O.pbe-n-rrkjus_psl.1.0.0.UPF
ATOMIC_POSITIONS crystal
W 0.000000 0.000000 0.000000
W 0.166667 0.166667 0.166667
W 0.000000 0.000000 0.333333
W 0.166667 0.166667 0.500000
W 0.000000 0.000000 0.666667
W 0.166667 0.166667 0.833333
W 0.000000 0.333333 0.000000
W 0.166667 0.500000 0.166667
W 0.000000 0.333333 0.333333
W 0.166667 0.500000 0.500000
W 0.000000 0.333333 0.666667
W 0.166667 0.500000 0.833333
W 0.000000 0.666667 0.000000
W 0.166667 0.833333 0.166667
W 0.000000 0.666667 0.333333
W 0.166667 0.833333 0.500000
W 0.000000 0.666667 0.666667
W 0.166667 0.833333 0.833333
W 0.333333 0.000000 0.000000
W 0.500000 0.166667 0.166667
W 0.333333 0.000000 0.333333
W 0.500000 0.166667 0.500000
W 0.333333 0.000000 0.666667
W 0.500000 0.166667 0.833333
W 0.333333 0.333333 0.000000
W 0.500000 0.500000 0.166667
W 0.333333 0.333333 0.333333
W 0.500000 0.500000 0.500000
W 0.333333 0.333333 0.666667
W 0.500000 0.500000 0.833333
W 0.333333 0.666667 0.000000
W 0.500000 0.833333 0.166667
W 0.333333 0.666667 0.333333
W 0.500000 0.833333 0.500000
W 0.333333 0.666667 0.666667
W 0.500000 0.833333 0.833333
W 0.666667 0.000000 0.000000
W 0.833333 0.166667 0.166667
W 0.666667 0.000000 0.333333
W 0.833333 0.166667 0.500000
W 0.666667 0.000000 0.666667
W 0.833333 0.166667 0.833333
W 0.666667 0.333333 0.000000
O 0.416666 0.500000 0.500000
W 0.833333 0.500000 0.166667
W 0.666667 0.333333 0.333333
W 0.833333 0.500000 0.500000
W 0.666667 0.333333 0.666667
W 0.833333 0.500000 0.833333
W 0.666667 0.666667 0.000000
W 0.833333 0.833333 0.166667
W 0.666667 0.666667 0.333333
W 0.833333 0.833333 0.500000
W 0.666667 0.666667 0.666667
W 0.833333 0.833333 0.833333
! this is a comment that the code will ignore
K_POINTS automatic
4 4 4 1 1 1