Ultrarelativistic Collisions Research Articles

Thermodynamics of a rotating hadron resonance gas with van der Waals interaction

Studying the thermodynamics of the systems produced in ultra-relativistic heavy-ion collisions is crucial in understanding the QCD phase diagram. Recently, a new avenue has opened regarding the implications of large initial angular momentum and subsequent vorticity in the medium evolution in high-energy collisions. This adds a new type of chemical potential into the partonic and hadronic systems, called the rotational chemical potential. We study the thermodynamics of an interacting hadronic matter under rotation, formed in an ultra-relativistic collision. We introduce attractive and repulsive interactions through the van der Waals equation of state. Thermodynamic properties like the pressure (P), energy density (ε\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\varepsilon $$\\end{document}), entropy density (s), trace anomaly ((ε-3P)/T4\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$(\\varepsilon - 3P)/T^{4}$$\\end{document}), specific heat (cv\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$c_\ extrm{v}$$\\end{document}) and squared speed of sound (cs2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$c_\ extrm{s}^{2}$$\\end{document}) are studied as functions of temperature (T) for zero and finite rotation chemical potential. The conserved charge fluctuations, which can be quantified by their respective susceptibilities, are also studied. The rotational (spin) density corresponding to the rotational chemical potential is explored. In addition, we explore the possible liquid–gas phase transition in the hadron gas with van der Waals interaction in the T – ω\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\omega $$\\end{document} phase space.

Impact of medium anisotropy on quarkonium dissociation and regeneration

Quarkonium production in ultra-relativistic collisions plays a crucial role in probing the existence of hot QCD matter. This study explores quarkonia states dissociation and regeneration in the hot QCD medium while considering momentum anisotropy. The net quarkonia decay width (ΓD\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\Gamma _{D}$$\\end{document}) arises from two essential processes: collisional damping and gluonic dissociation. The quarkonia regeneration includes the transition from octet to singlet states within the anisotropic medium. Our study utilizes a medium-modified potential that incorporates anisotropy via particle distribution functions. This modified potential gives rise to collisional damping for quarkonia due to the surrounding medium, as well as the transition of quarkonia from singlet to octet states due to interactions with gluons. Furthermore, we employ the detailed balance approach to investigate the regeneration of quarkonia within this medium. Our comprehensive analysis spans various temperature settings, transverse momentum values, and anisotropic strengths. Notably, we find that, in addition to medium temperatures and heavy quark transverse momentum, anisotropy significantly influences the dissociation and regeneration of various quarkonia states.

Quasi particle model vs lattice QCD thermodynamics: extension to Nf=2+1+1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$N_f=2+1+1$$\\end{document} flavors and momentum dependent quark masses

In the last decade a quasi-particle model (QPM) has supplied the basis for the study of heavy quark (HQ) production in ultra-relativistic collisions, allowing for a phenomenological estimate of the HQ diffusion coefficient Ds(T)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$D_s(T)$$\\end{document}. Using the new lattice QCD results for the equation of state (EoS) with 2+1+1 dynamical flavors, we extend the QPM from Nf=2+1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$N_f=2+1$$\\end{document} to Nf=2+1+1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$N_f=2+1+1$$\\end{document}, where the charm quark is included. Fixing the coupling g(T) by a fit to the lQCD energy density ϵ(T)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\epsilon (T)$$\\end{document}, we evaluate the impact of different temperature parametrizations of charm quark mass on EoS and susceptibilities χq(T)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\chi _q(T)$$\\end{document} of light, χs(T)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\chi _s(T)$$\\end{document} of strange and χc(T)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\chi _c(T)$$\\end{document} of charm quarks, the last favouring a charm quark mass increasing toward Tc\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$T_c$$\\end{document}. We also explore the extension of the QPM to a more realistic approach called QPMp\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_p$$\\end{document}, where quark and gluon masses explicitly depend on their momentum converging to the current quark mass at high momenta, as expected from asymptotic free dynamics. The QPMp\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_p$$\\end{document} allows for a simultaneous quantitative description not only of the EoS but also of the quark susceptibilities (χq(T)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\chi _q(T)$$\\end{document}, χs(T)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\chi _s(T)$$\\end{document}), which instead are underestimated in the simple QPM. Furthermore, evaluating the spatial diffusion coefficient 2πTDs(T)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$2\\pi T D_s(T)$$\\end{document} in the QPMp\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_p$$\\end{document}, we find it is also closer than QPM to the recent lQCD data performed including dynamical fermions. Finally, in a 1+1D expanding system, we evaluate the RAA(pT)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$R_{AA}(p_T)$$\\end{document} in the QPM and QPMp\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_p$$\\end{document}, finding a significant reduction at low momenta for QPMp\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_p$$\\end{document} which could lead in a realistic scenario to a better agreement to experimental data.

