Next Article in Journal
Two Sides of the Same Coin: Sterile Neutrinos and Dark Radiation, Status and Perspectives
Next Article in Special Issue
CP Violation for the Heavens and the Earth
Previous Article in Journal
Impact of AGB Stars on the Chemical Evolution of Neutron-Capture Elements
Previous Article in Special Issue
Fluctuations of Initial State and Event-by-Event Pseudo-Rapidity Correlations in High Energy Nuclear Collisions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multiplicity Dependencies of Midrapidity Transverse Momentum Distributions of Identified Charged Particles in Proton-Proton Collisions at (s)1/2 = 7 TeV at the LHC

by
Khusniddin K. Olimov
1,*,
Fu-Hu Liu
2,*,
Kobil A. Musaev
1 and
Maratbek Z. Shodmonov
1
1
Physical-Technical Institute of Uzbekistan Academy of Sciences, Chingiz Aytmatov Str. 2b, Tashkent 100084, Uzbekistan
2
State Key Laboratory of Quantum Optics and Quantum Optics Devices, Collaborative Innovation Center of Extreme Optics, Institute of Theoretical Physics, Shanxi University, Taiyuan 030006, China
*
Authors to whom correspondence should be addressed.
Universe 2022, 8(3), 174; https://doi.org/10.3390/universe8030174
Submission received: 21 October 2021 / Revised: 2 February 2022 / Accepted: 7 March 2022 / Published: 10 March 2022
(This article belongs to the Special Issue Universe: Feature Papers–High Energy Nuclear and Particle Physics)

Abstract

:
Dependencies of midrapidity pt distributions of the charged pions and kaons, protons and antiprotons on charged-particle multiplicity density (<dNch/>) in inelastic proton-proton collisions at (s)1/2 = 7 TeV at the LHC, measured by ALICE Collaboration, are investigated. The simultaneous minimum χ2 fits with the Tsallis function with thermodynamical consistence and the Hagedorn function with included transverse flow have well-described the pt spectra of the particle species in the ten studied groups of charged-particle multiplicity density. The effective temperatures, T, of the Tsallis function with thermodynamical consistence have shown a steady rise with increasing the charged-particle multiplicity in proton-proton collisions at (s)1/2 = 7 TeV, in agreement with the similar result obtained recently in proton-proton collisions at (s)1/2 = 13 TeV at the LHC. The respective T versus <dNch/> dependence in proton-proton collisions at (s)1/2 = 7 TeV is reproduced quite well by the simple power function with the same value (≈ 1/3) of the exponent parameter as that extracted in proton-proton collisions at (s)1/2 = 13 TeV. The identical power dependence T ~ ε 1 / 3 between the initial energy density and effective temperature of the system has been observed in proton-proton collisions at (s)1/2 = 7 and 13 TeV. We have observed that the transverse radial flow emerges at <dNch/> ≈ 6 and then increases, becoming substantial at larger multiplicity events in proton-proton collisions at (s)1/2 = 7 TeV. We have estimated, analyzing T0 and β t versus <dNch/> dependencies, that the possible onset of deconfinement phase transition in proton-proton collisions at (s)1/2 = 7 TeV occurs at <dNch/> ≈ 6.1 ± 0.3, which is close to the corresponding recent estimate (<dNch/> ≈ 7.1 ± 0.2) in proton-proton collisions at (s)1/2 = 13 TeV. The corresponding critical energy densities for probable onset of deconfinement phase transition in proton-proton collisions at (s)1/2 = 7 and 13 TeV at the LHC have been estimated to be 0.67 ± 0.03 and 0.76 ± 0.02 GeV/fm3, respectively.

1. Introduction

The widespread use [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25] of the Tsallis distribution function in high-energy proton-proton collisions is explained by its very good parameterization of the experimental pt spectra of hadrons with just a few parameters: the first one—effective temperature (T), the second one—parameter of non-extensivity, q, which accounts for deviation of pt distribution from the Boltzmann–Gibbs exponential distribution, and the third parameter—the fitting constant, assumed to be proportional to the system volume. It is important to note that the Tsallis function has the advantage of being connected through the entropy to thermodynamics, which is not true with other power law distributions [8]. The q and T parameters can also be employed for the identification of the system size scaling as well as initial conditions [25].
There are various modifications of the Tsallis function, which have equally well-described the pt distributions of final hadrons in proton-proton collisions up to the largest available pt values at the Relativistic Heavy Ion Collider (RHIC) and Large Hadron Collider (LHC) experiments [3,4,5,6,7,8,9,10,15]. The non-extensivity parameter, q, of the Tsallis function has shown quite noticeable sensitivity to the large pt region (pt > 3 GeV/c) of the invariant pt distributions of hadrons, suggesting the necessity of analyzing the longer pt intervals for extracting the more correct q values [22,23,24].
The different transverse flow models have been incorporated into Tsallis statistics to describe the pt distributions of hadrons in high-energy heavy-ion and proton-proton collisions at the RHIC and LHC. In most cases, the Blast-Wave model with Boltzmann–Gibbs statistics (the BGBW model) [26,27,28], the Blast-Wave model with Tsallis statistics (the TBW model) [29,30], the Tsallis distribution with transverse flow effect—improved Tsallis distribution [30,31,32], and the Hagedorn formula (function) with transverse flow [11,16,21,33] have been used to estimate the kinetic freeze-out temperature and transverse expansion velocity.
Most of the works analyzing high-energy collisions with various versions of the Tsallis function and its modified forms, incorporating other model functions, have used separate fits to pt distribution of particular particle species in a studied collision. However, as indicated in [34], it is not feasible to assign the physical meaning to collective properties, such as kinetic freeze-out temperature or radial transverse flow velocity, extracted from separate model fits to each particle type. The simultaneous model fits to pt distributions of various particle species in a studied collision system, performed by keeping the temperature and transverse flow velocity as the common (shared) fitted parameters for all particle species, can produce the physically meaningful collective parameters of a collision system [16,21,33,34]. The combined (global) fits have proved to be quite efficient for the extraction of collective properties and the comparison of different collision systems with the help of few parameters [16,21,33,34].
The modern LHC experiments have the goal to create the Quark-Gluon Plasma (QGP) in collisions of heavy ions at high energies to investigate the various properties of QGP matter, which is assumed to have been created a few microseconds after the so-called “Big Bang”, thought to be a starting point for the birth of our Universe. Nonetheless, analysis of tiny p + p collision systems at high energies is also important and interesting. Such analyses are needed not only because the results from proton-proton collisions are used as a baseline for the investigation of heavy-ion collisions, but also to study the collective properties of a system produced in proton-proton collisions at the largest energies achieved at modern collider experiments [15]. The observations [21,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50] of different QGP signals in high-multiplicity proton-proton collisions at the LHC, including strong similarity and resemblance of their collective properties to those of heavy ions, serve as a good motivation to further investigate the high-energy proton-proton collisions.
In the present article, we analyze the midrapidity (mid-y) pt spectra of given particle species at ten groups of the average charged-particle (pseudo-rapidity) multiplicity density (<dNch/>) in proton-proton collisions at (s)1/2 = 7 TeV at the LHC, measured by the ALICE collaboration and included in [45]. The goal of the present paper is to extract information on the evolution of collective parameters of the p + p system with changing <dNch/> through combined (simultaneous) fits of pt distributions of the charged pions and kaons, protons and antiprotons in each group of <dNch/>, applying the Tsallis distribution with thermodynamical consistence and the Hagedorn function with included transverse flow implemented over the entire measured long pt region. The present work follows and uses the same analysis methods as the recent paper [21], in which the mid-y pt distributions of the same species of particles at ten different groups of <dNch/> in proton-proton collisions at (s)1/2 = 13 TeV at the LHC, measured by ALICE collaboration and presented in [51], have been investigated. Therefore, the results obtained in the present work for proton-proton collisions at (s)1/2 = 7 TeV are compared systematically with the corresponding results in [21] extracted for proton-proton collisions at (s)1/2 = 13 TeV.

2. The Data and Methods

