スピントロニクス理論の基礎/8-1 の履歴(No.18)
更新>>>> [次の章へ]
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 \frac{e^{-\beta H(t_0)}}{\red{Z_0}}O(t_0) \bigr\rangle );
は
を含まず
は
を含む和になっている。
できる限り をあらわに書いてみた。
は分配関数。
この Σ の部分を trace と呼んでいるのは、 を基底として を行列に表した際の 要素が であり、 Σ の部分はちょうど対角要素の和になっているためである。
線形代数で学ぶように、trace は のような ユニタリ変換による基底の付け替えに対して値が変化しない量であり、 基底の選び方に依らず値の決まる物理量である。
がエネルギー固有関数の時は、
(8.2)
は が で持つエネルギーとして定義される。
このとき、
&math( &\bra{\alpha(t_0)}e^{-\beta H(t_0)} =\big[(e^{-\beta H(t_0)})^\dagger\ket{\alpha(t_0)}\big]^\dagger =\big[e^{-\beta H^\dagger(t_0)}\ket{\alpha(t_0)}\big]^\dagger\\ &=\big[e^{-\beta H(t_0)}\ket{\alpha(t_0)}\big]^\dagger =\big[e^{-\beta E_\alpha}\ket{\alpha(t_0)}\big]^\dagger =e^{-\beta E_\alpha}\bra{\alpha(t_0)} );
より、
(8.1A)
trace を取ることが量子力学で期待値を求めることに他ならず(量子平均)
そこに
の因子を含めることで同時に統計平均も取れることになる。
がエネルギーの固有関数でないときには、(8.1) 式は正しいのだろうか? エネルギー固有関数でなければ、一般に
となるような気がするのだが、どちらが目的にかなうかが明確ではない。
恐らく、 はエネルギーの固有関数でなければ (8.1) は意味をなさないのではないか?言い換えると、(8.1) ではなく (8.1A) が基本式なのではないか、と思える。
(2015-02-20 追記) 上でも述べているとおり、トレースの値は任意のユニタリー変換に対して保存するから、異なる物理量に対する固有ベクトルで書き直してもトレースの値が変わることはない。(8.1) は基底の取り方によらず意味を持つ(のかもしれない・・・)。
時間発展演算子の積分表示†
さて、 はハミルトニアンの固有関数なので、 ハミルトニアンが時間的に変化しない限り、 も時間に対して位相以外変化しない。
(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( U(t,t_0)\\ &=\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\left(\frac{-i}{\hbar}\right)^n H^n \left[\frac{1}{n!}(t-t_0)^n\right]\\ =&\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) より
したがって、
そこで時刻 と の間を に分割すれば、
(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)} );
を得る。
一方で、この式のエルミート共役を取ってみる。これは個々の行列のエルミート共役を取り、 行列の積を入れ替え、 を に付け替えることで得られる。 ハミルトニアンがエルミート行列 であることを用いると、次式が得られる。
(8.8A)
&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.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)
は一目瞭然であるし、 と途中の任意の時刻 ただし で分割して、
&math( U(t, t_0) &= \underbrace{ 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_{k})(t_{k+1}-t_{k})} }_{U(t,t_k)} \underbrace{ e^{\frac{-i}{\hbar}H(t_{k-1})(t_{k}-t_{k-1})} \cdots e^{\frac{-i}{\hbar}H(t_1)(t_2-t_1)} e^{\frac{-i}{\hbar}H(t_0)(t_1-t_0)} }_{U(t_{k},t_0)} \\&= U(t,t_k)U(t_k,t_0) );
と表せることも容易に理解できる。 (物理的には から までの時間発展が から までの時間発展と から までの時間発展と を続けて行うことにより得られるのは自然であるが、 教科書の (8.7) の式を見ていてもなかなか理解できないのではないか)
注)これらの表記では として、 個々の が無限小の時間間隔を表していると考えている。
(8.7A) の指数関数の積の部分を指数の和に、さらに積分に直すことで、
&math(U(t, t_0) &= e^{\frac{-i}{\hbar}H(t_n)(t-t_n)+\frac{-i}{\hbar}H(t_{n-1})(t_n-t_{n-1})+ \cdots +\frac{-i}{\hbar}H(t_1)(t_2-t_1)+\frac{-i}{\hbar}H(t_0)(t_1-t_0)}\\ &= e^{\frac{-i}{\hbar}[H(t_n)(t-t_n)+H(t_{n-1})(t_n-t_{n-1})+\cdots +H(t_1)(t_2-t_1)+H(t_0)(t_1-t_0)]}\\ &= e^{\frac{-i}{\hbar}[H(t_n)(t-t_n)+\sum_{k=1}^n H(t_{k-1})(t_k-t_{k-1})]}\\ &= e^{\frac{-i}{\hbar}T\int_{t_0}^t dt'H(t')}\\ &= Te^{\frac{-i}{\hbar}\int_{t_0}^t dt'H(t')}\\ );
のようにして (8.7) の最終式を理解できる。
この表記では、 や は のような項をどちらから順に並べるかを表す演算子であると自然に理解できる。
(8.48) などで が ばかりでなく も含めて並べ替える話になるが、このときには無数に並ぶ の間に が挟まる形を思い浮かべれば良い。
(8.10)
上記より、 であり、 がユニタリ演算子であることが証明された。
疑問点†
以下 2012/1/30 のセミナーで出た話の内容を追記:
がユニタリ演算子であることは、 シュレーディンガー方程式が時間反転に対して対称であることに由来している。
不純物散乱などを入れた結果、 系が時間反転に対して対称でなくなった場合におかしなことが起きないかどうか、 注意深く考えなければならない。
(2013-07-01 追記)→ ここで展開している理論では、 統計力学で言うところの微視的状態のそれぞれは時間反転に対して対称な、 純粋にシュレーディンガー方程式に沿って時間発展を行う波動関数であり、 量子統計平均後の物理量がその熱力学的な時間の矢に従って時間反転に対して非対称な時間発展を行う。
そのあたり、注意して読み進めるべき。
物理量の時間発展†
(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.5A)
さらに、(8.9) で得た を使った。そして、
(8.16)
は、Heisenberg 表示でのハミルトニアンである。
ということで、
(8.15)
が正しいような・・・
ただしこの場合の は演算子が陽に時刻に依存していない場合にはゼロになるため、次章で扱うように であったり、 であったりする場合には で、結局 (8.15) と同じ形になる。
(8.16) の脚注:
のように時間に依存しない場合には、(8.9B) でも見たように
(8.7B)
&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.7A) からも
(8.7C)
&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)} );
として確認できる。
質問・コメント†
私の勘違いでした。†
東北大CMPT M2 ()
下の質問は取り消させてください。
私の勘違いでした。
http://park.itc.u-tokyo.ac.jp/kato-yusuke-lab/nagai/note_070106_in.pdf
を読んで納得しました。
結局(8.7)は不適切な表現ですよね。。†
東北大CMPT M2 ()
(8.7)についてですが、結局、Hが時間発展してしまうとn!の逆数は現れないわけですから、Texp[\int^t_(t_0) dt' H(t')]と書くのは間違いだと思うのですが、どうでしょうか。