Next Article in Journal
Superdeterminism without Conspiracy
Next Article in Special Issue
Equation of State of Quark–Gluon Matter in the Clustering-of-Color-Sources Approach
Previous Article in Journal
Scaling Behaviour of dN/dy in High-Energy Collisions
Previous Article in Special Issue
Study of Transverse-Spherocity Biased pp Collisions at the LHC Energies Using the PYTHIA 8 Event Generator
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

String Interactions as a Source of Collective Behaviour

by
Christian Bierlich
Department of Physics, Lund University, P.O. Box 118, 221 000 Lund, Sweden
Universe 2024, 10(1), 46; https://doi.org/10.3390/universe10010046
Submission received: 8 December 2023 / Revised: 10 January 2024 / Accepted: 15 January 2024 / Published: 17 January 2024
(This article belongs to the Special Issue Relativistic Heavy Ion Collision)

Abstract

:
The discovery of collective effects in small collision systems has spurred a renewed interest in hadronization models, and is also a source for collective effects all the way to large collision systems, where they are usually ascribed to the creation of a Quark–Gluon Plasma. In this topical mini-review, the microscopic model for string interactions, based on the Lund string hadronization model, developed with exactly this aim in mind is reviewed, and some prospects for the future are presented.

1. Introduction

Collisions of heavy nuclei at high energy are most often analyzed under the assumption that a Quark–Gluon Plasma (QGP) is created fractions of a fm after the collision. As evidence for QGP creation, several experimental signatures are considered. These include observations1 of strangeness enhancement relative to the proton–proton baseline [1,2] and anisotropic flow [3,4]. Over the past decades, a rich research program has emerged to studying the properties of the QGP.
The QGP paradigm differs vastly from the dynamics traditionally assumed to govern high-energy collisions of protons. In the latter, models based on the formation of strings [5] or clusters [6], which subsequently decay into hadrons, have dominated the field. These models are implemented in the highly successful “General purpose Monte Carlo event generators” [7,8], such as Pythia [9], Herwig [10], and Sherpa [11]. In proton collisions, this includes the scattering of multiple partons in each collision event, which, in the case of strings, leads to multiple strings connecting the emerging partons with each other and the beam remnants. No formation of QGP has traditionally been assumed, and the models have done remarkably well in reproducing most features of  e + e , ep, and pp collisions.
Over the past decade, the distinction between small ( e + e , ep, and pp) and large (pA and  A A ) collision systems has become much less clear. Measurements from the LHC have shown that pp collisions exhibit similar qualitative features of strangeness enhancement [12] and flow [13], as observed in larger collision systems (as well as many more recent observations. See ref. [14] for a recent review of ALICE results). There are even indications that flow signals can potentially form in  e + e , as shown by re-analyzed data from ALEPH [15]. This raises the question of whether a QGP is also created in small collision systems, challenging the “general purpose” paradigm of the aforementioned event generators. This assumption is successfully pursued in, for example, core–corona models [16,17]. Another possibility is to consider whether dynamics not including QGP formation could be responsible for the same signatures. This latter possibility is the underlying assumption of the microscopic model of interacting Lund strings, which has been developed over the past years, and is the topic of this mini-review2.
The main aim of the model is to establish a QGP-free baseline for predictions and postdictions for experimental measurements in heavy ion collisions.
The mini-review is structured as follows. In Section 2, the underlying Angantyr model is briefly described, before diving into the string formalism in Section 3. Here, interactions of multi-string systems are treated through Section 3.2 before some concluding remarks and an outlooks are given in Section 4.

2. Multi-Parton Interactions and the Angantyr Model

Strings are produced by partons created in hard collisions or subsequently radiated off in the parton shower. For multi-parton interactions (MPIs) in proton collisions, the main model by Sjöstrand and van Zijl [23] treats all parton–parton scatterings as independent  2 2  processes, with a scattering cross section calculable in perturbative QCD, which has the approximate behaviour [9]:
d σ ^ d p 2 α s 2 ( p 2 ) p 4 .
The expression is evidently divergent as  p 0 , and as such, it will at some point exceed the total nucleon-nucleon cross section, which is clearly unphysical. Even a cut-off scale of 1 GeV, a reasonable value for when one would expect perturbative QCD to break down, is too small. Instead, a screening scale  p , 0  is introduced, entering as a dampening factor, modifying Eqution (1) as:
α s 2 ( p 2 ) p 4 α s 2 ( p 2 + p , 0 2 ) ( p 2 + p 0 , 2 ) 2 .
This cut-off is usually interpreted as a colour screening scale, though suggestions linking it to saturation in the spirit of Color Glass Condensate (CGC) [24] instead, have been put forth.
The extension of this model, Angantyr [25,26], underpins all heavy ion results presented in this review. Before it is briefly outlined in the following, it should be mentioned that it is not at all obvious why a model assuming that soft particle production can be correctly described with  2 2  parton scatterings as a starting point. Other models, such as the already mentioned CGC models, as well as the also already mentioned EPOS model, have other starting points. In the latter case, soft inelastic processes are represented directly as cut Pomerons.

Angantyr for Heavy Ion Collisions

Angantyr is initialized by a Glauber calculation [27,28] with particular emphasis on colour fluctuations in the initial state [29,30,31]. Through the Good–Walker formalism [32], these can be translated to semi-inclusive nucleon-nucleon cross sections, which can be used to fit the parameters of the initial state model from pp data only. Crucially, this treatment can distinguish not only which nucleons are hit and which are spectating, but also how the nucleons are hit, i.e., whether they participate elastically, diffractively or absorbed via a colour exchange.
Once it is determined which nucleons interact, multiple partons from each nucleon can interact, in the same spirit as the MPI model for pp. The ansatz for the Angantyr model, however, is the wounded nucleon model3 [34]. In this model, an ad hoc multiplicity function is applied to each wounded nucleon in the collision, which can then be stacked to construct a pA or an  A A  collision. The process is sketched in Figure 1 (left). Instead of a function, Angantyr uses a dynamical model, where a wounded nucleon is generated as normal MPIs, but with reduced phase space. Two wounded nucleons correspond to a pp collision, and parameters of the model can, therefore, be fitted to pp only. Interestingly, it is possible to interpret the addition of an additional wounded nucleon to a pp collision (making it a proton-Deuteron collision) in the cut Pomeron language [25]. This may point to the reasons why the seemingly very different models such as the CGC and EPOS models mentioned above, and the Angantyr model, can reproduce very similar physics from very different starting points. This is a direction of research that is currently left unexplored.
As shown in Figure 1 (right), this procedure reproduces centrality measures (here,  E  in the forward direction, which is the ATLAS choice of centrality measure) in pA well (the model performs similarly well for  A A  [25]). This is essential for providing direct comparisons to centrality-binned experimental data, where centrality is experimentally determined by cutting percentiles in this quantity.
Relating the forward production (centrality bins) to mid-rapidity production yields the comparisons shown for pPb in Figure 2 (left), and for PbPb in Figure 2 (right). This basic extension of the MPI model to heavy ion collisions performs very well for total multiplicities, as observed.
The primary purpose of the Angantyr model is to provide a “blank canvas” upon which microscopic, string-based models for collectivity can be constructed. This is the main topic of this mini-review, which will be expanded upon in the coming sections.

3. The Lund String Model and Microscopic Collectivity

The basis of the microscopic model for collectivity lies in the so-called Lund string. This is the physics model underlying string fragmentation in the Pythia Monte Carlo event generator. It represents a specific version of a more general class of “hadronic string” or “flux-tube” models, all of which make distinct phenomenological assumptions about the dynamics, particularly of string breakups. In this section, the main components4 of the Lund string model for a single string will be introduced in Section 3.1, proceeding to systems of multiple strings in Section 3.2. Colour reconnection is presented first in Section 3.3, with string shoving and rope hadronization presented in Section 3.4 and Section 3.5, respectively.

3.1. Physics of a Single Lund String

