化学や熱力学の学習を進めていくと、必ずと言っていいほど登場するのが活量係数という概念です。
モル濃度やモル分率だけで反応や平衡を説明しようとすると、実際の溶液の挙動とズレが生じてしまうことがあります。
そのズレを補正するために使われるのが、まさに活量係数というパラメータです。
とはいえ、活量係数の求め方や計算方法、そしてその本質的な意味については、なんとなく公式だけを覚えて終わってしまっている方も多いのではないでしょうか。
本記事では活量係数の求め方は？計算方法と意味も！というテーマで、化学・溶液・熱力学・非理想溶液・γといったキーワードを軸に、基礎から実践的な計算例までを丁寧に解説していきます。
デバイ-ヒュッケルの式や非理想溶液のモデルなど、実務や試験で役立つ内容も盛り込んでいますので、ぜひ最後までご覧ください。
活量係数の求め方の結論とは
それではまず活量係数の求め方の結論について解説していきます。
結論から言うと、活量係数γ（ガンマ）は溶液の種類や濃度域によって適した計算モデルを選び、実測値または理論式から逆算する形で求めます。
理想溶液からのズレを表す係数であるため、理想的な挙動をどれだけ基準として設定するかによって、求め方そのものが変わってくるのです。
希薄な電解質溶液であればデバイ-ヒュッケルの式が使われますし、非電解質の非理想溶液であればマーギュラス式やファンラール式、あるいはNRTLやUNIQUACといった活量係数モデルが用いられます。
活量係数は実測と理論の両輪で求める
活量係数を求める方法は、大きく分けて実験的アプローチと理論的アプローチの二つに分類できます。
実験的アプローチでは、蒸気圧測定や凝固点降下、浸透圧測定などのデータから逆算して活量係数を導き出します。
理論的アプローチでは、イオン強度や分子間相互作用パラメータをもとに、数式モデルを用いて予測的に算出します。
実務では両者を組み合わせ、理論式で概算した値を実験データで補正するという流れが一般的でしょう。
濃度域によって使うべき式が異なる
活量係数の求め方は、対象となる溶液の濃度域によって大きく変わってきます。
極めて希薄な溶液であれば理想溶液に近いためγはほぼ1に近づきますが、濃度が上がるにつれてイオン同士や分子同士の相互作用が無視できなくなります。
そのため低濃度域ではデバイ-ヒュッケルの極限則、中濃度域では拡張デバイ-ヒュッケル式やデービス式、高濃度域ではピッツァー式など、段階的に精密なモデルへ切り替える必要があるのです。
結論部分のまとめ
活量係数の求め方の結論は、溶液の濃度域と種類に応じて適切なモデル式を選択し、実測データと理論式をすり合わせながら算出するという点に尽きます。
単一の万能公式が存在するわけではなく、状況に応じた使い分けこそが活量係数を正しく求めるための最大のポイントと言えるでしょう。
活量係数とはそもそも何を意味するのか
続いては活量係数とはそもそも何を意味するのかを確認していきます。
活量係数とは、実在する溶液の挙動を理想溶液の挙動からの補正係数として数値化したものです。
化学ポテンシャルや化学平衡の計算では、本来はモル濃度やモル分率をそのまま使いたいところですが、実在溶液では分子間力やイオン間の相互作用によって理想的な挙動からズレが生じます。
このズレを補正するために導入される概念が活量（activity）であり、活量係数γは活量と濃度を結びつける橋渡し役を担っています。
活量と濃度の関係式
活量aと濃度cの関係は、活量係数γを用いて次のように表されます。
a = γ × c
ここでaは活量、γは活量係数、cはモル濃度やモル分率などの濃度指標です。
この式からも分かるように、γが1であれば活量と濃度は完全に一致し、溶液は理想溶液として振る舞っていることになります。
逆にγが1から離れるほど、実在溶液の非理想性が強く現れていると解釈できるでしょう。
理想溶液とのズレを表す指標
理想溶液とは、分子同士やイオン同士の相互作用がすべて均一であると仮定した仮想的な溶液のことです。
現実の溶液では、溶質同士の静電的な引力や斥力、水和構造の違いなどによって、理想からのズレが必ず発生します。
活量係数はこのズレの大きさと方向を数値として表現するものであり、γが1より小さければ引力的な相互作用が支配的、γが1より大きければ斥力的な相互作用が支配的と考えられます。
なぜ活量係数が重要視されるのか
化学平衡定数や電極電位、溶解度積などを正確に扱うためには、モル濃度ではなく活量を用いる必要があります。
特に電解質溶液や高濃度の非理想溶液を扱う分析化学、電気化学、材料工学の分野では、活量係数を無視すると計算結果に大きな誤差が生じてしまいます。
そのため活量係数は、単なる補正値ではなく、実在系を正しく理解するための本質的なパラメータとして位置づけられているのです。
活量係数の計算方法にはどのようなものがあるか
続いては活量係数の計算方法にはどのようなものがあるかを確認していきます。
活量係数の計算方法は、対象が電解質溶液か非電解質溶液かによって大きく異なります。
電解質溶液ではイオン強度をベースにしたモデルが中心となり、非電解質の非理想溶液では分子間相互作用パラメータをベースにしたモデルが用いられます。
デバイ-ヒュッケルの式による計算
電解質溶液における活量係数の代表的な計算方法が、デバイ-ヒュッケルの極限則です。
log γ± = −A × z+ × z− × √I
ここでγ±は平均活量係数、Aは定数（水溶液・25℃で約0.509）、z+とz−はイオンの価数、Iはイオン強度を表します。
この式は希薄溶液、目安としてイオン強度が0.01mol/kg程度以下の領域で高い精度を発揮します。
濃度が上がるにつれて誤差が大きくなるため、より濃い溶液では拡張デバイ-ヒュッケル式やデービス式への切り替えが必要になるでしょう。
非理想溶液モデルによる計算
非電解質の非理想溶液では、マーギュラス式やファンラール式、さらにはNRTLモデルやUNIQUACモデルといった過剰ギブスエネルギーに基づく計算方法が使われます。
これらのモデルは、分子間相互作用パラメータを実験データからフィッティングして決定し、そこから活量係数を算出するという流れです。
プロセス設計や蒸留計算のシミュレーションでは、こうしたモデルがソフトウェア上で標準的に組み込まれています。
計算モデルの比較表
ここで代表的な活量係数の計算モデルを、適用範囲とともに整理してみましょう。
| モデル名 | 対象溶液 | 適用濃度域の目安 |
|---|---|---|
| デバイ-ヒュッケル極限則 | 電解質溶液 | イオン強度0.01mol/kg以下 |
| 拡張デバイ-ヒュッケル式 | 電解質溶液 | イオン強度0.1mol/kg程度まで |
| デービス式 | 電解質溶液 | イオン強度0.5mol/kg程度まで |
| ピッツァー式 | 電解質溶液 | 高濃度域まで対応可能 |
| マーギュラス式・ファンラール式 | 非電解質の非理想溶液 | 中程度の非理想性 |
| NRTL・UNIQUAC | 非電解質の非理想溶液 | 強い非理想性・多成分系 |
このように、対象となる溶液の性質と濃度域を見極めた上で、適切な計算モデルを選ぶことが求められます。
非理想溶液における活量係数の求め方の実践
続いては非理想溶液における活量係数の求め方の実践について確認していきます。
非理想溶液とは、ラウールの法則から外れる挙動を示す溶液の総称であり、多くの実在溶液がこれに該当します。
非理想溶液における活量係数の求め方には、実験的に蒸気圧を測定するアプローチと、状態方程式や過剰ギブスエネルギーモデルを用いるアプローチが存在します。
蒸気圧測定からの逆算
非理想溶液の活量係数は、気液平衡における蒸気圧データから求めることができます。
γi = Pi / (xi × Pi*)
ここでPiは成分iの実測分圧、xiは液相のモル分率、Pi*は純成分iの蒸気圧です。
実測された分圧を理想溶液での理論値と比較することで、その差からγを直接算出できるのです。
この方法は精度が高い一方で、正確な気液平衡装置と分析機器が必要になる点には注意が必要でしょう。
過剰ギブスエネルギーからの算出
過剰ギブスエネルギー(GE)は、実在溶液が理想溶液からどれだけエネルギー的にズレているかを表す量です。
このGEを各成分のモル分率で偏微分することで、それぞれの成分の活量係数を求めることができます。
RT ln γi = (∂GE/∂ni)T,P,nj
この関係式はギブス-デュエムの式とも呼ばれる基礎関係の一部であり、多くの活量係数モデルの理論的な土台となっています。
マーギュラス式やNRTLモデルは、このGEの数学的表現の違いによって使い分けられていると理解すると分かりやすいでしょう。
共沸点や相分離との関係
非理想溶液の活量係数を求める際には、共沸点や相分離といった現象との関連にも注目する必要があります。
活量係数が大きく1から外れる系では、正の偏差を示す場合に最低共沸点、負の偏差を示す場合に最高共沸点が現れることがあります。
さらに非理想性が極端に強い場合には、二液相に分離する液液平衡が生じることもあるため、蒸留プロセスの設計では活量係数の挙動を事前に把握しておくことが欠かせません。
活量係数の求め方における具体的な計算例
続いては活量係数の求め方における具体的な計算例を確認していきます。
理論だけでは実感がわきにくいため、ここでは電解質溶液を例にした計算の流れを追ってみましょう。
イオン強度を求める
まずは対象となる電解質溶液のイオン強度Iを求める必要があります。
I = 1/2 × Σ(ci × zi^2)
例えば0.01mol/kgの塩化ナトリウム水溶液の場合、Na+とCl−はどちらも価数1なので、I = 1/2 × (0.01×1^2 + 0.01×1^2) = 0.01mol/kgとなります。
この値が、次のデバイ-ヒュッケル式に代入するための基礎データとなるのです。
デバイ-ヒュッケル式に代入する
先ほど求めたイオン強度をデバイ-ヒュッケルの極限則に代入してみましょう。
log γ± = −0.509 × 1 × 1 × √0.01 = −0.0509
これを計算するとγ± ≈ 0.889という値が得られます。
つまりこの塩化ナトリウム水溶液では、理想溶液と比較して活量がおよそ11パーセント低く見積もられる状態であると解釈できます。
濃度が非常に薄いにもかかわらず、すでに理想溶液からのズレが確認できる点は興味深いのではないでしょうか。
濃度が上がった場合の補正計算
次に、同じ塩化ナトリウムでも濃度が0.1mol/kgまで上がった場合を考えてみましょう。
この濃度域ではデバイ-ヒュッケルの極限則の精度が落ちるため、拡張デバイ-ヒュッケル式を使うのが一般的です。
log γ± = −A × |z+ × z−| × √I / (1 + B × a × √I)
ここでBとaはイオンの大きさなどに関する経験的パラメータであり、実測データに基づいて設定されます。
このように濃度が上がるほど計算式に補正項が追加され、より複雑なモデルへと移行していくことが分かります。
活量係数を求める際に注意すべきポイント
続いては活量係数を求める際に注意すべきポイントを確認していきます。
活量係数の計算では、公式に数値を当てはめるだけでなく、いくつかの前提条件を正しく理解しておくことが欠かせません。
基準状態の取り方に注意する
活量係数は、どのような基準状態を採用するかによって数値そのものが変わってきます。
代表的な基準にはラウール基準とヘンリー基準があり、溶媒側にはラウール基準、希薄な溶質側にはヘンリー基準が使われることが一般的です。
基準状態を取り違えると、同じ溶液を扱っていてもまったく異なるγの値が算出されてしまうため、モデルや文献を参照する際には必ず基準を確認しましょう。
適用できる濃度域を超えないようにする
デバイ-ヒュッケルの極限則は非常に有名な式ですが、万能ではありません。
適用範囲を超えた高濃度域で使用してしまうと、実測値との誤差が急激に拡大してしまいます。
濃度域に応じたモデルの使い分けこそが、活量係数計算の精度を左右する最大の分かれ道と言っても過言ではないでしょう。
温度依存性を見落とさない
活量係数は濃度だけでなく温度にも依存するパラメータです。
デバイ-ヒュッケル式の定数Aも温度によって変化するため、25℃以外の条件で計算する場合は温度補正済みの定数を用いる必要があります。
プロセス設計や研究データを扱う際には、測定温度と計算モデルの前提温度が一致しているかどうかを必ず確認するようにしましょう。
活量係数を求める際は、基準状態・適用濃度域・温度条件という三つの前提をそろえることが、正確な計算結果を得るための必須条件です。
どれか一つでもズレていると、モデル自体は正しくても得られる数値の意味が変わってしまう点には十分注意しましょう。
まとめ
今回は活量係数の求め方は？計算方法と意味も！というテーマで、化学・溶液・熱力学・非理想溶液・γに関する内容を解説してきました。
活量係数とは、実在する溶液が理想溶液からどれだけズレているかを表す補正係数であり、溶液の種類や濃度域に応じて求め方が変わってきます。
電解質溶液ではデバイ-ヒュッケルの式や拡張式、非電解質の非理想溶液ではマーギュラス式やNRTL、UNIQUACといったモデルが用いられ、いずれも実測データと理論式のすり合わせによって精度が高められています。
基準状態や適用濃度域、温度条件といった前提をそろえることを忘れなければ、活量係数の計算はそれほど難しいものではないでしょう。
ぜひ本記事の内容を参考に、ご自身が扱う溶液系に最適な活量係数の求め方を選択してみてください。
