Cavity squeezing by a quantum conductor

Hybrid architectures integrating mesoscopic electronic conductors with resonant microwave cavities have a great potential for investigating unexplored regimes of electron-photon coupling. In this context, producing nonclassical squeezed light is a key step towards quantum communication with scalable solid-state devices. Here we show that parametric driving of the electronic conductor induces a squeezed steady state in the cavity. We find that squeezing properties of the cavity are essentially determined by the electronic noise correlators of the quantum conductor. In the case of a tunnel junction, we predict that squeezing is optimized by applying a time-periodic series of quantized $\delta-$peaks in the bias voltage. For an asymmetric quantum dot, we show that a sharp Leviton pulse is able to achieve perfect cavity squeezing.


I. INTRODUCTION
Squeezed states of light [1] exhibit reduced noise below the vacuum level in one of their quadrature and amplification in the other quadrature. Their realization is a key step in the development of quantum communication. They are important tools for continuous variable quantum information protocols [2,3] where they serve as building blocks for generating non-classical states. Their enhanced sensitivity can also be used for quantum non-demolition measurements of position and force [4]. Easily produced in optical systems, squeezed states have been observed more recently in circuit quantum electrodynamics at microwave frequencies [5], either as single-mode [6,7], two-mode squeezing [8,9] or as Einstein-Podolsky-Rosen states [10].
The parametric driving used so far in experiments is limited to a half-squeezed quadrature for a cavity mode because of the inevitable coupling to the external vacuum fluctuations [11]. This limit however does not apply to dissipative squeezing, in which one steers the environment to stabilize the cavity into a non-classical state. In this case, perfect squeezing can be achieved, at least in principle, with minimum uncertainty [12,13].
A recent development in the field of superconducting quantum circuits is the realization of hybrid systems in which a quantum conductor is coupled to a microwave resonator. These systems offer an appealing platform for investigating fundamental matterlight interactions with an experimental control on both the electronic and photonic parts [14][15][16][17][18]. Experiments have been realized with metallic tunnel junctions connected to a resonating line [19], or with quantum dots, realized in carbon nanotubes, nanowires or two-dimensional electron gases [20][21][22][23][24][25][26] embedded in high-finesse coplanar cavities. The interplay of electron transport and emission of photons can lead to an electronic-induced lasing state in the cavity [27][28][29][30], and more generally produce bunched or antibunched photons [31][32][33][34][35], and nonclassicality in the light emitted by a quantum conductor [36,37]. Squeezed light emitted by a tunnel junction was recently demonstrated experimentally in the absence of a cavity [38][39][40].
In this paper, we describe dissipative squeezing of a cavity mode coupled to an ac driven electronic reservoir. The system, depicted in Fig. 1, is a quantum conductor coupled to a microwave resonator. In addition, a classical bias voltage, with an oscillating part at twice the resonator frequency, is applied to the conductor. The conductor plays the role of a nonlinear environment: photons from the ac modulation are broken into pairs and transmitted to the cavity, thereby producing squeezing. We show that the amount of cavity squeezing is determined solely by current noise fluctuations in the conductor. We focus our attention on ac excitations which optimizes squeezing. In the case of a tunnel junction, we find that the best solution consists of periodic and quantized voltage peaks occurring in phase with the compressed quadrature. We also discuss how squeezing improves with the number of harmonics in the ac signal. For a quantum dot, we identify the conditions of optimum squeezing: asymmetric coupling to the leads, a narrow single-level resonance, a far-detuned single-level energy and a dc bias voltage matching the resonator frequency. In addition, we show that a Leviton pulse results in a vacuum squeezed state with minimal uncertainty. Perfect squeezing is approached by narrowing the width of the voltage pulses in the Leviton, again in phase with the compressed quadrature. The center line is coupled to external electromagnetic modes (ports) via capacitive gaps delimiting the cavity. The output field of the cavity can be measured in these ports. A dc voltage can be applied to the center line at a voltage node in order to preserve the high-Q factor of the cavity [22,41]. (b) Focus on the galvanic coupling between the cavity and a quantum conductor. The electric field (arrows) is spatially constant at the scale of the conductor. The conductor is lead-contacted on the upper side to the center line and on the lower side to the ground plane. (c) The quantum conductor can be a tunnel junction or a quantum dot. The upper lead is connected to the center line where the voltage potential is the sum of dc [V ] and ac [Vac(t)] parts, as well as the cavity quantum field ∝ a + a † .

II. TUNNEL JUNCTION
We start by considering the case of a tunnel junction for the quantum conductor. The voltage bias across the junction is the sum of a driven classical part applied on the upper lead, and a quantum part associated to the cavity field. Using a convenient choice of electromagnetic gauge, it simply dresses the electron tunneling operators where the operator T transfers one electron from the upper to the lower lead while T † does the opposite [42,43].φ(t) is the total semi-classical phase accumulated during a tunneling event, decomposed asφ(t) = (eV / )t+φ ac (t)+i(g/ω 0 )(â † −â), with the ac phase φ ac (t) = (e/ ) t dt V ac (t ). The last term is the quantum part and g measures the junction-cavity coupling strength.
The cavity Hamiltonian is reduced for simplicity to a single mode, H cav = ω 0â †â with the cavity annihilation operatorâ. Eq. (1) contains both the excitation of the cavity state by electron tunneling events and photo-assisted transport phenomena triggered by the ac modulation [44,45]. H T (t) is thus responsible for an exchange of energy between three sub-systems: the cavity, the ac classical field and lead (free) electrons.
We assume weak junction-cavity coupling and therefore neglect the backaction-induced change in electron tunneling resulting from the cavity. We thus set g = 0 to examine the current fluctuations of the tunnel junction and latter reinstate a finite g when considering the dynamics of the cavity fieldâ. In the presence of ac voltage modulation, the photo-assisted noise properties of the tunnel junction are characterized by the correlator (at g = 0) [46] whereÎ(ω) is the Fourier transform of the Heisenberg current operatorÎ(t) of the junction. The n = 0 term gives the stationary part of the noise, meaning that, when Fourier transformed with times t 1 and t 2 , it depends only on the time difference t 1 − t 2 and not on the mean timet = (t 1 + t 2 )/2. S 0 (ω) is called the absorption (emission) noise of the tunnel junction for ω > 0 (ω < 0) [47]. It governs the rate of energy transfer between the junction and its environment, here the cavity, via single photons of energy |ω|. As shown in the Supplementary Note 1, the different noise terms can be calculated to leading order in the tunneling strength in terms of the (unsymmetrized) equilibrium Johnson-Nyquist noise of the tunnel junctionS(ω) where R T is the junction dc resistance. The coefficients c m are defined by the Fourier expansion of the ac phase φ ac . ϕ is the overall phase of the ac signal. For the lead connected to the voltage V ac (t) = ( /e)φ ac (t), c m gives the probability amplitude for an electron to absorb m energy quanta from the classical ac field when m > 0. m < 0 describes correspondingly photon emission to the ac field [48]. In the particular case of a sinusoidal excitation, V ac (t) = V 1 cos(2ω 0 t+ϕ), these coefficients are written in terms of Bessel functions The energy of absorbed photons, 2m ω 0 , can be either used in exciting energetic electron-hole pairs in the conductor, or transferred to the cavity. The non-stationary noise terms S n =0 in Eq. (2) oscillate with the mean timet. They do not conserve energy and can provide n quanta of energy 2 ω 0 (or absorb if n < 0) to the cavity. We will return below to the physical significance of these terms when analysing the cavity stationary state and squeezing effects. Now that we have detailed the possible transfers of energy between the ac-excited tunnel junction and its environment (the cavity), we study the cavity evolution under the dissipative influence of the electrons. Assuming weak junction-cavity coupling, we expand H T to first order in g/ω 0 , H T = H 0 T + iλÎ(â † −â) with the coupling constant λ = g/(eω 0 ). The cavity evolution is described by a Heisenberg-Langevin equationȧ justified, either by an input-output calculation [49] detailed in the Supplementary Note 2, or by a Keldysh path integral formulation, discussed in appendix B, assessing the Gaussian character of current and cavity field fluctuations. The tunnel junction currentÎ(t) in Eq. (6) plays the role of a quantum noise term, with fluctuations characterized by the correlator of Eq. (2). In deriving this equation, we have neglected the intrinsic (bare) damping of the cavity κ 0 , assuming that the cavity dissipation caused by the electrons dominates. The corresponding damping rate κ = λ 2 [S 0 (ω 0 ) − S 0 (−ω 0 )] balances absorption and emission noises, since absorption (emission) noise corresponds to photon loss (gain) from the cavity. Our calculation is also based on the use of the rotating-wave approximation where rapidly oscillating terms are averaged to zero. The validity of this approximation is controlled by the smallness of κ/ω 0 and is consistent with our first order conductor-cavity coupling and with the absence of cavity backaction. The first order differential Eq. (6) can be solved straightforwardly in time or frequency space, and yields the steady-state correlation functions for the cavity fieldâ. We find for the anomalous correlator, using that κ ω 0 , This result can be given a physical interpretation: S 1 (ω 0 ) describes the coherent emission by the junction of a quantum of energy 2 ω 0 to the cavity. This energy quantum breaks into a pair of cavity photons thereby contributing to the â 2 correlator. This effect is limited in the denominator by the rate at which cavity photons are absorbed by the tunnel junction. In the same way, the number of cavity photons â †â = λ 2 S 0 (−ω 0 )/κ is unsurprisingly governed by the electronic emission noise. We now investigate field squeezing more precisely and introduce the two quadraturesX 1 = i(â † e −iϕ/2 −âe iϕ/2 ) and X 2 =â † e −iϕ/2 +âe iϕ/2 , where we use the same phase ϕ as in Eq. (4). Their variance is readily obtained where we introduced the evenS e and oddS o parts of the Johnson-Nyquist noiseS. The two quadrature fluctuations are thus sensitive to the electronic temperature T , the dc bias voltage V and the pulse shape of the ac signal. The squeezing mechanism is optimized by taking the limit of vanishing temperature and by setting eV = ω 0 . We first consider a single-tone driving of the tunnel junction, V ac (t) = V 1 cos(2ω 0 t + ϕ), the experimentally most accessible situation. The photo-assisted coefficients c n are then given by Bessel functions, as detailed in appendix A. A numerical minimization of ∆X 2 1 in Eq. (8) gives ∆X 2 1 = 0.618 for eV 1 = 0.706 × 2 ω 0 , with ∆X 2 2 = 1.864. This optimal squeezing value coincides exactly with the squeezing of the emitted light predicted and measured in Refs. [37,38], see also the more recent Ref. [50]. This is explained by noting that the zero-temperature cavity damping κ/λ 2 = 2 ω 0 /R T (the denominator in Eq. (8)) is constant for a tunnel junction, regardless of the bias voltage shape. This independence no longer holds at finite temperature or in the case of a conductor with a non-linear I-V characteristic. We turn to an ac modulation with the same fundamental frequency 2ω 0 but including higher harmonics [51]. Fig. 2 shows the improvement in squeezing ∆X 2 1 by adding more and more harmonics while ∆X 2 2 is further amplified. Considering a general periodic signal, we find analytically, as shown in appendix C, that the minimum value ∆X 2 1 = 4/π 2 = 0.405 is reached when c n = (1/π)(n+1/2) −1 for n ∈ Z, in agreement with a numerical minimization. The corresponding ac phase across the junction is a periodic piecewise linear function φ ac (t) = π/2−ω 0 t for t ∈]0, π/ω 0 [, with a jump discontinuity of π at t = 0 and multiples of π/ω 0 . Adding the dc voltage V, we find for the optimal voltage applied to the junction a series of δ-peaks, It is useful to give an intuitive classical picture for this squeezing optimization: the bias potential V opt (t) acts on the conductor specifically at times where the amplitude of the squeezed quadrature is maximum and the other quadrature vanishes. This is in fact a strong perturbation, the emission and absorption noises are infinite and the second variance ∆X 2 2 = +∞. Despite its linear current response, the tunnel junction is able to squeeze the cavity state down to 40% of the zero-point level. This is because squeezing is not governed by the current itself but by current fluctuations, and the noise of a tunnel junction is only a piecewise linear function enabling rectification [38]. Eq. (8) nevertheless suggests that better squeezing can be achieved by using a genuine non-linear system.

III. ASYMMETRIC QUANTUM DOT
We consider a quantum dot for the conductor embedded in the cavity. The situation where the dot is symmetrically coupled to the two leads, discussed in the Supplementary Note 3, is not optimal for squeezing. We thus focus on the asymmetric case where the upper lead is more weakly coupled to the quantum dot than the lower lead. In this case, the voltage drop from the central strip to the ground mainly takes place at the upper dot-lead tunnel contact. The first order cavity-conductor coupling is then of the form iλÎ U (â † −â) whereÎ U denotes the electrical current of electrons incoming from the upper lead. The coupling to the lower lead current is neglected.
In practice, for quantum dot geometries, it may be important to also take into account the coupling of electronic transport to phonons. Following Ref. [52], it would lead to tunneling processes involving the excitation of phonon-photon pairs, degrading the quality of squeezing. Such study is nevertheless beyond the scope of this work and we neglect electron-phonon coupling in what follows.
The analysis developed above for the tunnel junction can be essentially carried over to the quantum dot, withÎ U replacingÎ in the Heisenberg-Langevin Eq. (6). The noise properties of the quantum dot are derived using scattering theory as discussed in appendix D. We retrieve noise factors similar to Eq. (3), here involving two different equilibrium noise terms with the Fermi function f (ε) = (1+e ε/k B T ) −1 . The broadening Γ of the dot single energy level, denoted ε d , can be decomposed according to its coupling to the upper and lower leads Γ = Γ U + Γ L with Γ U Γ L . The Lorentzian form in the integrand of Eq. (11) describes the Breit-Wigner resonance for transmitting electrons through the dot [4].S + (ω) describes electron-hole excitations with energy ω, where the electron, with energy ε, tunnels from the lower to the upper lead and has to meet the resonance condition of the dot single-level.S − is the same but with hole tunneling. Proceeding with the calculation of cavity properties based on the Heisenberg-Langevin equation of motion, we retrieve the quadrature variances of Eq. (8) if we define the even/odd parts as 2S e/o (ω) =S + (ω) ±S − (−ω).
We have studied numerically the minimization of the variance ∆X 2 1 . Quite generally, squeezing optimization requires zero temperature and the four-wave mixing condition eV = ω 0 [38], as well as a sharp resonance, Γ ω 0 . In this regime, for ε d > 0,S − becomes negligible andS + in Eq. (8) is either constant for n > n th = ε d /(2 ω 0 ) − 1, or vanishingly small below this threshold.S e andS o simplify in Eq. (8) and the summation involves only values of n above the threshold n th .
In the case of a single-frequency ac modulation, Fig. 3 shows that the squeezed variance ∆X 2 1 displays a series of local minima, where the values of the minima decrease with ε d . Large single-level energy ε d however implies stronger power in the ac excitation signal in order to meet the Breit-Wigner resonance condition. In practice, this requires an ac signal amplitude V 1 close to ε d such that electrons can tunnel through the dot. Perfect squeezing ∆X 2 1 = 0 is reached at very large ε d but only with a weak power law of coefficient −1/3.
Alternatively, a vacuum squeezed state can be reached in the cavity by using a Leviton ac signal with the fundamental frequency 2ω 0 . Leviton pulses were originally [54] proposed as voltage excitations designed to transfer a finite number of electrons through a coherent conductor with minimal noise, in analogy with coherent states minimizing quantum-mechanical uncertainty. They consist of sums of Lorentzian pulses with unit flux each. A unit flux represents the attempt to transmit a single electron. Mathematical details are briefly reviewed in appendix A. A Leviton can be periodized by having an infinite train of evenly spaced Lorentzian pulses [48,55]. Leviton pulses have recently been synthesized and used to perform Hong-Ou-Mandel electronic experiments [56] and electron quantum tomography [57].
The use of a Leviton pulse for squeezing is natural. A Leviton train with periodic short pulses addresses specifically one quadrature (the one to be squeezed), while producing a minimal disturbance (noise) on the quantum conductor. Taking the limit of a very sharp resonance Γ/ω 0 → 0 and the bias voltage eV = ω 0 , one obtains for the squeezed variance for 0 < ε d < 2 ω 0 , and a minimal Heisenberg uncertainty ∆X 1 ∆X 2 = 1 reflecting the minimized perturbation by the Leviton compared to other types of ac excitation. a = e −2ω0τ is a parameter related to the width τ of each Lorentzian pulse. A Leviton pulse is thus able to produce in optimal conditions an ideal squeezed state with arbitrary compression. As shown in Fig. 4, squeezing saturates at finite damping Γ. For a reasonable damping rate Γ/ω 0 = 0.1, we find that the squeezed variance can still be reduced down to ∆X 2 1 = 0.075 for a = 0.91. A Leviton excitation with parameter a close to 1 exhibits sharp voltage pulses of width τ (1 − a)/(2ω 0 ). Similarly to the tunnel junction, it corresponds to concentrate shorttime pulses of high voltage when the amplitude of the squeezed quadrature is maximum while the other quadrature vanishes.

IV. CAVITY READOUT
Cavity squeezing could be evidenced by an in-situ measurement using a qubit and its anisotropic radiative decay [58]. It is also possible to demonstrate squeezing in the cavity by measuring the output field. So far, we have neglected in our discussion the coupling of the cavity to the external electromagnetic modes, see Fig. 1. They introduce an additional damping rate κ 0 . As detailed in the Supplementary Note 2, this external damping adds two terms in the Heisenberg-Langevin Eq. (6), which now readsȧ b in (t) is the input field, it describes the quantum state of incident photons on the cavity. In the absence of an input drive, it corresponds to vacuum noise with b in (t)b † in (t ) = δ(t − t ) at zero temperature. The output field, describing photons escaping the cavity, is given byb out =b in + √ κ 0â and obtained by solving Eq. (13).
Squeezing properties of the cavity are revealed by homodyne detection of the output field [5,6], mixing it with a local oscillator with the cavity frequency ω 0 and phase θ. The power spectrum S D (ω) [1] of the homodyne detector signalÎ D (t) ∝ b out (t)e i(ω0t+θ) + h.c. exhibits a deep at zero frequency and θ = π/2 + ϕ/2, where the squeezing effect is most visible, with assuming weak external damping κ 0 κ. S 0 D denotes the vacuum level, measured in the absence of ac excitation. A value smaller than one in this equation indicates a squeezed output field. Interestingly, output squeezing is weak in this limit due to pollution by the input vacuum noise. As discussed in appendix E, the output field gets more squeezed as κ 0 increases, the best squeezing being obtained for equal external and electronic dampings. This is however at the price of weaker squeezing in the cavity. Perfect squeezing could even be reached in the output field, limiting in this case the cavity state to half-squeezing.

V. CONCLUSION AND OUTLOOK
We studied the squeezing generated in a resonant cavity by coupling it to a mesoscopic conductor under parametric excitation. We showed that the quality of squeezing can be improved by enhancing nonlinearities in the conductor and by concentrating the voltage excitation pulses at instants where the squeezed quadrature amplitude reaches its maximum. In optimal conditions, perfect squeezing can even be achieved. We remark that nonlinearities could also be enhanced in a tunnel junction by increasing the coupling to the cavity and the associated dynamical Coulomb blockade [60]. More generally, our results can be easily extended to other quantum conductors for which the photo-assisted noise spectra are known or can be computed. They also suggest the possibility to engineer squeezed light for quantum information using electronic quantum conductors. time-dependent voltage is given by a sum of quantized Lorentzian pulses of common width τ . The ac phase is obtained by time-integration, φ ac (t) = (e/ ) t dt V ac (t ). It is conveniently written using the cyclic variable z = e 2iω0t , namely after a few algebraic manipulation [55] where 0 < a = e −2ω0τ < 1 is related to the pulses width. The conventional Leviton pulse, shaped to minimize zero-frequency shot noise, has the dc level eV = 2 ω 0 , cancelling the first term in Eq. (A1). In the main text, the dc bias chosen to optimize squeezing is eV = ω 0 . It can be understood by noting that the goal here is to minimize finite frequency noise correlators at the frequency ω 0 . Expanding Eq. (A2) in powers of z, we obtain c −1 = −a, c n≥0 = a n (1 − a 2 ): there is a strong imbalance between absorption and emission of photons. a = 1 is no longer an ac modulation but corresponds to a shift of the dc voltage by −2( /e)ω 0 .

Appendix B: Keldysh formulation and Heisenberg-Langevin equation
The out-of-equilibrium physics of our system is conveniently described within the Keldysh path-integral formalism [61], enabling a systematic integration of the electronic degrees of freedom and yielding an effective action for the photons [62]. The action obtained, Eqs. (B3) and (B4), can be shown to be equivalent to the Heisenberg-Langevin equation of motion Eq. (6) used in the main text, in the limit of small κ/ω 0 1. For a direct derivation of the Heisenberg-Langevin equation in the spirit of input-output theory, see Supplementary Note 2.
The partition function (hereafter = 1) involves an integration over complex-valued fields, a,a * for photons and c,c * for electrons. S e is the action for the isolated conductor, S cav the photonic part corresponding to the Hamiltonian H cav = ω 0â †â and the electron-photon coupling is to first order where η denotes the Keldysh time branch. I is the quantum conductor current written with complex-valued fields. The electronic part can be rigorously integrated using the cumulant expansion e iSep e = e i Sep e−(1/2) δS 2 ep e +... , where we use the notations δS ep = S ep − S ep and D[c, c * ]e iSe A = A e . We assume S ep e = 0, a finite S ep e can be absorbed by a small shift ∝ g in the cavity fields. To summarize, we square S ep , take its quantum average restricted to electronic degrees of freedom, and thus obtain a self-energy kernel for the photons involving current noise correlators. For clarity, we switch to classical/quantum variables, a cl/q = (a + ± a − )/ √ 2 and write the action in frequency space in order to take advantage of the current noise correlators given in Eq. (2). Summing S cav and (i/2) S 2 ep e , we find the effective action where we used the notation ω = +∞ −∞ dω 2π , and a cl (t) = ω a cl (ω)e −iωt . The retarded photon Green's function G −1 R (ω) = ω − ω 0 + iκ/2 has a pole shifted by half the damping rate κ (see main text). The quantum-quantum self-energy part is Σ K = −iλ 2 [S 0 (ω 0 ) + S 0 (−ω 0 )]. The effective action S eff includes also an anomalous term, responsible for state squeezing, Note that the real part of the self-energy induces in general a cavity pull which has been absorbed into a redefinition of ω 0 . Computing this frequency shift consistently requires the second order term in the expansion of H T in powers of g/ω 0 , namely where V T = T + T † . The cavity frequency shift then vanishes for a tunnel junction. The derivation of Eqs. (B3) and (B4) relies on the rotating-wave approximation, valid for κ ω 0 , where the fields a cl/q (ω), a * cl/q (ω) take significant values only for ω ω 0 . In principle, the anomalous part of the action S a also contains terms with a cl (ω)a q (2ω 0 − ω), corresponding to the effect of non-stationary noise terms on the damping. Those terms are found to be proportional to S 1 (ω) − S 1 (2ω 0 − ω) and thus vanish for ω ω 0 , with the small parameter κ/ω 0 .
Finally, we discuss the connexion between the quadratic action of Eq. (B3) and the Heisenberg-Langevin equation (6). Quite generally, it is known that current fluctuations in a tunnel junction, or a quantum conductor, are not gaussian. Nevertheless, computing the non-gaussian current contributions to the statistics of photons, one finds that they are small compared to the dominant Wick-like contractions among the current operators. This is true in the limit of weak damping, κ ω 0 , where cavity correlation functions only involve current operatorsÎ at frequencies ±ω 0 . For example, the fourth-order cavity field correlator gives, implying photon gaussian statistics except for the specific case of eV close to ω 0 with no ac excitation. A related discussion can be found in Refs. [31,32]. The reason is that an electron-hole excitation with energy ω 0 created by a current operator I ω0 must be destroyed by another single current operator, the phase-space for alternative processes -where electron and hole are annihilated by two distinct current operators -being negligible for weak damping κ. This argument pertains to higher-order correlation functions such that, for the purpose of photon statistics, it is legitimate to keep only the gaussian part of electronic current fluctuations. The resulting cavity field statistics are obviously gaussian. The above cumulant expansion can be rigorously stopped after the second order and the gaussian action in Eq. (B3) becomes exact as long as κ is negligible with respect to ω 0 . Computing second-order cavity correlations functions, with different ordering ofâ andâ † , we find coinciding results for the gaussian action and the Heisenberg-Langevin evolution. This completes the proof of the equivalence of the two formulations. This comparison differs from the standard derivation of a classical Langevin equation using the Keldysh action [63], in which case information about operator ordering is lost.

Appendix C: Optimized squeezing for a tunnel junction
In this section, we set = 1 for simplicity. We focus on the zero-temperature case, relevant to maximize the cavity state compression. In this case,S o [eV + (2n + 1)ω 0 ] = 1 R T eV + (2n + 1)ω 0 (C1) and the denominator in Eq. (8) of the main text simplifies to n∈Z |c n | 2 − |c n+1 | 2 eV + (2n + 1)ω 0 = 2ω 0 , regardless of V and the c n coefficients. At the optimal dc voltage eV = ω 0 , the squeezed variance takes the simple form which we still need to minimize with respect to the distribution of Fourier coefficients c n , We now prove that the pulse shape of Eq. (9) in the main text, corresponding to the piecewise linear phase φ ac,opt (t) = π/2−ω 0 t, extremizes the variance ∆X 2 1 . We first differentiate Eq. (C4) to obtain which gives −(ω 0 /π)e −i(2n+1)ω0t when evaluated at φ ac,opt (t). Using this result, we can proceed with the derivative of ∆X 2 1 with respect to an arbitrary form of φ ac (t) We evaluate this derivative with φ ac,opt (t) and its coefficients and obtain which completes the proof. Inserting the coefficients Eq. (C7) into the quadrature variance Eq. (C3), we find (C9) We also checked numerically that φ ac,opt (t) reaches the global minimum of ∆X 2 1 .

Appendix D: Landauer-Büttiker calculation of the noise
The noise properties of the quantum dot are derived using the scattering, or Landauer-Büttiker, formalism [5]. The current operator is expanded over the basis of one-particle scattering states originating from both leads. The general case is reviewed in the Supplementary Note 4, we focus here on the asymmetric case where the probability for single-electron transmission at resonance 4Γ U /Γ 1 is small, and expressions simplify. Omitting spin, the current operator has the form with the Breit-Wigner resonant function g(ε) = (ε − ε d + i Γ/2) −1 . The operatorĉ U/L (ε) annihilates an electron in a scattering state of energy ε incoming from the upper/lower lead. The normalisation is fixed by the average with the Fermi function f (ε) = (1 + e ε/k B T ) −1 . Due to the small capacitance at the upper dot-lead contact, the dc and ac bias voltages are applied essentially across this tunnel contact, the voltage potentials on both the quantum dot and the lower lead are fixed to the ground. Apart from the Breit-Wigner function, the rest of the calculation is similar to the case of a tunnel junction. The two-current correlators have the form of Eq. (2) and Eq. (3), where the equilibrium noise terms are given Eq. (11) in the main text.
We briefly discuss the case of a bare cavity damping κ 0 comparable to the electronic damping κ, but still much smaller than the resonator frequency ω 0 . The complete Heisenberg-Langevin equation (13) is solved by considering both the input field and electronic current fluctuations. One obtains for the two cavity field quadratures where denote their variances in the absence of intrinsic damping κ 0 , also given by Eq. (8) and discussed in length in the main text. Additionally, one finds for the output field squeezing, characterized by the power spectrum, Whereas a vanishing ∆X 2 1 (0) clearly optimizes squeezing in both the cavity and output fields, there is no such choice for κ 0 . Increasing κ 0 from zero improves squeezing in the output field but degrades cavity squeezing. Perfect squeezing in the output field is reached for κ = κ 0 , with vanishing ∆X 2 1 (0), in which case the cavity field is only half-squeezed.

CAVITY SQUEEZING BY A QUANTUM CONDUCTOR (SUPPLEMENTARY MATERIAL)
Cited equations not preceded by S− refer to the main text. Unless indicated otherwise, we set = 1 throughout the Supplementary material.

S-I: Current-current correlation function
In this section, we derive the current-current correlation function of an ac driven tunnel junction in the absence of cavity under the general form of Eq. (2) in the main text. We use the definition (S-1) The current operator through the junction is given bŷ corresponding to the operator where we incorporated the phase terms in the Fourier coefficients c n e inϕ → c n and defined ω kq = ε k − ε q . f (ε) = (1 + e βε ) −1 andf (ε) = 1 − f (ε) are the electron and hole Fermi-Dirac distribution functions. Substituting Eq. (S-5) in Eq. (S-1) and integrating over t 1 and t 2 , we obtain In deriving this expression, we used the standard assumptions of weak momentum dependence of γ kq and of the density of state in the leads, taken to be constant and denoted ν 0 . We introduced the tunnel resistance R T = 1/(2π|γ| 2 ν 2 0 ). Before pursuing our derivation of current-current correlation function, we define the unsymmetrized equilibrium Johnson-Nyquist noise for a tunnel junction, in absence of both dc and ac bias, asS (S-7) Using the expression Eq. (S-2) of the current operator, we readily find Using the above definition, we are able to rewrite the current-current correlator [Eq. (S-6)] as where the different noise terms are defined as Computing C ± with Eq. (S-12), we find an imaginary principal part term, describing frequency pull of the cavity, which we absorb into the cavity resonant frequency. We apply the rotating-wave approximation to the real part discarding all fast oscillating terms (beating with frequencies multiple of ω 0 ) and obtain We use these results in Eq. (S-14) and go back to the original time frame for the cavity field a(t) → a(t)e iω0t to find the Heisenberg-Langevin equationȧ

S-III: Finite-frequency Noise of a Quantum dot under ac excitation
We consider a single-level quantum dot attached to two leads, denoted upper and lower leads, with tunnel couplings t U and t L . Ignoring weak Coulomb interaction on the dot and assuming spinless electrons, the Hamiltonian takes the form with the constant density of states ν 0 for the two leads. ε d is the dot single-level energy.ĉ † α annihilates an electron with energy ε in lead α = U/L and satisfies the anticommutation relation [3] {ĉ † α (ε)ĉ α (ε )} = δ α,α δ(ε − ε ).

(S-2)
The transport of electrons in this system is described by using scattering states. The scattering matrix is given by [4] S U U (ε) S U L (ε) with the function g(ε) = (ε − ε d + iΓ/2) −1 . We have introduced the dot escape rates towards the two leads Γ α = πν 0 |t α | 2 , and the total escape rate Γ = Γ L + Γ U . We write the current operator in the presence of an ac driving of the two leads following Ref. [5], where the scattering problem in the dot is disentangled from the potential modulation and the photo-excitation in the leads [6]. Operatorsĉ α (ε) are introduced, which describe scattering modes just after entering or leaving the dot. They are assumed to be adiabatically connected to the lead scattering modes througĥ Eq. (S-4) describes the fact that an electron with energy ε arriving at the dot-lead contact is composed by a coherent sum of lead electrons with energies ε − 2mω 0 which have absorbed or emitted m photons from the ac field with probability amplitudes c α m . Writting the current in terms of operatorsĉ α (ε), as if no ac modulation were present, and substituting with Eq. (S-4), one obtains the general expression for the current operator in lead α I α (t) = e h dε 1 dε 2 n1,n2 e i(ε1−ε2)t γ,γ c γ n1 c γ * n2 ×ĉ † γ (ε 1 − 2n 1 ω 0 )A γγ (α, ε 1 , ε 2 )ĉ γ (ε 2 − 2n 2 ω 0 ),

Asymmetric quantum dot
Expressions simplify in the limit of an asymmetric dot, Γ U Γ L , more strongly coupled to the lower lead than to the upper lead. Computing the following scattering probabilities where we defined the scattering phase δ(ε) = arctan[Γ/(2(ε d − ε))], we see that they both scale as (Γ U /Γ) 2 . In contrast to that, for Γ U Γ L , we find A U L (U, ε 1 , ε 2 ) ig(ε 2 ) Γ U Γ, (S-10a) A LU (U, ε 1 , ε 2 ) −ig * (ε 1 ) Γ U Γ, (S-10b) and their square modulus scale as Γ U /Γ. We therefore neglect the LL and U U terms in the current. In addition, we assume that the voltage drop occurs at the upper dot-lead contact so that c L m = δ m,0 and we note c U m ≡ c m . The dc voltage is V ≡ V U at the upper lead and V L = 0 in the lower lead. The potential of the dot is also to the ground. With these different notations and approximations, the currentÎ U (t) is given by Eq. (D1) in the main text (see Methods) where the dc voltage is included in the operator time-dependence.

Fully symmetric quantum dot
We briefly discuss the case of a quantum dot symmetrically coupled to the leads Γ U = Γ L and to the cavity. The first order cavity-dot then takes the form iλ(Î U −Î L )(â † −â) involving the upper and lower currents with equal weights. The general approach devised for the tunnel junction also applies in this case, with the currentÎ being replaced byÎ U −Î L in the Heisenberg-Langevin equation (Eq. (6) in the main text). Settinĝ I =Î U −Î L (S-11)