1　Introduction

According to the Big Bang theory, under the extremely high temperature and high energy density in the early stage of the universe, the quark–gluon plasma (QGP), a new form of matter, was generated by the release of quarks and gluons that had been bound in hadrons by the strong interaction. The masses of the heavy-flavor quarks, m_c ˜ 1.3 GeV/c² and m_b ˜ 4.8 GeV/c², are larger than those of light quarks and the quantum chromodynamics (QCD) energy scale (Λ_QCD). Therefore, the generation of heavy-flavor quarks requires sufficiently large energy and momentum transfer by initial hard scattering in heavy-ion collisions and can be calculated using perturbation QCD [1, 2]. Heavy-flavor quarks have a high probability of undergoing relatively complete evolution in QCD materials; thus, they are ideal probes to study the properties of the QGP in heavy-ion collisions. In particular, it is crucial to study the interaction between heavy-flavor quarks and the medium by measuring the nuclear modification factor, flow, and production yield of charmed mesons [3-12].

The transverse momentum (p_T) spectra of identified particles provide valuable data for discoveries based on relativistic heavy-ion collisions [13, 14]. The Boltzmann–Gibbs blast-wave model, which has a compact set of parameters, namely, temperature (T) and flow velocity (β), has shown that the spectral shape is sensitive to the dynamics of nucleus-nucleus (AA) collisions [13, 15] and can be used to describe the transverse momentum distributions of light-flavor hadrons with different masses. Moreover, the Tsallis blast-wave model, which includes the fluctuation of the initial conditions in the hydrodynamic evolution on an event-by-event basis, has been used to study the π^±, K^±, p($\overline{\text{p}}$), ϕ, $\Lambda (\overline{\Lambda })$, and ${\Xi }^{-}({\overline{\Xi }}^{+})$ spectra in Au-Au collisions at 200 GeV [16]. It was found that the average transverse flow velocity 〈β〉 and temperature at kinetic freeze-out T_fro increased with centrality, whereas the non-extensive parameter q showed the opposite behavior. The thermodynamic formulation of non-extensive statistics was recently applied to the data for all light hadrons, and a strong grouping phenomenon in the T–(q-1) parameter space was considered to be a probe for the study of QGP in small and large systems. The non-extensive temperature converges to T≈ 0.144 ± 0.010 GeV [17], which is consistent with the QCD results [18, 19], indicating that the collisional system may have undergone QGP evolution. However, the non-extensive parameter q ≠ 1, indicating that the conclusion is based entirely on the non-extensive properties. Although the implications and consequences of such an application to particle production are still under investigation, the function is relatively easy to understand [20-22]. Heavy flavors are expected to behave differently from light flavors and thus may have different non-extensive thermal parameters. The same concept with a non-extensive statistical distribution function is extended to fit the charmed meson spectra recently measured in the ALICE and STAR experiments.

We present the Tsallis–Pareto distribution formula and explain the physical origin of the relationship between the parameters T and q in Sect. 2 of this article. The analysis procedure and fitting results for charmed mesons and the thermal temperature after transverse flow correction are also described in this section. In Sect. 3, the results for the 0–10% (central) and 30–50% (semiperipheral) centrality bins in Pb-Pb collisions at 2.76 TeV and for the 0–10%, 30–50%, and 60–80% (peripheral) centrality bins in Pb-Pb collisions at 5.02 TeV are presented. In the same section, the results obtained in Pb-Pb collisions are compared with the values of T and q measured in proton-(anti)proton ［pp($\overline{\text{p}}$)］ collisions over a wide energy range. The results and conclusions are summarized in Sect. 4.

2　Tsallis–Pareto distribution and its fit to charmed meson spectra

Much work on high-energy particle collisions has focused on the study of the transverse momentum distributions of outgoing particles. In the low-pT regime of the spectra, the conventional exponential distributions can be used to describe the spectral shape, and the formula, assuming vanishing chemical potential at high energies, is given as

