应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 声子计算又负频的出现！！！
汕头大学海洋科学接受调剂
53/1返回列表
查看: 2360  |  回复: 24
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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
/
回复此楼

» 猜你喜欢
什么时候能算完！！！！
1楼 2011-10-26 15:05:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huazhorg

铁杆木虫 (正式写手)

【答案】应助回帖

★
贺仪(金币+1): 多谢指教！ 2011-11-02 21:05:01
引用回帖:
23楼: Originally posted by hanyanli0475 at 2011-11-02 19:10:05:
NC当时没有做收敛测试用了45ry，超软后来做了收敛测试用了70ry。

你这顺序完全颠倒了，超软的用小点的就够了，NC的用大的！难怪你的超软反而比NC时间长。以我的经验，超软一般30～40就够了，NC赝势本身会有个推荐值，一般会在80～90以上。你用45Ry的给NC，显然小了，结果要是正常那才怪。
赞一下(1人) 回复此楼
24楼2011-11-02 19:32:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 25 个回答

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★
贺仪(金币+2): 多谢指教！ 2011-10-26 21:59:43
4x4x4的k点也敢拿来做phonon计算啊？太大无畏了～
顺便，对于模守赝势，这点截断也忒低了～用超软做声子时这个截断还差不多。
另外，对于模守赝势，不要自己设置电荷密度截断能。
赞一下(1人) 回复此楼
思想重于技巧，内涵重于表象
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也有点小，也要做收敛测试。
赞一下(1人) 回复此楼
4楼2011-10-26 18:45:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 25 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 290调剂生物0860 +35 哇哈哈，。 2026-04-11 41/2050 2026-04-14 23:47 by Xurambo2014
[考研] 化学070300 求调剂 +23 哈哈哈^_^ 2026-04-12 23/1150 2026-04-14 16:30 by zhouxiaoyu
[考研] 297，工科调剂? +6 河南农业大学-能 2026-04-14 6/300 2026-04-14 16:04 by Art1977
[考研] 一志愿085502，267分求调剂 +19 再忙也要吃饭啊 2026-04-08 20/1000 2026-04-14 16:03 by zs92450
[考研] 296求调剂 +8 汪！？！ 2026-04-09 9/450 2026-04-14 15:15 by 逆水乘风
[考研] 085408光电信息工程专硕355一志愿长春光机所调剂 +6 王ymaa 2026-04-13 13/650 2026-04-14 11:33 by 王ymaa
[考研] 机械还有还有名额吗？太难了 +8 笑笑袁 2026-04-10 8/400 2026-04-14 08:44 by screening
[考研] 297工科，求调剂? +13 河南农业大学-能 2026-04-12 13/650 2026-04-13 14:12 by dingyanbo1
[考研] 求调剂 +16 张番茄不炒蛋 2026-04-10 17/850 2026-04-12 13:58 by 熬夜成！
[考研] 291求调剂 +8 关忆北. 2026-04-11 8/400 2026-04-12 09:32 by 逆水乘风
[考研] 331求调剂 +5 王国帅 2026-04-11 5/250 2026-04-11 22:56 by 溪涧流水
[考研] 086003调剂求助 +21 苏弋万 2026-04-09 22/1100 2026-04-11 20:25 by dongdian1
[考研] 282，求调剂 +12 jggshjkkm 2026-04-09 14/700 2026-04-11 09:39 by 猪会飞
[考研] 0854调剂 +4 长弓傲 2026-04-09 4/200 2026-04-11 09:18 by 猪会飞
[考研] 求调剂 +5 不会飞的鱼@ 2026-04-10 5/250 2026-04-10 19:07 by chemisry
[考研] 282，电气工程专业，求调剂，不挑专业 +9 jggshjkkm 2026-04-10 9/450 2026-04-10 14:55 by 逆水乘风
[考研] 求调剂 材料与工程 324分 专硕 +19 翩翩一书生 2026-04-10 21/1050 2026-04-10 11:41 by wp06
[考研] 296求调剂 +6 汪！？！ 2026-04-08 6/300 2026-04-10 11:02 by mattzhming
[考研] 本科西工大 0856 324求调剂 +10 wysyjs25 2026-04-09 11/550 2026-04-10 08:37 by 5268321
[考研] 复试调剂，一志愿郑州大学材料与化工289分 +31 硕星赴 2026-04-08 31/1550 2026-04-09 16:54 by Delta2012
信息提示
关闭
请填处理意见
关闭