# Quantum dynamics of two quantum dots coupled through localized plasmons: An intuitive and accurate quantum optics approach using quasinormal modes

## Abstract

We study the quantum dynamics of two quantum dots (QDs) or artificial atoms coupled through the fundamental localized plasmon of a gold nanorod resonator. We derive an intuitive and efficient time-local master equation, in which the effect of the metal nanorod is taken into consideration self-consistently using a quasinormal mode (QNM) expansion technique of the photon Green function. Our efficient QNM technique offers an alternative and more powerful approach over the standard Jaynes-Cummings model, where the radiative decay, nonradiative decay, and spectral reshaping effect of the electromagnetic environment is rigorously included in a clear and transparent way. We also show how one can use our approach to compliment the approximate Jaynes-Cummings model in certain spatial regimes where it is deemed to be valid. We then present a study of the quantum dynamics and photoluminescence spectra of the two plasmon-coupled QDs. We first explore the non-Markovian regime, which is found to be important only on the ultrashort time scale of the plasmon mode which is about 40fs. For the field free evolution case of excited QDs near the nanorod, we demonstrate how spatially separated QDs can be effectively coupled through the plasmon resonance and we show how frequencies away from the plasmon resonance can be more effective for coherently coupling the QDs. Despite the strong inherent dissipation of gold nanoresonators, we show that qubit entanglements as large as 0.7 can be achieved from an initially separate state, which has been limited to less than 0.5 in previous work for weakly coupled reservoirs. We also study the superradiance and subradiance decay dynamics of the QD pair. Finally, we investigate the rich quantum dynamics of QDs that are incoherently pumped, and study the polarization dependent behaviour of the emitted photoluminescence spectrum where a double-resonance structure is observed due to the strong photon exchange interactions. Our general quantum plasmonics formalism can easily be extended to include multiple QDs interacting through the QNMs of metallic resonator structures, fully accounting for radiative and non-radiative coupling, as well as nonlinear light-matter interaction processes.

###### pacs:

42.50.Nn, 73.20.Mf, 78.67.Bf## I introduction

Quantum nanophotonics is an active field of research, which is driven in part by fundamental ideas in light-plasmon-matter interactions, applications in nanophotonics, and by continued advances in nanofabrication technologies. In particular, various types of nanoscale structures have been designed and fabricated to manipulate the interaction between quantum emitters and local electric fields, which can be enhanced by tailoring the local density of optical states (LDOS) typically through some discrete cavity resonance (1); (2). For coupling to electric fields below the diffraction limit, metallic nanoparticles (MNPs) have been shown to yield an unprecedented confinement of light at the nm scale due to the spatial structure of the localized surface plasmon (LSP) resonances. For quantum dot (QD) emitters or artificial atoms placed sufficiently close to the MNP, the strong coupling regime has also been shown to be experimentally accessible (3); (4); (5); (6). As a consequence of the extreme spatial confinement of the LSP, the corresponding effective mode volume of the electric field is much smaller than more traditional dielectric cavity structures, which leads to a strong enhancement of the spontaneous emission (SE) rate in the weak-to-intermediate coupling regime (7); (8); (9). Moreover, fine spatial control of the QD dynamics at the single quantum excitation level and processing of the light signal at the nanoscale is possible (2); (10); (11), resulting in a broad range of applications in fields such as high precision quantum information processing and quantum computation, efficient solar cells (12), and high precision chemical or biological detection (13); (14); (15).

