In the realm of fault-tolerant quantum computing, stabilizer operations play a pivotal role, characterized by their remarkable efficiency in classical simulation. This efficiency sets them apart from nonstabilizer operations within the quantum computational theory. In this Letter, we investigate the limitations of classically simulable measurements in distinguishing quantum states. We demonstrate that any pure magic state and its orthogonal complement of odd prime dimensions cannot be unambiguously distinguished by stabilizer operations, regardless of how many copies of the states are supplied. We also reveal intrinsic similarities and distinctions between the quantum resource theories of magic states and entanglement in quantum state discrimination. The results emphasize the inherent limitations of classically simulable measurements and contribute to a deeper understanding of the quantum-classical boundary.

Digital quantum simulation relies on Trotterization to discretize time evolution into elementary quantum gates. On current quantum processors with notable gate imperfections, there is a critical trade-off between improved accuracy for finer time steps, and increased error rate on account of the larger circuit depth. We present an adaptive Trotterization algorithm to cope with time dependent Hamiltonians, where we propose a concept of piecewise “conserved” quantities to estimate errors in the time evolution between two (nearby) points in time; these allow us to bound the errors accumulated over the full simulation period. They reduce to standard conservation laws in the case of time independent Hamiltonians, for which we first developed an adaptive Trotterization scheme [H. Zhao et al., Making Trotterization adaptive and energy-self-correcting for NISQ devices and beyond, PRX Quantum 4, 030319 (2023).]. We validate the algorithm for a time dependent quantum spin chain, demonstrating that it can outperform the conventional Trotter algorithm with a fixed step size at a controlled error.

Pulsar timing arrays perform Bayesian posterior inference with expensive Markov chain Monte Carlo (MCMC) methods. Given a dataset of $\sim 10–100$ pulsars and $\mathcal{O}({10}^3)$ timing residuals each, producing a posterior distribution for the stochastic gravitational wave background (SGWB) can take days to a week. The computational bottleneck arises because the likelihood evaluation required for MCMC is extremely costly when considering the dimensionality of the search space. Fortunately, generating simulated data is fast, so modern simulation-based inference techniques can be brought to bear on the problem. In this Letter, we demonstrate how conditional normalizing flows trained on simulated data can be used for extremely fast and accurate estimation of the SGWB posteriors, reducing the sampling time from weeks to a matter of seconds.

Absolute angular rotation rate measurements with sensitivity better than prad/s would be beneficial for fundamental science investigations. In this regard, large frame Earth based ring laser gyroscopes are top instrumentation as far as bandwidth, long-term operation, and sensitivity are concerned. Here, we demonstrate that the GINGERINO active-ring laser upper limiting noise is close to $2\times {10}^{-15}\mathrm{rad}/\mathrm{s}$ for $\sim 2\times {10}^5\mathrm{s}$ of integration time, as estimated by the Allan deviation evaluated in a differential measurement scheme. This result is more than a factor of 10 better than the theoretical prediction so far accounted for ideal ring lasers shot-noise with the two beams counterpropagating inside the cavity considered as two independent propagating modes. This feature is related to the peculiarity of real ring laser system dynamics that causes phase crosstalking among the two counterpropagating modes. In this context, the independent beam model is, then, not applicable, and the measured noise limit falls below the expected one.

The magnetic skyrmions generated in a centrosymmetric crystal were recently first discovered in ${\mathrm{Gd}}_2{\mathrm{PdSi}}_3$. In light of this, we observe the electronic structure by angle-resolved photoemission spectroscopy and unveil its direct relationship with the magnetism in this compound. The Fermi surface and band dispersions are demonstrated to have a good agreement with the density functional theory calculations carried out with careful consideration of the crystal superstructure. Most importantly, we find that the three-dimensional Fermi surface has extended nesting which matches well the $\mathit{q}$ vector of the magnetic order detected by recent scattering measurements. The consistency we find among angle-resolved photoemission spectroscopy, density functional theory, and the scattering measurements suggests the Ruderman-Kittel-Kasuya-Yosida interaction involving itinerant electrons to be the formation mechanism of skyrmions in ${\mathrm{Gd}}_2{\mathrm{PdSi}}_3$.

