Quantum process learning is a fundamental primitive that draws inspiration from machine learning with the goal of better studying the dynamics of quantum systems. One approach to quantum process learning is quantum compilation, whereby an analog quantum operation is digitized by compiling it into a series of basic gates. While there has been significant focus on quantum compiling for discrete-variable systems, the continuous-variable (CV) framework has received comparatively less attention. We present an experimental implementation of a CV quantum compiler that uses two-mode squeezed light to learn a Gaussian unitary operation. We demonstrate the compiler by learning a parameterized linear phase unitary through the use of target and control phase unitaries to demonstrate a factor of 5.4 increase in the precision of the phase estimation and a 3.6-fold acceleration in the time-to-solution metric when leveraging quantum resources. We further show how our approach can be extended to higher-dimensional compilation tasks. Our results are enabled by the tunable control of our cost landscape via variable squeezing, thus providing a critical framework to simultaneously increase precision and reduce time-to-solution.

Bosonic codes in superconducting resonators are a hardware-efficient avenue for quantum error correction and benefit from the inherent bias toward relaxation errors provided by long-lived cavities compared to typical superconducting qubits. The required coupling to an ancillary control qubit, however, can negate these benefits by inducing highly detrimental effects such as excess decoherence and undesired nonlinearities. This is exemplified by the typically employed transmon qubit that imposes limits on cavity operations due to its finite lifetime as well as the inevitably induced self-Kerr nonlinearity. An important question is thus whether a cavity-qubit coupling can be realized that offers readout and control capabilities without spoiling the cavity. Here, motivated by its long lifetime and design flexibility of its Hamiltonian, we experimentally investigate the fluxonium as a control qubit for superconducting cavities. We couple a fluxonium qubit to a superconducting resonator in the strong-dispersive regime and use it to measure the coherence and flux-tunable inherited nonlinearities of the resonator. We then demonstrate quantum control of the oscillator by preparing and characterizing resonator Fock states and their superpositions, with fidelities limited by resonator decay in our planar prototype device. Finally, we show that the fluxonium can realize cavity-couplings that eliminate undesirable cavity nonlinearities and reach regimes of suppressed self-Kerr that are inaccessible to transmon-based devices. Our results thus demonstrate the potential of the fluxonium as a high-performance bosonic control qubit for superconducting cavities that can outperform existing transmon-based architectures.

In this paper, we combine the $k$-reduction map, the moment method, and randomized measurements into a practical protocol for certifying the entanglement dimensionality. Our approach is based on the observation that a state with entanglement dimensionality at most $k$ must stay positive under the action of the $k$-reduction map. The core of our protocol utilizes the moment method to determine whether the $k$-reduced operator, i.e., the operator obtained after applying the $k$-reduction map on a quantum state, contains a negative eigenvalue or not. Notably, we propose a systematic method for constructing reduction moment criteria that apply to a broader range of states, including $k$-unfaithful states, compared with fidelity-based methods. The performance of our approach gets better and better with the moment order employed, which is corroborated by extensive numerical simulations. To apply our approach, it suffices to implement a unitary 3-design instead of a 4-design, which is more feasible in practice than the correlation matrix method. In the course of study, we show that the $k$-reduction negativity, the absolute sum of the negative eigenvalues of the $k$-reduced operator, is monotonic under local operations and classical communication for pure states.

Quantum state tomography serves as a key tool for identifying quantum states generated in quantum computers and simulators, typically involving local operations on individual particles or qubits to enable independent measurements. However, this approach requires an exponentially larger number of measurement setups as quantum platforms grow in size, highlighting the necessity of more scalable methods to efficiently perform quantum state estimation. Here, we present a tomography scheme that scales far more efficiently and, remarkably, eliminates the need for local addressing of single constituents before measurements. Inspired by the “spin-spiral” structure in magnetic materials, our scheme combines a series of measurement setups, each with different spiraling patterns, with compressed-sensing techniques. The results of the numerical simulations demonstrate a high degree of tomographic efficiency and accuracy. Additionally, we show how this method is suitable for the measurement of specific entanglement properties of interesting quantum many-body states, such as entanglement entropy, under various realistic experimental conditions. This method offers a positive outlook across a wide range of quantum platforms, including those in which precise individual operations are challenging, such as optical lattice systems.

