MS计算加载轴压条件下吸附出现问题

hpuCA

         使用wiser模型构建煤的分子结构，并尝试计算加不同轴相压力条件下吸附甲烷量的变化情况，我的计算过程是首先绘制wiser分子模型，然后结构优化，然后退火处理（Dynamics 设置为NVT），然后再进行退火处理（Dynamics设置NPT 并使用S-M加压函数设置Z方向压力），最后在返回结果中挑选里面能量最小的一帧进行吸附计算。可是在使用设置不同压力条件下优化后的模型吸附甲烷时候，得到的吸附量结果没有规律且偏离实际结果较远。请问各位大神，这样的计算流程是否正确？模拟结果很差，中间问题大概出现在哪里？谢谢各位！