[論文レビュー] Inner Shell Definition and Absolute Hydration Free Energy of K+(aq) on the Basis of Quasi-chemical Theory and Ab Initio Molecular Dynamics

Susan B. Rempe, D. Asthagiri+1 more|ArXiv.org|Oct 30, 2003

Spectroscopy and Quantum Chemical Studies被引用数 7

ひとこと要約

本研究では、準化学理論とab initio分子動力学（AIMD）を組み合わせて、K+(aq)の内層水和殻を定義し、その絶対的水和ギブス自由エネルギーを計算した。4つの水分子が化学的に定義された内層と判明し、径分布関数に明確な肩が観察されたことから、より弱く結合した第二の水和集団が存在することが示され、水和ギブス自由エネルギーは−70.5 ± 2.1 kcal/molと算出された。これは、類縁種の実験的推定値と良好に一致している。

ABSTRACT

K+(aq) ion is an integral component of many cellular processes, amongst which the most important, perhaps, is its role in transmitting electrical impulses along the nerve. Understanding its hydration structure and thermodynamics is crucial in dissecting its role in such processes. Her we address these questions using both the statistical mechanical quasi-chemical theory of solutions and ab initio molecular dynamics simulations. Simulations predict an interesting hydration structure for K+(aq): the population of about six (6) water molecules within the initial minimum of the observed g_{KO}(r) at infinite dilution involves four (4) inner-most molecules that the quasi-chemical theory suggests should be taken as the theoretical inner shell. The contribution of the 5th and 6th closest water molecules is observable as a distinct shoulder on the principal maximum of the g_{KO}(r). The quasi-chemical estimate of solvation free energy for the neutral pair KOH is also in good agreement with experiments.

研究の動機と目的

• K+(aq)の内層水和殻を理論的に一貫した定義を確立すること。これは、径分布関数に明確な第一最小がないためである。
• 準化学理論とab initio分子動力学を組み合わせたハイブリッド手法を用いて、K+(aq)の絶対的水和ギブス自由エネルギーを計算すること。
• 中性種（例：KOH）の実験データと照合することで、先行研究における水和ギブス自由エネルギー推定値の不一致を解消すること。
• Li+ や Na+ などの小さなアルカリ金属イオンとは異なる、K+ の水和の構造的・熱力学的特徴を明確にすること。
• 拡散的な溶媒層構造を示すイオンに対して、理論的およびシミュレーション手法を統合した堅牢な内層定義フレームワークを提供すること。

提案手法

• K+ および水を用いた無限希釈条件下で、VASPコードを用いてab initio分子動力学（AIMD）シミュレーションを実施。K+ と水のためのPW91汎関数と超ソフト擬ポテンシャルを用いた。
• 32個の水分子を含むシミュレーションボックスを用い、周期的境界条件を適用。300 Kで体積密度を液体水と同等に維持し、20 ps間の平衡化を実施した。
• 準化学理論を用いて溶媒層を内層と外層に分割。内層は量子力学的に取り扱い、外層は連続体モデルまたは力場モデルを用いた。
• K+ と水分子酸素原子間の径分布関数 g(r) を分析し、構造的サブグループを特定し、第一協調殻の位置を同定した。
• 準化学理論による内層寄与と古典的MDによる外層寄与を組み合わせることで、水和ギブス自由エネルギーを推定した。
• 中性種KOHの計算された水和ギブス自由エネルギーを実験データと比較することで、結果の妥当性を検証し、標準状態の表記慣習に合わせて補正した。

実験結果

リサーチクエスチョン

• RQ1 g_{KO}(r) に明確な第一最小がないにもかかわらず、K+(aq) の内層水和殻の正しい理論的定義は何か？
• RQ2 K+(aq) の水和構造は、Li+(aq) や Na+(aq) とはどのように異なっているか。特に、径分布関数の特徴に注目して。
• RQ3 準化学理論とab initio分子動力学を組み合わせることで、K+(aq) の絶対的水和ギブス自由エネルギーを正確に推定できるか？
• RQ4 径分布関数に観察される特徴的要素（主ピークの肩）が、明確に区別できる水和集団に対応している程度はどの程度か？
• RQ5 KOH種の計算された水和ギブス自由エネルギーは、実験値とどの程度一致するか。これにより、全体の手法的アプローチの妥当性がどの程度裏付けられるか？

主な発見

• 準化学理論により、K+(aq) の内層水和殻として最も確率の高い水分子数は4つと特定された。これは、最も強く結合したサブグループに基づく。
• ab initio分子動力学シミュレーションにより、g_{KO}(r) の第一ピークが広がっており、r ≈ 3.0 Åに明確な肩が観察された。これは、内層の4つを除いた2つの水分子からなる第二の水和集団が存在することを示唆している。
• K+ から3.50 Å以内の水分子の平均数は5.9 ± 0.9であるが、そのうち化学的に定義された内層は、協調エネルギーおよび動的挙動の差異から4分子が最も適切に表している。
• K+(aq) の水和ギブス自由エネルギーは−70.5 ± 2.1 kcal/molと算出された。これは、標準状態の表記慣習に合わせた調整後、実験的推定値と整合的である。
• 中性KOH種の水和ギブス自由エネルギーは−194.3 ± 2.1 kcal/molと算出され、実験値の−191 kcal/molと良好に一致した。
• 結果から、Li+ や Na+ とは異なり、K+(aq) の内層は、第一最小が浅く、径分布関数における重なり合う集団のため、はっきりと定義されていないことが明らかになった。