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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » Vasp+phonopy 声子谱出现离奇虚频！
262/3返回列表上一页123下一页
查看: 4052  |  回复: 25
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

tinkie0411

银虫 (小有名气)
送红花一朵
引用回帖:
9楼: Originally posted by chuanghua304 at 2014-04-01 08:48:16
其实理论理论上两种晶胞做优化都可以，手册的例子里好像就是用的立方结构，但是要在band.conf里设置，primitive axial。用原胞也可以，那就不用设置primitive axial了。具体体系不同扩展的倍数也不同，你128个原子应 ...

是的 128原子是用立方基扩的 16x2x2x2=128 如果换成3x3x3的话16x3x3x3=432个原子 程序太大了 跑不起来 我也试过用primitive cell 的基矢扩成3x3x3 的4x3x3x3=108个原子 结果差别不大
我用的是PAW-GGA91 Ecut最大值是282 我取的500ev
好的 以后的计算我会把EDIFFG = -1.0e-06 和NPAR=4删除 我现在再用primitive cell扩成4x4x4的超胞 256个原子 已经跑的很慢很慢了 而且我
也不确定算这么长时间最后能否有个可以用的结果 很困惑
10楼dx0620说DFPT只能算Gamma点 你怎么看 还有你告诉我的IBRION=6的差分法 是不是优化完之后 再算超胞的时候加这个 这样的算法超胞计算只需要一次 我的理解对吗？非常感谢！
赞一下 回复此楼
11楼2014-04-01 16:15:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tinkie0411

银虫 (小有名气)
引用回帖:
10楼: Originally posted by dx0620 at 2014-04-01 09:19:29
我怎么感觉出现这样的虚频不是计算方法的问题啊，你用的什么赝势？如果出现这种问题，建议先不要用考虑半芯态赝势。你如果惯用胞的超胞只有128个原子的话，就用惯用胞吧，对原胞的设置如果不熟悉的话很容易出错。据 ...

感谢你的回复！我用的是PAW-GGA91，你是说是赝势的问题吗？那我试试PAW-PBE。惯用胞2x2x2的超胞是128个原子，但是我不确定我这种离奇的虚频是否是因为超胞不够大而引起的。我也试过用primitive cell的基矢扩胞扩成4x3x3x3=108个原子的，也有这种虚频的情况。我现在用的是直接位移法，也打算试试DFPT和6楼chuanghua说的差分法。我现在最大的问题就是不知道产生这种离奇虚频的原因，所以就只能不停地尝试各种方法但也不确定能否解决。
赞一下 回复此楼
12楼2014-04-01 16:32:17
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tinkie0411

银虫 (小有名气)
引用回帖:
10楼: Originally posted by dx0620 at 2014-04-01 09:19:29
我怎么感觉出现这样的虚频不是计算方法的问题啊，你用的什么赝势？如果出现这种问题，建议先不要用考虑半芯态赝势。你如果惯用胞的超胞只有128个原子的话，就用惯用胞吧，对原胞的设置如果不熟悉的话很容易出错。据 ...

我看phonopy manual 1.7.4手册上62页用DFPT方法画的NaCl声子谱与第4页直接位移法是一致的 是不是DFPT方法算的也不只是Gamma点了呢 期待您的回复 谢谢！
赞一下 回复此楼
13楼2014-04-01 18:00:34
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dx0620

木虫 (正式写手)
引用回帖:
12楼: Originally posted by tinkie0411 at 2014-04-01 16:32:17
感谢你的回复！我用的是PAW-GGA91，你是说是赝势的问题吗？那我试试PAW-PBE。惯用胞2x2x2的超胞是128个原子，但是我不确定我这种离奇的虚频是否是因为超胞不够大而引起的。我也试过用primitive cell的基矢扩胞扩成 ...

你就用惯用胞2x2x2的超胞就行，IBRION=-1。我还是觉得出现你这种虚频的话跟赝势有很大关系，我之前算过的就有这种情况，你可以看看，更奇怪。赝势你都先别用考虑半芯态的
Vasp+phonopy 声子谱出现离奇虚频！
赞一下 回复此楼

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

tinkie0411
14楼2014-04-01 19:48:35
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dx0620

木虫 (正式写手)
引用回帖:
13楼: Originally posted by tinkie0411 at 2014-04-01 18:00:34
我看phonopy manual 1.7.4手册上62页用DFPT方法画的NaCl声子谱与第4页直接位移法是一致的 是不是DFPT方法算的也不只是Gamma点了呢 期待您的回复 谢谢！...

