• We report on broadly based systematic investigations of the modeling components for open heavy-flavor diffusion and energy loss in strongly interacting matter in their application to heavy-flavor observables in high-energy heavy-ion collisions, conducted within an EMMI Rapid Reaction Task Force framework. Initial spectra including cold-nuclear-matter effects, a wide variety of space-time evolution models, heavy-flavor transport coefficients, and hadronization mechanisms are scrutinized in an effort to quantify pertinent uncertainties in the calculations of nuclear modification factors and elliptic flow of open heavy-flavor particles in nuclear collisions. We develop procedures for error assessments and criteria for common model components to improve quantitative estimates for the (low-momentum) heavy-flavor diffusion coefficient as a long-wavelength characteristic of QCD matter as a function of temperature, and for energy loss coefficients of high-momentum heavy-flavor particles.
  • The Future Circular Collider (FCC) Study is aimed at assessing the physics potential and the technical feasibility of a new collider with centre-of-mass energies, in the hadron-hadron collision mode, seven times larger than the nominal LHC energies. Operating such machine with heavy ions is an option that is being considered in the accelerator design studies. It would provide, for example, Pb-Pb and p-Pb collisions at sqrt{s_NN} = 39 and 63 TeV, respectively, per nucleon-nucleon collision, with integrated luminosities above 30 nb^-1 per month for Pb-Pb. This is a report by the working group on heavy-ion physics of the FCC Study. First ideas on the physics opportunities with heavy ions at the FCC are presented, covering the physics of the Quark-Gluon Plasma, of gluon saturation, of photon-induced collisions, as well as connections with other fields of high-energy physics.
  • A. Dainese, E. Scomparin, G. Usai, P. Antonioli, R. Arnaldi, A. Beraudo, E. Bruna, G.E. Bruno, S. Bufalino, P. Di Nezza, M.P. Lombardo, R. Nania, F.Noferini, C. Oppedisano, S. Piano, F. Prino, A. Rossi, M. Agnello, W.M. Alberico, B. Alessandro, A. Alici, G. Andronico, F. Antinori, S. Arcelli, A. Badala, A.M. Barbano, R. Barbera, F. Barile, M. Basile, F. Becattini, C. Bedda, F. Bellini, S. Beole, L. Bianchi, C. Bianchin, C. Bonati, F. Bossu, E. Botta, D. Caffarri, P. Camerini, F. Carnesecchi, E. Casula, P. Cerello, C. Cicalo, M.L. Cifarelli, F. Cindolo, F. Colamaria, D. Colella, M. Colocci, Y. Corrales Morales, P. Cortese, A. De Caro, G. De Cataldo, A. De Falco, D. De Gruttola, M. D'Elia, N. De Marco, S. De Pasquale, D. Di Bari, D. Elia, A. Fantoni, A. Feliciello, A. Ferretti, A. Festanti, F. Fionda, G. Fiorenza, E. Fragiacomo, G.G. Fronze, M. Fusco Girard, M. Gagliardi, M. Gallio, K. Garg, P. Giubellino, V. Greco, E. Grossi, B. Guerzoni, D. Hatzifotiadou, E. Incani, G.M. Innocenti, N. Jacazio, S. Kumar Das, P. La Rocca, R. Lea, L. Leardini, M. Leoncino, M. Lunardon, G. Luparello, V. Mantovani Sarti, V. Manzari, M. Marchisone, G.V. Margagliotti, M. Masera, A. Masoni, A. Mastroserio, M. Mazzilli, M.A. Mazzoni, E. Meninno, M. Mesiti, L. Milano, S. Moretto, V. Muccifora, E. Nappi, M. Nardi, M. Nicassio, P. Pagano, G.S. Pappalardo, C. Pastore, B. Paul, C. Petta, O. Pinazza, S. Plumari, R. Preghenella, M. Puccio, G. Puddu, L. Ramello, C. Ratti, I. Ravasenga, F. Riggi, F. Ronchetti, A. Rucci, M. Ruggieri, R. Rui, S. Sakai, E. Scapparone, F. Scardina, F. Scarlassara, G. Scioli, S. Siddhanta, M. Sitta, F. Soramel, M. Suljic, C. Terrevoli, S. Trogolo, G. Trombetta, R. Turrisi, E. Vercellin, G. Vino, T. Virgili, G. Volpe, M.C.S. Williams, C. Zampolli
    Feb. 12, 2016 nucl-ex, nucl-th
    This document was prepared by the community that is active in Italy, within INFN (Istituto Nazionale di Fisica Nucleare), in the field of ultra-relativistic heavy-ion collisions. The experimental study of the phase diagram of strongly-interacting matter and of the Quark-Gluon Plasma (QGP) deconfined state will proceed, in the next 10-15 years, along two directions: the high-energy regime at RHIC and at the LHC, and the low-energy regime at FAIR, NICA, SPS and RHIC. The Italian community is strongly involved in the present and future programme of the ALICE experiment, the upgrade of which will open, in the 2020s, a new phase of high-precision characterisation of the QGP properties at the LHC. As a complement of this main activity, there is a growing interest in a possible future experiment at the SPS, which would target the search for the onset of deconfinement using dimuon measurements. On a longer timescale, the community looks with interest at the ongoing studies and discussions on a possible fixed-target programme using the LHC ion beams and on the Future Circular Collider.
  • Predictions for heavy-flavour production in relativistic heavy-ion experiments provided by the POWLANG transport setup, including now also an in-medium hadronization model, are displayed, After showing some representative findings for the Au-Au and Pb-Pb cases, a special focus will be devoted to the results obtained in the small systems formed in proton(deuteron)-nucleus collisions, where recent experimental data suggest the possible formation of a medium featuring a collective behaviour.
  • This report reviews the study of open heavy-flavour and quarkonium production in high-energy hadronic collisions, as tools to investigate fundamental aspects of Quantum Chromodynamics, from the proton and nucleus structure at high energy to deconfinement and the properties of the Quark-Gluon Plasma. Emphasis is given to the lessons learnt from LHC Run 1 results, which are reviewed in a global picture with the results from SPS and RHIC at lower energies, as well as to the questions to be addressed in the future. The report covers heavy flavour and quarkonium production in proton-proton, proton-nucleus and nucleus-nucleus collisions. This includes discussion of the effects of hot and cold strongly interacting matter, quarkonium photo-production in nucleus-nucleus collisions and perspectives on the study of heavy flavour and quarkonium with upgrades of existing experiments and new experiments. The report results from the activity of the SaporeGravis network of the I3 Hadron Physics programme of the European Union 7th Framework Programme.
  • We present recent results for heavy-flavor observables in nucleus-nucleus collisions at LHC energies, obtained with the POWLANG transport setup. The initial creation of c-cbar and b-bbar pairs is simulated with a perturbative QCD approach (POWHEG+PYTHIA); their propagation in the medium (created in the nucleus-nucleus or in proton-nucleus collision) is studied with the relativistic Langevin equation, here solved using weak-coupling transport coefficients. Successively, the heavy quarks hadronize in the medium. We compute the nuclear modification factor and the elliptic flow parameter of the final D mesons both in nucleus-nucleus and in (for the first time, in the POWLANG setup) proton-nucleus collisions and compare our results to experimental data.
  • We present a quantitative study of vorticity formation in peripheral ultrarelativistic heavy ion collisions at sqrt(s)NN = 200 GeV by using the ECHO-QGP numerical code, implementing relativistic dissipative hydrodynamics in the causal Israel-Stewart framework in 3+1 dimensions with an initial Bjorken flow profile. We consider and discuss different definitions of vorticity which are relevant in relativistic hydrodynamics. After demonstrating the excellent capabilities of our code, which proves to be able to reproduce Gubser flow up to 8 fm/c, we show that, with the initial conditions needed to reproduce the measured directed flow in peripheral collisions corresponding to an average impact parameter b=11.6 fm and with the Bjorken flow profile for a viscous Quark Gluon Plasma with \eta/s=0.1 fixed, a vorticity of the order of some 10^{-2} c/fm can develop at freezeout. The ensuing polarization of Lambda baryons does not exceed 1.4% at midrapidity. We show that the amount of developed directed flow is sensitive to both the initial angular momentum of the plasma and its viscosity.
  • We present results for the quenching, elliptic flow and azimuthal correlations of heavy flavour particles in high-energy nucleus-nucleus collisions obtained through the POWLANG transport setup, developed in the past to study the propagation of heavy quarks in the Quark-Gluon Plasma and here extended to include a modeling of their hadronization in the presence of a medium. Hadronization is described as occurring via the fragmentation of strings with endpoints given by the heavy (anti-)quark Q(Qbar) and a thermal parton qbar(q) from the medium. The flow of the light quarks is shown to affect significantly the R_AA and v_2 of the final D mesons, leading to a better agreement with the experimental data. The approach allows also predictions for the angular correlation between heavy-flavour hadrons (or their decay electrons) and the charged particles produced in the fragmentation of the heavy-quark strings.
  • We present results for the quenching, elliptic flow and azimuthal correlations of heavy flavour particles in high-energy nucleus-nucleus collisions obtained through the POWLANG transport setup, developed in the past to study the propagation of heavy quarks in the Quark-Gluon Plasma and here extended to include a modeling of their hadronization in the presence of a medium. Hadronization is described as occurring via the fragmentation of strings with endpoints given by the heavy (anti-)quark Q(Qbar) and a thermal parton $qbar(q)$ from the medium. The flow of the light quarks is shown to affect significantly the R_AA} and v_2 of the final D mesons, leading to a better agreement with the experimental data.
  • We present realistic equations of state for QCD matter at vanishing net-baryon density which embed recent lattice QCD results at high temperatures combined with a hadron resonance gas model in the low-temperature, confined phase. In the latter, we allow an implementation of partial chemical equilibrium, in which particle ratios are fixed at the chemical freeze-out, so that a description closer to the experimental situation is possible. Given the present uncertainty in the determination of the chemical freeze-out temperature from first-principle lattice QCD calculations, we consider different values within the expected range. The corresponding equations of state can be applied in the hydrodynamic modeling of relativistic heavy-ion collisions at the LHC and at the highest RHIC beam energies. Suitable parametrizations of our results as functions of the energy density are also provided.
  • A snapshot of the results for heavy-flavour observables in heavy-ion (AA) collisions at RHIC and LHC obtained with our transport calculations is displayed. The initial charm and beauty production is simulated through pQCD tools (POWHEG+PYTHIA) and is validated through the comparison with data from pp collisions. The propagation of c and b quarks in the medium formed in heavy-ion collisions is studied through a transport setup based on the relativistic Langevin equation. With respect to past works we perform a more systematic study, providing results with different choices of transport coefficients, either from weak-coupling calculations or from lattice-QCD simulations. Our findings are compared to a rich set of experimental data (D-mesons, non-photonic electrons, non-prompt J/\psi's) which have meanwhile become accessible.
  • We present ECHO-QGP, a numerical code for $(3+1)$-dimensional relativistic viscous hydrodynamics designed for the modeling of the space-time evolution of the matter created in high energy nuclear collisions. The code has been built on top of the \emph{Eulerian Conservative High-Order} astrophysical code for general relativistic magneto-hydrodynamics [\emph{Del Zanna et al., Astron. Astrophys. 473, 11, 2007}] and here it has been upgraded to handle the physics of the Quark-Gluon Plasma. ECHO-QGP features second-order treatment of causal relativistic viscosity effects in both Minkowskian or Bjorken coordinates; partial or complete chemical equilibrium of hadronic species before kinetic freeze-out; initial conditions based on the optical Glauber model, including a Monte-Carlo routine for event-by-event fluctuating initial conditions; a freeze-out procedure based on the Cooper-Frye prescription. The code is extensively validated against several test problems and results always appear accurate, as guaranteed by the combination of the conservative (shock-capturing) approach and the high-order methods employed. ECHO-QGP can be extended to include evolution of the electromagnetic fields coupled to the plasma.
  • The quenching of jets (and high-pT particle spectra) observed in heavy-ion collisions is interpreted as due to the energy lost by hard partons crossing the Quark Gluon Plasma. Here we review recent efforts to include in its modeling important qualitative features of QCD, like the correlations in multiple gluon emissions and the color-flow pattern in parton branchings. In particular, the modification of color connections among the partons of a shower developing in the presence of a medium is a generic occurrence accompanying parton energy-loss. We show how this effect can leave its fingerprints at the hadronization stage, leading by itself to a softening of hadron spectra and to an enhanced production of soft particles in jet-fragmentation.
  • A refined version of a multi-step calculation of heavy-flavor observables in pp and AA collisions has been developed, based on pQCD at NLO accuracy followed by parton shower evolution to describe heavy-quark production and on the relativistic Langevin equation to describe their stochastic evolution in the QCD plasma. Then, hadronization is modeled through an implementation of fragmentation functions based on pQCD and constrained by $e^{+}e^{-}$ collider data. Results of our calculations can be compared with recent measurements performed at the LHC in Pb--Pb collisions at $\sqrt{s_{NN}}$=2.76 TeV: nuclear modification factor $R_{AA}$ of the $p_{T}$ spectra at mid-rapidity of heavy-flavor decay electrons and of exclusively reconstructed open-charm mesons at different centralities, as well as their elliptic-flow $v_{2}$($p_{T}$) in semi-central collisions. To test the validity of our setup for such studies, its predictions are also checked against the $p_{T}$ spectra measured in pp collisions at $\sqrt{s}$=7 TeV and 2.76 TeV.
  • Multiple interactions between parton showers and the surrounding QCD matter are expected to underlie the strong medium-modifications of jet observables in ultra-relativistic heavy ion collisions at RHIC and at the LHC. Here, we note that such jet-medium interactions alter generically and characteristically the color correlations in the parton shower. We characterize these effects in a color-differential calculation of the medium-induced gluon radiation spectrum to first and second order in opacity. By interfacing simple branching histories of medium-modified color flow with the Lund hadronization model, we analyze how the medium modification of color correlations can affect the distribution of hadronic fragments in jets. Importantly, we observe that jet-medium interactions give rise to the medium-induced color decoherence of gluons from the parton shower. Since hadronization respects color flow and since each color singlet in a parton shower is hadronized separately, this medium-induced color decoherence leaves characteristic signatures in the jet fragmentation pattern. In particular, it can contribute to the quenching of leading hadron spectra. Moreover, it can increase strongly the yield of soft hadronic fragments from a jet, while the distribution of more energetic hadrons follows naturally the shape of a vacuum-like fragmentation pattern of lower total energy.
  • Medium-induced parton energy loss, resulting from gluon exchanges between the QCD matter and partonic projectiles, is expected to underly the strong suppression of jets and high-$p_T$ hadron spectra observed in ultra-relativistic heavy ion collisions. Here, we present the first color-differential calculation of parton energy loss. We find that color exchange between medium and projectile enhances the invariant mass of energetic color singlet clusters in the parton shower by a parametrically large factor proportional to the square root of the projectile energy. This effect is seen in more than half of the most energetic color-singlet fragments of medium-modified parton branchings. Applying a standard cluster hadronization model, we find that it leads to a characteristic additional softening of hadronic spectra. A fair description of the nuclear modification factor measured at the LHC may then be obtained for relatively low momentum transfers from the medium.
  • The suppression of high-pT single-hadron spectra in heavy-ion collisions is usually interpreted as due to parton energy-loss of high-momentum quarks and gluons propagating in the plasma. Here, we discuss to what extent this partonic picture must be complemented by a picture of medium-modified hadronization. In particular, we show how color-exchange with the medium modifies the properties of color singlet-clusters arising from the parton branchings, producing a softening of the hadron spectra.
  • We present a study of the heavy-flavor dynamics in nucleus-nucleus collisions. The initial (hard) production of c and b quarks is taken from NLO pQCD predictions. The presence of a hot medium (a Quark Gluon Plasma described by hydrodynamics) affects the final spectra of open-charm (beauty) hadrons and their decay electrons with respect to what found in pp collisions. The propagation of c and b quarks in the plasma is based on a picture of multiple uncorrelated random collisions, described by a relativistic Langevin equation. A microscopic evaluation of the transport coefficients is provided within a pQCD approach (with proper resummation of medium effects). Results for the final spectra of heavy-flavor hadrons and decay-electrons are given, with particular emphasis on R_AA and v_2.
  • We study the stochastic dynamics of c and b quarks in the hot plasma produced in nucleus-nucleus collisions at RHIC and LHC, providing results for the nuclear modification factor R_AA and the elliptic flow coefficient v_2 of the single-electron spectra arising from their semi-leptonic decays. The initial QQbar pairs are generated using the POWHEG code, implementing pQCD at NLO. For the propagation in the plasma we develop a relativistic Langevin equation (solved in a medium described by hydrodynamics) whose transport coefficients are evaluated through a first-principle calculation. Finally, at T_c, the heavy quarks are made hadronize and decay into electrons: the resulting spectra are then compared with RHIC results. Predictions for LHC are also attempted.
  • We study the stochastic dynamics of c and b quarks, produced in hard initial processes, in the hot medium created after the collision of two relativistic heavy ions. This is done through the numerical solution of the relativistic Langevin equation. The latter requires the knowledge of the friction and diffusion coefficients, whose microscopic evaluation is performed treating separately the contribution of soft and hard collisions. The evolution of the background medium is described by ideal/viscous hydrodynamics. Below the critical temperature the heavy quarks are converted into hadrons, whose semileptonic decays provide single-electron spectra to be compared with the current experimental data measured at RHIC. We focus on the nuclear modification factor R_AA and on the elliptic-flow coefficient v_2, getting, for sufficiently large p_T, a reasonable agreement.
  • The stochastic dynamics of heavy quarks in the fireball produced in heavy-ion collisions is followed through numerical simulations based on the Langevin equation. The modification of the final p_T spectra (R_AA) of c and b quarks, hadrons and single-electrons with respect to pp collisions is studied. The transport coefficients are evaluated treating separately the contribution of soft and hard collisions. The initial heavy-quark spectra are generated according to NLO-pQCD, accounting for nuclear effects through recent nPDFs. The evolution of the medium is obtained from the output of two hydro-codes (ideal and viscous). The heavy-quark fragmentation into hadrons and their final semileptonic decays are implemented according to up to date experimental data. A comparison with RHIC data for non-photonic electron spectra is given.
  • We propose a model for the propagation of a heavy-quark in a hot plasma, to be viewed as a first step towards a full description of the dynamics of heavy quark systems in a quark-gluon plasma, including bound state formation. The heavy quark is treated as a non relativistic particle interacting with a fluctuating field, whose correlator is determined by a hard thermal loop approximation. This approximation, which concerns only the medium in which the heavy quark propagates, is the only one that is made, and it can be improved. The dynamics of the heavy quark is given exactly by a quantum mechanical path integral that is calculated in this paper in the Euclidean space-time using numerical Monte Carlo techniques. The spectral function of the heavy quark in the medium is then reconstructed using a Maximum Entropy Method. The path integral is also evaluated exactly in the case where the mass of the heavy quark is infinite; one then recovers known results concerning the complex optical potential that controls the long time behavior of the heavy quark. The heavy quark correlator and its spectral function is also calculated semi-analytically at the one-loop order, which allows for a detailed description of the coupling between the heavy quark and the plasma collective modes.
  • Quark Gluon Plasma transport coefficients for heavy quarks and quark-antiquark pairs are computed through an extension of the results obtained for a hot QED plasma by describing the heavy-quark propagation in the eikonal approximation and by weighting the gauge field configurations with the Hard Thermal Loop effective action. It is shown that such a model allows to correctly reproduce, at leading logarithmic accuracy, the results obtained by other independent approaches. The results are then inserted into a relativistic Langevin equation allowing to follow the evolution of the heavy-quark momentum spectra. Our numerical findings are also compared with the ones obtained in a strongly-coupled scenario, namely with the transport coefficients predicted (though with some limitations and ambiguities) by the AdS/CFT correspondence.
  • The heavy-quark spectral function in a hot plasma is reconstructed from the corresponding euclidean propagator. The latter is evaluated through a path-integral simulation. A weak-coupling calculation is also performed, allowing to interpret the qualitative behavior of the spectral function in terms of quite general physical processes.
  • We investigate the behavior of a pair of heavy fermions, denoted by $Q$ and $\bar{Q}$, in a hot/dense medium. Although we have in mind the situation where $Q$ and $\bar{Q}$ denote heavy quarks, our treatment will be limited to simplified models, which bear only some general similarities with QCD. We study in particular the limiting case where the mass of the heavy fermions is infinite. Then a number of results can be derived exactly: a Schr\"odinger equation can be established for the correlator of the heavy quarks; the interaction effects exponentiate, leading to a simple instantaneous effective potential for this Schr\"odinger equation. We consider simple models for the medium in which the $Q\bar Q$ pair propagates. In the case where the medium is a plasma of photons and light charged fermions, an imaginary part develops in this effective potential. We discuss the physical interpretation of this imaginary part in terms of the collisions between the heavy particles and the light fermions of the medium; the same collisions also determine the damping rate of the heavy fermions. Finally we study the connection between the real-time propagator of the heavy fermion pair and its Euclidean counterpart, and show that the real part of the potential entering the Schr\"odinger equation for the real-time propagator is the free energy calculated in the imaginary-time formalism.