Real-time broadening of bath-induced density profiles from closed-system correlation functions.

The Lindblad master equation is one of the main approaches to open quantum systems. While it has been widely applied in the context of condensed matter systems to study properties of steady states in the limit of long times, the actual route to such steady states has attracted less attention yet. Here, we investigate the nonequilibrium dynamics of spin chains with a local coupling to a single Lindblad bath and analyze the transport properties of the induced magnetization. Combining typicality and equilibration arguments with stochastic unraveling, we unveil for the case of weak driving that the dynamics in the open system can be constructed on the basis of correlation functions in the closed system, which establishes a connection between the Lindblad approach and linear response theory at finite times. In this way, we provide a particular example where closed and open approaches to quantum transport agree strictly. We demonstrate this fact numerically for the spin-1/2 XXZ chain at the isotropic point and in the easy-axis regime, where superdiffusive and diffusive scaling is observed, respectively.

Eigenstate Thermalization Hypothesis and Its Deviations from Random-Matrix Theory beyond the Thermalization Time.

The eigenstate thermalization hypothesis explains the emergence of the thermodynamic equilibrium in isolated quantum many-body systems by assuming a particular structure of the observable's matrix elements in the energy eigenbasis. Schematically, it postulates that off-diagonal matrix elements are random numbers and the observables can be described by random matrix theory (RMT). To what extent a RMT description applies, more precisely at which energy scale matrix elements of physical operators become truly uncorrelated, is, however, not fully understood. We study this issue by introducing a novel numerical approach to probe correlations between matrix elements for Hilbert-space dimensions beyond those accessible by exact diagonalization. Our analysis is based on the evaluation of higher moments of operator submatrices, defined within energy windows of varying width. Considering nonintegrable quantum spin chains, we observe that matrix elements remain correlated even for narrow energy windows corresponding to timescales of the order of thermalization time of the respective observables. We also demonstrate that such residual correlations between matrix elements are reflected in the dynamics of out-of-time-ordered correlation functions.

Nontrivial damping of quantum many-body dynamics.

Understanding how the dynamics of a given quantum system with many degrees of freedom is altered by the presence of a generic perturbation is a notoriously difficult question. Recent works predict that, in the overwhelming majority of cases, the unperturbed dynamics is just damped by a simple function, e.g., exponentially as expected from Fermi's golden rule. While these predictions rely on random-matrix arguments and typicality, they can only be verified for a specific physical situation by comparing to the actual solution or measurement. Crucially, it also remains unclear how frequent and under which conditions counterexamples to the typical behavior occur. In this work, we discuss this question from the perspective of projection-operator techniques, where exponential damping of a density matrix occurs in the interaction picture but not necessarily in the Schrödinger picture. We show that a nontrivial damping in the Schrödinger picture can emerge if the dynamics in the unperturbed system possesses rich features, for instance due to the presence of strong interactions. This suggestion has consequences for the time dependence of correlation functions. We substantiate our theoretical arguments by large-scale numerical simulations of charge transport in the extended Fermi-Hubbard chain, where the nearest-neighbor interactions are treated as a perturbation to the integrable reference system.

Simulating Hydrodynamics on Noisy Intermediate-Scale Quantum Devices with Random Circuits.

In a recent milestone experiment, Google's processor Sycamore heralded the era of "quantum supremacy" by sampling from the output of (pseudo-)random circuits. We show that such random circuits provide tailor-made building blocks for simulating quantum many-body systems on noisy intermediate-scale quantum (NISQ) devices. Specifically, we propose an algorithm consisting of a random circuit followed by a trotterized Hamiltonian time evolution to study hydrodynamics and to extract transport coefficients in the linear response regime. We numerically demonstrate the algorithm by simulating the buildup of spatiotemporal correlation functions in one- and two-dimensional quantum spin systems, where we particularly scrutinize the inevitable impact of errors present in any realistic implementation. Importantly, we find that the hydrodynamic scaling of the correlations is highly robust with respect to the size of the Trotter step, which opens the door to reach nontrivial time scales with a small number of gates. While errors within the random circuit are shown to be irrelevant, we furthermore unveil that meaningful results can be obtained for noisy time evolutions with error rates achievable on near-term hardware. Our work emphasizes the practical relevance of random circuits on NISQ devices beyond the abstract sampling task.

