---
title: Relativistic Flux Tube Model in QCD
url: https://www.emergentmind.com/topics/relativistic-flux-tube-model
type: topic
---

# Relativistic Flux Tube Model in QCD

The relativistic flux tube model denotes a family of QCD-inspired descriptions in which confinement or early-time color dynamics is represented by a string-like or tube-like distribution of color flux with relativistic degrees of freedom. In hadron spectroscopy, it typically models mesons or quark–diquark baryons as relativistic endpoints connected by a rotating flux tube of tension $\sigma$ or $\tau$, with spin splittings added perturbatively; in confining-gauge-theory studies, it appears as the effective string description of long closed flux tubes and their world-sheet $S$-matrix; and in heavy-ion applications, it describes longitudinal color-electric flux tubes that decay through the Schwinger mechanism and couple self-consistently to kinetic transport [2509.14699] [1906.08098] [1505.08081]. Across these settings, the common structure is the replacement of purely static confinement by a dynamical relativistic object that carries energy, momentum, angular momentum, and, in some formulations, explicit particle-production or resonance degrees of freedom.

## 1. Conceptual domain and physical content

In the spectroscopy literature, the model starts from a heavy–light hadron, quarkonium, or quark–diquark baryon represented as two color sources connected by a gluonic flux tube. The flux tube has constant tension, contributes directly to the total energy and orbital angular momentum, and enforces Regge-like behavior through relativistic rotation of the endpoints [2509.14699] [1910.06065]. In heavy–light implementations, the light endpoint is taken to be effectively ultrarelativistic while the heavy endpoint is slower; in bottomonium the two endpoints are equal-mass heavy quarks; in singly and doubly heavy baryons the system is recast as a heavy quark or heavy diquark bound to a light diquark or light quark [2407.01655] [2306.06349].

A second usage concerns confining flux tubes themselves as physical extended objects in Yang–Mills theory. There the relevant degrees of freedom are the transverse Goldstone modes of a long string, and the low-energy dynamics is a relativistic two-dimensional quantum field theory on the world-sheet. The effective action is the Nambu–Goto action plus irrelevant operators such as $K^4$ terms, and the corresponding $2\to2$ massless $S$-matrix can be constrained nonperturbatively by analyticity, crossing, and unitarity [1906.08098]. Lattice calculations of closed winding flux tubes in $D=3+1$ $SU(N)$ gauge theory show that many low-lying states are accurately described by Nambu–Goto even at surprisingly short lengths, while some $0^-$ channels deviate strongly and are plausibly associated with massive internal modes [1007.4720].

A third usage appears in non-equilibrium QCD. There the flux tube is not a bound-state string but a longitudinal color-electric field generated in the overlap region of ultrarelativistic nuclei. In that setting the field decays into quarks, antiquarks, and gluons through the Schwinger mechanism, and the backreaction of the produced plasma modifies the field through Maxwell-like equations in an abelianized approximation [1505.08081] [1605.04126]. This suggests that “relativistic flux tube model” is best understood as a unifying geometric and dynamical motif rather than as a single universally fixed formalism.

## 2. Rotating-string dynamics in hadron spectroscopy

The spectroscopy version is built from a straight rotating string with massive endpoints. In one heavy–light implementation, the spin-averaged mass $\bar{M}$ obeys the Regge-like relation
\[
(\bar{M} - m_h)^2 \;=\; \sigma\,\frac{L + \lambda\,n_r}{2}\;+\;\big(m_l + m_h\,v_h^2\big)^2,
\]
with
\[
r \;=\; (v_l + v_h)\,\sqrt{\frac{8\,(L + \lambda\,n_r)}{\sigma}},
\]
where $L$ is the orbital angular momentum, $n_r=n-1$ is the radial quantum number, $\lambda$ controls the spacing of radial daughters, and $v_l$, $v_h$ are the endpoint velocities [2509.14699]. An equivalent heavy–light formulation for open-charm hadrons writes the same structure as the working spin-averaged mass formula and derives it from a rigidly rotating Nambu–Goto-like string with endpoint masses $m_{0,h}$ and $m_{0,l}$ [2509.24772]. For singly charmed and singly bottom baryons in the quark–diquark picture, the same dependence is retained with the heavy constituent identified as a charm or bottom quark and the light constituent as a scalar or axial diquark [2306.06349] [2407.01655].

