Source author record

Mari Carmen Bañuls

Mari Carmen Bañuls appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2022arXiv

Finite time large deviations via matrix product states

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.

preprint2020arXiv

How much entanglement is needed to reduce the energy variance?

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 $δ^2$ for $N$ spins, with bond dimension scaling as $\sqrt{N} D_0^{1/δ}$, 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.

preprint2020arXiv

Prethermalization of quantum systems interacting with non-equilibrium environments

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.

preprint2020arXiv

Quantum East model: localization, non-thermal eigenstates and slow dynamics

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.

preprint2020arXiv

Simulating 2+1d $\mathbb{Z}_3$ lattice gauge theory with iPEPS

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.

preprint2019arXiv

Probing thermalization through spectral analysis with matrix product operators

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.

preprint2016arXiv

Chiral condensate in the Schwinger model with Matrix Product Operators

Tensor network (TN) methods, in particular the Matrix Product States (MPS) ansatz, have proven to be a useful tool in analyzing the properties of lattice gauge theories. They allow for a very good precision, much better than standard Monte Carlo (MC) techniques for the models that have been studied so far, due to the possibility of reaching much smaller lattice spacings. The real reason for the interest in the TN approach, however, is its ability, shown so far in several condensed matter models, to deal with theories which exhibit the notorious sign problem in MC simulations. This makes it prospective for dealing with the non-zero chemical potential in QCD and other lattice gauge theories, as well as with real-time simulations. In this paper, using matrix product operators, we extend our analysis of the Schwinger model at zero temperature to show the feasibility of this approach also at finite temperature. This is an important step on the way to deal with the sign problem of QCD. We analyze in detail the chiral symmetry breaking in the massless and massive cases and show that the method works very well and gives good control over a broad range of temperatures, essentially from zero to infinite temperature.

preprint2016arXiv

The multi-flavor Schwinger model with chemical potential - Overcoming the sign problem with Matrix Product States

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.

preprint2015arXiv

Non-Abelian string breaking phenomena with Matrix Product States

Using matrix product states, we explore numerically the phenomenology of string breaking in a non-Abelian lattice gauge theory, namely 1+1 dimensional SU(2). The technique allows us to study the static potential between external heavy charges, as traditionally explored by Monte Carlo simulations, but also to simulate the real-time dynamics of both static and dynamical fermions, as the latter are fully included in the formalism. We propose a number of observables that are sensitive to the presence or breaking of the flux string, and use them to detect and characterize the phenomenon in each of these setups.

preprint2015arXiv

Slowest local operators in quantum spin chains

We numerically construct slowly relaxing local operators in a nonintegrable spin-1/2 chain. Restricting the support of the operator to $M$ consecutive spins along the chain, we exhaustively search for the operator that minimizes the Frobenius norm of the commutator with the Hamiltonian. We first show that the Frobenius norm bounds the time scale of relaxation of the operator at high temperatures. We find operators with significantly slower relaxation than the slowest simple "hydrodynamic" mode due to energy diffusion. Then, we examine some properties of the nontrivial slow operators. Using both exhaustive search and tensor network techniques, we find similar slowly relaxing operators for a Floquet spin chain; this system is hydrodynamically "trivial", with no conservation laws restricting their dynamics. We argue that such slow relaxation may be a generic feature following from locality and unitarity.

preprint2015arXiv

Variational Matrix Product Operators for the Steady State of Dissipative Quantum Systems

We present a new variational method, based on the matrix product operator (MPO) ansatz, for finding the steady state of dissipative quantum chains governed by master equations of the Lindblad form. Instead of requiring an accurate representation of the system evolution until the stationary state is attained, the algorithm directly targets the final state, thus allowing for a faster convergence when the steady state is a MPO with small bond dimension. Our numerical simulations for several dissipative spin models over a wide range of parameters illustrate the performance of the method and show that indeed the stationary state is often well described by a MPO of very moderate dimensions.

preprint2014arXiv

The temperature dependence of the chiral condensate in the Schwinger model with Matrix Product States

We present our recent results for the tensor network (TN) approach to lattice gauge theories. TN methods provide an efficient approximation for quantum many-body states. We employ TN for one dimensional systems, Matrix Product States, to investigate the 1-flavour Schwinger model. In this study, we compute the chiral condensate at finite temperature. From the continuum extrapolation, we obtain the chiral condensate in the high temperature region consistent with the analytical calculation by Sachs and Wipf.

preprint2013arXiv

Matrix Product States for Lattice Field Theories

The term Tensor Network States (TNS) refers to a number of families of states that represent different ansätze for the efficient description of the state of a quantum many-body system. Matrix Product States (MPS) are one particular case of TNS, and have become the most precise tool for the numerical study of one dimensional quantum many-body systems, as the basis of the Density Matrix Renormalization Group method. Lattice Gauge Theories (LGT), in their Hamiltonian version, offer a challenging scenario for these techniques. While the dimensions and sizes of the systems amenable to TNS studies are still far from those achievable by 4-dimensional LGT tools, Tensor Networks can be readily used for problems which more standard techniques, such as Markov chain Monte Carlo simulations, cannot easily tackle. Examples of such problems are the presence of a chemical potential or out-of-equilibrium dynamics. We have explored the performance of Matrix Product States in the case of the Schwinger model, as a widely used testbench for lattice techniques. Using finite-size, open boundary MPS, we are able to determine the low energy states of the model in a fully non-perturbative manner and to accurately extract the mass spectrum. In this work we extend the analysis to the determination of the chiral condensate, both for massless and massive fermions. The method allows for accurate finite size and continuum limit extrapolations and produces remarkably precise results, thus showing the feasibility of these techniques for gauge theory problems.

preprint2012arXiv

Tensor network techniques for the computation of dynamical observables in 1D quantum spin systems

We analyze the recently developed folding algorithm [Phys. Rev. Lett. 102, 240603 (2009)] to simulate the dynamics of infinite quantum spin chains, and relate its performance to the kind of entanglement produced under the evolution of product states. We benchmark the accomplishments of this technique with respect to alternative strategies using Ising Hamiltonians with transverse and parallel fields, as well as XY models. Additionally, we evaluate its ability to find ground and thermal equilibrium states.

preprint2010arXiv

Strong and weak thermalization of infinite non-integrable quantum systems

When a non-integrable system evolves out of equilibrium for a long time, local observables are expected to attain stationary expectation values, independent of the details of the initial state. However, intriguing experimental results with ultracold gases have shown no thermalization in non-integrable settings, triggering an intense theoretical effort to decide the question. Here we show that the phenomenology of thermalization in a quantum system is much richer than its classical counterpart. Using a new numerical technique, we identify two distinct thermalization regimes, strong and weak, occurring for different initial states. Strong thermalization, intrinsically quantum, happens when instantaneous local expectation values converge to the thermal ones. Weak thermalization, well-known in classical systems, happens when local expectation values converge to the thermal ones only after time averaging. Remarkably, we find a third group of states showing no thermalization, neither strong nor weak, to the time scales one can reliably simulate.