计算化学公社

 找回密码 Forget password
 注册 Register
Views: 131|回复 Reply: 6
打印 Print 上一主题 Last thread 下一主题 Next thread

[GROMACS] 求助各位大佬,我的分子动力学模拟修复和不修复蛋白都跑不通AMBER力场

[复制链接 Copy URL]

6

帖子

0

威望

17

eV
积分
23

Level 1 能力者

现在问题是不修复蛋白,蛋白的拓扑文件构建失败,但是如果修复蛋白的话在能量最小化的时候会失败
这是失败的指令:D:\autodock\benfagromacs>gmx mdrun -v -deffnm em
                      :-) GROMACS - gmx mdrun, 2025.2 (-:

Executable:   D:\fenzidonglixuemoniruanjian\GROMACS-2025.2_GPU_Windows-AMD64-AVX2_CUDA1263_msvc\bin\gmx.exe
Data prefix:  D:\fenzidonglixuemoniruanjian\GROMACS-2025.2_GPU_Windows-AMD64-AVX2_CUDA1263_msvc
Working dir:  D:\autodock\benfagromacs
Command line:
  gmx mdrun -v -deffnm em

Reading file em.tpr, VERSION 2025.2 (single precision)
1 GPU selected for this run.
Mapping of GPU IDs to the 1 GPU task in the 1 rank on this node:
  PP:0
PP tasks will do (non-perturbed) short-ranged interactions on the GPU
PP task will update and constrain coordinates on the CPU
Using 1 MPI thread
Using 16 OpenMP threads


Steepest Descents:
   Tolerance (Fmax)   =  1.00000e+03
   Number of steps    =        50000
Step=    0, Dmax= 1.0e-02 nm, Epot=  3.39863e+16 Fmax=         inf, atom= 182
Step=   14, Dmax= 1.2e-06 nm, Epot=  3.39863e+16 Fmax=         inf, atom= 182
Energy minimization has stopped because the force on at least one atom is not
finite. This usually means atoms are overlapping. Modify the input
coordinates to remove atom overlap or use soft-core potentials with the free
energy code to avoid infinite forces.
You could also be lucky that switching to double precision is sufficient to
obtain finite forces.

writing lowest energy coordinates.

Steepest Descents converged to machine precision in 15 steps,
but did not reach the requested Fmax < 1000.
Potential Energy  =  3.3986272e+16
Maximum force     =            inf on atom 182
Norm of force     =            inf

GROMACS reminds you: "A little less conversation, a little more action, please." (Elvis Presley)

202601051836418533..png (56.85 KB, 下载次数 Times of downloads: 0)

202601051836418533..png

202601051837251429..png (30.99 KB, 下载次数 Times of downloads: 1)

202601051837251429..png

202601051837438002..png (31.51 KB, 下载次数 Times of downloads: 0)

202601051837438002..png

6

帖子

0

威望

17

eV
积分
23

Level 1 能力者

2#
 楼主 Author| 发表于 Post on 2026-1-5 18:40:53 | 只看该作者 Only view this author
我使用的是spdbv修复的,求各位大佬给小弟看看怎么解决,被卡了一周了快烦死了

11

帖子

0

威望

81

eV
积分
92

Level 2 能力者

3#
发表于 Post on 2026-1-5 21:35:59 | 只看该作者 Only view this author
你的能量最小化的提醒信息有说,This usually means atoms are overlapping,看起来就是有两个或以上原子离得太近了,可以写脚本,或者把构型所有原子的xyz坐标导入到excel里面排序,找到离得太近的原子,然后再根据具体情况来判断怎么处理。

6万

帖子

99

威望

5万

eV
积分
124668

管理员

公社社长

4#
发表于 Post on 2026-1-6 01:21:59 | 只看该作者 Only view this author
肉眼检查结构,确保修复完之后的结构合理性
并且认真看产生拓扑文件整个过程的所有程序提示信息,所有warning都要留意
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6

帖子

0

威望

17

eV
积分
23

Level 1 能力者

5#
 楼主 Author| 发表于 Post on 2026-1-6 20:49:17 | 只看该作者 Only view this author
Myth 发表于 2026-1-5 21:35
你的能量最小化的提醒信息有说,This usually means atoms are overlapping,看起来就是有两个或以上原子离 ...

是的非常感谢大佬地回复,错误信息指出原子182(ASP 182的N原子)受力无限大。这通常意味着原子182与附近某个原子距离太近,导致范德华力无穷大。大佬我后来用了charmm-gui网站里的蛋白修复解决这个问题了,但是charmm-gui这个网站里要求自己写二硫键,大佬有没有专用来修复蛋白做amber的软件啊(自动修复的那种)我是一个计算机小白

6

帖子

0

威望

17

eV
积分
23

Level 1 能力者

6#
 楼主 Author| 发表于 Post on 2026-1-6 20:50:11 | 只看该作者 Only view this author
sobereva 发表于 2026-1-6 01:21
肉眼检查结构,确保修复完之后的结构合理性
并且认真看产生拓扑文件整个过程的所有程序提示信息,所有warn ...

大佬这个肉眼怎么观察啊,有没有推荐的自动化的软件啊,专门做amber力场的

6

帖子

0

威望

17

eV
积分
23

Level 1 能力者

7#
 楼主 Author| 发表于 Post on 2026-1-6 21:04:04 | 只看该作者 Only view this author
sobereva 发表于 2026-1-6 01:21
肉眼检查结构,确保修复完之后的结构合理性
并且认真看产生拓扑文件整个过程的所有程序提示信息,所有warn ...

大佬你看182号原子说是和别的原子离得太近,这种情况怎么辨别啊

本版积分规则 Credits rule

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

GMT+8, 2026-1-23 18:47 , Processed in 0.172325 second(s), 23 queries , Gzip On.

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