Strangeness production in high-multiplicity events

Marat Siddikov, Iván Schmidt Departamento de Física, Universidad Técnica Federico Santa María, y Centro Científico - Tecnológico de Valparaíso, Casilla 110-V, Valparaíso, Chile (Dated: January 5, 2021) In this paper we analyze in detail the production of strangeness in -proton collisions in the kinematics of large transverse momenta pT of produced . Using the color dipole framework, we estimated the production cross-sections for and demonstrated that the shapes of the pT - dependence are in agreement with available experimental data. We also analyzed the self-normalized yields of strange hadrons as a function of multiplicity of co-produced hadrons, and found that the predictions are in agreement with the faster-than-linear growth seen in experimental data. Our description is largely parameter-free and complements our previous studies dedicated to the explanation of multiplicity enhancement of quarkonia, as well as D- and B-. Keywords: CGC approach, strangeness production, multiplicity dependence.

I. INTRODUCTION

Since the early experiments at RHIC and SPS [1–6], the production of hadrons containing strange has been used as one of the probes of - Plasma (QGP) formation, described in the framework of QGP-inspired hydrodynamic models [7–21]. The QGP manifests itself in different ways in experimental observables. For example, it might lead to an enhancement of strange yields in heavy ion collisions, compared to pA and pp production of the same strange hadrons. Another possibility to observe the effects of QGP is via the enhancement of strangeness in events with large multiplicity of co-produced hadrons. This observable might be studied independently in pp, pA or AA collisions. While early experiments confirmed the enhancement of strangeness in heavy ion collisions, similar multiplicity enhancements has been recently observed at the LHC not only in heavy ion, but also in pA [22, 23] and even in pp collisions [24], where QGP formation in significant amounts is highly unlikely even at TeV-range collision energies. For this reason it makes sense to understand better the microscopic mechanisms of this phenomenon, at least in pp collisions. In general, application of perturbation theory for strangeness production is challenging due to lack of the hard scale (like e.g. heavy mass of the quark). Threfore the phenomenological description of strangeness production has been mostly limited to studies in the framework of Monte-Carlo generators [25, 26], and inevitably includes additional model-dependent assumptions. Recently detailed studies of heavy quarkonia [27, 28] and open heavy flavor mesons [29] in pp collisions have discovered that similar enhancement with multiplicity also happens for the production of heavier quarks, charm and bottom. This enhancement has a quite complicated dependence on the rapidity separation of the bins used to collect quarkonia and light , on the existence of rapidity gaps between the heavy hadrons and colliding [30], as well as (possibly) on the quantum numbers of the produced quarkonia states [31]. While the production of these mesons can be described in the two- fusion picture [32–42], as was pointed out in [25], the description of the multiplicity dependence presents challenges for the established two-pomeron paradigm. For this reason a number of new mechanisms have been suggested for its description: e.g. the percolation approach [43], a modification of the slope of the elastic amplitude [44] or contributions of multipomeron diagrams [45–50]. In view of the similarity of multiplicity enhancements observed in strange, charm and bottom sectors, it is very desirable to describe the phenomenon for all flavors in the same framework. Since the strange quarks have very arXiv:2101.00461v1 [hep-ph] 2 Jan 2021 light mass, in general it is very challenging to apply the theoretical tools which rely on the perturbative QCD or heavy quark mass limit for their justification. Nevertheless, in the kinematics of very large transverse momenta pT of produced strange hadrons, the latter variable effectively plays the role of hard scale which partially justifies the use of such perturbative tools. In what follows we will use the color dipole (CGC/Saturation) approach [51–59] which was previously applied to D- and B- production in [60–64]. The generalization of this framework to high-multiplicity events is well-known from the literature [64–71] and allows to explain the multiplicity dependence observed in both the charm and bottom sectors. Our analysis will be mostly focused on the production of kaons and Λ-, due to lack of information about fragmentation functions of other strange hadrons. The paper is structured as follows. In the next Section II we describe a framework for strangeness production in the CGC/Sat approach. In Section III we make numerical estimates for the cross-sections and compare with available experimental data. In Section IV we discuss the multiplicity dependence of strange hadrons in the large-pT kinematics and demonstrate that our approach can describe the experimentally observed dependence for kaons and Λ-baryons. Finally, in Section V we draw conclusions. 2

p p p p

p p p p

Figure 1: Left plot: The leading order two-pomeron mechanism of quark pair production. The diagram includes two cut (upper and lower gluon ladders). Right plot: Example of 3-pomeron mechanisms which might be relevant in the small-pT kinematics (additional reggeons shown with gray color, not considered in this paper for the reasons discussed in the text). In all plots the vertical dashed line stands for unitarity cuts. A summation over all possible permutations of gluon vertices in the quark line (loop) is implied.

II. PRODUCTION OF STRANGE HADRONS VIA FRAGMENTATION

We assume that all strange hadrons are produced via a fragmentation mechanism, and we will perform our eval- uations within the framework developed earlier in [60–64]. In this approach the cross-section is related to the quark pair QQ¯ production cross-section by

1   dσpp→M+ X dz xQ(y) dσpp→Q¯iQi+X = Di (1) dy d2p ˆ z2 z dy∗d2p∗ T i xQ T

∗ where y is the rapidity of the produced strange , y = y ln z is the rapidity of the quark, pT is the transverse − momentum of the produced strange hadron, Di(z) is the fragmentation function which describes the formation of a ∗ given final state from a parton of flavor i, and dσpp→Q¯iQi+X /dy is the cross-section of quark pair production with ∗ ∗ quark rapidity y and transverse momentum pT = pT /z. For the fragmentation functions of kaons and Λ-baryons we will use the expressions available from the literature (see the Appendix A for details). Up to the best of our knowledge, currently there is no data for the fragmentation functions for Ω, Ξ-baryons and φ-mesons, for this reason we will not consider them in what follows. Naturally, the dominant contribution in the strange sector stems from the strange quarks, although there are also contributions from other flavors. In what follows we will focus on the ∗ 2 ∗ evaluation of the cross-section dσpp→Q¯iQi+X /dy d pT which appears in the integrand of (1). In high energy kinematics, the inclusive production gets its dominant contribution from the fusion of two pomerons, which for heavy quarkonia production is given by the diagram shown in the left panel of Figure 1. In the rest frame of one of the protons, this process might be viewed as a fluctuation of the incoming virtual gluon into a heavy QQ¯ pair, with subsequent scattering of the QQ¯ dipole on the target proton. In the kinematics of LHC experiments the average light-cone momentum fractions x1,2 carried by are very small ( 1), and the gluon densities are enhanced. This enhancement implies that there could be sizable corrections from multiple pomeron exchanges between the heavy dipole and the target, which are formally suppressed for small dipoles. For this reason, instead of a hard process on individual partons it is more appropriate to use the color dipole framework (also known as CGC/Sat) [51–59]. At high energies, the color dipoles are eigenstates of interaction, and therefore they can be used as universal elementary building blocks, automatically accumulating both the hard and soft fluctuations [72]. In fact, the light-cone color dipole framework has been successfully applied to phenomenological descriptions of both hadron-hadron and - hadron collisions [73–80]. Another advantage of the CGC/Sat framework is that it allows a relatively straightforward extension for the description of high-multiplicity events, as discussed in [64–71]. 3

In the dipole approach, the quark production cross-section is given by [63, 64]

1 1 ¯ dσpp→QiQi+X (y, √s) 2 0 2 = d kT x1 g (x1, pT kT ) dz dz (2) dy d pT ˆ − ˆ0 ˆ0 2 2 d r1 d r2 i(r1−r2)·kT † † e Ψ (r2, z, pT )Ψ (r1, z, pT ) × ˆ 4π ˆ 4π QQ¯ QQ¯ NM (x2(y); ~r1, ~r2) , × p 2 2 mM + p⊥M ±y x1,2 h ie (3) ≈ √s where y and pT are the rapidity and transverse momenta of the produced in the center-of-mass frame of the colliding protons; kT is the transverse momentum of the strange quark with respect to incident gluon; g (x1, pT ) ¯ in the first line of (2) is the unintegrated gluon PDF; Ψg→QQ¯ (r, z) is the light-cone wave function of the QQ pair with transverse separation between quarks r and the light-cone fraction carried by the quark z. In general this is a nonperturbative object, and there is no model-independent way to evaluate it [101]. For this reason, in what follows we will restrict our consideration to the kinematics of large transverse momenta of produced hadrons, which in view of (2) implies that typical sizes of the dipoles are also small, 1/pT . In this kinematics we may use standard perturbative expressions [81, 82] ∼

† 2 2 αsNc  2  iθ12 2 −iθ12 2 Ψ r2, z, Q ΨT r1, z, Q =  K1 (f r1) K1 (f r2) e z + e (1 z) (4) T 2π2 f − 2 +mf K0 (f r1) K0 (f r2) ,

† 2 2 αsNc  2 2 2 Ψ r2, z, Q ΨL r1, z, Q = 4Q z (1 z) K0 (f r1) K0 (f r2) , (5) L 2π2 −

2 = z (1 z) Q2 + m2 (6) f − f

2 2 2 (f) 2 (f) 2 (f) 2 Ψ r, z, Q = ΨT r, z, Q + ΨL r, z, Q (7)

The meson production amplitude NM depends on the mechanism of QQ¯ pair formation. For the case of the two- pomeron fusion, in leading order it is given by [60, 63]

1 1 9 NM (x, ~r1, ~r2) = N (x, ~r1 ~r2) [N (x, ~r1) + N (x, ~r2)] N (x, z¯(~r1 ~r2)) (8) −2 − − 16 − 8 − 9 + [N (x, z~r¯ 1 ~r2) + N (x, z~r¯ 2 ~r1) + N (x, z~r¯ 1) + N (x, z~r¯ 2)] , 16 − − where N (x, ~r) is the color singlet dipole scattering amplitude. In the LHC kinematics at large transverse momenta (our principal interest) the natural choice of the saturation scale is µF pT , which significantly exceeds the saturation ∼ scale Qs(x). This finding justifies the use of two-pomeron approximation. However, in the kinematics of smaller pT , potentially there could be multipomeron contributions, like those shown in the right panel of Figure 1. We will not consider such contributions since, as we mentioned earlier, we can’t describe the small-pT region, due to lack of the nonperturbative wave function ΨQQ¯ . The unintegrated gluon PDF, which appears in the prefactor of (2), can be related to the integrated PDF xG (x, µF ) as [83]

2 ∂ x g x, k = 2 xG (x, µF ) , (9) ∂µ 2 2 F µF =k and the latter is closely related to the dipole scattering amplitude N (y, r) = d2b N (y, r, b) via a set of identities [65, 84] ´

2 ! µF 2 CF d kT i~kT ·~r d kT 2 N (y, ~r) = 4 φ (y, kT ) 1 e ; x G (x, µF ) = 2 φ (x, kT ) , (10) 2π α¯S ˆ kT − ˆ0 kT 4 where y = ln(1/x). The Eq. (10) can be inverted and gives the gluon uPDF in terms of the dipole amplitude,

CF µF 2 J1 (r µF ) 2 xG (x, µF ) = 2 d r rN (y, ~r) . (11) 2π α¯S ˆ r ∇ This allows to rewrite the result in a symmetric and self-consistent form, which allows straightforward generalization for high-multiplicity events.

III. NUMERICAL RESULTS

For the sake of definiteness, in our numerical evaluations we will take the “CGC” parametrization of the dipole cross-section [85–87],

 2γeff (r)  r Qs(x)  2  N0 , r N (x, ~r) = σ0 2 Qs(x) , (12) ≤ 2 ×  1 exp ( ln ( r Qs)) , r > − −A B Qs(x) 2 2 1−N N0 γs 1 − 0 N0γs = 2 , = (1 N0) , (13) A −(1 N0) ln (1 N0) B 2 − − − λ/2   x0  1 2 Qs(x) = , γeff (r) = γs + ln , (14) x κλY r Qs(x) −5 γs = 0.762, λ = 0.2319, σ0 = 21.85 mb, x0 = 6.2 10 (15) × 0 which is widely used in the literature, and in what follows we will also focus on KS-meson production. The fragmen- ± 0 tation functions for K and KS are constrained by the isospin symmetry relation [88]

