小平坦波面光场的衍射与准直特性 下载: 843次
1 引言
1678年,惠更斯基于直观信念[1]提出惠更斯原理,该原理以球面子波的次级扰动中心和包络面解释波传播过程中两个波面之间的关系。1818年,菲涅耳用干涉理论补充了惠更斯原理,形成惠更斯-菲涅耳原理[1-2],该原理以振幅随方向变化的子波相干叠加解释衍射现象。在考察单色点光源对空间任意一点的光作用时,衍射源参考面为点源产生的波面,惠更斯-菲涅耳原理数学表达式中的倾斜因子假设[2]明确了子波的传播区域。
目前,常见的以场分布为基础的衍射积分公式有3个,分别为菲涅耳-基尔霍夫衍射公式、第一种和第二种瑞利-索末菲衍射公式[1-6]。它们由以亥姆霍兹-基尔霍夫积分定理[1-2]为基础,在球面曲率半径足够大时大球面部分对积分公式的贡献为零条件下[2](或在索末菲辐射条件下[1])简化衍射积分公式,并根据所选择的格林函数和基尔霍夫边界条件[1-2]的不同应用推导而得。菲涅耳-基尔霍夫衍射公式是第一种和第二种瑞利-索末菲衍射公式的数学平均值,目前,3个衍射积分公式同时存在于同一电磁场理论体系中,尽管开展了一些衍射积分公式的比较研究[4-7],但衍射过程的能量问题未被充分关注,3个衍射积分公式的计算精度优劣仍未明确[1]。其中,第一种瑞利-索末菲衍射公式的倾斜因子符合惠更斯-菲涅耳原理关于倾斜因子的假设,被认为是较为合理的衍射积分公式,被应用于光场特性分析中[8-12]。以平面波角谱(空间频谱)或场分布傅里叶变换为基础的衍射积分[1,13-18]也被用于衍射场表达。文献[ 1]以有限频带的空间滤波器解释光波的传播,以倏逝波解释横向空间频率大于1/λ的光场;文献[ 13]采用傅里叶变换公式和卷积定理论证了两种瑞利-索末菲衍射公式;文献[ 14-17]分别运用鞍点法求积分的渐进公式和空间脉冲响应表达式等公式论证第一种瑞利-索末菲衍射公式。
通常,单模光纤[8-9,19]的模场半径、强聚焦贝塞尔-高斯光束[20]和低阶拉盖尔-高斯光束[21-22]的束腰光斑半径与光波波长为同一数量级,光纤端面光场和各类强聚焦高斯光束束腰光场的波面为平面。这些径向尺寸较小且波面为平面的光场可称为小平坦波面光场,小平坦波面光场的衍射光束为非傍轴光束,衍射积分公式的倾斜因子将影响其衍射远场分布的表达。
光波为矢量波,波面的法线方向(波矢)是波的传播方向,本文以惠更斯原理为基础,将用于分析衍射问题的衍射源参考面和衍射远场观察参考面均设定为波面,并以惠更斯-菲涅耳原理的数学表达式为基础,分析波面为平面的简单直角坐标偏振和简单柱面坐标偏振[23]小尺寸圆形衍射源光场的远场衍射特性。根据空间频谱理论和能量守恒定律,推断出一种远场衍射积分公式的倾斜因子,给出一种用于分析小平坦波面远场衍射特性的衍射积分公式表达,并以虚拟余弦辐射体辐射的子波相干叠加解释衍射过程。同时,给出小平坦波面光场经准直物镜后的像方焦面光场表达式,并以一个简单直角坐标偏振LG00模和简单柱面坐标偏振LG10模拉盖尔-高斯光束变换特性为例说明该公式的应用。
2 倾斜因子的选择
受开孔衍射屏限制的衍射过程可分解为开孔衍射屏对入射光场的制约和从衍射屏通光孔出射的光场的传播。若将从衍射屏通光孔出射的光场称为衍射源,可将受限衍射过程简化为相干光场从衍射源参考面到观察参考面的传播过程。若不存在衍射屏,将光波传播过程中被选择作为分析对象起始参考位置的光场称为衍射源,相干光场的行波传播过程即为无受限衍射过程。
参照惠更斯原理[2]分析衍射过程,衍射源参考面和观察参考面均为光波的波面,设衍射源的光波波面Σ为平面,该平面与直角坐标系的xOy平面(或柱面坐标系的rOφ平面)重合,衍射源参考面上任一参考点P处的光波波矢ks均与z轴平行。若小平坦波面衍射源光场分布为旋转对称形式,其中心与坐标原点O重合,当衍射源参考点P与观察参考点P1的距离L足够大时,衍射远场的光波波面Σ1可以近似为以O为圆心、以R为半径的半球面。可设该半球形波面为观察参考面,在半球面Σ1上,观察参考点P1处的光波波矢ko与z轴的夹角为θ。设衍射源方倾斜角θs为参考子波波矢k与衍射源参考点P处微面元dσ法线的夹角,观察方倾斜角θo为参考子波波矢k与观察参考点P1处微面元dσ1法线的夹角。角度的符号规则是以波面法线为始边沿锐角方向转向参考线,逆时针则为正,顺时针则为负,如
设衍射过程位于无源的各向同性均匀介质空间区域,区域中的时谐电磁场在麦克斯韦方程组中表现为对偶性[24],其电场和磁场的分布形式和传播规律相同。场分布可采用相同函数形式,表达为U(P,t)=U(P)exp(-i2πνt),其中U(P)为P点位置的复数形式场分布函数,ν为光波频率,t为时间变量。根据惠更斯-菲涅耳原理[2],并参照3个衍射积分公式[1-2]的表达形式,假设衍射积分公式的倾斜因子K(θs,θo)由衍射源方倾斜角θs和观察方倾斜角θo共同决定,当衍射源参考点P与观察参考点P1之间的距离L远大于光波波长λ时,观察参考面Σ1上P1点的光场分布函数U(P1,t1)可由衍射源参考面Σ上的光场分布函数U(P,t)表达,公式为
图 1. 小平坦波面光场的远场衍射和准直过程示意图
Fig. 1. Schematic of collimation and far-field diffraction from a small light field with planar wavefront
式中:k为光波的波数;Δt=t1-t;dσ为衍射源的微面元。
设在z=0的衍射源参考平面上的光场相位为零,光场分布为旋转对称结构形式,参考平面上的电场偏振为简单直角坐标偏振或简单柱面坐标偏振[23],前者电场偏振方向为直角坐标系的x轴或y轴或二者组合,后者电场偏振方向为柱面坐标系的径向或角向或二者组合。两类偏振光的场分布函数均与柱面坐标系中的角向坐标φ无关,可由同一形式表达为Um(r,t)=Um(r)exp(-i2πνt),其中m为标注衍射源电场偏振类型的参量。当衍射源为简单直角坐标偏振光场时,m=0,Um(r)代表光场的轴向分量Ux(r)和Uy(r);当衍射源为简单柱面坐标偏振光场时,m=1,Um(r)代表光场的径向分量Ur(r)和角向分量Uφ(r)。
在衍射源平面内,若半径为rmax的圆斑内的光场功率接近于衍射源总功率,可将rmax称为衍射源光场的有效半径。若衍射源光场的有效半径rmax较小,根据瑞利判断,在衍射远场、观察参考半球面半径R≫2
式中:circ(ξ)为圆孔函数,当0≤ξ≤1时,circ(ξ)=1,当ξ>1时,circ(ξ)=0;fr=sin θ/λ为径向空间频率;fr,max为最大径向空间频率。
当衍射源为简单直角坐标偏振光时,数学概念上的径向空间频谱Sm(fr)为衍射源场分布Um(r)的零阶汉克尔变换[19];当衍射源为简单柱面坐标偏振光时,数学概念上的径向空间频谱Sm(fr)为衍射源场分布Um(r)的一阶汉克尔变换[23]。两类偏振光场在数学概念上的衍射源径向空间频谱Sm(fr)可由同一形式的计算式表达,公式为
当光波的参考波面尺寸远大于波长时,参考波面上微面元的能流密度可在几何光学的精度范围内参照平面波坡印廷矢量[2,24,26]时间平均值的计算公式计算,由于观察参考半球面Σ1上参考点P1处的微面元dσ1=R2sin θdθdφ1,场分布函数为Um(R,θ)的光场在微面元dσ1上的功率dPOH计算式为
式中:Um(R,θ)为场分布函数;C0为常数[7,11]。当函数Um(R,θ)代表电场强度分布时,C0=(1/2)
由于径向空间频率fr=sin θ/λ,有dfr=(cos θ/λ)dθ,由(2)~(4)式可知,在观察半球面Σ1上,由衍射源径向空间频谱Sm(fr)表达的衍射远场总功率[7]计算式为
当采用空间频谱分布表达衍射源光场时,由Parseval定理[1]知,在参考平面Σ上,由场分布函数Um(r)表达的衍射源总功率与数学概念上的空间频谱分布Sm(fr)表达的衍射源总功率等效,即衍射源总功率[11]计算式为
设倾斜因子K(θ)=
而当采用菲涅耳-基尔霍夫衍射公式、第一种或第二种瑞利-索末菲衍射公式计算小平坦波面衍射远场总功率时,由于3个衍射公式的倾斜因子可分别表达为KFK(θ)=(1+cos θ)/2,KRS-1(θ)=cos θ,KRS-2(θ)=1,由(5)、(6)式知,无论最大径向空间频率fr,max采用物理意义上的最大径向空间频率或数学概念上的最大径向空间频率,衍射远场总功率计算值POH均偏离衍射源总功率PDS,不满足能量守恒定律。其中,当采用第一种瑞利-索末菲衍射公式计算小平坦波面衍射远场总功率时,KRS-1(θ)=cos θ,无论最大径向空间频率fr,max采用物理意义上的最大径向空间频率或数学概念上的最大径向空间频率,衍射远场总功率计算值POH均小于衍射源总功率PDS。可见,对于波面为平面的小面积衍射源,在分析衍射远场特性时,由于观察方倾斜角θo≈0,衍射源方倾斜角θs≈θ,倾斜因子K(θs,θo)由衍射源方倾斜因子K(θs)表达,衍射源方倾斜因子的合理表达式为K(θs)=
作为一种解释,本文以惠更斯-菲涅耳原理为基础,认为衍射源波面上每个微面元dσ的虚拟辐射是各向同性的,其光度学亮度为常数,即将衍射源波面上的每个面元看成余弦辐射体,球面子波辐射强度随衍射源方倾斜角θs的变化规律遵从Lambert定律[2],即球面子波分布函数在子波波矢k方向的振幅us(θs)与球面子波分布函数在衍射源微面元dσ法线方向的振幅us(0)的关系为us(θs)=us(0)×
3 小平坦波面衍射光束的准直特性
若准直物镜的像方焦距为f',物镜的物、像双方介质均为空气,物、像双方介质的折射率均为1,小平坦波面衍射源与准直物镜的物方焦面重合,半球形参考波面Σ1的曲率半径R=f',如
式中:准直物镜的最大物方孔径角θNA=arcsin(NA),NA为准直物镜的数值孔径。
由于小平坦波面衍射源经过准直物镜准直之后的光场径向尺寸较大,准直物镜像方焦面位于像方主面的菲涅耳衍射深区,准直物镜像方焦面处的光场分布函数Um(r')与准直物镜像方主面处的光场分布函数Um(r'H)的关系为Um(r')=Um(r'H)×exp(ikf')。根据物、像双方的坐标变换关系和径向空间频率的定义,有r'=λfrf'。由(2)、(3)、(7)式知,当倾斜因子K(θ)=
式中:相位延迟Φ=k(2f'+s)-(m+1)π/2;准直物镜像方焦面光场的最大半径r'NA=f'NA。
当圆孔函数circ(r'/r'NA)=1时,将准直物镜像方焦面上的光场分布函数Um(r')记为U0,m(r'),当衍射源的最大径向空间频率为fr,ms时,若准直物镜的数值孔径NA<λfr,ms,则准直物镜的孔径拦截了空间频率为fr>NA/λ的小平坦波面衍射源的光场,准直物镜像方焦面光场分布函数Um(r')=circ(r'/r'NA)U0,m(r')。当准直物镜的数值孔径NA≥λfr,ms时,圆孔函数circ(r'/r'NA)=1,若小平坦波面衍射源为简单直角坐标偏振光场或简单柱面坐标偏振光场,准直物镜像方焦面处的光场分布函数Um(r')=U0,m(r')。U0,m(r')为准直物镜物方焦面处的光场分布函数Um(r)的零阶或一阶贝塞尔-傅里叶变换形式,其结论与Goodman在傅里叶光学导论[1]中基于傍轴条件或菲涅耳近似下表达的傅里叶变换透镜的光场变换特性表达式一致,可见,倾斜因子K(θ)=
若准直物镜物方焦面处的衍射源为小平坦波面的简单直角坐标偏振LG00模或简单柱面坐标偏振LG10模拉盖尔-高斯光束束腰光斑,小平坦波面衍射源光场分布表达为
式中:Cm,1为常数;ω0为基模高斯光束的束腰半径。准直物镜的焦距f'与基模高斯光束的瑞利距离zR=π
式中:ω'0=λf'/(πω0)为准直物镜的数值孔径NA≫θ0时的准直物镜像方焦面处的基模高斯光束的束腰半径,其中θ0=λ/(πω0 )为准直物镜物方基模高斯光束的远场发散角。由(10)式知,当准直物镜的孔径对变换的光束产生限制作用时,准直物镜像方焦面光场分布函数Um(r')由U0,m(r')与圆孔函数circ(r'/r'NA)的乘积表达。
当准直物镜的数值孔径NA≫θ0时,可以认为准直物镜的孔径没有对变换的光束产生限制作用,圆孔函数circ(r'/r'NA)=1,准直物镜像方光束为束腰位于像方焦面且基模高斯光束束腰半径为ω'0的简单直角坐标偏振LG00模或简单柱面坐标偏振LG10模拉盖尔-高斯光束,与傍轴条件下的高斯光束变换特性一致。
4 结论
根据惠更斯-菲涅耳原理、空间频谱理论及能量守恒定律,推断小平坦波面衍射远场积分公式的倾斜因子K(θs,θo)由衍射源方倾斜因子K(θs)表达,衍射源方倾斜因子为衍射源方倾斜角θs余弦的平方根,并给出小平坦波面衍射远场分布的表达式。结合准直物镜对光场的变换功能,给出小平坦波面衍射源经满足正弦条件的准直物镜准直后的光场表达式。当准直物镜的数值孔径大于衍射源的最大径向空间频率fr,ms与波长乘积结果时,准直物镜物、像双方焦面光场分布关系与傍轴条件下的准直物镜物、像双方焦面光场分布关系一致,进一步论证了倾斜因子表达式的合理性。
[1] Goodman JW. Introduction to Fourier optics[M]. 2nd Ed. New York: The McGraw-Hill Companies, 1996: 1- 106.
[2] BornM, WolfE. Principles of optics [M]. 7th Ed. Oxford: Cambridge University Press, 1999: 116-203, 412- 516.
[3] GuM. Advanced optical imaging theory[M]. Berlin: Springer, 2000: 143- 169.
[7] 郭福源, 李连煌. 标量衍射积分公式比较分析[J]. 光学学报, 2013, 33(2): 0226001.
[8] Artiglia M. Coppa G, di Vita P, et al. Mode field diameter measurements in single-mode optical fibers[J]. Journal of Lightwave Technology, 1989, 7(8): 1139-1152.
[9] 林斌, 郭福源, 陈钰清, 等. 弱导光纤的标量衍射光束特性分析[J]. 中国激光, 2003, 30(9): 809-813.
[10] 郭福源, 王明华. 星形光波导耦合器的耦合特性分析[J]. 光学学报, 2006, 26(12): 1797-1802.
[11] 郭福源, 李连煌. 非傍轴衍射光束的传输特性[J]. 中国激光, 2013, 40(1): 0102004.
[12] 朱志绿, 李连煌, 郭福源, 等. 非傍轴LG00模与LG10模高斯光束的远场分布特性[J]. 光学学报, 2017, 37(5): 0526002.
[13] Sherman G C. Application of the convolution theorem to Rayleigh's integral formulas[J]. Journal of the Optical Society of America, 1967, 57(4): 546-547.
[14] 羊国光, 宋菲君. 高等物理光学[M]. 2版. 合肥: 中国科学技术大学出版社, 2008: 21- 86.
Yang GG, Song FJ. Advanced physical optics[M]. 2nd Ed. Hefei: University of Science and Technology of China Press, 2008: 21- 86.
[15] Lucke R L. Rayleigh-Sommerfeld diffraction and Poisson's spot[J]. European Journal of Physics, 2006, 27(2): 193-204.
[16] Harvey J E. Fourier treatment of near-field scalar diffraction theory[J]. American Journal of Physics, 1979, 47(11): 974-980.
[17] 殷琢豫. 22(5): 74-[J]. . 建立严格衍射积分的简单方法. 上海交通大学学报, 1988, 82: 117-118.
. 22(5): 74-[J]. Yin Z Y. A simple method for obtaining strigent diffraction integral. Journal of Shanghai Jiaotong University, 1988, 82: 117-118.
[19] 王奇, 杨天龙, 沈文达. 弱波导光纤出射模在衍射远区的场分布[J]. 光学学报, 1987, 7(1): 62-66.
[20] Youngworth K S, Brown T G. Focusing of high numerical aperture cylindrical-vector beams[J]. Optics Express, 2000, 7(2): 77-87.
[21] 李维超, 杨艳芳, 何英, 等. 基于高阶角向偏振拉盖尔高斯涡旋光束强聚焦的三光链结构[J]. 光学学报, 2019, 39(8): 0826002.
[22] 李红艳, 杨艳芳, 何英, 等. 基于4pi聚焦系统实现特殊焦场结构之间的转换[J]. 光学学报, 2020, 40(12): 1226001.
[23] 郭福源, 李连煌, 郑华, 等. 简单柱面坐标偏振光束的传输特性[J]. 中国激光, 2013, 40(8): 0802003.
[24] 陈军. 光学电磁理论[M]. 北京: 科学出版社, 2005: 35- 74.
ChenJ. Optical electromagnetic theory[M]. Beijing: Science Press, 2005: 35- 74.
[26] 江月松, 徐洪武, 唐华, 等. 厄米-高斯模与拉盖尔-高斯模的图像比较[J]. 光学学报, 2013, 33(6): 0606001.
Article Outline
郭福源, 李连煌, 郑华. 小平坦波面光场的衍射与准直特性[J]. 光学学报, 2021, 41(4): 0426001. Fuyuan Guo, Lianhuang Li, Hua Zheng. Diffraction and Collimation Characteristics of Small Light Fields with Planar Wavefronts[J]. Acta Optica Sinica, 2021, 41(4): 0426001.