We present a simple algorithm that implements an arbitrary $n$-qubit unitary operator using a Clifford + T circuit with T-count $O(2^{4n/3}{n}^{2/3})$ and ancilla-count $O(2^{2n/3}{n}^{1/3})$. The previous best T-count for unitary synthesis was $O(2^{3n/2}n)$ with $O(2^{n/2})$ ancillae by [Low et al., Quantum 8, 1375 (2024)]. There is a fundamental tradeoff between T gates and ancillae, and our algorithm extends the achievable regime of this tradeoff, reducing the exponent in T-count from $3n/2$ to $4n/3$. The best known lower bound remains $\mathrm{\Omega }(2^n)$, so the optimal T-count for general unitary synthesis is still open. Our construction is based on a recursive application of the cosine-sine decomposition, together with a generalization of the optimal diagonal unitary synthesis method by [D. Gosset et al., Quantum 10, 2168 (2026)] to multi-controlled $k$-qubit unitaries.

Quantum batteries, microscopic devices designed to address energy demands in quantum technologies, promise high power during charging and discharging processes. However, their performance depends crucially on reliability, quantified by the inverse of noise-to-signal ratios (NSRs), i.e., normalized fluctuations of work and power. We establish fundamental limits to this reliability: both work and power NSRs are universally bounded from below by a function of charging speed, imposing a reliability limit inherent to any quantum battery. More strikingly, we find that a quantum mechanical uncertainty relation forbids the simultaneous suppression of work and power fluctuations, revealing a fundamental trade-off that also limits the reliability of quantum batteries. We analyze the trade-off and limits, as well as their scaling behavior, across parallel, collective, and hybrid charging schemes for many-body quantum batteries, finding that increasing power by exploiting stronger entanglement comes at the cost of diminished reliability of power. Similar trends are also observed in the charging of quantum batteries utilizing Ising-like interactions. These suggest that achieving both high power and reliability requires neither parallel nor collective charging, but a hybrid charging scheme with an intermediate range of interactions. Therefore, our analysis shapes the practical and efficient design of reliable and high-performance quantum batteries.

Tunnel-coupled optically active quantum dot molecules (QDMs) have the potential to operate as spin-photon interfaces with coupled spins that interact with two different photon frequencies at the same time. A prerequisite is to deterministically prepare two (electron or hole) spins in the QDM and be able to electrically tune the orbital state couplings. Here, we demonstrate the sequential optical charging of a single QDM with two-electron spins while simultaneously maintaining the ability to widely tune orbital couplings using static electric fields and optically drive the system for quantum light generation. We optically prepare one- and two-spin states, initialize via optical pumping and explore orbital and spin relaxation dynamics for one- and two-spin states as a function of the energy detuning and hybridization of orbital states. For two-spin states, remarkably long singlet-triplet relaxation times are observed, extending beyond $\sim 100\mathrm{\mu }\mathrm{s}$ with strong dependence on the relative energy of ground and excited two-spin states. Qualitative agreement is observed with $\mathit{k}·\mathit{p}$ calculations of phonon-mediated spin relaxation. Our results provide new quantitative understanding of the dynamics of one- and two-spin states and address a number of key prerequisites for using QDMs to create multidimensional photonic cluster states.

Benchmarking physical devices and verifying logical algorithms are important tasks for scalable fault-tolerant quantum computing. Numerous protocols exist for benchmarking devices before running actual algorithms. In this work, we show that both physical and logical errors of fault-tolerant circuits can even be characterized in situ using syndrome data. To achieve this, we map general fault-tolerant Clifford circuits to subsystem codes using the spacetime code formalism and develop a scheme for estimating Pauli noise in Clifford circuits using syndrome data. We give necessary and sufficient conditions for the learnability of physical and logical noise from given syndrome data, and show that we can accurately predict logical fidelities from the same data. Importantly, our approach requires only a polynomial sample size, even when the logical error rate is exponentially suppressed by the code distance, and thus gives an exponential advantage against methods that use only logical data such as direct fidelity estimation. We demonstrate the practical applicability of our methods in various scenarios using synthetic data as well as the experimental data from a recent demonstration of fault-tolerant circuits by Bluvstein et al. [Nature 626, 7997 (2024)]. Our methods provide an efficient, in situ way of characterizing a fault-tolerant quantum computer to help gate calibration, improve decoding accuracy, and verify logical circuits.

