Detecting Emergent Continuous Symmetries at Quantum Criticality.

New or enlarged symmetries can emerge at the low-energy spectrum of a Hamiltonian that does not possess the symmetries, if the symmetry breaking terms in the Hamiltonian are irrelevant under the renormalization group flow. We propose a tensor network based algorithm to numerically extract lattice operator approximation of the emergent conserved currents from the ground state of any quantum spin chains, without the necessity to have prior knowledge about its low-energy effective field theory. Our results for the spin-1/2 J-Q Heisenberg chain and a one-dimensional version of the deconfined quantum critical points demonstrate the power of our method to obtain the emergent lattice Kac-Moody generators. It can also be viewed as a way to find the local integrals of motion of an integrable model and the local parent Hamiltonian of a critical gapless ground state.

Contrasting pseudocriticality in the classical two-dimensional Heisenberg and RP

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.

Scaling Hypothesis for Matrix Product States.

We study critical spin systems and field theories using matrix product states, and formulate a scaling hypothesis in terms of operators, eigenvalues of the transfer matrix, and lattice spacing in the case of field theories. The critical point, exponents, and central charge are determined by optimizing them to obtain a data collapse. We benchmark this method by studying critical Ising and Potts models, where we also obtain a scaling Ansatz for the correlation length and entanglement entropy. The formulation of those scaling functions turns out to be crucial for studying critical quantum field theories on the lattice. For the case of λϕ^{4} with mass parameter µ^{2} and lattice spacing a, we demonstrate a double data collapse for the correlation length δξ(µ,λ,D)=ξ[over ˜]((α-α_{c})(δ/a)^{-1/ν}) with D the bond dimension, δ the gap between eigenvalues of the transfer matrix, and α_{c}=µ_{R}^{2}/λ the parameter which fixes the critical quantum field theory.

Approaching the Kosterlitz-Thouless transition for the classical XY model with tensor networks.

We apply variational tensor-network methods for simulating the Kosterlitz-Thouless phase transition in the classical two-dimensional XY model. In particular, using uniform matrix product states (MPS) with non-Abelian O(2) symmetry, we compute the universal drop in the spin stiffness at the critical point. In the critical low-temperature regime, we focus on the MPS entanglement spectrum to characterize the Luttinger-liquid phase. In the high-temperature phase, we confirm the exponential divergence of the correlation length and estimate the critical temperature with high precision. Our MPS approach can be used to study generic two-dimensional phase transitions with continuous symmetries.