计算化学公社

标题: top文件不引用posre.itp能量爆炸问题 [打印本页]

作者
Author: 仙人球没有球 时间: 2024-4-23 14:19
标题: top文件不引用posre.itp能量爆炸问题
我现在在做一个蛋白模拟的体系，它是由一个蛋白和若干水分子还有平衡电荷的NA和CL离子组成。蛋白的电荷是分开写的，主链部分是根据残基库里的电荷添加的，侧链部分是利用gaussian生成的，蛋白的拓扑是我在gromos力场中添加了rtp然后用gmx pdb2gmx生成的，缺失的键长键角二面角参数是我用antechamber一系列操作生成prmtop手动添加的。目前遇到的问题是能量最小化过程能进行，然后固定主链原子跑NVT也能正常跑，但是和posre里的力常数有关，如果力常数变小，他就会报错，显示段错误（核心已转储）。由于我需要对蛋白进行一个拉伸，后面放开它的约束，在NPT系综下模拟，所以后续的posre力常数需要为0。但是目前无法解决NVT阶段不用力常数约束就不能模拟的问题，还请各位老师伸出援助之手！附件是我的相关模拟文件~

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