计算化学公社

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

[GROMACS] 求助-PLUMED+GROMACS跑 funnel metadynamics 报错

[复制链接 Copy URL]

29

帖子

0

威望

768

eV
积分
797

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
我打算用funnel metadynamics计算一个小分子配体和受体之间的结合能,然后照着PLUMED的官方教程去做,然后一开始模拟就报错,但是普通MD可以正常跑完,表明我的结构和拓扑应该没问题。
-------------------------------------------------------
Program:     gmx mdrun, version 2024.3-plumed_2.9.3

Unknown exception:
(exception type: N4PLMD6Plumed14ExceptionErrorE)

(tools/Grid.cpp:168) PLMD::GridBase::index_t PLMD::GridBase::getIndex(const
std::vector<unsigned int>&) const
Looking for a value outside the grid along the 0 dimension (arg name: fps.lp)

The above exception was the direct cause of the following exception:

(core/PlumedMain.cpp:972) void PLMD::PlumedMain::justCalculate()
An error happened while calculating funnel

For more information and tips for troubleshooting, please check the GROMACS
website at https://manual.gromacs.org/curre ... un-time-errors.html
-------------------------------------------------------
这是我利用FMAP输出的PLUMED输入文件上修改的:
MOLINFO STRUCTURE=ref2.pdb
WHOLEMOLECULES ENTITY0=1-4811
WHOLEMOLECULES ENTITY0=4812-4853
c1: CENTER ATOMS=4812-4853
c2: CENTER ATOMS={@mda:{name CA and resid 24:29:254:259:262}}
lig: COM ATOMS=4812,4813,4814,4815,4816,4817,4818,4819,4820,4821,4822,4823,4824,4825,4826,4827,4828,4829,4830,4831,4832,4833,4834,4835,4836,4837,4838,4839,4840,4841,4842,4843,4844,4845,4846,4847,4848,4849,4850,4851,4852,4853
d1: DISTANCE ATOMS=c1,c2

fps: FUNNEL_PS LIGAND=lig REFERENCE=ref2.pdb ANCHOR=4171 POINTS=3.048,1.843,3.653,3.126,1.469,4.468

FUNNEL ARG=fps.lp,fps.ld ZCC=1.5 ALPHA=0.6 RCYL=0.2 MINS=0.00 MAXS=4.00 KAPPA=35000 NBINS=500 NBINZ=500 FILE=BIAS LABEL=funnel

#METAD ARG=d1 SIGMA=0.015 HEIGHT=0.7 PACE=500 TEMP=298.15 BIASFACTOR=15 GRID_MIN=0.00 GRID_MAX=6.00 GRID_BIN=500 CALC_RCT RCT_USTRIDE=1 LABEL=metad

LOWER_WALLS ARG=fps.lp AT=0.40 KAPPA=35000 EXP=2 OFFSET=0 LABEL=lwall

UPPER_WALLS ARG=fps.lp AT=3.80 KAPPA=35000 EXP=2 OFFSET=0 LABEL=uwall

UPPER_WALLS ARG=d1 AT=4.10 KAPPA=35000 EXP=2 OFFSET=0 LABEL=distwall

UPPER_WALLS ARG=fps.ld AT=0.09 KAPPA=35000 EXP=2 LABEL=ld_wall

PRINT STRIDE=250 ARG=d1,fps.* FILE=COLVAR
我也把结构拓扑文件上传了,望大佬帮助。

gromacs.zip

1.33 MB, 下载次数 Times of downloads: 3

487

帖子

1

威望

1136

eV
积分
1643

Level 5 (御坂)

A Student

2#
发表于 Post on 2025-4-24 17:27:18 | 只看该作者 Only view this author
Funnel metadynamics的nature protocol有常见问题的解决方案的,参考我在这帖的连结和截图 http://bbs.keinsci.com/forum.php ... 0&fromuid=64740
敬仰一针见血的指责,厌倦别有用心的赞美。

29

帖子

0

威望

768

eV
积分
797

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2025-4-24 20:15:11 | 只看该作者 Only view this author
student0618 发表于 2025-4-24 17:27
Funnel metadynamics的nature protocol有常见问题的解决方案的,参考我在这帖的连结和截图 http://bbs.kein ...

谢谢

本版积分规则 Credits rule

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

GMT+8, 2025-8-14 08:53 , Processed in 0.166831 second(s), 23 queries , Gzip On.

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