水モデル

計算化学において、古典的水モデル（みずモデル）は、水クラスターや液体の水、露な溶媒を用いた水溶液のシミュレーションのために用いられる。水モデルは量子力学や分子力学、実験結果、これらの組み合わせによって決定される。分子の特定の性質を模倣するため、多くの種類のモデルが開発されている。一般的に、これらは (i) 「サイト」と呼ばれる相互作用点の数、(ii) モデルが剛直なのか柔軟なのか、(iii) モデルが分極効果を含んでいるか、という3つの点によって分類することができる。

露な（明示的な）水モデルの代わりとしては非明示的溶媒和モデル（連続体モデルとも呼ばれる）がある。非明示的溶媒和モデルの例としては、COSMO溶媒和モデルや分極連続体モデル（PCM）、ハイブリッド溶媒和モデルがある^[1]。

単純な水モデル

剛直なモデルは非結合相互作用に依存する最も単純な水モデルとして知られる。これらのモデルにおいて、結合性相互作用はホロノミック拘束によって陰に扱われる。静電相互作用はクーロンの法則を用いてモデル化され、分散力および斥力はレナード-ジョーンズ・ポテンシャルを用いてモデル化される^[2]^[3]。TIP3PおよびTIP4Pといったモデルに対するポテンシャルは

E_{ab}=\sum _{i}^{{\text{on }}a}\sum _{j}^{{\text{on }}b}{\frac {k_{C}q_{i}q_{j}}{r_{ij}}}+{\frac {A}{{r_{\text{OO}}}^{12}}}-{\frac {B}{{r_{\text{OO}}}^{6}}}

によって表わされる。ここで、静電定数k_Cは、分子モデリングにおいて一般的に使用される単位で332.1 Å·kcal/mol の値を持つ^[要出典]。q_i は電子の電荷に対する部分電荷、r_ijは2つの原子または電荷サイトの距離、AおよびBはレナード＝ジョーンズパラメータである。荷電サイトは原子あるいは（孤立電子対といった）ダミーサイト上にある。ほとんどの水モデルにおいて、レナード＝ジョーンズ項は酸素原子間の相互作用にのみ適用される。

下図は3から6サイト水モデルの一般的形状を示している。正確な幾何パラメータ（OHの距離とHOHの角度）はモデルによって異なる。

2-サイト

よく知られた3-サイトSPCモデルに基づく水の2-サイトモデルは、サイトくりこみ（再規格化）分子流体理論を用いて水の誘電特性を予測することが示されている^[4]。

3-サイト

3-サイトモデルは水分子の3つの原子に対応した3つの相互作用点を持つ。それぞれのサイトは点電荷を持ち、酸素原子に対応するサイトはレナード＝ジョーンズパラメータも持つ。3-サイトモデルは高い計算効率を達成しているため、これらは分子動力学シミュレーションの多くの応用に広く用いられている。モデルのほとんどは実際の水分子のものと合致する剛体構造を使用している。例外はSPCモデルであり、このモデルはHOH角について観測値の104.5°ではなく、理想的な正四面体形状（109.47°）を仮定している。

以下の表は一部の3-サイトモデルに対するパラメータを載せている。

	TIPS^[5]	SPC^[6]	TIP3P^[7]	SPC/E^[8]
r(OH), Å	0.9572	1.0	0.9572	1.0
HOH, deg	104.52	109.47	104.52	109.47
A × 10⁻³, kcal Å¹²/mol	580.0	629.4	582.0	629.4
B, kcal Å⁶/mol	525.0	625.5	595.0	625.5
q(O)	−0.80	−0.82	−0.834	−0.8476
q(H)	+0.40	+0.41	+0.417	+0.4238

SPC/Eモデルはポテンシャルエネルギー関数に平均分極補正を加えている。

E_{pol}={\frac {1}{2}}\sum _{i}{\frac {(\mu -\mu ^{0})^{2}}{\alpha _{i}}}

ここで、μは有効に分極した水分子の双極子（SPC/Eモデルでは2.35 D）、 μ⁰は孤立した水分子の双極子モーメント（1.85 D、実験値から）、α_iは等方性分極率定数（1.608 × 10⁻⁴⁰ F m²）である。モデル中の電荷は一定であるため、この補正は全エネルギーに1.25 kcal/mol (5.22 kJ/mol) を加える結果となる。SPC/EモデルはSPCモデルよりも良い密度と拡散定数をもたらす。