ALICE collaboration has selected [45] inelastic proton-proton collision events at (s)1/2 = 7 TeV with at least one produced charged particle in the pseudo-rapidity region | η | < 1, which make up around 75% of the total inelastic cross-section. To study the dependence of light-flavor hadron production on multiplicity, the selected collision events have been divided into event classes (groups) according to the total charge accumulated in two V0 detectors (V0M amplitude). The V0M amplitude scales linearly with the total number of the charged particles produced in the pseudo-rapidity window corresponding to acceptance of V0 scintillators [45]. The respective average charged-particle pseudo-rapidity densities, <dNch/>, in each class of events have been measured at mid-pseudo-rapidity ( | η | < 0.5). The low pt part of the particle spectra has been reconstructed in the rapidity interval | y | < 0.5, while the high pt part has been extracted in the pseudo-rapidity range | η | < 0.8 in order to make use of the full statistics of inelastic proton-proton collisions [45]. The average charged-particle multiplicity density and corresponding fractions of inelastic cross-sections for different groups of event multiplicity are shown in Table 1. The pt ranges [45] of midrapidity pt distributions in inelastic proton-proton collisions at (s)1/2 = 7 TeV are as follows: [0.1–20.0] GeV/c for π + + π , [0.2–20.0] GeV/c for K + + K , and [0.3–20.0] GeV/c for p + p ¯ . These pt ranges coincide with those measured [21] by ALICE collaboration at midrapidity in ten groups of <dNch/> in inelastic proton-proton collisions at (s)1/2 = 13 TeV.
It is well-known that the high pt part of the pt spectra of particles in high-energy nucleon–nucleon collisions can be reproduced well by the QCD-inspired Hagedorn function [52]:
d 2 N 2 π N e v p t d p t d y = C ( 1 + m t p 0 ) n
where C is the normalization constant, p0 and n are free parameters, m t = p t 2 + m 0 2 is the transverse energy, and m0 is the rest mass of a particle. As mentioned above, the various forms of the Tsallis function [1,2] can well-reproduce the pt and mt spectra of particles in proton-proton collisions at high energies [3,4,5,6,7,8,9,10,11,12,13,22,53]. In its simplest version [10,11,20], the Tsallis function is given at mid-y (y ≈ 0) by:
d 2 N 2 π N e v p t d p t d y = C x ( 1 + ( q 1 ) m t T ) 1 / ( q 1 )
where C x is the normalization constant, T is the effective temperature, and q is the non-extensivity parameter, characterizing the deviation of pt distribution from the exponential Boltzmann–Gibbs distribution. The function in Equation (2) is called a non-extensive generalization of the Boltzmann–Gibbs distribution ( ~ exp ( E T )), with the q parameter added to the temperature parameter. The parameter q is assumed to also be a measure of non-thermalization [54]. When q approaches one (1), the Tsallis function tends to the usual Boltzmann–Gibbs distribution. The degree of system thermalization increases as q approaches one.
Compared to various versions of Tsallis distributions, the following one [3,9,10] at mid-y leads to consistent thermodynamics for the pressure, energy density, and particle number:
d 2 N 2 π N e v p t d p t d y = C q m t ( 1 + ( q 1 ) m t T ) q / ( q 1 )
which we call the Tsallis function with thermodynamical consistence in the present analysis. The values of T, extracted from Tsallis function (Equations (2) and (3)) fits, represent the effective temperatures, which include contributions of both the thermal motion and collective expansion effects. To separate these two factors, the transverse flow velocity is incorporated into the Tsallis distribution function. It has been demonstrated in [10,16] that the Tsallis distribution with thermodynamical consistence, provided in Equation (3), results in noticeably smaller values of T compared to T obtained using the simple non-consistent Tsallis function in Equation (2). This is due to the additional mt term in Equation (3).
The functions represented by Equations (1) and (2) are mathematically equivalent when making substitutions, n = 1 / ( q 1 ) and p 0 = n T . The larger n values correspond to lower q values. For the quark-quark point scattering, n ≈ 4, and the parameter n becomes larger in case of the multiple scattering centers involved [6,55,56]. To embed the effective temperature, T, the function in Equation (1) is modified, setting p 0 = n T :
d 2 N 2 π N e v p t d p t d y = C ( 1 + m t n   T ) n
In the present work, we include the transverse flow into Equation (4) by using the transformation m t γ t ( m t p t β t ) , as also performed in [11,15,16,21,33]. Then, Equation (4) is modified to:
d 2 N 2 π N e v p t d p t d y = C ( 1 + γ t ( m t p t β t ) n T 0 ) n
Here, γ t = 1 / 1 β t 2 , < β t > is the mean transverse flow velocity, and T0 is an estimated kinetic freeze-out temperature. We call the function in Equation (5) the Hagedorn function with included transverse flow in the present work.
Successful analysis, using the function in Equation (5), of pt distributions of particles in collisions at high energies has been performed in [11,15,16,21,33]. The model represented by the function in Equation (5) is a powerful tool, which probes the long pt ranges of particles, allowing to compare different collision systems using few parameters [15,16,21,33]. The global parameters < β t > and T0 obtained in [33] in central Cu + Cu, Au + Au, and Pb + Pb collisions at mid-y at (snn)1/2 = 200–2760 GeV at the RHIC and LHC, applying the modified Hagedorn function in Equation (5) over the long pt range, have reproduced all the observed dependencies of < β t > and T0 on <Npart> and (snn)1/2, extracted using three different transverse expansion (blast-wave) models in the low pt range. The values of < β t > in the most central (0–5%) Pb + Pb collisions at (snn)1/2 = 2.76 and 5.02 TeV, extracted using the function in Equation (5) over long pt ranges in [16], have agreed within fit errors with the respective values of β t , obtained by ALICE collaboration in the most central (0–5%) Pb + Pb collisions at (snn)1/2 = 2.76 and 5.02 TeV in [34] and [57], respectively, using the combined Boltzmann–Gibbs blast-wave fits applied over low pt regions. Simultaneous fits with the modified Hagedorn function in Equation (5), applied over long pt intervals in [16], have confirmed that β t grows and T0 decreases with increasing the collision centrality in Pb + Pb collisions at (snn)1/2 = 2.76 and 5.02 TeV, which agrees well with the similar result of ALICE collaboration obtained from combined Boltzmann–Gibbs blast-wave fits in the low pt ranges in the same collisions in [34,57].
In the present analysis, for the description of pt distributions, d2N/(Nevdptdy), of the charged particles in ten groups of event multiplicity in inelastic proton-proton collisions at (s)1/2 = 7 TeV, we used the Tsallis function with thermodynamical consistence (Equation (3)) in the following form:
d 2 N N e v d p t d y = 2 π C q p t m t ( 1 + ( q 1 ) m t T ) q q 1
Similarly, to describe the pt spectra, d2N/(Nevdptdy), of particles in the present work, we applied the Hagedorn function with included transverse flow (Equation (5)) in the following form:
d 2 N N e v d p t d y = 2 π C p t ( 1 + γ t ( m t p t β t ) n T 0   ) n
In the present work, the simultaneous fits by the model functions of pt distributions of all studied particle types in each group of <dNch/> have been conducted using the Nonlinear Curve Fitting of the Origin 9.1 Graphing and Data Analysis Software. The T/T0 and <βt> are extracted as global (common) parameters for all particle types in the fit procedures.
The combined statistical and systematic errors (added in quadrature) are shown in the figures for the experimental data points. The details for the calculation of systematic errors in the pt spectra of the particles are provided in [45]. The minimum χ2 fit processes have been conducted accounting for the combined statistical and systematic errors as the weights (1/(error)2) for data points. During fits of pt distributions of the charged pions, the range pt < 0.5 GeV/c, containing significant contribution from decays of baryon resonance, has been excluded, as is also the case in [15,16,21,33,34,57]. It is important to mention that in the present work, we used the same model functions, presented in Equations (6) and (7), and identical fitting pt ranges for the studied particle species as used in recent work [21] for proton-proton collisions at (s)1/2 = 13 TeV.

3. Analysis and Results

The results extracted from simultaneous minimum χ2 fits by the Tsallis function with thermodynamical consistence (Equation (6)) of pt distributions of the studied particle types in various groups of <dNch/> in proton-proton collisions at (s)1/2 = 7 TeV are shown in Table 2. The parameter T is extracted as a shared (common) fitted parameter for all the particle species during the combined fitting procedures by Equation (6) in each group of <dNch/>. The results of this global fitting procedure are shown on Figure 1 for four different groups of <dNch/>.
Table 3 shows the parameters extracted from combined minimum χ2 fits with the Hagedorn function with included transverse flow (Equation (7)) of pt distributions of particles in different groups of <dNch/> in proton-proton collisions at (s)1/2 = 7 TeV. The parameters β t and T0 are extracted as shared (common) fitted parameters for all the particle species during the combined fitting procedures by Equation (7) in each group of <dNch/>. Figure 2 illustrates the respective combined minimum χ2 fit curves by the Hagedorn function with included transverse flow of mid-y pt distributions of the studied charged particles in four different groups of <dNch/>. Figure 1 and the χ 2 /n.d.f. values in Table 2 show that the combined fits with the Tsallis function with thermodynamical consistence reproduce quite satisfactory pt distributions of hadrons, except that for (anti)protons at high pt values in the low-multiplicity events, in ten various groups of <dNch/> in proton-proton collisions at (s)1/2 = 7 TeV. Similarly, as observed from Figure 2 and the values of χ 2 /n.d.f. in Table 3, the simultaneous fits with the Hagedorn function with included transverse flow reproduce pt distributions of particles quite well, except that for (anti)protons at high pt values in the low-multiplicity events, in ten groups of <dNch/>.
As one can see from Figure 1 and Figure 2, the pt spectra of (anti)protons could not be described well only in the low-multiplicity p + p collisions in region pt > 6–8 GeV/c. Even in these low-multiplicity (with low <dNch/> < 6 values) p + p collision events at (s)1/2 = 7 TeV, the quality of the combined fits is quite satisfactory, as one can see from the obtained χ2/n.d.f. values in Table 2 and Table 3. We attribute these deviations from model fits in the pt spectra of (anti)protons in region pt > 6–8 GeV/c in low-multiplicity p + p collisions to the fact that the protons are already available (as compared to the newly produced pions and kaons) as initial (colliding) particles, and it requires significantly larger energy to be deposited to the collision system (and hence, a significantly higher degree of thermalization) for the production of new (anti)protons, consisting of three (anti)quarks, as compared to that required for the production of pions and kaons, consisting of two (anti)quarks. Practically no thermalized collective system with collective transverse flow is expected or observed (see Table 3) in peripheral p + p interactions (with low <dNch/> < 6 values) with low energy transferred to the system.
Figure 3 illustrates the <dNch/> dependencies of the obtained q and T parameters of the Tsallis function with thermodynamical consistence for the studied charged particles in proton-proton collisions at (s)1/2 = 7 TeV, extracted from simultaneous minimum χ2 fits and presented in Table 2. Figure 4 summarizes the <dNch/> dependencies of the extracted parameters β t , T0, and n of the Hagedorn function with included transverse flow for the studied particles in proton-proton collisions at (s)1/2 = 7 TeV, presented in Table 3. The corresponding data obtained in [21] for the respective particle species in proton-proton collisions at (s)1/2 = 13 TeV are also presented in Figure 3 and Figure 4 for comparison. Here, it is worth mentioning the importance of analyzing the behavior of β t and T0 and respective excitation functions [14] because of their relation to map the QCD phase diagram, even though the chemical freeze-out temperature (Tch) is normally used in such phase diagrams.
As observed in Figure 3a, the overall values of q for mesons (pions and kaons) and baryons (protons and antiprotons) are clearly separated with a q(baryons) < q(mesons) relation holding in the entire studied <dNch/> region in proton-proton collisions at (s)1/2 = 7 and 13 TeV. This is consistent with the q(baryons) < q(mesons) relation in minimum bias, high-energy proton-proton collisions obtained earlier in [15,16,22,25]. As observed in Figure 3a, the behavior of q versus <dNch/> dependence for all considered particle types in proton-proton collisions at (s)1/2 = 7 TeV is similar to that obtained in proton-proton collisions at (s)1/2 = 13 TeV in [21]. However, as observed in Figure 3a, on the whole, the values of q for the analyzed particle species in proton-proton collisions at (s)1/2 = 7 TeV are found to be noticeably lower compared to those in proton-proton collisions at (s)1/2 = 13 TeV in the whole analyzed <dNch/> range. This can suggest that the systems produced in proton-proton collisions at (s)1/2 = 7 TeV are characterized by the noticeably larger degree of equilibrium and thermalization as compared to that in proton-proton interactions at (s)1/2 = 13 TeV. This finding is consistent with that of the recent work [16], in which the non-extensivity parameter q has been shown to increase systematically for all studied particle species with increasing the energy (s)1/2 of proton-proton collisions from 2.76 to 5.02 TeV. These results suggest that higher-energy proton-proton collisions are characterized by the lower degree of thermalization (larger degree of non-equilibrium).
In [8], the combined pt spectra of the charged particles, consisting predominantly of pions, in minimum bias p + p interactions at wide energy regions (s)1/2 = 0.54–7 TeV have been analyzed using the thermodynamically consistent Tsallis distribution. It has been obtained [8] that q increases weakly but consistently with beam energy, reaching the highest value ≈ 1.15 in minimum bias proton-proton interactions at (s)1/2 = 7 TeV. This value agrees with q ≈ 1.15 for pions and kaons in Table 2 in the majority of multiplicity classes at (s)1/2 = 7 TeV in the present analysis. The effective Tsallis temperature of 82 ± 1 MeV obtained in [8] in minimum bias proton-proton interactions at (s)1/2 = 7 TeV agrees within uncertainties with the average effective Tsallis temperature <T> = 81 ± 1 MeV calculated in the present analysis. The average effective Tsallis temperature <T> for minimum bias proton-proton interactions at (s)1/2 = 7 TeV in the present work has been calculated using the extracted temperature, Ti, from Table 2 and the corresponding fraction of events, fi (with respect to the total inelastic cross-section), of the i-th group of <dNch/> from Table 1 as follows: < T > = i = 1 10 ( T i · f i ) .
As seen from Figure 3b, the extracted T values of the Tsallis function with thermodynamical consistence show a similar consistent growth with increasing <dNch/> in proton-proton interactions at (s)1/2 = 7 TeV as that observed in collisions of protons at (s)1/2 = 13 TeV. This result can be explained by the fact that the larger multiplicity collisions are consistent with harder proton-proton collisions, in which a larger amount of energy-momentum is deposited to a system [21]. We have fitted the data points in Figure 3b, extracted from proton-proton collisions at (s)1/2 = 7 TeV, employing the simple power function T = A · < d N c h d η > α , where A is the normalization constant and α is the exponent parameter, as was carried out in [21]. Similar to [21], we have conducted the minimum χ2 fits by this simple power function in the entire studied <dNch/> region, and in the entire <dNch/> region excluding the first data point, i.e., excluding the softest (ultra-peripheral) proton-proton collisions. The respective results from minimum χ2 fits are presented in Table 4. The corresponding data extracted from [21] for proton-proton interactions at (s)1/2 = 13 TeV are also shown in Table 4 for comparison. Figure 3b and χ 2 /n.d.f. values in Table 4 show that this power function is not able to well-reproduce the T versus <dNch/> dependence in the entire studied region. However, as seen from Figure 3b and values of χ 2 /n.d.f. from Table 4, the T versus <dNch/> dependence in proton-proton collisions at (s)1/2 = 7 TeV is described quite satisfactorily, with the function T = A · < d N c h d η > α having exponent parameter ≈ 1/3 in region <dNch/> > 3, when the first data point is excluded from the fit range. This agrees very well, as observed from Figure 3b and Table 4, with the similar result in [21], in which the T versus <dNch/> dependence in proton-proton interactions at (s)1/2 = 13 TeV has been described quite well with the function T = A · < d N c h d η > α with exponent parameter ≈ 1/3 in the whole <dNch/> region, excluding the first data point.
It was observed from Figure 4a that β t values, obtained employing the Hagedorn function with included transverse flow, are essentially zero in the range <dNch/> < 6 (the first three data points) in proton-proton interactions at (s)1/2 = 7 TeV. As observed in Figure 4a, beginning from <dNch/> ≈ 6, the transverse flow starts emerging and developing, with parameter β t growing steadily in region <dNch/> > 6 until the largest <dNch/> values. The corresponding parameter T0 grows systematically in the low-multiplicity region <dNch/> < 6 (the first three data points) and does not change within the fit uncertainties in the broad region <dNch/> > 6 in proton-proton interactions at (s)1/2 = 7 TeV. The broad plateau range of T0 starting at <dNch/> ≈ 6 in Figure 4b goes along with a simultaneous and systematic rise of β t with increasing <dNch/> in the wide region <dNch/> > 6 in Figure 4a. The observation that the temperature, T0, first rises in the range <dNch/> < 6, reaching a wide plateau beginning at <dNch/> ≈ 6, which matches the absence of transverse flow at <dNch/> < 6 and the emergence and systematic increase of β t at the wide range <dNch/> > 6, suggests the probable onset of deconfinement phase transition at <dNch/> ≈ 6 in proton-proton interactions at (s)1/2 = 7 TeV. The β t and T0 versus <dNch/> dependencies in collisions of protons at (s)1/2 = 7 TeV, as observed from Figure 4a,b, qualitatively match those in collisions of protons at (s)1/2 = 13 TeV obtained in [21], with some differences observed in the absolute values of parameters. It is seen in Figure 4a,b that the absolute values of β t are noticeably smaller and the <dNch/> value for the probable onset of phase transition and the value of T0 in the plateau region are larger in p + p collisions at (s)1/2 = 13 TeV as compared to those at (s)1/2 = 7 TeV. It has been estimated in [21] from the analysis of T0 and β t versus <dNch/> dependencies that the expected onset of deconfinement phase transition in proton-proton collisions at (s)1/2 = 13 TeV takes place at <dNch/> ≈ 7.1 ± 0.2. We have followed a similar procedure as that in [21] to better estimate the critical value of <dNch/> for the probable onset of deconfinement phase transition in collisions of protons at (s)1/2 = 7 TeV: The critical <dNch/> is estimated to be the middle value on the <dNch/> axis between the third and fourth points in Figure 4b, which correspond to <dNch/> = 5.40 ± 0.17 and <dNch/> = 6.72 ± 0.21, respectively, in Table 1. The respective estimate for the probable onset of deconfinement phase transition in p + p collisions at (s)1/2 = 7 TeV proved to be <dNch/> ≈ 6.1 ± 0.3, which is smaller than, but close to, the corresponding <dNch/> ≈ 7.1 ± 0.2 obtained in p + p collisions at (s)1/2 = 13 TeV in [21].
As observed in Figure 4c, the n versus <dNch/> dependencies, obtained by employing the Hagedorn function with included transverse flow, in proton-proton collisions at (s)1/2 = 7 and 13 TeV [21] approximately match the respective inverse dependencies of the q parameter for the studied particles in Figure 3a.

4. Further Analysis and Discussion

