多重散乱理論

出典: フリー百科事典『ウィキペディア（Wikipedia）』

散乱理論では単独のポテンシャルにおける電子（散乱するものは電子以外にも光や他の粒子など様々なものが存在）の散乱を扱ったが、現実の散乱は、多数のポテンシャル下でかつ散乱される対象も多数存在する。また一つの電子に限っても、散乱は一回限りでなく複数回散乱される。このような多重な散乱を扱う理論が多重散乱理論(Multiple scattering theory)である。

[編集] 格子上に配置したランダムなポテンシャル下での電子

多重散乱理論には扱う対象により様々なものが考えられるが、以下に一つの例として並進対称に配置した格子系において、各格子（サイト）上にポテンシャルがランダム（非周期的）に配置した場合を考える。以下、散乱されるのは電子としておく。

サイトnにあるポテンシャルをV_n、自由電子（または無摂動）のハミルトニアンをH₀として、系を記述するハミルトニアンHを、

$H =\, H_0 + \sum_n V_n$

とする。次にこれを以下のように変形する。

$H = [H_0 + \tilde{V}(z)] + \sum_n [V_n - \tilde{V}_n(z)] = \tilde{H}(z) + v(z)$

ここで、

$\tilde{V}(z) = \sum_n \tilde{V}_n (z)$

であり、 $\tilde{V}_n(z)$ は任意の周期ポテンシャル。つまりポテンシャルV_nを周期的部分 $\tilde{V}$ と非周期的部分 $\,v$ とに分けた訳である。zは複素エネルギー。上式で、v(z)は次のようv_n(z)の和になっている。

$v(z) = \sum_n [V_n - \tilde{V}_n (z)] = \sum_n v_n(z)$

更に、この系におけるグリーン関数をG(z)とすると、G(z)は、

$G(z) = {1 \over {z - \tilde{H} - v} } = {1 \over {(z - \tilde{H}) (1 - {v \over {z - \tilde{H}} }) }}$

であり、

$\tilde{G} = {1 \over {z - \tilde{H} }}$

とし、非周期ポテンシャル部分vに関して展開すると、

$G(z) = \tilde{G} \{1 + v \tilde{G} + v \tilde{G} v \tilde{G} + \cdot \cdot \cdot \} = \tilde{G} + \tilde{G} T \tilde{G}$

$T = v + v \tilde{G} v + v \tilde{G} v \tilde{G} v + \cdot \cdot \cdot$

となる。Tを総散乱行列と言う。総散乱行列Tをサイトの和の形で表すと、

$T = \sum_n v_n + \sum_n v_n \tilde{G} \sum_m v_m + \sum_n \tilde{G} \sum_m v_m \tilde{G} \sum_p v_p \cdot \cdot \cdot$

となる。サイトnのポテンシャルv_nのみを考え、散乱理論の場合と同じ要領でt行列が定義できる。

$t_n = v_n \{1 + \tilde{G} v_n + \tilde{G} v_n \tilde{G} v_n \cdot \cdot \cdot \} = v_n {1 \over {1 - v_n \tilde{G} }} = v_n [1 - v_n \tilde{G}]^{-1}$

加えて、

$t_n = v_n + v_n \tilde{G} \{ v_n + v_n \tilde{G} v_n + v_n \tilde{G} v_n \tilde{G} v_n \cdot \cdot \cdot \} = v_n + v_n \tilde{G} t_n$

である。総散乱行列TはサイトnでのT_nの和、

$T = \, \sum_n T_n$

と表現でき、各T_nは、

$\begin{matrix} T_n & = & v_n + v_n \tilde{G} \sum_m v_m + v_n \tilde{G} \sum_m v_m \tilde{G} \sum_p v_p + \cdot \cdot \cdot \\ \ & = & v_n \{ 1 + \tilde{G} [ \sum_m v_m + \sum_m v_m \tilde{G} \sum_p v_p + \cdot \cdot \cdot ] \} \\ \ & = & v_n [1 + \tilde{G} \sum_m T_m ] \end{matrix}$

