This shows you the differences between two versions of the page.
quantumespresso:dft_u [2020/03/09 20:22] koudai [DFT計算] |
quantumespresso:dft_u [2021/06/27 22:04] |
||
---|---|---|---|
Line 1: | Line 1: | ||
- | ====== 概要 ====== | ||
- | 密度汎関数理論(density functional theory, DFT)ではバンドギャップの大きさが過小評価されたり、絶縁体であるにもかかわらず金属になってしまう場合もある。 | ||
- | このような問題を回避するため、オンサイトの相互作用を取り込んだDFT+Uの計算方法を解説し、ギャップの大きさの改善を目指す。 | ||
- | |||
- | このページのタイトルはDFT+Uとしているが、LDA+UおよびGGA+Uのことを総称してそう呼んでいる。 | ||
- | どちらの場合でも計算手順は同じである。 | ||
- | |||
- | なお、Uの値は手で入れないといけないが、どの程度の大きさか見積もる方法は別のページで触れる(予定)。 | ||
- | ただし、Uの大きさを第一原理的に見積もるための計算コストは大きい上にそんなに信頼性も高くないので、普通は実験で観測されたギャップや光吸収スペクトル、磁性体であれば磁気モーメントに合うように決める。 | ||
- | |||
- | |||
- | ====== NiO ====== | ||
- | |||
- | 酸化ニッケル(NiO)は典型的な反強磁性体である。 ただし、電子相関が強いため、通常のDFTによるバンド計算ではバンドギャップの大きさを説明できない。 | ||
- | 反強磁性のスピン配置については以下の論文を参照のこと。 | ||
- | |||
- | * W. L. Roth, Phys. Rev. 110, 1333 (1958). | ||
- | ====== 計算手順 ====== | ||
- | |||
- | |||
- | |||
- | |||
- | ===== DFT計算 ===== | ||
- | |||
- | まずはUを入れずに計算を行う。 | ||
- | 擬ポテンシャルはウルトラソフト型のO.pbe-rrkjus.UPFとNi.pbe-nd-rrkjus.UPFを使う。 | ||
- | 単位胞は反強磁性秩序を示すので2倍となる。 | ||
- | アップスピンのNiとダウンスピンのNiが非等価になることに注意する。 | ||
- | ここではそれぞれをNi1とNi2で分けて指定している。 | ||
- | |||
- | |||
- | <file - NiO.scf.in> | ||
- | & | ||
- | | ||
- | | ||
- | | ||
- | | ||
- | / | ||
- | &system | ||
- | | ||
- | | ||
- | | ||
- | ntyp= 3, | ||
- | | ||
- | | ||
- | | ||
- | | ||
- | | ||
- | | ||
- | | ||
- | / | ||
- | & | ||
- | / | ||
- | CELL_PARAMETERS | ||
- | 0.50 0.50 1.00 | ||
- | 0.50 1.00 0.50 | ||
- | 1.00 0.50 0.50 | ||
- | ATOMIC_SPECIES | ||
- | | ||
- | | ||
- | | ||
- | ATOMIC_POSITIONS crystal | ||
- | | ||
- | | ||
- | Ni1 0.0 0.0 0.0 | ||
- | Ni2 0.5 0.5 0.5 | ||
- | K_POINTS automatic | ||
- | 4 4 4 1 1 1 | ||
- | </ | ||
- | |||
- | スピン分極させる計算ではoccupations, | ||
- | 実行は次の通り。 | ||
- | |||
- | $ pw.x < NiO.scf.in > NiO.scf.out | ||
- | |||
- | 出力ファイルを見れば全磁化がゼロで絶対磁化が有限の値を持つので、反強磁性秩序を起こしていることがわかる。 | ||
- | |||
- | <file - NiO.scf.out> | ||
- | (略) | ||
- | total magnetization | ||
- | | ||
- | (略) | ||
- | </ | ||
- | |||
- | ひとつのNiが持つ磁化は1.56μ_Bで、実験で得られている値1.90μ_Bよりも小さい。 | ||
- | |||
- | * A. K. Cheetham and D. A. O. Hope, Phys. Rev. B 27, 6964 (1983). | ||
- | |||
- | 次に、バンドギャップを見積もるために状態密度を計算する。 | ||
- | 以下のNiO.nscf.inとNiO.dos.inを用意する。 | ||
- | NiO.nscf.inはNiO.scf.inをコピーして、calculationを' | ||
- | このときk点に対称点が含まれていたほうが状態密度が精度よく出る。 | ||
- | |||
- | <file - nio.nscf.in> | ||
- | & | ||
- | | ||
- | | ||
- | | ||
- | | ||
- | / | ||
- | &system | ||
- | | ||
- | | ||
- | | ||
- | ntyp= 3, | ||
- | | ||
- | | ||
- | | ||
- | | ||
- | | ||
- | / | ||
- | & | ||
- | / | ||
- | CELL_PARAMETERS | ||
- | 0.50 0.50 1.00 | ||
- | 0.50 1.00 0.50 | ||
- | 1.00 0.50 0.50 | ||
- | ATOMIC_SPECIES | ||
- | | ||
- | | ||
- | | ||
- | ATOMIC_POSITIONS crystal | ||
- | | ||
- | | ||
- | Ni1 0.0 0.0 0.0 | ||
- | Ni2 0.5 0.5 0.5 | ||
- | K_POINTS automatic | ||
- | 8 8 8 0 0 0 | ||
- | </ | ||
- | |||
- | <file - NiO.dos.in> | ||
- | &dos | ||
- | | ||
- | | ||
- | | ||
- | / | ||
- | </ | ||
- | |||
- | 実行は次のようにする。 | ||
- | |||
- | $ pw.x < NiO.nscf.in > NiO.nscf.out | ||
- | $ dos.x < NiO.dos.in > NiO.dos.out | ||
- | |||
- | 状態密度はnio.dosに出力される。 | ||
- | |||
- | |||
- | ===== DFT+U計算 ===== | ||
- | |||
- | 状態密度をプロットするとわかるが、バンドギャップは1eV程度しか開いていない。 | ||
- | しかしながら、実験で観測されているバンドギャップのサイズは4.3eVである。 | ||
- | |||
- | * G. A. Sawatzky and J. W. Allen, Phys. Rev. Lett. 53, 2339 (1984). | ||
- | |||
- | 電子相関の効果を取り入れることで、ギャップの値を再現する。 | ||
- | NiO.scf.inをコピーしてNiO_U.scf.inを作成する(ファイルが多くなるとややこしいので、新しいディレクトリで作業することを推奨)。 | ||
- | |||
- | <file - NiO_U.scf.in> | ||
- | & | ||
- | | ||
- | | ||
- | | ||
- | | ||
- | / | ||
- | &system | ||
- | | ||
- | | ||
- | | ||
- | ntyp= 3, | ||
- | | ||
- | | ||
- | | ||
- | | ||
- | | ||
- | | ||
- | | ||
- | | ||
- | | ||
- | | ||
- | / | ||
- | & | ||
- | / | ||
- | CELL_PARAMETERS | ||
- | 0.50 0.50 1.00 | ||
- | 0.50 1.00 0.50 | ||
- | 1.00 0.50 0.50 | ||
- | ATOMIC_SPECIES | ||
- | | ||
- | | ||
- | | ||
- | ATOMIC_POSITIONS crystal | ||
- | | ||
- | | ||
- | Ni1 0.0 0.0 0.0 | ||
- | Ni2 0.5 0.5 0.5 | ||
- | K_POINTS automatic | ||
- | 4 4 4 1 1 1 | ||
- | </ | ||
- | |||
- | さきほどとの違いは& | ||
- | それぞれのオプションの意味は次の通り。 | ||
- | |||
- | ^変数^説明^ | ||
- | |lda_plus_u|.true.にすることでDFT+U計算を行う。(LDA+U、GGA+Uで共通)| | ||
- | |Hubbard_U(i)|i番目の原子のオンサイトクーロン相互作用の大きさを表す。iは1からntypまでの値を取る。指定しなければゼロとなる。単位はeV。| | ||
- | |||
- | ここでは酸素の電子には相関は入れていない。 実行は次のようにする。 | ||
- | |||
- | $ pw.x < nio_u.scf.in > nio_u.scf.out | ||
- | |||
- | モーメントの値は原子あたり1.72μ_Bとなり、Uを入れることで改善された。 | ||
- | |||
- | <file - NiO_U.scf.out> | ||
- | (略) | ||
- | total magnetization | ||
- | | ||
- | (略) | ||
- | </ | ||
- | |||
- | 状態密度の計算はさきほどと同様である。 Uを入れることでバンドギャップが4eV程度となった。 | ||
- | |||
- | {{: |