tecktom.com

この章のまとめ

2026年10月15日10分で読めますAI生成教科書『虎の巻』第2部 第10章 光線の解析力学

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第2部 第10章 §2.10.7 ★★★ 光線の力学は、光路長を作用に選んだ解析力学そのもの。Snell・アイコナール・ABCD・エタンデュは Noether・Hamilton–Jacobi・正準変換・Liouville の姿

光線は光路長 ℓ=∫n ds\ell=\int n\,\mathrm{d} s を停留させる経路を通る(Fermat の原理)。光が進む向き zz を時刻とみると Lagrangian は L=n1+x′2+y′2L=n\sqrt{1+x'^2+y'^2} で、Euler–Lagrange 方程式は光線の方程式 dds(n dr/ds)=∇n\frac{\mathrm{d} }{\mathrm{d} s}(n\,\mathrm{d}\bm{r}/\mathrm{d} s)=\nabla n になる。Legendre 変換で得る正準運動量は p=n dr/ds\bm{p}=n\,\mathrm{d}\bm{r}/\mathrm{d} s の横の成分、Hamiltonian は H=−n2−px2−py2H=-\sqrt{n^2-p_x^2-p_y^2} である。屈折率が層をなす媒質では並進の Noether の保存量 nsin⁡θn\sin\theta が一定で、これが Snell の法則と全反射・蜃気楼を与える。光路長を終点の関数とみたアイコナール S\mathcal S は ∣∇S∣2=n2\lvert \nabla\mathcal S\rvert ^2=n^2 を満たし、これは光線の HH の Hamilton–Jacobi 方程式である。Helmholtz 方程式に aeik0Sa\mathrm{e}^{\mathrm{i} k_0\mathcal S} を入れて λ0→0\lambda_0\to0 とすると同じ式が出て、ℏ↔1/k0\hbar\leftrightarrow1/k_0 の置き換えで Schrödinger 方程式の古典極限と重なる。光軸の近くでは HH が2次式になり、伝搬・屈折・レンズは行列式 11 の ABCD 行列(線形正準変換)で、結像の条件は B=0B=0、倍率と角の倍率の積は 11 である。屈折率が放物型の媒質は調和振動子で、振幅の大きい光線では周期が縮む。相空間の体積 n2cos⁡θ dA dΩn^2\cos\theta\,\mathrm{d} A\,\mathrm{d}\Omega(エタンデュ)は Liouville の定理で保たれ、明るさ Le/n2L_{\mathrm e}/n^2 と集光の上限 n2/sin⁡2θsn^2/\sin^2\theta_{\mathrm s} を決める。

光の巻(光学と量子エレクトロニクス 第2章)は、この章の原理の上に結像と光線追跡の実際を、光学と量子エレクトロニクス 第3章は HH を近軸から4次まで展開した残りとしての収差を扱う。光学と量子エレクトロニクス 第4章の近軸波動方程式は、2次の HH (2.10.17) を ℏ→1/k0\hbar\to1/k_0 で量子化した Schrödinger 方程式と同じ形をしている。重力による光の曲がりは、曲がった時空での Fermat の原理として 第2部第16章で扱う。

この章の道具

  • Fermat の原理:δ∫n ds=0\delta\int n\,\mathrm{d} s=0。zz を時刻にとると L=n1+x′2+y′2L=n\sqrt{1+x'^2+y'^2}。
  • 光線の方程式:dds(n dr/ds)=∇n\frac{\mathrm{d} }{\mathrm{d} s}(n\,\mathrm{d}\bm{r}/\mathrm{d} s)=\nabla n、径数 dσ=ds/n\mathrm{d}\sigma=\mathrm{d} s/n では d2r/dσ2=∇(n2/2)\mathrm{d}^2\bm{r}/\mathrm{d}\sigma^2=\nabla(n^2/2)。
  • 光線の Hamiltonian:p=n dr/ds\bm{p}=n\,\mathrm{d}\bm{r}/\mathrm{d} s、H=−n2−px2−py2=−pzH=-\sqrt{n^2-p_x^2-p_y^2}=-p_z。層状の媒質で nsin⁡θ=n\sin\theta= 一定(Snell)。
  • アイコナール方程式 ∣∇S∣2=n2\lvert \nabla\mathcal S\rvert ^2=n^2 は光線の Hamilton–Jacobi 方程式。光線は波面に直交。
  • 近軸:H≈−n0+p2/(2n0)+n2x2/2H\approx-n_0+p^2/(2n_0)+n_2x^2/2、ABCD 行列は AD−BC=1AD-BC=1。伝搬 (1d/n01)\bigl(\begin{smallmatrix}1&d/n\\0&1\end{smallmatrix}\bigr)、レンズ (10−1/f1)\bigl(\begin{smallmatrix}1&0\\-1/f&1\end{smallmatrix}\bigr)。
  • エタンデュ n2cos⁡θ dA dΩ=dx dy dpx dpyn^2\cos\theta\,\mathrm{d} A\,\mathrm{d}\Omega=\mathrm{d} x\,\mathrm{d} y\,\mathrm{d} p_x\,\mathrm{d} p_y は保存(Liouville)。Le/n2L_{\mathrm e}/n^2 が一定、集光の上限 n2/sin⁡2θsn^2/\sin^2\theta_{\mathrm s}。

