スポンサーサイト

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


AkaiKKRのewidth その2

AkaiKKRのewidth その1では、AkaiKKR(machikayenama)入力ファイルの中の ewidth, edelt と実際に計算される複素エネルギーの経路の関係をまとめました。
今回は ewidth が不適切な値に設定されているときにどのような問題が起こり得るかについて書きます。

ewidth が適切でない場合、計算が収束しない場合があるだけでなく、計算が間違った結果に収束してしまう場合も存在するので注意が必要です。


ewidthが小さすぎる場合


まず ewidth が小さすぎる場合について見てみます。多くの場合、go計算で収束しないだけですが、稀に正しくない結果に収束する場合があるようなので注意が必要です。

Fig.1に示したのは正しいダイヤモンドの状態密度で、下記の入力ファイルで計算したものです。更にこの入力ファイルを ewidth=1.0 Ryに変更した結果がFig.2です。

c------------------------------------------------------------
go data/diamond
c------------------------------------------------------------
c brvtyp a c/a b/a alpha beta gamma
fcc 6.74 , , , , , ,
c------------------------------------------------------------
c edelt ewidth reltyp sdftyp magtyp record
0.001 2.0 sra vwnasa nmag 2nd
c------------------------------------------------------------
c outtyp bzqlty maxitr pmix
update 8 200 0.02
c------------------------------------------------------------
c ntyp
2
c------------------------------------------------------------
c type ncmp rmt field mxl anclr conc
C 1 1 0.0 2 6 100
Vc 1 1 0.0 2 0 100
c------------------------------------------------------------
c natm
4
c------------------------------------------------------------
c atmicx atmtyp
0.00000 0.00000 0.00000 C
0.25000 0.25000 0.25000 C
0.75000 0.75000 0.75000 Vc
0.50000 0.50000 0.50000 Vc
c------------------------------------------------------------

c------------------------------------------------------------
dos data/diamond
c------------------------------------------------------------
c brvtyp a c/a b/a alpha beta gamma
fcc 6.74 , , , , , ,
c------------------------------------------------------------
c edelt ewidth reltyp sdftyp magtyp record
0.001 2.67 sra vwnasa nmag 2nd
c------------------------------------------------------------
c outtyp bzqlty maxitr pmix
update 20 200 0.02
c------------------------------------------------------------
c ntyp
2
c------------------------------------------------------------
c type ncmp rmt field mxl anclr conc
C 1 1 0.0 2 6 100
Vc 1 1 0.0 2 0 100
c------------------------------------------------------------
c natm
4
c------------------------------------------------------------
c atmicx atmtyp
0.00000 0.00000 0.00000 C
0.25000 0.25000 0.25000 C
0.75000 0.75000 0.75000 Vc
0.50000 0.50000 0.50000 Vc
c------------------------------------------------------------


Diamond.png
Fig.1: ダイヤモンドの正しい状態密度(ewidth = 2.0 Ry)


Diamond-ewidth.png
Fig.2: ダイヤモンドの誤った状態密度(ewidth = 1.0 Ry)


状態密度の形状を比べてみると明らかなように、これら二つの計算は、それぞれ異なった解に収束しています。あらかじめ正しい状態密度の形を知らずに間違った計算例だけを見せられても、結果が正しいのか間違っているのかは判断できないでしょう。従ってAkaiKKRで計算を行う場合は、あらかじめ先行研究の結果を調べておく、あるいは、別の計算コード(HiLAPWとかecaljとか)と同じ結果が得られるかを確認する必要があると思います(参考: AkaiKKRとecaljでCuGaTe2 その1その2)。

ewidth がコア近傍を横切る場合


Fig.3では ewidth が valence state を完全に含み、core state を含まないようにしていることが分かります。このようにするのが理想的なのですが、実際には core state が valence state に近くて ewidth の底がちょうどその間に来るようにすることが難しい場合もあります。そういうときには ewidth が core state を含むようにしても問題なく計算できます。(core state を valence state として扱うようになります。)

akaikkr-handson-2-140815051743-phpapp02.png

Fig.3: 状態密度と複素数エネルギーメッシュの関係。第5回CCMSハンズオン(ソフトウェア講習会): AkaiKKRチュートリアル 2. AkaiKKRの実習より


なお ewidth の底が core state に近い場合、下記のメッセージが出力されます。

***msg in cstate...corelevel near ebtm found for nclr=原子番号


これはメッセージなので、必ずしも問題が起きていることを意味しませんが、もしも上記メッセージが出ていて、かつ、収束しないときや、収束していても結果が怪しいようなら ewidth の値を変更してみると問題が解決するかもしれません。

関連エントリ




参考URL




参考文献/使用機器




フィードバック



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

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


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


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


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