⚛️ nuclear theory

Physics-informed interpretable forms for dipole resonances and quantities relevant for astrophysics

This paper utilizes Bayesian inference to model isovector giant dipole resonance properties with structural effects, deriving symmetry energy coefficients that correlate with neutron-skin thickness and provide critical inputs for nuclear astrophysics calculations.

Original authors: Gourab Banerjee, Jhilam Sadhukhan, Debasish Mondal, B. K. Agrawal, Chandrani Sen, S. Mukhopadhyay

Published 2026-10-09

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Gourab Banerjee, Jhilam Sadhukhan, Debasish Mondal, B. K. Agrawal, Chandrani Sen, S. Mukhopadhyay

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.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

Inside the heart of every atom lies a dense, churning core where protons and neutrons are packed together in a state of extreme pressure. These particles are not static; they move, vibrate, and interact in complex ways that define the very nature of matter. One of the most fascinating ways these particles behave is through a collective motion known as the isovector giant dipole resonance. Imagine the protons and neutrons inside the nucleus as two distinct fluids. In this specific type of vibration, the fluid of positively charged protons sloshes in one direction while the fluid of neutral neutrons sloshes in the opposite direction, creating a rhythmic back-and-forth oscillation. This motion is not just a curiosity of the atomic world; it acts as a sensitive probe for the forces that hold the nucleus together. By studying how these fluids move, scientists can measure the "symmetry energy," a fundamental property that dictates how the nucleus behaves when it has an imbalance of protons and neutrons. This property is crucial for understanding some of the most extreme environments in the universe, such as the dense interiors of neutron stars and the violent collisions where heavy elements are forged.

For decades, physicists have tried to describe this sloshing motion using mathematical models. Early attempts treated the nucleus like a simple drop of liquid, offering two different ways to calculate the energy of the vibration. One model suggested the energy depended on the size of the nucleus in a specific way, while another suggested a different relationship. Real-world data showed that actual nuclei did not fit perfectly into either of these simple boxes; instead, they displayed a mix of behaviors that depended on whether the nucleus was round or stretched out, and how many more neutrons it had than protons. While powerful computer simulations exist to predict these behaviors, they are often too complex to use for the thousands of calculations required to model the life cycles of stars. Scientists needed a simpler, clearer way to describe these vibrations that still captured the essential physics, one that could be applied broadly without needing a supercomputer for every single calculation.

A team of researchers from India has now developed a new, refined approach to solve this problem. They turned to a statistical method called Bayesian inference, which allows scientists to extract the most likely values for physical properties from noisy experimental data while accounting for uncertainties. The team analyzed a vast collection of experimental data, specifically looking at how different nuclei absorb gamma rays to trigger this proton-neutron sloshing. They examined data from 81 spherical nuclei, which are round like billiard balls, and 60 deformed nuclei, which are stretched like rugby balls. By feeding this data into their statistical framework, they were able to pinpoint the exact energy and width of the resonance for each nucleus with remarkable precision.

The key to their success was building a new formula that explicitly included the structural quirks of the nucleus. Previous models often ignored how the shape of the nucleus or the specific arrangement of its internal particles affected the vibration. The new approach, however, incorporated the effects of nuclear deformation and the "shell structure," which refers to how protons and neutrons fill up energy levels similar to seats in a theater. When the researchers tested their new formulas against the experimental data, they found that the predictions matched the measurements almost perfectly, even for nuclei that were far from the stable, round shapes usually studied. They also discovered that the two different historical models for the vibration energy were not mutually exclusive but were actually connected through a deeper physical mechanism, with the new formula successfully bridging the gap between them.

Beyond just describing the vibration, the researchers used these new findings to uncover deeper truths about the forces inside the nucleus. By analyzing the resonance data, they were able to calculate the coefficients of the nuclear symmetry energy, which describe how the energy of the nucleus changes as the balance between protons and neutrons shifts. Their results aligned closely with recent independent measurements of nuclear dipole polarizability, a related property that measures how easily the nucleus can be distorted by an electric field. This agreement provided a strong check on their work, confirming that their statistical extraction of the data was reliable.

