Source author record

Kari J. Eskola

Kari J. Eskola appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

42works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

42 published item(s)

preprint2022arXiv

EPPS21: A global QCD analysis of nuclear PDFs

We present an updated global analysis of collinearly factorized nuclear parton distribution functions (PDFs) at next-to-leading order in perturbative QCD. In comparison to our previous fit, EPPS16, the present analysis includes more data from proton-lead collisions at the Large Hadron Collider: 5TeV double-differential CMS dijet and LHCb D-meson data, as well as 8TeV CMS W$^\pm$ data. These new data lead to significantly better-constrained gluon distributions at small and intermediate values of the momentum fraction $x$, confirming the presence of shadowing and antishadowing for gluons in large nuclei. In addition, we include Jefferson Lab measurements of deeply inelastic scattering which probe nuclear PDFs at large $x$ and low virtualities. For the first time within the Hessian framework, we now also explore the uncertainties of nuclear PDFs due to the errors in the baseline proton PDFs. We release the results of our analysis as a new public parametrization of nuclear PDFs called EPPS21.

preprint2022arXiv

Next-to-leading order QCD predictions for exclusive $J/ψ$ photoproduction in ultraperipheral Pb+Pb collisions at the LHC

We present the first study of rapidity-differential cross sections of coherent exclusive $J/ψ$ photoproduction in heavy-ion ultraperipheral collisions (UPCs) at the LHC, $\text{d}σ(\text{Pb}+\text{Pb}\rightarrow \text{Pb}+J/ψ+\text{Pb})/\text{d}y$, to NLO in collinear factorisation. Taking the generalised parton distributions as their forward limit counterparts, the NLO contributions in the cross sections are quantified and we show that the real part of the amplitude and quark-PDF contributions must not be neglected. We further evaluate the uncertainties arising from the scale-choice at both LO and NLO, and compare our results with exclusive $J/ψ$ photoproduction data in Pb+Pb UPCs from LHC. The scale dependence in $\text{d}σ(\text{Pb}+\text{Pb}\rightarrow \text{Pb}+J/ψ+\text{Pb})/\text{d}y$ is sizable, but we can find a scale-choice at NLO that reproduces the Pb+Pb UPC data at both Run 1 and Run 2 collision energies. This process has traditionally been suggested to be a direct probe of nuclear gluon distributions. We show that the situation changes rather dramatically from LO to NLO, where the interplay of both gluon and quark contributions, and in particular their relative signs, play a significant role. We end with some concluding remarks about the outlook of this observable within a tamed collinear factorisation.

preprint2022arXiv

Proton-PDF uncertainties in extracting nuclear PDFs from $W^\pm$ production in p+Pb collisions

We discuss the recent CMS Collaboration measurement of $W^\pm$ boson production in p+Pb collisions at 8.16 TeV in terms of the constraining power on nuclear parton distribution functions (PDFs). The impact of the free-proton PDF uncertainties on the nuclear PDF extraction is quantified by using a theoretical covariance-matrix method and Hessian PDF reweighting. We discuss different ways to mitigate these theoretical uncertainties, including self-normalization, forward-to-backward ratios and nuclear-modification ratios. It is found that none of these methods offer perfect cancellation of the free-proton PDFs but, with the present data uncertainties, the residual free-proton-PDF dependence has, conveniently for the global analyses, little effect on the extraction of the nuclear modifications. Based on a simple estimate of obtainable statistics at the LHC Run 3, we argue that this will change in the near future and it becomes more important to propagate the proton-PDF uncertainties accordingly. Using the obtained information on the correlations of the free-proton uncertainties, we also identify a new charge asymmetry ratio, where the cancellation of the proton-PDF uncertainties is found to be extremely good.

preprint2020arXiv

A QCD analysis of LHCb D-meson data in p+Pb collisions

We scrutinize the recent LHCb data for D$^0$-meson production in p+Pb collisions within a next-to-leading order QCD framework. Our calculations are performed in the SACOT-$m_{\rm T}$ variant of the general-mass variable-flavour-number scheme (GM-VFNS), which has previously been shown to provide a realistic description of the LHC p+p data. Using the EPPS16 and nCTEQ15 nuclear parton distribution functions (PDFs) we show that a very good agreement is obtained also in the p+Pb case both for cross sections and nuclear modification ratios in the wide rapidity range covered by the LHCb data. Encouraged by the good correspondence, we quantify the impact of these data on the nuclear PDFs by the Hessian reweighting technique. We find compelling direct evidence of gluon shadowing at small momentum fractions $x$, with no signs of parton dynamics beyond the collinear factorization. We also compare our theoretical framework to a fixed-order calculation supplemented with a parton shower. While the two frameworks differ in the absolute cross sections, these differences largely cancel in the nuclear modification ratios. Thus, the constraints for nuclear PDFs appear solid.

