熱力学的積分法のソースを表示

'''熱力学的積分法'''（ねつりきがくてきせきぶんほう、{{lang-en-short|Thermodynamic integration}}、TI法）は、ポテンシャルエネルギーが空間座標に異なる依存性を持つ2つの任意の状態（ポテンシャルエネルギー<math>U_A</math>および<math>U_B</math>を持つ状態AおよびB）間の[[自由エネルギー]]における差を比較するために用いられる手法である。系の自由エネルギーは単に系の位相空間座標の関数ではなく、位相空間上の[[ボルツマン分布|ボルツマン重み]]をかけた積分の関数（すなわち[[分配関数]]）であるため、2つの状態間の自由エネルギー差は2つの座標セット（それぞれ状態AおよびBのもの）のポテンシャルエネルギーから直接的には計算できない。熱力学的積分法では、自由エネルギー差は状態間の熱力学的経路を定義し、この経路に沿ったアンサンブル平均エンタルピー変化を積分することによって計算される。こういった経路は実際の化学的（chemical）な過程でも錬金術的（alchemical）な過程でもよい。錬金術的（非科学的）過程の一例は、Kirkwoodの結合パラメーター法である<ref>{{Cite journal |doi = 10.1063/1.1749657|bibcode = 1935JChPh...3..300K|title = Statistical Mechanics of Fluid Mixtures|year = 1935|last1 = Kirkwood|first1 = John G.|journal = The Journal of Chemical Physics|volume = 3|issue = 5|pages = 300–313}}</ref>。

== 導出 ==
ポテンシャルエネルギー<math>U_A</math>および<math>U_B</math>を持つ状態AおよびBを考える。いずれの系におけるポテンシャルエネルギーも、[[分子動力学]]法あるいは適切なボルツマン重みを持つ[[モンテカルロ法|モンテカルロ]]シミュレーションからサンプリングされた配置のアンサンブル平均として計算できる。ここで、
:<math>U(\lambda) = U_A + \lambda(U_B - U_A)</math>
として定義される新たなポテンシャルエネルギー関数を考える。上式において、<math>\lambda</math>は0と1の間の値を持つ結合パラメーターとして定義され、したがって<math>\lambda</math>の関数としてのポテンシャルエネルギは<math>\lambda = 0</math>の系Aのエネルギーおよび<math>\lambda = 1</math>の系Bのエネルギーとは異なる。[[カノニカル分布|カノニカルアンサンブル]]において、系の分配関数は以下の式で書くことができる。
:<math>Q(N, V, T, \lambda) = \sum_{s} \exp [-U_s(\lambda)/k_{B}T]</math>

この表記法では、 <math>U_s(\lambda)</math>はポテンシャルエネルギー関数<math>U(\lambda)</math>を持つアンサンブルにおける状態<math>s</math>のポテンシャルエネルギーである。この系の自由エネルギーは以下の式で定義される。
:<math>F(N,V,T,\lambda)=-k_{B}T \ln Q(N,V,T,\lambda)</math>,

λに関してFの導関数を取ると、これがλに関するポテンシャルエネルギーの導関数のアンサンブル平均と等しいことが分かる。
:<math>\begin{align}
\Delta F(A \rightarrow B)
&= \int_0^1 \frac{\partial F(\lambda)}{\partial\lambda} d\lambda \\
&= -\int_0^1 \frac{k_{B}T}{Q} \frac{\partial Q}{\partial\lambda} d\lambda \\
&= \int_0^1 \frac{k_{B}T}{Q} \sum_{s} \frac{1}{k_{B}T} \exp[- U_s(\lambda)/k_{B}T ] \frac{\partial U_s(\lambda)}{\partial \lambda} d\lambda \\
&= \int_0^1 \left\langle\frac{\partial U(\lambda)}{\partial\lambda}\right\rangle_{\lambda} d\lambda \\
&= \int_0^1 \left\langle U_B(\lambda) - U_A(\lambda) \right\rangle_{\lambda} d\lambda
\end{align}
</math>

したがって、状態AとBとの間の自由エネルギー差の変化は、カップリングパラメータ<math>\lambda</math>上でのポテンシャルエネルギーのアンサンブル平均導関数の積分から計算することができる<ref>Frenkel, Daan and Smit, Berend. ''Understanding Molecular Simulation: From Algorithms to Applications.'' Academic Press, 2007</ref>。実際のところは、これはポテンシャルエネルギー関数<math>U(\lambda)</math>を定義し、一連の<math>\lambda</math>値における平均配置のアンサンブルをサンプリングし、個々の<math>\lambda</math>値において<math>\lambda</math>に関する<math>U(\lambda)</math>のアンサンブル平均導関数を計算し、最後にアンサンブル平均導関数上の積分を計算することによって行われる。

[[アンブレラサンプリング法]]は関連した自由エネルギー計算手法である。アンブレラサンプリング法はポテンシャルエネルギーにバイアスを加える。無限に強いバイアスの限界においては、これは熱力学的積分法と等価である<ref>{{cite journal | doi = 10.1021/ct050252w | pmid = 26626532| author = J Kästner| year = 2006 | title = QM/MM Free-Energy Perturbation Compared to Thermodynamic Integration and Umbrella Sampling: Application to an Enzymatic Reaction | journal = Journal of Chemical Theory and Computation | volume = 2 | issue = 2 | pages = 452–461 | url =  | format =  | accessdate = |display-authors=etal}}</ref>。

==出典==
<references/>

== 補足資料 ==
* {{cite journal|和書|journal=アンサンブル|title=自由エネルギー計算の理論: 創薬応用を実現する定量的予測への挑戦|author=山下 雄史|year=2015|volume= 17|issue= 2|pages= 83-91|doi=10.11436/mssj.17.83 }}

==関連項目==
* [[自由エネルギー摂動法]]
* [[ベネット受容比法]]
* [[クルックスの揺動定理]]

{{DEFAULTSORT:ねつりきかくてきせきふんほう}}
[[Category:計算化学]]
[[Category:統計力学]]