⚛️ nuclear theory

QCD thermodynamics through the crossover as a gas of confining strings with repulsive interactions

This paper demonstrates that QCD thermodynamics in the intermediate-temperature regime can be effectively described by a gas of confining strings with repulsive excluded-volume interactions, where an open-string spectrum characterized by a Hagedorn temperature of approximately 300 MeV and distinct mesonic and baryonic volume parameters successfully reproduces bulk lattice-QCD data above the pseudocritical temperature.

Original authors: Yuki Fujimoto, Volodymyr Vovchenko

Published 2026-09-22

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Yuki Fujimoto, Volodymyr Vovchenko

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

To understand the story of how matter behaves under extreme heat, one must first look at the building blocks of the universe. At the heart of this story are protons and neutrons, the particles that make up the atomic nucleus. Under normal conditions, these particles are held together by a powerful force that acts like a cosmic glue, binding smaller particles called quarks into tight, stable groups. However, if you heat this matter to temperatures far hotter than the center of the sun, this glue begins to loosen. The protons and neutrons melt apart, releasing their quarks to roam freely in a hot, dense soup known as a quark-gluon plasma. Scientists have long been fascinated by the moment this transition happens. It is not a sudden explosion, but a smooth crossover, a gradual shift where the rules of the game change. The big question is: what does the matter look like in the middle of this shift? Is it a chaotic mess of free particles, or does it retain some of its old, structured nature?

A team of researchers has taken a fresh look at this intermediate zone, the temperature range just above where protons and neutrons dissolve. They proposed a model that treats this hot matter not as a gas of free particles, but as a gas of tiny, vibrating strings. In this view, the particles we usually think of as solid objects are actually the ends of these strings. When the temperature rises, these strings stretch and vibrate, creating a vast spectrum of possible states. The researchers added a crucial twist to this idea: they accounted for the fact that these strings cannot occupy the same space at the same time. Just as people in a crowded room need personal space to move, these strings repel each other, preventing them from packing too tightly. By combining this string-like behavior with a simple rule about excluded space, the team was able to create a mathematical description that matches the data from supercomputer simulations of the early universe.

The researchers focused on a specific temperature range, starting just above the point where protons and neutrons break apart, which occurs at about 156 MeV. In their model, they replaced the traditional list of known particles with a continuous spectrum of string states. This spectrum grows exponentially as the mass increases, a feature that suggests there is a limiting temperature beyond which the system cannot be heated without changing its fundamental nature. They tested different values for this limiting temperature, which they call the Hagedorn temperature. While earlier studies of empty space suggested this limit might be around 340 MeV, the researchers found that a value of 300 MeV provided the best match for the behavior of matter in the hot, dense environment. At this temperature, the model successfully reproduced the pressure and energy density observed in the most advanced computer simulations of quantum physics.

A key part of their discovery involved how they handled the repulsion between these strings. In the simplest version of their model, they assumed that all strings, whether they formed mesons or baryons, pushed away from each other with the same strength. This worked well for describing the overall pressure of the system. However, when they looked closer at the fluctuations of electric charge and other conserved quantities, a discrepancy appeared. The model with equal repulsion overestimated how much the baryons—particles like protons and neutrons—were interacting. To fix this, the researchers allowed the baryons to have a larger "excluded volume" than the mesons. In plain terms, this means that in their model, baryons act as if they take up more space and push away from their neighbors more strongly than the lighter meson particles do. With this adjustment, the model's predictions for charge fluctuations aligned much better with the simulation data, particularly for the baryon sector.

The results suggest that even after the protons and neutrons have melted, the matter does not immediately become a completely free gas of quarks. Instead, it retains a structure that can be described as a gas of confining strings, where the strings still feel the pull of the force that usually binds quarks together. The study indicates that this string-like description remains valid and useful well above the transition point, up to temperatures of about 230 MeV. The researchers found that the repulsive interactions between these strings are essential to getting the physics right; without them, the model fails to capture the subtle details of how the matter behaves. The finding that baryons require a stronger effective repulsion than mesons hints that the internal structure of these heavier particles plays a unique role in the hot soup of the early universe.

While the model fits the data remarkably well, the authors are careful to note that it is a description of effective behavior rather than a final proof of what is happening at the microscopic level. The "strings" in their model are a mathematical tool to represent the complex interactions, and the parameters they adjusted, such as the size of the excluded volume, are effective values that capture the net result of many different forces. They acknowledge that there are still uncertainties, particularly regarding how these strings behave at even higher temperatures or in denser environments. Nevertheless, the work provides a compelling picture of a transitional phase where the universe is neither fully hadronic nor fully partonic, but a unique state where string-like correlations persist. This approach offers a new way to understand the thermodynamics of the early universe, bridging the gap between the world of stable particles and the world of free quarks.

Problem Statement
The paper addresses the challenge of describing Quantum Chromodynamics (QCD) thermodynamics in the intermediate-temperature regime, specifically above the chiral pseudo-critical temperature ( $T_c \simeq 156$ MeV). While the Hadron Resonance Gas (HRG) model successfully describes thermodynamics below $T_c$ , the transition to partonic degrees of freedom at higher temperatures remains nontrivial. Recent lattice studies suggest that confinement-related correlations persist above $T_c$ , motivating the investigation of whether a thermodynamic description based on confining strings can be extended into this regime. The specific problem is to determine if a gas of confining strings, supplemented by repulsive interactions, can reproduce bulk lattice-QCD thermodynamics and conserved-charge fluctuations without invoking a full partonic description.

Methodology
The authors construct a model combining a string-based description of the hadronic spectrum with van der Waals-type excluded-volume (EV) corrections to account for residual repulsive interactions.

String Density of States: The discrete hadron resonance spectrum is replaced by an exponentially growing open-string spectrum. The asymptotic density of states for open bosonic strings in 3+1 dimensions is given by $\rho_{str}(m) \sim m^{-3/2} e^{m/T_H}$ , where $T_H$ is the Hagedorn temperature related to the confining string tension $\sigma$ by $T_H = \sqrt{3\sigma/2\pi}$ . The model retains light pseudo-Nambu-Goldstone mesons ( $\pi, K, \eta$ ) explicitly, while higher-mass excitations are represented by a continuum string density above channel-dependent thresholds. Baryons are modeled as strings ending on a quark and a diquark.
Excluded-Volume Corrections: To address thermodynamic singularities and short-range repulsion, the authors implement EV corrections. They consider two scenarios:
- A common EV parameter ( $b$ ) applied uniformly to all mesonic and baryonic channels.
- A multi-component non-diagonal EV prescription allowing distinct parameters for mesons ( $b_M$ ) and baryons ( $b_B$ ). This uses coefficients derived from hard-sphere mixture virial expansions to define the available volume.
Calibration: The model parameters are calibrated against lattice-QCD data at vanishing baryon chemical potential. The common EV parameter is fitted to the pressure data ( $P/T^4$ ) in the range $110 \le T \le 280$ MeV. In the multi-component scenario, $b_M$ and $b_B$ are simultaneously fitted to bulk thermodynamic observables ( $P, \epsilon, s, c_s^2$ ) and conserved-charge susceptibilities ( $\chi_2^B, \chi_2^S, \chi_{11}^{QS}$ ).

Key Contributions and Results

Bulk Thermodynamics: The model provides a good description of bulk lattice-QCD thermodynamics (pressure, entropy density, speed of sound, and trace anomaly) up to $T \simeq 200\text{--}230$ $T ≃ 200 - 230$ MeV.
- The best agreement with lattice data is achieved with a Hagedorn temperature of $T_H \simeq 300$ MeV. This value is lower than the $T_H \simeq 320\text{--}340$ MeV inferred from fits to the vacuum hadron spectrum but aligns with the lower end of zero-temperature string-tension estimates.
- With a common EV parameter, the fit yields $b \simeq 0.33$ fm $^3$ (corresponding to an effective radius $r \simeq 0.27$ fm). This correction successfully reproduces the peak structure in the trace anomaly and the "softest point" in the speed of sound, features that the ideal string gas fails to capture.
Conserved-Charge Fluctuations:
- Using a common EV parameter, the model overestimates baryon-sensitive observables ( $\chi_2^B$ and $\chi_{11}^{BQ}$ ), indicating insufficient suppression in the baryon sector.
- Introducing separate EV parameters improves the description of these fluctuations. The simultaneous calibration yields $b_M \simeq 0.19$ fm $^3$ and $b_B \simeq 0.73$ fm $^3$ .
- This hierarchy ( $b_B > b_M$ ) suggests that stronger effective suppression is required in the baryon sector to match lattice data, consistent with recent observations that ideal string-gas descriptions underestimate repulsion in the light-baryon sector.
- Despite this improvement, noticeable tensions remain in certain observables, particularly the mixed susceptibility $\chi_{11}^{BQ}$ , indicating that a two-parameter EV prescription does not provide a complete description of charge-sensitive observables.

Significance and Claims
The paper claims that an open-string description supplemented by excluded-volume corrections offers a compact and effective model for lattice QCD thermodynamics in the intermediate-temperature regime. The results support the physical relevance of string-like degrees of freedom above the pseudocritical temperature, suggesting an intermediate region where effective string-like and quark-gluon descriptions may capture complementary aspects of the same thermodynamics.

The authors explicitly state that their findings do not prove the system remains confined, but rather that string-like correlations remain thermodynamically relevant. The fitted hierarchy of EV parameters ( $b_B > b_M$ ) demonstrates that the baryonic contribution requires stronger effective suppression within the model, though this does not necessarily establish a stronger microscopic baryon-baryon repulsive force. The authors acknowledge systematic uncertainties, including the use of asymptotic string density near thresholds and the assumption of mass-independent EV parameters, noting that further systematic studies are required for a quantitative determination of effective string interactions. The model is intended for the intermediate regime and is not expected to apply to asymptotically high temperatures where weak-coupling partonic methods are required.

More like this