本研究では、実験的に困難である占有率の圧力依存性の理諭的予測法を確立し、それを水素クラスレートに適用した。Grandcanonical Monte Carlo(GCMC)シミュレーションは、グランドカノニカル集団(温度T、体積V、化学ポテンシャルμが決定された、粒子の出入りがある開いた系)を扱ったシミュレーション法である。仮想的な、ゲスト分子の気体とクラスレート内のケージを占有するゲスト分子との平衡を考えると、ゲスト分子のケージ占有はケージを吸着点とみなすことで吸着現象として扱うことができる。このとき占有しているゲスト分子は、ゲスト流体と平衡にある学ポテンシャルμを持つグランドカノニカル集団としてみなす。一方、ホストの結晶格子をなす水分子は統計力学のモデル計算ではカノニカル集団(温度T、体積V、粒子数Nが決定された閉じた系)として扱うのが便利であるが、シミュレーションではホストである水に対しては(N、T、圧力P)一定のほうが現実的である、またゲストに対しては(μPT)一定とし、また水素の化学ポテンシャルμはPとTで規定される値に定める。このために、余分なガスなどを考慮することなく、クラスレートの体積の温度と圧力依存性が得られることになる。事前に決定した化学ポテンシャルを用いて、水素クラスレートの占有状態を明らかにするためのGCMCシミュレーションを行った。大ケージでは圧力に応じて水素分子の占有状態が大きく変化する。、高圧になるにつれてより多くの水素分子がケージに入っていき、600MPaでは大ケージの40%近くを最大4個の水素分子が占有するという結果が得られた。これに対して、小ケージは圧力をあげても占有状態はあまり変化せず、ほとんどのケージが単一の水素分子に占められるという結果となった。
|