基于连续域束缚态的可调谐手性响应超表面

,  ,  , 

上一篇

下一篇

吕淑媛, 张璐璐, 罗文峰, 安岚. 基于连续域束缚态的可调谐手性响应超表面[J]. 物理学报, 2026, 75(9): 090404. doi: 10.7498/aps.75.20251726
引用本文: 吕淑媛, 张璐璐, 罗文峰, 安岚. 基于连续域束缚态的可调谐手性响应超表面[J]. 物理学报, 2026, 75(9): 090404. doi: 10.7498/aps.75.20251726
Shuyuan LV, Lulu ZHANG, Wenfeng LUO, Lan AN. Tunable chiral response metasurface based on bound states in the continuum[J]. Acta Physica Sinica, 2026, 75(9): 090404. doi: 10.7498/aps.75.20251726
Citation: Shuyuan LV, Lulu ZHANG, Wenfeng LUO, Lan AN. Tunable chiral response metasurface based on bound states in the continuum[J]. Acta Physica Sinica, 2026, 75(9): 090404. doi: 10.7498/aps.75.20251726

基于连续域束缚态的可调谐手性响应超表面

    通讯作者: E-mail:1159955131@qq.com.; 

Tunable chiral response metasurface based on bound states in the continuum

    Corresponding author: E-mail:1159955131@qq.com.; 
  • 摘要: 具有高品质因子(Q)共振与可调谐手性响应的光学超表面是当前纳米光子学领域的研究热点, 对推动手性光学器件提升具有重要意义. 本研究提出了一种基于连续域束缚态(BIC)的手性锗超表面设计, 通过打破C2对称性, 将BIC转变为准BIC(q-BIC), 从而激发支持手性响应的共振模式, 并通过改变不对称参数对共振特性与手性响应实现精确调控, 在近红外波段实现了超高Q因子(Q = 6121.14)与强圆二色性(CD = –0.94). 此外, 通过调节结构参数, 可实现对手性响应的静态调控. 最后, 利用石墨烯的可调谐特性, 在不改变结构参数的情况下通过调控其费米能级, 实现了手性响应的动态调控. 本研究为可控平面光学手性的设计提供了一种可行路径, 也为手性光子器件提供了新的设计思路.
  • 加载中
  • 图 1  超表面结构示意图 (a)由锗超表面和布拉格反射镜组成的结构示意图, 超表面上方环境折射率为1.34, 以模拟其在水溶液; (b)结构的单元胞顶视图, 结构参数为la = 148 nm, lb = la + δ, wa = wb = 124 nm, d = 1500 nm, d1 = 85 nm, d2 = 204 nm, Δd = 92.25 nm

    Figure 1.  Schematic of the metasurface structure: (a) Metasurface consisting of germanium and Bragg reflector; the refractive index of the environment above the metasurface is 1.34, to simulate its environment in an aqueous solution; (b) top view of the unit cell structure, where the structural parameters la = 148 nm, lb = la + δ, wa = wb = 124 nm, d = 1500 nm, d1 = 85 nm, d2 = 204 nm, Δd = 92.25 nm.

    图 2  手性超表面的本征模态分析 (a), (b) Γ点附近超表面的Q因子与能带结构; (c)—(e) δ = 0 nm, δ = 20 nm, δ = 35 nm时的远场辐射极化矢量分布图, 其中红色与蓝色椭圆分别代表RCP和LCP圆极化状态

    Figure 2.  Eigenmodal analysis of chiral metasurfaces: (a), (b) The Q-factor and band structure of the metasurface near the Γ point, respectively; (c)–(e) far-field radiation polarization vector distributions for δ = 0 nm, δ = 20 nm, and δ = 35 nm, where the red and blue ellipses represent RCP and LCP circular polarization states, respectively.

    图 3  不同结构参数对CD光谱的影响, 即锗超表面随着(a) H, (b) P, (c) l变化时的CD谱线

    Figure 3.  Effects of different structural parameters on the CD spectra: The CD spectral lines of germanium metasurface with the parameters (a) H, (b) P, (c) l.

    图 4  手性q-BIC的光学特性分析 (a)随着不对称参数δ变化的CD谱线; (b) δ = 35 nm时吸收谱线和CD谱线; (c) Q因子随不对称参数δ的变化趋势; (d) δ = 35 nm时双凹槽中心间距变化时的吸收谱线; (e) CD光谱随Δd变化的演化图

    Figure 4.  Analysis of the optical properties of the chiral q-BIC: (a) CD spectral lines as a function of the asymmetry parameter δ; (b) absorption spectra and CD spectral lines for δ = 35 nm; (c) variation trend of the Q-factor as a function of the asymmetry parameter δ; (d) absorption spectra with varying center-to-center spacing of the double grooves for δ = 35 nm; (e) evolution of the CD spectra with Δd.

    图 5  q-BIC模式的散射功率 (a), (b) RCP和LCP入射下手性q-BIC模式不同多极矩的散射功率; (c) MD在x, y和z三个方向上的散射功率; (d), (e) RCP和LCP入射下x-y平面下手性q-BIC模式的磁场分布

    Figure 5.  Scattering power of the q-BIC modes: (a), (b) Scattering power of different multipole moments for the chiral q-BIC mode under RCP and LCP illumination, respectively; (c) scattering power of the MD along the x, y, and z-directions; (d), (e) magnetic field distributions of the chiral q-BIC mode in the x-y plane under RCP and LCP illumination, respectively.

    图 6  石墨烯的手性响应调控 (a)引入单层石墨烯的结构示意图; (b), (c)不同费米能级下石墨烯介电常数实部与虚部随波长的变化; (d) RCP和LCP的吸收谱线; (e)不同费米能级EF下结构的CD谱线

    Figure 6.  Chiral response modulation via graphene: (a) Schematic of the structure with a monolayer graphene introduced; (b), (c) variation of the real and imaginary parts of graphene permittivity with wavelength at different Fermi energy levels; (d) absorption spectra under RCP and LCP illumination; (e) CD spectra of the structure at different Fermi levels EF.

    表 1  与现有结构的性能比较

    Table 1.  Performance comparison with existing structures.

    文献 波段 物理机制 光谱类型 Q因子 |CD|
    [31] 近红外 对称保护型BIC 透射 3216 0.98
    [44] 太赫兹 对称保护型BIC 透射 2274
    3736
    0.996
    0.395
    [45] 近红外 对称保护型BIC 透射 350 0.89
    [46] 近红外 偶然型BIC 透射 4423 0.99
    [47] 近红外 对称保护型BIC 透射 1517.9 0.979
    [48] 可见光 对称保护型BIC 反射 2663 0.93
    本工作 近红外 对称保护型BIC 吸收 6121.14 0.94
    下载: 导出CSV
  • [1] Chen Y, Gao J, Yang X D 2018 Nano Lett. 18 520 doi: 10.1021/acs.nanolett.7b04515
    [2] Moss G P 1996 Pure Appl. Chem. 68 2193 doi: 10.1351/pac199668122193
    [3] Menzel C, Rockstuhl C, Lederer F 2010 Phys. Rev. A 82 053811 doi: 10.1103/PhysRevA.82.053811
    [4] Chela-Flores J 1991 Chirality 3 389 doi: 10.1002/chir.530030503
    [5] Hentschel M, Schäferling M, Duan X Y, Giessen H, Liu N 2017 Sci. Adv. 3 e1602735 doi: 10.1126/sciadv.1602735
    [6] Tang Y Q, Cohen A E 2011 Science 332 333 doi: 10.1126/science.1202817
    [7] Zhu Y Y, Xu L G, Ma W, Zhou X, Kuang H, Wang L B, Xu C L 2012 Chem. Commun. 48 11889 doi: 10.1039/c2cc36559f
    [8] Yeom J, Guimaraes P P G, Ahn H M, Jung B K, Hu Q, McHugh K, Mitchell M J, Yun C O, Langer R, Jaklenec A 2020 Adv. Mater. 32 1903878 doi: 10.1002/adma.201903878
    [9] Omatsu T, Miyamoto K, Toyoda K, Morita R, Arita Y, Dholakia K 2019 Adv. Opt. Mater. 7 1801672 doi: 10.1002/adom.201801672
    [10] Wade J, Hilfiker J N, Brandt J R, et al. 2020 Nat. Commun. 11 6137 doi: 10.1038/s41467-020-19951-y
    [11] Ouyang L, Rosenmann D, Czaplewski D A, Gao J, Yang X 2020 Nanotechnology 31 295203 doi: 10.1088/1361-6528/ab88ea
    [12] Chen Y, Du W, Zhang Q, Ávalos-Ovando O, Wang J, Xu Q H, Liu N, Okamoto H, Govorov A O, Xiong Q H, Qiu C W 2022 Nat. Rev. Phys. 4 113 doi: 10.1038/s42254-021-00391-6
    [13] Li M H, Guo L Y, Dong J F, Yang H L 2014 J. Phys. D: Appl. Phys. 47 185102 doi: 10.1088/0022-3727/47/18/185102
    [14] Zhang Y D, Yu S, Han B, Zhou Y L, Zhang X W, Gao X Q, Tang Z Y 2022 Matter 5 837 doi: 10.1016/j.matt.2022.01.001
    [15] Lininger A, Palermo G, Guglielmelli A, Nicoletta G, Goel M, Hinczewski M, Strangi G 2023 Adv. Mater. 35 2107325 doi: 10.1002/adma.202107325
    [16] Zhou J F, Dong J F, Wang B N, Koschny T, Kafesaki M, Soukoulis C M 2009 Phys. Rev. B 79 121104 doi: 10.1103/PhysRevB.79.121104
    [17] Zhang S, Zhou J F, Park Y S, Rho J, Singh R, Nam S, Azad A K, Chen H T, Yin X, Taylor A J, Zhang X 2012 Nat. Commun. 3 942 doi: 10.1038/ncomms1908
    [18] Hu G W, Wang M S, Mazor Y, Qiu C W, Alù A 2021 Trends Chem. 3 342 doi: 10.1016/j.trechm.2021.02.004
    [19] Tang Y H, Liang Y, Yao J, Chen M K, Lin S R, Wang Z, Zhang J C, Huang X G, Yu C Y, Tsai D P 2023 Laser Photonics Rev. 17 2200597 doi: 10.1002/lpor.202200597
    [20] Shi T, Deng Z L, Geng G Z, Zeng X Z, Zeng Y X, Hu G W, Overvig A, Li J J, Qiu C W, Alu A, Kivshar Y S, Li X P 2022 Nat. Commun. 13 4111 doi: 10.1038/s41467-022-31877-1
    [21] Hsu C W, Zhen B, Stone A D, Joannopoulos J D, Soljačić M 2016 Nat. Rev. Mater. 1 16048 doi: 10.1038/natrevmats.2016.48
    [22] Friedrich H, Wintgen D 1985 Phys. Rev. A 32 3231 doi: 10.1103/PhysRevA.32.3231
    [23] Marinica D C, Borisov A G, Shabanov S V 2008 Phys. Rev. Lett. 100 183902 doi: 10.1103/PhysRevLett.100.183902
    [24] Liu Z J, Xu Y, Lin Y, Xiang J, Feng T H, Cao Q T, Li J T, Lan S, Liu J 2019 Phys. Rev. Lett. 123 253901 doi: 10.1103/PhysRevLett.123.253901
    [25] Overvig A, Yu N, Alù A 2021 Phys. Rev. Lett. 126 073001 doi: 10.1103/PhysRevLett.126.073001
    [26] Tang R J, Shi Y L, Shang H P, Wu J H, Ma H, Wei M L, Luo Y, Chen Z Q, Ye Y T, Jian J L, Zheng X R, Lin H T, Li L 2023 Nanophotonics 12 3257 doi: 10.1515/nanoph-2023-0156
    [27] Zhang Z Y, Bulgakov E, Pichugin K, Sadreev A, Xu Y, Qin Y W 2023 Phys. Rev. Appl. 20 L011003 doi: 10.1103/PhysRevApplied.20.L011003
    [28] Wang Y L, Han Z H, Du Y, Qin J Y 2021 Nanophotonics 10 1295 doi: 10.1515/nanoph-2020-0582
    [29] Koshelev K, Lepeshov S, Liu M K, Bogdanov A, Kivshar Y 2018 Phys. Rev. Lett. 121 193903 doi: 10.1103/PhysRevLett.121.193903
    [30] Li S, Wu Z Y, Liu X S, Ye Q, Liu G Q, Chen J, Du W, Tang C J, Liu Z Q 2025 Laser Photonics Rev. 19 e00477 doi: 10.1002/lpor.202500477
    [31] Li J B, Duan Q L, Dong X, Liu T, Zhou H, Liu Y, Li C, Wang Z, Zhang Y 2024 Phys. Rev. Appl. 22 044035 doi: 10.1103/PhysRevApplied.22.044035
    [32] Lv S Y, Hu F, Luo W F, Xu H F, An L 2024 Opt. Express 32 30053 doi: 10.1364/OE.530418
    [33] Huang Z S, Wang J, Jia W, Zhang S W, Zhou C H 2024 Opt. Express 32 33029 doi: 10.1364/OE.533165
    [34] Zeng J H, Zhou Y, Fu X, Yang J R, Chen Y X, Hong W Y 2025 Opt. Express 33 27014 doi: 10.1364/OE.566464
    [35] Shen Y, Xiao Y F, Yang L L, Wu L, Cai B, Luo H, Cheng Y Z, Li X C 2025 Results Eng. 28 107432 doi: 10.1016/j.rineng.2025.107432
    [36] Li H J, Zhou H M, Wei G G, Xu H S, Qin M, Liu J Q, Wu F 2023 Nanoscale 15 6636 doi: 10.1039/D3NR00055A
    [37] Liu Q K, Li Y, Lu Z D, Zhou Y J, Liu W M, Luo X Q, Wang X L 2023 Phys. Rev. B 108 155410 doi: 10.1103/PhysRevB.108.155410
    [38] Zong S, Chen J F, Cai S J, Liu X S, Liu G Q, Chen J, Tang C J, Du W, Liu Z Q 2024 Appl. Phys. Lett. 124 031702 doi: 10.1063/5.0180635
    [39] Zhong H Z, Liu Z Q, Liu X S, Fu G L, Liu G Q, Chen J, Tang C J 2020 Opt. Express 28 37294 doi: 10.1364/OE.412861
    [40] Masoudian Saadabad R, Huang L J, Miroshnichenko A E 2021 Phys. Rev. B 104 235405 doi: 10.1103/PhysRevB.104.235405
    [41] Palik E D 1985 Handbook of Optical Constants of Solids (New York: Academic Press
    [42] Koshelev K, Favraud G, Bogdanov A, Kivshar Y, Fratalocchi A 2019 Nanophotonics 8 725 doi: 10.1515/nanoph-2019-0024
    [43] Kong X, Khorashad L K, Wang Z, Govorov A O 2018 Nano Lett. 18 2001 doi: 10.1021/acs.nanolett.7b05446
    [44] Ma T, Tian J K, Sang W, Wang D D, Li J 2024 J. Opt. 26 075101 doi: 10.1088/2040-8986/ad44ab
    [45] Kim K H, Kim J R 2021 Adv. Opt. Mater. 9 2101162 doi: 10.1002/adom.202101162
    [46] Li Z Y, Liu C, Wang Y X, Yao Q H, Peng Y X, Wang K J, Liu J Q, Lu D Y, He M D 2025 Phys. Lett. A 560 130965 doi: 10.1016/j.physleta.2025.130965
    [47] Wang P P, Gao Z, He Y, Xu Z G, Zhao T G 2024 IEEE Photonics J. 16 1 doi: 10.1109/JPHOT.2024.3401030
    [48] Chen Y, Deng H C, Sha X B, Chen W J, Wang R Z, Chen Y H, Wu D, Chu J R, Kivshar Y S, Xiao S M, Qiu C W 2023 Nature 613 474 doi: 10.1038/s41586-022-05467-6
    [49] Savinov V, Fedotov V A, Zheludev N I 2014 Phys. Rev. B 89 205112 doi: 10.1103/PhysRevB.89.205112
    [50] Kaelberer T, Fedotov V A, Papasimakis N, Tsai D P, Zheludev N I 2010 Science 330 1510 doi: 10.1126/science.1197172
    [51] Zhang J F, Zhu Z H, Liu W, Yuan X D, Qin S Q 2015 Nanoscale 7 13530 doi: 10.1039/C5NR03060A
    [52] Xiao S Y, Wang T, Liu Y B, Xu C, Han X, Yan X C 2016 Phys. Chem. Chem. Phys. 18 26661 doi: 10.1039/C6CP03731C
  • 加载中
图( 6) 表( 1)
计量
  • 文章访问数:  464
  • HTML全文浏览数:  464
  • PDF下载数:  11
  • 施引文献:  0
出版历程
  • 收稿日期:  2025-12-16
  • 刊出日期:  2026-05-05

基于连续域束缚态的可调谐手性响应超表面

    通讯作者: E-mail:1159955131@qq.com.; 
  • 西安邮电大学电子工程学院, 西安 710121

摘要: 具有高品质因子(Q)共振与可调谐手性响应的光学超表面是当前纳米光子学领域的研究热点, 对推动手性光学器件提升具有重要意义. 本研究提出了一种基于连续域束缚态(BIC)的手性锗超表面设计, 通过打破C2对称性, 将BIC转变为准BIC(q-BIC), 从而激发支持手性响应的共振模式, 并通过改变不对称参数对共振特性与手性响应实现精确调控, 在近红外波段实现了超高Q因子(Q = 6121.14)与强圆二色性(CD = –0.94). 此外, 通过调节结构参数, 可实现对手性响应的静态调控. 最后, 利用石墨烯的可调谐特性, 在不改变结构参数的情况下通过调控其费米能级, 实现了手性响应的动态调控. 本研究为可控平面光学手性的设计提供了一种可行路径, 也为手性光子器件提供了新的设计思路.

English Abstract

    • 手性是自然界普遍存在的一种基本几何属性, 描述物体无法通过旋转或平移操作使其与自身的镜像完全重合的性质[1–3]. 该性质广泛存在于从微观分子到宏观天体的各个尺度的客体上, 有机化合物[4]、蛋白质分子[5]、星系自转等都呈现出手性特征[6], 为科学研究提供了丰富的素材, 在多个学科领域的研究均具有重要意义[7–9]. 在光学领域, 手性现象的研究和应用价值尤为突出. 手性物质与手性偏振光相互作用时会产生手性光学效应现象, 可用旋光性[10]或圆二色性(circular dichroism, CD)[11]进行描述, 表现为对左旋圆偏振光(left cirularly polarized, LCP)和右旋圆偏振光(right cirularly polarized, RCP)的差异响应. 由于天然手性材料的手性吸收较弱, 需要足够的光程积累才能产生可观测的信号, 限制了其在偏振成像、生物传感等方面的应用. 因此, 研究者致力于发展人工手性超材料与超表面, 为操控手性光学响应提供了全新途径[12,13], 通过在亚波长尺度上精确调控电磁响应, 增强光与物质的相互作用, 为手性光学的发展提供新的机遇[14,15]. 目前, 基于等离子体的金属手性超表面虽能产生显著的光学手性, 但伴随着不可避免的金属损耗与热效应问题, 且激发的模式多为宽带模式[16–18], 在实际应用中受到一定的限制. 因此, 实现同时具备高品质因子(Q)和可调谐手性响应的光学超表面吸收体成为一个亟待解决的挑战.

      近年来, 连续域束缚态(bound states in the continuum, BIC)的研究为超表面的手性响应提供了新的解决方案[19,20]. BIC是一种独特的光学共振态, 虽处于辐射态的连续域中, 却能保持完美的局域特性[21–23]. 作为连续域中的奇异点, 理想BIC具有无限高的Q因子和极强的光场局域能力, 可显著增强光与物质的相互作用. BIC已广泛用于非线性[24,25]、激光[26,27]、传感[28]等领域. 在实际应用中, 理想的BIC通常会因对称性破缺而与辐射场耦合, 转化为具有类似超腔模式特征的准BIC(q-BIC)模式[29], 表现出超高Q因子共振. 具体而言, 可通过破坏结构对称性, 使得q-BIC模式与辐射场耦合, 并支持圆偏振光选择性共振, 进而在动量空间中形成手性点(C点). 高Q因子保证了共振模式的强场局域能力, 能够显著增强光与物质的相互作用, 从而提升CD响应. 因此, q-BIC机制通过结合高Q因子与对称性调控, 为实现强手性响应提供了物理基础.

      在手性光学响应的多种表现形式中, 手性吸收因其在偏振选择性光热转换[30]、控制手性光子器件[31]等应用中的潜力而备受关注. 通过引入动态可调材料(如石墨烯、二氧化钒(VO2)等)或改变结构参数, 可进一步实现对手性吸收响应的主动调控[32–35]. 2023年, Li等[36]通过引入金属基底构建了无传输通道的单端口系统超表面, 借助单端口临界耦合机制, 使手性超表面实现了对一种圆偏振光的完美吸收, 同时有效抑制了另一种圆偏振光的吸收, CD接近0.812. 同年, Liu等[37]将平面手性硅超表面与VO2结合形成 Si-VO2杂化超表面, 利用VO2的相变特性来主动调节耗散损耗, 进而实现对光学手性的精确控制. 2024年Zong等[38]利用布拉格反射镜来增强光吸收, 在垂直入射条件下, 在1185.35 nm波长处CD达到最大值0.91.

      本研究提出一种全介电手性超表面, 为增强锗超表面吸收能力, 借助布拉格反射镜, 利用其全反射机制增强光吸收效果[39], 实现超表面的超高Q因子和高CD信号强度. 超表面与布拉格反射镜之间设置适当间距, 以确保反射光不会激发结构中的导模共振[40]. 当间距小于某一临界值时, 超表面与布拉格反射镜之间会产生强耦合, 激发导模共振, 导致吸收谱出现分裂峰并降低手性响应, 随着间距d增大, 这种耦合效应减弱. 当d 增大至1300 nm以上, 反射光与超表面之间的耦合被有效抑制, 导模共振消失, 吸收谱仅表现为由超表面自身支持的单一q-BIC手性共振峰. 选取d = 1500 nm, 使系统处于弱耦合状态, 此时布拉格反射镜仅用于构建单端口系统并阻断透射, 不引入额外模式, 从而在保持高手性响应与高吸收率的同时避免了导模共振的干扰. 进一步探讨了在垂直入射下, 如何通过调节结构的面内C2对称性, 实现圆偏振光诱导的手性吸收, 观测到的CD值最高可达–0.94. 此外, 通过调节结构参数, 可实现对手性响应的静态调控, 最后, 利用石墨烯的可调谐特性, 在不改变结构参数的情况下实现了手性响应的动态调制. 本研究不仅丰富了手性超表面的设计, 还为光子器件提供了新的设计思路.

    • 图1(a)展示了所设计的手性超表面结构, 由带双凹槽的方形锗纳米柱周期性阵列组成, 其位于带有布拉格反射镜的反射腔上, 布拉格反射镜由硅(Si)与二氧化硅(SiO2)交替堆叠构成, 其厚度分别为d1和d2, 反射腔厚度为d. 单元胞顶视图如图1(b)所示, 周期P = 640 nm, 方形锗纳米柱的边长l = 350 nm, 顶部和底部设计两个尺寸相同的凹槽, 且沿x轴反向平移. 凹槽的宽度为wa = wb, 一个凹槽长度为la, 另一个凹槽长度为lb = la + δ. Δd表示两个凹槽中心点之间的距离, 当Δd ≠ 0时, 该超表面具有C2对称性但不存在面内镜像对称性. δ为不对称参数, 当δ ≠ 0时, 面内C2对称性被打破. 这些纳米柱采用高折射率锗材料(n = 4.27, nk = 8.5×10–4)[41]制成, 其厚度H = 210 nm, Si和SiO2的折射率分别为nH = 3.48和nL = 1.45. 在圆偏振光沿–z方向垂直入射的条件下, 本研究利用COMSOL Multiphysics软件, 对所提出手性超表面的光学特性进行数值模拟研究. 在x和y方向设置周期性边界条件, 以确保物理场在边界上的分布是连续的. 同时, 在z方向的两个界面处采用完美匹配层(PML)作为边界条件, 有效吸收入射波能量, 降低反射率, 模拟无限大的空间, 确保模拟结果精确和高效.

    • 在满足C2旋转对称性条件的前提下, 将重点分析TE模式, 因为其在特定波长范围内比TM模式具有更高的Q因子, 针对TE模式研究了手性超表面中BIC的物理特性. 图2(a), (b)展示了第一布里渊区沿M-Г-X方向分布的Q因子及其能带结构. 图2(a)插图展示了方形晶格的第一布里渊区, 可以观察到Q因子在Γ点处趋于无穷大, 呈现出理想对称保护性BIC的特征[42]. 当δ = 0 nm, Δd = 92.25 nm时, 镜像对称性被打破, BIC在Γ点处表现为一个特殊的涡旋偏振奇点(V点, 以黑点标注), 该奇点被具有非零螺旋度的椭圆偏振态所环绕, 如图2(c)所示. 在镜面对称性破缺的前提下, 进一步引入不对称参数, 可实现超表面的本征手性响应. 如图2(d), (e)所示, 通过改变凹槽长度lb, 引入不对称参数δ, 面内C2对称性被打破. 当δ = 20 nm, Δd = 92.25 nm时, V点可分裂为一对圆偏振态, 其中以红点标注的右旋C+点的位置, 如图2(d)所示. 进一步将δ增至35 nm时, 如图2(e)所示, C+点移动到Γ点, 从而形成本征手性q-BIC. 综上所述, 镜面对称性的破缺是构建非简并手性模式的基础, 而C2对称性的打破则实现了手性响应在动量空间的调控. 因此, 通过在超表面中构建具有可调k空间C点的手性q-BIC, 是实现高Q因子与可调谐的手性光学响应的关键.

      为了研究超表面对圆偏振光的不同光学响应, 本研究采用圆二色性(CD)进行表征, 并分析了结构参数对CD光谱的影响. 在其他参数与图1保持一致的情况下, 分别计算了锗超表面厚度H在200—220 nm、周期P在630—650 nm、锗纳米柱边长l在346—354 nm变化时的CD光谱, 结果分别如图3(a)—(c)所示. 其中CD的表达式定义为[43]

      其中ALCP和ARCP分别是该结构在LCP和RCP入射下的吸收.

      如图3(a)所示, 随着H的增大, CD谱线的峰值呈现先增后减趋势, 当H = 210 nm时, 在1186 nm波长处出现最小值–0.94, 同时CD谱线的共振波长发生显著红移. 考虑到设计手性吸收超表面的要求, 选用H = 210 nm作最优厚度参数. 如图3(b)所示, 随着P的增大, CD谱线的峰值无明显变化, 但CD谱线的共振波长发生红移. 表明周期P对共振波长存在弱调控作用. 如图3(c)所示, 随着l的增大, CD谱线的峰值呈现先增后减的变化规律, 当l = 350 nm时, CD信号强度达到最大值, 该值被选为边长l的最优取值. 此外, 本研究通过系统调控结构参数H, P, l, 实现了在确定波段内对CD响应的优化, 为后续研究高性能手性吸收奠定了可靠的参数基础.

      在∆d = 92.25 nm的条件下, 通过引入不对称参数δ, 打破了结构的面内C2对称性, 可使理想BIC转变为手性q-BIC. 为了探究手性q-BIC的响应特性, 我们分析了不对称参数δ对超表面光学响应的影响. 图4(a)展示了随不对称参数δ变化的CD谱线. 随着δ的增大, CD谱线的共振波长发生蓝移且谱线线宽逐渐展宽的同时其峰值呈波动性变化, 在δ = 35 nm处CD谱线的峰值达到最小值. 为了更直观地展示手性q-BIC超表面的光学特性, 图4(b)展示了当δ = 35 nm时, 在垂直入射条件下, 圆偏振光的吸收光谱和CD光谱. 在1186 nm波长处, RCP入射下吸收率高达0.97, 而LCP入射下吸收率仅为0.026; 同时可观察到, 在1186 nm的共振波长处产生了尖锐且最大值为–0.94的CD值. 该吸收峰的形成源于结构在共振波长处满足临界耦合状态, q-BIC模式的辐射损耗率与锗材料的本征吸收损耗率相互匹配, 使得反射波相干抵消, 从而实现近完美的吸收. 此外, 图4(c)展示了Q因子对不对称参数δ的依赖关系. 随着δ从0增至40 nm, 几何对称性被打破, 观察到Q因子从107量级下降到103量级. 在30 nm < δ < 40 nm的范围内, 保证高CD值的同时, Q因子可被调谐到9043—4322范围内. 上述结果证实了图2(d), (e)中改变δ可调节C2对称性破缺程度, 从而实现手性模式在动量空间中的演化. 此外, 表1通过系统比较本研究与近期文献在波段、Q因子和CD值等关键指标上的表现, 清晰展示了本结构在近红外波段兼顾高Q值与强CD响应的综合优势, 为实现二者的协同优化提供了有效途径, 也为发展高性能可调手性器件奠定了新思路.

      我们进一步研究了参数∆d对手性响应的调控作用. 在保持δ = 35 nm的条件下, 图4(d)展示了当∆d分别为0 nm, ±30 nm, ±60 nm和±92 nm时, LCP和RCP的吸收谱线. 如图4(d)所示, 当∆d = 0 nm时, 结构具备对称性, LCP和RCP吸收谱完全重合, 未表现出手性. 随着|∆d|逐渐增大, 镜像对称性被打破, 吸收谱出现显著的不对称性; 当∆d符号反转时, 可以观察到相反的吸收选择特性. CD谱随∆d的变化关系, 如图4(e)所示. 可以看出, 当∆d的符号反转时, 手性响应也会反转; 同时∆d符号的反转虽引起CD符号的反转, 却不影响CD的绝对值大小, 说明∆d在不影响共振模式本征特性的基础上, 是实现手性翻转的关键参数. 综上所述, 本超表面能够在确定共振波长附近实现接近理论极限的CD响应以及灵活的手性切换, 为手性光子器件的设计提供了有效途径.

      为了进一步探索q-BIC模式的机理, 在δ = 35 nm, Δd = 92.25 nm的条件下, 对笛卡尔坐标系中不同多极子的散射功率进行分析. 这些多极子包括电偶极子(ED)、磁偶极子(MD)、环形偶极子(TD)、电四极子(EQ)以及磁四极子(MQ), 相应的米氏共振公式如下[49,50]:

      式中, c是真空中的光速; ω是入射光的角频率; r是笛卡尔坐标系中从原点到点(x, y, z)的距离矢量; α, β是x, y, z三个坐标轴的方向; j是手性锗超表面单位晶胞内的位移电流密度. 由多极子激发产生的远场散射能量通过下式进行计算:

      基于多极散射理论, 在RCP波和LCP波的激励下, 在超表面中诱导的不同多极散射功率分别如图5(a), (b)所示. 对于RCP波, 从图5(a)可以看出, MD的贡献占主导地位, 这种主导性源于MD模式涡旋状电流分布特征, 能够产生高度局域化的磁场(如图5(d)). 因此, 该q-BIC共振主要由MD诱导产生, 强烈的MD共振产生了如图4(b)所示的尖锐吸收谱. 如图5(b)所示, 在LCP入射下, 多极子贡献在1186 nm处被抑制, 从而阻止了q-BIC共振的激发. 以上结果进一步证实, 该超表面的本征手性响应源于手性q-BIC共振, 且主要由MD模式主导. 图5(c)展示了分解后的MD散射功率分量. 可以看出MD的x和y分量显著受到抑制, 而MD的z分量则被共振激发, 其值几乎与MD散射功率相等, 证实了q-BIC中MD的z方向特性. 为研究MD产生的手性响应, 我们模拟了1186 nm谐振波长下z = H/2截面处RCP和LCP波入射下的磁场分布, 结果分别如图5(d), (e)所示, 其中白色箭头表示位移电流密度. 图5(d)表明, 在RCP激发下, 该结构的位移电流呈现反向涡旋特征, 是MD共振的典型特征. 相比之下, 如图5(e)所示, 在LCP激发下, 该结构仅表现出弱的电磁响应. 这一显著对比再次表明MD模式在LCP与RCP激发下的选择性耦合. 正是这种差异化耦合, 导致了工作波长处显著的手性光学效应.

      手性超表面的CD响应在实际加工后难以调控. 为此, 在结构表面引入单层石墨烯, 实现了对CD的动态调制, 如图6(a)所示. 石墨烯的表面电导率(σg)可通过Kubo公式得到, 包括带内和带间的贡献, 表述如下[51,52]:

      其中, e为电子电荷, $\hbar$为约化普朗克常数, kB为玻尔兹曼常数, T为温度, EF是费米能级, τ是弛豫时间, ω = 2πf是角频率. 此处设置T = 300 K, τ = 10–13 s. 通过结合电导率σg和石墨烯厚度, 可计算出石墨烯的介电常数:

      其中, ε0为真空介电常数, tg为单层石墨烯片厚度. 根据(9)式可知, 石墨烯的介电常数由其表面电导率σg决定. 在不同费米能级EF下, σg发生变化, 进而引起εg的实部与虚部相应改变. 图6(b), (c)分别展示了不同费米能级下, 介电常数实部与虚部随波长变化的曲线. 其中, 虚部反映石墨烯的光吸收特性, 在EF ≈ $\hbar $ω/2附近, imag(εg)的变化趋势出现明显拐点. 当EF>$\hbar $ω/2时, imag(εg)急剧下降并趋近于零, 表明石墨烯由有损材料转变为无损材料. 这是由于在近红外波段, 入射光子能量较高, 当费米能级EF低于光子能量的一半时, 带间跃迁贡献显著, 导致其具有较强的光吸收, 表现为较大的介电常数虚部, 石墨烯表现为有损材料, 会抑制超表面的本征手性共振; 当EF超过光子能量一半时, 泡利阻塞效应抑制了带间跃迁, 带内跃迁占主导, imag(εg)急剧下降, 石墨烯趋近于无损状态, 使得结构中原本被吸收所抑制的手性模式得以恢复, 从而在调控EF时实现了CD信号的显著变化. 因此, 通过调节EF可主动调控石墨烯的表面电导率, 进而实现对手性光学响应的动态调制, 为近红外波段可调谐手性器件的实现提供了极大便利.

      如图6(d)所示, 在保持Δd = 92.25 nm, δ = 25 nm时, 随着费米能级EF的增大, RCP吸收共振增强, 手性q-BIC的带宽逐渐变窄. 当EF接近1 eV时, 由于泡利封锁效应增强, RCP的吸收率达到极值. 同时, 在费米能级变化的过程中, LCP吸收率始终保持极低水平, 几乎可忽略不计. 图6(e)展示了不同费米能级下的CD光谱, 随着费米能级的升高, CD光谱的峰值逐渐增大并在EF = 0.7 eV后趋于稳定. 将费米能级从0.4 eV增至0.7 eV, 可实现CD谱线峰值从–0.230到–0.952的调控, 最大调控深度达到0.722. 此外, 随着费米能级的增大, CD谱线发生蓝移现象, 该现象源于石墨烯介电常数实部随费米能级变化的色散特性, 如图6(b)所示. 本研究在保持结构参数不变的条件下, 通过调控石墨烯费米能级, 实现了手性响应的动态调制.

    • 本研究提出了一种基于双凹槽锗纳米柱的可调谐手性响应超表面, 通过打破面内C2对称性, 将理想的BIC转变为q-BIC. 该设计在近红外波段实现了超高Q因子(Q = 6121.14)与强圆二色性(CD = –0.94)响应的协同优化, 证明了在单一结构中兼顾高Q因子与显著手性的可行性. 在调控方面, 通过调节C2对称性破缺程度实现了手性响应在动量空间的调控, 并创新性地利用双凹槽中心间距Δd实现了CD的符号反转, 为手性态的可控设计提供了清晰理论机制. 进一步, 通过引入石墨烯层, 阐明了通过调节费米能级改变材料损耗与色散, 从而动态调制q-BIC共振特性与CD响应的物理路径, 实现了CD值在–0.230—–0.952范围内的大幅度调谐(调谐深度0.722), 拓展了在可调谐手性器件中的应用潜力. 本研究为发展高性能、可重构手性光子器件提供了新思路. 所提出的结构兼具高Q因子带来的强场局域能力, 以及通过几何参数与电学调控实现的手性灵活操控, 在手性光电子器件领域具有重要的应用潜力.

    参考文献 (52)

目录

/

返回文章
返回