Eigenstate thermalization hypothesis beyond standard indicators: Emergence of random-matrix behavior at small frequencies.

Using numerical exact diagonalization, we study matrix elements of a local spin operator in the eigenbasis of two different nonintegrable quantum spin chains. Our emphasis is on the question to what extent local operators can be represented as random matrices and, in particular, to what extent matrix elements can be considered as uncorrelated. As a main result, we show that the eigenvalue distribution of band submatrices at a fixed energy density is a sensitive probe of the correlations between matrix elements. We find that, on the scales where the matrix elements are in a good agreement with all standard indicators of the eigenstate thermalization hypothesis, the eigenvalue distribution still exhibits clear signatures of the original operator, implying correlations between matrix elements. Moreover, we demonstrate that at much smaller energy scales, the eigenvalue distribution approximately assumes the universal semicircle shape, indicating transition to the random-matrix behavior, and in particular that matrix elements become uncorrelated.

Exponential damping induced by random and realistic perturbations.

Given a quantum many-body system and the expectation-value dynamics of some operator, we study how this reference dynamics is altered due to a perturbation of the system's Hamiltonian. Based on projection operator techniques, we unveil that if the perturbation exhibits a random-matrix structure in the eigenbasis of the unperturbed Hamiltonian, then this perturbation effectively leads to an exponential damping of the original dynamics. Employing a combination of dynamical quantum typicality and numerical linked cluster expansions, we demonstrate that our theoretical findings for random matrices can, in some cases, be relevant for the dynamics of realistic quantum many-body models as well. Specifically, we study the decay of current autocorrelation functions in spin-1/2 ladder systems, where the rungs of the ladder are treated as a perturbation to the otherwise uncoupled legs. We find a convincing agreement between the exact dynamics and the lowest-order prediction over a wide range of interchain couplings.

Relaxation of dynamically prepared out-of-equilibrium initial states within and beyond linear response theory.

We consider a realistic nonequilibrium protocol, where a quantum system in thermal equilibrium is suddenly subjected to an external force. Due to this force, the system is driven out of equilibrium and the expectation values of certain observables acquire a dependence on time. Eventually, upon switching off the external force, the system unitarily evolves under its own Hamiltonian and, as a consequence, the expectation values of observables equilibrate towards specific constant long-time values. Summarizing our main results, we show that, in systems which violate the eigenstate thermalization hypothesis (ETH), this long-time value exhibits an intriguing dependence on the strength of the external force. Specifically, for weak external forces, i.e., within the linear response regime, we show that expectation values thermalize to their original equilibrium values, despite the ETH being violated. In contrast, for stronger perturbations beyond linear response, the quantum system relaxes to some nonthermal value which depends on the previous nonequilibrium protocol. While we present theoretical arguments which underpin these results, we also numerically demonstrate our findings by studying the real-time dynamics of two low-dimensional quantum spin models.

Impact of eigenstate thermalization on the route to equilibrium.

The eigenstate thermalization hypothesis (ETH) and the theory of linear response (LRT) are celebrated cornerstones of our understanding of the physics of many-body quantum systems out of equilibrium. While the ETH provides a generic mechanism of thermalization for states arbitrarily far from equilibrium, LRT extends the successful concepts of statistical mechanics to situations close to equilibrium. In our work, we connect these cornerstones to shed light on the route to equilibrium for a class of properly prepared states. We unveil that, if the off-diagonal part of the ETH applies, then the relaxation process can become independent of whether or not a state is close to equilibrium. Moreover, in this case, the dynamics is generated by a single correlation function, i.e., the relaxation function in the context of LRT. Our analytical arguments are illustrated by numerical results for idealized models of random-matrix type and more realistic models of interacting spins on a lattice. Remarkably, our arguments also apply to integrable quantum systems where the diagonal part of the ETH may break down.

