Chemistry Labs

理論・計算化学

ポテンシャルエネルギー面と自由エネルギー計算

分子幾何とポテンシャルエネルギー地形を結びつけ、統計力学と拡張サンプリングで自由エネルギー差や障壁を推定します。

1 中学2 高校3 大学4 上級5 研究

直観エネルギー地形と到達可能な状態

ポテンシャルエネルギー面は各分子構造にエネルギーを割り当てます。傾きは力を示し、谷は安定構造を示唆しますが、温度によって統計的に占有される領域は変わります。平衡確率を決めるのはポテンシャルだけでなく自由エネルギーです。

分子エネルギー地形と射影自由エネルギープロファイルを比較します。低エネルギー経路も速度論的に到達困難な場合があり、自由エネルギープロファイルは選んだ集団変数と十分なサンプリングに依存します。

学校エネルギー図の読み方

定義: ポテンシャルエネルギー面

電子状態と核座標 R を固定したとき V(R) がポテンシャルエネルギーです。停留点では勾配がゼロです。局所極小は拘束されない全方向で正の曲率を持ち、一次鞍点には下り方向が一つあります。

F(R)=−∇RV(R)F(R) = −∇_R V(R)

例: 反応座標上の障壁

単純化した一次元座標で、反応物の極小は遷移状態の極大より12 kJ mol⁻¹低い。反応物から遷移状態までのポテンシャル障壁はいくらですか。

解答

障壁は V‡−V_反応物=12 kJ mol⁻¹ です。これはポテンシャルエネルギー差であり、有限温度での活性化自由エネルギーと自動的に同じではありません。

大学微視状態から自由エネルギーへ

定義: 平均力ポテンシャル

集団変数 ξ に対する自由エネルギープロファイル W(ξ) は平衡確率密度 p(ξ) により W(ξ)=−kBT ln p(ξ)+定数 と定義されます。ξ に明示されない他の座標は積分消去されています。

ΔG(A→B)=−kBTln(PB/PA)ΔG(A→B) = −k_B T ln(P_B/P_A)

二つの盆地の平衡存在確率が自由エネルギー差を決め、各極小点のエネルギーだけでは決まりません。エントロピー寄与は各盆地が占める配置空間の広さを反映します。

方法戦略注意点
アンブレラサンプリング重なり合う窓にバイアス窓の重なりと再重み付けが重要
メタダイナミクス集団変数に沿ってバイアスを蓄積不適切な変数は遅いモードを隠す
自由エネルギー摂動状態間で標本を再重み付け位相空間の重なりが必要

上級障壁・座標・不確かさ

反応座標は縮約記述です。不適切な選択は直交する遅い再配列を隠し、ヒステリシスを生むことがあります。自由エネルギー法では窓間の重なり、可逆性または収束と不確かさ評価を示す必要があり、滑らかな曲線だけでは証拠として不十分です。

研究最前線：希少事象と信頼できる自由エネルギー

参考文献