preprint2020arXiv

Evidence of shadowing in inelastic nucleon-nucleon cross section

The Glauber modeling plays a key role in centrality-dependent measurements of heavy-ion collisions. A central input parameter in Glauber models is the inelastic nucleon-nucleon cross section $σ_{\text{nn}}^{\text{inel}}$ which is nearly always taken from proton-proton measurements. At the LHC energies $σ_{\text{nn}}^{\text{inel}}$ depends on the QCD dynamics at small $x$ and low interaction scales where the shadowing/saturation phenomena are expected to become relatively more important for larger nuclei than for the proton. Thus, $σ_{\text{nn}}^{\text{inel}}$ e.g. in Pb+Pb collisions may well be lower than what is seen in proton-proton collisions. In this talk, we demonstrate how to use the recent $W^\pm$ and $Z$ measurements as a "standard candle" to extract $σ_{\text{nn}}^{\text{inel}}$ in Pb+Pb collisions. Our analysis -- built on the ATLAS data, state-of-the-art NNLO QCD calculations and nuclear PDFs -- indicate that at the LHC energies $σ_{\text{nn}}^{\text{inel}}$ in Pb+Pb collisions is suppressed relative to the proton-proton measurements by tens of percents. We demonstrate that this is in line with expectations from nuclear PDFs.

preprint2020arXiv

Impact of dijet and D-meson data from 5.02 TeV p+Pb collisions on nuclear PDFs

We discuss the new constraints on gluon parton distribution function (PDF) in lead nucleus, derivable with the Hessian PDF reweighting method from the 5.02 TeV p+Pb measurements of dijet (CMS) and $D^0$-meson (LHCb) nuclear modification ratios. The impact is found to be significant, placing stringent constraints in the mid- and previously unconstrained small-$x$ regions. The CMS dijet data confirm the existence of gluon anti-shadowing and the onset of small-$x$ shadowing, as well as reduce the gluon PDF uncertainties in the larger-$x$ region. The gluon constraints from the LHCb $D^0$ data, reaching down to $x \sim 10^{-5}$ and derived in a NLO perturbative QCD approach, provide a remarkable reduction of the small-$x$ uncertainties with a strong direct evidence of gluon shadowing. Furthermore, we find a good description of the data even down to zero $D^0$-meson transverse momentum within a purely DGLAP-based approach without a need for imposing any non-linear effects. Importantly, the constraints obtained from the dijet and $D^0$ data are mutually fully consistent, supporting the universality of nuclear PDFs in hard-scattering processes.

preprint2016arXiv

Neutron-skin effect and centrality dependence of high-$p_{\mathrm{T}}$ observables in nuclear collisions

We report on our studies of the neutron-skin effects in high-$p_{\mathrm{T}}$ observables at the LHC. We study the impact of the neutron-skin effect on the centrality dependence of inclusive direct photon, high-$p_{\mathrm{T}}$ hadron and $W^{\pm}$ production in nuclear collisions at the LHC. The neutron-skin effect refers to the observation that in spherical heavy nuclei, the tail of the neutron distribution extends farther than the distribution of protons, which can affect observables sensitive to electroweak phenomena in very peripheral collisions. We quantify this effect for direct photons, charged hadrons and W bosons as a function of the collision centrality. In the case of direct photons we find that it will be difficult to resolve the neutron-skin effect, given the uncertainties in the nuclear PDFs and their spatial dependence. With charged hadrons and W's, however, up to 20~\% unambiguous effects are expected for most peripheral collisions.

preprint2015arXiv

Nuclear PDF constraints from p+Pb collisions at the LHC

As the current nuclear PDF analyses are mainly constrained by fixed-target Drell-Yan and deeply inelastic scattering data only the quark nuclear modifications at fairly large $x$ values are in a good control. Inclusive pion production in d+Au collisions at RHIC provides some constraints for gluons but due to the limited kinematic reach of the data the gluon modifications remain uncertain especially at small values of $x$. In this talk, we discuss how the existing data from p+Pb collisions at the LHC can improve the nuclear PDF fits and which measurements would be sensitive to the small-$x$ gluons. In particular, we consider inclusive hadron production, compare this to direct photons, and show estimates of the effect of CMS dijet measurements to the uncertainty of nuclear gluon distributions.

preprint2014arXiv

Centrality and rapidity dependence of inclusive pion and prompt photon production in p+Pb collisions at the LHC with EPS09s nPDFs

The centrality dependencies of the inclusive neutral pion and prompt photon nuclear modification factors for p+Pb collisions at the LHC are studied using a spatially dependent set of nuclear PDFs, EPS09s. The calculations are performed at mid- and forward rapidities searching for an observable which would optimally probe the spatial dependence of the nuclear PDFs. In addition, we discuss to which $x$ values of the nucleus the different observables are sensitive.

