AkaiKKRの角運動量(方位量子数)のメモ

AkaiKKR(machikaneyama)を用いた第一原理計算でいつもいくつにするのだったか忘れがちな方位量子数の設定値に関してメモをまとめました。

d状態まで計算するなら mxl = 2 (入力ファイル)と mxlmx = 3 (specx.f)が必要で、f状態まで計算するなら mxl = 3, mxlmx = 4 が必要になります。

またspecx.fのmsizmxは、以下の値が必要と考えるのが簡単です。

msizmx ≧ natmmx * mxlmx2


軌道角運動量(方位量子数)


AkaiKKR(machikaneyama)の入力ファイル、及びspecx.fには軌道角運動量量子数(方位量子数)の設定に関するパラメータがあります。

  • mxl (入力ファイル)
  • mxlmx (specx.f)
  • msizmx (specx.f)

この方位量子数は l = 0, 1, 2, 3, ... と整数値を取り、順にs軌道、p軌道、d軌道、f軌道 ... と対応します(参考:電子配置 - Wikipedia)。

mxlとmxlmxの設定


AkaiKKRの入力ファイルの中では、どの軌道まで計算するかを入力ファイルの mxl で指定する必要があります。
更に、これを受けてspecx.fではmakeの前にあらかじめどの軌道まで計算できるかを mxlmx で指定しておく必要があります。この対応関係を以下に示します。


spdf
mxl0123
mxlmx1234
table.1: 計算する電子の軌道と設定すべきmxl,mxlmxの関係

ここで注意しなければならないのは、mxlmxのほうが常にmxlよりも1以上大きくなければならないと言う点です。このことはAkaiKKRのマニュアルや計算機マテリアルデザイン入門 (大阪大学新世紀レクチャー)では分かりにくいのですが、実を言うとspecx.fのなかにコメントとして記入されています。

note: mxlmx is l_max + 1 where l_max is the maximum angular momentum used in the calculation.

msizmxの設定


更にspecx.fのなかのmsizmxも、どの軌道まで計算するかに依存します。

AkaiKKRのマニュアルや計算機マテリアルデザイン入門 (大阪大学新世紀レクチャー)には以下の式が書かれています。

\sum_{i = 1}^{\mathrm{natm}}(l_{max,i}+1)^2

これは各原子位置において、異なる軌道まで計算するときのことを考えてシグマを使った和の形にしていますが、全て同じ軌道まで計算するのなら以下の様に簡単になります。

msizmx ≧ natmmx * mxlmx2

例えばhow to run a system of over 30 atoms?の例の様に単位格子の中の原子の数が30個でd軌道(l=2)まで計算する場合は msizemx ≧ 30 * (2 + 1)2 = 270 が必要になります。msizmxにいくつが必要になるのかの早見表を以下に示します。一列目が原子数です。

spdf
114916
2281832
33122748
44163664
55204580
66245496
772863112
883272128
993681144
10104090160
11114499176
121248108192
131352117208
141456126224
151560135240
161664144256
171768153272
181872162288
191976171304
202080180320
212184189336
222288198352
232392207368
242496216384
2525100225400
2626104234416
2727108243432
2828112252448
2929116261464
3030120270480
3131124279496
3232128288512
table.2: msizmxに必要な値の原子数・角運動量依存性


参考URL




参考文献/使用機器




フィードバック



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

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


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


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

tag: AkaiKKR machikaneyama specx.f 

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関数フィッティング

最新コメント
リンク

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