水モデルのソースを表示

[[ファイル:Water-2D-labelled.png|thumb|水モデルはその幾何構造と原子電荷やレナード＝ジョーンズパラメータといったその他のパラメータによって定義される。]]

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

露な（明示的な）水モデルの代わりとしては[[陰溶媒|非明示的溶媒和]]モデル（連続体モデルとも呼ばれる）がある。非明示的溶媒和モデルの例としては、[[COSMO溶媒和モデル]]や[[分極連続体モデル]]（PCM）、ハイブリッド溶媒和モデルがある<ref name = "Skyner et al">{{cite journal | last = Skyner | first = R. | coauthors = McDonagh, J. L., Groom, C. R., van Mourik, T., Mitchell, J. B. O.| title = A Review of Methods for the Calculation of Solution Free Energies and the Modelling of Systems in Solution  | publisher = RSC, PCCP | year = 2015 | doi = 10.1039/C5CP00288E | pmid=25660403 | volume=17 | journal=Phys Chem Chem Phys | pages=6174–91}}</ref>。

==単純な水モデル==
剛直なモデルは[[分子間力|非結合相互作用]]に依存する最も単純な水モデルとして知られる。これらのモデルにおいて、結合性相互作用は[[ホロノミック拘束]]によって陰に扱われる。静電相互作用は[[クーロンの法則]]を用いてモデル化され、分散力および斥力は[[レナード-ジョーンズ・ポテンシャル]]を用いてモデル化される<ref name=Allen>{{cite book | author=Allen and Tildesley. | title=Computer Simulation of Liquids | 
year=1989}}</ref><ref name=Kirby>{{cite book | author=Kirby BJ. | title=Micro- and Nanoscale Fluid Mechanics: Transport in Microfluidic Devices | url=http://www.kirbyresearch.com/textbook}}</ref>。TIP3PおよびTIP4Pといったモデルに対するポテンシャルは
:<math> E_{ab} = \sum_{i} ^{\text{on }a} \sum_{j} ^{\text{on }b}
  \frac {k_Cq_iq_j}{r_{ij}}
  + \frac {A}{{r_{\text{OO}}}^{12}}
  - \frac {B}{{r_{\text{OO}}}^6}
</math>
によって表わされる。ここで、[[クーロン定数|静電定数]]''k<sub>C</sub>''は、分子モデリングにおいて一般的に使用される単位で332.1 Å·kcal/mol の値を持つ{{Citation needed|date=November 2010}}。''q<sub>i</sub>'' は電子の電荷に対する[[部分電荷]]、''r<sub>ij</sub>''は2つの原子または電荷サイトの距離、''A''および''B''はレナード＝ジョーンズパラメータである。荷電サイトは原子あるいは（[[孤立電子対]]といった）ダミーサイト上にある。ほとんどの水モデルにおいて、レナード＝ジョーンズ項は酸素原子間の相互作用にのみ適用される。

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

[[Image:Water models.svg|700px]]

==2-サイト==
よく知られた3-サイトSPCモデルに基づく水の2-サイトモデルは、サイトくりこみ（再規格化）分子流体理論を用いて水の誘電特性を予測することが示されている<ref>{{cite journal | last1 = Dyer | first1 = KM | last2 = Perkyns | first2 = JS | last3 = Stell | first3 = G | last4 = Pettitt | first4 =  | year = 2009   | url = | journal = Mol. Phys. | volume = 107 | issue = | pages = 423–431 |title=An angle dependent site-renormalized theory for the conformations of n-butane in a simple fluid|pmc=2909681|pmid= 20664717}}</ref>。

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

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

{| class="wikitable"
! 
! TIPS<ref>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.</ref> 
! SPC<ref>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.</ref> 
! TIP3P<ref name="tip3p">{{cite journal | last1 = Jorgensen | first1 = W. L. | last2 = Chandrasekhar | first2 = J. | last3 = Madura | first3 = J. D. | last4 = Impey | first4 = R. W. | last5 = Klein | first5 =  | year = 1983 | title =  Comparison of simple potential functions for simulating liquid water| url = | journal = J. Chem. Phys | volume = 79 | issue = | pages = 926–935 | doi = 10.1063/1.445869 }}</ref> 
! SPC/E<ref>{{cite journal | last1 = Berendsen | first1 = H. J. C. | last2 = Grigera | first2 = J. R. | last3 = Straatsma | first3 = T. P. | last4 = The | first4 =  | last5 = Pair Potentials | first5 = Effective | year = 1987 | title =  The missing term in effective pair potentials| url = | journal = J. Phys. Chem | volume = 91 | issue = | pages = 6269–6271 | doi = 10.1021/j100308a038 }}</ref>
|-
!r(OH), Å
| 0.9572 || 1.0 || 0.9572 || 1.0
|-
!HOH, deg
| 104.52 || 109.47 || 104.52 || 109.47
|-
! A × 10<sup>−3</sup>, kcal Å<sup>12</sup>/mol
| 580.0 || 629.4 || 582.0 || 629.4
|-
! B, kcal Å<sup>6</sup>/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モデルはポテンシャルエネルギー関数に平均分極補正を加えている。
:<math>E_{pol} = \frac 1 2 \sum_{i}
  \frac {(\mu - \mu^0)^2}{\alpha_i}
</math>
ここで、μは有効に分極した水分子の双極子（SPC/Eモデルでは2.35 [[デバイ|D]]）、 μ<sup>0</sup>は孤立した水分子の双極子モーメント（1.85 D、実験値から）、α<sub>i</sub>は等方性分極率定数（1.608 × 10<sup>−40</sup> F m<sup>2</sup>）である。
モデル中の電荷は一定であるため、この補正は全エネルギーに1.25 kcal/mol (5.22 kJ/mol) を加える結果となる。SPC/EモデルはSPCモデルよりも良い密度と拡散定数をもたらす。

[[CHARMM]]力場に実装されているTIP3Pモデルは元のモデルからわずかに修正された版である。この差異はレナード＝ジョーンズパラメータにある。元のTIP3Pとは異なり、CHARMM版のTIP3Pは酸素原子上に加えて、2つの水素原子上にもレナード＝ジョーンズパラメータを置いている。電荷は修正されていない<ref>{{cite journal | last1 = MacKerell | first1 = A. D. | last2 = Jr | first2 =  | last3 = Bashford | first3 = D. | last4 = Bellott | first4 = R. L. | last5 = Dunbrack | first5 = R. L. | last6 = Jr | first6 =  | last7 = Evanseck | first7 = J. D. | last8 = Field | first8 = M. J. | last9 = Fischer | first9 = S. | last10 = Gao | first10 = J. | last11 = Guo | first11 = H. | last12 = Ha | first12 = S. | last13 = Joseph-McCarthy | first13 = D. | last14 = Kuchnir | first14 = L. | last15 = Kuczera | first15 = K. | last16 = Lau | first16 = F. T. K. | last17 = Mattos | first17 = C. | last18 = Michnick | first18 = S. | last19 = Ngo | first19 = T. | last20 = Nguyen | first20 = D. T. | last21 = Prodhom | first21 = B. | last22 = Reiher | first22 = W. E. | last23 = III | first23 =  | last24 = Roux | first24 = B. | last25 = Schlenkrich | first25 = M. | last26 = Smith | first26 = J. C. | last27 = Stote | first27 = R. | last28 = Straub | first28 = J. | last29 = Watanabe | first29 = M. | last30 = Wiorkiewicz-Kuczera | first30 = J. | last31 = Yin | first31 = D. | last32 = Karplus | first32 =  | last33 = Modeling | first33 = Molecular | last34 = Dynamics | first34 =  | year = 1998 | title =  All-Atom Empirical Potential for Molecular Modeling and Dynamics Studies of Proteins †| url = | journal = J. Phys. Chem. | volume = 102 | issue = | pages = 3586–3616 | doi = 10.1021/jp973084f }}</ref>。3-サイトモデル（TIP3P）は比熱の計算に優れた性能を示す<ref>{{cite journal | last1 = Mao | first1 = Y. | last2 = Zhang | first2 = Y. | last3 = Conductivity | first3 = Thermal | last4 = Viscosity | first4 = Shear | last5 = Specific | first5 =  | last6 = Water Models | first6 = Rigid | year = 2012 | title =  Thermal conductivity, shear viscosity and specific heat of rigid water models| url = | journal = Chemical Physics Letters | volume = 542 | issue = | pages = 37–41 | doi = 10.1016/j.cplett.2012.05.044 }}</ref>。

=== 柔軟なSPC水モデル===
[[ファイル:FlexibleSPC-E.png|thumb|柔軟なSPC水モデル]]

柔軟な単純点電荷水モデル（柔軟なSPC水モデル）は3-サイトSPC水モデルの再パラメータ版である<ref>
{{cite journal
 |author = K. Toukan and A. Rahman
 |year= 1985
 |title = Molecular-dynamics study of atomic motions in water
 |journal = [[Physical Review B]]
 |volume = 31 |pages = 2643–2648
 |doi=10.1103/PhysRevB.31.2643
|bibcode = 1985PhRvB..31.2643T }}</ref><ref>
{{cite journal
 |author = H. J. C. Berendsen, J. R. Grigera, and T. P. Straatsma
 |year=1987
 |title = The missing term in effective pair potentials
 |journal = Journal of Physical Chemistry
 |volume = 91 |pages = 6269–6271
 |doi=10.1021/j100308a038
}}</ref>。SPCモデルは剛直であるのに対して、「柔軟な」SPCモデルは柔軟である。ToukanとRahmanのモデルにおいて、O-H伸縮は非調和とされ、ゆえに動的挙動がよく記述されている。これは[[分極電荷密度|分極]]を考慮に入れない中では最も正確な3中心水モデルの1つである。分子動力学シミュレーションにおいて、柔軟なSPCモデルは水の正確な[[密度]]および[[誘電率]]を与える<ref>
{{cite journal
 |author = M. Praprotnik, D. Janezic and J. Mavri
 |year= 2004
 |title = Temperature Dependence of Water Vibrational Spectrum: A Molecular Dynamics Simulation Study
 |journal = Journal of Physical Chemistry A
 |volume = 108 |pages = 11056–11062
 |doi=10.1021/jp046158d
}}</ref>。

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

===その他のモデル===
*Ferguson（柔軟なSPC）
*CVFF（柔軟）
*MG（柔軟、解離性）<ref>{{cite journal|journal=J. Phys. Chem. B|year= 2007|volume=111|issue=30|pages=8919-8927|title=Dissociative water potential for molecular dynamics simulations|author=Mahadevan TS, Garofalini SH|pmid= 17604393 |doi=10.1021/jp072530o|url=http://glass.rutgers.edu/dir.papers/mahaJPCB07-H2Opot.pdf}}</ref>
*KKYポテンシャル（柔軟なモデル）<ref>{{cite journal|journal=Molecular Simulation|volume=12|issue= 3-6|year= 1994|title=An Interatomic Potential Model for H2O: Applications to Water and Ice Polymorphs|doi=10.1080/08927029408023028|author=Naoki Kumagai, Katsuyuki Kawamura & Toshio Yokokawa}}</ref>
*BLXL (smear charged potential)<ref>{{cite journal|title=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|author=Christian J. Burnham, Jichen Li, Sotiris S. Xantheas and Maurice Leslie||journal=J. Chem. Phys. |volume=110|pages= 4566 |year=1999|doi=10.1063/1.478797}}</ref>

==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<sub>2v</sub>[[分子対称性]]以外のいかなる構造拘束もかけられていない。代わりに、点電荷とそれらの位置が水分子の静電学を最も良く記述するように最適化されている。OPCは一般的に用いられる剛直n-サイト水モデルよりも正確に包括的な一連のバルク特性を再現する。OPCモデルは[[AMBER (分子動力学)|AMBER]]力場に実装されている。

{|class="wikitable"
! 
! BF<ref>{{cite journal | last1 = Bernal | first1 = J. D. | last2 = Fowler | first2 = R.H. | year = 1933 | title =  A Theory of Water and Ionic Solution, with Particular Reference to Hydrogen and Hydroxyl Ions| url = | journal = J. Chem. Phys. | volume = 1 | issue = | page = 515 | doi = 10.1063/1.1749327 }}</ref>
! TIPS2<ref>{{cite journal | last1 = Jorgensen | first1 =  | year = 1982 | title =  Revised TIPS for simulations of liquid water and aqueous solutions| url = | journal = J. Chem. Phys | volume = 77 | issue = | pages = 4156–4163 | doi = 10.1063/1.444325 }}</ref>
! TIP4P<ref name="tip3p" />
! TIP4P-Ew<ref>{{cite journal | last1 = Horn | first1 = H. W. | last2 = Swope | first2 = W. C. | last3 = Pitera | first3 = J. W. | last4 = Madura | first4 = J. D. | last5 = Dick | first5 = T. J. | last6 = Hura | first6 = G. L. | last7 = Head-Gordon | first7 = T. | year = 2004 | title = Development of an improved four-site water model for biomolecular simulations: TIP4P-Ew. | url = | journal = J. Chem. Phys. | volume = 120 | issue = | pages = 9665–9678 | doi = 10.1063/1.1683075 }}</ref>
! TIP4P/Ice<ref>{{cite journal | last1 = Abascal | first1 = J. L. F. | last2 = Sanz | first2 = E. | last3 = Fernández | first3 = R. García | last4 = Vega | first4 = C. | year = 2005 | title = A potential model for the study of ices and amorphous water: TIP4P/Ice | url = | journal = J. Chem. Phys. | volume = 122 | issue = | page = 234511 | doi = 10.1063/1.1931662 }}</ref>
! TIP4P/2005<ref>{{cite journal | last1 = Abascal | first1 = J. L. F. | last2 = Vega | first2 = C. | year = 2005 | title = A general purpose model for the condensed phases of water: TIP4P/2005 | url = | journal = J. Chem. Phys. | volume = 123 | issue = | page = 234505 | doi = 10.1063/1.2121687 }}</ref>
! OPC<ref>{{cite journal | last1 = Izadi | first1 = S. | last2 = Anandakrishnan | first2 = R. | last3 = Onufriev | first3 = A. V. | year = 2014 | title = Building Water Models: A Different Approach | url = | journal = The Journal of Physical Chemistry Letters | volume = 5 | issue = 21| pages = 3863–3871 | doi = 10.1021/jz501780a }}</ref>
! TIP4P-D<ref>{{cite journal|last1=Piana|first1=Stefano|last2=Donchev|first2=Alexander G.|last3=Robustelli|first3=Paul|last4=Shaw|first4=David E.|title=Water Dispersion Interactions Strongly Influence Simulated Structural Properties of Disordered Protein States|journal=The Journal of Physical Chemistry B|year=2015|pages=150312133835006|issn=1520-6106|doi=10.1021/jp508971m|volume=119}}</ref>
|-
!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<sup>−3</sup>, kcal Å<sup>12</sup>/mol
| 560.4 || 695.0 || 600.0 || 656.1 || 857.9 || 731.3||865.1||904.7
|-
! B, kcal Å<sup>6</sup>/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を再パラメータ化したものである。

{|class="wikitable"
!
! BNS<ref name="st2">{{cite journal | last1 = Stillinger | first1 = F.H. | last2 = Rahman | first2 = A. | last3 = Improved | first3 =  | year = 1974 | title =  Improved simulation of liquid water by molecular dynamics| url = | journal = J. Chem. Phys. | volume = 60 | issue = | pages = 1545–1557 | doi = 10.1063/1.1681229 }}</ref>
! ST2<ref name="st2" />
! TIP5P<ref>{{cite journal | last1 = Mahoney | first1 = M. W. | last2 = Jorgensen | first2 =  | year = 2000 | title =  A five-site model for liquid water and the reproduction of the density anomaly by rigid, nonpolarizable potential functions| url = | journal = J. Chem. Phys. | volume = 112 | issue = | pages = 8910–8922 | doi = 10.1063/1.481505 }}</ref>
! TIP5P-E<ref>{{cite journal | last1 = Rick | first1 = S. W. | year = 2004 | title = A reoptimization of the five-site water potential (TIP5P) for use with Ewald sums | url = | journal = J. Chem. Phys. | volume = 120 | issue = | pages = 6085–6093 | doi = 10.1063/1.1652434 }}</ref>
|-
!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<sup>−3</sup>, kcal Å<sup>12</sup>/mol
| 77.4 || 238.7 || 544.5 || 554.3
|-
! B, kcal Å<sup>6</sup>/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<sub>L</sub>'', Å
| 2.0379 || 2.0160
|-
! ''R<sub>U</sub>'', Å
| 3.1877 || 3.1287
|}

しかしながら、BNSおよびST2モデルは静電項についてクーロンの法則を直接的に使用しないことに注意する必要がある。これらのモデルは短距離においてスイッチング関数''S(r)''
:<math>
S(r_{ij}) = 
\begin{cases} 
  0,    & \mbox{if }r_{ij} \le R_L \\
  \frac{(r_{ij} - R_L)^2(3R_U - R_L - 2r_{ij})}{(R_U - R_L)^2},  & \mbox{if }R_L \le r_{ij} \le R_U \\
  1, & \mbox{if }R_U \le r_{ij}
\end{cases}
</math>
を掛けることで静電項を縮小する。したがって、''R<sub>L</sub>'' および''R<sub>U</sub>''パラメータはBNSおよびST2にのみ適用される。

==6-サイト==
4-および5-サイトモデルの全てのサイトを混合した6-サイトモデルはNadaとvan der Eerdenによって開発された<ref>{{cite journal | last1 = Nada | first1 = H. | year = 2003 | title =  An intermolecular potential model for the simulation of ice and water near the melting point: A six-site model of H[sub 2]O| url = | journal = J. Chem. Phys. | volume = 118 | issue = | page = 7401 | doi = 10.1063/1.1562610 }}</ref>。水/氷系の研究のために当初設計されたが、非常に高い融点を持つ<ref>{{cite journal|title=The melting temperature of the six site potential model of water|author=José L. F. Abascal, Ramón García Fernández, Carlos Vega1 and Marcelo A. Carignano||journal=J. Chem. Phys.|volume= 125|pages= 166101 |year=2006|doi=10.1063/1.2360276|url=http://cacharro.quim.ucm.es/publicaciones/JCP_125_166101_2006.pdf}}</ref>。

==その他==
*生体分子シミュレーションにおける溶質の挙動に対する明示的溶質モデルの効果も広く研究されている。明示的水モデルは折り畳まれていないペプチドの特異的溶媒和と動力学には影響するのに対して、折り畳まれたペプチドのコンホメーションの挙動と柔軟性は損なわれないままであることが明らかにされている<ref>{{cite journal | last1 = Florova | first1 = P. | last2 = Sklenovsky | first2 = P. | last3 = Banas | first3 = P. | last4 = Otyepka | first4 = M. | year = 2010 | title =  Explicit Water Models Affect the Specific Solvation and Dynamics of Unfolded Peptides While the Conformational Behavior and Flexibility of Folded Peptides Remain Intact| url = | journal = J. Chem. Theory Comput | volume = 6 | issue = | pages = 3569–3579 | doi = 10.1021/ct1003687 }}</ref>。
*MBモデル。[[メルセデス・ベンツ]]のロゴに似たより抽象的なモデルは2次元系における水の一部の特性を再現する。これは「現実の」（すなわち3次元の）系のシミュレーションには使われないが、定性的研究 や教育目的のためには有用である<ref>{{cite journal | last1 = Silverstein | first1 = K. A. T. | last2 = Haymet | first2 = A. D. J. | last3 = Dill | first3 = K. A. | last4 = Effect | first4 = Hydrophobic | year = 1998 | title =  A Simple Model of Water and the Hydrophobic Effect| url = | journal = J. Am. Chem. Soc. | volume = 120 | issue = | pages = 3166–3175 | doi = 10.1021/ja973029k }}</ref>。
*粗視化モデル。水の1-および2-サイトモデルも開発されている<ref>{{cite journal | last1 = Izvekov | first1 = S. | last2 = Voth | first2 = G. A. | last3 = Multiscale | first3 =  | year = | title = J. Chem. Phys. ' | url = | journal = 2005' | volume = 123 | issue = | page = 134105 | doi = 10.1063/1.2038787 }}</ref>。粗視化モデルでは、それぞれのサイトが複数の水分子を表わすことができる。
*多体モデル。量子力学的に解かれたトレーニングセットと、ポテンシャルエネルギー表面を抽出するための機械学習プロトコルを用いて作られた水モデル。これらのポテンシャルエントロピー表面はMDシミュレーションに送り込まれ、凝集相系の物理特性の計算において前例のない精度をもたらした<ref>{{cite journal | last1 = Medders | first1 = G.R. | last2 = Paesani | first2 = F. | year = 2015 | title = Infrared and Raman Spectroscopy of Liquid Water through "First Principles" Many-Body Molecular Dynamics | url = | journal = J. Chem. Theory Comput | volume = 11 | issue = | pages = 1145–1154 | doi = 10.1021/ct501131j }}</ref>。

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

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

==脚注==
{{reflist}}

==関連項目==
*[[力場 (化学)]]
*[[分子力学法]]
*[[分子動力学法]]
*[[溶媒モデル]]

{{DEFAULTSORT:みすもてるl}}
[[Category:水]]
[[Category:計算化学]]
[[Category:種類別のモデル]]