Contributor info: Prof. Tomaz Prosen

Alexios Christopoulos, Andrea De Luca, Dmitry L. Kovrizhin, Tomaž Prosen

SciPost Phys. 16, 049 (2024) · published 20 February 2024 | Toggle abstract · pdf

We propose a general exact method of calculating dynamical correlation functions in dual symplectic brick-wall circuits in one dimension. These are deterministic classical many-body dynamical systems which can be interpreted in terms of symplectic dynamics in two orthogonal (time and space) directions. In close analogy with quantum dual-unitary circuits, we prove that two-point dynamical correlation functions are non-vanishing only along the edges of the light cones. The dynamical correlations are exactly computable in terms of a one-site Markov transfer operator, which is generally of infinite dimensionality. We test our theory in a specific family of dual-symplectic circuits, describing the dynamics of a classical Floquet spin chain. Remarkably, expressing these models in the form of a composition of rotations leads to a transfer operator with a block diagonal form in the basis of spherical harmonics. This allows us to obtain analytical predictions for simple local observables. We demonstrate the validity of our theory by comparison with Monte Carlo simulations, displaying excellent agreement with the latter for a choice of observables.

João Costa, Pedro Ribeiro, Andrea De Luca, Tomaž Prosen, Lucas Sá

SciPost Phys. 15, 145 (2023) · published 9 October 2023 | Toggle abstract · pdf

We study spectral and steady-state properties of generic Markovian dissipative systems described by quadratic fermionic Liouvillian operators of the Lindblad form. The Hamiltonian dynamics is modeled by a generic random quadratic operator, i.e., as a featureless superconductor of class D, whereas the Markovian dissipation is described by $M$ random linear jump operators. By varying the dissipation strength and the ratio of dissipative channels per fermion, $m=M/(2N_F)$, we find two distinct phases where the support of the single-particle spectrum has one or two connected components. In the strongly dissipative regime, this transition occurs for $m=1/2$ and is concomitant with a qualitative change in both the steady-state and the spectral gap that rules the large-time dynamics. Above this threshold, the spectral gap and the steady-state purity qualitatively agree with the fully generic (i.e., non-quadratic) case studied recently. Below $m=1/2$, the spectral gap closes in the thermodynamic limit and the steady-state decouples into an ergodic and a nonergodic sector yielding a non-monotonic steady-state purity as a function of the dissipation strength. Our results show that some of the universal features previously observed for fully random Liouvillians are generic for a sufficiently large number of jump operators. On the other hand, if the number of dissipation channels is decreased the system can exhibit nonergodic features, rendering it possible to suppress dissipation in protected subspaces even in the presence of strong system-environment coupling.

Felix Fritzsch, Roopayan Ghosh, Tomaž Prosen

SciPost Phys. 15, 092 (2023) · published 13 September 2023 | Toggle abstract · pdf

We compute the dynamics of entanglement in the minimal setup producing ergodic and mixing quantum many-body dynamics, which we previously dubbed boundary chaos. This consists of a free, non-interacting brickwork quantum circuit, in which chaos and ergodicity is induced by an impurity interaction, i.e., an entangling two-qudit gate, placed at the system's boundary. We compute both the conventional bipartite entanglement entropy with respect to a connected subsystem including the impurity interaction for initial product states as well as the so-called operator entanglement entropy of initial local operators. Thereby we provide exact results in a particular scaling limit of both time and system size going to infinity for either very small or very large subsystems. We show that different classes of impurity interactions lead to very distinct entanglement dynamics. For impurity gates preserving a local product state forming the bulk of the initial state, entanglement entropies of states show persistent spikes with period set by the system size and suppressed entanglement in between, contrary to the expected linear growth in ergodic systems. We observe similar dynamics of operator entanglement for generic impurities. In contrast, for T-dual impurities, which remain unitary under partial transposition, we find entanglement entropies of both states and operators to grow linearly in time with the maximum possible speed allowed by the geometry of the system. The intensive nature of interactions in all cases causes entanglement to grow on extensive time scales proportional to system size.

Žiga Krajnik, Enej Ilievski, Tomaž Prosen, Vincent Pasquier

SciPost Phys. 11, 051 (2021) · published 9 September 2021 | Toggle abstract · pdf

We construct an integrable lattice model of classical interacting spins in discrete space-time, representing a discrete-time analogue of the lattice Landau-Lifshitz ferromagnet with uniaxial anisotropy. As an application we use this explicit discrete symplectic integration scheme to compute the spin Drude weight and diffusion constant as functions of anisotropy and chemical potential. We demonstrate qualitatively different behavior in the easy-axis and the easy-plane regimes in the non-magnetized sector. Upon approaching the isotropic point we also find an algebraic divergence of the diffusion constant, signaling a crossover to spin superdiffusion.