It has been observed that β t becomes considerable at larger multiplicities in proton-proton collisions at (s)1/2 = 7 TeV, attaining the maximal value (see Table 3) of β t = 0.29 ± 0.02 at the highest studied value of <dNch/> = 21.3 ± 0.6. The transverse flow has been absent in proton-proton collisions in the low-multiplicity region <dNch/> < 6 (see Table 3 and Figure 4a). These results are consistent with the similar emergence and development of transverse radial flow in higher <dNch/> events in collisions of protons at (s)1/2 = 13 TeV at the LHC recently deduced in [21]. As shown in Table 1, the β t ≈ 0 events with <dNch/> < 6 make up around 62% of the total statistics of inelastic proton-proton collisions at (s)1/2 = 7 TeV. In recent papers [15,16], almost zero transverse flow velocity has been extracted in minimum bias inelastic proton-proton collisions at (s)1/2 = 2.76, 5.02, and 7 TeV, analyzing mid-y pt distributions of particles with the Hagedorn function with embedded transverse flow. The practically zero β t values obtained in [15,16] in minimum bias inelastic proton-proton interactions at (s)1/2 = 2.76, 5.02, and 7 TeV are likely due to the dominant fraction of events with low-multiplicity with the absence of transverse flow in the total ensemble of inelastic proton-proton collisions.
The various features of QGP: the enhancement of strangeness [38], hardening of pt distributions [39,40], and other QGP characteristics [35,36,37,41,42,43,44,45,46,47,48,49,50,58], obtained in proton-proton collisions at the LHC have been published. The onset of collective radial expansion in proton-proton collisions at (s)1/2 = 900 and 7000 GeV at the LHC have been reported in [35]. The emergence and development of transverse radial flow in larger multiplicity proton-proton events at (s)1/2 = 13 TeV at the LHC have been reported in [21]. The signatures of an equilibrated and collective system in large-multiplicity proton-proton events at (s)1/2 = 7 TeV at the LHC were discovered in [45]. In [38], the integrated yields of strange and multi-strange particles, relative to pions, have been shown to increase significantly with increasing <dNch/> in proton-proton interactions at (s)1/2 = 7 TeV at the LHC. This has been the first reported finding [38] of strangeness enhancement in large-multiplicity p + p events. These measurements [38] of ALICE collaboration have proven to be in remarkable agreement with those of p + Pb collisions at (snn)1/2 = 5.02 TeV at LHC, implying that the observed phenomenon is due to the final system produced in a collision. In large-multiplicity p + p collisions [38], the strangeness production has reached values similar to those obtained in Pb + Pb collisions at (snn)1/2 = 2.76 TeV, in which QGP is created. In a recent analysis [44], pt distributions of hadrons as a function of charged-particle multiplicity and transverse sphericity have been investigated, using non-extensive Tsallis statistics and the Boltzmann–Gibbs blast-wave (BGBW) model in collisions of protons at (s)1/2 = 13 TeV, with the help of the PYTHIA8 event generator. It has been deduced that the isotropic proton-proton collisions approach to thermal equilibrium, while the jetty events remain far from equilibrium [44].
The recent analyses in [46,47,48,49,50] present convincing explanations regarding the source of thermalization in proton-proton collisions at high energies. Thermal abundances of particle yield in proton-proton collisions at high energies with an exponential component in pt distributions of hadrons are considered to be a clear signature of thermalization [46,50,58,59]. Since a few secondary interactions in high-energy proton-proton collisions do not support thermalization via interactions at the final state, the discovery of a thermal feature in proton-proton collisions has been quite surprising [46,50]. The thermalization in high-energy proton-proton interactions has been explained in [46,47,48,49] as the process occurring during the rapid quench caused by collisions due to the large extent of quantum entanglement inside the wave functions of partons of colliding protons. Therefore, the effective T extracted from the pt spectra of hadrons depends on momentum transfer, constituting an ultraviolet cutoff of quantum modes resolved by collision [46,50]. In [50], pt distributions at various multiplicities of Pb + Pb collisions at (snn)1/2 = 2.76 TeV and proton-proton collisions at (s)1/2 = 7 TeV have been fitted well by the sum of an exponential and power-like function, represented by thermal-like temperature, Tth, and hard temperature, Th. The thermalization caused by quantum entanglement has been confirmed in [50], analyzing the multiplicity dependence of proton-proton and Pb + Pb collision data at the LHC. No thermal radiation is expected in diffractive proton-proton events with a large rapidity gap, because in such events the whole proton wave function becomes involved and entanglement entropy does not occur [50]. It has been indeed found in [46] that the thermal component vanishes in diffractive proton-proton interactions at (s)1/2 = 13 TeV, in spite of many hadrons still being created. In the studied proton-proton and Pb + Pb collision data at LHC, the effective thermalization temperature, Tth, proved to be proportional [50] to the hard temperature, Th, determined by the average pt. The obtained coefficient of proportionality has been shown as being universal and independent from the type of collision [50]. The proportionality between Tth and Th has been interpreted as that arising from color source clustering [50]. The Tth and Th have been shown to rise with increasing multiplicity of both p + p and Pb + Pb collisions [50].
The systematic rise of parameter T (extracted employing the Tsallis function with thermodynamical consistence) with increasing particle multiplicity in proton-proton collisions at (s)1/2 = 7 and 13 TeV [21], as seen from Figure 3b, is in agreement with the increase of Tth with multiplicity in collisions of protons at the LHC in [50]. We have confirmed an interesting result of [21], namely that T versus <dNch/> dependence in collisions of protons at (s)1/2 = 13 TeV is described well with the function T = A · < d N c h d η > α with exponent parameter ≈ 1/3. In the present paper, we have obtained that T versus <dNch/> dependence in p + p collisions at (s)1/2 = 7 TeV is also described quite satisfactorily, with this power function having the exponent parameter ≈ 1/3 in the wide region <dNch/> > 3. Hence, assuming the proportionality of the charged-particle multiplicity density to the initial energy density (ε), it is deduced that T ~ ε 1 / 3 in proton-proton collisions at (s)1/2 = 7 and 13 TeV. The linear proportionality of the charged-particle multiplicity density to the initial energy density (ε) has been demonstrated in Figure 3 of [42], in which advanced calculations of the initial energy density of (s)1/2 = 7 and 8 TeV p + p collisions have been performed, using the accelerating, exact, and explicit solutions of relativistic hydrodynamics [42]. It is interesting to mention that in a simple model of an ideal gas of massless pions [60], the energy density as a function of temperature is described by the Stefan–Boltzmann form ε π = ( π 2 10 ) T 4 , and hence T ~ ε π 1 / 4 . Comparing T ~ ε 1 / 3 extracted in the present analysis and in [21] with the relation T ~ ε π 1 / 4 for the model of an ideal gas of massless pions, one can see that the energy density dependencies of the effective temperatures of the systems, formed in proton-proton collisions at (s)1/2 = 7 and 13 TeV [21], and of the temperature of an ideal gas of massless pions, are compatible with each other due to the closeness of the respective exponent parameters.
Employing the explicit and exact relativistic hydrodynamic solutions, the advanced calculations of initial energy density in minimum bias proton-proton interactions at mid-y at (s)1/2 = 7 and 8 TeV at the LHC were performed in [42]. The improved advanced hydrodynamic evaluations [42] of initial energy density in minimum bias proton-proton collisions at mid-y at the LHC have yielded ε(7 TeV) = 0.645 GeV/fm3 and ε(8 TeV) = 0.641 GeV/fm3, respectively, which are below the critical energy density (1 GeV/fm3) from lattice QCD calculations. These results were extracted using the mean values of the charged-particle multiplicity density at mid-y in minimum bias proton-proton collisions at (s)1/2 = 7 and 8 TeV at the LHC ALICE and CMS experiments. It has been concluded [42] that a large enough initial energy density to create a non-hadronic perfect fluid (QGP-like) is available in high-multiplicity proton-proton events at the LHC. The analogous advanced calculations, employing the relativistic hydrodynamic solutions in [43], have yielded ε(13 TeV) ≈ 0.69 GeV/fm3 using <dNch/> ≈ 6.5 in minimum bias proton-proton collisions at (s)1/2 = 13 TeV at mid-y at the LHC. Assuming the proportionality of the charged-particle multiplicity density to the initial energy density (ε) and employing the results of [43], the critical value of <dNch/> at mid-y in proton-proton collisions at (s)1/2 = 13 TeV, required to achieve the critical QCD energy density (1 GeV/fm3), has been calculated in [21] as follows:
< d N c h d η > ( p + p   at   13   TeV ) = 1   GeV / fm 3 0.69   GeV / fm 3 · 6.5 9.4
Similarly, using the results of [42], we can evaluate the critical <dNch/> at mid-y in proton-proton collisions at (s)1/2 = 7 TeV to achieve the critical QCD energy density (1 GeV/fm3):
< d N c h d η > ( p + p   at   7   TeV ) = 1   GeV / fm 3 0.645   GeV / fm 3 · 5.9 9.1
Our estimate of <dNch/> ≈ 6.1 ± 0.3 for the probable onset of deconfinement phase transition was found to be lower than the critical value of <dNch/> for achieving the QCD critical energy density (1 GeV/fm3), evaluated above (Equation (9)) in collisions of protons at (s)1/2 = 7 TeV at mid-y. It is an interesting result that <dNch/> ≈ 6.1 ± 0.3 for the probable deconfinement phase transition estimated in the present work in collisions of protons at (s)1/2 = 7 TeV proved to be smaller than, but close to, the corresponding <dNch/> ≈ 7.1 ± 0.2 estimated in [21] in collisions of protons at (s)1/2 = 13 TeV, which is also observed in Figure 4a,b. Using the results of Equations (8) and (9), in the present analysis we can estimate the respective critical energy densities for the probable onset of deconfinement phase transition in proton-proton collisions at (s)1/2 = 7 and 13 TeV corresponding to <dNch/> ≈ 6.1 ± 0.3 and <dNch/> ≈ 7.1 ± 0.2, respectively:
ε t r ( p + p   at   7   TeV ) = 6.1 9.1 × 1 GeV fm 3 0.67 ± 0.03   GeV fm 3
and
ε t r ( p + p   at   13   TeV ) = 7.1 9.4 × 1 GeV fm 3 0.76 ± 0.02   GeV fm 3
As seen from the above calculations, the estimated critical energy densities for the probable onset of deconfinement phase transitions in proton-proton collisions at (s)1/2 = 7 and 13 TeV have proven to be significantly smaller than the critical value of 1 GeV/fm3 from lattice QCD calculations. From the comparison of non-extensivity parameter q versus <dNch/> dependencies for two collision energies in Figure 3a, it has been deduced that the systems produced in collisions of protons at (s)1/2 = 7 TeV are characterized by a noticeably larger degree of equilibrium and thermalization as compared to that in proton-proton collisions at (s)1/2 = 13 TeV, which could be due to the faster and more violent proton-proton collisions at (s)1/2 = 13 TeV than those at (s)1/2 = 7 TeV. The deduced larger degree of equilibrium and thermalization at (s)1/2 = 7 TeV, compared to that at (s)1/2 = 13 TeV, can probably explain our finding that the probable onset of deconfinement phase transition in proton-proton collisions at (s)1/2 = 7 TeV takes place at a lower estimated critical energy density (0.67 ± 0.03 GeV/fm3) than that (0.76 ± 0.02 GeV/fm3) in collisions of protons at (s)1/2 = 13 TeV.
Here, it is worth mentioning the results of relevant works [61,62,63] of Campanini et al. In an early paper [61], the available experimental data on average pt (<pt>) versus charged-particle rapidity density, dn/dy, dependencies in proton-proton and proton–antiproton collisions at SPS collider and ISR have been examined. It has been suggested [61] that evidence may already be available of high-temperature phase transition to quark-gluon plasma in hadronic collision events at SPS and ISR with central rapidity density dn/dy within one or two units around 6. In [62], it has been shown that in various experiments in proton-proton and proton–antiproton collisions at (s)1/2 energies ranging from 22 GeV to 7 TeV, there is a well-defined slope change in <pt> versus charged-particle rapidity density, which may indicate a phase transition to a new production mechanism of particles at dn/ = 5.5 ± 1.2. In [63], the approximate equations of state (EOS) have been extracted from the analysis of <pt> versus dNch/ dependencies in experimental data on proton-proton and proton–antiproton interactions at (s)1/2 energies ranging from 31 GeV to 7 TeV. The obtained results [63] proved to be consistent with theory predictions in case of crossover from hadronic gas to quark gluon plasma, starting at (dNch/) ≈ 6 and ending at (dNch/) ≈ 24. Hence, the critical values <dNch/> ≈ 6.1 ± 0.3 and <dNch/> ≈ 7.1 ± 0.2 [21] for the probable onset of deconfinement phase transition in collisions of protons at (s)1/2 = 7 and 13 TeV, respectively, estimated from the analysis of the dependency of T0 and β t on <dNch/>, are consistent and in good agreement with the respective estimates of Campanini, deduced from detailed analyses of the average pt (<pt>) versus charged-particle (pseudo)rapidity density dependencies in [61,62,63].
For more precise estimations of the critical <dNch/> and corresponding critical energy density for the probable onset of deconfinement phase transition in collisions of protons at high LHC energies, it is necessary to conduct a further comprehensive analysis of dependencies of particle production on <dNch/> in proton-proton collisions, preferably with smaller binning in <dNch/>.

