Latest Papers

The simulation of large-scale quantum systems is one of the most sought-after applications of quantum computers. Of particular interest for near-term demonstrations of quantum computational advantage are analog quantum simulations, which employ analog controls instead of digitized gates. Most analog quantum simulations to date, however, have been performed using qubit-based processors, despite the fact that many physical systems are more naturally represented in terms of qudits (i.e., $d$-level systems). Motivated by this, we present an experimental realization of the Lipkin-Meshkov-Glick (LMG) model using an analog simulator based on a single superconducting transmon qudit with up to $d = 9$ levels. This is accomplished by moving to a rotated frame in which evolution under any time-dependent local field and one-axis twisting can be realized by the application of multiple simultaneous drives. Combining this analog drive scheme with universal control and single-shot readout of the qudit state, we provide a detailed study of five finite-size precursors of quantum criticality in the LMG model: dynamical phase transitions, closing of the energy gap, Kibble-Zurek-like dynamics, statistics of the order parameter, and excited-state phase transitions. For each experiment we devise a protocol for extracting the relevant properties which does not require any prior knowledge of the system eigenstates, and can therefore be readily extended to higher dimensions or more complicated models. Our results cement high-dimensional transmon qudits as an exciting path towards simulating many-body physics.

Scaling quantum computers to large sizes requires the implementation of many parallel qubit readouts. Here we present an ultrastable superconducting-qubit readout method using the multi-tone self-phase-referenced Pound-Drever-Hall (PDH) technique, originally developed for use with optical cavities. In this work, we benchmark PDH readout of a single transmon qubit, using room-temperature heterodyne detection of all tones to reconstruct the PDH signal. We demonstrate that PDH qubit readout is insensitive to microwave phase drift, displaying $0.73^\circ$ phase stability over 2 hours, and capable of single-shot readout in the presence of phase errors exceeding the phase shift induced by the qubit state. We show that the PDH sideband tones do not cause unwanted measurement-induced state transitions for a transmon qubit, leading to a potential signal enhancement of at least $14$~dB.

Motivated by the importance of dispersive readout in quantum technology, we study a prototypical dispersive readout setup that is probed by a squeezed vacuum in a time-reversal-symmetric fashion. To this end, we develop a full-counting-statistics framework for dispersive readout and analyze its measurement information, accompanied by a generalized mean-field approach suitable to deal with non-unitary dynamics. Distinct from conventional input-output theory, our full-counting-statistics approach enables the direct calculation of arbitrary-order cumulants for the measured cumulative (i.e., time-integrated) photonic distribution while maintaining applicability to nonlinear systems. The corresponding Fisher information exhibits an exponential dependence on the squeezing parameter and a robustness against residual nonlinearity, which can even approach the quantum Fisher information, setting an upper limit. This work introduces a conceptually streamlined and computationally efficient framework for continuous quantum measurements, making it well suited for widespread adoption in quantum technologies.

Recent studies find that even drives far below the superconducting gap frequency may cause drive-induced quasiparticle generation (QPG) across Josephson junctions (JJs) of superconducting qubits (SCQs), posing a serious concern for fault-tolerant superconducting quantum computing (FTSQC). Nonetheless, quantitative experimental estimation on QPG rates has remained vague. Here, we investigate QPG using strongly driven SCQs, reaching qubit drive amplitudes up to $2π\times$300 GHz by applying intense drive fields through the readout resonators. The resonator nonlinear responses enable quantification of the energy loss at SCQs, including the contribution from QPG. Surprisingly, the estimated total energy loss rates are far lower than those expected by the Floquet-Markov formalism with QPG as the sole loss mechanism. Meanwhile, calculations that incorporate high-frequency cutoffs (HFCs) in the QPG conductance at approximately 17-20 GHz effectively explain the experimental observations. These results suggest limitations in either the QPG conductance model or the Markovian treatment of the QPG processes. Both possibilities possess crucial implications for handling QPG problems toward FTSQC and for a more deeper understanding of Josephson junctions.

