tecktom.com

測地線偏差と潮汐力

2026年10月22日17分で読めますAI生成教科書『基礎理論』第11章 曲率テンソル
目次
  1. 測地線の族
  2. 2次元の Jacobi 方程式
  3. Newton の重力の潮汐は曲率である
  4. 点質量の潮汐
  5. 確認テスト
  6. 木下の層

AI 生成・木下が未勉強

物理のための微分幾何入門 第11章 §11.4 ★★★ 近くを通る2本の測地線の間隔の加速度が曲率で決まることを示し、Newton 力学の潮汐をその例として読む

測地線の族

近くを通る測地線の族を考える。2つの変数の写像 xμ(t,s)x^\mu(t,s) を考える。各 ss について、tt の関数として見たこの写像がアフィンパラメータの測地線になるものとする。その上の2つのベクトル場を次で定める。

Tμ=∂xμ∂t,Jμ=∂xμ∂s.(11.19)T^\mu=\frac{\partial x^\mu}{\partial t},\qquad J^\mu=\frac{\partial x^\mu}{\partial s}.\tag{11.19}

TT は測地線の速度で、JJ は隣の測地線への間隔を表す。JJ を Jacobi 場(偏差ベクトル)と呼ぶ。

曲線に沿う共変微分を D/dtD/\mathrm{d} t・D/dsD/\mathrm{d} s と書く。たとえば DZρ/ds=∂sZρ+ΓρνσJνZσDZ^\rho/\mathrm{d} s=\partial_sZ^\rho+\Gamma^\rho{}_{\nu\sigma}J^\nu Z^\sigma である。

補題 11.15

(1) DJρ/dt=DTρ/dsDJ^\rho/\mathrm{d} t=DT^\rho/\mathrm{d} s。(2) 族の上のベクトル ZZ について

DdtDdsZρ−DdsDdtZρ=RρσμνZσTμJν.(11.20)\frac{D}{\mathrm{d} t}\frac{D}{\mathrm{d} s}Z^\rho-\frac{D}{\mathrm{d} s}\frac{D}{\mathrm{d} t}Z^\rho=R^\rho{}_{\sigma\mu\nu}Z^\sigma T^\mu J^\nu .\tag{11.20}

証明

(1) ∂tJρ=∂t∂sxρ=∂sTρ\partial_tJ^\rho=\partial_t\partial_sx^\rho=\partial_sT^\rho である。接続の項は ΓρμνTμJν\Gamma^\rho{}_{\mu\nu}T^\mu J^\nu と ΓρνμJνTμ\Gamma^\rho{}_{\nu\mu}J^\nu T^\mu で、下の添字の対称性から等しい。

(2) DdsZρ=∂sZρ+ΓρνσJνZσ\dfrac{D}{\mathrm{d} s}Z^\rho=\partial_sZ^\rho+\Gamma^\rho{}_{\nu\sigma}J^\nu Z^\sigma に D/dtD/\mathrm{d} t を作用させる。接続の tt 微分は ∂tΓ=Tμ∂μΓ\partial_t\Gamma=T^\mu\partial_\mu\Gamma なので

DdtDdsZρ=∂t∂sZρ+Tμ∂μΓρνσJνZσ+Γρνσ(∂tJνZσ+Jν∂tZσ)+ΓρμλTμ(∂sZλ+ΓλνσJνZσ).\frac{D}{\mathrm{d} t}\frac{D}{\mathrm{d} s}Z^\rho=\partial_t\partial_sZ^\rho+T^\mu\partial_\mu\Gamma^\rho{}_{\nu\sigma}J^\nu Z^\sigma +\Gamma^\rho{}_{\nu\sigma}(\partial_tJ^\nu Z^\sigma+J^\nu\partial_tZ^\sigma) +\Gamma^\rho{}_{\mu\lambda}T^\mu\bigl(\partial_sZ^\lambda+\Gamma^\lambda{}_{\nu\sigma}J^\nu Z^\sigma\bigr).

この式から、2つのパラメータを入れ替えた式を引く(T↔JT\leftrightarrow J)。2階の偏微分 ∂t∂sZ\partial_t\partial_sZ は消える。(1) の ∂tJ=∂sT\partial_tJ=\partial_sT から、Γ ∂J\Gamma\,\partial J の項も消える。Γ ∂Z\Gamma\,\partial Z を含む4項も、2つずつ同じ形なので消える。

