⚛️ nuclear experiments

Fundamental geometric limitations on disentangling nuclear-surface properties in relativistic heavy ion collisions

This paper demonstrates that the extraction of nuclear surface diffuseness from relativistic heavy ion collisions is fundamentally limited by strong geometric degeneracies between the surface parameter and intrinsic nuclear deformations, though multiparticle triangular correlations offer a more independent pathway to constrain these properties.

Original authors: Hadi Mehrabpour, Behnaz Behzadmoghaddam, S. M. A. Tabatabaee Mehr, Oscar Garcia-Montero, Li Yan

Published 2026-09-02

📖 5 min read🧠 Deep dive

CC0 1.0

Original authors: Hadi Mehrabpour, Behnaz Behzadmoghaddam, S. M. A. Tabatabaee Mehr, Oscar Garcia-Montero, Li Yan

Original paper dedicated to the public domain under CC0 1.0 (http://creativecommons.org/publicdomain/zero/1.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Atoms are the fundamental building blocks of matter, but they are not solid, uniform spheres. Instead, they possess a dense core surrounded by a fuzzy cloud of particles, much like a peach with a pit and a soft, gradually thinning flesh. Physicists call this gradual thinning the "surface diffuseness," a property that describes how sharply the edge of the atom's core is defined. While scientists have long known the size of the core and the arrangement of protons inside, the exact nature of this fuzzy edge, particularly for neutrons, remains one of the most elusive puzzles in nuclear physics. Because neutrons carry no electric charge, they cannot be probed directly with the same tools used to map protons. To understand them, researchers must rely on indirect methods, often smashing heavy atoms together at nearly the speed of light and analyzing the spray of particles that results.

For decades, scientists have hoped that these high-speed collisions could reveal the hidden shape of the neutron cloud. The idea is that the initial geometry of the colliding atoms imprints itself on the final pattern of particles, allowing physicists to work backward and deduce the original structure. However, a new study suggests that this task is far more complicated than previously thought. The researchers found that the signals used to measure the fuzziness of the nuclear surface are deeply entangled with another property: the overall shape of the nucleus itself. Just as a slightly squashed ball and a ball with a fuzzier edge can produce similar ripples when dropped in water, the intrinsic deformation of an atom can mask the very surface details scientists are trying to measure.

The team, led by researchers from institutions in China, Iran, and Spain, set out to untangle this confusion using sophisticated computer simulations. They focused on two specific types of heavy atoms: uranium-238 and neon-20. Uranium is naturally elongated, resembling a rugby ball, while neon is more spherical but can still be distorted. The researchers simulated millions of collisions, systematically changing two key variables in their models: the surface diffuseness, which controls how soft the nuclear edge is, and the deformation parameters, which control how stretched or pear-shaped the nucleus is. They then watched how these changes affected the patterns of particles emerging from the collisions, specifically looking at how groups of two or three particles correlated with one another.

The results revealed a significant roadblock. The researchers discovered that the most common signals used to study nuclear collisions are overwhelmingly dominated by the overall shape of the nucleus. When they tried to measure the surface fuzziness using these standard signals, the changes in the nuclear shape acted like a loud noise drowning out a whisper. In the most central collisions, where the atoms hit head-on, the signal for surface fuzziness was so weak compared to the signal for nuclear shape that it was nearly impossible to separate the two. The simulations showed that a nucleus with a sharp edge but a specific shape could produce the exact same particle pattern as a nucleus with a fuzzy edge but a slightly different shape. This creates a "degeneracy," a situation where two different physical realities look identical to the observer.

However, the study did not end in a dead end. The researchers found that not all signals are equally confused. While the standard two-particle correlations were heavily biased toward the nuclear shape, other signals involving three particles or specific triangular patterns offered a different perspective. These less common correlations responded more independently to the surface fuzziness, retaining some of the information that the other signals lost. By combining these different types of measurements, the team showed that it is possible to gain a clearer picture, though the path remains difficult. They used advanced statistical tools to map out exactly how much information could be recovered, finding that even with the best combination of signals, the uncertainty regarding the surface fuzziness remains substantial.

The study also tested whether these limitations were specific to the type of atom or the computer model used. They compared their results for uranium collisions at one energy level with neon collisions at a much higher energy, and they used two different theoretical frameworks to describe the initial collision. In every case, the fundamental problem persisted: the intrinsic shape of the nucleus consistently obscured the surface details. The specific numbers changed depending on the model, but the core finding remained robust. The inability to cleanly separate the surface fuzziness from the nuclear shape is not an error in the simulation or a quirk of a specific element; it is a fundamental geometric limitation inherent to the way these collisions work.

Ultimately, this research clarifies the boundaries of what can be learned from heavy ion collisions. It demonstrates that while these high-energy experiments are powerful tools for probing nuclear structure, they cannot simply extract the surface diffuseness in isolation. The shape of the nucleus and the fuzziness of its edge are inextricably linked in the data. To truly understand the neutron skin of an atom, scientists must treat the nuclear shape and the surface properties as a single, interconnected system. The path forward requires using a diverse set of particle correlations and combining high-energy collision data with other, lower-energy measurements to break the geometric deadlock. The study does not close the door on this field, but it provides a crucial map of the terrain, showing exactly where the shadows lie and where the light of new information can still penetrate.

Technical Summary: Fundamental Geometric Limitations on Disentangling Nuclear-Surface Properties in Relativistic Heavy Ion Collisions

Problem Statement
The extraction of nuclear surface diffuseness ( $a_0$ ) from relativistic heavy ion collisions is fundamentally challenged by the intertwined responses of geometric observables to both surface diffuseness and intrinsic nuclear deformation (specifically quadrupole $\beta_2$ and octupole $\beta_3$ deformations). While ultrarelativistic collisions convert initial geometric fluctuations into collective final-state particle correlations via hydrodynamic response, these observables often lack specificity. Intrinsic deformation can dominate the same geometric observables that are sensitive to the nuclear surface, creating a "deformation–diffuseness degeneracy." This makes it difficult to determine whether a measured correlation signal arises from the nuclear surface profile or from the intrinsic shape of the nucleus. The central question addressed is whether multiparticle correlations can provide independent information on $a_0$ once intrinsic deformation is allowed to vary, or if the response remains obscured by the deformation parameters.

Methodology
The authors employ event-by-event Monte Carlo Glauber simulations to investigate this limitation, focusing on $^{238}\text{U} + ^{238}\text{U}$ and $^{20}\text{Ne} + ^{20}\text{Ne}$ collisions.

Initial State Modeling: The nuclear density is described using a deformed Woods–Saxon (WS) parameterization, where the half-density radius ( $R_0$ ) is fixed to experimental values, while the surface diffuseness ( $a_0$ ), quadrupole deformation ( $\beta_2$ ), and octupole deformation ( $\beta_3$ ) are systematically varied. The study utilizes the TrENTo framework for initial-state generation, with selected comparisons to the CGC-based McDipper model to assess model dependence regarding nucleon width and short-range correlations.
Observables: The analysis focuses on two-particle correlation measures ( $\varepsilon_2\{2\}$ , $\varepsilon_3\{2\}$ , and transverse-size correlation $c_d\{2\}$ ) and three-particle flow–transverse-size correlations ( $\text{cov}(2)$ and $\text{cov}(3)$ ). These are linked to final-state observables through linear hydrodynamic response.
Sensitivity and Information Analysis:
- Local Sensitivity: The authors define sensitivity vectors ( $S_X$ ) to quantify the local dependence of observables on scaled WS parameters. The angle between these vectors determines the complementarity of different observables.
- Degeneracy Directions: Local constant-observable contours are calculated to identify directions in parameter space where changes in $a_0$ can be compensated by changes in deformation parameters.
- Information Geometry: A Fisher-information matrix is constructed to quantify the combined information content and identify the least-constrained directions (weak eigenvectors) in the parameter space.
- Bayesian Analysis: A Bayesian response-model analysis is performed to evaluate posterior distributions for the nuclear parameters, incorporating both Monte Carlo statistical uncertainties and model-discrepancy contributions.

Key Results

Dominance of Quadrupole Deformation: Observables dominated by elliptic geometry, particularly $\varepsilon_2\{2\}$ and $c_d\{2\}$ , exhibit a strong response to quadrupole deformation ( $\beta_2$ ). This leads to a local degeneracy where variations in $a_0$ are largely compensated by changes in $\beta_2$ , substantially limiting the independent sensitivity of these observables to surface diffuseness.
Complementary Role of Triangular Correlations: Triangular correlations ( $\varepsilon_3\{2\}$ ) and three-particle covariance observables ( $\text{cov}(3)$ ) provide a more distinct response to the nuclear surface. Their sensitivity directions differ significantly from those of elliptic correlations, offering additional, albeit not complete, information on $a_0$ .
Local Nature of Degeneracy: The deformation–diffuseness degeneracy is not a global property of the parameter space but is local and configuration-dependent. The direction of the least-constrained parameter combination changes as the underlying nuclear configuration (values of $\beta_2$ and $\beta_3$ ) varies. For instance, in the presence of finite octupole deformation, the response of $\varepsilon_3\{2\}$ rotates to become predominantly sensitive to $\beta_3$ , while $\varepsilon_2\{2\}$ acquires a stronger $\beta_3$ component as $\beta_2$ decreases.
Persistence of Weak Constraints: Despite combining multiple observables, the joint Fisher-information contours remain strongly elongated along the $a_0$ direction. The least-constrained direction is nearly aligned with the surface-diffuseness axis across various initial-state model configurations (TrENTo and McDipper) and collision systems (U+U and Ne+Ne). This indicates that increasing the number of correlations does not automatically eliminate the ambiguity.
Robustness: The systematic dependence of correlations on nuclear deformation is robust against changes in collision energy, nuclear species, and initial-state prescriptions (nucleon width and short-range correlations). However, the quantitative magnitude of the response varies with the model, necessitating the inclusion of model uncertainties in any extraction.

Significance and Claims
The paper establishes that the extraction of nuclear surface diffuseness from relativistic heavy ion collisions is fundamentally limited by intrinsic multipole deformation. The authors claim that:

The limitation is not merely a statistical issue or a consequence of a specific model choice but arises from the intrinsic structure of the observable response to nuclear geometry.
No single multiparticle correlation can uniquely isolate $a_0$ ; instead, a simultaneous analysis of complementary correlations (elliptic and triangular) is required to partially break the degeneracy.
The identifiability of the Woods–Saxon parameters ( $a_0, \beta_2, \beta_3$ ) is intrinsically configuration-dependent. There is no universal hierarchy where a specific correlation maps uniquely to a single nuclear-structure parameter.
Reliable extraction of surface diffuseness requires the simultaneous treatment of intrinsic deformation, the use of complementary multiparticle correlations, and a rigorous accounting of initial-state model uncertainties. The authors suggest that combining these high-energy observables with independent low-energy nuclear-structure constraints will be essential for disentangling the surface profile from intrinsic deformation.

More like this