amber模擬kcl水溶液
最近剛開始學習amber軟體，看網上的教程勉強知道怎麼操作這個amber了。就暫時跑了個分子動力學，其他的啥也沒處理。先把我的操作過程記錄下來吧，免得日後忘記。
一、構建kcl.pdb結構
利用GaussView構建結構，其實很簡單，就分別選擇K和Cl，然後在窗口點兩下，保存成kcl.pdb。
二、利用xleap獲得sander的輸入.top拓撲文件和.crd坐標文件
①打開xleap，在命令行輸入xleap
②導入spce水的leaprc文件，這個文件包含一些基本的離子的參數，我這裡使用的是spce水。amber還自帶有其他的水模型文件，在$AMBERHOME/dat/leap/cmd下面。
source leaprc.water.spce
kcl = loadpdb kcl.pdb
edit kcl
③修改kcl的原子類型和電荷
滑鼠選擇兩個原子，然後Edit-Edit selected atoms，修改原子類型為K+、Cl-，這是因為第②導入的力場裡面K、Cl的原子類型就是K+和Cl-。修改完之後保存
④設置盒子以及溶劑化
set kcl box { 30 30 30 }
solvatebox kcl TIP3PBOX 0
⑤保存amber拓撲和坐標文件
saveamberparm kcl kcl-water-box.top kcl-water-box.crd
savepdb kcl kcl-water-box.pdb
三、跑能量最小化
最小化文件enmin.in如下，參數控制在&cntrl … /裡面，注釋符號是啥我暫時還沒搞清楚
Minimize
&cntrl
imin=1,
ntx=1,
irest=0,
maxcyc=2000,
ncyc=1000,
ntpr=100,
ntwx=0,
cut=8.0,
/
運行的指令，我寫了個腳本，提交到課題組的伺服器集群上面計並行算。MPIRUN和EXEC是相應的目錄
MPIRUN=/opt/intel/impi/5.0.2.044/intel64/bin/mpirun
EXEC=/home/lkq/mysoft/amber18/bin/sander.MPI
TOP=kcl-water-box.top
${MPIRUN} -np ${nprocs} ${EXEC} -O -i enmin.in -p ${TOP} -c kcl-water-box.crd -o kcl-min.out -r kcl-min.rst
四、跑npt系綜
npt系綜的文件npt.in如下
kcl: 100ps npt run
&cntrl
imin = 0,
irest = 0,
ntx = 1,
ntb = 2,
ntp = 1,
cut = 10,
ntc = 2,
ntf = 2,
tempi = 300.0,
temp0 = 300.0,
ntt = 3,
gamma_ln = 1.0,
nstlim = 50000, dt = 0.002
ntpr = 100, ntwx = 1000, ntwr = 1000
/
END
END
運行命令
${MPIRUN} -np ${nprocs} ${EXEC} -O -i npt.in -p ${TOP} -c kcl-min.rst -o kcl-npt.out -r kcl-npt.rst -x kcl-npt.mdcrd