The implications of this work extend far beyond the laboratory. The team used their newly determined symmetry energy values to predict the thickness of the "neutron skin" in heavy nuclei. This is a layer of excess neutrons that forms on the surface of a nucleus when there are many more neutrons than protons. The thickness of this skin is a critical piece of information for astrophysicists because it is directly linked to the size and structure of neutron stars. The researchers found a strong correlation between the dipole polarizability and the neutron-skin thickness, allowing them to estimate the skin's thickness with high confidence. When they applied their method to the lead-208 nucleus, their predictions overlapped significantly with recent measurements from the PREX collaboration, a major international experiment dedicated to measuring this very property.

This research provides a powerful new tool for the scientific community. By offering a compact, physically grounded formula that accurately describes how nuclei vibrate, the team has created a resource that can be easily used in large-scale calculations for nuclear astrophysics. These calculations are essential for modeling how stars explode and how heavy elements are created in the universe. The new formulas allow scientists to estimate the strength of the interaction between light and matter in unstable nuclei that cannot be studied directly in a lab. This capability is vital for understanding the conditions inside neutron stars and the processes that occur during their mergers. The work demonstrates that by combining rigorous statistical analysis with a deep understanding of nuclear structure, scientists can refine our picture of the atomic world and, in doing so, illuminate the most extreme corners of the cosmos.

Technical Summary: Physics-informed Interpretable Forms for Dipole Resonances and Quantities Relevant for Astrophysics

Problem Statement
The isovector giant dipole resonance (IVGDR) is a critical collective excitation mode providing insight into nuclear structure, dynamics, and the nuclear symmetry energy. While microscopic models (e.g., TDHF, RPA) can predict IVGDR centroid energies ( $E_G$ ) with reasonable accuracy, they often struggle to reproduce the resonance width ( $\Gamma_G$ ) and are computationally prohibitive for large-scale applications. Conversely, classical hydrodynamic models offer simple analytical expressions but lack the precision to capture structural effects like isospin asymmetry, nuclear deformation, and shell corrections. Furthermore, existing data-driven machine learning approaches often lack physical interpretability. There is a need for improved analytical formulations that retain clear physical interpretability while accurately reproducing IVGDR systematics across the nuclear chart to support large-scale astrophysical calculations, such as radiative capture processes.

Methodology
The authors employ a Bayesian inference framework to extract and optimize IVGDR parameters from experimental photoabsorption data.

Data Source: The analysis utilizes measured $\gamma$ -ray spectra for 81 spherical and 60 axially deformed nuclei from the EXFOR database.
Modeling: The photoabsorption cross-section is modeled using a Simplified Lorentzian (SLO) distribution. The authors propose new, physically motivated functional forms for the centroid energy ( $E_G$ $E_{G}$ ) and width ( $\Gamma_G$ $Γ_{G}$ ) that explicitly incorporate:
- Isospin Asymmetry: A term dependent on $I = |N-Z|/A$ is added to the centroid energy scaling, extending beyond the traditional Goldhaber-Teller ( $A^{-1/6}$ ) and Steinwedel-Jensen ( $A^{-1/3}$ ) models.
- Nuclear Deformation: For axially deformed nuclei, the centroid energy and width are adjusted based on the quadrupole deformation parameter $\beta_2$ , which is derived from the experimental splitting of the IVGDR strength distribution.
- Shell Corrections: The damping coefficient ( $k_s$ ) in the width formula is refined to include a dependence on the ground-state shell-correction energy ( $\delta W$ ), addressing local fluctuations near magic numbers.
Optimization: Bayesian inference is used to determine the posterior probability distributions of the model parameters ( $a_0, a_1, a_I$ for energy; $k_s, \delta$ for width) and to quantify correlations between them.
Astrophysical Application: The extracted parameters are used to constrain the volume ( $K_v$ ) and surface ( $K_s$ ) symmetry energy coefficients. These are subsequently employed to predict nuclear dipole polarizability ( $\alpha_D$ ) and investigate its correlation with neutron-skin thickness ( $r_{skin}$ ) using various Skyrme energy density functionals (EDFs).

