Smearing and charge balance for tungsten with interstitial oxygen in Quantum ESPRESSO

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

• I don't know much about QE, but W psuedopotentials tend to be difficult (in my experience). You may need more valence electrons to stabilize your calculation maybe? – Tristan Maxson Jan 2 at 19:23
• For studying effects due to an interstitial atom, I'm not sure as to why you've chosen a vc-relax calculation. Usually, the unit cell of the unaltered crystal will be subjected to vc-relax and with the optimized unit-cell parameters you would construct the supercell. Include the interstitial atom and then do a "relax calculation". Atleast this is how I've seen it being done. – Anoop A Nair Jan 2 at 20:13
• Since youre having convergence issues, could you try to change the default diagonalization form "david" to "cg" using diagonalization = "cg" in your &electrons namelist. It has helped me sometimes atleast. – Anoop A Nair Jan 2 at 20:17
• @Dr.Viper I think his main point was that you should relax the pure cell before addition, but also keep it fixed after addition. This assumes you want to model the dopant in its non-interacting dilute limit rather than relaxing the crystal to account for the dopant being incorporated uniformly. Think of surface calculations where the cell dimensions do not change, even after absorption of an adsorbate if this helps you think about it. The adsorbate in high coverage may play a slight role in the surface cell, but in the dilute limit it should definitely not influence things. – Tristan Maxson Jan 2 at 22:54
• @Dr.Viper Scratch my suggestion, you already are using a spn potential which means s and p states are involved as valence. – Tristan Maxson Jan 3 at 0:11