Relativistic Flux Tube Model in QCD

• The relativistic flux tube model is a QCD-inspired framework that represents confinement as a rotating string-like flux tube connecting color sources.
• It dynamically incorporates relativistic kinematics to produce Regge trajectories and mass systematics in hadron spectroscopy, validated by lattice and phenomenological studies.
• Applications span effective 1+1D QCD and heavy-ion collisions where flux tubes decay via the Schwinger mechanism, illustrating energy-to-particle conversion.

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 (Jakhad et al., 18 Sep 2025, Miro et al., 2019, Ruggieri et al., 2015). 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 (Jakhad et al., 18 Sep 2025, Chen et al., 2019). 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 (Jakhad et al., 2024, Jakhad et al., 2023).

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 (Miro et al., 2019). 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 (Athenodorou et al., 2010).

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 (Ruggieri et al., 2015, Ryblewski, 2016). 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

$\tau$0

with

$\tau$1

where $\tau$2 is the orbital angular momentum, $\tau$3 is the radial quantum number, $\tau$4 controls the spacing of radial daughters, and $\tau$5, $\tau$6 are the endpoint velocities (Jakhad et al., 18 Sep 2025). 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 $\tau$7 and $\tau$8 (Jakhad et al., 29 Sep 2025). 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 (Jakhad et al., 2023, Jakhad et al., 2024).

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

$\tau$9

after semiclassical elimination of the endpoint velocity and promotion of $S$0 to incorporate radial excitations (Chen et al., 2019). In that work the parameters are fixed from $S$1, $S$2, and $S$3, yielding

$S$4

for the spin-averaged bottomonium spectrum (Chen et al., 2019).

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 (Jakhad et al., 18 Sep 2025, Jakhad et al., 2023, Jakhad et al., 2024).

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

$S$5

with

$S$6

$S$7

and matrix elements evaluated in the $S$8 basis appropriate to $S$9–$K^4$0 coupling (Jakhad et al., 18 Sep 2025). The same operator structure is used for doubly charmed and doubly bottom baryons after replacing the heavy endpoint by an axial-vector diquark (Jakhad et al., 29 Sep 2025, Jakhad et al., 18 Sep 2025).

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, $K^4$1, and then coupling $K^4$2 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 (Jakhad et al., 2023, Jakhad et al., 2024). This machinery is used to assign several experimentally observed states. For example, the bottom-meson analysis identifies $K^4$3 as a candidate for $K^4$4, $K^4$5 as the first radial excitation $K^4$6, and $K^4$7 as a $K^4$8-wave orbital excitation, while the bottom-strange sector places $K^4$9 and $2\to2$0 in the $2\to2$1 multiplet (Jakhad et al., 18 Sep 2025). The singly bottom baryon study identifies $2\to2$2 as a $2\to2$3 state, $2\to2$4 as a $2\to2$5 state of $2\to2$6, and $2\to2$7 as a $2\to2$8 excitation of $2\to2$9 with $S$0, combining mass and strong-decay information (Jakhad et al., 2024).

In open charm, the same framework reproduces the established $S$1 and $S$2 low-lying states and is used to assign higher resonances such as $S$3, $S$4, $S$5, $S$6, and $S$7 (Jakhad et al., 29 Sep 2025). At the same time, the anomalously light $S$8 and $S$9 are explicitly noted as not aligning with conventional states in this framework, and coupled-channel, molecular, or tetraquark interpretations are discussed as plausible resolutions (Jakhad et al., 29 Sep 2025). A similar caveat appears for bottomonium: $D=3+1$0, $D=3+1$1, and $D=3+1$2 are interpreted respectively as predominantly $D=3+1$3, $D=3+1$4, and $D=3+1$5, while threshold effects are invoked to explain tensions in the higher vector spectrum (Chen et al., 2019).

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

$D=3+1$6

where $D=3+1$7 is the induced metric and the first non-universal corrections arise from two independent $D=3+1$8 structures with Wilson coefficients $D=3+1$9 and $SU(N)$0 (Miro et al., 2019). In two-dimensional left–right scattering, the low-energy phase shifts satisfy

$SU(N)$1

$SU(N)$2

$SU(N)$3

with universal $SU(N)$4 (Miro et al., 2019). The bootstrap then constrains $SU(N)$5 and $SU(N)$6 by analyticity and unitarity. In $SU(N)$7, for example, one finds

$SU(N)$8

and boundary solutions exhibit a sharp antisymmetric resonance whose mass and width match the world-sheet axion inferred on the lattice (Miro et al., 2019).

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 $SU(N)$9 is

$0^-$0

with $0^-$1 and level matching $0^-$2 (Athenodorou et al., 2010). For many low-lying states in $0^-$3 $0^-$4 gauge theories, the spectrum follows this Nambu–Goto form down to $0^-$5–$0^-$6, even where the naive $0^-$7 expansion no longer converges (Athenodorou et al., 2010). The major exception is a family of $0^-$8 states lying anomalously low. A simple relativistic massive-mode ansatz,

$0^-$9

with $\bar{M}$0, qualitatively describes these deviations and suggests that some flux-tube channels probe internal massive structure rather than only transverse Goldstone motion (Athenodorou et al., 2010). This is closely aligned with the bootstrap picture in which boundary $\bar{M}$1-matrices develop resonance structure, including the world-sheet axion and broader states in other irreducible channels (Miro et al., 2019).

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

A more field-theoretic realization derives an effective $\bar{M}$2-dimensional theory inside a flux-tube environment. Starting from the QCD$\bar{M}$3 action and imposing longitudinal dominance, transverse confinement, and a factorized ansatz for the quark field,

$\bar{M}$4

one obtains a compactified QCD$\bar{M}$5 action with effective coupling

$\bar{M}$6

and effective quark mass

$\bar{M}$7

where $\bar{M}$8 is the transverse eigenvalue (Koshelkin et al., 2012). 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 $\bar{M}$9 (Koshelkin et al., 2012). 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

$\tau$00

or, in Bjorken coordinates,

$\tau$01

while the particle distributions satisfy a relativistic Boltzmann equation with Lorentz force, source, and collision terms,

$\tau$02

(Ruggieri et al., 2015). Pair production is implemented through the Schwinger source, for example in the static-box frame,

$\tau$03

with

$\tau$04

and the total current contains both matter and displacement components (Ruggieri et al., 2015). For small $\tau$05 the model yields $\tau$06 fm/c and $\tau$07 fm/c, while the initial color-electric field decays within $\tau$08 fm/c almost irrespective of viscosity; at large $\tau$09, persistent plasma oscillations modulate the late-time pressure anisotropy (Ruggieri et al., 2015). A related review emphasizes that for $\tau$10 the transverse-momentum slopes are close to equilibrium within $\tau$11 fm/c, and for $\tau$12 viscous hydrodynamics already describes the system well, indicating rapid hydrodynamization (Ryblewski, 2016).

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 (Ruggieri et al., 2015, Ryblewski, 2016).

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 $\tau$13 pure gauge theory and $\tau$14-flavor QCD, the longitudinal chromoelectric field in the mid-plane is well fitted by the dual Abrikosov/Clem profile up to about $\tau$15 fm in pure gauge and $\tau$16 fm with dynamical quarks (Cea et al., 2017). The extracted penetration length is approximately $\tau$17–$\tau$18 fm, the coherence length is of order $\tau$19–$\tau$20 fm in pure gauge and $\tau$21–$\tau$22 fm in $\tau$23-flavor QCD, the Ginzburg–Landau parameter satisfies $\tau$24 throughout the measured range, and the transverse width is approximately $\tau$25–$\tau$26 fm (Cea et al., 2017). 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 (Cea et al., 2017).

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 $\tau$27 versus $\tau$28 become non-linear, with effective slopes reduced relative to the massless-end value

$\tau$29

and with curvature increasing for heavier endpoint configurations (Nandan et al., 2017). In heavy-ion phenomenology, flux-tube networks with $\tau$30-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 (Patel, 2011).

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 $\tau$31 or $\tau$32 (Jakhad et al., 18 Sep 2025, Chen et al., 2019). 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 $\tau$33, $\tau$34, $\tau$35, and some $\tau$36 excitations (Jakhad et al., 29 Sep 2025, Jakhad et al., 18 Sep 2025). In baryon applications, the quark–diquark approximation neglects internal diquark excitations and possible three-body flux configurations (Jakhad et al., 2023, Jakhad et al., 2024). In heavy-ion applications, the field is abelianized, chromomagnetic components are neglected, and transport is imposed at fixed $\tau$37 rather than derived from microscopic QCD (Ruggieri et al., 2015). World-sheet EFT and lattice descriptions, conversely, are rigorous for long confining strings but are not by themselves constituent models of finite hadrons (Miro et al., 2019, Athenodorou et al., 2010).

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 (Jakhad et al., 18 Sep 2025, Miro et al., 2019, Ruggieri et al., 2015).