quantumespresso:磁性体
Differences
This shows you the differences between two versions of the page.
| Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
| quantumespresso:磁性体 [2018/11/02 19:27] – [SCF計算] koudai | quantumespresso:磁性体 [2021/06/27 22:04] (current) – external edit 127.0.0.1 | ||
|---|---|---|---|
| Line 8: | Line 8: | ||
| Niは金属なので、Alのときと同様に状態に広がりを持たせる必要があります。 | Niは金属なので、Alのときと同様に状態に広がりを持たせる必要があります。 | ||
| - | 擬ポテンシャルとしてNi.pz-nd-rrkjus.UPFを使います。 | + | 擬ポテンシャルとしてウルトラソフト型の |
| 入力ファイルを作ります。 | 入力ファイルを作ります。 | ||
| Line 39: | Line 39: | ||
| |ecutrho |電荷密度の計算の際の運動エネルギーのカットオフ| | |ecutrho |電荷密度の計算の際の運動エネルギーのカットオフ| | ||
| |nspin | |nspin | ||
| - | |starting_magnetization(i)|自己無撞着な計算において、開始時点でのi番目の原子の磁化。1ならすべてアップスピンで、−1ならすべてダウンスピンで開始することを意味します。0から始めてしまうと磁化がでないことに注意します。iは1からnatまでの値をとります。| | + | |starting_magnetization(i)|自己無撞着な計算において、開始時点でのi番目の原子の磁化。1ならすべてアップスピンで、−1ならすべてダウンスピンで開始することを意味します。0から始めてしまうと磁化がでないことに注意します。iは1からntypまでの値をとります。| |
| - | k点の数も磁性体では大きい値が要求されますが、通常の2倍くらいにしておけば問題ありません。 | + | ecutrhoは何も指定しなければ自動的にecutwfcの4倍になりますが、ウルトラソフト型の擬ポテンシャルではさらに大きな値が要求されるので変更します。 |
| - | 計算は金属の場合よりもさらに時間がかかりますが、今回は簡単な構造なのですぐに終わります。 | + | |
| それではPWscfを実行します。 | それではPWscfを実行します。 | ||
| + | 計算は金属の場合よりもさらに時間がかかりますが、今回は簡単な構造なのですぐに終わります。 | ||
| - | $ pw.x < ni.scf.in > ni.scf.out | + | $ pw.x < Ni.scf.in > Ni.scf.out |
| 得られた磁気モーメントは以下のように確認します。 | 得られた磁気モーメントは以下のように確認します。 | ||
| - | $ grep -e magnetization | + | $ grep -e magnetization |
| atomic species magnetization | atomic species magnetization | ||
| total magnetization = 1.85 Bohr mag/cell | total magnetization = 1.85 Bohr mag/cell | ||
quantumespresso/磁性体.1541154459.txt.gz · Last modified: 2021/06/27 22:00 (external edit)