Squeezing with a flux-driven Josephson parametric amplifier

Josephson parametric amplifiers (JPA) are promising devices for applications in circuit quantum electrodynamics (QED) and for studies on propagating quantum microwaves because of their good noise performance. In this work, we present a systematic characterization of a flux-driven JPA at millikelvin temperatures. In particular, we study in detail its squeezing properties by two different detection techniques. With the homodyne setup, we observe squeezing of vacuum fluctuations by superposing signal and idler bands. For a quantitative analysis we apply dual-path cross-correlation techniques to reconstruct the Wigner functions of various squeezed vacuum and thermal states. At 10 dB signal gain, we find 4.9+-0.2 dB squeezing below vacuum. In addition, we discuss the physics behind squeezed coherent microwave fields. Finally, we analyze the JPA noise temperature in the degenerate mode and find a value smaller than the standard quantum limit for phase-insensitive amplifiers.


Introduction
The tremendous progress in the field of quantum electrodynamics (QED) with solidstate superconducting circuits [1][2][3][4] has recently triggered massive efforts aiming at the investigation of propagating quantum microwaves [5][6][7][8][9][10][11][12][13][14]. Towards this end, the analysis of propagating microwave light at frequencies of a few gigahertz and with power levels below a single photon on average has become an important task. However, due to the low signal energy of typically only a few attowatt per megahertz bandwidth, the measurement of such signals requires amplification. For a long time, phaseinsensitive high electron mobility transistor (HEMT) amplifiers are considered as a good choice. They feature a broad operation bandwidth, high gain, but still add 10-20 noise photons [15,16] to the signal. This noise can be significantly reduced by using Josephson parametric amplifiers (JPA) [17][18][19][20][21][22][23][24][25][26], which achieve amplification via a high-frequency modulation of a Josephson inductance. In the phase-insensitive or nondegenerate operation mode, the JPA noise temperatures approach the standard quantum limit dictated by the Heisenberg uncertainty relation [5,8,20,26,27]. Even better, JPAs can in principle amplify a single signal quadrature without adding any noise [15]. This property is intimately connected to the fact that a JPA can squeeze a single quadrature below the level of the vacuum fluctuations. Hence, in combination with lately developed state reconstruction methods for propagating quantum microwaves [6,7,9], a JPA is the ideal tool to study the rich physics of squeezed microwave states. Recently, the squeezed vacuum generated by the JPA investigated in this work was successfully used to generate distributable path entanglement [7]. In addition, squeezing is considered as a key ingredient for the creation of quantum memories for entangled continuous-variable states [28].
Here, we present a comprehensive study on the squeezing physics produced by a flux-driven JPA [24]. The particular advantage of this design is the large isolation between the pump and signal ports. For this reason, the elimination of the pump tone becomes obsolete [22] and, thus, experimental complexity is greatly reduced. After characterizing the basic properties of our JPA with a spectroscopy setup and a standard homodyne detection scheme, we reconstruct the Wigner functions of the JPA output states using dual-path cross-correlation techniques [6]. The latter allows us to investigate squeezed vacuum, squeezed thermal and squeezed coherent states. We find a maximum squeezing of 4.9 ± 0.2 dB and confirm that the dependence of photon number and displacement on the squeezed and coherent state angle follow the expectations from theory. In addition, a detailed noise analysis shows that our JPA is operated near the quantum limit in the phase-insensitive mode and that its noise temperature is below the standard quantum limit in the degenerate mode. All in all, our work illuminates in unprecedented detail the fundamental physics of single-mode squeezing in the microwave domain.
The paper is organized as follows: In section 2, we discuss the characterization of the flux-driven JPA using spectroscopy and homodyne detection. After introducing sample and measurement setups, we analyze signal and idler gain, the associated bandwidths, compression effects, and the noise properties in non-degenerate mode. Next, we investigate in section 3 the squeezing of vacuum and thermal states with both homodyne detection and the dual-path state reconstruction method. In section 4, we extend our discussion to squeezed coherent states by operating the JPA in the degenerate mode with coherent input signals. After discussing the JPA noise properties in the degenerate mode based on dual-path measurements in section 5, we summarize our main results and give a short outlook in section 6.

JPA sample
In this subsection, we describe the operation principle of a flux-driven JPA in general and our sample in particular. A parametric amplifier is an oscillator whose resonant frequency is modulated periodically in time. In the case of the JPA, the oscillating system is a quarter-wavelength transmission line resonator whose resonant frequency is determined by its capacitance and inductance (see figure 1(a)). The latter can be varied by a dc superconducting quantum interference device (SQUID), which consists of a superconducting loop interrupted by two Josephson junctions and acts as a fluxdependent non-linear inductor. Thus, by modifying the magnetic flux Φ dc threading the SQUID loop, the resonant frequency can be adjusted (see figure 1(b)). By fitting a physical model to the experimental data (black squares), we can estimate a Josephson coupling energy E J = h × 1305 GHz for each junction, where h is the Planck constant.  The slightly different values compared to those in the supplementary material of [7] result from the fact that, here, we employ a more sophisticated distributed-element model [29] instead of a simple lumped-element approach.
Periodically varying the resonant frequency with an ac flux (pump tone) at 2f 0 , where f 0 is the operation point frequency, results in parametric amplification: A signal at f 0 − f impinging at the signal port is amplified by the signal gain G and reflected back out of the signal port. At the same time, an idler mode at f 0 + f is created, whose amplitude is determined by the intermodulation gain M . This operation principle is depicted in figure 1(c). If the incoming signal consists of vacuum fluctuations, this process is the analogue of parametric downconversion in optics, where a pump photon is split into a signal and an idler photon. Therefore, strong quantum correlations between the signal and idler mode are established which finally result in squeezing.
An optical micrograph of the sample is shown in figure 2. The resonator and antenna are patterned from a sputtered 50 nm thick Nb film. At the contacts, 95 nm of gold is deposited on top of a 5 nm titanium bonding layer. As substrate we use thermally oxidized (300 nm) silicon with a thickness of 300 µm. In the last step, the aluminum dc SQUID is fabricated using shadow evaporation [30]. Each Al electrode has a thickness of 50 nm. The sample chip is placed between two alumina printed circuit boards inside a gold-plated copper box.

Spectroscopy and homodyne setup
The spectroscopy setup used to characterize the JPA and the homodyne setup to detect squeezing are sketched in figure 3. In JPA characterization and homodyne measurements, the same cryogenic setup (see figure 3(a)) is used. The sample is anchored to the sample stage of a dilution refrigerator and its temperature ranges between 90 mK and 130 mK. The signal, generated either by a microwave source or a vector network analyzer (VNA), passes a series of warm (-40 to -60 dB) and cold attenuators (−69 dB) as shown in figure 3(a). The signal power levels stated in this work are referred to the output of the 30 dB-attenuator in figure 3(a), while the pump power level is estimated at the input of the JPA sample box. The transmission from the 30 dB-attenuator output to the spectrum analyzer has been determined by sweeping the temperature of the 30 dB-attenuator and measuring the power of the emitted black body microwave radiation [31]. We subtract this value from the total transmission between the microwave source output and the spectrum analyzer to calibrate for the loss of the input line. For signal and idler gain measurements, a coherent signal is fed through the attenuated input line via the measurement circulator to the JPA. This circulator separates the outgoing from the incoming signal and protects the JPA from the noise generated by the amplification chain. The amplified output signal can be detected by a VNA or a spectrum analyzer (see figure 3(b)). The VNA allows to investigate the complex reflection coefficient of the JPA. The spectrum analyzer is used to investigate the idler gain and the degenerate operation of the JPA (see subsection 2.8). Figure 3(c) shows the homodyne receiver used in the squeezing measurements presented in section 3. It consists of a local oscillator microwave source and a mixer which downconverts the signal

Operation Point
The first step in characterizing the JPA is to determine a suitable operating point. To this end, we measure the signal gain with a VNA while synchronously sweeping the pump tone, fulfilling the relation f pump = 2f signal + 10 kHz. Since the VNA measurement bandwidth of 30 Hz is much smaller than 10 kHz, only the signal mode is detected preventing interference effects present in the degenerate mode. We emphasize that this measurement is different from the signal bandwidth measurements discussed in subsection 2.5, where the pump tone is at a fixed frequency and only the signal frequency is swept. In figure 4, we have plotted the measured signal gain when synchronously sweeping the pump and signal frequency for different flux values. For lower frequencies, the signal gain is increasing because the dependence of the resonant frequency on the flux becomes steeper (see also figure 1(b)). We choose our operating point f 0 between 5.634 GHz and 5.639 GHz, depending on the measurement. As it can be seen from figure 4, our operation point is located in the center of a region where the signal gain is appreciable and its frequency dependence is well behaved. At this operation frequency, the external quality factor is measured as Q ext = 312, and the isolation between antenna and resonator is at least 28 dB.

Non-degenerate gain
When the signal frequency is detuned from half the pump frequency, signal and idler modes are at different frequencies and can be observed individually. This mode of operation is therefore called "non-degenerate mode". Figure 5 shows the pump power dependence of the non-degenerate signal and idler gain for a detuning of 10 kHz. For low pump power, no significant signal gain is observed and the idler gain is small. For large pump power, the two curves converge and both idler and signal gain reach 19.5 dB.

Bandwidth
Besides the gain properties described above, bandwidth is an important feature of an amplifier. Therefore, we determine the instantaneous bandwidth of signal and idler modes at a fixed operation point (constant flux and pump frequency) by measuring the signal and idler gain for various detuning between half the pump frequency f pump /2 and the signal frequency f . As shown in figure 6, we observe a signal and idler bandwidth of 1.72 MHz for a pump power of −39 dBm. We define the gain-bandwidth-product (GBP), which is defined as the product of the voltage gain in linear units and the bandwidth of our amplifier. In the large gain limit, the GBP is nearly constant and close to the theoretical limit of f 0 /Q ext = 18 MHz [32] (see figure 6(c)). Going to low values of the signal gain, the idler gain vanishes and the signal gain approaches one because the signal gain is normalized to the pump-off condition. Therefore, we expect the idler GBP to decrease and signal GBP to increase. However, in the low gain limit the signal bandwidth measurement is very sensitive to the calibration data and fluctuations. Consequently, the signal GBP does not diverge in contrast to expectations from theory.

1 dB-compression point
Another important figure of merit for amplifiers is their 1 dB-compression point. It denotes the power where the signal gain is 1 dB below the value expected for a perfectly linear device [26]. In other words, at some point the signal gain starts to decrease as a function of input power due to the nonlinearity of the amplifier. In figure 7(a) bandwidth measurements of the signal gain are displayed. In contrast to figure 6, here the pump power is fixed at -39 dBm and the dependence on the signal power is studied. For small  signal powers (P signal ≤ −136 dBm) the curves overlap meaning that the gain is constant and that the amplifier is in the linear regime. For larger signal powers a reduction of the gain is observed. Figure 7(b) shows the maxima of the signal gain of figure 7(a) for several signal powers. The 1 dB-compression point occurs at −133 dBm. An analogous analysis for the idler gain (data not shown) reveals that the compression occurs for both the idler and signal gain at the same input power. For a circuit QED experiment with a cavity decay rate of 1 MHz, this power level is equivalent to 10 photons on average.

JPA noise properties in non-degenerate mode
A low noise temperature of an amplifier is very important when using it for measuring signals at the quantum level. Here, we use the spectroscopy setup for a rough characterization of the noise properties of our JPA when it is operated in the phaseinsensitive mode. For details on the noise temperature in the degenerate mode, we refer to section 5. In order to determine the JPA noise properties, we measure the amplified noise power emitted by the 30 dB-attenuator, whose temperature is controlled in the range T att = 50 − 800 mK. The total noise power P at the spectrum analyzer is given by where G denotes the total gain, B the detection bandwidth, h = 6.626 × 10 −34 J · s the Planck constant, k B = 1.38 × 10 −23 J/K the Boltzmann constant and T total the total noise temperature of the complete detection chain, which includes the measurement circulator, the JPA and the amplification chain. T total is related to the number of photons n total added by the complete detection chain by k B T total = n total hf 0 . The first term in (1) describes thermal fluctuations and vacuum fluctuations according to [33]. Possible deviations between the electronic temperature of the attenuator and the measured temperature are taken into account by δT . We set the signal gain to G signal (T →0) = 19 dB at the operation point f 0 . However, sample heating due to the pumping process and compression effects at high noise source temperatures may cause the signal gain to deviate from this value. Therefore, we measure the dependence of the gain-corrected power on the noise source temperature T att (see figure 8). To this end, we implement the following protocol for each temperature point: after measuring the signal gain with the VNA, we turn off the VNA and measure the total noise power using a spectrum analyzer. We obtain the gain corrected power by dividing the total noise power at each temperature point by the effective JPA gain for white uncorrelated noise, G eff (T ) = 2G signal (T ) − 1. Here, we consider the idler gain using G idler = G signal − 1 [18]. Taking into account the cable and connector losses between the noise source and the measurement circulator, we estimate T total = 167 mK, which corresponds to n total = 0.62 for the noise photons added by the whole detection chain referred to the input of the measurement circulator. This value is close to the standard quantum limit for phaseinsensitive amplifiers of 0.5 photons (135 mK). In addition, we find δT = − 19 mK to be reasonably small.

Degenerate gain
In order to demonstrate the usability of our JPA as a phase-dependent amplifier, we investigate the degenerate gain. In this mode of operation, the pump frequency is twice the signal frequency. Thus, the idler mode is created at the frequency of the amplified signal mode. This results in an interference of the two modes which is constructive or destructive depending on the phase between the idler and the signal modes, which can be controlled by shifting the phase difference between the probe signal and the pump tone. The degenerate gain is measured with a spectrum analyzer, where a measurement with zero pump power is used as a reference. In figure 9, the degenerate gain is plotted as a function of the phase between the probe signal and the pump for different pump power levels. We observe a maximum degenerate gain of 25.  Phase (degree) Figure 9. Degenerate gain as a function of the phase difference between probe and pump signal for different values of the pump power. For clarity, the curves are shifted in phase direction, so that the minima coincide. The temperature of the JPA is stabilized at 95 mK.

Squeezing of vacuum and thermal fluctuations
As mentioned in the discussion of the operation principle (subsection 2.1), the JPA creates quantum correlations between the signal and idler modes. In the degenerate operation mode, these correlations result in deamplification or amplification depending on the quadrature direction. In the case of deamplification, the quadrature fluctuations can be squeezed below those of the vacuum. In this section, we investigate squeezed vacuum fluctuations generated by the flux-driven JPA with two different detection techniques: homodyne detection and dual-path state reconstruction [6,7].
To define the squeezing level S in decibel, we compare the variance of the squeezed quadrature (∆X sq ) 2 with the quadrature variance of vacuum fluctuations, obtaining S = max{0, −10 lg (∆X sq ) 2 /0.25 }.
We note that (∆X sq ) 2 < 0.25 indicates squeezing and S is positive. Larger S corresponds to a higher squeezing level. (∆X sq ) 2 ≥ 0.25 indicates no squeezing and S equals zero. Hence, in this work we use the nomenclature that the term "squeezing" is equivalent to "squeezing below the vacuum level".

