磁気相互作用
拡張機能をご使用になりたい場合はアカデメイアまでお問い合わせ下さい。
概要
交換結合定数Jijからキュリー温度を計算する方法について、AkaiKKRには有限のrange内でJijの和を取る”jモード”と、逆格子空間でJijの和をとる(無限遠点までの和を取る)”tcモード”の2つのモードが実装されている。以下、それら2つのモードについて説明する。
jモードによるキュリー温度の計算
“jモード”でのキュリー温度の計算は、in/feファイル内の”go”を”j”に書き換えることで行える。
c----------------------Fe------------------------------------
c go file
j data/fe
c------------------------------------------------------------
c brvtyp a c/a b/a alpha beta gamma
bcc 5.27 , , , , , ,
c------------------------------------------------------------
c edelt ewidth reltyp sdftyp magtyp record
0.001 1.0 nrl mjw mag 2nd
c------------------------------------------------------------
c outtyp bzqlty maxitr pmix
update 4 50 0.023
c------------------------------------------------------------
c ntyp
1
c------------------------------------------------------------
c type ncmp rmt field mxl anclr conc
Fe 1 1 0.0 2
26 100
c------------------------------------------------------------
c natm
1
c------------------------------------------------------------
c atmicx type
0 0 0 Fe
c------------------------------------------------------------
入力ファイルの準備ができたら、goによる計算と同様の方法で計算を実行できる。
$ ./specx < in/fe
goによる計算と同様の結果が表示された後、Jijの出力とキュリー温度の出力が得られる。Jijについてはrange内の区別できるペアのリストが表示されている。dgnは等価なペアの数である。
J_ij
Fe-Fe
index site comp cell distance J_ij J_ij(meV) dgn
1 1 1 1 1 0.5000 0.5000 -0.5000 0.866025 0.001253 17.045871 8
2 1 1 1 1 1.0000 0.0000 0.0000 1.000000 0.000790 10.752440 6
3 1 1 1 1 0.0000 -1.0000 1.0000 1.414214 0.000154 2.101422 12
Tc (in mean field approximation) = 1749.175K
tcモードによるキュリー温度の計算
“tcモード”でのキュリー温度の計算は、in/feファイル内の”go”を”tc”に書き換えることで行える。
c----------------------Fe------------------------------------
c go file
tc data/fe
c------------------------------------------------------------
c brvtyp a c/a b/a alpha beta gamma
bcc 5.27 , , , , , ,
c------------------------------------------------------------
c edelt ewidth reltyp sdftyp magtyp record
0.001 1.0 nrl mjw mag 2nd
c------------------------------------------------------------
c outtyp bzqlty maxitr pmix
update 4 50 0.023
c------------------------------------------------------------
c ntyp
1
c------------------------------------------------------------
c type ncmp rmt field mxl anclr conc
Fe 1 1 0.0 2
26 100
c------------------------------------------------------------
c natm
1
c------------------------------------------------------------
c atmicx type
0 0 0 Fe
c------------------------------------------------------------
入力ファイルの準備ができたら、goによる計算と同様の方法で計算を実行できる。
$ ./specx < in/fe
goによる計算と同様の結果が表示された後、分子場の値(Hmf。mfはmolecular fieldの略。)とキュリー温度が得られる。
Tc calculation starts
Molecular field on each sublattice in
the assumed magnetic configuration
sublat cmpnt Hmf(Ry) Hmf(meV)
-------------------------------------
1 1 0.026992 367.241361
Tc (in mean field approximation) = 1420.551K