tecktom.com

出てくるもの1:第零法則と第一法則

2026年10月15日8分で読めますAI生成教科書『虎の巻』第2部 第19章 ブラックホール力学の4法則
目次
  1. 第零法則
  2. 質量を面積と角運動量で書く
  3. 第一法則
  4. 粒子を落として確かめる
  5. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第2部 第19章 §2.19.3 ★★★ 表面重力は地平面の上で一定で(第零法則)、質量の公式 M(A,J)M(A,J) の偏微分が κc2/(8πG)\kappa c^2/(8\pi G) と ΩH\Omega_{\mathrm H} を与える(第一法則)

第零法則

式 (2.19.11) は r±r_\pm と aa だけで書け、地平面の上の位置 θ\theta によらない。Kerr の地平面は赤道の方が膨らんだ回転楕円体のような形をしているが、表面重力は極でも赤道でも同じである。Newton の重力で回転する星の表面重力が極と赤道で異なるのと対照的である。

これは Kerr 解に限った偶然ではない。Bardeen・Carter・Hawking(1973、Communications in Mathematical Physics 31、161)は、定常なブラックホールの地平面で κ\kappa が一定であることを、Einstein 方程式と物質の優勢エネルギー条件(どの観測者が測ってもエネルギー密度が非負で、エネルギーの流れが光速を超えない)を仮定して示した。熱力学の第零法則「平衡にある系では温度がいたるところ等しい」と同じ形である。定常であることが熱平衡に、κ\kappa が温度に対応する。

質量を面積と角運動量で書く

MM と JJ を与えると地平面の面積 AA が決まる(式 (2.19.4))。逆に AA と JJ から MM を求めてみる。面積が同じになる Schwarzschild ブラックホールの質量を MirrM_{\mathrm{irr}} と書く:

A=16πG2Mirr2c4.(2.19.15)A=16\pi\frac{G^2M_{\mathrm{irr}}^2}{c^4}.\tag{2.19.15}

式 (2.19.4) の A=8πrgr+A=8\pi r_{\mathrm g}r_+ を使うと Mirr2=M2r+/(2rg)M_{\mathrm{irr}}^2=M^2r_+/(2r_{\mathrm g}) で、r+r_+ を消去して整理すると

M2=Mirr2+J2c24G2Mirr2(2.19.16)M^2=M_{\mathrm{irr}}^2+\frac{J^2c^2}{4G^2M_{\mathrm{irr}}^2}\tag{2.19.16}

を得る。Christodoulou(1970、Physical Review Letters 25、1596)が導いた式で、MirrM_{\mathrm{irr}} を既約質量と呼ぶ(呼び名の理由は次節で分かる)。導出を追うと、r+=rg+rg2−a2r_+=r_{\mathrm g}+\sqrt{r_{\mathrm g}^2-a^2} より (r+−rg)2=rg2−a2(r_+-r_{\mathrm g})^2=r_{\mathrm g}^2-a^2、すなわち r+2−2rgr++a2=0r_+^2-2r_{\mathrm g}r_++a^2=0 で、これに r+=2rgMirr2/M2r_+=2r_{\mathrm g}M_{\mathrm{irr}}^2/M^2 と a=J/(Mc)a=J/(Mc)、rg=GM/c2r_{\mathrm g}=GM/c^2 を代入して 4G2M2/c44G^2M^2/c^4 で割れば式 (2.19.16) になる。右辺は AA と JJ だけの関数なので、

Mc2=c2c4A16πG2+4πJ2c2A(2.19.17)Mc^2=c^2\sqrt{\frac{c^4A}{16\pi G^2}+\frac{4\pi J^2}{c^2A}}\tag{2.19.17}

と、質量エネルギーが (A,J)(A,J) の関数として書けた。

第一法則

式 (2.19.17) の全微分を取ると

c2δM=κc28πG δA+ΩH δJ(2.19.18)c^2\delta M=\frac{\kappa c^2}{8\pi G}\,\delta A+\Omega_{\mathrm H}\,\delta J\tag{2.19.18}

となる。係数がちょうど表面重力 (2.19.11) と地平面の角速度 (2.19.6) になることは、偏微分を計算して r±r_\pm で書き直すと確かめられる。Schwarzschild(J=0J=0)で確かめておくと、A=16πG2M2/c4A=16\pi G^2M^2/c^4 から δA=32πG2MδM/c4\delta A=32\pi G^2M\delta M/c^4 で、

κc28πGδA=c44GM⋅c28πG⋅32πG2Mc4δM=c2δM\frac{\kappa c^2}{8\pi G}\delta A=\frac{c^4}{4GM}\cdot\frac{c^2}{8\pi G}\cdot\frac{32\pi G^2M}{c^4}\delta M=c^2\delta M

となって合う。一般の aa では verify_ch2-19.py で偏微分 ∂(Mc2)/∂A∣J=κc2/(8πG)\partial(Mc^2)/\partial A|_J=\kappa c^2/(8\pi G)、∂(Mc2)/∂J∣A=ΩH\partial(Mc^2)/\partial J|_A=\Omega_{\mathrm H} を確かめた。