State-of-the-art single-qubit gates on superconducting transmon qubits can achieve the fidelities required for error-corrected computations. However, parameter fluctuations due to qubit instabilities, environmental changes, and control inaccuracies make it difficult to maintain this performance. To mitigate the effects of these parameter variations, we numerically derive gates robust to amplitude and frequency errors using gradient ascent pulse engineering (GRAPE). We analyze how fluctuations in qubit frequency, drive amplitude, and coherence affect gate performance over time. The robust pulses suppress coherent errors from drive amplitude drifts over 15 times more than a Gaussian pulse with derivative removal by adiabatic gate (DRAG) corrections. Furthermore, the robust gates, originally designed to compensate for quasi-static errors, also demonstrate resilience to stochastic, time-dependent noise, which is reflected in the dephasing time. They suppress added errors during increases in dephasing by up to 1.7 times more than DRAG.

We propose a quantum-enhanced picostrain sensor that achieves Heisenberg-limited strain sensing using superconducting qubits. A strain-sensitive qubit s Hamiltonian is coupled to the momentum quadrature of a microwave resonator, transducing mechanical strain $ε$ into amplified spatial displacements of the resonator s phase space. Using homodyne detection of the resonator field and multipartite entanglement of N qubits, the protocol achieves a strain sensitivity $Δε\sim pε$ (picostrain), two orders of magnitude better than classical sensors. The scheme integrates natively with superconducting processors, enabling in-situ diagnostic and nanoscale material characterization.

Quantum error correction (QEC) requires ancilla qubits to extract error syndromes from data qubits which store quantum information. However, ancilla errors can propagate back to the data qubits, introducing additional errors and limiting fault-tolerance. In superconducting quantum circuits, Kerr-cat qubits (KCQs), which exhibit strongly biased noise, have been proposed as ancillas to suppress this back-action and enhance QEC performance. Here, we experimentally demonstrate a beamsplitter interaction between a KCQ and a transmon, realizing an effective $\hat{Z}_{cat}\hat{X}_q$ coupling that can be employed for parity measurements in QEC protocols. We characterize the interaction across a range of cat sizes and drive amplitudes, confirming the expected scaling of the interaction rate. These results establish a step towards hybrid architectures that combine transmons as data qubits with noise-biased bosonic ancillas, enabling hardware-efficient syndrome extraction and advancing the development of fault-tolerant quantum processors.

Introducing new components and functionalities into quantum devices is critical in advancing state-of-the-art hardware. Here, we propose superconducting diodes (SDs) as a coherent nonreciprocal element in circuit quantum electrodynamics (cQED) architectures. In particular, we use an asymmetric SQUID as an SD controlled with a flux bias - nonreciprocal element with single control handle and on-chip modality. We spectroscopically characterize SD and show that flux bias acts cooperatively with the nonlinear diode response to induce direction-dependent resonance shifts in the transmission spectrum. We show that even with modest diode efficiency the isolation isolation ratio is sufficiently high, and scales with multiple SDs. We demonstrate the use of the SD as a coupler to realize coherent nonreciprocal qubit-qubit coupling. With a minimal two qubit system, we demonstrate nonreciprocal half-iSWAP, thereby showcasing the potential of intrinsic nonreciprocity as a tool to perform arbitrary two-qubit gates. Our work enables high-fidelity signal routing and entanglement generation in all-to-all connected microwave quantum networks, where nonreciprocity is embedded at the device level.

We present an emulator-based and hardware feasibility study of Aurora-DD, a phase-coherence compensation method that integrates a sign-based feedback update of a global phase offset (Delta phi) with a fixed-depth XY8 dynamical decoupling (DD) scaffold. The feedback optimization is performed offline on a calibrated emulator and the resulting Delta phi* is deployed as pre-calibrated phase compensation on hardware. This represents an "offline closed-loop, online open-loop" feasibility demonstration. Using an Aer-based emulator calibrated with ibm_fez device parameters, Aurora-DD achieves substantial reductions in mean-squared error of the measured expectation value <Z>, yielding 68-97% improvement across phase settings phi = 0.05, 0.10, 0.15, 0.20 over n=30 randomized trials. These large-n emulator results provide statistically stable evidence that the combined effect of XY8 and Delta phi* suppresses both dephasing and systematic phase bias. On real superconducting hardware (ibm_fez), we perform a small-sample (n=3) multi-phase validation campaign. Aurora-DD yields point estimates corresponding to approximately 99.2-99.6% reduction in absolute error relative to a no-DD baseline across all tested phase points. These hardware numbers are reported transparently as feasibility evidence under tight queue and credit constraints. In contrast, the auxiliary Aurora+ZNE branch exhibits instability: shallow two-point ZNE occasionally amplifies calibration inconsistencies and produces large error outliers. We therefore relegate ZNE analysis to the Appendix and position Aurora-DD (without ZNE) as the primary contribution. Overall, the combined results support pre-calibrated Aurora-DD as a practical, stable, and hardware-compatible phase-coherence compensator for NISQ devices in single-qubit settings.

