CN113504642A - 一种具有多中空的紧聚焦光场分布构造方法 - Google Patents

一种具有多中空的紧聚焦光场分布构造方法 Download PDF

Info

Publication number
CN113504642A
CN113504642A CN202110643219.3A CN202110643219A CN113504642A CN 113504642 A CN113504642 A CN 113504642A CN 202110643219 A CN202110643219 A CN 202110643219A CN 113504642 A CN113504642 A CN 113504642A
Authority
CN
China
Prior art keywords
matrix
light field
vortex
distribution
calculating
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110643219.3A
Other languages
English (en)
Other versions
CN113504642B (zh
Inventor
季开来
覃亚丽
徐梦华
刘小旭
郑欢
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University of Technology ZJUT
Original Assignee
Zhejiang University of Technology ZJUT
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang University of Technology ZJUT filed Critical Zhejiang University of Technology ZJUT
Priority to CN202110643219.3A priority Critical patent/CN113504642B/zh
Publication of CN113504642A publication Critical patent/CN113504642A/zh
Application granted granted Critical
Publication of CN113504642B publication Critical patent/CN113504642B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B27/00Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
    • G02B27/0012Optical design, e.g. procedures, algorithms, optimisation routines
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Optics & Photonics (AREA)
  • Polarising Elements (AREA)

Abstract

一种具有多中空的紧聚焦光场分布构造方法,包括以下步骤:步骤1,在紧聚焦系统中,设置入射平面的初始参数;步骤2,设计包含多涡旋奇点的相位分布矩阵P,对涡旋的数量为S、第n个涡旋的位置为和携带的拓扑电荷大小为ln进行调控;步骤3,以径向偏振与角向偏振矢量叠加为混合偏振状态;步骤4,将HPBGVB的透射光场矩阵Et,通过Deybe‑Wolf衍射积分方程得到聚焦光场矩阵Eout;步骤5,根据所得的聚焦光场矩阵Eout计算出聚焦光强矩阵。本发明可以灵活地调控紧聚焦光场分布中空的数量。

Description

一种具有多中空的紧聚焦光场分布构造方法
技术领域
本发明涉及一种具有多中空的紧聚焦光场分布构造方法,属于焦场调控的技术领域。
背景技术
焦场调控和构造以满足改善光学操纵、粒子捕获等领域的需要以及性能的提升,这是当前的重要研究课题。它指的是使用透镜在焦点附近(焦平面上)获得聚焦光场,可以在入射端对相关的参数(例如偏振形式、拓扑电荷等)调整,控制焦场的强度以及相位分布满足特定的需要。一个毫米级别的高斯光束通过高数值孔径(NA>0.7)物镜,产生一个大小可达到波长量级的光斑,可以有效地提高聚焦成像系统的空间分辨率。相较于普通聚焦,光束的紧聚焦会在聚焦中心产生一个折射率远小于背景折射率的能量分布,这将有利于粒子的捕获。因此光束的紧聚焦特性受到许多研究者的关注。
涡旋光束可以简单地看成是拉盖尔高斯模式的线性叠加。涡旋光束又称局域空心光束,在相位奇点处具有中心凹陷的强度分布特点。在1992年,Allen等人验证了涡旋光束具有相位并且每个光子携带携带轨道角动量(Orbital angular momentum,OAM)为,其中为拓扑电荷(Topological charge,TC)。具有OAM的光束会围绕相位奇点旋转沿着光轴向前传播。涡旋光束携带的TC幅值和符号,可以通过与球面波干涉产生的明亮条纹数量和旋转方向进行检测。利用涡旋光束的空间维度信息,可以极大提高通信焦场调控的自由度。
偏振是光束的一项重要特性,它描述的是光束在横向平面内的电场强度的振动方向。相较于常见的线性偏振和椭圆偏振属于均匀偏振,径向偏振和角向偏振属于非均匀偏振。其中径向偏振描述光场的横向分量的振动方向总是和径矢方向保持一致,而角向偏振的振动方向总是与径矢垂直。近些年,由于径向偏振光束和角向偏振光束在光纤通信、光学操纵、粒子捕获、光学数据存储、聚焦显微镜等领域宽泛的应用,它们已经被广泛的研究。
将涡旋特性加在矢量光束中,再经高数值孔径聚焦后,会有许多有趣的现象生成。当TC为零时,径向偏振光束紧聚焦后会产生一个很小的光斑,在焦点附近表现出一个很强的纵向光场分量;角向偏振光束在焦点附近表现为中心凹陷的圆环。当TC为1时,径向偏振和角向偏振光束都表现为一个光斑。但是这种聚焦程度有所不同,相较于径向偏振,角向偏振光束所得的光斑的半径要更小。此外,当TC为非整数时,矢量涡旋光束表现为更多聚焦特性,形成具有开口的涡旋光束。由此可见,偏振状态和涡旋特性对光束的聚焦产生巨大的影响。
发明内容
为了解决在粒子捕获中所需要提供多个中心暗洞的光笼结构的焦场调控问题,本发明提供一种具有多中空的紧聚焦光场分布构造方法,采用携带多拓扑电荷的双涡旋相位的方式,可以灵活地调控紧聚焦光场分布中空的数量。
本发明解决其技术问题所采用的技术方案是:
一种具有多中空的紧聚焦光场分布构造方法,包括以下步骤:
步骤1,在紧聚焦系统中,设置入射平面的初始参数:入射平面先采用在直角坐标系下N×N的矩阵X和Y,上下边界分别为R_xy和-R_xy,点间距为2R_xy/(N-1),计算出矩阵X和Y在极坐标下的极矢量长度矩阵ρ和角矩阵
Figure BDA0003107920560000021
根据透镜的数值孔径NA与光束折射角θ的关系,计算出折射角矩阵θ:
θ=arcsin[ρ·NA/(nt·R)]
其中,nt为透镜的折射率,R为透镜的最大半径;
根据标量波矢k=2π/λ,折射角θ和角矩阵
Figure BDA0003107920560000022
计算出波矢矩阵k=(kx,ky,kz):
Figure BDA0003107920560000023
再依据横向波矢矩阵
Figure BDA0003107920560000024
计算出横向波矢矩阵kρ
步骤2,设计包含多涡旋奇点的相位分布矩阵P,对涡旋的数量为S、第n个涡旋的位置为和携带的拓扑电荷大小为ln进行调控,涡旋相位分布矩阵P计算表达式为:
Figure BDA0003107920560000025
步骤3,以径向偏振与角向偏振矢量叠加为混合偏振状态,设置代表入射光束混合偏振状态的混合偏振角φ0,计算出混合偏振贝塞尔高斯涡旋光束(Hybrid PolarizedBessel Gaussian Vortex Beam,HPBGVB)在透射面的混合偏振状态矩阵V:
Figure BDA0003107920560000026
根据极矢量矩阵ρ、横向波矢矩阵kρ、设计的涡旋相位分布矩阵P和混合偏振状态矩阵V,计算出HPBGVB的透射光场矩阵Et
Figure BDA0003107920560000027
其中,入射光束携带的等效拓扑电荷数量为
Figure BDA0003107920560000031
w0为高斯腰宽,最大峰值1/e强度对应的半宽,HPBGVB的入射光束光场矩阵为
Figure BDA0003107920560000032
步骤4,将HPBGVB的透射光场矩阵Et,通过Deybe-Wolf衍射积分方程得到聚焦光场矩阵Eout
Figure BDA0003107920560000033
其中,r和Φ分别为焦平面的极矢量长度矩阵和角矩阵,α为透射光束的最大折射角,z为光束聚焦的传输距离;
步骤5,根据所得的聚焦光场矩阵Eout计算出聚焦光强矩阵I:
I=|Ex|2+|Ey|2+|Ez|2
矩阵I即为最后获得的具有多中空的紧聚焦光场分布,通过设计不同涡旋相位分布矩阵P,获得不同中空数量的紧聚焦光场分布。
鉴于入射光束的偏振状态和附加涡旋对聚焦特性的影响,本发明中选择贝塞尔高斯涡旋光束作为标量场,将每个涡旋引入位置信息,将径向偏振和角向偏振进行矢量组合,成为HPBGVB。紧聚焦光强分布中涡旋相位奇点所在区域常常表现为暗洞的结构。因此,可以通过采用携带多拓扑电荷的双离轴涡旋相位的方式,来调控聚焦光强中心暗洞的数量的生成。这些具有多中空的焦场结构可以同时捕获多个粒子分别在这些中空暗洞中,将在粒子捕获、光学操纵等领域具有一定的研究意义。此外,本发明中提出的混合偏振比例、涡旋的数量、涡旋位置的坐标信息,将为在焦平面上获得特定光强分布的聚焦光束提供更多的调制自由度提供可能。
本发明的有益效果主要表现在:相较于复梯度相位方法相比,本发明使用携带相同拓扑电荷数量关于光轴对称的双涡旋相位方法,这种加载多个离轴涡旋的方法将更方便地在焦平面上获得具有多中空的紧聚焦光场分布。获得焦场结构的中空数量,可以通过两个涡旋携带的拓扑电荷的大小进行准确调控,有利于提高中空焦场结构制备的效率。而且获得的紧聚焦光场分布的中空数量包含奇数和偶数,更具有普遍性。通过这种使用涡旋位置和携带拓扑电荷数量来设计涡旋相位分布的方式,可以有效降低涡旋相位设计的复杂度,提高数据处理效率。获得的多中空紧聚焦光场分布,可以实现同时捕获多个粒子的实际应用。
附图说明
图1为本发明中获得具有多中空的紧聚焦光场分布的流程图。
图2为本发明中获得具有多中空的紧聚焦光光场分布的紧聚焦系统框图。
图3为本发明中设计的多拓扑电荷的双离轴涡旋相位分布。
图4为本发明中获得具有多中空的紧聚焦光场分布结果。
具体实施方式
下面结合附图对本发明作进一步描述。
参考图1~图4,一种具有多中空的紧聚焦光场分布构造方法,包括以下步骤:
步骤1,在紧聚焦系统中,设置入射平面的初始参数:入射平面先采用在直角坐标系下N×N的矩阵X和Y,上下边界分别为R_xy和-R_xy,点间距为2R_xy/(N-1),计算出矩阵X和Y在极坐标下的极矢量长度矩阵ρ和角矩阵
Figure BDA0003107920560000041
根据透镜的数值孔径NA与光束折射角θ的关系,计算出折射角矩阵θ:
θ=arcsin[ρ·NA/(nt·R)]
其中,nt为透镜的折射率,R为透镜的最大半径;
根据标量波矢k=2π/λ,折射角θ和角矩阵
Figure BDA0003107920560000042
计算出波矢矩阵k=(kx,ky,kz):
Figure BDA0003107920560000043
再依据横向波矢矩阵
Figure BDA0003107920560000044
计算出横向波矢矩阵kρ
步骤2,设计包含多涡旋奇点的相位分布矩阵P,对涡旋的数量为S、第n个涡旋的位置为和携带的拓扑电荷大小为ln进行调控,涡旋相位分布矩阵P计算表达式为:
Figure BDA0003107920560000045
由于相位奇点所在的位置,光强分布展现为暗洞的结果。基于这个原理,可以将多个涡旋相位奇点在不同的位置上进行排列,将产生一些特定暗洞形状分布的紧聚焦光场分布。为了制备具有多中空的中空紧聚焦光场分布,本发明设计的多拓扑电荷的双离轴涡旋相位分布,两个涡旋关于原点(光轴)对称,且携带的拓扑电荷相同l1=l2=ln。通过调控每个涡旋的携带的拓扑电荷数量ln和涡旋位置的极矢量大小ρk,设计的多拓扑电荷的双离轴涡旋相位分布以及对应的入射光束光强分布(右上角),如图3所示;
步骤3,以径向偏振与角向偏振矢量叠加为混合偏振状态,设置代表入射光束混合偏振状态的混合偏振角φ0,计算出HPBGVB在透射面的混合偏振状态矩阵V:
Figure BDA0003107920560000051
根据极矢量矩阵ρ、横向波矢矩阵kρ、设计的涡旋相位分布矩阵P和混合偏振状态矩阵V,计算出HPBGVB的透射光场矩阵Et
Figure BDA0003107920560000052
其中,入射光束携带的等效拓扑电荷数量为
Figure BDA0003107920560000053
w0为高斯腰宽(最大峰值1/e强度对应的半宽),HPBGVB的入射光束光场矩阵为
Figure BDA0003107920560000054
步骤4,将HPBGVB的透射光场矩阵Et,通过Deybe-Wolf衍射积分方程得到聚焦光场矩阵Eout
Figure BDA0003107920560000055
其中,r和Φ分别为焦平面的极矢量长度矩阵和角矩阵,α为透射光束的最大折射角,z为光束聚焦的传输距离;
步骤5,根据所得的聚焦光场矩阵Eout计算出聚焦光强矩阵I:
I=|Ex|2+|Ey|2+|Ez|2
制备获得的多中空焦场结构的光强分布结果,如图4所示。图4所示了紧聚焦光强分布的中心凹陷区域中空数量分别为二、三、四和五的结果。为了更方便地获取具有多中空的紧聚焦光场分布,最大的中空数量与加载的双离轴涡旋中每个涡旋携带的具有相同符号拓扑电荷的大小有关。设定最大的中空数量为Sp,可以表述为:最大的中空数量=每个涡旋携带的拓扑电荷大小之和-1(Sp=l1+l2-1)。本发明中所提到的紧聚焦光场分布的中空数量可以通过在入射平面对所加载的双离轴涡旋相位携带的拓扑电荷大小进行调控,这将为调控紧聚焦光场分布的中空数量提供便利。
上述实施例仅例示性说明本发明的原理及其功效,而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下,对上述实施例进行修饰或改变。因此,举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变,仍应由本发明的权利要求所涵盖。

Claims (1)

1.一种具有多中空的紧聚焦光场分布构造方法,其特征在于,所述方法包括以下步骤:
步骤1,在紧聚焦系统中,设置入射平面的初始参数:入射平面先采用在直角坐标系下N×N的矩阵X和Y,上下边界分别为R_xy和-R_xy,点间距为2R_xy/(N-1),计算出矩阵X和Y在极坐标下的极矢量长度矩阵ρ和角矩阵
Figure FDA0003107920550000019
根据透镜的数值孔径NA与光束折射角θ的关系,计算出折射角矩阵θ:
θ=arcsin[ρ·NA/(nt·R)]
其中,nt为透镜的折射率,R为透镜的最大半径;
根据标量波矢k=2π/λ,折射角θ和角矩阵
Figure FDA00031079205500000110
计算出波矢矩阵k=(kx,ky,kz):
Figure FDA0003107920550000011
再依据横向波矢矩阵
Figure FDA0003107920550000012
计算出横向波矢矩阵kρ
步骤2,设计包含多涡旋奇点的相位分布矩阵P,对涡旋的数量为S、第n个涡旋的位置为和携带的拓扑电荷大小为ln进行调控,涡旋相位分布矩阵P计算表达式为:
Figure FDA0003107920550000013
步骤3,以径向偏振与角向偏振矢量叠加为混合偏振状态,设置代表入射光束混合偏振状态的混合偏振角φ0,计算出HPBGVB在透射面的混合偏振状态矩阵V:
Figure FDA0003107920550000014
根据极矢量矩阵ρ、横向波矢矩阵kρ、设计的涡旋相位分布矩阵P和混合偏振状态矩阵V,计算出HPBGVB的透射光场矩阵Et
Figure FDA0003107920550000015
其中,入射光束携带的等效拓扑电荷数量为
Figure FDA0003107920550000016
w0为高斯腰宽,最大峰值1/e强度对应的半宽,HPBGVB的入射光束光场矩阵为
Figure FDA0003107920550000017
步骤4,将HPBGVB的透射光场矩阵Et,通过Deybe-Wolf衍射积分方程得到聚焦光场矩阵Eout
Figure FDA0003107920550000018
其中,r和Φ分别为焦平面的极矢量长度矩阵和角矩阵,α为透射光束的最大折射角,z为光束聚焦的传输距离;
步骤5,根据所得的聚焦光场矩阵Eout计算出聚焦光强矩阵I:
I=|Ex|2+|Ey|2+|Ez|2
矩阵I即为最后获得的具有多中空的紧聚焦光场分布,通过设计不同涡旋相位分布矩阵P,获得不同中空数量的紧聚焦光场分布。
CN202110643219.3A 2021-06-09 2021-06-09 一种具有多中空的紧聚焦光场分布构造方法 Active CN113504642B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110643219.3A CN113504642B (zh) 2021-06-09 2021-06-09 一种具有多中空的紧聚焦光场分布构造方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110643219.3A CN113504642B (zh) 2021-06-09 2021-06-09 一种具有多中空的紧聚焦光场分布构造方法

Publications (2)

Publication Number Publication Date
CN113504642A true CN113504642A (zh) 2021-10-15
CN113504642B CN113504642B (zh) 2022-12-20

Family

ID=78009322

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110643219.3A Active CN113504642B (zh) 2021-06-09 2021-06-09 一种具有多中空的紧聚焦光场分布构造方法

Country Status (1)

Country Link
CN (1) CN113504642B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115047620A (zh) * 2022-07-08 2022-09-13 上海理工大学 一种在强聚焦条件下生成类超环偏振拓扑结构时空波包的方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050195480A1 (en) * 2004-01-12 2005-09-08 Brown Thomas G. System and method converting the polarization state of an optical beam into an inhomogeneously polarized state
US20100128268A1 (en) * 2007-05-02 2010-05-27 John Christopher Dainty A vectorial polarimetry method and apparatus for analyzing the three-dimensional electromagnetic field resulting from an interaction between a focused illuminating field and a sample to be observed
CN105589203A (zh) * 2016-03-13 2016-05-18 南京理工大学 产生径向偏振阵列光束的方法及装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050195480A1 (en) * 2004-01-12 2005-09-08 Brown Thomas G. System and method converting the polarization state of an optical beam into an inhomogeneously polarized state
US20100128268A1 (en) * 2007-05-02 2010-05-27 John Christopher Dainty A vectorial polarimetry method and apparatus for analyzing the three-dimensional electromagnetic field resulting from an interaction between a focused illuminating field and a sample to be observed
CN105589203A (zh) * 2016-03-13 2016-05-18 南京理工大学 产生径向偏振阵列光束的方法及装置

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
CHEN JIAN-NONG等: "Tight focus of a radially polarized and amplitudemodulated annular multi-Gaussiian beam", 《CHINESE PHYSICS B》 *
GILAD M. LERMAN, LIRON STERN AND URIEL LEVY: "Generation and Tight Focusing of Hybridly Polarized", 《OPTICAL SOCIETY OF AMERICA》 *
于娜娜等: "多涡旋相位径向矢量光束紧聚焦产生多聚焦光斑", 《河北工程大学学报》 *
曹重阳等: "紧聚焦角向偏振分数阶涡旋光诱导磁化场特性", 《物理学报》 *
李新忠等: "分数阶高阶贝塞尔涡旋光束的矢量波分析法研究", 《中国激光》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115047620A (zh) * 2022-07-08 2022-09-13 上海理工大学 一种在强聚焦条件下生成类超环偏振拓扑结构时空波包的方法
CN115047620B (zh) * 2022-07-08 2023-06-27 上海理工大学 一种在强聚焦条件下生成类超环偏振拓扑结构时空波包的方法

Also Published As

Publication number Publication date
CN113504642B (zh) 2022-12-20

Similar Documents

Publication Publication Date Title
Herrmann Optical theory of thermal velocity effects in cylindrical electron beams
CN106933027B (zh) 一种涡旋数目可控的环形涡旋阵列掩模板的设计方法
CN112147721B (zh) 偏振阶数可调且可连续变焦的柱矢量光束透镜及构造方法
CN105870604A (zh) 一种基于相位梯度超表面产生微波轨道角动量的阵列天线
CN105629463A (zh) 一种基于人工微结构超表面的圆偏振光分离器的设计方法
CN112946882B (zh) 一种空间指向和长度可定制的光管焦场产生方法
CN115359944B (zh) 一种实现任意空间指向光链焦场的方法
CN113759575A (zh) 一种调控焦场时空波包轨道角动量旋转轴指向的方法
CN115291386B (zh) 一种实现同一光焦斑任意指向直线阵列的方法
CN113504642B (zh) 一种具有多中空的紧聚焦光场分布构造方法
CN111722398B (zh) 一种强聚焦条件下生成亚波长时空涡旋的方法
Sheppard Electromagnetic field in the focal region of wide-angular annular lens and mirror systems
CN109683239B (zh) 光子集成芯片内矢量涡旋光束辐射器及其应用
CN112130231A (zh) 产生偏振阶数可调的柱矢量光束的超表面系统及构造方法
CN102967928B (zh) 一种柱偏振矢量光束紧聚焦光斑的产生方法及装置
CN112558297B (zh) 基于均匀线源辐射场产生任意空间指向光针焦场的方法
CN116165792B (zh) 一种产生携带任意指向轨道角动量涡旋焦场的方法
CN112946881B (zh) 一种产生任意指向光针立体阵的方法
CN112946883A (zh) 一种产生横向光针平面阵的方法
WO2023103112A1 (zh) 一种偏振传输不变光场的产生系统及方法
CN109375368A (zh) 一种基于空间偶极子阵的三维多焦斑阵列的产生方法
Song et al. Diversiform propagation of hollow Gaussian beam clusters in non-linear media with self-induced parabolic potentials
CN106896615A (zh) 非线性螺旋相位器件
Budhu et al. Synthesis of 3D-printed dielectric lens antennas via optimization of geometrical optics ray tracing
CN116661134B (zh) 基于超透镜产生聚焦杂化偏振矢量光束的光链生成方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant