The QCD scalar susceptibility and thermal scalar resonances in chiral symmetry restoration

Building upon recent results on the role of thermal resonances in chiral symmetry restoration, we show that a description of the QCD scalar susceptibility at finite temperature $T$ saturated by the thermal properties of the lightest scalar resonance, the $f_0(500)$, is compatible both with lattice QCD data at nonzero $T$ and with the $T=0$ light resonance properties coming from experimental data. The thermal $f_0(500)$ is generated within the framework of Unitarized Chiral Perturbation Theory. This method allows us to achieve a good description of lattice QCD results with a reliable pion mass dependence. In particular, we perform direct fits to the chiral susceptibility measured in lattice data at different pion masses and temperatures, obtaining a remarkable agreement for the susceptibility and for mass differences of the light quark condensate. In addition, the fitted low-energy constants are compatible with $T=0$ phenomenology. Our results confirm the role of unitarized approaches and thermal resonances in the dynamics of the QCD transition.

💡 Research Summary

In this work the authors investigate the QCD scalar susceptibility χ S(T), a key observable that quantifies fluctuations of the light‑quark condensate and signals the restoration of chiral symmetry at finite temperature. Starting from the definition χ S(T)=−∂⟨ \bar qq ⟩_l/∂m_l, they note that χ S is directly related to the two‑point correlator of the scalar operator σ_l∼\bar q_l q_l, whose quantum numbers coincide with those of the lightest scalar resonance f₀(500). The central hypothesis is that the temperature dependence of χ S can be saturated by the thermal properties of the f₀(500) pole.

The authors first review the limitations of the Linear Sigma Model (LSM), which, although it provides a simple analytic link between χ S and the σ self‑energy Σ(k=0;T), fails to generate the characteristic peak of χ S around the crossover temperature T_c. To overcome this, they employ Unitarized Chiral Perturbation Theory (UChPT). In UChPT the ππ scattering amplitude is unitarized via the Inverse Amplitude Method, generating a pole on the second Riemann sheet that moves with temperature: s_p(T)=M_p²(T)−i M_p(T)Γ_p(T). The real part of this pole, M_S²(T)=Re s_p(T), is identified with the thermal “sigma” mass that enters the susceptibility. The saturated form reads

 χ_U(T)=A B_phys² M_S²(0)/M_S²(T),

where B≡M_π²/(2 m_l) and A=χ_U(0)/B_phys² is a normalization constant to be fitted.

To connect with lattice QCD data, the authors use the subtracted condensate Δ_{l,s}=⟨\bar qq⟩l−(2 m_l/m_s)⟨\bar ss⟩ and the derived dimensionless order parameter M=−m_s^phys Δ{l,s}/F_K^4. Its derivative with respect to the light quark mass defines the “chiral susceptibility” χ_M, which, after a short algebraic manipulation, can be expressed as a linear combination of χ S, the strange condensate, and mixed derivatives. Using SU(3) ChPT at next‑to‑leading order (NLO), the authors show that the contributions from strange quarks are Boltzmann‑suppressed at temperatures near T_c, while the pion loops dominate. Consequently, χ_M≈χ S up to higher‑order corrections, a relation that is verified both analytically and numerically.

Armed with this approximation, the authors construct a fitting function that combines the saturated UChPT piece with the NLO ChPT thermal correction:

 χ_UM(T)=A(M_π) M_S²(0)/M_S²(T) + (1/4) B_phys² Δχ_ChPT^M(T).

Here Δχ_ChPT^M(T) contains the well‑known pion thermal function g₂(M_π,T) and small kaon/eta contributions. The low‑energy constants (LECs) l₁ʳ and l₂ʳ, which dominate the ππ scattering vertices and thus the pole position, are treated as fit parameters, while l₃ʳ and l₄ʳ are fixed to their standard ChPT values.

The numerical analysis uses lattice data from the HotQCD collaboration (N_τ=8) for three pion masses: M_π≈110, 140, and 160 MeV. For each mass the authors fit the normalization A(M_π) and, in separate fits, also l₁ʳ, l₂ʳ. The resulting A values are compatible with the theoretical expectation A_ChPT, confirming the internal consistency of the approach. The fitted LECs lie within the phenomenological ranges obtained from zero‑temperature ππ scattering analyses, demonstrating that the finite‑temperature fit does not spoil the low‑energy description.

The peak positions of the susceptibility are reproduced with impressive accuracy for the physical and slightly heavier pion masses: T_c≈163.8±5.8 MeV for M_π=140 MeV and T_c≈160.9±5.7 MeV for M_π=160 MeV, both in line with the lattice crossover temperature (≈155 MeV). For the lighter mass M_π=110 MeV the model predicts a much higher and sharper peak, reflecting the fact that the present N_f=2 framework does not fully capture the lattice procedure where the strange quark mass is held fixed while the light mass is varied. The authors acknowledge this limitation and suggest that the method is reliable down to pion masses of about 130 MeV.

Overall, the study provides strong quantitative evidence that the thermal f₀(500) resonance dominates the scalar susceptibility near the QCD crossover. The unitarized approach improves upon the Hadron Resonance Gas model, which lacks a peak structure, and offers a theoretically sound bridge between lattice results, experimental resonance properties, and chiral effective theory. The work also highlights the sensitivity of χ S to the LEC combination l₂ʳ+1.5 l₁ʳ, suggesting that future high‑precision lattice measurements of χ S could be used to refine these constants.

In the concluding section the authors outline possible extensions: incorporating the I=1/2 scalar channel (K*₀(700)), exploring the impact of an axial chemical potential μ₅, and studying the temperature dependence of the critical line T_c(μ₅). Such studies would further elucidate the interplay between unitarized resonances, Ward identities, and the dynamics of chiral symmetry restoration in QCD.