Charged-Particle Pseudorapidity Densities

• Charged-particle pseudorapidity densities are a metric that quantifies the number of primary charged particles produced as a function of pseudorapidity in high-energy collisions.
• Measurement techniques include silicon tracker hit counting, pixel-based tracklet reconstruction, and energy-deposit methods in forward detectors to ensure precise event characterization.
• Analyses of these densities uncover key trends in energy, centrality, and system-size scaling, which help constrain soft QCD mechanisms and probe quark–gluon plasma properties.

Charged-particle pseudorapidity densities, commonly denoted as $\frac{dN_\mathrm{ch}}{d\eta}$, quantify the number of primary charged particles produced per event in high-energy collisions as a function of pseudorapidity $\eta$. They serve as foundational observables for characterizing the longitudinal profile of particle production in hadronic, nuclear, and, more generally, multiparticle production processes. Analyses of $\frac{dN_\mathrm{ch}}{d\eta}$ in varied systems and energies enable constraints on soft and hard QCD mechanisms, inform energy density estimates, and probe the initial-state geometry and collective dynamics.

1. Experimental Determination of Pseudorapidity Densities

Charged-particle pseudorapidity densities are measured by counting primary charged particles in well-characterized detector regions as a function of $\eta$, where

$\eta = -\ln\left[\tan\left(\theta/2\right)\right]$

is defined by the polar angle $\theta$ with respect to the beam axis.

Techniques for extracting $\frac{dN_\mathrm{ch}}{d\eta}$ depend on detector architecture and include:

• Silicon tracker multiplicity arrays (e.g., PHOBOS at RHIC): Hit counting in single-layer silicon sensors, with analog corrections for multiply-charged fragments, over a broad pseudorapidity interval ($|\eta| < 5.4$). Centrality determination is based on auxiliary scintillator counters and pathlength-corrected calorimetric energy deposition, cross-calibrated with a Glauber model to obtain $\langle N_\mathrm{part} \rangle$ (0806.2803).
• Pixel detector-based tracklet reconstruction (e.g., ALICE, CMS): Short track segments (“tracklets”) are built from hits in consecutive pixel layers, accepted according to angular quality cuts, and corrected for combinatorial background, primary/secondary particle separation, and effective detector acceptance, yielding high efficiency and an effective $p_T$ threshold down to $\eta$0 MeV/$\eta$1 (Modak, 2024, Collaboration, 2010).
• Statistical/energy-deposit methods in forward detectors (e.g., FMD in ALICE, T2 in TOTEM): Multiplicities in forward/backward regions are estimated by relating energy deposition to mean charged-particle number using fitting or clustering techniques and statistical correction for backgrounds (Modak, 2022, Collaboration et al., 2014).
• Vertexing and event selection: Precise primary vertex determination optimizes acceptance, while trigger and offline selection criteria (using, e.g., FIT or V0 counters) ensure robust distinction of event classes (minimum-bias, inelastic, NSD, etc.) (Collaboration, 3 Apr 2025, Modak, 2024).

Correction factors account for acceptance, efficiency, bin migration, contamination by secondary particles, and possible multiplicity biases from triggers or centrality estimators.

2. Energy, Centrality, and System-Size Dependence

The energy, centrality, and system-size dependence of $\eta$2 encodes crucial information about soft QCD and initial geometric effects.

• Energy dependence: In proton–proton (pp) collisions, the midrapidity charged-particle density rises with energy as a power law,

$\eta$3

where $\eta$4–$\eta$5 in pp and increases to $\eta$6–$\eta$7 in central nucleus–nucleus (A–A) collisions, reflecting the increased phase-space and possible onset of additional particle production mechanisms (Collaboration, 2022, Modak, 2024).

• Centrality and system-size scaling: In A–A collisions, $\eta$8 increases with centrality (smaller impact parameter, larger overlap) and scales approximately with the number of participant nucleons $\eta$9, especially at low/intermediate energies and in peripheral collisions. The normalized quantity,

$\frac{dN_\mathrm{ch}}{d\eta}$0

