Search

State-of-the-art quantum simulators permit local temporal control of interactions and midcircuit readout. These capabilities open the way towards the exploration of intriguing nonequilibrium phenomena. We illustrate this with a kinetically constrained many-body quantum system that has a natural implementation on Rydberg quantum simulators. The evolution proceeds in discrete time and is generated by repeatedly entangling the system with an auxiliary environment that is monitored and reset after each time-step. Despite featuring an uncorrelated infinite-temperature average stationary state, the dynamics displays coexistence of fast and slow space-time regions in stochastic realizations of the system state. The time-record of measurement outcomes on the environment serves as natural probe for such dynamical heterogeneity, which we characterize using tools from large deviation theory. Our work establishes the large deviation framework for discrete-time open quantum many-body systems as a means to characterize complex dynamics and collective phenomena in quantum processors and simulators., Comment: 8+6 pages, 3+2 figures, comments welcome

Using tensor network methods, we simulate the real-time evolution of the lattice Thirring model quenched out of equilibrium in both the critical and massive phases, and study the appearance of dynamical quantum phase transitions, as non-analyticities in the Loschmidt rate. Whereas the presence of a dynamical quantum phase transition in the model does not correspond to quenches across the critical line of the equilibrium phase diagram at zero temperature, we identify a threshold in the energy density of the initial state, necessary for a dynamical quantum phase transition to be present. Moreover, in the case of the gapped quench Hamiltonian, we unveil a connection of this threshold to a transition between different regions in the finite temperature phase diagram.

We present a numerical method to simulate a system of multiple emitters coupled to a non-interacting bath, in any parameter regime. Our method relies on a Block Lanczos transformation that maps the whole system onto a strip geometry, whose width is given by the number of emitters. Utilizing the spatial symmetries of the problem and identifying the relevant range of energies of the bath we achieve a more efficient simulation, which we perform using tensor network techniques. As a demonstration, we study the collective emission from multiple emitters coupled to a square lattice of bosons and observe how the departure from Markovianity as coupling strength and emitter number is increased prevents collective radiation. We also simulate the dynamic preparation of an excitation in a bound state from a multi-excitation initial state. Our work opens new possibilities for the systematic exploration of non-Markovian effects in the dynamics and equilibrium properties of multi-emitter systems. Furthermore, it can easily be extended to other setups, including finite bath temperature or impurities coupled to fermionic environments.

We study the non-equilibrium dynamics of the dissipative quantum East model via numerical tensor networks. We use matrix product states to represent evolution under quantum-jump unravellings for sizes beyond those accessible to exact diagonalisation. This allows us to demonstrate that dynamical heterogeneity accompanies slow relaxation, in analogy with what is seen in classical glassy systems. Furthermore, using variational matrix product operators we: (i) compute the spectral gap of the Lindbladian, and show that glassiness is enhanced in the presence of weak quantum fluctuations compared to the pure classical case, and (ii) obtain the dynamical large deviations by calculating the leading eigenvector of the tilted Lindbladian, and find clear evidence for a first-order active-inactive dynamical phase transition. We also show how to directly sample the rare quantum trajectories associated to the large deviations., Comment: 5+11 pages, 4+6 figures, added SM

We study the dynamics and thermalization of the Fredkin spin chain, a system with local three-body interactions, particle conservation and explicit kinetic constraints. We consider deformations away from its stochastic point in order to tune between regimes where kinetic energy dominates and those where potential energy does. By means of exact diagonalisation, perturbation theory and variational matrix product states, we show that the stochastic point is where a transition occurs between a phase of fast thermalization to one of slow metastable (prethermal) dynamics. This change in relaxation is connected to the emergence of additional kinetic constraints which lead to the fragmentation of Hilbert space in the limit of a large potential energy. We also show that this transition can lead to thermalization being evaded for special initial conditions due to non-thermal eigenstates (akin to quantum many-body scars). We provide clear evidence for the existence of these non-thermal states for large system sizes even when far from the large-potential-energy limit, and explain their connection to the emergent kinetic constraints., Comment: 11+3 pages, 9+1 figures

Tensor network methods have demonstrated their suitability for the study of equilibrium properties of lattice gauge theories, even close to the continuum limit. We use them in an out-of-equilibrium scenario, much less explored so far, by simulating the real-time collisions of composite mesons in the lattice Schwinger model. Constructing wave-packets of vector mesons at different incoming momenta, we observe the opening of the inelastic channel in which two heavier mesons are produced and identify the momentum threshold. To detect the products of the collision in the strong coupling regime we propose local quantitites that could be measured in current quantum simulation platforms., Comment: 11 pages, 6 figures

