Long-time dynamics in quantum spin lattices: ergodicity and hydrodynamic projections at all frequencies and wavelengths

Obtaining rigorous and general results about the non-equilibrium dynamics of extended many-body systems is a difficult task. In quantum lattice models with short-range interactions, the Lieb-Robinson bound tells us that the spatial extent of operators grows at most linearly in time. But what happens within this light-cone? We discuss rigorous results on ergodicity and the emergence of the hydrodynamic scale, which establish fundamental principles at the root of non-equilibrium physics. One key idea of the present work is that general structures of hydrodynamics at the Euler scale follow independently from the details of the microscopic dynamics, and in particular do not necessitate chaos; they are consequences of"extensivity". Another crucial observation is that these apply at arbitrary frequencies and wavelengths. That is, long-time oscillatory behaviours can be reproduced from a natural extension of standard hydrodynamic notions, thus enlarging the hydrodynamic paradigm beyond the zero-frequency / infinite-wavelength point that it traditionally addresses.


Introduction
The dynamics of many-body, extended, isolated, interacting quantum systems has been the focus of much recent theoretical and experimental research [1,2,3]. Of particular interest is their universal, large-scale behaviours away from equilibrium. Given the panoply of phenomena observed, it is crucial to rigorously establish the physical principles that guide them, such as ergodicity and hydrodynamics. Can we prove that certain forms of ergodicity occur, and that hydrodynamics emerges whereby a large amount of information is lost? Can we characterise the remaining degrees of freedom? For reversible, Hamiltonian dynamics in extended quantum systems, these are some of the deepest questions and remain largely open, forming an important part of Hilbert's sixth problem [4].
On the one hand, ergodicity is the idea that, over long times, the system covers uniformly enough the manifold of states, or at least the part of it that is dynamically accessible, such as the energy shell. In the language of Gibbs' statistical mechanics, this translates into the equivalence of time averages with ensemble averages. In what situations does this happen in extended quantum systems? Conventional wisdom says that this requires the presence of chaos: the exponential separation of initially nearby trajectories, or, in quantum mechanics, certain conditions on the energies and matrix elements of macroscopic observables [5,6,7], or on outof-time-ordered correlators (OTOC) [8,9,10,11]. Is there such a link between short-time (chaotic) and long-time (ergodic) behaviours?
On the other hand, hydrodynamics purports that the dynamics of the system reduces to that of a few emergent, slowly-decaying degrees of freedom. One stark realisation of this is the Boltzmann-Gibbs principle [12,13,14]: in hydrodynamic linear response, this is the idea that correlations due to local perturbations are carried by long-lived modes over which local perturbations "project", and that propagate at hydrodynamic velocities. Indeed, the dominant correlations between a person's vocal chords and another's eardrum are carried by sound waves. The principle holds quite generally, including in extended Hamiltonian quantum models, and including with integrability [15,16,17]. What are the conditions for the emergence of such hydrodynamic projections, and what are the hydrodynamic modes?
In this paper, we report on progress in these directions. First, after reviewing basic aspects of ergodicity including von Neumann's ergodic theorem, we discuss notions of many-body ergodicity. We show that a "non-localisation condition" (there is a t such that no local observable returns to itself under evolution by time t) implies time-ergodicity for local observables. We also show that it implies a stability theorem of Kubo-Martin-Schwinger (KMS) states with correlations that vanish quickly enough in space: after local perturbations, over long periods of time, the state returns 100% of the time arbitrarily close, in a metric based on local observables, to the KMS state. We then discuss the authors' recent "almost-everywhere ergodicity" result [18], which show that for displacements along space-time rays at almost every velocity, these conclusions in fact hold in all generality in short-range hypercubic quantum lattice models (it is always true that, for almost every speed v, there is a t such that no local observable returns to itself under evolution by time t and displacement by distance vt). In particular, infinite-time averages of local observables along rays in space-time of almost every velocity are non-fluctuating, and, as operators, are "thin", i.e. approach the identity (thus become "classical").
Next, we discuss the phenomenon of hydrodynamic projections. Using a form of almosteverywhere ergodicity, we prove hydrodynamic projection theorems in arbitrary dimension, and explain how such results immediately extend to arbitrary frequencies f and wavenumbers k. We give a precise definition of the space of hydrodynamic modes, as the kernel of a the (in general foscillatory) unitary evolution operator on the Hilbert space of (in general k-)extensive quantities. In general, this space depends on the state and on the chosen frequency and wavenumber. Again, results apply to all short-range hypercubic quantum lattices, and a large family of states.
Finally, oscillatory hydrodynamic projections give rise to the "oscillatory linearised Euler equation", and we illustrate in the free-fermionic chain how this reproduces the oscillatory behaviours of two-point functions at large space-time separations. This gives a "proof of principle" for the use of hydrodynamics to address such oscillatory behaviours at long times.
From the physics perspective, our main messages are: 1. In contrast to the few-body case, the "right" notion of ergodicity in many-body systems does not rely on the conventional idea of covering the energy shell. Rather, in the thermodynamic limit, the fact that timelike averages of local observables become non-fluctuating on large times follows from the interplay between extensivity of the system and locality of the observables. It happens on timescales which are much smaller than the Poincaré recurrence scale, or that for covering the energy shell. This is many-body ergodicity (related to what is sometimes called "typicality"). Interestingly, despite these physical differences, it is nevertheless the same mathematical principles from the basic ergodicity theory that underpins many-body and few-body ergodicity, including von Neumann's ergodic theorem.
2. Euler-scale hydrodynamics hold for a wide family of systems, independently from their specific dynamics (integrable or chaotic, constrained or not, etc.). General structures, such as hydrodynamic projections, always emerge at large scales, and only specific properties that control the phenomenology are model dependent. The idea that hydrodynamics is a useful concept beyond its conventional field of applications, not requiring chaos, has come to the fore with the recent development of generalised hydrodynamics for integrable systems [19,20,21,22], and ties in with the idea that thermodynamic concepts should not rely on short-time behaviours, as emphasised recently in [23]. However, rigorous results about many-body dynamics, especially of such generality, are notoriously difficult to obtain.
3. Hydrodynamics at the Euler scale can be generalised to arbitrary frequencies and wavelengths, going beyond the conventional zero-frequency and large-wavelength region. Thus, in D dimensions, there is a (f, k)-hydrodynamic theory for every frequency f ∈ R and wavenumber k ∈ R D . The (f, k)-hydrodynamic modes control the (f, k)-oscillatory behaviours of correlation functions at large separations in space-time (for any given (f, k), if the set of modes is non-empty, then such (f, k)-oscillatory behaviours are observed). This shows that recent ideas on dynamical symmetries [24,25,26,27,28] can be extended to the hydrodynamic regime of correlation functions. In fact, the proofs we provide make it clear that the extension goes beyond simple phase oscillations, and include oscillations involving nontrivial internal transformations, such as spin rotations. This in principle allows hydrodynamics to be extended to describe large-scale behaviours involving nontrivial microscopic internal structures, such as helicoidal spin structures.
From a more philosophical perspective, an important message is that it is extremely fruitful to extract the relevant many-body physics out of equilibrium by analysing the operator algebra of local observables already in the thermodynamic limit. This is to be contrasted with methods based, in quantum systems, on how the spectrum of large Hamiltonians behave as the thermodynamic limit is taken, or, in classical models, on kinetic equations. The operator algebra framework allows one to directly concentrate on the small part of the Hilbert space (or phase space in classical mechanics) that is significant in the thermodynamic limit. In particular, we emphasise how kernels of naturally defined unitary operators on Hilbert spaces of observables encode the relevant aspects of large-scale physics.
We believe most of these notions are valid in a large class of many-body systems, quantum or classical. In order to be precise, however, we will concentrate on translation-invariant, infinitelyextended D-dimensional hyper-cubic quantum lattices with finite local spaces and short-range interactions (interaction strength decaying at least exponentially with distance).
The paper is organised as follows. First, in Section 2, we discuss notions of ergodicity, and in particular many-body ergodicity, and "almost-everywhere ergodicity" established recently [18] within the context of the C * algebraic formulation of quantum statistical mechanics. We contrast this with von Neumann's ergodic theorem, and discuss application to the problem of return to equilibrium after local quenches. We then prove rigorously, in Section 3, the principle of hydrodynamic projection onto conserved modes in the Euler scaling limit, within the same context. This generalises the result of [29] to any dimension, although the proof requires the introduction of quite a few new ideas. Crucially, mathematical formalism make it clear that the general principle can take into account oscillatory behaviours at any frequency and wavelength. We give an accurate definition of the full, complete space of hydrodynamic modes, which in general depends on the state, and on the frequency and wavelength. Finally, in Section 4 we review the linearised (oscillatory) Euler hydrodynamic equations established in [29] in one  [18], and of the hydrodynamic projection results of section 3. Within the Lieb-Robinson cone |x| < v LR t, very little was known until now. We find that ray-averagesĀ T of observables up to time T , with respect to every frequency ω and wavelength k and at every velocity v, vanish at large times T → ∞ (represented by the colour-gradient shaded region), except perhaps for a set of v of measure zero (represented by solid ray segments). Observables become thin as time evolve. The large-scale information that remains, after oscillatory fluid-cell means, is their projection onto a smaller space of degrees of freedom Q ω,k , the slowly decaying hydrodynamic modes. dimension, and discuss how they recover oscillatory behaviours in free fermion models. The structures uncovered are pictorially illustrated in Figure 1. 2 Ergodicity, many-body ergodicity, almost-everywhere ergodicity Ergodicity plays a fundamental role in physics. In this section, we pedagogically review some concepts of ergodicity in many-body systems, one of the goals being to put in context some recent results by the authors. We first review the basic rigorous ergodicity result in classical mechanics from Birkhoff's ergodic theorem. We then review von Neumann's fundamental work on the quantum ergodic theorem, recasting the nice discussion in [6] into the language of local observables in quantum spin lattices. We observe that the eigenstate thermalisation hypothesis (ETH) is in fact closely related to the conditions of von Neumann's quantum ergodic theorem. Both the von Neumann quantum ergodic theorem and the ETH rely heavily on the structure of the spectrum, and how it behaves in the thermodynamic limit. Lastly, we avoid such difficulties by passing to the C * -algebra formalism, in which context we discuss what many-body ergodicity may look like, and review the recent results [18], where almost-everywhere ergodicity is proven for every quantum spin model on (hyper-)cubic lattices with short-range interactions; we discuss some applications and relations with the problem of return to equilibrium. We denote throughout this section the time average as

