Collective-Flow Nuclear Shape Imaging

• The paper introduces an imaging-by-smashing strategy that maps initial nuclear deformation into final-state flow observables, providing a novel method to extract parameters such as β₂, γ, β₃, and β₄.
• Collective-flow-assisted nuclear shape imaging leverages multiparticle cumulant analysis and hydrodynamic simulations to link impact geometry with measurable anisotropies like v₂ and v₃.
• The method integrates high-energy QCD techniques with low-energy nuclear structure data, carefully managing model uncertainties and correlations to refine deformation inferences.

Searching arXiv for the core paper and closely related work to support the article. Searching for the 2025 nuclear-structure critique and the 2024/2025 STAR imaging papers. Collective-flow-assisted nuclear shape imaging is a program in which ultrarelativistic nuclear collisions are used to constrain nuclear deformation parameters by exploiting the mapping from initial-state geometry to final-state collective flow. In this framework, deformation of the colliding nuclei modifies the initial transverse energy density of the quark–gluon plasma, and the subsequent hydrodynamic expansion converts these spatial anisotropies into measurable momentum anisotropies $v_n$, mean-transverse-momentum fluctuations, and mixed cumulants. The method has been formulated as an “imaging-by-smashing” strategy for extracting quantities such as the quadrupole deformation $\beta_2$, triaxiality $\gamma$, and, in extensions, octupole $\beta_3$ and hexadecapole $\beta_4$ contributions (Collaboration, 2024). Its status is, however, intrinsically dual: collider flow is a well-founded probe of initial geometry, yet claims of direct laboratory-frame “precision imaging” of even-even $J^\pi=0^+$ ground states have been criticized as conceptually overstated, with robust inference requiring explicit attention to higher-body correlations and to pre-existing low-energy nuclear-structure constraints (Dobaczewski et al., 7 Jul 2025).

1. Historical emergence and domain of application

The method emerged from a convergence of relativistic heavy-ion phenomenology and nuclear-structure modeling. Event-shape engineering established that event-by-event flow fluctuations can be used to select collision ensembles with targeted initial anisotropies through the reduced flow vector $q_n$, thereby providing a practical route to geometry-conditioned measurements (Schukraft et al., 2012). In parallel, studies of small systems argued that collective flow in $p+A$, $d+A$, and light-ion collisions retains sensitivity to projectile geometry, including deuteron deformation, three-body triangularity, and $\alpha$-cluster substructure (Bozek et al., 2014).

A more explicit precursor of nuclear shape imaging was the proposal to use ultrarelativistic $\beta_2$0-$\beta_2$1 collisions to detect triangular $\beta_2$2-cluster correlations through the transmutation of initial triangularity into $\beta_2$3 (Broniowski et al., 2014). Subsequent AMPT studies of $\beta_2$4 argued that the ratio $\beta_2$5 can distinguish chain, triangular, and Woods–Saxon-like $\beta_2$6 configurations, especially in the most-central events (Zhang et al., 2017). At lower energies, EQMD simulations of $\beta_2$7 around $\beta_2$8–$\beta_2$9 MeV/nucleon found that directed flow $\gamma$0 of free protons is more sensitive than elliptic flow to four-$\gamma$1 cluster geometries, establishing an analogous shape-imaging logic outside the QGP regime (Guo et al., 2018).

The modern formulation was articulated for ultrarelativistic heavy-ion collisions at RHIC and the LHC, with benchmark applications to $\gamma$2 versus $\gamma$3, to the Ru/Zr isobars, to Xe+Xe, and to proposed light-ion systems such as Ne on Pb in fixed-target configurations (Zhang, 2024). In this expanded usage, the targeted shape descriptors are primarily $\gamma$4 and $\gamma$5, with increasing attention to $\gamma$6, $\gamma$7, neutron-skin effects, and cluster-driven non-ellipsoidal structure (Dobaczewski et al., 7 Jul 2025).

2. Geometric and hydrodynamic framework

Most implementations start from a deformed Woods–Saxon or Fermi density in which the nuclear surface is expanded in spherical harmonics. In its axial form, the radius is written

$\gamma$8

with density

$\gamma$9

For quadrupole triaxiality one uses

$\beta_3$0

A common low-energy comparison is the mapping

$\beta_3$1

and typical magnitudes quoted for well-deformed nuclei are $\beta_3$2–$\beta_3$3, $\beta_3$4 up to $\beta_3$5 in pear-shaped regions, and non-negligible $\beta_3$6 in actinides (Dobaczewski et al., 7 Jul 2025).

Event by event, initial anisotropies are quantified by participant eccentricities,

$\beta_3$7

with density-weighted averages in the transverse plane. Static quadrupole structure contributes primarily to $\beta_3$8, whereas $\beta_3$9 is often fluctuation-dominated at high multiplicity. Hydrodynamics then supplies the leading-order response

$\beta_4$0

with $\beta_4$1 determined by system size, multiplicity, viscosities, and freeze-out conditions. Higher harmonics receive nonlinear contributions such as

$\beta_4$2

This structure makes quadrupole imaging comparatively direct through $\beta_4$3, whereas triaxiality, octupole, and hexadecapole inference is entangled with nonlinear mode coupling and model dependence (Dobaczewski et al., 7 Jul 2025).

In the ultra-central uranium program, STAR further emphasized mixed size–shape observables. In that context, hydrodynamic simulations yielded the parametric forms

$\beta_4$4

which underlie the extraction of both $\beta_4$5 and $\beta_4$6 from central $\beta_4$7 data (Zhang, 2024).

3. Observables and inference strategy

The basic observables are integrated and differential flow harmonics $\beta_4$8, their multiparticle cumulants $\beta_4$9, $J^\pi=0^+$0, $J^\pi=0^+$1, symmetric cumulants

$J^\pi=0^+$2

event-plane correlations, and event-by-event mean-transverse-momentum fluctuations. In the STAR implementation, the shape-sensitive set for uranium imaging included $J^\pi=0^+$3, $J^\pi=0^+$4, and $J^\pi=0^+$5, while the isobar program used ratios of $J^\pi=0^+$6 and $J^\pi=0^+$7 versus multiplicity to constrain $J^\pi=0^+$8, $J^\pi=0^+$9, $q_n$0, and $q_n$1 (Zhang, 2024).

A central methodological device is the use of matched-system ratios, such as U/Au or Ru/Zr, to suppress common QGP-response uncertainties. In the STAR isobar analysis,

$q_n$2

while in the uranium case the ratio of U+U to Au+Au observables isolates deformation-driven changes in central collisions (Zhang, 2024). A later STAR analysis extended the same logic to $q_n$3, $q_n$4, and their covariances with $q_n$5, arguing that $q_n$6-based ratios provide sensitivity to a possible octupole component of $q_n$7 (Collaboration, 21 Jun 2025).

Nonflow suppression is handled primarily through multiparticle cumulants, $q_n$8 gaps, and subevent methods. Event-shape engineering provides a second layer of selection: by binning on the reduced flow vector

$q_n$9

one enriches classes with larger or smaller initial anisotropy, thereby increasing the resolving power of geometry-conditioned measurements (Schukraft et al., 2012). The critique literature argues that, for credible deformation inference, this observable set should be embedded in a global, multi-observable, multi-model fit with priors anchored to low-energy data, uncertainty propagation across initial conditions, hydrodynamics, and hadronic afterburners, and explicit reporting of model discrepancy (Dobaczewski et al., 7 Jul 2025).

4. Benchmark systems and reported results

The method has been applied or proposed across a broad set of systems. The table summarizes the principal benchmarks.

System or comparison	Main observables	Targeted structure
$p+A$0 vs $p+A$1	$p+A$2, $p+A$3, $p+A$4	$p+A$5, $p+A$6, with extensions to $p+A$7, $p+A$8
$p+A$9 vs $d+A$0	$d+A$1, $d+A$2, multiplicity dependence, ratios	$d+A$3, $d+A$4, neutron skin and radius differences
Xe+Xe	flow systematics	triaxiality claims
$d+A$5+heavy target	$d+A$6, multiplicity trends	triangular $d+A$7-cluster correlations
Pb+Ne vs Pb+Ar or Pb+O	$d+A$8, HBT anisotropies	light-nucleus deformation and clustering
$d+A$9 vs $\alpha$0	$\alpha$1, $\alpha$2, $\alpha$3 cumulants	neutron-skin thickness