Although long distance and large scale transmission of information using metallic structures is typically not practical because of the strong inherent Ohmic losses of metals at optical frequencies, improvements can be made by using hybrid semiconductor-metallic structures (16); (17), in which the transmission is carried out by semiconductor optical technology, while a nanoscale metallic resonator can be used as an effective modulator (10) and/or transistor (18). Thus it is of significant fundamental and applied interest to study the interaction between dipole emitters such as QDs and individual MNPs. Recently, there have been a few works studying how quantum emitters couple to LSPs, e.g., to describe enhanced SE (Purcell effect) (7); (8); (9), entanglement dynamics (19), and the fluorescence spectrum (20). In the classical or semiclassical regime, with the exception of a particular type of simple geometry such as a spherical cavity (21), for which classical analytical results are available, most of the nanoplasmonic studies are carried out by numerical analysis which is numerically cumbersome and not physically intuitive (22); (23); (24). In the quantum optics regime of cavity-QED (cQED), it has been common to exploit a standard cavitylike master equation with phenomenological decay rates that are implicitly Lorentzian in their decay dynamics (25); (26); (27); such an approach is useful and easy to understand, but it is ultimately limited since the general non-Lorentzian nature of the LSP is neglected, and it it not clear how to obtain the various coupling parameters, e.g., as a function of QD distance from the MNP resonator.

As a result of the fluctuation-dissipation theorem, in general a continuous mode theory instead of the simple single mode theory needs to be employed for a quantum optics description of an inhomogeneous lossy structure, and a quantum noise term can be included phenomenologically (28); (29) or by including a continuous reservoir at the level of a microscopic theory (30); (31); both of these approaches result in a powerful framework with the continuous response of the medium embedded in the medium Green function, which is obtained from an electric dipole source in Maxwell’s equations. On the other hand, it is highly desirable to be able to describe the physics of LSPs in terms of one or a few discrete modes, which has been the standard approach in dielectric cQED systems. Recently, it is shown that the LSP can be effectively described as the quasinormal modes (QNMs) of the MNP (32); (33), which are defined as the eigenfunctions of wave equation in the frequency domain with open boundary conditions (34); (35). A generalized mode expansion technique of the classical photon Green function based on QNMs has been shown to work extremely well for various shaped MNPs, and the SE enhancement of an electric dipole located both inside and outside of MNP shows excellent agreement with full numerical calculations over a broad range of frequencies around the the LSP resonance (37); (36); (38). The combination of an insightful QNM approach and a rigorous Green function approach to quantum optics is thus highly desired, as MNPs facilitate a coupling regime, in general with a non-Lorentzian spectral density, i.e., beyond a dissipative Jaynes-Cummings (JC) model. In certain limits, it can also be used to aid a JC model and justify when such a simpler model can work, with clearly identified coupling rates that can be obtained from QNM theory. Indeed, the MNP yields a rich mode coupling regime as a function of position and polarization, and allows one to explore a complex interplay of radiative and nonradiative dynamics that are unique to the metal environment.

In this paper we present a quantum optics framework to model the quantum dynamics between two QDs coupled to the LSP of a MNP system. The extension to model more than two QDs is straightforward and also described. While there have been several papers studying the dynamics of two QDs coupled by a MNP (39); (40), these approaches, similar to the methods mentioned above, start from assuming the system could be described by the standard cQED master equation by adding in phenomenological decay parameters by hand; by doing so, they neglect the possible non-Lorentzian features of the LDOS which is important in the case for QDs that are sufficiently close to the LSP resonator (20), and they do not incorporate the full electromagnetic response of the MNP environment, including both radiative and nonradiative coupling effects. Instead of assuming a standard Lorentzian decay rate of the LSP, we start from a microscopic model and derive a master equation that takes into consideration the electromagnetic response of the MNP in detail by exploiting a QNM expansion technique for the photon Green function (37); (36). As an example application of this theory, we consider two QDs in the vicinity of a gold nanorod, as shown in Fig. 1(a). While other MNP shapes can also be used in our theory, including metal dimers (36), the single nanorod is partly motivated by the following reasons: () the LSP resonance is around 1.4 eV, which is close to the wavelengths used in optical communication and for many QD emitters, () the nanorod is a nontrivial geometry for which analytic methods are not readily available, and it is dominated by a single cavity mode, polarized along the axis of the rod. While the technique we exemplify below is a single mode theory, it can easily be generalized to include multiple LSP modes if there are several QNMs in the frequency regime of interest, and it properly includes the QNM dissipation. Recently, Yang et al. have carried out a somewhat similar effort, to study the simple linear optical properties of a single dipole coupled to a metal resonator (41).

