read
  mol pqr _pid~water_com.pqr
  mol pqr _pid~water_pro.pqr
  mol pqr _pid~water_lig.pqr
end

ELEC name _pid~water_com
  mg-auto
  mol 1
  dime   161  129  161  # 格点数目
  cglen  78.72  68.1495  80.748  # 粗略格点长度
  fglen  57.48  50.433  58.832  # 细密格点长度
  fgcent 62.819  62.8575  29.631  # 细密格点中心
  cgcent 62.819  62.8575  29.631  # 粗略格点中心

  lpbe              # PB方程求解方法, lpbe(线性), npbe(非线性), smbpe(大小修正)

  temp  300         # 温度

  pdie  2           # 溶质介电常数
  sdie  80          # 溶剂介电常数, 真空1, 水78.54
  bcfl  mdh         # 粗略格点PB方程的边界条件, zero, sdh/mdh(single/multiple Debye-Huckel), focus, map
  srfm  smol        # 构建介质和离子边界的模型, mol(分子表面), smol(平滑分子表面), spl2/4(三次样条/7阶多项式)
  chgm  spl4        # 电荷映射到格点的方法, spl0/2/4, 三线性插值, 立方/四次B样条离散
  swin  0.3         # 立方样条的窗口值, 仅用于 srfm=spl2/4

  srad  1.4         # 溶剂探测半径
  sdens 10          # 表面密度, 每A^2的格点数, (srad=0)或(srfm=spl2/4)时不使用

  ion  1 0.15 0.95  # 阳离子的电荷, 浓度, 半径
  ion -1 0.15 1.81  # 阴离子

  calcforce  no
  calcenergy comps
end

ELEC name _pid~water_com~VAC
  mg-auto  mol 1
  dime   161  129  161  # 格点数目
  cglen  78.72  68.1495  80.748  # 粗略格点长度
  fglen  57.48  50.433  58.832  # 细密格点长度
  fgcent 62.819  62.8575  29.631  # 细密格点中心
  cgcent 62.819  62.8575  29.631  # 粗略格点中心

  lpbe              # PB方程求解方法, lpbe(线性), npbe(非线性), smbpe(大小修正)

  temp  300         # 温度

  pdie  2           # 溶质介电常数
  sdie  1          # 溶剂介电常数, 真空1, 水78.54
  bcfl  mdh         # 粗略格点PB方程的边界条件, zero, sdh/mdh(single/multiple Debye-Huckel), focus, map
  srfm  smol        # 构建介质和离子边界的模型, mol(分子表面), smol(平滑分子表面), spl2/4(三次样条/7阶多项式)
  chgm  spl4        # 电荷映射到格点的方法, spl0/2/4, 三线性插值, 立方/四次B样条离散
  swin  0.3         # 立方样条的窗口值, 仅用于 srfm=spl2/4

  srad  1.4         # 溶剂探测半径
  sdens 10          # 表面密度, 每A^2的格点数, (srad=0)或(srfm=spl2/4)时不使用

  ion  1 0.15 0.95  # 阳离子的电荷, 浓度, 半径
  ion -1 0.15 1.81  # 阴离子

  calcforce  no
  calcenergy comps
end

APOLAR name _pid~water_com~SAS
  mol 1

  temp  300  # 温度
  srfm  sacc # 构建溶剂相关表面或体积的模型
  swin  0.3  # 立方样条窗口(A), 用于定义样条表面

  # SASA模型
  srad  1.4    # 探测半径(A)
  gamma 1      # 表面张力(kJ/mol-A^2)

  #gamma const 0.027     0        # 表面张力, 常数
  #gamma const 0.0226778 3.84928  # 表面张力, 常数

  press  0     # 压力(kJ/mol-A^3)
  bconc  0     # 溶剂本体密度(A^3)
  sdens 10
  dpos  0.2
  grid  0.1 0.1 0.1

  # SAV模型 截距0
  #srad  1.29      # SAV探测半径(A)
  #press 0.234304  # 压力(kJ/mol-A^3)

  # WCA模型
  #srad   1.25           # 探测半径(A)
  #sdens  200            # 表面的格点密度(1/A)
  #dpos   0.05           # 表面积导数的计算步长
  #bconc  0.033428       # 溶剂本体密度(A^3)
  #grid   0.45 0.45 0.45 # 算体积分时的格点间距(A)

  calcforce no
  calcenergy total
end

