Dynaphopy 计算高温声子谱

Shana · 发表于 Post on 2025-2-3 21:12:48

mly 发表于 2025-1-24 23:14
想请教您一下，我想计算300K下的吉布斯自由能，我需要在AIMD的INCAR中设置300K，在使用phonopy提取热力学性 ...

phonopy的300k自由能可以通过自由能的定义去算，即考虑T*S部分，S 是熵，和声子有关，phonopy可以根据声子谱计算出熵。当然，考虑到温度对声子的影响，phonopy支持QHA方法来计算有限温的自由能。建议看看phonopy的例子部分https://phonopy.github.io/phonopy/examples.html，以及理论部分https://phonopy.github.io/phonop ... perties-expressions。

mly · 发表于 Post on 2025-2-14 09:46:22

Shana 发表于 2025-2-3 21:12
phonopy的300k自由能可以通过自由能的定义去算，即考虑T*S部分，S 是熵，和声子有关，phonopy可以根据声 ...

感谢

高阁 · 发表于 Post on 2025-2-17 17:48:34

请问非谐性声子谱的数据怎么像phonopy一样导出呢？我需要用它们计算后续内容，感谢！还有我觉得这个声子谱似乎有些奇怪，请问它正常吗？这是2000k下的

Shana · 发表于 Post on 2025-2-22 12:02:43

高阁发表于 2025-2-17 17:48
请问非谐性声子谱的数据怎么像phonopy一样导出呢？我需要用它们计算后续内容，感谢！还有我觉得这个声子谱 ...

dynaphopy input OUTCAR -sfc FORCE_CONSTANTS 这个指令会给出修正后的力常数，然后用phonopy画出修正后的声子谱

高阁 · 发表于 Post on 2025-2-23 11:09:22

Shana 发表于 2025-2-22 12:02
dynaphopy input OUTCAR -sfc FORCE_CONSTANTS 这个指令会给出修正后的力常数，然后用phonopy画出修正 ...

谢谢

yuanzhe · 发表于 Post on 2025-2-24 21:11:53

本帖最后由 yuanzhe 于 2025-2-25 14:56 编辑

shana老师，用vasp跑aimd需要用phonopy扩胞的POSCAR进行计算吗？kpoints是直接vaspkit还是手动设置吗？

Shana · 发表于 Post on 2025-2-25 17:28:27

yuanzhe 发表于 2025-2-24 21:11
shana老师，用vasp跑aimd需要用phonopy扩胞的POSCAR进行计算吗？kpoints是直接vaspkit还是手动设置吗？

扩胞系数要和phonopy的系数一样，不然结果会有问题。晶胞比较大的话（扩胞系数比较大），用单k点就很好了。

yuanzhe · 发表于 Post on 2025-3-3 15:47:09

Shana 发表于 2025-2-25 17:28
扩胞系数要和phonopy的系数一样，不然结果会有问题。晶胞比较大的话（扩胞系数比较大），用单k点就很好了 ...

谢谢shana老师解惑，我一般扩2*2*2，用vaspkit生成的就可以了嘛

Shana · 发表于 Post on 2025-3-5 15:22:35

yuanzhe 发表于 2025-3-3 15:47
谢谢shana老师解惑，我一般扩2*2*2，用vaspkit生成的就可以了嘛

可以，phonopy也可以扩胞。phonopy扩完胞后，会有一个SPOSCAR，这个是未做位移的supercell。

yuanzhe · 发表于 Post on 2025-3-11 15:19:42

Shana 发表于 2025-3-5 15:22
可以，phonopy也可以扩胞。phonopy扩完胞后，会有一个SPOSCAR，这个是未做位移的supercell。

谢谢shana老师！！！

Autumn899 · 发表于 Post on 2025-4-8 22:56:06

请问什么命令可以导出能在origin里画图的声子数据文件呀

Shana · 发表于 Post on 2025-4-11 09:32:30

Autumn899 发表于 2025-4-8 22:56
请问什么命令可以导出能在origin里画图的声子数据文件呀

先保存重整后的力常数，然后用phonopy计算声子谱，导出dat数据

Autumn899 · 发表于 Post on 2025-4-11 17:01:24

Shana 发表于 2025-4-11 09:32
先保存重整后的力常数，然后用phonopy计算声子谱，导出dat数据

好的！

cdx · 发表于 Post on 2025-6-16 08:59:06

本帖最后由 cdx 于 2025-6-16 09:00 编辑

老师请问利用Dynaphopy可以得到声子寿命吗？已完成非谐计算，需要如何导出声子寿命数据，求指导

cdx · 发表于 Post on 2025-6-16 15:03:13

Shana 发表于 2025-2-22 12:02
dynaphopy input OUTCAR -sfc FORCE_CONSTANTS 这个指令会给出修正后的力常数，然后用phonopy画出修正 ...

老师你好，请问可以正常获得非谐修正后的力常数，可以利用phonopy导出声子寿命数据吗

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

[综合讨论] Dynaphopy 计算高温声子谱

浏览过的版块