# QCD challenges from pp to AA collisions – 4th edition

Javira Altmann, Carlota Andres, Anton Andronic, Federico Antinori, Pietro Antonioli, Andrea Beraudo, Eugenio Berti, Livio Bianchi, Thomas Boettcher, Lorenzo Capriotti, et al.

► To cite this version:

Javira Altmann, Carlota Andres, Anton Andronic, Federico Antinori, Pietro Antonioli, et al., QCD challenges from pp to AA collisions – 4th edition. European Physical Journal C: Particles and Fields, 2024, 84 (4), pp.421. 10.1140/epjc/s10052-024-12650-8. hal-04425529

HAL Id: hal-04425529

https://hal.science/hal-04425529v1

Submitted on 15 Jan 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
QCD challenges from pp to AA collisions: 4th edition

Javira Altmann¹, Carlota Andres², Anton Andronic³, Federico Antinori⁴, Pietro Antonioli⁵, Andrea Beraudo⁶,
Eugenio Bertì⁷, Livio Bianchi⁶,⁸, Thomas Boettcher⁹, Lorenzo Capriotti¹⁰, Peter Christiansen¹¹,
Jesus Guillermo Contreras Nuño¹², Leticia Cunqueiro Mendez¹³, Cesar da Silva¹⁴, Andrea Dainese⁴,
Hans Peter Dembinski¹⁵, David Dobrigkeit Chinellato¹⁶,¹⁷, Andrea Dubla¹⁸, Mattia Faggin¹⁹,²⁰, Chris Flett²¹,
Vincenzo Greco²²,²³, Ilia Grishmanovskii²⁴, Jack Holguin², Yuuka Kanakubo²⁵, Dong Jo Kim²⁵, Ramona Lea²⁶,
Su Houng Lee²⁷, Saverio Mariani²⁸, Adam Matyja²⁹, Aleksas Mazeliauskas³⁰, Vincenzo Minissale²²,³¹,
Andreas Morsch²⁸, Lucia Oliva²²,³¹, Luca Orusa⁶,⁸, Petja Paakkinen²⁵, Daniel Pablos⁶,³²,³³, Guy Paić³⁴,
Tanguy Pierog³⁵, Salvatore Plumari²²,²³, Francesco Prino⁶, Andrea Rossi⁴,ª, Lorenzo Sestini⁴,ᵇ, Peter Skands¹,
Olga Soloveva³⁶, Francesca Soramel⁴,¹⁹, Alba Soto Ontoso²⁸, Martin Spousta³⁷, Andre Govinda Stahl Leiton²⁸,
Jiayin Sun³⁸, Adam Takacs³⁹, Stefano Trogolo⁸,⁴⁰, Rosario Turrisi⁴, Marta Verweij⁴¹, Vytautas Vislavicius¹¹,
Jing Wang⁴², Klaus Werner⁴³, Valentina Zaccolo²⁰, Mingyu Zhang⁴⁴, Jianhui Zhu⁴, Davide Zuliani⁴,¹⁹

¹ Monash University, Melbourne, Australia

² CPHT, CNRS, Ecole polytechnique, Institut Polytechnique de Paris, Palaiseau, France

³ Universität Münster, Münster, Germany

⁴ INFN, Sezione di Padova, Padua, Italy

⁵ Università and INFN Bologna, Bologna, Italy

⁶ INFN, Sezione di Torino, Turin, Italy

⁷ INFN, Sezione di Firenze, Florence, Italy

⁸ Università di Torino, Turin, Italy

⁹ University of Cincinnati, Cincinnati, USA

¹⁰ Università and INFN Ferrara, Ferrara, Italy

¹¹ Lund University, Lund, Sweden

¹² Czech Technical University in Prague, Prague, Czech Republic

¹³ Università Sapienza, Rome, Italy

¹⁴ Los Alamos National Lab, Santa Fe, USA

¹⁵ Technische Universität Dortmund, Dortmund, Germany

¹⁶ University of Campinas UNICAMP, Campinas, Brazil

¹⁷ Stefan Meyer Institute of the Austrian Academy of Sciences, Vienna, Austria

¹⁸ GSI Helmholtzzentrum für Schwerionenforschung, Darmstadt, Germany

¹⁹ Università di Padova, Padua, Italy

²⁰ Università and INFN Trieste, Trieste, Italy

²¹ Université Paris-Saclay, CNRS, IJCLab, Orsay, France

²² Università di Catania, Catania, Italy

²³ INFN-LNS, Catania, Italy

²⁴ Institut für Theoretische Physik, Frankfurt, Germany

²⁵ University of Jyväskylä and Helsinki Institute of Physics, Helsinki, Finland

²⁶ Universita di Brescia, Brescia, Italy

²⁷ Yonsei University, Seoul, Korea

²⁸ CERN, Geneva, Switzerland

²⁹ The Henryk Niewodniczanski Institute of Nuclear Physics, Polish Academy of Sciences, Kraków, Poland

³⁰ Ruprecht Karls Universitaet Heidelberg, Heidelberg, Germany

³¹ INFN Sezione di Catania, Catania, Italy

³² Departamento de Física, Universidad de Oviedo, Avda. Oviedo, Spain

³³ Instituto Universitario de Ciencias y Tecnologías Espaciales de Asturias (ICTEA), Oviedo, Spain

³⁴ Universidad Nacional Autonoma, Mexico City, Mexico

³⁵ Karlsruhe Institute of Technology, Institut für Astroteilchenphysik, Karlsruhe, Germany

³⁶ Helmholtz Research Academy Hesse for FAIR, Goethe University Frankfurt, Frankfurt, Germany

³⁷ Charles University, Prague, Czech Republic

³⁸ INFN Sezione di Cagliari, Cagliari, Italy

³⁹ University of Bergen, Bergen, Norway

⁴⁰ University of Houston, Houston, USA

⁴¹ Utrecht University, Utrecht, The Netherlands

⁴² Massachusetts Institute of Technology, Cambridge, USA

⁴³ SUBATECH, Nantes University-IN2P3/CNRS-IMT Atlantique, Nantes, France

⁴⁴ Central China Normal University, Wuhan, China
**Abstract** This paper is a write-up of the ideas that were presented, developed and discussed at the fourth International Workshop on QCD Challenges from pp to AA, which took place in February 2023 in Padua, Italy. The goal of the workshop was to focus on some of the open questions in the field of high-energy heavy-ion physics and to stimulate the formulation of concrete suggestions for making progresses on both the experimental and theoretical sides. The paper gives a brief introduction to each topic and then summarizes the primary results.

# 1 Introduction

In hadronic collisions complex many-body systems of strongly-interacting particles are produced. The strong interaction, described in the Standard Model by quantum chromodynamics (QCD), determines the properties of the systems formed in these collisions and the evolution of the observed phenomena with the system complexity. In collisions of heavy nuclei (AA) at ultra-relativistic energies a system of deconfined quarks and gluons, the quark-gluon plasma (QGP), is formed. The existence of a QGP state is expected from QCD calculations on the lattice [1–4]. The search of evidences of QGP formation and the study of its properties have shaped the high-energy heavy-ion physics program at accelerators in the last 30 years [5–14] and traced the road for future experiments [15–18]. Most of the observables and probes used to investigate the QGP were first studied in detail in more elementary, smaller, systems, like electron–positron ($e^−e^+$), proton–proton (pp), and proton–nucleus (pA) collisions. These studies were fundamental for the comprehension of many aspects of the strong force and the development of QCD. On top of providing a baseline for interpreting the measurements in AA collisions, pp and pA collisions can be exploited to study the onset of phenomena ascribed to QGP in AA collisions. The goal of the workshop series “QCD challenges from pp to AA collisions” is to bring together experimental and theoretical physicists involved in the study of strong-interaction phenomena, mainly at high-energy particle colliders, with the aim of identifying the main open questions in the field and proposing new ideas and directions of research for addressing them. The focus is on those phenomena in hadronic collisions that are possibly sensitive to the size of the interacting system.

The 4th edition of the workshop was organized in Padua, Italy, from 13 February to 17 February 2023. It consisted in plenary and round-table sessions related to the following six discussion tracks:

1. Initial state and ultra-peripheral collisions

2. Jet production and properties in pp collisions and in the medium

3. Event properties and hydro in small and large systems

4. Hadronization of light and heavy flavour across collision systems

5. Energy loss and transport in the medium and in small systems

6. QCD and astrophysics

In plenary sessions invited speakers gave an overview on their specific field of interest. During round-table sessions, one per track, small groups (7–8 persons) discussed in detail new ideas and studies from both the experimental and theoretical sides. Moreover, special round-table sessions where groups from two or more tracks discussed together about possible synergies have been organized. The conveners of each track wrapped up the main ideas emerged and presented the outcome of the discussion in the plenary session on the last day of the workshop. In this document the outcome of these discussions is reported.

A section is dedicated to each track. The choice of the tracks was determined by considering the latest advances and open questions in the field, briefly outlined in what follows. In recent years, an intensive experimental campaign of measurements of jet and heavy-flavour (charm and beauty) production, azimuthal anisotropy, and correlations was carried out, giving an unprecedented insight into the partonic interactions in the medium and setting important constraints for the understanding of partonic energy loss and of heavy-quark transport. This campaign will continue in the next years, also exploiting detector upgrades, with novel measurements of jet substructure and with new and precise measurements of several charm and beauty hadron species, in particular of baryons. A prerequisite for the interpretation of production measurements will be the characterisation of the nuclear parton distribution functions (nPDF) and of initial-state effects. These topics are discussed in Sects. 2, 3, and 6.

In the last decade, phenomena like long-range flow-like angular correlations, studied in heavy-ion collisions and typically interpreted as a consequence of the formation and expansion of a medium have been observed, with different magnitude, in pp and pA collisions systems. A major challenge for the future is understanding the onset of these phenomena across collision systems, their dependence on event properties like particle multiplicity, the role played by multiple parton interactions (MPI), colour reconnection, and fluctuations in the initial state, as well as explaining the absence of evidences for other phenomena like partonic energy loss that naturally accompany flow observations in heavy-ion collisions. This aspect, which was discussed in most sessions, is one of the main subjects of Sects. 4 and 6. In recent years, the paradigm that heavy-quark hadronisation should proceed

<footnote>a e-mail: andrea.rossi@pd.infn.it (corresponding author)</footnote>

<footnote>b e-mail: lorenzo.sestini@pd.infn.it (corresponding author)</footnote>
similarly in $e^−e^+$ and pp collisions, which motivated the usage of a factorisation approach for calculating charm- and beauty-hadron production cross sections, has been severely questioned by the observation that charm and beauty baryon production relative to that of mesons is significantly larger in pp collisions than in $e^−e^+$ collisions. Tracing the modification of the hadronisation process in different hadronic environments, from pp to AA collisions, and the possible emergence of quark coalescence as a relevant hadronisation process already in small collision systems, is a major goal for the experiments. This topic is the subject of Sect. 5.

Many QCD-related measurements performed at hadronic colliders can provide important pieces of information to answer open questions related to cosmic-ray physics and astrophysical objects. The outcome of the discussions emanating from this theme is reported in Sect. 7. Section 8 closes the document with a summary and an outlook.

We remark that the choice of topics addressed in the following sections does not aspire to systematically review all open questions in the field. It is rather driven by the most important points presented by the participants at the workshop and the results of discussions and brainstorming.

## 2 Initial state and ultraperipheral collisions

The colliding particles in high-energy hadron collisions are dense gluonic systems. As a result, describing the initial state of high-energy hadron collisions requires understanding QCD at high gluon densities. These conditions can be studied by probing the structure of heavy nuclei at low $x$ and $Q^2$. The longitudinal structure of nuclei can be described using nPDFs [19–22]. In nPDF analyses, the nucleon-to-proton PDF ratio is parameterized as a function of $x$, and evolution in $Q^2$ is governed by the linear DGLAP equation. At low $x$ and $Q^2$, however, the spatial separation between gluons may be smaller than the size scale of the probing interaction. In this low-$x$ regime, nonlinearity may affect the evolution of parton densities, and coherent interactions with multiple gluons could become important. These so-called saturation effects are described by the color-glass condensate (CGC) effective field theory [23]. Experimentally, the low-$x$ regime is most commonly accessed at hadron colliders using two classes of observables: forward inclusive particle production in pA collisions and photoproduction in ultraperipheral collisions (UPCs). Observing and studying nonlinear QCD effects is one of the primary goals of initial state physics.

### 2.1 Forward measurements and low-$x$ suppression

Forward detectors at high-energy hadron colliders are ideal tools for studying the low-$x$ regime. The LHCb detector, for example, can probe $x$ below $10^{-5}$ in charm production

in pPb collisions. LHCb measurements of the $D^0$ nuclear modification factor in pPb collisions at $\sqrt{s_{NN}} = 5$ TeV [24] have been used in state-of-the-art nPDF fits to obtain precise descriptions of the gluon distribution at low $x$ [19, 20]. The improvement in precision over previous-generation nPDF fits is illustrated in Fig. 1. The LHCb $D^0$ data offer clear evidence of suppression of the low-$x$ gluon density in heavy nuclei. The resulting suppression is consistent with both collinearly factorizable shadowing of the gluon nPDF and gluon saturation in the CGC framework. As a result, it is now clear that observing nonlinear QCD effects in forward particle production will require high experimental precision across a wide range of observables.

In order to push to lower $x$ and $Q^2$, the LHCb collaboration has studied both charged hadron [27] and neutral pion [28] production in pPb collisions. The forward results again agree with both nPDF and CGC predictions. The ALICE collaboration has also measured the nuclear modification of charged hadrons and $\pi^0$ mesons in pPb collisions at central rapidity [29–31]. The ALICE pPb results show signs of suppression at $p_T \lesssim 2$ GeV, as opposed to the Cronin-like peak observed in PbPb collisions. The CMS collaboration has also measured the nuclear modification of charged hadrons [32], observing a similar suppression to that observed by ALICE at low $p_T$. Altogether, the LHCb, ALICE, and CMS results probe the nuclear gluon density from $x \lesssim 10^{-5}$ to $x \gtrsim 10^{-1}$ over a wide range of $Q^2$.

In addition to nonlinear effects in the evolution of parton densities, gluon saturation can also result in multiple scattering with low-$x$ gluons, leading to modifications of particle correlations at forward rapidity in pPb collisions. Both the PHENIX and STAR experiments have studied di-pion correlations at forward rapidity, observing a suppression of back-to-back low-$p_T$ di-pion pairs in deuteron-gold (dAu) and pA collisions [33, 34]. The STAR collaboration also observed that this suppression scales with $A^{1/3}$, consistent with saturation predictions. The ATLAS collaboration sees a similar suppression of back-to-back dijets in pPb collisions [35].

An ideal probe of the nuclear gluon distribution is direct photon production in pA collisions [36]. Direct photons probe the gluon nPDF at leading order and could provide access to the saturation regime at low $p_T$ and forward rapidity. The nuclear modification of low-$p_T$ direct photons could reveal nonlinear QCD effects, and photon-hadron correlations could provide evidence for gluon saturation. The ALICE collaboration is planning to install a high-granularity forward electromagnetic calorimeter (FOCAL) that will allow for measurements of forward direct photon production in Run 4 of the LHC [37]. The LHCb collaboration is also planning on installing tracking stations inside of its dipole magnet before Run 4, dramatically increasing its acceptance for low-momentum converted photons [38].
<figure><img src="image_2.png" /><figcaption>**Fig. 1** Comparisons of current- and previous-generation gluon nuclear modification factors $R_{pPb}^g$ in $^{208}$Pb from the (left) EPPS [20,25] and (right) nNNPDF [19,26] collaborations. The shaded regions show the</figcaption></figure>

68% CL uncertainties. In both cases, the fit results shown in gray do not use LHCb D⁰ data [24], while the results shown in color do use these data

## 2.2 Probing low $x$ with vector-meson photoproduction in UPCs

Studies of photoproduction of vector mesons in UPCs have provided complementary probes of the gluonic structure of nucleons. In UPCs, the electric field of a nucleus produces quasi-real photons, which can interact coherently with the colliding nucleus. The photon flux scales with $Z^2$, where Z is the atomic number of the nucleus. Furthermore, hadronic interactions are suppressed in UPCs due to the large impact parameter of the collision, which exceeds the sum of the radii of the colliding nuclei. This results in an experimentally clean signal. Vector meson photoproduction is sensitive to the gluon density of the nucleus at leading order at $x = me^{\pm y}/\sqrt{s_{NN}}$, where m is the vector meson's mass and y is its rapidity. As a result, vector meson photoproduction in high-energy UPCs can probe gluon densities at extremely low $x$ [39,40].

Coherent vector meson photoproduction has been recently measured in heavy-ion collisions for $\rho^0$ [41–43], $J/\psi$ [44–47], and $\psi(2S)$ [44,46]. The results cover a wide range in rapidity and show partial agreement with theoretical models. The heavy-ion cross section can be used to calculate photonuclear cross section. However, this calculation suffers from an ambiguity because the photon could be emitted by either nucleus, resulting in two possible $x$ values for each data point. Two solutions have been proposed to resolve this ambiguity. First, vector meson photoproduction can be measured in multiple nuclear breakup classes, which can then be used to extract photoproduction cross sections at both possible values of $x$ [48]. This method has been used by both the ALICE and CMS collaborations to extract low-$x$ photoproduction cross sections [49,50], which are shown in Fig. 2 (left). Second, coherent production in peripheral and ultraperipheral collisions can be combined [51]. Coherent production in peripheral heavy-ion collisions has been measured

by both the ALICE and LHCb collaborations. Both methods allow for access of unique $x$ range down to $x \sim 10^{-5}-10^{-6}$. The experimental results show strong gluon suppression relative to the Impulse Approximation [52], which can be best described by models with either shadowing or saturation effects [53]. Additionally, the $|t|$-dependence of the photonuclear cross section, where $|t| \simeq p_T^2$, is directly sensitive to the spatial distribution of gluons in the nucleus [54–56].

The gluon distribution in the proton can be studied using vector meson photoproduction in pp or pPb collisions. The results of LHCb [58,59] and ALICE [60–62] extend the exclusive $J/\psi$ Bjorken-$x$ photonuclear cross section coverage down to $\sim 10^{-6}$. However, no clear indication of gluon saturation at low $x$ is visible, and data follow trends observed in ep collisions at HERA [63–65]. A similar trend is visible in $Y(nS)$ [66,67] and $\rho^0$ [68] meson data from the CMS collaboration. Deviation from this trend might be visible in the dissociative vector meson photoproduction cross section at higher collision energies reachable during Runs 3 and 4 of the LHC.

## 2.3 Progress in higher-order calculations for $J/\psi$ photoproduction

The theoretical interpretation of vector meson photoproduction data from UPCs is limited by large factorization scale uncertainties [69–72]. Furthermore, the LO and NLO results have opposite signs. This contrasts with the expectations of a forward $t=0$ elastic scattering amplitude based on Regge theory, and in general hints at poor perturbative stability of the amplitude. This, together with the sensitivity of the process to generalized parton distributions, poses a major obstacle to including exclusive heavy vector meson photoproduction data in global PDF analyses.

Recent theoretical work has aimed to overcome these challenges. PDFs can be related to GPDs using the so-called
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /><figcaption>**Fig. 2** J/ψ photoproduction cross sections for γPb (left) from ALICE [49] and CMS [50] measured using nuclear breakup categories to unambiguously determine the energy dependence and for γp (right) from HERA and LHCb measurements versus NLO predictions [57]</figcaption></figure>

Shuvaev integral transform at moderate-to-low values of the skewness parameter $\xi \sim x/2$, which quantifies the fraction of longitudinal momentum transfer from the initial state hadron to that in the final state [73–75]. The conventional collinear factorization coefficient functions for exclusive heavy vector meson photoproduction at NLO can then be supplemented with corrections of $O(Q_0^2/\mu^2)$, where $Q_0$ is the GPD input scale and $\mu$ the factorization scale [76]. For exclusive J/ψ photoproduction, $\mu = O(m_c)$, where $m_c$ is the charm quark mass, so such corrections are of $O(1)$ and are crucial ingredients in the description of the process. Moreover, the choice $\mu = m_c$ allows for the resummation of a particular class of large double logarithms $(\alpha_s \ln(1/\xi) \ln(m_c^2/\mu^2))^n$ that arise in the high-energy limit of the amplitude [77]. Collectively, the implementation of this low cut-off procedure and program of effective low-$x$ resummation provide for a more reliable and stable $\overline{\text{MS}}$ amplitude to NLO with a reduced factorization scale dependence, allowing for a sensible comparison to the experimental data [57].

The central cross section predictions based on input GPDs constructed from three global PDF analyses via the Shuvaev transform differ dramatically in the kinematic region accessible to the LHC but are compatible at HERA energies, see Fig. 2 (right). This observation, together with the relative size of the input PDF uncertainties in the low-$x$ and low-$\mu$ domain compared to the experimental uncertainties and the systematically tamed scale uncertainty discussed above, motivated an extraction of a low-$x \sim 10^{-3} - 3 \times 10^{-6}$ and low-$\mu^2 \sim 2.4 \text{ GeV}^2$ gluon PDF over both HERA and LHC kinematic regions via independent fitting and statistical reweighting approaches using the exclusive J/ψ photoproduction data gathered from HERA and the LHC [78]. More global impact studies of this data within a larger fitting framework would be desirable to allow for refined extractions of PDFs and provide novel constraints at low $x$ and low

$\mu$, which would be of practical use to the PDF fitting groups and to the wider particle physics community.

## 2.4 Challenges and opportunities

Substantial experimental and theoretical progress has been made towards developing a coherent picture of the gluonic structure of nuclei from forward particle production and vector meson photoproduction measurements. However, the proton PDF extractions using HERA and LHCb J/ψ photoproduction data discussed above differ dramatically from proton PDF determinations using LHCb $D^0$ meson production data [57]. This discrepancy shows that a great deal of work is needed to reconcile these two classes of observables through better theoretical understanding of the processes e.g. in terms of low-$x$ resummations and applied subtraction schemes. This reconciliation is becoming increasingly urgent as nPDFs achieve new levels of precision yet continue to produce predictions compatible with CGC calculations, demonstrating the difficulty of conclusively observing the effects of non-linear QCD. Vital new experimental inputs are expected in the near future as the LHC prepares to produce high-energy OO and pO collisions. Vector meson photoproduction data from OO UPCs could be used to further tame factorization scale uncertainties by constructing ratios with PbPb measurements [70]. Furthermore, forward particle production measurements in pO collisions will probe the A-dependence of gluon suppression in nuclei, potentially revealing the onset of saturation effects at low $x$.

# 3 Jet production and properties in pp and in the medium

The fragmentation of energetic quarks and gluons in relativistic heavy-ion collisions is modified with respect to its vacuum
counterpart due to the presence of a QGP. Hard-propagating partons scatter with the constituents of this color-deconfined medium both elastically and inelastically (see e.g. Ref. [79] for a recent review). Inelastic scattering processes trigger additional radiation governed by an emission probability distinct from the standard Altarelli-Parisi (AP) splitting kernels. One of the peculiarities of the medium-induced spectrum is the absence of a collinear singularity. Consequently, medium-induced emissions typically appear at large angles, sometimes even outside the cone of the reconstructed jet, leading to a degradation of the jet energy. This suppression of the jet $p_T$ spectrum in AA collisions with respect to pp collisions has been experimentally confirmed both at RHIC and LHC energies. An in-depth understanding of energy loss mechanisms requires an experimental program that goes significantly beyond inclusive observables.

During the workshop, we explored the potential of jet substructure measurements to pin down one of the fundamental scales of in-medium jet physics, namely the characteristic angle of medium-induced emissions, $\theta_c$. This angular scale also controls the interference pattern of multiple emitters [80] and, as a consequence, the amount of energy lost. An unambiguous determination of $\theta_c$ would not only inform about these fundamental properties of jet evolution in a dense medium but also about properties of the QGP itself, such as the effective length $L$ or transport properties often characterized by the coefficient $\hat{q}$, the mean squared $p_T$ transferred per unit length. In this report, we first review the current status of jet substructure measurements that target the determination of $\theta_c$. Next, we discuss the potential of new observables proposed in the talks by the participants and conclude with a brief overview of the follow-up discussions.

## 3.1 Experimental searches of color coherence effects

The quest for color coherence effects in jet observables has been an active field of research during the last decade. Here, we focus on the most recent developments and refer the reader to Ref. [84] for a complete review. One of the most recent measurements aiming at unveiling color coherence dynamics is that of the groomed jet radius, $\theta_g$, by the ALICE experiment [81] (Fig. 3, left). The distribution was measured to be narrower in PbPb than in pp, i.e. splittings with $\theta_g < 0.2$ are enhanced in the medium. The data points agree with a wide set of theoretical predictions based on quite different ingredients, some of them lacking a color coherence angle implementation.

The ATLAS collaboration extended the measurement of $\theta_g$ by studying the nuclear modification factor $R_{AA}$ as a function of $\theta_g$ [82] (Fig. 3, middle). For pp collisions, this double differential measurement revealed sizeable discrepancies (up to 50%) between general purpose Monte Carlo event generators and data, thus highlighting the importance of measur-

ing the pp baseline to guarantee a meaningful interpretation of heavy-ion-to-pp ratios. Regarding the PbPb data, ATLAS observes an increasing jet suppression with increasing $\theta_g$, with a possible inflection point at around $\theta_g = 0.05$. This result is in qualitative agreement with a very recent measurement [83] that does not rely on SoftDrop grooming but rather utilizes a re-clustering procedure in the spirit of trimming [85] in which hard subjets are found inside large-$R$ jets. In that measurement (Fig. 3, right), $R_{AA}$ was evaluated for configurations with different subjet multiplicity, and it was found that large-$R$ jets with single subjets are quenched significantly less than jets with a complex topology (a factor of 2 difference in terms of the jet $R_{AA}$). In addition, $R_{AA}$ remains flat as a function of the opening angle between subjets ($\Delta R_{12}$) for subjets with $\Delta R_{12} > 0.2$ which may be interpreted as a constraint on the maximal value of $\theta_c$. In fact, the value of $\theta_c < 0.2$ is consistent with previous measurements of radius ($r$) dependent fragmentation functions [86] where only minimal changes of the jet fragmentation pattern in PbPb collisions were observed for $r < 0.1$.

The aforementioned measurements point to a narrowing of the angular scale of jets in PbPb relative to pp collision. A possible caveat of such inclusive jet measurements is that the comparison is made at the same reconstructed jet energy, which does not map to the same scattered parton energy in pp and PbPb due to the additional energy loss that is present in the latter case. This can potentially lead to selection biases [87,88].

The theoretical interpretation of these results is not yet unambiguous and it is therefore important to discuss novel observables with $\theta_c$-sensitivity that can help disentangle between different effects.

## 3.2 Novel observables with $\theta_c$-sensitivity

The observables discussed during the meeting can be categorized into two classes: clustering-tree based and energy-flow based. In the former category two proposals rely on the dynamical grooming (DyG) method [89]. On the one hand, the dynamically groomed jet radius (akin to $\theta_g$ but probing a wider region of phase-space) was proposed as a potential candidate to detect color coherence effects [90]. This observable has been already measured in pp collisions [91] and benchmarked against a theoretical prediction at high-logarithmic accuracy in perturbative QCD [92]. In the medium case, an analytic calculation was also presented together with a dedicated MC study of non-perturbative effects such as hadronisation and medium response. The key difference between the SoftDrop and Dynamical Grooming versions of $\theta_g$ is that while the former leads to a shift in the distribution towards smaller values in PbPb, the latter predicts a genuine change in the shape of the distribution, i.e. the vacuum distribution is almost flat (for some DyG settings) and color coherence
<figure><img src="image_5.png" /><figcaption>**Fig. 3** A compilation of recent jet substructure measurements indicating a narrowing of the jet core in heavy-ion collisions. From left-to-right the figures have been extracted from Refs. [81–83]</figcaption></figure>

ence dynamics induces a peak around $\theta_c$. The discriminating power of this observable is therefore potentially enhanced.

Vacuum and medium-modified showers factorize in time to first approximation [93], with the harder vacuum part happening earlier. A broader early vacuum shower will undergo a more active medium-induced shower due to the larger number of emitters. This introduces a jet-width dependence of energy loss, which is naturally interwound with flavour dependence. These effects coexist with color coherence, which further regulates the effective number of emitters in the medium.

The color charge dependence of inclusive jet suppression was explicitly studied in several phenomenological works, see e.g. Refs. [94,95]. An experimental strategy to distinguish between these confounding effects was introduced in Ref. [96] which proposes to exploit the forward capabilities of LHC detectors.

Finally, a complementary set of observables that do not rely on clustering sequences are the so-called energy-energy correlators [97] (EECs) that describe n-point correlation functions of energy flux. In recent years, there has been a renewed interest on EECs both theoretically and experimentally. On the formal side, the logarithmic structure of these observables in vacuum entails a beautiful connection to light-ray operators and conformal field theory [98–100]. In a heavy-ion context, first steps towards calculating the two-point energy-energy correlator were presented in Refs. [101,102]. A leading-order calculation demonstrated an enhancement of this observable at large angles attributed to color decoherence dynamics. The authors also proposed an experimental measurement of this observable in $\gamma$+jet events to mitigate the selection bias effect. The quantitative impact of energy loss, underlying event and medium response on this observable, although formally power-suppressed by construction of the observable, remains to be studied. It would also be beneficial to perform an in-depth study of the com-

plementarity of clustering-tree based and energy-flow based observables to study jet quenching.

### 3.3 Future steps

The multiple scales involved in the evolution of jets in heavy-ion collisions makes any attempt of isolating $\theta_c$-dynamics an extremely challenging task. We identified a few basic requirements that any observable needs to meet in order to qualify for such a task: (i) calculable vacuum baseline, (ii) resilience to underlying event, and (iii) reduced sensitivity to medium response. We believe that jet substructure observables computed in high-$p_T$ jets could potentially meet these requirements since a sufficiently large separation of scales is guaranteed.

## 4 Event properties and hydro in small and large systems

The development of collective flow behavior [103–106] and the enhancement of strange hadron production [107,108] have been considered as signatures of the presence of the QGP in nuclear collisions. Such phenomena were considered exclusive to heavy-ion collisions until similar features were also observed in high-multiplicity pp and pA collisions. This raised the question of whether they are due to the formation of a strongly interacting QGP medium or could be explained by other mechanisms, e.g., initial-state momentum correlations. The former scenario is supported by the success of hydrodynamical frameworks in describing the anisotropic flow patterns seen in small collision systems [109]. However, collective-like effects have been observed also in ultraperipheral AA collisions [110] and recent ALICE results suggest that a non zero elliptic flow is measured even in low-multiplicity pp collisions. Therefore, small systems have grown into prominence in the study of the hot QCD medium and are no longer regarded solely as reference measurements for the studies
performed with heavy-ion collisions. This discussion track focused on the frameworks describing event properties and hydro in small and large systems.

## 4.1 Global description across systems: core-corona approaches

The open questions regarding the interpretation of QGP-like effects in small colliding systems require further theoretical development of the description of the medium created in relativistic nuclear collisions. In particular, there is a need for phenomenological approaches that take into account non-equilibrium effects that become more and more important going from large to small systems and that also include the description of both a ‘core’ (hydrodynamically evolving QGP) and a ‘corona’ (non-thermalized particles at low densities or high momentum) part of the produced medium [111].

Indeed, the relativistic hydrodynamic approaches describe the matter in local thermal equilibrium (ideal hydro) or with moderate deviations around it (viscous and anisotropic hydro). Still, other parts of the system, such as matter far out-of-equilibrium and propagating jets, are often not included in the hydrodynamic description. This limits the applicability of the hydro framework for small systems and at lower collision energies. The core-corona approaches try to solve this issue, allowing simulations of large and small systems, as well as a wide range of collision energies: from the highest LHC energies down to the lowest energies of the beam energy scan at RHIC. It is worth mentioning that this capability is also intrinsic in kinetic approaches, in which both the partonic and the hadronic phases are described through the Boltzmann equations (e.g., AMPT [112]) or generalized transport equations (e.g., PHSD [113]). Here the dynamical separation of a core and corona part is automatically achieved during the space-time evolution of the fireball created in relativistic nuclear collisions. In this sense, such approaches can be considered as belonging to the broader class of core-corona approaches.

Two examples of core-corona models with hydrodynamic formulation for the QGP phase are DCCI2 [114,115], EPOS2/3 [116,117] and EPOS4 [118,119]. The latter has been recently publicly released [120]. The experimental data on charged-particle multiplicity and $\Omega/\pi$ ratio from pp to AA collisions are used both in the DCCI2 model and EPOS to fix the model parameters.

As seen from Fig. 4 (left) the fractions of the core and corona parts exhibit a clear scaling with multiplicity. There is a smooth increase of the core part and a decrease of the corona part going from low to high multiplicity. The onset of core dominance within the DCCI2 approach happens at $\langle dN_{ch}/d\eta \rangle \sim 20$. The same scaling is present in observables sensitive to particle multiplicity, such as the $\Omega/\pi$ ratio, see Fig. 4(middle). EPOS can also describe the trend of the

experimental measurements for the mean transverse momentum $\langle p_T \rangle$ as a function of multiplicity shown in Fig. 4 (right), including the jump seen when passing from pp to PbPb collisions. Therefore, core-corona approaches are able to describe the hadrochemistry and kinematics across systems. It is interesting to note that the core part has an important role even in pp collisions.

Furthermore, an important feature of core-corona models is the ability to describe multiple observables across many collision systems simultaneously. This feature is crucial for an over-arching understanding of how similar phenomena can appear from small to large collision systems.

## 4.2 Deeper insights into nuclear collisions: realistic 3D view

While the majority of studies are related to observables at midrapidity, in the last decade, it has become more and more important to gain a realistic three-dimensional view of the entire interaction area of the collision. For example, describing the tilted fireball produced in non-central heavy-ion reactions [122] as well as the asymmetric particle rapidity distributions generated in pA and other asymmetric collisions [123,124]. In relativistic AA collisions, the tilt of the source is converted by the collective expansion of the produced hot matter into the directed flow $v_1$, which is the first harmonic of the particle azimuthal distribution in momentum space. The directed flow of light hadrons constitutes a probe of the asymmetric geometry of the matter distribution in the longitudinal direction in pA collisions [124]. Full 3D simulations are also crucial for the comparison to the elliptic and triangular flow data in small systems. In particular, in [123] it is estimated that approximately 50% of the difference between triangular-flow values measured by PHENIX and STAR originates from the different flow correlations between different rapidity regions. This highlights the importance of the correct 3D view of the fireball and requires theoretical approaches to have the initial condition modeled in ‘realistic’ ways and support further investigation of collectivity and correlations at large rapidity.

Surprisingly, it turned out recently that the $v_1$ of heavy-flavor particles is much larger than that of light charged hadrons in both large and small systems [125,126]; the origin of the $v_1$ of D mesons is different from the one of the bulk medium and its large value manifest only due to the non-perturbative interaction of heavy quarks with the QGP medium [127]. Therefore, the directed flow allows 3D access to the production of both soft and hard probes and a deeper look into their mutual interaction.

LHCb collaboration is a relative newcomer to collective flow measurements in heavy-ion collisions. However, thanks to upgrades in Run 3 LHCb will be able to study PbPb collisions up to 30% centrality. The large forward rapidity reach
<figure><img src="image_6.png" /><figcaption>**Fig. 4** Left: fraction of core and corona as a function of multiplicity from pp to PbPb. Middle and right: integrated $\Omega/\pi$ yield ratio (middle) and average transverse momentum (right) as a function of charged-particle multiplicity density at midrapidity from pp to PbPb compared to EPOS 4.0.0 predictions from various parts of the collision system [121].</figcaption></figure>

The “co+co” (core+corona) curve is an interpolation between the core and corona cases, with the core weight increasing continuously with multiplicity. The “full” case is equal to the co+co case with in addition the inclusion of hadron rescattering

of LHCb puts it in a unique position to contribute to the study of the full 3D picture of nuclear collisions. In parallel, LHCb will take data in a fixed target mode with SMOG2 for a wide range of nuclear targets at a smaller center of mass energy but with high luminosity, allowing to probe a variety of initial-state geometries.

## 4.3 The evergreen tale of collectivity

Multi-particle correlation measurements are arguably the most versatile and precise tool in studying collective phenomena in nuclear collisions. The centrality and momentum dependence of elliptic and triangular flow coefficients are the key inputs of the state-of-the-art Bayesian inference analyses of QGP properties like specific shear and bulk viscosities [128–130]. The momentum dependence of baryon and meson elliptic flow is used to benchmark different hadronization mechanisms (e.g. fragmentation or coalescence) at different momentum scales. More generally the observed similar mass ordering of collective flow in PbPb and pPb collisions points to the common origin of this effect.

The crucial outstanding question is whether the initial state effects, e.g. initial momentum correlations, could be resolved in flow measurements in small systems. One very promising observable is the correlation between fluctuations in mean $p_T$ and elliptic flow magnitude [131,132]. In central PbPb collisions, the elliptic and radial flows have a positive correlation. When inspecting this correlation as a function of multiplicity, it is important to consider that the correlation between particle multiplicity and impact parameter becomes weaker going towards peripheral collisions, and it is even more diluted in small collision systems (see e.g. [124]). For a given multiplicity interval, at a larger impact parameter (smaller overlap area) the initial energy density must

<figure><img src="image_7.png" /><figcaption>**Fig. 5** Pearson correlation coefficient $\rho(c_2\{2\}, [p_T])$ as a function of multiplicity as measured by ATLAS [134] and CMS [135]</figcaption></figure>

increase to maintain the same multiplicity. In very peripheral collisions, the collision geometry is no longer controlled by the impact parameter, but rather by the large fluctuations of the number of participating nucleons relative to the average. With the net entropy results, more deformed shapes have a lower density than more compact configurations leading to the anti-correlation of elliptic flow and mean $p_T$ fluctuations. In contrast, the initial state models predict positive correlations. Indeed the Pearson correlation coefficient measured by ATLAS and CMS collaborations shows the characteristic change in behavior at low and high multiplicities, as shown in Fig. 5. However, a firm conclusion about measuring the initial state effects is complicated by the measurement's sensitivity to experimental cuts and similar behavior seen in models without initial state correlations, e.g. AMPT [133]. Further studies of this fascinating behavior will be pursued.
The difficulty to link the collective behavior seen in small systems to initial-state effects or to identify it as a QGP signature has stimulated first attempts to study observables related to particle production and collectivity by means of multidifferential analysis, in which centrality selection is flanked by event-shape categorization of events. An example of such event-shape selection methods employed from small to large systems is the transverse spherocity analysis [136–140], which has been used to classify events based on different degree of the collective effects in heavy-ion collisions [139]. Other event-shape selection techniques are those based on the relative transverse activity classifier [141, 142] and on flattenicity [143, 144], mainly employed in pp collisions, as well as the event-shape engineering based on flow vector [145], until now exploited especially in heavy-ion experiments at LHC for studying correlations between flow harmonics of different order [146]. Various studies are exploring the potentiality of these event-shape classifiers. For instance, the capability of flattenicity to identify collisions with multiple lower transverse-momentum parton-parton scatterings and softer $p_T$ spectrum compared to the multiplicity estimator may be exploited during LHC Run 3 and 4 to isolate high-multiplicity pp collisions originated by soft partonic processes. Moreover, flattenicity is a promising variable to overcome the drawback of multiplicity-based event classifiers that, suffering from biases towards multijet final states, make challenging the quest for jet-quenching in high-multiplicity pp collisions [144]. Pursuing further studies based on multidifferential measurements in pp and in pA collisions may help to gain deeper insights into the origin of collectivity in small colliding systems.

The collective effects of QGP can manifest themselves not only in the kinematic distribution of final state particles but also leave imprints in their polarization. The short, but rapid rotation of the medium achieves very large values of vorticity. In a thermal system, this leads to spin alignment of quarks and consequently the spin of produced hadrons. By measuring the spin of produced hadrons with respect to the collision axis, the amount and direction of hadron polarization can be estimated. Polarization physics has attracted a lot of attention from both experimental and theoretical communities [147]. The Λ polarization puzzle and recently measured vector meson polarization have led to new theoretical developments. Various phenomena can be studied with polarization, ranging from global features due to the collectively expanding medium to more localized vorticity that appears due to high-energy particles propagating through the medium [148, 149]. Notably, vorticity - and subsequently also polarization - has been shown to encode information about medium properties like viscosity, potentially rendering polarization measurements a fundamental tool in characterizing the QGP.

## 5 Hadronization of light and heavy flavour across collision systems

By now, there is abundant evidence that many non-perturbative “constants” such as baryon-to-meson ratios, strangeness fractions, and (moments of) hadron kinematic distributions, can depend sensitively on the type of colliding beam particles (the collision system) and on event activity (the production environment).

Using rapidity density of charged particles as a proxy for event activity, many (though not all) salient observables exhibit continuous increases with progressively higher average multiplicity densities in pp, pA, and AA collisions, typically reaching approximately constant asymptotic values for (central) AA collisions.

Several physical models have been proposed that are capable of modelling these trends qualitatively, and increasingly even quantitatively. During the workshop the following types of models were discussed:

* General-purpose Monte-Carlo event generators (such as PYTHIA [150] and HERWIG [151]) mainly focus on $e^+e^-$, $e^\pm p$, and pp collision systems. Colour reconnections (CR) can increase baryon-to-meson ratios [152, 153] and non-perturbative colour interactions (CI) between strings or clusters can increase both baryon and strangeness fractions [154–157]. Both types of effects (CR and CI) depend (implicitly or explicitly) on the effective local density of coloured partons. Since there is no “medium” in these models, the main physical driver for the local density of coloured partons in hadron-hadron collisions is multi-parton interactions (MPI). These are absent in $e^+e^-$ and $e^\pm p$ collisions. Note however that $e^+e^- \to WW \to$ hadrons furnishes a case similar to that of double-parton interactions, and was instrumental to demonstrating the existence of CR at LEP 2 [158], albeit with relatively low statistics. Both mechanisms (CR and CI) can also produce flow-type effects [159–162].

* Models of quark coalescence are based on the idea that in the dense medium of quarks and gluons that is formed in high-energy AA collisions, but even in high multiplicity pA and pp ones, hadrons come from the coalescence of quarks with a similar velocity. Nearly twenty years ago it was realized that such an hadronization mechanism would be quite efficient in creating baryons, resulting in larger baryon yields w.r.t. those expected from standard fragmentation functions. Most of the coalescence approaches evaluate the hadrons yields as a convolution of the quark thermal (at least at low $p_T$) distribution functions at a temperature $T_H \simeq 155$ MeV and the hadron wave function [163, 164]. A statistical homogeneous color distribution is assumed in evaluating the
probability of quark coalescence. An alternative formulation of the idea of recombination instead evaluates the hadron formation as a convolution between the quark distribution function with a Breit–Wigner resonant scattering cross section instead of the hadron Wigner wave function [165]. Despite differences in the details, all these approaches show the general feature of an enhancement of baryon to meson ratio at intermediate $p_T$ with $p/\pi$, $\Lambda/K_S^0 \simeq 1$ [164]; furthermore they all agree that at $p_T > 6-8$ GeV the dominant process will be independent string fragmentation. Coalescence models share also the feature of predicting a quark-number scaling of the hadron elliptic flow, which has been also observed experimentally [163,164,166]. In the last decade these approaches have been extended to the charm sector and have successfully predicted large values of the $\Lambda_c^+/D^0$ ratio [167,168] with peak reaching about unity in AA collisions and about 0.5 in pp and pA collisions [169].

* Statistical hadronization models (SHM) are successful in describing hadron yields in (central) heavy-ion collisions over a broad range of collision energies [170]. The evolution of hadron ratios with multiplicity in pp, pPb and collisions is also well described [171], here with some assumptions on the canonical conservation volume for strangeness and baryon numbers [172].

* Core-corona models are designed to simulate relativistic pp, pA, and AA collisions across a wide range of energies. These models aim to provide a unified approach by considering collective effects in all systems and assuming the creation of a medium. In the core-corona procedure, high-multiplicity events are made by a dense core (bulk-matter) that thermalizes and expands collectively, and a corona near the surface. The separation between core and corona is considered a dynamical process based on density. The two main models considered are EPOS and DCCI and are explained in detail in Sect. 4.1.

* Another approach to in-medium hadronization, also based on the idea of a quark recombination, is the one developed on the background of a POWHEG+PYTHIA Monte Carlo event generators [173]. On a hadronization surface identified by the local temperature $T_H = 155$ MeV, a charm quark is recombined with a thermal quark or di-quark at finite mass values typical of those of the constituent quark picture. The recombination occurs if the invariant mass of the pair, $M_C$, is larger than the mass of the lightest charmed hadron. Then the cluster decays into a charmed hadron that has the same baryon number and strangeness of the cluster or, for large invariant masses, $M_C > M_{max} \simeq 4$ GeV, it hadronizes via string-fragmentation according to PYTHIA 6.4. A specific trait of this approach is to assume a thermal population of di-quarks that acquire also a hydrodynamical

radial and elliptic flows. In Sect. 5.2, this feature will be discussed in relation to the elliptic flow of $\Lambda_c^+$.

In Sects. 5.1 and 5.2 respectively, the current status in terms of available experimental measurements and their description by theoretical models is presented, focusing on light and heavy flavours, respectively. In Sect. 5.3, the current status of observations and physical models of exotic hadrons (carrying charm quarks) is discussed. The structure of these hadrons, whether tetraquarks or pentaquarks, is interesting and currently theoretically debated. These hadrons can also shed light on hadronization mechanisms, as suggested by the recent observation by LHCb [174] of an event-multiplicity dependence of the $\chi_{c1}$ (3872) hadron production in pp collisions. Finally, in Sect. 5.4, a summary is followed by a menu for next steps.

## 5.1 Light flavours

One of the most intensely studied observable is the multiplicity dependence of the yield ratio to pions for various strange particle species. Measured by ALICE [108] across pp, pA, and AA collisions, the ratio is observed to exhibit a strong multiplicity dependence for pp and pA collisions, progressively stronger for multi-strange hyperons.

It was discussed during the workshop that it would be advisable to indicate also the values measured in $e^+e^-$ collisions when showing these plots, for a clearer “vacuum” reference value. Interesting attempts were also made to separate this into quark and gluon jets [175] in $e^+e^-$ collisions: the interpretation of the data is however complex also because of the dependence of the results on the details of the method used to separate the classes. Non-trivial multiplicity dependence of these quantities in $e^+e^-$ collisions is also possible [176], though so far little explored.

An important conundrum that any physical model of these effects will need to address is the fact that the proton-to-pion ratio has been observed to exhibit very little (if any) dependence on multiplicity, combined with the intriguing fact that the values of this ratio observed in LHC collisions appears to be a bit *below* the value observed at LEP, see Fig. 6 (left). According to the PDG table of average identified-particle multiplicities in $e^+e^-$ collisions [179], $\langle N_p \rangle / \langle N_\pi \rangle$ at LEP is about 0.062±0.002, whereas the values measured by ALICE in pp collisions [108] range from about 0.048 ± 0.006 for the lowest charged-particle densities to about 0.055 ± 0.005 for the highest ones. This presents a challenge for many of the current dynamical models on the market, for which the LEP value typically acts as an effective lower bound. Even in models (or tunes) that do not use the LEP value as a direct constraint, one would still have to address at least in principle what physical mechanism accounts for the uni-
<figure><img src="image_8.png" /><figcaption>**Fig. 6** Left: Baryon-to-meson ratios measured by ALICE [108], with approximate indications of the values measured at LEP (DELPHI [177]) superimposed. Right: proton-Λ correlations in azimuth [178]</figcaption></figure>

versality breakdown between $e^+e^-$ and pp collisions. One potential mechanism that has been hypothesised to possibly reduce the number of proton-antiproton pairs is hadronic reannihilation [180] though it may be doubtful if that could be sufficiently active even at low multiplicities, as required by data. Another point that was made during the workshop is that it may be useful as a cross check to study how the numerator and denominator, i.e. the proton and pion yields, evolve separately, although the denominator is of course largely degenerate with the charged multiplicity itself. Other potential cross checks could include measurements of this (and other) ratios in high-$p_\perp$ (quark and gluon) jets. Hadrochemistry in jets is since recently being investigated with data already existing on hyperons [181] and on deuteron [182]. The former were important to demonstrate and quantify the different production rates of strange baryons relative to $K_S^0$ mesons in and out-of jets. The measured ratios also set constraints to diquark formation in jets to models like PYTHIA8, which cannot reproduce the data. The deuteron measurement instead supports the formation of deuteron in-jets via hadronic interactions occurring between hadrons close in phase-space data. Within uncertainties, the data can be reproduced both by modelling these interactions as a simplistic coalescence process and by a more realistic, four-momentum conserving, hadron rescattering model introduced in PYTHIA8.3.

A further interesting observation in baryon production in pp collisions is that there is a dip in both baryon-baryon and baryon-antibaryon correlations near $\Delta\phi \sim 0$ [178], which is not described by MC models, Fig. 6 (right), see also Ref. [183]. The suggestion was made that this might be interesting to follow up with similar measurements involving heavy-flavour baryons, though admittedly experimentally very challenging.

## 5.2 Heavy flavours

A rather surprising observation at the LHC is the large fragmentation fraction for charmed baryons in pp collisions [184], which is significantly larger than that measured at LEP. PYTHIA (CR Mode 2) [152] describes $\Lambda_c^+$ data but underpredicts $\Xi_c^{0,+}$ data. In a statistical hadronization approach, He and Rapp [185] assumed, under the guidance of the Relativistic Quark Model [186], the existence of a vast set of excited charm-baryon states, most of which yet unobserved [179], in order to describe the measured $\Lambda_c^+/D^0$ ratio. This approach underpredicts $\Xi_c^{0,+}$. SHM was recently applied to bottom hadronization in pp collisions [187].

The current status of the comparison of data and models is shown in Fig. 7, where the ratios of $D_S^+$ mesons and $\Lambda_c^+$ baryons to $D^0$ mesons are shown as a function of $p_T$ for pp and central PbPb collisions. The data exhibit a clear difference between pp and PbPb, which the models capture well, though overall the model description of the data needs further improvements.

The rather complete charm chemistry was recently measured by ALICE in central PbPb collisions [14], see Fig. 8 (left). The SHM describes (predicted [188]) the data very well, except for the $\Lambda_c$ baryon, where a good description is achieved only by the inclusion of many (tripled in number) excited charm-baryon states compared to the default (PDG) hadron spectrum. It is important to recall that SHM assumes a concurrent hadronization of all flavors in QGP at the crossover temperature of 156.5 MeV. A similarly quantitative description of charmed hadron yields in dynamical models remains challenging, but developments are under way. Measurements of flow for different charm-hadron species can provide an important test for the modelling of hadronization. Figure 8 (right) shows the $p_T$ dependence of the differ-
<figure><img src="image_9.png" /><figcaption>**Fig. 7** Charm hadron ratios vs. $p_T$, data in comparison to models [14] (Catania: coalescence; TAMU, GSI-Hd: SHM, TAMU with enhanced charm-baryon spectrum, GSI-Hd with the PDG one)</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Fig. 8** Left: charm hadron yields in central PbPb collisions, data [14] and SHM predictions. Note the scale factor of 30 for the $J/\psi$ meson. Right: The difference between the $v_2$ values of $\Lambda_c^+$ and $D^0$ as a function of $p_T$ as predicted by the Catania coalescence model and POWLANG</figcaption></figure>

ence between $v_2(\Lambda_c)$ and $v_2(D^0)$ predicted in a coalescence approach and in the POWLANG mechanism. In the coalescence model the $v_2$ of $\Lambda_c$ receives a contribution from two light quarks while the D meson only from one. In the POWLANG approach instead $\Lambda_c^+$ combines to a diquark evolving hydrodynamically with the bulk medium and has a milder difference w.r.t. D in particular at intermediate $p_T$.

In closing this section it is worth mentioning the possible presence of a persisting problem in the data on $\Lambda_c^+$ baryon, where the rapidity distribution with LHCb (forward rapidity) and ALICE (midrapidity) data appears too peaked at midrapidity. It is also worth recalling that hadronization in the charm sector is not only interesting per se, but crucial for the attempts to extract the heavy quark transport coefficients via (transport) model description of the data. Currently, the

spread in the model predictions due to hadronization is large [189].

## 5.3 Exotic hadrons

Exotic hadrons are those hadrons which do not fit in the standard meson and baryon classification [190]. We discuss briefly the charm-carrying tetraquark candidates $\chi_{c1}(3872)$, which is charmonium-like (hidden-charm) and the recently-discovered $T_{cc}^+$, a doubly-charmed tetraquark candidate with a mass of about 3875 MeV/c² and a ccūd quark content [191]. Both hadrons can be interpreted as (loosely-bound) mesonic molecular states, namely $D^{*0}\bar{D}^0$ and $D^{*+}D^0$, respectively, as their masses lie narrowly below the thresholds for the respective channels. In such an interpretation these hadrons
<figure><img src="image_11.png" /><figcaption>**Fig. 9** Multiplicity dependence of the cross section ratio of $\chi_{c1}(3872)$ and $\psi$ meson (left) and of the $T_{cc}^+$ measured yield [192] (right)</figcaption></figure>

are consequently extended objects (of size up to 10 fm), but the case of compact tetraquark states is not excluded.

The study of exotic-hadron production can possibly also shed light on hadronization mechanisms, as suggested by the recent observation by LHCb [174] of an event-multiplicity dependence of the $\chi_{c1}(3872)$ meson production in pp collisions, shown in Fig. 9 (left). The measurement is described in the comover model [193] as break-up via interaction with other particles, modeled with an effective cross section dependent on the size of the hadron. A relatively-compact tetraquark is preferred in this case, but it is currently discussed how conclusive this comparison is [194,195]. The cross section ratio was observed to be larger in pPb collisions [196]. The multiplicity dependence of $T_{cc}^+$ production [192] suggests an increasing trend (Fig. 9 right), which is surprising, given the comparable size of the two hadrons.

Another very interesting exotic state, the X(6900), discovered in 2020 by LHCb [197] and recently confirmed by ATLAS [198] and CMS [199], is a candidate for a tetraquark of charm quarks and antiquarks. Studying it in the future as a function of multiplicity, challenging as it certainly is, will be an important component towards understanding the hadronization of complex exotic hadrons.

## 5.4 Next steps

In this section the possible next steps discussed at the workshop are presented (essentially a wish list for the experiments, clearly a long-term program):

* Other non-strange baryon than protons vs. multiplicity (which points to the $\Delta$ resonances, experimentally a very challenging measurement).

* Correlations (at low momentum) between baryons, both with light and heavy quarks.

* $\Lambda_c^+/D^0$ ratio at low multiplicities (and at high $p_T$) in pp collisions, to check if the value in $e^+e^-$ is recovered.

* Clarifying (experimentally) the apparent inconsistency in the $\Lambda_c^+/D^0$ dependence on rapidity.

* The measurement of $\Lambda_c$ elliptic flow in PbPb and its difference w.r.t. the D one would allow to infer the relevance of diquark degrees of freedom in the QGP, see Fig. 8.

* Search for new excited baryon states, both for the charm and bottom sector.

* More differential $\chi_{c1}(3872)$ measurements, in yield rather than ratio to $\psi(2S)$.

# 6 Energy loss and transport from small to large systems

Despite the impressive progresses on both the experimental and theoretical side in the last twenty years, a fully satisfactory description of in-medium energy loss and related phenomena has not been reached yet. Still challenging is the comprehensive description of $R_{AA}$ and flow data of several particles, which provide sensitivity to the energy-loss process and to its dependence on parton colour-charge and mass. Furthermore, the observation of collective effects in small systems without sizeable signs of energy loss is a major puzzle yet to be solved.

## 6.1 Signals for energy loss in small systems

One of the most exciting challenges posed by the experimental results obtained in small collision systems is the so-called $v_2 - R_{ppb}$ puzzle. Finite elliptic flow, $v_2 > 0$, is observed for heavy flavour hadrons [200,201] and even for jets [202,203], which could be naturally explained by partonic energy loss, while $R_{ppb}$ is compatible with nuclear PDF model calculations [204–207]. While such a behavior is predicted by the CGC model, it also predicts bottonium $v_2 (\Upsilon) > 0$ [208], which is not observed [209]. Other possible solutions to the $v_2 - R_{ppb}$ puzzle, which do not include energy loss, exist. For instance, AMPT predicts $v_2 > 0$ and $R_{ppb} \approx 1$ by means
of a parton escape mechanism (see e.g. [202]). Moreover, a Glasma phase alone could give “diffusion” and no energy loss resulting in $R_{ppb} > 1$ which moves back to unity through energy loss in a medium [210,211].

Taking into account experimental uncertainties, $R_{ppb}$ provides only relatively weak constraints on possible energy loss effects, and thus they cannot be completely excluded. Employing jet-hadron correlation measurements in pPb collisions for charged jet energies in the range 15–60 GeV/c, ALICE has provided a limit on the energy lost outside a jet cone of $R = 0.4$ of 0.4 GeV at 90% CL, ATLAS has extended the measurements to jet $p_T > 60$ GeV/c without observing a clear energy loss pattern. With the larger data samples projected to be collected during the LHC Run 3 for pp, pPb and OO collisions this limit can be pushed below 0.1 GeV [212].

Energy loss effects are expected to be larger in collisions producing a large particle multiplicity, which is interpreted as due to the high number of initial parton-parton interactions, and thus, high multiplicity events correspond to high initial energy densities. However, the analysis of such events is complicated due to selection biases [213]. An example is the jet-hadron azimuth angular correlation measurements in pp, where the observed broadening can be reproduced by PYTHIA8 calculations which do not include any medium effects [214].

New ideas for the analysis of pp collisions discussed during the workshop are a better control of the multiplicity dispersion over the rapidity covered by the signal and multiplicity measurements regions and search for signals of redistributed energy in the underlying event of jets. From the theory side, it would be important that models predicting a finite $v_2 > 0$ in small systems provide also the minimum parton energy loss that could explain the effect.

## 6.2 Role of the pre-equilibrium stage

Most model calculations of parton energy loss in heavy-ion collisions assume no interactions in the first ~1 fm/c, i.e. before the QGP formation. However, recent studies in the CGC and in effective kinetic theory predict a substantial $\hat{q}$ in the pre-hydrodynamics phases [210,215], which may have consequences for the determination of $\hat{q}$ in the QGP phase and the understanding of jet quenching (see e.g. [216]).

Recent energy loss calculations in heavy-ion collisions consider three simplified scenarios for the initial stages: (1) parton production time, $\tau_p$, and medium production time, $\tau_m$ are equal and 0, (2) $\tau_p = \tau_m = 1$ fm/c and (3) $\tau_p = 0$ and $\tau_m = 1$ fm/c. While all the scenarios can reproduce the $p_T$-dependence of the $R_{PbPb}$, albeit with different coupling parameters, once the coupling parameter is fixed, they yield a significantly different high-$p_T$ particle $v_2$. This highlights the importance of the treatment of energy loss in the initial stages to properly understand the $R_{PbPb}$ and high-$p_T$ $v_2$ data.

While a large $\hat{q}$ is predicted in the region 0–1 fm/c, so far no realistic medium-interactions have been considered in this pre-hydrodynamics phase. We think that efforts from the theory community are needed to understand the consequences for $v_2$ and energy loss. Moreover, we wonder whether small systems can be considered a proxy for pre-equilibrium effects since in these systems this phase dominates.

## 6.3 Energy dependence of $\hat{q}$

Another interesting theoretical result was brought to our attention during the workshop: rigorous analytical estimations within the framework of kinetic theory have revealed significant jet energy dependence of $\hat{q}$. For the first time, all possible diagrams for $2 \to 2$ [217] and $2 \to 3$ [218,219] quark-gluon scatterings are being thoroughly utilized to evaluate $\hat{q}$ within the Dynamical Quasi-Particle Model, which describes the QGP phase in the PHSD transport approach [220]. At a high-temperature limit obtained $E_{jet}$-dependence is in accordance with the well-known LO-HTL leading-log expression [221,222]. This energy dependence is sensitive to the choice of the coupling constant and it can be especially useful for further improvements of jet-quenching models, which usually do not take $E_{jet}$-dependence into account or fail to reproduce asymptotic leading-log behavior.

We think it is important to quantify the consequences for the in-medium parton shower evolution and jet shapes. Can the energy dependence be constrained from experimental data?

## 6.4 Energy loss in quarkonium production

The observation that in pp collisions prompt $J/\psi$ mesons are part of jets and carry, on average, only 50% of the jet energy came as a surprise, since event generators predict them to be produced almost isolated. A possible explanation is that $J/\psi$ pre-states are produced as colour-octets which radiate gluons evolving into a colour singlet. Also $J/\psi$ mesons produced as part of a parton shower with a $g \to c\bar{c}$ splitting are non-isolated. In both cases, the directly production of $J/\psi$ mesons can be modified by the presence of a QGP. Indeed, CMS observes that in PbPb collisions high-$p_T$ non-isolated $J/\psi$ are more suppressed than the isolated ones [223].

So far nothing is known about the nature of the particles produced in association with the $J/\psi$. We think it is important to characterise them measuring their $p_T$ and angular distributions (or studying jet shapes) in order to constrain the production mechanism. One should also search for modifications of the jet structure comparing PbPb and pp collisions. Do the same phenomena occur for $\Upsilon$ production? Moreover, the mechanisms responsible for high-$p_T$ $J/\psi$ suppression might also play a role in charmed hadron suppression, since
the $c\bar{c}$ pair stays for some time in a colour octet state. We propose to study charm suppression as a function of the $D\bar{D}$ opening angle.

## 6.5 Better constraints on heavy flavor diffusion

In heavy-ion collisions, charm hadron spectra in the low-$p_T$ region and $v_2$ are used to constrain the charm diffusion coefficient. However, collisional energy loss is not the only key ingredient for the calculations. It has been shown that without radiative energy loss, the predicted $R_{AA}$ is too high and $v_2$ too low [218,224]. Another important ingredient for the interpretation of the experimental results is the relation between the c-quark and the hadron $p_T$. This relation is different for hadronisation via recombination and fragmentation. For recombination the momentum of the hadron is larger than that of the charm quark and it also inherits part of the flow of the light hadron with which it combines [225].

We think that in order to mitigate the influence of the hadronisation mechanism on the determination of the diffusion coefficient one should limit the analysis to models which can describe simultaneously $R_{AA}$, $v_2$ and the baryon-to-meson ratio $\Lambda_c/D$. Further improvements will be obtained by extending the measurements to lower $p_T$. In the future, additional observables such as $D^0$-hadron $v_2$ with event shape engineering (ESE), $D^0 v_2${4}, and $v_3$ as well as correlations of $D^0$ versus pion $v_n$ in ESE classes will improve the constraints on the HF diffusion [226,227]. Moreover, we emphasise the importance for models to provide uncertainty bands since they have a large effect on $\chi^2/ndf$.

Ultimately, one will use B mesons (LHC Run 3 and beyond, sPHENIX) to measure the heavy-flavour diffusion coefficient. For b-quarks, there are less uncertainties in the transport description (Boltzmann, Langevin) [228,229]. Furthermore, the larger mass brings the calculations closer to the infinite-mass approximation used in lattice QCD calculations of the diffusion coefficient [230].

## 6.6 Signals for heavy flavor thermalisation

Modifications of charm hadron spectra and elliptic flow represent only indirect signs of charm thermalisation. Thermalisation implies that the information about the initial conditions of the production gets blurred completely. Ultimately, this can be demonstrated with $D\bar{D}$ azimuthal angle correlations. The relative height of the nearside and back-to-back peaks is sensitive to energy loss and the width of these peaks to the extent of thermalisation, which would imply full isotropisation at low $p_T$ (see e.g. [231]).

These measurements are very statistics hungry and might be only possible with the ALICE3 detector in LHC Run 5 [18]. It would be interesting to see how much one can already learn from charm hadron-inclusive hadron correla-

tions. On the one hand, the decorrelation from the decay renders the measurement less sensitive, on the other hand one gains a large factor in sample size. Moreover, experimental performance studies need calculations with state-of-the-art models.

## 6.7 Challenges for future high precision measurements

We close this chapter with a section on more general considerations concerning the challenges for future high precision measurements. Today one witnesses a significant variability among models with different description of partonic interactions in the medium, medium evolution, hadronisation, hadronic phase, and nuclear PDFs. Important questions to be answered by theorists are: whether the choice and number of modelling parameters are under sufficient control, and whether the observables currently provided by the experiments are already optimal.

Results from Bayesian analyses for parameter constraints are crucially sensitive to uncertainties: wrong uncertainties equals to wrong results. Hence, experimentalists must assure that the uncertainties they estimate are under sufficient control for the high-precision era. For instance, the tendency to choose systematic uncertainties conservatively can lead to overestimated overall uncertainties entering the fits. A more complicated issue are correlations among systematic uncertainties, since in many cases they are not evaluated within the same measured distribution and very rarely among different measurements performed by different groups. Last but not least, models are sometimes used to evaluate the systematic uncertainties, which also introduces correlations between the experimental uncertainties and the models.

# 7 QCD and astrophysics

Along studies of QCD that cast light on the structure and behavior of nucleons and nuclei, various connections with astroparticle and astrophysics-related “puzzles” have been discussed during the workshop. A better understanding of phenomena in soft QCD is expected to lead to breakthroughs in these fields. We summarize the input needed from collider experiments for the applications in astroparticle and astrophysics in Table 1.

## 7.1 Hyperon puzzle for neutron stars

The interaction of hyperons (Y) with nucleons (N) is one of the keys to understand the composition of the most dense objects in our Universe: neutron stars (NS) [232,233]. These are characterised by large masses ($M \approx 1.2-2.2$ solar masses $M_{\odot}$) and small radii ($R \approx 9-13$ km) [234–236]. In the standard scenario, the gravitational pressure is counter-balanced
<table><caption>**Table 1** Desired input for astroparticle and astrophysics from colliders</caption><thead><tr><th>Motivation</th><th>Experiments</th><th>Desirables</th></tr></thead><tbody><tr><td>Anti-hyperon puzzle in neutron stars</td><td>ALICE, LHCb</td><td>Precision measurements, rapidity dependence?</td></tr><tr><td>Interactions of hyperons (Y) and nucleons (N); NΣ, Λφ; three-body interactions NNY, NYY, YYY</td><td>LHCb SMOG2, AMBER, NA61</td><td>LHC run at 450 GeV</td></tr><tr><td>Searches for dark matter and primordial antimatter</td><td>LHCb SMOG2, AMBER, NA61</td><td>LHC run at 450 GeV</td></tr><tr><td>In pH<sub>2</sub>, pHe, pD<sub>2</sub>: production of d̄/He, p̄, e<sup>+</sup>, π<sup>0</sup> → γγ, η → γγ</td><td>ALICE, LHCb, LHCf, FOCAL</td><td>Precision measurements over wide range in η and √s in p-ion systems, especially pO</td></tr><tr><td>Muon puzzle in air showers</td><td>ALICE, LHCb, LHCf, FOCAL</td><td>Precision measurements over wide range in η and √s in p-ion systems, especially pO</td></tr><tr><td>Charged particle spectra; pO and πO inelastic cross-section; multiplicity-dependent hadrochemistry: R = E(π<sup>0</sup>)/E<sub>tot</sub>, p/π, Λ/K<sub>S</sub><sup>0</sup>, Ξ/Λ, ...</td><td>LHCb, FASER</td><td>Forward acceptance at large rapidity required to reach small x, LHC run at 450 GeV</td></tr><tr><td>Atmospheric background to astro-neutrinos</td><td>LHCb, FASER</td><td>Forward acceptance at large rapidity required to reach small x, LHC run at 450 GeV</td></tr><tr><td>D<sup>0</sup>, D<sup>+</sup> production cross-sections at forward rapidity; multiplicity-dependent charm production; oxygen gluon nPDF at x ~ 10<sup>-7</sup></td><td>LHCb, FASER</td><td>Forward acceptance at large rapidity required to reach small x, LHC run at 450 GeV</td></tr></tbody></table>

by the Fermi pressure of neutrons in the core, which, along with electrons, are the only particles from the parent star remnants.

The high-density environment ($\rho \approx 4 \rho_0$, with $\rho_0$ being the nuclear density) assumed to occur in the interior of NS leads to an increase in the Fermi energy of the nucleons, resulting into the appearance of new degrees of freedom such as hyperons. Nonetheless, this energetically-favoured production of strange hadrons induces a softening of the Equation of State (EoS) incompatible with the current highest mass limit from experimental observations of close to $3 M_\odot$ [237–239]. For this reason, the presence of hyperons inside the inner cores of NS is still under debate, and this “hyperon puzzle” is far from being solved [240,241]. A possible way out is represented by two-body and three-body repulsive YN and YNN interactions. In both cases, a sufficiently strong YN or YNN repulsive interaction can push the appearance of hyperons to larger densities, limiting the possible presence of these particle species inside NS, stiffening the EoS and leading to larger star masses.

Currently, ALICE femtoscopic measurements in small systems (pp, pPb) for baryon–baryon pairs involving hyperons deliver the most precise data on the residual strong interaction between nucleons and strange hadrons [242–254].

The results obtained in recent years from femtoscopic measurements in small colliding systems have proven that femtoscopy can play a central role in understanding the dynamics between hyperons and nucleons in vacuum. For example, the ALICE measurements on pΞ<sup>-</sup> pairs [243,244] confirmed a strong attractive interaction between these two hadrons and provided a direct confirmation of lattice potentials [255], leading to a better understanding of the interaction among hyperons and nucleons.

A comparison between hadronic models and these data is necessary in order to constrain calculations at finite density and to pin down the behavior of hyperons in a dense matter environment. The great possibility to investigate, within the femtoscopy technique, different Y–N interactions and to extend the measurements to three-body forces, can finally provide quantitative input to the long-standing hyperon puzzle.

## 7.2 Muon puzzle in atmospheric showers

Cosmic rays (CRs) are fully ionized nuclei that arrive at Earth with relativistic kinetic energies. At energies above 100 TeV, CRs are observed indirectly by ground-based experiments through extensive air showers, particle showers produced in Earth’s atmosphere. Air shower features are used to determine the direction, energy, and mass of the cosmic ray. Two main features of an air shower are used to estimate the mass, its depth of shower maximum $X_{max}$, and the number of muons $N_\mu$ produced in the shower. The bulk of muons is produced at the end of a hadronic cascade dominated by interactions in the non-perturbative regime of QCD, which are modeled by effective theories like the core-corona model EPOS-LHC [258,259], and the parton shower model Pythia with string-string interactions [154,260,261]. Especially thanks to the latest generation of large air shower experiments that measure air showers in unprecedented detail, a significant muon deficit is observed in current state-of-the-art simulations based on $X_{max}$ in comparison to observations in real air showers (Fig. 10). This muon production anomaly
<figure><img src="image_12.png" /><figcaption>**Fig. 10** Compilation of muon measurements converted to the abstract z-scale for the hadronic interaction model EPOS-LHC. The z-scale is given by $z = (\ln(N_{\mu}) - \ln(N_{\mu,p}))/((\ln(N_{\mu}, \text{Fe}) - \ln(N_{\mu,p})))$, where $N_{\mu}$ is the measured muon content and $N_{\mu,p}$ and $N_{\mu, \text{Fe}}$ are predictions for proton and iron CRs, respectively. The energy scales of the experimental datasets shown here have been cross-calibrated as described in Ref. [256]. The expected value $z_{\text{mass}}$ based on the cosmic-ray mass composition estimated by the GSF model [257] is subtracted. Also shown are $z_{\text{mass}}$ values computed from $X_{\text{max}}$ measurements by the Pierre Auger Observatory (grey band)</figcaption></figure>

is called the “Muon Puzzle in air showers”, more information can be found in Ref. [262]. The produced number of muons is very sensitive to the energy fraction carried away by photons, which are primarily produced in $\pi^0$ decays [263–265]. Small changes of the $\pi^0$-fraction among all produced hadrons in the forward region have a large effect on muon production [264,266].

To solve the Muon Puzzle, forward measurements of the light-flavor hadron production at small transverse momentum in minimum bias events are needed, in particular of hadro-chemistry. Indirectly, precision measurements of the inelastic cross-section in pA collisions are also needed, and ideally also the proton-$\pi$ inelastic cross-section, which affect predictions for $X_{\text{max}}$. Particularly important are measurements by forward experiments like LHCf [267–269], CMS with CASTOR [270], and LHCb [27, 271, 272]. A key ingredient to resolve the muon puzzle could be multiplicity-dependent strangeness enhancement observed by ALICE [108] in the central collision region. LHCb found evidence for this effect also in the forward region [273]. In the near future, the study of soft hadron production in pOcollisions at the LHC is another key step forward, which provides the best reference for cosmic proton with air, and will greatly reduce the uncertainty coming from interpolations of pp and pPb data.

### 7.3 Production of secondary particles in the galaxy

With the last generation of particle detectors in space (PAMELA [274], AMS-02 [275], DAMPE [276], CALET [277], Fermi-LAT [278]) and with the future ones like GAPS [279], physics of CRs and gamma rays ($\gamma$ rays) is more and more a precision discipline. The fluxes of positrons ($e^+$) [275] and light antinuclei like antiprotons ($\bar{p}$) [277], and potentially in the near future antideuteron ($\bar{d}$) and antihelium ($\bar{He}$), are being measured with percent uncertainties in a wide energy range and compared to the expectations from secondary-only production, namely primary CR interactions with the interstellar medium (ISM), to evidence primary astrophysical sources or contributions from dark matter annihilation or decay [280–283]. To improve the secondary production description, and infer conclusions on the primary contributions, the uncertainties on hadronic cross-section, currently contributing as $\pm 20\%$ [284] to the secondary $\bar{p}$ flux, need to scale down to the experimental precision. Cross-section measurements in the pp channel, dominating the secondary production [284–286], or either with the CR projectile or the ISM target replaced by helium (Hep, pHe, and HeHe), are hence needed. Both direct and indirect production have to be addressed. For $\bar{p}$, with the dominant contributions coming from prompt emission, hyperon-induced and the possible isospin asymmetry for antineutron production have to be constrained on data [287]. For $\bar{d}$ and $\bar{He}$ nuclei, suppressed in the CRs-ISM production at low energies with respect to DM signal predictions [288,289], in addition to direct production via coalescence, a possible contribution from $\bar{\Lambda}_b$ baryon decays [290] has to be constrained. For $e^+$ above 1 GeV energy, the flux data by AMS-02 are higher than the current secondary-only predictions [291], that could be improved with measurements of the $\pi^\pm$ and $K^\pm$ spallation cross-sections from Hep, pHe, and HeHecollisions.

Most of the $\gamma$ rays produced by hadronic interactions detected by Fermi-LAT [292] originate instead from the $\pi^0 \to \gamma\gamma$ decay and new data on the Lorentz invariant cross section for $\pi^0$ production are needed [286] to bring down the $\sigma(pp \to \pi^0 X)$ uncertainty to the same level as the Fermi-LAT statistical errors. Particularly, measurements for $p_T \lesssim 1$ GeV, large coverage in $x_F = 2p_L/\sqrt{s}$ and beam energies in the laboratory frame from a few tens of GeV to few TeV, both for pp and pHecollisions, would contribute most in decreasing the uncertainties.

In parallel to the theoretical progresses, experiments from different facilities are providing relevant measurements. At the SPS accelerator at CERN, NA61 [293] and COMPASS [294] fixed-target experiments are being upgraded to extend and improve $\bar{p}$, $e^+$, and $\gamma$ production cross-section measurements from pp and pHeinteractions. Leveraging on the injection of gases in the LHC beam-pipe with a system called SMOG, the LHCb experiment has also been oper-
ated in fixed-target mode since 2015, providing the first mea-
surements for both direct [295] and hyperon-induced [296]
$\sigma(\text{pHe} \rightarrow \bar{p}X, \sqrt{s_{\text{NN}}} = 110 \text{ GeV})$. Measurements of $\bar{d}$,
$K^{\pm}$, $\pi^{0}$, $\eta$, $e^{+}$ in the same dataset and of $\bar{p}$ at lower beam
energies are also ongoing or feasible with the SMOG sys-
tem. With the upgrade of the LHCb fixed-target programme,
namely the installation of a gas storage cell upstream of the
LHCb nominal interaction point and a new gas feed sys-
tem [297], hydrogen and deuterium could be injected as
well, shedding light on the possible isospin asymmetry in the
antineutron-to-antiproton production. Moreover, concurrent
data-taking with the beam-beam collisions was demonstrated
to be feasible [298], resulting in larger collected data samples
and more precise cross-section measurements [299]. All of
this is resulting in a unique facility for QCD measurements
at the LHC and, notably, for its connection to astrophysics.

8 Summary and outlook

The 4th edition of the workshop QCD challenges from pp to
AA collisions took place in Padua, Italy, from 13th to 17th
February 2023. It was attended by 65 researchers, includ-
ing both “senior” experimental and theoretical physicists and
young researchers (about 15%). The agenda consisted of
plenary sessions alternated to parallel round-table sessions
organised in parallel on six tracks, with a few sessions ded-
icated to inter-track discussions. The workshop format was
very effective in stimulating discussions. In the round-table
sessions there was ample time for following up the topics
presented during plenary sessions, the questions raised, as
well as for deepening items and ideas that session conveners
had prepared in advance.

In general, the necessity emerging in the last years of shift-
ing away from the paradigm considering proton–proton and
heavy-ion collisions as systems featuring different and unre-
lated properties was further established during the workshop.
Understanding the observation of collective flow behavior as
well as the dependence of strangeness production on the event
charged-particle multiplicity is considered of pivotal impor-
tance. It was noted in particular the need for phenomenologi-
cal approaches that take into account non-equilibrium effects
that become more and more important going from large to
small collision systems. The performance of core-corona
approaches in describing the hadrochemistry and kinemat-
ics across collision systems was discussed, together with the
possible relevant role played by the core in pp collisions.
Gaining a realistic three-dimensional description of the entire
interaction area of the collision was highlighted as a neces-
sary step for interpreting rapidity-differential measurements
and observables like $v_1$ and particle asymmetries sensitive
to effects and correlations spanning over wide rapidity inter-
vals. Full 3D simulations are also crucial for the comparison

to the elliptic and triangular flow data in small systems. More-
over, studying the correlation between fluctuations in mean
$p_T$ and elliptic-flow magnitude could help resolving initial-
state effects like initial momentum correlations in small sys-
tems.

So far no evidence was found of energy loss in small sys-
tems while positive elliptic-flow values have been observed
for several high-energy probes, like heavy-flavour, quarko-
nia and jets. Also with jet-hadron correlations, expected to
be less affected by the experimental uncertainties, only upper
limits to the out-of-cone radiation could be set in pPb col-
lisions. A better control of the multiplicity dispersion over
large rapidity regions and searches for signals of the redis-
tributed energy in the underlying event were discussed as
possible directions to solve the $R_{pPb} - v_2$ puzzle. On the
theory side, it would be important that models predicting
positive $v_2$ specified also the minimum parton energy loss
required to generate it. Sensitivity to possible energy loss
in the final state requires understanding of the initial state
also for the production of high-energy partons. At high colli-
sion energy the colliding nuclei can be represented as dense
gluonic systems. Saturation effects, which can be described
by the Colour Glass Condensate effective field theory, can
become relevant when low Bjorken-$x$ values are involved.
As reviewed during the workshop, in recent years, substan-
tial progress has been made to constrain this low-$x$ regime
with measurements of particle production at forward rapidity
in pA collisions and of photoproduction of vector mesons in
ultraperipheral collisions. These data stimulated significant
theoretical developments, especially to overcome the obsta-
cles posed by the poor stability of perturbative calculations of
the amplitude to the inclusion of heavy-vector meson data in
global PDF analyses. However, attempts of extracting proton
PDF by relating PDF to GPD and using HERA and LHCb
J/ψ photoproduction data show large differences from PDF
determined using LHCb $D^0$-meson data. This discrepancy
indicates the need of additional theoretical work but also
that new experimental data including those from pO and OO
runs at the LHC and forward direct-photon production mea-
surements that detector upgrades will allow, will be vital to
observe and understand nonlinear QCD effects.

Similarly to the appearance of long-range correlations
and flow-like effects, the hadronisation process is a fur-
ther important case in which concepts originally introduced
to describe hadron production yields and momentum dis-
tributions in heavy-ion collisions turned out to be effec-
tive for proton–proton and proton–Pb collisions. In partic-
ular, models including quark-recombination as a hadroni-
sation process or based on a statistical approach agnostic
of any partonic phase describe heavy-flavour baryon pro-
duction significantly better than calculations that employ
fragmentation functions constrained to reproduce $e^+e^-$ data
and that strongly underestimate the baryon-to-meson ratios.
Models based on string-fragmentation expect an enhance-
ment of baryon-to-meson ratios in hadronic collisions only if
colour reconnection is implemented more realistically, going
beyond the leading-colour approximation that is sufficient for
describe $e^+e^-$ data. It is expected that the experiments will
make an effort in the incoming years to measure the produc-
tion yields of $\Sigma_c(2455)^{0,++}$ and of the charm-and-strange
baryons $\Xi_c^{0,+}$ and $\Omega_c^0$ with improved precision over a wider
transverse momentum range and better profiling their evo-
lution with multiplicity across collision systems. A precise
measurement of beauty baryon production at midrapidity,
which is currently missing, will allow to determine whether
baryons are effectively less produced at forward rapidity,
resolving the apparent tension between ALICE and LHCb
charm baryon data. On the theoretical side, a careful study
of the expected multiplicity dependence in small systems as
well as a study of the possible rapidity dependence of the
baryon-to-meson ratios is desirable. Particle angular corre-
lations are also envisaged as powerful probes of the hadro-
nisation dynamics. Aside “standard” heavy-flavour baryons,
studying the production yields and the properties of exotic
hadrons with heavy-quark content may add further insight
into hadronisation, as well as into the nature and internal
structures of these states.

For “large” collision systems and for the characterisation
of the QGP, polarisation measurements are seen as a very
promising way to investigate vorticity phenomena that can
arise in the medium both as a global effect related to the
medium rotation as well as localised ones in response to the
passage of a high-energy particle. The possibility to con-
nect quark-spin alignment to hadron polarisation is based on
the assumption that a sizeable fraction of hadrons is formed
via coalescence, inheriting the spin of the system of recom-
bining quarks. Polarization measurements can therefore also
probe the hadronisation process. Constraining the hadroni-
sation is also fundamental for making a step further in the
study of in-medium energy loss and heavy-quark transport.
In particular, for determining the spatial diffusion coeffi-
cient from charm-hadron data, it is mandatory to rely only
on models that can simultaneously describe $R_{AA}$, $v_2$ and
the $\Lambda_c^+/D^0$ ratio. Additional constraints could be set by
new and precise measurements of D-meson flow coefficients
with Event Shape Engineering or of direct correlations of D-
meson and pion flow coefficients. Eventually, thanks to the
detector upgrades and the expected new measurements of
beauty flow and $R_{AA}$, the spatial-diffusion coefficient could
be determined with beauty hadrons. This possibility has the
advantage that the uncertainties in both transport models and
lattice-QCD calculations are smaller on beauty observables
than on charm ones.

Concerning in-medium energy loss, the energy depen-
dence of the $\hat{q}$ coefficient, which encodes the scattering
power of the medium, was discussed in the workshop. Under-

standing its consequences for in-medium parton shower and
jet shapes and how experimental data could constrain it are
important challenges for the future. It was also highlighted
that treating energy loss in the first instants after the colli-
sion, within the pre-equilibrium stage of the QGP, an interval
often neglected in energy-loss calculations, could be impor-
tant for the simultaneous description of $R_{AA}$ and $v_2$ and for
the determination of $\hat{q}$. It was also proposed to investigate
whether small systems could provide a proxy for addressing
the pre-equilibrium stage. As a different direction to study
energy loss and in particular medium-induced radiation, the
potential of jet substructure measurements for determining
the characteristic angle of medium-induced emissions, $\theta_c$
was extensively debated in the workshop. Observables based
on clustering-tree, like the dynamically groomed jet radius,
and on energy-flow, like energy-energy correlators, were dis-
cussed. For the latter, measurements in $\gamma$-jet events were
advertised as promising thanks to mitigation of selection-
bias effects. In general, isolating the $\theta_c$-related dynamics
remains a very challenging task: jet-substructure observables
in high-$p_T$ jets could more naturally match the requirements
identified as necessary to fulfill this task, i.e. having a cal-
culable vacuum baseline, resilience to the underlying event,
and reduced sensitivity to medium response.

At the workshop some QCD aspects influencing the under-
standing of astrophysics and cosmic-ray physics were also
discussed. Far from being solved is the "hyperon puzzle" for
neutron stars. It was remarked that there is a need of fur-
ther femtoscopic measurements for clarifying whether two-
body or three-body repulsive interactions exist and limit the
formation of hyperons in the core of neutron stars, stiffen-
ing the equation of state and allowing larger star masses.
Another "puzzle" discussed in the workshop concerns the
deficit of muons in state-of-the-art simulations compared to
observations in air showers. Important inputs from particle-
accelerator experiments are forward measurements of light-
flavor hadron production at small transverse momentum,
in particular studies of hadrochemistry, including profiling
strangeness enhancement as a function of the event multiplic-
ity. It was noted that a future run of pO collisions at the LHC
could provide a better reference for cosmic rays than current
pp and pPb data. Concerning the search of dark-matter sig-
nals, the constantly improving precision of measurements
of charged cosmic rays and $\gamma$ rays from particle detectors
in space is opening a new era. However, in order to draw
conclusions on the primary sources which the description of
the secondary production in cosmic-ray interactions with the
interstellar medium must be improved. To this purpose, pre-
cise measurements of the production cross section of $\bar{p}, \bar{d}, \bar{\pi}, \bar{e}, \bar{\mu}, \bar{e}^+, \bar{\mu}^+$ in pp as well as in pHe, Hep, and HeHe collisions
possible with the NA61 upgrade, the COMPASS++/AMBER
experiment, and the LHCb SMOG program, are fundamental
to reduce the uncertainties on hadronic cross sections. Con-
straining light antinuclei production in beauty-hadron decay with LHCb or in the future with ALICE3 will also be important.

We would like to conclude this document with few final remarks on the effectiveness of the workshop structure. With respect to previous editions, an effort was made to enlarge the participation of physicists working on QCD-related topics outside the heavy-ion physics community. This was a major objective of this edition that we think should be pursued also in the future. The introduction of a session connecting QCD at colliders with astrophysics themes was also a new ingredient of this edition appreciated by the participants. For young researchers, the participation in the round-table sessions was a particularly valuable experience. During the organisation of the event, it was deliberately chosen to give session conveners ample freedom to define the topics they preferred to deepen during the sessions and, to a large extent, define the speakers. While this freedom assures lively and focused discussions, it also makes it more complex to assure some continuity between the topics discussed in different editions. However, we do not consider this as a significant limitation, rather as a feature characterising what should be the expectations and the outcome of a similar event, which, on top of stimulating detailed discussions on well-established open points in the field, is also meant to start new collaborations and directions of research. The organisation of a 1-day online “prequel” event a couple of months before the workshop could be useful to make new participants familiar with the workshop format and to sharpen the discussion topics, also defining points and questions to be addressed in the workshop.

**Acknowledgements** We acknowledge the financial support of the Italian Institute for Nuclear Physics (INFN) and the University of Padua. This project has received funding from the European Union's Horizon 2020 research and innovation programme under grant agreement N. 824093 (STRONG 2020 project). We acknowledge the support from the project "CLASH: Pinning down the origin of collective effects in small collision systems". The work was further supported by the Czech Science Foundation (GA ČR), project No. 22-27262S. We acknowledge also the support of the Ministry of Education and Science and National Science Centre (Poland) and the U.S. National Science Foundation. L.O. acknowledges the Next Generation action of the European Commission and the MUR funding (PNRR Missione 4) under the HEFESTUS project. P.P. acknowledges financial support from the Academy of Finland project 330448, the Center of Excellence in Quark Matter of the Academy of Finland and the European Research Council project ERC-2018-ADG-835105 YoctoLHC. C.A. has received funding from the European Union's Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement No 893021 (JQ4LHC). Finally we acknowledge the grant ERC-CZ LL2327.

**Data Availability Statement** This manuscript has no associated data or the data will not be deposited. [Authors' comment: Data sharing not applicable to this article as no datasets were generated during the current study].

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation,

distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Funded by SCOAP³.

## References

1. F. Karsch, Lattice simulations of the thermodynamics of strongly interacting elementary particles and the exploration of new phases of matter in relativistic heavy ion collisions. J. Phys. Conf. Ser. **46**, 122–131 (2006). https://doi.org/10.1088/1742-6596/46/1/017. arXiv:hep-lat/0608003 [hep-lat]

2. S. Borsanyi et al., Is there still any $T_c$ mystery in lattice QCD? Results with physical masses in the continuum limit III. JHEP **1009**, 073 (2010). https://doi.org/10.1007/JHEP09(2010)073. arXiv:1005.3508 [hep-lat]

3. S. Borsanyi, Z. Fodor, C. Hoelbling, S.D. Katz, S. Krieg et al., Full result for the QCD equation of state with 2+1 flavors. Phys. Lett. B **730**, 99–104 (2014). https://doi.org/10.1016/j.physletb.2014.01.007. arXiv:1309.5258 [hep-lat]

4. A. Bazavov, T. Bhattacharya, M. Cheng, C. DeTar, H.T. Ding et al., The chiral and deconfinement aspects of the QCD transition. Phys. Rev. D **85**, 054503 (2012). https://doi.org/10.1103/PhysRevD.85.054503. arXiv:1111.1710 [hep-lat]

5. M.C. Abreu et al., Evidence for deconfinement of quarks and gluons from the J/psi suppression pattern measured in Pb + Pb collisions at the CERN SPS. Phys. Lett. B **477**, 28–36 (2000). https://doi.org/10.1016/S0370-2693(00)00237-9

6. E. Andersen et al., Strangeness enhancement at mid-rapidity in Pb Pb collisions at 158-A-GeV/c. Phys. Lett. B **449**, 401–406 (1999). https://doi.org/10.1016/S0370-2693(99)00140-9

7. I. Arsene et al., Quark gluon plasma and color glass condensate at RHIC? The perspective from the BRAHMS experiment. Nucl. Phys. A **757**, 1–27 (2005). https://doi.org/10.1016/j.nuclphysa.2005.02.130. arXiv:nucl-ex/0410020

8. K. Adcox et al., Formation of dense partonic matter in relativistic nucleus-nucleus collisions at RHIC: experimental evaluation by the PHENIX collaboration. Nucl. Phys. A **757**, 184–283 (2005). https://doi.org/10.1016/j.nuclphysa.2005.03.086. arXiv:nucl-ex/0410003

9. B.B. Back et al., The PHOBOS perspective on discoveries at RHIC. Nucl. Phys. A **757**, 28–101 (2005). https://doi.org/10.1016/j.nuclphysa.2005.03.084. arXiv:nucl-ex/0410022

10. J. Adams et al., Experimental and theoretical challenges in the search for the quark gluon plasma: the STAR Collaboration's critical assessment of the evidence from RHIC collisions. Nucl. Phys. A **757**, 102–183 (2005). https://doi.org/10.1016/j.nuclphysa.2005.03.085. arXiv:nucl-ex/0501009

11. Y. Schutz, U.A. Wiedemann (eds.) Quark Matter. Proceedings, 22nd International Conference on Ultra-Relativistic Nucleus-Nucleus Collisions, Quark Matter 2011, Annecy, France, May 23–28, 2011, vol. 38, p. 120301 (2011). https://doi.org/10.1088/0954-3899/38/12/120301

12. G. Roland, K. Safarik, P. Steinberg, Heavy-ion collisions at the LHC. Prog. Part. Nucl. Phys. **77**, 70–127 (2014). https://doi.org/10.1016/j.ppnp.2014.05.001
13. P. Braun-Munzinger, V. Koch, T. Schäfer, J. Stachel, Properties of hot and dense matter from relativistic heavy ion collisions. Phys. Rep. **621**, 76–126 (2016). https://doi.org/10.1016/j.physrep.2015.12.003. arXiv:1510.00442 [nucl-th]

14. The ALICE experiment—a journey through QCD (2022). arXiv:2211.04384 [nucl-ex]

15. A. Adare et al., An Upgrade Proposal from the PHENIX Collaboration (2015). arXiv:1501.06197 [nucl-ex]

16. C. Ahdida et al., Letter of Intent: the NA60+ experiment (2022). arXiv:2212.14452 [nucl-ex]

17. R. Abdul Khalek et al., Science requirements and detector concepts for the electron-ion collider. EIC Yellow Report. Nucl. Phys. A **1026**, 122447 (2022). https://doi.org/10.1016/j.nuclphysa.2022.122447. arXiv:2103.05419 [physics.ins-det]

18. Letter of intent for ALICE 3: a next-generation heavy-ion experiment at the LHC (2022). arXiv:2211.02491 [physics.ins-det]

19. R. Abdul Khalek, R. Gauld, T. Giani, E. R. Nocera, T.R. Rabemananjara, J. Rojo, nNNPDF3.0: evidence for a modified partonic structure in heavy nuclei. Eur. Phys. J. C **82**(6), 507 (2022). https://doi.org/10.1140/epjc/s10052-022-10417-7. arXiv:2201.12363 [hep-ph]

20. K.J. Eskola, P. Paakkinen, H. Paukkunen, C.A. Salgado, EPPS21: a global QCD analysis of nuclear PDFs. Eur. Phys. J. C **82**(5), 413 (2022). https://doi.org/10.1140/epjc/s10052-022-10359-0. arXiv:2112.12462 [hep-ph]

21. P. Duwentäster, T. Ježo, M. Klasen, K. Kovařík, A. Kusina, K.F. Muzakka, F.I. Olness, R. Ruiz, I. Schienbein, J.Y. Yu, Impact of heavy quark and quarkonium data on nuclear gluon PDFs. Phys. Rev. D **105**(11), 114043 (2022). https://doi.org/10.1103/PhysRevD.105.114043. arXiv:2204.09982 [hep-ph]

22. I. Helenius, M. Walt, W. Vogelsang, NNLO nuclear parton distribution functions with electroweak-boson production data from the LHC. Phys. Rev. D **105**(9), 094031 (2022). https://doi.org/10.1103/PhysRevD.105.094031. arXiv:2112.11904 [hep-ph]

23. L.D. McLerran, R. Venugopalan, Gluon distribution functions for very large nuclei at small transverse momentum. Phys. Rev. D **49**, 3352–3355 (1994). https://doi.org/10.1103/PhysRevD.49.3352. arXiv:hep-ph/9311205

24. R. Aaij et al., Study of prompt D⁰ meson production in pPb collisions at $\sqrt{s_{NN}}$ = 5 TeV. JHEP **10**, 090 (2017). https://doi.org/10.1007/JHEP10(2017)090. arXiv:1707.02750 [hep-ex]

25. K.J. Eskola, P. Paakkinen, H. Paukkunen, C.A. Salgado, EPPS16: nuclear parton distributions with LHC data. Eur. Phys. J. C **77**(3), 163 (2017). https://doi.org/10.1140/epjc/s10052-017-4725-9. arXiv:1612.05741 [hep-ph]

26. R. Abdul Khalek, J. J. Ethier, J. Rojo, G. Weelden, nNNPDF2.0: quark flavor separation in nuclei from LHC data. JHEP **09**, 183 (2020). https://doi.org/10.1007/JHEP09(2020)183. arXiv:2006.14629 [hep-ph]

27. R. Aaij et al., Measurement of the nuclear modification factor and prompt charged particle production in pPb and pp collisions at $\sqrt{s_{NN}}$ = 5 TeV. Phys. Rev. Lett. **128**(14), 142004 (2022). https://doi.org/10.1103/PhysRevLett.128.142004. arXiv:2108.13115 [hep-ex]

28. R. Aaij et al., Nuclear modification factor of neutral pions in the forward and backward regions in p-Pb collisions. Phys. Rev. Lett. **131**(4), 042302 (2023). https://doi.org/10.1103/PhysRevLett.131.042302. arXiv:2204.10608 [nucl-ex]

29. S. Acharya et al., Transverse momentum spectra and nuclear modification factors of charged particles in pp, p-Pb and Pb-Pb collisions at the LHC. JHEP **11**, 013 (2018). https://doi.org/10.1007/JHEP11(2018)013. arXiv:1802.09145 [nucl-ex]

30. S. Acharya et al., Neutral pion and η meson production in p-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Eur. Phys. J. C **78**(8), 624 (2018). https://doi.org/10.1140/epjc/s10052-018-6013-8. arXiv:1801.07051 [nucl-ex]

31. S. Acharya et al., Nuclear modification factor of light neutral-meson spectra up to high transverse momentum in p-Pb collisions at $s_{NN}$ = 8.16 TeV. Phys. Lett. B **827**, 136943 (2022). https://doi.org/10.1016/j.physletb.2022.136943. arXiv:2104.03116 [nucl-ex]

32. V. Khachatryan et al., Charged-particle nuclear modification factors in PbPb and pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. JHEP **04**, 039 (2017). https://doi.org/10.1007/JHEP04(2017)039. arXiv:1611.01664 [nucl-ex]

33. M.S. Abdallah et al., Evidence for nonlinear gluon effects in QCD and their mass number dependence at STAR. Phys. Rev. Lett. **129**(9), 092501 (2022). https://doi.org/10.1103/PhysRevLett.129.092501. arXiv:2111.10396 [nucl-ex]

34. C. Aidala et al., Nonperturbative transverse momentum broadening in dihadron angular correlations in $\sqrt{s_{NN}}$ = 200 GeV proton-nucleus collisions. Phys. Rev. C **99**(4), 044912 (2019). https://doi.org/10.1103/PhysRevC.99.044912. arXiv:1809.09045 [hep-ex]

35. M. Aaboud et al., Dijet azimuthal correlations and conditional yields in pp and p+Pb collisions at $s_{NN}$=5.02 TeV with the ATLAS detector. Phys. Rev. C **100**(3), 034903 (2019). https://doi.org/10.1103/PhysRevC.100.034903. arXiv:1901.10440 [nucl-ex]

36. J. Jalilian-Marian, A.H. Rezaeian, Prompt photon production and photon-hadron correlations at RHIC and the LHC from the Color Glass Condensate. Phys. Rev. D **86**, 034016 (2012). https://doi.org/10.1103/PhysRevD.86.034016. arXiv:1204.1319 [hep-ph]

37. Letter of Intent: A Forward Calorimeter (FoCal) in the ALICE experiment (2020)

38. C.M. LHCb Collaboration, Framework TDR for the LHCb Upgrade II - Opportunities in flavour physics, and beyond, in the HL-LHC era. Technical report, CERN, Geneva (2021). https://cds.cern.ch/record/2776420

39. J.G. Contreras, J.D. Tapia Takaki, Ultra-peripheral heavy-ion collisions at the LHC. Int. J. Mod. Phys. A **30**, 1542012 (2015). https://doi.org/10.1142/S0217751X15420129

40. S.R. Klein, H. Mäntysaari, Imaging the nucleus with high-energy photons. Nat. Rev. Phys. **1**(11), 662–674 (2019). https://doi.org/10.1038/s42254-019-0107-6. arXiv:1910.10858 [hep-ex]

41. S. Acharya et al., Coherent photoproduction of $\rho^0$ vector mesons in ultra-peripheral Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. JHEP **06**, 035 (2020). https://doi.org/10.1007/JHEP06(2020)035. arXiv:2002.10897 [nucl-ex]

42. S. Acharya et al., First measurement of coherent $\rho^0$ photoproduction in ultra-peripheral Xe-Xe collisions at $s_{NN}$ = 5.44 TeV. Phys. Lett. B **820**, 136481 (2021). https://doi.org/10.1016/j.physletb.2021.136481. arXiv:2101.02581 [nucl-ex]

43. M. Abdallah et al., Tomography of ultrarelativistic nuclei with polarized photon-gluon collisions. Sci. Adv. **9**(1), 3903 (2023). https://doi.org/10.1126/sciadv.abq3903. arXiv:2204.01625 [nucl-ex]

44. S. Acharya et al., Coherent $J/\psi$ and $\psi'$ photoproduction at midrapidity in ultra-peripheral Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Eur. Phys. J. C **81**(8), 712 (2021). https://doi.org/10.1140/epjc/s10052-021-09437-6. arXiv:2101.04577 [nucl-ex]

45. R. Aaij et al., Study of coherent $J/\psi$ production in lead-lead collisions at $\sqrt{s_{NN}}$ = 5 TeV. JHEP **07**, 117 (2022). https://doi.org/10.1007/JHEP07(2022)117. arXiv:2107.03223 [hep-ex]

46. R. Aaij et al., Study of exclusive photoproduction of charmonium in ultra-peripheral lead-lead collisions. JHEP **06**, 146 (2023). https://doi.org/10.1007/JHEP06(2023)146. arXiv:2206.08221 [hep-ex]

47. Photon-nucleus energy dependence of coherent $J/\psi$ cross session in ultraperipheral PbPb collisions at 5.02 TeV with CMS. Technical report, CERN, Geneva (2022). https://cds.cern.ch/record/2843162

48. V. Guzey, M. Strikman, M. Zhalov, Disentangling coherent and incoherent quasielastic $J/\psi$ photoproduction on nuclei by
neutron tagging in ultraperipheral ion collisions at the LHC. Eur. Phys. J. C **74**(7), 2942 (2014). https://doi.org/10.1140/epjc/s10052-014-2942-z. arXiv:1312.6486 [hep-ph]

49. S. Acharya et al., Energy dependence of coherent photonuclear production of $J/\psi$ mesons in ultra-peripheral Pb-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV. JHEP **10**, 119 (2023). https://doi.org/10.1007/JHEP10(2023)119. arXiv:2305.19060 [nucl-ex]

50. A. Tumasyan et al., Probing small Bjorken-x nuclear gluonic structure via coherent $J/\psi$ photoproduction in ultraperipheral PbPb collisions at $\sqrt{s_{NN}} = 5.02$ TeV (2023). arXiv:2303.16984 [nucl-ex]

51. J.G. Contreras, Gluon shadowing at small $x$ from coherent $J/\psi$ photoproduction data at energies available at the CERN Large Hadron Collider. Phys. Rev. C **96**(1), 015203 (2017). https://doi.org/10.1103/PhysRevC.96.015203. arXiv:1610.03350 [nucl-ex]

52. S.R. Klein, J. Nystrand, J. Seger, Y. Gorbunov, J. Butterworth, STARLight: a Monte Carlo simulation program for ultraperipheral collisions of relativistic ions. Comput. Phys. Commun. **212**, 258-268 (2017). https://doi.org/10.1016/j.cpc.2016.10.016. arXiv:1607.03838 [hep-ph]

53. A. Accardi et al., Electron ion collider: the next QCD Frontier: understanding the glue that binds us all. Eur. Phys. J. A **52**(9), 268 (2016). https://doi.org/10.1140/epja/i2016-16268-9. arXiv:1212.1701 [nucl-ex]

54. J. Cepila, J.G. Contreras, J.D. Tapia Takaki, Energy dependence of dissociative $J/\psi$ photoproduction as a signature of gluon saturation at the LHC. Phys. Lett. B **766**, 186-191 (2017). https://doi.org/10.1016/j.physletb.2016.12.063. arXiv:1608.07559 [hep-ph]

55. Exclusive and dissociative $J/\psi$ photoproduction, and exclusive dimuon production, in p-Pb collisions at $\sqrt{s_{NN}} = 8.16$ TeV (2023). arXiv:2304.12403 [nucl-ex]

56. S. Acharya et al., First measurement of the $l^l$-dependence of coherent $J/\psi$ photonuclear production. Phys. Lett. B **817**, 136280 (2021). https://doi.org/10.1016/j.physletb.2021.136280. arXiv:2101.04623 [nucl-ex]

57. C.A. Flett, S.P. Jones, A.D. Martin, M.G. Ryskin, T. Teubner, How to include exclusive $J/\psi$ production data in global PDF analyses. Phys. Rev. D **101**(9), 094011 (2020). https://doi.org/10.1103/PhysRevD.101.094011. arXiv:1908.08398 [hep-ph]

58. R. Aaij et al., Updated measurements of exclusive $J/\psi$ and $\psi(2S)$ production cross-sections in pp collisions at $\sqrt{s} = 7$ TeV. J. Phys. G **41**, 055002 (2014). https://doi.org/10.1088/0954-3899/41/5/055002. arXiv:1401.3288 [hep-ex]

59. R. Aaij et al., Central exclusive production of $J/\psi$ and $\psi(2S)$ mesons in $pp$ collisions at $\sqrt{s} = 13$ TeV. JHEP **10**, 167 (2018). https://doi.org/10.1007/JHEP10(2018)167. arXiv:1806.04079 [hep-ex]

60. B.B. Abelev et al., Exclusive $J/\psi$ photoproduction off protons in ultra-peripheral p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV. Phys. Rev. Lett. **113**(23), 232504 (2014). https://doi.org/10.1103/PhysRevLett.113.232504. arXiv:1406.7819 [nucl-ex]

61. S. Acharya, Energy dependence of exclusive $J/\psi$ photoproduction off protons in ultra-peripheral p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV. Eur. Phys. J. C **79**(5), 402 (2019). https://doi.org/10.1140/epjc/s10052-019-6816-2. arXiv:1809.03235 [nucl-ex]

62. T. Herman, $J/\psi$ photoproduction results from ALICE. Acta Phys. Polon. Suppl. **16**(1), 98 (2023). https://doi.org/10.5506/APhysPolBSupp.16.1-A98

63. A. Aktas et al., Elastic $J/\psi$ production at HERA. Eur. Phys. J. C **46**, 585-603 (2006). https://doi.org/10.1140/epjc/s2006-02519-5. arXiv:hep-ex/0510016

64. C. Alexa et al., Elastic and proton-dissociative photoproduction of $J/\psi$ mesons at HERA. Eur. Phys. J. C **73**(6), 2466 (2013). https://doi.org/10.1140/epjc/s10052-013-2466-y. arXiv:1304.5162 [hep-ex]

65. S. Chekanov et al., Exclusive photoproduction of $J/\psi$ mesons at HERA. Eur. Phys. J. C **24**, 345-360 (2002). https://doi.org/10.1007/s10052-002-0953-7. arXiv:hep-ex/0201043

66. R. Aaij et al., Measurement of the exclusive $Υ$ production cross-section in $pp$ collisions at $\sqrt{s} = 7$ TeV and 8 TeV. JHEP **09**, 084 (2015). https://doi.org/10.1007/JHEP09(2015)084. arXiv:1505.08139 [hep-ex]

67. A.M. Sirunyan, et al., Measurement of exclusive $Υ$ photoproduction from protons in $pPb$ collisions at $\sqrt{s_{NN}} = 5.02$ TeV. Eur. Phys. J. C **79**(3), 277 (2019). https://doi.org/10.1140/epjc/s10052-019-6774-8. arXiv:1809.11080 [hep-ex] [Erratum: Eur. Phys. J. C 82, 343 (2022)]

68. A.M. Sirunyan, Measurement of exclusive $\rho(770)^0$ photoproduction in ultraperipheral $pPb$ collisions at $\sqrt{s_{NN}} = 5.02$ TeV. Eur. Phys. J. C **79**(8), 702 (2019). https://doi.org/10.1140/epjc/s10052-019-7202-9. arXiv:1902.01339 [hep-ex]

69. D.Y. Ivanov, A. Schafer, L. Szymanowski, G. Krasnikov, Exclusive photoproduction of a heavy vector meson in QCD. Eur. Phys. J. C **34**(3), 297-316 (2004) [Erratum: Eur. Phys. J. C 75, 75 (2015)]. https://doi.org/10.1140/epjc/s2004-01712-x. arXiv:hep-ph/0401131

70. K.J. Eskola, C.A. Flett, V. Guzey, T. Löytäinen, H. Paukkunen, Next-to-leading order perturbative QCD predictions for exclusive $J/\psi$ photoproduction in oxygen-oxygen and lead-lead collisions at energies available at the CERN Large Hadron Collider. Phys. Rev. C **107**(4), 044912 (2023). https://doi.org/10.1103/PhysRevC.107.044912. arXiv:2210.16048 [hep-ph]

71. K.J. Eskola, C.A. Flett, V. Guzey, T. Löytäinen, H. Paukkunen, Exclusive $J/\psi$ photoproduction in ultraperipheral Pb+Pb collisions at the CERN Large Hadron Collider calculated at next-to-leading order perturbative QCD. Phys. Rev. C **106**(3), 035202 (2022). https://doi.org/10.1103/PhysRevC.106.035202. arXiv:2203.11613 [hep-ph]

72. K.J. Eskola, C.A. Flett, V. Guzey, T. Löytäinen, H. Paukkunen, Predictions for exclusive $Υ$ photoproduction in ultraperipheral Pb + Pb collisions at the LHC at next-to-leading order in perturbative QCD. Eur. Phys. J. C **83**(8), 758 (2023). https://doi.org/10.1140/epjc/s10052-023-11927-8. arXiv:2303.03007 [hep-ph]

73. A.G. Shuvaev, K.J. Golec-Biernat, A.D. Martin, M.G. Ryskin, Off diagonal distributions fixed by diagonal partons at small $x$ and $\chi$. Phys. Rev. D **60**, 014015 (1999). https://doi.org/10.1103/PhysRevD.60.014015. arXiv:hep-ph/9902410

74. A. Shuvaev, Solution of the off forward leading logarithmic evolution equation based on the Gegenbauer moments inversion. Phys. Rev. D **60**, 116005 (1999). https://doi.org/10.1103/PhysRevD.60.116005. arXiv:hep-ph/9902318

75. K.J. Golec-Biernat, A.D. Martin, M.G. Ryskin, Diagonal input for the evolution of off diagonal partons. Phys. Lett. B **456**, 232-239 (1999). https://doi.org/10.1016/S0370-2693(99)00504-3. arXiv:hep-ph/9903327

76. S.P. Jones, A.D. Martin, M.G. Ryskin, T. Teubner, The exclusive $J/\psi$ process at the LHC tamed to probe the low $x$ gluon. Eur. Phys. J. C **76**(11), 633 (2016). https://doi.org/10.1140/epjc/s10052-016-4493-y. arXiv:1610.02272 [hep-ph]

77. S.P. Jones, A.D. Martin, M.G. Ryskin, T. Teubner, Exclusive $J/\psi$ and $Υ$ photoproduction and the low $x$ gluon. J. Phys. G **43**(3), 035002 (2016). https://doi.org/10.1088/0954-3899/43/3/035002. arXiv:1507.06942 [hep-ph]

78. C.A. Flett, A.D. Martin, M.G. Ryskin, T. Teubner, Very low $x$ gluon density determined by LHCb exclusive $J/\psi$ data. Phys. Rev. D **102**, 114021 (2020). https://doi.org/10.1103/PhysRevD.102.114021. arXiv:2006.13857 [hep-ph]

79. L. Apolinário, Y.-J. Lee, M. Winn, Heavy quarks and jets as probes of the QGP. Prog. Part. Nucl. Phys. **127**, 103990 (2022). https://doi.org/10.1016/j.ppnp.2022.103990. arXiv:2203.16352 [hep-ph]
80. J. Casalderrey-Solana, Y. Mehtar-Tani, C.A. Salgado, K. Tywoniuk, New picture of jet quenching dictated by color coherence. Phys. Lett. B 725, 357-360 (2013). https://doi.org/10.1016/j.physletb.2013.07.046. arXiv:1210.7765 [hep-ph]

81. ALICE Collaboration, Measurement of the groomed jet radius and momentum splitting fraction in pp and Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Phys. Rev. Lett. 128(10), 102001 (2022). https://doi.org/10.1103/PhysRevLett.128.102001. arXiv:2107.12984 [nucl-ex]

82. ATLAS Collaboration, Measurement of substructure-dependent jet suppression in Pb+Pb collisions at 5.02 TeV with the ATLAS detector. Phys. Rev. C 107(5), 054909 (2023). https://doi.org/10.1103/PhysRevC.107.054909. arXiv:2211.11470 [nucl-ex]

83. ATLAS Collaboration, Measurement of suppression of large-radius jets and its dependence on substructure in Pb+Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV with the ATLAS detector (2023). arXiv:2301.05606 [nucl-ex]

84. L. Cunqueiro, A.M. Sickles, Studying the QGP with Jets at the LHC and RHIC. Prog. Part. Nucl. Phys. 124, 103940 (2022). https://doi.org/10.1016/j.ppnp.2022.103940. arXiv:2110.14490 [nucl-ex]

85. D. Krohn, J. Thaler, L.-T. Wang, Jet trimming. JHEP 02, 084 (2010). https://doi.org/10.1007/JHEP02(2010)084. arXiv:0912.1342 [hep-ph]

86. ATLAS Collaboration, Measurement of angular and momentum distributions of charged particles within and around jets in Pb+Pb and $pp$ collisions at $\sqrt{s_{NN}}$ = 5.02 TeV with the ATLAS detector. Phys. Rev. C 100(6), 064901 (2019). https://doi.org/10.1103/PhysRevC.100.064901. arXiv:1908.05264 [nucl-ex] [Erratum: Phys. Rev. C 101, 059903 (2020)]

87. J. Brewer, Q. Brodsky, K. Rajagopal, Disentangling jet modification in jet simulations and in $Z^+\to jet$ data. JHEP 02, 175 (2022). https://doi.org/10.1007/JHEP02(2022)175. arXiv:2110.13159 [hep-ph]

88. Y.-L. Du, D. Pablos, K. Tywoniuk, Deep learning jet modifications in heavy-ion collisions. JHEP 21, 206 (2020). https://doi.org/10.1007/JHEP03(2021)206. arXiv:2012.07797 [hep-ph]

89. Y. Mehtar-Tani, A. Soto-Ontoso, K. Tywoniuk, Dynamical grooming of QCD jets. Phys. Rev. D 101(3), 034004 (2020). https://doi.org/10.1103/PhysRevD.101.034004. arXiv:1911.00375 [hep-ph]

90. P. Caucal, A. Soto-Ontoso, A. Takacs, Dynamically groomed jet radius in heavy-ion collisions. Phys. Rev. D 105(11), 114046 (2022). https://doi.org/10.1103/PhysRevD.105.114046. arXiv:2111.14768 [hep-ph]

91. ALICE Collaboration, Measurements of the groomed jet radius and momentum splitting fraction with the soft drop and dynamical grooming algorithms in $pp$ collisions at $\sqrt{s} = 5.02$ TeV. JHEP 05, 244 (2023). https://doi.org/10.1007/JHEP05(2023)244. arXiv:2204.10246 [nucl-ex]

92. P. Caucal, A. Soto-Ontoso, A. Takacs, Dynamical grooming meets LHC data. JHEP 07, 020 (2021). https://doi.org/10.1007/JHEP07(2021)020. arXiv:2103.06566 [hep-ph]

93. P. Caucal, E. Iancu, A.H. Mueller, G. Soyez, Vacuum-like jet fragmentation in a dense QCD medium. Phys. Rev. Lett. 120, 232001 (2018). https://doi.org/10.1103/PhysRevLett.120.232001. arXiv:1801.09703 [hep-ph]

94. M. Spousta, B. Cole, Interpreting single jet measurements in Pb+Pb collisions at the LHC. Eur. Phys. J. C 76(2), 50 (2016). https://doi.org/10.1140/epjc/s10052-016-3896-0. arXiv:1504.05169 [hep-ph]

95. J.-W. Qiu, F. Ringer, N. Sato, P. Zurita, Factorization of jet cross sections in heavy-ion collisions. Phys. Rev. Lett. 122(25), 252301 (2019). https://doi.org/10.1103/PhysRevLett.122.252301. arXiv:1903.01993 [hep-ph]

96. D. Pablos, A. Soto-Ontoso, Pushing forward jet substructure measurements in heavy-ion collisions. Phys. Rev. D 107(9), 094003 (2023). https://doi.org/10.1103/PhysRevD.107.094003. arXiv:2210.07901 [hep-ph]

97. C.L. Basham, L.S. Brown, S.D. Ellis, S.T. Love, Energy correlations in electron-positron annihilation: testing QCD. Phys. Rev. Lett. 41, 1585 (1978). https://doi.org/10.1103/PhysRevLett.41.1585

98. D.M. Hofman, J. Maldacena, Conformal collider physics: energy and charge correlations. JHEP 05, 012 (2008). https://doi.org/10.1088/1126-6708/2008/05/012. arXiv:0803.1467 [hep-th]

99. P. Kravchuk, D. Simmons-Duffin, Light-ray operators in conformal field theory. JHEP 11, 102 (2018). https://doi.org/10.1007/JHEP11(2018)102. arXiv:1805.00098 [hep-th]

100. L.J. Dixon, I. Moult, H.X. Zhu, Collinear limit of the energy-energy correlator. Phys. Rev. D 100(1), 014009 (2019). https://doi.org/10.1103/PhysRevD.100.014009. arXiv:1905.01310 [hep-ph]

101. C. Andres, F. Dominguez, R. Kunnawalkam Elayavalli, J. Holguin, C. Marquet, I. Moult, Resolving the scales of the quark-gluon plasma with energy correlators (2022). arXiv:2209.11236 [hep-ph]

102. C. Andres, F. Dominguez, J. Holguin, C. Marquet, I. Moult, A coherent view of the quark-gluon plasma from energy correlators (2023). arXiv:2303.03413 [hep-ph]

103. B.B. Abelev et al., Multiparticle azimuthal correlations in $p$-Pb and Pb-Pb collisions at the CERN Large Hadron Collider. Phys. Rev. C 90(5), 054901 (2014). https://doi.org/10.1103/PhysRevC.90.054901. arXiv:1406.2474 [nucl-ex]

104. M. Aaboud et al., Measurement of long-range multiparticle azimuthal correlations with the subevent cumulant method in $pp$ and $p+Pb$ collisions with the ATLAS detector at the CERN Large Hadron Collider. Phys. Rev. C 97(2), 024904 (2018). https://doi.org/10.1103/PhysRevC.97.024904. arXiv:1708.03559 [hep-ex]

105. A.M. Sirunyan et al., Observation of correlated azimuthal anisotropy Fourier harmonics in $pp$ and $p+Pb$ collisions at the LHC. Phys. Rev. Lett. 120(9), 092301 (2018). https://doi.org/10.1103/PhysRevLett.120.092301. arXiv:1709.09189 [nucl-ex]

106. C. Aidala et al., Creation of quark-gluon plasma droplets with three distinct geometries. Nat. Phys. 15(3), 214-220 (2019). https://doi.org/10.1038/s41567-018-0360-0. arXiv:1805.02973 [nucl-ex]

107. B.I. Abelev et al., Enhanced strange baryon production in Au+Au collisions compared to $p+p$ at $s(NN)^{1/2} = 200$-GeV. Phys. Rev. C 77, 044908 (2008). https://doi.org/10.1103/PhysRevC.77.044908. arXiv:0705.2511 [nucl-ex]

108. J. Adam et al., Enhanced production of multi-strange hadrons in high-multiplicity proton-proton collisions. Nat. Phys. 13, 535-539 (2017). https://doi.org/10.1038/nphys4111. arXiv:1606.07424 [nucl-ex]

109. J.L. Nagle, W.A. Zajc, Small system collectivity in relativistic hadronic and nuclear collisions. Annu. Rev. Nucl. Part. Sci. 68, 211-235 (2018). https://doi.org/10.1146/annurev-nucl-101916-123209. arXiv:1801.03477 [nucl-ex]

110. G. Aad et al., Two-particle azimuthal correlations in photonuclear ultraperipheral Pb+Pb collisions at 5.02 TeV with ATLAS. Phys. Rev. C 104(1), 014903 (2021). https://doi.org/10.1103/PhysRevC.104.014903. arXiv:2101.10771 [nucl-ex]

111. K. Werner, Core-corona separation in ultra-relativistic heavy ion collisions. Phys. Rev. Lett. 98, 152301 (2007). https://doi.org/10.1103/PhysRevLett.98.152301. arXiv:0704.1270 [nucl-th]

112. Z.-W. Lin, C.M. Ko, B.-A. Li, B. Zhang, S. Pal, A Multi-phase transport model for relativistic heavy ion collisions. Phys. Rev. C 72, 064901 (2005). https://doi.org/10.1103/PhysRevC.72.064901. arXiv:nucl-th/0411110

113. W. Cassing, E.L. Bratkovskaya, Parton-hadron-string dynamics: an off-shell transport approach for relativistic energies.
Nucl. Phys. A **831**, 215-242 (2009). https://doi.org/10.1016/j.nuclphysa.2009.09.007. arXiv:0907.5331 [nucl-th]

114. Y. Kanakubo, Y. Tachibana, T. Hirano, Interplay between core and corona components in high-energy nuclear collisions. Phys. Rev. C **105**(2), 024905 (2022). https://doi.org/10.1103/PhysRevC.105.024905. arXiv:2108.07943 [nucl-th]

115. Y. Kanakubo, Y. Tachibana, T. Hirano, Nonequilibrium components in the region of very low transverse momentum in high-energy nuclear collisions. Phys. Rev. C **106**(5), 054908 (2022). https://doi.org/10.1103/PhysRevC.106.054908. arXiv:2207.13966 [nucl-th]

116. K. Werner, I. Karpenko, T. Pierog, M. Bleicher, K. Mikhailov, Event-by-event simulation of the three-dimensional hydrodynamical evolution from flux tube initial conditions in ultrarelativistic heavy ion collisions. Phys. Rev. C **82**, 044904 (2010). https://doi.org/10.1103/PhysRevC.82.044904. arXiv:1004.0805 [nucl-th]

117. K. Werner, B. Guiot, I. Karpenko, T. Pierog, Analysing radial flow features in p-Pb and p-p collisions at several TeV by studying identified particle production in EPOS3. Phys. Rev. C **89**(6), 064903 (2014). https://doi.org/10.1103/PhysRevC.89.064903. arXiv:1312.1233 [nucl-th]

118. K. Werner, On a deep connection between factorization and saturation: new insight into modeling high-energy proton-proton and nucleus-nucleus scattering in the EPOS4 framework (2023). arXiv:2301.12517 [hep-ph]

119. K. Werner, Core-corona procedure and microcanonical hadronization to understand strangeness enhancement in proton-proton and heavy ion collisions in the EPOS4 framework. arXiv:2306.10277 [hep-ph]

120. EPOS4: A Monte Carlo Tool for Simulating High-energy Scattering. https://klaus/pages.in2p3.fr/epos4/

121. K. Werner, Core-corona procedure and microcanonical hadronization to understand strangeness enhancement in proton-proton and heavy ion collisions in the EPOS4 framework (2023). arXiv:2306.10277 [hep-ph]

122. P. Bozek, I. Wyskiel, Directed flow in ultrarelativistic heavy-ion collisions. Phys. Rev. C **81**, 054902 (2010). https://doi.org/10.1103/PhysRevC.81.054902. arXiv:1002.4999 [nucl-th]

123. W. Zhao, S. Ryu, C. Shen, B. Schenke, 3D structure of anisotropic flow in small collision systems at energies available at the BNL Relativistic Heavy Ion Collider. Phys. Rev. C **107**(1), 014904 (2023). https://doi.org/10.1103/PhysRevC.107.014904. arXiv:2211.16376 [nucl-th]

124. L. Oliva, P. Moreau, V. Voronyuk, E. Bratkovskaya, Influence of electromagnetic fields in proton-nucleus collisions at relativistic energy. Phys. Rev. C **101**(1), 014917 (2020). https://doi.org/10.1103/PhysRevC.101.014917. arXiv:1909.06770 [nucl-th]

125. S. Chatterjee, P. Božek, Large directed flow of open charm mesons probes the three dimensional distribution of matter in heavy ion collisions. Phys. Rev. Lett. **120**(19), 192301 (2018). https://doi.org/10.1103/PhysRevLett.120.192301. arXiv:1712.01189 [nucl-th]

126. M.R. Haque, S. Singha, B. Mohanty, Probing the profile of bulk matter in p+Pb collisions via directed flow of heavy quarks. Phys. Rev. C **104**(2), 024901 (2021). https://doi.org/10.1103/PhysRevC.104.024901. arXiv:2107.07457 [hep-ph]

127. L. Oliva, S. Plumari, V. Greco, Directed flow of D mesons at RHIC and LHC: non-perturbative dynamics, longitudinal bulk matter asymmetry and electromagnetic fields. JHEP **05**, 034 (2021). https://doi.org/10.1007/JHEP05(2021)034. arXiv:2009.11066 [hep-ph]

128. J.E. Bernhard, J.S. Moreland, S.A. Bass, Bayesian estimation of the specific shear and bulk viscosity of quark-gluon plasma. Nat. Phys. **15**(11), 1113-1117 (2019). https://doi.org/10.1038/s41567-019-0611-8

129. D. Everett et al., Multisystem Bayesian constraints on the transport coefficients of QCD matter. Phys. Rev. C **103**(5), 054904 (2021). https://doi.org/10.1103/PhysRevC.103.054904. arXiv:2011.01430 [hep-ph]

130. G. Nijs, W. Schee, Hadronic nucleus-nucleus cross section and the nucleon size. Phys. Rev. Lett. **129**(23), 232301 (2022). https://doi.org/10.1103/PhysRevLett.129.232301. arXiv:2206.13522 [nucl-th]

131. A. Mazeliauskas, D. Teaney, Fluctuations of harmonic and radial flow in heavy ion collisions with principal components. Phys. Rev. C **93**(2), 024913 (2016). https://doi.org/10.1103/PhysRevC.93.024913. arXiv:1509.07492 [nucl-th]

132. P. Bozek, Transverse-momentum-flow correlations in relativistic heavy-ion collisions. Phys. Rev. C **93**(4), 044908 (2016). https://doi.org/10.1103/PhysRevC.93.044908. arXiv:1601.04513 [nucl-th]

133. S.H. Lim, J.L. Nagle, Exploring origins for correlations between flow harmonics and transverse momentum in small collision systems. Phys. Rev. C **103**(6), 064906 (2021). https://doi.org/10.1103/PhysRevC.103.064906. arXiv:2103.01348 [nucl-th]

134. G. Aad, et al., Measurement of flow harmonics correlations with mean transverse momentum in lead-lead and proton-lead collisions at $\sqrt{s_{NN}}$ = 5.02 TeV with the ATLAS detector. Eur. Phys. J. C **79**(12) (2019). https://doi.org/10.1140/epjc/s10052-019-7489-6

135. Correlations between multiparticle cumulants and mean transverse momentum in small collision systems with the CMS detector. Technical report, CERN, Geneva (2022). https://cds.cern.ch/record/2805932

136. A. Banfi, G.P. Salam, G. Zanderighi, Phenomenology of event shapes at hadron colliders. JHEP **06**, 038 (2010). https://doi.org/10.1007/JHEP06(2010)038. arXiv:1001.4082 [hep-ph]

137. A. Ortiz, G. Paić, E. Cuautle, Mid-rapidity charged hadron transverse spherocity in pp collisions simulated with Pythia. Nucl. Phys. A **941**, 78-86 (2015). https://doi.org/10.1016/j.nuclphysa.2015.05.010. arXiv:1503.03129 [hep-ph]

138. S. Acharya et al., Charged-particle production as a function of multiplicity and transverse spherocity in pp collisions at $\sqrt{s} = 5.02$ and 13 TeV. Eur. Phys. J. C **79**(10), 857 (2019). https://doi.org/10.1140/epjc/s10052-019-7350-y. arXiv:1905.07208 [nucl-ex]

139. N. Mallick, R. Sahoo, S. Tripathy, A. Ortiz, Study of transverse spherocity and azimuthal anisotropy in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV using a multi-phase transport model. J. Phys. G **48**(4), 045104 (2021). https://doi.org/10.1088/1361-6471/abeb59. arXiv:2008.13616 [hep-ph]

140. L. Oliva, W. Fan, P. Moreau, S.A. Bass, E. Bratkovskaya, Nonequilibrium effects and transverse spherocity in ultrarelativistic proton-nucleus collisions. Phys. Rev. C **106**(4), 044910 (2022). https://doi.org/10.1103/PhysRevC.106.044910. arXiv:2204.04194 [nucl-th]

141. T. Martin, P. Skands, S. Farrington, Probing collective effects in hadronisation with the extremes of the underlying event. Eur. Phys. J. C **76**(5), 299 (2016). https://doi.org/10.1140/epjc/s10052-016-4135-4. arXiv:1603.05298 [hep-ph]

142. G. Bencedi, A. Ortiz, A. Paz, Disentangling the hard gluon bremsstrahlung effects from the relative transverse activity classifier in pp collisions. Phys. Rev. D **104**(1), 016017 (2021). https://doi.org/10.1103/PhysRevD.104.016017. arXiv:2105.04838 [hep-ph]

143. A. Ortiz, G. Paic, A look into the “hedgehog” events in pp collisions. Rev. Mex. Fis. Suppl. **3**(4), 040911 (2022). https://doi.org/10.31349/SuplRevMexFis.3.040911. arXiv:2204.13733 [hep-ph]

144. A. Ortiz, A. Khuntia, O. Vázquez-Rueda, S. Tripathy, G. Bencedi, S. Prasad, F. Fan, Unveiling the effects of multiple soft partonic interactions in pp collisions at $s=13.6$ TeV using a new event
classifier. Phys. Rev. D **107**(7), 076012 (2023). https://doi.org/
10.1103/PhysRevD.107.076012. arXiv:2211.06093 [hep-ph]

145. A.M. Poskanzer, S.A. Voloshin, Methods for analyzing anisotropic flow in relativistic nuclear collisions. Phys. Rev. C **58**, 1671–1678 (1998). https://doi.org/10.1103/PhysRevC.58.1671. arXiv:nucl-ex/9805001

146. G. Aad et al., Measurement of the correlation between flow harmonics of different order in lead-lead collisions at $\sqrt{s_{NN}}$=2.76 TeV with the ATLAS detector. Phys. Rev. C **92**(3), 034903 (2015). https://doi.org/10.1103/PhysRevC.92.034903. arXiv:1504.01289 [hep-ex]

147. F. Becattini, Spin and polarization: a new direction in relativistic heavy ion physics. Rep. Prog. Phys. **85**(12), 122301 (2022). https://doi.org/10.1088/1361-6633/ac97a9. arXiv:2204.01144 [nucl-th]

148. W.M. Serenone, J.G.P. Barbon, D.D. Chinellato, M.A. Lisa, C. Shen, J. Takahashi, G. Torrieri, $\lambda$ polarization from thermalized jet energy. Phys. Lett. B **820**, 136500 (2021). https://doi.org/10.1016/j.physletb.2021.136500

149. V.H. Ribeiro, D.D. Chinellato, M.A. Lisa, W.M. Serenone, C. Shen, J. Takahashi, G. Torrieri, A polarization from vortex ring as medium response for jet thermalization (2023)

150. C. Bierlich et al., A comprehensive guide to the physics and usage of PYTHIA 8.3 (2022). https://doi.org/10.21468/SciPostPhysCodeb.8. arXiv:2203.11601 [hep-ph]

151. J. Bellm et al., Herwig 7.2 release note. Eur. Phys. J. C **80**(5), 452 (2020). https://doi.org/10.1140/epjc/s10052-020-8011-x. arXiv:1912.06509 [hep-ph]

152. J.R. Christiansen, P.Z. Skands, String formation beyond leading colour. JHEP **08**, 003 (2015). https://doi.org/10.1007/JHEP08(2015)003. arXiv:1505.01681 [hep-ph]

153. S. Gieseke, P. Kirchgaeßer, S. Plätzer, Baryon production from cluster hadronisation. Eur. Phys. J. C **78**(2), 99 (2018). https://doi.org/10.1140/epjc/s10052-018-5585-7. arXiv:1710.10906 [hep-ph]

154. C. Bierlich, G. Gustafson, L. Lönnblad, A. Tarasov, Effects of overlapping strings in pp collisions. JHEP **03**, 148 (2015). https://doi.org/10.1007/JHEP03(2015)148. arXiv:1412.6259 [hep-ph]

155. N. Fischer, T. Sjöstrand, Thermodynamical string fragmentation. JHEP **01**, 140 (2017). https://doi.org/10.1007/JHEP01(2017)140. arXiv:1610.09818 [hep-ph]

156. C.B. Duncan, P. Kirchgaeßer, Kinematic strangeness production in cluster hadronization. Eur. Phys. J. C **79**(1), 61 (2019). https://doi.org/10.1140/epjc/s10052-019-6573-2. arXiv:1811.10336 [hep-ph]

157. C. Bierlich, S. Chakraborty, G. Gustafson, L. Lönnblad, Strangeness enhancement across collision systems without a plasma. Phys. Lett. B **835**, 137571 (2022). https://doi.org/10.1016/j.physletb.2022.137571. arXiv:2205.11170 [hep-ph]

158. S. Schael et al., Electroweak measurements in electron-positron collisions at W-boson-pair energies at LEP. Phys. Rep. **532**, 119–244 (2013). https://doi.org/10.1016/j.physrep.2013.07.004. arXiv:1302.3415 [hep-ex]

159. A. Ortiz Velasquez, P. Christiansen, E. Cuautle Flores, I. Maldonado Cervantes, G. Païc, Color reconnection and flowlike patterns in *pp* collisions. Phys. Rev. Lett. **111**(4), 042001 (2013). https://doi.org/10.1103/PhysRevLett.111.042001. arXiv:1303.6326 [hep-ph]

160. C. Bierlich, G. Gustafson, L. Lönnblad, A shoving model for collectivity in hadronic collisions (2016). arXiv:1612.05132 [hep-ph]

161. C.B. Duncan, P. Skands, Fragmentation of two repelling lund strings. SciPost Phys. **8**, 080 (2020). https://doi.org/10.21468/SciPostPhys.8.5.080. arXiv:1912.09639 [hep-ph]

162. C. Bierlich, S. Chakraborty, G. Gustafson, L. Lönnblad, Setting the string shoving picture in a new frame. JHEP **03**, 270 (2021).

https://doi.org/10.1007/JHEP03(2021)270. arXiv:2010.07595 [hep-ph]

163. V. Greco, C.M. Ko, P. Levai, Parton coalescence at RHIC. Phys. Rev. C **68**, 034904 (2003). https://doi.org/10.1103/PhysRevC.68.034904. arXiv:nucl-th/0305024

164. R.J. Fries, V. Greco, P. Sorensen, Coalescence models for hadron formation from quark gluon plasma. Annu. Rev. Nucl. Part. Sci. **58**, 177–205 (2008). https://doi.org/10.1146/annurev.nucl.58.110707.171134. arXiv:0807.4939 [nucl-th]

165. L. Ravagli, R. Rapp, Quark coalescence based on a transport equation. Phys. Lett. B **655**, 126–131 (2007). https://doi.org/10.1016/j.physletb.2007.07.043. arXiv:0705.0021 [hep-ph]

166. V. Minissale, F. Scardina, V. Greco, Hadrons from coalescence plus fragmentation in AA collisions at energies available at the BNL Relativistic Heavy Ion Collider to the CERN Large Hadron Collider. Phys. Rev. C **92**(5), 054904 (2015). https://doi.org/10.1103/PhysRevC.92.054904. arXiv:1502.06213 [nucl-th]

167. Y. Oh, C.M. Ko, S.H. Lee, S. Yasui, Heavy baryon/meson ratios in relativistic heavy ion collisions. Phys. Rev. C **79**, 044905 (2009). https://doi.org/10.1103/PhysRevC.79.044905. arXiv:0901.1382 [nucl-th]

168. S. Plumari, V. Minissale, S.K. Das, G. Coci, V. Greco, Charmed Hadrons from coalescence plus fragmentation in relativistic nucleus-nucleus collisions at RHIC and LHC. Eur. Phys. J. C **78**(4), 348 (2018). https://doi.org/10.1140/epjc/s10052-018-5828-7. arXiv:1712.00730 [hep-ph]

169. V. Minissale, S. Plumari, V. Greco, Charm hadrons in pp collisions at LHC energy within a coalescence plus fragmentation approach. Phys. Lett. B **821**, 136622 (2021). https://doi.org/10.1016/j.physletb.2021.136622. arXiv:2012.12001 [hep-ph]

170. A. Andronic, P. Braun-Munzinger, K. Redlich, J. Stachel, Decoding the phase structure of QCD via particle production at high energy. Nature **561**(7723), 321–330 (2018). https://doi.org/10.1038/s41586-018-0491-6. arXiv:1710.09425 [nucl-th]

171. J. Cleymans, P.M. Lo, K. Redlich, N. Sharma, Multiplicity dependence of (multi)strange baryons in the canonical ensemble with phase shift corrections. Phys. Rev. C **103**(1), 014904 (2021). https://doi.org/10.1103/PhysRevC.103.014904. arXiv:2009.04844 [hep-ph]

172. V. Vovchenko, B. Dönigus, H. Stoecker, Canonical statistical model analysis of p-p, p-Pb, and Pb-Pb collisions at energies available at the CERN Large Hadron Collider. Phys. Rev. C **100**(5), 054906 (2019). https://doi.org/10.1103/PhysRevC.100.054906. arXiv:1906.03145 [hep-ph]

173. A. Beraudo, A. De Pace, D. Pablos, F. Prino, M. Monteno, M. Nardi, Heavy-flavor transport and hadronization in pp collisions (2023). arXiv:2306.02152 [hep-ph]

174. R. Aaij et al., Observation of multiplicity dependent prompt $\chi_{c1}$ (3872) and $\psi(2S)$ production in *pp* collisions. Phys. Rev. Lett. **126**(9), 092001 (2021). https://doi.org/10.1103/PhysRevLett.126.092001. arXiv:2009.06619 [hep-ex]

175. K. Ackerstaff et al., Production of K0(S) and Lambda in quark and gluon jets from Z0 decay. Eur. Phys. J. C **8**, 241–254 (1999). https://doi.org/10.1007/s100529901058. arXiv:hep-ex/9805025

176. N. Hunt-Smith, P. Skands, String fragmentation with a time-dependent tension. Eur. Phys. J. C **80**(11), 1073 (2020). https://doi.org/10.1140/epjc/s10052-020-08654-9. arXiv:2005.06219 [hep-ph]

177. P. Abreu et al., Tuning and test of fragmentation models based on identified particles and precision event shape data. Z. Phys. C **73**, 11–60 (1996). https://doi.org/10.1007/s002880050295

178. J. Adam et al., Insight into particle production mechanisms via angular correlations of identified particles in pp collisions at $\sqrt{s} = 7$ TeV. Eur. Phys. J. C **77**(8), 569 (2017). https://doi.org/10.1140/epjc/s10052-017-5129-6. arXiv:1612.08975 [nucl-ex]
179. P.A. Zyla et al., Review of particle physics. PTEP 2020(8), 083-01 (2020). https://doi.org/10.1093/ptep/ptaa104

180. T. Sjöstrand, M. Utheim, A framework for hadronic rescattering in pp collisions. Eur. Phys. J. C 80(10), 907 (2020). https://doi.org/10.1140/epjc/s10052-020-8399-3. arXiv:2005.05658 [hep-ph]

181. S. Acharya et al., Production of $K_S^0$, $\Lambda(\bar{\Lambda})$, $\Xi^{\pm}$ and $\Omega^{\pm}$ in jets and in the underlying event in pp and p-Pb collisions. JHEP 07, 136 (2023). https://doi.org/10.1007/JHEP07(2023)136. arXiv:2211.08936 [nucl-ex]

182. S. Acharya et al., Enhanced deuteron coalescence probability in jets. Phys. Rev. Lett. 131(4), 042301 (2023). https://doi.org/10.1103/PhysRevLett.131.042301. arXiv:2211.15204 [nucl-ex]

183. S. Acharya, et al., Studying strangeness and baryon production mechanisms through angular correlations between charged $\Xi$ baryons and identified hadrons in pp collisions at $\sqrt{s} = 13$ TeV (2023). arXiv:2308.16706 [hep-ex]

184. S. Acharya et al., Charm-quark fragmentation fractions and production cross section at midrapidity in pp collisions at the LHC. Phys. Rev. D 105(1), 011103 (2022). https://doi.org/10.1103/PhysRevD.105.L011103. arXiv:2105.06335 [nucl-ex]

185. M. He, R. Rapp, Charm-baryon production in proton-proton collisions. Phys. Lett. B 795, 117-121 (2019). https://doi.org/10.1016/j.physletb.2019.06.004. arXiv:1902.08889 [nucl-th]

186. D. Ebert, R.N. Faustov, V.O. Galkin, Spectroscopy and Regge trajectories of heavy baryons in the relativistic quark-diquark picture. Phys. Rev. D 84, 014025 (2011). https://doi.org/10.1103/PhysRevD.84.014025. arXiv:1105.0583 [hep-ph]

187. M. He, R. Rapp, Bottom hadrochemistry in high-energy hadronic collisions. Phys. Rev. Lett. 131(1), 012301 (2023). https://doi.org/10.1103/PhysRevLett.131.012301. arXiv:2209.13419 [hep-ph]

188. A. Andronic, P. Braun-Munzinger, M.K. Köhler, A. Mazeliauskas, K. Redlich, J. Stachel, V. Vislavicius, The multiplet charm hierarchy in the statistical hadronization model. JHEP 07, 035 (2021). https://doi.org/10.1007/JHEP07(2021)035. arXiv:2104.12754 [hep-ph]

189. A. Beraudo et al., Extraction of heavy-flavor transport coefficients in QCD matter. Nucl. Phys. A 979, 21-86 (2018). https://doi.org/10.1016/j.nuclphysa.2018.09.002. arXiv:1803.03824 [nucl-th]

190. N. Brambilla, S. Eidelman, C. Hanhart, A. Nefediev, C.-P. Shen, C.E. Thomas, A. Vairo, C.-Z. Yuan, The XYZ states: experimental and theoretical status and perspectives. Phys. Rep. 873, 1-154 (2020). https://doi.org/10.1016/j.physrep.2020.05.001. arXiv:1907.07583 [hep-ex]

191. R. Aaij et al., Observation of an exotic narrow doubly charmed tetraquark. Nat. Phys. 18(7), 751-754 (2022). https://doi.org/10.1038/s41567-022-01614-y. arXiv:2109.01038 [hep-ex]

192. R. Aaij et al., Study of the doubly charmed tetraquark $T_{cc}^+$. Nat. Commun. 13(1), 3351 (2022). https://doi.org/10.1038/s41467-022-30206-w. arXiv:2109.01056 [hep-ex]

193. A. Esposito, E.G. Ferreiro, A. Pilloni, A.D. Polosa, C.A. Salgado, The nature of X(3872) from high-multiplicity pp collisions. Eur. Phys. J. C 81(7), 669 (2021). https://doi.org/10.1140/epjc/s10052-021-09425-w. arXiv:2006.15044 [hep-ph]

194. B. Wu, X. Du, M. Sibila, R. Rapp, X(3872) transport in heavy-ion collisions. Eur. Phys. J. A 57(4), 122 (2021). https://doi.org/10.1140/epja/s10050-021-00623-4 arXiv:2006.09945 [nucl-th] [Erratum: Eur. Phys. J. A 57, 314 (2021)]

195. S.H. Lee, S. Noh, D. Park, H.-O. Yoon, Exotic particles and heavy ion collisions. EPJ Web Conf. 276, 06015 (2023). https://doi.org/10.1051/epjconf/202327606015

196. Modification of $\chi_{c1}(3872)$ and $\psi(2S)$ production in $pPb$ collisions at $\sqrt{s_{NN}} = 8.16$ TeV. Technical report, CERN, Geneva (2022). https://cds.cern.ch/record/2807146

197. R. Aaij et al., Observation of structure in the $J/\psi$-pair mass spectrum. Sci. Bull. 65(23), 1983-1993 (2020). https://doi.org/10.1016/j.scib.2020.08.032. arXiv:2006.16957 [hep-ex]

198. G. Aad et al., Observation of an excess of di-charmonium events in the four-muon final state with the ATLAS detector (2023). arXiv:2304.08962 [hep-ex]

199. A. Hayrapetyan et al., Observation of new structure in the $J/\psi J/\psi$ mass spectrum in proton-proton collisions at $\sqrt{s} = 13$ TeV (2023). arXiv:2306.07164 [hep-ex]

200. A.M. Sirunyan et al., Studies of charm and beauty hadron long-range correlations in pp and $pPb$ collisions at LHC energies. Phys. Lett. B 813, 136036 (2021). https://doi.org/10.1016/j.physletb.2020.136036. arXiv:2009.07065 [hep-ex]

201. G. Aad et al., Measurement of azimuthal anisotropy of muons from charm and bottom hadrons in $pp$ collisions at $\sqrt{s} = 13$ TeV with the ATLAS detector. Phys. Rev. Lett. 124(8), 082301 (2020). https://doi.org/10.1103/PhysRevLett.124.082301. arXiv:1909.01650 [nucl-ex]

202. Azimuthal anisotropy of jet particles in $p-Pb$ and $Pb-Pb$ collisions at $\sqrt{s_{NN}} = 5.02$ TeV (2022). arXiv:2212.12609 [nucl-ex]

203. G. Aad et al., Measurement of the sensitivity of two-particle correlations in $pp$ collisions to the presence of hard scatterings (2023). arXiv:2303.17357 [nucl-ex]

204. S. Acharya et al., Measurement of inclusive charged-particle b-jet production in pp and $p-Pb$ collisions at $\sqrt{s_{NN}} = 5.02$ TeV. JHEP 01, 178 (2022). https://doi.org/10.1007/JHEP01(2022)178. arXiv:2110.06104 [nucl-ex]

205. V. Khachatryan et al., Measurement of inclusive jet production and nuclear modifications in $pPb$ collisions at $\sqrt{s_{NN}} = 5.02$ TeV. Eur. Phys. J. C 76(7), 372 (2016). https://doi.org/10.1140/epjc/s10052-016-4205-7. arXiv:1601.02001 [nucl-ex]

206. G. Aad et al., Centrality and rapidity dependence of inclusive jet production in $\sqrt{s_{NN}} = 5.02$ TeV proton-lead collisions with the ATLAS detector. Phys. Lett. B 748, 392-413 (2015). https://doi.org/10.1016/j.physletb.2015.07.023. arXiv:1412.4092 [hep-ex]

207. Measurement of the prompt $D^0$ nuclear modification factor in $pPb$ collisions at $\sqrt{s_{NN}} = 8.16$ TeV (2022). arXiv:2205.03936 [nucl-ex]

208. C. Zhang, C. Marquet, G.-Y. Qin, S.-Y. Wei, B.-W. Xiao, Elliptic flow of heavy quarkonia in $pA$ collisions. Phys. Rev. Lett. 122(17), 172302 (2019). https://doi.org/10.1103/PhysRevLett.122.172302. arXiv:1901.10320 [hep-ph]

209. Azimuthal anisotropy of $Y(1S)$ mesons in $pPb$ collisions at $\sqrt{s_{NN}} = 8.16$ TeV (2022)

210. A. Ipp, D.I. Müller, D. Schuh, Anisotropic momentum broadening in the 2+1D Glasma: analytic weak field approximation and lattice simulations. Phys. Rev. D 102(7), 074001 (2020). https://doi.org/10.1103/PhysRevD.102.074001. arXiv:2001.10001 [hep-ph]

211. J.H. Liu, S.K. Das, V. Greco, M. Ruggieri, Ballistic diffusion of heavy quarks in the early stage of relativistic heavy ion collisions at RHIC and the LHC. Phys. Rev. D 103(3) (2021). https://doi.org/10.1103/physrevd.103.034029

212. A. Collaboration, ALICE physics projections for a short oxygen-beam run at the LHC, ALICE-PUBLIC-2021-004 (2021)

213. Study of charged particle production at high $pT$ using event topology in pp, $pPb$ and $PbPb$ collisions at $\sqrt{s_{NN}} = 5.02$ TeV. Phys. Lett. B 843, 137649 (2023). https://doi.org/10.1016/j.physletb.2022.137649. arXiv:2204.10157 [nucl-ex]

214. F. Krizek, Search for jet quenching in high-multiplicity $pp$ collisions using inclusive and semi-inclusive jet production in ALICE. Acta Phys. Polon. Supp. 16(1), 69 (2023). https://doi.org/10.5506/APhysPolBSupp.16.1-A69

215. K. Boguslavski, A. Kurkela, T. Lappi, F. Lindenbauer, J. Peuron, Jet momentum broadening during initial stages in heavy-ion collisions (2023). arXiv:2303.12595 [hep-ph]

216. C. Andres, L. Apolinário, F. Dominguez, M.G. Martinez, C.A. Salgado, Medium-induced radiation with vacuum propagation in the pre-hydrodynamics phase. JHEP 03, 189 (2023). https://doi.org/10.1007/JHEP03(2023)189. arXiv:2211.10161 [hep-ph]
217. I. Grishmanovskii, T. Song, O. Soloveva, C. Greiner, E. Bratkovskaya, Exploring jet transport coefficients by elastic scattering in the strongly interacting quark-gluon plasma. Phys. Rev. C 106 (1), 014903 (2022). https://doi.org/10.1103/PhysRevC.106.014903. arXiv:2204.01561 [nucl-th]

218. T. Song, I. Grishmanovskii, O. Soloveva, Soft gluon emission from heavy quark scattering in strongly interacting quark-gluon plasma. Phys. Rev. D 107 (3), 036009 (2023). https://doi.org/10.1103/PhysRevD.107.036009. arXiv:2210.04010 [nucl-th]

219. I. Grishmanovskii, O. Soloveva, T. Song, C. Greiner, E. Bratkovskaya, Inelastic and elastic parton scatterings in the strongly interacting quark-gluon plasma (2023). arXiv:2308.03105 [hep-ph]

220. P. Moreau, O. Soloveva, L. Oliva, T. Song, W. Cassing, E. Bratkovskaya, Exploring the partonic phase at finite chemical potential within an extended off-shell transport approach. Phys. Rev. C 100 (1), 014911 (2019). https://doi.org/10.1103/PhysRevC.100.014911. arXiv:1903.10257 [nucl-th]

221. P.B. Arnold, W. Xiao, High-energy jet quenching in weakly-coupled quark-gluon plasmas. Phys. Rev. D 78, 125008 (2008). https://doi.org/10.1103/PhysRevD.78.125008. arXiv:0810.1026 [hep-ph]

222. S. Caron-Huot, O(g) plasma effects in jet quenching. Phys. Rev. D 79, 065039 (2009). https://doi.org/10.1103/PhysRevD.79.065039. arXiv:0811.1603 [hep-ph]

223. A. Tumasyan et al., Fragmentation of jets containing a prompt $J/\psi$ meson in PbPb and pp collisions at $\sqrt{s_{NN}} = 5.02$ TeV. Phys. Lett. B 825, 136842 (2022). https://doi.org/10.1016/j.physletb.2021.136842. arXiv:2106.13235 [hep-ex]

224. S. Acharya et al., Prompt $D^0$, $D^+$, and $D^{*+}$ production in Pb-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV. JHEP 01, 174 (2022). https://doi.org/10.1007/JHEP01(2022)174. arXiv:2110.09420 [nucl-ex]

225. T. Song, H. Berrehrah, D. Cabrera, J.M. Torres-Rincon, L. Tolos, W. Cassing, E. Bratkovskaya, Tomography of the quark-gluon-plasma by charm quarks. Phys. Rev. C 92(1), 014910 (2015). https://doi.org/10.1103/PhysRevC.92.014910. arXiv:1503.03039 [nucl-th]

226. S. Plumari, G. Coci, V. Minissale, S.K. Das, Y. Sun, V. Greco, Heavy-light flavor correlations of anisotropic flows at LHC energies within event-by-event transport approach. Phys. Lett. B 805, 135460 (2020). https://doi.org/10.1016/j.physletb.2020.135460. arXiv:1912.09350 [hep-ph]

227. M.L. Sambataro, Y. Sun, V. Minissale, S. Plumari, V. Greco, Event-shape engineering analysis of D meson in ultrarelativistic heavy-ion collisions. Eur. Phys. J. C 82(9), 833 (2022). https://doi.org/10.1140/epjc/s10052-022-10802-2. arXiv:2206.03160 [hep-ph]

228. F. Prino, R. Rapp, Open heavy flavor in QCD matter and in nuclear collisions. J. Phys. G Nucl. Part. Phys. 43(9), 093002 (2016). https://doi.org/10.1088/0954-3899/43/9/093002

229. J. Zhao, K. Zhou, S. Chen, P. Zhuang, Heavy flavors under extreme conditions in high energy nuclear collisions (2020)

230. M.L. Sambataro, V. Minissale, S. Plumari, V. Greco, B meson production in Pb+Pb at 5.02 ATeV at LHC: estimating the diffusion coefficient in the infinite mass limit. Phys. Lett. B 849, 138480 (2024). https://doi.org/10.1016/j.physletb.2024.138480. arXiv:2304.02953 [hep-ph]

231. M. Nahrgang, J. Aichelin, P.B. Gossiaux, K. Werner, Azimuthal correlations of heavy quarks in Pb + Pb collisions at $\sqrt{s} = 2.76$ TeV at the CERN Large Hadron Collider. Phys. Rev. C 90(2), 024907 (2014). https://doi.org/10.1103/PhysRevC.90.024907. arXiv:1305.3823 [hep-ph]

232. F. Özel, P. Freire, Masses, radii, and the equation of state of neutron stars. Annu. Rev. Astron. Astrophys. 54, 401–440 (2016). https://doi.org/10.1146/annurev-astro-081915-023322. arXiv:1603.02698 [astro-ph.HE]

233. T.E. Riley et al., A NICER view of PSR J0030+0451: millisecond pulsar parameter estimation. Astrophys. J. Lett. 887(1), 21 (2019). https://doi.org/10.3847/2041-8213/ab481c. arXiv:1912.05702 [astro-ph.HE]

234. P. Demorest, T. Pennucci, S. Ransom, M. Roberts, J. Hesels, Shapiro delay measurement of a two solar mass neutron star. Nature 467, 1081–1083 (2010). https://doi.org/10.1038/nature09466. arXiv:1010.5788 [astro-ph.HE]

235. J. Antoniadis et al., A massive pulsar in a compact relativistic binary. Science 340, 6131 (2013). https://doi.org/10.1126/science.1233232. arXiv:1304.6875 [astro-ph.HE]

236. H.T. Cromartie et al., Relativistic Shapiro delay measurements of an extremely massive millisecond pulsar. Nat. Astron. 4(1), 72–76 (2019). https://doi.org/10.1038/s41550-019-0880-2. arXiv:1904.06759 [astro-ph.HE]

237. B.P. Abbott et al., GW170817: observation of gravitational waves from a binary neutron star inspiral. Phys. Rev. Lett. 119(16), 161101 (2017). https://doi.org/10.1103/PhysRevLett.119.161101. arXiv:1710.05832 [gr-qc]

238. J.M. Lattimer, Neutron star mass and radius measurements. Universe 5(7), 159 (2019). https://doi.org/10.3390/universe5070159

239. J.M. Lattimer, Neutron stars and the nuclear matter equation of state. Annu. Rev. Nucl. Part. Sci. 71, 433–464 (2021). https://doi.org/10.1146/annurev-nucl-102419-124827

240. H. Djapo, B.-J. Schaefer, J. Wambach, On the appearance of hyperons in neutron stars. Phys. Rev. C 81, 035803 (2010). https://doi.org/10.1103/PhysRevC.81.035803. arXiv:0811.2939 [nucl-th]

241. L. Tolos, L. Fabbietti, Strangeness in nuclei and neutron stars. Prog. Part. Nucl. Phys. 112, 103770 (2020). https://doi.org/10.1016/j.ppnp.2020.103770. arXiv:2002.09223 [nucl-ex]

242. S. Acharya et al., p-p, p-Λ and Λ-Λ correlations studied via femtoscopy in pp reactions at $\sqrt{s} = 7$ TeV. Phys. Rev. C 99(2), 024001 (2019). https://doi.org/10.1103/PhysRevC.99.024001. arXiv:1805.12455 [nucl-ex]

243. S. Acharya et al., First observation of an attractive interaction between a proton and a cascade baryon. Phys. Rev. Lett. 123(11), 112002 (2019). https://doi.org/10.1103/PhysRevLett.123.112002. arXiv:1904.12198 [nucl-ex]

244. A. Collaboration et al., Unveiling the strong interaction among hadrons at the LHC. Nature 588, 232–238 (2020). https://doi.org/10.1038/s41586-020-3001-6. arXiv:2005.11495 [nucl-ex] [Erratum: Nature 590, E13 (2021)]

245. S. Acharya et al., Search for a common baryon source in high-multiplicity pp collisions at the LHC. Phys. Lett. B 811, 135849 (2020). https://doi.org/10.1016/j.physletb.2020.135849. arXiv:2004.08018 [nucl-ex]

246. S. Acharya et al., Investigation of the p-Σ0 interaction via femtoscopy in pp collisions. Phys. Lett. B 805, 135419 (2020). https://doi.org/10.1016/j.physletb.2020.135419. arXiv:1910.14407 [nucl-ex]

247. S. Acharya et al., Study of the Λ-Λ interaction with femtoscopy correlations in pp and p-Pb collisions at the LHC. Phys. Lett. B 797, 134822 (2019). https://doi.org/10.1016/j.physletb.2019.134822. arXiv:1905.07209 [nucl-ex]

248. S. Acharya et al., Exploring the NA-NΣ coupled system with high precision correlation techniques at the LHC. Phys. Lett. B 833, 137272 (2022). https://doi.org/10.1016/j.physletb.2022.137272. arXiv:2104.04427 [nucl-ex]

249. S. Acharya et al., First study of the two-body scattering involving charm hadrons. Phys. Rev. D 106(5), 052010 (2022). https://doi.org/10.1103/PhysRevD.106.052010. arXiv:2201.05352 [nucl-ex]

250. First measurement of the Λ-Ξ interaction in proton-proton collisions at the LHC (2022). https://doi.org/10.1016/j.physletb.2022.137223. arXiv:2204.10258 [nucl-ex]
251. Towards the understanding of the genuine three-body interaction for p–p–p and p–p–Λ (2022). arXiv:2206.03344 [nucl-ex]

252. S. Acharya et al., Scattering studies with low-energy kaon-proton femtoscopy in proton-proton collisions at the LHC. Phys. Rev. Lett. **124**(9), 092301 (2020). https://doi.org/10.1103/PhysRevLett.124.092301. arXiv:1905.13470 [nucl-ex]

253. S. Acharya et al., Kaon-proton strong interaction at low relative momentum via femtoscopy in Pb-Pb collisions at the LHC. Phys. Lett. B **822**, 136708 (2021). https://doi.org/10.1016/j.physletb.2021.136708. arXiv:2105.05683 [nucl-ex]

254. Constraining the $\overline{K}N$ coupled channel dynamics using femtoscopic correlations at the LHC (2022). arXiv:2205.15176 [nucl-ex]

255. K. Sasaki et al., $\Lambda\Lambda$ and $N\Xi$ interactions from lattice QCD near the physical point. Nucl. Phys. A **998**, 121737 (2020). https://doi.org/10.1016/j.nuclphysa.2020.121737. arXiv:1912.08630 [hep-lat]

256. D. Soldin, Update on the combined analysis of muon measurements from nine air shower experiments. PoS ICRC2021, 349 (2021). https://doi.org/10.22323/1.395.0349. arXiv:2108.08341 [astro-ph.HE]

257. H.P. Dembinski, R. Engel, A. Fedynitch, T. Gaisser, F. Riehn, T. Stanev, Data-driven model of the cosmic-ray flux and mass composition from 10 GeV to 10¹¹ GeV. PoS ICRC2017, 533 (2018). https://doi.org/10.22323/1.301.0533. arXiv:1711.11432 [astro-ph.HE]

258. T. Pierog, K. Werner, EPOS model and ultra high energy cosmic rays. Nucl. Phys. B Proc. Suppl. **196**, 102–105 (2009). https://doi.org/10.1016/j.nuclphysbps.2009.09.017. arXiv:0905.1198 [hep-ph]

259. T. Pierog, I. Karpenko, J.M. Katzy, E. Yatsenko, K. Werner, EPOS LHC: test of collective hadronization with data measured at the CERN Large Hadron Collider. Phys. Rev. C **92**(3), 034906 (2015). https://doi.org/10.1103/PhysRevC.92.034906. arXiv:1306.0121 [hep-ph]

260. C. Bierlich, J.R. Christiansen, Effects of color reconnection on hadron flavor observables. Phys. Rev. D **92**(9), 094010 (2015). https://doi.org/10.1103/PhysRevD.92.094010. arXiv:1507.02091 [hep-ph]

261. C. Bierlich, G. Gustafson, L. Lönnblad, Collectivity without plasma in hadronic collisions. Phys. Lett. B **779**, 58–63 (2018). https://doi.org/10.1016/j.physletb.2018.01.069. arXiv:1710.09725 [hep-ph]

262. J. Albrecht et al., The muon puzzle in cosmic-ray induced air showers and its connection to the large hadron collider. Astrophys. Space Sci. **367**(3), 27 (2022). https://doi.org/10.1007/s10509-022-04054-5. arXiv:2105.06148 [astro-ph.HE]

263. R. Ulrich, R. Engel, M. Unger, Hadronic multiparticle production at ultra-high energies and extensive air showers. Phys. Rev. D **83**, 054026 (2011). https://doi.org/10.1103/PhysRevD.83.054026. arXiv:1010.4310 [hep-ph]

264. S. Baur, H. Dembinski, M. Perlin, T. Pierog, R. Ulrich, K. Werner, Core-corona effect in hadron collisions and muon production in air showers (2019). arXiv:1902.09265 [hep-ph]

265. T. Pierog, S. Baur, H. Dembinski, M. Perlin, R. Ulrich, K. Werner, When heavy ions meet cosmic rays: potential impact of QGP formation on the muon puzzle. PoS ICRC2021, 469 (2021). https://doi.org/10.22323/1.395.0469

266. L.A. Anchordoqui, C. García Canal, S.J. Sciutto, J.F. Soriano, Through the looking-glass with ALICE into the quark-gluon plasma: a new test for hadronic interaction models used in air shower simulations. Phys. Lett. B **810**, 135837 (2020). https://doi.org/10.1016/j.physletb.2020.135837. arXiv:1907.09816 [hep-ph]

267. O. Adriani et al., Measurements of longitudinal and transverse momentum distributions for neutral pions in the forward-

rapidity region with the LHC detector. Phys. Rev. D **94**(3), 032007 (2016). https://doi.org/10.1103/PhysRevD.94.032007. arXiv:1507.08764 [hep-ex]

268. O. Adriani et al., Measurement of forward photon production cross-section in proton-proton collisions at $\sqrt{s}$ = 13 TeV with the LHC detector. Phys. Lett. B **780**, 233–239 (2018). https://doi.org/10.1016/j.physletb.2017.12.050. arXiv:1703.07678 [hep-ex]

269. O. Adriani et al., Measurement of energy flow, cross section and average inelasticity of forward neutrons produced in $\sqrt{s}$ = 13 TeV proton-proton collisions with the LHC Arm2 detector. JHEP **07**, 016 (2020). https://doi.org/10.1007/JHEP07(2020)016. arXiv:2003.02192 [hep-ex]

270. A.M. Sirunyan et al., Measurement of the average very forward energy as a function of the track multiplicity at central pseudorapidities in proton-proton collisions at $\sqrt{s}$ = 13 TeV. Eur. Phys. J. C **79**(11), 893 (2019). https://doi.org/10.1140/epjc/s10052-019-7402-3. arXiv:1908.01750 [hep-ex]

271. R. Aaij et al., Measurement of prompt charged-particle production in pp collisions at $\sqrt{s}$ = 13 TeV. JHEP **01**, 166 (2022). https://doi.org/10.1007/JHEP01(2022)166. arXiv:2107.10090 [hep-ex]

272. Nuclear modification factor of neutral pions in the forward and backward regions in $p\bar{p}$ collisions (2022). arXiv:2204.10608 [nucl-ex]

273. Evidence for modification of b quark hadronization in high-multiplicity $pp$ collisions at $\sqrt{s}$ = 13 TeV (2022). arXiv:2204.13042 [hep-ex]

274. O. Adriani, G.C. Barbarino, G.A. Bazilevskaya, R. Bellotti, M. Boezio, E.A. Bogomolov, L. Bonechi, M. Bongi, V. Bonvicini, S. Borisov, Pamela results on the cosmic-ray antiproton flux from 60 mev to 180 gev in kinetic energy. Phys. Rev. Lett. **105**(12) (2010). https://doi.org/10.1103/physrevlett.105.121101

275. M. Aguilar, L. Ali Cavasonza, G. Ambrosi et al., Towards understanding the origin of cosmic-ray positrons. Phys. Rev. Lett. **122**, 041102 (2019). https://doi.org/10.1103/PhysRevLett.122.041102

276. Direct detection of a break in the teraelectronvolt cosmic-ray spectrum of electrons and positrons. Nature **552**(7683), 63–66 (2017). https://doi.org/10.1038/nature24475

277. M. Aguilar et al., Antiproton flux, antiproton-to-proton flux ratio, and properties of elementary particle fluxes in primary cosmic rays measured with the alpha magnetic spectrometer on the international space station. Phys. Rev. Lett. **117**(9), 091103 (2016). https://doi.org/10.1103/PhysRevLett.117.091103

278. W.B. Atwood et al., The large area telescope on the fermi gamma-ray space telescope mission. Astrophys. J. **697**, 1071–1102 (2009). https://doi.org/10.1088/0004-637X/697/2/1071. arXiv:0902.1089 [astro-ph.IM]

279. T. Aramaki, C.J. Hailey, S.E. Boggs, P. Doetinchem, H. Fuke, S.I. Mognet, R.A. Ong, K. Perez, J. Zweerink, Antideuteron sensitivity for the GAPS experiment. Astropart. Phys. **74**, 6–13 (2016). https://doi.org/10.1016/j.astropartphys.2015.09.001

280. G. Giesen, M. Boudaud, Y. Génolini, V. Poulin, M. Cirelli, P. Salati, P.D. Serpico, Ams-02 antiprotons, at last! secondary astrophysical component and immediate implications for dark matter. J. Cosmol. Astropart. Phys. **2015**(09), 023 (2015). https://doi.org/10.1088/1475-7516/2015/09/023

281. A. Cuoco, J. Heisig, L. Klamt, M. Korsmeier, M. Krämer, Scrutinizing the evidence for dark matter in cosmic-ray antiprotons. Phys. Rev. D **99**(10) (2019). https://doi.org/10.1103/physrevd.99.103014

282. M. Di Mauro, M.W. Winkler, Multimessenger constraints on the dark matter interpretation of the fermi-lat galactic center excess. Phys. Rev. D **103**, 123005 (2021). https://doi.org/10.1103/PhysRevD.103.123005

283. L. Orusa, S. Manconi, F. Donato, M.D. Mauro, Constraining positron emission from pulsar populations with AMS-02 data. J.
Cosmol. Astropart. Phys. **2021**(12), 014 (2021). https://doi.org/10.1088/1475-7516/2021/12/014

284. M. Korsmeier, F. Donato, M.D. Mauro, Production cross sections of cosmic antiprotons in the light of new data from the NA61 and LHCb experiments. Phys. Rev. D **97**(10) (2018). https://doi.org/10.1103/physrevd.97.103019

285. L. Orusa, M.D. Mauro, F. Donato, M. Korsmeier, New determination of the production cross section for secondary positrons and electrons in the galaxy. Phys. Rev. D **105**(12) (2022). https://doi.org/10.1103/physrevd.105.123021

286. L. Orusa, M. Di Mauro, F. Donato, M. Korsmeier, New determination of the production cross section for γ rays in the galaxy. Phys. Rev. D **107**, 083031 (2023). https://doi.org/10.1103/PhysRevD.107.083031

287. M.W. Winkler, Cosmic ray antiprotons at high energies. J. Cosmol. Astropart. Phys. **2017**(02), 048 (2017). https://doi.org/10.1088/1475-7516/2017/02/048

288. F. Donato, N. Fornengo, P. Salati, Antideuterons as a signature of supersymmetric dark matter. Phys. Rev. D **62**(4) (2000). https://doi.org/10.1103/physrevd.62.043003

289. M. Cirelli, N. Fornengo, M. Taoso, A. Vittino, Anti-helium from dark matter annihilations. J. High Energy Phys. **2014**(8) (2014). https://doi.org/10.1007/jhep08(2014)009

290. M.W. Winkler, T. Linden, Dark matter annihilation can produce a detectable antihelium flux through $\Lambda_b$ decays. Phys. Rev. Lett. **126**, 101101 (2021). https://doi.org/10.1103/PhysRevLett.126.101101

291. M. Di Mauro, F. Donato, M. Korsmeier, S. Manconi, L. Orusa, Novel prediction for secondary positrons and electrons in the galaxy. Phys. Rev. D **108**, 063024 (2023). https://doi.org/10.1103/PhysRevD.108.063024

292. S. Abdollahi et al., Incremental fermi large area telescope fourth source catalog. Astrophys. J. Suppl. **260**(2), 53 (2022). https://doi.org/10.3847/1538-4365/ac6751. arXiv:2201.11184 [astro-ph.HE]

293. A. Aduszkiewicz et al., Measurements of $\pi^{\pm}$, $K^{\pm}$, $p$ and $\bar{p}$ spectra in proton-proton interactions at 20, 31, 40, 80 and 158 GeV/$c$ with the NA61/SHINE spectrometer at the CERN SPS. Eur. Phys. J. C **77**(10), 671 (2017). https://doi.org/10.1140/epjc/s10052-017-5260-4. arXiv:1705.02467 [nucl-ex]

294. B. Adams et al., Letter of intent: a new QCD facility at the M2 beam line of the CERN SPS (COMPASS++/AMBER) (2019)

295. R. Aaij et al., Measurement of antiproton production in p-He collisions at $\sqrt{s_{NN}}$ = 110GeV. Phys. Rev. Lett. **121**, 222001 (2018). https://doi.org/10.1103/PhysRevLett.121.222001

296. R. Aaij et al., Measurement of antiproton production from antihyperon decays in pHe collisions at $\sqrt{s_{NN}}$ = 110 GeV. Eur. Phys. J. C **83**, 543 (2023). https://doi.org/10.1140/epjc/s10052-023-11673-x. arXiv:2205.09009 [hep-ex]

297. L. Collaboration, LHCb SMOG Upgrade. Technical report, CERN, Geneva (2019). https://doi.org/10.17181/CERN.SAQC.EOWH. https://cds.cern.ch/record/2673690

298. S. Mariani, Fixed-target physics with the LHCb experiment at CERN Presented 08 Apr 2022. https://cds.cern.ch/record/2806641

299. A. Bursche, H.P. Dembinski, P. Di Nezza, M. Ferro-Luzzi, F. Fleuret, G. Graziani, G. Manca, E.A. Maurice, N. Neri, L.L. Papalardo, P. Robbe, M. Schmelling, M.A. Winn, V. Zhukov, Physics opportunities with the fixed-target program of the LHCb experiment using an unpolarized gas target. Technical report, CERN, Geneva (2018). https://cds.cern.ch/record/2649878