スピントロニクス理論の基礎/8-10

(782d) 更新


[前の章へ] <<<<  スピントロニクス理論の基礎(目次) >>>> [次の章へ]

8-10 不純物散乱のもとでの Green 関数

(8.113) の Fourier 変換表示を求める

(8.116)

&v_i(\bm q)\equiv \frac{1}{V}\int d^3re^{i\bm q\cdot\bm r}v_i(\bm r)\\ &=\frac{1}{V}\sum_k^{N_i}v_ka^3\int d^3re^{i\bm q\cdot\bm r}[\delta(\bm r-\bm R_k)-1/V]\\ &=\frac{a^3}{V}\sum_k^{N_i}v_k(e^{i\bm q\cdot\bm R_k}-\delta_{\bm q,0})

を不純物散乱ポテンシャルの Fourier 変換と定義する。

この定義は (8.74) 等に与えられた c(\bm r,t) の Fourier 変換とは V の配置や i の符号などが異なるので注意が必要。

これに対応するハミルトニアンは、

(8.116)

&V_i=\int d^3rv_i(\bm r)c^\dagger(\bm r)c(\bm r)\\ &=\int d^3r \left(\sum_{\bm q}e^{-i\bm q\cdot\bm r}v_i(\bm q)\right) \left(\frac{1}{\sqrt V}\sum_{\bm k}e^{-i\bm k\cdot\bm r}c_{\bm k}^\dagger\right) \left(\frac{1}{\sqrt V}\sum_{\bm k'}e^{i\bm k'\cdot\bm r}c_{\bm k'}\right)\\ &=\sum_{\bm q,\bm k,\bm k'}v_i(\bm q)c_{\bm k}^\dagger c_{\bm k'}\frac{1}{V}\int d^3re^{i(-\bm q-\bm k+\bm k')\cdot\bm r}\\ &=\sum_{\bm q,\bm k,\bm k'}v_i(\bm q)c_{\bm k}^\dagger c_{\bm k'}\delta_{\bm q+\bm k,\bm k'}\\ &=\sum_{\bm q,\bm k}v_i(\bm q)c^\dagger(\bm k)c(\bm k+\bm q)

となる。

今考えているようなポテンシャルでは電子の振動数が変わらないため、単一の \omega によって表せる。これは以下で自然に導出される。また実際には「不純物平均」により k も変わらないという結果になる。

g_{\bm k,\bm k',\omega,\omega'}= 2\pi \delta(\omega-\omega')g_{\bm k,\bm k',\omega}

これらは 8-7 で見た g_0 と同じ性質である。

g_{0\bm k,\bm k',\omega,\omega'}= 2\pi\delta(\omega-\omega')\delta_{\bm k,\bm k'}g_{0\bm k,\omega}

→ 単一の \omega \bm k で表せるという結果は実時間・実空間で見れば 系の運動が絶対的な時刻や座標に依存せず、相対値のみで記述できることに対応する。

(8.117)

係数を正しく付けるために次元について確認:

  • g_{\bm k,\bm k',\omega,\omega'} は (8.76) によれば Gd^3r で無次元、 d^3r'/V で無次元、 dtdt'/\hbar が残って (時間/エネルギー) の次元を持つ
  • g_{\bm k,\bm k',\omega} \delta(\omega-\omega') が外に出るので、(1/エネルギー) の次元を持つ
  • v_i(\bm q) は (8.116) によれば v_i(\bm r) と同じ次元で、(8.106) によればこれはポテンシャルではなく (エネルギー) の次元を持つ

(8.47)〜(8.77) および (8.116) に与えられたフーリエ係数の定義は係数の付け方などが数学の授業とは異なるので、 慎重に計算していかないと正しい答えが得られない。実際に計算を進めるとかなり大変。

&g_{\bm k,\bm k',\omega,\omega'}^r= \frac{1}{\hbar V}\int dt\int dt'\int d^3r\int d^3r' e^{i(\omega t-\omega't'-\bm k\cdot\bm r+\bm k'\cdot\bm r')}g^r(\bm r,t,\bm r',t')\\ &= \frac{1}{\hbar V}\int dt\int dt'\int d^3r\int d^3r' e^{i(\omega t-\omega't'-\bm k\cdot\bm r+\bm k'\cdot\bm r')} g_0^r(\bm r,t,\bm r',t')\\ &+ \frac{1}{\hbar V}\int dt\int dt'\int d^3r\int d^3r' e^{i(\omega t-\omega't'-\bm k\cdot\bm r+\bm k'\cdot\bm r')} \frac{1}{\hbar}\int dt_1\int d^3r_1 g_0^r(\bm r,t,\bm r_1,t_1)v_i(\bm r_1)g^r(\bm r_1,t_1,\bm r',t')\\ &= g_{0\bm k,\bm k',\omega,\omega'}^r+ \frac{1}{\hbar^2 V}\int dt_1\int d^3r_1 \int dt\int dt'\int d^3r\int d^3r' e^{i(\omega t_1-\omega't'-\bm k\cdot\bm r_1+\bm k'\cdot\bm r')} \\ &\hspace{1cm}\times e^{i\omega (t-t_1)}e^{-i\bm k\cdot(\bm r-\bm r_1)} g_0^r(\bm r-\bm r_1,t-t_1)v_i(\bm r_1)g^r(\bm r_1,t_1,\bm r',t')\\ &= 2\pi\delta(\omega-\omega')\delta_{\bm k,\bm k'}g_{0\bm k,\omega}^r+ \frac{1}{\hbar V}\int dt_1\int d^3r_1 \int dt'\int d^3r' e^{i(\omega t_1-\omega't'-\bm k\cdot\bm r_1+\bm k'\cdot\bm r')}\\ &\hspace{1cm}\times \left[\frac{1}{\hbar}\int d^3(r-r_1)\int d(t-t_1)e^{-i\bm k\cdot(\bm r-\bm r_)}e^{i\omega (t-t_1)} g_0^r(\bm r-\bm r_1,t-t_1)\right]v_i(\bm r_1)g^r(\bm r_1,t_1,\bm r',t')\\ &=2\pi\delta(\omega-\omega')\delta_{\bm k,\bm k'}g_{0\bm k,\omega}^r+ \frac{g_{0\bm k,\omega}^r}{\hbar V}\int dt_1\int d^3r_1 \int dt'\int d^3r' e^{i(\omega t_1-\omega't'-\bm k\cdot\bm r_1+\bm k'\cdot\bm r')} \Big[\sum_{\bm q}e^{-i\bm q\cdot\bm r_1}v_i(\bm q) \Big] g^r(\bm r_1,t_1,\bm r',t')\\ &=2\pi\delta(\omega-\omega')\delta_{\bm k,\bm k'}g_{0\bm k,\omega}^r+ \frac{g_{0\bm k,\omega}^r}{\hbar V} \sum_{\bm q}v_i(\bm q) \int d^3r_1\int d^3r'\int dt_1\int dt'e^{i(\omega t_1-\omega't'-i(\bm k+\bm q)\cdot\bm r_1+\bm k'\cdot\bm r')} g^r(\bm r_1,t_1,\bm r',t')\\ &=2\pi\delta(\omega-\omega')\delta_{\bm k,\bm k'}g_{0\bm k,\omega}^r+ g_{0\bm k,\omega}^r \sum_{\bm q}v_i(\bm q) g_{\bm k+\bm q,\bm k',\omega,\omega'}^r\\ &=2\pi\delta(\omega-\omega')\left[\delta_{\bm k,\bm k'}g_{0\bm k,\omega}^r+ g_{0\bm k,\omega}^r \sum_{\bm q}v_i(\bm q) g_{\bm k+\bm q,\bm k',\omega}^r\right]

4つ目の等号は、 g_0(\bm r,t,\bm r',t')=g_0(\bm r-\bm r',t-t') と表し、 \bm r''=\bm r-\bm r' t''=t-t' と変換することで、

&g_{0\bm k,\bm k',\omega,\omega'} \\ &= \frac{1}{V}\int d^3r\int d^3r' \int dt \int dt'\ e^{i\bm k\cdot \bm r} e^{-i\bm k'\cdot \bm r'} e^{-i\omega t} e^{i\omega' t'} g_0(\bm r-\bm r',t-t') \\ &= \frac{1}{V}\int d^3r' \int d^3r'' \int dt' \int dt''\, e^{i\bm k\cdot (\bm r'+\bm r'')} e^{-i\bm k'\cdot \bm r'} e^{-i\omega (t'+t'')} e^{i\omega' t'} g_0(\bm r'',t'') \\ &= \frac{1}{V} \underbrace{\int d^3r' e^{-i(\bm k'-\bm k)\cdot \bm r'}}_{V\delta_{\bm k,\bm k'}} \underbrace{\int dt' e^{i(\omega'-\omega) t'}}_{2\pi\delta(\omega-\omega')} \underbrace{\int d^3r'' \int dt'' e^{i\bm k\cdot \bm r''} e^{-i\omega t''} g_0(\bm r'',t'')}_{g_{0\bm k,\omega}} \\ &= 2\pi\delta(\omega-\omega') \delta_{\bm k,\bm k'} f_{0\bm k,\omega}

と表せるが、9/26 のセミナーで植田先生、佐野先生に、上記の変数変換に加えて \bm R=\frac{\bm r+\bm r'}{2} T=\frac{t+t'}{2} の変換も行った方が物理的に見通しが良くなることを教えていただいた。

この置き換えにより、

\bm r=\bm R+\bm r''/2

\bm r'=\bm R-\bm r''/2

t=T+t''/2

t'=T-t''/2

となるから・・・

(後でさらに追記予定)

したがって、