The Lund string model is built upon the observation from lattice QCD (and previously from onium spectra) that the confining potential between a  q q ¯  pair is linear, possessing a string tension of approximately  κ  1 GeV/fm. This is given a physical interpretation as flux tubes, akin to vortex lines in a superconductor. The Lund string model [5] is a 1 + 1 dimensional phenomenological model realization of this concept. It is based on the assumption that the fragmentation dynamics are insensitive to the width of the string, allowing it to be approximated by a massless relativistic string of infinitesimal width. Mesons are modelled as small  q q ¯  string pieces (“yo-yo” modes) produced by successively fragmenting pieces off the original string. The main components of the original string model include the breakup mechanism dictating the phase space distribution of produced hadrons, as well as the inclusion of gluons as kinks on the string.
The string breaking can be modelled as an iterative process, where each break produces a new hadron. With minimal constraints, amounting to the requirement that the same result should be obtained whether fragmenting from one end of the string or the other, the distribution of the momentum fraction (z) of the remaining light-cone momentum taken away by the hadron produced in each splitting is given as:
f ( z ) ( 1 z ) a z exp b m 2 z ,
where  m  is the hadron’s transverse mass ( m 2 = m 2 + p 2 ), and a and b are parameters determined from experiments. This formulation fixes the longitudinal degrees of freedom in the string-breaking process. The transverse momentum and flavour/mass of the quarks are determined from the probability of the  q q ¯  pair tunnelling through a region of size  m / κ , given by:
d P d 2 p κ exp ( π m , q 2 / κ ) = exp ( π m q 2 / κ ) × exp ( π p , q 2 / κ ) .
Due to the factorization in the above expression, the transverse momenta and mass of the quarks can be treated separately. In practice, this is carried out by generating Gaussian  p -kicks, with a parameter  σ p  governing the width of the distribution5. The suppression of strange quarks relative to u or d quarks is described by a factor  ρ = exp ( π ( m s 2 m u 2 ) / κ ) , which is also treated as a parameter6.

3.2. Interactions of Multiple Strings

In collisions of protons (and thus also in larger collision systems), strings will overlap with each other (see Figure 3, left). It is therefore reasonable to question whether the initial assumption of strings being infinitely thin should be relaxed. Discarding the assumption altogether, allowing strings to have a transverse extension, naturally leads to string interactions. Results from lattice QCD [39] can estimate the field shape, as shown in the black points in Figure 3. As indicated by the fitted curve, the field is well approximated by a Gaussian7. However, the transverse size of the string’s electrical field as shown in the figure should be taken with caution. Since lattice calculations are performed in arbitrary “lattice units”, a connection to a physical quantity is necessary to translate to physical units. Typically, the string tension is used for such calculations, but further assumptions are needed about the fraction of energy going into the field versus the fraction needed to break the condensate. Therefore, these types of calculations cannot provide a precise estimate of the string width, but they do give a good idea of the order—which is significant enough for string interactions to occur.
In the Lund string model, the string is produced by chaining together quarks and gluons, which were previously produced by the hard scattering and successive emissions from the parton shower. Due to causality, the field is initially very thin, and interactions between strings do not depend on overlapping fields. Instead, they are formulated in terms of “colour reconnections”, the term given to sub-leading colour corrections to the string topology provided by the parton shower. It is only after this phase that the string reaches a size where the (colour)-electric fields start interacting as fields, pushing each other apart. Eventually, the strings will hadronize, and if they are still overlapping, they may hadronize with a larger string tension than  κ = 1  GeV/fm. After hadronization, hadrons may interact with each other in a hadronic cascade.
Importantly, the time allowed for the strings to interact with each other is not a free parameter of the model. The derivation of the fragmentation function Eqution (3) also yields, as a by-product, the distribution of string breakup vertices in proper time ( τ ). The distribution is more conveniently written in terms of  Γ = ( κ τ ) 2 , given by:
P ( Γ ) d Γ = Γ a exp ( b Γ ) d Γ .
From this distribution, the average string break time is:
τ 2 = 1 + a b κ 2 ,
which, with reasonable values for a and b (determined from total multiplicities in  e + e  using Eqution (3)), gives  τ 2 2  fm, albeit with a tail to larger times. While this cannot be directly interpreted as the hadron formation time (it is unclear from the model whether a hadron can be considered produced once the string breaks or after, for example, one full oscillation of the yo-yo mode), this figure establishes some conceptual facts about the microscopic model. First and foremost, the time allowed for strings to interact before hadronization is significantly shorter (a factor of 5–10) than in hydrodynamic models assuming a QGP. Secondly, there is no single time defined for “string freezeout”. The process will inevitably result in some strings hadronizing earlier and some later8. Third, the earlier hadronization compared to QGP-based models allows the subsequent hadronic cascade to have a much larger impact on observable results. This latter effect is not discussed in this review, but interested readers are referred to refs. [42,43].
Ordered roughly by time after the initial collision, the microscopic model for collective behaviour consists of the following components:
  • τ  0–0.6 fm: Parton shower and colour reconnection. Emissions down to the parton shower cut-off scale and altered colour topologies, which can give rise to short-range collective effects [44,45] and enhanced baryon production [46,47].
  • τ  0.6–1.4 fm: String shoving. Strings reach a transverse size where they interact maximally before hadronizing, giving rise to soft collective effects [48], and affecting jets as well [49].
  • τ  1.4 fm: Rope hadronization. Strings can hadronize in a colour multiplet if they are still sufficiently close together, forming a so-called colour rope [50], which leads to enhanced strangeness and baryon production [41].
  • After  τ  1.4 fm: Hadronic cascade. Hadronic rescattering [51] results in significant corrections to both flavour and flow, particularly in the dense environment of  A A  collisions [43].

3.3. Colour Reconnections

The colour topology delivered by the matrix element calculation and parton shower uses the “leading colour” limit9, in which the string topology is uniquely determined by the fact that the probability of two random colours are the same, vanishes. As an example, two  q q ¯  pairs in the leading colour limit will never be allowed to form a singlet, but will always form an octet. Thus, in  3 3 ¯ = 8 1 , the  1  vanishes in the limit.
Colour reconnection models are an attempt to reintroduce the colour topologies which should have been generated, had the perturbative calculation been carried out in full colour. The simplest reconnection models aims at reintroducing the singlet. This is, for example, the case for the default Pythia reconnection model [23] (often referred to as “the MPI based model”). The main motivation for the model was the observation that with completely independent MPIs,  p  will not rise with multiplicity, as was observed at UA1 [53]. If MPIs are allowed to reconnect, and in essence, transfer some multiplicity to  p , this is achieved.
Since one cannot perform a full-colour correction from first principles, colour reconnection models contain some ad hoc modelling elements. In the case of string-based models, this is often centred around the notion of reducing total string length. This is quantified in the so-called  λ -measure, which measures how many hadrons of some reference mass  m 0 , a given string piece has room (in phase space) to produce. For a simple  q q ¯  system it can easily be written as, e.g.,  λ q q ¯ = ln ( m q q ¯ 2 / m 0 2 ) , but for more complicated topologies, approximate expressions must be used, e.g.,
λ = i = 0 n ln 1 + ( k i p i + k i + 1 p i + 1 ) 2 m 0 2 ,
where k is 1 for quarks and  1 / 2  for gluons, as the gluon momentum is shared between the two connecting string pieces.
It has been shown that even the simplest colour reconnection model in Pythia can give rise to some QGP-like effects, in particular baryon–meson ratios and flow [44]. The effects are, however, only short range in rapidity [45] and can therefore not explain observations such as the pp ridge [13] or  v 2  with rapidity gap between particles of interest.
More advanced colour reconnection models have been developed, in particular to further address baryon–meson ratios. Most noteworthy is the so-called sub-leading colour or “QCD-CR” model [46]. The model aims to address also the higher neglected multiplets when constructing string topologies, rather than just the singlet, for example, the anti-triplet in  3 3 = 6 3 ¯ , or the triplet and anti-sextet in  3 8 = 15 6 ¯ 3 . Importantly, the  3 ¯  and  6 ¯  in these states involves colour states10 which, in the string model, are interpreted as “junctions”: states where three string legs are connected to a single point, and which carry intrinsic baryon number [54]. Allowing the formation of such states in colour reconnection (in pairs of junction–anti-junction, so it does not violate baryon number conservation) increases the amount of produced baryons when a high number of MPIs increases the chance of a junction colour reconnection. Technically, one performs an approximate calculation of the weight factors based on the  SU ( 3 )  algebra, determining which configurations are more or less likely. The actual configuration is then chosen based on which of the possible configurations minimizes the  λ -measure the most.
The model was introduced to address discrepancies of baryon yields in the first years of LHC running wrt.  e + e  (see, e.g., ref. [55]), but crucially also allows for the production of double heavy flavour baryons, not allowed in the normal string model. In particular, production of heavy flavour baryons has been a large success for the model, where the  Λ c / D  ratio vs.  p  in pp collisions (see Figure 4 (left)) was found by ALICE [56] to differ drastically from expectation. The colour reconnection model gives a satisfactory description. In ion collisions, later refinements of the model [57], taking into account the non-negligible mass of the charm quark, allowed predictions also for pA, as shown in Figure 4 (right). Also, here, the description is drastically improved.
The last word is, however, not said for this type of model. Further studies of charm baryon yields have revealed that, for example, charmed  Ξ  is not equally well reproduced [57,59]. Some are on account of further enhancement of strangeness (see Section 3.5), but not all can be explained in this way. Colour reconnection remains to be an active area of study.

