计算化学公社

标题: 求助 nvt模拟co2在生物炭上的吸附解吸，计算压力为负值 [打印本页]

作者
Author: 茉莉花 时间: 2025-7-15 20:58
标题: 求助 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)

欢迎光临计算化学公社 (http://bbs.keinsci.com/)