抄録
本研究では,アルミノケイ酸塩ガラスの分極力場を開発し,その分子動力学計算を用いて,ガラスの架橋構造について考察した.分極力場のパラメータは,密度汎関数法に基づく第一原理計算の分極状態と力を再現するように,Force Matching法を用いて最適化した.本研究では,まず等方的な斥力を仮定したモデルを用いて分極力場を最適化したところ,Alの四面体構造が歪んで,実験値と異なる5配位構造を形成した.一方で,イオンの変形を考慮した異方的な斥力モデルで最適化したところ,5配位構造は解消され,NMRスペクトルも実験値と良い一致を示した.そこで,ガラスの架橋構造を詳しく調べたところ,Al-O-Alの架橋酸素の双極子モーメントだけが顕著に小さいこと,またその架橋角度は比較的広く分布し,不安定であることを明らかにした.またこの架橋酸素は密度が増加すると大きく折れ曲がることから,アルミノケイ酸塩ガラスの架橋構造はAl-O-Alの架橋酸素が優先的に壊れて,三配位酸素となることを明らかにした.