计算化学公社

标题: 求助，lammps中NVT系综下求解热容的问题 [打印本页]

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

作者
Author: Computistry 时间: 2023-3-27 01:04
楼主后来解决了吗？我也遇到这个问题了。。非常感谢！！！

作者
Author: 十五楼的钢板 时间: 2023-5-4 20:25
楼主和楼上后来解决这个问题了么？我用计算等压热容时也遇到这个问题了。。。测试了二氧化碳和水两种物质的热容与实验值差距非常大，可以交流一下么，非常感谢！！！

作者
Author: kimwyy 时间: 2024-6-17 12:09
nvt系综计算热容也是需要升温改变温度的，你是完全在恒温下nvt。

作者
Author: haitaier 时间: 2024-6-17 19:24

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

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

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