5. Summary and Conclusions

We have analyzed experimental pt distributions of particles at ten groups of <dNch/> in inelastic proton-proton collisions at (s)1/2 = 7 TeV at mid-y at the LHC, measured by ALICE collaboration. We have studied the change of collective characteristics of the collision system with varying <dNch/> through combined minimum χ2 model fits of pt distributions of identified charged particles, using the Tsallis distribution function with thermodynamical consistence and the Hagedorn function with included transverse (radial) flow. The combined minimum χ2 fits with the Tsallis function with thermodynamical consistence and the Hagedorn function with transverse (radial) flow quite well-describe the pt distributions of the studied particle species in ten different groups of <dNch/> in collisions of protons at (s)1/2 = 7 TeV. The findings of the present analysis for proton-proton collisions at (s)1/2 = 7 TeV were systematically compared with the respective results of recent work [21] for proton-proton collisions at (s)1/2 = 13 TeV.
The values of non-extensivity parameter q for pions and kaons, protons and antiprotons in proton-proton collisions at (s)1/2 = 7 TeV were found to be noticeably smaller compared to those in collisions of protons at (s)1/2 = 13 TeV in the whole analyzed <dNch/> range. This indicates that the systems produced in proton-proton collisions at (s)1/2 = 7 TeV have a noticeably larger degree of equilibrium and thermalization than those at (s)1/2 = 13 TeV.
The extracted effective temperatures, T, of the Tsallis function with thermodynamical consistence have shown a consistent increase with the increasing <dNch/> in proton-proton collisions at (s)1/2 = 7 TeV, in agreement with the similar result [21] obtained in proton-proton collisions at (s)1/2 = 13 TeV. The respective T versus <dNch/> dependence in proton-proton collisions at (s)1/2 = 7 TeV has been well-described with the simple power function T = A · < d N c h d η > α with the same value ≈ 1/3 of exponent parameter as that obtained recently [21] in proton-proton collisions at (s)1/2 = 13 TeV. Comparing the relation T ~ ε 1 / 3 extracted in the present analysis and in [21] with the relation T ~ ε π 1 / 4 for the simple model of an ideal gas of massless pions, it was found that dependencies on the initial energy density of the effective temperatures of the systems, created in collisions of protons at (s)1/2 = 7 and 13 TeV, and that of an ideal pion gas, are compatible to each other due to the closeness of the corresponding exponent parameters.
It was obtained that the transverse (radial) flow emerged at <dNch/> ≈ 6 and then increased, becoming significant at higher multiplicity events and reaching the maximal value of β t = 0.29 ± 0.02 at the highest studied value of <dNch/> = 21.3 ± 0.6 in proton-proton collisions at (s)1/2 = 7 TeV. These results are consistent with the similar emergence and development of transverse radial flow in larger multiplicity proton-proton events at (s)1/2 = 13 TeV at the LHC demonstrated in [21].
We have estimated, analyzing T0 and β t versus <dNch/> dependencies, extracted employing the Hagedorn function with included transverse flow, that the probable onset of deconfinement phase transition in proton-proton collisions at (s)1/2 = 7 TeV takes place at <dNch/> ≈ 6.1 ± 0.3, which is close to the corresponding estimate (<dNch/> ≈ 7.1 ± 0.2) obtained recently in collisions of protons at (s)1/2 = 13 TeV in [21]. The critical values <dNch/> ≈ 6.1 ± 0.3 and <dNch/> ≈ 7.1 ± 0.2 [21] for the probable onset of deconfinement phase transition in collisions of protons at (s)1/2 = 7 and 13 TeV, respectively, estimated by us from the analysis of T0 and β t versus <dNch/> dependencies, proved to be consistent and in good agreement with the respective estimates of Campanini et al. in [61,62,63] and theory predictions [63] in case of crossover from hadronic gas to quark gluon plasma starting at (dNch/) ≈ 6. We have also estimated the corresponding critical energy densities for the probable onset of deconfinement phase transitions in proton-proton collisions at (s)1/2 = 7 and 13 TeV at the LHC to be 0.67 ± 0.03 and 0.76 ± 0.02 GeV/fm3, respectively, being significantly lower than the critical QCD energy density (1 GeV/fm3). The noticeably larger degree of equilibrium and thermalization deduced at (s)1/2 = 7 TeV than that at (s)1/2 = 13 TeV could probably explain our finding that the possible onset of deconfinement phase transition in proton-proton collisions at (s)1/2 = 7 TeV takes place at the lower estimated critical energy density (0.67 ± 0.03 GeV/fm3), as compared to that (0.76 ± 0.02 GeV/fm3) in proton-proton collisions at (s)1/2 = 13 TeV.

Author Contributions

K.K.O., F.-H.L., K.A.M. and M.Z.S. have contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

The work of F.-H.L. has been supported by the National Natural Science Foundation of China under Grant Nos. 11575103, 12047571, and 11947418, and the Shanxi Provincial Natural Science Foundation under Grant No. 201901D111043. The work of Kh.K.O., K.A.M., and M.Z.S. has been supported by the Ministry of Innovative Development of Uzbekistan within the framework of fundamental project No. F3-20200929146 on analysis of open data on proton-proton and heavy-ion collisions at the LHC.

Data Availability Statement

The data analyzed in this work are included within the manuscript and cited as references at relevant places within the text of the article.

Conflicts of Interest

The authors declare no conflict of interest regarding the article. The funding agencies did not have any role in the design of the study; in collection and analysis of the data; in the interpretation of the results; in the writing of the article, or in the decision to publish the manuscript.

