ブラック2

第6章 ブラックホ−ルへの道2

逆説の相対性理論 第3部に戻る

第5章 アインシュタインの方程式に戻る

第7章 シュバルトシルトの解(怪?)へ行く


アインシュタインのテンソルが出たところでもう一度ブラックホ−ルへ挑戦します。前回(4章)はアフィン係数や曲率テンソルの所でポショちゃったのですが夢よモウ一度....と言うわけです。

ただアインシュタインのテンソルそのものを理解しながら進むのは少々骨が折れます。とりあえずブラックホ−ルの麓(フモト)から6合目(?)あたりまで まだ雲の下、シュバルツシルトの外部解まで登って行きます。

六根清浄


6−1)夜空はなぜ暗い? 簡単そうで実は困った問題

        本当の目的 重力ポテンシャルΦ

6−2)アインシュタインの方程式

        a)定理

         b)エネルギ−運動量テンソルμν

6−3)3次元空間、2次元空間の計量テンソルの解釈

      6-3-1)3次元空間の計算

        a)曲座標の計量テンソル  時間を無視した3次元空間での考察!

        b)曲座標のクリストッフェルのカッコ

        c)曲座標のリ−マンクリストッフェルのテンソル

      6-3-2)2次元空間の計算

        a)3次元球面上の2次元平面

        b)アインシュタインのテンソルの式を適用する。

        c)イロイロな無重力

6−4)アインシュタインのテンソル/方程式 ってナンだ?

       a)アインシュタインのテンソルの式  曲がる条件

                       シュバルツシルト 見参!

       b)クリストッフェルのカッコ

      c)リ−マンクリストッフェルのテンソル R10,01

     d)リ−マンクリストッフェルのテンソル R01,10

     e)どう曲がってる

6−5)この章の終りに...実は口上

          ちょっと一服 「漆黒洞


6−1)夜空はなぜ暗い? 簡単そうで実は困った問題

夜空はなぜ暗い? この答えは小学生も熟年の男女も大体同じである。

太陽が沈むから

しかし夜空には星が輝いている。太陽の数百倍の星がゴロゴロ(?)と何万個も輝いている。それなのになぜ暗い?

それは地球から う〜んと遠いので 光が弱まるから

しかし 星(恒星)がう〜んと数があれば星(恒星)がう〜んと遠くにあろうとやはり夜は明るくなる。つまり夜が暗くなるためには星(恒星)がう〜んと数があってもその効果がなくなるほどう〜んと遠くに配置されなければならない。つまりふたつのう〜んう〜んの兼ね合いで夜空は暗くなる!。トイレじゃないぞ!(^o^)

そこでこの兼ね合いを計算する...と言っても別に難しくはない。まず半径rの球の表面積は

πr     ;球の表面積

である。この式は電磁気学でもちょろちょろ出て来る。

球の中心に地球を置き、rを数百光年くらいにする。するとその内側に太陽クラスの恒星がいくつか入る。その密度をρとする。すると半径rと半径r+drの球のカラ(厚さdr)の中には

ρπrdr 個

太陽クラスの恒星が入ることになる。

太陽クラスの恒星ひとつからの地球に到達する光は1/rに比例して弱まる。よって ρπrdr 個の光は地球に到達する頃は合計すると

ρπrdr×(1/r)=ρπdrに比例している。

よって地球から半径rの球の中にある太陽クラスの恒星からの全ての光量は

ρπdr=ρπrに比例することになる。

式 ρπr を解釈すると「rを無限に大きくすると光量も無限に大きくなる。」つまり宇宙が無限に広がっており、その宇宙に太陽クラスの恒星が(密度が薄くとも)無限に有れば、地球の夜は無限に明るくなる!。(^^)

言葉を変えて言えば、恒星が遠ざかるために光が弱まる効果う〜んは恒星の数の増加で光が強くなる効果う〜んによって打ち消され、さらに強まるわけである。う〜んう〜んの兼ね合いは矛盾する。

この計算に具体的な数値を入れるともっと悲惨になるらしいが...ともかく夜が暗いのは地球に到達する光は1/rに比例して弱まる以上に弱まるためである。宇宙に限界があるとか、暗黒物質があるとかも、言えなくはないがメジャ−ではない。

具体的に言えば、遠くの星ほど波長が長くなる(=弱くなる)効果が加わるため、夜空は無限に明るくならない。遠くの星ほど波長が長くなるのは所謂ドップラ−効果、そのドップラ−効果がなぜ起きるかは宇宙が膨張しているから...よくよく言われているビッグバン理論である。夜が暗いのは宇宙が膨張している証拠である。

本当の目的 重力ポテンシャルΦ

この辺の話はこのHPを読んでいる読者にとって当たり前かもしれない。我々があえて一般相対性理論の真っ最中にこの話をするのは

ρπdr=ρπr と 式刄ウ=πGρ(X) ρ(X);物質の質量密度

の類似性に注目して 式刄ウ=πGρ(X) の意味を汲み取って欲しいからである。Φは重力ポテンシャルである。そろそろトイレを出よう!。(^ヘ^)トイレじゃないぞ!

式刄ウ=πGρ(X) を言葉で説明しようとすると結構苦しい。まず数学的な類推から類推する。

ρπdr=ρπr=SとすればdS/dr=ρπである。

 dS/drは地球を原点とした極座標のrで全光量を微分したものである。

この微分をSの増加率儡とすれば 儡=ρπである。

式刄ウ=πGρ(X)はその延長線で考えれば良い。

重力も光も距離の2乗に反比例して弱まる。重力の場合 ニュ−トンの万有引力の公式

F=GmM/r

G=6.67±0.005×10-8   cm3-1sec-2

=6.67±0.005×10-11  Nm2kg-2=m3kg-1S-2 

となるから恒星の質量をMとすると GM/rを質量Mの恒星が作る重力場Φと考えて良い。

Φ=GM/r

半径rと半径r+drの球のカラ(厚さdr)の中には

ρ(X)πrdr 個    ρ(X);物質の質量密度?

恒星が入る。これを上記Mとして代入すれば

Φ=GM/r

=G(ρ(X)πrdr)/r

Gρ(X)πdrになるので

増加率刄ウ=πGρ(X)になる。

