Observation of quantum many-body effects due to zero point fluctuations in superconducting circuits.

Electromagnetic fields possess zero point fluctuations which lead to observable effects such as the Lamb shift and the Casimir effect. In the traditional quantum optics domain, these corrections remain perturbative due to the smallness of the fine structure constant. To provide a direct observation of non-perturbative effects driven by zero point fluctuations in an open quantum system we wire a highly non-linear Josephson junction to a high impedance transmission line, allowing large phase fluctuations across the junction. Consequently, the resonance of the former acquires a relative frequency shift that is orders of magnitude larger than for natural atoms. Detailed modeling confirms that this renormalization is non-linear and quantum. Remarkably, the junction transfers its non-linearity to about thirty environmental modes, a striking back-action effect that transcends the standard Caldeira-Leggett paradigm. This work opens many exciting prospects for longstanding quests such as the tailoring of many-body Hamiltonians in the strongly non-linear regime, the observation of Bloch oscillations, or the development of high-impedance qubits.

Introduction

The realization of many-body effects in quantum matter, often associated with remarkable physical properties, hinges on strong interactions between constituents. Mechanisms to achieve strong interactions include the Coulomb interaction in narrow band electronic materials, and Feshbach resonances that can produce arbitrarily large scattering lengths in ultracold atomic gases. In contrast, while providing great design versatility, purely photonic platforms[1-3], are not easily amenable to realizing strong correlations, since they usually come with weak nonlinearity. To circumvent this, superconducting circuits, which operate in the microwave range and display high tunability, have been proposed[4] for the exploration of correlated states of light. Here, correlations originate from nonlinear elements, such as Josephson junctions, and the enhancement of nonlinearities is accompanied by large zero-point fluctuations (ZPF). This can be understood in the electronics language of impedance as follows. The dynamics of a Josephson junction is described by two conjugate variables: the number of transferred Cooper pairs and the superconducting phase difference . Despite being an anharmonic oscillator, a Josephson junction with Josephson energy and charging energy , can be associated to an impedance which sets the amplitude of the fluctuations of and . When is sufficiently smaller than unity, and , with the superconducting quantum of resistance. Consequently, at low , phase fluctuations are weak and the anharmonic Josephson cosine potential can be reduced to a quadratic potential plus a quartic perturbation, as is the case for the transmon qubit[5]. On the other hand, if is large, the full cosine potential is explored due to strong phase fluctuations. Anharmonicity then becomes important, as observed with the Cooper-pair box[6] or the fluxonium qubit[7], and as a result, the oscillation frequency can strongly deviate from the harmonic value . Thus, exploring many-body physics in circuit quantum electrodynamics must rely on a careful tailoring of ZPF.

The approach[8] that we follow to explore many-body effects originating from a single nonlinear superconducting element is to couple it to many harmonic modes. In the presence of such an environment, the degree of anharmonicity of a nonlinear junction will also depend on the external impedance, and three regimes can be identified. When does not match the environmental impedance at frequencies close to , the junction is accurately described as an almost isolated system, so that the effect of the environment only amounts to small perturbative corrections, similar to the Lamb shift[9]. At the same time, an impedance-mismatched environment remains weakly perturbed by the nonlinear junction, and this absence of back-action allows it to be described as a set of harmonic oscillators, following the Caldeira–Leggett approach[10]. This simplified description is at the core of the current understanding of open quantum systems, and was already verified experimentally in the early studies of macroscopic quantum tunneling[11]. The important role of ZPF in the damping effect that such an environment has on a Josephson junction was already noticed experimentally[12] and explained theoretically[13,14] three decades ago. In these early works, the effect of ZPF was to renormalize junction properties such as the critical DC current by about , which nonetheless had a large effect on macroscopic quantum tunneling rates. When , the junctions and its environment fully hybridize, since they are impedance matched, but the anharmonicity of the junction remains weak and can be treated pertubatively[15-18]. The case is much more challenging, both experimentally and theoretically since the strongly anharmonic junction hybridizes with many modes of its environment. In DC measurements, such effects result in the celebrated Schmid–Bulgadaev transition predicted more than 30 years ago[19,20], a localization phenomenon whose relevance for microwave AC measurements requires further experimental and theoretical investigations[21]. The environment provides a strong action on the junction, which itself induces a sizeable back-action on many modes of the environment, the combined circuit forming a complex many-body system reminiscent of quantum impurity problems encountered in condensed matter[22]. More specifically, the frequency shift induced by the environment on the junction can be comparable with , a nonperturbative effect due to a modification of the vacuum[23]. At the same time, the nonlinearity of the junction is transferred into the environmental modes, affecting for instance their broadening, and producing a physical regime that was not addressed so far.

In this work, we report on the effects of ZPF in a device consisting of a fully characterized multimode environment and a highly nonlinear single Josephson junction, acting as a weak link between two linear transmission lines, with all subsystems reaching the high-impedance regime. As a result, the transmission of single photons through our device is strongly affected by the interplay of nonlinearities and ZPF. We observe a 30% renormalization of the junction frequency as compared with the value that would have been obtained without ZPF—analogous to a giant Lamb shift—and we provide clear evidence for modifications of the environmental vacuum, which inherits strong nonlinear effects. A detailed temperature analysis of our system proves the quantum origin of these fluctuations and eliminates an explanation in terms of classical hybridization effects. Finally, our experimental findings are in quantitative agreement with a microscopic theory based on the self-consistent harmonic approximation (SCHA), embedded within a fully fledged microscopic description of our circuit using microwave simulation tools.

Results

Background

The many-body regime of a single nonlinear junction coupled to a high-impedance environment has remained largely unexplored experimentally, since obtaining at gigahertz frequencies is very challenging. One option is to use on-chip resistors[24]. However, this may lead to unwanted Joule heating[25]. Therefore, we rather pursue a solution that relies on superconducting (lossless) high-inductance materials such as Josephson junction arrays[26-28], noting that disordered superconductors[29] are also promising. In Josephson junction arrays, can reach given the large inductance of these materials, while maintaining good quality factors in the device.

Early experiments have embedded ultrasmall Josephson junctions between highly resistive leads, demonstrating the incoherent tunneling of Cooper pairs[24] in the framework of the theory[30]. In this case however, no supercurrent flows through the junction and no quantum coherent effects were observed. Later, the phase/charge duality in the regime ,  >  was explored using SQUID arrays as the environment[31-33]. Experimental results were explained by fluctuations due to the finite temperature of the electromagnetic environment and the effect of ZPF could not be investigated. Moreover, these two series of experiments relied on DC measurements. This has the disadvantage that nonequilibrium effects need to be taken into account when results are interpreted, while the system is not directly probed at the finite frequencies—around —that are of greatest interest.

It has since become possible to obtain a frequency-resolved picture of the environment of quantum systems such as Josephson junctions, thanks to the advent of circuit QED[34]. Here, microwave techniques allow a more accurate examination of the effects of ZPF on Josephson junctions[35], and observations of perturbative spectral shifts (below 1%) attributed to ZPF were reported[9,36,37]. Several bottom-up experiments explored nonperturbative effects of light–matter interaction at ultrastrong coupling between a qubit and a single-mode resonator (for a review see Forn-Dìaz et al.[38] or Kockum et al.[39]). An effect similar to the Lamb shift—a reduction of the effective Josephson energy—was also reported recently for a DC-biased Josephson junction coupled to a single mode high-impedance resonator[40]. Moving towards many-body territory, a nonperturbative renormalization of the frequency of a flux qubit was demonstrated[41,42]. However, in this experiment, fluctuations were mainly thermal, and in addition, the environment cutoff frequency could not be clearly measured. The resulting unknown parameters prevented a quantitative modeling of the experiments. Indeed, as pointed by various authors[43-46], it is necessary to account for all the microscopic details of the circuit to get rid of unphysical divergences in multimode models. Furthermore, a thorough modeling of such circuits is mandatory to discriminate the trivial effects of normal mode splitting (spectral shifts observed when two classical harmonic oscillators hybridize) from the dynamical ones associated to true vacuum fluctuations. With the exception of  Gely et al.[47], this important issue has received surprisingly little attention in the circuit-QED context.

Presentation of the experiment

Our system builds on recent advances in the fabrication and control of large-scale Josephson arrays[18,48]. It consists of a small Josephson junction of characteristic impedance on the order of (), which is embedded in the middle of two SQUID chains, each consisting of 1500 unit cells (Fig. 1), forming high-characteristic impedance transmission lines. We measure the characteristics of this environment precisely: its high-frequency cutoff—or plasma frequency—  17 GHz and its wave impedance (see Table 1 and Supplementary Note 10). The SQUID parameters were carefully chosen to maintain a negligible phase slip rate (), ensuring that these chains can be described as a linear environment. They are capacitively coupled to the measurement setup to suppress DC noise which could affect the small junction (Fig. 1c). In order to vary the degree of nonlinearity and hence the strength of the ZPF, we measured three samples with different small junction sizes, connected to nominally identical chains (see Table 1).

Fig. 1

SQUID chains coupled to a small Josephson junction (weak link). The upper part represents the spatial phase distribution of the two first standing waves or resonant modes of the total system (Josephson junction + chains). And odd (even) mode—which couples (does not couple) to the junction—is represented in purple (orange). The lower part is a schematic of the system. The SQUID chains, depicted as blue transmission lines, are capacitively coupled to the input and output 50 Ω coaxial cables and galvanically coupled to the small Josephson junction (in red). a Optical picture of the input and output capacitive couplings. b SEM picture of a few of the SQUIDs (1500 in total for each chain) that are coupled to the small Josephson junction (in red). c Equivalence between the transmission line effective picture and the SQUID chain characterized by three microscopic parameters and the inductance and capacitance per SQUID respectively and the ground capacitance

Table 1

Parameters of three samples. The bare Josephson energy is inferred using the Ambegaokar–Baratoff law. is the measured value of the renormalized Josephson energy. As a consistency check, the bare value is also extracted from the fit of using the SCHA. is the capacitance shunting the small Josephson junction (see Supplementary Note 9). , , and are obtained from the dispersion relation of the chain (see Supplementary Note 10)

Sample	A	B	C
Small junction
Area [μm^2]	315 × 195	370 × 190	440 × 185
C_J [fF]	2.7 ± 0.3	3.2 ± 0.3	3.7 ± 0.4
C_sh [fF]	3.0 ± 0.5	2.4 ± 0.4	5.1 ± 1.0
EJ*[GHz]	1.8 ± 0.1	3.1 ± 0.2	5.7 ± 0.3
EJ,bareAB [GHz]	3.7 ± 0.2	5.8 ± 0.3	6.8 ± 0.5
EJ,bareth [GHz]	3.7	5.5	8.2
Nonlinearity E_J,bare∕E_c	0.27	0.40	0.93
Renormalization EJ*∕E_J,bare	0.49	0.56	0.70
Chain
C [fF]	144	144	144
C_g [fF]	0.189	0.192	0.181
L [nH]	0.66	0.60	0.61
E_J∕E_c	460	506	498

The broadband microwave transmission of the full system shows a series of resonances (see Fig. 2b). A broadening of the modes in the array is expected since the SQUID chains are capacitively coupled to the measurement lines, hence forming very long microwave resonators. The transmission of the system is measured using very low microwave power, down to the single-photon regime. This prevents any power-induced broadening or frequency shift of these resonances (see Supplementary Note 9). A closer look at Fig. 2b reveals that resonances come in pairs. This is expected given the symmetry of the sample: our system can be decomposed into two subsystems (See Supplementary Note 2). One is made of even modes, which are decoupled from the small Josephson junction, while the other is composed of odd modes, with impedance , ultrastrongly coupled to the small Josephson junction[18,48]. A more surprising observation is that the odd modes are much more damped than the even ones. We interpret this as resulting from the nonlinearity that odd modes inherit from the small Josephson junction. This is experimental evidence of the strong back-action of the small Josephson junction on the many modes of the chain forming its linear environment.

Fig. 2

Inferring the renormalized resonant frequency of the small Josephson junction. a Amplitude of the microwave transmission versus frequency (sample A, 24 mK). The even–odd modes frequency splitting changes sign precisely at . Arrows are guides to the eye of the splitting sign. b Fit of the double peaks for three cases: well below the resonant frequency of the small Josephson junction (blue) its inductive part dominates, close to (orange) the impedance of the junction is large so that the two modes are almost decoupled, and well above (green) the capacitive part of the junction dominates. c Experimental normalized frequency splittings obtained from the previous fits (dots) and theoretical prediction (full line). The resonance frequency of the small Josephson junction corresponds to the vanishing value of the normalized splitting

