求助进行CO2的模拟，C-rescale压浴导致密度异常

raku

各位老师好！

我在在进行CO2的模拟，使用的力场为TraPPE-Small，来自论坛。http://bbs.keinsci.com/forum.php ... 52&highlight=trappe
T=293.15 K； P= 60 bar. 能量最小化（EM）之后直接使用NPT计算密度，其中pcoupl = C-rescale。dt = 0.001; nstep = 5000000; tau-p = 5;Pcoupltype = isotropic.
CO2的临界温度压力为：304K，73.8 bar。测试压力比较靠近临界点。
结算结果发现在该温度压力条件下，计算的密度仅为200 kg/m3；而NIST数据为700+。改用pcoupl = Parrinello-Rahman，结果可以与NIST结果吻合。
C-rescale在其他远离临界点温压条件下结果与NIST差别不大。

请教下这种情况可能是什么原因呢？