EPM刚体模型计算CO2粘度，与柔性模型的Greenkubo计算方法一致，但计算结果相差很大

CO2CCC

求助，利用EPM刚体模型在lammps软件中计算CO2粘度时，与柔性模型的Greenkubo计算方法一致，但计算结果相差很大，下图贴有输入脚本，请教我的刚体模型设置是否合理，fix指令有何错误，CO2分子模板文件是否错误

yuzc

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