残るのは ∂Γ\partial\Gamma の2項と ΓΓ\Gamma\Gamma の2項である。

(∂μΓρνσ−∂νΓρμσ+ΓρμλΓλνσ−ΓρνλΓλμσ)ZσTμJν.\bigl(\partial_\mu\Gamma^\rho{}_{\nu\sigma}-\partial_\nu\Gamma^\rho{}_{\mu\sigma}+\Gamma^\rho{}_{\mu\lambda}\Gamma^\lambda{}_{\nu\sigma}-\Gamma^\rho{}_{\nu\lambda}\Gamma^\lambda{}_{\mu\sigma}\bigr)Z^\sigma T^\mu J^\nu .

括弧の中は (11.3) そのものである。

定理 11.16(測地線偏差の方程式 ★★★)

測地線の族の Jacobi 場は

D2Jρdt2=RρσμνTσTμJν(11.21)\frac{D^2J^\rho}{\mathrm{d} t^2}=R^\rho{}_{\sigma\mu\nu}T^\sigma T^\mu J^\nu\tag{11.21}

を満たす。これを Jacobi 方程式(測地線偏差の方程式)と呼ぶ。

証明

補題 11.15 の (1) から D2Jdt2=DdtDdsT\dfrac{D^2J}{\mathrm{d} t^2}=\dfrac{D}{\mathrm{d} t}\dfrac{D}{\mathrm{d} s}T である。(2) で Z=TZ=T とすると、これは DdsDdtT+RρσμνTσTμJν\dfrac{D}{\mathrm{d} s}\dfrac{D}{\mathrm{d} t}T+R^\rho{}_{\sigma\mu\nu}T^\sigma T^\mu J^\nu に等しい。各曲線が測地線なので DT/dt=0DT/\mathrm{d} t=0 で、第1項は消える。

平らな空間では右辺が 00 で、間隔は tt の1次式でしか変わらない。曲率があると、初めは平行に出発した測地線どうしが加速度を持って近づいたり離れたりする。逆に、測地線の間隔の加速度を測れば、曲率を知ることができる。

2次元の Jacobi 方程式

2次元で、速度 TT は長さ1、Jacobi 場 JJ はそれに直交するとする。例題 11.2 の式を (11.21) の右辺に入れると、g(T,J)=0g(T,J)=0、g(T,T)=1g(T,T)=1 から次の形になる。

RρσμνTσTμJν=K(Tρg(T,J)−Jρg(T,T))=−KJρ.R^\rho{}_{\sigma\mu\nu}T^\sigma T^\mu J^\nu=K\bigl(T^\rho g(T,J)-J^\rho g(T,T)\bigr)=-KJ^\rho .

測地線に沿って平行移動した長さ1の法線 NN を使い、J=j(t)NJ=j(t)N と書く。NN は平行なので D2J/dt2=j′′ND^2J/\mathrm{d} t^2=j''N である。

命題 11.17(2次元の Jacobi 方程式 ★★★)

2次元で、測地線に直交する Jacobi 場の大きさ j(t)j(t) は

j′′(t)+K j(t)=0(11.22)j''(t)+K\,j(t)=0\tag{11.22}

を満たす。KK は測地線上の点での Gauss 曲率である。

KK が一定で、1点から角度をつけて出発する族 j(0)=0, j′(0)=1j(0)=0,\ j'(0)=1 を考える。解は曲率の符号で3通りに分かれる。

  • K=1/a2>0K=1/a^2>0:j=asin⁡(t/a)j=a\sin(t/a)。間隔は最大になったあと縮み、t=πat=\pi a で再び 00 になる。
  • K=0K=0:j=tj=t。間隔は一定の割合で広がる。
  • K=−1/ℓ2<0K=-1/\ell^2<0:j=ℓsinh⁡(t/ℓ)j=\ell\sinh(t/\ell)。間隔は指数関数的に広がる。
図

