解析力学から様々な理論分野を組み立てる 第5部 第12章 §5.12.3 ★★★ 加速度 a a a の観測者の正の振動数は Minkowski の正の振動数と指数関数でずれ、その Bogoliubov 係数が ∣ β / α ∣ = e − π c Ω / a \lvert \beta/\alpha\rvert =\mathrm{e}^{-\pi c\Omega/a} ∣ β / α ∣ = e − π c Ω/ a を満たす。Minkowski 真空はその観測者に温度 ℏ a / ( 2 π c k B ) \hbar a/(2\pi c k_{\mathrm B}) ℏ a / ( 2 π c k B ) の熱浴に見える
加速する観測者の正の振動数
右の楔に止まる観測者(ξ = \xi= ξ = 一定)の時計は t R t_{\mathrm R} t R で進む。この観測者が「正の振動数」と呼ぶのは e − i Ω t R \mathrm{e}^{-\mathrm{i}\Omega t_{\mathrm R}} e − i Ω t R で、右へ進むモードでは e − i Ω u \mathrm{e}^{-\mathrm{i}\Omega u} e − i Ω u である。これを Minkowski の光的座標で書くと、u = − c a ln ( − a U / c ) u=-\frac ca\ln(-aU/c) u = − a c ln ( − a U / c ) から
g Ω R = ℏ 4 π Ω e − i Ω u = ℏ 4 π Ω ( − a U c ) 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} g Ω R = 4 π Ω ℏ e − i Ω u = 4 π Ω ℏ ( − c a U ) i ν ( U < 0 ) , g Ω R = 0 ( U > 0 ) , ν = a c Ω ( 5.12.21 )
となる。ν \nu ν は無次元の数で、振動数を加速度の目盛り a / c a/c a / c で測ったものである。このモードは U → 0 − U\to0^- U → 0 − (楔の境界=観測者の地平面)で無限に速く振動し、U U U の関数として単一の振動数を持たない。観測者の時計で一定の振動数で振動する波は、Minkowski の時計では振動数が指数関数で変わる波に見える。
図 5.12.1 に楔の配置を示す。U = 0 U=0 U = 0 と V = 0 V=0 V = 0 の2本の光線が平面を4つに分け、右の楔 R と左の楔 L に加速する観測者の座標が張られる。右の楔の観測者には U > 0 U>0 U > 0 の領域から光が届かず、U = 0 U=0 U = 0 はこの観測者の地平面である。
図 5.12.1 Minkowski 時空の Rindler の楔。破線が光線 U = 0 U=0 U = 0 ・V = 0 V=0 V = 0 。太い双曲線は ξ = 0 \xi=0 ξ = 0 に止まる観測者(右の楔 R)と、その鏡像の観測者(左の楔 L)の世界線、細い双曲線は別の ξ \xi ξ 、原点から出る直線は t R = t_{\mathrm R}= t R = 一定の面。R の観測者には U = 0 U=0 U = 0 より上(領域 F と L)から光が届かない。
Bogoliubov 係数
式 (5.12.17) の係数を直接計算する。内積 (5.12.13) を部分積分して ( f ω , g ) = ( 2 ω / ℏ ) ∫ g f ω ∗ d U (f_\omega,g)=(2\omega/\hbar)\int g\,f_\omega^*\,\mathrm{d} U ( f ω , g ) = ( 2 ω /ℏ ) ∫ g f ω ∗ d U 、− ( f ω ∗ , g ) = ( 2 ω / ℏ ) ∫ g f ω d U -(f_\omega^*,g)=(2\omega/\hbar)\int g\,f_\omega\,\mathrm{d} U − ( f ω ∗ , g ) = ( 2 ω /ℏ ) ∫ g f ω d U を使い、x = − U x=-U x = − U とおくと
α Ω ω = 1 2 π ω Ω ∫ 0 ∞ ( a x c ) i ν e − i ω x d x , β Ω ω = 1 2 π ω Ω ∫ 0 ∞ ( a x c ) i ν e + i ω x d x . (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} α Ω ω = 2 π 1 Ω ω ∫ 0 ∞ ( c a x ) i ν e − i ω x d x , β Ω ω = 2 π 1 Ω ω ∫ 0 ∞ ( c a x ) i ν e + i ω x d x . ( 5.12.22 )
積分は Euler のガンマ関数の公式 ∫ 0 ∞ x s − 1 e − λ x d x = Γ ( s ) λ − s \int_0^\infty x^{s-1}\mathrm{e}^{-\lambda x}\mathrm{d} x=\Gamma(s)\lambda^{-s} ∫ 0 ∞ x s − 1 e − λ x d x = Γ ( s ) λ − s を λ = ± i ω \lambda=\pm\mathrm{i}\omega λ = ± i ω に解析接続したもので(e ∓ ϵ x \mathrm{e}^{\mp\epsilon x} e ∓ ϵ x を掛けて ϵ → 0 + \epsilon\to0^+ ϵ → 0 + とする)、( ± i ) − 1 − i ν = ∓ i e ± π ν / 2 (\pm\mathrm{i})^{-1-\mathrm{i}\nu}=\mp\mathrm{i}\,\mathrm{e}^{\pm\pi\nu/2} ( ± i ) − 1 − i ν = ∓ i e ± π ν /2 から
α Ω ω = − i 2 π ω Ω ( a c ) i ν Γ ( 1 + i ν ) ω 1 + i ν e π ν / 2 , β Ω ω = i 2 π ω Ω ( a c ) 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 π − i Ω ω ( c a ) i ν ω 1 + i ν Γ ( 1 + i ν ) e π ν /2 , β Ω ω = 2 π i Ω ω ( c a ) i ν ω 1 + i ν Γ ( 1 + i ν ) e − π ν /2 ( 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} β Ω ω = − e − π ν α Ω ω = − e − π c Ω/ a α Ω ω ( 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') ( 1 − e − 2 π ν ) ∫ d ω α Ω ω α Ω ′ ω ∗ = δ ( Ω − Ω ′ ) となり、式 (5.12.20) は
⟨ 0 M ∣ b ^ Ω R † b ^ Ω ′ R ∣ 0 M ⟩ = δ ( Ω − Ω ′ ) e 2 π 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 M ∣ b ^ Ω R † b ^ Ω ′ R ∣ 0 M ⟩ = e 2 π c Ω/ a − 1 δ ( Ω − Ω ′ ) ( 5.12.25 )
になる。δ ( 0 ) \delta(0) δ ( 0 ) は、連続の振動数で長い時間 t o b s t_{\mathrm{obs}} t obs 観測すると δ ( 0 ) → t o b s / ( 2 π ) \delta(0)\to t_{\mathrm{obs}}/(2\pi) δ ( 0 ) → t obs / ( 2 π ) になる量で、単位時間・単位角振動数あたり ( 2 π ) − 1 ( e 2 π c Ω / a − 1 ) − 1 (2\pi)^{-1}(\mathrm{e}^{2\pi c\Omega/a}-1)^{-1} ( 2 π ) − 1 ( e 2 π c Ω/ a − 1 ) − 1 個の粒子が観測者に届く。分母の形は Planck 分布 n B ( Ω ) = 1 / ( e ℏ Ω / ( k B T ) − 1 ) n_{\mathrm B}(\Omega)=1/(\mathrm{e}^{\hbar\Omega/( k_{\mathrm B} T)}-1) n B ( Ω ) = 1/ ( e ℏΩ/ ( k B T ) − 1 ) そのもので、温度は
ℏ Ω k B T U = 2 π c Ω a ⟹ T U = ℏ a 2 π c k B (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} k B T U ℏΩ = a 2 π c Ω ⟹ T U = 2 π c k B ℏ a ( 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) ∣ Γ ( 1 + i ν ) ∣ 2 = π ν / sinh ( π ν ) を使うと ∣ β Ω ω ∣ 2 = ν e − π ν / ( 4 π Ω ω sinh ( π ν ) ) \lvert \beta_{\Omega\omega}\rvert ^2=\nu\,\mathrm{e}^{-\pi\nu}/(4\pi\Omega\omega\sinh(\pi\nu)) ∣ β Ω ω ∣ 2 = ν e − π ν / ( 4 π Ω ω sinh ( π ν )) で、ω \omega ω について積分すると ∫ d ω / ω \int\mathrm{d}\omega/\omega ∫ d ω / ω が対数的に発散する。この発散が δ ( 0 ) \delta(0) δ ( 0 ) にあたり、時間が長いほど多く数えられることを表している。
解析性から見る理由
同じ結果は、積分をせずに解析性の議論で出すこともできる。Minkowski の正の振動数の関数 e − i ω U \mathrm{e}^{-\mathrm{i}\omega U} e − i ω U (ω > 0 \omega>0 ω > 0 )は、U U U を複素数に広げると下半面 Im U < 0 \operatorname{Im} U<0 Im U < 0 で有界で解析的である。逆に、下半面で有界で解析的な関数は正の振動数の成分だけを持つ。
右の楔のモード ( − a U / c ) i ν (-aU/c)^{\mathrm{i}\nu} ( − a U / c ) i ν を、U < 0 U<0 U < 0 から下半面を通って U > 0 U>0 U > 0 まで解析接続する。− U = ∣ U ∣ e i θ -U=\lvert U\rvert \mathrm{e}^{\mathrm{i}\theta} − U = ∣ U ∣ e i θ と書くと、下半面を回るとき偏角 θ \theta θ は 0 0 0 から π \pi π まで動き、U > 0 U>0 U > 0 に着いたとき ( − a U / c ) i ν = e − π ν ( a U / c ) i ν (-aU/c)^{\mathrm{i}\nu}=\mathrm{e}^{-\pi\nu}(aU/c)^{\mathrm{i}\nu} ( − a U / c ) i ν = e − π ν ( a U / c ) i ν となる。途中の絶対値は e − ν θ ≤ 1 \mathrm{e}^{-\nu\theta}\le1 e − ν θ ≤ 1 で有界である。したがって
h Ω ( U ) = ( − a U c ) i ν θ ( − U ) + e − π ν ( a U c ) 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} h Ω ( U ) = ( − c a U ) i ν θ ( − U ) + e − π ν ( c a U ) i ν θ ( U ) ( 5.12.27 )
は Minkowski の正の振動数だけでできている(θ ( ⋅ ) \theta(\cdot) θ ( ⋅ ) は階段関数)。U > 0 U>0 U > 0 の部分は左の楔のモードである。左の楔の観測者の未来向きの時間を t L t_{\mathrm L} t L とすると U = c a e a ( t L + ξ / c ) / c U=\frac ca\mathrm{e}^{a(t_{\mathrm L}+\xi/c)/c} U = a c e a ( t L + ξ / c ) / c で、左の観測者の正の振動数は ( a U / c ) − i ν (aU/c)^{-\mathrm{i}\nu} ( a U / c ) − i ν 、式 (5.12.27) の第2項はその複素共役(左の観測者の負の振動数)にあたる。左の楔のモード g Ω L g^{\mathrm L}_\Omega g Ω L を同じ規格化で定めると、h Ω ∝ g Ω R + e − π ν g Ω L ∗ h_\Omega\propto g^{\mathrm R}_\Omega+\mathrm{e}^{-\pi\nu}g^{\mathrm L*}_\Omega h Ω ∝ g Ω R + e − π ν g Ω L ∗ で、( h Ω , ϕ ^ ) (h_\Omega,\hat\phi) ( h Ω , ϕ ^ ) は Minkowski の消滅演算子だけの組み合わせだから ∣ 0 M ⟩ | 0_{\mathrm M} \rangle ∣ 0 M ⟩ を消す。( g ∗ , ϕ ^ ) = − b ^ † (g^*,\hat\phi)=-\hat b^\dagger ( g ∗ , ϕ ^ ) = − b ^ † を使うと、∣ 0 M ⟩ | 0_{\mathrm M} \rangle ∣ 0 M ⟩ は
( b ^ Ω R − e − π ν b ^ Ω L † ) ∣ 0 M ⟩ = 0 , ( b ^ Ω L − e − π ν b ^ Ω R † ) ∣ 0 M ⟩ = 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} ( b ^ Ω R − e − π ν b ^ Ω L † ) ∣ 0 M ⟩ = 0 , ( b ^ Ω L − e − π ν b ^ Ω R † ) ∣ 0 M ⟩ = 0 ( 5.12.28 )
を満たす(2つ目は U → − U U\to-U U → − U の鏡像で同じ議論をしたもの)。この2式から b ^ Ω R \hat b^{\mathrm R}_\Omega b ^ Ω R を解くと、⟨ 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}) ⟨ b ^ Ω R † b ^ Ω R ⟩ = e − 2 π ν / ( 1 − e − 2 π ν ) が出て、式 (5.12.25) と一致する(δ ( 0 ) \delta(0) δ ( 0 ) の因子を除く)。式 (5.12.28) は、Minkowski 真空が R と L のモードを対にした状態であることを示しており、その形は5.12.6 節で熱分布との対応に使う。
何が熱的なのか
Minkowski 真空は1つの純粋状態で、慣性系の観測者が見ても加速する観測者が見ても同じ状態である。変わるのは「粒子」の数え方で、加速する観測者は右の楔の中の観測量しか測れないので、左の楔の自由度について平均した状態を見ることになる。式 (5.12.28) の対の構造から、この平均が温度 T U T_{\mathrm U} T U のカノニカル分布になる(5.12.6 節)。熱的に見える理由は、地平面の向こう側の自由度を見られないことにある。
例題 5.12.3 地上の加速度と 1 K の Unruh 温度 ★★
オリジナル
加速度 a = 9.80665 m s − 2 a=9.80665\,\mathrm{m\,s^{-2}} a = 9.80665 m s − 2 の観測者の Unruh 温度を求めよ。
Unruh 温度が 1 K 1\,\mathrm K 1 K になる加速度を求めよ。
単位時間・単位角振動数あたりに届く粒子の数 ( 2 π ) − 1 ( e ℏ Ω / ( k B T U ) − 1 ) − 1 (2\pi)^{-1}(\mathrm{e}^{\hbar\Omega/( k_{\mathrm B} T_{\mathrm U})}-1)^{-1} ( 2 π ) − 1 ( e ℏΩ/ ( k B T U ) − 1 ) − 1 に ℏ Ω \hbar\Omega ℏΩ を掛けて積分し、右へ進む粒子が運ぶエネルギーの流れが π ( k B T U ) 2 / ( 12 ℏ ) \pi( k_{\mathrm B} T_{\mathrm U})^2/(12\hbar) π ( k B T U ) 2 / ( 12ℏ ) になることを示せ。
解答 (1) T U = ℏ a / ( 2 π c k B ) = 1.054571817 × 10 − 34 × 9.80665 / ( 2 π × 2.99792458 × 10 8 × 1.380649 × 10 − 23 ) = 3.98 × 10 − 20 K T_{\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 T U = ℏ a / ( 2 π c k B ) = 1.054571817 × 1 0 − 34 × 9.80665/ ( 2 π × 2.99792458 × 1 0 8 × 1.380649 × 1 0 − 23 ) = 3.98 × 1 0 − 20 K 。
(2) a = 2 π c k B × 1 K / ℏ = 2.47 × 10 20 m s − 2 a=2\pi c k_{\mathrm B}\times1\,\mathrm K/\hbar=2.47\times10^{20}\,\mathrm{m\,s^{-2}} a = 2 π c k B × 1 K /ℏ = 2.47 × 1 0 20 m s − 2 。地上の重力加速度の 2.5 × 10 19 2.5\times10^{19} 2.5 × 1 0 19 倍で、加速度と温度の比 ℏ / ( 2 π c k B ) = 4.06 × 10 − 21 K s 2 m − 1 \hbar/(2\pi c k_{\mathrm B})=4.06\times10^{-21}\,\mathrm{K\,s^2\,m^{-1}} ℏ/ ( 2 π c k B ) = 4.06 × 1 0 − 21 K s 2 m − 1 が小さいので、この効果を直接観測するのは難しい。
(3) x = ℏ Ω / ( k B T U ) x=\hbar\Omega/( k_{\mathrm B} T_{\mathrm U}) x = ℏΩ/ ( k B T U ) とおくと
∫ 0 ∞ d Ω 2 π ℏ Ω e ℏ Ω / ( k B T U ) − 1 = ( k B T U ) 2 2 π ℏ ∫ 0 ∞ x d x e x − 1 = ( k B T U ) 2 2 π ℏ ⋅ π 2 6 = π ( k B T U ) 2 12 ℏ . \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}. ∫ 0 ∞ 2 π d Ω e ℏΩ/ ( k B T U ) − 1 ℏΩ = 2 π ℏ ( k B T U ) 2 ∫ 0 ∞ e x − 1 x d x = 2 π ℏ ( k B T U ) 2 ⋅ 6 π 2 = 12ℏ π ( k B T U ) 2 .
これは1次元の黒体放射の片向きのエネルギーの流れで、3次元の Stefan–Boltzmann の法則の T 4 T^4 T 4 が、1次元では T 2 T^2 T 2 になる。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com