---
title: 'QuLTRA: Hybrid Lumped & Transmission Analyzer'
url: https://www.emergentmind.com/topics/qultra
type: topic
---

# QuLTRA: Hybrid Lumped & Transmission Analyzer

Searching arXiv for QuLTRA and closely related superconducting-circuit analysis tools.
First, locating the QuLTRA paper itself.
QuLTRA, short for **Quantum hybrid Lumped and TRansmission lines circuits Analyzer**, is an **open-source Python package** for simulating superconducting quantum circuits containing **both lumped elements** and **distributed microwave structures** such as **coplanar waveguide (CPW) transmission lines** and **multi-line CPW couplers**. Its defining feature is the direct treatment of distributed components, without discretizing them into lumped-element equivalents, combined with Hamiltonian extraction via the **energy participation ratio (EPR)** method. In the formulation presented in "QuLTRA: Quantum hybrid Lumped and TRansmission lines circuits Analyzer" [2509.03651], this enables the extraction of mode frequencies, anharmonicities, cross-Kerr interactions, and Purcell decay rates while naturally including higher-order modes of distributed structures.

## 1. Scope and design objective

QuLTRA addresses a practical bottleneck in superconducting circuit design: the need to map a target Hamiltonian onto a realizable physical circuit when that circuit contains both nonlinear lumped elements and distributed microwave components. The design loop is described as iterative: one guesses component values, extracts the Hamiltonian, compares it to the target, and revises the design accordingly [2509.03651].

Within that workflow, existing tools are divided into two broad categories. **Full electromagnetic solvers** such as **Ansys HFSS + pyEPR** are described as very accurate, but computationally expensive and difficult to use iteratively for large or complex layouts. **Lumped-element-only Hamiltonian tools** such as **QuCAT, SQcircuit, and CircuitQ** are described as faster, but unable to naturally represent distributed components like CPW resonators, feedlines, and couplers unless these are discretized into many LC sections [2509.03651].

The central purpose of QuLTRA is to bridge that gap. It retains a circuit-based modeling workflow while handling distributed microwave structures in a physically direct manner. This suggests that QuLTRA is positioned as an intermediate layer between full-wave electromagnetic verification and reduced lumped-network quantization, especially for circuit-QED architectures in which CPW resonators, couplers, and feedlines are structurally central.

## 2. Quantization framework and Hamiltonian extraction

QuLTRA uses the **energy participation ratio** method to extract effective Hamiltonian parameters from the linearized circuit [2509.03651]. The Hamiltonian is written as a linear multimode part plus Josephson nonlinearity:
$$
\hat{H}=\sum_{m=1}^M{\hbar \omega_m \hat{a}_m^\dagger \hat{a}_m} -\sum_i N_i^2 E_i \left[\cos(\hat{\phi}_i/N_i)+\hat{\phi}_i^2/(2N_i^2)\right],
$$
with
$$
E_i=\left(\frac{\hbar}{2e}\right)^2\frac{1}{L_i},
$$
where $\hat{\phi}_i$ is the reduced flux across the $i$-th Josephson array.

The reduced flux is expanded in normal modes as
$$
\phi_i = \sum_{m=1}^M \phi_{m,i}(a_m^\dagger + a_m),
$$
and the mode-dependent flux zero-point fluctuation $\phi_{m,i}$ is related to the energy participation ratio $p_{m,i}$ by
$$
\phi_{m,i}^2 = p_{m,i}\frac{\hbar \omega_m}{2E_i}.
$$
The participation ratio itself is defined as
$$
p_{m,i}= \frac{\text{linear inductive energy stored in the $i$-th array in mode $m$}}{\text{total inductive energy stored in mode $m$}}.
$$

After expanding the cosine to fourth order and applying the rotating-wave approximation, the effective Hamiltonian becomes
$$
\hat{H} \simeq \sum_{m=1}^M
\left[ \hbar(\omega_{m}- \Delta_{m}) \hat{a}_{m}^\dagger \hat{a}_{m} - \frac{1}{2} \hbar \alpha_{m} \hat{a}_{m}^{\dagger 2} \hat{a}_{m}^2 \right] - \sum_{n \ne m} \frac{1}{2} \hbar \chi_{mn} \hat{a}_{m}^\dagger \hat{a}_{m} \hat{a}_{n}^\dagger \hat{a}_{n},
$$
where $\omega_m$ is the bare mode frequency, $\Delta_m$ the Lamb shift, $\alpha_m$ the self-Kerr or anharmonicity, and $\chi_{mn}$ the cross-Kerr coupling [2509.03651]. The relations
$$
\Delta_m=\frac{1}{2}\sum_n \chi_{mn}, \qquad \alpha_m=\frac{\chi_{mm}}{2}
$$
and
$$
\chi_{mn} = \sum_i \frac{1}{\hbar}\frac{E_i}{N_i^2}\phi_{m,i}^2\phi_{n,i}^2 = \sum_i \frac{\hbar\,\omega_m\omega_n}{4E_i} \frac{p_{m,i}p_{n,i}}{N_i^2}
$$
show that once the linear mode frequencies $\omega_m$ and the EPRs $p_{m,i}$ are known, the nonlinear Hamiltonian parameters follow directly [2509.03651].

