⚛️ nuclear theory

A novel method to determine the 4th moment of the 208Pb charge density from low-momentum-transfer electron scattering

This paper proposes a novel method to precisely determine the fourth moment of the 208Pb nuclear charge-density distribution using low-momentum-transfer elastic electron scattering.

Original authors: Rika Danjo, Toshimi Suda, Toshio Suzuki

Published 2026-09-22

📖 4 min read🧠 Deep dive

CC0 1.0

Original authors: Rika Danjo, Toshimi Suda, Toshio Suzuki

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 the visible world, but inside them lies a complex, invisible architecture. At the heart of every atom sits the nucleus, a dense cluster of protons and neutrons. While protons carry a positive electric charge, neutrons are neutral, yet both contribute to the nucleus's overall size and shape. For decades, physicists have mapped the distribution of electric charge within these nuclei by firing streams of electrons at them. When electrons bounce off a nucleus, their scattering patterns reveal how the charge is arranged, much like how a sonar ping reveals the shape of a submerged object. This technique has allowed scientists to measure the average size of the nucleus with great precision. However, knowing the average size is only part of the story. To fully understand the internal structure, particularly how the invisible neutrons are distributed relative to the visible protons, scientists need to measure higher-order details of the charge distribution. These details act as a sensitive probe for the neutron cloud, which is crucial for understanding the forces that hold the nucleus together and the behavior of exotic, unstable atoms found in the cosmos.

A team of researchers at Tohoku University in Japan has now proposed a refined method to extract one of these elusive details, known as the fourth moment of the charge density, from the nucleus of lead-208. Lead-208 is a heavy, stable isotope that serves as a standard reference in nuclear physics. The researchers' goal is to determine the shape of the charge distribution with such precision that it reveals the subtle influence of the neutrons. Previously, obtaining this information required scattering electrons at a wide range of speeds, including very high speeds that are difficult to achieve with unstable, short-lived nuclei. The new approach, however, relies entirely on data from low-speed electron scattering. By focusing on this low-energy range, the method avoids the technical hurdles associated with high-speed collisions, opening the door to studying nuclei that cannot be held in a target for long.

The core challenge in this work is that heavy nuclei like lead-208 possess a massive electric charge, which significantly distorts the path of the incoming electrons. This distortion, caused by the Coulomb force, complicates the simple mathematical models usually used to interpret scattering data. In the past, scientists often had to rely on approximations that were not precise enough for the level of detail required here. The authors of this paper developed a new strategy to handle this distortion. Instead of using a rough approximation, they calculated the exact effect of the electric field on the electron waves using numerical methods. They then incorporated these precise calculations as a correction factor into their analysis. This allows them to strip away the distortion and recover the true underlying pattern of the charge distribution from the experimental data, even for such a heavy nucleus.

To test their method, the researchers used existing theoretical models of the lead-208 nucleus to simulate what the experimental data would look like. They applied their new correction technique to these simulated results and found that they could extract the fourth moment of the charge distribution with an uncertainty of about 0.8 percent. This is a significant improvement over the current best measurements, which carry an uncertainty of roughly 1.5 percent. The improvement is not merely a matter of better numbers; it fundamentally changes what can be learned. Because the fourth moment is sensitive to the distribution of neutrons, reducing the uncertainty in this measurement allows physicists to pin down the size of the neutron cloud with much greater accuracy. The study shows that by including specific higher-order terms in their mathematical expansion and focusing on a narrow range of low momentum transfers, the method becomes stable and reliable.

The researchers emphasize that this approach is particularly valuable because it does not require data from the high-speed, high-momentum range where the scattering signal becomes extremely weak. Since the scattering probability drops sharply as the speed of the electrons increases, measuring at high speeds requires an immense number of collisions, which is often impossible with rare, unstable isotopes. By proving that the necessary information can be extracted from the low-speed range alone, this method offers a practical path forward for studying neutron-rich nuclei that are difficult to produce. The team is currently conducting experiments at the Ultra-Low momentum transfer facility at Tohoku University to verify these findings with real data. If successful, this technique will provide a clearer window into the neutron skin of atomic nuclei, helping to refine our understanding of the fundamental forces that govern the structure of matter.

Technical Summary: A Novel Method to Determine the 4th Moment of the $^{208}$ Pb Charge Density