preprint2014arXiv

Confronting current NLO parton fragmentation functions with inclusive charged-particle spectra at hadron colliders

The inclusive spectra of charged particles measured at high transverse momenta ($p_T\gtrsim$2GeV/c) in proton-proton and proton-antiproton collisions in the range of center-of-mass energies $\sqrt{s}=200-7000$GeV are compared with next-to-leading order perturbative QCD calculations using seven recent sets of parton-to-hadron fragmentation functions (FFs). Accounting for the uncertainties in the scale choices and in the parton distribution functions, we find that most of the theoretical predictions tend to overpredict the measured LHC and Tevatron cross sections by up to a factor of two. We identify the currently too-hard gluon-to-hadron FFs as the probable source of the problem, and justify the need to refit the FFs using the available LHC and Tevatron data in a region of transverse momenta, $p_T\gtrsim$10GeV/c, which is supposedly free from additional non-perturbative contributions and where the scale uncertainty is only modest.

preprint2014arXiv

Constraining nPDFs with inclusive pions and direct photons at forward rapidities in p+Pb collisions at the LHC

In this talk, we present NLO pQCD predictions for inclusive pion and direct photon nuclear modifications in p+Pb collisions at mid- and forward rapidities at the LHC. In addition to the minimum bias predictions, we also address the centrality dependence with spatially dependent nuclear PDFs. To understand which regions of the nuclear momentum fraction $x_2$ these observables predominantly probe, we present also the underlying $x_2$ distributions at different rapidities. We are led to conclude that the isolated photons at forward rapidities are more sensitive to the small-$x_2$ dynamics than the inclusive pions.

preprint2014arXiv

Dijets in p+Pb collisions and their quantitative constraints for nuclear PDFs

We present a perturbative QCD analysis concerning the production of high-pT dijets in p+Pb collisions at the LHC. The next-to-leading order corrections, scale variations and free-proton PDF uncertainties are found to have only a relatively small influence on the normalized dijet rapidity distributions. Interestingly, however, these novel observables prove to retain substantial sensitivity to the nuclear effects in the PDFs. Especially, they serve as a more robust probe of the nuclear gluon densities at $x>0.01$, than e.g. the inclusive hadron production. We confront our calculations with the recent data by the CMS collaboration. These preliminary data lend striking support to the gluon antishadowing similar to that in the EPS09 nuclear PDFs.

preprint2014arXiv

LHC data challenges the contemporary parton-to-hadron fragmentation functions

We discuss the inclusive high-pT charged-particle production in proton-proton collisions at the LHC. The experimental data are compared to the NLO perturbative QCD calculations employing various sets of parton-to-hadron fragmentation functions. Most of the theoretical predictions are found to disastrously overpredict the measured cross sections, even if the scale variations and PDF errors are accounted for. The problem appears to arise from the presently too hard gluon-to-hadron fragmentation functions.

preprint2014arXiv

Probing the small-$x$ nuclear gluon distributions with isolated photons at forward rapidities in p+Pb collisions at the LHC

Inclusive direct photon production in p+Pb collisions at the LHC is studied within the NLO perturbative QCD. Our aim is to quantify the dominant $x$ regions probed at different rapidities and to identify the best conditions for testing the nuclear gluon parton distribution functions (nPDFs) at small $x$. A comparison to the inclusive pion production reveals that from these two processes the photons carry more sensitivity to the small-$x$ partons and that this sensitivity can be further increased by imposing an isolation cut for the photon events. The details of the isolation criteria, however, seem to make only a small difference to the studied $x$ sensitivity and have practically no effect on the expected nuclear modifications. We consider also the yield asymmetry between forward and backward rapidities which can be used to probe the nPDFs irrespectively of whether an accurate p+p baseline is available.

preprint2014arXiv

Zeroing in on the initial state -- tomography using bulk, jets and photons

One of the unsolved problems in the current 'standard model' of heavy ion physics is the apparent rapid thermalization of QCD matter in the pre-equilibrium stage. While it is challenging to probe this mechanism directly, there are now several observables available which allow tomographic imaging of the initial state geometry, which is expected to carry remnant information of the equilibration mechanism. On the fluid dynamics side, scaled fluctuations in the momentum space anisotropy parameters v_n image the initial eccentricity fluctuations epsilon_n almost directly with only a weak dependence on the details of the fluid dynamical evolution. From a different direction, due to the strong non-linear dependence of their emission rates on temperature, thermal photons and their v_n are very sensitive to the initial state graininess. Finally, the v_2 and v_3 of high P_T hadrons coming from hard processes reflect the attenuation pattern of partons propagating through the inhomogeneous matter density after some fluid dynamical evolution. Combining information from all these channels does not yet lead to a fully consistent picture, however intriguing trends pointing towards non-trivial initial state dynamics emerge.

