Local Scale Transformations on the Lattice with Tensor Network Renormalization.

Consider the partition function of a classical system in two spatial dimensions, or the Euclidean path integral of a quantum system in two space-time dimensions, both on a lattice. We show that the tensor network renormalization algorithm [G. Evenbly and G. Vidal Phys. Rev. Lett. 115, 180405 (2015)] can be used to implement local scale transformations on these objects, namely, a lattice version of conformal maps. Specifically, we explain how to implement the lattice equivalent of the logarithmic conformal map that transforms the Euclidean plane into a cylinder. As an application, and with the 2D critical Ising model as a concrete example, we use this map to build a lattice version of the scaling operators of the underlying conformal field theory, from which one can extract their scaling dimensions and operator product expansion coefficients.

Tensor Network Renormalization.

We introduce a coarse-graining transformation for tensor networks that can be applied to study both the partition function of a classical statistical system and the Euclidean path integral of a quantum many-body system. The scheme is based upon the insertion of optimized unitary and isometric tensors (disentanglers and isometries) into the tensor network and has, as its key feature, the ability to remove short-range entanglement or correlations at each coarse-graining step. Removal of short-range entanglement results in scale invariance being explicitly recovered at criticality. In this way we obtain a proper renormalization group flow (in the space of tensors), one that in particular (i) is computationally sustainable, even for critical systems, and (ii) has the correct structure of fixed points, both at criticality and away from it. We demonstrate the proposed approach in the context of the 2D classical Ising model.

Tensor Network Renormalization Yields the Multiscale Entanglement Renormalization Ansatz.

We show how to build a multiscale entanglement renormalization ansatz (MERA) representation of the ground state of a many-body Hamiltonian H by applying the recently proposed tensor network renormalization [G. Evenbly and G. Vidal, Phys. Rev. Lett. 115, 180405 (2015)] to the Euclidean time evolution operator e(-ßH) for infinite ß. This approach bypasses the costly energy minimization of previous MERA algorithms and, when applied to finite inverse temperature ß, produces a MERA representation of a thermal Gibbs state. Our construction endows tensor network renormalization with a renormalization group flow in the space of wave functions and Hamiltonians (and not merely in the more abstract space of tensors) and extends the MERA formalism to classical statistical systems.

Class of highly entangled many-body states that can be efficiently simulated.

We describe a quantum circuit that produces a highly entangled state of N qubits from which one can efficiently compute expectation values of local observables. This construction yields a variational ansatz for quantum many-body states that can be regarded as a generalization of the multiscale entanglement renormalization ansatz (MERA), which we refer to as the branching MERA. In a lattice system in D dimensions, the scaling of entanglement of a region of size L(D) in the branching MERA is not subject to restrictions such as a boundary law L(D-1), but can be proportional to the size of the region, as we demonstrate numerically.

Real-space decoupling transformation for quantum many-body systems.

We propose a real-space renormalization group method to explicitly decouple into independent components a many-body system that, as in the phenomenon of spin-charge separation, exhibits separation of degrees of freedom at low energies. Our approach produces a branching holographic description of such systems that opens the path to the efficient simulation of the most entangled phases of quantum matter, such as those whose ground state violates a boundary law for entanglement entropy. As in the coarse-graining transformation of Vidal [Phys. Rev. Lett. 99, 220405 (2007).

Frustrated antiferromagnets with entanglement renormalization: ground state of the spin-1/2 Heisenberg model on a kagome lattice.

Entanglement renormalization techniques are applied to numerically investigate the ground state of the spin-1/2 Heisenberg model on a kagome lattice. Lattices of N={36,144, infinity} sites with periodic boundary conditions are considered. For the infinite lattice, the best approximation to the ground state is found to be a valence bond crystal with a 36-site unit cell, compatible with a previous proposal. Its energy per site, E=-0.432 21, is an exact upper bound and is lower than the energy of any previous (gapped or algebraic) spin liquid candidate for the ground state.

Entanglement renormalization in two spatial dimensions.

We propose and test a scheme for entanglement renormalization capable of addressing large two-dimensional quantum lattice systems. In a translationally invariant system, the cost of simulations grows only as the logarithm of the lattice size; at a quantum critical point, the simulation cost becomes independent of the lattice size and infinite systems can be analyzed. We demonstrate the performance of the scheme by investigating the low energy properties of the 2D quantum Ising model on a square lattice of linear size L={6, 9, 18, 54, infinity} with periodic boundary conditions. We compute the ground state and evaluate local observables and two-point correlators. We also produce accurate estimates of the critical magnetic field and critical exponent beta. A calculation of the energy gap shows that it scales as 1/L at the critical point.