计算化学公社

标题: 如何在amber模擬時對特定的兩個原子做restraint?? [打印本页]

作者
Author: xroxzero 时间: 2017-9-28 16:24
标题: 如何在amber模擬時對特定的兩個原子做restraint??
各位老師好

我想在進行模擬時進行兩個原子的restraint，我在input檔裡加了 restraint_wt=10 restraintmask="@12,17"

我還想在這兩個指定原子加上距離上的限制例如:將這兩個原子restraint在3A

請問我該怎麼設定呢??

謝謝各位老師

作者
Author: agent99 时间: 2017-9-28 23:19
你用的是cartesian restraint，那么所有原子的平衡位置都是你的参考坐标（运行命令里 -ref 后面的文件）。对于你的情况，显然是个internal coordinate的restraint，推荐用nmropt，然后只需写一个两行的小文件定义一个harmonic restraint就行了，具体请看amber的tutorial

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