では,実際に解を求めてみましょう。3次元空間において,原点対称であるから,極座標表示を用いて考えることにします。線素の表現を得ることを試みましょう。
ds2=−c2dt2+dx2+dy2+dz2 に dt,dr,dθ,dϕ で表されたものを代入していくので,
極座標表示した線素はこれらの 2 次式のみを含みます。
ここで,drdθ,dθdϕ,dϕdr を含む項に関しては消えることがすぐわかります。なぜなら,θ→−θ としたとき
drdθ→−drdθ であり,対称性からこの変化に対して線素の変化があってはならないので,
drdθ−(−drdθ)=0
∴drdθ=0
となるからです。その他もこれと同様です。さらに,dtdθ,dtdϕ,dtdr ついても,時間に重力は無関係なことから, t を変化させて同様の議論をすれば
これらを含む項に関しては消えることがわかります。
よって,
ds2=−Uc2dt2+Vdr2+Wr2dθ2+Xr2sin2θdϕ2
という形で線素を表せるはずです。球対称であることから,U,V,W,X は r のみの関数です。
ここで,W=X であるとすると,θ と ϕ の回転により線素が異なることになってしまうのでおかしいです。よって,
ds2=−Uc2dt2+Vdr2+Wr2(dθ2+sin2θdϕ2)
となります。
いままで,r の取り決めの制約はなかったので,簡単のため Wr2=r2 つまり W=1 となるような r を
基準として系をとりなおすことにします。r の尺度を変えただけであり,たんなる計算の簡略化のためです。
よって,
ds2=−Uc2dt2+Vdr2+r2(dθ2+sin2θdϕ2)
で線素を表すことができます。これを
ds2=−e2νc2dt2+e2λdr2+r2(dθ2+sin2θdϕ2)
と表しても一般性を失わないのでこうおくことにします。
ここで,条件の3つめにより,r→∞ において,線素はMinkowski空間におけるものと同じになる必要があるので,
e2ν,e2λ→1
となる必要があることを注意しておきましょう。
ここから,gμν について,
g00=−e2ν,g11=e2λ,g22=r2,g33=r2sin2θ
であり,これ以外は0を成分に持ちます。また,これより,gμν について,
g00=−e−2ν,g11=e−2λ,g22=r−2,g33=r−2sin−2θ
Ricciテンソルを計算します。
Γνσμ=gμλΓλνσ=21gμλ(∂xσ∂gλν+∂xν∂gλσ−∂xλ∂gνσ)
を利用します。r による偏微分を,プライムをつけて表すことにすれば,
Γ001=g1λΓλ00=g11Γ100=21g11(∂x0∂g10+∂x0∂g10−∂x1∂g00)=ν′⋅e2ν−2λ
Γ111=g11Γ111=21g11(∂x1∂g11)=λ′
などと計算していくと,
Γ001Γ111Γ221Γ331Γ100Γ122Γ233Γ332=ν′e2ν−2λ=λ′=−re−2λ=−re−2λsin2θ=ν′=Γ133=r−1=sinθcosθ=−sinθcosθ
これら以外の成分は0です。ただし,Christoffelの記号の下2つの添字については対称であることに留意してください。つまり,
Γ100=Γ010
等が成立します。
さらにこれらを用いて,Ricciテンソル
を計算すれば,
R00R11R22R33=(ν′′−λ′ν′+ν′2+r2ν′)e2ν−2λ=−ν′′+λ′ν′−ν′2+r2ν′=(−1−rν′+rλ′)e−2λ+1=R22sin2θ
となります。これら以外の成分は0です。
さて,Einsteinの重力場方程式により,Rμν=0 が要求されます。
{R00=0R11=0
⟹⎩⎨⎧ν′′−λ′ν′+ν′2+r2ν′=0−ν′′+λ′ν′−ν′2+r2ν′=0
⟹λ′+ν′=0
⟹(λ+ν)′=0(1)
これより,λ+ν は定数です。ここで,式:
e2ν,e2λ→1
より,r→∞ では,λ,ν→0。
よって,
λ+ν=0(2)
さらに,
R22=0
より,
(−1−rν′+rλ′)e−2λ+1=0
ここで,式 (1),(2) より,
(1+2rν′)e2ν(re2ν)′re2ν=1=1=r−2m(3)
ここで,−2m は積分定数です。ここまでをまとめると,
Rμν=0⟺⎩⎨⎧R00=0R11=0R22=0R33=0⇒{re2ν=r−2mλ+ν=0
ちなみに,逆に最右辺を仮定しても,Rμν=0 を示すことができます。つまり,これらは全て同値であることを証明できます。
これらの式より,計量テンソルを求められます。式 (3) により,
g00=−e2ν=−(1−r2m)
また,
g11=e2λ=e−2ν=(1−r2m)−1
これらにより,Schwarzschild解:
ds2=−(1−r2m)c2dt2+(1−r2m)−1dr2+r2(dθ2+sin2θdϕ2)
を得ます。