Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Numerical Simulations of the Decaying Transverse Oscillations in the Cool Jet
Next Article in Special Issue
Probing Gluons with the Future Spin Physics Detector
Previous Article in Journal
Selected Topics of Social Physics: Equilibrium Systems
Previous Article in Special Issue
Photons as a Signal of Deconfinement in Hadronic Matter under Extreme Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Interacting Colour Strings Approach in Modelling of Rapidity Correlations

Laboratory of Ultra-High Energy Physics, St. Petersburg University, 7/9 Universitetskaya Nab., St. Petersburg 199034, Russia
*
Author to whom correspondence should be addressed.
Physics 2023, 5(2), 636-654; https://doi.org/10.3390/physics5020042
Submission received: 28 February 2023 / Revised: 14 May 2023 / Accepted: 18 May 2023 / Published: 20 June 2023
(This article belongs to the Special Issue From Heavy Ions to Astroparticle Physics)

Abstract

:
In this paper, using the concept of multi-pomeron exchange, we develope a Monte Carlo model of interacting quark–gluon strings acting as particle-emitting sources aimed at describing inelastic proton–proton interactions at high energies. The implemented 3D (three-dimensional) dynamics of colour string formation resulted in their finite length in the rapidity space and in the fluctuating event-by-event spatial density. Thus, this results in string cluster formation because of the fusion mechanism and the appearance of long-range multiplicity and mean transverse momentum (mean- p T ) correlations in rapidity. We study, via the pseudorapidity dependence, the sensitivity to the details of the 3D dynamical formation of strings for several observables such as the forward–backward correlation coefficient value, strongly intensive quantity, Σ , and the “almost” strongly intensive observable, the variance, σ C 2 , of the distribution of the asymmetry coefficient, C. The strongly intensive quantity Σ is used in this study to suppress trivial statistical fluctuations in the number of particles emitting similar types of sources and to reveal the intrinsic fluctuations of a single source. We demonstrate the connection between Σ and such often used observables as cumulants, factorial cumulants, and σ C 2 . We stress the importance of the contribution of “short” strings and the event asymmetry of the initial conditions on the long-range correlation measures. We argue that string cluster formation because of the fusion mechanism explains the collective effects seen in multiplicity and transverse momentum–multiplicity, p T N, long-range correlation functions.

1. Introduction

Living in the era of the LHC (Large Hadron Collider), one has access to an unprecedentedly huge amount of data recorded by the state-of-the-art experimental facilities. Many intriguing results have been obtained for both proton–proton and ion–ion collisions that require a deep perception from the theoretical point of view. For example, the two most-significant ones are (i) the appearance of a so-called ridge phenomena in proton-proton (pp) collisions at the LHC—the long-range rapidity structure [1]—and (ii) the onset of collectivity in small systems [2].
Since the calculations from the first principles are often not possible in the theory of strong interactions and in the QCD (quantum chromodynamics), then the most-popular way to understand the experimental results is to compare them to simulations performed with some widely used Monte Carlo event generators. They usually contain many versions and tunes, including various options for some physical processes to be selected to be switched on and off.
In this study, we probe the sensitivity of different experimental observables to the details of the particle source formation in 3D (three-dimensional) dynamics to see a particular effect of one or another model’s features on the final correlation measures.
We focus on the modelling of initial-state configurations of pp collisions, and we attempt to mimic soft processes that result in multiparticle production. To implement this, we follow the philosophy of the phenomenological colour string model that was initially based on the two-stage particle production scenario [3]. In the first stage, the longitudinally extended tubes of colour flux, so-called quark–gluon strings, are formed between wounded partons of the colliding hadrons. In the second stage, strings fragment into hadrons via the colour vacuum break-down because of the creation of quark–antiquark pairs. Therefore, in this framework, one naturally accounts for the initial conditions of hadronic collision via the quark–gluon string kinematics. The effective string breaking results later in the production of observed particles.
What in principle differentiates our Monte Carlo model from the aforementioned event generators and from the original string model approach is the mechanism of the string–string interaction taken in the form of the so-called string fusion [4]. In the case of the absence of any interactions, the factorisation occurs, and the resulting spectra can be considered as a convolution of the spectra originating from the independent strings and of the distribution in several stretched strings. In the case of a high density of strings formed, the interactions between the strings as a fusion will increase the colour field density inside the new colour flux tubes formed. This changes the strings’ characteristics, which affects the particle production.
In this paper, we argue that not only string fusion plays an important role at the initial stages, but another peculiarity, the consideration of the detailed dynamics of strings’ formation is essential. The latter complicates the inner structure of the model and deprives it of the translational invariance in the rapidity. Thus, the model is a powerful tool for studying long-range rapidity correlations, and in combination with the string fusion effect, it may serve to study different types of particle production sources.
The main object is to quantify the long-range correlations between some observables measured event-by-event in pp collisions in two separate pseudorapidity intervals that are usually selected symmetrically with respect to the midrapidity. Based on the causality principle, it was shown in Ref. [5] that, if the long-range rapidity correlation between particles exists, then the correlation must be formed at a proper and very early stage of collision. Therefore, in this study, we concentrate on the long-range rapidity correlations, as we are most interested in the role of the initial conditions and the sensitivity to the details of the evolution of the colliding system.
In this paper, the event-by-event fluctuations in the number of particle-emitting sources are for brevity called “volume” fluctuations. This reflects our approach to the study of the initial states of a system that has a complex spatial structure. We consider the so-called strongly intensive observable, Σ , which suppresses these trivial fluctuations. In addition, we look at the event asymmetry coefficient that reflects the asymmetrical distribution of the particle sources (strings and string clusters) in rapidity, which results in the forward–backward correlations as well.
In the present study, the resonance decays, jets, and other particle production mechanisms characterised by short-range correlations are not considered. This is because of the main interest in the search for observables being related to the string fusion phenomena and the appearance of long-range rapidity correlations measured in hadron collisions in well-separated forward and backward rapidity intervals (about a unit of rapidity). Moreover, as soon as we study soft processes of particle production, i.e., non-perturbative effects, the main role is played by the hadrons with the transverse momentum, p T , up to 1.5–2 GeV/ c (with c the speed of light), and our model does not include the simulation of hard processes.
The paper is organized as follows. In Section 2, we briefly present the developed model and describe its key features and parameters. In this Section, we introduce the longitudinal and transverse dynamics of strings and describe the mechanism of their fusion. In Section 3, we perform the fits to the global observables and spectra, i.e., we explain how the parameters listed in Section 2 are obtained. Further, in Section 3, we present the model predictions for the forward–backward rapidity correlations measured with the correlation coefficient, the strongly intensive quantity, Σ , and the “almost” strongly intensive observable, σ 2 ( C ) , where C is the event-by-event asymmetry coefficient of the multiplicity distribution in rapidity. In addition, we show the connection of Σ with such quantities like the often used cumulants and factorial cumulants. We also obtain the dependence of Σ on the selected forward-multiplicity-based classes. Finally, in Section 4, we discuss the results obtained.

2. Definitions and Model

In the model considered here, we assume an inelastic collision of two protons (called an event).

2.1. Multi-Pomeron Exchange in pp Collisions

The number of strings, n str , formed in the pp-collision event comes from the pomeron number, n pom , distribution [6] as n str = 2 n pom . This is due to the mechanism, in which colour strings appear from the unitarity cut of a cylindrical pomeron exchange diagram [7], while without the cut, the diagram represents an elastic scattering. We sample n pom from
P n pom = C z 1 z n pom 1 exp z l = 0 n pom 1 z l l ! ,
where z = 2 w γ s Δ R 2 + α ln s , w = 1.5 is the quasi-eikonal parameter related to the small-mass diffraction dissociation of incoming hadrons, Δ = α ( 0 ) 1 = 0.2 is the residue of the pomeron trajectory, γ = 1.035 GeV 2 and R 2 = 3.3 GeV 2 characterise the coupling of the pomeron trajectory with the initial hadrons, the slope of the pomeron trajectory α = 0.05 GeV 2 , s is the square of the collision centre-of-mass energy, and C ( z ) is a normalising coefficient.
The values of these parameters are taken from Ref. [8], where one assumes that the number of primary strings in an event (related to n pom from Equation (1)) is initially formed, and then, some of them can fuse, forming the string clusters with the modified fragmentation characteristics. This framework can be viewed as a three-step scenario: first, strings are produced; then, they interact (fusion); lastly, the strings hadronisation occurs.
One may note that this is in contrast to another approach used in Ref. [9], where the strings were considered to be produced of different types from the very beginning, so that the string fusion was effectively implemented already at the string formation moment. These two methods will necessarily differ in the values of the parameters in Equation (1) so that the final charged particle multiplicity distribution, P N ch (2), describes the data. However, the present scenario allows us to introduce the 3D dynamics of string formation and evolution. Formally, in our approach,
P N ch = n pom = 1 P n pom P n pom N ch ,
although the exact form of P n pom N ch is unknown because of the complexity of an event picture in our model of interacting strings of different lengths and positions in the rapidity space; thus, we treat P n pom N ch within a Monte Carlo simulation.
The mechanism of pp collision is described as a multi-pomeron exchange, and we consider only strings stretched between quarks and anti-quarks of the colliding protons (no diquarks as string end-points were implemented). Taking this into account, we suppose that the number of partons in each of the colliding protons is equal to the number of strings per event.

2.2. Event-by-Event Sets of Protons and Parton’s Permutations

The problem to solve is that for every collision event, characterised by a certain number of exchanged pomerons, one has to form two protons with a given number of partons. Besides, on the one hand, it is necessary, to respect the energy and momentum conservation laws and, on the other hand, not to spoil the parton distribution functions. To meet these requirements in the Monte Carlo simulation, we propose to use the following algorithm of partons’ permutations.
First, one prepares extensive sets of protons for each possible number of partons. Each parton momentum is sampled as proton momentum fraction, x i (where i denotes the ith parton), according to the parton distribution functions (PDFs) from CT10nnlo, next-to-next-to-leading order approximation [10], by the CTEQ (Coordinate Theoretical-Experimental Project on QCD) group based on CT10 PDFs [11], set 1, of LHAPDF (Les Houches Accord Parton Density Functions) [12] at the momentum transferred, Q 2 = 1 (GeV/c) 2 . A parton is assigned with the current quark mass of a certain flavour that is found from the probability distribution for a given x i . At this point, there is no restriction to the parton content in terms of x i .
Second, the generated partons are rearranged between protons in order to have i x i 1 and e i 1 for all protons, where e i = m i 2 m p 2 cosh 2 y beam + x i 2 tanh 2 y beam , m i is a mass of the ith parton, m p is the mass of the proton, and y beam is the proton beam rapidity, y beam = arccosh ( s / ( 2 m p ) ) .
Thus, it is ensured that the total energy and momenta of the partons are equal to the proton’s energy and momentum. To perform this, two protons (necessarily with the same number of partons) are selected and two randomly chosen partons, one taken from each proton, are exchanged between them. This permutation is accepted only if such a rearrangement decreases i x i and i e i (if the sums are larger than 1 in the previous step) or increases i x i and i e i (if the sums are smaller than 1 in the previous step) simultaneously for both considered protons.
The computation time of such an algorithm grows crucially with the number of considered partons. Therefore, one actually cannot reach exactly the value of 1 for i x i and i e i . Thus, the permutations stop at some reasonable number of iterations (i.e., when the improvement caused by the exchange of the pair of partons becomes negligible), which leads to the lack of i x i and i e i compared to 1. Therefore, to compensate, an object called a gluon cloud with the momentum fraction, x g cloud = 1 i x i and e g cloud = 1 i e i is introduced. It is assigned the mass, m g cloud = E g cloud 2 p g cloud 2 , with E g cloud the energy and p g cloud the momentum of the gluon cloud, and is considered a parton that later is also used to form strings.
Finally, a string is stretched between two randomly selected partons belonging to two randomly selected protons (but with the same number of partons). A string is accepted only if its energy is sufficient to decay at least into two pions at rest meaning that s x 1 x 2 2 m π [13] with m π denoting the pion mass. Thus, all partons from the two colliding protons should form such strings. Otherwise, another random proton pair to be sought.

2.3. Strings Motion in Transverse Plane

Initial transverse positions of the strings are sampled from the Gaussian distribution respecting the cross-section area of a proton. Then, the strings are put in motion in the transverse plane according to the system of differential equations of the second order,
r ¨ i = f i j = r i j r ˜ i j ( g N σ ) m σ 2 K 1 ( m σ r ˜ i j ) ,
defined by the two-dimensional Yukawa interaction [14]. Here, r i j is a 2D distance between the ith and jth strings, r ˜ i j = r i j 2 + s string 2 is a regularized 2D distance in which there is no Coulomb singularity at small r i j due to the introduction of s string = 0.176 fm [14] that is a genuine string width, unlike the effective string width, which is a result of quantum fluctuations; g N σ = 0.2 is the QCD string self-interaction coupling [15]; m σ = 0.6 GeV/c 2 [14] is the mass of the σ -meson that is a mediator of the force between strings; K 1 is a modified Bessel function due to the cylindrical symmetry of the problem. Following Equation (3), strings are considered moving as a whole, with no kinks.
It is worth mentioning that, in the case of many overlapping strings, their interaction will lead to the formation of string clusters, which at some threshold, eventually may cause the quark–gluon plasma formation. That is to say that some density of overlapping strings will be sufficient to eliminate the quark-antiquark chiral condensate, q q ¯ , because of its suppression in the vicinity of colour strings [14], which will lead to the restoration of the chiral symmetry.
Transverse strings’ evolution governed by Equation (3) can be terminated at some proper time, τ , which will affect the final string density. We consider two scenarios. In the first scenario, the system of strings is frozen at some conventional time (same for all events) before the start of string fragmentation at τ = 1.5 fm/c [14]. In the second case, the system is fixed at τ deepest time, obtained event-by-event. This is the time that takes for some initial configuration of strings in the transverse plane to reach the global minimum of the potential energy, changing according to Equation (3). Thus, τ deepest depends not only on the initial string transverse positions, but also on the number of strings, hence on the collision energy. It is also worth emphasising that, for τ = 1.5 fm/c, one obtains a more dilute system, while for τ deepest , the string density is highest.

2.4. Strings’ Longitudinal Evolution

In the longitudinal direction, initial strings’ end rapidities, defined by partons forming a string and expressed by the first term in
y q = ± arcsinh x q p beam m q arccosh τ 2 σ 2 2 m q 2 + 1 ,
are decreased [16] by the loss of string length in rapidity, y loss , expressed by the second term in Equation (4). It is the string tension σ that slows down the massive quarks flying outwards according to d p q / d t = σ , where p beam = s N N / 4 m p 2 is a proton beam momentum, s N N is the square of the collision centre-of-mass energy per nucleon pair, m p = 0.938 GeV/c 2 , m q is a current quark mass (we consider m u = 0.0022 GeV/c 2 , m d = 0.00479 GeV/c 2 , m s = 0.095 GeV/c 2 , and m c = 1.275 GeV/c 2 for u, d, s, and c quarks, respectively) or m g cloud .
Note that τ here is the same time as for the transverse dynamics; therefore, string evolution is synchronised in both dimensions.
An important remark is that the value of the evolution time changes not only the length of the strings, but also strings’ positions with respect to midrapidity. This is because y loss (the second term in Equation (4)) applies to both ends of the string and depends on the quark’s mass and momentum. Therefore, some strings can even lie entirely in one hemisphere if τ is large enough such that the second term in Equation (4) is greater than the first term in absolute value.

