光学学报, 2021, 41 (5): 0526001, 网络出版: 2021-04-07   

会聚波面光场的衍射特性 下载: 914次

Diffraction Characteristic of Light Field with Converging Wave Front
作者单位
福建师范大学光电与信息工程学院激光与光电子技术研究所, 光子技术福建省重点实验室, 医学光电科学与技术教育部重点实验室, 福建 福州 350007
摘要
以能量守恒定律为基础分析一种会聚光场的衍射积分公式,推导出以一个衍射源参考波面为半球面的衍射积分公式的倾斜因子,给出一个计算焦平面上光场分布的表达式,并以虚拟矢量子波干涉解释会聚光场的衍射过程。根据满足正弦条件的聚焦物镜对光波波面的变换特性,给出平坦波面光场经过物镜聚焦后的焦面光场计算表达式,并以简单直角坐标偏振和简单柱面坐标偏振的拉盖尔-高斯光束聚焦特性为例说明该公式的应用。
Abstract
On the basis of the law of energy conservation, we analyzed a diffraction integral formula for the converging light field, derived an inclination factor of the diffraction integral formula for a semi-spherical reference wave front of the diffraction source, and suggested a formula for calculating the light field distribution on the focal plane. In addition, we explained the diffraction process of the converging light field by the interference process of suppositional vector wavelets. Then, according to the wave front transformation properties of focusing objective lens satisfying the sine condition, a calculation formula for the light field in focal plane which focuses from a planar wave front by the focusing objective lens was proposed. Finally, the applications of the above formulas were discussed through the focusing characteristics of the Laguerre-Gaussian beams with simple rectangular coordinate polarization and simple cylindrical coordinate polarization.

1 引言

会聚波面为聚焦光学系统中的常见波面,在考察单色会聚球面波通过圆孔衍射后在焦点附近的光场特性时,由于忽略了衍射源波面上倾斜因子的影响,惠更斯-菲涅耳(Huggens- Fresnel)原理[1-2]由德拜(Debye)积分[2-4]或里查德-沃耳夫(Richards-Wolf)衍射积分公式[5]表达。

柱面坐标偏振光场由于具有一些如光镊子、粒子加速和高分辨率成像等潜在的特殊应用,光场表达和传输特性分析及其应用研究已成为热点课题[6-11]。德拜积分或里查德-沃耳夫矢量衍射积分公式常被用于各类偏振光束聚焦特性的数值计算和模拟[10-15],由于忽略了衍射积分公式中的观察方倾斜因子,基于数值计算的电场强度或光强分布图成为焦面光场的常用表达形式,但能够用于焦面光场特性分析的解析函数较少。

本文以矢量光场特性和行波场辐射能量守恒为基础分析适用于会聚光场特性分析的衍射积分公式,将分析衍射问题的衍射源参考面设定于半球形波面,将观察参考面设定于焦平面,推断出用于分析会聚光场焦平面特性的衍射积分公式中的观察方倾斜因子,观察方倾斜因子为观察方倾斜角余弦的平方根,并以虚拟矢量子波相干叠加解释会聚光场的衍射过程。

根据聚焦物镜对光波波面的变换功能,当入射平坦波面垂直于聚焦物镜的光轴时,满足正弦条件的聚焦物镜的变迹(或切趾)因子[3,5]为物镜像方光线倾斜角余弦的平方根,基于此给出一种平坦波面光场经过满足正弦条件的聚焦物镜聚焦时的焦面光场表达式,并以一个简单直角坐标偏振LG00模高斯(Gauss)光束和简单柱面坐标偏振LG10模拉盖尔(Laguerre)-高斯光束[9,16]的聚焦特性为例说明该公式的应用。

2 观察方倾斜因子的选择

根据惠更斯原理[2]分析衍射过程,衍射源参考面Σ为会聚光场的半球形波面,波面曲率半径为R,曲率中心为O',观察参考面Σ'为焦平面处的平坦波面,它与直角坐标系的x'O'y'平面重合,如图1所示。

图 1. 会聚波面衍射过程示意图

Fig. 1. Schematic of the diffraction from converging wave front

下载图片 查看所有图片

设衍射过程位于无源的各向同性均匀介质空间区域,波长为λ的光波在区域中的电场和磁场分布函数具有相同的表达形式,参考子波波矢k与衍射源参考点P处光波波矢ks的夹角为衍射源方倾斜角θs,参考子波波矢k与观察参考点P'处光波波矢k0的夹角为观察方倾斜角θo,并假设衍射积分公式的倾斜因子K(θs,θo)由衍射源方倾斜角θs和观察方倾斜角θo共同决定[16],根据惠更斯-菲涅耳原理[2],当衍射源参考点P与观察参考点P'之间的距离Lλ时,观察参考面Σ'P'点的光场分布函数U(P')可表达为

U(P')=ΣU(P)exp(ikL)iλLK(θs,θo)dσ,(1)

式中:U(P)为衍射源参考面Σ上的光场分布函数,i为虚数单位,k为光波的波数,dσ为衍射源参考面的微面元。

若焦平面Σ'上的光场有效径向尺寸r'max较小,衍射源波面Σ的曲率半径R≫2r'max2,衍射源参考点P处的光波波矢ksz轴的夹角为θ',可借鉴瑞利判断简化衍射源参考点P(Rsin θ'cos φ,Rsin θ'sin φ,-Rcos θ')与观察参考点P'(r'cos φ',r'sin φ',0)之间的距离L的表达式,即LR-r'sin θ'cos(φ-φ'),简化引起的误差远小于λ/4。当Rr'max时,衍射源方倾斜角θs≈0,观察方倾斜角θo≈-θ',可以认为倾斜因子K(θs,θo)几乎与衍射源方倾斜角θs无关,只由观察方倾斜角θo计算,假设K(θo)为偶函数,有K(θs,θo)≈K(θo)≈K(θ')。

由于衍射源微面元dσ=R2sin θ'dθ'dφ,由(1)式知,观察平面Σ'上的光场径向分布函数U(r',φ')可由衍射源参考半球面上的光场角向分布函数U(R,θ',φ)表达为

U(r',φ')=Rexp(ikR)iλ0π/202πK(θ')U(R,θ',φ)×exp[-ikr'sinθ'cos(φ'-φ)]sinθ'dφdθ'(2)

采用场分布函数Um(r')表达两类偏振光场,下标m标注偏振类型[16]m=0表示简单直角坐标偏振光场,Um(r')代表光场的x轴或y轴分量,或x轴分量与y轴分量的组合;m=1表示光场为简单柱面坐标偏振光场,Um(r')代表光场的径向或角向分量,或径向分量和角向分量的组合。根据m阶贝塞尔(Bessel)函数Jm(α)的积分形式表达式[17-19],观察参考平面Σ'上光场的径向空间频率表达式f'r=sin θ'/λ及其微分形式表达式df'r=(cos θ'/λ)dθ',观察平面Σ'上的两类偏振光场的径向分布函数Um(r')可由衍射源参考半球面Σ上的光场角向分布函数Um(R,θ')表达为

Um(r')=2πexp(iΦ)0K(θ')Um(R,θ')cosθ'circ(λf'r)×Jm(2πr'f'r)f'rdf'r,(3)

式中,Φ=kR-(m+1)π/2为相位延迟,θ'=arcsin(λf'r),用于θ'f'r之间的变量替换, circ(λf'r)为圆孔函数。

衍射源参考面为半球形波面,倾斜角θ'≤π/2,焦面光场的径向空间频率f'r≤1,(3)式中的圆孔函数circ(λf'r)=1。

根据数学概念上的空间频谱定义,焦面光场的径向空间频谱Sm(f'r)为场分布函数Um(r')的m阶汉克尔(Hankel)变换,根据汉克尔变换[9,16,20]的性质,由(3)式知,焦平面光场的径向空间频谱可表达为

Sm(f'r)=RλK(θ')Um(R,θ')exp(iΦ)cosθ'(4)

根据帕斯瓦尔(Parseval)定理[1],基于参考平面上光场径向分布函数Um(r')计算的焦面光场总功率与基于数学概念上的径向空间频谱分布函数Sm(f'r)计算的焦面光场总功率等效,由基于径向空间频谱分布函数计算的光场总功率表达式[16]和(4)式知,在观察参考平面Σ'上,焦面光场总功率POP的计算式为

POP=2πC00R2λ2K2(θ')Um(R,θ')2cos2θ'f'rdf'r,(5)

式中:C0常数[9,16], Um(R,θ')为光场分布函数Um(R,θ')的模。

衍射源波面为半球面,当波面曲率半径Rλ时,波面的尺寸远大于λ,衍射源微面元dσ上的能流密度可由坡印亭(Poynting)矢量[2,21]时间平均值计算,衍射源微面元上的光场功率dPDS=C0Um(R,θ')2dσ,由微面元dσ、空间频率f'r及其微分形式df'r的表达式知,衍射源参考面Σ上的光场总功率PDS可表达为

PDS=2πR2λ2C00Um(R,θ')2cosθ'f'rdf'r(6)

由(5)式和(6)式知,当且仅当衍射积分公式的倾斜因子K(θ')=cos1/2θ'时,焦平面上的光场总功率POP等于半球形衍射源参考面上的光场总功率PDS,满足能量守恒定律。

而当采用德拜积分或里查德-沃耳夫衍射积分公式分析焦面光场特性时,由于未考虑子波干涉叠加时的矢量特性,忽略了观察方倾斜角θo对倾斜因子K(θs,θo)的影响,认为倾斜因子K(θ')=1,由(5)式和(6)式知,基于K(θ')=1计算的焦平面上的光场总功率POP大于半球形衍射源参考面上的光场总功率PDS

可见,对于半球形波面衍射源,在分析焦面光场特性时,由于衍射源方倾斜角θs≈0,倾斜因子K(θs,θo)由观察方倾斜因子K(θo)表达,观察方倾斜因子的合理表达式为K(θo)=cos1/2θo,基于该倾斜因子表达的衍射场符合行波场辐射能守恒定律。

作为一种解释,本文以矢量光场干涉特性、惠更斯-菲涅耳原理和能量守恒定律为基础,认为在观察参考点P'处参加干涉叠加的虚拟子波为矢量形式的子波,参考子波波矢k方向、振幅为u(θo)的矢量子波对P'处观察参考波面微面元dσ'法线(波矢ko方向)上的能流贡献等效于振幅为u(θo)cos1/2θo标量子波的贡献。以观察参考点P'处振动方向为O'y'平面内的某个子波电场E为例,由于EHk三者互相垂直,当子波波矢kz轴的夹角为θo时,子波磁场Hx'O'z'平面内,Hz轴的夹角为π/2-θo,不同衍射源参考点P产生的子波波矢k方向不同,即使子波电场E方向相同,子波磁场H方向不同,在干涉叠加时的贡献也不同。若子波波矢k方向振幅为u(θo)的矢量子波能流密度为s(θo),则它对P'处观察参考波面微面元dσ'法线方向的能流贡献为能流密度在微面元法线上的分量s(θo)cos θo,与振幅为u(θo)cos1/2θo的标量子波的能流贡献等效,表现为观察方倾斜因子K(θo)=cos1/2θo

可见,衍射过程可以认为是虚拟矢量子波相干叠加的过程。

3 平坦波面光场的聚焦特性

设像方焦距为f'的聚焦物镜的物、像双方介质均为空气,光场径向分布函数为Um(r)的平坦波面简单直角坐标偏振光场或简单柱面坐标偏振光场与聚焦物镜的物方焦面ΣF重合,如图1所示。若Um(r)表达的平坦波面光场的径向尺寸较大,聚焦物镜物方主面位于物方焦面光场的菲涅耳衍射深区,聚焦物镜物方主面处的光场径向分布函数Um(rH)可表达为Um(rH)=Um(r)exp(ikf')。

若聚焦物镜的成像条件满足正弦条件,物方焦面径向坐标r与像方光线倾斜角θ'的关系为r=f'sin θ',聚焦物镜的变迹因子[3,5]a(θ')=cos1/2θ',根据透镜对波面的变换特性,平坦波面光场经过聚焦物镜变换后的光场波面为立体角小于2π球面度的球冠面,当参考波面Σ的顶点与聚焦物镜像方主点重合时,波面曲率半径R=f',波面曲率中心与聚焦物镜的像方焦点重合。若聚焦物镜的数值孔径为NA,聚焦物镜的最大像方孔径角θ'max=arcsin(NA),则参考波面Σ上的光场分布函数Um(R,θ')可由聚焦物镜的变迹因子a(θ')和物方焦面处的光场分布函数Um(r)表达为

Um(R,θ')=circθ'θ'maxUm(r)cos1/2θ'exp[ik(f'+s)],(7)

式中:s为聚焦物镜物、像双方主面之间的光程。

当聚焦物镜的焦距f'与像方焦面光场的有效径向尺寸r'max关系满足f'≫2r'max2时,聚焦物镜像方焦面的光场波面可近似为小平面,根据满足正弦条件的聚焦物镜物、像双方的坐标变换关系,像方焦面的空间频率f'r与物方焦面的径向坐标r的关系为f'r=r/(λf'),由于倾斜因子K(θ')=cos1/2θ',由(3)式和(7)式知,聚焦物镜像方焦面的光场径向分布函数Um(r')可表达为

Um(r')=2πλf'exp(iΦ1)0circrrNA×Um(r)Jm2πr'rλf'rdr,(8)

式中:Φ1=k(2f'+s)-(m+1)π/2为相位延迟,rNA=NA·f'为由聚焦物镜数值孔径NA和焦距f'决定的物方焦面光场最大径向尺寸。

当聚焦物镜物方焦面光场的最大径向尺寸大于rNA时,聚焦物镜的通光孔径限制参加变换的光场径向尺寸,聚焦物镜像方焦面的光场分布函数Um(r')为circ(r/rNA)Um(r)的m阶汉克尔变换形式。

当聚焦物镜物方焦面光场的最大径向尺寸小于rNA时,圆孔函数circ(r/rNA)=1,根据满足正弦条件的聚焦物镜物、像双方的坐标变换关系,物方焦面的空间频率fr与像方焦面的径向坐标r'的关系为fr=r'/(λf'),聚焦物镜像方焦面的光场径向分布函数Um(r')由物方焦面光场的径向频谱函数表达为

Um(r')=2πλf'exp(iΦ1)0Um(r)Jm(2πrfr)rdr(9)

由帕斯瓦尔(Parseval)定理[1]、(9)式和空间频率fr与径向坐标r'的关系式知,基于径向分布函数Um(r')计算的聚焦物镜像方焦面光场总功率等于基于径向分布函数Um(r)计算的聚焦物镜物方焦面光场总功率,满足能量守恒定律。

(8)式和(9)式与傅里叶光学导论[1]中基于菲涅耳或傍轴近似条件下表达的傅里叶变换透镜的传输特性[1]基本一致,(9)式与满足正弦条件的准直物镜在数值孔径大于物方焦面最大径向空间频率与波长乘积时的物、像双方焦面光场关系式[16]一致。可见,观察方倾斜因子K(θ')=cos1/2θ'具有合理性,它与聚焦物镜的变迹因子a(θ')=cos1/2θ'将物、像双方焦面光场的m阶汉克尔变换关系从傍轴领域推广到满足正弦条件的聚焦物镜领域。

若聚焦物镜物方焦面处的光场为平坦波面的简单直角坐标偏振LG00模或简单柱面坐标偏振LG10模拉盖尔-高斯光束[22]束腰光斑,光场分布函数表达为

Um(r)=Cm,11ω0rω0mexp-r2ω02,(10)

式中:Cm,1为常数,ω0为基模高斯光束束腰光斑半径。

由(8)式和(10)式知,聚焦物镜像方焦面的光场分布函数Um(r')可表达为

Um(r')=2πCm,1λf'exp(iΦ1)0rNArω0m+1×exp-r2ω02Jm2πr'rλf'dr(11)

当聚焦物镜的数值孔径NAω0/f',且聚焦物镜物方焦面的基模高斯光束瑞利距离zRf'时,由(9)式和(10)式知,聚焦物镜像方焦面处的光场分布Um(r')可表达为

Um(r')=Cm,1exp(iΦ1)1ω'0r'ω'0mexp-r'2ω'02,(12)

式中,ω'=λf'/ω0)为聚焦物镜像方焦面处的基模高斯光束束腰光斑半径。

(12)式与傍轴条件下的高斯光束变换特性一致,这表明观察方倾斜因子由观察方倾斜角余弦平方根表达具有合理性。

4 结论

根据矢量光场特性、能量守恒定律和空间频谱理论分析衍射积分公式,推断出用于计算半球形波面衍射源的焦面光场分布的衍射积分公式倾斜因子K(θs,θo),由观察方倾斜角θo表达为K(θs,θo)=cos1/2θo,并基于此给出计算会聚波面的焦面光场分布的表达式,表明衍射过程可以认为是虚拟的矢量子波相干叠加的过程。结合聚焦物镜对光场的变换功能,给出平坦波面光场经过满足正弦条件的聚焦物镜后的焦面光场表达式,当物方焦面上的光场最大径向尺寸小于NA·f'时,聚焦物镜像方焦面光场分布函数与傍轴条件下的聚焦物镜像方焦面光场分布函数一致,进一步论证了观察方倾斜因子为观察方倾斜角余弦平方根的合理性。

参考文献

[1] Goodman JW. Introduction to Fourier optics[M]. 2nd ed. New York: The McGraw-Hill Companies, Inc., 1996: 1- 106.

[2] BornM, WolfE. Principles of optics[M]. 7th ed. Oxford: Cambridge University Press, 1999: 116-203, 412- 516.

[3] MinG. Imaging with a high numerical-aperture objective[M] ∥Advanced Optical Imaging Theory. Berlin, Heidelberg: Springer, 2000: 143- 176.

[4] Wolf E. Electromagnetic diffraction in optical systems I. An integral representation of the image field[J]. Proceedings of the Royal Society of London Series A Mathematical and Physical Sciences, 1959, 253(1274): 349-357.

[5] Richards B, Wolf E. Electromagnetic diffraction in optical systems, II. Structure of the image field in an aplanatic system[J]. Proceedings of the Royal Society of London Series A Mathematical and Physical Sciences, 1959, 253(1274): 358-379.

[6] 陈建, 詹其文. 矢量光场与激光焦场定制[J]. 光学学报, 2019, 39(1): 0126002.

    Chen J, Zhan Q W. Tailoring laser focal fields withvectorial optical fields[J]. Acta Optica Sinica, 2019, 39(1): 0126002.

[7] Nieminen T A, Heckenberg N R, Rubinsztein-Dunlop H. Forces in optical tweezers with radially and azimuthally polarized trapping beams[J]. Optics Letters, 2008, 33(2): 122-124.

[8] Zhan Q W. Cylindrical vector beams: from mathematical concepts to applications[J]. Advances in Optics and Photonics, 2009, 1(1): 1-57.

[9] 郭福源, 李连煌, 郑华, 等. 简单柱面坐标偏振光束的传输特性[J]. 中国激光, 2013, 40(8): 0802003.

    Guo F Y, Li L H, Zheng H, et al. Propagation characteristics of simple cylindrical polarized beam[J]. Chinese Journal of Lasers, 2013, 40(8): 0802003.

[10] Wang F, Zhao C L, Dong Y, et al. Generation and tight-focusing properties of cylindrical vector circular Airy beams[J]. Applied Physics B, 2014, 117(3): 905-913.

[11] Janet C A P, Rajesh K B, Udhayakumar M, et al. Tight focusing properties of radially polarized Gaussian beams with pair of vortices[J]. Chinese Physics Letters, 2016, 33(12): 124206.

[12] Jiang Y S, Li X P, Gu M. Generation of sub-diffraction-limited pure longitudinal magnetization by the inverse Faraday effect by tightly focusing an azimuthally polarized vortex beam[J]. Optics Letters, 2013, 38(16): 2957-2960.

[13] Youngworth K, Brown T. Focusing of high numerical aperture cylindrical-vector beams[J]. Optics Express, 2000, 7(2): 77-87.

[14] 李维超, 杨艳芳, 何英, 等. 基于高阶角向偏振拉盖尔高斯涡旋光束强聚焦的三光链结构[J]. 光学学报, 2019, 39(8): 0826002.

    Li W C, Yang Y F, He Y, et al. Triple optical chain generated by tight focusing of azimuthally polarized higher-order Laguerre-Gauss vortex beams[J]. Acta Optica Sinica, 2019, 39(8): 0826002.

[15] 李浩然, 樊承锦, 党锦超, 等. 高次方涡旋光束经大数值孔径透镜的聚焦特性[J]. 强激光与粒子束, 2018, 30(1): 011002.

    Li H R, Fan C J, Dang J C, et al. Focusing properties of power-exponent-phase vortex beam focused by high numerical-aperture objective[J]. High Power Laser and Particle Beams, 2018, 30(1): 011002.

[16] 郭福源, 李连煌, 郑华. 小平坦波面光场的衍射与准直特性[J]. 光学学报, 2021, 41(4): 0426001.

    Guo F Y, Li L H, Zheng H. The characteristics of collimation and far field diffraction from a small light field with planar wave front[J]. Acta Optica Sinica, 2021, 41(4): 0426001.

[17] 刘庆伟, 李重光, 楼宇丽, 等. 双曲线顶点对轴锥镜聚焦柱矢量光束的影响[J]. 中国激光, 2020, 47(4): 0409001.

    Liu Q W, Li C G, Lou Y L, et al. Influence of hyperbolic vertex on axicon-focusing cylindrical vector beams[J]. Chinese Journal of Lasers, 2020, 47(4): 0409001.

[18] Saghafi S. Sheppard C J R. Near field and far field of elegant Hermite-Gaussian and Laguerre-Gaussian modes[J]. Journal of Modern Optics, 1998, 45(10): 1999-2009.

[19] 黎芳, 宫瑞婷, 刘慧. 非傍轴拉盖尔-高斯光束的轨道角动量密度特性[J]. 激光与光电子学进展, 2017, 54(3): 032601.

    Li F, Gong R T, Liu H. Orbital angular momentum density characteristics of nonparaxial laguerre-Gaussian beams[J]. Laser & Optoelectronics Progress, 2017, 54(3): 032601.

[20] 王奇, 杨天龙, 沈文达. 弱波导光纤出射模在衍射远区的场分布[J]. 光学学报, 1987, 7(1): 62-66.

    Wang Q, Yang T L, Shen W D. Far field patterns of the modes emitted from a weakly guiding optical fiber[J]. Acta Optica Sinica, 1987, 7(1): 62-66.

[21] 陈军. 光学电磁理论[M]. 北京: 科学出版社, 2005: 35- 74.

    ChenJ. Optical electromagnetic theory[M]. Beijing: Science Press, 2005: 35- 74.

[22] Deng D M. Nonparaxial propagation of radially polarized light beams[J]. Journal of the Optical Society of America B, 2006, 23(6): 1228-1234.

郭福源, 李连煌, 郑华. 会聚波面光场的衍射特性[J]. 光学学报, 2021, 41(5): 0526001. Fuyuan Guo, Lianhuang Li, Hua Zheng. Diffraction Characteristic of Light Field with Converging Wave Front[J]. Acta Optica Sinica, 2021, 41(5): 0526001.

本文已被 2 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!