3.4. The String-Shoving Mechanism

As explained at the beginning of Section 3.2, the transverse (colour)-electric field can be approximated well as a Gaussian. After the string has had time after its initial creation to expand to its full transverse size, strings will start “shoving” each other, with a force calculable from the field. If the field is  E = N exp ( ρ 2 / 2 R 2 ) , where  ρ  is the radius in cylindrical coordinates and R is the equilibrium radius. N is a normalization factor, which is determined by letting the energy in the field correspond to a fraction g of the total string tension.
It is now assumed that all multi-string colour configurations such as  3 3  have already been correctly colour reconnected. This means that only octet configurations are left, while singlets have vanished. The remaining string dipoles will therefore only repel each other and never attract. The energy per unit length of two strings overlapping is  d 2 ρ ( E 1 + E 2 ) 2 / 2 , from which the force between two strings transversely separated by  d  can be calculated as:
f ( d ) = g κ d R 2 exp d 2 4 R 2 .
This is the main ingredient of the string-shoving mechanism. Technical complications arise when implementing this in a Monte Carlo event generator. String pieces are not aligned in parallel, but a suitable Lorentz transformation can be found to align each pair in parallel planes [40]. Also, simply conducting a step-wise time evolution is not practical from a computational point of view. Instead, the shoving model is implemented in a similar spirit as a parton shower, ordered in the  p  pushes.

3.4.1. Soft Collective Effects

The main intent of the shoving model is to be able to generate anisotropic flow as a response to the spatial initial conditions. In proton collisions the initial conditions are currently just produced as the convolution of two Gaussian mass distributions, which is itself a Gaussian, i.e., fully symmetric. Higher order flow is therefore produced only from fluctuations in the initial conditions, and not from intrinsic asymmetry11. Overall observables, such as the long-range ridge first measured in pp by CMS [13], are well reproduced as shown in Figure 5. (Note that the parameter g in the figure is not identical to g in Eqution (8), but differs by a normalizing factor).
Since the model is implemented as part of the normal Pythia package, it can also be used to study the subtle effects of introducing other scales to the collision. ATLAS has measured [61] the ridge in collisions with a Z boson present, instead of the normal QCD minimum bias collision systems. The presence of the Z introduces the  m Z  scale as the largest hard scale of the collision, which could potentially alter the distribution (in  p ) of MPIs wrt. minimum bias. As shown in Figure 5 (right), the shoving model does a good job at reproducing the ridge also under these special circumstances.
Small collision systems are, however, not the main avenue for flow measurements. The shoving model has been applied, through Angantyr, to  A A  collisions as well12. Here, flow is quantified in flow coefficients ( v n s); see, e.g., ref. [14] for an overview. Currently, the model does not perform well in full collisions due to the presence of many soft gluons—leading to many strings with many kinks—which are difficult to handle technically [40]. Instead, toy approaches have been considered with some success.
Since the multiplicity generated by a single string is well known (roughly one hadron per unit of rapidity), a system consisting of straight strings, i.e., without soft gluons, corresponding to the multiplicity of  A A  collisions in a given centrality interval can be set up. In ref. [40], this was carried out along with spatial initial conditions resembling those of the given centrality interval. Furthermore, since the initial conditions are known in this case, the flow does not have to be calculated in the same way as the experiment13, but can be calculated wrt. the true reaction plane, as well as correlations with the initial geometry. In Figure 6 (left), results for  v 2  in this toy configuration, compared to data from ALICE [63], is shown. While the agreement is by no means perfect, it should be emphasized that this calculation is performed using the same set of parameters as used for smaller collision systems in the same paper. This suggests that flow can indeed have the same origin across collision systems, and that the origin may be string interactions.
The normal explanation for flow in heavy ion collisions, is a hydrodynamic origin, due to the deconfined QGP phase. In hydrodynamic models [64] as well as multiphase transport [65], it is known that the correlation between initial state eccentricity [66] and final state  v 2  is linear. In Figure 6 (right), the same response is shown for the shoving model, revealing that it exhibits similar behaviour in dense systems. This suggests, as is perhaps not surprising, that hydrodynamic behaviour is not limited to deconfined systems. From this, it follows that flow signals cannot be taken as proof of a QGP if the same signal can be generated from the shoving model without QGP.

3.4.2. Effects on Jets

As indicated above, the shoving model applies everywhere there are two or more strings available, irrespective of the collision system or hard process. Since jet quenching in  A A  collisions constitutes some of the most solid evidence for QGP production, the so far unsuccessful search for jet modification in pp (see, e.g., [22]) is understandably a topic of some interest. When discussing jet modifications from string shoving, it is important to first emphasize that the model makes no distinction between jets, and the rest of the collision event. In the string model, every coloured parton must be connected to the rest of the event through a string, and when there are strings, there will be shoving effects. However, the small size of pp collisions combined with the fast-moving jet parton makes jet modifications from string shoving very difficult to observe. However, string shoving is a final state interaction effect, and as such, jet modifications will show up given enough statistics/patience.
Consider a pp collision where a Z is produced together with a reasonably high  p  jet. With a reasonable lower bound on jet  p , here14 p , j e t > 80  GeV, and the jet cross section  σ j = p , 0 d p , j d σ d p , j  is calculable in perturbative QCD, but receives large corrections from MPI and hadronization. In Figure 7, the jet cross section in  s = 7  TeV collisions is shown as function of jet R (= Δ η 2 + Δ ϕ 2 ). As expected, the effect is not nearly as large as the effect from MPIs, which at large R blows up as the jet cone radius covers more and more of the underlying event. At large R, shoving contributes to the jet cross section at the same level as hadronization, i.e., at the level of a few percent. The main point of the figure is twofold. Firstly, one must be precise when asking for studies of jet modifications due to final state effects because, as shown, the jet definition can be doctored in such a way that it reaches way out into the underlying event, where even basic quantities like the jet cross section are affected. Second, measurements of jet modifications due to such miniscule final state effects are indeed—as expected—difficult to measure.
An example of a concrete observable, the jet mass of large  R = 0.7  jets is shown in Figure 7 (right) compared with measurements by CMS at  s = 7  TeV. Sizeable effects (over 10%) show up at low jet masses of less than 20 GeV, but current measurements have no chance of distinguishing.

3.5. Rope Hadronization