Ergodicity in classical mechanics
In the classical mechanics of interacting particles, the basic notion of ergodicity is the equivalence between time average and "ensemble average", or average over the uniform measure on the energy shell. Consider a system of N particles, with phase-space coordinates (x n , p n ) ∈ R 2D : n = 1, . . . , N and time evolution generated by a Hamiltonian H(x, p) (we denote (x, p) = ((x n , p n ) : n = 1, . . . , N ) ∈ R 2DN ). The microcanonical ensemble average of an observable A(x, p) on the shell Ω E = {(x, p) : H(x, p) = E} of energy E can be defined by a limit on the thickened shell Note that this naturally induces a measure µ E on Ω E . Now consider points (x, p) ∈ Ω E . The ergodicity statement is that, if there is no nontrivial constant of the motion 1 , then, for every observable A, every E ∈ im H, and almost every where (x(t), p(t)) is the point (x, p) time-evolved by time t. This is equivalent to a more intuitive statement: the fraction of time that a trajectory t → (x(t), p(t)) spends in a neighbourhood Γ ⊂ Ω E of a point (x , p ) ∈ Γ is, over long times, proportional to its "surface area" µ E (Γ) on the energy shell, for almost every initial point (x, p) and for every (x , p ). In particular, over long times, the trajectory comes as close as desired to the point (x , p ). Ergodicity in classical mechanics, as expressed above, follows immediately from Birkhoff's ergodic theorem [30] .

Von Neumann's quantum ergodic theorem
It is natural to ask about the equivalent ergodicity notion for quantum systems. A problem arises: in quantum mechanics, phase space is essentially "discretised" as dxdp /2. Now we can not ask a more precise question than if a trajectory will cross a neighbourhood of volume /2 of a point. Requiring that this happens for every point (in some precise way that can be worked out, see below) is a much weaker condition, that of "quantum recurrence". Of course, at large values of the classical action S , this quantisation is not seen, and one recovers classical mechanics. The above classical ergodicity notion, with the more stringent condition on the measure covered by trajectories, is then a good approximation. But how is this recovered directly from quantum mechanics? Large classical actions naturally occur in macroscopic systems. Can we formulate 1 A nontrivial constant of the motion is a measurable function Q on phase space (an observable) that is invariant under the dynamics (constant of the motion), and that breaks the energy shells (nontrivial) in that for at least one E, there exist two non-intersecting subsets Γ, Γ ⊂ ΩE of nonzero µE-measure such that Q(Γ) ∩ Q(Γ ) = ∅. ergodicity for macroscopic quantum systems, and find the quantum properties that guarantee it? In particular, timescales for quantum recurrence are prohibitively large, while those for ergodicity of the corresponding classical system at S are short; how to account for this? Von Neumann's idea is to look at "macroscopic observables". Intuitively, such observables should behave like classical ones, but one asks what the microscopic quantum evolution says about these observables, and in particular their ergodicity properties. Perhaps the original idea was to consider a classical object (such as a tennis ball) as a macroscopic combination of quantum particles, and macroscopic observables for such objects. But in fact, it turns out that the abstract formulation applies quite generally. Making connections with modern developments, we will recast von Neumann's line of thought to the context of many-body quantum systems with short-range interactions: what about ergodicity in large (macroscopic) quantum spin lattices?
Take a finite spin lattice: a quantum system with N -dimensional local spaces lying on where h x is an operator supported on a finite, x-independent number of sites near to x. For simplicity we take the system periodic and translation invariant: there is a representation y → ι y of the group (Z/(2L+1)Z) D on the group of automorphisms of the operator algebra Aut(End H), such that ι y h x = h x+y (and ι y A is supported on the y-translate of the support of A). The Hamiltonian has spectrum of eigenvalues {E n } and normalised eigenvectors |n . Consider a microcanonical shell S = {n : |Λ L |e − < E n < |Λ L |e + }, keeping the energy densities e − , e + implicit for readability. The microcanonical ensemble is represented by the density matrix It can be shown that [31], in the limit of large volume L → ∞, followed by the limit of small shell e − → e + → e, averages of local observables A (operators supported on finite numbers of sites) in this ensemble tend to averages in the Gibbs ensemble, where β is chosen such that Tr ρh 0 = e. A number of ergodicity-like results can be obtained depending on the strength of the assumptions on the underlying system.

Assumption of non-degeneracy and quantum recurrence
Before discussing von Neumann's QET, let us first consider the simplest ideas of ergodicity, keeping L finite and the shell of finite thickness. Let us assume that E n = E m if n = m. Take a state in the microcanonical shell, |ψ = n∈S c n |n .
Then clearly The result is within the so-called "diagonal ensemble": the ensemble of density matrices that are diagonal in the energy eigenbasis. This is certainly not ergodicity, even in the limit of large L and thin shell! The equivalent of the classical uniform shell-averaging is the uniform sum over n's (Tr ρ mc A), but here, the c n 's may behave wildely as functions of n.
One can nevertheless go further, and get a weak statement of ergodicity. Beyond the trivial stationarity statement n|A(t)|n = n|A|n , in any state |n , the observable A is non-fluctuating: This is "mean-square ergodicity", and in fact implies that the distribution of the random variable A(t) in the state |n is supported on A(t) = n|A|n . Mean-square ergodicity has a clear technical parallel with the classical ergodicity notion. Indeed, the condition of non-degeneracy is equivalent to the condition that the space of all operators invariant under time evolution is spanned by the projections P m = |m m|, and this parallels the classical condition that there be no nontrivial constants of the motion 2 . In evaluating time averages, it is von Neumann's ergodic theorem from functional analysis [32, Theorem II.11] (not to be confused with von Neumann's quantum ergodic theorem!) that replaces Birkhoff's theorem from measure theory. This theorem implies that, if U t is a unitary one-parameter group on a Hilbert space, then where P 1 the projector onto its unit-eigenvalue subspace. In order to use this, we must consider time evolution as a unitary operator acting on matrices (say under the Hilbert-Schmidt inner product), and then its long-time average projects onto span(P m ). Of course, in the finitedimensional context considered here, such more involved mathematics is not necessary: it is elementary to see that the long-time average of an observable reduces it to its diagonal, thus to a linear combination of projections P m . Then, because we look at a single state |n , this implies that observables are non-fluctuating. Thus, within a fix state |n , ergodicity holds, paralleling what happens within a fix classical energy shell (we fix the state |n , because the limit of an infinitesimally thin shell, for fixed L, is just a single state). As |n is stationary, the statement 2 In this respect, it is sometimes stated in the literature that integrability -in the sense of the existence of large enough amount of nontrivial constants of the motion -for finite quantum spin lattices does not make sense, because there always are as many projections Pm, which commute with H, as the dimension of the Hilbert space dim H. This is misleading, as a similar statement could be made for classical systems: functions of the Hamiltonian are constants of the motion. But they are not "nontrivial": they do not divide the energy shells. Similarly, Pm are not nontrivial constants of the motion. Only with degeneracies can new, nontrivial constants of the motion appear: the space of operators is of dimension (dim H) 2 , and integrability typically leads to enough degeneracies so that there are ∝ log(dim H) dim H linearly independent projections commuting with H; this is associated to the presence of a large enough amount of extensive conserved quantities. that A(t) be non-fluctuating and supported on n|A|n is the natural ergodicity statement.
But this parallel is essentially technical, as the above simple statements only relate to quantum recurrence, not many-body ergodicity. Indeed, the order of limits is simply incorrect for many-body systems: the limit in time T → ∞, in the time averaging (1), is taken before the limit of large volumes L → ∞ (which is in fact never taken in the mean-square ergodicity statement). This is thus valid, in realistic many-body systems, only at astronomically large times! As we will see in subsection 2.3, however, with the right space of operators, such a technical parallel becomes fully meaningful.

Assumptions on maximal deviations and von Neumann's QET
Assume that max n∈S n|A|n − Tr ρ mc A 2 is "small". A natural requirement would be that it is uniformly bounded in L (thus stays finite in the large volume limit), and this bound tends to zero in the limit of thin microcanonical shell. Of course, this cannot be expected to be valid for all observables A. But, at least for local observables, this should hold (certainly uniform boundedness in L is immediate by finiteness of the operator norm). A clear example is the energy density h x . The states |n may be assumed to diagonalise the translation operator, and therefore Thus and max n∈S n|h x |n − Tr ρ mc h x 2 ≤ (e + − e − ) 2 . The requirement is that such a smallness bound holds not just for the energy density, but for all (or at least for a large family of) local observables.
This looks much more like an ergodicity statement. Making connection with more modern ideas, we note that the above assumption is in fact equivalent to a part of the eigenstate thermalisation hypothesis (ETH) [2]: that in the thermodynamic limit L → ∞, averages of local observables in energy eigenstates vary continuously with the energy density.
One can get further with one more assumption: Assume that max n,m | n|A|m | 2 is small, and that the energy differences E m − E n are nondegenerate (this implies, but is stronger than, non-degeneracy of energies).
Then, von Neumann proved that ψ(t)|A|ψ(t) has small time variance, This is (a loose formulation of) von Neumann's quantum ergodic theorem.
The assumption that max n,m | n|A|m | 2 be small is implied by the full formulation of the ETH. Indeed, in the ETH, one equates off-diagonal matrix elements of local operators with an entropy factor that vanishes exponentially with the system size, times matrix elements from random matrix theory. Importantly, the assumption that max n,m | n|A|m | 2 be small is in fact weaker. The energy-difference non-degeneracy condition would be paralleled by the assumed random-matrix distribution of eigenvalues in quantum chaos; although the full connection is not obvious.
One sees immediately two potential drawbacks from von Neumann's ergodic theorem: • The assumptions on maximal deviations from the microcanonical averages are somewhat contrived; it is not clear, given a short-range interaction on a quantum lattice, how to show their validity. (The same can be said a fortiori of the ETH.) Is there something more universal that can be said? Are these assumptions really necessary to understand the structure of long-time dynamics in many-body quantum systems?
• Non-degeneracy assumptions are natural, but they are seemingly unrelated to the thermodynamic limit, where a continuum of states exists. Can something be said when there are degeneracies, most interestingly in the presence of a large number of nontrivial conserved quantities? Can this be formulated directly in the thermodynamic limit?
Remark 2.1. In von Neumann's work, instead of local observables, it is "macroscopic observables" that are the main focus. One should consider a family of macroscopic observables M i , with the property that they almost mutually commute, [M i , M j ] ≈ 0, and thus can be almost simultaneously diagonalised. The associated subspaces correspond to different "macrostates", intuitively different cells in the classical phase space. Von Neumann's theory can be formulated in terms of projections on these subspaces, and conditions on these. One may take spatial averages of local observables,Â := |Λ L | −1 x∈Λ L A(x), which indeed have vanishing commutators, ||[Â,B]|| ∼ |Λ L | −1 . We find that it is more transparent to forgo the abstract notion of macroscopic observables and their associated eigenspaces, and instead concentrate on local observables and their matrix elements.

Many-body and almost-everywhere ergodicity theorems
Recall that in von Neumann's quantum setup, we considered the "macroscopic", or many-body, limit L → ∞ in order to say something about ergodicity.
It is indeed commonly believed that, in typical interacting many-body systems and states · · · of physical importance, auto-correlation functions of local observables in the thermodynamic limit vanish at long times, A(0, t)B(0, 0) − A B → 0 (t → ∞). This is the notion of mixing 3 : the effect of the small localised perturbation B(0, 0) mixes with its surrounding and, locally, decays in time. A weaker notion, implied by mixing, is that of ergodicity: . This is equivalent to mean-square ergodicity: the vanishing of the variance of the observablē at long times, Ā 2 T − A 2 → 0. As explained above, mean-square ergodicity says that infinitetime averages of observables, in the state · · · , do not fluctuate, and are simply equal to ensemble averages.
But an important observation, which we believe was mostly overlooked in the literature, is that the resulting question of ergodicity, in fact, has less to do with quantum mechanics, and much more with the emergence of large-scale behaviours in many-body systems. This is about the passage from microscopic laws to macroscopic dynamics. As we mentioned, in manybody quantum systems, quantum recurrence, as discussed in paragraph 2.2.1, is irrelevant to this question. But also, in many-body classical systems, the uniform covering of the energy shell, as discussed in subsection 2.1, is also irrelevant to this question. In both cases, timescales are too large. The question of emergent ergodicity in many-body systems is of a different nature, and it brings out similar concepts both in quantum and classical many-body systems; "quantumness" does not play a fundamental role.
In this subsection, we take the C * algebra viewpoint on statistical mechanics. This, we argue, is much better adapted to the thermodynamic limit and the question of many-body ergodicity. We discuss this question, and in particular review results on almost-everywhere ergodicity established recently by the authors, which are extremely general and free from hardto-check conditions. For clarity we still concentrate on quantum spin lattices, however this viewpoint is much more easily adapted to other many-body systems, with only small differences between the quantum and classical cases.
Interestingly, as we will see, the basic mathematical ideas behind quantum recurrence (especially as "mean-square ergodicity") and classical ergodicity are still the right ones for many-body ergodicity, however only after the space of observables has been reduced to that relevant for local physics in the thermodynamic limit. It is the reduction to this space that allows these ideas to extract the correct form of ergodicity seen in many-body systems, associated to times that do not grow with the system size.
We note that the notions of ergodicity discussed in this subsection are not to be confused with what is sometimes called the "ergodic principle" (see e.g. [33]). This principle stipulates that relaxation must lead to a state of the Gibbs form. What we discuss here is about correlation functions of local observables, and is associated to conserved quantities that are "localised": time-invariant local observables. By contrast, thermalisation, the ergodic principle, and the Gibbs form, are concepts related to extensive conserved quantities. Extensive conserved quantities are involved in the hydrodynamic projection theorem proved in Section 3, and in the thermalisation theorems shown in [34], where it is proposed that the Hilbert space of extensive conserved quantities is in fact the tangent space to the manifold of maximal entropy states. Localised and extensive conserved quantities span the kernels of different evolution operators, acting on different spaces of observables (local observable, and extensive observables, repsectively).
We believe that analysing the structure of local observables in thermodynamic state is much more directly relevant to local physics than analysing the structure of spectra and how they behave as L → ∞. But, perhaps most importantly, in the thermodynamic limit, more can be said, giving almost-everywhere ergodicity, which ultimately is at the basis of the emergence of hydrodynamic structures.

C * algebra formulation and basic results
The strength of the C * algebra formulation is that we can directly study dynamics in the thermodynamic limit: the limit L → ∞ is already taken, and we concentrate on the remaining, relevant degrees of freedom.
We thus consider a quantum spin model on an infinite lattice, with local spaces C N on Z D (N, D ≥ 1). The Hamiltonian is homogeneous (space-translation invariant) and with shortrange interaction (exponentially decaying or faster, see Appendix A for the precise formulation); this generalises slightly the form (5), which was of finite-range interaction. Observables form a C * algebra U, the norm-completion of the algebra U loc of local observables (supported on finite numbers of sites). As above, we denote A(x, t), B(x, t), . . . observables A ∈ U translated to space-time points x, t, and sometimes omit the time argument if t = 0. A state 4 ω is a continuous, positive linear functional on U, interpreted as giving ensemble averages, normalised to ω(1) = 1. It is assumed homogeneous and stationary, Most importantly, the state is assumed to satisfy a property which we will simply refer to as "spatial mixing". We consider two different implementation: "factoriality", and a uniformenough spatial mixing, see Appendix A. Both include, in particular, the clustering of correlations at large separations, lim Results reviewed below on almost-everywhere ergodicity apply either if the state is factor ([35, p. 81, Def. 2.4.8]), or uniformly-enough spatially mixing (or both). Any thermal state that is in a single thermodynamic phase (i.e. Kubo-Martin-Schwinger (KMS) state satisfying (17), see below) is an example of a state with all above properties: for D = 1 this is every thermal state at nonzero temperature as there are no thermal phase transitions, for D > 1 this includes every thermal state above a certain model-dependent temperature [36,Section 5.3].
Fundamental theorems have been established that guarantee that the infinite-volume limit of finite-volume thermal states exist and gives a state ω on U; this is so at least at large enough temperatures if D > 1 (see e.g. [36,Section 5.3]), and at all temperatures if D = 1 [37]. For instance, for H of the form (5), there exists β c , with β c = ∞ if D = 1, such that for all β < β c and for all local observable A ∈ U loc , the limit exists and gives rise to a positive, continuous (with respect to the topology induced by the operator norm), normalised linear functional, that is invariant under space-time translations, The result is extended by continuity to U. The thermodynamic limit point in Equation (18) for all A, B in a dense (with respect to the norm topology) τ -invariant * -subalgebra of U τ , where U τ is the set of entire analytic elements for τ .  [37]. This unique state is factor, as an immediate consequence of [36,Theorem 5.3.30]. By the result of existence of space invariant KMS states [36, pg. 296], this unique state will also be space invariant. Note that time invariance is an immediate consequence of the KMS condition, Equation (20).
Further, it is guaranteed, again for β < β c , that the state ω from Eq. (18) is exponentially clustering: there exists γ > 0 such that for every local A, B ∈ U loc , there is c > 0 such that (the distance can be taken as the L 1 or L 2 distance -as chosen for convenience -on the square lattice between the supports of A and B). In particular, ω satisfies the spatial mixing condition. More generally, any Kubo-Martin-Schwinger (KMS) state satisfying (17) satisfies also (19) and the spatial mixing condition. Finally, by the Lieb-Robinson bound [38], this exponential decay can be extended in spacetime [18], to the full region outside of a Lieb-Robinson "light-cone": there exists v LR > 0 and γ > 0 such that for every local A, B and v > v LR , there is c > 0 such that The structure in space-time is depicted in Figure 2. We note that in the C * algebra formulation, the concept of microcanonical shell, Eq. (6), does not make immediate sense, especially in view of the ensemble equivalence (7). Nevertheless, for the purposes of interpretation of the results, it will be convenient to consider any state ψ that is "near to" ω as being part of the microcanonical shell in the thermodynamic limit, ψ ∈ S micro . Here, "near to" simply means "absolutely continuous" with respect to ω, that is, of the form These states are still positive, continuous, normalised linear functionals, but not space-time translation invariant (which is of course as expected for states in the microcanonical shell). Further, by spatial mixing, the density of any extensive observable in ψ is the same as it is in ω: In particular, the energy density is the same in ψ as in ω.
Finally, in the theory of C * algebras, a very important role is played by the Gelfand-Naimark-Segal (GNS) construction. In this construction, one associates to any state ω a Hilbert space H ω , constructed from the observables in U, with in particular the vector |Ω associated with the identity operator 1. This Hilbert space is such that ω(A) = Ω|A|Ω , and more generally ω(B * AC) = B|A|C [35, Chapter 2.3] (here we use the same notation for A ∈ U and its representation as a bounded operator on H ω ). More precisely, the Hilbert space is that arising from the completion of U/N with respect to the inner product induced by the sesquilinear form (A, B) → ω(A * B) on the equivalence classes A + N ∈ U/N, with null space N (satisfying ω(B * N) = 0 ∀ B ∈ U). Operators acting on H ω are sometimes called "superoperators" in the condensed matter literature. For ω being the infinite-temperature state or trace state (see above), the resulting norm is in fact (the infinite-dimensional completion of) the Hilbert-Schmidt norm, but the GNS construction is more general. Physically, in thermal states, H ω is the space of particle and hole excitations above the thermal "vacuum".
Space and time translations can be extended to unitary group actions on the GNS representation, see also [35, Theorem 2.3.16, Corollary 2.3.17]. Here we express the related theorem; below we will use the same notation τ t and ι x for time and space translation unitaries on H ω . Theorem 2.1 (GNS representation). Given a state ω ∈ E U of a unital C * algebra U there exists a (unique, up to unitary equivalence) triple (H ω , π ω , Ω ω ) where H ω is a Hilbert space with inner product ·, · , π ω is a representation of the C * algebra by bounded operators acting on H ω and Ω ω is a cyclic vector for π ω , i.e. the span{π ω (A)Ω ω : A ∈ U} is dense in H ω , such that If additionally we have a group G of automorphisms {τ g } g∈G of U and ω is τ -invariant, then there exists a representation of G by unitary operators U ω (g) acting on H ω . This representation is uniquely determined by and invariance of the cyclic vector See Appendix A for the full mathematical setup.

Many-body ergodicity
In order to look at ergodicity, in the classical realm we needed to consider the energy shell to which the initial condition belonged, and in the von Neumann quantum setup, we needed a small quantum energy shell. So, in the C * algebra context, it is natural to take a state ω assumed spacetime stationary and spatially mixing as described above, and states ψ in its microcanonical shell, as in (23).
In the C * -operator algebra context, one can still consider time-averaging (1) of operators (either as elements of the C * algebra, or bounded operators on H ω ), where time integration is now to be understood as a Bochner integral [39]. For clarity, we will write explicitly the long-time limits, hence we will use (16).
The basic formulation of many-body ergodicity is conceptually similar to that of subsection 2.1 (classical case), and paragraph 2.2.1 (quantum case): Theorem 2.2. (Mean-square ergodicity.) Let ω be spacetime stationary and spatially mixing. Suppose that the unitary time evolution operator τ t on the GNS space H ω has trivial kernel, ker(τ s − 1) = C|Ω for some s ∈ (0, ∞). Then for all integers n ≥ 1, That is, we have mean-square ergodicity, and in fact the "random variable" A is non-fluctuating within the state ω.
The theorem follows from von Neumann's ergodic theorem for unitary operators on Hilbert spaces, as in paragraph 2.2.1, but now with the GNS space H ω as the Hilbert space. Using von Neumann's ergodic theorem the time average gives a projection onto the subspce of H ω that is invariant under τ s , for some s ∈ (0, ∞), see also [40,Theorem 3.1]. If then ker(τ s − 1) = C|Ω , this projection is the rank one projection onto |Ω , and the result is straightforward, as a direct application of the techniques used in [18] (also discussed in [41,Section 2]). Relevant results exist in various forms, including [42, Proposition 6.3.5], [43], [44].
Clearly, as we have taken already the limit L → ∞, the timescales for convergence are not controlled by the volume of the system, contrary to the phenomenon of quantum recurrence as discussed in paragraph 2.2.1; in (28) convergence of long-time averages happen much before the recurrence time. Yet, the techniques are essentially the same. What has happened? The most important aspect, as mentioned, is that in the C * algebra formulation, we have reduced the set of degrees of freedom to the relevant ones for local physics. Thus, in contrast to the condition of non-degeneracy, the condition of triviality of the kernel of the time evolution properly encodes many-body ergodicity.
In a formal sense, this kernel triviality condition is akin to asking that there be no nontrivial constants of the motion in the classical case, or to asking that there be no degeneracies in the quantum case. However, in a more physical sense, it is starkly different. In the classical context, for instance, integrability implies that there are nontrivial constants of the motion (as many as degrees of freedom). But in extensive systems, integrability is not related to the presence of nontrivial elements in ker(τ t − 1). Instead, it is related to the presence of extensive conserved quantities in addition to the Hamiltonian. These span the kernel of a different evolution operator, that which acts not on the GNS space built from the local observables, but on a Hilbert space of extensive observables, see section 3. Extensive conserved quantities may exist, the system may be integrable, and yet many-body ergodicity may still hold. Thus, we have separated the condition for many-body ergodicity to hold, from the requirement on nontrivial extensive conserved quantities, such as those found in integrable systems.
The "nontrivial constants of the motion" that should not be present for many-body ergodicity to hold, are the time-invariant local observables (locality is here in the weak sense of elements of the C * algebra). Typically, except in special systems such as those with dynamical constraints, there are no such observables beyond (the trivial) multiples of 1. In fact, we note that timeinvariant local observables play a rôle in the phenomenon of localisation -thus many-body ergodicity occurs if there is no "localisation" (in the sense that there is no time-invariant local observables). As mentioned, extensive conserved quantities, instead, play a rôle in integrability. And in fact, most importantly for us, they are at the basis of the emergent hydrodynamic structures, which we explain in Section 3.
It turns out that one can go further if the state ω is also a KMS state, Eq. 20. In fact, it is not necessary for H in (18) to be the evolution Hamiltonian, or equivalently for τ t in (20) to be the dynamics we take for time averaging. Any other H , any other dynamics τ t , that keeps the state spacetime stationary will do -for instance, any other extensive conserved quantity, if available. Thus, in integrable model, any generalised Gibbs ensemble [45,46,47].
Again, the proof is omitted, but follows quite directly from the techniques developed in [18]. Recall that C T → 0 in the strong operator topology means ||C T |ψ || → 0 for all |ψ ∈ H ω . Now the result is much stronger: this is truly an ergodicity property, stating that the time average of any observable, in any state within a microcanonical shell surrounding ω, gives the average in ω. Time averaging gives back ensemble averaging. The operatorial equation, for A as an operator on the GNS space H ω , is in fact an equivalent statement.
Finally, an additional insight may be gained into many-body ergodicity by the following construction. Recall that in the classical case, ergodicity implies that every point in the energy shell is approached arbitrarily closely by the trajectory as time passes. Here, instead, with states playing the role of "points in phase space", a trajectory starting from ψ approaches arbitrarily closely the final state ω: for every > 0, there is a finite time t such that the timeevolved state ψ • τ t is a distance less than from ω, and in fact, this is true for a fraction of time that tends to 100% over long periods. This construction starts with the choice of a sequence of unit-norm, strictly local observables A n ∈ U loc , ||A n || = 1, n = 1, 2, 3, . . . which spans U loc (every element of U loc may be written as a finite linear combination of these). Such a sequence certainly exists. Then, we define the following distance function between states: We note that there is a lot of freedom in defining the distance function; in particular, any finite set of C * -algebra elements A n ∈ U, n = 1, 2, . . . , N may be adjoined to the sequence A n . Thus, the general result on the distance function can be used to deduce results for any finite-dimensional subspace of such elements.
As d(ψ, ψ ) = 0, then ψ(A n ) = ψ (A n ) for all n. As A n 's span U loc , then also ψ(A) = ψ (A). Therefore by the triangle inequality, |ψ(B) − ψ (B)| ≤ 2 . As this holds for every > 0, we and we are interested in the long time limit of the time average: It is easy to see, by the dominated convergence theorem, that we can move the lim T →∞ past the summation ∞ n=1 . We then want to apply Theorem 2.3. We can immediately do this for the second term: For the first term, 1 T T 0 ψ 2 (τ t A n ) dt , we notice that its limit coincides with the limit of ω(A) 1 T T 0 ψ(τ t A n ) dt, by the following argument: (35) and again by applying Theorem 2.3. Finally, applying the limits in Equation (33), we show: This implies the following. Proof. We prove the last statement as it implies for first. Let us denote µ + T and µ − T the measures for the sets {t ∈ [0, T ] : d(ψ t , ω) ≥ } and {t ∈ [0, T ] : d(ψ t , ω) < }, respectively. Note that , hence by (37) we must have In particular, from the above theorem we conclude that, for any A ∈ U, we have |ψ(τ t A) − ω(A)| < for 100% of the times: the average in a state in a microcanonical shell, is almost always as near as desired to the average in the stationary state ω.
Some remarks are in order: 1. Because the results apply to the C * algebra obtained from local observables, the results address local relaxation, a form of typicality (see e.g. [7]) which expresses many-body ergodicity.
2. Note how both Theorems 2.2 and 2.3 apply not just to thermal states, but to more general spacetime stationary states, such as generalised Gibbs ensembles. This therefore takes into account the possibility that the model possesses a large number of extensive conserved quantities similar to the Hamiltonian; these would lead, in the finite-volume setup, to a large number of degeneracies. Here, no condition arises on the presence or not of extensive conserved quantities.
3. The microcanonical shell represented by the states ψ should be understood as a being a shell around the microcanonical state with respect to all extensive conserved quantities, not just the energy: indeed all extensive variables have densities that are the same as in ω, Eq. (24). This is why ergodicity can hold in ψ without conditions on the presence or not of extensive conserved quantities.
4. Requiring the KMS condition is not too strong a requirement. It is expected that most, or maybe all, spacetime stationary spatially mixing states are in fact KMS states. Perhaps the Tomita-Takesaki theory [36] could be used.

5.
A similar C * algebra construction can be made for classical systems, such as a classical gas, where the limit of infinitely many particles is taken. We expect that similar ergodicity statements can be obtained. In the classical case, the GNS space can also be constructed, and time evolution is again -at least in well-behaved systems -a one-parameter unitary group. Thus, von Neumann's ergodic theorem for unitary operators applies. Again, as we have taken the limit of infinitely-many particles, timescales are not those that would be necessary to cover the full energy shell -it is not true that over long times, the many-particle trajectory comes as close (in the conventional many-particle metric) as desired to a given point in the many-particle phase space. Particle numbers are so large that the phase space is immense, and "long times" are never enough.
Thus we have solved the problem of obtaining ergodicity in many-body systems with timescales below the quantum recurrence timescale (in a different way than in von Neumann's QET, with different-looking conditions), or below the classical shell-covering timescales. We have also solved the problem of generalising to cases where nontrivial conservation laws exist (the extensive conserved quantities), and of having a formulation that is directly in the thermodynamic limit.
There remains the problem that the condition of triviality of the kernel of the time-evolution unitary one-parameter group on H ω is rather difficult to verify. It is also perhaps not too clear if the resulting ergodicity statement is the one that is most physically relevant. These two problems are solved by almost-everywhere ergodicity, to which we now turn.

Almost-everywhere ergodicity
We now obtain ergodicity results which are valid in every quantum spin lattices as described above, and in every spacetime stationary spatially mixing states. No special, hard-to-check condition is required.
Such general results however are not valid for the standard notion of ergodicity, with integration purely in the time direction. Instead, we must look in spacetime. On the lattice, it is natural to restrict to the set of rational spatial directions S D−1 Q = {x/|x| : x ∈ Z D } (this is dense on the sphere). Given T > 0, v ∈ R and n ∈ S D−1 Q , we consider the average of A ∈ U along the ray with velocity vector v = vn, where a = ( a i ) i is the vector of the integer parts.
In fact, the results apply more generally. That is, instead of considering the flat averages, we may modulate the average with oscillating factors. So, we consider instead (keeping the wavenumber k and frequency f implicit in the notation for the average, for lightness of notation). Let us also denote See Fig. 1 for every A ∈ U, m ∈ N. If, further, ω is a KMS state, then: (Ergodicity) for every ψ ∈ S micro , and for every A, B ∈ U, in the strong operator topology.
Additionally, under the same assumptions of the Theorem above, and similarly to the derivation around Equation (33), one can show that for any ψ ∈ S micro in the microcanonical shell of ω, the ray-averaged distance (with respect to the metric of Equation (30)) between ψ and ω tends to zero, for almost every ray. Denoting ψ υt := ψ • ι υt τ t , we have: Taking the language of probability, Eq. (41) means "Ā v T → E A in law". This is almosteverywhere ergodicity: along almost every velocity,Ā v T tends to a non-fluctuating quantity. As a consequence, by the Cauchy-Schwartz inequality, Eqs. (42) are starker statements. They hold in spatially mixing, space-time translation invariant KMS states, including Gibbs and generalised Gibbs states. They imply that the limits in (42) hold within any expectation value, multiplied by any other observables in any order. We thus obtain an extension of (a weak version of) the Lieb-Robinson bound to within the Lieb-Robinson cone: the ray-averaged operator A v T becomes "thin" as T → ∞, being un-observable by any B, at every frequency and wavelength, and at almost every velocity v Further, Theorem 2.5 means that time-averages, at every frequency and wavelength and almost every velocity v, converge to non-fluctuating (cluster out), classical (commutators vanish) variables.
In particular, we have [Ā v T , B] 2 → 0. At short times, the OTOC [A( vt , t), B] 2 is expected to grow exponentially in chaotic systems [10,48] (eventually reaching O(1) values by boundedness of the state). The velocities where this happens are bounded by a statedependent "butterfly velocity" v B , |v| < v B < v LR . We find that the time-averaged version of this quantity, , t ), B] , in fact decays for almost every velocity.
The results do not say anything about what happens in the pure time direction, v = 0. For instance, there are non-interacting Hamiltonians for which the velocity v = 0 is not ergodic; as mentioned, this is in fact related to localisation. However, all results hold for velocities as near as desired to 0.
Again, Theorem 2.5 is based on the von Neumann ergodic theorem [32, Theorem II.11], which relates time averages to projections onto invariant spaces. As far as we know, before our work there were no nontrivial results concerning the invariant subspace of H ω for unitary time evolution in quantum lattice models. We show it to be spanned by |Ω at almost every velocity. This happens because finite-dimensional local spaces (in fact, we believe that countabledimensionality is sufficient) are too small to allow for operators to extend thickly on large distances over time: the set of speeds where the operator remains supported under averaging should be "small enough".
Remarks 1 -5 of paragraph 2.3.2 apply to almost-everywhere ergodicity results as well. We emphasise again that the results of Theorem 2.5 do not require any nontrivial condition on the interaction or the dynamics (such as chaos, non-degenerate eigenvalues, absence of nontrivial constants of the motion, non-localisation condition, etc.). The results are purely consequences of extensivity of the system -that is, the fact that the interaction range is short enough, and that the thermodynamic, large-volume limit has been taken.
Thus we have now addressed the question of the conditions under which notions of ergodicity may hold: almost-everywhere ergodicity only requires extensivity. However, we have not yet addressed the question of the physical relevance of this notion. Surprisingly, it turns out that almost-everywhere ergodicity is in fact closely related to the hydrodynamic structure emerging on large scales. We now turn to this important result.

Hydrodynamic projections and spaces of extensive charges
The hydrodynamic projection theorem says that the Euler-scale connected correlations decompose into the conserved charges of the system. For its formulation, we will construct Hilbert spaces of extensive quantities: intuitively to every local observable A ∈ U loc (density) we will associate an element, formally, ΣA = x∈Z D A(x), defined as an appropriate equivalence class 5 . The Cauchy-completion of the space of such elements will give a Hilbert space of extensive quantities H 0 . Under appropriate clustering assumptions, we can define the time evolution as a unitary action on the Hilbert spaces of extensive quantities, and thus have an unambiguous definition of the conserved extensive charges Q 0 . This will be a closed subspace of H 0 , thus allowing us to define the orthogonal projection P : H 0 → Q 0 .
The objects of interest are the time-averaged, long-wavelength Fourier transform of two-point connected correlations: where is the sesquilinear connected correlation. We do not know how to show the existence of the limit lim T →∞ , but we can show that x∈Z D e iκ·x/t (A(x, t), B) is uniformly bounded. This allows us to use the notion of a Banach limit lim T →∞ , see [49,Chapter III.7] and [29, Appendix A]. We can also show that the result is indeed a function not merely of A, B, but in fact of the equivalence classes ΣA, ΣB, as the notation suggests; and that it is continuous, on both variables, with respect to the norm || · || 0 on H 0 . The essence of the hydrodynamic projection is that inside correlation functions every extensive observable will project to a conserved charge, irrespectively of the choice of a Banach limit: This is an expression of the Boltzmann-Gibbs principle [15,16,17]: the reduction of the number of degrees of freedom at large space-time separations by projections over hydrodynamic modes. The physical idea is that the initial, dynamically complicated disturbance quickly relaxes and projects, at the Euler scale of long time and large distances, onto the extensive conserved quantities, that then carry correlations. This occurs thanks almost-everywhere ergodicity, as is made apparent in the proof of Theorem 3.1 in Section 3.4 (see [29] for proof in D=1). We believe this is the first general, rigorous result concerning the Boltzmann-Gibbs principle in deterministic interacting systems of arbitrary dimensions. The principle is applicable to the large class of quantum lattices, going beyond interacting particle systems conventionally studied in statistical physics and hydrodynamics. It provides further support to the idea that the basic principles of hydrodynamics hold independently from the details of the microscopic dynamics.
Surprisingly, the hydrodynamic projection can be generalised to describe oscillatory behaviours: to any frequency f ∈ R, and wavenumber k ∈ R D . We show that the projection still takes place for k-extensive quantities, formally Σ k A = e −ik·x A(x), which will project to (f, k)-conserved charges, defined as those Σ k A such that their time evolution take the form Oscillatory hydrodynamic projection describes oscillatory behaviours that emerge at large space and time separations in correlation functions. One in general expects that (f, k)-conserved charges only exist for certain frequency-wavenumber pairs (f, k), depending on the specific model. A simple example is the free fermionic lattice, with some dispersion relation E(k).
In this system, creation and annihilation operators at momentum k are (E(k), k)-extensive conserved quantities. In Section 4, in order to illustrate the phenomemon we will explain (in a non-rigorous fashion) how oscillatory hydrodynamic projection recovers the oscillatory algebraic decay of fermion two-point functions in a free fermionic quantum chain, in agreement with a saddle-point analysis.
In the present section, we concentrate on the rigorous and general results of hydrodynamic projections in quantum lattice models.

Clustering and basic assumptions
Throughout all of Section 3 we consider a dynamical system (U, ι, τ ), see Definition A.1 and Appendix A for more details. We use the notation A(x, t) for space-time translated observables ι x τ t A. We also denote A(x) = A(x, 0), while we keep the notation τ t A for pure time translations. The system will be in a space-time invariant state ω, Definition A.2, which will be assumed to have clustering properties with respect to space translations. The almost-everywhere ergodicity theorems shown in [18], and discussed in Section 2.3, require the state ω to be clustering in space, in particular lim x→∞ ω(A(x)B) − ω(A)ω(B) = 0 for any A, B ∈ U. Almost-everywhere ergodicity holds no matter how fast the connected correlations of observables decay at large spatial separations. For the hydrodynamic projection to occur stronger clustering assumptions need to be imposed. In particular, clustering of correlations at large space separations will need to happen faster than |x| p , for p large enough: Additionally we will need to control the dependence on the size of the supports of the observables in the clustering bound; the support of a local observable A is defined as Λ A := {X ⊂ Z D : A ∈ U X }. That is because time evolving observables "spreads out" their support. A state satisfying this form of clustering will be called r-sizeably p-clustering: and C A,B is bounded as: A subset C ⊂ U is called uniformly p-clustering for p > D if there exist a C > 0 s.t. for all A, B ∈ C: Of course C A,B will, in general, depend on the pair of observables, i.e. on their norm and the size of their support. In D = 1 quantum spin chains, [29], the state is assumed to be r-sizeably p-clustering; this notion originates from [34,Definition 4.2]. The definition of p-clustering, for every p > 0, includes high temperature KMS states, which are exponentially clustering in space, as shown in [50,51]. and Uniform (exponential) clustering of a subset C ⊂ U is defined similarly.
In order for the time evolution to be well defined as a unitary group action on the Hilbert spaces of extensive quantities, we need to assume that the time translations of local observables cluster in a uniform enough manner. First, we require that any element A in span{τ t A : A ∈ U loc , t ∈ R} is approximated by a sequence σ n A of local elements. We will show that this follows from the Lieb-Robinson bound. Additionally, we require that for any local A, B the set of pairs {(σ n A, σ n B)} is uniformly clustering in space, as per Equation (49). Combining spatial clustering of the state ω with the Lieb-Robinson bound, one can obtain space-like p c clustering, defined as follows: Definition 3.3. The dynamical system (U, ι, τ ) in the state ω is called space-like p c -clustering with velocity υ c , if 1. ∀A ∈Û loc = span{τ t A : A ∈ U loc , t ∈ R} there exists a sequence σ n A ∈ U loc , n ∈ N, such that lim σ n A = A. For any A ∈ U loc we define σ n A = A, ∀n.
2. ∀A, B ∈Û loc the set of pairs {(σ n A, σ n B)} is uniformly p-clustering for some p > p c .
3. ∀A, B ∈ U loc there exist p > p c , 0 < V < υ c and C A,B > 0 such that for all x ∈ Z D , |x| ≥ V |t|.
Similarly, we define exponential space-like clustering by the same conditions, with Condition 2 replaced by uniform exponential clustering and Equation Space-like p c -clustering for p c > D, where D is the lattice dimension, is a sufficient condition for showing the hydrodynamic projection theorem. In turn, for space-like p c -clustering to hold, it is sufficient to have the following two conditions:  [29,Section 8] to arbitrary lattice dimension D, to show that whenever p > D(r + 1), a r-sizeably p-clustering state will be space-like p c -clustering for p c > D, in any system with interaction that satisfies Equation (127). This is slightly more technical in D > 1, but it largely follows the same ideas as in the D = 1 proof in [29], hence we omit it. Instead, in Appendix C we show how exponential space clustering and the Lieb-Robinson bound lead to exponential space-like clustering.

Spaces of extensive charges
We start by constructing Hilbert spaces of extensive quantities from the dynamical system (U, ι, τ ). For each wavenumber k ∈ R D we define the positive-semidefinite sesquilinear form where we recall that (A, B) = ω(A † B) − ω(A † )ω(B) is the sesquilinear connected correlation. We define the equivalence relation A ∼ k A on U loc by A − A, A − A k = 0. The Hilbert spaces of extensive quantities are the norm-completion of the quotient spaces formed of the set of equivalence classes of ∼ k . The equivalence class of A ∈ U loc is Σ k A := {A ∈ U loc : A ∼ k A} and it is to be understood as the k-extensive observable associated to the density A. Formally, we can associate it with the infinite series Σ k A = x∈Z D e −ik·x A(x), the "total" A of the full quantum lattice. Of course this series does not converge within U, but it does converge weakly in the GNS representation associated to the state: this weak limit coincides with the above definition of Σ k A. The set of equivalence classes is V k := U loc / ∼ k = {Σ k A : A ∈ U loc } and its Cauchy completion gives the Hilbert space of k-extensive observables H k = V k . In order to examine the long time dynamics of these extensive observables we have to extend the action of time evolution τ t , t ∈ R to the Hilbert spaces H k . This is done rigorously in [29, Section 5.3] for D = 1 (spin chains) and can be immediately extended to arbitrary dimension D. The key in order to be able to do this is the uniform clustering condition 2 in Definition 3.3. Thus, we can show that time evolution acts as a unitary operator τ k t on H k . For lightness of notation, we omit the superscript and simply write τ t for the unitary action of time evolution on these spaces.
With this construction, we can define the subspace of conserved extensive charges as those elements of H k that are invariant under τ t : In fact, we can go one step further and define the subspace of (f, k)-extensive conserved quantities (or f -oscillatory k-extensive charges) as: See Section 4 for a discussion of these charges and examples in a free fermion chain. It is immediate that Q (f,k) is a closed subspace of H k , hence there is an orthogonal projection

A hydrodynamic projection theorem
The correlation functions with (f, k)-fluid-cell averaging are given by Note how one extracts, thanks to the time integral and factor e −if t , the time-oscillatory behaviour of the correlation function with frequency f . Note also how the full wavenumber is k + κ/t, representing, in the large-time limit, a long-wavelength modulation of a k-oscillatory factor; this extracts the space-oscillatory behaviour with wavenumber k.
The limit on T is in general a Banach limit [49, Chapter III.7], [29, Appendix A]. The result will hold irrespectively of the choice of Banach limit, and all the aspects of the proof regarding the Banach limit are the same as in D = 1, as shown in [29, Section 6]. As the notation implies, the result only depends on the equivalence classes Σ k A, Σ k B, which we'll denote with the respective lowercase letters a, b. In fact this defines a continuous sesquilinear form on H k . The hydrodynamic projection theorem is rigorously stated as follows: For every frequency-wavenumber pair (f, k) ∈ R×R D , rational vector κ ∈ R D : ∃r ∈ R | rk ∈ Z D , and any k-extensive elements a, b ∈ H k , Specifically for the simple case (f, k) = (0, 0) we have The case κ = 0 is of special interest, and is known at the (oscillatory version of the) Drude weight: Here the limit is in fact an ordinary limit. The projection formula, when written explicitly in a basis decomposition, is exactly (an oscillatory version of) the so-called Mazur bound, here shown rigorously to be saturated. This is a very immediate extension to the oscillatory realm, and to arbitrary dimensions, of the result [29, Thm 6.1], and can be proven quite straightforwardly without the extensive machinery we develop below for dealing with the case κ = 0, see the proof of [29, Thm 6.1] and the discussion there for more details.
We proceed with the proof of the hydrodynamic projection theorem. We will first prove it for (f, k) = (0, 0), and then easily generalise the proof to arbitrary (f, k). We note again that for κ = 0, the theorem is a straightforward application of von Neumann's ergodic theorem for unitary operators. We thus concentrate on the case κ = 0.

