Temperature dependence of η/s of strongly interacting matter: Effects of the equation of state and the parametric form of (η/s)(T)

We investigate the temperature dependence of the shear viscosity to entropy density ratio $\ensuremath{\eta}/s$ using a piecewise linear parametrization. To determine the optimal values of the parameters and the associated uncertainties, we perform a global Bayesian model-to-data comparison on $\mathrm{Au}+\mathrm{Au}$ collisions at $\sqrt{{s}_{\mathrm{NN}}}=200$ GeV and $\mathrm{Pb}+\mathrm{Pb}$ collisions at 2.76 TeV and 5.02 TeV, using a $2+1\mathrm{D}$ hydrodynamical model with the Eskola-Kajantie-Ruuskanen-Tuominen (EKRT) initial state. We provide three new parametrizations of the equation of state (EoS) based on contemporary lattice results and hadron resonance gas, and use them and t…

Massive Quarks at One Loop in the Dipole Picture of Deep Inelastic Scattering

We calculate the light cone wave functions for a virtual photon to split into quark-antiquark states, including for the first time quark masses at one loop accuracy. These wave functions can be used to calculate cross sections for several precision probes of perturbative gluon saturation at the Electron-Ion Collider. Using these wave functions we derive, for the first time, the dipole picture DIS cross sections at one loop for longitudinal and transverse virtual photons including quark masses. The quark masses are renormalized in the pole mass scheme, satisfying constraints from the requirement of Lorentz invariance of the quark Dirac and Pauli form factors.

Pinning down QCD-matter shear viscosity in A + A collisions via EbyE fluctuations using pQCD + saturation + hydrodynamics

We compute the initial energy densities produced in ultrarelativistic heavy-ion collisions from NLO perturbative QCD using a saturation conjecture to control soft particle production, and describe the subsequent space-time evolution of the system with hydrodynamics, event by event. The resulting centrality dependence of the low-$p_T$ observables from this pQCD + saturation + hydro ("EKRT") framework are then compared simultaneously to the LHC and RHIC measurements. With such an analysis we can test the initial state calculation, and constrain the temperature dependence of the shear viscosity-to-entropy ratio $\eta/s$ of QCD matter. Using these constraints from the current RHIC and LHC measu…

Event-by-event fluctuations in a perturbative QCD + saturation + hydrodynamics model: Determining QCD matter shear viscosity in ultrarelativistic heavy-ion collisions

We introduce an event-by-event perturbative-QCD + saturation + hydro (“EKRT”) framework for ultrarelativistic heavy-ion collisions, where we compute the produced fluctuating QCD-matter energy densities from next-to-leading-order perturbative QCD using a saturation conjecture to control soft-particle production and describe the space-time evolution of the QCD matter with dissipative fluid dynamics, event by event. We perform a simultaneous comparison of the centrality dependence of hadronic multiplicities, transverse momentum spectra, and flow coefficients of the azimuth-angle asymmetries against the LHC and RHIC measurements. We compare also the computed event-by-event probability distribut…

Latest predictions from the EbyE NLO EKRT model

We present the latest results from the NLO pQCD + saturation + viscous hydrodynamics (EbyE NLO EKRT) model. The parameters in the EKRT saturation model are fixed by the charged hadron multiplicity in the 0-5 \% 2.76 TeV Pb+Pb collisions. The $\sqrt{s}$, $A$ and centrality dependence of the initial particle production follows then from the QCD dynamics of the model. This allows us to predict the $\sqrt{s}$ and $A$ dependence of the particle production. We show that our results are in an excellent agreement with the low-$p_T$ data from 2.76 TeV and 5.02 TeV Pb+Pb collisions at the LHC as well as with the data from the 200 GeV Au+Au collisions at RHIC. In particular, we study the centrality de…

Color charge correlations in the proton at NLO: Beyond geometry based intuition

Color charge correlators provide fundamental information about the proton structure. In this Letter, we evaluate numerically two-point color charge correlations in a proton on the light cone including the next-to-leading order corrections due to emission or exchange of a perturbative gluon. The non-perturbative valence quark structure of the proton is modelled in a way consistent with high-$x$ proton structure data. Our results show that the correlator exhibits startlingly non-trivial behavior at large momentum transfer or central impact parameters, and that the color charge correlation depends not only on the impact parameter but also on the relative transverse momentum of the two gluon pr…