Nonstabilizerness, also known as ``magic'', stands as a crucial resource for achieving a potential advantage in quantum computing. Its connection to many-body physical phenomena is poorly understood at present, mostly due to a lack of practical methods to compute it at large scales. We present a novel approach for the evaluation of nonstabilizerness within the framework of matrix product states (MPS), based on expressing the MPS directly in the Pauli basis. Our framework provides a powerful tool for efficiently calculating various measures of nonstabilizerness, including stabilizer R\'enyi entropies, stabilizer nullity, and Bell magic, and enables the learning of the stabilizer group of an MPS. We showcase the efficacy and versatility of our method in the ground states of Ising and XXZ spin chains, as well as in circuits dynamics that has recently been realized in Rydberg atom arrays, where we provide concrete benchmarks for future experiments on logical qubits up to twice the sizes already realized., Comment: 7+10 pages, 3+7 figures

We introduce a quantum algorithm to efficiently prepare states with a small energy variance at the target energy. We achieve it by filtering a product state at the given energy with a Lorentzian filter of width $\delta$. Given a local Hamiltonian on $N$ qubits, we construct a parent Hamiltonian whose ground state corresponds to the filtered product state with variable energy variance proportional to $\delta\sqrt{N}$. We prove that the parent Hamiltonian is gapped and its ground state can be efficiently implemented in $\mathrm{poly}(N,1/\delta)$ time via adiabatic evolution. We numerically benchmark the algorithm for a particular non-integrable model and find that the adiabatic evolution time to prepare the filtered state with a width $\delta$ is independent of the system size $N$. Furthermore, the adiabatic evolution can be implemented with circuit depth $\mathcal{O}(N^2\delta^{-4})$. Our algorithm provides a way to study the finite energy regime of many body systems in quantum simulators by directly preparing a finite energy state, providing access to an approximation of the microcanonical properties at an arbitrary energy., Comment: 10 pages, 8 figures

Energy filter methods in combination with quantum simulation can efficiently access the properties of quantum many-body systems at finite energy densities [Lu et al. PRX Quantum 2, 020321 (2021)]. Classically simulating this algorithm with tensor networks can be used to investigate the microcanonical properties of large spin chains, as recently shown in [Yang et al. Phys. Rev. B 106, 024307 (2022)]. Here we extend this strategy to explore the properties of off-diagonal matrix elements of observables in the energy eigenbasis, fundamentally connected to the thermalization behavior and the eigenstate thermalization hypothesis. We test the method on integrable and non-integrable spin chains of up to 60 sites, much larger than accessible with exact diagonalization. Our results allow us to explore the scaling of the off-diagonal functions with the size and energy difference, and to establish quantitative differences between integrable and non-integrable cases., Comment: Accepted version, 16 pages, 8 figures

Many quantum algorithms rely on the measurement of complex quantum amplitudes. Standard approaches to obtain the phase information, such as the Hadamard test, give rise to large overheads due to the need for global controlled-unitary operations. We introduce a quantum algorithm based on complex analysis that overcomes this problem for amplitudes that are a continuous function of time. Our method only requires the implementation of real-time evolution and a shallow circuit that approximates a short imaginary-time evolution. We show that the method outperforms the Hadamard test in terms of circuit depth and that it is suitable for current noisy quantum computers when combined with a simple error-mitigation strategy.

In the out-of-equilibrium evolution induced by a quench, fast degrees of freedom generate long-range entanglement that is hard to encode with standard tensor networks. However, local observables only sense such long-range correlations through their contribution to the reduced local state as a mixture. We present a tensor network method that identifies such long-range entanglement and efficiently transforms it into mixture, much easier to represent. In this way, we obtain an effective description of the time-evolved state as a density matrix that captures the long-time behavior of local operators with finite computational resources., Comment: 5 pages, 4 figures, comments are welcome!

Quantum simulation is one of the most promising scientific applications of quantum computers. Due to decoherence and noise in current devices, it is however challenging to perform digital quantum simulation in a regime that is intractable with classical computers. In this work, we propose an experimental protocol for probing dynamics and equilibrium properties on near-term digital quantum computers. As a key ingredient of our work, we show that it is possible to study thermalization even with a relatively coarse Trotter decomposition of the Hamiltonian evolution of interest. Even though the step size is too large to permit a rigorous bound on the Trotter error, we observe that the system prethermalizes in accordance with previous results for Floquet systems. The dynamics closely resemble the thermalization of the model underlying the Trotterization up to long times. We extend the reach of our approach by developing an error mitigation scheme based on measurement and rescaling of survival probabilities. To demonstrate the effectiveness of the entire protocol, we apply it to the two-dimensional XY model and numerically verify its performance with realistic noise parameters for superconducting quantum devices. Our proposal thus provides a route to achieving quantum advantage for relevant problems in condensed matter physics.

We explore the rich nature of correlations in the ground state of ultracold atoms trapped in state-dependent optical lattices. In particular, we consider interacting fermionic ytterbium or strontium atoms, realizing a two-orbital Hubbard model with two spin components. We analyze the model in one-dimensional setting with the experimentally relevant hierarchy of tunneling and interaction amplitudes by means of exact diagonalization and matrix product states approaches, and study the correlation functions in density, spin, and orbital sectors as functions of variable densities of atoms in the ground and metastable excited states. We show that in certain ranges of densities these atomic systems demonstrate strong density-wave, ferro- and antiferromagnetic, as well as antiferroorbital correlations., Comment: 8 pages, 5 figures

We use projected entangled-pair states (PEPS) to calculate the large deviations (LD) statistics of the dynamical activity of the two dimensional East model, and the two dimensional symmetric simple exclusion process (SSEP) with open boundaries, in lattices of up to 40x40 sites. We show that at long-times both models have phase transitions between active and inactive dynamical phases. For the 2D East model we find that this trajectory transition is of the first-order, while for the SSEP we find indications of a second order transition. We then show how the PEPS can be used to implement a trajectory sampling scheme capable of directly accessing rare trajectories. We also discuss how the methods described here can be extended to study rare events at finite times., Comment: 6+7 pages, 4+9 figures, added most updated version and SM

We propose to use quantum information notions to characterize thermally induced melting of nonperturbative bound states at high temperatures. We apply tensor networks to investigate this idea in static and dynamical settings within the Ising quantum field theory, where bound states are confined fermion pairs - mesons. An equilibrium signature of meson melting is identified in the temperature dependence of the thermal-state second R\'enyi entropy, which varies from exponential to power-law scaling. Out of equilibrium, we identify as the relevant signature the transition from an oscillatory to a linear growing behavior of reflected entropy after a thermal quench. These analyses apply more broadly, which brings new ways of describing in-medium meson phenomena in quantum many-body and high-energy physics., Comment: minor extended discussions and layout changes

The quantum simulation of fermionic gauge field theories is one of the anticipated uses of quantum computers in the NISQ era. Recently work has been done to simulate properties of the fermionic Z2 gauge field theory in (1+1) D and the pure gauge theory in (2+1) D. In this work, we investigate various options for simulating the fermionic Z2 gauge field theory in (2+1) D. To simulate the theory on a NISQ device it is vital to minimize both the number of qubits used and the circuit depth. In this work we propose ways to optimize both criteria for simulating time dynamics. In particular, we develop a new way to simulate this theory on a quantum computer, with minimal qubit requirements. We provide a quantum circuit, simulating a single first order Trotter step, that minimizes the number of 2-qubit gates needed and gives comparable results to methods requiring more qubits. Furthermore, variational approaches are investigated that allow to further decrease the circuit depth., Comment: 11 pages, 3 figures. Improved numerical results, added references, fixed typos

Tensor networks provide extremely powerful tools for the study of complex classical and quantum many-body problems. Over the last two decades, the increment in the number of techniques and applications has been relentless, and especially the last ten years have seen an explosion of new ideas and results that may be overwhelming for the newcomer. This short review introduces the basic ideas, the best established methods and some of the most significant algorithmic developments that are expanding the boundaries of the tensor network potential. The goal is to help the reader not only appreciate the many possibilities offered by tensor networks, but also find their way through state-of-the-art codes, their applicability and some avenues of ongoing progress., Comment: Review paper. To appear (after edition) in Annu. Rev. Condens. Matter Phys. v2 fixes a problem with the placement of figures

We investigate quantum inspired algorithms to compute physical observables of quantum many-body systems at finite energies. They are based on the quantum algorithms proposed in [Lu et al. PRX Quantum 2, 020321 (2021)], which use the quantum simulation of the dynamics of such systems, as well as classical filtering and sampling techniques. Here, we replace the quantum simulation by standard classical methods based on matrix product states and operators. As a result, we can address significantly larger systems than those reachable by exact diagonalization or by other algorithms. We demonstrate the performance with spin chains up to 80 sites.

The transverse folding algorithm [Phys. Rev. Lett. 102, 240603] is a tensor network method to compute time-dependent local observables in out-of-equilibrium quantum spin chains that can sometimes overcome the limitations of matrix product states. We present a contraction strategy that makes use of the exact light cone structure of the tensor network representing the observables. The strategy can be combined with the hybrid truncation proposed for global quenches in [Phys. Rev. A 91, 032306], which significantly improves the efficiency of the method. We demonstrate the performance of this transverse light cone contraction also for transport coefficients, and discuss how it can be extended to other dynamical quantities., Comment: Added the data as ancillary files

Recent work has shown the effectiveness of tensor network methods for computing large deviation functions in constrained stochastic models in the infinite time limit. Here we show that these methods can also be used to study the statistics of dynamical observables at arbitrary finite time. This is a harder problem because, in contrast to the infinite time case where only the extremal eigenstate of a tilted Markov generator is relevant, for finite time the whole spectrum plays a role. We show that finite time dynamical partition sums can be computed efficiently and accurately in one dimension using matrix product states, and describe how to use such results to generate rare event trajectories on demand. We apply our methods to the Fredrickson-Andersen (FA) and East kinetically constrained models, and to the symmetric simple exclusion process (SSEP), unveiling dynamical phase diagrams in terms of counting field and trajectory time. We also discuss extensions of this method to higher dimensions., Comment: 6+2 pages, 3+2 figures

We introduce a Metropolis-Hastings Markov chain for Boltzmann distributions of classical spin systems. It relies on approximate tensor network contractions to propose correlated collective updates at each step of the evolution. We present benchmarks for a wide variety of instances of the two-dimensional Ising model, including ferromagnetic, antiferromagnetic, (fully) frustrated and Edwards-Anderson spin glass cases, and we show that, with modest computational effort, our Markov chain achieves sizeable acceptance rates, even in the vicinity of critical points. In each of the situations we have considered, the Markov chain compares well with other Monte Carlo schemes such as the Metropolis or Wolff algorithm: equilibration times appear to be reduced by a factor that varies between 40 and 2000, depending on the model and the observable being monitored. We also present an extension to three spatial dimensions, and demonstrate that it exhibits fast equilibration for finite ferro and antiferromagnetic instances. Additionally, and although it is originally designed for a square lattice of finite degrees of freedom with open boundary conditions, the proposed scheme can be used as such, or with slight modifications, to study triangular lattices, systems with continuous degrees of freedom, matrix models, a confined gas of hard spheres, or to deal with arbitrary boundary conditions., Comment: 21 pages, 25 figures

Using a recently introduced tensor network method, we study the density of states of the lattice Schwinger model, a standard testbench for lattice gauge theory numerical techniques, but also the object of recent experimental quantum simulations. We identify regimes of parameters where the spectrum appears to be symmetric and displays the expected continuum properties even for finite lattice spacing and number of sites. However, we find that for moderate system sizes and lattice spacing of $ga\sim O(1)$, the spectral density can exhibit very different properties with a highly asymmetric form. We also explore how the method can be exploited to extract thermodynamic quantities., Comment: 10 pages, 5 figures

The large deviation (LD) statistics of dynamical observables is encoded in the spectral properties of deformed Markov generators. Recent works have shown that tensor network methods are well suited to compute the relevant leading eigenvalues and eigenvectors accurately. However, the efficient generation of the corresponding rare trajectories is a harder task. Here we show how to exploit the MPS approximation of the dominant eigenvector to implement an efficient sampling scheme which closely resembles the optimal (so-called "Doob") dynamics that realises the rare events. We demonstrate our approach on three well-studied lattice models, the Fredrickson-Andersen and East kinetically constrained models (KCMs), and the symmetric simple exclusion process (SSEP). We discuss how to generalise our approach to higher dimensions., Comment: 13 pages, 3 figures

We propose the construction of thermodynamic ensembles that minimize the R\'enyi free energy, as an alternative to Gibbs states. For large systems, the local properties of these R\'enyi ensembles coincide with those of thermal equilibrium, and they can be used as approximations to thermal states. We provide algorithms to find tensor network approximations to the 2-R\'enyi ensemble. In particular, a matrix-product-state representation can be found by using gradient-based optimization on Riemannian manifolds, or via a non-linear evolution which yields the desired state as a fixed point. We analyze the performance of the algorithms and the properties of the ensembles on one-dimensional spin chains.

For an isolated generic quantum system out of equilibrium, the long time average of observables is given by the diagonal ensemble, i.e. the mixed state with the same probability for energy eigenstates as the initial state but without coherences between different energies. In this work we present a method to approximate the diagonal ensemble using tensor networks. Instead of simulating the real time evolution, we adapt a filtering scheme introduced earlier in [Phys. Rev. B 101, 144305 (2020)] to this problem. We analyze the performance of the method on a non-integrable spin chain, for which we observe that local observables converge towards thermal values polynomially with the inverse width of the filter., Comment: 9 pages, 9 figures

