tecktom.com

密度行列を2本の経路で運ぶ

2026年10月17日11分で読めますAI生成教科書『虎の巻』第4部 第6章 経路積分との合流:Feynman–Vernon の影響汎関数と、密度行列の構成の限界
目次
  1. 伝播関数を2つ使う
  2. 重心と差の変数で書く
  3. 木下の層

AI 生成・木下が未勉強

量子力学の様々な構成方法 第4部 第6章 §4.6.1 ★★★ ρ^(t)=U^ρ^U^†\hat\rho(t)=\hat U\hat\rho\hat U^\dagger の位置表示は伝播関数 KK と K∗K^* の積で運ばれ、経路積分では前向きの経路 xx と後ろ向きの経路 x′x' の2本の積分 ∫Dx Dx′ ei(S[x]−S[x′])/ℏ\int\mathcal Dx\,\mathcal Dx'\,\mathrm{e}^{\mathrm{i}(S[x]-S[x'])/\hbar} になる

第4部第5章(第4部第5章)は、環境と触れている系の時間発展を CPTP 写像として書き、環境が過去を覚えていない Markov な場合に Lindblad 方程式を導いた。その章の橋の節は、環境の相関時間が長い場合には Lindblad 方程式が立てられず、別の書き方が要ると述べて終わっている。この章はその書き方の1つ、Feynman と Vernon の影響汎関数を扱う。部の型(出発点の要請 → 基本の対象 → 基準との同値 → 共通例題 → 得意な問題 → 限界)の中では、同値性の地図の「経路積分 × 密度行列」の橋と、第4部全体の「限界」を受け持つ。

出発点に新しい公理は置かない。使うのは、第4部第1章(第4部第1章)の時間発展の要請 (D3) ρ^↦U^ρ^U^†\hat\rho\mapsto\hat U\hat\rho\hat U^\dagger と部分トレースの要請 (D5)、それに第3部第2章(第3部第2章)で基準から導いた伝播関数の経路積分表示だけである。U^\hat U は左から、U^†\hat U^\dagger は右から密度行列に掛かるので、密度行列の時間発展は「前向きの経路」と「後ろ向きの経路」の2本の経路積分になる(第4.6.1節)。系と環境の全体をこの形に書いて環境の経路を先に積分すると、系の2本の経路だけの汎関数が1つ残る。これが影響汎関数で、環境が系の2本の経路をどれだけ区別して記録したかを表す(第4.6.2節)。環境が調和振動子の集まりで系と線形に結合している場合(Caldeira–Leggett 模型)、影響汎関数は厳密に計算でき、摩擦の核と雑音の核の2つで決まる。高温でこれが時間について局所的になると、Caldeira–Leggett の主方程式が出てくる(第4.6.3節)。共通例題①の振動子と共通例題②の2準位系は、どちらも影響汎関数で厳密に解くことができる。振動子の定常状態は量子の揺動散逸定理で決まり、結合が弱い極限で第4部第4章の熱平衡に戻る。2準位系の位相緩和は、温度 0 で指数関数でなく冪で減衰する(第4.6.4節)。第4.6.5節では、Caldeira–Leggett の主方程式が Lindblad 形でないために正値性を壊すこと、強い結合の熱平衡が系の Gibbs 状態でないこと、初めの相関と記憶のある環境が写像の半群の外にあることを順に示す。橋の節(第4.6.6節)では、影響汎関数を基準・密度行列・古典の Langevin 方程式・位相空間の Kramers 方程式の言葉に書き直して並べる。

この章では、系の座標は共通例題の約束どおり xx と書き、2本の経路 x(s)x(s)、x′(s)x'(s) から重心と差

X(s)=x(s)+x′(s)2,ξ(s)=x(s)−x′(s),Σ(s)=x(s)+x′(s)=2X(s)(4.6.1)X(s)=\frac{x(s)+x'(s)}{2},\qquad \xi(s)=x(s)-x'(s),\qquad \Sigma(s)=x(s)+x'(s)=2X(s)\tag{4.6.1}

を作る(XX・ξ\xi は第4部第2章の混合した Gauss 波束と同じ記号)。環境の振動子の座標は qkq_k、質量は mkm_k、角周波数は ωk\omega_k と書く。恒等演算子は 1\mathbf 1、反交換子は {A^,B^}+\{\hat A,\hat B\}_{+} である。

伝播関数を2つ使う

孤立した系の密度演算子は (D3) により ρ^(t)=U^(t)ρ^(0)U^(t)†\hat\rho(t)=\hat U(t)\hat\rho(0)\hat U(t)^\dagger、U^(t)=e−iH^t/ℏ\hat U(t)=\mathrm{e}^{-\mathrm{i}\hat Ht/\hbar} で発展する。位置表示の行列要素 ρ(xf,xf′,t)=⟨xf∣ρ^(t)∣xf′⟩\rho(x_f,x_f',t)=\langle x_f | \hat\rho(t) | x_f' \rangle を取り、ρ^(0)\hat\rho(0) の両側に位置の完全系を挟むと

ρ(xf,xf′,t)=∫dxi dxi′  K(xf,t;xi,0) ρ(xi,xi′,0) K(xf′,t;xi′,0)∗(4.6.2)\rho(x_f,x_f',t)=\int\mathrm{d} x_i\,\mathrm{d} x_i'\;K(x_f,t;x_i,0)\,\rho(x_i,x_i',0)\,K(x_f',t;x_i',0)^*\tag{4.6.2}

になる。K(xf,t;xi,0)=⟨xf∣U^(t)∣xi⟩K(x_f,t;x_i,0)=\langle x_f | \hat U(t) | x_i \rangle は伝播関数で、右側の ⟨xi′∣U^†∣xf′⟩=⟨xf′∣U^∣xi′⟩∗\langle x_i' | \hat U^\dagger | x_f' \rangle =\langle x_f' | \hat U | x_i' \rangle ^* から複素共役が出る。第3部第2章で示したとおり KK は経路積分 ∫Dx eiS[x]/ℏ\int\mathcal Dx\,\mathrm{e}^{\mathrm{i} S[x]/\hbar} に書けるので、(4.6.2) は

ρ(xf,xf′,t)=∫dxi dxi′∫x(0)=xix(t)=xf ⁣Dx∫x′(0)=xi′x′(t)=xf′ ⁣Dx′  ei(S[x]−S[x′])/ℏ ρ(xi,xi′,0)(4.6.3)\rho(x_f,x_f',t)=\int\mathrm{d} x_i\,\mathrm{d} x_i'\int_{x(0)=x_i}^{x(t)=x_f}\!\mathcal Dx\int_{x'(0)=x_i'}^{x'(t)=x_f'}\!\mathcal Dx'\; \mathrm{e}^{\mathrm{i}(S[x]-S[x'])/\hbar}\,\rho(x_i,x_i',0)\tag{4.6.3}

という2本の経路の積分になる。xx は U^\hat U から来る前向きの経路、x′x' は U^†\hat U^\dagger から来る経路で、作用の前の符号が逆になっている。x′x' を「時刻 tt から 00 へ戻る経路」と読めば、2本は時刻 tt でつながり、時刻 00 の密度行列を挟んで閉じた1本の時間の道(図 4.6.1)を往復する。この読み方は Schwinger (1961), J. Math. Phys. 2, 407 が量子振動子の Brown 運動で用い、Keldysh (1964), Zh. Eksp. Teor. Fiz. 47, 1515(英訳 Sov. Phys. JETP 20, 1018 (1965))が非平衡の摂動論の図形の規則に組み上げた。この章では名前だけ挙げ、その図形の規則には進まない。

図

図 4.6.1 密度行列の時間発展 (4.6.3) の2本の経路。上の経路は U^\hat U から、下の経路は U^†\hat U^\dagger から来て、作用の符号が逆になる。2本を時刻 tt でつなぐと閉じた時間の道になる。環境を積分して消すと、2本の経路を同じ時刻で結ぶ因子 F[x,x′]F[x,x'] が残る(破線)。

重心と差の変数で書く

2本の経路を重心 XX と差 ξ\xi に書き直すと、作用の差が見通しのよい形になる。x=X+ξ/2x=X+\xi/2、x′=X−ξ/2x'=X-\xi/2 を Lagrangian L=m2x˙2−V(x)L=\tfrac m2\dot x^2-V(x) に入れると

L(x,x˙)−L(x′,x˙′)=mX˙ξ˙−[V(X+ξ2)−V(X−ξ2)]=mX˙ξ˙−V′(X) ξ−V′′′(X)24ξ3−⋯(4.6.4)L(x,\dot x)-L(x',\dot x')=m\dot X\dot\xi-\Big[V\big(X+\tfrac\xi2\big)-V\big(X-\tfrac\xi2\big)\Big] =m\dot X\dot\xi-V'(X)\,\xi-\frac{V'''(X)}{24}\xi^3-\cdots\tag{4.6.4}

になる。運動エネルギーの差は m2[(X˙+ξ˙/2)2−(X˙−ξ˙/2)2]=mX˙ξ˙\tfrac m2[(\dot X+\dot\xi/2)^2-(\dot X-\dot\xi/2)^2]=m\dot X\dot\xi で、ちょうど ξ\xi について1次である。ポテンシャルの差は ξ\xi の奇数次だけを含み、VV が2次以下なら ξ\xi について厳密に1次になる。この「2次の系では作用の差が ξ\xi の1次で尽きる」ことが、第4.6.4節で振動子を厳密に解くときの要になる。

自由粒子では、古典の経路の作用 Scl=m(xf−xi)2/(2t)S_{\mathrm{cl}}=m(x_f-x_i)^2/(2t)(第3部第3章、第3部第3章)の差が

Scl(xf,xi)−Scl(xf′,xi′)=mt (Xf−Xi)(ξf−ξi)(4.6.5)S_{\mathrm{cl}}(x_f,x_i)-S_{\mathrm{cl}}(x_f',x_i')=\frac mt\,(X_f-X_i)(\xi_f-\xi_i)\tag{4.6.5}

になり、KK∗K K^* の前因子は ∣m/(2πiℏt)∣=m/(2πℏt)\lvert m/(2\pi\mathrm{i}\hbar t)\rvert =m/(2\pi\hbar t) である。対角要素 ξf=0\xi_f=0 では指数が ξi\xi_i について1次なので、ξi\xi_i の積分は ρ(⋅,⋅,0)\rho(\cdot,\cdot,0) の ξi\xi_i についての Fourier 変換になる。初めの状態の ξi\xi_i 依存は運動量の分布を表すから、自由粒子では「各点から運動量の分布どおりに粒子が飛び出す」という古典の像が、KK∗K K^* の積分から厳密に出てくる。

例題 4.6.1 混合した Gauss 波束を KK∗K K^* で運ぶ ★★

オリジナル第4部第2章(第4部第2章)の混合した Gauss 波束ρ(x,x′,0)=(2πσ02)−1/2exp⁡ ⁣(−X2/(2σ02)−Δp2ξ2/(2ℏ2))\rho(x,x',0)=(2\pi\sigma_0^2)^{-1/2}\exp\!\big(-X^2/(2\sigma_0^2)-\Delta p^2\xi^2/(2\hbar^2)\big) を、自由粒子の (4.6.2) で時刻 tt まで運び、位置の分散 σ(t)2\sigma(t)^2 を求めよ。Δp=ℏ/(2σ0)\Delta p=\hbar/(2\sigma_0) の純粋な場合に共通例題③の正答 σ02+(ℏt/(2mσ0))2\sigma_0^2+(\hbar t/(2m\sigma_0))^2 と一致することを確かめ、87^{87}Rb、σ0=1 μm\sigma_0=1\,\mu\mathrm m、t=1 mst=1\,\mathrm{ms} で、純粋な場合と温度 1 μK1\,\mu\mathrm K の混合(Δp=mkBT\Delta p=\sqrt{m k_{\mathrm B} T})の場合の σ(t)\sigma(t) を数値で出せ。

解答

対角要素 ξf=0\xi_f=0 を取る。(4.6.5) から指数は −im(Xf−Xi)ξi/(ℏt)-\mathrm{i} m(X_f-X_i)\xi_i/(\hbar t) で、

ρ(Xf,Xf,t)=m2πℏt∫dXi e−Xi2/(2σ02)2πσ02∫dξi exp⁡ ⁣(−im(Xf−Xi)ξiℏt−Δp2ξi22ℏ2)\rho(X_f,X_f,t)=\frac{m}{2\pi\hbar t}\int\mathrm{d} X_i\,\frac{\mathrm{e}^{-X_i^2/(2\sigma_0^2)}}{\sqrt{2\pi\sigma_0^2}}\int\mathrm{d}\xi_i\,\exp\!\Big(-\frac{\mathrm{i} m(X_f-X_i)\xi_i}{\hbar t}-\frac{\Delta p^2\xi_i^2}{2\hbar^2}\Big)

になる(dxi dxi′=dXi dξi\mathrm{d} x_i\,\mathrm{d} x_i'=\mathrm{d} X_i\,\mathrm{d}\xi_i)。ξi\xi_i の Gauss 積分は 2π (ℏ/Δp)exp⁡ ⁣(−m2(Xf−Xi)2/(2t2Δp2))\sqrt{2\pi}\,(\hbar/\Delta p)\exp\!\big(-m^2(X_f-X_i)^2/(2t^2\Delta p^2)\big) で、前因子とまとめると

ρ(Xf,Xf,t)=∫dXi e−Xi2/(2σ02)2πσ02⋅e−(Xf−Xi)2/(2(Δp t/m)2)2π(Δp t/m)2\rho(X_f,X_f,t)=\int\mathrm{d} X_i\,\frac{\mathrm{e}^{-X_i^2/(2\sigma_0^2)}}{\sqrt{2\pi\sigma_0^2}}\cdot\frac{\mathrm{e}^{-(X_f-X_i)^2/(2(\Delta p\,t/m)^2)}}{\sqrt{2\pi(\Delta p\,t/m)^2}}

となる。これは分散 σ02\sigma_0^2 の Gauss 分布と分散 (Δp t/m)2(\Delta p\,t/m)^2 の Gauss 分布の畳み込みなので、

σ(t)2=σ02+(Δp tm)2\sigma(t)^2=\sigma_0^2+\Big(\frac{\Delta p\,t}{m}\Big)^2

で、第4部第2章の結果と一致する。Δp=ℏ/(2σ0)\Delta p=\hbar/(2\sigma_0) を入れると σ02+(ℏt/(2mσ0))2\sigma_0^2+(\hbar t/(2m\sigma_0))^2 で、共通例題③の正答の形である。数値は、m=1.443×10−25 kgm=1.443\times10^{-25}\,\mathrm{kg} で純粋な場合 ℏt/(2mσ0)=0.365 μm\hbar t/(2m\sigma_0)=0.365\,\mu\mathrm m、σ(1 ms)=1.065 μm\sigma(1\,\mathrm{ms})=1.065\,\mu\mathrm m。1 μK1\,\mu\mathrm K の混合では Δp t/m=kBT/m t=9.79 μm\Delta p\,t/m=\sqrt{ k_{\mathrm B} T/m}\,t=9.79\,\mu\mathrm m、σ=9.84 μm\sigma=9.84\,\mu\mathrm m になる。混合の場合の ξ\xi の幅 ℏ/Δp=0.0747 μm\hbar/\Delta p=0.0747\,\mu\mathrm m は位置の幅 σ0\sigma_0 よりずっと小さく、ξ\xi の方向に細い密度行列ほど速く広がる。

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『虎の巻』の記事一覧へ