0 ± KS 2 1 K 2 D z, µ = D¯ z, µ , (16) i 2 i

0 ± and for this reason the cross-sections for KS and K production are proportional to each other with very good precision. For Λ-baryons we found only one parametrization for the fragmentation function only [88], and as shown in Appendix (A), at large-z its strange flavor component is approximately proportional to that of kaons; therefore the Λ- cross-sections are proportional to that of kaons in the large-pT kinematics. 0 In the left panel of Figure 2 we show the results for KS production cross-section, in the kinematics of ongoing and planned experiments. As can be seen in the literature, all the experimental papers contain results for the self- −1 normalized yields Nev dN/dy dpT , instead of cross-sections. The normalization parameter Nev is chosen as the total number of events Nev dpT dN/dpT or the number of non-single-diffractive events NNSD. The self-normalized yields are difficult to describe∼ ´ in our approach, because the normalization coefficient gets its dominant contribution from the small-pT region and thus cannot be evaluated reliably. This introduces a normalization uncertainty in our evaluations of such self-normalized yields. In Figure 3 we show the calculated yields in comparison with available experimental data from ALICE [89], CMS [90], CDF [91] and STAR [92] collaborations [102]. The STAR data [92] were included by authors of [88] in their global fit of fragmentation functions, for this reason the description of these data is nearly perfect. For data from LHC and the model provides a very reasonable description of the shape, although, as expected, there is a mismatch in the normalization by a factor of two. The theoretical shape of the pT -dependence starts deviating from the experimental data in the region pT .2 GeV, where nonperturbative effects become pronounced.

IV. MULTIPLICITY DEPENDENCE

Since the dipole approach (2, 8, 12) provides a reasonable description of the strangeness production in the large-pT kinematics, we can apply it to the study of the dependence on the number of charged particles Nch co-produced together with a given strange hadron. The extension of the color dipole (CGC/Sat) framework to the description of high-multiplicity events is quite straightforward, as was discussed in [30, 45, 60, 64–71]. In what follows we will briefly summarize the main results which will be used for the phenomenological estimates of the multiplicity dependence. In view of the Local Parton Hadron Duality (LPHD) hypothesis [93–95], the multiplicity of produced hadrons in a given event is directly proportional to the number of partons produced in a collision. For this reason, the study of high 5

101

= 100 s 7 TeV s=13 TeV 10-1 s=100 TeV , μ b / GeV T p

ⅆ -2

y 10

ⅆ σ / 2 -3 ⅆ 10 0 KS 10-4 5 10 20

pT , GeV

0 Figure 2: The pT -dependence of the cross-section dσ/dpT for KS -meson production at central rapidities, evaluated with the two-pomeron fusion mechanism. Results for Λ have similar shape, due to the similarity of their fragmentation functions (see details in Appendix A). multiplicity events in different channels allows to understand better the onset of the saturation regime in high energy collisions. The probability of multiplicity fluctuations decreases rapidly as function of the number of produced charged particles Nch [96], for this reason in the study of the multiplicity dependence it is more common to use a self-normalized ratio [89]

2  dNM /dy dσM (y, η, √s, n) /dy dσch η, √s, Q , n /dη = / 2 (17) dNM /dy dσM (y, η, √s, n = 1) /dy dσch (η, √s, Q , n = 1) /dη h i h i h i where n = Nch/ Nch is the relative enhancement of the number of charged particles in a given observation window h i (e.g. pseudorapidity bin (η ∆η/2, η+∆η/2)); dσM (y, √s, n) is the strange hadron M production cross-section, with − rapidity y and Nch = n Nch charged particles; dσch(y, √s, n) is the total production cross-section for Nch = n Nch charged particles in theh samei observation window. Since the cross-sections are proportional to the probabilityh toi produce a given final state, the ratio (17) might be interpreted as a conditional probability to produce a strange hadron M in a pp collision in which Nch charged particles are produced. We expect that even in high multiplicity events each gluon cascade (“pomeron”) should satisfy the nonlinear Balitsky- Kovchegov equation, and therefore the dipole amplitude (12) should keep its form, although the value of the saturation scale Qs might be modified. As was demonstrated in [65–67], the observed number of charged multiplicity dNch/dy of 2 soft hadrons in pp collisions is proportional to the saturation scale Qs (modulo logarithmic corrections), and therefore 2 in the dipole framework the events with large multiplicity might be described by simply rescaling Qs as a function of n [65–71],

Q2 (x, b; n) n Q2 (x, b) . (18) s ≈ The accuracy of the approximation (18) was tested in [64], and it was found that its error does not exceed 10 per cent in the region of interest (n . 10), on par with the precision of current evaluations. Therefore, in what follows we 6

100 100 ALICE, pp, s=13 TeV CMS, pp, s=7 TeV - 1 - 1 10-1 10-1 1 ⅆ2N , GeV , GeV

T ⅆ ⅆ T Nev y pT 2 p -2 p -2 1 ⅆ N 10 2 10 ⅆ 1 ⅆ N ⅆ Nev ⅆ y ⅆ pT y λinc⨯ y Nev ⅆ y ⅆ pT ⅆ ⅆ 1 ⅆ2N -3 -3 λ ⨯ N / 10 N / 10 inc ⅆ ⅆ 2 2 Nev y pT ) ⅆ ) ⅆ ev ev 10-4 |y|<0.5 10-4 |y|<0.5

( 1 / N 0 ( 1 / N 0 10-5 KS 10-5 KS 1 5 10 1 5 10

pT , GeV pT , GeV

100 100 CDF, pp, s=1.8 TeV STAR, pp, s=200 GeV - 1 - 1 10-1 10-1

, GeV , GeV 2 T 2 T 1 ⅆ N p -2 1 ⅆ N p -2 10 10 ⅆ ⅆ ⅆ ⅆ Nev y pT Nev ⅆ y ⅆ pT y y ⅆ 1 ⅆ2N ⅆ -3 λ ⨯ -3 N / 10 inc ⅆ ⅆ N / 10 2 Nev y pT 2 ) ⅆ ) ⅆ ev ev 10-4 |y|<0.5 10-4 |y|<0.5

( 1 / N 0 ( 1 / N 0 10-5 KS 10-5 KS 1 5 10 1 2 3 4 5 6

pT , GeV pT , GeV

0 Figure 3: Self-normalized pT -dependence of the KS -meson yields at central rapidities. Theoretical predictions (solid line) are compared with experimental data from ALICE [89], CMS [90], CDF [91] and STAR [92] collaborations. As we explained in the text, our approach is not very reliable at small pT , and thus the evaluation of the global normalization factor Nev might have large nonperturbative corrections. In order to demonstrate that the description of the shape is correct, we have plotted also the yields multiplied by a constant factor λinc ≈ 1.7 − 2.2 (dashed lines). will use (18) for our estimates. While at LHC energies it is expected that the typical values of the saturation scale Qs (x, b) fall into the range 0.5-1 GeV, from (18), we can see that in events with enhanced multiplicity this parameter might lead to an interplay of the large-Qs and large-pT limits. Since increasing multiplicity and increasing energy 2 (decreasing x) affect Qs in a similar way, the study of the high-multiplicity events allows to study a deeply saturated regime, which determines the dynamics of all processes at significantly higher energies. Since at high energies each pomeron hadronizes independently [97], the observed enhancement of multiplicity in the whole process must be shared between all pomerons which might contribute in a given rapidity window. For this reason, for phenomenological estimates of the multiplicity dependence it is important if the rapidity bin used to collect charged particles Nch overlaps with the bin used for the observation of strange large-pT hadrons, as explained in Figure 4. In the limit of large-pT the typical sizes of the dipoles are small, r 1/pT , so we may expect from (2, 12, 18) ∼ hγeff i that the contribution of each cut pomeron to the multiplicity dependence is given by the factor n , where ni is ∼ i the relative enhancement of multiplicity assigned to a given pomeron, and the value of the parameter γeff is given in (14, 15). For the configuration when the bins used to collect charged and strange particles are separated by rapidity (see the left panel of the Figure 4) the multiplicity is assigned to one of the pomerons, so the expected multiplicity dependence of the cross-section is nhγeff i. For the case when the strange and bins partially overlap (as shown in the right panel of the∼ Figure 4), we should average over all possible partitions of the observed number of charged particles. This evaluation technically is quite complicated, although we can assume with good precision that the multiplicity enhancement is shared equally between both pomerons [46]. For this reason in this case we expect that the multiplicity cross-section would be (n/2)2hγeff i. ∼ 7

p p

∆y ∆y (y , y + ) (y ∆y , y + ∆y ) − 2 2 − 2 2

∆η ∆η (η ,η + ) ∆η ∆η 2 2 (η ,η + ) − − 2 2

p p

Figure 4: (color online) Demonstration that sharing of the enhanced multiplicity of the whole process between individual pomerons depends on the experimental setup. Left plot: The experimental setup in which the rapidity bin used for collection of strange hadrons (blue box) does not overlap with the bin used for the collection of charged particles (red box). The elevated multiplicity in this case should be unambiguously attributed to the lower pomeron. Right plot: The experimental setup when the bins partially overlap. For the partons in the intersection region (magenta color) the assignment to upper or lower pomerons depends on the position of large-pT quark inside the bin (y − ∆y/2, y + ∆y/2). In the final result we should average over all possible rapidities of strange quark inside the bin.

Currently the data on multiplicity dependence of strange hadrons are available from the ALICE experiment [89]. As we can see from Figure 5, the theoretical curves can describe reasonably well the slope of the experimentally observed n-dependence (in logarithmic coordinates), although apparently overestimate all the experimental points by the same normalization factor 1.2. We would like to stress that by definition at the point n = 1 the self-normalized ratio (17) equals one. This condition∼ is fulfilled in our theoretical curves, and therefore we believe that the normalization of our curves is correct.

V. CONCLUSIONS

In this paper we studied the production of strange hadrons in the color dipole approach. We found that the CGC/Sat approach can describe the shapes of pT distributions in the large-pT kinematics, although it might not be very reliable for smaller pT . The latter restriction implies that the suggested approach cannot be applied to pT -integrated observables, which get its dominant contribution from the nonperturbative small-pT region. As a consequence, our predictions for experimentally measurable self-normalized yields suffer from a global normalization uncertainty, which complicates the direct comparison of model predictions with data, even for large pT . For this reason we call experimentalists to publish also the cross-sections, for which there is unambiguous separation of small-pT and large-pT . We made predictions for such cross-sections in the kinematics of ongoing and future experiments, and can provide further predictions on demand. We also applied the CGC/Sat approach to the description of the multiplicity dependence measured by ALICE [89]. Fortunately, all kaons and Λ-baryons were collected with sufficiently large transverse momenta pT & 4 GeV, where our approach is well justified. We found that the theoretical predictions are in reasonable agreement with experimental data. Our evaluation is largely parameter-free and relies only on the choice of the parametrization for the dipole cross-section (12) and fragmentation functions of strange hadrons. This study complements our previous analysis of the multiplicity dependence of heavier charm and bottom pro- duction [60] and demonstrates that at sufficiently large pT it is possible to describe all of them within the same framework. 8

ALICE, pp, s=13 TeV 10 0 4 12 GeV KS (

Λ+Λ (4

K K 1 dN  dN

0.1 0.1 0.5 1 5 10

dNch /dη

〈dNch /dη〉

0 Figure 5: Comparison of the theoretical multiplicity dependence for KS meson (production solid curve) and Λ baryons (dashed curve) with experimental data from ALICE [89]. For the sake of reference we have also shown a dotted line, which corresponds to a linear dependence. The charged particles and strange hadrons are collected at central rapidities.

Acknowledgements

We thank our colleagues at UTFSM University for encouraging discussions. This research was partially supported by Proyecto Basal FB 0821(Chile) and Fondecyt (Chile) grant 1180232. Also, we thank Yuri Ivanov for technical support of the USM HPC cluster, where part of the evaluations were performed.

Appendix A: Fragmentation functions

In this section we would like to summarize briefly the fragmentation functions used in our evaluations. These functions are nonperturbative objects, which cannot be evaluated from first principles. For this reason currently their parametrization is extracted from the phenomenological fits of experimental data. For the sake of definiteness, for our evaluations we used the fragmentation functions for kaons and Λ from [88] (so-called AKK08 parametrization). ± 0 The fragmentation functions for K and KS are constrained by the isospin symmetry relation

0 ± KS 2 1 K 2 D z, µ = D¯ z, µ , (A1) i 2 i 0 therefore in what follows we will consider only the fragmentation function of neutral kaons KS. For kaons we checked that the alternative parametrizations of fragmentation functions DSS17 [98], NNPDF [99] and JAM [100] give similar results in the region of interest. We have not found parametrizations for fragmentation functions of strange baryons 0∗ Ω, Ξ, and neither of KS and φ-mesons, and for this reason we do not consider them in this paper. 9

4 K ±, AKK08 0.08 3 Λ, AKK08 0.06 ± ( z )

( z ) K , DSS17 2 s / H 0.04 s / H D D 3

1 z 0.02

0 0. 0. 0.2 0.4 0.6 0.8 1. 0. 0.2 0.4 0.6 0.8 1. z z

Figure 6: Left: The fragmentation function Ds/H of kaons and Λ-baryon (s-quark component), evaluated in AKK08 [88] and DSS17 [98] parametrizations. Right: The same function multiplied by z3. As explained in the text, for the physical cross-section 3 in the large-pT kinematics the fragmentation function contributes multiplied by an additional factor ∼ z , so the difference of the fragmentation functions in the small-z domain has a minor effect on physical observables.

In the AKK08 parametrization [88] it is assumed that the fragmentation function is given by

i/H ai bi  di  D (z) = Niz (1 z) 1 + ci(1 z) , (A2) − − where Ni, ai, bi, ci, di are some numerical coefficients which depend on the hadron and quark flavor i. We expect that for strange hadrons the largest contribution comes from the fragmentation of the strange quark, thus we will discuss below the fragmentation function Ds/H . As we can see from Figure 6, the parametrizations for K± mesons and Λ-baryons differ quite substantially in the region of small z . 0.3, although become comparable for all hadrons at larger values of z. In this paper we are mostly interested in the large-pT kinematics, and it is possible to show that this region has stronger sensitivity to the region of large z. Indeed, as we can see from the structure of (2), at large pT the cross- n QQ¯  QQ¯  QQ¯ section dσ ¯ /dp is suppressed as 1/p with n 5. The momentum of the quark pair p is related to the QQ T ∼ T & T QQ¯ momentum of the strange hadron as pT = pT /z, so this implies that in the integral over the fragmentation fraction z in (1) effectively we get an additional prefactor zn−2, which suppresses the contribution of the small-z domain. As we can see from the right panel of Figure 6, the∼ dominant contribution comes from the region z 0.6 0.8, where the difference between fragmentation functions does not exceed a factor of two. ∼ −

[1] T. Alber, et al. Z. Phys. C 64 (1994), 195. [2] E. Anderson, et al. Phys. Lett. B294 (1992), 127; Phys. Lett. B 327 (1994), 433. [3] S. Abatzis S, et al. Nucl. Phys. A566 (1994), 499c; Phys. Lett. B316 (1993), 615. [4] D. DiBari et al. Nucl. Phys. A590 (1995) 307c. [5] J.B. Kinson, et al. Nucl. Phys. A590 (1995) 317c. [6] Gazdzicki M, et al. Nucl. Phys. A590 (1995), 197c. [7] J. Rafelski, B. Müller, Phys. Rev. Lett. 48 (1982), 1066 [Erratum: ibid. 56 (1986), 2334]. [8] J. Rafelski, Phys. Rep. 88 (1982), 331. [9] P. Koch, B. Müller, J. Rafelski. Phys. Rep. 142 (1986), 167. [10] J. Letessier, A. Tounsi, U. W. Heinz, J. Sollfrank and J. Rafelski, Phys. Rev. Lett. 70 (1993), 3530. [11] J. Letessier, J. Rafelski, A. Tounsi, Phys. Lett. B328 (1994), 499. [12] J. Rafelski, M. Danos, Phys. Rev. C50 (1994), 1684. [13] Sollfrank J, Gazdzicki M, Heinz U, Rafelski J. Z. Phys. C61 (1994), 659. [14] J. Letessier, A. Tounsi, U. Heinz, J. Sollfrank, J. Rafelski, Phys. Rev. D 51 (1995), 3408. [15] P. Braun-Munzinger, J. Stachel, J.P. Wessels, N. Xu, Phys. Lett. B344 (1995), 43. [16] J. W. Harris and B. Muller, “The Search for the quark - gluon plasma,” Ann. Rev. Nucl. Part. Sci. 46 (1996), 71-107 [arXiv:hep-ph/9602235 [hep-ph]]. 10

[17] P. Koch, B. Muller and J. Rafelski, “From strangeness enhancement to quark-gluon plasma discovery,” Int. J. Mod. Phys. A 32 (2017) no.31, 1730024 [arXiv:1708.08115 [nucl-th]]. [18] J. Rafelski and B. Muller, “Strangeness Production in the Quark - Gluon Plasma,” Phys. Rev. Lett. 48 (1982), 1066 [erratum: Phys. Rev. Lett. 56 (1986), 2334]. [19] J. Rafelski, “Discovery of Quark-Gluon-Plasma: Strangeness Diaries,” Eur. Phys. J. ST 229 (2020) no.1, 1-140 [arXiv:1911.00831 [hep-ph]]. [20] I. Schmidt, M. Siddikov and M. Musakhanov, “Charmonia suppression in heavy ion collisions: Hot versus cold effects,” Phys. Rev. C 98, no.2, 025207 (2018). [21] B. Z. Kopeliovich, I. Schmidt and M. Siddikov, “Suppression versus enhancement of heavy quarkonia in pA collisions,” Phys. Rev. C 95, no.6, 065203 (2017) [arXiv:1701.07134 [hep-ph]]. [22] Abelev, B. et al. (ALICE Collaboration), Phys. Lett. B 728, 25–38 (2014). [23] Adam, J. et al. (ALICE Collaboration), Phys. Lett. B 758, 389–401 (2016). [24] J. Adam et al. [ALICE Collaboration], Nature Phys. 13, 535 (2017) [arXiv:1606.07424 [nucl-ex]]. [25] N. Fischer and T. Sjöstrand, JHEP 1701, 140 (2017) [arXiv:1610.09818 [hep-ph]]. [26] H. J. Pirner, B. Z. Kopeliovich and K. Reygers, “Strangeness Enhancement due to Fluctuations,” Phys. Rev. D 101 (2020) no.11, 114010 [arXiv:1810.04736 [hep-ph]]. [27] D. Thakur [ALICE Collaboration], “J/ψ production as a function of charged-particle multiplicity with ALICE at the LHC,” arXiv:1811.01535 [hep-ex]. [28] B.√ Abelev et al. [ALICE Collaboration], “J/ψ production as a function of charged particle multiplicity in pp collisions at s=7 TeV”, Phys. Lett. B 712 (2012), 165. [29] J. Adam et al. [ALICE Collaboration], JHEP 1509, 148 (2015) [arXiv:1505.00664 [nucl-ex]]. [30] M. Siddikov and I. Schmidt, “Single diffractive production of open heavy flavor mesons,” Phys. Rev. D 102 (2020) no.7, 076020 [arXiv:2008.12446 [hep-ph]]. [31] M. Siddikov and I. Schmidt, “Multiplicity dependence of χc and χb meson production,” [arXiv:2012.08284 [hep-ph]]. [32] G. T. Bodwin, E. Braaten and G. P. Lepage, Phys. Rev. D 51, 1125 (1995) Erratum: [Phys. Rev. D 55, 5853 (1997)] [hep-ph/9407339]. [33] F. Maltoni, M. L. Mangano and A. Petrelli, Nucl. Phys. B 519, 361 (1998) [hep-ph/9708349]. [34] N. Brambilla, E. Mereghetti and A. Vairo, Phys. Rev. D 79, 074002 (2009) Erratum: [Phys. Rev. D 83, 079904 (2011)] [arXiv:0810.2259 [hep-ph]]. [35] Y. Feng, J. P. Lansberg and J. X. Wang, Eur. Phys. J. C 75, no. 7, 313 (2015) [arXiv:1504.00317 [hep-ph]]. [36] N. Brambilla et al., Eur. Phys. J. C 71, 1534 (2011) [arXiv:1010.5827 [hep-ph]]. [37] S. P. Baranov, A. V. Lipatov and N. P. Zotov, Eur. Phys. J. C 75, no. 9, 455 (2015) [arXiv:1508.05480 [hep-ph]]. [38] S. P. Baranov and A. V. Lipatov, Phys. Rev. D 96, no. 3, 034019 (2017) [arXiv:1611.10141 [hep-ph]]. [39] R. Baier and R. Ruckl, Phys. Lett. 102B (1981) 364. [40] E. L. Berger and D. L. Jones, Phys. Rev. D 23, 1521 (1981). [41] C. H. Chang, Nucl. Phys. B 172, 425 (1980). [42] R. Maciula and A. Szczurek, Phys. Rev. D 87, no. 9, 094022 (2013) [arXiv:1301.3033 [hep-ph]]. [43] E.G.Ferreiro and C.Pajares, Phys.Rev.C86,034903(2012)[arXiv:1203.5936[hep-ph]]. [44] B. Z. Kopeliovich, H. J. Pirner, I. K. Potashnikova, K. Reygers and I. Schmidt, Phys. Rev. D 88, no. 11, 116002 (2013) [arXiv:1308.3638 [hep-ph]] . [45] M. Siddikov, E. Levin and I. Schmidt, Eur. Phys. J. C 80 (2020), 560 [arXiv:1910.13579 [hep-ph]]. [46] E. Levin and M. Siddikov, “J/ψ production in hadron scattering: three-pomeron contribution,” Eur. Phys. J. C 79 (2019) no.5, 376 [arXiv:1812.06783 [hep-ph]]. [47] V. A. Khoze, A. D. Martin, M. G. Ryskin and W. J. Stirling, “Inelastic J/ψ and υ hadroproduction,” Eur. Phys. J. C 39, 163 (2005), [hep-ph/0410020]. [48] L. Motyka and M. Sadzikowski, “On relevance of triple gluon fusion in J/ψ hadroproduction,” Eur. Phys. J. C 75 (2015) no.5, [arXiv:1501.04915 [hep-ph]]. [49] I. Schmidt and M. Siddikov, “Contribution of digluons to charmonia production,” J. Phys. G 46 (2019) no.6, 065002 [arXiv:1801.09974 [hep-ph]]. [50] E. Levin, L. Lipatov and M. Siddikov, “BFKL pomeron with massive gluons and running coupling,” Phys. Rev. D 94, no.9, 096004 (2016) [arXiv:1608.03816 [hep-ph]]. [51] L. V. Gribov, E. M. Levin and M. G. Ryskin, “Semihard processes in QCD", Phys. Rep. 100 (1983) 1. [52] L. D. McLerran and R. Venugopalan, Phys. Rev. D 49, 2233 (1994) [hep-ph/9309289]. [53] L. D. McLerran and R. Venugopalan, Phys. Rev. D 49, 3352 (1994) [hep-ph/9311205]. [54] L. D. McLerran and R. Venugopalan, Phys. Rev. D 50, 2225 (1994) [hep-ph/9402335]. [55] A. H. Mueller and J. Qiu, Nucl. “Gluon recombination and shadowing at small values of x",Phys. B268 (1986) 427 [56] L. McLerran and R. Venugopalan,“Gluon distribution functions for very large nuclei at small transverse momentum", Phys. Rev. D49 (1994) 3352;‘Green’s function in the color field of a large nucleus" D50 (1994) 2225;“ Fock space distributions, structure functions, higher twists, and small x" , D59 (1999) 09400. [57] K. J. Golec-Biernat and M. Wusthoff, Phys. Rev. D 60, 114023 (1999) [hep-ph/9903358]. [58] B. Z. Kopeliovich and A. V. Tarasov, Nucl. Phys. A 710, 180 (2002) [hep-ph/0205151]. [59] B. Kopeliovich, A. Tarasov and J. Hufner, Nucl. Phys. A 696, 669 (2001) [hep-ph/0104256]. [60] I. Schmidt and M. Siddikov, Phys. Rev. D 101 (2020) no.9, 094020 [arXiv:2003.13768 [hep-ph]]. [61] J. Binnewies, B. A. Kniehl and G. Kramer, Phys. Rev. D 58, 034016 (1998) [hep-ph/9802231]. 11

[62] B. A. Kniehl and G. Kramer, Phys. Rev. D 60, 014006 (1999) [hep-ph/9901348]. [63] V. P. Goncalves, B. Kopeliovich, J. Nemchik, R. Pasechnik and I. Potashnikova, Phys. Rev. D 96, no. 1, 014010 (2017) [arXiv:1704.04699 [hep-ph]]. [64] Y. Q. Ma, P. Tribedy, R. Venugopalan and K. Watanabe, Phys. Rev. D 98, no. 7, 074025 (2018) [arXiv:1803.11093 [hep-ph]]. [65] Y. V. Kovchegov and E. Levin, at high energy Vol. 33 (Cambridge University Press, 2012). [66] D. Kharzeev and M. Nardi, “Hadron production in nuclear collisions at RHIC and high density QCD,” Phys. Lett. B 507 (2001) 121; [nucl-th/0012025]. D. Kharzeev and E. Levin, ‘ ‘Manifestations of high density QCD in the first RHIC data,” Phys. Lett. B 523 (2001) 79, [nucl-th/0108006]; D. Kharzeev, E. Levin and M. Nardi, “The Onset of classical QCD dynamics in relativistic heavy ion collisions,” Phys. Rev. C 71 (2005) 054903, [hep-ph/0111315]; “Hadron multiplicities at the LHC,” J. Phys. G 35 (2008) no.5, 054001.38 [arXiv:0707.0811 [hep-ph]]. [67] A. Dumitru, D. E. Kharzeev, E. M. Levin and Y. Nara, “ “Gluon Saturation in pA Collisions at the LHC: KLN Model Predictions For Hadron Multiplicities,” Phys. Rev. C 85 (2012) 044920 [arXiv:1111.3031 [hep-ph]]. [68] D. Kharzeev and M. Nardi, Phys. Lett. B 507, 121 (2001) [nucl-th/0012025]. [69] Y. V. Kovchegov, Nucl. Phys. A 692, 557 (2001) [hep-ph/0011252]. [70] E. Levin and A. H. Rezaeian, “Gluon saturation and inclusive hadron production at LHC,” Phys. Rev. D 82 (2010) 014022, [arXiv:1005.0631 [hep-ph]]. [71] T. Lappi, Eur. Phys. J. C 71, 1699 (2011) [arXiv:1104.3725 [hep-ph]]. [72] N.N. Nikolaev, B.G. Zakharov; J. Exp. Theor. Phys. 78, 598 (1994). [73] Y. V. Kovchegov, “Small x F(2) structure function of a nucleus including multiple pomeron exchanges,” Phys. Rev. D 60 (1999), 034008 [arXiv:hep-ph/9901281 [hep-ph]]. [74] Y. V. Kovchegov and H. Weigert, “Triumvirate of Running Couplings in Small-x Evolution,” Nucl. Phys. A 784 (2007), 188-226 [arXiv:hep-ph/0609090 [hep-ph]]. [75] I. Balitsky and G. A. Chirilli, “Next-to-leading order evolution of color dipoles,” Phys. Rev. D 77 (2008), 014019 [arXiv:0710.4330 [hep-ph]]. [76] Y. V. Kovchegov and E. Levin, “Quantum chromodynamics at high energy,” Camb. Monogr. Part. Phys. Nucl. Phys. Cosmol. 33 (2012). [77] I. Balitsky, “Effective field theory for the small x evolution,” Phys. Lett. B 518 (2001), 235-242 [arXiv:hep-ph/0105334 [hep-ph]]. [78] F. Cougoulic and Y. V. Kovchegov, “Helicity-dependent generalization of the JIMWLK evolution,” Phys. Rev. D 100 (2019) no.11, 114020 [arXiv:1910.04268 [hep-ph]]. [79] C. A. Aidala, E. Aschenauer, F. Aslan, A. Bacchetta, I. Balitsky, S. Benic, S. Bhattacharya, M. Boglione, M. Burkardt, J. Cammarota, G. A. Chirilli, C. Cocuzza, A. Courtoy, D. de Florian, P. Di Nezza, A. Dumitru, S. Fucini, K. Fukushima, Y. Furletova, L. Gamberg, O. Garcia-Montero, F. Gelis, V. Guzey, Y. Hatta, F. Hautmann, T. J. Hobbs, T. Horn, E. Iancu, S. Joosten, Z. B. Kang, R. Kishore, Y. V. Kovchegov, P. Kroll, K. Kumerički, K. Kutak, T. Lappi, H. W. Lin, X. Liu, S. Liuti, C. Lorcé, H. Mäntysaari, C. Marquet, Y. Makris, K. A. Mamo, Y. Mehtar-Tani, A. Metz, Z. E. Meziani, G. A. Miller, J. Miller, A. Mukherjee, P. M. Nadolsky, F. I. Olness, B. Pasquini, B. Pire, C. Pisano, D. Pitonyak, M. V. Polyakov, A. Prokudin, J. W. Qiu, M. Radici, A. Rajan, S. Rajesh, M. Rinaldi, K. Roy, C. Royon, N. Sato, M. Schlegel, G. Schnell, P. Schweitzer, S. Scopetta, R. Seidl, K. Semenov-Tian-Shansky, A. Signori, D. Sokhan, A. M. Stasto, L. Szymanowski, A. Tarasov, D. Triantafyllopoulos, T. Ullrich, R. Venugopalan, I. Vitev, W. Vogel- sang, A. Vossen, B. T. Wang, S. Wallon, K. Watanabe, C. Weiss, B. W. Xiao, H. U. Yee and Y. Zhao, “Probing and Nuclei in High Energy Collisions,” [arXiv:2002.12333 [hep-ph]]. [80] Y. Q. Ma and R. Venugopalan, ‘‘Comprehensive Description of J/ψ Production in Proton-Proton Collisions at Collider Energies,” Phys. Rev. Lett. 113 (2014) no.19, 192301 [arXiv:1408.4075 [hep-ph]]. [81] H. G. Dosch, T. Gousset, G. Kulzinger and H. J. Pirner, “ leptoproduction and nonperturbative gluon fluctuations in QCD,” Phys. Rev. D 55 (1997), 2602-2615 [arXiv:hep-ph/9608203 [hep-ph]]. [82] J. D. Bjorken, J. B. Kogut and D. E. Soper, “ at Infinite Momentum: Scattering from an External ,” Phys. Rev. D 3 (1971), 1382. [83] M. A. Kimber, A. D. Martin and M. G. Ryskin, Phys. Rev. D 63, 114027 (2001), [hep-ph/0101348]. [84] R. S. Thorne, “Gluon distributions and fits using dipole cross-sections,” AIP Conf. Proc. 792 (2005) no.1, 324. [85] H. Kowalski and D. Teaney, Phys. Rev. D 68, 114005 (2003) [hep-ph/0304189]. [86] H. Kowalski, L. Motyka and G. Watt, Phys. Rev. D 74, 074016 (2006) [hep-ph/0606272]. [87] A. H. Rezaeian and I. Schmidt, Phys. Rev. D 88 (2013) 074016, [arXiv:1307.0825 [hep-ph]]. [88] S. Albino, B. A. Kniehl and G. Kramer, “AKK Update: Improvements from New Theoretical Input and Experimental Data,” Nucl. Phys. B 803 (2008), 42-104 [arXiv:0803.2768 [hep-ph]]. [89] S.√ Acharya et al. [ALICE], “Multiplicity dependence of (multi-)strange hadron production in proton-proton collisions at s = 13 TeV,” Eur. Phys. J. C 80, no.2, 167 (2020) [arXiv:1908.01861 [nucl-ex]].√ [90] V. Khachatryan et al. [CMS], ‘‘ Production in pp Collisions at s = 0.9 and 7 TeV,” JHEP 05, 064 (2011) [arXiv:1102.4282 [hep-ex]]. 0 0 √ [91] D. Acosta et al. [CDF], “KS and Λ production studies in pp¯ collisions at s = 1800-GeV and 630-GeV,” Phys. Rev. D 72 (2005), 052001 [arXiv:hep-ex/0504048 [hep-ex]]. [92] B. I. Abelev et al. [STAR], “Strange particle production in p+p collisions at s**(1/2) = 200-GeV,” Phys. Rev. C 75 (2007), 064901 [arXiv:nucl-ex/0607033 [nucl-ex]]. [93] Y. L. Dokshitzer, V. A. Khoze and S. I. Troian, “On the concept of local parton hadron duality,” J. Phys. G17 (1991) 12

1585. 36. [94] V. A. Khoze, W. Ochs and J. Wosiek, “Analytical QCD and multiparticle production,” arXiv:hep-ph/0009298. [95] V. A. Khoze and W. Ochs, “Perturbative QCD approach to multiparticle production,” Int. J. Mod. Phys. A12 (1997) 2949 [arXiv:hep-ph/9701421]. [96] B. Abelev et al. [ALICE Collaboration], Phys. Lett. B 712, 165 (2012) [arXiv:1202.2816 [hep-ex]]. [97] G. P. Korchemsky, J. Kotanski and A. N. Manashov, Phys. Rev. Lett. 88, 122002 (2002) [hep-ph/0111185]. [98] D. de Florian, M. Epele, R. J. Hernandez-Pinto, R. Sassot and M. Stratmann, “Parton-to- Fragmentation Revisited,” Phys. Rev. D 95 (2017) no.9, 094019 [arXiv:1702.06353 [hep-ph]]. [99] V. Bertone et al. [NNPDF], “A determination of the fragmentation functions of , kaons, and protons with faithful uncertainties,” Eur. Phys. J. C 77 (2017) no.8, 516 [arXiv:1706.07049 [hep-ph]]. [100] N. Sato, J. J. Ethier, W. Melnitchouk, M. Hirai, S. Kumano and A. Accardi, “First Monte Carlo analysis of fragmentation functions from single-inclusive e+e− ,” Phys. Rev. D 94 (2016) no.11, 114004 [arXiv:1609.00899 [hep-ph]]. [101] We would like to remind that in case of D- and B-meson production studied earlier in [60–64] the scale in the small-pT kinematics was set by the heavy quark mass, so the use of the framework was justified up to pT ≈ 0. In case of the strangeness production this is no longer true, and we have to restrict our consideration to the large-pT domain only. [102] We also would like to mention that the CMS and STAR data were not normalized to unity. For example, the CMS data instead of Nev used the number of non-single-diffractive events NNSD, so we also corrected the normalization of experimental data.