AkaiKKRでルチル構造SnO2 その1

AkaiKKR(machikaneyama)を用いてルチル構造のSnO2のバンド構造を計算しました。

001_20160710125630051.png
Fig.1: ルチル構造SnO2のバンド構造



ルチル構造半導体


Fig.2に示したのがルチル型の結晶構造です。

002_20160710125630226.jpg
Fig.2: ルチル型結晶構造


今回は、ルチル構造のSnO2の状態密度とバンド構造を計算します。
ルチル構造は比較的格子間領域の広い結晶構造なので、マフィンティン半径の比や追加のempty sphereなど考えるべきものは色々あるかと思いますが、今回は単純に下記のような入力ファイルを用いることにします。

c------------------------------------------------------------
go data/SnO2
c------------------------------------------------------------
c brvtyp a c/a b/a alpha beta gamma
st 8.95, 0.673, , , , ,
c------------------------------------------------------------
c edelt ewidth reltyp sdftyp magtyp record
0.001 2.1 sra vwn nmag 2nd
c------------------------------------------------------------
c outtyp bzqlty maxitr pmix
update 4 200 0.023
c------------------------------------------------------------
c ntyp
2
c------------------------------------------------------------
c type ncmp rmt field mxl anclr conc
Sn 1 1 0.0 2
50 100
O 1 1 0.0 2
8 100
c------------------------------------------------------------
c natm
6
c------------------------------------------------------------
c atmicx atmtyp
0.0a 0.0b 0.0c Sn
1/2a 1/2b 1/2c Sn
0.307a 0.307b 0.0c O
0.693a 0.693b 0.0c O
0.193a 0.807b 1/2c O
0.807a 0.193b 1/2c O
c------------------------------------------------------------


ポスト処理として、状態密度とバンド構造の計算も行いました。ブロッホスペクトル関数の計算(spc計算)では、AkaiKKRのspc計算の入出力フォーマットの通りに22 May 2015のバージョンと同じ入出力フォーマットにしました。k点の座標は、AkaiKKRでSrTiO3ペロフスカイトのScilabスクリプトを利用しました。

結果


Fig.3に示したのがSnO2の状態密度の計算結果です。

003_20160710125629bb2.png
Fig.3: SnO2の状態密度


今回の計算では ewidth を上手に選ぶ必要がありました(参考: AkaiKKRのewidth その1)。今回の例では、-1.8 Ry 付近にあるコアの状態も ewidth の範囲に含むことで収束させました。

Fig.4に示したのがSnO2のバンド構造です。

004_201607101256294af.png
Fig.4: SnO2のバンド構造


これだと価電子帯や伝導帯の構造が分かりにくいので、さらに ewidth を小さくして計算した結果が、冒頭のFig.1です。
この結果はSvane and Antoncik (1987)による計算結果と調和的です。

関連エントリ




参考URL




付録


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


参考文献/使用機器




フィードバック



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

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


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


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


tag: AkaiKKR machikaneyama KKR ルチル構造 

comment

Secret

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

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

最新コメント
リンク

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