Showing 1–50 of 84 results for author: Elfner, H

SMASH as an event generator for heavy-ion collisions

Authors: A. Sciarra, H. Elfner

Abstract: In this article we present an overview of the SMASH hadronic transport approach that is applied for non-equilibrium dynamics of hadrons in \heavyion\ collisions. We will give an overview about the ingredients of the approach and the applications for the dynamical description of \heavyion\ collisions and for calculations of fundamental properties of the hadron gas. The main emphasis of the article… ▽ More In this article we present an overview of the SMASH hadronic transport approach that is applied for non-equilibrium dynamics of hadrons in \heavyion\ collisions. We will give an overview about the ingredients of the approach and the applications for the dynamical description of \heavyion\ collisions and for calculations of fundamental properties of the hadron gas. The main emphasis of the article will be the infrastructure for sustainable software development that we have developed over the last 10 years including extensive unit tests and continuous integration. We will also provide one section about the performance of the code and how it can be analyzed and improved in the future. △ Less

Submitted 30 September, 2024; originally announced September 2024.

Comments: 9 pages, 4 figures, invited contribution to Frontiers on 'Promoting Green Computing in High-Energy Physics and Astrophysics'

Constraints on the Equation of State of Nuclear Matter from Systematically Comparing SMASH Calculations to HADES Data

Authors: Justin Mohs, Simon Spies, Hannah Elfner

Abstract: We aim to constrain the equation of state of nuclear matter by comparing calculations with the SMASH transport model to directed and elliptic flow measurements for protons and deuterons performed by the HADES collaboration in a systematic way. A momentum-dependent term is included in the potential for which we show that it is needed to describe flow data. We further incorporate a simple symmetry p… ▽ More We aim to constrain the equation of state of nuclear matter by comparing calculations with the SMASH transport model to directed and elliptic flow measurements for protons and deuterons performed by the HADES collaboration in a systematic way. A momentum-dependent term is included in the potential for which we show that it is needed to describe flow data. We further incorporate a simple symmetry potential in the transport model and present constraints on the stiffness of the equation of state of nuclear matter at saturation density and on the symmetry potential. The constraints are obtained by performing a Bayesian analysis such that we can also provide an uncertainty for the estimated parameters. The posterior distribution is obtained by Markov chain Monte Carlo sampling for which we emulate the transport model with a Gaussian process to lower computational costs. We find that a relatively stiff equation of state is favoured in our analysis with a small uncertainty whereas the constraints obtained for the symmetry potential are rather loose. △ Less

Submitted 25 September, 2024; originally announced September 2024.

Comments: 11 pages, 11 figures

Deuteron, triton, helium-3 and hypertriton production in relativistic heavy-ion collisions via stochastic multi-particle reactions

Authors: Martha Ege, Justin Mohs, Jan Staudenmaier, Hannah Elfner

Abstract: The production of light nuclei in heavy -ion collisions is an excellent probe for studying the phase diagram of quantum chromodynamics and for the search of a critical end point. In this work we apply a hybrid approach in which we study the light nuclei production in the afterburner stage of central Au+Au collisions at $\sqrt{s}_{NN}=7.7$, 14.5 and 19.6 GeV. In this stage, light nuclei are produce… ▽ More The production of light nuclei in heavy -ion collisions is an excellent probe for studying the phase diagram of quantum chromodynamics and for the search of a critical end point. In this work we apply a hybrid approach in which we study the light nuclei production in the afterburner stage of central Au+Au collisions at $\sqrt{s}_{NN}=7.7$, 14.5 and 19.6 GeV. In this stage, light nuclei are produced dynamically in $4\leftrightarrow 2$ catalysis reactions. A comparison of the dynamic production and a coalescence approach is presented for transverse momentum spectra of deuterons, tritons, $^3\rm He$ nuclei and hypertritons and ratios of light nuclei yields. A good agreement with the experimentally measured yield of nuclei is found and we proceed to further investigate the production mechanisms of light nuclei by calculating the rates of the important channels for the formation and disintegration. We find that the afterburner stage is essential for the description of light nuclei formation in heavy-ion collisions, as light nuclei undergo a large number of interactions. △ Less

Submitted 6 September, 2024; originally announced September 2024.

Comments: 15 pages, 25 figures

Distinguishing the sources of dielectron anisotropic flow at low beam energies

Authors: Renan Hirayama, Hannah Elfner

Abstract: We present calculations of dielectron anisotropic flow in heavy-ion collisions at HADES beam energies from a hadronic transport approach. The ongoing experimental analysis employs the traditional reaction plane method to evaluate the flow coefficients $v_n$ and claims to see isotropic radiation from the thermal quark-gluon plasma. We show in this work, that in the region above the pion mass, the d… ▽ More We present calculations of dielectron anisotropic flow in heavy-ion collisions at HADES beam energies from a hadronic transport approach. The ongoing experimental analysis employs the traditional reaction plane method to evaluate the flow coefficients $v_n$ and claims to see isotropic radiation from the thermal quark-gluon plasma. We show in this work, that in the region above the pion mass, the dilepton flow measurement might suffer from cancellation effects that mask the complicated underlying dynamics. Contributions from different baryonic and mesonic resonances show collective behaviour with different signs and lead to an overall vanishing elliptic flow. To differentiate the different contributions, we propose to employ the scalar product method, which exploits the previously measured hadronic flow to create different reference planes. As a proof of concept, we calculate the $v_2$ of dielectrons for Ag+Ag collisions at $\sqrt{s_{NN}}=2.55$ $\mathrm{GeV}$ with both methods and investigate the contribution of each source, concluding that the scalar product method provides the proton and pion tagged flow coefficients as two distinct measurements, disentangling the various dilepton sources. △ Less

Submitted 29 August, 2024; originally announced August 2024.

Comments: 10 pages, 7 figures

Bayesian Inference analysis of jet quenching using inclusive jet and hadron suppression measurements

Authors: R. Ehlers, Y. Chen, J. Mulligan, Y. Ji, A. Kumar, S. Mak, P. M. Jacobs, A. Majumder, A. Angerami, R. Arora, S. A. Bass, R. Datta, L. Du, H. Elfner, R. J. Fries, C. Gale, Y. He, B. V. Jacak, S. Jeon, F. Jonas, L. Kasper, M. Kordell II, R. Kunnawalkam-Elayavalli, J. Latessa, Y. -J. Lee , et al. (28 additional authors not shown)

Abstract: The JETSCAPE Collaboration reports a new determination of the jet transport parameter $\hat{q}$ in the Quark-Gluon Plasma (QGP) using Bayesian Inference, incorporating all available inclusive hadron and jet yield suppression data measured in heavy-ion collisions at RHIC and the LHC. This multi-observable analysis extends the previously published JETSCAPE Bayesian Inference determination of… ▽ More The JETSCAPE Collaboration reports a new determination of the jet transport parameter $\hat{q}$ in the Quark-Gluon Plasma (QGP) using Bayesian Inference, incorporating all available inclusive hadron and jet yield suppression data measured in heavy-ion collisions at RHIC and the LHC. This multi-observable analysis extends the previously published JETSCAPE Bayesian Inference determination of $\hat{q}$, which was based solely on a selection of inclusive hadron suppression data. JETSCAPE is a modular framework incorporating detailed dynamical models of QGP formation and evolution, and jet propagation and interaction in the QGP. Virtuality-dependent partonic energy loss in the QGP is modeled as a thermalized weakly-coupled plasma, with parameters determined from Bayesian calibration using soft-sector observables. This Bayesian calibration of $\hat{q}$ utilizes Active Learning, a machine--learning approach, for efficient exploitation of computing resources. The experimental data included in this analysis span a broad range in collision energy and centrality, and in transverse momentum. In order to explore the systematic dependence of the extracted parameter posterior distributions, several different calibrations are reported, based on combined jet and hadron data; on jet or hadron data separately; and on restricted kinematic or centrality ranges of the jet and hadron data. Tension is observed in comparison of these variations, providing new insights into the physics of jet transport in the QGP and its theoretical formulation. △ Less

Submitted 28 August, 2024; v1 submitted 15 August, 2024; originally announced August 2024.

Comments: 20 pages, 10 figures, 2 tables, submitted to PRC; updated acknowledgements

A soft-hard framework with exact four momentum conservation for small systems

Authors: I. Soudi, W. Zhao, A. Majumder, C. Shen, J. H. Putschke, B. Boudreaux, A. Angerami, R. Arora, S. A. Bass, Y. Chen, R. Datta, L. Du, R. Ehlers, H. Elfner, R. J. Fries, C. Gale, Y. He, B. V. Jacak, P. M. Jacobs, S. Jeon, Y. Ji, L. Kasper, M. Kelsey, M. Kordell II, A. Kumar , et al. (28 additional authors not shown)

Abstract: A new framework, called x-scape, for the combined study of both hard and soft transverse momentum sectors in high energy proton-proton ($p$-$p$) and proton-nucleus ($p$-$A$) collisions is set up. A dynamical initial state is set up using the 3d-Glauber model with transverse locations of hotspots within each incoming nucleon. A hard scattering that emanates from two colliding hotspots is carried ou… ▽ More A new framework, called x-scape, for the combined study of both hard and soft transverse momentum sectors in high energy proton-proton ($p$-$p$) and proton-nucleus ($p$-$A$) collisions is set up. A dynamical initial state is set up using the 3d-Glauber model with transverse locations of hotspots within each incoming nucleon. A hard scattering that emanates from two colliding hotspots is carried out using the Pythia generator. Initial state radiation from the incoming hard partons is carried out in a new module called I-matter, which includes the longitudinal location of initial splits. The energy-momentum of both the initial hard partons and their associated beam remnants is removed from the hot spots, depleting the energy-momentum available for the formation of the bulk medium. Outgoing showers are simulated using the matter generator, and results are presented for both cases, allowing for and not allowing for energy loss. First comparisons between this hard-soft model and single inclusive hadron and jet data from $p$-$p$ and minimum bias $p$-$Pb$ collisions are presented. Single hadron spectra in $p$-$p$ are used to carry out a limited (in number of parameters) Bayesian calibration of the model. Fair comparisons with data are indicative of the utility of this new framework. Theoretical studies of the correlation between jet $p_T$ and event activity at mid and forward rapidity are carried out. △ Less

Submitted 24 July, 2024; originally announced July 2024.

Comments: 18 pages, 15 figures

Flow and Equation of State of nuclear matter at $\mathbf{E_{\mathrm{kin}}}$/A=0.25-1.5 GeV with the SMASH transport approach

Authors: Lucia Anna Tarasovičová, Justin Mohs, Anton Andronic, Hannah Elfner, Karl-Heinz Kampert

Abstract: We present a comparison of directed and elliptic flow data by the FOPI collaboration in Au-Au, Xe-CsI, and Ni-Ni collisions at beam kinetic energies from 0.25 to 1.5 GeV per nucleon to simulations using the SMASH hadronic transport model. The Equation of State is parameterized as a function of nuclear density and momentum dependent potentials are newly introduced in SMASH. With a statistical analy… ▽ More We present a comparison of directed and elliptic flow data by the FOPI collaboration in Au-Au, Xe-CsI, and Ni-Ni collisions at beam kinetic energies from 0.25 to 1.5 GeV per nucleon to simulations using the SMASH hadronic transport model. The Equation of State is parameterized as a function of nuclear density and momentum dependent potentials are newly introduced in SMASH. With a statistical analysis, we show that the collective flow data at lower energies is in best agreement with a soft momentum dependent potential, while the elliptic flow at higher energies requires a harder momentum dependent Equation of State. △ Less

Submitted 16 May, 2024; originally announced May 2024.

BHAC-QGP: three-dimensional MHD simulations of relativistic heavy-ion collisions, II. Application to Au-Au collisions

Authors: Markus Mayer, Ashutosh Dash, Gabriele Inghirami, Hannah Elfner, Luciano Rezzolla, Dirk H. Rischke

Abstract: We present BHAC-QGP, a new numerical code to simulate the evolution of matter created in heavy-ion collisions. BHAC-QGP is based on the Black Hole Accretion Code (BHAC), which has been designed to model astrophysical processes through the solution of the equations of general-relativistic magnetohydrodynamics. Like the mother code, BHAC-QGP uses Adaptive Mesh Refinement (AMR), which allows for a dy… ▽ More We present BHAC-QGP, a new numerical code to simulate the evolution of matter created in heavy-ion collisions. BHAC-QGP is based on the Black Hole Accretion Code (BHAC), which has been designed to model astrophysical processes through the solution of the equations of general-relativistic magnetohydrodynamics. Like the mother code, BHAC-QGP uses Adaptive Mesh Refinement (AMR), which allows for a dynamic adjustment of the resolution in regions of the computational domain where a particularly high accuracy is needed. We here discuss a number of applications of BHAC-QGP to Au-Au collisions at Relativistic Heavy-Ion Collider (RHIC) energies and show that the code is able to reproduce results of other simulations of these scenarios, but with much higher accuracy. △ Less

Submitted 4 October, 2024; v1 submitted 13 March, 2024; originally announced March 2024.

Comments: 19 pages, 14 figures, new sections and figures added

BHAC-QGP: three-dimensional MHD simulations of relativistic heavy-ion collisions, I. Methods and tests

Authors: Markus Mayer, Ashutosh Dash, Gabriele Inghirami, Hannah Elfner, Luciano Rezzolla, Dirk H. Rischke

Abstract: We present BHAC-QGP, a new numerical code to simulate the evolution of matter created in heavy-ion collisions in the presence of electromagnetic fields. It is derived from the Black Hole Accretion Code (BHAC), which has been designed to model astrophysical processes in a general-relativistic magnetohydrodynamical description. As the original Black Hole Accretion Code, BHAC-QGP benefits from the us… ▽ More We present BHAC-QGP, a new numerical code to simulate the evolution of matter created in heavy-ion collisions in the presence of electromagnetic fields. It is derived from the Black Hole Accretion Code (BHAC), which has been designed to model astrophysical processes in a general-relativistic magnetohydrodynamical description. As the original Black Hole Accretion Code, BHAC-QGP benefits from the use of Adaptive Mesh Refinement (AMR), which allows us to dynamically adjust the resolution where necessary, and makes use of time-dependent Milne coordinates and the ultrarelativistic equation of state, $P = e/3$. We demonstrate that BHAC-QGP accurately passes a number of systematic and rigorous tests. △ Less

Submitted 4 October, 2024; v1 submitted 13 March, 2024; originally announced March 2024.

Comments: 20 pages, 9 figures, new sections and figures added

Photon-triggered jets as probes of multi-stage jet modification

Authors: C. Sirimanna, Y. Tachibana, A. Angerami, R. Arora, S. A. Bass, S. Cao, Y. Chen, L. Du, R. Ehlers, H. Elfner, W. Fan, R. J. Fries, C. Gale, Y. He, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, Y. Ji, L. Kasper, M. Kordell II, A. Kumar, R. Kunnawalkam-Elayavalli, J. Latessa, S. Lee , et al. (28 additional authors not shown)

Abstract: Prompt photons are created in the early stages of heavy ion collisions and traverse the QGP medium without any interaction. Therefore, photon-triggered jets can be used to study the jet quenching in the QGP medium. In this work, photon-triggered jets are studied through different jet and jet substructure observables for different collision systems and energies using the JETSCAPE framework. Since t… ▽ More Prompt photons are created in the early stages of heavy ion collisions and traverse the QGP medium without any interaction. Therefore, photon-triggered jets can be used to study the jet quenching in the QGP medium. In this work, photon-triggered jets are studied through different jet and jet substructure observables for different collision systems and energies using the JETSCAPE framework. Since the multistage evolution used in the JETSCAPE framework is adequate to describe a wide range of experimental observables simultaneously using the same parameter tune, we use the same parameters tuned for jet and leading hadron studies. The same isolation criteria used in the experimental analysis are used to identify prompt photons for better comparison. For the first time, high-accuracy JETSCAPE results are compared with multi-energy LHC and RHIC measurements to better understand the deviations observed in prior studies. This study highlights the importance of multistage evolution for the simultaneous description of experimental observables through different collision systems and energies using a single parameter tune. △ Less

Submitted 30 January, 2024; originally announced January 2024.

Comments: 4 pages, 5 figures. Proceedings of Quark Matter 2023 - XXXth International Conference on Ultra-relativistic Nucleus-Nucleus Collisions, Houston, TX, 3-9 September 2023

Measuring jet quenching with a Bayesian inference analysis of hadron and jet data by JETSCAPE

Authors: R. Ehlers, A. Angerami, R. Arora, S. A. Bass, S. Cao, Y. Chen, L. Du, H. Elfner, W. Fan, R. J. Fries, C. Gale, Y. He, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, Y. Ji, L. Kasper, M. Kordell II, A. Kumar, R. Kunnawalkam-Elayavalli, J. Latessa, S. Lee, Y. -J. Lee, D. Liyanage , et al. (28 additional authors not shown)

Abstract: The JETSCAPE Collaboration reports the first multi-messenger study of the QGP jet transport parameter $\hat{q}$ using Bayesian inference, incorporating all available hadron and jet inclusive yield and jet substructure data from RHIC and the LHC. The theoretical model utilizes virtuality-dependent in-medium partonic energy loss coupled to a detailed dynamical model of QGP evolution. Tension is obse… ▽ More The JETSCAPE Collaboration reports the first multi-messenger study of the QGP jet transport parameter $\hat{q}$ using Bayesian inference, incorporating all available hadron and jet inclusive yield and jet substructure data from RHIC and the LHC. The theoretical model utilizes virtuality-dependent in-medium partonic energy loss coupled to a detailed dynamical model of QGP evolution. Tension is observed when constraining $\hat{q}$ for different kinematic cuts of the inclusive hadron data. The addition of substructure data is shown to improve the constraint on $\hat{q}$, without inducing tension with the constraint due to inclusive observables. These studies provide new insight into the mechanisms of jet interactions in matter, and point to next steps in the field for comprehensive understanding of jet quenching as a probe of the QGP. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: 4 pages, 2 figures. Proceedings of Quark Matter 2023 - XXXth International Conference on Ultra-relativistic Nucleus-Nucleus Collisions, Houston, TX, 3-9 September 2023

3D Multi-system Bayesian Calibration with Energy Conservation to Study Rapidity-dependent Dynamics of Nuclear Collisions

Authors: Andi Mankolli, Aaron Angerami, Ritu Arora, Steffen Bass, Shanshan Cao, Yi Chen, Lipei Du, Raymond Ehlers, Hannah Elfner, Wenkai Fan, Rainer J. Fries, Charles Gale, Yayun He, Ulrich Heinz, Barbara Jacak, Peter Jacobs, Sangyong Jeon, Yi Ji, Lauren Kasper, Michael Kordell II, Amit Kumar, R. Kunnawalkam-Elayavalli, Joseph Latessa, Sook H. Lee, Yen-Jie Lee , et al. (26 additional authors not shown)

Abstract: Considerable information about the early-stage dynamics of heavy-ion collisions is encoded in the rapidity dependence of measurements. To leverage the large amount of experimental data, we perform a systematic analysis using three-dimensional hydrodynamic simulations of multiple collision systems -- large and small, symmetric and asymmetric. Specifically, we perform fully 3D multi-stage hydrodynam… ▽ More Considerable information about the early-stage dynamics of heavy-ion collisions is encoded in the rapidity dependence of measurements. To leverage the large amount of experimental data, we perform a systematic analysis using three-dimensional hydrodynamic simulations of multiple collision systems -- large and small, symmetric and asymmetric. Specifically, we perform fully 3D multi-stage hydrodynamic simulations initialized by a parameterized model for rapidity-dependent energy deposition, which we calibrate on the hadron multiplicity and anisotropic flow coefficients. We utilize Bayesian inference to constrain properties of the early- and late- time dynamics of the system, and highlight the impact of enforcing global energy conservation in our 3D model. △ Less

Submitted 31 December, 2023; originally announced January 2024.

Role of initial transverse momentum in a hybrid approach

Authors: Niklas Götz, Lucas Constantin, Hannah Elfner

Abstract: The purpose of this work is to study the effect of exchanging initial condition models in a modular hybrid approach. The focus lies on the event-by-event correlations of elliptic and triangular flow. This study is performed in the hybrid approach SMASH-vHLLE, composed of the hadronic transport approach SMASH and the (3+1)d viscous hydrodynamic code vHLLE. The initial condition models investigated… ▽ More The purpose of this work is to study the effect of exchanging initial condition models in a modular hybrid approach. The focus lies on the event-by-event correlations of elliptic and triangular flow. This study is performed in the hybrid approach SMASH-vHLLE, composed of the hadronic transport approach SMASH and the (3+1)d viscous hydrodynamic code vHLLE. The initial condition models investigated are SMASH IC, Trento and IP-Glasma. Correlations are calculated on an event-by-event basis between the eccentricities and momentum anisotropies of the initial state as well as the momentum anisotropies in the final state, both for ultra-central and off-central collisions for AuAu collisions at $\sqrt{s_{NN}} = 200$ GeV. This work demonstrates that, although averaged values for the eccentricities of these models are very similar, substantial differences exist both in the distributions of eccentricities, the correlations amongst the initial state properties as well as in the correlations between initial state and final state properties. Notably, whereas initial state momentum anisotropy is shown to not affect the final state flow, the presence of radial flow affects the emergence of final state momentum anisotropies. Inclusion of radial flow in the linear fit improves the prediction of final state flow from initial state properties. The presence of momentum in the initial state has an effect on the emergence of flow and is therefore a relevant part of initial state models, challenging the common understanding of final state momentum anisotropies being a linear response to initial state eccentricity only. △ Less

Submitted 8 December, 2023; originally announced December 2023.

Comments: 17 pages, 12 figures. Raw data of figures is provided in the supplemental material

Momentum anisotropy generation in a hybrid approach

Authors: Niklas Götz, Lucas Constantin, Hannah Elfner