often exhibits a weak centrality dependence from central to semi-peripheral collisions and only a modest decrease (factor $\frac{dN_\mathrm{ch}}{d\eta}$1–$\frac{dN_\mathrm{ch}}{d\eta}$2) from most central to most peripheral events (Collaboration, 3 Apr 2025, Collaboration, 2015).

• System-size hierarchy: Systematic comparisons across pp, p–A, and A–A at the same $\frac{dN_\mathrm{ch}}{d\eta}$3 using one detector (ALICE) reveal a smooth increase of midrapidity $\frac{dN_\mathrm{ch}}{d\eta}$4 with system size. Central Pb–Pb yields reach $\frac{dN_\mathrm{ch}}{d\eta}$5 per unit $\frac{dN_\mathrm{ch}}{d\eta}$6 at 5.02–5.36 TeV; in central p–Pb, $\frac{dN_\mathrm{ch}}{d\eta}$7; and in pp, $\frac{dN_\mathrm{ch}}{d\eta}$8 (Collaboration, 2022).

Tabulated energy and system-size scaling:

System	Energy (TeV)	$\frac{dN_\mathrm{ch}}{d\eta}$9 (midrapidity)
pp (INEL>0)	13.6	$\eta$0
Pb–Pb (0–5%)	5.36	$\eta$1
p–Pb (NSD)	8.16	$\eta$2

3. Geometric and Participant Scaling: Nuclear Overlap and Implications

The geometry of the nuclear overlap zone and the participant scaling variables define the scaling of the charged-particle production:

• Participant number vs geometry: The overall multiplicity and shape (height and width) of $\eta$3 in A–A and smaller systems are driven more by the geometry of the collision zone (quantified via $\eta$4, where $\eta$5 is the nuclear mass number) than by $\eta$6 alone. Collisions with the same participant fraction $\eta$7 from different nuclei yield more similar pseudorapidity distributions over the full range of $\eta$8 (0806.2803).
• Longitudinal scaling: Plotting $\eta$9 as a function of $\eta = -\ln\left[\tan\left(\theta/2\right)\right]$0 reveals longitudinal (“limiting fragmentation”) scaling at forward/backward rapidities across a broad energy range, consistent with an invariance of the fragmentation region as $\eta = -\ln\left[\tan\left(\theta/2\right)\right]$1 increases. This behavior persists in inclusive heavy-ion collisions from RHIC to LHC energies (Collaboration, 2013, Basu et al., 2020).
• Violations at high energy: At LHC energies, the per-participant-pair density at midrapidity grows with centrality, violating simple participant scaling and indicating increasing contributions from small-$\eta = -\ln\left[\tan\left(\theta/2\right)\right]$2 gluon production/saturation (Basu et al., 2020).

4. Theoretical Interpretation and Model Comparisons

Charged-particle pseudorapidity densities serve as critical constraints for QCD-based models of multiparticle production and medium formation:

• Soft–hard two-component models: Particle production is often modeled as a superposition of “soft” (participant-dependent) and “hard” (binary collision, $\eta = -\ln\left[\tan\left(\theta/2\right)\right]$3-dependent) components. The centrality and energy trends favor a dominantly participant-driven (“soft”) scaling, especially at RHIC and lower LHC energies ($\eta = -\ln\left[\tan\left(\theta/2\right)\right]$4 per participant pair approximately independent of centrality), but deviations appear as energy increases (Basu et al., 2020).
• Color glass condensate (CGC)/gluon saturation: Enhanced centrality dependence and the observed shape are more faithfully captured in models involving initial-state gluon saturation, with pseudorapidity densities parameterized as $\eta = -\ln\left[\tan\left(\theta/2\right)\right]$5 where $\eta = -\ln\left[\tan\left(\theta/2\right)\right]$6 is larger for A–A than for pp (Basu et al., 2020, Collaboration, 2012, Collaboration, 2018).
• Hydrodynamic models: Hydrodynamical approaches, especially when incorporating initial-state conditions from color domain models (e.g., IP-Glasma+MUSIC+UrQMD), reproduce the absolute level and broad shape of the most central Pb–Pb multiplicities. Nevertheless, many models tend to underestimate the observed multiplicities at the highest energies and centralities (Collaboration, 3 Apr 2025, Modak, 2024).
• Rapidity width and Landau model: The rapidity distribution width, $\eta = -\ln\left[\tan\left(\theta/2\right)\right]$7, extracted from a Gaussian parameterization of $\eta = -\ln\left[\tan\left(\theta/2\right)\right]$8, increases with energy and system size, exceeding expectations from Landau hydrodynamics. Its proportionality to the beam rapidity $\eta = -\ln\left[\tan\left(\theta/2\right)\right]$9 reflects the dominance of available phase-space constraints at high energy (Collaboration, 2016, Collaboration, 2013).

5. Pseudorapidity Densities as Probes of the Initial State and Medium Properties

• Initial energy density estimates: Through the Bjorken energy density formula,

$\theta$0

where $\theta$1 is the overlap area and $\theta$2 the formation time, $\theta$3 constrains $\theta$4 via measured charged-particle multiplicities and mean $\theta$5. The rise in $\theta$6 with system size and centrality reflects a tenfold boost in initial energy density from pp to Pb–Pb at fixed $\theta$7, consistent with the formation of a strongly coupled quark–gluon plasma (QGP) phase (Collaboration, 2022).

• Collective phenomena: Particle production at midrapidity in central A–A collisions, alongside observed collective flow patterns, signals rapid thermalization and isotropization. Comparisons of $\theta$8 per participant across systems further inform whether collective effects extend to small systems (p–A, high-multiplicity pp) (Modak, 2022).

6. Scaling, Extrapolation, and Future Directions

• Phenomenological parameterizations: Three functional forms are commonly used for $\theta$9 extrapolation: trapezoidal, sum of two Gaussians, and a difference of Gaussians. Across experiment types, the difference of Gaussians achieves the most robust simultaneous description for pp, p–A, and A–A (Basu et al., 2020, Collaboration, 2013).
• Extrapolation and projections: By fitting the energy dependence of central A–A $\frac{dN_\mathrm{ch}}{d\eta}$0 per participant pair as $\frac{dN_\mathrm{ch}}{d\eta}$1, anticipated densities at midrapidity reach $\frac{dN_\mathrm{ch}}{d\eta}$2 for Pb–Pb at 39 TeV (FCC) (Basu et al., 2020).
• Instrumental advances: Upgrades to vertexing (ALPIDE MAPS), tracking (GEM-based TPC), and forward acceptance (MFT, FMD) allow extension of $\frac{dN_\mathrm{ch}}{d\eta}$3 coverage, enhanced centrality resolution, and systematic reduction of uncertainties, supporting the next generation of precision measurements (Tripathy, 17 Jun 2025, Modak, 2024).

7. Summary Table: Core Quantities and Scaling Expressions

Observable	Definition / Formula	Observational Context
Pseudorapidity $\frac{dN_\mathrm{ch}}{d\eta}$4	$\frac{dN_\mathrm{ch}}{d\eta}$5	Universally adopted for longitudinal kinematics
Density $\frac{dN_\mathrm{ch}}{d\eta}$6	Primary charged particles per unit $\frac{dN_\mathrm{ch}}{d\eta}$7	Central to all multiplicity analyses
Normalized density	$\frac{dN_\mathrm{ch}}{d\eta}$8	System size, centrality scaling
Longitudinal shift	$\frac{dN_\mathrm{ch}}{d\eta}$9	For limiting fragmentation studies
Energy scaling	$|\eta| < 5.4$0	$|\eta| < 5.4$1 (pp), $|\eta| < 5.4$2 (A–A, central)
Energy density	$|\eta| < 5.4$3 via $|\eta| < 5.4$4	Centrality, system-size dependence

In sum, charged-particle pseudorapidity densities function as a principal diagnostic for multiparticle production mechanisms, initial geometry, energy scaling, and collective effects in high-energy nuclear and hadronic collisions. The scaling behavior with energy, system size, and centrality—combined with new advancements in detector technology and robust, model-agnostic extrapolation schemes—continues to provide critical constraints on the dynamics of QCD matter from RHIC to the LHC and into the future collider era.