24小时热门版块排行榜    

CyRhmU.jpeg
查看: 4508  |  回复: 23
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

gemucai

木虫 (正式写手)

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

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




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


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



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

刘小伟110

木虫 (小有名气)

引用回帖:
23楼: Originally posted by gemucai at 2014-01-02 17:45:56
至少在QE中,我们还暂时没有用到零点能校正这个方法。另外,你说的INCAR,是VASP里面的输入文件,但是VASP是不能直接算声子的,必须外接PHONON程序。

f/i是虚频吧?这虚频确实很大,零点能不会那么大的,建议检 ...

gemucai好,首先感谢您的回复,我有所理解,我计算了个水的频率,存在很多f/i,照理说,单个分子气体频率计算,应该没有f/i的,而且按E=1/2hv公式,把所有频率的能量相加,所得能量比较大,零点能的话有可能比本身能量还大,还有就是零点能校正,它是和谁校正啊,单独加上个这么大的数,岂不是很糟糕。
还望您给解释一下。。。。我用这个公式是不是不对啊,用频率之后的能量相加是不是不对啊?
我如果要对气态分子吸附态进行零点能校正的话,比如H2O在Cu上的吸附,要计算频率,我是不是应该把H2O和Cu的第一层原子弛豫,把所有的频率的能量相加,就算是零点能呢?是不是应该把Cu原子弛豫的频率去掉,只加H2O的频率呢,但是那样的话,我根本就不能从所有的频率中区分H2O和Cu的频率啊?
您能不能给个INCAR呢?
最近甚是苦恼啊,还望您能指点迷津,期待您的回复,在此,万分感谢。。。
慢慢的思考。。。
24楼2014-01-03 22:31:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 24 个回答

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的回帖

gemucai

木虫 (正式写手)

问题解决了。

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

希望这对大家以后的计算起到帮助的作用,也不妨讨论下为何虚频对degauss如此敏感。
3楼2012-05-25 18:59:36
已阅   回复此楼   关注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的回帖
信息提示
请填处理意见