The superconductor-insulator-superconductor Josephson junction is the fundamental nonlinear element of superconducting circuits. Connecting two junctions in series gives rise to higher-harmonic content in the total energy-phase relation, enabling new design opportunities in multimode circuits. However, the double-junction element hosts an internal mode whose spectrum is set by the finite capacitances of the individual junctions. Using a Born-Oppenheimer approximation that treats the additional mode as fast compared to the qubit mode, we analyze the double-junction circuit element shunted by a large capacitor. Here, we derive an effective single-mode model of the qubit containing a correction term owing to the presence of the internal mode. We explore experimentally relevant parameter regimes and find that our model accurately describes the low-energy spectrum of the qubit. We further discuss how eliminating the internal degree of freedom affects the system's periodic boundary conditions and how this leads to non-uniqueness in performing the Born-Oppenheimer approximation. Finally, we analyze the harmonic content of the double-junction element and discuss its sensitivity to charge noise.

Superconducting transmon qubits based on hybrid superconductor-semiconductor Josephson junctions (gatemons) offer gate tunability, but their relaxation times remain well below those of state-of-the-art transmons, and the origin of this discrepancy is not fully understood. Here, we co-fabricate gatemons and SIS-junction transmons with nominally identical circuit layouts, gate dielectrics, and control lines, so that the Josephson element is the only intentional distinction. Across multiple chips, transmons in this architecture reach relaxation times in the tens of microseconds, whereas gatemons saturate in the few-microsecond range. Using the transmons as on-chip references, we construct a loss budget including Purcell decay, spontaneous emission through the control line, and internal dielectric loss, and find that the corresponding T1 limits exceed all measured gatemon values by more than an order of magnitude. Temperature-dependent T1 measurements follow a common quasiparticle-activation model and yield similar superconducting gaps for S-Sm-S and SIS junctions, indicating that the reduced gatemon coherence is dominated by additional temperature-independent, junction-intrinsic dissipation.

Superconducting quantum circuits derive their nonlinearity from the Josephson energy-phase relation. Besides the fundamental cos ϕ term, this relation can also contain higher Fourier harmonics cos ( k ϕ ) corresponding to correlated tunneling of k Cooper pairs. The parity of the dominant tunneling process, i.e. whether an odd or even number of Cooper pairs tunnel, results in qualitatively different properties, and controlling this opens up a wide range of applications in superconducting technology. However, access to even-dominated regimes has remained challenging and has so far relied on complex multi-junction or all-hybrid architectures. Here, we demonstrate a simple "harmonic parity qubit"; an element that combines two aluminum-oxide tunnel junctions in parallel to a gate-tunable InAs/Al nanowire junction forming a SQUID, and use spectroscopy versus flux to reconstruct its energy-phase relation at 85 gate voltage points. At half flux quantum, the odd harmonics of the Josephson potential can be suppressed by up to two orders of magnitude relative to the even harmonics, producing a double-well potential dominated by even harmonics with minima near  ± π/2. The ability to control harmonic parity enables supercurrent carried by pairs of Cooper pairs and provides a new building block for Fourier engineering in superconducting circuits.

Accurate models of quantum processors are essential for understanding, calibrating, and improving their performance. In practice, model construction must balance physical detail against the experimental and computational effort required to reliably learn parameters. Compact descriptions therefore often rely on assumptions about which interactions, noise processes, or hidden degrees of freedom are relevant. Here we introduce LIMINAL, a data-driven framework for testing such assumptions and selecting minimal adequate Lindblad models. LIMINAL fits nested candidate models to time-resolved tomographic data and uses likelihood-ratio tests to decide when added physical mechanisms are warranted. We apply LIMINAL to a five-qubit superconducting processor, identifying an idling model with three-local Hamiltonian terms and two-local dissipation, while finding no support for three-local dissipation. We further apply it to recover driven single-qubit Hamiltonians, reconstruct a shaped-pulse Hamiltonian without assuming an analytic pulse model, and test hidden-qubit extensions in coupler-mediated dynamics, demonstrating the applicability of the framework for a wide range of tasks.