#print elecEnergy _pid~water_com - _pid~water_com~VAC end
#print apolEnergy _pid~water_com~SAS end
#quit
 ELEC name _pid~water_pro
  mg-auto
  mol 2
  dime   161  129  161  # 格点数目
  cglen  78.72  68.1495  80.748  # 粗略格点长度
  fglen  57.48  50.433  58.832  # 细密格点长度
  fgcent 62.819  62.8575  29.631  # 细密格点中心
  cgcent 62.819  62.8575  29.631  # 粗略格点中心

  lpbe              # PB方程求解方法, lpbe(线性), npbe(非线性), smbpe(大小修正)

  temp  300         # 温度

  pdie  2           # 溶质介电常数
  sdie  80          # 溶剂介电常数, 真空1, 水78.54
  bcfl  mdh         # 粗略格点PB方程的边界条件, zero, sdh/mdh(single/multiple Debye-Huckel), focus, map
  srfm  smol        # 构建介质和离子边界的模型, mol(分子表面), smol(平滑分子表面), spl2/4(三次样条/7阶多项式)
  chgm  spl4        # 电荷映射到格点的方法, spl0/2/4, 三线性插值, 立方/四次B样条离散
  swin  0.3         # 立方样条的窗口值, 仅用于 srfm=spl2/4

  srad  1.4         # 溶剂探测半径
  sdens 10          # 表面密度, 每A^2的格点数, (srad=0)或(srfm=spl2/4)时不使用

  ion  1 0.15 0.95  # 阳离子的电荷, 浓度, 半径
  ion -1 0.15 1.81  # 阴离子

  calcforce  no
  calcenergy comps
end

ELEC name _pid~water_pro~VAC
  mg-auto  mol 2
  dime   161  129  161  # 格点数目
  cglen  78.72  68.1495  80.748  # 粗略格点长度
  fglen  57.48  50.433  58.832  # 细密格点长度
  fgcent 62.819  62.8575  29.631  # 细密格点中心
  cgcent 62.819  62.8575  29.631  # 粗略格点中心

  lpbe              # PB方程求解方法, lpbe(线性), npbe(非线性), smbpe(大小修正)

  temp  300         # 温度

  pdie  2           # 溶质介电常数
  sdie  1          # 溶剂介电常数, 真空1, 水78.54
  bcfl  mdh         # 粗略格点PB方程的边界条件, zero, sdh/mdh(single/multiple Debye-Huckel), focus, map
  srfm  smol        # 构建介质和离子边界的模型, mol(分子表面), smol(平滑分子表面), spl2/4(三次样条/7阶多项式)
  chgm  spl4        # 电荷映射到格点的方法, spl0/2/4, 三线性插值, 立方/四次B样条离散
  swin  0.3         # 立方样条的窗口值, 仅用于 srfm=spl2/4

  srad  1.4         # 溶剂探测半径
  sdens 10          # 表面密度, 每A^2的格点数, (srad=0)或(srfm=spl2/4)时不使用

  ion  1 0.15 0.95  # 阳离子的电荷, 浓度, 半径
  ion -1 0.15 1.81  # 阴离子

  calcforce  no
  calcenergy comps
end

APOLAR name _pid~water_pro~SAS
  mol 2

  temp  300  # 温度
  srfm  sacc # 构建溶剂相关表面或体积的模型
  swin  0.3  # 立方样条窗口(A), 用于定义样条表面

  # SASA模型
  srad  1.4    # 探测半径(A)
  gamma 1      # 表面张力(kJ/mol-A^2)

  #gamma const 0.027     0        # 表面张力, 常数
  #gamma const 0.0226778 3.84928  # 表面张力, 常数

  press  0     # 压力(kJ/mol-A^3)
  bconc  0     # 溶剂本体密度(A^3)
  sdens 10
  dpos  0.2
  grid  0.1 0.1 0.1

  # SAV模型 截距0
  #srad  1.29      # SAV探测半径(A)
  #press 0.234304  # 压力(kJ/mol-A^3)

  # WCA模型
  #srad   1.25           # 探测半径(A)
  #sdens  200            # 表面的格点密度(1/A)
  #dpos   0.05           # 表面积导数的计算步长
  #bconc  0.033428       # 溶剂本体密度(A^3)
  #grid   0.45 0.45 0.45 # 算体积分时的格点间距(A)

  calcforce no
  calcenergy total
end

