计算化学公社

标题: 求助卡那霉素+水在NVT平衡时,分子结构出现问题,多次尝试仍然报错 [打印本页]

作者
Author:
HYUAN    时间: 2023-7-5 15:14
标题: 求助卡那霉素+水在NVT平衡时,分子结构出现问题,多次尝试仍然报错
老师好,我用高斯09的opt freq rb3lyp/6-311+g(d,p),优化卡那霉素的结构,优化成功后通过multiwfn计算RESP电荷(计算了真空和水溶液下两个,取其平均值),然后通过sobtop生成其拓扑结构,用于GROMACS中进行MD模拟,一开始单独将卡那霉素在真空中进行em,nvt,npt,通过逐步升温的,以及减小时间步长等方法,卡那霉素在真空中可以成功平衡,之后将平衡后的结构提取出来,重新建立盒子,加入水分子,再进行能量最小化是可以运行成功,但后续的NVT平衡一直报错,从报错文件来看是卡那霉素的氨基与其两侧的羟基有交联成键,尝试了多次能量最小化、减小时间步长、逐步升温等方法,都不能使其平衡,一直出现类似这样的报错:
Fatal error:
1 particles communicated to PME rank 1 are more than 2/3 times the cut-off out
of the domain decomposition cell of their charge group in dimension x.
This usually means that your system is not well equilibrated.

附件给出了卡那霉素的拓扑文件,以及通过高斯优化出的结构,还有进行Nvt平衡时的报错文件,已经被困扰了好多天了,一直平衡不成功,直接跑MD也跑不了,希望老师可以帮忙看一下,这个是由什么问题造成的,非常感谢!
(, 下载次数 Times of downloads: 2)
(, 下载次数 Times of downloads: 2)
(, 下载次数 Times of downloads: 1) (, 下载次数 Times of downloads: 3) (, 下载次数 Times of downloads: 2) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0)

作者
Author:
sobereva    时间: 2023-7-6 10:43
确保给sobtop用的mol2文件里的成键关系是正确的,不正确就在GaussView里调好了再保存mol2
作者
Author:
Arkin    时间: 2023-7-6 10:54
不一定是 mol2 里成键的问题。把完整的 itp 文件或 top 文件发出来,用 `gmx grompp -pp ...`,我大概知道是怎么回事,但是需要足够的数据支持。
作者
Author:
HYUAN    时间: 2023-7-6 14:52
老师好,附件里给出了原始从sobtop中得到的拓扑文件,麻烦老师可以再帮忙看一下,谢谢您。
(, 下载次数 Times of downloads: 1) (, 下载次数 Times of downloads: 1) (, 下载次数 Times of downloads: 1)

作者
Author:
HYUAN    时间: 2023-7-6 16:39
补充一下.mol2文件,我自己看的话成键结构好像显示都很正常,麻烦老师了 (, 下载次数 Times of downloads: 1)
作者
Author:
Arkin    时间: 2023-7-6 18:29
虽然看不出来你的这个分子是怎么摆的,但是能确定 MD 过程中,重合的是 1-5 位的氮和氢原子。



看以下这篇文献:https://doi.org/10.1002/jcc.20090:A biomolecular force field based on the free enthalpy of hydration and solvation: The GROMOS force-field parameter sets 53A5 and 53A6

具体的解释牵涉到力场的开发和非键能的计算。

其中的引用是:

In principle, all atom pairs should be included in this sum, but generally, the sum is restricted to a subset of atom pairs. First, covalently bound neighboring atoms (first neighbors) and second neighbors are standardly excluded from this sum. They are already directly interacting through the bonded interactions, and are thus considered to be excluded atom pairs for the nonbonded interactions. In addition, the third or 1– 4 covalently bound neighbor atoms that are part of or bound to aromatic rings are also excluded from the nonbonded interactions. This makes it easier for the improper dihedral angle interactions to keep the atoms of and bound to a planar aromatic ring in one plane. Furthermore, there are a few exceptions to this general rule. These involve hydrogen atoms that have a repulsive van der Waals interaction equal to zero. In cases where the hydrogen is in a 1– 4 or 1–5 covalently bound neighbor position to an oppositely charged atom, X45, and the heavy atom, X1, to which the hydrogen is bonded is an excluded atom with respect to atom X45 an additional exclusion of the hydrogen with atom X45 is sometimes needed. This is because the bonded interactions may fail to prevent the hydrogen atom from collapsing onto atom X45. For example, in the nucleotide adenosine atoms H61 and H62 are excluded from atom N7 for this reason.

所以,很好的一种解决方式是,你在 top 文件里额外添加 `[exclusions] ` (用对应的原子序号),把这个 LJ 为零的氢和和氮之间的非键相互左右给排除掉就可以了。

作者
Author:
HYUAN    时间: 2023-7-6 19:28
Arkin 发表于 2023-7-6 18:29
虽然看不出来你的这个分子是怎么摆的,但是能确定 MD 过程中,重合的是 1-5 位的氮和氢原子。

好的,非常感谢您的认真解疑,我去尝试一下这个办法!




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3