preprint2013arXiv

A perturbative QCD study of dijets in p+Pb collisions at the LHC

Inspired by the recent measurements of the CMS collaboration, we report a QCD study of dijet production in proton+lead collisions at the LHC involving large-transverse-momentum jets, $p_T \gtrsim 100$ GeV. Examining the inherent uncertainties of the next-to-leading order perturbative QCD calculations and their sensitivity to the free proton parton distributions (PDFs), we observe a rather small, typically much less than 5% clearance for the shape of the dijet rapidity distribution within approximately 1.5 units around the midrapidity. Even a more stable observable is the ratio between the yields in the positive and negative dijet rapidity, for which the baseline uncertainty can be made negligible by imposing a symmetric jet rapidity acceptance. Both observables prove sensitive to the nuclear modifications of the gluon distributions, the corresponding uncertainties clearly exceeding the estimated baseline uncertainties from the free-proton PDFs and scale dependence. From a theoretical point of view, these observables are therefore very suitable for testing the validity of the collinear factorization and have a high potential to provide precision constraints for the nuclear PDFs.

preprint2013arXiv

Centrality dependence of inclusive prompt photon production in d+Au, Au+Au, p+Pb, and Pb+Pb collisions

We calculate the centrality dependence of the midrapidity nuclear modification for inclusive prompt photon production in d+Au and Au+Au collisions at RHIC and in p+Pb and Pb+Pb collisions at the LHC. Our results, using the recent spatially dependent nuclear PDF set EPS09s, are consistent with the existing high-pT data from the PHENIX and CMS collaborations. The good agreement even in the case of nucleus+nucleus collisions suggests that the high-pT direct photon production is not significantly altered by the strongly interacting medium produced in such collisions. We find the centrality dependence of the nuclear modifications generally rather weak but perhaps measurable at low pT.

preprint2013arXiv

Elliptic flow of thermal photons from event-by-event hydrodynamic model

Elliptic flow of direct photons in relativistic heavy ion collisions is believed to be dominated by contribution from thermal radiation of quark gluon plasma up to $p_T$ $\sim 5$ GeV/$c$, although other sources start outshining the thermal contribution at already smaller values of $p_T$ in the direct photon spectrum. The elliptic flow of thermal photons from ideal hydrodynamics considering a smooth initial density distribution under-predicts the PHENIX direct photon data from 200A GeV Au+Au collisions at RHIC by a large margin in the range $1 < p_T < 5$ GeV/$c$. However, a significant enhancement of thermal photon production due to fluctuations in the initial QCD matter density distributions is expected. We show that such fluctuations result in substantially larger photon elliptic flow for $p_T > 2.5 $ GeV/$c$ compared to a smooth initial-state-averaged density profile. The results from event-by-event hydrodynamics are found to be sensitive to the fluctuation size parameter. However, the effects of initial state fluctuations are insufficient to account for the discrepancy to the PHENIX data for direct photon elliptic flow. Furthermore, the photon $v_2$ is reduced even more when we include the NLO pQCD prompt photon component. We also calculate the spectra and elliptic flow of thermal photons for 2.76A TeV Pb+Pb collisions at LHC and for the 0--40% centrality bin. Thermal photons from event-by-event hydrodynamics along with prompt photons from NLO pQCD calculations explain the ALICE preliminary direct photon data well in the region $p_T \ge 2.5$ GeV/c. Similar to RHIC, the elliptic flow results at LHC are again found to be much smaller than the ALICE preliminary $v_2$ data.

preprint2012arXiv

Centrality and initial formation time dependence of the emission of thermal photons from fluctuating initial conditions at RHIC and LHC

Event-by-event fluctuating initial conditions (IC) in the ideal hydrodynamic calculation are known to enhance the production of thermal photons significantly compared to a smooth initial state averaged profile in the range $p_T >$ 1 GeV/$c$ for 200A GeV Au+Au collisions at RHIC and 2.76A TeV Pb+Pb collisions at LHC. The 'hotspots' or the over-dense regions in the fluctuating IC produce more high $p_T$ photons compared to the smooth IC due to the strong temperature dependent emission of the thermal radiation. This enhancement is expected to be more pronounced for peripheral collisions, for lower beam energies, and for larger values of plasma formation time. A suitably normalized ratio of central to peripheral yield of thermal photons ($R_{cp}^γ$) is a potential probe to study the density fluctuations and their size in the initial conditions.

preprint2012arXiv

Collision centrality and $τ_0$ dependence of the emission of thermal photons from fluctuating initial state in ideal hydrodynamic calculation

