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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP接FROPHO计算H2时声子谱有虚频怎么办
212/3返回列表上一页123下一页
查看: 3211  |  回复: 20
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

rxd337

木虫 (正式写手)

养殖场经理

【答案】应助回帖

★ ★
liliangfang: 金币+1, 谢谢交流 2012-04-13 08:39:40
liu108: 金币+1, ★有帮助 2012-05-08 11:53:05
你用的也是个晶胞啊，和一般的晶胞计算一样。只是把晶胞放大而已，但一定要放的足够大，你可以试一下ISIF=3,如果晶胞发生变化，那么就说明你的晶胞取的还不够大。
赞一下 回复此楼
都是假设！
11楼2012-04-12 20:01:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rex881026

木虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
liliangfang: 金币+2, 谢谢指教 2012-04-14 10:29:24
这应该是属于没有优化好。。可以cp CONTCAR POSCAR再多迭代两次优化。。
赞一下 回复此楼
12楼2012-04-13 09:56:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liu108

木虫 (正式写手)
引用回帖:
6楼: Originally posted by rxd337 at 2012-04-12 16:58:12:
你核查下吗，我只是觉得有这个可能。计算软件这东西，可调的参数就那几个，出现问题我们的方法就是挨个排查。虚频的存在就表示结构可能不稳定吗，只要保证结构最优了应该就不会有这个虚频，我觉得。

ISIF=3的试过了，还是有虚频啊。
赞一下 回复此楼
会当凌绝顶，一览众山小。
13楼2012-04-13 11:39:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liu108

木虫 (正式写手)
引用回帖:
8楼: Originally posted by 铁昂笑 at 2012-04-12 17:26:06:
个人建议：楼主计算的是氢气分子，无论在0K还是其他温度，都不需要考虑是晶体的情形。

但晶格参数的设置不要排除人为引入的对称，建议a,b,c取不同的数值，不要都取10埃。K点一定是只有GAMMA点。

你的建议我试过了，声子谱里有4个小的虚频值，怎么也消除不了，是什么个情况呢？
赞一下 回复此楼
会当凌绝顶，一览众山小。
14楼2012-04-13 11:40:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rxd337

木虫 (正式写手)

养殖场经理

【答案】应助回帖

★ ★
liliangfang: 金币+1, 多谢分享 2012-04-14 10:29:35
liu108: 金币+1 2012-05-08 11:53:18
氢气的体系可能是比较特别的系统，以为VASP一般都是处理晶体。这里有个文件是说这个的希望对你会有帮助：http://wenku.baidu.com/view/52414b8c680203d8ce2f24ae.html
赞一下 回复此楼
都是假设！
15楼2012-04-13 15:52:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★
感谢参与，应助指数 +1
liliangfang: 金币+2, 专家考核, 谢谢指教 2012-04-14 10:29:45
liu108: 金币+2 2012-05-08 11:53:29
话说，像H2这类的无需周期条件的分子计算声子，本来就是只计算Gamma点的就可以了的。所以用Phonopy接到VASP上，实在不如直接用VASP自带的DFPT计算来得省事啊~
赞一下 回复此楼
思想重于技巧，内涵重于表象
16楼2012-04-13 20:36:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liu108

木虫 (正式写手)
引用回帖:
16楼: Originally posted by souledge at 2012-04-13 20:36:20:
话说，像H2这类的无需周期条件的分子计算声子，本来就是只计算Gamma点的就可以了的。所以用Phonopy接到VASP上，实在不如直接用VASP自带的DFPT计算来得省事啊~

你说的VASP自带的DFPT没有用过呀，在哪里呢？
赞一下 回复此楼
会当凌绝顶，一览众山小。
17楼2012-04-16 09:32:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
liu108: 金币+2 2012-05-08 11:53:53
优化的时候isif=2就可以了，格子可以适当加大，另外你的k点太多了，也不知道你取到了gamma点没有。一般分子只用一个gamma点就够了，多了反而不准确
赞一下 回复此楼

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

liu108
嘟啦啦嘟啦啦嘟
18楼2012-04-18 20:59:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liu108

木虫 (正式写手)
送鲜花一朵
引用回帖:
18楼: Originally posted by 贺仪 at 2012-04-18 20:59:07:
优化的时候isif=2就可以了，格子可以适当加大，另外你的k点太多了，也不知道你取到了gamma点没有。一般分子只用一个gamma点就够了，多了反而不准确

对啊，我优化时isif=2，格子我用了10*11*12的，k点我也换成1*1*1的了，但是我用fropho算的PDOS总是在波长为-3，-2，-1时有3个小虚频。这个情况遇到过么，怎么解决呢？
赞一下 回复此楼
会当凌绝顶，一览众山小。
19楼2012-04-20 09:28:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)
引用回帖:
19楼: Originally posted by liu108 at 2012-04-20 09:28:59:
对啊，我优化时isif=2，格子我用了10*11*12的，k点我也换成1*1*1的了，但是我用fropho算的PDOS总是在波长为-3，-2，-1时有3个小虚频。这个情况遇到过么，怎么解决呢？

如果不大的话，就没什么问题吧，我没做过声子，对于虚频的消除，论坛里面有很多讨论的，你看搜下能不能解决你的问题
赞一下 回复此楼

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

liu108
嘟啦啦嘟啦啦嘟
20楼2012-04-20 10:07:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liu108 的主题更新
212/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 ELBZ6WwEaRj6 2026-11-05 4/200 2026-11-06 02:37 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-05 7/350 2026-11-06 02:35 by qL1cU46YIdK0
[基金申请] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 02:24 by qL1cU46YIdK0
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 02:24 by qL1cU46YIdK0
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-04 9/450 2026-11-06 02:15 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 ttH3f1oA6ace 2026-11-04 8/400 2026-11-06 02:02 by qL1cU46YIdK0
[公派出国] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 ttH3f1oA6ace 2026-11-04 13/650 2026-11-06 02:02 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +6 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-06 01:51 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 15/750 2026-11-06 01:51 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 10/500 2026-11-06 01:51 by qL1cU46YIdK0
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 01:42 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 11/550 2026-11-06 01:29 by qL1cU46YIdK0
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 01:29 by qL1cU46YIdK0
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-05 4/200 2026-11-05 23:39 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +5 ncTVBM6AOJHF 2026-11-04 13/650 2026-11-05 22:44 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +6 ncTVBM6AOJHF 2026-11-04 15/750 2026-11-05 22:44 by qL1cU46YIdK0
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-05 11:10 by PYZmcHr5kSaq
信息提示
关闭
请填处理意见
关闭