スピントロニクス理論の基礎/9-3 の履歴(No.2)
更新[前の章へ] <<<< スピントロニクス理論の基礎(目次) >>>> [次の章へ]
9-3 拡散型因子†
や に現われた という因子は、 低エネルギー&長波長の極限で発散する。と教科書に書いてあるが、(9.56)〜(9.58) を教科書の形ではなく スピントロニクス理論の基礎/9-2#wea394fb のようにまとめれば、 その極限で単純に発散するわけでないようにも思える。
というのは置いておいて、(9.58) を変形すると以下の様に表せる。
(9.61)
&math( \rho_\phi^{(D)}(\bm r,t) &= \frac{\textcolor{red}{e^2}\nu(0)}{\textcolor{red}{a^3}} \sum_{\bm q} \int_{-\infty}^\infty \frac{d\Omega}{2\pi} \, e^{-i\bm q\cdot \bm r} e^{i\Omega t}\ \phi(\bm q,\Omega) \frac{ i\frac{\Omega}{Dq^2} }{1+i\frac{\Omega}{Dq^2}} \\&= \frac{\textcolor{red}{e^2}\nu(0)}{\textcolor{red}{a^3}} \sum_{\bm q} \int_{-\infty}^\infty \frac{d\Omega}{2\pi} \, e^{-i\bm q\cdot \bm r} e^{i\Omega t}\ \phi(\bm q,\Omega) \frac{ i\Omega }{Dq^2+i\Omega} \\&= \frac{\textcolor{red}{e^2}\nu(0)}{\textcolor{red}{a^3}} \int d^3r'\int_{-\infty}^\infty dt' \underbrace{\delta(\bm r'-\bm r)\delta(t'-t)}\, \left[ \sum_{\bm q} \int_{-\infty}^\infty \frac{d\Omega}{2\pi} \, e^{-i\bm q\cdot \bm r'} e^{i\Omega t'} \ \phi(\bm q,\Omega) \frac{ i\Omega }{Dq^2+i\Omega} \right] \\&= \frac{\textcolor{red}{e^2}\nu(0)}{\textcolor{red}{a^3}} \int d^3r'\int_{-\infty}^\infty dt' \overbrace{ \frac{1}{V}\sum_{\bm q'}e^{-i\bm q'\cdot(\bm r-\bm r')} \int\frac{d\Omega'}{2\pi}e^{i\Omega'(t-t')} } \ \left[ \sum_{\bm q} \int_{-\infty}^\infty \frac{d\Omega}{2\pi} \, e^{-i\bm q\cdot \bm r'} e^{i\Omega t'} \ \phi(\bm q,\Omega) \frac{ i\Omega }{Dq^2+i\Omega} \right] \\&= \frac{\textcolor{red}{e^2}\nu(0)}{\textcolor{red}{a^3V}} \sum_{\bm q'}e^{-i\bm q'\cdot\bm r} \int\frac{d\Omega'}{2\pi}e^{i\Omega't} \ \Bigg[ \sum_{\bm q} \int_{-\infty}^\infty \frac{d\Omega}{2\pi} \, \underbrace{\int d^3r'e^{-i(\bm q-\bm q')\cdot \bm r'}}_{\propto \, \delta_{\bm q\bm q'}} \ \underbrace{\int_{-\infty}^\infty dt' e^{i(\Omega-\Omega') t'} }_{\propto \, \delta(\Omega-\Omega')} \ \phi(\bm q,\Omega) \underbrace{\frac{ i\Omega }{Dq^2+i\Omega}} \Bigg] \\&= \frac{\textcolor{red}{e^2}\nu(0)}{\textcolor{red}{a^3V}} \int d^3r'\int_{-\infty}^\infty dt' \left[ \sum_{\bm q'} \int\frac{d\Omega'}{2\pi} \, e^{-i\bm q'\cdot(\bm r-\bm r')}e^{i\Omega'(t-t')} \overbrace{\frac{ 1 }{Dq'^2+i\Omega'}\phantom{\bigg|}} \right] \left[ \sum_{\bm q} \int_{-\infty}^\infty \frac{d\Omega}{2\pi} \, e^{-i\bm q\cdot \bm r'} i\Omega e^{i\Omega t'} \ \phi(\bm q,\Omega) \right] \\&= \frac{\textcolor{red}{e^2}\nu(0)}{\textcolor{red}{a^3V}} \int d^3r'\int_{-\infty}^\infty dt' D(\bm r-\bm r',t-t')\frac{\PD}{\PD t'}\phi(\bm r',t') );
教科書で負号が現われる理由、 の因子が落ちている理由が分からない?
ここで、
(9.62)
&math( D(\bm r,t) \equiv \sum_{\bm q} \int\frac{d\Omega}{2\pi} \, e^{-i\bm q\cdot\bm r}e^{i\Omega t} \frac{ 1 }{Dq^2+i\Omega} );
であり、 は
(9.64)
&math( \left(\frac{\PD}{\PD t}-D\nabla^2\right)D(\bm r,t) &= \sum_{\bm q} \int\frac{d\Omega}{2\pi} \, e^{-i\bm q\cdot\bm r} e^{i\Omega t} \frac{ Dq^2 + i\Omega }{Dq^2+i\Omega} \\&= \sum_{\bm q} \int\frac{d\Omega}{2\pi} \, e^{-i\bm q\cdot\bm r} e^{i\Omega t} \\&= \textcolor{red}{V}\delta(\bm r-\bm r')\delta(t-t') );
で定義される拡散型の微分方程式の Green 関数になっている。 (もしかして、(9.61) で抜けていた の因子は の定義に入れておくべき物なのか?つまり、(9.61) ではなく (9.62) に が抜けているのかも)
したがって、
(9.65)
&math( \left(\frac{\PD}{\PD t}-D\nabla^2\right)\rho_\phi^{(D)}(\bm r,t) &= \frac{\textcolor{red}{e^2}\nu(0)}{\textcolor{red}{a^3V}} \int d^3r'\int_{-\infty}^\infty dt' \left(\frac{\PD}{\PD t}-D\nabla^2\right) D(\bm r-\bm r',t-t')\frac{\PD}{\PD t'}\phi(\bm r',t') \\&= \frac{\textcolor{red}{e^2}\nu(0)}{\textcolor{red}{a^3V}} \int d^3r'\int_{-\infty}^\infty dt' \delta(\bm r-\bm r')\delta(t-t')\frac{\PD}{\PD t'}\phi(\bm r',t') \\&= \frac{\textcolor{red}{e^2}\nu(0)}{\textcolor{red}{a^3V}}\frac{\PD}{\PD t}\phi(\bm r,t) );
となって、電荷密度 はポテンシャルの時間変化により生み出されている。 (変化の方法を考えると、教科書の負号が正しい気がする。どこで間違えたんだろう???)
(9.66) は、
&math( D(\bm r-\bm r') &\equiv \sum_{\bm q} \, e^{-i\bm q\cdot(\bm r-\bm r')} \frac{1}{Dq^2} \\&= \frac{V}{D}\int\frac{d^3q}{(2\pi)^3} \, e^{-i\bm q\cdot(\bm r-\bm r')} \frac{1}{\,q^2\,} \\&= \frac{V}{8\pi^3D}\int dq \int qd\theta \int q\sin\theta d\varphi \, e^{-iq|\bm r-\bm r'|\cos \theta} \frac{1}{\,q^2\,} \\&= \frac{V}{4\pi^2D}\int dq \int d\theta \sin\theta \, e^{-iq|\bm r-\bm r'|\cos \theta} \\&= \frac{V}{4\pi^2D}\int dq \int d\theta \, \frac{1}{iq|\bm r-\bm r'|} \big( iq|\bm r-\bm r'|\sin \theta \big) e^{-iq|\bm r-\bm r'|\cos \theta} \\&= \frac{V}{4\pi^2D}\int dq \, \frac{1}{iq|\bm r-\bm r'|} \Big[e^{-iq|\bm r-\bm r'|\cos \theta}\Big]_0^{\pi} \\&= \frac{V}{2\pi^2D}\int dq \, \frac{1}{q|\bm r-\bm r'|} \frac{e^{iq|\bm r-\bm r'|}-e^{-iq|\bm r-\bm r'|}}{2i} \\&= \frac{V}{2\pi^2D}\int dq \, \frac{\sin \big(q|\bm r-\bm r'|\big)}{q|\bm r-\bm r'|} \\&= \frac{V}{2\pi^2D}\int \frac{dq'}{|\bm r-\bm r'|} \, \frac{\sin q'}{q'} \\&= \frac{V}{2\pi^2D} \frac{\pi}{|\bm r-\bm r'|} \, \\&= \frac{\red{V}}{\red{2}\pi D|\bm r-\bm r'|} );
ではないだろうか?
(やはり
に
の因子が足りないのかもしれない)
(以下勉強中)