In a
D=2+1 quantum critical system, the entanglement entropy across a boundary with a corner contains a subleading logarithmic scaling term with a universal coefficient. It has been conjectured that this coefficient is, to leading order, proportional to the number of field components
N in the associated
O(N) continuum ${\ensuremath{\phi}}^{4}$ field theory. Using density matrix renormalization group calculations combined with the powerful numerical linked cluster expansion technique, we confirm this scenario for the
O(2) Wilson-Fisher fixed point in a striking way, through direct calculation at the quantum critical points of two very different microscopic models. The value of this corner coefficient is, to within our numerical precision, twice the coefficient of the Ising fixed point. Our results add to the growing body of evidence that this universal term in the R\'enyi entanglement entropy reflects the number of low-energy degrees of freedom in a system, even for strongly interacting theories.
Entanglement scaling provides a powerful probe of universal properties at quantum critical points. In two dimensions, contributions originating from a corner-shaped bipartition exhibit a universal scaling, which is determined by the underlying conformal field theory. We develop a method to extract this corner entanglement entropy from projected entangled pair states directly in the thermodynamic limit. When applied to models at a quantum critical point, we show that the corner contribution exhibits scaling with the effective correlation length, in agreement with the hypothesis of finite-entanglement scaling. Our results for the corner coefficients are consistent with other methods, demonstrating the efficiency of our method for diagnosing strongly-correlated quantum critical points in two dimensions.
The subleading corner logarithmic corrections in entanglement entropy (EE) are crucial for revealing the universal characteristics of quantum critical points (QCPs), but they are challenging to detect. Motivated by recent developments in the stable computation of EE in (
2+1)-dimensional [(
2+1)D] quantum many-body systems, we have developed a method for directly measuring the corner contribution in EE with lower computational cost. The cornerstone of our approach is to measure the subtracted corner entanglement entropy (SCEE), defined as the difference between the EEs of subregions with the same boundary length for smooth and cornered boundaries during a sign-problem-free quantum Monte Carlo simulation. Our improved method inherently eliminates not only the area-law term of EE but also the subleading log corrections arising from Goldstone modes, leaving the universal corner contribution as the leading term of the SCEE with greatly improved data quality. Utilizing this advanced approach, we calculate the SCEE of the bilayer Heisenberg model on both square and honeycomb lattices across their (
2+1)D
O(3) QCPs with different opening angles on the entanglement boundary and obtain accurate values of the corresponding universal corner log coefficients. These findings will encourage further theoretical investigations to access controlled universal information for interacting conformal field theories at (
2+1) dimensions.