We introduce a family of nonintegrable 1D lattice models that feature robust periodic revivals under a global quench from certain initial product states, thus generalizing the phenomenon of many-body scarring recently observed in Rydberg atom quantum simulators. Our construction is based on a systematic embedding of the single-site unitary dynamics into a kinetically constrained many-body system. We numerically demonstrate that this construction yields new families of models with robust wave-function revivals, and it includes kinetically constrained quantum clock models as a special case. We show that scarring dynamics in these models can be decomposed into a period of nearly free clock precession and an interacting bottleneck, shedding light on their anomalously slow thermalization when quenched from special initial states.

In the presence of conservation laws, superpositions of eigenstates of the corresponding conserved quantities cannot be generated by quantum dynamics. Thus, any such coherence represents a potentially valuable resource of asymmetry, which can be used, for example, to enhance the precision of quantum metrology or to enable state transitions in quantum thermodynamics. Here we ask if such superpositions, already present in a reference system, can be broadcast to other systems, thereby distributing asymmetry indefinitely at the expense of creating correlations. We prove a no-go theorem showing that this is forbidden by quantum mechanics in every finite-dimensional system. In doing so, we also answer some open questions in the quantum information literature concerning the sharing of timing information of a clock and the possibility of catalysis in quantum thermodynamics. We also prove that even weaker forms of broadcasting, of which Åberg’s “catalytic coherence” is a particular example, can only occur in the presence of infinite-dimensional reference systems. Our results set fundamental limits on the creation and manipulation of quantum coherence and shed light on the possibilities and limitations of quantum reference frames to act catalytically without being degraded.

Symmetries of both closed- and open-system dynamics imply many significant constraints. These generally have instantiations in both classical and quantum dynamics (Noether’s theorem, for instance, applies to both sorts of dynamics). We here provide an example of such a constraint which has no counterpart for a classical system, that is, a uniquely quantum consequence of symmetric dynamics. Specifically, we demonstrate the impossibility of broadcasting asymmetry (symmetry breaking) relative to a continuous symmetry group, for bounded-size quantum systems. The no-go theorem states that if two initially uncorrelated systems interact by symmetric dynamics and asymmetry is created at one subsystem, then the asymmetry of the other subsystem must be reduced. We also find a quantitative relation describing the trade-off between the subsystems. These results cannot be understood in terms of additivity of asymmetry, because, as we show here, any faithful measure of asymmetry violates both subadditivity and superadditivity. Rather, they must be understood as a consequence of an (intrinsically quantum) information-disturbance principle. Our result also implies that if a bounded-size quantum reference frame for the symmetry group, or equivalently, a bounded-size reservoir of coherence (e.g., a clock with coherence between energy eigenstates in quantum thermodynamics) is used to implement any operation that is not symmetric, then the quantum state of the frame or reservoir is necessarily disturbed in an irreversible fashion, i.e., degraded.

Computing marginal distributions of discrete or semidiscrete Markov random fields (MRFs) is a fundamental, generally intractable problem with a vast number of applications in virtually all fields of science. We present a new family of computational schemes to approximately calculate the marginals of discrete MRFs. This method shares some desirable properties with belief propagation, in particular, providing exact marginals on acyclic graphs, but it differs with the latter in that it includes some loop corrections; i.e., it takes into account correlations coming from all cycles in the factor graph. It is also similar to the adaptive Thouless-Anderson-Palmer method, but it differs with the latter in that the consistency is not on the first two moments of the distribution but rather on the value of its density on a subset of values. The results on finite-dimensional Isinglike models show a significant improvement with respect to the Bethe-Peierls (tree) approximation in all cases and with respect to the plaquette cluster variational method approximation in many cases. In particular, for the critical inverse temperature ${\beta }_c$ of the homogeneous hypercubic lattice, the expansion of ${(d{\beta }_c)}^{-1}$ around $d=\infty$ of the proposed scheme is exact up to $d^{-4}$ order, whereas the latter two are exact only up to $d^{-2}$ order.

We revisit a family of integrals that delude intuition and that recently appeared in mathematical literature in connection with computer algebra package verification. We show that the remarkable properties displayed by these integrals become transparent when formulated in the language of random walks. In turn, the random walk view naturally leads to a plethora of nontrivial generalizations that are worked out. Related complex identities are also derived, without the need of explicit calculation. The crux of our treatment lies in a causality argument where a message that travels at finite speed signals the existence of a boundary.

We introduce quantum hypercube states, a class of continuous-variable quantum states that are generated as orthographic projections of hypercubes onto the quadrature phase space of a bosonic mode. In addition to their interesting geometry, hypercube states display phase-space features much smaller than Planck’s constant, and a large volume of Wigner negativity. We theoretically show that these features make hypercube states sensitive to displacements at extremely small scales in a way that is surprisingly robust to initial thermal occupation and to small separation of the superposed state components. In a high-temperature proof-of-principle optomechanics experiment we observe, and match to theory, the signature outer-edge vertex structure of hypercube states.

