String interactions as a source of collective behaviour

The discovery of collective effects in small collision systems, have spurred a renewed interest in hadronization models, also as 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 is reviewed, and some prospects for the future presented.


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 observations 1 of strangeness enhancement relative to the proton-proton baseline [1,2] and anisotropic flow [3,4].Over the past decades, a rich research programme has emerged, 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 AA) 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-review 2 .
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 sec. 2 the underlying Angantyr model is briefly described, before diving into the string formalism in sec.3.Here interactions of multi-string systems is treated through sec.3.2, before some concluding remarks and an outlook is given in sec.4.
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.[18,19,20,21,22,14] 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.

Projectile
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 AA collision.This treatment, along with emphasis on colour fluctuations in the initial state, provides a good description of forward production (right), important to determine centrality classes of collisions in the same way as experiments does.In this figure, and throughout the paper, data comparison is done using RIVET [23,24].(Right figure from ref. [25], data from ATLAS [26])

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 [27] 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]: 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 with 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 eq. ( 1) as: ( This cut-off is usually interpreted as colour screening scale, though suggestions linking it to saturation in the spirit of Color Glass Condensate (CGC) [28] instead, have been put forth.
The extension of this model, Angantyr [29,25], 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 [30,31] with particular emphasis on colour fluctuations in the initial state [32,33,34].Through the Good-Walker formalism [35] 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 through 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 [37].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 AA collision.The process is sketched in fig. 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 a very different starting point.This is a direction of research currently left very unexplored.
As shown in fig. 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 AA [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 fig. 2 (left), and for PbPb in fig. 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 components 4 of the Lund string model for a single string will be introduced in sec.3.1, proceeding to systems of multiple strings in sec.3.2.Colour reconnection is presented first in sec.3.3, with string shoving and rope hadronization presented in sec.3.4 and 3.5, respectively.

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.√ 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]) 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: 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: dP Due to the factorization in the above expression, the transverse momenta and mass of the quarks can be treated separately.In practice, this is done 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 2 s − m 2 u )/κ), which is also treated as a parameter6 .

Interactions of multiple strings
In collisions of protons (and thus also in larger collision systems), strings will overlap with each other (see fig. 3, left).It is therefore reasonable to question whether the initial assumption of strings being infinitely thin should be relaxed.Discarding the assumption alltogether, 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 fig. 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 eq. ( 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: which, with reasonable values for a and b (determined from total multiplicities in e + e − using eq.( 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 later 8 .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].
• τ ≈ 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 [40].
• After τ ≈ 1.4 fm: Hadronic cascade.Hadronic rescattering [51] results in significant corrections to both flavour and flow, particularly in the dense environment of AA collisions [43].

Colour reconnections
The colour topology delivered by the matrix element calculation and parton shower, uses the "leading colour" limit 9 , 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 done in full colour.The simplest reconnection models aims at reintroducing the singlet.This is for example the case for the default Pythia reconnection model [27] (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 [56].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 element.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 ALICE Data Pythia8 (def) + QCDCR (mode-0) + SC-CR (all)   [53], data from ALICE [54,55]) be written as e.g.λ q q = ln(m 2 q q/m 2 0 ), but for more complicated topologies approximate expressions must be used, as e.g.: 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 states 10 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 [57].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. [58]), but crucially also allows for production of double heavy flavour baryons, not allowed in the normal string model.In particular production of heavy flavour baryons have been a large success for the model, where the Λ c /D ratio vs p ⊥ in pp collisions (see fig. 4 (left)) was found by ALICE [54] to differ drastically from expectation.The colour reconnection model gives a satisfactory description.In ion collisions later refinements of the model [53], taking into account the non-negligible mass of the charm quark, allowed predictions also for pA, as shown in fig. 4 (right).Also here the description is drastically improved.
The last word is, however, not said for this type of models.Further studies of charm baryon yields have revealed that for example charmed Ξ is not equally well reproduced [53,59].Some is on the account of further enhancement of strangeness (see sec.3.5), but not all can be explained this way.Colour reconnection remains to be an active area of study.

The string shoving mechanism
As explained in the beginning of sec.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 /2R 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.Remaining string dipoles will therefore only repel each other, 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: This is the main ingredient of the string shoving mechanism.Technicals complications arise when implementing this in a Monte Carlo.String pieces are not aligned in parallel, but a suitable Lorentz transformation can be found to align each pair in parallel planes [41].Also, simply doing a step-wise time evolution, is not practical from a computational point of view.Instead the shoving model is implemented in similar spirit as a parton shower, ordered in the p ⊥ pushes.

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,   48,49], data from CMS and ATLAS [13,61].) i.e. fully symmetric.Higher order flow is therefore produced only from fluctuations in the initial conditions, and not from intrinsic asymmetry 11 .Overall observables, such as the long-range ridge first measured in pp by CMS [13], is well reproduced as shown in fig. 5.
(Note that the parameter g in the figure is not identical to g in eq. ( 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 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 fig. 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 AA collisions as well 12 .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 [41].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 AA collisions in a given centrality interval, can be set up.In ref. [41] this was done, along with spatial initial conditions resembling those of the 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. [41], data from ALICE [63]) given centrality interval.Furthermore, since the initial conditions are known in this case, flow does not have to be calculated in the same way as experiment 13 , but can be calculated wrt. the true reaction plane, as well as correlations with the initial geometry.In fig.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 fig.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 follows that flow signals cannot be taken as proof of a QGP, indeed not if the same signal can be generated from the shoving model without QGP.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 [67])

Effects on jets
As indicated above, the shoving model applies everywhere there are two or more strings available, irrespective of collision system or hard process.Since jet quenching in AA collisions constitute 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 ⊥,jet > 80 GeV, the jet cross section σ j = ∞ p ⊥,0 dp ⊥,j dσ dp ⊥,j is calculable in perturbative QCD, but receives large corrections from MPI and hadronization.In fig.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 few percent level.The main point of the figure is twofold.First of old, 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 in the underlying event, where even basic quantities like the jet cross section is affected.Second, that measurements of jet modifications due to so 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 fig.7 (right), compared with measurements by CMS at √ s = 7 TeV.Sizeable effects (over 10%) show up at low jet masses, less than 20 GeV, but current measurements have no chance of distinguishing.

Rope hadronization
After the shoving phase, 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" [40].Connecting back to the explanations on colour multiplets given in sec.3.3, this is the effect that the string tension will be enhanced for higher multiplets 15SU(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: 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: The energy density in a rope can quickly grow large enough to make a dramatic impact on final state observables (see sec.3.5.1).It will, however, not exceed more than around a factor 2 above the vacuum string tension, except in very central AA 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 (dE/d 3 x), which for example in CGC based models can reach over 100 GeV/fm 3 in PbPb collisions at LHC [83,84,41], are so large that plasma creation is inevitable.Crucially the energy density in a rope scenario never reaches such levels, as shown in fig.8. Other models of string interactions, crucially the string perculation model [85], has the opposite conclusion: that the percolation of strings is potentially the mechanism responsible for QGP formation [86].

Effects on inclusive production
When the string tension rises, the main effect is that suppression of strangeness in eq. ( 4) decreases.The derived quantity ρ is replaced by an effective quantity ρ, which can be obtained directly as: where the exponent can be obtained directly from (p, q) through eq. ( 10).Similar expressions exist for the other parameters of the hadronization model, see refs.[40,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 AA collision is 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 fig. 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.

Conclusion 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 eq. ( 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 proofs-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.

Figure 3 :
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 4 :
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.(Figuresfrom ref.[53], data from ALICE [54, 55])

Figure 5 :
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].)

1 2 v 2 (Figure 6 :
Figure 6: Shoving model results in AA collisions, for toy initial 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.[41], data from ALICE [63])

Figure 7 :
Figure7: 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[67])

Figure 8 :
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].)

1
Average ratio of Λ + c and D 0