|
|
1.INTRODUCTIONStrong coupling of localised surface plasmons in metal-dielectric structures with excitons in dye molecules or semiconductors attracts intense interest due to emerging applications such as ultrafast reversible switching,1–3 quantum computing,4, 5 and light harvesting.6 An extremely strong field confinement on the length scale well below the diffraction limit provided by plasmonic structures gives rise to exciton-plasmon coupling that is much stronger than the coupling to cavity modes in semiconductor microcavites. Strong coupling regime between two systems is established when the energy exchange between them takes place faster than the energy dissipation, which leads to an anticrossing gap (Rabi splitting) in the dispersion of mixed states.7 While relatively weak Rabi splittings ~ 1 meV were reported for semiconductor quantum dots (QDs) coupled to a cavity mode,8–10 much greater splittings (up to 500 meV) were observed for surface plasmons coupled to excitons in J-aggregates,11–20 various dye molecules,21–25 semiconductors QDs,26–28 or two-dimensional atomic crystalls.32–35 In the classical picture of coupled oscillators,7 the interaction between a quantum emitter (QE), modeled here by a two-level system with frequency ω0, and a cavity (or plasmon) mode with frequency wm, gives rise to mixed states that show up in optical spectra through splitting of scattering or emission peaks into two polaritonic bands separated (for ωm = ω0) by , where γ0 and γm are, respectively, the QE’s and mode’s decay rates, and g is the coupling parameter,9 Here μ is the QE dipole matrix element and is the mode volume characterizing field confinement at the QE position. Since typically , the onset of transition to the strong coupling regime is , which demands sufficiently strong field confinement (or small mode volume. For a QD placed inside a microcavity, the cavity mode volume is at least , where λ is the wavelength, resulting in Δ ~ 1 meV.9 However, if a QE is located at a “hot spot” near a plasmonic structure characterized by extremely strong field confinement, e.g., near a sharp metal tip or in a gap between two metal structure, the Rabi splitting can reach much greater values even though the plasmon decay rate is much higher than that for cavity modes.36–39 If an ensemble of N QEs is coupled to a cavity mode, they form a collective state comprised of all QEs oscillating in sync with each other.40 For such a collective state, the coupling scales as with the ensemble size, while the rest of collective states are “dark”, i.e., they interact only weakly with the electromagnetic field. Importantly, such scaling implies that individual QEs are coupled to the cavity mode with approximately equal strength, so that the mode, in fact, interacts with one “giant oscillator” whose amplitude is enhanced by the factor N relative to individual QE. This picture is affected dramatically if an ensemble of N QEs is placed near a plasmonic nanostructure, which is normally characterized by strongly varying local fields. For example, in a typical experimental setup, the QEs (excitons in J-aggregates or QDs) are embedded in a dielectric shell enclosing a metal core supporting localized surface plasmon (see Fig. 1). Strong exciton-plasmon coupling requires large plasmon local density of states (LDOS) facilitating an efficient energy transfer (ET) between a plasmon and a QE placed in a strong local field region. In plasmonic structures, the local fields can be very strong close to metal surface, especially near sharp features such as narrow tips or surface irregularities, but fall off rapidly away from the surface. In this case, the coupling between individual QEs and resonant plasmon mode can vary in a wide range, so that the classical “giant oscillator” picture is no longer useful and, instead, one has to resort to the underlying mechanism of energy exchange between the plasmon and QEs. Furthermore, since for sufficiently remote QEs, the individual QE-plasmon ET rates are small, it is evident that the ensemble QE-plasmon coupling should saturate as the region QEs are distributed in expands. Saturation of Rabi splitting with an increasing system size was recently reported for molecular excitons in J-aggregates embedded in dielectric shell enclosing a gold nanoprism.20 In this Letter, we develop a model for exciton-pasmon coupling based upon microscopic picture of energy exchange between system components. We establish an explicit relation between the ensemble QE-plasmon coupling and the rate of cooperative energy transfer (CET) between a collective state and resonant plasmon mode (see Fig. 1), and estimate the energy-exchange frequency in the strong coupling regime. By defining carefully the plasmon mode volume in terms of plasmon LDOS, we show that, for a single QE at a hot spot near sharp metal tip, the QE-plasmon coupling scales as where Vmet is the metal volume fraction that confines the plasmon field, and provide an analytical expression for the enhancement factor relative to exciton coupling to a cavity mode. For an ensemble of QEs near a plasmonic structure, we show that the ensemble coupling parameter scales as , where is the average plasmon mode volume in the QE region. If QEs are placed in a region with nearly constant, however large, plasmon LDOS, the coupling exhibits the usual cavity-like scaling . However, in open plasmonic systems, if QEs are uniformly distributed, with concentration n, in an extended region saturating the plasmon mode volume (see Fig. 1), the ensemble QE-plasmon coupling saturates to some value that, for system size below the diffraction limit, has a universal form where is the metal dielectric function and is the dielectric constant in the QE region. The saturated coupling gs is independent of system’s size and shape, except indirectly via the plasmon frequency wm. Surprisingly, gs is comparable with the coupling between QEs, uniformly distributed inside a microcavity with the same concentration, and a cavity mode, suggesting that large Rabi splittings, observed in open plasmonic systems, are likely due to high QE concentrations. 2.PLASMON LDOS, MODE VOLUME AND ENERGY TRANSFER RATETo establish a relation between cavity-like coupling (1) and microscopic picture of energy exchange in plasmonic systems, we employ a classical approach for plasmons and describe excitons by two-level systems, reffered to as QEs. We assume that QEs are placed near a metal-dielectric structure with characteristic size smaller than the radiation wavelength described by a complex dielectric function . The interaction of a QE, located at a position r with dipole orientation n, with eloctromagnetic environmentby is characterized by projected LDOS7 where is the rescaled electromagnetic dyadic Green function in the presence of plasmonic structure ( is the standard Green dyadic for Maxwell equations and k is the wave vector). In the near field limit, the Green function splits into free-space and plasmon contributions, , where the plasmon Green function, for frequency ω near the plasmon mode frequency ωm, has the form:41,42 Here, E(r) is the plasmon field, which we chose to be real, satisfying the Gauss’s law, , and tensor product of fields is implied. In Eq. (4), Um is the plasmon mode energy defined in the standard way as 43 while the factor ωm/4Um, is the plasmon pole residue in the complex frequency plane.42 Using Eqs. (3) and (4), we obtain the projected plasmon LDOS as a Lorentzian centered at the plasmon frequency, where Q = ωm/γm, is the plasmon quality factor. The plasmon LDOS describes the distribution of plasmon states in unit volume and frequency interval. Accordingly, its frequency integral, , defines the projected plasmon mode density as,41,42 where is the projected plasmon mode volume that characterizes the plasmon field confinement at a point r in the direction n. We stress that, for characteristic system size smaller than the radiation wavelength, the plasmon mode volume is a real function of r defined in terms of real plasmon fields at that position.42 For larger systems, where plasmon modes are strongly hybridized with the radiation continuum, field confinement is not a well-defined notion but, in this case, the QE-plasmon interaction can still be treated accurately using the expansion over quasi-normal modes characterized by complex mode volumes.44–48 Within our analytical approach, the above connection between the plasmon mode volume and plasmon LDOS allows relating the QE-plasmon coupling with the rate of energy transfer (ET) between them. In order to relate to the QE-plasmon ET rate , we note that the latter can be expressed via the plasmon Green function as . Using Eq. (4), we obtain the frequency-dependent ET rate as a Lorentzian centered at the plasmon frequency, At resonance frequency the QE-plasmon ET rate takes the form where is given by Eq. (7). Normalizing by the free-space radiative decay rate , we obtain the Purcell factor49 for a QE coupled to a plasmonic resonator: . To illustrate ET rate’s sensitivity to the QE position and system geometry, in Fig. 2 we show near the tip of Au nanorod in water. Nanorod was modeled by a prolate speroid with semimajor and semiminor axes a and b, respectively, and the experimental Au dielectric function was used in all calculations. To highlight the dependence of on system geometry, we normalize it by its value for a sphere with radius a and plot the result against the normalized distance d/a for several aspect ratios a/b. With increasing aspect ratio, i.e., reducing the nanorod volume relative to the sphere, the normalized ET rates increase dramatically near the nanorod tip, which translates to analogous behavior of the QE-plasmon coupling, as we discuss in detail later in the paper. Let us now turn to an ensemble of N QEs with dipole moments pi = μηi situated at positions ri near a plasmonic structure. Each QE interacts, via the coupling , with the common electric field generated by all QEs, where ) is the near-field electromagnetic Green’s function in the presence of plasmonic structure. Due to electromagnetic correlations between QEs, collective states are formed that are described by the eigenstates of the Green’s function matrix at QEs’ positions projected onto QEs’ dipole moments:50 . Note that this coupling matrix includes direct dipole-dipole interactions between the QEs which are causing random shifts of the QEs’ energies. However, for random dipole orientations, the direct dipole coupling vanishes on average,51,52 while its fluctuations contributes, among other factors, to inhomogeneous broadenning of the QE energies (we will return to this point later). Furthermore, for system size smaller than the radiation wavelength, the radiative (superradiant) coupling between the QEs is relatively weak as well and is not considered here. Therefore, near the plasmon resonance, the coupling matrix is dominated by the plasmon Green function Eq. (4), and we obtain53 The diagonal elements of plasmon coupling matrix Eq. (10) are complex, , and their real and imaginary parts describe, respectively, the frequency shifts δω(i) and decay rates of individual QEs due to interaction with a plasmon mode; the former contribute to inhomogeneous broadening of the QE energies, while the latter are QE-plasmon ET rates, given by Eq. (8). Note that at resonance, we have δω(i) = 0, while are given by Eq. (9). The off-diagonal elements of plasmon coupling matrix Eq. (10) describe plasmonic correlations between QEs which give rise to the collective states. These states are represented by the eigenstates ψ of the coupling matrix (10) satisfying , where the complex eigenvalues λ characterize the energy shifts and decay rates of the collective states. It is now easy to see that the vector is an eigenstate of the matrix Eq. (10), and the corresponding eigenvalue has the form , where is frequency shift of the collective state and is its decay rate. At resonance ω = ωη, we again have δωΝ = 0, while the decay rate of the collective state ψΝ takes the form [compare to Eq. (9)] where is the projected plasmon mode volume (7) at the ith QE position. Thus, the collective state described by ψΝ exchanges its energy with the plasmon mode cooperatively, i.e., at a rate equal to the sum of individual QE-plasmon ET rates.41,53 As we show below, the CET rate Eq. (11) determines onset of the transition to strong coupling regime for QE ensembles coupled to plasmonic resonators. Since the imaginary part of the eigenvalue λΝ saturates the ET rates from the QEs to resonant plasmon mode, the rest of the collective states are not coupled to that mode but, in principle, can still couple to the off-resonant modes and radiation field. Note, however, that for large ensembles, coupling of collective states to the off-resonant modes is relatively weak,54,55 while large plasmonic Purcell factors ensure that coupling to the radiation field is relatively weak as well, implying that the exciton-plasmon energy exchange at the CET rate (11) is the dominant energy flow channel in the system. We stress that, in plasmonic systems, the collective states are formed due to QEs’ correlations via the local fields that can vary strongly near a plasmonic structure. Therefore, these states are distinct from the superradiant and subradiant states, which emerge due to QEs’ coupling to the common radiation field that is nearly uniform on the system scale. 3.EXCITON-PLASMON COUPLING AND ENERGY EXCHANGEConsider first a single QE situated near a plasmonic resonator with its frequency ω0 close to the plasmon frequency ωm. We consider systems with characteristic size smaller the radiation wavelength and, therefore, consider only the near-field coupling between a QE and resonant plasmon mode; note, however, that for larger systems, radiative (superradiant) coupling between them can significatly affect the optical spectra in the strong coupling regime.18 Using Eq. (8), we can now relate the QE-plasmon coupling (1) to the QE-plasmon ET rate as The above expression for g implies a relation between the QE-plasmon coupling and the Purcell factor: g2 = . To describe the transition to strong coupling regime, we note that, within classical model of coupled oscillators, the spectrum splits into upper and lower polaritonic bands with complex frequencies7 , where and . To simplify the analysis, we assume for now that QE and plasmon are in resonance, ω0 = ωm. Typically, in plasmonic systems, the intrinsic QE decay rate γ0 is much smaller than either γm or γet, and so we disregard it here as well; however small, γ0 does play a crucial role in the interference effects, such as Fano resonances,39 but such phenomena are beyond our work’s scope. In the weak coupling regime, both polaritonic bands are centered at the same frequency ωm but are characterized by different spectral widths, where we used the relation (12). For , the expansion of Eq. (13) over the small parameter yields and , which represent, respectively, the decay rates of weakly-mixed plasmon and QE states. In the strong coupling regime, both polariton bands are characterized by constant spectral width γm/2, while their central frequencies are separated by the Rabi splitting: Thus, in terms of QE-plasmon ET rate, the transition to strong coupling regime take place at Note that the polaritonic bands become spectrally distinct when Rabi splitting exceeds their linewidth γm/2, which corresponds to the condition . The onset condition (15) implies that, for plasmonic resonators with high quality factor , the QE-plasmon energy exchange takes place on a much longer time scale than the optical period. Qualitatively, the energy exchange dynamics can be described by the damped oscillator equation, , where x characterizes the QE-plasmon energy balance while frequency is the two-way QE-plasmon ET rate. After a pulsed excitation, the energy balance undergoes damped oscillations with energy-exchange frequency and damping rate γm/2, which is similar to that for optical field intensity in the strong coupling regime [see Eq. (14)]. In this model, the transition from strong to weak coupling regime is analogous to the transition from underdamped to overdamped regime taking place at , which coincides with Eq. (15). Let us now turn to an ensemble of N QEs near a plasmonic structure. As discussed in the previous section, plasmonic correlations between the QEs give rise to a collective state that exchanges its energy with resonant plasmon mode at the CET rate , given by (11). Accordingly, the ensemble QE-plasmon coupling gN is related to the CET rate as [compare to Eq. (12)] where is the average plasmon mode volume for QE ensemble, defined as The ensemble QE-plasmon coupling (17) is related as to the cooperative Purcell factor characterizing enhancement of the CET-based cooperative emission rate.53,56 For a QE ensemble, the transition onset to strong coupling regime is obtained with replacing the ET rate in Eq. (15) by the CET rate: . In the strong coupling regime, the energy balance between the collective state and resonant plasmon mode oscillates with energy-exchange frequency , which vanishes at the transition point, and the same damping rate γm/2 [compare to Eq. (16)]. The ensemble QE-plasmon coupling gN is related to individual QE couplings , and, therefore, it depends sensitively on the plasmon LDOS variations over the region QEs are distributed in. For example, if the plasmon LDOS (and, accordingly, the mode volume) is nearly constant in the QE region, e.g., within dielectric core enclosed by a metallic shell, the individual couplings gi are approximately equal and, hence, the ensemble QE-plasmon coupling exhibits a cavity-like scaling, . However, in open plasmonic systems with large number of QEs distributed outside the plasmonic structure (see Fig. 1), the ensemble QE-plasmon coupling saturates to a universal value, as we show later in this paper. In the above analysis, we assumed same excitation frequencies for all QEs in the ensemble. In the experiment, however, the QE frequencies are distributed within some interval due to, e.g., direct dipole-dipole interactions between the QEs or, in the case of semiconductor QDs, their size variations. Here we note that the ET rate between a donor and an acceptor is determined by the spectral overlap of their respective emission and absorption bands.7 Therefore, as long as the ensemble inhomogeneous broadening stays within the broad plasmon resonance band, the QE-plasmon energy exchange mechanism remains largely unaffected. 4.EXCITON-PLASMON COUPLING NEAR SHARP METAL TIPThe largest values of exciton-plasmon coupling are achieved for QEs placed in a region with large plasmon LDOS that provides an efficient ET between a QE and a plasmon mode. Here, we consider a single QE situated near a sharp tip of a small plasmonic structure, such as a metal nanorod, where the field confinement can be extremely strong (hot spot). To estimate the QE-plasmon coupling, we note that the Gauss’s law implies , so the coupling (12) has the form where we assumed that only in the metallic region the dielectric function is dispersive. The integral over metallic region in Eq. (19) depends on the characteristic size of that region lm relative to the skin penetration length ls (about 20 nm for Au in plasmon frequency domain). For small structures with lm < ls, the QE-plasmon coupling scales as , while in opposite case lm < ls, the metal volume Vmet should be replaced by the effective volume the plasmon field is largely confined to. In either case, there is a significant enhancement of g relative to the QE coupling to a cavity mode, , so in the analysis below we use the notation Vmet throughout. Importantly, in addition to this geometric volume effect, there is also strong field enhancement due to “lightning rod” effect near sharp metal tips or surface irregularities. To elucidate the relative importance of these two enhancement sources, below we estimate the QE’s coupling g to a plasmon mode oscillating along the tip of a small metal structure (see Fig. 3). For small systems, the plasmon field does not significantly change inside the metallic structure while falling off rapidly outside of it, so the largest field enhancement takes place near the tip. For QE polarized along the tip (see Fig. 3), the coupling (19) at the system symmetry axis (z-axis) can be estimated as where subscripts in and out refer, respectively, to the fields at the interface on metal and dielectric sides. Using the boundary condition , we obtain where is the plasmon field outside the structure normalized by its value at the tip. For good Drude metals in the plasmonic frequency range, we have , and the coupling (21) takes the form The maximal value of gtip is achieved in a close proximity to the tip, where . Comparing this value to the coupling gcav for a QE inside a cavity, given by Eq. (1) with , we obtain the enhancement factor, indicating that the plasmonic enhancement is due to both the geometric volume effect and the plasmon field enhancement characterized by . Note that for metal nanostructures with characteristic size larger than the skin penetration length 1s, the metal volume in Eq. (23) should be replaced by , as discussed above. In Fig. 3, we plot distance dependence of the QE-plasmon coupling for a QE near the tip of Au nanorod in water. Note that for a prolate spheroid used here to model the nanorod, the expression (21) is exact. The curves show calculated , normalized by , plotted against normalized distance d/a for different values of aspect ratio a/b at fixed nanorod length a, i.e., for different metal volumes Vmet = ab2. It is clearly seen that for a/b > 1, i.e., when the system possesses hot spots near nanorod tips, the QE-plasmon coupling scales as with its overall magnitude falling off rapidly away from the tip. Note that while the coupling is largest at the tip , it is expected to saturate below QE-tip distances as the nonlocal effects become important.57–59 In noble metals, this length scale is < 1 nm in the plasmonic frequency range and, therefore, the distance dependence in Fig. 3 is cut off below d/a = 0.05. 5.SATURATION OF EXCITON-PLASMON COUPLING IN OPEN SYSTEMSIn open plasmonic systems, the QEs are distributed within some region outside the metal structure, e.g., within dielectric shell, with dielectric constant εd, enclosing a metallic core [see Fig. (1)]. In such systems, the plasmon field E(r) can vary substantially within the QE region and, in particular, falls off rapidly away from the metal surface. If QEs are uniformly distributed, with concentration n, within some volume V0, the sum in Eq. (18) is replaced by the integral over V0, and the corresponding QE-plasmon coupling g0 takes the form where the factor 1/3 comes from orientational averaging. In contrast to individual QE-plasmon coupling (19), the ensemble coupling (24) is determined by the ratio of integrated field intensities over QE and metallic regions. If the QE region is sufficiently extended, the remote QEs do not interact with the plasmon mode and, therefore, the integral in Eq. (24) is independent of V0. In this case, the QEs saturate the entire plasmon mode volume, i.e., beyond the QE region, the plasmon LDOS is negligibly small, so that, with expanding QE region, the ensemble QE-plasmon coupling (24) saturates to some value gs. Note that, since the coupling (24) is independent of losses, such saturation is due to reduction of QE-plasmon ET rates away from the plasmonic structure, rather than any dephasing processes. Let us now show that the actual value of saturated coupling gs does not depend directly on the size or shape of the metal structure as well and, remarkably, can be obtained explicitly for any plasmonic system with characteristic size below the diffraction limit. Indeed, due to the Gauss’s law, both volume integrals in Eq. (24) reduce to surface integrals over the system interfaces, including the common metal-dielectric interface. After matching the normal field components at this interface and disregarding, in the saturated case, the outer boundary of the QE region, the ratio of integrated intensities in Eq. (24) is found as , and we arrive at a universal saturated QE-plasmon coupling, Using the relation for good Drude metals, the coupling (25) simplifies to and the transition onset to strong coupling regime for large ensembles, , can be presented as Note that saturation of the QE-plasmon coupling in the strong coupling regime was recently reported for photolumenescence of molecular excitons in J-aggregates embedded in dielectric shell enclosing an Au nanoprism.20 In Fig. 4, we show the calculated QE-plasmon coupling g0, given by Eq. (24), for QEs distributed uniformly within SiO2 shell enclosing an Au nanorod. This core-shell system is modeled by two confocal spheroids with major semi-axes a and a1 corresponding to the Au/SiO2 and Si2O/H2O interfaces, respectively. With expanding QE region, the coupling g0 saturates to the value gs given by Eq. (25). Saturation is faster for more elongated particles possessing hot spots near their tips, in which case the QE-plasmon coupling reaches 90% of its saturated value for a1/a = 2. Finally, let us compare the coupling parameters for large QE ensembles coupled to plasmons and to cavity modes in the case when QEs saturate their respective mode volumes. While in plasmonic cavities, the QE-plasmon coupling is strongly enhanced relative to QEs’ coupling to a cavity mode, this is not the case for open plasmonic systems, where only a small fraction of QEs close to metal surface interacts strongly with the plasmon mode. Indeed, let us assume that QEs distributed uniformly, with concentration n, inside a microcavity. The QE-cavity coupling is also given by Eq. (17), but with cavity mode volume Vcav given by9 where is the cavity dielectric function. After averaging Eq. (28) over the cavity volume Vcav, we obtain , and the saturated coupling parameter for microcavities takes the form Comparing to Eq. (26), we conclude that, for same QE concentrations, the saturated coupling for open plasmonic system is of the same order as for microcavities, implying that large Rabi splitting observed in open plasmonic systems are due to high QE concentrations. Note also that, due to much higher microcavity quality factors, the transition onset is lower as well. 6.CONCLUSIONSIn conclusion, we developed a model for exciton-plasmon coupling based on a microscopic picture of energy exchange that drives the transition to strong coupling regime. Plasmonic correlations between the QEs give rise to a collective state that transfers its energy cooperatively to a resonant plasmon mode at a rate equal to the sum of individual QE-plasmon ET rates. It is this CET rate that, along with plasmon decay rate, determines the QE-plasmon coupling as well as the energy-exchange frequency in the strong coupling regime. By defining accurately the plasmon mode volume, we have shown that for a QE near a sharp metal tip, the QE-plasmon coupling is enhanced by the factor relative to the QE coupling to a cavity mode, where Vmet is metal volume that confines the plasmon field. For an ensemble of N QEs placed in a region with weakly varying plasmon LDOS, the QE-plasmon coupling exhibits cavity-like scaling . However, in open plasmonic systems, the ensemble QE-plasmon coupling saturates to a universal value gs that does not depend on plasmonic system’s size and shape except indirectly via the plasmon frequency. Finally, we compared coupling parameters for QEs coupled plasmon and cavity modes to find that, for same QE concentrations, the open plasmonic systems offer no significant enhancement, implying that the observed large Rabi splittings are likely due to high QE concentrations in such systems. ACKNOWLEDGMENTSThis work was supported in part by the National Science Foundation grants DMR-2000170, DMR-1856515, DMR-1826886 and HRD-1547754. REFERENCEST. Schwartz, J. A. Hutchison, C. Genet, and T. W. Ebbesen, Phys. Rev. Lett., 106 196405
(2011). https://doi.org/10.1103/PhysRevLett.106.196405 Google Scholar
A.-L. Baudrion, A. Perron, A. Veltri, A. Bouhelier, P.-M. Adam, and R. Bachelot, Reversible Strong Coupling in Silver Nanoparticle Arrays Using Photochromie Molecules, Nano Lett., 13 282
–286
(2013). Google Scholar
L. Lin, M. Wang, X. Wei, X. Peng, C. Xie, and Y. Zheng, Nano Lett., 16 7655
–7663
(2016). https://doi.org/10.1021/acs.nanolett.6b03702 Google Scholar
S. Sun, H. Kim, G. S. Solomon, and E. Waks, Nat. Nanotechnol., 11 539
–544
(2016). https://doi.org/10.1038/nnano.2015.334 Google Scholar
L. De Santis, C. Anton, B. Reznychenko, N. Somaschi, G. Coppola, J. Senellart, C. Gomez, A. Lemaitre, I. Sagnes, A. G. White, L. Lanco, A. Auffeves, and P. Senellart, Nat. Nanotech, 12 663
–667
(2017). https://doi.org/10.1038/nnano.2017.85 Google Scholar
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
–6856
(2016). https://doi.org/10.1021/acs.nanolett.6b02661 Google Scholar
L. Novotny and B. Hecht, Principles of Nano-Optics, CUP, New York
(2012). https://doi.org/10.1017/CBO9780511794193 Google Scholar
J. P. Reithmaier, G. Sek, A. Löffler, C. Hofmann, S. Kuhn, S. Reitzenstein, L. V. Keldysh, V. D. Kulakovskii, T. L. Reinecke, and A. Forchel, Nature, 432 197
(2004). https://doi.org/10.1038/nature02969 Google Scholar
G. Khitrova, H. M. Gibbs, M. Kira, S. W. Koch, and A. Scherer, Nature Phys., 2 81
(2006). https://doi.org/10.1038/nphys227 Google Scholar
K. Hennessy, A. Badolato, M. Winger, D. Gerace, M. Atatüre, S. Gulde, S. Falt, E. L. Hu, and A. Imamoglu, Nature, 445 896
(2006). https://doi.org/10.1038/nature05586 Google Scholar
J. Bellessa, C. Bonnand, J. C. Plenet, and J. Mugnier, Phys. Rev. Lett., 93 036404
(2004). https://doi.org/10.1103/PhysRevLett.93.036404 Google Scholar
Y. Sugawara, T. A. Kelf, J. J. Baumberg, M. E. Abdelsalam, and P. N. Bartlett, Phys. Rev. Lett., 97 266808
(2006). https://doi.org/10.1103/PhysRevLett.97.266808 Google Scholar
G. A. Wurtz, P. R. Evans, W. Hendren, R. Atkinson, W. Dickson, R. J. Pollard, A. V. Zayats, W. Harrison, and C. Bower, Nano Lett., 7 1297
(2007). https://doi.org/10.1021/nl070284m Google Scholar
N. T. Fofang, T.-H. Park, O. Neumann, N. A. Mirin, P. Nordlander, and N. J. Halas, Nano Lett., 8 3481
(2008). https://doi.org/10.1021/nl8024278 Google Scholar
P. Vasa, R. Pomraenke, S. Schwieger, Yu. I. Mazur, Vas. Kunets, P. Srinivasan, E. Johnson, J. E. Kihm, D. S. Kim, E. Runge, G. Salamo, and C. Lienau, Phys. Rev. Lett., 101 116801
(2008). https://doi.org/10.1103/PhysRevLett.101.116801 Google Scholar
J. Bellessa, C. Symonds, K. Vynck, A. Lemaitre, A. Brioude, L. Beaur, J. C. Plenet, P. Viste, D. Felbacq, E. Cambril, and P. Valvin, Phys. Rev. B, 80 033303
(2009). https://doi.org/10.1103/PhysRevB.80.033303 Google Scholar
A. E. Schlather, N. Large, A. S. Urban, P. Nordlander, and N. J. Halas, Nano Lett., 13 3281
(2013). https://doi.org/10.1021/nl4014887 Google Scholar
W. Wang, P. Vasa, R. Pomraenke, R. Vogelgesang, A. De Sio, E. Sommer, M. Maiuri, C. Manzoni, G. Cerullo, and C. Lienau, ACS Nano, 8 1056
(2014). https://doi.org/10.1021/nn405981k Google Scholar
G. Zengin, M. Wersäll, S. Nilsson, T. J. Antosiewicz, M. Käll, T. Shegai, Phys. Rev. Lett., 114 157401
(2015). https://doi.org/10.1103/PhysRevLett.114.157401 Google Scholar
M. Wersäll, J. Cuadra, T. J. Antosiewicz, S. Balci, and T. Shegai, Nano Lett., 17 551
–558
(2017). https://doi.org/10.1021/acs.nanolett.6b04659 Google Scholar
T. K. Hakala, J. J. Toppari, A. Kuzyk, M. Pettersson, H. Tikkanen, H. Kunttu, and P. Torma, Phys. Rev. Lett., 103 053602
(2009). https://doi.org/10.1103/PhysRevLett.103.053602 Google Scholar
A. Berrier, R. Cools, C. Arnold, P. Offermans, M. Crego-Calama, S. H. Brongersma, and J. Gomez-Rivas, ACS Nano, 5 6226
(2011). https://doi.org/10.1021/nn201077r Google Scholar
A. Salomon, R. J. Gordon, Y. Prior, T. Seideman, and M. Sukharev, Phys. Rev. Lett., 109 073002
(2012). https://doi.org/10.1103/PhysRevLett.109.073002 Google Scholar
A. De Luca, R. Dhama, A. R. Rashed, C. Coutant, S. Ravaine, P. Barois, M. Infusino, and G. Strangi, Appl. Phys. Lett., 104 103103
(2014). https://doi.org/10.1063/1.4868105 Google Scholar
V. N. Peters, T. U. Tumkur, Jing Ma, N. A. Kotov, and M. A. Noginov, Opt. Express, 24 25653
(2016). https://doi.org/10.1364/OE.24.025653 Google Scholar
D. E. Gomez, K. C. Vernon, P. Mulvaney, and T. J. Davis, Nano Lett., 10 274
(2010). https://doi.org/10.1021/nl903455z Google Scholar
D. E. Gomez, S. S. Lo, T. J. Davis, and G. V. Hartland, J. Phys. Chem. B, 117 4340
(2013). https://doi.org/10.1021/jp306830s Google Scholar
A. Manjavacas, F. J. Garcia de Abajo, and P. Nordlander, Nano Lett., 11 2318
(2011). https://doi.org/10.1021/nl200579f Google Scholar
P. Vasa, R. Pomraenke, S. Schwieger, Y. I. Mazur, V. Kunets, P. Srinivasan, E. Johnson, J. E. Kihm, D. S. Kim, E. Runge, G. Salamo, and C. Lienau, Phys. Rev. Lett., 101 116801
(2008). https://doi.org/10.1103/PhysRevLett.101.116801 Google Scholar
B. J. Lawrie, K.-W. Kim, D. P. Norton, and R. F. Haglund Jr., Nano Lett., 12 6152
(2012). https://doi.org/10.1021/nl3029784 Google Scholar
D. G. Baranov, M. Wersall, J. Cuadra, T. J. Antosiewicz, and T. Shegai, ACS Phot., 5 24
–42
(2018). https://doi.org/10.1021/acsphotonics.7b00674 Google Scholar
J. Wen, H. Wang, W. Wang, Z. Deng, C. Zhuang, Y. Zhang, F. Liu, J. She, J. Chen, H. Chen, S. Deng, and N. Xu, Nano Lett., 17 4689
–4697
(2017). https://doi.org/10.1021/acs.nanolett.7b01344 Google Scholar
A. Krasnok, S. Lepeshov, and A. Alu, Optics Express, 26 15972
(2018). https://doi.org/10.1364/OE.26.015972 Google Scholar
C. Schneider, M. M. Glazov, T. Korn, S. Häfling, and B. Urbaszek, Nat. Comm., 9 2695
(2018). https://doi.org/10.1038/s41467-018-04866-6 Google Scholar
M. Stährenberg, B. Munkhbat, D. G. Baranov, J. Cuadra, A. B. Yankovich, T. J. Antosiewicz, E. Olsson, and T. Shegai, Nano Lett., 18 5938
–5945
(2018). https://doi.org/10.1021/acs.nanolett.8b02652 Google Scholar
K. Santhosh, O. Bitton, L. Chuntonov, and G. Haran, Nat. Commun., 7 11823
(2016). https://doi.org/10.1038/ncomms11823 Google Scholar
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
–130
(2016). https://doi.org/10.1038/nature17974 Google Scholar
P. Vasa and C. Lienau, ACS Phot., 5 2
–23
(2018). https://doi.org/10.1021/acsphotonics.7b00650 Google Scholar
H. Leng, B. Szychowski, M.-C. Daniel, and M. Pelton, Nat. Comm., 9 4012
(2018). https://doi.org/10.1038/s41467-018-06450-4 Google Scholar
G. Khitrova, H. M. Gibbs, F. Jahnke, M. Kira, and S. W. Koch, Rev. Mod. Phys., 71 1591
–1639
(1999). https://doi.org/10.1103/RevModPhys.71.1591 Google Scholar
T. V. Shahbazyan, Phys. Rev. Lett., 117 207401
(2016). https://doi.org/10.1103/PhysRevLett.117.207401 Google Scholar
T. V. Shahbazyan, Phys. Rev. B, 98 115401
(2018). https://doi.org/10.1103/PhysRevB.98.115401 Google Scholar
L. D. Landau and E. M. Lifshitz, Electrodynamics of Continuous Media, Elsevier, Amsterdam
(2004). Google Scholar
C. Sauvan, J. P. Hugonin, I. S. Maksymov, and P. Lalanne, Phys. Rev. Lett., 110 237401
(2013). https://doi.org/10.1103/PhysRevLett.110.237401 Google Scholar
R.-C. Ge, P. T. Kristensen, J. F. Young, and S. Hughes, New J. Phys., 16 113048
(2014). https://doi.org/10.1088/1367-2630/16/11/113048 Google Scholar
P. T. Kristensen and S. Hughes, ACS Phot., 1 2
(2014). https://doi.org/10.1021/ph400114e Google Scholar
E. A. Muljarov and W. Langbein, Phys. Rev. B, 94 235438
(2016). https://doi.org/10.1103/PhysRevB.94.235438 Google Scholar
P. Lalanne, W. Yan, K. Vynck, C. Sauvan, and J.-P. Hugonin, Laser Photon. Rev., 12 1700113
(2018). https://doi.org/10.1002/lpor.v12.5 Google Scholar
E. M. Purcell, Phys. Rev., 69 681
(1946). https://doi.org/10.1103/PhysRev.69.37 Google Scholar
M. Gross and S. Haroche, Phys. Rep., 93 301
(1982). https://doi.org/10.1016/0370-1573(82)90102-8 Google Scholar
R. Friedberg, S. R. Hartmann, and J. T. Manassah, Phys. Rep. C, 7 101
(1973). https://doi.org/10.1016/0370-1573(73)90001-X Google Scholar
T. V. Shahbazyan, M. E. Raikh, and Z. V. Vardeny, Phys. Rev. B, 61 13266
(2000). https://doi.org/10.1103/PhysRevB.61.13266 Google Scholar
T. V. Shahbazyan, Phys. Rev. B, 99 125143
(2019). https://doi.org/10.1103/PhysRevB.99.125143 Google Scholar
A. Delga, J. Feist, J. Bravo-Abad, and F. J. Garcia-Vidal, Phys. Rev. Lett., 112 253601
(2014). https://doi.org/10.1103/PhysRevLett.112.253601 Google Scholar
L. S. Petrosyan and T. V. Shahbazyan, Phys. Rev. B, 96 075423
(2017). https://doi.org/10.1103/PhysRevB.96.075423 Google Scholar
M. ElKabbash, E. Miele, A. K. Fumani, M. S. Wolf, A. Bozzola, E. Haber, T. V. Shahbazyan, J. Berezovsky, F. De Angelis, and G. Strangi, Google Scholar
C. Ciraci, 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). https://doi.org/10.1126/science.1224823 Google Scholar
N. A. Mortensen, S. Raza, M. Wubs, T. Søndergaard, and S. I. Bozhevolnyi, Nat. Commun., 5 3809
(2014). https://doi.org/10.1038/ncomms4809 Google Scholar
C. Tserkezis, M. Wubs, N. A. Mortensen, ACS Phot., 5 133
(2018). https://doi.org/10.1021/acsphotonics.7b00538 Google Scholar
|