エッセンシャルテキスト 光学(第1~第3章)

左貝潤一著「エッセンシャルテキスト 光学」森北出版社(2019.7)
www.amazon.co.jp/dp/4627776314

読後メモ(行間の補充等)

第1章 光学の基礎事項

p1 波数と角波数

波数κ、角波数kκ=1λ        [1m]k=2πλ        [radm]

p8 光線の単位ベクトル

位置ベクトルをr=(x,y,z)、点Aの位置ベクトルをrA、点Bの位置ベクトルをrB、この位置ベクトルの変化をdrとすると、rA=(xA,yA,zA)rB=(xB,yB,zB) よって、dr=rBrA=(xBxA,yByA,zBzA) この点AからBへの変化が光線の経路上で起きる場合の光線に沿った点Aと点Bとの間の幾何学的距離をds、点Aにおける光線の単位ベクトル(光線の方向を示す大きさ1のベクトル)をsとすると、s ds=drよって、s=drdsdrの方向余弦成分(x,y,z要素)を用いれば、s=(dxds,dyds,dzds)=(xBxAds,yByAds,zBzAds)
sは単位ベクトルであるから、方向余弦成分として、α2+β2+γ2=1となる(α,β,γ)を用いて、s=(α,β,γ)と現わすこともできる。


p8 アイコナール

アイコナールφ
光路長φは、空間上の位置によって値が異なるスカラー値であるから、位置ベクトルをrとすると、φ(r)と表現できる。
いま、同じ光路長Cφ(ra)=Cφ(rb)=Cφ(rc)=Cを持つ空間上の各点ra=(xa,ya,za)rb=(xb,yb,zb)rc=(xc,yc,zc)(x,y,z)グラフ上に集めると三次元上の曲面となる。これは、特定の光路長(光源からの光路長が同じ場合、波動の位相値も等しい)φに対応した電磁波の波面Cであるが、他の特定の光路長φの集団についても、同様に対応する波面Cを考えることができるので、結局、光路長φ(r)の等位面としてφ(r)=C1φ(r)=C2φ(r)=C3対応する各波面C1,C2,C3,を(x,y,z)座標で表すことができる。

p8 波面の法線

波面の法線grad φ
空間上の位置によって値が異なるスカラー値たる光路長φについて、空間内の変化率(勾配)grad φ=(φx,φy,φz)と波面Cとの関係を考える。
波面Cに微小距離だけ離れた二つの点を考えると、両点を結ぶ直線ベクトルΔr=(dx,dy,dz)は、波面Cの接線となる。両ベクトルの内積を計算すると、grad φΔr=φxdx+φydy+φzdz=Δφとなる。Δφは、具体的には、両点間におけるφの変化量を意味するが、両点は、いずれも波面C上の点なのであるから、φの値に変動はない。よって、Δφの値は、0である。従って、この内積の値も0となるので、grad φΔr=0|grad φ||Δr|cosπ2=0と変形できる。すなわち、grad φは、波面Cの接線Δrと垂直の関係にある。そして、波面Cの接線Δrは、波面Cの接面上で360度どの方向にも想定できるから、結局、grad φは、波面Cに対して垂直な方向を示すベクトルであると言える。
なお、同ベクトルを同ベクトルの絶対値で除すると、大きさが1のベクトル(単位ベクトル)となるところ、これは光線の方向(波面Cに垂直)を示す大きさ1のベクトルと等しいからs=grad φ|grad φ|となる。

p8 アイコナール方程式