The layout of our paper is as follows. In Sec. II, we present our main theoretical technique and derive a quantum master equation based on a rigorous quantum optics approach for the medium in terms of the photon Green function, which is obtained from the QNM of the LSP. In Sec. III, we compare the QNM technique to the JC model, and present the improvements over the standard JC model with the help of the QNM technique. We also discuss how our approach could be used in conjunction with the driven JC model in certain regimes, providing a rigorous definition for the various coupling parameters. In Sec. IV.1, our first example studies the simple SE dynamics from a single QD located around the metal nanorod, and shows that the non-Markovian dynamics is important on a time scale of around the lifetime of the LSP. In Sec. IV.2, we study the free-field dynamics of two QDs in a homogeneous background, coupled by the nanorod LSP and show that two qubit (QD) entanglement can be established within a few picoseconds for separate states with a peak value larger than 0.7, despite the strong Ohmic losses; we also study the affect of QD pure dephasing on the peak value of the entanglement evolution, and investigate the concurrence (as a measure of entanglement) for different QD distances from both sides of the nanorod. In Sec. IV.3, we study the incoherent spectrum for the excited two-QD system; in particular, we show explicitly how the real and imaginary part of the Green function contributes to the coupling between spatially separated QDs, and find a rich polarization-dependent behavior of the spectra, including a double-resonance feature which is mediated by the strong photon exchange effects. We present our conclusions in Sec. V.

## Ii Theory: quantization scheme for a general medium, master equation, and Green function expansion in terms of quasinormal modes

For our MNP we consider a 3D gold nanorod as shown in Fig. 1(a) with two QDs (dark brown/light grey ellipsoids), QD and , located around both ends of the nanorod. We use parameters for a metal nanorod with length nm and radius nm, with the Drude model for the dielectric constant, , where rad/s (bulk plasmon frequency) and rad/s (collision rate), similar to the parameters for gold. The background refractive index is . We will also allow for the possibility of an incoherent pump field on the QDs.

For the medium quantization scheme, we start from the formalism developed by Scheel/Dung et al. (42); (29); (43), which has been widely used to describe the quantum electrodynamics of a quantum emitter around a spherical metallic nanoresonator (3); (20); (44). This approach, however, is completely general and can be employed for any lossy inhomogeneous structure as long as the corresponding Green function can be calculated. The photonic Green function is defined through , for the position dependent complex permittivity , where inside the nanorod and elsewhere; here and is the unit dyadic. The imaginary part of the Green function with the same position arguments, gives the projected LDOS (45); (46), while the Green function with different position arguments gives the propagator of the electric field. For a homogeneous dielectric, the imaginary part of the homogeneous Green function is given by .

By treating the QDs as two-level systems, and using the dipole and rotating-wave approximations, the total QD-MNP system is described by the Hamiltonian,

(1) |

where (with ) are the Pauli matrices of the two QDs excitons (electron-hole pairs), and and are the resonance frequency and dipole moment of the -th QD, respectively; are the boson field operators (43), where the electric field operator is given by , with the imaginary part of .

In a rotating frame at the frequency of the QD, , the total Hamiltonian becomes , where the system, reservoir, and the interaction terms are respectively defined by

(2a) | ||||

(2b) | ||||

(2c) |

where . Transforming into the interaction picture, and using the second-order Born-Markov approximation, the master equation for the reduced operator for QD pair is obtained from

(3) |

where is the state of the reservoir; here we have assumed a second-order Born approximation, which is valid in the weak-to-intermediate coupling regime. We assume the temperature of the reservoir is K, which is a good approximation at optical frequencies. The field operators satisfy the following relations: , . After calculating the integrand on the right hand side of Eq. (3) explicitly, we transform back to the Schrödinger picture, and the generalized master equation for the reduced system is obtained as

(4) |

Here , with the photon reservoir function defined through

(5) |

