スポンサーサイト

上記の広告は1ヶ月以上更新のないブログに表示されています。
新しい記事を書く事で広告が消せます。


AkaiKKRで反強磁性クロム

AkaiKKR(machikaneyama)を用いて整合反強磁性クロムの計算を行いました。反強磁性の計算の正しいやり方は分かりませんが、格子定数を大きくして無理矢理に強磁性状態にしたクロムのポテンシャルファイルを初期ポテンシャルとして計算を行うと、強磁性状態とも常磁性状態とも異なる反強磁性状態の解が得られませした。
しかしながら、全エネルギーが最小となる格子定数においては、反強磁性状態が安定とはなりませんでした。

CrDOS.gif
Fig.1: 反強磁性クロム(赤)と常磁性クロム(緑)の部分状態密度の格子体積依存性。格子体積が小さくなるとともに反強磁性的な磁気モーメントが消える。



反強磁性クロム


体心立方構造(bcc)のクロムは、反強磁性体の金属として知られています。Wikipediaの反強磁性の項目には以下のように、クロムを広義の反強磁性体として分類しています。

体心立方構造の頂点位置と体心位置のスピンが反対方向を向いているだけでなく、スピンの大きさが単位胞ごとに正弦関数的に変化している。

しかしながら、第一原理計算では、しばしば単純に体心と頂点のスピンが反対を向いているだけの反強磁性クロムの計算がなされます。(例えば第一原理計算ソフトウェア Advance/PHASEの応用機能と解析事例計算機マテリアルデザイン入門 (大阪大学新世紀レクチャー)のABCAPによる計算例など。)このような反強磁性を整合(commensurate)反強磁性と呼びます。



今回はAkaiKKR(machikaneyama)を用いて上記の単純な整合反強磁性クロムの計算を行いうことを目標にします。

反強磁性体のポテンシャルファイル


クロムの入力ファイルを作成し、そのままmagで計算を行っても反強磁性状態の解は得られません。

計算機マテリアルデザイン入門 (大阪大学新世紀レクチャー)のAkaiKKRの項目には、強磁性状態の鉄のポテンシャルファイルからfmgを用いてスピングラス状態のポテンシャルファイルを作成する方法が書かれています。さらにそのページには反強磁性状態の計算をする場合にも同じようにしてデータファイルを作ることができるとあります。

スピングラス状態の計算では、強磁性状態の計算で得られたポテンシャルファイルのスピンの向きを反転させたポテンシャルファイルを作成し、それらをCPAする事でスピングラス状態を実現していました。いうなればスピンの異なる原子の不規則合金です。ねがてぃぶろぐでもAkaiKKRで鉄のキュリー温度にて計算を行いました。
このことから類推するに、反強磁性状態も、スピンの異なる原子を交互に並べることで実現できそうです。つまりスピンの異なる原子の規則合金という事です。

そのためにはまず、強磁性状態のポテンシャルファイルを作成します。bcc構造のクロムの格子定数はa=2.88Å(a=5.45bohr)程度ですが、この格子定数で磁性状態を含む計算(mag)を行っても、常磁性状態の解が得られるのみです。
そこで意図的に大きな格子定数(例えばa=8.0bohrなど)を指定して強磁性解を得ておきます。

その後、utilにあるfmgを用いて、反転させたポテンシャルファイルを初期ポテンシャルとしてgo計算を行います。
シェルスクリプトやfmgのための入力ファイル、計算のテンプレートなどは以下のようになりました。実行ファイルfmgはパスの通ったところにおいてあるとし、入力ファイルのテンプレートはtemplateという名前のフォルダを作成して入れておきます。

#!/bin/csh -f

set A_LIST=( 5.60 5.55 5.50 5.45 5.40 5.35 5.30 5.25 5.20 )

rm analysis/CrAFM.dat
rm analysis/CrNMG.dat

specx < in/CrFMG.in > out/CrFMG.out
fmg < crafm

foreach A ( ${A_LIST} )
sed 's/'ABOHR'/'${A}'/g' template/CrAFM0.in > in/CrAFM_${A}.in
specx < in/CrAFM_${A}.in > out/CrAFM_${A}.out
sed 's/'ABOHR'/'${A}'/g' template/CrNMG0.in > in/CrNMG_${A}.in
specx < in/CrNMG_${A}.in > out/CrNMG_${A}.out

sed -e '$d' data/crafm.info | sed -n '$p' >> analysis/CrAFM.dat
sed -e '$d' data/crnmg.info | sed -n '$p' >> analysis/CrNMG.dat
end

grep "spin moment" out/CrAFM_*.out | awk 'NR % 4 == 1;NR % 4 == 2'




結果


ある一つのクロム原子に着目した部分状態密度(PDOS)の格子体積依存性がFig.1のgif動画です。格子体積が大きく、反強磁性的な磁気モーメントが存在する間はDOSの形状が非磁性状態のものと異なっていますが、体積が小さくなり磁気モーメントが消えるとともにDOSの形状も非磁性状態のものと同じになります。

DOS.png
Fig.3: 格子定数a=5.45bohrにおける反強磁性クロム(赤)と常磁性クロム(緑)の部分状態密度


Fig.3に示したのが実験により得られた格子定数付近(a=5.45bohr)における反強磁性クロムと常磁性クロムの部分状態密度です。反強磁性状態ではフェルミ準位近傍に特徴的なふくらみを見ることができます。

CrAFM-Energy.png
Fig.4: 全エネルギーの格子体積依存性


Fig.4に示したのが全エネルギーの格子体積依存性です。磁気モーメントが大きい間は、反強磁性相の方がエネルギー的に安定であることがわかります。磁気モーメントが小さくなるにしたがって全エネルギーの差が小さくなることも読み取れます。
しかし、全エネルギーが最小化するのはほとんど磁気モーメントが消失しているa=5.40bohr付近でエネルギー的にはほとんど差が見られません。それどころか、残念なことにわずかながら常磁性状態の方がエネルギー的に安定であるという結果になっています。

これが整合反強磁性を仮定したためなのか、単純にbzqlty(k点の数)が低いためなのか、数値計算の限界なのかは私には判断できません。

関連エントリ




参考URL




付録


このエントリで使用したAkaiKKR用のファイルを添付します。ファイル名末尾の".txt"を削除して、"_"を"."に変更すれば使えるはずです。(参考:ねがてぃぶろぐの付録)


参考文献/使用機器




フィードバック



にほんブログ村 その他趣味ブログ 電子工作へ

 ↑ 電子工作ブログランキング参加中です。1クリックお願いします。


コメント・トラックバックも歓迎です。 ↓      


 ↓ この記事が面白かった方は「拍手」をお願いします。


tag: AkaiKKR machikaneyama KKR 反強磁性 DOS 状態密度 

comment

Secret

FC2カウンター
カテゴリ
ユーザータグ

LTspiceAkaiKKRmachikaneyamaScilabKKRPSoC強磁性CPAPICOPアンプecalj状態密度モンテカルロ解析常微分方程式トランジスタodeDOSインターフェースPDS5022定電流スイッチング回路分散関係半導体シェルスクリプトレベルシフト乱数HP6632A可変抵抗トランジスタ技術R6452AI2C温度解析ブレッドボードバンドギャップ確率論反強磁性セミナーバンド構造数値積分偏微分方程式非線形方程式ソルバ熱設計絶縁三端子レギュレータISO-I2CA/DコンバータシュミットトリガフォトカプラカオスPWscfGW近似LM358LEDマフィンティン半径発振回路USB数値微分TL431PC817Cサーボアナログスイッチ直流動作点解析補間カレントミラー74HC4053bzqltyチョッパアンプFFT2ちゃんねる開発環境量子力学単振り子電子負荷VESTAQuantumESPRESSO標準ロジックパラメトリック解析ブラべ格子イジング模型アセンブラLDA基本並進ベクトルBSchSMPTLP621失敗談六方最密充填構造コバルト位相図QSGWGGAスイッチト・キャパシタewidth状態方程式VCAキュリー温度繰り返し最適化仮想結晶近似不規則合金熱伝導gfortran相対論抵抗FETMaximaQuantum_ESPRESSOcygwinランダムウォークラプラス方程式スピン軌道相互作用スレーターポーリング曲線マントルシュレディンガー方程式ZnO自動計測QNAP固有値問題ダイヤモンドデータロガー井戸型ポテンシャルTLP552CIFxcrysdenゼーベック係数熱力学条件分岐MCU最小値UPS格子比熱最大値ガイガー管平均場近似過渡解析Writer509スーパーセルFXA-7020ZR差し込みグラフ第一原理計算テスタ起電力OpenMP三角波ubuntuLM555NE555ブラウン運動詰め回路ハーフメタルawkfsolveUbuntuフェルミ面TLP521トランスMAS830LPGACK1026OPA2277フィルタトレーナーバトルEAGLEノコギリ波負帰還安定性ナイキスト線図MBEP-10LMC6622SC1815CapSenseAACircuitPIC16F785入出力固定スピンモーメントFSMTeX結晶磁気異方性全エネルギーc/a合金multiplotgnuplot非線型方程式ソルバL10構造正規分布等高線ジバニャン方程式ヒストグラム確率論初期値interp1fcc面心立方構造ウィグナーザイツ胞半金属デバイ模型磁気モーメント電荷密度重積分SIC不純物問題擬ポテンシャル状態図cif2cellPWgui二相共存ウルツ鉱構造edeltquantumESPRESSOフォノンリジッドバンド模型スワップ領域BaO岩塩構造ルチル構造マテリアルデザインspecx.fフラクタルマンデルブロ集合キーボードRealforceクーロン散乱三次元疎行列縮退化学反応関数フィッティング最小二乗法Excel直流解析PCTS-110TS-112日本語パラメータ・モデル等価回路モデル文字列ポケモンGO熱拡散方程式HiLAPW両対数グラフ片対数グラフ陰解法Crank-Nicolson法ifort境界条件連立一次方程式グラフの分割軸ラベルヒストグラム不規則局所モーメントスーパーリーグ円周率Gimp凡例線種シンボルトラックボール

最新コメント
リンク

にほんブログ村 その他趣味ブログ 電子工作へ
上記広告は1ヶ月以上更新のないブログに表示されています。新しい記事を書くことで広告を消せます。