A key goal of quantum chaos is to establish a relationship between widely observed universal spectral fluctuations of clean quantum systems and random matrix theory (RMT). Most prominent features of such RMT behavior with respect to a random spectrum, both encompassed in the spectral pair correlation function, are statistical suppression of small level spacings (correlation hole) and enhanced stiffness of the spectrum at large spectral ranges. For single-particle systems with fully chaotic classical counterparts, the problem has been partly solved by Berry [Proc. R. Soc. A 400, 229 (1985)] within the so-called diagonal approximation of semiclassical periodic-orbit sums, while the derivation of the full RMT spectral form factor $K(t)$ (Fourier transform of the spectral pair correlation function) from semiclassics has been completed by Müller et al. [Phys. Rev. Lett. 93, 014103 (2004)]. In recent years, the questions of long-time dynamics at high energies, for which the full many-body energy spectrum becomes relevant, are coming to the forefront even for simple many-body quantum systems, such as locally interacting spin chains. Such systems display two universal types of behaviour which are termed the “many-body localized phase” and “ergodic phase.” In the ergodic phase, the spectral fluctuations are excellently described by RMT, even for very simple interactions and in the absence of any external source of disorder. Here we provide a clear theoretical explanation for these observations. We compute $K(t)$ in the leading two orders in $t$ and show its agreement with RMT for nonintegrable, time-reversal invariant many-body systems without classical counterparts, a generic example of which are Ising spin-$1/2$ models in a periodically kicking transverse field. In particular, we relate $K(t)$ to partition functions of a class of twisted classical Ising models on a ring of size $t$; hence, the leading-order RMT behavior $K(t)\simeq 2t$ is a consequence of translation and reflection symmetry of the Ising partition function.

Verification is a task to check whether a given quantum state is close to an ideal state or not. In this paper, we show that a variety of many-qubit quantum states can be verified with only sequential single-qubit measurements of Pauli operators. First, we introduce a protocol for verifying ground states of Hamiltonians. We next explain how to verify quantum states generated by a certain class of quantum circuits. We finally propose an adaptive test of stabilizers that enables the verification of all polynomial-time-generated hypergraph states, which include output states of the Bremner-Montanaro-Shepherd-type instantaneous quantum polynomial time (IQP) circuits. Importantly, we do not make any assumption that the identically and independently distributed copies of the same states are given: Our protocols work even if some highly complicated entanglement is created among copies in any artificial way. As applications, we consider the verification of the quantum computational supremacy demonstration with IQP models, and verifiable blind quantum computing.

Quantum systems can be exquisite sensors thanks to their sensitivity to external perturbations. This same characteristic also makes them fragile to external noise. Quantum control can tackle the challenge of protecting a quantum sensor from environmental noise, while strongly coupling the sensor with the field to be measured. As the compromise between these two conflicting requirements does not always have an intuitive solution, optimal control based on a numerical search could prove very effective. Here, we adapt optimal control theory to the quantum-sensing scenario by introducing a cost function that, unlike the usual fidelity of operation, correctly takes into account both the field to be measured and the environmental noise. We experimentally implement this novel control paradigm using a nitrogen vacancy center in diamond, finding improved sensitivity to a broad set of time-varying fields. The demonstrated robustness and efficiency of the numerical optimization, as well as the sensitivity advantage it bestows, will prove beneficial to many quantum-sensing applications.

To implement fault-tolerant quantum computation with continuous variables, the Gottesman-Kitaev-Preskill (GKP) qubit has been recognized as an important technological element. However, it is still challenging to experimentally generate the GKP qubit with the required squeezing level, 14.8 dB, of the existing fault-tolerant quantum computation. To reduce this requirement, we propose a high-threshold fault-tolerant quantum computation with GKP qubits using topologically protected measurement-based quantum computation with the surface code. By harnessing analog information contained in the GKP qubits, we apply analog quantum error correction to the surface code. Furthermore, we develop a method to prevent the squeezing level from decreasing during the construction of the large-scale cluster states for the topologically protected, measurement-based, quantum computation. We numerically show that the required squeezing level can be relaxed to less than 10 dB, which is within the reach of the current experimental technology. Hence, this work can considerably alleviate this experimental requirement and take a step closer to the realization of large-scale quantum computation.

We introduce the “displacemon” electromechanical architecture that comprises a vibrating nanobeam, e.g., a carbon nanotube, flux coupled to a superconducting qubit. This platform can achieve strong and even ultrastrong coupling, enabling a variety of quantum protocols. We use this system to describe a protocol for generating and measuring quantum interference between trajectories of a nanomechanical resonator. The scheme uses a sequence of qubit manipulations and measurements to cool the resonator, to apply two effective diffraction gratings, and then to measure the resulting interference pattern. We demonstrate the feasibility of generating a spatially distinct quantum superposition state of motion containing more than $10^6$ nucleons using a vibrating nanotube acting as a junction in this new superconducting qubit configuration.

Stabilizer codes are among the most successful quantum error-correcting codes, yet they have important limitations on their ability to fault tolerantly compute. Here, we introduce a new quantity, the disjointness of the stabilizer code, which, roughly speaking, is the number of mostly nonoverlapping representations of any given nontrivial logical Pauli operator. The notion of disjointness proves useful in limiting transversal gates on any error-detecting stabilizer code to a finite level of the Clifford hierarchy. For code families, we can similarly restrict logical operators implemented by constant-depth circuits. For instance, we show that it is impossible, with a constant-depth but possibly geometrically nonlocal circuit, to implement a logical non-Clifford gate on the standard two-dimensional surface code.

Isolated quantum many-body systems with integrable dynamics generically do not thermalize when taken far from equilibrium. As one perturbs such systems away from the integrable point, thermalization sets in, but the nature of the crossover from integrable to thermalizing behavior is an unresolved and actively discussed question. We explore this question by studying the dynamics of the momentum distribution function in a dipolar quantum Newton’s cradle consisting of highly magnetic dysprosium atoms. This is accomplished by creating the first one-dimensional Bose gas with strong magnetic dipole-dipole interactions. These interactions provide tunability of both the strength of the integrability-breaking perturbation and the nature of the near-integrable dynamics. We provide the first experimental evidence that thermalization close to a strongly interacting integrable point occurs in two steps: prethermalization followed by near-exponential thermalization. Exact numerical calculations on a two-rung lattice model yield a similar two-timescale process, suggesting that this is generic in strongly interacting near-integrable models. Moreover, the measured thermalization rate is consistent with a parameter-free theoretical estimate, based on identifying the types of collisions that dominate thermalization. By providing tunability between regimes of integrable and nonintegrable dynamics, our work sheds light on the mechanisms by which isolated quantum many-body systems thermalize and on the temporal structure of the onset of thermalization.

The quantum Rabi model, involving a two-level system and a bosonic field mode, is arguably the simplest and most fundamental model describing quantum light-matter interactions. Historically, due to the restricted parameter regimes of natural light-matter processes, the richness of this model has been elusive in the lab. Here, we experimentally realize a quantum simulation of the quantum Rabi model in a single trapped ion, where the coupling strength between the simulated light mode and atom can be tuned at will. The versatility of the demonstrated quantum simulator enables us to experimentally explore the quantum Rabi model in detail, including a wide range of otherwise unaccessible phenomena, as those happening in the ultrastrong and deep strong-coupling regimes. In this sense, we are able to adiabatically generate the ground state of the quantum Rabi model in the deep strong-coupling regime, where we are able to detect the nontrivial entanglement between the bosonic field mode and the two-level system. Moreover, we observe the breakdown of the rotating-wave approximation when the coupling strength is increased, and the generation of phonon wave packets that bounce back and forth when the coupling reaches the deep strong-coupling regime. Finally, we also measure the energy spectrum of the quantum Rabi model in the ultrastrong-coupling regime.

The eigenstate thermalization hypothesis (ETH) posits that the reduced density matrix for a subsystem corresponding to an excited eigenstate is “thermal.” Here we expound on this hypothesis by asking: For which class of operators, local or nonlocal, is ETH satisfied? We show that this question is directly related to a seemingly unrelated question: Is the Hamiltonian of a system encoded within a single eigenstate? We formulate a strong form of ETH where, in the thermodynamic limit, the reduced density matrix of a subsystem corresponding to a pure, finite energy density eigenstate asymptotically becomes equal to the thermal reduced density matrix, as long as the subsystem size is much less than the total system size, irrespective of how large the subsystem is compared to any intrinsic length scale of the system. This allows one to access the properties of the underlying Hamiltonian at arbitrary energy densities (or temperatures) using just a single eigenstate. We provide support for our conjecture by performing an exact diagonalization study of a nonintegrable 1D quantum lattice model with only energy conservation. In addition, we examine the case in which the subsystem size is a finite fraction of the total system size, and we find that, even in this case, many operators continue to match their canonical expectation values, at least approximately. In particular, the von Neumann entanglement entropy equals the thermal entropy as long as the subsystem is less than half the total system. Our results are consistent with the possibility that a single eigenstate correctly predicts the expectation values of all operators with support on less than half the total system, as long as one uses a microcanonical ensemble with vanishing energy width for comparison. We also study, both analytically and numerically, a particle-number conserving model at infinite temperature that substantiates our conjectures.

We address the question of whether the super-Heisenberg scaling for quantum estimation is indeed realizable. We unify the results of two approaches. In the first one, the original system is compared with its copy rotated by the parameter-dependent dynamics. If the parameter is coupled to the one-body part of the Hamiltonian, the precision of its estimation is known to scale at most as $N^{-1}$ (Heisenberg scaling) in terms of the number of elementary subsystems used $N$. The second approach compares the overlap between the ground states of the parameter-dependent Hamiltonian in critical systems, often leading to an apparent super-Heisenberg scaling. However, we point out that if one takes into account the scaling of time needed to perform the necessary operations, i.e., ensuring adiabaticity of the evolution, the Heisenberg limit given by the rotation scenario is recovered. We illustrate the general theory on a ferromagnetic Heisenberg spin chain example and show that it exhibits such super-Heisenberg scaling of ground-state fidelity around the critical value of the parameter (magnetic field) governing the one-body part of the Hamiltonian. Even an elementary estimator represented by a single-site magnetization already outperforms the Heisenberg behavior providing the $N^{-1.5}$ scaling. In this case, Fisher information sets the ultimate scaling as $N^{-1.75}$, which can be saturated by measuring magnetization on all sites simultaneously. We discuss universal scaling predictions of the estimation precision offered by such observables, both at zero and finite temperatures, and support them with numerical simulations in the model. We provide an experimental proposal of realization of the considered model via mapping the system to ultracold bosons in a periodically shaken optical lattice. We explicitly derive that the Heisenberg limit is recovered when the time needed for preparation of quantum states involved is taken into account.

We study the dynamics of a supersonically expanding, ring-shaped Bose-Einstein condensate both experimentally and theoretically. The expansion redshifts long-wavelength excitations, as in an expanding universe. After expansion, energy in the radial mode leads to the production of bulk topological excitations—solitons and vortices—driving the production of a large number of azimuthal phonons and, at late times, causing stochastic persistent currents. These complex nonlinear dynamics, fueled by the energy stored coherently in one mode, are reminiscent of a type of “preheating” that may have taken place at the end of inflation.

Nonlocality and contextuality are at the root of conceptual puzzles in quantum mechanics, and they are key resources for quantum advantage in information-processing tasks. Bell nonlocality is best understood as the incompatibility between quantum correlations and the classical theory of causality, applied to relativistic causal structure. Contextuality, on the other hand, is on a more controversial foundation. In this work, I provide a common conceptual ground between nonlocality and contextuality as violations of classical causality. First, I show that Bell inequalities can be derived solely from the assumptions of no signaling and no fine-tuning of the causal model. This removes two extra assumptions from a recent result from Wood and Spekkens and, remarkably, does not require any assumption related to independence of measurement settings—unlike all other derivations of Bell inequalities. I then introduce a formalism to represent contextuality scenarios within causal models and show that all classical causal models for violations of a Kochen-Specker inequality require fine-tuning. Thus, the quantum violation of classical causality goes beyond the case of spacelike-separated systems and already manifests in scenarios involving single systems.