#print elecEnergy _pid~water_pro - _pid~water_pro~VAC end
#print apolEnergy _pid~water_pro~SAS end
#quit
 ELEC name _pid~water_lig
  mg-auto
  mol 3
  dime   65  33  65  # 格点数目
  cglen  20.3025  15.651  29.172  # 粗略格点长度
  fglen  18.535  15.434  24.448  # 细密格点长度
  fgcent 64.2885  60.25  28.627  # 细密格点中心
  cgcent 64.2885  60.25  28.627  # 粗略格点中心

  lpbe              # PB方程求解方法, lpbe(线性), npbe(非线性), smbpe(大小修正)

  temp  300         # 温度

  pdie  2           # 溶质介电常数
  sdie  80          # 溶剂介电常数, 真空1, 水78.54
  bcfl  mdh         # 粗略格点PB方程的边界条件, zero, sdh/mdh(single/multiple Debye-Huckel), focus, map
  srfm  smol        # 构建介质和离子边界的模型, mol(分子表面), smol(平滑分子表面), spl2/4(三次样条/7阶多项式)
  chgm  spl4        # 电荷映射到格点的方法, spl0/2/4, 三线性插值, 立方/四次B样条离散
  swin  0.3         # 立方样条的窗口值, 仅用于 srfm=spl2/4

  srad  1.4         # 溶剂探测半径
  sdens 10          # 表面密度, 每A^2的格点数, (srad=0)或(srfm=spl2/4)时不使用

  ion  1 0.15 0.95  # 阳离子的电荷, 浓度, 半径
  ion -1 0.15 1.81  # 阴离子

  calcforce  no
  calcenergy comps
end

ELEC name _pid~water_lig~VAC
  mg-auto  mol 3
  dime   65  33  65  # 格点数目
  cglen  20.3025  15.651  29.172  # 粗略格点长度
  fglen  18.535  15.434  24.448  # 细密格点长度
  fgcent 64.2885  60.25  28.627  # 细密格点中心
  cgcent 64.2885  60.25  28.627  # 粗略格点中心

  lpbe              # PB方程求解方法, lpbe(线性), npbe(非线性), smbpe(大小修正)

  temp  300         # 温度

  pdie  2           # 溶质介电常数
  sdie  1          # 溶剂介电常数, 真空1, 水78.54
  bcfl  mdh         # 粗略格点PB方程的边界条件, zero, sdh/mdh(single/multiple Debye-Huckel), focus, map
  srfm  smol        # 构建介质和离子边界的模型, mol(分子表面), smol(平滑分子表面), spl2/4(三次样条/7阶多项式)
  chgm  spl4        # 电荷映射到格点的方法, spl0/2/4, 三线性插值, 立方/四次B样条离散
  swin  0.3         # 立方样条的窗口值, 仅用于 srfm=spl2/4

  srad  1.4         # 溶剂探测半径
  sdens 10          # 表面密度, 每A^2的格点数, (srad=0)或(srfm=spl2/4)时不使用

  ion  1 0.15 0.95  # 阳离子的电荷, 浓度, 半径
  ion -1 0.15 1.81  # 阴离子

  calcforce  no
  calcenergy comps
end

APOLAR name _pid~water_lig~SAS
  mol 3

  temp  300  # 温度
  srfm  sacc # 构建溶剂相关表面或体积的模型
  swin  0.3  # 立方样条窗口(A), 用于定义样条表面

  # SASA模型
  srad  1.4    # 探测半径(A)
  gamma 1      # 表面张力(kJ/mol-A^2)

  #gamma const 0.027     0        # 表面张力, 常数
  #gamma const 0.0226778 3.84928  # 表面张力, 常数

  press  0     # 压力(kJ/mol-A^3)
  bconc  0     # 溶剂本体密度(A^3)
  sdens 10
  dpos  0.2
  grid  0.1 0.1 0.1

  # SAV模型 截距0
  #srad  1.29      # SAV探测半径(A)
  #press 0.234304  # 压力(kJ/mol-A^3)

  # WCA模型
  #srad   1.25           # 探测半径(A)
  #sdens  200            # 表面的格点密度(1/A)
  #dpos   0.05           # 表面积导数的计算步长
  #bconc  0.033428       # 溶剂本体密度(A^3)
  #grid   0.45 0.45 0.45 # 算体积分时的格点间距(A)

  calcforce no
  calcenergy total
end

#print elecEnergy _pid~water_lig - _pid~water_lig~VAC end
#print apolEnergy _pid~water_lig~SAS end
#quit

