硅烯声子谱G点附近虚频
请教各位大神,小弟入手pwscf不久,计算硅烯的声子谱,可发现两个问题实在难以解决:1、G点附近总有虚频,该怎么消除?2、在一个高对称点处两条光学支本该简并的,但是我计算出来的结果却没有,是什么原因?如下图
(硅烯的结构是组里已经优化好的结构,算声子谱直接做自洽计算吧?!)
跪求大神指导,万分感谢。
搜狗截图20160802150307.png 返回小木虫查看更多
今日热帖
请教各位大神,小弟入手pwscf不久,计算硅烯的声子谱,可发现两个问题实在难以解决:1、G点附近总有虚频,该怎么消除?2、在一个高对称点处两条光学支本该简并的,但是我计算出来的结果却没有,是什么原因?如下图
(硅烯的结构是组里已经优化好的结构,算声子谱直接做自洽计算吧?!)
跪求大神指导,万分感谢。
搜狗截图20160802150307.png 返回小木虫查看更多
scf.in参数设置
&control
calculation='scf',
restart_mode='from_scratch',
prefix='si'
pseudo_dir = '/home/bxu/zjp/2sixi/4/si/',
outdir='/home/bxu/zjp/2sixi/4/si/tmp/'
/
&system
ibrav = 4, celldm(1) =7.3052690583, celldm(3) =3.0593725672, nat= 2, ntyp= 1,
ecutwfc = 50,ecutrho=450
occupations='smearing',degauss=0.02,smearing='gaussian'
/
&electrons
electron_maxstep=200
mixing_beta = 0.5
conv_thr = 1.0d-10
diagonalization='cg'
/
ATOMIC_SPECIES
Si 28.086 Si.pbe-n-rrkjus_psl.1.0.0.UPF
ATOMIC_POSITIONS {crystal}
Si 0.333330005 0.666660011 0.480724089
Si 0.666660011 0.333330005 0.519265927
K_POINTS {automatic}
20 20 1 0 0 0,
ph.in参数设置
&inputph
tr2_ph=1.0d-16,
prefix='si',
amass(1)=28.086,
outdir='/home/bxu/zjp/2sixi/4/si/tmp',
fildyn='si.dyn',
elph=.false.,
trans=.true.,
ldisp=.true.
nq1=4, nq2=4,nq3=1
/
q2r.in 参数设置
&input
zasr='simple',fildyn='si.dyn',flfrc='si.fc',la2F=.false.
gamma 点本来就有虚频,对应长波极限。一般直接忽略不计
重新优化结构,力收敛判据调至0.0001eV/A。
改用MP型展宽,做严格收敛测试。
K点改成18*18*1.
diagonalization改成david。cg只在空带数很大的时候才有优势,速度上并不占优。
大神,你应该搞错了。。。问题不是在gamma点,是在gamma点附近的出现较大虚频,这跟文献上明显不对应啊。其次,在另一个好对称点处两光学支没有简并。是什么原因呢