After the shoving phase, the string will hadronize. Hadrons will have properties as determined by Equations (3) and (4). If the strings are, at this point, still overlapping, they may, however, form a “rope” [41]. Connecting back to the explanations on colour multiplets given in Section 3.3, this is the effect that the string tension will be enhanced for higher multiplets15.
SU ( 3 )  multiplets can conveniently be characterized by the two quantum numbers  ( p , q ) , which can be given the heuristic interpretation as the number of parallel and anti-parallel strings acting together coherently when overlapping. The simple example of  3 3 = 6 3 ¯  can as such be written as  ( 1 , 0 ) ( 1 , 0 ) = ( 2 , 0 ) ( 0 , 1 ) . The quadratic Casimir ( C 2 ) of a given  ( p , q )  multiplet can be written in terms of p and q, normalized by  C 2  of a triplet. Since the string tension scales like  C 2  in lattice calculations [81], this can be used to calculate the effective string tension ( κ ˜ ) of any multiplet as:
κ ˜ ( p , q ) κ ( 1 , 0 ) = C 2 ( p , q ) C 2 ( 1 , 0 ) = 1 4 ( p 2 + p q + q 3 + 3 p + 3 q ) .
In the Pythia rope model, the rope hadronizes by breaking up one “thread” at the time. One such breaking will therefore take the multiplet from the state  ( p , q ) ( p 1 , q ) . The effective string tension in the string break therefore becomes:
κ ˜ κ = 2 p + q + 2 4 .
The energy density in a rope can quickly grow large enough to make a dramatic impact on final state observables (see Section 3.5.1). It will, however, not exceed more than around a factor 2 above the vacuum string tension, except in very central  A A  collisions, where a factor 2 more can be reached. In QGP based models it is often argued (see, e.g., ref. [82]) that initial state energy densities ( d E / d 3 x ), which, for example, in CGC-based models can reach over 100 GeV/fm 3  in PbPb collisions at LHC [40,83,84], are so large that plasma creation is inevitable. Crucially, the energy density in a rope scenario never reaches such levels, as shown in Figure 8. Other models of string interactions, crucially the string percolation model [85], have the opposite conclusion: that the percolation of strings is potentially the mechanism responsible for QGP formation [86].

3.5.1. Effects on Inclusive Production

When the string tension rises, the main effect is that suppression of strangeness in Eqution (4) decreases. The derived quantity  ρ  is replaced by an effective quantity  ρ ˜ , which can be obtained directly as:
ρ ˜ = exp π ( m s 2 m u 2 ) κ ˜ = ρ κ / κ ˜ ,
where the exponent can be obtained directly from  ( p , q )  through Eqution (10). Similar expressions exist for the other parameters of the hadronization model; see refs. [41,87].
Crucially, the base value of  ρ  is not touched, but fixed to  e + e  data once and for all. Furthermore, the model does not require any fitting to heavy ion data at all. Once the geometry of the  A A  collision is set up, everything else is fixed by smaller collision systems.
Production of multi-strange baryons is now affected both by the introduction of junction in the colour reconnection mechanism explained above, as well as the enhancement of strange quark production. In Figure 9, a summary of strange multi-strange hadron production across pp, pPb and PbPb is shown. While the model does not fully capture the data, it exhibits the correct trend across all collision systems. As was the case for string shoving, work on extending this facet of the model to all collision systems is still an ongoing endeavour.

4. Conclusions and Outlook

The discovery of collective effects in small collision systems has spurred renewed interest in hadronization models, a field that has otherwise been quite dormant since the late 1980s and early 1990s. This mini-review discusses ongoing attempts to construct a QGP-free model to describe behaviour in both small and large collision systems, traditionally ascribed to QGP formation. However, developments are also occurring in other areas, both connected and unrelated to these endeavours. For the string hadronization model, this includes attempts to modify the usual Eqution (4) to a thermal one [88], colour reconnection corrections to the top quark mass [89], calculations of spatial vertices in the fragmentation [90] underlying hadronic rescattering, hyperfine splitting effects [91], and more theoretical attempts to investigate the single string entanglement entropy [92], with some phenomenological considerations also made [93]. Similarly, for the Herwig cluster model, efforts such as cluster reconnections [94,95] have made progress in the production of multi-strange baryons, akin to string interactions.
A very recently emerging direction involves the formulation of hadronization models using machine learning-based empirical models, with a special focus on deep learning techniques. Complementary approaches are being developed for both the string and the cluster models. The HadML initiative has demonstrated a proof of principle that the cluster model can be parametrized using GANs [96], and how such models can be tuned to real data [97], providing a path forward for real applications. The MLHad initiative [98] has similarly shown how the string model can be parametrized using normalizing flows [99], as well as how string hadronization can be re-parametrized [100] to allow for more efficient parameter variation and ultimately more robust comparisons to data. While it is still unclear whether such efforts will contribute directly to a deeper physics understanding of hadronization, it is evident that they will aid in making better quantitative comparisons to data, and ultimately in ruling out models that cannot globally describe data.
Ahead of us lies the high luminosity phase of the LHC. This phase is expected to bring new insights into hadronization, not least through improved statistics on heavy flavour production, particularly in the baryon sector. This advancement is expected to impact the colour reconnection models as new knowledge emerges. The upcoming electron–ion collider, as well as further reanalysis of LEP data until the FCC-ee, will hopefully shed further light on collectivity in the smallest systems.

Funding

Support from Vetenskapsrådet, contract number 2023-04316 and the Knut and Alice Wallenberg foundation, contract number 2017.003 is gratefully acknowledged.

Acknowledgments

I wish to thank Manuel Szewc for a well placed comment on Deep Learning techniques.

Conflicts of Interest

The author declare no conflicts of interest.

Notes

1
QGP formation predicts several other types of modifications with respect to proton collisions, such as jet quenching, heavy quark modifications, altered resonance production, and several more. The focus in this mini-review will, however, be on strangeness and flow.
2
This review is thus not a general review on collectivity in small systems. The reader is referred to excellent reviews in refs. [14,18,19,20,21,22] for other perspectives, both from an experimental and theoretical starting point. In particular ref. [20] offers a RHIC perspective, which is not covered in this review.
3
In this sense, Angantyr is more of a successor to the FRITIOF model [33], which had the same starting point.
4
For a more comprehensive, recent review, see ref. [9].
5
In principle,  σ p = κ / π ( 0.25  GeV ) 2 , but fits to LEP data suggest this number to be higher, indicating that a non-negligible factor of the  p  comes from another source.
6
In principle, this number could be obtained if a suitable value for quark masses were inserted. However, current quark masses lead to too little strangeness suppression, and constituent quark masses to too much.
7
The other curve, labelled “Clem profile”, is derived from Landau–Ginzburg theory and is given by  E = K 0 ( ρ 2 + ξ v 2 / λ )  (where  K 0  is a modified Bessel function,  ρ  the cylindrical radius coordinate,  ξ v  the condensate coherence length, and  λ  the penetration depth.) See ref. [40] for more details.
8
This, in principle, ought to give rise to a mixed phase of strings and hadrons. However, this possibility has not yet been pursued in this framework.
9
Formally obtained [52] by letting  N c  ( N c  is the number of colours), while keeping  α s N c  fixed. This eliminates coherence effects which are suppressed by a factor  1 / N c 2  and allows for a simple representation of gluons as colour–anti-colour.
10
Such states have no analogy in  N c  limit, which corresponds only to dipole-like connections. These states represent  ϵ i j k  structures in colour space, which are explicitly  N c = 3 .
11
It is possible to calculate initial geometries in other types of approaches like BFKL evolution in impact parameter space [60]. While the first steps have been made to interface this, it is not yet mature enough for real calculations.
12
As a historical note, it should be mentioned that the string-shoving model is not the first attempt to generate anisotropic flow from string interactions. Abramovsky et al. introduced a simple but similar idea already in 1988 [62], before the idea of flow had been coupled to QGP formation at all.
13
In the experiment, the reaction plane angle is not known event-by-event, and one must construct flow coefficients from two -or multiparticle correlations. In calculations, in particular in  A A  collisions where the non-flow contributions are small, the true reaction plane can be used.
14
As well as requiring a Z with  p , Z > 40  GeV, the jet is the leading anti- k  [67] jet with  Δ ϕ Z , j > 3 π / 4 .
15
The rope hadronization idea dates back to the 1980s [50], and has since been pursued in similar forms by many authors. Some very similar to the Pythia ropes, others in the same general class of string fusion models. Some implemented in event generators [69,70,71,72,73,74,75], and others with purely theoretical work [76,77,78,79,80].