where m_T = $\sqrt{p_{\text{T}}^2+m^2}$ is the transverse mass, m is the rest mass of the hadron species, and T is the corresponding temperature. Instead of Eq.(1), a power-law distribution has been used in high-energy physics to better characterize the spectra with $p_{\text{T}}\underset{˜}{>}$3 GeV/c. The Tsallis–Pareto distribution derived from non-extensive statistics, which is a generalization of the conventional Boltzmann–Gibbs theory, is thought to describe the full p_T spectrum of hadrons. The transverse momentum dependence of identified hadron production from recently published data has been successfully described by non-extensive statistics [23]. In addition, for the theoretical limit $q\to 1$, exponential and logarithmic functions can be obtained [24-26]. The exploration of this new research direction in statistical mechanics began when Tsallis proposed a non-extensive entropy in 1988 [24]. Non-extensive entropy (the Tsallis entropy) is a generalization of the conventional Boltzmann–Gibbs entropy, and it formed the basis of non-extensive statistics.

The Tsallis–Pareto distribution within the non-extensive statistics in [17] is also used in the current analyses:

where A is a normalizing factor and can reflect the production yield of the hadron p_T spectrum. Tq is the temperature in the non-extensive statistical theory; the subscript q is omitted for brevity in the following discussion. Note that Tq can differ from the temperature T in Eq.(1), but its physical meaning should be the same in the limiting case $q\to 1$.

The correlations between the parameters T and q have been presented in earlier studies [27-30]. Furthermore, the charged particle multiplicity can be derived from the Tsallis–Pareto-distributed transverse momentum, and concrete application to experimental data yields a negative binomial distribution parameter k ˜ $\mathcal{O}$(10) [31-38]. In addition, many studies have revealed the importance of measuring the event-by-event multiplicity and its fluctuation. They found that the yield of strange hadrons is positively related to the multiplicity, and a long-range correlation can be observed in small collisional systems as the multiplicity increases [39-43]. This study focuses on the relationship between the parameters of heavy-flavor hadrons in the T–(q-1) parameter space using the non-extensive statistics described above. For simplicity, the fluctuations in the number of produced particles can be explained in a one-dimensional relativistic gas model [44], and the Tsallis parameters under consideration are given as

where M is the number of particles at energy E.

In the thermodynamic picture, the relationship between T and q can be obtained from Eqs.(3) and (4), assuming that the relative size of the multiplicity fluctuations is constant as in [17]:

Thus, we have the following formula in the current analyses:

This formula is used to measure the relationship between the Tsallis parameters and event multiplicity in charmed meson production for both small and large systems over a wide range of collision energies and hadron transverse momenta, and the results are compared with the corresponding results for light hadrons.

2.1　 Application of Tsallis–Pareto distribution to charmed meson spectra

In this study, we analyze the transverse momentum dependence of charmed meson production in pp($\overline{\text{p}}$) and AA collisions measured by the ALICE, CDF, and STAR collaborations [45-52]. We find that the charmed meson data from small and large collision systems with broad selection criteria can be investigated simultaneously. D⁰ and D^* spectra were measured by the STAR Collaboration in pp collisions at 200 GeV for 0.4 < p_T < 8.8 GeV/c and at 500 GeV for 1.5 < p_T < 17.7 GeV/c. D⁰ spectra were measured by the CDF Collaboration in p $\overline{\text{p}}$ collisions at 1.96 TeV for 5.8 < p_T < 16.0 GeV/c and by the ALICE Collaboration in pp collisions at 7 TeV for 1.5 < p_T < 13.8 GeV/c. Figure 1 shows these pT spectra of D mesons, together with our fit results at four energies (200 GeV, 500 GeV, 1.96 TeV, and 7 TeV) in pp($\overline{\text{p}}$) collisions. The fit parameters and χ²/ndf are listed in Table 1 . The solid curves from the Tsallis–Pareto distributions describe the data well. Error bars denote the quadratic sums of the statistical and systematic uncertainties. Data are scaled by factors of 10ⁿ for better visibility.