Squeezing detected with the homodyne setup
First, we detect the output signal of the JPA with the homodyne detector shown in figure 3(c). By sweeping the temperature of the 30 dB-attenuator from 50 mK to 800 mK with inactive JPA, we calibrate the gain and the noise temperature of the detection chain using (1). Then, we cool the attenuator to 33 mK and investigate the squeezed vacuum state. In figure 10(a), the ratio dS ≡ S on /S off is plotted as a function of the relative phase between the pump and the local oscillator and the detected intermediate frequency f IF . The latter is obtained by downconverting the signal using a fixed local oscillator frequency f LO = f 0 = f pump /2. For these settings, signals at f LO + f IF and f LO − f IF , representing the signal and idler modes, are downconverted to the same intermediate frequency f IF . Therefore, the homodyne detector is sensitive to the correlations between the two modes created by the JPA. Here, S on and S off are the power spectral densities recorded with the JPA pump power on and off, respectively. In figure 10(b), the average dS avg of dS calculated in the frequency range 118 kHz < f IF < 487 kHz is plotted. Whenever dS avg < 0 dB, the noise detected at the spectrum analyzer referred to the input of the amplification chain is smaller than the vacuum noise emitted by the attenuator. This demonstrates the effect of vacuum squeezing. From the minimal value dS min indicated in the inset of figure 10(b), we calculate the squeezing level [27] S = −10 lg 1 − T n 1 − 10 dS min /10 0.5hf 0 /k B .
Here, T n is the noise temperature of the detection chain with the JPA off. Taking into account the cable losses, reference spectrum fluctuation and thermal population at the input of the JPA, we retrieve a lower bound of 2.8 dB of squeezing at the input of the amplification chain.

Dual-path setup
The estimation of the squeezing level with the homodyne detection setup only yields a lower bound for the amount of squeezing. Rather than improving this setup, we choose to fully reconstruct the squeezed state emitted by the JPA with the dual-path setup, which is based on cross-correlation techniques and realtime data processing [6,7]. First, we introduce the dual-path setup in figure 11. We split the signal under study described by the bosonic annihilation and creation operatorsâ andâ † using a microwave beam splitter and feed them into two amplification and detection paths. During the splitting process vacuum fluctuations are added to the split signals. While the signal is emitted from our JPA sample, the vacuum fluctuations are realized by terminating the other input port of the beam splitter with a broadband 50 Ω load. At the two output ports, we first amplify the signals using cold HEMT and room temperature amplifiers and then downconvert them to an intermediate frequency (IF) of 11 MHz using IQ-mixers. The resulting orthogonal quadrature signals I 1,2 and Q 1,2 are then digitized by four analogdigital-converters (ADCs), and processed in real time by a field programmable gate array (FPGA) logic. Based on the beam splitter relations and the fact that the noise contributions from two detection paths are independent, we get access to all quadrature moments of signal and the two detection paths up to fourth order. In the setup, we still have the same temperature-controllable 30 dB-attenuator at the end of the JPA input line to calibrate the photon number conversion factors which relate the autocorrelations measured in our detector (in units of V 2 ) to photon numbers at the input of the beam splitter.
After post-processing the ensemble averages of the noisy quadrature moments, the reconstructed signal and noise moments have the form of â † lâ m and V r and r + s ≤ 4. Following [7], with the third and fourth moments we calculate the third and fourth order cumulants to verify the Gaussianity of the state. Furthermore, we use the first two moments to reconstruct the signal Wigner functions at the input of the beam splitter. At the same time, the noise moments of the two detection paths are obtained. We find that the detection chains add 24.22 ± 0.02 and 27.32 ± 0.03 photons referred to the beam splitter input which corresponds noise temperatures of 6.55 ± 0.01 K and 7.39 ± 0.01 K for the detection path 1 and path 2, respectively. Again, the error bars we provide are of purely statistical nature. In all dual-path experiments, the temperature of the JPA sample is stabilized at 50 mK.

Squeezing detected with the dual-path setup
Selected Wigner function reconstructions of the state at the input of the beam splitter are plotted in figure 12. When the 30 dB-attenuator is stabilized at 50 mK, vacuum fluctuations are present at the input of the JPA. As shown in figure 12(a), the increase in pump power, which corresponds to an increase of signal gain, causes an increase of the squeezing level and an increase of the photon number (see also figure 13(a)). We achieve a maximum squeezing level of 4.9 ± 0.2 dB below vacuum at 10 dB signal gain. However, if we further increase the signal gain, the squeezing level decreases again. This behavior is expected [34] because the squeezing becomes suppressed when the JPA enters the bifurcation regime. In this regime, also the higher order cumulants do not vanish anymore. Indeed, we observe this effect in our data for signal gains larger than 10 dB. In addition, from the JPA input output relation [35] we identify the signal voltage gain in the non-degenerate mode in linear units as G signal,V = cosh (r). Applying this relation, we obtain from the photon number of a squeezed state, n = sinh 2 (r), the expression n = G signal,P − 1, where G signal,P = G 2 signal,V is the signal power gain in linear units. Therefore, we expect that the photon number increases linearly with G signal,p in the non-degenerate mode with a slope of one. Figure 13(b) confirms this behavior for small signal gains below the bifurcation regime. Next, we fix the signal gain at 1 dB. When the temperature of the 30 dB-attenuator is increased (see figure 12(b) and figure 13(c)), more and more thermal photons are incident at the input port of the JPA. Thus the squeezing level decreases and at some point the output state of the JPA is not squeezed below vacuum any more. In addition to the squeezed vacuum and squeezed thermal states discussed in the previous section, we here analyze a more general squeezed state: the squeezed coherent state. There are two ways to obtain a squeezed coherent state. First, one can squeeze the vacuum stateŜ (ξ) |0 and apply a displacement operationD (α)Ŝ (ξ) |0 . Here, ξ denotes the complex squeeze parameter ξ = r exp(iϕ) with squeeze factor r and squeezing angle ϕ. Second, one can apply a displacement operator on the vacuum stateD (α) |0 , and subsequently squeeze this displaced vacuumŜ (ξ)D (α) |0 . We define the coherent phase θ as the angle between displacement direction and p-axis, and the anti-squeezed angle γ = − ϕ/2 as the angle between the anti-squeezed quadrature and the p-axis. We illustrate the difference of the two methods in figure 14. For the former method, the final displacement of the squeezed coherent state only depends on the displacement operation and is independent of the squeeze factor r (figures 14(a)-(c)). However, the displacement of the squeezed coherent state obtained from the second method depends on both the displacement and squeezing operations. When the anti-squeezed quadrature is parallel to the displacement direction of a coherent state, the final displacement of the squeezed coherent state is maximum ( figure 14(d)). Contrary, the final displacement reaches its minimal value when the anti-squeezed quadrature is perpendicular to the displacement direction ( figure 14(e)). In this section, we present an analysis on squeezed coherent states generated with the second method and detected with the dual-path setup. To this end, we pump the JPA in the presence of a weak coherent signal. One measurement trace always records four regions: vacuum (V) region with both JPA pump and coherent signal off, coherent signal (C) region with JPA pump off and coherent signal on, squeezed coherent (SC) region with both JPA pump and coherent signal on, and squeezed vacuum (SV) region with JPA pump on and coherent signal off. Figure 15 shows typical averaged time traces of selected first and second moments from a single measurement, which corresponds to the situation in figure 14(f). Four Wigner function reconstructions are placed on top of the time trace. When both JPA pump and coherent signal are off, we can clearly identify the vacuum Wigner function with its maximum centered at the origin. Once we turn on the coherent signal, which is a displacement operation, the vacuum state is equally displaced along p-and qquadratures and becomes a coherent state with phase Θ = 45 • referred to the beam  Figure 17. Experimentally obtained photon numbers for coherent states, squeezed vacuum states and squeezed coherent states as a function of the coherent state phase Θ. The statistical uncertainty is smaller than the symbol size. The dashed curves are fits using (3). We note that the thermal contribution δN ≈ 0.005 is vanishingly small. Furthermore, although this fit is performed independently from the one in figure 16, the obtained fit parameters are consistent. The data are from the same measurements as shown in figure 16. splitter input. Next, keeping the coherent signal on, we turn on the JPA pump which gives a squeezing operation with the p-quadrature as the squeezed quadrature and the q-quadrature as the anti-squeezed quadrature (γ = 90 • ). This results in a suppression of the displacement of the state along the p-quadrature and a corresponding amplification of the displacement along the q-quadrature. Also the p-quadrature variance is squeezed and the q-quadrature variance is amplified, turning the circular profile of the vacuum and coherent state Wigner function into an elliptic profile. Thus, the state becomes a squeezed coherent state with squeezing level S = 4.3 dB. Finally, keeping the pump on, we turn off the coherent signal to generate a squeezed vacuum state with squeezing level S = 4.7 dB. In this context, we would like to point to the following experimental aspect. Compared with a coherent state, whose phase Θ is referred to the beam splitter input, in the case of a squeezed coherent state, we use θ as the phase of the coherent state before applying the squeezing operator, and θ is referred to the JPA input. During the propagation of the coherent state from the input of the JPA to the input of the beam splitter, the phase of the coherent state evolves. We account for this effect by a constant phase difference ∆Θ ≡ Θ − θ. We note that the angle γ always refers to the squeezed vacuum state recorded in each measurement trace. From theory [35], we expect the displacement of a squeezed coherent state after the squeeze operation to depend on the angles as â = α cosh r − α * e iϕ sinh r. ( Here, α = |α| exp [iπ (90 • − θ) /180 • ] is the complex amplitude of the coherent state before the squeeze operation and ϕ = − 2γ the angle of the complex squeeze parameter. By fixing the anti-squeezing angle of squeezed vacuum states γ at 0 • , 45 • , 90 • , and 135 • and rotating the phase Θ of the coherent signal, we map out the dependence of the displacement of a squeezed coherent state on γ and θ. In figure 16, we display the displacement, which is the center of the individual states in phase space given by their first moment â for various values of Θ and γ. The squeezed vacuum states are centered at the origin, and the coherent states are located on a circle around the origin. If we turn on the JPA pump and rotate the phase of a coherent signal, the squeezed coherent state moves mainly along the γ direction. The displacement of the squeezed coherent states reaches its maximum when γ = θ + 2n × 90 • , where n ∈ Z. Geometrically, this means that the anti-squeezed direction is collinear to the displacement vector pointing from the origin to the center of the state. For our choice of γ, we obtain a characteristic star-shaped pattern.
Similar to the center of the Wigner function, which represents the displacement, the photon number of a squeezed coherent state varies when we rotate the phase of the coherent signal Θ while keeping the anti-squeezed angle γ constant. Following [35], we obtain: â †â = |α| 2 + δN cosh 2 r + sinh 2 r − (α * ) 2 e iϕ cosh r sinh r −α 2 e −iϕ cosh r sinh r + sinh 2 r, where δN describes the thermal photons present in the vacuum state at the JPA input and |α| 2 is the number of photons in the coherent state. As we see from figure 17, the photon number oscillates and reaches a maximum when γ = θ + 2n × 90 • , n ∈ Z. Thus, the photon number is maximal when the displacement of the SC state is maximal. We From theory, we expect that the only difference between a SC state and the corresponding SV state is a displacement in phase space without any rotation or deformation. Therefore, we analyze the statistics of the variances (∆X anti ) 2 and (∆X sq ) 2 of the anti-squeezed and squeezed quadratures and that of the angle γ for the SC and SV states. As displayed in table 1, our data show that no significant rotations or deformations are present.

JPA noise properties in degenerate mode
A JPA operated in the degenerate mode can not only generate vacuum squeezing, it can also be used as a low-noise phase-sensitive amplifier which, in principle, does not need to add any noise to the amplified quadrature [15]. With the dual-path setup, we therefore study the noise properties of our JPA in the degenerate mode. More precisely, we perform a temperature sweep of the 30 dB-attenuator shown in figure 11. The variance of the fluctuations at the frequency f 0 generated with this procedure is where (∆X therm ) 2 has the unit of photon number. At each temperature, the JPA pump is operated in the pulsed mode. Consequently, a single time trace in our measurement always contains a region corresponding to a non-squeezed thermal state and a region corresponding to a squeezed thermal state. For any quadrature, the variance (∆X out ) 2 at the output of the JPA is related to the variance (∆X in ) 2 at the input via the relation [15] (∆X out ) 2 = G X (∆X in ) 2 + (∆X noise ) 2 .
Here, G X is the gain for this quadrature and (∆X noise ) 2 is the noise added by the amplifier referred to the output. In principle, we could determine the variance of the thermal state at the input of the JPA using the dual-path reconstructed signal moments at the input of the beam splitter taking into account the cable, circulator and JPA losses. However, the dual-path reconstruction detects a thermal population of 0.1 photons in the vacuum [7] which would result in a significant underestimation of the JPA noise. For this reason, we calculate (∆X in ) 2 based on (4), and model the cable loss between the 30 dB-attenuator output and the measurement circulator input with beam splitters, and account for the temperature gradients. In this way, we model an equivalent amplifier consisting of the measurement circulator, JPA, and cables to the beam splitter input.
The noise contributions of all these components are represented by the noise temperature of the equivalent amplifier. Therefore, the latter is a pessimistic estimate for the noise properties of the JPA itself. In figure 18, we show the anti-squeezed and squeezed quadrature variance as a function of the noise source temperature. Obviously, there is good agreement between theory and experiment. The corresponding numerical results are displayed in table 2. Most importantly, we observe a variance of 0.14 ± 0.01 photons (referred to the input) for the noise added by our JPA to the anti-squeezed quadrature. This value is clearly below the standard quantum limit of 0.25 photons for a single quadrature of a phaseinsensitive amplifier. The relevant noise number for the use of the JPA as a squeezer is the noise it adds to the squeezed quadrature at the JPA output. We do not refer the added noise from the squeezed quadrature to the input of the equivalent JPA, because systematic uncertainties in the setup [7] dominate the squeezed quadrature noise variance at the JPA output and would be amplified by 1/G sq when referring to the JPA input. However, these uncertainties are negligible for the anti-squeezed quadrature. Instead, following [15] we calculate a lower bound for the squeezed quadrature noise variance at the JPA input from the experimentally more robust quadrature gains G sq and G anti of the squeezed and anti-squeezed quadrature and antisqueezed quadrature noise variance. Using the values from table 2, we obtain the relation (∆X sq,in ) 2 ≥ 1 16 1 − (G sq G anti ) −1/2 2 / (∆X anti,in ) 2 = 0.02, where subscript "in" indicates the JPA input.

Conclusions and Outlook
In summary, we present a detailed analysis of the physics of squeezed microwave light generated with a flux-driven JPA. We first determine the operation point of the JPA and characterize its basic amplification properties, including non-degenerate gain, bandwidth, 1 dB-compression point and degenerate gain. We then use the JPA to squeeze vacuum fluctuations and find 4.9 ± 0.2 dB of squeezing at 10 dB signal gain. Furthermore, we investigate displacement and photon number of squeezed coherent microwave fields and find excellent agreement with theoretical expectations. In the degenerate mode, we verify that our JPA, as a phase-sensitive device, adds less noise to the amplified quadrature than an ideal phase-insensitive amplifier. This property is of utmost importance for high efficiency detection, state tomography and quantum communication applications in the microwave domain. Furthermore, with the setup used in this work, i.e., squeezed states incident at a linear beam splitter, the generation of path entanglement between continuous-variable propagating quantum microwaves was demonstrated in [7]. Considering recent work on the engineering of tunable beam splitter Hamiltonians [36][37][38], our setup could be extended to the interesting case of dynamical switching between two separable single-mode squeezed states and a pathentangled two-mode squeezed state.