Quantum criticality and entanglement for two dimensional long-range Heisenberg bilayer

The study of quantum criticality and entanglement in systems with long-range (LR) interactions is still in its early stages, with many open questions remaining. In this work, we investigate critical exponents and scaling of entanglement entropies (EE) in the LR bilayer Heisenberg model using large-scale quantum Monte Carlo (QMC) simulations and the recently developed nonequilibrium increment algorithm for measuring EE. By applying modified (standard) finite-size scaling (FSS) above (below) the upper critical dimension and field theory analysis, we obtain precise critical exponents in three regimes: the LR Gaussian regime with a Gaussian fixed point, the short-range (SR) regime with Wilson-Fisher (WF) exponents, and a LR non-Gaussian regime where the critical exponents vary continuously from LR Gaussian to SR values. We compute the R\'enyi EE both along the critical line and in the N\'eel phase and observe that as the LR interaction is enhanced, the area-law contribution in EE gradually vanishes both at quantum critical points (QCPs) and in the N\'eel phase. The log-correction in EE arising from sharp corners at the QCPs also decays to zero as LR interaction grows, whereas the log-correction for N\'eel states, caused by the interplay of Goldstone modes and restoration of the symmetry in a finite system, is enhanced as LR interaction becomes stronger. We also discuss relevant experimental settings to detect these nontrivial properties in critical behavior and entanglement information for quantum many-body systems with LR interactions.