tecktom.com

出てくるもの1:静的な遮蔽、Thomas–Fermi と Friedel 振動

2026年10月15日11分で読めますAI生成教科書『虎の巻』第5部 第1章 遮蔽と Fermi 液体
目次
  1. Thomas–Fermi の遮蔽
  2. Lindhard 関数の全体の形
  3. Friedel 振動
  4. 木下の層

AI 生成・木下が未勉強

ハミルトニアンの選び方で固体を考える 第5部 第1章 §5.1.3 ★★★ 金属に入れた電荷は、格子間隔より短い距離で完全に遮蔽される。ただし Fermi 面が鋭いので、遮蔽の雲の外側に波長 π/kF\pi/k_{\mathrm F} で振動して r−3r^{-3} で減る電荷の揺れ(Friedel 振動)が残る

Thomas–Fermi の遮蔽

電荷 QQ の点電荷を電子気体に入れる。ϕext(q)=Q/(ε0q2)\phi_{\mathrm{ext}}(\bm{q})=Q/(\varepsilon _0q^2) を式 (5.1.9) の εr\varepsilon _{\mathrm r} で割ると ϕ(q)=Q/(ε0(q2+qTF2))\phi(\bm{q})=Q/\bigl(\varepsilon _0(q^2+q_{\mathrm{TF}}^2)\bigr) で、実空間に戻すと

ϕ(r)=Q4πε0r e−qTFr(5.1.10)\phi(r)=\frac{Q}{4\pi\varepsilon _0r}\,\mathrm{e}^{-q_{\mathrm{TF}}r}\tag{5.1.10}

になる。裸の Coulomb ポテンシャルが、距離 1/qTF1/q_{\mathrm{TF}} を超えると指数関数的に消える。これが第4部第2章で不純物の散乱に使った遮蔽された Coulomb ポテンシャルで、そこで使った遮蔽の波数がこの qTFq_{\mathrm{TF}} である。遮蔽された相互作用の Fourier 成分 e2/(ε0(q2+qTF2))e^2/\bigl(\varepsilon _0(q^2+q_{\mathrm{TF}}^2)\bigr) は q→0q\to0 で有限になる。5.1.1節の Hartree–Fock で速さが発散したのは q→0q\to0 の 1/q21/q^2 のためだったので、遮蔽された相互作用で交換の自己エネルギーを計算し直すと、発散は消える。

遮蔽を担う電荷の分布は、Poisson 方程式 ρind=−ε0∇2ϕ−ρext\rho_{\mathrm{ind}}=-\varepsilon _0\nabla^2\phi-\rho_{\mathrm{ext}} から

ρind(r)=−QqTF24πre−qTFr\rho_{\mathrm{ind}}(r)=-\frac{Qq_{\mathrm{TF}}^2}{4\pi r}\mathrm{e}^{-q_{\mathrm{TF}}r}

で、全体を積分すると −Q-Q になる。外の電荷はちょうど打ち消される(完全遮蔽)。金属の内部に静電場が入り込めないことの、微視的な言い方である。

Lindhard 関数の全体の形

長波長の近似を外し、静的な Lindhard 関数を全ての qq で求める。式 (5.1.5) の2つの項を k→−k−q\bm{k}\to-\bm{k}-\bm{q} でまとめると、絶対零度ではχ0(q,0)=4V∑∣k∣<kF1εk−εk+q\chi_0(q,0)=\frac{4}{V}\sum_{\lvert \bm{k}\rvert <k_{\mathrm F}}\frac{1}{\varepsilon _{\bm{k}}-\varepsilon _{\bm{k}+\bm{q}}}(主値)になる。εk+q−εk=ℏ22m(2kqu+q2)\varepsilon _{\bm{k}+\bm{q}}-\varepsilon _{\bm{k}}=\frac{\hbar^2}{2m}(2kqu+q^2)(uu は角の余弦)の角の積分を先に行うと ∫−11du2kqu+q2=12kqln⁡∣q+2kq−2k∣\int_{-1}^{1}\frac{\mathrm{d} u}{2kqu+q^2}=\frac{1}{2kq}\ln\lvert \frac{q+2k}{q-2k}\rvert で、残りの kk の積分から

χ0(q,0)=−D(εF) L(q2kF)(5.1.11)\chi_0(q,0)=-D(\varepsilon _{\mathrm F})\,L\Bigl(\frac{q}{2k_{\mathrm F}}\Bigr)\tag{5.1.11}

を得る。LL は Hartree–Fock の式 (5.1.4) と同じ関数である。どちらも「Fermi 球の中の状態について 1/(波数の差)1/(\text{波数の差}) の重みで足す」という同じ和なので、同じ関数になる。