The spin degrees of freedom is crucial for the understanding of any condensed matter system. Knowledge of spin-mixing mechanisms is not only essential for successful control and manipulation of spin qubits, but also uncovers fundamental properties of investigated devices and material. For electrostatically defined bilayer graphene quantum dots, in which recent studies report spin-relaxation times $T_1$ up to 50 ms with strong magnetic field dependence, we study spin-blockade phenomena at charge configuration $(1,2)\leftrightarrow (0,3)$. We examine the dependence of the spin-blockade leakage current on interdot tunnel coupling and on the magnitude and orientation of externally applied magnetic field. In out-of-plane magnetic field, the observed zero-field current peak could arise from finite-temperature cotunneling with the leads; though involvement of additional spin- and valley-mixing mechanisms are necessary for explaining the persistent sharp side peaks observed. In in-plane magnetic field, we observe a zero-field current dip, attributed to the competition between the spin Zeeman effect and the Kane-Mele spin-orbit interaction. Details of the line shape of this current dip, however, suggest additional underlying mechanisms are at play.

Nonstabilizerness, also known as “magic,” stands as a crucial resource for achieving a potential advantage in quantum computing. Its connection to many-body physical phenomena is poorly understood at present, mostly due to a lack of practical methods to compute it at large scales. We present a novel approach for the evaluation of nonstabilizerness within the framework of matrix product states (MPSs), based on expressing the MPS directly in the Pauli basis. Our framework provides a powerful tool for efficiently calculating various measures of nonstabilizerness, including stabilizer Rényi entropies, stabilizer nullity, and Bell magic, and enables the learning of the stabilizer group of an MPS. We showcase the efficacy and versatility of our method in the ground states of Ising and XXZ spin chains, as well as in circuits dynamics that has recently been realized in Rydberg atom arrays, where we provide concrete benchmarks for future experiments on logical qubits up to twice the sizes already realized.

We present a novel classical algorithm designed to learn the stabilizer group—namely, the group of Pauli strings for which a state is a $\pm 1$ eigenvector—of a given matrix product state (MPS). The algorithm is based on a clever and theoretically grounded biased sampling in the Pauli (or Bell) basis. Its output is a set of independent stabilizer generators whose total number is directly associated with the stabilizer nullity, notably a well-established nonstabilizer monotone. We benchmark our method on $T$-doped states randomly scrambled via Clifford unitary dynamics, demonstrating very accurate estimates up to highly entangled MPS with bond dimension $\chi \sim {10}^3$. Our method, thanks to a very favorable scaling $\mathcal{O}({\chi }^3)$, represents the first effective approach to obtain a genuine magic monotone for MPS, enabling systematic investigations of quantum many-body physics out of equilibrium.

Recent advances in quantum simulation based on neutral atoms have largely benefited from high-resolution, single-atom sensitive imaging techniques. A variety of approaches have been developed to achieve such local detection of atoms in optical lattices or optical tweezers. For alkaline-earth and alkaline-earth-like atoms, the presence of narrow optical transitions opens up the possibility of performing novel types of Sisyphus cooling, where the cooling mechanism originates from the capability to spatially resolve the differential optical level shifts in the trap potential. Up to now, it has been an open question whether high-fidelity imaging could be achieved in a “repulsive Sisyphus” configuration, where the trap depth of the ground state exceeds that of the excited state involved in cooling. Here, we demonstrate high-fidelity (99.971(1)%) and high-survival (99.80(5)%) imaging of strontium atoms using repulsive Sisyphus cooling. We use an optical lattice as a pinning potential for atoms in a large-scale tweezer array with up to 399 tweezers and show repeated, high-fidelity lattice-tweezer-lattice transfers. We furthermore demonstrate loading the lattice with approximately 10 000 atoms directly from the MOT and scalable imaging over $>10000$ lattice sites with a combined survival probability and classification fidelity better than 99.2%. Our lattice thus serves as a locally addressable and sortable reservoir for continuous refilling of optical tweezer arrays in the future.

We highlight a noncanonical yet natural choice of variables for an efficient derivation of a kinetic equation for the energy density in nonisotropic systems, including internal gravity waves on a vertical plane, inertial, and Rossby waves. The existence of a second quadratic invariant simplifies the kinetic equation and leads to extra conservation laws for resonant interactions. We analytically determine the scaling of the radial turbulent energy spectrum. Our findings suggest the existence of an inverse energy cascade of internal gravity waves, from small to large scales, in practically relevant scenarios.

The state of an open quantum system undergoing an adiabatic process evolves by following the instantaneous stationary state of its time-dependent generator. This observation allows one to characterize, for a generic adiabatic evolution, the average dynamics of the open system. However, information about fluctuations of dynamical observables, such as the number of photons emitted or the time-integrated stochastic entropy production in single experimental runs, requires controlling the whole spectrum of the generator and not only the stationary state. Here, we show how such information can be obtained in adiabatic open quantum dynamics by exploiting tools from large deviation theory. We prove an adiabatic theorem for deformed generators, which allows us to encode, in a biased quantum state, the full counting statistics of generic time-integrated dynamical observables. We further compute the probability associated with an arbitrary “rare” time history of the observable and derive a dynamics which realizes it in its typical behavior. Our results provide a way to characterize and engineer adiabatic open quantum dynamics and to control their fluctuations.

We demonstrate how to incorporate a catalyst to enhance the performance of a heat engine. Specifically, we analyze efficiency in one of the simplest engine models, which operates in only two strokes and comprises of a pair of two-level systems, potentially assisted by a $d$-dimensional catalyst. When no catalysis is present, the efficiency of the machine is given by the Otto efficiency. Introducing the catalyst allows for constructing a protocol which overcomes this bound, while new efficiency can be expressed in a simple form as a generalization of Otto’s formula: $1-(1/d)({\omega }_c/{\omega }_h)$. The catalyst also provides a bigger operational range of parameters in which the machine works as an engine. Although an increase in engine efficiency is mostly accompanied by a decrease in work production (approaching zero as the system approaches Carnot efficiency), it can lead to a more favorable trade-off between work and efficiency. The provided example introduces new possibilities for enhancing performance of thermal machines through finite-dimensional ancillary systems.

Twin-field quantum key distribution (TFQKD) overcomes the linear rate-loss limit, which promises a boost of secure key rate over long distance. However, the complexity of eliminating the frequency differences between the independent laser sources hinders its practical application. We analyzed and determined the frequency stability requirements for implementing TFQKD using frequency-stabilized lasers. Based on this analysis, we proposed and demonstrated a simple and practical approach that utilizes the saturated absorption spectroscopy of acetylene as an absolute reference, eliminating the need for fast frequency locking to achieve TFQKD. Adopting the 4-intensity sending-or-not-sending TFQKD protocol, we experimentally demonstrated the TFQKD over 502, 301, and 201 km ultralow-loss optical fiber, respectively. We expect this high-performance scheme will find widespread usage in future intercity and free-space quantum communication networks.

The superconformal index of half-BPS states in $\mathcal{N}=4$ supersymmetric Yang-Mills with gauge group $\mathrm{U}(N)$ admits an expansion in terms of giant gravitons, ${\mathcal{I}}_N(q)={\mathcal{I}}_{\infty }(q){\sum }_{m=0}^{\infty }{q}^{mN}{\widehat{\mathcal{I}}}_m(q)$, where $m$ is the number of giant gravitons and ${\mathcal{I}}_{\infty }(q)$ is the graviton index. The expansion can be viewed as the implementation of trace relations for finite $N$. We derive this expansion directly in supergravity from the class of half-BPS solutions due to Lin, Lunin, and Maldacena in type IIB supergravity. The moduli space of these configurations can be quantized using covariant quantization methods. We show how this quantization leads to the precise expression for the expansion in terms of giant gravitons. Our proposal provides a derivation of the giant graviton expansion directly in terms of quantized supergravity degrees of freedom, and it recovers discrete data via quantum geometries that are classically nonsmooth.

Resonant positron annihilation on atomic electrons provides a powerful method to search for light new particles coupled to $e^{+}{e}^{-}$. Reliable estimates of production rates require a detailed characterization of electron momentum distributions. We describe a general method that harnesses the target material Compton profile to properly include electron velocity effects in resonant annihilation cross sections. We additionally find that high-$Z$ atoms can efficiently act as particle physics accelerators, providing a density of relativistic electrons that allows one to extend by several times the experimental mass reach.