Numerical optimal control (GRAPE) can in principle discover pulse shapes that suppress all coherent gate error to machine precision. But when does that capability actually matter? We present a systematic comparison of Gaussian, DRAG, and GRAPE pulses for single-qubit gates on a three-level transmon model parameterized by IQM Garnet hardware ($T_1 = 37\,μ$s, $T_2 = 9.6\,μ$s,$α/2π= -200$ MHz), with the explicit goal of identifying the regimes where numerical optimization provides genuine practical advantage over analytical methods. Our central finding is that properly calibrated DRAG already operates near the decoherence floor. At 20 ns gate time, GRAPE eliminates all coherent error ($1 - F < 10^{-15}$), but DRAG achieves $1 - F = 4.9 \times 10^{-4}$ in coherent error alone,and $8.4 \times 10^{-4}$ under full decoherence -- only $1.2\times$ above GRAPE's decoherence-limited performance. More surprisingly,DRAG is \emph{more robust} than GRAPE to qubit frequency detuning (minimum fidelity 0.990 vs.\ 0.931 over $\pm 5$ MHz), the dominant calibration uncertainty in charge-noise-limited transmons. GRAPE retains superior amplitude robustness (minimum fidelity 0.994 vs.\ 0.990) and provides the only route to guaranteed zero coherent error, which matters at short gate times ($\lesssim 15$ ns) where perturbative corrections break down. These results lead to concrete calibration guidance: (1) properly calibrated DRAG is sufficient for gate times $\gtrsim 20$ ns on hardware with $T_2/T \gtrsim 500$, (2) GRAPE is necessary at short gate times or when targeting error rates well below the decoherence floor, and (3) robust optimal control incorporating frequency uncertainty should be used when detuning is the dominant noise source. We decompose the full error budget (coherent, $T_1$, $T_2$, control noise) and provide the open-source QubitPulseOpt framework for reproducing all results.

The development of large-scale quantum processors benefits from superconducting qubits that can operate at elevated temperatures and be fabricated with scalable, foundry-compatible processes. Atomic layer deposition (ALD) is increasingly being adopted as an industrial standard for thin-film growth, particularly in applications requiring precise control over layer thickness and composition. Here, we report superconducting qubits based on NbN/AlN/NbN trilayers deposited entirely by ALD. By varying the number of ALD cycles used to form the AlN barrier, we achieve Josephson tunneling through barriers of different thicknesses, with critical current density spanning seven orders of magnitude, demonstrating the uniformity and versatility of the process. Owing to the high critical temperature of NbN, transmon qubits based on these all-nitride trilayers exhibit microsecond-scale relaxation times, even at temperatures above 300 mK. These results establish ALD as a viable low-temperature deposition technique for superconducting quantum circuits and position all-nitride ALD qubits as a promising platform for operation at elevated temperatures.

Superconducting transmon qubits are a promising platform for quantum computation, yet they face significant fidelity degradation due to connectivity noise, particularly in the intermediate coupling regime where noise levels are substantial. While prior works suggest that high fidelity requires operating in regimes with strongly suppressed noise, maintaining such conditions under practical experimental constraints remains challenging. To address this, we investigate quantum gate operations in fully connected transmon rings, examining both SWAP and general circuits. Our study reveals that fidelity can be significantly enhanced by tuning gate operation durations, with local maxima emerging even under strong noise conditions. These fidelity enhancements occur consistently across different qubit numbers and operation types, and for specific initial states -- particularly those with favorable symmetry or entanglement properties -- the achieved fidelities approach quantum error correction thresholds. Furthermore, we develop a supervised machine learning model that accurately predicts the optimal operation durations for new devices, enabling efficient optimization without extensive experimental simulations. These results provide a pathway toward robust quantum circuit design in noisy experimental environments.

