应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 哭了，各种虚频！
131/2返回列表12下一页
查看: 4592  |  回复: 23
查看全部回帖 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

gemucai

木虫 (正式写手)

[求助] 哭了，各种虚频！

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




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


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

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

VASP and MS condensed matter physics

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-05-24 09:26:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
问题解决了。

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

希望这对大家以后的计算起到帮助的作用，也不妨讨论下为何虚频对degauss如此敏感。
赞一下(2人) 回复此楼
3楼2012-05-25 18:59:36
已阅   回复此楼   关注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的回帖

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

木虫 (正式写手)
刚用occupation=fixed算了一遍，跟取比较大的smearing效果是一样一样滴，这是怎么回事呢？
赞一下 回复此楼
11楼2012-05-26 11:17:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
引用回帖:
12楼: Originally posted by 姚海峰456 at 2012-08-01 21:42:17
已经用第一性原理优化好的单胞，用gulp计算声子谱会出现虚频，怎么解决？——求大神帮忙啊！！！

我记得gulp是经典的紧束缚方法吧？朝出现虚频的方向改变原子位置试试。
赞一下 回复此楼
13楼2012-08-02 18:56:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
内容已删除
回复此楼
16楼2012-08-03 09:12:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
引用回帖:
14楼: Originally posted by souledge at 2012-08-02 20:39:39
GULP是紧束缚么……我怎么记得GULP是类似于经典的分子动力学的东西呢……记得实验室的人用这货跑上千个原子的分子动力学，也就一天的事～还是在一台三年前的四核电脑上……不过，至于他参数是不是足够收敛，我就不 ...

其实这不是真正的紧束缚，只是一种经典的紧束缚，类似于phonon里面的小位移方法。它通过算势能对位移的二阶偏导得出近邻次近邻相互作用的动力学矩阵，解久期方程得到振动频率。貌似这个势能确定了之后算起来确实相当快，我去年手算一简单体系也只花了一个下午。
赞一下 回复此楼
17楼2012-08-03 09:16:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gemucai 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