Extensive theoretical and experimental work has established high-fidelity electron shuttling in Si/SiGe systems, whereas demonstrations in $\text{Si}/{\text{SiO}}_2$, also known as silicon metal-oxide-semiconductor (SiMOS), devices remain at an early stage. To help address this, we perform 3D simulations of conveyor-belt charge shuttling in a realistic SiMOS device, building on earlier 2D modeling. We solve the Poisson and time-dependent Schrödinger equations for varying shuttling speeds and gate voltages, focusing on potential pitfalls of typical SiMOS devices such as oxide-interface roughness, gate fabrication imperfections, and charge defects along the transport path. The simulations reveal that for low clavier-gate voltages, the additional oxide screening in multi-layer gate architectures causes conveyor-belt shuttling to collapse to the bucket-brigade mode, inducing considerable orbital excitation in the process. Increasing the confinement restores conveyor-belt operation, which we find to be robust against interface roughness, gate misalignment, and charge defects buried in the oxide. However, our results indicate that defects located at the $\text{Si}/{\text{SiO}}_2$ interface can induce significant orbital excitation. For lower conveyor gate biases, positive defects in the transport channel can even capture passing electrons. Hence we identify key challenges and find operating regimes for reliable charge transport in SiMOS architectures.

As fault-tolerant quantum computers scale, certifying the accuracy of computations performed with encoded logical qubits will soon become classically intractable. This creates a critical need for scalable, device-independent certification methods. In this work, we introduce logical accreditation, a framework for efficiently certifying the correctness of quantum computations performed on logical qubits. Our protocol is robust against general noise models, far beyond those typically considered in performance analyses of quantum error-correcting codes. Through numerical simulations, we demonstrate that logical accreditation can scalably certify quantum advantage experiments and indicate the crossover point where encoded computations begin to outperform physical computations. The framework also enables evaluation of whether logical error rates are sufficiently low that error mitigation can be efficiently performed, extends entropy benchmarking to the regime of fault-tolerant computation, and upper bounds the infidelity of the logical output state of a computation. Underlying the framework is a novel randomized compilation scheme that converts arbitrary logical circuit noise into stochastic Pauli noise. This scheme includes a method for twirling non-transversal logical gates beyond the standard $T$ gate, resolving an open problem posed by [Piveteau et al. PRL 127, 200505 (2021)]. By bridging fault-tolerant computation and computational certification, logical accreditation offers a scalable, practical means of certifying the accuracy of quantum computations performed using encoded logical qubits.

Quantum phases at zero temperature can be characterized as equivalence classes under local unitary transformations: two ground states within a gapped phase can be transformed into each other via a local unitary circuit. We generalize this circuit-based characterization of phases to systems at finite-temperature thermal equilibrium described by Gibbs states. We construct a channel circuit that approximately transforms one Gibbs state into another provided the two are connected by a path in parameter space along which a certain correlation-decay condition holds. For finite-dimensional systems of linear size $L$ and approximation error $ϵ$, the locality of the circuit is $\text{Polylog}(\text{Poly}(L)/ϵ)$. The correlation-decay condition, which we specify, is expected to be satisfied in the interior of many noncritical thermal phases, including those displaying discrete symmetry breaking and topological order. As an application, we show that any system in the same thermal phase as a zero-temperature topological code coherently preserves quantum information for a macroscopically long time, establishing self-correction as a universal property of thermal phases. As part of the proof, we provide explicit encoding and decoding channel circuits to encode information into, and decode it from, a system in thermal equilibrium.

Finite-size general security proofs for quantum key distribution based on Rényi entropies have recently been introduced. These approaches are more flexible and provide tighter bounds on the secret key rate than traditional formulations based on the von Neumann entropy. However, deploying them requires minimizing the conditional Rényi entropy, a difficult optimization problem that has hitherto been tackled using ad hoc techniques based on the Frank-Wolfe algorithm, which are unstable and can only handle particular cases. In this work, we introduce a method based on non-symmetric conic optimization for solving this problem. Our technique is fast, reliable, and completely general. We illustrate its performance on several protocols and demonstrate a dramatic runtime improvement over the state of the art.

