Tunable photonic cavity coupled to a voltage-biased double quantum dot system: Diagrammatic NEGF approach

Tunable photonic cavity coupled to a voltage-biased double quantum dot system: Diagrammatic NEGF approach

Bijay Kumar Agarwalla Chemical Physics Theory Group, Department of Chemistry, and Centre for Quantum Information and Quantum Control, University of Toronto, 80 Saint George St., Toronto, Ontario, Canada M5S 3H6    Manas Kulkarni Department of Physics, New York City College of Technology, The City University of New York, Brooklyn, New York 11201, USA    Shaul Mukamel Department of Chemistry, University of California, Irvine, California 92697, USA    Dvira Segal Chemical Physics Theory Group, Department of Chemistry, and Centre for Quantum Information and Quantum Control, University of Toronto, 80 Saint George St., Toronto, Ontario, Canada M5S 3H6
August 22, 2019

We investigate gain in microwave photonic cavities coupled to voltage-biased double quantum dot systems with an arbitrary strong dot-lead coupling and with a Holstein-like light-matter interaction, by adapting the diagrammatic Keldysh nonequilibrium Green’s function approach. We compute out-of-equilibrium properties of the cavity : its transmission, phase response, mean photon number, power spectrum, and spectral function. We show that by the careful engineering of these hybrid light-matter systems, one can achieve a significant amplification of the optical signal with the voltage-biased electronic system serving as a gain medium. We also study the steady state current across the device, identifying elastic and inelastic tunnelling processes which involve the cavity mode. Our results show how recent advances in quantum electronics can be exploited to build hybrid light-matter systems that behave as single-atom amplifiers and photon source devices. The diagrammatic Keldysh approach is primarily discussed for a cavity-coupled double quantum dot architecture, but it is generalizable to other hybrid light-matter systems.

I Introduction

Recent years have seen a significant progress in probing and controlling hybrid light-matter systems at the interface of quantum optics and condensed matter physics Xiang et al. (2012); Le-Hur et al. (2015); Kurizki et al. (2015); Hincks et al. (2015). Few examples of hybrid quantum systems include cavity-Quantum Electrodynamics (c-QED) arrays Houck et al. (2012); Le-Hur et al. (2015); Underwood et al. (2012); Schmidt and Koch (2013), cold atoms coupled to light Baumann et al. (2010); Kulkarni et al. (2013); Krinner et al. (2015), optomechanical devices Aspelmeyer et al. (2014); Ludwig and Marquardt (2013) and cavity-coupled quantum dots Petersson et al. (2012); Delbecq et al. (2011, 2013); Liu et al. (2014); Kulkarni et al. (2014); Le-Hur et al. (2015); Stockklauser et al. (2015). The motivation for this paper is a class of recent experiments where quantum dots have been integrated with superconducting resonators, accomplishing sufficiently strong charge-cavity coupling of MHz Frey et al. (2012); Petersson et al. (2012); Toida et al. (2013); Viennot et al. (2014); Deng et al. (2015a, b). Such quantum-dot circuit QED systems (QD-cQED) offer a rich platform for studying non-equilibrium open quantum systems at the interface of quantum optics and mesoscopic solid-state physics. Experiments are versatile, with a highly tunable window of parameters. Recent breakthroughs in such devices include the observation of photon emission proceeding via the DC transport of electrons Liu et al. (2014), and the realization of microwave lasers (masers) Gullans et al. (2015); Liu et al. (2015).

