Interplay between core and corona from small to large systems

We present new results in $p$+$p$ and Pb+Pb collisions at the LHC energies from the updated dynamical core--corona initialization framework (DCCI2). The fractions of final hadron yields originating from equilibrated and non-equilibrated components are extracted as functions of multiplicity. We find that the contributions from non-equilibrated components are non-negligible even in Pb+Pb collisions and affect $p_T$-integrated multi-particle correlations. These suggest the importance of non-equilibrated components for the sophisticated extraction of properties of the quark gluon plasma from comparisons between dynamical frameworks and experimental data.


Introduction
The properties of the quark-gluon plasma (QGP) have been investigated through dynamical frameworks based on relativistic hydrodynamics.Recently, several attempts from dynamical frameworks are reported [1][2][3][4][5][6][7][8][9][10] to quantitatively extract the QGP transport coefficients from the data-to-model comparison.This indicates that sophisticated description with dynamical frameworks is needed for a further investigation of transport coefficients.On the other hand, the possibility of the QGP formation in small colliding systems has been discussed following the report of strangeness enhancement in p+p collisions by the ALICE Collaboration [11].Thus, the extension of the applicability of those dynamical frameworks to small colliding systems is indispensable.
Motivated by these backgrounds, we build the dynamical core-corona initialization as a state-of-the-art dynamical framework that is capable of describing from small to large colliding systems [12][13][14].Under the core-corona picture [15], the system generated in highenergy nuclear collisions is described with two components: equilibrated matter (core) and non-equilibrated matter (corona).In our framework, we realize the dynamical description of the separation of the system into the core and corona by incorporating the picture with the novel dynamical initialization approach [16].

Dynamical Core-Corona Initialization framework
Under the dynamical initialization approach, initial conditions of QGP fluids are generated dynamically via source terms of relativistic hydrodynamic equations.We assume that the sys-tem is described with the two components, QGP fluids as core and non-equilibrated partons as corona, respectively, and consider the energy and momentum conservation for the entire system generated in a collision.Thus, the depositions of energy and momentum from nonequilibrated partons become sources of QGP fluids, and which makes it possible to model the generation of initial conditions of QGP fluids with keeping the energy and momentum conservation of the entire system.The dynamical realization of the core-corona picture is attained by modeling the energy and momentum deposition rates of non-equilibrated partons.Based on the core-corona picture, non-equilibrated partons with low transverse momentum, p T , and in high-density regions are likely to deposit their energy and momentum while those with high p T and in low-density regions are not.
A brief summary of the latest version of our framework [14,17] is as follows: Initial conditions are phase-space distributions of partons after the final-state radiation obtained from Pythia8 [18] and Pythia8 Angantyr [19,20] for p+p and A+A collisions, respectively.The dynamical initialization is performed with those initially produced partons and the system is dynamically separated into the core and corona.The space-time evolution of the core, QGP fluids, is described with (3+1)-D ideal hydrodynamics [21] by incorporating the equation of state s95p-v.1.1 [22].Once temperatures of fluids become lower than a switching temperature, fluids are particlized with iS3D [23], a Monte-Carlo sampler of hadrons from hydrodynamic fields.The non-equilibrium partons are hadronized via the Lund string fragmentation with Pythia8.Finally, those directly emitted hadrons are handed to Jam [24], a hadronic transport model that performs hadronic rescatterings and resonance decays.

Results
We simulate p+p and Pb+Pb collisions at the LHC energies with our framework.Some major parameters control the fraction of generated core and corona components.We determine them so that the full simulation results can reasonably describe the multiplicity dependence of the yield ratios of omega baryons to charged pions in p+p and Pb+Pb collisions [14] 1 reported by the ALICE Collaboration [11].By fixing those parameters, we extract the fractions of hadron yields originating from core and corona as functions of multiplicity near mid-pseudo-rapidity (|η| < 0.5) in p+p collisions at √ s = 7 and 13 TeV and Pb+Pb collisions at TeV as shown in Fig. 1.One sees that the fractions of core and corona components show clear multiplicity scaling regardless of the system size or collision energies.Thus, within our framework, final hadron yields can be one of the reasonable estimators to evaluate the fraction of equilibrated and non-equilibrated components.We also find that the fraction of core components becomes dominant when multiplicity exceeds ⟨dN ch /dη⟩ ≈ 20.It should be also noted that the fraction of core (corona) components does not reach ≈ 100% (≈ 0%) even at the highest multiplicity bin in Pb+Pb collisions.For a further investigation, we show p T spectra of charged pions obtained in 20-40% centrality Pb+Pb collisions at √ s NN = 2.76 TeV in Fig. 2 (left).As an overall tendency, the core (corona) component dominates low (high) p T regions, which is exactly the consequence of implementing the core-corona picture.However, against the naive expectation of the above tendency, there is a slight enhancement of the corona components at the very low p T regions (p T < 1 GeV).This can be intuitively understood just as a result of the interplay between the two different shapes of p T spectra: exponential for the core components and power law for the corona components.The hadronic productions at such low p T regions in the p T spectra of corona components originate from the feed-down from string fragmentation of non-  equilibrated partons that lose their initial energy and momentum in dynamical initialization [14].
To see the influence of the existence of the corona components in the very low p T on observables, four-particle cumulant of charged particles, c 2 {4}, as a function of the number of charged particles, N ch , produced near midrapidity obtained from Pb+Pb collisions at √ s NN = 2.76 TeV is shown in Fig. 2 (right).For the kinematic ranges, c 2 {4} and N ch are obtained from charged particles with |η| < 0.8 and 0.2 < p T < 3.0 GeV which is the same ranges used in Ref. [25].The results from simulations switching off hadronic rescatterings and from the core components show the negative values of c 2 {4}, while those from the corona components are zero-consistent in the entire range of N ch .Notably, from the comparisons among those three results, one sees that the absolute values of c 2 {4} solely from the core components are diluted due to the existence of the corona contributions.

Figure 1 .
Figure 1.Fractions of hadron yields originating from the core and corona as functions of the number of charged particles produced at midrapidity |η| < 0.5.Results in p+p collisions at √ s = 7 and 13 TeV (open squares and triangles, respectively) and Pb+Pb collisions at √ s NN = 2.76 TeV (closed circles) are shown.

Figure 2 .
Figure 2. (Left) Charged pion p T spectra obtained from 20-40% centrality class of Pb+Pb collisions at √ s NN = 2.76 TeV are shown in the upper panel.Fractions of core and corona contributions in each p T bin are shown in the lower panel.(Right) Four-particle cumulant of charged particles as a function of the number of charged particles produced at midrapidity obtained in Pb+Pb collisions at √ s NN = 2.76 TeV.The solid lines are guides to the eyes.In both left and right figures, results with switching off hadronic rescatterings are shown in orange crosses and their breakdowns into the core and corona contributions are shown in red diagonal crosses and blue triangles, respectively.