Contributor info: Dr Marko Medenjak

Denis Bernard, Fabian H. L. Essler, Ludwig Hruza, Marko Medenjak

SciPost Phys. 12, 042 (2022) · published 28 January 2022 | Toggle abstract · pdf

We consider the dynamics of fluctuations in the quantum asymmetric simple exclusion process (Q-ASEP) with periodic boundary conditions. The Q-ASEP describes a chain of spinless fermions with random hoppings that are induced by a Markovian environment. We show that fluctuations of the fermionic degrees of freedom obey evolution equations of Lindblad type, and derive the corresponding Lindbladians. We identify the underlying algebraic structure by mapping them to non-Hermitian spin chains and demonstrate that the operator space fragments into exponentially many (in system size) sectors that are invariant under time evolution. At the level of quadratic fluctuations we consider the Lindbladian on the sectors that determine the late time dynamics for the particular case of the quantum symmetric simple exclusion process (Q-SSEP). We show that the corresponding blocks in some cases correspond to known Yang-Baxter integrable models and investigate the level-spacing statistics in others. We carry out a detailed analysis of the steady states and slow modes that govern the late time behaviour and show that the dynamics of fluctuations of observables is described in terms of closed sets of coupled linear differential-difference equations. The behaviour of the solutions to these equations is essentially diffusive but with relevant deviations, that at sufficiently late times and large distances can be described in terms of a continuum scaling limit which we construct. We numerically check the validity of this scaling limit over a significant range of time and space scales. These results are then applied to the study of operator spreading at large scales, focusing on out-of-time ordered correlators and operator entanglement.

Marko Medenjak, Jacopo De Nardis, Takato Yoshimura

SciPost Phys. 9, 075 (2020) · published 19 November 2020 | Toggle abstract · pdf

We introduce non-trivial contributions to diffusion constant in generic many-body systems arising from quadratic fluctuations of ballistically propagating, i.e. convective, modes. Our result is obtained by expanding the current operator in the vicinity of equilibrium states in terms of powers of local and quasi-local conserved quantities. We show that only the second-order terms in this expansion carry a finite contribution to diffusive spreading. Our formalism implies that whenever there are at least two coupled modes with degenerate group velocities, the system behaves super-diffusively, in accordance with the non-linear fluctuating hydrodynamics theory. Finally, we show that our expression saturates the exact diffusion constants in quantum and classical interacting integrable systems, providing a general framework to derive these expressions.

Marko Medenjak, Tomaz Prosen, Lenart Zadnik

SciPost Phys. 9, 003 (2020) · published 7 July 2020 | Toggle abstract · pdf

Dynamical response functions are standard tools for probing local physics near the equilibrium. They provide information about relaxation properties after the equilibrium state is weakly perturbed. In this paper we focus on systems which break the assumption of thermalization by exhibiting persistent temporal oscillations. We provide rigorous bounds on the Fourier components of dynamical response functions in terms of extensive or local dynamical symmetries, i.e. extensive or local operators with periodic time dependence. Additionally, we discuss the effects of spatially inhomogeneous dynamical symmetries. The bounds are explicitly implemented on the example of an interacting Floquet system, specifically in the integrable Trotterization of the Heisenberg XXZ model.

Marko Medenjak, Vladislav Popkov, Tomaž Prosen, Eric Ragoucy, Matthieu Vanicat

SciPost Phys. 6, 074 (2019) · published 24 June 2019 | Toggle abstract · pdf

In this paper we study the statistical properties of a reversible cellular automaton in two out-of-equilibrium settings. In the first part we consider two instances of the initial value problem, corresponding to the inhomogeneous quench and the local quench. Our main result is an exact matrix product expression of the time evolution of the probability distribution, which we use to determine the time evolution of the density profiles analytically. In the second part we study the model on a finite lattice coupled with stochastic boundaries. Once again we derive an exact matrix product expression of the stationary distribution, as well as the particle current and density profiles in the stationary state. The exact expressions reveal the existence of different phases with either ballistic or diffusive transport depending on the boundary parameters.