tecktom.com

出てくるもの1:Unruh 効果

2026年10月15日13分で読めますAI生成教科書『虎の巻』第5部 第12章 曲がった時空の場の量子論:Unruh 効果と Hawking 放射
目次
  1. 加速する観測者の正の振動数
  2. Bogoliubov 係数
  3. 解析性から見る理由
  4. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第5部 第12章 §5.12.3 ★★★ 加速度 aa の観測者の正の振動数は Minkowski の正の振動数と指数関数でずれ、その Bogoliubov 係数が ∣β/α∣=e−πcΩ/a\lvert \beta/\alpha\rvert =\mathrm{e}^{-\pi c\Omega/a} を満たす。Minkowski 真空はその観測者に温度 ℏa/(2πckB)\hbar a/(2\pi c k_{\mathrm B}) の熱浴に見える

加速する観測者の正の振動数

右の楔に止まる観測者(ξ=\xi= 一定)の時計は tRt_{\mathrm R} で進む。この観測者が「正の振動数」と呼ぶのは e−iΩtR\mathrm{e}^{-\mathrm{i}\Omega t_{\mathrm R}} で、右へ進むモードでは e−iΩu\mathrm{e}^{-\mathrm{i}\Omega u} である。これを Minkowski の光的座標で書くと、u=−caln⁡(−aU/c)u=-\frac ca\ln(-aU/c) から

gΩR=ℏ4πΩ  e−iΩu=ℏ4πΩ(−aUc)iν(U<0),gΩR=0(U>0),ν=cΩa(5.12.21)g^{\mathrm R}_\Omega=\sqrt{\frac{\hbar}{4\pi\Omega}}\;\mathrm{e}^{-\mathrm{i}\Omega u}=\sqrt{\frac{\hbar}{4\pi\Omega}}\Bigl(-\frac{aU}{c}\Bigr)^{\mathrm{i}\nu}\quad(U<0),\qquad g^{\mathrm R}_\Omega=0\quad(U>0), \qquad \nu=\frac{c\Omega}{a}\tag{5.12.21}

となる。ν\nu は無次元の数で、振動数を加速度の目盛り a/ca/c で測ったものである。このモードは U→0−U\to0^-(楔の境界=観測者の地平面)で無限に速く振動し、UU の関数として単一の振動数を持たない。観測者の時計で一定の振動数で振動する波は、Minkowski の時計では振動数が指数関数で変わる波に見える。

図 5.12.1 に楔の配置を示す。U=0U=0 と V=0V=0 の2本の光線が平面を4つに分け、右の楔 R と左の楔 L に加速する観測者の座標が張られる。右の楔の観測者には U>0U>0 の領域から光が届かず、U=0U=0 はこの観測者の地平面である。

図

図 5.12.1 Minkowski 時空の Rindler の楔。破線が光線 U=0U=0・V=0V=0。太い双曲線は ξ=0\xi=0 に止まる観測者(右の楔 R)と、その鏡像の観測者(左の楔 L)の世界線、細い双曲線は別の ξ\xi、原点から出る直線は tR=t_{\mathrm R}= 一定の面。R の観測者には U=0U=0 より上(領域 F と L)から光が届かない。

Bogoliubov 係数

式 (5.12.17) の係数を直接計算する。内積 (5.12.13) を部分積分して (fω,g)=(2ω/ℏ)∫g fω∗ dU(f_\omega,g)=(2\omega/\hbar)\int g\,f_\omega^*\,\mathrm{d} U、−(fω∗,g)=(2ω/ℏ)∫g fω dU-(f_\omega^*,g)=(2\omega/\hbar)\int g\,f_\omega\,\mathrm{d} U を使い、x=−Ux=-U とおくと

αΩω=12πωΩ∫0∞(axc)iνe−iωx dx,βΩω=12πωΩ∫0∞(axc)iνe+iωx dx.(5.12.22)\alpha_{\Omega\omega}=\frac{1}{2\pi}\sqrt{\frac{\omega}{\Omega}}\int_0^\infty\Bigl(\frac{ax}{c}\Bigr)^{\mathrm{i}\nu}\mathrm{e}^{-\mathrm{i}\omega x}\,\mathrm{d} x,\qquad \beta_{\Omega\omega}=\frac{1}{2\pi}\sqrt{\frac{\omega}{\Omega}}\int_0^\infty\Bigl(\frac{ax}{c}\Bigr)^{\mathrm{i}\nu}\mathrm{e}^{+\mathrm{i}\omega x}\,\mathrm{d} x .\tag{5.12.22}

