计算化学公社

 找回密码 Forget password
 注册 Register

两相系统(有机分子/水)分子动力学模拟力场选择问题

查看数: 5532 | 评论数: 7 | 收藏 Add to favorites 1
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2021-11-19 16:00

正文摘要:

各位老师好,请教一下做有机分子(醇和中链脂肪酸)与水组成的两相系统(如图)的分子动力学模拟选择什么力场比较合适?目前使用的是amber99sb-ildn,通过acpype获取有机分子GROMACS拓扑文件,用RESP2(0.5)分配原子 ...

回复 Reply

xryao 发表于 Post on 2021-11-21 14:39:53
sobereva 发表于 2021-11-21 06:15
这没有体现任何没达到平衡的信息
要考察平衡没有,通常应当绘制粒子数分布、氢键数目、盒子尺寸等等随时 ...

好的 谢谢老师
sobereva 发表于 Post on 2021-11-21 06:15:27
xryao 发表于 2021-11-20 10:34
sob老师,请问如何判断两相系统是否收敛呢,或者说怎么判断该体系达到稳定状态了?对该两相系统进行了100 ...

这没有体现任何没达到平衡的信息
要考察平衡没有,通常应当绘制粒子数分布、氢键数目、盒子尺寸等等随时间的变化曲线,平坦了就说明收敛了
xryao 发表于 Post on 2021-11-20 09:10:11

好的,谢谢sob老师。
sobereva 发表于 Post on 2021-11-20 09:07:00
xryao 发表于 2021-11-20 09:04
好的,谢谢sob老师,每种分子的原子电荷加和是为0的,体系总电荷有0.02左右,这个是否可以忽略,我去尝试 ...

基本可以忽略
xryao 发表于 Post on 2021-11-20 09:04:23
sobereva 发表于 2021-11-20 08:07
脂肪酸用GAFF结合RESP2(0.5)
水用OPC,用比较流行的SPC/E也凑合。
当前不牵扯生物分子,根本不涉及AMBER ...

好的,谢谢sob老师,每种分子的原子电荷加和是为0的,体系总电荷有0.02左右,这个是否可以忽略,我去尝试一下您推荐的力场。
sobereva 发表于 Post on 2021-11-20 08:07:42
脂肪酸用GAFF结合RESP2(0.5)
水用OPC,用比较流行的SPC/E也凑合。
当前不牵扯生物分子,根本不涉及AMBER力场
弄清楚净电荷不为0是怎么回事。先看看每种分子的原子电荷加和是否为0(Multiwfn输出的RESP/RESP2电荷的小数位数已经很多了,不应当会偏离0明显)。如果每种分子都为0,体系总电荷只是零点零零零零几那种程度(加和时数值误差累积所致),直接当成0即可。

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-23 16:51 , Processed in 0.166451 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list