Variational minimization of tensor network states enables the exploration of low energy states of lattice gauge theories. However, the exact numerical evaluation of high-dimensional tensor network states remains challenging in general. In [E. Zohar, J. I. Cirac, Phys. Rev. D 97, 034510 (2018)] it was shown how, by combining gauged Gaussian projected entangled pair states with a variational Monte Carlo procedure, it is possible to efficiently compute physical observables. In this paper we demonstrate how this approach can be used to investigate numerically the ground state of a lattice gauge theory. More concretely, we explicitly carry out the variational Monte Carlo procedure based on such contraction methods for a pure gauge Kogut-Susskind Hamiltonian with a $\mathbb{Z}_3$ gauge field in two spatial dimensions. This is a first proof of principle to the method, which provides an inherent way to increase the number of variational parameters and can be readily extended to systems with physical fermions., Comment: 13 pages, 7 figures; v2: as published

We simulate a zero-temperature pure $\mathbb{Z}_3$ Lattice Gauge Theory in 2+1 dimensions by using an iPEPS (Infinite Projected Entangled-Pair State) ansatz for the ground state. Our results are therefore directly valid in the thermodynamic limit. They clearly show two distinct phases separated by a phase transition. We introduce an update strategy that enables plaquette terms and Gauss-law constraints to be applied as sequences of two-body operators. This allows the use of the most up-to-date iPEPS algorithms. From the calculation of spatial Wilson loops we are able to prove the existence of a confined phase. We show that with relatively low computational cost it is possible to reproduce crucial features of gauge theories. We expect that the strategy allows the extension of iPEPS studies to more general LGTs., Comment: 7 pages, 7 figures

We study a one-dimensional classical stochastic kinetically constrained model (KCM) inspired by Rydberg atoms in their "facilitated" regime, where sites can flip only if a single of their nearest neighbours is excited. We call this model "XOR-FA" to distinguish it from the standard Fredrickson-Andersen (FA) model. We describe the dynamics of the XOR-FA model, including its relation to simple exclusion processes in its domain wall representation. The interesting relaxation dynamics of the XOR-FA is related to the prominence of large dynamical fluctuations that lead to phase transitions between active and inactive dynamical phases as in other KCMs. By means of numerical tensor network methods we study in detail such transitions in the dynamical large deviation regime., Comment: 13+2 pages, 7+1 figures

We introduce two kinds of quantum algorithms to explore microcanonical and canonical properties of many-body systems. The first one is a hybrid quantum algorithm that, given an efficiently preparable state, computes expectation values in a finite energy interval around its mean energy. This algorithm is based on a filtering operator, similar to quantum phase estimation, which projects out energies outside the desired energy interval. However, instead of performing this operation on a physical state, it recovers the physical values by performing interferometric measurements without the need to prepare the filtered state. We show that the computational time scales polynomially with the number of qubits, the inverse of the prescribed variance, and the inverse error. In practice, the algorithm does not require the evolution for long times, but instead a significant number of measurements in order to obtain sensible results. Our second algorithm is a quantum-assisted Monte Carlo sampling method to compute other quantities which approach the expectation values for the microcanonical and canonical ensembles. Using classical Monte Carlo techniques and the quantum computer as a resource, this method circumvents the sign problem that is plaguing classical Quantum Monte Carlo simulations, as long as one can prepare states with suitable energies. All algorithms can be used with small quantum computers and analog quantum simulators, as long as they can perform the interferometric measurements. We also show that this last task can be greatly simplified at the expense of performing more measurements., Comment: 21 pages, 9 figures; v3: a new schematic and other minor improvements and corrections

The usual paradigm of open quantum systems falls short when the environment is actually coupled to additional fields or components that drive it out of equilibrium. Here we explore the simplest such scenario, by considering a two level system coupled to a first thermal reservoir that in turn couples to a second thermal bath at a different temperature. We derive a master equation description for the system and show that, in this situation, the dynamics can be especially rich. In particular, we observe prethermalization, a transitory phenomenon in which the system initially approaches thermal equilibrium with respect to the first reservoir, but after a longer time converges to the thermal state dictated by the temperature of the second environment. Using analytical arguments and numerical simulations, we analyze the occurrence of this phenomenon, and how it depends on temperatures and coupling strengths. The phenomenology gets even richer if the system is placed between two such non-equilibrium environments. In this case, the energy current through the system may exhibit transient features and even switch direction, before the system eventually reaches a non-equilibrium steady state., Comment: 14 pages, 8 figures

