tecktom.com

出てくるもの1:測地線方程式と Christoffel 記号

2026年10月15日11分で読めますAI生成教科書『虎の巻』第2部 第12章 測地線は固有時の作用の停留
目次
  1. Euler–Lagrange 方程式を計算する
  2. 径数を固有時に選ばなかったとき
  3. 平方根を外した
  4. 平坦な時空でも は 0 でない
  5. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第2部 第12章 §2.12.3 ★★★ Euler–Lagrange 方程式から測地線方程式と Christoffel 記号が出る。固有時を径数にすると平方根を外した L=12gx˙x˙L=\frac12g\dot x\dot x に置き換えられる

Euler–Lagrange 方程式を計算する

w=gαβ(x)x˙αx˙βw=g_{\alpha\beta}(x)\dot x^\alpha\dot x^\beta とおくと L=−mcwL=-mc\sqrt w で、

∂L∂x˙μ=−mcw gμνx˙ν,∂L∂xμ=−mc2w ∂μgαβ x˙αx˙β\frac{\partial L}{\partial \dot x^\mu}=-\frac{mc}{\sqrt w}\,g_{\mu\nu}\dot x^\nu,\qquad \frac{\partial L}{\partial x^\mu}=-\frac{mc}{2\sqrt w}\,\partial_\mu g_{\alpha\beta}\,\dot x^\alpha\dot x^\beta

である(∂μ=∂/∂xμ\partial_\mu=\partial/\partial x^\mu)。Euler–Lagrange 方程式 (2.12.7) を −mc-mc で割ると

ddλ(gμνx˙νw)−12w ∂μgαβ x˙αx˙β=0.(2.12.9)\frac{\mathrm{d} }{\mathrm{d} \lambda}\Bigl(\frac{g_{\mu\nu}\dot x^\nu}{\sqrt w}\Bigr)-\frac{1}{2\sqrt w}\,\partial_\mu g_{\alpha\beta}\,\dot x^\alpha\dot x^\beta=0 .\tag{2.12.9}

ここで 2.12.2節の順序どおり、方程式を書いたあとで λ=τ\lambda=\tau と置く。式 (2.12.3) から w=c\sqrt w=c は定数になり、x˙μ=uμ\dot x^\mu=u^\mu として

ddτ(gμνuν)−12∂μgαβ uαuβ=0.\frac{\mathrm{d} }{\mathrm{d} \tau}\bigl(g_{\mu\nu}u^\nu\bigr)-\frac12\partial_\mu g_{\alpha\beta}\,u^\alpha u^\beta=0 .

第1項は gμνg_{\mu\nu} が世界線に沿って変わる分を含み、dgμνdτ=∂αgμνuα\frac{\mathrm{d} g_{\mu\nu}}{\mathrm{d} \tau}=\partial_\alpha g_{\mu\nu}u^\alpha なので

gμνduνdτ+∂αgμβ uαuβ−12∂μgαβ uαuβ=0.g_{\mu\nu}\frac{\mathrm{d} u^\nu}{\mathrm{d} \tau}+\partial_\alpha g_{\mu\beta}\,u^\alpha u^\beta-\frac12\partial_\mu g_{\alpha\beta}\,u^\alpha u^\beta=0 .

∂αgμβuαuβ\partial_\alpha g_{\mu\beta}u^\alpha u^\beta は α\alpha と β\beta を入れ替えても同じなので 12(∂αgμβ+∂βgμα)uαuβ\frac12(\partial_\alpha g_{\mu\beta}+\partial_\beta g_{\mu\alpha})u^\alpha u^\beta と書け、まとめると

gμνduνdτ+Γμαβ uαuβ=0,Γμαβ=12(∂αgμβ+∂βgμα−∂μgαβ)(2.12.10)g_{\mu\nu}\frac{\mathrm{d} u^\nu}{\mathrm{d} \tau}+\Gamma_{\mu\alpha\beta}\,u^\alpha u^\beta=0,\qquad \Gamma_{\mu\alpha\beta}=\frac12\bigl(\partial_\alpha g_{\mu\beta}+\partial_\beta g_{\mu\alpha}-\partial_\mu g_{\alpha\beta}\bigr)\tag{2.12.10}

となる。両辺に gρμg^{\rho\mu} を掛ける。

定理 2.12.2(測地線方程式 ★★★)

作用 (2.12.4) の停留点を固有時 τ\tau で表した世界線は

d2xρdτ2+Γρμνdxμdτdxνdτ=0,Γρμν=12gρσ(∂μgσν+∂νgσμ−∂σgμν)(2.12.11)\frac{\mathrm{d}^2x^\rho}{\mathrm{d}\tau^2}+\Gamma^\rho{}_{\mu\nu}\frac{\mathrm{d} x^\mu}{\mathrm{d} \tau}\frac{\mathrm{d} x^\nu}{\mathrm{d} \tau}=0,\qquad \Gamma^\rho{}_{\mu\nu}=\frac12g^{\rho\sigma}\bigl(\partial_\mu g_{\sigma\nu}+\partial_\nu g_{\sigma\mu}-\partial_\sigma g_{\mu\nu}\bigr)\tag{2.12.11}

を満たす。Γρμν\Gamma^\rho{}_{\mu\nu} を Christoffel 記号と呼ぶ。下の2つの添字について対称(Γρμν=Γρνμ\Gamma^\rho{}_{\mu\nu}=\Gamma^\rho{}_{\nu\mu})である。

測地線方程式の第2項は、Newton の運動方程式で「力 /m/m」があった場所にある。ただし Γ\Gamma は質量 mm を含まないので、どの粒子も同じ世界線を進む。これは第2部第11章の等価原理(重力による加速度は質量によらない)を式にしたものである。mm は作用 (2.12.4) の全体に掛かる定数にすぎず、停留の条件から消える。

径数を固有時に選ばなかったとき

式 (2.12.9) を λ=τ\lambda=\tau と置かずに計算すると、s(λ)=w=c dτdλs(\lambda)=\sqrt w=c\,\frac{\mathrm{d} \tau}{\mathrm{d} \lambda} として

x¨ρ+Γρμνx˙μx˙ν=s˙s x˙ρ(2.12.12)\ddot x^\rho+\Gamma^\rho{}_{\mu\nu}\dot x^\mu\dot x^\nu=\frac{\dot s}{s}\,\dot x^\rho\tag{2.12.12}

になる(ドットは λ\lambda 微分)。右辺は世界線に沿った向き x˙ρ\dot x^\rho への加速度で、径数の進み方の不揃いを表しているだけである。同じ世界線を径数だけ取り替えた解なので、世界線の形は変わらない。右辺が 00 になるのは s˙=0\dot s=0、つまり λ=aτ+b\lambda=a\tau+b(a,ba,b は定数)のときに限る。このように測地線方程式を右辺 00 の形にする径数をアフィン径数と呼ぶ。

平方根を外した L2=12gμνx˙μx˙νL_2=\frac12g_{\mu\nu}\dot x^\mu\dot x^\nu

アフィン径数だけを使うなら、平方根の付いた LL の代わりに

L2=12gμν(x) x˙μx˙ν(2.12.13)L_2=\frac12g_{\mu\nu}(x)\,\dot x^\mu\dot x^\nu\tag{2.12.13}

を使ってよい。L2L_2 の Euler–Lagrange 方程式は、∂L2∂x˙μ=gμνx˙ν\frac{\partial L_2}{\partial \dot x^\mu}=g_{\mu\nu}\dot x^\nu と ∂L2∂xμ=12∂μgαβx˙αx˙β\frac{\partial L_2}{\partial x^\mu}=\frac12\partial_\mu g_{\alpha\beta}\dot x^\alpha\dot x^\beta から、上と同じ計算で

gμν(x¨ν+Γναβx˙αx˙β)=0g_{\mu\nu}\bigl(\ddot x^\nu+\Gamma^\nu{}_{\alpha\beta}\dot x^\alpha\dot x^\beta\bigr)=0

となり、径数を固有時に選んだ式 (2.12.11) がそのまま出る。

右辺に s˙/s\dot s/s の項が出ないのは、L2L_2 が径数の取り替えで変わらないという性質を持たないからである。L2L_2 は λ\lambda をあらわに含まないので、第1部第10章の Noether の定理(時刻の並進)から H2=x˙μ∂L2∂x˙μ−L2H_2=\dot x^\mu\frac{\partial L_2}{\partial \dot x^\mu}-L_2 が保存し、L2L_2 は x˙\dot x の2次の同次関数なので H2=L2H_2=L_2 である(第1部第10章)。つまり L2L_2 の解では w=2L2w=2L_2 が自動的に一定で、径数は自動的にアフィン径数になる。逆に、LL の解のうちアフィン径数で書いたものは L2L_2 の方程式を満たす。2つの LL は同じ世界線を与え、違うのは径数の選び方の自由度だけである。

L2L_2 には LL に無い利点が2つある。1つ目は、Hessian が gμνg_{\mu\nu} そのもので逆行列を持つことである。第1部第7章の Legendre 変換がそのまま使え、pμ=gμνx˙νp_\mu=g_{\mu\nu}\dot x^\nu から

H2=12gμν(x) pμpν(2.12.14)H_2=\frac12g^{\mu\nu}(x)\,p_\mu p_\nu\tag{2.12.14}

を得る(第1部第7章)。λ=τ\lambda=\tau のとき H2=c2/2H_2=c^2/2 で、これは第2部第3章の質量殻の拘束 pμpμ=m2c2p_\mu p^\mu=m^2c^2 を m2m^2 で割ったものにあたる。平方根の付いた LL では HH が恒等的に 00 で拘束が別に残ったが、L2L_2 ではその拘束が保存する H2H_2 の値になっている。2つ目は、m=0m=0 の粒子(光)を扱えることである。光の世界線は ds=0\mathrm{d} s=0 なので LL は使えないが、L2L_2 では w=0w=0 の解として

x¨ρ+Γρμνx˙μx˙ν=0,gμνx˙μx˙ν=0(2.12.15)\ddot x^\rho+\Gamma^\rho{}_{\mu\nu}\dot x^\mu\dot x^\nu=0,\qquad g_{\mu\nu}\dot x^\mu\dot x^\nu=0\tag{2.12.15}

を満たす世界線(光的な測地線)を考えることができる。光の場合は固有時が無いので、径数はアフィン径数としか言えない。

平坦な時空でも Γ\Gamma は 0 でない

Γρμν\Gamma^\rho{}_{\mu\nu} は計量の1階微分で書けているので、平坦な時空でも曲がった座標を使えば 00 でない。次の例題がその例である。重力の無い時空で Γ≠0\Gamma\neq0 になり、重力のある時空でも1点では Γ=0\Gamma=0 にできる(自由落下する系を選ぶ。第2部第11章)ので、Γ\Gamma そのものは「重力があるかどうか」を判定する量ではない。判定する量は計量の2階微分を含む曲率で、第2部第14章で扱う。Γ\Gamma がある座標系で 00 で別の座標系で 00 でないことは、Γ\Gamma が座標変換でテンソルのように変換しないことを意味し、そのことは第2部第13章で式で確かめる。

例題 2.12.3 平坦な時空の極座標の Christoffel 記号 ★★

オリジナル平坦な時空を極座標で ds2=c2dt2−dr2−r2dφ2\mathrm{d} s^2=c^2\mathrm{d} t^2-\mathrm{d} r^2-r^2\mathrm{d}\varphi^2 と書く(zz 方向の運動は考えない)。

  1. (1) 0 でない Christoffel 記号をすべて求めよ。
  2. (2) 時刻 tt を径数にして測地線方程式を書け(遅い粒子では tt は固有時に比例する)。
  3. (3) 直線 x=bx=b, y=V(t−t0)y=V(t-t_0) を極座標で書いたものが (2) を満たすことを確かめよ。
解答

(1) 0 でない計量の微分は ∂rgφφ=−2r\partial_rg_{\varphi\varphi}=-2r だけである。grr=−1g^{rr}=-1、gφφ=−1/r2g^{\varphi\varphi}=-1/r^2 なので

Γrφφ=12grr(−∂rgφφ)=−r,Γφrφ=Γφφr=12gφφ∂rgφφ=1r\Gamma^r{}_{\varphi\varphi}=\frac12g^{rr}(-\partial_rg_{\varphi\varphi})=-r,\qquad \Gamma^\varphi{}_{r\varphi}=\Gamma^\varphi{}_{\varphi r}=\frac12g^{\varphi\varphi}\partial_rg_{\varphi\varphi}=\frac1r

で、ほかはすべて 00 である。

(2) 計量は tt によらず g00=1g_{00}=1 なので Γ0μν=0\Gamma^0{}_{\mu\nu}=0 で、dtdτ\frac{\mathrm{d} t}{\mathrm{d} \tau} は一定になる。tt を径数にしてよく、

r¨−rφ˙2=0,φ¨+2rr˙φ˙=0.\ddot r-r\dot\varphi^2=0,\qquad \ddot\varphi+\frac{2}{r}\dot r\dot\varphi=0 .

第1式の −rφ˙2-r\dot\varphi^2 は、極座標で習う加速度の動径成分 r¨−rφ˙2\ddot r-r\dot\varphi^2 の第2項である。第2式に r2r^2 を掛けると ddt(r2φ˙)=0\frac{\mathrm{d} }{\mathrm{d} t}(r^2\dot\varphi)=0 で、面積速度の保存になる。

(3) r=b2+V2(t−t0)2r=\sqrt{b^2+V^2(t-t_0)^2}、tan⁡φ=V(t−t0)/b\tan\varphi=V(t-t_0)/b である。φ˙=bV/r2\dot\varphi=bV/r^2 なので r2φ˙=bVr^2\dot\varphi=bV は一定で第2式を満たす。r˙=V2(t−t0)/r\dot r=V^2(t-t_0)/r、r¨=V2/r−V4(t−t0)2/r3=V2b2/r3\ddot r=V^2/r-V^4(t-t_0)^2/r^3=V^2b^2/r^3 で、rφ˙2=b2V2/r3r\dot\varphi^2=b^2V^2/r^3 と等しいので第1式も満たす。重力の無い時空の測地線は直線であり、極座標で Γ≠0\Gamma\neq0 になったのは座標が曲がっているためである。

木下の勉強ステータス

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

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

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