Source author record

K. J. Eskola

K. 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

41works
6topics
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

41 published item(s)

preprint2016arXiv

Heavy ions at the Future Circular Collider

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.

preprint2016arXiv

Predictions for $p+$Pb Collisions at $\sqrt{s_{NN}} = 5$ TeV: Comparison with Data

Predictions made in Albacete {\it et al} prior to the LHC $p+$Pb run at $\sqrt{s_{NN}} = 5$ TeV are compared to currently available data. Some predictions shown here have been updated by including the same experimental cuts as the data. Some additional predictions are also presented, especially for quarkonia, that were provided to the experiments before the data were made public but were too late for the original publication are also shown here.

preprint2015arXiv

Event-by-event fluctuations in perturbative QCD + saturation + hydro model: pinning down 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 distributions of relative fluctuations of elliptic flow, and event-plane angle correlations, with the experimental data from Pb+Pb collisions at the LHC. We show how such a systematic multi-energy and multi-observable analysis tests the initial state calculation and the applicability region of hydrodynamics, and in particular how it constrains the temperature dependence of the shear viscosity-to-entropy ratio of QCD matter in its different phases in a remarkably consistent manner.

preprint2015arXiv

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 $η/s$ of QCD matter. Using these constraints from the current RHIC and LHC measurements we then predict the charged hadron multiplicities and flow coefficients for the 5.023 TeV Pb+Pb collisions.

preprint2015arXiv

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 data. Our systematic multi-energy and -observable analysis not only tests the initial state calculation and applicability of hydrodynamics, but also makes it possible to constrain the temperature dependence of the shear viscosity-to-entropy ratio, eta/s(T), of QCD matter in its different phases. Remarkably, we can describe all these different flow observables and correlations consistently with eta/s(T) that is independent of the collision energy.

preprint2015arXiv

Predictions for 5.023 TeV Pb+Pb collisions at the LHC

We compute predictions for various low-transverse-momentum bulk observables in $\sqrt{s_{NN}} = 5.023$ TeV Pb+Pb collisions at the 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 coefficients may, however, offer slightly more discriminating power on the temperature dependence of QCD matter viscosity than the 2.76 TeV measurements. Our prediction for the multiplicity in the 0-5\% centrality class, obtained using the two temperature-dependent shear-viscosity-to-entropy ratios that give the best overall fit to RHIC and LHC data is $dN_{\rm ch}/dη\big|_{|η|\le 0.5} =1876\dots2046$. We also predict a power-law increase from 200 GeV Au+Au collisions at RHIC to 2.76 and 5.023 TeV Pb+Pb collisions at the LHC, $dN_{\rm ch}/dη\big|_{|η|\le 0.5} \propto s^{0.164\dots0.174}$.

preprint2014arXiv

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.

preprint2013arXiv

A Large Hadron Electron Collider at CERN

This document provides a brief overview of the recently published report on the design of the Large Hadron Electron Collider (LHeC), which comprises its physics programme, accelerator physics, technology and main detector concepts. The LHeC exploits and develops challenging, though principally existing, accelerator and detector technologies. This summary is complemented by brief illustrations of some of the highlights of the physics programme, which relies on a vastly extended kinematic range, luminosity and unprecedented precision in deep inelastic scattering. Illustrations are provided regarding high precision QCD, new physics (Higgs, SUSY) and electron-ion physics. The LHeC is designed to run synchronously with the LHC in the twenties and to achieve an integrated luminosity of O(100) fb$^{-1}$. It will become the cleanest high resolution microscope of mankind and will substantially extend as well as complement the investigation of the physics of the TeV energy scale, which has been enabled by the LHC.

preprint2013arXiv

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 constrained in this new framework simultaneously by all these data.

preprint2013arXiv

On the Relation of the LHeC and the LHC

The present note relies on the recently published conceptual design report of the LHeC and extends the first contribution to the European strategy debate in emphasising the role of the LHeC to complement and complete the high luminosity LHC programme. The brief discussion therefore focuses on the importance of high precision PDF and $α_s$ determinations for the physics beyond the Standard Model (GUTs, SUSY, Higgs). Emphasis is also given to the importance of high parton density phenomena in nuclei and their relevance to the heavy ion physics programme at the LHC.

preprint2012arXiv

