| 查看: 3767 | 回复: 10 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
xhzha铁杆木虫 (正式写手)
|
[求助]
怎么消除虚频，细节
|
|||
大家好，用siesta算声子振动谱，出现虚频，如附件图所示，怎么消除虚频呢，我看到很多人讨论过，说提高收敛精度，还有就是用虚频对应的简正坐标加到体系的坐标上去，重新优化结构，再计算声子。我精度设到10-8，已经很高，但结果还是不好，所以我想试试用虚频对应简正坐标加到体系坐标上去重新优化试试，但是这个具体操作过程我不是很了解，一个虚频的振动模式不是体系所有的原子的振动结果么？虚频的简正坐标是根据对应gamma点的虚频的频率值，解附件中的常微分方程得到的么，得到的简正坐标又加到体系的哪个原子上呢，是x,y,z的哪一个分量上呢，很多东西不是很明白，还望大家不腻赐教，谢谢!  |
回复此楼» 收录本帖的淘帖专辑推荐
 ATK和电子输运 | 碳纳米管声子 |
» 猜你喜欢
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有72人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
-
招聘 金属平磨液，抛光液研发工程师
已经有3人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何查看计算精度？
已经有16人回复
- 虚频 是怎么回事？
已经有7人回复
- 声子谱图（速度自相关函数的傅里叶变换之类的问题）
已经有9人回复
- 消除七个虚频的方法
已经有20人回复
- 碳化硅细粉注浆成型
已经有15人回复
- 过渡态中怎样消除不必要的虚频
已经有3人回复
- 哭了，各种虚频！
已经有23人回复
- MS中如何调整键长？
已经有6人回复
- 如何改变能带结构图的精度
已经有16人回复
- 急求细节！学院的副院长用assistant dean 对不对吗？急着封装推荐信。。。
已经有8人回复
- 【求助】为什么有虚频的结构能量比没有虚频的结构能量低？
已经有6人回复
- 【求助】关于虚频
已经有8人回复
- 【求助】请教CASTEP过渡态里面如何消除虚频？
已经有16人回复
- 【求助】siesta计算声子
已经有16人回复
- 【求助】再怎么消虚频啊？
已经有10人回复
- 【求助】虚频消除问题--虚频振动坐标能找到，怎么找不到相应的原坐标
已经有6人回复
- 【求助】消除虚频问题
已经有8人回复
xhzha
铁杆木虫 (正式写手)
- 应助: 37 (小学生)
- 金币: 9454.6
- 散金: 30
- 红花: 7
- 帖子: 409
- 在线: 577.7小时
- 虫号: 1101711
- 注册: 2010-09-17
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
10楼2011-09-21 16:03:19
yinwei_li
木虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 2079.1
- 散金: 1532
- 红花: 10
- 帖子: 896
- 在线: 194.2小时
- 虫号: 443243
- 注册: 2007-10-27
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
2楼2011-09-20 08:49:11
xhzha
铁杆木虫 (正式写手)
- 应助: 37 (小学生)
- 金币: 9454.6
- 散金: 30
- 红花: 7
- 帖子: 409
- 在线: 577.7小时
- 虫号: 1101711
- 注册: 2010-09-17
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
3楼2011-09-20 09:52:38
yinwei_li
木虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 2079.1
- 散金: 1532
- 红花: 10
- 帖子: 896
- 在线: 194.2小时
- 虫号: 443243
- 注册: 2007-10-27
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
【答案】应助回帖
★ ★
zzy870720z(金币+2): 鼓励讨论 2011-09-20 12:39:51
xhzha(金币+5): 2011-09-20 16:45:32
zzy870720z(金币+2): 鼓励讨论 2011-09-20 12:39:51
xhzha(金币+5): 2011-09-20 16:45:32
|
通过原子移动来找稳定结构的方法是有一定限制的，如果你的声子谱只是在布里渊区某个特殊点有虚频，找到此点对应的本征矢（本征矢可表征原子如何振动），沿着本征矢逐步移动原子并计算每移动一步体系的能量的变化，你会发现能量的变化会出现一个低谷，那么这个能量低谷的位置就对应新的稳定结构。计算这个稳定结构的声子谱应该是没有虚频的。但是你给出的声子谱在整个布里渊区都有虚频，用这种方法可能就不适应了，因为你没办法同时按照所有布里渊区点的本征矢移动原子。 或许大家有别的方法，这个我就不太清楚了！ |
4楼2011-09-20 11:41:18