References

  1. Tsallis, C. Possible generalization of Boltzmann-Gibbs statistics. J. Stat. Phys. 1988, 52, 479. [Google Scholar] [CrossRef]
  2. Tsallis, C. Nonadditive entropy: The concept and its use. Eur. Phys. J. A 2009, 40, 257. [Google Scholar] [CrossRef] [Green Version]
  3. Cleymans, J.; Worku, D. The Tsallis distribution in proton-proton collisions at (snn)1/2 = 0.9 TeV at the LHC. J. Phys. G 2012, 39, 025006. [Google Scholar] [CrossRef]
  4. Sena, I.; Deppman, A. Systematic analysis of pT -distributions in p + p collisions. Eur. Phys. J. A 2013, 49, 17. [Google Scholar] [CrossRef]
  5. PHENIX Collab. Measurement of neutral mesons in p+p collisions at (snn)1/2 = 200 GeV and scaling properties of hadron production. Phys. Rev. D 2011, 83, 052004. [Google Scholar] [CrossRef] [Green Version]
  6. Khandai, P.K.; Sett, P.; Shukla, P.; Singh, V. Hadron spectra in p+p collisions at RHIC and LHC energies. Int. J. Mod. Phys. A 2013, 28, 1350066. [Google Scholar] [CrossRef]
  7. Wong, C.Y.; Wilk, G. Tsallis Fits to pT Spectra for pp Collisions at LHC. Acta Phys. Polon. B 2012, 43, 2047. [Google Scholar] [CrossRef]
  8. Cleymans, J.; Lykasov, G.; Parvan, A.; Sorin, A.; Teryaev, O.; Worku, D.S. Systematic properties of the Tsallis distribution: Energy dependence of parameters in high energy pp collisions. Phys. Lett. B 2013, 723, 351. [Google Scholar] [CrossRef] [Green Version]
  9. Cleymans, J. On the Use of the Tsallis Distribution at LHC Energies. J. Phys. Conf. Ser. 2017, 779, 012079. [Google Scholar] [CrossRef]
  10. Zheng, H.; Zhu, L. Comparing the Tsallis Distribution with and without Thermodynamical Description in p-p Collisions. Adv. High Energy Phys. 2016, 2016, 9632126. [Google Scholar] [CrossRef] [Green Version]
  11. Khandai, P.K.; Sett, P.; Shukla, P.; Singh, V. System size dependence of hadron pT spectra in p+p and Au+Au collisions at (snn)1/2 = 200 GeV. J. Phys. G 2014, 41, 025105. [Google Scholar] [CrossRef] [Green Version]
  12. Tsallis, C.; Mendes, R.; Plastino, A.R. The role of constraints within generalized nonextensive statistics. Phys. A 1998, 261, 534–554. [Google Scholar] [CrossRef]
  13. Biró, T.S.; Purcsel, G.; Ürmössy, K. Non-extensive approach to quark matter. Eur. Phys. J. A 2009, 40, 325. [Google Scholar] [CrossRef]
  14. Li, L.-L.; Liu, F.-H.; Olimov, K.K. Excitation Functions of Tsallis-Like Parameters in High-Energy Nucleus–Nucleus Collisions. Entropy 2021, 23, 478. [Google Scholar] [CrossRef] [PubMed]
  15. Olimov, K.K.; Iqbal, A. Systematic analysis of midrapidity transverse momentum spectra of identified charged particles in p+p collisions at (snn)1/2 = 2.76, 5.02, and 7 TeV at the LHC. Int. J. Mod. Phys. A 2020, 35, 2050167. [Google Scholar] [CrossRef]
  16. Olimov, K.K.; Kanokova, S.Z.; Olimov, A.K.; Umarov, K.I.; Tukhtaev, B.J.; Gulamov, K.G.; Yuldashev, B.S.; Lutpullaev, S.L.; Saidkhanov, N.S.; Kosim Olimov, K.; et al. Combined analysis of midrapidity transverse momentum spectra of the charged pions and kaons, pro-tons and antiprotons in p+p and Pb+Pb collisions at (snn)1/2 = 2.76 and 5.02 TeV at the LHC. Mod. Phys. Lett. A 2020, 35, 2050237. [Google Scholar] [CrossRef]
  17. STAR Collab. Strange particle production in p+p collisions at (s)1/2 = 200 GeV. Phys. Rev. C 2007, 75, 064901. [Google Scholar] [CrossRef] [Green Version]
  18. PHENIX Collab. Identified charged hadron production in p+p collisions at (s)1/2 = 200 and 62.4 GeV. Phys. Rev. C 2011, 83, 064903. [Google Scholar] [CrossRef] [Green Version]
  19. ALICE Collab. Production of pions, kaons and protons in pp collisions at (s)1/2 = 900 GeV with ALICE at the LHC. Eur. Phys. J. C 2011, 71, 1655. [Google Scholar] [CrossRef] [Green Version]
  20. CMS Collab. Study of the inclusive production of charged pions, kaons, and protons in pp collisions at (s)1/2 = 0.9, 2.76, and 7 TeV. Eur. Phys. J. C 2012, 72, 2164. [Google Scholar] [CrossRef] [Green Version]
  21. Olimov, K.K.; Liu, F.-H.; Musaev, K.A.; Olimov, K.; Tukhtaev, B.J.; Yuldashev, B.S.; Saidkhanov, N.S.; Umarov, K.I.; Gulamov, K.G. Multiplicity dependencies of midrapidity transverse momentum spectra of identified charged particles in p+p collisions at (s)1/2 = 13 TeV at LHC. Int. J. Mod. Phys. A 2021, 36, 2150149. [Google Scholar] [CrossRef]
  22. Bíró, G.; Barnaföldi, G.G.; Biró, T.S.; Ürmössy, K.; Takács, Á. Systematic Analysis of the Non-Extensive Statistical Approach in High Energy Particle Collisions—Experiment vs. Theory. Entropy 2017, 19, 88. [Google Scholar] [CrossRef]
  23. Shen, K.; Barnaföldi, G.G.; Biró, T.S. Hadron Spectra Parameters within the Non-Extensive Approach. Universe 2019, 5, 122. [Google Scholar] [CrossRef] [Green Version]
  24. Grigoryan, S. Using the Tsallis distribution for hadron spectra in pp collisions: Pions and quarkonia at (snn)1/2 = 5–13000 GeV. Phys. Rev. D 2017, 95, 056021. [Google Scholar] [CrossRef] [Green Version]
  25. Bíró, G.; Barnaföldi, G.G.; Biró, T.S.; Shen, K. Mass hierarchy and energy scaling of the Tsallis–Pareto parameters in hadron productions at RHIC and LHC energies. EPJ Web Conf. 2018, 171, 14008. [Google Scholar] [CrossRef]
  26. Schnedermann, E.; Solfrank, J. Thermal phenomenology of hadrons from 200A GeV S+S collisions. Phys. Rev. C 1993, 48, 2462. [Google Scholar] [CrossRef] [Green Version]
  27. STAR Collab. Systematic measurements of identified particle spectra in pp, d+Au, and Au+Au collisions at the STAR detector. Phys. Rev. C 2009, 79, 034909. [Google Scholar] [CrossRef]
  28. STAR Collab. Identified particle production, azimuthal anisotropy, and interferometry measurements in Au+Au collisions at (snn)1/2 = 9.2 GeV. Phys. Rev. C 2010, 81, 024911. [Google Scholar] [CrossRef]
  29. Tang, Z.; Xu, Y.; Ruan, L.; van Buren, G.; Wang, F.; Xu, Z. Spectra and radial flow in relativistic heavy ion collisions with Tsallis statistics in a blast-wave description. Phys. Rev. C 2009, 79, 051901. [Google Scholar] [CrossRef] [Green Version]
  30. Lao, H.-L.; Liu, F.-H.; Lacey, R.A. Extracting kinetic freeze-out temperature and radial flow velocity from an improved Tsallis distribution. Eur. Phys. J. A 2017, 53, 44. [Google Scholar] [CrossRef] [Green Version]
  31. Bhattacharyya, T.; Cleymans, J.; Khuntia, A.; Pareek, P.; Sahoo, R. Radial flow in non-extensive thermodynamics and study of particle spectra at LHC in the limit of small (q-1). Eur. Phys. J. A 2016, 52, 30. [Google Scholar] [CrossRef]
  32. Thakur, D.; Tripathy, S.; Garg, P.; Sahoo, R.; Cleymans, J. Indication of a Differential Freeze-Out in Proton-Proton and Heavy-Ion Collisions at RHIC and LHC Energies. Adv. High Energy Phys. 2016, 2016, 4149352. [Google Scholar] [CrossRef]
  33. Olimov, K.K.; Kanokova, S.Z.; Olimov, K.; Gulamov, K.G.; Yuldashev, B.S.; Lutpullaev, S.L.; Umarov, F.Y. Average transverse expansion velocities and global freeze-out temperatures in central Cu+Cu, Au+Au, and Pb+Pb collisions at high energies at RHIC and LHC. Mod. Phys. Lett. A 2020, 35, 2050115. [Google Scholar] [CrossRef]
  34. ALICE Collab. Centrality dependence of π, K, and p production in Pb-Pb collisions at (snn)1/2 = 2.76 TeV. Phys. Rev. C 2013, 88, 044910. [Google Scholar] [CrossRef] [Green Version]
  35. Jiang, K.; Zhu, Y.; Liu, W.; Chen, H.; Li, C.; Ruan, L.; Tang, Z.; Xu, Z. Onset of radial flow in p+p collisions. Phys. Rev. C 2015, 91, 024910. [Google Scholar] [CrossRef] [Green Version]
  36. Bashir, I.U.; Bhat, R.A.; Uddin, S. Evidence of collective flow in p-p collisions at LHC. arXiv 2015, arXiv:1502.04185v2. [Google Scholar]
  37. Bashir, I.U.; Parra, R.A.; Bhat, R.A.; Uddin, S. Particle Transverse Momentum Distributions in p-p Collisions at (s)1/2 = 0.9 TeV. Adv. High Energy Phys. 2019, 2019, 8219567. [Google Scholar] [CrossRef] [Green Version]
  38. ALICE Collab. Enhanced production of multi-strange hadrons in high-multiplicity proton-proton collisions. Nat. Phys. 2017, 13, 535. [Google Scholar] [CrossRef] [Green Version]
  39. Tripathy, S.; et al.; [ALICE Collaboration] Energy dependence of φ(1020) production at mid-rapidity in pp collisions with ALICE at the LHC. Nucl. Phys. A 2019, 982, 180. [Google Scholar] [CrossRef]
  40. Dash, A.K. Multiplicity dependence of strangeness and hadronic resonance production in pp and p-Pb collisions with ALICE at the LHC. Nucl. Phys. A 2019, 982, 467–470. [Google Scholar] [CrossRef]
  41. Khuntia, A.; Sharma, H.; Kumar Tiwari, S.; Sahoo, R.; Cleymans, J. Radial flow and differential freeze-out in proton-proton collisions at (s)1/2 = 7 TeV at the LHC. Eur. Phys. J. A 2019, 55, 3. [Google Scholar] [CrossRef]
  42. Csanad, M.; Csörgő, T.; Jiang, Z.-F.; Yang, C.-B. Initial energy density of (snn)1/2 = 7 and 8 TeV p-p collisions at LHC. Universe 2017, 3, 9. [Google Scholar] [CrossRef] [Green Version]
  43. Jiang, Z.; Csanád, M.; Kasza, G.; Yang, C.; Csörgő, T. Pseudorapidity and Initial-energy Densities in p+p and Heavy-ion Collisions at RHIC and LHC. Acta Phys. Pol. B Proc. Suppl. 2019, 12, 261. [Google Scholar] [CrossRef]
  44. Tripathy, S.; Bisht, A.; Sahoo, R.; Khuntia, A.; Panikkassery Salvan, M. Event shape and multiplicity dependence of freeze-out scenario and system thermodynamics in pro-ton+proton collisions at (snn)1/2 = 13 TeV using PYTHIA8. Adv. High Energy Phys. 2021, 2021, 8822524. [Google Scholar] [CrossRef]
  45. ALICE Collab. Multiplicity dependence of light flavor hadron production in pp collisions at (snn)1/2 = 7 TeV. Phys. Rev. C 2019, 99, 024906. [Google Scholar] [CrossRef] [Green Version]
  46. Baker, O.K.; Kharzeev, D.E. Thermal radiation and entanglement in proton-proton collisions at energies available at the CERN Large Hadron Collider. Phys. Rev. D 2018, 98, 054007. [Google Scholar] [CrossRef] [Green Version]
  47. Kharzeev, D.E.; Levin, E.M. Deep inelastic scattering as a probe of entanglement. Phys. Rev. D 2017, 95, 114008. [Google Scholar] [CrossRef] [Green Version]
  48. Berges, J.; Floerchinger, S.; Venugopalan, R. Dynamics of entanglement in expanding quantum fields. J. High Energy Phys. 2018, 2018, 145. [Google Scholar] [CrossRef] [Green Version]
  49. Berges, J.; Floerchinger, S.; Venugopalan, R. Thermal excitation spectrum from entanglement in an expanding quantum string. Phys. Lett. B 2018, 778, 442. [Google Scholar] [CrossRef]
  50. Feal, X.; Pajares, C.; Vazquez, R.A. Thermal behavior and entanglement in Pb-Pb and p–p collisions. Phys. Rev. C 2019, 99, 015205. [Google Scholar] [CrossRef] [Green Version]
  51. ALICE Collab. Multiplicity dependence of π, K, and p production in pp collisions at (snn)1/2 = 13 TeV. Eur. Phys. J. C 2020, 80, 693. [Google Scholar] [CrossRef]
  52. Hagedorn, R. Multiplicities, pT Distributions and the Expected Hadron→Quark-Gluon Phase Transition. Riv. Nuovo Cim. 1983, 6, 1. [Google Scholar] [CrossRef]
  53. Wilk, G.; Wlodarczyk, Z. Power laws in elementary and heavy-ion collisions. A story of fluctuations and nonextensivity? Eur. Phys. J. A 2009, 40, 299. [Google Scholar] [CrossRef] [Green Version]
  54. Wilk, G.; Wlodarczyk, Z. Interpretation of the Nonextensivity Parameter q in Some Applications of Tsallis Statistics and Lévy Distributions. Phys. Rev. Lett. 2000, 84, 2770. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Blankenbecler, R.; Brodsky, S.J.; Gunion, J. Analysis of particle production at large transverse momentum. Phys. Rev. D 1975, 12, 3469. [Google Scholar] [CrossRef] [Green Version]
  56. Brodsky, S.; Pirner, H.; Raufeisen, J. Scaling properties of high pT inclusive hadron production. Phys. Lett. B 2006, 637, 58. [Google Scholar] [CrossRef] [Green Version]
  57. ALICE Collab. Production of charged pions, kaons and (anti-)protons in Pb-Pb and inelastic pp collisions at (s)1/2 = 5.02 TeV. Phys. Rev. C 2020, 101, 044907. [Google Scholar] [CrossRef]
  58. Mishra, A.N.; Paić, G.; Pajares, C.; Scharenberg, R.P.; Srivastava, B.K. Deconfinement and degrees of freedom in pp and A-A collisions at LHC energies. Eur. Phys. J. A 2021, 57, 245. [Google Scholar] [CrossRef]
  59. Becattini, F. An introduction to the statistical hadronization model. In Proceedings of the Quark-Gluon Plasma and Heavy Ion Collisions: Past, Present, Future, Torino, Italy, 8–14 December 2008. [Google Scholar]
  60. Sarkar, S.; Satz, H.; Sinha, B. (Eds.) The Physics of the Quark-Gluon Plasma: Introductory Lectures; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar] [CrossRef] [Green Version]
  61. Campanini, R. Quark gluon plasma and multiplicity dependence of transverse momentum in hadronic collisions. Lett. Al Nuovo Cimento 1985, 44, 343. [Google Scholar] [CrossRef]
  62. Campanini, R. Possible Signals of new phenomena in hadronic interactions at dn/dη = 5.5 ± 1.2. arXiv 2010, arXiv:1012.5219. [Google Scholar]
  63. Campanini, R.; Ferri, G. Experimental equation of state in proton-proton and proton-antiproton collisions and phase transition to quark gluon plasma. Phys. Lett. B 2011, 703, 237. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The combined minimum χ 2 fits (solid curves) with the function in Equation (6) of the experimental mid-y pt distributions of the charged pions (●) and kaons (▲), protons and antiprotons (■) in proton-proton collisions at (s)1/2 = 7 TeV in various groups of <dNch/>: <dNch/> = 21.3 ± 0.6 (a), <dNch/> = 11.5 ± 0.3 (b), <dNch/> = 6.72 ± 0.21 (c), and <dNch/> = 2.26 ± 0.12 (d). The vertical error bars are combined statistical and systematic errors, dominated by the systematic ones.