A Large Hadron Electron Collider at CERN: Report on the Physics and Design Concepts for Machine and Detector

The physics programme and the design are described of a new collider for particle and nuclear physics, the Large Hadron Electron Collider (LHeC), in which a newly built electron beam of 60 GeV, up to possibly 140 GeV, energy collides with the intense hadron beams of the LHC. Compared to HERA, the kinematic range covered is extended by a factor of twenty in the negative four-momentum squared, $Q^2$, and in the inverse Bjorken $x$, while with the design luminosity of $10^{33}$ cm$^{-2}$s$^{-1}$ the LHeC is projected to exceed the integrated HERA luminosity by two orders of magnitude. The physics programme is devoted to an exploration of the energy frontier, complementing the LHC and its discovery potential for physics beyond the Standard Model with high precision deep inelastic scattering measurements. These are designed to investigate a variety of fundamental questions in strong and electroweak interactions. The physics programme also includes electron-deuteron and electron-ion scattering in a $(Q^2, 1/x)$ range extended by four orders of magnitude as compared to previous lepton-nucleus DIS experiments for novel investigations of neutron's and nuclear structure, the initial conditions of Quark-Gluon Plasma formation and further quantum chromodynamic phenomena. The LHeC may be realised either as a ring-ring or as a linac-ring collider. Optics and beam dynamics studies are presented for both versions, along with technical design considerations on the interaction region, magnets and further components, together with a design study for a high acceptance detector. Civil engineering and installation studies are presented for the accelerator and the detector. The LHeC can be built within a decade and thus be operated while the LHC runs in its high-luminosity phase. It thus represents a major opportunity for progress in particle physics exploiting the investment made in the LHC.

preprint2012arXiv

Initial state anisotropies and their uncertainties in ultrarelativistic heavy-ion collisions from the Monte Carlo Glauber model

In hydrodynamical modeling of heavy-ion collisions, the initial-state spatial anisotropies are translated into momentum anisotropies of the final-state particle distributions. Thus, understanding the origin of the initial-state anisotropies and their uncertainties is important before extracting specific QCD matter properties, such as viscosity, from the experimental data. In this work we review the wounded nucleon approach based on the Monte Carlo Glauber model, charting in particular the uncertainties arising from modeling of the nucleon-nucleon interactions between the colliding nucleon pairs and nucleon-nucleon correlations inside the colliding nuclei. We discuss the differences between the black disk model and a probabilistic profile function approach for the inelastic nucleon-nucleon interactions, and investigate the influence of initial-state correlations using state-of-the-art modeling of these.

preprint2012arXiv

Initial state anisotropies in ultrarelativistic heavy-ion collisions from the Monte Carlo Glauber model

In hydrodynamicalmodeling of heavy-ion collisions the initial state spatial anisotropies translate into momentum anisotropies of the final state particle distributions. Thus, understanding the origin of the initial anisotropies and quantifying their uncertainties is important for the extraction of specific QCD matter properties, such as viscosity, from the experimental data. In this work we study the wounded nucleon approach in the Monte Carlo Glauber model framework, focusing especially on the uncertainties which arise from the modeling of the nucleon-nucleon interactions between the colliding nucleon pairs and nucleon-nucleon correlations inside the colliding nuclei. We compare the black disk model and a probabilistic profile function approach for the inelastic nucleon-nucleon interactions, and study the effects of initial state correlations using state-of-theart modeling of these.

preprint2012arXiv

Multiplicities and $p_T$ spectra 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 computed $E_T$, saturation scales and the final-state multiplicities. The key parameters, which need to be fixed using the measurements, are identified, and their correlation is discussed. We convert the saturated minijet $E_T$ into QCD-matter initial conditions for longitudinally boost-invariant ideal hydrodynamics. We compute the charged-particle multiplicities and identified bulk hadron $p_T$ spectra in 5% most central Au+Au collisions at RHIC and Pb+Pb at the LHC. We obtain an encouragingly good agreement with the experimental data, simultaneously at RHIC and LHC, showing that the approach has a definite predictive power.

preprint2010arXiv

Nuclear PDFs at NLO - status report and review of the EPS09 results