Fluctuations in the initial QCD matter density distribution are found to enhance the production of thermal photons significantly in the range 2 \leq pT \leq 4 GeV/c compared to a smooth initial state averaged profile in ideal hydrodynamic calculation for 200 AGeV Au+Au collisions at the Relativistic Heavy Ion Collider (RHIC) and 2.76 ATeV Pb+Pb collisions at the Large Hadron Collider (LHC). The thermal emission of photons is strongly dependent on the initial temperature of the system where the presence of 'hotspots' in the initial state translates into enhanced production of photons compared to a smooth profile. The effect of fluctuations in the initial state is found to be stronger for peripheral collisions and for lower beam energies. The pT spectra are found to be quite sensitive to the value of the initial formation time of the plasma which is not known unambiguously and which may vary with collision centralities at a particular beam energy. Increase in the value of the formation time lowers the production of thermal photons compared to the results from a shorter formation time. However, the relative enhancement from fluctuating initial tates (compared to a smooth initial state) is found to be stronger for the larger values of formation time. The pT spectra alone are found to be insufficient to quantify the fluctuations in the initial density distribution due to the uncertainties in the initial conditions. A suitably normalized ratio of central-to-peripheral yield as a function of collision centrality and pT can be a useful measure of the fluctuation size scale.

preprint2012arXiv

EPS09s and EKS98s: Impact parameter dependent nPDF sets

In our recent study we have determined two new spatially dependent nuclear PDF (nPDF) sets, EPS09s and EKS98s. With these, the hard-process cross-sections can be calculated in different centrality classes consistently with the globally analyzed nPDFs for the first time. The sets were determined by exploiting the $A$-systematics of the globally fitted nPDF sets, EPS09 and EKS98. For the spatial dependence of the nPDFs we used a power series ansatz in the nuclear thickness function $T_A$. In this flash talk we introduce the framework, and present our NLO EPS09s-based predictions for the nuclear modification factor in four centrality classes for inclusive neutral pion production in p+Pb collisions at the LHC and for inclusive prompt photon production in d+Au collisions at RHIC at midrapidity.

preprint2012arXiv

Extracting the impact parameter dependence of the nPDFs from the EKS98 and EPS09 global fits

As all the globally fitted nuclear PDFs (nPDFs) have been so far impact parameter independent, it has not been possible to calculate the hard process cross sections in different centrality classes consistently with the global analyses. In \cite{Helenius:2012wd} we have offered a solution to this problem by determining two spatially dependent nPDF sets, \texttt{EPS09s} and \texttt{EKS98s}, using the $A$-systematics of the earlier global fits EPS09 and EKS98 and an assumption that the spatial dependence can be written as a power series of the nuclear thickness function. For a data comparison, we have calculated the nuclear modification factor of inclusive neutral pion production in d+Au collisions at RHIC in four centrality bins at midrapidity and compared these to a PHENIX measurement. In addition, we have also performed a similar calculation for inclusive photon production in d+Au collisions at RHIC.

preprint2012arXiv

Global analysis of nuclear PDFs - latest developments

In this review talk I discuss the latest developments in the DGLAP-based global analysis of nuclear parton distribution functions (nPDFs), focusing on the recent EPS09, nCTEQ and DSSZ global fits. I also briefly review the recent analysis for assigning a spatial dependence to the globally analysed nPDFs, resulting in the new sets EPS09s and EKS98s ("s" for spatial). With these, one can now compute nuclear hard-process cross sections and estimate their nPDF-originating uncertainties in different centrality classes for the first time consistently with the global nPDF fits.

preprint2012arXiv

Impact-parameter dependent nuclear parton distribution functions: EPS09s and EKS98s and their applications in nuclear hard processes

We determine the spatial (impact parameter) dependence of nuclear parton distribution functions (nPDFs) using the $A$-dependence of the spatially independent (averaged) global fits EPS09 and EKS98. We work under the assumption that the spatial dependence can be formulated as a power series of the nuclear thickness functions $T_A$. To reproduce the $A$-dependence over the entire $x$ range we need terms up to $[T_A]^4$. As an outcome, we release two sets, EPS09s (LO, NLO, error sets) and EKS98s, of spatially dependent nPDFs for public use. We also discuss the implementation of these into the existing calculations. With our results, the centrality dependence of nuclear hard-process observables can be studied consistently with the globally fitted nPDFs for the first time. As an application, we first calculate the LO nuclear modification factor $R^{1jet}_{AA}$ for primary partonic-jet production in different centrality classes in Au+Au collisions at RHIC and Pb+Pb collisions at LHC. Also the corresponding central-to-peripheral ratios $R_{CP}^{1jet}$ are studied. We also calculate the LO and NLO nuclear modification factors for single inclusive neutral pion production, $R_{dAu}^{π^0}$, at mid- and forward rapidities in different centrality classes in d+Au collisions at RHIC. In particular, we show that our results are compatible with the PHENIX mid-rapidity data within the overall normalization uncertainties given by the experiment. Finally, we show our predictions for the corresponding modifications $R_{pPb}^{π^0}$ in the forthcoming p+Pb collisions at LHC.