更に、

$T_n = v_n + v_n \tilde{G} T_n + v_n \tilde{G} \sum_{m \ne n} T_m = t_n [1 + \tilde{G} \sum_{m \ne n} T_m ]$

である。ここで、

$(1 - v_n \tilde{G})T_n = v_n + v_n \tilde{G} \sum_{m \ne n} T_m, \quad \quad t_n = v_n [1 - v_n \tilde{G}]^{-1}$

よりt_nが出てくる。以上から総散乱行列Tは、t行列により次のように表される。

$T = \sum_n t_n + \sum_n t_n \tilde{G} \sum_{m \ne n} t_m + \sum_n t_n \tilde{G} \sum_{m \ne n} t_m \tilde{G} \sum_{p \ne m} t_p + \cdot \cdot \cdot$

[編集] 形式解の提示

ここでポテンシャルが全て同じであると考える。そして総散乱行列Tを次のように分解する。

$T = \, \sum_{n,n'}T_{nn'}$

分解されたT_nn'は、

$T_{nn'} = t_n \delta_{nn'} + t_n G_0 \sum_{m \ne n} T_{mn'}$

となる。G₀は自由電子のグリーン関数とする（ $\tilde{G} \to G_0$ ）。これにより、厳密な形式解を得ることができる。T_nn'は更に、

$\begin{matrix} T_{nn'} & = & t_n \delta_{nn'} + t_n G_0 \sum_{m \ne n} (t_m \delta_{mn'} + t_m G_0 \sum_{p \ne m} T_{pn'}) \\ \ & = & t_n \delta_{nn'} + t_n G_0 t_{n'} + t_n G_0 \sum_{m \ne n} t_m G_0 t_{n'} + \cdot \cdot \cdot \end{matrix}$

となる。T_nn'はサイトnから始まって、サイトn'で終わる全ての散乱過程を記述していることとなる。一方T_nは、

$T_n = \, \sum_{n'} T_{nn'}$

であり、これはサイトnは考慮されるが、終点としてのサイトn'を考えていない。そして、T_nn'の形式解は（但し、ここでr→kへのフーリエ変換及び、角運動量表示を導入している）、

$T_{nn'}^{LL'}(\kappa) = \tau_{n}^{l}(\kappa) [ \delta_{nn'}^{LL'} + \sum_{n_1,L_1} B_{nn_1}^{LL_1}(\kappa) \cdot T_{{n_1}n'}^{{L_1}L'}(\kappa)]$

$T_{nn'} (\mathbf{k}, \mathbf{k'}) = (4 \pi)^2 \sum_{L,L'} Y_L (\mathbf{k}) T_{nn'}^{LL'} (k, k')Y_{L'} (\mathbf{k'})$ 　：　角運動量表示

となる（形式解導出の詳細は省略）。 $B_{nn_1}^{LL_1}$ は構造定数と言われるもので、結晶格子の種類にのみ依存する定数である。 $\kappa = k = \sqrt{E}$ であり。L,L',lなどは軌道角運動量に関しての指標である。τ_n(κ)はt行列t_nに相当する。ここで構造定数は具体的には、

