Search

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.

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.

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.

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.

We use the formalism of strange correlators to construct a critical classical lattice model in two dimensions with the Haagerup fusion category H_{3} as input data. We present compelling numerical evidence in the form of finite entanglement scaling to support a Haagerup conformal field theory (CFT) with central charge c=2. Generalized twisted CFT spectra are numerically obtained through exact diagonalization of the transfer matrix, and the conformal towers are separated in the spectra through their identification with the topological sectors. It is further argued that our model can be obtained through an orbifold procedure from a larger lattice model with input Z(H_{3}), which is the simplest modular tensor category that does not admit an algebraic construction. This provides a counterexample for the conjecture that all rational CFT can be constructed from standard methods.

Just as matrix product states represent ground states of one-dimensional quantum spin systems faithfully, continuous matrix product states (cMPS) provide faithful representations of the vacuum of interacting field theories in one spatial dimension. Unlike the quantum spin case, however, for which the density matrix renormalization group and related matrix product state algorithms provide robust algorithms for optimizing the variational states, the optimization of cMPS for systems with inhomogeneous external potentials has been problematic. We resolve this problem by constructing a piecewise linear parameterization of the underlying matrix-valued functions, which enables the calculation of the exact reduced density matrices everywhere in the system by high-order Taylor expansions. This turns the variational cMPS problem into a variational algorithm from which both the energy and its backwards derivative can be calculated exactly and at a cost that scales as the cube of the bond dimension. We illustrate this by finding ground states of interacting bosons in external potentials and by calculating boundary or Casimir energy corrections of continuous many-body systems with open boundary conditions.

We provide a generalization of the matrix product operator formalism for string-net projected entangled pair states (PEPS) to include nonunitary solutions of the pentagon equation. These states provide the explicit lattice realization of the Galois conjugated counterparts of (2+1)-dimensional topological quantum field theories, based on tensor fusion categories. Although the parent Hamiltonians of these renormalization group fixed point states are gapless, these states can still be the topological ground states of a gapped non-Hermitian Hamiltonian. We show by example that the topological sectors of the Yang-Lee theory (the nonunitary counterpart of the Fibonacci fusion category) can be constructed, even in the absence of closure under Hermitian conjugation of the basis elements of the Ocneanu tube algebra. The topological sector construction is demonstrated by applying the concept of strange correlators to the Yang-Lee model, giving rise to a nonunitary version of the classical hard hexagon model in the Yang-Lee universality class and obtaining all generalized twisted boundary conditions on a finite cylinder of the Yang-Lee edge singularity. Finally, we construct the PEPS transfer matrix and show that taking the Hermitian conjugate changes the topological phase for these nonunitary string-net models.

Although artificial neural networks have recently been proven to provide a promising new framework for constructing quantum many-body wave functions, the parametrization of a quantum wave function with non-abelian symmetries in terms of a Boltzmann machine inherently leads to biased results due to the basis dependence. We demonstrate that this problem can be overcome by sampling in the basis of irreducible representations instead of spins, for which the corresponding ansatz respects the non-abelian symmetries of the system. We apply our methodology to find the ground states of the one-dimensional antiferromagnetic Heisenberg (AFH) model with spin-1/2 and spin-1 degrees of freedom, and obtain a substantially higher accuracy than when using the s_{z} basis as an input to the neural network. The proposed ansatz can target excited states, which is illustrated by calculating the energy gap of the AFH model. We also generalize the framework to the case of anyonic spin chains.

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.

We extend the concept of strange correlators, defined for symmetry-protected phases in You et al. [Phys. Rev. Lett. 112, 247202 (2014)PRLTAO0031-900710.1103/PhysRevLett.112.247202], to topological phases of matter by taking the inner product between string-net ground states and product states. The resulting two-dimensional partition functions are shown to be either critical or symmetry broken, since the corresponding transfer matrices inherit all matrix product operator symmetries of the string-net states. For the case of critical systems, these nonlocal matrix product operator symmetries are the lattice remnants of topological conformal defects in the field theory description. Following Aasen et al. [J. Phys. A 49, 354001 (2016)JPAMB51751-811310.1088/1751-8113/49/35/354001], we argue that the different conformal boundary conditions can be obtained by applying the strange correlator concept to the different topological sectors of the string net obtained from Ocneanu's tube algebra. This is demonstrated on the lattice by calculating the conformal field theory spectra in the different topological sectors for the Fibonacci (hard-hexagon) and Ising string net. Additionally, we provide a complementary perspective on symmetry-preserving real-space renormalization by showing how known tensor network renormalization methods can be understood as the approximate truncation of an exactly coarse-grained strange correlator.

