スピントロニクス理論の基礎/8-10 の履歴(No.23)
更新- 履歴一覧
- 差分 を表示
- 現在との差分 を表示
- ソース を表示
- スピントロニクス理論の基礎/8-10 へ行く。
[前の章へ] <<<< スピントロニクス理論の基礎(目次) >>>> [次の章へ]
8-10 不純物散乱のもとでの Green 関数†
(8.113) の Fourier 変換表示を求める†
(8.116)
&math( &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) 等に与えられた の Fourier 変換とは の配置や の符号などが異なるので注意が必要。
これに対応するハミルトニアンは、
(8.116)
&math( &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) );
となる。
今考えているようなポテンシャルでは電子の振動数が変わらないため、単一の によって表せる。これは以下で自然に導出される。また実際には「不純物平均」により も変わらないという結果になる。
&math(g_{\bm k,\bm k',\omega,\omega'}= 2\pi \delta(\omega-\omega')g_{\bm k,\bm k',\omega});
これらは 8-7 で見た と同じ性質である。
&math(g_{0\bm k,\bm k',\omega,\omega'}= 2\pi\delta(\omega-\omega')\delta_{\bm k,\bm k'}g_{0\bm k,\omega});
→ 単一の や で表せるという結果は実時間・実空間で見れば 系の運動が絶対的な時刻や座標に依存せず、相対値のみで記述できることに対応する。
(8.117)
係数を正しく付けるために次元について確認:
- は (8.76) によれば で無次元、 で無次元、 が残って (時間/エネルギー) の次元を持つ
- は が外に出るので、(1/エネルギー) の次元を持つ
- は (8.116) によれば と同じ次元で、(8.106) によればこれはポテンシャルではなく (エネルギー) の次元を持つ
(8.47)〜(8.77) および (8.116) に与えられたフーリエ係数の定義は係数の付け方などが数学の授業とは異なるので、 慎重に計算していかないと正しい答えが得られない。実際に計算を進めるとかなり大変。
&math( &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つ目の等号は、 と表し、 、 と変換することで、
&math( &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 のセミナーで植田先生、佐野先生に、上記の変数変換に加えて 、 の変換も行った方が物理的に見通しが良くなることを教えていただいた。
この置き換えにより、
となるから・・・
(後でさらに追記予定)
したがって、
&math(&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) のフーリエ変換の定義で係数はいろいろ工夫してあった。 を埋め込む場所が と で異なったり、 数学的には意味のない が埋め込まれていたり。
これらの工夫は上式で余計な係数が出てこないようにするための、 細心の注意を払った定義であったことがここへ来て始めて分かる。
もう一度次元を確認しておくと、
- および は (1/エネルギー) の次元
- は (エネルギー) の次元を持っている
不純物ポテンシャルの次数で展開する†
右辺の を (8.117) で繰り返し展開すると、
(8.118)
&math( &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\\ );
特徴を見ておくと、
- 次の項には 個の が現われ、
- その間に 個の が入っている
- 一番右の は 関数により に等しくなるため、 個の のうち実質的に自由に動かせるのは 個である。
不純物ポテンシャルの位置平均†
不純物の位置平均を次のように定義する。
つまり、不純物が系の体積 の中のどこにある確率も一定であるという条件である。
(10/3 のセミナー後追記)
佐野先生からこの積分に不純物密度の重み付けを入れることで、
不純物密度に偏りがある場合に対応可能であるとの指摘があった。
この場合、後に出てくる などと合わせて、 δ関数ではなく不純物密度のフーリエ係数 成分が残ることになる。 不純物存在確率のフーリエ変換はキュムラント母関数などとも深い関係があるとのこと。
追記ここまで。
や の中に現われる和について不純物散乱の位置平均を取る。
(8.119A)
&math(&\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});
&math( &\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} );
&math( &\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} );
ここで、 の成分は不純物平均を取る前からゼロなので、 を含む項はすべて消えて、
などとなる・・・というのが教科書の論法だけれど、10/17 のセミナーで佐野先生からこれだと落ちている項があるとの指摘があった。
実際、上の第2式で と置けば、 左辺は平均を取る前から明らかにゼロであるにも関わらず、 右辺には が残ってしまう。
ちゃんと計算すると、上記の式には、
など、 の積が残るのが正しい。
これらの項があると、グリーン関数には教科書で説明されている項の他に
&math( \frac{n_iv_i^{n-1}}{N^n}g_{0\bm k,\omega}\left[\sum_{\bm q} \delta_{\bm q,\bm 0} g_{0\bm k+\bm q,\omega}\right]^n );
のような項が余計に現われることになる。
しかしこの項は教科書で扱われている項に比べて少なくとも
程度小さいことになる。
は格子点の数で非常に大きいため、この項は無視できるのだと思う。
← 教科書で扱う項では
の評価結果に
が残らないところに注意
(よく考えると の分母が非常に小さくなるときも無視して良いのかどうか、 完全に自信を持って言い切れないところがある? の次数が上がっている分、分母が小さい時はより大きな値を取るはず。もう少し検討が必要かも)
(10/17 ここまで追記)
次の項では、 個の のうち、そのすべてが、 どれか他の と等しくなる項以外は消えることになる。
どの とどの とが等しいかを表すために使われているのが あの不思議な点線の三角形。たとえば の時の図として、
が与えられれば、 、 、 となることを表している。 これに対応して、δ関数部分は となる。 また、和を取る際に独立に動かせるのは の3つだけであり、しかもそれらは互いに等しくなってはいけないため、 その項数は となる。 のとき、 と見なしてしまって構わない。
したがって、上記の図に対応する因子は となる。
であれば、
- →
- →
- →
- →
の4つの場合以外の項は消えて、これらが (8.135)〜(8.139) で評価されている4つの場合に相当する。
運動量の保存 = 波数成分の回復†
このようにして現われる 関数の部分に注目すると、 たとえば は の時以外ゼロとなるため、 この項を評価するに当たっては、 と仮定して良い。
同様に、 次に現われるすべての項で を仮定できる。
すると (8.118) に見るように、すべての 次の項には が含まれているため、 これを と書き換えることができて、
これは、不純物平均により電子の運動量の保存が回復したことに対応している。 (平均により系の並進対称性が回復したと言い換えても良い、らしい)
そこで、
(8.121)
と書けて、
&math(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次の項に現われる に対する和は、 を新しい変数と考えると独立に評価することができて、
(8.122)
&math( \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} );
に対する積分を、 に対する積分に置き換えると、 ( は 分散関係)
(8.123)
&math( \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} );
ここでの は運動エネルギーではなく、 運動エネルギーからフェルミエネルギーを引いた物であることに注意。
この積分、本当は δ関数の性質 を使うとまともに計算できるんだけど、 まずはしばらく教科書の通りやってみる。
(8.124)
&math( \sum_{\bm k}\ ?=V\int \frac{d^3k}{(2\pi)^3}\ ?=N\int_{-\varepsilon_F}^\infty\nu(\varepsilon)d\varepsilon\ ? );
から始めればいいのかな?
&math( \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} );
などを使った。
(8.125)
は 付近でのみ大きな値を取る。そこで、
&math( \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} );
どうして教科書で でなく になっているのか分からない・・・
→ この疑問への答えは 9-1 で出てくる。 という近似は、 系の温度が低いと仮定し、フェルミ分布関数 をステップ関数と見なす近似に対応している。
&math(\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次の項については出てこなかった
&math( &\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)\\ );
この第一項、実数成分の被積分関数は の時に と見なすことができるから、 この積分は発散する。しかし、今は件の積分の虚数成分のみに興味があって、 実数成分は無視してしまって良い・・・と (8.126) の直後の 注11) に書いてある。
しかし、
虚数成分にはやはり ではなく が出てくることが確かめられた。が、以降では(低温近似を先取りして)教科書通り で進めてみる。
(8.127)
&math( &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)
この は電子エネルギーのぼやけであることが後に分かる。
そこで、系が金属であるための条件
(8.129)
を先取りして用いていく。
さて、実はここから先は面倒なだけで遠回りかつ不正確な記述が続く、ちゃんと解く方法は この問題、実はとても簡単に解けるのでは? を参照。
ただ、9章を理解するにはやはり以下の内容についてもある程度分かっていないとダメなので、完全に無駄というわけでもない。実際、練習問題と思えば十分に価値がある。
↑ これ、勘違いだったみたい
(10/3 セミナー後の追記)
(8.128) の を左辺に移すと、
となるが、これとフェルミの黄金律との関係に注目するよう佐野先生から指摘があった。
高次の項 (n>2)†
(8.130)
を計算するに当たり、
- 0次には1つも無いが
- 1次は を1つ
- 2次は を2つ
- ・・・
含んでいる。
不純物平均により、複数の に共通の がかかり、 が現われる項以外が消えてしまうため、
(8.131), (8.132)
のように、 個の を必ず2本以上束ねた項だけが残る。
(8.134)
(8.121) の3次の項に出てくる和は、 、 と置けば、
&math( &\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)
すなわち の項で、 とすれば、
&math( &\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)
すなわち の項で、
&math( &\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)
すなわち の項で、
&math( &\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\\ &=(*) );
の部分が、
&math( &\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} );
であれば、教科書の通り、
&math( (*)=-\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) で解説されていた。 この手の計算はちゃんと複素解析的に考えないとダメみたい(汗
教科書は誤植もあるようなので、この積分のやり方を δ関数の性質 にまとめた。
恐らく正しくは、
&math( &\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} );
であり、 との近似の下、
&math( (*)=\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)
&math(&\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} \\);
と置いた。
最後の式の と とを入れ替えると教科書の式になる。 入れ替えずに書いておいてくれれば分かりやすいのに。
足のクロスする山の評価について†
この次に出てくる説明はいろいろおかしいと思う。たぶん正しいのは以下。
- は 付近でのみ大きな値を取る。なぜなら分母 は低温近似で と見なせることから、 となり、 を満たす においてピークを持つ。
- が大きな値を取るのは となるときのみ
-
について考えると、
- 与えられた に対して任意の を選ぶことができるが、
- は となる球面上に来るように選ばなければならない
- これは原点を中心とした半径 の球殻と、 を中心とした半径 の球殻との交線上に が来ることと同義である
- フェルミ面のぼやけを 程度と仮定すると、
- を完全に自由に取れば が球殻の体積であるが、
-
のために上記の円弧状の交線に限られると、体積はどれほど減るだろうか?
- と とが為す角を とすると、
- 円弧の半径は で、円周は
- 円弧の太さ(断面積)は
- したがって、円弧の体積は
-
が動いたことを想定し、
に対して平均を取ると、
&math(&\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); - この最終結果は の結果と(偶然?)一致する。
- したがって、完全に自由に を取ったときに比べて値は だけ小さくなる
- フェルミレベル近傍に於いて を使うと、この比を と表せる。
- 教科書に現われる という比はこの値。
ということで、ちょっと意味は無いけれど本文中の誤植のみ指摘:
閑話休題して、
(8.139)
「自由に取った場合」を考えるには、 を数えるときに が2つ掛かっていることを考慮して (8.137) に帰着する。
&math( &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 のセミナーで植田先生から教えていただいた内容:
や
に比べて、
が大きい物理的理由は次のように考えると納得できる。
は と 、 と の2組の散乱がそれぞれ独立して生じたという項であるのに対して、
や
は と の散乱が終わる前に の散乱が起きたと見なせる。
散乱頻度が低いときには、複数の散乱が同時に起きる確率は それらが独立に起きる確率に比べて十分低い。
散乱頻度( に相当)が大きくなると、 そのような確率も無視できなくなるが、それを判別する条件が
である。
ν(0) の大きさ†
各項の大きさを評価する前に、 の大きさを評価しておく。
がフェルミエネルギーであることから(実際にはこの分を引いてあるため フェルミレベルは 0 である)
&math( &N\int_{-\epsilon_F}^0\nu(\varepsilon)d\varepsilon=N_{\mathrm{electrons}}\\ );
の関係が成り立つ。 すなわち、フェルミレベル以下の状態数が電子の粒子数に等しい。
現在の所、「格子点数 」と「電子の粒子数 」との関係を 良く理解できていないのだけれど、これらがおよそ等しい
とすれば、
&math( &\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) );
を得る。
したがって、
で、教科書にあるとおり を得る。 金属であるためには、
が必要である。
&math( &3\pi n_iv_i^2\ll\varepsilon_F^2\\ &(3\pi n_i)^\frac{1}{2}v_i\ll\varepsilon_F\\ );
と書き直せば、これは不純物のポテンシャルがフェルミエネルギーに比べて小さく、 さらに不純物密度が高すぎないことを示している。
これに加えて、
すなわち
どの項が支配的か†
上記を仮定すると、
(8.140)
&math( 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)
&math( \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} );
&math( 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}} );
となる。
本文中で などと書いてあるが、 これは の間違いか?
結局、元々の偶数次項から出てくる「2次の寄与と同様な物が繰り返される項」、 すなわち となる項が支配的な寄与を及ぼし、他の項は無視できることになる。
支配的とされる項に赤丸を付けつつ、6次の途中まで列挙してみた。 6次はまだまだたくさんあり、7次も非常にたくさん出てくるが、 次の支配項は8次の項になる。
以下、6次の項がまだまだ続く・・・
(8.142)
&math( g_{\bm k,\omega}^a =\frac{1}{\hbar\omega-\varepsilon_{\bm k}-\frac{i\hbar}{2\tau}} );
不純物散乱がない場合には のために非常に急峻に変化する 関数であったが、不純物散乱によってピーク幅が広がっている、 すなわちフェルミレベルがぼやけている、ということになる。
自己エネルギー†
「自己エネルギー」についてはここだけ読んでも物理的意味がよく分からないので後で勉強が必要。
(8.143)
と置いて、これらを自己エネルギーと呼べば、
(8.144)
&math( g_{\bm k,\omega}^r =\frac{1}{\hbar\omega-\varepsilon_{\bm k}-\Sigma^r} );
&math( g_{\bm k,\omega}^a =\frac{1}{\hbar\omega-\varepsilon_{\bm k}-\Sigma^a} );
の部分が虚数であることが一目で分からないので注意が必要。
計算時には、
&math( g_{\bm k,\omega}^r =\frac{1}{\hbar\omega-\varepsilon_{\bm k}+i\delta_\varepsilon} );
&math( g_{\bm k,\omega}^a =\frac{1}{\hbar\omega-\varepsilon_{\bm k}-i\delta_\varepsilon} );
ただし、
という表示が分かりやすい。
高次項に対する独自考察†
一般的にn次を考える†
次の項で、重なり合いながら存在する 個の山のそれぞれが 本の点線を束ねたものであったとする。
(8-10.1)
たとえば下図は の場合である。
(8.119) にあたる式において、
- Σを取るべき の個数は 個である。
- それぞれの は異なっていなければならないため、現われる係数は である。
- では と見なせる。
したがって、その項は
(8-10.2)
&math( &\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 );
のような形になる。
連山ごとに和が取れる†
この項が完全に 個の連山に分かれているとき、 (山と山との間が標高ゼロになるとき、完全に山が分かれていると呼ぶ) それら連山ごとに別々に和が取れて、 間に が挟まる形になる。
それぞれの連山に含まれる単山の個数を
それぞれの連山に含まれる点線の本数を
とすると
、
であり、
(8-10.3)
&math( &= 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)
&math( &=\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 );
したがって、この項が本当に支配的であることを言うためには、 他の を出す項すべてを加えても、この項が十分に大きいと見なせる ことを確認しなければならない。
を出すのは、 個の連山となる項なので、
連山ごとの比較に帰着する†
上記の 次、 個の山、 個の連山を持つ項と比べると、 すべての単連山について以下の式が成り立てば良いことになる。 (連山に含まれる単山の個数を 、点線の本数を とする)
(8-10.5)
&math( \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本にまとめた 本足の単山より ずっと小さいことを示した上で、さらに 本足の単山が 本足の単山より小さいことを言えば良い。 あ、クロスせずに (8.137) のように入れ子になっている場合も考える必要がありますね。
単山からなる場合†
この後半は (8.134) や (8.136) のようにすれば簡単で、 本足の単山 ( ) は、
(8-10.6)
&math( &\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} );
となって、 との前提を置けば支配項に比べて無視できることが分かる。
前記手順の前半は、クロスしている部分については (8.138) から (8.139) の部分でやったような考え方で、入れ子になっている部分については (8.137) でやったような考え方で、評価すれば良いのだと思うけれど、 計算手順について、まだあまり納得できていない。
入れ子も個別に和が取れる†
入れ子について:外側の連山のある足の直後に、入れ子になった小さな連山が 個あったとすると、内側の連山は個別に和を取ることができて、
(8-10.7)
&math( &\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] );
などとなる。
ここで、 は外側の連山の、 は内側の連山の、山の数および線の数。 は入れ子連山の直前の外側の山の足に対応する である。
通常ならば連山と連山の間には がはさまるが、 入れ子の場合には代わりに がはさまる。このため、外側連山の和を取る際に の因子が (内側連山の数 + 1) 乗になる。それが全体の和をどれほど小さくするのか、 (8.137) の評価方法を学んでからもう一度考えたい。
項の大きさだけでなく数も数えないといけないのでは?†
もう1つ、本当は 次の項の大きさが「小さい」ことだけでなく、 その数がそれほど「多くない」ことも言わなければならないはずだけれど・・・ 次の項から 連山の項が何個出るかを数えるのは かなり大変。
加えて、ちょっと考えただけでも そんなに自信を持って「多くない」とは言い切れない感じがしている。
例えば 13次の項から出る6連山の項は 12次の項から出る6連山の支配項に比べて、 1つであれば支配項より だけ小さく、無視できるかもしれない。しかし下図の様に、 3本足の山をどこに置くかで同じ大きさの項が6つ出てくることから、 すべてを合わせると比率は まで縮まってしまう。
同様にして、 次の項から出る 連山の項を全て合わせれば 対応する支配項の 倍になるから、 が有限である限り が大きいところで大小は逆転する。
教科書はこの辺曖昧なままになっているけれど、 実際には高次の項は「小さいから無視できる」という話では無く、 以下に見るように「きれいに打ち消し合って出てこない」 のが本当のところなのだと思う。
ここで用いた近似の意味†
上で用いた近似の意味を考えてみる。
(この項は 2011/10/17 まで間違えた内容を書いてあったのを、
植田先生に指摘していただいて直している最中です)
(8.117) を変形していく。
(8-10.8)
&math( &\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'}\left\langle v_i(\bm q) g_{0\bm k+\bm q,\omega}^\alpha v_i(\bm q') g_{\bm k+\bm q+\bm q',\omega}^\alpha\right\rangle_i );
この式の &math(\left\langle v_i(\bm q) g_{0\bm k+\bm q,\omega}^\alpha v_i(\bm q') g_{\bm k+\bm q+\bm q',\omega}^\alpha\right\rangle_i); について、 の部分と の部分とでそれぞれ別々に不純物平均を取れるとすると、
&math( &\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)
&math( &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} );
が得られる。
9章でやるのと同じに考えると、上での近似は先頭の2つの の項と、 以降の、 に含まれる との間で等しくなる項がない、すなわちまたがる「点線」がないということを表していて、 2つの からなる単純な山の繰り返しの項のみを取り出したという、 教科書で紹介された近似そのものである。
んー、つまり・・・あんまり面白い話にはならなかった感じ。