2.5. String Fusion

Considering strings density evolution in the rapidity and transverse plane dimensions, we study the fusion of strings in the string’s final configuration. Having a finite size in the transverse plane caused by the colour confinement, strings may overlap with the different degrees of intersection.
As a reference scenario to string fusion, we consider the case of independent sources, so, at this step, the possibility of strings interacting is ignored, even when overlapped.
Our main interest here is to consider string–string interaction effects in a form of fusion leading to a new type of string with higher tension. There are several options on how to account for this; see [17] for corresponding studies.
First, one may assume that the colour field changes only in the areas where strings overlap and stays the same in the rest of the transverse strings’ area. Another possibility is that, if strings overlap, even on an arbitrarily small area, the strings form a cluster, which has the uniformly modified colour field on the area of the union of all these strings. However, so far, the attempts were in vain to find the observables sensitive to the type of overlap mode. Therefore, we follow a cellular approach in the transverse plane developed in Ref. [17]. This means that we only look after strings’ centre position on the transverse grid with the constant cell size, which reflects the string transverse area. Thereby, the only strings to interact are the strings whose centres lie in the same transverse cell. This simplification is beneficial from the point of view of computing resources, and besides, it gives the same results as a full-scale simulation that should account for the particular areas of overlap.
One has to mention that the fusion problem becomes three-dimensional since in the rapidity space there are the fluctuations of both the lengths of strings and the strings’ positions with respect to midrapidity. Thus, for strings occupying the same cell in the transverse plane, one has to search for their overlaps in the rapidity space to perform fusion in slices in rapidity.
Technically, k overlapping strings that have centres in some transverse cell and overlap over some rapidity interval are replaced by a cluster of k strings, which can be seen as a new string with the modified particle production characteristics [18]:
μ k = μ 0 k ,
while μ 0 is the mean particle multiplicity per rapidity unit for an independent string. The parts of the initial strings that do not overlap in the rapidity space (even though string centres may lie in the same transverse cell) maintain the same colour field characteristics, and are considered also as separate strings. Thus, the process of string fusion makes overlapping strings shorter in the rapidity space, although the strings obtain the modification factor k as stated in Equation (5), which results in an overall decrease of mean multiplicity for interacting sources.
One can qualitatively explain Equation (5) as follows. Given the case of the overlap of a few strings, their colour field vectors could be randomly oriented with respect to each other. To find the total colour field inside the cluster of overlapping strings, one can consider this problem as a random walk in the vector colour space, thus giving not the factor of k, but k for k overlapping strings compared to a single independent string. The first consideration of this behaviour was proposed in Ref. [19], while a rigorous comparison of the string percolation approach to the colour–glass condensate model can be found in Ref. [20].
The modification of the transverse momentum, p T , spectrum was supposed [21] to be performed according to the mean, p T k = p 0 k 4 , where p 0 stands for the mean transverse momentum of particles produced by independent sources. Thus, strings’ interaction causes the increase in particle mean- p T . However, in Refs. [9,22], instead p 0 is made dependent on the total number, M, of strings in an event as power, M β . The β = 1.16 [ 1 ( ln s 2.52 ) 0.19 ] depends on the collision energy, and for LHC energies, β 0.25 . This gives us the idea to combine these two approaches and to introduce the factor k β instead of k 4 . However, let us stress that the k here is a number of strings in a string cluster in some rapidity interval in some transverse cell, while M from Refs. [9,22] is the overall number of particle sources in an event. Eventually, we apply the following modification for the mean- p T from a cluster:
p T k = p 0 k β ,
which would be of special interest in our future studies of pp interactions at lower collision energies, where β changes sign according to the fit in Ref. [22].
In the current implementation of the model, at this moment of simulation, the string system is already frozen, and no transverse position is assigned to the cluster of strings (new strings), as soon as there is no need of it in the further calculations. However, one can introduce the transverse position as an average value over transverse positions of k fused strings.
Thus, after the initial strings’ dynamics in transverse and longitudinal spaces is considered and after the string fusion procedure is completed, one has a mixture of particle production sources: (i) strings (or parts of strings in rapidity) that did not undergo the process of fusion (let us call them single strings) so that one has to use for them k = 1 in Equations (5) and (6) and (ii) new strings (which we call clusters of strings) with k = 2 , 3, or other values in Equations (5) and (6) that reflect the degree of the strings’ overlap for these clusters.

2.6. Effective String Hadronisation

We perform effective string hadronisation by dividing the string in the y-direction into units of length, ε , to relate particles’ rapidities with the corresponding pieces of strings. This gives the mean multiplicity, N ε = μ 0 ε k , for ε unit and the actual multiplicity, N ε , sampled from the Poisson distribution with the mean N ε . Thus, multiplicity from the jth string in an event is the sum of its ε pieces: N j = l N ε l , and event multiplicity is the sum over multiplicities from all strings in event: N = j N j . Particles’ rapidities are sampled from the Gauss distribution for each ε unit with the mean equal to the centre of the ε unit and with variance equal to ε . It is important to mention that this ε -division is only a technical solution to approximate the string fragmentation. Therefore, this procedure is the same for the single strings and clusters of strings. The substantial changes when considering string fusion appear in the factor k for mean multiplicities (5) and mean transverse momenta (6), as well as for the Schwinger-like probabilities for the production of particle species.
Particles’ transverse momentum is sampled from the distribution,
f ( p T ) = π p T 2 p T k 2 exp π p T 2 4 p T k 2 ,
which corresponds to the Schwinger mechanism of particle production [23,24,25], with p T k from Equation (6). Here, particle species, i, are assigned according to Schwinger-like probabilities, exp ( π m i 2 / σ k 2 β ) , which is consistent with Equation (7), since we set σ = 4 p 0 2 . In the model, we consider π , K, and p particles and ρ -resonance, with the latter decaying into two charged pions. With particle masses, the longitudinal component of momentum, p z , is found and, thus, the pseudorapidity,
η = 1 2 ln | p | + p z | p | p z ,
is calculated for each particle, where | p | = p T 2 + p z 2 is an absolute value of particle total momentum.

2.7. Model Tuning

To tune our model, we simultaneously fit the ALICE experiment pp inelastic interaction data at s = 900 GeV on the charged particle multiplicity distribution, η -spectrum [26], and p T N correlation function [27], where the denote the averaging over events. Global model parameters not affecting the result are: the hadronisation parameter, ε = 0.1 and the transverse grid cell size of 0.3 fm. The best fit results are obtained for the model free parameters μ 0 = 0.87 and p 0 = 0.38 GeV/c. The values of the parameters for multi-pomeron distribution and for the transverse dynamics of strings are fixed to those from Section 2.1, Section 2.2, Section 2.3, Section 2.4 and Section 2.5.

3. Results

We consider a few cases of the evolution of strings’ system:
  • no transverse string dynamics, but for longitudinal dynamics τ = 1.5 fm/c is used;
  • both the transverse and longitudinal string dynamics occur till τ = 1.5 fm/c;
  • both the transverse and longitudinal dynamics occur till τ = τ deepest (which changes event to event depending on the initial string configuration in the transverse plane).
