24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 4741  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂,一志愿 南京航空航天大学 ,080500材料科学与工程学硕,总分289分 +10 @taotao 2026-03-29 10/500 2026-04-01 09:30 by oooqiao
[考研] 08工科275求调剂,可跨考。 +3 AaAa7420 2026-03-31 3/150 2026-04-01 09:27 by JourneyLucky
[考博] 26申博 +4 加油冲啊! 2026-03-26 4/200 2026-03-31 22:42 by greychen00
[考研] 一志愿北化085600材料专硕275|有文章专利|求调剂 +16 Micky11223 2026-03-25 17/850 2026-03-31 22:04 by 544594351
[考研] 一志愿西交大080500材料学硕349 +6 jqx1258 2026-03-31 7/350 2026-03-31 21:08 by yuq
[考研] 318求调剂 +3 笃行致远. 2026-03-31 3/150 2026-03-31 20:27 by 求调剂zz
[考研] 考研调剂 +9 小蜡新笔 2026-03-29 10/500 2026-03-31 19:52 by Dyhoer
[考研] 化学308分调剂 +11 你好明天你好 2026-03-30 12/600 2026-03-31 19:14 by 杨金金金
[考研] 070300一志愿211,312分求调剂院校 +12 小黄鸭宝 2026-03-30 12/600 2026-03-31 17:53 by 544594351
[考研] 086000生物与医药 初试274求调剂 +4 小叮当来了 2026-03-30 4/200 2026-03-31 16:48 by shengliu165
[考研] 317求调剂 +12 蛋黄咸肉粽 2026-03-26 13/650 2026-03-31 16:44 by asdfzly
[考研] 理学07化学 303求调剂 +8 睿08 2026-03-27 8/400 2026-03-31 16:31 by 690616278
[考研] 08工科,295,接受跨专业调剂 +4 lmnlzy 2026-03-31 4/200 2026-03-31 15:25 by 西京学院招办
[考研] 270求调剂 +3 小杰pp 2026-03-31 4/200 2026-03-31 12:59 by wxiongid
[考博] 材料专业申博 +5 杜雨婷dyt 2026-03-29 5/250 2026-03-31 11:19 by oooqiao
[考研] 08工科,295,接受跨专业调剂 +6 lmnlzy 2026-03-30 6/300 2026-03-31 10:04 by cal0306
[考研] 085602化工求调剂(331分) +8 111@127 2026-03-30 8/400 2026-03-30 21:23 by 研究僧导导
[考研] 329求调剂 +10 钮恩雪 2026-03-25 10/500 2026-03-29 13:32 by peike
[有机交流] 高温高压反应求助 10+4 chibby 2026-03-25 4/200 2026-03-27 21:08 by BT20230424
[考研] 341求调剂 +7 青柠檬1 2026-03-26 7/350 2026-03-27 00:19 by wxiongid
信息提示
请填处理意见