総合問題 屈折率が放物型のファイバーの中の光線を5通りに見る ★★★

オリジナル屈折率が n(x)2=n02(1−α2x2)n(x)^2=n_0^2(1-\alpha^2x^2) の媒質(n0=1.5n_0=1.5, α=1.0 mm−1\alpha=1.0\,\mathrm{mm^{-1}})の中を、xzxz 面内で zz の正の向きに進む光線を考える。

  1. (1) (使う道具:2.10.4節・第1部第12章)近軸の Hamiltonian (2.10.22) で、長さ zz 進む ABCD 行列を求め、行列式が 11 であることを確かめよ。長さ π/(2α)\pi/(2\alpha)(4分の1周期)の行列が、平行な光の束を1点に集めるレンズとして働くことを示せ。
  2. (2) (使う道具:2.10.2節・2.10.4節)光軸上から出て振幅 xmax⁡=0.2 mmx_{\max}=0.2\,\mathrm{mm} に達する光線について、近軸の周期と厳密な周期 (2.10.23) を求めよ。光軸での光線の傾き θ0\theta_0(sin⁡θ0=p/n0\sin\theta_0=p/n_0)も求めよ。
  3. (3) (使う道具:第1部第14章の作用変数と断熱不変量)近軸の光線の作用変数 J=12π∮p dxJ=\frac{1}{2\pi}\oint p\,\mathrm{d} x を xmax⁡x_{\max} で表せ。長い距離をかけて α\alpha を 1.0 mm−11.0\,\mathrm{mm^{-1}} から 0.25 mm−10.25\,\mathrm{mm^{-1}} までゆっくり変えたファイバーの出口で、振幅 xmax⁡x_{\max} と最大の運動量 pmax⁡p_{\max} を求めよ。
  4. (4) (使う道具:2.10.5節)(3) の光線の相空間の軌道が囲む面積が 2πJ2\pi J に等しいことを確かめ、(3) の結果をエタンデュの保存として読め。
  5. (5) (使う道具:2.10.3節・第1部第14章の Bohr–Sommerfeld の条件)波の位相が k0Sk_0\mathcal S で ℏ\hbar の役を 1/k01/k_0 がすることから、J=(m+12)/k0J=(m+\frac12)/k_0(m=0,1,2,…m=0,1,2,\dots)の光線だけが波として許されると考える。真空中の波長 λ0=1.0 μm\lambda_0=1.0\,\mu\mathrm m のとき、(2) の光線はおよそ何番目の mm に当たるか。mm 番目の波が zz の向きに進むときの位相の増え方 βm\beta_m(伝搬定数)を求めよ。
解答

(1) 式 (2.10.22) の Hamilton 方程式は x′=p/n0x'=p/n_0, p′=−n0α2xp'=-n_0\alpha^2x で、調和振動子の流れなので

M(z)=(cos⁡αzsin⁡αzn0α−n0αsin⁡αzcos⁡αz),det⁡M=cos⁡2αz+sin⁡2αz=1.M(z)=\begin{pmatrix}\cos\alpha z&\dfrac{\sin\alpha z}{n_0\alpha}\\[1ex]-n_0\alpha\sin\alpha z&\cos\alpha z\end{pmatrix},\qquad\det M=\cos^2\alpha z+\sin^2\alpha z=1 .

z=π/(2α)z=\pi/(2\alpha) では M=(01/(n0α)−n0α0)M=\bigl(\begin{smallmatrix}0&1/(n_0\alpha)\\-n_0\alpha&0\end{smallmatrix}\bigr) で、p1=0p_1=0(平行な光の束)なら x2=0x_2=0 になり、高さによらずすべての光線が出口の光軸上の1点に集まる。出口の面に焦点のあるレンズとして働き、p2=−n0αx1p_2=-n_0\alpha x_1 から焦点距離に当たる量は 1/(n0α)1/(n_0\alpha)(−1/C-1/C)である。

