# Tunable Hybrid Quantum Electrodynamics from Non-Linear Electron Transport

###### Abstract

Recent advances in quantum electronics have allowed to engineer hybrid nano-devices comprising on chip a microwave electromagnetic resonator coupled to an artificial atom, a quantum dot. These systems realize novel platforms to explore non-equilibrium quantum impurity physics with light and matter. Coupling the quantum dot system to reservoir leads (source and drain) produces an electronic current as well as dissipation when applying a bias voltage across the system. Focusing on a standard model of biased quantum dot coupled to a photon mode which gives rise to an Anderson-Holstein Hamiltonian, we elucidate the signatures of the electronic correlations in the phase of the transmitted microwave signal. In addition, we illustrate the effect of the electronic degrees of freedom on the photon field, giving rise to non-linearities, damping and dissipation, and discuss how to control these effects by means of gate and bias voltages.

###### pacs:

73.23.-b,42.50.-p,72.10.Fk## I Introduction

Recent years have seen tremendous experimental progress in probing and controlling quantum mechanics at the level of single constituents. A variety of physical systems, at the interface between quantum optics, atomic and solid state physics, have emerged as natural platforms to store and process information at the genuine quantum level. Examples include systems of entangled atoms and photons Raimond et al. (2001); Cohen-Tannoudji et al. (1997) and their solid-state analogue Yamamoto et al. (2000); Schoelkopf and Girvin (2008); Akkermans and Montambaux (2007) or systems of electrons and spins confined in semiconducting nanostructures Petta et al. (2005), to name a few.

More recently the experimental effort moved toward the realization of hybrid nanosystems obtained interfacing different kinds of quantum information platforms, such as quantum dots, ultracold atoms in cavities, mechanical resonators or superconducting circuits Mekhov et al. (2007); Majer et al. (2007); Suffczynski et al. (2009); Delbecq et al. (2012); Verdú et al. (2009); Parigi et al. (2012); Ritsch et al. (2013); Marquardt and Girvin (2009); Pirkkalainen et al. (2013); Kubo et al. (2010); Dumeige et al. (2011); Petersson et al. (2012); Frey et al. (2012); Jin et al. (2012, 2013); Miguel-Sánchez et al. (2013); Bergenfeldt et al. (2013); Viennot et al. (2013); Bocquillon et al. (2012); Dubois et al. (2013). These novel hybrid architectures are not only of great practical relevance for quantum information but also of fundamental interest as platforms to explore many body physics in novel and unconventional regimes, particularly away from thermal equilibrium. Indeed when small interacting quantum mechanical systems involving inherently different degrees of freedom, such as electrons and photons, are coupled to their fluctuating quantum environments a plethora of interesting many body phenomena are expected to emerge Goldhaber-Gordon et al. (1998); Cronenwett et al. (1998); Türeci et al. (2011); Latta et al. (2011); Mebrahtu et al. (2012); Jezouin et al. (2012); Hofheinz et al. (2011); Ingold et al. (1992); Safi and Saleur (2004); Le Hur (2004); Le Hur and Li (2005); Le Hur et al. (2007); Chung et al. (2013); LeClair et al. (1997); Le Hur (2012); Orth et al. (2013); Goldstein et al. (2013); Bera et al. (2013); Sbierski et al. (2013); Beugnon et al. (2006); Henriet et al. (2014). For example, in a recent pioneering experiment Delbecq et al. (2011) the coupling between a quantum dot, a microwave resonator and two biased electronic reservoirs was achieved on chip and signatures of many body correlations observed in the phase of the microwave transmitted photon as a function of the voltage bias across the leads. Motivated by these rapid developments Delbecq et al. (2011); Frey et al. (2012), in this paper we discuss the non-equilibrium physics of an hybrid quantum impurity model made by electrons and photons.

The paper is organized as follows. In section II we introduce our hybrid quantum impurity model. Section III is devoted to study transmission of photons through the microwave cavity and to elucidate the signature of many body correlations in the phase of the microwave signal. In section IV we derive an effective theory for the cavity photon coupled to the density fluctuations of the dot electron, while section V we present a discussion of the results.

## Ii Hybrid Impurity Model with Electrons and Photons

We model our system as a single level quantum dot tunnel coupled to biased source-drain leads and capacitively coupled to a single mode electromagnetic resonator (see Fig. 1). The Anderson-Holstein type of Hamiltonian for the system reads (setting )

(1) | |||||

where are annihilation operators of, respectively, leads electrons, correlated dot and photon field, and while . The above Anderson-Holstein model Holstein (1959); Anderson (1975) has been studied extensively in the context of non-linear electron transport, with most of the emphasis focused on signatures of many body correlations through electronic properties such as current-voltage characteristics and conductance Mitra et al. (2004); Cornaglia et al. (2004); Vinkler et al. (2012); Laakso et al. (2014). Here our scope is to probe the many-body Kondo-type correlations of this system through photon transport.

In order to probe this hybrid system we capacitively couple the resonator to two long transmission lines Clerk et al. (2010). Without any loss of generality we consider a symmetric coupling to the two bosonic baths, such that the Hamiltonian of the full system takes the form

(2) |

where describes photonic modes of the transmission lines () with dispersion and spectrum , with being an high-energy cut off. In the next section, using exact equations of motion and the input-output relations Clerk et al. (2010); Le Hur (2012) we will obtain an expression for the transmission/reflection coefficients which allows to characterize photon transport in the system.

## Iii Phase Spectroscopy of a Cavity-Emitted Microwave Photon

### iii.1 Input-Output Theory

Let’s start by deriving an input-output relation for this problem, following the treatment of Refs.[Clerk et al., 2010; Le Hur, 2012]. We consider the full Hamiltonian of the problem and we write the equation of motion of the operator which reads

(3) |

After integration we get for

(4) |

while if we integrate backward in time from we get

(5) |

We introduce the left input voltage at time

(6) |

as well as the output voltage at time

(7) |

Using Eq. (4-5) for we can express the average output voltage in terms of the input field. Indeed we have first

(8) |

then by taking the quantum average from both sides and introducing a susceptibility defined as

(9) |

we end up with

(10) |

where we have introduced the reflection coefficient

(11) |

with

Using the same approach we can obtain an analog result for the transmission coefficient, , defined as

(12) |

The susceptibility is nothing but the response of the photon displacement to an input signal, carried by , and it is given by the Kubo formula

(13) |

where the average is taken over the system hamiltonian with no transmission line but fully including the fermionic leads. As a result, the expression for the transmission and reflection coefficients have to be interpreted as perturbative in the coupling to the photonic bath (). Indeed parameters in figure 2 are chosen such that the scattering matrix is almost unitary, i.e., the inelastic contribution to the transmission of light is small compared to the elastic part.

### iii.2 Phase Spectroscopy

Given the previous result for the transmission coefficient, by defining we can obtain the phase of the transmitted microwave signal,

(14) |

The retarded photon Green’s function can be written in Fourier space in terms of the photon self-energy as

(15) |

where includes both the effects of frequency renormalization and the damping due to the fermionic environment. From this expression we immediately get the frequency dependence of real and imaginary part of the transmission coefficient

and from them we can extract the phase

(16) |

We start evaluating the phase to the lowest order in the electron-photon coupling . To this extent it is convenient to compute directly the photon retarded self-energy , defined in Eq. (15). Within a Keldysh approach this turns to be related to the density response function of a purely electronic Anderson Impurity Model (AIM),

(17) |

with the electronic charge susceptibility. For an Anderson Impurity Model which exhibits a Fermi-Liquid type of ground state this must satisfy the Korringa-Shiba relation Filippone et al. (2011) which implies

(18) |

However much less is known in general about this quantity at finite frequency in presence of both strong interaction and finite bias (but see references [Werner et al., 2009; Schiró and Fabrizio, 2009; Dutt et al., 2011] for recent developments).

### iii.3 Resonant Level Model

To gain some analytical insight we start evaluating the response function for . In this non interacting limit, the dot electrons form a resonance of width centered at , their spectral function in the wideband limit becomes a simple lorentzian, and the associated Green’s functions read

(19) |

and

(20) |

Evaluating the density response function reduces then to a convolution of single particle Green’s functions

(21) |

which can be evaluated in closed form. Here we discuss the resulting behavior for the transmission coefficient. We remind here that we are evaluating the transmission coefficient according to Eq. (12), which is valid in the linear response regime with respect to the coupling to the photonic bath (external transmission line).

In figure 2 we plot the behaviour of the real and imaginary part of the transmission coefficient, as well as its phase, as a function of incoming photon frequency and for different bias voltages . The real part of shows a peak at a renormalized frequency , which shifts and sharpens at larges voltages. Electrons induce a renormalization of photon frequency, which can be seen as a many body light-shift.

The imaginary part of the transmission vanishes at and similarly does the phase , as it follows immediately from Eq.(14). More interesting and rich is the behavior of the phase at the bare resonator frequency, , which is also what typically measured in experiments Delbecq et al. (2011). As we see from figure 3 at large voltages the phase smoothly approaches , quite independently from photon frequency and other system parameters and also from the coupling to the transmission line. This universal behaviour of the phase in the strongly non linear regime can be understood at first by looking at the structure of photon self-energy . Indeed one can write the phase of the transmission coefficient on resonance as

(22) |

We see how a phase of crucially requires the photon damping induced by conduction electrons to vanish sufficiently fast at large bias. Perturbatively we find both the real part and the imaginary part of the self-energy to vanish at voltages larger than a crossover scale as

(23) |

hence, . Corrections to the universal large bias limit depend from the photon frequency (bottom right panel), stronger quantum effects (larger ) shift the crossover scale to larger values, suggesting an effective classical picture to explain this effect as we will see later in the paper.

### iii.4 The Role of Photonic Dissipation

It is interesting to discuss at this point how the above results change in presence of an extra source of photonic dissipation, such as for example that one coming from a small coupling to the transmission lines. In order to include this effect in our treatment we modify the expression of the (bare) photon Green’s function to include some coupling to a dissipative bosonic bath

(24) |

where is the effective dissipation induced by a coupling to the modes in Eq. (2). Here we assumed the photonic bath to be Markovian, with a cutoff frequency . In principle, the computation of the photon self-energy to the lowest order is unchanged, since it only involves electronic propagators. As a result we get for the phase of the transmission coefficient the result

(25) |

In figure 5 we plot the phase as a function of the bias voltage for different values of the photonic dissipation. We see that at finite the phase crossesover away from the value to . A simple calculation using the scaling of the photonic self-energy (23) gives the crossover voltage scale .

### iii.5 Infinite Anderson Impurity Model