We study quasiparticle excitations for quantum spin chains with long-range interactions using variational matrix product state techniques. It is confirmed that the local quasiparticle ansatz is able to capture those excitations very accurately, even when the correlation length becomes very large and in the case of topological nontrivial excitation such as spinons. It is demonstrated that the breaking of the Lieb-Robinson bound follows from the appearance of cusps in the dispersion relation, and evidence is given for a crossover between different quasiparticles as the long-range interactions are tuned.

Edge theories of symmetry-protected topological phases are well known to possess global symmetry anomalies. In this Letter we focus on two-dimensional bosonic phases protected by an on-site symmetry and analyze the corresponding edge anomalies in more detail. Physical interpretations of the anomaly in terms of an obstruction to orbifolding and constructing symmetry-preserving boundaries are connected to the cohomology classification of symmetry-protected phases in two dimensions. Using the tensor network and matrix product state formalism we numerically illustrate our arguments and discuss computational detection schemes to identify symmetry-protected order in a ground state wave function.

We demonstrate that perturbative expansions for quantum many-body systems can be rephrased in terms of tensor networks, thereby providing a natural framework for interpolating perturbative expansions across a quantum phase transition. This approach leads to classes of tensor-network states parametrized by few parameters with a clear physical meaning, while still providing excellent variational energies. We also demonstrate how to construct perturbative expansions of the entanglement Hamiltonian, whose eigenvalues form the entanglement spectrum, and how the tensor-network approach gives rise to order parameters for topological phase transitions.

We study the entanglement structure of lattice gauge theories from the local operational point of view, and, similar to Soni and Trivedi [J. High Energy Phys. 1 (2016) 1], we show that the usual entanglement entropy for a spatial bipartition can be written as the sum of an undistillable gauge part and of another part corresponding to the local operations and classical communication distillable entanglement, which is obtained by depolarizing the local superselection sectors. We demonstrate that the distillable entanglement is zero for pure Abelian gauge theories at zero gauge coupling, while it is in general nonzero for the non-Abelian case. We also consider gauge theories with matter, and show in a perturbative approach how area laws-including a topological correction-emerge for the distillable entanglement. Finally, we also discuss the entanglement entropy of gauge fixed states and show that it has no relation to the physical distillable entropy.

The matrix product state formalism is used to simulate Hamiltonian lattice gauge theories. To this end, we define matrix product state manifolds which are manifestly gauge invariant. As an application, we study (1+1)-dimensional one flavor quantum electrodynamics, also known as the massive Schwinger model, and are able to determine very accurately the ground-state properties and elementary one-particle excitations in the continuum limit. In particular, a novel particle excitation in the form of a heavy vector boson is uncovered, compatible with the strong coupling expansion in the continuum. We also study full quantum nonequilibrium dynamics by simulating the real-time evolution of the system induced by a quench in the form of a uniform background electric field.

We use the matrix product state formalism to construct stationary scattering states of elementary excitations in generic one-dimensional quantum lattice systems. Our method is applied to the spin-1 Heisenberg antiferromagnet, for which we calculate the full magnon-magnon S matrix for arbitrary momenta and spin, the two-particle contribution to the spectral function, and higher order corrections to the magnetization curve. As our method provides an accurate microscopic representation of the interaction between elementary excitations, we envisage the description of low-energy dynamics of one-dimensional spin chains in terms of these particlelike excitations.

We prove an upper bound on the maximal rate at which a Hamiltonian interaction can generate entanglement in a bipartite system. The scaling of this bound as a function of the subsystem dimension on which the Hamiltonian acts nontrivially is optimal and is exponentially improved over previously known bounds. As an application, we show that a gapped quantum many-body spin system on an arbitrary lattice satisfies an area law for the entanglement entropy if and only if any other state with which it is adiabatically connected (i.e., any state in the same phase) also satisfies an area law.

For quantum lattice systems with local interactions, the Lieb-Robinson bound serves as an alternative for the strict causality of relativistic systems and allows the proof of many interesting results, in particular, when the energy spectrum exhibits an energy gap. In this Letter, we show that for translation invariant systems, simultaneous eigenstates of energy and momentum with an eigenvalue that is separated from the rest of the spectrum in that momentum sector can be arbitrarily well approximated by building a momentum superposition of a local operator acting on the ground state. The error satisfies an exponential bound in the size of the support of the local operator, with a rate determined by the gap below and above the targeted eigenvalue. We show this explicitly for the Affleck-Kennedy-Lieb-Tasaki model and discuss generalizations and applications of our result.

We construct a class of projected entangled pair states which is exactly the resonating valence bond wave functions endowed with both short range and long range valence bonds. With an energetically preferred resonating valence bond pattern, the wave function is simplified to live in a one-parameter variational space. We tune this variational parameter to minimize the energy for the frustrated spin-1/2 J(1)-J(2) antiferromagnetic Heisenberg model on the square lattice. Taking a cylindrical geometry, we are able to construct four topological sectors with an even or odd number of fluxes penetrating the cylinder and an even or odd number of spinons on the boundary. The energy splitting in different topological sectors is exponentially small with the cylinder perimeter. We find a power law decay of the dimer correlation function on a torus, and a lnL correction to the entanglement entropy, indicating a gapless spin-liquid phase at the optimum parameter.

We introduce a variational method for calculating dispersion relations of translation invariant (1+1)-dimensional quantum field theories. The method is based on continuous matrix product states and can be implemented efficiently. We study the critical Lieb-Liniger model as a benchmark and excellent agreement with the exact solution is found. Additionally, we observe solitonic signatures of Lieb's type II excitation. In addition, a nonintegrable model is introduced where a U(1)-symmetry breaking term is added to the Lieb-Liniger Hamiltonian. For this model we find evidence of a nontrivial bound-state excitation in the dispersion relation.

We show how to construct renormalization group (RG) flows of quantum field theories in real space, as opposed to the usual Wilsonian approach in momentum space. This is achieved by generalizing the multiscale entanglement renormalization ansatz to continuum theories. The variational class of wave functions arising from this RG flow are translation invariant and exhibits an entropy-area law. We illustrate the construction for a free nonrelativistic boson model, and argue that the full power of the construction should emerge in the case of interacting theories.

We present a quantum algorithm to prepare injective projected entangled pair states (PEPS) on a quantum computer, a class of open tensor networks representing quantum states. The run time of our algorithm scales polynomially with the inverse of the minimum condition number of the PEPS projectors and, essentially, with the inverse of the spectral gap of the PEPS's parent Hamiltonian.

The numerical renormalization group (NRG) is rephrased as a variational method with the cost function given by the sum of all the energies of the effective low-energy Hamiltonian. This allows us to systematically improve the spectrum obtained by NRG through sweeping. The ensuing algorithm has a lot of similarities to the density matrix renormalization group (DMRG) when targeting many states, and this synergy of NRG and DMRG combines the best of both worlds and extends their applicability. We illustrate this approach with simulations of a quantum spin chain and a single impurity Anderson model where the accuracy of the effective eigenstates is greatly enhanced as compared to the NRG, especially in the transition to the continuum limit.

We develop a new algorithm based on the time-dependent variational principle applied to matrix product states to efficiently simulate the real- and imaginary-time dynamics for infinite one-dimensional quantum lattices. This procedure (i) is argued to be optimal, (ii) does not rely on the Trotter decomposition and thus has no Trotter error, (iii) preserves all symmetries and conservation laws, and (iv) has low computational complexity. The algorithm is illustrated by using both an imaginary-time and a real-time example.

We consider the manifold of all quantum many-body states that can be generated by arbitrary time-dependent local Hamiltonians in a time that scales polynomially in the system size, and show that it occupies an exponentially small volume in Hilbert space. This implies that the overwhelming majority of states in Hilbert space are not physical as they can only be produced after an exponentially long time. We establish this fact by making use of a time-dependent generalization of the Suzuki-Trotter expansion, followed by a well-known counting argument. This also demonstrates that a computational model based on arbitrarily rapidly changing Hamiltonians is no more powerful than the standard quantum circuit model.

We show how continuous matrix product states of quantum fields can be described in terms of the dissipative nonequilibrium dynamics of a lower-dimensional auxiliary boundary field by demonstrating that the spatial correlation functions of the bulk field correspond to the temporal statistics of the boundary field. This equivalence (1) illustrates an intimate connection between the theory of continuous quantum measurement and quantum field theory, (2) gives an explicit construction of the boundary field allowing the extension of real-space renormalization group methods to arbitrary dimensional quantum field theories without the introduction of a lattice parameter, and (3) yields a novel interpretation of recent cavity QED experiments in terms of quantum field theory, and hence paves the way toward observing genuine quantum phase transitions in such zero-dimensional driven quantum systems.