Flux Tube Entanglement Entropy (FTE²)

• Flux Tube Entanglement Entropy (FTE²) is a gauge-invariant measure quantifying the excess entanglement from the color flux tube connecting a static quark–antiquark pair in Yang–Mills theory.
• It decomposes into internal (color) entropy from Gauss’s law constraints and vibrational entropy from transverse string fluctuations, computed via lattice gauge techniques.
• FTE² scales with entangling region geometry and flux tube topology, offering insights into confinement, nonlocal quantum correlations, and glueball mass relationships.

Flux Tube Entanglement Entropy (FTE$^2$) quantifies the excess entanglement entropy attributed solely to the color flux tube connecting a static quark–antiquark pair in Yang–Mills gauge theory, relative to the vacuum. Distinguished by its gauge-invariant formulation and independence from center-algebra ambiguities, FTE$^2$ isolates nonlocal quantum correlations intrinsic to the confining gauge sector. Its construction and interpretation leverage lattice gauge theory, the replica trick, and operator-based approaches, revealing two fundamental contributions: internal (color) entropy tied to Gauss’s-law constraints, and vibrational (string) entropy associated with low-energy transverse fluctuations. The interplay of entangling region geometry, color group rank, and flux-tube topology governs the behavior of FTE$^2$, rendering it a probe of both the physical and topological structure of non-Abelian flux tubes.

1. Formal Definition and Gauge-Invariant Construction

For a spatial region $V$ within a Yang–Mills lattice containing a static $Q\bar Q$ pair, let $\rho_{V|Q\bar Q}$ denote the reduced density matrix on $V$ with the flux tube present, and $\rho_V^0$ that in the vacuum. The von Neumann FTE$^2$ is defined by

$$\mathrm{FTE}^2 = S_V|Q\bar Q - S_V^0 = -\mathrm{Tr}(\rho_{V|Q\bar Q}\ln\rho_{V|Q\bar Q}) + \mathrm{Tr}(\rho_V^0\ln\rho_V^0).$$

Similarly, the $q$th Rényi FTE$^2$ is

$$\mathrm{FTE}^{2(q)} = S_V^{(q)}|Q\bar Q - S_V^{(q)}, \quad S_V^{(q)} = \frac{1}{1-q}\ln\mathrm{Tr}[\rho_V^q].$$

On the lattice, using the replica trick, FTE$^2$ is expressed in terms of ratios of Polyakov-loop correlators embedded on a $q$-sheeted manifold: $$\mathrm{FTE}^{2(q)} = -\frac{1}{q-1}\ln\left( \frac{ \langle\prod_{r=1}^q P^{(r)}(x_Q)P^{(r)\dagger}(x_{\bar Q})\rangle }{ [\langle P(x_Q)P^\dagger(x_{\bar Q})\rangle ]^q }\right),$$ with $P(x)$ the traced Polyakov line at $x$. The continuum (von Neumann) limit is recovered via $\lim_{q\to1}$.

This construction is manifestly gauge-invariant and eliminates UV-divergent boundary contributions, yielding a finite result in the continuum limit and removing dependence on the choice of center algebra (Amorosso et al., 2024, Amorosso et al., 23 Jan 2026, Amorosso et al., 12 Feb 2025).

2. Decomposition: Internal (Color) and Vibrational Entropy

FTE$^2$ naturally splits into two additive components: $$\mathrm{FTE}^2 \equiv S_{\mathrm{internal}} + S_{\mathrm{vibrational}}.$$

Internal Color Entropy arises from non-Abelian gauge constraints (Gauss’s law) at boundaries where the flux tube is severed by the entangling surface. For $SU(N_c)$ with sources in representation $R$, this takes the universal form

$S_{\mathrm{internal}} = k\ln d_R,$

where $k$ is the number of full boundary crossings of the flux tube, and $d_R$ is the dimension of $R$ (Amorosso et al., 2024, Amorosso et al., 12 Feb 2025).

Vibrational Entropy encodes entanglement from transverse string fluctuations. In (2+1)D, this is modeled by the ground-state entropy of a thin relativistic string: $$S_{\mathrm{vibrational}}\sim \frac{1}{3}\ln(L/\epsilon),$$ for the von Neumann entropy and string of length $L$, with UV cutoff $\epsilon$. For second Rényi entropy,

$$S_{\mathrm{vibrational}}^{(2)}\sim \frac{1}{4}\ln(L/\epsilon),$$

where the fractional subsegment length governs the detailed dependence (Amorosso et al., 2024).

In (1+1)D, transverse vibrations are absent, and FTE$^2$ reduces strictly to $S_{\mathrm{internal}}$ (Amorosso et al., 2024).

3. Geometry, Topology, and Scaling

The entangling region $V$ determines the FTE$^2$ profile through its topology and intersection with the flux tube:

