Back to all posts
Hardware Engineering

Q3D Extractor Output Parameters: RLGC Matrices & Parasitic Extraction for Superconducting Qubits

Silicofeller Engineering TeamPublished by Silicofeller · Technical Reference · June 2026

Q3D Extractor by ANSYS is the industry-standard tool for parasitic extraction in superconducting quantum circuits. It computes the RLGC matrices—Resistance, Inductance, Conductance, and Capacitance—for arbitrary conductor geometries on layered substrates. These matrices are the primary inputs to circuit Hamiltonian models that predict qubit frequency, anharmonicity, coupling strengths, and coherence times.

This reference catalogues 47+ output parameters across 11 categories, with design rules derived from ANSYS Q3D documentation, IBM Quantum, Google AI Quantum, Krantz 2019, and Blais RMP 2021.

"Even small parasitic entries in the RLGC matrices shift qubit frequencies by 10–100 MHz. At 4K, aluminium and niobium become superconducting (R → 0), and kinetic inductance dominates. Normal-metal values from Q3D need temperature-dependent correction before use in the quantum Hamiltonian."

Parameter Category Overview

CategoryParametersKey Focus
1. Resistance Matrix (R) 5Ohmic loss, ground plane, via resistance
2. Inductance Matrix (L) 5Self/mutual inductance, kinetic inductance, Josephson L_J
3. Conductance Matrix (G) 4Substrate leakage, surface conductance
4. Capacitance Matrix (C) 6Qubit C_Σ, coupling caps, pad-to-ground parasitics
5. Parasitic Resistance 4JJ series R, shunt R, wirebond/bump R
6. Parasitic Inductance 4Wirebond L, lead L, slot L, package L
7. Parasitic Capacitance 5Trace-to-ground, substrate, inter-layer, fringe
8. Electromagnetic Coupling 6Q_ext, Q_int, Q_L, Z₀, ε_eff, coupling k²
9. Substrate & Dielectric Loss 5tan δ_bulk, tan δ_MA, tan δ_SA, tan δ_MS, SPR
10. Skin Effect & Frequency-Dependent 7Skin depth, R_ac/R_dc, propagation, RLGC per-unit-length
11. Post-Processing Derived Outputs 7E_C, E_J, g, χ, ζ, α, Γ_P

1. Resistance Matrix (R)

The resistance matrix captures ohmic losses in the conductor network. At room temperature these values are used for fabrication quality assessment; at cryogenic temperatures Al becomes superconducting (R → 0) and the normal-state resistance R_N serves as a junction quality proxy.

IDParameterIdealGood RangeKey Design Note
1DC Self-Resistance R_ii< 1 mΩ 0.5 – 5 mΩAl superconducting at 4K → R→0; use R_N·A as fab quality proxy
2AC Self-Resistance (5–6 GHz)< 5 mΩ (SC at 4K) 5 – 20 mΩSkin effect: R_ac ∝ √f for bulk; &approx; R_dc for thin film < δ_s
3Mutual Resistance R_ij&approx; 0< 50 μΩNon-zero R_ij reveals overlapping ground return paths
4Contact / Via Resistance< 1 mΩ 1 – 5 mΩCritical for multi-chip modules; oxidation is the main failure mode
5Ground Plane Sheet Resistance< 0.1 mΩ/sq 0.1 – 0.5 mΩ/sqPerforated ground planes add ~5–10 pH/sq for vortex pinning

2. Inductance Matrix (L)

The inductance matrix sets the Josephson energy , which directly determines qubit frequency. Both geometric and kinetic inductance contributions must be accounted for in the design model.

