# Meson Production in High Energy p+p Collisions at the RHIC Energies

Authors: Bao-Chun Li; Ya-Zhou Wang; Er-Qin Wang
Journal: Advances in High Energy Physics
 (2013)

---

## Abstract

Transverse momentum spectra of mesons produced in p+p collisions are studied in the framework of a thermalized cylinder model. In the region of high transverse momentum, the considered distributions have a tail part at the maximum energy of RHIC. A two-component distribution based upon the improved cylinder model is used to fit the experimental data of the PHENIX Collaboration. It is found that the improved approach can describe the meson production in the wider range of transverse momenta.


---

## Body

## 1. Introduction

Relativistic heavy ion collisions are performed to study nuclear matter under extreme temperature and density conditions [1]. Proton-proton collisions are conventionally used as a reference to compare with nuclear collisions and to understand the observed collective effects. And the Large Hadron Collider (LHC) [2] was originally designed to rise center-of-mass (cm) energy to 14 TeV, which is almost 8 times the record of 1.8 TeV achieved by the Fermilab in the United States. After the LHC data obtained in p+p interactions at sNN=900 GeV and 7 TeV [3, 4], the new interest in general features of p+p collisions at ultrarelativistic energies appeared. In p+p collisions, the meson spectra provide insight into the particle production mechanism and interaction in the hadronic and quark gluon plasma (QGP) phases. Furthermore, the detailed study of meson spectra is important because it acts as an ingredient for estimating the hadronic decay backgrounds in the photon, single lepton, and dilepton spectra, which are the penetrating probes of QGP. To isolate phenomena related to the dense and hot medium created in such collisions, it is also important to measure particle production in smaller collision systems like p-p and d-Au. Measurements of transverse momentum spectra for particles produced in p+p collisions are used as a baseline, to which similar measurements for heavy ion collisions are compared. In addition, the nuclear modification factor RAA of several identified hadrons with high transfer momentum is used to probe jet quenching.The information about the production process is retained by the final-state particle distributions in the collisions [5–9]. Single-hadron production at large transverse momenta in high-energy hadronic and nuclear collisions results from the fragmentation of quarks and gluons issuing from parton-parton scatterings with large momentum transfer. To explain the abundant experimental data, different phenomenological mechanisms of initial coherent multiple interactions and particle transports were proposed and extended in recent years [10–15]. The comparison to model calculations can provide valuable information of the collision evolution and help better understand properties of the QGP. In particular, several theory models of high energy collisions were reported in a workshop held at the CERN Theory Institute [16]. Recently, systematic studies on the production of final-state particles performed by the NA49 collaboration were discussed [17]. Hadronic transport models fail to describe the production of final-state particles, while the results of statistical models are generally in good agreement with the measured particle yields at available energies.Based on the one-dimensional string model [19] and the fireball model [20], we have developed a thermalized cylinder model, which successfully describes the particle production in heavy ion collisions over an energy range from the Alternating Gradient Synchrotron (AGS) to the relativistic heavy ion collider (RHIC) [21–25]. In our previous work [25], the transverse momentum distributions of strange hadrons produced in Cu+Cu and Au+Au collisions at RHIC energies were explained by a single component distribution. The excitation degree of the emission source is allowed to be determined by studying the transverse momentum spectra, anisotropic flow effects, and their correlations. We found that the single component distribution only describes a narrow transverse momentum rang. Recently, the invariant differential cross-section for production of neutral mesons in p+p collisions at sNN=200 GeV were published by the PHENIX Collaboration [18]. The high pT reach of the transverse momentum is helpful to characterize the mechanisms of truly perturbative parton-parton scatterings and parton fragmentation in different QCD environments. Particle yields observed in these experiments inspired our work. It is interesting for us to analyze the results of pp collisions at the RHIC energies. In order to verify the thermalized cylinder model and describe the broad distribution range, in this paper, we will develop the single component spectra into a two-component distribution.


## 2. The Formula