Predictions for multiplicities and flow harmonics in 5.44 TeV Xe+Xe collisions at the CERN Large Hadron Collider

We present the next-to-leading-order event-by-event EKRT model predictions for the centrality dependence of the charged hadron multiplicity in the pseudorapidity interval $|\eta|\le 0.5$, and for the centrality dependence of the charged hadron flow harmonics $v_n\{2\}$ obtained from 2-particle cumulants, in $\sqrt{s_{NN}}=5.44$ TeV Xe+Xe collisions at the CERN Large Hadron Collider. Our prediction for the 0-5 \% central charged multiplicity is $dN_{\rm ch}/d\eta =1218\pm 46$. We also predict $v_n\{2\}$ in Xe+Xe collisions to increase more slowly from central towards peripheral collisions than those in a Pb+Pb system. We find that at $10 \dots 50$\% centralities $v_2\{2\}$ is smaller and $v_…

Extracting qˆ in event-by-event hydrodynamics and the centrality/energy puzzle

In our analysis, we combine event-by-event hydrodynamics, within the EKRT formulation, with jet quenching -ASW Quenching Weights- to obtain high- for charged particles at RHIC and LHC energies for different centralities. By defining a K-factor that quantifies the departure of from an ideal estimate, , we fit the single-inclusive experimental data for charged particles. This K-factor is larger at RHIC than at the LHC but, surprisingly, it is almost independent of the centrality of the collision. peerReviewed

Temperature dependence of η / s of strongly interacting matter: Effects of the equation of state and the parametric form of ( η / s ) ( T )

We investigate the temperature dependence of the shear viscosity to entropy density ratio η/s using a piecewise linear parametrization. To determine the optimal values of the parameters and the associated uncertainties, we perform a global Bayesian model-to-data comparison on Au+Au collisions at √sNN=200 GeV and Pb+Pb collisions at 2.76 TeV and 5.02 TeV, using a 2+1D hydrodynamical model with the Eskola-Kajantie-Ruuskanen-Tuominen (EKRT) initial state. We provide three new parametrizations of the equation of state (EoS) based on contemporary lattice results and hadron resonance gas, and use them and the widely used s95p parametrization to explore the uncertainty in the analysis due to the c…

Impact parameter dependence of color charge correlations in the proton

The impact parameter dependence of color charge correlators in the proton is obtained from the light front formalism in light cone gauge. We include NLO corrections due to the $|qqqg\rangle$ Fock state via light-cone perturbation theory. Near the center of the proton, the $b$-dependence of the correlations is very different from a "transverse profile function". The resulting $t$-dependence of exclusive $J/\Psi$ photoproduction transitions from exponential to power law at $|t| \approx 1$ GeV$^2$. This prediction could be tested at upcoming DIS facilities or in nucleus-proton ultraperipheral collisions (UPCs).

Event-by-event fluctuations in a perturbative QCD plus saturation plus hydrodynamics model : Determining QCD matter shear viscosity in ultrarelativistic heavy-ion collisions

We introduce an event-by-event perturbative-QCD + saturation + hydro ("EKRT") framework for ultrarelativistic heavy-ion collisions, where we compute the produced fluctuating QCD-matter energy densities from next-to-leading-order perturbative QCD using a saturation conjecture to control soft-particle production and describe the space-time evolution of the QCD matter with dissipative fluid dynamics, event by event. We perform a simultaneous comparison of the centrality dependence of hadronic multiplicities, transverse momentum spectra, and flow coefficients of the azimuth-angle asymmetries against the LHC and RHIC measurements. We compare also the computed event-by-event probability distribut…

Stronger C -odd color charge correlations in the proton at higher energy

The non-forward eikonal scattering matrix for dipole-proton scattering at high energy obtains an imaginary part due to a $C$-odd three gluon exchange. We present numerical estimates for the perturbative Odderon amplitude as a function of dipole size, impact parameter, their relative azimuthal angle, and light-cone momentum cutoff $x$. The proton is approximated as $\psi_\mathrm{qqq}|qqq\rangle + \psi_\mathrm{qqqg}|qqqg\rangle$, where $\psi_\mathrm{qqq}$ is a non-perturbative three quark model wave function while the gluon emission is computed in light-cone perturbation theory. We find that the Odderon amplitude increases as $x$ decreases from 0.1 to 0.01. At yet lower $x$, the reversal of t…

