求助：利用phonopy生成声子态密度时报错

XYL123 · 发表于 Post on 2022-11-7 15:05:35

对结构进行精优化之后，利用DFPT微扰法，计算声子态密度时有问题。
1.cp CONTCAR POSCAR

2.phonopy -d --dim="1 1 1"

3.cp SPOSCAR POSCAR 然后进行一次vasp静态计算，得到vasprun.xml

4.phonopy --fc vasprun.xml，生成FORCE_CONSTANTS力常数矩阵文件

5.phonopy -p -s mesh.conf mesh.conf需要提前准备mesh.conf 输出partial_dos.pdf文件(声子态密度)

执行完第五步之后会出现这样的报错。求助怎么解决？（扩胞的方法试过了，也重新精优化过几次。）

有的结构可以，有的结构不行

dockery · 发表于 Post on 2022-11-7 22:38:48

mesh.conf 的原子数是否有问题

XYL123 · 发表于 Post on 2022-11-8 08:41:04

dockery 发表于 2022-11-7 22:38
mesh.conf 的原子数是否有问题

mesh.conf里，因为没扩胞dim设置为111，K点取样密度MP设置为30 30 30
试过扩胞，修改dim也不行

get-it · 发表于 Post on 2022-11-9 04:29:33

检查一下原胞里面有多少个原子

XYL123 · 发表于 Post on 2022-11-9 08:58:44

get-it 发表于 2022-11-9 04:29
检查一下原胞里面有多少个原子

原胞里面只有五个原子

get-it · 发表于 Post on 2022-11-9 17:35:27

XYL123 发表于 2022-11-9 08:58
原胞里面只有五个原子

所以你不能设6

XYL123 · 发表于 Post on 2022-11-9 18:46:38

get-it 发表于 2022-11-9 17:35
所以你不能设6

那请问需要修改哪个文件，谢谢！

get-it · 发表于 Post on 2022-11-12 18:58:34

XYL123 发表于 2022-11-9 18:46
那请问需要修改哪个文件，谢谢！

改你的mesh.conf

乐平 · 发表于 Post on 2023-8-12 13:05:34

get-it 发表于 2022-11-9 04:29
检查一下原胞里面有多少个原子

看到您回答了很多关于声子计算的问题，请问 phonopy 计算声子谱的 band.config，声子态密度的 mesh.config 文件是如何生成的？

我看了 phonopy 的官网，似乎没找到如何生成这些文件的方法，只是罗列了一些参数解释，感觉很支离破碎……

get-it · 发表于 Post on 2023-8-16 00:58:08

乐平发表于 2023-8-12 13:05
看到您回答了很多关于声子计算的问题，请问 phonopy 计算声子谱的 band.config，声子态密度的 mesh.confi ...

一般从phonopy的例子中得到，再自己改改。
phoonpy的手册里该有的都有，但不适合入门

乐平 · 发表于 Post on 2023-8-16 09:57:08

get-it 发表于 2023-8-16 00:58
一般从phonopy的例子中得到，再自己改改。
phoonpy的手册里该有的都有，但不适合入门

谢谢回复。
的确 phonopy 手册里的参数很权威，但是太零散，对初学者不友好……

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[VASP] 求助：利用phonopy生成声子态密度时报错