LL の3つの性質がそのまま物理になる。第1に、L(0)=1L(0)=1 で式 (5.1.9) の Thomas–Fermi に戻る。第2に、q≫2kFq\gg2k_{\mathrm F} では L≈13(2kF/q)2L\approx\frac13(2k_{\mathrm F}/q)^2 と小さくなり、εr→1\varepsilon _{\mathrm r}\to1 になる。電子は Fermi 波長より細かい凹凸のポテンシャルには集まりきれないので、短い距離の相互作用はほとんど遮蔽されない。第3に、q=2kFq=2k_{\mathrm F} で LL 自身は連続(値 1/21/2)だが、その傾きが対数的に発散する。q=2kFq=2k_{\mathrm F} は、Fermi 球の一方の端の電子を反対側の端へ移す波数で、この qq を境に「占有された状態から空いた状態へ、エネルギーをほとんど使わずに移れるかどうか」が切り替わる。Fermi 面が鋭いことが、qq の関数の角として残る。

図

図 5.1.1 銅の密度(qTF/kF=1.33q_{\mathrm{TF}}/k_{\mathrm F}=1.33)の電子気体で、遮蔽された相互作用と裸の相互作用の比 1/εr(q,0)1/\varepsilon _{\mathrm r}(q,0)。実線が Thomas–Fermi の近似 (5.1.9)、点線が Lindhard 関数 (5.1.11) による RPA の値。長い波長では両者が一致し、短い波長では Lindhard のほうが早く 11 に近づく(遮蔽が弱まる)。点線は q=2kFq=2k_{\mathrm F} で傾きが不連続に見えるが、実際には傾きが対数的に発散する角になっている。

Friedel 振動

q=2kFq=2k_{\mathrm F} の角は、実空間では遠方まで続く振動になる。最も簡単な場合として、原点に置いた短距離のポテンシャル u(r)=u0δ(3)(r)u(\bm{r})=u_0\delta^{(3)}(\bm{r}) に対する、相互作用しない電子の密度の変化 δn(r)=u0χ0(r)\delta n(r)=u_0\chi_0(r) を計算する。χ0(r)=∫d3q(2π)3eiq⋅rχ0(q,0)\chi_0(r)=\int\frac{\mathrm{d}^3q}{(2\pi)^3}\mathrm{e}^{\mathrm{i}\bm{q}\cdot\bm{r}}\chi_0(q,0) を式 (5.1.11) から求めると、

χ0(r)=−2mkF4π3ℏ2 sin⁡x−xcos⁡xx4,x=2kFr(5.1.12)\chi_0(r)=-\frac{2mk_{\mathrm F}^4}{\pi^3\hbar^2}\,\frac{\sin x-x\cos x}{x^4},\qquad x=2k_{\mathrm F}r\tag{5.1.12}

になる(数値の Fourier 変換と照合した)。遠方 x≫1x\gg1 では

δn(r)≈u0 2mkF4π3ℏ2 cos⁡(2kFr)(2kFr)3(5.1.13)\delta n(r)\approx u_0\,\frac{2mk_{\mathrm F}^4}{\pi^3\hbar^2}\,\frac{\cos(2k_{\mathrm F}r)}{(2k_{\mathrm F}r)^3}\tag{5.1.13}

で、密度は波長 π/kF\pi/k_{\mathrm F}(Fermi 波長の半分)で振動しながら r−3r^{-3} で減る。これを Friedel 振動と呼ぶ。Fourier 変換の遠方の振る舞いは、qq の関数のいちばん強い特異点で決まる。χ0(q,0)\chi_0(q,0) は滑らかでない点を q=2kFq=2k_{\mathrm F} にだけ持つので、遠方の振動の波数は 2kF2k_{\mathrm F} になる。Thomas–Fermi の指数関数的な遮蔽 (5.1.10) には振動が無いが、それは LL を L(0)=1L(0)=1 で置き換えて角を捨てたからである。遮蔽を自己無撞着に入れても、εr(q,0)\varepsilon _{\mathrm r}(q,0) の滑らかでない点は同じ q=2kFq=2k_{\mathrm F} にあるので、振動の波長と r−3r^{-3} の減り方は変わらない。

式 (5.1.12) は、スピンを持つ磁性の不純物が伝導電子のスピンを偏らせる場合にもそのまま使え、離れた2つの磁性の不純物の間に符号が振動する交換相互作用を作る(RKKY 相互作用。第5部第2章)。また、イオンの間の相互作用も電子に遮蔽されるので、フォノンの分散 ω(q)\omega(\bm{q}) も q=2kFq=2k_{\mathrm F} で傾きが変わる(Kohn 異常)。

図

