CN104777538A - 衍射光相位可控的二维达曼光栅 - Google Patents

衍射光相位可控的二维达曼光栅 Download PDF

Info

Publication number
CN104777538A
CN104777538A CN201510174137.3A CN201510174137A CN104777538A CN 104777538 A CN104777538 A CN 104777538A CN 201510174137 A CN201510174137 A CN 201510174137A CN 104777538 A CN104777538 A CN 104777538A
Authority
CN
China
Prior art keywords
diffraction
span
order
light
diffraction time
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
CN201510174137.3A
Other languages
English (en)
Other versions
CN104777538B (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.)
Shanghai Institute of Optics and Fine Mechanics of CAS
Original Assignee
Shanghai Institute of Optics and Fine Mechanics of CAS
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 Shanghai Institute of Optics and Fine Mechanics of CAS filed Critical Shanghai Institute of Optics and Fine Mechanics of CAS
Priority to CN201510174137.3A priority Critical patent/CN104777538B/zh
Publication of CN104777538A publication Critical patent/CN104777538A/zh
Application granted granted Critical
Publication of CN104777538B publication Critical patent/CN104777538B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B5/00Optical elements other than lenses
    • G02B5/18Diffraction gratings
    • G02B5/1814Diffraction gratings structurally combined with one or more further optical elements, e.g. lenses, mirrors, prisms or other diffraction gratings
    • G02B5/1819Plural gratings positioned on the same surface, e.g. array of gratings
    • G02B5/1823Plural gratings positioned on the same surface, e.g. array of gratings in an overlapping or superposed manner

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Holo Graphy (AREA)

Abstract

一种衍射光相位可控的二维达曼光栅,将一个周期分为许多相同的矩形单元,通过优化各个矩形单元的相位延迟来获得所需要的远场衍射光分布。特别的,所述单元的相位延迟是二值分布的,所述的衍射光相位可控的二维达曼光栅能使单色平面波在远场形成特定衍射级次的衍射光光强相等。本发明一方面可以大幅度提高衍射光的效率,另一方面还可以实现分束不同级次之间的二元相位控制,在相干通信和相干检测等领域具有重要的应用价值。

Description

