User Tools

Site Tools


quantumespresso:磁性体

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

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