An isolated quantum system evolves reversibly, yet a small subsystem can lose nearly all information about its initial state. Conserved quantities may retain some of this information, depending on how strongly they overlap with local observables. Here we consider a subsystem in an arbitrary initial state coupled to a maximally mixed bath. We show that the worst-input mean-square equilibration error, maximised over window origins, is bounded from above and below by the mean square of a complete local autocorrelation, with the stationary part subtracted. The two mean-square quantities are linearly comparable. A correlation record over a finite interval therefore controls the equilibration error in every interval of the same length, and the result extends to nonnegative observation weights with nonnegative Fourier transform. The result allows arbitrary degeneracies of energies and gaps. We apply it to two static Hamiltonian families. For fixed spin subsystems in the one-particle sector of a bounded-range periodic clock, fixed-accuracy thermalisation is guaranteed on windows that grow only slightly faster than linearly with system size, despite exponentially small positive energies. A random commuting model with logarithmic interaction range retains an extensive set of conserved observables. Its expected mean-square error is bounded by a finite-size term and an inverse-time term. The inverse-time dependence is unavoidable because of the initial transient. These examples show how local relaxation can coexist with exact conservation laws and a finely spaced spectrum.
Comments
comments welcome; 13 pages main plus supplementary, 2 figures, v2: substantially rewritten version with new results
An isolated quantum system evolves reversibly, yet a small subsystem can lose nearly all information about its initial state. Conserved quantities may retain some of this information, depending on how strongly they overlap with local observables. Here we consider a subsystem in an arbitrary initial state coupled to a maximally mixed bath. We show that the worst-input mean-square equilibration error, maximised over window origins, is bounded from above and below by the mean square of a complete local autocorrelation, with the stationary part subtracted. The two mean-square quantities are linearly comparable. A correlation record over a finite interval therefore controls the equilibration error in every interval of the same length, and the result extends to nonnegative observation weights with nonnegative Fourier transform. The result allows arbitrary degeneracies of energies and gaps. We apply it to two static Hamiltonian families. For fixed spin subsystems in the one-particle sector of a bounded-range periodic clock, fixed-accuracy thermalisation is guaranteed on windows that grow only slightly faster than linearly with system size, despite exponentially small positive energies. A random commuting model with logarithmic interaction range retains an extensive set of conserved observables. Its expected mean-square error is bounded by a finite-size term and an inverse-time term. The inverse-time dependence is unavoidable because of the initial transient. These examples show how local relaxation can coexist with exact conservation laws and a finely spaced spectrum.
Markovian open-system dynamics have widespread applications throughout quantum information science, including algorithmic state preparation. Their convergence is commonly quantified using the worst case global trace distance between the evolving and stationary states. However, this criterion can be unnecessarily stringent when only physically relevant observables are of interest. Here we introduce and study observable-specific mixing times. We prove that, for quasi-local, rapidly mixing Lindbladians, sums of geometrically local observables equilibrate in a time independent of system size, in contrast to the logarithmic dependence of global state mixing. This separation reduces the runtime of dissipative quantum algorithms, including quantum Gibbs samplers, for estimating quantities such as the Gibbs state energy and local order parameters, yielding an overall scaling that is linear in system size. Complementing this quantum result, we develop a quantum-inspired classical algorithm for estimating the same quantities. Its runtime is likewise linear in system size, but scaling exponentially in O(log(1/ε)D), where D denotes the spatial dimension of the lattice. We further analyse non-interacting Lindbladians over qudits, fermions, and bosons, demonstrating that locality of observables is not always necessary for a qualitatively faster mixing. Small-scale simulations of quantum Gibbs samplers reveal no large hidden constants in our asymptotic analysis and show that the theoretical predictions closely capture the finite-size dynamics.
Quantum systems governed by unitary and reversible microscopic dynamics may nevertheless exhibit equilibration, in the sense that some effective description becomes time-independent. Standard equilibration results usually consider two separate situations: system-environment structures, in which the composite system evolves unitarily while the system of interest equilibrates, and restricted measurements, such as coarse-grained POVMs and observables, in which the measurement statistics equilibrate. Here, we bring these descriptions into a common state-level framework using the concept of generalized subsystems, where the accessible effective state appears as the output of a quantum channel acting on the microscopic state. We derive bounds showing that generalized subsystems equilibrate when their dimension is small compared with the effective dimension of the discarded microscopic information. We further show that this condition is met for typical initial states in large subspaces and that the resulting equilibrium description is largely insensitive to microscopic initial details. The framework recovers the usual equilibration bounds for ordinary subsystems and finite families of POVMs. As an example, we also introduce a finite-resolution energy channel that maps unresolved microscopic energy levels into effective energy levels, thereby making residual effective coherences explicit and showing how spectral multiplicities constrain those coherences while strengthening equilibration. Our results provide a unified state-level formulation of quantum equilibration under general forms of limited accessible information.
Using a Krylov-subspace time evolution algorithm, we simulate the real-time dynamics of translation invariant non-integrable finite spin rings to quite long times with high accuracy. We systematically study the finite-size deviation between the resulting equilibrium state and the thermal state, and we highlight the importance of the energy variance on the deviations. We find that the deviations are well described by the eigenstate thermalization hypothesis, and that the von Neumann entropy correction scaling is the square of the local operator scaling. We reveal also an area law contribution to the relaxed von Neumann entropy, which we connect to the mutual information between the considered subsystem and its immediate environment. We also find that local observables relax towards equilibrium exponentially with a relaxation time scale that grows linearly with system length and is somewhat independent of the local operator, but depends strongly on the energy of the initial state, with the fastest relaxation times found towards one end of the overall energy spectrum. To contrast this behaviour we also study domain wall initial states, which exhibit clear diffusive behaviour, with a Thouless time scaling quadratically with the systems size, leading to a rather precise estimate for the diffusion constant for states in the vicinity of the middle of the energy spectrum.
Temperature is one of the central concepts of thermodynamics, yet its meaning far from equilibrium remains unclear. The problem is especially challenging in isolated quantum many-body systems, whose states evolve unitarily, may remain far from equilibrium, and retain energy coherence, a genuinely quantum feature with no direct classical counterpart. Specifically, energy fluctuations in a nonstationary quantum state have two distinct origins. Part of these fluctuations comes from uncertainty in the energy populations and has the usual thermodynamic meaning. The rest comes from quantum coherence between energy sectors and is responsible for the state's time dependence. We propose that, even away from equilibrium, temperature identifies the state within the family of regular states sharing the same energy-coherence structure. The resulting inverse temperature is obtained from a principle of minimum discrimination information, which selects the least distinguishable state compatible with the prescribed energy and coherence structure. This replaces the usual maximum-entropy principle, which, applied naively, can yield a biased estimator of the inverse temperature. We also lay the groundwork for an extension of the construction to subsystems and show that, although the inverse temperature of a subsystem is not determined by its reduced state alone, under an alignment assumption, its instantaneous rate of change is a local quantity for sufficiently large subsystems.
We consider time evolution in models close to integrable points with hidden symmetries that generate infinitely many local conservation laws that do not commute with one another. The system is expected to (locally) relax to a thermal ensemble if integrability is broken, or to a so-called generalised Gibbs ensemble if unbroken. In some circumstances expectation values exhibit quasi-stationary behaviour long before their typical relaxation time. For integrability-breaking perturbations, these are also called pre-thermalisation plateaux, and emerge e.g. in the strong coupling limit of the Bose–Hubbard model. As a result of the hidden symmetries, quasi-stationarity appears also in integrable models, for example in the Ising limit of the XXZ model. We investigate a weak coupling limit, identify a time window in which the effects of the perturbations become significant and solve the time evolution through a mean-field mapping. As an explicit example we study the XYZ spin- chain with additional perturbations that break integrability. One of the most intriguing results of the analysis is the appearance of persistent oscillatory behaviour. To unravel its origin, we study in detail a toy model: the transverse-field Ising chain with an additional nonlocal interaction proportional to the square of the transverse spin per unit length (2013 Phys. Rev. Lett . 111 197203). Despite being nonlocal, this belongs to a class of models that emerge as intermediate steps of the mean-field mapping and shares many dynamical properties with the weakly interacting models under consideration.
One of the main methods for protecting quantum information against decoherence is to encode information in the ground subspace (or the low energy sector) of a Hamiltonian with a large energy gap which penalizes errors from environment. The protecting Hamiltonian is chosen such that its degenerate ground subspace is an error detecting code for the errors caused by the interaction with the environment. We consider environments with arbitrary number of local sites, e.g. spins, whose interactions among themselves are local and bounded. Then, assuming the system is interacting with a finite number of sites in the environment, we prove that, up to second order with respect to the coupling constant, decoherence and relaxation are suppressed by a factor which grows exponentially fast with the ratio of energy penalty to the norm of local interactions in the environment. The state may, however, still evolve unitarily inside the code subspace due to the Lamb shift effect. In the context of adiabatic quantum computation, this means that the evolution inside the code subspace is effectively governed by a renormalized Hamiltonian. The result is derived from first principles, without use of master equations or their assumptions, and holds even in the infinite temperature limit. We also prove that unbounded or non-local interactions in the environment at sites far from the system do not considerably modify the exponential suppression. Our main technical tool is a new bound on the decay of power spectral density at high frequencies for local observables and many-body Hamiltonians with bounded and local interactions in a neighborhood around the support of the observable.