Fig. 1Full-screen View

Transverse momentum distributions dN/dp_T of D⁰/D^* in pp($\overline{\text{p}}$) collisions at 200 GeV, 500 GeV, 1.96 TeV, and 7 TeV, from bottom to top. Solid curves are results of Tsallis–Pareto fit. Error bars are quadratic sums of statistical and systematic uncertainties, and data are scaled by factors of 10ⁿ for better visibility.

Table 1Full-screen View

Values of parameters from Tsallis–Pareto fit to charmed meson spectra in Pb-Pb ［pp($\overline{\text{p}}$)］ collisions. The uncertainties are from the fit.

$\sqrt{s_{\text{NN}}}$(GeV)	Centrality	Charmed meson	A	T	q	χ2/ndf
Pb-Pb, 2760	0–10%	D0	1.279± 0.516	0.239±0.030	1.166± 0.014	5.15/6
		D+	1.052± 0.350	0.201± 0.024	1.180± 0.013	9.13/7
		D*+	0.515± 0.212	0.240±0.031	1.179± 0.012	2.9/7
	30–50%	D0	0.178± 0.076	0.278±0.041	1.169±0.020	0.88/5
		D+	0.052± 0.025	0.322± 0.055	1.151±0.025	0.51/5
		D*+	0.084± 0.045	0.250±0.049	1.205± 0.025	0.34/5
Pb-Pb, 5020	0–10%	D0	1.957± 0.551	0.240± 0.020	1.187±0.008	14.95/9
		D+	0.768±0.277	0.245± 0.025	1.190± 0.008	5.16/9
		D*+	0.675± 0.209	0.258± 0.026	1.184±0.010	6.12/9
	30–50%	D0	0.145± 0.040	0.328± 0.026	1.175± 0.008	2.12/9
		D+	0.088± 0.024	0.311± 0.024	1.179± 0.007	5.35/9
		D*+	0.057± 0.020	0.331±0.034	1.185± 0.008	3.76/9
	60–80%	D0	0.010± 0.003	0.427± 0.042	1.151± 0.012	2.52/8
		D+	0.005± 0.001	0.402± 0.043	1.170± 0.013	2.21/8
		D*+	0.004± 0.002	0.430± 0.069	1.156± 0.017	4.48/8
pp, 200		D0+D*	3.8e-4 ± 5.3e-5	0.322± 0.022	1.081± 0.011	4.39/7
pp, 500		D0+D*	6.8e-4 ± 1.7e-4	0.310± 0.020	1.132± 0.006	4.13/10
p($\overline{\text{p}}$), 1960		D0	1.8e-3±1.7e-3	0.386± 0.058	1.143± 0.011	0.87/3
pp, 7000		D0	1.6e-3 ± 7.6e-4	0.494± 0.062	1.139± 0.023	1.75/6

Show more

In Fig. 2, the fits are applied to the transverse momentum distributions of charmed mesons for different centrality bins in AA collisions [49-52]. The fits of the transverse momentum distributions of prompt D⁰, D⁺, and D^*+ mesons in Pb-Pb collisions at 2.76 TeV are shown in Fig. 2(a). Solid circles, diamonds, and triangles represent D⁰, D⁺, and D^*+, respectively. The 0–10% (solid) and 30–50% (open) centrality bins are shown in Fig. 2(a) [49], where the D⁰ and D^*+ production yields at 0–10% and 30–50% are scaled by factors of 10 and 0.05, respectively. For the 0–10% centrality bin, fits are performed in the range 0 < p_T < 20.0 GeV/c for D⁰ and in the range 0 < p_T < 30.0 GeV/c for D⁺ and D^*+. For the 30–50% centrality bin, the same fitting procedures were performed in the range 0 < p_T < 14.0 GeV/c for D⁰, D⁺, and D^*+. The vertical bars represent the sums of the statistical and systematic uncertainties. Table 1 shows that the fit parameters of the charmed meson spectra in the 0–10% centrality bin have smaller uncertainties with more data points than those in the 30–50% centrality bin at 2.76 TeV.