References

  1. Rafelski, J.; Muller, B. Strangeness Production in the Quark—Gluon Plasma. Phys. Rev. Lett. 1982, 48, 1066, Erratum in Phys. Rev. Lett. 1986, 56, 2334. [Google Scholar] [CrossRef]
  2. Stock, R.; Eschke, J.; Heck, W.; Kabana, S.; Kuehmichel, A.; Lahanas, M.; Lee, Y.; Margetis, S.; Renfordt, R.; Roehrich, D.; et al. Strangeness enhancement in central S + S collisions at 200-GeV/nucleon. Nucl. Phys. A 1991, 525, 221C–226C. [Google Scholar] [CrossRef]
  3. Ollitrault, J.Y. Anisotropy as a signature of transverse collective flow. Phys. Rev. D 1992, 46, 229–245. [Google Scholar] [CrossRef]
  4. Ackermann, K.H. Elliptic flow in Au + Au collisions at (S(NN))**(1/2) = 130 GeV. Phys. Rev. Lett. 2001, 86, 402–407. [Google Scholar] [CrossRef]
  5. Andersson, B.; Gustafson, G.; Ingelman, G.; Sjöstrand, T. Parton Fragmentation and String Dynamics. Phys. Rep. 1983, 97, 31–145. [Google Scholar] [CrossRef]
  6. Webber, B.R. A QCD Model for Jet Fragmentation Including Soft Gluon Interference. Nucl. Phys. B 1984, 238, 492–528. [Google Scholar] [CrossRef]
  7. Buckley, A.; Butterworth, J.; Gieseke, S.; Grellsheid, D.; Höche, S.; Hoeth, H.; Krauss, F.; Lönnblad, L.; Nurse, E.; Richardson, P.; et al. General-purpose event generators for LHC physics. Phys. Rep. 2011, 504, 145–233. [Google Scholar] [CrossRef]
  8. Campbell, J.M.; Diefenthaler, M.; Hobbs, T.J.; Höche, S.; Isaacson, J.; Kling, F.; Mrenna, S.; Reuter, J.; Alioli, S.; Andersen, J.R.; et al. Event Generators for High-Energy Physics Experiments. In Proceedings of the Snowmass 2021, Seattle, WA, USA, 17–26 July 2022; p. 3. [Google Scholar]
  9. 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. SciPost Phys. Codeb. 2022, 8. [Google Scholar] [CrossRef]
  10. Bellm, J.; Gieseke, S.; Grellscheid, D.; Plätzer, S.; Rauch, M.; Reuschle, C.; Richardson, P.; Schichtel, P.; Seymour, M.H.; Siódmok, A.; et al. Herwig 7.0/Herwig++ 3.0 release note. Eur. Phys. J. C 2016, 76, 196. [Google Scholar] [CrossRef]
  11. Bothmann, E.; Chahal, G.S.; Höche, S.; Krause, J.; Krauss, F.; Kuttimalai, S.; Liebschner, S.; Napoletano, D.; Schönherr, M.; Schulz, H.; et al. Event Generation with Sherpa 2.2. SciPost Phys. 2019, 7, 034. [Google Scholar] [CrossRef]
  12. ALICE Collaboration. Enhanced production of multi-strange hadrons in high-multiplicity proton–proton collisions. Nat. Phys. 2017, 13, 535–539. [Google Scholar] [CrossRef]
  13. Khachatryan, V.; Sirunyan, A.M.; Tumasyan, A.; Adam, W.; Bergauer, T.; Dragicevic, M.; Erö, J.; Fabjan, C.; Friedl, M.; Frühwirth, R.; et al. Observation of Long-Range Near-Side Angular Correlations in Proton-Proton Collisions at the LHC. J. High Energy Phys. 2010, 2010, 91. [Google Scholar] [CrossRef]
  14. The ALICE experiment—A journey through QCD. arXiv 2022, arXiv:2211.04384.
  15. Chen, Y.C.; Lee, Y.J.; Chen, Y.; Chang, P.; McGinn, C.; Sheng, T.A.; Innocenti, G.M.; Maggi, M. Analysis note: Two-particle correlation in e+e collisions at 91-209 GeV with archived ALEPH data. arXiv 2023, arXiv:2309.09874. [Google Scholar]
  16. Werner, K. Core-corona separation in ultra-relativistic heavy ion collisions. Phys. Rev. Lett. 2007, 98, 152301. [Google Scholar] [CrossRef] [PubMed]
  17. Kanakubo, Y.; Tachibana, Y.; Hirano, T. Unified description of hadron yield ratios from dynamical core–corona initialization. Phys. Rev. C 2020, 101, 024912. [Google Scholar] [CrossRef]
  18. Li, W. Observation of a ‘Ridge’ correlation structure in high multiplicity proton–proton collisions: A brief review. Mod. Phys. Lett. A 2012, 27, 1230018. [Google Scholar] [CrossRef]
  19. Schlichting, S.; Tribedy, P. Collectivity in Small Collision Systems: An Initial-State Perspective. Adv. High Energy Phys. 2016, 2016, 8460349. [Google Scholar] [CrossRef]
  20. Nagle, J.L.; Zajc, W.A. Small System Collectivity in Relativistic Hadronic and Nuclear Collisions. Ann. Rev. Nucl. Part. Sci. 2018, 68, 211–235. [Google Scholar] [CrossRef]
  21. Citron, Z.; Dainese, A.; Grosse-Oetringhaus, J.F.; Jowett, J.M.; Lee, Y.-J.; Wiedemann, U.A.; Winn, M.; Andronic, A.; Bellini, F.; Bruna, E.; et al. Report from Working Group 5: Future physics opportunities for high-density QCD at the LHC with heavy-ion and proton beams. CERN Yellow Rep. Monogr. 2019, 7, 1159–1410. [Google Scholar] [CrossRef]
  22. Adolfsson, J.; Andronic, A.; Bierlich, C.; Bozek, P.; Chakraborty, S.; Christiansen, P.; Chinellato, D.D.; Fries, R.J.; Gustafson, G.; van Hees, H.; et al. QCD challenges from pp to A—A collisions. Eur. Phys. J. A 2020, 56, 288. [Google Scholar] [CrossRef]
  23. Sjöstrand, T.; van Zijl, M. A Multiple Interaction Model for the Event Structure in Hadron Collisions. Phys. Rev. D 1987, 36, 2019. [Google Scholar] [CrossRef] [PubMed]
  24. Gelis, F.; Iancu, E.; Jalilian-Marian, J.; Venugopalan, R. The Color Glass Condensate. Ann. Rev. Nucl. Part. Sci. 2010, 60, 463–489. [Google Scholar] [CrossRef]
  25. Bierlich, C.; Gustafson, G.; Lönnblad, L.; Shah, H. The Angantyr model for Heavy-Ion Collisions in PYTHIA8. J. High Energy Phys. 2018, 2018, 134. [Google Scholar] [CrossRef]
  26. Bierlich, C.; Gustafson, G.; Lönnblad, L. Diffractive and non-diffractive wounded nucleons and final states in pA collisions. J. High Energy Phys. 2016, 2016, 139. [Google Scholar] [CrossRef]
  27. Miller, M.L.; Reygers, K.; Sanders, S.J.; Steinberg, P. Glauber modeling in high energy nuclear collisions. Ann. Rev. Nucl. Part. Sci. 2007, 57, 205–243. [Google Scholar] [CrossRef]
  28. Rybczynski, M.; Stefanek, G.; Broniowski, W.; Bozek, P. GLISSANDO 2: GLauber Initial-State Simulation AND mOre…, ver. 2. Comput. Phys. Commun. 2014, 185, 1759–1772. [Google Scholar] [CrossRef]
  29. Gribov, V.N. Glauber corrections and the interaction between high-energy hadrons and nuclei. Sov. Phys. JETP 1969, 29, 483–487. [Google Scholar]
  30. Heiselberg, H.; Baym, G.; Blaettel, B.; Frankfurt, L.L.; Strikman, M. Color transparency, color opacity, and fluctuations in nuclear collisions. Phys. Rev. Lett. 1991, 67, 2946–2949. [Google Scholar] [CrossRef]
  31. Blaettel, B.; Baym, G.; Frankfurt, L.L.; Heiselberg, H.; Strikman, M. Hadronic cross-section fluctuations. Phys. Rev. D 1993, 47, 2761–2772. [Google Scholar] [CrossRef]
  32. Good, M.L.; Walker, W.D. Diffraction disssociation of beam particles. Phys. Rev. 1960, 120, 1857–1860. [Google Scholar] [CrossRef]
  33. Andersson, B.; Gustafson, G.; Nilsson-Almqvist, B. A Model for Low p(t) Hadronic Reactions, with Generalizations to Hadron - Nucleus and Nucleus-Nucleus Collisions. Nucl. Phys. B 1987, 281, 289–309. [Google Scholar] [CrossRef]
  34. Bialas, A.; Bleszynski, M.; Czyz, W. Multiplicity Distributions in Nucleus-Nucleus Collisions at High-Energies. Nucl. Phys. B 1976, 111, 461–476. [Google Scholar] [CrossRef]
  35. Bierlich, C.; Buckley, A.; Butterworth, J.; Christensen, C.H.; Corpe, L.; Grellscheid, D.; Grosse-Oetringhaus, J.F.; Gutschow, C.; Karczmarczyk, P.; Klein, J.; et al. Robust Independent Validation of Experiment and Theory: Rivet version 3. SciPost Phys. 2020, 8, 026. [Google Scholar] [CrossRef]
  36. Bierlich, C.; Buckley, A.; Christensen, C.H.; Christiansen, P.H.L.; Duncan, C.B.; Grosse-Oetringhaus, J.F.; Karczmarczyk, P.; Kirchgaeßer, P.; Klein, J.; Lönnblad, L.; et al. Confronting experimental data with heavy-ion models: RIVET for heavy ions. Eur. Phys. J. C 2020, 80, 485. [Google Scholar] [CrossRef]
  37. ATLAS Collaboration. Measurement of the centrality dependence of the charged-particle pseudorapidity distribution in proton–lead collisions at  s NN  = 5.02 TeV with the ATLAS detector. Eur. Phys. J. C 2016, 76, 199. [Google Scholar] [CrossRef] [PubMed]
  38. ALICE Collaboration. Centrality dependence of the pseudorapidity density distribution for charged particles in Pb-Pb collisions at  s NN  = 5.02 TeV. Phys. Lett. B 2017, 772, 567–577. [Google Scholar] [CrossRef]
  39. Cea, P.; Cosmai, L.; Cuteri, F.; Papa, A. Flux tubes in the SU(3) vacuum: London penetration depth and coherence length. Phys. Rev. D 2014, 89, 094505. [Google Scholar] [CrossRef]
  40. Bierlich, C.; Chakraborty, S.; Gustafson, G.; Lönnblad, L. Setting the string shoving picture in a new frame. J. High Energy Phys. 2021, 2021, 270. [Google Scholar] [CrossRef]
  41. Bierlich, C.; Gustafson, G.; Lönnblad, L.; Tarasov, A. Effects of Overlapping Strings in pp Collisions. J. High Energy Phys. 2015, 2015, 148. [Google Scholar] [CrossRef]
  42. da Silva, A.V.; Serenone, W.M.; Chinellato, D.D.; Takahashi, J.; Bierlich, C. Studies of heavy-ion collisions using PYTHIA Angantyr and UrQMD. arXiv 2020, arXiv:2002.10236. [Google Scholar]
  43. Bierlich, C.; Sjöstrand, T.; Utheim, M. Hadronic rescattering in pA and AA collisions. Eur. Phys. J. A 2021, 57, 227. [Google Scholar] [CrossRef]
  44. Ortiz Velasquez, A.; Christiansen, P.; Cuautle Flores, E.; Maldonado Cervantes, I.; Paić, G. Color Reconnection and Flowlike Patterns in pp Collisions. Phys. Rev. Lett. 2013, 111, 042001. [Google Scholar] [CrossRef] [PubMed]
  45. Bierlich, C. Microscopic collectivity: The ridge and strangeness enhancement from string–string interactions. Nucl. Phys. A 2019, 982, 499–502. [Google Scholar] [CrossRef]
  46. Christiansen, J.R.; Skands, P.Z. String Formation Beyond Leading Colour. J. High Energy Phys. 2015, 2015, 3. [Google Scholar] [CrossRef]
  47. Bierlich, C.; Christiansen, J.R. Effects of color reconnection on hadron flavor observables. Phys. Rev. D 2015, 92, 094010. [Google Scholar] [CrossRef]
  48. Bierlich, C.; Gustafson, G.; Lönnblad, L. Collectivity without plasma in hadronic collisions. Phys. Lett. B 2018, 779, 58–63. [Google Scholar] [CrossRef]
  49. Bierlich, C. Soft modifications to jet fragmentation in high energy proton–proton collisions. Phys. Lett. B 2019, 795, 194–199. [Google Scholar] [CrossRef]
  50. 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]
  51. Sjöstrand, T.; Utheim, M. A Framework for Hadronic Rescattering in pp Collisions. Eur. Phys. J. C 2020, 80, 907. [Google Scholar] [CrossRef]
  52. ’t Hooft, G. A Planar Diagram Theory for Strong Interactions. Nucl. Phys. B 1974, 72, 461. [Google Scholar] [CrossRef]
  53. UA1 Collaboration. A Study of the General Characteristics of p p ¯  Collisions at  s  = 0.2-TeV to 0.9-TeV. Nucl. Phys. B 1990, 335, 261–287. [Google Scholar] [CrossRef]
  54. Sjöstrand, T.; Skands, P.Z. Baryon number violation and string topologies. Nucl. Phys. B 2003, 659, 243. [Google Scholar] [CrossRef]
  55. CMS Collaboration. Strange Particle Production in pp Collisions at  s  = 0.9 and 7 TeV. J. High Energy Phys. 2011, 2011, 64. [Google Scholar] [CrossRef]
  56. ALICE Collaboration. Λ c +  Production and Baryon-to-Meson Ratios in pp and p-Pb Collisions at  s NN  = 5.02 TeV at the LHC. Phys. Rev. Lett. 2021, 127, 202301. [Google Scholar] [CrossRef] [PubMed]
  57. Bierlich, C.; Gustafson, G.; Lönnblad, L.; Shah, H. The dynamic hadronization of charm quarks in heavy-ion collisions. arXiv 2023, arXiv:2309.12452. [Google Scholar]
  58. ALICE Collaboration. Measurement of Prompt  D 0 Λ c + , and  Σ c 0 , + + (2455) Production in Proton–Proton Collisions at  s  = 13 TeV. Phys. Rev. Lett. 2022, 128, 012001. [Google Scholar] [CrossRef]
  59. The ALICE Collaboration. Measurement of the production cross section of prompt  Ξ c 0  baryons at midrapidity in pp collisions at  s  = 5.02 TeV. J. High Energy Phys. 2021, 2021, 159. [Google Scholar] [CrossRef]
  60. Bierlich, C.; Rasmussen, C.O. Dipole evolution: Perspectives for collectivity and γ*A collisions. J. High Energy Phys. 2019, 2019, 26. [Google Scholar] [CrossRef]
  61. Aaboud, M.; Aad, G.; Abbott, B.; Abbott, D.C.; Abdinov, O.; Abeloos, B.; Abhayasinghe, D.K.; Abidi, S.H.; AbouZeid, O.S.; Abraham, N.L.; et al. Measurement of long-range two-particle azimuthal correlations in Z-boson tagged pp collisions at  s  = 8 and 13 TeV. Eur. Phys. J. C 2020, 80, 64. [Google Scholar] [CrossRef]
  62. Abramovsky, V.A.; Gedalin, E.V.; Gurvich, E.G.; Kancheli, O.V. Long Range Azimuthal Correlations in Multiple Production Processes at High-energies. JETP Lett. 1988, 47, 337–339. [Google Scholar]
  63. The ALICE Collaboration. Anisotropic flow of charged particles in Pb-Pb collisions at  s NN  = 5.02 TeV. Phys. Rev. Lett. 2016, 116, 132302. [Google Scholar] [CrossRef] [PubMed]
  64. Noronha-Hostler, J.; Yan, L.; Gardim, F.G.; Ollitrault, J.Y. Linear and cubic response to the initial eccentricity in heavy-ion collisions. Phys. Rev. C 2016, 93, 014909. [Google Scholar] [CrossRef]
  65. Wei, D.X.; Huang, X.G.; Yan, L. Hydrodynamic response in simulations within a multiphase transport model. Phys. Rev. C 2018, 98, 044908. [Google Scholar] [CrossRef]
  66. Alver, B.; Back, B.B.; Baker, M.D.; Ballintijn, M.; Barton, D.S.; Betts, R.R.; Bickley, A.A.; Bindel, R.; Busza, W.; Carroll, A. System size, energy, pseudorapidity, and centrality dependence of elliptic flow. Phys. Rev. Lett. 2007, 98, 242302. [Google Scholar] [CrossRef] [PubMed]
  67. Cacciari, M.; Salam, G.P.; Soyez, G. The anti-kt jet clustering algorithm. J. High Energy Phys. 2008, 2008, 63. [Google Scholar] [CrossRef]
  68. CMS Collaboration. Studies of Jet Mass in Dijet and W/Z + Jet Events. J. High Energy Phys. 2013, 2013, 90. [Google Scholar] [CrossRef]
  69. Merino, C.; Pajares, C.; Ranft, J. Effects of interaction of strings in the Dual Parton Model. Phys. Lett. B 1992, 276, 168–172. [Google Scholar] [CrossRef]
  70. Mohring, H.J.; Ranft, J.; Merino, C.; Pajares, C. String fusion in the dual parton model and the production of anti-hyperons in heavy ion collisions. Phys. Rev. D 1993, 47, 4142–4145. [Google Scholar] [CrossRef]
  71. Sorge, H.; Berenguer, M.; Stoecker, H.; Greiner, W. Color rope formation and strange baryon production in ultrarelativistic heavy ion collisions. Phys. Lett. B 1992, 289, 6–11. [Google Scholar] [CrossRef]
  72. Bleicher, M.; Greiner, W.; Stoecker, H.; Xu, N. Strangeness enhancement from strong color fields at RHIC. Phys. Rev. C 2000, 62, 061901. [Google Scholar] [CrossRef]
  73. Soff, S.; Randrup, J.; Stoecker, H.; Xu, N. Effects of strong color fields on baryon dynamics. Phys. Lett. B 2003, 551, 115–120. [Google Scholar] [CrossRef]
  74. Amelin, N.S.; Braun, M.A.; Pajares, C. String fusion and particle production at high-energies: Monte Carlo string fusion model. Z. Phys. C 1994, 63, 507–516. [Google Scholar] [CrossRef]
  75. Armesto, N.; Braun, M.A.; Ferreiro, E.G.; Pajares, C. Strangeness enhancement and string fusion in nucleus-nucleus collisions. Phys. Lett. B 1995, 344, 301–307. [Google Scholar] [CrossRef]
  76. Bialas, A.; Czyz, W. Chromoelectric Flux Tubes and the Transverse Momentum Distribution in High-energy Nucleus-nucleus Collisions. Phys. Rev. D 1985, 31, 198. [Google Scholar] [CrossRef] [PubMed]
  77. Kerman, A.K.; Matsui, T.; Svetitsky, B. Particle Production in the Central Rapidity Region of Ultrarelativistic Nuclear Collisions. Phys. Rev. Lett. 1986, 56, 219. [Google Scholar] [CrossRef]
  78. Gyulassy, M.; Iwazaki, A. Quark and gluon pair production in SU(N) covariant constant fields. Phys. Lett. B 1985, 165, 157–161. [Google Scholar] [CrossRef]
  79. Braun, M.; Pajares, C. A Probabilistic model of interacting strings. Nucl. Phys. B 1993, 390, 542–558. [Google Scholar] [CrossRef]
  80. Braun, M.; Pajares, C. Cross-sections and multiplicities in hadron nucleus collisions with interacting color strings. Phys. Rev. D 1993, 47, 114–122. [Google Scholar] [CrossRef]
  81. Bali, G.S. Casimir scaling of SU(3) static potentials. Phys. Rev. D 2000, 62, 114503. [Google Scholar] [CrossRef]
  82. Chen, G.; Fries, R.J.; Kapusta, J.I.; Li, Y. Early Time Dynamics of Gluon Fields in High Energy Nuclear Collisions. Phys. Rev. C 2015, 92, 064912. [Google Scholar] [CrossRef]
  83. Schenke, B.; Tribedy, P.; Venugopalan, R. Event-by-event gluon multiplicity, energy density, and eccentricities in ultrarelativistic heavy-ion collisions. Phys. Rev. C 2012, 86, 034908. [Google Scholar] [CrossRef]
  84. Schenke, B.; Tribedy, P.; Venugopalan, R. Fluctuating Glasma initial conditions and flow in heavy ion collisions. Phys. Rev. Lett. 2012, 108, 252301. [Google Scholar] [CrossRef]
  85. Bautista, I.; Pajares, C.; Ramírez, J.E. String percolation in AA and p+p collisions. Rev. Mex. Fis. 2019, 65, 197–223. [Google Scholar] [CrossRef]
  86. Ramírez, J.E.; Díaz, B.; Pajares, C. Interacting color strings as the origin of the liquid behavior of the quark-gluon plasma. Phys. Rev. D 2021, 103, 094029. [Google Scholar] [CrossRef]
  87. Bierlich, C.; Chakraborty, S.; Gustafson, G.; Lönnblad, L. Strangeness enhancement across collision systems without a plasma. Phys. Lett. B 2022, 835, 137571. [Google Scholar] [CrossRef]
  88. Fischer, N.; Sjöstrand, T. Thermodynamical String Fragmentation. J. High Energy Phys. 2017, 2017, 140. [Google Scholar] [CrossRef]
  89. Argyropoulos, S.; Sjöstrand, T. Effects of color reconnection on tt¯ final states at the LHC. J. High Energy Phys. 2014, 2014, 43. [Google Scholar] [CrossRef]
  90. Ferreres-Solé, S.; Sjöstrand, T. The space–time structure of hadronization in the Lund model. Eur. Phys. J. C 2018, 78, 983. [Google Scholar] [CrossRef]
  91. Bierlich, C.; Chakraborty, S.; Gustafson, G.; Lönnblad, L. Hyperfine splitting effects in string hadronization. Eur. Phys. J. C 2022, 82, 228. [Google Scholar] [CrossRef]
  92. Berges, J.; Floerchinger, S.; Venugopalan, R. Thermal excitation spectrum from entanglement in an expanding quantum string. Phys. Lett. B 2018, 778, 442–446. [Google Scholar] [CrossRef]
  93. Hunt-Smith, N.; Skands, P. String fragmentation with a time-dependent tension. Eur. Phys. J. C 2020, 80, 1073. [Google Scholar] [CrossRef]
  94. Gieseke, S.; Kirchgaeßer, P.; Plätzer, S. Baryon production from cluster hadronisation. Eur. Phys. J. C 2018, 78, 99. [Google Scholar] [CrossRef]
  95. Duncan, C.B.; Kirchgaeßer, P. Kinematic strangeness production in cluster hadronization. Eur. Phys. J. C 2019, 79, 61. [Google Scholar] [CrossRef]
  96. Ghosh, A.; Ju, X.; Nachman, B.; Siodmok, A. Towards a deep learning model for hadronization. Phys. Rev. D 2022, 106, 096020. [Google Scholar] [CrossRef]
  97. Chan, J.; Ju, X.; Kania, A.; Nachman, B.; Sangli, V.; Siodmok, A. Fitting a deep generative hadronization model. J. High Energy Phys. 2023, 2023, 84. [Google Scholar] [CrossRef]
  98. Ilten, P.; Menzo, T.; Youssef, A.; Zupan, J. Modeling hadronization using machine learning. SciPost Phys. 2023, 14, 027. [Google Scholar] [CrossRef]
  99. Bierlich, C.; Ilten, P.; Menzo, T.; Mrenna, S.; Szewc, M.; Wilkinson, M.K.; Youssef, A.; Zupan, J. Towards a data-driven model of hadronization using normalizing flows. arXiv 2023, arXiv:2311.09296. [Google Scholar]
  100. Bierlich, C.; Ilten, P.; Menzo, T.; Mrenna, S.; Szewc, M.; Wilkinson, M.K.; Youssef, A.; Zupan, J. Reweighting Monte Carlo Predictions and Automated Fragmentation Variations in Pythia 8. arXiv 2023, arXiv:2308.13459. [Google Scholar]