&g_{\bm k,\bm k',\omega}^r=\delta_{\bm k,\bm k'}g_{0\bm k,\omega}^r+ g_{0\bm k,\omega}^r \sum_{\bm q}v_i(\bm q) g_{\bm k+\bm q,\bm k',\omega}^r

このように見直してみると、 (8.74)〜(8.76) や (8.116) のフーリエ変換の定義で係数はいろいろ工夫してあった。 V を埋め込む場所が c v_i で異なったり、 数学的には意味のない \hbar が埋め込まれていたり。

これらの工夫は上式で余計な係数が出てこないようにするための、 細心の注意を払った定義であったことがここへ来て始めて分かる。

もう一度次元を確認しておくと、

  • g_{\bm k,\bm k',\omega} および g_{\bm k,\omega} は (1/エネルギー) の次元
  • v_i(\bm q) は (エネルギー) の次元を持っている

不純物ポテンシャルの次数で展開する

右辺の g^r を (8.117) で繰り返し展開すると、

(8.118)

&g_{\bm k,\bm k',\omega}^r=\delta_{\bm k,\bm k'}g_{0\bm k,\omega}^r+ g_{0\bm k,\omega}^r \sum_{\bm q}v_i(\bm q) \big(\delta_{\bm k+\bm q,\bm k'}g_{0\bm k+\bm q,\omega}^r+ g_{0\bm k+\bm q,\omega}^r \sum_{\bm q'}v_i(\bm q')g_{\bm k+\bm q+\bm q',\bm k',\omega}^r \big)\\ &=g_{0\bm k,\omega}^r\delta_{\bm k,\bm k'}+ \sum_{\bm q}g_{0\bm k,\omega}^r v_i(\bm q) g_{0\bm k+\bm q,\omega}^r \delta_{\bm k+\bm q,\bm k'}+ \sum_{\bm q,\bm q'} g_{0\bm k,\omega}^r v_i(\bm q) g_{0\bm k+\bm q,\omega}^r v_i(\bm q')g_{\bm k+\bm q+\bm q',\bm k',\omega}^r\\ &=g_{0\bm k,\omega}^r\delta_{\bm k,\bm k'}+ \sum_{\bm q_1}g_{0\bm k,\omega}^r v_i(\bm q_1) g_{0\bm k+\bm q_1,\omega}^r \delta_{\bm k+\bm q_1,\bm k'}+ \sum_{\bm q_1,\bm q_2} g_{0\bm k,\omega}^r v_i(\bm q_1)g_{0\bm k+\bm q_1,\omega}^r v_i(\bm q_2)g_{0\bm k+\bm q_1+\bm q_2,\omega}^r \delta_{\bm k+\bm q_1+\bm q_2,\bm k'}+\dots\\ &+\sum_{\bm q_1,\bm q_2,\dots,\bm q_n} g_{0\bm k,\omega}^r v_i(\bm q_1)g_{0\bm k+\bm q_1,\omega}^r v_i(\bm q_2)g_{0\bm k+\bm q_1+\bm q_2,\omega}^r \dots v_i(\bm q_{n-1})g_{0\bm k+\bm q_1+\dots+\bm q_{n-1},\omega}^r v_i(\bm q_n)g_{0\bm k+\bm q_1+\dots+\bm q_n,\omega}\, \delta_{\bm k+\bm q_1+\dots+\bm q_n,\bm k'}\\ &+\dots\\

特徴を見ておくと、

  • n 次の項には n+1 個の g_0 が現われ、
  • その間に n 個の v_i(\bm q_?) が入っている
  • 一番右の \bm k+\bm q_1+\dots+\bm q_n \delta 関数により \bm k' に等しくなるため、 n 個の \sum のうち実質的に自由に動かせるのは n-1 個である。

不純物ポテンシャルの位置平均

不純物の位置平均を次のように定義する。

&\langle \hspace{5 mm}\rangle_i\equiv \frac{1}{V^{N_i}}\int_V d^3R_{1}\int_V d^3R_{2}\dots\int_V d^3R_{N_i}

つまり、不純物が系の体積 V の中のどこにある確率も一定であるという条件である。

(10/3 のセミナー後追記)
佐野先生からこの積分に不純物密度の重み付けを入れることで、 不純物密度に偏りがある場合に対応可能であるとの指摘があった。

この場合、後に出てくる e^{i\bm q_?\cdot\bm R_{k_?}} などと合わせて、 δ関数ではなく不純物密度のフーリエ係数 \bm q_? 成分が残ることになる。 不純物存在確率のフーリエ変換はキュムラント母関数などとも深い関係があるとのこと。

追記ここまで。

v_i(\bm q) v_i(\bm q)v_i(\bm q') の中に現われる和について不純物散乱の位置平均を取る。

(8.119A)

&\langle \sum_{k=1}^{N_i}e^{\textcolor{red}{i}\bm q\cdot\bm R_k}\rangle_i \equiv \frac{1}{V^{N_i}}\int d^3R_{1}\int d^3R_{2}\dots\int d^3R_{N_i} \sum_{k=1}^{N_i}e^{\textcolor{red}{i}\bm q\cdot\bm R_k}\\ &=\sum_{k=1}^{N_i} \frac{1}{V}\int d^3R_k e^{\textcolor{red}{i}\bm q\cdot\bm R_k} = N_i\delta_{\bm q,\bm 0}

&\langle \sum_{k_1,k_2}e^{\textcolor{red}{i}\bm q_1\cdot\bm R_{k_1}}e^{\textcolor{red}{i}\bm q_2\cdot\bm R_{k_2}} \rangle_i\\ &=\frac{1}{V^{N_i}}\int d^3R_{1}\int d^3R_{2}\dots\int d^3R_{N_i} \sum_{k_1,k_2}e^{\textcolor{red}{i}\bm q_1\cdot\bm R_{k_1}}e^{\textcolor{red}{i}\bm q_2\cdot\bm R_{k_2}} \\ &= \sum_{k_1\ne k_2} \frac{1}{V^2} \int d^3R_{k_1}\int d^3R_{k_2} e^{\textcolor{red}{i}\bm q_1\cdot\bm R_{k_1}}e^{\textcolor{red}{i}\bm q_2\cdot\bm R_{k_2}} + \sum_{k} \frac{1}{V} \int d^3R_{k} e^{\textcolor{red}{i}\bm q_1\cdot\bm R_{k}}e^{\textcolor{red}{i}\bm q_2\cdot\bm R_{k}}\\ &= N_i(N_i-1)\delta_{\bm q_1,\bm 0}\delta_{\bm q_2,\bm 0} + N_i\delta_{\bm q_1+\bm q_2,\bm 0}

&\langle \sum_{k_1,k_2,k_3}e^{\textcolor{red}{i}\bm q_1\cdot\bm R_{k_1}}e^{\textcolor{red}{i}\bm q_2\cdot\bm R_{k_2}}e^{\textcolor{red}{i}\bm q_2\cdot\bm R_{k_3}} \rangle_i\\ &= \sum_{k_1\ne k_2 \ne k_3} \delta_{\bm q_1,0}\delta_{\bm q_2,0}\delta_{\bm q_3,0} + \sum_{k_1 = k_2 \ne k_3} \delta_{\bm q_1+\bm q_2,0}\delta_{\bm q_3,0} + \sum_{k_1 \ne k_2 = k_3} \delta_{\bm q_1,0}\delta_{\bm q_2+q_3,0} + \sum_{k_1 = k_3 \ne k_2} \delta_{\bm q_1+\bm q_3,0}\delta_{\bm q_2,0} + \sum_{k_1 = k_3 = k_2} \delta_{\bm q_1+\bm q_2+\bm q_3,0}\\ &= N_i(N_i-1)(N_i-2)\delta_{\bm q_1,\bm 0}\delta_{\bm q_2,\bm 0} +N_i(N_i-1)(\delta_{\bm q_1+\bm q_2,0}\delta_{\bm q_3,0}+\delta_{\bm q_1,0}\delta_{\bm q_2+q_3,0}+\delta_{\bm q_1+\bm q_3,0}\delta_{\bm q_2,0}) +N_i\delta_{\bm q_1+\bm q_2+\bm q_3,\bm 0}

ここで、 v_i(\bm q=\bm 0) の成分は不純物平均を取る前からゼロなので、 \delta_{\bm q_?,\bm 0} を含む項はすべて消えて、

\langle v_i(\bm q_1) \rangle_i=0

\langle v_i(\bm q_1)v_i(\bm q_2) \rangle_i=N_iv_i^2\left(\frac{a^3}{V}\right)^2\delta_{\bm q_1+\bm q_2,\bm 0}=n_iv_i^2\frac{1}{N}\delta_{\bm q_1+\bm q_2,\bm 0}

\langle v_i(\bm q_1)v_i(\bm q_2)v_i(\bm q_3) \rangle_i=n_iv_i^3\frac{1}{N^2}\delta_{\bm q_1+\bm q_2+\bm q_3,\bm 0}

などとなる・・・というのが教科書の論法だけれど、10/17 のセミナーで佐野先生からこれだと落ちている項があるとの指摘があった。

実際、上の第2式で \bm q_1=\bm q_2=\bm o と置けば、 左辺は平均を取る前から明らかにゼロであるにも関わらず、 右辺には n_iv_i^2\frac{1}{N} が残ってしまう。

ちゃんと計算すると、上記の式には、

\langle v_i(\bm q_1)v_i(\bm q_2) \rangle_i=N_iv_i^2\left(\frac{a^3}{V}\right)^2\delta_{\bm q_1+\bm q_2,\bm 0}=n_iv_i^2\frac{1}{N}(\delta_{\bm q_1+\bm q_2,\bm 0}\red{-\delta_{\bm q_1,\bm 0}\delta_{\bm q_2,\bm 0})}

\langle v_i(\bm q_1)v_i(\bm q_2)v_i(\bm q_3) \rangle_i=n_iv_i^3\frac{1}{N^2}(\delta_{\bm q_1+\bm q_2+\bm q_3,\bm 0}\red{-\delta_{\bm q_1,\bm 0}\delta_{\bm q_2,\bm 0}\delta_{\bm q_3,\bm 0})}

など、 \delta_{\bm q_k,\bm 0} の積が残るのが正しい。

これらの項があると、グリーン関数には教科書で説明されている項の他に余計な項が出る。 下の (8-10.4) の結果を先取りすると、 2M 次の主要項が

(g_{0\bm k,\omega})^{M+1}\left(-\pi in_iv_i^2\nu(0)\right)^M

であったところに、

-(g_{0\bm k,\omega})^{M+1}\Big(\frac{n_iv_i^2}{N^{2}}g_{0\bm k,\omega}\Big)^M

の項が余計に現われる。

両者の比を考えると、

&\left(-\pi i\nu(0)\right)^M \>> \Big(\frac{1}{N^{2}}g_{0\bm k,\omega}\Big)^M

が成立すれば右辺を無視できるが、 N は格子点の数で非常に大きいため、 g_{0\bm k,\omega} が発散する点の周辺以外では十分成り立ちそう。

ただ、 g_{0\bm k,\omega} の分母が非常に小さくなるときも無視して良いのかどうか、 完全に自信を持って言い切れないところがある?

(10/17 + 10/24 ここまで追記)

n 次の項では、 n 個の k_? のうち、そのすべてが、 どれか他の k_{?'} と等しくなる項以外は消えることになる。

どの k_? とどの k_? とが等しいかを表すために使われているのが あの不思議な点線の三角形。たとえば n=9 の時の図として、

8-119A.png

が与えられれば、 k_1=k_4=k_7 k_2=k_3=k_5=k_8 k_6=k_9 となることを表している。 これに対応して、δ関数部分は \delta_{\bm q_1+\bm q_4+\bm q_7,\bm 0}\delta_{\bm q_2+\bm q_3+\bm q_5+\bm q_8,\bm 0}\delta_{\bm q_6+\bm q_9,\bm 0} となる。 また、和を取る際に独立に動かせるのは k_1,k_2,k_6 の3つだけであり、しかもそれらは互いに等しくなってはいけないため、 その項数は _{N_i}C_3=N_i(N_i-1)(N_i-2) となる。 3\ll N_i のとき、 _{N_i}C_3=N_i^3 と見なしてしまって構わない。

したがって、上記の図に対応する因子は N_i^3\delta_{\bm q_1+\bm q_4+\bm q_7,\bm 0}\delta_{\bm q_2+\bm q_3+\bm q_5+\bm q_8,\bm 0}\delta_{\bm q_6+\bm q_9,\bm 0} となる。

n=4 であれば、

  • k_1=k_2,k_3=k_4 N_i^2\delta_{\bm q_1+\bm q_2,0}\delta_{\bm q_3+\bm q_4,0}
  • k_1=k_2=k_3=k_4 N_i\delta_{\bm q_1+\bm q_2+\bm q_3+\bm q_4,0}
  • k_1=k_4,k_2=k_3 N_i^2\delta_{\bm q_1+\bm q_4,0}\delta_{\bm q_2+\bm q_3,0}
  • k_1=k_3,k_2=k_4 N_i^2\delta_{\bm q_1+\bm q_3,0}\delta_{\bm q_2+\bm q_4,0}

の4つの場合以外の項は消えて、これらが (8.135)〜(8.139) で評価されている4つの場合に相当する。

運動量の保存 = 波数成分の回復

このようにして現われる \delta 関数の部分に注目すると、 たとえば \delta_{\bm q_1+\bm q_2,0}\delta_{\bm q_3+\bm q_4,0} \bm q_1+\bm q_2=\bm q_3+\bm q_4=0 の時以外ゼロとなるため、 この項を評価するに当たっては、 \bm q_1+\bm q_2+\bm q_3+\bm q_4=0 と仮定して良い。

同様に、 n 次に現われるすべての項で \bm q_1+\bm q_2+\dots+\bm q_n=0 を仮定できる。

すると (8.118) に見るように、すべての n 次の項には \delta_{\bm k+\bm q_1+\bm q_2+\dots+\bm q_n,\bm k'} が含まれているため、 これを \delta_{\bm k,\bm k'} と書き換えることができて、

g_{\bm k,\bm k',\omega}^r\propto\delta_{\bm k,\bm k'}

これは、不純物平均により電子の運動量の保存が回復したことに対応している。 (平均により系の並進対称性が回復したと言い換えても良い、らしい)

そこで、

(8.121)

g_{\bm k,\bm k',\omega}^r\equiv\delta_{\bm k,\bm k'}g_{\bm k,\omega}^r

と書けて、

g_{\bm k,\omega}^r=g_{0\bm k,\omega}^r+ n_iv_i^2\frac{1}{N}\sum_{\bm q_1}g_{0\bm k,\omega}^r g_{0\bm k+\bm q_1,\omega}^r g_{0\bm k,\omega}^r+ n_iv_i^3\frac{1}{N^2}\sum_{\bm q_1,\bm q_2}g_{0\bm k,\omega}^r g_{0\bm k+\bm q_1,\omega}^r g_{0\bm k+\bm q_1+\bm q_2,\omega}^r g_{0\bm k,\omega}^r+ \cdots

などとなる。

2次の項

1次の項は消えるので、2次の項が最低次項になる。

2次の項に現われる \bm q_1 に対する和は、 \bm k_1=\bm k+\bm q_1 を新しい変数と考えると独立に評価することができて、

(8.122)

8-122.png

\frac{1}{N}\sum_{\bm q_1}g_{\textcolor{red}{0}\bm k+\bm q_1,\omega}^r= \frac{1}{N}\sum_{\bm k_1}g_{\textcolor{red}{0}\bm k_1,\omega}^r= a^3\int\frac{d^3k}{(2\pi)^3}g_{\textcolor{red}{0}\bm k_1,\omega}^r= a^3\int\frac{d^3k}{(2\pi)^3}\frac{1}{\hbar\omega-\varepsilon_{\bm k}+i0}

\bm k に対する積分を、 \varepsilon\equiv\varepsilon_{\bm k}=\frac{\hbar^2\bm k^2}{2m}-\varepsilon_F に対する積分に置き換えると、 ( \nu(\varepsilon) は 分散関係)

(8.123)

\frac{1}{N}\sum_{\bm q_1}g_{\textcolor{red}{0}\bm k+\bm q_1,\omega}^r= \int_{-\epsilon_F}^\infty d\varepsilon\nu(\varepsilon)\frac{1}{\hbar\omega-\varepsilon+i0}

ここでの \varepsilon は運動エネルギーではなく、 運動エネルギーからフェルミエネルギーを引いた物であることに注意。

この積分、本当は δ関数の性質 を使うとまともに計算できるんだけど、 まずはしばらく教科書の通りやってみる。

(8.124)

\sum_{\bm k}\ ?=V\int \frac{d^3k}{(2\pi)^3}\ ?=N\int_{-\varepsilon_F}^\infty\nu(\varepsilon)d\varepsilon\ ?

から始めればいいのかな?

\nu(\varepsilon)\equiv\frac{V}{N(2\pi)^3}\frac{d^3k}{d\varepsilon}=\frac{a^3}{2\pi^3}\frac{k^2dk}{d\varepsilon} =\frac{mk}{2\pi^2\hbar^2}a^3 =\frac{m^{3/2}a^3}{\sqrt 2\pi^2\hbar^3}\sqrt{\varepsilon+\varepsilon_F}

N=V/a^3
d^3k=4\pi k^2dk
\frac{d\epsilon}{dk}=\frac{d}{dk}\left(\frac{\hbar^2}{2m}k^2-\varepsilon_F\right)=\frac{\hbar^2}{m}k
k=\frac{1}{\hbar}\sqrt{2m(\varepsilon+\varepsilon_F)}

などを使った。

(8.125)

\frac{1}{\hbar\omega-\varepsilon+i0} \varepsilon\sim\hbar\omega 付近でのみ大きな値を取る。そこで、

\int_{-\epsilon_F}^\infty d\varepsilon\nu(\varepsilon)\frac{1}{\hbar\omega-\varepsilon+i0}\sim \nu(\textcolor{red}{\hbar\omega})\int_{-\infty}^\infty d\varepsilon\frac{1}{\hbar\omega-\varepsilon+i0}

どうして教科書で \nu(\hbar\omega) でなく \nu(0) になっているのか分からない・・・

→ この疑問への答えは 9-1 で出てくる。 \hbar\omega=0 という近似は、 系の温度が低いと仮定し、フェルミ分布関数 f(\omega) をステップ関数と見なす近似に対応している。

\frac{1}{N}\sum_{\bm k}g_{\bm k,\omega}^r\sim \nu(0)\int_{-\infty}^\infty d\varepsilon\frac{1}{-\varepsilon+i0} =\nu(0)\int_{-\infty}^\infty d\varepsilon\frac{-\varepsilon-i0}{\varepsilon+(0)^2} \ \textcolor{red}{\stackrel{?}{=}}\,-\pi i\nu(0)

うーん、これって説明になっているんだろうか?なんだかよく分からない。

(9.18) あたりでもう少しまじめにやると書いてあるので重複するかもしれないけれど、 δ関数の性質 を使って少しちゃんとやってみる。

→ (9.18) は2次の項だけで、1次の項については出てこなかった

&\int_{-\epsilon_F}^\infty d\varepsilon\nu(\varepsilon)\frac{1}{\hbar\omega-\varepsilon+i0}\\ &=-\int_{\epsilon_F+\hbar\omega}^{-\infty} d\varepsilon'\nu(\hbar\omega-\varepsilon')\frac{1}{\varepsilon'+i0}\\ &=\int_{-\infty}^{\epsilon_F+\hbar\omega} d\varepsilon'\nu(\hbar\omega-\varepsilon')\frac{1}{\varepsilon'+i0}\\ &=\int_{-\infty}^{\epsilon_F+\hbar\omega} d\varepsilon'\nu(\hbar\omega-\varepsilon') \left[\frac{1}{\varepsilon'}-i\pi\delta(\varepsilon')\right]\\ &=\left[\ \frac{m^{3/2}a^3}{\sqrt 2\pi^2\hbar^3}\int_{-\infty}^{\epsilon_F+\hbar\omega} d\varepsilon'\frac{\sqrt{\hbar\omega-\varepsilon'+\varepsilon_F}}{\varepsilon'}\ \right] -i\pi\nu(\hbar\omega)\\

この第一項、実数成分の被積分関数は \varepsilon'\rightarrow -\infty の時に \varepsilon'^{1/2} と見なすことができるから、 この積分は発散する。しかし、今は件の積分の虚数成分のみに興味があって、 実数成分は無視してしまって良い・・・と (8.126) の直後の 注11) に書いてある。

しかし、

虚数成分にはやはり \nu(0) ではなく \nu(\hbar\omega) が出てくることが確かめられた。が、以降では(低温近似を先取りして)教科書通り \nu(0) で進めてみる。

(8.127)

&g_{\bm k,\omega}^r=g_{0\bm k,\omega}^r- n_iv_i^2 \pi i\nu(0)(g_{0\bm k,\omega}^r)^2+\cdots\\ &=g_{0\bm k,\omega}^r- \frac{i\hbar}{2\tau}(g_{0\bm k,\omega}^r)^2+\cdots\\

ここで、

(8.128)

\frac{\hbar}{\tau}\equiv 2 n_iv_i^2\pi\nu(0)

この \frac{\hbar}{\tau} は電子エネルギーのぼやけであることが後に分かる。

そこで、系が金属であるための条件

(8.129)

\frac{\hbar}{\tau}= 2\pi n_iv_i^2\nu(0)\ll \varepsilon_F

を先取りして用いていく。

さて、実はここから先は面倒なだけで遠回りかつ不正確な記述が続く、ちゃんと解く方法は この問題、実はとても簡単に解けるのでは? を参照。

ただ、9章を理解するにはやはり以下の内容についてもある程度分かっていないとダメなので、完全に無駄というわけでもない。実際、練習問題と思えば十分に価値がある。

↑ これ、勘違いだったみたい

(10/3 セミナー後の追記)

(8.128) の \hbar を左辺に移すと、

\frac{\hbar}{\tau}\equiv 2 n_iv_i^2\pi\nu(0)

となるが、これとフェルミの黄金律との関係に注目するよう佐野先生から指摘があった。

高次の項 (n>2)

(8.130)

g_{\bm k,\bm k'} を計算するに当たり、

  • 0次には1つも無いが
  • 1次は e^{i\bm q\cdot\bm R_k} を1つ
  • 2次は e^{i\bm q\cdot\bm R_k} を2つ
  • ・・・

含んでいる。

不純物平均により、複数の \bm q_{k_1},\bm q_{k_2},\dots,\bm q_{k_m} に共通の R_{k} がかかり、 \delta_{\bm q_{k_1}+\bm q_{k_2}+\dots+\bm q_{k_m},\bm 0} が現われる項以外が消えてしまうため、

(8.131), (8.132)

のように、 n 個の e^{i\bm q\cdot\bm R_k} を必ず2本以上束ねた項だけが残る。

(8.134)

8-133.png

(8.121) の3次の項に出てくる和は、 \bm k_1=\bm k+\bm q_1 \bm k_2=\bm k+\bm q_1+\bm q_2 と置けば、

&\sum_{\bm q_1,\bm q_2}g_{0\bm k,\omega}^rv_i(\bm q_1)g_{0\bm k+\bm q_1,\omega}^rv_i(\bm q_2)g_{0\bm k+\bm q_1+\bm q_2,\omega}^r v_i(\bm q_3=-\bm q_1-\bm q_2)g_{0\bm k,\omega}^r\\ &=n_iv_i^3\textcolor{red}{\frac{1}{N^2}}\sum_{\bm q_1,\bm q_2}g_{0\bm k,\omega}^r g_{0\bm k+\bm q_1,\omega}^r g_{0\bm k+\bm q_1+\bm q_2,\omega}^r g_{0\bm k,\omega}^r\\ &=n_iv_i^3g_{0\bm k,\omega}^r \Big(\textcolor{red}{\frac{1}{N}}\sum_{\bm k_1}g_{0\bm k_1,\omega}^r\Big)\Big(\textcolor{red}{\frac{1}{N}}\sum_{\bm k_2}g_{0\bm k_2,\omega}^r\Big) g_{0\bm k,\omega}^r\\ &=n_iv_i^3g_{0\bm k,\omega}^r \Big(\textcolor{red}{\frac{1}{N}}\sum_{\bm k_1}g_{0\bm k_1,\omega}^r\Big)^2 g_{0\bm k,\omega}^r\\ &=n_iv_i^3g_{0\bm k,\omega}^r \big(-\pi i\nu(0)\big)^2 g_{0\bm k,\omega}^r\\ &=-\frac{\textcolor{red}{\hbar}}{2\tau}v_i\pi \nu(0)\big(g_{0\bm k,\omega}^r \big)^2\\

(8.135)

8-135.png

\delta_{\bm q_1+\bm q_2,0}\delta_{\bm q_3+\bm q_4,0} すなわち \bm q_1+\bm q_2=\bm q_3+\bm q_4=0 の項で、 N_i(N_i-1)\sim N_i^2 とすれば、

&\sum_{\bm q_1,\bm q_2}g_{0\bm k,\omega}^rv_i(\bm q_1)g_{0\bm k+\bm q_1,\omega}^rv_i(-\bm q_1)g_{0\bm k,\omega}^rv_i(\bm q_3)g_{0\bm k+\bm q_3,\omega}^rv_i(-\bm q_3)g_{0\bm k,\omega}^r\\ &=n_i^2v_i^4\frac{1}{N^2}\sum_{\bm q_1,\bm q_3}g_{0\bm k,\omega}^r g_{0\bm k+\bm q_1,\omega}^r g_{0\bm k,\omega}^r g_{0\bm k+\bm q_3,\omega}^r g_{0\bm k,\omega}^r\\ &=n_i^2v_i^4\Big(\frac{1}{N}\sum_{\bm k_1}g_{0\bm k+\bm k_1,\omega}^r\Big)^2 \big(g_{0\bm k,\omega}^r\big)^3\\ &=\big(-n_iv_i^2\pi i\nu(0)\big)^2 \big(g_{0\bm k,\omega}^r\big)^3\\ &=\Big(\frac{-i\textcolor{red}{\hbar}}{2\tau}\Big)^2\big(g_{0\bm k,\omega}^r\big)^3\\

(8.136)

8-136.png

\delta_{\bm q_1+\bm q_2+\bm q_3+\bm q_4,0} すなわち \bm q_1+\bm q_2+\bm q_3+\bm q_4=0 の項で、

&\sum_{\bm q_1,\bm q_2,\bm q_3}g_{0\bm k,\omega}^rv_i(\bm q_1)g_{0\bm k+\bm q_1,\omega}^rv_i(\bm q_2)g_{0\bm k+\bm q_1+\bm q_2,\omega}^rv_i(\bm q_3)g_{0\bm k+\bm q_1+\bm q_2+\bm q_3,\omega}^rv_i(-\bm q_1-\bm q_2-\bm q_3)g_{0\bm k,\omega}^r\\ &=n_iv_i^4\frac{1}{N^3}\sum_{\bm q_1,\bm q_2,\bm q_3}g_{0\bm k,\omega}^r g_{0\bm k+\bm q_1,\omega}^r g_{0\bm k+\bm q_1+\bm q_2,\omega}^r g_{0\bm k+\bm q_1+\bm q_2+\bm q_3,\omega}^r g_{0\bm k,\omega}^r\\ &=n_iv_i^4\Big(\frac{1}{N}\sum_{\bm k_1}g_{0\bm k_1,\omega}^r\Big)^3 \big(g_{0\bm k,\omega}^r\big)^2\\ &=n_iv_i^4\big(-\pi i\nu(0)\big)^3 \big(g_{0\bm k,\omega}^r\big)^2\\ &=\frac{i\textcolor{red}{\hbar}}{2\tau}\big(v_i\textcolor{red}{\pi} \nu(0)\big)^2 (g_{0\bm k,\omega}^r\big)^2

(8.137)

8-137.png

\delta_{\bm q_1+\bm q_4,0}+\delta_{\bm q_2+\bm q_3,0} すなわち \bm q_1+\bm q_4=\bm q_2+\bm q_3=0 の項で、

&\sum_{\bm q_1,\bm q_2}g_{0\bm k,\omega}^rv_i(\bm q_1)g_{0\bm k+\bm q_1,\omega}^rv_i(\bm q_2)g_{0\bm k+\bm q_1+\bm q_2,\omega}^rv_i(-\bm q_2)g_{0\bm k+\bm q_1,\omega}^rv_i(-\bm q_1)g_{0\bm k,\omega}^r\\ &=n_i^2v_i^4\frac{1}{N^2}\sum_{\bm q_1,\bm q_2}g_{0\bm k,\omega}^r g_{0\bm k+\bm q_1,\omega}^r g_{0\bm k+\bm q_1+\bm q_2,\omega}^r g_{0\bm k+\bm q_1,\omega}^r g_{0\bm k,\omega}^r\\ &=n_i^2v_i^4\Big(\frac{1}{N}\sum_{\bm k_1}(g_{0\bm k_1,\omega}^r)^2\Big)\Big(\frac{1}{N}\sum_{\bm k_2}g_{0\bm k_2,\omega}^r\Big) \big(g_{0\bm k,\omega}^r\big)^2\\ &=n_i^2v_i^4\big(\ ?\ \big)\big(-\pi i\nu(0)\big) \big(g_{0\bm k,\omega}^r\big)^2\\ &=n_iv_i^2\big(\ ?\ \big)\big(n_iv_i^2\cdot-\pi i\nu(0)\big) \big(g_{0\bm k,\omega}^r\big)^2\\ &=-\frac{i\hbar}{2\tau}n_iv_i^2\big(\ ?\ \big)\big(g_{0\bm k,\omega}^r\big)^2\\ &=(*)

(\ ?\ ) の部分が、

&\frac{1}{N}\sum_{\bm k_1}(g_{0\bm k_1,\omega}^r)^2=-\frac{\pi}{2}i\nu(0)\frac{1}{\varepsilon_F}

であれば、教科書の通り、

(*)=-\frac{1}{2}\Big(\frac{\hbar}{2\tau}\Big)^2\frac{1}{\varepsilon_F}\big(g_{0\bm k,\omega}^r\big)^2

となる。

この計算は (9.18) から (9.22) で解説されていた。 この手の計算はちゃんと複素解析的に考えないとダメみたい(汗

教科書は誤植もあるようなので、この積分のやり方を δ関数の性質 にまとめた。

恐らく正しくは、

&\frac{1}{N}\sum_{\bm k_1}(g_{0\bm k_1,\omega}^r)^2 =\textcolor{red}{+}\frac{\pi}{2}i\nu(\hbar\omega)\frac{1}{\hbar\omega+\varepsilon_F}

であり、 \hbar\omega\sim 0 との近似の下、

(*)=\textcolor{red}{+}\frac{1}{2}\Big(\frac{\hbar}{2\tau}\Big)^2\frac{1}{\varepsilon_F}\big(g_{0\bm k,\omega}^r\big)^2

となるような?

(8.138)

8-138A.png

&\sum_{\bm q_1,\bm q_2}g_{0\bm k,\omega}^rv_i(\bm q_1)g_{0\bm k+\bm q_1,\omega}^rv_i(\bm q_2)g_{0\bm k+\bm q_1+\bm q_2,\omega}^rv_i(-\bm q_1)g_{0\bm k+\bm q_2,\omega}^rv_i(-\bm q_2)g_{0\bm k,\omega}^r\\ &=n_i^2v_i^4\frac{1}{N^2}\sum_{\bm q_1,\bm q_2}g_{0\bm k,\omega}^r g_{0\bm k+\bm q_1,\omega}^r g_{0\bm k+\bm q_1+\bm q_2,\omega}^r g_{0\bm k+\bm q_2,\omega}^r g_{0\bm k,\omega}^r\\ &=n_i^2v_i^4\textcolor{red}{\frac{1}{N^2}}\big(g_{0\bm k,\omega}^r\big)^2 \sum_{\bm k_1,\bm k_2} g_{0\bm k_1,\omega}^r \textcolor{red}{g_{0\bm k_1+\bm k_2-\bm k,\omega}^r g_{0\bm k_2,\omega}^r} \\

\bm k_1=\bm k+\bm q_1,\bm k_2=\bm k+\bm q_2 と置いた。

最後の式の g_{0\bm k_1+\bm k_2-\bm k,\omega}^r g_{0\bm k_2,\omega}^r とを入れ替えると教科書の式になる。 入れ替えずに書いておいてくれれば分かりやすいのに。

足のクロスする山の評価について

この次に出てくる説明はいろいろおかしいと思う。たぶん正しいのは以下。

  • g_{0\bm k} |\bm k|=k_F 付近でのみ大きな値を取る。なぜなら分母 \hbar\omega-\varepsilon+i\delta_\varepsilon は低温近似で \hbar\omega=0 と見なせることから、 0-\varepsilon+i\delta_\varepsilon=\textstyle \frac{\hbar^2k^2}{2m}-\varepsilon_F+i\delta_\varepsilon となり、 \textstyle \frac{\hbar^2k_F^2}{2m}=\varepsilon_F を満たす k_F においてピークを持つ。
  • g_{0\bm k} g_{0\bm k_1} g_{0\bm k_2} g_{0\bm k+\bm k_2-\bm k_1} が大きな値を取るのは |\bm k|=|\bm k_1|=|\bm k_2|=|\bm k_1+\bm k_2-\bm k_1|=k_f となるときのみ
  • |\bm k_1+\bm k_2-\bm k|=k_f について考えると、
    • 与えられた -\bm k に対して任意の \bm k_1 を選ぶことができるが、
    • -\bm k_2 |\bm k_1+\bm k_2-\bm k|=k_F となる球面上に来るように選ばなければならない
    • これは原点を中心とした半径 k_F の球殻と、 -\bm k+\bm k_1 を中心とした半径 k_F の球殻との交線上に |\bm k_1+\bm k_2-\bm k|=k_F が来ることと同義である

8-138.png

  • フェルミ面のぼやけを \delta k 程度と仮定すると、
  • \bm k_2 を完全に自由に取れば 4\pi k_F^2\delta k が球殻の体積であるが、
  • g_{0\bm k+\bm k_2-\bm k_1,\omega}^r のために上記の円弧状の交線に限られると、体積はどれほど減るだろうか?
    • -\bm k \bm k_1 とが為す角を \theta とすると、
    • 円弧の半径は k_F\sin(\theta/2) で、円周は 2\pi k_F\sin(\theta/2)
    • 円弧の太さ(断面積)は \frac{\delta k}{\sin\theta}\cdot \delta k=\frac{\delta k^2}{\sin\theta}
    • したがって、円弧の体積は \frac{2\pi k_F\delta k^2\sin(\theta/2)}{\sin\theta}
    • \bm k_1 が動いたことを想定し、 \theta,\psi に対して平均を取ると、
      &\frac{1}{4\pi}\int_0^\pi d\theta\int_0^{2\pi}\sin\theta d\psi \frac{2\pi k_F\delta k^2\sin(\theta/2)}{\sin\theta}\\ &=\frac{2\pi k_F\delta k^2}{2}\int_0^\pi d\theta\sin(\theta/2)\\ &=\frac{2\pi k_F\delta k^2}{2}[-2\cos(\theta/2)]_0^\pi\\ &=2\pi k_F\delta k^2
    • この最終結果は \theta = 90^\circ の結果と(偶然?)一致する。
  • したがって、完全に自由に \bm k_2 を取ったときに比べて値は \frac{\delta k}{2k_F} だけ小さくなる
  • フェルミレベル近傍に於いて \frac{\delta \varepsilon}{\varepsilon_F}=\frac{2\delta k}{k_F} を使うと、この比を \frac{\delta k}{2k_F}=\frac{\delta \varepsilon}{4\varepsilon_F}=\frac{\hbar}{4\varepsilon_F\tau} と表せる。
  • 教科書に現われる \frac{\hbar}{4\varepsilon_F\tau} という比はこの値。

ということで、ちょっと意味は無いけれど本文中の誤植のみ指摘:

\frac{\textcolor{red}{\hbar^2}k_F}{m}\delta k\sim \frac{\hbar}{\tau}

\ell\equiv \frac{\textcolor{red}{\hbar}k_F}{m}\tau

(\hbar/\varepsilon_F^{\red 2}\tau)^{\red 1}

閑話休題して、

(8.139)

「自由に取った場合」を考えるには、 \bm k_2 を数えるときに g が2つ掛かっていることを考慮して (8.137) に帰着する。

&n_i^2v_i^4\textcolor{red}{\frac{1}{N^2}}\big(g_{0\bm k,\omega}^r\big)^2 \sum_{\bm k_1,\bm k_2} g_{0\bm k_1,\omega}^r g_{0\bm k_2,\omega}^r g_{0\bm k+\bm k_2-\bm k_1,\omega}^r \\ &\sim \left(\frac{\hbar}{4\varepsilon_F\tau}\right)n_i^2v_i^4\textcolor{red}{\frac{1}{N^2}}\big(g_{0\bm k,\omega}^r\big)^2 \sum_{\bm k_1,\bm k_2} g_{0\bm k_1,\omega}^r \big(g_{0\bm k_2,\omega}^r\big)^2 \\ &=\left(\frac{\hbar}{4\varepsilon_F\tau}\right) \frac{1}{2}\Big(\frac{\hbar}{2\tau}\Big)^2\frac{1}{\varepsilon_F}\big(g_{0\bm k,\omega}^r\big)^2\\ &\propto \frac{\hbar}{\tau}\left(\frac{\hbar}{\varepsilon_F\tau}\right)^2\big(g_{0\bm k,\omega}^r\big)^2

4次の項の大小関係に関する物理的解釈

10/17 のセミナーで植田先生から教えていただいた内容:

8-137.png

8-138A.png

に比べて、

8-135.png

が大きい物理的理由は次のように考えると納得できる。

8-135.png

\bm q -\bm q \bm q' -\bm q' の2組の散乱がそれぞれ独立して生じたという項であるのに対して、

8-137.png

8-138A.png

\bm q -\bm q の散乱が終わる前に \bm q' の散乱が起きたと見なせる。

散乱頻度が低いときには、複数の散乱が同時に起きる確率は それらが独立に起きる確率に比べて十分低い。

散乱頻度( 1/\tau に相当)が大きくなると、 そのような確率も無視できなくなるが、それを判別する条件が

\hbar/\varepsilon_F \tau \ll 1

である。

ν(0) の大きさ

各項の大きさを評価する前に、 \nu(0) の大きさを評価しておく。

\varepsilon_F がフェルミエネルギーであることから(実際にはこの分を引いてあるため フェルミレベルは 0 である)

&N\int_{-\epsilon_F}^0\nu(\varepsilon)d\varepsilon=N_{\mathrm{electrons}}\\

の関係が成り立つ。 すなわち、フェルミレベル以下の状態数が電子の粒子数に等しい。

現在の所、「格子点数 N 」と「電子の粒子数 N_{\mathrm{electrons}} 」との関係を 良く理解できていないのだけれど、これらがおよそ等しい

N\sim N_{\mathrm{electrons}}

とすれば、

&\int_{-\epsilon_F}^0\nu(\varepsilon)d\varepsilon=1\\ &=\int_{-\epsilon_F}^0 d\varepsilon \frac{m^{3/2}a^3}{\sqrt 2\pi^2\hbar^3}\sqrt{\varepsilon+\varepsilon_F}\\ &=\frac{m^{3/2}a^3}{\sqrt 2\pi^2\hbar^3}\int_0^{\epsilon_F} d\varepsilon' \sqrt{\varepsilon'}\\ &=\frac{m^{3/2}a^3}{\sqrt 2\pi^2\hbar^3}\left[\frac{2}{3}\varepsilon'^{3/2}\right]_0^{\epsilon_F}\\ &=\frac{2\varepsilon_F}{3}\cdot\frac{m^{3/2}a^3}{\sqrt 2\pi^2\hbar^3}\sqrt{\varepsilon_F}\\ &=\frac{2\varepsilon_F}{3}\nu(0)

を得る。

したがって、

\nu(0)=\frac{3}{2\varepsilon_F}

で、教科書にあるとおり \nu(0)\sim 1/\varepsilon_F を得る。 金属であるためには、

\frac{\hbar}{\tau}=2\pi n_iv_i^2\nu(0)=\frac{3\pi n_iv_i^2}{\varepsilon_F}\ll\varepsilon_F

が必要である。

&3\pi n_iv_i^2\ll\varepsilon_F^2\\ &(3\pi n_i)^\frac{1}{2}v_i\ll\varepsilon_F\\

と書き直せば、これは不純物のポテンシャルがフェルミエネルギーに比べて小さく、 さらに不純物密度が高すぎないことを示している。

(8.136) ではこれに加えて、

v_i\ll\varepsilon_F\sim 1/\nu(0)

すなわち

v_i\nu(0)\ll 1

が要求されている。

(以下 10/24 追記)

\nu(0)\sim 1/\varepsilon_F を考えるとこれは

v_i\ll \varepsilon_F

を表しており、上記と比べると (3\pi n_i)^\frac{1}{2} のファクターだけ厳しい条件になっている。

すなわち、 (3\pi n_i)^\frac{1}{2}\ll 1 であるから、 v_i\ll \varepsilon_F が成立するならば (3\pi n_i)^\frac{1}{2}v_i\ll\varepsilon_F は当然成立する。

金属である条件として与えられている 1/εFτ << 1 について

2章に与えられている値を参照すると(P17 表2.1)、

\rho\ (\mathrm{\mu\Omega\,cm}) \textstyle\frac{\hbar}{\varepsilon_F\tau}
Fe 7.0 10^{-2}
Ni 9.8 10^{-2}
Au 2.2 5\times 10^{-3}
Fe81Ni19 50\sim 60 10^{-1}

注)Fe81Ni19 はパーマロイ

となって、それほど極端に小さいわけでもないことには注意が必要。

となると・・・

3\pi n_iv_i^2/\varepsilon_F^2= 10^{-1}\sim 10^{-3}

より、

\frac{v_i}{\nu (0)}=\sqrt{\frac{10^{-1}\sim 10^{-3}}{3\pi n_i}}

ここで n_i=\frac{N_i}{N}\sim 10^{-4} あたりなので、

\frac{v_i}{\nu (0)}\ll 1

は、とてもじゃないけど成り立ちそうにないけれど?!

N 〜 Nelectrons について

N\sim N_{\mathrm{electrons}} が成り立つ意味は何だろうか?

まず、 \bm k 空間における格子点の密度は (8.77) あたりの式にも出てくるように、

\sum_{\bm k}\Big[\ \ ?\ \ \Big] = \left(\frac{V}{(2\pi)^3}\right)\int d^3k\Big[\ \ ?\ \ \Big]

であり、これは1粒子の波動関数が体積 V を持つ系の大きさに 対応する境界条件から、取り得る \bm k の値が離散化することに対応している。

この式には a^3 N も現われないことに注意が必要。

すなわち、 N の表す「格子点数」と \bm k 空間の「格子点数」とは関係がない。

では N とは何だったか?

a^3 は (8.106) で始めて現われた体積で、 本来δ関数状の不純物散乱ポテンシャルを a^3 にならしたときの高さとして v_i を定義したのであった。この段階では a^3 には物理的意味が与えられておらず、 v_ia^3 を一定に保つ限り a^3 の取り方には任意性が残っていた。

教科書ではその後、(8.121) にて N=V/a^3 として N を導入し、 これを「全格子点数」と呼んでいる。つまり、 a^3 は実空間における結晶格子の基本単位胞体積を考えていたことが読み取れる。

さかのぼって考えると、 v_i は不純物散乱の強さを基本単位胞体積にならしたときの 平均高さとして定義される。この点は v_i \varepsilon との大きさを比べる際に重要となりそう。

すなわち、 N の定義は「実空間において系に金属の単位胞がいくつ含まれるか」である。

すると例えば n 価の面心立方金属であれば単位胞当たり 4n 個の自由電子を生じ、

N_\mathrm{electrons}=4n\left(\frac{V}{a^3}\right)=4nN

である。体心であれば

N_\mathrm{electrons}=2nN

である。実際の数は考える金属の価数と結晶構造に依るけれど、オーダー的には

N_\mathrm{electrons}\sim N

と見なせることが分かる。

(追記ここまで)

どの項が支配的か

上記を仮定すると、

(8.140)

g_{\bm k,\omega}^r= g_{0\bm k,\omega}^r+\left(-\frac{i\hbar}{2\tau}\right)(g_{0\bm k,\omega}^r)^2 \left[ 1\textcolor{red}{-}i\pi v_i\nu(0)-(\textcolor{red}{\pi}v_i\nu(0))^2 -\frac{i}{4}\frac{\hbar}{\varepsilon_F\tau}+O\left(\frac{\hbar}{\varepsilon_F\tau}\right)^2+\cdots \right]+ \left(-\frac{i\hbar}{2\tau}\right)^2(g_{0\bm k,\omega}^r)^3+\cdots

の [ ] の中は1を残して消し去ることができる。

同様にして、高次項を評価すれば2次の項やそれ以降にも小さな項が付け加わるが、 上記の2つの条件下ではやはり無視できて、

(8.141)

\frac{g_{\bm k,\omega}^r}{g_{0\bm k,\omega}^r}= 1+\left(-\frac{i\hbar}{2\tau}g_{0\bm k,\omega}^r\right)+ \left(-\frac{i\hbar}{2\tau}g_{0\bm k,\omega}^r\right)^2+\cdots =\frac{1}{1+\frac{i\hbar}{2\tau}g_{0\bm k,\omega}^r}

g_{\bm k,\omega}^r= \frac{1}{(g_{0\bm k,\omega}^r)^{-1}+\frac{i\hbar}{2\tau}} =\frac{1}{\hbar\omega-\varepsilon_{\bm k}+\frac{i\hbar}{2\tau}}

となる。

本文中で \frac{i\hbar}{2\tau}\ll 1 などと書いてあるが、 これは \textcolor{red}{\frac{\hbar}{2\varepsilon_F\tau}}\ll 1 の間違いか?

結局、元々の偶数次項から出てくる「2次の寄与と同様な物が繰り返される項」、 すなわち \bm q_1+\bm q_2=\bm q_3+\bm q_4=\dots=\bm q_{2n-1}+\bm q_{2n}=0 となる項が支配的な寄与を及ぼし、他の項は無視できることになる。

支配的とされる項に赤丸を付けつつ、6次の途中まで列挙してみた。 6次はまだまだたくさんあり、7次も非常にたくさん出てくるが、 次の支配項は8次の項になる。

8-140.png

以下、6次の項がまだまだ続く・・・

(8.142)

g_{\bm k,\omega}^a =\frac{1}{\hbar\omega-\varepsilon_{\bm k}-\frac{i\hbar}{2\tau}}

不純物散乱がない場合には i0 のために非常に急峻に変化する 関数であったが、不純物散乱によってピーク幅が広がっている、 すなわちフェルミレベルがぼやけている、ということになる。

自己エネルギー

「自己エネルギー」についてはここだけ読んでも物理的意味がよく分からないので後で勉強が必要。

(8.143)

\Sigma^r=-\frac{i\hbar}{2\tau}

\Sigma^a=+\frac{i\hbar}{2\tau}

と置いて、これらを自己エネルギーと呼べば、

(8.144)

g_{\bm k,\omega}^r =\frac{1}{\hbar\omega-\varepsilon_{\bm k}-\Sigma^r}

g_{\bm k,\omega}^a =\frac{1}{\hbar\omega-\varepsilon_{\bm k}-\Sigma^a}

\Sigma の部分が虚数であることが一目で分からないので注意が必要。

計算時には、

g_{\bm k,\omega}^r =\frac{1}{\hbar\omega-\varepsilon_{\bm k}+i\delta_\varepsilon}

g_{\bm k,\omega}^a =\frac{1}{\hbar\omega-\varepsilon_{\bm k}-i\delta_\varepsilon}

ただし、 \delta_\varepsilon=\frac{\hbar}{2\tau}

という表示が分かりやすい。

10/31 のセミナーで「自己エネルギー」という言葉の意味について教えていただいた。

g_{\bm k,\omega}^\alpha =\frac{1}{\hbar\omega-\varepsilon_{\bm k}-\Sigma^\alpha}

の形を見ると、自由な場合の Green 関数と比べて

\varepsilon_{\bm k}\rightarrow\varepsilon_{\bm k}+\Sigma^\alpha

とした形になっている。

すなわち、不純物散乱の効果を取り込むことによって系のエネルギーが \Sigma^\alpha だけ増加しており、これを系が自ら持つエネルギーという意味で 自己エネルギーと呼んでいるのだと思う。

高次項に対する独自考察

一般的にn次を考える

n 次の項で、重なり合いながら存在する m 個の山のそれぞれが l_j 本の点線を束ねたものであったとする。

(8-10.1)

\sum_{j=1}^m l_j=n

たとえば下図は m=6,n=21 の場合である。

8-10-1.png

(8.119) にあたる式において、

  • Σを取るべき k_i の個数は m 個である。
  • それぞれの k_i は異なっていなければならないため、現われる係数は _{N_i}C_m=N_i(N_i-1)(N_i-2)\dots(N_i-(m-1)) である。
    (コメントで指摘をいただいたのですが、ここは _{N_i}C_m=N_i(N_i-1)(N_i-2)\dots(N_i-(m-1))/m! としなければならないのに間違っていたようです。以降の議論にどう影響するか未確認なのでお気を付け下さい。すみません。)
  • m<<N_i では _{N_i}C_m \sim N_i^m と見なせる。

したがって、その項は

(8-10.2)

&\frac{N_i^m}{N^n}v_i^n \sum_{\underbrace{^{\bm q,\bm q,\dots,\bm q}}_{n-m}} g_{0\bm k}^\alpha \underbrace{ g_{0\bm k+\bm q}^\alpha \dots g_{0\bm k+\bm q\dots}^\alpha }_{n-1} g_{0\bm k}^\alpha \\&= g_{0\bm k}^\alpha \Bigg[ \frac{n_i^mv_i^n}{N^{n-m}} \sum_{\underbrace{^{\bm q,\bm q,\dots,\bm q}}_{n-m}} \underbrace{ g_{0\bm k+\bm q}^\alpha \dots g_{0\bm k+\bm q\dots}^\alpha }_{n-1}\Bigg]g_{0\bm k}^\alpha

のような形になる。

連山ごとに和が取れる

この項が完全に M 個の連山に分かれているとき、 (山と山との間が標高ゼロになるとき、完全に山が分かれていると呼ぶ) それら連山ごとに別々に和が取れて、 間に g_{0\bm k}^\alpha が挟まる形になる。

それぞれの連山に含まれる単山の個数を m_j
それぞれの連山に含まれる点線の本数を n_j
とすると \sum_j m_j=m \sum_j n_j=n であり、

8-10-3.png

(8-10.3)

&= g_{0\bm k}^\alpha \Bigg[\frac{n_i^{m_1}v_i^{n_1}}{N^{n_1-m_1}}\sum_{\underbrace{^{\bm q,\bm q,\dots,\bm q}}_{n_1-m_1}} \underbrace{ g_{0\bm k+\bm q}^\alpha \dots g_{0\bm k+\bm q\dots}^\alpha }_{n_1-1}\Bigg] g_{0\bm k}^\alpha \Bigg[\frac{n_i^{m_2}v_i^{n_2}}{N^{n_2-m_2}}\sum_{\underbrace{^{\bm q,\bm q,\dots,\bm q}}_{n_2-m_2}} \underbrace{ g_{0\bm k+\bm q}^\alpha \dots g_{0\bm k+\bm q\dots}^\alpha }_{n_2-1}\Bigg] g_{0\bm k}^\alpha \dots \\&= g_{0\bm k}^\alpha \underbrace{ \Big[\hspace{4mm}\hspace{4mm}\Big]g_{0\bm k}^\alpha \Big[\hspace{4mm}\hspace{4mm}\Big]g_{0\bm k}^\alpha \cdots \Big[\hspace{4mm}\hspace{4mm}\Big] }_{M}g_{0\bm k}^\alpha \\&= \left(g_{0\bm k}^\alpha\right)^{M+1} \underbrace{ \Big[\hspace{4mm}\hspace{4mm}\Big] \Big[\hspace{4mm}\hspace{4mm}\Big] \dots \Big[\hspace{4mm}\hspace{4mm}\Big] }_{M}

支配項との比較

すべての連山が最も単純な2本の線からなる山でできている場合、 その項は支配的な項としてカウントされているタイプになる。

8-10-4.png

このとき、 n=2M であるから、その項は以下の様に評価できる。

(8-10.4)

&=\left(g_{0\bm k}^\alpha\right)^{M+1} \underbrace{ \Big[\hspace{4mm}(m_1=1,n_1=2)\hspace{4mm}\Big] \dots \Big[\hspace{4mm}(m_M=1,n_M=2)\hspace{4mm}\Big] }_{M} \\ &=\left(g_{0\bm k}^\alpha\right)^{M+1} \Big[\frac{n_iv_i^2}{N}\sum_{\bm q'}g_{0\bm q'\omega}^\alpha\Big]^M\\ &=\left(g_{0\bm k}^\alpha\right)^{M+1} \Big[n_iv_i^2\big(-\pi i\nu(0)\big)\Big]^M\\ &=\left(g_{0\bm k}^\alpha\right)^{M+1} \Big[-\frac{i\hbar}{2\tau}\ \Big]^M

したがって、この項が本当に支配的であることを言うためには、 他の \left(g_{0\bm k}^\alpha\right)^{M+1} を出す項すべてを加えても、この項が十分に大きいと見なせる ことを確認しなければならない。

\left(g_{0\bm k}^\alpha\right)^{M+1} を出すのは、 M 個の連山となる項なので、

8-10-5.png

連山ごとの比較に帰着する

上記の n 次、 m 個の山、 M 個の連山を持つ項と比べると、 すべての単連山について以下の式が成り立てば良いことになる。 (連山に含まれる単山の個数を m' 、点線の本数を n' とする)

8-10-5B.png

(8-10.5)

\Bigg|\ \frac{n_i^{m'}v_i^{n'}}{N^{n'-m'}}\sum_{\underbrace{^{\bm q,\bm q,\dots,\bm q}}_{n'-m'}} \underbrace{ g_{0\bm k+\bm q}^\alpha \dots g_{0\bm k+\bm q\dots}^\alpha }_{n'-1}\Bigg| \stackrel{?}{\ll} \frac{n_iv_i^2}{N}\sum_{\bm q'}g_{0\bm q'\omega}^\alpha =\Bigg|-\frac{i\hbar}{2\tau}\ \Bigg|

これを証明するには、(8.138) で見た様な足のクロスした連山が、 クロスした部分の足を1本にまとめた n'' 本足の単山より ずっと小さいことを示した上で、さらに n''>2 本足の単山が 2 本足の単山より小さいことを言えば良い。 あ、クロスせずに (8.137) のように入れ子になっている場合も考える必要がありますね。

8-10-6.png

単山からなる場合

この後半は (8.134) や (8.136) のようにすれば簡単で、 n''>2 本足の単山 ( m'=1 ) は、

(8-10.6)

&\frac{n_iv_i^{n''}}{N^{n''-1}}\sum_{\underbrace{^{\bm q,\bm q,\dots,\bm q}}_{n''-1}} \underbrace{ g_{0\bm k+\bm q}^\alpha \dots g_{0\bm k+\bm q\dots}^\alpha }_{n''-1} \\&= \Big[\frac{n_iv_i^2}{N}\sum_{\bm q'}g_{0\bm q'}\Big] \Big[\frac{v_i}{N}\sum_{\bm q'}g_{0\bm q'}\Big]^{n''-2} \\&= \Big[-\frac{i\hbar}{2\tau}\Big] \Big[-\pi iv_i\nu(0)\Big]^{n''-2}

となって、 v_i\nu(0)\ll 1 との前提を置けば支配項に比べて無視できることが分かる。

前記手順の前半は、クロスしている部分については (8.138) から (8.139) の部分でやったような考え方で、入れ子になっている部分については (8.137) でやったような考え方で、評価すれば良いのだと思うけれど、 計算手順について、まだあまり納得できていない。

入れ子も個別に和が取れる

入れ子について:外側の連山のある足の直後に、入れ子になった小さな連山が M' 個あったとすると、内側の連山は個別に和を取ることができて、

8-10-7.png

(8-10.7)

&\frac{n_i^{m'}v_i^{n'}}{N^{n'-m'}}\sum_{\underbrace{^{\bm q,\bm q,\dots,\bm q}}_{n'-m'}} \underbrace{ g_{0\bm k+\bm q}^\alpha \dots g_{0\bm k+\bm q\dots}^\alpha }_{n'-1} \\&= \Bigg[\frac{n_i^{m'_o}v_i^{n'_o}}{N^{n'_o-m'_o}}\sum_{\underbrace{^{\bm q,\bm q,\dots,\bm q}}_{n'_o-m'_o}} \underbrace{ g_{0\bm k+\bm q}^\alpha \dots g_{0\bm k+\bm q\dots}^\alpha }_{n'_o-1}\Bigg] g_{0\bm k_{\mathrm{prev}}}^\alpha \Bigg[\frac{n_i^{m'_{i1}}v_i^{n'_{i1}}}{N^{n'_{i1}-m'_{i1}}}\sum_{\underbrace{^{\bm q,\bm q,\dots,\bm q}}_{n'_{i1}-m'_{i1}}} \underbrace{ g_{0\bm k+\bm q}^\alpha \dots g_{0\bm k+\bm q\dots}^\alpha }_{n'_{i1}-1}\Bigg]\\ &\hspace{2cm} g_{0\bm k_{\mathrm{prev}}}^\alpha \Bigg[\frac{n_i^{m'_{i2}}v_i^{n'_{i2}}}{N^{n'_{i2}-m'_{i2}}}\sum_{\underbrace{^{\bm q,\bm q,\dots,\bm q}}_{n'_{i2}-m'_{i2}}} \underbrace{ g_{0\bm k+\bm q}^\alpha \dots g_{0\bm k+\bm q\dots}^\alpha }_{n'_{i2}-1}\Bigg]\dots \\&= \Bigg[\frac{n_i^{m'_o}v_i^{n'_o}}{N^{n'_o-m'_o}}\sum_{\underbrace{^{\bm q,\bm q,\dots,\bm q}}_{n'_o-m'_o}} \underbrace{ g_{0\bm k+\bm q}^\alpha \dots (g_{0\bm k_{\mathrm{prev}} q\dots}^\alpha)^{M'+1} \dots g_{0\bm k+\bm q\dots}^\alpha }_{n'_o-1}\Bigg] \Bigg[\ \ (i1)\ \ \Bigg]\dots\Bigg[\ \ (iM')\ \ \Bigg]

などとなる。

ここで、 m'_o,n'_o は外側の連山の、 m'_{i1},n'_{i1},m'_{i2},n'_{i2},\dots は内側の連山の、山の数および線の数。 g_{0\bm k_{\mathrm{prev}}}^\alpha は入れ子連山の直前の外側の山の足に対応する g である。

通常ならば連山と連山の間には g_{0\bm k}^\alpha がはさまるが、 入れ子の場合には代わりに g_{0\bm k_{\mathrm{prev}}}^\alpha がはさまる。このため、外側連山の和を取る際に g_{0\bm k_{\mathrm{prev}}}^\alpha の因子が (内側連山の数 + 1) 乗になる。それが全体の和をどれほど小さくするのか、 (8.137) の評価方法を学んでからもう一度考えたい。

項の大きさだけでなく数も数えないといけないのでは?

もう1つ、本当は n 次の項の大きさが「小さい」ことだけでなく、 その数がそれほど「多くない」ことも言わなければならないはずだけれど・・・ n 次の項から M 連山の項が何個出るかを数えるのは かなり大変。

加えて、ちょっと考えただけでも そんなに自信を持って「多くない」とは言い切れない感じがしている。

例えば 13次の項から出る6連山の項は 12次の項から出る6連山の支配項に比べて、 1つであれば支配項より v_i\nu(0) だけ小さく、無視できるかもしれない。しかし下図の様に、 3本足の山をどこに置くかで同じ大きさの項が6つ出てくることから、 すべてを合わせると比率は 6v_i\nu(0) まで縮まってしまう。

8-10-6B.png

同様にして、 2M+1 次の項から出る M 連山の項を全て合わせれば 対応する支配項の Mv_i\nu(0) 倍になるから、 v_i\nu(0) が有限である限り M が大きいところで大小は逆転する。

教科書はこの辺曖昧なままになっているけれど、 実際には高次の項は「小さいから無視できる」というだけの話では無く、 「互いに打ち消し合ってくれるために表に出てこない」 のが本当のところなのだと思う。

ここで用いた近似の意味

上で用いた近似の意味を考えてみる。
(この項は 2011/10/17 まで間違えた内容を書いてあったのを、 植田先生に指摘していただいて直している最中です)

(8.117) を変形していく。

(8-10.8)

&\langle g_{\bm k,\omega}^\alpha\rangle_i = g_{0\bm k,\omega}^\alpha+ g_{0\bm k,\omega}^\alpha\sum_{\bm q}\langle v_i(\bm q) g_{\bm k+\bm q,\omega}^\alpha \rangle_i \\&= g_{0\bm k,\omega}^\alpha+ g_{0\bm k,\omega}^\alpha\sum_{\bm q}\left\langle v_i(\bm q) \left[ g_{0\bm k+\bm q,\omega}^\alpha+ g_{0\bm k+\bm q,\omega}^\alpha\sum_{\bm q'}v_i(\bm q') g_{\bm k+\bm q+\bm q',\omega}^\alpha \right]\right\rangle_i \\&= g_{0\bm k,\omega}^\alpha+ g_{0\bm k,\omega}^\alpha\sum_{\bm q}\left\langle v_i(\bm q)\right\rangle_i g_{0\bm k+\bm q,\omega}^\alpha + g_{0\bm k,\omega}^\alpha\sum_{\bm q,\bm q'} g_{0\bm k+\bm q,\omega}^\alpha \left\langle v_i(\bm q) v_i(\bm q') g_{\bm k+\bm q+\bm q',\omega}^\alpha\right\rangle_i

この式の \left\langle v_i(\bm q) v_i(\bm q') g_{\bm k+\bm q+\bm q',\omega}^\alpha\right\rangle_i について、 v_i(\bm q) v_i(\bm q') の部分と g_{\bm k+\bm q+\bm q',\omega}^\alpha の部分とでそれぞれ別々に不純物平均を取れるとすると、

&\sim g_{0\bm k,\omega}^\alpha+ g_{0\bm k,\omega}^\alpha\sum_{\bm q}\langle v_i(\bm q)\rangle_i g_{0\bm k+\bm q,\omega}^\alpha + g_{0\bm k,\omega}^\alpha\sum_{\bm q,\bm q'}\langle v_i(\bm q)v_i(\bm q')\rangle_i g_{0\bm k+\bm q,\omega}^\alpha \langle g_{\bm k+\bm q+\bm q',\omega}^\alpha\rangle_i \\&= g_{0\bm k,\omega}^\alpha+ g_{0\bm k,\omega}^\alpha\sum_{\bm q} 0 \cdot g_{0\bm k+\bm q,\omega}^\alpha + g_{0\bm k,\omega}^\alpha\sum_{\bm q,\bm q'}\frac{N_iv_i^2}{N^2}\delta_{\bm q+\bm q',\bm 0} g_{0\bm k+\bm q,\omega}^\alpha \langle g_{\bm k+\bm q+\bm q',\omega}^\alpha\rangle_i \\&= g_{0\bm k,\omega}^\alpha+ g_{0\bm k,\omega}^\alpha\frac{n_iv_i^2}{N}\sum_{\bm q} g_{0\bm k+\bm q,\omega}^\alpha \langle g_{\bm k+\bm q+\bm q',\omega}^\alpha\rangle_i \\&= g_{0\bm k,\omega}^\alpha+ g_{0\bm k,\omega}^\alpha \left[\frac{n_iv_i^2}{N}\sum_{\bm k'} g_{0\bm k',\omega}^\alpha\right] \langle g_{\bm k+\bm q+\bm q',\omega}^\alpha\rangle_i \\&= g_{0\bm k,\omega}^\alpha + \Sigma^\alpha g_{0\bm k,\omega}^\alpha \langle g_{\bm k+\bm q+\bm q',\omega}^\alpha\rangle_i

ゴールが見えてきた。

(8-10.9)

&1/g_{0\bm k,\omega}^\alpha = 1/\langle g_{\bm k,\omega}^\alpha\rangle_i + \Sigma^\alpha \\ &1/\langle g_{\bm k,\omega}^\alpha\rangle_i = 1/g_{0\bm k,\omega}^\alpha - \Sigma^\alpha \\ &\langle g_{\bm k,\omega}^\alpha\rangle_i = \frac{1}{1/g_{0\bm k,\omega}^\alpha - \Sigma^\alpha} = \frac{1}{\hbar\omega-\varepsilon_{\bm k} - \Sigma^\alpha}

が得られる。

このことは、(8-10.8) の途中の式で右辺の \langle g_{\bm k,\omega}^\alpha\rangle_i に対してく繰り返し左辺を代入することで、

(8-10.10)

\langle g_{\bm k,\omega}^\alpha\rangle_i =& g_{0\bm k,\omega}^\alpha+ g_{0\bm k,\omega}^\alpha\Big\{\sum_{\bm q}\langle v_i(\bm q)g_{0\bm k+\bm q,\omega}^\alpha v_i(-\bm q)\rangle_i\Big\} \langle g_{\bm k,\omega}^\alpha\rangle_i \\=& g_{0\bm k,\omega}^\alpha+ g_{0\bm k,\omega}^\alpha\Big\{\sum_{\bm q}\langle v_i(\bm q)g_{0\bm k+\bm q,\omega}^\alpha v_i(-\bm q)\rangle_i\Big\} g_{0\bm k,\omega}^\alpha \\ &+ g_{0\bm k,\omega}^\alpha\Big\{\sum_{\bm q}\langle v_i(\bm q)g_{0\bm k+\bm q,\omega}^\alpha v_i(-\bm q)\rangle_i\Big\} g_{0\bm k,\omega}^\alpha\Big\{\sum_{\bm q}\langle v_i(\bm q)g_{0\bm k+\bm q,\omega}^\alpha v_i(-\bm q)\rangle_i\Big\}g_{0\bm k,\omega}^\alpha \\ &+ g_{0\bm k,\omega}^\alpha\Big\{\sum_{\bm q}\langle v_i(\bm q)g_{0\bm k+\bm q,\omega}^\alpha v_i(-\bm q)\rangle_i\Big\} g_{0\bm k,\omega}^\alpha \Big\{\sum_{\bm q}\langle v_i(\bm q)g_{0\bm k+\bm q,\omega}^\alpha v_i(-\bm q)\rangle_i\Big\}g_{0\bm k,\omega}^\alpha \Big\{\sum_{\bm q}\langle v_i(\bm q)g_{0\bm k+\bm q,\omega}^\alpha v_i(-\bm q)\rangle_i\Big\}g_{0\bm k,\omega}^\alpha \\ &+\dots

となって、確かに教科書で扱った支配項のみをすべて含んでいることからも確認できる。

んー、つまり・・・あんまり面白い話にはならなかった感じ?

高次の項を考える場合

以下は 11/7 のゼミで出た話

近似を上のようにとらえると、2次までではなく高次の項を入れた近似について考えることも容易になる・・・かと思いきや、そんなこともないみたい???

すなわち、

上記では (8.117) 2次まで展開した後、「不純物平均の分割」という近似を用いた。

そこで同様に、3次や4次まで展開した後に近似を用いれば、より高次の項を取り込むことができる・・・と思ったのだけれど、やってみるとそうはうまく行かない。

たとえば4次まで取り込んだ場合:

(8-10.11)

\langle g_{\bm k,\omega}^\alpha\rangle_i =& \, g_{0\bm k,\omega}^\alpha+ g_{0\bm k,\omega}^\alpha\sum_{\bm q}\left\langle v_i(\bm q)\right\rangle_i g_{0\bm k+\bm q,\omega}^\alpha + g_{0\bm k,\omega}^\alpha\sum_{\bm q,\bm q'} g_{0\bm k+\bm q,\omega}^\alpha g_{0\bm k+\bm q+\bm q',\omega}^\alpha \left\langle v_i(\bm q) v_i(\bm q') \right\rangle_i \\ &+ g_{0\bm k,\omega}^\alpha\sum_{\bm q,\bm q',\bm q''} g_{0\bm k+\bm q,\omega}^\alpha g_{0\bm k+\bm q+\bm q',\omega}^\alpha g_{0\bm k+\bm q+\bm q'+\bm q'',\omega}^\alpha \left\langle v_i(\bm q) v_i(\bm q') v_i(\bm q'') \right\rangle_i \\ &+ g_{0\bm k,\omega}^\alpha\sum_{\bm q,\bm q',\bm q'',\bm q'''} g_{0\bm k+\bm q,\omega}^\alpha g_{0\bm k+\bm q+\bm q',\omega}^\alpha g_{0\bm k+\bm q+\bm q'+\bm q'',\omega}^\alpha \left\langle v_i(\bm q) v_i(\bm q') v_i(\bm q'') v_i(\bm q''') g_{\bm k+\bm q+\bm q'+\bm q''+\bm q''',\omega}^\alpha \right\rangle_i \\ \sim& \, g_{0\bm k,\omega}^\alpha+ g_{0\bm k,\omega}^\alpha\sum_{\bm q}\left\langle v_i(\bm q)\right\rangle_i g_{0\bm k+\bm q,\omega}^\alpha + g_{0\bm k,\omega}^\alpha\sum_{\bm q,\bm q'} g_{0\bm k+\bm q,\omega}^\alpha g_{0\bm k+\bm q+\bm q',\omega}^\alpha \left\langle v_i(\bm q) v_i(\bm q') \right\rangle_i \\ &+ g_{0\bm k,\omega}^\alpha\sum_{\bm q,\bm q',\bm q''} g_{0\bm k+\bm q,\omega}^\alpha g_{0\bm k+\bm q+\bm q',\omega}^\alpha g_{0\bm k+\bm q+\bm q'+\bm q'',\omega}^\alpha \left\langle v_i(\bm q) v_i(\bm q') v_i(\bm q'') \right\rangle_i \\ &+ g_{0\bm k,\omega}^\alpha\sum_{\bm q,\bm q',\bm q'',\bm q'''} g_{0\bm k+\bm q,\omega}^\alpha g_{0\bm k+\bm q+\bm q',\omega}^\alpha g_{0\bm k+\bm q+\bm q'+\bm q'',\omega}^\alpha \left\langle v_i(\bm q) v_i(\bm q') v_i(\bm q'') v_i(\bm q''')\right\rangle_i \left\langle g_{\bm k+\bm q+\bm q'+\bm q''+\bm q''',\omega}^\alpha \right\rangle_i

これを (8-10.10) と同様に順次展開すると、

  • 0次
  • 2次
  • 3次
  • 4次
  • 4+2次
  • 4+3次
  • 4+4次
  • 4+4+2次
  • 4+4+3次
  • 4+4+4次
  • ・・・

といった項が現われるが、2+3次 や 3+2次 など、非常に多くの項を取りこぼしてしまうことになる。

うーん、やはり植田先生のおっしゃった通り、 この形で理解しようとしても、一筋縄ではいかないみたいだ。。。

4次以下の単山からなる項をすべて取り込むには、上記ではなく、

(8-10.12)

\langle g_{\bm k,\omega}^\alpha\rangle_i \sim& \, g_{0\bm k,\omega}^\alpha\sum_{\bm q,\bm q'} g_{0\bm k+\bm q,\omega}^\alpha \left\langle v_i(\bm q) v_i(\bm q') \right\rangle_i \left\langle g_{\bm k+\bm q+\bm q',\omega}^\alpha \right\rangle_i \\ &+ g_{0\bm k,\omega}^\alpha\sum_{\bm q,\bm q',\bm q''} g_{0\bm k+\bm q,\omega}^\alpha g_{0\bm k+\bm q+\bm q',\omega}^\alpha \left\langle v_i(\bm q) v_i(\bm q') v_i(\bm q'') \right\rangle_i \left\langle g_{\bm k+\bm q+\bm q'+\bm q'',\omega}^\alpha \right\rangle_i \\ &+ g_{0\bm k,\omega}^\alpha\sum_{\bm q,\bm q',\bm q'',\bm q'''} g_{0\bm k+\bm q,\omega}^\alpha g_{0\bm k+\bm q+\bm q',\omega}^\alpha g_{0\bm k+\bm q+\bm q'+\bm q'',\omega}^\alpha \left\langle v_i(\bm q) v_i(\bm q') v_i(\bm q'') v_i(\bm q''')\right\rangle_i \left\langle g_{\bm k+\bm q+\bm q'+\bm q''+\bm q''',\omega}^\alpha \right\rangle_i

の形が真っ先に思い浮かぶが、これだと数え落としはない物の、今度は数えすぎが生じる。

なぜなら (8.135) で扱った 2+2 の形の 4次 の項は、 文字通り 2次 の項を2回重ねることでも生じるため、 このような項を取り込みすぎてしまうことになる。

したがって正しく求めるには、4次 の形から 2+2次 を抜いておく必要がある。

(8-10.13)

\langle g_{\bm k,\omega}^\alpha\rangle_i \sim& \, g_{0\bm k,\omega}^\alpha+ g_{0\bm k,\omega}^\alpha\sum_{\bm q,\bm q'} g_{0\bm k+\bm q,\omega}^\alpha \left\langle v_i(\bm q) v_i(\bm q') \right\rangle_i \left\langle g_{\bm k+\bm q+\bm q',\omega}^\alpha \right\rangle_i \\ &+ g_{0\bm k,\omega}^\alpha\sum_{\bm q,\bm q',\bm q''} g_{0\bm k+\bm q,\omega}^\alpha g_{0\bm k+\bm q+\bm q',\omega}^\alpha \left\langle v_i(\bm q) v_i(\bm q') v_i(\bm q'') \right\rangle_i \left\langle g_{\bm k+\bm q+\bm q'+\bm q'',\omega}^\alpha \right\rangle_i \\ &+ g_{0\bm k,\omega}^\alpha\sum_{\bm q,\bm q',\bm q'',\bm q'''} g_{0\bm k+\bm q,\omega}^\alpha g_{0\bm k+\bm q+\bm q',\omega}^\alpha g_{0\bm k+\bm q+\bm q'+\bm q'',\omega}^\alpha \left\langle v_i(\bm q) v_i(\bm q') v_i(\bm q'') v_i(\bm q''')\right\rangle_i \big( 1 - \delta_{\bm q+\bm q',\bm o}\delta_{\bm q''+\bm q''',\bm o} \big) \left\langle g_{\bm k+\bm q+\bm q'+\bm q''+\bm q''',\omega}^\alpha \right\rangle_i

4次の項が2つの2次の項からなるのは \bm q+\bm q'=\bm o \bm q''+\bm q'''=\bm o の時なので、これを除くために ( 1 - \delta_{\bm q+\bm q',\bm o}\delta_{\bm q''+\bm q''',\bm o} ) という係数を導入した。

不純物平均から出てくるδ関数により右辺の \left\langle g_{\bm k+\bm q+\bm q'+\bm q'',\omega}^\alpha \right\rangle_i などはすべて \left\langle g_{\bm k,\omega}^\alpha \right\rangle_i に等しくなるため、

(8-10.14)

\langle g_{\bm k,\omega}^\alpha\rangle_i \sim& \, g_{0\bm k,\omega}^\alpha+ g_{0\bm k,\omega}^\alpha\Big\{\sum_{\bm q,\bm q'} g_{0\bm k+\bm q,\omega}^\alpha \left\langle v_i(\bm q) v_i(\bm q') \right\rangle_i \Big\} \left\langle g_{\bm k,\omega}^\alpha \right\rangle_i \\ &+ g_{0\bm k,\omega}^\alpha\Big\{\sum_{\bm q,\bm q',\bm q''} g_{0\bm k+\bm q,\omega}^\alpha g_{0\bm k+\bm q+\bm q',\omega}^\alpha \left\langle v_i(\bm q) v_i(\bm q') v_i(\bm q'') \right\rangle_i \Big\} \left\langle g_{\bm k,\omega}^\alpha \right\rangle_i \\ &+ g_{0\bm k,\omega}^\alpha\Big\{\sum_{\bm q,\bm q',\bm q'',\bm q'''} g_{0\bm k+\bm q,\omega}^\alpha g_{0\bm k+\bm q+\bm q',\omega}^\alpha g_{0\bm k+\bm q+\bm q'+\bm q'',\omega}^\alpha \left\langle v_i(\bm q) v_i(\bm q') v_i(\bm q'') v_i(\bm q''')\right\rangle_i \big( 1 - \delta_{\bm q+\bm q',\bm o}\delta_{\bm q''+\bm q''',\bm o} \big) \Big\} \left\langle g_{\bm k,\omega}^\alpha \right\rangle_i \\ =& g_{0\bm k,\omega}^\alpha+g_{0\bm k,\omega}^\alpha \Sigma' \left\langle g_{\bm k,\omega}^\alpha \right\rangle_i

したがって、

(8-10.15)

&\langle g_{\bm k,\omega}^\alpha\rangle_i = \frac{1}{\hbar\omega-\varepsilon_{\bm k} - \Sigma'}

ただし、

(8-10.16)

\Sigma' =& \sum_{\bm q,\bm q'} g_{0\bm k+\bm q,\omega}^\alpha \left\langle v_i(\bm q) v_i(\bm q') \right\rangle_i +\sum_{\bm q,\bm q',\bm q''} g_{0\bm k+\bm q,\omega}^\alpha g_{0\bm k+\bm q+\bm q',\omega}^\alpha \left\langle v_i(\bm q) v_i(\bm q') v_i(\bm q'') \right\rangle_i \\ &+ \sum_{\bm q,\bm q',\bm q'',\bm q'''} g_{0\bm k+\bm q,\omega}^\alpha g_{0\bm k+\bm q+\bm q',\omega}^\alpha g_{0\bm k+\bm q+\bm q'+\bm q'',\omega}^\alpha \left\langle v_i(\bm q) v_i(\bm q') v_i(\bm q'') v_i(\bm q''')\right\rangle_i \big( 1 - \delta_{\bm q+\bm q',\bm o}\delta_{\bm q''+\bm q''',\bm o} \big)

となる(んじゃないかな)

この右辺のうち第一項は2次までを取り入れた際の自己エネルギーそのものであり、 第二項、第三項がそれぞれ3次および4次を取り入れたことによる修正項になっている。

すなわち高次の項を取り入れたとしても、自己エネルギーの値が少し(?)変わるだけで Green 関数の関数形そのものは変化しない。(どこかに発散する項がない限り)

もともと自己エネルギーの値は正確な値の不明な v_i なる値に依存しており、 また実数成分はフェルミエネルギーの値に吸収されてしまうこともあり、 上で疑問に思った「近似の精度」についてはそれほど目くじらを立てる必要も無いのかもしれない?
→ んー、 g^r g^a g^< などの相対的な大きさは重要なので、やはりそう雑には扱えない。

どこまでで打ち切るにせよ、「一貫性のある近似」になっていないとダメだ。

8-11 で見るように、 g^< を2次で打ち切ったり、 4次で打ち切ったりした結果を上記のような手法で求めると、 「一貫性を持った近似」を得られるようだ。

質問・コメント




無題

kt? ()

この本を読むうえで非常に参考にさせていただいております。勉強になります。

質問なのですが、不純物ポテンシャルの位置平均のところや、高次項に対する独自考察のところで、NiC3、NiCmなどコンビネーションを計算しているところで、3!やm!で割らなくてよい理由はなぜなのでしょうか。

コンビネーションを用いるのはわかるのですが、その場合普通のコンビネーションの計算だとそれぞれ3!、m!で割らなければいけない気がします。

  • 記憶が曖昧なところがあるのですが、今見る限りおっしゃるとおり m! などで割らないといけないように思えます。単にミスをしていた可能性が高いです。ちょっとその影響を調べる余力がないので、本文中にコメントのみ追加することにいたします。 -- 武内(管理人)?
  • お忙しいところ早速のご返信ありがとうございます。承知しました。 -- kt?

添付ファイル: file8-138A.png 1013件 [詳細] file8-137.png 937件 [詳細] file8-136.png 975件 [詳細] file8-135.png 1000件 [詳細] file8-133.png 972件 [詳細] file8-122.png 912件 [詳細] file8-10-7.png 863件 [詳細] file8-119A.png 994件 [詳細] file8-10-6B.png 776件 [詳細] file8-10-6.png 806件 [詳細] file8-10-4.png 825件 [詳細] file8-10-5B.png 780件 [詳細] file8-10-5.png 816件 [詳細] file8-10-3.png 789件 [詳細] file8-10-1.png 978件 [詳細] file8-140.png 1030件 [詳細] file8-138.png 1100件 [詳細]

Counter: 7324 (from 2010/06/03), today: 2, yesterday: 3