アイコナール方程式(grad φ)2=n(r)2
光線上の点A及び点Bとの間の光路長(光学的距離)の差をφBAとし、点Aと光線上にない点C(点Bから微小距離離れた点)との光路長の差をφCAとした場合の両者の差φCAφBAは、φCBに等しい。φCAφBA=φCB他方で、点Bと点Cとの距離差は微小であること(=点Bや点Cが、点Aからは相対的に十分に遠い距離に位置すること)を考えると、線分ABと線分ACとは平行であり、また、点Bの属する波面と点Cの属する波面とは平行であると近似できる。この場合、光線ABが属する直線と点Cが属する波面の交点をDとすれば、点B、C、Dの組合せからなる各線分は、幾何学的に、以下の関係にある(θは両線分のなす角度)。BC cosθ=BDこれを光路長に置き換えると、屈折率をn として、nBC cosθ=nBDベクトルとその内積で左辺を表記すると、nBCs=nBD右辺を光路長φの差分として表現すると、nBCs=Δφいま、点Bと点Cとの距離差は微小であることを前提とするので、BCdrΔφdφで表現しなおすと、n drs=dφベクトル部分について、成分表示dr=(dxdydz) ds=(αβγ) に置き換えたうえで、内積を演算すると、dφ=n(dxdydz)(αβγ)=n(α dx+β dy+γ dz)が得られる。
この式を、各々、x,y,zにて偏微分すると、φx=nαφ y=nβφz=nγとなるので、各々を二乗して合計した値は、sが単位ベクトルであることを踏まえれば(φx)2+(φy)2+(φz)2=n2α2+n2β2+n2γ2=n2(α2+β2+γ2)=n2|s|=n2左辺を勾配(grad)で表現し、右辺の屈折率nが位置ベクトルrの関数であることを踏まえて表現すると、(grad φ)2=n(r)2が得られる。

p8 光線方程式

s=grad φn(r)drds=grad φn(r)n(r)drds=grad φ両辺をsで微分すると、dds[n(r)drds]=ddsgrad φ=grad ddsφ=grad ddsn(r)ds=grad n(r)

p11 合成波

進行波1(最大振幅(電界)A、角速度ω1、角波数k1)の時刻t、位置zの瞬時振幅(電界)u1u1=Acos(ω1tk1z)進行波2(最大振幅(電界)A、角速度ω2、角波数k2)の時刻t、位置zの瞬時振幅(電界)u2u2=Acos(ω2tk2z)の合成波Uの瞬時振幅は、U=u1+u2=Acos(ω1tk1z)+Acos(ω2tk2z)となるところ、三角関数の公式cosα+cosβ=2cosα+β2cosαβ2を用いると、U=2Acosω1tk1z+ω2tk2z2cosω1tk1zω2t+k2z2=2Acos(ω1+ω2)t(k1+k2)z2cos(ω1ω2)t(k1k2)z2=2Acos(ω1ω2)t(k1k2)z2cos(ω1+ω2)t(k1+k2)z2=2Acos(ω1ω22tk1k22z)cos(ω1+ω22tk1+k22z)となり、波長の異なる二種類の波の積の形に変形できる。
波1の波長をλ1、波2の波長をλ2とすると、λ1=2πk1λ2=2πk2であり、合成波を構成する二種類の波のうち波長の長い方(包絡線)の波長λgは、λg=2πkg=2π(k1k22)=2×2π(2πλ12πλ2)=2(1λ11λ2)=2λ1λ2λ2λ1となる(p11図1.6の包絡線の矢印間は半波長)。

p12 群速度

群速度vg
一般に波速度v、波長λ、振動数(周波数)ν、角波数k、角速度ωとの間には、v=λν=λ2π2πν=1kω=ωkの関係があるところ、p11の式(1.26)U=2Acos(ω1ω2)t(k1k2)z2cos(ω1+ω2)t(k1+k2)z2のうち、cos(ω1ω2)t(k1k2)z2の部分が包絡線の要素であるから、同部分の速度(群速度)は、vg=ωgkg=(ω1ω22)(k1k22)=ω1ω2k1k2となる。ここで、平均値として、ω=ω1+ω22k=k1+k22を考えると、ロピタルの定理limxaf(x)f(a)g(x)g(a)=limxaf(x)g(x)より、ω1ω2k1k2が微小である場合には、ω1ω2k1k2dωdkであるから、vgは、vgdωdkとなる。
これを前提に、周波数ωによって屈折率nが変化する媒質中における群速度を検討すると、vg=1dkdω=1ddω(n(ω)ωc)=cddω(n(ω)ω)=cddω(n(ω))×ω+n(ω)ddωω=cωddω(n(ω))+n(ω)となる(式中のωは平均値)。