衍射光相位可控的二维达曼光栅
技术领域
本发明涉及一种衍射光学元器件,特别是一种衍射光相位可控的二维达曼光栅。
背景技术
达曼光栅是一种能将单一光束分成若干等强度光束的衍射光学元件,由于其制作简便衍射效率高,在并行刻写、图像复制等领域发挥着重大作用。然而,相干通信及相干检测技术的发展对分束的光波提出了更高的要求,不仅要求强度相等,而且相位值也相同或呈二值分布。
专利(CN102360091B)针对上述要求提出了等位相等光强分束达曼光栅及制备方法,该技术在传统达曼光栅设计的基础上对衍射光的相位进行了约束,实现了将单一波长的入射平面波在远场分成若干等相位等光强的光束。但该技术只针对一维的达曼光栅进行设计,当分束的数目较大时,由于衍射光只在一维方向上排列,占用空间较大,使得系统结构较大,整体装置笨重。虽然可以通过叠加两个一维的光栅形成可分离二维达曼光栅以解决上述问题,但此类二维光栅的衍射效率较低,因为其效率为两个一维光栅效率的乘积。
发明内容
本发明的目的在于提供一种衍射光相位可控的达曼光栅,该达曼光栅在单色平面光波的照射下可使得二维分布的若干束等强度衍射光相位至多为2阶,极大地方便了相干通信及相干检测等技术的相位匹配过程。
本发明的技术解决方案如下:
一种衍射光相位可控的二维达曼光栅,其特点在于构成是把每一个光栅周期分为N*M个矩形单元,其中N和M均为整数;该光栅分束阵列定义为一个矩形的光束阵列分布,其中列方向的光束数目为该光栅在x方向上的分束比,行方向光束数目为该光栅在y方向上的分束比;
当分束阵列为(2L+1)*(2K+1)时,所需的衍射级次的光强I(fl,fk)均相等,其中fl表示沿x轴方向上的衍射级次,fk表示沿y轴方向上的衍射级次,I(fl,fk)表示沿x轴衍射级次为fl同时沿y轴衍射级次为fk的光强,l的取值范围为[-L,L]的整数,k的取值范围为[-K,K]的整数;
当分束阵列为(2L+1)*(2K)时,所需的衍射级次的光强I(fl,fk)均相等,l的取值范围为[-L,L]的整数,k的取值范围为[-2K+1,2K-1]的奇数;
当分束阵列为(2L)*(2K+1)时,所需的衍射级次的光强I(fl,fk)均相等,l的取值范围为[-2L+1,2L-1]的奇数,k的取值范围为[-K,K]的整数;
当分束阵列为(2L)*(2K)时,所需的衍射级次的光强I(fl,fk)均相等,l的取值范围为[-2L+1,2L-1]的奇数,k的取值范围为[-2K+1,2K-1]的奇数;
所需的衍射级次的相位为0或π,该二维光栅两个方向上的周期均已归一化,一个周期的相位延迟A(n,m)的值采用全局优化算法进行编程优化算,其中(n,m)表示矩形单元在该光栅周期中的坐标,n表示x轴上的坐标,m表示y轴上的坐标,n的取值范围为[1,N]的整数,m的取值范围为[1,M]的整数,每个单元的相位延迟A(n,m)的值二值化,为1或为exp(iθ),当该单元值为1时表示该矩形单元对入射光相位延迟为0,当该矩形单元值为exp(iθ)时,该单元对入射光相位延迟为θ;远场衍射光的复振幅表达式为:
T ( f l , f k ) = 1 NM Σ n = 1 N Σ m = 1 M A ( n , m ) · exp ( j 2 π f l n - 1 / 2 N ) · exp ( j 2 π f k m - 1 / 2 M ) · sin c ( f l N ) · sin c ( f k M )
由此得出每个衍射级次的光强为:
I(fl,fk)=T(fl,fk)·T*(fl,fk)
其中,T*(fl,fk)表示T(fl,fk)的复共轭;
每个衍射级次的相位为:
Φ = ( f l , f k ) = ar tan { Imag { T ( f l , f k ) } Real { T ( f l , f k ) } }
其中,Imag{.}表示取虚部,Real{.}表示取实部。
一个光栅周期内A(n,m)分布的不同对应着不同的远场分布,通过全局优化算法优化A(n,m)的分布可得到所需衍射级次的强度和相位分布,特别地,使所需衍射级次拥有相等的强度分布,衍射级次的相位值至多为2个。当分束比为奇数(2L+1)时,所需的衍射级次为0,±1,±2,…,±L。当分束比为偶数2L时,所需的衍射级次为±1,±3,…,±(2L-1)。
衍射光相位可控的二维达曼光栅的描述参数包括:
1.衍射光相位可控的二维达曼光栅的分束比,即在二维方向上所需的衍射级次的个数。
2.一个光栅周期被分为完全相同的单元的数目N*M。
3.对于衍射光相位可控的二维达曼光栅,衍射效率的定义为:
分束阵列为(2L+1)*(2K+1)时,
分束阵列为(2L+1)*(2K)时,
分束阵列为(2L)*(2K+1)时,
分束阵列为(2L*2K)时, Eff = Σ l = - L L - 1 Σ k = - K K - 1 I ( f 2 l + 1 , f 2 k + 1 )
4.光强均匀度
Unif Int = I max - I min I max + I min
式中Imax和Imin分别表示所需衍射级次中光强最大值和最小值。
5.相位均匀度
Unif Phs = max { abs { Φ π - round { Φ π } } · π }
式中Φ为任意一个所需衍射级次的相位值,round{.}表示四舍五入求整,abs{.}表示取绝对值,max{.}表示取最大值。
在全局优化算法中,优化A(n,m)的分布使得衍射效率尽可能大,光强均匀度参量和相位均匀度参量尽可能小。
本发明的技术效果如下:
单位单色平面波经本发明衍射光相位可控的二维达曼光栅衍射时,所需衍射级次的衍射光远场光强相等,且相位值至多只有2个,极大地方便了相干通信和相干检测中的相位匹配过程,同时由于采用的光栅结构,使得光栅的衍射效率相对于分离二维达曼光栅的衍射效率高。
附图说明
图1是分束阵列为4*4的衍射光相位可控的二维达曼光栅一个周期的示意图。
图2是分束阵列为4*4的衍射光相位可控的二维达曼光栅结构示意图。
图3是分束阵列为5*5的衍射光相位可控的二维达曼光栅一个周期的示意图。
图4是分束阵列为5*5的衍射光相位可控的二维达曼光栅结构示意图。
具体实施方式
下面结合附图对和实施例进一步详细说明发明,但不应以此限制本发明的保护范围。
把光栅的一个周期分为N*M个完全相同的矩形单元,每个矩形单元的坐标用(n,m)表示,其中1≤n≤N,1≤m≤M。每个矩形单元对入射光的相对相位延迟为0和π,对应A(n,m)的取值分别为1和-1。采用模拟退火算法优化A(n,m)的分布,使得目标函数取值最低,目标函数可表示为,
costF=1/Eff+UnifInt+UnifPhs
本发明实例给出了分束阵列为4*4和5*5的衍射光相位可控的二维达曼光栅,其优化结果如图1和图3所示。
图1是分束阵列为4*4的衍射光相位可控的二维达曼光栅一个周期的示意图。所述的一个周期由128*128个完全相同的矩形单元组成,示意图中用一个像素代表一个矩形单元,像素的颜色代表不同的相对相位延迟,黑色代表相位延迟为0,白色代表相位延迟为π。
图2是分束阵列为4*4的衍射光相位可控的二维达曼光栅结构示意图。该结构由10*10个光栅周期组成,在实际应用中将会由更多的周期组成。单位单色光垂直入射到实施例2的光栅结构的远场衍射光的强度分布见表1。单位单色光垂直入射到实施例2的光栅结构的远场衍射光的相位分布见表2。表格列方向的衍射级次表示x方向的衍射级次,表格行方向的衍射级次表示y方向的衍射级次。光栅衍射效率为72.5%,光强均匀度为0.0108,相位均匀度为0.0108。
表1.分束阵列4*4的二维达曼光栅远场衍射光的光强分布
衍射级次 -3 -1 +1 +3
-3 0.045738 0.044934 0.045308 0.044759
-1 0.045494 0.045449 0.045465 0.045296
+1 0.045296 0.045465 0.045449 0.045494
+3 0.044759 0.045308 0.044934 0.045738
表2.分束阵列4*4的二维达曼光栅远场衍射光的相位分布
衍射级次 -3 -1 +1 +3
-3 3.147085 0.0011656 3.151785 3.145185
-1 0.01077 0.00037042 3.1346 3.1377
+1 3.145485 3.148585 0.00037042 0.01077
+3 3.138 3.1314 0.0011656 3.1361
图3是分束阵列为5*5的衍射光相位可控的二维达曼光栅一个周期的示意图。所述的一个周期由64*64的完全相同的矩形单元组成,示意图中用一个像素代表一个矩形单元,像素的颜色代表不同的相对相位延迟,黑色代表相位延迟为0,白色代表相位延迟为π。
图4为本发明实施例2的光栅结构示意图,该结构由10*10个光栅周期组成,在实际应用中将会由更多的周期组成。单位单色光垂直入射到实施例2的光栅结构的远场衍射光的强度分布见表3。单位单色光垂直入射到实施例2的光栅结构的远场衍射光的相位分布见表4。表格列方向的衍射级次表示x方向的衍射级次,表格行方向的衍射级次表示y方向的衍射级次。光栅衍射效率为73.4%,光强均匀度为0.0088,相位均匀度为0.0160。
表3.实施例2远场衍射光的光强分布
衍射级次 -2 -1 0 1 2
-2 0.0295 0.0296 0.0295 0.0292 0.0294
-1 0.0295 0.0294 0.0291 0.0296 0.0291
0 0.0291 0.0293 0.0294 0.0293 0.0291
+1 0.0291 0.0296 0.0291 0.0294 0.0295
+2 0.0294 0.0292 0.0295 0.0296 0.0295
表4.实施例2远场衍射光的相位分布
衍射级次 -2 -1 0 1 2
-2 0.0112 3.1374 0.0123 0.0047 0.016
-1 0.0135 3.156785 3.149185 0.0062 0.0136
0 3.1275 3.1383 0 3.144885 3.155685
+1 0.0136 0.0062 3.134 3.1264 0.0135
+2 0.016 0.0047 0.0123 3.145785 0.0112

