计算化学公社

标题: CHARMM-GUI生成的文件夹之后怎么用gromacs模拟 [打印本页]

作者
Author: shanwo 时间: 2024-7-22 12:48
标题: CHARMM-GUI生成的文件夹之后怎么用gromacs模拟
各位大佬，我在做一个膜蛋白体系的模拟，用CHARMM-GUI生成了文件夹，里面的gromacs文件夹有一个readme文件，里面的教程命令是这样的：
#!/bin/bash

# Equilibration
init="step5_input"
rest_prefix="step5_input"
mini_prefix="step6.0_minimization"
equi_prefix="step6.%d_equilibration"
prod_prefix="step7_production"
prod_step="step7"

cnt=1
cntmax=6

while [ $cnt -le $cntmax ]
do
pcnt=$(($cnt - 1))
istep=$(printf "$equi_prefix" $cnt)
if [ $cnt -eq 1 ]; then
      pstep=$mini_prefix
else
      pstep=$(printf "$equi_prefix" $pcnt)
fi

gmx grompp -f ${istep}.mdp -o ${istep}.tpr -c ${pstep}.gro -r ${rest_prefix}.gro -p topol.top -n index.ndx
gmx mdrun -v -deffnm ${istep}

cnt=$(($cnt + 1))
done

我打包好的sh文件运行总是报错：

Error in user input:Invalid input values  In option s Required option was not provided, and the default file 'topol' does not exist or is not accessible. The following extensions were tried to complete the file name:    .tpr

和

Error in user input:
Invalid command-line options
  In command-line option -c
File 'step6.5_equilibration.gro' does not exist or is not accessible.
The file could not be opened.
   Reason: No such file or directory
   (call to fopen() returned error code 2)

作者
Author: Biny 时间: 2024-8-30 19:42
关于第一个问题，我也遇到了，参考GROMACS Tutorials的命令，我发现平衡和模拟过程的命名是不含有-v的，所以我就删掉了，删掉之后可以正常允运行，没有报错

作者
Author: student0618 时间: 2024-8-31 02:16
看上去minimization時出问题了，导致后来跑不了。
检查一下最后跑到的Log (step6.0_minimization.log)，大概会见到很多warnings。
也用例如VMD检查stepxxx.pdb，那些是出错前的座标。

作者
Author: HNUST 时间: 2024-12-14 18:41
readme是一个脚本

作者
Author: CrysW555 时间: 2024-12-24 14:21
本帖最后由 CrysW555 于 2024-12-24 14:36 编辑

可以把能量最小化的log发出来看一看

作者
Author: 许正东 时间: 3 day ago

HNUST 发表于 2024-12-14 18:41
readme是一个脚本

请问readme脚本怎么运行呢？

作者
Author: HNUST 时间: 3 day ago

许正东发表于 2025-8-19 15:56
请问readme脚本怎么运行呢？

csh README
要修改一点点代码

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