计算化学公社

标题: EPM刚体模型计算CO2粘度，与柔性模型的Greenkubo计算方法一致，但计算结果相差很大 [打印本页]

作者
Author: CO2CCC 时间: 2025-12-28 13:52
标题: EPM刚体模型计算CO2粘度，与柔性模型的Greenkubo计算方法一致，但计算结果相差很大
求助，利用EPM刚体模型在lammps软件中计算CO2粘度时，与柔性模型的Greenkubo计算方法一致，但计算结果相差很大，下图贴有输入脚本，请教我的刚体模型设置是否合理，fix指令有何错误，CO2分子模板文件是否错误
(, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0)

作者
Author: yuzc 时间: 2025-12-28 20:02
先预言一下，你模拟的值可能是偏高很多。
1. NVE下Green-Kubo理论上是最干净的，但是对温度偏移敏感，尤其粘度还是温度的函数。所以可以用一个弱耦合的NVT下采样。调高一下timecon的值。
2. trap(f_SS[i])那块，如果你不控制积分上下限，虽然原理上在噪声区统计值正负抵消，但是更常见的情况是只会随着时间而越来越大。建议你先找到关联函数衰减为0并且刚开始振荡的t，从0积分到这个t，再看下结果。
3. 具体问题具体分析了，听听其他同行的建议。

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