| 查看: 3301 | 回复: 10 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
xhzha铁杆木虫 (正式写手)
|
[求助]
怎么消除虚频，细节
|
|
大家好，用siesta算声子振动谱，出现虚频，如附件图所示，怎么消除虚频呢，我看到很多人讨论过，说提高收敛精度，还有就是用虚频对应的简正坐标加到体系的坐标上去，重新优化结构，再计算声子。我精度设到10-8，已经很高，但结果还是不好，所以我想试试用虚频对应简正坐标加到体系坐标上去重新优化试试，但是这个具体操作过程我不是很了解，一个虚频的振动模式不是体系所有的原子的振动结果么？虚频的简正坐标是根据对应gamma点的虚频的频率值，解附件中的常微分方程得到的么，得到的简正坐标又加到体系的哪个原子上呢，是x,y,z的哪一个分量上呢，很多东西不是很明白，还望大家不腻赐教，谢谢!  |
回复此楼» 猜你喜欢
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有7人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有8人回复
-
申请2026年博士
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有7人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
z841888077
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 730.9
- 散金: 21
- 红花: 1
- 帖子: 93
- 在线: 76.5小时
- 虫号: 2330591
- 注册: 2013-03-08
- 专业: 无机纳米化学
11楼2018-05-02 20:08:16
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 (小学生)
- 金币: 9466.1
- 散金: 30
- 红花: 7
- 帖子: 407
- 在线: 577.2小时
- 虫号: 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