Abstract
Channeling surface plasmon–polaritons to control their propagation direction is of the utmost importance for future optoelectronic devices. Here, we develop an effective-index method to describe and characterize the properties of a 2D material’s channel plasmon–polaritons (CPPs) guided along a V-shaped channel. Focusing on the case of graphene, we derive a universal Schrödinger-like equation from which one can determine the dispersion relation of graphene CPPs and corresponding field distributions at any given frequency, since they depend on the geometry of the structure alone. The results are then compared against more rigorous theories, having obtained very good agreement. Our calculations show that CPPs in graphene and other 2D materials are attractive candidates to achieve deep subwavelength waveguiding of light, holding potential as active components for the next generation of tunable photonic devices.
© 2017 Optical Society of America
1. INTRODUCTION
Surface plasmon–polaritons (SPPs) [1] have been considered a viable tool to tailor light–matter interactions at the nanoscale owing to their ability to squeeze electromagnetic fields below the diffraction limit [2,3]. Such property makes SPP modes particularly attractive in the design of photonic components that rival their electronic counterparts in terms of miniaturization, while delivering larger bandwidths and operational speeds [4]. Among the vast number of applications in optoelectronics where SPPs may be put into use, one of the most prominent is subwavelength waveguiding of electromagnetic radiation. In this context, gap-SPPs (SPP modes sustained at dielectric gaps separating two metal surfaces) are appealing candidates due to their favorable balance between losses and field confinement [5–7]. The width of the dielectric gap can either remain constant or vary in the direction perpendicular to it. An example of the latter configuration is a V-groove carved into a metallic substrate. Within this geometry, the corresponding propagating SPP modes are generally referred to as channel plasmon–polaritons (CPPs) [8,9]. Contrarily to elementary (planar) dielectric–metal interfaces, where the field is confined only along one direction, CPP modes render a two-dimensional (2D) confinement of the electromagnetic field in the plane normal to its propagation (the latter being along the translational invariant direction). As of today, a plethora of fundamental explorations [10–12] and proof-of-concept experiments have been carried out demonstrating the usage of noble metal CPPs in plasmonic interferometers [13,14], waveguides [15–19], and ring-resonators [13], and for nanofocusing [20]. In the context of quantum plasmonics [21], CPPs have also been explored for the control of quantum emitters [22]. Indeed, such realizations unambiguously announce that the field is reaching maturity [9].
In parallel, the rise of graphene as a novel plasmonic medium emerged only about half a decade ago [23–26]. In spite of the nanophotonics of 2D materials [27] being in its early stages, we have been witnessing extraordinarily rapid progress in graphene plasmonics, fueled by the prospect of long-lived, gate-tunable graphene plasmons (GPs) that yield large field confinements in the terahertz and mid-infrared [24–26]. At the time of writing, the majority of studies on graphene plasmons have dealt with nanostructured graphene (ribbons, disks, etc.) [23,25,28–31]. Nonplanar graphene plasmons have recently gained interest [32–35], while the investigation of graphene channel plasmons still remains largely unexplored. Only very recently has research on the plasmonic properties of graphene-covered triangular wedges and grooves been conducted [36–38].
In this paper, we derive a universal description to characterize channel plasmons guided along the apex of a triangularly shaped sheet of a conductive 2D material [graphene, 2D electron gas, doped 2D semiconductor transition metal dichalcogenides (TMDCs), etc.]; see Fig. 1 with the 2D material exemplified by graphene. Our approach follows in part the effective-index framework previously applied to CPPs occurring in traditional metals [6,39,40]. The concept behind the effective-index method (EIM) is that, at each height (with respect to the bottom of the channel), one can identify a one-dimensional (1D) dielectric/graphene/dielectric/graphene/dielectric (DGDGD) waveguide in which the core dielectric has a width given by . The combination of such 1D waveguide configurations—i.e., DGDGD slabs stacked together—can then be used to construct the whole 2D channel waveguide akin to graphene CPPs; see Fig. 1. Here, we show that this conceptually intuitive picture can be used to determine the dispersion of CPPs and corresponding electric field distributions by solving a Schrödinger-like equation whose eigenvalues depend only on the angle of the channel. Finally, the results obtained using this unified description are then compared against the ones calculated using more rigorous theories [36], having obtained good agreement within the EIM regime.
2. METHODS
We consider a V-shaped geometry making an angle as illustrated in Fig. 1, where it is assumed henceforth that in order for the EIM to be valid. Within this limit, one can make use of the transverse magnetic (TM) and electric (TE) representations. Therefore, we shall look for TM-like solutions in which the main component of the electric field lies along the axis. Recall that the underlying idea behind the EIM applied to the description of CPPs in V-grooves is to model the inner region as two coupled 1D waveguides. This foundational principle does not change when considering a dielectric V-channel covered with graphene (or any other 2D material). However, the corresponding equations will be quite different, asserting the natural differences between the two materials (i.e., 3D metal and 2D graphene). Therefore, we write the field as [6,39,40]
where denotes the propagation constant of the 2D material channel plasmon, and the coupled waveguide equation [stemming from the Helmholtz equation for the electric field (1)] reads with , and where acts as a separation constant; see Supplement 1. We note that in writing the above equations, we have neglected the and derivatives, whose contribution is small in the regime where [40]. Notice that Eq. (2a) is essentially an equation for a 1D waveguide, parameterized by . The specific form of is, in general, non-trivial. Therefore, our strategy here consists in developing from the GP condition at a planar graphene double-interface with the same dielectric media—DGDGD—which, for the mode with the sought-after symmetry, is given by [24] with , and where is the gap-width as a function of the coordinate. Here, is the dynamical conductivity of graphene (or, more generally, any conducting 2D crystal). A closed-form expression for may be obtained by expanding the previous equation for small widths [41], which leads to [42] where . In the case of graphene, materializes the assumption that the conductivity is given by its Drude-like expression with negligible damping [24], with denoting the fine-structure constant. In possession of an explicit relation for , a solution for the coupled waveguide Eq. (2) may be fetched. The solution to the first equation is trivial, simply being in the inner region, where is a constant. The remaining differential Eq. (2b) is somewhat less trivial. However, we can make such an equation more affordable by defining the wavevector , denoting plasmons propagating in the planar 2D material [24], and identifying , while at the same time introducing the dimensionless variable . Performing these transformations, we arrive at a dimensionless eigenvalue equation [43] resembling a Schrödinger equation with an “energy potential” where the 2D CPP dispersion relation is given in terms of the dimensionless eigenvalue . Interestingly, notice that mimics an “effective dielectric function” that is entirely determined by the geometry—namely the angle —of the channel. We stress that the solution of Eq. (5) yields the eigenfunctions and propagation constants of the CPP modes, while holding the advantage of portraying a clearer picture of the physics of the problem, which is epitomized by the universal scaling behavior for the CPP wavevector depending entirely on the geometry of the V-channel.We bring to the reader’s attention that the above differential equation for is substantially more intricate in our present case with a 2D material than it is for conventional CPPs in metallic V-grooves, where a simple differential equation and corresponding analytic solution can be straightforwardly derived [40]. In what follows, we solve the differential Eq. (5) numerically using the shooting method [44].
3. RESULTS AND DISCUSSION
The eigenfunctions resulting from the numerical solution of Eq. (5) are depicted in Fig. 2, for 2D CPP modes in a V-shaped channel with . As portrayed in the figure, these behave like bound modes trapped in a potential well, even displaying features familiar to the ones encountered in typical quantum mechanics problems, such as the decrease in energy difference as the mode order increases and the number of nodes given by (where labels the mode order). Moreover, Fig. 2 clearly shows that the fundamental CPP mode is highly confined within the bottom of the groove (i.e., near the apex of the channel), while the successively higher-order modes tend to be concomitantly more spread along the vertical direction of the channel [i.e., the () direction].
Having computed the eigenvalues and field profiles along the coordinate, one may readily construct the 2D electric field distributions, , corresponding to the aforementioned CPP modes using Eq. (1). The outcome of such a procedure is shown in Fig. 3, for V-shaped graphene with different opening angles, and excited at different resonant frequencies. The top panels illustrate the modal distributions for the fundamental [Fig. 3(a)] and second-order [Fig. 3(b)] CPP modes, whose main features echo our above description of Fig. 2, namely the higher level of localization achieved by the fundamental mode when compared with its leading order ones. In addition, the influence of the groove angle is reflected in the two bottom panels of the figure, with smaller opening angles yielding correspondingly higher field confinements near the apex of the channel. This behavior is a consequence of the stronger “confining potential” for CPPs in sharper V-structures. Furthermore, let us stress that we expect the propagation constants and corresponding electric field profiles to remain accurate even for finite-sized channels, as long as the truncation of the V-shaped structure occurs at sufficiently large distances away from the apex (i.e., at coordinates where the field is essentially zero). In this regard, the computed field patterns may be used to estimate the required channel depth/height a device based on 2D CPPs should have.
A characteristic of particular physical importance, both from a fundamental perspective and also from a device-engineering viewpoint, is the dispersion relation exhibited by the 2D CPPs. The spectrum of graphene CPPs guided along triangular channels stems from the solution of Eq. (5), which allows the retrieval of the propagation constant in terms of its planar graphene plasmon counterpart. Thus, with a single calculation one may construct the entire plasmon dispersion in the system of interest, for any frequency or any other parameter, since the solution depends uniquely on the angle of the structure. In Fig. 4 we show the obtained dispersion relations for graphene CPPs in representative triangular channels. Figures 4(a) and 4(b) exhibit the spectrum for the fundamental CPP in two dielectric environments with different relative permittivities, keeping the angle constant. As an eye-guide, we have also plotted the dispersion of graphene plasmons in flat graphene, along with the dispersion of light propagating freely inside the dielectric. It is apparent from the figure that the dispersion curve akin to graphene CPPs lies to the right of their corresponding flat GPs. This is indicative of the higher amount of field localization attained by the former when compared with the latter. The degree of subwavelength localization of the electromagnetic field is even more dramatic for smaller angles, as a comparison between Figs. 4(b) and 4(c) plainly shows. Additionally, in order to gauge the level of fidelity of our model we have also superimposed the dispersion of the graphene CPP modes obtained using a more rigorous technique [36] (dashed lines in Fig. 4). Performing a comparison of results, very good agreement can be observed regarding the dispersion of the fundamental CPP mode in Figs. 4(a) and 4(b), and both the fundamental and second-order modes in Fig. 4(c). This clearly illustrates that, despite its inherent simplicity, the EIM framework can still provide consistent results, thus being a very valuable tool in judging the main characteristics of a 2D material’s channel plasmons. Unfortunately, the EIM does not capture the higher-order modes in the case of (not shown); the justification for this lies in the breakage of the small-argument expansion performed in the cotangent figuring in Eq. (3), as the field distributions arising from these (higher-order) modes extend to regions relatively far away from the apex [cf. Fig. 3(b)]. Nevertheless, note that for smaller angles (e.g., ) the second-order mode is already well described by the EIM owing to the smaller gap. A more all-encompassing perspective can be laid out by taking advantage of the universal scaling of the graphene (or any 2D material) CPP spectrum with the angle , given by the eigenvalues . Figure 5 shows the scaling factor (since ) as a function of the opening angle of the system. The plot demonstrates quite remarkably the ability of our model to correctly determine the behavior of 2D channel plasmons. In the case of the fundamental mode, the agreement is maintained even beyond the small angle () regime. Naturally, and as already anticipated above, the description of the higher-order modes is limited to smaller angles. Figure 5 plays the role of an epilogue: it clearly displays the universal scaling law for the 2D CPP’s propagation constant as a function of the angle—valid independently of the frequency, 2D conductivity, or dielectric constant—while at the same time outlining the regime of validity of the model. Moreover, we have found that the curves plotted in the figure follow a simple analytical expression, in the form of , to which we have fitted the constants and ; cf. Supplement 1. Such a feature may be of extreme relevance in the experimental design of graphene-based CPPs, functioning as a “ruler” or “recipe” to fabricate graphene CPPs with tailored plasmonic properties. In fact, this relation is entirely general irrespective of the 2D material that is employed, since the specific properties of each material are all encapsulated within the flat , which is known; check Supplement 1. Therefore, although we have mainly focused on graphene throughout this paper, our results can be made completely general, as in Fig. 5, for any plasmonic 2D medium.
4. CONCLUSIONS
In conclusion, we have developed a universal framework to characterize plasmonic excitations propagating along the apex of a V-shaped graphene sheet forming a 1D channel. Using the EIM applied to the case of graphene, we derived the dispersion relation for graphene CPPs and their corresponding electric field distributions. As our calculations demonstrate, the EIM not only endows us with an intuitive picture of the physics governing these guided modes, but also can provide results consistent with the ones obtained using more elaborated theories, in particular for the fundamental channel plasmon mode (which in any case should be the most relevant from an applications perspective owing to its ability to reach deeper subwavelength regimes). This is courtesy of a simple yet powerful universal equation and corresponding scaling law that depends only on the system’s geometry (namely the angle ). Therefore, we believe our work may be applied as a first step—since the solution is essentially obtained instantaneously—in the design of novel devices for extreme subwavelength waveguiding of electromagnetic radiation [21] based either on graphene CPPs or on any other plasmon-supporting 2D material belonging to the now vast “2D library” of atomically thin materials [45].
Furthermore, in principle, one could envision the inclusion of a gain medium for loss-compensation in order to obtain longer propagation lengths. The effect of gain media could be incorporated in our model by replacing with an adequate dielectric function , so that it balances the Ohmic losses of the 2D material [46]. Note, however, that this is possible only provided that we are below the lasing threshold and for ; cf. Eq. (4). Naturally, a rigorous and complete description of gain can only be achieved using a quantum optical framework; see, for example, Ref. [47].
Finally, as a last remark, we note that our notion of a 2D material may in principle also include ultra-thin metallic films. For instance, a planar Drude-metal film of finite, but small, thickness has a symmetric mode with a low-frequency dependence of the form [48], thus being analogous to our generic 2D materials. In this way, CPPs in ultra-thin wedge/groove geometries [49,50] can also be covered by the universality outlined above.
Funding
Danmarks Grundforskningsfond (DNRF) (DNRF103); H2020 European Research Council (ERC) (341054 (PLAQNAP)); Fundação para a Ciência e a Tecnologia (FCT) (UID/FIS/04650/2013); European Commission (EC) (696656); Villum Fonden.
Acknowledgment
N. M. R. P. acknowledges funding from the European Commission within the project “Graphene-Driven Revolutions in ICT and Beyond” and the Portuguese Foundation for Science and Technology (FCT) in the framework of the Strategic Financing UID/FIS/04650/2013. S. I. B. acknowledges financial support by the European Research Council, and N. A. M. is a Villum Investigator supported by Villum Fonden.
See Supplement 1 for supporting content.
REFERENCES AND NOTES
1. A. A. Maradudin, W. L. Barnes, and J. R. Sambles, Modern Plasmonics (Elsevier, 2014).
2. D. K. Gramotnev and S. I. Bozhevolnyi, “Plasmonics beyond the diffraction limit,” Nat. Photonics 4, 83–91 (2010). [CrossRef]
3. W. L. Barnes, A. Dereux, and T. W. Ebbesen, “Surface plasmon subwavelength optics,” Nature 424, 824–830 (2003). [CrossRef]
4. E. Ozbay, “Plasmonics: merging photonics and electronics at nanoscale dimensions,” Science 311, 189–193 (2006). [CrossRef]
5. S. I. Bozhevolnyi, Plasmonic Nanoguides and Circuits (Pan Stanford, 2008).
6. S. I. Bozhevolnyi, “Effective-index modeling of channel plasmon polaritons,” Opt. Express 14, 9467–9476 (2006). [CrossRef]
7. E. Moreno, F. J. García-Vidal, S. G. Rodrigo, L. Martín-Moreno, and S. I. Bozhevolnyi, “Channel plasmon-polaritons: modal shape, dispersion, and losses,” Opt. Lett. 31, 3447–3449 (2006). [CrossRef]
8. T. W. Ebbesen, C. Genet, and S. I. Bozhevolnyi, “Surface-plasmon circuitry,” Phys. Today 61(5), 44–50 (2008). [CrossRef]
9. C. L. C. Smith, N. Stenger, A. Kristensen, N. A. Mortensen, and S. I. Bozhevolnyi, “Gap and channeled plasmons in tapered grooves: a review,” Nanoscale 7, 9355–9386 (2015). [CrossRef]
10. S. I. Bozhevolnyi, V. S. Volkov, E. Devaux, and T. W. Ebbesen, “Channel plasmon-polariton guiding by subwavelength metal grooves,” Phys. Rev. Lett. 95, 046802 (2005). [CrossRef]
11. T. Søndergaard, S. M. Novikov, T. Holmgaard, R. L. Eriksen, J. Beermann, Z. Han, K. Pedersen, and S. I. Bozhevolnyi, “Plasmonic black gold by adiabatic nanofocusing and absorption of light in ultra-sharp convex grooves,” Nat. Commun. 3, 969 (2012). [CrossRef]
12. S. Raza, N. Stenger, A. Pors, T. Holmgaard, S. Kadkhodazadeh, J. B. Wagner, K. Pedersen, M. Wubs, S. I. Bozhevolnyi, and N. A. Mortensen, “Extremely confined gap surface-plasmon modes excited by electrons,” Nat. Commun. 5, 4125 (2014). [CrossRef]
13. S. I. Bozhevolnyi, V. S. Volkov, E. Devaux, J.-Y. Laluet, and T. W. Ebbesen, “Channel plasmon subwavelength waveguide components including interferometers and ring resonators,” Nature 440, 508–511 (2006). [CrossRef]
14. C. Haffner, W. Heni, Y. Fedoryshyn, J. Niegemann, A. Melikyan, D. L. Elder, B. Baeuerle, Y. Salamin, A. Josten, U. Koch, C. Hoessbacher, F. Ducry, L. Juchli, A. Emboras, D. Hillerkuss, M. Kohl, L. R. Dalton, C. Hafner, and J. Leuthold, “All-plasmonic Mach-Zehnder modulator enabling optical high-speed communication at the microscale,” Nat. Photonics 9, 525–528 (2015). [CrossRef]
15. R. B. Nielsen, I. Fernandez-Cuesta, A. Boltasseva, V. S. Volkov, S. I. Bozhevolnyi, A. Klukowska, and A. Kristensen, “Channel plasmon polariton propagation in nanoimprinted V-groove waveguides,” Opt. Lett. 33, 2800–2802 (2008). [CrossRef]
16. D. K. Gramotnev and D. F. P. Pile, “Single-mode subwavelength waveguide with channel plasmon-polaritons in triangular grooves on a metal surface,” Appl. Phys. Lett. 85, 6323–6325 (2004). [CrossRef]
17. X. Li, T. Jiang, L. Shen, and X. Deng, “Subwavelength guiding of channel plasmon polaritons by textured metallic grooves at telecom wavelengths,” Appl. Phys. Lett. 102, 031606 (2013). [CrossRef]
18. D. F. P. Pile and D. K. Gramotnev, “Channel plasmon-polariton in a triangular groove on a metal surface,” Opt. Lett. 29, 1069–1071 (2004). [CrossRef]
19. O. Lotan, C. L. C. Smith, J. Bar-David, N. A. Mortensen, A. Kristensen, and U. Levy, “Propagation of channel plasmons at the visible regime in aluminum V-groove waveguides,” ACS Photon. 3, 2150–2157 (2016). [CrossRef]
20. V. S. Volkov, S. I. Bozhevolnyi, S. G. Rodrigo, L. Martín-Moreno, F. J. García-Vidal, E. Devaux, and T. W. Ebbesen, “Nanofocusing with channel plasmon polaritons,” Nano Lett. 9, 1278–1282 (2009). [CrossRef]
21. S. I. Bozhevolnyi and N. A. Mortensen, “Plasmonics for emerging quantum technologies,” Nanophotonics6 (2017), https://doi.org/10.1515/nanoph-2016-0179. [CrossRef]
22. E. Bermúdez-Ureña, C. Gonzalez-Ballestero, M. Geiselmann, R. Marty, I. P. Radko, T. Holmgaard, Y. Alaverdyan, E. Moreno, F. J. García-Vidal, S. I. Bozhevolnyi, and R. Quidant, “Coupling of individual quantum emitters to channel plasmons,” Nat. Commun. 6, 7883 (2015). [CrossRef]
23. L. Ju, B. Geng, J. Horng, C. Girit, M. Martin, Z. Hao, H. A. Bechtel, X. Liang, A. Zettl, Y. R. Shen, and F. Wang, “Graphene plasmonics for tunable terahertz metamaterials,” Nat. Nanotechnol. 6, 630–634 (2011). [CrossRef]
24. P. A. D. Gonçalves and N. M. R. Peres, An Introduction to Graphene Plasmonics (World Scientific, 2016).
25. S. Xiao, X. Zhu, B.-H. Li, and N. A. Mortensen, “Graphene-plasmon polaritons: from fundamental properties to potential applications,” Front. Phys. 11, 117801 (2016). [CrossRef]
26. F. J. García de Abajo, “Graphene plasmonics: challenges and opportunities,” ACS Photon. 1, 135–152 (2014). [CrossRef]
27. F. Xia, H. Wang, D. Xiao, M. Dubey, and A. Ramasubramaniam, “Two-dimensional material nanophotonics,” Nat. Photonics 8, 899–907 (2014). [CrossRef]
28. H. Yan, T. Low, W. Zhu, Y. Wu, M. Freitag, X. Li, F. Guinea, P. Avouris, and F. Xia, “Damping pathways of mid-infrared plasmons in graphene nanostructures,” Nat. Photonics 7, 394–399 (2013). [CrossRef]
29. Z. Fang, S. Thongrattanasiri, A. Schlather, Z. Liu, L. Ma, Y. Wang, P. M. Ajayan, P. Nordlander, N. J. Halas, and F. J. García de Abajo, “Gated tunability and hybridization of localized plasmons in nanostructured graphene,” ACS Nano 7, 2388–2395 (2013). [CrossRef]
30. H. Yan, F. Xia, Z. Li, and P. Avouris, “Plasmonics of coupled graphene micro-structures,” New J. Phys. 14, 125001 (2012). [CrossRef]
31. X. Zhu, W. Wang, W. Yan, M. B. Larsen, P. Bøggild, T. G. Pedersen, S. Xiao, J. Zi, and N. A. Mortensen, “Plasmon-phonon coupling in large-area graphene dot and antidot arrays fabricated by nanosphere lithography,” Nano Lett. 14, 2907–2913 (2014). [CrossRef]
32. A. Kumar, K. H. Fung, M. T. H. Reid, and N. X. Fang, “Transformation optics scheme for two-dimensional materials,” Opt. Lett. 39, 2113–2116 (2014). [CrossRef]
33. T. Christensen, A.-P. Jauho, M. Wubs, and N. A. Mortensen, “Localized plasmons in graphene-coated nanospheres,” Phys. Rev. B 91, 125414 (2015). [CrossRef]
34. M. Riso, M. Cuevas, and R. A. Depine, “Tunable plasmonic enhancement of light scattering and absorption in graphene-coated subwavelength wires,” J. Opt. 17, 075001 (2015). [CrossRef]
35. W. Wang, B.-H. Li, E. Stassen, N. A. Mortensen, and J. Christensen, “Localized surface plasmons in vibrating graphene nanodisks,” Nanoscale 8, 3809–3815 (2016). [CrossRef]
36. P. A. D. Gonçalves, E. J. C. Dias, S. Xiao, M. I. Vasilevskiy, N. A. Mortensen, and N. M. R. Peres, “Graphene plasmons in triangular wedges and grooves,” ACS Photon. 3, 2176–2183 (2016). [CrossRef]
37. D. Smirnova, S. H. Mousavi, Z. Wang, Y. S. Kivshar, and A. B. Khanikaev, “Trapping and guiding surface plasmons in curved graphene landscapes,” ACS Photon. 3, 875–880 (2016). [CrossRef]
38. P. Liu, X. Zhang, Z. Ma, W. Cai, L. Wang, and J. Xu, “Surface plasmon modes in graphene wedge and groove waveguides,” Opt. Express 21, 32432–32440 (2013). [CrossRef]
39. S. I. Bozhevolnyi and K. V. Nerkararyan, “Channel plasmon polaritons guided by graded gaps: closed-form solutions,” Opt. Express 17, 10327–10334 (2009). [CrossRef]
40. S. I. Bozhevolnyi and K. V. Nerkararyan, “Analytic description of channel plasmon polaritons,” Opt. Lett. 34, 2039–2041 (2009). [CrossRef]
41. In fact, more rigorously, for .
42. We remark that if the dielectric material within the V-shaped region differs from the cladding dielectric media, it is not possible to derive a simple, closed-form expression for . For this reason, here we consider only the case of a homogeneous dielectric environment.
43. In writing this equation we have explicitly approximated the width (for small angles) as , for the sake of clarity. The expression using the exact form is recoverable upon the replacement (we have used the exact expression in all our results).
44. W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical Recipes: the Art of Scientific Computing (Cambridge University, 2007).
45. A. K. Geim and I. V. Grigorieva, “Van der Waals heterostructures,” Nature 499, 419–425 (2013). [CrossRef]
46. S. Xiao, V. P. Drachev, A. V. Kildishev, X. Ni, U. K. Chettiar, H.-K. Yuan, and V. M. Shalaev, “Loss-free and active optical negative-index metamaterials,” Nature 466, 735–738 (2010). [CrossRef]
47. E. Amooghorban, N. A. Mortensen, and M. Wubs, “Quantum optical effective-medium theory for loss-compensated metamaterials,” Phys. Rev. Lett. 110, 153602 (2013). [CrossRef]
48. S. Raza, T. Christensen, M. Wubs, S. I. Bozhevolnyi, and N. A. Mortensen, “Nonlocal response in thin-film waveguides: loss versus nonlocality and breaking of complementarity,” Phys. Rev. B 88, 115401 (2013). [CrossRef]
49. S. Lee and S. Kim, “Long-range channel plasmon polaritons in thin metal film V-grooves,” Opt. Express 19, 9836–9847 (2011). [CrossRef]
50. X. Chen, N. C. Lindquist, D. J. Klemme, P. Nagpal, D. J. Norris, and S.-H. Oh, “Split-wedge antennas with sub-5 nm gaps for plasmonic nanofocusing,” Nano Lett. 16, 7849–7856 (2016). [CrossRef]