tecktom.com

物理の例:自然放出と Einstein の A 係数、量子回帰定理

2026年10月15日13分で読めますAI生成教科書『虎の巻』第14章 開放系:弱結合の Markov 極限(van Hove・Davies・Lindblad)
目次
  1. 2準位原子の Davies 生成子
  2. Davies の前提は満たされるか
  3. 量子回帰定理
  4. 木下の層

AI 生成・木下が未勉強

近似の数理 第14章 §14.5 ★★★ 2準位原子と真空の電磁場に Davies の定理を当てはめると A 係数と光 Bloch 方程式の緩和が出る。前提の可積分性は高い周波数の切断で満たされる。量子回帰定理は別の主張

2準位原子の Davies 生成子

2準位原子 H^S=ℏω0σ+σ−\hat H_S=\hbar\omega_0\sigma_+\sigma_- を、体積 VV の箱で量子化した電磁場に双極子相互作用 −d^⋅E^-\hat{\bm{d}}\cdot\hat{\bm{E}} で結合させる(光学と量子エレクトロニクス 第34章)。双極子の行列要素を実数にとって d^=deg(σ++σ−)\hat{\bm{d}}=\bm{d}_{eg}(\sigma_++\sigma_-) と書くと、(14.1) の形で A^=σx\hat A=\sigma_x、B^=−deg⋅E^\hat B=-\bm{d}_{eg}\cdot\hat{\bm{E}} で、Bohr 周波数は ±ω0\pm\omega_0、A^(ω0)=σ−\hat A(\omega_0)=\sigma_-、A^(−ω0)=σ+\hat A(-\omega_0)=\sigma_+ である。場の各モードとの結合を ℏgks=ℏωk/(2ε0V) (deg⋅eks)\hbar g_{\bm{k}s}=\sqrt{\hbar\omega_k/(2\varepsilon_0V)}\,(\bm{d}_{eg}\cdot\bm{e}_{\bm{k}s}) とおき、

J(ω)=∑k,s∣gks∣2δ(ω−ωk)=ω2ε0ℏV⋅Vω22π2c3⋅23deg2=ω3deg26π2ε0ℏc3\mathcal J(\omega)=\sum_{\bm{k},s}\lvert g_{\bm{k}s}\rvert ^2\delta(\omega-\omega_k)=\frac{\omega}{2\varepsilon_0\hbar V}\cdot\frac{V\omega^2}{2\pi^2c^3}\cdot\frac23d_{eg}^2=\frac{\omega^3d_{eg}^2}{6\pi^2\varepsilon_0\hbar c^3}

を得る(2つ目の因子は1つの偏光あたりのモードの密度、23deg2\frac23d_{eg}^2 は2つの偏光と向きについての ∣deg⋅e∣2\lvert \bm{d}_{eg}\cdot\bm{e}\rvert ^2 の和の平均)。温度 TT の場では (14.10) の速さが

γ(ω0)=Γ(nB(ω0)+1),γ(−ω0)=ΓnB(ω0),Γ=2πJ(ω0)=ω03deg23πε0ℏc3(14.19)\gamma(\omega_0)=\Gamma\bigl(n_{\mathrm B}(\omega_0)+1\bigr),\qquad \gamma(-\omega_0)=\Gamma n_{\mathrm B}(\omega_0),\qquad \Gamma=2\pi\mathcal J(\omega_0)=\frac{\omega_0^3d_{eg}^2}{3\pi\varepsilon_0\hbar c^3}\tag{14.19}

となり(λ\lambda は degd_{eg} に含めた)、Davies の生成子 (14.14) は例題 14.4 で γ↓=Γ(nB+1)\gamma_\downarrow=\Gamma(n_{\mathrm B}+1)、γ↑=ΓnB\gamma_\uparrow=\Gamma n_{\mathrm B}、γϕ=0\gamma_\phi=0 とおいたものになる。比 γ(−ω0)/γ(ω0)=nB/(nB+1)=e−βℏω0\gamma(-\omega_0)/\gamma(\omega_0)=n_{\mathrm B}/(n_{\mathrm B}+1)=\mathrm{e}^{-\beta\hbar\omega_0} は (14.16) そのものである。

