24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

leila3409

捐助贵宾 (正式写手)

应助: 0 (幼儿园)
金币: 616.2
散金: 406
红花: 7
帖子: 310
在线: 71.5小时
虫号: 985885
注册: 2010-03-30
性别: MM
专业: 爆炸与冲击动力学

[求助] 声子谱虚频求助已有4人参与

本人用SIESTA算Nb的声子谱，扩散晶胞为原胞的5*5*5倍，单边长度为3.2985*5=16.5 Ang，共125个原子。EngeryCutoff 选1200 Ry，K点用6*6*6 ，MP展宽用0.3 eV，电子收敛标准为 1.d-5，原子移动位移为0.04 Bohr。

计算的结果总在GAMMA点附近出现虚频。如下图及与实验值的对比。正确的结果应该是在Gamma附近的声学支出现一个小节，kink，但不应该是虚频。。。

本人所用晶格常数为SIESTA优化后的结果3.2985 Ang,实验值为3.3Ang.

求指教！

计算值

实验值

回复此楼

» 猜你喜欢

论文终于录用啦！满足毕业条件了已经有13人回复
2025年遐想已经有5人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复
求个博导看看已经有18人回复

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

1楼 2013-12-20 18:31:55

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

xhzha

铁杆木虫 (正式写手)

应助: 37 (小学生)
金币: 9454.6
散金: 30
红花: 7
帖子: 409
在线: 577.7小时
虫号: 1101711
注册: 2010-09-17
性别: GG
专业: 凝聚态物性I:结构、力学和

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
leila3409: 金币+50, ★★★很有帮助 2013-12-23 08:00:46
franch: 金币+2, 谢谢回帖交流 2013-12-23 21:34:33

EnergyCutoff 的值好大，有验证需要那么大么，不浪费时间么？没有看到你力的收敛标准，可以在优化的时候将力的收敛标准设高一些。如若再出现虚频，可以考虑将你的虚频的对应的坐标（对应在Eigenvector里)加到原子坐标优化再计算。还有一种是优化比较大的超胞，然后用超胞算声子谱。如果都不行，我也没辙了~