スポンサーサイト

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


AkaiKKRで点欠陥の形成エネルギー(未完)

Wang et al. (2004) (PDF)を参考に、AkaiKKR (machikaneyama)を用いて面心立方構造(fcc)のアルミニウムの点欠陥の生成エネルギーをスーパーセル法を用いて計算しました。
結果は、Wang et al. (2004) (PDF)の計算結果に対して、一桁程度の過大評価となってしまいました。


点欠陥の形成エネルギー


結晶中に点欠陥を作るために必要な欠陥の形成エネルギーは、スーパーセル法を用いた第一原理計算から計算することができます。
具体的には、まず、スーパーセルで完全結晶の全エネルギーを計算します(E0)。次にスーパーセルから1つの原子を取り除いた系の計算を行い、全エネルギー(E1)を求めます。スーパーセルの原子数がN個のとき、欠陥形成エネルギーは、以下の式から求めることができます(参考: Wang et al. (2004) (PDF))。
\begin{equation}
\Delta E_{f} = E_{1} - \frac{N-1}{N} E_{0}
\end{equation}

同様の計算は、スーパーセルではなくコヒーレントポテンシャル近似(CPA)でもできる可能性はあるのでしょうか?このような質問がAkaiKKR (machikaneyama)の掲示板に投稿されて・・・いましたが、現在は存在していないようです(#6678)。質問者の計算は、うまく行っておらず、形成エネルギーが一桁程度過大評価されているようです。

そこで今回は、AkaiKKR (machikaneyama)でもスーパーセル法を用いて点欠陥の形成エネルギーを計算してみます。

計算手法


Wang et al. (2004) (PDF)では、面心立方構造(fcc)のアルミニウムとニッケル、体心立方構造(bcc)のモリブデンとタンタルの点欠陥の形成エネルギーの計算が行われています。今回は、この中で面心立方構造のアルミニウムに関して計算を行います。スーパーセルのサイズは 2*2*2=32 原子としました。
Wang et al. (2004) (PDF)によると、交換相関汎関数はGGAよりもLDAのほうが良いであるとか、構造緩和はしないほうがむしろ良いだとか、色々議論があるようです。とりあえず今回は、LDA(mjw)で構造緩和もしない(というか大変なのでやりたくない)という方針で行きます。

c--------------------Al--------------------------------------
go data/SuperAlVc
c------------------------------------------------------------
c brvtyp a c/a b/a alpha beta gamma
sc 15.3 , , , , , ,
c------------------------------------------------------------
c edelt ewidth reltyp sdftyp magtyp record
0.001 1.0 sra mjw nmag 2nd
c------------------------------------------------------------
c outtyp bzqlty maxitr pmix
update 4 200 0.035
c------------------------------------------------------------
c ntyp
32
c------------------------------------------------------------
c type ncmp rmt field mxl anclr conc
Al1 1 1 0.0 2 0 100
Al2 1 1 0.0 2 13 100
Al3 1 1 0.0 2 13 100
Al4 1 1 0.0 2 13 100
Al5 1 1 0.0 2 13 100
Al6 1 1 0.0 2 13 100
Al7 1 1 0.0 2 13 100
Al8 1 1 0.0 2 13 100
Al9 1 1 0.0 2 13 100
Al10 1 1 0.0 2 13 100
Al11 1 1 0.0 2 13 100
Al12 1 1 0.0 2 13 100
Al13 1 1 0.0 2 13 100
Al14 1 1 0.0 2 13 100
Al15 1 1 0.0 2 13 100
Al16 1 1 0.0 2 13 100
Al17 1 1 0.0 2 13 100
Al18 1 1 0.0 2 13 100
Al19 1 1 0.0 2 13 100
Al20 1 1 0.0 2 13 100
Al21 1 1 0.0 2 13 100
Al22 1 1 0.0 2 13 100
Al23 1 1 0.0 2 13 100
Al24 1 1 0.0 2 13 100
Al25 1 1 0.0 2 13 100
Al26 1 1 0.0 2 13 100
Al27 1 1 0.0 2 13 100
Al28 1 1 0.0 2 13 100
Al29 1 1 0.0 2 13 100
Al30 1 1 0.0 2 13 100
Al31 1 1 0.0 2 13 100
Al32 1 1 0.0 2 13 100
c------------------------------------------------------------
c natm
32
c------------------------------------------------------------
c atmicx atmtyp
0 0 0 Al1
1/4 1/4 0 Al2
1/4 0 1/4 Al3
0 1/4 1/4 Al4

1/2 0 0 Al5
3/4 1/4 0 Al6
3/4 0 1/4 Al7
1/2 1/4 1/4 Al8

0 1/2 0 Al9
1/4 3/4 0 Al10
1/4 1/2 1/4 Al11
0 3/4 1/4 Al12

0 0 1/2 Al13
1/4 1/4 1/2 Al14
1/4 0 3/4 Al15
0 1/4 3/4 Al16

1/2 1/2 0 Al17
3/4 3/4 0 Al18
3/4 1/2 1/4 Al19
1/2 3/4 1/4 Al20

1/2 0 1/2 Al21
3/4 1/4 1/2 Al22
3/4 0 3/4 Al23
1/2 1/4 3/4 Al24

0 1/2 1/2 Al25
1/4 3/4 1/2 Al26
1/4 1/2 3/4 Al27
0 3/4 3/4 Al28

1/2 1/2 1/2 Al29
3/4 3/4 1/2 Al30
3/4 1/2 3/4 Al31
1/2 3/4 3/4 Al32
c------------------------------------------------------------


AkaiKKRでスーパーセル その1で書いたとおり、AkaiKKRでスーパーセルの計算を行うためには、それに適したパラメータをspecx.fに設定して再コンパイルする必要があります。計算するコンピュータのメモリが少ない場合、スワップ領域を使う必要があるかもしれません。今回はAkaiKKRとUbuntu 12.04 のスワップ領域で指定した下記のパラーメータをspecx.fに設定しました。

     & (natmmx=32, ncmpmx=32, msizmx=288, mxlmx=3, nk1x=500, nk3x=701,


結果と議論


面心立方構造のアルミニウムの点欠陥の形成エネルギーを計算するために、計算セルの中に32個のアルミニウム原子を置いた完全結晶の全エネルギー(E0)とスーパーセルからひとつの原子を点欠陥に置き換えたスーパーセルの全エネルギー(E1)の計算を行いました。得られた全エネルギーは、以下のようになりました。

E0 = -15482.440171667 (Ry)
E1 = -14998.335584355 (Ry)

よって点欠陥の生成エネルギーは
ΔEf = 0.2783319 (Ry) = 3.786906 (eV)
となりました。

この値はWang et al. (2004) (PDF)で報告されている 0.568 (eV) @ N=31, 0.511 (eV) @ N=108 とくらべて一桁程度の過大評価となってしまいました。したがって、AkaiKKRのCPA計算で点欠陥の生成エネルギーを過大評価してしまうのは、必ずしもCPAの問題ではないかもしれません。

関連エントリ




参考URL




付録


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


参考文献/使用機器




フィードバック



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

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


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


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


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

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ヶ月以上更新のないブログに表示されています。新しい記事を書くことで広告を消せます。