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.
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.
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.
Superconducting circuits are being used for large-scale quantum devices, and a major challenge is to perform accurate numerical simulations of device parameters. One of the most advanced methods for analyzing superconducting circuit designs is the energy-participation-ratio (EPR) method, which constructs quantum Hamiltonians based on the energy distribution extracted from classical electromagnetic simulations. In the EPR approach, we extract linear terms from finite-element simulations and add nonlinear terms using the energy participation ratio extracted from the classical simulations. However, the EPR method relies on a low-order expansion of nonlinear terms, which is prohibitive for accurately describing highly anharmonic circuits. An example of such a circuit is the fluxonium qubit, which has recently attracted increasing attention due to its high lifetimes and low error rates. In this work, we extend the EPR approach to effectively address highly nonlinear superconducting circuits, and, as a proof of concept, we apply our approach to a fluxonium qubit. Specifically, we design, fabricate, and experimentally measure a fluxonium qubit coupled to a readout resonator. We compare the measured frequencies of both the qubit and the resonator to those extracted from the EPR analysis, and we find an excellent agreement. Furthermore, we compare the dispersive shift as a function of external flux obtained from experiments with our EPR analysis and a simpler lumped-element model. Our findings reveal that the EPR results closely align with the experimental data, providing more accurate estimations compared to the simplified lumped-element simulations.
Fluxonium qubits, designed to mitigate charge noise and enhance anharmonicity, are among the most promising superconducting platforms for quantum computing. To understand and exploit their quantum properties and design novel fluxonium-based architectures with improved functionalities, these systems require an accurate Hamiltonian formulation to capture their energy level structure and quantum dynamics. This work presents a systematic method for constructing the Hamiltonian for multi-loop circuits that partitions the system into a set of uncoupled harmonic oscillators and a coupled anharmonic part originating from the Josephson circuit elements, allowing clear identification of independent modes and isolating the nonlinearity in the Josephson terms. While demonstrated for fluxonium-type multi-loop circuits, this method can be generalized to other superconducting qubit architectures within the broader context of circuit QED, making it a versatile tool for exploring different circuit configurations. Our systematic and flexible modeling approach forms the theoretical basis for the qubit measurement and control experiments validating multi-loop fluxonium architectures.
Superconducting quantum circuits is one of the leading candidates for a universal quantum computer. Designing novel qubit and multiqubit superconducting circuits requires the ability to simulate and analyze the properties of a general circuit. In particular, going outside the transmon approach, we cannot make assumptions on anharmonicity, thus precluding blackbox quantization approaches and necessitating the formal circuit quantization approach. We consider and solve two issues involved in simulating general superconducting circuits. One of the issues is the handling of free modes in the circuit, that is, circuit modes with no potential term in the Hamiltonian. Another issue is circuit size, namely, the challenge of simulating strongly coupled multimode circuits. The main mathematical tool we use to address these issues is the linear canonical transformation in the setting of quantum mechanics. We address the first issue by giving a provably correct algorithm for removing free modes by performing a linear canonical transformation to completely decouple the free modes from other circuit modes. We address the second by giving a series of different linear canonical transformations to reduce intermode couplings, thereby reducing the problem to the weakly coupled case and greatly mitigating the overhead for classical simulation. We benchmark our decoupling methods by applying them to the circuit of two inductively coupled fluxonium qubits, obtaining several orders of magnitude reduction in the size of the Hilbert space that needs to be simulated.
Quantum computers utilize qubits to store and process quantum information. In superconducting quantum computers, qubits are implemented as quantum superconducting resonant circuits. The circuits are operated only at the two energy states, which form the computational basis for the qubit. To suppress leakage to uncomputational states, superconducting qubits are designed to be anharmonic oscillators, which is achieved using one or more Josephson junctions, a nonlinear superconducting element. One of the main challenges in developing quantum computers is minimizing the decoherence caused by environmental noise. Decoherence is characterized by two coherence times, T1 for depolarization processes and T2 for dephasing. This thesis reviews and investigates the decoherence properties of superconducting qubits. The main goal of the thesis is to analyze the tradeoff between anharmonicity and dephasing in a qubit unimon. Recently developed unimon incorporates a single Josephson junction shunted by a linear inductor and a capacitor. Unimon is tunable by external magnetic flux, and at the half flux quantum bias, the Josephson energy is partially canceled by the inductive energy, allowing unimon to have relatively high anharmonicity while remaining fully protected against low-frequency charge noise. In addition, at the sweet spot with respect to the magnetic flux, unimon becomes immune to first-order perturbations in the flux. The sweet spot, however, is relatively narrow, making unimon susceptible to dephasing through the quadratic coupling to the flux noise. In the first chapter of this thesis, we present a comprehensive look into the basic theory of superconducting qubits, starting with two-state quantum systems, followed by superconductivity and superconducting circuit elements, and finally combining these two by introducing circuit quantum electrodynamics (cQED), a framework for building superconducting qubits. We follow with a theoretical discussion of decoherence in two-state quantum systems, described by the Bloch-Redfield formalism. We continue the discussion by estimating decoherence using perturbation theory, with special care put into the dephasing due to the low-frequency 1/f noise. Finally, we review the theoretical model of unimon, which is used in the numerical analysis. As a main result of this thesis, we suggest a design parameter regime for unimon, which gives the best ratio between anharmonicity and T2.