求助：动力学模拟得到CO2密度与实际不一致

姜酱将 · 发表于 Post on 2023-10-27 19:24:10

在近期工作中遇到一个问题：经过refprop的计算，在418K，30MPa条件下，CO2的密度为0.50605g/cm3。但是在MS中，我初始利用ac模块建立了密度为0.1g/cm3，想着用NPT给它压到0.50605，但是最后的结果才大概0.233左右，时间是200ps，NPT系综。请问这个是怎么回事？

Qing-958 · 发表于 Post on 2024-6-24 15:03:24

请问你解决了吗？

姜酱将 · 发表于 Post on 2024-6-24 16:16:05

Qing-958 发表于 2024-6-24 15:03
请问你解决了吗？

没，直接把正常密度输进去用的我

Qing-958 · 发表于 Post on 2024-6-24 16:25:08

姜酱将发表于 2024-6-24 16:16
没，直接把正常密度输进去用的我

好吧，我做的煤吸附CO2，分子密度比较大，输入正常密度会导致AC盒子无法建立。所以只能跑NPT，但是跑出来的密度相比实际值会低0.2左右。。。

chenzhe · 发表于 Post on 2024-8-19 08:19:51

Qing-958 发表于 2024-6-24 16:25
好吧，我做的煤吸附CO2，分子密度比较大，输入正常密度会导致AC盒子无法建立。所以只能跑NPT，但是跑出来 ...

初始模型密度低这么多，npt肯定无法正常控制到合理密度。npt要达到合理密度只能初始模型尽量接近，起码不能太离谱。从0.1控压到0.5的密度，盒子要缩到原来1/5这绝对不可能。另外无论多大密度ac模块都能建立，报错和这个无关。

至于实际能不能控制到你还得看MD力场和其他参数。对于forcite，估计八成无脑选了compass，compass算二氧化碳远不如专门的二氧化碳力场。还有就是压缩系数，液态和超临界态二氧化碳的压缩系数远比水大得多，而forcite或者其他MD程序应该是默认是水的压缩系数，所以也得修改

Huschein · 发表于 Post on 2024-8-21 01:18:08

盲猜用的默认的等温可压缩系数，密度问题90%的问题都出在这

Qing-958 · 发表于 Post on 2024-9-9 14:40:42

本帖最后由 Qing-958 于 2024-9-9 14:42 编辑

Huschein 发表于 2024-8-21 01:18
盲猜用的默认的等温可压缩系数，密度问题90%的问题都出在这

那我问一下大佬如何修改MS中默认的压缩系数，不清楚在哪里修改这个参数

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Material Studio] 求助：动力学模拟得到CO2密度与实际不一致

评分 Rate

浏览过的版块