计算化学公社

标题: 求助：动力学模拟得到CO2密度与实际不一致 [打印本页]

作者
Author: 姜酱将 时间: 2023-10-27 19:24
标题: 求助：动力学模拟得到CO2密度与实际不一致
在近期工作中遇到一个问题：经过refprop的计算，在418K，30MPa条件下，CO2的密度为0.50605g/cm3。但是在MS中，我初始利用ac模块建立了密度为0.1g/cm3，想着用NPT给它压到0.50605，但是最后的结果才大概0.233左右，时间是200ps，NPT系综。请问这个是怎么回事？
(, 下载次数 Times of downloads: 22)

作者
Author: Qing-958 时间: 2024-6-24 15:03
请问你解决了吗？

作者
Author: 姜酱将 时间: 2024-6-24 16:16

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

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

作者
Author: Qing-958 时间: 2024-6-24 16:25

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

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

作者
Author: chenzhe 时间: 2024-8-19 08:19

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程序应该是默认是水的压缩系数，所以也得修改

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

作者
Author: Qing-958 时间: 2024-9-9 14:40
本帖最后由 Qing-958 于 2024-9-9 14:42 编辑

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

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

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