Line shapes

The line shape of a given even–odd pair of resonances can be obtained by associating with it two effective LC oscillator[49] connected via the small Josephson junction (see insets in Fig. 2b). In the regime of interest (See Supplementary Note 5 and 6) this junction can be treated as a ZPF-dependent inductance in parallel with a capacitance , with resonance frequency . The odd and even modes mentioned earlier are characterized by respective frequencies and , with and . Then for modes at frequencies such that , the capacitance of the small junction can be neglected () leading to . In the opposite case () the inductance can be neglected, giving . The most interesting regime is when the system is probed close to . In that case, the impedance of the small junction diverges and consequently the two effective oscillators are uncoupled leading to . In the Supplementary Note 7, we confirm that a fully microscopic model of the whole circuit also predicts that the frequency splitting between even and odd modes changes sign at the renormalized frequency of the junction . The frequencies of each even–odd pair of modes is extracted by fitting the peaks Fig. 2 to line shapes of an input–output formalism based on the simple model just described (see Supplementary Fig. 3 and “Methods” section).

Renormalized Josephson energy

The effective resonance frequency of the junction, , depends on its environment due to the interplay of strong anharmonicity and many-body ZPF, and can be inferred by tracking the evolution of the normalized frequency splitting , between even (uncoupled) and odd (coupled) modes, where refers to mode number. As shown in Puertas-Martìnez et al.[18], in a long chain, this quantity equals the phase shift difference between even and odd modes. It vanishes when the left and right halves of the device decouple, so that even and odd modes become degenerate. Figure 2c shows the experimentally obtained for one of our samples, from which we extract . As we show in the Supplementary Notes 4 and 5, the ZPF-dependent effective inductance of the weak link is related to a renormalized Josephson energy aswhere , with the intrinsic capacitance of the junction and a shunting capacitance due to the surrounding circuitry. Note that we define in terms of , and not in terms of the DC critical current as is done in for instance[23]. We use Eq. (1) to infer experimentally. is given by the junction size measured from an SEM picture. The way is extracted is explained in the Supplementary Note 9. Values for sample A, B, and C are reported in Table 1. To see the effect of vacuum fluctuations, we compare with the bare Josephson energy of the weak link, which was obtained as follows. We fabricated many nominally identical Josephson junctions on the same chip and measured their room temperature resistances. The expected bare Josephson energy of the small Josephson junction (see Table 1) was then inferred using the Ambegaokar–Baratoff law. We observe a systematic shift between this bare energy and the renormalized one we inferred from measurements, a shift that is more pronounced for sample A that shows a high nonlinearity. This points towards a large renormalization induced by the strong zero-point phase fluctuations of the hybridized junction-chain modes, as expected since the small junction is impedance matched to the chains.

We now show that this renormalization is quantitatively captured by a microscopic model based on the SCHA. Its success in accounting for nonlinearities introduced by Josephson junctions is well-established[13,14,50,51]. More recently it was employed in detailed microscopic models in the field of circuit QED[18,52]. The idea behind the SCHA is that the strong phase fluctuations allowed by the environment average the nonlinear potential of the small Josephson junction, lowering its effective Josephson energy from the bare value to the renormalized one . This is valid, provided the phase , though strongly fluctuating, is still sufficiently localized. In this regard we note the following. Though large, the effective environmental impedance seen by the weak link, is still less than . Under this condition, the environment is known to produce spontaneous symmetry breaking of the 2 periodicity in the phase difference across the weak link[19,20,22]. It is therefore reasonable to approximate the system’s full wave function with a Gaussian that is fairly well localized in the direction, which is the essence of the SCHA. At zero temperature, the interplay of many-body ZPF and nonlinearity can be described is approximated by replacing the cosine Josephson potential by an effective quadratic term , where the renormalized Josephson energy is given by the self-consistent equation:

Here, the total phase fluctuation across the junction is given bywhere is the contribution to the small junction ZPF coming from odd mode . Importantly, in the strong ZPF regime, the expectation value must be taken with respect to the modified vacuum of the hybridized modes, which means that the normal modes of the systems has to be updated during the numerical iteration of Eq. (2). This is in contrast to familiar examples of ZPF induced phenomena, such as the Lamb shift in hydrogen, where the perturbative nature of the effect allows one to calculate fluctuations with respect to the bare vacuum of the environment. We independently extracted the parameters of the whole circuit (junction + chains), and then used Eq. (2) to determine the theoretical bare Josephson energy required to find back the measured renormalized (see next section for more details). The agreement between experimentally and theoretically estimated (see Table 1) provides strong evidence that our system displays large ZPF, which leads to a renormalization of up to of the Josephson energy of the small junction (or equivalently of its resonant frequency ). Moreover, as expected, this renormalization increases when the ratio decreases, or equivalently when the nonlinearity of the small Josephson junction increases.

Quantum versus thermal fluctuations

As is renormalized by phase fluctuations across the weak link, one expects a crossover from quantum to thermally driven fluctuations as temperature increases. Extending the SCHA to nonzero temperatures (see Supplementary Note 4 and 5), we find that the fluctuations of mode contain a Bose factor contribution:with the frequency of mode , and its zero temperature ZPF. Therefore, at low temperature, fluctuations saturate to a finite ZPF value (a hallmark of quantum uncertainty), while at high temperature they increase linearly with temperature (4). According to Eq. (2), should decrease when the system is heated up. Consequently, odd modes’ frequencies are shifted to lower values when temperature increases, while the even modes stay put. This striking experimental signature of nonlinearity can clearly be seen in 3a. This constitutes smoking gun evidence of the back-action of the Josephson junction on its environment: the shift of to smaller values at increasing temperatures indicates that fluctuations are thermally enhanced.

Fig. 3

Temperature-induced renormalization. a Zoom on a even–odd pair of transmission peaks for sample A at temperature ranging from 23 to 150 mK. The even mode (gray) does not move while the odd mode (blue is at 25 mK, red at 130 mK) shifts down in frequency when warming up, showing a downward renormalization of the junction frequency . b ZPF of the small junction as a function of the temperature for three samples (A, B, and C ranging from dark to light blue), extracted from Eq. (2). ZPF are stronger in sample A, which is associated to a smaller ratio (large nonlinearity). The measured quantum to classical crossover is in good agreement with theory (full lines). The inset displays the corresponding renormalized junction frequency of the three samples. The full lines are the SCHA predictions while the dashed lines represent what would be the temperature evolution of these frequencies if ZPF were omitted from , using the same values of

The recipe to extract is the following: is obtained from measurements at different temperatures. Since all the other parameters (, , , , and ) are known, we can fit using Eqs. (2) and (4), taking as the (only) fitting parameter. Then, being determined, we can compute the phase fluctuations across the small Josephson junction using Eq. (1) and (2):We checked that at the lowest temperature of our cryostat, the phase fluctuations experienced by the small Josephson junction are fully in the quantum regime, by measuring from 25 to 130 mK. Results are shown in Fig. 3b. We observe that the quantum to classical crossover appears at decreasing temperatures from sample C to A. This is because decreases from sample C to A. Therefore the junction is coupled to modes with lower and lower frequencies, which are thermally occupied at lower temperatures. The inset in Fig. 3b shows the corresponding fit of for the three samples. The dashed lines represent obtained using the value of extracted from the previous fit but including only thermal renormalization of i.e. disregarding ZPF. Consequently, is given by:The discrepancy between the dashed lines and the fit clearly shows that the fluctuations have mainly quantum origin. At increasing temperatures, thermal fluctuations add to the quantum ZPF, and cause a rise in , witnessed both in the experimental extraction and the predictions from SCHA, see Fig. 3b. It is likely that the extracted for sample A is systematically underestimated due to sizeable errors in the SCHA that rapidly set in after , leading to a mismatch with the theory at high temperatures.

Many-body nature of the ZPF

In order to confirm the many-body character of this renormalization, we can estimate how many modes are affecting the small junction simultaneously. The ZPF are quantitatively determined by how the full vacuum of the whole circuit is dressed by the coupling through the weak link. Within the SCHA, the number of modes contributing a finite amount of provides a measure of the number of interacting particles in the system. In Fig. 4 we compare the experimentally extracted with various calculated values. In each calculation, the full system was truncated to a finite number of modes in a window around . If the window is too narrow, important contributions to the ZPF are neglected, and is underestimated. The comparison unambiguously shows that, in sample B, around 30 modes contribute to the total phase fluctuations. In circuit-QED language, the full width at half maximum (FWHM) of the environmental ZPF – labeled – is about 7 GHz for our samples (see inset of Fig. 4). Therefore, our device operates in a regime where due to the impedance matching to the transmission line. Moreover, our device is strongly nonlinear. Consequently, it is not possible to treat perturbatively the nonlinearity as is usually done in the field for the Transmon qubit or other weakly nonlinear circuits[5,15-17] (a detailed analysis is given in Supplementary Note 11).