For each of the above cases, the results are compared for independent and interacting sources, as the evolution time determined the three-dimensional event string density, which stimulated string fusion.
The model results for different sets of string dynamics parameters are compared with PYTHIA8.3 Monte Carlo generator simulations [28] with a default list of settings (the Monash pp tune for non-diffractive events), which are also plotted with the colour reconnection option switched off and on. A qualitative comparison with the ALICE data on pp collisions at s N N = 900 GeV is also performed.

3.1. Global Observables

Starting from quite a complicated string structure in the rapidity space (strings’ lengths and positions of strings’ ends significantly vary event-by-event), a symmetrical inclusive η -distribution is obtained in the whole available η -range.
We tune our model with interacting particle production sources, evolving till τ deepest , to describe the ALICE data [26] at midrapidity (Figure 1). To describe the charged particle multiplicity distribution, the events with N ch = 0 are removed since these events are highly influenced by diffractive processes that are not considered in the model. For comparison, we also present here the results for the case of the absence of string interactions with the parameters obtained from the fit. One can see from the multiplicity distribution (Figure 1, left) and pseudorapidity spectrum (Figure 1, right) that the inclusion of the string fusion effect resulted into a relatively smaller average multiplicity.
One should keep in mind that, in the case of the multiplicity variable, the model can be retuned so that the no-interaction scenario would describe these spectra. One can see how important is the string fusion solely at the level of correlations, as described in the following Subsections.
Finally, it should be noted that p 0 = 0.38 GeV/c is selected in order to describe p T = 0.489 ± 0.001 GeV/c measured by ALICE [27].

3.2. p T –N Correlation Function

We study p T N correlations in the acceptance of 0.15 < p T < 4 GeV/c, | η | < 0.8 . The model results with interacting particle production sources, evolving till τ deepest , tuned to the ALICE data [27], are compared to the PYTHIA simulations and to the findings with other model options in Figure 2.
From Figure 2 one can see that the model results for independent strings (empty markers, except the squares) almost coincide, and they exhibit no dependence of p T on N regardless the transverse string dynamics. Thus, one can conclude that there is no collective behaviour without strings’ interaction.
As for the interacting particle production sources, one can see more significant correlations in Figure 2 for the cases with the larger strings’ density formed in the event configuration. Namely, the first case is the weak p T N dependence for string fusion and fixed strings’ positions in the transverse plane (full grey crosses): fusion occurres quite rarely. Next case is the correlation function for the system of strings evolving till τ = 1.5 fm/c (full blue circles). The largest p T N correlation in the model is seen for the largest density of strings at τ deepest , where the string density is the highest one.
One could also see that the colour reconnection option in PYTHIA (full green squares) plays formally a similar role as our string fusion mechanism, although PYTHIA without the colour reconnection included (empty green squares) still has some background correlations.
In brief, even though one may try to tune global observables without collective effects of some sort, it is not possible to explain the growth of the event mean- p T with the multiplicity.

3.3. Multiplicity Correlations

Now, we study the multiplicity correlations in pseudorapidity in terms of the correlation coefficient, b B F , which represents the slope of the correlation function, N B ( N F ) [29],
b B F = d N B ( N F ) d N F | N F = N F ,
defined for multiplicities N F and N B in two pseudorapidity intervals, (so-called “forward” and “backward”) separated by a Δ η gap.
Correlation function, N B ( N F ) , indicates how average value of N B depends on the selected value of N F . In case of a linear dependence the correlation coefficient reads [30]:
b corr [ N F , N B ] = N F N B N F N B N B 2 N B 2 .
It is important to mention that, in the case of a boost invariance in rapidity, i.e., when the string longitudinal dynamics is ignored, the fluctuations in the number of strings would provide the same input to the long-range correlations at all positions of forward and backward rapidity windows. This is not the case in our model, as with the loss of translational invariance, one might obtain different “volume” fluctuations (fluctuating number of strings acting as particle-emitting sources) for different Δ η .
We study b corr [ N F , N B ] as a function of the distance Δ η between forward and backward pseudorapidity acceptances of width δ η = 0.2 each in order to access the information about particle sources and their distribution over η . As one can see from Figure 3, the behaviour of b corr [ N F , N B ] versus Δ η is not trivial for a purely symmetric η -distribution, shown in Figure 1, right.
Figure 3 shows two groups of results. The first is a set of plots for b corr [ N F , N B ] versus Δ η , which for τ deepest , exhibit almost no dependence on Δ η —both with (full red triangles) and without (empty red triangles) string fusion. This is because τ deepest is not that large (the average τ deepest at s = 900 GeV is about 0.73) and, therefore, the colour strings stay long enough to produce particles into both forward and backward η -windows. This is why one sees a strong correlation which does not weaken with Δ η . Note that the magnitudes of b corr [ N F , N B ] are different for various model assumptions, and the calculations for τ deepest with string fusion are closer to the experimental values.
The second group shows the decreasing tendency of b corr [ N F , N B ] for all other model cases. The values of b corr fall down with the increase of Δ η since the colour strings are significantly shrunk by y loss with a τ = 1.5 fm/c evolution time. Therefore, the strings impact more independently in forward and backward η windows the more the windows are separated. Thus, the forward–backward correlation weakens with Δ η . As a cross-check, one can notice that the results for τ = 1.5 fm/c and for the case without transverse string dynamics almost coincide for independent strings (open circles and crosses). This is because in these both cases, τ = 1.5 fm/c for the longitudinal dynamics.
For the interacting strings, the results for τ = 1.5 fm/c lie below the case without transverse string dynamics because of the higher frequency of string fusion and, consequently, lower multiplicity in the former case then in the latter one. This is also correct for any system evolution time: the results for interacting strings lie below the corresponding plots for independent strings because of the lower multiplicity caused by string fusion in the former case. The same prediction on the suppression of b corr [ N F , N B ] due to string fusion was obtained earlier for strings that are infinite in rapidity [32].
A similar behaviour one can observe for the PYTHIA simulations: colour reconnection gives an effect similar to that for the string fusion mechanism—an overall decrease of multiplicity and that of b corr [ N F , N B ] .
The model results are also compared with the ALICE measurements [31]. The simultaneous fit to the experimental d N / d η and P ( N ch ) distributions allows us to reproduce the experimental value of b corr [ N F , N B ] only in the midrapidity (interacting strings evolving till τ deepest at Δ η = 0 ). However, for the moment, the model does not follow the steady decrease of the data b corr [ N F , N B ] : the full red triangles versus the black line.
To conclude here, we note that the evolution time and the relevant string density are found playing the major role in the magnitude of long-range correlation coefficient and in its behaviour with the distance, Δ η , separating forward and backward pseudorapidity acceptance intervals.

3.4. Multiplicity Fluctuations with Strongly Intensive Observable

To avoid the dependence of the results on the system “volume”, we study joint multiplicity fluctuations in forward and backward pseudorapidity intervals in terms of the special fluctuation measure, Σ [ N F , N B ] [33],
Σ [ N F , N B ] = N F ω [ N B ] + N B ω [ N F ] 2 ( N F N B N F N B ) N F + N B ,
where ω [ A ] = A 2 A 2 / A is a scaled variance of an extensive event variable A distribution.
By construction, in the case of similar types of strings, Σ [ N F , N B ] depends neither on the system “volume” (the number of particle-emitting sources), nor on its event-by-event fluctuations. As was proposed for the family of strongly intensive quantities [34], the quantity Σ [ N F , N B ] is normalised to unity in the case of the models of independent particle production.
Figure 4 demonstrates an increase from the unity of Σ [ N F , N B ] values with the distance between forward and backward η -acceptances. One can argue that the results for independent and interacting strings split, which made Σ [ N F , N B ] dependent on the types of particle production sources. The results for τ = 1.5 fm/c (blue circles) tended to decrease for the case of string clusters for the largest Δ η , which resembled also the difference for the PYTHIA plots (green squares). The results for τ deepest lied even lower (red triangles), although it is complicated to distinguish the two model regimes here.
A considerable difference with the PYTHIA predictions can be explained here by the absence of short-range effects in our model that was currently aimed at the study of long-range correlation effects.

