分子の反応性や安定性を調べるためにはギブズエネルギーなどの自由エネルギーの情報は欠かせません。
じつは量子化学計算パッケージを使うことで、分子軌道だけではなく、自由エネルギーを簡便に調べることができます。
このページではフリーソフトであるGAMESSの計算結果を使って、分子の自由エネルギーがどのようにして得られるかを紹介します。
- 自由エネルギーへは電子エネルギーの寄与が最も大きい
- 自由エネルギーへの補正値は統計力学の知見を使うと得られる
- 並進運動は理想気体を仮定する
- 回転運動は対称性の考慮が重要
- 振動運動は全ての振動モードを取り扱う
自由エネルギーの寄与
このページは、量子化学計算パッケージのGAMESSを使った水分子の計算結果をもとに作成しています。計算ログは下からダウンロードできます。(量子化学計算条件の詳細はログを確認願います。)
統計力学を使えば自由エネルギーがわかる
まず重要な知識ですが、自由エネルギーには電子エネルギーがほとんどの寄与を占めます。
例えば基底状態の水分子に対して量子化学計算を行うと、電子状態のエネルギー値として
-----------------
ENERGY COMPONENTS
-----------------
(中略)
TOTAL ENERGY = -76.0257177698

電子状態のエネルギーが大きいため、化学反応熱が大きくなります。
自由エネルギーはこの電子エネルギーに補正を加えることで得ることができます。
ギブズエネルギー

GAMESSのデフォルト設定では、温度

つまり、内部エネルギー
では、どのようにして
答えは統計力学にあります。
統計力学に基づいた分配関数
Q LN Q
ELEC. 1.00000E+00 0.000000
TRANS. 3.00431E+06 14.915558
ROT. 4.02909E+01 3.696126
VIB. 1.00029E+00 0.000289
TOT. 1.21081E+08 18.611973
これはGAMESSによる水分子の計算結果です。Qが分配関数

通常、電子エネルギーでは熱励起はほとんどないので
自由エネルギーは並進・回転・振動に分けて考える
では、どのようにして分子の分配関数
「複雑なものは分けて考えよう」ということで、並進運動・回転運動・振動運動のそれぞれの寄与を合算することで、分配関数や自由エネルギーを計算することができます。

例えば
統計力学を利用すると分子の並進、回転、振動の分配関数やエントロピーは次の式を利用することで計算できます。

GAMESSでも各寄与ごとに足しあわされて内部エネルギーやギブズエネルギーなどが計算されています。(下は水分子の例)
E H G CV CP S
KJ/MOL KJ/MOL KJ/MOL J/MOL-K J/MOL-K J/MOL-K
ELEC. 0.000 0.000 0.000 0.000 0.000 0.000
TRANS. 3.718 6.197 -36.975 12.472 20.786 144.800
ROT. 3.718 3.718 -9.162 12.472 12.472 43.203
VIB. 59.288 59.288 59.281 0.160 0.160 0.022
TOTAL 66.725 69.204 13.144 25.103 33.417 188.025
例えば、

以下では、熱力学量に対する並進、回転、振動による寄与についてそれぞれ解説していきます。
並進運動による寄与
背景
並進運動はエネルギーとしては大きな値は持ちませんが、エントロピーには大きな寄与があります。これは並進運動はほとんど連続したエネルギー準位を持つことに寄ります。
質量


水でもベンゼンでも
内部エネルギーについては、理想気体として取り扱うので、
エントロピーは分配関数を利用することで、
実際の計算例
GAMESSによる水分子の計算結果を見てみましょう。
E H G CV CP S
KJ/MOL KJ/MOL KJ/MOL J/MOL-K J/MOL-K J/MOL-K
ELEC. 0.000 0.000 0.000 0.000 0.000 0.000
TRANS. 3.718 6.197 -36.975 12.472 20.786 144.800
ROT. 3.718 3.718 -9.162 12.472 12.472 43.203
VIB. 59.288 59.288 59.281 0.160 0.160 0.022
TOTAL 66.725 69.204 13.144 25.103 33.417 188.025
まず、並進は「TRANS.」で表されており、
また、エントロピーについても、


単位換算がややこしければSI単位系で常に考えましょう
この得られた
回転運動による寄与
背景
3つの慣性モーメント
まず、内部エネルギーについては、並進運動と同様にエネルギー等分配の法則から、
一方、エントロピーについては分配関数として

例えば水分子だと
慣性モーメント
これを利用すると、エントロピーは
実際の計算例
GAMESSによる水分子の計算結果を見てみましょう。
E H G CV CP S
KJ/MOL KJ/MOL KJ/MOL J/MOL-K J/MOL-K J/MOL-K
ELEC. 0.000 0.000 0.000 0.000 0.000 0.000
TRANS. 3.718 6.197 -36.975 12.472 20.786 144.800
ROT. 3.718 3.718 -9.162 12.472 12.472 43.203
VIB. 59.288 59.288 59.281 0.160 0.160 0.022
TOTAL 66.725 69.204 13.144 25.103 33.417 188.025
回転は「ROT.」で表されており、
次にエントロピーについてですが、
THE ROTATIONAL CONSTANTS ARE (IN GHZ)
988.63138 406.32750 287.97130
つまり、
振動運動による寄与
背景
内部エネルギーに対する振動モードによる寄与は並進・回転に比べるとかなり大きくなります。それはゼロ点振動によって書く振動モードが必ず一定以上のエネルギーを持つためです。
振動運動のエネルギー準位幅は並進・回転に比べて大きいため、内部エネルギーについてエネルギー等分配の法則を利用して
ここで、振動数
エントロピーについては、
実際の計算例
GAMESSによる水分子の計算結果を見てみましょう。
E H G CV CP S
KJ/MOL KJ/MOL KJ/MOL J/MOL-K J/MOL-K J/MOL-K
ELEC. 0.000 0.000 0.000 0.000 0.000 0.000
TRANS. 3.718 6.197 -36.975 12.472 20.786 144.800
ROT. 3.718 3.718 -9.162 12.472 12.472 43.203
VIB. 59.288 59.288 59.281 0.160 0.160 0.022
TOTAL 66.725 69.204 13.144 25.103 33.417 188.025
振動は「VIB.」で表されており、


水分子だと
水の振動モードはGAMESSによってすべて出力されており、次の通りです。
MODE FREQ(CM**-1) SYMMETRY RED. MASS IR INTENS.
(中略)
7 1688.567 A1 1.086338 2.259042
8 4033.275 A1 1.041776 0.471534
9 4189.330 B1 1.087619 2.144479
振動モードの波数
最後に振動のエントロピーについても
まとめ
このページではフリーソフトであるGAMESSの計算結果を参考にしながら、分子の自由エネルギーがどのようにして得られるかを紹介しました。
分子の反応性や安定性を調べるためにはギブズエネルギーなどの自由エネルギーの情報は欠かせません。
この自由エネルギーは量子化学計算パッケージを使うことで簡便に調べることができます。また、これらの自由エネルギーの計算結果が統計力学の知見を使っていることを知っていただければと思います。