1. Introduction
Long ago, in experiments [
1,
2] on the interaction of high-energy protons with nuclei, an unexpectedly large yield of particles into the rear hemisphere, in a region kinematically forbidden for nucleon-nucleon interaction, was discovered. This served as the basis for putting forward the hypothesis [
3] about the presence of nuclear density fluctuations in nuclei, called “fluctons”. It was assumed that in a target nucleus periodically two or even more nucleons can approach each other at a short distance, forming a cluster of several nucleons. In this case, the interaction of an incident proton with such a multinucleon cluster (flucton) made it possible to explain the observed particle production in the backward direction - outside the region permitted by nucleon-nucleon kinematics. Later this area of momenta was called cumulative.
Over time, when the first beams of relativistic nuclei were obtained, it became possible to study cumulative production in the region of projectile nucleus fragmentation [
4,
5]. In this case, it corresponds to the production of particles with a longitudinal momentum greater than that per nucleon of the incident nucleus, what explains the name “cumulative production”. In both cases, in the regions of fragmentation of the target nucleus and the projectile nucleus, experimental studies of this phenomenon are limited to the region of rather small transverse momenta (less than 2 GeV) [
6,
7,
8,
9,
10].
It was suggested a variety of models to describe the process of the formation of particles in the cumulative region. They can be conditionally divided into two large groups. The first group [
11,
12,
13,
14,
15,
16,
17,
18,
19,
20,
21] assumes the presence of fluctons already in the initial state of the nucleus, and the second group of models [
22,
23,
24,
25] suggests that the dense nuclear matter clusters are formed later, in the process of a nuclear collision.
Currently, the construction of the NICA collider at JINR in Dubna [
26,
27,
28], designed for relatively low energies of colliding nuclei compared to the LHC and RHIC colliders and high luminosity, opens up the possibility of studying the production of particles in a new cumulative region of central rapidities and large transverse momenta. These investigations are of great interest, since from the modern point of view, multinucleon fluctons in nuclei are clumps of dense cold baryon-enriched quark-gluon matter. Studies of the clusters of dense cold quark-gluon matter intrinsic to the nuclei are also possible (see [
29]) in future experiments at FAIR (Darmstadt).
In the present paper, we study the possibility of observing the particle production in the new cumulative region of central rapidities and large transverse momenta in Au-Au collision using MPD installation of the NICA complex. We estimate the yield of pions and protons in this new cumulative region due to the process of interaction of a nucleon of one nucleus with a flucton of another. To describe the dependence of cumulative particles on transverse momentum, we use a microscopic (at the quark level) approach developed earlier [
16,
17,
18,
19,
20,
21] for description of the cumulative particle production in the fragmentation region of one of the colliding nuclei.
2. Features of cumulative production in the nucleus fragmentation region
Most of the experimental data on cumulative production were obtained in the rest frame of a fragmenting nucleus. In this case, when incident protons are scattered on a fixed target nucleus, cumulative particles are emitted into the rear hemisphere, which is convenient for the experimental study of this process.
For the inclusive cross section for particle production in the cumulative region in pA collisions:
the so-called nuclear scaling was experimentally established:
As a scaling variable
x, it was proposed to use the so-called cumulative number
x, defined as the minimum target mass (measured in nucleon masses) that allows the production of a cumulative particle c with a given momentum
k. This variable is suitable because for an integer
it determines the kinematic boundaries for the production of a particle c with momentum
k when an incident proton collides with a flucton consisting of
x nucleons (see
Figure 1 and
Figure 2).
The variable
x is also relativistically invariant and can be calculated in an arbitrary frame of reference from the equation:
Here
is the 4-momentum of the incident proton and
is the 4-momentum per nucleon in nucleus A. The
is equal to
for production of pions and
for production of protons. Here
is the mass of the produced particle c,
m is the nucleon mass and
(for more details see [
34]).
For pA interaction at relativistic energies, as follows from the data (see, for example, the analysis in the articles [
6,
8,
9,
10]), the function
in (
1) does not depend either on the initial energy (
) or on the atomic number of the nucleus A (at least , for heavy nuclei). It can be presented in the following form:
where
The value of the parameter
in (
5) is 0.139 for pions and 0.135 for protons, and the exponential dependence is well satisfied at
for pions and
for protons [
10]. Please note that for the production of pions and protons
=4.0·10
2 and
=2.64·10
6 mb/Gev
2, so, in the region of nuclear fragmentation (region F in
Figure 1 and 2) the yield of cumulative protons is almost 10
4 times higher compared to pions.
As for the dependence of the yields of cumulative particles on the transverse momentum, in the works [
9,
10], when analyzing experimental data, it was parameterized by the Gaussian dependence:
and the parameter
increases with the cumulative number
x and depends on the sort of the particle. These dependencies can be parameterized for pion and proton production as follows [
32,
34]:
where the square of the transverse momentum is measured in
.
This Gaussian approximation (
7), used in [
9,
10], does not seem quite realistic to us, especially in the light of extrapolation of this dependence to the region of sufficiently large transverse momenta (region C in the
Figure 1 and
Figure 2). Therefore, in the works [
32,
34] we proposed an alternative parameterization of the dependence of the data [
9,
10] on the transverse momentum by a simple exponent, which usually gives a more adequate description of the data for large transverse momenta:
where
Unfortunately, as shown in [
32,
34], although both approximations provide an adequate description of the experimental data [
9,
10] in the nucleus fragmentation region (region F in the
Figure 1 and
Figure 2) they give very different results, especially at higher initial energy (
=8 GeV), when we use them to calculate particle yields in the new cumulative region of central rapidities and large transverse momenta (region C in the
Figure 1 and
Figure 2) available for study at the NICA MPD. In these figures we depictured the pseudorapidity acceptance (
) of the MPD detector in the NN collision center of mass frame. The region F corresponds to the production in the rear hemisphere in the rest frame of a fragmenting nucleus (
).
Comparing
Figure 1 and
Figure 2, we also see that in the new cumulative region the transverse momenta of particles increase with the initial energy. This explaines why the study of particle production in this region are absolutely impossible at the LHC and RHIC energies.
3.
Theoretical description of the dependence on transverse momentum in the cumulative region
In the present paper, to eliminate the uncertainty described at the end of the previous section and to estimate the yield of pions and protons with large transverse momenta in the new cumulative region of mid-rapidities at the NICA energies, we use the quark-parton model of the cumulative particle production from a flucton, developed earlier [
16,
17,
18,
19,
20,
21].
As shown in this model, the formation of cumulative pions and protons is dominated by two different mechanisms. In the case of pion production, the fragmentation of one flucton quark into a pion predominates [
16,
17,
18] (see left panel in
Figure 3), whereas in the case of proton production, the mechanism of coherent coalescence (recombination) of three flucton quarks into a proton is dominant [
19,
20,
21] (right panel in
Figure 3).
As shown in [
21] in the framework of this approach the dependence on transverse momentum for the production of cumulative pions is given by the expression:
where
is the number of donor quarks transferring their momentum to the active quark forming a cumulative pion (
denotes the integer part) and
Here
is the Bessel function, and
is the modified Bessel function (McDonald function). Note that for
the integral (
12) is calculated explicitly:
The dependence on transverse momentum for the production of cumulative protons is given by the expression:
where for the considered interval of
x for protons (
):
where
is a step function. The total number of donor quarks transferring their momentum to three active quarks forming a cumulative proton is equal to
.
Note that
, which is the constituent quark mass, is the only parameter in these formulas (see (
11) and (
14)). With a natural value of this parameter
= 310 MeV, we obtain a simultaneous description of the dependence of the production of cumulative pions and protons on the transverse momentum (see
Figure 4).
Really in
Figure 4 we see that with the natural value of a single parameter
this approach correctly describes the broadening of transverse momentum distributions with increasing cumulative number, which is observed for both pion and proton yields. Moreover, at the same time, it gives correct wider transverse momentum distributions for pions compared to protons for the same value of the cumulative variable. In this approach, it arises due to different mechanisms of the formation of particles with momenta outside the pp-kinematics - fragmentation of one flucton quark for a pion (left panel in
Figure 3) and coherent coalescence (recombination) of three flucton-quarks for a proton (right panel in
Figure 3).
4. Yield of cumulative particles at mid-rapidities in pA colllisions
In this section, we estimate the yield of cumulative pions and protons with high transverse momentum at central rapidities in pA collisions, using the theoretical description of the dependence of cumulative production on transverse momentum obtained in the previous section. The idea is to use the description of the inclusive cross section of the cumulative particle production
(
1) obtained in the nucleus fragmentation region (region F in
Figure 1 and
Figure 2), for to estimate the particle yields in the new cumulative region of central rapidities (region C in
Figure 1 and
Figure 2), available for study at the NICA collider.
For region C (see
Figure 1 and
Figure 2), the rapidity interval
is selected. We mean that in the center of mass frame the incident proton moves along the
z axis, and the fragmenting nucleus moves in the opposite direction. We also exclude the region near zero rapidity, since the theoretical approach [
16,
17,
18,
19,
20,
21] (
Figure 3), used in the previous section to describe the transverse momentum dependence of cumulative production, is not valid at small
values.
Compairing the
Figure 1 and
Figure 2, we expect that the results obtained for the initial energy
=4 GeV will be more reliable than for
=8 GeV, since in the latter case it is necessary to continue the cross section
to significantly higher values of transverse momenta.
Taking into account the definition of
(
1) the multiplicity of particles in acceptance
in pA collision is determined by the expression:
Using the relativistic invariance of
we can write (
18) in the center-of-mass system of NN collision and move from
to rapidity
y
where
According to (
20), for given
y and
, the values of
and
are given by the formulas:
Then we can calculate
solving the equation (
2) with (
3) and (
4). We find for pions:
and for protons:
Here
and
are the energy and momentum of the incident proton in the center-of-mass system of NN collision:
In our calculations, we first performed integration over the transverse momentum
, and then over the rapidity
y:
The limits of integration over the transverse momentum
were specified by setting the minimum and maximum values of the cumulative number under consideration:
where
is defined by the formulas
The
is given by (
3) and (
4) for pions and for protons.
These limits of integration over transverse momentum
are of practical interest, since they determine the interval of transverse momenta
for a given rapidity
y, in which particles with given values of the cumulative number must be registered in the experiment. In
Table 1 we present them for the production of cumulative pions and protons for two values of the cumulative number
x=1.6 and 3.0 and values of rapidity
-1 and -0.5 for two initial energies
4 and 8 GeV. From this table we see that at an initial energy of 8 GeV the transverse momenta are approximately two times greater than at 4 GeV (see
Figure 1 and
Figure 2).
The results of our calculations of pion and proton multiplicities using the formula (
25) in the new cumulative region at
and
(region C in the
Figure 1 and
Figure 2) for pAu collisions are presented in
Table 2. For the dependence of particle yields in the cumulative region on transverse momentum, we used dependences (
11) and (
14), obtained within the framework of the theoretical approach [
16,
17,
18,
19,
20,
21] presented in the previous section.
For comparison, we also carried out similar calculations with the Gaussian dependence on transverse momentum (
7), used in the original experimental works [
9,
10], and with the exponential dependence (
7), proposed in works [
32,
34]. The value
was chosen as the beginning of the cumulative region, since, starting from this value, the fits used in these works are valid for both pions and protons. For the total cross section of the pAu interaction we used the value
=2 bn.
From
Table 2 we see that at an initial energy of 4 GeV, the predictions for particle multiplicities obtained within the framework of the described theoretical model are approximately in order of magnitude consistent with the results obtained using both Gaussian and exponential fits of experimental data. At an initial energy of 8 GeV, the results of present theoretical calculations support the results obtained using the more natural exponential fits for the transverse momentum dependence proposed in [
32,
34], while the use of Gaussian type fits [
9,
10] predicts extremely low particle multiplicities into this new cumulative region at this energy.
The reason for this is, of course, that at an initial energy of 8 GeV the transverse momenta in the new cumulative region C (see
Figure 1 and
Figure 2 and
Table 1) are approximately two times greater than at 4 GeV .
5.
Estimates of cumulative production in the region available for study with NICA MPD.
Using the estimates of the cumulative production in pAu collisions in the rapidity region
(region C in
Figure 1 and
Figure 2), made in the previous section, we will now try to make a rough estimate of the production of cumulative pions and protons for the symmetric AuAu reaction in the region
available for study with NICA MPD.
Recall that we have excluded the region near zero rapidity, since the theoretical approach [
16,
17,
18,
19,
20,
21] (
Figure 3), used in the present work to describe the transverse momentum dependence of cumulative production (see
Section 3), is not valid at small
values. Another reason for excluding this region from the present consideration is the fact that in this region it may be important to take into account the contribution of the rarer flucton-flucton scattering process [
29,
30,
31]. This process, which is of great physical interest, requires a separate special study, which we leave for our future studies. Let us only note that such a process can be studied experimentally only in the new cumulative region of central rapidities and large transverse momenta, available for research using MPD and SPD facilities of NICA collider [
26,
27,
28], and cannot be studied in the traditional cumulative region of fragmentation of one of the nuclei.
In order to obtain estimates for cumulative particles in the region in the AuAu reaction, based on estimates of their yields in this region in the pAu collision, it is necessary to take into account the increased effective flux of nucleons compared to protons, which will interact with the flucton in the gold nucleus. Of course, it is also necessary to take into account that in this case there is also a symmetrical contribution to the rapidity region , which comes from the interaction of the nucleons of the second nucleus with the flucton of the first nucleus.
Replacing an incident proton with a nucleus increases the number of projectile nucleons interacting with a flucton in another nucleus. To take this into account, we introduce an effective factor . The magnitude of this factor can be estimated through the ratio of the number of participating nucleons or NN collisions in pAu and AuAu reactions. Clearly, the result will strongly depend on the centrality of the AuAu collision.
At high energies, the values of
and
are significantly different. So from the ALICE experiment [
36,
37,
38] we know that at LHC energies for pPb and PbPb collisions
=15 for min.bias events and increases to a value of 24 for 0-5% of the most central events, while
= 55 for min.bias of events, increasing to 115 for central ones. As is known, the number of participating nucleons,
, is determined mainly by the collision geometry (the value of the impact parameter [
39]) and weakly depends on the initial energy, while the number of nucleon-nucleon collisions,
, increases significantly with increasing initial energy, therefore, taking into account the relatively low energies of the NICA collider, for further rough estimates we chose the value of
equal to
=15 for the case of min.bias collisions.
The obtained etimates of integral multiplicities and yields of pions and protons in min.bias AuAu collisions at NICA collider in the cumulative region
and
due to the process of interaction of a nucleon of one nucleus with a flucton of another are presented in
Table 3. The contribution of nucleon-flucton interactions was calculated using the dependences on the transverse momentum (
11) and (
14) obtained in the theoretical approach [
16,
17,
18,
19,
20,
21] (see
Section 3).
Estimates of cumulative particle yields,
in
Table 3 are presented for one hour (
t) of operation of the NICA collider. When performing them, we take into account that the design luminosity of the NICA collider for AuAu collisions at an energy of 8 GeV will be 100 times higher than at an energy of 4 GeV -
=
and
respectively, resulting in interaction rates,
of 7 kHz and 70 Hz with
7 bn [
26,
27]. When analyzing the yields of cumulative particles presented in
Table 3, it is necessary also to keep in mind that in a real experiment the final number of “good” events that have been selected according to various criteria, such as the position of the interaction vertex, the activation of various triggers, etc., usually turns out to be significant (approximately an order of magnitude) less than that given by general estimates.
From the
values in
Table 3 we see that the increase in the NICA collider luminosity with increasing initial energy from 4 to 8 GeV practically compensates the overall drop of the integral multiplicities
, arising due to the general increase of transverse momenta in the cumulative region with energy. This allows us to conclude that it is possible to observe the production of particles in a new cumulative region of central rapidities and high transverse momenta in Au-Au collisions using the MPD installation of the NICA complex both at initial energies of 4 and 8 GeV.
Conducting research at higher energies of the NICA collider (for example, at
= 11 GeV) will be more difficult, because a further drop in integral multiplicities
will no longer be compensated by an increase in luminosity, which remains at the same level at 11 GeV as at 8 GeV [
26,
27].
To study the dependence of the production of pions and protons on the cumulative number within the described theoretical approach, we also calculated the inclusive cross section
characterizing the distribution of cumulative particles in
x. The results of the calculations are presented in the
Figure 5.
From
Table 3 and
Figure 5 we see that with increasing initial energy, the proton yield decreases much faster than the pion yield for the same fixed value of the cumulative number
x. At 4 GeV the ratio of proton to pion yields is about 10 for the same value of
x, but at 8 GeV the pion yield already dominates the proton yield by more than 10 times. Let us recall that in the traditional cumulative region of nucleus fragmentation at low values of transverse momenta, the ratio of proton and pion yields was about
for the same value of
x (see the values of the constants
and
in the paragraph after the formula (
6)).
This effect occurs due to different mechanisms of the formation of particles with momentum outside the pp-kinematics - coherent coalescence (recombination) of three flucton quarks for a proton and fragmentation of one flucton quark for a pion in the theoretical approach [
16,
17,
18,
19,
20,
21] described in
Section 3.
There is some experimental indication that such an effect does occur. Results of the SPIN collaboration on the production of protons and pions with large transverse momenta at an angle of
on stationary nuclear targets by protons with an energy of 50 GeV/c (
= 9.8), which corresponds to the cumulative number
x up to 1.2 show [
40] that in this region the ratio p/
for the same value of the cumulative number
x is of the order of 100, which is significantly less than
, which was the case in traditional cumulative region of nucleus fragmentation.
It is assumed that the obtained dependences will be studied in MPD and SPD experiments at the NICA collider [
26,
27,
28] using existing and new ultrathin pixel detector systems. It is important that for reliable registration of very rare events of particle creation in the cumulative region and reliable separation of their tracks from various kinds of false background tracks, it is necessary to have a signal simultaneously from several types of detectors used by the installation. In this regard, information from the internal tracking system is especially important, making it possible to reliably confirm the exit of the track of a cumulative particle from the vertex of the primary interaction, and thus isolate it from the inevitable noise background.
6. Conclusion
A consequence of the presence of nuclear density fluctuations in colliding nuclei - the so-called fluctons [
3] - is the production of particles with momentum in the region outside the nucleon-nucleon kinematics, called cumulative [
4]. From the modern point of view, these multinucleon fluctons occasionally appearing in nuclei are clumps of dense cold baryon-enriched quark-gluon matter. Therefore, studying the process of nucleon scattering on such nuclear density fluctuations with the production of a particle in the cumulative region is of great interest.
In this work, we estimate the yields of pions and protons, due to the interaction of one of the nucleons of a nucleus with a flucton in another nucleus, in the new cumulative region of central rapidities and large transverse momenta, available for study at the MPD and SPD facilities of the NICA collider [
26,
27,
28].
Calculations were carried out using a previously developed microscopic approach [
16,
17,
18,
19,
20,
21], which allows us to describe dependence of particle yields on transverse momentum at different values of the cumulative number simultaneously for both pions and protons using a single parameter - the mass of the constituent quark
= 310 MeV.
It is shown that the found value of pion and proton yields in this new cumulative region indicates the possibility of studying this phenomenon in collisions of heavy nuclei at the MPD facility of the NICA complex at low initial energies of 4 - 8 GeV.
Theoretical calculations also predict that in this region of initial energies, when going from 4 to 8 GeV, the dominance of proton yields over pions is replaced by dominance of pion yields. This effect arises due to different mechanisms of particle formation in the interaction of a nucleon with a flucton within the framework of the used theoretical approach [
16,
17,
18,
19,
20,
21] - fragmentation of one flucton quark for pion production and coherent fusion (recombination) of three flucton quarks for proton production (see
Section 3).
However, when making these conclusions, it is necessary to keep in mind that in this work only the contribution of the interaction of a nucleon with a flucton was taken into account and the rarer process of flucton-flucton scattering [
29,
30,
39] was not taken into account, which, however, can make a significant contribution especially in the region of the most central rapidities
. The process of interaction of fluctons is of great physical interest and requires a separate special study, which we leave for our future research. It is important that this process can be studied experimentally only in the new cumulative region of central rapidities and large transverse momenta, available for research at the MPD and SPD facilities of the NICA collider, and cannot be studied in the traditional cumulative region of fragmentation of one of the nuclei.