IDParameterIdealGood RangeKey Design Note
6Self-Inductance L_ii 1 – 3 nH 0.5 – 5 nHL = L_geo + L_kinetic; target E_J/E_C ~ 50–80
7Mutual Inductance M_ij< 5 pH (idle) 5 – 50 pHIntentional M_ij in flux-tunable couplers; unintentional sets ZZ floor
8Geometric Inductance< 0.5 pH/μm 0.5 – 1 pH/μmSlot cuts in ground plane drastically increase L_geo
9Kinetic Inductance L_k< 2 pH/sq (Al) 1 – 10 pH/sqHigh-L_k materials (NbTiN, TiN) used for KI qubit designs
10Josephson Inductance L_J 8 – 12 nH 5 – 20 nHThe ONLY nonlinear element; L_J/C_Σ ratio sets anharmonicity

3. Conductance Matrix (G)

Conductance matrix entries represent leakage current paths through the substrate and along surfaces. On high-resistivity silicon at 4K, these should be negligible—any measurable conductance indicates substrate quality or surface contamination issues.

IDParameterIdealGood RangeKey Design Note
11Self-Conductance G_ii< 0.1 nS 0.1 – 1 nSG_ii = 1/R_leak; appears as parallel resistance in resonator model
12Substrate Bulk Conductance< 0.01 nS 0.01 – 0.5 nSUse HR-Si (> 10 kΩ·cm) or sapphire; resistivity ↑ 100× at 4K
13Surface / Interface Conductance< 0.001 nS/μm 0.001 – 0.05 nS/μmAdsorbed water and organics increase G_surf; HF dip before cooldown
14Mutual Conductance G_ij< 1 pS 1 – 50 pSNon-zero in presence of surface water; fixed by clean or guard rings

4. Capacitance Matrix (C)

The capacitance matrix is arguably the most important Q3D output for transmon design. The total qubit self-capacitance sets the charging energy , which determines the ratio, anharmonicity, and charge noise sensitivity.

IDParameterIdealGood RangeKey Design Note
15Qubit Self-Capacitance C_Σ 60 – 100 fF 40 – 200 fFC_Σ = Σ|C_ij| from Maxwell matrix; target E_J/E_C = 50–80
16Readout Resonator Cap C_r 200 – 500 fF 100 – 600 fFSets ω_r = 1/√(L_r C_r); target 6.5–8 GHz window
17Qubit–Resonator Coupling C_g 1 – 10 fF 0.5 – 15 fFg/2π target 50–150 MHz; Purcell decay ∝ (g/Δ)² × κ
18Qubit–Qubit Coupling C_J 0.5 – 5 fF 0.2 – 10 fFModern heavy-hex lattice uses tunable couplers to cancel residual ZZ
19Pad-to-Ground Parasitic Cap< 5 fF 1 – 20 fFEach 1 fF of parasitic shifts f_qubit by ~10–30 MHz
20Trace Mutual Capacitance C_ij< 1 fF 1 – 5 fFOverlapping traces on adjacent layers is the primary source

5. Parasitic Resistance

Parasitic resistances in the qubit circuit cause energy dissipation that directly limits T₁. At millikelvin temperatures, the dominant sources are quasiparticle conductance in JJ leads and substrate leakage through lithography residues.

IDParameterIdealGood RangeKey Design Note
21Series JJ Parasitic R< 0.01 Ω 0.01 – 0.1 ΩQuasiparticle poisoning transiently raises R_ser; shielding critical
22Shunt Parasitic R (R_p)> 1 MΩ 100 kΩ – 1 MΩSubstrate residues from lithography most common cause; O₂ plasma clean
23Wirebond / Bump R< 5 mΩ 2 – 20 mΩAu–Au bonds have lower R than Al wedge bonds
24Metal Interface Contact R< 1 mΩ 1 – 10 mΩNative Al₂O₃ must be removed by Ar ion milling before deposition

6. Parasitic Inductance

Parasitic inductance creates impedance discontinuities that shift resonator frequencies, cause reflections, and introduce AC flux errors in bias lines. Flip-chip integration reduces wirebond inductance by 10–20× compared to traditional wire bonding.

