-
本文中提出的非周期三角HCG如图 1所示。它由硅(Si)材料和空气构成,Si的折射率为3.48,空气折射率为1。在图中, Λn和Sn(n=0,1,2…)分别为非周期三角HCG光栅条的周期和底边长,非周期三角HCG光栅条占空比ηn=Sn/Λn(n=0, 1, 2…),d=1μm为非周期三角光栅条的高度,xn和ψn(n=0,1,2…)分别为某一个非周期三角HCG光栅条的中心位置和对应的相位,w是非周期三角HCG整体宽度。在实际器件制备中,可以先在Si层上生长一层牺牲层,然后利用电子束光刻和干法刻蚀等工艺程序将Si层刻蚀成非周期三角HCG,最后通过腐蚀液将牺牲层腐蚀掉,这就可以形成悬空、由空气层包围的高对比度非周期三角HCG[13]。
对于矩形HCG,通常是基于等效介质理论直接将其等效为一层介质薄膜,该介质薄膜的等效折射率由矩形光栅的参量和入射光偏振所决定[14]。但是对于三角HCG,无法利用等效介质原理直接将其等效为一层薄膜,因为其宽度由下而上逐渐变窄,可以将三角HCG等效成由多个不同占空比矩形HCG叠加而成,这就可以将整个三角HCG等效为多层等效折射率不同的薄膜,具体如图 2所示。通过计算多个矩形HCG的折射率,即可得到三角HCG的等效折射率。
矩形HCG参量和不同偏振态对其等效折射率的影响,可由下式表示[15-16]:
$ n_{\mathrm{TE}, 1}=\left[(1-\eta) n_{1}{}^{2}+\eta n_{\mathrm{h}}{}^{2}\right]^{\frac{1}{2}} $
(1) $ n_{\mathrm{TM}, 1}=\left[\frac{1-\eta}{n_{1}{}^{2}}+\frac{\eta}{n_{\mathrm{h}}{}^{2}}\right]^{-\frac{1}{2}} $
(2) $ \begin{array}{c} {n_{{\rm{TE}},2}} = \left[ {{n_{{\rm{TE}},1}}^2 + \frac{{{{\rm{ \mathsf{ π} }}^2}}}{3}{{\left( {\frac{\mathit{\Lambda }}{\lambda }} \right)}^2}{\eta ^2}(1 - } \right.\\ {\left. {\eta {)^2}{{\left( {{n_{\rm{h}}}^2 - {n_1}^2} \right)}^2}} \right]^{\frac{1}{2}}} \end{array} $
(3) $ \begin{array}{*{20}{c}} {{n_{{\rm{TM}},2}} = \left[ {{n_{{\rm{TM}},1}}^2 + \frac{{{{\rm{ \mathsf{ π} }}^2}}}{3}{{\left( {\frac{\mathit{\Lambda }}{\lambda }} \right)}^2}{\eta ^2}(1 - } \right.}\\ {{{\left. {\eta {)^2}{{\left( {\frac{1}{{{n_{\rm{h}}}^2}} - \frac{1}{{{n_1}^2}}} \right)}^2}{n_{{\rm{TE}},1}}^6{n_{{\rm{TM}},1}}^2} \right]}^{\frac{1}{2}}}} \end{array} $
(4) 式中, nl和nh分别是光栅的低折射率材料和高折射率材料,η是光栅占空比(光栅条宽/光栅周期),Λ是光栅周期,nTE, 1和nTM, 1分别是不同偏振下光栅等效成薄膜的1阶等效折射率。从(1)式和(2)式可知,在1阶等效中,光栅等效薄膜折射率只与光栅占空比和光栅材料有关,nTE, 2和nTM, 2分别是不同偏振下光栅等效成薄膜的2阶等效折射率,从(3)式和(4)式中可以发现,在2阶等效中,薄膜折射率还与入射光波长λ和Λ有关。
η和Λ/λ对光栅等效薄膜折射率的影响如图 3所示。此时nh=3.48,nl=1。图 3a是当Λ/λ=0.5时,占空比对光栅等效折射率的影响。从图中可以看到,随着占空比的增加,TE偏振下的折射率要大于TM偏振的折射率,并且当η等于0或1时,两个偏振下的折射率相等。因为三角HCG可看为由几个不同占空比的矩形HCG构成,如图 2所示,所以三角HCG从上到下的折射率分布与图 3a中的曲线相似。图 3b是当η=0.5时,周期与入射波长比值Λ/λ对光栅等效折射率的影响。从图中可看到, 随着Λ/λ的增加,TM偏振下的光栅等效折射率变化要比TE偏振下的显著。
Figure 3. The influence of grating parameters and incident wavelength on the equivalent refractive index of grating
为了让非周期三角HCG实现透射光束偏转,需让非周期三角HCG相位呈线性分布[11],本文中设计的非周期三角HCG其整体相位分布是从左往右按照线性关系依次升高,具体见图 1。为设计出相位线性分布且具有较高透射率的三角HCG,首先需计算周期性三角HCG参量对其透射率和透射相位的影响,然后选择合适的光栅参量设计出透射光束可偏转非周期三角HCG。周期性三角HCG参量对其透射率和透射相位的影响可用严格耦合波法计算[7],为了保证计算结果的准确性,本文中将三角HCG细分为50个高度相同的小矩形来计算[17],具体如图 4所示。此时Λ=0.5μm, η=1。
-
Λ和η对周期性三角HCG相位和透射率的影响如图 5所示(入射光波长(TM)λ=1.55μm,光栅高度d=1μm)。之所以随着Λ和η变化,周期性三角HCG透射率会发生变化,是因为根据等效介质原理可知,Λ和η变化,则三角HCG的等效折射率会发生变化,等效折射率变化,则会影响光栅的透射率,因为薄膜透射和反射与介质薄膜的折射率有关[18]。而三角HCG相位随着Λ和η变化,也是因为Λ和η变化会影响三角HCG的等效薄膜折射率,根据(5)式可知,薄膜的相位变化是与入射波长、薄膜折射率、薄膜厚度有关的,所以三角HCG的周期和占空比发生变化,则相位也会发生变化。
$ \Delta \mathit{\Phi } = \frac{{2{\rm{ \mathsf{ π} }}}}{\lambda }nh $
(5) 式中,ΔΦ是薄膜的相位变化,h是薄膜厚度,n是薄膜折射率。
根据Λ和η对周期性三角HCG相位和透射率的影响,选择出一组透射率大于85%的光栅周期和占空比,根据(6)式对它们进行位置排列,相邻光栅条相位关系应与(7)式相符[9]:
$ {x_{n + 1}} = {x_n} + \left( {{\mathit{\Lambda }_n} + {\mathit{\Lambda }_{n + 1}}} \right)/2,(n = 0,1,2,3 \cdots ) $
(6) $ \begin{array}{l} \psi \left( {{x_{n + 1}}} \right) = \psi \left( {{x_n} + \left( {{\mathit{\Lambda }_n} + {\mathit{\Lambda }_{n + 1}}} \right)/2} \right) = \\ {k_0}{x_{n + 1}}\sin \theta \end{array} $
(7) 式中,k0是波数,λ=1.55μm(TM偏振),θ是光栅透射光束偏转角,其大小由光栅整体相位差Δψ, k0和w决定,如下式所示[8]:
$ \theta=\arcsin \left(\frac{\Delta \psi}{w k_{0}}\right) $
(8) 结构参量和排列顺序如图 6所示,具体参量数值可见表 1。从表 1可以发现,非周期三角HCG能实现大角度光束偏转是因为相位累积的结果[19],另外,本文中选择的三角HCG结构参量在电子束光刻的精度范围内。图 6a是不同位置选择的Λn和ηn(n=0,1,2…),图 6b是所设计的非周期三角HCG整体相位分布和理论相位分布。从图 6b可知, 本文中设计的非周期三角HCG总相位差Δψ=20.25rad,本文中w=10μm,根据(8)式计算可得本文中设计的光束偏转角为30°。
Table 1. Parameters and phases of different elements of non-periodic triangular HCG
period Λ/μm duty cycle η phase/rad 0.5 0.66 0.029 0.5 0.52 1.025 0.8 0.156 2.027 0.63 0.884 3.038 0.57 0.852 4.045 0.455 0.888 5.060 0.32 0.896 6.076 0.335 0.68 0.810 0.46 0.296 1.820 0.68 0.832 2.833 0.65 0.752 3.855 0.48 0.88 4.864 0.365 0.88 5.879 0.25 0.792 0.608 0.26 0.44 1.623 0.68 0.86 2.635 0.59 0.872 3.643 0.495 0.744 5.673 0.365 0.74 0.046 0.35 0.496 1.418 依据上述选择的光栅参量利用FDTD方法建立非周期三角HCG 2维模型并对其进行模拟计算,模型总共包含有20个不同周期和占空比的三角HCG,四周边界条件设置的是完美匹配层(perfectly matched layer, PML),其主要作用是解决在仿真区域边界上的反射问题,光源为高斯光束(TM偏振)。通过模拟计算可知, 非周期等腰三角HCG光束偏转角可达30.3°,如图 7所示。图 7a是由FDTD方法建立的非周期三角HCG模型, 图 7b是非周期等腰三角HCG整体的光强分布图,图 7c是在距离光栅透射面20μm和30μm处的光强分布。从图 7c可以看到, 随着离透射面的距离从20μm增加到30μm,光强的波峰也从x=-7.0086μm左移到x=-12.8664μm,据此得出透射光束偏转角为θ=arctan(5.8578/10)=30.3°,非常接近由(8)式计算得到的光束偏转角度。最终模拟得到的偏转角稍微偏差目标的30°,主要因为理想的相位分布是连续的,而实际设计的相位分布是离散的,因而导致实际设计的非周期三角HCG偏转角度与理论有略有偏差[9-10]。
Figure 7. a—HCG model of non-periodic triangle established by FDTD b—global light intensity distribution of non-periodic triangulated HCG c—light intensity at different height from transmission surface
另外根据计算发现,非周期三角HCG对1.55μm波长透射率为82.6%,略低于85%。这是因为在非周期性的三角形HCG中,是由多个不同周期和占空比的三角HCG组合,不是无限大,而周期三角HCG在模拟计算时是无限大(即无限个周期),所以它们总的透射率会略低于周期性HCG[9],具体如图 8所示。
-
图 9是不同波长入射光照射到第2.1节中设计的非周期三角HCG后的光束偏转角。从图中可以看出, 当入射光波长在1.5μm~1.6μm之间变化时,非周期等腰三角HCG依然能够实现大角度偏转。当入射光波长从1.5μm增加到1.6μm,透射光束偏转角从28.323°增加到31.85°,可实现透射光束3.527°的调谐。这种能够实现透射光束偏转角可调谐特性对在雷达系统中的应用具有重要意义[20-22]。
入射光波长在λ=1.55μm附近变化时, 非周期等腰三角HCG能够保持大角度透射光束偏转,主要是因为非周期等腰三角HCG每一个光栅条相位变化较小,而且整体的光栅相位依然能够呈大致的线性关系,如图 10所示。从图 10可知,虽然波长变化使光栅整体相位偏离了原先的线性关系,因为根据(5)式可知波长会影响光栅等效薄膜相位,但随着位置的增加,相位依然随之增加,保持着准线性分布。另外随着波长的增加光栅透射光束偏转角增加可由(8)式解释,因为总相位差Δψ随着波长增加变化不大,w一直保持10μm,而k0却随着波长增大而减小,这就导致透射光束偏转角随着波长增加而增加,实现了透射光束偏转角的调谐。
-
不仅通过改变入射波长可以实现光束偏转角的调谐,还可以通过改变光栅低折射率材料nl实现光束偏转角的调谐,如图 11和图 12所示。图 11是不同nl下的电场分布,图 12是不同nl下的非周期三角HCG光束偏转角。从图 11、图 12可以发现, 随着nl的增加,在入射波长为1.55μm的条件下,非周期等腰三角HCG的光束偏转角逐渐从30.3°下降到19.3°,实现了11°的光束偏转角调谐,这为将来利用液晶作为低折射率材料实现电控光束偏转提供了理论基础。
Figure 12. The influence of low refractive index materials on the deflection angle of transmitted light
折射率nl对光束偏转角的影响可以解释如下。由(2)式可知,随着折射率nl的增大,三角HCG的等效折射率增大,这会引起非周期三角HCG的相位差减小,从而导致偏角减小,具体如图 13所示。图 13为不同nl条件下用严格耦合波法计算得到的非周期三角形HCG的相位分布。可以发现, 随着nl的增大,非周期三角HCG相位差Δψ在下降。nl=1时,Δψ=20.25rad; nl=1.4时,Δψ降到16.87rad。根据(8)式可知,相位差的减小会使偏转角度减小。此外,从图 13可以发现,nl分别是1.1和1.2时,非周期三角HCG相位在x= 0.5μm和x=10μm之间线性分布; n1是1.3和1.4时,非周期三角HCG相位在x=1μm和x=10μm之间线性分布; nl=1.4时,非周期三角HCG有效宽度是9μm,代入(8)式计算可得光束偏转角为19.27°,这一结果与图 12中的模拟值高度一致,可以表明该模型是准确的。
利用高对比度光栅实现光束大角度偏转
Large angle deflection of beam using high contrast grating
-
摘要: 为了分析高折射率对比度光栅(HCG)参量和入射波长对光束偏转角的影响,采用严格耦合波法设计了透射光束可偏转非周期三角HCG,并通过时域有限差分法证明了所设计的非周期三角HCG可实现30.3°的透射光束偏转。结果表明,当低折射率介质材料折射率从1增加到1.4时,透射光束偏转角可实现11°的调谐;当入射波长从波长1.5μm增加到1.6μm时,该非周期三角HCG可实现3.527°的透射光束偏转角调谐。这一研究结果可对将来制备高性能光束偏转光栅提供理论指导。Abstract: In order to analyze the influence of high contrast grating(HCG) parameters and incident wavelength on beam deflection angle, the rigorous coupled wave method was used to design a deflectable non-periodic triangular HCG, and it is proved that 30.3° beam deflection can be achieved in the designed non-periodic triangular HCG by the finite-difference time-domain (FDTD). At the same time, it is found that when the refractive index of the low refractive index material increases from 1 to 1.4, the deflection angle of the transmission beam can achieve 11° tuning, and when the incident wavelength increases from 1.5μm to 1.6μm, the aperiodic triangular HCG can achieve 3.527° transmission beam deflection angle tuning. The results can provide theoretical guidance for the fabrication of high-performance beam deflection gratings in the future.
-
Key words:
- gratings /
- beam deflection /
- transmission beam /
- non-periodic
-
Table 1. Parameters and phases of different elements of non-periodic triangular HCG
period Λ/μm duty cycle η phase/rad 0.5 0.66 0.029 0.5 0.52 1.025 0.8 0.156 2.027 0.63 0.884 3.038 0.57 0.852 4.045 0.455 0.888 5.060 0.32 0.896 6.076 0.335 0.68 0.810 0.46 0.296 1.820 0.68 0.832 2.833 0.65 0.752 3.855 0.48 0.88 4.864 0.365 0.88 5.879 0.25 0.792 0.608 0.26 0.44 1.623 0.68 0.86 2.635 0.59 0.872 3.643 0.495 0.744 5.673 0.365 0.74 0.046 0.35 0.496 1.418 -
[1] MATEUS C F R, HUANG M C Y, DENG Y, et al. Ultrabroadband mirror using low-index cladded subwavelength grating[J]. IEEE Photonics Technology Letters, 2004, 16(2): 518-520. doi: 10.1109/LPT.2003.821258 [2] TSUNEMI Y, YOKOTA N, MAJIMA S, et al. 1.55μm VCSEL with polarization-independent HCG mirror on SOI[J]. Optics Express, 2013, 21(23): 28685-28692. doi: 10.1364/OE.21.028685 [3] HUANG M C Y, ZHOU Y, CHANG C J. Single mode high-contrast subwavelength grating vertical cavity surface emitting lasers[J]. A-pplied Physics Letters, 2008, 92(17): 171108. doi: 10.1063/1.2917447 [4] ZHOU Y, HUANG M C, CHANG H C J. Tunable VCSEL with ultra-thin high contrast grating for high-speed tuning[J]. Optics Express, 2008, 16(18): 14221-14226. doi: 10.1364/OE.16.014221 [5] LI K, RAO Y, CHASE C, et al. Monolithic high-contrast metastructure for beam-shaping VCSELs[J]. Optica, 2018, 5(1): 10-13. doi: 10.1364/OPTICA.5.000010 [6] GAO S, PARK C S, LEE S S, et al. A highly efficient bifunctional dielectric metasurface enabling polarization-tuned focusing and deflection for visible light[J]. Advanced Optical Materials, 2019, 7(9): 1801337. doi: 10.1002/adom.201801337 [7] FANG W J, HUANG Y Q, DUAN X F et al. Non-periodic high-index contrast gratings reflector with large-angle beam forming ability[J]. Optics Communications, 2016, 367: 6-11. doi: 10.1016/j.optcom.2016.01.025 [8] CARLETTI L, MALUREANU R, JESPER M, et al. High-index-contrast grating reflector with beam steering ability for the transmitted beam[J]. Optics Express, 2011, 19(23): 23567. doi: 10.1364/OE.19.023567 [9] MA Ch L, HUANG Y Q, DUAN X F, et al. High-transmissivity non-periodic sub-wavelength high-contrast grating with large-angle beam-steering ability[J]. Chinese Optics Letters, 2014, 12(12): 120501. doi: 10.3788/COL201412.120501 [10] FANG W J, HUANG Y Q, FEI J, et al. Concentric circular focusing reflector realized using high index contrast gratings[J]. Optics Communications, 2017, 402: 572-576. doi: 10.1016/j.optcom.2017.06.096 [11] GORO T, YOSUKE T, MOE T, et al. Thermally controlled Si photonic crystal slow light waveguide beam steering device[J]. Optics Express, 2018, 26(9): 11529. doi: 10.1364/OE.26.011529 [12] HUANG Y W, ZHANG X, ZANG J W, et al. Non-periodic wide-angle beam steering HCG array for application in VCSEL[J]. Journal of Infrared and Millimeter Waves, 2018, 37(1): 20-29(in Ch-inese). [13] CHUNG I S, IAKOVLEV V, SIRBU A, et al. Broadband MEMS-tunable high-index-contrast subwavelength grating long-wavelength VCSEL[J]. IEEE Journal of Quantum Electronics, 2010, 46(9): 1245-1253(in Chinese). doi: 10.1109/JQE.2010.2047494 [14] XU P, TANG Sh T, YUAN X, et al. Design of an embedded trico-lor-shifting device[J]. Acta Physica Sinica, 2018, 67(2): 60-71(in Chinese). [15] JOSEPH N M, DENNIS W P, MARK S M. Design of binary subwavelength diffractive lenses by use of zeroth-order effective-medium theory[J]. Journal of the Optical Society of America, 1999, A16(5): 1157-1167. [16] JIANG X W, WU H. Research of 1-D sub-wavelength grating anti-reflection[J]. Chinese Journal of Luminescence, 2017, 38(2): 177-181(in Chinese). doi: 10.3788/fgxb20173802.0177 [17] JING X F, ZHANG J Ch, JIN Sh Zh, et al. Design of highly efficient transmission gratings with deep etched triangular grooves[J]. Applied Optics, 2012, 51(33): 7920-7933. doi: 10.1364/AO.51.007920 [18] GRUNHUT V, DEPING R A. Influence of the sign of the refractive index in the reflectivity of a metamaterial surface with localized roughness[J]. European Physical Journal, 2011, D62(2): 227-236. doi: 10.1140/epjd/e2011-10438-4 [19] SHI W Y, DENG W M, LIU W N, et al. Rectangular dielectric metagrating for high-efficiency diffraction with large-angle deflection[J]. Chinese Optics Letters, 2020, 18(7): 073601. doi: 10.3788/COL202018.073601 [20] WANG Ch M, PENG Z H, LIU Y G, et al. Radial sub-aperture coherence method to achieve beam steering with high precision and stability[J]. Optics Express, 2019, 27(5): 6331-6347. doi: 10.1364/OE.27.006331 [21] LI Zh S. Design of liquid crystal polarization grating and study of its beam deflection characteristics[D]. Changchun: Changchun Institute of Optics, Fine Mechanics and Physics, Chinese Academy of Sciences, 2019: 89-90(in Chinese). [22] PAN W J, WU Zh Y, ZHANG X L. Identification of aircraft wake vortex based on k-nearest neighbor[J]. Laser Technology, 2020, 44(4): 471-477(in Chinese).