We present a next-to-leading order (NLO) global DGLAP analysis of nuclear parton distribution functions (nPDFs) and their uncertainties. Carrying out an NLO nPDF analysis for the first time with three different types of experimental input -deep inelastic ℓ+A scattering, Drell-Yan dilepton production in p+A collisions, and inclusive pion production in d+Au and p+p collisions at RHIC -we find that these data can well be described in a conventional collinear factorization framework. Although the pion production has not been traditionally included in the global analyses, we find that the shape of the nuclear modification factor R dAu of the pion p T -spectrum at midrapidity retains sensitivity to the gluon distributions, providing evidence for shadowing and EMC-effect in the nuclear gluons. We use the Hessian method to quantify the nPDF uncertainties which originate from the uncertainties in the data. In this method the sensitivity of χ 2 to the variations of the fitting parameters is mapped out to orthogonal error sets which provide a user-friendly way to calculate how the nPDF uncertainties propagate to any factorizable nuclear cross-section. The obtained NLO and LO nPDFs and the corresponding error sets are collected in our new release called EPS09. These results should find applications in precision analyses of the signatures and properties of QCD matter at the LHC and RHIC.

Abstract. This White Paper presents the science case of an Electron-Ion Collider (EIC), focused on the structure and interactions of gluon-dominated matter, with the intent to articulate it to the broader nuclear science community. It was commissioned by the managements of Brookhaven National Laboratory (BNL) and Thomas Jefferson National Accelerator Facility (JLab) with the objective of presenting a summary of scientific opportunities and goals of the EIC as a follow-up to the 2007 NSAC Long Range plan. This document is a culmination of a community-wide effort in nuclear science following a series of workshops on EIC physics over the past decades and, in particular, the focused ten-week program on "Gluons and quark sea at high energies" at the Institute for Nuclear Theory in Fall 2010. It contains a brief description of a few golden physics measurements along with accelerator and detector concepts required to achieve them. It has been benefited profoundly from inputs by the users' communities of BNL and JLab. This White Paper offers the promise to propel the QCD science program in the US, established with the CEBAF accelerator at JLab and the RHIC collider at BNL, to the next QCD frontier. Preamble Editors' note for the second editionThe first edition of this White Paper was released in 2012. In the current (second) edition, the science case for the EIC is further sharpened in view of the recent data from BNL, CERN and JLab experiments and the lessons learnt from them. Additional improvements were made by taking into account suggestions from the larger nuclear physics community including those made at the EIC Users Group meeting at Stony Brook University in July 2014, and the QCD Town Meeting at Temple University in September 2014.Abhay Deshpande, Zein-Eddine Meziani and Jian-Wei Qiu November 2014 Editors' note for the third edition Since the 2nd release of this White Paper, the NSAC's Long Range Plan (2015) was successfully completed. The EIC is a major recommendation of the US nuclear science community. In the current release (version 3) we have fixed some minor remaining errors in the text, and have added a few new references. While the core science case for the EIC remains the same, the machine designs of both options, the eRHIC at BNL and the JLEIC at JLab keep evolving. In this 3rd release of the EIC White Paper instead of making substantial changes to the machine design sections (5.1 and 5.2), we give references to the most recent machine design documents.

We introduce a global analysis of collinearly factorized nuclear parton distribution functions (PDFs) including, for the first time, data constraints from LHC proton–lead collisions. In comparison to our previous analysis, EPS09, where data only from charged-lepton–nucleus deep inelastic scattering (DIS), Drell–Yan (DY) dilepton production in proton–nucleus collisions and inclusive pion production in deuteron–nucleus collisions were the input, we now increase the variety of data constraints to cover also neutrino–nucleus DIS and low-mass DY production in pion–nucleus collisions. The new LHC data significantly extend the kinematic reach of the data constraints. We now allow much more freedom for the flavor dependence of nuclear effects than in other currently available analyses. As a result, especially the uncertainty estimates are more objective flavor by flavor. The neutrino DIS plays a pivotal role in obtaining a mutually consistent behavior for both up and down valence quarks, and the LHC dijet data clearly constrain gluons at large momentum fraction. Mainly for insufficient statistics, the pion–nucleus DY and heavy-gauge-boson production in proton–lead collisions impose less visible constraints. The outcome – a new set of next-to-leading order nuclear PDFs called EPPS16 – is made available for applications in high-energy nuclear collisions.