-
Microcanonical Free Cumulants in lattice systems
Authors:
Felix Fritzsch,
Tomaž Prosen,
Silvia Pappalardi
Abstract:
Recently, the full version of the Eigenstate Thermalization Hypothesis (ETH) has been systematized using Free Probability. In this paper, we present a detailed discussion of the Free Cumulants approach to many-body dynamics within the microcanonical ensemble. Differences between the later and canonical averages are known to manifest in the time-dependent fluctuations of extensive operators. Thus,…
▽ More
Recently, the full version of the Eigenstate Thermalization Hypothesis (ETH) has been systematized using Free Probability. In this paper, we present a detailed discussion of the Free Cumulants approach to many-body dynamics within the microcanonical ensemble. Differences between the later and canonical averages are known to manifest in the time-dependent fluctuations of extensive operators. Thus, the microcanonical ensemble is essential to extend the application of Free Probability to the broad class of extensive observables. We numerically demonstrate the validity of our approach in a non-integrable spin chain Hamiltonian for extensive observables at finite energy density. Our results confirm the full ETH properties, specifically the suppression of crossing contributions and the factorization of non-crossing ones, thus demonstrating that the microcanonical free cumulants encode ETH smooth correlations for both local and extensive observables.
△ Less
Submitted 2 September, 2024;
originally announced September 2024.
-
Bethe-ansatz diagonalization of steady state of boundary driven integrable spin chains
Authors:
Vladislav Popkov,
Xin Zhang,
Carlo Presilla,
Tomaz Prosen
Abstract:
We find that the non-equilibrium steady state (NESS) of integrable spin chains undergoing boundary dissipation, can be described in terms of quasiparticles, with renormalized -- dissipatively dressed -- dispersion relation. The spectrum of the NESS is then fully accounted for by Bethe ansatz equations for a related coherent system, described by a dissipation-projected Hamiltonian of the original s…
▽ More
We find that the non-equilibrium steady state (NESS) of integrable spin chains undergoing boundary dissipation, can be described in terms of quasiparticles, with renormalized -- dissipatively dressed -- dispersion relation. The spectrum of the NESS is then fully accounted for by Bethe ansatz equations for a related coherent system, described by a dissipation-projected Hamiltonian of the original system. We find explicit analytic expressions for the dressed energies of $XXX$ and $XXZ$ models with effective, i.e., induced by the dissipation, diagonal boundary fields, which are U(1) invariant, as well as XXZ and XYZ models with effective non-diagonal boundary fields. In all cases, the dissipative dressing generates an extra singularity in the dispersion relation, which strongly modifies the nonequilibrium steady state spectrum with respect to the spectrum of the corresponding coherent model. This leads, in particular, to a dissipation-assisted entropy reduction, due to the suppression in the NESS spectrum of plain wave-type Bethe states in favor of Bethe states localized at the boundaries.
△ Less
Submitted 17 August, 2024;
originally announced August 2024.
-
Symmetry classes of classical stochastic processes
Authors:
Lucas Sá,
Pedro Ribeiro,
Tomaž Prosen,
Denis Bernard
Abstract:
We perform a systematic symmetry classification of the Markov generators of classical stochastic processes. Our classification scheme is based on the action of involutive symmetry transformations of a real Markov generator, extending the Bernard-LeClair scheme to the arena of classical stochastic processes and leading to a set of up to fifteen allowed symmetry classes. We construct families of sol…
▽ More
We perform a systematic symmetry classification of the Markov generators of classical stochastic processes. Our classification scheme is based on the action of involutive symmetry transformations of a real Markov generator, extending the Bernard-LeClair scheme to the arena of classical stochastic processes and leading to a set of up to fifteen allowed symmetry classes. We construct families of solutions of arbitrary matrix dimensions for five of these classes with a simple physical interpretation of particles hopping on multipartite graphs. In the remaining classes, such a simple construction is prevented by the positivity of entries of the generator particular to classical stochastic processes, which imposes a further requirement beyond the usual symmetry classification constraints. We partially overcome this difficulty by resorting to a stochastic optimization algorithm, finding specific examples of generators of small matrix dimensions in six further classes, leaving the existence of the final four allowed classes an open problem. Our symmetry-based results unveil new possibilities in the dynamics of classical stochastic processes: Kramers degeneracy of eigenvalue pairs, dihedral symmetry of the spectra of Markov generators, and time reversal properties of stochastic trajectories and correlation functions.
△ Less
Submitted 25 June, 2024;
originally announced June 2024.
-
Integrability of open boundary driven quantum circuits
Authors:
Chiara Paletta,
Tomaž Prosen
Abstract:
In this paper, we address the problem of Yang-Baxter integrability of doubled quantum circuit of qubits (spins 1/2) with open boundary conditions where the two circuit replicas are only coupled at the left or right boundary. We investigate the cases where the bulk is given by elementary six vertex unitary gates of either the free fermionic XX type or interacting XXZ type. By using the Sklyanin's c…
▽ More
In this paper, we address the problem of Yang-Baxter integrability of doubled quantum circuit of qubits (spins 1/2) with open boundary conditions where the two circuit replicas are only coupled at the left or right boundary. We investigate the cases where the bulk is given by elementary six vertex unitary gates of either the free fermionic XX type or interacting XXZ type. By using the Sklyanin's construction of reflection algebra, we obtain the most general solutions of the boundary Yang-Baxter equation for such a setup. We use this solution to build, from the transfer matrix formalism, integrable circuits with two step discrete time Floquet (aka brickwork) dynamics. We prove that, only if the bulk is a free-model, the boundary matrices are in general non-factorizable, and for particular choice of free parameters yield non-trivial unitary dynamics with boundary interaction between the two chains. Then, we consider the limit of continuous time evolution and we give the interpretation of a restricted set of the boundary terms in the Lindbladian setting. Specifically, for a particular choice of free parameters, the solutions correspond to an open quantum system dynamics with the source terms representing injecting or removing particles from the boundary of the spin chain.
△ Less
Submitted 17 July, 2024; v1 submitted 18 June, 2024;
originally announced June 2024.
-
Quantum many-body spin ratchets
Authors:
Lenart Zadnik,
Marko Ljubotina,
Žiga Krajnik,
Enej Ilievski,
Tomaž Prosen
Abstract:
Introducing a class of SU(2) invariant quantum unitary circuits generating chiral transport, we examine the role of broken space-reflection and time-reversal symmetries on spin transport properties. Upon adjusting parameters of local unitary gates, the dynamics can be either chaotic or integrable. The latter corresponds to a generalization of the space-time discretized (Trotterized) higher-spin qu…
▽ More
Introducing a class of SU(2) invariant quantum unitary circuits generating chiral transport, we examine the role of broken space-reflection and time-reversal symmetries on spin transport properties. Upon adjusting parameters of local unitary gates, the dynamics can be either chaotic or integrable. The latter corresponds to a generalization of the space-time discretized (Trotterized) higher-spin quantum Heisenberg chain. We demonstrate that breaking of space-reflection symmetry results in a drift in the dynamical spin susceptibility. Remarkably, we find a universal drift velocity given by a simple formula which, at zero average magnetization, depends only on the values of SU(2) Casimir invariants associated with local spins. In the integrable case, the drift velocity formula is confirmed analytically based on the exact solution of thermodynamic Bethe ansatz equations. Finally, by inspecting the large fluctuations of the time-integrated current between two halves of the system in stationary maximum-entropy states, we demonstrate violation of the Gallavotti-Cohen symmetry, implying that such states cannot be regarded as equilibrium ones. We show that the scaled cumulant generating function of the time-integrated current instead obeys a generalized fluctuation relation.
△ Less
Submitted 3 June, 2024;
originally announced June 2024.
-
Spectral form factor in chaotic, localized, and integrable open quantum many-body systems
Authors:
Jiachen Li,
Stephen Yan,
Tomaž Prosen,
Amos Chan
Abstract:
We numerically study the spectral statistics of open quantum many-body systems (OQMBS) as signatures of quantum chaos (or the lack thereof), using the dissipative spectral form factor (DSFF), a generalization of the spectral form factor to complex spectra. We show that the DSFF of chaotic OQMBS generically displays the $\textit{quadratic}$ ramp-plateau behaviour of the Ginibre ensemble from random…
▽ More
We numerically study the spectral statistics of open quantum many-body systems (OQMBS) as signatures of quantum chaos (or the lack thereof), using the dissipative spectral form factor (DSFF), a generalization of the spectral form factor to complex spectra. We show that the DSFF of chaotic OQMBS generically displays the $\textit{quadratic}$ ramp-plateau behaviour of the Ginibre ensemble from random matrix theory, in contrast to the linear ramp-plateau behaviour of the Gaussian ensemble in closed quantum systems. Furthermore, in the presence of many-body interactions, such RMT behaviour emerges only after a time scale $τ_{\mathrm{dev}}$, which generally increases with system size for sufficiently large system size, and can be identified as the non-Hermitian analogue of the $\textit{many-body Thouless time}$. The universality of the random matrix theory behavior is demonstrated by surveying twelve models of OQMBS, including random Kraus circuits (quantum channels) and random Lindbladians (Liouvillians) in several symmetry classes, as well as Lindbladians of paradigmatic models such as the Sachdev-Ye-Kitaev (SYK), XXZ, and the transverse field Ising models. We devise an unfolding and filtering procedure to remove variations of the averaged density of states which would otherwise hide the universal RMT-like signatures in the DSFF for chaotic OQMBS. Beyond chaotic OQMBS, we study the spectral statistics of non-chaotic OQMBS, specifically the integrable XX model and a system in the many-body localized (MBL) regime in the presence of dissipation, which exhibit DSFF behaviours distinct from the ramp-plateau behaviour of random matrix theory. Lastly, we study the DSFF of Lindbladians with the Hamiltonian term set to zero, i.e. only the jump operators are present, and demonstrate that the results of RMT universality and scaling of many-body Thouless time survive even without coherent evolution.
△ Less
Submitted 2 May, 2024;
originally announced May 2024.
-
Integrability of Goldilocks quantum cellular automata
Authors:
Logan E. Hillberry,
Lorenzo Piroli,
Eric Vernier,
Nicole Yunger Halpern,
Tomaž Prosen,
Lincoln D. Carr
Abstract:
Goldilocks quantum cellular automata (QCA) have been simulated on quantum hardware and produce emergent small-world correlation networks. In Goldilocks QCA, a single-qubit unitary is applied to each qubit in a one-dimensional chain subject to a balance constraint: a qubit is updated if its neighbors are in opposite basis states. Here, we prove that a subclass of Goldilocks QCA -- including the one…
▽ More
Goldilocks quantum cellular automata (QCA) have been simulated on quantum hardware and produce emergent small-world correlation networks. In Goldilocks QCA, a single-qubit unitary is applied to each qubit in a one-dimensional chain subject to a balance constraint: a qubit is updated if its neighbors are in opposite basis states. Here, we prove that a subclass of Goldilocks QCA -- including the one implemented experimentally -- map onto free fermions and therefore can be classically simulated efficiently. We support this claim with two independent proofs, one involving a Jordan--Wigner transformation and one mapping the integrable six-vertex model to QCA. We compute local conserved quantities of these QCA and predict experimentally measurable expectation values. These calculations can be applied to test large digital quantum computers against known solutions. In contrast, typical Goldilocks QCA have equilibration properties and quasienergy-level statistics that suggest nonintegrability. Still, the latter QCA conserve one quantity useful for error mitigation. Our work provides a parametric quantum circuit with tunable integrability properties with which to test quantum hardware.
△ Less
Submitted 3 April, 2024;
originally announced April 2024.
-
Loss-induced quantum information jet in an infinite temperature Hubbard chain
Authors:
Patrik Penc,
Cătălin Paşcu Moca,
Örs Legeza,
Tomaž Prosen,
Gergely Zaránd,
Miklós Antal Werner
Abstract:
Information propagation in the one-dimensional infinite temperature Hubbard model with a dissipative particle sink at the end of a semi-infinite chain is studied. In the strongly interacting limit, the two-site mutual information and the operator entanglement entropy exhibit a rich structure with two propagating information fronts and superimposed interference fringes. A classical reversible cellu…
▽ More
Information propagation in the one-dimensional infinite temperature Hubbard model with a dissipative particle sink at the end of a semi-infinite chain is studied. In the strongly interacting limit, the two-site mutual information and the operator entanglement entropy exhibit a rich structure with two propagating information fronts and superimposed interference fringes. A classical reversible cellular automaton model quantitatively captures the transport and the slow, classical part of the correlations, but fails to describe the rapidly propagating information jet. The fast quantum jet resembles coherent free particle propagation, with the accompanying long-ranged interference fringes that are exponentially damped by short-ranged spin correlations in the many-body background.
△ Less
Submitted 29 February, 2024;
originally announced February 2024.
-
Structural Stability Hypothesis of Dual Unitary Quantum Chaos
Authors:
Jonathon Riddell,
Curt von Keyserlingk,
Tomaž Prosen,
Bruno Bertini
Abstract:
Having spectral correlations that, over small enough energy scales, are described by random matrix theory is regarded as the most general defining feature of quantum chaotic systems as it applies in the many-body setting and away from any semiclassical limit. Although this property is extremely difficult to prove analytically for generic many-body systems, a rigorous proof has been achieved for du…
▽ More
Having spectral correlations that, over small enough energy scales, are described by random matrix theory is regarded as the most general defining feature of quantum chaotic systems as it applies in the many-body setting and away from any semiclassical limit. Although this property is extremely difficult to prove analytically for generic many-body systems, a rigorous proof has been achieved for dual-unitary circuits -- a special class of local quantum circuits that remain unitary upon swapping space and time. Here we consider the fate of this property when moving from dual-unitary to generic quantum circuits focussing on the \emph{spectral form factor}, i.e., the Fourier transform of the two-point correlation. We begin with a numerical survey that, in agreement with previous studies, suggests that there exists a finite region in parameter space where dual-unitary physics is stable and spectral correlations are still described by random matrix theory, although up to a maximal quasienergy scale. To explain these findings, we develop a perturbative expansion: it recovers the random matrix theory predictions, provided the terms occurring in perturbation theory obey a relatively simple set of assumptions. We then provide numerical evidence and a heuristic analytical argument supporting these assumptions.
△ Less
Submitted 27 August, 2024; v1 submitted 29 February, 2024;
originally announced February 2024.
-
The spectral boundary of the Asymmetric Simple Exclusion Process (ASEP) -- free fermions, Bethe ansatz and random matrix theory
Authors:
Goran Nakerst,
Tomaž Prosen,
Masudul Haque
Abstract:
In non-equilibrium statistical mechanics, the Asymmetric Simple Exclusion Process (ASEP) serves as a paradigmatic example. We investigate the spectral characteristics of the ASEP, focusing on the spectral boundary of its generator matrix. We examine finite ASEP chains of length $L$, under periodic (pbc) and open boundary conditions (obc). Notably, the spectral boundary exhibits $L$ spikes for pbc…
▽ More
In non-equilibrium statistical mechanics, the Asymmetric Simple Exclusion Process (ASEP) serves as a paradigmatic example. We investigate the spectral characteristics of the ASEP, focusing on the spectral boundary of its generator matrix. We examine finite ASEP chains of length $L$, under periodic (pbc) and open boundary conditions (obc). Notably, the spectral boundary exhibits $L$ spikes for pbc and $L+1$ spikes for obc. Treating the ASEP generator as an interacting non-Hermitian fermionic model, we extend the model to have tunable interaction. In the non-interacting case, the analytically computed many-body spectrum shows a spectral boundary with prominent spikes. For pbc, we use the coordinate Bethe ansatz to interpolate between the noninteracting case to the ASEP limit, and show that these spikes stem from clustering of Bethe roots. The robustness of the spikes in the spectral boundary is demonstrated by linking the ASEP generator to random matrices with trace correlations or, equivalently, random graphs with distinct cycle structures, both displaying similar spiked spectral boundaries.
△ Less
Submitted 1 February, 2024;
originally announced February 2024.
-
Boundary Chaos: Spectral Form Factor
Authors:
Felix Fritzsch,
Tomaž Prosen
Abstract:
Random matrix spectral correlations is a defining feature of quantum chaos. Here, we study such correlations in a minimal model of chaotic many-body quantum dynamics where interactions are confined to the system's boundary, dubbed \textit{boundary chaos}, in terms of the spectral form factor and its fluctuations. We exactly calculate the latter in the limit of large local Hilbert space dimension…
▽ More
Random matrix spectral correlations is a defining feature of quantum chaos. Here, we study such correlations in a minimal model of chaotic many-body quantum dynamics where interactions are confined to the system's boundary, dubbed \textit{boundary chaos}, in terms of the spectral form factor and its fluctuations. We exactly calculate the latter in the limit of large local Hilbert space dimension $q$ for different classes of random boundary interactions and find it to coincide with random matrix theory, possibly after a non-zero Thouless time. The latter effect is due to a drastic enhancement of the spectral form factor, when integer time and system size fulfill a resonance condition. We compare our semiclassical (large $q$) results with numerics at small local Hilbert space dimension ($q=2,3$) and observe qualitatively similar features as in the semiclassical regime.
△ Less
Submitted 14 December, 2023;
originally announced December 2023.
-
Dual symplectic classical circuits: An exactly solvable model of many-body chaos
Authors:
Alexios Christopoulos,
Andrea De Luca,
D L Kovrizhin,
Tomaž Prosen
Abstract:
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 correlati…
▽ More
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 different choices of observables.
△ Less
Submitted 24 January, 2024; v1 submitted 4 July, 2023;
originally announced July 2023.
-
Localised Dynamics in the Floquet Quantum East Model
Authors:
Bruno Bertini,
Pavel Kos,
Tomaz Prosen
Abstract:
We introduce and study the discrete-time version of the Quantum East model, an interacting quantum spin chain inspired by simple kinetically constrained models of classical glasses. Previous work has established that its continuous-time counterpart displays a disorder-free localisation transition signalled by the appearance of an exponentially large (in the volume) family of non-thermal, localised…
▽ More
We introduce and study the discrete-time version of the Quantum East model, an interacting quantum spin chain inspired by simple kinetically constrained models of classical glasses. Previous work has established that its continuous-time counterpart displays a disorder-free localisation transition signalled by the appearance of an exponentially large (in the volume) family of non-thermal, localised eigenstates. Here we combine analytical and numerical approaches to show that: i) The transition persists for discrete times, in fact, it is present for any finite value of the time step apart from a zero measure set; ii) It is directly detected by following the non-equilibrium dynamics of the fully polarised state. Our findings imply that the transition is currently observable in state-of-the-art platforms for digital quantum simulation.
△ Less
Submitted 24 February, 2024; v1 submitted 21 June, 2023;
originally announced June 2023.
-
Kardar-Parisi-Zhang scaling in the Hubbard model
Authors:
Cătălin Paşcu Moca,
Miklós Antal Werner,
Angelo Valli,
Gergely Zaránd,
Tomaž Prosen
Abstract:
We explore the Kardar-Parisi-Zhang (KPZ) scaling in the one-dimensional Hubbard model, which exhibits global $SU_c(2)\otimes SU_s(2)$ symmetry at half-filling, for the pseudo-charge and the total spin. We analyze dynamical scaling properties of high temperature charge and spin correlations and transport. At half-filling, we observe a clear KPZ scaling in both charge and spin sectors. Away from hal…
▽ More
We explore the Kardar-Parisi-Zhang (KPZ) scaling in the one-dimensional Hubbard model, which exhibits global $SU_c(2)\otimes SU_s(2)$ symmetry at half-filling, for the pseudo-charge and the total spin. We analyze dynamical scaling properties of high temperature charge and spin correlations and transport. At half-filling, we observe a clear KPZ scaling in both charge and spin sectors. Away from half-filling, the $SU_c(2)$ charge symmetry is reduced to $U_c(1)$, while the $SU_s(2)$ symmetry for the total spin is retained. Consequently, transport in the charge sector becomes ballistic, while KPZ scaling is preserved in the spin sector. These findings confirm the link between non-abelian symmetries and KPZ scaling in the presence of integrability. We study two settings of the model: one involving a quench from a bi-partitioned state asymptotically close to the $T=\infty$ equilibrium state of the system, and another where the system is coupled to two markovian reservoirs at the two edges of the chain.
△ Less
Submitted 20 June, 2023;
originally announced June 2023.
-
Quantum Chaos and Coherence: Random Parametric Quantum Channels
Authors:
Apollonas S. Matsoukas-Roubeas,
Tomaž Prosen,
Adolfo del Campo
Abstract:
The survival probability of an initial Coherent Gibbs State (CGS) is a natural extension of the Spectral Form Factor (SFF) to open quantum systems. To quantify the interplay between quantum chaos and decoherence away from the semi-classical limit, we investigate the relation of this generalized SFF with the corresponding $l_1$-norm of coherence. As a working example, we introduce Parametric Quantu…
▽ More
The survival probability of an initial Coherent Gibbs State (CGS) is a natural extension of the Spectral Form Factor (SFF) to open quantum systems. To quantify the interplay between quantum chaos and decoherence away from the semi-classical limit, we investigate the relation of this generalized SFF with the corresponding $l_1$-norm of coherence. As a working example, we introduce Parametric Quantum Channels (PQC), a discrete-time model of unitary evolution mixed with the effects of measurements or transient interactions with an environment. The Energy Dephasing (ED) dynamics arises as a specific case in the Markovian limit. We demonstrate our results in a series of random matrix models.
△ Less
Submitted 18 August, 2024; v1 submitted 30 May, 2023;
originally announced May 2023.
-
On two non-ergodic reversible cellular automata, one classical, the other quantum
Authors:
Tomaz Prosen
Abstract:
We propose and discuss two variants of kinetic particle models - cellular automata in 1+1 dimensions, which have some appeal due to their simplicity and intriguing properties which could warrant further research and applications. The first model is a deterministic and reversible automaton describing two species of quasiparticles: stable massless matter particles moving with velocity $\pm 1$ and un…
▽ More
We propose and discuss two variants of kinetic particle models - cellular automata in 1+1 dimensions, which have some appeal due to their simplicity and intriguing properties which could warrant further research and applications. The first model is a deterministic and reversible automaton describing two species of quasiparticles: stable massless matter particles moving with velocity $\pm 1$ and unstable, standing (zero velocity) field particles. We discuss two distinct continuity equations for three conserved charges of the model. While the first two charges and the corresponding currents have support three (3) lattice sites and represent a lattice analogue of conserved energy-momentum tensor, we find an additional conserved charge and current with support of nine (9) sites, implying non-ergodic behaviour and potentially signalling integrability of the model with a highly nested R-matrix structure. The second model represents a quantum (or stochastic) deformation of a recently introduced and studied charged hardpoint lattice gas, where particles of different binary charge ($\pm 1$) and binary velocity ($\pm 1$) can nontrivially mix upon elastic collisional scattering. We show that while the unitary evolution rule of this model does not satisfy the full Yang-Baxter equation, it still satisfies an intriguing related identity which gives birth to an infinite set of local conserved operators, the-so-called glider operators.
△ Less
Submitted 3 April, 2023;
originally announced April 2023.
-
Dynamical criticality of magnetization transfer in integrable spin chains
Authors:
Žiga Krajnik,
Johannes Schmidt,
Enej Ilievski,
Tomaž Prosen
Abstract:
Recent studies have found that fluctuations of magnetization transfer in integrable spin chains violate the central limit property. Here we revisit the problem of anomalous counting statistics in the Landau-Lifshitz field theory by specializing to two distinct anomalous regimes featuring a dynamical critical point. By performing optimized numerical simulations using an integrable space-time discre…
▽ More
Recent studies have found that fluctuations of magnetization transfer in integrable spin chains violate the central limit property. Here we revisit the problem of anomalous counting statistics in the Landau-Lifshitz field theory by specializing to two distinct anomalous regimes featuring a dynamical critical point. By performing optimized numerical simulations using an integrable space-time discretization we extract the algebraic growth exponents of time-dependent cumulants which attain their threshold values. The distinctly non-Gaussian statistics of magnetization transfer in the easy-axis regime is found to converge towards the universal distribution of charged single-file systems. At the isotropic point we infer a weakly non-Gaussian distribution, corroborating the view that superdiffusive spin transport in integrable spin chains does not belong to any known dynamical universality class.
△ Less
Submitted 22 December, 2023; v1 submitted 29 March, 2023;
originally announced March 2023.
-
Ergodicity breaking and deviation from Eigenstate Thermalisation in relativistic QFT
Authors:
Miha Srdinšek,
Tomaž Prosen,
Spyros Sotiriadis
Abstract:
The validity of the ergodic hypothesis in quantum systems can be rephrased in the form of the eigenstate thermalisation hypothesis (ETH), a set of statistical properties for the matrix elements of local observables in energy eigenstates, which is expected to hold in any ergodic system. We test the ETH in a nonintegrable model of relativistic quantum field theory (QFT) using the numerical method of…
▽ More
The validity of the ergodic hypothesis in quantum systems can be rephrased in the form of the eigenstate thermalisation hypothesis (ETH), a set of statistical properties for the matrix elements of local observables in energy eigenstates, which is expected to hold in any ergodic system. We test the ETH in a nonintegrable model of relativistic quantum field theory (QFT) using the numerical method of Hamiltonian truncation in combination with analytical arguments based on Lorentz symmetry and renormalisation group theory. We find that there is an infinite sequence of eigenstates with the characteristics of quantum many-body scars, that is, exceptional eigenstates with observable expectation values that lie far from thermal values, and we show that these states are one-quasiparticle states. We argue that in the thermodynamic limit the eigenstates cover the entire area between two diverging lines: the line of one-quasiparticle states, whose direction is dictated by relativistic kinematics, and the thermal average line. Our results suggest that the strong version of the ETH is violated in any relativistic QFT whose spectrum admits a quasiparticle description.
△ Less
Submitted 16 January, 2024; v1 submitted 27 March, 2023;
originally announced March 2023.
-
Full Eigenstate Thermalization via Free Cumulants in Quantum Lattice Systems
Authors:
Silvia Pappalardi,
Felix Fritzsch,
Tomaž Prosen
Abstract:
The Eigenstate-Thermalization-Hypothesis (ETH) has been established as the general framework to understand quantum statistical mechanics. Only recently has the attention been paid to so-called full ETH, which accounts for higher-order correlations among matrix elements, and that can be rationalized theoretically using the language of Free Probability. In this work, we perform the first numerical i…
▽ More
The Eigenstate-Thermalization-Hypothesis (ETH) has been established as the general framework to understand quantum statistical mechanics. Only recently has the attention been paid to so-called full ETH, which accounts for higher-order correlations among matrix elements, and that can be rationalized theoretically using the language of Free Probability. In this work, we perform the first numerical investigation of the full ETH in physical many-body systems with local interactions by testing the decomposition of higher-order correlators into thermal free cumulants for local operators. We perform exact diagonalization on two classes of local non-integrable (chaotic) quantum many-body systems: spin chain Hamiltonians and Floquet brickwork unitary circuits. We show that the dynamics of four-time correlation functions are encoded in fourth-order free cumulants, as predicted by ETH. Their dependence on frequency encodes the physical properties of local many-body systems and distinguishes them from structureless, rotationally invariant ensembles of random matrices.
△ Less
Submitted 2 September, 2024; v1 submitted 1 March, 2023;
originally announced March 2023.
-
Boundary Chaos: Exact Entanglement Dynamics
Authors:
Felix Fritzsch,
Roopayan Ghosh,
Tomaž Prosen
Abstract:
We compute the dynamics of entanglement in the minimal setup producing ergodic and mixing quantum many-body dynamics, which we previously dubbed {\em 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 convent…
▽ More
We compute the dynamics of entanglement in the minimal setup producing ergodic and mixing quantum many-body dynamics, which we previously dubbed {\em 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 cause entanglement to grow on extensive time scales proportional to system size.
△ Less
Submitted 1 June, 2023; v1 submitted 19 January, 2023;
originally announced January 2023.
-
Localization challenges quantum chaos in the finite two-dimensional Anderson model
Authors:
Jan Šuntajs,
Tomaž Prosen,
Lev Vidmar
Abstract:
It is believed that the two-dimensional (2D) Anderson model exhibits localization for any nonzero disorder in the thermodynamic limit and it is also well known that the finite-size effects are considerable in the weak disorder limit. Here we numerically study the quantum-chaos to localization transition in the finite 2D Anderson model using standard indicators used in the modern literature such as…
▽ More
It is believed that the two-dimensional (2D) Anderson model exhibits localization for any nonzero disorder in the thermodynamic limit and it is also well known that the finite-size effects are considerable in the weak disorder limit. Here we numerically study the quantum-chaos to localization transition in the finite 2D Anderson model using standard indicators used in the modern literature such as the level spacing ratio, spectral form factor, variances of observable matrix elements, participation entropy and the eigenstate entanglement entropy. We show that many features of these indicators may indicate emergence of robust single-particle quantum chaos at weak disorder. However, we argue that a careful numerical analysis is consistent with the single-parameter scaling theory and predicts the breakdown of quantum chaos at any nonzero disorder value in the thermodynamic limit. Among the hallmarks of this breakdown are the universal behavior of the spectral form factor at weak disorder, and the universal scaling of various indicators as a function of the parameter $u = \left(W \ln V\right)^{-1}$ where $W$ is the disorder strength and $V$ is the number of lattice sites.
△ Less
Submitted 26 February, 2023; v1 submitted 20 December, 2022;
originally announced December 2022.
-
Anomalous diffusion in the Long-Range Haken-Strobl-Reineker model
Authors:
Alberto Giuseppe Catalano,
Francesco Mattiotti,
Jérôme Dubail,
David Hagenmüller,
Tomaž Prosen,
Fabio Franchini,
Guido Pupillo
Abstract:
We analyze the propagation of excitons in a $d$-dimensional lattice with power-law hopping $\propto 1/r^α$ in the presence of dephasing, described by a generalized Haken-Strobl-Reineker model. We show that in the strong dephasing (quantum Zeno) regime the dynamics is described by a classical master equation for an exclusion process with long jumps. In this limit, we analytically compute the spatia…
▽ More
We analyze the propagation of excitons in a $d$-dimensional lattice with power-law hopping $\propto 1/r^α$ in the presence of dephasing, described by a generalized Haken-Strobl-Reineker model. We show that in the strong dephasing (quantum Zeno) regime the dynamics is described by a classical master equation for an exclusion process with long jumps. In this limit, we analytically compute the spatial distribution, whose shape changes at a critical value of the decay exponent $α_{\rm cr} = (d+2)/2$. The exciton always diffuses anomalously: a superdiffusive motion is associated to a Lévy stable distribution with long-range algebraic tails for $α\leqα_{\rm cr}$, while for $α> α_{\rm cr}$ the distribution corresponds to a surprising mixed Gaussian profile with long-range algebraic tails, leading to the coexistence of short-range diffusion and long-range Lévy-flights. In the many-exciton case, we demonstrate that, starting from a domain-wall exciton profile, algebraic tails appear in the distributions for any $α$, which affects thermalization: the longer the hopping range, the faster equilibrium is reached. Our results are directly relevant to experiments with cold trapped ions, Rydberg atoms and supramolecular dye aggregates. They provide a way to realize an exclusion process with long jumps experimentally.
△ Less
Submitted 23 August, 2023; v1 submitted 15 December, 2022;
originally announced December 2022.
-
Symmetry Classification of Many-Body Lindbladians: Tenfold Way and Beyond
Authors:
Lucas Sá,
Pedro Ribeiro,
Tomaž Prosen
Abstract:
We perform a systematic symmetry classification of many-body Lindblad superoperators describing general (interacting) open quantum systems coupled to a Markovian environment. Our classification is based on the behavior of the many-body Lindbladian under antiunitary symmetries and unitary involutions. We find that Hermiticity preservation reduces the number of symmetry classes, while trace preserva…
▽ More
We perform a systematic symmetry classification of many-body Lindblad superoperators describing general (interacting) open quantum systems coupled to a Markovian environment. Our classification is based on the behavior of the many-body Lindbladian under antiunitary symmetries and unitary involutions. We find that Hermiticity preservation reduces the number of symmetry classes, while trace preservation and complete positivity do not, and that the set of admissible classes depends on the presence of additional unitary symmetries: in their absence or in symmetry sectors containing steady states, many-body Lindbladians belong to one of ten non-Hermitian symmetry classes; if however, there are additional symmetries and we consider non-steady-state sectors, they belong to a different set of 19 classes. In both cases, it does not include classes with Kramer's degeneracy. Remarkably, our classification admits a straightforward generalization to the case of non-Markovian, and even non-trace-preserving, open quantum dynamics. While the abstract classification is completely general, we then apply it to general (long-range, interacting, spatially inhomogeneous) spin-$1/2$ chains. We explicitly build examples in all ten classes of Lindbladians in steady-state sectors, describing standard physical processes such as dephasing, spin injection and absorption, and incoherent hopping, thus illustrating the relevance of our classification for practical physics applications. Finally, we show that the examples in each class display unique random-matrix correlations. To fully resolve all symmetries, we employ the combined analysis of bulk complex spacing ratios and the overlap of eigenvector pairs related by symmetry operations. We further find that statistics of levels constrained onto the real and imaginary axes or close to the origin are not universal due to spontaneous breaking of Lindbladian PT symmetry.
△ Less
Submitted 22 August, 2023; v1 submitted 1 December, 2022;
originally announced December 2022.
-
Spectral and steady-state properties of fermionic random quadratic Liouvillians
Authors:
João Costa,
Pedro Ribeiro,
Andrea de Luca,
Tomaž Prosen,
Lucas Sá
Abstract:
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…
▽ More
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.
△ Less
Submitted 13 October, 2023; v1 submitted 14 October, 2022;
originally announced October 2022.
-
Universal anomalous fluctuations in charged single-file systems
Authors:
Žiga Krajnik,
Johannes Schmidt,
Vincent Pasquier,
Tomaž Prosen,
Enej Ilievski
Abstract:
Introducing a general class of one-dimensional single-file systems (meaning that particle crossings are prohibited) of interacting hardcore particles with internal degrees of freedom (called charge), we exhibit a novel type of dynamical universality reflected in anomalous statistical properties of macroscopic fluctuating observables such as charge transfer. We find that stringent dynamical constra…
▽ More
Introducing a general class of one-dimensional single-file systems (meaning that particle crossings are prohibited) of interacting hardcore particles with internal degrees of freedom (called charge), we exhibit a novel type of dynamical universality reflected in anomalous statistical properties of macroscopic fluctuating observables such as charge transfer. We find that stringent dynamical constraints lead to universal anomalous statistics of cumulative charge currents manifested both on the timescale characteristic of typical fluctuations and also in the rate function describing rare events. By computing the full counting statistics of net transferred charge between two extended subsystems, we establish a number of unorthodox dynamical properties in an analytic fashion. Most prominently, typical fluctuations in equilibrium are governed by a universal distribution that markedly deviates from the expected Gaussian statistics, whereas large fluctuations are described by an exotic large-deviation rate function featuring an exceptional triple critical point. Far from equilibrium, competition between dynamical phases leads to dynamical phase transitions of first and second order and spontaneous breaking of fluctuation symmetry of the univariate charge large-deviation function. The rich phenomenology of the outlined dynamical universality is exemplified on an exactly solvable classical cellular automaton of charged hardcore particles. We determine the dynamical phase diagram in the framework of Lee-Yang's theory of phase transitions and exhibit a hyper-dimensional diagram of distinct dynamical regimes. Our findings lead us to conclude that the conventional classification of dynamical universality classes based on the algebraic dynamical exponents and asymptotic scaling functions that characterize hydrodynamic relaxation of the dynamical structure factor is incomplete and calls for refinement.
△ Less
Submitted 1 February, 2024; v1 submitted 2 August, 2022;
originally announced August 2022.
-
Dual unitary circuits in random geometries
Authors:
Yusuf Kasim,
Tomaž Prosen
Abstract:
Recently introduced dual unitary brickwork circuits have been recognised as paradigmatic exactly solvable quantum chaotic many-body systems with tunable degree of ergodicity and mixing. Here we show that regularity of the circuit lattice is not crucial for exact solvability. We consider a circuit where random 2-qubit dual unitary gates sit at intersections of random arrangements of straight lines…
▽ More
Recently introduced dual unitary brickwork circuits have been recognised as paradigmatic exactly solvable quantum chaotic many-body systems with tunable degree of ergodicity and mixing. Here we show that regularity of the circuit lattice is not crucial for exact solvability. We consider a circuit where random 2-qubit dual unitary gates sit at intersections of random arrangements of straight lines in two dimensions (mikado) and analytically compute the variance of the spatio-temporal correlation function of local operators. Note that the average correlator vanishes due to local Haar randomness of the gates. The result can be physically motivated for two random mikado settings. The first corresponds to the thermal state of free particles carrying internal qubit degrees of freedom which experience interaction at kinematic crossings, while the second represents rotationally symmetric (random euclidean) space-time.
△ Less
Submitted 20 June, 2022;
originally announced June 2022.
-
Spectral form factor in a minimal bosonic model of many-body quantum chaos
Authors:
Dibyendu Roy,
Divij Mishra,
Tomaž Prosen
Abstract:
We study spectral form factor in periodically-kicked bosonic chains. We consider a family of models where a Hamiltonian with the terms diagonal in the Fock space basis, including random chemical potentials and pair-wise interactions, is kicked periodically by another Hamiltonian with nearest-neighbor hopping and pairing terms. We show that for intermediate-range interactions, random phase approxim…
▽ More
We study spectral form factor in periodically-kicked bosonic chains. We consider a family of models where a Hamiltonian with the terms diagonal in the Fock space basis, including random chemical potentials and pair-wise interactions, is kicked periodically by another Hamiltonian with nearest-neighbor hopping and pairing terms. We show that for intermediate-range interactions, random phase approximation can be used to rewrite the spectral form factor in terms of a bi-stochastic many-body process generated by an effective bosonic Hamiltonian. In the particle-number conserving case, i.e., when pairing terms are absent, the effective Hamiltonian has a non-abelian $SU(1,1)$ symmetry, resulting in universal quadratic scaling of the Thouless time with the system size, irrespective of the particle number. This is a consequence of degenerate symmetry multiplets of the subleading eigenvalue of the effective Hamiltonian and is broken by the pairing terms. In the latter case, we numerically find a nontrivial systematic system-size dependence of the Thouless time, in contrast to a related recent study for kicked fermionic chains.
△ Less
Submitted 29 August, 2022; v1 submitted 10 March, 2022;
originally announced March 2022.
-
Exact anomalous current fluctuations in a deterministic interacting model
Authors:
Žiga Krajnik,
Johannes Schmidt,
Vincent Pasquier,
Enej Ilievski,
Tomaž Prosen
Abstract:
We analytically compute the full counting statistics of charge transfer in a classical automaton of interacting charged particles. Deriving a closed-form expression for the moment generating function with respect to a stationary equilibrium state, we employ asymptotic analysis to infer the structure of charge current fluctuations for a continuous range of timescales. The solution exhibits several…
▽ More
We analytically compute the full counting statistics of charge transfer in a classical automaton of interacting charged particles. Deriving a closed-form expression for the moment generating function with respect to a stationary equilibrium state, we employ asymptotic analysis to infer the structure of charge current fluctuations for a continuous range of timescales. The solution exhibits several unorthodox features. Most prominently, on the timescale of typical fluctuations the probability distribution of the integrated charge current in a stationary ensemble without bias is distinctly non-Gaussian despite diffusive behavior of dynamical charge susceptibility. While inducing a charge imbalance is enough to recover Gaussian fluctuations, we find that higher cumulants grow indefinitely in time with different exponents, implying singular scaled cumulants. We associate this phenomenon with the lack of a regularity condition on moment generating functions and the onset of a dynamical critical point. In effect, the scaled cumulant generating function does not, irrespectively of charge bias, represent a faithful generating function of the scaled cumulants, yet the associated Legendre dual yields the correct large-deviation rate function. Our findings hint at novel types of dynamical universality classes in deterministic many-body systems.
△ Less
Submitted 4 April, 2022; v1 submitted 13 January, 2022;
originally announced January 2022.
-
Simulating Lindbladian evolution with non-abelian symmetries: Ballistic front propagation in the $SU(2)$ Hubbard model with a localized loss
Authors:
Cătălin Paşcu Moca,
Miklós Antal Werner,
Örs Legeza,
Tomaž Prosen,
Márton Kormos,
Gergely Zaránd
Abstract:
We develop a non-Abelian time evolving block decimation (NA-TEBD) approach to study of open systems governed by Lindbladian time evolution, while exploiting an arbitrary number of abelian or non-abelian symmetries. We illustrate this method in a one-dimensional fermionic $SU(2)$ Hubbard model on a semi-infinite lattice with localized particle loss at one end. We observe a ballistic front propagati…
▽ More
We develop a non-Abelian time evolving block decimation (NA-TEBD) approach to study of open systems governed by Lindbladian time evolution, while exploiting an arbitrary number of abelian or non-abelian symmetries. We illustrate this method in a one-dimensional fermionic $SU(2)$ Hubbard model on a semi-infinite lattice with localized particle loss at one end. We observe a ballistic front propagation with strongly renormalized front velocity, and a hydrodynamic current density profile. For large loss rates, a suppression of the particle current is observed, as a result of the quantum Zeno effect. Operator entanglement is found to propagate faster than the depletion profile, preceding the latter.
△ Less
Submitted 31 December, 2021;
originally announced December 2021.
-
Lindbladian dissipation of strongly-correlated quantum matter
Authors:
Lucas Sá,
Pedro Ribeiro,
Tomaž Prosen
Abstract:
We propose the Sachdev-Ye-Kitaev Lindbladian as a paradigmatic solvable model of dissipative many-body quantum chaos. It describes $N$ strongly coupled Majorana fermions with random all-to-all interactions, with unitary evolution given by a quartic Hamiltonian and the coupling to the environment described by $M$ quadratic jump operators, rendering the full Lindbladian quartic in the Majorana opera…
▽ More
We propose the Sachdev-Ye-Kitaev Lindbladian as a paradigmatic solvable model of dissipative many-body quantum chaos. It describes $N$ strongly coupled Majorana fermions with random all-to-all interactions, with unitary evolution given by a quartic Hamiltonian and the coupling to the environment described by $M$ quadratic jump operators, rendering the full Lindbladian quartic in the Majorana operators. Analytical progress is possible by developing a dynamical mean-field theory for the Liouvillian time evolution on the Keldysh contour. By disorder-averaging the interactions, we derive an (exact) effective action for two collective fields (Green's function and self-energy). In the large-$N$, large-$M$ limit, we obtain the saddle-point equations satisfied by the collective fields, which determine the typical timescales of the dissipative evolution, particularly, the spectral gap that rules the relaxation of the system to its steady state. We solve the saddle-point equations numerically and find that, for strong or intermediate dissipation, the system relaxes exponentially, with a spectral gap that can be computed analytically, while for weak dissipation, there are oscillatory corrections to the exponential relaxation. In this letter, we illustrate the feasibility of analytical calculations in strongly correlated dissipative quantum matter.
△ Less
Submitted 29 June, 2022; v1 submitted 22 December, 2021;
originally announced December 2021.
-
Boundary driven $XYZ$ chain: Exact inhomogeneous triangular matrix product ansatz
Authors:
Vladislav Popkov,
Xin Zhang,
Tomaz Prosen
Abstract:
We construct an explicit matrix product ansatz for the steady state of a boundary driven $XY\!Z$ spin-$\tfrac{1}{2}$ chain for arbitrary local polarizing channels at the chain's ends. The ansatz, where the Lax operators are written explicitly in terms of infinite-dimensional bidiagonal (triangular) site-dependent matrices, becomes exact either in the (Zeno) limit of infinite dissipation strength,…
▽ More
We construct an explicit matrix product ansatz for the steady state of a boundary driven $XY\!Z$ spin-$\tfrac{1}{2}$ chain for arbitrary local polarizing channels at the chain's ends. The ansatz, where the Lax operators are written explicitly in terms of infinite-dimensional bidiagonal (triangular) site-dependent matrices, becomes exact either in the (Zeno) limit of infinite dissipation strength, or thermodynamic limit of infinite chain length. The solution is based on an extension of the newly discovered family of separable eigenstates of the model.
△ Less
Submitted 10 December, 2021;
originally announced December 2021.
-
Boundary Chaos
Authors:
Felix Fritzsch,
Tomaž Prosen
Abstract:
Scrambling in many-body quantum systems causes initially local observables to spread uniformly over the whole available Hilbert space under unitary dynamics, which in lattice systems causes exponential suppression of dynamical correlation functions with system size. Here, we present a perturbed free quantum circuit model, in which ergodicity is induced by an impurity interaction placed on the syst…
▽ More
Scrambling in many-body quantum systems causes initially local observables to spread uniformly over the whole available Hilbert space under unitary dynamics, which in lattice systems causes exponential suppression of dynamical correlation functions with system size. Here, we present a perturbed free quantum circuit model, in which ergodicity is induced by an impurity interaction placed on the system's boundary, that allows for demonstrating the underlying mechanism. This is achieved by mapping dynamical correlation functions of local operators acting at the boundary to a partition function with complex weights defined on a two dimensional lattice with a helical topology. We evaluate this partition function in terms of transfer matrices, which allow for numerically treating system sizes far beyond what is accessible by exact diagonalization and whose spectral properties determine the asymptotic scaling of correlations. Combining analytical arguments with numerical results we show that for impurities which remain unitary under partial transpose correlations are exponentially suppressed with system size in a particular scaling limit. In contrast for generic impurities or generic locations of the local operators correlations show persistent revivals with a period given by the system size.
△ Less
Submitted 16 June, 2022; v1 submitted 9 December, 2021;
originally announced December 2021.
-
Exact Spectral Statistics in Strongly Localised Circuits
Authors:
Bruno Bertini,
Pavel Kos,
Tomaz Prosen
Abstract:
Since the seminal work of Anderson, localisation has been recognised as a standard mechanism allowing quantum many-body systems to escape ergodicity. This idea acquired even more prominence in the last decade as it has been argued that localisation -- dubbed many-body localisation (MBL) in this context -- can sometimes survive local interactions in the presence of sufficiently strong disorder. A c…
▽ More
Since the seminal work of Anderson, localisation has been recognised as a standard mechanism allowing quantum many-body systems to escape ergodicity. This idea acquired even more prominence in the last decade as it has been argued that localisation -- dubbed many-body localisation (MBL) in this context -- can sometimes survive local interactions in the presence of sufficiently strong disorder. A conventional signature of localisation is in the statistical properties of the spectrum -- spectral statistics -- which differ qualitatively from those in the ergodic phase. Although features of the spectral statistics are routinely used as numerical diagnostics for localisation, they have never been derived from first principles in the presence of non-trivial interactions. Here we fill this gap and provide the example of a simple class of quantum many-body systems -- which we dub strongly localised quantum circuits -- that are interacting, localised, and where the spectral statistics can be characterised exactly. Furthermore, we show that these systems exhibit a cascade of three different regimes for spectral correlations depending on the energy scale: at small, intermediate, and large scales they behave as disconnected patches of three decreasing sizes. We argue that these features appear in generic MBL systems, with the difference that only at the smallest scale they do become Poissonian.
△ Less
Submitted 28 March, 2022; v1 submitted 29 October, 2021;
originally announced October 2021.
-
Exact solution of the Rule 150 reversible cellular automaton
Authors:
Joseph W. P. Wilkinson,
Tomaž Prosen,
Juan P. Garrahan
Abstract:
We study the dynamics of the Rule 150 reversible cellular automaton (RCA). This is a one-dimensional lattice system of binary variables with synchronous (Floquet) dynamics, corresponding to a bulk deterministic and reversible discrete version of the kinetically constrained XOR-Fredrickson-Andersen model, whereby the local dynamics is restricted: a site flips if and only if the states of its neighb…
▽ More
We study the dynamics of the Rule 150 reversible cellular automaton (RCA). This is a one-dimensional lattice system of binary variables with synchronous (Floquet) dynamics, corresponding to a bulk deterministic and reversible discrete version of the kinetically constrained XOR-Fredrickson-Andersen model, whereby the local dynamics is restricted: a site flips if and only if the states of its neighbouring sites are different from each other. Like other RCA which have been studied recently, such as Rule 54 and Rule 201, Rule 150 is integrable, however, in contrast is noninteracting. In particular, the emergent quasiparticles - the domain walls - behave as free fermions. This then allows us to solve the model by means of matrix product ansätze. We find the exact equilibrium and nonequilibrium stationary states for systems with closed (periodic) and open (stochastic) boundaries, respectively, resolve the full spectrum of the time evolution operator and, therefore, gain access to the relaxation dynamics, and obtain the exact large deviation statistics of dynamical observables in the long time limit.
△ Less
Submitted 3 December, 2021; v1 submitted 28 October, 2021;
originally announced October 2021.
-
Quantum chaos in triangular billiards
Authors:
Črt Lozej,
Giulio Casati,
Tomaž Prosen
Abstract:
We present an extensive numerical study of spectral statistics and eigenfunctions of quantized triangular billiards. We compute two million consecutive eigenvalues for six representative cases of triangular billiards, three with generic angles with irrational ratios with $π$, whose classical dynamics is presumably mixing, and three with exactly one angle rational with $π$, which are presumably onl…
▽ More
We present an extensive numerical study of spectral statistics and eigenfunctions of quantized triangular billiards. We compute two million consecutive eigenvalues for six representative cases of triangular billiards, three with generic angles with irrational ratios with $π$, whose classical dynamics is presumably mixing, and three with exactly one angle rational with $π$, which are presumably only weakly mixing or even only non-ergodic in case of right-triangles. We find excellent agreement of short and long range spectral statistics with the Gaussian orthogonal ensemble of random matrix theory for the most irrational generic triangle, while the other cases show small but significant deviations which are attributed either to scarring or super-scarring mechanism. This result, which extends the quantum chaos conjecture to systems with dynamical mixing in the absence of hard (Lyapunov) chaos, has been corroborated by analysing distributions of phase-space localisation measures of eigenstates and inspecting the structure of characteristic typical and atypical eigenfunctions.
△ Less
Submitted 22 January, 2022; v1 submitted 8 October, 2021;
originally announced October 2021.
-
Absence of Normal Fluctuations in an Integrable Magnet
Authors:
Žiga Krajnik,
Enej Ilievski,
Tomaž Prosen
Abstract:
We investigate dynamical fluctuations of transferred magnetization in the one-dimensional lattice Landau--Lifshitz magnet with uniaxial anisotropy, representing an emblematic model of interacting spins. We demonstrate that the structure of fluctuations in thermal equilibrium depends radically on the characteristic dynamical scale. In the ballistic regime, typical fluctuations are found to follow a…
▽ More
We investigate dynamical fluctuations of transferred magnetization in the one-dimensional lattice Landau--Lifshitz magnet with uniaxial anisotropy, representing an emblematic model of interacting spins. We demonstrate that the structure of fluctuations in thermal equilibrium depends radically on the characteristic dynamical scale. In the ballistic regime, typical fluctuations are found to follow a normal distribution and scaled cumulants are finite. In stark contrast, on the diffusive and superdiffusive timescales, relevant respectively for the easy-axis and isotropic magnet at vanishing total magnetization, typical fluctuations are no longer Gaussian and, remarkably, scaled cumulants are divergent. The observed anomalous features disappear upon breaking integrability, suggesting that the absence of normal fluctuations is intimately tied to the presence of soliton modes. In a nonequilibrium setting of the isotropic magnet with weakly polarized step-profile initial state we find a slow drift of dynamical exponent from the superdiffusive towards the diffusive value.
△ Less
Submitted 30 December, 2021; v1 submitted 27 September, 2021;
originally announced September 2021.
-
On two reversible cellular automata with two particle species
Authors:
Katja Klobas,
Tomaž Prosen
Abstract:
We introduce a pair of time-reversible models defined on the discrete space-time lattice with 3 states per site, specifically, a vacancy and a particle of two flavours (species). The local update rules reproduce the rule 54 reversible cellular automaton when only a single species of particles is present, and satisfy the requirements of flavour exchange (C), space-reversal (P), and time-reversal (T…
▽ More
We introduce a pair of time-reversible models defined on the discrete space-time lattice with 3 states per site, specifically, a vacancy and a particle of two flavours (species). The local update rules reproduce the rule 54 reversible cellular automaton when only a single species of particles is present, and satisfy the requirements of flavour exchange (C), space-reversal (P), and time-reversal (T) symmetries. We find closed-form expressions for three local conserved charges and provide an explicit matrix product form of the grand canonical Gibbs states, which are identical for both models. For one of the models this family of Gibbs states seems to be a complete characterisation of equilibrium (i.e. space and time translation invariant) states, while for the other model we empirically find a sequence of local conserved charges, one for each support size larger than 2, hinting to its algebraic integrability. Finally, we numerically investigate the behaviour of spatio-temporal correlation functions of charge densities, and test the hydrodynamic prediction for the model with exactly three local charges. Surprisingly, the numerically observed 'sound velocity' does not match the hydrodynamic value. The deviations are either significant, or they decay extremely slowly with the simulation time, which leaves us with an open question for the mechanism of such a glassy behaviour in a deterministic locally interacting system.
△ Less
Submitted 19 November, 2021; v1 submitted 3 September, 2021;
originally announced September 2021.
-
Thermalisation Dynamics and Spectral Statistics of Extended Systems with Thermalising Boundaries
Authors:
Pavel Kos,
Tomaz Prosen,
Bruno Bertini
Abstract:
We study thermalisation and spectral properties of extended systems connected, through their boundaries, to a thermalising Markovian bath. Specifically, we consider periodically driven systems modelled by brickwork quantum circuits where a finite section (block) of the circuit is constituted by arbitrary local unitary gates while its complement, which plays the role of the bath, is dual-unitary. W…
▽ More
We study thermalisation and spectral properties of extended systems connected, through their boundaries, to a thermalising Markovian bath. Specifically, we consider periodically driven systems modelled by brickwork quantum circuits where a finite section (block) of the circuit is constituted by arbitrary local unitary gates while its complement, which plays the role of the bath, is dual-unitary. We show that the evolution of local observables and the spectral form factor are determined by the same quantum channel, which we use to characterise the system's dynamics and spectral properties. In particular, we identify a family of fine-tuned quantum circuits -- which we call strongly non-ergodic -- that fails to thermalise even in this controlled setting, and, accordingly, their spectral form factor does not follow the random matrix theory prediction. We provide a set of necessary conditions on the local quantum gates that lead to strong non-ergodicity, and in the case of qubits, we provide a complete classification of strongly non-ergodic circuits. We also study the opposite extreme case of circuits that are almost dual-unitary, i.e., where thermalisation occurs with the fastest possible rate. We show that, in these systems, local observables and spectral form factor approach respectively thermal values and random matrix theory prediction exponentially fast. We provide a perturbative characterisation of the dynamics and, in particular, of the time-scale for thermalisation.
△ Less
Submitted 15 November, 2021; v1 submitted 17 August, 2021;
originally announced August 2021.
-
Absence of thermalization of free systems coupled to gapped interacting reservoirs
Authors:
Marko Ljubotina,
Dibyendu Roy,
Tomaž Prosen
Abstract:
We study the thermalization of a small $XX$ chain coupled to long, gapped $XXZ$ leads at either side by observing the relaxation dynamics of the whole system. Using extensive tensor network simulations, we show that such systems, although not integrable, appear to show either extremely slow thermalization or even lack thereof since the two can not be distinguished within the accuracy of our numeri…
▽ More
We study the thermalization of a small $XX$ chain coupled to long, gapped $XXZ$ leads at either side by observing the relaxation dynamics of the whole system. Using extensive tensor network simulations, we show that such systems, although not integrable, appear to show either extremely slow thermalization or even lack thereof since the two can not be distinguished within the accuracy of our numerics. We show that the persistent oscillations observed in the spin current in the middle of the $XX$ chain are related to eigenstates of the entire system located within the gap of the boundary chains. We find from exact diagonalization that some of these states remain strictly localized within the $XX$ chain and do not hybridize with the rest of the system. The frequencies of the persistent oscillations determined by numerical simulations of dynamics match the energy differences between these states exactly. This has important implications for open systems, where the strongly interacting leads are often assumed to thermalize the central system. Our results suggest that if we employ gapped systems for the leads, this assumption does not hold; this finding is particularly relevant to any potential future experimental studies of open quantum systems.
△ Less
Submitted 15 June, 2021;
originally announced June 2021.
-
Reversible Cellular Automata as Integrable Interactions Round-a-Face: Deterministic, Stochastic, and Quantized
Authors:
Tomaz Prosen
Abstract:
A family of reversible deterministic cellular automata, including the rules 54 and 201 of [Bobenko et al., Commun. Math. Phys. 158, 127 (1993)] as well as their kinetically constrained quantum (unitary) or stochastic deformations, is shown to correspond to integrable Floquet circuit models with local interactions round-a-face. Using inhomogeneous solutions of the star-triangle relation with a one…
▽ More
A family of reversible deterministic cellular automata, including the rules 54 and 201 of [Bobenko et al., Commun. Math. Phys. 158, 127 (1993)] as well as their kinetically constrained quantum (unitary) or stochastic deformations, is shown to correspond to integrable Floquet circuit models with local interactions round-a-face. Using inhomogeneous solutions of the star-triangle relation with a one or two dimensional spectral parameter, changing their functional form depending on the orientation, we provide an explicit construction of the transfer matrix and establish its conservation law and involutivity properties. Integrability is independently demonstrated by numerically exploring the spectral statistics via the Berry-Tabor conjecture. Curiously, we find that the deformed rule 54 model generically possesses no other local conserved quantities besides the net soliton current.
△ Less
Submitted 3 June, 2021; v1 submitted 2 June, 2021;
originally announced June 2021.
-
Many Body Quantum Chaos and Dual Unitarity Round-a-Face
Authors:
Tomaz Prosen
Abstract:
We propose a new type of locally interacting quantum circuits which are generated by unitary interactions round-a-face (IRF). Specifically, we discuss a set (or manifold) of dual-unitary IRFs with local Hilbert space dimension $d$ (DUIRF$(d)$) which generate unitary evolutions both in space and time directions of an extended 1+1 dimensional lattice. We show how arbitrary dynamical correlation func…
▽ More
We propose a new type of locally interacting quantum circuits which are generated by unitary interactions round-a-face (IRF). Specifically, we discuss a set (or manifold) of dual-unitary IRFs with local Hilbert space dimension $d$ (DUIRF$(d)$) which generate unitary evolutions both in space and time directions of an extended 1+1 dimensional lattice. We show how arbitrary dynamical correlation functions of local observables can be evaluated in terms of finite dimensional completely positive trace preserving unital maps, in complete analogy to recently studied circuits made of dual unitary brick gates (DUBG). In fact, we show that the simplest non-trivial (non-vanishing) local correlation functions in dual-unitary IRF circuits involve observables non-trivially supported on at least two sites. We completely characterise the 10-dimensional manifold of DUIRF$(2)$ for qubits ($d=2$) and provide, for $d=3,4,5,6,7$, empirical estimates of its dimensionality based on numerically determined dimensions of tangent spaces at an ensemble of random instances of dual-unitary IRF gates. In parallel, we apply the same algorithm to determine ${\rm dim}\,{\rm DUBG}(d)$ and show that they are of similar order though systematically larger than ${\rm dim}\,{\rm DUIRF}(d)$ for $d=2,3,4,5,6,7$. It is remarkable that both sets have rather complex topology for $d\ge 3$ in the sense that the dimension of the tangent space varies among different randomly generated points of the set. Finally, we provide additional data on dimensionality of the chiral extension of DUBG circuits with distinct local Hilbert spaces of dimensions $d\neq d'$ residing at even/odd lattice sites.
△ Less
Submitted 17 May, 2021;
originally announced May 2021.
-
Anisotropic Landau-Lifshitz Model in Discrete Space-Time
Authors:
Žiga Krajnik,
Enej Ilievski,
Tomaž Prosen,
Vincent Pasquier
Abstract:
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 demonstra…
▽ More
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.
△ Less
Submitted 8 November, 2021; v1 submitted 28 April, 2021;
originally announced April 2021.
-
Rule 54: Exactly solvable model of nonequilibrium statistical mechanics
Authors:
Berislav Buča,
Katja Klobas,
Tomaž Prosen
Abstract:
We review recent results on an exactly solvable model of nonequilibrium statistical mechanics, specifically the classical Rule 54 reversible cellular automaton and some of its quantum extensions. We discuss the exact microscopic description of nonequilibrium dynamics as well as the equilibrium and nonequilibrium stationary states. This allows us to obtain a rigorous handle on the corresponding eme…
▽ More
We review recent results on an exactly solvable model of nonequilibrium statistical mechanics, specifically the classical Rule 54 reversible cellular automaton and some of its quantum extensions. We discuss the exact microscopic description of nonequilibrium dynamics as well as the equilibrium and nonequilibrium stationary states. This allows us to obtain a rigorous handle on the corresponding emergent hydrodynamic description, which is treated as well. Specifically, we focus on two different paradigms of Rule 54 dynamics. Firstly, we consider a finite chain driven by stochastic boundaries, where we provide exact matrix product descriptions of the nonequilibrium steady state, most relevant decay modes, as well as the eigenvector of the tilted Markov chain yielding exact large deviations for a broad class of local and extensive observables. Secondly, we treat the explicit dynamics of macro-states on an infinite lattice and discuss exact closed form results for dynamical structure factor, multi-time-correlation functions and inhomogeneous quenches. Remarkably, these results prove that the model, despite its simplicity, behaves like a regular fluid with coexistence of ballistic (sound) and diffusive (heat) transport. Finally, we briefly discuss quantum interpretation of Rule 54 dynamics and explicit results on dynamical spreading of local operators and operator entanglement.
△ Less
Submitted 13 May, 2021; v1 submitted 30 March, 2021;
originally announced March 2021.
-
Eigenstate thermalization in dual-unitary quantum circuits: Asymptotics of spectral functions
Authors:
Felix Fritzsch,
Tomaž Prosen
Abstract:
The eigenstate thermalization hypothesis provides to date the most successful description of thermalization in isolated quantum systems by conjecturing statistical properties of matrix elements of typical operators in the (quasi-)energy eigenbasis. Here we study the distribution of matrix elements for a class of operators in dual-unitary quantum circuits in dependence of the frequency associated w…
▽ More
The eigenstate thermalization hypothesis provides to date the most successful description of thermalization in isolated quantum systems by conjecturing statistical properties of matrix elements of typical operators in the (quasi-)energy eigenbasis. Here we study the distribution of matrix elements for a class of operators in dual-unitary quantum circuits in dependence of the frequency associated with the corresponding eigenstates. We provide an exact asymptotic expression for the spectral function, i.e., the second moment of this frequency resolved distribution. The latter is obtained from the decay of dynamical correlations between local operators which can be computed exactly from the elementary building blocks of the dual-unitary circuits. Comparing the asymptotic expression with results obtained by exact diagonalization we find excellent agreement. Small fluctuations at finite system size are explicitly related to dynamical correlations at intermediate times and the deviations from their asymptotical dynamics. Moreover, we confirm the expected Gaussian distribution of the matrix elements by computing higher moments numerically.
△ Less
Submitted 18 May, 2021; v1 submitted 22 March, 2021;
originally announced March 2021.
-
Spectral properties of three-dimensional Anderson model
Authors:
Jan Šuntajs,
Tomaž Prosen,
Lev Vidmar
Abstract:
The three-dimensional Anderson model represents a paradigmatic model to understand the Anderson localization transition. In this work we first review some key results obtained for this model in the past 50 years, and then study its properties from the perspective of modern numerical approaches. Our main focus is on the quantitative comparison between the level sensitivity statistics and the level…
▽ More
The three-dimensional Anderson model represents a paradigmatic model to understand the Anderson localization transition. In this work we first review some key results obtained for this model in the past 50 years, and then study its properties from the perspective of modern numerical approaches. Our main focus is on the quantitative comparison between the level sensitivity statistics and the level statistics. While the former studies the sensitivity of Hamiltonian eigenlevels upon inserting a magnetic flux, the latter studies the properties of unperturbed eigenlevels. We define two versions of dimensionless conductance, the first corresponding to the width of the level curvature distribution relative to the mean level spacing, and the second corresponding to the ratio of the Heisenberg and the Thouless time obtained from the spectral form factor. We show that both conductances look remarkably similar around the localization transition, in particular, they predict a nearly identical critical point consistent with other measures of the transition. We then study some further properties of those quantities: for level curvatures, we discuss particular similarities and differences between the width of the level curvature distribution and the characteristic energy studied by Edwards and Thouless in their pioneering work [J. Phys. C. 5, 807 (1972)]. In the context of the spectral form factor, we show that at the critical point it enters a broad time-independent regime, in which its value is consistent with the level compressibility obtained from the level variance. Finally, we test the scaling solution of the average level spacing ratio in the crossover regime using the cost function minimization approach introduced in [Phys. Rev. B. 102, 064207 (2020)]. We find that the extracted transition point and the scaling coefficient agree with those from the literature to high numerical accuracy.
△ Less
Submitted 14 December, 2021; v1 submitted 9 March, 2021;
originally announced March 2021.
-
Spectral Statistics of Non-Hermitian Matrices and Dissipative Quantum Chaos
Authors:
Jiachen Li,
Tomaž Prosen,
Amos Chan
Abstract:
We propose a measure, which we call the dissipative spectral form factor (DSFF), to characterize the spectral statistics of non-Hermitian (and non-Unitary) matrices. We show that DSFF successfully diagnoses dissipative quantum chaos, and reveals correlations between real and imaginary parts of the complex eigenvalues up to arbitrary energy (and time) scale. Specifically, we provide the exact solut…
▽ More
We propose a measure, which we call the dissipative spectral form factor (DSFF), to characterize the spectral statistics of non-Hermitian (and non-Unitary) matrices. We show that DSFF successfully diagnoses dissipative quantum chaos, and reveals correlations between real and imaginary parts of the complex eigenvalues up to arbitrary energy (and time) scale. Specifically, we provide the exact solution of DSFF for the GinUE and for a Poissonian random spectrum (Poisson) as minimal models of dissipative quantum chaotic and integrable systems respectively. For dissipative quantum chaotic systems, we show that DSFF exhibits an exact rotational symmetry in its complex time argument $τ$. Analogous to the spectral form factor (SFF) behaviour for GUE, DSFF for GinUE shows a ``dip-ramp-plateau'' behavior in $|τ|$: DSFF initially decreases, increases at intermediate time scales, and saturates after a generalized Heisenberg time which scales as the inverse mean level spacing. Remarkably, for large matrix size, the ``ramp'' of DSFF for GinUE increases quadratically in $|τ|$, in contrast to the linear ramp in SFF for Hermitian ensembles. For dissipative quantum integrable systems, we show that DSFF takes a constant value except for a region in complex time whose size and behavior depends on the eigenvalue density. Numerically, we verify the above claims and show that DSFF for real and quaternion real Ginibre ensembles coincides with the GinUE behaviour except for a region in complex time plane of measure zero in the limit of large matrix size. As a physical example, we consider the quantum kicked top model with dissipation, and show that it falls under the Ginibre universality class and Poisson as the `kick' is switched on or off. Lastly, we study spectral statistics of ensembles of random classical stochastic matrices, and show that these models fall under the Ginibre universality class.
△ Less
Submitted 1 October, 2021; v1 submitted 8 March, 2021;
originally announced March 2021.
-
Random Matrix Spectral Form Factor of Dual-Unitary Quantum Circuits
Authors:
Bruno Bertini,
Pavel Kos,
Tomaz Prosen
Abstract:
We investigate a class of brickwork-like quantum circuits on chains of $d-$level systems (qudits) that share the so-called `dual unitarity' property. Namely, these systems generate unitary dynamics not only when propagating in the time direction, but also when propagating in the space direction. We consider space-time homogeneous (Floquet) circuits and perturb them with a quenched single-site diso…
▽ More
We investigate a class of brickwork-like quantum circuits on chains of $d-$level systems (qudits) that share the so-called `dual unitarity' property. Namely, these systems generate unitary dynamics not only when propagating in the time direction, but also when propagating in the space direction. We consider space-time homogeneous (Floquet) circuits and perturb them with a quenched single-site disorder, i.e. by applying independent single site random unitaries drawn from arbitrary non-singular distribution over ${\rm SU}(d)$, e.g. one concentrated around the identity, after each layer of the circuit. We identify the spectral form factor at time $t$ in the limit of long chains as the dimension of the commutant of a finite set of operators on a qudit ring of $t$ sites. For general dual unitary circuits of qubits $(d=2)$ and a family of their extensions to higher $d>2$, we provide explicit construction of the commutant and prove that spectral form factor exactly matches the prediction of circular unitary ensemble for all $t$, if only the local 2-qubit gates are different from a SWAP (non-interacting gate). We discuss and partly prove possible extensions of our results to a weaker (more singular) forms of disorder averaging, as well as to quantum circuits with time-reversal symmetry, and to computing higher moments of the spectral form factor.
△ Less
Submitted 16 July, 2021; v1 submitted 22 December, 2020;
originally announced December 2020.
-
Signatures of Chaos in Non-integrable Models of Quantum Field Theory
Authors:
Miha Srdinsek,
Tomaz Prosen,
Spyros Sotiriadis
Abstract:
We study signatures of quantum chaos in (1+1)D Quantum Field Theory (QFT) models. Our analysis is based on the method of Hamiltonian truncation, a numerical approach for the construction of low-energy spectra and eigenstates of QFTs that can be considered as perturbations of exactly solvable models. We focus on the double sine-Gordon, also studying the massive sine-Gordon and ${φ^4}$ model, all of…
▽ More
We study signatures of quantum chaos in (1+1)D Quantum Field Theory (QFT) models. Our analysis is based on the method of Hamiltonian truncation, a numerical approach for the construction of low-energy spectra and eigenstates of QFTs that can be considered as perturbations of exactly solvable models. We focus on the double sine-Gordon, also studying the massive sine-Gordon and ${φ^4}$ model, all of which are non-integrable and can be studied by this method with sufficiently high precision from small to intermediate perturbation strength. We analyze the statistics of level spacings and of eigenvector components, both of which are expected to follow Random Matrix Theory predictions. While level spacing statistics are close to the Gaussian Orthogonal Ensemble as expected, on the contrary, the eigenvector components follow a distribution markedly different from the expected Gaussian. Unlike in the typical quantum chaos scenario, the transition of level spacing statistics to chaotic behaviour takes place already in the perturbative regime. On the other hand, the distribution of eigenvector components does not appear to change or approach Gaussian behaviour, even for relatively large perturbations. Moreover, our results suggest that these features are independent of the choice of model and basis.
△ Less
Submitted 1 April, 2021; v1 submitted 15 December, 2020;
originally announced December 2020.
-
Integrable nonunitary open quantum circuits
Authors:
Lucas Sá,
Pedro Ribeiro,
Tomaž Prosen
Abstract:
We explicitly construct an integrable and strongly interacting dissipative quantum circuit via a trotterization of the Hubbard model with imaginary interaction strength. To prove integrability, we build an inhomogeneous transfer matrix, from which conserved superoperator charges can be derived, in particular, the circuit's dynamical generator. After showing the trace preservation and complete posi…
▽ More
We explicitly construct an integrable and strongly interacting dissipative quantum circuit via a trotterization of the Hubbard model with imaginary interaction strength. To prove integrability, we build an inhomogeneous transfer matrix, from which conserved superoperator charges can be derived, in particular, the circuit's dynamical generator. After showing the trace preservation and complete positivity of local maps, we reinterpret them as the Kraus representation of the local dynamics of free fermions with single-site dephasing. The integrability of the map is broken by adding interactions to the local coherent dynamics or by removing the dephasing. In particular, even circuits built from convex combinations of local free-fermion unitaries are nonintegrable. Moreover, the construction allows us to explicitly build circuits belonging to different non-Hermitian symmetry classes, which are characterized by the behavior under transposition instead of complex conjugation. We confirm all our analytical results by using complex spacing ratios to examine the spectral statistics of the dissipative circuits.
△ Less
Submitted 19 March, 2021; v1 submitted 12 November, 2020;
originally announced November 2020.
-
Chaos and Ergodicity in Extended Quantum Systems with Noisy Driving
Authors:
Pavel Kos,
Bruno Bertini,
Tomaž Prosen
Abstract:
We study the time evolution operator in a family of local quantum circuits with random fields in a fixed direction. We argue that the presence of quantum chaos implies that at large times the time evolution operator becomes effectively a random matrix in the many-body Hilbert space. To quantify this phenomenon we compute analytically the squared magnitude of the trace of the evolution operator --…
▽ More
We study the time evolution operator in a family of local quantum circuits with random fields in a fixed direction. We argue that the presence of quantum chaos implies that at large times the time evolution operator becomes effectively a random matrix in the many-body Hilbert space. To quantify this phenomenon we compute analytically the squared magnitude of the trace of the evolution operator -- the generalised spectral form factor -- and compare it with the prediction of Random Matrix Theory (RMT). We show that for the systems under consideration the generalised spectral form factor can be expressed in terms of dynamical correlation functions of local observables in the infinite temperature state, linking chaotic and ergodic properties of the systems. This also provides a connection between the many-body Thouless time $τ_{\rm th}$ -- the time at which the generalised spectral form factor starts following the random matrix theory prediction -- and the conservation laws of the system. Moreover, we explain different scalings of $τ_{\rm th}$ with the system size, observed for systems with and without the conservation laws.
△ Less
Submitted 18 March, 2021; v1 submitted 23 October, 2020;
originally announced October 2020.