Figure 1. The main concept of Angantyr (left) is to build up a heavy ion collision from contributions from each wounded nucleon in the projectile and target. The figure shows how this can be used to build a pp, a pA and an  A A  collision. This treatment, along with emphasis on colour fluctuations in the initial state, provides a good description of forward production (right) (yellow bars in the ratio denotes error ratio), important to determine centrality classes of collisions in the same way as experiments does. In this figure, and throughout the paper, data comparison is carried out using RIVET [35,36]. (Right figure from ref. [25], data from ATLAS [37]).
Figure 1. The main concept of Angantyr (left) is to build up a heavy ion collision from contributions from each wounded nucleon in the projectile and target. The figure shows how this can be used to build a pp, a pA and an  A A  collision. This treatment, along with emphasis on colour fluctuations in the initial state, provides a good description of forward production (right) (yellow bars in the ratio denotes error ratio), important to determine centrality classes of collisions in the same way as experiments does. In this figure, and throughout the paper, data comparison is carried out using RIVET [35,36]. (Right figure from ref. [25], data from ATLAS [37]).
Universe 10 00046 g001
Figure 2. Particle production at mid-rapidity from the Angantyr model, compared to data in pPb collisions (left) and PbPb (right) in centrality bins. (Figures from ref. [25], data from ATLAS and ALICE [37,38]).
Figure 2. Particle production at mid-rapidity from the Angantyr model, compared to data in pPb collisions (left) and PbPb (right) in centrality bins. (Figures from ref. [25], data from ATLAS and ALICE [37,38]).
Universe 10 00046 g002
Figure 3. In pp collisions at  s  = 7 TeV (left), many strings can overlap in transverse space and rapidity as shown, giving rise to interactions and coherence effects. Note that the string radius is here put at 0.2 fm for illustrative purposes, but can in reality be much larger. The string’s transverse (colour)-electric field itself is not a cylinder, but can rather be approximated as a Gaussian (right), with black points being lattice calculations [39]. (Figures from refs. [40,41]).
Figure 3. In pp collisions at  s  = 7 TeV (left), many strings can overlap in transverse space and rapidity as shown, giving rise to interactions and coherence effects. Note that the string radius is here put at 0.2 fm for illustrative purposes, but can in reality be much larger. The string’s transverse (colour)-electric field itself is not a cylinder, but can rather be approximated as a Gaussian (right), with black points being lattice calculations [39]. (Figures from refs. [40,41]).
Universe 10 00046 g003
Figure 4. Results from QCD-CR colour reconnection in pp (left) and pPb (right), showing that the junction production mechanism drastically improves the description of  Λ c / D  ratios in both collision systems.(Figures from ref. [57], data from ALICE [56,58]).
Figure 4. Results from QCD-CR colour reconnection in pp (left) and pPb (right), showing that the junction production mechanism drastically improves the description of  Λ c / D  ratios in both collision systems.(Figures from ref. [57], data from ALICE [56,58]).
Universe 10 00046 g004
Figure 5. The string-shoving model reproducing the ridge in pp collisions in minimum bias collisions in different multiplicity bins (left), and in Z-tagged, high multiplicity collisions (right). (Figures from ref. [48,49], data from CMS and ATLAS [13,61]).
Figure 5. The string-shoving model reproducing the ridge in pp collisions in minimum bias collisions in different multiplicity bins (left), and in Z-tagged, high multiplicity collisions (right). (Figures from ref. [48,49], data from CMS and ATLAS [13,61]).
Universe 10 00046 g005
Figure 6. Shoving model results in  A A  collisions for initial toy conditions (see text). Comparison to  v 2  in PbPb collisions at  s N N = 5.02  TeV (left) shows a remarkable agreement with data, and studies of the correlation of  v 2  vs. the eccentricity  ϵ 2  (right) shows a similar response as one would expect from hydrodynamic calculations. (Figures from ref. [40], data from ALICE [63]).
Figure 6. Shoving model results in  A A  collisions for initial toy conditions (see text). Comparison to  v 2  in PbPb collisions at  s N N = 5.02  TeV (left) shows a remarkable agreement with data, and studies of the correlation of  v 2  vs. the eccentricity  ϵ 2  (right) shows a similar response as one would expect from hydrodynamic calculations. (Figures from ref. [40], data from ALICE [63]).
Universe 10 00046 g006
Figure 7. The effect on string shoving on jet cross section (details in text) vs. jet radius (left) shows a small effect, similar to that of hadronization. On observables such as jet mass (right) the effect can reach the 10–20% level for small jet masses, out of reach of current experiments. (Figures from ref. [49], data from CMS [68]).
Figure 7. The effect on string shoving on jet cross section (details in text) vs. jet radius (left) shows a small effect, similar to that of hadronization. On observables such as jet mass (right) the effect can reach the 10–20% level for small jet masses, out of reach of current experiments. (Figures from ref. [49], data from CMS [68]).
Universe 10 00046 g007
Figure 8. Production points of hadrons in impact parameter space for representative collisions of pp at 7 TeV (top row), pPb at 5.02 TeV (middle row) and PbPb at 2.76 TeV (bottom row). The hadron production points are coloured according to the rope tension, giving an image of the energy density reached in the ropes. (Figure from ref. [87]).
Figure 8. Production points of hadrons in impact parameter space for representative collisions of pp at 7 TeV (top row), pPb at 5.02 TeV (middle row) and PbPb at 2.76 TeV (bottom row). The hadron production points are coloured according to the rope tension, giving an image of the energy density reached in the ropes. (Figure from ref. [87]).
Universe 10 00046 g008
Figure 9. Production of multi-strange hadrons with rope hadronization in pp, pPb and PbPb. (Figure from ref. [87]. Data from ALICE [12]).
Figure 9. Production of multi-strange hadrons with rope hadronization in pp, pPb and PbPb. (Figure from ref. [87]. Data from ALICE [12]).
Universe 10 00046 g009
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

Bierlich, C. String Interactions as a Source of Collective Behaviour. Universe 2024, 10, 46. https://doi.org/10.3390/universe10010046

AMA Style

Bierlich C. String Interactions as a Source of Collective Behaviour. Universe. 2024; 10(1):46. https://doi.org/10.3390/universe10010046

Chicago/Turabian Style

Bierlich, Christian. 2024. "String Interactions as a Source of Collective Behaviour" Universe 10, no. 1: 46. https://doi.org/10.3390/universe10010046

APA Style

Bierlich, C. (2024). String Interactions as a Source of Collective Behaviour. Universe, 10(1), 46. https://doi.org/10.3390/universe10010046

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop