リンドハード理論のソースを表示

'''リンドハード理論'''<ref>{{cite journal |last=Lindhard |first=J. |date=1954 |title=On the properties of a gas of charged particles |url=http://gymarkiv.sdu.dk/MFM/kdvs/mfm%2020-29/mfm-28-8.pdf |journal=Danske Matematisk-fysiske Meddeleiser |publisher=Det Kongelige Danske Videnskabernes Selskab |volume=28 |issue=8 |pages=1–57 |access-date=2016-09-28 }}</ref><ref name=Ashcroft>N. W. Ashcroft and N. D. Mermin, ''Solid State Physics'' (Thomson Learning, Toronto, 1976)</ref>（リンドハードりろん）とは、固体中の電子による[[遮蔽]]効果を計算する手法である。リンドハード理論は、量子力学（第一原理摂動論）と[[乱雑位相近似]]に基づいている。

[[トーマス-フェルミ遮蔽]]は、より一般的なリンドハード公式の特別な場合として導出される。具体的に言うと波数ベクトルがフェルミ波数ベクトルよりはるかに小さいときのリンドハード公式、すなわち長距離極限がトーマス-フェルミ遮蔽である<ref name=Ashcroft/>。

この記事は[[CGSガウス単位系]]を用いることにする。

== 公式 ==
縦[[誘電関数]]のリンドハード公式は次式で与えられる。
::{|cellpadding="2" style="border:2px solid #ccccff"
| <math>\epsilon(q,\omega) = 1 - V_q \sum_k{\frac{f_{k-q}-f_k}{\hbar(\omega+i\delta)+E_{k-q}-E_k}}.</math>
|}
ここで<math>V_q</math>は<math>V_{eff}(q) - V_{ind}(q)</math>、<math>f_k</math>は熱力学的な平衡にある電子の[[フェルミ分布関数]]であるキャリア分布関数である。しかしこのリンドハード公式は非平衡分布関数においても有効である。

== リンドハード公式の解析 ==
リンドハード公式を理解するために、2次元と3次元における極限状態を考える。1次元の場合も、別のやり方で考える。

=== 3次元 ===
==== 長波長極限 ====
まず長波長極限(<math>q\to0</math>)を考える。
リンドハード公式の分母について、
: <math>E_{k-q} - E_k = \frac{\hbar^2}{2m}(k^2-2\vec{k}\cdot\vec{q}+q^2) - \frac{\hbar^2 k^2}{2m} \simeq -\frac{\hbar^2 \vec{k}\cdot\vec{q}}{m}</math>,
また分子について、
: <math>f_{k-q} - f_k = f_k - \vec{q}\cdot\nabla_k f_k + \cdots - f_k \simeq - \vec{q}\cdot\nabla_k f_k</math>.
これらをリンドハード公式に代入し、極限<math>\delta \to 0</math>をとると、
: <math>
\begin{alignat}{2}
\epsilon(0,\omega_0) & \simeq 1 + V_q \sum_{k,i}{ \frac{q_i \frac{\partial f_k}{\partial k_i}}{\hbar \omega_0 - \frac{\hbar^2 \vec{k}\cdot\vec{q}}{m}} }\\
& \simeq 1 + \frac{V_q}{\hbar \omega_0} \sum_{k,i}{q_i \frac{\partial f_k}{\partial k_i}}(1+\frac{\hbar \vec{k}\cdot\vec{q}}{m \omega_0})\\
& \simeq 1 + \frac{V_q}{\hbar \omega_0} \sum_{k,i}{q_i \frac{\partial f_k}{\partial k_i}}\frac{\hbar \vec{k}\cdot\vec{q}}{m \omega_0}\\
& = 1 - V_q \frac{q^2}{m \omega_0^2} \sum_k{f_k}\\
& = 1 - V_q \frac{q^2 N}{m \omega_0^2} \\
& = 1 - \frac{4 \pi e^2}{\epsilon q^2 L^3} \frac{q^2 N}{m \omega_0^2} \\
& = 1 - \frac{\omega_{pl}^2}{\omega_0^2}
\end{alignat}
</math>
ここで<math>E_k = \hbar \omega_k</math>、<math>V_q = \frac{4 \pi e^2}{\epsilon q^2 L^3}</math>、<math>\omega_{pl}^2 = \frac{4 \pi e^2 N}{\epsilon L^3 m}</math>を用いた。

（SI単位では因子<math>4\pi </math>を<math>1/\epsilon_{0}</math>に置き換わる）

この結果は古典的な誘電関数と同様である。

==== 静止限界 ====
次に静止極限(<math>\omega + i\delta \to 0</math>)を考える。
リンドハード公式は次のようになる。
: <math>\epsilon(q,0) = 1 - V_q \sum_k{\frac{f_{k-q}-f_k}{E_{k-q}-E_k}}</math>.
分母と分子に上記の式を代入すると、
: <math>\epsilon(q,0) = 1 - V_q \sum_{k,i}{\frac{-q_i \frac{\partial f}{\partial k_i} }{ -\frac{\hbar^2 \vec{k}\cdot\vec{q}}{m} }}
 = 1 - V_q \sum_{k,i}{\frac{q_i \frac{\partial f}{\partial k_i} }{\frac{\hbar^2 \vec{k}\cdot\vec{q}}{m} }}</math>
熱平衡におけるフェルミ-ディラックキャリア分布を仮定すると、
: <math>\sum_{i}{ q_i \frac{\partial f_k}{\partial k_i} } = -\sum_{i}{ q_i \frac{\partial f_k}{\partial \mu} \frac{\partial \epsilon_k}{\partial k_i} } = -\sum_{i}{ q_i k_i \frac{\hbar^2}{m} \frac{\partial f_k}{\partial \mu}}
</math>
ここで<math>\epsilon_k = \frac{\hbar^2 k^2}{2m}</math>、<math>\frac{\partial \epsilon_k}{\partial k_i} = \frac{\hbar^2 k_i}{m} </math>を用いた。

よって、
: <math>
\begin{alignat}{2}
\epsilon(q,0) & = 
1 + V_q \sum_{k,i}{\frac{ q_i k_i \frac{\hbar^2}{m} \frac{\partial f_k}{\partial \mu} }{\frac{\hbar^2 \vec{k}\cdot\vec{q}}{m} }} = 
1 + V_q\sum_k{\frac{\partial f_k}{\partial \mu}} = 
1 + \frac{4 \pi e^2}{\epsilon q^2} \frac{\partial}{\partial \mu} \frac{1}{L^3} \sum_k{f_k} \\
& = 1 + \frac{4 \pi e^2}{\epsilon q^2} \frac{\partial}{\partial \mu} \frac{N}{L^3} =
1 + \frac{4 \pi e^2}{\epsilon q^2} \frac{\partial n}{\partial \mu} \equiv
1 + \frac{\kappa^2}{q^2}.
\end{alignat}
</math>
ここで<math>\kappa</math>は<math>\kappa = \sqrt{ \frac{4\pi e^2}{\epsilon} \frac{\partial n}{\partial \mu} }</math>で定義される3次元遮蔽波数（3次元遮蔽長の逆数）である。

ここで3次元での静的に遮蔽されたクーロンポテンシャルは次のように与えられる。
: <math>V_s(q,\omega=0) \equiv \frac{V_q}{\epsilon(q,\omega=0)} = \frac {\frac{4 \pi e^2}{\epsilon q^2 L^3} }{ \frac{q^2 + \kappa^2}{q^2} }  = \frac{4 \pi e^2}{\epsilon L^3} \frac{1}{q^2 + \kappa^2}</math>
またこの結果のフーリエ変換は、
: <math>V_s(r) = \sum_q{ \frac{4\pi e^2}{\epsilon L^3 (q^2+\kappa^2)} e^{i \vec{q} \cdot \vec{r}} } = \frac{e^2}{\epsilon r} e^{-\kappa r}</math>
これは[[湯川ポテンシャル]]として知られる。
ここで、このフーリエ変換では基本的に「全ての」<math>\vec{q}</math>についての和をとり、「各」<math>\vec{q}</math>における小さな<math>|\vec{q}|</math>についての表現を使用するのは正しくないことに注意。

[[File:Screening.png|500px|thumb|3次元における静的に遮蔽されたポテンシャル（上の曲面）とクーロンポテンシャル（下の曲面）。]]

縮退したガス(T=0)において、[[フェルミエエネルギー]]は次式で与えられる。
: <math>E_f = \frac{\hbar^2}{2m}(3\pi^2 n)^{\frac{2}{3}} </math>,
よって密度は、
: <math>n = \frac{1}{3\pi^2}  \left(\frac{2m}{\hbar^2} E_f\right)^{\frac{3}{2}} </math>.
T=0では<math>E_f \equiv \mu</math>、よって<math>\frac{\partial n}{\partial \mu} = \frac{3}{2}\frac{n}{E_f}</math>。

これを上述の3次元遮蔽波数の式に代入すると、
::{|cellpadding="2" style="border:2px solid #ccccff"
| <math>\kappa = \sqrt{ \frac{4\pi e^2}{\epsilon} \frac{\partial n}{\partial \mu} } = \sqrt{ \frac{6\pi e^2 n}{\epsilon E_f} }</math>.
|}
これは3次元における[[トーマス-フェルミ遮蔽]]波数である。

なお、[[デバイ遮蔽]]は非縮退極限の場合を記述する。結果は<math>\kappa = \sqrt{ \frac{4\pi e^2 n \beta}{\epsilon} }</math>であり、3次元のデバイ遮蔽波数である。

=== 2次元 ===
==== 長波長極限 ====
まず長波長極限 (<math>q\to0</math>)を考える。リンドハード公式の分母について、
: <math>E_{k-q} - E_k = \frac{\hbar^2}{2m}(k^2-2\vec{k}\cdot\vec{q}+q^2) - \frac{\hbar^2 k^2}{2m} \simeq -\frac{\hbar^2 \vec{k}\cdot\vec{q}}{m}</math>,
また分子について、
: <math>f_{k-q} - f_k = f_k - \vec{q}\cdot\nabla_k f_k + \cdots - f_k \simeq - \vec{q}\cdot\nabla_k f_k</math>.
これらをリンドハード公式に代入し、極限<math>\delta \to 0</math>をとると、
: <math>
\begin{alignat}{2}
\epsilon(0,\omega) & \simeq 1 + V_q \sum_{k,i}{ \frac{q_i \frac{\partial f_k}{\partial k_i}}{\hbar \omega_0 - \frac{\hbar^2 \vec{k}\cdot\vec{q}}{m}} }\\
& \simeq 1 + \frac{V_q}{\hbar \omega_0} \sum_{k,i}{q_i \frac{\partial f_k}{\partial k_i}}(1+\frac{\hbar \vec{k}\cdot\vec{q}}{m \omega_0})\\
& \simeq 1 + \frac{V_q}{\hbar \omega_0} \sum_{k,i}{q_i \frac{\partial f_k}{\partial k_i}}\frac{\hbar \vec{k}\cdot\vec{q}}{m \omega_0}\\
& = 1 + \frac{V_q}{\hbar \omega_0} 2 \int d^2 k (\frac{L}{2 \pi})^2 \sum_{i,j}{q_i \frac{\partial f_k}{\partial k_i}}\frac{\hbar k_j q_j}{m \omega_0}\\
& = 1 + \frac{V_q L^2}{m \omega_0^2} 2 \int \frac{d^2 k}{(2 \pi)^2}  \sum_{i,j}{q_i q_j k_j \frac{\partial f_k}{\partial k_i}}\\
& = 1 + \frac{V_q L^2}{m \omega_0^2} \sum_{i,j}{ q_i q_j 2 \int \frac{d^2 k}{(2 \pi)^2} k_j \frac{\partial f_k}{\partial k_i}}\\
& = 1 - \frac{V_q L^2}{m \omega_0^2} \sum_{i,j}{ q_i q_j 2 \int \frac{d^2 k}{(2 \pi)^2} k_k \frac{\partial f_j}{\partial k_i}}\\
& = 1 - \frac{V_q L^2}{m \omega_0^2} \sum_{i,j}{ q_i q_j n \delta_{ij}}\\
& = 1 - \frac{2 \pi e^2}{\epsilon q L^2} \frac{L^2}{m \omega_0^2} q^2 n\\
& = 1 - \frac{\omega_{pl}^2(q)}{\omega_0^2},
\end{alignat}
</math>
ここで<math>E_k = \hbar \epsilon_k</math>、<math>V_q = \frac{2 \pi e^2}{\epsilon q L^2}</math>、<math>\omega_{pl}^2(q) = \frac{2 \pi e^2 n q}{\epsilon m}</math>を用いた。

==== 静止限界 ====
次に静止極限(<math>\omega + i\delta \to 0</math>)を考える。リンドハード公式は次のようになる。
: <math>\epsilon(q,0) = 1 - V_q \sum_k{\frac{f_{k-q}-f_k}{E_{k-q}-E_k}}</math>.
上述の式を分母と分子に代入すると、
: <math>\epsilon(q,0) = 1 - V_q \sum_{k,i}{\frac{-q_i \frac{\partial f}{\partial k_i} }{ -\frac{\hbar^2 \vec{k}\cdot\vec{q}}{m} }}
 = 1 - V_q \sum_{k,i}{\frac{q_i \frac{\partial f}{\partial k_i} }{\frac{\hbar^2 \vec{k}\cdot\vec{q}}{m} }}</math>
熱平衡でのフェルミ-ディラックキャリア分布を仮定すると、
: <math>\sum_{i}{ q_i \frac{\partial f_k}{\partial k_i} } = -\sum_{i}{ q_i \frac{\partial f_k}{\partial \mu} \frac{\partial \epsilon_k}{\partial k_i} } = -\sum_{i}{ q_i k_i \frac{\hbar^2}{m} \frac{\partial f_k}{\partial \mu}}
</math>
ここで<math>\epsilon_k = \frac{\hbar^2 k^2}{2m}</math>、<math>\frac{\partial \epsilon_k}{\partial k_i} = \frac{\hbar^2 k_i}{m} </math>を用いた。

よって、
: <math>
\begin{alignat}{2}
\epsilon(q,0) & = 
1 + V_q \sum_{k,i}{\frac{ q_i k_i \frac{\hbar^2}{m} \frac{\partial f_k}{\partial \mu} }{\frac{\hbar^2 \vec{k}\cdot\vec{q}}{m} }} = 
1 + V_q\sum_k{\frac{\partial f_k}{\partial \mu}} = 
1 + \frac{2 \pi e^2}{\epsilon q L^2} \frac{\partial}{\partial \mu} \sum_k{f_k} \\
& = 1 + \frac{2 \pi e^2}{\epsilon q} \frac{\partial}{\partial \mu} \frac{N}{L^2} =
1 + \frac{2 \pi e^2}{\epsilon q} \frac{\partial n}{\partial \mu} \equiv
1 + \frac{\kappa}{q}.
\end{alignat}
</math>

<math>\kappa</math>は<math>\kappa = \frac{2\pi e^2}{\epsilon} \frac{\partial n}{\partial \mu}</math>で定義される2次元での遮蔽波数（2次元での遮蔽長の逆数）である。

ここで2次元での静的に遮蔽されたクーロンポテンシャルは次式で与えられる。
: <math>V_s(q,\omega=0) \equiv \frac{V_q}{\epsilon(q,\omega=0)} = \frac{2 \pi e^2}{\epsilon q L^2} \frac{q}{q + \kappa} = \frac{2 \pi e^2}{\epsilon L^2} \frac{1}{q + \kappa}</math>.
2次元フェルミ気体の化学ポテンシャルは次式で与えられることが知られている。
: <math>\mu (n,T) = \frac{1}{\beta} \ln{(e^{\hbar^2 \beta \pi n/m}-1)}</math>,
また<math>\frac{\partial \mu}{\partial n} = \frac{\hbar^2 \pi}{m} \frac{1}{1-e^{-\hbar^2 \beta \pi n / m}}</math>である。

よって2次元での遮蔽波数は、
::{|cellpadding="2" style="border:2px solid #ccccff"
| <math>\kappa = \frac{2\pi e^2}{\epsilon}\frac{\partial n}{\partial \mu} = \frac{2\pi e^2}{\epsilon} \frac{m}{\hbar^2 \pi} (1-e^{-\hbar^2 \beta \pi n / m}) = \frac{2 m e^2}{ \hbar^2\epsilon} f_{k=0} .</math>
|}
これはnに依存しないことに注意。

=== 1次元 ===
ここでは次元を下げたいくつかの限られた場合を考える。次元を下げると、遮蔽効果は弱くなる。低次元では、一部の力線が遮蔽効果が無い物質を貫く。1次元の場合、ワイヤ軸に非常に近い力線にのみ遮蔽が影響を与えると考えられる。

==== 実験 ====
実際の実験では、単一フィラメントのような1次元の場合を扱っていたとしても、3次元バルクの遮蔽効果も考慮する必要がある。D. Davisは、フィラメントと同軸円筒に閉じ込められた電子ガスにトーマス-フェルミ遮蔽を適用した<ref>D. Davis ''[http://prola.aps.org/abstract/PRB/v7/i1/p129_1 Thomas-fermi screening in one dimension]'', Phys. Rev. B, 7(1), 129, (1973)</ref>。K<sub>2</sub>Pt(CN)<sub>4</sub>Cl<sub>0.32</sub>·2.6H<sub>2</sub>0では、フィラメントと同軸円筒との間の領域内のポテンシャルは<math>e^{-k_{eff}r}/r</math>で変化し、有効遮蔽長は金属[[白金]]の約10倍であることが分かっている。

== 関連項目 ==
* [[遮蔽]]

== 参考文献 ==
{{Reflist}}
*{{cite book |author1=Haug, Hartmut |author2=W. Koch, Stephan | title=Quantum Theory of the Optical and Electronic Properties of Semiconductors (4th ed.) | publisher=World Scientific Publishing Co. Pte. Ltd. | year=2004 | isbn=981-238-609-2}}

{{DEFAULTSORT:りんとはあとりろん}}
[[Category:物性物理学]]
[[Category:理論]]
[[Category:物理学のエポニム]]