Continuous tensor networks are variational wavefunctions proposed in recent years to efficiently simulate quantum field theories (QFTs). Prominent examples include the continuous matrix product state (cMPS) and the continuous multi-scale entanglement renormalization ansatz (cMERA). While the cMPS can approximate ground states of a class of QFT Hamiltonians that are both local and interacting, cMERA is only well-understood for QFTs that are quasi-local and non-interacting. In this paper we propose the magic cMERA, a concrete realization of cMERA for a free boson QFT that simultaneously satisfies four remarkable properties: (i) it is the exact ground state of a strictly local Hamiltonian; (ii) in the massless case, its spectrum of scaling operators is exactly soluble in real space; (iii) it has the short-distance structure of a cMPS; (iv) it is generated by a quasi-local entangler that can be written as a continuous matrix product operator. None of these properties is fulfilled by previous cMERA proposals. Properties (iii)-(iv) establish a firm connection between cMERA and cMPS wavefunctionals, opening the path to applying powerful cMPS numerical techniques, valid for interacting QFTs, also to cMERA calculations.
We numerically study the one-dimensional long-range transverse field Ising model (TFIM) in the antiferromagnetic (AFM) regime at zero temperature using generalized Hartree-Fock (GHF) theory. The spin-spin interaction extends to all spins in the lattice and decays as $1/{r}^{\ensuremath{\alpha}}$, where
r denotes the distance between two spins and $\ensuremath{\alpha}$ is a tunable exponent. We map the spin operators to Majorana operators and approximate the ground state of the Hamiltonian with a fermionic Gaussian state (FGS). Using this approximation, we calculate the ground-state energy and the entanglement entropy, which allows us to map the phase diagram for different values of $\ensuremath{\alpha}$. In addition, we compute the scaling behavior of the entanglement entropy with the system size to determine the central charge at criticality for the case of $\ensuremath{\alpha}>1$. For $\ensuremath{\alpha}<1$ we find a logarithmic divergence of the entanglement entropy even far away from the critical point, a feature of systems with long-range interactions. We provide a detailed comparison of our results to outcomes of density matrix renormalization group (DMRG) and the linked cluster expansion (LCE) methods. In particular, we find excellent agreement of GHF with DMRG and LCE in the weak long-range regime $\ensuremath{\alpha}\ensuremath{\ge}1$, and qualitative agreement with DMRG in the strong-long range regime $\ensuremath{\alpha}\ensuremath{\le}1$. Our results highlight the power of the computationally efficient GHF method in simulating interacting quantum systems.
Tensor network methods are playing a central role in multiple disciplines of modern quantum physics, such as condensed matter physics, quantum information and quantum gravity. Based on renormalization group (RG) ideas, they provide a large class of variational ansatz for quantum many-body systems. There are various tensor network structures which are suited for different systems based on their entanglement structures, such as matrix product states (MPS), projected entangled pair states (PEPS), multi-scale entanglement renormalization ansatz (MERA). Numerous techniques involving coarse-graining tensor networks have been developed, such as tensor renormalization group (TRG), tensor network renormalization (TNR) and its variants. \n \nFirst, we study the generalization of MERA to continuous systems, or cMERA, which is expected to become a powerful variational ansatz for the ground state of strongly interacting quantum field theories. We investigate, in the simpler context of Gaussian cMERA for free theories, the extent to which the cMERA state
ketPsicM with finite UV cut-off
Lambda can capture the spacetime symmetries of the ground state
ketPsi. For a free boson conformal field theory (CFT) in 1+1 dimensions as a concrete example, we build a quasi-local unitary transformation
V that maps
ketPsi into
ketPsicM and show two main results. (i) Any spacetime symmetry of the ground state
ketPsi is also mapped by
V into a spacetime symmetry of the cMERA
ketPsicM. However, while in the CFT the stress-energy tensor
Tmunu(x) (in terms of which all the spacetime symmetry generators are expressed) is local, the corresponding cMERA stress-energy tensor
TmunucM(x)=VTmunu(x)Vdagger is quasi-local. (ii) From the cMERA, we can extract quasi-local scaling operators
OalphacM(x) characterized by the exact same scaling dimensions
Deltaalpha, conformal spins
salpha, operator product expansion coefficients
Calphabetagamma, and central charge
c as the original CFT. We argue that these results should also apply to interacting theories. \n \nSecond we extend TNR to field theories in the continuum. A short-distance length scale
1/Lambda is introduced in the continuum partition function by smearing the fields. The resulting object is still defined in the continuum but has no fluctuations at distances shorter than
1/Lambda. An infinitesimal coarse-graining step is then generated by the combined action of a \\textit{rescaling} operator
mathfrakL and a \\textit{disentangling} operator
mathfrakK that implements a quasilocal field redefinition. As demonstrated for a free boson in two dimensions, continuous TNR exactly preserves translation and rotation symmetries and can generate a proper RG flow. Moreover, from a critical fixed point of this RG flow one can then extract the conformal data of the underlying conformal field theory. \n \nLeaving behind tensor network methods, we explore the universal structure that emerges from the eigenvectors of the reduced density matrix (RDM) for 1+1 dimensional critical lattice systems. Using the Ising model as an example, we demonstrate that the large-weight eigenvectors are related to the low-lying states of the underlying CFT with suitable boundary conditions by a conformal map, and they transform accordingly under the Virasoro generators. \n \nFinally, using the quantum renormalization group (QRG), we derive the bulk geometry that emerges in the holographic dual of the fermionic U(
N) vector model at a nonzero charge density. The obstruction that prohibits the metallic state from being smoothly deformable to the direct product state under the renormalization group flow gives rise to a horizon at a finite radial coordinate in the bulk. The region outside the horizon is described by the Lifshitz geometry with a higher-spin hair determined by microscopic details of the boundary theory. On the other hand, the interior of the horizon is not described by any Riemannian manifold, as it exhibits an algebraic non-locality. The non-local structure inside the horizon carries the information on the shape of the filled Fermi sea.
At a quantum critical point, the low-energy physics of a quantum spin chain is described by conformal field theory (CFT). Given the Hamiltonian of a critical spin chain, we propose a variational method to build an approximate lattice representation ϕ_{α} of the corresponding primary CFT operators ϕ_{α}^{CFT}. We then show how to numerically compute the operator product expansion coefficients C_{αβγ}^{CFT} governing the fusion of two primary fields. In this way, we complete the implementation of Cardy's program, outlined in the 1980s, which aims to extract the universality class of a phase transition, as encoded in the so-called conformal data of the underlying CFT, starting from a microscopic description. Our approach, demonstrated here for the critical quantum Ising model, only requires a generic (i.e., in general, nonintegrable) critical lattice Hamiltonian as its input.
Ising model has been successful in describing ferromagnetism and its phase transition to paramagnet. At the critical point, the free energy density function and its derivatives diverge. Their behaviour near the critical point are described by power-laws with associated critical exponents. In many cases, the critical exponents can be determined from analytic solutions via conventional renormalization groups methods or from Monte Carlo simulations. However, for quantum many-body systems, very few are tractable to analytical solutions. The quantum many-body wavefunction belongs to large dimensional Hilbert space that increases exponentially with system size. If the Hamiltonian is gapped and only local interaction is considered, then the wavefunction can be efficiently truncated. Tensor Network formalism provides a scheme to truncate the less important degrees of freedom via Singular Value Decomposition (SVD) of the density matrix. In this study, we investigated the thermodynamic properties and phase transition of one-dimensional transverse-field quantum Ising model (1D-tQIM) under the finite-size effect and random coupling strength. Starting with Matrix Product States (MPS) as a wavefunction ansatz, the Density Matrix Renormalization Group algorithm is applied to the MPS. The variational algorithm, which iteratively performs SVDs and truncation at each bond, approximates the ground state MPS wavefunction. All quantum observables are calculated from the contraction of the resultant ground state MPS. Although theoretically, divergence at critical points only happen in an infinite system, one can obtain the critical exponents through simulation of finite-size 1D-tQIM. Using the analytic solution as a benchmark, we compared the finite-size effects of the system using finite-size scaling analysis and MPS methods. The critical exponents of 1D-tQIM are independently calculated and compared with the analytical results. Thermodynamic quantities such as magnetization, susceptibility and correlation function are calculated for system sizes of 20, 40, 60, 80, 100 and 120 spins. We determined the respective critical exponents: