24小时热门版块排行榜    

查看: 5136  |  回复: 23

gemucai

木虫 (正式写手)

[求助] 哭了,各种虚频!

我以前计算声子谱很少出现虚频,但最近各种虚频都出现了,十分的崩溃!
首先是这个非高对称点的虚频:




然后是无限长的烷烃声子谱也出现了虚频,不做asr的话,gamma点也有很大虚频:


两个结构的受力我优化到0.00001 ry/bohr,stress优化到0.5kbar以内,且ecut和k点也做过收敛测试的。我已经束手无策了,谁能帮帮我?



回复此楼

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

VASP and MS condensed matter physics

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

gemucai

木虫 (正式写手)

问题解决了。

跟大家分享一下我的解决方法,以后碰到有虚频的情况不妨增加smearing的宽度,即degauss的大小,效果很灵的。其实有时候根本不需要过度地优化,我优化到0.0000001ry/bohr也有很大虚频,但是degauss增加了0.01虚频就立刻消失。

希望这对大家以后的计算起到帮助的作用,也不妨讨论下为何虚频对degauss如此敏感。
3楼2012-05-25 18:59:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hanyanli0475

铜虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
WDD880227: 金币+1, 多谢交流提示 2012-05-24 16:54:03
gemucai: 金币+5, ★★★很有帮助, 超软和模守恒我都用过,会不会一维结构更容易失稳? 2012-05-25 08:01:54
以前我算声子也遇到虚频过,尤其是用超软贗势时,虚频非常大,我当时,的解决办法就是,从头做收敛测试,然后优化,都是重新做了一遍,最后就没负频了,所以,个人认为,声子对对称性要求非常高,所以,收敛测试精度高一些,优化时,步长密一些,保证结构非常稳定吧。个人建议,请指教。
什么时候能算完!!!!
2楼2012-05-24 15:58:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

姚海峰456

银虫 (初入文坛)

已经用第一性原理优化好的单胞,用gulp计算声子谱会出现虚频,怎么解决?——求大神帮忙啊!!!
12楼2012-08-01 21:42:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

souledge

专家顾问 (著名写手)

引用回帖:
3楼: Originally posted by gemucai at 2012-05-25 18:59:36
问题解决了。

跟大家分享一下我的解决方法,以后碰到有虚频的情况不妨增加smearing的宽度,即degauss的大小,效果很灵的。其实有时候根本不需要过度地优化,我优化到0.0000001ry/bohr也有很大虚频,但是degauss增 ...

LZ做的是导电体系么?至于smearing值的大小,也要进行收敛测试的……
思想重于技巧,内涵重于表象
4楼2012-05-25 21:53:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

condensed

木虫 (著名写手)

显然,含氢材料
5楼2012-05-26 03:16:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

引用回帖:
4楼: Originally posted by souledge at 2012-05-25 21:53:28
LZ做的是导电体系么?至于smearing值的大小,也要进行收敛测试的……...

绝缘体哦。。。
6楼2012-05-26 07:56:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

我怀疑是这样的,smearing取得过小会导致电子密度n(r)出现不该有的尖峰,本来是线性的分布结果变成了spiky的形状,那么求泛函和偏导的时候就会出现不该有的震荡,甚至是很大的虚频。
7楼2012-05-26 08:00:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

just_play

至尊木虫 (正式写手)

请教楼主绝缘体系为什么不用fixed占据方法?效果不好吗?
So Trivial !
8楼2012-05-26 10:15:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

引用回帖:
8楼: Originally posted by just_play at 2012-05-26 10:15:12
请教楼主绝缘体系为什么不用fixed占据方法?效果不好吗?

我一直以为这个occupation不重要,现在发现还是需要在上面花点力气的。什么情况下用fixed,有哪些优点,能告诉我吗?
9楼2012-05-26 10:26:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

刚刚看了PW的输入文件,里面提到fixed适用于绝缘体,我可能是弄混了。以前我设置过别的参数也是fixed,结果他们要我输入占据数,我输入了好几次系统都提示我说是错的,所以我后来都不用fixed了。但今天看来好像不是这个参数。我后来算声子都用smearing,也没发现有什么异样。。。
10楼2012-05-26 10:33:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gemucai 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 岩土工程学报什么时候才能终审完呐 25+3 yeager111 2026-08-10 4/200 2026-08-16 20:42 by tfang
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +9 Ldrop2023 2026-08-13 9/450 2026-08-16 13:24 by wukun_007
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +11 家与远方 2026-08-10 16/800 2026-08-16 10:28 by ray43
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
信息提示
请填处理意见