利用Gromacs模拟CO2与辛烷混相问题

ocean6

最近研究了一个CO2和辛烷的混相问题，发现计算的体系压力结果总是偏大，参考着相关的文献设置，辛烷力场用的是OPLASS AA（电荷用的是1.2*CM，之前看的社长发的帖子，比较推荐的)，CO2力场用的是TrappE Small力场（五点的，也是公社之前发布的），这两个力场分别单独测试了辛烷和CO2在温压下的的密度，结果误差很小；截断半径为1.4nm，NVT系综，初始模拟见下图，具体的模拟设置见附件
对于很多文献中所说，气体分子数量控制体系的压力这一说法我不太明白，我自己想的是根据加入的CO2的数量和区域体积，算出密度，再结合温度值，反算出压力；根据这个思路，我在20*5*5nm3的一半体积的盒子中加入了1045个CO2分子，温度是333K，计算出来的CO2密度应该是305.47kg/m^3，压力应该是10.23MPa左右，发现没有混相（见下图所示，辛烷的二维密度分布图），要加入2000个左右CO2才差不多混相，这远远大于已报道的最小混相压力了，不知道哪里有问题，有没有大佬帮忙解答一下