Bottomonium uses a related but distinct analytic formula. Starting from the exact energy and angular momentum of a straight rotating string with equal endpoint masses, one obtains for the spin-averaged mass
\[
M_{nL}\equiv \bar{\mathrm{M}(nL)}=2m_b+\left(\frac{\sigma^2}{2\pi^2m_b}\right)^{1/3}(\lambda n+L)^{2/3},
\]
after semiclassical elimination of the endpoint velocity and promotion of $L\to \lambda n+L$ to incorporate radial excitations [1910.06065]. In that work the parameters are fixed from $\bar{\mathrm{M}(1S)}$, $\bar{\mathrm{M}(2S)}$, and $h_b(1P)$, yielding
\[
m_b=4.7224~\mathrm{GeV},\qquad \sigma=2.96~\mathrm{GeV}^2,\qquad \lambda=1.41
\]
for the spin-averaged bottomonium spectrum [1910.06065].

These formulations share two structural consequences. First, the flux tube carries angular momentum, so confinement is not reduced to a static linear potential. Second, orbital and radial excitations organize into approximately linear, parallel, and equidistant Regge patterns when plotted in suitable variables, a property repeatedly emphasized for bottom mesons, singly charmed baryons, singly bottom baryons, and doubly heavy baryons [2509.14699] [2306.06349] [2407.01655].

## 3. Spin splittings, coupling schemes, and spectroscopic applications

Physical masses are obtained by adding spin-dependent corrections to the spin-averaged flux-tube result. In heavy–light mesons the correction is written as
\[
\Delta M \;=\; H_{so} \;+\; H_t \;+\; H_{ss},
\]
with
\[
\begin{aligned}
H_{so} &= \Bigg[\left(\frac{2\,\alpha}{3\,r^3} - \frac{b_0}{2\,r}\right)\frac{1}{m_l^2}
+ \frac{4\,\alpha}{3\,r^3}\frac{1}{m_h\,m_l}\Bigg]\;\mathbf{L}\cdot\mathbf{S}_l \\
&\quad+ \Bigg[\left(\frac{2\,\alpha}{3\,r^3} - \frac{b_0}{2\,r}\right)\frac{1}{m_h^2}
+ \frac{4\,\alpha}{3\,r^3}\frac{1}{m_h\,m_l}\Bigg]\;\mathbf{L}\cdot\mathbf{S}_h,
\end{aligned}
\]
\[
H_t \;=\; \frac{4\,\alpha}{3\,r^3\,m_h\,m_l}\;\hat{\mathbf{B}},
\qquad
H_{ss} \;=\; \frac{32\pi\,\alpha\,\sigma_0^3}{9\sqrt{\pi}\,m_h\,m_l}\;e^{-\sigma_0^2 r^2}\;\mathbf{S}_l\cdot\mathbf{S}_h,
\]
and matrix elements evaluated in the $|J,j\rangle$ basis appropriate to $j$–$j$ coupling [2509.14699]. The same operator structure is used for doubly charmed and doubly bottom baryons after replacing the heavy endpoint by an axial-vector diquark [2509.24772] [2509.14699].