One of the most interesting directions in theoretical high-energy and condensed-matter physics is understanding dynamical properties of collective states of quantum field theories. The most elementary tool in this quest is retarded equilibrium correlators governing the linear response theory. In this article we examine tensor networks as a way of determining them in a fully ab initio way in a class of (1+1)-dimensional quantum field theories arising as infrared descriptions of quantum Ising chains.We show that, complemented with signal analysis using the Prony method, tensor network calculations for intermediate times provide a powerful way to explore the structure of singularities of the correlator in the complex frequency plane and to make predictions about the thermal response to perturbations in a class of nonintegrable interacting quantum field theories., Comment: v2: minor layout changes; added appendix C on ground state quench

We explore the relation between the entanglement of a pure state and its energy variance for a local one dimensional Hamiltonian, as the system size increases. In particular, we introduce a construction which creates a matrix product state of arbitrarily small energy variance $\delta^2$ for $N$ spins, with bond dimension scaling as $\sqrt{N} D_0^{1/\delta}$, where $D_0>1$ is a constant. This implies that a polynomially increasing bond dimension is enough to construct states with energy variance that vanishes with the inverse of the logarithm of the system size. We run numerical simulations to probe the construction on two different models, and compare the local reduced density matrices of the resulting states to the corresponding thermal equilibrium. Our results suggest that the spatially homogeneous states with logarithmically decreasing variance, which can be constructed efficiently, do converge to the thermal equilibrium in the thermodynamic limit, while the same is not true if the variance remains constant., Comment: small changes to fix typos and bibliographic references

We present concluding results from our study for zero-temperature phase structure of the massive Thirring model in 1+1 dimensions with staggered regularisation. Employing the method of matrix product states, several quantities, including two types of correlators, are investigated, leading to numerical evidence of a Berezinskii-Kosterlitz-Thouless phase transition. Exploratory results for real-time dynamics pertaining to this transition, obtained using the approaches of variational uniform matrix product state and time-dependent variational principle, are also discussed., Comment: 7 pages, 4 figures, contribution to the 37th International Symposium on Lattice Field Theory (LATTICE2019), 16-22 June 2019, Wuhan, China

We study in detail the properties of the quantum East model, an interacting quantum spin chain inspired by simple kinetically-constrained models of classical glasses. Through a combination of analytics, exact diagonalization and tensor-network methods we show the existence of a transition, from a fast to a slow thermalization regime, which manifests itself throughout the spectrum. On the slow side, by exploiting the localization of the ground state and the form of the Hamiltonian, we explicitly construct a large (exponential in size) number of non-thermal states which become exact finite-energy-density eigenstates in the large-size limit, as expected for a true phase transition. A ``super-spin'' generalization allows us to find a further large class of area-law states proved to display very slow relaxation. These states retain memory of their initial conditions for extremely long times. Our numerical analysis reveals that the localization properties are not limited to the ground state and that} many eigenstates have large overlap with product states and can be approximated well by matrix product states at arbitrary energy densities. The mechanism that induces localization to the ground state, and hence the non-thermal behavior of the system, can be extended to a wide range of models including a number of simple spin chains. We discuss implications of our results for slow thermalization and non-ergodicity more generally in disorder-free systems with constraints and we give numerical evidence that these results may be extended to two dimensional systems., Comment: 21 pages, 16 figures

Formulating gauge theories on a lattice offers a genuinely non-perturbative way of studying quantum field theories, and has led to impressive achievements. In particular, it significantly deepened our understanding of quantum chromodynamics. Yet, some very relevant problems remain inherently challenging, such as real time evolution, or the presence of a chemical potential, cases in which Monte Carlo simulations are hindered by a sign problem. In the last few years, a number of possible alternatives have been put forward, based on quantum information ideas, which could potentially open the access to areas of research that have so far eluded more standard methods. They include tensor network calculations, quantum simulations with different physical platforms and quantum computations, and constitute nowadays a vibrant research area. Experts from different fields, including experimental and theoretical high energy physics, condensed matter, and quantum information, are turning their attention to these interdisciplinary possibilities, and driving the progress of the field. The aim of this article is to review the status and perspectives of these new avenues for the exploration of lattice gauge theories., Comment: 28 pages (including 15 figures) plus bibliography. Comments are most welcome

We combine matrix product operator techniques with Chebyshev polynomial expansions and present a method that is able to explore spectral properties of quantum many-body Hamiltonians. In particular, we show how this method can be used to probe thermalization of large spin chains without explicitly simulating their time evolution, as well as to compute full and local densities of states. The performance is illustrated with the examples of the Ising and PXP spin chains. For the non-integrable Ising chain, our findings corroborate the presence of thermalization for several initial states, well beyond what direct time-dependent simulations have been able to achieve so far.