We review the current status of the global DGLAP analysis of nuclear parton distribution functions, nPDFs, focusing on the recent EPS09 analysis, whose output, EPS09NLO, is the best-constrained NLO nPDF set on the market. Collinear factorization is found to work very well in the kinematical region studied. With the error sets released in the EPS09 package one can compute how the nPDF-related uncertainties propagate into factorizable nuclear hard-process cross sections. A comparison with the other existing NLO nPDF sets is shown, and the BRAHMS forward-$η$ hadron data from d+Au collisions are discussed in the light of the EPS09 nPDFs and their error sets.

preprint2009arXiv

EPS09 - Global NLO analysis of nuclear PDFs and their uncertainties

In this talk, we introduce our recently completed next-to-leading order (NLO) global analysis of the nuclear parton distribution functions (nPDFs) called EPS09 - a higher order successor to the well-known leading-order (LO) analysis EKS98 and also to our previous LO work EPS08. As an extension to similar global analyses carried out by other groups, we complement the data from deep inelastic $l+A$ scattering and Drell-Yan dilepton measurements in p+$A$ collisions by inclusive midrapidity pion production data from d+Au collisions at RHIC, which results in better constrained gluon distributions than before. The most important new ingredient, however, is the detailed error analysis, which employs the Hessian method and which allows us to map out the parameter-space vicinity of the best-fit to a collection of nPDF error sets. These error sets provide the end-user a way to compute how the PDF-uncertainties will propagate into the cross sections of his/her interest. The EPS09 package to be released soon, will contain both the NLO and LO results for the best fits and the uncertainty sets.

preprint2008arXiv

An improved global analysis of nuclear parton distribution functions including RHIC data

We present an improved leading-order global DGLAP analysis of nuclear parton distribution functions (nPDFs), supplementing the traditionally used data from deep inelastic lepton-nucleus scattering and Drell-Yan dilepton production in proton-nucleus collisions, with inclusive high-$p_T$ hadron production data measured at RHIC in d+Au collisions. With the help of an extended definition of the $χ^2$ function, we now can more efficiently exploit the constraints the different data sets offer, for gluon shadowing in particular, and account for the overall data normalization uncertainties during the automated $χ^2$ minimization. The very good simultaneous fit to the nuclear hard process data used demonstrates the feasibility of a universal set of nPDFs, but also limitations become visible. The high-$p_T$ forward-rapidity hadron data of BRAHMS add a new crucial constraint into the analysis by offering a direct probe for the nuclear gluon distributions -- a sector in the nPDFs which has traditionally been very badly constrained. We obtain a strikingly stronger gluon shadowing than what has been estimated in previous global analyses. The obtained nPDFs are released as a parametrization called EPS08.

preprint2005arXiv

Transverse Spectra of Hadrons in Central $AA$ Collisions at RHIC and LHC from pQCD+Saturation+Hydrodynamics and from pQCD+Energy Losses

We study the transverse spectra of hadrons in nearly central $AA$ collisions at RHIC and LHC in a broad transverse momentum range Low-$p_T$ spectra are calculated by using boost-invariant hydrodynamics with initial energy and net-baryon densities from the EKRT pQCD+saturation model. High-$p_T$ spectra are obtained from pQCD jet calculation including the energy loss of the parton in the matter prior to its fragmentation to final hadrons.

preprint2004arXiv

D meson enhancement in pp collisions at the LHC due to nonlinear gluon evolution

When nonlinear effects on the gluon evolution are included with constraints from HERA, the gluon distribution in the free proton is enhanced at low momentum fractions, x < 0.01, and low scales, Q^2 < 10 GeV^2, relative to standard, DGLAP-evolved, gluon distributions. Consequently, such gluon distributions can enhance charm production in pp collisions at center of mass energy 14 TeV by up to a factor of five at midrapidity, y \sim 0, and transverse momentum p_T -> 0 in the most optimistic case. We show that most of this enhancement survives hadronization into D mesons. Assuming the same enhancement at leading and next-to-leading order, we show that the D enhancement may be measured by D^0 reconstruction in the K^-π^+ decay channel with the ALICE detector.

preprint2004arXiv

Enhanced charm hadroproduction due to nonlinear corrections to the DGLAP equations

We have studied the effects of nonlinear scale evolution of the parton distribution functions to charm production in $pp$ collisions at center-of-mass energies of 5.5, 8.8 and 14 TeV. We find that the differential charm cross section can be enhanced up to a factor of 4-5 at low $p_T$. The enhancement is quite sensitive to the charm quark mass and the renormalization/factorization scales.

preprint2003arXiv

Enhancement of charm quark production due to nonlinear corrections to the DGLAP equations

We have studied how parton distributions based on the inclusion of nonlinear scale evolution and constraints from HERA data affect charm production in $pp$ collisions at center-of-mass energies of 5.5, 8.8 and 14 TeV. We find that, while the resulting enhancement can be substantial, it is very sensitive to the charm quark mass and the scale entering the parton densities and the strong coupling constant.

preprint2003arXiv

Nonlinear corrections to the DGLAP equations in view of the HERA data

The effects of the first nonlinear corrections to the DGLAP evolution equations are studied by using the recent HERA data for the structure function $F_2(x,Q^2)$ of the free proton and the parton distributions from CTEQ5L and CTEQ6L as a baseline. By requiring a good fit to the H1 data, we determine initial parton distributions at $Q_0^2=1.4$ GeV$^2$ for the nonlinear scale evolution. We show that the nonlinear corrections improve the agreement with the $F_2(x,Q^2)$ data in the region of $x\sim 3\cdot 10^{-5}$ and $Q^2\sim 1.5$ GeV$^2$ without paying the price of obtaining a worse agreement at larger values of $x$ and $Q^2$. For the gluon distribution the nonlinear effects are found to play an increasingly important role at $x\lsim 10^{-3}$ and $Q^2\lsim10$ GeV$^2$, but rapidly vanish at larger values of $x$ and $Q^2$. Consequently, contrary to CTEQ6L, the obtained gluon distribution at $Q^2=1.4$ GeV$^2$ shows a power-like growth at small $x$. Relative to the CTEQ6L gluons, an enhancement up to a factor $\sim6$ at $x=10^{-5}$, $Q_0^2=1.4$ GeV$^2$ reduces to a negligible difference at $Q^2\gsim 10$ GeV$^2$.

preprint2002arXiv

A perturbative QCD analysis of charged-particle distributions in hadronic and nuclear collisions

We compute the distributions of charged particles at large transverse momenta in $p\bar p(p)$, $pA$ and $AA$ collisions in the framework of perturbative QCD, by using collinear factorization and the modern PDFs and fragmentation functions. At the highest cms-energies the shape of the spectra measured in $p\bar p(p)$ collisions at large $q_T$ can be well explained. The difference between the data and the lowest-order computation is quantified in terms of a constant $K$-factor for each energy. The $K$-factor is found to systematically decrease with growing $\sqrt s$. Also a lower limit for the partonic transverse momentum, $p_0$, is extracted for each $\sqrt s$ based on the comparison with the measurements. A systematic increase of $p_0$ as a function of $\sqrt s$ is found. Nuclear effects in the charged-particle spectra in $pA$ and $AA$ collisions at RHIC and LHC are studied in the framework of collinear factorization by applying the EKS98 nuclear corrections to the parton distributions. The nuclear effects are shown to mostly enhance the computed spectra. A comparison with the recent PHENIX data from central and peripheral Au+Au collisions at RHIC is done.

preprint2002arXiv

Rapidity dependence of particle production in ultrarelativistic nuclear collisions

We compute the rapidity dependence of particle and transverse energy production in ultrarelativistic heavy ion collisions at various beam energies and atomic numbers using the perturbative QCD + saturation model. The distribution is a broad gaussian near $y=0$ but the rapid increase of particle production with the beam energy will via energy conservation strongly constrain the rapidity distribution at large $y$.

preprint2001arXiv

Heavy ion collision multiplicities and gluon distribution functions

Atomic number ($A$) and energy ($\roots$) scaling exponents of multiplicity and transverse energy in heavy ion collisions are analytically derived in the perturbative QCD + saturation model. The exponents depend on the small-$x$ behaviour of gluon distribution functions at an $x$-dependent scale. The relation between initial state and final state saturation is also discussed.

preprint2001arXiv

Multiplicities and Transverse Energies in Central AA Collisions at RHIC and LHC from pQCD, Saturation and Hydrodynamics

