Relaxation Dynamics of a Quantum Emitter Resonantly Coupled to a Coherent State of a Localized Surface Plasmon

We investigate the relaxation dynamics of a quantum dipole emitter (QDE), e.g., a molecule or quantum dot, located near a metal nanoparticle (MNP) exhibiting a dipolar localized surface plasmon (LSP) resonance at the frequency of the QDE radiative transition. A generic three-level QDE, which is pumped with an external laser pulse and thereby brought into an optically active excited state, is considered to be weakly coupled to the resonant LSP described by a coherent state. It is shown that, under the condition of the QDE–MNP characteristic relaxation time being much shorter than that of the QDE in free space but much longer than the LSP lifetime, the QDE relaxation dynamics can be described analytically and feature, in general, non-exponential decay with complicated transient behaviour. The main physical consequence of this relaxation process is that the emission, being largely determined by the MNP, comes out with a substantial delay. It is also shown that energy dissipation in the QDE–MNP system is relatively weak with the probability of the photon emission being $0.75, a number which, rather surprisingly, does not explicitly depend on the metal absorption characteristics. A large number of QDE–MNP system parameters in our analytical description open new possibilities for controlling quantum emitter dynamics.


Introduction
The interaction of quantum dipole emitters (QDEs), such as molecules or quantum dots, with metal nanoparticles (MNPs) at optical frequencies allows control over the ow of electromagnetic energy and lies at the core of an explosively growing eld of quantum plasmonics. 1 The coupling of electromagnetic elds in dielectrics and free electron oscillations in metals, which results in surface plasmon excitations, enables the localization of light in subwavelengthsized volumes and, in general, command over light-matter interactions at the nanoscale.One of the very fundamental effects occurring as a result of these interactions is the modication of QDE relaxation dynamics due to the presence of metal nanostructures, e.g.MNPs, strongly perturbing local electromagnetic elds.Recent advances in nano-optics, especially experiments with single molecules interacting with well-dened metal nanostructures, [2][3][4] oen referred to as nanoantennas, serve as a strong impetus for further developments in this direction. 5,68][9] It is also expected the QDE-MNP interaction can enter the regime of strong coupling, where excitation energy is coherently transferred between the QDE and MNP in the form of Rabi oscillations. 1,10he strong coupling between a QDE and a resonator is a key element of the quantum-optics toolbox.In the absence of energy dissipation, two identical oscillators are strongly coupled even for small coupling constants.In practice, their coupling should be sufficiently strong so that the energy level splitting becomes larger than the width of the oscillator resonances or, alternatively, that the period of Rabi oscillations becomes smaller than the oscillator lifetimes.The QDE-MNP system is characterized by an enormous mismatch in the lifetimes of its components due to inevitable light absorption in metals, with the LSP lifetime being of the order of a few femtoseconds in the visible region. 11It is clear that, for the strong QDE-MNP coupling regime to be realized, the normal mode splitting should be very large in order to exceed the LSP linewidth or, alternatively, for the Rabi frequency to exceed the LSP damping rate.Considering this requirement from the viewpoint of the QDE decay rate enhanced near a resonant MNP, 3,4,[7][8][9][10] the MNP-enhanced QDE decay rate should be increased by many orders of magnitude, reaching the LSP damping rate. 11The situation changes when considering the interaction of an ensemble of QDEs with plasmonic resonators due to the square-root scaling of the coupling constant with respect to the QDE concentration. 12Several successful demonstrations of the strong-coupling regime in various resonant plasmonic congurations interacting with QDE ensembles have been reported during the last 10 years, 1 including recent real-time observations of ultrafast Rabi oscillations between J-aggregate excitons and surface plasmon polaritons supported by a periodic metal nanostructure. 13However, for the QDE-MNP system, the strong-coupling regime requires very large QDE dipole moments ([1 a.u.) and very small ((10 nm) QDE-MNP separations, 10 and is yet to be experimentally realized.
Numerous investigations considering the modication of spontaneous emission in various QDE-MNP congurations have focused on the effect of dramatic enhancements of the QDE (radiative and nonradiative) decay rates in the vicinity of the MNP, 3,4,[7][8][9]14,15 always implicitly assuming that the relaxation dynamics is purely exponential as obtained in the Weisskopf-Wigner treatment of an individual two-level atom. 16 Atthe same time, the occurrence of transient effects at early times was noted in the consideration of molecular dynamics modied by the presence of an MNP, stressing the following (at later times) exponential decay in the case of weak coupling but without further analysis of transient behaviour.10 Here, we argue that, given an extremely large difference between the relaxation times of the excited QDE and the LSP mode of the MNP, the experimentally accessible regime of weak coupling, in which the QDE-MNP relaxation time is much shorter than that of the QDE in free space but much longer than the LSP lifetime, occurs via coherent interaction between an optically active excited QDE state and a coherent state of the resonant LSP.Recently, we considered a similar system of a generic three-level QDE located near an MNP exhibiting a dipolar resonance (that coincides with the frequency of the QE radiative transition) using a quasi-classical description, in which the QDE was represented by the coherent superposition of the excited and ground states obeying the time-dependent Schrödinger equation (for a two-level system), while the LSP eld was considered to be induced in the MNP by the classical electromagnetic eld created by the oscillating QDE dipole.17 For this conguration under pulsed excitation, we have found an intermediate regime of relaxation (which we associated with the selfstimulated QDE transition 17 ) that occurs in a coherent fashion (like in the strongcoupling regime) with a substantial delay of the emission followed by exponential relaxation dynamics akin to that of modied spontaneous emission in the weakcoupling regime.3,4,[7][8][9][14][15][16] In this work, we develop further our treatment by elaborating the description of the LSP eld and introducing the density-matrix formalism for essentially the same conguration as was treated previously.17 The considered MNP-QDE conguration exhibits the following three main features: (1) Resonant excitation of the LSP is realized as a result of free electron oscillations in the MNP that, for nm-size MNPs, can be regarded as classical current oscillations, since a large number of free electrons ($100 nm À3 ) are involved and their energy spectrum can be considered continuous.This classical oscillating current is then represented by a quantum coherent state 18 of the LSP.Note that the coherent LSP state is fundamentally different from (oen considered) LSP states with a denite number of quantized plasmons. 12) The LSP decay rate is larger than the QDE spontaneous emission (decay) rate by ve orders of magnitude.Therefore, even for relatively strong QDE-MNP interactions, the relaxation of the QDE-MNP system is much slower than the LSP decay, a feature that allows one to disregard the LSP dynamics and consider the MNP response as being instantaneous.
(3) The magnitude of a dipole moment associated with the radiative QDE transition is one order of magnitude smaller than that of an LSP dipole moment induced by the QDE, a feature that allows one to consider the MNP acting as an antenna of the QDE-MNP system.
In our opinion, these distinct features of the considered MNP-QDE conguration open the possibility for the realization of unique (coherent) relaxation dynamics that exhibit characteristic signatures of both weak and strong coupling regimes.
Under these conditions and in the absence of external illumination, it becomes crucial to properly take into account self-action of the excited QDE, in which its dipole eld generates an LSP mode that acts back on the QDE, thus providing a feedback in the QDE-MNP system.Considering the resonant QDE-MNP coupling to be strong enough to signicantly decrease the excited QDE lifetime but weak enough so that the emission rate is still far smaller than the LSP dissipation rate, we nd that the QDE-MNP emission comes with a delay producing a single emission pulse with characteristics that can be controlled by adjusting the system parameters, such as the QDE-MNP separation and MNP dimensions.We further consider energy dissipation in the QDE-MNP system due to the light absorption by the MNP and obtain a general result showing that the probability of the photon emission by the system is $0.75, a number which, rather surprisingly, does not explicitly depend on the system parameters, including the metal absorption.We also discuss several publications, whose theoretical and experimental results might have already (albeit indirectly) supported our ndings.

Theoretical framework
The QDE-MNP system under consideration is schematically presented in Fig. 1 and consists of a generic three-level QDE 9,10 and a spherical MNP.It is assumed that an external pump laser brings the QDE from the ground state g into the state s, where it decays nonradiatively into the optically active excited state e, and that the spherical MNP exhibits a dipolar LSP resonance at the frequency u 0 of the radiative (dipole-allowed) transition e / g [Fig.1(b)].This allows us to separate the excitation dynamics, which are not inuenced by the presence of the MNP, from the relaxation dynamics of state e, whose modication due to the QDE-MNP coupling is the main subject of this work.Note that the shape of an MNP is not important in this context and can be chosen specically in order to produce a dipolar resonance at a given frequency, 10,19 for instance, to coincide with the QDE radiative transition frequency.
The Hamiltonian of the system can be represented as: (2) Here, u 0 is the frequency of the radiative QDE transition, ħ is the reduced Planck constant, b + e and b e are the creation and annihilation operators of the excited QDE state e, b + g and b g are the creation and annihilation operators of the ground QDE state g, u p is the frequency of the resonant LSP excitation, a + and a are the creation and annihilation operators of the LSP, and q is the coupling constant characterizing the interaction between the QDE and LSP.By using the unitary transformation: U 0 ¼ exp(ÀiH 0 t/ħ), we transform the system Hamiltonian into  one that, within the rotating wave approximation and under the condition of strict resonance (u 0 ¼ u p ), has the following form: The rst two of the QDE-MNP system features described in the previous section allow one to signicantly simplify the solution procedure.Considering free electron oscillations in the resonantly excited MNP as classical current oscillations (due to the very large number of electrons involved and the continuity of their energy spectrum), we further make use of the concept of coherent states in quantum optics 18 for the description of this classical current.We also assume that, due to an extremely large difference in the decay rates of an LSP and an isolated QDE, it is possible and, indeed, highly probable that the relaxation of the QDE-MNP system is much slower than that of the LSP, but much faster than that of the isolated QDE.In such a situation, one can neglect the QDE relaxation due to its spontaneous emission and disregard the LSP dynamics, considering the MNP response as instantaneous.Under these conditions, the wave function of the full QDE-MNP system can be represented as follows: where the LSP relaxation rate G is introduced.Here, we operate under the assumption of instantaneous MNP response: va/vt ( Ga, hence Finally, one can work out the following solution of the system of eqn ( 6)-( 9): and r (s) ee and r (s) gg are the probability density matrix elements at the initial moment t ¼ s.
It is seen that the relaxation dynamics of the QDE-MNP conguration is characterized by the parameter m that inuences both the time delay in reaching the maximum QDE transition rate [eqn (11)] and the relaxation rate at later moments of time [eqn (10)].This important system characteristic is related to the basic parameters of the considered conguration.Let us assume that the QDE dipole moment is collinear with the QDE-MNP axis and that the MNP center-to-QDE distance R is considerably larger than the MNP radius r [Fig.1(a)], with all dimensions being much smaller than the wavelength l of light, i.e., that l [ R [ r.Then, by making use of the electrostatic approximation we obtain for the coupling constant q, which represents the energy of the QDE-MNP (dipoledipole) interaction, the following expression: where 3 0 is the vacuum permittivity, 3 d is the relative permittivity of the dielectric environment, d eg is the dipole moment of the QDE radiative transition, and d p is the dipole moment associated with the LSP transition.In the vicinity of the LSP resonance, the polarizability d of the MNP can be expressed as follows: where Finally, the characteristic relaxation rate m of the QDE-MNP system can be expressed via the system parameters by substituting eqn ( 13) and ( 15) into eqn (12): The QDE relaxation dynamics are strongly inuenced by the excitation of the LSP mode that opens a very efficient relaxation channel.The QDE-induced dipole moment of the MNP, which plays the role of an antenna in the QDE-MNP system, can be represented in the following form: d eg e Àiu 0 t coshfmðt À sÞ À 4g : The above expressions [eqn (10), ( 11), ( 16) and ( 17)] constitute the main theoretical outcome of our work, providing simple analytical formulas for the QDE-MNP relaxation and emission dynamics and demonstrating that the dynamics are in general quite complicated.Note that these expressions are also consistent with those obtained previously by us using the quasi-classical approach, in which the QDE was represented by the coherent superposition of the excited and ground states obeying the time-dependent Schrödinger equation (for a two-level system), while the LSP eld was considered to be induced in the MNP by the classical electromagnetic eld created by the oscillating QDE dipole. 17In this respect, the above derivation can be considered as a justication of the previously used approach.

Results and discussion
One of the most important assumptions made in our theoretical framework described in the previous section is related to the strength of the QDE-MNP coupling, which should ensure considerably larger relaxation rates m than that for the QDE in free space g 0 .Their ratio can be evaluated now with the help of eqn ( 16) and the Weisskopf-Wigner result 16 as follows: with l 0 being the vacuum wavelength corresponding to the QDE transition frequency u 0 .For a typical dielectric environment with 3 d ¼ 2.25 (e.g., glass or polymer), the resonance condition (i.e., 3 0 m ¼ À4.5) is met, for gold, at a wavelength of $530 nm with 3 00 m(g) ¼ 2.35 and, for silver, at $400 nm with 3 00 m(s) ¼ 0.22. 22onsidering an MNP with a radius of 5 nm and the QDE distance to the MNP center being 15 nm (in order to be within the electrostatic dipole description), one obtains the ratio b z 17 for gold and b z 77 for silver, thereby justifying the above assumption, m [ g 0 .It is interesting that the effect is already pronounced at relatively large ($10 nm) distances between the QDEs and the MNP surface, which are in the range of distances explored in the recent experiments with 10 nm-size gold nanoparticles. 8It is also apparent that even larger ratios can be achieved by exploiting the LSP shape dependence 10,19 and red-shiing the MNP resonance towards smaller metal absorption. 22he QDE relaxation described by eqn ( 10) and ( 11) begins at some moment in time when the QDE optically active state e, which is created by nonradiative decay of the excited state s, is partially relaxed into the ground state g, so that r (s) gg > 0. This starting process can occur due to other inducements always found in an open system, for example, due to the free-space spontaneous emission, i.e., without the inuence of the QDE-MNP interaction, because the MNP dipole moment [eqn (17)] is still negligibly small.Let us now assume that the QDE relaxation at initial moments of time can be described via an exponential decay, so that r r ee ¼ exp(Àgt).It is reasonable to expect that this initial QDE relaxation is very similar to that of the free-space spontaneous emission, i.e., that g $ g 0 .Applying the continuity condition at the transition between these two processes to both functions, r r ee (t) and r ee (t), and their derivatives, one can determine the Note that the time s does not depend on the initial QDE relaxation rate g, aer which the investigated process starts.Therefore, the present consideration allows one to analyze the whole relaxation process.The QDE relaxation dynamics depend strongly on the efficiency of the QDE-MNP interaction, which is characterized by the relaxation rate ratio m/g.Sharp step-like behavior, which is observed for very large ratios m/g, changes to more gradual population decay for smaller ratios m/g [Fig.2(a)].Emission of radiation from the QDE-MNP system is determined primarily by the square magnitude of the (induced) MNP dipole moment [eqn (17)], whose maximum is attained aer a certain delay time t d ¼ [1 + ln(2m/g)](2m) À1 .For strong QDE-MNP interactions (very large m/g), the emission peak is narrow and occurs close to the initial moment of time [Fig.2(b)].Note that, in practice, m/g $ b [eqn (18)] that can attain rather large values in realistic conditions as discussed above.It should also be noted that the delay time t d is of the same order of magnitude as the width (at half maximum) of the emission peak.This interesting feature might be found useful when conducting and analyzing the corresponding experiments.
The described QDE relaxation process can be considered as the self-stimulated QDE transition from the excited coherent superposition state into the ground state because it is stimulated not by an external (to the QDE) monochromatic light but by the feedback eld from the LSP, which is in turn excited in the MNP by the QDE dipole moment.It can be shown that this eld represents in fact a p pulse, which ensures the QDE transition into the ground state.
We would like to mention that the occurrence of transient effects at early times was noted in the theoretical consideration of molecular dynamics modied by the presence of an MNP, stressing the following (at later times) exponential decay in the case of weak coupling but without further analysis of transient behavior. 10owever, a numerical analysis based on the SPP quantum-mechanical description indicated the occurrence of the emission peak with a certain delay in time in a fashion similar to our results (cf.Fig. 2 here and Fig. 2 in ref. 10).It is also relevant to mention the reported observation of non-Markovian dynamics of a quantum dot resonantly coupled to a micropillar cavity, resulting in a  nonexponential decay in time. 23Finally, we believe that such a delay could have actually been present (but not elucidated) in the recent experiments with 10 nmsize gold nanoparticles connected by DNA to individual uorophores (see Fig. 1(d) in ref. 8).
As already mentioned, the MNP serves as a radiative antenna in the QDE-MNP conguration.The energy dissipation in this conguration (that causes photon loss) occurs due to inevitable ohmic losses in the MNP.Then the probability of non-radiative decay of the system can be determined by relating the energy dissipated in the system 24 to the photon energy: where E is the electric eld induced at the MNP, which is characterized by the dipole moment p [eqn (17)], by the QDE dipole moment d eg : In the discussed case, when r (s) ee [ r (s) gg , aer performing integration in eqn ( 19) one obtains a very simple formula for the loss probability: P y 0.25[1 À (3g/8m)] z 0.25.Under the condition of the resonant QDE-MNP coupling, the transition from the QDE optically active state e to the ground state g occurs during the time interval of $2t d (Fig. 2), implying the transition rate of $0.5/t d .Assuming that the excitation of the upper state s by the pump [Fig.1(b)] together with the subsequent relaxation to the QDE optically active state e takes considerably less time (also that the pump rate is smaller than any of the system relaxation rates), we can evaluate the uorescence enhancement due to the QDE-MNP resonant coupling: which is substantial [eqn (18)] and can easily reach several orders of magnitude.The most interesting physical nding of our work, viz., the spontaneous emission delay under pulsed excitation, is in fact quite general.Indeed, it is only required that a QDE is placed near a resonator that, at the frequency of QDE radiative transition, features a well-dened dipolar resonance with the damping rate, which is substantially larger than the QDE relaxation rate in free space.The former justies the derivation of eqn (4) and ensures the p/2 phase delay in the QDE-MNP system (electromagnetic) feedback, while the latter is needed to realize the desirable weak-coupling regime.The appropriate interaction can be realized at practically any wavelength with non-spherical MNPs 19 or by using low absorbing dielectric (semiconductor) nanoparticles having large permittivity values and supporting strong Mie resonances that can be chosen propitiously by adjusting the particle shapes and sizes. 25Another possibility would be to place a QDE near a metal surface, a conguration that is resonant if 3 0 m (u 0 ) z À3 d , with the strong-coupling regime requiring sub-nanometer QDE-surface distances. 26ur approach can also be applied in this case, provided that |3 with R being in this case the QDE-surface distance, and considering that the QDE dipole moment is perpendicular to the metal surface.Note, that the R À6 scaling in eqn ( 16) is transformed, for this conguration, into the R À3 scaling [eqn (22)], which is also expected to be the case for small separations with the dipolar MNP response to the homogeneous eld becoming strongly multipolar and approaching that of a at metal surface. 27The condition 3 0 m (u 0 ) z À3 d denes the physical situation, where the predicted relaxation dynamics can be observed.In the visible region of the spectrum, where the magnitude of the real part of the dielectric constant of metals, such as gold and silver, is large enough, semiconductor substrates can be used.In this case, excitons, which are typically localized in the vicinity of the interface, can be playing the role of individual QDEs.Indeed, very recent time-resolved spontaneous emission measurements in the conguration, consisting of a semiconductor (CdS) nanocrystal separated from a metal (silver) surface by a 8 nm-thin dielectric (MgF 2 ) layer, revealed that the spontaneous emission reaches its maximum with a signicant delay in time. 28This delay (a few hundred ps) is signicantly larger than the characteristic time of reaching the equilibrium state (a few ps), with the overall time dependence being quite similar to our results (cf.Fig. 2 here and Fig. 3 in ref. 28).

Conclusions
In summary, we have considered the relaxation dynamics of a generic QDE excited with short pump pulses and located near a MNP that exhibits a dipolar LSP resonance at the frequency of the QDE radiative transition.It has been shown that, under the condition of the QDE-MNP characteristic relaxation time being much shorter than that of the QDE in free space but much longer than the LSP lifetime, the QDE relaxation dynamics can be described analytically and feature, in general, non-exponential decay with complicated transient behaviour.Our theoretical analysis resulted in the following main conclusions: (i) the relaxation dynamics in the resonantly coupled QDE-MNP system exhibits step-like behaviour thereby deviating signicantly from the generally accepted exponential decay, 3,4,[7][8][9] (ii) the QDE-MNP radiation emission reaches its maximum with a signicant delay in time, and (iii) energy dissipation in the QDE-MNP system is relatively weak with the probability of the photon emission being $0.75, a number which, rather surprisingly, does not explicitly depend on the metal absorption characteristics.A large number of system parameters in our analytical description opens new possibilities for controlling the QDE relaxation and emission dynamics.Given the variety of resonant plasmonic 19 and semiconductor 25 nanoparticles, the experimental observation of the predicted effect seems feasible, 28 while the possibility of tuning the delay time by changing the QDE-MNP separation can be exploited in many applications, e.g., for practical implementation of a nanoscopic ruler, 29 or optimization of scanning single QDE uorescence lifetime imaging, 30 as well as in fundamental studies within quantum plasmonics.This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.

Fig. 1
Fig. 1 Schematic of a system with a QDE placed near an MNP, indicating (a) system parameters and (b) QDE energetic levels along with an oscillating current associated with the LSP excitation.
n¼0 a n ðtÞ ffiffiffiffi n! p jni: (5) Here, b e (t) and |ei are the probability amplitude and the wave function of the excited QDE state, b g (t) and |gi are the probability amplitude and the wave function of the ground QDE state, a(t) is the eigenvalue of the operator a and |ni is the LSP wave function corresponding to the energy eigenvalue nħu p .Eqn (4) and (5) represent the starting point of our theoretical framework, allowing us to nally obtain relations for the probability density matrix elements of the QDE transitions and for the eigenvalue a(t) of the LSP operators:

00m
is the metal permittivity.Using the condition of strict resonance u 0 ¼ u p and, alternatively, |3 0 m + 23 d | (

Fig. 2
Fig. 2 Relaxation dynamics of the QDE-MNP configuration for different ratios of the established and initial relaxation rates m/g, showing (a) the population decay of the excited state e [eqn (10)] and (b) the squared magnitude of the MNP dipole moment [eqn (17)] normalized to its maximum value.Insets display the same dependencies on a logarithmic scale for the case of m/g ¼ 20.