Figure 1. The combined minimum χ 2 fits (solid curves) with the function in Equation (6) of the experimental mid-y pt distributions of the charged pions (●) and kaons (▲), protons and antiprotons (■) in proton-proton collisions at (s)1/2 = 7 TeV in various groups of <dNch/>: <dNch/> = 21.3 ± 0.6 (a), <dNch/> = 11.5 ± 0.3 (b), <dNch/> = 6.72 ± 0.21 (c), and <dNch/> = 2.26 ± 0.12 (d). The vertical error bars are combined statistical and systematic errors, dominated by the systematic ones.
Universe 08 00174 g001
Figure 2. The combined minimum χ 2 fits (solid curves) with the function in Equation (7) of the experimental mid-y pt distributions of the charged pions (●) and kaons (▲), protons and antiprotons (■) in proton-proton collisions at (s)1/2 = 7 TeV in various groups of <dNch/>: <dNch/ > = 21.3 ± 0.6 (a), <dNch/> = 11.5 ± 0.3 (b), <dNch/> = 6.72 ± 0.21 (c), and <dNch/> = 2.26 ± 0.12 (d). The vertical error bars are combined statistical and systematic errors, dominated by the systematic ones.
Figure 2. The combined minimum χ 2 fits (solid curves) with the function in Equation (7) of the experimental mid-y pt distributions of the charged pions (●) and kaons (▲), protons and antiprotons (■) in proton-proton collisions at (s)1/2 = 7 TeV in various groups of <dNch/>: <dNch/ > = 21.3 ± 0.6 (a), <dNch/> = 11.5 ± 0.3 (b), <dNch/> = 6.72 ± 0.21 (c), and <dNch/> = 2.26 ± 0.12 (d). The vertical error bars are combined statistical and systematic errors, dominated by the systematic ones.
Universe 08 00174 g002
Figure 3. (a)—The <dNch/> dependencies of the obtained q values of the Tsallis function with thermodynamical consistence, presented in Table 2, for the charged pions (●) and kaons (▲), protons and antiprotons (■) in proton-proton collisions at (s)1/2 = 7 TeV. (b)—The same for the obtained effective temperatures T (●) of the Tsallis function with thermodynamical consistence, presented in Table 2. The corresponding results extracted from [21] for the respective charged particles in proton-proton collisions at (s)1/2 = 13 TeV are shown by the corresponding open (hollow) symbols for comparison. The data points in panel (b) are fitted with the simple power function T = A · < d N c h d η > α , where A is the normalization constant, and α is the exponent parameter. The dashed and solid curves are the simple power function fits of the whole range (10 data points), and of the whole range excluding the first data point (9 data points), respectively.
Figure 3. (a)—The <dNch/> dependencies of the obtained q values of the Tsallis function with thermodynamical consistence, presented in Table 2, for the charged pions (●) and kaons (▲), protons and antiprotons (■) in proton-proton collisions at (s)1/2 = 7 TeV. (b)—The same for the obtained effective temperatures T (●) of the Tsallis function with thermodynamical consistence, presented in Table 2. The corresponding results extracted from [21] for the respective charged particles in proton-proton collisions at (s)1/2 = 13 TeV are shown by the corresponding open (hollow) symbols for comparison. The data points in panel (b) are fitted with the simple power function T = A · < d N c h d η > α , where A is the normalization constant, and α is the exponent parameter. The dashed and solid curves are the simple power function fits of the whole range (10 data points), and of the whole range excluding the first data point (9 data points), respectively.
Universe 08 00174 g003
Figure 4. The <dNch/> dependencies of the obtained β t (a) and T0 (b) parameters (●) of the Hagedorn function with included transverse flow, presented in Table 3, in proton-proton collisions at (s)1/2 = 7 TeV. (c)—The same for the obtained n values of the Hagedorn function with included transverse flow, presented in Table 3, for the charged pions (●) and kaons (▲), protons and antiprotons (■). The corresponding results extracted from [21] for the respective particles in proton-proton collisions at (s)1/2 = 13 TeV are shown by the corresponding open (hollow) symbols for comparison.
Figure 4. The <dNch/> dependencies of the obtained β t (a) and T0 (b) parameters (●) of the Hagedorn function with included transverse flow, presented in Table 3, in proton-proton collisions at (s)1/2 = 7 TeV. (c)—The same for the obtained n values of the Hagedorn function with included transverse flow, presented in Table 3, for the charged pions (●) and kaons (▲), protons and antiprotons (■). The corresponding results extracted from [21] for the respective particles in proton-proton collisions at (s)1/2 = 13 TeV are shown by the corresponding open (hollow) symbols for comparison.
Universe 08 00174 g004
Table 1. The mean charged-particle multiplicity density and fraction of inelastic cross-section in ten multiplicity groups in inelastic proton-proton collisions at (s)1/2 = 7 TeV.
Table 1. The mean charged-particle multiplicity density and fraction of inelastic cross-section in ten multiplicity groups in inelastic proton-proton collisions at (s)1/2 = 7 TeV.
V0M Mult. Class<dNch/>σ/σINEL>0 (%)
I21.3 ± 0.60–0.95
II16.5 ± 0.50.95–4.7
III13.5 ± 0.44.7–9.5
IV11.5 ± 0.39.5–14
V10.1 ± 0.314–19
VI8.45 ± 0.2519–28
VII6.72 ± 0.2128–38
VIII5.40 ± 0.1738–48
IX3.90 ± 0.1448–68
X2.26 ± 0.1268–100
Table 2. Parameters obtained using combined minimum χ2 fits with the Tsallis function with thermodynamical consistence (Equation (6)) of pt spectra of particles in various groups of <dNch/> in proton-proton collisions at (s)1/2 = 7 TeV. The fitted pt ranges are [0.5–20.0] GeV/c for π + + π , [0.2–20.0] GeV/c for K + + K , and [0.3–20.0] GeV/c for p + p ¯ . Here, n.d.f. denotes the number of degrees of freedom.
Table 2. Parameters obtained using combined minimum χ2 fits with the Tsallis function with thermodynamical consistence (Equation (6)) of pt spectra of particles in various groups of <dNch/> in proton-proton collisions at (s)1/2 = 7 TeV. The fitted pt ranges are [0.5–20.0] GeV/c for π + + π , [0.2–20.0] GeV/c for K + + K , and [0.3–20.0] GeV/c for p + p ¯ . Here, n.d.f. denotes the number of degrees of freedom.
<dNch/> q   ( π + + π ) q   ( K + + K ) q   ( p + p ¯ ) T (MeV) χ 2 / n . d . f .   ( n . d . f . )
21.3 ± 0.61.143 ± 0.0011.153 ± 0.0011.133 ± 0.001131 ± 21.85 (114)
16.5 ± 0.51.145 ± 0.0011.154 ± 0.0011.131 ± 0.001121 ± 21.54 (114)
13.5 ± 0.41.146 ± 0.0011.155 ± 0.0011.130 ± 0.001114 ± 11.06 (114)
11.5 ± 0.31.147 ± 0.0011.155 ± 0.0011.128 ± 0.001109 ± 10.79 (114)
10.1 ± 0.31.148 ± 0.0011.155 ± 0.0011.128 ± 0.001104 ± 10.65 (114)
8.45 ± 0.251.149 ± 0.0011.156 ± 0.0011.128 ± 0.00198 ± 10.61 (114)
6.72 ± 0.211.150 ± 0.0011.156 ± 0.0011.125 ± 0.00191 ± 10.45 (114)
5.40 ± 0.171.151 ± 0.0011.157 ± 0.0011.123 ± 0.00185 ± 10.39 (114)
3.90 ± 0.141.151 ± 0.0011.156 ± 0.0011.122 ± 0.00175 ± 10.79 (114)
2.26 ± 0.121.149 ± 0.0011.153 ± 0.0011.114 ± 0.00157 ± 11.00 (114)
Table 3. Parameters extracted using combined minimum χ2 fits with the Hagedorn function with included transverse flow (Equation (7)) of pt distributions of particles in various groups of <dNch/> in proton-proton collisions at (s)1/2 = 7 TeV. The fitted pt ranges are [0.5–20.0] GeV/c for π + + π , [0.2–20.0] GeV/c for K + + K , and [0.3–20.0] GeV/c for p + p ¯ .
Table 3. Parameters extracted using combined minimum χ2 fits with the Hagedorn function with included transverse flow (Equation (7)) of pt distributions of particles in various groups of <dNch/> in proton-proton collisions at (s)1/2 = 7 TeV. The fitted pt ranges are [0.5–20.0] GeV/c for π + + π , [0.2–20.0] GeV/c for K + + K , and [0.3–20.0] GeV/c for p + p ¯ .
<dNch/> n   ( π + + π ) n   ( K + + K ) n   ( p + p ¯ ) β t   ( in   c   Units ) T0 (MeV) χ 2 / n . d . f .   ( n . d . f . )
21.3 ± 0.66.80 ± 0.056.61 ± 0.048.11 ± 0.060.29 ± 0.02113 ± 30.60 (113)
16.5 ± 0.56.77 ± 0.056.57 ± 0.048.13 ± 0.060.23 ± 0.02113 ± 30.72 (113)
13.5 ± 0.46.78 ± 0.046.55 ± 0.048.17 ± 0.060.18 ± 0.02114 ± 30.63 (113)
11.5 ± 0.36.76 ± 0.046.52 ± 0.048.24 ± 0.060.15 ± 0.02113 ± 30.54 (113)
10.1 ± 0.36.77 ± 0.046.52 ± 0.048.24 ± 0.060.12 ± 0.02115 ± 30.58 (113)
8.45 ± 0.256.74 ± 0.046.48 ± 0.048.15 ± 0.060.07 ± 0.02115 ± 40.74 (113)
6.72 ± 0.216.74 ± 0.046.47 ± 0.048.29 ± 0.060.04 ± 0.02112 ± 30.60 (113)
5.40 ± 0.176.73 ± 0.046.46 ± 0.038.40 ± 0.060 ± 0.02108 ± 30.56 (113)
3.90 ± 0.146.70 ± 0.056.49 ± 0.048.48 ± 0.090 ± 0.0395 ± 40.97 (113)
2.26 ± 0.126.75 ± 0.066.59 ± 0.059.06 ± 0.130 ± 0.0369 ± 41.19 (113)
Table 4. Parameters extracted from minimum χ2 fits using the simple power function T = A · < d N c h d η > α of the T versus <dNch/> dependencies in Figure 3b in the whole <dNch/> region (I), and in the whole region excluding the first data point (II) with <dNch/> = 2.26 ± 0.12 and <dNch/> = 2.55 ± 0.04 in proton-proton collisions at (s)1/2 = 7 and 13 TeV, respectively. The corresponding data extracted from [21] for proton-proton collisions at (s)1/2 = 13 TeV are shown for comparison.
Table 4. Parameters extracted from minimum χ2 fits using the simple power function T = A · < d N c h d η > α of the T versus <dNch/> dependencies in Figure 3b in the whole <dNch/> region (I), and in the whole region excluding the first data point (II) with <dNch/> = 2.26 ± 0.12 and <dNch/> = 2.55 ± 0.04 in proton-proton collisions at (s)1/2 = 7 and 13 TeV, respectively. The corresponding data extracted from [21] for proton-proton collisions at (s)1/2 = 13 TeV are shown for comparison.
<dNch/>
Fit Range
Collision Type, (s)1/2A (MeV)α χ 2 / n . d . f .   ( n . d . f . )
p + p, 7 TeV45.8 ± 1.40.353 ± 0.0133.67 (8)
Ip + p, 13 TeV [21]49.9 ± 1.70.33 ± 0.013.66 (8)
p + p, 7 TeV48.7 ± 0.60.327 ± 0.0050.37 (7)
IIp + p, 13 TeV [21]51.5 ± 0.90.31 ± 0.011.05 (7)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Olimov, K.K.; Liu, F.-H.; Musaev, K.A.; Shodmonov, M.Z. Multiplicity Dependencies of Midrapidity Transverse Momentum Distributions of Identified Charged Particles in Proton-Proton Collisions at (s)1/2 = 7 TeV at the LHC. Universe 2022, 8, 174. https://doi.org/10.3390/universe8030174

AMA Style

Olimov KK, Liu F-H, Musaev KA, Shodmonov MZ. Multiplicity Dependencies of Midrapidity Transverse Momentum Distributions of Identified Charged Particles in Proton-Proton Collisions at (s)1/2 = 7 TeV at the LHC. Universe. 2022; 8(3):174. https://doi.org/10.3390/universe8030174

Chicago/Turabian Style

Olimov, Khusniddin K., Fu-Hu Liu, Kobil A. Musaev, and Maratbek Z. Shodmonov. 2022. "Multiplicity Dependencies of Midrapidity Transverse Momentum Distributions of Identified Charged Particles in Proton-Proton Collisions at (s)1/2 = 7 TeV at the LHC" Universe 8, no. 3: 174. https://doi.org/10.3390/universe8030174

APA Style

Olimov, K. K., Liu, F. -H., Musaev, K. A., & Shodmonov, M. Z. (2022). Multiplicity Dependencies of Midrapidity Transverse Momentum Distributions of Identified Charged Particles in Proton-Proton Collisions at (s)1/2 = 7 TeV at the LHC. Universe, 8(3), 174. https://doi.org/10.3390/universe8030174

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop