7
$\begingroup$

I am working on $\ce{GeHfTe}$ and in the case of a non-collinear calculation, I am getting the DOS in the attached photo. I have already relaxed my structure with VC-relax and it is converged with respect to k-points, EOS parameter, and degauss value.

If any one can guide me what is wrong in the DOS obtained it will be very helpful. I am not quite sure about the two sharp peaks around $\pu{-25eV}$.

DOS of GeHfTe

&CONTROL
  calculation = 'scf'
  restart_mode  = 'from_scratch'
  etot_conv_thr =   1.0000000000d-13
  forc_conv_thr =   1.0000000000d-09
  outdir = './'
  prefix = 'GeHfTe'
  pseudo_dir = '/home/yash/project/qe-6.7-ReleasePack/qe-6.7/pseudo'
  tprnfor = .true.
  tstress = .true.
  verbosity = 'high'
/
&SYSTEM
  degauss =   0.001
  ecutwfc =   150
  ibrav = 0
  nat = 6
  ntyp = 3
  occupations = 'smearing'
  smearing = 'cold'
 noncolin=.true.
  lspinorb=.true.
/
&ELECTRONS
  conv_thr =   1.2000000000d-15
  mixing_beta =   4.0000000000d-01
/
&IONS
/
&CELL
/
ATOMIC_SPECIES
Ge     72.64  Ge.rel-pz-dn-rrkjus_psl.0.2.2.UPF
Hf     178.49 Hf.rel-pz-spn-rrkjus_psl.1.0.0.UPF
Te     127.6 Te.rel-pz-dn-rrkjus_psl.0.2.2.UPF

CELL_PARAMETERS (angstrom)
   3.799849132   0.000000000   0.000000000
   0.000000000   3.799849132   0.000000000
   0.000000000   0.000000000   8.387787463

ATOMIC_POSITIONS (crystal)
Hf            0.0000000000        0.5000000000        0.7536986859
Hf            0.5000000000        0.0000000000        0.2463013141
Ge            0.0000000000        0.0000000000       -0.0000000000
Ge            0.5000000000        0.5000000000        0.0000000000
Te            0.5000000000        0.0000000000        0.6218959934
Te            0.0000000000        0.5000000000        0.3781040066

K_POINTS (automatic)
10 10 10 0 0 0
$\endgroup$
9
  • 2
    $\begingroup$ Perhaps you could show what input file you used to get this DOS? It may be that there is nothing wrong but you have a low-lying state that isn't dispersed in energy. It all depends ;) $\endgroup$
    – nickpapior
    Oct 26, 2021 at 7:00
  • $\begingroup$ Thank you so much here is the input scf file and out put [link] (easyupload.io/m/5g6uln) and in dos i am using guassian smearing. $\endgroup$
    – Yash Ghoda
    Oct 26, 2021 at 7:07
  • $\begingroup$ I think you should just put it in the question (surrounded by triple back-ticks) $\endgroup$
    – nickpapior
    Oct 26, 2021 at 10:47
  • 5
    $\begingroup$ Frequently occurs with semi-core states. They'll be localized and well below the Fermi energy. $\endgroup$ Oct 26, 2021 at 18:43
  • 1
    $\begingroup$ I have a few questions/suggestions for troubleshooting this. Your SCF input looks "OK" (tight thresholds, you said your degauss was optimized) but I wonder... Were you able to get the DOS without spin orbit coupling (SOC)? Try to calculate the SCF, NSCF, and DOS with scalar-relativistic pseudopotentials instead of fully relativistic. It will be easier to solve if we simplify your calculation. And, also, I think spin-polarized should be fine for DOS, and SOC effect for the band structure. $\endgroup$
    – epalos
    Oct 30, 2021 at 9:27

1 Answer 1

8
$\begingroup$

I suspect that these are the $5s$ and $5p$ states of $\ce{Hf}$, which are treated as valence in the pseudopotential you're using. These states will be fairly tightly bound to the nucleus, which always show up as sharp peaks in the DOS.

The reason why tightly-bound states show up like this can be understood by noting that if the potential $V$ is large compared to the kinetic term $T$, then the Hamiltonian $H$ can be approximated as:

$$H\approx V$$

The energy $E$ of a wavefunction $\psi_{k}$ at k-point $k$ is,

\begin{align} E & \approx \int d^3 r \psi^\ast_k(r)V(r)\psi_k(r)\\ &= \int d^3 r V(r)\psi^\ast_k(r)\psi_k(r)\\ &= \int d^3 r V(r)\left\vert \psi_k(r)\right\vert^2\\ \end{align}

Now, recall that k-points are introduced when we apply Bloch's theorem, and write,

$$\psi_k(r)=u_k(r)e^{ik\cdot r}$$

where $u_k$ is the Bloch function, which does not explicitly depend on $k$. Substituting into the energy expression, we have

\begin{align} E & \approx \int d^3 r V(r)\left\vert u_k(r)\right\vert^2\\ \end{align}

In other words, the phase factor is irrelevant and the dependence on the phase $k$ has disappeared! The energy does not depend on $k$ at all, and the state has the same energy regardless of $k$. This means that the density of states for this state is a delta-function at $E$ -- which is essentially what you're seeing in your DOS plots.

$\endgroup$
3
  • $\begingroup$ How can we solve this problem please and avoid getting a similar DOS? I'm facing a similar issue $\endgroup$
    – Camilla
    Apr 12, 2023 at 13:27
  • 1
    $\begingroup$ @camilla if it's because of semi-core states, like the example here, then there's nothing wrong - they are correct. The DOS of a localised state really is a sharp spike. $\endgroup$ Apr 12, 2023 at 15:06
  • $\begingroup$ @Camilla ,My guide told me that, this is fine and could be a reason of Pseudo potential that i have used here. I have stopped working on this system as of now. so i do not know much about this as of now, although if you are working on same system i would like to know if you get any other explanation regarding the issue. $\endgroup$
    – Yash Ghoda
    Apr 13, 2023 at 16:00

You must log in to answer this question.

Not the answer you're looking for? Browse other questions tagged .