
We report on broadly based systematic investigations of the modeling
components for open heavyflavor diffusion and energy loss in strongly
interacting matter in their application to heavyflavor observables in
highenergy heavyion collisions, conducted within an EMMI Rapid Reaction Task
Force framework. Initial spectra including coldnuclearmatter effects, a wide
variety of spacetime evolution models, heavyflavor 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 heavyflavor particles in nuclear collisions. We develop
procedures for error assessments and criteria for common model components to
improve quantitative estimates for the (lowmomentum) heavyflavor diffusion
coefficient as a longwavelength characteristic of QCD matter as a function of
temperature, and for energy loss coefficients of highmomentum heavyflavor
particles.

The Future Circular Collider (FCC) Study is aimed at assessing the physics
potential and the technical feasibility of a new collider with centreofmass
energies, in the hadronhadron 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, PbPb and pPb collisions at sqrt{s_NN} = 39 and 63 TeV, respectively,
per nucleonnucleon collision, with integrated luminosities above 30 nb^1 per
month for PbPb. This is a report by the working group on heavyion physics of
the FCC Study. First ideas on the physics opportunities with heavy ions at the
FCC are presented, covering the physics of the QuarkGluon Plasma, of gluon
saturation, of photoninduced collisions, as well as connections with other
fields of highenergy physics.

This document was prepared by the community that is active in Italy, within
INFN (Istituto Nazionale di Fisica Nucleare), in the field of
ultrarelativistic heavyion collisions. The experimental study of the phase
diagram of stronglyinteracting matter and of the QuarkGluon Plasma (QGP)
deconfined state will proceed, in the next 1015 years, along two directions:
the highenergy regime at RHIC and at the LHC, and the lowenergy 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 highprecision 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 fixedtarget programme using the LHC ion beams and on
the Future Circular Collider.

Predictions for heavyflavour production in relativistic heavyion
experiments provided by the POWLANG transport setup, including now also an
inmedium hadronization model, are displayed, After showing some representative
findings for the AuAu and PbPb 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 heavyflavour and quarkonium production
in highenergy 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 QuarkGluon 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 protonproton, protonnucleus and nucleusnucleus
collisions. This includes discussion of the effects of hot and cold strongly
interacting matter, quarkonium photoproduction in nucleusnucleus 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 heavyflavor observables in nucleusnucleus
collisions at LHC energies, obtained with the POWLANG transport setup. The
initial creation of ccbar and bbbar pairs is simulated with a perturbative
QCD approach (POWHEG+PYTHIA); their propagation in the medium (created in the
nucleusnucleus or in protonnucleus collision) is studied with the
relativistic Langevin equation, here solved using weakcoupling 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 nucleusnucleus and in (for the first time, in the
POWLANG setup) protonnucleus 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
ECHOQGP numerical code, implementing relativistic dissipative hydrodynamics in
the causal IsraelStewart 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 highenergy nucleusnucleus
collisions obtained through the POWLANG transport setup, developed in the past
to study the propagation of heavy quarks in the QuarkGluon 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 heavyflavour hadrons (or their decay
electrons) and the charged particles produced in the fragmentation of the
heavyquark strings.