Problem Statement
The charge density distribution of nuclei, $\rho_c(r)$ , and its moments, $\langle r^n \rangle_c$ , are fundamental to understanding nuclear structure. While the mean-square radius ( $\langle r^2 \rangle_c$ ) is well-determined, the 4th moment ( $\langle r^4 \rangle_c$ ) is critical because it provides sensitivity to the point-neutron density distribution ( $\langle r^2 \rangle_n$ ). Current determinations of $\langle r^4 \rangle_c$ for stable nuclei like $^{208}$ Pb rely on extracting $\rho_c(r)$ from elastic electron scattering data over a wide momentum transfer ( $q$ ) range using "model-independent" methods (e.g., Fourier-Bessel or Sum of Gaussians). However, this approach faces two primary limitations:

Precision: The current uncertainty for $^{208}$ Pb is approximately 1.5%, which, while sensitive to the small neutron contribution (estimated at 2–3%), limits the precision of extracted neutron radii.
Experimental Constraints: Extracting $\rho_c(r)$ requires data at high $q$ , where the scattering cross section drops as $1/q^4$ . This necessitates high luminosity, making such measurements difficult or impossible for short-lived, unstable nuclei.

Methodology
The authors propose a new method to determine $\langle r^4 \rangle_c$ using only low-momentum-transfer data ( $q < 0.315 \text{ fm}^{-1}$ ), avoiding the need for high- $q$ measurements. The approach involves three key technical components:

Taylor Expansion in PWBA: Within the Plane Wave Born Approximation (PWBA), the squared charge form factor $|F_{PWBA}(q)|^2$ is expanded as a Taylor series in powers of $q^2$ . The coefficients of this series depend on the charge moments $\langle r^n \rangle_c$ . By truncating this series, $\langle r^4 \rangle_c$ can be isolated and extracted from low- $q$ data.
Convergence Criterion: To ensure the validity of the truncated series, the authors establish a convergence condition based on the Leibniz criterion for alternating series. For $^{208}$ Pb, this restricts the analysis to $q < 0.315 \text{ fm}^{-1}$ .
Numerical Coulomb Correction: For heavy nuclei like $^{208}$ Pb, Coulomb distortion significantly affects electron scattering, rendering the standard PWBA insufficient. The authors reject the Effective Momentum Approximation (EMA) for this specific low- $q$ regime, as it fails to reproduce the "bump-like" structures observed in Phase-Shift Analysis (PSA). Instead, they numerically calculate a multiplicative correction factor, $f_{\rho_c(r)}(E, \theta)$ , by comparing PSA results with PWBA results. They demonstrate that this factor is insensitive to the specific charge density model (FB vs. SOG) and incident energy within the relevant low- $q$ range.

Key Contributions and Results
The paper validates this method using pseudo-experimental data generated from the Fourier-Bessel (FB) charge density distribution of $^{208}$ Pb. The study systematically evaluates the precision and accuracy of the extraction by varying the truncation order of the Taylor expansion (from $FF_2$ to $FF_6$ ) and the $q$ range.

Truncation Order: The analysis reveals that lower-order truncations (e.g., up to $q^4$ ) result in a strong dependence of the extracted $\langle r^4 \rangle_c$ on $q$ , indicating insufficient accuracy. Convergence is achieved only when terms up to at least order $q^{10}$ (corresponding to $FF_5$ ) are included.
Precision Improvement: By utilizing data near the upper convergence limit ( $q \approx 0.3 \text{ fm}^{-1}$ ) and including terms up to $q^{10}$ , the proposed method reduces the uncertainty in $\langle r^4 \rangle_c$ from the current $\sim 1.5\%$ to approximately 0.8%.
Robustness: The extraction of $\langle r^4 \rangle_c$ is shown to be weakly sensitive to the specific values used for higher-order moments ( $\langle r^6 \rangle_c$ and above), allowing the use of phenomenological models (like the 2-parameter Fermi distribution) for these higher terms without compromising the result.

Significance
The paper claims that this method offers a pathway to determine $\langle r^4 \rangle_c$ with significantly improved precision without requiring high-luminosity facilities or high- $q$ data.

Neutron Radius Constraints: The reduction in uncertainty for $\langle r^4 \rangle_c$ is projected to reduce the uncertainty in the root-mean-square radius of the point-neutron density ( $\sqrt{\langle r^2 \rangle_n}$ ) by approximately 30% compared to previous studies.
Applicability to Unstable Nuclei: Because the method relies on the low- $q$ region where cross sections are larger, it is presented as a promising approach for future studies of unstable, neutron-rich nuclei. The authors note that such measurements are currently in progress at the Ultra-Low momentum transfer $Q^2$ (ULQ2) facility at RARiS, Tohoku University, and suggest potential applicability at the SCRIT facility at RIKEN.

The work does not claim to solve the fundamental theoretical uncertainties in the relation between charge and neutron densities but provides a more precise experimental input to constrain these theoretical frameworks.

More like this