$B_{nn_1}^{LL_1}(\kappa) = - [(4 \pi i \kappa) \sum_{L_1} i^{l-l'-l_1} C_{LL'L_1} Y_{L_1}(\mathbf{R}_n - \mathbf{R}_{n'}) h_l^{+}(\kappa |\mathbf{R}_n - \mathbf{R}_{n'}|) ](1 - \delta_{nn'})$

$C_{LL'L_1} = \int Y_L(\mathbf{q}) Y_{L'} (\mathbf{q}) Y_{L_1}(\mathbf{q}) d \Omega_q$

となる。 $h_l^{+}$ は球ハンケル関数、Y_Lは球面調和関数である。尚、形式解は次のようにも表される。

$T_{nn'}^{LL'}(\kappa) = \{ [\tau^{-1}(\kappa) - B(\kappa)]^{-1} \}_{nn'}^{LL'}$

この形式解から、状態密度をD(E)の表式を得ることができる。この時、上式左辺を $T_{nn'}^{LL'}(\kappa) = T(\kappa)$ と略して表示。

$\begin{matrix} D(E) - D_0(E) & = & {2 \over {N \pi} } \mathrm{Im Tr} {d \over {dE} } \ln [T(\kappa)] \\ \ & = & - {2 \over {N \pi} } \mathrm{Im Tr} {d \over {dE} } \ln [\tau^{-1} - B(\kappa)] \\ \ & = & - {2 \over {N \pi} } \mathrm{Im Tr} [T(\kappa) ({ d \tau^{-1} \over {dE} } - { d B(\kappa) \over {dE} } ) ] \\ \ & = & - {2 \over {N \pi} } \mathrm{Im} \sum_{n,L} \sum_{n_1, L_1} T_{nn_1}^{LL_1} [ \delta_{{n_1}n}^{{L_1}L} {d (\tau_n^l (\kappa) )^{-1} \over {dE} } - {d \over {dE} } B_{{n_1}n}^{{L_1}L} (\kappa) ] \end{matrix}$

ここで、係数2はスピンの縮重度、Nは全サイト数、Imは虚数部分、Trはトレース（跡）を取ることを意味する。D₀(E)は自由電子の状態密度。

[編集] D₀(E)の起源

グリーン関数から状態密度を求める式は（エネルギーは全てEとする）、

$D(E) = - {1 \over {\pi}} \mathrm{Im} G(E) = - {1 \over {\pi}} \mathrm{Im} \{ \tilde{G}(E) + \tilde{G}(E) T(E) \tilde{G} (E) \}$

であり（スピン縮重度などの係数は省略）、ここで $\tilde{G}(E) \to G_0(E)$ とすると、

$\begin{matrix} D(E) = - {1 \over {\pi}} \mathrm{Im} G(E) & = & - {1 \over {\pi}} \mathrm{Im} \{ G_0(E) + G_0(E) T(E) G_0 (E) \} \\ \ & = & D_0 (E) - {1 \over {\pi}} \mathrm{Im} G_0(E) T(E) G_0(E) \end{matrix}$

となりD₀(E)を移項すると、D(E) - D₀(E)が出てくる。

以上は、ポテンシャルを全て同一とみなしたが、最初の前提であるポテンシャルがランダムである場合、その扱いは難しくなる。ランダムさ（乱れ）には構造的な乱れ、配置の乱れなど多様な状況を考えることができるが、ここでは先にあるように原子の配置のみが乱れた系である置換型の不規則二元合金を考えるのが比較的扱いが楽である。このランダムな問題を解くものとして、平均化によってランダムさを一様なものとして扱うアプローチがある。これに関係する近似手法として単サイト近似、平均場近似（有効媒質近似）がある。多重散乱理論を出発点として、このようなランダムな系を扱うバンド計算手法として、ATAやCPAがある。ランダムでない通常の周期的な系を、多重散乱理論を利用して解くバンド計算手法にKKR法がある。

"http://ja.wikipedia.org../../../%E5%A4%9A/%E9%87%8D/%E6%95%A3/%E5%A4%9A%E9%87%8D%E6%95%A3%E4%B9%B1%E7%90%86%E8%AB%96.html" より作成

カテゴリ: 量子力学

多重散乱理論

出典: フリー百科事典『ウィキペディア（Wikipedia）』

[編集] 格子上に配置したランダムなポテンシャル下での電子

[編集] 形式解の提示

[編集] D₀(E)の起源

Views

ナビゲーション

ヘルプ

検索

多重散乱理論

出典: フリー百科事典『ウィキペディア（Wikipedia）』

[編集] 格子上に配置したランダムなポテンシャル下での電子

[編集] 形式解の提示

[編集] D0(E)の起源

Views

ナビゲーション

ヘルプ

検索

[編集] D₀(E)の起源