We present a semi-classical method for determining the effective low-energy quantum Hamiltonian of weakly anharmonic superconducting circuits containing mesoscopic Josephson junctions coupled to electromagnetic environments made of an arbitrary combination of distributed and lumped elements. A convenient basis, capturing the multi-mode physics, is given by the quantized eigenmodes of the linearized circuit and is fully determined by a classical linear response function. The method is used to calculate numerically the low-energy spectrum of a 3D-transmon system, and quantitative agreement with measurements is found.
We present a semi-classical method for determining the effective low-energy quantum Hamiltonian of weakly anharmonic superconducting circuits containing mesoscopic Josephson junctions coupled to electromagnetic environments made of an arbitrary combination of distributed and lumped elements. A convenient basis, capturing the multi-mode physics, is given by the quantized eigenmodes of the linearized circuit and is fully determined by a classical linear response function. The method is used to calculate numerically the low-energy spectrum of a 3D-transmon system, and quantitative agreement with measurements is found.
We present a semiclassical method for determining the effective low-energy quantum Hamiltonian of weakly anharmonic superconducting circuits containing mesoscopic Josephson junctions coupled to electromagnetic environments made of an arbitrary combination of distributed and lumped elements. A convenient basis, capturing the multimode physics, is given by the quantized eigenmodes of the linearized circuit and is fully determined by a classical linear response function. The method is used to calculate numerically the low-energy spectrum of a 3D transmon system, and quantitative agreement with measurements is found.
We present a method for calculating the energy levels of superconducting circuits that contain highly anharmonic, inductively shunted modes with arbitrarily strong coupling. Our method starts by calculating the normal modes of the linearized circuit and proceeds with numerical diagonalization in this basis. As an example, we analyze the Hamiltonian of a fluxonium qubit inductively coupled to a readout resonator. While elementary, this simple example is nontrivial because it cannot be efficiently treated by the method known as ``black-box quantization,'' numerical diagonalization in the bare harmonic oscillator basis, or perturbation theory. Calculated spectra are compared to measured spectroscopy data, demonstrating excellent quantitative agreement between theory and experiment.
Abstract 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.
Building on the first-order circuit quantization method [arXiv:2304.12252, arXiv:2401.09120], we provide simple formulas to construct exact Hamiltonians for Josephson-junction-based superconducting qudits capacitively, inductively, or galvanically coupled to passive linear environments. These environments may be multiport, multimode, discrete or continuous, reciprocal or nonreciprocal, and are characterized directly by their impedance or admittance matrices. In the weak-coupling regime, we further derive \emph{divergence-free} dispersive Hamiltonians for mode-resolved environments and transition-resolved weak-coupling master equations for dissipative continua. Mode structure, frequency renormalizations, environment-mediated interactions, decay rates, and directional cross couplings then follow from the same causal immittance response, while spurious Lamb-shift divergences arising from uncontrolled approximations in previous treatments are made explicit and avoided. We apply the theory to a set of illustrative circuits comprising a discrete resonator filter, finite-band metamaterial environments, nonreciprocal waveguide-QED systems, and superconducting giant atoms, for which analytical response matrices can be obtained, although the method is particularly well suited to numerical responses from electromagnetic solvers or experimental characterization. We thereby extend the black-box quantization framework to multiport, dissipative, and nonreciprocal settings, establishing a simple and scalable route toward optimized and automated electromagnetic design of large-scale superconducting quantum hardware.
We propose a new quantization method for superconducting electronic circuits involving a Josephson-junction device coupled to a linear microwave environment. The method is based on an exact impedance synthesis of the microwave environment considered as a blackbox with impedance function Z(s). The synthesized circuit captures dissipative dynamics of the system with resistors coupled to the reactive part of the circuit in a nontrivial way. We quantize the circuit and compute relaxation rates following previous formalisms for lumped element circuit quantization. Up to the errors in the fit our method gives an exact description of the system and its losses.
Superconducting circuits are one of the leading candidates for scalable quantum computation, owing to their long coherence time and compatibility with fabrication techniques. Despite their experimental success, the conventional theoretical framework for superconducting circuit quantisation relies on a semiclassical procedure that requantises the phenomenological Ginzburg-Landau order parameter, introducing the postulated canonical relations and leaving the microscopic origins of circuit parameters hidden. In the present thesis, we develop a rigorous and first-principles approach for the quantisation of superconducting circuits, directly grounded in Bardeen-Cooper-Schrieffer (BCS) theory for superconductivity, and extend this framework to investigate the dynamics of the superconducting Higgs amplitude mode. The central methodological contribution of this work is a projective approach, in which the full many-body electronic Hamiltonian of a superconducting system is projected onto a carefully constructed low-energy BCS Hilbert space, HBCS, spanned by phase parameterised BCS ground states. We systematically characterise the kinematic properties of this subspace, quantify the orthogonality and completeness of the phase basis via the BCS overlap and the projection error, and show that the effective dimensionality of the space is scaled as the square root of the number of single electronic modes. This approach does not impose a spontaneously broken symmetry so that it consistently describes quantised circuits that support superpositions of phases, and the canonical commutation relations between phase and number are derived from the underlying fermionic commutation properties, rather than imposed. This establishes the fully quantum mechanical foundation upon which the standard circuit quantisation approach implicitly rests. The projective formalism is subsequently applied to derive the effective Hamiltonian for superconducting circuit elements. For point-like elements, \eg, capacitive islands and Josephson junctions, projection of the fermionic Hamiltonian onto HBCS yields analytical expressions for the charging energy and the Josephson energy in terms of microscopic material parameters. The resulting Hamiltonian reproduces the standard Cooper-pair-box and transmon Hamiltonians, but now on an unambiguous microscopic footing. We further extend the approach to spatially distributed elements. Inductance is derived from the constitutive response to spatially varying phase torsion in a superconducting wire, yielding the lumped-element inductive energy from first principles. This analysis also identifies the emergent low-energy Hilbert space of the inductor. Combining the results, a complete microscopic derivation of the standard superconducting inductance-capacitance-junction (LCJ) circuit Hamiltonian is recovered. The final part of this thesis extends the projective methodology to the superconducting Higgs mode associated to the amplitude dynamics of the superconducting order parameter. A general BCS space parameterised by phase and gap amplitude is constructed, onto which we project the microscopic Hamiltonian. This yields a generalised eigenvalue problem governing the dynamics of the gap amplitude. The Higgs mode of a superconducting island is thereby shown to behave as a quantum harmonic oscillator. We derive the analytical expressions for the superconducting Higgs mass, spring constant, and oscillation frequency of the gap dynamics. While our results reveal a logarithmic correction to the Higgs frequency reported in previous literature, we numerically validate these results, and compute quantum anharmonic corrections to the Higgs frequency for higher excitations. This thesis proposes a novel, formal and quantitative circuit quantisation theory for superconducting circuits, based on which we explore the Higgs amplitude mode dynamics. The approach is capable of resolving the fundamental question about the compactness of superconducting phase and is extendable to cover quasiparticle contributions. This microscopic circuit theory equips us with a new tool to understand superconducting devices at a level of microscopic detail inaccessible to the standard semiclassical framework, paving the way for a principled treatment of quasiparticle losses and the exploration of exotic superconducting circuit elements.