====== 概要 ======


入門編の最後に、Niを例に磁性体の磁気モーメントを計算してみましょう。


====== SCF計算 ======

Niは金属なので、Alのときと同様に状態に広がりを持たせる必要があります。
擬ポテンシャルとしてウルトラソフト型の Ni.pz-nd-rrkjus.UPF を使います。

入力ファイルを作ります。

<file - Ni.scf.in >
&control
  prefix = 'nickel',
  pseudo_dir = './',
  outdir = './tmp/',
/
&system
  ibrav = 2, celldm(1) = 6.48, nat = 1, ntyp = 1,
  ecutwfc = 24.0, ecutrho = 192.0,
  occupations = 'smearing', smearing = 'mv', degauss = 0.02
  nspin = 2, starting_magnetization(1) = 0.7,
/
&electrons
/
ATOMIC_SPECIES
 Ni 58.69 Ni.pz-nd-rrkjus.UPF
ATOMIC_POSITIONS
 Ni 0.0 0.0 0.0
K_POINTS automatic
 8 8 8 1 1 1
</file>

新しく出てきたオプションの意味は次のとおりです。

^変数^説明^
|ecutrho |電荷密度の計算の際の運動エネルギーのカットオフ|
|nspin   |スピン分極しているかどうか。2でz軸方向に分極していることを表します。もし何も指定しなければ1となり磁化は起こりません。|
|starting_magnetization(i)|自己無撞着な計算において、開始時点でのi番目の原子の磁化。1ならすべてアップスピンで、−1ならすべてダウンスピンで開始することを意味します。0から始めてしまうと磁化がでないことに注意します。iは1からntypまでの値をとります。|

ecutrhoは何も指定しなければ自動的にecutwfcの4倍になりますが、ウルトラソフト型の擬ポテンシャルではさらに大きな値が要求されるので変更します。 

それではPWscfを実行します。
計算は金属の場合よりもさらに時間がかかりますが、今回は簡単な構造なのですぐに終わります。

  $ pw.x < Ni.scf.in > Ni.scf.out

得られた磁気モーメントは以下のように確認します。

  $ grep -e magnetization Ni.scf.out
      atomic species magnetization
      total magnetization = 1.85 Bohr mag/cell
      absolute magnetization = 1.87 Bohr mag/cell
      total magnetization = 0.66 Bohr mag/cell
      absolute magnetization = 0.78 Bohr mag/cell
      total magnetization = 1.00 Bohr mag/cell
      absolute magnetization = 1.10 Bohr mag/cell
      total magnetization = 0.62 Bohr mag/cell
      absolute magnetization = 0.74 Bohr mag/cell
      total magnetization = 0.62 Bohr mag/cell
      absolute magnetization = 0.71 Bohr mag/cell
      total magnetization = 0.62 Bohr mag/cell
      absolute magnetization = 0.69 Bohr mag/cell
      total magnetization = 0.62 Bohr mag/cell
      absolute magnetization = 0.69 Bohr mag/cell

一番下の2つの値が、SCF計算の収束によって得られた磁化です。
全磁化（total magnetization）$m_t$と絶対磁化（absolute magnetization）$m_a$の違いは、$n_{\uparrow}(r)$を位置$r$におけるアップスピンの電子数密度、$n_{\downarrow}(r)$をダウンスピンの電子数密度として

\begin{align}
m_t &= \int dr (n_{\uparrow}(r) − n_{\downarrow}(r)) \\
m_a &= \int dr |n_{\uparrow}(r) − n_{\downarrow}(r)|
\end{align}
で定義されます。 ただし積分範囲は単位胞内です。

強磁性体であれば全磁化と絶対磁化は同じになります。 一方で反強磁性体の場合は全磁化は0になりますが、絶対磁化は1つの原子が持つ磁化の2倍の値になります。

Niの場合は全磁化が実験で0.606 µBであることが知られていて、これに近い値が得られました。