In classical physics, events follow a definite causal order: the past influences the future, but not the reverse. Quantum theory, however, permits superpositions of causal orders—the so-called indefinite causal orders (ICOs)—which can provide operational advantages over classical scenarios. Verifying such phenomena has sparked significant interest, much like earlier efforts devoted to refuting local realism and confirming quantum entanglement. To date, demonstrations of ICO have all been based a process called the quantum switch and have relied on device-dependent or semi-device-independent protocols. Achieving a device-independent verification of ICO would imply that nature allows for correlations that do not respect causality, independent of any experimental assumptions or underlying theoretical description of the experiment. To this end, a recent theoretical development introduced a Bell-like inequality that allows for fully device-independent verification of ICO in a quantum switch. Here we implement this verification by experimentally violating this inequality. In particular, we measure a value of $1.8328\pm 0.0045$, which is 18 standard deviations above the definite causal order bound of $1.75$. Our work presents the first implementation of a device-independent protocol to verify ICO, albeit in the presence of experimental loopholes. This represents an important step toward the device-independent verification of an ICO and provides a context in which to identify loopholes specifically related to the verification of ICO.

Magic state cultivation is a leading approach for generating the resource states required for fault-tolerant quantum computation. Here we present a cultivation protocol that increases the success probability of magic state generation in platforms with flexible and nonlocal connectivity. Our method implements cultivation directly on the surface code, avoiding the detour through alternative, less efficient error-correcting codes used in prior approaches. This both improves the acceptance rate and preserves compatibility with the geometry of the code and the hardware. Numerical simulations show that our protocols improve success probabilities and reduce output error rates compared with protocols tailored to locally connected platforms. Under realistic noise models for cold-atom and trapped-ion systems, we improve the rate of magic state generation by more than a factor of 20. Finally, we study qubit loss and erasure and show that very low error rates can be achieved with minimal overhead, reaching below ${10}^{-6}$ infidelity using only nine physical erasure qubits.

Kinetically constrained models were originally introduced to capture slow relaxation in glassy systems, where dynamics are hindered by local constraints instead of energy barriers. Their quantum counterparts have recently drawn attention for exhibiting highly degenerate eigenstates at zero energy—known as zero modes—stemming from chiral symmetry. Yet, the structure and implications of these zero modes remain poorly understood. In this work, we focus on the properties of the zero mode subspace in quantum kinetically constrained models with a $U(1)$ particle-conservation symmetry. We use the $U(1)$ East, which lacks inversion symmetry, and the inversion-symmetric $U(1)$ East-West models to illustrate our two main results. First, we observe that the simultaneous presence of constraints and chiral symmetry generally leads to a parametric increase in the number of zero modes due to the fragmentation of the many-body Hilbert space into disconnected sectors. Second, we generalize the concept of compact localized states from single-particle physics and introduce the notion of collective bound states, a special kind of nonergodic eigenstates that are robust to enlarging the system size. We formulate sufficient criteria for their existence, arguing that the degenerate zero mode subspace plays a central role, and demonstrate bound states in both example models and in a two-dimensional model, the $U(1)$ North-East, and in the pair-flip model, a system without particle conservation. Our results motivate a systematic study of bound states and their relation to ergodicity breaking, transport, and other properties of quantum kinetically constrained models.

Consider the problem of deciding, for a particular multipartite quantum state, whether or not it is realizable in a quantum network with a particular causal structure. This is a fully quantum version of what causal inference researchers refer to as the problem of causal discovery. In this work, we introduce a fully quantum version of the inflation technique for causal inference, which leverages the quantum marginal problem. The primary example by which we illustrate the utility of this method is testing compatibility of tripartite quantum states with the quantum network known as the triangle scenario. We show, in particular, how the method yields a complete classification of pure three-qubit states into those that are and those that are not compatible with the triangle scenario. We also provide some illustrative examples involving mixed-states and some where one or more of the systems is higher dimensional. Furthermore, we examine the question of when the incompatibility of a multipartite quantum state with a causal structure can be inferred from the incompatibility of a joint probability distribution induced by implementing measurements on each subsystem. Finally, we present a family of networks, which includes the triangle scenario as a special case, for which causal compatibility constraints can be derived.

This work introduces an approach rooted in quantum thermodynamics to enhance sampling efficiency in quantum machine learning (QML). We propose conceptualizing quantum supervised learning as a thermodynamic cooling process. Building on this concept, we develop a quantum refrigerator protocol that enhances sample efficiency during training and prediction without the need for Grover iterations or quantum phase estimation. Inspired by heat-bath algorithmic cooling protocols, our method alternates entropy compression and thermalization steps to decrease the entropy of qubits, increasing polarization toward the dominant bias. This technique minimizes the computational overhead associated with estimating classification scores and gradients, presenting a practical and efficient solution for QML algorithms compatible with noisy intermediate-scale quantum devices.

Quantum speed limits are usually regarded as fundamental restrictions, constraining the amount of computation that can be achieved within some given time and energy. Complementary to this intuition, here we show that these limitations are also of operational value: they enable the secure generation of certified randomness. We consider a prepare-and-measure scenario with some (experimentally determined or promised) upper bound on the energy uncertainty $\mathrm{\Delta }E$ of the average prepared quantum state, but without any further assumptions on the devices, Hilbert space or Hamiltonian. Given that we can freely choose the time at which to apply the untrusted preparation procedure, we show that this scenario admits the generation of randomness that is secure against adversaries with additional classical information. We show how to determine the amount of certified randomness given the observed correlations, discuss how interactions with the environment are taken into account, and sketch a conceivable experimental implementation. In particular, we show that single-mode coherent states admit this kind of certification of nonzero randomness in some parameter regimes, reinforcing existing demonstrations of nonclassicality in the simple harmonic oscillator. Our results extend earlier efforts to devise semi-device-independent protocols grounded in reasonable physical assumptions, and they contribute to the understanding of time-energy uncertainty relations via their operational consequences.

Achieving robust and scalable remote quantum entanglement is one of the fundamental challenges for the development of distributed quantum networks and modular quantum computing systems. In this context, perfect state transfer (PST) and fractional state transfer (FST) have emerged as promising schemes for quantum state transfer and remote entanglement generation using only nearest-neighbor couplings. In this work, we propose a quantum state transfer scheme based on a zig-zag configuration, which introduces a control parameter for PST and FST. We show that this new parameter can suppress the population in the intermediate qubits, thereby reducing losses and enhancing state transfer. In certain limiting cases, our new scheme reduces to the conventional PST scheme, revealing an elegant mathematical structure inherent in the design. We experimentally demonstrate our proposed scheme on a superconducting quantum processor, achieving an 18% reduction in error for remote Bell state generation in a 1D ($1\times 5$) qubit chain with enhanced robustness to specific noise channels. Furthermore, we extend our approach to a 2D ($3\times 3$) network and successfully generate a $\mathcal{W}$ state among the four corner qubits. These results highlight the potential of our new enhanced quantum state transfer scheme for scalable and noise-resilient quantum communication and computing.

Recent advances have demonstrated that $\mathcal{O}(\log M)$ measurements suffice to predict $M$ properties of arbitrarily large quantum many-body systems. However, these remarkable findings assume that the properties to be predicted are chosen independently of the data. This assumption can be violated in practice, where scientists adaptively select properties after looking at previous predictions. This work investigates the adaptive setting for three classes of observables: local, Pauli, and bounded-Frobenius-norm observables. We prove that $\mathrm{\Omega }(\sqrt{M})$ samples of an arbitrarily large unknown quantum state are necessary to predict expectation values of $M$ adaptively chosen local and Pauli observables, where the system size scales exponentially and polynomially in $M$, respectively. We also present computationally efficient algorithms that achieve this information-theoretic lower bound. In contrast, for bounded-Frobenius-norm observables, we devise an algorithm requiring only $\mathcal{O}(\log M)$ samples, independent of system size. These results highlight the potential pitfalls of adaptivity in analyzing data from quantum experiments and provide algorithmic tools to safeguard against erroneous predictions in quantum experiments.