In this section we briefly consider the limit of the Anderson Impurity Model (AIM) and discuss in particular the computation of the density-density response which is relevant to obtain the photon self-energy. In the limit we can project out all the states with double occupancy of the fermionic level. This can be done formally by introducing a representation in terms of an auxiliary fermion and a slave boson, , living in a constrained Hilbert space

(26) |

In terms of these new degrees of freedom the AIM Hamiltonian becomes Coleman (1987); Houghton et al. (1987)

(27) |

with the Hamiltonian of a the two free leads. In order to account for non-equilibrium effects due to the voltage bias we wil use a Keldysh slave boson approach Rosch et al. (2001); Ratiani and Mitra (2009) and consider the Non-Crossing-Approximation(NCA), corresponding to the lowest order self-consistent perturbation theory in the hybridization and an exact treatment of the constraint. While in general the NCA equations have to be solved numerically, the limit of large voltage biases has been studied in detail in the context of quantum transport and it is now well understood even at the analytical level Rosch et al. (2001). Here, due to bias-induced decoherence effects the pseudo fermion acquires a finite damping, , which is given for by Rosch et al. (2001)

(28) |

being the Kondo temperature. This immediately translates into a narrow peak in the pseudo-fermion () spectral function, while the slave boson remains incoherent. Since we are interested in computing the electron density response function, which is directly related to the density response of the auxiliary fermion Coleman (1987); Houghton et al. (1987), we can immediately borrow the results obtained previously for the non-interacting case and use them for the case as well, with the difference that now the lifetime of the fermionic excitation is strongly bias dependent. From this we conclude that since for one has , i.e. we are always in the regime where the phase saturates to .

## Iv Integrating Out the Fermions - Effective Theory for the Photon Field

To get a further insight, we now derive an effective theory for the photon in the resonator after integrating out all electronic degrees of freedom. We will show that (i) fermionic correlations induce photon non-linearities as well as noise and dissipation and (ii) that the photon long-time dynamics can be described effectively by a Langevin equation for a noisy Duffing oscillator Kovacic and Brennan (2011). Finally, we will use this effective theory to interpret the result on the phase spectroscopy obtained in previous section.

In order to proceed we formulate the problem in terms of Keldysh action Kamenev (2011).Let’s start considering the photon in the resonator that we describe as an harmonic oscillator with coordinate and momentum and with Keldysh action

(29) |

where is a field defined on the two branches Keldysh contour , or equivalently in the right hand side we have introduced the two components which are defined on the forward/backward branch, respectively. It is useful for what will follow to perform a change of variables and introduce the classical and quantum components of the field defined as

(30) | |||

(31) |

in terms of which the action becomes

(32) |

A few things are worth to be noticed from this simple result. First, the Keldysh action in this basis vanishes when evaluated at , i.e.

This is a general result that follow essentially by the unitarity of the evolution and that will hold also for the more complicated problem of electrons coupled to photons we will consider now. Second important observation is that, performing the functional integral over the quantum field we obtain a functional Dirac-delta which constraint the dynamics on the classical trajectory

(33) |

This result is in general true only perturbatively, in other words expanding the Keldysh action to linear order in will gives us the effective quantum dynamics. We will discuss later the meaning of the higher powers of .

We now consider the full problem, including the fermionic degrees of freedom, and we couple the photon field to the electronic density so that the full Keldysh action reads

(34) |

where and are, respectively, the photon and the AIM action while the last term describes electron-photon coupling

(35) |

Here we introduced the so called classical (c) and quantum (q) basis for other operators as well, such the electron density.

We can formally define an effective action for the photon by averaging over electronic degrees of freedom

(36) |

with

(37) |

This effective theory describes a photon field coupled to the non-equilibrium environment generated by the density fluctuations of dot electrons. While the functional cannot be evaluated exactly, due to the form of the electron-photon coupling, the theory admits a simple analytical description in the large bias regime Mozyrsky et al. (2006); Nocera et al. (2012).

We start noticing that, again, the effective action (36) vanishes when evaluated at . This is because by construction

(38) |

the last identity follows from the fact that in absence of any quantum component for the field we have by construction

(39) |

due to causality. Alternatively, this can be seen explicitly from the expression of the effective action. This suggests to expand the effective action in power of the quantum field around , while keeping for the time being still fluctuating and to stop the theory to the lowest non trivial orders. By proceeding in this way we obtain

(40) |

where is a term which is purely local in time while

(41) |

In the next two sections we will discuss separately the physical meaning of both terms. We notice that there are no linear terms nor quadratic terms in the classical field in this expansion, since

(42) | |||

(43) |

again by construction due to the causality of the theory. As we mentioned in the previous paragraph the solution of this non-interacting problem for a generic time dependent field is still a challenging task to be accomplished in fully generality. However we can take advantage of the fact that for very large bias the dynamics of the electronic degrees of freedom is much faster than the dynamics of the photon, with typical time scales

As a result we can to leading order assume an adiabatic approximation and assume that the electronic degrees of freedom just see a static classical field . In other words we can evaluate the coefficient of the effective action by considering a AIM which is in its instantaneous ground state corresponding to the configuration of the classical field .

### iv.1 Electron Induced Photon Non-Linearities

The first term is local in time and reads

(44) |

with . We notice that, in absence of other terms in , the functional integral over would result in the photon equation of motion from which we can interpret as an effective potential for the photon, that renormalizes the bare harmonic one. Expanding it around we get the general form

(45) |

where represent respectively the renormalized frequency, the anisotropy term and the effective anharmonicity. We can obtain analytical expressions for these coefficients in terms of the parameters of the electronic problem, in particular the external gate and bias voltages. Their expressions can be obtained by differentiating in Eq. (45) and read for , respectively

(46) |

and

(47) |

and finally the anharmonicity

(48) |

The anisotropy term is particularly interesting to discuss. As we see from Eq. (47) this term at zero gate voltage reduces to

(49) |

namely it vanishes for symmetric bias voltages, . This is indeed not surprising since the anisotropy term breaks the inversion symmetry , which is a direct consequence of particle hole symmetry of the Holstein model that holds at and . The results discussed here concern the non-interacting , but can be extended to the interacting case in the large bias regime, as discussed in previous section upon substituting in Eqns. (46,47,48)

### iv.2 Electron Induced Photon Damping and Dissipation

Let’s now consider the effects of the other terms in the effective action which are non local in time. We can distinguish two contributions, and , which include retardation effects and read, respectively

(50) | |||

(51) |

with being respectively the retarded (R) and Keldysh (K) component of the density-density response function of the AIM. We now focus on the low frequency behaviour of the kernels .

We start from the retarded component, encoding the spectrum of density fluctuations. At low frequency, , its imaginary part vanishes linearly, see figure 6

(52) |

where the low frequency slope can be computed analytically for the RLM case, and reads

(53) |

More generally, this result can be also understood by noticing that, upon integrating out the driven cavity, the problem becomes reminiscent of an AC driven quantum dot coupled to fermionic leads Buttiker et al. (1993); Gabelli et al. (2006); Fève et al. (2007); Gabelli et al. (2012) which has been recently shown to possess an effectively ohmic dynamical charge response Mora and Le Hur (2010); Nigg et al. (2006); BÃ¼ttiker and Nigg (2009); Etzioni et al. (2011); Hamamoto et al. (2010); Filippone et al. (2011, 2013); Dutt et al. (2013). A finite bias voltage does not change qualitatively this linear behavior and only renormalize slightly the slope.

As opposite, the Keldysh component is strongly affected by a finite bias voltage which turns the characteristic linear in frequency equilibrium behaviour, imposed by fluctuation-dissipation theorem at , into a finite value at zero frequency (see Fig. 6). We can obtain a full analytical understanding of this in the case of a biased RLM. Using the expression for the lesser and greater RLM Green’s function we can write the Keldysh component as

The low frequency behaviour of this kernel controls the nature of the noise. It is easy to see already from this expression that for the behaviour of depends strongly from the bias voltage. While at zero bias the kernel vanishes as a power law, see fig 2, a finite voltage acts as an infrared cutoff and results into a finite value at zero frequency, , with

(55) |

which plays the role of an effective diffusion coefficient for the dynamics of the photon. A simple calculation gives

Quite interestingly we see that the diffusion coefficient vanishes linearly with the bias voltage as

(56) |

which suggests in some sense that the bias voltage is playing the role of an effective temperature, . For larger values of the applied bias the diffusion coefficient saturates to a constant value, independent on the oscillator position , which reads

(57) |

Again as in the previous section, one can extend the above discussion to the interacting case in the large bias regime upon substituting in the above results.

Finally, if we translate the above results back in the time domain we can write the effective action as

(58) |

This immediately shows that the retared component of the electron density fluctuation gives rise to a friction term in the photon effective action with damping coefficient , while the Keldysh term plays the role of a noise source with effective diffusion coefficient . To make explicit the connection with a classical stochastic Langevin dynamics we introduce the noise field through an Hubbard-Stratonovich decoupling Kamenev (2011) and then perform the integral over the quantum component which gives the classical dynamics

(59) |

with a white multiplicative noise. If now one couples this system to a coherent mode of the photonic transmission line playing the role of an explicit driving term, see Eq.(2), and with the force field deriving from the effective anharmonic potential in Eq. (45), the above dynamics describes a Duffing oscillator Kovacic and Brennan (2011) driven by noise Mallick and Marcq (2004), which has recently attracted interest in the context of circuit QED Larson and O’Dell (2013); André et al. (2012) and nanomechanical systems Zou et al. (2012).

### iv.3 Effective Temperature, Photon Occupation and Dissipated Electronic Power

We conclude this analysis by discussing in this section the nature of the environment generated by fermionic density fluctuations. In order to do so we start recalling that for a generic bosonic Green’s function in thermal equilibrium the following fluctuation dissipation relation (FDR) holds between the Keldysh and the retarded/advanced components

(60) |

where is the temperature of the system. In a generic out of equilibrium condition this relation does not hold, however it is still interesting to define an effective distribution function

(61) |

and to discuss its frequency dependence. In particular, a singular behaviour at low frequency

(62) |

it is often used to define an effective temperature for the out of equilibrium system, a concept that has emerged in many different contexts both classical and quantum Mitra and Millis (2005); Mitra et al. (2006); Foini et al. (2011); Dalla Torre et al. (2012); Ribeiro et al. (2013); Schiró and Mitra (2013). In the case of interest here the photon resonator is coupled to the effective bath induced by the electronic density fluctuations, encoded in the functions . We define therefore a photon effective distribution function as

(63) |

We plot in figure 7 the behaviour of this distribution function. It is easy to check analytically that, at zero bias , the bosonic self-energy satisfies the FDR, i.e.

(64) |

However, at finite bias a characteristic bosonic singularity emerges at low frequency, that allows to define a bias-dependent effective temperature . In particular, using the definition (62) and the results for the diffusion coefficient and the damping

(65) | |||

(66) |

we obtain an effective equilibrium result for the effective temperature

(67) |

that we plot in the inset of figure 7 as a function of the bias voltage. We see the small and large bias behaviours (compared to the electronic lifetime ) are characterized by two different power laws, at small bias when is almost set by the noise while at large voltage when the noise as we have seen saturates while the dissipation decays fast .

Quite interestingly, not only the bath distribution function but also the photon occupation is effectively thermal at large voltages (see right panel). In order to clarify this point we now compute the photon occupation as a function of the bias voltage. Let’s start considering the photon occupation, . This can be obtained from the Keldysh component of the photon Green’s function as

(68) |

In absence of electron-photon interaction, , we have

(69) |

from which we immediately get the thermal occupation .

In presence of a finite interaction the expression for the Keldysh component generally reads

(70) |

where the retarded/advanced component reads

(71) |

while are respectively the retarded/advanced/keldysh component of the photon self-energy, which we computed to lowest order in perturbation theory in as already discussed in previous sections.

It is now easy to see that for finite the first term in Eq. (70) gives a vanishing contribution to the photon occupation. Therefore, for the sake of computing the photon occupation we may write

(72) |

For weak coupling the retarded photon Green’s function is significantly different from zero only close to (neglecting small frequency renormalization effects), while the Keldysh component of the self-energy varies smoothly around . This suggests to write the retarded and advanced photon Green’s function as

(73) |

and to evaluate the above integral to get

(74) |

In figure 7 we plot the behaviour of the the photon number as a function of the voltage bias respectively for weak and large biases and for different photon frequencies . We see that the small bias regime features a clear threshold behavior, in accordance with previous results Mitra et al. (2004) namely the photon occupation sets in only as . This has to be understood in terms of energetic balance to excite a photon in the resonator by means of conduction electrons with energy above the Fermi sea. As opposite, the large bias regime shows a steep increase of the photon number, with a power-law behaviour that can be understood in term of an effective temperature description. Indeed as we have seen can be identified with the noise (or diffusion coefficient) while so that we have

(75) |

as we discussed in the previous section. Finally, we would like to clarify the relation between the photon occupation as a function of the bias voltage and the dissipated electronic power. In order to do this we consider the weak coupling case and proceed by perturbative evaluation of the photon and electron Green’s functions. It is tempting to compare the photon dissipated power with the electronic power due to the current flow, that generally reads . However this identification reveals a rather different behaviour both at small and at large voltage biases. Indeed for vanishing bias we expect

(76) |

while for large biases, when the current saturates we would expect

(77) |

These different power laws, that should survive also in presence of electron-photon coupling are difficult to reconcile with the result for the photon occupation. We finally remark that at higher frequencies departures from this effective equilibrium appear and the resulting semiclassical stochastic dynamics in this non-equilibrium environment Bradde and Biroli (2012) deserves further investigations.

## V Discussion

We can finally interpret the phase spectroscopy in light of the effective photon model valid at large bias. In this regime the non-linearity is only subleading and the main effects due to fermions are a many-body light frequency shift and a damping, so that one can write the photon propagator as

(78) |

resembling the spin susceptibility of an Ohmic spin-boson model for Le Hur (2012). This renormarlized frequency can be identified with a many-body Lamb shift. The crucial observation is now that, quite differently from conventional regimes, a large bias pushes the system into the regime , from which a phase of immediately follows on resonance.

We stress that this result is sensitive to external sources of decoherence for the photonic system, such as for example those provided by the auxiliary transmission lines. Indeed the presence of a finite dissipation would contribute to the photon damping with a finite value even in the large bias regime and then diminish the phase of the transmitted photon. As a result one should expect the phase to reach a maximum close to and then to cross over to zero for sufficiently large voltage biases, . Yet, for small values of the photon dissipation we expect a sizeable window of voltages where the phase stays close to , as confirmed by a direct evaluation of the transmission coefficient in presence of dissipation.

## Vi Conclusions

In this paper we discussed the physics of an hybrid quantum impurity model with interacting electrons and photons. We have shown that the phase of the transmitted photon at large voltage biases approaches the value of and we interpreted this result in light of an effective photon theory obtained integrating out the electronic degrees of freedom. It would be interesting to explore the connection between this result and the electronic (Friedel) phase shift at low bias, related to the Kondo effect. Finally we have discussed the robustness of the above result with respect to photon dissipation and shown that this hybrid system can be used as a simulator of the stochastic Duffing oscillator in relation with chaos.

## Vii Acknowledgment

It is a pleasure to thank M. Büttiker, A. Houck, T. Kontos, A. Browaeys, L. Henriet, C. Mora, J. Petta, P. Dutt, M. Filippone, L. Messio, O. Parcollet, Z. Ristivojevic, P. Simon, J. Gabelli, J. Esteve, and H. Tureci for interesting discussions and comments related to this work. This work has been funded by the LABEX PALM Paris-Saclay, grant Quantum Dyna and by DOE through DE-FG02-08ER46541 (KLH).

## References

- Raimond et al. (2001) J. M. Raimond, M. Brune, and S. Haroche, Rev. Mod. Phys. 73, 565 (2001).
- Cohen-Tannoudji et al. (1997) C. Cohen-Tannoudji, J. Dupont-Roc, and G. Grynberg, Photons and Atoms: Introduction to Quantum Electrodynamics (Wiley, 1997).
- Yamamoto et al. (2000) Y. Yamamoto, F. Tassone, and H. Cao, Semiconductor Cavity Quantum Electrodynamics (Springer, 2000), 1st ed., ISBN 3540675205.
- Schoelkopf and Girvin (2008) R. J. Schoelkopf and S. M. Girvin, Nature 451, 664 (2008).
- Akkermans and Montambaux (2007) E. Akkermans and G. Montambaux, Mesoscopic Physics of Electrons and Photons (Cambridge University Press, USA, 2007), 1st ed.
- Petta et al. (2005) J. R. Petta, A. C. Johnson, J. M. Taylor, E. A. Laird, A. Yacoby, M. D. Lukin, C. M. Marcus, M. P. Hanson, and A. C. Gossard, Science 309, 2180 (2005).
- Mekhov et al. (2007) I. B. Mekhov, C. Maschler, and H. Ritsch, Nature Physics 3, 319 (2007).
- Majer et al. (2007) J. Majer, J. M. Chow, J. M. Gambetta, J. Koch, B. R. Johnson, J. A. Schreier, L. Frunzio, D. I. Schuster, A. A. Houck, A. Wallraff, et al., Nature 449, 443 (2007).
- Suffczynski et al. (2009) J. Suffczynski, A. Dousse, K. Gauthron, A. Lemaitre, I. Sagnes, L. Lanco, J. Bloch, P. Voisin, and P. Senellart, Phys. Rev. Lett. 103, 027401 (2009).
- Delbecq et al. (2012) M. R. Delbecq, L. E. Bruhat, J. J. Viennot, S. Datta, A. Cottet, and T. Kontos, Nature Communication 4, 1400 (2012).
- Verdú et al. (2009) J. Verdú, H. Zoubi, C. Koller, J. Majer, H. Ritsch, and J. Schmiedmayer, Phys. Rev. Lett. 103, 043603 (2009).
- Parigi et al. (2012) V. Parigi, E. Bimbard, J. Stanojevic, A. J. Hilliard, F. Nogrette, R. Tualle-Brouri, A. Ourjoumtsev, and P. Grangier, Phys. Rev. Lett. 109, 233602 (2012).
- Ritsch et al. (2013) H. Ritsch, P. Domokos, F. Brennecke, and T. Esslinger, Rev. Mod. Phys. 85, 553 (2013).
- Marquardt and Girvin (2009) F. Marquardt and S. M. Girvin, ArXiv e-prints (2009), eprint 0905.0566.
- Pirkkalainen et al. (2013) J.-M. Pirkkalainen, S. U. Cho, J. Li, G. S. Paraoanu, P. J. Hakonen, and M. A. A. Sillanpää, Nature 494, 211 (2013).
- Kubo et al. (2010) Y. Kubo, F. R. Ong, P. Bertet, D. Vion, V. Jacques, D. Zheng, A. Dréau, J.-F. Roch, A. Auffeves, F. Jelezko, et al., Phys. Rev. Lett. 105, 140502 (2010).
- Dumeige et al. (2011) Y. Dumeige, R. AllÃ©aume, P. Grangier, F. Treussart, and J.-F. Roch, New Journal of Physics 13, 025015 (2011).
- Petersson et al. (2012) K. D. Petersson, L. W. McFaul, M. D. Schroer, M. Jung, J. M. Taylor, A. H. Houck, and J. R. Petta, Nature 490, 380 (2012).
- Frey et al. (2012) T. Frey, P. J. Leek, M. Beck, J. Faist, A. Wallraff, K. Ensslin, T. Ihn, and M. Büttiker, Phys. Rev. B 86, 115303 (2012).
- Jin et al. (2012) P.-Q. Jin, M. Marthaler, A. Shnirman, and G. Schön, Phys. Rev. Lett. 108, 190506 (2012).
- Jin et al. (2013) J. Jin, M. Marthaler, P.-Q. Jin, D. Golubev, and G. Schön, New Journal of Physics 15, 025044 (2013).
- Miguel-Sánchez et al. (2013) J. Miguel-Sánchez, A. Reinhard, E. Togan, T. Volz, A. Imamoglu, B. Besga, J. Reichel, and J. Estéve, New Journal of Physics 15, 045002 (2013).
- Bergenfeldt et al. (2013) C. Bergenfeldt, P. Samuelsson, B. Sothmann, C. Flindt, and M. Büttiker, ArXiv e-prints (2013), eprint 1307.4833.
- Viennot et al. (2013) J. J. Viennot, M. R. Delbecq, M. C. Dartiailh, A. Cottet, and T. Kontos, ArXiv e-prints (2013), eprint 1310.4363.
- Bocquillon et al. (2012) E. Bocquillon, F. D. Parmentier, C. Grenier, J.-M. Berroir, P. Degiovanni, D. C. Glattli, B. Plaçais, A. Cavanna, Y. Jin, and G. Fève, Phys. Rev. Lett. 108, 196803 (2012).
- Dubois et al. (2013) J. Dubois, T. Jullien, F. Portier, P. Roche, A. Cavanna, Y. Jin, W. Wegscheider, P. Roulleau, and D. C. Glattli, Nature (2013).
- Goldhaber-Gordon et al. (1998) D. Goldhaber-Gordon, J. Göres, M. A. Kastner, H. Shtrikman, D. Mahalu, and U. Meirav, Phys. Rev. Lett. 81, 5225 (1998).
- Cronenwett et al. (1998) S. M. Cronenwett, T. H. Oosterkamp, and L. P. Kouwenhoven, Science 281, 540 (1998).
- Türeci et al. (2011) H. E. Türeci, M. Hanl, M. Claassen, A. Weichselbaum, T. Hecht, B. Braunecker, A. Govorov, L. Glazman, A. Imamoglu, and J. von Delft, Phys. Rev. Lett. 106, 107402 (2011).
- Latta et al. (2011) C. Latta, F. Haupt, M. Hanl, A. Weichselbaum, M. Claassen, W. Wuester, P. Fallahi, S. Faelt, L. Glazman, J. Von Delft, et al., Nature 474, 627 (2011).
- Mebrahtu et al. (2012) H. T. Mebrahtu, I. V. Borzenets, D. E. Liu, X. Zheng, Y. V. Bomze, A. I. Smirnov, H. U. Baranger, and G. Finkelstein, Nature 488 (2012).
- Jezouin et al. (2012) S. Jezouin, M. Albert, F. D. Parmentier, A. Anthore, U. Gennser, A. Cavanna, I. Safi, and F. Pierre, Nature Communications 4 (2012).
- Hofheinz et al. (2011) M. Hofheinz, F. Portier, Q. Baudouin, P. Joyez, D. Vion, P. Bertet, P. Roche, and D. Esteve, ArXiv e-prints (2011), eprint 1102.0131.
- Ingold et al. (1992) G. L. Ingold, D. V. Averin, and Y. V. Nazarov, in Single Charge Tunneling, NATO ASI Series B 294 (Plenum, New York, 1992).
- Safi and Saleur (2004) I. Safi and H. Saleur, Phys. Rev. Lett. 93, 126602 (2004).
- Le Hur (2004) K. Le Hur, Phys. Rev. Lett. 92, 196804 (2004).
- Le Hur and Li (2005) K. Le Hur and M.-R. Li, Phys. Rev. B 72, 073305 (2005).
- Le Hur et al. (2007) K. Le Hur, P. Doucet-Beaupré, and W. Hofstetter, Phys. Rev. Lett. 99, 126801 (2007).
- Chung et al. (2013) C.-H. Chung, K. Le Hur, G. Finkelstein, M. Vojta, and P. Wölfle, Phys. Rev. B 87, 245310 (2013).
- LeClair et al. (1997) A. LeClair, F. Lesage, S. Lukyanov, and H. Saleur, Physics Letters A 235, 203 (1997), ISSN 0375-9601.
- Le Hur (2012) K. Le Hur, Phys. Rev. B 85, 140506 (2012).
- Orth et al. (2013) P. P. Orth, A. Imambekov, and K. Le Hur, Phys. Rev. B 87, 014305 (2013).
- Henriet et al. (2014) L. Henriet, Z. Ristivojevic, P. P. Orth, and K. Le Hur, ArXiv e-prints (2014), eprint 1401.4558.
- Goldstein et al. (2013) M. Goldstein, M. H. Devoret, M. Houzet, and L. I. Glazman, Phys. Rev. Lett. 110, 017002 (2013).
- Bera et al. (2013) S. Bera, S. Florens, H. Baranger, N. Roch, A. Nazir, and A. Chin, ArXiv e-prints (2013), eprint 1301.7430.
- Sbierski et al. (2013) B. Sbierski, M. Hanl, A. Weichselbaum, H. E. Türeci, M. Goldstein, L. I. Glazman, J. von Delft, and A. Imamoglu, Phys. Rev. Lett. 111, 157402 (2013).
- Beugnon et al. (2006) J. Beugnon, M. P. A. Jones, J. Dingjan, B. Darquié, G. Messin, A. Browaeys, and P. Grangier, Nature 440, 779 (2006).
- Delbecq et al. (2011) M. R. Delbecq, V. Schmitt, F. D. Parmentier, N. Roch, J. J. Viennot, G. Fève, B. Huard, C. Mora, A. Cottet, and T. Kontos, Phys. Rev. Lett. 107, 256804 (2011).
- Holstein (1959) T. Holstein, Annals of Physics 8, 325 (1959), ISSN 0003-4916.
- Anderson (1975) P. W. Anderson, Phys. Rev. Lett. 34, 953 (1975).
- Mitra et al. (2004) A. Mitra, I. Aleiner, and A. J. Millis, Phys. Rev. B 69, 245302 (2004).
- Cornaglia et al. (2004) P. S. Cornaglia, H. Ness, and D. R. Grempel, Phys. Rev. Lett. 93, 147201 (2004).
- Mravlje and Ramsak (2008) J. Mravlje and A. Ramsak, Phys. Rev. B 78, 235416 (2008).
- Vinkler et al. (2012) Y. Vinkler, A. Schiller, and N. Andrei, Phys. Rev. B 85, 035411 (2012).
- Laakso et al. (2014) M. A. Laakso, D. M. Kennes, S. G. Jakobs, and V. Meden, New Journal of Physics 16, 023007 (2014).
- Clerk et al. (2010) A. A. Clerk, M. H. Devoret, S. M. Girvin, F. Marquardt, and R. J. Schoelkopf, Rev. Mod. Phys. 82, 1155 (2010).
- Filippone et al. (2011) M. Filippone, K. Le Hur, and C. Mora, Phys. Rev. Lett. 107, 176601 (2011).
- Werner et al. (2009) P. Werner, T. Oka, and A. J. Millis, Phys. Rev. B 79, 035320 (2009).
- Schiró and Fabrizio (2009) M. Schiró and M. Fabrizio, Phys. Rev. B 79, 153302 (2009).
- Dutt et al. (2011) P. Dutt, J. Koch, J. E. Han, and K. Le Hur, Annals of Physics 326, 2693 (2011).
- Coleman (1987) P. Coleman, Phys. Rev. B 35, 5072 (1987).
- Houghton et al. (1987) A. Houghton, N. Read, and H. Won, Phys. Rev. B 35, 5123 (1987).
- Rosch et al. (2001) A. Rosch, J. Kroha, and P. Wölfle, Phys. Rev. Lett. 87, 156802 (2001).
- Ratiani and Mitra (2009) Z. Ratiani and A. Mitra, Phys. Rev. B 79, 245111 (2009).
- Kovacic and Brennan (2011) I. Kovacic and M. J. Brennan, The Duffing Equation: Non Linear Oscillators and their Behaviour (John Wiley and Sons, 2011).
- Kamenev (2011) A. Kamenev, Field Theory of Non-Equilibrium Systems (Cambridge University Press, 2011).
- Mozyrsky et al. (2006) D. Mozyrsky, M. B. Hastings, and I. Martin, Phys. Rev. B 73, 035104 (2006).
- Nocera et al. (2012) A. Nocera, C. A. Perroni, V. Marigliano Ramaglia, and V. Cataudella, Phys. Rev. B 86, 035420 (2012).
- Buttiker et al. (1993) M. Büttiker, H. Thomas, and A. Pretre, Physics Letters A 180, 364 (1993), ISSN 0375-9601.
- Gabelli et al. (2006) J. Gabelli, G. Fève, J.-M. Berroir, B. Plaçais, A. Cavanna, B. Etienne, Y. Jin, and D. C. Glattli, Science 313, 499 (2006).
- Fève et al. (2007) G. Fève, A. MahÃ©, J.-M. Berroir, T. Kontos, B. Plaçais, D. C. Glattli, A. Cavanna, B. Etienne, and Y. Jin, Science 316, 1169 (2007).
- Mora and Le Hur (2010) C. Mora and K. Le Hur, Nature Physics 6, 697 (2010).
- Gabelli et al. (2012) J. Gabelli, G. Fève, J.-M. Berroir, and B. Plaçais, Reports on Progress in Physics 75, 126504 (2012).
- Nigg et al. (2006) S. E. Nigg, R. López, and M. Büttiker, Phys. Rev. Lett. 97, 206804 (2006).
- BÃ¼ttiker and Nigg (2009) M. Büttiker and S. E. Nigg, The European Physical Journal Special Topics 172, 247 (2009), ISSN 1951-6355.
- Etzioni et al. (2011) Y. Etzioni, B. Horovitz, and P. Le Doussal, Phys. Rev. Lett. 106, 166803 (2011).
- Hamamoto et al. (2010) Y. Hamamoto, T. Jonckheere, T. Kato, and T. Martin, Phys. Rev. B 81, 153305 (2010).
- Filippone et al. (2013) M. Filippone, K. Le Hur, and C. Mora, Phys. Rev. B 88, 045302 (2013).
- Dutt et al. (2013) P. Dutt, T. L. Schmidt, C. Mora, and K. Le Hur, Phys. Rev. B 87, 155134 (2013).
- Mallick and Marcq (2004) K. Mallick and P. Marcq, The European Physical Journal B - Condensed Matter and Complex Systems 38, 99 (2004), ISSN 1434-6028.
- Larson and O’Dell (2013) J. Larson and D. H. J. O’Dell, ArXiv e-prints (2013), eprint 1304.3582.
- André et al. (2012) S. André, L. Guo, V. Peano, M. Marthaler, and G. Schön, Phys. Rev. A 85, 053825 (2012).
- Zou et al. (2012) J. Zou, S. Buvaev, M. Dykman, and H. B. Chan, Phys. Rev. B 86, 155420 (2012).
- Mitra and Millis (2005) A. Mitra and A. J. Millis, Phys. Rev. B 72, 121102 (2005).
- Mitra et al. (2006) A. Mitra, S. Takei, Y. B. Kim, and A. J. Millis, Phys. Rev. Lett. 97, 236808 (2006).
- Foini et al. (2011) L. Foini, L. F. Cugliandolo, and A. Gambassi, Phys. Rev. B 84, 212404 (2011).
- Dalla Torre et al. (2012) E. G. Dalla Torre, E. Demler, T. Giamarchi, and E. Altman, Phys. Rev. B 85, 184302 (2012).
- Ribeiro et al. (2013) P. Ribeiro, Q. Si, and S. Kirchner, EPL (Europhysics Letters) 102, 50001 (2013).
- Schiró and Mitra (2013) M. Schiró and A. Mitra, ArXiv e-prints (2013), eprint 1310.0117.
- Bradde and Biroli (2012) S. Bradde and G. Biroli, ArXiv e-prints (2012), eprint 1204.6027.