The fate of quantum entanglement at finite-temperature phase transitions remains an open question, particularly for continuous symmetry breaking where zero-temperature Goldstone modes generate long-range correlations. Using large-scale quantum Monte Carlo simulations, we investigate the third Rényi negativity across the O(3) transition in the three-dimensional Heisenberg antiferromagnet. The first such study for a thermal critical point with continuous symmetry. We uncover two fundamental results. First, the negativity exhibits a pure area law at the critical point, with the subleading constant term vanishing within statistical uncertainty. This demonstrates that thermal fluctuations completely destroy the long-range entanglement present at zero temperature. The divergent classical correlation length leaves no imprint on quantum entanglement itself. Second, despite this absence of singular behavior in the negativity, its temperature derivative follows the exact scaling of the specific heat, yielding critical exponents -α/ν=0.190(1) and 1/ν=1.350(5) in precise agreement with the O(3) universality class. Our work establishes that while quantum entanglement is blind to thermal criticality, its thermodynamic derivatives encode the full universal scaling, revealing an unexpected connection between entanglement and classical phase transitions. Furthermore, this also provides further evidence that Rényi negativity can still effectively shield classical correlations in systems with continuous symmetry.
Quantum entanglement is fragile to thermal fluctuations, which raises the question whether finite temperature phase transitions support long-range entanglement similar to their zero temperature counterparts. Here we use quantum Monte Carlo simulations to study the third Renyi negativity, a generalization of entanglement negativity, as a proxy of mixed-state entanglement in the 2D transverse field Ising model across its finite temperature phase transition. We find that the area-law coefficient of the Renyi negativity is singular across the transition, while its subleading constant is zero within the statistical error. This indicates that the entanglement is short-range at the critical point despite a divergent correlation length. Renyi negativity in several exactly solvable models also shows qualitative similarities to that in the 2D transverse field Ising model.
Abstract We investigate several entanglement-related quantities at finite-temperature criticality in the three-dimensional quantum spherical model, both as a function of temperature T and of the quantum parameter g , which measures the strength of quantum fluctuations. While the von Neumann and the Rényi entropies exhibit a volume-law for any g and T , the mutual information obeys an area law. The prefactors of the volume-law and of the area-law are regular across the transition, reflecting that universal singular terms vanish at the transition. This implies that the mutual information is dominated by nonuniversal contributions. This hampers its use as a witness of criticality, at least in the spherical model. We also study the logarithmic negativity. For any value of g , T , the negativity exhibits an area-law. The negativity vanishes deep in the paramagnetic phase, it is larger at small temperature, and it decreases upon increasing the temperature. For any g , it exhibits the so-called sudden death, i.e. it is exactly zero for large enough T . The vanishing of the negativity defines a ‘death line’, which we characterise analytically at small g . Importantly, for any finite T the area-law prefactor is regular across the transition, whereas it develops a cusp-like singularity in the limit . Finally, we consider the single-particle entanglement and negativity spectra. The vast majority of the levels are regular across the transition. Only the larger ones exhibit singularities. These are related to the presence of a zero mode, which reflects the symmetry breaking. This implies the presence of sub-leading singular terms in the entanglement entropies. Interestingly, since the larger levels do not contribute to the negativity, sub-leading singular corrections are expected to be suppressed for the negativity.
We investigate several entanglement-related quantities at finite-temperature criticality in the three-dimensional quantum spherical model, both as a function of temperature
T and of the quantum parameter
g, which measures the strength of quantum fluctuations. While the von Neumann and the R\\'enyi entropies exhibit the volume-law for any
g and
T, the mutual information obeys an area law. The prefactors of the volume-law and of the area-law are regular across the transition, reflecting that universal singular terms vanish at the transition. This implies that the mutual information is dominated by nonuniversal contributions. This hampers its use as a witness of criticality, at least in the spherical model. We also study the logarithmic negativity. For any value of
g,T, the negativity exhibits an area-law. The negativity vanishes deep in the paramagnetic phase, it is larger at small temperature, and it decreases upon increasing the temperature. For any
g, it exhibits the so-called sudden death, i.e., it is exactly zero for large enough
T. The vanishing of the negativity defines a "death line", which we characterise analytically at small
g. Importantly, for any finite
T the area-law prefactor is regular across the transition, whereas it develops a cusp-like singularity in the limit
Tto0. Finally, we consider the single-particle entanglement and negativity spectra. The vast majority of the levels are regular across the transition. Only the larger ones exhibit singularities. These are related to the presence of a zero mode, which reflects the symmetry breaking. This implies the presence of sub-leading singular terms in the entanglement entropies. Interestingly, since the larger levels do not contribute to the negativity, sub-leading singular corrections are expected to be suppressed for the negativity.
Establishing the nature of a quantum phase transition in finite-size simulations -- whether continuous, first-order, or weak first-order -- is a fundamental challenge in quantum many-body computation. Especially, the weak first-order phase transition is affected by a super large correlation length and always displays as a continuous critical point in simulated finite-sizes. The core difficulty lies in the fact that there is no effective finite-size theory to distinguish these phase transitions in the realistic simulations limited by the computational resource. In this work, we have fixed this problem by introducing a unified finite-size framework based on the Renyi thermal entropy (RTE) and its derivative (DRTE) to detect and characterize quantum phase transitions. We derive complete scaling theories for the RTE and DRTE at second-order, first-order, and weak first-order transitions, showing that the DRTE naturally isolates the singular part of the free energy and strengthens the characteristics of various phase transitions in finite sizes. Using quantum Monte Carlo simulations, we demonstrate accurate data collapse and extraction of critical exponents at (2+1)-dimensional O(
N) critical points. More importantly, the DRTE provides a smoking-gun signature of weak first-order transitions through a clear double-peak structure and a crossing at zero, which we unambiguously observe in debated deconfined quantum criticality candidates such as the
J--
Q models. Our approach offers a general, unbiased, and numerically efficient tool for probing the universal properties of quantum phase transitions, resolving long-standing ambiguities between continuous and weak first-order scenarios.
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ényi EE both along the critical line and in the Néel 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éel 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éel 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.
At a quantum critical point, the universal scaling behavior of Rényi entanglement entropy is controlled by the universality class of the codimension-two Rényi (or conical) defects in the infrared theory. In this work we perform a systematic study of critical correlations along Rényi defect lines in (2+1)
D quantum spin models realizing quantum phase transitions described by the O(3) Wilson-Fisher universality class, using large-scale quantum Monte Carlo simulations. We present numerical evidence that, for a fixed Rényi index
n, there exist multiple Rényi defect universality classes, with distinct critical exponents for the O(3) order parameter on the defect. These universality classes are realized by choosing microscopically different entanglement cuts in lattice models, which we classify as ordinary, special and extraordinary according to their relation to surface criticality. For the extraordinary entanglement cut, we further find evidence for a phase transition on the defect as a function of the Rényi index. Our results highlight the key role of defect universality classes in determining the universal scaling of Rényi entropy, and provide a framework for understanding the previously observed dependence of Rényi entropy scaling on microscopic lattice details.