この記事は,東京大学航空宇宙工学科/専攻 Advent Calendar 2017向けの記事です.
今年の4~6月あたりに得た波動光学の基礎をまとめておく.
波動光学,和書があまりなく,毎回毎回洋書を参照するのもだるいので.
トップ画像の出典は[2].
「【幾何光学】カメラのレンズ性能の式などの導出」などの記事で,レンズの基本的な物理量についてみてきた.
F値や被写界深度は幾何光学,つまり光が進む経路(反射や屈折など)を考えるだけで十分だった.
しかしながら,MTFなどは幾何光学では説明不可能である.
([1]のようにレンズメーカーのホームページにはMTF特性図が載っている.)
そこで,ここでは光を波動として考える波動光学の基礎について,備忘録的にまとめたい.
(備忘録的,というのは,波動光学の考え方を卒業研究で多用し修士研究でも用いる予定なのだが,現在は卒業設計に時間をとられ,せっかく学んだ波動光学の基礎を忘れそう,という事情による.)
波動光学とは,光を波動として考える光学であり,干渉・回折・偏光などが扱える.
中学,高校で学ぶレンズの公式などは幾何光学であり,ヤングの実験などが波動光学である.
さらに発展として,量子力学も考慮した量子光学という分野も存在する.
はじめに「2.波動光学の基本的な回折の定式化」で波動光学の基礎式となる回折現象を定式化する.
次に「3.像形成」で薄肉レンズでの像をそのインパルス応答を定式化することによって理解する.
さらに「4.光学系のシステム応答解析」で一般の結像系について,そのシステム応答を考える.
そして続編の記事「【波動光学】波動光学入門2 - 適当な数値計算例」では,波面収差の取り扱いと実際の数値計算を行っている.
最初は式変形ばかりでわかりにくいが,続編で可視化したものを見れば雰囲気はわかるだろう.
幾何光学では,点光源をレンズを介して再び完全な点像へと集光できた.
しかし,実際にはレンズを通ったときの回折により一点へと集光させることは不可能である.
そこで,まずはじめに回折を簡単に記述するための波動光学の基礎近似式を導出する.
ここでは簡単のために完全な単色光を仮定する.
位置rでの正弦波複素振幅U(r,t)は,波数ベクトルkを用いて,
U(r,t)=Aexp[i(k⋅r−ωt)]となるのは流石にいいよね?(卒業設計後にこれを読み返すであろう自分を煽っておく.)
(まあ,ここではもっぱら時間平均しか考えないので,時間項は落ちるが.)
上の座標系を用いて,ホイヘンス - フレネル原理に基づく回折式を導入する.
ホイヘンス - フレネル原理とは,高校の教科書ではホイヘンスの原理と呼ばれているあれで,
“前進波の波面の各点が二次波とよばれる新しい波の波源となり,全体としての前進波は全ての二次波を重ね合わせたものとなる.”
というもの.
このとき,点P0での複素振幅U(P0)は,点P0,P1間距離r01を用いて
U(P0)=1iλ∬ΣU(P1)eikr01r01cosθdsとなる.さらに
cosθ=zr01より,ホイヘンス - フレネル原理に基づく回折式は
U(x,y)=ziλ∬ΣU(ξ,η)eikr01r201dξdηとなる.
...と,まあいきなり(2-2)式を持ち出して変形したわけだが,一般的な光学入門書だと,(2-4)式でz≈r01より
U(x,y)=∬ΣU(ξ,η)eikriλrdξdηから議論がスタートする.スケール因子1/iλの説明もなしで.
波動光学を初めて勉強したときはこの(2-5)式をいきなり出されて??となっていたが,定性的には,
U(ξ,η)exp[ikr]/rがホイヘンス - フレネル原理の二次波の球面波であり,それを回折開口で積分している,と理解できる.
(球面波とは,点波源より伝播する波面が球となる波.点波源(波面の中心)からの距離rを用いるとその振幅はexp[ikr]/rに比例.)
(2-2)式の導出はとても煩雑であるので,[4]あたりを参照されたいが,概要は,波動方程式を開口によって決まる境界条件(キルヒホッフ境界条件)の下で解けば求まる.
(2-4)式を近似していく.
r01=√z2+(x−ξ)2+(y−η)2≈z[1+12(x−ξz)2+12(y−ηz)2]となるが,(2-4)式において,分母のr01はこの二次の項を落とすことによる影響は小さいが,exp内のr01でこれを落とすのは影響が無視できない.
したがって,(2-4)式は
U(x,y)=eikziλz∞∬−∞U(ξ,η)eik2z[(x−ξ)2+(y−η)2]dξdηとなり,畳み込みの形で表される.これがフレネル回折積分と呼ばれるものである.
ここで積分範囲が無限だが,開口が有限範囲であるという条件はU(ξ,η)に含まれている.
さらにξ,ηがかかわらない項を外に出して,
U(x,y)=eikziλzeik2z(x2+y2)∞∬−∞{U(ξ,η)eik2z(ξ2+η2)}e−j2πλz(xξ+yη)dξdη=eikziλzeik2z(x2+y2)F{U(ξ,η)eik2z(ξ2+η2)}fξ=xλz,fη=yλzとフーリエ変換F{}で表すことができる.
さらに厳しい近似(フラウンホーファー近似)
z≫k(ξ2+η2)max2を(2-8)式にほどこすと,
U(x,y)=eikzeik2z(x2+y2)iλz∞∬−∞{U(ξ,η)}e−j2πλz(xξ+yη)dξdη=eikzeik2z(x2+y2)iλzF{U(ξ,η)}fξ=xλz,fη=yλzとなり,これをフラウンホーファー回折積分といい,zが大きい,つまり無限遠で観測させる場合に有効である.
今までは,有限の開口を光が通過することによって生じる回折を定量的に記述した.
次は下のような座標系で像形成を考えよう.
なお,以下では定数係数は省略する.
一般論として,線形時不変なシステム応答を考えた場合,出力は単位インパルス応答の畳み込みで表示される.
これを光学系へあてはめると,単位インパルス応答は単位振幅の点光源像が相当し,それをhとする.
するとある物体Uoが作る像Uiは
Ui(u,v)=h(u,v)∗Uo(u,v)=∞∬−∞h(u−ξ,v−η)Uo(ξ,η)dξdηとなる.
後の式変形からわかるが,実際に求めると
Ui(u,v)=∞∬−∞h(u−Mξ,v−Mη)Uo(ξ,η)dξdηとスケール因子がかかってしまうので,ここではh(u−ξ,v−η)を一般的にh(u,v;ξ,η)と書くことにする.
光学系に限った話ではなく,線形時不変システムに一般的な話だが,結像系の特性はインパルス応答hを特定することによって完全に記述できるので,このhを求める.
インパルス応答は点光源像となるため,物体を点(ξ,η)においたデルタ関数とする.
この場合,レンズに入射するのは点(ξ,η)を中心とする発散球面波であり,球面波は波源からの距離rを用いて
eikrrとなるが,フレネル近似と同様に近軸近似をほどこし
Ul(x,y)=1iλz1eik2z1[(x−ξ)2+(y−η)2]となる((2-7)式の開口を無限小点にした場合と同値).
また,焦点距離fのレンズを通過した後の光波場分布は,導出は省略するが,
U′l(x,y)=Ul(x,y)P(x,y)e−ik2f(x2+y2)となる(凸レンズは発散球面波を収束球面波へ変換すると考えれば自明).
なお,Pとは瞳関数と呼ばれるものであり,レンズ開口内で1を,開口外で0をとる関数であり,レンズの形状を規定する.
最後に,距離z2の伝播をフレネル回折の式(2-7)を用いると,
h(u,v;ξ,η)=1iλz2∞∬−∞U′l(x,y)eik2z2[(u−x)2+(v−y)2]dxdyが得られる.
上記3式を組み合わせると,
h(u,v;ξ,η)=1λ2z1z2eik2z2(u2+v2)eik2z1(ξ2+η2)×∞∬−∞P(x,y)eik2(1z1+1z2−1f)(x2+y2)×e−jk[(ξz1+uz2)x+(ηz1+vz2)y]dxdyとなるが,とても煩雑である.
詳細な議論は[4]などに任せるが,結像面をレンズの公式をみたす距離に配置すれば
1z1+1z2−1f=0とでき,
1. | 像面の強度分布のみに注目する場合で,雑煮関連する位相分布は重要ではない. |
2. | 像の光波場分布に注目するが,像の観察は,光軸が薄肉レンズを貫く点を中心とする半径z2の球面上で行う. |
のどちらかをみたす場合,
eik2z2(u2+v2)≈1の近似が成立し,さらに
1. | 光軸が薄肉レンズを貫く点を中心とする半径z1の球面上に物体がある. |
2. | 光軸がレンズを貫く点に向かって収束する球面波で物体が照射される. |
3. | この2次位相因子の位相が,特定の像点(u,v)の光波場に大きく寄与する物体部分の内側ではほとんど変化しない. |
のいずれかをみたす場合,
eik2z1(ξ2+η2)≈1の近似が成立する.
さらにシステムの倍率
M=−z2z1を定義すると,(3-7)式は
h(u,v;ξ,η)≈1λ2z1z2∞∬−∞P(x,y)e−i2πλz2[(u−Mξ)x+(v−Mη)y]dxdyとなる.
したがって,レンズの公式がみたされる場合,インパルス応答は(スケール因子1/λz1も含め),像座標(u=Mξ,v=Mη)を中心に位置するレンズ開口のフラウンホーファー回折像と等しくなっている.
これは少し考えてみればそれほど意外な結果ではない.
つまりレンズの公式をみたすようにz2を設定したので,像平面はレンズ通過後の球面波が収束する点であり,この収束点近傍の光波場分布は,球面波の拡がりを制限するレンズ開口のまさにフラウンホーファー回折像である.
さらに結像系を一般化し,結像系を線形時不変システムとして描写していく.
これ以降の議論は和書ではほとんど触れられておらず,仕方がないので[2][5]などで学んだ.
ここではまず光学系の一般化モデルを定義し,次に回折限界(幾何光学的には点光源(発散球面波)が点像(収束球面波)へと移されるシステム.波動光学的には回折の影響により,完全な点像とはならない.)における結像システムを定式化する.
さらにその後にその周波数特性についてみていく.
再びCanonのレンズを例に取るが,下図の用に光学系は多数のレンズ,鏡群によって構成されている.
しかし物体上の点からは球面波が放射され,そして結像系を経て(球面波となり)最終的に像として結像すると仮定すれば,結像系を一つのブラックボックスにまとめてしまっても,そのシステム特性を記述する上では問題ない.
ここでは,下図のように物体,入射瞳,射出瞳,像を定義し,また座標系も下図のように設定する.
空間的にコヒーレントな光学系の場合,像の複素振幅Uiは物体の複素振幅Uo と点像振幅分布関数 (Amplitude Spread Function / ASF) hの畳み込み
Ui(u,v)=∞∬−∞h(u−ξ,v−η)Uo(ξ,η)dξdη=h(u,v)∗Uo(u,v)で表される.ASFは,「3.像形成」での議論により,理想像(幾何光学的に求まる無回折時の像)点付近の光振幅は射出瞳のフラウンホーファー回折像であるため,
h(u,v;ξ,η)=Aλzi∞∬−∞P(x,y)e−i2πλzi[(u−Mξ)x+(v−Mη)y]dxdyとなった.
ただしziは射出瞳―像平面間距離,Pは瞳関数(開口内側で1,外側で0)である.
ここで
ξ=Mξ,η=Mηと座標を取り直すとASFは,
h(u−ξ,v−η)=Aλzi∞∬−∞P(x,y)e−i2πλzi[(u−ξ)x+(v−η)y]dxdyすなわち
h(u,v)=Aλzi∞∬−∞P(x,y)e−i2πλzi(ux+vy)dxdy=AλziF{P(x,y)}fx=uλzi,fy=vλziと瞳関数のフーリエ変換で表される.
以上より,一般的な回折限界の光学系の場合の像は,幾何光学的像と,射出瞳のフラウンホーファー回折像であるインパルス応答との畳み込みであるとみなせる.
つまり,コヒーレント結像系では,複素振幅に対して線形である.
空間的にインコヒーレントな光学系の場合,観測されるのは像の複素振幅Uiではなく強度Iiである.
像の強度Iiは物体の強度Ioと点拡がり関数 (Point Spread Function / PSF) sの畳み込み
Ii(u,v)=∞∬−∞s(u−ξ,v−η)Io(ξ,η)dξdη=s(u,v)∗Io(u,v)で表される.
Iは時間平均演算子⟨⟩を用いて
I(u,v)=⟨|U(u,v,t)|2⟩=⟨U(u,v,t)U∗(u,v,t)⟩となることより,PSFはASFの絶対値の二乗
s(u,v)=|h(u,v)|2=|AλziF{P(x,y)}fx=uλzi,fy=vλzi|2となる.
以上より,インコヒーレント結像系は強度に関して線形であり,そのような結像系のインパルス応答は,複素インパルス応答の大きさの二乗で表される.
(4-1)式のフーリエ変換をとると,畳み込みの公式より,
F{Ui(u,v)}=F{h(u,v)}F{Uo(u,v)}となり,
H(fx,fy)=F{h(u,v)}は振幅伝達関数 (amplitude transfer function / ATF) と呼ばれ,システムの周波数応答を表す.
ATFは
H(fx,fy)=F{h(u,v)}=F{AλziF{P(x,y)}fx=uλzi,fy=vλzi}=AλziP(−λzifx,−λzify)となる.
したがって,コヒーレント結像系ではその周波数特性は瞳関数と一致し,その通過周波数では振幅や位相を歪めることなく通過させる.
インコヒーレント結像系でもコヒーレント結像系の場合と同様にシステムの周波数応答Hを(0周波数の値で正規化して)
H(fx,fy)=F{s(u,v)}∞∬−∞s(u,v)dudv=F{|h(u,v)|2}∞∬−∞|h(u,v)|2dudvと求めることができる.
Hは光伝達関数 (optical transfer function / OTF) といい,その絶対値|H|を変調伝達関数 (modulation transfer function / MTF) という.
なお,OTFはパーセバルの定理を用いて
H(fx,fy)=F{|h(u,v)|2}∞∬−∞|h(u,v)|2dudv=∞∬−∞H∗(p−fx,q−fy)H(p,q)dpdq∞∬−∞|H(p,q)|2dpdq=∞∬−∞P∗(x−λzifx,y−λzify)P(x,y)dxdy∞∬−∞|P(x,y)|2dxdqとなり,瞳関数の自己相関関数であることがわかる.
「【波動光学】波動光学入門2 - 適当な数値計算例」で,この記事のつづきとして,
・ | 波面収差が存在する場合 |
・ | 具体的な数値計算結果 |
を掲載する.
[1] | Canon. EF24-70mm F4L IS USM 仕様. Retrieved December 24, 2017, from http://cweb.canon.jp/ef/lineup/standard-zoom/ef24-70-f4l/spec.html |
[2] | Jason D. Schmidt. Numerical Simulation of Optical Wave Propagation with example in MATLAB®. SPIE PRESS, 3 edition, 2013. |
[3] | Joseph W. Goodman. Introduction to Fourier Optics. W. H. Freeman, 3 edition, 2004. |
[4] | Joseph W. Goodman. Trans. by 尾崎 義治, 朝倉 利光. フーリエ光学. 森北出版株式会社, 第3版, 2012. |
[5] | James B. Breckinridge. Basic Optics for the ASTRONOMICAL SCIENCES. SPIE PRESS, 2012. |
1. | 【幾何光学】カメラのレンズ性能の式などの導出 [2017-08-15] |
2. | 【波動光学】波動光学入門2 - 適当な数値計算例 [2017-12-24] |
おにぎり [2019/12/11 14:26]
分かりやすい説明ですごくためになりました。
自分は4f系の顕微鏡(焦点距離の違う2枚レンズで観察対象をフーリエ変換して逆フーリエ変換する)をシミュレーションをしたい思っているのですが、まず最初にコヒーレントな場を考えようとしています。そこで質問なのですが、レンズでフーリエ変換をしたいとき(z1=z2=f)インパルス応答などを考える必要があるんでしょうか?
(3-8)式が成り立たないのでこの場合不適切なのでしょうか?
溶けかけてるうさぎ(管理者) [2019/12/14 00:58]
おにぎりさん.
光学の専門家ではないので,はっきりしたことは言えませんが,レンズが何枚だろうと基本的には1つの光学系(システム)として取り扱えるはずなので,結像させるのであれば(3-8)式が成立しないことはないと思います.
そして,(線形時不変を仮定できる場合は)インパルス応答が求まっていれば,どんな入力に対しても(任意の形状の入力に対しても)その出力を求められます.
おにぎり [2019/12/14 01:27]
返信ありがとうございます。
勉強になります!
名前
Email (※公開されることはありません)
コメント