24小时热门版块排行榜    

查看: 3472  |  回复: 10

xhzha

铁杆木虫 (正式写手)

[求助] 怎么消除虚频,细节

大家好,用siesta算声子振动谱,出现虚频,如附件图所示,怎么消除虚频呢,我看到很多人讨论过,说提高收敛精度,还有就是用虚频对应的简正坐标加到体系的坐标上去,重新优化结构,再计算声子。我精度设到10-8,已经很高,但结果还是不好,所以我想试试用虚频对应简正坐标加到体系坐标上去重新优化试试,但是这个具体操作过程我不是很了解,一个虚频的振动模式不是体系所有的原子的振动结果么?虚频的简正坐标是根据对应gamma点的虚频的频率值,解附件中的常微分方程得到的么,得到的简正坐标又加到体系的哪个原子上呢,是x,y,z的哪一个分量上呢,很多东西不是很明白,还望大家不腻赐教,谢谢!


回复此楼

» 收录本帖的淘帖专辑推荐

ATK和电子输运 碳纳米管声子

» 猜你喜欢

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

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

yinwei_li

木虫 (正式写手)

【答案】应助回帖

★ ★
xhzha(金币+5): 2011-09-20 09:53:08
zzy870720z(金币+2): 谢谢指点 2011-09-20 12:39:38
引用回帖:
1楼: Originally posted by xhzha at 2011-09-18 22:31:59:
大家好,用siesta算声子振动谱,出现虚频,如附件图所示,怎么消除虚频呢,我看到很多人讨论过,说提高收敛精度,还有就是用虚频对应的简正坐标加到体系的坐标上去,重新优化结构,再计算声子。我精度设到10-8,已 ...

有虚频说明你的结构动力学不稳定,原则上按照虚频的模式的原子振动方式移动原子(冷冻声子方法)可以找到动力学稳定的结构,但是你这个结构虚的太离谱,基本类似能带了,无法用冷冻声子计算。
虚频就虚频,为什么要消除呢?只有那些由于计算精度不够引起的虚频才能想办法消除,要是真是由于结构不稳定引起的虚频,如果你人为的消除虚频是不是有改数据的嫌疑呢?
与其临渊羡鱼,不如退而结网!
2楼2011-09-20 08:49:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xhzha

铁杆木虫 (正式写手)

引用回帖:
2楼: Originally posted by yinwei_li at 2011-09-20 08:49:11:
有虚频说明你的结构动力学不稳定,原则上按照虚频的模式的原子振动方式移动原子(冷冻声子方法)可以找到动力学稳定的结构,但是你这个结构虚的太离谱,基本类似能带了,无法用冷冻声子计算。
虚频就虚频,为什 ...

谢谢您的回复,DFT优化是在初始模型的基础上小范围移动原子找到局域内稳定结构,可能并不是体系的最稳定结构,所以我就想按照虚频的原子振动方式移动原子,改变初始构型,优化找到动力学最稳定的结构,所以用这种方法消除虚频应该还是合理的吧。但是也就是怎样去移动原子这点我比较模糊,虚频是个别原子造成的么?siesta虚频的振动方向可以可视化么,希望您能够讲述移动原子的细节。非常感谢您,我算的是纳米管的周期性结构。
3楼2011-09-20 09:52:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yinwei_li

木虫 (正式写手)

【答案】应助回帖

★ ★
zzy870720z(金币+2): 鼓励讨论 2011-09-20 12:39:51
xhzha(金币+5): 2011-09-20 16:45:32
引用回帖:
3楼: Originally posted by xhzha at 2011-09-20 09:52:38:
谢谢您的回复,DFT优化是在初始模型的基础上小范围移动原子找到局域内稳定结构,可能并不是体系的最稳定结构,所以我就想按照虚频的原子振动方式移动原子,改变初始构型,优化找到动力学最稳定的结构,所以用这种 ...

通过原子移动来找稳定结构的方法是有一定限制的,如果你的声子谱只是在布里渊区某个特殊点有虚频,找到此点对应的本征矢(本征矢可表征原子如何振动),沿着本征矢逐步移动原子并计算每移动一步体系的能量的变化,你会发现能量的变化会出现一个低谷,那么这个能量低谷的位置就对应新的稳定结构。计算这个稳定结构的声子谱应该是没有虚频的。但是你给出的声子谱在整个布里渊区都有虚频,用这种方法可能就不适应了,因为你没办法同时按照所有布里渊区点的本征矢移动原子。

或许大家有别的方法,这个我就不太清楚了!
与其临渊羡鱼,不如退而结网!
4楼2011-09-20 11:41:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yinwei_li

木虫 (正式写手)

【答案】应助回帖


zzy870720z(金币+1): 谢谢补充 2011-09-20 12:42:53
引用回帖:
3楼: Originally posted by xhzha at 2011-09-20 09:52:38:
谢谢您的回复,DFT优化是在初始模型的基础上小范围移动原子找到局域内稳定结构,可能并不是体系的最稳定结构,所以我就想按照虚频的原子振动方式移动原子,改变初始构型,优化找到动力学最稳定的结构,所以用这种 ...

补充:我做的一般都是晶体,三维的。你做的是纳米管,可能会有区别!你构造的纳米结构是不是本身就是一个近似结构?如果只是个近似结构,有虚频可能也很正常吧!
与其临渊羡鱼,不如退而结网!
5楼2011-09-20 11:43:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xhzha

铁杆木虫 (正式写手)

引用回帖:
5楼: Originally posted by yinwei_li at 2011-09-20 11:43:53:
补充:我做的一般都是晶体,三维的。你做的是纳米管,可能会有区别!你构造的纳米结构是不是本身就是一个近似结构?如果只是个近似结构,有虚频可能也很正常吧!

谢谢您的回复,我是用Material Studio 建的模型,然后用siesta去优化,但siesta只是局域内的稳定结构,所以我才想移动原子,调整结构试试。
6楼2011-09-20 16:43:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yinwei_li

木虫 (正式写手)

【答案】应助回帖


xhzha(金币+2): 2011-09-21 08:08:01
zzy870720z(金币+1): 鼓励交流 2011-09-21 13:02:56
引用回帖:
6楼: Originally posted by xhzha at 2011-09-20 16:43:13:
谢谢您的回复,我是用Material Studio 建的模型,然后用siesta去优化,但siesta只是局域内的稳定结构,所以我才想移动原子,调整结构试试。

MS建的模型是不是在一个大的框架里模拟纳米结构啊?这样的话有虚频可能很正常!因为你建的结构实际上还是有周期性的。不知道是否可以这样理解。
与其临渊羡鱼,不如退而结网!
7楼2011-09-20 19:55:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xhzha

铁杆木虫 (正式写手)

引用回帖:
7楼: Originally posted by yinwei_li at 2011-09-20 19:55:34:
MS建的模型是不是在一个大的框架里模拟纳米结构啊?这样的话有虚频可能很正常!因为你建的结构实际上还是有周期性的。不知道是否可以这样理解。

恩,我建的是周期性的结构,我的体系的键长等也是根据文献报导的结果得到的,所以初始模型可能会有点偏离。弱弱的问一句,虚频的本征矢是怎么得到的?
8楼2011-09-21 08:07:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yinwei_li

木虫 (正式写手)

【答案】应助回帖


bingmou(金币+1): 感谢交流! 2011-09-21 10:04:07
xhzha(金币+8): 2011-09-21 16:03:33
引用回帖:
8楼: Originally posted by xhzha at 2011-09-21 08:07:27:
恩,我建的是周期性的结构,我的体系的键长等也是根据文献报导的结果得到的,所以初始模型可能会有点偏离。弱弱的问一句,虚频的本征矢是怎么得到的?

siesta没用过,具体不知道该怎么找,你可以找所有输出文件,看看有没有矩阵形势的,那个可能就是。建议咨询用过siesta的专家!
与其临渊羡鱼,不如退而结网!
9楼2011-09-21 09:47:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xhzha

铁杆木虫 (正式写手)

引用回帖:
9楼: Originally posted by yinwei_li at 2011-09-21 09:47:49:
siesta没用过,具体不知道该怎么找,你可以找所有输出文件,看看有没有矩阵形势的,那个可能就是。建议咨询用过siesta的专家!

哦,呵呵,谢谢你的耐心讲解,金币都给你了~
10楼2011-09-21 16:03:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xhzha 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 307求调剂 +6 冷笙123 2026-03-17 6/300 2026-03-19 15:14 by peike
[考研] 085601材料工程专硕求调剂 +9 慕寒mio 2026-03-16 9/450 2026-03-19 15:07 by jhhcooi
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +11 yangfz 2026-03-17 11/550 2026-03-19 15:06 by houyaoxu
[考研] 0703化学调剂 +3 18889395102 2026-03-18 3/150 2026-03-19 13:53 by houyaoxu
[考研] 求材料调剂 +10 隔壁陈先生 2026-03-12 10/500 2026-03-19 09:15 by 无际的草原
[考研] 材料专硕英一数二306 +5 z1z2z3879 2026-03-18 5/250 2026-03-19 07:43 by BruceLiu320
[考研] 085700资源与环境308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-18 22:35 by bingxueer79
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +7 步川酷紫123 2026-03-13 7/350 2026-03-18 17:12 by 尽舜尧1
[考研] 295求调剂 +3 一志愿京区211 2026-03-18 5/250 2026-03-18 17:03 by zhaoqian0518
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 0854,计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[基金申请] 被我言中:新模板不强调格式了,假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 301求调剂 +4 A_JiXing 2026-03-16 4/200 2026-03-17 17:32 by ruiyingmiao
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 321求调剂(食品/专硕) +3 mxcz321 2026-03-12 6/300 2026-03-13 08:45 by xc321
信息提示
请填处理意见