解析力学から様々な理論分野を組み立てる 第2部 第16章 §2.16.5 ★★ 同じ屈折率 n = 1 + 2 G M / ( c 2 r ) n=1+2GM/(c^2r) n = 1 + 2 GM / ( c 2 r ) から、光の到着の遅れ(Shapiro 遅延)と、像の位置と増光(重力レンズ)を出す
光の到着の遅れ
(2.16.24) の屈折率は、光の進路を曲げるだけでなく、光が届くまでの時間も延ばす。直線の経路(曲がりの効果は時間に対して2次)に沿ってd t = n d z / c \mathrm{d} t=n\,\mathrm{d} z/c d t = n d z / c を積分すると、平坦な時空より余計にかかる座標時は
Δ t = ∫ − z 1 z 2 n − 1 c d z = 2 G M c 3 ∫ − z 1 z 2 d z b 2 + z 2 = 2 G M c 3 ln ( z 1 + z 1 2 + b 2 ) ( z 2 + z 2 2 + b 2 ) b 2 (2.16.26) \Delta t=\int_{-z_1}^{z_2}\frac{n-1}{c}\,\mathrm{d} z=\frac{2GM}{c^3}\int_{-z_1}^{z_2}\frac{\mathrm{d} z}{\sqrt{b^2+z^2}}
=\frac{2GM}{c^3}\ln\frac{\bigl(z_1+\sqrt{z_1^2+b^2}\bigr)\bigl(z_2+\sqrt{z_2^2+b^2}\bigr)}{b^2}\tag{2.16.26} Δ t = ∫ − z 1 z 2 c n − 1 d z = c 3 2 GM ∫ − z 1 z 2 b 2 + z 2 d z = c 3 2 GM ln b 2 ( z 1 + z 1 2 + b 2 ) ( z 2 + z 2 2 + b 2 ) ( 2.16.26 )
である(z 1 z_1 z 1 , z 2 z_2 z 2 は最接近点から発信点・受信点までの距離)。両端が太陽から十分遠く z 1 , 2 ≈ r 1 , 2 ≫ b z_{1,2}\approx r_{1,2}\gg b z 1 , 2 ≈ r 1 , 2 ≫ b なら
Δ t ≈ 2 G M c 3 ln 4 r 1 r 2 b 2 = r s c ln 4 r 1 r 2 b 2 (2.16.27) \Delta t\approx\frac{2GM}{c^3}\ln\frac{4r_1r_2}{b^2}=\frac{r_{\mathrm s}}{c}\ln\frac{4r_1r_2}{b^2}\tag{2.16.27} Δ t ≈ c 3 2 GM ln b 2 4 r 1 r 2 = c r s ln b 2 4 r 1 r 2 ( 2.16.27 )
となる。対数の中の量は座標の取り方で O ( 1 ) \mathcal{O}(1) O ( 1 ) だけ変わるが、b b b を変えたときの Δ t \Delta t Δ t の変化、つまり惑星が太陽の向こう側を通る前後の時間変化は座標によらない。I. I. Shapiro はこれを「一般相対論の第4の検証」として、金星や水星にレーダー波を当てて往復時間を測ることを提案した1 1 I. I. Shapiro (1964), Physical Review Letters 13, 789–791. 。最も精密な測定は、土星へ向かう探査機 Cassini との電波の往復で、γ − 1 = ( 2.1 ± 2.3 ) × 10 − 5 \gamma-1=(2.1\pm2.3)\times10^{-5} γ − 1 = ( 2.1 ± 2.3 ) × 1 0 − 5 を得ている2 2 B. Bertotti, L. Iess, P. Tortora (2003), Nature 425, 374–376. 。
例題 2.16.5 金星のレーダーの往復の遅れ ★★
オリジナル 地球(太陽から 1 a u 1\,\mathrm{au} 1 au )から金星(0.723332 a u 0.723332\,\mathrm{au} 0.723332 au )へレーダー波を送り、反射波を受ける。金星は太陽のほぼ真裏にあり(外合)、電波は太陽の縁 b = R ⊙ = 6.957 × 10 8 m b=R_\odot=6.957\times10^8\,\mathrm m b = R ⊙ = 6.957 × 1 0 8 m をかすめる。往復で余計にかかる時間を (2.16.27) で求めよ。
解答 片道の遅れは r s c ln 4 r ⊕ r V b 2 \frac{r_{\mathrm s}}{c}\ln\frac{4r_\oplus r_{\mathrm V}}{b^2} c r s ln b 2 4 r ⊕ r V 。r s / c = 2953 / 2.998 × 10 8 = 9.85 × 10 − 6 s r_{\mathrm s}/c=2953/2.998\times10^8=9.85\times10^{-6}\,\mathrm s r s / c = 2953/2.998 × 1 0 8 = 9.85 × 1 0 − 6 s 、4 r ⊕ r V / b 2 = 4 × 1.496 × 10 11 × 1.0821 × 10 11 / ( 6.957 × 10 8 ) 2 = 1.338 × 10 5 4r_\oplus r_{\mathrm V}/b^2=4\times1.496\times10^{11}\times1.0821\times10^{11}/(6.957\times10^8)^2=1.338\times10^5 4 r ⊕ r V / b 2 = 4 × 1.496 × 1 0 11 × 1.0821 × 1 0 11 / ( 6.957 × 1 0 8 ) 2 = 1.338 × 1 0 5 で、その対数は 11.80 11.80 11.80 。片道 1.16 × 10 − 4 s 1.16\times10^{-4}\,\mathrm s 1.16 × 1 0 − 4 s 、往復でその2倍の 2.33 × 10 − 4 s 2.33\times10^{-4}\,\mathrm s 2.33 × 1 0 − 4 s 、つまり約 233 μ s 233\,\mu\mathrm s 233 μ s である。
往復の時間そのもの(約 1.7 × 10 3 s 1.7\times10^3\,\mathrm s 1.7 × 1 0 3 s )に比べると 10 − 7 10^{-7} 1 0 − 7 の大きさだが、レーダーの時間の分解能(μ s \mu\mathrm{s} μ s 以下)で十分に測ることができる。外合をはさんで b b b が変わると遅れが対数的に増えて減るので、往復時間の変化の形として取り出すことができる。
重力レンズ
光の曲がり (2.16.23) を、遠い光源と観測者の間に置かれた質量(レンズ)に当てる。観測者からレンズまでの距離を D l D_{\mathrm l} D l 、光源までを D s D_{\mathrm s} D s 、レンズから光源までを D l s D_{\mathrm{ls}} D ls とする。角度はすべて小さいとし、曲がりはレンズの位置で一度に起きるとみなす(薄いレンズの近似)。光源の本当の方向を β \beta β 、像の見える方向を θ \theta θ (どちらもレンズの方向から測る)とすると、光はレンズの位置で衝突径数 b = D l θ b=D_{\mathrm l}\theta b = D l θ を通って ϑ = 4 G M / ( c 2 D l θ ) \vartheta=4GM/(c^2D_{\mathrm l}\theta) ϑ = 4 GM / ( c 2 D l θ ) だけ曲がる。光源の面での位置のずれを比べて
β = θ − θ E 2 θ , θ E 2 ≔ 4 G M c 2 D l s D l D s (2.16.28) \beta=\theta-\frac{\theta_{\mathrm E}^2}{\theta},\qquad \theta_{\mathrm E}^2\coloneqq\frac{4GM}{c^2}\frac{D_{\mathrm{ls}}}{D_{\mathrm l}D_{\mathrm s}}\tag{2.16.28} β = θ − θ θ E 2 , θ E 2 : = c 2 4 GM D l D s D ls ( 2.16.28 )
を得る。これをレンズ方程式と呼ぶ。θ \theta θ の2次方程式 θ 2 − β θ − θ E 2 = 0 \theta^2-\beta\theta-\theta_{\mathrm E}^2=0 θ 2 − β θ − θ E 2 = 0 なので、像は
θ ± = β ± β 2 + 4 θ E 2 2 (2.16.29) \theta_\pm=\frac{\beta\pm\sqrt{\beta^2+4\theta_{\mathrm E}^2}}{2}\tag{2.16.29} θ ± = 2 β ± β 2 + 4 θ E 2 ( 2.16.29 )
の2つできる。θ + \theta_+ θ + はレンズから見て光源と同じ側、θ − \theta_- θ − (負)は反対側に見える。光源がレンズの真後ろ(β = 0 \beta=0 β = 0 )にあると、軸のまわりの対称性から像は半径 θ E \theta_{\mathrm E} θ E の輪になる(Einstein リング)。1979 年には、約 6 ′ ′ 6'' 6 ′′ 離れて並ぶ2つのクエーサーのスペクトルが同じであることから、1つのクエーサーの2つの像であることが見つかった3 3 D. Walsh, R. F. Carswell, R. J. Weymann (1979), Nature 279, 381–384. 。
光の曲がりは光の面積あたりの明るさ(表面輝度)を変えない。これは Liouville の定理から出るエタンデュの保存(第2部第10章 )である。したがって像の明るさの比は像の立体角の比に等しく、増光率は μ ± = ∣ θ ± β d θ ± d β ∣ \mu_\pm=\lvert \frac{\theta_\pm}{\beta}\frac{\mathrm{d} \theta_\pm}{\mathrm{d} \beta}\rvert μ ± = ∣ β θ ± d β d θ ± ∣ で与えられる。x ≔ β / θ E x\coloneqq\beta/\theta_{\mathrm E} x : = β / θ E とおいて2つの像の和をとると
μ t o t = μ + + μ − = x 2 + 2 x x 2 + 4 (2.16.30) \mu_{\mathrm{tot}}=\mu_++\mu_-=\frac{x^2+2}{x\sqrt{x^2+4}}\tag{2.16.30} μ tot = μ + + μ − = x x 2 + 4 x 2 + 2 ( 2.16.30 )
となる。銀河系の中の星がレンズになるとき、θ E \theta_{\mathrm E} θ E はミリ秒角ほどで2つの像は分けて見えないが、レンズの星が光源の前を横切ると x x x が時間とともに変わり、光源が一時的に明るくなる。これが重力マイクロレンズで、B. Paczyński が銀河系の暗い天体を探す方法として提案した4 4 B. Paczyński (1986), The Astrophysical Journal 304, 1–5. 。
図 2.16.3 薄いレンズの配置(角度は誇張した)。光源から出た光は衝突径数 b = D l θ b=D_{\mathrm l}\theta b = D l θ でレンズの横を通り、角度 ϑ = 4 G M / ( c 2 b ) \vartheta=4GM/(c^2b) ϑ = 4 GM / ( c 2 b ) 曲がって観測者に届く。観測者は像を、光源の本当の方向 β \beta β でなく θ \theta θ の方向に見る。点線は、曲がりが無いときの光の延長と、光源の本当の方向。
例題 2.16.6 銀河系の中心方向の星によるマイクロレンズ ★★
オリジナル 太陽質量の星(G M ⊙ = 1.32712440018 × 10 20 m 3 s − 2 GM_\odot=1.32712440018\times10^{20}\,\mathrm{m^3\,s^{-2}} G M ⊙ = 1.32712440018 × 1 0 20 m 3 s − 2 )が、地球から D l = 4 k p c D_{\mathrm l}=4\,\mathrm{kpc} D l = 4 kpc の所にあり、その真後ろ付近の D s = 8 k p c D_{\mathrm s}=8\,\mathrm{kpc} D s = 8 kpc に光源の星がある(1 k p c = 3.0857 × 10 19 m 1\,\mathrm{kpc}=3.0857\times10^{19}\,\mathrm m 1 kpc = 3.0857 × 1 0 19 m )。
(1) Einstein 半径 θ E \theta_{\mathrm E} θ E をミリ秒角で求めよ。
(2) 光源がレンズから β = θ E / 2 \beta=\theta_{\mathrm E}/2 β = θ E /2 だけずれているとき、2つの像の位置を θ E \theta_{\mathrm E} θ E を単位に求め、全体の増光率 μ t o t \mu_{\mathrm{tot}} μ tot を求めよ。
(3) (2.16.30) を μ ± \mu_\pm μ ± の和として導け。
解答 (1) D l s = 4 k p c D_{\mathrm{ls}}=4\,\mathrm{kpc} D ls = 4 kpc なので θ E 2 = 4 G M c 2 1 D s = 2 r s 8 k p c = 5906 m 2.469 × 10 20 m = 2.39 × 10 − 17 \theta_{\mathrm E}^2=\frac{4GM}{c^2}\frac{1}{D_{\mathrm s}}=\frac{2r_{\mathrm s}}{8\,\mathrm{kpc}}=\frac{5906\,\mathrm m}{2.469\times10^{20}\,\mathrm m}=2.39\times10^{-17} θ E 2 = c 2 4 GM D s 1 = 8 kpc 2 r s = 2.469 × 1 0 20 m 5906 m = 2.39 × 1 0 − 17 。θ E = 4.89 × 10 − 9 r a d = 1.01 \theta_{\mathrm E}=4.89\times10^{-9}\,\mathrm{rad}=1.01 θ E = 4.89 × 1 0 − 9 rad = 1.01 ミリ秒角である。
(2) x = 1 / 2 x=1/2 x = 1/2 を (2.16.29) に入れると θ ± / θ E = ( 0.5 ± 4.25 ) / 2 \theta_\pm/\theta_{\mathrm E}=(0.5\pm\sqrt{4.25})/2 θ ± / θ E = ( 0.5 ± 4.25 ) /2 、つまり θ + = 1.28 θ E \theta_+=1.28\,\theta_{\mathrm E} θ + = 1.28 θ E 、θ − = − 0.78 θ E \theta_-=-0.78\,\theta_{\mathrm E} θ − = − 0.78 θ E 。2つの像は 2.06 θ E ≈ 2 2.06\,\theta_{\mathrm E}\approx2 2.06 θ E ≈ 2 ミリ秒角しか離れていないので、望遠鏡では1つに見える。増光率は μ t o t = ( 0.25 + 2 ) / ( 0.5 4.25 ) = 2.18 \mu_{\mathrm{tot}}=(0.25+2)/(0.5\sqrt{4.25})=2.18 μ tot = ( 0.25 + 2 ) / ( 0.5 4.25 ) = 2.18 で、光源は約2倍明るく見える。
(3) θ ± \theta_\pm θ ± を θ E \theta_{\mathrm E} θ E を単位に y ± = ( x ± x 2 + 4 ) / 2 y_\pm=(x\pm\sqrt{x^2+4})/2 y ± = ( x ± x 2 + 4 ) /2 と書く。d y ± d x = 1 2 ( 1 ± x x 2 + 4 ) \frac{\mathrm{d} y_\pm}{\mathrm{d} x}=\frac12\bigl(1\pm\frac{x}{\sqrt{x^2+4}}\bigr) d x d y ± = 2 1 ( 1 ± x 2 + 4 x ) なので
μ ± = ∣ y ± x d y ± d x ∣ = ( x 2 + 4 ± x ) 2 4 x x 2 + 4 \mu_\pm=\lvert \frac{y_\pm}{x}\frac{\mathrm{d} y_\pm}{\mathrm{d} x}\rvert =\frac{(\sqrt{x^2+4}\pm x)^2}{4x\sqrt{x^2+4}} μ ± = ∣ x y ± d x d y ± ∣ = 4 x x 2 + 4 ( x 2 + 4 ± x ) 2
(μ − \mu_- μ − では y − < 0 y_-<0 y − < 0 なので絶対値をとった)。和をとると分子は ( x 2 + 4 + x ) 2 + ( x 2 + 4 − x ) 2 = 2 ( x 2 + 4 ) + 2 x 2 = 4 ( x 2 + 2 ) (\sqrt{x^2+4}+x)^2+(\sqrt{x^2+4}-x)^2=2(x^2+4)+2x^2=4(x^2+2) ( x 2 + 4 + x ) 2 + ( x 2 + 4 − x ) 2 = 2 ( x 2 + 4 ) + 2 x 2 = 4 ( x 2 + 2 ) で、μ t o t = x 2 + 2 x x 2 + 4 \mu_{\mathrm{tot}}=\frac{x^2+2}{x\sqrt{x^2+4}} μ tot = x x 2 + 4 x 2 + 2 を得る。x = 1 x=1 x = 1 (光源が Einstein 半径の所)では μ t o t = 3 / 5 = 1.34 \mu_{\mathrm{tot}}=3/\sqrt5=1.34 μ tot = 3/ 5 = 1.34 になる。
近似の裏付け
型 B(極限) | 小さいもの 角度 θ , β ≪ 1 \theta,\beta\ll1 θ , β ≪ 1 と r s / b → 0 r_{\mathrm s}/b\to0 r s / b → 0 、レンズの厚み/距離 → 0 \to0 → 0 (薄いレンズ) | 証明 条件付き(曲がり角を1次の (2.16.23) で置き換え、曲がりをレンズの面に集める。Shapiro 遅延の式も同じ1次の屈折率で、座標によらないのは b b b への依存だけ) | 行き先 近似の数理 第7章
1 I. I. Shapiro (1964), Physical Review Letters 13, 789–791.
2 B. Bertotti, L. Iess, P. Tortora (2003), Nature 425, 374–376.
3 D. Walsh, R. F. Carswell, R. J. Weymann (1979), Nature 279, 381–384.
4 B. Paczyński (1986), The Astrophysical Journal 304, 1–5.
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com