Third-Order Negativity in Quantum Entanglement

• Third-order negativity is a diagnostic measure derived from the third moment of a partially transposed density operator that quantifies mixed-state and tripartite entanglement.
• It utilizes constructions like the third Rényi negativity and the invariant I5 to analyze finite-temperature phase transitions and separability criteria in multipartite settings.
• The approach reveals universal scaling laws, boundary-local singularities, and multifractal behavior in both disordered systems and conformal field theory contexts.

Third-order negativity denotes a class of third-moment diagnostics built from a partial transpose and used to characterize mixed-state entanglement. In the recent literature, two closely related constructions are prominent. The first is the third Rényi negativity,

$$R_3 = -\log \left( \frac{\mathrm{Tr}\left[ (\rho^{T_A})^3 \right]}{\mathrm{Tr} \rho^3} \right),$$

used as a proxy of mixed-state entanglement in finite-temperature many-body systems. The second is the tripartite invariant

$$I_5 = \operatorname{Tr}\left[ (\rho_{BC}^{\Gamma})^3 \right],$$

a permutation-invariant polynomial function of a reduced density matrix that serves as a separability criterion in multipartite settings. Both constructions are based on the third moment of a partially transposed density operator, but they are deployed for different tasks: diagnosing mixed-state entanglement structure in extended systems, resolving disorder-averaged negativity spectra, and providing necessary and sufficient criteria for full separability in tripartite pure and mixed states (Wu et al., 2019, Ma et al., 3 May 2026).

1. Definitions and formal structure

The standard entanglement negativity for a bipartite density matrix $\rho$ on $\mathcal{H}_A \otimes \mathcal{H}_B$ is

$E_N = \log \| \rho^{T_A}\|_1,$

with $\rho^{T_A}$ the partial transpose and $\|\cdot\|_1$ the trace norm. Rényi negativities generalize this through moments of the partial transpose,

$$R_n = -\log \left( \frac{\mathrm{Tr}\left[ (\rho^{T_A})^n \right]}{\mathrm{Tr} \rho^n} \right),$$

and the third Rényi negativity is the $n=3$ case. For pure states, $R_n$ is related to the usual Rényi entanglement entropy; for even $$I_5 = \operatorname{Tr}\left[ (\rho_{BC}^{\Gamma})^3 \right],$$0, analytic continuation as $$I_5 = \operatorname{Tr}\left[ (\rho_{BC}^{\Gamma})^3 \right],$$1 provides the ordinary negativity. Because $$I_5 = \operatorname{Tr}\left[ (\rho_{BC}^{\Gamma})^3 \right],$$2 is odd, $$I_5 = \operatorname{Tr}\left[ (\rho_{BC}^{\Gamma})^3 \right],$$3 is typically used as a computable proxy rather than as the analytic continuation itself (Wu et al., 2019).

In multipartite quantum information, the third-order negativity is formulated as

$$I_5 = \operatorname{Tr}\left[ (\rho_{BC}^{\Gamma})^3 \right],$$4

for a three-qubit or three-qudit pure state $$I_5 = \operatorname{Tr}\left[ (\rho_{BC}^{\Gamma})^3 \right],$$5 with $$I_5 = \operatorname{Tr}\left[ (\rho_{BC}^{\Gamma})^3 \right],$$6. Using the replica trick, the same invariant can be written as

$$I_5 = \operatorname{Tr}\left[ (\rho_{BC}^{\Gamma})^3 \right],$$7

where the $$I_5 = \operatorname{Tr}\left[ (\rho_{BC}^{\Gamma})^3 \right],$$8 are permutation operators acting on replicas of subsystem $$I_5 = \operatorname{Tr}\left[ (\rho_{BC}^{\Gamma})^3 \right],$$9 (Ma et al., 3 May 2026).

A third line of usage appears in the negativity spectrum of disordered systems, where one studies

$\rho$0

The third negativity moment is then $\rho$1, and its scaling encodes information about the full spectrum of the partially transposed reduced density matrix rather than only the logarithmic negativity (Turkeshi et al., 2019).

Quantity	Definition	Primary setting
Third Rényi negativity	$\rho$2	Bipartite mixed states at finite temperature
Third-order negativity	$\rho$3	Tripartite separability of pure and mixed states
Third negativity moment	$\rho$4	Negativity spectrum in disordered systems

2. Finite-temperature many-body entanglement and the third Rényi negativity

In the two-dimensional transverse field Ising model, the third Rényi negativity was computed with quantum Monte Carlo simulations across the finite-temperature phase transition. For local models, the Rényi negativity obeys an area law,

$\rho$5

