Rope Hadronization in Dense QCD

• Rope hadronization is a phenomenological extension of the Lund string model that combines overlapping color strings into higher SU(3) multiplets with enhanced effective tension.
• It modifies tunneling probabilities to reduce suppression for s‑sbar pairs and diquarks, thereby explaining strangeness and baryon enhancement in high-multiplicity collisions.
• Implemented in DIPSY, PYTHIA8, and PYTHIA/Angantyr, RH offers a coherent framework to interpret collective-like dynamics and near‐side ridge phenomena without invoking a thermalized medium.

Searching arXiv for recent and foundational papers on rope hadronization and string shoving. arXiv search: rope hadronization string shoving PYTHIA 8 Rope hadronization (RH) is a phenomenological extension of the Lund string hadronization model for dense QCD environments in which several nearby color strings overlap in transverse space and act coherently as a higher SU(3) color multiplet, or color rope, with an enhanced effective string tension. In this framework, the fragmentation field strength is no longer that of an isolated $q\bar q$ string; reduced tunneling suppression for $s\bar s$ pairs and diquarks modifies hadrochemistry, while companion mechanisms such as string shoving can convert overlap energy into collective-like transverse dynamics. RH has been implemented in DIPSY, PYTHIA 8, and PYTHIA/Angantyr, and has been used to interpret multiplicity-dependent strangeness enhancement, baryon production, resonance systematics, and near-side ridge-like correlations in small and large collision systems (Bierlich et al., 2014, Bierlich, 2017, Bierlich, 2018).

1. String-fragmentation basis and color-rope formalism

In the baseline Lund picture, a color flux tube between separating colored partons carries approximately constant tension,

$\kappa \approx 1~\text{GeV/fm},$

and string breaking follows a Schwinger-like tunneling law,

$$\frac{d\mathcal{P}}{dp_\perp} \propto \kappa \exp\!\left(-\frac{\pi m_\perp^2}{\kappa}\right).$$

The corresponding strange-quark suppression factor is written as

$$\rho = \frac{\mathcal{P}_s}{\mathcal{P}_{u,d}} = \exp\!\left[-\frac{\pi\left(m_s^2-m_u^2\right)}{\kappa}\right].$$

RH modifies this independent-string approximation in events where many strings are produced close together and overlap in transverse space, so that the combined color field behaves coherently rather than as separate singlets (Bierlich, 2016, Bierlich, 2017).

The rope is represented as an SU(3) multiplet labeled by $\{p,q\}$. In the random-walk formulation, adding a triplet to a multiplet gives

$$\{p,q\}\oplus\{1,0\}= \left\{ \{p+1,q\},\; \{p,q-1\},\; \{p-1,q+1\} \right\},$$

with probabilities proportional to the multiplet dimensions,

$N_{\{p,q\}}=\frac{(p+1)(q+1)(p+q+2)}{2}.$

The full rope tension scales with the quadratic Casimir,

$$\frac{\kappa^{\{p,q\}}}{\kappa^{\{1,0\}}} = \frac{1}{4}\left(p^2+pq+q^2+3p+3q\right),$$

while the tension relevant for a single breakup is the reduction in rope tension after removing one color charge. In the formulations used in DIPSY and PYTHIA, this gives

$\tilde{\kappa}=\frac{2p+q+2}{4}\,\kappa.$

Replacing $s\bar s$0 by $s\bar s$1 in the tunneling probabilities yields

$s\bar s$2

so increasing overlap drives $s\bar s$3 toward unity and enhances $s\bar s$4 and diquark production. This is the central microscopic origin of the predicted enhancement of strange hadrons and baryons in RH (Bierlich et al., 2014, Bierlich et al., 2022).

2. Relation to color reconnection and string shoving

RH is closely related to, but distinct from, two other nonperturbative ingredients used in contemporary event-generator studies: color reconnection (CR) and string shoving. The distinction is structural. CR changes the topology of color connections, typically by minimizing total string length or potential energy; RH changes the effective fragmentation field strength by merging overlapping strings into ropes; string shoving adds a dynamical transverse pressure generated by the excess energy in overlapping strings (Hushnud et al., 2023, Bierlich et al., 2016).

Mechanism	What it changes	Typical phenomenology
Color reconnection	String topology	Junction production, resonance suppression, $s\bar s$5 trends
Rope hadronization	Effective string tension	Strangeness and baryon enhancement
String shoving	Transverse pressure between strings	Mass ordering and ridge-like correlations

This separation is operationally important. CR is described as a rearrangement of color topology, while RH is explicitly described as not just a rearrangement of color topology. In the PYTHIA8 studies of strange-particle production, the complementarity is stated directly: CR changes the string topology, RH changes the effective fragmentation strength (Hushnud et al., 2023). In the shoving model, the transverse field profile is taken as

$s\bar s$6

leading to a repulsive force per unit length

$s\bar s$7

so the same spatial overlap that later produces ropes can earlier generate collective transverse kicks (Bierlich, 2018).

The implementations reflect this decomposition. Early pp studies used DIPSY, which provides transverse-coordinate information and interfaces to PYTHIA 8 for fragmentation. Later work introduced a direct PYTHIA8 rope-hadronization plug-in, and heavy-ion extensions were built in PYTHIA/Angantyr. For nuclei and nontrivial string topologies, the parallel-frame formalism constructs a common Lorentz frame for each pair of string pieces and estimates overlap from spacetime location, relative angle, and Gaussian transverse string profiles; string pieces with gluon kinks are included. The same geometric framework is also used in newer PYTHIA8 developments to evaluate both shoving and rope formation for all string pairs, including strings from high-$s\bar s$8 partons and jets (Bierlich, 2017, Bierlich et al., 2022, Chakraborty, 2021).

3. Strangeness enhancement and baryochemistry across systems

The most established RH phenomenology concerns strange and multi-strange hadrons. The generic prediction is a monotonic rise of strange-hadron-to-pion ratios with event activity, with stronger enhancement for hadrons containing more strange quarks. In pp collisions at $s\bar s$9 and $\kappa \approx 1~\text{GeV/fm},$0 TeV, PYTHIA8 studies with RH reported successful descriptions of $\kappa \approx 1~\text{GeV/fm},$1, $\kappa \approx 1~\text{GeV/fm},$2, $\kappa \approx 1~\text{GeV/fm},$3, $\kappa \approx 1~\text{GeV/fm},$4, $\kappa \approx 1~\text{GeV/fm},$5, and integrated yields $\kappa \approx 1~\text{GeV/fm},$6, together with multiplicity-dependent hardening of the strange-hadron $\kappa \approx 1~\text{GeV/fm},$7 spectra. Those studies also emphasized that the enhancement is controlled primarily by event multiplicity rather than by the center-of-mass energy, and that the multiplicity dependence appears to saturate at higher multiplicities (Nayak et al., 2018).

A more differential comparison at $\kappa \approx 1~\text{GeV/fm},$8 TeV separated the roles of CR and RH. In that study, the MPI-based CR model described $\kappa \approx 1~\text{GeV/fm},$9 reasonably well but underestimated strange baryons; QCD-based CR improved $$\frac{d\mathcal{P}}{dp_\perp} \propto \kappa \exp\!\left(-\frac{\pi m_\perp^2}{\kappa}\right).$$0 through junction formation but still failed for $$\frac{d\mathcal{P}}{dp_\perp} \propto \kappa \exp\!\left(-\frac{\pi m_\perp^2}{\kappa}\right).$$1 and $$\frac{d\mathcal{P}}{dp_\perp} \propto \kappa \exp\!\left(-\frac{\pi m_\perp^2}{\kappa}\right).$$2; the combination of RH with QCD-based CR gave the best description of $$\frac{d\mathcal{P}}{dp_\perp} \propto \kappa \exp\!\left(-\frac{\pi m_\perp^2}{\kappa}\right).$$3, $$\frac{d\mathcal{P}}{dp_\perp} \propto \kappa \exp\!\left(-\frac{\pi m_\perp^2}{\kappa}\right).$$4, and $$\frac{d\mathcal{P}}{dp_\perp} \propto \kappa \exp\!\left(-\frac{\pi m_\perp^2}{\kappa}\right).$$5, while overestimating $$\frac{d\mathcal{P}}{dp_\perp} \propto \kappa \exp\!\left(-\frac{\pi m_\perp^2}{\kappa}\right).$$6 at high multiplicity. The same analysis found that RH was essential for yields and baryon-to-meson ratios, but that the multiplicity dependence of $$\frac{d\mathcal{P}}{dp_\perp} \propto \kappa \exp\!\left(-\frac{\pi m_\perp^2}{\kappa}\right).$$7 was described best by MPI-based CR with reconnection range $$\frac{d\mathcal{P}}{dp_\perp} \propto \kappa \exp\!\left(-\frac{\pi m_\perp^2}{\kappa}\right).$$8; RH alone underestimated $$\frac{d\mathcal{P}}{dp_\perp} \propto \kappa \exp\!\left(-\frac{\pi m_\perp^2}{\kappa}\right).$$9 and was reported to “totally fail” for its multiplicity dependence. A plausible implication is that RH is most robust for hadrochemistry, whereas collective hardening of spectra is more configuration-dependent (Hushnud et al., 2023).

The cross-system extension in PYTHIA/Angantyr found a coherent increase of strange-hadron and baryon yields with average charged central multiplicity in pp, pPb, and PbPb, in qualitative agreement with LHC data. The same mechanism therefore provides a continuous multiplicity-based description across systems without switching to a plasma picture. The principal limitation in that study was that high-multiplicity PbPb events showed an overshoot of baryonic yields, and in some cases too steep a rise for kaons, indicating that RH alone does not saturate the enhancement as strongly as the data (Bierlich et al., 2022).

Event-topology studies have refined the multiplicity picture. In PYTHIA 8 with CR and RH at $$\rho = \frac{\mathcal{P}_s}{\mathcal{P}_{u,d}} = \exp\!\left[-\frac{\pi\left(m_s^2-m_u^2\right)}{\kappa}\right].$$0 TeV, strange-to-pion ratios were analyzed versus $$\rho = \frac{\mathcal{P}_s}{\mathcal{P}_{u,d}} = \exp\!\left[-\frac{\pi\left(m_s^2-m_u^2\right)}{\kappa}\right].$$1, midrapidity and forward multiplicity, transverse sphericity, transverse spherocity, unweighted transverse spherocity, relative transverse activity, and charged-particle flattenicity. The reported hierarchy was

$$\rho = \frac{\mathcal{P}_s}{\mathcal{P}_{u,d}} = \exp\!\left[-\frac{\pi\left(m_s^2-m_u^2\right)}{\kappa}\right].$$2

with $$\rho = \frac{\mathcal{P}_s}{\mathcal{P}_{u,d}} = \exp\!\left[-\frac{\pi\left(m_s^2-m_u^2\right)}{\kappa}\right].$$3 nearly flat and $$\rho = \frac{\mathcal{P}_s}{\mathcal{P}_{u,d}} = \exp\!\left[-\frac{\pi\left(m_s^2-m_u^2\right)}{\kappa}\right].$$4 the most sensitive. Charged-particle flattenicity was identified as the most suitable classifier because it reproduced a quantitative enhancement similar to that obtained with the number of multi-parton interactions, whereas some other selections were strongly affected by autocorrelation bias or limited multiplicity reach (Prasad et al., 2024).

4. Collective-like dynamics and the near-side ridge

Original RH mainly modified hadrochemistry through $$\rho = \frac{\mathcal{P}_s}{\mathcal{P}_{u,d}} = \exp\!\left[-\frac{\pi\left(m_s^2-m_u^2\right)}{\kappa}\right].$$5. The string-shoving extension added a dynamically generated transverse pressure arising from the excess energy in the overlap region of nearby strings. In that picture, strings first overlap and repel each other, then later hadronize as ropes. The combined model was designed to reproduce soft features usually associated with hydrodynamics—higher $$\rho = \frac{\mathcal{P}_s}{\mathcal{P}_{u,d}} = \exp\!\left[-\frac{\pi\left(m_s^2-m_u^2\right)}{\kappa}\right].$$6 for heavier particles and long-range azimuthal correlations—without assuming a thermalized medium (Bierlich et al., 2016).

The ridge application was made explicit in a PYTHIA 8 study of two-particle $$\rho = \frac{\mathcal{P}_s}{\mathcal{P}_{u,d}} = \exp\!\left[-\frac{\pi\left(m_s^2-m_u^2\right)}{\kappa}\right].$$7–$$\rho = \frac{\mathcal{P}_s}{\mathcal{P}_{u,d}} = \exp\!\left[-\frac{\pi\left(m_s^2-m_u^2\right)}{\kappa}\right].$$8 correlations in pp collisions at $$\rho = \frac{\mathcal{P}_s}{\mathcal{P}_{u,d}} = \exp\!\left[-\frac{\pi\left(m_s^2-m_u^2\right)}{\kappa}\right].$$9 and $\{p,q\}$0 TeV. The analysis used same-event and mixed-event pair distributions, with 10 mixed events, and employed the ZYAM subtraction procedure for $\{p,q\}$1 projections. In high-multiplicity events, an additional long-range near-side enhancement extending over roughly three units in $\{p,q\}$2 was observed only when RH with shoving was enabled. The structure was absent in low-multiplicity events, and it was also absent in high-multiplicity events when rope formation was not included. The study further reported that the correlation strength was larger for CR-1 and CR-2 than for CR-0, but that color reconnection alone was not sufficient to generate the long-range near-side ridge. The stated interpretation was that microscopic final-state string dynamics can reproduce the qualitative ridge feature without requiring a deconfined thermal medium (Chakraborty et al., 2020).

This division of labor is consistent across the RH literature. RH at hadronization primarily changes hadrochemistry, especially strange-particle yields, whereas shoving acts earlier and generates transverse kicks that can produce mass ordering in $\{p,q\}$3 and long-range azimuthal structure. The two mechanisms are therefore usually treated as coupled but non-identical components of a string-interaction picture of small-system collectivity (Bierlich, 2018).

5. Resonances, heavy flavor, and multiplicity-dependent yields

RH has also been tested outside light-flavor strange-hadron ratios. In heavy flavor, PYTHIA 8 studies of pp collisions at $\{p,q\}$4 and $\{p,q\}$5 TeV examined the relative yields of $\{p,q\}$6 mesons, $\{p,q\}$7, and predicted $\{p,q\}$8-meson yields as functions of relative charged-particle multiplicity for several $\{p,q\}$9 intervals. The reported pattern was a stronger-than-linear rise of open-charm yields with multiplicity, with the best qualitative agreement obtained when RH was combined with CR. The same framework reproduced the increasing trend of inclusive $$\{p,q\}\oplus\{1,0\}= \left\{ \{p+1,q\},\; \{p,q-1\},\; \{p-1,q+1\} \right\},$$0 yields but underestimated the magnitude, and it predicted a similar non-linear rise for $$\{p,q\}\oplus\{1,0\}= \left\{ \{p+1,q\},\; \{p,q-1\},\; \{p-1,q+1\} \right\},$$1 mesons in $$\{p,q\}\oplus\{1,0\}= \left\{ \{p+1,q\},\; \{p,q-1\},\; \{p-1,q+1\} \right\},$$2. The interpretation emphasized microscopic partonic production mechanisms—especially MPI and hard heavy-quark production—rather than collective thermal-medium effects (Tripathy et al., 2022).

Resonance production reveals a more differentiated interplay between CR and RH. In PYTHIA 8 studies at $$\{p,q\}\oplus\{1,0\}= \left\{ \{p+1,q\},\; \{p,q-1\},\; \{p-1,q+1\} \right\},$$3 and $$\{p,q\}\oplus\{1,0\}= \left\{ \{p+1,q\},\; \{p,q-1\},\; \{p-1,q+1\} \right\},$$4 TeV, mesonic resonance-to-stable-hadron ratios such as $$\{p,q\}\oplus\{1,0\}= \left\{ \{p+1,q\},\; \{p,q-1\},\; \{p-1,q+1\} \right\},$$5, $$\{p,q\}\oplus\{1,0\}= \left\{ \{p+1,q\},\; \{p,q-1\},\; \{p-1,q+1\} \right\},$$6, and $$\{p,q\}\oplus\{1,0\}= \left\{ \{p+1,q\},\; \{p,q-1\},\; \{p-1,q+1\} \right\},$$7 were suppressed at high event activity primarily by QCD-based CR, which forms shorter strings and biases fragmentation toward lighter hadrons. The main RH signatures appeared in channels with explicit strange content: $$\{p,q\}\oplus\{1,0\}= \left\{ \{p+1,q\},\; \{p,q-1\},\; \{p-1,q+1\} \right\},$$8 and $$\{p,q\}\oplus\{1,0\}= \left\{ \{p+1,q\},\; \{p,q-1\},\; \{p-1,q+1\} \right\},$$9 were enhanced at high multiplicity when RH was included, reflecting increased production of strange and anti-strange quark pairs. In the baryon sector, ratios such as $N_{\{p,q\}}=\frac{(p+1)(q+1)(p+q+2)}{2}.$0, $N_{\{p,q\}}=\frac{(p+1)(q+1)(p+q+2)}{2}.$1, and $N_{\{p,q\}}=\frac{(p+1)(q+1)(p+q+2)}{2}.$2 were enhanced especially with CR(1), with RH further strengthening the effect in strange baryonic channels. An important caveat in that study was that RH without CR had no significant effect on the studied resonance ratios, implying that CR dominated the suppression and enhancement patterns while RH supplied the strangeness-driven component (Goswami et al., 2019).

These extensions broaden the scope of RH beyond standard strange-hadron systematics. They also show that its phenomenological role is not uniform: in some observables RH is the decisive ingredient, while in others it is subleading to color-topology effects or to the hard-production sector. This suggests that RH is best viewed as one component of a coupled MPI–CR–hadronization framework rather than as a standalone explanation for all multiplicity-dependent phenomena.

6. Interpretation, limitations, and ongoing issues

RH is often discussed because it challenges a common inference: that strangeness enhancement and some collectivity-like signals uniquely imply quark–gluon plasma formation. Across pp, pPb, and PbPb studies, RH has been shown to generate rising strange-hadron and baryon yields as a function of multiplicity without invoking a deconfined thermal medium. The same line of work argues that ridge-like structures in small systems can arise from microscopic string interactions when shoving is included. At the same time, the RH papers do not present this as a proof that QGP-like dynamics are absent; one study states explicitly that RH shows QGP is not required to explain particular observables, but does not exclude other collective effects (Nayak et al., 2018, Bierlich et al., 2022).

The principal limitations are empirical and implementation-dependent. No single PYTHIA8 configuration was found to simultaneously and quantitatively describe both strangeness enhancement and the multiplicity-dependent rise of $N_{\{p,q\}}=\frac{(p+1)(q+1)(p+q+2)}{2}.$3; RH helps yields and ratios, whereas MPI-based CR with large $N_{\{p,q\}}=\frac{(p+1)(q+1)(p+q+2)}{2}.$4 gives the best $N_{\{p,q\}}=\frac{(p+1)(q+1)(p+q+2)}{2}.$5 description in the 13 TeV strange-hadron study (Hushnud et al., 2023). In high-multiplicity PbPb, RH overshoots baryonic yields and does not saturate the enhancement as strongly as the data, motivating combinations with color reconnection between sub-collisions, string shoving, and eventually hadronic rescattering (Bierlich et al., 2022).

Model dependence also enters through overlap geometry and tuning. A review of hadronisation models and color reconnection emphasized that rope effects are very sensitive to the modeling of transverse proton structure. Later PYTHIA8/Angantyr work noted that default fragmentation parameters tuned to LEP data are no longer fully consistent once RH and the strange-quark magnetic-moment contribution are included, and that some observables, notably $N_{\{p,q\}}=\frac{(p+1)(q+1)(p+q+2)}{2}.$6 in pp, can be overshot in current implementations (Bierlich, 2016, Chakraborty, 2021). Related phenomenological work on thermodynamical string fragmentation with environment-dependent effective tension or temperature explicitly described those modifications as proxies for close-packed strings rather than as a full microscopic RH implementation, and concluded that such proxies were not sufficient to generate ridge physics or full collective flow (Fischer et al., 2016).

Within contemporary event-generator phenomenology, RH therefore occupies a specific position. It is a microscopic final-state string-interaction mechanism that robustly modifies hadrochemistry in dense events, especially for strange and multi-strange hadrons, and that can participate in a broader explanation of collectivity-like observables when coupled to shoving and CR. Its successes are substantial but selective, and its limitations are sufficiently well documented that RH is best treated as a nontrivial alternative—and sometimes complement—to hydrodynamic interpretations rather than as a complete substitute for them.