24小时热门版块排行榜    

查看: 5297  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 20:58 by Equinoxhua
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +3 zju2000 2026-09-26 3/150 2026-09-27 16:26 by lala1008
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:26 by 4IW0sJvtEqX8
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:21 by 4IW0sJvtEqX8
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:03 by 4IW0sJvtEqX8
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 oeyPlfmMUpOC 2026-09-26 3/150 2026-09-27 09:13 by 4IW0sJvtEqX8
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:44 by XgCC8uTcwILl
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:39 by XgCC8uTcwILl
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 GxESXJmptOtQ 2026-09-25 3/150 2026-09-27 06:32 by XgCC8uTcwILl
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 06:08 by XgCC8uTcwILl
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 06:05 by XgCC8uTcwILl
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 05:56 by XgCC8uTcwILl
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:56 by XgCC8uTcwILl
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 vQZoDrrm7mUF 2026-09-25 5/250 2026-09-27 05:53 by XgCC8uTcwILl
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:50 by XgCC8uTcwILl
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:44 by XgCC8uTcwILl
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 5/250 2026-09-27 05:44 by XgCC8uTcwILl
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 6/300 2026-09-27 05:42 by XgCC8uTcwILl
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 PGTSwSC3F6nQ 2026-09-24 5/250 2026-09-25 11:19 by yhq9807
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 PGTSwSC3F6nQ 2026-09-24 7/350 2026-09-25 11:18 by yhq9807
信息提示
请填处理意见