ここで2つのことが同時に出ている。第一に、A^=σx\hat A=\sigma_x には σ+b^†\sigma_+\hat b^\dagger のような回転しない項が含まれるが、永年近似で A^(ω)ρ^A^(ω′)†\hat A(\omega)\hat\rho\hat A(\omega')^\dagger(ω≠ω′\omega\ne\omega')を落とすので、結合を最初から回転波近似で書いた場合と同じ生成子になる。原子と連続体の結合での回転波近似は、van Hove の極限の中に含まれている。第二に、温度 00 の速さ Γ\Gamma は自然放出の Einstein の A 係数 A21A_{21} で、温度 TT での吸収の速さ ΓnB\Gamma n_{\mathrm B} を、角周波数あたりのエネルギー密度u(ω)=ℏω3π2c3nB(ω)u(\omega)=\frac{\hbar\omega^3}{\pi^2c^3}n_{\mathrm B}(\omega) を使って B12u(ω0)B_{12}u(\omega_0) と書くと、

A21=ℏω03π2c3B21,B12=B21=πdeg23ε0ℏ2(14.20)A_{21}=\frac{\hbar\omega_0^3}{\pi^2c^3}B_{21},\qquad B_{12}=B_{21}=\frac{\pi d_{eg}^2}{3\varepsilon_0\hbar^2}\tag{14.20}

が出る(光学と量子エレクトロニクス 第20章の流儀。準位の縮退はない場合)。Einstein が熱平衡から要請した関係は、KMS の浴に対する Davies の生成子の性質として導かれる。

Davies の前提は満たされるか

温度 00 の相関関数は C(u)=ℏ2∫0∞J(ω)e−iωudωC(u)=\hbar^2\int_0^\infty\mathcal J(\omega)\mathrm{e}^{-\mathrm{i}\omega u}\mathrm{d}\omega である。J∝ω3\mathcal J\propto\omega^3 のままではこの積分は発散し、相関関数そのものが定義されない。双極子近似は波長が原子の大きさ(a0a_0 の程度)より長いことを仮定しており、ω≳c/a0\omega\gtrsim c/a_0 の高い周波数では結合が ω3\omega^3 より弱くなる。これを J(ω)∝ω3e−ω/ωc\mathcal J(\omega)\propto\omega^3\mathrm{e}^{-\omega/\omega_{\mathrm c}} の切断で表すと

∫0∞ω3e−ω/ωce−iωudω=6(ωc−1+iu)4\int_0^\infty\omega^3\mathrm{e}^{-\omega/\omega_{\mathrm c}}\mathrm{e}^{-\mathrm{i}\omega u}\mathrm{d}\omega=\frac{6}{(\omega_{\mathrm c}^{-1}+\mathrm{i} u)^4}

で、絶対値は u−4u^{-4} で減衰し、時間について可積分である。相関時間は τB∼1/ωc\tau_{\mathrm B}\sim1/\omega_{\mathrm c} で、Davies の定理の前提 (iii) はこの切断のもとで満たされる。定理の結論は切断の細部によらない Γ=2πJ(ω0)\Gamma=2\pi\mathcal J(\omega_0) を与えるが、Lamb のずれ S(ω0)S(\omega_0) は主値積分で ωc\omega_{\mathrm c} に依存し、非相対論的な扱いの範囲を超える(ここでは計算しない)。

例題 14.5 水素の 2p→1s の自然放出の速さ ★★

オリジナル水素原子の 2p→1s2\mathrm p\to1\mathrm s 遷移を2準位とみなす。行列要素は ∣⟨1s∣z∣2p0⟩∣=1282243a0\lvert \langle 1\mathrm s | z | 2\mathrm p_0 \rangle \rvert =\frac{128\sqrt2}{243}a_0、遷移のエネルギーは Rydberg エネルギーの 34\frac34 である。

  • (1) (14.19) の Γ\Gamma と寿命 1/Γ1/\Gamma を求めよ。
  • (2) Davies の極限の小さいパラメータ Γ/ω0\Gamma/\omega_0 を見積もれ。
  • (3) 室温(300 K300\,\mathrm K)の熱放射による励起の速さ ΓnB(ω0)\Gamma n_{\mathrm B}(\omega_0) は無視できるか。
解答