Claims (1)

1.一种衍射光相位可控的二维达曼光栅,其特征在于构成是把每一个光栅周期分为N*M个矩形单元,其中N和M均为整数;该光栅分束阵列定义为一个矩形的光束阵列分布,其中列方向的光束数目为该光栅在x方向上的分束比,行方向光束数目为该光栅在y方向上的分束比;
当分束阵列为(2L+1)*(2K+1)时,所需的衍射级次的光强I(fl,fk)均相等,其中fl表示沿x轴方向上的衍射级次,fk表示沿y轴方向上的衍射级次,I(fl,fk)表示沿x轴衍射级次为fl同时沿y轴衍射级次为fk的光强,l的取值范围为[-L,L]的整数,k的取值范围为[-K,K]的整数;
当分束阵列为(2L+1)*(2K)时,所需的衍射级次的光强I(fl,fk)均相等,l的取值范围为[-L,L]的整数,k的取值范围为[-2K+1,2K-1]的奇数;
当分束阵列为(2L)*(2K+1)时,所需的衍射级次的光强I(fl,fk)均相等,l的取值范围为[-2L+1,2L-1]的奇数,k的取值范围为[-K,K]的整数;
当分束阵列为(2L)*(2K)时,所需的衍射级次的光强I(fl,fk)均相等,l的取值范围为[-2L+1,2L-1]的奇数,k的取值范围为[-2K+1,2K-1]的奇数;
所需的衍射级次的相位为0或π,该二维光栅两个方向上的周期均已归一化,一个周期的相位延迟A(n,m)的值采用全局优化算法进行编程优化算,其中(n,m)表示单元在该光栅周期中的坐标,n表示x轴上的坐标,m表示y轴上的坐标,n的取值范围为[1,N]的整数,m的取值范围为[1,M]的整数,每个单元的相位延迟A(n,m)的值二值化,为1或为exp(iθ),当该单元值为1时表示该单元对入射光相位延迟为0,当该单元值为exp(iθ)时,该单元对入射光相位延迟为θ;远场衍射光的复振幅表达式为:
T ( f l , f k ) = 1 NM Σ n = 1 N Σ m = 1 M A ( n , m ) · exp ( j 2 π f l m - 1 / 2 N ) · exp ( j 2 π f k m - 1 / 2 M ) · sin c ( f l N ) · sin c ( f k M )
由此得出每个衍射级次的光强为:
I(fl,fk)=T(fl,fk)·T*(fl,fk)
其中,T*(fl,fk)表示T(fl,fk)的复共轭;
每个衍射级次的相位为:
Φ ( f l , f k ) = arctan { Imag { T ( f l , f k ) } Real { T ( f l , f k ) } }
其中,Imag{.}表示取虚部,Real{.}表示取实部。
CN201510174137.3A 2015-04-14 2015-04-14 衍射光相位可控的二维达曼光栅 Active CN104777538B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510174137.3A CN104777538B (zh) 2015-04-14 2015-04-14 衍射光相位可控的二维达曼光栅

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510174137.3A CN104777538B (zh) 2015-04-14 2015-04-14 衍射光相位可控的二维达曼光栅