Nonextensivity and temperature fluctuations of the Higgs boson production

We determine the temperature fluctuations associated with the Higgs boson pT spectrum through the derivation of the string tension distribution corresponding to the QCD-based Hagedorn function, frequently used to fit the transverse momentum distribution (TMD). The identified string tension fluctuations are heavy tailed, behaving similarly to the q-Gaussian distribution. After the convolution with the Schwinger mechanism, both approaches correctly describe the entire TMD. This approach leads to the nonthermal description of the particle production in ultrarelativistic pp collisions. By analyzing the data of pp collisions at s=13 TeV, we found that the average temperature associated with the Higgs boson differential cross section is around 85 times greater than the estimated value for the charged particle TMD. Our results show that the Higgs boson production exhibits the largest deviation from the thermal description. Published by the American Physical Society 2024

J/ψ and ψ(2S) polarization in proton-proton collisions at energies available at the CERN Large Hadron Collider using PYTHIA8

The production mechanisms of charmonium states in both hadronic and heavy-ion collisions hold great significance for investigating the hot and dense QCD matter. Studying charmonium polarization in ultrarelativistic collisions can also provide insights into the underlying production mechanisms. With this motivation, we explore the J/ψ and ψ(2S) polarization in proton+proton collisions at s=7, 8, and 13 TeV using a pQCD-inspired Monte Carlo event generator called 8. This work considers reconstructed quarkonia through their dimuons decay channel in the ALICE forward rapidity acceptance range of 2.5<yμμ<4. Further, we calculate the polarization parameters λθ, λϕ, λθϕ from the polar and azimuthal angular distributions of the dimuons in helicity and Collins-Soper frames. This study presents a comprehensive measurement of the polarization parameters as a function of transverse momentum, charged-particle multiplicity, and rapidity at the LHC energies. Our findings of charmonium polarization are in qualitative agreement with the corresponding experimental data. Published by the American Physical Society 2024

Femtoscopy with Lévy Sources from SPS through RHIC to LHC

Femtoscopy is a unique tool to investigate the space-time geometry of the matter created in ultra-relativistic collisions. If the probability density distribution of hadron emission is parametrized, then the dependence of its parameters on particle momentum, collision energy, and collision geometry can be given. In recent years, several measurements came to light that indicated the adequacy of assuming a Lévy-stable shape for the mentioned distribution. In parallel, several new phenomenological developments appeared, aiding the interpretation of the experimental results or providing tools for the measurements. In this paper, we discuss important aspects of femtoscopy with Lévy sources in light of some of these advances, including phenomenological and experimental ones.

Measurements of (anti)(hyper)nuclei with ALICE

The investigation of the production mechanism of light (anti)(hyper)nuclei in ultrarelativistic hadronic collisions is one of the main topics in modern nuclear physics. The ALICE Collaboration has significantly contributed to this specific field of research with systematic measurements of the production of (anti)(hyper)nuclei in different collision systems and centre-ofmass energies provided by the Large Hadron Collider. Measurements of the properties of hypernuclei, such as their lifetimes and binding energies, provide information on the hadronic interaction between hyperons and nucleons, which is complementary to that obtained from correlation measurements. In this contribution, recent results on the production of (anti)(hyper)nuclei measured with ALICE during the LHC Run 2 in different collision systems will be presented. These results will be discussed within the context of the statistical hadronization model and baryon coalescence. For the first time, the observations of the (anti)hyperhydrogen-4 and (anti)hyperhelium-4 in Pb–Pb collisions at 5.02 TeV will be shown. Moreover, new results on (anti)(hyper)nuclei measurements obtained using the LHC Run 3 data will be presented.

Impact of vorticity and viscosity on the hydrodynamic evolution of hot QCD medium

The strongly interacting transient state of quark-gluon plasma (QGP) medium created in ultra-relativistic collisions survives for a duration of a few fm/c. The spacetime evolution of QGP crucially depends on the equation of state (EoS), vorticity, viscosity, and external magnetic field. In the present study, we obtain the lifetime of a vortical QGP fluid within the ambit of relativistic second-order viscous hydrodynamics. We observe that the coupling of vorticity and viscosity significantly increases the lifetime of vortical QGP. The inclusion of a static magnetic field, vorticity, and viscosity makes the evolution slower. However, the static magnetic field slightly decreases the QGP lifetime by accelerating the evolution process for a non-rotating medium. We also report the rate of change of vorticity in the QGP, which will be helpful in studying the behavior of the medium in detail.

