lammps中如何用pimd计算比热容

想问一下:
lammps中用路径积分函数(pimd)计算比热容有人会么?
本人主要研究MOFS材料的分子动力学的,查看文献发现pimd计算比热容,相较于md计算的可靠性要更高一些,但是不知道具体的实现方式。
想问一下有没有朋友会,能够交流分享一下。