计算化学公社

标题: 用VASP5.4.4 DFPT的方法计算声子谱出现错误 [打印本页]

作者
Author: ShuaiCui 时间: 2021-11-20 20:40
标题: 用VASP5.4.4 DFPT的方法计算声子谱出现错误
尊敬的各位老师您好，我再用VASP5.4.4版本计算声子谱的时候出现了这样的错误（见图片），请问这样的错误是什么原因导致的呢？谢谢

作者
Author: 卡开发发 时间: 2021-11-20 21:05
VASP应该没做DFPT部分对IVDW的支持，你要不换有限差分做（IBRION=5或6）。

作者
Author: ShuaiCui 时间: 2021-11-20 21:13

卡开发发发表于 2021-11-20 21:05
VASP应该没做DFPT部分对IVDW的支持，你要不换有限差分做（IBRION=5或6）。

那我还需要重新结构优化嘛？

作者
Author: 卡开发发 时间: 2021-11-21 02:33

ShuaiCui 发表于 2021-11-20 21:13
那我还需要重新结构优化嘛？

声子谱计算与优化过程的计算方法和精度一致即可。

作者
Author: sobereva 时间: 2021-11-21 05:31
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并清楚、准确反映出帖子具体内容，避免有任何歧义，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题“用DFPT的方法计算声子谱”改了，以后务必注意

作者
Author: ShuaiCui 时间: 2021-11-21 08:06
谢谢老师

作者
Author: 松鼠 时间: 2024-10-24 15:21
本帖最后由松鼠于 2024-10-24 15:30 编辑

请问一下有限差分法用IBRION=-1和用IBRION=5或6具体的区别是啥啊？

作者
Author: 松鼠 时间: 2024-10-24 15:37

卡开发发发表于 2021-11-20 21:05
VASP应该没做DFPT部分对IVDW的支持，你要不换有限差分做（IBRION=5或6）。

请问一下有限差分法用IBRION=-1和用IBRION=5或6具体的区别是啥吗？

作者
Author: 卡开发发 时间: 2024-10-24 17:46

松鼠发表于 2024-10-24 15:37
请问一下有限差分法用IBRION=-1和用IBRION=5或6具体的区别是啥吗？

IBRION=-1本身没办法算有限差分，只是通过外部程序例如Pymatgen或者ASE产生一批扩超晶格原子位移的结构，然后通过IBRION=-1+NSW=0进行单点计算获得力，然后再通过后处理得到Hessian矩阵。IBRION=5其实相当于上面的过程全部在VASP内部进行，IBRION=6用DFPT替代有限差分。

作者
Author: 松鼠 时间: 2024-10-25 21:29

卡开发发发表于 2024-10-24 17:46
IBRION=-1本身没办法算有限差分，只是通过外部程序例如Pymatgen或者ASE产生一批扩超晶格原子位移的结构， ...

好的，感谢！

作者
Author: going6667 时间: 2024-12-26 23:08

卡开发发发表于 2024-10-24 17:46
IBRION=-1本身没办法算有限差分，只是通过外部程序例如Pymatgen或者ASE产生一批扩超晶格原子位移的结构， ...

卡卡老师，您好。我是一名刚刚接触声子色散计算的学生，我从很多途径中了解到声子色散计算需要构建超胞，从vaspwiki了解到声子色散需要构建超胞的原因是：半导体和绝缘体的电子屏蔽作用比较弱，需要考虑长程相互作用对力常数矩阵的影响。目前我有几个比较简单的问题一直想不出答案，能麻烦您指导一下吗：
(1)声子色散计算需要构建超胞的原因只是因为vaspwiki上说的这个原因吗？那么金属体系是不是不需要进行超胞？
(2)vasp+phonopy计算动力学矩阵有有限位移法和DFPT两种方法，我了解到有限位移法要求声子波矢与原胞边界supersize正交，而DFPT原则上不需要扩胞，但因为vasp原因，计算DFPT还是需要扩胞，请问我了解的地方有错误的地方吗？
(3)上述两种方法扩胞大小有经验条件吗？我从一些DFT数据库网站上了解到声子谱计算扩胞后原则上不应小于12埃，到20埃比较保险。是不是还得去测试超胞大小直到声子谱变化收敛？

作者
Author: 卡开发发 时间: 2024-12-27 00:45

going6667 发表于 2024-12-26 23:08
卡卡老师，您好。我是一名刚刚接触声子色散计算的学生，我从很多途径中了解到声子色散计算需要构建超胞， ...

1、金属仍然可能需要，甚至可以说对于小晶胞的情况绝大多数情况都得扩。只是有可能不见得要那么大的超胞，其他情况有可能需要的超胞更大，按到也需要通过足够的测试，不过因为声子的计算量太大大家没人愿意这样做。除此之外，即便对DFPT来说VASP也没办法处理声子BZ采样q点为Gamma以上的情况，得通过超晶胞的声子谱进行折叠。
2、有限位移没有这方面限制吧，但是确实需要超晶胞上面已经提到，因为原子位移导致的密度变化也可以看成是一种响应。DFPT原则上可以不扩胞，这方面大概是VASP自身算法的原因，具体如何实现我没有深究，如果有人能讨论得清楚可以进行补充。
3、严格说得测试，科学的角度说主观经验不太可取，但碍于算不动，没人愿意这么做。

作者
Author: going6667 时间: 2024-12-27 10:07

卡开发发发表于 2024-12-27 00:45
1、金属仍然可能需要，甚至可以说对于小晶胞的情况绝大多数情况都得扩。只是有可能不见得要那么大的超胞 ...

感谢卡卡老师百忙之中抽空回答！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)