光学と量子エレクトロニクス 第9章 §9.4 ★★★ 広がった光源からの光の2点の間の相関は、光源の強度分布の Fourier 変換で決まる。式の形は Fraunhofer 回折と同じで、光源が小さいほど遠くの2点まで相関する
2点の相関と Young の縞
時間の相関が「同じ点の、時刻の違う2つの場」の相関だったのに対し、空間の相関は「同じ時刻の、場所の違う2つの場」の相関である。光の進む向きに垂直な面の上の2点 x 1 \bm{x}_1 x 1 , x 2 \bm{x}_2 x 2 の場の相関を
G 12 ( 1 ) = ⟨ E ∗ ( x 1 , t ) E ( x 2 , t ) ⟩ , g 12 ( 1 ) = G 12 ( 1 ) ⟨ ∣ E ( x 1 ) ∣ 2 ⟩ ⟨ ∣ E ( x 2 ) ∣ 2 ⟩ (9.16) G^{(1)}_{12}=\langle E^*(\bm{x}_1,t)E(\bm{x}_2,t) \rangle ,\qquad g^{(1)}_{12}=\frac{G^{(1)}_{12}}{\sqrt{\langle \lvert E(\bm{x}_1)\rvert ^2 \rangle \langle \lvert E(\bm{x}_2)\rvert ^2 \rangle }}\tag{9.16} G 12 ( 1 ) = ⟨ E ∗ ( x 1 , t ) E ( x 2 , t )⟩ , g 12 ( 1 ) = ⟨∣ E ( x 1 ) ∣ 2 ⟩ ⟨∣ E ( x 2 ) ∣ 2 ⟩ G 12 ( 1 ) ( 9.16 )
と書く。この面に2つの小さな穴をあけて後ろの幕に Young の縞を作ると、2つの穴からの光の強度が等しいとき、幕の中心付近の縞の可視度は ∣ g 12 ( 1 ) ∣ \lvert g^{(1)}_{12}\rvert ∣ g 12 ( 1 ) ∣ になる。計算は式 (9.12) と同じで、穴から幕の点までの2本の経路の差が、時間差 τ \tau τ の代わりに位相差を作る(経路の差がコヒーレンス長より十分短い範囲で)。
van Cittert–Zernike の定理
定理 9.3(van Cittert–Zernike の定理)
面 z = 0 z=0 z = 0 に広がった光源があり、次の3つを仮定する。(i) 光は準単色で、問題にする経路の差はコヒーレンス長より十分短い(波数を中心の値 k k k で代表させてよい)。(ii) 光源の異なる点は互いに無相関で、光源の面の上の場が ⟨ s ∗ ( ξ ) s ( ξ ′ ) ⟩ = I ( ξ ) δ ( 2 ) ( ξ − ξ ′ ) \langle s^*(\bm{\xi})s(\bm{\xi}') \rangle =I(\bm{\xi})\,\delta^{(2)}(\bm{\xi}-\bm{\xi}') ⟨ s ∗ ( ξ ) s ( ξ ′ )⟩ = I ( ξ ) δ ( 2 ) ( ξ − ξ ′ ) を満たす(I ( ξ ) I(\bm{\xi}) I ( ξ ) は光源の強度分布)。(iii) 光源の大きさと観測点の横の位置は距離 z z z より十分小さい(近軸)。このとき距離 z z z の面の上の2点の相関は
g 12 ( 1 ) = e i ψ ∫ I ( ξ ) e − i k ξ ⋅ ( x 2 − x 1 ) / z d 2 ξ ∫ I ( ξ ) d 2 ξ , ψ = k ( ∣ x 2 ∣ 2 − ∣ x 1 ∣ 2 ) 2 z (9.17) g^{(1)}_{12}=\mathrm{e}^{\mathrm{i}\psi}\,\frac{\displaystyle\int I(\bm{\xi})\,\mathrm{e}^{-\mathrm{i} k\bm{\xi}\cdot(\bm{x}_2-\bm{x}_1)/z}\,\mathrm{d}^2\xi}{\displaystyle\int I(\bm{\xi})\,\mathrm{d}^2\xi},\qquad
\psi=\frac{k\bigl(\lvert \bm{x}_2\rvert ^2-\lvert \bm{x}_1\rvert ^2\bigr)}{2z}\tag{9.17} g 12 ( 1 ) = e i ψ ∫ I ( ξ ) d 2 ξ ∫ I ( ξ ) e − i k ξ ⋅ ( x 2 − x 1 ) / z d 2 ξ , ψ = 2 z k ( ∣ x 2 ∣ 2 − ∣ x 1 ∣ 2 ) ( 9.17 )
となる。可視度 ∣ g 12 ( 1 ) ∣ \lvert g^{(1)}_{12}\rvert ∣ g 12 ( 1 ) ∣ は、光源の強度分布の2次元 Fourier 変換を、空間周波数 k ( x 2 − x 1 ) / z k(\bm{x}_2-\bm{x}_1)/z k ( x 2 − x 1 ) / z で読んだものである。
証明 光源の点 ξ \bm{\xi} ξ から観測点 x \bm{x} x へは球面波 e i k r / r \mathrm{e}^{\mathrm{i} kr}/r e i k r / r (r r r は2点の距離)が伝わるので、E ( x ) = ∫ s ( ξ ) e i k r ( x , ξ ) / r d 2 ξ E(\bm{x})=\int s(\bm{\xi})\mathrm{e}^{\mathrm{i} kr(\bm{x},\bm{\xi})}/r\,\mathrm{d}^2\xi E ( x ) = ∫ s ( ξ ) e i k r ( x , ξ ) / r d 2 ξ である(共通の定数は省く)。(ii) を使うと二重積分が一重になり、G 12 ( 1 ) = ∫ I ( ξ ) e i k ( r 2 − r 1 ) / ( r 1 r 2 ) d 2 ξ G^{(1)}_{12}=\int I(\bm{\xi})\mathrm{e}^{\mathrm{i} k(r_2-r_1)}/(r_1r_2)\,\mathrm{d}^2\xi G 12 ( 1 ) = ∫ I ( ξ ) e i k ( r 2 − r 1 ) / ( r 1 r 2 ) d 2 ξ 。(iii) から r j ≈ z + ∣ x j − ξ ∣ 2 / ( 2 z ) r_j\approx z+\lvert \bm{x}_j-\bm{\xi}\rvert ^2/(2z) r j ≈ z + ∣ x j − ξ ∣ 2 / ( 2 z ) で、
r 2 − r 1 ≈ ∣ x 2 ∣ 2 − ∣ x 1 ∣ 2 2 z − ξ ⋅ ( x 2 − x 1 ) z r_2-r_1\approx\frac{\lvert \bm{x}_2\rvert ^2-\lvert \bm{x}_1\rvert ^2}{2z}-\frac{\bm{\xi}\cdot(\bm{x}_2-\bm{x}_1)}{z} r 2 − r 1 ≈ 2 z ∣ x 2 ∣ 2 − ∣ x 1 ∣ 2 − z ξ ⋅ ( x 2 − x 1 )
となる(∣ ξ ∣ 2 / ( 2 z ) \lvert \bm{\xi}\rvert ^2/(2z) ∣ ξ ∣ 2 / ( 2 z ) は差をとると消える)。分母の r 1 r 2 r_1r_2 r 1 r 2 は z 2 z^2 z 2 で置き換えてよく、規格化の分母は ∫ I d 2 ξ / z 2 \int I\,\mathrm{d}^2\xi/z^2 ∫ I d 2 ξ / z 2 である。
式 (9.17) は、開口の透過関数を光源の強度分布に置き換えた Fraunhofer 回折の式(第6章 )と同じ形をしている。ただし Fraunhofer 回折では開口の上の場の振幅を Fourier 変換したのに対し、ここでは光源の強度を Fourier 変換している。光源の各点が無相関なので、振幅でなく強度が足し合わされるからである。また、光源の面の上の位相 k ∣ ξ ∣ 2 / ( 2 z ) k\lvert \bm{\xi}\rvert ^2/(2z) k ∣ ξ ∣ 2 / ( 2 z ) が差をとると消えるので、Fraunhofer 回折のような「光源が λ z \sqrt{\lambda z} λ z より小さい」という条件は要らない。
一様な明るさの円盤(半径 a a a )の光源では、円の Fourier 変換(第6章の Airy 円盤と同じ計算)から
∣ g 12 ( 1 ) ∣ = ∣ 2 J 1 ( v ) v ∣ , v = k θ s d = π θ D d λ (9.18) \lvert g^{(1)}_{12}\rvert =\Bigl\lvert\frac{2J_1(v)}{v}\Bigr\rvert,\qquad v=k\theta_{\mathrm s}d=\frac{\pi\theta_{\mathrm D}d}{\lambda}\tag{9.18} ∣ g 12 ( 1 ) ∣ = v 2 J 1 ( v ) , v = k θ s d = λ π θ D d ( 9.18 )
となる。d = ∣ x 2 − x 1 ∣ d=\lvert \bm{x}_2-\bm{x}_1\rvert d = ∣ x 2 − x 1 ∣ は2点の間隔、θ s = a / z \theta_{\mathrm s}=a/z θ s = a / z は光源の見込む角の半径、θ D = 2 θ s \theta_{\mathrm D}=2\theta_{\mathrm s} θ D = 2 θ s は角直径、J 1 J_1 J 1 は1次の Bessel 関数である。J 1 J_1 J 1 の最初の零点 v = 3.832 v=3.832 v = 3.832 で可視度が 0 0 0 になるので、相関が最初に消える間隔は
d 0 = 3.832 π λ θ D = 1.22 λ θ D d_0=\frac{3.832}{\pi}\frac{\lambda}{\theta_{\mathrm D}}=1.22\,\frac{\lambda}{\theta_{\mathrm D}} d 0 = π 3.832 θ D λ = 1.22 θ D λ
である(図 9.4 )。同じ 1.22 1.22 1.22 が Airy 円盤の半径にも出るのは、どちらも一様な円の Fourier 変換だからである。
図 9.4 一様な円盤の光源(角直径 θ D \theta_{\mathrm D} θ D )からの光の、間隔 d d d の2点の相関の大きさ(式 (9.18) )。d = 1.22 λ / θ D d=1.22\lambda/\theta_{\mathrm D} d = 1.22 λ / θ D で最初に 0 0 0 になり、その先は小さな山が続く。
例題 9.7 太陽の光の相関の幅 ★★
オリジナル 太陽を、角直径 θ D = 0.53 ∘ \theta_{\mathrm D}=0.53^\circ θ D = 0.5 3 ∘ の一様な円盤とみなす。波長 550 n m 550\,\mathrm{nm} 550 nm の光(緑のフィルタを通した光)について、
(1) 地上で相関が最初に消える2点の間隔 d 0 d_0 d 0 を求めよ。
(2) 太陽の光で Young の縞を見るには、2つの穴の間隔をどうすればよいか。可視度が 0.88 0.88 0.88 以上になる間隔を、式 (9.18) を v v v の2次まで展開して見積もれ。
解答 (1) θ D = 0.53 × π / 180 = 9.25 × 10 − 3 r a d \theta_{\mathrm D}=0.53\times\pi/180=9.25\times10^{-3}\,\mathrm{rad} θ D = 0.53 × π /180 = 9.25 × 1 0 − 3 rad で、d 0 = 1.22 × 550 × 10 − 9 / ( 9.25 × 10 − 3 ) = 72.5 μ m d_0=1.22\times550\times10^{-9}/(9.25\times10^{-3})=72.5\,\mu\mathrm m d 0 = 1.22 × 550 × 1 0 − 9 / ( 9.25 × 1 0 − 3 ) = 72.5 μ m である。
(2) J 1 ( v ) = v / 2 − v 3 / 16 + ⋯ J_1(v)=v/2-v^3/16+\cdots J 1 ( v ) = v /2 − v 3 /16 + ⋯ から 2 J 1 ( v ) / v ≈ 1 − v 2 / 8 2J_1(v)/v\approx1-v^2/8 2 J 1 ( v ) / v ≈ 1 − v 2 /8 。1 − v 2 / 8 ≥ 0.88 1-v^2/8\ge0.88 1 − v 2 /8 ≥ 0.88 は v ≤ 0.98 v\le0.98 v ≤ 0.98 で、d ≤ 0.98 λ / ( π θ D ) = 0.31 λ / θ D = 18.5 μ m d\le0.98\lambda/(\pi\theta_{\mathrm D})=0.31\lambda/\theta_{\mathrm D}=18.5\,\mu\mathrm m d ≤ 0.98 λ / ( π θ D ) = 0.31 λ / θ D = 18.5 μ m である(v = 0.98 v=0.98 v = 0.98 での厳密な値は 0.884 0.884 0.884 )。太陽の光をそのまま使うなら、2つの穴を 20 μ m 20\,\mu\mathrm m 20 μ m 程度まで近づける必要がある。穴の間隔を広げたいなら、太陽の光をいったん小さな穴(角直径を小さくした光源)に通してから使う。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com