CN114762191A - 基于高阶旋转对称单元的非线性几何相位超表面 - Google Patents
基于高阶旋转对称单元的非线性几何相位超表面 Download PDFInfo
- Publication number
- CN114762191A CN114762191A CN202280000514.XA CN202280000514A CN114762191A CN 114762191 A CN114762191 A CN 114762191A CN 202280000514 A CN202280000514 A CN 202280000514A CN 114762191 A CN114762191 A CN 114762191A
- Authority
- CN
- China
- Prior art keywords
- nonlinear
- order
- rotational symmetry
- super
- geometric phase
- 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
Links
- 230000010287 polarization Effects 0.000 claims abstract description 46
- 239000000463 material Substances 0.000 claims abstract description 40
- 239000000758 substrate Substances 0.000 claims abstract description 23
- 230000009021 linear effect Effects 0.000 claims abstract description 9
- 239000013078 crystal Substances 0.000 claims abstract description 4
- 238000005086 pumping Methods 0.000 claims abstract description 4
- 229910052751 metal Inorganic materials 0.000 claims description 36
- 239000002184 metal Substances 0.000 claims description 36
- 230000000737 periodic effect Effects 0.000 claims description 7
- GQYHUHYESMUTHG-UHFFFAOYSA-N lithium niobate Chemical compound [Li+].[O-][Nb](=O)=O GQYHUHYESMUTHG-UHFFFAOYSA-N 0.000 claims description 5
- 229910052782 aluminium Inorganic materials 0.000 claims description 4
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 claims description 4
- 239000011521 glass Substances 0.000 claims description 4
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 claims description 4
- 229910052737 gold Inorganic materials 0.000 claims description 4
- 239000010931 gold Substances 0.000 claims description 4
- AMGQUBHHOARCQH-UHFFFAOYSA-N indium;oxotin Chemical compound [In].[Sn]=O AMGQUBHHOARCQH-UHFFFAOYSA-N 0.000 claims description 4
- 239000010980 sapphire Substances 0.000 claims description 4
- 229910052594 sapphire Inorganic materials 0.000 claims description 4
- 230000001105 regulatory effect Effects 0.000 abstract description 8
- 230000033228 biological regulation Effects 0.000 description 9
- 230000005684 electric field Effects 0.000 description 8
- 230000003287 optical effect Effects 0.000 description 8
- 230000004044 response Effects 0.000 description 8
- 238000000034 method Methods 0.000 description 7
- 238000004088 simulation Methods 0.000 description 7
- RKTYLMNFRDHKIL-UHFFFAOYSA-N copper;5,10,15,20-tetraphenylporphyrin-22,24-diide Chemical compound [Cu+2].C1=CC(C(=C2C=CC([N-]2)=C(C=2C=CC=CC=2)C=2C=CC(N=2)=C(C=2C=CC=CC=2)C2=CC=C3[N-]2)C=2C=CC=CC=2)=NC1=C3C1=CC=CC=C1 RKTYLMNFRDHKIL-UHFFFAOYSA-N 0.000 description 6
- 230000001808 coupling effect Effects 0.000 description 6
- 230000009022 nonlinear effect Effects 0.000 description 6
- 238000009826 distribution Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000005502 phase rule Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000002708 enhancing effect Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000001737 promoting effect Effects 0.000 description 2
- 239000005361 soda-lime glass Substances 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000019771 cognition Effects 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01Q—ANTENNAS, i.e. RADIO AERIALS
- H01Q15/00—Devices for reflection, refraction, diffraction or polarisation of waves radiated from an antenna, e.g. quasi-optical devices
- H01Q15/0006—Devices acting selectively as reflecting surface, as diffracting or as refracting device, e.g. frequency filtering or angular spatial filtering devices
- H01Q15/0086—Devices acting selectively as reflecting surface, as diffracting or as refracting device, e.g. frequency filtering or angular spatial filtering devices said selective devices having materials with a synthesized negative refractive index, e.g. metamaterials or left-handed materials
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B5/00—Optical elements other than lenses
- G02B5/30—Polarising elements
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Optics & Photonics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Optical Modulation, Optical Deflection, Nonlinear Optics, Optical Demodulation, Optical Logic Elements (AREA)
Abstract
本公开提供一种基于高阶旋转对称单元的非线性几何相位超表面，自下而上依次包括：衬底(1)；非线性材料层(2)；超表面结构层(3)，包括多个具有局域高阶旋转对称性的超构单元(4)，超构单元(4)呈周期性排列形成具有全局高阶旋转对称性的晶格结构；非线性几何相位超表面在圆偏振态基波的泵浦下，能够同时产生与基波相同和相反圆偏振态的非线性谐波，且所产生非线性谐波的非线性几何相位受超构单元局域旋转对称阶数和晶格结构的调控。
Description
技术领域
本公开涉及电磁波调制技术领域，具体涉及一种基于高阶旋转对称单元的非线性几何相位超表面。
背景技术
在线性光学中，圆偏振光与各向异性超构单元作用可以产生具有相反手性的圆偏光，且转换后的光携带几何相位Φ＝2σθ，称为Pancharatnam-Berry几何相位。其中θ为超构单元的旋转角，σ＝±1代表入射光处于左旋或右旋圆偏振态(LCP或RCP)。
类似的概念也被引入到非线性光学中，当偏振态为σ的基波光入射到非线性超表面时，产生的+σ或-σ偏振态n次谐波具有几何相位(n-1)σθ或(n+1)σθ，称为非线性几何相位。非线性几何相位不依赖于结构尺寸和材料色散，且可以通过旋转超构单元实现0～2π的连续调控，因此被广泛应用于光的非线性调控。
在非线性超表面中，圆偏光n次谐波的产生需要满足选择定则：n＝Nm±1，其中N为整数，m为超构单元旋转对称阶数，符号+和-代表谐波与基波具有相同或相反的圆偏振态。以倍频为例，具有1重旋转对称(C1)结构中可以产生与基波相同和相反圆偏振的倍频光，而三重旋转对称(C3)结构中只能产生相反圆偏振的倍频光，其余对称结构则不能发生圆偏光倍频。
因此，对于具体某一非线性过程，非线性几何相位只能在少数特定的旋转对称结构中产生。尽管同一过程可能在不同的旋转对称结构中产生，但非线性几何相位规律并不因此改变，比如在C1和C3结构中，偏振态与基波相反的圆偏倍频光几何相位都为3σθ，这些限制关系制约了非线性几何相位的产生和调控灵活度。
发明内容
(一)要解决的技术问题
针对上述问题，本公开提供了一种基于高阶旋转对称单元的非线性几何相位超表面，用于解决传统超表面对非线性几何相位的调控受限制较大等技术问题。
(二)技术方案
本公开提供了一种基于高阶旋转对称单元的非线性几何相位超表面，自下而上依次包括：衬底；非线性材料层；超表面结构层，包括多个具有局域高阶旋转对称性的超构单元，超构单元呈周期性排列形成具有全局高阶旋转对称性的晶格结构；非线性几何相位超表面在圆偏振态基波的泵浦下，能够同时产生与基波相同和相反圆偏振态的非线性谐波。
进一步地，同时产生与基波相同和相反圆偏振态的非线性谐波的谐波阶数n为2或3。
进一步地，超构单元的局域旋转对称阶数m、晶格结构的全局旋转对称阶数q、谐波阶数n满足以下条件：n+1和n-1均能够被局域旋转对称阶数m与全局旋转对称阶数q的最大公因数整除。
进一步地，非线性谐波的谐波阶数n为2时，非线性谐波的几何相位Φ＝±mθ或Φ＝±3mθ，其中，θ为超构单元的旋转角；非线性谐波的谐波阶数n为3时，非线性谐波的几何相位Φ＝±2mθ、±mθ或0。
进一步地，具有局域高阶旋转对称性的超构单元的局域旋转对称阶数m大于4。
进一步地，超构单元呈周期性排列包括采用四方点阵或六角点阵的排布形成全局旋转对称阶数q为4或6的晶格结构。
进一步地，超构单元为亚波长金属结构，亚波长金属结构的材料包括但不限于金、铝中的一种。
进一步地，亚波长金属结构具有m阶旋转对称性即可，包括但不限于正m边形、正m星形，其中，m为超构单元的局域旋转对称阶数。
进一步地，非线性材料层包括但不限于氧化铟锡、铌酸锂中的一种；衬底包括但不限于玻璃、蓝宝石中的一种。
进一步地，超构单元的厚度h＜λ0/2、周期p＜λo；非线性材料层的厚度t＜λ0/2；衬底的厚度为d＞＞λ0，其中，λ0为谐波波长。
(三)有益效果
本公开提供的基于高阶旋转对称单元的非线性几何相位超表面，利用超构单元和晶格结构的耦合效应，能在具有局域高阶旋转对称性的超构单元中同时产生与基波相同和相反圆偏振态的非线性谐波，且所产生非线性谐波的非线性几何相位受超构单元局域旋转对称阶数和晶格结构的调控；这打破了传统非线性超表面中的选择定则与几何相位规律，同时也丰富了非线性光场调控的方式。
附图说明
图1示意性示出了根据本公开实施例中基于高阶旋转对称单元的非线性几何相位超表面的结构示意图；
图2示意性示出了根据本公开实施例中三种不同亚波长金属结构的非线性几何相位超表面的三维视图和俯视图；
图3示意性示出了根据本公开实施例1中设计的非线性几何相位超表面的倍频仿真结果；
图4示意性示出了根据本公开实施例1中利用超构单元结构设计的偏折器件仿真结果；
图5示意性示出了根据本公开实施例2中非线性几何相位超表面的结构示意图以及三倍频仿真结果；
附图标记说明：
1.衬底，2.非线性材料层，3.超表面结构层，4.超构单元。
具体实施方式
为使本公开的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本公开进一步详细说明。
在此使用的术语仅仅是为了描述具体实施例，而并非意在限制本公开。在此使用的术语“包括”、“包含”等表明了所述特征、步骤、操作和/或部件的存在，但是并不排除存在或添加一个或多个其他特征、步骤、操作或部件。
在此使用的所有术语(包括技术和科学术语)具有本领域技术人员通常所理解的含义，除非另外定义。应注意，这里使用的术语应解释为具有与本说明书的上下文相一致的含义，而不应以理想化或过于刻板的方式来解释。
本公开的实施例提供了一种基于高阶旋转对称单元的非线性几何相位超表面，请参见图1，自下而上依次包括：衬底1；非线性材料层2；超表面结构层3，包括多个具有局域高阶旋转对称性的超构单元4，超构单元4呈周期性排列形成具有全局高阶旋转对称性的晶格结构；该非线性几何相位超表面在圆偏振态基波的泵浦下，能够同时产生与基波相同和相反圆偏振态的非线性谐波，且所产生非线性谐波的非线性几何相位受超构单元局域旋转对称阶数和晶格结构的调控。
超表面结构层3的对称性由单个超构单元4的局域对称性和多个超构单元4的全局对称性共同决定，因此，本公开通过设计超构单元4的结构、尺寸以及其晶格结构的全局对称性来调节其非线性光学响应。具体地，本公开包括超表面结构层3、具有垂直方向高非线性响应的非线性材料层2以及非线性响应可忽略的衬底1，非线性响应主要在非线性材料层2中发生，超构单元4呈周期性排列形成的阵列主要用于调控电场对称性和增强非线性材料层2中z方向电场，并促进非线性谐波耦合到自由空间；衬底1则主要用于支撑非线性材料层2和超表面结构层3。
其中，超表面结构层3采用具有局域高阶旋转对称性的超构单元4，且该超构单元4排列形成具有全局高阶旋转对称性的晶格结构，利用超构单元4与晶格结构的耦合效应，在非线性超表面中可同时产生与基波相同和相反两种圆偏振态的非线性谐波，且所产生非线性谐波的非线性几何相位受超构单元4的旋转对称阶数和晶格结构调控。这打破了传统非线性超表面中的选择定则与几何相位规律，也丰富了非线性光场调控的手段。
在上述实施例的基础上，同时产生与基波相同和相反圆偏振态的非线性谐波的谐波阶数n为2或3。
本公开的非线性几何相位超表面既适用于二次谐波的调控也适用于三次谐波的调控，对应的非线性材料层2相应地选择二阶或三阶非线性材料，同时，利用超构单元4调控电场对称性，即可得到谐波阶数n为2或3的非线性谐波。
在上述实施例的基础上，为了实现上述技术效果，本公开中超构单元的局域旋转对称阶数m、晶格结构的全局旋转对称阶数q、非线性谐波的谐波阶数n需满足以下条件：n+1和n-1均能够被局域旋转对称阶数m与全局旋转对称阶数q的最大公因数整除。
当超构单元4的局域高阶旋转对称阶数m与晶格结构的全局高阶旋转对称阶数q不匹配时，即q不能被m整除时，传统选择定则可以被打破，这时选择定则应考虑整体结构的旋转对称阶数，其可以表示为超构单元4的局域旋转对称阶数m与晶格结构的全局旋转对称阶数q的最大公因数。在设计超构单元4时，按照以下规则进行超构单元4的旋转对称阶数m的具体设计，即m需满足：(n+1)和(n-1)均能够被局域旋转对称阶数m与全局旋转对称阶数q的最大公因数整除，且m为大于4的正整数。其中，n为非线性谐波的谐波阶数。
本公开设计得到的非线性几何相位超表面可以同时产生与基波相同和相反两种圆偏振的非线性谐波，且非线性几何相位由光场的偏振态和超构单元的旋转对称阶数，以及晶格结构共同决定。其中，所得二次谐波(倍频光)的几何相位为Φ＝±mθ或Φ＝±3mθ，其中，θ为超构单元4的旋转角，所得三次谐波(三倍频光)的几何相位为Φ＝±2mθ、±mθ或0。
下面说明本公开中谐波产生和几何相位调控的原理，如下：
以二阶非线性效应为例，由于选择定则，人们认为非线性超表面中圆偏光的倍频只能发生在C1和C3旋转对称结构中，其中C1结构中可以产生与基波相同和相反两种圆偏振的倍频光，而C3结构中只能产生与基波相反圆偏振的倍频光。与基波相同和相反圆偏振倍频光的非线性几何相位分别为Φ＝±θ和Φ＝±3θ。
因此，在四方晶格结构中，对于二阶非线性效应，m为4以上的奇数阶旋转对称超构单元都能同时产生两种圆偏振的倍频光。且由于超构单元与晶格结构的耦合效应，两种圆偏振倍频光都可以携带Φ＝±mθ的几何相位；其中，θ为超构单元4的旋转角。比如对于C5结构，倍频光的几何相位可以为Φ＝±5θ，而C7结构中倍频光几何相位可以达到Φ＝±7θ。这表明非线性几何相位受超构单元旋转对称阶数调控。
对于三阶非线性效应，旋转对称阶数m为4以上且不为4的整数倍的超构单元也能同时产生两种圆偏振的三倍频光，且偏振态与基波相同的三倍频光几何相位为Φ＝±2mθ(m为奇数)或Φ＝±mθ(m为偶数)，而偏振态与基波相反的三倍频光几何相位为0。
在六角晶格结构中，对于二阶非线性效应，旋转对称阶数m为4以上的奇数且不为3的整数倍的超构单元中都能产生两种圆偏振倍频光，其中与基波具有相同圆偏振态的倍频光可以携带Φ＝±mθ的几何相位，与基波具有相反圆偏振态的倍频光可以携带Φ＝±3mθ的几何相位，与四方晶格的情况稍有不同。这表明非线性几何相位也受晶格结构调控。
对于三阶非线性效应，旋转对称阶数m为4以上且不为3的整数倍的超构单元中都能产生两种圆偏振三倍频光，且两种圆偏振三倍频光都可以携带Φ＝±2mθ(m为奇数)或Φ＝±mθ(m为偶数)的几何相位。
以上总结如下表1：
表1
| 谐波阶数 | 四方晶格排列(q＝4) | 六角晶格排列(q＝6) |
| n＝2 | m为4以上的奇数 | m为4以上的奇数且不为3的整数倍 |
| n＝3 | m为4以上且不为4的整数倍 | m为4以上且不为3的整数倍 |
在上述实施例的基础上，本公开具有局域高阶旋转对称性的超构单元4的局域旋转对称阶数m＞4。
其中，超构单元4的局域旋转对称阶数为3或4时虽然也能在特定晶格中同时产生与基波相同和相反圆偏振态的非线性谐波，但此时超构单元4与晶格结构的耦合效应还不足够明显；当局域旋转对称阶数为3或4时，相反圆偏振态的非线性谐波的强度远大于相同圆偏振态的非线性谐波的强度，所以本公开选择超构单元4的局域旋转对称阶数m＞4。
在上述实施例的基础上，超构单元4呈周期性排列包括采用四方点阵或六角点阵的排布形成全局旋转对称阶数q为4或6的晶格结构。
超构单元4采用四方点阵或六角点阵的周期性排列方式形成具有特定排列方式的晶格结构，即该晶格结构的晶格旋转对称阶数q为4或6，晶格结构的高阶旋转对称性能与超构单元4的耦合效果更好。
在上述实施例的基础上，超构单元4为亚波长金属结构，本公开亚波长金属结构的材料可以包括但不限于金、铝中的一种。
亚波长金属结构以周期性排列方式形成晶格结构，主要用于调控电场对称性和增强非线性材料层2中z方向电场，并促进非线性谐波耦合到自由空间；亚波长金属结构以自身的中心旋转一定角度θ可以带来相应的几何相位变化。
在上述实施例的基础上，亚波长金属结构具有m阶旋转对称性即可，包括但不限于正m边形、正m星形，其中，m为超构单元4的局域旋转对称阶数。
亚波长金属结构例如具有五重(C5)、七重(C7)局域高阶旋转对称性，其具体结构如图2所示，图2(a)为正五星形，图2(b)为正五边形，图2(c)为正七边形，当然，这里并不限定于以上三种亚波长金属结构，其它具有局域高阶旋转对称性的结构也适用于本公开的技术方案。需要说明的是，当亚波长金属结构为具有五重(C5)、七重(C7)等5以上的奇数阶旋转对称性的结构时，在n＝2，q＝4时效果最优。
在上述实施例的基础上，非线性材料层2包括但不限于氧化铟锡、铌酸锂中的一种；衬底1包括但不限于玻璃、蓝宝石中的一种。
非线性材料层2的材料可以为氧化铟锡(ITO)、铌酸锂等，其中铌酸锂具有高二阶非线性响应，ITO同时具有较强的二阶和三阶非线性响应；衬底1的材料可以为玻璃、蓝宝石等，衬底1的材料只需要满足非线性响应低、能起支撑作用即可。
在上述实施例的基础上，超构单元4的厚度h＜λ0/2、周期p＜λ0；非线性材料层2的厚度t＜λ0/2；衬底1的厚度为d＞＞λ0，其中，λ0为谐波波长。
如图2所示，图2上半部分为非线性几何相位超表面的三维视图，其中，超构单元4的厚度和周期分别为h和p，非线性材料层2和衬底1的厚度分别为t和d；超构单元4的周期p需要小于谐波波长λ0，是为了避免产生高阶衍射；超构单元4的厚度h经验取值小于谐波半波长λ0/2，有利于效率优化。其中，二次谐波的波长是基波波长的一半，三次谐波的波长是基波波长的三分之一。
图2下半部分为非线性几何相位超表面的俯视图，θ为亚波长金属结构的旋转角，将亚波长金属结构中垂直于x轴方向的初始对称轴作为对称轴并将垂直于对称轴的x’轴(初始位置时与x轴重合)作为参考轴，此时θ＝0；当旋转亚波长金属结构时，对称轴随之旋转且参考轴x’也随之旋转，则参考轴x’与x轴之间的夹角即为旋转角θ，可以通过设计调整每个亚波长金属结构的旋转角θ并得到具有不同旋转角θ的亚波长金属结构单元，继而可以通过周期性排列得到亚波长金属结构阵列。
本公开通过设计具有五重(C5)、七重(C7)等高阶旋转对称性的超构单元、具有垂直方向高非线性响应的非线性材料层，并利用超构单元与晶格结构的耦合效应，可以同时产生与基波相同和相反两种圆偏振的非线性谐波，且非线性几何相位由光场的偏振态和超构单元的旋转对称阶数，以及晶格结构共同决定；颠覆了人们对非线性超表面中谐波产生和几何相位的认知，对非线性光场调控的原理和应用具有重大意义。
下面通过具体实施方式对本公开作进一步说明。在以下实施例中对上述基于高阶旋转对称单元的非线性几何相位超表面进行具体说明。但是，下述实施例仅用于对本公开进行例示，本公开的范围不限于此。
实施例1：
本实施例选择具有C5旋转对称性的非线性超表面为例进行说明。
该非线性几何相位超表面由具有周期性阵列排布的超构单元4、非线性材料层2和衬底1构成，其中，超构单元4为亚波长金属结构，该亚波长金属结构单元为正五边形，按四方点阵周期性排列形成晶格结构；亚波长金属结构单元的材料为金，其介电常数可以从Palik光学手册中获得。
非线性材料层2的材料为ITO，ITO是一种介电常数实部具有零点的材料，在零介电常数波长附近非线性相互作用可以得到大大增强。ITO的相对介电常数可以用德鲁德模型描述为ε＝ε∞-ωp 2/(ω2+iΓω)，其中ω为光的角频率。公式里其它参数取值可以通过实验制备条件进行控制，这里采用商用ITO薄膜的实验测试结果，即ε∞＝3.0404，ωp＝2.5695×1013Hz，Γ＝1.4832×106Hz，对应介电常数实部零点位于1300nm波长附近。ITO的最大二阶非线性系数矩阵元为d33＝0.18pm/V。
衬底1的材料为钠钙玻璃，其介电常数可以从Palik光学手册中获得。采用COMSOLMultiphysics软件模拟该结构中的倍频，基波波段设置为1250nm～1350nm，优化的亚波长金属结构单元的周期和厚度分别为p＝600nm，h＝60nm，正五边形的边长为200nm。非线性材料层2的厚度t为15nm，衬底1的厚度d为1mm。
图3为设计的非线性几何相位超表面的倍频仿真结果，其中图3(a)、图3(b)分别为左旋和右旋圆偏振基波入射下的倍频光(即二次谐波)不同圆偏分量强度，图3(c)、图3(d)分别为左旋圆偏基波光入射下产生的左旋和右旋圆偏振倍频光(即二次谐波)的几何相位。
倍频光强随波长变化如图3(a)和图3(b)所示，在左旋和右旋基波光入射下，该结构中都可以产生两种圆偏振的倍频光，且在ITO介电常数零点附近，倍频光强度明显得到了增强。这意味着非线性超表面中的常规选择定则被打破。
图3(c)和图3(d)展示了左旋基波光入射下，左旋和右旋圆偏振倍频光的几何相位随超构单元旋转角的变化。可以看出，亚波长金属结构单元每旋转9°，左旋和右旋圆偏振倍频光的相位分别变化+45°和-45°，对应几何相位分别满足Φ＝5θ和Φ＝-5θ。
接着，利用旋转亚波长金属结构单元设计了一个非线性光束偏折器件来验证这种相位梯度。选取24个旋转角依次增加3°的金属结构单元沿x方向排列，组成超构单元；再将超构单元周期性排布。当左旋圆偏振的平面基波光从下方入射时，产生的左旋和右旋倍频光将分别向+x和-x方向偏转，偏转角分别为ψ＝±arcsin(λ/P)。这里λ为倍频光波长，P＝24p为超级单元长度，p为600nm。
图4为利用上述超构单元结构设计的偏折器件仿真结果，其中图4(a)为截面倍频光不同圆偏振分量的归一化电场分布，图4(b)为远场倍频光强分布。
具体地，图4(a)展示了1300nm左旋基波光入射下，倍频光两种圆偏振分量的归一化电场分布仿真结果。可以看出，不同偏振的倍频光向不同的方向偏转。远场分析表明，左旋和右旋圆偏振倍频光的衍射角分别为+5.2°和-5.2°，与理论计算结果一致，如图4(b)所示。
实施例2：
本实施例选择具有C7旋转对称性的非线性超表面为例进行说明。
该非线性几何相位超表面由具有周期性阵列排布的超构单元4、非线性材料层2和衬底1构成，如图5(a)所示。其中，超构单元4为亚波长金属结构，该亚波长金属结构单元为正七星形，按六角点阵周期性排列形成晶格结构；亚波长金属结构单元的材料为铝，其介电常数可以从Palik光学手册中获得。
非线性材料层2的材料仍然为ITO，与实施例1相同，不过这里用到的是ITO的三阶非线性效应。ITO的相对介电常数在不同制备条件下具有不同的取值，这里用德鲁德模型描述为ε＝ε∞-ωp 2/(ω2+iΓω)，并取ε∞＝1.0878，ωp＝1.0419×1013Hz，Γ＝4.0950×106Hz，对应介电常数实部零点位于2000nm波长附近。ITO的最大三阶非线性系数矩阵元为χ1111＝3.5×10-18m2/V2。
衬底1的材料为钠钙玻璃，其介电常数可以从Palik光学手册中获得。采用COMSOLMultiphysics软件模拟该结构中的三倍频，基波波段设置为1900nm～2100nm，优化的亚波长金属结构单元的周期和厚度分别为p＝600nm，h＝100nm，正七星形的棱长和宽度分别为l＝120nm和w＝130nm。非线性材料层2的厚度t为15nm，衬底1的厚度d为1mm。
图5(b)～(d)为设计的非线性几何相位超表面的三倍频仿真结果，其中图5(b)为左旋和右旋圆偏振基波入射下的三倍频光(即三次谐波)不同圆偏分量强度，图5(c)、图5(d)分别为左旋圆偏基波光入射下产生的左旋和右旋圆偏振三倍频光(即三次谐波)的几何相位。
三倍频光强随波长变化如图5(b)所示，在左旋和右旋基波光入射下，该结构中都可以产生两种圆偏振的三倍频光。由于非线性过程同时受电场强度和交叠积分影响，三倍频最强位置并不完全位于ITO的介电常数零点波长。这意味着非线性超表面中的常规选择定则被打破。
图5(c)和图5(d)展示了左旋基波光入射下，左旋和右旋圆偏振三倍频光的几何相位随超构单元旋转角的变化。可以看出，亚波长金属结构单元每旋转4.3°，左旋和右旋圆偏振三倍频光的相位分别变化+60°和-60°，对应几何相位分别满足Φ＝14θ和Φ＝-14θ。
以上所述的具体实施例，对本公开的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本公开的具体实施例而已，并不用于限制本公开，凡在本公开的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本公开的保护范围之内。
Claims (10)
1.一种基于高阶旋转对称单元的非线性几何相位超表面，其特征在于，自下而上依次包括：
衬底(1)；
非线性材料层(2)；
超表面结构层(3)，包括多个具有局域高阶旋转对称性的超构单元(4)，所述超构单元(4)呈周期性排列形成具有全局高阶旋转对称性的晶格结构；
所述非线性几何相位超表面在圆偏振态基波的泵浦下，能够同时产生与所述基波相同和相反圆偏振态的非线性谐波。
2.根据权利要求1所述的基于高阶旋转对称单元的非线性几何相位超表面，其特征在于，所述同时产生与所述基波相同和相反圆偏振态的非线性谐波的谐波阶数n为2或3。
3.根据权利要求2所述的基于高阶旋转对称单元的非线性几何相位超表面，其特征在于，所述超构单元(4)的局域旋转对称阶数m、所述晶格结构的全局旋转对称阶数q、所述谐波阶数n满足以下条件：
n+1和n-1均能够被所述局域旋转对称阶数m与所述全局旋转对称阶数q的最大公因数整除。
4.根据权利要求3所述的基于高阶旋转对称单元的非线性几何相位超表面，其特征在于，所述非线性谐波的谐波阶数n为2时，所述非线性谐波的几何相位Φ＝±mθ或Φ＝±3mθ，其中，θ为所述超构单元(4)的旋转角；
所述非线性谐波的谐波阶数n为3时，所述非线性谐波的几何相位Φ＝±2mθ、±mθ或0。
5.根据权利要求3所述的基于高阶旋转对称单元的非线性几何相位超表面，其特征在于，所述具有局域高阶旋转对称性的超构单元(4)的局域旋转对称阶数m大于4。
6.根据权利要求3所述的基于高阶旋转对称单元的非线性几何相位超表面，其特征在于，所述超构单元(4)呈周期性排列包括采用四方点阵或六角点阵的排布形成全局旋转对称阶数q为4或6的晶格结构。
7.根据权利要求1所述的基于高阶旋转对称单元的非线性几何相位超表面，其特征在于，所述超构单元(4)为亚波长金属结构，所述亚波长金属结构的材料包括金、铝中的一种。
8.根据权利要求7所述的基于高阶旋转对称单元的非线性几何相位超表面，其特征在于，所述亚波长金属结构包括正m边形、正m星形，其中，m为所述超构单元(4)的局域旋转对称阶数。
9.根据权利要求1所述的基于高阶旋转对称单元的非线性几何相位超表面，其特征在于，所述非线性材料层(2)包括氧化铟锡、铌酸锂中的一种；所述衬底(1)包括玻璃、蓝宝石中的一种。
10.根据权利要求1所述的基于高阶旋转对称单元的非线性几何相位超表面，其特征在于，所述超构单元(4)的厚度h＜λ0/2、周期p＜λ0；所述非线性材料层(2)的厚度t＜λ0/2；所述衬底(1)的厚度为d＞＞λ0，其中，λ0为谐波波长。
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/CN2022/079528 WO2023168555A1 (zh) | 2022-03-07 | 2022-03-07 | 基于高阶旋转对称单元的非线性几何相位超表面 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114762191A true CN114762191A (zh) | 2022-07-15 |
| CN114762191B CN114762191B (zh) | 2023-05-26 |
Family
ID=82336685
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202280000514.XA Active CN114762191B (zh) | 2022-03-07 | 2022-03-07 | 基于高阶旋转对称单元的非线性几何相位超表面 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN114762191B (zh) |
| WO (1) | WO2023168555A1 (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117075409A (zh) * | 2023-10-16 | 2023-11-17 | 安徽大学 | 一种可增强二次谐波产生效率的bic超表面 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102683866A (zh) * | 2011-03-15 | 2012-09-19 | 深圳光启高等理工研究院 | 一种控制电磁波传播方向的装置 |
| JP2017092178A (ja) * | 2015-11-06 | 2017-05-25 | 株式会社豊田中央研究所 | 電磁波装置 |
| CN108107498A (zh) * | 2017-12-15 | 2018-06-01 | 中国科学院光电技术研究所 | 一种基于旋转晶向的非线性位相梯度超表面 |
| CN109683334A (zh) * | 2019-01-16 | 2019-04-26 | 南方科技大学 | 一种透射式光分束器及其制造方法 |
| US20220011646A1 (en) * | 2018-10-01 | 2022-01-13 | William Marsh Rice University | Harmonic light-generating metasurface |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113946009B (zh) * | 2021-10-27 | 2023-10-03 | 中国科学院光电技术研究所 | 一种超构表面器件及其设计、制备方法 |
-
2022
- 2022-03-07 CN CN202280000514.XA patent/CN114762191B/zh active Active
- 2022-03-07 WO PCT/CN2022/079528 patent/WO2023168555A1/zh unknown
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102683866A (zh) * | 2011-03-15 | 2012-09-19 | 深圳光启高等理工研究院 | 一种控制电磁波传播方向的装置 |
| JP2017092178A (ja) * | 2015-11-06 | 2017-05-25 | 株式会社豊田中央研究所 | 電磁波装置 |
| CN108107498A (zh) * | 2017-12-15 | 2018-06-01 | 中国科学院光电技术研究所 | 一种基于旋转晶向的非线性位相梯度超表面 |
| US20220011646A1 (en) * | 2018-10-01 | 2022-01-13 | William Marsh Rice University | Harmonic light-generating metasurface |
| CN109683334A (zh) * | 2019-01-16 | 2019-04-26 | 南方科技大学 | 一种透射式光分束器及其制造方法 |
Non-Patent Citations (2)
| Title |
|---|
| XIN XIE: "Generalized Pancharatnam-Berry Phase in Rotationally Symmetric Meta-Atoms" * |
| 赵文超: "基于Au-WS2复合超表面的高效率二次谐波微纳功能器件研究" * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117075409A (zh) * | 2023-10-16 | 2023-11-17 | 安徽大学 | 一种可增强二次谐波产生效率的bic超表面 |
| CN117075409B (zh) * | 2023-10-16 | 2023-12-26 | 安徽大学 | 一种可增强二次谐波产生效率的bic超表面 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114762191B (zh) | 2023-05-26 |
| WO2023168555A1 (zh) | 2023-09-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Li et al. | Tailoring terahertz wavefront with state switching in VO2 Pancharatnam–Berry metasurfaces | |
| Zhao et al. | Conjugated gammadion chiral metamaterial with uniaxial optical activity and negative refractive index | |
| WO2017020792A1 (zh) | 一种左旋圆偏振转换的超材料薄膜 | |
| CN114089539B (zh) | 一种基于复合相位调控的双完美涡旋光束超表面设计方法 | |
| Nakata et al. | Reconfigurable terahertz quarter-wave plate for helicity switching based on babinet inversion of an anisotropic checkerboard metasurface | |
| CN114762191A (zh) | 基于高阶旋转对称单元的非线性几何相位超表面 | |
| Kwon et al. | Transmission magnitude and phase control for polarization-preserving reflectionless metasurfaces | |
| Wang et al. | Tunable annular acoustic metasurface for transmitted wavefront modulation | |
| CN109901257A (zh) | 一种可见光超材料偏振转换器 | |
| Dong et al. | Dynamic chiroptical responses in transmissive metamaterial using phase-change material | |
| Wu et al. | Near-infrared chirality of plasmonic metasurfaces with gold rectangular holes | |
| Liu et al. | Boosting second harmonic generation by merging bound states in the continuum | |
| CN117148603A (zh) | 基于二氧化钒-石墨烯的太赫兹波前调控器件 | |
| Yang et al. | Polarization conversion in terahertz planar metamaterial composed of split-ring resonators | |
| Luo et al. | Chiral metasurface design with highly efficient and controllable asymmetric transmission and perfect polarization conversion of linearly polarized electromagnetic waves in the THz range | |
| Amin et al. | Metasurface supporting broadband circular dichroism for reflected and transmitted fields simultaneously | |
| Li et al. | Few-layer metasurfaces with arbitrary scattering properties | |
| Zhang et al. | A vortex-focused beam metalens array in the visible light range based on computer-generated holography | |
| Dong et al. | Full-space terahertz metasurface based on thermally tunable InSb | |
| Guo et al. | Functional metasurfaces based on metallic and dielectric subwavelength slits and stripes array | |
| Li et al. | Reflection-type transparent metamaterial polarization rotator with ultra-wide bandwidth | |
| Wang et al. | Broadband asymmetric transmission by rotated bilayer cross-shaped metamaterials | |
| Pourali et al. | Tunable enhanced Faraday rotation in a defected plasma photonic crystal under external magnetic field with different declinations | |
| Nemati et al. | Electrically tunable polarization-insensitive MIM plasmonic metasurface operating in transmission mode | |
| Liu et al. | Theoretical design of multilayered VO2-based switchable terahertz metasurfaces as broadband functional devices |
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 |