<figure><img src="image_1.png" /></figure>

HAL
open science

# Energy-participation quantization of Josephson circuits

Zlatko Minev, Zaki Leghtas, Shantanu Mundhada, Lysander Christakis, Ioan Pop, Michel Devoret

► To cite this version:

Zlatko Minev, Zaki Leghtas, Shantanu Mundhada, Lysander Christakis, Ioan Pop, et al., Energy-participation quantization of Josephson circuits. npj Quantum Information, 2021, 7 (1), pp.131. 10.1038/s41534-021-00461-8. [hal-03467417](https://hal-INRAE-03467417)

HAL Id: hal-03467417

https://minesparis-psl.hal.science/hal-03467417v1

Submitted on 10 Sep 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
ARTICLE OPEN

<figure><img src="image_3.png" /></figure>

# Energy-participation quantization of Josephson circuits

Zlatko K. Minev<sup>1,4,✉</sup>, Zaki Leghtas<sup>1,2</sup>, Shantanu O. Mundhada<sup>1,5</sup>, Lysander Christakis<sup>1,6</sup>, Ioan M. Pop<sup>1,3</sup> and Michel H. Devoret<sup>1</sup>

Superconducting microwave circuits incorporating nonlinear devices, such as Josephson junctions, are a leading platform for emerging quantum technologies. Increasing circuit complexity further requires efficient methods for the calculation and optimization of the spectrum, nonlinear interactions, and dissipation in multi-mode distributed quantum circuits. Here we present a method based on the energy-participation ratio (EPR) of a dissipative or nonlinear element in an electromagnetic mode. The EPR, a number between zero and one, quantifies how much of the mode energy is stored in each element. The EPRs obey universal constraints and are calculated from one electromagnetic-eigenmode simulation. They lead directly to the system quantum Hamiltonian and dissipative parameters. The method provides an intuitive and simple-to-use tool to quantize multi-junction circuits. We experimentally tested this method on a variety of Josephson circuits and demonstrated agreement within several percents for nonlinear couplings and modal Hamiltonian parameters, spanning five orders of magnitude in energy, across a dozen samples.

npj Quantum Information (2021)7:131; https://doi.org/10.1038/s41534-021-00461-8

## INTRODUCTION

Quantum information processing based on the control of microwave electromagnetic fields in Josephson circuits is a promising platform for both fundamental physics experiments and emerging quantum technologies<fnref n="1" />⁻³. Key to the success of this platform is the ability to quantitatively model the distributed quantized electromagnetic modes of the system, their nonlinear interactions, and their dissipation (see Fig. 1). This challenge is the subject of intensifying interest⁴⁻²⁶, as experimental architectures²⁷⁻³⁶ and nonlinear devices³⁷⁻⁴⁴ scale in both complexity and diversity.

In this paper, we introduce a circuit quantization method based on the concept of the energy-participation ratio (EPR). We reduce the quantization problem to answering the simple question: what fraction of the energy of mode *m* is stored in element *j*? This leads to a constrained number between zero and one, the EPR, denoted *p*<sub>mj</sub><sup>45</sup>. This ratio is the key quantity that bridges classical and quantum circuit analysis; we show that it plays the primary role in the construction of the system many-body Hamiltonian. Furthermore, dissipation in the system is treated on equal footing by calculating the EPR *p*<sub>ml</sub> of lossy element *l* in mode *m*.

The EPR method deviates from previous black-box quantization work<sup>4,6,7</sup>, which uses the impedance-response matrix, denoted $Z_{jj'}(\omega)$, where *j* and *j'* are index ports associated with nonlinear elements. For all pairs of ports, the complex function $Z_{jj'}(\omega)$ is calculated from a finite-element (FE) driven simulation in the vicinity of the eigenfrequency of every mode. Our method replaces these steps with a more economical FE eigenmode simulation, from which one extracts the energy participations $p_{mj}$ and $p_{mj}$ needed to fully characterize both the dissipative and Hamiltonian properties of the circuit.

To test the method, we compared EPR calculations of circuit parameters to experimentally measured ones for 8 superconducting devices designed with the EPR method, comprising a total of 15 qubits, 8 readout and storage resonator modes, and 1 waveguide system. The results demonstrate agreement for

Hamiltonian parameters spanning over five orders of magnitude in energy. Resonance frequencies were calculated to 1% accuracy, large nonlinear interactions, such as anharmonicities and cross-Kerr frequencies, to 5%, and small, nonlinear interactions to 10%. This level of accuracy is sufficient for most current quantum information experiments.

## RESULTS AND DISCUSSION

### To quantize a simple circuit: qubit coupled to a cavity

In this section, we introduce the EPR method of quantum circuit design on a modest, yet informative, example: a transmon qubit coupled to a cavity mode (see Fig. 2). The transmon<sup>46</sup> consists of a Josephson junction shunted by a capacitance. It is embedded in the cavity, which we will consider as a black-box distributed structure. The Hamiltonian of this system $\hat{H}_{\text{full}}$ can be conceptually separated into two contributions (see Supplementary Section A2),

$$ \hat{H}_{\text{full}} = \hat{H}_{\text{lin}} + \hat{H}_{\text{nl}}, \quad (1) $$

where $\hat{H}_{\text{lin}}$ consists of all terms associated with the linear response of the junction and the resonator structure, and $\hat{H}_{\text{nl}}$ consists of terms associated with the nonlinear response of the junction. Restricting our attention to the cavity and qubit modes of the otherwise black-box structure, the analytical form of the Hamiltonian follows from standard circuit quantization<sup>47,48</sup> (see Supplementary Information):

$$ \hat{H}_{\text{lin}} = \hbar\omega_c \hat{a}_c^\dagger \hat{a}_c + \hbar\omega_q \hat{a}_q^\dagger \hat{a}_q, \quad (2) $$

$$ \hat{H}_{\text{nl}} = -E_j[\cos(\hat{\phi}_j) + \hat{\phi}_j^2/2], \quad (3) $$

$$ \hat{\phi}_j = \varphi_q (\hat{a}_q + \hat{a}_q^\dagger) + \varphi_c (\hat{a}_c + \hat{a}_c^\dagger), \quad (4) $$

where $\omega_c$ and $\omega_q$ are the angular frequencies of the cavity and qubit eigenmodes defined associated with $\hat{H}_{\text{lin}}$, respectively, and where $\hat{a}_c$ and $\hat{a}_q$ are their annihilation operators, respectively. The Josephson energy $E_j$ can be computed from the

<footnote n="1">¹Department of Applied Physics, Yale University, New Haven, CT, USA. ²Centre Automatique et Systèmes, Mines-ParisTech, PSL Research University, Paris, France. ³IQMT, Karlsruhe Institute of Technology, Eggenstein-Leopoldshafen, Germany. ⁴Present address: IBM T.J. Watson Research Center, Yorktown Heights, NY, USA. ⁵Present address: Quantum Circuit Incorporated (QCI), New Haven, CT, USA. ⁶Present address: Department of Physics, Princeton University, Princeton, NJ, USA. ✉email: zlatko.minev@aya.yale.edu</footnote>
<figure><img src="image_4.png" /><figcaption>**Fig. 1** Conceptual overview. **a** Illustration of the physical model of an example quantum device, which comprises two three-dimensional (3D) cavities (gray enclosures), each housing several qubit chips (green boxes). A close-up view of one of the chips is depicted in the inset. The dotted box in the center of the chip schematically outlines a nonlinear inductive sub-circuit, referred to as a Josephson dipole. **b** Results of a finite-element eigenmode analysis (FE<sub>e</sub>) of the Josephson circuit linearized about its equilibrium. The mth mode eigenfrequency and electric and magnetic fields are ω<sub>m</sub>, ⃗E<sub>m</sub>(r), and ⃗H<sub>m</sub>(r), respectively, where r denotes spatial position. Center inset: |⃗E<sub>m</sub>| profile (red: high; blue: low) for the fundamental mode of one of the 3D cavities. Additional FE driven simulations (FE<sub>d</sub>) are unnecessary; i.e., the impedance matrix Z<sub>jj'</sub>(ω) is not calculated. **c** The Hamiltonian ⃗H<sub>full</sub>, which includes nonlinear interactions to arbitrary order (see "Results"), is computed directly from the eigenanalysis via the EPRs p<sub>mj</sub> and EPR signs s<sub>mj</sub> = ±1 of the junctions, j. Dissipative contributions due to a lossy element l are similarly computed from the loss EPRs p<sub>ml</sub> for linear dissipation, and the EPR signs s<sub>ml</sub> are unnecessary. Direct extraction of Hamiltonian and dissipative parameters from eigensolutions is unique to the EPR method. The geometry of the classical model is modified in an iterative search for the desired dissipative and Hamiltonian parameters (left-pointing arrow).</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Fig. 2** Quantizing a simple circuit. **a** Illustration of a 3D cavity enclosing a transmon qubit chip. The cross symbol marks the location of a Josephson junction. Vertical blue arrows depict the electric field ⃗E<sub>m</sub>(r) of the fundamental cavity mode, TE<sub>101</sub>. **b** Equivalent two-mode lumped-element representation of the distributed circuit. Operators ˆa<sub>q</sub> and ˆa<sub>c</sub> denote the qubit and cavity mode operators, respectively.</figcaption></figure>

Ambegaokar-Baratoff formula adapted to the measured room-
temperature resistance of the junction⁴⁹. The junction-reduced
generalized flux ˆφ<sub>J</sub> corresponds to the classical variable
φ<sub>J</sub>(t) := ∫<sup>t</sup><sub>−∞</sub> v<sub>J</sub>(t) dt/φ<sub>0</sub>, where v<sub>J</sub>(t) is the instantaneous voltage
across the junction<sup>47,48</sup>, and φ<sub>0</sub> := ℏ/2e is the reduced flux
quantum. The junction flux operator [Eq. (4)] is a linear, real-
valued, and non-negative combination of the mode operators (see
Supplementary Section A5), and in its expression, φ<sub>c</sub> and φ<sub>q</sub> are
the quantum zero-point fluctuations of junction flux in the cavity
and qubit mode, respectively. It is worth stating that the linear
coupling between the cavity and qubit, commonly denoted g<sup>46</sup>, is
fully factored in our analysis and is implicitly handled in the
extraction of the operators from the electromagnetic simulation.

Our principal aim is to determine the unknown quantities: ω_c, ω_q,
φ_q, and φ_c. As we will show, we extract and compute these
quantities from an eigenanalysis of the classical distributed circuit
corresponding to H<sub>lin</sub>. This includes the qubit-cavity layout,
materials, electromagnetic boundary conditions, and a model of
the junction as a lumped-element, linear inductor. The eigensolver
returns the requested set of eigenmodes and their frequencies,
quality factors, and field solutions. By running the eigensolver in
the frequency range of interest, we obtain the hybridized cavity
and qubit modes, whose eigenfrequencies ω_c and ω_q fully
determine Ĥ<sub>lin</sub> (see Supplementary Section C for FE methodology).

To determine H<sub>nl</sub>, we need the quantum zero-point fluctuations φ<sub>q</sub> and φ<sub>c</sub>, which are calculated from the participation of the junction in the eigenfield solutions. The participation p<sub>m</sub> of the

junction in mode $m \in \{c, q\}$ is defined to be the fraction of
inductive energy stored in the junction relative to the total
inductive energy stored in the entire circuit,

$$
p_m := \frac{\text{Inductive energy stored in the junction}}{\text{Total inductive energy stored in mode } m'} \quad (5)
$$

evaluated when only mode *m* is excited. Thus, *p*<sub>*m*</sub> can be
computed from the electric ⃗*E*<sub>*m*</sub>(*r*) and magnetic ⃗*H*<sub>*m*</sub>(*r*)
eigenfields as detailed in Supplementary Section C2; ⃗*r* denotes
spatial position. In the quantum setting, Eq. (5) links *p*<sub>*m*</sub>, φ<sub>*J*</sub>, and
the state of the circuit,

$$
p_m = \frac{\langle \psi_m | \frac{1}{2} E_J \dot{\phi}_J^2 | \psi_m \rangle}{\langle \psi_m | \frac{1}{2} \hat{H}_{\text{lin}} | \psi_m \rangle}, \quad (6)
$$

where |ψ<sub>m</sub>⟩ denotes a coherent state or a Fock excitation of mode m. Note that normal ordering must be used in Eq. (6); this correct treatment of vacuum fluctuations is detailed in Supplementary Section A6. Simplifying Eq. (6), one expresses the variance of the quantum zero-point fluctuations φ<sub>c</sub> and φ<sub>q</sub> as a function of the classical energy participations p<sub>m</sub>,

$$
\phi_c^2 = p_c \frac{\hbar \omega_c}{2E_J} \quad \text{and} \quad \phi_q^2 = p_q \frac{\hbar \omega_q}{2E_J}, \qquad (7)
$$

which completely determines ⃗H<sub>nl</sub> and thus completes the
description of the system Hamiltonian ⃗H<sub>full</sub>. Here φ<sub>c</sub> and φ<sub>q</sub> can
be taken as positive numbers. As we will see in the next section, in
the presence of multiple junctions, this is not always true.

Designing experiments with the EPR requires one to further
extract from $\hat{H}_{\text{full}}$ the transition frequencies and nonlinear
couplings between modes. Depending on the case, this can be
done approximately or exactly using numerical or analytical
techniques²⁰. This task is easily achieved if $\hat{H}_{\text{nl}}$ is a perturbation to
$\hat{H}_{\text{lin}}$⁴⁶. In this limit, $\hat{H}_{\text{full}}$ for our qubit–cavity example can be
approximated by the effective, excitation-number-conserving
Hamiltonian, see Supplementary Section B1,

$$
\hat{H}_{\text{eff}} = (\omega_q - \Delta_q) \hat{n}_q + (\omega_c - \Delta_c) \hat{n}_c - \chi_{qc} \hat{n}_q \hat{n}_c \\
- \frac{1}{2} a_q \hat{n}_q (\hat{n}_q - \hat{1}) - \frac{1}{2} a_c \hat{n}_c (\hat{n}_c - \hat{1}), \qquad (8)
$$

where $\hat{n}_q = \hat{a}_q^\dagger \hat{a}_q$ and $\hat{n}_c = \hat{a}_c^\dagger \hat{a}_c$ denote the qubit and cavity excitation-number operators, respectively, $\Delta_q$ denotes the "Lamb shift" of the qubit frequency due to the dressing of this nonlinear
<figure><img src="image_6.png" /><figcaption>Fig. 3 Schematic representation of the Josephson circuit and its nonlinear elements. a A simple example of a Josephson dipole—a Josephson tunnel junction. b An example of a composite junction, comprising four Josephson junctions in a ring, frustrated by an external magnetic flux $\Phi_j^{\text{ext}}$ threading the loop. c Conceptual decomposition of a general Josephson dipole, denoted $j$. For convenience, its potential energy function $\mathcal{E}_j(\Phi_j; \Phi_j^{\text{ext}})$ can be Taylor expanded in a sum of nonlinear inductive contributions of increasing order $\Phi_j^p$, with relative amplitude $c_{jp}$, where $p$ denotes the index in the series. The energy function can be subjected to external bias parameters $\Phi_j^{\text{ext}}$, such as flux or voltages. d Schematic diagram of a general Josephson circuit conceptually resolved into a purely dissipative (left, $\kappa_m^l$), linear (middle, $\hat{H}_{\text{lin}}$), and nonlinear (right, $\hat{H}_{\text{nl}}$) constitutions.</figcaption></figure>

mode by quantum fluctuations of the fields, $a_q$ ($a_c$) is the qubit (cavity) anharmonicity, and $\chi_{qc}$ is the qubit–cavity dispersive shift (cross-Kerr coupling). The Hamiltonian parameters can be calculated directly from the EPR, see Supplementary Section B2,

$$a_q = \frac{1}{2} \chi_{qq} = p_q^2 \frac{\hbar \omega_q^2}{8 E_j}, \quad (9)$$

$$a_c = \frac{1}{2} \chi_{cc} = p_c^2 \frac{\hbar \omega_c^2}{8 E_c}, \quad (10)$$

$$\chi_{qc} = p_q p_c \frac{\hbar \omega_q \omega_c}{4 E_j}. \quad (11)$$

Experimentally, the qubit Lamb shift can be obtained as $\Delta_q = a_q - \chi_{qc}/2$. Since a single EPR $p_m$ determines the nonlinear interaction for each mode, the parameters $\chi_{qc}$ and $a_q$ are interdependent,

$$\chi_{qc} = \sqrt{\chi_{qq}\chi_{cc}} = 2\sqrt{a_q a_c}. \quad (12)$$

As shown in Supplementary Section A7, the EPRs $p_c$ and $p_q$ obey the constraints

$$0 \le p_q, p_c \le 1 \text{ and } p_q + p_c = 1. \quad (13)$$

These relations together with Eqs. (9) and (12) are useful to budget the dilution of the nonlinearity of the junction (see Supplementary Section B2) and to provide insight on the limits of accessible parameters (see “Methods”). Further, Eq. (13) is used to validate the convergence of the FE simulation.

## Quantizing the general Josephson system

The simple results obtained in the preceding section will now be generalized to arbitrary nonlinear devices enclosed in a black-box, distributed, electromagnetic structure. While such structures are frequently classified as planar<sup>27,50–52</sup> (two-dimensional), quasi-planar<sup>28,29,32,53</sup> (2.5-dimensional (2.5D)), or three-dimensional<sup>30,54–56</sup>

(3D), we will treat all classes on equal footing. The electromagnetic structure is assumed to be linear in the absence of the enclosed nonlinear devices. For simplicity of discussion, we can consider these devices to be inductive and lumped; distributed nonlinear devices, such as kinetic-inductance transmission lines<sup>57–60</sup>, can be thought of as a series of lumped ones.

The simplest nonlinear device comprises a single element, such as a Josephson tunnel junction [see Fig. 3a], an atomic-point contact<sup>61,62</sup>, a nanobridge<sup>38,63</sup>, a semiconducting nanowire<sup>40,41,64–67</sup>, or another hybrid structure<sup>68</sup>. A multi-element device, such as a SQUID<sup>69,70</sup>, a SNAIL<sup>71</sup> [see Fig. 3b], a superinductance<sup>72–74</sup>, or a junction array<sup>5,75–80</sup> refers to a subcircuit composed of purely inductive lumped elements. This subcircuit can also be subjected to external controls, such as voltage or flux biases.

The general nonlinear device that we now consider, referred to as a Josephson dipole, is any lumped, purely inductive, nonlinear subcircuit with two terminals. The key characteristic of the Josephson dipole is that it possesses a characteristic energy function, which encapsulates all details of its constitution. For example, the two-terminal nonlinear device known as the symmetric SQUID<sup>70</sup> is described by the energy function $\mathcal{E}_j(\Phi_j; \Phi_j^{\text{ext}}) = -E_j(\Phi_j^{\text{ext}}) \cos(\Phi_j/\phi_0)$, where $\Phi_j$ is the generalized flux across the device terminals<sup>47,48</sup>, $E_j$ is the effective Josephson energy, $\Phi_j^{\text{ext}}$ is the external flux bias, and the subscript $j$ denotes the $j$th Josephson dipole in the circuit. The flux $\Phi_j$ is defined as the deviation away from the value in equilibrium, as discussed below. To ease the notation, parameters such as $\Phi_j^{\text{ext}}$ will be implicit hereafter. Similarly to the example of the single-junction transmon, the energy of a Josephson dipole can be separated in two parts. One part $\mathcal{E}_j^{\text{lin}}$ accounts for the linear response of the dipole, while the other $\mathcal{E}_j^{\text{nl}}$ accounts for the nonlinear response,

$$\mathcal{E}_j(\Phi_j) = \mathcal{E}_j^{\text{lin}}(\Phi_j) + \mathcal{E}_j^{\text{nl}}(\Phi_j), \quad (14)$$

where

$$\mathcal{E}_j^{\text{lin}}(\Phi_j) := \frac{1}{2} E_j \left( \frac{\Phi_j}{\phi_0} \right)^2, \quad (15)$$

and where the constant $E_j$ sets the scale of the junction energy. This energy scale can be represented by the linear inductance $L_j := \phi_0/E_j$ presented by the Josephson dipole when submitted to a small excitation about its equilibrium.

**Frustrated equilibrium.** External biases can set up persistent currents in the circuit. These can alter the static (direct current) equilibrium of the Josephson system. For example, frustrating a superconducting ring with a magnetic flux sets up a persistent circulating current in the ring. For a Josephson dipole in such a loop, the definition of the flux $\Phi_j$ will differ in Eqs. (14) and (15) as a function of the equilibrium. Due to this adjustment, terms linear in $\Phi_j$ are absent from Eq. (15) by construction. Supplementary Sections A8 and A9 discuss these equilibrium considerations in detail.

**Quantum Hamiltonian.** Having conceptually carved out non-linear contributions from the system Hamiltonian $\hat{H}_{\text{full}}$ and collected them in the set of $\mathcal{E}_j^{\text{nl}}$ functions, we define the linearized Josephson circuit to correspond to everything left over in the system. This linear circuit consists of the electromagnetic circuit external to the Josephson dipoles, combined with their linear inductances $L_j$. We will use the eigenmodes of the linearized circuit to explicitly construct $\hat{H}_{\text{full}}$. The eigenmode frequencies and field distributions are readily obtained using a conventional FE solver (see Supplementary Section C). The Hamiltonian of the linearized Josephson circuit can thus be expressed as (see Supplementary Section A5)

$$\hat{H}_{\text{lin}} = \sum_{m=1}^{M} \hbar \omega_m \hat{a}_m^{\dagger} \hat{a}_m, \quad (16)$$
where M is the number of modes addressed by the numerical simulation, $\omega_m$ is the solution eigenfrequency of mode m, and $\hat{a}_m$ the corresponding mode amplitude (annihilation operator), defined by the mode eigenvector. We emphasize that the frequencies $\omega_m$ will be significantly perturbed by the Lamb shifts $\Delta_m$ and should be seen as an intermediate parameter entering in the calculation of the rest of the nonlinear Hamiltonian,

$$ \hat{H}_{nl} = \sum_{j=1}^{J} \varepsilon_j^{nl} = \sum_{j=1}^{J} E_j (c_{j3} \hat{\phi}_j^3 + c_{j4} \hat{\phi}_j^4 + \dots) \quad (17) $$

$$ = \sum_{j=1}^{J} E_j \sum_{p=3}^{\infty} c_{jp} \hat{\phi}_j^p, \quad (18) $$

$$ \hat{\phi}_j = \sum_{m=1}^{M} \phi_{mj} (\hat{a}_m^\dagger + \hat{a}_m), \quad (19) $$

where J is the total number of junctions and $\hat{\phi}_j := \hat{\Phi}_j/\phi_0$. In Eq. (17), we have introduced a Taylor expansion of $\varepsilon_j^{nl}$, where the energy $E_j$ and expansion coefficients $c_{jp}$ are known from the fabrication of the Josephson circuit, see Fig. 3c. For example, for a Josephson junction, the constant $E_j$ is just the Josephson energy, while $c_{jp}$ are the coefficients of the cosine expansion; i.e., $c_{jp}$ is 0 for odd p and $(-1)^{p/2+1}/p!$ for even p. The expansion is helpful for analytics but does not need to be used in the numerical analysis of $\hat{H}_{nl}$, see Supplementary Section A.

The Hamiltonian $\hat{H}_{full}$ is specified since the operators $\hat{\phi}_j$ are expressed in terms of the mode amplitudes as a linear combination (see Supplementary Section A5). Here, $\phi_{mj}$ are the dimensionless, real-valued, quantum zero-point fluctuations of the reduced flux of junction j in mode m. Determination of $\hat{H}_{full}$ is now reduced to computing $\phi_{mj}$. We achieve this by employing a generalization of the EPR.

The EPR $p_{mj}$ of junction j in eigenmode m is defined to be the fraction of inductive energy stored in the junction when only that mode is excited,

$$ \begin{aligned} p_{mj} &:= \frac{\text{Inductive energy stored in junction } j}{\text{Inductive energy stored in mode } m} \\ &= \frac{\langle \psi_m | \frac{1}{2} E_j \hat{\phi}_j^2 | \psi_m \rangle}{\langle \psi_m | \frac{1}{2} \hat{H}_{\text{lin}} | \psi_m \rangle}, \end{aligned} \quad (20) $$

which is a straightforward extension of Eq. (5), and is similarly computed using normal ordering (see Supplementary Section A6). The EPR $p_{mj}$ is computed from the eigenfield solutions $\hat{E}_m(\vec{r})$ and $\hat{H}_m(\vec{r})$ as explained in Supplementary Section C3. It is a bounded, non-negative, real number, $0 \le p_{mj} \le 1$. A zero EPR $p_{mj} = 0$ means that junction j is not excited in mode m. A unity EPR $p_{mj} = 1$ means that junction j is the only inductive element excited in the mode.

From the EPR $p_{mj}$, one directly computes the variance of the quantum zero-point fluctuations,

$$ \phi_{mj}^2 = p_{mj} \frac{\hbar \omega_m}{2E_j} \quad (21) $$

Equation (21) constitutes the bridge between the classical solution of the linearized Josephson circuit and the quantum Hamiltonian $\hat{H}_{full}$ of the full Josephson system, up to the sign of $\phi_{mj}$.

**Universal EPR properties:** The quantum fluctuations $\phi_{mj}$ are not independent of each other, since the EPRs are submitted to three types of universal constraints—valid regardless of the circuit topology and nature of the Josephson dipoles. These are of practical importance, as they are useful guides in evaluating the performance of possible designs and assessing their limitations. As shown in Supplementary Section A7, the EPRs obey one sum rule

per junction j and one set of inequalities per mode m,

$$ \sum_{m=1}^{M} p_{mj} = 1 \quad \text{and} \quad 0 \le \sum_{j=1}^{J} p_{mj} \le 1. \quad (22) $$

The total EPR of a Josephson dipole is a quantity that is independent of the number of modes—it is precisely unity for all circuits in which the dipole is embedded. It can only be diluted among the modes. On the other hand, a given mode can accept at most a total EPR of unity from all the dipoles. In practice, this sum rule can be fully exploited only if the bound M reaches the total number of relevant modes of the system.

The next fundamental property concerns the orthogonality of the EPRs. Rewriting Eq. (21) in terms of the amplitude of the zero-point fluctuation, we have

$$ \varphi_{mj} = s_{mj} \sqrt{p_{mj} \hbar \omega_m / 2E_j}, \quad (23) $$

where the EPR sign $s_{mj}$ of junction j in mode m is either +1 or -1. The EPR sign encodes the relative direction of current flowing across the junction. Only the relative value between $s_{mj}$ and $s_{mj'}$ for $j \ne j'$ has physical significance (see Supplementary Fig. 8). The EPR sign $s_{mj}$ is calculated in parallel with the process of calculating $p_{mj}$ from the field solution $\hat{H}(\vec{r})$, see Supplementary Section C3. We now obtain the EPR orthogonality relationship

$$ \sum_{m=1}^{M} s_{mj} s_{mj'} \sqrt{p_{mj} p_{mj'}} = 0, \quad (24) $$

valid when the sum from 1 to M covers all the relevant modes, see Supplementary Section A7.

**Excitation-number-conserving interactions:** Thus, as announced, knowledge of the EPRs completely specifies $\hat{H}_{nl}$, through Eqs. (17), (19), and (21). The Hamiltonian can now be analytically or numerically diagonalized using various computational techniques²⁰. In this section, our focus will be now to explicitly handle the effect of the nonlinear interactions $\hat{H}_{nl}$ on the eigenmodes. Before treating the case of a general nonlinear interaction, we focus on the leading-order effect of $\hat{H}_{nl}$ in the case of the “canonical” Josephson system. In this case, the J Josephson dipoles are all Josephson tunnel junctions and the dispersive regime is satisfied for all pairs of modes k and m; i.e., $\omega_k - \omega_m \gg E_j c_{jp} < \hat{\phi}_j^p$ for $p \ge 3$ and in the absence of strong drives. The leading-order nonlinear terms are the subset of p = 4 terms that conserve excitation number. After normal ordering, see Supplementary Section B1, one finds the effective Hamiltonian

$$ \hat{H}_4 = -\hbar \sum_{m=1}^{M} \Delta_m \hat{a}_m^\dagger \hat{a}_m + \frac{a_m}{2} \hat{a}_m^{\dagger 2} \hat{a}_m^2 + \sum_{n<m} \chi_{mn} \hat{a}_m^\dagger \hat{a}_n^\dagger \hat{a}_m \hat{a}_n, \quad (25) $$

which is a generalization of the one found in Eq. (8). In Eq. (25), we have introduced the Lamb shift $\Delta_m$ of mode m, the anharmonicity $a_m$ of the mode, and its total dispersive shift $\chi_{mn}$ (so-called cross-Kerr term) with a different mode, labeled n. Each of these parameters is directly calculated from the EPRs. As shown in Supplementary Section B, for arbitrary m and n,

$$ \chi_{mn} = \sum_{j=1}^{J} \frac{\hbar \omega_m \omega_n}{4E_j} p_{mj} p_{nj}, \quad (26) $$

while $a_m = \chi_{mm}/2$ and $\Delta_m = \sum_{n=1}^{M} \chi_{mn}/2$. Equation (26) implements, mathematically, the idea that the amplitude of these nonlinear couplings is the result of a spatial-mode scalar product of the EPRs. Remarkably, from Eq. (26) it is seen that the EPRs are essentially the only free parameters subject to design, when determining the nonlinear couplings, since $\omega_m, \omega_n$, and $E_j$ are generally tightly constrained by experimental considerations.
Equation (26) can be cast in matrix form by introducing the EPR matrix

$$ \mathbf{P} := \begin{pmatrix} p_{11} & \cdots & p_{1J} \\ \vdots & \ddots & \vdots \\ p_{M1} & \cdots & p_{MJ} \end{pmatrix}, \qquad (27) $$

which we have found useful in handling large circuits, especially for those in excess of 100 modes. We also introduce the diagonal matrices of eigenfrequencies $\Omega := \langle \omega_1, \dots, \omega_M \rangle$ and junction energies $\mathbf{E}_J := \langle E_1, \dots, E_J \rangle$, which lead to the matrix form of Eq. (26),

Kerr matrix: $\chi = \frac{\hbar}{4} \Omega \mathbf{PE}_J^{-1} \mathbf{P} \Omega$.

$$ \text{Anharmonicity: } a_m = \frac{1}{2} [\mathbf{X}]_{mm'}, \qquad (28) $$

$$ \text{Lamb shift: } \Delta_m = \frac{1}{2} \sum_{m'=1}^{M} [\mathbf{X}]_{mm'}. $$

We have defined the symmetric matrix of dispersive shifts $\chi$, with elements $[\mathbf{X}]_{mm'} = \chi_{mm'}$. Further discussion of the matrix approach and applications to pth-order corrections is deferred to Supplementary Section B2, and the amplitude of an arbitrary multi-photon interaction stemming from the full $\hat{H}_{nl}$ is calculated in Supplementary Section B3.

**EPR for dissipation in the circuit:** The EPR method treats the calculation of Hamiltonian and dissipation parameters on equal footing. Unlike in the impedance method⁴, one can completely characterize both $\hat{H}_{full}$ and the effect of dissipative elements in the circuit from the eigenfield solutions, $\hat{E}_m(\vec{r})$ and $\hat{H}_m(\vec{r})$. The list of dissipative elements include bulk and surface dielectrics⁸¹⁻⁸³, thin-film metals⁵³,⁸⁴, surface interfaces⁸⁵⁻⁸⁹, and metal seams⁹⁰. The EPR of a dissipative element $l$ in mode $m$ will be denoted $p_{ml}$. It is calculated similarly to $p_{mj}$ as summarized in Supplementary Section D. The participation $p_{ml}$ and the quality factor $Q_l$ of the material of this element are used to estimate the total quality factor of mode $m$ in the standard way when the fields are not greatly altered by the dissipation ($Q_m \gg 1)^{30,91-93}$,

$$ Q_m^{-1} = \sum_l p_{ml} Q_l^{-1}. \qquad (29) $$

Experimental values of $Q_l$ are found in the literature, and some are provided in Supplementary Section D. Equation (29) and the dissipative EPR $p_{ml}$ provide a dissipation budget for the individual influence of each dissipation mechanism in the system, providing a useful tool to optimize design layout for quantum coherence⁹⁴.

## Comparison between theory and experiment

Applying the EPR method, we designed eight superconducting samples to test the agreement between the EPR theory and experimental results. We tested several sample configurations, comprising 15 qubits, 8 cavity modes, and 1 waveguide in 3 different circuit quantum electrodynamics (cQED) architectures. The samples were measured in a standard cQED set-up, see "Methods," at the 15 mK stage of a dilution unit, over multiple cool downs.

Six of the samples were each composed of two qubits and one 3D cavity, one sample was composed of two qubits and a waveguide, and one sample was a flip-chip, 2.5D system²⁸ consisting of a flip-chip qubit embedded in a two-mode whispering gallery mode resonator⁵³ (WGMR). The specifics of each sample are discussed in the "Methods" section.

For each sample, we measured the circuit parameters of interest: dressed mode frequencies $\omega_m - \Delta_{m}$, anharmonicities of qubits and high-Q cavities $a_m$, cross-Kerr frequencies $\chi_{mn}$, and input-output (I-O) quality factors $Q_c$ for any readout modes. Our measurement methodology is detailed in the "Methods" section.

The measured parameters were compared to those calculated using the energy-participation method. The linearized Josephson circuit of each sample was modeled in Ansys High-Frequency Electromagnetic-Field Simulator (HFSS). Junctions were modeled as lumped inductors, whose nominal energy $E_j$ was inferred from room-temperature resistance measurements⁴⁹. To account for the error bars of the measurement and the drift in resistance over time, $E_j$ was adjusted by no more than 10% to fit the measured qubit frequency. To minimize the number of free parameters, we neglect the small junction intrinsic capacitance $C_j$ in our modeling. The tradeoff is a small and estimable systematic offset of the bare simulated mode anharmonicities. We estimate this correction to be on the order of 4% for a $C_j = 4$ ff. From the eigenfield solutions, we calculated the EPRs $p_{mj}$ and the sign $s_{mj}$ to construct $H$ and extract its parameters. Detailed steps of the procedure can be found in Supplementary Section C. The results are presented in Tables 1–3.

**Figure 4** summarizes the agreement of the measured and calculated sample parameters, which span five orders of magnitude in frequency. Accounting for $C_j$, we find that mode frequencies are calculated to 1% accuracy, large nonlinear interaction energies (namely, anharmonicity and cross-Kerr frequencies >10 MHz) are calculated at the 5% level, and small nonlinear interaction energies agree at the 10% level. We highlight that we have used minimal, coarse adjustment to account for shifts in $E_j$ and otherwise, by neglecting $C_j$, the calculation is free from adjustable parameters.

The results of Fig. 4 demonstrate the accuracy and applicability of the EPR method. For each device, the EPR results are obtained from a single eigenmode simulation, using full automation of the analysis, provided by our open-source package *pyEPR* (see the *pyEPR*⁹⁵ code repository at http://github.com/zlatko-minev/pyEPR). For current standard applications, we find the agreement sufficient. Further improvements in accuracy would require improved ability to estimate the Josephson dipole energy $E_j$ and its intrinsic capacitance $C_j$. At the same level of accuracy, improvements in the precision and reproducibility of the implementation and assembly of the Josephson circuit design are needed, such as in chip-clamping techniques, precision machining of the device sample holder, and I-O couplers.

**Conclusion.** An intuitive, easy-to-use, and efficient method is needed to design and analyze Josephson microwave quantum circuits. We have described in this article such a method, based on the distribution of the electromagnetic energy in the circuit and its participation in nonlinear and dissipative elements. This so-called EPR method offers physical insight helping the design process and provides a simple link between the classical circuit and its quantum properties. By comparing our theory to eight experimental devices incorporating Josephson junctions, we have shown that our method is accurate and applicable to a large range of quantum circuit architectures. It is directly applicable to a broader class of nonlinear inductive elements, such as weak-link nanobridges³⁸,⁶³, nanowires⁴⁰,⁴¹,⁶⁴,⁶⁶,⁶⁷, and kinetic-inductance thin films⁴³,⁵⁸,⁵⁹. While best suited for weakly nonlinear systems, the EPR method is derived within circuit theory without approximations. It can be seen as arising from a change of basis adapted to nonlinear elements, as detailed in Supplementary Section A. In practice, the useful reach of the method is set by the numerical ability to include all relevant electromagnetic modes and to compute the spectrum of the extracted Hamiltonian²⁰. We contribute an open-source package *pyEPR* (see the *pyEPR*⁹⁵ code repository at http://github.com/zlatko-minev/pyEPR), which automates the EPR method, and was tested in the design of several further experiments²⁸,⁷⁴,⁹⁶⁻¹⁰⁵.

## METHODS

### Methods of the experiment

**Device fabrication.** Unless otherwise noted, samples were fabricated according to the following methodology. Sample patterns, both large and fine features, were defined by a 100 kV electron-beam pattern generator
<table><caption>**Table 1.** Two-qubit, one-cavity devices.</caption><thead><tr><th rowspan="2">Device</th><th colspan="3">Frequency (MHz)</th><th colspan="2">Anharmonicity (MHz)</th><th colspan="3">Cross-Kerr (MHz)</th><th rowspan="2">I-O coupling<br>Q<sub>C</sub></th></tr><tr><th>ω<sub>D</sub>/2π</th><th>ω<sub>B</sub>/2π</th><th>ω<sub>C</sub>/2π</th><th>α<sub>D</sub>/2π</th><th>α<sub>B</sub>/2π</th><th>χ<sub>DB</sub>/2π</th><th>χ<sub>BC</sub>/2π</th><th>χ<sub>DC</sub>/2π</th></tr></thead><tbody><tr><td rowspan="3">R9C1</td><td>4951</td><td>5664</td><td>9158</td><td>138</td><td>170</td><td>92</td><td>4.7</td><td>0.4</td><td>5.20 × 10<sup>3</sup></td></tr><tr><td>4866</td><td>5691</td><td>9154</td><td>150</td><td>185</td><td>99</td><td>4.2</td><td>0.55</td><td>7.40 × 10<sup>3</sup></td></tr><tr><td>-1.7%</td><td>0.5%</td><td>-0.04%</td><td>8%</td><td>8%</td><td>7%</td><td>-12%</td><td>27%</td><td>29%</td></tr><tr><td rowspan="3">R2C1</td><td>4823</td><td>5567</td><td>8947</td><td>150</td><td>192</td><td>64.5</td><td>4.8</td><td>0.3</td><td>4.97 × 10<sup>3</sup></td></tr><tr><td>4770</td><td>5640</td><td>8950</td><td>161</td><td>211</td><td>67.7</td><td>5.88</td><td>0.46</td><td>5.44 × 10<sup>3</sup></td></tr><tr><td>-1.1%</td><td>1.3%</td><td>0.03%</td><td>6.8%</td><td>9%</td><td>4.7%</td><td>18%</td><td>35%</td><td>9%</td></tr><tr><td rowspan="3">R7C1</td><td>4726</td><td>5475</td><td>8999</td><td>156</td><td>189</td><td>67</td><td>4.8</td><td>0.34</td><td>2.68 × 10<sup>3</sup></td></tr><tr><td>4770</td><td>5640</td><td>8950</td><td>161</td><td>211</td><td>67.7</td><td>5.88</td><td>0.46</td><td>3.07 × 10<sup>3</sup></td></tr><tr><td>0.9%</td><td>2.9%</td><td>-0.55%</td><td>3.1%</td><td>10%</td><td>1%</td><td>18%</td><td>26%</td><td>13%</td></tr><tr><td rowspan="3">R3C2</td><td>4845</td><td>5620</td><td>8979</td><td>152</td><td>195</td><td>61</td><td>5.1</td><td>0.3</td><td>2.11 × 10<sup>3</sup></td></tr><tr><td>4770</td><td>5640</td><td>8950</td><td>161</td><td>211</td><td>67.7</td><td>5.88</td><td>0.46</td><td>1.78 × 10<sup>3</sup></td></tr><tr><td>-1.5%</td><td>0.4%</td><td>-0.3%</td><td>5.6%</td><td>7.6%</td><td>9.9%</td><td>13%</td><td>35%</td><td>-19%</td></tr><tr><td rowspan="3">R3C1</td><td>4688</td><td>5300</td><td>9003</td><td>148</td><td>174</td><td>85</td><td>5</td><td>0.33</td><td>2.43 × 10<sup>3</sup></td></tr><tr><td>4745</td><td>5265</td><td>8922</td><td>159</td><td>198</td><td>73</td><td>5.1</td><td>0.37</td><td>5.65 × 10<sup>3</sup></td></tr><tr><td>1.2%</td><td>-0.7%</td><td>-0.9%</td><td>6.9%</td><td>12.1%</td><td>-16%</td><td>2%</td><td>9%</td><td>57%</td></tr><tr><td rowspan="3">DT3</td><td>6160</td><td>7110</td><td>9170</td><td>130</td><td>150</td><td>278</td><td>3</td><td>2.5</td><td>9.17 × 10<sup>3</sup></td></tr><tr><td>6100</td><td>7141</td><td>9155</td><td>140</td><td>177</td><td>312</td><td>3.9</td><td>3.1</td><td>7.33 × 10<sup>3</sup></td></tr><tr><td>-1.0%</td><td>0.4%</td><td>-0.15%</td><td>7%</td><td>15%</td><td>11%</td><td>23%</td><td>19%</td><td>-25%</td></tr></tbody></table>

Summary of measured and calculated Hamiltonian and input-output (I-O) coupling parameters for the six devices described in “Methods.” Indices D, B, and C denote the dark, bright, and cavity modes, respectively. The input-output quality factor to the readout cavity is denoted $Q_C$. For each device, the first (second) row quantifies the measured, *m*, (bare calculated, *c*) values. The third row quantifies the bare agreement, i.e., $(c - m)/c$. In the anharmonicity column, the bare agreement should be corrected by the systematic shift due to our choice to neglect the junction intrinsic capacitance in our modeling (see “Methods”). We evaluate the correction to be of order 4%, estimated by taking a nominal junction $C_f = 4$ fF; hence, an overall corrected agreement of 4.3% for this column.

<table><caption>**Table 2.** Flip-chip (2.5D), one-qubit, one-storage-cavity, one-readout-cavity devices.</caption><thead><tr><th rowspan="2">Device</th><th colspan="3">Frequency (MHz)</th><th colspan="2">Anharmonicity (MHz)</th><th colspan="2">Cross-Kerr (MHz)</th><th rowspan="2">I-O coupling<br>Q<sub>C</sub></th></tr><tr><th>ω<sub>O</sub>/2π</th><th>ω<sub>S</sub>/2π</th><th>ω<sub>C</sub>/2π</th><th>α<sub>D</sub>/2π</th><th>α<sub>B</sub>/2π</th><th>χ<sub>QS</sub>/2π</th><th>χ<sub>QC</sub>/2π</th></tr></thead><tbody><tr><td rowspan="3">WG1</td><td>4890</td><td>7070</td><td>7267</td><td>310</td><td></td><td>0.25</td><td>0.30</td><td>20 × 10<sup>3</sup></td></tr><tr><td>4820</td><td>7020</td><td>7340</td><td>325</td><td></td><td>0.29</td><td>0.33</td><td>16 × 10<sup>3</sup></td></tr><tr><td>-1.4%</td><td>-0.7%</td><td>1.0%</td><td>4.6%</td><td></td><td>13%</td><td>9%</td><td>-22%</td></tr></tbody></table>

Summary of measured and calculated Hamiltonian and input-output (I-O) coupling parameters for the device described in “Methods.” Indices Q, S, and C denote the qubit, storage, and readout cavity modes, respectively. The input-output quality factor to the readout cavity is denoted $Q_C$. For each device, the first (second) row quantifies the measured, *m*, (bare calculated, *c*) values. The third row quantifies the bare agreement, i.e., $(c - m)/c$.

<table><caption>**Table 3.** Two-qubit, one-waveguide devices.</caption><thead><tr><th>Frequency (MHz)</th><th colspan="2">Anharmonicity (MHz)</th><th>Cross-Kerr (MHz)</th></tr><tr><th>ω<sub>D</sub>/2π</th><th>ω<sub>B</sub>/2π</th><th>α<sub>D</sub>/2π</th><th>α<sub>B</sub>/2π</th><th>χ<sub>DB</sub>/2π</th></tr></thead><tbody><tr><td>6010</td><td>8670</td><td>85</td><td>180</td><td>278</td></tr><tr><td>5824</td><td>8878</td><td>97</td><td>206</td><td>281</td></tr><tr><td>-3.2%</td><td>2.3%</td><td>12%</td><td>13%</td><td>1.1%</td></tr></tbody></table>

Summary of measured and calculated Hamiltonian parameters for the device described in “Methods.” Indices D and B denote the dark and bright modes, respectively. For each device, the first (second) row summarizes the measured, *m*, (bare calculated, *c*) values. The third row quantifies the bare agreement, i.e., $(c - m)/c$.

(Raith EBPG 5000+) in a single step on a PMAA/MAA (Microchem A-4/ Microchem EL-13) resist bilayer coated on a 430 µm thick, double-side-polished, c-plane sapphire wafer, grown with the edge-defined film-fed growth technique. Using the bridge-free fabrication technique$^{106-108}$, the

Al/AlOx/Al Josephson tunnel junctions were formed by a double-angle aluminum evaporation under ultra-high vacuum in a multi-chamber Plassys UMS300 UHV. The two depositions were interrupted by a thermal oxidation step, static 100 Torr environment of 85% argon and 15% oxygen, to form the thin AlOx barrier of the tunnel junction. Prior to the first deposition, to reduce junction aging$^{108}$, the exposed wafer surfaces were exposed to 1 min oxygen-argon plasma cleaning, under a pressure of 3 × 10⁻³ mbar. After wafer dicing (ADT ProVecturs 7100) and chip cleaning, the normal-state resistance $R_N$ of the Josephson junctions was measured to provide an estimate of the Josephson energy, $E_J$, of the device junctions. The junction energy was to first order estimated by an extrapolation of $R_N$ from room temperature to the operating sample temperature, at approximately 15 mK, using the Ambegaokar-Baratoff relation$^{109}$,

$$E_J = \frac{h\Delta}{2(2e)^2} R_N^{-1}, \qquad (30)$$

where Δ is the superconducting gap of aluminum, *e* is the elementary charge, and *h* is Planck's constant.

**Sample holder.** Sample holders were machined in aluminum alloy 6061; seams were formed using thin indium gaskets placed in machined grooves
<figure><img src="image_7.png" /></figure>

**Fig. 4** Theory vs. experiment. Comparison between theory and experiment over five-orders of magnitude in energy scale of the system Hamiltonian $\hat{H}_{\text{full}}$ for eight distinct, multi-mode device samples, described in detail in the “Methods” section, including 3D, flip-chip (2.5D), and 3D waveguide architectures incorporating readout and storage resonators and qubit modes. For each device, the dominant parameters in $\hat{H}_{\text{full}}$, dressed frequencies $\omega_m$, bare anharmonicities $a_m$, and cross-Kerr interactions $\chi_{mn}$ were measured and calculated using the EPR method with our open-source pyEPR package (see the pyEPR⁹⁵ code repository at http://github.com/zlatko-minev/pyEPR). Gray line is of slope one, representing ideal agreement between theory and experiment.

in one of the mating surfaces. Only non-magnetic components were used in proximity to the samples, molybdenum washers, aircraft-alloy 7075 screws (McMaster/Fastener Express) with <1% iron impurities, and non-magnetic SubMiniature version A (SMA) connectors.

**Cryogenic set-up.** Samples were thermally anchored to the 15 mK stage of a cryogen-free dilution refrigerator (Oxford Triton 200) and were measured using a standard cQED measurement set-up⁴⁵,⁵¹,⁹⁶. High-magnetic-permeability, μ-metal (Amumetal A4K) shields together with aluminum superconducting shields enclosed all samples. Microwave input and output lines were filtered with Eccosorb CR-110 infrared-frequency filters⁵⁵,⁹², thermally anchored at the 15 mK stage. Output lines were additionally filtered with cryogenic isolators (Quinstar CWJ1019-K414) and 12 GHz K&L multi-section lowpass filters. Output lines leading up to the high-electron-mobility transistor (HEMT) amplifier (Low Noise Factory), anchored at 4 K, were superconducting (CoaxCo Ltd. SC-086/50-NbTi-NbTi PTFE).

**Quantum amplifier.** The output signal of a sample was processed by a Josephson parametric converter (JPC) anchored at the 15 mK stage and operated in amplification mode¹¹⁰,¹¹¹, before routing to the HEMT. The JPC provide a typical gain of 21 dB with a typical noise-visibility ratio of 6 dB. See ref. ¹¹² for a review of the parametric amplification.

**Frequency and I-O coupling measurements.** Spectroscopic measurements were used to determine the frequencies of the resonator modes. Anharmonicities were determined in two-tone spectroscopy⁹²,¹¹³. Cross-Kerr energies were determined from dressed dephasing measurements¹¹⁴,¹¹⁵. In particular, the dressed-dephasing measurement sequence consisted of first preparing the qubit in the ground state, then exciting it to the equator by a $\pi/2$ pulse. Subsequently, a weak readout tone excited the readout cavity of the qubit for a fixed duration, 10 times the readout cavity lifetime $\kappa_r$, after which we measure the qubit X and Y Bloch vectors, after waiting for a time $5/\kappa_r$ for any photons in the cavity to leak out. By varying the amplitude and frequency of the applied weak-readout tone, we could calibrate both the strength of our readout, in steady-state photon number in the readout cavity, and the value of the cross-Kerr frequency shift between the qubit and readout resonator. The values could be obtained

<figure><img src="image_8.png" /></figure>

**Fig. 5** Two-qubit, one-cavity devices. **a, b** Not-to-scale diagram illustrating chip configurations A and B, respectively. Vertical blue arrow indicates cavity electric field orientation. Crosses mark the location of Josephson tunnel junctions. **c** Optical photograph of sample R2C1. Bottom half of aluminum sample holder is visible; top half is removed. The two-qubit chip (outlined by the dashed green box) is housed in the middle of the readout cavity (highlighted in blue). Cavity fundamental mode electric field profile E depicted by arrows. *Input-output* SMA pin coupler labeled I-O.

from fits of the X and Y quadratures. For each sample, the coupling quality factor of the readout-cavity mode, denoted $Q_c$, was extracted from the spectroscopic response of the readout cavity at low photon numbers⁹²,¹¹³, by measuring the scattering parameters, $S_{21}$ or $S_{11}$.

To test EPR's robustness to experimental variability and its applicability over wide range of experimental conditions, the presented samples were fabricated in multiple runs and measured in different cooldowns. Some devices were subjected to as many as six thermal cycles.

The Hamiltonian parameters and coupling energies for each sample were also calculated, following the EPR method presented in section on the general approach. In particular, we modeled the sample geometry and materials in a FE electromagnetic simulation, as explicated in Supplementary Section C. Our aim in writing this supplementary section has been to provide an easy access point to the practical use of the EPR method, which we hope will benefit the reader, and allow them to adopt it easily. Our choice of simulation software was the Ansys HFSS, although we emphasize that the EPR ideas translate to any standard EM eigenmode simulation package. Further, we modeled the loss due to the I-O couplers in the simulation as 50 Ω resistive sheets, see Supplementary Section D4. The eigenmode analysis provided the calculated I-O quality factors and Purcell limits. All electromagnetic and quantum analyses, including the extraction of participations from the eigenfields and the numerical diagonalization of the Hamiltonian to extract its quantum spectrum, were performed in a fully automated manner using the freely available pyEPR package (see the pyEPR⁹⁵ code repository at http://github.com/zlatko-minev/pyEPR).

The mode quality due to the I-O coupling, $Q_c$, was set by the length of the I-O SMA-coupler pin. Its length inside the sample-holder box was measured at room temperature using calipers. This nominal length was used then used in the HFSS model to create a 3D model of the pin inside the sample holder. The quality factor $Q_c$ was then obtained from the eigenmode eigenvalue. We remark that the measurement of the pin length is accurate to no more than 20%; further, it can be affected by various idiosyncrasies, such as bending of the thin SMA center pin. Nonetheless, the predictions of the quality factors for low-Q modes were observed to be very reasonable estimates, and similarly, the predicted Purcell limits for qubit and high-Q cavity modes were consistent with estimates from measurements.

## Devices: two-qubit, one-cavity devices

**Device description.** We measured six samples that were each comprised of two qubits and one cavity. The cavity was a standard, machined aluminum cavity⁵⁴. It housed either one or two sapphire chips, which were either patterned with transmon qubits or simply blank. Each transmon consisted of two thin-film aluminum pads connected by a Josephson junction. We tested two configurations of chips and patterns. Configuration A consisted of one chip with two orthogonal qubits, as depicted in Fig. 5a. Similarly, configuration B consisted of one chip with two parallel qubits, depicted in Fig. 5b. The two qubits were aligned parallel to each other; however, unlike configuration B, there was no galvanic connection between them. The results of the measurements are presented in Table 1.
Samples R1C9, R2C1, R7C1, R2C1, and R3C2 were fabricated in configuration A; sample DT3 was fabricated in configuration B. Three of the sample (R2C1, R7C1, R3C1) were fabricated simultaneously on the same sapphire wafer, all with nominally identical dimensions. Additionally, R2C1 and R7C1 were designed to have nominally the same Josephson junctions energy, $E_J$. The rest of the samples (R1C9, DT3, and R3C2) were fabricated at different times and on different wafers. The dimensions of their transmons and the inductance of the junctions were designed to be different. Only sample R3C2 was designed to be very similar to the nominally identical sample R2C1 and R7C1 but with adjusted $E_J$. For samples R2C1, R7C1, R3C1, and R3C2 a second, un-processed, un-patterned, blank sapphire chip was placed in parallel with the qubit carrying chip [see Fig. 5c] to purposefully lower the readout cavity frequency, thus bringing it within the JPC amplification band.

Configurations A and B were designed to test the ability of the EPR method to calculate the mixing between strongly coupled modes. The strong coupling was achieved in two distinct ways. First, configuration A used the spatial proximity of the two qubits to yield a strong capacitive coupling between them, which resulted in large qubit-qubit mixing. Second, instead of spatial proximity, configuration B used a galvanic connection between the qubits to yield strong hybridization. Our two-qubit designs share some similar-in-spirit characteristics with the promising recent developments reported in refs. 116–121, but our implementation is distinct and is designed to provide several unique advantages.

**Mode structure and interesting physical insights.** Configuration A is characterized by strong capacitive coupling between the two transmons, which have different pad sizes, see Fig. 5c, and hence different normal-mode frequencies. Due to the strong hybridization, each qubit normal mode consists of some excitation in the vertical and some in the horizontal transmon. With some foresight, we will label the vertical mode *bright* (B) and the horizontal *dark* (D). The bright-mode resonance is higher in frequency and thus is closer to the resonance of the readout cavity mode (C). This smaller detuning made it a natural choice for designing stronger coupling between it, (B), and the readout mode (C). This was implemented by orienting the transmon design that participates in mode (B) vertical.

To understand this design choice, let us first consider the popular analogy46,122 between cQED and cavity QED, often used to discuss mode couplings. In this atomic analogy, the transmon qubit is analogous to a real atom inside the cavity. Thus, it can be described by an electric-dipole moment $d_B$. Meanwhile, its coupling, cross-Kerr, etc. to the cavity mode are derived from the electric-dipole coupling interaction. In particular, the coupling amplitude is proportional to $d_B \cdot E$, where $E$ is the cavity electric field at the transmon junction. From this analogy, one can infer that the coupling is maximized when the two are parallel, $d_B \parallel E$, and one could hope to measure a strong cross-Kerr between the bright qubit and the cavity. This successful conclusion is true, but a coincidence. We will shortly discuss how this popular analogy fails spectacularly for the dark mode in configuration B. Instead, we will argue that a correct way to understand the nonlinear coupling between the two modes is through the participation ratio, which will provide the correct coupling for both configurations A and B.

Before proceeding to configuration B, we note one further useful features that configuration A exhibits. In particular, while the bright qubit mode can be Purcell limited116,123, the dark mode is simultaneously Purcell protected. Thus, one can potentially achieve a high ratio in the I–O bath coupling of the two qubits.

Configuration B has two qubit modes, which we will also label *dark* (D) and *bright* (B). Since both transmons are designed with the exact same transmon pad geometry and junction energy $E_J$, see Fig. 5, we can expect that no single junction is preferred, due to the symmetry of the sample. This is in sharp contrast to the asymmetric energy distribution in configuration A. Returning to configuration B, we can estimate that, in each qubit mode, both junctions participate equally and with near maximal allowed participation,

$$p_{D1} = p_{D2} = p_{B1} = p_{B2} \approx \frac{1}{2}. \qquad (31)$$

If the two transmons were well-separated spatially and not connected, they would be uncoupled. However, the galvanic connection between the two lower pads, see Fig. 5a, results in a very strong hybridization and splitting between the nominally identical transmons. The result of the strong hybridization is a symmetric and antisymmetric combination of the two bare transmons. In other words, the hybridization results in a common mode, namely, (B), where both junctions oscillate in phase,

and a differential mode, namely, (D), where both junctions oscillate out of phase. These phase relationships are captured by the signs:

$$s_{D1} = 1, \quad s_{D2} = -1, \qquad (32)$$

$$s_{B1} = 1, \quad s_{B2} = +1. \qquad (33)$$

In an attempt to understand how these hybridized qubit modes will couple to the cavity mode (C), let us first consider the atomic analogy again. When the two junctions oscillate in phase, in the (B) mode, the net dipole moment of the bright mode, $\bar{d}_B$, must be large, since it is the sum of the two junction dipole contributions. Second, $\bar{d}_B$ must be oriented in the vertical direction, parallel to the cavity electric field $\bar{E}$. Hence, we would conclude that the bright mode coupling is large, $\bar{d}_B \cdot \bar{E} \gg 0$, and there should be a strong cross-Kerr interaction between the cavity and bright qubit. Continuing the analogy in the case of the dark mode, we would deduce that the net dipole moment of the dark mode is zero, since the two junctions oscillate out of phase and cancel each other’s contribution, $\bar{d}_D = 0$. Thus, we should not expect any coupling between the dark qubit and the cavity mode, $\bar{d}_D \cdot \bar{E} = 0$. To the contrary of this conclusion, as can be seen in the measured results in Table 1, the nonlinear coupling of the dark and bright qubit to the cavity is nearly equal. The atomic analogy and the dipole argument have failed completely. We can understand the origin of this failure and how to arrive at the correct conclusion by using the EPR. As embodied in Eq. (26), in the dispersive regime, the nonlinear coupling between two modes, in this case a qubit and cavity, is given by the overlap of the EPR distribution. In particular, the cross-Kerr amplitude between the dark qubit and the readout cavity mode is given by

$$\chi_{DC} = \frac{\hbar\omega_D\omega_C}{4E_J}(p_{D1}p_{C1} + p_{D2}p_{C2}), \qquad (34)$$

where both junctions have the same junction energy $E_J$, and $\omega_C$ (resp: $\omega_C$) denotes the dark qubit (resp: cavity) mode frequency. The signs, used in the atomic dipole logic, do not factor into the coupling, because the Josephson mechanics is fundamentally different. To obtain $\chi_{BC}$, one can replace the label “D” with “B” in Eq. 34. Then, it is easy to use Eq. (31) to show that the ratio of two Kerr couplings is not zero, but rather of order unity,

$$\chi_{BC}/\chi_{DC} = \omega_B/\omega_D. \qquad (35)$$

**Failure of the conventional dipole approach.** We showed that, although the heuristic atomic analogy seems seductively accurate, it fails completely in some cases to predict the nonlinear couplings. Instead, one can use the intuition and calculation method provided by the EPRs.

As an added note, we observe that Eqs. (32) and (33) embodies the orthogonality of the participations, see Eq. (24). We also remark that, although the atomic analogy fails in the case of the nonlinear couplings, it can yield some guidance when considering the linear mixing of the modes, useful for discussing the Purcell effect. To illustrate, let us briefly extend the atomic analogy. The dipole-like coupling between the bright mode and the cavity suggests that the bright mode will inherit some coupling to the environment, mediated by the cavity. Thus, since $\bar{d}_B \cdot \bar{E} \gg 0$, we can expect the bright qubit to potentially be Purcell limited. In contrast, since $\bar{d}_D \cdot \bar{E} = 0$, we could expect the dark qubit to be Purcell protected. Both of these qualitative Purcell predictions are valid, but to quantify them, we will use the EPR method and FE eigenmode simulation of the sample, as will be discussed shortly.

**Experimental results.** Table 1 summarizes the results of the agreement between the measured and calculated Hamiltonian parameters for all two-qubit, one-cavity samples. The three modes in each sample are labeled dark (D), bright (B), and cavity (C); the reason for this convention is described above. In all samples, the qubits were designed to be in the dispersive regime with respect to the cavity, which was detuned by 2–4 GHz. However, in a large number of the samples, the two qubits were strongly hybridized, often necessitating higher-order nonlinear corrections to be included in the calculation. This strong hybridization was used as a test of the theory in this more challenging and fickle regime.

In total, for each sample we measured and calculated eight frequency parameters and one dimensionless, coupling quality factor, $Q_C$, of the readout cavity mode. In particular, in the low-excitation limit, the nonlinear interactions among the modes were characterized by the effective
<figure><img src="image_9.png" /><figcaption>Fig. 6 Schematic representation (not to scale) of superconducting waveguide device DTW1. Two-qubit chip is placed λ/4 from a short termination in the waveguide. Guided input-output waves are launched and monitored through the I-O port connecting to a standard SMA adapter (not shown).</figcaption></figure>

dispersive Hamiltonian

$$ \begin{aligned} \hat{H}/\hbar ={}& \omega_D \hat{n}_D + \omega_B \hat{n}_B + \omega_C \hat{n}_C \\ & -\frac{1}{2} \alpha_D \hat{n}_D (\hat{n}_D - \hat{1}) - \frac{1}{2} \alpha_B \hat{n}_B (\hat{n}_B - \hat{1}) \\ & -\chi_{DB} \hat{n}_D \hat{n}_B - \chi_{DC} \hat{n}_D \hat{n}_C - \chi_{BC} \hat{n}_B \hat{n}_C, \end{aligned} \quad (36) $$

where $\hat{n}_D$, $\hat{n}_B$, and $\hat{n}_C$ denote the dark, bright, and cavity photon-number operator, respectively. The coupling of the resonator mode to the bath is given by the Lindblad superoperator term $\kappa_C D[\hat{a}_C]\rho$, where $\kappa_C = \omega_C/Q_C$, and $\rho$ is the density operator.

We remark that all samples in configuration A demonstrated a large asymmetry in the Kerr coupling between the bright-to-cavity and dark-to-cavity coupling, $\chi_{BC} \gg \chi_{DC}$. In contrast, samples in configuration B demonstrated near equal coupling, $\chi_{BC} \approx \chi_{DC}$. In both configurations A and B, the dark mode was Purcell protected, we calculated a Purcell coupling factor of $Q_{\text{Purcell}}^D \gg 10^7$, using the eigenmode method described in Supplementary Section D4. On the other hand, the bright mode was somewhat Purcell limited, $Q_{\text{Purcell}}^B \approx 10^6$. From the relative Rabi amplitudes of the dark and bright qubits, we could verify the order of magnitude scaling calculated for the Purcell effect.

### Devices: two-qubit, single-waveguide devices

We measured a two-qubit sample inside of a waveguide. Figure 6 presents the set-up and depicts the sample, which was of the configuration B type presented in Fig. 5. The sample chip was positioned inside an aluminum WR90 waveguide. The waveguide was terminated in a short at one side and attached to an impedance-matched SMA coupler port on the input-launcher side, which was used to drive and measure the waveguide. The chip was centered inside the cross-section of the waveguide and placed λ/4 away from the termination wall, at the measurement frequency. The rest of the experimental set-up was identical to that described in section “Methods of the experiment.” The two qubit modes were labeled dark and bright, similarly to the samples discussed in the section “Two-qubit, one-cavity devices.”

Table 3 presents the agreement between the measured and calculated key Hamiltonian parameters of the sample. These consist of the two-mode frequencies, two-qubit anharmonicities, and the strong cross-Kerr interaction between the two qubits.

### Devices: flip-chip (2.5D), one-qubit, one-storage-cavity, one-readout-cavity devices

We also designed a multilayer planar²⁸ (2.5D) cQED sample, depicted in Fig. 7, with the EPR method. It consisted of high-Q storage mode (S), one low-Q readout cavity (C), and one control transmon qubit (Q). The two cavity modes were formed in the footprint of a single WGMR⁵³. The three modes were in the dispersive regime, and the storage mode was used to encode and decode quantum information, as well as to observe parity revivals. Details of the sample design have been reported in ref. ²⁸. The agreement between the measured parameters of the sample and those obtained by the EPR calculation methods are presented in Table 2.

<figure><img src="image_10.png" /><figcaption>Fig. 7 Illustration of flip-chip (2.5D) device WG1. a Depiction (not to scale) of chip stack consisting of two chips separated by a 100 µm vacuum gap. The inner face of each chip supports part of the pattern of a multi-layer whispering gallery-mode resonator (WGMR) resonator⁵³. The lower layer contains a 2.5D, aperture transmon qubit²⁸ embedded in the WGMR. b Zoomed-in view of the lower layer of the aperture transmon. The cross marks the location of the Josephson tunnel junction device, which connects the lower center trace to an island embedded in the lower aperture. The capacitance to the top layer significantly participates in determining qubit parameters.</figcaption></figure>

## DATA AVAILABILITY

Data are available from the authors on reasonable request.

## CODE AVAILABILITY

The source code for the EPR method is open-sourced and can be found at http://github.com/zlatko-minev/pyEPR.

Received: 22 October 2020; Accepted: 17 June 2021; Published online: 25 August 2021

## REFERENCES

1. Devoret, M. H. & Schoelkopf, R. J. Superconducting circuits for quantum information: an outlook. Science 339, 1169–1174 (2013).

2. Arute, F. et al. Quantum supremacy using a programmable superconducting processor. Nature 574, 505–510 (2019).

3. Blais, A., Grimsmo, A. L., Girvin, S. M. & Wallraff, A. Circuit quantum electrodynamics. https://journals.aps.org/rmp/abstract/10.1103/RevModPhys.93.025005 (2020).

4. Nigg, S. E. et al. Black-box superconducting circuit quantization. Phys. Rev. Lett. 108, 240502 (2012).

5. 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).

6. Solgun, F., Abraham, D. W. & DiVincenzo, D. P. Blackbox quantization of superconducting circuits using exact impedance synthesis. Phys. Rev. B 90, 134504 (2014).

7. Solgun, F. & DiVincenzo, D. P. Multiport impedance quantization. Ann. Phys. 361, 605–669 (2015).

8. Smith, W. C. et al. Quantization of inductively shunted superconducting circuits. Phys. Rev. B 94, 144507 (2016).

9. Gely, M. F. et al. Convergence of the multimode quantum Rabi model of circuit quantum electrodynamics. Phys. Rev. B 95, 245115 (2017).

10. Malekakhlagh, M., Petrescu, A. & Türeci, H. E. Cutoff-free circuit quantum electrodynamics. Phys. Rev. Lett. 119, 073601 (2017).

11. Pechal, M. & Safavi-Naeini, A. H. Millimeter-wave interconnects for microwave-frequency quantum machines. Phys. Rev. A 96, 042305 (2017).

12. Parra-Rodriguez, A., Egusquiza, I. L., DiVincenzo, D. P. & Solano, E. Canonical circuit quantization with linear nonreciprocal devices. Phys. Rev. B 99, 014514 (2019).

13. Parra-Rodriguez, A., Rico, E., Solano, E. & Egusquiza, I. L. Quantum networks in divergence-free circuit QED. Quantum Sci. Technol. 3, 024012 (2018).

14. Ansari, M. H. Superconducting qubits beyond the dispersive regime. Phys. Rev. B 100, 024509 (2019).

15. Krupko, Y. et al. Kerr nonlinearity in a superconducting Josephson metamaterial. Phys. Rev. B 98, 094516 (2018).

16. Malekakhlagh, M., Petrescu, A. & Türeci, H. E. Lifetime renormalization of weakly anharmonic superconducting qubits. I. Role of number nonconserving terms. Phys. Rev. B 101, 134509 (2020).

17. Solgun, F., DiVincenzo, D. P. & Gambetta, J. M. Simple impedance response formulas for the dispersive interaction rates in the effective Hamiltonians of low
anharmonicity superconducting qubits. *IEEE Trans. Microwave Theory Tech.* **67**, 928–948 (2019).

18. Petrescu, A., Malekakhlagh, M. & Türeci, H. E. Lifetime renormalization of driven weakly anharmonic superconducting qubits: II. The readout problem. *Phys. Rev. B* **101**, 134510 (2019).

19. You, X., Sauls, J. A. & Koch, J. Circuit quantization in the presence of time-dependent external flux. *Phys. Rev. B* **99**, 174512 (2019).

20. Di Paolo, A., Baker, T. E., Foley, A., Sénéchal, D. & Blais, A. Efficient modeling of superconducting quantum circuits with tensor networks. Preprint at https://www.nature.com/articles/s41534-020-00352-4 (2019).

21. Menke, T. et al. Automated design of superconducting circuits and its application to 4-local couplers. *npj Quantum Inf.* **7**, 49 (2021).

22. Gely, M. F. & Steele, G. A. QuCAT: quantum circuit analyzer tool in Python. *N. J. Phys.* **22**, 013025 (2020).

23. Kyaw, T. H. et al. Quantum computer-aided design: digital quantum simulation of quantum processors. Preprint at http://arxiv.org/abs/2006.03070 (2020).

24. Yan, F. et al. Engineering framework for optimizing superconducting qubit designs (2020). Preprint at http://arxiv.org/abs/2006.04130. 2006.04130.

25. Minev, Z. K., McConkey, T. G., Takita, M., Corcoles, A. D. & Gambetta, J. M. Circuit quantum electrodynamics (cQED) with modular quasi-lumped models (2021). Preprint at http://arxiv.org/abs/2103.10344.

26. Minev, Z. K. et al. Qiskit metal: an open-source framework for quantum device design & analysis (Q-EDA). zenodo. https://zenodo.org/record/4618154 (2021).

27. Barends, R. et al. Coherent Josephson qubit suitable for scalable quantum integrated circuits. *Phys. Rev. Lett.* **111**, 080502 (2013).

28. Minev, Z. et al. Planar multilayer circuit quantum electrodynamics. *Phys. Rev. Appl.* **5**, 044021 (2016).

29. Brecht, T. et al. Multilayer microwave integrated quantum circuits for scalable quantum computing. *npj Quantum Inf.* **2**, 16002 (2016).

30. Reagor, M. et al. Quantum memory with millisecond coherence in circuit QED. *Phys. Rev. B* **94**, 014506 (2016).

31. Gambetta, J. M., Chow, J. M. & Steffen, M. Building logical qubits in a superconducting quantum computing system. *npj Quantum Inf.* **3**, 2 (2017).

32. Rosenberg, D. et al. 3D integrated superconducting qubits. *npj Quantum Inf.* **3**, 42 (2017).

33. Versluis, R. et al. Scalable quantum circuit and control for a superconducting surface code. *Phys. Rev. Appl.* **8**, 034021 (2017).

34. Naik, R. K. et al. Random access quantum information processors using multimode circuit quantum electrodynamics. *Nat. Commun.* **8**, 1904 (2017).

35. Krantz, P. et al. A quantum engineer's guide to superconducting qubits. *Appl. Phys. Rev.* **6**, 021318 (2019).

36. Kjaergaard, M. et al. Superconducting qubits: current state of play. *Annu. Rev. Condens. Matter Phys.* **11**, 369–395 (2020).

37. Josephson, B. Possible new effects in superconductive tunnelling. *Phys. Lett.* **1**, 251–253 (1962).

38. Vijay, R., Levenson-Falk, E. M., Slichter, D. H. & Siddiqi, I. Approaching ideal weak link behavior with three dimensional aluminum nanobridges. *Appl. Phys. Lett.* **96**, 223112 (2010).

39. Kerman, A. J. Metastable superconducting qubit. *Phys. Rev. Lett.* **104**, 027002 (2010).

40. Larsen, T. W. et al. Semiconductor-nanowire-based superconducting qubit. *Phys. Rev. Lett.* **115**, 127001 (2015).

41. De Lange, G. et al. Realization of microwave quantum circuits using hybrid superconducting-semiconducting nanowire Josephson elements. *Phys. Rev. Lett.* **115**, 127002 (2015).

42. Janvier, C. et al. Coherent manipulation of Andreev states in superconducting atomic contacts. *Science* **349**, 1199–1202 (2015).

43. Maleeva, N. et al. Circuit quantum electrodynamics of granular aluminum resonators. *Nat. Commun.* **9**, 3889 (2018).

44. Wang, J. I.-J. et al. Coherent control of a hybrid superconducting circuit made with graphene-based van der Waals heterostructures. *Nat. Nanotechnol.* **14**, 120–125 (2019).

45. Minev, Z. K. *Catching and Reversing a Quantum Jump Mid-Flight*. PhD thesis, Yale Univ. (2019).

46. Koch, J. et al. Charge-insensitive qubit design derived from the Cooper pair box. *Phys. Rev. A* **76**, 42319 (2007).

47. Yurke, B. & Denker, J. S. Quantum network theory. *Phys. Rev. A* **29**, 1419–1437 (1984).

48. Devoret, M. H. Quantum fluctuations in electrical circuits. In *A UJFG/a NATO Advanced Study Institute Conference* (eds Reynaud, S., Giacobino, E. & Zinn-Justin, J.) 351 (Elsevier, 1997).

49. Gloos, K., Poikolainen, R. S. & Pekola, J. P. Wide-range thermometer based on the temperature-dependent conductance of planar tunnel junctions. *Appl. Phys. Lett.* **77**, 2915 (2000).

50. Blais, A., Huang, R.-S., Wallraff, A., Girvin, S. M. & Schoelkopf, R. J. Cavity quantum electrodynamics for superconducting electrical circuits: an architecture for quantum computation. *Phys. Rev. A* **69**, 062320 (2004).

51. Wallraff, A. et al. Strong coupling of a single photon to a superconducting qubit using circuit quantum electrodynamics. *Nature* **431**, 162–167 (2004).

52. Yan, F. et al. The flux qubit revisited to enhance coherence and reproducibility. *Nat. Commun.* **7**, 12964 (2016).

53. Minev, Z., Pop, I. M. & Devoret, M. H. Planar superconducting whispering gallery mode resonators. *Appl. Phys. Lett.* **103**, 142604 (2013).

54. Paik, H. et al. Observation of high coherence in Josephson junction qubits measured in a three-dimensional circuit QED architecture. *Phys. Rev. Lett.* **107**, 240501 (2011).

55. Rigetti, C. et al. Superconducting qubit in a waveguide cavity with a coherence time approaching 0.1 ms. *Phys. Rev. B* **86**, 100506 (2012).

56. Axline, C. et al. An architecture for integrating planar and 3D cQED devices. *Appl. Phys. Lett.* **109**, 42601 (2016).

57. Yurke, B. & Buks, E. Performance of cavity-parametric amplifiers, employing kerr nonlinearites, in the presence of two-photon loss. *J. Lightwave Technol.* **24**, 5054–5066 (2006).

58. Ho Eom, B., Day, P. K., LeDuc, H. G. & Zmuidzinas, J. A wideband, low-noise superconducting amplifier with high dynamic range. *Nat. Phys.* **8**, 623–627 (2012).

59. Vissers, M. R. et al. Frequency-tunable superconducting resonators via nonlinear kinetic inductance. *Appl. Phys. Lett.* **107**, 062601 (2015).

60. Mortensen, H. L., Mølmer, K. & Andersen, C. K. Normal modes of a superconducting transmission-line resonator with embedded lumped element circuit components. *Phys. Rev. A* **94**, 053817 (2016).

61. Koops, M. C., van Duyneveldt, G. V. & de Bruyn Ouboter, R. Direct observation of the current-phase relation of an adjustable superconducting point contact. *Phys. Rev. Lett.* **77**, 2542–2545 (1996).

62. Bretheau, L., Girit, Ç., Ö., Pothier, H., Esteve, D. & Urbina, C. Exciting Andreev pairs in a superconducting atomic contact. *Nature* **499**, 312–315 (2013).

63. Peltonen, J. T. et al. Coherent dynamics and decoherence in a superconducting weak link. *Phys. Rev. B* **94**, 180508 (2016).

64. Mooij, J. E. & Nazarov, Y. V. Superconducting nanowires as quantum phase-slip junctions. *Nat. Phys.* **2**, 169–172 (2006).

65. Ku, J., Manucharyan, V. & Bezryadin, A. Superconducting nanowires as nonlinear inductive elements for qubits. *Phys. Rev. B* **82**, 134518 (2010).

66. Abay, S. et al. Charge transport in InAs nanowire Josephson junctions. *Phys. Rev. B* **89**, 214508 (2014).

67. Casparis, L. et al. Gatemon benchmarking and two-qubit operations. *Phys. Rev. Lett.* **116**, 150505 (2016).

68. Shim, Y.-P. & Tahan, C. Bottom-up superconducting and Josephson junction devices inside a group-IV semiconductor. *Nat. Commun.* **5**, 4225 (2014).

69. Zimmerman, J. E. & Silver, A. H. Macroscopic quantum interference effects through superconducting point contacts. *Phys. Rev.* **141**, 367–375 (1966).

70. Clarke, J. & Braginski, A. I. (eds *The SQUID Handbook* (Wiley-VCH Verlag, 2004).

71. Frattini, N. E. et al. 3-Wave mixing Josephson dipole element. *Appl. Phys. Lett.* **110**, 222603 (2017).

72. Manucharyan, V. E. et al. Evidence for coherent quantum phase slips across a Josephson junction array. *Phys. Rev. B* **85**, 024521 (2012).

73. Pop, I. M. et al. Coherent suppression of electromagnetic dissipation due to superconducting quasiparticles. *Nature* **508**, 369–372 (2014).

74. Muppalla, P. R. et al. Bistability in a mesoscopic Josephson junction array resonator. *Phys. Rev. B* **97**, 024518 (2018).

75. 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).

76. Hutter, C., Tholén, E. A., Stannigel, K., Lidmar, J. & Haviland, D. B. Josephson junction transmission lines as tunable artificial crystals. *Phys. Rev. B* **83**, 014511 (2011).

77. Masluk, N. A., Pop, I. M., Kamal, A., Minev, Z. K. & Devoret, M. H. Microwave characterization of Josephson junction arrays: implementing a low loss superinductance. *Phys. Rev. Lett.* **109**, 137002 (2012).

78. Bell, M. T., Sadovskyy, I. A., Ioffe, L. B., Kitaev, A. Y. & Gershenson, M. E. Quantum superinductor with tunable nonlinearity. *Phys. Rev. Lett.* **109**, 137003 (2012).

79. Weißl, T. et al. Kerr coefficients of plasma resonances in Josephson junction chains. *Phys. Rev. B* **92**, 104508 (2015).

80. Macklin, C. et al. A near-quantum-limited Josephson traveling-wave parametric amplifier. *Science* **350**, 307–310 (2015).

81. Martinis, J. M. et al. Decoherence in Josephson qubits from dielectric loss. *Phys. Rev. Lett.* **95**, 210503 (2005).

82. Patel, U. et al. Coherent Josephson phase qubit with a single crystal silicon capacitor. *Appl. Phys. Lett.* **102**, 012602 (2013).

83. Dial, O. et al. Bulk and surface loss in superconducting transmon qubits. *Supercond. Sci. Technol.* **29**, 044001 (2016).
84. Vissers, M. R., Weides, M. P., Kline, J. S., Sandberg, M. & Pappas, D. P. Identifying capacitive and inductive loss in lumped element superconducting hybrid titanium nitride/aluminum resonators. *Appl. Phys. Lett.* **101**, 022601 (2012).

85. Wenner, J. et al. Surface loss simulations of superconducting coplanar waveguide resonators. *Appl. Phys. Lett.* **99**, 113513 (2011).

86. Geerlings, K. et al. Improving the quality factor of microwave compact resonators by optimizing their geometrical parameters. *Appl. Phys. Lett.* **100**, 192601 (2012).

87. Sandberg, M. et al. Long-lived, radiation-suppressed superconducting quantum bit in a planar geometry. *Appl. Phys. Lett.* **102**, 072601 (2012).

88. Wang, C. et al. Surface participation and dielectric loss in superconducting qubits. *Appl. Phys. Lett.* **107**, 162601 (2015).

89. Bruno, A. et al. Reducing intrinsic loss in superconducting resonators by surface treatment and deep etching of silicon substrates. *Appl. Phys. Lett.* **106**, 182601 (2015).

90. Brecht, T. et al. Demonstration of superconducting micromachined cavities. *Appl. Phys. Lett.* **107**, 192603 (2015).

91. Gao, J. *The Physics of Superconducting Microwave Resonators*. PhD thesis, California Institute of Technology (2008).

92. Geerlings, K. L. *Improving Coherence of Superconducting Qubits and Resonators*. PhD thesis, Yale Univ. (2013).

93. Brecht, T. et al. Micromachined integrated quantum circuit containing a superconducting qubit. *Phys. Rev. Appl.* **7**, 044018 (2017).

94. Martinis, J. M. & Megrant, A. UCSB final report for the CSQ program: review of decoherence and materials physics for superconducting qubits. Preprint at http://arxiv.org/abs/1410.5793 (2014).

95. Minev, Z. K. et al. pyEPR: the energy-participation-ratio (EPR) open-source framework for quantum device design. zenodo https://doi.org/10.5281/zenodo.4744447 (2021).

96. Leghtas, Z. et al. Confining the state of light to a quantum manifold by engineered two-photon loss. *Science* **347**, 853–857 (2015).

97. Mundhada, S. O. et al. Generating higher-order quantum dissipation from lower-order parametric processes. *Quantum Sci. Technol.* **2**, 024005 (2017).

98. Touzard, S. et al. Coherent oscillations inside a quantum manifold stabilized by dissipation. *Phys. Rev. X* **8**, 021005 (2018).

99. Campagne-Ibarcq, P. et al. Deterministic remote entanglement of superconducting circuits through microwave two-photon transitions. *Phys. Rev. Lett.* **120**, 200501 (2018).

100. Wang, Z. et al. Cavity attenuators for superconducting qubits. *Phys. Rev. Appl.* **11**, 014031 (2019).

101. Grimm, A. et al. The Kerr-Cat qubit: stabilization, readout, and gates. https://www.nature.com/articles/s41586-020-2587-z (2019).

102. Minev, Z. K. et al. To catch and reverse a quantum jump mid-flight. *Nature* **570**, 200–204 (2019).

103. Campagne-Ibarcq, P. et al. A stabilized logical quantum bit encoded in grid states of a superconducting cavity. https://www.nature.com/articles/s41586-020-2603-3 (2019).

104. Winkel, P. et al. Nondegenerate parametric amplifiers based on dispersion-engineered Josephson-junction arrays. *Phys. Rev. Appl.* **13**, 024015 (2020).

105. Winkel, P. et al. Implementation of a transmon qubit using superconducting granular aluminum. *Phys. Rev. X* **10**, 031032 (2020).

106. Rigetti, C. T. *Quantum Gates for Superconducting Qubits*. PhD thesis, Yale Univ. (2009).

107. Lecocq, F. et al. Junction fabrication by shadow evaporation without a suspended bridge. *Nanotechnology* **22**, 315302 (2011).

108. Pop, I. M. et al. Fabrication of stable and reproducible submicron tunnel junctions. *J. Vacuum Sci. Technol. B* **30**, 010607 (2012).

109. Ambegaokar, V. & Baratoff, A. Tunneling between superconductors. *Phys. Rev. Lett.* **11**, 104–104 (1963).

110. Bergeal, N. et al. Phase-preserving amplification near the quantum limit with a Josephson ring modulator. *Nature* **465**, 64–68 (2010).

111. Abdo, B., Kamal, A. & Devoret, M. Nondegenerate three-wave mixing with the Josephson ring modulator. *Phys. Rev. B* **87**, 014508 (2013).

112. Roy, A. & Devoret, M. H. Introduction to parametric amplification of quantum signals with Josephson circuits. *Comptes Rendus Phys.* **17**, 740–755 (2016).

113. Reagor, M. J. *Superconducting Cavities for Circuit Quantum Electrodynamics*. PhD thesis, Yale Univ. (2016).

114. Gambetta, J. et al. Qubit-photon interactions in a cavity: measurement-induced dephasing and number splitting. *Phys. Rev. A* **74**, 042318 (2006).

115. Gambetta, J. et al. Quantum trajectory approach to circuit QED: quantum jumps and the Zeno effect. *Phys. Rev. A* **77**, 012112 (2008).

116. Gambetta, J. M., Houck, A. A. & Blais, A. Superconducting qubit with Purcell protection and tunable coupling. *Phys. Rev. Lett.* **106**, 030502 (2011).

117. Srinivasan, S. J., Hoffman, A. J., Gambetta, J. M. & Houck, A. A. Tunable coupling in circuit quantum electrodynamics using a superconducting charge qubit with a V-shaped energy level diagram. *Phys. Rev. Lett.* **106**, 083601 (2011).

118. Diniz, I., Dumur, E., Buisson, O. & Auffèves, A. Ultrafast quantum nondemolition measurements based on a diamond-shaped artificial atom. *Phys. Rev. A* **87**, 033837 (2013).

119. Dumur, E. et al. V-shaped superconducting artificial atom based on two inductively coupled transmons. *Phys. Rev. B* **92**, 020515 (2015).

120. Zhang, G., Liu, Y., Raftery, J. J. & Houck, A. A. Suppression of photon shot noise dephasing in a tunable coupling superconducting qubit. *npj Quantum Inf.* **3**, 1 (2017).

121. Roy, T. et al. Implementation of pairwise longitudinal coupling in a three-qubit superconducting circuit. *Phys. Rev. Appl.* **7**, 054025 (2017).

122. Devoret, M., Girvin, S. & Schoelkopf, R. Circuit-QED: how strong can the coupling between a Josephson junction atom and a transmission line resonator be? *Ann. Phys.* **16**, 767–779 (2007).

123. Houck, A. A. et al. Controlling the spontaneous emission of a superconducting transmon qubit. *Phys. Rev. Lett.* **101**, 080502 (2008).

## ACKNOWLEDGEMENTS

We thank S.M. Girvin, R.J. Schoelkopf, S. Nigg, H. Paik, A. Blais, H.E. Türeci, F. Solgun, V. Sivak, S. Touzard, N. Frattini, S. Shankar, C. Axline, V.V. Albert, K. Chou, A. Petrescu, D. Cody, A. Eickbusch, and E. Flurin for valuable discussions and the I. Siddiqi and B. Huard groups for using the early versions of pyEPR. This research was supported by the US Army Research Office (ARO) Grant No. W911NF-18-1-0212. Z.M. acknowledges partial support from the ARO (W911NF-16-1-0349). M.H.D. acknowledges partial support from the ARO (W911NF-18-1-0020) and the Air Force Office of Scientific Research (FA9550-19-1-0399). The view and conclusions contained in this document are those of the authors and should not be interpreted as representing the official policies, either expressed or implied, of the Army Research Office or the US Government. The US Government is authorized to reproduce and distribute reprints for Government purposes notwithstanding any copyright notation herein.

## AUTHOR CONTRIBUTIONS

Z.K.M. conceived and developed the method and the theory, performed the experiments, and analyzed the data. Z.L. and M.H.D. contributed to the theory. I.M.P. and S.O.M. contributed to the measurement and fabrication of the devices L.C. assisted with the simulations. Z.K.M. and M.H.D. wrote the manuscript. All authors provided suggestions, discussed the results, and contributed to the manuscript.

## COMPETING INTERESTS

The authors declare no competing interests.

## ADDITIONAL INFORMATION

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41534-021-00461-8.

**Correspondence** and requests for materials should be addressed to Z.K.M.

**Reprints and permission information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_11.png" /></figure>

**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/.

© The Author(s) 2021