物理のための微分幾何入門 第10章 §10.7
10.1-1 ϕ ( λ ) = 2 λ + 1 \phi(\lambda)=2\lambda+1 ϕ ( λ ) = 2 λ + 1 と ϕ ( λ ) = − λ \phi(\lambda)=-\lambda ϕ ( λ ) = − λ はどちらも1次関数である。命題 10.4 の右辺は0になるので、x ( 2 t + 1 ) x(2t+1) x ( 2 t + 1 ) と x ( − t ) x(-t) x ( − t ) は測地線の方程式をそのまま満たす。ϕ ′ < 0 \phi'<0 ϕ ′ < 0 の場合も、証明の式変形は ϕ ′ ≠ 0 \phi'\ne0 ϕ ′ = 0 しか使わないのでそのまま通る。
ϕ ( λ ) = λ 3 \phi(\lambda)=\lambda^3 ϕ ( λ ) = λ 3 では ϕ ′ = 3 λ 2 , ϕ ′ ′ = 6 λ \phi'=3\lambda^2,\ \phi''=6\lambda ϕ ′ = 3 λ 2 , ϕ ′′ = 6 λ である。(10.4) の右辺の係数は次のようになる。
ϕ ′ ′ ϕ ′ = 6 λ 3 λ 2 = 2 λ \frac{\phi''}{\phi'}=\frac{6\lambda}{3\lambda^2}=\frac2\lambda ϕ ′ ϕ ′′ = 3 λ 2 6 λ = λ 2
これは0でないので、x ( t 3 ) x(t^3) x ( t 3 ) は測地線の方程式を満たさない。曲線としては同じ測地線の道をたどるが、径数がアフィンでない。λ = 0 \lambda=0 λ = 0 では ϕ ′ = 0 \phi'=0 ϕ ′ = 0 となり、速度が0になる。2. 10.1-2 1. (a) g x x = g y y = y − 2 g_{xx}=g_{yy}=y^{-2} g xx = g y y = y − 2 で、0でない微分は ∂ y g x x = ∂ y g y y = − 2 y − 3 \partial_yg_{xx}=\partial_yg_{yy}=-2y^{-3} ∂ y g xx = ∂ y g y y = − 2 y − 3 である。(10.1) から次の3つ(と対称な成分)が0でない。
Γ x x y = 1 2 y 2 ∂ y g x x = − 1 y , Γ y x x = − 1 2 y 2 ∂ y g x x = 1 y , Γ y y y = 1 2 y 2 ∂ y g y y = − 1 y \Gamma^x{}_{xy}=\frac12y^2\,\partial_yg_{xx}=-\frac1y,\qquad
\Gamma^y{}_{xx}=-\frac12y^2\,\partial_yg_{xx}=\frac1y,\qquad
\Gamma^y{}_{yy}=\frac12y^2\,\partial_yg_{yy}=-\frac1y Γ x x y = 2 1 y 2 ∂ y g xx = − y 1 , Γ y xx = − 2 1 y 2 ∂ y g xx = y 1 , Γ y y y = 2 1 y 2 ∂ y g y y = − y 1
測地線の方程式は次の2本になる。
x ¨ − 2 y x ˙ y ˙ = 0 , y ¨ + x ˙ 2 − y ˙ 2 y = 0 \ddot x-\frac2y\,\dot x\dot y=0,\qquad \ddot y+\frac{\dot x^2-\dot y^2}{y}=0 x ¨ − y 2 x ˙ y ˙ = 0 , y ¨ + y x ˙ 2 − y ˙ 2 = 0
2. (b) x = x 0 , y = e s x=x_0,\ y=\mathrm{e}^s x = x 0 , y = e s では x ˙ = 0 \dot x=0 x ˙ = 0 で、1本めは満たされる。y ˙ = y ¨ = e s \dot y=\ddot y=\mathrm{e}^s y ˙ = y ¨ = e s なので、2本めの左辺は e s − e 2 s / e s = 0 \mathrm{e}^s-\mathrm{e}^{2s}/\mathrm{e}^s=0 e s − e 2 s / e s = 0 である。速さは ( x ˙ 2 + y ˙ 2 ) / y 2 = e 2 s / e 2 s = 1 (\dot x^2+\dot y^2)/y^2=\mathrm{e}^{2s}/\mathrm{e}^{2s}=1 ( x ˙ 2 + y ˙ 2 ) / y 2 = e 2 s / e 2 s = 1 なので、s s s は弧長である。y = 1 y=1 y = 1 から y = e y=\mathrm{e} y = e までの長さは次のとおりである。
∫ 1 e d y y = 1 \int_1^{\mathrm{e}}\frac{\mathrm{d} y}{y}=1 ∫ 1 e y d y = 1
3. (c) x ˙ = R / cosh 2 s \dot x=R/\cosh^2s x ˙ = R / cosh 2 s 、y ˙ = − R sinh s / cosh 2 s \dot y=-R\sinh s/\cosh^2s y ˙ = − R sinh s / cosh 2 s である。速さを調べると、次のとおり1になる。
x ˙ 2 + y ˙ 2 y 2 = R 2 ( 1 + sinh 2 s ) / cosh 4 s R 2 / cosh 2 s = 1 \frac{\dot x^2+\dot y^2}{y^2}=\frac{R^2(1+\sinh^2s)/\cosh^4s}{R^2/\cosh^2s}=1 y 2 x ˙ 2 + y ˙ 2 = R 2 / cosh 2 s R 2 ( 1 + sinh 2 s ) / cosh 4 s = 1
1本めの式は d d s ( x ˙ / y 2 ) = ( x ¨ y − 2 x ˙ y ˙ ) / y 3 \frac{\mathrm{d}}{\mathrm{d} s}(\dot x/y^2)=(\ddot x\,y-2\dot x\dot y)/y^3 d s d ( x ˙ / y 2 ) = ( x ¨ y − 2 x ˙ y ˙ ) / y 3 と同値で、x ˙ / y 2 = 1 / R \dot x/y^2=1/R x ˙ / y 2 = 1/ R が定数なので満たされる。2本めの式を調べる。y ¨ = − R ( 1 − sinh 2 s ) / cosh 3 s \ddot y=-R(1-\sinh^2s)/\cosh^3s y ¨ = − R ( 1 − sinh 2 s ) / cosh 3 s と計算でき、次のとおり0になる。
y ¨ + x ˙ 2 − y ˙ 2 y = − R ( 1 − sinh 2 s ) + R ( 1 − sinh 2 s ) cosh 3 s = 0 \ddot y+\frac{\dot x^2-\dot y^2}{y}=\frac{-R(1-\sinh^2s)+R(1-\sinh^2s)}{\cosh^3s}=0 y ¨ + y x ˙ 2 − y ˙ 2 = cosh 3 s − R ( 1 − sinh 2 s ) + R ( 1 − sinh 2 s ) = 0
3. 10.1-3 x = t 2 x=t^2 x = t 2 では速さは ∣ x ˙ ∣ = 2 t \left\vert \dot x\right\vert =2t ∣ x ˙ ∣ = 2 t である。長さとエネルギーは次のとおりである。
L = ∫ 0 1 2 t d t = 1 , E = 1 2 ∫ 0 1 4 t 2 d t = 2 3 L=\int_0^12t\,\mathrm{d} t=1,\qquad E=\frac12\int_0^14t^2\,\mathrm{d} t=\frac23 L = ∫ 0 1 2 t d t = 1 , E = 2 1 ∫ 0 1 4 t 2 d t = 3 2
L 2 = 1 ≤ 2 ( b − a ) E = 4 3 L^2=1\le2(b-a)E=\frac43 L 2 = 1 ≤ 2 ( b − a ) E = 3 4 で、不等式は成り立つ。x = t x=t x = t と走り直すと速さは1で一定になり、E = 1 2 E=\frac12 E = 2 1 である。L 2 = 1 = 2 ⋅ 1 ⋅ 1 2 L^2=1=2\cdot1\cdot\frac12 L 2 = 1 = 2 ⋅ 1 ⋅ 2 1 で等号になる。4. 10.1-4 y ( t ) = x ( c t ) y(t)=x(ct) y ( t ) = x ( c t ) とおくと、y ˙ ( t ) = c x ˙ ( c t ) , y ¨ ( t ) = c 2 x ¨ ( c t ) \dot y(t)=c\,\dot x(ct),\ \ddot y(t)=c^2\ddot x(ct) y ˙ ( t ) = c x ˙ ( c t ) , y ¨ ( t ) = c 2 x ¨ ( c t ) である。測地線の方程式に入れると c 2 c^2 c 2 がくくり出せる。
y ¨ k + Γ k i j ( y ) y ˙ i y ˙ j = c 2 ( x ¨ k + Γ k i j x ˙ i x ˙ j ) ∣ c t = 0 \ddot y^k+\Gamma^k{}_{ij}(y)\dot y^i\dot y^j=c^2\left(\ddot x^k+\Gamma^k{}_{ij}\dot x^i\dot x^j\right)\Big|_{ct}=0 y ¨ k + Γ k ij ( y ) y ˙ i y ˙ j = c 2 ( x ¨ k + Γ k ij x ˙ i x ˙ j ) c t = 0
よって y y y は測地線で、初速度は y ˙ ( 0 ) = c x ˙ ( 0 ) = c v \dot y(0)=c\,\dot x(0)=cv y ˙ ( 0 ) = c x ˙ ( 0 ) = c v である。命題 10.2 により、出発点と初速度が同じ測地線は1つしかない。したがって初速度 c v cv c v の測地線は y ( t ) = x ( c t ) y(t)=x(ct) y ( t ) = x ( c t ) そのものである。これはもとの測地線と同じ道を、∣ c ∣ \left\vert c\right\vert ∣ c ∣ 倍の速さで通る。5. 10.2-1 系 10.9 により、測地線になる緯線は f ′ ( u ) = 0 f'(u)=0 f ′ ( u ) = 0 の所である。f ′ ( u ) = − sin ( u / a ) f'(u)=-\sin(u/a) f ′ ( u ) = − sin ( u / a ) が0になるのは、u = 0 , π a u=0,\ \pi a u = 0 , π a の2つである。u = 0 u=0 u = 0 では f = R + a f=R+a f = R + a で、外側の赤道(軸から最も遠い円)である。u = π a u=\pi a u = π a では f = R − a f=R-a f = R − a で、内側の赤道(穴のふちの円)である。6. 10.2-2 1. (a) f = a sin θ f=a\sin\theta f = a sin θ (θ \theta θ は北極からの角)である。赤道では f = a , ψ = 60 ∘ f=a,\ \psi=60^\circ f = a , ψ = 6 0 ∘ なので、Clairaut の定数は J = a cos 60 ∘ = a / 2 J=a\cos60^\circ=a/2 J = a cos 6 0 ∘ = a /2 である。最も北の点では ψ = 0 \psi=0 ψ = 0 である。よって a sin θ min = a / 2 a\sin\theta_{\min}=a/2 a sin θ m i n = a /2 から θ min = 30 ∘ \theta_{\min}=30^\circ θ m i n = 3 0 ∘ となる。最高緯度は 90 ∘ − 30 ∘ = 60 ∘ 90^\circ-30^\circ=60^\circ 9 0 ∘ − 3 0 ∘ = 6 0 ∘ である。2. (b) 緯度 30 ∘ 30^\circ 3 0 ∘ の緯線では f = a sin 60 ∘ = a cos 30 ∘ f=a\sin60^\circ=a\cos30^\circ f = a sin 6 0 ∘ = a cos 3 0 ∘ である。f cos ψ = a / 2 f\cos\psi=a/2 f cos ψ = a /2 から次の値を得る。
cos ψ = 1 / 2 cos 30 ∘ = 1 3 = 0.577 , ψ = 54.7 ∘ \cos\psi=\frac{1/2}{\cos30^\circ}=\frac{1}{\sqrt3}=0.577,\qquad \psi=54.7^\circ cos ψ = cos 3 0 ∘ 1/2 = 3 1 = 0.577 , ψ = 54. 7 ∘
7. 10.2-3 1. (a) f ′ = u / c 2 + u 2 f'=u/\sqrt{c^2+u^2} f ′ = u / c 2 + u 2 なので、次のとおり1になる。
f ′ 2 + h ′ 2 = u 2 + c 2 c 2 + u 2 = 1 f'^2+h'^2=\frac{u^2+c^2}{c^2+u^2}=1 f ′2 + h ′2 = c 2 + u 2 u 2 + c 2 = 1
2. (b) 出発点では f = 1 + 4 = 5 f=\sqrt{1+4}=\sqrt5 f = 1 + 4 = 5 、ψ = 60 ∘ \psi=60^\circ ψ = 6 0 ∘ である。Clairaut の定数は J = 5 cos 60 ∘ = 1.118 J=\sqrt5\cos60^\circ=1.118 J = 5 cos 6 0 ∘ = 1.118 である。首では f = c = 1 < J f=c=1<J f = c = 1 < J なので、(10.12) から首には達しない。折り返すのは f ( u ) = J f(u)=J f ( u ) = J の所である。
1 + u 2 = 1.118 ⟹ u 2 = 1.25 − 1 = 0.25 , u = 0.5 \sqrt{1+u^2}=1.118\quad\Longrightarrow\quad u^2=1.25-1=0.25,\qquad u=0.5 1 + u 2 = 1.118 ⟹ u 2 = 1.25 − 1 = 0.25 , u = 0.5
首の手前の u = 0.5 u=0.5 u = 0.5 の緯線に接して折り返し、もとの側へ戻っていく。8. 10.2-4 接している点 P P P では、測地線の向きが緯線の向きと一致するので u ˙ = 0 \dot u=0 u ˙ = 0 である。速さが1なので f φ ˙ = ± 1 f\dot\varphi=\pm1 f φ ˙ = ± 1 、つまり φ ˙ 2 = 1 / f ( u 0 ) 2 \dot\varphi^2=1/f(u_0)^2 φ ˙ 2 = 1/ f ( u 0 ) 2 である。(10.10) の1本めから、P P P での u ¨ \ddot u u ¨ は次のように決まる。
u ¨ = f f ′ φ ˙ 2 = f ′ ( u 0 ) f ( u 0 ) \ddot u=ff'\dot\varphi^2=\frac{f'(u_0)}{f(u_0)} u ¨ = f f ′ φ ˙ 2 = f ( u 0 ) f ′ ( u 0 )
緯線が測地線でないので、系 10.9 により f ′ ( u 0 ) ≠ 0 f'(u_0)\ne0 f ′ ( u 0 ) = 0 で、u ¨ ≠ 0 \ddot u\ne0 u ¨ = 0 である。点 P P P では u ˙ = 0 , u ¨ ≠ 0 \dot u=0,\ \ddot u\ne0 u ˙ = 0 , u ¨ = 0 なので、母線方向の座標はそこで極値をとる。
f ′ ( u 0 ) > 0 f'(u_0)>0 f ′ ( u 0 ) > 0 なら極小で、測地線は u ≥ u 0 u\ge u_0 u ≥ u 0 の側(軸から遠い側)にある。f ′ ( u 0 ) < 0 f'(u_0)<0 f ′ ( u 0 ) < 0 なら極大で、測地線は u ≤ u 0 u\le u_0 u ≤ u 0 の側にある。どちらの場合も測地線は f f f の大きい側にあり、(10.12) の f ≥ ∣ J ∣ f\ge\left\vert J\right\vert f ≥ ∣ J ∣ と合っている。9. 10.3-1 定数 m m m を外して Euler–Lagrange 方程式を書く。
r ¨ − r φ ˙ 2 = 0 , d d t ( r 2 φ ˙ ) = r 2 φ ¨ + 2 r r ˙ φ ˙ = 0 , z ¨ = 0 \ddot r-r\dot\varphi^2=0,\qquad \frac{\mathrm{d}}{\mathrm{d} t}\left(r^2\dot\varphi\right)=r^2\ddot\varphi+2r\dot r\dot\varphi=0,\qquad \ddot z=0 r ¨ − r φ ˙ 2 = 0 , d t d ( r 2 φ ˙ ) = r 2 φ ¨ + 2 r r ˙ φ ˙ = 0 , z ¨ = 0
2本めを r 2 r^2 r 2 で割ると φ ¨ + 2 r r ˙ φ ˙ = 0 \ddot\varphi+\frac2r\dot r\dot\varphi=0 φ ¨ + r 2 r ˙ φ ˙ = 0 である。r ˙ φ ˙ \dot r\dot\varphi r ˙ φ ˙ の係数2は、Γ φ r φ + Γ φ φ r \Gamma^\varphi{}_{r\varphi}+\Gamma^\varphi{}_{\varphi r} Γ φ r φ + Γ φ φ r にあたる。0でない Christoffel 記号は、Γ r φ φ = − r \Gamma^r{}_{\varphi\varphi}=-r Γ r φφ = − r と Γ φ r φ = Γ φ φ r = 1 / r \Gamma^\varphi{}_{r\varphi}=\Gamma^\varphi{}_{\varphi r}=1/r Γ φ r φ = Γ φ φ r = 1/ r だけである。10. 10.3-2 1. (a) 円柱面の線素は R 2 d φ 2 + d z 2 R^2\mathrm{d}\varphi^2+\mathrm{d} z^2 R 2 d φ 2 + d z 2 である。座標 ( R φ , z ) (R\varphi,z) ( R φ , z ) をとれば、これは平面のデカルト座標の線素と同じ形になる。定理 10.10 により、質点はこの平面で等速直線運動をする。
初速度の成分は、軸方向が v cos 30 ∘ v\cos30^\circ v cos 3 0 ∘ 、周方向が v sin 30 ∘ v\sin30^\circ v sin 3 0 ∘ である。方位角も軸方向の位置も時刻に比例して増えるので、軌道は螺旋である。1周にかかる時間は 2 π R / ( v sin 30 ∘ ) 2\pi R/(v\sin30^\circ) 2 π R / ( v sin 3 0 ∘ ) なので、ピッチは次のとおりである。
v cos 30 ∘ × 2 π R v sin 30 ∘ = 2 π R cot 30 ∘ = 2 π × 0.10 × 1.732 m = 1.09 m v\cos30^\circ\times\frac{2\pi R}{v\sin30^\circ}=2\pi R\cot30^\circ=2\pi\times0.10\times1.732\ \mathrm m=1.09\ \mathrm m v cos 3 0 ∘ × v sin 3 0 ∘ 2 π R = 2 π R cot 3 0 ∘ = 2 π × 0.10 × 1.732 m = 1.09 m
2. (b) 空間での速さは v v v で一定なので、加速度は法線成分だけで、大きさは v 2 κ v^2\kappa v 2 κ である(第6章 §6.4 )。周方向の速さ v φ = v sin 30 ∘ v_\varphi=v\sin30^\circ v φ = v sin 3 0 ∘ の円運動の向心加速度は v φ 2 / R v_\varphi^2/R v φ 2 / R で、軸方向の加速度は0である。
κ = v φ 2 R v 2 = sin 2 30 ∘ R = 0.25 0.10 m = 2.5 m − 1 \kappa=\frac{v_\varphi^2}{Rv^2}=\frac{\sin^230^\circ}{R}=\frac{0.25}{0.10\ \mathrm m}=2.5\ \mathrm{m^{-1}} κ = R v 2 v φ 2 = R sin 2 3 0 ∘ = 0.10 m 0.25 = 2.5 m − 1
円柱面の上では測地線でも、空間の曲線としては曲がっている。主法線は軸を向き、曲面に垂直である。11. 10.3-3 1. (a) 単位質量あたりの運動エネルギーは T = 1 2 l 2 ( θ ˙ 2 + sin 2 θ φ ˙ 2 ) T=\frac12l^2(\dot\theta^2+\sin^2\theta\,\dot\varphi^2) T = 2 1 l 2 ( θ ˙ 2 + sin 2 θ φ ˙ 2 ) である。計量は g = l 2 d i a g ( 1 , sin 2 θ ) g=l^2\,\mathrm{diag}(1,\sin^2\theta) g = l 2 diag ( 1 , sin 2 θ ) で、0でない Christoffel 記号は次の2つである。
Γ θ φ φ = − sin θ cos θ , Γ φ θ φ = cot θ \Gamma^\theta{}_{\varphi\varphi}=-\sin\theta\cos\theta,\qquad \Gamma^\varphi{}_{\theta\varphi}=\cot\theta Γ θ φφ = − sin θ cos θ , Γ φ θ φ = cot θ
右辺は − g θ θ ∂ θ V = − 1 l 2 g l sin θ = − g l sin θ -g^{\theta\theta}\partial_\theta V=-\frac{1}{l^2}\,gl\sin\theta=-\frac gl\sin\theta − g θ θ ∂ θ V = − l 2 1 g l sin θ = − l g sin θ と、− g φ φ ∂ φ V = 0 -g^{\varphi\varphi}\partial_\varphi V=0 − g φφ ∂ φ V = 0 である。運動方程式は次の2本になる。
θ ¨ − sin θ cos θ φ ˙ 2 = − g l sin θ , φ ¨ + 2 cot θ θ ˙ φ ˙ = 0 \ddot\theta-\sin\theta\cos\theta\,\dot\varphi^2=-\frac gl\sin\theta,\qquad \ddot\varphi+2\cot\theta\,\dot\theta\dot\varphi=0 θ ¨ − sin θ cos θ φ ˙ 2 = − l g sin θ , φ ¨ + 2 cot θ θ ˙ φ ˙ = 0
2. (b) θ = θ 0 \theta=\theta_0 θ = θ 0 を保つなら θ ˙ = θ ¨ = 0 \dot\theta=\ddot\theta=0 θ ˙ = θ ¨ = 0 である。2本めから、角速度 φ ˙ = ω \dot\varphi=\omega φ ˙ = ω は一定になる。1本めは次の式になる。
sin θ 0 cos θ 0 ω 2 = g l sin θ 0 ⟹ ω = g l cos θ 0 \sin\theta_0\cos\theta_0\,\omega^2=\frac gl\sin\theta_0\quad\Longrightarrow\quad\omega=\sqrt{\frac{g}{l\cos\theta_0}} sin θ 0 cos θ 0 ω 2 = l g sin θ 0 ⟹ ω = l cos θ 0 g
l = 1.0 m , θ 0 = 30 ∘ l=1.0\ \mathrm m,\ \theta_0=30^\circ l = 1.0 m , θ 0 = 3 0 ∘ では ω = 9.8 / 0.866 s − 1 = 3.36 s − 1 \omega=\sqrt{9.8/0.866}\ \mathrm{s^{-1}}=3.36\ \mathrm{s^{-1}} ω = 9.8/0.866 s − 1 = 3.36 s − 1 である。周期は 2 π / ω = 1.87 s 2\pi/\omega=1.87\ \mathrm s 2 π / ω = 1.87 s となる。12. 10.3-4 命題 10.3 の証明の Leibniz 則で、2つのベクトル場をともに q ˙ \dot q q ˙ にとる。D q ˙ k / d t D\dot q^k/\mathrm{d} t D q ˙ k / d t は (10.15) の左辺なので、− g k l ∂ l V -g^{kl}\partial_lV − g k l ∂ l V に等しい。
d d t 1 2 g ( q ˙ , q ˙ ) = g ( D q ˙ d t , q ˙ ) = g k j ( − g k l ∂ l V ) q ˙ j = − ∂ j V q ˙ j = − d V d t \frac{\mathrm{d}}{\mathrm{d} t}\,\frac12g(\dot q,\dot q)=g\!\left(\frac{D\dot q}{\mathrm{d} t},\dot q\right)=g_{kj}\left(-g^{kl}\partial_lV\right)\dot q^j=-\partial_jV\,\dot q^j=-\frac{\mathrm{d} V}{\mathrm{d} t} d t d 2 1 g ( q ˙ , q ˙ ) = g ( d t D q ˙ , q ˙ ) = g k j ( − g k l ∂ l V ) q ˙ j = − ∂ j V q ˙ j = − d t d V
最後の等号は連鎖律である。移項すると d d t ( 1 2 g ( q ˙ , q ˙ ) + V ) = 0 \frac{\mathrm{d}}{\mathrm{d} t}\bigl(\frac12g(\dot q,\dot q)+V\bigr)=0 d t d ( 2 1 g ( q ˙ , q ˙ ) + V ) = 0 となる。13. 10.4-1 運動エネルギーの計量は g = m ( d x 2 + d y 2 ) g=m(\mathrm{d} x^2+\mathrm{d} y^2) g = m ( d x 2 + d y 2 ) なので、Jacobi 計量は次のとおりである。
g ~ = 2 m ( E − m g y ) ( d x 2 + d y 2 ) \tilde g=2m(E-mgy)\left(\mathrm{d} x^2+\mathrm{d} y^2\right) g ~ = 2 m ( E − m g y ) ( d x 2 + d y 2 )
計量として使えるのは係数が正の領域で、y < E / ( m g ) y<E/(mg) y < E / ( m g ) である。y = E / ( m g ) y=E/(mg) y = E / ( m g ) は速さが0になる高さで、そこで計量はつぶれる。14. 10.4-2 F ( y ) = 2 m ( E − m g y ) F(y)=2m(E-mgy) F ( y ) = 2 m ( E − m g y ) と書く。Jacobi 計量の弧長 s s s での微分を ′ ' ′ で表す。関数 F F F は x x x によらない。そのため、Jacobi 計量のエネルギー汎関数の x x x についての Euler–Lagrange 方程式から、J = F x ′ \mathcal J=F\,x' J = F x ′ が保存する。
弧長の条件 F ( x ′ 2 + y ′ 2 ) = 1 F(x'^2+y'^2)=1 F ( x ′2 + y ′2 ) = 1 を x ′ 2 x'^2 x ′2 で割る。そこへ x ′ = J / F x'=\mathcal J/F x ′ = J / F を入れると、次の式になる。
( d y d x ) 2 = 1 F x ′ 2 − 1 = F J 2 − 1 \left(\frac{\mathrm{d} y}{\mathrm{d} x}\right)^2=\frac{1}{Fx'^2}-1=\frac{F}{\mathcal J^2}-1 ( d x d y ) 2 = F x ′2 1 − 1 = J 2 F − 1
両辺を x x x で微分して 2 d y / d x 2\,\mathrm{d} y/\mathrm{d} x 2 d y / d x で割る。∂ y F = − 2 m 2 g \partial_yF=-2m^2g ∂ y F = − 2 m 2 g を使うと、次の式を得る。
d 2 y d x 2 = ∂ y F 2 J 2 = − m 2 g J 2 \frac{\mathrm{d}^2y}{\mathrm{d} x^2}=\frac{\partial_yF}{2\mathcal J^2}=-\frac{m^2g}{\mathcal J^2} d x 2 d 2 y = 2 J 2 ∂ y F = − J 2 m 2 g
(10.18) により d s = 2 ( E − V ) d t = ( F / m ) d t \mathrm{d} s=2(E-V)\mathrm{d} t=(F/m)\,\mathrm{d} t d s = 2 ( E − V ) d t = ( F / m ) d t なので、x ′ = m x ˙ / F x'=m\dot x/F x ′ = m x ˙ / F である。よって J = F x ′ = m x ˙ \mathcal J=F\,x'=m\dot x J = F x ′ = m x ˙ となる。
これを入れると d 2 y / d x 2 = − g / x ˙ 2 \mathrm{d}^2y/\mathrm{d} x^2=-g/\dot x^2 d 2 y / d x 2 = − g / x ˙ 2 である。放物運動 x = v x t , y = v y t − 1 2 g t 2 x=v_xt,\ y=v_yt-\frac12gt^2 x = v x t , y = v y t − 2 1 g t 2 を x x x の関数として書くと、次のようになる。
y = v y v x x − g 2 v x 2 x 2 y=\frac{v_y}{v_x}x-\frac{g}{2v_x^2}x^2 y = v x v y x − 2 v x 2 g x 2
2階微分は − g / v x 2 -g/v_x^2 − g / v x 2 で、x ˙ = v x \dot x=v_x x ˙ = v x なので両者は一致する。15. 10.4-3 Jacobi 計量は F ( r ) ( d x 2 + d y 2 ) F(r)(\mathrm{d} x^2+\mathrm{d} y^2) F ( r ) ( d x 2 + d y 2 ) で、e 2 σ = F \mathrm{e}^{2\sigma}=F e 2 σ = F である。σ = 1 2 ln F \sigma=\frac12\ln F σ = 2 1 ln F は動径だけの関数なので、ラプラシアンは Δ σ = σ ′ ′ + σ ′ / r \Delta\sigma=\sigma''+\sigma'/r Δ σ = σ ′′ + σ ′ / r である。F = A + B / r F=A+B/r F = A + B / r とおき、A = 2 m E , B = 2 m k A=2mE,\ B=2mk A = 2 m E , B = 2 mk として計算すると、次の結果になる。
K = − 1 F ( σ ′ ′ + σ ′ r ) = − A B 2 ( A r + B ) 3 = − E k 4 m ( E r + k ) 3 K=-\frac{1}{F}\left(\sigma''+\frac{\sigma'}{r}\right)=-\frac{AB}{2(Ar+B)^3}=-\frac{Ek}{4m(Er+k)^3} K = − F 1 ( σ ′′ + r σ ′ ) = − 2 ( A r + B ) 3 A B = − 4 m ( E r + k ) 3 E k
運動できる領域では F > 0 F>0 F > 0 、つまり E r + k > 0 Er+k>0 E r + k > 0 である。したがって K K K の符号は − E -E − E の符号と同じになる。
E < 0 E<0 E < 0 (楕円軌道)なら K > 0 K>0 K > 0 である。
E = 0 E=0 E = 0 (放物線)なら K = 0 K=0 K = 0 である。
E > 0 E>0 E > 0 (双曲線)なら K < 0 K<0 K < 0 である。
E = 0 E=0 E = 0 で平らになることは、本文の exbox の結果と合っている。16. 10.4-4 実際の運動ではエネルギーの保存から 1 2 g i j q ˙ i q ˙ j = E − V \frac12g_{ij}\dot q^i\dot q^j=E-V 2 1 g ij q ˙ i q ˙ j = E − V である。これを使うと次の式を得る。
p i q ˙ i = g i j q ˙ j q ˙ i = 2 ( E − V ) p_i\dot q^i=g_{ij}\dot q^j\dot q^i=2(E-V) p i q ˙ i = g ij q ˙ j q ˙ i = 2 ( E − V )
Jacobi 計量での長さの要素は d s = g ~ i j q ˙ i q ˙ j d t \mathrm{d} s=\sqrt{\tilde g_{ij}\dot q^i\dot q^j}\,\mathrm{d} t d s = g ~ ij q ˙ i q ˙ j d t である。g ~ = 2 ( E − V ) g \tilde g=2(E-V)g g ~ = 2 ( E − V ) g を入れると次のようになる。
d s = 2 ( E − V ) ⋅ 2 ( E − V ) d t = 2 ( E − V ) d t = p i q ˙ i d t \mathrm{d} s=\sqrt{2(E-V)\cdot2(E-V)}\,\mathrm{d} t=2(E-V)\,\mathrm{d} t=p_i\dot q^i\,\mathrm{d} t d s = 2 ( E − V ) ⋅ 2 ( E − V ) d t = 2 ( E − V ) d t = p i q ˙ i d t
これは (10.18) とも一致している。17. 10.5-1 n n n が定数なら、光学的な計量 n 2 δ i j n^2\delta_{ij} n 2 δ ij はユークリッド計量の定数倍である。定数倍は Christoffel 記号を変えないので、測地線は直線である。光路長は n l nl n l 、かかる時間は n l / c nl/c n l / c である。n = 1.5 n=1.5 n = 1.5 、l = 2.0 m l=2.0\ \mathrm m l = 2.0 m では次のとおりである。
n l = 3.0 m , n l c = 3.0 3.00 × 10 8 s = 1.0 × 10 − 8 s nl=3.0\ \mathrm m,\qquad \frac{nl}{c}=\frac{3.0}{3.00\times10^8}\ \mathrm s=1.0\times10^{-8}\ \mathrm s n l = 3.0 m , c n l = 3.00 × 1 0 8 3.0 s = 1.0 × 1 0 − 8 s
18. 10.5-2 境界面を水平にとり、法線からの入射角を i i i とすると θ = 90 ∘ − i \theta=90^\circ-i θ = 9 0 ∘ − i である。(10.22) は n 1 sin i 1 = n 2 sin i 2 n_1\sin i_1=n_2\sin i_2 n 1 sin i 1 = n 2 sin i 2 となる。
sin i 2 = 1.00 × sin 45 ∘ 1.33 = 0.532 , i 2 = 32.1 ∘ \sin i_2=\frac{1.00\times\sin45^\circ}{1.33}=0.532,\qquad i_2=32.1^\circ sin i 2 = 1.33 1.00 × sin 4 5 ∘ = 0.532 , i 2 = 32. 1 ∘
19. 10.5-3 光線の方程式 (10.21) の x x x 成分を、光線の全長にわたって積分する。曲がりが小さいので d l ≈ d z \mathrm{d} l\approx\mathrm{d} z d l ≈ d z で、n ≈ 1 n\approx1 n ≈ 1 とみなせる。横向きの方向余弦の変化が曲がりの角 δ \delta δ になり、それは ∂ x n \partial_xn ∂ x n の積分で決まる。x = b x=b x = b の直線の上では r = b 2 + z 2 r=\sqrt{b^2+z^2} r = b 2 + z 2 である。
δ = ∫ − ∞ ∞ ∂ b ( 2 G M c 2 b 2 + z 2 ) d z = − 2 G M b c 2 ∫ − ∞ ∞ d z ( b 2 + z 2 ) 3 / 2 = − 2 G M b c 2 ⋅ 2 b 2 = − 4 G M c 2 b \delta=\int_{-\infty}^{\infty}\partial_b\left(\frac{2GM}{c^2\sqrt{b^2+z^2}}\right)\mathrm{d} z=-\frac{2GMb}{c^2}\int_{-\infty}^{\infty}\frac{\mathrm{d} z}{(b^2+z^2)^{3/2}}=-\frac{2GMb}{c^2}\cdot\frac{2}{b^2}=-\frac{4GM}{c^2b} δ = ∫ − ∞ ∞ ∂ b ( c 2 b 2 + z 2 2 GM ) d z = − c 2 2 GM b ∫ − ∞ ∞ ( b 2 + z 2 ) 3/2 d z = − c 2 2 GM b ⋅ b 2 2 = − c 2 b 4 GM
負号は、光線が x x x の小さい側、つまり質量の側へ曲がることを表す。太陽の縁をかすめる光では次の値になる。
∣ δ ∣ = 4 × 1.327 × 10 20 ( 3.00 × 10 8 ) 2 × 6.96 × 10 8 = 8.49 × 10 − 6 r a d = 1.75 ′ ′ \left\vert \delta\right\vert =\frac{4\times1.327\times10^{20}}{(3.00\times10^8)^2\times6.96\times10^8}=8.49\times10^{-6}\ \mathrm{rad}=1.75'' ∣ δ ∣ = ( 3.00 × 1 0 8 ) 2 × 6.96 × 1 0 8 4 × 1.327 × 1 0 20 = 8.49 × 1 0 − 6 rad = 1.7 5 ′′
20. 10.5-4 位置ベクトルと n t n\bm t n t の外積 M = x × n t \bm M=\bm x\times n\bm t M = x × n t を考える。∣ M ∣ = n r sin ϕ \left\vert \bm M\right\vert =n\,r\sin\phi ∣ M ∣ = n r sin ϕ である。l l l で微分すると、d x / d l = t \mathrm{d}\bm x/\mathrm{d} l=\bm t d x / d l = t と光線の方程式から次の式を得る。
d M d l = t × n t + x × d d l ( n t ) = 0 + x × ∇ n \frac{\mathrm{d}\bm M}{\mathrm{d} l}=\bm t\times n\bm t+\bm x\times\frac{\mathrm{d}}{\mathrm{d} l}(n\bm t)=\bm 0+\bm x\times\nabla n d l d M = t × n t + x × d l d ( n t ) = 0 + x × ∇ n
屈折率が動径だけの関数なら ∇ n = n ′ ( r ) x / r \nabla n=n'(r)\,\bm x/r ∇ n = n ′ ( r ) x / r で、これは x \bm x x に平行である。よって x × ∇ n = 0 \bm x\times\nabla n=\bm 0 x × ∇ n = 0 となり、M \bm M M は一定である。とくにその大きさ n r sin ϕ n\,r\sin\phi n r sin ϕ も一定になる。
この保存量は、回転の対称性から出たもので、力学の角運動量にあたる。M \bm M M の向きも一定なので、光線は原点を含む1つの平面の中を進む。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com