アンブレラサンプリング法

アンブレラサンプリング法（アンブレラサンプリングほう 、傘サンプリングとも）は、計算物理学および計算化学の手法の一つであり、エネルギー地形の形状によってが妨げられる系のサンプリングを改善するために用いられる。1977年に、TorrieおよびValleauによって初めて提案された。アンブレラサンプリング法は統計学においてより一般的なの特別な物理学的応用である。

エネルギー障壁が配置空間の2つの領域を隔てている系では、不十分なサンプリングに悩まされる。メトロポリス・モンテカルロシミュレーションでは、ポテンシャル障壁を乗り越える低い確率は、シミュレーションによって不十分にサンプリングされる、または全くサンプリングされない到達しにくい配置が残る可能性がある。容易に思い浮かべられる例は融点における固体で起こる。秩序パラメータ  を持つ系の状態を考えると、液相（低い）と固相 （高い ）のどちらもエネルギー的に低いが、 の中間の値での自由エネルギー障壁によって隔てられている。これによって両方の相の十分なサンプリングをシミュレーションで行うことができない。

アンブレラサンプリング法はこの状況における困難を乗り越える方法の一つである。この方法では、モンテカルロサンプリングで標準的に用いられるボルツマン重み付けを、エネルギー障壁の影響を打ち消すよう選ばれた重みに置き換える。生成されたマルコフ連鎖は、以下の式で与えられる確率分布を持つ。

上式において、 は真のポテンシャルエネルギー、 はボルツマン重み付けモンテカルロシミュレーションでは近づき難い配置を促進するよう選ばれた関数である。上記の例において、 は  が中間の  で高い値、低い/高い  で低い値を取り、障壁の乗り越えを容易にするように選ばれる。

このようなやり方で行われるサンプリング実行から推定される熱力学特性  は以下の式によってカノニカルアンサンブル値へと変換できる。

添字  はアンブレラサンプリングされたシミュレーションからの値であることを示している。

重み関数  を導入することの効果は系のポテンシャルエネルギーにバイアスポテンシャル  を加えるのと等価である。

バイアスポテンシャルが厳密に反応座標または秩序パラメータ  の関数ならば、次に反応座標の（バイアスのない）自由エネルギープロファイルはバイアスを受けた自由エネルギープロファイルからバイアスポテンシャルを差し引くことで計算できる。

上式において、  はバイアスのない系の自由エネルギープロファイル、  はバイアスを受けた（アンブレラサンプリングされた）系について計算された自由エネルギープロファイルである。

一連のアンブレラサンプリングシミュレーションは重み付きヒストグラム解析法（）あるいはその一般化を用いて解析できる。WHAMは最尤法を使って導くことができる。

平均力ポテンシャルまたは反応速度を計算するためのアンブレラサンプリング法に代わる方法としては、自由エネルギー摂動法とがある。完全な非平衡において機能する別の方法としてS-PRES（Stochastic process rare event sampling）法がある。

脚注

関連項目 
 熱力学的積分法

モンテカルロ法
分子動力学
計算化学
理論化学