|
|
如题,我的自由能计算参数设置如下:
free-energy = yes
couple-moltype = BIP
couple-lambda0 = none
couple-lambda1 = vdw-q
couple-intramol = no
separate-dhdl-file = yes
sc-alpha = 0.5
sc-power = 1
sc-sigma = 0.3
init-lambda-state = $LAMBDA$
coul-lambdas = 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 0.65 0.70 0.75 0.80 0.85 0.90 0.95 1.00
vdw-lambdas = 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 0.65 0.70 0.75 0.80 0.85 0.90 0.95 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00
nstdhdl = 100
calc-lambda-neighbors = -1
用alchemical_analysis分析得到采样窗口的重叠分布图和相应的重叠矩阵如下,这样的采样窗口该如何改进
|
-
cfm.png
(249.52 KB, 下载次数 Times of downloads: 46)
重叠分布图
-
O_MBAR.png
(253.8 KB, 下载次数 Times of downloads: 29)
重叠矩阵
|