p12 位相速度

位相速度vp
式(1.26)U=2Acos(ω1ω2)t(k1k2)z2cos(ω1+ω2)t(k1+k2)z2のうち、cos(ω1+ω2)t(k1+k2)z2の部分が包絡線内部の周期的変化を示す要素であるから、同部分の速度(位相速度)は、vp=ωpkp=ω1+ω22k1+k22=ω1+ω2k1+k2=wkとなる。

p12 群速度と位相速度の関係

群速度vgと位相速度vpの関係
vg=dωdkvp=ωkから、vg=ddk(vp×k)=vpdkdk+dvpdkk=vp+kdvpdkであり、更に変形をすると、vg=vp+kdvpdλ×dλdk=vp+kdvpdλ×ddk(2πk)=vp+kdvpdλ×(2πk2)=vpdvpdλ×2πk=vpdvpdλ×λ=vpλdvpdλ

第2章 屈折と反射


p16 反射の法則

(誤)AE=AB sinθr
(正)AE=AB sin(πθr)

p17 多層被膜の屈折角

光線が、屈折率n0の空中から、入射角θiで、多数被膜を通過する場合における最終層の屈折角θtを求める。第1層の屈折率をn1、第2層の屈折率をn2、…第m1層の屈折率をnm1、第m層(最終層)の屈折率をnmとした場合、被膜層厚さが十分に小さい場合、各層は平行平面とみなせるので、各層での屈折角は次層への入射角に等しくなる。これらを順に、θ1, θ2, , θm1として各層境界面でスネルの法則を適用すると、{n0sinθi=n1sinθ1n1sinθ1=n2sinθ2nm1sinθm1=nmsinθtとなるので、結局、n0sinθi=nm sinθtと整理される。
よって、最初の層での入射角θi、最初の層の屈折率n0、最終層の屈折率nmのみで、最終層からの屈折角θtは求められる。

p19 マリュスの定理(ラグランジュの積分不変量)

n(r)を位置ベクトルrで定まる屈折率、光線の単位ベクトルをsとすると、
s=grad φn(r)n(r)s=grad φrot[n(r)s]=rot grad φrot[n(r)s]=0ここで、閉曲線Cで囲まれた曲面Sについてベクトルの回転(rot)を面積分したものと、閉曲線Cに従ってベクトルを線積分したものとが等しいというストークスの定理Srot AdS=CAdrをベクトルであるn(r)sに用いると、Srot [n(r)s]dS=Cn(r)sdr0=Cn(r)sdr

p21 プリズムの偏角

プリズムの偏角\(\delta\)
入射光線と出射光線との交点をDとすると、P21図の三角形DBCにおいて、B=θ1θ1C=θ2θ2であるから、δ=B+C=(θ1θ1)+(θ2θ2)=θ1+θ2(θ1+θ2)となり、三角形ABCにおいて、、B=π2θ1C=π2θ2であるから、α=A=π(B+C)=θ1+θ2となるので、偏角δは、δ=θ1+θ2αとまとめられる。

p21 最小偏角

最小偏角δm

スネルの法則1sinθ1=nsinθ1nsinθ2=1sinθ2即ちθ1=arcsin(n sinθ1)θ2=arcsin(n sinθ2)δ式に代入するとδ=arcsin(n sinθ1)arcsin(n sinθ2)α=arcsin(n sinθ1)arcsin(n sin(αθ1))αとなる。最小偏角δmは、この式の極値が候補となるから、ddθ1δ=0となるような条件(θ1の値)を求めてみる。0=ddθ1δ=ddθ1[arcsin(n sinθ1)arcsin(n sin(αθ1))α]=ddθ1[arcsin(n sinθ1)]ddθ1[arcsin(n sin(αθ1))]ここで、arcsinの微分公式ddxarcsinx=11x2を用いると、与式のうちddθ1[arcsin(n sinθ1)]の部分は、ddθ1[arcsin(n sinθ1)]=ddθ1(n sinθ1)×dd(n sinθ1)[arcsin(n sinθ1)]=n cosθ1×11(n sinθ1)2=n cosθ11(n sinθ1)2また、与式のうちddθ1[arcsin(n sin(αθ1))]の部分は、d[arcsin(n sin(αθ1))]dθ1=d(n sin(αθ1))dθ1×d[arcsin(n sin(αθ1))]d(n sin(αθ1))=d(αθ1)dθ1×d(n sin(αθ1))d(αθ1)×11(n sin(αθ1))2=1×n cos(αθ1)×11(n sin(αθ1))2=n cos(αθ1)1(n sin(αθ1))2となるので、結局、与式は、0=ddθ1δ=n cosθ11(n sinθ1)2n cos(αθ1)1(n sin(αθ1))2n cos(αθ1)1(n sin(αθ1))2=n cosθ11(n sinθ1)2となる。この方程式の分母どおし、分子どおしを比較すると、{cos (αθ1)=cos θ1sin (αθ1)=sin θ1であれば、方程式が満たされる。よって、αθ1=θ1α=2θ1θ1=α2がこの方程式の解となる(δが極値をとるための条件)。このとき、θ2=αθ1=αα2=α2であり、スネルの法則より、θ1=arcsin(n sinθ1)=arcsin(n sinα2)θ2=arcsin(n sinθ2)=arcsin(n sinα2)となり、θ1=θ2が成立し、最小偏角δmは、δm=θ1+θ2α=arcsin(n sinα2)+arcsin(n sinα2)α=2arcsin(n sinα2)αとなる(屈折率nが既知のとき頂角αのプリズムで最小偏角δmを求める式)。
最小偏角δmで、入射角を表現すれば、δm=θ1+θ2αδm=2 θ1αδm+α=2 θ1θ1=δm+α2となる。
逆に、頂角αのプリズムを入射光線に対して回転させて最小偏角δmを測定すれば、以下のとおり、プリズムの材質の屈折率nを求めることができる(最小偏角法)。1sinθ1=n sinθ1n=sinθ1sinθ1=sinδm+α2sinα2

p22 入射角等が微小の場合の近似

ラジアン角θに関して、limnsinθθ=1が成立することから、微小角θについて、sinθθと近似できる。この場合、空気から媒質(屈折率n)への光線入射にかかるスネルの法則は、入射角をθ、屈折角をθとして、1sinθ1=n sinθ1θ1n θ1と近似できる。出射角も同様にθ2n θ2である。この場合(入射角θ1及び出射角θ2が微小の場合。なお、α=θ1+θ2なので、αも微小角である)、空気中のプリズムの偏角δは、δ=(θ1θ1)+(θ2θ2)(nθ1θ1)+(nθ2θ2)=(n1)θ1+(n1)θ2=(n1)(θ1+θ2)=(n1)αと近似できる。

p23 ガウスのレンズ公式

物体位置s1、後側焦点距離f、像の位置s2として、s2s1=fs2fs2s1=1s2f1s1=1s21f1f=1s1+1s2

第3章 屈折・反射を用いた結像作用

p29 単一球面での屈折(式3.14)

p29図3.3において、三角形の内角と外角の関係から、以下の二式が成立する。ζ2+θt=AOPAOP+ζ1=θiよって、ζ2+θt+ζ1=θiζ2=θiθtζ1

p29 正弦定理

一般に三角形(頂点A,B,C。各頂点の対辺a,b,c)において、外接円の半径をRとすると、以下の等式が成り立つ(正弦定理)。asinA=bsinB=csinC=2Rこれを、p29図3.3の三角形APOに適用すると、AOsinP=PAsinO=POsinA対応する値を入れると、Rsinζ1=η1sin(ζ2+θt)=s1+Rsin(πθi)から、Rsinζ1=η1sin(θiζ1)=s1+Rsinθiを得る。左辺と右辺から、R sinθi=(s1+R)sinζ1sinθi=s1+RRsinζ1ここで、近軸光線近似として、sinθiθisinζ1ζ1s1η1を用いると、与式は、θiη1+RRζ1

p36 行列法での基本式

球面上の光線入射点での接平面におけるスネルの法則から、njsinθj=nj+1sinθj+1njsin(ζj+θ0)=nj+1sin(ζj+1+θ0)ここで、ζjζj+1は微小角なので、sinζjζjかつsinζj+1ζj+1である。また、sinθ0=sinxjRj1であるから、θ0も微小角なので、sinθ0θ0これらの合計であるζj+θ0ζj+1+θ0も微小角と言えるから、sin(ζj+θ0)ζj+θ0かつsin(ζj+1+θ0)ζj+1+θ0となるので、冒頭のスネルの法則の適用結果は、以下のように近似できる。nj(ζj+θ0)=nj+1(ζj+1+θ0)ここで、ζj+ξj=π2かつζj+1+ξj+1=π2であるから、sinζj=cosξjかつsinζj+1=cosξj+1を上記式に各々代入して、nj(cosξj+θ0)=nj+1(cosξj+1+θ0)nj(cosξj+xjRj)=nj+1(cosξj+1+xjRj)njcosξj+njxjRjnj+1xjRj=nj+1cosξj+1njcosξj+xj(njnj+1)Rj=nj+1cosξj+1nj+1cosξj+1=njcosξj+xj(njnj+1)Rjとの近似結果を得る(屈折率nと光線方向余弦cosξとの積の屈折面における変換則)。

p36 屈折行列

屈折行列Rj
レンズ光軸をz軸、光軸からの距離をxj、光線がx軸となす角をxij、屈折率をnjとした場合、曲率半径Rjの屈折面前後での基底ベクトル(光軸距離、伝播方向)について、Qj+1=RjQj (xj+1nj+1cosξj+1 )= (10nj+1njRj1 )(xjnjcosξj )=(xjxj(nj+1nj)Rjnjcosξj )であり、屈折面前後で、光軸からの距離が不変であること(行列1行目)xj+1=xjまた、屈折率nと光線方向余弦cosξとの積ncosξにかかる変換則(行列2行目)nj+1cosξj+1=xj(nj+1nj)Rjnjcosξjnj+1cosξj+1=njcosξj+xj(njnj+1)Rjが行列法で再現記述されている。

p36 転送行列

転送行列Tj+1
曲面間の光線伝播の際の光軸からの距離xjについて、光軸上の球面間隔をsj+1とすると、以下の近似式が成立しているところ、xj+1xjsj+1cosξj+1転送行列Tj+1を用いた記述法でも、
Qj+1=Tj+1Qj (xj+1nj+1cosξj+1 )= (1sj+1nj+101 )(xjnjcosξj )=(xj+sj+1nj+1×(njcosξj)njcosξj )=(xj+sj+1njnj+1cosξjnj+1cosξj )となるところ、2曲面間では屈折率に変動はなく(nj+1=nj)、また、光線のx軸角度にも変動がないこと(ξj=ξj+1)を踏まえると、行列1行目は、xj+1=xj+sj+1nj+1nj+1cosξj+1xj+1xj=sj+1 cosξj+1となり、上記近似式が再現記述されている(行列2行目は左辺・右辺ともにnj+1cosξj+1となり常に成立)。

p37 システム行列の行列式

屈折行列Rj=(10nj+1njRj1)の行列式の値は、|Rj|=1×1+0×nj+1njRj=1転送行列Tj+1=(1sj+1nj+101)の行列式の値は、|Tj+1|=1×1+sj+1nj+1×0=1一般に、正方行列ABの積ABの行列式は各行列の行列式の積に等しい|AB|=|A||B|から、システム行列S=R2TLR1の行列式についても|S|=|R2TLR1|=|R2||TL||R1|=1×1×1=1となる。

p38 物体から像までの光線伝播式

Qim=DQobにおける行列D=TimSTobにおいて、システム行列Sの行列式の値は1であり、転送行列TobTimの行列式の値も計算をするといずれも1であるから、Dの行列式の値も1となり、|D|=1(cas2n2)(b+as1n1)+0×a=1cas2n2=1b+as1n1

p38 横倍率

横倍率M
式3.43 を展開した行列の1行目は、xim=xob(cas2n2)であるので、横倍率Mは、Mximxob=cas2n2となる。ここで、、Dの行列式の値が1であることから、 cas2n2=1b+as1n1であるので、 b+as1n1=1cas2n2=1Mよって、D=(M0a1M)

p39 角倍率

角倍率γ
γtanζ2tanζ1=1tanξim1tanξob=cotξimcotξob近軸光線では、p29図3.2において、cotξob=|H1A1||O1H1||H1A1||O1A1|=cosξob同様に、cotξimcosξimなので、γcosξimcosξob
なお、横倍率と角倍率の積は、レンズ前後の屈折率の比に等しく(p34式3.35)Mγ=n1n2なので、γ=n1n21M=n1n2(b+as1n1)

p40 焦点位置

式3.52 より、n1a=fAなので、焦点位置は、sF1=bn1a=b(fA)=(1+nLn2R2dLnL)fA

p41 薄肉レンズのシステム行列

主点とレンズ頂点が一致する(p37図3.7 でV1=H1V2=H2)ので、{sH1=0sH2=0また、レンズは空気中にあるので、{n1=1n2=1これを式3.46{sH1=n1(1b)asH2=n2(c1)aに代入すると{0=1(1b)a0=1(c1)aから、{b=1c=1が得られる。よって、薄肉レンズのシステム行列は、Sj=(10aj1)

p42 色収差

合成系の後側焦点距離fcについての公式3.601fc=j=1N1fjの両辺を波長λで微分することを考える(変数はλnjfの入れ子構造)。
逆数の微分の公式δ(1f)=δff2より、δ(1fc)=δ(j=1N1fj)δfcfc2=j=1Nδfjfj2ここで、右辺のの中身を考えると、(fjfと簡略表記して)δff2=1fδff=1f(1f)(1δf)=1f(nj1)(1R1j1R2j)δ((nj1)(1R1j1R2j))=1fnj1δ(nj1)=1f(nj1)dnjdλddnj(1nj1)=1f(nj1)dnjdλ(ddnj(nj1)(nj1)2)=1f(nj1)dnjdλ(1(nj1)2)=1f1nj1dnjdλ=1f1nj1δnjよって、fjfと簡略表記していたのを戻すと、右辺のの中身は、1fjδnjnj1である。従って、色収差δfcについて、以下の式が成立する。δfcfc2=j=1N1fjδnjnj1=j=1N1fjδnjnj1

p43 薄肉レンズ2枚の色消し条件

レンズ2枚であるので、j=1,2を式3.60 及び式3.63に用いると、以下の式が得られる。{1fc=1f1+1f21f1ν1+1f2ν2=0第2式より、1f1ν1=1f2ν2f1ν1=f2ν2f1f2=ν2ν1となるので、これを第1式に用いると、1fc=f2+f1f1f2=1+f1f2f1=1ν2ν1f1f1=fc(1ν2ν1)=fc(ν1ν2ν1)第2式より、f2=f1ν1ν2=fc(ν1ν2ν1)ν1ν2=fc(ν1ν2ν2)

p45 球面反射鏡における結像式

p30式3.21において、n2=nとおくと、ns1ns2=nnR1s11s2=2R=1fとなる。

p46 横倍率

式3.65 より、1s11s2=1f(1s1+1f)=1s2f+s1s1f=1s2s1ff+s1=s2ff+s1=s2s1よって、M=ximxob=|PV||QV|=s2s1=ff+s1

p46 放物線

放物線は焦点Fと準線lからの距離が等しい点の軌跡であるところ、放物線の形状を仮に、x=ay2と置いた場合、原点Oと焦点Fとの間の距離はfであり、原点Oと準線lすなわちx=lへの垂線距離もfである。ここで、x=fのときの放物線上の対応点をAとすると、点Ay軸からfだけ離れているので、y軸から更にfだけ離れた準線lとの距離は、2fとなり、放物線の定義から、点Aと焦点Fとの距離も2fである。よって、点Aの座標は、(x,y)=(f,2f)となる。これをx=ay2に代入して、f=a(2f)2f=4af2a=14fしたがって、放物線の方程式は、x=14fy2y2=4fx