11.8 fe.bcc(Crystalline Magnetism：結晶磁性)

The test calculation of body-centered cubic lattice(bcc)] in the magnetic state is performed.

体心立方格子(bcc)での磁性状態での計算チェックを行う。

Property（性質）
crystal structure	body centered cubic lattic IBRAV=3
space group	No.229 $F \frac{4}{m} \bar{3} \frac{2}{m}, O_{h}^{9}$
lattice constant	5.23 a.u.
xc potentail	CA-PZ(LSDA)
energy cut-offs	40 Ry, 250 Ry
k point sampling	Monkhost & Pack $12\times12\times12$, NKPT=68

It is specified that calculation of the bulk modulus is able to be performed as an optional extra. If you change the pseudopotential from fe_sr.app.fmt to uspp.Fe-ca-sp-oda-irel-2.fmt, you can confirm that pressure value and, an equilibrium lattice constant and the bulk modulus estimated from it almost correspond to them estimated from the energy change. The former pseudopotential seemes to be used for the calculation, and the latter pseudopotinal seems to be used for the crystal calculation. Both are different in the energy which decides the logarithmic derivative of components of $d$-orbital. The former has a narrow energy width, and the latter has a wide energy width. Additionally, to obtain exact value, the variation width of the volume is set to a little small value.

オプションとして、体積弾性率の計算ができるように設定されている。 擬ポテンシャルファイルをfe_sr.app.fmtからuspp.Fe-ca-sp-oda-irel-2.fmtへ変更して、 圧力値とそこから見積もられる平衡格子定数と体積弾性率が、エネルギーの変化から見積もられるものとほぼ一致する ことが確かめられた。前者の擬ポテンシャルファイルは、クラスタ用として使用し、後者は結晶用として使用するのが よいと思われる。両者は、$d$軌道成分の対数微分を決めるエネルギー値が異なっており、前者の場合はエネルギー幅 が狭くなっており、後者の場合はエネルギー幅が広くなっている。また、エネルギーから求める体積弾性率の値を正確 に求めるため、体積の変化分を少し小さくとっている。