|
|
本帖最后由 五月雨 于 2022-5-17 16:24 编辑
1、跑了一个纯MEA的模拟体系,在能量最小化之后首先进行了4ns的NPT,平衡密度和文献一致,之后进行了10ns的NVT。力场使用的是Amber03,文献(da Silva)中使用的是Amber,具体指哪个没提。MEA的拓扑文件是根据GAFF产生的,和文献相同。但是不知道为什么,itp文件中各种力常数和文献有不小的差异。
最后得到平衡密度和文献一致,自扩散系数为0.41×10-5 cm2/s ,文献值(da Silva)为0.61×10-5 cm2/s,但是最后da Silva说考虑误差建议值在(0.61±0.3)×10-5 cm2/s。
请问各位老师,0.41这个值可不可以接受,来证明力场参数的准确性?
2、我在相同理论水平上进行了关于CO2的参数验证,模型选择的是EPM2,最终的得到的平衡密度为2.0kg/m3,实验值为1.624±0.10kg/m3;扩散系数为0.15 cm2/s ,实验为0.16 cm2/s。请问这个数据能不能作为力场参数的验证呢?
|
-
Moosavi.png
(35.97 KB, 下载次数 Times of downloads: 13)
-
Dhawal Shah.png
(117.72 KB, 下载次数 Times of downloads: 21)
-
-
box_MEA.top
132 Bytes, 下载次数 Times of downloads: 2
-
-
em.mdp
384 Bytes, 下载次数 Times of downloads: 1
-
-
npt.mdp
815 Bytes, 下载次数 Times of downloads: 4
-
-
nvt.mdp
723 Bytes, 下载次数 Times of downloads: 2
|