NaCl对称性那么高，算的声子谱跟前面的一样很正常的，想出问题都难。1.7.4手册上第80页说的很清楚，只能计算Gamma点上的。
赞一下 回复此楼
15楼2014-04-01 19:55:03
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tinkie0411

银虫 (小有名气)
送红花一朵
引用回帖:
14楼: Originally posted by dx0620 at 2014-04-01 19:48:35
你就用惯用胞2x2x2的超胞就行，IBRION=-1。我还是觉得出现你这种虚频的话跟赝势有很大关系，我之前算过的就有这种情况，你可以看看，更奇怪。赝势你都先别用考虑半芯态的
...

非常感谢你的回复！你说的半芯态赝势是什么意思，是不是指赝势文件中带—sv —pv 和—d的赝势呢？
我找到的解释
有些轨道不属于价轨道，但是能量又不是太低（通常位于HOMO以下5-10eV的范围）。如果在post-HF计算中做冻结处理，会对分子结构和特性通常有不太大也不太小的影响（对有些量，影响是非常大的，比如四极矩）。这些轨道就是semi-core轨道，其中的电子称为semi-core电子。例如3d金属的3s3p电子，碘的4d电子，金的4f5s5p电子。
如果追求较高的精度，需要对这部分电子进行相关能处理，而不能冻结。当然对于HF计算，这么做是没有任何意义的。 是这样理解吗？
我现在做的是Ni2MnGa，画的是fcc奥氏体的，空间群225，按理说对称性也是相当好啊，不知道为什么会出现这种情况。POSCAR我在1楼贴出来了，POTCAR我用的是PAW-GGA91，Ni10Mn7Ga13 就是取的最简单的那种（Ga考虑了d），但是我也做过Ga3的，也出现这种情况，这应该没用你说的考虑半芯态的赝势吧？ 顺带一句，我也试过USPP-GGA91，结果一样。我现在打算试一下PAW-GGAPBE，如果还是有这种问题，我不知道还能换什么赝势了。另外我最近做了四方结构马氏体相的声子谱，就没有这种离奇虚频出现了。虽说奥氏体相低温不是稳定相，但是我现在做的是重复很多已经发表的文章，人家的都没有我这种情况，我到底错哪了呢？？？
赞一下 回复此楼
16楼2014-04-01 21:25:51
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tinkie0411

银虫 (小有名气)
引用回帖:
15楼: Originally posted by dx0620 at 2014-04-01 19:55:03
NaCl对称性那么高，算的声子谱跟前面的一样很正常的，想出问题都难。1.7.4手册上第80页说的很清楚，只能计算Gamma点上的。...

好的，看到了，多谢！
回复此楼
17楼2014-04-01 21:27:20
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dx0620

木虫 (正式写手)
引用回帖:
16楼: Originally posted by tinkie0411 at 2014-04-01 21:25:51
非常感谢你的回复！你说的半芯态赝势是什么意思，是不是指赝势文件中带—sv —pv 和—d的赝势呢？
我找到的解释
有些轨道不属于价轨道，但是能量又不是太低（通常位于HOMO以下5-10eV的范围）。如果在post-HF计算 ...

差不多这么理解，简单来说就是不单单只考虑最外侧电子。那你先试试PBE的赝势看看吧
赞一下 回复此楼
18楼2014-04-01 21:35:13
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tinkie0411

银虫 (小有名气)
引用回帖:
14楼: Originally posted by dx0620 at 2014-04-01 19:48:35
你就用惯用胞2x2x2的超胞就行，IBRION=-1。我还是觉得出现你这种虚频的话跟赝势有很大关系，我之前算过的就有这种情况，你可以看看，更奇怪。赝势你都先别用考虑半芯态的
...

这是文献中的声子谱，我画的是Gamma到X那段。奥氏体是不稳定，但是虚频我不是我画出来的那种。他们用的软件是PHONON software，我用的phonopy，虽说前者是收费的可能更好用，但是phonopy也发了很多好文章。我现在的情况总不会是因为软件的问题吧？很困惑！
Vasp+phonopy 声子谱出现离奇虚频！-1
QQ图片20140401153253.jpg
赞一下 回复此楼
19楼2014-04-01 21:35:16
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tinkie0411

银虫 (小有名气)
引用回帖:
18楼: Originally posted by dx0620 at 2014-04-01 21:35:13
差不多这么理解，简单来说就是不单单只考虑最外侧电子。那你先试试PBE的赝势看看吧...

回复的好快！实在抱歉跟你抱怨一大堆，这件事情真是让我困惑了好久！我会试试PBE的，有结果了我还会贴上来，衷心感谢！！！
赞一下 回复此楼
20楼2014-04-01 21:40:46
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tinkie0411 的主题更新
262/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