Tanaka, H., Matsumoto, M. & Yagasaki, T. Cage occupancies of CH4, CO2, and Xe hydrates: Mean field theory and grandcanonical Monte Carlo simulations. J. Chem. Phys. 160 (2024) doi:10.1063/5.0188679
「クラシックな氷」とは異なり、水分子がまるで「かご（ケージ）」のような立体枠組みを作り、その内部にメタンや二酸化炭素、キセノンといった気体分子を取り込んだ物質を「クラスレート水和物（ガスハイドレート）」と呼びます。将来のエネルギー資源や二酸化炭素の貯蔵担体として期待されるこの物質ですが、各ケージにどれくらいの割合で気体分子が閉じ込められるか（占有率）を正確に予測することは、その熱力学的な安定性を知る上で極めて重要です。
これまでの理論では、ケージの中に入った気体分子（ゲスト分子）同士の相互作用は単純化されがちでした。本研究では、隣り合うケージに入ったゲスト分子同士が及ぼし合う力を考慮できる新しい統計力学理論を提案しました(TMY2022, TMY2023)。各分子が周囲から受ける影響を平均的な場として扱う「平均場近似」を取り入れることで、ゲスト分子同士の相互作用がケージ占有率に与える影響を正確に見積もることに成功しました。
理論の正しさを確かめるため、分子の出し入れを計算できる「大正準モンテカルロ（GCMC）シミュレーション」を行いました。さらに研究グループは、ケージの中心のみにゲスト分子が吸着すると仮定して計算を大幅に効率化する「格子GCMC」という新手法も開発しました。これらの計算結果は提案した理論と見事に一致し、温度が低い場合や分子間の引き合う力が強い場合には、ケージ内の気体密度が大きくゆらぐ臨界現象のような挙動が起きることも明らかになりました。
Abstract
本研究では、クラスレート水和物の熱力学的安定性に関する統計力学理論を提案し、ケージの占有率に対するゲスト分子間相互作用の影響を考察した。この影響の大きさを調べるため、平均場近似を展開した。我々の新手法は非常に良好に機能し、2種類のグランドカノニカル・モンテカルロ（GCMC）シミュレーションによってその有効性が示された。1つは通常のGCMCであり、もう1つは本研究でクラスレート水和物のために設計された「格子GCMC」である。後者のシミュレーションでは、各ゲスト分子はケージの中心のいずれかに吸着されるとし、他のゲストが存在しない状態でのホスト-ゲスト相互作用をケージ占有の自由エネルギーに組み込むことで、ゲスト分子間の相互作用のみを明示的に扱う。温度が低い場合やゲスト分子間の相互作用が強い場合には、大きな密度ゆらぎのようなゲスト種の臨界現象が観察された。(Geminiによる概要の機械翻訳)
We propose a statistical mechanical theory for the thermodynamic stability of clathrate hydrates, considering the influence of the guest-guest interaction on the occupancies of the cages. A mean field approximation is developed to examine the magnitude of the influence. Our new method works remarkably well, which is manifested by two sorts of grandcanonical Monte Carlo (GCMC) simulations. One is full GCMC, and the other is designed in the present study for clathrate hydrates, called lattice-GCMC, in which each guest can be adsorbed at one of the centers of the cage. In the latter simulation, only the guest-guest interaction is explicitly treated, incorporating the host-guest interaction into the free energy of the cage occupation without other guests. Critical phenomena for guest species, such as large density fluctuations, are observed when the temperature is low or the guest-guest interaction is strong.
research papers paper2024 clathratehydrate simulation GCMC thermodynamics