Moments of ensembles of unitaries play a central role in quantum information theory as they capture the statistical properties of dynamics of systems with some form of randomness. Indeed, concepts such as approximate $t$-designs arise when comparing how close an associated moment operator of a given unitary ensemble is to that of another, reference ensemble. Despite the importance of moment operators, their properties have not been as explored for quantum channels. In this work we develop a theoretical framework to compute moment operators for ensembles of quantum channels, for all moment orders $t$, with a special focus on determining ensembles that can be used as points of reference. By deriving hierarchies between ensembles, via inequalities of their moment operator norms, we give them operational meaning and define useful concepts such as that of channel $t$-designs. Finally, we perform theoretical and numerical studies which show that different types of noise can decrease the norm of the moment operators (e.g., depolarizing noise), as well as increase it (e.g., amplitude damping noise), and generalize noise-induced concentration phenomena to channel-design-induced phenomena. Along the way, we find a block-orthogonal basis for permutations, which greatly simplifies our analyses, and may be of independent interest in moment calculations.

Open quantum systems can be described by unraveling their Lindblad master equation into ensembles of quantum trajectories. Here we investigate how the complexity of such trajectories is affected by conservation laws and other dynamical constraints of the underlying Lindblad evolution. We characterize this complexity using a data-driven approach based on the intrinsic dimension, defined as the minimal number of variables required to encode the information contained in a data set. Applying this framework to several systems, including variants of the quantum top and of the XXZ chain, we find that the intrinsic dimension is sensitive to the structure of their dynamics. The Lindblad evolution in these systems is typically chaotic, but additional constraints arise at specific parameter values, where the dynamics becomes integrable, exhibits Hilbert-space fragmentation, or develops a closed BBGKY hierarchy, leading to pronounced minima in the intrinsic dimension. Our approach results in an unsupervised probe of the complexity of dissipative quantum systems that is sensitive to chaos and ergodicity-breaking phenomena beyond the initial transient regime.

Out-of-equilibrium dynamics of non-integrable Hamiltonian many-body quantum systems are characterized by highly entangled wave functions. Near-maximal entanglement arises in systems exhibiting thermalization or pre-thermalization, where the system converges to a steady state with a fixed energy density. Classical simulation of the time dependence of such wave functions requires exponential resources. However, typical computations aim to estimate expectation values of local operators and correlation functions to some expected precision. For thermalizing systems at sufficiently high-energy densities, such computations can be done without storing the full wave function by instead simulating the evolution of the local operator, which requires significantly fewer resources. Nonetheless, constructing such resource-efficient classical algorithms remains a challenge for intermediate energy densities, where simulating both the wave function and operator evolution is costly. In this paper, we propose a heuristic approach to accelerate the convergence of matrix product state (MPS) simulations of expectation values, applicable across a broad range of energy densities. We estimate the desired observables by rescaling the MPS results at low bond dimensions with a factor that depends on the fidelity of the MPS wave function. Using this technique, we simulated the dynamics of the two-dimensional transverse-field Ising model (TFIM) on a $7\times 8$ grid with periodic boundary conditions, using a maximum bond dimension of $\chi =4096$ on a single A100 GPU, as well as the dynamics of the two-dimensional XY model on grids of size up to $9\times 9$. We compare our TFIM results to similar simulations on a digital quantum processor [R. Haghshenas et al. Nature (London) 653, 56 (2026)], demonstrating excellent agreement and confirming the predictive power of our method.

Random many-body states and matrices are both a useful tool to model certain physical systems and an important asset for quantum computation. Realizing them, however, generally requires an exponential (in system size) amount of resources. Recent research has presented a way out by showing that one can generate random states and matrices, or more precisely a controlled approximation of them, by applying a quantum circuit built in terms of few-body unitary gates. Most of this research has been focussed on the case of quantum circuits composed by fully random unitary gates. Here we consider what happens for circuits that, instead, involve a minimal degree of randomness. Specifically, we concentrate on two different settings: (a) brickwork quantum circuits with a single one-qudit random matrix at a boundary; (b) brickwork quantum circuits with fixed interactions but random one-qudit gates everywhere. We show that, for any given initial state, (a) and (b) produce a distribution of states approaching the Haar distribution in the limit of large circuit depth. More precisely, we show that the moments of the distribution produced by our circuits can approximate the ones of the Haar distribution in a depth proportional to the system size. Interestingly we find that in both Cases (a) and (b) the relaxation to the Haar distribution occurs in two steps—this is in contrast with what happens in fully random circuits. Moreover, we show that choosing appropriately the fixed interactions, for example taking the local gate to be a dual-unitary gate with high enough entangling power, minimally random circuits produce a Haar random distribution more rapidly than fully random circuits. In particular, dual-unitary circuits with maximal entangling power—i.e. perfect tensors—appear to provide the optimal quantum state design preparation for any design number. Our results also suggest that the minimally random settings (a) and (b) give an optimal framework to prepare unitary designs, i.e., random unitary matrices.

