|
シュレジンガー方程式の数値解法を研究した。
1) 2原子分子に限らず、直線分子の電子波動関数を数値積分法で求めるために、円筒座標を導入し、さらに、等間かく有限差分法の適用に合った座標変換を行なった。その結果、ハートリー・ホツク方程式は、2次元のシュレジンガー方程式となり、さらに、数4次元の固有値問題となる。本年度は、一電子問題の解法が完成し、【H(^+_2)】イオンについて、厳密解と比較することができた。【H(^(2+)_3)】イオンの数値解がこの方法によって始めて得られた。
2) シュレジンガー方程式の数値解法を、分子中の原子核の運動の方程式に応用した。この場合、有限要素法を用いた。2次元素の方程式に応用すると解くべき方程式は、大次元の対称一般固有値問題となる。本研究では数十方次元の一般固有値問題を解くアルゴリズムを堤出した。この研究ではスーパーコンピューターの機能を十分引きだすことに成功し、いくつかの応用計算を行なった。しかしながら、有限要素法において、一次の補間関数を使ったため、精度が十分でなかった。
3) 精度の向上を、有限要素法で企るには、きざみ(要素の大きさ)を小さくするのも一方法であるが、その結果は、解くべき一般固有値問題が巨大になる。他の方法は、補間関数の次数をあげることである。このことによる精度の向上の仕方を調べるため、一次元のシュレジンガー方程式の解法に応用した。4次の補間関数を用い、かつ、位置エネルギー項の部分の積分にも適当な補間をほどこすと、多くの例題において、10根以上が、実用上、厳密解と等しい値を高速でえられることが明らかになった。
|
