Scaling Hypothesis for Matrix Product States

Tensor renormalization group for fermions

We review the basic ideas of the tensor renormalization group method and show how they can be applied for lattice field theory models involving relativistic fermions and Grassmann variables in arbitrary dimensions. We discuss recent progress for entanglement filtering, loop optimization, bond-weighting techniques and matrix product decompositions for Grassmann tensor networks. The new methods are tested with two-dimensional Wilson-Majorana fermions and multi-flavor Gross-Neveu models. We show that the methods can also be applied to the fermionic Hubbard model in 1+1 and 2+1 dimensions.

Emergent Conformal Boundaries from Finite-Entanglement Scaling in Matrix Product States

The use of finite entanglement scaling with matrix product states (MPS) has become a crucial tool for studying one-dimensional critical lattice theories, especially those with emergent conformal symmetry. We argue that finite entanglement introduces a relevant deformation in the critical theory. As a result, the bipartite entanglement Hamiltonian defined from the MPS can be understood as a boundary conformal field theory with a physical and an entanglement boundary. We are able to exploit the symmetry properties of the MPS to engineer the physical conformal boundary condition. The entanglement boundary, on the other hand, is related to the concrete lattice model and remains invariant under this relevant perturbation. Using critical lattice models described by the Ising, Potts, and free compact boson conformal field theories, we illustrate the influence of the symmetry and the relevant deformation on the conformal boundaries in the entanglement spectrum.

Finite-Entanglement Scaling of 2D Metals

We extend the study of finite-entanglement scaling from one-dimensional gapless models to two-dimensional systems with a Fermi surface. In particular, we show that the entanglement entropy of a contractible spatial region with linear size L scales as S∼Llog[ξf(L/ξ)] in the optimal tensor network, and hence area-law entangled, state approximation to a metallic state, where f(x) is a scaling function which depends on the shape of the Fermi surface and ξ is a finite correlation length induced by the restricted entanglement. Crucially, the scaling regime can be realized with numerically tractable bond dimensions. We also discuss the implications of the Lieb-Schultz-Mattis theorem at fractional filling for tensor network state approximations of metallic states.

Extracting the Speed of Light from Matrix Product States

We provide evidence that the spectrum of the local effective Hamiltonian and the transfer operator in infinite-system matrix product state simulations are identical up to a global rescaling factor, i.e., the speed of light of the system, when the underlying system is described by a 1+1 dimensional CFT. We provide arguments for this correspondence based on a path integral point of view. This observation turns out to yield very precise estimates for the speed of light in practice, confirming exact results to high precision where available, but also allowing us to finally determine the speed of light of the non-integrable, critical SU(2) Heisenberg chains with half-integer spin S>1/2 with unprecedented accuracy. We also show that the same technology applied to doped Hubbard ladders provides highly accurate velocities for a range of dopings. Combined with measurements of compressibilities we present new results for the Luttinger liquid parameter in the Luther-Emery regime of doped Hubbard ladders, outperforming earlier approaches based on the fitting of real-space correlation functions.

Contrasting pseudocriticality in the classical two-dimensional Heisenberg and RP^{2} models: Zero-temperature phase transition versus finite-temperature crossover

Tensor-network methods are used to perform a comparative study of the two-dimensional classical Heisenberg and RP^{2} models. We demonstrate that uniform matrix product states (MPSs) with explicit SO(3) symmetry can probe correlation lengths up to O(10^{3}) sites accurately, and we study the scaling of entanglement entropy and universal features of MPS entanglement spectra. For the Heisenberg model, we find no signs of a finite-temperature phase transition, supporting the scenario of asymptotic freedom. For the RP^{2} model we observe an abrupt onset of scaling behavior, consistent with hints of a finite-temperature phase transition reported in previous studies. A careful analysis of the softening of the correlation length divergence, the scaling of the entanglement entropy, and the MPS entanglement spectra shows that our results are inconsistent with true criticality, but are rather in agreement with the scenario of a crossover to a pseudocritical region which exhibits strong signatures of nematic quasi-long-range order at length scales below the true correlation length. Our results reveal a fundamental difference in scaling behavior between the Heisenberg and RP^{2} models: Whereas the emergence of scaling in the former shifts to zero temperature if the bond dimension is increased, it occurs at a finite bond-dimension independent crossover temperature in the latter.

Scaling Hypothesis for Projected Entangled-Pair States

We introduce a new paradigm for scaling simulations with projected entangled-pair states (PEPS) for critical strongly correlated systems, allowing for reliable extrapolations of PEPS data with relatively small bond dimensions D. The key ingredient consists of using the effective correlation length ξ for inducing a collapse of data points, f(D,χ)=f(ξ(D,χ)), for arbitrary values of D and the environment bond dimension χ. As such we circumvent the need for extrapolations in χ and can use many distinct data points for a fixed value of D. Here, we need that the PEPSs have been optimized using a fixed-χ gradient method, which can be achieved using a novel tensor-network algorithm for finding fixed points of 2D transfer matrices, or by using the formalism of backwards differentiation. We test our hypothesis on the critical 3D dimer model, the 3D classical Ising model, and the 2D quantum Heisenberg model.

Tensor Networks Can Resolve Fermi Surfaces

We demonstrate that projected entangled-pair states are able to represent ground states of critical, fermionic systems exhibiting both 1d and 0d Fermi surfaces on a 2D lattice with an efficient scaling of the bond dimension. Extrapolating finite size results for the Gaussian restriction of fermionic projected entangled-pair states to the thermodynamic limit, the energy precision as a function of the bond dimension is found to improve as a power law, illustrating that an arbitrary precision can be obtained by increasing the bond dimension in a controlled manner. In this process, boundary conditions and system sizes have to be chosen carefully so that nonanalyticities of the Ansatz, rooted in its nontrivial topology, are avoided.