A unified classification-quantification framework for bubble-like nuclei within the extended quantum molecular dynamics model

This paper presents a unified classification-quantification framework using the dimensionless $BHTU$ parameters to systematically characterize bubble-like nuclear morphologies across the AME2020 database within the extended quantum molecular dynamics model, revealing widespread bubble and toroidal structures in medium-mass, heavy, and superheavy nuclei.

The Gist

Imagine the atomic nucleus not as a smooth, solid marble, but as a dynamic, shifting cloud of tiny particles (protons and neutrons) that can arrange themselves in all sorts of strange shapes. For a long time, scientists thought these clouds were mostly uniform balls. But this paper suggests that under certain conditions, these clouds can puff up in the middle, creating hollow spaces, or even form ring shapes, much like a donut.

Here is a simple breakdown of what the researchers did and what they found, using everyday analogies:

The Problem: Mapping a Shape-Shifting Universe

Think of the "Periodic Table" of elements as a giant map. Scientists have known about some weird shapes on this map (like "bubbles" where the center is empty), but they only knew about a few specific islands. They didn't have a complete map of where these strange shapes appear, nor did they have a standard ruler to measure exactly how hollow or thick a nucleus is.

The Tool: A "Frictional Cooling" Simulation

The researchers used a computer model called EQMD (Extended Quantum Molecular Dynamics).

• The Analogy: Imagine you have a bowl full of marbles (the protons and neutrons) that are vibrating wildly. If you just let them be, they bounce around chaotically. To see their natural, resting shape, you need to slow them down.
• The Method: The researchers added a "frictional cooling" mechanism to their simulation. Think of this like putting the vibrating marbles in a thick, cold syrup. It slows them down gently until they settle into their most stable, relaxed arrangement. This allowed them to see the "true" shape of the nucleus without the noise of constant shaking.

The Discovery: Three Main Shapes

After cooling down thousands of different nuclei, the researchers found that the nuclei generally fell into three categories, which they named based on their shape:

1. The Droplet (B = 0):

   • What it is: A standard, solid ball. The density is highest in the center and fades out toward the edge, just like a drop of water.
   • Where they are: Mostly found in light nuclei (small atoms).

2. The Bubble (B = 1):

   • What it is: A hollow ball. The center is empty or very thin, and the matter is packed into a shell around the outside.
   • Where they are: Mostly found in medium-sized nuclei. The researchers highlighted a specific area around the element Calcium-40 and neutron-rich areas as "prime candidates" where these bubbles are most likely to be found.

3. The Toroidal Bubble (B = 2):

   • What it is: A donut or a ring. The density dips in the very center, rises up in a ring in the middle, and then dips again before the outer edge.
   • Where they are: These start appearing in heavier nuclei (around atomic number 25) and become common in very heavy, super-heavy elements.

The New "Ruler": The B-H-T-U Framework

To stop guessing and start measuring, the team created a unified classification system using four "factors" (like a scorecard for nuclear shapes):

• B (The Shape Score): This counts the "bumps" in the density curve.
  • 0 bumps = Droplet.
  • 1 bump = Bubble.
  • 2 bumps = Toroidal Bubble.
• H (The Hollow Score): This measures how empty the center is. A high score means a very hollow center; a low score means a solid center.
• T (The Thickness Score): This measures how thick the "skin" or outer layer of the nucleus is.
• U (The Bubble Size Score): This measures how big the empty hole in the middle is compared to the whole nucleus.

What They Found on the Map

By applying this new ruler to the entire map of known elements (from the AME2020 database), they created a visual guide:

• Light elements are mostly solid droplets.
• Medium elements (like the area around Calcium) are the "bubble capital," showing the most significant hollow centers.
• Heavy elements start turning into "donuts" (toroidal bubbles).
• Super-heavy elements also show widespread bubble structures.

Why This Matters (According to the Paper)

The paper claims this work does two main things:

1. It reveals the "richness" of nuclear shapes: It shows that nuclei are much more varied than just solid balls; they can be hollow, ring-shaped, and everything in between.
2. It provides a predictive tool: By using this B-H-T-U framework, scientists now have a standardized way to predict which specific atoms might have these exotic shapes. This gives experimentalists a "treasure map" to know exactly where to look in future experiments to find these bubble-like structures.

In short, the researchers built a new way to sort and measure the shapes of atomic nuclei, discovering that "hollow" and "ring" shapes are much more common in nature than previously mapped, especially in medium and heavy elements.

Technical Summary

Technical Summary: A Unified Classification-Quantification Framework for Bubble-Like Nuclei within the Extended Quantum Molecular Dynamics Model

Problem Statement
The structural diversity of atomic nuclei, particularly the existence of exotic shapes such as bubble, toroidal, and halo configurations, remains a significant challenge in nuclear physics. While traditional semiclassical liquid drop models assume a dense, nearly spherical system with uniform central density, theoretical and experimental evidence suggests that under specific conditions—such as large isospin asymmetry or high angular momentum—nuclei can exhibit significantly reduced central densities. Previous investigations into these "bubble" structures have largely relied on static mean-field theories (e.g., Self-Consistent Mean Field, Hartree-Fock-Bogoliubov, Energy Density Functional). However, these approaches often focus on specific nuclides (e.g., $^{34}$Si, $^{36}$Ar) or localized regions of the nuclear chart, lacking a systematic census across all known nuclides. Furthermore, existing frameworks often exhibit inconsistencies in predicting bubble nuclei when accounting for pairing correlations and deformation. There is a distinct need for a unified classification and quantification framework capable of mapping bubble-like structures across the entire nuclear chart, particularly within the context of cluster configurations where local density variations are naturally pronounced.

