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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP计算稳定体系的零点能为什么也会出现虚频？
51/1返回列表
查看: 3774  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

cindywei

木虫 (初入文坛)

[求助] VASP计算稳定体系的零点能为什么也会出现虚频？

最近用VASP计算Gamma点的零点能做校正，原本热力学及力学上都能稳定存在的钯氢化合物有些成分点出现了虚频，是哪里出现问题了吗？
INCAR设置
ISTART = 0
ICHARG = 2
ENCUT = 450
ISMEAR = 0
PREC = A
EDIFF = 1E-6
EDIFFG = -0.001
NSW = 1
IBRION = 5
POTIM = 0.015
NFREE = 2

OUTCAR
1 f  =   20.970164 THz   131.759427 2PiTHz  699.489360 cm-1    86.725659 meV
2 f/i=   15.743741 THz    98.920841 2PiTHz  525.154652 cm-1    65.110902 meV
3 f/i=   15.775374 THz    99.119597 2PiTHz  526.209815 cm-1    65.241726 meV
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-01-02 17:15:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘小伟110

木虫 (小有名气)
送红花一朵
引用回帖:
4楼: Originally posted by cindywei at 2014-01-03 08:19:33
计算零点能的方法有很多种吧，如果你只考虑轻质的H及O的零点能的话，可以在POSCAR选择H及O不固定而Cu固定，这种计算不完全等同于晶格优化的弛豫，是每个不固定原子按特定几个方向有限的弛豫来获得其零点振动能。我 ...

cindywei您好，首先感谢您的回复，我有所理解，我计算了个水的频率，存在很多f/i，照理说，单个分子气体频率计算，应该没有f/i的，而且按E=1/2hv公式，把所有频率的能量相加，所得能量比较大，零点能的话有可能比本身能量还大，还有就是零点能校正，它是和谁校正啊，单独加上个这么大的数，岂不是很糟糕。
还望您给解释一下。。。。
VASP计算稳定体系的零点能为什么也会出现虚频？
3.3.jpg
赞一下 回复此楼
慢慢的思考。。。
5楼2014-01-03 21:55:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

tider

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
franch: 金币+2, 谢谢回帖交流，， 2013-04-22 20:41:32
cindywei: 金币+5, ★★★很有帮助 2013-04-27 11:44:45
http://wolverton.northwestern.ed ... phonon-calculations

One drawback of the frozen-phonon method is that large supercells are needed to accurately calculate the force constant matrix. The periodic boundary conditions used in DFT calculations cause problems for the frozen phonon method. Displacing one atom in a small unit cell creates forces on all the atoms in the same unit cell, but also on the periodic images of these atoms.
赞一下 回复此楼
2楼2013-04-18 10:59:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘小伟110

木虫 (小有名气)
cindywei您好，我想问您一个问题就是我要做零点能校正的话，按E=1/2hv公式，把所有频率的能量相加，大约就是0.7eV，按理说，零点能校正的话，不会是这么大的值，也就是0.02eV，是不是这中计算方法有误呢？
我如果要对气态分子吸附态进行零点能校正的话，比如H2O在Cu上的吸附，要计算频率，我是不是应该把H2O和Cu的第一层原子弛豫，把所有的频率的能量相加，就算是零点能呢？是不是应该把Cu原子弛豫的频率去掉，只加H2O的频率呢，但是那样的话，我根本就不能从所有的频率中区分H2O和Cu的频率啊？
最近甚是苦恼啊，还望您能指点迷津，期待您的回复，在此，万分感谢。。。
赞一下 回复此楼
慢慢的思考。。。
3楼2014-01-02 16:16:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cindywei

木虫 (初入文坛)
引用回帖:
3楼: Originally posted by 刘小伟110 at 2014-01-02 16:16:15
cindywei您好，我想问您一个问题就是我要做零点能校正的话，按E=1/2hv公式，把所有频率的能量相加，大约就是0.7eV，按理说，零点能校正的话，不会是这么大的值，也就是0.02eV，是不是这中计算方法有误呢？
我如果要 ...

计算零点能的方法有很多种吧，如果你只考虑轻质的H及O的零点能的话，可以在POSCAR选择H及O不固定而Cu固定，这种计算不完全等同于晶格优化的弛豫，是每个不固定原子按特定几个方向有限的弛豫来获得其零点振动能。我跟个INCAR的设置，你可以通过参数来理解一下，希望对你有用
ISTART = 0
ICHARG = 2
ENCUT = 375
ISMEAR = 0
PREC = A
EDIFF = 1E-7
EDIFFG = -0.001
NSW = 1
IBRION = 5
POTIM = 0.015
NFREE = 2
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

刘小伟110
4楼2014-01-03 08:19:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 03:55 by OEbVnUOu01ol
[论文投稿] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:43 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 03:40 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 00:55 by OEbVnUOu01ol
[基金申请] 2026国自然函评费到账 +11 羊腰板 2026-08-21 11/550 2026-08-22 22:03 by lsbin3733
[基金申请] 什么时候开奖？ +7 CrisMessi 2026-08-18 8/400 2026-08-22 18:29 by 淀粉搬运工
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[教师之家] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 时间戳变了，能看出什么问题？ +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
关闭
请填处理意见
关闭