where, for ease of notation, we have introduced with the positions of -th-th QD; in addition, with the unit vector of -th dipole moment, and we have included a pure dephasing term, , on the right hand side of Eq. (4), with a dephasing rate , where is the standard Lindblad superoperator; finally, the last term allows for the possibility of an incoherent pump term on each QD with pump rate . In the following calculations, we assume and D 0.62 e-nm. As can be seen from the time-dependent integral, the the LSP bath sampling depends on the system Hamiltonian. The time-dependent Pauli matrices are given by , and for resonant QDs, so we have . In order to derive Eq. (4), we have used the identity (42). In a single QNM picture, note that can be divergent; however, in our master equation, we calculate , which is convergent. In a practical QNM calculation, we compute the frequency integral over a finite bandwidth that covers the QNM resonance, with limits at approximately 0.5 eV from the QNM resonance frequency. Furthermore, we have checked that this leads to the correct decay rate from a full dipole numerical simulation.

In Eq. (3), some non-Markovian effects are captured through the time integration over the photon reservoir. Applying a second Markov approximation (i.e., ) for the reservoir bath sampling, then we obtain the following Markovian master equation

(6) |

where we have introduced the QD coupling term , a LSP-induced SE rate , a photonic Lamb shift , and a LSP coupling strength between the QDs, . When the QDs are resonant with each other, the coupling term can be simplified to , where describes the coherent coupling between the QDs, and describes the incoherent coupling. The master equation (Eq. (6)) then simplifies to

(7) |

where the effective Hamiltonian term is defined as .

From the above master equations, it is clear that the dynamics of the coupled QDs will show a strong positional dependence through the Green function terms, and this is fully captured in the theory. Unfortunately, the calculation of the Green functions (apart from very simple geometries) is generally a very difficult and a time consuming process, even with computations carried out on clustered computers. In some previous studies, the coupling to the LSP was treated phenomenologically without taking the full detail of the geometry and electromagnetic response into consideration (39); (40). However, for the MNP, recently it has been shown that the Green function can be accurately obtained in terms of an expansion of the QNM, and for the gold nanorod (and indeed many MNP geometries), a single QNM expansion represents an accurate description of the Green function over broadband frequencies and spatial positions. For any two spatial points near the MNP, but outside the regime of Ohmic heating, the dyadic Green function (47) is accurately described through (37); (36)

(8) |

where and are the QNM of interest and the correspondent complex eigenfrequency, respectively. The QNMs are normalized through (34); (35)

(9) |

where in practise we use a computational volume of about 1.5 micron cubed. Alternative QNM normalization schemes are presented in Refs. [(32); (48)], which have been shown to be equivalent (49) to the one above. For our MNP resonator, we have assumed, and verified, that there is only one mode in the regime of interest (near the LSP), and for the gold nanorod, the resonance of the LSP is calculated to be THz (1.344 - i0.0684 eV) (37) with quality factor ; in order to obtain the QNM numerically, a 6-fs length (Gaussian shape in time domain) spatial plane wave near 325 THz with polarization along the axis of the nanorod is injected, and a run-time Fourier transform with a time window 60 fs is employed; a nonuniform conformal mesh scheme is used, with a mesh size of 1 nm cubed is employed around the nanorod. The spatial dependence of the mode profile, , is shown around the nanorod in Fig. 1(b).

In the calculation of the propagator, we use the regularized mode , since it allows one to model spatial regimes from the near to far field regimes (37),

(10) |

with the regularized field given by , where the volume of the integral is now confined to the nanorod volume and . The regularized mode has a simple physical interpretation: it is the solution to a scattering problem when the nanorod is excited by the QNM, which ensures the correct output characteristics in the far field. Usually Eq. (8), which uses the divergent QNM field, gives an excellent approximation to the full Green function as long as the distance between the two positions is no more than a few hundred nm away from the surface of the nanorod; but as the distance becomes sufficiently large, Eq. (10) should be employed to calculate both the propagator and enhancement of LDOS (see Ref. [(37)] for more details).

The enhancement of the projected LDOS, in direction , is defined as

(11) |

