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 $\sim 50$ MeV/$c$ (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 $\frac{dN_\mathrm{ch}}{d\eta}$ 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,

$$\left.\frac{dN_\mathrm{ch}}{d\eta}\right|_{\eta=0} \propto s^\beta,$$

where $\beta \sim 0.10$–$0.12$ in pp and increases to $\sim 0.15$–$0.16$ 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, $\frac{dN_\mathrm{ch}}{d\eta}$ increases with centrality (smaller impact parameter, larger overlap) and scales approximately with the number of participant nucleons $N_\mathrm{part}$, especially at low/intermediate energies and in peripheral collisions. The normalized quantity,

$$\frac{dN_\mathrm{ch}/d\eta}{\langle N_\mathrm{part}/2 \rangle},$$

often exhibits a weak centrality dependence from central to semi-peripheral collisions and only a modest decrease (factor $\sim1.7$–$1.8$) 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 $\sqrt{s_{NN}}$ using one detector (ALICE) reveal a smooth increase of midrapidity $\left.\frac{dN_\mathrm{ch}}{d\eta}\right|_{\eta=0}$ with system size. Central Pb–Pb yields reach $\sim 2000$ per unit $\eta$ at 5.02–5.36 TeV; in central p–Pb, $\sim 60$; and in pp, $\sim 5.7$ (Collaboration, 2022).

Tabulated energy and system-size scaling:

System	Energy (TeV)	$\left.\frac{dN_\mathrm{ch}}{d\eta}\right|_{\eta=0}$ (midrapidity)
pp (INEL>0)	13.6	$7.12^{+0.12}_{-0.09}$
Pb–Pb (0–5%)	5.36	$2004 \pm 52$
p–Pb (NSD)	8.16	$19.1 \pm 0.7$

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 $\frac{dN_\mathrm{ch}}{d\eta}$ in A–A and smaller systems are driven more by the geometry of the collision zone (quantified via $N_\mathrm{part}/2A$, where $A$ is the nuclear mass number) than by $N_\mathrm{part}$ alone. Collisions with the same participant fraction $N_\mathrm{part}/2A$ from different nuclei yield more similar pseudorapidity distributions over the full range of $\eta$ (0806.2803).
• Longitudinal scaling: Plotting $\frac{dN_\mathrm{ch}}{d\eta}$ as a function of $\eta' = \eta - y_\mathrm{beam}$ reveals longitudinal (“limiting fragmentation”) scaling at forward/backward rapidities across a broad energy range, consistent with an invariance of the fragmentation region as $\sqrt{s_{NN}}$ 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-$x$ 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, $N_\mathrm{coll}$-dependent) components. The centrality and energy trends favor a dominantly participant-driven (“soft”) scaling, especially at RHIC and lower LHC energies ($\frac{dN_\mathrm{ch}}{d\eta}$ 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 $N_\mathrm{part}^\alpha s_{NN}^\gamma$ where $\gamma$ 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, $\sigma$, extracted from a Gaussian parameterization of $dN_\mathrm{ch}/dy$, increases with energy and system size, exceeding expectations from Landau hydrodynamics. Its proportionality to the beam rapidity $y_\mathrm{beam}$ 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,

$$\epsilon_\mathrm{Bj} = \frac{1}{c\tau A}\,\langle dE_T/dy \rangle,$$

where $A$ is the overlap area and $\tau$ the formation time, $dN_\mathrm{ch}/d\eta$ constrains $\langle dE_T/dy \rangle$ via measured charged-particle multiplicities and mean $p_T$. The rise in $\frac{dN_\mathrm{ch}}{d\eta}$ with system size and centrality reflects a tenfold boost in initial energy density from pp to Pb–Pb at fixed $\sqrt{s_{NN}}$, 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 $\frac{dN_\mathrm{ch}}{d\eta}$ 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 $\frac{dN_\mathrm{ch}}{d\eta}$ 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}$ per participant pair as $s_{NN}^{0.153}$, anticipated densities at midrapidity reach $\sim 3600$ 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 $\eta$ 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 $\eta$	$-\ln \tan(\theta/2)$	Universally adopted for longitudinal kinematics
Density $\frac{dN_\mathrm{ch}}{d\eta}$	Primary charged particles per unit $\eta$	Central to all multiplicity analyses
Normalized density	$$\frac{dN_\mathrm{ch}/d\eta}{\langle N_\mathrm{part}/2 \rangle}$$	System size, centrality scaling
Longitudinal shift	$\eta^\prime = \eta - y_\mathrm{beam}$	For limiting fragmentation studies
Energy scaling	$$\left.\frac{dN_\mathrm{ch}}{d\eta}\right|_{\eta=0} \propto s^\beta$$	$\beta\sim0.115$ (pp), $0.156$ (A–A, central)
Energy density	$\epsilon_\mathrm{Bj}$ via $dN_\mathrm{ch}/d\eta$	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.