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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP接FROPHO计算H2时声子谱有虚频怎么办
111/2返回列表12下一页
查看: 2196  |  回复: 20
查看全部回帖@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

liu108

木虫 (正式写手)

[求助] VASP接FROPHO计算H2时声子谱有虚频怎么办

我是用VASP优化氢气后接FROPHO来计算H2的声子谱，
VASP优化中的INCAR设置为：
SYSTEM = H2
ENCUT = 550
ISTART = 0
ICHARG = 2
PREC = High
EDIFFG = 1E-5
EDIFFG = -0.001
ISMEAR = 0
GGA = 91
NSW = 100
IBRION = 2
ISIF = 2
ISYM = 0
赝势用的PW91， K点为4*4*4，盒子长度为10Å。
FROPHO计算中用的1*1*1的supercell，计算的声子谱（PDOS）存在虚频
2*2*2的supercell试过，也有虚频。望各位大仙指点。

H2的PDOS
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
会当凌绝顶，一览众山小。
1楼2012-04-12 09:12:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liu108

木虫 (正式写手)
引用回帖:
2楼: Originally posted by rxd337 at 2012-04-12 11:43:16:
ISIF=2是只弛豫离子，你check过ISIF=3吗？

ISIF=3明显不对啊，氢气不用优化格子的啊。
赞一下 回复此楼
会当凌绝顶，一览众山小。
3楼2012-04-12 14:10:34
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liu108

木虫 (正式写手)
引用回帖:
4楼: Originally posted by rxd337 at 2012-04-12 15:23:46:
那就很难说了，格子也是一种自由度，不优化格子怎样保证结构是在最稳定状态。

氢气不是晶体啊，不能按照优化晶体的方法做。
赞一下 回复此楼
会当凌绝顶，一览众山小。
5楼2012-04-12 16:07:28
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liu108

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

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

对氢气的INCAR要怎样设置呢？还有就是K点该怎样设置？你说的不排除人为引入的对称是什么意思呢？谢谢啦。
赞一下 回复此楼
会当凌绝顶，一览众山小。
9楼2012-04-12 19:44:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liu108

木虫 (正式写手)
引用回帖:
7楼: Originally posted by rxd337 at 2012-04-12 16:59:52:
氢气在0k下应该也是晶体

氢气的计算不能当做晶体来计算啊
赞一下(1人) 回复此楼
会当凌绝顶，一览众山小。
10楼2012-04-12 19:45:58
 已阅   回复此楼   关注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的回帖

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

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

liu108

木虫 (正式写手)
送鲜花一朵
引用回帖:
20楼: Originally posted by 贺仪 at 2012-04-20 10:07:43:
如果不大的话，就没什么问题吧，我没做过声子，对于虚频的消除，论坛里面有很多讨论的，你看搜下能不能解决你的问题

我还有一个问题，就是在这三个小虚频存在的情况下，算出来氢气的熵值很奇怪，从10 K到1000 K的熵值变化为102--165， 0 K时为零，我是每隔10 K记一次数值。
赞一下 回复此楼
会当凌绝顶，一览众山小。
21楼2012-04-20 10:17:18
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liu108 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