24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2350  |  回复: 24
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

hanyanli0475

铜虫 (小有名气)

[求助] 声子计算又负频的出现!!!

今天计算氮化铟的声子,输出后低频有-26cm-1的负频,而且其他低频的走势和文献也不太一样,我把输入文件粘出来,请大家帮我找找错误呗,看看有没有需要修改的地方,谢谢了!!!
自洽计算:
&CONTROL
calculation = 'scf' ,
restart_mode = 'from_scratch' ,
outdir = '/home/zhao/10.21/tmp/' ,
pseudo_dir = '/opt/espresso-4.2.1/pseudo/' ,
prefix = 'InN' ,
tprnfor= .true. ,
tstress= .true. ,
/
&SYSTEM
ibrav = 4,celldm(1) = 6.893,celldm(3) = 1.605,nat = 4,ntyp = 2,
ecutwfc = 45.0 ,ecutrho = 180 ,
/
&ELECTRONS
conv_thr = 1.0d-8 ,
mixing_beta = 0.7 ,
/
ATOMIC_SPECIES
In   114.820  In.pbe-d-rrkjus.UPF
N    14.007   N.pbe-rrkjus.UPF
ATOMIC_POSITIONS crystal
In       0.333333487   0.666666513   0.003058286    1   1   1
In       0.666666513   0.333333487   0.503058286    1   1   1
N        0.333333439   0.666666561   0.381941714    1   1   1
N        0.666666561   0.333333439   0.881941714    1   1   1
K_POINTS automatic
4 4 4      1 1 1
声子计算
phonon for InN
&inputph
outdir ='/home/zhao/10.21/tmp/',
tr2_ph =1.0d-12,
prefix ='InN',
fildvscf ='InNdv',
amass(1) =114.8200,
amass(2) =14.00700,
fildyn ='InN.dyn',
elph =.false.,
trans =.true.,
ldisp =.true.,
nq1=4,nq2=4,nq3=4
/
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

什么时候能算完!!!!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hanyanli0475

铜虫 (小有名气)

引用回帖:
4楼: Originally posted by goldenfisher at 2011-10-26 18:45:29:
他用的貌似是超软的。
超软赝势ecutrho至少要是ecutwfc的8-10倍。对于rrkjus类型的,ecutwfc=45确实有些小。建议做一个能量收敛测试。
此外,4*4*4也有点小,也要做收敛测试。

请问如何做个测试呢,测试的标准怎么看呢?我以前没做过测试,请您把测试的过程说的仔细点呗,谢谢了。
什么时候能算完!!!!
5楼2011-10-26 19:45:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 25 个回答

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★
贺仪(金币+2): 多谢指教! 2011-10-26 21:59:43
4x4x4的k点也敢拿来做phonon计算啊?太大无畏了~
顺便,对于模守赝势,这点截断也忒低了~用超软做声子时这个截断还差不多。
另外,对于模守赝势,不要自己设置电荷密度截断能。
思想重于技巧,内涵重于表象
2楼2011-10-26 16:38:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hanyanli0475

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by souledge at 2011-10-26 16:38:30:
4x4x4的k点也敢拿来做phonon计算啊?太大无畏了~
顺便,对于模守赝势,这点截断也忒低了~用超软做声子时这个截断还差不多。
另外,对于模守赝势,不要自己设置电荷密度截断能。

谢谢您的回答,我现在就去改改,重新跑一下。
今天下午把ph的源程序大概看看,真难啊,呵呵!
什么时候能算完!!!!
3楼2011-10-26 17:24:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

【答案】应助回帖

★ ★
贺仪(金币+2): 多谢指教! 2011-10-26 22:00:00
他用的貌似是超软的。
超软赝势ecutrho至少要是ecutwfc的8-10倍。对于rrkjus类型的,ecutwfc=45确实有些小。建议做一个能量收敛测试。
此外,4*4*4也有点小,也要做收敛测试。
4楼2011-10-26 18:45:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 279学硕食品专业求调剂院校 20+3 孤独的狼爱吃羊 2026-04-12 14/700 2026-04-12 16:32 by 钰璞
[考研] 0831一轮调剂失败求助 +9 小熊睿睿_s 2026-04-11 9/450 2026-04-12 14:18 by 昆工高欣
[考研] 347求调剂 +4 mhyqyy 2026-04-06 4/200 2026-04-12 02:27 by 秋豆菜芽
[考研] 085404 293求调剂 +9 勇远库爱314 2026-04-08 9/450 2026-04-12 02:24 by 秋豆菜芽
[考研] 一志愿新疆大学085401,314分 +4 咔咔咔咔9 2026-04-05 4/200 2026-04-12 02:21 by 秋豆菜芽
[教师之家] 请问地理、遥感方面,可以做哪些横向项目啊,纵向完不成考核啊 +3 锦衣卫寒战 2026-04-07 5/250 2026-04-11 20:51 by 豫椒
[考研] 化学308分求调剂 +22 你好明天你好 2026-04-07 24/1200 2026-04-11 11:14 by ChemPharm
[考研] 298求调剂 +5 残荷新柳 2026-04-07 5/250 2026-04-11 11:02 by 紫曦紫棋
[考研] 调剂 +4 电气300求调剂不 2026-04-08 7/350 2026-04-11 10:44 by 紫曦紫棋
[考研] 085500求调剂材料 +10 易11122 2026-04-09 10/500 2026-04-11 10:39 by maddjdld
[考研] 机械专硕270求调剂,接受跨专业 +12 老师看看我吧aba 2026-04-09 14/700 2026-04-11 10:21 by laoshidan
[考研] 281求调剂 +11 觉得好的吧 2026-04-10 11/550 2026-04-11 09:35 by 逆水乘风
[考研] 本科211 工科085400 280分求调剂 可跨专业 +3 LZH(等待调剂中 2026-04-09 3/150 2026-04-09 21:29 by wutongshun
[考研] 332,085601求调剂 +12 ydfyh 2026-04-09 14/700 2026-04-09 17:28 by wp06
[考研] 085600材料与化工,求调剂 +11 won_qii 2026-04-07 11/550 2026-04-09 17:03 by luoyongfeng
[考研] 招收有机化学、化工,药学,食品灯专业学生 +3 yrfhjgdj 2026-04-08 3/150 2026-04-09 10:15 by QYQX_123
[考研] 0703化学调剂 348分 +14 唉我超真没招了 2026-04-06 15/750 2026-04-08 19:16 by 我减肥1
[考研] 327求调剂 +12 Xxjc1107. 2026-04-06 12/600 2026-04-08 16:46 by luoyongfeng
[考研] 电子信息346 +4 zuoshaodian 2026-04-08 4/200 2026-04-08 11:54 by zzucheup
[考研] 材料调剂 +13 汉123456 2026-04-07 14/700 2026-04-07 22:53 by 来看流星雨10
信息提示
请填处理意见