We compute the particle multiplicities and transverse energies at central and nearly central AA collisions at RHIC and LHC. The initial state is computed from perturbative QCD supplemented by the conjecture of saturation of produced partons. The expansion stage is described in terms of hydrodynamics assuming longitudinal boost invariance and azimuthal symmetry. Transverse flow effects, a realistic list of hadrons and resonance decays are included. Comparison with the data of the multiplicities at $\sqrt s=56$ AGeV and 130 AGeV from RHIC is done and predictions for the full RHIC energy and LHC energy are made for the multiplicities and transverse energies. The reduction from the initially released minijet transverse energy to the $E_T$ in the final state is less than in the one-dimensional case but still dramatic: a factor of 2.7 at RHIC, and 3.6 at the LHC.

preprint2000arXiv

Production of transverse energy from minijets in next-to-leading order perturbative QCD

We compute in next-to-leading order (NLO) perturbative QCD the transverse energy carried into the central rapidity unit of hadron or nuclear collisions by the partons freed in the few-GeV subcollisions. The formulation is based on a rapidity window and a measurement function of a new type. The behaviour of the NLO results as a function of the minimum transverse momentum and as a function of the scale choice is studied. The NLO results are found to be stable relative to the leading-order ones even in the few-GeV domain.

preprint2000arXiv

Transverse energy from minijets in ultrarelativistic nuclear collisions: a next-to-leading order analysis

We compute in next-to-leading order (NLO) perturbative QCD the amount of transverse energy produced into a rapidity region $ΔY$ of a nuclear collision from partons created in the few-GeV subcollisions. The NLO formulation assumes collinear factorization and is based on the subtraction method. We first study the results as a function of the minimum transverse momentum scale and define and determine the associated $K$-factors. The dependence of the NLO results on the scale choice and on the size of $ΔY$ is also studied. The calculations are performed for GRV94 and CTEQ5 sets of parton distributions. Also the effect of nuclear shadowing to the NLO results is investigated. The main results are that the NLO results are stable relative to the leading-order (LO) ones even in the few-GeV domain and that there exists a new kinematical region not present in the LO, which is responsible for the magnitude of the $K$-factors. The dependence on the size of $ΔY$ is found to be practically linear within two central units of rapidity, and the scale dependence present at RHIC energies is seen to vanish at the LHC energies. Finally, the effect of NLO corrections to the determination of transverse energies and multiplicities in nuclear collisions is discussed within the framework of parton saturation.

preprint1999arXiv

Nuclear Parton Distributions - a DGLAP Analysis

Nuclear parton distributions $f_A(x,Q^2)$ are studied within a framework of the DGLAP evolution. Measurements of $F_2^A/F_2^D$ in deep inelastic $lA$ collisions, and Drell--Yan dilepton cross sections measured in $pA$ collisions are used as constraints. Also conservation of momentum and baryon number is required. It is shown that the calculated $Q^2$ evolution of $F_2^{\rm Sn}/F_2^{\rm C}$ agrees very well with the recent NMC data, and that the ratios $R_f^A=f_A/f$ are only moderately sensitive to the choice of a specific modern set of free parton distributions. For general use, we offer a numerical parametrization of $R_f^A(x,Q^2)$ for all parton flavours $f$ in $A>2$, and at $10^{-6}\le x \le 1$ and $2.25 {\rm GeV}^2\le Q^2\le 10^4 {\rm GeV}^2$.

preprint1999arXiv

Scaling of transverse energies and multiplicities with atomic number and energy in ultrarelativistic nuclear collisions

We compute how the initial energy density and produced gluon, quark and antiquark numbers scale with atomic number and beam energy in ultrarelativistic heavy ion collisions. The computation is based on the argument that the effect of all momentum scales can be estimated by performing the computation at one transverse momentum scale, the saturation momentum. The initial numbers are converted to final ones by assuming kinetic thermalisation and adiabatic expansion. The main emphasis of the study is at LHC and RHIC energies but it is observed that even at SPS energies this approach leads to results which are not unreasonable: what is usually described as a completely soft nonperturbative process can also be described in terms of gluons and quarks. The key element is the use of the saturation scale.

preprint1998arXiv

Scale evolution of nuclear parton distributions