For singly charmed and singly bottom baryons in the quark–diquark picture, the heavy-quark limit motivates coupling the orbital motion first to the light-diquark spin, $j_\ell=\mathbf{L}+\mathbf{S}_{\mathcal D}$, and then coupling $j_\ell$ to the heavy-quark spin. The spin splittings are again decomposed into spin–orbit, tensor, and contact terms, but the relevant matrix elements differ between scalar-diquark antitriplets and axial-diquark sextets [2306.06349] [2407.01655]. This machinery is used to assign several experimentally observed states. For example, the bottom-meson analysis identifies $B_J^*(5732)$ as a candidate for $1P(1^+,j=1/2)$, $B_J(5840)$ as the first radial excitation $2S$, and $B_J(5970)$ as a $1D$-wave orbital excitation, while the bottom-strange sector places $B_{sJ}(6063)$ and $B_{sJ}(6114)$ in the $1D$ multiplet [2509.14699]. The singly bottom baryon study identifies $\Sigma_b(6097)$ as a $1P(3/2^-)$ state, $\Xi_b(6100)$ as a $1P(1/2^-)$ state of $\Xi_b$, and $\Xi_b(6227)$ as a $1P$ excitation of $\Xi_b'$ with $J^P=3/2^-$, combining mass and strong-decay information [2407.01655].

In open charm, the same framework reproduces the established $D$ and $D_s$ low-lying states and is used to assign higher resonances such as $D_2(2740)^0$, $D_3^*(2750)$, $D_{s1}^*(2860)^\pm$, $D_{s3}^*(2860)^\pm$, and $D_{sJ}(3040)^+$ [2509.24772]. At the same time, the anomalously light $D_{s0}(2317)$ and $D_{s1}(2460)$ are explicitly noted as not aligning with conventional states in this framework, and coupled-channel, molecular, or tetraquark interpretations are discussed as plausible resolutions [2509.24772]. A similar caveat appears for bottomonium: $\Upsilon(10750)$, $\Upsilon(10860)$, and $\Upsilon(11020)$ are interpreted respectively as predominantly $3^3D_1$, $5^3S_1$, and $4^3D_1$, while threshold effects are invoked to explain tensions in the higher vector spectrum [1910.06065].

## 4. Flux tubes as strings: closed spectra and world-sheet bootstrap

In confining gauge theories, the relativistic flux tube also appears as a long string whose low-energy degrees of freedom are world-sheet Goldstone bosons. The effective action is written schematically as
\[
A = \int d^2\sigma \sqrt{-h}\, [ \ell_s^{-2} + \mathcal{R}(h) + K^2 + \ell_s^2 K^4 + \dots ],
\]
where $h_{\alpha\beta}$ is the induced metric and the first non-universal corrections arise from two independent $K^4$ structures with Wilson coefficients $\alpha_3$ and $\beta_3$ [1906.08098]. In two-dimensional left–right scattering, the low-energy phase shifts satisfy
\[
2 \delta_{sym}  =  s/4  + \alpha_2 s^2 + \alpha_3 s^3 + O(s^4),
\]
\[
2 \delta_{anti} =  s/4 - \alpha_2 s^2 + (\alpha_3 + 2 \beta_3) s^3 + O(s^4),
\]
\[
2 \delta_{sing} =  s/4 - (D-3) \alpha_2 s^2 + (\alpha_3 - (D-2)\beta_3) s^3 + O(s^4),
\]
with universal $\alpha_2=(D-26)/(384\pi)$ [1906.08098]. The bootstrap then constrains $\alpha_3$ and $\beta_3$ by analyticity and unitarity. In $D=4$, for example, one finds
\[
\alpha_3 \ge -\frac{1}{768} + \frac{121}{9216\pi^2},
\qquad
\alpha_3 \ge -\frac{1}{768} + |\beta_3|,
\]
and boundary solutions exhibit a sharp antisymmetric resonance whose mass and width match the world-sheet axion inferred on the lattice [1906.08098].

Lattice calculations of closed flux tubes winding around a spatial torus provide a complementary nonperturbative test. The Nambu–Goto prediction for a closed string of length $l$ is
\[
E^2(l) = (K l)^2 + 8\pi K \left[N_L + N_R - \frac{D-2}{24}\right] + \left(\frac{2\pi q}{l}\right)^2,
\]
with $K\equiv \sigma$ and level matching $N_L-N_R=q$ [1007.4720]. For many low-lying states in $D=3+1$ $SU(N)$ gauge theories, the spectrum follows this Nambu–Goto form down to $l\sqrt{\sigma}\approx 2$–$3$, even where the naive $1/l$ expansion no longer converges [1007.4720]. The major exception is a family of $0^-$ states lying anomalously low. A simple relativistic massive-mode ansatz,
\[
E(p) \approx E_0(l) + \sqrt{m^2 + p^2},
\]
with $m/\sqrt{\sigma}\approx 1.85$, qualitatively describes these deviations and suggests that some flux-tube channels probe internal massive structure rather than only transverse Goldstone motion [1007.4720]. This is closely aligned with the bootstrap picture in which boundary $S$-matrices develop resonance structure, including the world-sheet axion and broader states in other irreducible channels [1906.08098].

## 5. QCD realizations beyond spectroscopy: compactification and heavy-ion flux tubes

A more field-theoretic realization derives an effective $1+1$-dimensional theory inside a flux-tube environment. Starting from the QCD$_4$ action and imposing longitudinal dominance, transverse confinement, and a factorized ansatz for the quark field,
\[
\Psi(x)=\frac{1}{\sqrt{2}}
\begin{pmatrix}
G_1(\mathbf{r}_\perp)\,[f_+(X)+f_-(X)]\\
-\,G_2(\mathbf{r}_\perp)\,[f_+(X)-f_-(X)]\\
G_1(\mathbf{r}_\perp)\,[f_+(X)-f_-(X)]\\
G_2(\mathbf{r}_\perp)\,[f_+(X)+f_-(X)]
\end{pmatrix},
\qquad X=(x^0,x^3),
\]
one obtains a compactified QCD$_2$ action with effective coupling
\[
g_{2D}=\int dx^1dx^2\,g_{4D}\,\big(|G_1|^2+|G_2|^2\big)^{3/2},
\]
and effective quark mass
\[
m_{q,T}=\lambda,
\]
where $\lambda$ is the transverse eigenvalue [1212.3301]. The gauge field acquires a transverse-confinement mass term as well, and the resulting one-dimensional gluodynamics yields linear confinement with string tension proportional to $g_{2D}^2$ [1212.3301]. This provides a first-principles bridge between flux-tube geometry and the effective parameters that spectroscopy models usually treat phenomenologically.

In ultrarelativistic heavy-ion collisions, the “color-flux-tube model” describes the early-time system by a longitudinal color-electric field which decays by Schwinger tunneling. In the transport formulation, the field obeys
\[
\partial_\mu F^{\mu\nu}=J^\nu,
\]
or, in Bjorken coordinates,
\[
\frac{dE}{d\tau} = \rho \sinh\eta_s - j \cosh\eta_s,
\]
while the particle distributions satisfy a relativistic Boltzmann equation with Lorentz force, source, and collision terms,
\[
p^\mu \partial_\mu f_{jc}(x,p) + g Q_{jc} F^{\mu\nu} p_\nu \frac{\partial f_{jc}}{\partial p^\mu}
= S_{jc}(x,p) + C_{jc}[f]
\]
[1505.08081]. Pair production is implemented through the Schwinger source, for example in the static-box frame,
\[
\frac{dN_{jc}}{d^4x\, d^2p_T\, dp_z} = R_{jc}(p_T)\,\delta(p_z)\,p_0,
\]
with
\[
R_{jc}(p_T) = (\mathcal{E}_{jc}/4\pi^3)\, \left| \ln(1 \pm e^{-\pi p_T^2/\mathcal{E}_{jc}}) \right|,
\]
and the total current contains both matter and displacement components [1505.08081]. For small $\eta/s$ the model yields $\tau_{\text{isotropization}}\approx 0.8$ fm/c and $\tau_{\text{thermalization}}\approx 1$ fm/c, while the initial color-electric field decays within $\sim 1$ fm/c almost irrespective of viscosity; at large $\eta/s$, persistent plasma oscillations modulate the late-time pressure anisotropy [1505.08081]. A related review emphasizes that for $\eta/s=1/(4\pi)$ the transverse-momentum slopes are close to equilibrium within $\tau\lesssim 1$ fm/c, and for $\eta/s\lesssim 3/(4\pi)$ viscous hydrodynamics already describes the system well, indicating rapid hydrodynamization [1605.04126].

These heavy-ion formulations are physically distinct from bound-state RFT spectroscopy, but they retain the same central idea: a relativistic color-electric tube stores energy, accelerates color charges, and transfers field energy into matter through explicitly relativistic dynamics [1505.08081] [1605.04126].

## 6. Lattice calibration, phenomenological scope, and limitations

Lattice studies of static quark–antiquark flux tubes at zero temperature provide quantitative scales for the transverse structure that underlies the flux-tube picture. In both $SU(3)$ pure gauge theory and $(2+1)$-flavor QCD, the longitudinal chromoelectric field in the mid-plane is well fitted by the dual Abrikosov/Clem profile up to about $d\simeq 1.5$ fm in pure gauge and $d\simeq 1.1$ fm with dynamical quarks [1702.06437]. The extracted penetration length is approximately $\lambda\simeq 0.15$–$0.17$ fm, the coherence length is of order $0.4$–$1.0$ fm in pure gauge and $0.4$–$0.74$ fm in $(2+1)$-flavor QCD, the Ginzburg–Landau parameter satisfies $\kappa<1/\sqrt{2}$ throughout the measured range, and the transverse width is approximately $0.42$–$0.53$ fm [1702.06437]. The study therefore supports a type-I dual-superconductor description of the static tube over the distances probed and supplies concrete scales for any phenomenological flux-tube model [1702.06437].

The model’s phenomenological reach extends beyond conventional hadrons. A pentaquark study treats exotic states as two compact clusters connected by a rotating string, with finite current-quark masses at the endpoints. The resulting Regge trajectories in $J$ versus $M^2$ become non-linear, with effective slopes reduced relative to the massless-end value
\[
\alpha' = \frac{1}{2\pi\sigma},
\]
and with curvature increasing for heavier endpoint configurations [1706.07223]. In heavy-ion phenomenology, flux-tube networks with $3$-point junctions are also used to argue that baryon-number correlations on the last scattering surface may preserve the pattern of flux-tube vertices near deconfinement, yielding short-range sign-alternating baryon–antibaryon correlation structures [1111.0177].

Several limitations recur across the literature. Spectroscopic implementations often do not solve a full dynamical wave equation for the flux tube and instead encode radial excitations by a fitted replacement $L\to L+\lambda n_r$ or $\lambda n + L$ [2509.14699] [1910.06065]. Coupled-channel effects, threshold mixing, and strong-decay mass shifts are generally omitted, although they are explicitly identified as plausible sources of discrepancies for states such as $D_{s0}(2317)$, $D_{s1}(2460)$, $B_J(5970)$, and some $B_s$ excitations [2509.24772] [2509.14699]. In baryon applications, the quark–diquark approximation neglects internal diquark excitations and possible three-body flux configurations [2306.06349] [2407.01655]. In heavy-ion applications, the field is abelianized, chromomagnetic components are neglected, and transport is imposed at fixed $\eta/s$ rather than derived from microscopic QCD [1505.08081]. World-sheet EFT and lattice descriptions, conversely, are rigorous for long confining strings but are not by themselves constituent models of finite hadrons [1906.08098] [1007.4720].

Taken together, these strands indicate that the relativistic flux tube model is most successful when used as a controlled representation of long-range color dynamics with explicit relativistic kinematics. Its most robust outputs are Regge organization, qualitative and often quantitative mass systematics, nontrivial world-sheet constraints, and dynamical field-to-particle conversion in strong longitudinal fields. Its least controlled sector remains the treatment of short-distance physics, threshold effects, and internal structure beyond the leading string degree of freedom [2509.14699] [1906.08098] [1505.08081].

Source: https://www.emergentmind.com/topics/relativistic-flux-tube-model