⚛️ nuclear theory

Sensitivity of fission barrier heights to nuclear matter properties

Using covariant density functional theory and Gaussian-process emulators, this study demonstrates a strong correlation between the surface energy coefficient and fission barrier heights in $^{240}$ Pu and $^{236}$ U, identifying the isoscalar-vector derivative coupling constant as the primary factor governing the outer barrier height.

Original authors: X. Zhang, K. Hagino, J. M. Yao

Published 2026-10-06

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: X. Zhang, K. Hagino, J. M. Yao

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

At the heart of the atomic nucleus lies a delicate balance of forces that determines whether an atom remains stable or splits apart. Inside the nucleus, protons and neutrons are held together by a powerful attraction, but the protons also push each other apart because they all carry a positive electric charge. This tug-of-war creates a kind of energetic landscape. For a nucleus to split, or undergo fission, it must climb over a hill of energy known as a barrier. If the hill is high, the nucleus stays intact; if the hill is low, the nucleus can easily break apart. Understanding the height of these hills is crucial for everything from generating nuclear power to explaining how heavy elements are forged in the violent explosions of stars. However, predicting exactly how high these barriers are is difficult because it requires calculating the behavior of thousands of particles interacting in complex ways.

A new study by researchers at Kyoto University and Sun Yat-sen University has taken a fresh look at this problem by connecting the behavior of individual nuclei to the properties of nuclear matter in general. The team focused on two heavy elements, plutonium-240 and uranium-236, which are classic examples of atoms that can split. They wanted to know how the "surface energy" of a nucleus—essentially the cost of creating a new surface area when the nucleus stretches and deforms—relates to the height of the fission barrier. In the liquid drop model, a long-standing way of thinking about nuclei, this surface energy is a key factor. The researchers used a sophisticated computer framework called covariant density functional theory, which treats the nucleus as a cloud of particles governed by specific mathematical rules, to simulate how these atoms behave.

To make sense of the vast number of possible outcomes, the scientists employed a statistical tool known as a Gaussian-process emulator. Instead of running a single calculation, they generated thousands of variations of their theoretical model, each with slightly different settings for the fundamental forces inside the nucleus. They then used the known properties of plutonium and uranium, such as their size and binding energy, to filter out the unrealistic variations and focus on the most probable ones. This allowed them to map out the likely range of fission barrier heights and see how they correlated with the surface energy coefficient. The results revealed a strong, direct link: a higher surface energy coefficient leads to higher fission barriers. This makes physical sense, as a stiffer surface requires more energy to stretch and deform, effectively raising the hill the nucleus must climb to split.

The study went further to identify exactly which microscopic settings in their theory were responsible for these results. They found that the height of the outer fission barrier—the second, larger hill the nucleus must cross—is governed almost entirely by a single specific setting in their model, known as the isoscalar-vector derivative coupling. This setting controls how the nuclear force changes as the density of the nucleus varies. In contrast, the inner barrier, which is the first hill the nucleus encounters, is influenced by a combination of several different settings working together. This distinction is important because it tells physicists that different parts of the fission process are sensitive to different aspects of the nuclear force. While the outer barrier is a straightforward reflection of one specific interaction, the inner barrier is a more complex result of multiple factors.

The researchers also examined how the surface energy relates to other properties, such as the fissility parameter, which measures how likely a nucleus is to split. They found that the surface energy coefficient is the dominant driver of this likelihood, rather than the Coulomb energy or other symmetry effects. This confirms that the ability of the nucleus to resist deformation is primarily a matter of its surface tension. The study concludes that by measuring fission properties in the laboratory, scientists can now place tighter constraints on the surface energy coefficient, a value that is difficult to measure directly. This improved understanding could refine our models of nuclear astrophysics, particularly in the extreme environments of neutron stars, where the surface properties of matter dictate how the star's crust behaves and evolves.

Technical Summary: Sensitivity of Fission Barrier Heights to Nuclear Matter Properties

Problem Statement
Nuclear fission is a critical process in energy production, superheavy element synthesis, and astrophysical nucleosynthesis. A central quantity in describing fission is the fission barrier height, which requires a consistent microscopic treatment of large-amplitude collective deformation and single-particle shell structure. While Nuclear Density Functional Theory (DFT) has successfully described fission across the nuclear chart, the connection between bulk nuclear matter properties and fission observables remains less explored than other correlations (e.g., neutron-skin thickness and symmetry energy). Specifically, the dependence of fission barrier heights on the surface energy coefficient ( $a_{\text{surf}}$ ) and the underlying microscopic couplings of the Energy Density Functional (EDF) lacks a comprehensive statistical analysis. Previous studies established a linear correlation between $a_{\text{surf}}$ and barrier heights using Skyrme parametrizations, but a rigorous quantification of parametric uncertainty and the identification of the specific microscopic drivers of this correlation were missing.

Methodology
The authors employ Covariant Density Functional Theory (CDFT) using the point-coupling functional PC-PK1 to investigate the microscopic origin of correlations between surface energy and fission properties. The study focuses on two representative actinide nuclei: $^{240}\text{Pu}$ and $^{236}\text{U}$ .

