pyHMA一个新的Python软件包，它与VASP接口，通过对NVTBorn-Oppenheimerabinitio分子动力学(AIMD)模拟的数据进行后处理来计算晶体系统的(经典)非谐特性。它以最近开发的谐波映射平均法(HMA)为基础，利用分析上已知的谐波行为来重新表述直接/传统的集合平均，以便在一定的CPU时间内显著提高精度。该软件包包括两个阶段：从vasprun.xml文件读取AIMD数据，然后计算非谐属性。第一阶段是与MD包相关的，而第二阶段是通用的，因为它接收的是所需格式的数据。尽管pyHMA目前与VASP接口计算HMA非谐能量和压力，但它的修改方式允许通过添加新的阅读器与其他代码(如LAMMPS)接口，并且一旦有相关数据，可以通过添加新方法计算其他HMA非谐波特性(如热容量)。
pyHMA官方文档：