Azimuthal correlations in diffractive scattering at the Electron-Ion Collider

We calculate azimuthal correlations between the exclusively produced vector meson and the scattered electron in Deep Inelastic Scattering processes at the future Electron-Ion Collider (EIC). We identify "kinematical" and "intrinsic" contributions to these correlations, and show that the correlations are sensitive to the non-trivial correlations in the gluon distribution of the target. Realistic predictions at the EIC kinematics are provided using two different approaches to describe the dipole-proton interaction at relatively small xx. peerReviewed

Jet quenching as a probe of the initial stages in heavy-ion collisions

Jet quenching provides a very flexible variety of observables which are sensitive to different energy- and time-scales of the strongly interacting matter created in heavy-ion collisions. Exploiting this versatility would make jet quenching an excellent chronometer of the yoctosecond structure of the evolution process. Here we show, for the first time, that a combination of jet quenching observables is sensitive to the initial stages of heavy-ion collisions, when the approach to local thermal equilibrium is expected to happen. Specifically, we find that in order to reproduce at the same time the inclusive particle production suppression, $R_{AA}$, and the high-$p_T$ azimuthal asymmetries, $v…

One-loop corrections to light cone wave functions: the dipole picture DIS cross section

We develop methods needed to perform loop calculations in light cone perturbation theory using a helicity basis, refining the method introduced in our earlier work. In particular this includes implementing a consistent way to contract the four-dimensional tensor structures from the helicity vectors with d-dimensional tensors arising from loop integrals, in a way that can be fully automatized. We demonstrate this explicitly by calculating the one-loop correction to the virtual photon to quark-antiquark dipole light cone wave function. This allows us to calculate the deep inelastic scattering cross section in the dipole formalism to next-to-leading order accuracy. Our results, obtained using …

Pinning down QCD-matter shear viscosity in ultrarelativistic heavy-ion collisions via EbyE fluctuations using pQCD + saturation + hydrodynamics

We introduce an event-by-event pQCD + saturation + hydro ("EKRT") framework for high-energy heavy-ion collisions, where we compute the produced fluctuating QCD-matter energy densities from next-to-leading order (NLO) perturbative QCD (pQCD) using saturation to control soft particle production, and describe the space-time evolution of the QCD matter with viscous hydrodynamics, event by event (EbyE). We compare the computed centrality dependence of hadronic multiplicities, p_T spectra and flow coefficients v_n against LHC and RHIC data. We compare also the computed EbyE probability distributions of relative fluctuations of v_n, as well as correlations of 2 and 3 event-plane angles, with LHC d…

Massive quarks in NLO dipole factorization for DIS : Longitudinal photon

In this work, we will present the first complete calculation of the one-loop longitudinal photon-to-quark-antiquark light cone wave function, with massive quarks. The quark masses are renormalized in the pole mass scheme. The result is used to calculate the next-to-leading order correction to the high energy Deep Inelastic Scattering longitudinal structure function on a dense target in the dipole factorization framework. For massless quarks the next-to-leading order correction was already known to be sizeable, and our result makes it possible to evaluate it also for massive quarks.

Cubic color charge correlator in a proton made of three quarks and a gluon

The three point correlation function of color charge densities is evaluated explicitly in light cone gauge for a proton on the light cone. This includes both $C$-conjugation even and odd contributions. We account for perturbative corrections to the three-quark light cone wave function due to the emission of an internal gluon which is not required to be soft. We verify the Ward identity as well as the cancellation of UV divergences in the sum of all diagrams so that the correlator is independent of the renormalization scale. It does, however, exhibit the well known soft and collinear singularities. The expressions derived here provide the $C$-odd contribution to the initial conditions for hi…

Fluid dynamics with saturated minijet initial conditions in ultrarelativistic heavy-ion collisions

Using next-to-leading order perturbative QCD and a conjecture of saturation to suppress the production of low-energy partons, we calculate the initial energy densities and formation times for the dissipative fluid dynamical evolution of the quark-gluon plasma produced in ultrarelativistic heavy-ion collisions. We identify the framework uncertainties and demonstrate the predictive power of the approach by a good global agreement with the measured centrality dependence of charged particle multiplicities, transverse momentum spectra and elliptic flow simultaneously for the Pb+Pb collisions at the LHC and Au+Au at RHIC. In particular, the shear viscosity in the different phases of QCD matter is…

The one loop gluon emission light cone wave function