図 11.3 1点から角度をつけて出発した2本の測地線の間隔 j(t)j(t)(j(0)=0j(0)=0、j′(0)=1j'(0)=1、∣K∣=1\lvert K\rvert=1)。正の曲率では t=πt=\pi で再び会い、負の曲率では指数関数的に離れる((11.22))。

正の曲率で再び会う点を共役点と呼ぶ。球面では北極から出たすべての大円が南極で会い、北極と南極は互いに共役点である。共役点を越えた測地線は、もう最短の道ではない(端点を少しずらした近くの道の方が短くなる)。

例題 11.4 球面の経線の間隔 ★★

オリジナル半径 aa の球面で、北極から角 δφ\delta\varphi をなして出発する2本の経線を考える。

  1. 北極からの弧長 tt での2本の間隔 j(t)j(t) を、(11.22) を解いて求めよ。
  2. 地球(a=6.371×106 ma=6.371\times10^6\ \mathrm m)で δφ=1∘\delta\varphi=1^\circ のとき、赤道での間隔を求めよ。
  3. 同じ結果を、緯線の半径から幾何で確かめよ。
解答

(1) 北極で2本は一致するので j(0)=0j(0)=0 である。出発の角が δφ\delta\varphi なので、北極のごく近くでは平面と同じく j≈t δφj\approx t\,\delta\varphi である。したがって j′(0)=δφj'(0)=\delta\varphi である。K=1/a2K=1/a^2 で (11.22) を解くと

j(t)=a δφ sin⁡(t/a)j(t)=a\,\delta\varphi\,\sin(t/a)

を得る。

(2) 赤道は t=πa/2t=\pi a/2 なので j=a δφj=a\,\delta\varphi である。δφ=π/180=0.017453\delta\varphi=\pi/180=0.017453 を入れると、次のようになる。

j=6.371×106×0.017453 m=1.112×105 m≈111 km.j=6.371\times10^6\times0.017453\ \mathrm m=1.112\times10^5\ \mathrm m\approx111\ \mathrm{km}.

(3) 北極からの弧長 tt の点は余緯度 t/at/a にある。そこの緯線の半径は asin⁡(t/a)a\sin(t/a) である。角 δφ\delta\varphi の経線どうしの間隔はその半径に δφ\delta\varphi を掛けたもので、(1) と一致する。t=πat=\pi a の南極で j=0j=0 となり、2本は共役点で再び会う。

Newton の重力の潮汐は曲率である

Newton 力学で、重力ポテンシャル Φ(x)\Phi(\bm x) の中を自由に運動する粒子は x¨=−∇Φ\ddot{\bm x}=-\nabla\Phi に従う。近くを運動する2粒子の位置の差 ξ\bm\xi は、右辺を1次まで展開して

ξ¨i=−∂i∂jΦ ξj(11.23)\ddot\xi^i=-\partial_i\partial_j\Phi\,\xi^j\tag{11.23}

に従う。行列 ∂i∂jΦ\partial_i\partial_j\Phi を潮汐テンソルと呼ぶ。

この式は (11.21) と同じ形をしている。実際、時刻 tt と位置 xix^i をまとめた4次元の空間に、0でない成分が

Γitt=∂iΦ(11.24)\Gamma^i{}_{tt}=\partial_i\Phi\tag{11.24}

だけの対称な接続を置くと、測地線の方程式 (11.2) は t¨=0\ddot t=0 と x¨i+∂iΦ=0\ddot x^i+\partial_i\Phi=0 になる。tt をパラメータに選べば、これは Newton の運動方程式そのものである。

この接続はどんな計量の Levi-Civita 接続でもない。

(t,x)(t,\bm x) の空間には計量を勝手に置かない(第5章 §5.5)が、定義 11.1 は接続だけで曲率を決めるので、そのまま計算できる。(11.3) に入れると、接続どうしの積はすべて 00 になる。積 ΓρμλΓλνσ\Gamma^\rho{}_{\mu\lambda}\Gamma^\lambda{}_{\nu\sigma} が0でないには、前の因子から λ=t\lambda=t が要る。ところが後ろの因子は上の添字が空間のときしか0でないので、両立しないからである。

命題 11.18(Newton 力学の曲率 ★★)

接続 (11.24) の曲率の0でない成分は

Ritjt=−Rittj=∂i∂jΦ(11.25)R^i{}_{tjt}=-R^i{}_{ttj}=\partial_i\partial_j\Phi\tag{11.25}

で、Ricci テンソルの0でない成分は Rtt=∇2Φ=4πGρR_{tt}=\nabla^2\Phi=4\pi G\rho だけである。

証明

Ritjt=∂jΓitt−∂tΓijt=∂j∂iΦR^i{}_{tjt}=\partial_j\Gamma^i{}_{tt}-\partial_t\Gamma^i{}_{jt}=\partial_j\partial_i\Phi である。上の添字が tt の成分は、Γtμν=0\Gamma^t{}_{\mu\nu}=0 から消える。2番目の添字が空間の添字 kk の成分も、Γiνk=0\Gamma^i{}_{\nu k}=0 から消える。Ricci は Rtt=Ritit=∑i∂i2ΦR_{tt}=R^i{}_{tit}=\sum_i\partial_i^2\Phi で、Poisson 方程式 ∇2Φ=4πGρ\nabla^2\Phi=4\pi G\rho を使えば最後の等号になる。

この接続で (11.21) を書く。同時刻の2粒子を比べるので J=(0,ξ)J=(0,\bm\xi)、速度は T=(1,v)T=(1,\bm v) である。右辺で0でないのは σ=μ=t, ν=j\sigma=\mu=t,\ \nu=j の項だけである。

RittjTtTtξj=−∂i∂jΦ ξj.R^i{}_{ttj}T^tT^t\xi^j=-\partial_i\partial_j\Phi\,\xi^j.

Jt=0J^t=0 なので D2ξi/dt2=ξ¨iD^2\xi^i/\mathrm{d} t^2=\ddot\xi^i で、Newton の (11.23) がそのまま出る。

この見方では、重力は「接続が平坦でないこと」で、潮汐はその曲率である。一様な重力(∂i∂jΦ=0\partial_i\partial_j\Phi=0)は曲率を持たず、自由落下する箱の中では消せる。一方、潮汐は座標の取り替えでは消せない。一般相対論では同じ役割を時空の計量の曲率 c2Ri0j0≃∂i∂jΦc^2R^i{}_{0j0}\simeq\partial_i\partial_j\Phi が担い、その導出は『相対論入門』第10章 §10.2 で行う。

点質量の潮汐

点の質量を MM とし、ポテンシャル Φ=−GM/r\Phi=-GM/r を考える。単位ベクトル n=x/r\bm n=\bm x/r を使うと、潮汐テンソルは次のとおりである。

∂i∂jΦ=GMr3(δij−3ninj)(11.26)\partial_i\partial_j\Phi=\frac{GM}{r^3}\bigl(\delta_{ij}-3n_in_j\bigr)\tag{11.26}

である。固有値は、質量の方向に −2GM/r3-2GM/r^3、それに垂直な2方向に +GM/r3+GM/r^3 である(確認テスト 11.4-4)。(11.23) の負号から、質量の方向には引き伸ばされ、垂直な方向には押し縮められる。跡は 00 で、物質の無い所では Rtt=0R_{tt}=0 であることに対応する。

図

図 11.4 地球の周りに輪に並べた自由な試験粒子が月の潮汐で受ける相対加速度(赤の矢印、(11.26))。月の方向の成分は +2GMx/d3+2GMx/d^3、垂直な成分は −GMy/d3-GMy/d^3 で、輪は月の方向に伸びた楕円(青)になる。満潮が地球の両側に1つずつできる理由である。

例題 11.5 地球と月の潮汐の見積もり ★★★

オリジナル次の値を使う(MmM_{\mathrm m} は月の質量、dd は地球と月の距離)。

G=6.674×10−11 m3 kg−1 s−2,g=9.81 m/s2,Mm=7.342×1022 kg,ME=5.972×1024 kg,d=3.844×108 m,RE=6.371×106 m.\begin{aligned} &G=6.674\times10^{-11}\ \mathrm{m^3\,kg^{-1}\,s^{-2}},\qquad g=9.81\ \mathrm{m/s^2},\qquad M_{\mathrm m}=7.342\times10^{22}\ \mathrm{kg},\\ &M_{\mathrm E}=5.972\times10^{24}\ \mathrm{kg},\qquad d=3.844\times10^8\ \mathrm m,\qquad R_{\mathrm E}=6.371\times10^6\ \mathrm m . \end{aligned}

  1. 地球の中心と、月に最も近い地表の点との、月の引力による相対加速度(潮汐加速度)を (11.26) から求めよ。月の引力そのものの大きさとも比べよ。
  2. 海面が潮汐のポテンシャルの等ポテンシャル面に沿うとして(平衡潮)、満潮と干潮の海面の高さの差を見積もれ。
解答

(1) 月の方向に間隔 ξ=RE\xi=R_{\mathrm E} をとると、(11.23) と (11.26) から相対加速度は 2GMmRE/d32GM_{\mathrm m}R_{\mathrm E}/d^3 である。数値を入れると

atide=2×6.674×10−11×7.342×1022×6.371×106(3.844×108)3 m/s2=1.10×10−6 m/s2a_{\mathrm{tide}}=\frac{2\times6.674\times10^{-11}\times7.342\times10^{22}\times6.371\times10^{6}}{(3.844\times10^8)^3}\ \mathrm{m/s^2} =1.10\times10^{-6}\ \mathrm{m/s^2}

で、gg の 1.1×10−71.1\times10^{-7} 倍である。月の引力そのものは GMm/d2=3.32×10−5 m/s2GM_{\mathrm m}/d^2=3.32\times10^{-5}\ \mathrm{m/s^2} で、その約 3%3\% にすぎない。地球全体が月に向かって自由落下しているので、引力そのものは打ち消され、場所による差だけが残る。

(2) (11.23) の右辺は、潮汐のポテンシャル Φtide=12∂i∂jΦ ξiξj\Phi_{\mathrm{tide}}=\tfrac12\partial_i\partial_j\Phi\,\xi^i\xi^j の勾配の −1-1 倍である。(11.26) を入れ、月の方向から測った角を ϑ\vartheta とすると、地表(∣ξ∣=RE\lvert\bm\xi\rvert=R_{\mathrm E})で

Φtide=GMmRE22d3(1−3cos⁡2ϑ)\Phi_{\mathrm{tide}}=\frac{GM_{\mathrm m}R_{\mathrm E}^2}{2d^3}\bigl(1-3\cos^2\vartheta\bigr)

となる。

月の真下 ϑ=0\vartheta=0 と真横 ϑ=π/2\vartheta=\pi/2 を比べる。ポテンシャルの差は 32GMmRE2/d3\tfrac32GM_{\mathrm m}R_{\mathrm E}^2/d^3 である。海面の高さの差 hh は、これを重力加速度 g=GME/RE2g=GM_{\mathrm E}/R_{\mathrm E}^2 で割ったものになる。

h=32 MmME(REd)3RE=1.5×0.01229×4.553×10−6×6.371×106 m=0.53 m.h=\frac32\,\frac{M_{\mathrm m}}{M_{\mathrm E}}\Bigl(\frac{R_{\mathrm E}}{d}\Bigr)^3R_{\mathrm E} =1.5\times0.01229\times4.553\times10^{-6}\times6.371\times10^6\ \mathrm m=0.53\ \mathrm m .

実際の潮位差は海岸の地形と海水の運動で大きく変わるが、外洋での大きさの目安はこの値である。

確認テスト

オリジナル

  1. 11.4-1 平面上の Jacobi 方程式 (11.22) の一般解を書き、それが2本の直線の間隔の式になっていることを確かめよ。
  2. 11.4-2 K=−1/ℓ2K=-1/\ell^2 の面で、j(0)=0, j′(0)=1j(0)=0,\ j'(0)=1 の Jacobi 場を求めよ。それを使い、1点から測地線の長さ rr にある点がなす円の周の長さを求め、平面の 2πr2\pi r と比べよ。
  3. 11.4-3 太陽の質量 Ms=1.989×1030 kgM_{\mathrm s}=1.989\times10^{30}\ \mathrm{kg}、地球と太陽の距離 ds=1.496×1011 md_{\mathrm s}=1.496\times10^{11}\ \mathrm m として、太陽による地表の潮汐加速度を求め、月によるものとの比を求めよ。太陽の引力そのものは月の引力の何倍か。また、月と太陽が一直線に並ぶとき(大潮)と直角に並ぶとき(小潮)の潮汐加速度の比を見積もれ。
  4. 11.4-4 (11.26) を導き、その固有値と固有ベクトルを求めよ。跡が 00 になることも確かめよ。

木下の勉強ステータス

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

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

← AI生成教科書『基礎理論』の記事一覧へ