Abstract: Anisotropic flow emerges in all three of hybrid approaches: initial conditions, viscous relativistic hydrodynamics as well as hadronic transport. Previous works focus mainly on a constant or temperature dependent shear viscosity $η/s$. Here instead, we study qualitatively the effect of a generalized $η/s(T,μ_B)$ in the hybrid approach SMASH-vHLLE-hybrid. The parameterization takes into account the… ▽ More Anisotropic flow emerges in all three of hybrid approaches: initial conditions, viscous relativistic hydrodynamics as well as hadronic transport. Previous works focus mainly on a constant or temperature dependent shear viscosity $η/s$. Here instead, we study qualitatively the effect of a generalized $η/s(T,μ_B)$ in the hybrid approach SMASH-vHLLE-hybrid. The parameterization takes into account the constraints of matching to the transport coefficients in the hadronic phase, as well as of recent Bayesian analysis results. We compare the effect of the different parameterizations in the intermediate energy region of $\sqrt{s_{NN}}$=7.7 - 39.0 GeV. We observe that using the energy density dependent parameterization decreases the effect of the point of particlization. In addition, we quantify the uncertainty due to different initial state profiles, including the SMASH initial conditions as well as TrENTo and IP-Glasma profiles. It can be shown that the initial state transverse momentum impacts final state momentum anisotropy. △ Less

Submitted 6 December, 2023; originally announced December 2023.

Comments: 4 pages, 3 figures. Presented at Hot Quarks 2022 Workshop for Young Scientists on the Physics of Ultrarelativistic Nucleus-Nucleus Collisions, Estes Park, USA

Charge and energy deposition in the McDIPPER framework

Authors: Oscar Garcia-Montero, Sören Schlichting, Hannah Elfner

