熱力学的積分法

熱力学的積分法（ねつりきがくてきせきぶんほう、、TI法）は、ポテンシャルエネルギーが空間座標に異なる依存性を持つ2つの任意の状態（ポテンシャルエネルギーformula_1およびformula_2を持つ状態AおよびB）間の自由エネルギーにおける差を比較するために用いられる手法である。系の自由エネルギーは単に系の位相空間座標の関数ではなく、位相空間上のボルツマン重みをかけた積分の関数（すなわち分配関数）であるため、2つの状態間の自由エネルギー差は2つの座標セット（それぞれ状態AおよびBのもの）のポテンシャルエネルギーから直接的には計算できない。熱力学的積分法では、自由エネルギー差は状態間の熱力学的経路を定義し、この経路に沿ったアンサンブル平均エンタルピー変化を積分することによって計算される。こういった経路は実際の化学的（chemical）な過程でも錬金術的（alchemical）な過程でもよい。錬金術的（非科学的）過程の一例は、Kirkwoodの結合パラメーター法である。
導出.
ポテンシャルエネルギーformula_1およびformula_2を持つ状態AおよびBを考える。いずれの系におけるポテンシャルエネルギーも、分子動力学法あるいは適切なボルツマン重みを持つモンテカルロシミュレーションからサンプリングされた配置のアンサンブル平均として計算できる。ここで、
として定義される新たなポテンシャルエネルギー関数を考える。上式において、formula_6は0と1の間の値を持つ結合パラメーターとして定義され、したがってformula_6の関数としてのポテンシャルエネルギはformula_8の系Aのエネルギーおよびformula_9の系Bのエネルギーとは異なる。カノニカルアンサンブルにおいて、系の分配関数は以下の式で書くことができる。
この表記法では、 formula_11はポテンシャルエネルギー関数formula_12を持つアンサンブルにおける状態formula_13のポテンシャルエネルギーである。この系の自由エネルギーは以下の式で定義される。
λに関してFの導関数を取ると、これがλに関するポテンシャルエネルギーの導関数のアンサンブル平均と等しいことが分かる。
したがって、状態AとBとの間の自由エネルギー差の変化は、カップリングパラメータformula_6上でのポテンシャルエネルギーのアンサンブル平均導関数の積分から計算することができる。実際のところは、これはポテンシャルエネルギー関数formula_12を定義し、一連のformula_6値における平均配置のアンサンブルをサンプリングし、個々のformula_6値においてformula_6に関するformula_12のアンサンブル平均導関数を計算し、最後にアンサンブル平均導関数上の積分を計算することによって行われる。
アンブレラサンプリング法は関連した自由エネルギー計算手法である。アンブレラサンプリング法はポテンシャルエネルギーにバイアスを加える。無限に強いバイアスの限界においては、これは熱力学的積分法と等価である。