Fig. 2Full-screen View

Transverse momentum distributions dN/dp_T of D⁰, D⁺, and D^*+ for different centrality bins in Pb-Pb collisions at 2.76 TeV (a) and 5.02 TeV (b), (c), (d) , where the production yields are scaled by various factors for visibility. Vertical bars represent quadratic sums of statistical and systematic errors; symbols are placed at the center of the bin. Detailed descriptions are presented in Sect. 2.1.

Figure 2(b), (c), and (d) show the transverse momentum distributions dN/dp_T of D⁰ mesons (solid circles), D⁺ mesons (diamonds), and D^*+ mesons (triangles) in the 0–10%, 30–50%, and 60–80% centrality bins, respectively, in Pb-Pb collisions at 5.02 TeV [50]. The vertical bars represent the quadratic sums of the statistical and systematic uncertainties, and symbols are placed at the center of the bin. The solid curves representing the Tsallis–Pareto distributions describe the data well. For visibility, the D⁰ and D^*+ distributions in the three centrality bins are scaled by factors of 10 and 1/10, respectively. To more physically constrain the D⁺ and D^*+ yields at p_T = 1.5 GeV/c, we applied a D⁺/D⁰ and D^*+/D⁰ ratio of approximately 0.5 from PYTHIA and performed the fit. The ratio obtained from PYTHIA is consistent with the experimental data [50]. The fitted T, q, and A parameters and χ²/ndf values are listed in Table 1 . The T and q parameters after transverse flow correction are shown in Fig. 5. In addition, we also applied the same Tsallis–Pareto fits to the transverse momentum spectra of π^±, K^±, and p($\overline{\text{p}}$) in 0–5%, 5–10%, 10–20%, 20–30%, 30–40%, 40–50%, 50–60%, 60–70%, 70–80%, and 80–90% centrality bins in Pb-Pb collisions at 2.76 and 5.02 TeV [51, 52]. The fit parameters after transverse flow correction are also shown in Fig. 5. we studied ϕ, Λ⁰, and Ξ for different centralities at 2.76 TeV using the same method. The point-by-point statistical and systematic uncertainties were added as a quadratic sum when we performed these fits.

Fig. 5Full-screen View

Thermal temperature T versus q-1 and values of parameters from the Tsallis–Pareto fit of the identified particle spectra at different centralities in Pb-Pb ［pp($\overline{\text{p}}$)］ collisions at 2.76 and 5.02 TeV (200 GeV, 500 GeV, 1.96 TeV, and 7 TeV) after transverse flow correction. Shaded vertical band marks the saturated value of q-1 = 0.142 ± 0.010 in pp($\overline{\text{p}}$) collisions with increasing energy. The solid and dotted lines are from Eq.(6); the parameters are listed in Table 1.

2.2　 Thermal temperature with flow correction

The phenomenological model can describe almost all hadronic spectra by beginning with thermalization and collective flow as basic assumptions [15]. The mass dependence of the effective temperature T has been described by introducing a Gaussian parameterization [53-56] and can be interpreted as the presence of a radial flow. The velocity of the radial flow, which is generated by violent nucleon-nucleon collisions in two colliding nuclei and developed in both the QGP phase and hadronic rescattering, increases the transverse momentum of particles in proportion to their mass [15, 57, 58]. Many models have been used to investigate the radial flow [59]; a radial flow model [15, 58] we can use in this analysis is written as

$T=T_{\text{fro}}+m{\langle u_{\text{t}}\rangle }^2\mathrm{,}$ (7)

where T_fro is the hadron kinetic freeze-out temperature, and 〈u_t〉 is a measure of the strength of the (average radial) transverse flow. The relationship between the average transverse velocity 〈β_t〉 and 〈u_t〉 is given as