Employing matrix product states as an ansatz, we study the non-thermal phase structure of the (1+1)-dimensional massive Thirring model in the sector of vanishing total fermion number with staggered regularization. In this paper, details of the implementation for this project are described. To depict the phase diagram of the model, we examine the entanglement entropy, the fermion bilinear condensate and two types of correlation functions. Our investigation shows the existence of two phases, with one of them being critical and the other gapped. An interesting feature of the phase structure is that the theory with non-zero fermion mass can be conformal. We also find clear numerical evidence that these phases are separated by a transition of the Berezinskii-Kosterlitz-Thouless type. Results presented in this paper establish the possibility of using the matrix product states for probing this type of phase transition in quantum field theories. They can provide information for further exploration of scaling behaviour, and serve as an important ingredient for controlling the continuum extrapolation of the model., Comment: 31 pages, 18 figures; minor changes to the text, typos corrected, references added; version published in Physical Review D

Here we demonstrate that tensor network techniques - originally devised for the analysis of quantum many-body problems - are well suited for the detailed study of rare event statistics in kinetically constrained models (KCMs). As concrete examples we consider the Fredrickson-Andersen and East models, two paradigmatic KCMs relevant to the modelling of glasses. We show how variational matrix product states allow to numerically approximate - systematically and with high accuracy - the leading eigenstates of the tilted dynamical generators which encode the large deviation statistics of the dynamics. Via this approach we can study system sizes beyond what is possible with other methods, allowing us to characterise in detail the finite size scaling of the trajectory-space phase transition of these models, the behaviour of spectral gaps, and the spatial structure and "entanglement" properties of dynamical phases. We discuss the broader implications of our results., Comment: 15 pages (including supplemental material), 3+18 figures. Extended Supplemental material

Tensor Network States are ans\"atze for the efficient description of quantum many-body systems. Their success for one dimensional problems, together with the fact that they do not suffer from the sign problem and can address the simulation of real time evolution, have turned them into one of the most promising techniques to study strongly correlated systems. In the realm of Lattice Gauge Theories they can offer an alternative to standard lattice Monte Carlo calculations, which are suited for static properties and regimes where no sign problem appears. The application of Tensor Networks to this kind of problems is a young but rapidly evolving research field. This paper reviews some of the recent progress in this area, and how, using one dimensional models as testbench, some fundamental milestones have been reached that may pave the way to more ambitious goals., Comment: Contribution to the proceedings of the 36th Annual International Symposium on Lattice Field Theory - LATTICE2018, 22-28 July, 2018 Michigan State University, East Lansing, Michigan, USA. Bibliography updated in second version

We present preliminary results of a study on the non-thermal phase structure of the (1+1) dimensional massive Thirring model, employing the method of matrix product states. Through investigating the entanglement entropy, the fermion correlators and the chiral condensate, it is found that this approach enables us to observe numerical evidence of a Kosterlitz-Thouless phase transition in the model., Comment: 7 pages, 4 figures; contribution to the proceedings of Lattice 2018 conference

Building on a previously introduced block Lanczos method, we demonstrate how to approximate any operator function of the form Trf (A) when the argument A is given as a Hermitian matrix product operator. This gives access to quantities that, depending on the full spectrum, are difficult to access for standard tensor network techniques, such as the von Neumann entropy and the trace norm of an MPO. We present a modified, more efficient strategy for computing thermal properties of short- or long-range Hamiltonians, and illustrate the performance of the method with numerical results for the thermal equilibrium states of the Lipkin-Meshkov-Glick and Ising Hamiltonians., Comment: 8 pages, comments very welcome

We provide a detailed formulation of the recently proposed variational approach [Y. Ashida et al., Phys. Rev. Lett. 121, 026805 (2018)] to study ground-state properties and out-of-equilibrium dynamics for generic quantum spin-impurity systems. Motivated by the original ideas by Tomonaga, Lee, Low, and Pines, we construct a canonical transformation that completely decouples the impurity from the bath degrees of freedom. By combining this transformation with a Gaussian ansatz for the fermionic bath, we obtain a family of variational many-body states that can efficiently encode the strong entanglement between the impurity and fermions of the bath. We give a detailed derivation of equations of motions in the imaginary- and real-time evolutions on the variational manifold. We benchmark our approach by applying it to investigate ground-state and dynamical properties of the anisotropic Kondo model and compare results with those obtained using matrix-product state (MPS) ansatz. We show that our approach can achieve an accuracy comparable to MPS-based methods with several orders of magnitude fewer variational parameters than the corresponding MPS ansatz. Comparisons to the Yosida ansatz and the exact solution from the Bethe ansatz are also discussed. We use our approach to investigate the two-lead Kondo model and analyze its long-time spatiotemporal behavior and the conductance behavior at finite bias and magnetic fields. The obtained results are consistent with the previous findings in the Anderson model and the exact solutions at the Toulouse point., Comment: 17 pages, 10 figures