preprint2012arXiv

Influence of initial state fluctuations on the production of thermal photons

Inhomogeneities in the initial QCD matter density distribution increase the production of thermal photons significantly compared to a smooth initial-state-averaged profile in the region $p_T > 1$ GeV/$c$ in an ideal hydrodynamic calculation. This relative enhancement is more pronounced for peripheral collisions, for smaller size systems as well as for lower beam energies. A suitably normalized ratio of central-to-peripheral yield of thermal photons reduce the uncertainties in the hydrodynamical initial conditions and can be a useful parameter to study the density fluctuations and their size. The fluctuations in the initial density distribution also lead to a larger elliptic flow of thermal photons for $p_T >$ 2.0 GeV/$c$ compared to the flow from a smooth profile.

preprint2012arXiv

Modeling the Impact Parameter Dependence of the nPDFs With EKS98 and EPS09 Global Fits

So far the nuclear PDFs (nPDFs) in the global DGLAP fits have been taken to be spatially independent. In this work, using the $A$-dependence of the globally fitted sets EPS09 and EKS98, we have determined the spatial dependence of the nPDFs in terms of powers of the nuclear thickness functions. New spatially dependent nPDF sets EPS09s (NLO, LO, error sets) and EKS98s (LO) are released. As an application, we consider the nuclear modification factor $R_{dAu}^{π^0}$ at midrapidity for neutral pion production in deuteron-gold collisions at RHIC in NLO. Comparison with the PHENIX data in different centrality classes is also shown. In addition, predictions for the corresponding nuclear modification factor $R_{pPb}^{π^0}$ in proton-lead collisions at the LHC are discussed.

preprint2011arXiv

Elliptic flow from event-by-event hydrodynamics

We present an event-by-event hydrodynamical framework which takes into account the initial density fluctuations arising from a Monte Carlo Glauber model. The elliptic flow is calculated with the event plane method and a one-to-one comparison with the measured event plane $v_2$ is made. Both the centrality- and $p_T$-dependence of the $v_2$ are remarkably well reproduced. We also find that the participant plane is a quite good approximation for the event plane.

preprint2011arXiv

Elliptic flow of thermal photons in heavy-ion collisions at Relativistic Heavy Ion Collider and Large Hadron Collider

We calculate the thermal photon transverse momentum spectra and elliptic flow in $\sqrt{s_{NN}} = 200$ GeV Au+Au collisions at RHIC and in $\sqrt{s_{NN}} = 2.76$ TeV Pb+Pb collisions at the LHC, using an ideal-hydrodynamical framework which is constrained by the measured hadron spectra at RHIC and LHC. The sensitivity of the results to the QCD-matter equation of state and to the photon emission rates is studied, and the photon $v_2$ is discussed in the light of the photonic $p_T$ spectrum measured by the PHENIX Collaboration. In particular, we make a prediction for the thermal photon $p_T$ spectra and elliptic flow for the current LHC Pb+Pb collisions.

preprint2011arXiv

Energy loss in a fluctuating hydrodynamical background

Recently it has become apparent that event-by-event fluctuations in the initial state of hydrodynamical modelling of ultrarelativistic heavy-ion collisions are crucial in order to understand the full centrality dependence of the elliptic flow coefficient v_2. In particular, in central collisions the density fluctuations play a major role in generating the spatial eccentricity in the initial state. This raises the question to what degree high P_T physics, in particular leading-parton energy loss, which takes place in the background of an evolving medium, is sensitive to the presence of the event-by-event density fluctuations in the background. In this work, we report results for the effects of fluctuations on the nuclear modification factor R_AA in both central and noncentral sqrt(s_NN) = 200 GeV Au+Au collisions at RHIC. Two different types of energy-loss models, a radiative and an elastic, are considered. In particular, we study the dependence of the results on the assumed spatial size of the density fluctuations, and discuss the angular modulation of R_AA with respect to the event plane.

preprint2011arXiv

Enhancement of thermal photon production in event-by-event hydrodynamics

Thermal photon emission is widely believed to reflect properties of the earliest, hottest evolution stage of the medium created in ultra-relativistic heavy-ion collisions. Previous computations of photon emission have been carried out using a hydrodynamical medium description with smooth, averaged initial conditions. Recently, more sophisticated hydrodynamical models which calculate observables by averaging over many evolutions with event-by-event fluctuating initial conditions (IC) have been developed. Given their direct connection to the early time dynamics, thermal photon emission appears an ideal observable to probe fluctuations in the medium initial state. In this work, we demonstrate that including fluctuations in the IC may lead to an enhancement of the thermal photon yield of about a factor of 2 in the region $2 < p_T < 4$ GeV/$c$ (where thermal photon production dominates the direct photon yield) compared to a scenario using smooth, averaged IC. Consequently, a much better agreement with PHENIX data is found. This can be understood in terms of the strong temperature dependence of thermal photon production, translating into a sensitivity to the presence of 'hotspots' in an event and thus establishing thermal photons as a suitable probe to characterize IC fluctuations.