Interactions that mix conjugate variables, such as the flux through a circuit element and the charge across it, lie outside the reach of the elementary couplings of superconducting circuits. Capacitors connect charge to charge, and inductors connect flux to flux, while no two-terminal element couples flux to charge directly. A native flux-charge coupling would thus serve as a circuit primitive in its own right, opening direct routes to non-reciprocity, protected modes, and unconventional readout. In this work, we demonstrate a flux-charge coupling by harnessing a voltage-tunable Josephson junction with parametrically modulated critical current, which mediates the interaction between a classical charge variable and a quantum flux operator. Relying on parity-selection rules in a hybrid superconducting-semiconductor fluxonium, we isolate the flux-charge coupling from other parasitic capacitive contributions and perform cross-quadrature-activated coherent control of states. Critically, we realize a flux-charge coupling that scales linearly with driving amplitude while keeping the transition energy first-order-insensitive to gate voltage. Such unconventional interaction broadens the toolbox of superconducting circuits with a critical missing component that enables the coherent coupling of conjugate variables.

The fidelity of operations on a solid-state quantum processor is fundamentally bounded by environmental decoherence. Characterizing environmental fluctuations is challenging because the acquisition time of nonadaptive experimental protocols limits temporal precision and can average out rapid features of the underlying dynamics. Here, we overcome this temporal-resolution limit by 2 orders of magnitude using a field-programmable gate-array powered classical controller that adaptively and continuously tracks the relaxation-time fluctuations of two fixed-frequency superconducting transmon qubits, which exhibit average relaxation times of approximately 0.17 ms and occasionally exceed 0.5 ms. We report events in which the relaxation time switches by nearly an order of magnitude over timescales of just tens of milliseconds, rather than minutes or hours as previously reported. Our real-time Bayesian estimation protocol estimates relaxation times within a few milliseconds, close to the decoherence timescale itself. Our statistical analysis further suggests that some of these fast fluctuations arise from two-level systems switching at rates up to 10 Hz, 4 orders of magnitude faster than earlier reports. These results redefine the timescales relevant for calibration in superconducting quantum processing units, establish a reference for rapid relaxation-rate characterization in device screening, and improve our understanding of fast relaxation dynamics.

Superconducting qubit parameters drift on sub-second timescales, motivating calibration and benchmarking techniques that can be executed on millisecond timescales. We demonstrate an on-FPGA workflow that co-locates pulse generation, data acquisition, analysis, and feed-forward, eliminating CPU round trips. Within this workflow, we introduce sparse-sampling and on-FPGA inference tools, including computationally efficient methods for estimation of exponential and sine-like response functions, as well as on-FPGA implementations of Nelder-Mead optimization and golden-section search. These methods enable low-latency primitives for readout calibration, spectroscopy, pulse-amplitude calibration, coherence estimation, and benchmarking. We deploy this toolset to estimate T_1 in 10 ms, optimize readout parameters in 100 ms, optimize pulse amplitudes in 1 ms, and perform Clifford randomized gate benchmarking in 107 ms on a flux-tunable superconducting transmon qubit. Running a closed-loop on-FPGA recalibration protocol continuously for 6 hours enables more than 74,000 consecutive recalibrations and yields gate errors that consistently retain better performance than the baseline initial calibration. Correlation analysis shows that recalibration suppresses coupling of gate error to control-parameter drift while preserving a coherence-linked performance. Finally, we quantify uncertainty versus time-to-decision under our sparse sampling approaches and identify optimal parameter regimes for efficient estimation of qubit and pulse parameters.

Spurious couplings and decoherence degrade the performance of solid-state quantum processors, demanding careful design, calibration, and mitigation protocols. These strategies often rely on characterization of the idling processor, but tomographic recovery of (time-independent) Lindblad dynamics scales exponentially with qubit count. Here, we experimentally benchmark and demonstrate that randomized ("shadow") measurements accelerate Lindblad tomography on a superconducting transmon processor. We first implement extensible Lindblad tomography, which estimates Lindblad parameters using a complete tomographic dataset, and use it as a baseline to benchmark a shadow tomography approach, shadow Lindblad tomography. The shadow approach recycles randomized configurations to estimate the same Lindblad parameters using far fewer resources under physically motivated locality assumptions. We experimentally verify these assumptions in our processor by implementing the protocols on one- and three-qubit subsystems; here, shadow Lindblad tomography reproduces extensible Lindblad tomography within uncertainties while using exponentially fewer configurations. Leveraging this efficiency, we apply shadow Lindblad tomography to the full five-qubit processor and recover all single qubit dissipation and two-qubit coupling parameters in 9 hours of acquisition time compared to an estimated 58 hours for extensible Lindblad tomography. Additionally, our shadow implementation is compatible with conventional Gaussian error propagation, avoiding the use of median-of-means estimators. Together, these results demonstrate how randomized shadow tomography protocols can be practically implemented to learn quantum processor dynamics at an increasing qubit count.

Tunable Josephson harmonics open new avenues for qubit design. We demonstrate a superconducting circuit element consisting of a tunnel junction in series with a superconducting quantum interference device (SQUID) loop, yielding a Josephson potential whose harmonic content is strongly tunable by magnetic flux. Through spectroscopy of the first four qubit transitions, together with an effective single-mode model renormalized by the internal mode, we resolve a second harmonic with an amplitude up to ∼ 10 % of the fundamental. We identify a flux sweet spot where the ground-state dispersive shift vanishes, achieved by balancing the dispersive couplings to the internal and qubit modes. This highly tunable element provides a route toward protected qubits and customizable nonlinear microwave devices.

Tunable Josephson harmonics open up for new qubit design. We demonstrate a superconducting circuit element with a tunnel junction in series with a SQUID loop, yielding a highly magnetic-flux tunable harmonic content of the Josephson potential. We analyze spectroscopy of the first four qubit transitions with a circuit model which includes the internal mode, revealing a second harmonic up to ∼10% of the fundamental harmonic. Interestingly, a sweet spot where the dispersive shift vanishes is achieved by balancing the dispersive couplings to the internal and qubit modes. The highly tunable set-up provides a route toward protected qubits, and customizable nonlinear microwave devices.

The development of quantum circuits based on hybrid superconductor-semiconductor Josephson junctions holds promise for exploring their mesoscopic physics and for building novel superconducting devices. The gate-tunable superconducting transmon qubit (gatemon) is the paradigmatic example of such a superconducting circuit. However, gatemons typically suffer from unstable and hysteretic qubit frequencies with respect to the applied gate voltage and reduced coherence times. Here we develop methods for characterizing these challenges in gatemons and deploy these methods to compare the impact of shunt capacitor designs on gatemon performance. Our results indicate a strong frequency- and design-dependent behavior of the qubit stability, hysteresis, and dephasing times. Moreover, we achieve highly reliable tuning of the qubit frequency with 1 MHz precision over a range of several GHz, along with improved stability in grounded gatemons compared to gatemons with a floating capacitor design.

If a quantum dot is coupled to a topological superconductor via tunneling contacts, each contact hosts a Majorana zero mode in the limit of zero transmission. Close to a resonance and at a finite contact transparency, the resonant level in the quantum dot couples the Majorana modes, but a ground state degeneracy per fermion parity subspace remains if the number of Majorana modes coupled to the dot is five or larger. Upon varying shape-defining gate voltages while remaining close to resonance, a nontrivial evolution within the degenerate ground-state manifold is achieved. We characterize the corresponding non-Abelian holonomy for a quantum dot with chaotic classical dynamics using random matrix theory and discuss measurable signatures of the non-Abelian time-evolution.

