光学学报, 2018, 38 (10): 1005002, 网络出版: 2019-05-09   

正六边形多孔阵列的涡旋光衍射 下载: 1132次

Diffraction of Vortex Beam by Regular Hexagonal Multi-Hole Array
作者单位
北京交通大学理学院光信息科学与技术研究所发光与光信息技术教育部重点实验室, 北京 100044
摘要
基于菲涅耳衍射理论,研究了涡旋光经正六边形排布的多孔阵列衍射后的光强分布,分析了正六边形多孔阵列的结构参数对蜂窝状光场的影响。研究结果表明,衍射光场会随着相位结构发生周期性变化,得到蜂窝状或花瓣状的光强分布;圆孔半径会影响衍射光场的范围,正六边形的边长对条纹的宽度和间距存在影响。
Abstract
Based on the Fresnel diffraction theory, the intensity distribution of a vortex beam diffracted by a regular hexagonal multi-hole array is investigated and the effects of the structural parameters of this regular hexagonal multi-hole array on the cellular optical field are also analyzed. The research results show that the diffracted field periodically changes with the phase structure and thus a cellular or petal shaped intensity distribution can be obtained. The radius of the circular hole has an effect on the diffracted field range and the side length of the regular hexagon has influences on the width and the separation distance of the diffraction fringes.

1 引言

涡旋光是一种携带螺旋相位波前的光束,其每个光子都携带有轨道角动量(OAM)。Allen等[1-2]通过实验发现了具有螺旋相位因子的光束,随后涡旋光得到了广泛关注,被应用于光学操控[3]、空间光通信[4]、旋转测量[5]及光学微刻[6-7]等领域。涡旋光因其携带的轨道角动量及螺旋相位结构,在干涉衍射方面具有特别的现象。Sztul等[8]进行了涡旋光的杨氏双缝干涉实验。Ghai等[9]对涡旋光进行了单缝衍射实验,发现通过观察干涉衍射条纹的弯曲变化等现象可确定涡旋光的拓扑荷符号及部分其他信息。Hickmann等[10]研究了三角光阑的涡旋光衍射,衍射图样为均匀排布成三角形阵列的光斑点阵,点阵中的光斑个数与拓扑荷数直接相关。因此,可以利用干涉衍射实验来确定涡旋光拓扑荷的方向及数值大小等信息。

本文主要讨论了涡旋光经多孔光阑后的衍射特性,分析了涡旋光经多孔光阑衍射后所得到的蜂窝状光场特性,研究了光场与涡旋光相位结构及光阑结构参数之间的关系。

2 涡旋光多孔衍射

拉盖尔-高斯(LG)光束是一种比较常用的涡旋光,其复振幅可表示为

uml(r1,θ1,z)=2m!π(m+|l|)!1w(z)r12w(z)lexp-r12w2zLml2r12w2zexp(ilθ1)×expik'r12z2(z2+zR2)exp-i(2m+l+1)arctanzzR,(1)

式中(r1,θ1)为极坐标,r1为极坐标半径,θ1为方位角;z为传播距离;m为径向指数;k'为波数;w(z)=w0[(z2+zR2)/zR2]1/2z处的光斑半径;w0为束腰半径;zR为瑞利距离;(2m+|l|+1)tan-1(z/zR)为Guoy相移项;l为拓扑荷数; Lml(x)为广义拉盖尔多项式。设衍射光阑的窗函数为T(x0, y0),根据菲涅耳衍射理论,LG光束垂直照射衍射光阑中心时,取径向指数m=0,衍射光场的极坐标表示为