A versatile and efficient variational approach is developed to solve in- and out-of-equilibrium problems of generic quantum spin-impurity systems. Employing the discrete symmetry hidden in spin-impurity models, we present a new canonical transformation that completely decouples the impurity and bath degrees of freedom. Combining it with Gaussian states, we present a family of many-body states to efficiently encode nontrivial impurity-bath correlations. We demonstrate its successful application to the anisotropic and two-lead Kondo models by studying their spatiotemporal dynamics and universal behavior in the correlations, relaxation times and the differential conductance. We compare them to previous analytical and numerical results. In particular, we apply our method to study new types of nonequilibrium phenomena that have not been studied by other methods, such as long-time crossover in the ferromagnetic easy-plane Kondo model. The present approach will be applicable to a variety of unsolved problems in solid-state and ultracold-atomic systems., Comment: 6 pages, 4 figures

Tensor Network methods have been established as a powerful technique for simulating low dimensional strongly-correlated systems for over two decades. Employing the formalism of Matrix Product States, we investigate the phase diagram of the massive Thirring model. We also show the possibility of studying soliton dynamics and topological phase transition via the Thirring model., Comment: 8 pages, 3 figures; contribution to the proceedings of the 35th International Symposium on Lattice Field Theory, 18 - 24 June 2017, Granada, Spain

We construct almost conserved local operators, that possess a minimal commutator with the Hamiltonian of the system, near the many-body localization transition of a one-dimensional disordered spin chain. We collect statistics of these slow operators for different support sizes and disorder strengths, both using exact diagonalization and tensor networks. Our results show that the scaling of the average of the smallest commutators with the support size is sensitive to Griffiths effects in the thermal phase and the onset of many-body localization. Furthermore, we demonstrate that the probability distributions of the commutators can be analyzed using extreme value theory and that their tails reveal the difference between diffusive and sub-diffusive dynamics in the thermal phase., Comment: 9 pages, 7 figures

We propose an explicit formulation of the physical subspace for a (1+1)-dimensional SU(2) lattice gauge theory, where the gauge degrees of freedom are integrated out. Our formulation is completely general, and might be potentially suited for the design of future quantum simulators. Additionally, it allows for addressing the theory numerically with matrix product states. We apply this technique to explore the spectral properties of the model and the effect of truncating the gauge degrees of freedom to a small finite dimension. In particular, we determine the scaling exponents for the vector mass. Furthermore, we also compute the entanglement entropy in the ground state and study its scaling towards the continuum limit., Comment: 19 pages, 16 figures, version 2: published version

The study of lattice gauge theories with Monte Carlo simulations is hindered by the infamous sign problem that appears under certain circumstances, in particular at non-zero chemical potential. So far, there is no universal method to overcome this problem. However, recent years brought a new class of non-perturbative Hamiltonian techniques named tensor networks, where the sign problem is absent. In previous work, we have demonstrated that this approach, in particular matrix product states in 1+1 dimensions, can be used to perform precise calculations in a lattice gauge theory, the massless and massive Schwinger model. We have computed the mass spectrum of this theory, its thermal properties and real-time dynamics. In this work, we review these results and we extend our calculations to the case of two flavours and non-zero chemical potential. We are able to reliably reproduce known analytical results for this model, thus demonstrating that tensor networks can tackle the sign problem of a lattice gauge theory at finite density., Comment: 10 pages, 5 figures; proceedings of XIIth Quark Confinement and the Hadron Spectrum, Thessaloniki 2016

During recent years there has been an increasing interest in the application of matrix product states, and more generally tensor networks, to lattice gauge theories. This non-perturbative method is sign problem free and has already been successfully used to compute mass spectra, thermal states and phase diagrams, as well as real-time dynamics for Abelian and non-Abelian gauge models. In previous work we showed the suitability of the method to explore the zero-temperature phase structure of the multi-flavor Schwinger model at non-zero chemical potential, a regime where the conventional Monte Carlo approach suffers from the sign problem. Here we extend our numerical study by looking at the spatially resolved chiral condensate in the massless case. We recover spatial oscillations, similar to the theoretical predictions for the single-flavor case, with a chemical potential dependent frequency and an amplitude approximately given by the homogeneous zero density condensate value., Comment: 7 pages, 3 figures, proceedings of the 34th International Symposium on Lattice Field Theory, 24-30 July 2016, Southampton, UK; extension of the results for the chiral condensate from arXiv:1611.00705