愛因斯坦的曲率公式和 光經過太陽的偏折角度
張海潮
一、1911年愛因斯坦首次計算光經過太陽的偏折角度
愛因斯坦在 1907 年1和 1911 年2 分別發表了兩篇重要的論文, 對廣義相對論進行了初步 的探索。 在論文中, 愛氏指出 : 因為重力場不是慣性系統, 因此光速會受到重力場的影響。 這個 看法推廣了 1905 年 6 月論文3 中所說, 光速在慣性系中是恆定的, 與光源的速度無關。 後來就 以重力是否納入作為狹義和廣義相對論的區分。
在 1907 和 1911 的論文, 愛氏均論證在重力位能是 φ 的情形下, 光速是 c
1 + φ
c2
; (1.1)
式中的 c 指慣性系中的光速, 等於 299,792,458 公尺/秒, φ 是重力位能。 以太陽為例, 若是忽 略其他星球, φ(r) = −GMr , 其中 G 是萬有引力常數, M 是太陽的質量, r 是與太陽 (中心) 的距離。 當 r 以 ∞ 代入, φ = 0, 光速就是慣性系中的光速4 c。
在 1911 論文中, 愛氏對光經過太陽的偏折現象推論如下 : 考慮光在一個平面上行進, 在 此平面的法線方向, 重力位能 φ 與光速皆為常數。
1德文本發表於 《放射學和電子學年鑑》 (Jahrbuch fur Radioakt und Elektronik, 1907), 英譯 : Einstein’s comprehensive 1907 essay on relativity (part III) H. M. Schwartz American Journal of Physics. Vol.45, No.10, October 1977.
中譯 : 紀念愛因斯坦文集第二卷 : 學術論文, 150-209, 台灣凡異出版社 《關於相對性原理和由此得出的結論》
2德文本發表於 《物理年鑑》 (Annalen der Physik, 1911), 英譯 : On the Influence of Gravitation on the Propagation of Light. The Principle of Relativity, 99-108, A collection of Original memoirs on the Special and General Theory of Relativity, Dover Books, 1920. 中譯 : 《關於引力對光傳播路徑的影響》, 同註 1, 台灣凡異出版社, 212-223.
3德文本《論動體的電動力學》 發表於德國 《物理年鑑》 (Annalen der Physik, 1905), 英譯 : On the Electrodynamics of Moving Bodies, 同註 2, The principle of Relativity, 37-65, 中譯 : 《論動體的電動力學》, 同註 1, 台灣凡異出版社, 83-115.
4現在已經通用記號 c 來代表慣性系統在真空中量得的光速, 此處, 當 r = ∞ 時, 代表光速不受重力影響, 因此光速是 c, 其值為 299,792,458 公尺/秒。 (1.1) 式中的光速公式首見於 1907 的論文, 並在 1911 論文中以此計算光經過太陽時偏折的角度。 本文第 一節重現了 1911 年的計算。 然後在 1916 年愛因斯坦修正了 (1.1) 式, 重新做了計算, 請見本文第三節。
圖一
如圖一, 若太陽表面是一平台, 不妨設平面波所在的平面與太陽表面垂直, 平面上的光線 行進時, 始終維持在同一個平面。
圖二
如圖二, 若 x 軸是平面波的波前線, 假設 P1, P2 接近, 例如距離 1 公尺, 則因 P2 比 P1
更接近太陽, 在 P2 的光速 c2 小於在 P1 的光速 c1。 在 ∆t 的時間, 分別以 P2, P1 為圓心, c2∆t, c1∆t 為半徑畫圓, 在 P1, P2 靠近的情形, 假設 c2 到 c1 的變化是線性的, 則此兩圓的 公切線, 就是新的波前線。
圖三說明新的波前線與舊的波前線的夾角 ∆φ 及 c2∆t, c1∆t 的關係。
圖三
因為 P1P2 = 1, 所以
∆φ ≈ sin(∆φ) = (c1−c2)∆t P1P2
≈ ∂c2
∂x∆t, (1.2)
注意到 ∆φ 也是波行進的方向改變的角度, 此一角度必須除以弧長 ∆s = c2∆t 才是曲率, 所 以波前在 P2 位置的曲率 κ 是
κ = lim
∆s→0
∆φ
∆s = 1 c2
∂c2
∂x = ∂
∂x(ln c2). (1.3) 此即愛因斯坦曲率公式, 式中 (由 (1.1) 式知)c2 = c 1 − GMrc2, r 是 P2 到太陽 (中心) 的距離, x 軸是波前線。
我們將在第二節說明愛因斯坦計算的曲率和一般的平面曲線的曲率公式是一樣的。 下面將 以曲率對弧長的積分來計算光偏折的角度。
圖四
如圖四, 太陽在 (0, 0), 光經過太陽偏折很小, 計算時不妨假設波前線永遠與 x 軸平行, 圖 中光的路徑貼近直線 x = R, R 是太陽的半徑。
因為 c2 = c 1 − GMrc2, 則 ln c2 = ln c + ln 1 − GMrc2 , 故曲率 κ= ∂
∂xln c2 = ∂
∂x
ln c + ln
1 −GM rc2
= 1
1 −GMrc2 GM
c2
∂
∂x
−1 r
≈GM c2
∂
∂x
−1 r
. (1.4)
又因為 x2+ y2 = r2 對 x 偏微分有
2x = 2r∂r
∂x, (1.5)
即
∂r
∂x = x
r. (1.6)
將 (1.6) 式代入 (1.4) 式可得 κ ≈ GM
c2 1 r2
∂r
∂x = GM c2
1 r2
x r x=R
= GM c2
R
r3. (1.7)
又因偏折很小, 計算時不妨設弧長 ds = dy, y = R tan θ, dy = R sec2θdθ, r = R cos θ, 則偏折角為 R∞
−∞κds=R∞
−∞κdy, 其中 κ = GM c2
R
r3 = GM c2
cos3θ R2 故 Z ∞
−∞
κdy= GM Rc2
Z π2
−π2
cos3θsec2θdθ = GM Rc2
Z π2
−π2
cos θdθ = GM
Rc2 ·2, (1.8) 將 G, M, R, c 的值代入, 得偏折角是 0.83 弧秒。5
以上是愛因斯坦在 1911 年論文的計算, 到了 1916 年, 他以完整的廣義相對論再計算了 一次, 答案 4GM
Rc2 ≈1.7 弧秒, 是 1911 年計算的兩倍, 我們將在第三節重現 1916 年的計算6。
二、 檢討愛因斯坦的曲率公式
圖五
如圖五, 時間 t = 0 時的波前線是 x 軸, 從 (x, 0) 出發的路徑是 (p(x, t), q(x, t)), p(x, 0) = x, q(x, 0) = 0, 並且 ˙p(x, 0) = 0。 (此處 ˙p = ∂p∂t, ˙q = ∂q∂t)。
現計算曲率。 假設時間 t 的波前線是 y = m(t)x + b(t), 且 ( ˙p(x, t), ˙q(x, t)) 與波前線
5試算
2GM
Rc2 =2 × 6.67 × 10−11m3/kg · s2 × 2 × 1030kg
(0.7 × 109m) × (3 × 108m/s)2 ≈4 × 10−6弧度 (無單位),
其中 4 × 10−6弧度 = 4 × 10−6×180π ×3600 ≈ 0.83 弧秒。 這個數值出現在愛因斯坦 1911 年論文的最後一段, 愛氏的符號 是2kM
c2∆, 其中 k = G, ∆ = R。 另外注意到, c2= c
1 −GMrc2 中的 GM
rc2 ≤GM
Rc2 ≈2 × 10−6, 也就是 c2與 c 的差異最多 差 3 × 108×2 × 10−6≈600m/s。
61916 年愛氏在德國物理年鑑 Annalen der Physik 發表 《廣義相對論的基礎》。 英譯 : The foundation of the General Theory of Relativity, 同註二 The Principle of Relativity, 109-163。 中譯 : 同註一, 台灣凡異出版社 278-334。
y = m(t)x + b(t) 垂直, 則 m(t) = tan θ(t), 兩邊對 t 微分得 m′(t) = sec2θdθ
dt. (2.1)
注意到, 平面曲線曲率 κ 的定義是 dθds, 其中 ds 代表弧長, 因爲 dsdt 代表瞬時速率, 所以 有 ds 和 dt 的關係式 ds =p ˙p2+ ˙q2dt, 故將 (2.1) 式左右同除瞬時速率p ˙p2 + ˙q2 得
m′(t)
p ˙p2+ ˙q2 = sec2θ dθ
p ˙p2+ ˙q2dt = sec2θdθ
ds = sec2θκ(t). (2.2) 現令 t = 0, 觀察到因為 ˙p(x, 0) = 0, 且 sec2θ(0) = sec20 = 1, 則曲率
κ(0) = m′(0)
˙q(x, 0). (2.3)
但
q(x, t) = m(t)p(x, t) + b(t), (2.4) 將 (2.4) 式對 x 偏微分有
qx = mpx. (2.5)
再對 t 偏微分得
qxt = m′px+ mpxt. (2.6) 且 m(0) = 0, p(x, 0) = x, px(x, 0) = 1, 所以當 t = 0 時
m′(0) = qxt(x, 0)
px(x, 0) = qxt(x, 0). (2.7) 因此在 t = 0 時, 曲率
κ(0) = qxt(x, 0)
˙q(x, 0) . (2.8)
回到 (1.4) 式愛因斯坦公式
∂lnp ˙p2+ ˙q2
∂x = ˙pptx+ ˙qqtx
˙p2+ ˙q2 . (2.9) 當 t = 0 時, ˙p = 0, 上式等於 (2.8)
˙qqtx
˙q2 = qtx(x, 0)
˙q(x, 0) = κ(0). (2.10) 愛氏的計算與平面曲線曲率的計算相同。
三、 愛因斯坦在 1916 年再次計算光經過太陽的偏折角度
在 1916 年的論文 (同註 6) 愛因斯坦改進了 1907 年與 1911 年的計算方法, 引入微分 幾何和 (愛因斯坦) 場方程式。 在假設空間中只有一個太陽位於原點的情形, 他得到度規張量 (metric tensor) 的解 c2dτ2 如下 (現在通稱 Schwarzschild 解)
c2dτ2 = c2
1 −2GM rc2
dt2−
1 − 2GM rc2
−1
dr2−r2dθ2−r2sin2θdφ2; (3.1) 式中 G 是萬有引力常數, M 是太陽質量, c 是慣性系統下真空中的光速, r, θ, φ 是球座標, 與 直角座標的關聯是 :
x= r sin θ cos φ y = r sin θ sin φ z = r cos θ
, (3.2)
其中 π
2 −θ 是北緯的緯度, 且滿足 dr2+ r2dθ2+ r2sin2θdφ2 = dx2 + dy2+ dz2, 如圖 :
圖六
Schwarzchild 解 (3.1) 式中用 r, θ, φ 座標是尊重解對原點的球形對稱7。 在 1916 論文
7關於 (3.1) 式的解, 現在通稱 Schwarzschild solution, 它是一個滿足愛因斯坦場方程式的度規張量, 並且要求 : (1) 度規張量與時間無關。
(2) 度規張量球形對稱於原點 (0, 0, 0)。
(3) 在原點之外, 無質量分佈, 即滿足 Rij= 0, Rij是此度規張量的 Ricci 張量。
(4) 當原點的質量趨近於 0 時, 此一度規張量趨近
c2dt2−dr2−r2dθ2−r2sin2θdφ2= c2dt2−dx2−dy2−dz2
請參考 Foster 和 Nightingale 寫的教科書 《A Short Course in General Relativity》 第二版 3.7 The Schwarzschild Solution.
我們應注意, 第一、 當 M 代表太陽質量而 r 大於太陽半徑時, 2GMrc2 ≤ 4 × 10−6 (見註五), 第二、 當 M = 0 時, c2dt2− dr2−r2dθ2−r2sin2θdφ2 稱為 Minkowski flat metric, 它首次出現在 1908 年 9 月 21 日, Minkowski 在德國科隆 (Cologne) 的演講, 講題為 《Space and Time》, 見 The Principle of Relativity p.p 73-91. (同註二), 文中度規張量 c2dτ2= c2dt2−dx2−dy2−dz2 出現在該篇講稿的 p85., 第三、 愛因斯坦在 1916 年論文的前言中提到 : 「下面所要論述的 理論, 是對今天一般稱之為 『相對論』 的理論所做的可能想像得到的最為廣泛的推廣; 為便於區別起見, 以後我把上述 『相對論』 稱為
『狹義相對論』, 並且假定它已為大家所熟悉。 用了 Minkowski 所給予狹義相對論的形式, 相對論的這種推廣就變得很容易; 這位數學 家首先清楚地認識到空間座標和時間座標形式上的等價性, 並把它應用在建立這一理論方面。」 見台灣凡異出版社 p.278 (同註一)。
中, 愛因斯坦所用的符號是
x4 = ct x1 = x x2 = y x3 = z
,
並且用 a 代表 2GMc2 , 我們在本文中, 維持 ct, x, y, z 的記號。 因此 (3.1) 式 Schwarzschild solution 是8 :
c2dτ2 = c2 1 −a
r
dt2−
1 −a r
−1
dr2−r2dθ2−r2sin2θdφ2
≈c2 1 −a
r
dt2 −
1 + a
r
dr2−r2dθ2−r2sin2θdφ2
= c2 1 −a
r
dt2− a
rdr2−dr2−r2dθ2 −r2sin2θdφ2
= c2 1 −a
r
dt2− a
rdr2−dx2−dy2−dz2.
(3.3)
由於
d(r2) = 2rdr = d(x2+ y2+ z2) = 2xdx + 2ydy + 2zdz, (3.4) 所以
dr = x
rdx+ y
rdy+z
rdz. (3.5)
因此
c2dτ2 = c2 1 −a
r
dt2−
a r
x
rdx+y
rdy+z rdz2
−dx2 −dy2−dz2. (3.6) 根據廣義相對論對 c2dτ2 的詮釋, 光的路徑會落在 c2dτ2 = 0 中9。 現在進一步假設平面波在 x − y 平面進行, 如圖
圖七
8根據註五, 當 r 大於太陽半徑時, ar =2GMrc2 < 4 × 10−6, 因此 (3.1) 式的 Schwarzschild metric 與 flat metric 相當接近。
9在 Minkowski flat metric 的情形, c2dt2−dx2−dy2−dz2= 0 表示 (dxdt)2+ (dydt)2+ (dzdt)2= c2, 亦即光的路徑滿足 c2dt2−dx2−dy2−dz2= 0。 在一般的情況, 光的路徑亦必須滿足 c2dτ2= 0, 請參考 Foster 和 Nightingale 寫的教科書 A Short Course in General Relativity 第二版 p.88 的解釋。
因此, 光路徑滿足 : z = 0, x = R 即 dx = 0, 代入上式得 c2dτ2 =
1 − a r
c2dt2− a r
y rdy2
−dy2. (3.7) 此時 c2dτ2 = 0 代表光行路徑, 所以光行滿足
1 − a r
c2dt2 =a
r3y2+ 1
dy2, (3.8)
或
dy dt
2
= c2 1 −a
r
1 + a
r3y2−1
≈c2 1 − a
r
1 − a r3y2
≈c2
1 −2GM
rc2 − 2GM r3c2 y2
.
(3.9)
故
dy dt ≈c
1 −2GM
rc2 − 2GM r3c2 y2
12
≈c
1 − GM
rc2 − GM r3c2y2
. (3.10) 式 (3.10) 中, c2 ≈ dy
dt 就是在太陽重力場下 (地球觀察到) 的光速。 而愛因斯坦曲率公式 :
κ = ∂ln c2
∂x = ∂
∂x
ln
1 −GM rc2 −
GM y2 r3c2
=
1 − GM
rc2 − GM y2 r3c2
−1
∂
∂x
1 − GM
rc2 − GM y2 r3c2
≈ −GM c2
∂(1r)
∂x − GM y2 c2
∂ r13
∂x
= −GM c2
−1 r2 · x
r
−GM y2 c2
−3 r4 ·x
r
= GM x
r3c2 + 3GMy2x r5c2 .
(3.11)
當 x = R = 太陽半徑, 曲率
κ= ∂ln c2
∂x ≈ GM R
r3c2 +3GMRy2
r5c2 . (3.12) 現將 κ 對 ds 積分, 不妨以 dy 代替 ds, 如圖四。 則積分之第一部分
Z ∞
−∞
GM R
r3c2 dy, (3.13)
與圖四之積分相同, 答案是 2GMRc2 。
而積分之第二部分為
Z ∞
−∞
3GMRy2
r5c2 dy. (3.14)
注意到, 此處 r 與 y 有關, 故不可直接積分, 令 y = R tan θ, 則有 dy = R sec2θdθ, 且 cos θ = Rr, 故
y2
r5dy = R2tan2θ · Rsec2θdθ
R5 cos5θ
= 1
R2 sin2θcos θdθ. (3.15) 則積分之第二部分為
Z ∞
−∞
3GMRy2 r5c2 dy=
Z π2
−π2
3GMR c2 · 1
R2 sin2θcos θdθ
= 3GMR c2 · 1
R2 Z θ=π2
θ=−π2
sin2θd(sin θ)
= 3GM
Rc2 · sin3θ 3
θ=π2
θ=−π2
= 2GM Rc2 .
(3.16)
所以 Z ∞
−∞
κds ≈ 2GM
Rc2 +2GM
Rc2 = 4GM
Rc2 ≈1.7 弧秒。 (3.17)
—本文作者為台大數學系退休教授—