|
|
请问大家,谁做过umbrella或者constant-force
我对参考点(参考组)和弹簧点有点不理解,虽然我已经将仔细阅读过李老师的博文和sob老师的博文了,也pull成功了,但是我还是处于一知半解的状态。
; Pull code
pull = yes ; Pull start
pull_ngroups = 2 ; The number of pull groups
pull_ncoords = 1 ; The number of pull coordinates
pull_group1_name = A ; The name of the pull group
pull_group2_name = B ; The name of the pull group (It is frozen here)
pull_coord1_type = constant-force ; Center of mass pulling using a linear potential and therefore a constant force
pull-coord1_geometry = direction-periodic ; As direction, but allows the distance to be larger than half the box size
pull_coord1_dim = Y N N ; the distance components to be used with geometry distance
pull_coord1_vec = 1.0 0.0 0.0 ; Pull direction
pull_coord1_k = 1000 ; pulling force kJ mol^-1 nm^-1
pull-coord1-groups = 1 2 ; The two groups indices should be given on which this pull coordinate will operate
pull-nstxout = 10 ; frequency for writing out the COMs of all the pull group
pull-nstfout = 1 ; frequency for writing out the force of all the pulled group
这是我模仿教程的pull code写的。使用的拉方式为constant-force所以没有弹簧点,有A、B两个拉伸组,但是我限制了B的运动,所以只能拉伸A组。
请问以上代码没有定义参考组吧?
其实就像想问,拉伸方式是umbrella或constant-force怎么设置参考组和弹簧点。
|
|