Relation between far-from-equilibrium dynamics and equilibrium correlation functions for binary operators.

Linear response theory (LRT) is one of the main approaches to the dynamics of quantum many-body systems. However, this approach has limitations and requires, e.g., that the initial state is (i) mixed and (ii) close to equilibrium. In this paper, we discuss these limitations and study the nonequilibrium dynamics for a certain class of properly prepared initial states. Specifically, we consider thermal states of the quantum system in the presence of an additional static force which, however, become nonequilibrium states when this static force is eventually removed. While for weak forces the relaxation dynamics is well captured by LRT, much less is known in the case of strong forces, i.e., initial states far away from equilibrium. Summarizing our main results, we unveil that, for high temperatures, the nonequilibrium dynamics of so-called binary operators is always generated by an equilibrium correlation function. In particular, this statement holds true for states in the far-from-equilibrium limit, i.e., outside the linear response regime. In addition, we confirm our analytical results by numerically studying the dynamics of local fermionic occupation numbers and local energy densities in the spin-1/2 Heisenberg chain. Remarkably, these simulations also provide evidence that our results qualitatively apply in a more general setting, e.g., in the anisotropic XXZ model where the local energy is a non-binary operator, as well as for a wider range of temperature. Furthermore, exploiting the concept of quantum typicality, all of our findings are not restricted to mixed states, but are valid for pure initial states as well.

Pull or Push? Octopuses Solve a Puzzle Problem.

Octopuses have large brains and exhibit complex behaviors, but relatively little is known about their cognitive abilities. Here we present data from a five-level learning and problem-solving experiment. Seven octopuses (Octopus vulgaris) were first trained to open an L shaped container to retrieve food (level 0). After learning the initial task all animals followed the same experimental protocol, first they had to retrieve this L shaped container, presented at the same orientation, through a tight fitting hole in a clear Perspex partition (level 1). This required the octopuses to perform both pull and release or push actions. After reaching criterion the animals advanced to the next stage of the test, which would be a different consistent orientation of the object (level 2) at the start of the trial, an opaque barrier (level 3) or a random orientation of the object (level 4). All octopuses were successful in reaching criterion in all levels of the task. At the onset of each new level the performance of the animals dropped, shown as an increase in working times. However, they adapted quickly so that overall working times were not significantly different between levels. Our findings indicate that octopuses show behavioral flexibility by quickly adapting to a change in a task. This can be compared to tests in other species where subjects had to conduct actions comprised of a set of motor actions that cannot be understood by a simple learning rule alone.

Octopus arm movements under constrained conditions: adaptation, modification and plasticity of motor primitives.

The motor control of the eight highly flexible arms of the common octopus (Octopus vulgaris) has been the focus of several recent studies. Our study is the first to manage to introduce a physical constraint to an octopus arm and investigate the adaptability of stereotypical bend propagation in reaching movements and the pseudo-limb articulation during fetching. Subjects (N=6) were placed inside a transparent Perspex box with a hole at the center that allowed the insertion of a single arm. Animals had to reach out through the hole toward a target, to retrieve a food reward and fetch it. All subjects successfully adjusted their movements to the constraint without an adaptation phase. During reaching tasks, the animals showed two movement strategies: stereotypical bend propagation reachings, which were established at the hole of the Perspex box and variant waving-like movements that showed no bend propagations. During fetching movements, no complete pseudo-joint fetching was observed outside the box and subjects pulled their arms through the hole in a pull-in like movement. Our findings show that there is some flexibility in the octopus motor system to adapt to a novel situation. However, at present, it seems that these changes are more an effect of random choices between different alternative motor programs, without showing clear learning effects in the choice between the alternatives. Interestingly, animals were able to adapt the fetching movements to the physical constraint, or as an alternative explanation, they could switch the motor primitive fetching to a different motor primitive 'arm pulling'.