Geometric Poisson distribution of photons produced in the ultrarelativistic hadronic collisions

We show that the multiplicity distribution of photons produced with enhanced void probability in inelastic proton-proton collisions at sqrt{s} = 900 GeV, 2.76 TeV, and 7 TeV, measured at forward rapidities by the ALICE experiment at LHC, can be described by the geometric Poisson distribution. The traditionally used negative binomial distribution fails to reproduce the enhanced void probability and the shape of the modified combinants simultaneously. Our findings are relevant for the theoretical modeling of photon production processes in high-energy hadronic collisions.

Correlation of heavy-flavour hadron production and charged-particle multiplicity in pp collisions measured by ALICE

Heavy-flavour (HF) quarks, i.e. charm and beauty, are produced in the early stages of ultra-relativistic collisions via hard scattering processes. The measurement of heavy flavour as a function of charged-particle multiplicity not only provides information on how the production mechanisms are influenced by the event activity, but can also be considered as a tool to test the role of multiple parton interactions. Furthermore, charged-particle multiplicity studies are essential for reference measurements, to tune theoretical models and to provide information on the global characteristics of events. In ALICE, charged-particle multiplicity and heavy-flavour production in the hadronic and electronic decay channels are measured at central rapidity while HF-decay muons and J/ψ → μ + μ − are measured at forward rapidity. In this contribution, results will be presented on relative J/ψ, average D-mesons and leptons from HF hadron decay yields as a function of charged-particle multiplicity in proton-proton (pp) collisions at , 7, 8 and 13 TeV. These results will also be compared to theoretical model calculations.

Inferring nuclear structure from heavy isobar collisions using Trajectum

Nuclei with equal number of baryons but varying proton number (isobars) have many commonalities, but differ in both electric charge and nuclear structure. Relativistic collisions of such isobars provide unique opportunities to study the variation of the magnetic field, provided the nuclear structure is well understood. In this work we simulate collisions using several state-of-the-art parametrizations of the ^{96}_{40}4096Zr and ^{96}_{44}4496Ru isobars and show that a comparison with the exciting STAR measurement [arXiv:2109.00131] of ultrarelativistic collisions can uniquely identify the structure of both isobars. This not only provides an urgently needed understanding of the structure of the Zirconium and Ruthenium isobars, but also paves the way for more detailed studies of nuclear structure using relativistic heavy ion collisions.

Understanding the QCD medium by the diffusion of charm quarks using a color string percolation model

We study the drag and diffusion coefficients of the charm quark in the deconfined matter produced in the ultra-relativistic collisions by taking the Color String Percolation Model (CSPM) approach. CSPM, being a QCD-inspired model, can give us essential information about the hot and dense system produced in ultra-relativistic collisions. With the information on initial percolation temperature and percolation density, we estimate the relaxation time ($\tau_{c}$), drag coefficient ($\gamma$), transverse momentum diffusion coefficient ($B_{0}$), and spatial diffusion coefficient ($D_{s}$) of charm quark inside a deconfined medium. Finally, we compare the obtained results with lattice QCD and with various other theoretical models. A good agreement can be observed between the results obtained from CSPM and lattice QCD.

Nonequilibrium effects and transverse spherocity in ultrarelativistic proton-nucleus collisions

We investigate the effects of nonequilibrium dynamics in small colliding systems by comparing a nonequilibrium transport approach, Parton-Hadron-String Dynamics (PHSD), with a $(2+1)$-dimensional viscous hydrodynamic model, VISHNew. Focusing on $p+\mathrm{Pb}$ collisions at energies available at the Large Hadron Collider, we extract the initial conditions for the hydrodynamic model from PHSD, in order to reduce the impact of the early out-of-equilibrium dynamics and focus on the traces of nonequilibrium in the ensuing medium evolution. We find that in the transport approach quantities like energy density and bulk viscous pressure are highly inhomogeneous on the transverse plane during the whole evolution, whereas the hydrodynamic simulations dissolve more efficiently the initial spatial irregularities, still keeping a high degree of inhomogeneity due to the smaller size and lifetime of the medium produced in $p+\mathrm{Pb}$ collisions with respect to heavy-ion reactions. As a first step that will help to identify the impact of these nonequilibrium effects on final observables in proton-nucleus collisions, we perform an analysis of the transverse spherocity, an event-shape observable able to distinguish between jetty and isotropic configurations of transverse momenta. We found that the spherocity distribution in PHSD is slightly shifted towards the isotropic limit with respect to the hydrodynamic result. Even though this dissimilarity is partially due to the difference in the final charged particle production, it mainly comes from the different description within the two frameworks of the medium produced in small colliding systems. This finding supports the idea that multidifferential measurements, such as those based on event categorization according to multiplicity and spherocity, are useful to study final-state observables in ultrarelativistic proton-nucleus collisions.