Proof of hydrodynamic projection
The main idea for the proof of the hydrodynamic projection formula in quantum spin lattices is this: By using an appropriate geometric construction, we recast the D-dimensional problem of hydrodynamic projection into the a 1-dimensional problem. This is done by identifying the summation over space coordinates in a plane perpendicular to the wavenumber direction, with a sesquilinear form and its new associated Hilbert space, which is to play the role of the sesquilinear correlation in an effective one-dimensional problem. Once this is done, the proof follows that of [29] done for D = 1. The schematics of our proof is as follows: We assume a dynamical system in a p-clustering state ω, Definition 3.1 and with sufficiently fast decaying interactions, Equation (127), so that the Lieb-Robinson bounds, Equation (129), holds. We combine p-clustering with the Lieb-Robinson bound (proof in Appendix C) in order to obtain space-like p c -clustering, Definition 3.3. This leads to the special Property 3.3, below, for the dynamical system. Using this property we can prove Lemma 3.4, which forms the basis for the hydrodynamic projection formula, Theorem 3.1.
The most important difference from the D = 1 case is the new geometric constructon, which we now explain. We define the rational unit sphere in D-dimensional space as S D−1 Q = {κ/|κ| : κ ∈ Z D }. We also denote the set of vectors in rational directions in R D \ {0} by For any such vector κ ∈ R D Q , there is a unique unit vector κ = ±κ/|κ| ∈ S D−1 Q that lies on the half-unit sphere: such that κ 1 > 0, or κ 1 = 0, κ 2 > 0, or · · · , or κ 1 = 0, . . . , κ D−1 = 0, κ D > 0. Then, given such a unit vector κ , we denote by κ = r κ the "smallest" integer vector associated to it, where r > 0 is the smallest positive number such that κ ∈ Z D ( κ exists because κ ∈ R D Q ). Clearly rκ = κ for every r ∈ R, and in fact, κ is an element of Z D that identifies the "rational ray" on which κ lies: there is a bijection between the set of κ's and the set of {rκ : r ∈ R}'s for κ ∈ R D Q . Below we will refer to κ as a rational ray, and it will be important that it lies in Z D . For any rational ray κ, we may build its perpendicular plane in Z D : a sublattice of Z D that is isomorphic to Z D−1 , and whose vectors have vanishing vector dot product with κ. Given κ, there is a set (possibly empty) of zero components ζ = {i : κ i = 0} and its complement (always non-empty)ζ = {1, . . . , D} \ ζ = {j 1 , . . . , j d } for unique 1 ≤ d = |ζ| ≤ D and j 1 < . . . < j d . We construct the vectors (with implied summation over repeated indices). Of course, we could have taken away the greatest common divisor of the components of h (i) in order to have "denser" planes, but this is not necessary in the following construction. By stacking the perpendicular planes parallely, we obtain a sublattice of Z D which is itself isomorphic to Z D , defined as the Z-module freely generated by the set {h (i) , κ}: In-between the points of this lattice lie "foundamental cells". That is, consider the region R κ = {x i h (i) + z κ : x i ∈ [0, 1)∀i, z ∈ [0, 1)} ⊂ R D and the cell in Z D given by Then the union of all unit cells shifted by the sublattice gives back Z D : Indeed, by linear independence and the fact that the Jacobian J is finite and nonzero, every u ∈ R D may be written in a unique way as u = x i h (i) + z κ for some (x, z) ∈ R D , and thus u = x i h (i) + z κ + λ where λ ∈ R κ . This holds in particular for u ∈ Z D , and since x i h (i) + z κ ∈ Z D κ it gives (66). We define for any wavenumber κ ∈ R D Q and κ ∈ R D , some positive-semidefinite sesquiliniear forms, and associated Hilbert spaces, in a way similar to that done above: first by summation over the perpendicular planes H D−1 κ , and second by summation over the full sublattices Z D κ := Note that these are well defined by clustering of the connected correlator (A, Similarly to the construction of the spaces H κ , we define the Hilbert spaces H κ,⊥ and H κ κ as the norm completions of V κ,⊥ := U loc / ∼ κ,⊥ and V κ The equivalence class of A ∈ U loc is denoted by the respective lowercase a and it is immediate that (·, ·) κ,⊥ and (·, ·) κ κ are inner products on their respective Hilbert spaces, and thus satisfies the Cauchy-Schwarz inequality. See Appendix D for proofs of basic properties. It can be established that ι x := ι x κ (x ∈ Z) and τ t (t ∈ R) act unitarily on H κ,⊥ as representations of the groups Z and R respectively, and that τ t (t ∈ R) act unitarily on H κ κ as a representation of the group R. Details and other basic results are discussed in Appendix E. We denote space-time translations of elements a ∈ Hκ ,⊥ by ι z τ t a, time translations of a ∈ Hκ κ by τ t a, while keeping the notation A(x, t) for A ∈ U. Property 3.3. Consider a dynamical system (U, ι, τ ) (Definition A.1) that is space-like D + 1clustering (Definition 3.3). For every κ ∈ R D Q , it follows that: 2. ∀ a, b ∈ V κ,⊥ , there exist T > 0 (independent of κ) and a Lebesgue measurable function f : R → R + such that: Proof. Property 3.3.1 is a consequence of the almost-everywhere ergodicity theorem, [18, Theorem 3.2]. One can see this by proving that (ι vt τ t a, b) κ,⊥ is space-like ergodic, [18, Definition 3.1], see Appendix B for the general idea. That is, first, note that our assumptions immediately imply that the state ω is space-like ergodic: Consider a, b ∈ Vκ ,⊥ and class representatives A, B ∈ U loc respectively. Then, Using space-like clustering we can see for any z ∈ R, υ > υ c and A, B ∈ U loc where the summand is N -independent, x-summable, and uniformly bounded by a m-summable function, allowing us to apply the dominated convergence theorem in order to move lim N inside the sum x∈H D−1 κ in Eq. (71). Thus, Applying the almost-everywhere ergodicity theorem for the system (H κ,⊥ , ι, τ ) (see Appendix B) we get Property 1. First consider the case |v| ≥ V /| κ| + T −1 . Note that (|v|t − 1)| κ| ≥ V t and (|v|t − 1)| κ| ≥ (|v|T − 1)| κ| ≥ V T . We have, for some p > D + 1 (each step is explained below): Note that the right-hand side never gets infinite, sinceκ is fixed and (|v|t − 1)| κ| ≥ V T > 0. In the second line we used the fact that which allowed us to use space-like p-clustering. We now use the fact that x∈Z D−1 F (x) = d D−1 x F ( x ), hence we'll have to deal with x i h (i) . The idea is to do the following transformation: The finite Jacobian of the change of variable above is (nonzero) We have by the triangle inequality: where we used x i = x i − ε i for some ε i ∈ [0, 1), and where h = i |h (i) |. Hence, doing the transformation x → y, defined in Eq. (76): In the third line we used the fact that (|v|T − 1)| κ| > 0. In the final line, we used (|v|T − 1)| κ| ≥ V T and, for every ∈ (0, h/V T ], where I, as defined by the right-hand side of the inequality, is finite, as p > D + 1 > D − 1, and only depends on h/(V T ) and D.
For |v| < V /| κ| + T −1 , we instead use the fact that by the Cauchy-Schwartz inequality and space-time translation invariance.
Thus we set We see that, for p > D + 1, this indeed satisfies the right properties, and in particular this lower bound on p is necessary for the integral of (|v| + 1)f (v) to exists on R.
From this, we obtain the following crucial lemma: Lemma 3.4. If the dynamical system satisfies Property 3.3, then ∀A, B ∈ U loc , ∀s ∈ R, there exists T 0 > 0 such that for every κ ∈ R D Q , the following holds: Proof. First we write Let us first simplify G s (t) for arbitrary s. Writing (in a unique fashion) x = r + λ where r ∈ Z D κ and λ ∈ Λ κ , we get where (recall that Λ κ is finite)Ã and the lowercaseã u , b ∈ V k κ/(t+s) are the respective equivalence classes ofÃ u , B (see Equation (67) and discussion below). We note that where we used the fact that lim κ →0 ||c|| κ κ = ||c|| κ 0 exists for any c. Therefore, using the form At this point, we have managed to recast the problem into an effectively one-dimensional problem. The proof now broadly follows that in D = 1, [29]. In order to show the Lemma, i.e show lim T →∞ 1 T T T 0 dtg(t) = 0, we want to commute the limit and the time integral past the summation z∈Z in g(t). The result would then follow by applying Property 3.3.1. We proceed to uniformly (in t) bound the summand in g(t). As in the proof of Lemma 6.6 in [29], we write (denoting κ = |κ||κ|): where υ t = t −1 υt . Having established Property 3.3, the rest of the proof follows exactly as in [29,Lemma 6.6].
Finally, having shown Lemma 83, we have for all A, B ∈ U loc , any rational vector κ ∈ R D , s ∈ R and any choice of Banach Limit lim t→∞ : where a, b ∈ V κ/t are the respective equivalence classes (see Equation (53) and below). We then proceed exactly as in the proof [29, Theorem 6.7], which completes the proof of the Hydrodynamic projection theorem for any dimension D and for zero wavenumber and frequency.
The proofs can easily be generalised to all frequencies and wavelengths, as all the bounds remain essentially unchanged. We consider different representations of the groups of space and time translations on the C * -algebra U, for each frequency-wavenumber pair, defined as (91) These form representations of the groups Z D , R respectively, by bijective linear maps on U. They are not * -automorphisms like ι x and τ t . However, linearity and unitarity on the Hilbert spaces defined from the state are all that is used in our proofs (and of course the group representation properties). Hence we can relax the Definition of a dynamical system, Definition A.1, to a triplet (U,ι,τ ) whereτ is a strongly continuous representation of the group R by bijective linear maps {τ t : U ∼ − → U} t∈R , andι is a representation of the translation group Z D by bijective linear maps {ι x : U ∼ − → U} x∈Z D , and the state is required to have the invariance property (ι xτt (A),ι xτt (B)) = (A, B). It is clear that if uniform p-clustering holds for a subset of elements in U loc under ι x , then it also holds underι x . Under this new represenations of space translations we can construct the Hilbert space of extensive quantities as Further, it is immediate that the Lieb-Robinson bound also holds forτ t = e if t τ t . Thus, by the Lieb-Robinson bound, in this caseτ t is also a one-parameter unitary group onH 0 . We then can define the space ofτ t -invariantsQ, and the orthogonal projectionP Finally, it is easily seen (by p-clustering and the Lieb-Robison bound) that if the dynamical system (U, ι, τ ) and state ω is space-like p c clustering according to Definition 3.3, then so is (U,ι,τ ), ω. From there on, the results we proved above also hold for these new quantities:

Linearised oscillatory Euler equation in one dimension
In [29], it is shown, in one-dimensional quantum lattices (quantum chains), that the theorem of hydrodynamic projection can be used to rigorously obtain linearised Euler equations. These are equations for two-point functions of local conserved densities, in the limit of large times and large wavelengths. In the previous sections, we have seen how both the concepts of ergodiciy and that of hydrodynamic projections hold as well with respect to any wavenumber k and frequency f . In [29], it was commented that indeed also the linearised Euler equation will hold at arbitrary wavenumber and frequency. Thus, this allows us to extend the principles of hydrodynamic correlations to predict oscillatory behaviours. This is perhaps the most important consequence of the observation that large-scale concepts in fact hold with oscillatory factors. In essence, the hydrodynamic expansion becomes an expansion near a nonzero point in the (f, k)-plane.
In this section, we illustrate this principle in a simple quantum chain, in a non-rigorous fashion. Recall Q f,k , Eq. (55): this is the f -oscillatory closed subspace of the Hilbert space H k generated by the k-extensive observables. One may consider the corresponding (f, k)conserved densities and currents A, B ∈ U loc : the density gives rise to the k-extensive charge Such densities and currents should satisfy the oscillatory continuity equation In terms of the infinitesimal generator i[H, ·] for τ t , the k-extensive charge satisfies [H, Σ k A] + f Σ k A = 0. If such a density and current exist, then two-point correlation functions of observables that have nonzero overlap with this density (or more precisely, with Σ k A), should present oscillatory behaviours at large scales of space and time. Indeed, under the (f, k)-fluid-cell mean, defined, as on the r.h.s. of (57), with a factor e ik·x−if t , correlation functions at large scales project onto Q f,k , as per the projection Theorem 3.1. As this holds for the oscillatory fluid-cell mean, that is under averaging with oscillatory factors, this extracts the oscillatory behaviour of the correlation function. For any given model, one expects only certain values of (f, k) that would give nontrivial Q f,k = {0}. That is, only certain type of oscillatory behaviour can be observed (if any). A simple example is the free fermionic lattice, with some dispersion relation E(k). In this system, creation and annihilation operators at momentum k are (E(k), k)-extensive conserved quantities, and oscillatory hydrodynamic projection corresponds to the oscillatory algebraic decay of fermion two-point functions conventionally obtained by a saddle-point analysis. In fact, it is convenient to go beyond just the hydrodynamic projection principle, and get the oscillatory linearised Euler hydrodynamic equations, as in [29]. This gives the explicit oscillatory behaviours of correlation functions of conserved densities. Below we discuss this, and illustrate the phenomenon in a free fermion model. We concentrate on the one-dimensional case D = 1.
(97) (f, k)-hydrodynamic projection is a hydrodynamic expansion, of linear-response type, in the neighbourhood of arbitrary (f, k), which extends the paradigm of hydrodynamics to oscillatory behaviours. It generalises ideas in recent works on time crystals, where dynamical symmetries are used to explain persistent oscillations in Drude weights [24,26].

Oscillatory linearised Euler hydrodynamics
As mentioned, in one-dimensional models, the results of [29] are stronger than those reported in Section 3, as the linearised Euler equation is also proven. This makes the asymptotic of correlation functions more explicit. At a non-rigorous level, the formulae obtained are reviewed in [17] in the usual (non-oscillatory) case. For our purpose, in the non-oscillatory case, it is sufficient to recall that if A i , B i are pairs of conserved densities and currents, such that the extensive charges a i = ΣA i form a basis for the space of conserved quantities Q, then the connected correlation function behaves as where the flux Jacobian is the statistic susceptibility matrix is and it is expected that, for the results to hold, the fluid-cell mean can be taken in the more intuitive form of an average over a cell in space-time: Here the mesoscopic length can be taken as L = L( ) with L → ∞ fast enough as → ∞, and L/ → 0 (or as L = , then taking → 0 + on the asymptotic large-result). Three remarks are in order, see the explanations in [17]: (i) The right-hand side of eq. (98) is obtained from the right-hand side of eq. (45) by (1) expressing the continuity equation relating (A i (x, t), A j (0, 0)) and (B i (x, t), A j (0, 0)); (2) using the projection formula, eq. (45), expressing S ΣB i ,ΣA j (κ) in terms of S ΣA l ,ΣA j (κ), Fourier transforming back to real space to get (B i (x, t), A j (0, 0)) in terms of (A l (x, t)A j (0, 0)); and (3) solving the resulting continuity equation for the matrix of correlators (A i (x, t)A j (0, 0)), with appropriate initial condition. The solution is (98).
(ii) The fluid-cell mean (101) is different from that in (44). The latter is used to express the rigorous projection result, but the result is expected to hold for a variety of possible definitions of fluid-cell means. The fluid-cell mean chosen above is more physically transparent and more convenient for our purposes.
(iii) If there is a finite set of basis charges (a finite set of indices i), the right-hand side of (98) is a "generalised function". Its meaning is that forx/t equal to an eigenvalue of A, the largeasypmtotic decays more slowly than 1/ , while for other velocities, it decays more rapidly. The delta-function comes about from the normalisation condition: where we use 1 = dx and the fact that x∈Z A i (x, t) is independent of time. If there is a continuous set of basis charges, as is typical in integrable models, then the large-asymptotics is exactly 1/ , and the right-hand side of (98) is an ordinary function. See [16,17]. The observation made in [29] is that all formulae apply equally well in the oscillatory case: with the space-time translation group defined as This is because η ω,k x,t still is a unitary operator on the GNS space. Thus we have a space of k-extensive observables, the equivalence classes Σ k A built using the positive semidefinite sesquilinear form and from it we may construct, formally, the (f, k)-conserved charges, Q f,k . Thus, for frequency-wavenumber f, k, if A i , B i are pairs of oscillatory conserved densities and currents, eq. (96), such that a i = Σ k A i form a basis for Q f,k , then where the oscillatory flux Jacobian is the statistic susceptibility matrix is (note that these depend on f because the basis A i 's does), and the fluid-cell mean is

Oscillatory behaviour from a saddle-point analysis
For simplicity, and in order to verify the ideas, we consider a one-dimensional quadratic model, where asymptotics of correlation functions can be obtained by an elementary saddle-point analysis. Such models are integrable, and the usual hydrodynamic projection principle has been studied widely in the context of generalised hydrodynamics [15,16,22,17], and verified to agree with a saddle-point analysis [52].
The mathematical results reported in the main text apply to quantum lattice models with finite local spaces, and are based on the bosonic version of the C * algebra formulation of quantum statistical mechanics, where local operators commute with each other. Hence, they do not cover free bosonic chains (infinite-dimensional local space), nor free fermionic chains (fermionic formulation). Nevertheless, all results are expected to hold in both cases, as they only rely on general properties of correlation functions; in particular, for free fermionic chains, all such properties are well established [35]. In this section, we give the example of the free fermionic chain, with Hamiltonian We are looking to evaluate the asymptotic form of the thermal correlation function as with energy spectrum E(k) = −J cos k and a(k, t) = e iHt a(k)e −iHt = e −iE(k)t a(k).
Therefore we obtain the usual expression for the correlator A similar analysis as that below can be performed in a generalised Gibbs ensemble (GGE), where n(k) takes an arbitrary form that characterises the GGE. The asymptotic behaviour is easily obtained by a saddle point analysis. With k ± solving v(k ± ) = x/t where v(k) = E (k) = J sin k and k + ∈ [−π/2, π/2], k − = sgn(k + )(π − k + ), we obtain as t → ∞ with x/t = ξ fixed. The decay in 1/ √ t, which is slower than 1/t, indicates that the hydrodynamic projection formula should give a generalised function. As all factors explicitly written as functions of k ± are slowly varying, they may be assumed to be constant within the fluid cell. Then we see that the support of the (f, k)-fluid-cell mean (108) is on k = k ± and f = E(k ± ), which is In fact, this saddle point analysis does not provide the full information that we need about the shape of the correlation function around the velocity ξ = v(k). However, because we know the support, we may immediately write, as a consequence of this saddle-point result, The normalisation R(k) is obtained by evaluating where again we use 1 = dx. The left-hand side is found to be e if t n(k) from (114), while the right-hand side is found to be e if t R(k) from (117), giving

Oscillatory behaviour from oscillatory linearised Euler hydrodynamics
We note that e iE(k)t a(k, t) is independent of time. Further, It is a simple matter to verify that the associated current has the form in that A(x, t) and B(x, t) satisfy eq. (96). Then we evaluate and thus Hence, the prediction from hydrodynamic projections is in agreement with (117) and (119).

Conclusion
We have analysed the large-time, long-wavelength behaviours of short-range quantum spin models on Z D .
We have first considered many-body ergodicity, which is ergodicity as viewed from the viewpoint of local physics in the thermodynamic limit. We contrasted it with the standard notion of ergodicity in systems with finite degrees of freedom, and von Neumann's quantum ergodic theorem. In particular, we have obtained strong results in KMS states. We have shown that, under a condition of non-localisation, time averages of local observables become classical, taking the values of their ensemble averages -this is a natural notion of ergodicity. We have also shown, under the same condition, that a time-evolved state within a "microcanonical shell" of the KMS state -defined as perturbations by elements of the C * algebra -is arbitrarily closed to the KMS state 100% of the times over long periods, under a natural metric based on local observables. These results show that the long-time evolution locally looses a large amount of information, in agreement with expected ergodicity.
The non-localisation condition is that the kernel of the time-evolution operator on the GNS space should be "trivial" (one-dimensional): there should not be non-trivial observables (elements of the C * algebra different from 1) that come back to themselves after some time. This is in general hard to check, and it would be interesting to study this condition is explicit models.
Only the new aspects of the proofs of the above many-body ergodicity statements were presented, as the main parts follow directly from the techniques presented in [18].
We have then reviewed the notion of almost-everywhere ergodicity, introduced in [29,18]. This is ergodicity, as above, for almost every ray in space-time. It is shown in [18] to hold for all short-range quantum lattices on Z D . That is, when considering displacements not just in time, but in space-time, there is no need for non-localisation conditions in order for ergodicity to hold almost everywhere. This includes the parallel result on the distance of states in the microcanonical shell, under space-time displacements. This strongly constrains the structure of time-evolved local observables at large times: essentially, local observables become "thin" over time. Intuitively, this is because the finite-dimensional local space of quantum spin lattices does not allow to have operators nontrivially supported over growing regions -only a countable set of space-time rays may remain where the observable is nontrivial.
We have also observed that all these ergodicity results in fact hold with arbitrary oscillatory factors, going beyond conventional ergodicity to strong constraints on how observables may oscillate over long times.
Almost everywhere ergodicity was then used to obtain new rigorous results on hydrodynamic projections. We have shown that the principle of hydrodynamic projection -by which correlation functions, at large wavelengths and long times, project onto the extensive conserved quantities admitted by the model -hold in every short-range quantum spin model on Z D , and at every frequency f and wavenumber k. Projection occurs on f, k-dependent spaces, which are kernels of the time evolution operator on f, k-dependent Hilbert spaces, different from, but constructed in a similar way to, the GNS space.
The hydrodynamic projection theorem naturally generalises that shown in [29] to arbitrary dimension D > 1. However, the proof is a nontrivial extension of that presented there for D = 1, requiring the construction of new Hilbert spaces and the analysis of their properties. The hydrodynamic projection theorem was presented in its most general form, involving arbitrary frequency f and wavenumber k.
The almost-everywhere ergodicity and hydrodynamic projection results underline the large universality of general ergodic and hydrodynamic principles: they emerge solely from the sepa-ration of scales, the large gap between extensivity and locality. Further, the mathematical construction based on Hilbert spaces built from local observables and their correlations shows the inherent flexibility in defining notions of extensivity (or homogeneity) and stationarity, allowing for oscillations.
Finally, using this flexibility, we have considered the oscillatory linearised Euler equation: the linearised Euler equation, but obtained from oscillatory hydrodynamic projections instead of the conventional ones. This was proposed in [29], where the linearised Euler equation was shown for finite-range quantum spin chains. The equation describes the oscillatory behaviours of correlation functions at large space-time separations. We have illustrated it, in a non-rigorous fashion, on a free-fermion chain, showing that the oscillatory hydrodynamic principle correctly reproduces the oscillatory behaviour that can be inferred from a simple saddle-point evaluation of two-point fermion correlation functions at large space-time separations. This gives a proof-ofprinciple for the idea that conventional hydrodynamics can be generalised to describe not only averaged, smoothed-out behaviours at large space-time, but also oscillatory behaviours, with frequencies that are far from 0 and wavelengths that are microscopic.
The results presented open many doors for further studies. Applying the results to more examples, including specific structures inside the LR light-cone, the kernel of the GNS evolution operator and the related many-body ergodicity, and oscillatory hydrodynamic projections in interacting systems generalising the free-fermion example, would be very interesting. It is likely that, in integrable systems, oscillatory hydrodynamic projection would connect with the finitedensity form factor expansions of quantum models; this may also provide a way of defining and studying such form factor expansions in classical models.
The results illustrate how the choice of Hilbert space relates to the large-scale physics of interest; the important elements always lie within the kernel of the evolution operator, but the latter can be taken as acting on local observables, or extensive observables. The idea that different Hilbert spaces relate to different scales was introduced in [53], and it would be interesting to extend the rigorous analysis presented here to the "diffusive" Hilbert space constructed there.
The (f, k)-hydrodynamic projection result paves the way for a full (f, k)-hydrodynamics, a subject which should help uncover new universal dynamics and which we hope to investigate in the near future. In fact, going beyond simple oscillatory phases, it is also possible to re-define space-and time-translation operators by adjoining internal motions of local observables (such as spin rotations). It is clear that all proofs provided immediately generalise to this case whenever the internal motion is an extensive internal symmetry transformation (it commutes with time evolution and acts tensorially on individual spatial cells). It would be extremely interesting to study the related hydrodynamics.
Finally, it would be interesting to go beyond clean spin lattices, towards quantum and classical gases of particles, and also systems with disorder. We believe that the general methods of operator algebras, and in particular the principles presented here, give an alternative set of tools to the traditional ones (such as kinetic or Boltzmann equations), which have the potential to provide new fruitful ways of understanding the emergence of large-scale behaviours from the microscopic physics. current collaboration on related aspects. We also thank Abhishek Dhar for sharing a draft of [23] and for related discussions on ergodicity. BD was supported by EPSRC under the grant "Emergence of hydrodynamics in many-body systems: new rigorous avenues from functional analysis", ref. EP/W000458/1. DA is supported by a studentship from the Engineering and Physical Sciences Research Council.
Definition A.2 (States, Invariance). A state of a dynamical system (U, ι, τ ) is a positive linear functional ω : U → C such that ω = 1. The set of states is denoted by E U . A state ω ∈ E U is called space invariant if ω(ι n (A)) = ω(A), ∀A ∈ U, n ∈ Z D and time invariant if ω(τ t (A)) = ω(A), ∀A ∈ U, t ∈ R. We will refer to a space and time invariant state simply as invariant.
Of great importance in all our results is the Lieb-Robinson bound, which we state as in [54,Corollary 4.3

Appendix B Ergodicity in the various completions of the space of local observables
Consider any positive-semidefinite sesquiliniear form ·, · defined on elements of U loc . For example, this can be the connected correlator (A, B) = ω(A † B) − ω(A † )ω(B), or other forms defined in the main text, such as (A, B) k,⊥ . We define, by the usual construction, a Hilbert space H from ·, · . On it, we assume that we have space-time translation group actions, ι x (for x ∈ Z d ) and τ t . For the case of (A, B) k,⊥ , as summation over the perpendicular plane is done, there remain only a one-dimensional translation group, d = 1. We say that we have the system (H, ι, τ ). We can generalise the definitions of clustering of Section 3.1 for any such forms, so that we define: a subset C ⊂ U loc to be called uniformly p-clustering, with respect to ·, · , for p > D if there exist a C > 0 s.t. for all A, B ∈ C: Likewise, we can adapt Definition 3.3 of space-like p c clustering to ·, · , and we define, as per [18, Def 3.1, Theorem 3.2], ·, · to be space-like ergodic when there exists a υ c such that for every A, B ∈ U loc and υ > υ c : The proof of the almost everywhere ergodicity theorem [18,Theorem 3.2] can now be adapted to the system (H, ι, τ ), as long as space-like ergodicity holds.
Further, with space-like p c -clustering for p c > d, we can define an equivalence relation and similarly to the construction of the spaces of extensive quantities, we can construct a Hilbert space by the completion of the space of equivalence classes. We can also extend the action of ι, τ on the Hilbert space and it will be unitary. This is a matter of replicating the methods in [29,Section 5].

Appendix C Space-like clustering proof
Throughout this appendix we use the notation ι x τ t A for space-time translations of observables A ∈ U, in order to maintain visual clarity of expressions. We prove the following proposition: Theorem C.1. Consider a dynamical system (U, ι, τ ) with interaction satisfying Equation (127), in an exponentially clustering state (Definition 3.2). It follows that the system is space-like exponentially clustering with respect to ω, Definition 3.3, for the Lieb-Robison velocity υ LR .
First, we show the first condition of Definition 3.3. The Lieb-Robinson bound allows us to approximate the time evolved observables by local ones, by projecting the time evoluted τ t (A), A ∈ U loc onto local ones σ Λ (τ t (A)) supported on finite Λ ⊂ Z D . This is done by using the result [55,Corollary 4.4] and satisfies the first condition of Definition 3.3: Lemma C.2. Let A ∈ U and consider a finite Λ ⊂ Z D . If there is an ε > 0 such that then we can approximate A by a strictly local σ Λ (A) ∈ U Λ : This Lemma states that if an observable A ∈ U almost commutes with every B ∈ U Z D \Λ supported outside a finite Λ, then it can be well approximated by a local σ Λ (A) ∈ U Λ . Combined with the Lieb-Robinson bound Equation (129) we can prove (this is also described in [54,Chapter 4.3]): Theorem C.3. Consider a dynamical system with exponentially decaying interactions, the time evolution τ t A of a local A ∈ U Λ and the finite sets Λ r = ∪ x∈Λ B x (r), r = 1, 2, 3, . . . , where B x (r) is the ball of radius r > 0 around x, i.e. Λ r is Λ extended by a distance r around all of its points. Then, we can approximate τ t A by the local σ Λr (τ t A) ∈ U Λr : with ε r = C|Λ|N 2|Λ| exp −λ(r − υ LR |t|), r = 1, 2, 3 . . . .
Let r ∈ N and by Theorem C.3 consider the local approximation of τ t B by σ Λr (τ t B), supported in Λ r = ∪ x∈Λ B x (r): Using this we can approximate the time evolution of observables by a local sequence σ r (τ t A) := σ Λr (τ t A) so that lim r σ r (τ t A) = τ t A.
We show the uniformity condition 2 of the elements σ r (τ t A) in the end of the appendix, we first proceed to show the third condition of Definition 3.3.
It remains to prove the uniformity condition, i.e. that the set {(σ n τ t a, σ m τ s b) : n, m ∈ N} is uniformly clustering, Equation (49) in Definition 3.1.
Proof of 3.3.3. We start with I := |ω(ι n (A)σ Λr (τ t B)) − ω(A)ω(σ Λr (τ t B))| and choose l = dist(Λ A + n, Λ). If r < l then we bound I by the space clustering property, as in Equation (139): and e −λ dist(Λ A +n,Λr) ≤ e −λ|n| e λ diam(Λ A ∪Λ)−λl = e −2λ|n| e 2λ diam(Λ A ∪Λ) . Thus, we get Now, if r > l we proceed as in Equation (137) to get Hence, we can overall obtain a uniform in r exponential clustering. Similarly we can repeat the process for the replacement ι n A → σ r τ s ι n A, which will finally yield uniform exponential clustering for the set {(σ r τ s a, σ r τ t b) : r, r ∈ N} for any t, s ∈ R. This obviously implies uniform p-clustering for all p.

Appendix D Sesquilinear forms properties
Here we prove some basic properties for the sesquilinear forms used throughout the main text.
Lemma D.1. For any k ∈ R D and A ∈ U loc in a p-clustering with p > D dynamical system, it follows A, A k ≥ 0.
We divide both sides of the inequality by (L + 1) D and the Lemma follows by taking the limit L → ∞. That is because the first term becomes A, A k , while the term o((L + 1) D−1 ) tends to 0 for p > D, after diving by (L + 1) D and taking L → ∞. This can be seen when writing o((L + 1) D ) as a sum over all permutations σ of z 1 , z 2 , ..., z D and all k = 1, 2, ..., D, of terms of the form: where the product is taken over a choice of k elements of z 1 , z 2 , ..., z D . We can bound the absolute value of any such term as follows, where we denote S D (d) = {x ∈ Z D : |x| = d} the number of lattice points on a circle of radius d: