Electromagnetic simulations form an indispensable part of the design and optimization process for superconducting quantum devices. Although several commercial platforms exist, open-source alternatives optimized for high-performance computing remain limited. To address this gap, we introduce SQDMetal, a Python-based API that integrates Qiskit Metal (IBM), Gmsh, Palace (AWS), and Paraview (Kitware) into an open-source, highly parallel simulation workflow for superconducting quantum circuits. SQDMetal enables accurate, efficient, and scalable simulations while remaining community-driven and free from commercial constraints. In this work, we validate SQDMetal through mesh convergence studies which benchmark SQDMetal against COMSOL Multiphysics and Ansys, demonstrating excellent agreement for both eigenmode and electrostatic (capacitance) simulations. Furthermore, we simulate superconducting resonators and transmon qubits, showing reasonable agreement with experimental measurements. SQDMetal also supports advanced capabilities, including Hamiltonian extraction via the energy participation ratio (EPR) method, incorporation of kinetic inductance effects, and full 3D modelling of device geometry for improved predictive accuracy. By unifying open-source tools into a single framework, SQDMetal lowers the barriers to entry for community members seeking to access high-performance simulations to assist in the design and optimization of their devices.
Electromagnetic simulations form an indispensable part of the design and optimization process for superconducting quantum devices. Although several commercial platforms exist, open-source alternatives optimized for high-performance computing remain limited. To address this gap, we introduce SQDMetal, a Python-based API that integrates Qiskit Metal (IBM), Gmsh, Palace (AWS), and Paraview (Kitware) into an open-source, highly parallel simulation workflow for superconducting quantum circuits. SQDMetal enables accurate, efficient, and scalable simulations while remaining community-driven and free from commercial constraints. In this work, we validate SQDMetal through mesh convergence studies which benchmark SQDMetal against COMSOL Multiphysics and Ansys, demonstrating excellent agreement for both eigenmode and electrostatic (capacitance) simulations. Furthermore, we simulate superconducting resonators and transmon qubits, showing reasonable agreement with experimental measurements. SQDMetal also supports advanced capabilities, including Hamiltonian extraction via the energy participation ratio (EPR) method, incorporation of kinetic inductance effects, and full 3D modelling of device geometry for improved predictive accuracy. By unifying open-source tools into a single framework, SQDMetal lowers the barriers to entry for community members seeking to access high-performance simulations to assist in the design and optimization of their devices.
Accurate electromagnetic (EM) feature extraction, including element characterization, eigenmodes, and field distributions, is essential for superconducting quantum circuit design. To streamline this process, we present a workflow built around Palace, an open-source, high-performance finite element method solver tailored for quantum applications. Starting from circuit layouts, the workflow automates mesh generation, multiple EM solver processing, and EM-to-Hamiltonian post-processing. We benchmark the workflow on a chip with bare resonators and qubits coupled with readout resonators, achieving resonator frequencies prediction within 0.3% and 3 out of 4 external couplings within 16% of cryogenic measurements. These results demonstrate open-source EM tools can match commercial accuracy while offering scalable, license-free analysis. Our Palace-based workflow provides an accessible and extensible foundation for rapid superconducting circuit development and materials-loss studies.
Quantum computing is advancing rapidly, with superconducting quantum computing (SQC) emerging as one of the most promising platforms for scalable quantum information processing. In contrast to fabrication-focused approaches prevalent in current literature, this work presents a fully simulation-based design and tuning methodology for a small-scale superconducting quantum chip - addressing practical limitations posed by fabrication resource constraints - and consequently excludes physical chip fabrication or real-time performance evaluation (e.g., fidelity and error rate measurements). A meticulous, step-by-step electromagnetic simulation was conducted on the chip, designed using Qiskit Metal and Ansys Electronics, with qubit parameters tuned and verified using the Lumped Oscillator Model (LOM) and Energy Participation Ratio (EPR) methods. The design achieves good parameter convergence with relative tuning error rates (Δ%) below 0.05% after 15 adaptive passes and approximately 50,000 solved mesh elements, yielding an Energy Participation Ratio (EPR) of 91.8 % for the chip substrate (Silicon). Furthermore, the study includes a theoretical comparison between two fundamental qubit-qubit coupling schemes. A key contribution of this work lies in systematically analyzing how geometric variations in qubit design affect critical quantum parameters such as qubit frequencies, anharmonicities, coupling strengths, and dispersive shifts - thereby offering a comprehen-sive, simulation-driven framework for early-stage quantum chip prototyping.
We present an open-source database of superconducting quantum device designs that may be used as the starting point for customized devices. Each design can be generated programmatically using the open-source Qiskit Metal package, and simulated using finite-element electromagnetic solvers. We present a robust workflow for achieving high accuracy on design simulations. Many designs in the database are experimentally validated, showing excellent agreement between simulated and measured parameters. Our database includes a front-end interface that allows users to generate ``best-guess'' designs based on desired circuit parameters. This project lowers the barrier to entry for research groups seeking to make a new class of devices by providing them a well-characterized starting point from which to refine their designs.
An accurate and efficient numerical electromagnetic model for superconducting qubits is essential for characterizing and minimizing design-dependent dielectric losses. The energy participation ratio (EPR) is the commonly adopted metric used to evaluate these losses, but its calculation presents a severe multiscale computational challenge. Conventional finite element method (FEM) requires 3D volumetric meshing, leading to prohibitive computational costs and memory requirements when attempting to capture singular electric fields at nanometer-thin material interfaces. To address this bottleneck, we propose SesQ, a surface integral equation simulator tailored for the precise simulation of the EPR. By applying discretization on 2D surfaces, deriving a semi-analytical multilayer Green's function, and employing a dedicated non-conformal boundary mesh refinement scheme, SesQ accurately resolves singular edge fields without an explosive growth in the number of unknowns. Validations with analytically solvable models demonstrate that SesQ accelerates capacitance extraction by roughly two orders of magnitude compared to commercial FEM tools. While achieving comparable accuracy for capacitance extraction, SesQ delivers superior precision for EPR calculation. Simulations of practical transmon qubits further reveal that FEM approaches tend to significantly underestimate the EPR. Finally, the high efficiency of SesQ enables rapid iteration in the layout optimization, as demonstrated by minimizing the EPR of the qubit pattern, establishing the simulator as a powerful tool for the automated design of low-loss superconducting quantum circuits.
We propose SPICE-Q, a SPICE-inspired design-technology co-optimization framework for superconducting quantum processors. Rather than replacing tools such as HFSS, Qiskit Metal, pyEPR, SQcircuit, SQuADDS, scqubits, or QuTiP, SPICE-Q aims to connect them through a unified, traceable data chain spanning process rules, layout, electromagnetic simulation, energy-participation-ratio and circuit quantization, Hamiltonian extraction, noise analysis, cryogenic test, and manufacturing feedback. The central mapping is from process and PDK constraints to layout geometry, electromagnetic modes, equivalent circuit parameters, effective Hamiltonians, and finally metrics such as frequency, coupling, anharmonicity, decoherence, readout performance, and yield. This flow must capture Josephson-junction variability, transmon frequency allocation, resonator and Purcell constraints, coupler crosstalk, microwave routing, 3D interconnects, material/interface loss, package modes, and wafer-scale process statistics. By introducing standardized model interfaces, statistical parameter models, model cards, version governance, and closed-loop calibration from cryogenic and fabrication data, SPICE-Q frames superconducting quantum-chip design as an engineering workflow rather than a collection of isolated simulations. We argue that scalable and fault-tolerant quantum processors will require such a continuous model chain from device physics and electromagnetic fields to quantum dynamics, noise, manufacturability, and system-level yield.
Abstract Designing superconducting quantum circuits involves optimizing the layout to achieve certain target parameters. This optimization process usually depends on iterative electromagnetic simulations, which are computationally expensive and require manual intervention to adjust the layout parameters. Here, we present a method to efficiently automate the optimization of superconducting circuits, which significantly reduces the need for manual intervention. The method’s efficiency arises from approximate nonlinear model-driven (ANMod) parameter updates, which are constructed from the user’s physical knowledge. Additionally, we provide a full implementation using the ANMod-method as an open-source Python package, QDesignOptimizer. The package automates the design workflow by combining high-accuracy electromagnetic simulations in ansys HFSS and energy participation ratio (pyEPR) analysis integrated with the design tool quantum-metal (formerly known as Qiskit-Metal). Our implementation supports modular and flexible subsystem-level analysis and is easily extensible to optimize for additional parameters. The ANMod-method is not specific to superconducting circuits; as such, it can be applied to a range of nonlinear optimization problems across science and technology.