Collective Monte Carlo updates through tensor network renormalization

Abstract

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.

21 pages, 25 figures

References in corpus (10)

• The density-matrix renormalization group in the age of matrix product states
• Tensor renormalization group approach to 2D classical lattice models
• Criticality, the area law, and the computational power of PEPS
• Minimally Entangled Typical Thermal State Algorithms
• Algorithms for finite Projected Entangled Pair States
• Variational quantum Monte Carlo simulations with tensor-network states
• Strings, Projected Entangled Pair States, and variational Monte Carlo methods
• Critical behavior of lattice Schwinger model with topological term at $\theta =\pi$ using Grassmann tensor renormalization group
• Phase Transitions of Ferromagnetic Potts Models on the Simple Cubic Lattice
• Do tensor renormalization group methods work for frustrated spin systems?

Cited by in corpus (7)

• Nonstabilizerness via matrix product states in the Pauli basis
• Critical line of the triangular Ising antiferromagnet in a field from a $C_3$-symmetric corner transfer matrix algorithm
• Efficient Representation of Minimally Entangled Typical Thermal States in two dimensions via Projected Entangled Pair States
• Tensor network Monte Carlo simulations for the two-dimensional random-bond Ising model
• Limitations of tensor network approaches for optimization and sampling: A comparison to quantum and classical Ising machines
• Hyperoptimized approximate contraction of tensor networks for rugged-energy-landscape spin glasses on periodic square and cubic lattices
• Phase transitions and crtical exponents in the six-vertex model on kagome lattices