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 was discovered in a region where nucleon–nucleon interaction is kinematically forbidden. This served as the basis for putting forward a hypothesis [
3] about the presence of nuclear density fluctuations in nuclei, called “fluctons”. It was assumed that, in a target nucleus, two or even more nucleons can periodically approach each other at short distances, forming a cluster of several nucleons. In this case, the interaction of an incident proton with such a multinucleon cluster (flucton) makes 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 referred to as cumulative.
In 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 the per-nucleon momentum of the incident nuclei, which explains the name of “cumulative production”. In both cases, experimental studies of this phenomenon in the regions of fragmentation of the target nucleus and the projectile nucleus are limited to the region of rather small transverse momenta (less than 2 GeV) [
6,
7,
8,
9,
10].
A variety of models has been suggested to describe the process of of particle formation in the cumulative region. These 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 in the initial state of the nucleus, while the second [
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,
29], designed for relatively low energies of colliding nuclei compared to the LHC and RHIC colliders as well as high luminosity, has opened 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, as from the contemporary 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 (see [
30]) are expected to be possible in future experiments at FAIR (Darmstadt).
In the present paper, we study the possibility of observing particle production in the new cumulative region of central rapidities and large transverse momenta in Au–Au collisions using the MPD installation at the NICA complex. We estimate the yield of pions and protons in this new cumulative region caused by the process in which a nucleon of one nucleus interacts with a flucton from 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] to describe the cumulative particle production in the fragmentation region of one of these 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 of particle production in the cumulative region in pA collisions,
and the so-called nuclear scaling was experimentally established as follows:
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 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 [
31]).
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 on either the initial energy (
) or 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
and
= 2.64·10
mb/Gev
; thus, the yield of cumulative protons in the region of nuclear fragmentation (region F in
Figure 1 and
Figure 2) is almost 10
times higher compared to that of pions.
As for the dependence of the yields of cumulative particles on the transverse momentum, in [
9,
10], when analyzing experimental data, it was parameterized by the Gaussian dependence:
the parameter
increases with the cumulative number
x and depends on the kind of particle under investigation. These dependencies can be parameterized for pion and proton production as follows [
31]:
where the square of the transverse momentum is measured in GeV
.
This Gaussian approximation (
7), used in [
9,
10], does not seem quite realistic to us, especially in light of the extrapolation of this dependence to the region of sufficiently large transverse momenta (region C in
Figure 1 and
Figure 2). Therefore, in the works [
31,
32] we proposed an alternative parameterization of the dependence of the data [
9,
10] on the transverse momentum by a simple exponent, which usually provides a more adequate description of the data for large transverse momenta:
where
Unfortunately, as shown in [
31,
32], although both approximations provide an adequate description of the experimental data [
9,
10] in the nucleus fragmentation region (region F in
Figure 1 and
Figure 2), they provide 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
Figure 1 and
Figure 2) available for study at the NICA MPD. In these figures, we depict the pseudorapidity acceptance (
) of the MPD detector in the NN collision center of the mass frame. Region F corresponds to the production in the rear hemisphere in the rest frame of a fragmenting nucleus (
).
Comparing
Figure 1 and
Figure 2, it is apparent that in the new cumulative region the transverse momenta of particles increase with the initial energy. This explains why the study of particle production in this region is 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, in order to eliminate the uncertainty described at the end of the previous section and estimate the yield of pions and protons with large transverse momenta in the new cumulative region of mid-rapidities at NICA energies, we use the quark–parton model of the cumulative particle production from a flucton, which was 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 the 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] (the 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 provided 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 provided 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 provides 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, that is 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), to estimate 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. This means 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 exclude the region near zero rapidity, as 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.
Comparing the
Figure 1 and
Figure 2, we expect the results obtained for the initial energy
= 4 GeV to be more reliable than those for
= 8 GeV, as in the latter case it is necessary to extrapolate the cross-section
to significantly higher values of transverse momenta.
Taking into account the definition of
(
1), the multiplicity of particles in acceptance
in a 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 provided by the formulas
Then, we can calculate
by solving the Equation (
2) with (
3) and (
4). For pions, we find the following:
and for protons
Here,
and
are the energy and momentum, respectively, of the incident proton in the center-of-mass system of NN collisions:
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 value of
is provided by (
3) and (
4) for pions and protons, respectively.
These limits of integration over transverse momentum
are of practical interest, as 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. The difficulty of experimentally measuring the predicted effects lies in the fact that a given interval of change in the cumulative number
x at a fixed value of the particle’s rapidity corresponds to a certain interval of its transverse momenta, which depends on both the initial energy and the type of particle being detected (see the yellow region C in
Figure 1 and
Figure 2). In the general case, the boundaries of the transverse momentum interval in which particles with given values of the cumulative number should be detected in the experiment are easily calculated using Formulas (
26) and (
27).
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
and
for two initial energies
4 and 8 GeV. From this table, it can be seen that the transverse momenta at an initial energy of 8 GeV 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
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 dependencies (
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 carried out similar calculations using the Gaussian dependence on the transverse momentum (
7), as used in the original experimental works [
9,
10], and with the exponential dependence (
7) proposed in [
31,
32]. The value
was chosen as the beginning of the cumulative region because the fits used in these works are valid for both pions and protons when starting from this value. For the total cross section of the pAu interaction, we used the value
= 2 barns.
From
Table 2 it can be seen that, at an initial energy of 4 GeV, the predictions for particle multiplicities obtained within the framework of the described theoretical model are consistent to within approximately an order of magnitude with the results obtained using both Gaussian and exponential fits of the experimental data. At an initial energy of 8 GeV, the results of the present theoretical calculations support the results obtained using the more natural exponential fits for the transverse momentum dependence proposed in [
31,
32], 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 now try to make a rough estimate of the production of cumulative pions and protons for the symmetric AuAu reaction in the region
that is available for study with NICA MPD.
Recall that we have excluded the region near zero rapidity, as 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 [
30,
33,
34]. This process, which is of great physical interest, requires a separate special study, which we leave for our future studies. We only note that such a process can be studied experimentally only in the new cumulative region of central rapidities and large transverse momenta that is available for research using the MPD and SPD facilities of the NICA collider [
26,
27,
28,
29], 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, as they will interact with the flucton in the gold nucleus. Of course, it is necessary to take into account that in this case there is 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. From the ALICE experiment [
35,
36,
37] we know that, for pPb and PbPb collisions at LHC energies,
= 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 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), 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 a value of
equal to
= 15 for the case of min. bias collisions.
The obtained estimates of integral multiplicities and yields of pions and protons in min. bias AuAu collisions at the 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).
The estimates of cumulative particle yields
in
Table 3 are presented for one hour (
t) of operation of the NICA collider. When performing them, we took 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 barns, which corresponds to the values in [
26].
When analyzing the yields of cumulative particles presented in
Table 3, it is necessary to keep in mind that they refer to the case of an “ideal” detector, and provide only an upper limit on the number of detected particles. Despite the fairly high overall efficiency of detecting charged particles at the NICA MPD (about 80–90% [
28]), 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 significantly less (about an order of magnitude) than that provided by general estimates.
From the
values in
Table 3, it is apparent that the increase in the NICA collider luminosity with increasing initial energy from 4 to 8 GeV practically compensates for 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 the new cumulative region of central rapidities and high transverse momenta in Au–Au collisions using the MPD installation of the NICA complex at initial energies of both 4 and 8 GeV.
Conducting research at higher energies of the NICA collider, for example, at
= 11 GeV, is more difficult, as 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,
28].
To study the dependence of the production of pions and protons on the cumulative number within the described theoretical approach, we calculated the inclusive cross-section
characterizing the distribution of cumulative particles in
x. The results of the calculations are presented in the
Figure 5. Note that the errors presented in the figure reflect only the uncertainties arising within the framework of the model.
From
Table 3 and
Figure 5, it can be seen that when the energy increases, 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 yields to pion yields is about 10 for the same value of
x, while at 8 GeV the pion yield already dominates the proton yield more than ten-fold. Recall that in the traditional cumulative region of nucleus fragmentation at low values of transverse momenta, the ratio of proton and pion yields is about
for the same value of
x (see the values of the constants
and
in the paragraph after Formula (
6)).
This effect occurs due to different mechanisms of particle formation with momentum outside the pp-kinematics, that is, 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 have been some experimental indications 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 40
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 [
38] that in this region the ratio p/
for the same value of the cumulative number
x is on the order of 100, which is significantly less than 10
in the 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,
29] using existing and new ultra-thin pixel detector systems. It is important to note 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 simultaneously obtain a signal from several different types of detectors used by the MPD installation [
28]. In this regard, information from the internal tracking system [
39] 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 thereby isolate it from the inevitable background noise.
6. Conclusions
A consequence of the presence of nuclear density fluctuations in colliding nuclei (so-called fluctons [
3]) is the production of particles with momentum in the region outside the nucleon–nucleon kinematics, called cumulative [
4]. From the contemporary point of view, these multi-nucleon fluctons that occasionally appear 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 caused by the interaction of one of the nucleons of a nucleus with a flucton in another nucleus within the new cumulative region of central rapidities and large transverse momenta now available for study at the MPD and SPD facilities of the NICA collider [
26,
27,
28,
29].
Calculations were carried out using a previously developed microscopic approach [
16,
17,
18,
19,
20,
21], allowing us to describe the dependence of particle yields on transverse momentum at different values of the cumulative number simultaneously for both pions and protons using a single parameter, namely, the mass of the constituent quark
= 310 MeV.
It is shown that the values of pion and proton yields found in this new cumulative region indicate 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 predict that in this region of initial energies, the dominance of proton yields over pions is replaced by dominance of pion yields when changing from 4 to 8 GeV. This effect arises due to different mechanisms of particle formation in the interaction of a nucleon with a flucton within the framework of the 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 drawing 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; the rarer process of flucton–flucton scattering [
30,
33,
34], which was not taken into account, can have a significant contribution, especially in the region of the most central rapidities
. The process of flucton interaction 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.