| 研究概要 |
Mg_2SiO_4オリビン、変型スピネル、及びスピネルは地球上部マントル、マントル遷移層の最重要構成鉱物である。平成10年度は、それら3種の多形の高温高圧物性計算に適用できる高精度な有効原子間ポテンシャルを求めることを主目標に置いた。
結晶のポテンシャルエネルギーを、クーロン項、ファンデァワールス引力項、反発項から成る二体間相互作用の和で表した。加えて、結晶中における多体相互作用を取り扱うべく、酸素イオンについてはbreathingシェルモデルを考慮した。必要なエネルギーパラメータは、3多形の実測の結晶構造、熱膨張、圧縮率、弾性定数を可能な限り再現するとの条件を用いて経験的に求めた。
3多形のそれぞれについて、求めた物性が非常に多様で、且つ取り扱った温度圧力範囲が極めて広範であるにもかかわらず、実測値と分子動力学シミュレーションによる計算値の一致は極めて満足ゆくものであった。例えば、300Kにおける静水圧縮については、実測データが存在する常圧から17GPa(オリビン)、27GPa(変型スピネル)、31GPa(スピネル)に至る圧力範囲で、3多形のそれぞれについて、実測のモル体積の圧力依存を極めて高精度で再現することに成功した。常圧下における体積熱膨張についても300Kから2100K(オリビン)、1100K(変型スピネル)、1000K(スピネル)に至る温度範囲でそれぞれの実測値を高精度で再現した。300K,OGPaにおける弾性定数(独立な弾性定数は、オリビンと変型スピネルでは共に9個、スピネルでは3個)とそれらの温度圧力依存についても、3多形の各々について、計算値が高精度で実測値を再現することを見出した。
|