図 5.1.2 Friedel 振動。式 (5.1.12) に x3=(2kFr)3x^3=(2k_{\mathrm F}r)^3 を掛けた (sin⁡x−xcos⁡x)/x(\sin x-x\cos x)/x を kFrk_{\mathrm F}r の関数として描いた(厳密な式なので点線)。r−3r^{-3} の減衰を掛け戻したので振幅は一定に見え、波数 2kF2k_{\mathrm F}(kFrk_{\mathrm F}r で周期 π\pi)の振動が遠方まで続くことが分かる。

例題 5.1.3 銅の中の点電荷を遮蔽する雲の大きさ ★★

オリジナル銅の伝導電子(n=8.47×1028 m−3n=8.47\times10^{28}\,\mathrm{m^{-3}}、εF=7.03 eV\varepsilon _{\mathrm F}=7.03\,\mathrm{eV}、kF=1.36×1010 m−1k_{\mathrm F}=1.36\times10^{10}\,\mathrm{m^{-1}})に、電荷 QQ の不純物を入れる。

  1. (1) qTFq_{\mathrm{TF}} と遮蔽長 1/qTF1/q_{\mathrm{TF}} を求めよ。qTF/(2kF)q_{\mathrm{TF}}/(2k_{\mathrm F}) はいくらか。
  2. (2) 半径 RR の球の中に入る誘起電荷の割合を qTFRq_{\mathrm{TF}}R の式で求め、90%90\% が入る RR を求めよ。銅の最近接原子の距離 0.256 nm0.256\,\mathrm{nm} と比べよ。
  3. (3) 遮蔽の雲が不純物の位置に作る静電ポテンシャル ϕind(0)\phi_{\mathrm{ind}}(0) を求めよ。
解答

(1) D(εF)=3n/(2εF)D(\varepsilon _{\mathrm F})=3n/(2\varepsilon _{\mathrm F}) から qTF2=3ne22ε0εF=3.27×1020 m−2q_{\mathrm{TF}}^2=\frac{3ne^2}{2\varepsilon _0\varepsilon _{\mathrm F}}=3.27\times10^{20}\,\mathrm{m^{-2}}、qTF=1.81×1010 m−1q_{\mathrm{TF}}=1.81\times10^{10}\,\mathrm{m^{-1}}、1/qTF=0.0552 nm1/q_{\mathrm{TF}}=0.0552\,\mathrm{nm}。qTF/(2kF)=0.665q_{\mathrm{TF}}/(2k_{\mathrm F})=0.665 で、Thomas–Fermi が当てはまる q≪2kFq\ll2k_{\mathrm F} の範囲の端にある。遮蔽長が Fermi 波長と同じ大きさなので、遮蔽の雲の形の細部は Lindhard 関数で計算し直す必要がある(図 5.1.1)。

(2) ∫0R4πr2ρind dr=−Q[1−(1+qTFR)e−qTFR]\int_0^R4\pi r^2\rho_{\mathrm{ind}}\,\mathrm{d} r=-Q\bigl[1-(1+q_{\mathrm{TF}}R)\mathrm{e}^{-q_{\mathrm{TF}}R}\bigr]。1−(1+t)e−t=0.91-(1+t)\mathrm{e}^{-t}=0.9 を解くと t=3.89t=3.89 で、R=3.89/qTF=0.215 nmR=3.89/q_{\mathrm{TF}}=0.215\,\mathrm{nm}。最近接距離 0.256 nm0.256\,\mathrm{nm}(1/qTF1/q_{\mathrm{TF}} の 4.64.6 倍)より小さい。不純物の電荷は、隣の原子の位置ではほとんど打ち消されている。

(3) ϕ(r)=Q4πε0re−qTFr\phi(r)=\frac{Q}{4\pi\varepsilon _0r}\mathrm{e}^{-q_{\mathrm{TF}}r} から裸の Q4πε0r\frac{Q}{4\pi\varepsilon _0r} を引いて r→0r\to0 とすると ϕind(0)=−QqTF4πε0\phi_{\mathrm{ind}}(0)=-\frac{Qq_{\mathrm{TF}}}{4\pi\varepsilon _0}。雲の電荷 −Q-Q が、半径 1/qTF1/q_{\mathrm{TF}} の殻に集まっているのと同じポテンシャルを作る。Q=eQ=e なら −eqTF/(4πε0)=−26 V-eq_{\mathrm{TF}}/(4\pi\varepsilon _0)=-26\,\mathrm V で、電子のエネルギーにして 26 eV26\,\mathrm{eV} と大きい。この値は Thomas–Fermi の雲を点の近くまで信じた見積もりで、実際には Lindhard 関数の q≫2kFq\gg2k_{\mathrm F} での弱い遮蔽のために点の近くの遮蔽はもっと弱い。

ほかの巻では

  • 近似の数理 第6章:Thomas–Fermi 理論を多体の極限として正当化した所で、遮蔽の波数 qTFq_{\mathrm{TF}} の数学の側が読める

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『虎の巻』の記事一覧へ