Methodology
This study employs the Extended Quantum Molecular Dynamics (EQMD) model to perform a systematic investigation of relaxed low-energy cluster configurations for all nuclides listed in the AME2020 database. Unlike static mean-field approaches, EQMD treats nucleons as Gaussian wave packets with complex width parameters that evolve dynamically under a Hamiltonian comprising kinetic energy, center-of-mass corrections, and effective interactions (Skyrme, Coulomb, symmetry, and Pauli terms).

Key methodological features include:

1. Frictional Cooling: To stabilize the system and obtain relaxed low-energy states, a frictional cooling mechanism is integrated into the equations of motion. This reduces statistical fluctuations and allows the system to settle into local energy minima, favoring the emergence of clustered structures.
2. Initialization and Relaxation: Nucleon positions are sampled uniformly within a sphere, and momenta are drawn via Monte Carlo using the local Fermi gas approximation. The system evolves via momentum dissipation equations until the time derivatives of coordinates, momenta, and width parameters vanish, indicating a relaxed state.
3. Density Analysis: The radial density distribution is derived by integrating the phase-space distribution function. To characterize these distributions, the authors introduce a unified framework based on dimensionless parameters derived from the radial density profile.

Key Contributions and Framework
The primary contribution of this work is the establishment of a unified classification-quantification framework based on four dimensionless parameters: B, H, T, and U.

• Classification Factor ($B$): Determined by the number of inflection points in the radial density profile, $B$ categorizes nuclei into three morphological types:
  • $B = 0$: Droplet (monotonically decreasing density).
  • $B = 1$: Bubble (pronounced central density depletion).
  • $B = 2$: Toroidal Bubble (density decreases then increases with radius, indicating a local minimum at intermediate radii).
• Quantification Factors:
  • $H$ (Degree of Central Depletion): Defined as $H = (\rho_{max} - \rho_{center}) / \rho_{max}$, quantifying the relative reduction of central density.
  • $T$ (Relative Surface Thickness): Characterizes the thickness of the surface liquid-like layer relative to the root-mean-square radius ($R_{rms}$).
  • $U$ (Relative Bubble Size): Characterizes the relative size of the internal low-density region.

This framework allows for a systematic mapping of nuclear shapes across the nuclide chart, moving beyond simple classification to a quantitative description of density distributions.

Results
The systematic calculations reveal distinct distribution patterns for different nuclear morphologies:

• Light Nuclei: Predominantly exhibit droplet-like structures ($B=0, H=0, T=1, U=0$). The limited number of nucleons and strong overlap of wave packets prevent the formation of distinct low-density central regions.
• Medium-Mass Nuclei: Bubble structures ($B=1$) are prevalent, particularly in the vicinity of $^{40}$Ca and in neutron-rich regions. These nuclei display pronounced central hollowing with large $H$ and $U$ values. The region near $^{40}$Ca is identified as a prime candidate for experimental searches due to high stability and significant central depletion.
• Heavy and Superheavy Nuclei: Toroidal bubble structures ($B=2$) emerge around $Z \approx 25$ and become increasingly prevalent in heavier systems ($Z > 28$). These configurations feature a local density minimum at intermediate radii and a shell-like low-density region.
• Superheavy Region: Bubble structures are found to be widespread, consistent with earlier theoretical predictions.
• Proton vs. Neutron Distributions: The appendix details that proton bubbles are primarily distributed in symmetric nuclei and regions with slight neutron excess, while neutron bubbles appear on both sides of the chart, reflecting isospin-symmetric pairing within alpha clusters and surface delocalization effects.

Significance and Claims
The authors claim that this work overcomes the limitations of traditional mean-field theories, which rely on smooth, homogeneous density distributions, by capturing cluster configurations through the EQMD framework. The study asserts that:

1. It provides the first systematic census of bubble-like structures across the entire nuclear chart using a dynamic cluster model.
2. The proposed BHTU parametric framework offers a robust theoretical tool for identifying and predicting bubble-like nuclei, distinguishing between different density characteristics and revealing systematic trends as functions of mass number and isospin asymmetry.
3. By focusing on relaxed low-energy cluster states, the study unveils a rich landscape of nuclear shapes that may be experimentally accessible through dynamical nonequilibrium states (e.g., giant resonances) in heavy-ion collisions.
4. The results provide theoretical guidance for future experimental efforts at radioactive ion beam facilities, specifically targeting the regions identified as having high probabilities of bubble formation (e.g., near $^{40}$Ca and in the superheavy region).

The paper maintains a modest stance, noting that while the global scan is influenced by specific model assumptions and cluster-structure characteristics, the results serve as a useful reference and a foundation for exploring the complexity of nuclear density distributions.