Research Article

## Broadband Negative Refraction of Highly Squeezed Hyperbolic Polaritons in 2D Materials

## Abstract

Negative refraction of highly squeezed polaritons is a fundamental building block for nanophotonics, since it can enable many unique applications, such as deep-subwavelength imaging. However, the phenomenon of all-angle negative refraction of highly squeezed polaritons, such as graphene plasmons with their wavelength squeezed by a factor over 100 compared to free-space photons, was reported to work only within a narrow bandwidth (<1 THz). Demonstrating this phenomenon within a broad frequency range remains a challenge that is highly sought after due to its importance for the manipulation of light at the extreme nanoscale. Here we show the broadband all-angle negative refraction of highly squeezed hyperbolic polaritons in 2D materials in the infrared regime, by utilizing the naturally hyperbolic 2D materials or the hyperbolic metasurfaces based on nanostructured 2D materials (e.g., graphene). The working bandwidth can vary from several tens of THz to over a hundred of THz by tuning the chemical potential of 2D materials.

## 1. Introduction

Realizing negative refraction of highly squeezed polaritons, especially that supported by two-dimensional (2D) materials [1–6], such as graphene plasmon polaritons, is an important step toward the active manipulation of light at the extreme nanoscale [7–9] and can promise many photonic and optoelectronic applications [10–17]. In 2017, the phenomenon of all-angle negative refraction between highly squeezed* isotropic* graphene plasmons and hexagonal boron nitride’s (BN) phonon polaritons, with their in-plane polaritonic wavelengths squeezed by a factor over 100, is theoretically shown possible in the graphene-BN heterostructures [18]. By following [18], henceforth we use the parameter of squeezing factor to define the ratio between the photon wavelength in free space and the in-plane polaritonic wavelength. The extreme spatial confinement of these highly squeezed isotropic polaritons, however, also limits our ability to tailor their dispersion relations [18, 19]; in other words, their effective negative refractive index exists only within a certain frequency range. Consequently, the phenomenon of all-angle negative refraction of these highly squeezed isotropic polaritons is restricted to only a narrow frequency range (near 23 THz) with a bandwidth of less than 1 THz [18]. The narrow bandwidth unavoidably hinders the demonstration of this exotic phenomenon in experiments and thus the potential applications. It remains a challenge to realize the all-angle negative refraction of highly squeezed polaritons within a broad frequency range.

Here we predict the phenomenon of broadband all-angle negative refraction of highly squeezed* hyperbolic* polaritons, with their squeezing factors over 100, in the infrared regime. The highly squeezed hyperbolic polaritons can be supported by nanostructured 2D materials such as graphene nanoribbon arrays [20, 21] or BN nanoribbon arrays [22, 23], by 2D materials supported by nanostructured substrates [13, 24], or by naturally existing hyperbolic 2D materials such as black phosphorous [2, 25]. In addition, we note that the hyperbolic polaritons can also be supported by metal based metasurfaces at visible [26, 27] and even microwave [28] regimes. The peculiar isofrequency contour of hyperbolic polaritons [2, 20, 29] gives us an extra freedom to tailor their in-plane propagation direction and thus the flexibility to realize the all-angle negative refraction of highly squeezed polaritons in a wide frequency range. As an example, for the all-angle negative refraction of hyperbolic polaritons shown below, the working bandwidth can vary from several tens of THz to over a hundred of THz by simply tuning the chemical potential of 2D materials. Our work thus indicates that 2D materials are a versatile platform for the design of advanced nanometasurfaces and nanoimaging elements, which are of fundamental importance for the active control of light at the nanoscale.

We note that there are other types of negative refraction studied in the platform of 2D materials, including the negative refraction of electrons [30, 31], negative refraction of light propagating through a monolayer graphene [32], and plasmonic [33, 34] and optical [35, 36] negative refraction in 3D bulk materials (i.e., 2D material-based periodic structures). However, these electromagnetic refractions [32–36] occur out of the plane of 2D materials and not in plane, thus without taking advantage of the high spatial confinement of polaritons that is critical in this work. In addition, the all-angle negative refraction of graphene plasmons was studied in the hybrid graphene-photonic crystal structures [37], but the working bandwidth is still far less than 1 THz for a given chemical potential of graphene, along with the squeezing factor less than 20.

## 2. Results and Discussion

To highlight the underlying physics, we begin with the dispersion of hyperbolic polaritons supported by a uniaxial metasurface. The metasurface, such as that in the left region of Figure 1(a), can be modelled by an effective anisotropic surface conductivity of . Contrary to the isotropic metasurface (i.e., ) which supports the propagation of either transverse-magnetic (TM) [38, 39] or transverse-electric (TE) polaritons [40–43], the uniaxial metasurface supports the hybrid TM-TE polaritons [20, 24, 44–46]. From the electromagnetic theory [47], the dispersion relation for hybrid polaritons [44–46] is derived asIn the above, is the relative permittivity of region above/below the metasurface; and are the out-of-plane and in-plane wavevectors of the hybrid polariton, respectively; is the angle between and ; is the angular frequency; , , and are the permittivity, permeability, and light speed in free space, respectively. It is worth noting that the left side of (1) is only relevant to the TM field components of hybrid polaritons, while the denominator of the right side of (1) only corresponds to the TE field components. To some extent, the numerator of the right side of (1) denotes the coupling strength between TM and TE field components.

Importantly, for the highly squeezed hybrid polaritons (i.e., ), the right side of (1) is shown to be much smaller than 1, due to the difference of the spatial confinement between pure TE and TM polaritons; see Supplementary Materials. This way, the dispersion relation of hybrid polaritons in (1) can be compactly reduced to [20, 24]In other words, the highly squeezed hybrid TM-TE polaritons are dominant by the TM field components. The isofrequency contour governed by (2) is hyperbolic, when .

Figure 1(b) shows the hyperbolic isofrequency contour of highly squeezed polaritons supported by a graphene metasurface at 15 THz. As a conceptual demonstration, here we use the nanostructured graphene, i.e., graphene nanoribbon arrays, to create the hyperbolic metasurface; see the schematic structure in Figure 1(a). According to [20], when the pitch of nanoribbons is much smaller than the polaritonic wavelength , i.e., , the effective medium theory can be applied to describe the graphene metasurface. Then the effective anisotropic surface conductivity of graphene metasurface can be described by and , where is the width of nanoribbon, is an equivalent conductivity associated with the near-field coupling between adjacent nanoribbons, and is the surface conductivity of monolayer graphene modelled by the Kubo formula [20, 48]. Here the nanostructured graphene has a chemical potential of 0.1 eV, a conservative electron mobility of 10000 cm^{2}V^{−1}s^{−1} [38, 39], a pitch of nm, and a width of nm. The region below/above the graphene metasurface is the dielectric substrate (e.g., SiO_{2}) and air, respectively.

For the emergence of refraction phenomenon, the graphene metasurface in the right region in Figure 1(a) shall be different from the left region. One simple way is to rotate the graphene metasurface in the right region with a certain angle with respect to the left region, such as 90^{0} shown in Figure 1(a). This way, the graphene metasurface in the right region can be characterized by a surface conductivity of , where and . By applying the conservation law for wavevectors parallel to the boundary between left and right regions (i.e., parallel to ), we can determine the eigenmodes for the incident polariton in the left region and the transmitted polariton in the right region; see Figure 1(b). Through judicious design (e.g., by rotating the right region 90° with respect to the left region), we can create the case that the directions of the component of group velocity parallel to the interface (i.e., the y-component) are opposite for the eigenmodes of incident and transmitted polaritons, as shown in Figure 1(b). Since the group velocity (instead of the phase velocity) determines the direction of power flow, the directions of the component of power flow parallel to the interface are also opposite for the incident and transmitted polaritons. This enables the negative refraction of hyperbolic polaritons at the interface between left and right regions. The above underlying mechanisms for negative refraction of polaritons are the same as [18, 49]. In addition, we note that the hyperbolic polaritons in the left region can only propagate within a certain range of directions. Figure 1(b) indicates that, for these hyperbolic polaritons (i.e., incident from the left region with arbitrary angles), negative refraction can always happen at the boundary. Here we denote this phenomenon as the all-angle negative refraction [18, 49, 50] of hyperbolic polaritons (i.e., graphene plasmons).

Figure 1(c) numerically demonstrates the all-angle negative refraction of hyperbolic polaritons at 15 THz, by using the finite-element method (COMSOL Multiphysics). Hyperbolic polaritons are excited by a dipole source in the left region and propagate directionally towards the boundary. At the boundary, the polaritonic beams are negatively refracted. Moreover, these beams converge to form an image in the right region, which computationally validates the all-angle negative refraction. We note that the full width at half maximum (FWHM) of the image is only 0.035 *μ*m [Fig. S2], which is less than 1/100 of the working wavelength (i.e., 20 *μ*m) in free space and can enable deep-subwavelength imaging. Here the squeezing factor [see Figure 1(b)] indicates that when compared with the wavelength in free space, the polaritonic wavelength is squeezed at least by a factor over 100. It is worth emphasizing that the all-angle negative refraction of highly squeezed hyperbolic polaritons with directional propagation, as an important advantage over that of isotropic polaritons which propagate omnidirectionally [18, 46, 47], might facilitate the design of novel compact guidance. We note that the imaging mechanism here is not a perfect image recovery. This is because the reflection at the boundary is unavoidable due to the impedance mismatch (caused by the modal mismatch) between the two hyperbolic metasurfaces, and the propagation loss of hyperbolic polaritons will degrade the imaging quality; see Supplementary Materials. To optimize the quality of formed image in the right region, one shall consider both the reflection and the propagation loss of polaritons [18, 46, 47].

Figure 2 shows the effective anisotropic conductivity of graphene metasurface, which can help to infer the working frequency range of the all-angle negative refraction of hyperbolic polaritons. As shown in Figure 2(a), the frequency range that has and thus supports the hyperbolic polaritons spans from 0 to 47 THz, when the chemical potential of graphene is 0.1 eV. It shall be noted that the phenomenon of all-angle negative refraction of hyperbolic polaritons can happen at arbitrary frequency within this frequency range, as long as the effective medium theory for graphene metasurface is valid (i.e., when ). To guarantee the validity of the effective medium theory for metasurfaces based on 2D materials, a small value of pitch , which although might increase the complexity in structural fabrication [22, 23], can be adopted. This way, the working frequency range for the all-angle negative refraction of highly squeezed polaritons revealed here is no longer limited by the frequency range supporting negative-index polaritons (such as [18]). Consequently, the working frequency range can be broadband and actively controllable via tuning the chemical potential of 2D materials. In addition, the relative bandwidth for negative refraction of hyperbolic polaritons can be up to 2, if we define the relative bandwidth as , where and () are the maximum and minimum frequencies supporting the negative refraction of hyperbolic polaritons.

Figure 2(b) shows the working bandwidth of all-angle negative refraction in graphene metasurfaces. The bandwidth of graphene metasurface having changes from 40 THz to over 100 THz, by increasing the chemical potential from 0.1 eV to 0.5 eV. Therefore, in principle, the bandwidth of all-angle negative refraction of hyperbolic polaritons can vary from several tens of THz to even over a hundred THz, by simply increasing the chemical potential of graphene. Interestingly, Figure 2(b) shows that the loss can increase the bandwidth when is smaller than 0.18 eV; see analysis in Supplementary Materials.

To numerically validate the all-angle negative refraction of hyperbolic polaritons in a broad bandwidth, Figure 3 demonstrates this phenomenon at other frequencies, i.e., at 10 THz in Figure 3(a) and at 20 THz in Figure 3(b). The squeezing factors for hyperbolic polaritons at these two frequencies are both over 100; see Supplementary Materials. Therefore, by using hyperbolic metasurfaces based on 2D materials, we can extend the working bandwidth of all-angle negative refraction of highly squeezed polaritons to at least several tens of THz, which is favored for practical applications.

To further extend the bandwidth of all-angle negative refraction of highly squeezed polaritons, one may adopt the naturally anisotropic 2D materials to support tunable hyperbolic polaritons, such as those described in [2] including black phosphorous. From [2], these 2D materials can be directly characterized by an anisotropic surface conductivity, i.e.,where is the concentration of electrons, is the electron’s effective mass along the direction, is the relaxation time, is the frequency of the onset of interband transitions for the component of conductivity, accounts for the strength of interband component, and is a step function. Equation (3) circumvents the requirement (i.e., discussed in the above) for the validity of effective medium theory for metasurfaces based on nanostructured 2D materials. For these anisotropic 2D materials, the value of can be flexibly tunable via electrostatic gating, just like graphene, and can be different from . These, along with the broad class of anisotropic materials, give us the flexibility to realize the all-angle negative refraction of highly squeezed polaritons in a broad frequency range.

As a concrete example, Figure 4 shows the all-angle negative refraction of hyperbolic polaritons supported by the naturally anisotropic 2D materials. By following [2], here we set , , , , , , and , where and is the free-electron mass. The frequency range, which has and thus supports the hyperbolic polaritons, spans from 24 to 143 THz in Figure 4(a), with a bandwidth of 119 THz. This indicates that the negative refraction of hyperbolic polaritons can happen within the above frequency range; see the negative refraction of hyperbolic polaritons at 60 THz in Figure 4(b), for example.

In conclusion, we have revealed a viable way to realize the all-angle negative refraction of highly squeezed polariton in a broadband infrared regime, by utilizing hyperbolic metasurfaces based on 2D materials or naturally anisotropic 2D materials. Due to the combined advantages of highly directional propagation, active tunability, low loss, and ultrahigh confinement provided by hyperbolic polaritons in 2D materials, the broad class of 2D materials can provide a versatile platform for the manipulation of light-matter interaction at the extreme nanoscale and for the design of highly compact nanodevices and circuits.

## 3. Materials and Methods

The finite element simulation is implemented via the frequency domain simulation in the commercial software of COMSOL Multiphysics. To enable the high calculation accuracy, the 2D material is modeled as a surface, where an anisotropic surface conductivity is used to fulfill the conditions for discontinuities in the electromagnetic fields. The meshing resolution in the plane of 2D material is 5 nm. For Figures 1(c) and 3, a z-polarized dipole source and for Figure 4 a y-polarized electric dipole is placed in the left region at 5 nm above the 2D material; the phenomenon of all-angle negative refraction is pronounced, independent of the vertical position of the source. The fields in Figures 1(c), 3, and 4 are obtained at the plane with 2 nm above the plane of 2D material. For the clarity of conceptual demonstration, the value of is artificially set to be equal to in Figures 1(c), 3, and 4. For the cases considering the realistic material loss, the phenomenon of all-angle negative refraction is shown in Fig. S3. The phenomenon of all-angle negative refraction is also independent of the relative permittivity of the dielectric substrate , as shown in Fig. S6; in the main text, we set .

## Data Availability

All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.

## Conflicts of Interest

The authors declare that they have no conflicts of interest.

## Authors’ Contributions

Xiao Lin conceived the research. Jing Jiang performed the main calculation. Xiao Lin and Baile Zhang contributed insight and discussion on the results. Jing Jiang, Xiao Lin, and Baile Zhang wrote the paper. Xiao Lin and Baile Zhang supervised the project. Jing Jiang and Xiao Lin contributed equally to this work.

## Acknowledgments

The authors thank Y. Yang, I. Kaminer and M. Soljačić for useful discussions. This work was sponsored by the Singapore Ministry of Education (Grant nos. MOE2015-T2-1-070, MOE2016-T3-1-006, and Tier 1 RG174/16 (S)).

## Supplementary Materials

Section S1: Dispersion of hybrid polaritons supported by anisotropic metasurfaces. Section S2: All-angle negative refraction of hyperbolic graphene plasmons. Section S3: Loss influence on the bandwidth having . Fig. S1. Real part of effective surface conductivity of graphene metasurface. Fig. S2. Full width at half maximum of the image for the point source in Figure 1(c). Fig. S3. All-angle negative refraction of hyperbolic polaritons when the real material loss is considered. Fig. S4. Isofrequency contours of hyperbolic graphene plasmons at 10 THz, 15 THz and 20 THz. Fig. S5. Isofrequency contours of hyperbolic graphene plasmons at 1 THz and 40 THz. Fig. S6. Substrate influence on the all-angle negative refraction of hyperbolic graphene plasmons at 15 THz. Fig. S7. All-angle negative refraction of hyperbolic polaritons in nanostructures of patterned 2D materials with a pitch of nm. Fig. S8. Loss influence on the bandwidth having * (Supplementary Materials)*

## References and notes

- D. N. Basov, M. M. Fogler, and F. J. García de Abajo, “Polaritons in van der Waals materials,”
*Science*, vol. 354, no. 6309, 2016. View at Scopus · View at Google Scholar - A. Nemilentsau, T. Low, and G. Hanson, “Anisotropic 2D Materials for Tunable Hyperbolic Plasmonics,”
*Physical Review Letters*, vol. 116, no. 6, 2016. View at Publisher · View at Google Scholar - T. Low, A. Chaves, J. D. Caldwell et al., “Polaritons in layered two-dimensional materials,”
*Nature Materials*, vol. 16, no. 2, pp. 182–194, 2017. View at Publisher · View at Scopus · View at Google Scholar - X. Lin, I. Kaminer, X. Shi et al., “Splashing transients of 2D plasmons launched by swift electrons,”
*Science Advances*, vol. 3, no. 1, 2017. View at Scopus · View at Google Scholar - F. H. L. Koppens, D. E. Chang, and F. J. García de Abajo, “Graphene plasmonics: a platform for strong light–matter interactions,”
*Nano Letters*, vol. 11, no. 8, pp. 3370–3377, 2011. View at Publisher · View at Scopus · View at Google Scholar - L. J. Wong, I. Kaminer, O. Ilic, J. D. Joannopoulos, and M. Soljačić, “Towards graphene plasmon-based free-electron infrared to X-ray sources,”
*Nature Photonics*, vol. 10, no. 1, pp. 46–52, 2016. View at Publisher · View at Scopus · View at Google Scholar - G. X. Ni, A. S. McLeod, Z. Sun et al., “Fundamental limits to graphene plasmonics,”
*Nature*, vol. 557, no. 7706, pp. 530–533, 2018. View at Publisher · View at Google Scholar - D. Alcaraz Iranzo, S. Nanot, E. J. Dias et al., “Probing the ultimate plasmon confinement limits with a van der Waals heterostructure,”
*Science*, vol. 360, no. 6386, pp. 291–295, 2018. View at Publisher · View at Google Scholar - X. Shi, X. Lin, I. Kaminer et al., “Superlight inverse Doppler effect,”
*Nature Physics*, vol. 14, no. 10, pp. 1001–1005, 2018. View at Publisher · View at Google Scholar - E. Yoxall, M. Schnell, A. Y. Nikitin et al., “Direct observation of ultraslow hyperbolic polariton propagation with negative phase velocity,”
*Nature Photonics*, vol. 9, no. 10, pp. 674–678, 2015. View at Publisher · View at Scopus · View at Google Scholar - S. Dai, Q. Ma, T. Andersen et al., “Subdiffractional focusing and guiding of polaritonic rays in a natural hyperbolic material,”
*Nature Communications*, vol. 6, 2015. View at Scopus · View at Google Scholar - S. Chakraborty, O. P. Marshall, T. G. Folland, Y.-J. Kim, A. N. Grigorenko, and K. S. Novoselov, “Applied optics: Gain modulation by graphene plasmons in aperiodic lattice lasers,”
*Science*, vol. 351, no. 6270, pp. 246–248, 2016. View at Publisher · View at Scopus · View at Google Scholar - A. Vakil and N. Engheta, “Transformation optics using graphene,”
*Science*, vol. 332, no. 6035, pp. 1291–1294, 2011. View at Publisher · View at Scopus · View at Google Scholar - S. Dai, Z. Fei, Q. Ma et al., “Tunable phonon polaritons in atomically thin van der Waals crystals of boron nitride,”
*Science*, vol. 343, no. 6175, pp. 1125–1129, 2014. View at Publisher · View at Scopus · View at Google Scholar - P. Li, M. Lewin, A. V. Kretinin et al., “Hyperbolic phonon-polaritons in boron nitride for near-field optical imaging and focusing,”
*Nature Communications*, vol. 6, 2015. View at Scopus · View at Google Scholar - G. X. Ni, L. Wang, M. D. Goldflam et al., “Ultrafast optical switching of infrared plasmon polaritons in high-mobility graphene,”
*Nature Photonics*, vol. 10, no. 4, pp. 244–247, 2016. View at Publisher · View at Scopus · View at Google Scholar - C. Qian, X. Lin, Y. Yang et al., “Multifrequency Superscattering from Subwavelength Hyperbolic Structures,”
*ACS Photonics*, vol. 5, no. 4, pp. 1506–1511, 2018. View at Publisher · View at Google Scholar - X. Lin, Y. Yang, N. Rivera et al., “All-angle negative refraction of highly squeezed plasmon and phonon polaritons in graphene–boron nitride heterostructures,”
*Proceedings of the National Acadamy of Sciences of the United States of America*, vol. 114, pp. 6717–6721, 2017. View at Publisher · View at Google Scholar - Y. Jiang, X. Lin, T. Low, B. Zhang, and H. Chen, “Group-velocity-controlled and gate-tunable directional excitation of polaritons in graphene-boron nitride heterostructures,”
*Laser & Photonics Reviews*, vol. 12, no. 5, Article ID 1870024, 2018. View at Publisher · View at Google Scholar - J. S. Gomez-Diaz, M. Tymchenko, and A. Alù, “Hyperbolic Plasmons and Topological Transitions Over Uniaxial Metasurfaces,”
*Physical Review Letters*, vol. 114, no. 23, 2015. View at Publisher · View at Google Scholar - H. Yan, T. Low, W. Zhu et al., “Damping pathways of mid-infrared plasmons in graphene nanostructures,”
*Nature Photonics*, vol. 7, no. 5, pp. 394–399, 2013. View at Publisher · View at Scopus · View at Google Scholar - P. Li, I. Dolado, F. J. Alfaro-Mozaz et al., “Infrared hyperbolic metasurface based on nanostructured van der Waals materials,”
*Science*, vol. 359, no. 6378, pp. 892–896, 2018. View at Publisher · View at Scopus · View at Google Scholar - J. J. López, A. Ambrosio, S. Dai et al., “Photothermal Effect: Large Photothermal Effect in Sub-40 nm h-BN Nanostructures Patterned Via High-Resolution Ion Beam (Small 22/2018),”
*Small*, vol. 14, no. 22, p. 1870101, 2018. View at Publisher · View at Google Scholar - E. Forati, G. W. Hanson, A. B. Yakovlev, and A. Alù, “Planar hyperlens based on a modulated graphene monolayer,”
*Physical Review B: Condensed Matter and Materials Physics*, vol. 89, no. 8, Article ID 081410, 2014. View at Publisher · View at Scopus · View at Google Scholar - X. Huang, Y. Cai, X. Feng et al., “Black Phosphorus Carbide as a Tunable Anisotropic Plasmonic Metasurface,”
*ACS Photonics*, vol. 5, no. 8, pp. 3116–3123, 2018. View at Publisher · View at Google Scholar - A. A. High, R. C. Devlin, A. Dibos et al., “Visible-frequency hyperbolic metasurface,”
*Nature*, vol. 522, no. 7555, pp. 192–196, 2015. View at Publisher · View at Scopus · View at Google Scholar - J. S. T. Smalley, F. Vallini, S. A. Montoya et al., “Luminescent hyperbolic metasurfaces,”
*Nature Communications*, vol. 8, 2017. View at Scopus · View at Google Scholar - Y. Yang, L. Jing, L. Shen et al., “Hyperbolic spoof plasmonic metasurfaces,”
*NPG Asia Materials*, vol. 9, no. 8, p. e428, 2017. View at Publisher · View at Google Scholar - A. Poddubny, I. Iorsh, P. Belov, and Y. Kivshar, “Hyperbolic metamaterials,”
*Nature Photonics*, vol. 7, pp. 948–957, 2013. View at Google Scholar - G.-H. Lee, G.-H. Park, and H.-J. Lee, “Observation of negative refraction of Dirac fermions in graphene,”
*Nature Physics*, vol. 11, no. 11, pp. 925–929, 2015. View at Publisher · View at Scopus · View at Google Scholar - S. Chen, Z. Han, M. M. Elahi et al., “Electron optics with p-n junctions in ballistic graphene,”
*Science*, vol. 353, no. 6307, pp. 1522–1525, 2016. View at Publisher · View at Scopus · View at Google Scholar - H. Harutyunyan, R. Beams, and L. Novotny, “Controllable optical negative refraction and phase conjugation in graphite thin films,”
*Nature Physics*, vol. 9, no. 7, pp. 423–425, 2013. View at Publisher · View at Scopus · View at Google Scholar - B. Wang, X. Zhang, F. J. García-Vidal, X. Yuan, and J. Teng, “Strong coupling of surface plasmon polaritons in monolayer graphene sheet arrays,”
*Physical Review Letters*, vol. 109, no. 7, Article ID 073901, 2012. View at Publisher · View at Scopus · View at Google Scholar - H. Huang, B. Wang, H. Long, K. Wang, and P. Lu, “Plasmon-negative refraction at the heterointerface of graphene sheet arrays,”
*Optics Expresss*, vol. 39, no. 20, pp. 5957–5960, 2014. View at Publisher · View at Scopus · View at Google Scholar - K. V. Sreekanth, A. De Luca, and G. Strangi, “Negative refraction in graphene-based hyperbolic metamaterials,”
*Applied Physics Letters*, vol. 103, no. 2, Article ID 023107, 2013. View at Publisher · View at Scopus · View at Google Scholar - A. A. Sayem, M. M. Rahman, M. R. C. Mahdy, I. Jahangir, and M. S. Rahman, “Negative Refraction with Superior Transmission in Graphene-Hexagonal Boron Nitride (hBN) Multilayer Hyper Crystal,”
*Scientific Reports*, vol. 6, 2016. View at Scopus · View at Google Scholar - S. Zhong, Y. Lu, C. Li, H. Xu, F. Shi, and Y. Chen, “Tunable plasmon lensing in graphene-based structure exhibiting negative refraction,”
*Scientific Reports*, vol. 7, no. 1, 2017. View at Publisher · View at Google Scholar - M. Jablan, H. Buljan, and M. Soljačić, “Plasmonics in graphene at infrared frequencies,”
*Physical Review B: Condensed Matter and Materials Physics*, vol. 80, no. 24, Article ID 245435, 2009. View at Publisher · View at Scopus · View at Google Scholar - X. Lin, N. Rivera, J. J. López, I. Kaminer, H. Chen, and M. Soljačić, “Tailoring the energy distribution and loss of 2D plasmons,”
*New Journal of Physics*, vol. 18, no. 10, p. 105007, 2016. View at Publisher · View at Google Scholar - S. A. Mikhailov and K. Ziegler, “New electromagnetic mode in graphene,”
*Physical Review Letters*, vol. 99, no. 1, Article ID 016803, 2007. View at Publisher · View at Scopus · View at Google Scholar - F. Gargiulo and O. V. Yazyev, “Structural and electronic transformation in low-angle twisted bilayer graphene,”
*2D Materials*, vol. 5, no. 1, 2018. View at Scopus · View at Google Scholar - X. Lin, Y. Shen, I. Kaminer, H. Chen, and M. Soljačić, “Transverse-electric Brewster effect enabled by nonmagnetic two-dimensional materials,”
*Physical Review A: Atomic, Molecular and Optical Physics*, vol. 94, no. 2, 2016. View at Publisher · View at Google Scholar - M. Y. Musa, M. Renuka, X. Lin et al., “Confined transverse electric phonon polaritons in hexagonal boron nitrides,”
*2D Materials*, vol. 5, no. 1, p. 015018, 2018. View at Publisher · View at Google Scholar - J. O. Island, M. Buscema, M. Barawi et al., “Ultrahigh photoresponse of few-layer TiS3 nanoribbon transistors,”
*Advanced Optical Materials*, vol. 2, no. 7, pp. 641–645, 2014. View at Publisher · View at Scopus · View at Google Scholar - J. Dai and X. C. Zeng, “Titanium Trisulfide Monolayer: Theoretical Prediction of a New Direct-Gap Semiconductor with High and Anisotropic Carrier Mobility,”
*Angewandte Chemie*, vol. 127, no. 26, pp. 7682–7686, 2015. View at Publisher · View at Google Scholar - J. Kang, H. Sahin, H. D. Ozaydin, R. T. Senger, and F. M. Peeters, “TiS
^{3}nanoribbons: Width-independent band gap and strain-tunable electronic properties,”*Physical Review B: Condensed Matter and Materials Physics*, vol. 92, no. 7, Article ID 075413, 2015. View at Publisher · View at Google Scholar - J. A. Kong,
*Electromagnetic Wave Theory*, EMW Publishing, Cambridge, Mass, USA, 2008. View at MathSciNet - G. W. Hanson, “Dyadic green's functions for an anisotropic, non-local model of biased graphene,”
*IEEE Transactions on Antennas and Propagation*, vol. 56, no. 3, pp. 747–757, 2008. View at Publisher · View at Scopus · View at Google Scholar - H. Shin and S. Fan, “All-angle negative refraction for surface plasmon waves using a metal-dielectric-metal structure,”
*Physical Review Letters*, vol. 96, Article ID 073907, 2006. View at Google Scholar - H. J. Lezec, J. A. Dionne, and H. A. Atwater, “Negative refraction at visible frequencies,”
*Science*, vol. 316, no. 5823, pp. 430–432, 2007. View at Publisher · View at Scopus · View at Google Scholar