Theoretical Framework:
- Fission Barriers: Potential Energy Surfaces (PES) are calculated within the relativistic mean-field (RMF) theory, including pairing correlations (BCS approximation), microscopic center-of-mass corrections, and rotational corrections. The PES is generated as a function of quadrupole deformation ( $\beta_2$ ) using axially deformed harmonic oscillator bases.
- Surface Energy Extraction: The isoscalar ( $a_{\text{surf}}$ ) and isovector ( $a_{\text{sym}}$ ) surface energy coefficients are extracted from Semi-Infinite Nuclear Matter (SINM) using the relativistic Thomas–Fermi approximation. This allows the characterization of surface properties distinct from infinite nuclear matter (INM) bulk properties.
- Statistical Analysis: To explore the parameter space, the authors vary nine EDF coupling constants ( $C = [\alpha_S, \beta_S, \gamma_S, \delta_S, \alpha_V, \gamma_V, \delta_V, \alpha_{TV}, \delta_{TV}]$ ) using Latin hypercube sampling.
- Emulators and Bayesian Inference: Gaussian-process (GP) emulators, specifically PCA-GP for potential energy curves, are trained on 200 EDF parameter sets to predict observables for a larger ensemble (6000 sets). Bayesian inference is used to construct posterior distributions of these quantities, constrained by the ground-state properties (binding energy and charge radius) of $^{240}\text{Pu}$ and $^{236}\text{U}$ .
- Sensitivity Analysis: A global sensitivity analysis (GSA) based on the Shapley effect is performed to quantify the contribution of each coupling constant to the variance of observables. This method accounts for correlations among the posterior EDF parameters.

Key Contributions and Results

Correlation between $a_{\text{surf}}$ and Fission Barriers: The analysis confirms a strong positive correlation between the surface energy coefficient $a_{\text{surf}}$ and both inner ( $E_A$ ) and outer ( $E_B$ ) fission barrier heights. This is physically consistent: a larger $a_{\text{surf}}$ implies a stiffer nuclear surface, increasing the energy cost of deformation.
Differential Sensitivity of Barriers:
- The inner barrier height ( $E_A$ ) is more tightly correlated with $a_{\text{surf}}$ than the outer barrier. The authors attribute this to the fact that the inner barrier occurs at smaller deformations where the competition between surface and Coulomb energies dominates.
- The outer barrier height ( $E_B$ ) exhibits a weaker correlation with $a_{\text{surf}}$ compared to the inner barrier. This is attributed to the large deformation and neck formation at the outer barrier, where curvature effects not captured by the simple $a_{\text{surf}}$ term become significant.
Microscopic Origin (Sensitivity Analysis):
- Outer Barrier ( $E_B$ ): The height of the outer barrier is governed predominantly by the isoscalar-vector derivative coupling constant ( $\delta_V$ ). The Shapley share analysis shows that $\delta_V$ accounts for the vast majority of the variance in $E_B$ , regardless of whether correlations among other coupling constants are retained.
- Inner Barrier ( $E_A$ ) and $a_{\text{surf}}$ : These quantities are jointly governed by $\delta_V$ and a set of mutually correlated coupling constants ( $\beta_S, \gamma_S, \alpha_V$ ). The variance in these quantities is explained by the combined effect of these four parameters.
- Other Observables: The Coulomb energy ( $E_{\text{Coul}}$ ) shows negligible sensitivity to $\delta_V$ , while the surface symmetry coefficient ( $a_{\text{sym}}$ ) is dominated by the isovector derivative coupling ( $\delta_{TV}$ ).
Uncertainty Quantification: The study quantifies the parametric uncertainty of fission observables. The relative uncertainty of the outer barrier ( $E_B$ ) is found to be approximately 25%, significantly larger than that of the inner barrier or ground-state properties. Constraining the EDF with ground-state properties alone provides limited additional information on fission-related quantities compared to the prior, suggesting that experimental fission data (e.g., half-lives) would be necessary to further constrain these parameters.

Significance and Claims
The paper claims to provide the first comprehensive statistical analysis linking fission barrier heights to the surface energy coefficient and the specific microscopic couplings within a covariant point-coupling framework.

Microscopic Identification: The work identifies $\delta_V$ as the primary driver for the outer fission barrier, offering a specific microscopic target for refining EDFs.
Constraint Strategy: The strong correlation between fission properties and $a_{\text{surf}}$ suggests that experimental fission data could be used as a constraint to determine the surface energy coefficient more precisely.
Astrophysical Relevance: The authors note that a more precise determination of $a_{\text{surf}}$ is relevant for nuclear astrophysics, particularly in modeling shape transitions in the inner crust of neutron stars where "nuclear pasta" phases may emerge.
Future Directions: The study is limited to axially symmetric deformations. The authors suggest that extending the framework to include triaxial deformations and a broader set of fissioning nuclei is a necessary future step.

The authors maintain a modest tone regarding the scope, acknowledging that while ground-state constraints reduce uncertainties in bulk properties, they are insufficient to fully constrain fission barriers without incorporating specific fission observables.

Technical Summary: Sensitivity of Fission Barrier Heights to Nuclear Matter Properties

More like this