AkaiKKRでスーパーセル その1

AkaiKKRではコヒーレントポテンシャル近似(CPA)を用いて不純物を含んだ系の第一原理計算を行うことが出来ます。しかしながら、これは(他のコードでしばしば行われる)スーパーセルの計算が出来ないことを意味しません。
今回は計算機マテリアルデザイン入門 に載っている例であるFe87.5Li12.5合金の計算をスーパーセル法で行いました。

Chaparral_Supercell_520.jpg

Fig.1: スーパーセル(気象)。なお今回のエントリとは関係ない。



不純物の計算


多くの第一原理計算パッケージでは、不純物や格子欠陥の計算を行う際に、単位格子を複数個並べてその中の少数の原子を他の組成や格子欠陥に置き換えると言うことを行います。このように単位格子をたくさん並べたものをスーパーセルと、スーパーセルを使う計算手法をスーパーセル法と呼びます。

他方、AkaiKKRではコヒーレントポテンシャル近似(CPA)を用いることで、単位格子を用いながら不純物や格子欠陥の計算ができます。
CPAでは不規則性の効果によりバンド構造がにじむ(電子が有限の寿命を持つ)効果がはっきりと確認できます。
これに対して、スーパーセル法では、不純物周りの構造緩和を取り入れることが出来ます。

AkaiKKRはCPAが使えるという点が特徴的であるのですが、当然ながら、スーパーセルの計算ができないと言うわけではありません。(実際に構造緩和までやろうとすると大変だとは思いますが。)

今回は計算機マテリアルデザイン入門 に載っている体心立方構造(bcc)の鉄を2×2×2のスーパーセルとして、その中の1原子をリチウムに置き換えた計算を行います。

入力ファイル


以下に8個の鉄原子のうち1個をリチウムに置換したスーパーセルの入力ファイルを示します。
この他にスーパーセル法のテストのために純鉄の計算もスーパーセルで行います。

 go    data/fe_li_super
bcc 10.85, , , , , ,
0.001 1.2 nrl mjw mag 2nd
update 10 200 0.024
3
Fe1 1 0 0 2 26 100
Fe2 1 0 0 2 26 100
Li 1 0 0 1 3 100
8
0.00000 0.00000 0.00000 Li
0.50000 0.00000 0.00000 Fe1
0.00000 0.50000 0.00000 Fe1
0.00000 0.00000 0.50000 Fe1
0.25000 0.25000 0.25000 Fe2
0.25000 0.25000 0.75000 Fe2
0.25000 0.75000 0.25000 Fe2
0.75000 0.25000 0.25000 Fe2


比較のためにCPAを用いた不規則合金の計算も行います。

 go    data/fe_li
bcc 5.43, , , , , ,
0.001 1.2 nrl mjw mag 2nd
update 10 200 0.024
1
FeLi 2 0 0 2 26 87.5
3 12.5
1
0.00000 0.00000 0.00000 FeLi


大きなセルでのspecx.fの設定


今回のスーパーセルは2×2×2なので、原子の数は高々8個なのですが、計算する原子の数が増えるとspecx.fを編集して再makeする事が必要になる可能性があります。

パラメータの設定にはhow to run a system of over 30 atoms?AkaiKKRの角運動量(方位量子数)のメモが参考になるはずです。
以下に、赤井先生の回答とそれを私が意訳したものを示します。

It is still possible to run for a system with over 30 atoms.
Take care about the following:

1) Decrease the number of k-points in order that the memory addressing does not overflow (in the case of 32bit addressing system, such as Pentium 4 with 32bit addressing). This can be done by decreasing "nk1x". The test run should be done bu "bzqlty=0" or "bzqlty=1" for which nk1x=1~20 would be enough.

2) Increase the number aoms that can be handled:
For that increase "natmmx", "ncmpmx", and "msizmx".

For example if you have 30 atoms and for each atom you would like to calculate up to d states (l=2), all atoms have different types, each site has a single component (i.e. not a random alloy) the following parameters are suitable:
natmmx=30, ncmpmx=30, msizmx=270, mxlmx=3

30原子以上の系に関しても計算を行うことが出来ます。以下の点に注意してください。

1) メモリのオーバーフローを避けるためk点の数を減らします。そのためには"nk1x"を小さくします。"bzqlty=0"や"bzqlty=1"のテスト計算ではnk1x=1~20で充分です。

2) 原子数を増やします。これには"natmmx", "ncmpmx" 及び "msizmx"を大きくします。

例えば30原子をd状態までの計算で、全ての原子が異なる組成で、各サイト(原子位置)には1種類の組成しか持たない(つまり不規則合金でない)場合は、下記のパラメータとなります。
natmmx=30, ncmpmx=30, msizmx=270, mxlmx=3


計算結果


まずテスト計算としての純鉄の状態密度をスーパーセル法で計算したものを示します。
AkaiKKRでニッケル・鉄・コバルトなどの純鉄の状態密度と一致することが確認できます。

FeSuper.png

Fig.2: スーパーセル法により計算した鉄の状態密度


次にスーパーセル法を用いたFeLi合金とコヒーレントポテンシャル近似(CPA)を用いて計算したFeLi合金の状態密度を示します。

FeLiCPA.png

Fig.3: スーパーセル法(赤)とコヒーレントポテンシャル近似(緑)によるFe87.5Li12.5合金の状態密度


スーパーセル法による結果とCPAによる結果は、状態密度のおおよその形状に関しては非常に良く似ています。しかしながらCPAでは不規則性の効果で状態密度の鋭さが減じているのに対して、スーパーセル法ではとがったままです。

関連エントリ




参考URL




付録


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


参考文献/使用機器




フィードバック



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

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


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


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


tag: AkaiKKR machikaneyama KKR CPA スーパーセル 

comment

Secret

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

LTspiceAkaiKKRmachikaneyamaScilabKKRPSoCOPアンプPICCPA強磁性常微分方程式モンテカルロ解析odeトランジスタ状態密度インターフェーススイッチング回路ecaljPDS5022DOS定電流半導体シェルスクリプト乱数レベルシフトHP6632Aブレッドボード分散関係温度解析R6452Aトランジスタ技術I2C可変抵抗反強磁性セミナー数値積分確率論偏微分方程式バンド構造非線形方程式ソルババンドギャップ絶縁熱設計シュミットトリガLEDA/Dコンバータ三端子レギュレータLM358ISO-I2CGW近似カオスフォトカプラマフィンティン半径TL431数値微分PC817Cアナログスイッチ直流動作点解析発振回路USBサーボカレントミラー74HC4053パラメトリック解析LDAbzqltyチョッパアンプ量子力学FFT2ちゃんねるアセンブラBSch開発環境電子負荷ブラべ格子イジング模型補間基本並進ベクトル標準ロジック単振り子キュリー温度繰り返しMaxima状態方程式失敗談相対論スピン軌道相互作用FETランダムウォーク熱伝導六方最密充填構造コバルトewidthTLP621GGAQSGW不規則合金位相図抵抗SMPcygwinラプラス方程式スレーターポーリング曲線gfortranスイッチト・キャパシタ詰め回路TLP552三角波格子比熱TLP521条件分岐LM555MCUNE555QNAPマントルテスタ過渡解析FXA-7020ZRダイヤモンドデータロガーガイガー管自動計測Writer509UPSシュレディンガー方程式ブラウン運動awk差し込みグラフ熱力学平均場近似仮想結晶近似VCAfsolve井戸型ポテンシャルVESTA起電力スーパーセルOpenMP第一原理計算ubuntu固有値問題L10構造OPA2277interp12SC1815fccウィグナーザイツ胞面心立方構造フィルタジバニャン方程式ヒストグラム確率論マテリアルデザインspecx.f等高線正規分布PGAフェルミ面非線型方程式ソルバ初期値固定スピンモーメントスワップ領域ルチル構造リジッドバンド模型edeltquantumESPRESSO岩塩構造BaOSIC二相共存ZnOウルツ鉱構造フォノンデバイ模型c/aノコギリ波全エネルギーFSMTeXgnuplotmultiplotハーフメタルCapSense半金属合金結晶磁気異方性Ubuntu文字列入出力TS-110TS-112疎行列Excel直流解析ヒストグラム円周率不規則局所モーメントトラックボールPC等価回路モデルパラメータ・モデルキーボードRealforce三次元マンデルブロ集合フラクタル化学反応重積分縮退日本語最小二乗法関数フィッティングGimpMAS830LHiLAPW熱拡散方程式両対数グラフナイキスト線図負帰還安定性陰解法Crank-Nicolson法P-10クーロン散乱境界条件連立一次方程式片対数グラフEAGLEPIC16F785LMC662トランスシンボルCK1026線種凡例MBEAACircuitグラフの分割軸ラベルifort

最新コメント
リンク

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