求助：在能量最小化时提示错误信息One or more water molecules can not be settled

芮啦

各位老师，我定义了一个尺寸为4nm的立方体盒子，在盒子的中心放置适配体-铅离子，在CHARMM27力场下用Tip3p水解。进行到能量最小化这一步一直报错如下：
step 15: One or more water molecules can not be settled.
Check for bad contacts and/or reduce the timestep if appropriate.
Wrote pdb files with previous and current coordinates

step 21: One or more water molecules can not be settled.
Check for bad contacts and/or reduce the timestep if appropriate.
Wrote pdb files with previous and current coordinates

step 23: One or more water molecules can not be settled.
Check for bad contacts and/or reduce the timestep if appropriate.
Wrote pdb files with previous and current coordinates

......
step 492: One or more water molecules can not be settled.
Check for bad contacts and/or reduce the timestep if appropriate.
Wrote pdb files with previous and current coordinates

我只用金属离子进行分子动力学模拟，只用适配体进行分子动力学模拟，缩小盒子的尺寸等操作后依然出现同样的信息，提示One or more water molecules can not be settled. Check for bad contacts and/or reduce the timestep if appropriate. Wrote pdb files with previous and current coordinates。我检查了我的pdb文件和坐标文件没有问题，麻烦各位老师能够指明一下解决的方向。

libo371324

这个我以前遇到过，我将盒子的尺寸增加就解决啦

芮啦

libo371324 发表于 2021-4-6 17:14
这个我以前遇到过，我将盒子的尺寸增加就解决啦

谢谢您的回复，我把盒子的尺寸分别增加到6nm、7nm和8nm，但是还是出现了同样的错误提示。
我在mdp文件中增加了一项  define    =-DFLEXIBLE，但是又出现了新的错误提示：
Step 10, time 0.01 (ps)  LINCS WARNING
relative constraint deviation after LINCS:
rms 0.000003, max 0.000063 (between atoms 406 and 407)
bonds that rotated more than 30 degrees:
atom 1 atom 2  angle  previous, current, constraint length
411    412   30.6    0.1010   0.1010      0.1010
和错误提示：
Energy minimization has stopped, but the forces have not converged to the
requested precision Fmax < 1000 (which may not be possible for your system).
It stopped because the algorithm tried to make a new step whose size was too
small, or there was no change in the energy since last step. Either way, we
regard the minimization as converged to within the available machine
precision, given your starting configuration and EM parameters.
请问您有遇到过这种情况吗

sobereva

这和盒子大小没有任何直接关系。你的体系盒子根本没必要那么大，白浪费计算量。

mdp这么重要的文件都不提供

芮啦

sobereva 发表于 2021-4-7 00:50
这和盒子大小没有任何直接关系。你的体系盒子根本没必要那么大，白浪费计算量。

mdp这么重要的文件都不 ...

谢谢老师，我在上面上传了我的mdp文件
我在mdp文件中添加了  define    =-DFLEXIBLE之后没有再出现one or more water molecules can not be settled的报错，但是出现了新的错误提示：
Steepest Descents:
   Tolerance (Fmax)   =  1.00000e+03
   Number of steps    =         5000
Energy minimization has stopped, but the forces have not converged to the
requested precision Fmax < 1000 (which may not be possible for your system).
It stopped because the algorithm tried to make a new step whose size was too
small, or there was no change in the energy since last step. Either way, we
regard the minimization as converged to within the available machine
precision, given your starting configuration and EM parameters.

Double precision normally gives you higher accuracy, but this is often not
needed for preparing to run molecular dynamics.

writing lowest energy coordinates.

Back Off! I just backed up em.gro to ./#em.gro.2#

Steepest Descents converged to machine precision in 1722 steps,
but did not reach the requested Fmax < 1000.
Potential Energy  =  7.1431397e+10
Maximum force     =  4.5007404e+07 on atom 116075
Norm of force     =  2.0984681e+05

能不能麻烦老师再指点一下，谢谢您

牧生

芮啦 发表于 2021-4-7 09:00
谢谢老师，我在上面上传了我的mdp文件
我在mdp文件中添加了  define    =-DFLEXIBLE之后没有再出现one o ...

已经达到最大步数，但能量仍然没有达到设定最小值。

解决办法:反复进行几次能量最小化，或者改用双精度的gromacs进行能量最小化。

芮啦

牧生 发表于 2021-4-7 13:44
已经达到最大步数，但能量仍然没有达到设定最小值。

解决办法:反复进行几次能量最小化，或者改用双精 ...

谢谢您的回复，我反复进行了很多次能量最小化，依然提示相同的错误信息，我看到在能量最小化后生成的gro文件中，我的DNA链完全散开了，请问您能不能提示一下这可能是什么原因造成的呢？
打扰您啦，非常感谢

sobereva

芮啦 发表于 2021-4-7 14:20
谢谢您的回复，我反复进行了很多次能量最小化，依然提示相同的错误信息，我看到在能量最小化后生成的gro ...

原因和4L说的一样，存在内在不合理性

芮啦

sobereva 发表于 2021-4-8 01:13
原因和4L说的一样，存在内在不合理性

谢谢sob老师，我的初始结构和拓扑文件应该没有问题，在建模的过程中系统也没有警告和注释，这个问题困扰我一周了，一直没有解决，我想再请教一下老师，可能是力场参数有问题吗，哪些参数可能是我需要修改的地方
打扰老师了，谢谢您的回复

木林森

芮啦 发表于 2021-4-8 09:41
谢谢sob老师，我的初始结构和拓扑文件应该没有问题，在建模的过程中系统也没有警告和注释，这个问题困扰 ...

请问楼主你这个问题解决了吗？我也遇到了相同的问题，可以分享一下吗？

lonemen

牧生 发表于 2021-4-7 13:44
已经达到最大步数，但能量仍然没有达到设定最小值。

解决办法:反复进行几次能量最小化，或者改用双精 ...

请问您说的反复进行能量最小化，是如何操作的呢？用上次最小化的结构，继续跑下一次最小化吗？如此反复进行？

牧生

lonemen 发表于 2021-6-19 14:16
请问您说的反复进行能量最小化，是如何操作的呢？用上次最小化的结构，继续跑下一次最小化吗？如此反复进 ...

是的。如果反复多次还不能最小化，就考虑用双精度，
如果还不行，就考虑是自己做的模型是否合理。

lonemen

牧生 发表于 2021-6-19 14:45
是的。如果反复多次还不能最小化，就考虑用双精度，
如果还不行，就考虑是自己做的模型是否合理。

好的，多谢指点和建议

芮啦

木林森 发表于 2021-6-19 11:05
请问楼主你这个问题解决了吗？我也遇到了相同的问题，可以分享一下吗？

我在windows下用gromacs计算出现了上面的提到的问题。但在linux下，在所有条件都相同的情况下试着计算没有出现问题了

木林森

芮啦 发表于 2021-6-21 18:32
我在windows下用gromacs计算出现了上面的提到的问题。但在linux下，在所有条件都相同的情况下试着计算没 ...

我现在就是在linux下，同样的问题一直报错