In the past decades, a number of bulk observables about charged particles, such as the Fourier coefficients of azimuth-angle distributions   , transverse momentum spectra  -  and pseudorapidity distributions   , have experienced a series of extensive investigations in both nucleus and hadron collisions. These investigations have confirmed two facts. One is that the quark-gluon plasma (QGP) might have come into being in nucleus collisions at current RHIC or LHC energies      . This also might be true even in hadron collisions at the early lower energies of Intersecting Storage Rings (ISR) and Super Proton Synchrotron (SPS) at CERN      . The other is that the matter created in nucleus or hadron collisions shows a clear feature of collective flow, expanding nearly like a perfect fluid with very low viscosity  -  .
Apart from collective movement, the quanta of produced matter also have the components of thermal motion. The evolution of the produced matter is the convolution of the collective movement and thermal motion. To clarify the role of thermal motion in the expansions of the produced matter in p-p collisions at LHC energies is the major subject of this paper. To this end, we may as usual ignore the minor collective flow in the transverse directions. The transverse movement of the produced matter is therefore only induced by the thermal motion.
The collective movement of produced matter in the longitudinal direction can be solved analytically. There are a number of schemes in dealing with such precise calculations  -  . In this paper, the hydrodynamic model proposed by N. Suzuki is employed  . Besides the analytical characteristics, the other typical feature of this theoretical model is that it incorporates the effects of phase transition into solutions. This coincides with the current experimental observations as mentioned above. Hence, the employed model is more in line with the realistic situations. In addition, the model is related to the initial temperature of QGP, the sound speed in both partonic and hadronic media, the baryochemical potential, as well as the critical temperature of phase transition. This work may therefore help us understand various transport coefficients of expanding system.
In Section 2, a brief introduction is given to the theoretical model  , presenting its analytical solutions. The solutions are then used in Section 3 to formulate the invariant multiplicity distributions of charged particles produced in collisions which are in turn compared with the measurements performed by LHC-CMS Collaboration in p-p collisions at LHC energies of , 2.76 and 7 TeV  , respectively. The last Section 4 is about conclusions.
2. A Brief Introduction to the Model
Here, for the purpose of completion and application, we shall list the key ingredients of the adopted model  .
1) The movement of fluid meets continuity equation
where , t is the time and z is the longitudinal coordinate along beam direction. is the energy-momentum tensor, which, for a perfect fluid, takes the form
where is the metric tensor.
is the 4-velocity of fluid, is its rapidity. and in Equation (2) are the energy density and pressure of fluid, which fit thermodynamic relations
where and are the temperature and entropy density of fluid, respectively. To close Equation (1), another relation, namely the equation of state
is needed, where is the sound speed of fluid, which takes different values in QGP and in hadronic phase.
2) Project Equation (1) to the direction of and the direction perpendicular to , respectively. This leads to equations
The first formula is the continuity equation for entropy conservation. The second one means the existence of a scalar function satisfying relations
From and Legendre transformation, Khalatnikov potential is introduced via
In terms of , the variables and can be expressed as
where is the initial temperature of fluid and . Through above equations, the coordinate base of is transformed to that of , and the first formula of Equation (6) is translated into the so called telegraphy equation
3) Along with the expansions of matter created in collisions, it becomes cooler and cooler. As its temperature drops from the initial to the critical , phase transition occurs. The matter transforms from QGP state to hadronic state. The produced hadrons are initially in the violent and frequent collisions. The major part of these collisions is inelastic. Hence, the abundances of identified hadrons are in changing. Furthermore, the mean free paths of these primary hadrons are very short. The movement of them is still like that of a fluid meeting Equation (10) with only difference being the value of . In QGP, , which is the sound speed of a massless perfect fluid, being the maximum of . In the hadronic state, . At the point of phase transition, that is as , is discontinuous.
4) The solution of Equation (10) for the sector of QGP is 
where is a constant determined by tuning the theoretical results to experimental data. is the 0th order modified Bessel function of the first kind.
In the sector of hadrons, the solution of Equation (10) is 
It is evident that, if , then , , , and thus . At the point of phase transition, , , , . Then
. That is, the potential is discontinuous at point of .
3. The Thermal Motion Induced Transverse Momentum Distributions of Identified Charged Particles
With Khalatnikov potential , the rapidity distributions of fluid read as 
where is the cross area of overlap region of collisions. Inserting Equation (9) into above equation, the part in the round brackets becomes
Along with the expansions of hadronic matter, its temperature continues becoming lower. According to the prescription of Cooper-Frye  , as the temperature drops to the so called chemical freeze-out temperature , the inelastic collisions among hadrons cease. The yields of identified hadrons maintain unchanged becoming the measured results in experiments. Convoluting Equation (15) with the contributions from thermal motion, the invariant multiplicity distributions of charged particles equal   
where is the transverse mass of produced charged particle with rest mass . in Equation (17) is the baryochemical potential. For Fermi charged particles, in the denominator of Equation (17), and for Bosons, . The meaning of Equation (17) is evident. It is the convolution of with the energy of the charged particles in the state with temperature . The former comes from the longitudinal collective flow. The latter stems from the thermal motion.
The right hand side of Equation (17) is evaluated with . At this moment, the fluid freezes out into the charged particles. Replacing in Equation (16) by of Equation (12), it turns into
where is the 1st order modified Bessel function of the first kind.
The integral interval of in Equation (17) is . By using Equations (15), (17), and (18)-(20), together with the definitions in Equation (13), we can get the thermal motion induced transverse momentum distributions of identified charged particles as shown Figure 1. This figure shows the transverse momentum distributions of , and at mid-rapidity of in p-p collisions at , 2.76 and 7 TeV, respectively. The solid dots are the experimental data  . The solid curves are the results of Equation (17). The left and right columns are for the positive and negative charged particles, respectively.
Figure 1. The transverse momentum spectra of identified charged particles produced at mid-rapidity of in p-p collisions at , 2.76 and 7 TeV (from top to bottom). The solid dots are the experimental data  . The solid curves are the results of Equation (17). The left and right columns are for the positive and negative charged particles, respectively.
It can be seen that the theoretical results are in good accordance with experimental data for transverse momentum up to about .
In calculations, in Equation (13) takes the well-recognized value of . in Equation (13) takes the value of from the investigations of Refs.     . Referring to the measurements of Ref.  , the chemical freeze-out temperature takes the value of 160 MeV for in the three different incident energies of 0.9, 2.76 and 7 TeV. For and , takes 170 and 175 MeV, respectively, in the three cases of collisions. The baryochemical potential in Equation (17) takes the value of since in the considered collisions, the abundances of particle and antiparticle are approximately equal to each other  . This can be clearly seen in Figure 1. As for the initial temperature in Equation (13), there is no widely accepted value so far. In this paper, consulting the results presented in Ref.  , takes the values of 3.0, 6.5 and 7.8 GeV for incident energies equaling 0.9, 2.76 and 7 TeV, respectively.
In order to see the importance of thermal motion in the expansions of the matter produced in p-p collisions, the minor transverse collective flow is as usual ignored. The movement of matter in the transverse directions is therefore only induced by the thermal motion. Convoluted it with longitudinal collective movement, we can get the invariant multiplicity distributions of charged particles.
The longitudinal collective flow follows the relativistic hydrodynamics. The most striking feature of the model used in this paper is that it takes the phase transition into consideration, which might have appeared in the current experiments in nuclear or even in hadronic collisions. Hence, the adopted model might reflect a more realistic process of spatiotemporal evolution of matter created in collisions.
The model contains a rich information about the transport coefficients of fluid, such as the sound speed in QGP , the sound speed in hadronic state , the phase transition temperature , the chemical freeze-out temperature , the baryochemical potential and the initial temperature . Exception for , the other five coefficients take the values either from the well-known theoretical results or from experimental measurements. As for , there is no much knowledge about it at present. In our calculations, is estimated referring to other researches.
The present investigations show that the distributions of charged particles measured in p-p collisions at LHC energies of , 2.76 and 7 TeV can be well described by the theoretical model for up to about . Certainly, the transverse collective flow objectively exists even though it is minor. It should affect the motion of the matter. The existing investigations have shown that if the transverse collective flow is taken into account, the describable region is extended to about    . Furthermore, to describe the experimental data in the region of , some other mechanisms, such as the method of blast-wave and perturbative QCD  , should be introduced. Considering that the overwhelming majority of charged particles situates in the region of , we may conclude that the thermal motion plays a non-ignorable role in the transverse expansions of the matter produced in p-p collisions at LHC energies.
This work is supported by the Shanghai Key Lab of Modern Optical System.