造岩物質の構造と物性の理論

• 研究課題/領域番号: 60121013
• 研究種目: 特定研究
• 配分区分: 補助金
• 研究機関: 岡山大学
• 研究代表者: 松井 義人 岡山大学, 温研, 教授 (20033128)
• 研究期間 (年度): 1985
• 研究課題ステータス: 完了 (1985年度)
• 配分額: 4,300千円 (直接経費: 4,300千円); 1985年度: 4,300千円 (直接経費: 4,300千円)
• キーワード: ケイ酸塩 / ケイ酸塩メルト / 計算機シミュレーション / 高圧力 / 固体物性 / 鉱物化学 / 鉱物物理学

研究概要

本年度は、計画研究の初年度として、分子動力学計算、固体量子論(第一原理)計算についてこれまでに提示された成果の全面再評価を実行した。同時に、地球科学に無経験の物理系分担者と、問題の所在を討論し、将来最も有効な戦略は何かについて討議を進めた。これらの目的のため、地球内部研究センターにおいて8月25日〜27日にわた1)研究会を開くなど、情報・意見の交換に努力した。
分子動力学計算について得られた結果は、1)2体中心力モデルにおいては、形式電荷がそのままクーロン力を構成すると仮定した場合、圧力および剛性率などを再現することはほとんど不可能であること。したがって部分電荷を仮定するほかに道がないこと、2)従来の格子エネルギー最小化法によって推定された2体中心力のパラメタは、分子動力学法に適用できない(原子の変位が無限小のとき以外は全く無効)ものがほとんどであること。3)以上の結論をふまえた上で、本研究のために結晶構造が精密に決定されたMgSi【O_3】ペロフスカイトについて弾性率を強制条件とし、空間群の制約をとり払ったWMIN法を実行し、広い適用範囲を持つポテンシャルパラメタが得られたこと、4)同様に部分電荷を仮定したポテンシャルにより、Si【O_2】の多形転移の計算機実験に成功したこと、などである。
固体量子論計算から、実効ポテンシャルを導く試みも進行中であるが、今年度内に結論は出なかった。第一原理計算から分子動力学計算に利用できるポテンシャル関数を導き出せるか否かは、現在理論鉱物物理学の中心をなす問題であって、国際的競争が激烈であるので、特に重点的な補強が必要である。