Abstract We present quantitative predictions for quantum simulator experiments on Ising models from trapped ions to Rydberg chains and show how the thermalization, and thus decoherence times, can be controlled by considering common, independent, and end-cap couplings to the bath. We find (i) independent baths enable more rapid thermalization in comparison to a common one; (ii) the thermalization timescale depends strongly on the position in the Ising phase diagram; (iii) for a common bath larger system sizes show a significant slow down in the thermalization process; and (iv) finite-size scaling indicates a subradiance effect slowing thermalization rates toward the infinite spin chain limit. We find it is necessary to treat the full multi-channel Lindblad master equation rather than the commonly used single-channel local Lindblad approximation to make accurate predictions on a classical computer. This method reduces the number of qubits one can practically classical simulate by at least a factor of 4, in turn showing a quantum advantage for such thermalization problems at a factor of 4 smaller qubit number for open quantum systems as opposed to closed ones. Thus, our results encourage open quantum system exploration in noisy intermediate-scale quantum technologies.
Abstract We present quantitative predictions for quantum simulator experiments on Ising models from trapped ions to Rydberg chains and show how the thermalization, and thus decoherence times, can be controlled by considering common, independent, and end-cap couplings to the bath. We find (i) independent baths enable more rapid thermalization in comparison to a common one; (ii) the thermalization timescale depends strongly on the position in the Ising phase diagram; (iii) for a common bath larger system sizes show a significant slow down in the thermalization process; and (iv) finite-size scaling indicates a subradiance effect slowing thermalization rates toward the infinite spin chain limit. We find it is necessary to treat the full multi-channel Lindblad master equation rather than the commonly used single-channel local Lindblad approximation to make accurate predictions on a classical computer. This method reduces the number of qubits one can practically classical simulate by at least a factor of 4, in turn showing a quantum advantage for such thermalization problems at a factor of 4 smaller qubit number for open quantum systems as opposed to closed ones. Thus, our results encourage open quantum system exploration in noisy intermediate-scale quantum technologies.
We study decoherence induced on a two-level system coupled to a one-dimensional quantum spin chain. We consider the cases where the dynamics of the chain is determined by the Ising, XY, or Heisenberg exchange Hamiltonian. This model of quantum baths can be of fundamental importance for the understanding of decoherence in open quantum systems, since it can be experimentally engineered by using atoms in optical lattices. As an example, here we show how to implement a pure dephasing model for a qubit system coupled to an interacting spin bath. We provide results that go beyond the case of a central spin coupled uniformly to all the spins of the bath, in particular showing what happens when the bath enters different phases, or becomes critical; we also study the dependence of the coherence loss on the number of bath spins to which the system is coupled and we describe a coupling-independent regime in which decoherence exhibits universal features, irrespective of the system-environment coupling strength. Finally, we establish a relation between decoherence and entanglement inside the bath. For the Ising and the XY models we are able to give an exact expression for the decay of coherences, while for the Heisenberg bath we resort to the numerical time-dependent density matrix renormalization group.
We address the out-of-equilibrium dynamics arising from quantum-quench (QQ) protocols (instantaneous changes of the Hamiltonian parameters) in many-body systems within their quantum critical regime and in contact with (homogeneously coupled) thermal baths. We consider two classes of QQ protocols. In one of them, the thermal bath is used to prepare the initial finite-temperature Gibbs state; then, after quenching, the thermal bath is removed and the dynamics of the system is unitary. We also address a more complex QQ protocol where the thermal bath is not removed after quenching, thus the quantum evolution is also driven by the interaction with the bath, which may be described by appropriate master equations for the density matrix of the system, where a further relevant time scale, or inverse decay rate, characterizes the system-bath coupling. Under these QQ protocols, the critical system develops out-of-equilibrium scaling behaviors, which extend those for isolated critical systems, by introducing further scaling variables proportional to the temperature of the thermal bath and the decay rate of the system-bath interactions. These out-of-equilibrium scaling behaviors are checked by analyzing QQ protocols within fermionic Kitaev wires, or equivalently quantum Ising chains, supplemented with a particular modelization of thermal bath that guarantees the asymptotic thermalization within the Lindblad master equation for the dynamics of open systems.
Exactly solvable models provide an opportunity to study different aspects of reduced quantum dynamics in detail. We consider the reduced dynamics of a single spin in finite XX and XY spin 1/2 chains. First we introduce a general expression describing the evolution of the reduced density matrix. This expression proves to be tractable when the combined closed system (i.e. open system plus environment) is integrable. Then we focus on comparing decoherence and thermalization timescales in the XX chain. We find that for a single spin these timescales are comparable, in contrast to what should be expected for a macroscopic body. This indicates that the process of quantum relaxation of a system with few accessible states can not be separated in two distinct stages - decoherence and thermalization. Finally, we turn to finite-size effects in the time evolution of a single spin in the XY chain. We observe three consecutive stages of the evolution: regular evolution, partial revivals, irregular (apparently chaotic) evolution. The duration of the regular stage is proportional to the number of spins in the chain. We observe a "quiet and cold period" in the end of the regular stage, which breaks up abruptly at some threshold time.
Many physical phenomena, including thermalization in open quantum systems and quantum Gibbs sampling, are modeled by Lindbladians approximating a system weakly coupled to a bath. Understanding the convergence speed of these Lindbladians to their steady states is crucial for bounding algorithmic runtimes and thermalization timescales. We study two such families of processes: one characterizing a repeated-interaction Gibbs sampling algorithm, and another modeling open many-body quantum thermalization. We prove that both converge in polynomial time for several non-commuting systems, including high-temperature local lattices, weakly interacting fermions, and 1D spin chains. These results demonstrate that simple dissipative quantum algorithms can prepare complex Gibbs states and that Lindblad dynamics accurately capture thermal relaxation. Our proofs rely on a novel technical result that extrapolates spectral gap lower bounds from quasi-local Lindbladians to the non-local generators governing these dynamics.
The occupation number is a key observable for diagnosing thermalization, as it connects directly to standard statistical laws such as Fermi-Dirac, Bose-Einstein, and Boltzmann distributions. In the context of spin systems, it represents the population of the sublevels of the magnetization in the z direction. We use this quantity to probe the onset of thermalization in an isolated one-dimensional quantum spin-1 Ising model with transverse and longitudinal fields and in its classical counterpart. Thermalization is achieved when the long-time average of the occupation number converges to the microcanonical prediction as the chain length L increases, consistent with the emergence of ergodicity. However, the finite-size scaling analysis in the quantum model is challenged by the exponential growth of the Hilbert space with L. To overcome this limitation, we turn to the classical model, which enables access to much larger system sizes. By tracking the dynamics of individual spins on their three-dimensional Bloch spheres and employing tools from random matrix theory, we establish a quantitative criterion for classical ergodicity in interacting spin systems. We find that deviations from classical ergodicity decay algebraically with system size. This power-law scaling then provides a quantitative bound on the approach to thermal equilibrium in the quantum model.
Abstract Tensor network methods as presented in our open source matrix product states software have opened up the possibility to study many-body quantum physics in one and quasi-one-dimensional systems in an easily accessible package similar to density functional theory codes but for strongly correlated dynamics. Here, we address methods which allow one to capture the full entanglement without truncation of the Hilbert space. Such methods are suitable for validation of and comparisons to tensor network algorithms, but especially useful in the case of new kinds of quantum states with high entanglement violating the truncation in tensor networks. Quantum cellular automata are one example for such a system, characterized by tunable complexity, entanglement, and a large spread over the Hilbert space. Beyond the evolution of pure states as a closed system, we adapt the techniques for open quantum systems simulated via the Lindblad master equation. We present three algorithms for solving closed-system many-body time evolution without truncation of the Hilbert space. Exact diagonalization methods have the advantage that they not only keep the full entanglement but also require no approximations to the propagator. Seeking the limits of a maximal number of qubits on a single core, we use Trotter decompositions or Krylov approximation to the exponential of the Hamiltonian. All three methods are also implemented for open systems represented via the Lindblad master equation built from local channels. We show their convergence parameters and focus on efficient schemes for their implementations including Abelian symmetries, e.g. symmetry used for number conservation in the Bose–Hubbard model or discrete symmetries in the quantum Ising model. We present the thermalization timescale in the long-range quantum Ising model as a key example of how exact diagonalization contributes to novel physics.