Abstract
The study of lightâmatter interaction has led to important advances in quantum optics and enabled numerous technologies. Over recent decades, progress has been made in increasing the strength of this interaction at the single-photon level. More recently, a major achievement has been the demonstration of the so-called strong coupling regime1,2, a key advancement enabling progress in quantum information science. Here, we demonstrate lightâmatter interaction over an order of magnitude stronger than previously reported, reaching the nonperturbative regime of ultrastrong coupling (USC). We achieve this using a superconducting artificial atom tunably coupled to the electromagnetic continuum of a one-dimensional waveguide. For the largest coupling, the spontaneous emission rate of the atom exceeds its transition frequency. In this USC regime, the description of atom and light as distinct entities breaks down, and a new description in terms of hybrid states is required3,4. Beyond lightâmatter interaction itself, the tunability of our system makes it a promising tool to study a number of important physical systems, such as the well-known spin-boson5 and Kondo models6.
Similar content being viewed by others
Main
Light propagating in a one-dimensional (1D) waveguide is described by a 1D electromagnetic field with a continuous spectrum of frequencies. The strong coupling regime7 between an atom and such an electromagnetic continuum is defined as the regime in which the atom emits radiation predominantly into the waveguide with a rate ÎG that significantly exceeds the decoherence rate of the atom as well as emission into any other channel. In this regime, the atomic transition frequency Î far exceeds the emission rate ÎG ⪠Î. Achieving strong coupling to a continuum is a recent achievement in quantum optics8. Strong atomâwaveguide coupling has numerous applications, such as the development of quantum networks9 for quantum communication10 and quantum simulation11. This technology, first demonstrated with superconducting qubits in open transmission lines8,10,12,13, has also been implemented with both neutral atoms14 and quantum dots15 in photonic crystal waveguides. The distinctive signature of strong coupling is a decrease below 50% of the amplitude of transmitted signals due to coherent atomic scattering of photons.
A distinct regime of lightâmatter interaction is reached when ÎG becomes comparable to the atomic transition frequency ÎG/Î â¼ 0.1, the ultrastrong coupling (USC) regime. Most studies involving atom-field interactions are in the regime ÎG ⪠Î, where the common rotating-wave approximation (RWA) applies. In the USC regime, the RWA breaks down but perturbative treatments still allow an effective atom-field description when ÎG/Î â¼ 0.1 (refs 16,17). A novel, unexplored regime of lightâmatter interaction is the nonperturbative USC regime, where ÎG approaches or exceeds the atomic transition frequency ÎG/Î â¼ 1 and perturbation theory breaks down. This is a general definition also applicable to the case of discrete modes in cavity quantum electrodynamics systems18. We note that the nonperturbative USC regime has also been referred to in the literature as the deep strong coupling regime19. In the nonperturbative USC regime, the atomâphoton system is described by photons dressing the atom even in the ground state3,4,18. In this regime, the Markovian approximation also breaks down because the broad qubit linewidth ÎG implies that the spectral density of the environment seen by the atom is not independent of frequency. The presence of a continuum of modes ultrastrongly coupled to an atom has the additional effect of renormalizing the atomic frequency from the bare value Î0, which is a generalization of the well-known Lamb shift to arbitrary coupling strengths. These renormalization effects are also central to the well-known spin-boson model5, which has been used to describe, for example, open quantum systems20, quantum stochastic resonance21 and phase transitions in Josephson junctions22. Reaching the nonperturbative USC regime allows the exploration of the ultimate limits in lightâmatter interaction strength and relativistic quantum information phenomena23. In addition, ultrastrong couplings may have technological applications, such as single-photon nonlinearities24 and broadband single-photon sources3.
Superconducting qubits are artificial atoms with transitions in the microwave range of frequencies. Recently, flux-type superconducting qubits have been put forward as candidates to reach the nonperturbative USC regime3,25, having demonstrated large galvanic couplings to resonators17 and a large anharmonicity that allows them to remain an effective two-level system when ÎG â¼ Î. This is in contrast to other more weakly anharmonic qubits whose transitions would overlap for large enough ÎG.
Here, we demonstrate nonperturbative ultrastrong coupling of a superconducting flux qubit26 coupled to an open 1D transmission line via a shared Josephson junction. As predicted3,25, we observe that ÎG increases with the inverse of the coupling junction size. For devices with a small-enough coupling junction we measure ÎG â¼ Î, indicating that we reach the nonperturbative USC regime. Our flux qubit has four Josephson junctions. Two reference junctions are designed with the same area, while the areas of the other two junctions are scaled by the factors α â¼ 0.6 and β > 1 with respect to the area of the reference junctions27. The flux qubit is galvanically attached to the centre line of a 1D coplanar waveguide transmission line (Fig. 1a). To achieve ultrastrong couplings, we place the β-junction in parallel to the other three (Fig. 1b). The coupling to the line is then mainly determined3,25 by the matrix element between ground |0ã and excited |1ã qubit states of the superconducting phase operator across the β-junction , which is the dominant contribution to the coupling for β < 4. Further, we make the coupling tunable by turning the β-junction into a superconducting quantum interference device (SQUID) threaded by a flux Ïβ, as shown in Fig. 1b (Methods).
The experiments are performed by applying a probe field with a variable frequency and recording the transmitted field amplitude and phase on a vector network analyzer. For emission rates Î1/Π⪠1, where Î1 is the total emission rate and, in the presence of thermal excitations, the transmitted coherent scattering amplitude at low driving power is given by8,28:
Here Î2 â¡ ÎÏ + (Î1/2)(1 + 2nth) is the total decoherence rate, ÎÏ is the pure dephasing rate, Î´Ï = Ï â Î is the detuning of the probe field, and nth is the thermal photon occupation number at the qubit frequency (Supplementary Information). The maximum reflection amplitude is r0 = Î1/[2Î2(1 + 2nth)]. As in other experiments on superconducting quantum circuits8,10, relaxation into channels other than the waveguide is negligible. Therefore, we assume Î1 = ÎG. We note that equation (1) applies in the RWA. However, it has recently been shown4 that the scattering lineshapes are approximately Lorentzian in the USC regime up to Î1/Î â¼ 1 if we consider Î and Î1 to be renormalized parameters. This can be shown using a polaron transformation, allowing us to interpret the scattering centre as an atom dressed by a cloud of photons.
We first show measurements on a device with a fixed coupling junction with β â 3.5 (Fig. 1b). The transmission spectrum as a function of applied magnetic field (Fig. 2a) shows a maximum extinction at the symmetry point of 95%, indicating strong coupling. By fitting equation (1) (dashed line), we infer Î1/2Ï = 88 ± 11âMHz (see Methods), Î/2Ï = 3.996 ± 0.001âGHz, giving Î1/Î = 0.02, which is not in the USC regime. Flux qubit spectra in transmission lines similar to this one have previously been reported8,29.
To enhance the coupling strength, we designed a second device where the size of the β-junction was decreased to β â 1.8. The resulting qubit spectrum in Fig. 2b shows striking differences compared to the previous device with β â 3.5. The qubit linewidth at the symmetry point is very large, comparable to the total measurement bandwidth of 3â11âGHz. The deviations from a Lorentzian lineshape are due to bandwidth limitations of our set-up, still allowing us to infer a full-width at half-maximum of 2Î2/2Ï â 10.90 ± 0.44âGHz (see Methods). The extracted qubit emission rate Î1/2Ï â 9.24 ± 0.52âGHz exceeds the qubit splitting Î/2Ï = 7.68 ± 0.08âGHz, giving Î1/Î = 1.20 ± 0.07, a clear indication that this device reaches the nonperturbative USC regime.
Having observed two devices with Î1 ⪠Πand Î1 > Î, we now explore the intermediate region using a device with tunable coupling (Fig. 1b) designed with a tunable range of β â¼ 1.6â3.6. In Fig. 3aâc, spectroscopy of the tunable coupling device is shown at three different values of Ïβ. Using scanning-electron microscope (SEM) images of the measured device, we identify Fig. 3aâc as effectively having, respectively, β(a) â 3.6, β(b) â 2.3, β(c) â 1.6. Figure 3a corresponds to the highest effective β-junction size, therefore the lowest coupling strength. A flux qubit spectrum can be identified with Î/2Ï = 5.20 ± 0.02âGHz and 2Î2/2Ï â 2.40 ± 0.07âGHz. The maximum extinction at the symmetry point is over 95%. The quality of the signal below 4âGHz degrades due to the measurement taking place outside the optimal bandwidth of our amplifier and circulators (4â8âGHz, Supplementary Information). In Fig. 3b, the qubit gap decreases to Î/2Ï â 2.90 ± 0.05âGHz, as expected for a smaller β-junction. The width 2Î2/2Ï = 5.90 ± 0.22âGHz is clearly enhanced, with the extinction decreasing to 30%. In Fig. 3c, the qubit spectrum is barely discernible. The extinction is only 10%, with a response that appears featureless in our frequency range. Figure 3d, e shows the extracted values of r0 and Î2 using equation (1). The value of 2Î2/2Ï â 13 ± 3âGHz from Fig. 3c is an inferred bound due to the difficulty in fitting the transmission at this value of flux.
To understand the spectrum of the tunable coupling device and extract the corresponding emission rates Î1, we need to take into account finite temperature effects. We can set an upper bound on nth, which is (Methods). Figure 3f shows that the values of nmax for β > 2 are consistent with a unique maximum effective temperature of Teff = 90 mK, comparable to other superconducting qubit experiments. Using 0 < nth < nmax, we then put bounds on Î1: . Using these bounds, we plot Î1/Î in Fig. 4a. The plot clearly shows that we can tune the device from the regime of strong coupling all the way into the nonperturbative USC regime. The curve in Fig. 4a corresponds to the theoretical value of the normalized coupling strength (Supplementary Information)
with RQ = h/(2e)2 = 6.5 kΩ the resistance quantum and Z0 the characteristic impedance of the line. The matrix element values of the phase operator across the coupling junction β, |Ïβ|2, are calculated using the methods of ref. 3. The observed values of Î1/Î agree very well with the calculated values based on our circuit3 for an impedance close to the nominal 50 Ω. Above Î1/Î â Ï/2, equation (2) becomes a lower bound (Supplementary Information). This is consistent with data in the range β < 2 lying above equation (2). Including renormalization effects5 in equation (2) might further improve the agreement with the measurements for β < 2.
Our system allows us to explore the spin-boson (SB) model in an ohmic bath. According to the SB model, the high-frequency modes of the transmission line renormalize the bare qubit splitting Î0 to4,5
αSB is the SB normalized coupling strength that is related to the spectral density of the environment J(Ï). For an ohmic system such as our transmission line, αSB = J(Ï)/ÏÏ. ÏC â« Î0 is the cutoff frequency of the environment and p is a constant of order 1. Up to αSB â¼ 0.5, we identify αSB = Î1/ÏÎ. Above αSB â 0.5 (or Î1/Î â Ï/2) this relation becomes a lower bound for αSB (Supplementary Information). In Fig. 4b we plot the experimental qubit splittings Î (circles). Using qubit junction dimensions extracted from SEM images of the device, we diagonalize the qubit Hamiltonian at each flux Ïβ (triangles) to give the bare qubit gaps Î0. We then renormalize the calculated Î0 using equation (3) and a value of p = exp(1 + γ) â 4.8, which is derived using an exponential cutoff model4,5. γ is the Euler constant. We find the best fit to the measured Î using a cutoff of ÏC/2Ï = 50âGHz, which is consistent with characteristic system frequencies such as the plasma frequency of the qubit junctions and the superconducting gap. The agreement between the observed qubit splittings Î and our estimates of the renormalized gaps is clear3,4,5.
As a prelude to future work, we can place our results in the context of the SB model. The SB model defines three dynamical regimes for the qubit: underdamped (αSB < 0.5), overdamped (1 > αSB > 0.5) and localized (αSB > 1). The connection between Î1/Î and αSB allows us to draw the boundaries between these regimes in Fig. 4a. We see that our tunable device enters well into the overdamped regime, and very possibly into the localized regime for β < 2. More detailed measurements of the dynamics of the device in these regimes could further confirm the predictions of the SB model. Suggestively, the strong reduction of the qubit response seen in Fig. 3c (leftmost data points in Fig. 4a) with a flat response as a function of frequency is consistent with simulations of classical double-well dynamics in the overdamped regime (P. Forn-DÃaz, manuscript in preparation).
We have presented measurements of superconducting flux qubits in 1D open transmission lines in regimes of interaction starting at strong coupling and ranging deeply into the ultrastrong coupling regime. In particular, we observed qubits with emission rates exceeding their own frequency, a clear indication of nonperturbative ultrastrong coupling. These results are very relevant for the study of open systems in the USC regime, opening the door to the development of a new generation of quantum electronics with ultrahigh bandwidth for quantum and nonlinear optics applications. The tunability of our system also makes it well-suited to the simulation of other quantum systems. In particular, we showed that the device can span the various transition regions of the SB model. With further development of our quantum circuit, the structure of the photon dressing cloud could also be directly detected, allowing the study of the physics of the Kondo model6 in a well-controlled setting. The ultrastrong coupling regime has other interesting intrinsic properties on its own, such as the entangled nature of the ground state.
Note added in proof: After acceptance of our paper, a related manuscript was published30 showing similar results to this work using a resonator instead of a transmission line.
Methods
Device details and fabrication.
We made the device with tunable coupling by replacing the β-junction with a SQUID threaded by a flux Ïβ. The tunable coupling device then consists of two loops, the main loop that changes primarily the qubit magnetic energy through the flux Ïε and the β-loop that changes the effective coupling to the transmission line through Ïβ. Changing β also modifies the minimum qubit splitting Î. To minimize this effect, we make the SQUID junctions asymmetric, which lowers the sensitivity of Î to Ïβ. Similar tunable coupling architectures were already suggested in ref. 31. In the experiment, we sweep the global magnetic field, therefore simultaneously changing Ïε and Ïβ. The qubit spectrum shows minima near Ïε â Ï0(1/2 + n), with Ï0 = h/2e the quantum of flux, n being an integer (Supplementary Information). Here, different n will correspond to different Ïβ, leading to different coupling strengths. The loop areas Aε/Aβ are designed to have a large, incommensurate ratio, allowing the exploration of many different values of β.
The fabrication methods used are based on those of ref. 27. The fabrication of devices starts by patterning the transmission line using optical lithography followed by an evaporation of 200ânm of aluminium. A gap in the transmission line is left to place the qubit in a second lithography stage. We pattern the qubit using an electron beam writer. Prior to the second aluminium evaporation an Ar milling step is applied to remove the native oxide on the first aluminium layer, guaranteeing optimal conduction between the two aluminium layers. The qubit is evaporated using double-angle shadow mask evaporation, resulting in a total thickness of 105ânm. After the first shadow evaporation step, we oxidize the film with dynamical flow at â¼0.01âmbar for 7âmin, yielding critical current densities of â¼12âμAâμmâ2. The chip is then diced and the transmission line is wire-bonded to a printed circuit board connecting to the rest of the circuitry in our cryostat.
The transmission line consists of a 6.5âmm long on-chip coplanar waveguide with a centre line and gaps 8âμm and 4âμm wide, respectively, resulting in a 50 Ω characteristic impedance. Numerical simulations are run to verify the impedance of the circuit. We use a squared webbed ground to reduce superconducting vortex motion on the ground plane.
Bounds on qubit emission rate.
The dependence of r0 and Î2 on nth shown below equation (1) does not allow the independent extraction of all parameters, Î1, ÎÏ, nth at each value of β. However, we can set bounds on nth. The lower bound case assumes no thermal excitations, therefore nth = 0. If we instead set ÎÏ = Î2 (1 â r0(1 + 2nth)2) ⥠0, we identify an upper bound on the photon occupation number . In Fig. 3f, the values of nmax were extracted assuming ÎÏ = 0. If we were to assume ÎÏ/2Ï = 17âMHz as the nonthermal dephasing rate, extracted from the narrower linewidth of the device in Fig. 2a assuming nth = 0, the resulting nth would not differ significantly from nmax. Now, bounds on Î1 = 2Î2r0(1 + 2nth) can be set as Î1(nth = 0) and Î1(nth = nmax), giving . The lower bound, nth = 0, is close to the calculated value of nth at the cryostat temperature of 10âmK for all qubit frequencies.
Spectroscopic analysis.
In all data shown, we use equation (1) to simultaneously fit the real and imaginary parts of the transmission. Supplementary Section 3 shows the full set of fitted resonances used in Figs 3 and 4 of the main text. Note that the baseline is fixed to a normalized value of 1 and is not adjusted. The baseline value is itself determined by measuring the transmitted background when the qubit is flux-tuned away from the frequency band of interest.
Data availability.
The data that support the plots within this paper and other findings of this study are available from the corresponding authors upon request.
References
Haroche, S. & Raimond, J.-M. Exploring the Quantum: Atoms, Cavities, and Photons (Oxford Univ. Press, 2006).
Wallraff, A. et al. Circuit quantum electrodynamics: coherent coupling of a single photon to a Cooper pair box. Nature 431, 162â167 (2004).
Peropadre, B., Zueco, D., Porras, D. & GarcÃa-Ripoll, J. J. Nonequilibrium and nonperturbative dynamics of ultrastrong coupling in open lines. Phys. Rev. Lett. 111, 243602 (2013).
DÃaz-Camacho, G., Bermudez, A. & GarcÃa-Ripoll, J. J. Dynamical polaron ansatz: a theoretical tool for the ultra-strong coupling regime of circuit QED. Phys. Rev. A 93, 043843 (2016).
Leggett, A. J. et al. Dynamics of the dissipative two-state system. Rev. Mod. Phys. 59, 1â85 (1987).
Le Hur, K. Kondo resonance of a microwave photon. Phys. Rev. B 85, 140506 (2012).
Shen, J.-T. & Fan, S. Coherent single photon transport in a one-dimensional waveguide coupled with superconducting quantum bits. Phys. Rev. Lett. 95, 213001 (2005).
Astafiev, O. et al. Resonance fluorescence of a single artificial atom. Science 327, 840â843 (2010).
Kimble, H. J. The quantum internet. Nature 453, 1023â1030 (2008).
Hoi, I.-C. et al. Demonstration of a single-photon router in the microwave regime. Phys. Rev. Lett. 107, 073601 (2011).
GarcÃa-Ãlvarez, L. et al. Fermionâfermion scattering in quantum field theory with superconducting circuits. Phys. Rev. Lett. 114, 070502 (2015).
Hoi, I.-C. et al. Giant cross-Kerr effect for propagating microwaves induced by an artificial atom. Phys. Rev. Lett. 111, 053601 (2013).
van Loo, A. F. et al. Photon-mediated interactions between distant artificial atoms. Science 342, 1494â1496 (2013).
Goban, A. et al. Superradiance for atoms trapped along a photonic crystal waveguide. Phys. Rev. Lett. 115, 063601 (2015).
Arcari, M. et al. Near-unity coupling efficiency of a quantum emitter to a photonic crystal waveguide. Phys. Rev. Lett. 113, 093603 (2014).
Niemczyk, T. et al. Circuit quantum electrodynamics in the ultrastrong-coupling regime. Nat. Phys. 6, 772â776 (2010).
Forn-DÃaz, P. et al. Observation of the BlochâSiegert shift in a qubit-oscillator system in the ultrastrong coupling regime. Phys. Rev. Lett. 105, 237001 (2010).
Ciuti, C., Bastard, G. & Carusotto, I. Quantum vacuum properties of the intersubband cavity polariton field. Phys. Rev. B 72, 115303 (2005).
Casanova, J. et al. Deep strong coupling regime of the JaynesâCummings model. Phys. Rev. Lett. 105, 263603 (2010).
Grifoni, M., Paladino, E. & Weiss, U. Dissipation, decoherence and preparation effects in the spin-boson system. Eur. Phys. J. B 10, 719â729 (1999).
Gammaitoni, L., Hänggi, P., Jung, P. & Marchesoni, F. Stochastic resonance. Rev. Mod. Phys. 70, 223â287 (1998).
Penttilä, J. S., Parts, U., Hakonen, P. J., Paalanen, M. A. & Sonin, E. B. Superconductorâinsulator transition in a single Josephson junction. Phys. Rev. Lett. 82, 1004â1007 (1999).
SabÃn, C., Peropadre, B., Rey, M. D. & MartÃn-MartÃnez, E. Extracting past-future vacuum correlations using circuit QED. Phys. Rev. Lett. 109, 033602 (2012).
Sánchez-Burillo, E., Zueco, D., GarcÃa-Ripoll, J. J. & MartÃn-Moreno, L. Scattering in the ultrastrong regime: nonlinear optics with one photon. Phys. Rev. Lett. 113, 263604 (2014).
Bourassa, J. et al. Ultrastrong coupling regime of cavity QED with phase-biased flux qubits. Phys. Rev. A 80, 032109 (2009).
Mooij, J. E. et al. Josephson persistent-current qubit. Science 285, 1036â1039 (1999).
Orgiazzi, J. L. et al. Flux qubits in a planar circuit quantum electrodynamics architecture: quantum control and decoherence. Phys. Rev. B 93, 104518 (2016).
Peropadre, B. et al. Scattering of coherent states on a single artificial atom. New J. Phys. 15, 035009 (2013).
Haeberlein, M. et al. Spin-boson model with an engineered reservoir in circuit quantum electrodynamics. Preprint at https://arXiv.org/abs/1506.09114 (2015).
Yoshihara, F. et al. Superconducting qubitâoscillator circuit beyond the ultrastrong-coupling regime. Nat. Phys. http://dx.doi.org/10.1038/nphys3906 (2016).
Peropadre, B., Forn-DÃaz, P., Solano, E. & GarcÃa-Ripoll, J. J. Switchable ultrastrong coupling. Phys. Rev. Lett. 105, 023601 (2010).
Acknowledgements
We acknowledge financial support from NSERC of Canada, the Canadian Foundation for Innovation, the Ontario Ministry of Research and Innovation, Industry Canada, Canadian Microelectronics Corporation, EU FP7 FET-Open project PROMISCE, Spanish Mineco Project FIS2012-33022 and CAM Network QUITEMAD+. B.P. acknowledges the Air Force of Scientific Research for support under award FA9550-12-1-0046. The University of Waterlooâs Quantum NanoFab was used for this work. We thank A. J. Leggett and A. Garg for fruitful discussions, and M. Otto, S. Chang, A. M. Vadiraj and C. Deng for help with device fabrication and with the measurement set-ups.
Author information
Authors and Affiliations
Contributions
P.F.-D., C.M.W. and A.L. designed the experiment. P.F.-D. designed the devices and fabricated them. P.F.-D., C.M.W. and A.L. conducted the experiments. J.-L.O. provided input to device design and fabrication. M.A.Y. and R.B. assisted in numerical modelling of devices. J.J.G.-R. and B.P. provided theoretical support to interpret the measurements. P.F.-D., C.M.W. and A.L. performed the data analysis and wrote the manuscript with feedback from all authors. C.M.W. and A.L. supervised the project.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary information
Supplementary information (PDF 894 kb)
Rights and permissions
About this article
Cite this article
Forn-DÃaz, P., GarcÃa-Ripoll, J., Peropadre, B. et al. Ultrastrong coupling of a single artificial atom to an electromagnetic continuum in the nonperturbative regime. Nature Phys 13, 39â43 (2017). https://doi.org/10.1038/nphys3905
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/nphys3905
This article is cited by
-
Emergent quantum phase transition of a Josephson junction coupled to a high-impedance multimode resonator
Nature Communications (2024)
-
Probing the symmetry breaking of a lightâmatter system by an ancillary qubit
Nature Communications (2023)
-
Enhancing the direct charging performance of an open quantum battery by adjusting its velocity
Scientific Reports (2023)
-
Extremely large Lamb shift in a deep-strongly coupled circuit QED system with a multimode resonator
Scientific Reports (2023)
-
Noisy intermediate-scale quantum computers
Frontiers of Physics (2023)