Robustness of the Rabi splitting under nonlocal corrections in plexcitonics
We explore theoretically how nonlocal corrections in the description of the metal affect the strong coupling between excitons and plasmons in typical examples where nonlocal effects are anticipated to be strong, namely small metallic nanoparticles, thin metallic nanoshells or dimers with narrow separations, either coated with or encapsulating an excitonic layer. Through detailed simulations based on the generalised nonlocal optical response theory, which simultaneously accounts both for modal shifts due to screening and for surface-enhanced Landau damping, we show that, contrary to expectations, the influence of nonlocality is rather limited, as in most occasions the width of the Rabi splitting remains largely unaffected and the two hybrid modes are well distinguishable. We discuss how this behaviour can be understood in view of the popular coupled-harmonic-oscillator model, while we also provide analytic solutions based on Mie theory to describe the hybrid modes in the case of matryoshka-like single nanoparticles. Our analysis provides an answer to a so far open question, that of the influence of nonlocality on strong coupling, and is expected to facilitate the design and study of plexcitonic architectures with ultrafine geometrical details.
The strong coupling of excitons in organic molecules and quantum dots with surface plasmons has been attracting increasing interest in recent years (1); (2). Both propagating surface plasmons in thin metallic films (3); (4); (5); (6); (7); (8); (9) or localised surface plasmons in metallic nanoparticles (10); (11); (12); (13); (14); (15); (16) have been explored, not only from the point of view of theory, for which it constitutes a new class of light-matter interactions characterised by unique, hybrid optical modes (17); (18); (19); (20); (21), but also because of its technological prospects in quantum optics and the emerging field of quantum plasmonics (22); (23) as well as the design of novel optical components (24); (25); (26); (27); (28).
Originally, strong coupling effects were studied for atomic and solid-state systems inside optical or photonic cavities (29); (30); (31); (32). -aggregates of organic molecules are however increasingly more frequently preferred, as they are characterised by large dipole moments and narrow transition lines (11); (33). At the same time, plasmonic architectures are excellent templates as cavities for the experimental realisation of hyprid exciton-photon systems, because they provide nanoscale confinement and small modal volumes (34), thus permitting even single-molecule strong coupling at room temperature (35); (36). The interaction of excitons with plasmons is usually traced through the Rabi splitting in optical spectra (4); (11), but recent elaborate experiments showed that it is also possible to observe the coherent energy exchange between the optical states of the two components in real time through the corresponding Rabi oscillations(37); (38); (39). Tailoring the plasmonic environment is gradually departing from the regime of single nanoparticles (12); (13); (40); (41), and dimers (42); (43); (44); (45); (46); (47); (48) or nanoparticle arrays (49); (50); (51); (52); (53) are exploited to benefit from their stronger field confinement and enhancement. A similar coupling has also been observed for the interaction of emitters with the phononic modes of SiC antennae (54), while engineering the electromagnetic vacuum through evanescent plasmon modes was recently proposed as a promising route towards increased light-matter interactions (55).
As novel plasmonic architectures are explored in the quest for stronger, possibly even ultrastrong coupling (56); (57); (58), reducing characteristic lengths such as nanoparticle sizes and/or separations is a natural step that goes hand-in-hand with advances in nanofabrication. In such systems however, the excitation of a richness of hybrid optical modes characterised by large field intensities (59); (60) is accompanied by the increased influence of nonclassical effects that require a description beyond classical electrodynamics (61). Nonlocal screening due to the spatial dispersion of induced charges (62); (63); (64); (65), enhanced Landau damping near the metal surface (66); (67); (68), electron spill-out (69); (70); (71); (72); (73); (74) and quantum tunnelling (75); (76); (77) are different manifestations of the quantum-mechanical nature of plasmonic nanostructures in the few-nm regime, which need to be considered for an accurate evaluation of the emitter-plasmon coupling. For instance, nonlocal effects strongly affect spontaneous emission rates for quantum emitters in plasmonic environments (78), and lead to reduced near-fields as compared to the common local-response approximation (LRA) (79). Fluorescence of molecules in the vicinity of single metallic nanoparticles (80) or inside plasmonic cavities (81) is also affected by a combination of nonlocal induced-charge screening and surface-enhanced Landau damping. In the context of strong coupling, Marinica et al. have reported quenching of the plexcitonic strength in dimers, attributed to quantum tunnelling (82). Even before entering the tunnelling regime, one anticipates that experimentally observed nonlocal modal shifts and broadening (83); (60); (84) can manifest themselves as deviations in the dispersion diagram and changes in the width of the Rabi splitting in the optical spectra. Nevertheless, related studies are still missing and these expectations are yet to be confirmed.
In the recent review by Törmä and Barnes (1), one section was devoted on how smaller dimensions of the plasmon-supporting structures become increasingly more important. In particular, it was anticipated that one consequence of nonlocal effects that is of great interest in the context of strong coupling is that of a reduction in the field enhancement that can be achieved when light is confined to truly nm dimensions. Motivated by this discussion, and also by our recent studies of single emitters in plasmonic environments (80); (81), we explore here the influence of nonlocality on plexcitonics through theoretical calculations based on the generalised nonlocal optical response (GNOR) theory (85). GNOR has proven particularly efficient in recent years in simultaneously describing both nonlocal screening and Landau damping through a relatively simple correction of the wave equation that is straightforward to implement to any plasmonic geometry model (86). Here we take advantage of its flexibility to describe the optical response of typical plexcitonic architectures: metallic nanoshells either coated by an excitonic layer or containing it as nanoparticle core, dimers of such nanoshells and homogeneous nanosphere dimers encapsulated in an excitonic matrix. In all situations it is shown that, while nonlocal modal shifts introduce a detuning for fixed geometrical details, the width of the Rabi splitting remains in practice unaffected if one follows the more practical approach of modifying nanoparticle sizes and/or separations so as to tune the plasmon mode to the exciton energy. Broadening of the plasmon modes due to Landau damping is also shown to be weak enough, so as to secure that the two hybrid modes remain distinguishable, and linewidth-based criteria for strong coupling are satisfied. This somehow unexpected result is interpreted in view of a popular coupled-harmonic-oscillator (CHO) model and the strength of the coupling as related to the modal volume, while analytic solutions for the hybrid modes are obtained based on Mie theory. Our analysis relaxes concerns about nonlocal effects in the coupling of plasmons with excitons, and provides additional flexibility to the experimental realisation of novel plexcitonic architectures with sizes in the few-nm regime.
Ii Results and discussion
In order to facilitate a strong influence of nonlocality, we consider throughout this paper small nanoparticles with radii that do not exceed 20 nm, and thin metallic shells of 1–2 nm width. While such thin nanoshells are still experimentally challenging, successful steps towards this direction have been presented recently (87). The choice of nanoshells over homogeneous metallic spheres serves a dual purpose. On the one hand, thin metallic layers ensure the reduced lengths required for a strong nonlocal optical response, even though the far-field footprint of nonlocality can often prove negligible (88); (89). On the other hand, and more importantly, modifying the shell thickness provides the required plasmon tuning (90) to match the exciton energy and observe the anticrossing of the two hybrid modes in dispersion diagrams (12); (58). As discussed by Törmä and Barnes (1), metals with low loss are beneficial for strong-coupling applications, because the plasmon modes should have linewidths comparable to those of the excitonic material. In our case, low intrinsic loss also ensures that the role of Landau damping will be better illustrated. For this reason, silver is adopted as the metallic material throughout the paper.
As a first example of a nonlocal plexcitonic system, we study a silver nanoshell of total radius 5 nm in air, encapsulating an organic-dye core of variable radius , so that the metal thickness is , as shown in the schematics of Fig. 1. The dye is modelled as a homogeneous excitonic layer with a Drude–Lorentz dielectric function as described in Sec. IV, with its transition energy at 2.7 eV. The extinction () spectrum for a homogeneous sphere of radius 5 nm described by such a permittivity is shown in Fig. 1(a) by the black line (magnified in the inset), and is indeed characterised by a weak resonance at 2.7 eV. With a red (blue) dashed line is depicted the corresponding LRA (GNOR) extinction spectrum for a hollow silver nanoshell (in the absence of the excitonic core) with 5 nm, and 0.94 nm, the thickness for which the symmetric particle-like plasmon mode of the nanoshell (91) matches in the local description. As expected, when nonlocal effects are taken into account through GNOR, for the same geometrical parameters the plasmon mode experiences both a blueshift, of about 0.023 eV, and additional broadening, so that the full width at half maximum of the resonance increases from 0.05 eV to 0.067 eV. Once the dye core is included, its excitonic mode interacts with the localised surface plasmon, and two hybrid modes are formed, as shown by solid red and blue lines in Fig. 1(a), with a Rabi splitting of energy 0.156 eV. While both peaks experience a nonlocal blueshift, its strength is practically the same, so that remains unaffected. However, one should be cautious and not draw conclusions just from Fig. 1(a), since the spectra have been calculated for a shell thickness for which the plasmon mode is detuned from in the nonlocal description.
In order to get a better visualisation of the strong coupling of Fig. 1(a), we plot in Fig. 1(b) the dispersion diagram of the two modes as a function of the plasmon energy , calculated by modifying the shell thickness and obtaining the resonance peak for hollow nanoshells. We then repeat the calculation in the presence of the dye, and obtain the frequencies of the hybrid modes, , from the extinction spectra. Plotting the energies of the two hybrid modes as a function of leads to dispersion lines that are practically indistinguishable for the LRA (red lines) and GNOR (blue lines) models. The same dispersion diagram is plotted in the inset of Figure 1b as a function of shell thickness. It is straightforward to see that the GNOR dispersion lines are just horizontally shifted with respect to the corresponding LRA ones, as the nonlocal blueshift makes thicker shells necessary to achieve the same plasmon energy.
Together with the numerical results, in Fig. 1(b) we also present the solutions obtained through two different analytic approaches: solid red circles represent the hybrid mode frequencies obtained from the CHO model, while open triangles correspond to analytic solutions based on Mie theory – both are in excellent agreement with the full simulations. Due to its simplicity, the oscillator model is frequently adopted in literature to describe the modes in strongly-coupled plexcitonic geometries (7); (42); (45). In the general case, the (complex) energies of the two hybrid modes are given by (1); (16)
where and are the uncoupled plasmon and exciton damping rates, respectively. The frequency of the Rabi splitting is usually obtained by the simulated or experimental spectra at the crossing point, and it is directly related to the coupling strength through . When absorptive losses are low, and the plasmon and exciton linewidths are similar, Eq. (II) simplifies to (92)
We note here that, since the width of the Rabi splitting is itself often obtained from the dispersion diagrams, instead of more elaborate, quantum-mechanical analyses based on the coupling strength, it is not that surprising that Eq. (2) reproduces the same data so well. On the other hand, for spherical particles one can obtain the position of the modes through the poles of the scattering matrix in the Mie solution (81), an approach which is the nanoparticle equivalent to introducing the metal and excitonic dielectric functions into the surface plasmon polariton dispersion relation (1). This was done for example by Fofang et al. (12) in the simplified case where many of the (background) dielectric constants involved could be taken equal to unity. To generalise this, the dipolar Mie eigenfrequencies of a core-shell particle can be obtained through (81)
where is the dielectric function of the core, the corresponding one of the shell, and describes the environment, as shown in the schematics of Fig. 1. Using the Drude–Lorentz model discussed in the Methods section for , and a simple Drude model for (with plasma frequency equal to 9 eV and background dielectric constant equal to 3.65 to model silver), and disregarding absorptive losses in both materials as we are interested in the real resonance frequencies, leads to a relation in which is raised to the power of 6, and therefore six eigenfrequencies are obtained. Three of them are negative and have no physical meaning; the other three are the two hybrid plasmon-exciton modes, and , shown by open triangles in Fig. 1(b), and the antisymmetric, cavity-like plasmon mode of the nanoshell (91), which is at much higher frequencies and does not interact with the exciton.
Apart from the resonance frequency of the two hybrid modes and the width of the Rabi splitting, of particular interest is also the linewidth of the initial uncoupled modes, and that of the hybrid ones. In the description of mechanical harmonic oscillators, the condition for reaching the strong coupling regime is usually defined as , where is the largest frequency between the two of the original uncoupled modes. In other contexts, this condition describes the ultrastrong coupling regime. In plasmonics these requirements are usually relaxed, and strong coupling is defined in a more pragmatic way, through (1)
This is the condition required for the Rabi splitting to be observable in the spectra, and it indeed holds true in the example studied in Fig. 1, as can be verified by introducing the linewidths of the uncoupled modes discussed above into Eq (4). To better illustrate the broadening and damping of the modes predicted by nonlocal theories, in Figs. 1(c)-(d) we show contour plots of extinction versus energy, obtained with the LRA and GNOR model, respectively. For all shell thicknesses (and therefore the corresponding plasmon energies) studied here, the two hybrid modes are clearly distinguishable, and the relaxed strong-coupling condition of Eq. (4) holds. We have repeated the same analysis for several nanoparticle sizes and radius–thickness combinations, even reproducing results in the ultra-strong coupling regime (58), without observing qualitative or quantitative differences from the above discussion.
As a second example, it is natural to explore the inverse topology, where the excitonic layer grows around a metallic nanoshell, a geometry which is easier to achieve experimentally and has been shown to be preferable for obtaining large Rabi splittings and even entering the ultrastrong coupling regime (58). In Fig. 2 we replace the excitonic core of the nanoshell of Fig. 1 with a SiO one (described by a dielectric constant ). We assume a constant nanoshell radius 5 nm, covered by an external excitonic layer of thickness 5 nm and exciton energy 2.7 eV as previously, and modify the thickness of the nanoshell, as shown in the schematics of Fig. 2. The extinction spectra of the uncoupled components are shown in Fig. 2(a), with a black line for the spectrum of a homogeneous excitonic sphere of 10 nm, and with red (blue) dashed lines for the LRA (GNOR) spectra of a silver nanoshell with 1.22 nm. A thicker shell than in Fig. 1(a) is required to bring the plasmon resonance at 2.7 eV, since a higher-index dielectric is now used as the core material. The corresponding coupled spectra are shown by solid lines, and are characterised by a Rabi splitting 0.206 eV. In addition to the two hybrid modes, a third peak of much lower intensity is visible at 2.7 eV, which does not shift at all under nonlocal corrections. This is the surface exciton polariton mode excited at the outer surface of the dye layer, which indeed has no reason to be affected by the plasmon modal shifts (93). It is always excited for this kind of nanoparticle topology, but it is usually masked by the much stronger peaks of the hybrid modes due to the large size of the plasmonic particles. The small nanoshell dimensions explored here ensure that it is clearly discernible in the extinction spectra. Such modes were recently explored as substitutes for plasmons in metals for applications in nanophotonics (94).
In Fig. 2(b) we present the corresponding dispersion diagram, as a function of the plasmon energy, similarly to Fig. 1(b). Given our previous analysis in relation to Fig. 1(b), and the spectra presented in Fig. 2(a), it is not surprising that the dispersion lines for the LRA and GNOR models are again indistinguishable. The CHO model captures excellently the two hybrid modes, as shown by the red solid dots, but it lacks a description of the uncoupled excitonic mode. This weak resonance is excellently reproduced however by the corresponding Mie-based analytic solutions. Since in this case we study a three-layer nanoparticle, one should in principle obtain the corresponding scattering matrix (95); (96) and use appropriate asymptotic expressions (97) to obtain the resonances from its poles. A simpler approximation, which we successfully adopt here, is to separate the nanoparticle into two components: a metallic nanoshell embedded in an infinite excitonic environment described by a dielectric function , for which the hybrid modes can be obtained from Eq. (3) (replacing with ), and an excitonic sphere of in an infinite environment described by , for which the resonance condition is (81)
As can be seen by the open triangles in Fig. 2(b), due to the small sizes involved, this approach works extremely well to accurately describe not only the hybrid modes , but the uncoupled excitonic mode as well. Finally, in Figs. 2(c)-(d) it is shown through full extinction contours that the additional nonlocal broadening of the modes is not significant in this case either. Considering the initial, uncoupled linewidths of the plasmon resonances shown in Fig. 2(a) (0.049 eV for LRA and 0.072 eV for GNOR), Eq. (4) again implies that strong coupling is still achievable under nonlocal corrections.
In search for a geometry where the footprint of nonlocality might be stronger, we depart from the description of single nanoparticles, and study in Figure 3 an exciton core–silver shell dimer. The particles have a total radius 20 nm and a metallic shell of thickness 2 nm, and are separated by a variable gap of width . Both the thin shell thickness and the narrow dimer gap are expected to enhance nonlocal effects in this case. Instead of tuning the plasmon through the shell thickness, as in Figs 1 and 2, here we change the width of the gap. Plasmon hybridisation (91) causes the bonding dimer plasmon modes to redshift as the gap shrinks (59), allowing for efficient tuning of the optical response. Extinction contours obtained with the LRA and GNOR models as a function of are shown in Figs. 3(a) and (b), respectively, for an excitonic core with 1.7 eV. The vertical dashed lines denote the gap for which the bonding dimer plasmon mode in the absence of a dye core is at 1.7 eV, and indicate a 0.22 nm difference between the two models. In Fig. 3(c) extinction spectra for each model are presented both for this gap width within the corresponding model (solid lines), and also for the gap which tunes the plasmon at 1.7 eV within the other model (dashed lines). Examining the tuned situation for each model, it is evident that no change in the width of the Rabi splitting and no significant broadening can be observed for this geometry either. This is more clear in the dispersion diagram of Fig. 3(d), where the energy of the hybrid modes is plotted versus the plasmon energy in each model, leading once more to almost identical lines, also in good agreement with the CHO model.
Finally, we turn again to the inverse topology, of a dimer enclosed in an excitonic matrix, as shown schematically on the right-hand side of Fig. 4. To avoid repetitions and make things more interesting, instead of nanoshells we consider homogeneous silver spheres of radius 15 nm separated by a gap of 1 nm, covered with an excitonic layer made of two spheres of radius 20 nm, overlapping around the gap as they are concentric with the corresponding silver particles. In the absence of the excitonic layer the plasmon resonance is found to be at 2.92 eV within LRA, and at 2.97 eV within GNOR. Instead of tuning the dimer plasmon through some geometrical parameter, we keep in this case the geometry fixed and modify instead, an approach often adopted in theoretical calculations due to its convenience (45). Extinction contour maps obtained with LRA and GNOR are shown in Figs. 4(a) and (b) respectively, together with the uncoupled excitonic and plasmonic lines (white dashed lines). The corresponding extinction spectra in the presence (solid lines) or absence (dashed lines) of the excitonic layer are shown in Fig. 4(c). One again observes that the width of the Rabi splitting appears to be nearly unaffected, but the hybrid modes are significantly broadened. Furthermore, the uncoupled excitonic peak, which is strong within LRA, nearly vanishes within GNOR, and can only be traced through the corresponding scattering spectra. Nevertheless, this strong damping, originating mainly from the presence of a much larger amount of metallic material (homogeneous spheres instead of shells), is still not enough to bring the system outside the strong coupling regime as defined by Eq. (4). In addition, careful examination of the dispersion diagrams in Fig. 4(d) shows that the Rabi splitting is in fact narrower by 0.012 eV in the case of GNOR, finally identifying a small impact of nonlocality on the plexcitonic response.
In all the examples explored above, we focused on the dominant, dipolar plasmon modes, which are more relevant from an experimental point of view. While we have identified situations in which our conclusions hold for higher-order modes as well, the coupling of excitons with e.g. quadrupolar single-particle or bonding-dimer modes, and how it is affected by nonlocal effects, relies on several parameters, and a more systematic study is required. Efficient excitation of higher-order modes in single nanoparticles calls for larger particle sizes, for which nonlocal effects tend to become negligible (60). On the other hand, for very thin metallic shells, nonlocality can wash out higher-order modes (80), so that a classically predicted Rabi splitting could disappear completely. Similarly, the efficient excitation of quadrupolar bonding-dimer modes might require to enter the truly sub-nm regime (59), where quantum tunnelling prevails, restricting nonlocality to a secondary role. Further theoretical work in the future should shed more light on such issues.
To better understand the reason why nonlocal frequency shifts do not usually affect the width of the Rabi splitting, one can recall once more the CHO model of Eq. (2), according to which the energy difference between the two hybrid modes is
where . Nonlocal effects introduce a small correction to the plasmon, . In the strong-coupling regime, where , and assuming at a first step that remains constant for simplicity, Eq. (6) becomes
For Rabi splittings around 0.2 eV, and nonlocal blueshifts of the order of 0.03 eV, typical values in the examples explored here, this introduces a change in of about 1%, which is hard to observe. According to the above analysis, one needs to identify systems with much larger nonlocal frequency shifts, accompanied with narrow Rabi splittings, to obtain a strong effect. What is still missing however is the influence of nonlocality on the coupling strength, and therefore, through , on the Rabi splitting itself. In a semi-classical approach, and as long as the exact positions and orientations of individual emitters in the excitonic layer can be disregarded, the coupling strength is given by (15)
where is the number of emitters characterised by a transition dipole moment , is the transition frequency, the vacuum permittivity and the mode volume. For spherical particles, it has been shown analytically that, for the dipolar mode which is of interest here, the mode volume depends only on the the geometrical volume and the environment through (98). Nevertheless, this result was obtained assuming local response functions for the metal, and nonlocal corrections to should be evaluated through the general definition (99)
where is the electromagnetic energy density at position . Introducing nonlocal corrections into Eq. (9) has been discussed in detail by Toscano et al. (100). Our calculations have shown that within the GNOR model typically increases by no more than 20%, which translates into a decrease in the coupling strength. Since this modification enters Eq. (6) through a correction, the nonlocal changes in coupling strength and plasmon energy tend, to first order, to cancel each other out. This tendency is further strengthened once absorptive losses are taken into account, which, according to Eq. (II) introduce an extra damping correction , working towards the same direction as in Eq. (7), to further counterbalance the (relatively stronger) effect of increasing the mode volume. In view of the above discussion, interpreting the results of Figs. 1–4 and understanding why nonlocal effects do not usually play an important role in the exciton-plasmon coupling can be achieved in a simple and intuitive manner.
In summary, we explored the influence of nonlocal effects in the description of the metal on the coupling of plasmonic nanoparticles with dye layers characterised by an excitonic transition. Through detailed simulations, in conjunction with analytical modelling, we showed that, contrary to expectations based mainly on results for single emitters in plasmonic environments, neither nonlocal frequency shifts due to screening, nor surface-enhanced Landau damping produces strong deviations from a description within classical electrodynamics. Apart from extreme situations of dramatic nonlocal blueshifts in situations of already narrow Rabi splittings, nonlocality affects plasmon-exciton coupling only incrementally, as long as the excitonic material is sufficiently described as a homogeneous layer. For simplicity, we have used here the GNOR model, leaving quantum corrections based on the Feibelman parameters for the centroid of charge or on surface-dipole moments for future investigation (77); (74). By analysing how nonlocal corrections in plasmonics enter into the common coupled-oscillator model, we provided a simple, intuitive interpretation of our findings. Novel plexcitonic architectures involving reduced, few-nm length scales, can be analysed without the necessity to resort to elaborate models going beyond classical electrodynamics, thus providing additional flexibility in the design and optimisation of systems with strong light-matter interactions.
Homogeneous silver nanoparticles and thin silver shells are described by the experimental dielectric function of Johnson and Christy (101). When nonlocal corrections are introduced in the metal, the Drude part is subtracted from to produce the background contribution to the metal dielectric function, , according to
where for silver we use 8.99 eV, 0.025 eV (65).
In the case of single nanoparticles, we include nonlocal corrections in the scattering matrix of the Mie solution as described in Ref. (102). For a homogeneous metallic sphere of radius described by a dielectric function in a homogeneous medium of the electric-type nonlocal Mie coefficients of multipole order , , are given by
where and are the spherical Bessel and first-type Hankel functions, respectively, while with being the (transverse) wave number in medium and denotes derivation with respect to the argument. The nonlocal correction to the Mie coefficients is given as
where and is the longitudinal wave number in the sphere, associated with the longitudinal dielectric function , which is frequency- and wave vector-dependent (103); (62). The dispersion of longitudinal waves is given by . In the limiting case where we retrieve the local result of standard Mie theory. The corresponding analytical solution is lengthy and can be found in the Supplementary Information of Ref. (80).
In the case of dimers, we use a commercial finite-element solver (Comsol Multiphysics 5.0) (104), appropriately adapted to include the description of nonlocal effects (105), to solve the system of coupled equations (85)
where and are the electric field and the induced current density respectively; is the Drude conductivity, and is the vacuum permeability, related to the velocity of light in vacuum through .
In both approaches, the hydrodynamic parameter is taken equal to , where m s is the Fermi velocity of silver (65), while for the diffusion constant we use m s (102). Values of this order of magnitude reproduce well the experimentally observed size-dependent broadening in small metallic nanoparticles (106). As the additional boundary condition we adopt the usual condition of zero normal component of the induced current density at the metal boundary (65), which implies a hard-wall description of the metal, an approach which might prove inefficient in the case of good jellium metals (107); (108), but is reasonable for a noble metal with high work function like silver. When describing the interface between silver and the excitonic layer, we assume that carriers from the one medium are not allowed to enter the other, an assumption justified by the fact that the excitonic layer consists of an assembly of dye molecules, and is only effectively described as a homogeneous layer.
For the excitonic material we use a Drude–Lorentz model according to
where is the reduced oscillator strength. Throughout the paper is allowed to vary as stated in each case, 0.02, and 0.052 eV (12).
Acknowledgements.C. T. was supported by funding from the People Programme (Marie Curie Actions) of the European Union’s Seventh Framework Programme (FP7/2007-2013) under REA grant agreement number 609405 (COFUNDPostdocDTU). We gratefully acknowledge support from VILLUM Fonden via the VKR Centre of Excellence NATEC-II and from the Danish Council for Independent Research (FNU 1323-00087). The Center for Nanostructured Graphene (CNG) was financed by the Danish National Research Council (DNRF103). N.A.M. is a VILLUM Investigator supported by VILLUM Fonden.
- P. Törmä and W. L. Barnes, Rep. Prog. Phys. 78, 013901 (2015).
- R. J. Moerland, T. K. Hakala, J.-P. Martikainen, H. T. Rekola, A. I. Väkeväinen, and P. Törmä, Strong coupling between organic molecules and plasmonic nanostructures, in Quantum Plasmonics, Eds. S. I. Bozhevolnyi, L. Martín-Moreno, and F. J. García-Vidal (Springer, Cham, 2017).
- I. Pockrand, A. Brillante, and D. Möbius, J. Chem. Phys. 77, 6289 (1982).
- J. Bellessa, C. Bonnand, J. C. Plenet, and J. Mugnier, Phys. Rev. Lett. 93, 036404 (2004).
- J. Dintinger, S. Klein, F. Bustos, W. L. Barnes, and T. W. Ebbesen, Phys. Rev. B 71, 035424 (2005).
- T. K. Hakala, J. J. Toppari, A. Kuzyk, M. Pettersson, H. Tikkanen, H. Kunttu, and P. Törmä, Phys. Rev. Lett. 103, 053602 (2009).
- D. E. Gómez, K. C. Vernon, P. Mulvaney, and T. J. Davis, Nano Lett. 10, 274 (2010).
- A. González-Tudela, P. A. Huidobro, L. Martín-Moreno, C. Tejedor, and F. J. García-Vidal, Phys. Rev. Lett. 110, 126801 (2013).
- T. U. Tumkur, G. Zhu, and M. A. Noginov, Opt. Express 24, 3921 (2016).
- G. P. Wiederrecht, G. A. Wurtz, and J. Hranisavljevic, Nano Lett. 4, 2121 (2004).
- Y. Sugawara, T. A. Kelf, J. J. Baumberg, M. E Abdelsalam, and P. N. Bartlett, Phys. Rev. Lett. 97, 266808 (2006).
- N. T. Fofang, T.-H. Park, O. Neumann, N. A. Mirin, P. Nordlander, and N. J. Halas, Nano Lett. 8, 3481 (2008).
- A. Trügler and U. Hohenester, Phys. Rev. B 77, 115403 (2008).
- G. Zengin, G. Johansson, P. Johansson, T. J. Antosiewicz, M. Käll, and T. Shegai, Sci. Rep. 3, 3074 (2013).
- G. Zengin, M. Wersäll, S. Nilsson, T. J. Antosiewicz, M. Käll, and T. Shegai, Phys. Rev. Lett. 114, 157401 (2015).
- D. Melnikau, R. Esteban, D. Savateeva, A. Sánchez-Iglesias, M. Grzelczak, M. K. Schmidt, L. M. Liz-Marzán, J. Aizpurua, and Y. P. Rakovich, J. Phys. Chem. Lett. 7, 354 (2016).
- D. G. Lidzey, D. D. C. Bradley, M. S. Skolnick, T. Virgili, S. Walker, and D. M. Whittaker, Nature 395, 53 (1998).
- D. E. Chang, A. S. Sørensen, P. R. Hemmer, and M. D. Lukin, Phys. Rev. B 76, 035420 (2007).
- A. Delga, J. Feist, J. Bravo-Abad, and F. J. García-Vidal, Phys. Rev. Lett. 112, 253601 (2014).
- J. George, S. Wang, T. Chervy, A. Canaguier-Durand, G. Schaeffer, J.-M. Lehn, J. A. Hutchison, C. Genet, and T. W. Ebbesen, Faraday Discuss. 178, 281 (2015).
- E. K. Tanyi, H. Thuman, N. Brown, S. Koutsares, V. A. Podolskiy, and M. A. Noginov, Adv. Opt. Mater. 5, 1600941 (2017).
- Quantum Plasmonics, Eds. S. I. Bozhevolnyi, L. Martín-Moreno, and F. J. García-Vidal (Springer, Cham, 2017).
- S. I. Bozhevolnyi and N. A. Mortensen, Nanophotonics, doi: 10.1515/nanoph-2016-0179 (2017).
- P. A. Hobson, S. Wedge, J. A. E. Wasey, I. Sage, and W. L. Barnes, Adv. Mater. 14, 1393 (2002).
- J. Dintinger, I. Robel, P. V. Kamat, C. Genet, and T. W. Ebbesen, Adv. Mater. 18, 1645 (2006).
- D. E. Chang, A. S. Sørensen, E. A. Demler, and M. D. Lukin, Nature Phys. 3, 807 (2007).
- A. Auffèves, J.-M. Gérard, and J.-P. Poizat, Phys. Rev. A 79, 053838 (2009).
- T. Schwartz, J. A. Hutchison, C. Genet, and T. W. Ebbesen, Phys. Rev. Lett. 106, 196405 (2011).
- R. J. Thompson, G. Rempe, and H. J. Kimble, Phys. Rev. Lett. 68, 1132 (1992).
- T. Yoshie, A. Scherer, J. Hendrickson, G. Khitrova, H. M. Gibbs, G. Rupper, C. Ell, O. B. Shchekin, and D. G. Deppe, Nature 432, 200 (2004).
- T. Aoki, B. Dayan, E. Wilcut, W. P. Bowen, A. S. Parkins, T. J. Kippenberg, K. J. Vahala, and H. J. Kimble, Nature 443, 671 (2006).
- A. Faraon, I. Fushman, D. Englund, N. Stoltz, P. Petroff, and J. Vučković, Nature Phys. 4, 859 (2008).
- D. G. Lidzey, D. D. C. Bradley, A. Armitage, S. Walker, and M. S. Skolnick, Science 288, 1620 (2000).
- J. A. Schuller, E. S. Barnard, W. Cai, Y. C. Jun, J. S. White, and M. L. Brongersma, Nature Mater. 9, 193 (2010).
- K. Santhosh, O. Bitton, L. Chuntonov, and G. Haran, Nature Commun. 7, 11823 (2016).
- R. Chikkaraddy, B. de Nijs, F. Benz, S. J. Barrow, O. A. Scherman, E. Rosta, A. Demetriadou, P. Fox, O. Hess, and J. J. Baumberg, Nature 535, 127 (2016).
- P. Vasa, W. Wang, R. Pomraenke, M. Lammers, M. Maiuri, C. Manzoni, G. Cerullo, and C. Lienau, Nature Photon. 7, 128 (2013).
- P. Vasa, W. Wang, R. Pomraenke, M. Maiuri, C. Manzoni, G. Cerullo, and C. Lienau, Phys. Rev. Lett. 114, 036802 (2015).
- H. Wang, H.-Y. Wang, A. Bozzola, A. Toma, S. Panaro, W. Raja, A. Alabastri, L. Wang, Q.-D. Chen, H.-L. Xu, F. De Angelis, H.-B. Sun, and R. P. Zaccaria, Adv. Funct. Mater. 26, 6198 (2016).
- E. Waks and D. Sridharan, Phys. Rev. A 82, 043945 (2010).
- R. Liu, Z.-K. Zhou, Y.-C. Yu, T. Zhang, H. Wang, G. Liu, Y. Wei, H. Chen, and X.-H. Wang, Phys. Rev. Lett. 118, 237401, (2017).
- S. Savasta, R. Saija, A. Ridolfo, O. Di Stefano, P. Denti, and F. Borghese, ACS Nano 4, 6369 (2010).
- A. Manjavacas, F. J. García de Abajo, and P. Nordlander, Nano Lett. 11, 2318 (2011).
- A. E. Schlather, N. Large, A. S. Urban, P. Nordlander, and N. J. Halas, Nano Lett. 13, 3281 (2013).
- O. Pérez-González, N. Zabala, and J. Aizpurua, Nanotechnology 25, 035201 (2014).
- E.-M. Roller, C. Argyropoulos, A. Högele, T. Liedl, and M. Pilo-Pais, Nano Lett. 16, 5962 (2016).
- , R.-Q. Li, D. Hernángomez-Pérez, F. J. García-Vidal, and A. I. Fernández-Domínguez, Phys. Rev. Lett. 117, 107401 (2016).
- X. Chen, Y.-H. Chen, J. Qin, D. Zhao, B. Ding, R. J. Blaikie, and M. Qiu, Nano Lett. 17, 3246 (2017).
- G. A. Wurtz, P. R. Evans, W. Hendren, R. Atkinson, W. Dickson, R. J. Pollard, and A. V. Zayats, Nano Lett. 7, 1297 (2007).
- E. Eizner, O. Avayu, R. Ditcovski, T. Ellenbogen, Nano Lett. 15, 6215 (2015).
- A. Tsargorodska, M. L. Cartron, C. Vasilev, G. Kodali, O. A. Mass, J. J. Baumberg, P. L. Dutton, C. N. Hunter, P. Törmä, and G. J. Leggett, Nano Lett. 16, 6850 (2016).
- F. Todisco, M. Esposito, S. Panaro, M. De Giorgi, L. Dominici, D. Ballarini, A. I. Fernández-Domínguez, V. Tasco, M. Cuscunà, A. Passaseo, C. Ciracì, G. Gigli, and D. Sanvitto, ACS Nano 10, 11360 (2016).
- J. Li, K. Ueno, H. Uehara, J. Guo, T. Oshikiri, and H. Misawa, J. Phys. Chem. Lett. 7, 2786 (2016).
- R. Esteban, J. Aizpurua, and G. W. Bryant, New J. Phys. 16, 013052 (2014).
- J. Ren, Y. Gu, D. Zhao, F. Zhang, T. Zhang, and Q. Gong, Phys. Rev. Lett. 118, 073604 (2017).
- M. Geiser, F. Castellano, G. Scalari, M. Beck, L. Nevou, and J. Faist, Phys. Rev. Lett. 108, 106402 (2012).
- S. Balci, Opt. Lett. 38, 4498 (2013).
- A. Cacciola, O. Di Stefano, R. Stassi, R. Saija, and S. Savasta, ACS Nano 8, 11483 (2014).
- I. Romero, J. Aizpurua, G. W. Bryant, and F. J. García de Abajo, Opt. Express 14, 9988 (2006).
- S. Raza, S. Kadkhodazadeh, T. Christensen, M. Di Vece, M. Wubs, N. A. Mortensen, and N. Stenger, Nature Commun. 6, 8788 (2015).
- W. Zhu, R. Esteban, A. G. Borisov, J. J. Baumberg, P. Nordlander, H. J. Lezec, J. Aizpurua, and K. B. Crozier, Nature Commun. 7, 11495 (2016).
- F. J. García de Abajo, J. Phys. Chem. C 112, 17983 (2008).
- J. M. McMahon, S. K. Gray, and G. C. Schatz, Phys. Rev. Lett. 103, 097403 (2009).
- S. Raza, G. Toscano, A.-P. Jauho, M. Wubs, and N. A. Mortensen, Phys. Rev. B 84, 121412(R), (2011).
- S. Raza, S. I. Bozhevolnyi, M. Wubs, and N. A. Mortensen, J. Phys.: Condens. Matter 27, 183204 (2015).
- X. Li, D. Xiao, and Z. Zhang, New J. Phys. 15, 023011 (2013).
- J. B. Khurgin, Faraday Discuss. 178, 109 (2015).
- T. V. Shahbazyan, Phys. Rev. B 94, 235431 (2016).
- G. Weick, G.-L. Ingold, R. A. Jalabert, and D. Weinmann, Phys. Rev. B 74, 165421 (2006).
- J. A. Scholl, A. L. Koh, and J. A. Dionne, Nature 483, 421 (2012).
- R. C. Monreal, T. J. Antosiewicz, and S. P. Apell, New J. Phys. 15, 083044 (2015).
- G. Toscano, J. Straubel, A. Kwiatkowski, C. Rockstuhl, F. Evers, H. Xu, N. A. Mortensen, and M. Wubs, Nature Commun. 6, 7132 (2015).
- C. Ciracì and F. Della Sala, Phys. Rev. B 93, 205405 (2016).
- T. Christensen, W. Yan, A.-P. Jauho, M. Soljačić, and N. A. Mortensen, Phys. Rev. Lett. 118, 157402 (2017).
- K. J. Savage, M. M. Hawkeye, R. Esteban, A. G. Borisov, J. Aizpurua, and J. J. Baumberg, Nature 491, 574 (2012).
- J. A. Scholl, A. García-Etxarri, A. L. Koh, and J. A. Dionne, Nano Lett. 13, 564 (2013).
- W. Yan, M. Wubs, and N. A. Mortensen, Phys. Rev. Lett. 115, 137403 (2015).
- R. Filter, C. Bösel, G. Toscano, F. Lederer, and C. Rockstuhl, Opt. Lett. 39, 6118 (2014).
- S. Raza, M. Wubs, S. I. Bozhevolnyi, and N. A. Mortensen, Opt. Lett. 40, 839 (2015).
- C. Tserkezis, N. Stefanou, M. Wubs, and N. A. Mortensen, = Nanoscale 8, 17532 (2016).
- C. Tserkezis, N. A. Mortensen, and M. Wubs, arXiv 1703.0072 (2017).
- D. C. Marinica, H. Lourenço-Martins, J. Aizpurua, and A. G. Borisov, Nano Lett. 13, 5972 (2013).
- C. Ciracì, R. T. Hill, J .J. Mock, Y. Urzhumov, A. I. Fernández-Domínguez, S. A. Maier, J. B. Pendry, A. Chilkoti, and D. R. Smith, Science 337, 1072 (2012).
- H. Shen, L. Chen, L. Ferrari, M.-H. Lin, N. A. Mortensen, S. Gwo, and Z. Liu, Nano Lett. 17, 2234 (2017).
- N. A. Mortensen, S. Raza, M. Wubs, T. Søndergaard, and S. I. Bozhevolnyi, Nature Commun. 5, 3809 (2014).
- M. Wubs and N. A. Mortensen, Nonlocal response in plasmonic nanostructures in Quantum Plasmonics, Eds. S. I. Bozhevolnyi, L. Martín-Moreno, and F. J. García-Vidal, (Springer, Cham, 2017).
- X. Zhang, L. Guo, J. Luo, X. Zhao, T. Wang, Y. Li, and Y. Fu, ACS Appl. Mater. Interfaces 8, 9889 (2016).
- C. Tserkezis, G. Gantzounis, and N. Stefanou, J. Phys.: Condens. Matter 20, 075232 (2008).
- S. Raza, G. Toscano, A.-P. Jauho, N. A. Mortensen, and M. Wubs, Plasmonics 8, 193 (2013).
- S. J. Oldenburg, R. D. Averitt, S. L. Westcott, and N. J. Halas, Chem. Phys. Lett. 288, 243 (1998).
- E. Prodan, C. Radloff, N. J. Halas, and P. Nordlander, Science 302, 419 (2003).
- S. Rudin, and T. L. Reinecke, Phys. Rev. B 59, 10227 (1999).
- M. R. Philpott, A. Brillante, I. R. Pockrand, and J. D. Swalen, Mol. Cryst. Liq. Cryst. 50, 139 (1979).
- M. J. Gentile and W. L. Barnes, J. Opt. 19, 035003 (2017).
- Absorption and Scattering of Light by Small Particles, C. F. Bohren and D. R. Huffman (Wiley, New York, 1983).
- N. Stefanou, C. Tserkezis, and G. Gantzounis, Proc. SPIE 6989, 698910 (2008).
- Mathematical Methods for Physicist, Fifth Edition, G. B. Arfken and H. J. Weber (Academic Press, San Diego, 2000).
- J. B. Khurgin and G. Sun, J. Opt. Soc. Am. B 26, B83 (2009).
- S. A. Maier, Opt. Express 14, 1957 (2006).
- G. Toscano, S. Raza, W. Yan, C. Jeppesen, S. Xiao, M. Wubs, A.-P. Jauho, S. I. Bozhevolnyi, N. A. Mortensen, Nanophotonics 2, 161 (2013).
- P. B. Johnson and R. W. Christy, Phys. Rev. B 6, 4370 (1972).
- C. Tserkezis, J. R. Maack, Z. Liu, M. Wubs, and N. A. Mortensen, Sci. Rep. 6, 28441 (2016).
- R. Ruppin, Phys. Rev. Lett. 31, 1434 (1973).
- Comsol Multiphysics, www.comsol.com
- Nanoplasmonics Lab, www.nanopl.org
- U. Kreibig and L. Genzel, Surf. Sci. 156 678 (1985).
- T. V. Teperik, P. Nordlander, J. Aizpurua, and A. G. Borisov, Phys. Rev. Lett. 110, 263901 (2013).
- C. Ciracì, Phys. Rev. B 95, 245434 (2017).