E(r,θ,z)=-iλzexp(ik'z)T(r1,θ1)u0l(r1,θ1,0)expik'2zr2+r12-2rr1cos(θ-θ1)rdr1dθ1,(2)

式中r为极径,θ为极角,λ为光波长。

将衍射光阑设置为图1所示的结构,孔径函数可表示为

  T(x1,y1)=1,  |x1-r'22+|y1-3r'22r02 or (|x1|-r')2+y12r020,  else,(3)

图 1. 正六边形多孔阵列示意图及其对具有不同拓扑荷数的LG光束的影响。(a)示意图;(b) l=1; (c) l=2; (d) l=3; (e) l=4; (f) l=5; (g) l=6

Fig. 1. Schematic of regular hexagonal multi-hole array and its effect on LG beams with different topological charge numbers. (a) Diagram; (b) l=1; (c) l=2; (d) l=3; (e) l=4; (f) l=5; (g) l=6

下载图片 查看所有图片

式中(x1,y1)为孔径函数的直角坐标,r'为六边形边长,r0为圆孔半径,r0r'/2。取r0=50 μm,r=0.9 mm,λ=1.55 μm,w0=1 mm,z=1.0 m,将(3)式代入(2)式计算得到LG光经此光阑后的衍射图样,如图1(a)所示。图1(b)、(f)所示为l=1,5时的衍射光场,衍射图样为六边形环状光斑,暗纹衔接成蜂窝结构。图1(c)、(e)所示为l=2,4时的衍射光场,衍射图样为花瓣状光斑,光斑组合成三角形和六边形;图1(d)所示为l=3时的衍射光场,衍射图样为花瓣状,光斑组合成六边形;图1(g)所示为l=6时的衍射光场,场中出现圆环状暗纹,环内存在点状光斑,环外侧较弱光强的条纹彼此衔接,形如莲蓬。

当继续增大拓扑荷数l时,l=6k+n(拓扑荷周期数k=0,1,2,…,n为整数且n<6)的衍射光场与l=n的光场类似,图2所示为拓扑荷数l=7,8,9时的衍射光场,与图1l=1,2,3的光场相似,其中xy为衍射平面的横、纵坐标。将l=1,7的衍射光场彩图转换为灰度图,对图像中的所有像素按照灰度值的大小统计其出现的概率,得到图3(a)、(b)所示的灰度分布直方图,计算其归一化系数巴氏距离DB,计算公式为

DB=xXp(x)q(x),(4)

式中X为灰度值的集合;p(x)、q(x)为两幅图像的灰度值分布概率函数。l=1,7两衍射图间的巴氏距离DB=0.9814(0≤DB≤1),DB越接近1表示这两幅图像越相似[11]。由图1可知,l=1,5时衍射光场均为蜂窝状。分别计算样本ls=1,2,3,4,5,6的衍射图样与高阶拓扑荷l=61(k=0,1,2,…)间的巴氏距离DBlsDBls与拓扑荷周期数k之间的关系如图3(c)所示,其中同种标识的两条曲线分别代表l=6k+1(k=0,1,2,…)与l=6k+5(k=0,1,2,…)的情况。在前四个周期内(k≤4), DB1DB5的值要大于0.95,其他四条曲线的值均远小于0.95,可以认为l=6k±1(k=0,1,2,3,4)的衍射光场与l=1,5的具有较高的相似性。当k>4时,由于拓扑荷数l的增大,光场的蜂窝状结构变差, DB1DB5的数值下降,DB1DB5两条曲线近乎重叠,间接地反映了两者衍射光场分布的近似性,DB2DB4亦是相同的情况。

图 2. 具有不同拓扑荷数l的高阶涡旋光的衍射光场。(a) l=7; (b) l=8; (c) l=9

Fig. 2. Diffraction fields of high order vortex beams with different topological charge numbers l . (a) l=7; (b) l=8; (c) l=9

下载图片 查看所有图片

图 3. 不同拓扑荷数衍射光场的灰度直方图及巴氏距离。(a) l=1;(b) l=7;(c)巴氏距离的变化曲线

Fig. 3. Gray histograms of diffraction fields with different topological charge numbers l and Bhattacharyya distances. (a) l=1;(b) l=7;(c) curves of Bhattacharyya distance

下载图片 查看所有图片

3 蜂窝状光场分析

3.1 拓扑荷数l的影响

图1(b)、(f)所示蜂窝状光场的x方向光强分布,如图4(a)、(d)所示。将中心六边形亮环定义为0级亮条纹,两侧亮环依次为±N级,x=0两侧的主极大值代表0级条纹两条对边的光强值,往外相邻两个主极大对应同一个环状条纹两条对边的光强值,如图4(a)所示。取λ=1.55 μm,w0=1 mm,r0=50 μm,r=0.9 mm,z=1 m。图4(a)~(c)所示为拓扑荷数l=6k+1(k=0,2,4)时的光强分布。随着拓扑荷数l的增大,0级亮环的光强I0逐渐减小,非0级亮环光强值增大,光强包络线的中心逐渐下凹;非0级亮环内外边存在光强差,内侧的光强 IN-in低于外侧的光强 IN-out(N≠0),衍射级数越大,差值越大,且同一衍射级l越大,光强差值也越大。图4(d)~(e)所示为拓扑荷l=6k+5(k=0,2)时的光强分布,也有同样的规律,但此时非0级亮环的外侧光强低于内侧光强。令 L0=I0/IMax,LN = IN-in/ IN-out(N≠0)(计算时均采用归一化光强,最大光强IMax=1)。图4(f)、(g)所示分别为L0L1与拓扑荷数l=6k+1、l=6k+5间的关系,可以看出,拓扑荷数l越大,L0越小;l=6k+1时,L1值恒小于1,且随着拓扑荷数l的增大而逐渐减小;l=6k+5时,L1值恒大于1,且随着拓扑荷数l的增大而逐渐增大。

当拓扑荷数l增大到一定数值时,0级亮环的光强很小,且非0级亮环内外侧光强的差值较大,亮环光强分布不均,因此才会出现图3(c)所示k>4时 DB1DB5数值减小的现象。

图 4. 不同拓扑荷数下衍射光场的横向光强分布及L0和L1。(a) 光强,l=1,k=0; (b) 光强,l=13,k=2; (c) 光强,l=25,k=4; (d) 光强,l=5,k=0; (e) 光强,l=17,k=2; (f) L0和L1, l=6k+1;(g) L0和L1,l=6k+5

Fig. 4. Transverse intensity distributions of diffraction fields and L0 and L1 under different topological charge numbers.(a) Light intensity, l=1,k=0; (b) light intensity, l=13,k=2; (c) light intensity, l=25,k=4; (d) light intensity, l=5,k=0; (e) light intensity, l=17,k=2; (f) L0 and L1, l=6k+1; (g) L0 and L1, l=6k+5

下载图片 查看所有图片

图 5. 不同r0下衍射场的横向光强分布及w1/2与r0间的关系。(a)r0=50 μm; (b)r0=100 μm; (c)r0=150 μm; (d) w1/2与r0的关系

Fig. 5. Transverse intensity distributions of diffraction fields under different r0 and relationship between w1/2 and r0.(a)r0=50 μm; (b)r0=100 μm; (c)r0=150 μm; (d) relationship between w1/2 and r0

下载图片 查看所有图片

3.2 圆孔半径的影响

λ=1.55 μm,w0=1 mm,l=1,r=0.9 mm,z=1.0 m,改变r0得到的衍射光场如图5所示,可以看出,当r0增大时,光场范围逐渐收缩,亮环的数量明显减少。光强分布包络线的半峰全宽w1/2(单位为mm)与r0之间的关系如图5(d)所示,可以看出,随着r0的增大,w1/2减小,但减小的趋势逐渐放缓。

3.3 六边形边长的影响

r0=50 μm,其他参数不变,改变r得到的衍射光场如图6所示,其中D、ΔL分别表示条纹宽度及相邻条纹间隔。随着r的增大,光场范围不会改变,但亮环的尺寸会逐渐减小,条纹变密集。D、ΔLr之间的关系如图6(d)所示,可以看出,随着r的增大,D、ΔL逐渐减小,条纹近乎重叠,说明条纹宽度及条纹间隔数值相等。

图 6. 不同r下衍射场的横向光强分布及条纹宽度和密度与六边形边长间的关系。(a)r=0.4 mm; (b)r=0.6 mm; (c)r=0.8 mm; (d)关系曲线

Fig. 6. Transverse intensity distributions of diffraction fields under different r and relationship of fringe width and density with side length of regular hexagon. (a) r=0.4 mm; (b)r=0.6 mm; (c)r=0.8 mm; (d) relationship curves

下载图片 查看所有图片

4 实验

实验装置图如图7所示。激光器输出激光,经过单模光纤(SMF)传输,利用偏振控制器(PC)控制输出光的偏振方向,再利用准直镜A1和物镜A2将光斑准直扩束,并照射到反射式空间光调制器(SLM)上,SLM上加载计算得到叉形光栅,衍射光中±1级衍射光即为LG光。利用物镜A3将LG光聚焦后照射到衍射光阑A4上,利用电荷耦合器件(CCD)观察LG光经过光阑A4后所得到的衍射图样。

图 7. 涡旋光衍射实验。(a)示意图;(b)实物图

Fig. 7. Diffraction experiment of vortex beam. (a) Diagram; (b) physical map

下载图片 查看所有图片

图 8. 叉型光栅。(a)示意图;(b)l=1和(c) l=2时的衍射图

Fig. 8. Fork grating. (a) Schematic; diffraction diagrams when (b) l=1 and (c) l=2

下载图片 查看所有图片

SLM上加载的叉型光栅如图8(a)所示,所展示的是产生l=2的LG光的叉型光栅。图8(b)、(c)所示为叉型光栅衍射图,两张图的±1级衍射光分别为l=±1, ±2的LG光,可以看出,产生的LG光为环形光束,0级光斑为反射的高斯光。利用l=1的LG光进行衍射实验,得到的实验结果如图9所示,可以看出,衍射图样中有许多六边形环形光斑,但产生的LG光圆环上的光强均匀性较差,导致六边形亮环的光强分布不均。

图 9. l=1时的衍射光场

Fig. 9. Diffraction field under l=1

下载图片 查看所有图片

5 结论

理论分析了涡旋光经多孔的衍射情况,研究发现,LG光经正六边形排布的多孔衍射后,衍射光场会随着拓扑荷l的变化出现周期性变化。l=6k±1(k=0,1,2,…)时,会出现蜂窝状光场,其他拓扑荷的衍射会得到花瓣状的图像。对于l=6k±1的蜂窝状光场,随着拓扑荷数l的增大,中心亮环的光强逐渐减小而外侧的光强逐渐增大,且非0级亮环上的光强分布不均匀;当圆孔半径r0增大时,衍射光场的范围减小;当六边形边长r增大时,亮环的条纹宽度D及条纹间距ΔL减小,且条纹宽度与条纹间距数值相等。涡旋光束经多孔衍射后得到的光场可运用于涡旋光检测等,根据不同的衍射光场对比模板便可得到涡旋光的拓扑荷数值。

参考文献

[1] Allen L, Beijersbergen M W. Spreeuw R J C, et al. Orbital angular momentum of light and the transformation of Laguerre-Gaussian laser modes[J]. Physical Review A, 1992, 45(11): 8185-8189.

[2] Barnett S M, Allen L. Orbital angular momentum and nonparaxial light beams[J]. Optics Communications, 1994, 110(5/6): 670-678.

[3] He H. Friese M E J, Heckenberg N R, et al. Direct observation of transfer of angular momentum to absorptive particles from a laser beam with a phase singularity[J]. Physical Review Letters, 1995, 75(5): 826-829.

[4] Zhang D, Feng X, Huang Y. Encoding and decoding of orbital angular momentum for wireless optical interconnects on chip[J]. Optics Express, 2012, 20(24): 26986-26995.

[5] Lavery M P J, Speirits F C, Barnett S M, et al. . Detection of a spinning object using light’s orbital angular momentum[J]. Science, 2013, 341(6145): 537-540.

[6] Stankevicius E, Gertus T, Rutkauskas M, et al. Fabrication of micro-tube arrays in photopolymer SZ2080 by using three different methods of a direct laser polymerization technique[J]. Journal of Micromechanics and Microengineering, 2012, 22(6): 065022.

[7] Ni J, Wang C, Zhang C, et al. Three-dimensional chiral microstructures fabricated by structured optical vortices in isotropic material[J]. Light: Science & Applications, 2017, 6(7): e17011.

[8] Sztul H I, Alfano R R. Double-slit interference with Laguerre-Gaussian beams[J]. Optics Letters, 2006, 31(7): 999-1001.

[9] Ghai D P, Senthilkumaran P, Sirohi R S. Single-slit diffraction of an optical beam with phase singularity[J]. Optics and Lasers in Engineering, 2009, 47(1): 123-126.

[10] Hickmann J M. Fonseca E J S, Soares W C, et al. Unveiling a truncated optical lattice associated with a triangular aperture using light’s orbital angular momentum[J]. Physical Review Letters, 2010, 105(5): 053904.

[11] Bhattacharyya A. On a measure of divergence between two multinomial populations[J]. The Indian Journal of Statistics, 1946, 7(4): 401-406.

李润泉, 王智, 崔粲, 王慧莹, 李颖. 正六边形多孔阵列的涡旋光衍射[J]. 光学学报, 2018, 38(10): 1005002. Runquan Li, Zhi Wang, Can Cui, Huiying Wang, Ying Li. Diffraction of Vortex Beam by Regular Hexagonal Multi-Hole Array[J]. Acta Optica Sinica, 2018, 38(10): 1005002.

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

相关论文

加载中...

关于本站 Cookie 的使用提示

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