Continuously monitored quantum systems are emerging as promising platforms for quantum metrology, where a central challenge is to identify measurement strategies that optimally extract information about unknown parameters encoded in the complex quantum state of emitted radiation. Different measurement strategies effectively access distinct temporal modes of the emitted field, and the resulting choice of mode can strongly impact the information available for parameter estimation. While a ubiquitous approach in quantum optics is to select frequency modes through spectral filtering, the metrological potential of this technique has not yet been systematically quantified. We develop a theoretical framework to assess this potential by modeling spectral detection as a cascaded quantum system, allowing us to reconstruct the full density matrix of frequency-filtered photonic modes and to compute their associated Fisher information. This framework provides a minimal yet general method to benchmark the performance of spectral measurements in quantum optics, allowing to identify optimal filtering strategies in terms of frequency selection, detector linewidth, and metrological gain accessible through higher-order frequency-resolved correlations and mean-field engineering. These results lay the groundwork for identifying and designing optimal sensing strategies in practical quantum-optical platforms.

In this work we demonstrate the use of adapted classical phase retrieval algorithms to perform control-free quantum phase estimation. We eliminate the costly controlled time evolution and Hadamard test commonly required to access the complex time-series needed to reconstruct the spectrum. This significant reduction of the number of coherent controlled-operations lowers the circuit depth and considerably simplifies the implementation of statistical quantum phase estimation in near-term devices. This seemingly impossible task can be achieved by extending the problem that one wishes to solve to one with a larger set of input signals while exploiting natural constraints on the signal and/or the spectrum. We leverage well-established algorithms that are widely used in the context of classical signal processing, demonstrating two complementary methods to do this, vectorial phase retrieval and two-dimensional phase retrieval. We numerically investigate the feasibility of both approaches for estimating the spectrum of the Fermi-Hubbard model and discuss their resilience to inherent statistical noise.

We introduce magic-augmented Clifford circuits—architectures in which Clifford circuits are preceded and/or followed by constant-depth circuits of non-Clifford (“magic”) gates—as a resource-efficient way to realize approximate $k$-designs, with reduced circuit depth and usage of magic. We prove that shallow-Clifford circuits, when augmented with constant-depth circuits of magic gates, can generate approximate unitary and state $k$-designs with $\epsilon$ relative error. The total circuit depth for these constructions on $N$ qubits is $O(\log (N/\epsilon ))+2^{O(k\log k)}$ in one dimension and $O(\log \log (N/\epsilon ))+2^{O(k\log k)}$ in all-to-all circuits using ancillas, which improves upon previous results for small $k\ge 4$. Furthermore, our construction of relative-error state $k$-designs only involves states with strictly local magic. The required number of magic gates is parametrically reduced when considering $k$-designs with bounded additive error. As an example, we show that shallow-Clifford circuits followed by $O(k^2)$ single-qubit magic gates, independent of system size, can generate an additive-error state $k$-design. We develop a classical statistical mechanics description of our random circuit architectures, which provides a quantitative understanding of the required depth and number of magic gates for additive-error state $k$-designs. We also prove no-go theorems for various architectures to generate designs with bounded relative error.

We provide practical simulation methods for scalar field theories on a quantum computer that yield improved asymptotics as well as concrete gate estimates for the simulation and physical qubit estimates using the surface code. We achieve these improvements through two optimizations. First, we consider a finite volume approach for estimating the elements of the S-matrix. This approach is appropriate in general for 1+1D and for certain low-energy elastic collisions in higher dimensions. Second, we implement our approach using a series of different fault-tolerant simulation algorithms for Hamiltonians formulated both in the field occupation basis and field amplitude basis. Our algorithms are based on either second-order Trotterization or qubitization. The cost of Trotterization in occupation basis scales as $O(\lambda N^7|\mathrm{\Omega }{|}^3/(M^{5/2}{ϵ}^{3/2}))$ where $\lambda$ is the coupling strength, $N$ is the occupation cutoff, $|\mathrm{\Omega }|$ is the volume of the spatial lattice, $M$ is the mass of the particles and $ϵ$ is the uncertainty in the energy calculation used for the $S$-matrix determination. Qubitization in the field basis scales as $O(|\mathrm{\Omega }{|}^2(k^2\mathrm{\Lambda }+k{M}^2)/ϵ)$, where $k$ is the cutoff in the field and $\mathrm{\Lambda }$ is a scaled coupling constant. We find in both cases that the bounds suggest physically meaningful simulations can be performed using on the order of $4\times {10}^6$ physical qubits and ${10}^{12}T$-gates which corresponds to roughly one day on a superconducting quantum computer with surface code and a cycle time of 100 ns. This places the simulation of scalar field theory within striking distance of the gate counts for the best available chemistry simulation results.