式 (2.19.18) は熱力学の第一法則 dU=TdS+ΩdJ\mathrm{d} U=T\mathrm{d} S+\Omega\mathrm{d} J(回転する物体の場合)と同じ形をしている。第3部第4章では内部エネルギー U(S,V,N)U(S,V,N) の偏微分が温度・圧力・化学ポテンシャルを与えることを見るが、ここでは Mc2(A,J)Mc^2(A,J) の偏微分が κc2/(8πG)\kappa c^2/(8\pi G) と ΩH\Omega_{\mathrm H} を与えている。AA がエントロピーに、κ\kappa が温度に対応する。第1部第7章の言葉では、(κ,ΩH)(\kappa,\Omega_{\mathrm H}) は (A,J)(A,J) と Legendre 変換で入れ替えられる共役な変数である。

粒子を落として確かめる

第一法則を、2.19.2節の保存量で読み直す。エネルギー EE、角運動量 LzL_z の試験粒子を落とすと、式 (2.19.8) と (2.19.18) から

κc28πG δA=E−ΩHLz(2.19.19)\frac{\kappa c^2}{8\pi G}\,\delta A=E-\Omega_{\mathrm H}L_z\tag{2.19.19}

となる。右辺の符号は、粒子が地平面を越えられるという条件で決まる。粒子の4元運動量 pμp^\mu は未来向きの時間的(光子なら光的)で、地平面の生成子 χ\chi も未来向きの光的なベクトルである。(+,−,−,−)(+,-,-,-) の流儀では、未来向きの時間的または光的なベクトル2つの内積は非負なので、地平面の上で

χμpμ=p0+ΩHcpφ=1c(E−ΩHLz)≥0(2.19.20)\chi^\mu p_\mu=p_0+\frac{\Omega_{\mathrm H}}{c}p_\varphi=\frac{1}{c}\bigl(E-\Omega_{\mathrm H}L_z\bigr)\ge0\tag{2.19.20}

が成り立つ。EE と LzL_z は保存量なので、地平面で評価した値は遠方で準備したときの値と同じである。したがって δA≥0\delta A\ge0 で、粒子を落とすどんな過程でも地平面の面積は減らない。等号 E=ΩHLzE=\Omega_{\mathrm H}L_z は、粒子が地平面の生成子とちょうど同じ向きに回りながら接して入る極限で、このとき面積は1次の範囲で変わらない。これを可逆な変化と呼ぶ。熱力学で準静的な過程がエントロピーを変えないのと同じ位置づけである。

例題 2.19.3 スピン 0.8 のブラックホールに粒子を落とす ★★

オリジナル例題 2.19.1 のブラックホール(M=10M⊙M=10M_\odot、a∗=0.8a_\ast=0.8)に、エネルギー E=10−4Mc2E=10^{-4}Mc^2 の粒子を落とす。

  1. (1) Lz=0L_z=0 のとき、第一法則から地平面の面積の相対変化 δA/A\delta A/A を求めよ。
  2. (2) 面積を1次で変えない角運動量 LzL_z を求め、JJ との比を答えよ。
解答

(1) 式 (2.19.11) で r+=1.6 rgr_+=1.6\,r_{\mathrm g}、r−=0.4 rgr_-=0.4\,r_{\mathrm g}、r+2+a2=3.2 rg2r_+^2+a^2=3.2\,r_{\mathrm g}^2 から κ=(3/16)c2/rg\kappa=(3/16)c^2/r_{\mathrm g}(Schwarzschild の値 c2/(4rg)c^2/(4r_{\mathrm g}) の 0.75 倍)。式 (2.19.19) から

δA=8πGκc2E=8πG⋅10−4M(3/16)c2/rg=8π⋅10−4⋅163rg2,\delta A=\frac{8\pi G}{\kappa c^2}E=\frac{8\pi G\cdot10^{-4}M}{(3/16)c^2/r_{\mathrm g}}=\frac{8\pi\cdot10^{-4}\cdot16}{3}r_{\mathrm g}^2,

A=12.8πrg2A=12.8\pi r_{\mathrm g}^2 で割って δA/A=3.33×10−4\delta A/A=3.33\times10^{-4}。MM と JJ を変えた Kerr 解の面積を厳密に計算し直しても 3.33×10−43.33\times10^{-4} で、違いは 10−810^{-8} の桁の2次の項だけである。

(2) E=ΩHLzE=\Omega_{\mathrm H}L_z とすればよい。ΩH=c/(4rg)\Omega_{\mathrm H}=c/(4r_{\mathrm g})、J=Mac=0.8MrgcJ=Mac=0.8Mr_{\mathrm g}c から

LzJ=EΩHJ=10−4Mc2(c/4rg)(0.8Mrgc)=5×10−4.\frac{L_z}{J}=\frac{E}{\Omega_{\mathrm H}J}=\frac{10^{-4}Mc^2}{(c/4r_{\mathrm g})(0.8Mr_{\mathrm g}c)}=5\times10^{-4}.

厳密に計算し直した面積の相対変化は −6.7×10−8-6.7\times10^{-8} で、(1) の値より4桁小さい2次の量になる。回転の向きにたっぷり角運動量を持たせて落とすと、エネルギーを足しても面積はほとんど増えない。

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『虎の巻』の記事一覧へ