We present results for the quenching, elliptic flow and azimuthal
correlations of heavy flavour particles in highenergy nucleusnucleus
collisions obtained through the POWLANG transport setup, developed in the past
to study the propagation of heavy quarks in the QuarkGluon 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
netbaryon density which embed recent lattice QCD results at high temperatures
combined with a hadron resonance gas model in the lowtemperature, confined
phase. In the latter, we allow an implementation of partial chemical
equilibrium, in which particle ratios are fixed at the chemical freezeout, so
that a description closer to the experimental situation is possible. Given the
present uncertainty in the determination of the chemical freezeout temperature
from firstprinciple 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 heavyion 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 heavyflavour observables in heavyion (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 heavyion
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 weakcoupling calculations or from latticeQCD simulations. Our
findings are compared to a rich set of experimental data (Dmesons,
nonphotonic electrons, nonprompt J/\psi's) which have meanwhile become
accessible.

We present ECHOQGP, a numerical code for $(3+1)$dimensional relativistic
viscous hydrodynamics designed for the modeling of the spacetime evolution of
the matter created in high energy nuclear collisions. The code has been built
on top of the \emph{Eulerian Conservative HighOrder} astrophysical code for
general relativistic magnetohydrodynamics [\emph{Del Zanna et al., Astron.
Astrophys. 473, 11, 2007}] and here it has been upgraded to handle the physics
of the QuarkGluon Plasma. ECHOQGP features secondorder treatment of causal
relativistic viscosity effects in both Minkowskian or Bjorken coordinates;
partial or complete chemical equilibrium of hadronic species before kinetic
freezeout; initial conditions based on the optical Glauber model, including a
MonteCarlo routine for eventbyevent fluctuating initial conditions; a
freezeout procedure based on the CooperFrye prescription. The code is
extensively validated against several test problems and results always appear
accurate, as guaranteed by the combination of the conservative
(shockcapturing) approach and the highorder methods employed. ECHOQGP can be
extended to include evolution of the electromagnetic fields coupled to the
plasma.

The quenching of jets (and highpT particle spectra) observed in heavyion
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 colorflow 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 energyloss. 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 jetfragmentation.

A refined version of a multistep calculation of heavyflavor observables in
pp and AA collisions has been developed, based on pQCD at NLO accuracy followed
by parton shower evolution to describe heavyquark 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 PbPb collisions at $\sqrt{s_{NN}}$=2.76 TeV: nuclear
modification factor $R_{AA}$ of the $p_{T}$ spectra at midrapidity of
heavyflavor decay electrons and of exclusively reconstructed opencharm mesons
at different centralities, as well as their ellipticflow $v_{2}$($p_{T}$) in
semicentral 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 mediummodifications of jet observables in
ultrarelativistic heavy ion collisions at RHIC and at the LHC. Here, we note
that such jetmedium interactions alter generically and characteristically the
color correlations in the parton shower. We characterize these effects in a
colordifferential calculation of the mediuminduced gluon radiation spectrum
to first and second order in opacity. By interfacing simple branching histories
of mediummodified 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 jetmedium
interactions give rise to the mediuminduced 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 mediuminduced 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 vacuumlike fragmentation pattern of lower total energy.

Mediuminduced 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
ultrarelativistic heavy ion collisions. Here, we present the first
colordifferential 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 colorsinglet fragments of
mediummodified 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 highpT singlehadron spectra in heavyion collisions is
usually interpreted as due to parton energyloss of highmomentum quarks and
gluons propagating in the plasma. Here, we discuss to what extent this partonic
picture must be complemented by a picture of mediummodified hadronization. In
particular, we show how colorexchange with the medium modifies the properties
of color singletclusters arising from the parton branchings, producing a
softening of the hadron spectra.

We present a study of the heavyflavor dynamics in nucleusnucleus
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 opencharm (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 heavyflavor hadrons and decayelectrons 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 nucleusnucleus collisions at RHIC and LHC, providing results for the
nuclear modification factor R_AA and the elliptic flow coefficient v_2 of the
singleelectron spectra arising from their semileptonic 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 firstprinciple 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 singleelectron spectra to be compared with the
current experimental data measured at RHIC. We focus on the nuclear
modification factor R_AA and on the ellipticflow coefficient v_2, getting, for
sufficiently large p_T, a reasonable agreement.

The stochastic dynamics of heavy quarks in the fireball produced in heavyion
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 singleelectrons with respect to pp collisions is studied. The
transport coefficients are evaluated treating separately the contribution of
soft and hard collisions. The initial heavyquark spectra are generated
according to NLOpQCD, accounting for nuclear effects through recent nPDFs. The
evolution of the medium is obtained from the output of two hydrocodes (ideal
and viscous). The heavyquark fragmentation into hadrons and their final
semileptonic decays are implemented according to up to date experimental data.
A comparison with RHIC data for nonphotonic electron spectra is given.

We propose a model for the propagation of a heavyquark in a hot plasma, to
be viewed as a first step towards a full description of the dynamics of heavy
quark systems in a quarkgluon 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 spacetime 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
semianalytically at the oneloop 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
quarkantiquark pairs are computed through an extension of the results obtained
for a hot QED plasma by describing the heavyquark 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 heavyquark momentum
spectra. Our numerical findings are also compared with the ones obtained in a
stronglycoupled scenario, namely with the transport coefficients predicted
(though with some limitations and ambiguities) by the AdS/CFT correspondence.

The heavyquark spectral function in a hot plasma is reconstructed from the
corresponding euclidean propagator. The latter is evaluated through a
pathintegral simulation. A weakcoupling 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
realtime 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 realtime propagator is the free energy calculated in the
imaginarytime formalism.