Fast, high fidelity control and readout of protected superconducting qubits are fundamentally challenging due to their inherent insensitivity. We propose a flux qubit variation which enjoys a tunable level of protection against relaxation to resolve this outstanding issue. Our qubit design, the double-shunted flux qubit (DSFQ), realizes a generic double-well potential through its three junction ring geometry. One of the junctions is tunable, making it possible to control the barrier height and thus the level of protection. We analyze single- and two-qubit gate operations that rely on lowering the barrier. We show that this is a viable method that results in high fidelity gates as the non-computational states are not occupied during operations. Further, we show how the effective coupling to a readout resonator can be controlled by adjusting the externally applied flux while the DSFQ is protected from decaying into the readout resonator. Finally, we also study a double-loop gradiometric version of the DSFQ which is exponentially insensitive to variations in the global magnetic field, even when the loop areas are non-identical.

Reading out the parity degree of freedom of Majorana bound states is key to demonstrating their nonabelian exchange properties. Here, we present a low-energy model describing localized edge states in a two-arm device. We study parity-to-charge conversion based on coupling the superconductor bound states to a quantum dot whose charge is read out by a sensor. The dynamics of the system, including the readout device, is analyzed in full using a quantum-jump approach. We show how the resulting signal and noise ratio differentiates between local Majorana and Andreev bound states.

Demonstration of Majorana non-Abelian properties is a major challenge in the field of topological superconductivity. In this work, we propose a minimal device and protocol for testing non-Abelian properties using charge-transfer operations between a quantum dot and two Majorana bound states combined with reading the parity state using a second dot. We use an adiabatic perturbation theory to find fast adiabatic paths to perform operations and to account for nonadiabatic errors. We find the ideal parameter sweep and a region in parameter space which reduces the charge-transfer operation time 1-2 orders of magnitude with respect to constant velocity driving. Using realistic parameters, we estimate that the lower bound for the time scale can be reduced to ∼10 ns. Deviations from the ideal parameters lead to the accumulation of an undesired dynamical phase, affecting the outcome of the proposed protocol. We furthermore suggest to reduce the influence from the dynamical phase using a flux echo. The echo protocol is based on the 4π-periodicity of the topological state, absent for trivial bound states.

Demonstration of Majorana nonabelian exchange properties is a major challenge in the field which would show the existence of a new topological quasiparticle. In this work, we propose a minimal device and protocol for testing nonabelian exchange properties using charge-transfer operations between a quantum dot and two Majorana bound states. We use an adiabatic perturbation theory to find fast adiabatic paths to perform operations and to account for non-adiabatic errors. We find the ideal parameter sweep and a region in parameter space where MBS nonabelian properties can be shown. The ideal parameter sweep reduces the charge-transfer operation time in 1-2 orders of magnitude with respect to constant velocity driving. Using realistic parameters, we estimate that the lower bound for the time scale can be reduced to ∼ 10 ns. Deviations from the ideal parameters lead to the accumulation of an undesired dynamical phase, affecting the outcome of the proposed protocol. We reduce the influence from the dynamical phase using an angular echo.

Taro V. Brown, Maarten Golterman, Svend Krøjer, Yigal Shamir, K. Splittorff The Niels Bohr Institute, University of Copenhagen, Blegdamsvej 17, DK-2100, Copenhagen Ø, Denmark Department of Physics and Astronomy, San Francisco State University, San Francisco, CA 94132, USA Raymond and Beverly Sackler School of Physics and Astronomy, Tel Aviv University, 69978 Tel Aviv, Israel (Dated: September 25, 2019)

The ϵ-regime of dilaton chiral perturbation theory is introduced. We compute the dilaton mass, the chiral condensate and the topological susceptibility in the ϵ-regime, as a function of the fermion mass. The microscopic spectral density of the Dirac operator is obtained from dilaton chiral perturbation theory. Our main result is that the chiral condensate and the spectral density are related to their counterparts from ordinary chiral perturbation theory via a simple scaling relation. This relation originates from the mass dependence of the dilaton potential, and is valid in both the ϵ-regime and the p-regime. In the ϵ-regime, moreover, all results agree with the universal predictions to leading order in ϵ.