Abstract
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.
Similar content being viewed by others
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 platforms1,2,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 proposed4 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 \(\hat{n}\) and the superconducting phase difference \(\hat{\phi }\). Despite being an anharmonic oscillator, a Josephson junction with Josephson energy \({E}_{{\rm{J}}}\) and charging energy \({E}_{{\rm{c}}}\), can be associated to an impedance \({Z}_{{\rm{J}}}\,=\,{\hbar} /{(2e)}^{2}\sqrt{2{E}_{{\rm{c}}}/{E}_{{\rm{J}}}},\) which sets the amplitude of the fluctuations of \(\hat{n}\) and \(\hat{\phi }\). When \(\langle {\hat{\phi }}^{2}\rangle\) is sufficiently smaller than unity, \(\langle {\hat{\phi }}^{2}\rangle\,\propto\,{Z}_{{\rm{J}}}/{R}_{{\rm{Q}}}\) and \(\langle {\hat{n}}^{2}\rangle \propto {R}_{{\rm{Q}}}/{Z}_{{\rm{J}}}\), with \({R}_{{\rm{Q}}}=h/{(2e)}^{2}\simeq \ 6.5\ {\rm{k}}\Omega\) the superconducting quantum of resistance. Consequently, at low \({Z}_{{\rm{J}}}\), phase fluctuations are weak and the anharmonic Josephson cosine potential \({E}_{{\rm{J}}}(1-\cos \hat{\phi })\) can be reduced to a quadratic potential plus a quartic perturbation, as is the case for the transmon qubit5. On the other hand, if \({Z}_{{\rm{J}}}\) is large, the full cosine potential is explored due to strong phase fluctuations. Anharmonicity then becomes important, as observed with the Cooper-pair box6 or the fluxonium qubit7, and as a result, the oscillation frequency \({\omega }_{{\rm{J}}}\) can strongly deviate from the harmonic value \(\sqrt{2{E}_{{\rm{J}}}{E}_{{\rm{C}}}}\). Thus, exploring many-body physics in circuit quantum electrodynamics must rely on a careful tailoring of ZPF.
The approach8 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 \({Z}_{{\rm{J}}}\) does not match the environmental impedance \({Z}_{{\rm{env}}}(\omega )\) at frequencies \(\omega\) close to \({\omega }_{{\rm{J}}}\), 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 shift9. 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 approach10. 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 tunneling11. The important role of ZPF in the damping effect that such an environment has on a Josephson junction was already noticed experimentally12 and explained theoretically13,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 \(1 \%\), which nonetheless had a large effect on macroscopic quantum tunneling rates. When \({Z}_{{\rm{env}}} \sim {Z}_{{\rm{J}}}\ll {R}_{{\rm{Q}}}\), the junctions and its environment fully hybridize, since they are impedance matched, but the anharmonicity of the junction remains weak and can be treated pertubatively15,16,17,18. The case \({Z}_{{\rm{env}}} \sim {Z}_{{\rm{J}}} \sim {R}_{{\rm{Q}}}\) 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 ago19,20, a localization phenomenon whose relevance for microwave AC measurements requires further experimental and theoretical investigations21. 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 matter22. More specifically, the frequency shift induced by the environment on the junction can be comparable with \({\omega }_{{\rm{J}}}\), a nonperturbative effect due to a modification of the vacuum23. 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 \({Z}_{{\rm{env}}} \sim {R}_{{\rm{Q}}}\) at gigahertz frequencies is very challenging. One option is to use on-chip resistors24. However, this may lead to unwanted Joule heating25. Therefore, we rather pursue a solution that relies on superconducting (lossless) high-inductance materials such as Josephson junction arrays26,27,28, noting that disordered superconductors29 are also promising. In Josephson junction arrays, \({Z}_{{\rm{env}}}=\sqrt{L/C}\) can reach \({R}_{{\rm{Q}}}\) 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 pairs24 in the framework of the \(P(E)\) theory30. 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 \({Z}_{{\rm{J}}}\), \({Z}_{{\rm{env}}}\) Â >Â \({R}_{{\rm{Q}}}\) was explored using SQUID arrays as the environment31,32,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 \({\omega }_{{\rm{J}}}\)â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 QED34. Here, microwave techniques allow a more accurate examination of the effects of ZPF on Josephson junctions35, and observations of perturbative spectral shifts (below 1%) attributed to ZPF were reported9,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 resonator40. Moving towards many-body territory, a nonperturbative renormalization of the frequency of a flux qubit was demonstrated41,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 authors43,44,45,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 arrays18,48. It consists of a small Josephson junction of characteristic impedance on the order of \({R}_{{\rm{Q}}}\) (\({E}_{{\rm{J}}}/{E}_{{\rm{c}}}\ \lesssim \ 1\)), 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â\({\omega }_{{\rm{plasma}}}\ \simeq\) 17âGHz and its wave impedance \({Z}_{{\rm{chain}}}=\sqrt{L/{C}_{{\rm{g}}}}\ \simeq \ 1.8\ {\rm{k}}\Omega\) (see Table 1 and Supplementary Note 10). The SQUID parameters were carefully chosen to maintain a negligible phase slip rate (\({E}_{{\rm{J}}}/{E}_{{\rm{c}}}\ \lesssim \ 500\)), 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).
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 \(L\) and \(C\) the inductance and capacitance per SQUID respectively and \({C}_{g}\) the ground capacitance
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 \(50\ \Omega\) 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 \({Z}_{{\rm{env}}}=2{Z}_{{\rm{chain}}}\), ultrastrongly coupled to the small Josephson junction18,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.
Inferring the renormalized resonant frequency \({\omega }_{{\rm{J}}}^{* }\) of the small Josephson junction. a Amplitude of the microwave transmission \(| {S}_{21}|\) versus frequency (sample A, 24âmK). The evenâodd modes frequency splitting \(S\) changes sign precisely at \({\omega }_{{\rm{J}}}^{* }\). 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 \({\omega }_{{\rm{J}}}^{* }\) (orange) the impedance of the junction is large so that the two modes are almost decoupled, and well above \({\omega }_{{\rm{J}}}^{* }\) (green) the capacitive part of the junction dominates. c Experimental normalized frequency splittings \(S\) obtained from the previous fits (dots) and theoretical prediction (full line). The resonance frequency \({\omega }_{{\rm{J}}}^{* }\) of the small Josephson junction corresponds to the vanishing value of the normalized splitting \(S\)
Line shapes
The line shape of a given evenâodd pair of resonances can be obtained by associating with it two effective LC oscillator49 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 \({L}_{{\rm{J}}}^{* }\) in parallel with a capacitance \({C}_{\parallel }\), with resonance frequency \({\omega }_{{\rm{J}}}^{* }=1/\sqrt{{L}_{{\rm{J}}}^{* }{C}_{\parallel }}\). The odd and even modes mentioned earlier are characterized by respective frequencies \({\omega }_{{\rm{even}}}=1/\sqrt{LC}\) and \({\omega }_{{\rm{odd}}}=1/\sqrt{{L}_{\Sigma }{C}_{\Sigma }}\), with \(1/{L}_{\Sigma }=1/2L+1/{L}_{{\rm{J}}}^{* }\) and \({C}_{\Sigma }=C/2+{C}_{\parallel }\). Then for modes at frequencies such that \({\omega }_{{\rm{odd}}},{\omega }_{{\rm{even}}}\ll {\omega }_{{\rm{J}}}^{* }\), the capacitance of the small junction can be neglected (\({C}_{\parallel } \sim 0\)) leading to \({\omega }_{{\rm{odd}}}\ > \ {\omega }_{{\rm{even}}}\). In the opposite case (\({\omega }_{{\rm{odd}}},{\omega }_{{\rm{even}}}\gg {\omega }_{{\rm{J}}}^{* }\)) the inductance can be neglected, giving \({\omega }_{{\rm{odd}}}\ <\ {\omega }_{{\rm{even}}}\). The most interesting regime is when the system is probed close to \({\omega }_{{\rm{J}}}^{* }\). In that case, the impedance of the small junction diverges and consequently the two effective oscillators are uncoupled leading to \({\omega }_{{\rm{odd}}}={\omega }_{{\rm{even}}}\). 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 \({\omega }_{{\rm{J}}}^{* }\). 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 \({E}_{{\rm{J}}}^{* }\)
The effective resonance frequency of the junction, \({\omega }_{{\rm{J}}}^{* }\), 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 \(S=({\omega }_{{\rm{odd,k}}}-{\omega }_{{\rm{even,k}}})/({\omega }_{{\rm{even}},{\rm{k}}+1}-{\omega }_{{\rm{even,k}}})\), between even (uncoupled) and odd (coupled) modes, where \({\rm{k}}=0\ldots M\) 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 \(S\) for one of our samples, from which we extract \({\omega }_{{\rm{J}}}^{* }\). 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 \({E}_{{\rm{J}}}^{* }={({\hbar} /2e)}^{2}/{L}_{{\rm{J}}}^{* }\) as
where \({E}_{{\rm{c}}}={(2e)}^{2}/(2({C}_{{\rm{J}}}+{C}_{{\rm{sh}}}))\), with \({C}_{{\rm{J}}}\) the intrinsic capacitance of the junction and \({C}_{{\rm{sh}}}\) a shunting capacitance due to the surrounding circuitry. Note that we define \({E}_{{\rm{J}}}^{* }\) in terms of \({L}_{{\rm{J}}}^{* }\), and not in terms of the DC critical current as is done in for instance23. We use Eq. (1) to infer \({E}_{{\rm{J}}}^{* }\) experimentally. \({C}_{{\rm{J}}}\) is given by the junction size measured from an SEM picture. The way \({C}_{{\rm{sh}}}\) 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 \({E}_{{\rm{J}}}^{* }\) 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 \({E}_{{\rm{J,bare}}}^{{\rm{AB}}}\) (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 \(| {S}_{21}|\) 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-established13,14,50,51. More recently it was employed in detailed microscopic models in the field of circuit QED18,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 \({E}_{{\rm{J}}}\) to the renormalized one \({E}_{{\rm{J}}}^{* }\). This is valid, provided the phase \({\phi }_{{\rm{J}}}\), though strongly fluctuating, is still sufficiently localized. In this regard we note the following. Though large, the effective environmental impedance \(2{Z}_{{\rm{chain}}}\simeq \ 3.8\ {\rm{k}}\Omega\) seen by the weak link, is still less than \({R}_{{\rm{Q}}}\). Under this condition, the environment is known to produce spontaneous symmetry breaking of the 2\(\pi\) periodicity in the phase difference \({\phi }_{{\rm{J}}}\) across the weak link19,20,22. It is therefore reasonable to approximate the systemâs full wave function with a Gaussian that is fairly well localized in the \({\phi }_{{\rm{J}}}\) 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 \({E}_{{\rm{J}}}^{* }{\phi }_{{\rm{J}}}^{2}/2\), where the renormalized Josephson energy \({E}_{{\rm{J}}}^{* }\) is given by the self-consistent equation:
Here, the total phase fluctuation across the junction \(\langle {\phi }_{{\rm{J}}}^{2}\rangle\) is given by
where \({\phi }_{k}^{2}\) is the contribution to the small junction ZPF coming from odd mode \(k\). 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 \({E}_{{\rm{J}}}^{* }\) (see next section for more details). The agreement between experimentally and theoretically estimated \({E}_{{\rm{J,bare}}}\) (see Table 1) provides strong evidence that our system displays large ZPF, which leads to a renormalization of up to \(50 \%\) of the Josephson energy of the small junction (or equivalently \(30 \%\) of its resonant frequency \({\omega }_{{\rm{J}}}\)). Moreover, as expected, this renormalization increases when the ratio \({E}_{{\rm{J,bare}}}/{E}_{{\rm{c}}}\) decreases, or equivalently when the nonlinearity of the small Josephson junction increases.
Quantum versus thermal fluctuations
As \({\omega }_{{\rm{J}}}^{* }\) 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 \(k\) contain a Bose factor contribution:
with \({\omega }_{k}\) the frequency of mode \(k\), and \({\phi }_{k}^{2}\) 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), \({\omega }_{{\rm{J}}}^{* }\) 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 \({\omega }_{{\rm{J}}}^{* }\) to smaller values at increasing temperatures indicates that fluctuations are thermally enhanced.
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 \({\omega }_{{\rm{J}}}^{* }\). b ZPF of the small junction \(\langle {\phi }_{{\rm{J}}}^{2}\rangle\) 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 \({E}_{{\rm{J,bare}}}/{E}_{{\rm{c}}}\) (large nonlinearity). The measured quantum to classical crossover is in good agreement with theory (full lines). The inset displays the corresponding renormalized junction frequency \({f}_{{\rm{J}}}^{* }={\omega }_{{\rm{J}}}^{* }/2\pi\) 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 \(\langle {\phi }_{{\rm{J}}}^{2}\rangle\), using the same values of \({E}_{{\rm{J,bare}}}\)
The recipe to extract \({E}_{{\rm{J,bare}}}^{{\rm{th}}}\) is the following: \({E}_{{\rm{J}}}^{* }(T)\) is obtained from \({S}_{21}\) measurements at different temperatures. Since all the other parameters (\(L\), \(C\), \({C}_{{\rm{g}}}\), \({C}_{{\rm{J}}}\), and \({C}_{{\rm{sh}}}\)) are known, we can fit \({E}_{{\rm{J}}}^{* }(T)\) using Eqs. (2) and (4), taking \({E}_{{\rm{J,bare}}}^{{\rm{th}}}\) as the (only) fitting parameter. Then, \({E}_{{\rm{J,bare}}}^{{\rm{th}}}\) 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 \(| {S}_{21}|\) 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 \({\omega }_{{\rm{J}}}^{* }\) 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 \({\omega }_{{\rm{J}}}^{* }\) for the three samples. The dashed lines represent \({\omega }_{{\rm{J}}}^{* }\) obtained using the value of \({E}_{{\rm{J,bare}}}\) extracted from the previous fit but including only thermal renormalization of \({\omega }_{{\rm{J}}}^{* }\) i.e. disregarding ZPF. Consequently, \(\langle {\phi }_{{\rm{J}}}^{2}\rangle\) 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 \(\langle {\phi }_{{\rm{J}}}^{2}\rangle\), witnessed both in the experimental extraction and the predictions from SCHA, see Fig. 3b. It is likely that the extracted \(\langle {\phi }_{{\rm{J}}}^{2}\rangle\) for sample A is systematically underestimated due to sizeable errors in the SCHA that rapidly set in after \(\langle {\phi }_{{\rm{J}}}^{2}\rangle \ \gtrsim \ 1\), 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 \({\phi }_{k}^{2}\) provides a measure of the number of interacting particles in the system. In Fig. 4 we compare the experimentally extracted \(\langle {\phi }_{{\rm{J}}}^{2}(T)\rangle\) with various calculated values. In each calculation, the full system was truncated to a finite number of modes in a window around \({\omega }_{{\rm{J}}}^{* }\). If the window is too narrow, important contributions to the ZPF are neglected, and \(\langle {\phi }_{{\rm{J}}}^{2}(T)\rangle\) 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 \({\phi }_{{\rm{k}}}^{2}({\omega }_{{\rm{k}}})\) â labeled \({\Gamma }_{{\rm{J}}}\) â is about 7âGHz for our samples (see inset of Fig. 4). Therefore, our device operates in a regime where \({\Gamma }_{{\rm{J}}}/{\omega }_{{\rm{J}}}^{* } \sim 1\) 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 circuits5,15,16,17 (a detailed analysis is given in Supplementary Note 11).
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 \({\Gamma }_{{\rm{J}}}\) 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 QED53,54,55,56,57,58,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 reads
with \({\hat{n}}_{0}\equiv {\hat{n}}_{J}\) and \({\hat{\phi }}_{{\rm{J}}}\) referring to the charge the phase drop across the small junction while \({\hat{n}}_{i}\) and \({\hat{\phi }}_{i}\), \(i\in [1..N]\) refer to the charge and phase operators on chain site \(i\in [1..N]\). Charge and phase operators obey the commutation rules \(\left[\right.{\hat{\phi }}_{k},{\hat{n}}_{p}\left]\right.=i{\delta }_{k,p}\). The microscopic parameters are \({E}_{{\rm{J,S}}}\), the Josephson energy of the SQUIDs, \({E}_{{\rm{J,bare}}}\) 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 \({\hat{H}}_{{\rm{t}}}\) the trial harmonic Hamiltonian that will approximate \(\hat{H}\), optimized with respect to the renormalized Josephson energy \({E}_{{\rm{J}}}^{* }\). The variational principle gives
with \(\left|{\Psi }_{{\rm{t}}}\right\rangle\) the many-body ground state of \({\hat{H}}_{{\rm{t}}}\). Because of the harmonic character of \({\hat{H}}_{{\rm{t}}}\), we have
Inserting (14) into (13) we end up with the self-consistent equation:
The physical interpretation is the following: when ZPF are negligible, \(\langle {\hat{\phi }}_{{\rm{J}}}^{2}\rangle \ \simeq \ 0\) and \({E}_{{\rm{J}}}^{* }={E}_{{\rm{J,bare}}}\), so in its ground state the junction behaves as an harmonic oscillator of frequency \({\omega }_{{\rm{J,bare}}}=\sqrt{2{E}_{{\rm{J,bare}}}{E}_{{\rm{c}}}}\). For weak nonlinearity, fluctuations increase but remain such that \(\langle {\hat{\phi }}_{{\rm{J}}}^{2}\rangle \ll 1\), resulting in \({E}_{{\rm{J}}}^{* }\simeq {E}_{{\rm{J,bare}}}(1-\langle {\phi }_{{\rm{J}}}^{2}\rangle /2)\), so that the junction behaves as a weakly anharmonic oscillator with fundamental frequency \({\omega }_{{\rm{J,bare}}}(1-\langle {\hat{\phi }}_{{\rm{J}}}^{2}\rangle /4)\). For an isolated junction \(\langle {\hat{\phi }}_{{\rm{J}}}^{2}\rangle =\sqrt{{E}_{{\rm{c}}}/2{E}_{{\rm{J,bare}}}}\), and the frequency becomes \({\omega }_{{\rm{J,bare}}}-{E}_{{\rm{c}}}/4\), a well-known result for the Transmon qubit5. 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 \(S\) is linked to the phase shift difference \(\theta\) between even and odd modes in the thermodynamic limit18:
The analytical formula of the phase shift differenceâderived in the Supplementary 6 and 7 reads
with
\({\omega }_{{\rm{p}}}=1/\sqrt{L(C+{C}_{{\rm{g}}}/4)}\) being the plasma frequency of the chain and \({L}_{{\rm{J}}}^{* }={\hbar}^{2}/{(2e)}^{2}{E}_{{\rm{J}}}^{* }\) 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 \(\alpha\) and \(\beta\) with mutual coupling rate \(g\), external coupling \({\kappa }_{{\rm{ext}}}\) and internal loss \({\kappa }_{{\rm{in}}}\), with Hamiltonian:
Here \({\hat{a}}_{{{\rm{in}}}_{{\rm{L}}}}\), \({\hat{a}}_{{{\rm{out}}}_{{\rm{L}}}}\) are the left input and output signals and \({\hat{a}}_{{{\rm{out}}}_{{\rm{R}}}}\) is the right output signal. Thus the inputâoutput relations are :
The equations of motion are:
The complex transmission is defined as \({S}_{21}={\hat{a}}_{{{\rm{out}}}_{{\rm{R}}}}/{\hat{a}}_{{{\rm{in}}}_{{\rm{L}}}}\), and can be calculated using Eqs. (21)â(24). We define the even \({\omega }_{{\rm{e}}}={\omega }_{{\rm{r}}}+g\) and odd \({\omega }_{{\rm{o}}}={\omega }_{{\rm{r}}}-g\) frequencies, and add phenomenologically losses in the odd modes \({\kappa }_{{\rm{o}}}={\kappa }_{{\rm{in}}}+{\kappa }_{{\rm{add}}}\) (we keep \({\kappa }_{{\rm{e}}}={\kappa }_{{\rm{in}}}\)), 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 \({\Pi }_{\delta \omega }(\omega )=1/\delta \omega\) if \(\omega \in [\omega -\delta \omega /2,\omega +\delta \omega /2]\). 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.
Data availability
The data that support the findings of this study as well the treatment scripts are available at the open access repository Zenodo [https://doi.org/10.5281/zenodo.3520349].
References
Greentree, A. D., Tahan, C., Cole, J. H. & Hollenberg, L. C. L. Quantum phase transitions of light. Nat. Phys. 2, 856â861 (2006).
Carusotto, I. & Ciuti, C. Quantum fluids of light. Rev. Mod. Phys. 85, 299â366 (2013).
Hur, K. L. et al. Many-body quantum electrodynamics networks: non-equilibrium condensed matter physics with light. Comptes Rendus Phys. 17, 808â835 (2016).
Houck, A. A., Türeci, H. E. & Koch, J. On-chip quantum simulation with superconducting circuits. Nat. Phys. 8, 292â299 (2012).
Koch, J. et al. Charge-insensitive qubit design derived from the Cooper pair box. Phys. Rev. A 76, 042319 (2007).
Vion, D. et al. Manipulating the quantum state of an electrical circuit. Science 296, 886â889 (2002).
Manucharyan, V. E., Koch, J., Glazman, L. I. & Devoret, M. H. Fluxonium: single cooper-pair circuit free of charge offsets. Science 326, 113â116 (2009).
U. Weiss, Quantum Dissipative Systems (World Scientific., 1992).
Fragner, A. et al. Resolving vacuum fluctuations in an electrical circuit by measuring the lamb shift. Science 322, 1357â1360 (2008).
Leggett, A. J. et al. Dynamics of the dissipative two-state system. Rev. Mod. Phys. 59, 1â85 (1987).
Clarke, J., Cleland, A. N., Devoret, M. H., Esteve, D. & Martinis, J. M. Quantum-mechanics of a macroscopic variableâthe phase difference of a Josephson junction. Science 239, 992â997 (1988).
Schwartz, D. B., Sen, B., Archie, C. N. & Lukens, J. E. Quantitative study of the effect of the environment on macroscopic quantum tunneling. Phys. Rev. Lett. 55, 1547â1550 (1985).
Zaikin, A. D. & Panyukov, S. V. Lifetime of macroscopic current states. JETP Lett. 43, 670 (1986).
Panyukov, S. V. & Zaikin, A. D. Quantum fluctuations and the current-phase relation in Josephson junctions and squids. Phys. B 152, 162 (1988).
Nigg, S. E. et al. Black-box superconducting circuit quantization. Phys. Rev. Lett. 108, 240502 (2012).
Bourassa, J., Beaudoin, F., Gambetta, J. M. & Blais, A. Josephson-junction-embedded transmission-line resonators: From Kerr medium to in-line transmon. Phys. Rev. A 86, 013814 (2012).
Weissl, T. et al. Kerr coefficients of plasma resonances in Josephson junction chains. Phys. Rev. B 92, 104508 (2015).
Puertas MartÃnez, J. et al. A tunable Josephson platform to explore many-body quantum optics in circuit-QED. npj Quantum Inf. 5, 19 (2019).
Schmid, A. Diffusion and localization in a dissipative quantum system. Phys. Rev. Lett. 51, 1506â1509 (1983).
Bulgadaev, S. A. Phase diagram of a dissipative quantum system. JETP Lett. 39, 315â319 (1984).
A. Murani, et al. âAbsence of a dissipative quantum phase transition in Josephson junctions.â Preprint at https://arxiv.org/abs/1905.01161 (2019).
Schön, G. & Zaikin, A. D. Quantum coherent effects, phase transitions, and the dissipative dynamics of ultra small tunnel junctions. Phys. Rep. 198, 237â412 (1990).
Hekking, F. W. J. & Glazman, L. I. Quantum fluctuations in the equilibrium state of a thin superconducting loop. Phys. Rev. B 55, 6551â6558 (1997).
Kuzmin, L. S., Nazarov, Y. V., Haviland, D. B., Delsing, P. & Claeson, T. Coulomb blockade and incoherent tunneling of Cooper pairs in ultrasmall junctions affected by strong quantum fluctuations. Phys. Rev. Lett. 67, 1161â1164 (1991).
Huard, B., Pothier, H., Esteve, D. & Nagaev, K. E. Electron heating in metallic resistors at sub-Kelvin temperature. Phys. Rev. B 76, 165426 (2007).
V. E. Manucharyan, âSuperinductance,â Thesis, Yale Univ. (2012).
Masluk, N., Pop, I., Kamal, A., Minev, Z. & Devoret, M. Microwave characterization of Josephson junction arrays: implementing a low loss superinductance. Phys. Rev. Lett. 109, 137002 (2012).
Bell, M., Sadovskyy, I., Ioffe, L., Kitaev, A. & Gershenson, M. Quantum superinductor with tunable nonlinearity. Phys. Rev. Lett. 109, 137003 (2012).
Maleeva, N. et al. Circuit quantum electrodynamics of granular aluminum resonators. Nat. Commun. 9, 3889 (2018).
G. L. Ingold and Y. V. Nazarov, âCharge Tunneling Rates in Ultrasmall Junctions,â in Single Charge Tunneling (Springer, Boston, MA, Boston, MA, 1992) pp. 21â107.
Corlevi, S., Guichard, W., Hekking, F. W. J. & Haviland, D. B. Phase-charge duality of a Josephson junction in a fluctuating electromagnetic environment. Phys. Rev. Lett. 97, 096802 (2006).
Ergül, A. et al. Localizing quantum phase slips in one-dimensional Josephson junction chains. New J Phys. 15, 095014 (2013).
Weissl, T. et al. Bloch band dynamics of a Josephson junction in an inductive environment. Phys. Rev. B 91, 014507 (2015).
Wallraff, A. et al. Strong coupling of a single photon to a superconducting qubit using circuit quantum electrodynamics. Nature 431, 162â166 (2004).
Hoi, I. C. et al. Probing the quantum vacuum with an artificial atom in front of a mirror. Nat. Phys. 11, 1045â1049 (2015).
Silveri, M. et al. Broadband Lamb shift in an engineered quantum system. Nat. Phys. 15, 533â537 (2019).
P. Y. Wen, et al. âLarge collective Lamb shift of two distant superconducting artificial atoms.â Preprint at https://arxiv.org/abs/1904.12473 (2019).
P. Forn-DÃaz, L. Lamata, E. Rico, J. Kono, and E. Solano. Ultrastrong coupling regimes of lightâmatter interaction. Reviews of Modern Physics 91, 025005 (2019).
Kockum, A. F., Miranowicz, A., Liberato, S. D., Savasta, S. & Nori, F. Ultrastrong coupling between light and matter. Nat. Rev. Phys. 1, 19â40 (2019).
Rolland, C. et al. Antibunched photons emitted by a dc-biased Josephson junction. Phys. Rev. Lett. 122, 186804 (2019).
Forn-DÃaz, P. et al. Ultrastrong coupling of a single artificial atom to an electromagnetic continuum in the nonperturbative regime. Nat. Phys. 13, 39â43 (2017).
Magazzù, L. et al. Probing the strongly driven spin-boson model in a superconducting quantum circuit. Nat. Commun. 9, 1403 (2018).
Garcia-Ripoll, J. J., Peropadre, B. & Liberato, S. D. Lightâmatter decoupling and \({A}^{2}\) term detection in superconducting circuits. Sci. Rep. 5, 16055 (2015).
Malekakhlagh, M., Petrescu, A. & Türeci, H. E. Cutoff-free circuit quantum electrodynamics. Phys. Rev. Lett. 119, 073601 (2017).
Gely, M. F. et al. Convergence of the multimode quantum Rabi model of circuit quantum electrodynamics. Phys. Rev. B 95, 245115 (2017).
Parra-Rodriguez, A., Rico, E., Solano, E. & Egusquiza, I. L. Quantum networks in divergence-free circuit QED. Quantum Sci. Technol. 3, 024012 (2018).
Gely, M. F., Steele, G. A. & Bothner, D. Nature of the Lamb shift in weakly anharmonic atoms: from normal-mode splitting to quantum fluctuations. Phys. Rev. A 98, 053808 (2018).
Kuzmin, R., Mehta, N., Grabon, N. & Mencia, R. Superstrong coupling in circuit quantum electrodynamics. npj Quantum Inf. 5, 20 (2019).
D. M. Pozar, Microwave engineering (John Wiley & Sons, 2009).
Kampf, A. & Schön, G. Quantum effects and the dissipation by quasiparticle tunneling in arrays of Josephson junctions. Phys. Rev. B 36, 3651â3660 (1987).
Chakravarty, S., Ingold, G. L., Kivelson, S. & Zimanyi, G. Quantum statistical mechanics of an array of resistively shunted Josephson junctions. Phys. Rev. B 37, 3283â3294 (1988).
Joyez, P. Self-consistent dynamics of a Josephson junction in the presence of an arbitrary environment. Phys. Rev. Lett. 110, 217003 (2013).
Garcia-Ripoll, J. J., Solano, E. & Martin-Delgado, M. A. Quantum simulation of Anderson and Kondo lattices with superconducting qubits. Phys. Rev. B 77, 024522 (2008).
Hur, K. L. Kondo resonance of a microwave photon. Phys. Rev. B 85, 140506 (2012).
Goldstein, M., Devoret, M. H., Houzet, M. & Glazman, L. I. Inelastic microwave photon scattering off a quantum impurity in a Josephson-junction array. Phys. Rev. Lett. 110, 017002 (2013).
Peropadre, B., Zueco, D., Porras, D. & GarcÃa-Ripoll, J. Nonequilibrium and nonperturbative dynamics of ultrastrong coupling in open lines. Phys. Rev. Lett. 111, 243602 (2013).
Snyman, I. & Florens, S. Robust Josephson-Kondo screening cloud in circuit quantum electrodynamics. Phys. Rev. B 92, 085131 (2015).
Gheeraert, N., Bera, S. & Florens, S. Spontaneous emission of Schrödinger cats in a waveguide at ultrastrong coupling. New J Phys. 19, 023036 (2017).
Gheeraert, N. et al. Particle production in ultrastrong-coupling waveguide QED. Phys. Rew. A 98, 043816 (2018).
Acknowledgements
The authors would like to thank F. Balestro, L. Del Rey, D. Dufeu, E. Eyraud, J. Jarreau, T. Meunier, and W. Wernsdorfer, for early support with the experimental setup. Very fruitful discussions with K. R. Amin, P. Forn-Diaz, J.-J. Garcia-Ripoll, D. B. Haviland, M. Houzet, P. Joyez, V. E. Manucharyan, F. Portier, and H. E. Tureci are acknowledged. The sample was fabricated in the Nanofab clean room. This research was supported by the ANR under contracts CLOUD (project number ANR-16-CE24-0005), GEARED (project number ANR-14-CE26-0018), by the National Research Foundation of South Africa (Grant No. 90657), and by the PICS contract FERMICATS. J.P.M. acknowledges support from the Laboratoire dâexcellence LANEF in Grenoble (ANR-10-LABX-51-01). R.D. and S.L. acknowledge support from the CFM foundation and the âInvestisements dâavenirâ (ANR-15-IDEX-02) programs of the French National Research Agency. K.B. and J.D acknowledge the European Unionâs Horizon 2020 research and innovation programme under the Marie SkÅodowska-Curie grant agreement No 754303.
Author information
Authors and Affiliations
Contributions
S.L., J.P.M., S.F. and N.R. designed the experiment. S.L. fabricated the device. S.L. performed the experiment and analysed the data with help from S.F., N.R. and I.S., while S.F. and I.S. provided the theoretical support. S.L., J.P.M., K.B., R.D., J.D., F.F., V.M., L.P., O.B., C.N., W.H.G., S.F., I.S. and N.R. participated in setting up the experimental platform, and took part in writing the paper.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review information Nature Communications thanks Gary Steele and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Léger, S., Puertas-MartÃnez, J., Bharadwaj, K. et al. Observation of quantum many-body effects due to zero point fluctuations in superconducting circuits. Nat Commun 10, 5259 (2019). https://doi.org/10.1038/s41467-019-13199-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-019-13199-x
This article is cited by
-
Emergent quantum phase transition of a Josephson junction coupled to a high-impedance multimode resonator
Nature Communications (2024)
-
Soliton confinement in a quantum circuit
Nature Communications (2023)
-
Bolometric detection of Josephson inductance in a highly resistive environment
Nature Communications (2023)
-
Evidence of dual Shapiro steps in a Josephson junction array
Nature Physics (2023)
-
Kerr reversal in Josephson meta-material and traveling wave parametric amplification
Nature Communications (2022)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.