We present a practical approach for interfacing light with a two-dimensional atomic tweezer array. Typical paraxial fields are poorly matched to the array’s multi-diffraction-order radiation pattern, thus severely limiting the interface coupling efficiency. Instead, we propose to design a field mode that naturally couples to the array: it consists of a unique superposition of multiple beams corresponding to the array’s diffraction orders. This composite mode can be generated from a single Gaussian beam using standard free-space optics, including spatial light modulators and a single objective lens. For a triangular array with lattice spacing about twice the wavelength, all diffraction angles remain below ${35}^{\circ }$, making the scheme compatible with standard objectives of numerical aperture $\mathrm{NA}\le 0.7$. Our analytical theory and scattering simulations reveal that the interface efficiency $r_0$ for quantum information tasks scales favorably with the array atom number $N$: reaching $>0.99$ ($>0.9999$) for $N=149$ ($N\sim 1000$) and scaling as $1-r_0\sim 1/N$ for large $N$. The scheme is robust to optical imperfections and atomic-position errors, offering a viable path for quantum light-matter applications in current tweezer-array platforms.

We develop a framework for constructing quantum error-correcting codes and logical gates for three types of spaces—composite permutation-invariant spaces of many qubits or qudits, composite constant-excitation Fock state spaces of many bosonic modes, and monolithic nuclear state spaces of atoms, ions, and molecules. By identifying all three spaces with discrete simplices and representations of the Lie group $SU(q)$, we prove that many codes and their gates in $SU(q)$ can be interconverted between the three state spaces. We construct code instances for all three spaces using classical ${\ell }_1$ codes and Tverberg’s theorem, a classic result from convex geometry. We obtain families of quantum codes with distance that scales almost linearly with the code length $N$ by constructing ${\ell }_1$ codes based on combinatorial patterns called Sidon sets and utilizing their Tverberg partitions. This compares favorably with the existing designs for all the state spaces. We present explicit constructions of codes with shorter length or lower total spin/excitation than known codes with similar parameters, bosonic codes with exotic Gaussian gates, as well as examples of short codes with distance larger than the known constructions.

Divisibility of dynamical maps is a central notion in the study of quantum non-Markovianity, providing a natural framework to characterize memory effects via time-local master equations. In this work, we generalize the notion of divisibility of quantum dynamical maps from the Schrödinger to the Heisenberg picture. While the two pictures are equivalent at the level of physical predictions, we show that the divisibility properties of the corresponding dual maps are, in general, not equivalent. This inequivalence originates from the distinction between left and right generators of time-local master equations, which interchange roles under duality. We demonstrate that Schrödinger and Heisenberg divisibility are distinct concepts by constructing explicit dynamics divisible only in one picture. Furthermore, we introduce a quantifier for the violation of Heisenberg P-divisibility, analogous to the trace-distance-based measure of non-Markovianity, and provide it with an operational interpretation in terms of the guessing probability between effects. Our results show that Heisenberg divisibility is an independent witness of memory effects and highlight the need to consider both pictures when characterizing non-Markovian quantum dynamics.

Recent demonstrations indicate that silicon-spin quantum processing unit will be able to shuttle physical qubits rapidly and with high fidelity—a desirable feature for maximizing logical connectivity, supporting new codes, and routing around damage. However it may seem that shuttling at the logical level is unwise: static defects in the device may “scratch” a logical qubit as it passes, causing correlated errors to which the code is highly vulnerable. Here we explore an architecture where logical qubits are 1D strings (“snakes”) which can be moved freely over a planar latticework. Possible scratch events are inferred via monitor qubits and the complementary gap; if deemed a risk, remarkably the shuttle process can be undone in a way that negates any corruption. This leads to high levels of tolerance against shuttling-related imperfections and enables logical operations between snakes by a semitransversal method. We conclude that this approach is suitable for fault-tolerant computing in both near-term and long-term, mature-era silicon devices.