我目前在用lammps例子中给的green-kubo方法在求粘度，我先在NVT下驰豫了50万步（500ps），然后跑NVE下的100万步，求粘度，跟例子中的设置是一样的，可是我发现在NVE计算过程中，温度一直在升高，比如我想求4484K的粘度，NVT驰豫完后，体系温度大约在4500K左右，是合理的，可是NVE下计算完成后，温度就升到了5000K，老师，我想问一下，我这样求出的粘度还准确吗？如果这种做法不对的话，我应该怎么办？（我也尝试了在NPT和NVT下求粘度，可是都出错了） 返回小木虫查看更多