and in terms of the QNM contribution, one simply simply replaces by [i.e., Eq. (8)]. Figure 1(c) shows the comparison between the enhancement of the -projected (axis of nanorod) LDOS, , at 10 nm [nm] away from the nanorod [as in shown in Fig. 1(a) by QD], calculated via Eq. (8) (orange solid line) and with a full numerical dipole calculation using finite-difference time domain method (FDTD) (50) (blue dashed line). Clearly the mode expansion technique gives an excellent agreement with the full-dipole FDTD calculation, and thus includes the LSP reservoir function accurately for use in the presented quantum master equation. The total SE rate induced by the QNM, including radiative and nonradiative coupling, is given by

(12) |

where is obtained from Eq. (8).

From the analysis above, it is clear that two spatially separated QDs could be coupled to each other by the coupling term as shown in Eq. (6). In the following, we will also give the emitted spectrum that can be measured at the detector position . For a system containing QDs, the spectrum at the position of the detector, , is given by , with (in the rotating wave approximation). For continuous wave excitation (e.g., from the incoherent pump field), in the presence of just one QD (e.g., QD ), the spectrum is given by

(13) |

However for more than one QD, the the power spectrum is derived to be

(14) |

where , and are the generalized propagator factors from the position of the QDs to the detector, with (or in the present case of two QDs) and ; the incoherent spectrum due to the th QD is defined as ; and the cross term due interference effects between the th and th QD is given by . For convenience, we define the -polarized incoherent spectrum as

(15) |

which we will use later to help explain the polarization features of the emitted spectrum.

## Iii Quasinormal mode model compared to a Jaynes-Cummings model

In the standard JC model both the cavity field and the quantum emitters (e.g., two level atoms) are treated as system operators, which makes the model suitable for studying the physics of strong coupling between the cavity mode and the quantum emitters. To include dissipation into the cavity, the JC model assumes an electromagnetic environment with a Lorentzian spectral density, and this works well for many dielectric cavities. However, the reservoir function from metals resonators can be highly non-Lorentzian; moreover, it is well known that the plasmonic resonance/spectral density of metallic nanoresonators can change as a function of position around the resonator, which can be probed experimentally by measuring the near field electromagnetic response at different positions (51). As is shown clearly through Eqs. (5), (8) and (10), the spectral function of the LSP, which is given by the imaginary part of the Green function at the same spatial point, also depends on the phase of the QNM. Our general model can actually be used to assess when the JC may work, with a rigorous definition of the coupling parameters, and it can go beyond the Lorentzian lineshape model as well. The JC model, if in a regime of validity, can explore effects beyond the 2nd-order Born approximation, e.g., in the strong coupling regime.

The Green function in Eq. (8) can be used to obtain the photon reservoir function [Eq. (5)], e.g., for some dipole position (assumed to be near the resonator), . In a rotating wave approximation, the imaginary part of this function can be written as

(16) |

where and are the real and imaginary parts of the QNM function at the dipole position, and we have introduced a non-Lorentzian coupling factor defined through

(17) |

To better quantify the QNM lineshape, consider a dipole position nm, 10 nm away from the right side of the metal nanorod, as is shown in Fig. 1(b) by the white arrow; the non-Lorentzian coupling factor is , which mainly leads to a small frequency shift of the resonance frequency (which is easily captured in a Lorentzian function by just moving the resonance frequency). Thus for this example, the -projected spectral density obtained from Eq. (5) (orange solid) is well described by a Lorentzian line shape (blue dashed) as shown in Fig. 2(a). In general, however, a position dependent non-Lorentzian spectral density will be obtained around the nanorod; this effect is shown in Fig. 2(b) for at the example position nm, by the solid (orange) line, and a slight blueshift of the resonance peak is also seen which is consistent with the observation in Ref. [(51)]; the dashed (blue) line is the same Lorentzian fit used in the previous case, which clearly shows that the line shape changes as a function of dipole position. The corresponding enhancement of the LDOS/SE is shown in Fig. 2(c). At position nm, the non-Lorentzian shape of the spectral density influences the enhanced LDOS/SE as is shown in Fig. 2(d); the non-Lorentzian coupling factor is now , which has a more dramatic effect on the spectral line shape. We stress that all the information of the resonance shift and non-Lorentzian spectral function is included in the spatial dependence of the phase factor of the QNM; and this information naturally comes into the calculations below through the QNM normalization condition the analytical Green function. Although the Lorentzian fit of the spectral function is valid at certain spatial locations, the non-Lorentzian spectrum becomes important when is large enough, and one then requires the imaginary part of the QNM as well as the real part [(32)].