IDParameterIdealGood RangeKey Design Note
25Wirebond / Bump L< 1 nH 0.3 – 3 nHFlip-chip In bumps reduce to ~0.1 nH; key for 3D scaling
26Control Line Lead L< 100 pH 50 – 500 pHLong coax adds 1–10 nH; de-embed by calibration
27Ground Plane Slot L< 1 pH/sq 1 – 5 pH/sqVortex-pinning holes add ~5–10 pH/sq but necessary at B > 0
28Package / Board L< 0.5 nH 0.5 – 2 nHSMP connectors (< 0.3 nH) preferred over SMA (~0.5–1 nH) for cryo

7. Parasitic Capacitance

Parasitic capacitances shift qubit frequencies from design targets and create unwanted coupling paths. Fringe capacitance at conductor edges accounts for 30–50% of total coupling capacitance in typical transmon designs.

IDParameterIdealGood RangeKey Design Note
29Trace-to-Ground C 0.1 – 0.4 fF/μm 0.1 – 1 fF/μmSets CPW Z₀ = √(L/C); target 50 Ω
30Pad-to-Substrate C< 10 fF 5 – 30 fFThinning substrate 500→200 μm reduces by ~2.5×
31Inter-Layer Cap (3D)< 5 fF 1 – 20 fFBump height variation (σ ~ 1–2 μm) causes ~0.5–1 fF spread
32Fringe Capacitance 0.02 – 0.1 fF/μm 0.05 – 0.5 fF/μm~30–50% of C_g comes from fringe; underestimating shifts f by 50+ MHz
33Wirebond Pad Parasitic Cap< 50 fF 20 – 150 fFReducing pad 150→80 μm cuts C by ~2.5× with no yield penalty

8. Electromagnetic Coupling

Coupling parameters bridge the gap between Q3D parasitic extraction and the quantum Hamiltonian. External quality factor Q_ext sets readout bandwidth, while Z₀ and ε_eff determine transmission line geometry for target frequencies.

IDParameterIdealGood RangeKey Design Note
34External Quality Factor Q_ext 5×10³ – 2×10⁴ 10³ – 10⁵T₁_Purcell = Q_ext/ω_r × (Δ/g)²; Purcell filter relaxes trade-off
35Internal Quality Factor Q_int> 10⁶ 10⁵ – 10⁶Requires HR-Si or sapphire, clean deposition, minimal surface TLS
36Loaded Quality Factor Q_L 10³ – 10⁴ 500 – 2×10⁴In practice Q_L &approx; Q_ext when Q_int ≫ Q_ext (under-coupled limit)
37CPW Characteristic Impedance Z₀ 50 Ω ± 1 Ω 45 – 55 ΩOn 500 μm Si: 10 μm trace / 6 μm gap → Z₀ &approx; 50 Ω
38Effective Permittivity ε_eff 6.0 – 6.5 5.5 – 7.0ε_eff &approx; (1 + εr)/2 for CPW in air on substrate
39Coupling Coefficient k² 1 – 10 ×10⁻³ 0.5 – 20 ×10⁻³Etch depth variation of 0.1 μm → δk²/k² ~ 5%

9. Substrate & Dielectric Loss

Dielectric loss is the dominant T₁ limiter in planar transmon designs. The loss tangent at each interface (metal-air, substrate-air, metal-substrate) multiplied by the surface participation ratio determines the contribution to qubit decay.

IDParameterIdealGood RangeKey Design Note
40Substrate Bulk Loss tan δ_bulk< 10⁻⁶ 10⁻⁶ – 10⁻⁵tan δ improves by 10–100× on cooling from 300K to 4K
41Metal-Air Interface tan δ_MA< 10⁻³ 10⁻³ – 5×10⁻³Dominant T₁ source in planar transmons; HF vapor clean reduces by 10×
42Substrate-Air Interface tan δ_SA< 5×10⁻⁴ 5×10⁻⁴ – 5×10⁻³H-passivated Si surface (HF dip) shows 5× lower tan δ_SA
43Metal-Substrate Interface tan δ_MS< 5×10⁻³ 5×10⁻³ – 10⁻²Amorphous SiOx layer (1–2 nm) is primary TLS host; HF clean removes it
44Surface Participation Ratio (SPR)< 5 ppm 5 – 50 ppm 1/Q_TLS = Σ pᵢ × tan δᵢ; SPR is the design lever, tan δ is the material lever

10. Skin Effect & Frequency-Dependent Parameters

These parameters characterise the frequency-dependent behaviour of the conductor network. At cryogenic temperatures the skin depth concept is replaced by the London penetration depth λ_L, but room-temperature Q3D values remain essential for pre-cooldown verification.

IDParameterIdealGood RangeKey Design Note
45Skin Depth at 5 GHz (δ_s) 0.5 – 2 μm 0.5 – 3 μmAt 4K SC: replaced by London λ_L &approx; 60–163 nm for thin-film Al
46AC/DC Resistance Ratio&approx; 1.0 (thin film) 1.0 – 2.0Thin-film qubits (100–200 nm) operate below skin-depth limit
47Propagation Constant γα < 0.1 dB/mα 0.1 – 1 dB/mFor interconnects > 10 mm even 0.1 dB/m causes measurable loss
48Phase Velocity v_ph 1.2 – 1.4 ×10⁸ m/s 1.0 – 1.6 ×10⁸ m/sλ/4 at 7 GHz on Si &approx; 4.25 mm
49Per-Unit-Length R'< 0.1 mΩ/mm 0.1 – 2 mΩ/mmAt 4K Al is SC: R' → 0 below T_c
50Per-Unit-Length L' 0.3 – 0.5 nH/mm 0.2 – 0.8 nH/mmL'_kinetic small for Al (~0.01–0.05 nH/mm)
51Per-Unit-Length C' 0.1 – 0.2 pF/mm 0.05 – 0.3 pF/mmCheck: Z₀ = √(L'/C') &approx; 50 Ω as consistency verification

11. Post-Processing Derived Outputs

These parameters are computed from the raw RLGC matrices and represent the quantum Hamiltonian parameters that ultimately determine device performance. The standard design loop is: Q3D → Hamiltonian → Optimise → Iterate.

IDParameterIdealGood RangeKey Design Note
52Charging Energy 200 – 350 MHz 150 – 400 MHz; critical Hamiltonian input
53Josephson Energy 10 – 30 GHz 5 – 50 GHz tunable via flux in split-junction transmons
54Qubit–Resonator Coupling 50 – 150 MHz 20 – 300 MHz ensures dispersive limit;
55Dispersive Shift 1 – 5 MHz 0.5 – 10 MHzSingle-shot readout SNR ; Purcell filter allows larger
56ZZ Coupling Rate < 10 kHz 10 – 50 kHzCentral challenge of transmon scaling; tunable coupler pushes
57Anharmonicity −300 to −150 MHz−350 to −100 MHzGate BW to avoid leakage;
58Purcell Decay Rate < 1 kHz 1 – 10 kHzPurcell filter reduces by 10–100× without affecting readout

Key Takeaways

  • RLGC Matrices from Q3D are the primary inputs to circuit Hamiltonian models. Even small parasitic entries shift qubit frequencies by 10–100 MHz.
  • Superconducting Regime: At 4K, Al and Nb become superconducting → R → 0, L_kinetic dominates. Normal-metal values from Q3D need temperature-dependent correction.
  • Surface Participation: SPR × tan δ controls T₁. Design rule: minimise p_MA below 5 ppm via thick metal, wider CPW gaps, and clean interfaces.
  • Derived Outputs: , , , , , are all computed from Q3D matrices. Iterating Q3D → Hamiltonian → Optimise is the standard qubit design loop.

Sources: ANSYS Q3D · IBM Quantum · Google AI · Krantz 2019 · Blais RMP 2021 · Q3D Analysis — QuantumChipGen Reference Deck

About the Authors

SF

Silicofeller Engineering Team

The Silicofeller team specialises in superconducting quantum chip design automation, electromagnetic simulation, and VLSI-grade layout tooling. Our mission is to make quantum hardware design accessible, reproducible, and physics-grounded.