Fig. 4

Many-body nature of the ZPF. Total phase fluctuations across the small Josephson junction in sample B, taking into account in our model (full lines) different numbers of modes of the environment, ranging from one (light blue) to the total number (dark blue). The inset shows the relative contribution of the different modes to the total fluctuations, with being the FWHM of this quantity

This work provides a direct observation of several quantum many-body effects driven by ZPF in an open quantum system. This was achieved by developing a spectroscopic setup where the high-impedance environment of a single nonlinear Josephson junction was monitored mode by mode, and compared with a detailed microscopic model. A strong quantum renormalization (up to 50%) of the Josephson energy of the single junction was demonstrated, analogous to a nonperturbative Lamb shift. In addition, the back-action of the small Josephson junction causes nonlinear broadening and strong temperature dependence of the environmental modes, providing the most striking signature of the many-body effects that take place in our circuit. The measured temperature dependence of the phase fluctuation across the Josephson junction indicates that our device remains quantum coherent at cryogenic temperatures. As many as 30 modes are involved in the renormalization of the small junction. Our superconducting circuit thus behaves as a fully fledged quantum many body simulator, paving the way for the further observation of various many-body nonlinear effects in circuit QED[53-59].

Methods

Full model

The Hamiltonian of the full system can be decomposed into odd and even parts—containing respectively the modes coupled and not coupled to the junction (see Supplementary Note 2). The odd Hamiltonian readswith and referring to the charge the phase drop across the small junction while and , refer to the charge and phase operators on chain site . Charge and phase operators obey the commutation rules . The microscopic parameters are , the Josephson energy of the SQUIDs, the bare Josephson energy of the small junction, and the capacitance matrix:with:

Self-consistent harmonic approximation (SCHA)

Because of the cosine term in Eq. (7), we are dealing with an interacting many-body problem that cannot be solved analytically. To study the best variational harmonic approximation we use the SCHA:with the trial harmonic Hamiltonian that will approximate , optimized with respect to the renormalized Josephson energy . The variational principle giveswith the many-body ground state of . Because of the harmonic character of , we haveInserting (14) into (13) we end up with the self-consistent equation:The physical interpretation is the following: when ZPF are negligible, and , so in its ground state the junction behaves as an harmonic oscillator of frequency . For weak nonlinearity, fluctuations increase but remain such that , resulting in , so that the junction behaves as a weakly anharmonic oscillator with fundamental frequency . For an isolated junction , and the frequency becomes , a well-known result for the Transmon qubit[5]. For larger fluctuations, the principle remains the same but no analytical formula can be derived, so that one should solve the self-consistent equation numerically. A more detailed derivation—including thermal fluctuations—is presented in Supplementary Notes 4 and 5.

Frequency splitting S between odd and even modes

The splitting is linked to the phase shift difference between even and odd modes in the thermodynamic limit[18]:The analytical formula of the phase shift difference—derived in the Supplementary 6 and 7 readswith being the plasma frequency of the chain and the effective inductance of the small junction.

Fitting formula for the peaks

Input–output theory is used to fit the parameters associated to the double resonances observed in the transmission spectrum. These are mapped to two coupled harmonic modes and with mutual coupling rate , external coupling and internal loss , with Hamiltonian:Here ,  are the left input and output signals and is the right output signal. Thus the input–output relations are :The equations of motion are:The complex transmission is defined as , and can be calculated using Eqs. (21)–(24). We define the even and odd frequencies, and add phenomenologically losses in the odd modes (we keep ), so that:For some of the odd modes, we found a signature of inhomogeneous broadening, that we modeled by a convolution of their frequency with a gate function defined as if . Understanding microscopically this additional broadening, possibly due to offset charges, is beyond the scope of the description using the SCHA, and will require additional theoretical developments.