Light cone perturbation theory has become an essential tool to calculate cross sections for various small-$x$ dilute-dense processes such as deep inelastic scattering and forward proton-proton and proton-nucleus collisions. Here we set out to do one loop calculations in an explicit helicity basis in the four dimensional helicity scheme. As a first process we calculate light cone wave function for one gluon emission to one-loop order in Hamiltonian perturbation theory on the light front. We regulate ultraviolet divergences with transverse dimensional regularization and soft divergences with using a cut-off on longitudinal momentum. We show that when all the renormalization constants are comb…

Massive quarks in NLO dipole factorization for DIS : Transverse photon

We calculate the light cone wave functions for the QCD Fock components in a transverse virtual photon necessary for applications at next-to-leading order in the QCD coupling, including quark masses. We present a detailed calculation of both the one loop wave function for the quark-antiquark Fock component and the tree level wave function for the quark-antiquark-gluon Fock component. The quark masses are renormalized in the pole mass scheme, satisfying constraints from the requirement of Lorentz invariance. In particular the quark Pauli form factor at NLO is recovered from the on-shell limit of the quark-antiquark Fock component. We use our result to calculate the next-to-leading order corre…

Next-to-leading order improved perturbative QCD + saturation + hydrodynamics model for A + A collisions

We calculate initial conditions for the hydrodynamical evolution in ultrarelativistic heavy-ion collisions at the LHC and RHIC in an improved next-to-leading order perturbative QCD + saturation framework. Using viscous relativistic hydrodynamics, we show that we obtain a good simultaneous description of the centrality dependence of charged particle multiplicities, transverse momentum spectra and elliptic flow at the LHC and at RHIC. In particular, we discuss how the temperature dependence of the shear viscosity is constrained by these data.

Temperature dependence of η/s : uncertainties from the equation of state

We perform a global model-to-data comparison on Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV and Pb+Pb collisions at $2.76$ TeV and $5.02$ TeV, using a 2+1D hydrodynamics model with the EKRT initial state and a shear viscosity over entropy density ratio $(\eta/s)(T)$ with a linear $T$ dependence. To quantify the amount of uncertainty due to the choice of the equation of state (EoS), we compare analysis results based on four different EoSs: the well known $s95p$ parametrisation, an updated parametrisation based on the same list of particles in hadron resonance gas, but using recent lattice results for the partonic part of the EoS, and two new parametrisations based on the Particle Data Group …

Temperature dependence of eta/s of strongly interacting matter : Effects of the equation of state and the parametric form of (eta/s)(T)

We investigate the temperature dependence of the shear viscosity to entropy density ratio $\eta/s$ using a piecewise linear parametrization. To determine the optimal values of the parameters and the associated uncertainties, we perform a global Bayesian model-to-data comparison on Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV and Pb+Pb collisions at $2.76$ TeV and $5.02$ TeV, using a 2+1D hydrodynamical model with the EKRT initial state. We provide three new parametrizations of the equation of state (EoS) based on contemporary lattice results and hadron resonance gas, and use them and the widely used $s95p$ parametrization to explore the uncertainty in the analysis due to the choice of the eq…

From minijet saturation to global observables in A + A collisions at the LHC and RHIC

We review the recent results from the computation of saturated next-to-leading order perturbative QCD minijet intial conditions combined with viscous hydrodynamical evolution of ultrarelativistic heavy-ion collisions at the LHC and RHIC. Comparison with experimental data is shown.

Proton Structure Functions at Next-to-Leading Order in the Dipole Picture with Massive Quarks

We predict heavy quark production cross sections in deep inelastic scattering at high energy by applying the color glass condensate effective theory. We demonstrate that, when the calculation is performed consistently at next-to-leading order accuracy with massive quarks, it becomes possible, for the first time in the dipole picture with perturbatively calculated center-of-mass energy evolution, to simultaneously describe both the light and heavy quark production data at small xBj. Furthermore, we show how the heavy quark cross section data provides additional strong constraints on the extracted nonperturbative initial condition for the small-xBj evolution equations.

Predictions for 5.023 TeV Pb + Pb collisions at the CERN Large Hadron Collider

