7
$\begingroup$

I used VASP as my DFT software to do a calculation on an anti-perovskite crystal. I try to reproduce the band structure from literature which was calculated based on HSE + SOC calculation. Their result was shown below:

]2

Source: https://www.researchgate.net/publication/264084992_Topological_Crystalline_Insulators_and_Dirac_Octets_in_Anti-perovskites

10.1103/PhysRevB.90.081112

I try to redo the calculation and plot the band along the same high symmetry line, here is what I got ( My K-point mesh is a bit low, so it may not be so smooth) :

enter image description here

It seems that the two graphs are qualitatively the same. But the gap at Gamma point was open up. I have no idea what kind of mistake I've made. And here is my INCAR setting in my calculation:

#Calculation Parameter

 ISMEAR = 0 ; SIGMA = 0.05

 PREC = Normal

 ENCUT = 350

 ISPIN = 2

 LWAVE = .TRUE.

 LCHARG = .TRUE.


 LSORBIT = .TRUE.

 MAGMOM = 15*0

 SAXIS = 0 0 1

 GGA_COMPAT = .FALSE.

#Electronic Relaxation

 EDIFF = 1E-05

 ALGO = Normal

 NELM = 500

#Ionic Relaxation

 IBRION = -1

 NSW = 0

#HSE06

 LHFCALC = .TRUE.

 HFSCREEN = 0.2

 ALGO = Damped

 TIME = 0.4

 #AEXX = 0.25
$\endgroup$
4
  • $\begingroup$ You should give the link to the reference paper. $\endgroup$ – Jack Dec 16 '20 at 4:31
  • $\begingroup$ I’ve already added the reference ! Thanks for your reminder $\endgroup$ – JensenPang Dec 16 '20 at 5:12
  • $\begingroup$ You are not using the same parameters. $\endgroup$ – Jack Dec 16 '20 at 5:26
  • $\begingroup$ I gave you +1 long ago, but does there need to be so much space between each line in what is now the code block? $\endgroup$ – Nike Dattani Feb 28 at 5:08

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.