This organization is consequential. It separates the problem into a linear eigenmode computation and a nonlinear parameter extraction stage, which is the same broad philosophy associated with EPR-based quantization. A plausible implication is that the method is especially attractive when the linear microwave environment is structurally complicated but the nonlinearity remains localized in Josephson elements.

## 3. Linear-mode computation through admittance matrices

QuLTRA represents the circuit as a network of multi-port subcircuits described by **admittance matrices** [2509.03651]. After selecting a reference node, Kirchhoff’s current law yields
$$
\mathbf{Y}(z)\mathbf{V}=0,
$$
and the mode frequencies are obtained from the roots of
$$
\det(\mathbf{Y}(z))=0.
$$

For a lossless circuit, $z=j\omega$. For a lossy circuit, the paper uses
$$
z=\kappa/2 + j\omega,
$$
where $\omega$ is the mode frequency and $\kappa$ is the mode linewidth or dissipation rate [2509.03651].

The root-finding procedure is split by circuit class. For **lossless circuits**, QuLTRA scans the target frequency interval in overlapping sub-intervals and uses **Brent’s method** to find roots of the real or imaginary part of the determinant, depending on matrix-size parity. Candidate roots are filtered by checking that $|\det(\mathbf{Y}(z_0))| < 10^{-3}$ [2509.03651]. For **lossy circuits**, it first finds the lossless roots by shorting resistors and then refines each root with **Newton–Raphson** using the lossy circuit as the target. This yields both resonant frequencies and linewidths.

After obtaining a root $z_m$, QuLTRA computes the corresponding node-voltage eigenvector $\mathbf{V}_m$ from the null space of $\mathbf{Y}(z_m)$. Those node voltages are then used to compute inductive energies and therefore participation ratios:
$$
p_{m,i}=\frac{E_{L,i}^{(m)}}{E_L^{(m)}}.
$$
For lumped inductors and Josephson junctions, the inductive energy uses
$$
E_L=\frac{1}{2}L|I|^2, \qquad I = YV.
$$
For distributed CPW lines and couplers, the current distribution is computed analytically and the inductive energy is integrated over the structure [2509.03651].

This admittance-matrix formalism is one of the package’s general design principles. The paper states that the method is not fundamentally limited to CPW lines; it can be generalized to any component described by an admittance matrix [2509.03651].

## 4. Direct modeling of CPW transmission lines and multi-line couplers

A central contribution of QuLTRA is the direct treatment of **straight, lossless CPW lines** as genuine transmission-line elements rather than as LC ladders [2509.03651]. For a line, the voltage and current profiles are written as
$$
V(x)=V^+e^{-jkx}+V^-e^{jkx}, \qquad I(x)=I^+e^{-jkx}+I^-e^{jkx},
$$
with
$$
k=\omega/v, \qquad v=\frac{1}{\sqrt{LC}}=\frac{c}{\sqrt{\varepsilon_{\mathrm{eff}}}}, \qquad Z_0=\sqrt{\frac{L}{C}},
$$
and
$$
V^\pm = \pm Z_0 I^\pm.
$$

A line of length $\ell$ is represented by the two-port admittance matrix
$$
\mathbf{Y}_\mathrm{tml}(z)= \frac{1}{Z_0 \sinh(z\ell /v)} \begin{pmatrix} \cosh(z\ell /v) & -1 \\ -1 & \cosh(z\ell /v) \end{pmatrix}.
$$
Given the endpoint voltages for a mode, QuLTRA reconstructs $V^\pm$, then $I^\pm$, then the full current profile $I(x)$, and computes the line inductive energy analytically:
$$
E_\mathrm{tml} = \frac{L}{2}\int_0^\ell |I(x)|^2\,dx.
$$
Because the CPW is treated as a continuum, QuLTRA automatically captures the **fundamental mode**, **higher-order harmonics**, and the correct mode structure of distributed resonators [2509.03651].