Unitary transformations are routinely modeled and implemented in the field of quantum optics. In contrast, nonunitary transformations, which can involve loss and gain, require a different approach. In this work, we present a universal method to deal with nonunitary networks. An input to the method is an arbitrary linear transformation matrix of optical modes that does not need to adhere to bosonic commutation relations. The method constructs a transformation that includes the network of interest and accounts for full quantum optical effects related to loss and gain. Furthermore, through a decomposition in terms of simple building blocks, it provides a step-by-step implementation recipe, in a manner similar to the decomposition by Reck et al. [Experimental Realization of Any Discrete Unitary Operator, Phys. Rev. Lett. 73, 58 (1994)] but applicable to nonunitary transformations. Applications of the method include the implementation of positive-operator-valued measures and the design of probabilistic optical quantum information protocols.

Random quantum circuits yield minimally structured models for chaotic quantum dynamics, which are able to capture, for example, universal properties of entanglement growth. We provide exact results and coarse-grained models for the spreading of operators by quantum circuits made of Haar-random unitaries. We study both $1+1\mathrm{D}$ and higher dimensions and argue that the coarse-grained pictures carry over to operator spreading in generic many-body systems. In $1+1\mathrm{D}$, we demonstrate that the out-of-time-order correlator (OTOC) satisfies a biased diffusion equation, which gives exact results for the spatial profile of the OTOC and determines the butterfly speed $v_B$. We find that in $1+1\mathrm{D}$, the “front” of the OTOC broadens diffusively, with a width scaling in time as $t^{1/2}$. We address fluctuations in the OTOC between different realizations of the random circuit, arguing that they are negligible in comparison to the broadening of the front within a realization. Turning to higher dimensions, we show that the averaged OTOC can be understood exactly via a remarkable correspondence with a purely classical droplet growth problem. This implies that the width of the front of the averaged OTOC scales as $t^{1/3}$ in $2+1\mathrm{D}$ and as $t^{0.240}$ in $3+1\mathrm{D}$ (exponents of the Kardar-Parisi-Zhang universality class). We support our analytic argument with simulations in $2+1\mathrm{D}$. We point out that, in two or higher spatial dimensions, the shape of the spreading operator at late times is affected by underlying lattice symmetries and, in general, is not spherical. However, when full spatial rotational symmetry is present in $2+1\mathrm{D}$, our mapping implies an exact asymptotic form for the OTOC, in terms of the Tracy-Widom distribution. For an alternative perspective on the OTOC in $1+1\mathrm{D}$, we map it to the partition function of an Ising-like statistical mechanics model. As a result of special structure arising from unitarity, this partition function reduces to a random walk calculation which can be performed exactly. We also use this mapping to give exact results for entanglement growth in $1+1\mathrm{D}$ circuits.

We generate and characterize entangled states of a register of 20 individually controlled qubits, where each qubit is encoded into the electronic state of a trapped atomic ion. Entanglement is generated amongst the qubits during the out-of-equilibrium dynamics of an Ising-type Hamiltonian, engineered via laser fields. Since the qubit-qubit interactions decay with distance, entanglement is generated at early times predominantly between neighboring groups of qubits. We characterize entanglement between these groups by designing and applying witnesses for genuine multipartite entanglement. Our results show that, during the dynamical evolution, all neighboring qubit pairs, triplets, most quadruplets, and some quintuplets simultaneously develop genuine multipartite entanglement. Witnessing genuine multipartite entanglement in larger groups of qubits in our system remains an open challenge.