• Half-slab geometry: $V$ is a “slab” of width $w$ crossing or parallel to the flux tube. Translating $V$ across the flux tube generates a broadened step function in $\mathrm{FTE}^2(x_0)$, accurately fit by an error function whose width encodes the intrinsic flux-tube thickness. The asymptotic value at full cross-cutting is $\mathrm{FTE}^2\to 2\ln N_c$ (for $SU(N_c)$ fundamental sources) (Amorosso et al., 2024, Amorosso et al., 12 Feb 2025).
• Multiple/slab and staggered geometries: For $k$ full intersections (e.g., two disconnected slabs), $S_{\mathrm{internal}}\rightarrow k\ln N_c$, with sharp jumps as regions merge to sever the flux tube topologically. Partial crossings—where the region does not fully disconnect the flux tube—do not contribute (Amorosso et al., 23 Jan 2026, Amorosso et al., 12 Feb 2025).
• Area-law scaling: Lattice studies in quenched QCD confirm that the excess (Rényi) FTE$^2$ scales linearly with the boundary area overlapping the flux tube, $\Delta S_2 \propto L_z\,w$, with proportionality constant determined empirically (Amorosso et al., 2023).

This dependence on topology certifies that $S_{\mathrm{internal}}$ reflects bulk, not boundary-local, entanglement, distinguishing it from UV-sensitive area-law terms in traditional entropy measures.

4. Physical Scales: Entanglement Radius and Intrinsic Width

Recent investigations have identified two important scales characterizing the flux tube's spatial structure:

• Entanglement Radius ($\xi_0$): The minimal transverse width required for the entangling region to register a full color entropy contribution is nonzero, reflecting the intrinsic finite width of the flux tube. $\xi_0$ grows approximately linearly with $N_c$ and is independent of Rényi order and interquark separation. (Amorosso et al., 23 Jan 2026)
• Intrinsic Width ($\lambda$): The tails of the FTE$^2$ transverse profile decay exponentially with characteristic scale $\lambda$, matching the inverse mass of the lightest $0^{++}$ glueball, i.e., $\lambda \sim 1/m_{0^{++}}$. The table below summarizes extracted values for several $SU(N_c)$ groups:

$N_c$	$\xi_0\sqrt{\sigma_0}$	$\lambda\sqrt{\sigma_0}$	$\sqrt{\sigma_0}/m_{0^{++}}$
2	0.185(6)	0.223(15)	0.212(2)
3	0.269(8)	0.199(45)	0.231(2)
4	0.321(8)	0.218(31)	0.236(3)
5	0.393(9)	0.191(29)	0.239(3)

This establishes the entanglement radius as a new, universal physical scale governing topological color entanglement.

5. Dimensional Reduction and Analytic Results in (1+1)D

In (1+1)D Yang–Mills theory, the absence of string vibrational modes enables analytic evaluation of FTE$^2$: $S_{\mathrm{FTE}^2} = k\,\ln d_R,$ with $k$ the number of boundary crossings and $d_R$ the source representation’s dimension. This result is independent of flux-tube length, lattice spacing, or placement of replica branch points (Amorosso et al., 2024, 1705.01549). The ground-state entanglement structure splits naturally into a classical (Shannon) term, a color entanglement term, and a Bell-pair (genuine) entanglement term, with the color piece universal and robust in the continuum limit (1705.01549).

6. Lattice Implementation and Computational Strategies

FTE$^2$ is evaluated through lattice simulations using the Polyakov-loop–based replica construction, exploiting periodicity and boundary conditions to reconstruct the reduced density matrices for the chosen region $V$ (Amorosso et al., 2024, Amorosso et al., 23 Jan 2026, Amorosso et al., 12 Feb 2025). Calculation steps include:

• Implementation of heatbath and over-relaxation updates for the gauge fields (Wilson plaquette action), supplemented by multilevel algorithms for high-precision Polyakov-loop correlator measurement.
• Formulation of entangling geometries (slab, half-slab, double-slab) to probe spatial dependence and topology.
• UV divergence elimination via vacuum subtraction.
• Control over the number and location of region–complement intersections to isolate individual $S_{\mathrm{internal}}$ contributions.

Numerical results confirm UV finiteness, scaling with $N_c$, the topological nature of the color contribution, and quantitative agreement with analytic (1+1)D and string-theory-based vibrational models (Amorosso et al., 2024, Amorosso et al., 12 Feb 2025).

7. Topological Nature, Open Problems, and Extensions

The color entanglement revealed by FTE$^2$ displays a distinctly topological character: only when the entangling region $V$ fully severs the confining flux tube does $S_{\mathrm{internal}}$ accrue a quantized contribution per crossing; partial or incomplete intersections yield little or no entropy increase (Amorosso et al., 23 Jan 2026, Amorosso et al., 12 Feb 2025). This marks FTE$^2$ as a probe of bulk, distillable entanglement in gauge theories.

Unresolved issues include:

• The precise behavior of FTE$^2$ for general, possibly non-simply connected regions.
• Refinements required for the thin-string vibrational approximation, especially in the presence of finite intrinsic width.
• The form of $S_{\mathrm{internal}}$ in the presence of dynamical quarks and in higher dimensions, for which string breaking and richer topological possibilities may play a role.
• The quantitative relationship between FTE$^2$ and glueball masses, and its potential connection to other entanglement measures in gauge theories.

Developments in FTE$^2$ thus offer new quantitative and topological diagnostics for understanding confinement, flux-tube structure, and quantum correlations in non-Abelian gauge theory (Amorosso et al., 2024, Amorosso et al., 2024, Amorosso et al., 23 Jan 2026, Amorosso et al., 12 Feb 2025, Amorosso et al., 2023).