QuLTRA also directly models **multi-line CPW couplers**, including two-line and three-line structures. For an $n$-conductor geometry,
$$
\mathbf{V}(x)=\mathbf{V}^+e^{-jkx}+\mathbf{V}^-e^{jkx}, \qquad \mathbf{I}(x)=\mathbf{I}^+e^{-jkx}+\mathbf{I}^-e^{jkx},
$$
with
$$
\mathbf{V}^\pm = \pm \mathbf{Z}\mathbf{I}^\pm,
$$
where
$$
\mathbf{Z}=\sqrt{\mathbf{L}\mathbf{C}^{-1}},
$$
and the matrices satisfy
$$
\mathbf{L}\mathbf{C}=\frac{1}{v^2}\mathbf{I}_n.
$$

The capacitance matrix is obtained via **conformal mapping / Schwarz–Christoffel transformations**. The paper gives the mapping
$$
w(z)=\int_0^z \frac{ \prod_{j\in\{0,\dots,n\}\setminus\{m,m+1\}}(z'-c_j) }{ \prod_{k=0}^{n}(z'-a_k)^{1/2} \prod_{l=0}^{n}(z'-b_l)^{1/2} } \,dz',
$$
with the $c_j$ chosen so that grounded conductors map to a continuous boundary, and then obtains the capacitance matrix elements from the mapped geometry as
$$
C_{im}=(\varepsilon_r+1)\varepsilon_0 \frac{\operatorname{Re}(b'_i)-\operatorname{Re}(a'_{i+1})}{\operatorname{Im}(b'_m)}.
$$
For a two-line or three-line coupler of length $\ell$, QuLTRA constructs a **4-port admittance matrix** $\mathbf{Y}_\mathrm{cpl}(z)$ by exciting one port at a time with 1 V while grounding the others, then calculating the resulting port currents from the propagated-wave solution. The stored inductive energy is
$$
E_\mathrm{cpl}=\frac{1}{2}\int_0^\ell \mathbf{I}(x)^\dagger \mathbf{L}\mathbf{I}(x)\,dx.
$$

The significance of this construction is methodological rather than merely implementational. It removes the need to approximate distributed structures with a large number of lumped sections and therefore avoids the associated growth in node count and the convergence issues of LC discretization. This suggests that the package is especially advantageous when higher distributed modes, coupler geometry, or line-mediated environmental effects are central to the device physics.

## 5. Higher-order modes, accuracy, and validation

QuLTRA treats distributed components as objects with **infinitely many modes**, and includes as many as needed by finding additional roots of $\det(\mathbf{Y}(z))=0$ [2509.03651]. The paper emphasizes that this is especially important for $\lambda/4$ and $\lambda/2$ resonators, ultra-strong coupling systems, Purcell filters using higher resonator poles, and multiplexed readout circuits where harmonics matter.

The validation program compares QuLTRA against **Ansys HFSS + pyEPR**, **QuCAT**, and **literature benchmarks** [2509.03651]. The reported examples are summarized below.

| Validation case | Reported comparison | Reported outcome |
|---|---|---|
| Transmon + $\lambda/2$ resonator | HFSS vs QuLTRA | Largest discrepancies below about 8% |
| $\lambda/4$ resonator inductively coupled to a feedline | HFSS vs QuLTRA | Very good linewidth agreement; QuLTRA a few seconds on a laptop, HFSS almost one hour on a server |
| Lumped-element test cases | QuCAT vs QuLTRA | Almost perfect agreement |
| Notch Purcell filter | Literature/HFSS-inspired benchmark | $Z_{12}$, notch behavior, resonance peaks, and Purcell suppression reproduced |

For a transmon capacitively coupled to a $\lambda/2$ readout resonator, the reported values are: qubit mode frequency **HFSS 5.76 GHz vs QuLTRA 5.85 GHz**, resonator mode frequency **9.36 GHz vs 9.34 GHz**, qubit anharmonicity **348 MHz vs 321 MHz**, and cross-Kerr **1.8 MHz vs 1.93 MHz** [2509.03651]. The largest discrepancies are stated to be below about **8%**, mainly in anharmonicity and cross-Kerr, with the likely explanation that QuLTRA uses a capacitance-matrix model for the qubit rather than a full electromagnetic simulation.

For a $\lambda/4$ resonator inductively coupled to a feedline, QuLTRA is reported to match HFSS very well for resonator linewidth versus coupler length, while requiring only **a few seconds on a laptop**, compared with **almost one hour on a server** for HFSS [2509.03651]. For lumped-element test cases, the comparison with QuCAT is described as **almost perfect agreement**.

A literature-inspired notch Purcell filter benchmark is used to compare impedance and dissipation-related quantities. At the optimal resonator length, the reported values are a qubit dissipation rate of $\kappa_q/2\pi = 0.0166$ Hz, a Purcell decay time of **9.59 s**, a readout mode linewidth of $\kappa_r/2\pi = 3.14$ MHz, and a dispersive shift of $\chi_{qr}/2\pi = 2.75$ MHz [2509.03651].

Taken together, these comparisons support two narrower conclusions stated in the source material: first, QuLTRA preserves agreement with lumped-circuit tools on problems those tools already cover; second, it extends direct simulation to distributed structures while avoiding the computational cost of full electromagnetic solvers.

## 6. Application domains in circuit QED

The paper discusses several application classes that illustrate the package’s intended operational regime [2509.03651].

### Purcell filters

QuLTRA is used to model a notch Purcell filter formed by two $\lambda/4$ CPW resonators coupled through a multi-line CPW structure. In that setting, it computes impedance spectra, predicts the placement of impedance notches, and evaluates qubit Purcell decay rates directly from the full circuit. The example also shows extreme sensitivity of Purcell suppression to qubit-frequency detuning from the notch.

A second Purcell-filter example exploits higher-frequency poles of a $\lambda/4$ resonator as an intrinsic notch filter. By tuning the resonator length to align a pole with the qubit frequency, QuLTRA yields a qubit dissipation rate of $\kappa_q/2\pi = 0.0166$ Hz and a Purcell time of **9.59 s**, while maintaining the reported readout linewidth and dispersive shift [2509.03651].

### Multimode ultra-strong coupling

In a transmon coupled to a $\lambda/4$ resonator in the **multimode ultra-strong coupling regime**, QuCAT is described as requiring a Foster ladder with at least **10 LC sections** for accurate qubit frequency and anharmonicity [2509.03651]. QuLTRA instead treats the resonator as a true CPW line and directly obtains a qubit frequency of **8.02 GHz** and an anharmonicity of **352 MHz**. It also computes the Lamb shift as higher resonator modes are included, showing the expected gradual saturation.

### Multiplexed qubit readout

The package is also applied to a shared CPW readout line connected to multiple qubit-resonator subsystems. The analysis demonstrates coupling-position-dependent linewidths, standing-wave effects on a non-terminated line, and maxima in linewidth at positions corresponding to integer multiples of $\lambda_r/2$ for capacitive coupling and integer multiples of $\lambda_r/4$ for inductive coupling [2509.03651]. The paper presents this for two structures: **QSTR1**, comprising qubit + readout resonator + Purcell filter with capacitive coupling, and **QSTR2**, comprising qubit + readout resonator with inductive coupling through a CPW coupler.

These examples indicate the intended niche of QuLTRA: architectures in which distributed microwave engineering and Hamiltonian-level quantum modeling cannot be cleanly separated. A plausible implication is that the package is most useful when one needs both layout-aware microwave structure and rapid iteration at the circuit-model level.

## 7. Advantages, assumptions, and present limitations

The paper enumerates several advantages of QuLTRA [2509.03651]. It models CPW transmission lines and couplers directly instead of approximating them as many lumped elements; it supports fast iterative design relative to full electromagnetic simulation; it extracts mode frequencies, anharmonicities, cross-Kerr couplings, Lamb shifts, and Purcell decay rates through EPR; it includes higher-order resonator modes naturally; and it is built on a general multi-port admittance-matrix formalism.

Those advantages, however, are paired with explicit constraints. The current implementation supports only **ideal**, **lossless**, **straight** CPW transmission lines and straight multi-line couplers [2509.03651]. The paper also notes that bends or curvature may introduce small deviations, that the substrate permittivity is currently hard-coded for silicon in the interface, and that future work may extend the method to more general distributed components and stronger nonlinearity.

These limitations are important for interpretation. QuLTRA is not presented as a replacement for final full-wave validation in arbitrary geometries. Rather, it occupies a specific methodological role in the superconducting-circuit design stack: a circuit-QED analysis tool for hybrid lumped-distributed networks that is substantially more direct than lumped discretization and substantially lighter than full electromagnetic simulation [2509.03651]. In that sense, its principal contribution is to make distributed microwave structures first-class objects in Hamiltonian extraction workflows for superconducting quantum hardware.

Source: https://www.emergentmind.com/topics/qultra