大家好,最近我在重复石墨烯纳米带的电子结构,我已经得到了能带和文献的差不多,但是在计算态密度时,发现PDOS之和不等于DOS,我已经把所有得到的每个原子的s和p态的态密度分别加起来了,但是,还是相差的太多,是不是我输入文件又问题啊,请高手们指点一下啊,真是无语啊,往下的工作就卡在这了,谢谢大家啊。下面是输入文件:
1.scf
&CONTROL
calculation = 'scf' ,
restart_mode = 'from_scratch' ,
outdir = '/' ,
pseudo_dir = '/' ,
prefix = 'GNRs' ,
etot_conv_thr = 1.0D-5 ,
forc_conv_thr = 1.0D-4 ,
tstress = .true. ,
tprnfor = .true. ,
/
&SYSTEM
ibrav = 8,
celldm(1) = 12.85444,celldm(2) = 2.64706,celldm(3) =0.62131,
nat = 18,ntyp = 2,
ecutwfc = 125.0 ,ecutrho = 500 ,nosym = .false.,
occupations = 'smearing' ,degauss = 0.02 ,smearing = 'mp' ,
/
&ELECTRONS
conv_thr = 1.0d-8 ,
mixing_beta = 0.7 ,
mixing_mode = 'plain' ,
diagonalization = 'david',
/
ATOMIC_SPECIES
C 12.00000 C.pz-vbc.UPF
H 1.00000 H.pz-vbc.UPF
ATOMIC_POSITIONS crystal
C 0.866030546 0.704408575 0.695292364
C 0.866030398 0.300318030 0.695284939
C 0.866028084 0.434371782 0.702120006
C 0.866028204 0.570355384 0.702121677
C 0.866028670 0.637974084 0.202842982
……………………………………
H 0.866033838 0.247278886 0.250636660
K_POINTS automatic
1 1 21 0 0 0
2.nscf
&CONTROL
calculation = 'nscf' ,
restart_mode = 'from_scratch' ,
outdir = './' ,
pseudo_dir = '。/' ,
prefix = 'GNRs' ,
tprnfor= .true. ,
tstress= .true. ,
etot_conv_thr= 1.0e-5,
forc_conv_thr= 1.0e-4,
/
&SYSTEM
ibrav = 8,celldm(1) = 12.85444,celldm(2) =2.64706,celldm(3) =0.62131,
nat = 18,ntyp = 2,
ecutwfc = 125.0 ,ecutrho = 500 ,nosym = .false., nbnd=80,
/
&ELECTRONS
conv_thr = 1.0d-8 ,
mixing_beta = 0.7 ,
mixing_mode = 'plain' ,
diagonalization = 'david',
/
ATOMIC_SPECIES
C 12.00000 C.pz-vbc.UPF
H 1.00000 H.pz-vbc.UPF
ATOMIC_POSITIONS crystal
C 0.866030546 0.704408575 0.695292364
C 0.866030398 0.300318030 0.695284939
C 0.866028084 0.434371782 0.702120006
C 0.866028204 0.570355384 0.702121677
…………
H 0.866033669 0.247279466 0.819441547
H 0.866033841 0.757447327 0.250639886
H 0.866033838 0.247278886 0.250636660
K_POINTS automatic
1 1 300 0 0 0
3.dos.x
&inputpp
outdir='./'
prefix='GNRs'
fildos='GNRs.dos',
Emin=-15.0, Emax=25.0, DeltaE=0.01
/
4.pdos.x
&inputpp
outdir='./'
prefix='GNRs'
Emin=-20.0, Emax=25.0, DeltaE=0.01 ngauss=1, degauss=0.02
/
图已给出,请打家帮帮忙呗,谢谢哈!
![]()
DOS.jpg
![]()
pdos.jpg |