The Rényi entanglement entropy in quantum many-body systems can be viewed as the difference in free energy between partition functions with different trace topologies. We introduce an external field λ that controls the partition function topology, allowing us to define a notion of nonequilibrium work as λ is varied smoothly. Nonequilibrium fluctuation theorems of the work provide us with statistically exact estimates of the Rényi entanglement entropy. This framework also naturally leads to the idea of using quench functions with spatially smooth profiles, providing us a way to average over lattice scale features of the entanglement entropy while preserving long distance universal information. We use these ideas to extract universal information from quantum Monte Carlo simulations of SU(N) spin models in one and two dimensions. The vast gain in efficiency of this method allows us to access unprecedented system sizes up to 192×96 spins for the square lattice Heisenberg antiferromagnet.
The Rényi entanglement entropy in quantum many-body systems can be viewed as the difference in free energy between partition functions with different trace topologies. We introduce an external field λ that controls the partition function topology, allowing us to define a notion of nonequilibrium work as λ is varied smoothly. Nonequilibrium fluctuation theorems of the work provide us with statistically exact estimates of the Rényi entanglement entropy. This framework also naturally leads to the idea of using quench functions with spatially smooth profiles, providing us a way to average over lattice scale features of the entanglement entropy while preserving long distance universal information. We use these ideas to extract universal information from quantum Monte Carlo simulations of SU(N) spin models in one and two dimensions. The vast gain in efficiency of this method allows us to access unprecedented system sizes up to 192×96 spins for the square lattice Heisenberg antiferromagnet.
Abstract We develop a nonequilibrium increment method in quantum Monte Carlo simulations to obtain the Rényi entanglement entropy of various quantum many-body systems with high efficiency and precision. To demonstrate its power, we show the results on a few important yet difficult (2 + 1)dquantum lattice models, ranging from the Heisenberg quantum antiferromagnet with spontaneous symmetry breaking, the quantum critical point with O(3) conformal field theory (CFT) to the toric code Z2 Z2 topological ordered state and the Kagome Z2 Z2 quantum spin liquid model with frustration and multi-spin interactions. In all these cases, our method either reveals the precise CFT data from the logarithmic correction or extracts the quantum dimension in topological order, from the dominant area law in finite-size scaling, with very large system sizes, controlled errorbars, and minimal computational costs. Our method, therefore, establishes a controlled and practical computation paradigm to obtain the difficult yet important universal properties in highly entangled quantum matter.
Our understanding of physical phenomena hinges on finding universal core mechanisms that unite them. The concept of universality is deeply ingrained in the study of quantum many-body systems. At zero temperature, microscopically different systems with long-range order collapse into their universal state described by a handful of universal parameters. Establishing those parameters implies identification of the universal theory that is effectively describing the system and ultimately providing the desired understanding. We take up this task with the help of Renyi entanglement entropy. Defined with respect to a system bipartition, this measure quantifies information shared between the subsystems. Understanding what insights into universality are encoded within this information-theoretic quantity as well as developing numerical tools to efficiently estimate the Renyi entanglement entropy are the subjects of this thesis. \n \nOn the computational end of this far-reaching goal, we develop a novel theoretical framework for constructing improved Renyi entanglement entropy estimators in the context of d+1 quantum Monte Carlo methods. The discovery of a connection of this methodology to the well-established Kandel-Domany formalism provides a clear path towards generalization. Additionally, we embrace a data-driven approach towards learning the ground state wavefunction. We demonstrate how a restricted Boltzmann machine can be used to reconstruct the Renyi entanglement entropy from projective measurements of a quantum ground state. Furthermore, we extend this classical generative architecture to a quantum analogue that we call the quantum Boltzmann machine. \n \nOn the theoretical side of this endeavour, we study the Renyi entanglement entropy scaling terms for two quantum lattice models embedded in two dimensional space via extensive quantum Monte Carlo simulations. For the ground state of the XY model, we provide conclusive numerical evidence for a logarithmic contribution that uniquely characterizes the continuous symmetry of the emerging order parameter. Moreover, we confirm the form of the subleading universal geometric contribution arising due to the bosonic nature of low-energy degrees of freedom in this model. For the critical ground state of the transverse field Ising model, we develop a novel scaling procedure to extract a universal number κ revealed via a cylindrical entangling bipartition in the thin-slice limit. The combined product of our work sheds new light on the entanglement-based classification of universality and brings a suite of new powerful numerical tools to continue illuminating this theoretical program in the future.
Quantum magic, or nonstabilizerness, provides a crucial characterization of quantum systems, regarding the classical simulability with stabilizer states. In this work, we propose an alternative and efficient algorithm for computing stabilizer R\'enyi entropy, one of the measures for quantum magic, in spin systems with sign-problem free Hamiltonians. This algorithm is based on the quantum Monte Carlo simulation of the path integral of the work between two partition function ensembles and it applies to all spatial dimensions and temperatures. We demonstrate this algorithm on the one- and two-dimensional transverse field Ising model at both finite and zero temperatures and show the quantitative agreements with tensor-network based algorithms. We analyze the computational cost and provide analytical and numerical evidences for it to be polynomial in system size. This work also suggests a unifying framework for calculating various types of entropy quantities including entanglement R\'enyi entropy and entanglement R\'enyi negativity.
In recent years entanglement measures, such as the von Neumann and the Rényi entropies, provided a unique opportunity to access elusive features of quantum many-body systems. However, extracting entanglement properties analytically, experimentally, or in numerical simulations can be a formidable task. Here, by combining the replica trick and the Jarzynski equality we devise an alternative effective out-of-equilibrium protocol for measuring the equilibrium Rényi entropies. The key idea is to perform a quench in the geometry of the replicas. The Rényi entropies are obtained as the exponential average of the work performed during the quench. We illustrate an application of the method in classical Monte Carlo simulations, although it could be useful in different contexts, such as in quantum Monte Carlo, or experimentally in cold-atom systems. The method is most effective in the quasistatic regime, i.e., for a slow quench. As a benchmark, we compute the Rényi entropies in the Ising universality class in 1+1 dimensions. We find perfect agreement with the well-known conformal field theory predictions.
Entanglement entropy encodes important features of strongly interacting quantum many-body systems and gauge theories, but its analytical study is still limited to systems with high level of symmetry. This motivates the search for efficient techniques to investigate this quantity numerically, through Monte Carlo calculations on the lattice. In this contribution, we discuss the computation of the entropic c-function by means of an algorithm based on Jarzynski's equality, which is an exact theorem from non-equilibrium statistical mechanics. After presenting benchmark results for the Ising model in two dimensions, where our algorithm successfully reproduces the analytical predictions from conformal field theory, we discuss its generalization to the three-dimensional Ising model, for which we were able to extract universal terms beyond the area law. Finally we point out some future generalizations of this calculation.
The large-scale behaviour of entanglement entropy in finite-density states, in and out of equilibrium, can be understood using the physical picture of particle pairs. However, the full theoretical origin of this picture is not fully established yet. In this work, we clarify this picture by investigating entanglement entropy using its connection with the large-deviation theory for thermodynamic and hydrodynamic fluctuations. We apply the universal framework of Ballistic Fluctuation Theory (BFT), based the Euler hydrodynamics of the model, to correlation functions of branch-point twist fields, the starting point for computing Rényi entanglement entropies within the replica approach. Focusing on free fermionic systems in order to illustrate the ideas, we show that both the equilibrium behavior and the dynamics of Rényi entanglement entropies can be fully derived from the BFT. In particular, we emphasise that long-range correlations develop after quantum quenches, and accounting for these explain the structure of the entanglement growth. We further show that this growth is related to fluctuations of charge transport, generalising to quantum quenches the relation between charge fluctuations and entanglement observed earlier. The general ideas we introduce suggest that the large-scale behaviour of entanglement has its origin within hydrodynamic fluctuations.