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.

求指教！

计算值

实验值

回复此楼

» 猜你喜欢

读博已经有5人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有13人回复
博士申请都是内定的吗？已经有6人回复
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？已经有5人回复
博士读完未来一定会好吗已经有29人回复
投稿精细化工已经有4人回复
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！已经有4人回复
导师想让我从独立一作变成了共一第一已经有9人回复
心脉受损已经有5人回复
Springer期刊投稿求助已经有4人回复

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

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

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

tinkie0411

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 387.7
散金: 13
帖子: 95
在线: 77.9小时
虫号: 1074695
注册: 2010-08-14
专业: 金属功能材料

【答案】应助回帖

引用回帖:

7楼: Originally posted by leila3409 at 2013-12-23 08:11:20
嗯，很可能是晶胞不够大，我昨天也想到了。我用原胞在各方向扩大5倍，但因为基失不正交，所以由周期性引起的晶胞间相互作用可能还是比较大，我已经将扩展胞正交，现在有250个原子了。
energycutoff 会对别的性质有 ...

请问楼主你将扩展包变为正交后与之前的声子谱有差别吗？谢谢