Abstract: In this short note we present aspects of the energy and charge deposition within the McDIPPER, a novel 3D resolved model for the initial state of ultrarelativistic Heavy-Ion collisions based on the $k_\perp$-factorized Color Glass Condensate hybrid approach. This framework is a initial-state Monte Carlo event generator which deposits the relevant conserved charges (energy, charge and baryon densit… ▽ More In this short note we present aspects of the energy and charge deposition within the McDIPPER, a novel 3D resolved model for the initial state of ultrarelativistic Heavy-Ion collisions based on the $k_\perp$-factorized Color Glass Condensate hybrid approach. This framework is a initial-state Monte Carlo event generator which deposits the relevant conserved charges (energy, charge and baryon densities) both in the midrapidity and forward/backward regions of the collision. The event-by-event generator computes the gluon and (anti-) quark phase-space densities using the IP-Sat model, from where the conserved charges can be extracted directly. In this work we present the centrality and collision energy dependence for the deposited conserved quantities at midrapidity and the full event, the so-called $4π$ solid angle range. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Hybrid Hadronization of Jet Showers from $e^++e^-$ to $A+A$ with JETSCAPE

Authors: Cameron Parker, Aaron Angerami, Ritu Arora, Steffen Bass, Shanshan Cao, Yi Chen, Raymond Ehlers, Hannah Elfner, Wenkai Fan, Rainer J. Fries, Charles Gale, Yayun He, Ulrich Heinz, Barbara Jacak, Peter Jacobs, Sangyong Jeon, Yi Ji, Lauren Kasper, Michael Kordell II, Amit Kumar, Joseph Latessa, Yen-Jie Lee, Roy Lemmon, Dananjaya Liyanage, Arthur Lopez , et al. (26 additional authors not shown)

Abstract: In this talk we review jet production in a large variety of collision systems using the JETSCAPE event generator and Hybrid Hadronization. Hybrid Hadronization combines quark recombination, applicable when distances between partons in phase space are small, and string fragmentation appropriate for dilute parton systems. It can therefore smoothly describe the transition from very dilute parton syst… ▽ More In this talk we review jet production in a large variety of collision systems using the JETSCAPE event generator and Hybrid Hadronization. Hybrid Hadronization combines quark recombination, applicable when distances between partons in phase space are small, and string fragmentation appropriate for dilute parton systems. It can therefore smoothly describe the transition from very dilute parton systems like $e^++e^-$ to full $A+A$ collisions. We test this picture by using JETSCAPE to generate jets in various systems. Comparison to experimental data in $e^++e^-$ and $p+p$ collisions allows for a precise tuning of vacuum baseline parameters in JETSCAPE and Hybrid Hadronization. Proceeding to systems with jets embedded in a medium, we study in-medium hadronization for jet showers. We quantify the effects of an ambient medium, focusing in particular on the dependence on the collective flow and size of the medium. Our results clarify the effects we expect from in-medium hadronization of jets on observables like fragmentation functions, hadron chemistry and jet shape. △ Less

Submitted 7 November, 2023; v1 submitted 31 October, 2023; originally announced October 2023.

Comments: 6 pages, 4 figures, Hard Probes 2023 conference, accepted for publication in Proceedings of Science: version 2, references added, typos fixed

Photon momentum anisotropies from the late stages of relativistic heavy-ion collisions

Authors: Hannah Elfner, Niklas Götz, Oscar Garcia-Montero, Jean-Francois Paquet, Charles Gale

Abstract: The photon emission from the late stages of the dynamical evolution of heavy-ion reactions at the highest RHIC and LHC energies is investigated. A comparison between a calculation from hadronic rates from a fluid dynamic evolution down to temperatures of 120 MeV and a full non-equilibrium hadronic transport approach is performed. The photon yields are very similar in both cases while the elliptic… ▽ More The photon emission from the late stages of the dynamical evolution of heavy-ion reactions at the highest RHIC and LHC energies is investigated. A comparison between a calculation from hadronic rates from a fluid dynamic evolution down to temperatures of 120 MeV and a full non-equilibrium hadronic transport approach is performed. The photon yields are very similar in both cases while the elliptic flow is slightly smaller in the non-equilibrium scenario. This study is important, since it is crucial to apply the same dynamical evolution model for hadronic and electromagnetic observables. △ Less

Submitted 27 October, 2023; originally announced October 2023.

Comments: 6 pages, 4 figures, accepted in PoS for Hard Probes 2023

Fate of critical fluctuations in an interacting hadronic medium using maximum entropy distributions

Authors: Jan Hammelmann, Marcus Bluhm, Marlene Nahrang, Hannah Elfner

Abstract: We study the evolution of critical fluctuations in an expanding system within a hadronic transport approach. The initialization of the system with critical fluctuations is achieved by coupling the ideal hadron resonance gas cumulants to the ones from the 3d Ising model and generating the net and total particle number distribution from the principle of maximum entropy. These distributions are then… ▽ More We study the evolution of critical fluctuations in an expanding system within a hadronic transport approach. The initialization of the system with critical fluctuations is achieved by coupling the ideal hadron resonance gas cumulants to the ones from the 3d Ising model and generating the net and total particle number distribution from the principle of maximum entropy. These distributions are then evolved using realistic hadronic interactions. We systematically investigate the evolution of the critical fluctuations initialized at various temperatures and chemical potentials along a freeze-out line. We find that resonance regeneration and isospin randomization processes have the strongest influence on the evolution of the fluctuations. Additionally, the sets of particles coupled to the critical mode are modified to assess the strength of the propagation of correlations through interactions. We find that in the scaling region of the critical point correlations are propagated through the whole collisional history and are still present after the kinetic freeze-out of the matter if the coupling strength is large enough. △ Less

Submitted 10 October, 2023; originally announced October 2023.

Thermodynamics of an updated hadronic resonance list and influence on hadronic transport

Authors: Jordi Salinas San Martín, Renan Hirayama, Jan Hammelmann, Jamie M. Karthein, Paolo Parotto, Jacquelyn Noronha-Hostler, Claudia Ratti, Hannah Elfner

Abstract: Hadron lists based on experimental studies summarized by the Particle Data Group (PDG) are a crucial input for the equation of state and thermal models used in the study of strongly-interacting matter produced in heavy-ion collisions. Modeling of these strongly-interacting systems is carried out via hydrodynamical simulations, which are followed by hadronic transport codes that also require a hadr… ▽ More Hadron lists based on experimental studies summarized by the Particle Data Group (PDG) are a crucial input for the equation of state and thermal models used in the study of strongly-interacting matter produced in heavy-ion collisions. Modeling of these strongly-interacting systems is carried out via hydrodynamical simulations, which are followed by hadronic transport codes that also require a hadronic list as input. To remain consistent throughout the different stages of modeling of a heavy-ion collision, the same hadron list with its corresponding decays must be used at each step. It has been shown that even the most uncertain states listed in the PDG from 2016 are required to reproduce partial pressures and susceptibilities from Lattice Quantum Chromodynamics with the hadronic list known as the PDG2016+. Here, we update the hadronic list for use in heavy-ion collision modeling by including the latest experimental information for all states listed in the Particle Data Booklet in 2021. We then compare our new list, called PDG2021+, to Lattice Quantum Chromodynamics results and find that it achieves even better agreement with the first principles calculations than the PDG2016+ list. Furthermore, we develop a novel scheme based on intermediate decay channels that allows for only binary decays, such that PDG2021+ will be compatible with the hadronic transport framework SMASH. Finally, we use these results to make comparisons to experimental data and discuss the impact on particle yields and spectra. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Comments: 17 pages, 16 figures, 2 tables

The McDIPPER: A novel saturation-based 3+1D initial state model for Heavy Ion Collisions

Authors: Oscar Garcia-Montero, Hannah Elfner, Sören Schlichting

Abstract: We present a new 3D resolved model for the initial state of ultrarelativistic heavy-ion collisions, based on the $k_\perp$-factorized Color Glass Condensate hybrid approach. The McDIPPER framework responds to the need for a rapidity-resolved initial-state Monte Carlo event generator which can deposit the relevant conserved charges (energy, charge and baryon densities) both in the midrapidity and f… ▽ More We present a new 3D resolved model for the initial state of ultrarelativistic heavy-ion collisions, based on the $k_\perp$-factorized Color Glass Condensate hybrid approach. The McDIPPER framework responds to the need for a rapidity-resolved initial-state Monte Carlo event generator which can deposit the relevant conserved charges (energy, charge and baryon densities) both in the midrapidity and forward/backward regions of the collision. This event-by-event generator computes the gluon and (anti-) quark phase-space densities using the IP-Sat model, from where the relevant conserved charges can be computed directly. In the present work we have included the leading order contributions to the light flavor parton densities. As a feature, the model can be systematically improved in the future by adding next-to-leading order calculations (in the CGC hybrid framework), and extended to lower energies by including sub-eikonal corrections the channels included. We present relevant observables, such as the eccentricities and flow decorrelation, as tests of this new approach. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: 18 pages, 11 figures

A multistage framework for studying the evolution of jets and high-$p_T$ probes in small collision systems

Authors: Abhijit Majumder, Aaron Angerami, Ritu Arora, Steffen Bass, Shanshan Cao, Yi Chen, Raymond Ehlers, Hannah Elfner, Wenkai Fan, Rainer J. Fries, Charles Gale, Yayun He, Ulrich Heinz, Barbara Jacak, Peter Jacobs, Sangyong Jeon, Yi Ji, Lauren Kasper, Michael Kordell II, Amit Kumar, Joseph Latessa, Yen-Jie Lee, Roy Lemmon, Dananjaya Liyanage, Arthur Lopez , et al. (26 additional authors not shown)

Abstract: Understanding the modification of jets and high-$p_T$ probes in small systems requires the integration of soft and hard physics. We present recent developments in extending the JETSCAPE framework to build an event generator, which includes correlations between soft and hard partons, to study jet observables in small systems. The multi-scale physics of the collision is separated into different stag… ▽ More Understanding the modification of jets and high-$p_T$ probes in small systems requires the integration of soft and hard physics. We present recent developments in extending the JETSCAPE framework to build an event generator, which includes correlations between soft and hard partons, to study jet observables in small systems. The multi-scale physics of the collision is separated into different stages. Hard scatterings are first sampled at binary collision positions provided by the Glauber geometry. They are then propagated backward in space-time following an initial-state shower to obtain the initiating partons' energies and momenta before the collision. These energies and momenta are then subtracted from the incoming colliding nucleons for soft-particle production, modeled by the 3D-Glauber + hydrodynamics + hadronic transport framework. This new hybrid approach (X-SCAPE) includes non-trivial correlations between jet and soft particle productions in small systems. We calibrate this framework with the final state hadrons' $p_T$-spectra from low to high $p_T$ in $p$-$p$, and and then compare with the spectra in $p$-$Pb$ collisions from the LHC. We also present results for additional observables such as the distributions of event activity as a function of the hardest jet $p_T$ in forward and mid-rapidity for both $p$-$p$ and $p$-$Pb$ collisions. △ Less

Submitted 1 November, 2023; v1 submitted 4 August, 2023; originally announced August 2023.

Comments: 6 pages, 6 figures, Proceedings of Hard Probes 2023, 26-31 March 2023 Aschaffenburg, Germany

Journal ref: Proceedings of Science (Hard Probes 2023) 2023

Collision term dependence of the hadronic shear viscosity and diffusion coefficients

Authors: Jan Hammelmann, Jan Staudenmaier, Hannah Elfner

Abstract: In this work the shear viscosity $η$ and the diffusion coefficients of conserved charges $κ_{ij}$ with $i,j\in\{B,Q,S\}$ of hadronic matter are investigated within the hadronic transport approach SMASH. We systematically study the effect of multi-particle reactions, angular distributions and additional elastic cross sections via the additive quark model description, the numbers of degrees of freed… ▽ More In this work the shear viscosity $η$ and the diffusion coefficients of conserved charges $κ_{ij}$ with $i,j\in\{B,Q,S\}$ of hadronic matter are investigated within the hadronic transport approach SMASH. We systematically study the effect of multi-particle reactions, angular distributions and additional elastic cross sections via the additive quark model description, the numbers of degrees of freedom and the baryon density on the transport coefficients using the Green-Kubo formalism. We find that multi-particle reactions decrease the shear viscosity in a simplified hadron gas whereas the electric charge diffusion coefficient is not modified. Furthermore, additional elastic cross sections have a strong impact on both $η$ and $κ_{ij}$ whereas anisotropic scatterings enhance the shear viscosity in the full hadron gas. When increasing the number of degrees of freedom the shear viscosity is only slightly modified in comparison to the diffusion coefficients. Finally, the calculation within a finite baryon chemical potential reveals that the shear viscosity itself does not depend on $μ_B$ but on the ratio $η/ s$. The diffusion coefficients show a strong dependency which we compare to Chapman-Enskog calculations. △ Less

Submitted 28 July, 2023; originally announced July 2023.

A new metric improving Bayesian calibration of a multistage approach studying hadron and inclusive jet suppression

Authors: W. Fan, G. Vujanovic, S. A. Bass, A. Angerami, R. Arora, S. Cao, Y. Chen, T. Dai, L. Du, R. Ehlers, H. Elfner, R. J. Fries, C. Gale, Y. He, M. Heffernan, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, Y. Ji, L. Kasper, M. Kordell II, A. Kumar, J. Latessa, Y. -J. Lee , et al. (30 additional authors not shown)

Abstract: We study parton energy-momentum exchange with the quark gluon plasma (QGP) within a multistage approach composed of in-medium DGLAP evolution at high virtuality, and (linearized) Boltzmann Transport formalism at lower virtuality. This multistage simulation is then calibrated in comparison with high $p_T$ charged hadrons, D-mesons, and the inclusive jet nuclear modification factors, using Bayesian… ▽ More We study parton energy-momentum exchange with the quark gluon plasma (QGP) within a multistage approach composed of in-medium DGLAP evolution at high virtuality, and (linearized) Boltzmann Transport formalism at lower virtuality. This multistage simulation is then calibrated in comparison with high $p_T$ charged hadrons, D-mesons, and the inclusive jet nuclear modification factors, using Bayesian model-to-data comparison, to extract the virtuality-dependent transverse momentum broadening transport coefficient $\hat{q}$. To facilitate this undertaking, we develop a quantitative metric for validating the Bayesian workflow, which is used to analyze the sensitivity of various model parameters to individual observables. The usefulness of this new metric in improving Bayesian model emulation is shown to be highly beneficial for future such analyses. △ Less

Submitted 27 October, 2023; v1 submitted 18 July, 2023; originally announced July 2023.

Comments: 24 pages, 16 figures

Multiscale evolution of heavy flavor in the QGP

Authors: G. Vujanovic, A. Angerami, R. Arora, S. A. Bass, S. Cao, Y. Chen, T. Dai, L. Du, R. Ehlers, H. Elfner, W. Fan, R. J. Fries, C. Gale, Y. He, M. Heffernan, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, Y. Ji, L. Kasper, M. Kordell II, A. Kumar, J. Latessa, Y. -J. Lee , et al. (30 additional authors not shown)

Abstract: Shower development dynamics for a jet traveling through the quark-gluon plasma (QGP) is a multiscale process, where the heavy flavor mass is an important scale. During the high virtuality portion of the jet evolution in the QGP, emission of gluons from a heavy flavor is modified owing to heavy quark mass. Medium-induced radiation of heavy flavor is sensitive to microscopic processes (e.g. diffusio… ▽ More Shower development dynamics for a jet traveling through the quark-gluon plasma (QGP) is a multiscale process, where the heavy flavor mass is an important scale. During the high virtuality portion of the jet evolution in the QGP, emission of gluons from a heavy flavor is modified owing to heavy quark mass. Medium-induced radiation of heavy flavor is sensitive to microscopic processes (e.g. diffusion), whose virtuality dependence is phenomenologically explored in this study. In the lower virtuality part of shower evolution, i.e. when the mass is comparable to the virtuality of the parton, scattering and radiation processes of heavy quarks differ from light quarks. The effects of these mechanisms on shower development in heavy flavor tagged showers in the QGP is explored here. Furthermore, this multiscale study examines dynamical pair production of heavy flavor (via virtual gluon splittings) and their subsequent evolution in the QGP, which is not possible otherwise. A realistic event-by-event simulation is performed using the JETSCAPE framework. Energy-momentum exchange with the medium proceeds using a weak coupling recoil approach. Using leading hadron and open heavy flavor observables, differences in heavy versus light quark energy-loss mechanisms are explored, while the importance of heavy flavor pair production is highlighted along with future directions to study. △ Less

Submitted 27 October, 2023; v1 submitted 18 July, 2023; originally announced July 2023.

Comments: 6 pages, 3 figures, contribution to the Hard Probes 2023 proceedings

Effects of multi-scale jet-medium interactions on jet substructures

Authors: JETSCAPE Collaboration, Y. Tachibana, A. Angerami, R. Arora, S. A. Bass, S. Cao, Y. Chen, T. Dai, L. Du, R. Ehlers, H. Elfner, W. Fan, R. J. Fries, C. Gale, Y. He, M. Heffernan, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, Y. Ji, K. Kauder, L. Kasper, W. Ke, M. Kelsey , et al. (35 additional authors not shown)

Abstract: We utilize event-by-event Monte Carlo simulations within the JETSCAPE framework to examine scale-dependent jet-medium interactions in heavy-ion collisions. The reduction in jet-medium interaction during the early high-virtuality stage, where the medium is resolved at a short distance scale, is emphasized as a key element in explaining multiple jet observables, particularly substructures, simultane… ▽ More We utilize event-by-event Monte Carlo simulations within the JETSCAPE framework to examine scale-dependent jet-medium interactions in heavy-ion collisions. The reduction in jet-medium interaction during the early high-virtuality stage, where the medium is resolved at a short distance scale, is emphasized as a key element in explaining multiple jet observables, particularly substructures, simultaneously. By employing the MATTER+LBT setup, which incorporates this explicit reduction of medium effects at high virtuality, we investigate jet substructure observables, such as Soft Drop groomed observables. When contrasted with existing data, our findings spotlight the significant influence of the reduction at the early high-virtuality stages. Furthermore, we study the substructure of gamma-tagged jets, providing predictive insights for future experimental analyses. This broadens our understanding of the various contributing factors involved in modifying jet substructures. △ Less

Submitted 16 July, 2023; originally announced July 2023.

Comments: 6 pages, 4 figures, contribution to the Hard Probes 2023 proceedings

Hot QCD White Paper

Authors: M. Arslandok, S. A. Bass, A. A. Baty, I. Bautista, C. Beattie, F. Becattini, R. Bellwied, Y. Berdnikov, A. Berdnikov, J. Bielcik, J. T. Blair, F. Bock, B. Boimska, H. Bossi, H. Caines, Y. Chen, Y. -T. Chien, M. Chiu, M. E. Connors, M. Csanád, C. L. da Silva, A. P. Dash, G. David, K. Dehmelt, V. Dexheimer , et al. (149 additional authors not shown)

Abstract: Hot QCD physics studies the nuclear strong force under extreme temperature and densities. Experimentally these conditions are achieved via high-energy collisions of heavy ions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC). In the past decade, a unique and substantial suite of data was collected at RHIC and the LHC, probing hydrodynamics at the nucleon scale, the… ▽ More Hot QCD physics studies the nuclear strong force under extreme temperature and densities. Experimentally these conditions are achieved via high-energy collisions of heavy ions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC). In the past decade, a unique and substantial suite of data was collected at RHIC and the LHC, probing hydrodynamics at the nucleon scale, the temperature dependence of the transport properties of quark-gluon plasma, the phase diagram of nuclear matter, the interaction of quarks and gluons at different scales and much more. This document, as part of the 2023 nuclear science long range planning process, was written to review the progress in hot QCD since the 2015 Long Range Plan for Nuclear Science, as well as highlight the realization of previous recommendations, and present opportunities for the next decade, building on the accomplishments and investments made in theoretical developments and the construction of new detectors. Furthermore, this document provides additional context to support the recommendations voted on at the Joint Hot and Cold QCD Town Hall Meeting, which are reported in a separate document. △ Less

Submitted 30 March, 2023; originally announced March 2023.

Comments: 190 pages, 69 figures

Theoretical and Experimental Constraints for the Equation of State of Dense and Hot Matter