$\langle {\beta }_{\text{t}}\rangle =\frac{\langle u_{\text{t}}\rangle }{\sqrt{1+{\langle u_{\text{t}}\rangle }^2}}\mathrm{.}$ (8)

Note that although the T value according to the non-extensive statistical theory can differ from the usual temperature in Eq.(1), the flow correction of the spectral temperatures is independent of the statistical model. In addition, the following function is used to study the collectivity of charmed mesons produced in heavy-ion collisions [5]:

$\frac{{\text{d}}^{2}N}{2\pi m_{\text{T}}\text{d}{m}_{\text{T}}\text{d}y}=\frac{\text{d}N/\text{d}y}{2\pi T(m_0+T)}{e}^{-(m_{\text{T}}-m0)/T}\mathrm{,}$ (9)

where m₀ is the rest mass of the hadron species. This method is used to analyze the charmed meson spectra and to understand the collective velocity of the radial flow from the data in Pb-Pb collisions at 2.76 TeV. The function in Eq.(9) is applied to the π, K, p($\overline{p}$), ϕ, Λ⁰, Ξ, Ω, and D⁰ spectra for 0–10% and 30–50% centrality at 2.76 TeV. For Λ⁰, Ξ, and Ω, the semiperipheral results are obtained at 20–40% centrality. The fitted parameters are shown in Fig. 3. The solid and open circles represent the results for 0–10% and 30–50% centrality, respectively. The slope parameters obtained for D⁰ mesons and strange hadrons at 0–10% and 30–50% centrality are 0.113± 0.037 and 0.107± 0.036, respectively. In addition, the linear fits of the ϕ, Λ⁰, Ξ, Ω, and D⁰ data points in Fig. 3 show a smaller slope than those of the π, K, and p($\overline{\text{p}}$) data points, indicating that the former may freeze out earlier and gain less transverse collectivity during system evolution. All fits were performed up to m_T - m₀ values of < 1 GeV/c² for π, K, and p($\overline{p}$), <2 GeV/c² for ϕ, Λ⁰, Ξ, and Ω, and <3 GeV/c² for D⁰.

Fig. 3Full-screen View

Effective temperature as a function of hadron mass for 0–10% (solid circles) and 30–50% (open circles) centrality in Pb-Pb collisions at 2.76 TeV. The lines are fits from Eq.(7).

The fitted Tsallis–Pareto T parameters in Sect. 2.1 are plotted as a function of the corresponding masses of π^±, K^±, and p($\overline{p}$) in the 0–5%, 5–10%, 10–20%, 20–30%, 30–40%, 40–50%, 50–60%, 60–70%, 70–80%, and 80–90% centrality bins in Pb-Pb collisions at 2.76 TeV and 5.02 TeV. The 〈u_t〉 values from Eq.(7) were extracted, and the 〈β_t〉 distributions as a function of 〈dN_ch/dη〉 in Pb-Pb collisions at 2.76 and 5.02 TeV are indicated by open squares and open diamonds, respectively, in Fig. 4. The D⁰ meson is added to the above process at 0–10% and 30–50% centrality in Pb-Pb collisions at 2.76 TeV and is extrapolated to a lower centrality with the same centrality dependence as light-flavor hadrons. We plot the relationship between 〈β_t〉 and 〈dN_ch/dη〉 as solid circles in Fig. 4. We finally obtained the thermal temperatures after flow correction of the spectral temperatures for charmed mesons, as shown in Sect. 3.

Figure 4 reveals that the multiplicity dependence of the average radial flow velocity is linear as a function of log 〈dN_ch/dη〉. Furthermore, the average radial flow velocity increases with increasing multiplicity, which is consistent with an earlier blast-wave analysis [17]. In addition, the values can be extracted for π^±, K^±, and p($\overline{\text{p}}$) using the non-extensive statistical theory:

