Abstract
Coherent perfect absorber (CPA), a resonator with critical losses that can perfectly absorb all incident light, has been observed at various frequency regimes (from microwave to visible light). Besides the functional frequency, the bandwidth is also an important parameter in characterizing the performance of CPA. Here, we explore the bandwidth of CPA in a kind of weakly-coupled-resonance metasurfaces with , where is the near-field coupling between the radiative and non-radiative resonant modes, and is the scattering loss rate of the radiative resonant mode. Based on the coupled mode theory, we analytically derive the upper and lower bounds of the bandwidth, and show that they are determined by the dissipation loss rates of the composed modes. To narrow the bandwidth, it is better to increase the radiative loss rate when designing a weakly coupled resonator. We also show that CPA is associated with a robust phase singularity with a winding number of ± 1. The conclusions are numerically verified in a designed resonant metasurface and could perform as a guideline for designing CPA in various resonant systems.
© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
1. Introduction
Absorption is an unavoidable phenomenon, besides the scattering, in various systems of wave. In general, absorption loss is believed to be the main factor impeding the utilizations of metamaterials and metasurfaces supporting localized resonances [1]. However, recently a class of synthetic systems with tailored absorption has recently generated huge interest in the optical and physical societies because they can reveal the virtue of absorption [2–8]. Among them the most representative one is the coherent perfect absorber (CPAs) [3–5]. A CPA refers to a resonator with critical losses that can perfectly absorb the entire incident light, and is believed to be the time-reversed counterpart of a laser operating at the threshold. Mathematically, the existence of CPA coincides with a zero eigenvalue of the scattering matrix S sitting at the real frequency axis. Physically, perfect absorption is due to a destructive interference in a critical-tuned two-port resonator [9–13]. Since the theoretical prediction by Chong et. al [3], most efforts were focused on the possible realization of CPAs in various systems such as photonic crystals and acoustic/optical metasurfaces [9–23].
The realization of CPA requires the presence of critical losses. As a critical underdamped resonator, the quality factor (Q factor) defined by the frequency-to-bandwidth ratio is a very important dimensionless parameter [24] of CPA. The bandwidth δf usually refers to the full-width half-maximum (FWHM) of the absorption peak. It is desirable to have a large Q factor to enhance the interaction between light and the resonator for potential applications such as in sensors and nonlinear optics. There are many approaches to increase the Q factor in photonic systems. For example, we can make use of avoided resonance crossing in a pair of coupled modes, where one of the eigen-modes has a reduced bandwidth [25–27]. On the contrary, a small Q factor, corresponding to a larger bandwidth, is also necessary to achieve an effective absorption within a broader frequency range. The detailed values and mechanisms of bounds in Q factor would help us to solve this tradeoff and provide us with a guideline in improving the performance of CPA toward a pre-defined application. Although there have seen a lot of theoretical and experimental verifications of CPAs at various frequency regimes, the bandwidth bond of CPA is rarely mentioned.
Here, we theoretically explore the bandwidth of CPAs within the scope of resonant metasurfaces by using the coupled mode theory (CMT), which is a very concise model to characterize resonant phenomena in various physical systems [28,29]. We show that in a typical resonant metasurface with weakly coupled bright and dark modes, there exist upper and lower bandwidth bounds that are determined by the dissipation loss rate. If the dissipation loss is very weak, it is better to increase the radiation loss rate of the bright resonant mode to narrow the bandwidth of CPA. In addition, the position of CPA in such a hybridized system is a robust phase singularity with a winding number of . The phase singularity implies an anomalous phase delay. To verify our theoretical explorations, we propose a metallic metasurface design that bears a sharp transition from positive to negative phase delay at the resonance of CPA. The numerical result is in good agreement with CMT. Our theoretical exploration can provide a guideline for designing CPAs in various resonant systems.
2. Theoretical conclusions from CMT
To derive the bounds on the bandwidth of CPA in a resonant metasurface, we adopt the simple CMT that can grasp main physical behaviors of various coupled resonators not limited to metasurfaces [30–35]. Without loss of generality, let us consider a metasurface exhibiting a mirror symmetry along the propagation direction. We further assume that the resonance of the metasurface is composed of a radiative resonant mode with a scattering loss rate of , and a non-radiative mode with a zero scattering. The complex amplitude of the coupled resonance is represented by, and usually corresponds to the energy in each resonant mode, where the subscripts and indicate the radiative and non-radiative modes, respectively. Such a kind of system with coupled bright and dark modes is a representative test-bed for electromagnetically induced transparency (EIT), Fano resonance and Autler-Townes effect [35]. To simplify our discussion on CPA, the resonant frequencies of the two modes are assumed to be equal, i.e. . This assumption has been used to investigate slow light phenomena in the classical analogue of EIT effect [30,31,33]. Different from former investigations, here we focus on CPA associated with a zero eigenvalue of the so-called scattering matrix. Considering a near-field coupling between these two resonant modes, and indicating the radiative coupling between the radiative mode and the input (output) light by, we get two coupled equations [28–31]:
The response matrix, directive coupling matrix between the input and the output fields, and the radiative coupling matrix () between the output (input) field and the resonance in the metasurface arerespectively. Here γd is the intrinsic dissipation rate by the absorption in the resonator, which can be made extremely small, e.g. by working in the microwave regime or by using proper lossless media. The relationship between the input and output fields is then given by , where the scattering matrix is formulated byMathematically, the position of CPA corresponds to a zero eigenvalue of . Due to the mirror symmetry along the propagation direction, CPA only occurs for a symmetric eigenvector of when the damping and coupling rates satisfy
To be physically observable, the solution of should be exactly at the real frequency axis by setting constraint conditions on, , and. In the strong coupling scenario of , by solving Eq. (4) it is readily to show that should be satisfied, and the frequencies of CPA are given by . CPA in such a strong coupling scenario has been explored to imitate the phase transition in parity-time-symmetric Hamiltonian systems [36]. The strict requirement on and constrains the realization of CPA in a limited frequency range.Here we are interested in the weak-coupling scenario of , and would like to discuss the exact bandwidth and explore the bounds on the existing CPA because it is desirable to maximize\minimize the operational bandwidth δf by optimizing the resonator design in many potential applications. In this weak-coupling scenario, the null of Eq. (4) requires that both the real and imaginary parts are zero. It is then clear that the frequency of CPA is always because . Starting from the extreme case of where the resonance of the metasurface is only dominated by the radiative dipolar mode, from Eq. (4) we can see that CPA at requires either a balanced scattering condition of [37] or a zero intrinsic dissipation of γd = 0. The latter case is physically impossible and should be neglected. Now, FWHM of the coherent absorption can be shown to be , which is only determined by the dissipation loss rate γd. The Q factor is then .
When is not zero, the bright and dark modes would couple with each other. Now CPA occurs at a real frequency of under the satisfaction of a critical coupling of
According to Eqs. (3) and (5), we can find the frequencies where the output intensity of output field b is half of the incident field a = (1, 1)T. The difference between them is defined as FWHM of CPA. Under the simple assumption of , FWHM of the coherent absorption is found to beWe can see that δf is determined not only by the dissipation loss rate , but also by the scattering one . For a selected system working at a desired frequency, the dissipation loss is generally fixed by the geometry of the structure and the constituent media forming it. Unlike the isolated radiative resonance of , we could engineer FWHM in weak-coupling scenario by optimizing the radiation loss rate with the help of resonator design, i.e. to modify the brightness of the radiative mode. When we select the scattering loss rate as the varied parameter, according to Eq. (6), there exists an upper bound in FWHMwhenWhen , FWHM is approximated given by. Together with the scenario of κ = 0 as discussed above, we can see that CPA only occurs at , and the lower bound of FWHM according to Eq. (6) is
Since the resonant frequency of CPA is fixed at, we can also utilize the Q factor to characterize the bandwidth. Figure 1 shows the effective Q factor of f0/δf versus and is confined by the upper and lower bounds predicted by Eqs. (7) and (8). For potential applications we can engineer the value of to obtain the desired functionality within the bounded region.Note that in the two-dimensional (2D) parameters space formed by and the solution of CPA satisfying Eq. (4) is a phase singularity with a winding number of [see the example shown in Fig. 1(b)]. The sign of the winding number is determined by the sign of the coupling strength. The existence of this robust phase singularity, which is immune to γd, γs and κ, implies an anomalous phase delay around CPA [38–40]. The phase delay for the symmetric eigenvalue of the scattering matrix at the resonant frequency is given by
The second term of the denominator is the one that gives the critical condition of Eq. (5), and consequently, τd indeed changes its sign when f crosses the frequency f0 of CPA.3. Metasurface design and numerical verification
In order to put the theoretical conclusion of CMT made in the former section into realistic systems, a resonant metallic metasurface as schematically shown in Fig. 2(a) is proposed. The dissipation loss rate of the metal, such as copper used here, is very small in the microwave region. Now it is difficult to satisfy the condition of in the metallic metasurface to approach the upper bound of Q factor. However, instead of searching for a very small radiative loss, we can switch to the scenario of a large value in our metasurface design according to the analysis present in Section 2. A large value could provide a FWHM of that can still approach the lower bound of Eq. (8). For example, if the value of equals 1000, δf is about.
In our structure design shown in Fig. 2(a), the metal of the resonant metasurface is copper with a deep subwavelength thickness of 0.5 mm and an electric conductivity ofS/m. The normally incident plane wave is linearly polarized with electric field E oscillating in the direction. Each unit cell is a rectangular lattice with a dimension of mm2. Width of each metallic strip is mm. The non-radiative mode is supported by the middle two metallic strips oriented along the direction, where the length of each strip is mm with a gap of mm. The radiative mode is generally from the outer two metallic strips oriented along the direction, where the length of each strip is mm with a gap size of mm. The tunable parameter here is the transverse displacement of the center metallic strips. When is varied, not only the coupling strength but also some other parameters such as the resonant frequency and the scattering loss rate are modified.
To analyze the performance of CPA, we firstly numerically calculate the complex transmission and reflection coefficients by using COMSOL Multiphysics. The frequency regime is from 18.5 to 22.5 GHz, and the transverse displacements varies from 0 to 0.48 mm with a step of 0.02 mm. Then, we retrieve a set of CMT parameters (red spheres in Fig. 2) by fitting the numerical reflection results with the analytical formulae of Eqs. (1)-(3). Finally, the discrete CMT parameters are fitted by using a polynomial fitting method to reveal simple relation between them and . The fitted polynomial curves (in the unit of GHz) versus tx (in the unit of mm) are displayed in Fig. 2 as solid lines, as: , , , and .
The fitting shows that an increased transverse displacement always leads to a blue (red) shift in the resonance of the radiative (non-radiative) mode. The scattering loss rate γs shows neglectable deviation from 5.24 GHz, and the dissipation loss rate γd is a constant of 0.005GHz, three orders smaller than γs.
Based on the above fitted function for each parameter and Eqs. (1)-(3), the amplitude of reflection is in good agreement with that from numerical simulation, see Fig. 3. Furthermore, we can obtain the coherent absorption spectra,, which is shown in Fig. 4. Around mm we get GHz, GHz, GHz, andGHz, respectively. Based on these discrete fitted parameters, the operational frequency of CPA is approximately given by GHz. These fitted parameters verify the assumptions ofand the weak coupling criterion of made in the former section. To clearly indicate the position of CPA, we plot the numerical and fitted coherent absorption spectra at mm in Fig. 5(a). The value of the coherent absorption from fitting (numerical simulation) at 20.55GHz (20.554 GHz) is 0.9978 (0.9977), and can be assumed to be unity. Based on the coherent absorption spectra, we can obtain the fitted (numerical) FWHM of CPA. The effective Q factor of CPA, defined by , can reach 1027 (1027), which is high enough for applications in sensors and nonlinear optics. The obtained Q factor is almost the same as that from our theoretical prediction, i.e., and approaches the upper bound of .
Besides the large Q factor close to the upper bound, we also verify that the eigenvalue of the scattering matrix at the CPA in such hybridized system is a phase singularity, see Eq. (4). Around such a phase singularity, the phase delay shows a transition from positive to negative when we change the coupling strength, see Eq. (9). Examples about the phase delays below ( mm) and above (mm) the critical coupling condition of CPA are shown in Fig. 5(b). The peaks of phase delay in these two conditions are opposite with each other. Note that when we vary the transverse displacement to tune the coupling strength κ, the resonant frequency of CPA also changes. The frequency of CPA is very close to the resonant frequency of the non-radiative mode, so the resonant peak\dip of the phase delay below\above the critical coupling also exhibits a weak red shift.
4. Conclusion
In summary, we derive the bandwidth of CPA in a resonant metasurface by using CMT, and show that the upper and lower bounds are determined by the dissipation loss rate of the resonances inside. Weakly-coupled resonators with a giant radiative loss rate can be utilized to get a narrow bandwidth approaching the lower bound of FWHM. We numerically verify our analysis in a designed metallic metasurface structure in the microwave frequency range, where the dissipation loss rate is very low. The numerical results are in good agreement with CMT. Our investigation can be transferred to higher frequency regimes such as infrared and visible light.
Funding
National Natural Science Foundation of China (NSFC) (11574162, 11674244, 11874228); Program for Innovative Research in University of Tianjin (TD13-5077).
References
1. W. Cai and V. Shalaev, Optical metamaterials: fundamentals and applications (Springer, 2009).
2. C. M. Bender and S. Boettcher, “Real spectra in non-Hermitian Hamiltonians having PT symmetry,” Phys. Rev. Lett. 80(24), 5243–5246 (1998). [CrossRef]
3. Y. D. Chong, L. Ge, H. Cao, and A. D. Stone, “Coherent perfect absorbers: time-reversed lasers,” Phys. Rev. Lett. 105(5), 053901 (2010). [CrossRef] [PubMed]
4. S. Longhi, “PT-symmetric laser absorber,” Phys. Rev. A 82, 031801 (2010).
5. W. Wan, Y. Chong, L. Ge, H. Noh, A. D. Stone, and H. Cao, “Time-reversed lasing and interferometric control of absorption,” Science 331(6019), 889–892 (2011). [CrossRef] [PubMed]
6. Y. D. Chong, L. Ge, and A. D. Stone, “PT-symmetry breaking and laser-absorber modes in optical scattering systems,” Phys. Rev. Lett. 106(9), 093902 (2011). [CrossRef] [PubMed]
7. L. Feng, R. El-Ganainy, and L. Ge, “Non-Hermitian photonics based on parity–time symmetry,” Nat. Photonics 11(12), 752–762 (2017). [CrossRef]
8. W. Wang, L. Q. Wang, R. D. Xue, H. L. Chen, R. P. Guo, Y. Liu, and J. Chen, “Unidirectional excitation of radiative-loss-free surface plasmon polaritons in parity-time symmetric systems,” Phys. Rev. Lett. 119(7), 077401 (2017). [CrossRef] [PubMed]
9. S. Dutta-Gupta, R. Deshmukh, A. Venu Gopal, O. J. F. Martin, and S. Dutta Gupta, “Coherent perfect absorption mediated anomalous reflection and refraction,” Opt. Lett. 37(21), 4452–4454 (2012). [CrossRef] [PubMed]
10. M. Pu, Q. Feng, M. Wang, C. Hu, C. Huang, X. Ma, Z. Zhao, C. Wang, and X. Luo, “Ultrathin broadband nearly perfect absorber with symmetrical coherent illumination,” Opt. Express 20(3), 2246–2254 (2012). [CrossRef] [PubMed]
11. S. Feng and K. Halterman, “Coherent perfect absorption in epsilon-near-zero metamaterials,” Phys. Rev. B Condens. Matter Mater. Phys. 86(16), 165103 (2012). [CrossRef]
12. J. W. Yoon, G. M. Koh, S. H. Song, and R. Magnusson, “Measurement and modeling of a complete optical absorption and scattering by coherent surface plasmon-polariton excitation using a silver thin-film grating,” Phys. Rev. Lett. 109(25), 257402 (2012). [CrossRef] [PubMed]
13. J. Zhang, K. F. MacDonald, and N. I. Zheludev, “Controlling light-with-light without nonlinearity,” Light Sci. Appl. 1(7), e18 (2012). [CrossRef]
14. S. Zanotto, F. P. Mezzapesa, F. Bianco, G. Biasiol, L. Baldacci, M. Vitiello, L. Sorba, R. Colombelli, and A. Tredicucci, “Perfect energy-feeding into strongly coupled systems and interferometric control of polariton absorption,” Nat. Phys. 10(11), 830–834 (2014). [CrossRef]
15. Y. Sun, W. Tan, H. Q. Li, J. Li, and H. Chen, “Experimental demonstration of a coherent perfect absorber with PT phase transition,” Phys. Rev. Lett. 112(14), 143903 (2014). [CrossRef] [PubMed]
16. M. Kang, Y. D. Chong, H.-T. Wang, W. Zhu, and M. Premaratne, “Critical route for coherent perfect absorption in a Fano resonance plasmonic system,” Appl. Phys. Lett. 105(13), 131103 (2014). [CrossRef]
17. P. Wei, C. Croenne, S. T. Chu, and J. Li, “Symmetrical and anti-symmetrical coherent perfect absorption for acoustic waves,” Appl. Phys. Lett. 104(12), 121902 (2014). [CrossRef]
18. J. Z. Song, P. Bai, Z. H. Hang, and Y. Lai, “Acoustic coherent perfect absorbers,” New J. Phys. 16(3), 033026 (2014). [CrossRef]
19. N. Arju, T. Ma, A. Khanikaev, D. Purtseladze, and G. Shvets, “Optical realization of double-continuum Fano interference and coherent control in plasmonic metasurfaces,” Phys. Rev. Lett. 114(23), 237403 (2015). [CrossRef] [PubMed]
20. G. Pirruccio, M. Ramezani, S. R. Rodriguez, and J. G. Rivas, “Coherent control of the optical absorption in a plasmonic lattice coupled to a luminescent layer,” Phys. Rev. Lett. 116(10), 103002 (2016). [CrossRef] [PubMed]
21. D. G. Baranov, A. E. Krasnok, T. Shegai, A. Alù, and Y. D. Chong, “Coherent perfect absorbers: linear control of light with light,” Nat. Rev. Mater. 2(12), 17064 (2017). [CrossRef]
22. S. Longhi, “Coherent virtual absorption for discretized light,” Opt. Lett. 43(9), 2122–2125 (2018). [CrossRef] [PubMed]
23. Y. Li and C. Argyropoulos, “Tunable nonlinear coherent perfect absorption with epsilon-near-zero plasmonic waveguides,” Opt. Lett. 43(8), 1806–1809 (2018). [CrossRef] [PubMed]
24. M. Tooley, Electronic circuits: fundamentals and applications, 3rd ed. (Elsevier, 2006).
25. J. Wiersig, “Formation of long-lived, scarlike modes near avoided resonance crossings in optical microcavities,” Phys. Rev. Lett. 97(25), 253901 (2006). [CrossRef] [PubMed]
26. Q. H. Song and H. Cao, “Improving optical confinement in nanostructures via external mode coupling,” Phys. Rev. Lett. 105(5), 053902 (2010). [CrossRef] [PubMed]
27. M. V. Rybin, K. L. Koshelev, Z. F. Sadrieva, K. B. Samusev, A. A. Bogdanov, M. F. Limonov, and Y. S. Kivshar, “High-Q supercavity modes in subwavelength dielectric resonators,” Phys. Rev. Lett. 119(24), 243901 (2017). [CrossRef] [PubMed]
28. H. A. Haus, Waves and fields in optoelectronics (Prentice-Hall, 1984).
29. W. Suh, Z. Wang, and S. Fan, “Temporal coupled mode theory and the presence of non-orthogonal modes in lossless multimode cavities,” IEEE J. Quantum Electron. 40(10), 1511–1518 (2004). [CrossRef]
30. S. Zhang, D. A. Genov, Y. Wang, M. Liu, and X. Zhang, “Plasmon-induced transparency in metamaterials,” Phys. Rev. Lett. 101(4), 047401 (2008). [CrossRef] [PubMed]
31. N. Liu, L. Langguth, T. Weiss, J. Kästel, M. Fleischhauer, T. Pfau, and H. Giessen, “Plasmonic analogue of electromagnetically induced transparency at the Drude damping limit,” Nat. Mater. 8(9), 758–762 (2009). [CrossRef] [PubMed]
32. A. Artar, A. A. Yanik, and H. Altug, “Multispectral plasmon induced transparency in coupled meta-atoms,” Nano Lett. 11(4), 1685–1689 (2011). [CrossRef] [PubMed]
33. P. Tassin, L. Zhang, R. Zhao, A. Jain, T. Koschny, and C. M. Soukoulis, “Electromagnetically induced transparency and absorption in metamaterials: the radiating two-oscillator model and its experimental confirmation,” Phys. Rev. Lett. 109(18), 187401 (2012). [CrossRef] [PubMed]
34. L. Verslegers, Z. Yu, Z. Ruan, P. B. Catrysse, and S. Fan, “From electromagnetically induced transparency to superscattering with a single structure: a coupled-mode theory for doubly resonant structures,” Phys. Rev. Lett. 108(8), 083902 (2012). [CrossRef] [PubMed]
35. M. F. Limonov, M. V. Rybin, A. N. Poddubny, and Y. S. Kivshar, “Fano resonances in photonics,” Nat. Photonics 11(9), 543–554 (2017). [CrossRef]
36. M. Kang, F. Liu, and J. Li, “Effective spontaneous PT-symmetry breaking in hybridized metamaterials,” Phys. Rev. A 87(5), 053824 (2013). [CrossRef]
37. M. Kang, F. Liu, T. F. Li, Q. H. Guo, J. Li, and J. Chen, “Polarization-independent coherent perfect absorption by a dipole-like metasurface,” Opt. Lett. 38(16), 3086–3088 (2013). [CrossRef] [PubMed]
38. D. R. Solli, C. F. McCormick, R. Y. Chiao, S. Popescu, and J. M. Hickmann, “Fast light, slow light, and phase singularities: a connection to generalized weak values,” Phys. Rev. Lett. 92(4), 043601 (2004). [CrossRef] [PubMed]
39. M. Kang, Y. Li, and H. T. Wang, “Fingerprints of topological defects in a metasurface,” Opt. Lett. 39(16), 4879–4882 (2014). [CrossRef] [PubMed]
40. M. Asano, K. Y. Bliokh, Y. P. Bliokh, A. G. Kofman, R. Ikuta, T. Yamamoto, Y. S. Kivshar, L. Yang, N. Imoto, S. K. Özdemir, and F. Nori, “Anomalous time delays and quantum weak measurements in optical micro-resonators,” Nat. Commun. 7(1), 13488 (2016). [CrossRef] [PubMed]