计算化学公社

标题: 求助：模拟气液界面系综的选择 [打印本页]

作者
Author: Lisy 时间: 2025-12-5 19:35
标题: 求助：模拟气液界面系综的选择
我的初始体系是一个气液界面，我想模拟气体分子逐渐扩散到液相中的过程，请问我应该选择NVT系综还是NPT系综来模拟比较好呢

作者
Author: sobereva 时间: 2025-12-5 20:17
通常NVT，恰当设置控压时NPT亦可

作者
Author: 牧生 时间: 2025-12-5 20:18
得用npt

nvt的体积是不变，不符合你的需求。

作者
Author: Donalaren 时间: 2026-4-11 22:48
本帖最后由 Donalaren 于 2026-4-11 22:49 编辑

牧生发表于 2025-12-5 20:18
得用npt

nvt的体积是不变，不符合你的需求。

请问大佬，我也想做气体（CO2分子）在水中的溶解。
在用npt系综模拟气液界面时，用半各项同性控压，x/y压缩系数为0，z压缩系数4.5e-5，控压5bar，控温293K。
这样跑出来的气体密度非常大，已经超过了真实情况下的液态CO2密度。
(TRAPPE力场和EPM2力场都试过)

作者
Author: 牧生 时间: 2026-4-12 08:17

Donalaren 发表于 2026-4-11 22:48
请问大佬，我也想做气体（CO2分子）在水中的溶解。
在用npt系综模拟气液界面时，用半各项同性控压，x/y ...

http://bbs.keinsci.com/thread-33597-1-1.html

作者
Author: Donalaren 时间: 2026-4-12 17:32
本帖最后由 Donalaren 于 2026-4-12 17:38 编辑

牧生发表于 2026-4-12 08:17
http://bbs.keinsci.com/thread-33597-1-1.html

谢谢您的回复和分享，TRAPPE力场和EPM2力场我都试过，控压后的密度都很大。
在我观察看来，NPT系综半各项同性控压模拟气液界面，气相会被一直压缩，密度压缩到接近800kg/m3，这已经是液态CO2的密度了。而我的z轴控压才5bar。
不知道是否是NPT控压不能用于高压下气-液界面？
D:\下载\471b6126-1840-4580-a33d-13505274cff6.jfif

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