[Paper Review] arXiv : Precision Monte Carlo Glauber predictions at present and future nuclear colliders

C. Loizides, D. d’Enterria+1 more|arXiv (Cornell University)|Oct 19, 2017

High-Energy Particle Collisions Research37 citations

TL;DR

This paper presents an improved Monte Carlo Glauber model for nuclear collisions at RHIC, LHC, and future hadron colliders, using data-driven inelastic pp cross sections and a refined nuclear density description with separate Fermi distributions for protons and neutrons. The model reduces uncertainties (mostly below 5%) and shows up to 7% changes in central Glauber quantities at √sNN = 5.02 TeV due to enhanced realism in nucleon spatial distribution and overlap calculations.

ABSTRACT

We present the results of an improved Monte Carlo Glauber (MCG) model of relevance for collisions involving nuclei at center-of-mass energies of BNL RHIC ($\sqrt{s_{ m NN}}=0.2$ TeV), CERN LHC ($\sqrt{s_{ m NN}}=2.76$-$8.8$ TeV), and proposed future hadron colliders ($\sqrt{s_{ m NN}}\approx 10$-$63$ TeV). The inelastic pp cross sections as a function of $\sqrt{s_{ m NN}}$ are obtained from a precise data-driven parametrization that exploits the many available measurements at LHC collision energies. We describe the nuclear transverse profile with two separated 2-parameter Fermi distributions for protons and neutrons to account for their different densities close to the nuclear periphery. Furthermore, we model the nucleon degrees of freedom inside the nucleus through a lattice with a minimum nodal separation, combined with a recentering and reweighting procedure, that overcomes some limitations of previous MCG approaches. The nuclear overlap function, number of participant nucleons and binary nucleon-nucleon collisions, participant eccentricity and triangularity, overlap area and average path length are presented in intervals of percentile centrality for lead-lead (PbPb) and proton-lead (pPb) collisions at all collision energies. We demonstrate for collisions at $\sqrt{s_{ m NN}}=5.02$ TeV that the central values of the Glauber quantities change by up to 7%, in a few bins of reaction centrality, due to the improvements implemented, though typically remain within the previously assigned systematic uncertainties, while their associated uncertainties are generally smaller (mostly below 5%) at all centralities than for earlier calculations. Tables for all quantities versus centrality at present and foreseen collision energies involving Pb nuclei, as well as for collisions of XeXe at $\sqrt{s_{ m NN}}=5.44$, and AuAu and CuCu at $\sqrt{s_{ m NN}}=0.2$ TeV, are provided.

Motivation & Objective

• To develop a more precise Monte Carlo Glauber model for nuclear collisions across current and future collider energies, including RHIC, LHC, and proposed 10–63 TeV hadron colliders.
• To reduce systematic uncertainties in Glauber-derived quantities such as participant nucleons, binary collisions, and eccentricity by improving the modeling of nuclear and nucleon density distributions.
• To incorporate data-driven parametrization of inelastic pp cross sections at LHC energies to ensure consistency with experimental measurements.
• To provide comprehensive tables of centrality-dependent Glauber observables for Pb, Xe, Au, and Cu nuclei at multiple energies for future experimental and theoretical reference.

Proposed method

• Utilizes a data-driven parametrization of inelastic pp cross sections as a function of √sNN, informed by LHC measurements across 2.76–8.8 TeV.
• Models the nuclear transverse profile using two separate 2-parameter Fermi distributions for protons and neutrons to reflect their differing densities near the nuclear periphery.
• Implements a lattice-based nucleon distribution with minimum nodal separation, followed by recentering and reweighting to better represent nucleon positions and avoid unphysical clustering.
• Calculates key Glauber observables—overlap function, participant nucleons, binary collisions, eccentricity, triangularity, overlap area, and average path length—across percentile centrality bins.
• Validates improvements by comparing results at √sNN = 5.02 TeV to previous calculations, quantifying changes and uncertainty reductions.
• Generates and provides detailed tables of all Glauber quantities for PbPb, pPb, XeXe, AuAu, and CuCu collisions at multiple energies.

Experimental results

Research questions

• RQ1 How do improved nucleon density modeling and data-driven pp cross sections affect centrality-dependent Glauber observables at √sNN = 5.02 TeV?
• RQ2 To what extent do the new model components reduce uncertainties in participant counts, binary collisions, and geometric eccentricities compared to previous Monte Carlo Glauber approaches?
• RQ3 What are the quantitative changes in Glauber quantities—such as participant eccentricity and overlap area—when using separate Fermi distributions for protons and neutrons?
• RQ4 How do the model predictions for PbPb and pPb collisions vary across the full range of present and future collider energies, from 0.2 TeV to 63 TeV?
• RQ5 What is the impact of the lattice-based nucleon placement with recentering and reweighting on the stability and accuracy of the Glauber calculations?

Key findings

• At √sNN = 5.02 TeV, the central values of Glauber quantities change by up to 7% in certain centrality bins due to model improvements, though most remain within prior systematic uncertainties.
• The uncertainties in all Glauber quantities are generally reduced to below 5% across all centralities, representing a significant improvement over earlier calculations.
• The use of separate 2-parameter Fermi distributions for protons and neutrons better captures the differing densities near the nuclear surface, improving the realism of the transverse nuclear profile.
• The lattice-based nucleon placement with nodal separation and reweighting procedure effectively mitigates clustering artifacts present in earlier MCG models.
• The model provides precise, data-constrained predictions for 146 unique combinations of nucleus, energy, and centrality, with full tables of observables for Pb, Xe, Au, and Cu systems.
• The improved model maintains consistency with experimental data while offering a more reliable and reduced-uncertainty baseline for interpreting heavy-ion collision data at RHIC and LHC.