Publications (2)

Publication Number Publication Date
CN104777538A true CN104777538A (zh) 2015-07-15
CN104777538B CN104777538B (zh) 2018-03-06

Family

ID=53619114

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510174137.3A Active CN104777538B (zh) 2015-04-14 2015-04-14 衍射光相位可控的二维达曼光栅

Country Status (1)

Country Link
CN (1) CN104777538B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106990529A (zh) * 2017-06-05 2017-07-28 南开大学 一种计算一维反射光栅远场分布的标量衍射理论算法
CN107490869A (zh) * 2017-08-24 2017-12-19 华天科技(昆山)电子有限公司 空间结构光发射装置
CN112558314A (zh) * 2020-12-01 2021-03-26 中国空间技术研究院 激光阵列分束器和光栅设计方法
CN112859215A (zh) * 2020-12-31 2021-05-28 中国科学院光电技术研究所 一种红外波段的准连续的超表面分束器件
CN112880844A (zh) * 2021-01-18 2021-06-01 中国科学院上海光学精密机械研究所 一种基于分束元件的单次曝光复振幅测量装置和方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003014915A (ja) * 2001-07-03 2003-01-15 Japan Science & Technology Corp Dammann型グレーティングをつけた光学素子
JP2003329822A (ja) * 2002-05-14 2003-11-19 Japan Science & Technology Corp 集光型光スプリッタ及びその作製方法
CN102360091A (zh) * 2011-09-28 2012-02-22 中国科学院上海光学精密机械研究所 等位相等光强分束达曼光栅及其制备方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003014915A (ja) * 2001-07-03 2003-01-15 Japan Science & Technology Corp Dammann型グレーティングをつけた光学素子
JP2003329822A (ja) * 2002-05-14 2003-11-19 Japan Science & Technology Corp 集光型光スプリッタ及びその作製方法
CN102360091A (zh) * 2011-09-28 2012-02-22 中国科学院上海光学精密机械研究所 等位相等光强分束达曼光栅及其制备方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
H.DAMMANN.ETC: "HIGH-EFFICIENCY IN-LINE MULTIPLE IMAGING BY MEANS OF MULTIPLE PHASE HOLOGRAMS", 《OPTICS COMMUNICATIONS》 *
谷雨: "对达曼光栅进行二维扩展的矩形孔径光栅设计", 《激光杂志》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106990529A (zh) * 2017-06-05 2017-07-28 南开大学 一种计算一维反射光栅远场分布的标量衍射理论算法
CN106990529B (zh) * 2017-06-05 2021-12-10 南开大学 一种计算一维反射光栅远场分布的标量衍射理论算法
CN107490869A (zh) * 2017-08-24 2017-12-19 华天科技(昆山)电子有限公司 空间结构光发射装置
CN112558314A (zh) * 2020-12-01 2021-03-26 中国空间技术研究院 激光阵列分束器和光栅设计方法
CN112859215A (zh) * 2020-12-31 2021-05-28 中国科学院光电技术研究所 一种红外波段的准连续的超表面分束器件
CN112880844A (zh) * 2021-01-18 2021-06-01 中国科学院上海光学精密机械研究所 一种基于分束元件的单次曝光复振幅测量装置和方法

Also Published As

Publication number Publication date
CN104777538B (zh) 2018-03-06

Similar Documents

Publication Publication Date Title
CN104777538A (zh) 衍射光相位可控的二维达曼光栅
Francesca et al. Three-term commutator estimates and the regularity of 1/2-harmonic maps into spheres
CN101701847B (zh) 基于光栅和ccd成像探测器的宽动态范围成像系统
CN104375277B (zh) 基于改进型迈克尔逊干涉仪的多通道空间结构光场生成装置
CN105490033A (zh) 一种l形子阵运用方法
CN105589203A (zh) 产生径向偏振阵列光束的方法及装置
CN101592783B (zh) 相干阵列激光逆达曼光栅合束孔径装填装置
CN102360091B (zh) 等位相等光强分束达曼光栅及其制备方法
Liang et al. Optimizing illumination’s complex coherence state for overcoming Rayleigh’s resolution limit
CN101187709A (zh) 消除零级衍射谱点的等光强分束光栅
Zeng et al. Three-dimensional vectorial multifocal arrays created by pseudo-period encoding
CN105242413A (zh) 一种六角阵列螺旋相位板及制作方法
Liu et al. Discrete emitters as a source of orbital angular momentum
CN108333784A (zh) 一种基于双光栅衍射产生规则涡旋光阵列的方法及系统
Han et al. Interaction structures of multi localized waves within the Kadomtsev–Petviashvili I equation
Hassan et al. Efficient beamforming and radiation pattern control using stacked intelligent metasurfaces
Abedin et al. The iterative Fresnel integrals method for Fresnel diffraction from tilted rectangular apertures: Theory and simulations
CN100468205C (zh) 一种消色差光变图像的制作方法
CN105372623A (zh) 一种基于l型阵列的信源仰角和方位角估计方法
CN102122089A (zh) 基于液晶光学相控阵的光束焦点三维独立控制的方法
Abedin et al. Computer simulation of Fresnel diffraction from double rectangular apertures in one and two dimensions using the iterative Fresnel integrals method
CN114355743B (zh) 一种多衍射级独立光波场调控的全息编码方法
Sabatyan et al. Considerable diversity in generation light-arm beams using multi-twisted phase structure zone plate
CN115308914A (zh) 基于渐变费马螺旋缝产生涡旋光束的系统及方法
Cruz et al. On achieving a desired flux distribution on the receiver of a solar power tower plant

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant