计算化学公社

标题: phonopy进行热力学计算，在室温下F与E为什么相差那么大？ [打印本页]

作者
Author: 松鼠 时间: 2025-4-29 21:49
标题: phonopy进行热力学计算，在室温下F与E为什么相差那么大？
我想请问，其中熵值那么大，是因为计算的是超胞的熵值吗？dim=222.
其中亥姆赫兹自由能F与简谐近似下的声子能量E的关系，应该是F=E-TS，TS这么大明显不合理，问题出在哪呢？
(, 下载次数 Times of downloads: 19)

作者
Author: sobereva 时间: 2025-4-29 23:38
熵是按照J/mol/K给出的而不是kJ/mol/K给出的，看起来数量级大没什么异常

作者
Author: 松鼠 时间: 2025-4-30 08:36

sobereva 发表于 2025-4-29 23:38
熵是按照J/mol/K给出的而不是kJ/mol/K给出的，看起来数量级大没什么异常

可是F为-14.77kJ/mol和E为168.57kJ/mol，在不考虑非谐的情况下，200kJ/mol会不会有点大？

作者
Author: 松鼠 时间: 2025-4-30 08:42
而且free_energy不是deta F，我看过从0 K 到293 K的数据，是从83kJ/mol的零点能慢慢降低的。

作者
Author: sobereva 时间: 2025-4-30 14:42

松鼠发表于 2025-4-30 08:36
可是F为-14.77kJ/mol和E为168.57kJ/mol，在不考虑非谐的情况下，200kJ/mol会不会有点大？

不会

作者
Author: 松鼠 时间: 2025-4-30 15:03
好的，谢谢

作者
Author: 松鼠 时间: 2025-4-30 15:05

sobereva 发表于 2025-4-30 14:42
不会

是否也就是说是计算的2*2*2超胞的数据？

作者
Author: sobereva 时间: 2025-5-1 12:04

松鼠发表于 2025-4-30 15:05
是否也就是说是计算的2*2*2超胞的数据？

扩胞前的

并且倘若你给的是conventional cell，自动被phonopy约化成了更小的primitive cell，则给出的是primitive cell的热力学数据

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