计算化学公社

标题: 求助：利用phonopy生成声子态密度时报错 [打印本页]

作者
Author: XYL123 时间: 2022-11-7 15:05
标题: 求助：利用phonopy生成声子态密度时报错
对结构进行精优化之后，利用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文件(声子态密度)

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

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

(, 下载次数 Times of downloads: 24)

作者
Author: dockery 时间: 2022-11-7 22:38
(, 下载次数 Times of downloads: 21)

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

作者
Author: XYL123 时间: 2022-11-8 08:41

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

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

作者
Author: get-it 时间: 2022-11-9 04:29
检查一下原胞里面有多少个原子

作者
Author: XYL123 时间: 2022-11-9 08:58

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

原胞里面只有五个原子

作者
Author: get-it 时间: 2022-11-9 17:35

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

所以你不能设6

作者
Author: XYL123 时间: 2022-11-9 18:46

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

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

作者
Author: get-it 时间: 2022-11-12 18:58

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

改你的mesh.conf

作者
Author: 乐平 时间: 2023-8-12 13:05

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

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

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

作者
Author: get-it 时间: 2023-8-16 00:58

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

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

作者
Author: 乐平 时间: 2023-8-16 09:57

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

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

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