スピントロニクス理論の基礎/8-1 の履歴(No.4)
更新培風館 多々良源 「スピントロニクス理論の基礎」 を勉強していく上で、気付いた誤植や難しい部分の注釈を書いておこうと思います。
8-1 物理量†
(8.1)
&math(\overline O(t_0) \equiv \frac{1}{Z_0} \trace[e^{-\beta H(t_0)}O(t_0)] = \frac{1}{Z_0} \sum_{\alpha(t_0)} \braket{\alpha(t_0)|e^{-\beta H(t_0)}O(t_0)|\alpha(t_0)} \equiv \bigl\langle e^{-\beta H(t_0)}O(t_0) \bigr\rangle \equiv \bigllangle O(t_0) \bigrrangle);
は
を含まない
は
を含む
できる限り をあらわに書いてみた。
は分配関数。
この Σ の部分を trace と呼んでいるのは、 を基底として を行列に表した際の 要素が であり、 Σ の部分はちょうど対角要素の和になっているためである。
線形代数で学ぶように、trace は のような ユニタリ変換による基底の付け替えに対して値が変化しない量であり、 基底の選び方に依らず値の決まる物理量である。
trace を取ることが量子力学で期待値を求めることに他ならず(量子平均)
そこに
の因子を含めることで同時に統計平均も取れることになる。
(8.2)
は が で持つエネルギーとして定義される。
はハミルトニアンの固有関数なので、 ハミルトニアンが時間的に変化しない限り、 も時間に対して位相以外変化しない。
(8.3)
において が変化すると、この式に沿って量子状態も時間変化を始める。
(8.4)
このようにして定義した がユニタリ演算子となることが (8.10) で示される。
(8.5)
(8.3) より
(8.6)
両辺を積分して を使う。
(8.7)
右辺の を分解していくと、
&math(&U(t,t_0)=1-\frac{i}{\hbar}\int_{t_0}^t dt_1 H(t_1) \left[ 1-\frac{i}{\hbar}\int_{t_0}^{t_1} dt_2 H(t_2) U(t_2,t_0) \right]\\ &=1+\frac{-i}{\hbar}\int_{t_0}^t dt_1 H(t_1)+\left(\frac{-i}{\hbar}\right)\int_{t_0}^t dt_1 \int_{t_0}^{t_1} dt_2 H(t_1)H(t_2) U(t_2,t_0)\\ &=1+\frac{-i}{\hbar}\int_{t_0}^t dt_1 H(t_1)
- \left(\frac{-i}{\hbar}\right)^2\int_{t_0}^t dt_1 \int_{t_0}^{t_1} dt_2 H(t_1)H(t_2)\\ &+\left(\frac{-i}{\hbar}\right)^3\int_{t_0}^t dt_1 \int_{t_0}^{t_1} dt_2 \int_{t_0}^{t_2} dt_3 H(t_1)H(t_2)H(t_3) + \cdots\\ &=\sum_{n=0}^\infty\left(\frac{-i}{\hbar}\right)^n\int_{t_0}^t dt_1 \int_{t_0}^{t_1} dt_2 \cdots \int_{t_0}^{t_{n-1}} dt_n H(t_1)H(t_2)\cdots H(t_n)\\ &\equiv Te^{\frac{-i}{\hbar}\int_{t_0}^tdt'H(t')} );
この式から が 「時間順序づけ演算子」 だと言われてもちょっとピンと来ない。 確かに積分の中では が時間の古い方が右側になるように並んでいるけれど・・・
この演算子について、この教科書では示されていないが (8.7A) や (8.8A) の形で理解すべきな気がしている。
(8.9) は唐突なので、以下ちょっと話の流れを変える。
(8.11)
&math(&\left[U(t,t_0)\right]^\dagger\\ &=\sum_{n=0}^\infty\left(\frac{i}{\hbar}\right)^n\int_{t_0}^t dt_1 \int_{t_0}^{t_1} dt_2 \cdots \int_{t_0}^{t_{n-1}} dt_n H(t_n)H(t_{n-1})\cdots H(t_1)\\ &\equiv \overline Te^{\frac{-i}{\hbar}\int_{t_0}^tdt'H(t')});
を考える。
のように、エルミート共役により演算子のかけ算の順が入れ替わること、 はエルミート演算子であること 、エルミート共役で虚数単位の符号が反転すること から、 のエルミート共役が上記のように書けることが分かる。
これと (8.7) との積を作ってみると、
(8.9)
&math(&\left[U(t,t_0)\right]^\dagger U(t,t_0)\\ &=\overline Te^{\frac{i}{\hbar}\int_{\textcolor{red}{t_0}}^tdt'H(t')}\, Te^{\frac{-i}{\hbar}\int_{\textcolor{red}{t_0}}^tdt'H(t')}\\ &=\left(1+\frac{i}{\hbar}\int_{t_0}^tdt_1H(t_1)+\left(\frac{i}{\hbar}\right)^2\int_{\textcolor{red}{t_0}}^t dt_1\int_{\textcolor{red}{t_0}}^{t_1}dt_2H(t_2)H(t_1)+\dots\right)\times\\ &\left(1+\frac{-i}{\hbar}\int_{t_0}^tdt_1H(t_1)+\left(\frac{-i}{\hbar}\right)^2\int_{\textcolor{red}{t_0}}^tdt_1\int_{\textcolor{red}{t_0}}^{t_1}dt_2H(t_1)H(t_2)+\dots\right)\\ &=1);
実際に最後の等号を示そうとすると大変であるが、ちゃんと計算すれば示せないことはない。
様子を見るため、 の2次の項を抜き出してみると、
(8.9A)
&math(
&\int_{t_0}^tdt_1\int_{t_0}^{t_1}dt_2 H(t_2)H(t_1)
- \int_{t_0}^tdt_1\int_{t_0}^{t_1}dt_2 H(t_1)H(t_2)
- \int_{t_0}^tdt_1\int_{t_0}^{t}dt_2 H(t_1)H(t_2)\\ =&\int_{t_0}^tdt_2\int_{t_0}^{t_2}dt_1 H(t_1)H(t_2)
- \int_{t_0}^tdt_1\int_{t_0}^{t_1}dt_2 H(t_1)H(t_2)
- \int_{t_0}^tdt_1\int_{t_0}^{t}dt_2 H(t_1)H(t_2)\\ =&\int_{t_0}^tdt_2\int_{t_0}^tdt_1 \theta(t_2-t_1)H(t_1)H(t_2)
- \int_{t_0}^tdt_1\int_{t_0}^tdt_2 \theta(t_1-t_2)H(t_1)H(t_2)
- \int_{t_0}^tdt_1\int_{t_0}^{t}dt_2 H(t_1)H(t_2)\\ =&\int_{t_0}^tdt_2\int_{t_0}^tdt_1 [\theta(t_2-t_1)+\theta(t_1-t_2)-1]H(t_1)H(t_2) = 0 );
となる。
- 2行目では第1項の と の名前を付け替えた
- 3行目では積分範囲を
に拡大し、
もとの範囲を
関数で表した
- 4行目では を用いた
少し脱線するが、(8.7) に現れる Σ が指数関数の展開を表しているにも関わらず、 の因子が現れない理由が上記の 関数の部分に 表された積分範囲にある。
これを見るために が時間に依存しない場合を計算しておく。
(8.9B)
&math( &\sum_{n=0}^\infty\left(\frac{-i}{\hbar}\right)^n\int_{t_0}^t dt_1 \int_{t_0}^{t_1} dt_2 \cdots \int_{t_0}^{t_{n-1}} dt_n H(t_1)H(t_2)\cdots H(t_n)\\ =&\sum_{n=0}^\infty\left(\frac{-i}{\hbar}\right)^n H^n \int_{t_0}^t dt_1 \int_{t_0}^{t_1} dt_2 \cdots \int_{t_0}^{t_{n-1}} dt_n 1\\ =&\sum_{n=0}^\infty\left(\frac{-i}{\hbar}\right)^n H^n \int_{t_0}^t dt_1 \int_{t_0}^{t_1} dt_2 \cdots \int_{t_0}^{t_{n-2}} dt_{n-1} (t_{n-1}-t_0)\\ =&\sum_{n=0}^\infty\left(\frac{-i}{\hbar}\right)^n H^n \int_{t_0}^t dt_1 \int_{t_0}^{t_1} dt_2 \cdots \int_{t_0}^{t_{n-3}} dt_{n-2} \left[\frac{1}{2}(t_{n-1}-t_0)^2\right]_{t_0}^{t_{n-2}}\\ =&\sum_{n=0}^\infty\left(\frac{-i}{\hbar}\right)^n H^n \int_{t_0}^t dt_1 \int_{t_0}^{t_1} dt_2 \cdots \int_{t_0}^{t_{n-4}} dt_{n-3} \left[\frac{1}{3!}(t_{n-2}-t_0)^3\right]_{t_0}^{t_{n-3}}\\ \vdots\\ =&\sum_{n=0}^\infty\left(\frac{-i}{\hbar}\right)^n H^n \left[\frac{1}{n!}(t_1-t_0)^n\right]_{t_0}^t\\ =&\sum_{n=0}^\infty\frac{1}{n!}\left(\frac{-i}{\hbar} H (t-t_0)\right)^n\\ =&e^{-\frac{i}{\hbar}H\cdot(t-t_0)} );
このように 番目の項の積分範囲は の になっている。 また、この形から分かるように、 が時間に依らないときには 自体が と可換になっている。
閑話休題で、(8.9) について: 次の項でも上と同様にゼロになることを示せるが、 この話は次のように考えることで、より直感的に理解できる。
(8.5) より
&math(U(t, t_0) = e^{\frac{-i}{\hbar}H(t_n)(t-t_n)} e^{\frac{-i}{\hbar}H(t_{n-1})(t_n-t_{n-1})} \cdots e^{\frac{-i}{\hbar}H(t_1)(t_2-t_1)} e^{\frac{-i}{\hbar}H(t_0)(t_1-t_0)} ); (8.7A)
一方で、
&math([U(t, t_0)]^\dagger = e^{\frac{i}{\hbar}H(t_0)(t_1-t_0)} e^{\frac{i}{\hbar}H(t_1)(t_2-t_1)} \cdots e^{\frac{i}{\hbar}H(t_{n-1})(t_n-t_{n-1})} e^{\frac{i}{\hbar}H(t_n)(t-t_n)}); (8.8A)
(8.7) や (8.8) の意味をこのように理解していると物理的な理解が深まる。
これらを用いれば、
&math(&[U(t, t_0)]^\dagger U(t, t_0) =\\ & e^{\frac{i}{\hbar}H(t_0)(t_1-t_0)} e^{\frac{i}{\hbar}H(t_1)(t_2-t_1)} \cdots e^{\frac{i}{\hbar}H(t_{n-1})(t_n-t_{n-1})} e^{\frac{i}{\hbar}H(t_n)(t-t_n)}\\ &e^{\frac{-i}{\hbar}H(t_n)(t-t_n)} e^{\frac{-i}{\hbar}H(t_{n-1})(t_n-t_{n-1})} \cdots e^{\frac{-i}{\hbar}H(t_1)(t_2-t_1)} e^{\frac{-i}{\hbar}H(t_0)(t_1-t_0)}\\ &=1); (8.9A)
は一目瞭然である。
注)これらの表記では として、 個々の が無限小の時間間隔を表していると考えている。
この表記では、 や は のような項をどちらから順に並べるかを表す演算子であると理解できる。
(8.48) などで が ばかりでなく も含めて並べ替える話になるが、このときには無数に並ぶ の間に が挟まる形を思い浮かべれば良い。
(8.10)
上記より、 であり、 がユニタリ演算子であることが証明された。
(8.12)
&math(\overline O(t)&=\frac{1}{Z_0}\sum_\alpha e^{-\beta E_\alpha} \braket{\alpha(t)|O(t)|\alpha(t)}\\ &=\frac{1}{Z_0}\textcolor{red}{\sum_\alpha} \braket{\alpha(t_0)|e^{-\beta H(t_0)}[U(t,t_0)]^\dagger O(t)U(t,t_0)|\alpha(t_0)} );
この式は前半部分の最重要の式である。 特に、 の部分に現れるエネルギーが時刻 での値ではなく の時の値であることに注目したい。
において系は熱平衡にあり、各エネルギー固有状態 に の確率で存在していた。量子力学的に見れば、それぞれの量子状態は系のハミルトニアンが変化しない限り安定で、お互いの間を行き来することはない。
その後ハミルトニアンが変化して系が熱平衡からずれると、 における固有状態 は (8.3) に従って時間発展する。 この時間発展が (8.4) のようにユニタリ変換で表せることから分かるとおり、 この時間発展はどんな に対してさえ決定論的であり、 また未来永劫、 が における別の固有状態 が時間発展した と混じり合うことはない。 これは &math(\braket{\alpha(t)|\alpha'(t)}=\braket{\alpha|U^\dagger U|\alpha'} =\braket{\alpha|\alpha'}=\delta_{\alpha,\alpha'});
したがって、時刻 における物理量 の統計平均は、 時刻 における状態 に対する の期待値 に、 において系が状態 にいた確率 を掛けて、 (8.12) のようにして求められるのである。
もちろん、系が時刻 において熱平衡にある必要は無いし、 ハミルトニアンの変化がどれだけ急峻であったとしても (8.12) は成り立つ。 ここまでの議論には断熱近似等は入っていない(と私は理解しています)。
(8.13)
(8.14)
&math(\overline O(t)&=\frac{1}{Z_0}\sum_\alpha \braket{\alpha(t_0)|e^{-\beta H(t_0)}O_H(t)|\alpha(t_0)}\\ &=\frac{1}{Z_0}\sum_\alpha e^{-\beta E_\alpha} \braket{\alpha(t_0)|O_H(t)|\alpha(t_0)} );
(8.15) 以下の議論は分かりにくい上に間違っている?
(8.17) は (8.13) そのものであるが、これを微分することにより、
(8.18)
&math(&\dot O_H \\ &= \dot U^\dagger O U\textcolor{red}{ + U^\dagger \dot O U} + U^\dagger O \dot U\\ &= \left[\frac{i}{\hbar}U^\dagger H\right] O U \textcolor{red}{+ \left(\frac{\PD O(t)}{\PD t}\right)_H} + U^\dagger O \left[\frac{-i}{\hbar}HU\right]\\ &=\frac{i}{\hbar}U^\dagger(HO-OH)U \textcolor{red}{+ \left(\frac{\PD O(t)}{\PD t}\right)_H}\\ &=\frac{i}{\hbar}(U^\dagger HU U^\dagger OU-U^\dagger OUU^\dagger HU) \textcolor{red}{+ \left(\frac{\PD O(t)}{\PD t}\right)_H}\\ &=\frac{i}{\hbar}(H_HO_H-O_HH_H) \textcolor{red}{+ \left(\frac{\PD O(t)}{\PD t}\right)_H}\\ &=\frac{i}{\hbar}\left[H_H(t), O_H(t)\right] \textcolor{red}{+ \left(\frac{\PD O(t)}{\PD t}\right)_H});
ここで、(8.5) およびそのエルミート共役である、
さらに、(8.9) で得た を使った。そして、
(8.16)
は、Heisenberg 表示でのハミルトニアンである。
ということで、
(8.15)
が正しいような・・・
ただしこの場合の は演算子が陽に時刻に依存していない場合 にはゼロになるため、次章で扱うように であったり、 であったりする場合には で、結局 (8.15) と同じ形になります。
(8.16) の脚注:
のように時間に依存しない場合には、
&math(&Te^{\textcolor{red}{\frac{-i}{\hbar}}\int_{\textcolor{red}{t_0}}^{\textcolor{red}{t}}dt'H(t')}\\ &=e^{\textcolor{red}{\frac{-i}{\hbar}}H\int_{\textcolor{red}{t_0}}^{\textcolor{red}{t}}dt'}\\ &=e^{\textcolor{red}{\frac{-i}{\hbar}}H(t-t_0)}); (8.7B)
となる。
これは、(8.7A) からも
&math(&U(t, t_0) = e^{\frac{-i}{\hbar}H(t_n)(t-t_n)} e^{\frac{-i}{\hbar}H(t_{n-1})(t_n-t_{n-1})} \cdots e^{\frac{-i}{\hbar}H(t_1)(t_2-t_1)} e^{\frac{-i}{\hbar}H(t_0)(t_1-t_0)}\\ &=e^{\frac{-i}{\hbar}H[(t-t_n)+(t_n-t_{n-1})+ \cdots
- (t_2-t_1)+(t_1-t_0)]}\\ &=e^{\frac{-i}{\hbar}H(t-t_0)} ); (8.7C)
として確認できる。