preprint2011arXiv

Event-by-event hydrodynamics and elliptic flow from fluctuating initial state

We develop a framework for event-by-event ideal hydrodynamics to study the differential elliptic flow which is measured at different centralities in Au+Au collisions at Relativistic Heavy Ion Collider (RHIC). Fluctuating initial energy density profiles, which here are the event-by-event analogues of the eWN profiles, are created using a Monte Carlo Glauber model. Using the same event plane method for obtaining $v_2$ as in the data analysis, we can reproduce both the measured centrality dependence and the $p_T$ shape of charged-particle elliptic flow up to $p_T\sim2$~GeV. We also consider the relation of elliptic flow to the initial state eccentricity using different reference planes, and discuss the correlation between the physical event plane and the initial participant plane. Our results demonstrate that event-by-event hydrodynamics with initial state fluctuations must be accounted for before a meaningful lower limit for viscosity can be obtained from elliptic flow data.

preprint2011arXiv

Inclusive prompt photon production in nuclear collisions at RHIC and LHC

Nuclear modification factors of inclusive prompt photon production in d-Au collisions at RHIC and p-Pb collisions at the LHC are provided at different rapidities. The calculations are performed at NLO accuracy using the EPS09 NLO nuclear parton distribution functions (nPDFs) and their error sets. The results are compared to the ones obtained with the nDS and HKN07 NLO nPDFs, and to the corresponding nuclear modification factors of neutral pion production in these collisions. The sensitivity of these results to the scale choice is also investigated. Interestingly, the predictions using the different nPDF sets differ from each other to the extent that this observable can be expected to become very useful for probing nPDFs over a wide range of Bjorken-x. In order to obtain a perturbative QCD baseline in heavy-ion collisions, calculations are carried out for minimum bias Au-Au collisions at RHIC and Pb-Pb collisions at the LHC. We also estimate the maximal possible suppression which the produced QCD matter can be expected to have on inclusive prompt photon production due to the quenching of the fragmentation component. The nuclear modification factor for prompt photon production is thus suggested to be used for gauging both the cold and the hot nuclear matter effects on other hard processes which are expected to be affected by quark-gluon plasma formation, such as large-pT hadron and jet production.

preprint2011arXiv

Monte Carlo Simulation for Elastic Energy Loss of Hard Partons in a Hydrodynamical Background

We have developed a Monte Carlo simulation describing the $2 \rightarrow 2$ scatterings of perturbatively produced, non-eikonally propagating high-energy partons with the quarks and gluons of the expanding QCD medium created in ultrarelativistic heavy ion collisions. The partonic scattering rates are computed in leading-order perturbative QCD (pQCD), while three different hydrodynamical scenarios are used to model the strongly interacting medium. We compare our results and tune the model with the neutral pion suppression observed in $\sqrt{s_{NN}}=200$ GeV Au+Au collisions at the BNL-RHIC. We find the incoherent nature of elastic energy loss incompatible with the measured angular dependence of the suppression. The effects of the initial state density fluctuations of the bulk medium are found to be small. Also the extrapolation from RHIC to the LHC is discussed.

preprint2011arXiv

Monte Carlo Simulation for Elastic Energy Loss of High-Energy Partons in Quark-Gluon Plasma

We examine the significance of $2 \rightarrow 2$ partonic collisions as the suppression mechanism of high-energy partons in the strongly interacting medium formed in ultrarelativistic heavy ion collisions. For this purpose, we have developed a Monte Carlo simulation describing the interactions of perturbatively produced, non-eikonally propagating high-energy partons with the quarks and gluons from the expanding QCD medium. The partonic collision rates are computed in leading-order perturbative QCD (pQCD), while three different hydrodynamical scenarios are used to model the medium. We compare our results with the suppression observed in $\sqrt{s_{NN}}=200$ GeV Au+Au collisions at the BNL-RHIC. We find the incoherent nature of elastic energy loss incompatible with the measured data and the effect of the initial state fluctuations small.

preprint2011arXiv

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 from 2.76 AGeV collisions, the available data probe much deeper into the model dynamics. A simultaneous fit of the nuclear suppression at both RHIC and LHC energies thus has great potential for discriminating between various models that yield equally good descriptions of RHIC data alone.

preprint2011arXiv