Despite ongoing theoretical advances in describing QD-cQED systems Le-Hur et al. (2015); Kulkarni et al. (2014); Schiró and Le-Hur (2014); Dmytruk et al. (2016); Härtle and Kulkarni (2015), there is a compelling need for adapting well-established techniques of non-equilibrium and condensed matter physics, e.g, the diagrammatic non-equilibrium Green’s function (NEGF) method, to explore the rich physics of these highly tunable and versatile hybrid light-matter systems. Existing literature typically treats electron-lead coupling in a perturbative manner (such as the Born-Markov approximation) Kulkarni et al. (2014); Xu and Vavilov (2013a, b); Jin et al. (2011); André et al. (2009); Karlewski et al. (2016), further enforcing the source-drain voltage to be very high, thereby incorporating only sequential, uni-directional electron transfer across the dots. Using such approaches one potentially misses important features in the optical and electronic signals, the result of finite bias voltage and strong dot-lead couplings. Moreover, approximate methods such as the Markovian-secular quantum master equation or mean field calculations are often uncontrolled and non-transparent. It is therefore of a great importance to introduce a systematic approach that allows for an arbitrary dot-lead coupling, (especially since experiments allow tunability from weak to strong dot-lead coupling), handles finite source-drain bias voltage, and treats light-matter coupling in a systematic (even if perturbative) manner. The diagrammatic NEGF approach Ventra (2008); Rammer (2007); Wang et al. (2014); Harbola and Mukamel (2008) is perfectly suited for this purpose. It allows us to simulate present cutting-edge experimental realizations of quantum-dot circuit-QED systems, and furthermore foresee new effects.

Our model includes two electronic levels corresponding to two quantum dots (DQD), each coupled to a primary microwave photon mode (cavity photons). This primary mode is coupled to left and right transmission lines, mimicking the openness of the cavity. A source-drain bias is applied across the DQD system, inducing DC electric current. This system can serve as a testbed for understanding the intricate interplay between light (cavity) and matter (voltage-biased DQD) degrees of freedom, specifically, in a non-equilibrium situation. For a schematic representation, see Fig. 1.

Figure 1: (Left) Schematics of an implementation of a biased double-quantum-dot circuit-QED setup: The electronic-mesoscopic system is integrated with a superconducting transmission line resonator (TLR). The coupling between the microwave photons of the TLR and the electronic levels of each quantum dot (QD), embedded in the nanowire, is described by a Holstein interaction of strength (). The DQD is driven out of equilibrium by the application of a finite source-drain bias . Tunnelling rates between the dots and the electron leads () and in between the dots () can be tuned via gate-controlled tunnel-barriers. Typical experimentally relevant values are given in the Table. (Right) A schematic representation of the model. The mesoscopic system is effectively housed in the microwave cavity. Cavity photons are coupled to the input and output ports with decay rates . The arrows inside the cavity represent tunnelling processes. Wavy lines indicate the light-matter coupling.

The paper in organized as follows. In Sec. II we introduce our model. In Sec. III we study the optical properties of the cavity, namely, the mean photon number, power spectrum and the spectral function, transmission coefficient, and phase response. We use the Keldysh NEGF method while relying on the random-phase-approximation (RPA), which is crucial for respecting symmetry conditions. Numerical simulations demonstrate the potential to use the system and develop novel quantum devices such as microwave amplifiers, photon sources and diodes/rectifiers. In Sec. IV we focus on the electronic part of the model, and demonstrate the influence of the cavity mode on the electric charge current. We summarize our work and provide an outlook in Sec. V.

Ii Model

The double quantum dot setup is placed between two metal leads composed of non-interacting electrons. Electron transfer between the dots takes place via direct tunnelling. Each dot is further coupled to a microwave cavity mode, designated as the “primary photon”. This mode is coupled to two transmission lines, namely input and output ports. The total Hamiltonian is (we set throughout the paper),




where are the site energies of the DQD, coupled to the left and right metal leads by real-valued hopping elements and , respectively. and are fermionic creation and annihilation operators for the respective dots. is the Hamiltonian for the photonic degrees of freedom. It consists of the primary photon of frequency and the secondary photon baths as two long transmission lines () with a symmetric coupling


Here, and are bosonic annihilation (creation) operators for the cavity mode and the two transmission lines.

The interaction between electrons in the junction and the primary optical mode is given by


with as the level number operator and the coupling strength, . For , describes the interaction between the microwave photon and the dipole moment of excess electrons in the DQD. Experimentally relevant parameters are given in the Table. is the decay rate of the cavity mode per port. In the wide-band limit we define it as , where is the bath density of states and being the average coupling between the cavity and the bath modes.

Table: Typical Parameter values from experiments (Refs. Liu et al., 2014; Kulkarni et al., 2014) Cavity loss rate eV MHz Light-matter coupling eV MHz Cavity frequency eV GHz Elastic tunneling eV GHz Drain tunneling rate eV GHz Source tunneling rate eV GHz Temperature mK 0.16 GHz

Iii Properties of cavity-emitted microwave photon

In this section, we compute various experimentally-measurable properties of the cavity, namely its average photon occupation, emission (power) spectrum, spectral response function, as well as the transmission amplitude and phase response of the cavity-emitted microwave photons. It is important to mention that, while performing the phase spectroscopy (transmission amplitude and phase) we relate incoming bosonic modes of the left () transmission lines with the input microwave signal. While outgoing bosonic modes of the right () transmission lines construct the output signal Le-Hur et al. (2015). For other types of photonic and electronic measurements the ports act as a source for dissipation.

iii.1 Average photon number

In recent years, the mean photon number became an experimentally accessible quantityEichler et al. (2012); Liu et al. (2015); Viennot et al. (2014) for QD-cQED setups. We compute the mean photon number in the cavity using the Keldysh NEGF technique. This method allows us to perform a perturbative expansion (second-order) in the electron-photon and cavity-photon bath coupling Hamiltonian, while capturing dot-lead interaction effects to all orders (non-perturbative). We consider the contour-ordered photon Green’s function,


It hands over all components required to calculate various optical signals. Here, is the contour-ordered operator (see Fig. 2) responsible for the rearrangement of operators according to their contour time. The earlier (later) contour time places operators to the right (left). In the second line of Eq. (5), operators are written in the interaction picture with respect to the non-interacting (quadratic) part of the Hamiltonian , for which both electronic and photonic Green’s functions are known exactly.

Figure 2: The complex-time contour in Keldysh formalism. are complex-time parameters. The contour path starts at some initial time , goes to observation time , then comes back to time .

The perturbative expansion of Eq. (5) generates terms of different orders in the electron-photon coupling . A naive perturbative calculation with diagrams up to a particular order leads to the violation of different symmetry-preserving physical processes, such as the conservation of charge and energy currents. In order to restore basic symmetries, one has to sum over an infinite-subclass of diagrams, taking into account all electron scattering events which are facilitated by the emission or absorption of a single photon quanta . This can be achieved by employing the so-called random phase approximation (RPA) Altland and Simons (2010); Utsumi et al. (2006); Agarwalla et al. (2015) where a particular type of ring diagrams are summed over, see Fig. 3. We can represent this infinite summation in a closed Dyson-like (kinetic) equation for ,


is the Green’s function of the primary photon which also includes the effect of the secondary photon modes (transmission lines). corresponds to the bubble diagrams involving the left and right dots’ Green’s function, see Fig. 3. It describes elastic and inelastic (energy exchange) processes, where electrons in the dots interact with the cavity mode. We will later identify the bubble diagrams, in other words, the photon self-energy (connected part), as the density-density correlation function of electrons. In terms of the contour variables, this Green’s function can be written as


This function is symmetric under the exchange of the contour time parameters and . are the non-interacting electronic Green’s functions (dressed by the arbitrarily strong electron-lead tunnelling Hamiltonian), defined as . The average is performed over the current-carrying steady state, determined by the inverse temperatures and chemical potentials of the electronic leads. Components of the non-interacting electron Green’s functions are given in Appendix A. In the third line of Eq. (7) we organize in a matrix form, with and as matrices with . Expressions for different components of and various relations among them such as the Korringa-Shiba relation are explained in Appendix B. It is important to mention that, if the DQDs are further coupled to a phononic environment Gullans et al. (2015), in the bubble diagrams should be replaced by the interacting , dressed by the phononic interaction (assuming Wick’s theorem).

In the steady state limit, different real-time components of can be obtained. The convolution in time domain results in a multiplicative form in the frequency domain. This gives

Figure 3: (a) Dyson equation for the photon Green’s function in contour time. (b) The bubble diagram consists of four dressed electronic Green’s functions, represented by dark solid line (), dashed line (), dotted line () and grey solid line ().

where are the time ordered, anti-time ordered, lesser and greater components of the Green’s function. The primary photon retarded Green’s function is given by


with as different components of the self-energy, materializing due to the coupling of the cavity photon to the input and output ports. In the wide-band limit we approximate (recall that is the decay rate of the cavity mode per port). We also receive and where stands for the Bose-Einstein distribution function, evaluated at temperature .

We now compute various components of the photon Green’s function by inverting the matrix in Eq. (8). We receive the greater and lesser components,


as well as the retarded and advanced Green’s functions


Here, and . In what follows we show that play a central role in enhancing gain in the cavity mode. For later use we also define the total self-energy which is additive in the electronic and transmission lines induced self-energies, i.e.,


With this at hand, we identify the mean photon number in the steady-state limit as


The integration can be performed to include terms to the second order in the electron-phonon coupling (order ) and in . We employ the residue theorem to perform the integration. Upto the second order the poles are located at . Assuming the integration in Eq. (13) then results in


The mean photon number is obtained as


where we alternatively express it in terms of the total self-energy . At equilibrium, the metallic leads are maintained at the same chemical potential () and at the same temperature, equal to the temperature of the photonic environment (). The detailed balance condition is then satisfied for and , i.e., , see Appendix B. This ensures the onset of the Bose-Einstein distribution for the cavity photon mode at equilibrium.

Figure 4: (Color online) Average photon number as a function of (a) applied bias voltage , (b) energy detuning , (c) dot energies and , and (d) dot-lead coupling . Numerical parameters are , , , (dashed), (solid). We use in panel (a), in panel (b), and in panels (c)-(d). Other parameters are the same as in the Table.

From Eq. (15) we can further obtain the mean-square displacement of the cavity mode as


In Fig. 4 we show the average cavity occupation number as a function of applied bias, energy detuning and dot-lead coupling. Unless otherwise stated we define here and below the detuning and enforce , with matter-light coupling . For the spectral functions of the metallic leads we make the wide-band approximation, and fix . We also set the equilibrium Fermi energy of the metal leads at zero and change the bias symmetrically with . The temperature of the two metals and the transmission lines are chosen to be identical.

In Fig. 4(a) we study the average photon occupation as a function of bias voltage. We find that increases in two steps. At the first step, , tunnelling electrons acquire sufficient energy to interact with the cavity mode and generate photons. The second step arises due to the additional resonance situation at . At this bias electrons arriving from the left metal at deposit energy () to the cavity mode, allowing them to resonantly cross the junction. In the positive detuning case , examined here, saturates at lower values for reverse (negative) bias, in comparison to that in the forward (positive) bias. This cavity-number asymmetry with respect to bias reflects the structural asymmetry in the DQD system. With increasing tunnelling strength , fast interdot charge transfer results in effectively weak interaction between electrons and the cavity mode, therefore showing low values for the average photon number.

In Fig. 4(b) we plot as a function of detuning and observe a significant enhancement in the photon number. For the bare system Hamiltonian, strong photon emission into the cavity is expected when , satisfying energy conservation. The position of the peaks in Fig. 4(b) is renormalized with respect to bare values due to dot-lead coupling. Close to the resonance condition a sharp increase in the photon number is observed, potentially creating a lasing function. For large detuning the DQD system does not well interact with the cavity mode, resulting in a vanishing photon occupation. To elucidate this behavior further we display as a function of dot energy levels and in Fig. 4(c). We see similar trends with large photon generation () for .

The nonlinear and non-monotonic behavior of as a function of dot-lead coupling is demonstrated in Fig. 4(d). At weak coupling and finite detuning, increases with , as the photon number in the cavity is amplified by charge transfer through the DQD system. In contrast, when the dot-leads coupling is strong, the renormalization and broadening of the dot energy levels allow electrons to tunnel through the DQD on a short timescale, only briefly interacting with the cavity photons, thus resulting in limited photon generation.

Figure 5: (Color online) (a) Power spectrum as a function of frequency for different dot-lead coupling. (dashed), (solid), (dashed-dotted). (b) Absolute value of transmission as a function of dot-lead coupling for a fixed bias. Numerical parameters are , , , , . Other parameters are the same as in the Table. Circles represent the values of at and at finite but very weak .

iii.2 Power spectrum and spectral function for the cavity mode

Next we look at the power spectrum and the spectral function of the cavity mode. The power spectrum has recently been measured for a similar setup Liu et al. (2015). In the stationary limit, we can immediately obtain the emission (power) spectrum in terms of the lesser Green’s function


It is obvious from this definition that . At the cavity frequency , is given by


Similarly, we obtain an expression for the spectral response function of the resonator, defined as the difference between the retarded and advanced photon Green’s function


with the normalization condition (sum rule) . This can be proved as well from the above equation invoking the residue theorem as explained before. In the absence of the light-matter interaction the cavity spectral function trivially satisfies the sum rule. The amplitude of the spectral function at is related to the emission spectrum as . In Fig. 5(a) we plot the power spectrum for different dot-lead coupling. It shows a nonmonotonic behavior with respect to tunnelling strength, with the maximum value taking place at an intermediate value for the tunnelling. The brodening, which is of the order of several MHz (1-3 MHz), results from the interplay between and , the two different sources of dissipation.

iii.3 Phase spectroscopy: Transmission and Phase

We calculate the transmission amplitude and the phase response of the emitted microwave photons by following the input-output relations in Refs. Clerk et al., 2010; Dmytruk et al., 2016; Schiró and Le-Hur, 2014. Recall that for phase spectroscopy measurements (performed via heterodyne detection Liu et al. (2014))the bosonic modes of the left () and the right () transmission lines are related with the input and output microwave signal, respectively Le-Hur et al. (2015). The transmission function reads as Clerk et al. (2010); Dmytruk et al. (2016); Schiró and Le-Hur (2014)


where can be shown to be proportional to the electronic charge susceptibility as Dmytruk et al. (2016); Schiró and Le-Hur (2014). Here


stands for the electron density response function, are indices for the dots, refers to the average over the electronic degrees of freedom in the nonequilibrium steady state. In our formulation, following Eq. (11) we identify the transmission to be proportional to the retarded Green’s function of the cavity photon mode, which in the time-domain is precisely the photon response function, given by


Therefore, we note that stands for the retarded component of the bubble diagram . Here represents average over the combined photonic-electronic steady state density operator.

Figure 6: (Color online) Spectroscopy of cavity-emitted photons. (a) Transmission (absolute value) and (b) phase response as a function of incoming photons frequency . Parameters are , , , , (solid) and (dashed). The red dashed-dotted line was obtained at . Other parameters are same as in the Table. The sum rule for transmission is .

We further write and identify the real and imaginary parts of the transmission function,


with the sum rule . Specifically, we get and . We further note that the real part of the transmission function provides a direct measure for the spectral function of the cavity photon. The phase of emitted photons is


Since we are mainly interested in the absolute value of the transmission function and the value of the phase response at the frequency of the cavity mode , we evaluate


Note that both the real and imaginary parts of show nontrivial dependence on bias voltage through the Keldysh component of the electronic Green’s function , see Appendix B.

In Fig. 6 we plot the absolute value of the transmission and the phase response as a function of the incoming photon frequency under different bias voltages . When the cavity is decoupled from the DQD, (dashed-dotted line), the transmission reaches unity at , and the broadening is determined by . As well, the phase response is zero at the resonant frequency, and it approaches in the off-resonant regime. For finite —yet at zero bias—charge fluctuations in the dots introduce shift in the transmission peak, further reducing the maximum amplitude. The frequency shift depends on the real part of the charge susceptibility , whereas the broadening reflects the difference between and . The phase response is zero when the transmission is at maximum. Most interestingly, we find that the absolute value of the transmission coefficient can be greatly enhanced—beyond unity—at finite bias, once . This situation is elaborated below.

Fig. 7 displays one of the central results of our work: The transmitted photon signal can be significantly enhanced at finite bias voltage, once the electronic system is fine-tuned to counteract dissipation from the (photonic) transmission lines. We study the behavior of the transmission coefficient and phase response at the bare cavity resonance frequency as a function of bias, for different cavity decay rates . We find that the transmission, or photon gain, increases for , and it saturates at high (positive and negative) biases. This behavior agrees with our observations for in Fig. 4(a). However, at a certain finite voltage (here around 30), the transmission jumps above the asymptotic value (panel a), while the phase response shows a sudden dip (dashed line in panel b). This sudden jump takes place precisely as the electronic part of the system acts to cancel out relaxation effects due to the tranmission lines, , see panel (d).

Figure 7: (Color online) Controlling transmitted photons by a DC bias voltage across the dots. (a) Transmission and (b) phase, as a function of bias voltage at the bare cavity frequency . In panels (c) and (d) we display the real and imaginary parts of the bubble electronic Green’s function , respectively, illustrating that a jump in transmission (and a dip in phase) precisely occurs once (panel d). Parameters are or 1.21 MHz (solid), or 0.24 MHz (dashed), , , , . Other parameters are given in the Table.
Figure 8: (Color online) Effect of detuning the DQD levels on the cavity signal. (a) and (c) Absolute value of transmission coefficient , and (b) and (d) phase response of the transmitted signal, displayed as a function of detuning at high bias using , , cavity loss rate , . in panels (a)-(b), in panels (c)-(d). Other parameters are same as in the Table.

We turn to the high bias limit and study in Fig. 8 the role of energy detuning on the transmission and phase response. In particular, we examine the dependence of and on the electron-photon coupling strength and on the incoherent tunnelling rate . Panels (a) and (c) display the transmission amplitude, showing gain (dip, ) on the positive (negative) side of detuning as a result of coupled electron photon transport processes: For positive detuning and positive bias, electron transport through the DQD system proceeds via inter-dot tunnelling, thereby reducing the energy of electrons via photon emission. In contrast, for negative detuning electron transport proceeds assisted by photon absorption, reflected by the dip in the transmission coefficient. This gain mechanism can be also corroborated with the amplification of . Finally, note that at zero detuning (), direct elastic tunnelling dominates over photon-induced contributions. For very large detuning, marginal charge flow through the dots results in an effectively minuscule electron-photon interaction. These two limits lead to unit transmission amplitude and zero phase response.

The relative strengths of and determines the gain and loss values in the transmission amplitude. Plots of the real and imaginary components of , displayed as a function of , are included in Appendix B (Fig. 13). Maximum gain is achieved when . Since , and typically , increasing shows a significant enhancement in gain and similarly loss.

Figure 9: (Color online) Control of cavity response by tuning the DQD site energies. Plot of (a) absolute value of transmission coefficient and (b) phase response of the transmitted signal as a function of and at high bias . Here , . Other parameters are the same as in the Table.

The dependence of and on the system-lead coupling strength is examined in Fig. 5(b), showing a nonliner behaviour, and in Figs. 8(c) and (d). The transmission is high in the sequential tunnelling regime (intermediate dot-lead coupling), whereas for large coupling renormalization and broadening of peaks lead to reduced gain. The reason is that the dwelling time of electrons in the dots is long ) at weak coupling, realizing an effectively significant electron-photon interaction. Indeed, figure 13 in Appendix B demonstrates that small returns large values for , thus a significant enhancement in as a whole. Upon increasing the coupling strength , the dots energies become broadened, thus electrons flow across the device without interacting with the cavity mode. This scenario shows small gain and loss. At very weak coupling, , the electronic medium introduces only dissipation, responsible for the sharp drop in transmission, see Fig. 5(b)

Plots of and , at the bare cavity frequency, as a function of and , reveal that degenerate quantum dots do not influence the cavity, with and , see Fig. 9. In contrast, at positive (negative) detuning, approximately satisfying , gain (loss) is observed.

In Fig. 10 we study the effect of a finite bias voltage, with energy detuning lying within the voltage bias window, on the transmission coefficient and phase response. Comparing finite-intermediate voltage results to Fig. 8, where high bias was employed, we note here an additional dip at , reflecting photon-assisted charge transfer processes from the right dot to the left dot. Correspondingly, a jump in phase is detected at the same value of detuning.

Figure 10: (Color online) Signatures of finite voltage bias on (a) and (b) , studied as a function of . Here, , , . Other parameters are the same as in Fig. 8.

iii.4 Special limits, scaling and universality

We discuss here scaling relations for measures related to the cavity, in different parameter regimes. We begin by considering the high bias and zero temperature limit. We further assume that the dot-lead coupling is strong, larger than detuning, . This situation is experimentally relevant and thereby potentially testable. We use Eq. (LABEL:bubble-GF) along with the expressions for the non-interacting electronic Green’s functions (see Appendix A). In this limit we found that


with . Further assuming that , the average photon number in Eq. (15) reduces to


with an effective temperature . It is remarkable to note that this linear scaling, , is universal to the Holstein-like class of models Mozyrsky and Martin (2002); Utsumi et al. (2013); Jin et al. (2015); Agarwalla and Segal (2016). In the opposite limit (though keeping ), the electronic part effectively decouples from the cavity. As a result, the cavity equilibrates with the secondary photon bath (ports), , and the transmission amplitude goes to unity.

Another interesting limit is the large detuning and high bias case, where we obtain


In this case, , the electronic current is negligible, and electron-photon coupling is effectively small. Again we find that when , , and in the opposite limit, , the cavity occupation number is thermal.

Iv Electronic properties: Steady state charge current

We study the electronic properties of the DQD system, focusing on the steady state charge current at the left contact. It is given by the powerful Meir-Wingreen formula Meir and Wingreen (1992); Haug and Jauho (1996), valid for an arbitrary large light-matter interaction and dot-lead coupling

Figure 11: Hartree () and Fock () diagrams for electrons in contour time. For expressions see Eq. (31)


Here are the lesser and greater components of the electronic Green’s function, fully dressed by the leads and the electron-photon interaction. As before, we obtain these components by relying on the Dyson equation, performing second order perturbation expansion in the electron-photon interaction Harbola et al. (2006); Lü and Wang (2007). We write


with as the nonlinear electronic self-energy arising due to the photonic degrees of freedom. Calculating it up to we receive the Hartree (H) and Fock (F) terms Park and Galperin (2011); Galperin et al. (2004) (see Fig. 11),


with . Following the Keldysh equation, we gather the lesser and greater components in the frequency domain as , with the total self-energy as the sum of left and right-lead self-energies, as well as the nonlinear component, .

Figure 12: (Color online) (a), (c) Charge current as a function of detuning . The nonlinear contribution to the current is shown separately in panels (b) and (d). Parameters are , , , , with in panels (a) and (b), in panels (c)-(d).

Substituting these expressions into the Meir-Wingreen formula we organize the charge current formula, written as a sum of elastic and inelastic contributions i.e., , with




Note that in the above equations the retarded and advanced components of the Green’s functions are renormalized by the electron-photon interaction. We expand these functions following the Dyson equation, , and organize the lowest order expression for the charge current in terms of the non-interacting . The elastic components of the current becomes , with