(2) 近軸の周期は 2π/α=6.283 mm2\pi/\alpha=6.283\,\mathrm{mm}。厳密には αxmax⁡=0.2\alpha x_{\max}=0.2 なので 2π1−0.04/α=6.156 mm2\pi\sqrt{1-0.04}/\alpha=6.156\,\mathrm{mm} で、約 2%2\% 短い。光軸での運動量は厳密な解で p=dx/dσ=n0αxmax⁡=0.30p=\mathrm{d} x/\mathrm{d}\sigma=n_0\alpha x_{\max}=0.30 で、sin⁡θ0=0.30/1.5=0.20\sin\theta_0=0.30/1.5=0.20、θ0=0.201 rad\theta_0=0.201\,\mathrm{rad}(11.5∘11.5^\circ)。

(3) 近軸では x=xmax⁡sin⁡αzx=x_{\max}\sin\alpha z, p=n0x′=n0αxmax⁡cos⁡αzp=n_0x'=n_0\alpha x_{\max}\cos\alpha z で、相空間の軌道は半軸 xmax⁡x_{\max} と n0αxmax⁡n_0\alpha x_{\max} の楕円である。∮p dx\oint p\,\mathrm{d} x は楕円の面積 π xmax⁡⋅n0αxmax⁡\pi\,x_{\max}\cdot n_0\alpha x_{\max} なので J=n0αxmax⁡2/2J=n_0\alpha x_{\max}^2/2(HH から定数 −n0-n_0 を除いた値 EE を使えば J=E/αJ=E/\alpha とも書ける)。はじめは J=1.5×1.0×0.22/2=0.030 mmJ=1.5\times1.0\times0.2^2/2=0.030\,\mathrm{mm}。α\alpha を変える長さが周期 2π/α2\pi/\alpha よりずっと長ければ JJ は断熱不変量で保たれるので、xmax⁡∝α−1/2x_{\max}\propto\alpha^{-1/2}、pmax⁡=n0αxmax⁡∝α1/2p_{\max}=n_0\alpha x_{\max}\propto\alpha^{1/2}。α\alpha が 14\frac14 倍になると xmax⁡=0.40 mmx_{\max}=0.40\,\mathrm{mm}、pmax⁡=0.30×12=0.15p_{\max}=0.30\times\frac12=0.15。光の束は太くなり、角の広がりは半分になる。

(4) 楕円の面積は πxmax⁡pmax⁡=πn0αxmax⁡2=2πJ\pi x_{\max}p_{\max}=\pi n_0\alpha x_{\max}^2=2\pi J。はじめ π×0.2×0.30=0.06π\pi\times0.2\times0.30=0.06\pi、出口で π×0.40×0.15=0.06π mm\pi\times0.40\times0.15=0.06\pi\,\mathrm{mm} で等しい。相空間の面積の保存(Liouville の定理)はゆっくり変える必要のない厳密な性質で、断熱不変量は「1本の光線の楕円が、形を変えながらも別の楕円に移る」ことを加えて主張している。同じ面積の楕円の中で太さと角の広がりを交換したのが (3) の結果である。

(5) k0=2π/λ0=6283 mm−1k_0=2\pi/\lambda_0=6283\,\mathrm{mm^{-1}} なので k0J=6283×0.030=188.5k_0J=6283\times0.030=188.5、m+12=188.5m+\frac12=188.5 から m≈188m\approx188。Hamiltonian の値は H=−n0+αJ=−n0+(m+12)α/k0H=-n_0+\alpha J=-n_0+(m+\frac12)\alpha/k_0(調和振動子の E=αJE=\alpha J)で、時間発展の位相 e−iHt/ℏ\mathrm{e}^{-\mathrm{i} Ht/\hbar} の t→zt\to z, ℏ→1/k0\hbar\to1/k_0 が e−ik0Hz\mathrm{e}^{-\mathrm{i} k_0Hz} なので、位相は zz について

βm=−k0H=k0n0−(m+12)α\beta_m=-k_0H=k_0n_0-\Bigl(m+\frac12\Bigr)\alpha

の割合で増える。隣り合う mm の差は α=1.0 mm−1\alpha=1.0\,\mathrm{mm^{-1}} で、どの mm の波も zz が 2π/α2\pi/\alpha 進むごとに位相の差が 2π2\pi の整数倍に戻る。近軸の光線が周期 2π/α2\pi/\alpha で同じ形に戻ることの、波の側の言い換えである。

木下の勉強ステータス

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

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

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