The coherence of superconducting quantum computers is severely limited by material defects that create parasitic two-level-systems (TLS). Progress is complicated by lacking understanding how TLS are created and in which parts of a qubit circuit they are most detrimental. Here, we present a method to determine the individual positions of TLS at the surface of a transmon qubit. We employ a set of on-chip gate electrodes near the qubit to generate local DC electric fields that are used to tune the TLS' resonance frequencies. The TLS position is inferred from the strengths at which TLS couple to different electrodes and comparing them to electric field simulations. We found that the majority of detectable surface-TLS was residing on the leads of the qubit's Josephson junction, despite the dominant contribution of its coplanar capacitor to electric field energy and surface area. This indicates that the TLS density is significantly enhanced near shadow-evaporated electrodes fabricated by lift-off techniques. Our method is useful to identify critical circuit regions where TLS contribute most to decoherence, and can guide improvements in qubit design and fabrication methods.

We investigate the decay dynamics of a three-level artificial atom, a superconducting transmon qubit, weakly coupled to a continuum of modes in a broadband, one-dimensional cavity. Using the resolvent formalism, we derive analytical expressions for the resonance frequency shifts and widths, which are then evaluated numerically for a Gaussian density of states. We identify two distinct dynamical regimes, differentiated by the ratio of the qubit's coupling strength to the continuum bandwidth. When this ratio is much less than one, the system exhibits a Markovian regime in which the resonance width is practically independent of energy within the continuum band. As the ratio increases, the system transitions to a non-Markovian regime where the resonance width becomes strongly energy-dependent. In this regime, the qubit interacts with the continuum faster than the continuum can erase the information from the qubit's past. Furthermore, we demonstrate that the coupling between the transmon's second level and its ground state significantly influences the decay dynamics of the third level. The interaction between these two levels opens a fast two-photon decay channel, which broadens the transmon's second level.

Precise control of superconducting qubits is essential for advancing both quantum simulation and quantum error correction. Recently, transmon qubit systems employing the single-transmon coupler (STC) scheme have demonstrated high-fidelity single- and two-qubit gate operations by dynamically tuning the effective coupling between qubits. However, the integration of STCs increases the number of control lines, thereby posing a significant bottleneck for chip routing and scalability. To address this challenge, we propose a robust control line multiplexing scheme based on a double-transmon coupler (DTC) architecture, which enables shared coupler control lines to substantially reduce wiring complexity. Moreover, we experimentally verify that this multiplexed configuration efficiently suppresses undesirable static $ZZ$ coupling while maintaining accurate control over two-qubit gate operations. We further demonstrate the feasibility of the architecture through two distinct gate implementations: a fast coupler $Z$-control-based CZ gate and a parametric iSWAP gate. To validate the practical applicability of this multiplexing approach in quantum circuits, we prepare Bell and three-qubit GHZ states using the proposed scheme with fidelity exceeding 99% and 96%, respectively. This multiplexed DTC architecture offers significant potential to minimize wiring overhead in two-dimensional qubit arrays, thereby greatly enhancing the scalability of superconducting quantum processors.

Eliminating residual ZZ interactions in a two-qubit system is essential for reducing coherent errors during quantum operations. In a superconducting circuit platform, coupling two transmon qubits via a transmon coupler has been shown to effectively suppress residual ZZ interactions. However, in such systems, perfect cancellation usually requires the qubit-qubit detuning to be smaller than the individual qubit anharmonicities, which exacerbates frequency crowding and microwave crosstalk. To address this limitation, we introduce TFT (Transmon-Fluxonium-Transmon) architecture, wherein two transmon qubits are coupled via a fluxonium qubit. The coupling mediated by the fluxonium eliminates residual ZZ interactions even for transmons detuned larger than their anharmonicities. We experimentally identified zero-ZZ interaction points at qubit-qubit detunings of 409 MHz and 616 MHz from two distinct TFT devices. We then implemented an adiabatic, coupler-flux-biased controlled-Z gate on both devices, achieving CZ gate fidelities of 99.64(6)% and 99.68(8)%.

