求助 nvt模拟co2在生物炭上的吸附解吸 ，计算压力为负值

茉莉花

求助各位老师 我是用LAMMPS软件小白 ，想要计算co2 ch4 饱和蒸汽压力计算吸附等温线 ，下面是我的一些问题：
1、以下脚本计算压力用ovito处理数据 但算出来的co2压力是负值怎么办？2、我觉得饱和蒸汽压力是在最后吸附平衡时整个区域co2的压力是否有误？
compute         1 all pe
compute         2 mobile_atoms ke/atom
compute         3 mobile_atoms cluster/atom 5.0
compute         peratom mobile_atoms stress/atom NULL
compute         p mobile_atoms reduce sum c_peratom[1] c_peratom[2] c_peratom[3]
variable        press equal -(c_p[1]+c_p[2]+c_p[3])/(55*55*40)