Mari Carmen Bañuls

Strong and Weak Thermalization of Infinite Nonintegrable Quantum Systems

When a nonintegrable system evolves out of equilibrium for a long time, local observables are in general expected to attain stationary expectation values, independent of the details of the initial state. But the thermalization of a closed quantum system is not yet well understood. Here we show that it presents indeed a much richer phenomenology than its classical counterpart. Using a new numerical technique, we identify two distinct 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, shows convergence to thermal values 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.

Observation of Correlated Particle-Hole Pairs and String Order in Low-Dimensional Mott Insulators

Parity correlations in a one-dimensional Bose gas in an optical lattice reveal a hidden “string order.” Quantum phases of matter are characterized by the underlying correlations of the many-body system. Although this is typically captured by a local order parameter, it has been shown that a broad class of many-body systems possesses a hidden nonlocal order. In the case of bosonic Mott insulators, the ground state properties are governed by quantum fluctuations in the form of correlated particle-hole pairs that lead to the emergence of a nonlocal string order in one dimension. By using high-resolution imaging of low-dimensional quantum gases in an optical lattice, we directly detect these pairs with single-site and single-particle sensitivity and observe string order in the one-dimensional case.

The mass spectrum of the Schwinger model with Matrix Product States

A bstractWe show the feasibility of tensor network solutions for lattice gauge theories in Hamiltonian formulation by applying matrix product states algorithms to the Schwinger model with zero and non-vanishing fermion mass. We introduce new techniques to compute excitations in a system with open boundary conditions, and to identify the states corresponding to low momentum and different quantum numbers in the continuum. For the ground state and both the vector and scalar mass gaps in the massive case, the MPS technique attains precisions comparable to the best results available from other techniques.

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.

Matrix Product States for Dynamical Simulation of Infinite Chains

We propose a new method for computing the ground state properties and the time evolution of infinite chains based on a transverse contraction of the tensor network. The method does not require finite size extrapolation and avoids explicit truncation of the bond dimension along the evolution. By folding the network in the time direction prior to contraction, time-dependent expectation values and dynamic correlation functions can be computed after much longer evolution time than with any previous method. Moreover, the algorithm we propose can be used for the study of some noninvariant infinite chains, including impurity models.

Thermal evolution of the Schwinger model with matrix product operators

We demonstrate the suitability of tensor network techniques for describing the thermal evolution of lattice gauge theories. As a benchmark case, we have studied the temperature dependence of the chiral condensate in the Schwinger model, using matrix product operators to approximate the thermal equilibrium states for finite system sizes with non-zero lattice spacings. We show how these techniques allow for reliable extrapolations in bond dimension, step width, system size and lattice spacing, and for a systematic estimation and control of all error sources involved in the calculation. The reached values of the lattice spacing are small enough to capture the most challenging region of high temperatures and the final results are consistent with the analytical prediction by Sachs and Wipf over a broad temperature range.

Unifying projected entangled pair state contractions

The approximate contraction of a tensor network of projected entangled pair states (PEPS) is a fundamental ingredient of any PEPS algorithm, required for the optimization of the tensors in ground state search or time evolution, as well as for the evaluation of expectation values. An exact contraction is in general impossible, and the choice of the approximating procedure determines the efficiency and accuracy of the algorithm. We analyze different previous proposals for this approximation, and show that they can be understood via the form of their environment, i.e. the operator that results from contracting part of the network. This provides physical insight into the limitation of various approaches, and allows us to introduce a new strategy, based on the idea of clusters, that unifies previous methods. The resulting contraction algorithm interpolates naturally between the cheapest and most imprecise and the most costly and most precise method. We benchmark the different algorithms with finite PEPS, and show how the cluster strategy can be used for both the tensor optimization and the calculation of expectation values. Additionally, we discuss its applicability to the parallelization of PEPS and to infinite systems.

Adiabatic Preparation of a Heisenberg Antiferromagnet Using an Optical Superlattice

We analyze the possibility to prepare a Heisenberg antiferromagnet with cold fermions in optical lattices, starting from a band insulator and adiabatically changing the lattice potential. The numerical simulation of the dynamics in 1D allows us to identify the conditions for success, and to study the influence that the presence of holes in the initial state may have on the protocol. We also extend our results to two-dimensional systems.

Non-Abelian string breaking phenomena with matrix product states

A bstractUsing 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.

Algorithms for finite projected entangled pair states

Projected Entangled Pair States (PEPS) are a promising ansatz for the study of strongly correlated quantum many-body systems in two dimensions. But due to their high computational cost, developing and improving PEPS algorithms is necessary to make the ansatz widely usable in practice. Here we analyze several algorithmic aspects of the method. On the one hand, we quantify the connection between the correlation length of the PEPS and the accuracy of its approximate contraction, and discuss how purifications can be used in the latter. On the other, we present algorithmic improvements for the update of the tensor that introduce drastic gains in the numerical conditioning and the efficiency of the algorithms. Finally, the state-of-the-art general PEPS code is benchmarked with the Heisenberg and quantum Ising models on lattices of up to