关于vasp+phonopy计算声子谱，光学支有虚频应如何解决的问题

876449830

各位老师，我用vasp+phonopy（密度泛函微扰理论）计算BaTiO4的声子谱，但是光学支有很大的虚频，而声学支却没有，这是为什么呢？应如何处理呢？谢谢各位的指点！

876449830

使用vasp继续优化后，虚频还剩三个，还是很难消除。更换QE计算后，显示虚频还有5个，不知如何操作了，这帖子莫沉呀，请论坛大佬指点一下呀

wzkchem5

你是想说BaTiO3吧？计算用的是什么空间群？消虚频的时候有没有降低空间对称性？

876449830

wzkchem5 发表于 2024-5-2 17:50
你是想说BaTiO3吧？计算用的是什么空间群？消虚频的时候有没有降低空间对称性？

谢谢您的回复。对的，就是钛酸钡，用的是Pm3̅m，我用的消虚频的方法是提高截断能优化，vasp计算时没考虑对称性的。

wzkchem5

876449830 发表于 2024-5-2 11:06
谢谢您的回复。对的，就是钛酸钡，用的是Pm3̅m，我用的消虚频的方法是提高截断能优化，vasp计算时没 ...

钛酸钡的空间群是P4mm（https://en.wikipedia.org/wiki/Barium_titanate），因为它是铁电体，低温下会发生自发的对称性破缺。
你用的空间群的对称性太高了，把结构往P4mm的方向微调一下，再优化。

876449830

wzkchem5 发表于 2024-5-2 18:45
钛酸钡的空间群是P4mm（https://en.wikipedia.org/wiki/Barium_titanate），因为它是铁电体，低温下会发 ...

啊，我的结构找错了嘛？

get-it

立方钛酸钡是高温相，简谐近似下就应该有虚频，需要考虑温度和非谐效应：
https://journals.aps.org/prb/abstract/10.1103/PhysRevB.103.094108

876449830

get-it 发表于 2024-5-2 19:07
立方钛酸钡是高温相，简谐近似下就应该有虚频，需要考虑温度和非谐效应：
https://journals.aps.org/prb/a ...

谢谢回复！学到了，我仔细看看这个工作。