i use qe to calculate the bands of mgb2 ,and get a semiconductor result with gap.but from the paper i know it is a superconductor, i dont know whats wrong with the parameter i choose?

 &control

calculation = 'scf'

restart_mode='from_scratch',

prefix='mgb2',

tstress = .true.

tprnfor = .true.

pseudo_dir = '****/qe-6.1/pseudo/',

outdir='*****/tempdir/'

/

&system

ibrav= 4, celldm(1) =5.826027964, celldm(3)=1.142  nat= 3, ntyp= 2,

nbnd =6,

ecutwfc =30.0,

ecutrho = 72,

occupations='smearing',

degauss=0.01,

smearing='mv'

/

&electrons

diagonalization='david'

mixing_mode = 'plain'

mixing_beta = 0.7

conv_thr = 1.0d-8

/

ATOMIC_SPECIES

mg 24.305 Mg.pz-n-vbc.UPF

b 10.811 B.pz-vbc.UPF

ATOMIC_POSITIONS {crystal}

mg 0.000000 0.000000 0.000000

b 0.333333 0.666667 0.500000

b 0.666667 0.333333 0.500000

K_POINTS {automatic}

24 24 24 0 0 0

thanx

More Cheedge Lee's questions See All
Similar questions and discussions