Authors: Rajesh Kumar, Veronica Dexheimer, Johannes Jahan, Jorge Noronha, Jacquelyn Noronha-Hostler, Claudia Ratti, Nico Yunes, Angel Rodrigo Nava Acuna, Mark Alford, Mahmudul Hasan Anik, Debarati Chatterjee, Katerina Chatziioannou, Hsin-Yu Chen, Alexander Clevinger, Carlos Conde, Nikolas Cruz-Camacho, Travis Dore, Christian Drischler, Hannah Elfner, Reed Essick, David Friedenberg, Suprovo Ghosh, Joaquin Grefa, Roland Haas, Alexander Haber , et al. (35 additional authors not shown)

Abstract: This review aims at providing an extensive discussion of modern constraints relevant for dense and hot strongly interacting matter. It includes theoretical first-principle results from lattice and perturbative QCD, as well as chiral effective field theory results. From the experimental side, it includes heavy-ion collision and low-energy nuclear physics results, as well as observations from neutro… ▽ More This review aims at providing an extensive discussion of modern constraints relevant for dense and hot strongly interacting matter. It includes theoretical first-principle results from lattice and perturbative QCD, as well as chiral effective field theory results. From the experimental side, it includes heavy-ion collision and low-energy nuclear physics results, as well as observations from neutron stars and their mergers. The validity of different constraints, concerning specific conditions and ranges of applicability, is also provided. △ Less

Submitted 12 June, 2024; v1 submitted 29 March, 2023; originally announced March 2023.

Journal ref: Living Reviews in Relativity (2024) 27:3

Dense Nuclear Matter Equation of State from Heavy-Ion Collisions

Authors: Agnieszka Sorensen, Kshitij Agarwal, Kyle W. Brown, Zbigniew Chajęcki, Paweł Danielewicz, Christian Drischler, Stefano Gandolfi, Jeremy W. Holt, Matthias Kaminski, Che-Ming Ko, Rohit Kumar, Bao-An Li, William G. Lynch, Alan B. McIntosh, William G. Newton, Scott Pratt, Oleh Savchuk, Maria Stefaniak, Ingo Tews, ManYee Betty Tsang, Ramona Vogt, Hermann Wolter, Hanna Zbroszczyk, Navid Abbasi, Jörg Aichelin , et al. (111 additional authors not shown)

Abstract: The nuclear equation of state (EOS) is at the center of numerous theoretical and experimental efforts in nuclear physics. With advances in microscopic theories for nuclear interactions, the availability of experiments probing nuclear matter under conditions not reached before, endeavors to develop sophisticated and reliable transport simulations to interpret these experiments, and the advent of mu… ▽ More The nuclear equation of state (EOS) is at the center of numerous theoretical and experimental efforts in nuclear physics. With advances in microscopic theories for nuclear interactions, the availability of experiments probing nuclear matter under conditions not reached before, endeavors to develop sophisticated and reliable transport simulations to interpret these experiments, and the advent of multi-messenger astronomy, the next decade will bring new opportunities for determining the nuclear matter EOS, elucidating its dependence on density, temperature, and isospin asymmetry. Among controlled terrestrial experiments, collisions of heavy nuclei at intermediate beam energies (from a few tens of MeV/nucleon to about 25 GeV/nucleon in the fixed-target frame) probe the widest ranges of baryon density and temperature, enabling studies of nuclear matter from a few tenths to about 5 times the nuclear saturation density and for temperatures from a few to well above a hundred MeV, respectively. Collisions of neutron-rich isotopes further bring the opportunity to probe effects due to the isospin asymmetry. However, capitalizing on the enormous scientific effort aimed at uncovering the dense nuclear matter EOS, both at RHIC and at FRIB as well as at other international facilities, depends on the continued development of state-of-the-art hadronic transport simulations. This white paper highlights the essential role that heavy-ion collision experiments and hadronic transport simulations play in understanding strong interactions in dense nuclear matter, with an emphasis on how these efforts can be used together with microscopic approaches and neutron star studies to uncover the nuclear EOS. △ Less

Submitted 25 January, 2024; v1 submitted 30 January, 2023; originally announced January 2023.

Comments: White paper prepared for the 2023 Long Range Plan. v3: Updated version as published in Progress in Particle and Nuclear Physics. Note: the published version does not include the executive summary; in the updated arXiv version, the executive summary is included as an appendix. v4: Corrected list of authors

Report number: INT-PUB-23-001, LA-UR-23-20514, LLNL-TR-844629

Journal ref: Prog.Part.Nucl.Phys. 134 (2024) 104080

Hard jet substructure in a multistage approach

Authors: Y. Tachibana, A. Kumar, A. Majumder, A. Angerami, R. Arora, S. A. Bass, S. Cao, Y. Chen, T. Dai, L. Du, R. Ehlers, H. Elfner, W. Fan, R. J. Fries, C. Gale, Y. He, M. Heffernan, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, Y. Ji, K. Kauder, L. Kasper, W. Ke , et al. (34 additional authors not shown)

Abstract: We present predictions and postdictions for a wide variety of hard jet-substructure observables using a multistage model within the JETSCAPE framework. The details of the multistage model and the various parameter choices are described in [A. Kumar et al., arXiv:2204.01163]. A novel feature of this model is the presence of two stages of jet modification: a high virtuality phase [modeled using the… ▽ More We present predictions and postdictions for a wide variety of hard jet-substructure observables using a multistage model within the JETSCAPE framework. The details of the multistage model and the various parameter choices are described in [A. Kumar et al., arXiv:2204.01163]. A novel feature of this model is the presence of two stages of jet modification: a high virtuality phase [modeled using the modular all twist transverse-scattering elastic-drag and radiation model (MATTER)], where modified coherence effects diminish medium-induced radiation, and a lower virtuality phase [modeled using the linear Boltzmann transport model (LBT)], where parton splits are fully resolved by the medium as they endure multiple scattering induced energy loss. Energy-loss calculations are carried out on event-by-event viscous fluid dynamic backgrounds constrained by experimental data. The uniform and consistent descriptions of multiple experimental observables demonstrate the essential role of modified coherence effects and the multistage modeling of jet evolution. Using the best choice of parameters from [A. Kumar et al., arXiv:2204.01163], and with no further tuning, we present calculations for the medium modified jet fragmentation function, the groomed jet momentum fraction $z_g$ and angular separation $r_g$ distributions, as well as the nuclear modification factor of groomed jets. These calculations provide accurate descriptions of published data from experiments at the Large Hadron Collider. Furthermore, we provide predictions from the multistage model for future measurements at the BNL Relativistic Heavy Ion Collider. △ Less

Submitted 16 October, 2024; v1 submitted 6 January, 2023; originally announced January 2023.

Comments: 20 pages, 12 figures

Journal ref: Phys. Rev. C 110, 044907 (2024)

Global angular momentum generation in heavy-ion reactions within a hadronic transport approach

Authors: Nils Sass, Marco Müller, Oscar Garcia-Montero, Hannah Elfner

Abstract: In 2017, the STAR collaboration at the Relativistic Heavy Ion Collider (RHIC) has measured finite global angular momentum in heavy-ion collisions through a spin polarization measurement of $Λ$ hyperons. This measurement revealed a high angular momentum of the heavy ions and provided experimental evidence for vorticity in the quark-gluon plasma (QGP) for the first time. In order to investigate the… ▽ More In 2017, the STAR collaboration at the Relativistic Heavy Ion Collider (RHIC) has measured finite global angular momentum in heavy-ion collisions through a spin polarization measurement of $Λ$ hyperons. This measurement revealed a high angular momentum of the heavy ions and provided experimental evidence for vorticity in the quark-gluon plasma (QGP) for the first time. In order to investigate the underlying mechanisms, a dynamic description of the transfer of angular momentum is required. In this work, the microscopic non-equilibrium transport approach SMASH (Simulating Many Accelerated Strongly-interacting Hadrons) is applied to study the generation of global angular momentum by the interaction of two nuclei. As SMASH provides access to the whole phase-space evolution of every particle at any given time, it allows to assess the fraction of angular momentum generated in the fireball by all participants. We confirm the previous modeling by Becattini \textit{et al} within a geometric Glauber model approach, which found that the angular momentum transfer reaches a unique maximum in mid-central collisions during time evolution. The corresponding impact parameter is around $b=4-6$ fm for all beam energies from $\sqrt{s_{\rm NN}}=2.41-200$ GeV. Even though angular momentum is not conserved locally in the transport approach a priori, we identify the contributions to the conservation violation and propose optimal setups for different energy regimes that recover conservation, based upon the test particle method and the treatment of Fermi motion. Furthermore, the system size and centrality dependence are investigated. △ Less

Submitted 12 January, 2024; v1 submitted 29 December, 2022; originally announced December 2022.

Comments: 8 pages, 6 figures

Journal ref: Phys.Rev.C 108 (2023) 4, 044903

Comprehensive Study of Multi-scale Jet-medium Interaction

Authors: Y. Tachibana, A. Angerami, R. Arora, S. A. Bass, S. Cao, Y. Chen, T. Dai, L. Du, R. Ehlers, H. Elfner, W. Fan, R. J. Fries, C. Gale, Y. He, M. Heffernan, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, Y. Ji, L. Kasper, W. Ke, M. Kelsey, M. Kordell II, A. Kumar , et al. (33 additional authors not shown)

Abstract: We explore jet-medium interactions at various scales in high-energy heavy-ion collisions using the JETSCAPE framework. The physics of the multi-stage modeling and the coherence effect at high virtuality is discussed through the results of multiple jet and high-$p_{\mathrm{T}}$ particle observables, compared with experimental data. Furthermore, we investigate the jet-medium interaction involved in… ▽ More We explore jet-medium interactions at various scales in high-energy heavy-ion collisions using the JETSCAPE framework. The physics of the multi-stage modeling and the coherence effect at high virtuality is discussed through the results of multiple jet and high-$p_{\mathrm{T}}$ particle observables, compared with experimental data. Furthermore, we investigate the jet-medium interaction involved in the hadronization process. △ Less

Submitted 23 December, 2022; originally announced December 2022.

Comments: 6 pages, 5 figures, contribution to the Quark Matter 2022 proceedings

Journal ref: Acta Phys. Pol. B Proc. Suppl. 16, 1-A50 (2023)

Long Range Plan: Dense matter theory for heavy-ion collisions and neutron stars

Authors: Alessandro Lovato, Travis Dore, Robert D. Pisarski, Bjoern Schenke, Katerina Chatziioannou, Jocelyn S. Read, Philippe Landry, Pawel Danielewicz, Dean Lee, Scott Pratt, Fabian Rennecke, Hannah Elfner, Veronica Dexheimer, Rajesh Kumar, Michael Strickland, Johannes Jahan, Claudia Ratti, Volodymyr Vovchenko, Mikhail Stephanov, Dekrayat Almaalol, Gordon Baym, Mauricio Hippert, Jacquelyn Noronha-Hostler, Jorge Noronha, Enrico Speranza , et al. (39 additional authors not shown)

Abstract: Since the release of the 2015 Long Range Plan in Nuclear Physics, major events have occurred that reshaped our understanding of quantum chromodynamics (QCD) and nuclear matter at large densities, in and out of equilibrium. The US nuclear community has an opportunity to capitalize on advances in astrophysical observations and nuclear experiments and engage in an interdisciplinary effort in the theo… ▽ More Since the release of the 2015 Long Range Plan in Nuclear Physics, major events have occurred that reshaped our understanding of quantum chromodynamics (QCD) and nuclear matter at large densities, in and out of equilibrium. The US nuclear community has an opportunity to capitalize on advances in astrophysical observations and nuclear experiments and engage in an interdisciplinary effort in the theory of dense baryonic matter that connects low- and high-energy nuclear physics, astrophysics, gravitational waves physics, and data science △ Less

Submitted 7 November, 2022; v1 submitted 3 November, 2022; originally announced November 2022.

Comments: 70 pages, 3 figures, White Paper for the Long Range Plan for Nuclear Science

Report number: LA-UR-22-31648

The exploration of hot and dense nuclear matter: Introduction to relativistic heavy-ion physics

Authors: Hannah Elfner, Berndt Müller

Abstract: This article summarizes our present knowledge about nuclear matter at the highest energy densities and its formation in relativistic heavy ion collisions. We review what is known about the structure and properties of the quark-gluon plasma and survey the observables that are used to glean information about it from experimental data. This article summarizes our present knowledge about nuclear matter at the highest energy densities and its formation in relativistic heavy ion collisions. We review what is known about the structure and properties of the quark-gluon plasma and survey the observables that are used to glean information about it from experimental data. △ Less

Submitted 21 October, 2022; originally announced October 2022.

Comments: 110 pages

Bayesian analysis of QGP jet transport using multi-scale modeling applied to inclusive hadron and reconstructed jet data

Authors: R. Ehlers, A. Angerami, R. Arora, S. A. Bass, S. Cao, Y. Chen, L. Du, T. Dai, H. Elfner, W. Fan, R. J. Fries, C. Gale, Y. He, M. Heffernan, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, Y. Ji, L. Kasper, W. Ke, M. Kelsey, M. Kordell II, A. Kumar, J. Latessa , et al. (33 additional authors not shown)

Abstract: The JETSCAPE Collaboration reports a new determination of jet transport coefficients in the Quark-Gluon Plasma, using both reconstructed jet and hadron data measured at RHIC and the LHC. The JETSCAPE framework incorporates detailed modeling of the dynamical evolution of the QGP; a multi-stage theoretical approach to in-medium jet evolution and medium response; and Bayesian inference for quantitati… ▽ More The JETSCAPE Collaboration reports a new determination of jet transport coefficients in the Quark-Gluon Plasma, using both reconstructed jet and hadron data measured at RHIC and the LHC. The JETSCAPE framework incorporates detailed modeling of the dynamical evolution of the QGP; a multi-stage theoretical approach to in-medium jet evolution and medium response; and Bayesian inference for quantitative comparison of model calculations and data. The multi-stage framework incorporates multiple models to cover a broad range in scale of the in-medium parton shower evolution, with dynamical choice of model that depends on the current virtuality or energy of the parton. We will discuss the physics of the multi-stage modeling, and then present a new Bayesian analysis incorporating it. This analysis extends the recently published JETSCAPE determination of the jet transport parameter $\hat{q}$ that was based solely on inclusive hadron suppression data, by incorporating reconstructed jet measurements of quenching. We explore the functional dependence of jet transport coefficients on QGP temperature and jet energy and virtuality, and report the consistency and tensions found for current jet quenching modeling with hadron and reconstructed jet data over a wide range in kinematics and $\sqrt{s_{\text{NN}}}$. This analysis represents the next step in the program of comprehensive analysis of jet quenching phenomenology and its constraint of properties of the QGP. △ Less

Submitted 16 August, 2022; originally announced August 2022.

Comments: 6 pages, 2 figures, contribution to the Quark Matter 2022 proceedings

Multi-scale evolution of charmed particles in a nuclear medium

Authors: JETSCAPE collaboration, W. Fan, G. Vujanovic, S. A. Bass, A. Majumder, A. Angerami, R. Arora, S. Cao, Y. Chen, T. Dai, L. Du, R. Ehlers, H. Elfner, R. J. Fries, C. Gale, Y. He, M. Heffernan, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, Y. Ji, K. Kauder, L. Kasper, W. Ke , et al. (35 additional authors not shown)

Abstract: Parton energy-momentum exchange with the quark gluon plasma (QGP) is a multi-scale problem. In this work, we calculate the interaction of charm quarks with the QGP within the higher twist formalism at high virtuality and high energy using the MATTER model, while the low virtuality and high energy portion is treated via a (linearized) Boltzmann Transport (LBT) formalism. Coherence effect that reduc… ▽ More Parton energy-momentum exchange with the quark gluon plasma (QGP) is a multi-scale problem. In this work, we calculate the interaction of charm quarks with the QGP within the higher twist formalism at high virtuality and high energy using the MATTER model, while the low virtuality and high energy portion is treated via a (linearized) Boltzmann Transport (LBT) formalism. Coherence effect that reduces the medium-induced emission rate in the MATTER model is also taken into account. The interplay between these two formalisms is studied in detail and used to produce a good description of the D-meson and charged hadron nuclear modification factor RAA across multiple centralities. All calculations were carried out utilizing the JETSCAPE framework. △ Less

Submitted 13 May, 2023; v1 submitted 1 August, 2022; originally announced August 2022.

Journal ref: 107(5), 2023, 054901

Influence of heavy resonances in SMASH

Authors: J. Salinas San Martin, J. Noronha-Hostler, H. Elfner, J. Hammelmann, R. Hirayama

Abstract: Recent lattice QCD results, comparing to a hadron resonance gas model, have shown the need for hundreds of particles in hadronic models. These extra particles influence both the equation of state and hadronic interactions within hadron transport models. Here, we introduce the PDG21+ particle list, which contains the most up-to-date database of particles and their properties. We then convert all pa… ▽ More Recent lattice QCD results, comparing to a hadron resonance gas model, have shown the need for hundreds of particles in hadronic models. These extra particles influence both the equation of state and hadronic interactions within hadron transport models. Here, we introduce the PDG21+ particle list, which contains the most up-to-date database of particles and their properties. We then convert all particles decays into 2 body decays so that they are compatible with SMASH in order to produce a more consistent description of a heavy-ion collision. △ Less

Submitted 19 July, 2022; originally announced July 2022.

Comments: Proceedings for the 37th Winter Workshop on Nuclear Dynamics

Temperature and net baryochemical potential dependence of $η/s$ in a hybrid approach

Authors: Niklas Götz, Hannah Elfner

Abstract: In this work, the qualitative impact of the net baryochemical potential dependence of the shear viscosity to entropy density ratio $η/s$ in hydrodynamical simulations is studied. The effect of a predicted non-constant $η/s$($μ_B$) is largely unexplored in hydrodynamic simulations. Previous studies focus only on a temperature dependence or even only a constant effective shear viscosity. This work a… ▽ More In this work, the qualitative impact of the net baryochemical potential dependence of the shear viscosity to entropy density ratio $η/s$ in hydrodynamical simulations is studied. The effect of a predicted non-constant $η/s$($μ_B$) is largely unexplored in hydrodynamic simulations. Previous studies focus only on a temperature dependence or even only a constant effective shear viscosity. This work addresses this issue by studying qualitatively the effect of a generalized $η/s$($T,μ_B$) in the hybrid approach SMASH-vHLLE, composed of the hadronic transport approach SMASH and the (3+1)d viscous hydrodynamic code vHLLE. In order to reduce the bias of the result on the equation of state used in the hydrodynamic part of the model, $η/s$ is parameterized directly in the energy density and net baryon number density. The parameterization takes into account the constraints of matching to the transport coefficients in the hadronic phase, as well as pQCD results. This work compares the impact of the density dependence for different system sizes and energies and compares the observables with experimental results in the RHIC - BES region $\sqrt{s_{NN}}$ =7.7 - 39.0 GeV, as the effect of this generalization is especially relevant for intermediate collision energies, for which the system is in equilibrium for a relevant amount of time, but the net baryochemical potential does not vanish. It is shown that the effect of an explicit net baryon number dependence on the elliptic flow is negligible and only relevant in the early stages of the collision. Additionally, we find that the proposed parameterization could be a good proxy for the shear viscosity in the non-equilibrium hadronic transport stage. △ Less

Submitted 30 October, 2022; v1 submitted 12 July, 2022; originally announced July 2022.

Comments: 11 pages, 12 figures

Effective spectral function of vector mesons via lifetime analysis

Authors: Renan Hirayama, Jan Staudenmaier, Hannah Elfner

Abstract: Effective spectral functions of the $ρ$ meson are reconstructed by considering the lifetimes inside different media using the hadronic transport SMASH (Simulating Many Accelerated Strongly-interacting Hadrons). Due to inelastic scatterings, resonance lifetimes are dynamically shortened (collisional broadening), even though the employed approach assumes vacuum resonance properties. Analyzing the… ▽ More Effective spectral functions of the $ρ$ meson are reconstructed by considering the lifetimes inside different media using the hadronic transport SMASH (Simulating Many Accelerated Strongly-interacting Hadrons). Due to inelastic scatterings, resonance lifetimes are dynamically shortened (collisional broadening), even though the employed approach assumes vacuum resonance properties. Analyzing the $ρ$ meson lifetimes allows to quantify an effective broadening of the decay width and spectral function, which is important in order to distinguish dynamical effects from additional genuine medium modifications to the spectral functions, indicating e.g. an onset of chiral symmetry restoration. The broadening of the spectral function in a thermalized system is shown to be consistent with other theoretical calculations. The effective $ρ$ meson spectral function is also presented for the dynamical evolution of heavy-ion collisions, finding a clear correlation of the broadening to system size, which is explained by an observed dependence of the width on the local hadron density. Furthermore, the difference in the results between the thermal system and full collision dynamics is explored, which may point to non-equilibrium effects. △ Less

Submitted 4 November, 2022; v1 submitted 30 June, 2022; originally announced June 2022.

Comments: 9 pages, 12 figures. Ancillary files divided by system

Inclusive jet and hadron suppression in a multistage approach

Authors: A. Kumar, Y. Tachibana, C. Sirimanna, G. Vujanovic, S. Cao, A. Majumder, Y. Chen, L. Du, R. Ehlers, D. Everett, W. Fan, Y. He, J. Mulligan, C. Park, A. Angerami, R. Arora, S. A. Bass, T. Dai, H. Elfner, R. J. Fries, C. Gale, F. Garza, M. Heffernan, U. Heinz, B. V. Jacak , et al. (35 additional authors not shown)

Abstract: We present a new study of jet interactions in the quark-gluon plasma created in high-energy heavy-ion collisions, using a multistage event generator within the JETSCAPE framework. We focus on medium-induced modifications in the rate of inclusive jets and high transverse momentum (high-$p_{\mathrm{T}}$) hadrons. Scattering-induced jet energy loss is calculated in two stages: A high virtuality stage… ▽ More We present a new study of jet interactions in the quark-gluon plasma created in high-energy heavy-ion collisions, using a multistage event generator within the JETSCAPE framework. We focus on medium-induced modifications in the rate of inclusive jets and high transverse momentum (high-$p_{\mathrm{T}}$) hadrons. Scattering-induced jet energy loss is calculated in two stages: A high virtuality stage based on the MATTER model, in which scattering of highly virtual partons modifies the vacuum radiation pattern, and a second stage at lower jet virtuality based on the LBT model, in which leading partons gain and lose virtuality by scattering and radiation. Coherence effects that reduce the medium-induced emission rate in the MATTER phase are also included. The TRENTo model is used for initial conditions, and the (2+1)dimensional VISHNU model is used for viscous hydrodynamic evolution. Jet interactions with the medium are modeled via 2-to-2 scattering with Debye screened potentials, in which the recoiling partons are tracked, hadronized, and included in the jet clustering. Holes left in the medium are also tracked and subtracted to conserve transverse momentum. Calculations of the nuclear modification factor ($R_{\mathrm{AA}}$) for inclusive jets and high-$p_{\mathrm{T}}$ hadrons are compared to experimental measurements at the BNL Relativistic Heavy Ion Collider (RHIC) and the CERN Large Hadron Collider (LHC). Within this framework, we find that with one extra parameter which codifies the transition between stages of jet modification -- along with the typical parameters such as the coupling in the medium, the start and stop criteria etc. -- we can describe these data at all energies for central and semicentral collisions without a rescaling of the jet transport coefficient $\hat{q}$. △ Less

Submitted 16 April, 2023; v1 submitted 3 April, 2022; originally announced April 2022.

Comments: 34 pages, 23 figures

Journal ref: Phys. Rev. C 107, no.3, 034911 (2023)

Role of bulk viscosity in deuteron production in ultrarelativistic nuclear collisions

Authors: D. Everett, D. Oliinychenko, M. Luzum, J. -F. Paquet, G. Vujanovic, S. A. Bass, L. Du, C. Gale, M. Heffernan, U. Heinz, L. Kasper, W. Ke, D. Liyanage, A. Majumder, A. Mankolli, C. Shen, D. Soeder, J. Velkovska, A. Angerami, R. Arora, S. Cao, Y. Chen, T. Dai, R. Ehlers, H. Elfner , et al. (31 additional authors not shown)

Abstract: We use a Bayesian-calibrated multistage viscous hydrodynamic model to explore deuteron yield, mean transverse momentum and flow observables in LHC Pb-Pb collisions. We explore theoretical uncertainty in the production of deuterons, including (i) the contribution of thermal deuterons, (ii) models for the subsequent formation of deuterons (hadronic transport vs coalescence) and (iii) the overall sen… ▽ More We use a Bayesian-calibrated multistage viscous hydrodynamic model to explore deuteron yield, mean transverse momentum and flow observables in LHC Pb-Pb collisions. We explore theoretical uncertainty in the production of deuterons, including (i) the contribution of thermal deuterons, (ii) models for the subsequent formation of deuterons (hadronic transport vs coalescence) and (iii) the overall sensitivity of the results to the hydrodynamic model -- in particular to bulk viscosity, which is often neglected in studies of deuteron production. Using physical parameters set by a comparison to only light hadron observables, we find good agreement with measurements of the mean transverse momentum $\langle p_T \rangle$ and elliptic flow $v_2$ of deuterons; however, tension is observed with experimental data for the deuteron multiplicity in central collisions. The results are found to be sensitive to each of the mentioned theoretical uncertainties, with a particular sensitivity to bulk viscosity, indicating that the latter is an important ingredient for an accurate treatment of deuteron production. △ Less

Submitted 15 March, 2022; originally announced March 2022.

Comments: 11 pages, 8 figures

Impact of hadronic interactions and conservation laws on cumulants of conserved charges in a dynamical model

Authors: Jan Hammelmann, Hannah Elfner

Abstract: Understanding the phase diagram of QCD by measuring fluctuations of conserved charges in heavy-ion collision is one of the main goals of the beam energy scan program at RHIC. Within this work, we calculate the role of hadronic interactions and momentum cuts on cumulants of conserved charges up to fourth order in a system in equilibrium within a hadronic transport approach (SMASH). In our model the… ▽ More Understanding the phase diagram of QCD by measuring fluctuations of conserved charges in heavy-ion collision is one of the main goals of the beam energy scan program at RHIC. Within this work, we calculate the role of hadronic interactions and momentum cuts on cumulants of conserved charges up to fourth order in a system in equilibrium within a hadronic transport approach (SMASH). In our model the net baryon, net charge and net strangeness is perfectly conserved on an event-by-event basis and the cumulants are calculated as a function of subvolume sizes and compared to analytic expectations. We find a modification of the kurtosis due to charge annihilation processes in systems with simplified degrees of freedom. Furthermore the result of the full SMASH hadron gas for the net baryon and net proton number fluctuations is presented for systems with zero and finite values of baryon chemical potential. Additionally we find that due to dynamical correlations the cumulants of the net baryon number cannot easily be recovered from the net protons. Finally the influence of deuteron cluster formation on the net proton and net baryon fluctuations in a simplified system is shown. This analysis is important to better understand the relation between measurements of fluctuations in heavy-ion collisions and theoretical calculation which are often performed in a grand canonical ensemble. △ Less

Submitted 23 February, 2022; originally announced February 2022.

Transport Model Comparison Studies of Intermediate-Energy Heavy-Ion Collisions

Authors: Hermann Wolter, Maria Colonna, Dan Cozma, Pawel Danielewicz, Che Ming Ko, Rohit Kumar, Akira Ono, ManYee Betty Tsang, Jun Xu, Ying-Xun Zhang, Elena Bratkovskaya, Zhao-Qing Feng, Theodoros Gaitanos, Arnaud Le Fèvre, Natsumi Ikeno, Youngman Kim, Swagata Mallik, Paolo Napolitani, Dmytro Oliinychenko, Tatsuhiko Ogawa, Massimo Papa, Jun Su, Rui Wang, Yong-Jia Wang, Janus Weil , et al. (27 additional authors not shown)

Abstract: Transport models are the main method to obtain physics information from low to relativistic-energy heavy-ion collisions. The Transport Model Evaluation Project (TMEP) has been pursued to test the robustness of transport model predictions in reaching consistent conclusions from the same type of physical model. Calculations under controlled conditions of physical input and set-up were performed with… ▽ More Transport models are the main method to obtain physics information from low to relativistic-energy heavy-ion collisions. The Transport Model Evaluation Project (TMEP) has been pursued to test the robustness of transport model predictions in reaching consistent conclusions from the same type of physical model. Calculations under controlled conditions of physical input and set-up were performed with various participating codes. These included both calculations of nuclear matter in a box with periodic boundary conditions, and more realistic calculations of heavy-ion collisions. In this intermediate review, we summarize and discuss the present status of the project. We also provide condensed descriptions of the 26 participating codes, which contributed to some part of the project. These include the major codes in use today. We review the main results of the studies completed so far. They show, that in box calculations the differences between the codes can be well understood and a convergence of the results can be reached. These studies also highlight the systematic differences between the two families of transport codes, known as BUU and QMD type codes. However, when the codes were compared in full heavy-ion collisions using different physical models, as recently for pion production, they still yielded substantially different results. This calls for further comparisons of heavy-ion collisions with controlled models and of box comparisons of important ingredients, like momentum-dependent fields, which are currently underway. We often indicate improved strategies in performing transport simulations and thus provide guidance to code developers. Results of transport simulations of heavy-ion collisions from a given code will have more significance if the code can be validated against benchmark calculations such as the ones summarized in this review. △ Less

Submitted 4 May, 2022; v1 submitted 14 February, 2022; originally announced February 2022.

Comments: 114 pages, 14 figures, 479 references, accepted for publication in Progress of Particle and Nuclear Phsics

Journal ref: Prog. Part. Nucl. Phys. 125 (2022) 103962

The applicability of hydrodynamics in heavy ion collisions at $\sqrt{s_{\rm NN}}$= 2.4-7.7 GeV

Authors: Gabriele Inghirami, Hannah Elfner

Abstract: To assess the degree of equilibration of the matter created in heavy-ion reactions at low to intermediate beam energies, a hadronic transport approach (SMASH) is employed. By using a coarse-graining method, we compute the energy momentum tensor of the system at fixed time steps and evaluate the degree of isotropy of the diagonal terms and the relative magnitude of the off-diagonal terms. This stud… ▽ More To assess the degree of equilibration of the matter created in heavy-ion reactions at low to intermediate beam energies, a hadronic transport approach (SMASH) is employed. By using a coarse-graining method, we compute the energy momentum tensor of the system at fixed time steps and evaluate the degree of isotropy of the diagonal terms and the relative magnitude of the off-diagonal terms. This study focuses mostly on Au+Au collisions in the energy range $\sqrt{s_{\rm NN}}$ = 2.4-7.7 GeV, but central collisions of lighter ions like C+C, Ar+KCl and Ag+Ag are considered as well. We find that the conditions concerning local equilibration for a hydrodynamic description are reasonably satisfied in a large portion of the system for a significant amount of time (several fm/c) when considering the average evolution of many events, yet they are rarely fulfilled on an event by event basis. This is relevant for the application of hybrid approaches at low beam energies as they are or will be reached by the HADES experiment at GSI, the future CBM experiment at FAIR as well as the beam energy scan program at RHIC. △ Less

Submitted 8 September, 2022; v1 submitted 15 January, 2022; originally announced January 2022.

Comments: Ancillary files with plots data and workflow provided. Version 2 corresponding to version accepted for publication by EPJC

CORSIKA 8 -- Contributions to the 37th International Cosmic Ray Conference in Berlin Germany (ICRC 2021)

Authors: Jean-Marco Alameddine, Johannes Albrecht, Jaime Alvarez-Muniz, Antonio Augusto Alves Jr, Luisa Arrabito, Dominik Baack, Konrad Bernlöhr, Marcus Bleicher, Johan Bregeon, Mathieu Carrere, Hans Dembinski, Hannah Elfner, Dominik Elsässer, Ralph Engel, Hu Fan, Anatoli Fedynitch, Dieter Heck, Tim Huege, Karl-Heinz Kampert, Nikolaos Karastathis, Lukas Nellen, Maximilian Nöthe, David Parello, Tanguy Pierog, Maria Pokrandt , et al. (17 additional authors not shown)

Abstract: Compilation of the six contributions to the ICRC conference 2021 by the CORSIKA 8 Collaboration. The status of the project is illustrated. In particular, the secondary hadron as well as the electromagnetic cascades are being validated individually, and current results are reviewed. A novel framework for radio emission simulations is presented, which is designed given the modular nature of CORSIKA… ▽ More Compilation of the six contributions to the ICRC conference 2021 by the CORSIKA 8 Collaboration. The status of the project is illustrated. In particular, the secondary hadron as well as the electromagnetic cascades are being validated individually, and current results are reviewed. A novel framework for radio emission simulations is presented, which is designed given the modular nature of CORSIKA 8 to support, both, the CoREAS as well as the ZHS formalism. At the same time, first Cherenkov emission calculations are shown which are based on CORSIKA 8 coupled with a GPU Cherenkov emission code. Finally, a new powerful feature of CORSIKA 8 is illustrated, where the entire genealogy of air shower particles can be studied in all details. △ Less

Submitted 22 December, 2021; originally announced December 2021.

Comments: ICRC 2021 conference, 6 proceedings, 51 pages

Particle production in a hybrid approach for a beam energy scan of Au+Au/Pb+Pb collisions between $\sqrt{s_\mathrm{NN}}$ = 4.3 GeV and $\sqrt{s_\mathrm{NN}}$ = 200.0 GeV

Authors: Anna Schäfer, Iurii Karpenko, Xiang-Yu Wu, Jan Hammelmann, Hannah Elfner

Abstract: Heavy-ion collisions at varying collision energies provide access to different regions of the QCD phase diagram. In particular collisions at intermediate energies are promising candidates to experimentally identify the postulated first order phase transition and critical end point. While heavy-ion collisions at low and high collision energies are theoretically well described by transport approac… ▽ More Heavy-ion collisions at varying collision energies provide access to different regions of the QCD phase diagram. In particular collisions at intermediate energies are promising candidates to experimentally identify the postulated first order phase transition and critical end point. While heavy-ion collisions at low and high collision energies are theoretically well described by transport approaches and hydrodynamics+transport hybrid approaches, respectively, intermediate energy collisions remain a challenge. In this work, a modular hybrid approach, the SMASH-vHLLE-hybrid coupling 3+1D viscous hydrodynamics (vHLLE) to hadronic transport (SMASH), is introduced. It is validated and subsequently applied in Au+Au/Pb+Pb collisions between $\sqrt{s_\mathrm{NN}}$ = 4.3 GeV and $\sqrt{s_\mathrm{NN}}$ = 200.0 GeV to study the rapidity and transverse mass distributions of identified particles as well as excitation functions for $\mathrm{dN}/\mathrm{d}y|_{y = 0}$ and $\langle p_\mathrm{T} \rangle$. A good agreement with experimental measurements is obtained, including the baryon stopping dynamics. The transition from a Gaussian rapidity spectrum of protons at lower energies to the double-hump structure at high energies is reproduced. The centrality and energy dependence of charged particle $v_2$ is also described reasonably well. This work serves as a basis for further studies, e.g. systematic investigations of different equations of state or transport coefficients. △ Less

Submitted 22 January, 2023; v1 submitted 16 December, 2021; originally announced December 2021.

Out-of-Equilibrium Photon Production in the Late Stages of Relativistic Heavy-Ion Collisions

Authors: Anna Schäfer, Oscar Garcia-Montero, Jean-François Paquet, Hannah Elfner, Charles Gale

Abstract: In this work, we assess the importance of non-equilibrium dynamics in the production of photons from the late stages of relativistic heavy-ion collisions. The p$_\mathrm{T}$-differential spectra and $v_2$ of photons from the late hadronic stage are computed within a non-equilibrium hadron transport approach, and compared to the results of a local equilibrium evolution using ideal relativistic hydr… ▽ More In this work, we assess the importance of non-equilibrium dynamics in the production of photons from the late stages of relativistic heavy-ion collisions. The p$_\mathrm{T}$-differential spectra and $v_2$ of photons from the late hadronic stage are computed within a non-equilibrium hadron transport approach, and compared to the results of a local equilibrium evolution using ideal relativistic hydrodynamics. It is found that non-equilibrium dynamics enhance the late-stage photon production at low p$_\mathrm{T}$ and decreases it at higher p$_\mathrm{T}$ compared to the estimate from hydrodynamics. This same comparison points to a significant increase in the momentum anisotropies of these photons due to non-equilibrium dynamics. Once combined with photons produced above the particlization temperature in the hydrodynamics evolution, the differences between the two approaches appear modest in what concerns the p$_\mathrm{T}$ differential spectra, but are clearly noticeable at low p$_\mathrm{T}$ for the elliptic flow: non-equilibrium dynamics enhance the photon $v_2$ below p$_\mathrm{T}$ $\approx 1.4$ GeV. △ Less

Submitted 14 June, 2022; v1 submitted 26 November, 2021; originally announced November 2021.

Comments: Data is available in the anc directory

Conservation laws in a novel hybrid approach

Authors: Anna Schäfer, Iurii Karpenko, Hannah Elfner

Abstract: Heavy-ion collisions covering a wide range of collision energies provide a vast amount of observables characterizing the properties of strongly-interacting matter. In particular collisions towards the high baryon-density regime of the QCD phase-diagram have become of interest to study the postulated first order phase transition and to locate a possible critical end point. In this work, the SMASH-v… ▽ More Heavy-ion collisions covering a wide range of collision energies provide a vast amount of observables characterizing the properties of strongly-interacting matter. In particular collisions towards the high baryon-density regime of the QCD phase-diagram have become of interest to study the postulated first order phase transition and to locate a possible critical end point. In this work, the SMASH-vHLLE-hybrid is presented as a novel hybrid model to theoretically describe such heavy-ion collisions. In addition, the SMASH hadron resonance gas equation of state is introduced. The accuracy of the latter is shown to be of fundamental importance in order to conserve energy, baryon number and electric charge throughout the different stages of the hybrid model. Furthermore, the impact of an inaccurate equation of state on final state observables is discussed. This work constitutes a first validation of the SMASH-vHLLE-hybrid in terms of conservation laws and excitation functions. It is expected to be applied to a broader range of observables in the future. △ Less

Submitted 17 September, 2021; originally announced September 2021.

Comments: Proceedings CPOD 2021

The BEST framework for the search for the QCD critical point and the chiral magnetic effect

Authors: Xin An, Marcus Bluhm, Lipei Du, Gerald V. Dunne, Hannah Elfner, Charles Gale, Joaquin Grefa, Ulrich Heinz, Anping Huang, Jamie M. Karthein, Dmitri E. Kharzeev, Volker Koch, Jinfeng Liao, Shiyong Li, Mauricio Martinez, Michael McNelis, Debora Mroczek, Swagato Mukherjee, Marlene Nahrgang, Angel R. Nava Acuna, Jacquelyn Noronha-Hostler, Dmytro Oliinychenko, Paolo Parotto, Israel Portillo, Maneesha Sushama Pradeep , et al. (18 additional authors not shown)

Abstract: The Beam Energy Scan Theory (BEST) Collaboration was formed with the goal of providing a theoretical framework for analyzing data from the Beam Energy Scan (BES) program at the relativistic heavy ion collider (RHIC) at Brookhaven National Laboratory. The physics goal of the BES program is the search for a conjectured QCD critical point as well as for manifestations of the chiral magnetic effect. W… ▽ More The Beam Energy Scan Theory (BEST) Collaboration was formed with the goal of providing a theoretical framework for analyzing data from the Beam Energy Scan (BES) program at the relativistic heavy ion collider (RHIC) at Brookhaven National Laboratory. The physics goal of the BES program is the search for a conjectured QCD critical point as well as for manifestations of the chiral magnetic effect. We describe progress that has been made over the previous five years. This includes studies of the equation of state and equilibrium susceptibilities, the development of suitable initial state models, progress in constructing a hydrodynamic framework that includes fluctuations and anomalous transport effects, as well as the development of freezeout prescriptions and hadronic transport models. Finally, we address the challenge of integrating these components into a complete analysis framework. This document describes the collective effort of the BEST Collaboration and its collaborators around the world. △ Less

Submitted 22 November, 2021; v1 submitted 31 August, 2021; originally announced August 2021.

Comments: 103 pages

Journal ref: Nuclear Physics A, Volume 1017, January 2022, 122343

The role of proton-antiproton regeneration in the late stages of heavy-ion collisions

Authors: Oscar Garcia-Montero, Jan Staudenmaier, Anna Schäfer, Juan M. Torres-Rincon, Hannah Elfner

Abstract: We investigate the long-standing question of the effect of proton-antiproton annihilation on the (anti-)proton yield, while respecting detailed balance for the 5-body back-reaction for the first time in a full microscopic description of the late stages of heavy-ion collisions. This is achieved by employing a stochastic collision criterion in a hadronic transport approach (SMASH), which allows to t… ▽ More We investigate the long-standing question of the effect of proton-antiproton annihilation on the (anti-)proton yield, while respecting detailed balance for the 5-body back-reaction for the first time in a full microscopic description of the late stages of heavy-ion collisions. This is achieved by employing a stochastic collision criterion in a hadronic transport approach (SMASH), which allows to treat arbitrary multi-particle reactions. It is used to account for the regeneration of (anti-)protons via $5π\rightarrow p\bar{p}$. Our results show that a back-reaction happens for a fraction of 15-20\% of all annihilations. Within a viscous hybrid approach Au+Au/Pb+Pb collisions from $\sqrt{s_{NN}}=17.3$ GeV$-5.02$ TeV are investigated and the quoted fraction is independent of the beam energy or centrality of the collision. Taking the back-reaction into account results in regeneration of half of the (anti-)proton yield that is lost due to annihilations at midrapidity. We also find that, concerning the multiplicities, treating the back-reaction as a chain of 2-body reactions is equivalent to a single 5-to-2 reaction. △ Less

Submitted 19 July, 2021; originally announced July 2021.

Comments: 6 pages, 4 figures

"Smashing more than two": Deuteron production in relativistic heavy ion collisions via stochastic multi-particle reactions

Authors: Jan Staudenmaier, Dmytro Oliinychenko, Juan M. Torres-Rincon, Hannah Elfner

Abstract: We study the deuteron production via the deuteron pion and nucleon catalysis reactions, $πp n \leftrightarrow πd$ and $N p n \leftrightarrow N d$, by employing stochastic multi-particle reactions in the hadronic transport approach SMASH for the first time. This is an improvement compared to previous studies, which introduced an artificial fake resonance $d'$ to simulate these… ▽ More We study the deuteron production via the deuteron pion and nucleon catalysis reactions, $πp n \leftrightarrow πd$ and $N p n \leftrightarrow N d$, by employing stochastic multi-particle reactions in the hadronic transport approach SMASH for the first time. This is an improvement compared to previous studies, which introduced an artificial fake resonance $d'$ to simulate these $3 \leftrightarrow 2$ reactions as a chain of $2\leftrightarrow 2$ reactions. The derivation of the stochastic criterion for multi-particle reactions is presented in a comprehensive fashion and its implementation is tested against an analytic expression for the scattering rate and the equilibrating particle yields in box calculations. We then study Au + Au collisions at $\sqrt{s_{\mathrm{NN}}} = 7.7$ GeV, where we find that multi-particle collisions substantially reduce the time required for deuterons to reach partial chemical equilibrium with nucleons. Subsequently, the final yield of $d$ is practically independent from the number of $d$ at particlization, confirming the results of previous studies. The mean transverse momentum and the integrated elliptic flow as a function of centrality are rather insensitive to the exact realization of the $2\leftrightarrow 3$ reactions. △ Less

Submitted 26 November, 2021; v1 submitted 27 June, 2021; originally announced June 2021.

Comments: 14 pages, 10 figures; included text files with the data to reproduce the figures; replaced with published version; additional appendix about grid setup

Journal ref: Phys. Rev. C 104, 034908 (2021)