Energy Participation Ratio Method

• Energy Participation Ratio Method is a universal metric that quantifies the fraction of electromagnetic energy stored in each circuit element of superconducting quantum devices.
• It employs finite-element eigenmode simulation to extract energy integrals, enabling precise Hamiltonian construction and dissipation rate predictions.
• The methodology facilitates design optimization by automating the integration of nonlinear and lossy elements in large-scale, multimode quantum circuits.

The energy participation ratio (EPR) method underpins contemporary modeling, quantization, and design of superconducting quantum circuits, including transmons, fluxonia, and coupled resonator networks. The EPR quantifies how much of a given electromagnetic mode’s energy resides in each discrete nonlinear, lossy, or dissipative element. This universal, geometry-agnostic metric transforms the complex problem of circuit quantization—especially for multimode, multi-junction systems—into a simulation-driven and automatable workflow, linking electromagnetic eigenmode analysis directly to system Hamiltonians and dissipation parameters.

1. Formal Definition and Physical Interpretation

For a linear eigenmode $m$ (frequency $\omega_m$) and device element $j$ (e.g., Josephson junction, lumped inductor, or lossy dielectric region), the EPR $p_{mj}$ is defined as the fraction of the mode's inductive or electric energy stored in that element:

$$p_{mj} = \frac{\text{inductive energy in } j \text{ when mode } m \text{ is excited}}{\text{total energy in mode } m}$$

In quantum terms, for each eigenstate $|\psi_m\rangle$ and inductive element $j$ (with energy $E_j$), this is

$$p_{mj} = \frac{\langle \psi_m | \tfrac{1}{2} E_j \hat{\varphi}_j^2 | \psi_m \rangle}{\langle \psi_m | \tfrac{1}{2} \hat{H}_{\text{lin}} | \psi_m \rangle}$$

where $\hat{\varphi}_j$ is the gauge-invariant phase drop across element $j$ and $\hat{H}_{\text{lin}}$ is the quadratic linearized Hamiltonian (Minev et al., 2020, Yilmaz et al., 2024).

For dielectric losses, the participation ratio $P_i$ is given by integrating the spatially localized energy density:

$$P_i = \frac{\int_{V_i} \epsilon(\mathbf{r}) | \mathbf{E}(\mathbf{r}) |^2 dV}{\int_{\text{all}} \epsilon(\mathbf{r}) | \mathbf{E}(\mathbf{r}) |^2 dV}$$

A high $p_{mj}$ or $P_i$ signifies a large fraction of energy stored in element $j$ or region $i$, which correspondingly exerts a strong influence on nonlinearity or dissipation.

2. Computational Methodology and Eigenmode Simulation

The EPR framework relies on a single electromagnetic finite-element eigenmode simulation of the entire linearized circuit layout, with nonlinear (or lossy) elements represented as parameterized boundary conditions:

• Linearization: Josephson junctions and superinductors are replaced by surface inductance sheets (e.g., $L_j = \Phi_0^2 / E_j$).
• Simulation: The solver (e.g., HFSS) produces mode profiles, frequencies $\omega_m$, and field distributions, from which energy integrals are extracted.
• For quasi-localized interfaces (e.g., metal-dielectric surfaces), a two-step simulation protocol is adopted: coarse-grained global mode simulation and fine-mesh local analyses near edges or junctions, allowing accurate reconstruction of divergent fields and, thus, local participations (Wang et al., 2015).

Table: Key Quantities Extracted in EPR Workflows | Quantity | Description | Method of Extraction | | ---------- | -------------------------------------------------- | ------------------------------------- | | $p_{mj}$ | Mode-$m$ inductive EPR of element $j$ | Sheet current and energy in FE sim. | | $P_i$ | Electric participation in dielectric region $i$ | Field energy density spatial integral | | $s_{mj}$ | EPR sign for relative phase/interference | Line integral of current/field |

3. Hamiltonian Construction and Dissipation Rates

Using all $p_{mj}$, the full circuit Hamiltonian is assembled as:

$H = H_{\text{lin}} + H_{\text{nl}}$

with

$$H_{\text{lin}} = \sum_m \hbar \omega_m a^\dagger_m a_m$$

and the nonlinear part for each element $j$:

$$H_{\text{nl},j} = -E_j \left[ \cos(\varphi_j) + \frac{1}{2} \varphi_j^2 \right]$$

Here, $\varphi_j$ decomposes as $\sum_m \varphi_{mj}(a_m + a^\dagger_m)$, with $$\varphi_{mj} = \sqrt{p_{mj} \hbar \omega_m / 2 E_j}$$ (Yilmaz et al., 2024). For highly anharmonic systems such as fluxonium, the expansion retains the full cosine form and is numerically exponentiated in the truncated Fock basis.

