ナノチューブの電子状態の研究 (1998) ミーティング報告

10/20 の内容

・overlap matrix S を考慮した時のリカージョン法

u_nSu_m = δ_nm である。

Η_mn と S は与えられているので、

次のように計算しておく必要がある。

H = S^-1[ Η_mn ]

必要に応じて、 H を用いる。

・炭素クラスタのT.B.ハミルトニアンの計算と状態密度計算プログラム

C₆₀ の状態密度の、計算結果が MOPAC の計算と似ないのは mopac 計算の性質かも知れない。

フェルミエネルギーの付近の状態密度を調べるべき。

・周期系の構造最適化プログラム

グラフファイトの２次元構造、ダイヤモンド構造を再現出来た。

最近接

このプログラムでは、グラフファイト AA スタック構造( ３次元) は、安定だった。

AB スタック構造はダイヤモンド構造に... とおもったら初期構造がABC スタックだった。

※後の計算ではどちらも層間 2.1Åで安定。

・分子軌道の数が少なくなっているのに計算が成り立つのは何故か?

系のハミルトニアンの対称性が関係している。

リカージョン法では、最初の波動関数の組を対称性が高くなるようにとることで、計算が良くなる。

課題: 1. C₆₀ の状態密度をリカージョン法で計算; 2. C₆₀ のフェルミエネルギーの付近の状態密度; 3. 疎行列を帯状列にするプログラム; 4. グラファイトの面間を再現できそうなポテンシャルを考える; 5. 正方格子の周期系の状態密度の再計算

コメントまたはアドバイスなどがあれば以下のアドレスへどうぞ。