where $\rho$6 is the boundary length between subsystems and $\rho$7 is the area-law coefficient. The principal numerical result is that the area-law coefficient is singular across the transition, specifically at $\rho$8 for the third Rényi negativity because of the replica construction (Wu et al., 2019).

The singularity is visible in the temperature derivative of the boundary-density contribution. The quantity $\rho$9 displays a cusp at $\mathcal{H}_A \otimes \mathcal{H}_B$0 and scales as $\mathcal{H}_A \otimes \mathcal{H}_B$1 at criticality, mirroring the scaling of the specific heat in the two-dimensional Ising universality class with $\mathcal{H}_A \otimes \mathcal{H}_B$2 and $\mathcal{H}_A \otimes \mathcal{H}_B$3. This identifies a sharply localized entanglement response at the entangling boundary, even though the thermal transition itself has a divergent correlation length (Wu et al., 2019).

The subleading constant $\mathcal{H}_A \otimes \mathcal{H}_B$4 was isolated using a Levin-Wen-type subtraction scheme,

$\mathcal{H}_A \otimes \mathcal{H}_B$5

Quantum Monte Carlo results found $\mathcal{H}_A \otimes \mathcal{H}_B$6 to be zero within statistical error for all studied system sizes and across the phase transition temperature. Each individual $\mathcal{H}_A \otimes \mathcal{H}_B$7 is singular at the transition, but the subtraction yields no net long-range or non-local quantum entanglement contribution. The corresponding interpretation in the paper is explicit: entanglement at the finite-temperature critical point is short-ranged, despite the divergent classical correlation length (Wu et al., 2019).

3. Solvable models, universality, and the absence of long-range entanglement

The same qualitative picture appears in several exactly solvable models. In the quantum spherical model and the two-dimensional Gaussian (free boson) model, the temperature derivative of the Rényi negativity is singular at $\mathcal{H}_A \otimes \mathcal{H}_B$8, in agreement with the finite-temperature Monte Carlo results for the non-integrable two-dimensional transverse field Ising model. In both solvable and non-integrable cases, the long-range part $\mathcal{H}_A \otimes \mathcal{H}_B$9 remains strictly zero in the thermodynamic limit, which the paper identifies as a universal feature of finite-temperature phase transitions in these systems (Wu et al., 2019).

In the solvable Gaussian and mean-field models, the subleading term vanishes exponentially with system size even at the critical point,

$E_N = \log \| \rho^{T_A}\|_1,$0

where $E_N = \log \| \rho^{T_A}\|_1,$1 is a finite quantum correlation length unrelated to the diverging classical correlation length. This supports a separation between classical criticality and the range of quantum entanglement in Gibbs states. A plausible implication, stated in the paper’s interpretive discussion, is that the Gibbs state near or at the transition can be decomposed into minimally entangled typical thermal states (METTS), that is, into a mixture of area-law pure states (Wu et al., 2019).

This framework also addresses a common misconception: long-range correlations at a thermal critical point do not, by themselves, imply long-range quantum entanglement. The singularity in the area-law coefficient is boundary-local, whereas the non-local term vanishes. The paper further argues that earlier linked-cluster claims regarding the absence of area-law singularity are likely artifacts of small system sizes, since the larger-scale quantum Monte Carlo simulations observe the singularity unambiguously (Wu et al., 2019).

4. Third negativity moments in the random singlet phase

In the random singlet phase, the third moment of the partially transposed reduced density matrix is part of the negativity spectrum. For a single disorder realization, the moments depend on the numbers of singlets crossing subsystem boundaries, and for $E_N = \log \| \rho^{T_A}\|_1,$2 one obtains

$E_N = \log \| \rho^{T_A}\|_1,$3

with $E_N = \log \| \rho^{T_A}\|_1,$4 the number of singlets shared between $E_N = \log \| \rho^{T_A}\|_1,$5 and $E_N = \log \| \rho^{T_A}\|_1,$6, and $E_N = \log \| \rho^{T_A}\|_1,$7 the singlets between $E_N = \log \| \rho^{T_A}\|_1,$8 and $E_N = \log \| \rho^{T_A}\|_1,$9 (Turkeshi et al., 2019).

The disordered problem distinguishes two inequivalent averages: $\rho^{T_A}$0 For adjacent intervals of size $\rho^{T_A}$1, the average of the log for the third moment obeys

$\rho^{T_A}$2

whereas the log of the average obeys

$\rho^{T_A}$3

The two coefficients are therefore distinct. This is one of the paper’s central results: negativity and logarithmic negativity are not trivially related after the average over the disorder, and higher negativity moments reveal finer information about the spectrum than the standard negativity alone (Turkeshi et al., 2019).

The scaling coefficients are universal in the sense used by the random-singlet literature: they do not depend on microscopic disorder details once the system flows to the infinite-randomness fixed point. Analytic predictions from strong disorder renormalization group were checked against SDRG numerics and exact computations for the random XX chain, with excellent agreement for both $\rho^{T_A}$4 and $\rho^{T_A}$5. The paper characterizes the inequivalence between different disorder averages as genuine multifractal, or “multiscaling,” behavior (Turkeshi et al., 2019).

5. Tripartite separability and the invariant $\rho^{T_A}$6

For normalized tripartite pure states, the principal separability statement is exact: $\rho^{T_A}$7 This is the main proposition in the 2026 work on multipartite measures. If $\rho^{T_A}$8, the state is a product across all three parties $\rho^{T_A}$9; if the state is fully separable, then direct computation gives $\|\cdot\|_1$0. The proof uses properties of the partial transpose and the spectrum of $\|\cdot\|_1$1, showing that only product states saturate the maximal value (Ma et al., 3 May 2026).

For mixed states, the paper introduces convex-roof extensions. One measure is based on

$\|\cdot\|_1$2

with

$\|\cdot\|_1$3

A second measure is

$\|\cdot\|_1$4

where the $\|\cdot\|_1$5 exponent ensures monotonicity as an entanglement measure. The vanishing criteria are again exact: $\|\cdot\|_1$6 and

$\|\cdot\|_1$7

Within the paper’s framing, third-order negativity thus extends the role played by bipartite negativity and the PPT criterion into a compact, permutation-invariant criterion for full separability in three-party systems (Ma et al., 3 May 2026).

This separates two notions that are often conflated. Bipartite negativity is a criterion on a chosen cut; third-order negativity in the $\|\cdot\|_1$8 sense is a multipartite criterion for full product structure. The latter is designed to bypass the combinatorial complexity of testing all bipartitions individually in multipartite systems (Ma et al., 3 May 2026).

6. Multipartite generalization, four-qubit structure, and conformal field theory

The replica-based formulation of $\|\cdot\|_1$9 generalizes to a hierarchy of multipartite invariants. For an $$R_n = -\log \left( \frac{\mathrm{Tr}\left[ (\rho^{T_A})^n \right]}{\mathrm{Tr} \rho^n} \right),$$0-partite system, the paper constructs

$$R_n = -\log \left( \frac{\mathrm{Tr}\left[ (\rho^{T_A})^n \right]}{\mathrm{Tr} \rho^n} \right),$$1

where the $$R_n = -\log \left( \frac{\mathrm{Tr}\left[ (\rho^{T_A})^n \right]}{\mathrm{Tr} \rho^n} \right),$$2 are distinct permutation operators on the replicas associated with each subsystem. The general proposition is

$$R_n = -\log \left( \frac{\mathrm{Tr}\left[ (\rho^{T_A})^n \right]}{\mathrm{Tr} \rho^n} \right),$$3

The paper states that this yields a permutation-symmetric necessary and sufficient criterion for full separability in arbitrary multipartite systems, including arbitrary qudit dimensions, and that a convex-roof extension preserves its role as an entanglement monotone and separability criterion for mixed states (Ma et al., 3 May 2026).

A concrete organizational example is the four-qubit pure state, which has 18 independent parameters and a corresponding set of 18 entanglement measures: 6 two-tangles, 8 tripartite measures, and 4 quadripartite measures. Among the tripartite measures are four $$R_n = -\log \left( \frac{\mathrm{Tr}\left[ (\rho^{T_A})^n \right]}{\mathrm{Tr} \rho^n} \right),$$4-type quantities derived from $$R_n = -\log \left( \frac{\mathrm{Tr}\left[ (\rho^{T_A})^n \right]}{\mathrm{Tr} \rho^n} \right),$$5, one for each three-party reduction. Their vanishing is equivalent to full separability of the corresponding three-party mixed state. If all four vanish, only the four-tangle may remain nonzero, corresponding to generalized GHZ-type states (Ma et al., 3 May 2026).

The same paper also places third-order negativity in quantum field theory. In CFT$$R_n = -\log \left( \frac{\mathrm{Tr}\left[ (\rho^{T_A})^n \right]}{\mathrm{Tr} \rho^n} \right),$$6, the third-order negativity corresponds to a three-point function of twist operators, and in the AdS/CFT correspondence this three-point function can be evaluated using geometric data and OPE coefficients. The paper further states that the universal value for adjacent intervals in CFT$$R_n = -\log \left( \frac{\mathrm{Tr}\left[ (\rho^{T_A})^n \right]}{\mathrm{Tr} \rho^n} \right),$$7, at large central charge, can be determined. This links third-order negativity to field-theoretic observables rather than treating it solely as an abstract invariant (Ma et al., 3 May 2026).