研究課題/領域番号 |
12640349
|
研究機関 | 大阪大学 |
研究代表者 |
鈴木 直 大阪大学, 基礎工学研究科, 教授 (40029559)
|
研究分担者 |
小田 竜樹 金沢大学, 理学部, 講師 (30272941)
|
キーワード | 局所密度汎関数近似 / 分子動力学 / ウルトラソフト擬ポテンシャル法 / ノンコリニア磁性 / 第一原理格子振動計算 / 線形応答理論 / 液体酸素 / 固体酸素 |
研究概要 |
(1)ノンコリニアー磁性の第一原理分子動力学計算コードの開発 昨年までに開発してきたノンコリニアー磁性の第一原理分子動力学計算コードについて次のような改良・拡張を行った:(A)計算コードの並列処理高速化、(B)一般化密度勾配近似計算への拡張(C)結晶系への拡張、(D)金属系への拡張 (C)では、基本並進部分を除いた空間群の対称操作をあらかじめ計算コードに入力として与えることにより電荷密度の対称化を行い、ほとんどのブラベ格子に対応するセルを計算可能にした。また、(D)では、擬固有値行列を導入することによりフェルミ準位付近の波動関数の更新をスムースに行うことに成功した。これにより、アルミニウムやシリコンの融点直上の液体に適用して分子動力学を行うことも可能となった。 (2)線形応答理論に基づく格子振動計算コード開発 擬ポテンシャルと平面波基底を用いて、変分的密度汎関数摂動論に基づく第一原理格子振動計算プログラムの実装を行った。我々の方法の特徴は、まず第1に、原子変位に対する波動関数と電荷密度の応答を線形応答の範囲内で全エネルギーの二階微分を最小化することにより最適化したことと、第2に、一次の波動関数の最適化(エネルギーの二階微分の最小化)の方法として、初めてdamped molecular dynamics法とLagrangeの未定乗数法を採用したプログラムコードを開発したことである。この方法は単純ではあるが、非常に安定した波動関数の収束が得られると同時に、計算にかかる労力などは他のエネルギー最小化の方法と同様、無摂動系の計算と同程度で済むことが確認された。また、この方法は、非常に単純で、既存のCar-Parrinello分子動力学コードに密度汎関数摂動論を実装するのに非常に適したものであり、今後原子変位に限らない様々な摂動に対する応答関数の計算コードの開発の基礎となり得るものである。
|