tecktom.com

出てくるもの(1):Fraunhofer 回折は開口の Fourier 変換

2026年10月15日8分で読めますAI生成教科書『虎の巻』第6章 回折
目次
  1. 遠方の場は角スペクトルそのもの
  2. スリットと回折格子
  3. 木下の層

AI 生成・木下が未勉強

光学と量子エレクトロニクス 第6章 §6.3 ★★★ NF≪1N_{\mathrm F}\ll1 の遠方では、向き θ\theta に届く場は横の波数 q=ksin⁡θ\bm{q}=k\sin\theta のモードの振幅そのもの。スリットと回折格子の像がここから出る

遠方の場は角スペクトルそのもの

Fresnel 近似の式 (6.10) の指数の中を ∣x−x′∣2=∣x∣2−2x⋅x′+∣x′∣2\lvert \bm{x}-\bm{x}'\rvert ^2=\lvert \bm{x}\rvert ^2-2\bm{x}\cdot\bm{x}'+\lvert \bm{x}'\rvert ^2 と分ける。開口の大きさを aa として、開口の Fresnel 数が小さい(NF=a2/(λz)≪1N_{\mathrm F}=a^2/(\lambda z)\ll1、つまり k∣x′∣2/(2z)≪1k\lvert \bm{x}'\rvert ^2/(2z)\ll1)なら、∣x′∣2\lvert \bm{x}'\rvert ^2 の項を落としてよい。残りは

U(x,z)≈eikz eik∣x∣2/(2z)iλz U~(q=kxz;0),I(x)I0=1λ2z2∣U~(kxz;0)∣2(6.13)U(\bm{x},z)\approx\frac{\mathrm{e}^{\mathrm{i} kz}\,\mathrm{e}^{\mathrm{i} k\lvert \bm{x}\rvert ^2/(2z)}}{\mathrm{i}\lambda z}\,\tilde U\Bigl(\bm{q}=\frac{k\bm{x}}{z};0\Bigr),\qquad \frac{I(\bm{x})}{I_0}=\frac{1}{\lambda^2z^2}\Bigl\lvert\tilde U\Bigl(\frac{k\bm{x}}{z};0\Bigr)\Bigr\rvert^2\tag{6.13}

となる。これを Fraunhofer 回折と呼ぶ。遠方の像は、開口の面の上の場の Fourier 変換(角スペクトル)を、横の波数 q=kx/z\bm{q}=k\bm{x}/z で見たものである。

この結果の意味は、モードの言葉で読むとはっきりする。遠くの点 x\bm{x} は開口から向き x/z≈sin⁡θ\bm{x}/z\approx\sin\theta に見えるので、そこには開口から sin⁡θ\sin\theta の向きに進む平面波だけが届く。その平面波は横の波数 q=ksin⁡θ\bm{q}=k\sin\theta のモードで、振幅は角スペクトル U~(q;0)\tilde U(\bm{q};0) である。遠くに離れると、重ね合わせていた平面波のモードが向きごとに分かれて見える。近軸の極限 ∣x∣≪z\lvert \bm{x}\rvert \ll z では x/z=tan⁡θ≈sin⁡θ\bm{x}/z=\tan\theta\approx\sin\theta で、この差は6.5節で数値で見る。

全パワーの保存も確かめておく。d2x=(z/k)2d2q\mathrm{d}^2x=(z/k)^2\mathrm{d}^2q なので、式 (6.13) の強度を面の上で積分すると∫d2x I/I0=∫d2q(2π)2∣U~∣2=∫d2x′ ∣U(x′,0)∣2\int\mathrm{d}^2x\,I/I_0=\int\frac{\mathrm{d}^2q}{(2\pi)^2}\lvert \tilde U\rvert ^2=\int\mathrm{d}^2x'\,\lvert U(\bm{x}',0)\rvert ^2(最後は Parseval の等式)となり、開口を通ったパワーがそのまま遠方に届く。レンズの焦点面でも、距離 zz を焦点距離 ff に置き換えた同じ式が、遠方まで行かずに成り立つ。レンズの位相が ∣x′∣2\lvert \bm{x}'\rvert ^2 の項を打ち消すからで、この性質を使う光学(4f 系・空間フィルタ)は第7章(第7章)で扱う。

スリットと回折格子

yy の向きに長い幅 bb のスリット(∣x′∣<b/2\lvert x'\rvert <b/2 で透過率 11)では、xx の向きの変換だけが要り、

t~(qx)=∫−b/2b/2e−iqxx′ dx′=b sin⁡(qxb/2)qxb/2,I(θ)∝(sin⁡uu)2,u=πbsin⁡θλ(6.14)\tilde t(q_x)=\int_{-b/2}^{b/2}\mathrm{e}^{-\mathrm{i} q_xx'}\,\mathrm{d} x'=b\,\frac{\sin(q_xb/2)}{q_xb/2},\qquad I(\theta)\propto\Bigl(\frac{\sin u}{u}\Bigr)^2,\quad u=\frac{\pi b\sin\theta}{\lambda}\tag{6.14}

となる。暗い所は bsin⁡θ=mλb\sin\theta=m\lambda(m=±1,±2,…m=\pm1,\pm2,\dots)で、中央の明るい帯の半分の角度幅 λ/b\lambda/b は、スリットが狭いほど広い。Δq∼1/b\Delta q\sim1/b の見積もりがそのまま正確な形になった。

同じスリットを間隔 dd で NN 本並べると、nn 番目のスリットの透過の関数は x′x' を ndnd ずらしたものなので、Fourier 変換は e−inqxd\mathrm{e}^{-\mathrm{i} nq_xd} 倍になる。全部を足すと

t~N(qx)=t~(qx)∑n=0N−1e−inqxd,∣∑n=0N−1e−2inγ∣2=sin⁡2Nγsin⁡2γ,γ=qxd2=πdsin⁡θλ(6.15)\tilde t_N(q_x)=\tilde t(q_x)\sum_{n=0}^{N-1}\mathrm{e}^{-\mathrm{i} nq_xd},\qquad \Bigl\lvert\sum_{n=0}^{N-1}\mathrm{e}^{-2\mathrm{i} n\gamma}\Bigr\rvert^2=\frac{\sin^2N\gamma}{\sin^2\gamma},\quad \gamma=\frac{q_xd}{2}=\frac{\pi d\sin\theta}{\lambda}\tag{6.15}

を得る。強度は1本のスリットの像(包絡線)に、並べたことによる因子 sin⁡2Nγ/sin⁡2γ\sin^2N\gamma/\sin^2\gamma を掛けたものである。この因子は γ=mπ\gamma=m\pi、つまり

dsin⁡θ=mλ(m=0,±1,±2,… )(6.16)d\sin\theta=m\lambda\qquad(m=0,\pm1,\pm2,\dots)\tag{6.16}

で最大値 N2N^2 をとる(回折格子の式)。mm を回折の次数と呼ぶ。主な最大の山から最初の 00 までは Δγ=π/N\Delta\gamma=\pi/N、角度で Δ(sin⁡θ)=λ/(Nd)\Delta(\sin\theta)=\lambda/(Nd) で、本数が多いほど山は細い。山の高さが N2N^2 で幅が 1/N1/N なので、各山のパワーは NN に比例し、本数に比例して光が通ることと合っている。

格子は波長を分ける道具になる。波長 λ\lambda と λ+Δλ\lambda+\Delta\lambda の2本の線を mm 次で見ると、山の位置は Δ(sin⁡θ)=mΔλ/d\Delta(\sin\theta)=m\Delta\lambda/d だけずれる。このずれが山の半分の幅 λ/(Nd)\lambda/(Nd) に等しいとき、一方の山の頂上がもう一方の最初の 00 に重なり、2本がちょうど分かれて見える(Rayleigh の基準)。したがって分解能は

λΔλ=mN(6.17)\frac{\lambda}{\Delta\lambda}=mN\tag{6.17}

で、照らされた本数と次数の積で決まる。格子の間隔 dd そのものは分解能に入らない。

図

図 6.2 幅 bb のスリットを間隔 d=3bd=3b で N=5N=5 本並べたときの Fraunhofer 回折の強度(実線)と、1本のスリットの像 (sin⁡u/u)2(\sin u/u)^2(破線)。主な山は dsin⁡θ=mλd\sin\theta=m\lambda に立ち、山の間に N−2=3N-2=3 個の小さな山がある。33 次の山は1本のスリットの像の 00(bsin⁡θ=λb\sin\theta=\lambda)と重なって消える

例題 6.3 回折格子でナトリウムの D 線を分ける ★★★

オリジナル1 mm あたり 600 本の回折格子に、ナトリウムの D 線(589.0 nm589.0\,\mathrm{nm} と 589.6 nm589.6\,\mathrm{nm}。平均 589.3 nm589.3\,\mathrm{nm})を垂直に入れる。

  1. (1) 1次の回折の角度を求めよ。何次までの光が出るか。
  2. (2) 格子を幅 25 mm25\,\mathrm{mm} にわたって照らすとき、1次と2次の分解能を求めよ。D 線を分けるのに要る分解能と、1次で分けられる最小の波長差を求めよ。
  3. (3) 照らす幅を 1 mm1\,\mathrm{mm} に絞ると、1次と2次で D 線は分かれるか。
解答

(1) d=1/600 mm=1.667 μmd=1/600\,\mathrm{mm}=1.667\,\mu\mathrm m。sin⁡θ=λ/d=0.3536\sin\theta=\lambda/d=0.3536 で θ=20.71∘\theta=20.71^\circ。∣sin⁡θ∣≤1\lvert \sin\theta\rvert \le1 なので次数は ∣m∣≤d/λ=2.83\lvert m\rvert \le d/\lambda=2.83、つまり2次までである。3次の「モード」は横の波数が kk を超え、6.1節のエバネッセント波になって格子の面に張り付いている。

(2) 照らされた本数は N=25×600=15000N=25\times600=15000 で、式 (6.17) から1次で λ/Δλ=15000\lambda/\Delta\lambda=15000、2次で 3000030000。D 線を分けるには 589.3/0.6=982589.3/0.6=982 あればよいので、どちらでも十分に分かれる。1次で分けられる最小の波長差は 589.3/15000=0.0393 nm589.3/15000=0.0393\,\mathrm{nm} である。

(3) N=600N=600 で、1次の分解能 600600 は 982982 に足りず分かれない。2次では 1200>9821200>982 で分かれる。照らす幅を狭めても格子の式 (6.16) の山の位置は変わらないが、山の幅が 1/N1/N に比例して太くなり、隣の線と重なる。

ほかの巻では

木下の勉強ステータス

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

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

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