Simulating open quantum systems is an essential technique for understanding complex physical phenomena and advancing quantum technologies. Some quantum algorithms simulate Lindblad dynamics exponentially accurately, i.e., they achieve logarithmically short circuit depth in terms of accuracy, but they need to coherently encode all possible jump operators with a large ancilla consumption. Minimizing the gate and ancilla counts while achieving such a logarithmic scaling in accuracy remains an important challenge. In this work, we present two randomized quantum algorithms for simulating general Lindblad dynamics with multiple jump operators aimed at an observable estimation that achieves a circuit depth with not only logarithmic scaling in accuracy but also either partial or complete independence from the parameters specifying the Lindbladian. This is based on a novel random circuit compilation method that leverages dissipative processes with only a single jump operator, leading to the proposed methods using minimal ancilla qubits—$4+\lceil {\log }_{2}M\rceil$ in the first case and seven in the other, where each single jump operator has at most $M$ Pauli strings. In addition, we numerically demonstrate the practical advantage over existing approaches by providing a detailed analysis of the required gate and ancilla counts. This work represents a significant step toward making open quantum system simulations more feasible on early fault-tolerant quantum computing devices.

To successfully perform quantum computations, it is often necessary to first accurately characterize the noise in the underlying hardware. However, it is well known that fundamental limitations prevent the unique identification of the noise model. This raises the question of whether these limitations impact the ability to predict noisy dynamics and mitigate errors. Here, we show, both theoretically and experimentally, that when learnable parameters are self-consistently characterized, the unlearnable (gauge) degrees of freedom do not impact predictions of noisy dynamics or error mitigation. We use the recently introduced framework of gate set Pauli noise learning to efficiently and self-consistently characterize and mitigate noise of a complete gate set, including state preparation, measurements, single-qubit gates and multi-qubit entangling Clifford gates. We validate our approach through experiments with up to 92 qubits and show that, while the gauge choice does not affect error-mitigated observables, optimizing it reduces sampling overhead.

Entanglement purification protocols (EPPs) are essential for generating high-fidelity entangled states in noisy quantum systems, enabling robust quantum networking and computation. Building on the circuit of the foundational recurrence protocol, we generalize two-way EPPs to arbitrary stabilizer codes. Through analytical derivations and noisy-circuit simulations incorporating circuit-level noise, we demonstrate enhanced purification performance, with fidelity improvements and finite distillation rates for distillable input states. We propose efficient circuit designs for EPPs tailored to dual-species Rydberg atom arrays, leveraging species-specific laser control and interspecies Rydberg interactions. Introducing a low-overhead operation set, the dual-species atom convenient operation set, we facilitate straightforward compilation of EPP circuits without the need for ancillary atoms or complex atom rearrangements. Our framework provides practical guidance for near-term implementations on dual-species platforms, advancing toward scalable entanglement distribution in neutral atom systems and paving the way for fault-tolerant quantum technologies.

$\mathrm{Ge}/\mathrm{SiGe}$ quantum well heterostructures confining a high-mobility two-dimensional hole gas (2DHG) have emerged as a compelling platform for hybrid superconductor(S)-semiconductor(Sm) quantum devices. Here, we investigate the low-temperature transport properties of split-gate quantum point contacts (QPC) defined in one such heterostructure and positioned at different distances from an aluminum superconducting contact. We observe ballistic one-dimensional transport evidenced by conductance quantization with at least four clearly visible plateaus. Andreev reflection at the S/Sm interface induces a 40% enhancement of the conductance steps relative to the normal-state conductance staircase measured under a 100-mT out-of-plane magnetic field. This result is in excellent agreement with the theoretical expectation for an interface transparency of 0.88. By operating the QPCs in the tunneling regime, we probe the local density of states of the proximitized 2DHG. We report direct experimental evidence of an induced superconducting gap, demonstrating that its magnitude can be tuned by a gate voltage acting on the carrier density in the 2DHG.

Measurement and readout speed impose limits for logical qubit and qudit performance in quantum computing. Here, we experimentally demonstrate a universal analytical method for designing bandwidth-efficient dispersive readout pulses that intrinsically reset the resonator and extend it to multiple quantum states. The approach is based on universal analytical pulses and requires knowledge of the qubit and resonator parameters, but needs no direct optimization of the pulse shape, even when accounting for the nonlinearity of the system. When correcting for the Kerr nonlinearity, we find even more precise solutions to the speed-versus-photon-number trade-off, achieving ${10}^2$ photon readout with reset to ${10}^{-2}$ within $3{\kappa }^{-1}$ duration and amplifier efficiency- and $T_1$-limited assignment error around 0.6%. Results on qubits and qutrits validate the approach.

We introduce a holographic framework for analyzing the steady states of repeated quantum channels with strong symmetries. Using channel-state duality, we show that the steady state of a $d$-dimensional quantum channel is holographically mapped to the boundary reduced density matrix of a $(d+1)$-dimensional wavefunction generated by a sequential unitary circuit. From this perspective, strong-to-weak spontaneous symmetry breaking (SW-SSB) in the steady state arises from the anyon condensation on the boundary of a topological order in one higher dimension. The conditional mutual information (CMI) associated with SW-SSB is then inherited from the bulk topological entanglement entropy. We make this duality explicit using isometric tensor network states (isoTNS) by identifying the channel’s time evolution with the transfer matrix of a higher-dimensional isoTNS. Built on isoTNS, we further construct continuously tunable quantum channels that exhibit distinct mixed-state phases and transitions in the steady states.

Achieving ultrafast single-qubit gates that approach the limits set by decoherence requires operating in the strong-driving regime, where conventional semi-classical descriptions and single-leakage models break down, and multi-photon transitions emerge as dominant error channels that grow rapidly as the gate time is reduced. This makes it essential to develop systematic methods to characterize and suppress these errors using pulse shapes that remain simple and experimentally straightforward to calibrate. In contrast to standard derivative removal by adiabatic gate (DRAG) approaches for single-qubit gates, which suppress leakage within an effective three-level description, our framework systematically treats multi-photon transitions that arise in the strong-driving regime due to the participation of higher energy levels. In particular, we identify and suppress the two-photon $|1⟩\leftrightarrow |3⟩$ channel, incorporate time-ordering corrections, and derive a recursive analytical construction that captures the dominant higher-order error processes. We term the resulting pulse families R1D and R2D. These constructions extend and unify DRAG-type approaches by incorporating higher-order processes and time-ordering effects, yielding analytical expressions for optimal quadrature corrections, detunings, and prefactors. In addition, we clarify the optimal choice of DRAG prefactors and constant detuning when time-ordering effects are properly taken into account. This provides a systematic prescription for calibrating additional correction terms to further improve gate performance. Using this framework, we numerically demonstrate gate infidelities below ${10}^{-5}$ for a $7\mathrm{ns}$ full $\pi$-rotation under realistic decoherence rates.

Time-dependent electromagnetic drives are fundamental for controlling complex quantum systems, including superconducting Josephson circuits. In these devices, accurate time-dependent Hamiltonian models are imperative for predicting their dynamics and designing high-fidelity quantum operations. Existing numerical methods, such as black-box quantization (BBQ) and energy-participation ratio (EPR), excel at modeling the static Hamiltonians of Josephson circuits. However, these techniques do not fully capture the behavior of driven circuits stimulated by external microwave drives, nor do they include a generalized approach to account for the inevitable noise and dissipation that enter through microwave ports. Here, we introduce numerical techniques that leverage classical microwave simulations, efficiently executable in finite-element solvers, to obtain the time-dependent Hamiltonian of microwave-driven superconducting circuits with arbitrary geometries under charge, flux, or mixed electromagnetic modulation. Importantly, our techniques do not rely on a lumped-element description of the superconducting circuit, in contrast to previous approaches to tackling this problem. We demonstrate the versatility of our approach by characterizing the driven properties of realistic circuit devices in complex electromagnetic environments, including coherent dynamics due to charge and flux modulation, as well as drive-induced relaxation and dephasing. Our techniques offer a powerful toolbox for optimizing circuit designs and advancing practical applications in superconducting quantum computing.