分子軌道法および分子動力学法の複合的計算による氷結晶における粒界拡散の研究

利用分子轨道和分子动力学方法联合计算研究冰晶中的晶界扩散

• 批准号: 15760001
• 负责人: 堀 彰
• 金额: $ 2.3万
• 依托单位: Hokkaido University
• 依托单位国家: 日本
• 项目类别: Grant-in-Aid for Young Scientists (B)
• 财政年份: 2003
• 资助国家: 日本
• 起止时间: 2003 至 2004
• 项目状态: 已结题
• 来源: https://kaken.nii.ac.jp/grant/KAKENHI-PROJECT-15760001/
• 关键词: 氷; 拡散; 粒界; 分子軌道法; 分子動力学法

氷中の物質の粒界拡散に関しては、氷床コアの中の物質移動を考える上で重要であるが、実験的にも理論的にもほぼ手つかずの状態である。本研究では氷の粒界拡散の理論的な取り扱いのため、分子軌道法と分子動力学法の複合的な計算を行うものである。実際には(1)構造モデルの作成と(2)計算に取り組んだ。(1)構造モデルの作成に関して、氷の構造にはアイス・ルールと呼ばれる水分子の配列の規則があるため、周期的境界条件を満たすモデルの作成は困難である。そこで、アイス・ルールおよび周期的境界条件を満たす、任意の大きさの構造モデルの作成プログラムの開発に取り組んだ。氷結晶中の個々の水分子は6通りの配置を取り得るので、N個の系では6^N通りの配置から上記の条件を満たす配置を探索するプログラムを作成したが、ごく小さな系を除いて、意図した構造を得ることができなかった。また、初期構造の周期的境界条件下での構造最適化を非経験的分子軌道法プログラムのGaussian03を用いて試みたが、プログラム上の問題で収束しなかったため、結晶の構造データをそのまま利用した。(2)(001)面どうしを組み合わせて擬似的に粒界層を作成し、その間にメタン分子を配置してその拡散の障壁エネルギーを求めた。周期的な境界条件では問題が生じたため課さなかった。水分子24個の層2層(計48個)を組み合わせた96個の水分子からなる界面構造で、メタンとメタンスルホン酸(MSA)の生成熱を半経験的分子軌道法プログラムのMOPACで計算したところ、後者の方が300kJ/mol以上も低く、MSAが大きな拡散係数を持つ原因と関連していると考えられる。分子動力学法による計算に関しては、周期的境界条件を満たす構造モデルの作成がその途上であるので、拡散係数の前指数因子を求めるための「熱力学的積分法」による自由エネルギー計算を実行する段階まで到達できなかった。

The material boundary of the material is dispersed, and the material movement in the bed is examined. It is important to understand the theoretical state of the material. In this study, the molecular orbital method and the molecular dynamics method are used to calculate the particle dispersion theory. In practice,(1) the structure of the structure is made (2) the calculation is made. (1)The structure of the structure is related to the structure of the structure The boundary conditions of the system, the structure of the system and the development of the system are determined. The water molecules in the crystal are arranged in 6 ways, and the N systems are arranged in 6 ways, and the upper conditions are explored. The structural optimization under the condition of initial structural period is carried out by non-theoretical molecular orbital method. The Gaussian03 method is used to solve the problem of structural optimization under the condition of initial structural period. (2) The (001) plane is composed of two layers, one layer is composed of two molecules, the other layer is composed of two molecules, and the other layer is composed of two molecules. The periodic boundary condition is not a problem. Water molecules 24 layers 2 layers (48 in total) are grouped into 96 water molecules. The interface structure, heat of formation of MSA, and the molecular orbital method are used to calculate the dispersion coefficient of MSA. The latter is more than 300kJ/mol. The calculation of molecular dynamics method is related to the boundary condition of the period. The structure of the calculation method is related to the calculation of the dispersion coefficient and the former exponential factor. The calculation method of thermodynamic integration is related to the calculation of the stage.