Precise computation of universal corner entanglement entropy at 2+1 dimension: from Ising to Gaussian quantum...

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 (Zhouet al2024Phys. Rev. B109165106; Zhanget al2024Phys. Rev. B109205147; Liaoet al2024Phys. Rev. B110235111), we develop abubble basisprojector 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 (Casini and Huerta 2007Nucl. Phys. B764183). 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.