ここでρ(X);物質の質量密度 としたのだが、実際には1辺数百光年のスペ−スに入っている恒星の全質量をスペ−スで割ったものだから(気分として)水や空気のような質量密度ではない。ただ物理的、数学的にはタダの質量密度である。ρ(X)のXは位置により密度が変わることを意味する。

 刄ウ=πGρ(X) 自体は重力ポテンシャルと宇宙的質量密度(?)ρ(X)の関係にしか過ぎない。これを言葉で説明しようとすると稀有壮大なオオボラになりそうである。

稀有壮大なオオボラは未だに成功しないのだが、大体こうなる。

地球を宇宙にポツンと浮かんでいる点と考え、それを原点とする座標を考える。そして適当な恒星が適当な存在濃度ρ(X)で地球を取り巻くとした時、地球は重力ポテンシャルΦの環境に置かれる。

重力ポテンシャルΦは地球を取り巻く恒星が作るポテンシャルでとりあえずスカラ−である。スカラ−であるからポテンシャルそのものを地球をどこかに引っ張る引力(これは方向と量が要るからベクトル)だとしてはならない。Φがいくら大きくともそれは引力が強いことを意味しない。

ただし、地球を引っ張る引力は地球を取り巻く恒星が作るポテンシャルを通して形成されることまで否定してはいない。その引力は−(gradΦ)であたえられる。mは地球の質量である。もし地球を取り巻く恒星が偏った分布をしているならそっちの方向に引力−(gradΦ)で引っ張られることになる。

−−−−−−−−−−−−−−−−−

このΦをニュ−トンの重力ポテンシャルΦと言う。ニュ−トン自身が作ったかどうかは知らないが、ニュ−トン力学の範囲では重力ポテンシャルΦはスカラ−でOKなのである。

このΦ(ニュ−トンの重力ポテンシャル)をもう少し考える。おなじみの位置エネルギ−

E=mgh  g;9.8m/s2

E=mΦ  Φ=ghと書替えても良い。

Φはまさに重力ポテンシャルである。地球によって作られたΦである。この場合hは「物体と地球の中心までの距離を意味しない」。Φの大きさは引力を意味しないからである。Φによる引力は−m(gradΦ)で与えられ、この場合h方向のみであるから

gradΦ=∂Φ/∂h=gであり、引力は−mgになる。

この重力ポテンシャルΦは後で出て来る。引力とは(gradΦ)である。

ただ この議論は後から気付いたのだが、質量Mの恒星だけが作る重力場をΦとした場合、その恒星による引力F=GmM/r=−(gradΦ)ならば

gradΦ=−GM/r Φ=GM/rとなる。

最初の定義 Φ=GM/r と異なる?。確かに!。...分かるような分からないような謎だが、今後質量Mの恒星だけが作る重力場Φは Φ=GM/rとする。(この項 未考慮)

6−2)アインシュタインの方程式

        a)定理

  アインシュタインのテンソルμνには次のような定理がある。

μμν=0 で Tμνはgμν,gμνの1階微分、gμνの2階微分を成分とするなら

μν=c1μν+cμνである。

これを証明するのがまた面倒なので来年まで先送りすることとして(^0^)、この定理を使うとアインシュタインの方程式を導ける。

μν=c1μν+cμν 

μν+(c/c1)gμν=(1/c1)Tμν ;*(1/c1)

μν+λgμν=κTμν  ;c/c1=λ (1/c1)=κ とした。

ここでアインシュタインのテンソルをGμν=Rμν−(1/2)gμνRとすれば

μν−(1/2)gμνR+λgμν=κTμνである。

このλgμνがアインシュタインをして「生涯最大の過ち」と言わしめた宇宙項である。このλが0なら宇宙は膨張し、マイナスなら収縮する。適当な値であれば定常宇宙になる。

実際には宇宙は膨張しているのでλ=0で良いと言うのがメジャ−だが、数学的には宇宙項が入っても良いわけである。ともかく教科書では

μν−(1/2)gμν=κTμνをアインシュタインの方程式としている。

        b)エネルギ−運動量テンソルμν

上記アインシュタインの方程式に登場したμν,gμνの1階微分、gμνの2階微分を成分とするμνをエネルギ−運動量テンソルと言う。このμνの成分が6−1)のΦと結びつく。

ここで普通なら

裳華房;「一般相対性及び重力の理論」(教科書)  Rμν−(1/2)gμνR=κμν

が出たわけであるからブラックホ−ルへ行く道すがらアインシュタインの方程式を考えれば良いのだが、「アインの細道第3部 イロイロあるアインシュタインの方程式」で列挙した下記の式が気になる。

学研;「最新アインシュタイン論」  Gμν=8πGTμν

図解雑学;「重力と一般相対性理論」  Rμν−(1/2)gμνR=(8πG/C)Tμν

通商産業研究社;「一般相対性理論の直観的方法」  Rμν−(1/2)gμνRμν

これらの式はエネルギ−運動量テンソルの係数の問題であり、エネルギ−運動量テンソルをどう取り扱っているかの問題である。幸いなことに係数(8πG/C)を呼び出す呪法(?)があるのでそれで考えて行く....つもりだったがこの呪法(?)は九十九(ツクモ)神まで召還するらしい。(^ι^)夢枕獏の「陽陰師」は面白い!

簡単に言えば重力場を作る物体が静止していてしかも非常に弱い場合、つまり重力場そのものも静止しておりほぼニュ−トンの重力理論に合う場合で、しかも重力場を作る物体の流れが完全流体に近い場合においては、多くの近似が成り立つので(8πG/C)が出て来るわけである。 いま少し 短兵急に結論を求めないで、アインシュタインのテンソルや方程式を眺めたほうが良いようである。H13.10.10 ハヤイハナシ 負けた!

我々はアインシュタインの方程式を数学的に(機械的に)導き出しただけであり、物理的な意味は考慮してない。少々遠回りだが、アインシュタインのテンソルと方程式を考察してからブラックホ−ルへ向かう。

 

6−3)3次元空間、2次元空間の計量テンソルの解釈

第5章を俗っぽく解釈する。

      6-3-1)3次元空間の計算

        a)曲座標の計量テンソル  時間を無視した3次元空間での考察!

       曲座標を下記に定める。 ψ;X-Y平面に対する仰角 θ;X-Y平面でのX軸に対する回転角

     X=Rsin(θ)cos(ψ)   Y=Rsin(θ)sin(ψ)   Z=Rcos(θ)

一応 S^2=X^2+Y^2+Z^2

^2sin(θ)^2cos(ψ)^2+^2sinθ)^2sin(ψ)^2+^2cos(θ)^2

^2{sin(θ)^2(cos(ψ)^2+sin(ψ)^2)+cosθ)^2}

^2{sin(θ)^2+cos(θ)^2}=^2

つまり S^2=^2 である。このS^2=X^2+Y^2+Z^2は気取って

S^2=ημνXμXν  μ,ν=1,2,3 X1X,XY,XZ

η11=1 η22=1 η33=1 ημν=0 μν と書けるし

dS^2=ημνdXμdXν  ;微小変化の式 とも書ける。

ところでこれは同じものを直行座標(X,Y,Z)で見てるか、曲座標(R,θ,ψ)で見てるかの違いである。

よって微小変化の式もR,θ,ψで表され

dS^2=gμνdXμdXν  μ,ν=1,2,3  X=R、X=θ、X=ψ

の形式にすれば g11=1   g22=(X1^2  g33=(X1sin(θ)^2  gμν=0 μν

となる。ピンとこないかもしれないが実は単純なことだ。

曲座標(R,θ,ψ)で表されたベクトルをR+dR,θ+dθ,ψ+dψ)とするベクトルにすることを考える。

最初θだけ-dθ動かした時、ベクトルの先が動く距離をS(dθ)とする。同じようにψだけdψ動かした距離をS(dψ)、RだけdR動かした距離をS(dR)とするなら、もとのベクトルをdθ,dψ,dRだけ動かした距離dSは

dS^2=S(dR)^2+S(dθ)^2+S(dψ)^2となる。

つまり 立体的な幾何学でしかない(下図参照)

動く距離 S(dR),S(dθ),S(dψ) のうちで一番考えやすいのがS(dR)=dR 次がS(dθ)=−Rdθだろう。

なお最初にdRだけ動かして次にdθ動かすと、本当はdθの変化で(R+dR)dθ=Rdθ+dRdθとすべきだろうが、dRdθ=0として無視することにする。

S(dψ)は注意が要る。考え方を図に示す。S(dψ)=sin(θ)dψになる。

よって

dS^2=S(dR)^2+S(dθ)^2+S(dψ)^2

=(dR)^2+(Rdθ)^2+(Rsin(θ)dψ)^2

テンソル風に書きなおすと

dS^2=gμνdXμdXν  μ,ν=1,2,3  X=R、X=θ、X=ψ

g11=1   g22=(X1^2  g33=(X1sin(X2)^2  他は0

一方直行座標(X,Y,Z)でのテンソル風の式は

dS^2=ημνdXμdXν  μ,ν=1,2,3 X1X,XY,XZ

 η11η22η331 他は0

ここで我々は初めてη**以外の計量テンソルを手に入れた。同じものでも測定方法で計量テンソルは異なる。この計量テンソルを元にしたリ−マンクリストッフェルのテンソルRαβ,μυは0になる。ハズ.....

        b)曲座標のクリストッフェルのカッコ

とりあえず、gμν gμν=E  となるgμν

g11=1   g22=1/(X1^2  g33=1/(X1sin(X2)^2   他は0

とする。要するに g11g11=1 g22g22=1  g33g33=1 なのである。

 

するとクリストッフェルのカッコΓσλν={λσν}=(∂λμν+∂νλμ−∂μνλσμ/2 が計算できる。 第5章 アインシュタインの方程式5−5) 参照

計算と言っても左辺は添え字νσλμで右辺はνσλだ。μをどうする? だがσμ ならσμ=0だ。だから実際に計算する価値のあるクリストッフェルのカッコはμ=σとしての

Γσλν={λσν}=(∂λσν+∂νλσ−∂σνλσσ/2 だけである。

ここでσ=1とする。するとσν,gλσは0か1であるから偏微分すれば0である。よって

Γ1λν−∂1νλ/2 だけが残る。  ;σσ11=1 

νλについては、ν=λ以外は0である。よってΓ1λνΓ122Γ133以外は0である。

Γ122−∂122/2=−(1/2)∂22X1-X1  :λ=ν=2 g22=(X1^2

 Γ133−∂133/2=−(1/2)∂33X1-X1sin(X2)^2  :λ=ν=3 g33=(X1sin(X2)^2

またσ=2とする。この場合、一端下記のように書いておき、ν,λを代入する前に0になるか否か判断すると早い。

Γ2λν={λ2ν}=(∂λ2ν+∂νλ2−∂2νλ22/2

ν=λ=1では青字が0になる。(下記)よってΓ211=0である。

Γ211(∂λ2ν+∂νλ22νλ22/2=0 

ν=1 λ=2では青字が0になり、式が1/3に短くなる。

Γ212(∂λ2ν+∂νλ2−∂2νλ22/2 

(∂122)g22/2=(1/2)(∂(X1^2/X1)g22

=X122=1/X1  ;g22=1/(X1^2

以下、同じ事の繰り返しである。面倒なので結論だけ列挙する。

 Γ213=0   Γ221=1/X1 ;上式Γ212と同じ つまりΓ221Γ212=1/X1

Γ222=0  Γ223=0  

Γ231=0   Γ232=0

Γ233(∂λ2ν+∂νλ2−∂2νλ22/2  ;ν=λ=3では青字が0

=−(∂23322/2=−(1/2)(∂g33/X2)g22  ;g33=(X1sin(X2)^2

sin(X2)cos(X2) 

またσ=3とする。

Γ3λν={λ3ν}=(∂λ3ν+∂νλ3−∂3νλ33/2

Γ311=0  Γ312=0  

Γ313(∂λ3ν+∂νλ3−∂3νλ33/2 ;ν=1 λ=3では青字が0

(∂13333/2=(1/2)(∂g33/X1)g33

 X1(sin(X2)^233  ;g33=(X1sin(X2)^2

=1/X1  ;g33=1/(X1sin(X2)^2 

Γ321=0  Γ322=0

Γ323(∂λ3ν+∂νλ3−∂3νλ33/2 ;ν=2 λ=3では青字が0

=(∂23333/2=(1/2)(∂g33/∂X2)g33

X1^2sin(X2)cos(X2)33   ;g33=(X1sin(X2)^2

     cos(X2)/sin(X2)  ;g33=1/(X1sin(X2)^2

Γ3311/X1  Γ332cos(X2)/sin(X2)  Γ333=0 

結局、曲座標における0でないリ−マンクリストッフェルのカッコは

Γ122=-X1   Γ133-X1sin(X2)^2  

Γ212=Γ2211/X1   Γ233sin(X2)cos(X2)

Γ313Γ3311/X1   Γ323Γ332cos(X2)/sin(X2) 

でもこの計算は合ってるのカイナ!→  ラシイ!

 

        c)曲座標のリ−マンクリストッフェルのテンソル

リ−マンクリストッフェルのテンソルRαβ,μυは(第5章 5−6)参照 ゲェッ〜

Rαβ,μυ=∂μβαυ−∂υβαμβλυ}{λαμβλμ}{λαυとなる。5-6-1)

この式も左辺はα,β,μ,υの4文字なのに右辺はα,β,μ,υ,λの5文字で落ち着かない。

教科書に従いα=1,β=3,μ=1,υ=3とすると

R13,13=∂1313−∂33113λ3}{λ113λ1}{λ13

=∂13133λ3}{λ113λ1}{λ13}  ;{311}=Γ131=0

−−−−−−−−−−−−−−−−−−−−

ここで初心(?)に戻り「Xυ,αの仮面をかぶっている式は一族全員が参加している。Συ=0が省略されている。」とする。<第3章 テンソルの散策/3−1)反変ベクトルと共変ベクトル/a)dXυ=(Xυα)dxαの解釈> 参照

つまり 

3λ3}{λ113λ1}{λ13}=λ=13[{3λ3}{λ113λ1}{λ13}]

333}{311+{323211+{313111311}{113321}{213}−{331}{313

=−{331}{313}  ;青字が0

=−Γ331Γ133 

=−(1/X1)(-X1sin(X2)^2)=(sin(X2)^2

となる。(小声)っと考えるべきだろうなぁ〜

よって

R13,13=∂13133λ3}{λ113λ1}{λ13

∂Γ133X1(sin(X2)^2

(-X1sin(X2)^2)X1+(sin(X2)^2

-(sin(X2)^2+sin(X2)^2=0  となる。

−−−−−−−−−−−−−−−−−−−−−

なおμ=υならRαβ,μμ=0 λはなんでも良い。下記で打ち消し合う!。

Rαβ,μμμβαμ−∂μβαμβλμ}{λαμβλμ}{λαμ

μβαμ−∂μβαμβλμ}{λαμβλμ}{λαμ

このようにユ−クリッド空間ではRαβ,μμ=0になる。(原語)(^0^)oh〜

なお ここで注意したいのだがRαβ,μμの変わりに

Rαβ,μυ=(1/2){∂μβυα−∂υβμα+∂υαμβ−∂μαυβ} 5-6-6)

を計算してRαβ,μμααRαβ,μυとするとまったく別のものになる。これはΓ(アフィン係数)を0としてRαβ,μυを導いているためである。第5章 アインシュタインの方程式 5−6)参照

     6-3-2)2次元空間の計算

        a)3次元球面上の2次元平面

早い話、地球の表面における平面の話である。幾何学の「三角形の内角の和は180゜である」というのは普通のサイズなら地表にかいた三角形でも成り立つ。だが三角形の中にロシアや中国がすっぽり入るほど大きい三角形なら「内角の和は180゜以上」になる。

曲がった2次元平面と言っても良いだろう。

曲座標を下記に定める。 ψ;X-Y平面に対する仰角 θ;X-Y平面でのX軸に対する回転角

     X=Rsin(θ)cos(ψ)   Y=Rsin(θ)sin(ψ)   Z=Rcos(θ)

この曲座標の原点を地球の中心に置き、Rを地球の半径とする。すると地上の任意の点はθとψで表される。これは2次元である。(ただし 曲がった!)

この平面でも

dS^2=μνdXμdXν  μ,ν=1,2 X1=θ,X=ψ

となり、g11=R^2  g22=(Rsin(X1)^2  他は0

および g11=1/R^2  g22=1/(Rsin(X1)^2  他は0  となる。

するとクリストッフェルのカッコΓσλν={λσν}=(∂λμν+∂νλμ−∂μνλσμ/2 が計算できる。 第5章 アインシュタインの方程式5−5) 参照

例によってσμ ならσμ=0だ。だから実際に計算する価値のあるクリストッフェルのカッコはμ=σとして

Γσλν(∂λσν+∂νλσ−∂σνλσσ/2となり

Γ1λν(∂λ1ν+∂νλ1−∂1νλ11/2

以下 青文字は0を示す

Γ111(11111111111/2=0

Γ112(∂112211−∂11211/2=0

Γ121Γ112=0

Γ122(∂212+∂221−∂12211/2

−(∂22/X111/2  

−sin(X1)cos(X1) g22=(Rsin(X1)^2  g11=1/R^2

Γ2λν(∂λ2ν+∂νλ2−∂2νλ22/2となり

Γ211(∂121+∂ν1221122/2=0

Γ212(∂122+∂212−∂22122/2

(∂22/X122/2

cos(X1)/sin(X1)=Γ221  g22=(Rsin(X1)^2 g22=1/(Rsin(X1)^2

Γ222(222222−∂22222/2=0

結局   Γ122−sin(X1)cos(X1)  Γ212Γ221cos(X1)/sin(X1)  他は0 となる。

Γσλνも{λσν}もΓ122もクリストッフェルのカッコと言う。(ラシイ!)....でそのクリスが出たので例によってリ−マンクリストッフェルのテンソルRαβ,μυを求めるわけだが「μ=υならRαβ,μμ=0 λはなんでも良い。」わけだから 0でないリ−クリのテンソル(?)は

Rαβ,12 ,Rαβ,21のみである。

Rαβ,121βα2−∂2βα1βλ2}{λα1βλ1}{λα2}とすれば

計算すべきRαβ,12R11,12 R12,12 R21,12 R22,12 

R11,121112−∂21111λ2}{λ111λ1}{λ12 ;青字が0

λ1λ2}{λ111λ1}{λ12

112}{111111}{112+{122211121212}=0

つまりR11,12=0

R12,121212−∂22112λ2}{λ112λ1}{λ12

1212212111211}{112222}{211221}{212} ;同時に書く

1212221}{212

∂(sin(X1)cos(X1))/X1{cos(X1)/sin(X1)}{−sin(X1)cos(X1)}

=−(cos(X1))^2+(sin(X1))^2+(cos(X1))^2

(sin(X1))^2

面倒になった! 結果だけ書くと

 R12,12(sin(X1))^2  R12,21−(sin(X1))^2  他は0

それで

R12,21−(sin(X1))^2=R22;リッチのテンソル

g2222  ;R22によるスカラ−曲率

(sin(X1))^2/(Rsin(X1))^2   ;g22=1/(Rsin(X1)^2

=−1/R^2

一方

R12,12g11R12,12R^2(sin(X1))^2 ;g11=R^2

R21,12=−R12,12−R^2(sin(X1))^2 

R21,12=g22R21,12 ;g22=1/(Rsin(X1)^2

    これも11;リッチのテンソル 

g1111  ;R11によるスカラ−曲率

=−1/R^2  ;g11=1/R^2

ただしスカラ−曲率自体は

gμνμνμνgμνμν

g1111g2222

=−2/R^2 である。

この関係式=-2/R^2より をスカラ−曲率と呼ぶ どうだ! 分かったかとなるのだが........今ひとつ!

        b)アインシュタインのテンソルの式を適用する。

まぁ ともかく...これは3次元的に見たら曲がっている2次元空間の話である。その空間μν

μν; g11=R^2  g22=(Rsin(X1)^2  他は0

と 書かれ、リッチのテンソル、スカラ−曲率

Rμν R22−(sin(X1))^2  R11−1  他は0  =−2/R^2

である。これらが

Rμν−(1/2)gμν=0  ;アインシュタインのテンソルの式

を満足していることを示そう。どうも教科書を除いてほとんどの本が添え字を下につけている。多分上記の形にも導けるのであろう。(^^)   Rμν,gμνを行列として書くと

となる。よって

行列=−Rμν=gμν/R^2

書き変えると

Rμν+gμν/R^2=0

Rμν−(1/2)gμν=0  ;=−2/R^2 より

−−−−−−−−−−−−−−−−−−−−−

       c)イロイロな無重力

 これは曲がった2次元空間の話の話である。結果としてそれはアインシュタインのテンソルGμν=0を満足する。

  3次元的に見たら曲がっている2次元空間の話??   簡単に考えよう。我々は地球の表面(=曲がった2次元空間)に縛り付けられているアワレ(?)な存在である。そこではミクロ的(実生活上)ではユ-クリッド幾何学が成り立つ。土地を買うとか売る場合、その面積は四角四面の土地ならば縦×横(m^2)として計算すれば良い。どこの不動産でもそれで良い。だけどその面積がベラボ−に大きくなり中国やロシアがすっぽり入るくらいになると×横(m^2)としての計算では成り立たない。「三角形の内角の和は180゜」ではないからであり、曲がった2次元空間だからである。

ではこれを4次元的に見たら曲がっている3次元空間の話に書きなおせるか?

これはミクロ的に見たら問題ないがマクロ的(総括的)に見たら問題となる事項があるかと言う問いでもある。その実例を所謂「自由落下のエレベ−タ−」では次のように記述する。

最初 エレベ−タ−内にふたつのリンゴを適当な距離において自由落下をさせる。このリンゴは地球の中心に近づくにつれ、徐々に接近するはずである。なぜなら各々のリンゴは地球の中心に向かって落ちて行くだけであるから。

これはリンゴ間に引力が働いたためではなく、空間がそうなっているから接近すると言うわけだ。

そうすると(無重力)のアインシュタインの方程式もイロイロあることになる。「自由落下のエレベ−タ−」での無重力も無重力である。だけど本当の無重力なら「リンゴが徐々に接近する」ことはない。「自由落下のエレベ−タ−」での無重力は時間的空間的に局所の無重力である。

そんなわけでμν−(1/2)gμν=0は完全なる(?)無重力ではなくバ-チャルな無重力の式である。完全なる(?)無重力ではRμν=0,=0だから0−(1/2)gμν×0=0となる。「一切が空?」とでも言おうか? ともかく

我々の空間は4次元的に見たら曲がっている3次元空間なのである。

我々が概念する宇宙はバ-チャルな無重力はあるが本当の無重力はない。宇宙にそっと置いた物体はそのうちどこかの星に落っこちるか、楕円の衛星軌道にハマルしかないわけである。(小声)ココまで言いきって良いのか?

そう言うわけだ。..お待たせしました。いよいよ ブラックホ−ルの本質に突入しよう。

 

6−4)アインシュタインのテンソル/方程式 ってナンだ?

まずアインシュタインのテンソルの式と アインシュタインの方程式を並べる。

Rμν−(1/2)gμν=0  ;アインシュタインのテンソルの式

Rμν−(1/2)gμν=κTμν  ;アインシュタインの方程式

アインシュタインのテンソルの式はμν=0でのアインシュタインの方程式であるからアインシュタインのテンソルの式はμν=0(無重力)でのアインシュタインの方程式であると言えそうである。←早口言葉?

そしてこの式Rμν−(1/2)gμν=0自体は計量テンソルgμνの条件である。スカラ−曲率の空間の「自由落下のエレベ−タ−」での無重力(=バ-チャルな無重力)での式である。

それでまずアインシュタインのテンソルの式から行く!

        a)アインシュタインのテンソルの式  曲がる条件

シュバルツシルト 見参!

Rμν−(1/2)gμν=0  ;アインシュタインのテンソルの式

これに  a)曲座標の計量テンソル  時間を無視した3次元空間での考察! の式

dS^2=S(dR)^2+S(dθ)^2+S(dψ)^2

=(dR)^2+(Rdθ)^2+(Rsin(θ)dψ)^2

にCTを考慮する。つまり

dS^2=−(CdT)^2+S(dR)^2+S(dθ)^2+S(dψ)^2

=−(CdT)^2+(dR)^2+(Rdθ)^2+(Rsin(θ)dψ)^2

テンソル風に書きなおすと

dS^2=gμνdXμdXν  μ,ν=0,1,2,3  X=CT、X=R、X=θ、X=ψ

g00=−1 g11=1   g22=(X1^2  g33=(X1sin(X2)^2  他は0

そして、リ−マンクリストッフェルのテンソル→リッチのテンソル→スカラ−曲率と追い詰めるわけであるが、残念ながらこのままではリ−マンクリストッフェルのテンソルRαβ,μυは全て0になる。その結果スカラ−曲率も0になる。つまり 曲がらない。

ではこれが曲がるためにはどうすれば良いか? 「a)3次元球面上の2次元平面」  「b)アインシュタインのテンソルの式を適用する。」にヒントが隠されている。

μν; g11=R^2  g22=(Rsin(X1)^2  他は0

Rμν R22−(sin(X1))^2  R11−1  他は0

この式からRμν/gμν=−(1/R^2)となる。これが=−2/R^2≠0に繋がるわけだ。

...と言うことはRμν−(1/2)gμν=0を満足することはRμν/gμν=(1/2)≠0 を満足させることである。ただしgμν=0ならRμν=0でこれは計算しない。

そうなるとR00,R11ともR00/g00=R11/g11≠0であるようにg00g11 を決めなければならない。

R00から検討する。

R00=Rα0,0α=∂00αα−∂α0α00λα}{λα00λ0}{λαα

上式にα=0を代入すれば簡単に計算できるが、R00,00=0である。だからαは1以上である。ところで上の式は全て0}と書ける項を含んでいる。これらの{0}が全部0であれば、R00は0になるからいくつかは0ではない。

ところでクリストッフェルのカッコ Γσλν={λσν}=(∂λμν+∂νλμ−∂μνλσμ/2  より

0}=(∂0μb+∂b−∂μb0aμ/2 

となるがこの{0}が0でないためにはaμ≠0 つまりaμでなければならないから

0}=(∂0ab+∂b0a−∂b0aa/2となる。

ここで0abに注目する。これは計量テンソルの成分をCTで微分することである。もしこれが0でないなら計量テンソルが時間的に変化することである。逆に0**=0なら時間経過によって空間が変化しないことである。今 どんな条件だったか? 空間がニョロニョロ変わるなら無重力どころではない。よって今の条件下では0ab=0である。....なら

0}=(0ab+∂b0a−∂b0aa/2

(∂b0a−∂b0aa/2であり、

この0}が0でないためには、a,bが同時に0であってはダメだが、aかbどちらかはかならず0でなければならないことになる。

言葉を変えて言えば、00X(=CT)を含んではならないが、X以外のものはなにかひとつはかならず含まなければならない

同様にして

R11Rα1,1α=∂11αα−∂α1α11λα}{λα11λ1}{λαα

各項は{1**}がかならず入っている。よってこれらの{1**}が全部0であれば、R11は0になるからいくつかは0ではない。

1}=Γ1b(∂1ab+∂b1a−∂b1aa/2 になる。

=0の場合10}=Γ01b(∂10b+∂b100b100/2 であるからb=0以外は0である。

Γ10={100}={001}=(∂10000/2  他のa=0のΓ0**は0

=0の場合10}=Γ10(∂1a001a−∂01aa/2であるからa=0以外は0である。

Γ10={100}={001}=(∂10000/2  他のb=0のΓ0**は0

要するにaかbが0なら、a=b=0以外の1}は0である。これらは少なくとも00X1(=R)を含んでいなければならないことになる。

時間経過によって空間が変化しない条件とは00X(=CT)を含んではならない かつX1(=R)を含んでいなければならない。

g11 も同じような条件がでるはずだが、うまく行かない。トホホ....

*****************

教科書に従ってこのg00

g00=−1−(2Φ/C^2)とする。  Φ;重力ポテンシャル −−−−−??

断っておくがこの章は教科書からかなり逸脱している。教科書通りだとまず分からない。さてΦは重力ポテンシャルである。引力とは(gradΦ)である。よって力をF=mMG/R^2とした時、Φ=−MG/Rとなる。つまり 座標の原点に質量Mを置いた時、質量の外側でかつ原点からRだけ離れたところのΦはRに反比例する。

g00=−1+(2MG/C^2)/R  テンソル表現ではg00=−1+(2MG/C^2)/X1である。

通常考えられる(2MG/C^2)/X1は非常に小さい値である。この値が次第に大きくなりついに1となると、g00=0となるのだが......その結果はどうなるかはさておき... (2MG/C^2)/X1=1となるX1の値をRに戻して

R=2MG/C^2  とする時、このRをシュバルツシルト半径 と言う。

このシュバルツシルト半径とは何か....それはまたまた さておき......(2MG/C^2)=aとしておく。aは普通 非常に小さい値である。(「この章の終りに」参照)

よって

g00=−(1−a/X1)   g00=1/g00=−1/(1−a/X1

これに対して

g11=1/(1−a/X1)   g11=1/g11=(1−a/X1

わけだが 今後 面倒な計算をまだ少しするので

g11=−(1−a(X1)^-1)^-1   g11=1/g11=−(1−a(X1)^-1)

とする。

早い話

g00=−(1−a/X1) g11=1/(1−a/X1  g22=(X1^2  g33=(X1sin(X2)^2  他は0

g00=−1/(1−a/X1) g11=(1−a/X1  g22=1/(X1^2  g33=1/(X1sin(X2)^2 

なら Rμν−(1/2)gμν=0 を満足する(ハズ)....教科書では!!!

そして再度クリストッフェルのカッコΓσλν={λσν}=(∂λμν+∂νλμ−∂μνλσμ/2 を計算する。....のだが、.......

       b)クリストッフェルのカッコ

何回もニャンカイも、指摘したと思うが、σμσμ=0(σμ)であるから

Γσλν(∂λσν+∂νλσ−∂σνλσσ/2である。

例によって青文字は0として

σ=3の場合Γ3λν(∂λ3ν+∂νλ33νλ33/2=(∂λ3ν+∂νλ333/2

これはλ,νのどちらかは3で残りは3以外である。λ=3とすると

Γ3(33ν+∂ν3333/2=(∂ν3333/2

Γ330(03333/2=0   Γ333(33333/2=0

Γ331(∂13333/2=(∂1(X1sin(X2)^2)/{2(X1sin(X2)^2}

2(X1sin(X2)/{2(X1sin(X2)^2}=1/(X1sin(X2)

結論  Γ331Γ3131/(X1sin(X2)

Γ332(∂23333/2=(∂2(X1sin(X2)^2)/{2(X1sin(X2)^2}

2(X1sin(X2))cos(X2)/{2(X1sin(X2)^2}=cos(X2)/(X1sin(X2)

結論  Γ332Γ323cos(X2)/(X1sin(X2)

σ=2の場合Γ2λν(∂λ2ν+∂νλ2−∂2νλ22/2

2νλ≠0はλ=ν=3 のみであるから先に計算する。

Γ233(323+∂332−∂23322/2

=−(∂23322/2=−(∂2(X1sin(X2)^2)/{2(X1^2}

2(X1sin(X2))cos(X2)/{2(X1^2}=sin(X2)cos(X2)/X1

λνなら2νλ=0であるからΓ2λν(∂λ2ν+∂νλ222/2 λ,νのどちらかは2で残りは2以外である。となるが、残りが0か3では結局0になるから残りは1しかない。

Γ221(221+∂12222/2=(∂12222/2

(∂1(X1^2)/{2(X1^2}=1/X1

結論  Γ221Γ212=1/X1

σ=1の場合Γ1λν(∂λ1ν+∂νλ1−∂1νλ11/2

1νλλ=νならばλ=0,1,2,3すべて0以外の解がある。参考 g11=(1−a/X1

Γ100(010+∂001−∂10011/2=(−∂100(1−a/X1/2

Γ111(∂111111−∂11111/2=(∂111(1−a/X1/2

Γ122(212+∂221−∂12211/2=(−∂122(1−a/X1/2

Γ133(313+∂331−∂13311/2=(−∂133(1−a/X1/2

計算する前に(忘れる前に)λνの場合も考慮しておく。これは1νλ=0だからΓ1λν(∂λ1ν+∂νλ111/2であり、λ,νのどちらかは1である。ただしそうすれば0でないgは11となるが、11X1しか含んでないので結局0になる。

では順番に行こう。下記の式をみて

g00=−(1−a/X1) g11=1/(1−a/X1  g22=(X1^2  g33=(X1sin(X2)^2 

1g00=−1(1−a/X1)=−a/X1^2

 1g111(1−a/X1^-1 ={−1/(1−a/X1^2}{/X1^2}

  1g22=2X1

  1g332sin(X2)^2X1

 よって(結論

Γ100(−∂100(1−a/X1/2=(1/2)(a/X1^2(1−a/X1

Γ111(∂111(1−a/X1/2=(−1/2)(a/X1^2)/(1−a/X1

Γ122(−∂122(1−a/X1/2=X1(1−a/X1

Γ133(−∂133(1−a/X1/2=−sin(X2)^2X1(1−a/X1

σ=0の場合Γ0λν(∂λ0ν+∂νλ00νλ00/2=(∂λ0ν+∂νλ000/2

λ,νのどちらかは0で残りは0以外である。となるが、残りが2か3では結局0になるから残りは1しかない。

Γ001Γ010(∂10001000/2

(∂10000/2=(∂100/2)/(1−a/X1

(1/2)(a/X1^2)/(1−a/X1

結論  Γ001Γ010(1/2)(a/X1^2)/(1−a/X1

−−−−−−−−−−−−−−−−−−−

これで全部だ。結論をまとめると

Γ001Γ010(1/2)(a/X1^2)/(1−a/X1{001}={100}

Γ100(1/2)(a/X1^2(1−a/X1{010}

Γ111(−1/2)(a/X1^2)/(1−a/X1{111}

Γ122X1(1−a/X1{212}

Γ133−sin(X2)^2X1(1−a/X1{313}

Γ221Γ212=1/X1{221}={122}

Γ331Γ3131/(X1sin(X2){331}={133}

Γ332Γ323cos(X2)/(X1sin(X2){332}={233}

なんとも スゴイ!

−−−−−−−−−−−−−−−−−−−−−−−−

    c)リ−マンクリストッフェルのテンソル R10,01

リ−マンクリストッフェルのテンソルRαβ,μυは(第5章 5−6)参照 ゲェッ〜その2

Rαβ,μυ=∂μβαυ−∂υβαμβλυ}{λαμβλμ}{λαυとなる。5-6-1)

υ=α、μ=βとして

Rαβ,βα=∂ββαα−∂αβαββλα}{λαββλβ}{λαα

β=0として

Rα0,0α00αα−∂α0α00λα}{0αλ0λ0}{λαα

=−∂α0α00λα}{0αλ0λ0}{λαα

α=1の場合

R10,01=−∂10100λ1}{01λ0λ0}{λ11

右辺をふたつに分けて計算する。まず−∂1010

−∂1010}=−∂1(1/2)(a/X1^2(1−a/X1

(−a/2)∂1{(1^-2−a(X1^-3)}

(−a/2){-2(1^-3+3a(X1^-4)}

a(1^-3(3/2)a^2(X1^-4)}

 

次に0λ1}{01λ0λ0}{λ11}   ;青字は0

001}{010000}{011011}{011010}{111

021}{012020}{211031}{013030}{311

={001}{010010}{111010}({001111})

(1/2)(a/X1^2(1−a/X1){(1/2)(a/X1^2)/(1−a/X1)−(−1/2)(a/X1^2)/(1−a/X1

(1/2)(a/X1^2(1−a/X1(a/X1^2)/(1−a/X1

(1/2)(a/X1^2(1−a/X1(a/X1^2)/(1−a/X1

(1/2)(a/X1^2^2

よって

R10,01=−∂1010}+{010}({001111})

a(1^-3(3/2)^2(1^-4+(1/2)a^2(X1^4)

=a(1^-3−a^2(X1^4)

=(a/1^3)(1−/X1)

=−(a/1^3)g00   ;g00=−(1−a/X1

つまりR00R10,01=−(a/1^3)g00である。

     d)リ−マンクリストッフェルのテンソル R01,10

R01,10も同じように計算できる。

β=1として

Rα1,1α=∂11αα−∂α1α11λα}{λα11λ1}{λαα

α=0として

R01,10=∂110001011λ0}{λ011λ1}{λ00

=∂11001λ0}{λ011λ1}{λ00

ところで{100101}=(1/2)(a/X1^2)/(1−a/X1)であるので、

1100=∂1101なのだがこれはR10,01を出す時計算している。

1100=∂1101=−a(1^-3)+(3/2)a^2(X1^-4

残りは

1λ0}{λ011λ1}{λ00

={100}{001101}{000110}{101111}{100

120}{201121}{200130}{301131}{300

100}({001111})

=(1/2)(a/X1^2)/(1−a/X1(a/X1^2)/(1−a/X1

=(1/2)(a^2)(1^-4)/(1−a/X1^2

結局

R01,10=−a(1^-3)+(3/2)a^2(X1^-4)+(1/2)(a^2)(1^-4)/(1−a/X1^2

となるのだが....少しややこしいので a/X1=Qとする。

R01,10

=−Q(1^-2)+(3/2)QQ(X1^-2)+(1/2)QQ(1^-2)/(1-Q)(1-Q)

=−Q(1^-2){1-(3/2)Q+(1/2)Q/((1-Q)(1-Q))}

=−Q(1^-2(1/2){2-3Q+Q/((1-Q)(1-Q))}

=−Q(1^-2(1/2)/((1-Q)(1-Q)){2(1-Q)(1-Q)-3Q(1-Q)(1-Q)+Q}

{ }内

2(1-Q)(1-Q)-3Q(1-Q)(1-Q)+Q

2-4Q+QQ-3Q+3QQ+Q

2-6Q+4QQ=2(1-Q)(1-2Q)

改めて書きなおすと

R01,10=−Q(1^-2(1/2)/((1-Q)(1-Q)){2(1-Q)(1-Q)-3Q(1-Q)(1-Q)+Q}

=−Q(1^-2(1/2)2(1-Q)(1-2Q)/((1-Q)(1-Q))

=−Q(1^-2(1-2Q)/(1-Q)

=−a/X1(1^-2(1-2a/X1)/(1-a/X1)

=−(1^-3/(1-a/X1)−(2^2/X1^4)/(1-a/X1)

ここでaは非常に小さいからa^2は0とみなし、(2^2/X1^4)/(1-a/X1)=0とすると

R01,10=−(1^-3/(1-a/X1)となる。

つまりR11R01,10=−(1^-3/(1-a/X1)=−(a/1^3g11である。; g11=1/(1−a/X1

(^-^)(^-^)(^-^)

      e)どう曲がってる?

「c)リ−マンクリストッフェルのテンソル R10,01」の結論 R00R10,01=−(a/1^3)g00 と

「d)リ−マンクリストッフェルのテンソル R01,10」の結論 R11R01,10=−(a/1^3g11 は

Rμν−(1/2)gμν=0  ;アインシュタインのテンソルの式

を満足しそうである。少なくとも「6-3-2)2次元空間の計算」の「b)アインシュタインのテンソルの式を適用する。」ではテンソルRμνとgμνを行列と見立てたならばその各々の成分 R00とg00 , R11νとg11に対して

R00−(1/2)g00=0 R11−(1/2)g11=0

をも満足する.....まで言いきってないかも知れないけど....そうなのだ。

そして c),d)における−(a/1^3)は(1/2)に該当する。つまりスカラ−曲率は−2a/(1^3)なのである。でもそれは結局 ナニがどうなることか?

 まず曲がってない空間では計量テンソルηは、直行座標、曲座標によって多少違うが、成分(00)、(11)については-1、1となる言わば不動の値を取る。だが曲がった空間では計量テンソルgの(00)成分、(11)成分が変化する。

今回のように重力の元での曲座標による計量テンソルの変化はg00は小さくなりg11成分は大きくなる。(上 マンガ)

空間の曲がりが非常に大きくなりついにg00=0となる場所の原点からの半径をシュバルツシルト半径と言う。その場所ではg11=∞となる。

00=0となるとは計量テンソルに寄与するCT つまり時間がなくなることを意味する。すなわち次元がひとつ減る!。すなわち時間が経過しない!。すなわち時間が凍りつく!。ブラックホ−ルに落ちて行く時、シュバルツシルト半径に到達するのに無限の時間がかかるというのは最後には時間がなくなることを意味する。時間が無いから因果律自体が成り立たない。それでだろうと思われるが、シュバルツシルト半径による地平線を「因果の地平線」と言う因果な名前で呼ぶ。

 

6−5)この章の終りに...実は口上

 普通 口上と言うものはなにかを始める前に言うものだけど、今回は逆!。いきなりシュバルツシルト半径とか因果の地平線なんかが唐突に出て来るから驚かれたかもしれない。普通、あまり科学的とは言い難い科学本でなくとも、「アインシュタインの重力場の方程式を一切の近似ナシに完全に解いたのはシュバルツシルトが初めてである」くらいは述べる。

 実際問題として重力場の方程式をちゃんと解くのはメチャクチャ難しい。唯一 曲座標でなら簡単な設定ではいたる所に計算上0が多発するのでなんとか解けるわけである。

 ただし 今回は教科書を大幅に逸脱している。この章は完全にオリジナルである。それでも結構あれこれ計算できるし、あれこれ述べることが可能らしい。

 まぁ 以上のことを口上として述べてから始めるべきだった。スミマセン!!

ところで計算に使った a=(2MG/C^2) である。別名シュバルツシルト半径である。 今現在、我々は重力g(=9.8m/S^2)のもとにいるのだがaの値に直すとどのくらいか?

F=mg  F=MG/r  M;地球の質量 G;重力定数  r;地球の半径

となるから  MG=gとなり

a=2g(r/C)^2 である。

ここで「光は1秒間に地球を7.5周する」と言う格言(?)を使う。

式に直すと2πr×7.5=C →r/C=1/(2π×7.5)=0.02

a=2g(r/C)^2=0.0078 m   7.8 mm 結構大きいんだよねぇ〜??....

 ところで5次元空間について少し話す!。4次元空間は空間座標X,Y,ZにCTが加わり4個の次元を作る。それに対して5次元では...まぁ 空間座標X,Y,Zに空間座標Wを加え、CTも合わせて5次元とするなら

S^2=−(CT)^2+X^2+Y^2+Z^2+^2

として後を考えれば良い。良いのだが、これを曲座標でどう表せるのか?

X=Rsin(θ)cos(ψ)   Y=Rsin(θ)sin(ψ)   Z=Rcos(θ)  W=??

 このWの曲座標形式を固定しない限り、5次元空間でアインシュタインの方程式が満足されるか否か、判断できない。判断できないことは真偽も判断できない。 真偽は別として5次元空間、6次元空間理論を主張するものはあるのだが、そこまで書き込まれた理論は まだ知らない。

 さてHPの行方だが、最低限 ブラックホ−ルの内部解までは探ってみる。タイトルはシュバルツシルトを全面に出して「第7章 シュバルツシルトの解(怪?)」としておく。

 多分アップは来年になると思う。 それでは良いお年を....って ちょっと早すぎるか!.....

チョンチョン

トップに戻る