A new distribution for systems of particles in equilibrium obeying the exclusion of correlated states is presented following Haldane’s state counting. It relies upon an ansatz to deal with the multiple exclusion that takes place when the states accessible to single particles are spatially correlated and it can be simultaneously excluded by more than one particle. Haldane’s statistics and Wu’s distribution are recovered in the limit of noncorrelated states of the multiple exclusion statistics. In addition, an exclusion spectrum function $\mathcal{G}(n)$ is introduced to account for the dependence of the state exclusion on the occupation number $n$. The results of thermodynamics and state occupation are shown for ideal lattice gases of linear particles of size $k$ ($k$-mers) where the multiple exclusion occurs. Remarkable agreement is found with grand-canonical Monte Carlo simulations from $k=2$ to 10 where the multiple exclusion dominates as $k$ increases.

We investigate how the presence of a single-particle mobility edge in a system can generate strong energy current rectification. Specifically, we study a quadratic bosonic chain subject to a quasiperiodic potential and coupled at its boundaries to spin baths of differing temperature. We find that rectification increases by orders of magnitude depending on the spatial position in the chain of localized eigenstates above the mobility edge. The largest enhancements occur when the coupling of one bath to the system is dominated by a localized eigenstate, while the other bath couples to numerous delocalized eigenstates. By tuning the parameters of the quasiperiodic potential it is thus possible to vary the amplitude, and even invert the direction, of the rectification.

In both thermodynamics and quantum mechanics, the arrow of time is characterized by the statistical likelihood of physical processes. We characterize this arrow of time for the continuous quantum measurement dynamics of a superconducting qubit. By experimentally tracking individual weak measurement trajectories, we compare the path probabilities of forward and backward-in-time evolution to develop an arrow of time statistic associated with measurement dynamics. We compare the statistics of individual trajectories to ensemble properties showing that the measurement dynamics obeys both detailed and integral fluctuation theorems, thus establishing the consistency between microscopic and macroscopic measurement dynamics.

We propose a new cellular automaton (CA), the sweep rule, which generalizes Toom’s rule to any locally Euclidean lattice. We use the sweep rule to design a local decoder for the toric code in $d\ge 3$ dimensions, the sweep decoder, and rigorously establish a lower bound on its performance. We also numerically estimate the sweep decoder threshold for the three-dimensional toric code on the cubic and body-centered cubic lattices for phenomenological phase-flip noise. Our results lead to new CA decoders with provable error-correction thresholds for other topological quantum codes including the color code.

We construct efficient deterministic dynamical decoupling schemes protecting continuous-variable degrees of freedom. Our schemes target decoherence induced by quadratic system-bath interactions with analytic time dependence. We show how to suppress such interactions to $N$th order using only $N$ pulses. Furthermore, we show how to homogenize a $2^m$-mode bosonic system using only $(N+1{)}^{2m+1}$ pulses, yielding—up to the $N$th order—an effective evolution described by noninteracting harmonic oscillators with identical frequencies. The decoupled and homogenized system provides natural decoherence-free subspaces for encoding quantum information. Our schemes only require pulses which are tensor products of single-mode passive Gaussian unitaries and swap gates between pairs of modes.

Chaotic dynamics in quantum many-body systems scrambles local information so that at late times it can no longer be accessed locally. This is reflected quantitatively in the out-of-time-ordered correlator of local operators, which is expected to decay to 0 with time. However, for systems of finite size, out-of-time-ordered correlators do not decay exactly to 0 and in this paper we show that the residual value can provide useful insights into the chaotic dynamics. When energy is conserved, the late-time saturation value of the out-of-time-ordered correlator of generic traceless local operators scales as an inverse polynomial in the system size. This is in contrast to the inverse exponential scaling expected for chaotic dynamics without energy conservation. We provide both analytical arguments and numerical simulations to support this conclusion.

Computing the spectra of operators is a fundamental problem in the sciences, with wide-ranging applications in condensed-matter physics, quantum mechanics and chemistry, statistical mechanics, etc. While there are algorithms that in certain cases converge to the spectrum, no general procedure is known that (a) always converges, (b) provides bounds on the errors of approximation, and (c) provides approximate eigenvectors. This may lead to incorrect simulations. It has been an open problem since the 1950s to decide whether such reliable methods exist at all. We affirmatively resolve this question, and the algorithms provided are optimal, realizing the boundary of what digital computers can achieve. Moreover, they are easy to implement and parallelize, offer fundamental speed-ups, and allow problems that before, regardless of computing power, were out of reach. Results are demonstrated on difficult problems such as the spectra of quasicrystals and non-Hermitian phase transitions in optics.

In a many-body quantum system, local operators in the Heisenberg picture $O(t)=e^{iHt}O{e}^{-iHt}$ spread as time increases. Recent studies have attempted to find features of that spreading which could distinguish between chaotic and integrable dynamics. The operator entanglement—the entanglement entropy in operator space—is a natural candidate to provide such a distinction. Indeed, while it is believed that the operator entanglement grows linearly with time $t$ in chaotic systems, we present evidence that it grows only logarithmically in generic interacting integrable systems. Although this logarithmic growth has been previously established for noninteracting fermions, there has been no progress on interacting integrable systems to date. In this Letter we provide an analytical upper bound on operator entanglement for all local operators in the “Rule 54” qubit chain, a cellular automaton model introduced in the 1990s [Bobenko et al., CMP 158, 127 (1993)], and recently advertised as the simplest representative of interacting integrable systems. Physically, the logarithmic bound originates from the fact that the dynamics of the models is mapped onto the one of stable quasiparticles that scatter elastically. The possibility of generalizing this scenario to other interacting integrable systems is briefly discussed.