For $\alpha$4, STAR’s proceedings reported, from IP-Glasma+MUSIC+UrQMD, $\alpha$5 and $\alpha$6, and after combination with Trajectum, $\alpha$7 and $\alpha$8 (Zhang, 2024). A later STAR analysis based on anisotropic and radial flow reported $\alpha$9 and $\beta_2$00 from the joint use of $\beta_2$01 and $\beta_2$02, while $\beta_2$03 alone yielded a lower bound $\beta_2$04 (Collaboration, 21 Jun 2025). Related RHIC proceedings quoted a combined result $\beta_2$05, $\beta_2$06, emphasizing a large deformation with a small but nonzero departure from axial symmetry (Zhang, 3 Oct 2025).

In the isobar system, STAR reported $\beta_2$07, $\beta_2$08, and neutron-distribution differences $\beta_2$09 fm and $\beta_2$10 fm for Ru–Zr, inferred from $\beta_2$11 and $\beta_2$12 ratio systematics (Zhang, 2024). For octupole deformation of uranium, hydrodynamic calculations predicted that in $\beta_2$13–$\beta_2$14 ultra-central events the ratio $\beta_2$15 should follow a linear dependence on $\beta_2$16, with $\beta_2$17 to $\beta_2$18 for $\beta_2$19–$\beta_2$20, while $\beta_2$21 should show a characteristic suppression (Zhang et al., 21 Apr 2025).

The method has also migrated to fixed-target LHC kinematics. In PbNe and PbAr collisions at $\beta_2$22 GeV, LHCb measured a significantly larger elliptic flow in central PbNe than in PbAr, with

$\beta_2$23

in the $\beta_2$24–$\beta_2$25 centrality bin, interpreting the result as qualitative confirmation of the distinctive bowling-pin shape of $\beta_2$26 (collaboration, 15 Sep 2025). A complementary AMPT study of Pb+Ne and Pb+O found that the freeze-out eccentricity

$\beta_2$27

is significantly larger for NLEFT-like $\beta_2$28 than for Woods–Saxon neon or either oxygen configuration, indicating that azimuthally sensitive femtoscopy can function as a geometric discriminator alongside flow (Kincses, 6 Jun 2025).

For clustered light nuclei, GLISSANDO and AMPT studies of $\beta_2$29 with heavy targets showed that triangular configurations produce increasing $\beta_2$30 and hence increasing $\beta_2$31 with multiplicity, while chain-like or smooth configurations favor larger $\beta_2$32; the discriminating observable is the ratio $\beta_2$33, which satisfies $\beta_2$34 in central $\beta_2$35 collisions (Zhang et al., 2017).

5. Nuclear-structure interpretation and conceptual controversy

The central conceptual dispute concerns the meaning of “shape” in an even-even $\beta_2$36 ground state. The critique literature stresses that such a state is rotationally invariant in the laboratory frame and therefore spherical in that frame; it does not possess a fixed orientation or an “instantaneous shape” that can be directly photographed before impact (Dobaczewski et al., 7 Jul 2025). On this view, deformed Woods–Saxon sampling of one-body densities is a practical phenomenological ansatz, not a literal representation of the laboratory-frame many-body wave function.

The same critique sharpens the information-theoretic content of the inference problem. According to that analysis, axial shape information in a $\beta_2$37 state resides in two-body densities, or conditional probabilities for finding one nucleon given another, whereas triaxiality requires three-body densities. Standard initial-condition generators, which sample nucleons independently from one-body deformed densities, do not encode this correlation structure. This is the basis for the claim that flow–geometry phenomenology is sound, but laboratory-frame “precision imaging” language is not (Dobaczewski et al., 7 Jul 2025).

