This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
wien2k:反強磁性 [2020/07/03 17:17] koudai [LDA+U] |
wien2k:反強磁性 [2021/07/22 18:38] (current) koudai [WIEN2kを使う場合] |
||
---|---|---|---|
Line 17: | Line 17: | ||
(反強磁性相転移をしても単位胞が拡大しない場合はこの操作は不要です) | (反強磁性相転移をしても単位胞が拡大しない場合はこの操作は不要です) | ||
+ | |||
+ | ==== WIEN2kを使う場合 ==== | ||
- supercellというプログラムを使って超格子を作成します。ここでは単位胞を各方向に2倍して、面心立方格子から単純立方格子に変更しています< | - supercellというプログラムを使って超格子を作成します。ここでは単位胞を各方向に2倍して、面心立方格子から単純立方格子に変更しています< | ||
Line 57: | Line 59: | ||
</ | </ | ||
* 超格子のstructファイルである NiO_super.struct ができます。 | * 超格子のstructファイルである NiO_super.struct ができます。 | ||
+ | |||
+ | |||
+ | ==== VESTAを使う場合 ==== | ||
+ | |||
+ | - structファイルをVESTAで開きます | ||
+ | - メニューから [Edit] -> [Edit Data] -> [Unit cell] -> [Transform] | ||
+ | - 開いたウィンドウの右に書いてある説明を参考に、変換行列を指定します | ||
+ | * 今回の場合は対角成分をすべて2にしてください | ||
+ | - [Search Atoms ...] を選択します | ||
===== スピン配置の指定 ===== | ===== スピン配置の指定 ===== | ||
Line 129: | Line 140: | ||
これは実験で得られているもの(手法にもよりますが4eV前後)よりもだいぶ小さいです。 | これは実験で得られているもの(手法にもよりますが4eV前後)よりもだいぶ小さいです。 | ||
* 参考: https:// | * 参考: https:// | ||
- | * ギャップに関してはk点数に関して収束していないので、増やしたほうが良いです。どのみち正確ではないのであまり意味がありませんが... | + | * このページの例だとギャップの大きさはk点数に対して収束していないです。どのみち正確ではないので増やしてもあまり意味がありませんが... |
実は第一原理計算は電子相関を無視した計算になっているため、電子相関の強いNiOではこのように実験に合わない結果となります。 | 実は第一原理計算は電子相関を無視した計算になっているため、電子相関の強いNiOではこのように実験に合わない結果となります。 | ||
Line 141: | Line 152: | ||
事前にスピン分極のSCF計算を済ませてください。 | 事前にスピン分極のSCF計算を済ませてください。 | ||
- | - LDA+U計算用の入力ファイルを作成します。今回はNiのd軌道に U = 6 eV (6/13.6 = 0.441...) のポテンシャルを入れます。Jはゼロにしました。< | + | - LDA+U計算用の入力ファイルを作成します。今回はNiのd軌道(l=2)に U = 6 eV (6/13.606 = 0.4409...) のポテンシャルを入れます。Jはゼロにしました。< |
$ init_orb_lapw -orb | $ init_orb_lapw -orb | ||
You have the following atoms: | You have the following atoms: | ||
Line 148: | Line 159: | ||
3 : O | 3 : O | ||
3 : O | 3 : O | ||
- | Enter the name, l, U(Ry) and J of the Atoms (eg. Fe 2 0.3 0.0; exit with RETURN):Ni 2 0.44 0 | + | Enter the name, l, U(Ry) and J of the Atoms (eg. Fe 2 0.3 0.0; exit with RETURN):Ni 2 0.441 0 |
Enter the name, l, U(Ry) and J of the Atoms (eg. Fe 2 0.3 0.0; exit with RETURN): | Enter the name, l, U(Ry) and J of the Atoms (eg. Fe 2 0.3 0.0; exit with RETURN): | ||
</ | </ | ||
- | - 入力ファイルが表示されるので、正しく設定されたか確認します。割と間違って出力されることが多いです。 | + | - 入力ファイルが表示されるので、正しく設定されたか確認します |
* NiO_AFM.indm | * NiO_AFM.indm | ||
* 2行目: +Uを入れる原子の数 | * 2行目: +Uを入れる原子の数 | ||
Line 176: | Line 187: | ||
2 1 2 index of atom, number of l, l | 2 1 2 index of atom, number of l, l | ||
1 nsic 0..AMF, 1..SIC, 2..HFM | 1 nsic 0..AMF, 1..SIC, 2..HFM | ||
- | 0.44 0 U J (Ry) | + | 0.441 0 U J (Ry) |
- | 0.44 0 U J (Ry) | + | 0.441 0 U J (Ry) |
</ | </ | ||
* https:// | * https:// | ||
* ブロイデン法の履歴を削除して、オプション-orbをつけて実行します< | * ブロイデン法の履歴を削除して、オプション-orbをつけて実行します< | ||
$ rm *.broyd* | $ rm *.broyd* | ||
- | $ runsp_lapw -cc 0.0001 -ec 0.00001 -i 100 -orb -p | + | $ runsp_lapw -cc 0.0001 -ec 0.00001 -i 100 -orb -p -NI |
</ | </ | ||
* k点数やUをかえて再度実行するときは、ブロイデン法の履歴を削除しなくてもオプション -NI をつけるだけで十分です。 | * k点数やUをかえて再度実行するときは、ブロイデン法の履歴を削除しなくてもオプション -NI をつけるだけで十分です。 |