欢迎来到天天文库
浏览记录
ID:9213886
大小:229.72 KB
页数:10页
时间:2018-04-23
《gromacs中mdp文件注解小结》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、GROMACS中mdp文件注解小结——踏雪无痕;预处理title=OPLSLysozymeMD;标题,可任意定义(最长64个字,简单点好)cpp=/lib/cpp;预处理器,与C/C++的预处理器一样,默认为(/lib/cpp)include=;引用文件,即拓扑文件中引用其他文件的路径,引用方式与C/C++引用一样:格式:-I/home/../include-L/home/../libdefine=-DFLEXIBLE;预定义。默认无预定义。可以使用预定义方法控制模拟进程。;-DFLEXIBLE会告诉grompp,采用柔性的水代替刚性水到拓扑结构,这对于共轭
2、梯度或l-bfgs是必要的,而且可以使最陡下降法效果更好。;-DPOSRES会告诉grompp包括posre.itp到拓扑结构,用于位置限制动力学模拟。;模拟控制integrator=md;动力学模拟方法,即整合牛顿力学定理的方法,根据模拟目的不同选择不同。选项:;md:使用跳蛙算法(leap-frog)整合牛顿定律。;sd:另外一种跳蛙法统计整合(stochasticdynamics)。使用这个选项时,某个或者某些原子组(tc_grps)的温度设定为某特定温度(ref_t[K]),这些组运动反方向的摩擦常数可以设定为某一个值(tau_t[ps])。tcou
3、pl参数在这个选项中被忽略。这个参数的随机算子由ld_seed设定。(NOTE:这个方法中温度偏差的回原要比使用Berendesen热浴方法快一倍,即使使用相同的tau_t值。);bd:使用Euler整合方法处理Brownian或者坐标Langevin动力学模拟,模拟中的粒子的速度为所受力除以摩擦因子(bd_fric[amups-1),加上一个随机的热力学噪音(bd_temp[K])。当bd_fric=0时,模拟粒子的摩擦因子为其质量除以tau_t,这与sd方法一致。随机算子由ld_seed指定。;以下几种算法不是integrators,但是也在此处指定,主
4、要用于能量最优化模拟等。;steep:使用最速下降法进行能量优化,能量优化最大位置移动用emstep[nm]设定,能量最大容忍度由emtol[kJmol-1nm-1]决定。;cg:使用共轭梯度法进行能量优化,能量最大容忍度为emtol[kJmol-1nm-1]决定。在进行最速下降法能量优化之后再进行一次共轭梯度法能量最优化是十分有效的能量最优化综合方法,可以使用nstcgsteep设定。在要对能量优化进行常态分析时,最好使用双精度的GROMACS,以保证较高的精确度。;l-bfgs:根据low-memoryBroyden-Fletcher-Goldfarb-
5、Shannoapproach进行的准牛顿算法能量优化。实际中似乎比cg收敛更快,但由于必须进行纠正步骤,还没有实现并行。;nm:对tpr文件中的系统结构进行常态分析。GROMACS必须为双精度。;tpi:插入测试粒子(Testparticleinsertion)。拓扑结构的最后一个原子为测试粒子。;Starttimeandtimestepinpstinit=0;模拟开始时间。默认为0,单位为[ps]。(能量最小化用不用,只针对integrator为md、sd和bd的情况)dt=0.002;时间步长(2fs)。默认为0.001,单位为[ps]。(能量最小化用不
6、用,只针对integrator为md、sd和bd的情况)nsteps=50000;最大模拟步数。默认为0。;Forexactruncontinuationorredoingpartofaruninit_step=0;开始步骤,默认为0。第i步的时间由下式计算:t=tinit+dt*(init_step+i)。自由能lambda由此计算:lambda=init_lambda+delta_lambda*(init_step+i)。非平衡模拟参数依赖于模拟步数。所以重启停止的模拟或重做某部分模拟时,必须设定init_step为重启那一步。tpbconv自动完成此工
7、作。;modeforcenterofmassmotionremoval---质心移动comm_mode=Linear;对系统或者系统中各个组质心的操作,有三种选项:;Linear:移动质心的平动;;Angular:移动质心的转动和平动;;No:不对质心进行任何操作。;numberofstepsforcenterofmassmotionremovalnstcomm=10;对质心进行操作的频率,默认为10,单位为[steps]。;group(s)forcenterofmassmotionremovalcomm_grps=;对质心进行操作的组,可以是索引文件中的一
8、个,或者多个组。默认为整个系统。;Langevin动
此文档下载收益归作者所有