Search

Knowledge of the pressure of cold and dense quark matter (QM) is known to significantly constrain the equation of state of neutron-star matter, a quantity playing a key role in deciphering the stars' internal structure. In this work, we make important progress towards determining the last unknown contribution to the pressure at Next-to-Next-to-Next-to-Leading Order (N3LO) in the strong coupling constant $\alpha_s$, available through the sum of all four-loop vacuum diagrams of dense Quantum Chromodynamics. We demonstrate the cancellation of both the covariant gauge parameter and all infrared (IR) divergences in the sum, showcasing the effective-field-theory paradigm in action. For the remaining IR-finite four-loop integrals, we demonstrate the efficacy of the dense Loop Tree Duality method -- a novel numerical framework for multiloop calculations in thermal field theory. Together, our results show that completing the N3LO pressure of cold QM is no longer merely possible but for the first time concretely within reach., Comment: 5 pages, 2 figures, 3 tables, appendices

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 correction to the high energy deep inelastic scattering (DIS) transverse structure function on a dense target in the dipole factorization framework. Together with our earlier result for longitudinal photons, this completes the calculation of the total deep inelastic scattering cross section in the dipole picture with massive quarks at next-to-leading order, enabling a comparison with experimental data., Comment: 58 pages, 8 figures

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., Comment: 10 pages, 7 figures. v2: references added/updated, typos fixed

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., Comment: v2: 42 pages, 5 figures, published version. v3: typos in eqs. 111 and 166 fixed

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_3\{2\}$ is larger than in the Pb+Pb system while $v_4\{2\}$ is of the same magnitude in both systems. We also find that the ratio of flow harmonics in Xe+Xe collisions and in Pb+Pb collisions shows a slight sensitivity to the temperature dependence of the shear-viscosity-to-entropy ratio. As we discuss here, the new nuclear mass-number systematics especially in the flow harmonics serves as a welcome further constraint for describing the space-time evolution of a heavy-ion system and for determining the shear viscosity and other transport properties of strongly interacting matter., Comment: 7 pages, 7 figures, revised version. The predicted flow coefficients vn with deformed Xe, and their ratios to 5.023 TeV Pb+Pb vn's are included as ancillary files

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 the four dimensional helicity scheme, agree with the recent calculation by Beuf using conventional dimensional regularization, confirming the regularization scheme independence of this cross section., Comment: v2: Improved version in response to referee's comments, matches the published version; 47 pages, 21 figures

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 measurements we then predict the charged hadron multiplicities and flow coefficients for the 5.023 TeV Pb+Pb collisions., Comment: 4 pages, Quark Matter 2015 proceedings

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\eta\big|_{|\eta|\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\eta\big|_{|\eta|\le 0.5} \propto s^{0.164\dots0.174}$., Comment: 8 pages, 7 figures

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., Comment: 5 pages, 4 figures, contributed talk by K.J.E. at the 7th International Conference on Hard and Electromagnetic Probes of High Energy Nuclear Collisions (Hard Probes 2015), Montreal, Canada, 29 June - 3 July, 2015

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., Comment: 32 pages, 23 figures

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., Comment: Prepared for the Proceedings of "Hard Probes 2013", the 6th International Conference on Hard and Electromagnetic Probes of High-Energy Nuclear Collisions, Stellenbosch, South Africa, November 4-8, 2013

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., Comment: International Conference on the Initial Stages of High-Energy Nuclear Collisions 2013, conference proceeding

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., Comment: 5 pages, 8 figures

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., Comment: 12 pages, 5 figures