As discussed above, and shown in Fig. 2(a), at some positions the spectral density could be well described by a Lorentzian line shape (for certain MNPs), so for QDs at these positions the quantum optical interactions could be approximately described by a dissipative JC model with the following QD-cavity coherent interaction terms (in a rotating wave approximation): ; however, we see that these parameters (and the parameters needed to describe QD-QD interactions) actually require the QNM technique in order to have a rigorous definition of these coupling parameters and the normalized field. The single QD-cavity coupling rate will be given by the usual rate , where , while dissipation from the cavity mode is then usually added through a Lindblad operator that describes only Lorentzian decay. For a dissipative Lorentzian model to work, we find that the QD positions must be located around high symmetry points within the vicinity of the field antinode points, but far enough away from the metal surface. Even when the JC model approximately works, then the decay rates still have to be obtained as a function of position in general. This is precisely what the QNM can provide, if the QD position is in a valid Lorentzian decay regime.

We also caution that the JC model still neglects some essential dissipative coupling processes from the metal environment. For example, the standard JC model does not provide an effective description of the nonradiative and radiative decay processes; such a description of the nonradiative/radiative decay will be necessary in order to compute important properties such as the quantum yield (or beta factor), e.g., of a single photon source. Below we demonstrate how one can use the QNM technique to achieve the separation of the total decay rate into radiative and nonradiative decay channels. Moreover, we will also show how one can add in Ohmic losses in a clear and simple way, which is needed for dipole positions very near the resonator (e.g., a few nm from the surface) (36). Importantly, in our approach, all of these physical rates can be computed analytically using the QNM theory, as a function of space and frequency. We describe and exemplify these scattering rates below.

Without the metal nanorod, the background decay rate is simply . While the nonradiative decay rate from the QNM is obtained from (52)

(18) |

where is the induced current density in the nanorod (MNP) at position . Thus the radiative decay rate from the QNM is given by

(19) |

In addition, the quasistatic decay rate can be obtained from

(20) |

with (37); (53) ( is for /-polarized dipoles, respectively). Consequently, the total nonradiative decay is given by . As is shown above, all of the decay rates are highly position dependent, but once the QNM is calculated, the decay rates at different positions can be computed immediately.

Figure 3 shows the decay rates as a function of distance, , away from the surface of the metal nanorod along the -axis; where we see that, at extremely small distances , the system is in quasi-static regime where the Ohmic heating effect due to (magenta dashed) is strong; as becomes larger, all of the decay rates decrease, but the quasi-static decay rate decreases much faster than the others; the inset of Fig. 3 shows the radiative coupling factor, (lower red solid), and nonradiative coupling factor, (upper blue solid), as function of in the near field regime. We see that the radiative output coupling efficiency is below 50%, though this can be increased to about 60% or greater using a dimer configuration (36). While it is not clear how to include such processes in a JC model, which would be further complicated by having different parameters at different QD positions, they could certainly help improve and guide such simpler models. More details on such an approach will be reported in a future publication when we will also explore effects beyond a second-order Born approximation.

With regards to computing the spectrum in a JC model, the spatially integrated far-field spectrum from the cavity operator is usually given by (assuming a rotating frame as the cavity frequency) ; however this assumes that the output coupling rate via the cavity (in this case the LSP) is purely radiative. For a metal resonator, one must include quenching effects by solving the input/output scattering problem, which is exactly what the Green function solution has done. In this way one can compliment the JC model by computing the spatially dependent output spectrum from the QD system operator dynamics and the medium electric field operators, so that . Furthermore, one could obtain the spatially averaged output spectrum (e.g., in the far-field) from , where and the radiative output coupling factor associated with the QNM is obtained from .