We compute predictions for various low-transverse-momentum bulk observables in √sNN = 5.023 TeV Pb+Pb collisions at the CERN Large Hadron Collider (LHC) from the event-by-event next-to-leading-order perturbative-QCD + saturation + viscous hydrodynamics (“EKRT”) model. In particular, we consider the centrality dependence of charged hadron multiplicity, flow coefficients of the azimuth-angle asymmetries, and correlations of event-plane angles. The centrality dependencies of the studied observables are predicted to be very similar to those at 2.76 TeV, and the magnitudes of the flow coefficients and event-plane angle correlations are predicted to be close to those at 2.76 TeV. The flow coeffic…

Extractingqˆin event-by-event hydrodynamics and the centrality/energy puzzle

Abstract In our analysis, we combine event-by-event hydrodynamics, within the EKRT formulation, with jet quenching -ASW Quenching Weights- to obtain high- p T R AA for charged particles at RHIC and LHC energies for different centralities. By defining a K -factor that quantifies the departure of q ˆ from an ideal estimate, K = q ˆ / ( 2 ϵ 3 / 4 ) , we fit the single-inclusive experimental data for charged particles. This K -factor is larger at RHIC than at the LHC but, surprisingly, it is almost independent of the centrality of the collision.

Multiplicities andpTspectra in ultrarelativistic heavy ion collisions from a next-to-leading order improved perturbative QCD+saturation+hydrodynamics model

We bring the EKRT framework, which combines perturbative QCD (pQCD) minijet production with gluon saturation and hydrodynamics, to next-to-leading order (NLO) in pQCD as rigorously as possible. We chart the model uncertainties, and study the viability and predictive power of the model in the light of the RHIC and LHC measurements in central $A+A$ collisions. In particular, we introduce a new set of measurement functions to define the infrared- and collinear-safe minijet transverse energy, ${E}_{T}$, in terms of which we formulate the saturation. We update the framework with the EPS09 NLO nuclear parton distributions (nPDFs), and study the propagation of the nPDF uncertainties into the compu…

Constraining energy loss from high-pT azimuthal asymmetries

The nuclear modification factor $R_{\rm AA}$ has been satisfactorily described by various jet quenching models. Nonetheless, all these formalisms, until very recently, underpredicted the high-$p_{\rm T}$ (> 10 GeV) elliptic flow $v_2$. We find that the simultaneous description of these observables requires to strongly suppress the quenching for the first $\sim 0.6$ fm after the collision. This shows the potential of jet quenching observables to constrain the dynamics of the initial stages of the evolution.

Extracting $\hat{q}$ in event-by-event hydrodynamics and the centrality/energy puzzle

In our analysis, we combine event-by-event hydrodynamics, within the EKRT formulation, with jet quenching -ASW Quenching Weights- to obtain high-$p_T$ $R_{\rm AA}$ for charged particles at RHIC and LHC energies for different centralities. By defining a $K$-factor that quantifies the departure of $\hat{q}$ from an ideal estimate, $K = \hat{q}/(2\epsilon^{3/4})$, we fit the single-inclusive experimental data for charged particles. This $K$-factor is larger at RHIC than at the LHC but, surprisingly, it is almost independent of the centrality of the collision.

Latest results from the EbyE NLO EKRT model

We review the results from the event-by-event next-to-leading order perturbative QCD + saturation + viscous hydrodynamics (EbyE NLO EKRT) model. With a simultaneous analysis of LHC and RHIC bulk observables we systematically constrain the QCD matter shear viscosity-to-entropy ratio eta/s(T), and test the initial state computation. In particular, we study the centrality dependences of hadronic multiplicities, pT spectra, flow coefficients, relative elliptic flow fluctuations, and various flow-correlations in 2.76 and 5.02 TeV Pb+Pb collisions at the LHC and 200 GeV Au+Au collisions at RHIC. Overall, our results match remarkably well with the LHC and RHIC measurements, and predictions for the…

Minijet initial state of heavy-ion collisions from next-to-leading order perturbative QCD

Systematics of parton-medium interaction from RHIC to LHC

Despite a wealth of experimental data for high-P_T processes in heavy-ion collisions, discriminating between different models of hard parton-medium interactions has been difficult. A key reason is that the pQCD parton spectrum at RHIC is falling so steeply that distinguishing even a moderate shift in parton energy from complete parton absorption is essentially impossible. In essence, energy loss models are effectively only probed in the vicinity of zero energy loss and, as a result, at RHIC energies only the pathlength dependence of energy loss offers some discriminating power. At LHC however, this is no longer the case: Due to the much flatter shape of the parton p_T spectra originating fr…