(1) ℏω0=34×13.606 eV=10.20 eV\hbar\omega_0=\frac34\times13.606\,\mathrm{eV}=10.20\,\mathrm{eV}、ω0=1.550×1016 s−1\omega_0=1.550\times10^{16}\,\mathrm{s^{-1}}。deg=e×0.7449 a0d_{eg}=e\times0.7449\,a_0 を入れるとΓ=ω03deg23πε0ℏc3=6.27×108 s−1\Gamma=\frac{\omega_0^3d_{eg}^2}{3\pi\varepsilon_0\hbar c^3}=6.27\times10^{8}\,\mathrm{s^{-1}}、寿命 1/Γ=1.60 ns1/\Gamma=1.60\,\mathrm{ns}(原子核を無限に重いとした値。2p±1\mathrm{2p}_{\pm1} からも同じ値になる)。

(2) Γ/ω0=4.0×10−8\Gamma/\omega_0=4.0\times10^{-8}。原子が1回振動する間に失う励起の割合がこの程度で、弱結合の極限の条件をよく満たす。切断の ωc∼c/a0≈6×1018 s−1\omega_{\mathrm c}\sim c/a_0\approx6\times10^{18}\,\mathrm{s^{-1}} を使えば、相関時間と寿命の比 Γ/ωc\Gamma/\omega_{\mathrm c} はさらに小さい。

(3) ℏω0/kBT=10.20 eV/0.0259 eV=395\hbar\omega_0/ k_{\mathrm B} T=10.20\,\mathrm{eV}/0.0259\,\mathrm{eV}=395 で、nB≈e−395n_{\mathrm B}\approx\mathrm{e}^{-395}。励起は完全に無視でき、緩和は温度 00 の自然放出だけで決まる。光の領域の遷移では室温でも浴はほぼ温度 00 であり、マイクロ波の遷移(次節の総合問題)とは事情が違う(検算スクリプト)。

近似の裏付け

型 B(極限) | 小さいもの Γ/ω0→0\Gamma/\omega_0\to0 と Γ/ωc→0\Gamma/\omega_{\mathrm c}\to0(原子と真空のモードの弱結合。時間を Γt\Gamma t 一定で伸ばす) | 証明 条件付き(Davies 1974 の前提のうち、浴の相関関数の可積分性は双極子近似の破れる周波数 ωc\omega_{\mathrm c} での切断で満たされる。回転波近似は極限に含まれる。Lamb のずれは切断に依存し定理の外) | 行き先 近似の数理 第7章

量子回帰定理

光の相関関数 g(1)g^{(1)}・g(2)g^{(2)}(光学と量子エレクトロニクス 第32章)には、2つの時刻の系の相関関数 ⟨X^(t+τ)Y^(t)⟩\langle \hat X(t+\tau)\hat Y(t) \rangle が要る。量子回帰定理は、これを1時刻の生成子で

⟨X^(t+τ)Y^(t)⟩=Tr⁡[X^ eτL(Y^ρ^S(t))](14.21)\langle \hat X(t+\tau)\hat Y(t) \rangle =\operatorname{Tr}\Bigl[\hat X\,\mathrm{e}^{\tau\mathcal L}\bigl(\hat Y\hat\rho_S(t)\bigr)\Bigr]\tag{14.21}

と計算してよい、と主張する。これは1時刻の Markov 性(ρ^S(t)=etLρ^S(0)\hat\rho_S(t)=\mathrm{e}^{t\mathcal L}\hat\rho_S(0))とは別の主張である。左辺は時刻 tt で系に Y^\hat Y を作用させた後の、系と浴を合わせた時間発展を含み、時刻 tt に系と浴の間にできた相関を引き継ぐ。右辺はその相関を捨て、時刻 tt で浴を ρ^B\hat\rho_B に戻したのと同じことをしている。G. Lindblad は、時間発展が半群であることだけからは多時刻の相関の Markov 性は出ないことを指摘した(Lindblad (1979), Commun. Math. Phys. 65, 281。二次文献で照合)。(14.21) が弱結合の極限で正しくなること(多時刻の相関関数の収束)は R. Dümcke が扱った(Dümcke (1983), J. Math. Phys. 24, 311)が、本書は原論文の前提を確かめていない。

例題 14.6 厳密に解ける模型で量子回帰定理を試す ★★

オリジナル例題 14.3 の模型(温度 00、∣e,0⟩| e,0 \rangle から出発、回転座標で c(t)c(t) は実数)で、⟨σ+(t+τ)σ−(t)⟩\langle \sigma_+(t+\tau)\sigma_-(t) \rangle を考える。

  • (1) 厳密な値が eiω0τc(t+τ)c(t)\mathrm{e}^{\mathrm{i}\omega_0\tau}c(t+\tau)c(t) であることを示せ。
  • (2) (14.21) の eτL\mathrm{e}^{\tau\mathcal L} の代わりに厳密な1時刻の写像 Φτ\Phi_\tau を使うと eiω0τ∣c(t)∣2c(τ)\mathrm{e}^{\mathrm{i}\omega_0\tau}\lvert c(t)\rvert ^2c(\tau) になることを示し、両者が一致する条件を述べよ。
  • (3) ΓτB=0.25\Gamma\tau_{\mathrm B}=0.25、t=τ=2τBt=\tau=2\tau_{\mathrm B} で2つの値を比べよ。
解答

(1) Heisenberg 描像の定義から ⟨σ+(t+τ)σ−(t)⟩=⟨ψ0∣U^(t+τ)†σ+U^(τ)σ−U^(t)∣ψ0⟩\langle \sigma_+(t+\tau)\sigma_-(t) \rangle =\langle \psi_0 | \hat U(t+\tau)^\dagger\sigma_+\hat U(\tau)\sigma_-\hat U(t) | \psi_0 \rangle。σ−U^(t)∣ψ0⟩\sigma_-\hat U(t)| \psi_0 \rangle は∣g,1k⟩| g,1_k \rangle の成分を消し、c(t)e−iω0t∣g,0⟩c(t)\mathrm{e}^{-\mathrm{i}\omega_0t}| g,0 \rangle だけを残す。∣g,0⟩| g,0 \rangle はエネルギー 00 の固有状態なので U^(τ)\hat U(\tau) で変わらない。左から (σ−U^(t+τ)∣ψ0⟩)†=c(t+τ)eiω0(t+τ)⟨g,0∣\bigl(\sigma_-\hat U(t+\tau)| \psi_0 \rangle \bigr)^\dagger=c(t+\tau)\mathrm{e}^{\mathrm{i}\omega_0(t+\tau)}\langle g,0 | をかけて eiω0τc(t+τ)c(t)\mathrm{e}^{\mathrm{i}\omega_0\tau}c(t+\tau)c(t)。

(2) σ−ρ^S(t)\sigma_-\hat\rho_S(t) のうち Tr⁡[σ+⋯ ]\operatorname{Tr}[\sigma_+\cdots] に残るのは ρee(t)∣g⟩⟨e∣\rho_{ee}(t)| g \rangle \langle e | の成分で、厳密な写像は振幅減衰の形 Φτ(∣g⟩⟨e∣)=c(τ)eiω0τ∣g⟩⟨e∣\Phi_\tau(| g \rangle \langle e |)=c(\tau)\mathrm{e}^{\mathrm{i}\omega_0\tau}| g \rangle \langle e | を持つので、eiω0τ∣c(t)∣2c(τ)\mathrm{e}^{\mathrm{i}\omega_0\tau}\lvert c(t)\rvert ^2c(\tau)。(1) と一致するのは c(t+τ)=c(t)c(τ)c(t+\tau)=c(t)c(\tau)、すなわち cc が指数関数のときに限る。指数関数は Markov の極限の答え e−Γt/2\mathrm{e}^{-\Gamma t/2} である。

(3) c(2τB)=0.863c(2\tau_{\mathrm B})=0.863、c(4τB)=0.665c(4\tau_{\mathrm B})=0.665 から、厳密な値(位相を除く)は 0.665×0.863=0.5740.665\times0.863=0.574、回帰させた値は 0.8632×0.863=0.6430.863^2\times0.863=0.643 で、12% 違う。1時刻の写像を厳密にしても、2時刻の相関は正しく出ない。時刻 tt に浴へ出た光子と系の相関(記憶)を、回帰の式は捨てているからである。ΓτB→0\Gamma\tau_{\mathrm B}\to0 では両者とも Markov の値に近づく(検算スクリプト)。

近似の裏付け

型 B(極限) | 小さいもの ΓτB→0\Gamma\tau_{\mathrm B}\to0(2つの時刻の両方を Γt\Gamma t 一定で伸ばす) | 証明 要確認(1時刻の主方程式は Davies 1974。多時刻の相関関数の弱結合極限での収束は Dümcke 1983 が扱ったとされるが原論文の前提を確かめていない。半群であることだけからは出ない〔Lindblad 1979〕) | 行き先 近似の数理 第7章

木下の勉強ステータス

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

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

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