求助：GROMACS用TI算吸附自由能时lambda设置问题

wuzhiyi

用
vdw-lambdas  = 0.00 0.05 0.075 0.10 0.0.125 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50
coul-lambdas  = 0.00 0.00 0.0 0.0 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
来跑rerun

活性炭

wuzhiyi 发表于 2021-7-5 03:21
用
vdw-lambdas  = 0.00 0.05 0.075 0.10 0.0.125 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50
coul-lambda ...

谢谢您的回复！我试下

王家小姐

您好，请问您的rerun跑成功了吗?我用rerun重新得到的dhld.xvg后，用alchemical_analysis.py计算结合能，老是提示我“File "/home/test2/K80-test/AFEG/Bina/solvent/dHdl_files/alchemical-analysis-master/alchemical_analysis/alchemical_analysis.py", line 1296, in <module>    main()
  File "/home/test2/K80-test/AFEG/Bina/solvent/dHdl_files/alchemical-analysis-master/alchemical_analysis/alchemical_analysis.py", line 1225, in main
    nsnapshots, lv, dhdlt, u_klt = parser_gromacs.readDataGromacs(P)
  File "/home/test2/K80-test/AFEG/Bina/solvent/dHdl_files/alchemical-analysis-master/alchemical_analysis/parser_gromacs.py", line 324, in readDataGromacs
    nsnapshots[nf,nf] += unixlike.wcPy(f.filename) - f.skip_lines - 1*bLenConsistency
IndexError: index 17 is out of bounds for axis 1 with size 17”

活性炭

王家小姐 发表于 2021-8-1 10:16
您好，请问您的rerun跑成功了吗?我用rerun重新得到的dhld.xvg后，用alchemical_analysis.py计算结合能，老 ...

你这个报错是将两组不同的lambda值一起分析引起的，你用同一组lambda值重新rerun，再分析

王家小姐

活性炭 发表于 2021-8-30 09:44
你这个报错是将两组不同的lambda值一起分析引起的，你用同一组lambda值重新rerun，再分析

请问具体要怎么做呀？我第一次接触不是很懂，您的命令和准备文件方便不方便给我参考一下呀？

活性炭

王家小姐 发表于 2022-2-19 11:23
请问具体要怎么做呀？我第一次接触不是很懂，您的命令和准备文件方便不方便给我参考一下呀？

单独跑完新增的采样窗口，不用输出dhdl文件，只需要得到轨迹文件即可，   然后根据轨迹用rerun把所有采样窗口用新的mdp重新跑一遍得到dhdl文件
最后，假设你增加采样窗口之后，现在一共有4个窗口 lambda.1 ~ lambda.4
gmx mdrun -deffnm prod -multidr lambda.1 lambda.2 lambda.3 lambda.4 -replex 1000 -nex 1000000

王家小姐

活性炭 发表于 2022-3-14 15:48
单独跑完新增的采样窗口，不用输出dhdl文件，只需要得到轨迹文件即可，   然后根据轨迹用rerun把所有采样 ...

请问“gmx mdrun -deffnm prod -multidr lambda.1 lambda.2 lambda.3 lambda.4 -replex 1000 -nex 1000000”是不是就是跑的rerun，还是需要用到rerun命令呢？如果用rerun命令的话，您能给个命令用法吗？我看教程尝试的命令都不对，多谢您

活性炭

王家小姐 发表于 2022-5-16 17:09
请问“gmx mdrun -deffnm prod -multidr lambda.1 lambda.2 lambda.3 lambda.4 -replex 1000 -nex 100000 ...

需要把所有采样窗口一起rerun一下
比如对第一个窗口，rerun的命令可以参考：
gmx_mpi grompp -f prod.1.mdp -c prod.gro -p topol.top -o prod.tpr
gmx_mpi mdrun -deffnm prod -rerun prod.xtc -dhdl dhdl
第二个窗口用prod.2.mdp，其它参数不变
其余窗口以此类推，用对应的mdp文件即可

ZINCxxy

你好，我想请教一下，我在使用TI算吸附自由能的时候，算完之后用生成的xvg和edr分析后提示我没有delta H这一项，这是什么原因所致呢，是因为mdp中设置的一些错误吗？此问题困扰颇久，非常感谢解答