One of the main aims in the field of quantum simulation is to achieve a quantum speedup, often referred to as “quantum computational supremacy,” referring to the experimental realization of a quantum device that computationally outperforms classical computers. In this work, we show that one can devise versatile and feasible schemes of two-dimensional, dynamical, quantum simulators showing such a quantum speedup, building on intermediate problems involving nonadaptive, measurement-based, quantum computation. In each of the schemes, an initial product state is prepared, potentially involving an element of randomness as in disordered models, followed by a short-time evolution under a basic translationally invariant Hamiltonian with simple nearest-neighbor interactions and a mere sampling measurement in a fixed basis. The correctness of the final-state preparation in each scheme is fully efficiently certifiable. We discuss experimental necessities and possible physical architectures, inspired by platforms of cold atoms in optical lattices and a number of others, as well as specific assumptions that enter the complexity-theoretic arguments. This work shows that benchmark settings exhibiting a quantum speedup may require little control, in contrast to universal quantum computing. Thus, our proposal puts a convincing experimental demonstration of a quantum speedup within reach in the near term.

Manipulating the state of a logical quantum bit (qubit) usually comes at the expense of exposing it to decoherence. Fault-tolerant quantum computing tackles this problem by manipulating quantum information within a stable manifold of a larger Hilbert space, whose symmetries restrict the number of independent errors. The remaining errors do not affect the quantum computation and are correctable after the fact. Here we implement the autonomous stabilization of an encoding manifold spanned by Schrödinger cat states in a superconducting cavity. We show Zeno-driven coherent oscillations between these states analogous to the Rabi rotation of a qubit protected against phase flips. Such gates are compatible with quantum error correction and hence are crucial for fault-tolerant logical qubits.

Single-photon detection is an essential component in many experiments in quantum optics, but it remains challenging in the microwave domain. We realize a quantum nondemolition detector for propagating microwave photons and characterize its performance using a single-photon source. To this aim, we implement a cavity-assisted conditional phase gate between the incoming photon and a superconducting artificial atom. By reading out the state of this atom in a single shot, we reach an external (internal) photon-detection fidelity of 50% (71%), limited by transmission efficiency between the source and the detector (75%) and the coherence properties of the qubit. By characterizing the coherence and average number of photons in the field reflected off the detector, we demonstrate its quantum nondemolition nature. We envisage applications in generating heralded remote entanglement between qubits and for realizing logic gates between propagating microwave photons.

The noncommutativity of position and momentum observables is a hallmark feature of quantum physics. However, this incompatibility does not extend to observables that are periodic in these base variables. Such modular-variable observables have been suggested as tools for fault-tolerant quantum computing and enhanced quantum sensing. Here, we implement sequential measurements of modular variables in the oscillatory motion of a single trapped ion, using state-dependent displacements and a heralded nondestructive readout. We investigate the commutative nature of modular variable observables by demonstrating no-signaling in time between successive measurements, using a variety of input states. Employing a different periodicity, we observe signaling in time. This also requires wave-packet overlap, resulting in quantum interference that we enhance using squeezed input states. The sequential measurements allow us to extract two-time correlators for modular variables, which we use to violate a Leggett-Garg inequality. Signaling in time and Leggett-Garg inequalities serve as efficient quantum witnesses, which we probe here with a mechanical oscillator, a system that has a natural crossover from the quantum to the classical regime.

We investigate the transport properties of neutral, fermionic atoms passing through a one-dimensional quantum wire containing a mesoscopic lattice. The lattice is realized by projecting individually controlled, thin optical barriers on top of a ballistic conductor. Building an increasingly longer lattice, one site after another, we observe and characterize the emergence of a band insulating phase, demonstrating control over quantum-coherent transport. We explore the influence of atom-atom interactions and show that the insulating state persists as contact interactions are tuned from moderately to strongly attractive. Using bosonization and classical Monte Carlo simulations, we analyze such a model of interacting fermions and find good qualitative agreement with the data. The robustness of the insulating state supports the existence of a Luther-Emery liquid in the one-dimensional wire. Our work realizes a tunable, site-controlled lattice Fermi gas strongly coupled to reservoirs, which is an ideal test bed for nonequilibrium many-body physics.