求助，lammps中NVT系综下求解热容的问题

AngYi

我想利用lammps求解一种有机材料的定容比热容，利用以下公式。首先输出NVT系综下总能量变化，每五千步输出一次，一共输出一千万步，取后两百万步的数据做计算，但是计算数据出现了数量级的差距。这样计算的结果为1.89E6 J/mol/K, 这类材料的比热容应该在几百到一千的范围内，详细计算过程如下，希望各位老师指导一下。
计算公式：
lammps输出文件计算示例：
由于是real体系，将输出的总能量单位由kal/mol转化为J/每单位，再对能量进行平方，然后对温度、能量、能量的平方做时间平均。带入上面公式，求解出来为3.15973102904048E-18 J/K/单位，1 mol的物质其包含6.022E23个单位，将其相乘，得到最终的比热容为1.89E6 J/mol/K，明显偏于正常值。请问各位老师，我的求解过程是否合理呢？
此外，NVT过程中能量应该是收敛的。 ，最终能量在62500 Kal/mol +-500的范围内波动。

Computistry

楼主后来解决了吗？我也遇到这个问题了。。非常感谢！！！

十五楼的钢板

楼主和楼上后来解决这个问题了么？我用计算等压热容时也遇到这个问题了。。。测试了二氧化碳和水两种物质的热容与实验值差距非常大，可以交流一下么，非常感谢！！！

kimwyy

nvt系综计算热容也是需要升温改变温度的，你是完全在恒温下nvt。

haitaier

十五楼的钢板 发表于 2023-5-4 20:25
楼主和楼上后来解决这个问题了么？我用计算等压热容时也遇到这个问题了。。。测试了二氧化碳和水两种物质的 ...

请问你解决了这个问题了吗？我也遇到了，和实际的比热差了好几个数量级