Fig. 4Full-screen View

Average radial flow velocity from Eq.(8) as a function of average event multiplicity. Solid circles represent ϕ, Λ⁰, Ξ, Ω, and D⁰, which are combined to extract 〈u_t〉 at 0–10% and 30–50% centrality. The results for light hadrons at 2.76 and 5.02 TeV are indicated by open squares and open diamonds, respectively. The fit result is shown at the bottom of the panel.

$\langle {\beta }_{\text{t}}\rangle =\left(0.286\pm 0.013\right)+\left(0.065\pm 0.006\right)\log \langle \text{d}{N}_{\text{ch}}/\text{d}\eta \rangle \mathrm{.}$ (10)

The linear dependence for charmed mesons is

$\langle {\beta }_{\text{t}}\rangle =\left(0.129\pm 0.037\right)+\left(0.065\pm 0.000\right)\log \langle \text{d}{N}_{\text{ch}}/\text{d}\eta \rangle \mathrm{.}$ (11)

3　Results and discussion

The measured transverse momentum distributions of charmed mesons in pp($\overline{\text{p}}$) collisions from the ALICE, CDF, and STAR collaborations are fitted with Tsallis–Pareto distributions in Fig. 1. The fit results are in good agreement with the data points overall. The same Tsallis–Pareto fits were also applied to the transverse momentum distributions of charmed mesons in Pb-Pb collisions at 2.76 and 5.02 TeV. Figure 2 shows the measured charmed meson spectra and the fitted Tsallis–Pareto curves for different centrality bins in Pb-Pb collisions at 2.76 TeV ［Fig. 2(a)］ and 5.02 TeV ［Fig. 2(b)-(d)］.

The non-extensive feature q ≠ 1 appears in Table 1 . To obtain a more physical understanding of this result, the dependence of the parameter q on the size of the collisional system is represented diagrammatically in Fig. 5. In addition, the spectral temperature T obtained from the Tsallis–Pareto fit of the hadron spectrum is larger than the original temperature T_thermal by a blue-shift factor owing to the presence of a radial flow [15]. The ordinate in Fig. 5 shows the thermal temperature values obtained by flow correction of the spectral temperatures for different hadron species. The flow correction formalism is given as

where β_t is given by Eq.(10) for π^±, K^±, and p($\overline{p}$) and by Eq.(11) for charmed mesons.

Figure 5 shows the T–(q-1) parameter space for π^±, K^±, p($\overline{p}$), and charmed mesons in Pb-Pb ［pp($\overline{\text{p}}$)］ collisions at 2.76 and 5.02 TeV (200 GeV, 500 GeV, 1.96 TeV, and 7 TeV). To simplify this figure, the results for strange particles (ϕ, Λ⁰, and Ξ) are not shown. The open and solid symbols represent the Tsallis–Pareto parameter values at 2.76 and 5.02 TeV, respectively. In addition, 〈dN_ch/dη〉 = 2047 and 〈dN_ch/dη〉 = 19.5 represent the central and peripheral centrality bins for light hadrons, respectively. There are five points around T = 0.25 at q - 1 = 0. This result indicates that at such high energies, the light-flavor particles are completely thermalized in most central Pb-Pb collisions; thus, their q values do not deviate from 1 at high temperature. The dotted and solid lines show the fitting results from Eq.(6). The fits for charmed mesons were performed in the range 0.14 < q-1 < 0.21 for 2.76 and 5.02 TeV, and the fit parameters are tabulated in Table 2 . The charmed meson results for pp($\overline{\text{p}}$) collisions at different energies are indicated by diamonds. The vertical shaded band marks the saturated value of q-1 = 0.142 ± 0.010 in pp($\overline{\text{p}}$) collisions with increasing energy. This value was obtained by constant fitting in the saturated region. To further investigate the T–(q-1) correlations between the pp($\overline{\text{p}}$) and Pb-Pb systems, we extended the fit lines along the decreasing q-1 direction for charmed mesons in Pb-Pb collisions.

We find that the Tsallis–Pareto distributions can provide satisfactory descriptions in a wide range of transverse momentum dependence of charmed meson production in pp($\overline{\text{p}}$) and AA collisions in a wide energy range. The values of the corresponding ζ² parameters in Table 2 can be used to test the strength of the correlation between T and q. This correlation is inversely proportional to the hadron mass overall. The following conclusions can be drawn from Table 2 and Fig. 5.

(i) There is a significant linear relationship between the thermal temperature and Tsallis q parameter for π^±, K^±, p($\overline{\text{p}}$), and charmed mesons in Pb-Pb collisions at 2.76 TeV. In addition, the slope of the T–(q-1) parameter is positively correlated with the hadron mass. The same conclusion can be drawn for Pb-Pb collisions at 5.02 TeV, where it is better supported because of the larger statistics and more complete centrality bins.

(ii) Charmed mesons have a significantly higher slope than light hadrons. The temperature of charmed mesons is found to be higher than that of light hadrons at the same q-1, indicating that heavy flavor requires a higher temperature to reach the same degree of non-extensivity as light flavors in heavy-ion collisions. The slope of a given hadron is smaller at 2.76 TeV than at 5.02 TeV. In Fig. 6, the slopes of the T–(q-1) correlations are plotted as a function of hadron mass for π^±, K^±, p($\overline{\text{p}}$), ϕ, Λ⁰, Ξ, and charmed mesons in Pb-Pb collisions at 2.76 and 5.02 TeV (open and solid circles, respectively). The data are fitted by a quadratic polynomial function, and a deeper theoretical explanation is needed.

Fig. 6Full-screen View

Slope of T–(q-1) correlations as a function of hadron mass in Pb-Pb collisions at 2.76 TeV (open circles) and 5.02 TeV (solid circles). The curves are quadratic polynomial fits.

(iii) The charmed meson results for pp($\overline{\text{p}}$) collisions at different energies show that the thermal temperature increases with system energy, whereas the q parameter becomes saturated at the pp($\overline{\text{p}}$) limit, q-1 = 0.142 ± 0.010. In addition, the results for most peripheral Pb-Pb collisions are found to approach the pp($\overline{\text{p}}$) limit, which suggests that more peripheral heavy-ion collisions are less affected by the medium and more similar to pp($\overline{\text{p}}$) collisions.

4　Summary

We presented fits of the transverse momentum spectra of D⁰, D⁺, and D^*+ mesons at mid-rapidity in Pb-Pb collisions at 2.76 and 5.02 TeV. A similar analysis with non-extensive statistics was performed to identify light hadron spectra for different centrality bins in Pb-Pb collisions at 2.76 TeV and 5.02 TeV after flow correction. Charmed meson production can be well described by the Tsallis–Pareto distributions. We observed that in the T–(q-1) parameter space, the slope has a positive dependence on hadron mass. In addition, the temperature of charmed mesons was found to be higher than that of light hadrons at the same q-1, indicating that heavy flavor requires a higher temperature to reach the same degree of non-extensivity as light flavors in heavy-ion collisions. In addition, the slope distribution of the T–(q-1) correlations (Fig. 6) and the anti-correlation between the thermal temperature and centrality for charmed mesons require a deeper theoretical explanation.

For the pp($\overline{\text{p}}$) collision system as a reference, we found that the thermal temperature increases with system energy, whereas the q parameter becomes saturated at the pp($\overline{\text{p}}$) limit, q-1 = 0.142 ± 0.010. Moreover, the results of most peripheral Pb-Pb collisions were found to approach the pp($\overline{\text{p}}$) limit, which suggests that more peripheral heavy-ion collisions are less affected by the medium and more similar to pp($\overline{\text{p}}$) collisions. In addition, uniform descriptions of both small and large systems over a wide range of collision energies and hadron transverse momenta were presented.