Holographic study of higher-order baryon number susceptibilities at finite temperature and density

The cumulants of baryon number fluctuations serve as a good probe for experimentally exploring the QCD phase diagram at finite density, giving rise to characteristic fluctuation patterns associated with a possible critical endpoint (CEP). We compute the higher-order baryon number susceptibilities at finite temperature and baryon chemical potential using a holographic QCD model to address the non-perturbative aspect of strongly coupled QCD matter. The model can accurately confront lattice QCD data on a quantitative level and the location of the CEP is found to fall within the range accessible to upcoming experimental measurements. The baryon number susceptibilities up to the twelfth order are computed, and the collision energy dependence of different ratios of these susceptibilities is examined along the chemical freeze-out line. The holographic results show quantitative agreement with experimental data and the functional renormalization group results in a large collision energy range, with all ratios exhibiting a peak structure around 5-10 GeV. The mismatching between our holographic results with experimental data for sufficiently low collision energy is possibly due to non-equilibrium effects and complex experimental environments. The future experiments with measurements in the low collision energy range $\sqrt{S_{NN}}\approx 1-10~\text{GeV}$ and reduced experimental uncertainty could reveal more non-monotonic behavior signals which can be used to locate the CEP.