Incremental SWAP Operator for Entanglement Entropy: Application for Exponential Observables in Quantum Monte Carlo Simulation

Abstract

We propose a new method to efficiently compute the entanglement entropy (EE) of quantum many-body systems. Our approach, called the incremental SWAP operator method , combines the simplicity of the SWAP operator used in projector quantum Monte Carlo simulations with recent advances in precisely computing exponential observables using incremental algorithms. We apply this technique to obtain accurate EE data at reduced computational cost for 1d and 2d antiferromagnetic Heisenberg models with different bipartition schemes. Using the computed EE data, we extract the area law coefficient, universal logarithmic corrections from Goldstone modes, and the geometric constant, finding quantitative agreement with analytical predictions. Moreover, for the first time in an unbiased numerical simulation of 2d antiferromagnetic Heisenberg model, we successfully obtain reliable universal logarithmic corrections from sharp corners that match expected theoretical values. The consistency between our numerical results and theoretical calculations demonstrates the power of our approach for accessing challenging universal entanglement properties. The extensions of our method to other quantum spin/boson models and the interacting fermion models, are outlined.

References in corpus (9)

• Universal terms for the entanglement entropy in 2+1 dimensions
• Universality of corner entanglement in conformal field theories
• Entanglement Entropy of the Two-Dimensional Heisenberg Antiferromagnet
• Finite Size Scaling of Mutual Information: A Scalable Simulation
• Entanglement entropy scaling in the bilayer Heisenberg spin system
• Improving entanglement and thermodynamic Rényi entropy measurements in quantum Monte Carlo
• Entanglement entropy of the $\nu =1/2$ composite fermion non-Fermi liquid state
• A Wang-Landau method for calculating Renyi entropies in finite-temperature quantum Monte Carlo simulations
• Dirac fermions with plaquette interactions. III. SU(N) phase diagram with Gross-Neveu criticality and first-order phase transition