In recent years, there have been many interests in studying electromagnetic superscattering because of its special applications[
All of the above-mentioned superscattering phenomena exhibit a single-frequency regime. Multifrequency superscattering may enable improved sensitivity and accuracy for extensive applications ranging from sensing, bioimaging, and optical tagging to spectroscopy[
In this work, we consider a two-dimensional subwavelength cylindrical structure. The structure consists of six layers of dielectric material, each of which is coated with a graphene shell, as shown in Fig. 1(a). Such a six-layer subwavelength structure is the main focus of this analysis, because six layers is sufficient for multifrequency superscattering and superscattering shaping. The conductivity of graphene can be expressed as the Kubo formula: [
Sign up for Chinese Optics Letters TOC Get the latest issue of Advanced Photonics delivered right to you！Sign up now
Figure 1.(a) Illustration of a subwavelength multilayered cylindrical structure consisting of six layers of dielectric materials, each of which is coated with a graphene shell. Three-dimensional and two-dimensional views of the structure are shown in the left and right inset, respectively. Graphene shells are represented by black dotted lines. Here, the plane wave is TM-polarized and incident from the left side of the cylinder structure. (b) The scattering spectra (normalized scattering cross section) for the same structure in (a) under the ideal lossless assumption. Both the total scattering spectrum and the partial spectra of the first five scattering terms are shown. The inset in (b) shows the details at the resonance peak circled by a dotted blue line.
Here, is the surface current of graphene satisfying Ohm’s law. To derive the scattering coefficient , dynamical matrix and are defined as
Then, from the recurrence relation of the coefficients of the adjacent layers, the relation between the coefficients of the innermost layer and the outermost layer can be described by the following formula:
Here, NSCS is normalized by , where is the wavelength of the incident wave. According to the electromagnetic multipole analysis of a cylindrical structure[
3. Results and Discussions
To highlight the underlying physical properties, we first consider a lossless subwavelength cylindrical structure by setting , where is the relaxation time of charge carriers of graphene. The scattering spectra (NSCS) as functions of the frequency for the structure in Fig. 1(a) are shown in Fig. 1(b). Both the total scattering spectrum and the partial spectra from channels with , , , , and are plotted. The geometrical parameters are and . Relative permittivities of dielectrics are , , , , , and . For all graphene shells, the chemical potential is at temperature , and the thickness can be negligible. In the system, resonances from channels with , , , and are mainly considered because the scattering cross sections from other channels are too small to be ignored. At 3.79715 THz, the total scattering cross section mostly comes from resonant channels with and , and the scattering cross sections from these four channels almost all reach the signal-channel limit, which makes the total reach four times the signal-channel limit. The total scattering cross section at 9.77464 THz exceeds four times the signal-channel limit, which originates from a combination of the resonance peaks from channels with and . Scattering cross sections of these peaks almost all reach signal-channel limit, and the contributions from channels with are negligible. Similarly, the resonance peaks at 12.65616 THz from channels with and are all nearly equal to the signal-channel limit, and they are overlapped well, which results in the total scattering cross section reaching four times the signal-channel limit. The resonance peaks with and are overlapped at 14.74908 THz, and these peaks are almost all equal to the signal-channel limit. Obviously, the contributions from channels with cannot be neglected, so the total scattering cross section is greater than four times the signal-channel limit at 14.74908 THz. This implies that multifrequency superscattering can be realized by overlapping electric and magnetic resonances with different orders.
The underlying mechanism of the resonant overlapping of multipoles in the subwavelength cylinder structure can be qualitatively understood from the dispersion of plasmon polaritons supported by a corresponding planar structure in Fig. 2. Starting from resonance modes supported in the corresponding planar structure, the resonances of the subwavelength cylinder structure can be understood using the whispering gallery condition based on the Bohr model[
Figure 2.Dispersion of plasmon polaritons versus |RTM| (in logarithmic scale) under the ideal lossless assumption is shown. RTM is the reflection coefficient of the TM-polarized plane wave for the corresponding planar structure. The parameters of the corresponding planar structure are the same as those in Fig.
Figure 3 intuitively shows the performance of the superscattering pattern shaping in the form of far and near fields for the cylinder structure in Fig. 1(a) at the four superscattering resonant frequencies. Far-field angular scattering amplitude can be analyzed by the derived formula[
Figure 3.Scattering angular distributions for the subwavelength cylinder structure in Fig.
Theoretically, various scattering pattern shaping can also be realized by resonantly overlapping different electromagnetic multipoles because these multipoles show different phase symmetries and parities[
Figure 3(a) shows the superscattering pattern shaping by resonantly overlapping ED and EQ at the resonant frequency of 3.79715 THz. Far-field scattering angular distribution is shown by a blue solid curve, and the scattering near field is shown in the background pattern, which is in accord with the blue solid curve. As shown in Fig. 3(a), the scattering is enhanced in the forward direction but suppressed in the backward direction, because the scatterings of ED and EQ exhibit odd and even parities, respectively. The two multipoles are in phase in the forward direction and out of phase in the backward direction. In fact, in the forward direction, the two resonantly overlapped multipoles are always in phase according to the optical theorem[
Figure 3(c) illustrates the superscattering pattern shaping with resonantly overlapped EQ and MO by showing the scattering-field distribution at the frequency of 9.77464 THz. The far-field scattering angular distribution is plotted by a blue solid curve, which matches well with the scattering near field shown in the background pattern. Like Fig. 3(a), the scattering is enhanced in the forward direction but suppressed in the backward direction due to opposite parities of EQ and MO (even and odd, respectively). The two overlapped multipoles are in phase and out of phase, respectively, in the forward and backward directions. It is observed that the scattering is also suppressed from 30° to 45° and 315° to 330° in the forward half-scattering circle, which can be understood analytically by when EQ and MO are only overlapping. The red crosses plot the ideal scattering pattern with the two modes purely overlapped. The slight discrepancy between blue solid curve and red crosses originates from the fact that the contribution of the ED is negligible. For the near field, the total field at 9.77464 THz is shown in Fig. 3(d). A remarkable “shadow” can also be seen in the forward direction, and the “shadow” is much larger than the diameter of the cylindrical structure, which agrees well with the scattering field in Fig. 3(c).
Figure 3(e) illustrates the superscattering pattern shaping with resonantly overlapped ED and MO by showing the scattering-field distribution at the frequency of 12.65616 THz. The scattering near and far fields correspond to the background pattern and the blue solid curve, respectively, and they are well matched. Clearly, the scattering is enhanced symmetrically both in the forward and backward directions due to the same odd parity of the scatterings of ED and MO. Additionally, the interferences of overlapped ED and MO can also suppress the scattering at , , , , , and due to . When only considering the overlapped ED and MO, the ideal scattering angular distribution plotted by red crosses matches well with the blue solid curve because other contributions except for ED and MO are too small to be neglected. For the near field, the total magnetic field at the resonant frequency of 12.65616 THz is shown in Fig. 3(f), and the cylinder structure leaves a large “shadow” larger than the structure diameter, not only in the forward direction but also in backward direction, which is consistent with the scattering field in Fig. 3(e).
Figure 3(g) shows the superscattering pattern shaping by resonantly overlapping EQ and EO at the frequency of 14.74908 THz in the form of the scattering far field plotted by the blue solid curve and the scattering near field shown in the bottom pattern. The scattering amplitudes are symmetrically enhanced both in forward and backward directions, which is caused by the same even parity of the EQ and EO. It can be also seen that the scattering is suppressed from 22.5° to 45°, 67.5° to 112.5°, 135° to 157.5°, 202.5° to 225°, 247.5° to 292.5°, and 315° to 337.5° according to . The ideal scattering pattern with purely overlapped ED and EO plotted by red crosses does not coincide well with the blue solid curve because the contribution of the ED cannot be neglected at the frequency of 14.74908 THz. For the near field, at 14.74908 THz, the remarkable “shadow”, which is much larger than the structure diameter, can be seen in both the forward and backward directions in Fig. 3(h), and the total magnetic field matches well with the scattering field shown in Fig. 3(g).
To facilitate the realization of the experiment and practical application, the effects of the material losses and the structural changes are analyzed in detail in Fig. 4. First, considering the optical losses of graphene, as shown in Fig. 4(a), the superscattering with resonantly overlapped multipoles still exists at corresponding resonant frequencies when changing the relaxation time , and the NSCSs of corresponding resonance peaks are slightly reduced due to the addition of graphene’s losses. Then, to illustrate the high tolerance to structural design, the total scattering spectra are plotted in Fig. 4(b) by varying thickness of the innermost layer of the dielectric and the corresponding scattering spectra of , , and are shown by the green, blue, and red solid curves, respectively. Except for , the other parameters are the same as those in the structure of Fig. 1(a). As the thickness increases gradually, though the positions of the resonant frequencies will move to the lower frequencies, superscattering with resonantly overlapped multipoles can still be obtained, and the NSCSs at the corresponding resonant frequencies are nearly unchanged. Compared with metal materials, graphene has the advantage of adjustability. Next, the tunability of graphene shells is illustrated by changing the chemical potential . The corresponding scattering spectra of , , and are shown by the green, blue, and red solid curves in Fig. 4(c), respectively. When the chemical potential increases, the positions of resonant frequencies will move towards the higher frequencies, and the superscattering phenomenon still exists with the NSCSs of corresponding resonance peaks almost unchanged. From Fig. 4, when considering the material losses and changing the structural parameters, superscattering with resonantly overlapped multipoles can be still obtained, which indicates that superscattering for a graphene-based subwavelength structure has a good tolerance to material losses and variations of structural parameters.
Figure 4.(a) Influence of graphene’s losses on multifrequency superscattering is shown through the total scattering spectra for different lossy cases. Except τ, other parameters are the same as those in the structure of Fig.
In conclusion, we propose that a graphene-based subwavelength structure can be a promising and versatile platform for the demonstration of simultaneous multifrequency superscattering and superscattering shaping with different engineered scattering patterns. Based on Bohr’s model and dispersion engineering, it is further demonstrated that such efficient superscattering pattern shaping is induced by the interferences of multipoles resonantly overlapped. In contrast with multilayer plasmonic structures, the performance of the proposed structure can be easily refined for the desired frequency by adjusting the chemical potential and loss of the graphene layer. Owing to multiple polaritonic dispersion lines in a graphene-based subwavelength structure, the phenomena of multifrequency superscattering are found to have a high tolerance to some structural variations and realistic graphene loss. In addition, more exotic scattering phenomena are expected in such a platform, such as dynamically tunable multifrequency cloaking and other multifrequency scattering-based devices.
 A. Kinkhabwala, Z. Yu, S. Fan, Y. Avlasevich, K. Müllen, W. E. Moerner. Large single-molecule fluorescence enhancements produced by a bowtie nanoantenna. Nat. Photon., 3, 654(2009).
 S. Arslanagic, R. W. Ziolkowski. Highly subwavelength, superdirective cylindrical nanoantenna. Phys. Rev. Lett., 120, 237401(2018).
 C. Loo, A. Lowery, N. Halas, J. West, R. Drezek. Immunotargeted nanoshells for integrated cancer imaging and therapy. Nano Lett., 5, 709(2005).
 J. Sheng, J. Xie, J. Liu. Multiple super-resolution imaging in the second band of gradient lattice spacing photonic crystal flat lens. Chin. Opt. Lett., 18, 120501(2020).
 J. Zhao, X. Y. Zhang, C. R. Yonzon, A. J. Haes, R. P. Van Duyne. Localized surface plasmon resonance biosensors. Nanomedicine, 1, 219(2006).
 J. Wang, X. Wang, M. Zeng. Broadband transverse displacement sensing of silicon hollow nanodisk under focused radial polarization illumination in the near-infrared region. Chin. Opt. Lett., 18, 063602(2020).
 A. Schliesser, N. Picqúe, T. W. Hänsch. Mid-infrared frequency combs. Nat. Photon., 6, 440(2012).
 M. A. Green, S. Pillai. Harnessing plasmonics for solar cells. Nat. Photon., 6, 130(2012).
 W. Xu, L. Xie, Y. Ying. Mechanisms and applications of terahertz metamaterial sensing: a review. Nanoscale, 9, 13864(2017).
 L. Novotny, B. Hecht. Principles of Nano-optics(2012).
 Y. Zhai, Y. Ma, S. N. David, D. Zhao, R. Lou, G. Tan, R. Yang, X. Yin. Scalable-manufactured randomized glass-polymer hybrid metamaterial for daytime radiative cooling. Science, 355, 1062(2017).
 T. Yang, H. Chen, X. Luo, H. Ma. Superscatterer: enhancement of scattering with complementary media. Opt. Express, 16, 18545(2008).
 W. H. Wee, J. B. Pendry. Shrinking optical devices. New J. Phys., 11, 073033(2009).
 Z. Ruan, S. Fan. Superscattering of light from subwavelength nanostructures. Phys. Rev. Lett., 105, 013901(2010).
 R. J. Li, X. Lin, S. S. Lin, X. Liu, H. S. Chen. Tunable deep-subwavelength superscattering using graphene monolayers. Opt. Lett., 40, 1651(2015).
 R. J. Li, X. Lin, S. S. Lin, X. Liu, H. S. Chen. Atomically thin spherical shell-shaped superscatterers based on Bohr model. Nanotechnology, 26, 505201(2015).
 R. Li, B. Zheng, X. Lin, R. Hao, S. Lin, W. Yin, E. Li, H. Chen. Design of ultracompact graphene-based superscatterers. IEEE J. Sel. Top. Quantum Electron., 23, 4600208(2017).
 W. Liu. Ultra-directional super-scattering of homogenous spherical particles with radial anisotropy. Opt. Express, 23, 14734(2015).
 W. Liu, B. Lei, J. Shi, H. Hu. Unidirectional superscattering by multilayered cavities of effective radial anisotropy. Sci. Rep., 6, 34775(2016).
 W. Liu. Superscattering pattern shaping for radially anisotropic nanowires. Phys. Rev. A, 96, 023854(2017).
 W. Wan, W. Zheng, Y. Chen, Z. Liu. From Fano-like interference to superscattering with a single metallic nanodisk. Nanoscale, 6, 9093(2014).
 Y. Huang, L. Gao. Superscattering of light from core-shell nonlocal plasmonic nanoparticles. J. Phys. Chem. C, 118, 30170(2014).
 A. Mirzaei, A. Miroshnichenko, I. Shadrivov, Y. Kivshar. Superscattering of light optimized by a genetic algorithm. Appl. Phys. Lett., 105, 011109(2014).
 C. Wang, C. Qian, H. Hu, L. Shen, Z. Wang, H. Wang, Z. Xu, B. Zhang, H. Chen, X. Lin. Superscattering of light in refractive-index near-zero environments. Prog. Electromagn. Res., 168, 15(2020).
 M. Zhou, L. Shi, J. Zi, Z. Yu. Extraordinarily large optical cross section for localized single nanoresonator. Phys. Rev. Lett., 115, 023903(2015).
 P. Del’Haye, T. Herr, E. Gavartin, M. Gorodetsky, R. Holzwarth, T. J. Kippenberg. Octave spanning tunable frequency comb from a microresonator. Phys. Rev. Lett., 107, 063901(2011).
 F. C. Cruz, D. L. Maser, T. Johnson, G. Ycas, A. Klose, F. R. Giorgetta, I. Coddington, S. A. Diddams. Mid-infrared optical frequency combs based on difference frequency generation for molecular spectroscopy. Opt. Express, 23, 26814(2015).
 C. Qian, X. Lin, Y. Yang, F. Gao, Y. Shen, J. Lopez, I. Kaminer, B. Zhang, E. Li, M. Soljačić, H. Chen. Multifrequency superscattering from subwavelength hyperbolic structures. ACS Photon., 5, 1506(2018).
 H. W. Wu, Y. Fang, J. Q. Quan, Y. Z. Han, Y. Q. Yin, Y. Li, Z. Q. Sheng. Multifrequency superscattering with high Q factors from a deep-subwavelength spoof plasmonic structure. Phys. Rev. B, 100, 235443(2019).
 C. Qian, X. Lin, Y. Yang, X. Xiong, H. Wang, E. Li, I. Kaminer, B. Zhang, H. Chen. Experimental observation of superscattering. Phys. Rev. Lett., 122, 063901(2019).
 V. I. Shcherbinin, V. I. Fesenko, T. I. Tkachova, V. R. Tuz. Superscattering from subwavelength corrugated cylinders. Phys. Rev. Appl., 13, 024081(2020).
 S. H. Raad, C. J. Zapata-Rodríguez, Z. Atlasbaf. Multi-frequency super-scattering from sub-wavelength graphene-coated nanotubes. J. Opt. Soc. Am. B, 36, 2292(2019).
 R. Kumar, K. Kajilawa. Superscattering from cylindrical hyperbolic metamaterials in the visible region. Opt. Express, 28, 1507(2020).
 R. Kumar, K. Kajikawa. Comparison of cylinder- and planar-effective medium approximations on calculation of scattering properties of cylindrical hyperbolic metamaterials. J. Opt. Soc. Am. B, 36, 559(2019).
 J. Chai, P. Hu, L. Ge, H. Xiang, D. Han. Tunable terahertz cloaking and lasing by the optically pumped graphene wrapped on a dielectric cylinder. J. Phys. Commun., 3, 035016(2019).
 W. Liu, A. E. Miroshnichenko, Y. S. Kivshar. Q- factor enhancement in all-dielectric anisotropic nanoresonators. Phys. Rev. B, 94, 195436(2016).
 K. Koshelev, A. Bogdanov, Y. Kivshar. Meta-optics and bound states in the continuum. Sci. Bull., 64, 836(2019).
 T. Zhan, X. Shi, Y. Dai, X. Liu, J. Zi. Transfer matrix method for optics in graphene layers. J. Phys. Condens. Matter., 25, 215301(2013).
 W. Liu, Y. S. Kivshar. Generalized Kerker effects in nanophotonics and meta-optics. Opt. Express, 26, 13085(2018).
 W. Liu. Generalized magnetic mirrors. Phys. Rev. Lett., 119, 123902(2017).
 H. C. van de Hulst. Light Scattering by Small Particles(1981).
Set citation alerts for the article
Please enter your email address