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カウンター
カテゴリ
ユーザータグ

LTspiceAkaiKKRmachikaneyamaScilabKKRPSoC強磁性CPAPICOPアンプecalj状態密度常微分方程式モンテカルロ解析odeトランジスタインターフェースDOSPDS5022スイッチング回路定電流半導体分散関係シェルスクリプト乱数レベルシフトHP6632A可変抵抗温度解析トランジスタ技術ブレッドボードR6452AI2C確率論セミナー数値積分反強磁性バンドギャップ熱設計非線形方程式ソルバ絶縁バンド構造偏微分方程式三端子レギュレータフォトカプラカオスマフィンティン半径ISO-I2CGW近似LM358A/DコンバータシュミットトリガLEDUSB数値微分サーボアナログスイッチ補間発振回路カレントミラー直流動作点解析TL43174HC4053PC817C単振り子FFTVESTA開発環境bzqlty電子負荷量子力学基本並進ベクトルパラメトリック解析標準ロジックチョッパアンプBSchLDAアセンブラブラべ格子2ちゃんねるイジング模型PWscf状態方程式仮想結晶近似キュリー温度Quantum_ESPRESSO熱伝導VCAスイッチト・キャパシタewidth最適化QSGWTLP621GGASMPMaxima失敗談位相図六方最密充填構造繰り返しスピン軌道相互作用相対論ランダムウォークFETgfortranコバルトスレーターポーリング曲線ラプラス方程式抵抗cygwin不規則合金格子比熱熱力学マントル条件分岐MCU井戸型ポテンシャルダイヤモンドQNAPUPS固有値問題シュレディンガー方程式自動計測ガイガー管詰め回路OpenMPTLP521ハーフメタルLM555ubuntufsolveブラウン運動平均場近似NE555ZnOTLP552QuantumESPRESSOxcrysdenCIF最小値最大値awkフェルミ面テスタ第一原理計算Ubuntu差し込みグラフFXA-7020ZR三角波過渡解析Writer509データロガースーパーセル起電力CK1026AACircuitMAS830LフィルタMBEP-10PGAトランスナイキスト線図ノコギリ波負帰還安定性EAGLEOPA2277PIC16F785CapSenseLMC6622SC1815入出力固定スピンモーメントFSMTeX結晶磁気異方性全エネルギーc/a合金multiplotgnuplot非線型方程式ソルバL10構造正規分布等高線ジバニャン方程式初期値interp1fcc面心立方構造ウィグナーザイツ胞半金属デバイ模型磁気モーメント電荷密度重積分SIC不純物問題ゼーベック係数cif2cellPWgui擬ポテンシャル二相共存ウルツ鉱構造edeltquantumESPRESSOフォノンリジッドバンド模型スワップ領域BaO岩塩構造ルチル構造ヒストグラム確率論マテリアルデザインフラクタルマンデルブロ集合キーボードRealforceクーロン散乱三次元疎行列縮退化学反応関数フィッティング最小二乗法Excel直流解析PCTS-110TS-112日本語パラメータ・モデル等価回路モデル文字列状態図陰解法熱拡散方程式HiLAPW両対数グラフCrank-Nicolson法連立一次方程式specx.fifort境界条件片対数グラフグラフの分割円周率ヒストグラム不規則局所モーメントGimpシンボル軸ラベル凡例線種トラックボール

最新コメント
リンク

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