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.
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.
Abstract We study the entanglement entropy in lattice field theory using a simulation algorithm based on Jarzynski’s theorem. We focus on the entropic c-function for the Ising model in two and in three dimensions: after validating our algorithm against known analytical results from conformal field theory in two dimensions, we present novel results for the three-dimensional case. We show that our algorithm, which is highly parallelized on graphics processing units, allows one to precisely determine the subleading corrections to the area law, which have been investigated in many recent works. Possible generalizations of this study to other strongly coupled theories are discussed.
A bstract The study of entanglement in gauge theories is expected to provide insights into many fundamental phenomena, including confinement. However, calculations of quantities related to entanglement in gauge theories are limited by ambiguities that stem from the non-factorizability of the Hilbert space. In this work we study lattice gauge theories that admit a dual description in terms of spin models, for which the replica trick and Rényi entropies are well defined. In the first part of this work, we explicitly perform the duality transformation in a replica geometry, deriving the structure of a replica space for a gauge theory. Then, in the second part, we calculate, by means of Monte Carlo simulations, the entropic c-function of the ℤ 2 gauge theory in three spacetime dimensions, exploiting its dual description in terms of the three-dimensional Ising model.
One of the most important predictions of quantum mechanics is the existence of quantum entanglement. Since the inception of quantum theory, the study of entanglement has become ubiquitous in several branches of physics, from high energy to condensed matter physics and quantum information. To this avail, different measures of entanglement in quantum systems have been introduced. The most successful ones are the von Neumann and the Rényi entanglement entropies, which have been heavily studied both in lattice models and in quantum field theories (QFTs). A different, more complete characterization of entanglement is provided by the entanglement (or modular) Hamiltonian, defined as the logarithm of the reduced density matrix. One of the key property of this operator is that, under generic conditions, in 1+1-dimensional conformal field theories (CFTs) presents a local and few-body structure. Beyond just being of theoretical interest, this property has been recently used to efficiently construct the ground state of spin chains in quantum simulators. A recently introduced generalisation of the entanglement Hamiltonian is the negativity Hamiltonian, a non-Hermitian operator which is more useful for mixed states and tripartite entanglement. In this thesis we study these quantifications in many-body quantum systems and in QFTs. In the first part of the thesis, we focus on fermionic systems at equilibrium, and we obtain analytical expressions for the entanglement and the negativity Hamiltonians in the massless Dirac fermion QFT, both in presence of boundaries and at finite temperature. We compare these predictions in the continuum QFT with numerical calculations on the lattice, showing that a proper comparison requires a non-trivial continuum limit of the lattice expression, which we detail. We also investigate a lattice realisation of a non-unitary CFT, the bc-ghost theory, showing that the entanglement Hamiltonian presents an additional term with respect to the unitary case. In the second part, we study the entanglement entropies in an integrable QFT, the massless renormalisation group flow from the tricritical to the critical Ising CFT, demonstrating that the entropy is monotonically decreasing under the renormalisation group flow, as expected from the entropic c-theorem. In the third part we consider out-of-equilibrium problems. We study the melting of a domain wall configuration, computing analytically the entanglement Hamiltonian and showing that the insertion of a defect dramatically modifies the entanglement properties, changing the growth of the entropies from logarithmic to linear in time. Finally, we study generic quantum quenches in free fermionic theories, using the celebrated quasiparticle picture to obtain a prediction for the evolution of the entanglement Hamiltonian during the quench.
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.
Entanglement is a fundamental feature of quantum mechanics and plays an important role in characterizing many-body systems and quantum phase transitions. A quantitative measure of quantum correlations between subsystems is entanglement entropy. In this thesis, we study the entanglement entropy of an interacting quantum field theory, the O(N) model. In local quantum field theories, the leading-order contribution in the scaling of entanglement entropy obeys an area law. This dominant term is governed by short-distance physics near the entangling surface and depends on the ultraviolet regularization, making it non-universal. Beyond this lead- ing scaling behavior, entanglement entropy can contain subleading universal contributions, which encode intrinsic information about the underlying quantum field theory. In this thesis, we calculate one such universal term: the correlation length coefficient r, which appears as a correction to the entanglement entropy evaluated slightly away from the quantum critical point. The analysis is carried out in D = 2 + 1 spacetime dimensions for an infinite cylindrical geometry. We study the model in the vicinity of criticality using the large-N expansion, which enables controlled analytical calculations. Our main result is that the coefficient r scales linearly with the number of field components N . This behavior is consistent with the interpretation of r as a universal quantity associated with the degrees of freedom of the theory. The calculation is based on the replica method and renormalization group arguments, which are presented in detail. Finally, we discuss the physical meaning of the result and possible methods for verifying it, and we review related results for the O(N) model in different geometries.
Computing the subleading logarithmic term in the entanglement entropy (EE) of (2+1)d quantum many-body systems remains a significant challenge, despite its central role in revealing universal information about quantum states and quantum critical points (QCPs). Building on recent algorithmic advances that enable the stable calculation of EE as an exponential observable~\cite{zhouIncremental2024,zhangIntegral2024,liaoExtracting2024}, we develop a {\it bubble basis} projector quantum Monte Carlo (QMC) algorithm to precisely and efficiently compute the universal corner of EE at QCPs in a (2+1)d square-lattice transverse-field Ising model augmented with a four-body interaction. Turning on this interaction allows us to trace an Ising critical line, reaching the tricritical point, and then a line of first-order phase transition. In (2+1)d, the tricritical point is described by the Gaussian theory, where a theoretical calculation of the corner logarithmic term in the 2nd Rényi entropy term is available~\cite{UniversalCasini2007}. Our QMC results are in quantitative agreement with this theoretical value, providing a highly nontrivial benchmark of the algorithm. Furthermore, we also study the Rényi EE at the Ising critical line and on the first-order transition line, obtaining results consistent with theoretical expectations. These findings establish the long-sought connection between the universal values of an exactly solvable limit and those of a strongly correlated regime at (2+1)d.