| 研究概要 |
計算機の進展に伴い、αb initio動力学が近年行われるようになってきているが、計算量を減らすため、高精度の計算手法は利用できないのが現状であることを考えると、ポテンシャル面を簡便にαb initio計算の結果から構成する方法が望まれる。CollinsらはShepard内挿法とトラジェクトリ計算を組み合わせたポテンシャル面の逐次改善法を用い、ポテンシャル面を決定する応用を行ってきた。しかしながら、この手法では、ポテンシャル面の二次微分を用いなければならないため、最新の量子化学的手法と組み合わせて用いることが困難である。そこでわれわれは、エネルギーに関する座標微分を全く使わない方法として、内挿移動最小二乗法(interpolant moving least squares,IMLS)法とShepard内挿法を併用した方法(IMLS/Shepard法)を開発し、O+H_2系、OH+H_2系に応用した。O+H_2系では、IMLS/Shepard法、Shepard法のいずれも反応断面積をサンプル点が増えるにつれ、再現していく傾向があったが、サンプル点が同数のときにrms誤差を比べると、IMLS/Shepard法はShepard法より小さかった。とくに、あるエネルギー以下の時にサンプルを行う方法の誤差が小さい。Shepard法では、ポテンシャル面の2次微分までを要するのに対し、IMLS/Shepard法では、ポテンシャル面の微分を全く必要としないことを考えると驚くべきことである。OH+H_2系でも、IMLS/Shepard法の結果がShepard法の結果より誤差がはるかに小さく、実用的な誤差の値(3kcal/mol程度)を示した。さらなる精度の改善、自動化法の確立が課題であるが、最新の量子化学的手法と局所内挿法の結合による、ポテンシャル面の自動生成の可能性が4原子系以上の多体系に開けたと考えている。ab intio法については、来年度における融合を目指してEOM法に基づく方法を開発中である。
|