Žiga Krajnik, Enej Ilievski, Tomaž Prosen

SciPost Phys. 9, 038 (2020) · published 14 September 2020 | Toggle abstract · pdf

We introduce a class of integrable dynamical systems of interacting classical matrix-valued fields propagating on a discrete space-time lattice, realized as many-body circuits built from elementary symplectic two-body maps. The models provide an efficient integrable Trotterization of non-relativistic $\sigma$-models with complex Grassmannian manifolds as target spaces, including, as special cases, the higher-rank analogues of the Landau-Lifshitz field theory on complex projective spaces. As an application, we study transport of Noether charges in canonical local equilibrium states. We find a clear signature of superdiffusive behavior in the Kardar-Parisi-Zhang universality class, irrespectively of the chosen underlying global unitary symmetry group and the quotient structure of the compact phase space, providing a strong indication of superuniversal physics.

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.

Katja Klobas, Tomaž Prosen

SciPost Phys. Core 2, 010 (2020) · published 26 June 2020 | Toggle abstract · pdf

In this paper we study the space evolution in the Rule 54 reversible cellular automaton, which is a paradigmatic example of a deterministic interacting lattice gas. We show that the spatial translation of time configurations of the automaton is given in terms of local deterministic maps with the support that is small but bigger than that of the time evolution. The model is thus an example of space-time dual reversible cellular automaton, i.e. its dual is also (in general different) reversible cellular automaton. We provide two equivalent interpretations of the result; the first one relies on the dynamics of quasi-particles and follows from an exhaustive check of all the relevant time configurations, while the second one relies on purely algebraic considerations based on the circuit representation of the dynamics. Additionally, we use the properties of the local space evolution maps to provide an alternative derivation of the matrix product representation of multi-time correlation functions of local observables positioned at the same spatial coordinate.

Bruno Bertini, Pavel Kos, Tomaz Prosen

SciPost Phys. 8, 068 (2020) · published 28 April 2020 | Toggle abstract · pdf

We provide exact results for the dynamics of local-operator entanglement in quantum circuits with two-dimensional wires featuring ultralocal solitons, i.e. single-site operators which, up to a phase, are simply shifted by the time evolution. We classify all circuits allowing for ultralocal solitons and show that only dual-unitary circuits can feature moving ultralocal solitons. Then, we rigorously prove that if a circuit has an ultralocal soliton moving to the left (right), the entanglement of local operators initially supported on even (odd) sites saturates to a constant value and its dynamics can be computed exactly. Importantly, this does not bound the growth of complexity in chiral circuits, where solitons move only in one direction, say to the left. Indeed, in this case we observe numerically that operators on the odd sublattice have unbounded entanglement. Finally, we present a closed-form expression for the local-operator entanglement entropies in circuits with ultralocal solitons moving in both directions. Our results hold irrespectively of integrability.

Bruno Bertini, Pavel Kos, Tomaz Prosen

SciPost Phys. 8, 067 (2020) · published 28 April 2020 | Toggle abstract · pdf

The entanglement in operator space is a well established measure for the complexity of the quantum many-body dynamics. In particular, that of local operators has recently been proposed as dynamical chaos indicator, i.e. as a quantity able to discriminate between quantum systems with integrable and chaotic dynamics. For chaotic systems the local-operator entanglement is expected to grow linearly in time, while it is expected to grow at most logarithmically in the integrable case. Here we study local-operator entanglement in dual-unitary quantum circuits, a class of "statistically solvable" quantum circuits that we recently introduced. We identify a class of "completely chaotic" dual-unitary circuits where the local-operator entanglement grows linearly and we provide a conjecture for its asymptotic behaviour which is in excellent agreement with the numerical results. Interestingly, our conjecture also predicts a "phase transition" in the slope of the local-operator entanglement when varying the parameters of the circuits.

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.

Marko Ljubotina, Spyros Sotiriadis, Tomaž Prosen

SciPost Phys. 6, 004 (2019) · published 10 January 2019 | Toggle abstract · pdf

We study quantum transport after an inhomogeneous quantum quench in a free fermion lattice system in the presence of a localised defect. Using a new rigorous analytical approach for the calculation of large time and distance asymptotics of physical observables, we derive the exact profiles of particle density and current. Our analysis shows that the predictions of a semiclassical approach that has been extensively applied in similar problems match exactly with the correct asymptotics, except for possible finite distance corrections close to the defect. We generalise our formulas to an arbitrary non-interacting particle-conserving defect, expressing them in terms of its scattering properties.