---
title: Charged-Particle Pseudorapidity Densities
url: https://www.emergentmind.com/topics/charged-particle-pseudorapidity-densities
type: topic
---

# Charged-Particle Pseudorapidity Densities

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$ [2409.18560, 1002.0621].
- **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 [2211.07576, 1411.4963].
- **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.) [2504.02505, 2409.18560].

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 [2211.15364, 2409.18560].
- **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 [2504.02505, 1509.07299].
- **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$ [2204.10210].

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 [1304.0347, 2008.07802].
- **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 [2008.07802].

## 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 [2008.07802].
- **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 [2008.07802, 1210.3615, 1812.01312].
- **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 [2504.02505, 2409.18560].
- **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 [1612.08966, 1304.0347].

## 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 [2204.10210].
- **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) [2211.07576].

## 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 [2008.07802, 1304.0347].
- **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) [2008.07802].
- **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 [2506.14989, 2409.18560].

## 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.

Source: https://www.emergentmind.com/topics/charged-particle-pseudorapidity-densities