Using the NMC and E665 nuclear structure function ratios $F_2^A/F_2^D$ and $F_2^A/F_2^{C}$ from deep inelastic lepton-nucleus collisions, and the E772 Drell--Yan dilepton cross sections from proton-nucleus collisions, and incorporating baryon number and momentum sum rules, we determine nuclear parton distributions at an initial scale $Q_0^2$. With these distributions, we study QCD scale evolution of nuclear parton densities. The emphasis is on small values of $x$, especially on scale dependence of nuclear shadowing. As the main result, we show that a consistent picture can be obtained within the leading twist DGLAP evolution, and in particular, that the calculated $Q^2$ dependence of $F_2^{Sn}/F_2^{C}$ agrees very well with the recent NMC data.

preprint1998arXiv

The scale dependent nuclear effects in parton distributions for practical applications

The scale dependence of the ratios of parton distributions in a proton of a nucleus $A$ and in the free proton, $R_i^A(x,Q^2)=f_{i/A}(x,Q^2)/f_i(x,Q^2)$, is studied within the framework of the lowest order leading-twist DGLAP evolution. By evolving the initial nuclear distributions obtained with the GRV-LO and CTEQ4L sets at a scale $Q_0^2$, we show that the ratios $R_i^A(x,Q^2)$ are only moderately sensitive to the choice of a specific modern set of free parton distributions. We propose that to a good first approximation, this parton distribution set-dependence of the nuclear ratios $R_i^A(x,Q^2)$ can be neglected in practical applications. With this result, we offer a numerical parametrization of $R_i^A(x,Q^2)$ for all parton flavours $i$ in any $A>2$, and at any $10^{-6}\le x \le 1$ and any $Q^2\ge 2.25$ GeV$^2$ for computing cross sections of hard processes in nuclear collisions.

preprint1997arXiv

Hydrodynamics of nuclear collisions with initial conditions from perturbative QCD

We compute the longitudinal hydrodynamic flow in ultrarelativistic heavy ion collisions at $\sqrt{s}$ = 5500 GeV by using boost non-invariant initial conditions following from perturbative QCD. The transfer of entropy and energy from the central region to larger rapidities caused by boost non-invariance is determined and the associated decrease in the lifetime of the system is estimated.

preprint1996arXiv

Baryon-to-entropy ratio in very high energy nuclear collisions

We compute as a function of rapidity $y$ the baryon number carried by quarks and antiquarks with $p_T > p_0 \approx$ 2 GeV produced in Pb+Pb collisions at TeV energies. The computation is carried out in lowest order QCD perturbation theory using structure functions compatible with HERA results. At $p_0=2$ GeV the initial gluon density is both transversally saturated and thermalised in the sense that the energy/gluon equals to that of an ideal gas with the same energy density. Even at these high energies the initial net baryon number density at $y=0$ at $τ=0.1$ fm will be more than the normal nuclear matter density but the baryon-to-entropy ratio is only $(B-\bar B)/S\sim 1/5000$. Further evolution of the system is discussed and the final baryon-to-entropy ratio is estimated.

preprint1995arXiv

Formation and Evolution of Quark-Gluon Plasma at RHIC and LHC

Initial conditions for quark-gluon plasma formation at $τ=0.1\fm$ are considered in lowest order perturbative QCD. Chemical composition, thermalization of the formed semihard quark-gluon system and especially implications of the new HERA parton distributions with the enhancement at small $x$ are studied. The plasma at $τ=0.1\fm$ is shown to be strongly gluon dominated both at RHIC and LHC, and a possibility for rapid thermalization of gluons at LHC is pointed out. Uncertainties in the calculations, particularly shadowing corrections to the parton distributions, are discussed. Free streaming and ideal hydro limits for the evolution of the gluonic plasma with the calculated minijet initial conditions are demonstrated, and a lower limit for final multiplicities obtained for the LHC nuclear collisions.

preprint1993arXiv

Perturbative Gluon Shadowing in Heavy Nuclei

We study how much gluon shadowing can be perturbatively generated through the modified QCD evolution in heavy nuclei. The evolution of small-$x$ gluons is investigated within the semiclassical approximation. The method of characteristics is used to evaluate the shadowed distributions in low-$Q$ and small-$x$ region. In solving the modified evolution equation, we model in simultaneously fusions from independent constituents and from the same constituent, both in a proton and in a large loosely bound nucleus of $A\sim 200$. In addition to the actual distributions at small $x$, we study the ratios of the distributions at an initial scale $Q_0 = 2$ GeV, and show that a strong nuclear shadowing can follow from the modified QCD evolution.