To summarize this section, we have discussed how our model can go well beyond the standard JC model while facilitating the simpler JC models in certain spatial regimes. To the extent that the approximate JC could be valid, one still has to obtain the coupling parameters from a model such as ours, and then carefully include quenching effects into any calculation of emitted fields far away from the system resonator. Thus our model can be used to guide and help the simpler JC models in certain regimes as well.

## Iv Results and applications

### iv.1 Localized plasmon induced SE from a single excited QD near the nanorod

Metal nanoparticles enhance the SE rate of excited single QDs due to the coupling with the LSP (QNM). The LDOS at positions around the nanorod changes rapidly in space compared to the homogeneous dielectric structure, and thus the SE rate of a QD around the nanorod can be significantly changed (32); (33); (37); (36), as is shown in Fig. 1(c). In this section, we present an analysis of the SE dynamics of a single QD on resonance with the LSP (; without loss of generality, we take the case of a QD polarized along the -axis at position , 10 nm away from the nanorod as shown schematically in Fig. 1(a) (QD). For the single-QD nanorod system, without an incoherent pump field (i.e., ), the non-Markovian master (Eq. (4)) becomes

(21) |

where we have explicitly used the result since and the kernel function is given by . Note the population decay is not affected by pure dephasing here so we can neglect pure dephasing for this single QD radiative decay study.

We assume here that the QD is initially excited. The QD population decay, , using Eq. (4), is shown in Fig. 4 by the light (green) solid line; the dashed (magenta) line shows the result of a Markovian exponential decay with the rate, given by Fermi’s golden rule with the on-resonant projected LDOS. The inset to Fig. 4 shows that the SE dynamics is recovered by Fermi’s golden rule after a characteristic timescale of about 40 fs (shown in the light gray region); this time scale agrees very well with the corresponding lifetime of the LSP, .

### iv.2 Localized plasmon induced coupling between two spatially separated QDs in homogeneous background

For two spatially separated QDs located around the nanorod (as is shown in Fig. 1(a) by QD and QD), these can be effectively coupled to each other by exchanging photons via the LSP; in the absence of a pump field (i.e., ), the non-Markovian master equation becomes

(22) |

Unless stated otherwise, we will assume the two QDs are resonant with each other (), but may be off resonant with the LSP, where ; however, later we also study the case with different QD resonance frequencies [e.g., in Fig. 11(d)]. Here the intercoupling between the QDs depends on the projected cross density of optical states (CDOS), via , which gives one part of the characteristic coupling strength between the QDs mediated by the electromagnetic environment of the nanorod. Since the Green function in use is the retarded Green function, the real and imaginary parts are related to each other through the Kramers-Kronig relation. The real part of the Green function between the two QDs yields the coherent coupling, , while the imaginary part gives the incoherent coupling, ; the relevant coupling strengths are shown in Fig. 5 with the coherent coupling () and the incoherent coupling strength () given by the chain (orange) and dashed (blue) lines, respectively, and the solid (cyan) line shows . It can be seen from Fig. 5(a), for -polarized QDs at positions nm, that may dominate over and when , where is some offset frequency from the real part of the LSP resonance; however, when QD is -polarized, there is almost no coupling between the QDs; but due to the complex position-dependent polarization characteristics of the LSP, the off-diagonal element of both the projected CDOS, , and the real part of the Green function are non-zero; consequently, QDs with different polarization can be effectively coupled to each other for certain QD positions. Figure 5(b) shows that an -polarized QD at nm could be effectively coupled to a -polarized QD at nm. Below, we will also look at the effect of the coherent exchange interactions in the presence of QD pure dephasing.

We first assume that the two QDs are initially in a separable state (with the first argument for QD and the second one for QD); as a result of the coherent coupling, non-classical correlations will be induced between the QDs, and the quantum correlations approach a maximum value at some characteristic interaction time, which eventually decays to zero due to the decoherence caused by the strong dissipation and the pure dephasing of the system. As a measure of the nonlocal quantum correlations between the separated QDs, we use the concurrence , which is obtained from the eigenvalues of the flipped density matrix