Key Contributions

Interpretable Parametrizations: The paper presents compact, analytical expressions for IVGDR centroid energy and width that successfully reproduce experimental systematics for both spherical and deformed nuclei.
Structural Integration: Unlike previous global fits, these forms explicitly integrate isospin asymmetry, deformation, and shell effects, offering a more robust description across the nuclear chart.
Bayesian Constraints: The study provides rigorous Bayesian constraints on the nuclear symmetry energy coefficients ( $K_v, K_s$ ) derived directly from IVGDR data.
Predictive Tool for Neutron Skins: The authors establish a linear correlation between the product of dipole polarizability and symmetry energy ( $\alpha_D J$ ) and the neutron-skin thickness ( $r_{skin}$ ), providing a practical formula to estimate $r_{skin}$ for heavy nuclei.

Results

Centroid Energy: The proposed expression, $E_G(A, I) = (a_0 A^{-1/3} + a_1 A^{-1/6})(1 - a_I I^2)$ , captures the subtle isospin dependence of centroid energies, showing better agreement with experimental data than traditional forms, particularly for nuclei with $A < 50$ . The analysis reveals a strong anti-correlation ( $r \approx -0.99$ ) between the $A^{-1/3}$ and $A^{-1/6}$ coefficients, suggesting a common underlying physical mechanism.
Width and Shell Effects: The inclusion of shell corrections in the damping coefficient $k_s$ significantly improves the reproduction of experimental widths, particularly near magic numbers where $k_s$ is suppressed. For deformed nuclei, distinct power-law exponents ( $\delta_\parallel \approx 3.66, \delta_\perp \approx -1.06$ ) were found necessary to describe the width anisotropy, outperforming single-exponent models.
Symmetry Energy and Polarizability: The optimized symmetry energy coefficients ( $K_v \approx 36.26$ MeV, $K_s \approx 57.64$ MeV) are in good agreement with recent measurements. The predicted dipole polarizability ( $\alpha_D$ ) matches experimental data well, with the model based on Eq. (18) providing an upper bound and the droplet model (Eq. 20) providing a lower bound.
Neutron-Skin Thickness: A strong linear correlation ( $r > 0.80$ for Eq. 18, $r > 0.90$ for Eq. 20) is observed between $\alpha_D J$ and $r_{skin}$ for heavy nuclei ( $A > 150$ ). Applying this relation to $^{208}\text{Pb}$ yields a neutron-skin thickness that overlaps with measurements from the PREX collaboration and electric dipole response studies.

Significance and Claims
The authors claim that their findings provide "valuable input for large-scale calculations of radiative capture processes relevant to nuclear astrophysics." Specifically, the derived analytical forms offer a computationally efficient and physically motivated description of the IVGDR strength function, which is essential for modeling photon-induced and radiative capture reaction rates in astrophysical environments (e.g., neutron star mergers).

The paper emphasizes that the extracted symmetry energy parameters are consistent with independent measurements of dipole polarizability, thereby strengthening the reliability of the inferred nuclear equation of state. Furthermore, the established correlation between $\alpha_D J$ and $r_{skin}$ offers a "practical pathway" to estimate neutron-skin thickness across a broad mass range, which is crucial for constraining the density dependence of the symmetry energy and understanding the structure of neutron stars. The work is presented as a bridge between IVGDR observables and fundamental isovector properties, improving the understanding of astrophysical scenarios under extreme conditions.

Technical Summary: Physics-informed Interpretable Forms for Dipole Resonances and Quantities Relevant for Astrophysics

More like this