The practical consequence is not that collider data are irrelevant, but that they must be interpreted as model-dependent constraints cross-validated against low-energy structure. Decades of spectroscopy, static quadrupole moments in odd-$\beta_2$38 nuclei, $\beta_2$39 values, charge radii, and electron scattering already define deformation systematics with high precision; collider analyses should therefore incorporate these as priors and benchmarks rather than attempt unconstrained extraction in isolation (Dobaczewski et al., 7 Jul 2025).

A later microscopic development moved the discussion closer to the critique’s requirements by formulating shape information directly in terms of transverse two-body densities. For the reduced azimuthal correlator

$\beta_2$40

the harmonic coefficients $\beta_2$41 furnish a pair-density fingerprint of intrinsic collectivity, and the associated moments define deformation measures $\beta_2$42 and $\beta_2$43 in the small-deformation limit (Blaizot et al., 21 Dec 2025). This suggests a microscopic route for reconciling collider observables with the requirement that axial and octupole information be formulated through correlated many-body densities rather than through one-body snapshots alone.

6. Systematics, best practice, and current directions

The method is highly sensitive to systematic uncertainty in both initial conditions and medium response. The main confounders include the choice of initial-condition model, such as MC-Glauber versus IP-Glasma; nucleon–nucleon correlations beyond independent sampling; subnucleonic fluctuations and finite nucleon size; neutron–proton density differences and surface diffuseness; multiplicity-based centrality calibration; residual nonflow from jets and resonance decays; transport coefficients $\beta_2$44 and $\beta_2$45; pre-equilibrium dynamics; hadronic afterburners; and the random orientation of the colliding nuclei (Dobaczewski et al., 7 Jul 2025). A separate hydrodynamic study of Au+Au, Cu+Au, and O+O emphasized that while $\beta_2$46 tracks $\beta_2$47 in large systems over suitable centralities, hadronization, nonlinear response, and hadronic rescattering become much more intrusive in small systems (Chen et al., 2024).

Best practice therefore consists of anchoring nuclear inputs to validated low-energy structure data, performing multi-observable global fits across centrality and collision species, scanning multiple initial-condition and hydrodynamic frameworks, and reporting credible intervals that include model discrepancy (Dobaczewski et al., 7 Jul 2025). Multi-particle cumulants and matched-system ratios are preferred over absolute two-particle observables, and consistent centrality calibration across species is essential. For claims involving triaxiality, the critique literature further recommends explicit tests of sensitivity to three-body correlations, since current one-body initialization practices do not encode them (Dobaczewski et al., 7 Jul 2025).

Several current directions extend the method beyond quadrupole imaging. Quantum constraints on collective fluctuations have been incorporated in ultra-central Pb+Pb through an energy-weighted sum-rule restriction on quadrupole variance, which reduces $\beta_2$48 and resolves the long-standing $\beta_2$49-to-$\beta_2$50 puzzle in central Pb collisions (Zakharov, 2020). Femtoscopy has been proposed as an additional imaging channel, with azimuthally sensitive HBT radii supplying direct information on freeze-out source eccentricity (Kincses, 6 Jun 2025). Calcium-isotope studies within AMPT indicate that neutron-skin thickness in $\beta_2$51 affects $\beta_2$52, $\beta_2$53, and especially the variance and skewness of $\beta_2$54, suggesting an “imaging-by-smashing” route to $\beta_2$55 (Vitsos et al., 27 Nov 2025). At the microscopic end, harmonic analysis of two-body correlations in ab initio calculations of $\beta_2$56 and $\beta_2$57 proposes a bridge from collider flow observables to the intrinsic many-body correlation structure itself (Blaizot et al., 21 Dec 2025).

Collective-flow-assisted nuclear shape imaging is thus best understood not as a literal camera for laboratory-frame ground-state shapes, but as a constrained inverse problem linking nuclear densities, many-body correlations, initial-state geometry, and final-state collective observables. In that more precise sense, it has become a distinct interface between nuclear structure and high-energy QCD phenomenology: strongest when used comparatively, statistically, and in conjunction with low-energy information, and most credible when its limitations are made as explicit as its geometric sensitivity.