According to the cylinder model [25], we assume the projectile and target cylinders to be formed in A+A, d+A and p+p collisions at high energies. The cylinders are wounded sources of particles. The idea comes from the observation that the process of particle production is not instantaneous, which was noted in [26, 27]. In the reference frame where the longitudinal momentum of a produced particle vanishes, the minimal time necessary for its creation is t0⩾1/MT, where MT=pT2+m2. In the laboratory frame, the particle in question acquires some longitudinal momentum and we have
(1)t⩾γt0=Emt2=coshymt,
where γ and E are the Lorentz factor and the energy of the particle, respectively. So, the resolving power of the longitudinal distance (the uncertainty of the distance from the collision point to those at which the particle is created) is
(2)l=vt=sinhymt.
When the rapidity of the produced particle is large enough so that l>Z(b), where Z(b) is the size of the nucleus at a given impact parameter, the particle cannot resolve individual collisions. It is natural to suggest that its creation may be insensitive to the number of collisions in the source. This is the origin of the idea of wounded sources. The concept of the wounded source may practically be applied in the whole rapidity region.In AA collisions the cylinder is thick, and in p+p collisions the cylinder is thin. Final-state particles are randomly produced from emission sources in the cylinder(s). The excitation degree of the side-surface region is naturally lower than that of the central axis region in the cylinders. From central axis region to side-surface region of the concerned emission source, the excitation degree is assumed to decrease linearly along the transverse axis direction [25]. When the excitation degree increases, the value of the distribution width σ becomes bigger and bigger. So, the excitation degree can be characterized by the momentum distribution width σ. The emission points with the same excitation degree form an emission circle in the transverse momentum space. Therefore, a given σ corresponds to wounded sources which emits a fixed density of particles. Let σs and σc denote the widths of the transverse momentum distribution of particles produced in the side-surface and central axis regions, respectively. The distribution of σ can be given by
(3)f(σ)=C(σc-σ),
with the normalization constant
(4)C=2(σc-σs)2.
According to the single component distribution in the Monte Carlo calculation, the transverse momentum distribution of final-state particles can be given by
(5)pT=(σc-(σc-σs)2-2R1C)-2lnR2,
where R1 and R2 are the random variable distributed in [0,1].In our previous work [25], transverse momentum spectra of strange particles produced in high-energy collisions were investigated by the above thermalized cylinder model. Considering the simply cylinder shape, we obtained the emission source location dependence of the exciting degree specifically. It is shown that the single component distribution, that is, (5), is successful in describing the experimental results measured by the STAR and PHOBOS Collaborations. Moreover, we found that the single component distribution can describe a narrow transverse momentum range. To explain the wider transverse momentum spectra of identified particles produced in p+p collisions, we need to consider the relative importance of hard versus soft processes in the particle production mechanisms at different energies. Hard parton-parton scatterings with large momentum transfer occur on a short time scale and are governed by perturbative QCD. The bulk of particle production occurs via soft processes (with low momentum transfer and consequently longer time scales) which are described with phenomenological models.In order to describe the invariant cross-section of hadrons as a function ofpT over a wide range, Hagedorn proposed an empirical formula [28], which is described as “inspired by QCD” and is given by
(6)Ed3Ndp3=B(1+pT/λ)n≃B{(-npTλ)forpT→0,(λpT)nforpT→∞,
where B, λ, and n are fit parameters. The limiting cases show that the distribution is an exponential form at low transverse momenta, but it is a power law at large transverse momenta, which came from “QCD inspired” quark interchange model [29],
(7)Ed3Ndp3~(mT2)-4~1(pT)8.The dominant contribution to the inclusive cross section is lowpT particles. However, with center of mass energy increasing, the hardening of the pT distribution implies an increase of the transverse momentum. At both low pT and high pT, UA1 collaboration [30] gives a hybrid form,
(8)Ed3Ndp3=Aexp(-bmT),forpT<pTc,=B(1+pT/p0)n,forpT>pTc,
where pTc is a free parameter. In order to give a good description of pion spectra in wide pT range, the PHENIX Collaboration [18, 31, 32] found a single form referred as the modified Hagedorn formula. The modification is to better describe the π0 spectrum for wider pT range, in particular at high pT where the spectrum behaves close to a simple power law function. The single formula has been used to successfully describe the hadron spectra measured in p-p collisions at different energies [18, 31, 32]. The formula has been extensively used in terms of mT,
(9)Ed3Ndp3=B[exp(-amT-bmT2+mT/λ)]n,
where a and b are the fit parameters. The distribution function is close to an exponential form at low pT and a pure power law form at high pT. Similarly, our model should be improved to a two-component distribution,
(10)PT=k(σ1c-(σ1c-σ1s)2-2R1/C1)-2lnR2+(1-k)(σ2c-(σ2c-σ2s)2-2R3/C2)-2lnR4,
where R3 and R4 are the random variable distributed in [0,1], k(1-k) and C1(C2) are the contribution coefficient and normalization constant of the first (second) component respectively, and σ1s(σ2s) and σ1c(σ2c) are the momentum distribution widths of particles produced in side-surface region and central axis region of the first (second) component, respectively. Our calculation shows that C1≠C2. The first and second items in (10) correspond to the contributions of soft production and hard emission, respectively.The basic model assumes the cylinders with the emitted momentum distribution width changing from the side-surface to the central axis. In fact, except for the most central collisions, the active (overlap) region is certainly not cylindrical. We could understand the “cylinder” as a nonperfect cylinder along the beam direction in the momentum space but not a perfect cylinder in the coordinate space. The emission sources with the same excitation degree stay at the same (sub)surface. This situation is similar to the equipotential surface in electromagnetism.


## 3. Comparison with PHENIX Results

The invariant differential cross-sections of neutral mesons produced in p-p collisions at nucleon-nucleon center-of-mass energysNN=200 GeV in various decay modes are presented in Figure 1. The symbols represent the experimental data of the PHENIX Collaboration [18] and the curves are our calculated results by using (10). The χ2 testing provides statistical indication of the most probable value of corresponding parameters. The values of k, σ1c-σ1s, and σ2c-σ2s obtained by fitting the experimental data are given in Table 1 with χ2 per degree of freedom (dof). The maximum value of the observed pT reaches about 13.0 GeV. We see that the calculated results approximately agree with the experimental data for neutral mesons in the region. The values of k, σ1c, and σ2c are taken to be 99.98±0.02, 1.70±0.12 and 4.64±0.35, respectively. Both σ1s, and σ2s values do not change obviously.Table 1
Values of the parametersk, σ1s~σ1c, and σ2s~σ2c obtained by fitting the experimental pT spectra in Figures 1, 2, and 3. The units of σ1s, σ1c, σ2s, and σ2c are GeV/c.
Particle
k (%)
σ
1
s
~
σ
1
c (×0.1)
σ
2
s
~
σ
2
c (×0.1)
χ
2/dof
ω
99.97
±
0.01
0.09
±
0.02
~
1.66
±
0.10
1.28
±
0.10
~
4.55
±
0.24
1.000
K
99.97
±
0.01
0.06
±
0.01
~
1.60
±
0.12
1.22
±
0.10
~
4.50
±
0.21
1.310
η
′
99.98
±
0.01
0.07
±
0.01
~
1.68
±
0.10
1.28
±
0.10
~
4.65
±
0.27
0.472
ϕ
99.99
±
0.00
0.05
±
0.01
~
1.72
±
0.10
1.30
±
0.10
~
4.75
±
0.24
0.407
π
99.96
±
0.01
0.08
±
0.01
~
1.65
±
0.10
1.26
±
0.10
~
4.52
±
0.22
0.210
η
99.97
±
0.01
0.06
±
0.01
~
1.61
±
0.11
1.24
±
0.11
~
4.48
±
0.26
0.240
p
+
p
-
100.00
±
0.01
0.10
±
0.02
~
1.78
±
0.10
—
0.278
J
/
ψ
100.00
±
0.01
0.16
±
0.03
~
2.25
±
0.15
—
0.447
ψ
′
100.00
±
0.01
0.18
±
0.04
~
2.10
±
0.12
—
0.470Figure 1
Invariant differential cross-section of neutral mesons measured in p+p collisions atsNN = 200 GeV in various decay modes. Experimental data measured by the PHENIX Collaboration [18] are shown by the scattered symbols. Our results calculated from the cylinder model are shown by the curves.To compare these results to other particles, we show the invariant differential cross sections of different particles produced in p+p collisions atsNN = 200 GeV in various decay modes in Figure 2. The different symbols are the experimental data of the PHENIX Collaboration [18]. The solid curves are our results calculated by the model. By fitting the experimental data, the parameter values are given in Table 1 with the values of χ2/dof. Similar to Figure 1, the values of σ1s(σ2s) do not change significantly. We see again that the model describes approximately the invariant differential cross-sections of final-state mesons produced in p+p collisions at the highest RHIC energy. Compared with the case of Figure 1, the values of parameters do not change significantly.Figure 2
Invariant differential cross-sections of different particles measured in p+p collisions atsNN = 200 GeV in various decay modes. Experimental data measured by the PHENIX Collaboration [18] are shown by the scattered symbols. Our results calculated from the cylinder model are shown by the curves.In Figure3, we present the invariant cross section at midrapidity for η in p+p collisions at sNN = 200 GeV. The symbols are the experimental data of the PHENIX Collaboration [33, 34]. The solid curves are our results calculated by the model. The dotted and dashed curves are our results corresponding to the contributions of the first and second items in (10). The parameter values and corresponding χ2/dof are given in Table 1. One can see the particle production of soft production and hard emission intuitively in the figure.Figure 3
Invariantη cross section at midrapidity as a function of transverse momentum in p+p at sNN = 200 GeV. The error bars are the statistical uncertainties. Experimental data taken from the PHENIX Collaboration [33, 34] are shown by the scattered symbols. Our results calculated from the cylinder model are shown by the curves.In order to testify the validity of the model, Figure4 shows the invariant π0 cross section, Ed3σ/d3p, in p+p collisions at different RHIC energies. The symbols are the experimental data of the PHENIX Collaboration [35]. The solid curves are our results calculated by the model. The dotted and dashed curves are our results corresponding to the contributions of the first and second items in (10). The values of k, σ1s~σ1c, and σ2s~σ2c are given in Table 2 with the corresponding χ2/dof. The calculated results are in agreement with the experimental data in p-p collisions at the RHIC energies. The k values decrease with the increase of the collision energy, and σs has no obvious change, equals only about 0.60 GeV/c. The values of σ1c (0.99–1.81 GeV/c) and σ2c (3.94–4.83 GeV/c) increase with the collision energy increasing. As discussed in detail in [25], the excitation degree of emission source on the cylinder central axis becomes higher due to the higher collision energy.Table 2
Values of the parametersk, σ1s~σ1c, and σ2s~σ2c obtained by fitting the experimental pT spectra Figure 4. The units of σ1s, σ1c, σ2s, and σ2c are GeV/c.
Particle
Energy (GeV)
k (%)
σ
1
s
~
σ
1
c (×0.1)
σ
2
s
~
σ
2
c (×0.1)
χ
2/dof
π
0
22.4
100.00
±
0.01
0.05
±
0.01
~
1.05
±
0.06
—
0.600
62.4
99.95
±
0.01
0.05
±
0.01
~
1.26
±
0.09
1.10
±
0.07
~
4.14
±
0.20
1.104
200
99.69
±
0.02
0.07
±
0.01
~
1.69
±
0.12
1.29
±
0.10
~
4.58
±
0.25
0.422Figure 4
Invariantπ0 cross-section as a function of transverse momentum in p+p collisions at different RHIC energies. The error bars are the quadratic sum of statistical and systematic uncertainties. Experimental data taken from the PHENIX Collaboration [35] are shown by the scattered symbols. Our results calculated from the cylinder model are shown by the curves.We would like to point out that the calculated results in Figures1–4 are obtained by (10), which is a general representation in the Monte Carlo method. Equations (6)–(9) are only used to explain why the two components need to be considered in (10). The curves in the figures are numerical results from (10). Some of them are not smooth due to low statistics.


## 4. Discussion

The above comparison shows that the improved cylinder model can be used in the description of the meson production in the wider range of transverse momenta. The values ofχ2/dof for all fits are shown in Tables 1 and 2. The maximum value is 1.420, and the minimum value is 0.246. We demonstrate that the calculated results are in good agreement with the available experimental data. In our calculation, k is the most important parameter which determines the distributive slope and range. In most cases, the second item cannot be neglected due to the contribution of the hard process. A large k (≃1) gives a precipitous and narrow distribution, whereas a smaller k (<1) gives a subdued and wide distribution. The second item in (10) is very small (<0.01) and does not contribute to the spectra at the lowest RHIC energy. In p+p collisions for different energies, we see that 1-k increases slightly with the increasing of collision energies. This means that the contribution of hard emission increases slightly. For the asymmetric system of d+Au collisions, it is predicted that the values of k do not change significantly due to the spatial asymmetry of the collision nucleons. Scattering processes at high energy hadron colliders can be classified as either hard or soft. Quantum Chromodynamics (QCD) is the underlying theory for all such processes, but, the approach and level of understanding arevery different for the two cases. For hard processes (e.g., Higgs boson or high pT jet production), the rates and event properties can be predicted by using perturbation theory. In addition, for soft processes (e.g., the total cross section, the underlying event, etc.), the rates and properties are dominated by nonperturbative QCD effects, which are less well understood. For many hard processes, soft interactions are occurring along with the hard interactions and their effects must be understood for comparisons to be made to perturbative predictions. An understanding of the rates and characteristics of predictions for hard processes, both signals and backgrounds, using perturbative QCD (pQCD) is crucial for both the Tevatron and LHC. We are in a position to evaluate the soft and hard contributions to the observed spectra by using k in the statistical model.The particles produced in high-energy nuclear collisions have attracted much attention since people are trying to understand the properties of strongly interacting quark-gluon plasma by studying the possible production mechanisms [36, 37]. The final-state particles are emitted isotropically in the rest frame of emission sources with the different excitation degree in collisions. If we consider that the local emission source has a motion in the transverse direction, there are interactions among the emission sources, then the transverse flow (directed flow and elliptic flow) can be explained by the model. Thermal-statistical models have been successful in describing particle yields in various systems at different energies. The cylinder model is developed from the fireball model, which is suggested in heavy-ion collisions [20]. The excitation degree varies with location in the cylinder. In the previous work [25], we obtained the emission source location dependence of the exciting degree specifically. From central axis to side-surface of the cylinder, the excitation degree of the emission source decreases linearly with the direction of radius. The excitation degree can easily be characterized by the corresponding distribution width, which is written as (4) by using the Monte Carlo calculation. In the model, because of the influence of the emission source temperature, the values of σs are much smaller than those of the σc. It is consist with the conclusion of [15], where the temperature of the fireball decrease linearly with time.In the present work, we improve the method by considering the difference between the soft and hard emission. The parametersσ1c(σ1s) and σ2c(σ2s) are used to describe the excitation degrees of emission sources closing to the side-surface and central axis of the cylinders, respectively. They are in fact reflections of the excitation degree of the soft emission process and the hard emission process, respectively. Our results agree well with the considered distributions of mesons with high transverse momenta produced in nuclear collisions at RHIC energies. At finite temperature, the stronger the collision strength, the larger the excitation degree. Therefore, on the central axis of the cylinder, the interaction between particles is strongest, and the excitation degree is highest. For different collision energies, our results in the improved cylinder model show that the parameters σ1c and σ2c slowly increase with the increasing of energy, and no obvious change can be observed in the values of σs1(σs2). The dependence of parameters on energies renders that intranuclear cascade collisions in the central axis regions play a more important role at higher energy. In the side-surface regions, a high multiplicity at high energy does not contribute to a high excitation degree due to weak intranuclear cascade collisions. In our opinion, high pT hadrons come primarily from the central region due to their large emission angles. And high momentum hadrons come primarily from the side-surface region due to their weak cascade collisions.


## 5. Conclusions

In summary, the transverse momentum distributions of mesons produced in p+p collisions at the RHIC energies have been studied in the framework of the improved cylinder model. The model is successful in description of the high transverse momentum meson production. Based on our phenomenological approach, we evaluated the soft and hard contributions to the observed spectra. At the same time, it can offer some information about soft and hard interactions in the collisions. In our previous work, a rudimentary investigation of the azimuthal anisotropy has been carried basing upon the cylinder model. Combined with the present work, the model can be used to describe uniformly the momentum distributions of final-state particles produced in high energy collisions.



---

*Source: 486476-2013-06-19.xml*