Harmonic flow correlations in Au+Au reactions at 1.23 AGeV: a new testing ground for the equation-of-state and expansion geometry

Correlations between the harmonic flow coefficients v_1, v_2, v_3 and v_4 of nucleons in semi-peripheral Au+Au collisions at a beam energy of 1.23 AGeV are investigated within the hadronic transport approach ultra-relativistic quantum molecular dynamics (UrQMD). In contrast to ultra-relativistic collision energies (where the flow coefficients are evaluated with respect to the respective event plane), we predict strong correlations between the flow harmonics with respect to the reaction plane. Based on an event-by-event selection of the midrapidity final state elliptic flow of nucleons we show that as a function of rapidity, (I) the sign of the triangular flow changes, (II) that the shape of v_4 changes from convex to concave, and (III) that v_3propto v_1v_2 and v_4propto v_2^2 for all different event classes, indicating strong correlations between all investigated harmonic flow coefficients.

Jet Transport Coefficient at the Large Hadron Collider Energies in a Color String Percolation Approach

Within the color string percolation model (CSPM), jet transport coefficient, q^, is calculated for various multiplicity classes in proton-proton and centrality classes in nucleus-nucleus collisions at the Large Hadron Collider energies for a better understanding of the matter formed in ultra-relativistic collisions. q^ is studied as a function of final state charged particle multiplicity (pseudorapidity density at midrapidity), initial state percolation temperature and energy density. The CSPM results are then compared with different theoretical calculations from the JET Collaboration those incorporate particle energy loss in the medium.

New Capabilities of the FLUKA Multi-Purpose Code