CHARMM力場に実装されているTIP3Pモデルは元のモデルからわずかに修正された版である。この差異はレナード＝ジョーンズパラメータにある。元のTIP3Pとは異なり、CHARMM版のTIP3Pは酸素原子上に加えて、2つの水素原子上にもレナード＝ジョーンズパラメータを置いている。電荷は修正されていない^[9]。3-サイトモデル（TIP3P）は比熱の計算に優れた性能を示す^[10]。

柔軟なSPC水モデル

柔軟な単純点電荷水モデル（柔軟なSPC水モデル）は3-サイトSPC水モデルの再パラメータ版である^[11]^[12]。SPCモデルは剛直であるのに対して、「柔軟な」SPCモデルは柔軟である。ToukanとRahmanのモデルにおいて、O-H伸縮は非調和とされ、ゆえに動的挙動がよく記述されている。これは分極を考慮に入れない中では最も正確な3中心水モデルの1つである。分子動力学シミュレーションにおいて、柔軟なSPCモデルは水の正確な密度および誘電率を与える^[13]。

柔軟なSPCモデルはMDynaMixおよびAbaloneに実装されている。

その他のモデル

Ferguson（柔軟なSPC）
CVFF（柔軟）
MG（柔軟、解離性）^[14]
KKYポテンシャル（柔軟なモデル）^[15]
BLXL (smear charged potential)^[16]

4-サイト

4-サイトモデルは、3-サイトモデルのHOH角の二等分線に沿った酸素の近くに1つのダミー原子（図中ではMと表示されている）を加えることによって4つの相互作用点を持つ。ダミー原子は負電荷のみを有している。このモデルは水分子の周りの静電分布を改善している。この手法を用いた最初のモデルは1933年に発表されたBernal-Fowlerモデルであった。このモデルは最も初期の水モデルでもある。しかしながら、BFモデルは密度や蒸発熱といった水のバルク特性をよく再現しない。したがって、歴史的に興味が持たれるのみである。これはパラメータ化の手法の結果である。現代コンピュータが利用可能になった後に開発されたより新しいモデルは、メトロポリス・モンテカルロ法あるいは分子動力学シミュレーションを行い、水のバルク特性を十分に再現するまでパラメータを調節することによってパラメータ化された。

1983年に初めて発表されたTIP4Pモデルは計算化学ソフトウェアパッケージに広く実装されており、生体分子系のシミュレーションにしばしば使用されている。TIP4Pモデルを特殊な用途のために再パラメータ化したモデルも存在する。TIP4P-Ewモデルはエバルトの方法と共に用いるため、TIP3P/Iceは固体の氷のシミュレーションのため、TIP4P/2--5は凝集水の全相図をシミュレーションするためにパラメータ化されている。

4-サイト水モデルのほとんどは、気相中の水分子のものと一致するOH距離とHOH角を用いている。例外はOPCモデルであり、これは水分子の根本的なC_2v分子対称性以外のいかなる構造拘束もかけられていない。代わりに、点電荷とそれらの位置が水分子の静電学を最も良く記述するように最適化されている。OPCは一般的に用いられる剛直n-サイト水モデルよりも正確に包括的な一連のバルク特性を再現する。OPCモデルはAMBER力場に実装されている。

	BF^[17]	TIPS2^[18]	TIP4P^[7]	TIP4P-Ew^[19]	TIP4P/Ice^[20]	TIP4P/2005^[21]	OPC^[22]	TIP4P-D^[23]
r(OH), Å	0.96	0.9572	0.9572	0.9572	0.9572	0.9572	0.8724	0.9572
HOH, deg	105.7	104.52	104.52	104.52	104.52	104.52	103.6	104.52
r(OM), Å	0.15	0.15	0.15	0.125	0.1577	0.1546	0.1594	0.1546
A × 10⁻³, kcal Å¹²/mol	560.4	695.0	600.0	656.1	857.9	731.3	865.1	904.7
B, kcal Å⁶/mol	837.0	600.0	610.0	653.5	850.5	736.0	858.1	900.0
q(M)	−0.98	−1.07	−1.04	−1.04844	−1.1794	−1.1128	-1.3582	-1.16
q(H)	+0.49	+0.535	+0.52	+0.52422	+0.5897	+0.5564	+0.6791	+0.58

その他:

TIP4PF（柔軟）

5-サイト

5-サイトモデルは、四面体様配置を持つ酸素原子の孤立電子対を表現する2つのダミー原子（Lと表示される）上に負電荷を置く。これらの種類の初期のモデルは1971年に提唱されたBen-NaimとStillingerのBNSモデルであり、すぐに1974年のStillingerとRahmanのST2モデルに継承された。主にそれらの高い計算コストのため、5-サイトモデルは、MahoneyとJorgensenのTIP5Pモデルが発表された2000年まではあまり開発されなかった。それ以前のモデルと比較した時、TIP5Pモデルは水二量体の構造（中性子散乱から得られる実験的動径分布関数をよく再現する、より「四面体型」な水構造）と水の密度が最大となる温度という点が改善している。TIP5P-Eモデルはエバルト和と共に用いるためにTIP5Pを再パラメータ化したものである。

	BNS^[24]	ST2^[24]	TIP5P^[25]	TIP5P-E^[26]
r(OH), Å	1.0	1.0	0.9572	0.9572
HOH, deg	109.47	109.47	104.52	104.52
r(OL), Å	1.0	0.8	0.70	0.70
LOL, deg	109.47	109.47	109.47	109.47
A × 10⁻³, kcal Å¹²/mol	77.4	238.7	544.5	554.3
B, kcal Å⁶/mol	153.8	268.9	590.3	628.2
q(L)	−0.19562	−0.2357	−0.241	−0.241
q(H)	+0.19562	+0.2357	+0.241	+0.241
R_L, Å	2.0379	2.0160
R_U, Å	3.1877	3.1287

しかしながら、BNSおよびST2モデルは静電項についてクーロンの法則を直接的に使用しないことに注意する必要がある。これらのモデルは短距離においてスイッチング関数S(r)

S(r_{ij})={\begin{cases}0,&{\mbox{if }}r_{ij}\leq R_{L}\\{\frac {(r_{ij}-R_{L})^{2}(3R_{U}-R_{L}-2r_{ij})}{(R_{U}-R_{L})^{2}}},&{\mbox{if }}R_{L}\leq r_{ij}\leq R_{U}\\1,&{\mbox{if }}R_{U}\leq r_{ij}\end{cases}}

を掛けることで静電項を縮小する。したがって、R_L およびR_UパラメータはBNSおよびST2にのみ適用される。

6-サイト

4-および5-サイトモデルの全てのサイトを混合した6-サイトモデルはNadaとvan der Eerdenによって開発された^[27]。水/氷系の研究のために当初設計されたが、非常に高い融点を持つ^[28]。

その他

生体分子シミュレーションにおける溶質の挙動に対する明示的溶質モデルの効果も広く研究されている。明示的水モデルは折り畳まれていないペプチドの特異的溶媒和と動力学には影響するのに対して、折り畳まれたペプチドのコンホメーションの挙動と柔軟性は損なわれないままであることが明らかにされている^[29]。
MBモデル。メルセデス・ベンツのロゴに似たより抽象的なモデルは2次元系における水の一部の特性を再現する。これは「現実の」（すなわち3次元の）系のシミュレーションには使われないが、定性的研究や教育目的のためには有用である^[30]。
粗視化モデル。水の1-および2-サイトモデルも開発されている^[31]。粗視化モデルでは、それぞれのサイトが複数の水分子を表わすことができる。
多体モデル。量子力学的に解かれたトレーニングセットと、ポテンシャルエネルギー表面を抽出するための機械学習プロトコルを用いて作られた水モデル。これらのポテンシャルエントロピー表面はMDシミュレーションに送り込まれ、凝集相系の物理特性の計算において前例のない精度をもたらした^[32]。

計算コスト

水シミュレーションの計算コストは、水モデルの持つ相互作用サイトの数と共に増加する。CPU時間は計算すべき原子間距離の数におおむね比例する。3-サイトモデルでは、水分子のそれぞれの対について9つの距離が必要である。4-サイトモデルでは10の距離、5-サイトモデルでは17の距離、6-サイトモデルでは26の距離が必要である。

分子動力学法において剛直な水モデルを用いる時は、拘束アルゴリズムを用いて構造を拘束し続けるための追加のコストが存在する（しかし拘束された結合長を用いると、タイムステップの増加が可能になることが多い）。

脚注

^ Skyner, R.; McDonagh, J. L., Groom, C. R., van Mourik, T., Mitchell, J. B. O. (2015). “A Review of Methods for the Calculation of Solution Free Energies and the Modelling of Systems in Solution”. Phys Chem Chem Phys (RSC, PCCP) 17: 6174–91. doi:10.1039/C5CP00288E. PMID 25660403.
^ Allen and Tildesley. (1989). Computer Simulation of Liquids
^ Kirby BJ.. Micro- and Nanoscale Fluid Mechanics: Transport in Microfluidic Devices. http://www.kirbyresearch.com/textbook
^ Dyer, KM; Perkyns, JS; Stell, G; Pettitt (2009). “An angle dependent site-renormalized theory for the conformations of n-butane in a simple fluid”. Mol. Phys. 107: 423–431. PMC 2909681. PMID 20664717.
^ Jorgensen, W. L. Quantum and statistical mechanical studies of liquids. 10. Transferable intermolecular potential functions for water, alcohols, and ethers. Application to liquid water. J. Am. Chem. Soc. 1981, 103, 335-340.
^ H.J.C. Berendsen, J.P.M. Postma, W.F. van Gunsteren, and J. Hermans, In Intermolecular Forces, edited by B. Pullman (Reidel, Dordrecht, 1981), p. 331.
^ ^a ^b Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; Impey, R. W.; Klein (1983). “Comparison of simple potential functions for simulating liquid water”. J. Chem. Phys 79: 926–935. doi:10.1063/1.445869.
^ Berendsen, H. J. C.; Grigera, J. R.; Straatsma, T. P.; The; Pair Potentials, Effective (1987). “The missing term in effective pair potentials”. J. Phys. Chem 91: 6269–6271. doi:10.1021/j100308a038.
^ MacKerell, A. D.; Jr; Bashford, D.; Bellott, R. L.; Dunbrack, R. L.; Jr; Evanseck, J. D.; Field, M. J. et al. (1998). “All-Atom Empirical Potential for Molecular Modeling and Dynamics Studies of Proteins †”. J. Phys. Chem. 102: 3586–3616. doi:10.1021/jp973084f.
^ Mao, Y.; Zhang, Y.; Conductivity, Thermal; Viscosity, Shear; Specific; Water Models, Rigid (2012). “Thermal conductivity, shear viscosity and specific heat of rigid water models”. Chemical Physics Letters 542: 37–41. doi:10.1016/j.cplett.2012.05.044.
^ K. Toukan and A. Rahman (1985). “Molecular-dynamics study of atomic motions in water”. Physical Review B 31: 2643–2648. Bibcode: 1985PhRvB..31.2643T. doi:10.1103/PhysRevB.31.2643.
^ H. J. C. Berendsen, J. R. Grigera, and T. P. Straatsma (1987). “The missing term in effective pair potentials”. Journal of Physical Chemistry 91: 6269–6271. doi:10.1021/j100308a038.
^ M. Praprotnik, D. Janezic and J. Mavri (2004). “Temperature Dependence of Water Vibrational Spectrum: A Molecular Dynamics Simulation Study”. Journal of Physical Chemistry A 108: 11056–11062. doi:10.1021/jp046158d.
^ Mahadevan TS, Garofalini SH (2007). “Dissociative water potential for molecular dynamics simulations”. J. Phys. Chem. B 111 (30): 8919-8927. doi:10.1021/jp072530o. PMID 17604393.
^ Naoki Kumagai, Katsuyuki Kawamura & Toshio Yokokawa (1994). “An Interatomic Potential Model for H2O: Applications to Water and Ice Polymorphs”. Molecular Simulation 12 (3-6). doi:10.1080/08927029408023028.
^ Christian J. Burnham, Jichen Li, Sotiris S. Xantheas and Maurice Leslie (1999). “The parametrization of a Thole-type all-atom polarizable water model from first principles and its application to the study of water clusters (n=2–21) and the phonon spectrum of ice Ih”. J. Chem. Phys. 110: 4566. doi:10.1063/1.478797.
^ Bernal, J. D.; Fowler, R.H. (1933). “A Theory of Water and Ionic Solution, with Particular Reference to Hydrogen and Hydroxyl Ions”. J. Chem. Phys. 1: 515. doi:10.1063/1.1749327.
^ Jorgensen (1982). “Revised TIPS for simulations of liquid water and aqueous solutions”. J. Chem. Phys 77: 4156–4163. doi:10.1063/1.444325.
^ Horn, H. W.; Swope, W. C.; Pitera, J. W.; Madura, J. D.; Dick, T. J.; Hura, G. L.; Head-Gordon, T. (2004). “Development of an improved four-site water model for biomolecular simulations: TIP4P-Ew.”. J. Chem. Phys. 120: 9665–9678. doi:10.1063/1.1683075.
^ Abascal, J. L. F.; Sanz, E.; Fernández, R. García; Vega, C. (2005). “A potential model for the study of ices and amorphous water: TIP4P/Ice”. J. Chem. Phys. 122: 234511. doi:10.1063/1.1931662.
^ Abascal, J. L. F.; Vega, C. (2005). “A general purpose model for the condensed phases of water: TIP4P/2005”. J. Chem. Phys. 123: 234505. doi:10.1063/1.2121687.
^ Izadi, S.; Anandakrishnan, R.; Onufriev, A. V. (2014). “Building Water Models: A Different Approach”. The Journal of Physical Chemistry Letters 5 (21): 3863–3871. doi:10.1021/jz501780a.
^ Piana, Stefano; Donchev, Alexander G.; Robustelli, Paul; Shaw, David E. (2015). “Water Dispersion Interactions Strongly Influence Simulated Structural Properties of Disordered Protein States”. The Journal of Physical Chemistry B 119: 150312133835006. doi:10.1021/jp508971m. ISSN 1520-6106.
^ ^a ^b Stillinger, F.H.; Rahman, A.; Improved (1974). “Improved simulation of liquid water by molecular dynamics”. J. Chem. Phys. 60: 1545–1557. doi:10.1063/1.1681229.
^ Mahoney, M. W.; Jorgensen (2000). “A five-site model for liquid water and the reproduction of the density anomaly by rigid, nonpolarizable potential functions”. J. Chem. Phys. 112: 8910–8922. doi:10.1063/1.481505.
^ Rick, S. W. (2004). “A reoptimization of the five-site water potential (TIP5P) for use with Ewald sums”. J. Chem. Phys. 120: 6085–6093. doi:10.1063/1.1652434.
^ Nada, H. (2003). “An intermolecular potential model for the simulation of ice and water near the melting point: A six-site model of H[sub 2]O”. J. Chem. Phys. 118: 7401. doi:10.1063/1.1562610.
^ José L. F. Abascal, Ramón García Fernández, Carlos Vega1 and Marcelo A. Carignano (2006). “The melting temperature of the six site potential model of water”. J. Chem. Phys. 125: 166101. doi:10.1063/1.2360276.
^ Florova, P.; Sklenovsky, P.; Banas, P.; Otyepka, M. (2010). “Explicit Water Models Affect the Specific Solvation and Dynamics of Unfolded Peptides While the Conformational Behavior and Flexibility of Folded Peptides Remain Intact”. J. Chem. Theory Comput 6: 3569–3579. doi:10.1021/ct1003687.
^ Silverstein, K. A. T.; Haymet, A. D. J.; Dill, K. A.; Effect, Hydrophobic (1998). “A Simple Model of Water and the Hydrophobic Effect”. J. Am. Chem. Soc. 120: 3166–3175. doi:10.1021/ja973029k.
^ Izvekov, S.; Voth, G. A.; Multiscale. “J. Chem. Phys. '”. 2005' 123: 134105. doi:10.1063/1.2038787.
^ Medders, G.R.; Paesani, F. (2015). “Infrared and Raman Spectroscopy of Liquid Water through "First Principles" Many-Body Molecular Dynamics”. J. Chem. Theory Comput 11: 1145–1154. doi:10.1021/ct501131j.