Systematics of the charged-hadron P_T spectrum and the nuclear suppression factor in heavy-ion collisions from sqrt{s}=200 GeV to sqrt{s} =2.76 TeV

In this paper, our goal is to make a simultaneous analysis of the high- and low-P_T parts of the charged-hadron P_T spectrum measured by the ALICE collaboration in central Pb-Pb collisions at sqrt{s}=2.76 TeV at the Large Hadron Collider (LHC), based on models which have been successfully applied and constrained in Au-Au collisions at Relativistic Heavy Ion Collider (RHIC). For the hydrodynamical modeling with which we obtain the low-P_T spectrum, we have computed the initial conditions based on perturbative QCD (pQCD) minijet production and saturation. The sensitivity of the obtained charged-hadron P_T spectrum on the hydrodynamic model parameters is studied. For the high-P_T part, we apply a number of parton-medium interaction models, which are tuned to describe the nuclear suppression factor R_AA measured at the RHIC in central Au-Au collisions at sqrt{s}=200 GeV. We find that the higher kinematic reach of the LHC, manifest in the hardening of the pQCD parton spectral slope, is in principle very efficient in discriminating the various models. However, due to the uncertainties in the p-p baseline, none of the tested models can be firmly ruled out with the present ALICE data. Comparison with the LHC data in this approach also shows that the matching of the hydrodynamic and pQCD+jet quenching components leaves fairly little room for other hadron production mechanisms in the cross-over region P_T=4-5 GeV.

preprint2011arXiv

Thermal photons from fluctuating initial conditions

Event-by-event fluctuations of initial QCD-matter density produced in heavy-ion collisions at RHIC enhance the production of thermal photons significantly in the region $2 \le p_T \le 4$ GeV/$c$ compared to a smooth initial-state averaged profile in the ideal hydrodynamic calculation. This enhancement is a an early time effect due to the presence of hotspots or over-dense regions in the fluctuating initial state. The effect of fluctuations is found to be stronger in peripheral than in central collisions.

preprint2010arXiv

A Monte-Carlo model for elastic energy loss in a hydrodynamical background

We present a computation of elastic energy loss of hard partons traversing the bulk hydrodynamical medium created in ultrarelativistic heavy-ion collisions. The model is based on perturbative Quantum Chromodynamics (pQCD) cross sections for $2\to 2$ processes in which a hard incoming parton is assumed to interact with a thermal parton from the medium. We model the interactions of this type in a Monte-Carlo framework to account properly for exact energy-momentum conservation, non-eikonal parton propagation, parton conversion reactions and the possibility to create additional hard recoiling partons from the medium. For the thermodynamical properties of the medium we use a hydrodynamical evolution model. We do not aim at a full description of high transverse momentum ($P_T$) observables at this point. Rather, we view the model as a starting point in obtaining a baseline of what to expect under the assumptions that the medium is describable by thermal quasifree partons and that their pQCD interactions with the high-energy partons are independent. Deviations from this baseline then call for more sophisticated medium description, as well as inclusion of higher-order processes and coherence effects in the pQCD scatterings.

preprint2010arXiv

Elastic energy loss with respect to the reaction plane in a Monte-Carlo model

We present a computation of $π^0$ nuclear modification factor with respect to the reaction plane in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV, based on a Monte-Carlo model of elastic energy loss of hard partons traversing the bulk hydrodynamical medium created in ultrarelativistic heavy-ion collisions. We find the incoherent nature of elastic energy loss incompatible with the measured data.

preprint2007arXiv

A global reanalysis of nuclear parton distribution functions

We determine the nuclear modifications of parton distribution functions of bound protons at scales $Q^2\ge 1.69$ GeV$^2$ and momentum fractions $10^{-5}\le x\le 1$ in a global analysis which utilizes nuclear hard process data, sum rules and leading-order DGLAP scale evolution. The main improvements over our earlier work {\em EKS98} are the automated $χ^2$ minimization, simplified and better controllable fit functions, and most importantly, the possibility for error estimates. The resulting 16-parameter fit to the N=514 datapoints is good, $χ^2/{\rm d.o.f}=0.82$. Within the error estimates obtained, the old {\em EKS98} parametrization is found to be fully consistent with the present analysis, with no essential difference in terms of $χ^2$ either. We also determine separate uncertainty bands for the nuclear gluon and sea quark modifications in the large-$x$ region where they are not stringently constrained by the available data. Comparison with other global analyses is shown and uncertainties demonstrated. Finally, we show that RHIC-BRAHMS data for inclusive hadron production in d+Au collisions lend support for a stronger gluon shadowing at $x<0.01$ and also that fairly large changes in the gluon modifications do not rapidly deteriorate the goodness of the overall fits, as long as the initial gluon modifications in the region $x\sim 0.02-0.04$ remain small.