Dissipation rates and device lifetime predictions follow directly:

$$\frac{1}{T_1} = \omega \sum_i P_i\, \tan\delta_i + \Gamma_0$$

where $P_i$ is the surface participation ratio, $\tan\delta_i$ the loss tangent, and $\Gamma_0$ encapsulates non-geometric loss channels (Wang et al., 2015).

4. Algorithmic Extensions: IEPR and Highly Anharmonic Circuits

The inductive-energy participation ratio (IEPR) generalizes EPR for mixed capacitive-inductive networks and explicit normal/bare mode decompositions (Yu et al., 2023). IEPR directly connects energy localization to unitary transformations between bare and normal basis and is systematically compatible with quantum electronic design automation workflows. In IEPR, for normal mode $m$ and circuit element $n$:

$$r_{mn} = \frac{\mathcal{E}_{mn}^{I}}{\mathcal{E}_m^{I}}$$

Hamiltonian extraction, Kerr and cross-Kerr calculations, and coupling constants are thus accessible in a unitary-invariant and computationally minimal procedure.

For very anharmonic circuits (e.g., fluxonium), recent advancements retain the full nonlinear Hamiltonian structure in the EPR formalism. Rather than truncating to quartic order, the complete cosine term is retained, with external flux bias seamlessly incorporated. Numerically, the resulting Hamiltonian matrix is constructed in truncated photon number bases and diagonalized for observable calculations (Yilmaz et al., 2024).

5. Practical Workflow and Design Optimization

Applying EPR techniques proceeds as follows:

1. Draw the full device geometry, including distributed and lumped elements, in an eigenmode solver.
2. Assign boundary conditions that linearize nonlinear elements for simulation.
3. Run eigenmode simulation to extract frequencies, mode profiles, and sub-element energy integrals.
4. Compute all required EPRs ($p_{mj}$ for nonlinear/dissipative elements; $P_i$ for dielectrics).
5. Assemble quantum Hamiltonian, extract Kerr parameters, cross-Kerr couplings, and Lamb shifts.
6. For lossy regions, determine $1/T_1$ contributions via respective $P_i$ and material $\tan\delta$.
7. Iterate device geometry, material choices, or layout to optimize relevant participations and minimize primary limiting losses (Wang et al., 2015, Minev et al., 2020, Yilmaz et al., 2024).

Design heuristics target minimal $P_i$ for deleterious interfaces, achieved via geometric strategies such as increasing electrode area, pad spacing, substrate etching, or using air-bridge architectures. Material cleanliness and in-situ processing are equally critical.

6. Experimental Validation and Comparative Performance

The EPR method and its variants have been validated across a diverse spectrum of superconducting circuits, including 15 transmon qubits, resonator arrays, and fluxonium devices, with frequency predictions matching measured spectra at the $\lesssim 1\%$ level, and nonlinear coupling rates (e.g., cross-Kerr $\chi_{mn}$, anharmonicities) typically agreeing within $5\%$ for large values (Minev et al., 2020, Yilmaz et al., 2024). Dissipation rate predictions align with measured quality factors, provided that dominant loss channels are faithfully represented. For highly nonlinear circuits, EPR with full-cosine treatment yields dispersive shifts $\chi(\Phi_{\text{ext}})$ and transition frequencies that closely track experiment, outperforming lumped-element models in accuracy, particularly outside the weak-anharmonicity regime (Yilmaz et al., 2024).

In comparison to traditional approaches, such as lumped LC network quantization, black-box quantization (BBQ), and frequency-sweep-based normal-mode search, the EPR framework offers:

• Compatibility with distributed multimode field simulations
• Extraction of both bare and dressed (normal) Hamiltonians
• Explicit dissipation budgeting for arbitrary loss channels
• Direct automation within quantum electronic design automation pipelines (Yu et al., 2023)

7. Limitations and Validity Regimes

While EPR-based quantization is universal in scope, several care points apply:

• Finite-element mesh convergence and careful mode normalization are essential for accurate $p_{mj},P_i$ extraction.
• For nanoscopic or spatially inhomogeneous interfaces, statistical absence of resonant two-level systems (TLS) requires judicious cutoffs in participation summation.
• Assumptions include the locality of nonlinear elements (junctions/inductors) and their weak energy participation (justifying perturbative treatment beyond quadratic order, unless the full cosine method is used).
• Extremely strong multimode coupling or extended, non-local nonlinear elements may require alternative theoretical treatments (Wang et al., 2015, Yilmaz et al., 2024).

Notwithstanding these caveats, the energy participation ratio method constitutes the state-of-the-art for first-principles, simulation-driven quantization and dissipation analysis in superconducting quantum circuits, facilitating the systematic and quantitative design of large-scale, fault-tolerant quantum processors.