3.5. Connection of Σ [ N F , N B ] with Cumulants and Factorial Cumulants

To illustrate the beauty of the mathematics standing behind the observable Σ [ N F , N B ] , let us briefly discuss the Σ [ N F , N B ] connection with the other often used objects, namely the cumulants and factorial cumulants. In the case of the joint probability distribution, P N F , N B , the cumulants and factorial cumulants can be introduced for any linear combination of N F and N B . Following the notations of Ref. [35], let us define a selected linear combination as a q-vector, q ( x ) = x 1 · N F + x 2 · N B . Let us consider these two specific q-vectors: q ( a ) = N F + N B and q ( b ) = N F N B , where a = 1 , 1 and b = 1 , 1 . The auxiliary q-vectors that one needs in order to calculate cumulants and factorial cumulants are: q ( a 2 ) = q ( a ) , q ( b 2 ) = q ( a ) , and q ( a b ) = q ( b ) .
Then, following Ref. [35], one can immediately express the first-order and second-order cumulants for the joint probability distribution, P N F , N B , in terms of the moments of the same distribution:
q ( a ) c = N F + N B ,
q ( b ) c = N F N B ,
q ( a ) 2 c = N F 2 N F 2 + N B 2 N B 2 + 2 · N F N B N F N B ,
q ( b ) 2 c = N F 2 N F 2 + N B 2 N B 2 2 · N F N B N F N B ,
q ( a ) · q ( b ) c = N F 2 N F 2 N B 2 + N B 2 .
In turn, the factorial cumulants can be straightforwardly expressed in terms of the cumulants:
q ( a ) fc = q ( a ) c ,
q ( b ) fc = q ( b ) c ,
q ( a ) 2 fc = q ( a ) 2 c q ( a 2 ) c = q ( a ) 2 c q ( a ) c ,
q ( b ) 2 fc = q ( b ) 2 c q ( b 2 ) c = q ( b ) 2 c q ( a ) c ,
q ( a ) · q ( b ) fc = q ( a ) · q ( b ) c q ( a b ) c = q ( a ) · q ( b ) c q ( b ) c .
One could immediately notice the connection of the presented cumulants and factorial cumulants with the strongly intensive quantity in the case of symmetrical pseudorapidity intervals, i.e., when N F = N B :
Σ [ N F , N B ] = q ( b ) 2 c q ( a ) c = 1 + q ( b ) 2 fc q ( a ) fc .
This connection is not surprising as soon as the goal of the construction of these special observables (or ratios of observables) is to suppress trivial statistical fluctuations.

3.6. Connection of Σ [ N F , N B ] with Asymmetry Coefficient

The main difference between the model results (Figure 4) for τ = 1.5 fm/c and τ deepest comes from the geometry of the strings in an event (considering transverse dynamics, longitudinal loss, and cluster formation via string fusion). Indeed, as was already mentioned in Section 3.3, the average τ deepest appears to be less than τ = 1.5 fm/c for this collision energy what results into different sets of strings and/or string clusters event-by-event producing particles in the forward and backward η -acceptances. Thus, another way to look at event-by-event multiplicity fluctuations is to construct an event observable,
C = N F N B N F + N B ,
that measures an event asymmetry in pseudorapidity [36] and to study the variance, σ C 2 , of its distribution over a set of events.
Let us consider the analytical form of σ C 2 :
σ C 2 Var ( C ) = ( N F N B ) 2 N F + N B N F N B N F + N B 2 .
This observable has been extensively studied in different models [37,38,39], e.g., in the cluster model of particle production, this observable helped to extract the average cluster size or, in other words, effective cluster multiplicity. Moreover, the extracted value both from pp and AuAu reactions was larger than the one predicted for the hadron-resonance gas model, indicating that correlations cannot be explained purely by the hadronic degrees of freedom considered in some statistical ensembles. In most of these studies, the two-stage scenario was adopted in one way or the other; some of the models take into account string interaction, while others provide a non-trivial non-uniform in rapidity fragmentation function of particle-emitting sources, but none of them consider simultaneously longitudinal dynamics and collectivity as we did in our model.
Let us show that σ C 2 is closely related to the strongly intensive quantity, Σ [ N F , N B ] . Using the assumptions,
( N F N B ) 2 N F + N B = ( N F N B ) 2 N F + N B ,
N F N B N F + N B 2 = ( N F N B ) 2 N F + N B ,
for the first and the second term of σ C 2 (24), respectively, one arrives at
σ C 2 = Var ( N F N B ) N F + N B = N F 2 2 N F N B + N B 2 N F 2 + 2 N F N B N B 2 N F + N B = Var ( N F ) + Var ( N B ) 2 cov ( N F , N B ) N F + N B = N F ω [ N B ] + N B ω [ N F ] 2 cov ( N F , N B ) N F + N B ,
where the last equality holds under the assumption of N F = N B ; “cov” stays for covariance. Thus,
σ C 2 Σ [ N F , N B ] .
The assumptions (25) and (26) are in accordance with the approximate calculations of Ref. [39], where deviations of N F and N B from N F and N B are considered to be small enough.
An important remark concerning σ C 2 needs to be made. Despite the demonstrated connection with Σ [ N F , N B ] , strictly speaking, σ C 2 is not a strongly intensive observable and may receive additional contributions coming from the “volume” fluctuations, as the latter are, in general, not cancelled out without the assumptions (25) and (26). Nevertheless, σ C 2 is considered as an important quantity since it shows the magnitude of the fluctuations in terms of event asymmetry in the longitudinal dimension, i.e., σ C 2 indirectly provides us with some information about the initial state.
Figure 5, left, shows that the variance, σ C 2 , taken from the distribution of the asymmetry coefficient C (23), indeed exhibits a behaviour similar to Σ [ N F , N B ] (Figure 4). Although both quantities gradually rise with Δ η , σ C 2 is not reaching the unity for the close η -acceptances ( Δ η = 0 ). As seen, in this region, there are the largest fluctuations of N F , B , and respectively, the assumption of small deviations of event-by-event multiplicities from their average values does not hold.
Figure 5, right, displays the values of σ C 2 above unity in a more detailed way. For two quantities that follow the Poissonian distribution, their difference would follow the Skellam distribution. Therefore, typically, the Skellam function is considered as a baseline in the studies of net-charge fluctuations. In our case, the numerator of σ C 2 can be treated as a net-hadronic charge in rapidity. The values of the model results of Var ( N F N B ) , in Figure 5, right, lie slightly above the corresponding Skellam baselines calculated as N F + N B . Therefore, the ratios of those model values over the lines are equal to σ C 2 , which is slightly above the unity.
Following Ref. [36], we may suggest that, at Δ η = 0, strings and/or string clusters predominated, which produce particles in both forward and backward η -acceptances, which suppressed the value of Σ [ N F , N B ] and σ C 2 . With the growth of Δ η , one expects more and more rare contributions of such sources. Instead, there were “short” quark–gluon strings that appeared only in the forward or in the backward η -intervals.
One has to mention that the same qualitative behaviour of Σ [ N F , N B ] with Δ η was obtained in Ref. [40], where the authors considered only strings infinite in rapidity, but included the short-range correlations. Thus, one may argue that both the absence of short-range correlations and the presence of short strings in our model gave the same results for Σ [ N F , N B ] that indicates the importance of the coherent application of the two approaches.

3.7. Forward Multiplicity Studies

In the analysis of the experimental data, along with the choice of stable variables, such as Σ [ N F , N B ] , another possibility to control an event’s initial conditions (such as unavoidable fluctuations due to, for instance, Fermi motion) is to select for analyses some sort of similar events.
In the context of the present study, it is necessary to make an important remark. We followed the picture of multi-pomeron exchange, which is needed to restore the unitarity of the scattering amplitudes for a pomeron intercept larger than one. Each cut pomeron resulted in the formation of two colour strings, which in total led to the formation of the number of strings in an event. In turn, each string had to be attached to either valence or sea quarks (at the moment, diquarks are not considered), which defines the number of partons “participating” in the interaction. Therefore, events differed by the number of participating partons, e.g., the number of colour strings that form the effective “volume” of the system radiating particles.
In this way, here, we repeat the procedure used in the ALICE experiment to define multiplicity-based classes of events in proton–proton collisions following Ref. [41].
Knowing η for each particle, we plot the multiplicity of particles that fly in the V0-detector of ALICE [42], whose acceptance covers 2.8 < η < 5.1 and 3.7 < η < 1.7 . This was to avoid auto-correlations in the TPC (Time-Projection Chamber) acceptance, where one measures particles for the analysis. We divide the distribution into five classes, each of which contained 20% of the events. For the analysis, the 0–20%, 20–40%, and 40–60% classes are taken and the corresponding strongly intensive quantity (11) is calculated for each set of a specific class of events.
It is important to mention that the acceptance of the ALICE V0-detector covers the region where quasi-diffractive processes play a role. Furthermore in this region, it is crucial to properly model the string fragmentation close to the string endpoints, e.g., to take into account the possibility to have a diquark as an endpoint, which can produce additional baryon yields, etc. Finally, these detectors have non-ideal resolution, and hence, this factor has to be included in the simulations as well. All these effects are not taken into account in this version of the model. Because of this, we do not perform any direct comparison with the available experimental data [41] obtained in very narrow multiplicity-based event classes. We only study whether Σ [ N F , N B ] depends on an event “type” and whether the property of strong intensity is violated for Σ [ N F , N B ] .
In Figure 6, one can see that, for all multiplicity-based event classes, Σ [ N F , N B ] (11) exhibits the same trend of an increase with Δ η . Moreover, one can see that, for more central events, this quantity reaches larger values. The same effect was observed by the ALICE experiment [41]. This is another manifestation of the string fusion effect because, for more “central” events, i.e., events with the higher forward multiplicity, one has a larger number of cut pomerons and, consequently, a larger number of strings. This means that one has a higher probability of strings interaction and possibly could produce string clusters of higher string tension. A similar analysis with the same conclusion was performed recently in Ref. [40] for interacting strings that are infinite in rapidity. Again, let us to emphasise the importance of taking into account both the short-range correlations that were modelled in Ref. [40] and the initial conditions that are considered in this model, e.g., strings of short length and/or shift with respect to midrapidity, as these two effects give similar results and cannot be easily disentangled by analysing Σ [ N F , N B ] only.

4. Discussion

Within the developed Monte Carlo multi-pomeron exchange model, we have obtained results on multiplicity and transverse momentum correlations for pp inelastic interactions at s = 900 GeV. The key elements of the model based on the event-by-event formation and fusion of longitudinally extended colour strings are the long-range correlations naturally introduced. The correlations appear from the fluctuations in the number of “long” strings that produced particles simultaneously in the forward and backward pseudorapidity intervals.
Tuned to some global observables in the ALICE data, such as the rapidity distribution of charged particles for inelastic pp interactions at s = 900 GeV and a p T N correlation function, the model with string fusion and the evolution time τ deepest gives results in qualitative agreement with the experiment in b corr [ N F , N B ] versus the distance Δ η between the forward and backward pseudorapidity acceptance intervals.
Strings’ evolution time and the relevant string densities are found to play a major role in the magnitude of long-range correlation coefficient and in its slope with the distance Δ η .
Results are found in trend with PYTHIA 8.3 event generator simulations, where the Colour Reconnection option plays a formally similar role, although the underlying physics of PYTHIA 8.3 is different.
Although it might seem, from Figure 3, that our model contains also short-range correlations that are visible at small Δ η in the b corr [ N F , N B ] plots, we would like to stress that this is just a result of the modification of the long-range correlation background because of the appearance of short strings (with respect to forward and backward η -windows). Thus, in the present study, the short-range correlations have been effectively implemented in the model via the presence of short strings that independently impact the forward and backward pseudorapidity intervals.
A similar phenomenon has exhibited itself in the case of the application of the strongly intensive quantity Σ aimed at eliminating the trivial effects of the so-called volume fluctuations, as it drives Σ increase with Δ η , which, in the case of the absence of longitudinal dynamics, may have originated only from the short-range correlations, as was shown in Ref. [43]. In future, we plan to introduce genuine short-range correlations related to string hadronisation, resonance decays, or mini-jets and profit from their interplay with the string fusion mechanism.
The dependence of Σ on the multiplicity-based event classes is found in accordance with the preliminary ALICE data: the model predicted larger values of Σ for pp interactions with higher forward multiplicity.
We have also shown that Σ is connected with such often used observables as cumulants, factorial cumulants, and the variance, σ C 2 , of the asymmetry coefficient distribution, which is also proven in the model simulation.
In this study, we have also shown that the introduction of the string fusion effect in the model results in the enhanced p T N correlation as the fused strings, on average, decay to the smaller number of particles with higher p T s.

Author Contributions

Conceptualisation, D.P., E.A. and G.F.; software, D.P. and E.A.; validation, D.P. and E.A.; formal analysis, D.P.; writing—original draft preparation, D.P.; writing—review and editing, E.A.; visualisation, D.P.; supervision, G.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by St. Petersburg State University, Russia (project No. 94031112).

Data Availability Statement

Not applicable.

Acknowledgments

The authors are grateful to Vladimir Kovalenko for fruitful discussions.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of the data; in the writing of the manuscript; nor in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
ALICE  A Large Ion Collider Experiment
CTEQ  Coordinate Theoretical-Experimental Project on QCD
CT10nnlo  CTEQ version 10 next-to-next-to-leading order
LHAPDF  Les Houches Accord Parton Density Functions
LHC  Large Hadron Collider
pp  proton-proton
QCD  Quantum ChromoDynamics
TPC  Time-Projection Chamber
1D, 2D, 3D  1-, 2-, 3-dimensional

References

  1. Khachatryan, V. et al. [CMS Collaboration]. Observation of long-range, near-side angular correlations in proton–proton collisions at the LHC. J. High Energy Phys. 2017, 2010, 91. [Google Scholar] [CrossRef] [Green Version]
  2. Adam, J. et al. [ALICE Collaboration]. Enhanced production of multi-strange hadrons in high-multiplicity proton–proton collisions. Nat. Phys. 2017, 13, 535–539. [Google Scholar] [CrossRef] [Green Version]
  3. Artru, X. Classical string phenomenology. How strings work. Phys. Rep. 1983, 97, 147–171. [Google Scholar] [CrossRef]
  4. Braun, M.; Pajares, C. A Probabilistic model of interacting strings. Nucl. Phys. B 1993, 390, 542–558. [Google Scholar] [CrossRef]
  5. Dumitru, A.; Gelis, F.; McLerran, L.; Venugopalan, R. Glasma flux tubes and the near side ridge phenomenon at RHIC. Nucl. Phys. A 2008, 810, 91–108. [Google Scholar] [CrossRef] [Green Version]
  6. Kaidalov, A.B.; Ter-Martirosyan, K.A. Pomeron as quark–gluon strings and multiple hadron production at SPS-collider energies. Phys. Lett. B 1982, 117, 247–251. [Google Scholar] [CrossRef]
  7. Capella, A.; Sukhatme, U.; Tan, C.-I.; Tran Thanh Van, J. Dual parton model. Phys. Rep. 1994, 236, 225–329. [Google Scholar] [CrossRef]
  8. Vechernin, V.V.; Belokurova, S.N. The strongly intensive observable in pp collisions at LHC energies in the string fusion model. J. Phys. Conf. Ser. 2020, 1690, 012088. [Google Scholar] [CrossRef]
  9. Armesto, N.; Derkach, D.A.; Feofilov, G.A. pt-multiplicity correlations in a multi-Pomeron-exchange model with string collective effects. Phys. Atom. Nucl. 2008, 71, 2087–2095. [Google Scholar] [CrossRef]
  10. Gao, J.; Guzzi, M.; Huston, J.; Lai, H.-L.; Li, Z.; Nadolsky, P.; Pumplin, J.; Stump, D.; Yuan, C.-P. CT10 next-to-next-to-leading order global analysis of QCD. Phys. Rev. D 2014, 89, 033009. [Google Scholar] [CrossRef] [Green Version]
  11. Lai, H.-L.; Guzzi, M.; Huston, J.; Li, Z.; Nadolsky, P.M.; Pumplin, J.; Yuan, C.-P. New parton distributions for collider physics. Phys. Rev. D 2010, 82, 074024. [Google Scholar] [CrossRef] [Green Version]
  12. Buckley, A.; Ferrando, J.; Lloyd, S.; Nordström, K.; Page, B.; Rüfenacht, M.; Schönherr, M.; Watt, G. LHAPDF6: Parton density access in the LHC precision era. Eur. Phys. J. C 2015, 75, 132. [Google Scholar] [CrossRef] [Green Version]
  13. Kovalenko, V.N. Modeling of exclusive parton distributions and long-range rapidity correlations in proton–proton collisions at the LHC energies. Phys. Atom. Nucl. 2013, 76, 1189–1195. [Google Scholar] [CrossRef] [Green Version]
  14. Kalaydzhyan, T.; Shuryak, E. Collective interaction of QCD strings and early stages of high-multiplicity pA collisions. Phys. Rev. C 2014, 90, 014901. [Google Scholar] [CrossRef] [Green Version]
  15. Kalaydzhyan, T.; Shuryak, E. Self-interacting QCD strings and string balls. Phys. Rev. D 2014, 90, 025031. [Google Scholar] [CrossRef] [Green Version]
  16. Shen, C.; Schenke, B. Dynamical initial-state model for relativistic heavy-ion collisions. Phys. Rev. C 2018, 97, 024907. [Google Scholar] [CrossRef] [Green Version]
  17. Braun, M.A.; Kolevatov, R.S.; Pajares, C.; Vechernin, V.V. Correlations between multiplicities and average transverse momentum in the percolating color strings approach. Eur. Phys. J. C 2004, 32, 535–546. [Google Scholar] [CrossRef]
  18. Braun, M.A.; Pajares, C.; Ranft, J. Fusion of strings vs. percolation and the transition to the quark–gluon plasma. Int. J. Mod. Phys. A 1999, 14, 2689–2704. [Google Scholar] [CrossRef] [Green Version]
  19. Biro, T.S.; Nielsen, H.B.; Knoll, J. Color rope model for extreme relativistic heavy ion collisions. Nucl. Phys. B 1984, 245, 449–468. [Google Scholar] [CrossRef]
  20. Braun, M.A.; Dias de Deus, J.; Hirsch, A.S.; Pajares, C.; Scharenberg, R.P.; Srivastava, B.K. De-confinement and clustering of color sources in nuclear collisions. Phys. Rep. 2015, 599, 1–50. [Google Scholar] [CrossRef] [Green Version]
  21. Braun, M.A.; Del Moral, F.; Pajares, C. Percolation of strings and the first RHIC data on multiplicity and tranverse momentum distributions. Phys. Rev. C 2002, 65, 024907. [Google Scholar] [CrossRef] [Green Version]
  22. Kovalenko, V.; Feofilov, G.; Puchkov, A.; Valiev, F. Multipomeron model with collective effects for high-energy hadron collisions. Universe 2022, 8, 246. [Google Scholar] [CrossRef]
  23. Schwinger, J. On gauge invariance and vacuum polarization. Phys. Rev. 1951, 82, 664–679. [Google Scholar] [CrossRef]
  24. Gurvich, E.G. The quark–antiquark pair production mechanism in a quark jet. Phys. Lett. B 1979, 87, 386–388. [Google Scholar] [CrossRef]
  25. Casher, A.; Neuberger, H.; Nussinov, S. Chromoelectric-flux-tube model of particle production. Phys. Rev. D 1979, 20, 179–188. [Google Scholar] [CrossRef]
  26. Aamodt, K. et al. [ALICE Collaboration]. Charged-particle multiplicity measurement in proton–proton collisions at s = 7 TeV with ALICE at LHC. Eur. Phys. J. C 2010, 68, 345–354. [Google Scholar] [CrossRef] [Green Version]
  27. Aamodt, K. et al. [ALICE Collaboration]. Transverse momentum spectra of charged particles in proton–proton collisions at s = 900 GeV with ALICE at the LHC. Phys. Lett. B 2010, 693, 53–68. [Google Scholar] [CrossRef] [Green Version]
  28. Bierlich, C.; Chakraborty, S.; Desai, N.; Gellersen, L.; Helenius, I.; Ilten, P.; Lönnblad, L.; Mrenna, S.; Prestel, S.; Preuss, C.T.; et al. A comprehensive guide to the physics and usage of PYTHIA 8.3. arXiv 2022, arXiv:2203.11601. [Google Scholar] [CrossRef]
  29. Uhlig, S.; Derado, I.; Meinke, R.; Preissner, H. Observation of charged particle correlations between the forward and backward hemispheres in pp collisions at ISR energies. Nucl. Phys. B 1978, 132, 15–28. [Google Scholar] [CrossRef]
  30. Capella, A.; Tran Thanh Van, J. Long-range rapidity correlations in hadron-nucleus interactions. Phys. Rev. D 1984, 29, 2512–2516. [Google Scholar] [CrossRef]
  31. Adam, J. et al. [The ALICE Collaboration]. Forward-backward multiplicity correlations in pp collisions at s = 0.9, 2.76 and 7 TeV. J. High Energy Phys. 2015, 5, 97. [Google Scholar] [CrossRef] [Green Version]
  32. Belokurova, S.; Vechernin, V. Long-range correlations between observables in a model with translational invariance in rapidity. Symmetry 2020, 12, 1107. [Google Scholar] [CrossRef]
  33. Andronov, E.V. Influence of the quark–gluon string fusion mechanism on long-range rapidity correlations and fluctuations. Theor. Math. Phys. 2015, 185, 1383–1390. [Google Scholar] [CrossRef]
  34. Gorenstein, M.I.; Gaździcki, M. Strongly intensive quantities. Phys. Rev. C 2011, 84, 014904. [Google Scholar] [CrossRef]
  35. Kitazawa, M.; Luo, X. Properties and uses of factorial cumulants in relativistic heavy-ion collisions. Phys. Rev. C 2017, 96, 024910. [Google Scholar] [CrossRef] [Green Version]
  36. Back, B.B. et al. [PHOBOS Collaboration]. Forward-backward multiplicity correlations in s NN = 200 GeV Au+Au collisions. Phys. Rev. C 2006, 74, 011901(R). [Google Scholar] [CrossRef] [Green Version]
  37. Haussler, S.; Abdel-Aziz, M.; Bleicher, M. Forward-backward charge fluctuations at RHIC energies. Nucl. Phys. A 2007, 785, 253–256. [Google Scholar] [CrossRef] [Green Version]
  38. Brogueira, P.; Dias de Deus, J. Forward-backward rapidity correlations in a two-step scenario. Phys. Lett. B 2007, 653, 202–205. [Google Scholar] [CrossRef] [Green Version]
  39. Bzdak, A.; Woźniak, K. Forward-backward multiplicity fluctuations in heavy nuclei collisions in the wounded nucleon model. Phys. Rev. C 2010, 81, 034908. [Google Scholar] [CrossRef] [Green Version]
  40. Belokurova, S.; Vechernin, V. Using a strongly intense observable to study the formation of quark-gluon string clusters in pp collisions at LHC energies. Symmetry 2022, 14, 1673. [Google Scholar] [CrossRef]
  41. Erokhin, A. Forward-backward multiplicity correlations with strongly intensive observables in pp collisions. In Proceedings of the VI International Conference on the Initial Stages of High-Energy Nuclear Collisions (IS2021), Rehovot, Israel, 10–15 January 2021; Available online: https://indico.cern.ch/event/854124/contributions/4134683/ (accessed on 15 May 2023).
  42. Aamodt, K. et al. [The ALICE Collaboration]. The ALICE experiment at the CERN LHC. J. Instrum. 2008, 3, S08002. [Google Scholar] [CrossRef] [Green Version]
  43. Andronov, E.; Vechernin, V. Strongly intensive observable between multiplicities in two acceptance windows in a string model. Eur. Phys. J. A 2019, 55, 14. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Comparison of model results for independent particle production sources (blue lines) and for interacting strings (red lines) with the ALICE data [26] (black squares) for pp inelastic interactions at s = 900 GeV: (left) charged particle multiplicity distribution, P ( N ch ) = P full ( N ch ) / ( 1 P full ( 0 ) ) in the pseudorapidity | η | < 1 acceptance for the full p T range of particles, and (right) η -distribution for the full p T range. Here, P full denotes the charged particle multiplicity distribution that includes events without registered particles and N ch denotes the charged particle multiplicity.