積分は Euler のガンマ関数の公式 ∫0∞xs−1e−λxdx=Γ(s)λ−s\int_0^\infty x^{s-1}\mathrm{e}^{-\lambda x}\mathrm{d} x=\Gamma(s)\lambda^{-s} を λ=±iω\lambda=\pm\mathrm{i}\omega に解析接続したもので(e∓ϵx\mathrm{e}^{\mp\epsilon x} を掛けて ϵ→0+\epsilon\to0^+ とする)、(±i)−1−iν=∓i e±πν/2(\pm\mathrm{i})^{-1-\mathrm{i}\nu}=\mp\mathrm{i}\,\mathrm{e}^{\pm\pi\nu/2} から

αΩω=−i2πωΩ(ac)iνΓ(1+iν)ω1+iνeπν/2,βΩω=i2πωΩ(ac)iνΓ(1+iν)ω1+iνe−πν/2(5.12.23)\alpha_{\Omega\omega}=\frac{-\mathrm{i}}{2\pi}\sqrt{\frac{\omega}{\Omega}}\Bigl(\frac{a}{c}\Bigr)^{\mathrm{i}\nu}\frac{\Gamma(1+\mathrm{i}\nu)}{\omega^{1+\mathrm{i}\nu}}\mathrm{e}^{\pi\nu/2},\qquad \beta_{\Omega\omega}=\frac{\mathrm{i}}{2\pi}\sqrt{\frac{\omega}{\Omega}}\Bigl(\frac{a}{c}\Bigr)^{\mathrm{i}\nu}\frac{\Gamma(1+\mathrm{i}\nu)}{\omega^{1+\mathrm{i}\nu}}\mathrm{e}^{-\pi\nu/2}\tag{5.12.23}

を得る。2つは前の因子が同じで、

βΩω=−e−πναΩω=−e−πcΩ/aαΩω(5.12.24)\beta_{\Omega\omega}=-\mathrm{e}^{-\pi\nu}\alpha_{\Omega\omega}=-\mathrm{e}^{-\pi c\Omega/a}\alpha_{\Omega\omega}\tag{5.12.24}

が、すべての ω\omega について成り立つ。β\beta は0でない。加速する観測者のモードには、Minkowski の負の振動数が混ざっている。

規格化の条件 (5.12.19) に式 (5.12.24) を入れると (1−e−2πν)∫dω αΩωαΩ′ω∗=δ(Ω−Ω′)(1-\mathrm{e}^{-2\pi\nu})\int\mathrm{d}\omega\,\alpha_{\Omega\omega}\alpha_{\Omega'\omega}^*=\delta(\Omega-\Omega') となり、式 (5.12.20) は

⟨0M∣b^ΩR†b^Ω′R∣0M⟩=δ(Ω−Ω′)e2πcΩ/a−1(5.12.25)\langle 0_{\mathrm M} | \hat b_\Omega^{\mathrm R\dagger}\hat b^{\mathrm R}_{\Omega'} | 0_{\mathrm M} \rangle =\frac{\delta(\Omega-\Omega')}{\mathrm{e}^{2\pi c\Omega/a}-1}\tag{5.12.25}

になる。δ(0)\delta(0) は、連続の振動数で長い時間 tobst_{\mathrm{obs}} 観測すると δ(0)→tobs/(2π)\delta(0)\to t_{\mathrm{obs}}/(2\pi) になる量で、単位時間・単位角振動数あたり (2π)−1(e2πcΩ/a−1)−1(2\pi)^{-1}(\mathrm{e}^{2\pi c\Omega/a}-1)^{-1} 個の粒子が観測者に届く。分母の形は Planck 分布 nB(Ω)=1/(eℏΩ/(kBT)−1)n_{\mathrm B}(\Omega)=1/(\mathrm{e}^{\hbar\Omega/( k_{\mathrm B} T)}-1) そのもので、温度は

ℏΩkBTU=2πcΩa⟹TU=ℏa2πckB(5.12.26)\frac{\hbar\Omega}{ k_{\mathrm B} T_{\mathrm U}}=\frac{2\pi c\Omega}{a}\quad\Longrightarrow\quad T_{\mathrm U}=\frac{\hbar a}{2\pi c k_{\mathrm B}}\tag{5.12.26}

である。これが Unruh 温度である。この結論は Fulling(1973、Physical Review D 7、2850)が量子化の一意性の問題として、Davies(1975、Journal of Physics A 8、609)が Rindler 計量の粒子の生成として、Unruh(1976、Physical Review D 14、870)が加速する検出器の応答として、それぞれ示した。

∣Γ(1+iν)∣2=πν/sinh⁡(πν)\lvert \Gamma(1+\mathrm{i}\nu)\rvert ^2=\pi\nu/\sinh(\pi\nu) を使うと ∣βΩω∣2=ν e−πν/(4πΩωsinh⁡(πν))\lvert \beta_{\Omega\omega}\rvert ^2=\nu\,\mathrm{e}^{-\pi\nu}/(4\pi\Omega\omega\sinh(\pi\nu)) で、ω\omega について積分すると ∫dω/ω\int\mathrm{d}\omega/\omega が対数的に発散する。この発散が δ(0)\delta(0) にあたり、時間が長いほど多く数えられることを表している。

解析性から見る理由

同じ結果は、積分をせずに解析性の議論で出すこともできる。Minkowski の正の振動数の関数 e−iωU\mathrm{e}^{-\mathrm{i}\omega U}(ω>0\omega>0)は、UU を複素数に広げると下半面 Im⁡U<0\operatorname{Im} U<0 で有界で解析的である。逆に、下半面で有界で解析的な関数は正の振動数の成分だけを持つ。

右の楔のモード (−aU/c)iν(-aU/c)^{\mathrm{i}\nu} を、U<0U<0 から下半面を通って U>0U>0 まで解析接続する。−U=∣U∣eiθ-U=\lvert U\rvert \mathrm{e}^{\mathrm{i}\theta} と書くと、下半面を回るとき偏角 θ\theta は 00 から π\pi まで動き、U>0U>0 に着いたとき (−aU/c)iν=e−πν(aU/c)iν(-aU/c)^{\mathrm{i}\nu}=\mathrm{e}^{-\pi\nu}(aU/c)^{\mathrm{i}\nu} となる。途中の絶対値は e−νθ≤1\mathrm{e}^{-\nu\theta}\le1 で有界である。したがって

hΩ(U)=(−aUc)iνθ(−U)+e−πν(aUc)iνθ(U)(5.12.27)h_\Omega(U)=\Bigl(-\frac{aU}{c}\Bigr)^{\mathrm{i}\nu}\theta(-U)+\mathrm{e}^{-\pi\nu}\Bigl(\frac{aU}{c}\Bigr)^{\mathrm{i}\nu}\theta(U)\tag{5.12.27}

は Minkowski の正の振動数だけでできている(θ(⋅)\theta(\cdot) は階段関数)。U>0U>0 の部分は左の楔のモードである。左の楔の観測者の未来向きの時間を tLt_{\mathrm L} とすると U=caea(tL+ξ/c)/cU=\frac ca\mathrm{e}^{a(t_{\mathrm L}+\xi/c)/c} で、左の観測者の正の振動数は (aU/c)−iν(aU/c)^{-\mathrm{i}\nu}、式 (5.12.27) の第2項はその複素共役(左の観測者の負の振動数)にあたる。左の楔のモード gΩLg^{\mathrm L}_\Omega を同じ規格化で定めると、hΩ∝gΩR+e−πνgΩL∗h_\Omega\propto g^{\mathrm R}_\Omega+\mathrm{e}^{-\pi\nu}g^{\mathrm L*}_\Omega で、(hΩ,ϕ^)(h_\Omega,\hat\phi) は Minkowski の消滅演算子だけの組み合わせだから ∣0M⟩| 0_{\mathrm M} \rangle を消す。(g∗,ϕ^)=−b^†(g^*,\hat\phi)=-\hat b^\dagger を使うと、∣0M⟩| 0_{\mathrm M} \rangle は

(b^ΩR−e−πνb^ΩL†)∣0M⟩=0,(b^ΩL−e−πνb^ΩR†)∣0M⟩=0(5.12.28)\bigl(\hat b^{\mathrm R}_\Omega-\mathrm{e}^{-\pi\nu}\hat b^{\mathrm L\dagger}_\Omega\bigr)| 0_{\mathrm M} \rangle =0,\qquad \bigl(\hat b^{\mathrm L}_\Omega-\mathrm{e}^{-\pi\nu}\hat b^{\mathrm R\dagger}_\Omega\bigr)| 0_{\mathrm M} \rangle =0\tag{5.12.28}

を満たす(2つ目は U→−UU\to-U の鏡像で同じ議論をしたもの)。この2式から b^ΩR\hat b^{\mathrm R}_\Omega を解くと、⟨b^ΩR†b^ΩR⟩=e−2πν/(1−e−2πν)\langle \hat b^{\mathrm R\dagger}_\Omega\hat b^{\mathrm R}_\Omega \rangle =\mathrm{e}^{-2\pi\nu}/(1-\mathrm{e}^{-2\pi\nu}) が出て、式 (5.12.25) と一致する(δ(0)\delta(0) の因子を除く)。式 (5.12.28) は、Minkowski 真空が R と L のモードを対にした状態であることを示しており、その形は5.12.6節で熱分布との対応に使う。

何が熱的なのか

Minkowski 真空は1つの純粋状態で、慣性系の観測者が見ても加速する観測者が見ても同じ状態である。変わるのは「粒子」の数え方で、加速する観測者は右の楔の中の観測量しか測れないので、左の楔の自由度について平均した状態を見ることになる。式 (5.12.28) の対の構造から、この平均が温度 TUT_{\mathrm U} のカノニカル分布になる(5.12.6節)。熱的に見える理由は、地平面の向こう側の自由度を見られないことにある。

例題 5.12.3 地上の加速度と 1 K の Unruh 温度 ★★

オリジナル

  1. 加速度 a=9.80665 m s−2a=9.80665\,\mathrm{m\,s^{-2}} の観測者の Unruh 温度を求めよ。
  2. Unruh 温度が 1 K1\,\mathrm K になる加速度を求めよ。
  3. 単位時間・単位角振動数あたりに届く粒子の数 (2π)−1(eℏΩ/(kBTU)−1)−1(2\pi)^{-1}(\mathrm{e}^{\hbar\Omega/( k_{\mathrm B} T_{\mathrm U})}-1)^{-1} に ℏΩ\hbar\Omega を掛けて積分し、右へ進む粒子が運ぶエネルギーの流れが π(kBTU)2/(12ℏ)\pi( k_{\mathrm B} T_{\mathrm U})^2/(12\hbar) になることを示せ。
解答

(1) TU=ℏa/(2πckB)=1.054571817×10−34×9.80665/(2π×2.99792458×108×1.380649×10−23)=3.98×10−20 KT_{\mathrm U}=\hbar a/(2\pi c k_{\mathrm B})=1.054571817\times10^{-34}\times9.80665/(2\pi\times2.99792458\times10^8\times1.380649\times10^{-23})=3.98\times10^{-20}\,\mathrm K。

(2) a=2πckB×1 K/ℏ=2.47×1020 m s−2a=2\pi c k_{\mathrm B}\times1\,\mathrm K/\hbar=2.47\times10^{20}\,\mathrm{m\,s^{-2}}。地上の重力加速度の 2.5×10192.5\times10^{19} 倍で、加速度と温度の比 ℏ/(2πckB)=4.06×10−21 K s2 m−1\hbar/(2\pi c k_{\mathrm B})=4.06\times10^{-21}\,\mathrm{K\,s^2\,m^{-1}} が小さいので、この効果を直接観測するのは難しい。

(3) x=ℏΩ/(kBTU)x=\hbar\Omega/( k_{\mathrm B} T_{\mathrm U}) とおくと

∫0∞dΩ2π ℏΩeℏΩ/(kBTU)−1=(kBTU)22πℏ∫0∞x dxex−1=(kBTU)22πℏ⋅π26=π(kBTU)212ℏ.\int_0^\infty\frac{\mathrm{d}\Omega}{2\pi}\,\frac{\hbar\Omega}{\mathrm{e}^{\hbar\Omega/( k_{\mathrm B} T_{\mathrm U})}-1}=\frac{( k_{\mathrm B} T_{\mathrm U})^2}{2\pi\hbar}\int_0^\infty\frac{x\,\mathrm{d} x}{\mathrm{e}^x-1} =\frac{( k_{\mathrm B} T_{\mathrm U})^2}{2\pi\hbar}\cdot\frac{\pi^2}{6}=\frac{\pi( k_{\mathrm B} T_{\mathrm U})^2}{12\hbar}.

これは1次元の黒体放射の片向きのエネルギーの流れで、3次元の Stefan–Boltzmann の法則の T4T^4 が、1次元では T2T^2 になる。

ほかの巻では

木下の勉強ステータス

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

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

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