Entangling gates between neighboring physical qubits are essential for quantum error correction. Implementing them in an all-microwave manner simplifies signal routing and control apparatus of superconducting quantum processors. We propose and experimentally demonstrate an all-microwave controlled-Z (CZ) gate that achieves high fidelity while suppressing residual ZZ interactions. Our approach utilizes a fixed-frequency transmon coupler and multi-path coupling, thereby sufficiently reducing the net transverse interaction between data transmons to suppress residual ZZ interactions. The controlled phase arises from the dispersive frequency shift of the $|gf\rangle$$\unicode{x2013}$$|eg\rangle$ transition between the coupler and one of the data transmons conditioned on the state of the other data transmon. Driving the transitions at the midpoint of two dispersively shifted resonance frequencies induces state-dependent geometric phases to achieve the CZ gate. Crucially, with this scheme, we can maintain a small net transverse interaction between two data transmons while increasing the coupling between the data and coupler transmons to accelerate the CZ-gate speed. Additionally, we measure the coupler state after the gate to detect a subset of decoherence-induced failures that occur during the gate operation. These events constitute erasure errors with known locations, enabling erasure-aware quantum error-correcting codes to improve future logical qubit performance.

Simulating the time-dynamics of quantum many-body systems was the original use of quantum computers proposed by Feynman, motivated by the critical role of quantum interactions between electrons in the properties of materials and molecules. Accurately simulating such systems remains one of the most promising applications of general-purpose digital quantum computers, in which all the parameters of the model can be programmed and any desired physical quantity output. However, performing such simulations on today's quantum computers at a scale beyond the reach of classical methods requires advances in the efficiency of simulation algorithms and error mitigation techniques. Here we demonstrate programmable digital quantum simulation of the dynamics of the 2D Fermi-Hubbard model -- one of the best-known simplified models of electrons in crystalline solids -- at a scale beyond exact classical state-vector simulation. We implement simulations of this model on lattice sizes up to ${6\times 6}$ using 72 qubits on Google's Willow quantum processor, across a range of physical parameters, including different on-site electron-electron interaction strengths and magnetic flux values, and study phenomena including formation of magnetic polarons (charge carriers surrounded by local magnetic polarisation), dynamical symmetry-breaking in stripe-ordered states, attraction of charge carriers on an entangled background state known as a valence bond solid, and the approach to equilibrium through thermalisation. We validate our results against exact calculations in parameter regimes where these are feasible, and compare them to approximate classical simulations performed using tensor network and operator propagation methods. Our results demonstrate that meaningful programmable digital quantum simulation of many-body interacting electron models is now feasible on state-of-the-art quantum hardware.

Unitary and dissipative models of quantum dynamics are linear maps on the space of states or density matrices. This linearity encodes the superposition principle, a key feature of quantum theory. However, this principle can break down in effective non-Hermitian dynamics arising from postselected quantum evolution. We theoretically characterize and experimentally investigate this breakdown in a dissipative superconducting transmon circuit. Within the circuit's three-level manifold, no-jump postselection generates an effective non-Hermitian Hamiltonian governing the excited two-level subspace and an anti-Hermitian nonlinearity. We prepare different initial states and use quantum state tomography to track their evolution under this effective, nonlinear Hamiltonian. By comparing the evolution of a superposition-state to a superposition of individually-evolved basis states, we test linearity and observe clear violations which we quantify across the exceptional-point (EP) degeneracy of the non-Hermitian Hamiltonian. We extend the analysis to density matrices, revealing a breakdown in linearity for the two-level subspace while demonstrating that linearity is preserved in the full three-level system. These results provide direct evidence of nonlinearity in non-Hermitian quantum evolution, highlighting unique features that are absent in classical non-Hermitian systems.

Decoherence of quantum bits arises primarily from the parasitic resistance within the qubit. This study presents the analysis of the decoherence process due to physical interactions between the qubit photons and parasitic resistance atoms, utilizing exclusively the Caldeira-Leggett electrical model, without relying on external Hamiltonians. The analysis shows a good agreement between the model of the electrical noise and the Johnson-Nyquist noise. The emission and absorption rates of the qubit's coherent loss, required for the Lindblad master equation that approximates the decoherence, are obtained. A numerical substitution in the analysis result yields a strong correlation with previous measurements. The present analysis enables also the derivation of the appropriate circuit characteristics for future simulations.