Figure 1. Comparison of model results for independent particle production sources (blue lines) and for interacting strings (red lines) with the ALICE data [26] (black squares) for pp inelastic interactions at s = 900 GeV: (left) charged particle multiplicity distribution, P ( N ch ) = P full ( N ch ) / ( 1 P full ( 0 ) ) in the pseudorapidity | η | < 1 acceptance for the full p T range of particles, and (right) η -distribution for the full p T range. Here, P full denotes the charged particle multiplicity distribution that includes events without registered particles and N ch denotes the charged particle multiplicity.
Physics 05 00042 g001
Figure 2. p T N correlation function calculated in | η | < 0.8 pseudorapidity acceptance with 0.15 < p T < 4 GeV/c for inelastic pp interactions at s = 900 GeV. Model results for string system evolving till τ deepest , represented by full red triangles, are tuned to follow ALICE data [27] (the black line: data points are connected to guide the eye). See text for details.
Figure 2. p T N correlation function calculated in | η | < 0.8 pseudorapidity acceptance with 0.15 < p T < 4 GeV/c for inelastic pp interactions at s = 900 GeV. Model results for string system evolving till τ deepest , represented by full red triangles, are tuned to follow ALICE data [27] (the black line: data points are connected to guide the eye). See text for details.
Physics 05 00042 g002
Figure 3. b corr [ N F , N B ] (10) as a function of the distance, Δ η , between forward and backward pseudorapidity acceptance intervals, where the forward, N F , and backward, N B , multiplicities, are calculated for inelastic pp interactions at s = 900 GeV (see text for details). Particles are selected with 0.3 GeV/c < p T < 1.5 GeV/c. The black line is drawn to guide the eye trough the ALICE data [31].
Figure 3. b corr [ N F , N B ] (10) as a function of the distance, Δ η , between forward and backward pseudorapidity acceptance intervals, where the forward, N F , and backward, N B , multiplicities, are calculated for inelastic pp interactions at s = 900 GeV (see text for details). Particles are selected with 0.3 GeV/c < p T < 1.5 GeV/c. The black line is drawn to guide the eye trough the ALICE data [31].
Physics 05 00042 g003
Figure 4. Σ [ N F , N B ] (11) as a function of the distance, Δ η , between forward and backward pseudorapidity acceptance intervals, where N F and N B multiplicities are calculated for inelastic pp interactions at s = 900 GeV (see text for details). Particles are selected with: 0.3 GeV/c < p T < 1.5 GeV/c.
Figure 4. Σ [ N F , N B ] (11) as a function of the distance, Δ η , between forward and backward pseudorapidity acceptance intervals, where N F and N B multiplicities are calculated for inelastic pp interactions at s = 900 GeV (see text for details). Particles are selected with: 0.3 GeV/c < p T < 1.5 GeV/c.
Physics 05 00042 g004
Figure 5. Left: the variance, σ C 2 (24), of the asymmetry coefficient, C (23), distribution. Right: the Var ( N F N B ) versus its Skellam baseline as a function of the distance, Δ η , between forward and backward pseudorapidity acceptance intervals, where N F and N B multiplicities are calculated event-by-event for inelastic pp interactions at s = 900 GeV (see text for details). Particles are selected with 0.3 GeV/c < p T < 1.5 GeV/c.
Figure 5. Left: the variance, σ C 2 (24), of the asymmetry coefficient, C (23), distribution. Right: the Var ( N F N B ) versus its Skellam baseline as a function of the distance, Δ η , between forward and backward pseudorapidity acceptance intervals, where N F and N B multiplicities are calculated event-by-event for inelastic pp interactions at s = 900 GeV (see text for details). Particles are selected with 0.3 GeV/c < p T < 1.5 GeV/c.
Physics 05 00042 g005
Figure 6. Centrality dependence of Σ [ N F , N B ] (11) for 0–20%, 20–40%, and 40–60% multiplicity-based centrality classes as a function of the distance Δ η between forward and backward pseudorapidity acceptance intervals, where N F and N B multiplicities are calculated for inelastic pp interactions at s = 900 GeV in the model configuration for interacting strings without transverse dynamics (see text for details). Particles are selected with 0.3 GeV/c < p T < 1.5 GeV/c.
Figure 6. Centrality dependence of Σ [ N F , N B ] (11) for 0–20%, 20–40%, and 40–60% multiplicity-based centrality classes as a function of the distance Δ η between forward and backward pseudorapidity acceptance intervals, where N F and N B multiplicities are calculated for inelastic pp interactions at s = 900 GeV in the model configuration for interacting strings without transverse dynamics (see text for details). Particles are selected with 0.3 GeV/c < p T < 1.5 GeV/c.
Physics 05 00042 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Prokhorova, D.; Andronov, E.; Feofilov, G. Interacting Colour Strings Approach in Modelling of Rapidity Correlations. Physics 2023, 5, 636-654. https://doi.org/10.3390/physics5020042

AMA Style

Prokhorova D, Andronov E, Feofilov G. Interacting Colour Strings Approach in Modelling of Rapidity Correlations. Physics. 2023; 5(2):636-654. https://doi.org/10.3390/physics5020042

Chicago/Turabian Style

Prokhorova, Daria, Evgeny Andronov, and Grigory Feofilov. 2023. "Interacting Colour Strings Approach in Modelling of Rapidity Correlations" Physics 5, no. 2: 636-654. https://doi.org/10.3390/physics5020042

Article Metrics

Back to TopTop