FLUKA is a general purpose Monte Carlo code able to describe the transport and interaction of any particle and nucleus type in complex geometries over an energy range extending from thermal neutrons to ultrarelativistic hadron collisions. It has many different applications in accelerator design, detector studies, dosimetry, radiation protection, medical physics, and space research. In 2019, CERN and INFN, as FLUKA copyright holders, together decided to end their formal collaboration framework, allowing them henceforth to pursue different pathways aimed at meeting the evolving requirements of the FLUKA user community, and at ensuring the long term sustainability of the code. To this end, CERN set up the FLUKA.CERN Collaboration1. This paper illustrates the physics processes that have been newly released or are currently implemented in the code distributed by the FLUKA.CERN Collaboration2 under new licensing conditions that are meant to further facilitate access to the code, as well as intercomparisons. The description of coherent effects experienced by high energy hadron beams in crystal devices, relevant to promising beam manipulation techniques, and the charged particle tracking in vacuum regions subject to an electric field, overcoming a former lack, have already been made available to the users. Other features, namely the different kinds of low energy deuteron interactions as well as the synchrotron radiation emission in the course of charged particle transport in vacuum regions subject to magnetic fields, are currently undergoing systematic testing and benchmarking prior to release. FLUKA is widely used to evaluate radiobiological effects, with the powerful support of the Flair graphical interface, whose new generation (Available at http://flair.cern) offers now additional capabilities, e.g., advanced 3D visualization with photorealistic rendering and support for industry-standard volume visualization of medical phantoms. FLUKA has also been playing an extensive role in the characterization of radiation environments in which electronics operate. In parallel, it has been used to evaluate the response of electronics to a variety of conditions not included in radiation testing guidelines and standards for space and accelerators, and not accessible through conventional ground level testing. Instructive results have been obtained from Single Event Effects (SEE) simulations and benchmarks, when possible, for various radiation types and energies. The code has reached a high level of maturity, from which the FLUKA.CERN Collaboration is planning a substantial evolution of its present architecture. Moving towards a modern programming language allows to overcome fundamental constraints that limited development options. Our long term goal, in addition to improving and extending its physics performances with even more rigorous scientific oversight, is to modernize its structure to integrate independent contributions more easily and to formalize quality assurance through state-of-the-art software deployment techniques. This includes a continuous integration pipeline to automatically validate the codebase as well as automatic processing and analysis of a tailored physics-case test suite. With regard to the aforementioned objectives, several paths are currently envisaged, like finding synergies with Geant4, both at the core structure and interface level, this way offering the user the possibility to run with the same input different Monte Carlo codes and crosscheck the results.

Thermodynamic and transport properties of matter formed in pp, p-Pb, Xe–Xe and Pb–Pb collisions at the Large Hadron Collider using color string percolation model

To have a better understanding of the matter formed in ultra-relativistic collisions, using the color string percolation model, we have estimated initial energy density, mean free path, squared speed of sound, shear viscosity to entropy density ratio (η/s), bulk viscosity to entropy density ratio and bulk modulus for pp, p-Pb, Xe–Xe and Pb–Pb collisions at the Large Hadron Collider (LHC). These observables are studied as a function of final state charged particle multiplicity density normalized by nuclear overlap area and initial state percolation temperature. The results from this work are compared with the results from well known QCD models and with the transport properties of matter found in day-to-day life. These observables allow us to conclude an array of important findings about the matter formed in high-energy collisions at the LHC. For the high-multiplicity events, the estimated η/s are comparable with that is observed for heavy-ion collisions, which is an indication of a strongly coupled quark gluon plasma. This makes LHC high-multiplicity pp collision events as probable candidates for the search of quark-gluon plasma droplets.

Model calculations of φ-meson production in small collision systems

Ultrarelativistic ion collisions provide the unique possibility to study the quark-gluon plasma, a state of matter formed in the universe at the very first moments after the Big Bang. The minimal temperature and baryon density for the quark-gluon plasma formation requires scrutiny, since the signatures of the quark-gluon plasma formation are observed in large systems (such as Au+Au) at , whereas collective effects in p+p collisions are not revealed. The φ-meson production measurements are considered to be a convenient tool to investigate the collision dynamics, as it is sensitive to the quark-gluon plasma effects. To interpret the nuclear modification effects and to study the process of the possible QGP formation the comparison with different theoretical models predictions is needed. This paper presents the comparison of the obtained experimental results on φ-meson production in small collision systems (p+Al, p+Au) at to default and string melting versions of the AMPT model and PYTHIA model predictions. The results indicate that the minimal conditions (temperature and baryon density) for a QGP formation may lie in between in p+Al and p+Au collisions.

Deep Learning for the classification of quenched jets

An important aspect of the study of Quark-Gluon Plasma (QGP) in ultrarelativistic collisions of heavy ions is the ability to identify, in experimental data, a subset of the jets that were strongly modified by the interaction with the QGP. In this work, we propose studying Deep Learning techniques for this purpose. Samples of Z+jet events were simulated in vacuum (pp collisions) and medium (PbPb collisions) and used to train Deep Neural Networks with the objective of discriminating between medium- and vacuum-like jets within the medium (PbPb) sample. Dedicated Convolutional Neural Networks, Dense Neural Networks and Recurrent Neural Networks were developed and trained, and their performance was studied. Our results show the potential of these techniques for the identification of jet quenching effects induced by the presence of the QGP.

The signature of charge dependent directed flow observables by electromagnetic fields in heavy ion collisions

We discuss the generation of the directed flow v_1(p_T,y_z) induced by the electromagnetic field as a function of p_T and y_z. Despite the complex dynamics of charged particles due to strong interactions generating several anisotropies in the azimuthal angle, it is possible at p_T > m to directly correlate the splitting in v_1 of heavy quarks with different charges to some main features of the magnetic field, and in particular its values at formation and freeze-out time. We further found that the slope of the splitting dDelta v_1/dy_z|_{y_z=0} of positively and negatively charged particles at high p_T can be formulated as dDelta v_1/dy_z|_{y_z=0}=-alpha frac{partial ln f}{partial p_T}+frac{2alpha -beta }{p_T}, where f is the p_T spectra of the charged particles and the constants alpha and beta (order of MeV) are constrained by the y component of magnetic fields and the sign of alpha is simply determined by the difference Delta [tB_y(t)] in the center of colliding systems at the formation time of particles and at the time when particles leave the effective range of electromagnetic fields or freeze out. The formula is derived from general considerations and is confirmed by several related numerical simulations; it supplies a useful guide to quantify the effect of different magnetic field configurations and provides an evidence of why the measurement of Delta v_1 of charm, bottom and leptons from Z^0 decay and their correlations are a powerful probe of the initial e.m. fields in ultra-relativistic collisions.