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アンプCPA強磁性PICモンテカルロ解析常微分方程式odeトランジスタecalj状態密度DOSインターフェース定電流スイッチング回路PDS5022半導体シェルスクリプト乱数レベルシフトHP6632A温度解析分散関係I2Cトランジスタ技術R6452A可変抵抗ブレッドボードセミナーバンドギャップ数値積分確率論反強磁性偏微分方程式バンド構造絶縁熱設計非線形方程式ソルバフォトカプラシュミットトリガLEDLM358カオスISO-I2C三端子レギュレータGW近似A/Dコンバータカレントミラーアナログスイッチ数値微分マフィンティン半径TL431発振回路サーボPC817CUSB直流動作点解析74HC4053補間FFTBSch開発環境パラメトリック解析2ちゃんねるチョッパアンプ量子力学bzqlty電子負荷イジング模型LDA標準ロジックアセンブラ基本並進ベクトルブラべ格子単振り子熱伝導位相図TLP621キュリー温度繰り返し状態方程式MaximaVESTAスイッチト・キャパシタ相対論FETランダムウォークスピン軌道相互作用SMP六方最密充填構造抵抗不規則合金ewidthスレーターポーリング曲線GGAラプラス方程式cygwingfortranQSGW失敗談コバルト条件分岐TLP521テスタLM555Writer509TLP552格子比熱マントルデータロガー自動計測詰め回路ガイガー管ダイヤモンドQNAPMCUFXA-7020ZR過渡解析三角波UPSNE555固有値問題熱力学ブラウン運動フェルミ面awk起電力第一原理計算OpenMPfsolveubuntu最大値xcrysden最小値最適化仮想結晶近似VCA差し込みグラフスーパーセル井戸型ポテンシャル平均場近似シュレディンガー方程式FSMフラクタルOPA2277固定スピンモーメント2SC1815全エネルギー合金multiplotgnuplotc/aTeX結晶磁気異方性interp1ウィグナーザイツ胞初期値マンデルブロ集合疎行列面心立方構造fcc不純物問題非線型方程式ソルバフィルタL10構造PGA半金属二相共存SICZnOウルツ鉱構造BaO重積分クーロン散乱磁気モーメント電荷密度三次元CIF岩塩構造CapSenseノコギリ波デバイ模型ハーフメタル正規分布フォノンquantumESPRESSOルチル構造スワップ領域リジッドバンド模型edelt縮退キーボード軸ラベルグラフの分割凡例トラックボールPC不規則局所モーメント片対数グラフトランス両対数グラフCK1026MAS830L直流解析Excel円周率パラメータ・モデルヒストグラム日本語最小二乗法等価回路モデルGimp線種シンボルTS-110TS-112PIC16F785LMC662化学反応文字列specx.f入出力ifortマテリアルデザインヒストグラム確率論Realforce等高線ジバニャン方程式P-10Ubuntuナイキスト線図Crank-Nicolson法陰解法熱拡散方程式HiLAPWAACircuit連立一次方程式負帰還安定性境界条件EAGLEMBE関数フィッティング

最新コメント
リンク

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