Constraining the nuclear symmetry energy from electric dipole polarizability and neutron skin in $^{208}\mathrm{Pb}$ within antisymmetrized molecular dynamics

Using the antisymmetrized molecular dynamics model, this study constrains the nuclear symmetry energy at subsaturation densities (0.2$\rho_0$ to 0.57$\rho_0$) by analyzing the electric dipole polarizability and neutron skin thickness of $^{208}\mathrm{Pb}$, favoring effective interaction parameters of $S_0\sim34$ MeV and $L=66-75$ MeV.

The Gist

Imagine the atomic nucleus not as a static marble, but as a bustling, crowded dance floor inside a tiny ball. On this floor, you have two types of dancers: Protons (who are positively charged and push each other away) and Neutrons (who are neutral and act as the glue).

In a perfect, balanced ball (like Lead-208, the subject of this study), there are roughly equal numbers of both. But in the universe, things aren't always balanced. In the cores of neutron stars, for example, there are way more neutrons than protons.

The big question physicists have been asking for decades is: How does the "glue" behave when the dance floor gets crowded with too many neutrons? This "glue" is called the Nuclear Symmetry Energy. It's the force that tries to keep the protons and neutrons happy and balanced. If this force is too weak, the star might collapse; if it's too strong, it might explode.

The Problem: We Can't See Inside

The trouble is, we can't stick a ruler inside a neutron star or a tiny atomic nucleus to measure this force directly. It's like trying to guess how stiff a mattress is just by looking at the box it came in.

Scientists have been using two different "probes" to guess the stiffness of this glue:

1. The Neutron Skin ($\Delta R_{np}$): Imagine the nucleus as a peach. The "skin" is the layer of neutrons sticking out on the outside because there are so many of them. A thicker skin suggests the glue is "stiff" (pushing neutrons out).
2. The Electric Dipole Polarizability ($\alpha_D$): Imagine poking the nucleus with a gentle electric stick. If the nucleus is "squishy," the protons and neutrons will wobble apart easily. If it's "stiff," they resist the wobble. This wobble is called the Giant Dipole Resonance.

The Conflict

For a long time, different theories gave different answers. Some said the glue was soft (thin skin, easy wobble), while others said it was stiff (thick skin, hard wobble). It was like two mechanics looking at the same car engine and arguing about whether the springs were soft or hard.

The New Approach: The "Anti-Symmetrized Molecular Dynamics" (AMD)

The authors of this paper decided to run a super-accurate computer simulation to settle the argument. They used a model called AMD.

Think of AMD as a high-definition, slow-motion video game of the nucleus.

• The Old Models: Some older models were like a blurry, low-resolution game where the dancers (nucleons) were just dots moving around. They sometimes forgot the rules of quantum mechanics (the "Pauli Exclusion Principle," which says no two dancers can stand in the exact same spot).
• The AMD Model: This model treats every dancer as a fuzzy, wiggly cloud (a Gaussian wave packet) and ensures they strictly follow the quantum rules. It's like a high-fidelity simulation where every dancer has a personal space bubble that never overlaps with another's.

What They Did

They simulated the Lead-208 nucleus thousands of times, tweaking the "rules of the game" (the parameters of the symmetry energy) to see which version matched reality.

They looked at two things in their simulation:

1. How thick the neutron skin was.
2. How the nucleus wobbled when poked (the Giant Dipole Resonance).

The Discovery: Finding the "Sweet Spot"

By running their high-definition simulation, they found a specific set of rules that made the simulation match the real-world data perfectly.

• The Skin: They found that to get the right thickness of the neutron skin (matching the famous PREX-II experiment), the "glue" needs to be moderately stiff.
• The Wobble: They found that to get the right amount of wobble (matching the RCNP experiment), the glue needs to be stiff in a specific way at low densities.

The Analogy: Imagine you are trying to tune a guitar string.

• If the string is too loose (soft symmetry energy), the skin is too thick, but the wobble is too slow.
• If the string is too tight (stiff symmetry energy), the skin is too thin, but the wobble is too fast.
• The authors found the perfect tuning. They determined that the "glue" has a specific strength at normal density ($S_0 \approx 34$ MeV) and a specific rate of change as you move away from normal density ($L \approx 66-75$ MeV).

Why This Matters

This isn't just about a lead atom in a lab.

• Neutron Stars: These findings tell us how neutron stars behave. A stiffer glue means neutron stars can be heavier and larger before collapsing into black holes.
• Supernovas: It helps explain how stars explode.
• The Universe: It refines our understanding of how matter behaves under extreme conditions, from the center of an atom to the edge of the observable universe.

The Bottom Line

The authors used a super-precise, quantum-mechanical computer simulation to act as a "universal translator" between two different experiments. They successfully narrowed down the properties of the nuclear "glue," resolving a long-standing tension in physics.

They essentially said: "We've measured the wobble and the skin, and our high-definition simulation tells us exactly how the nuclear glue must behave to make both of those things happen at the same time."

The result? A much clearer picture of the universe's most extreme matter.

Technical Summary

1. Problem Statement

The nuclear symmetry energy, $S(\rho)$, characterizes the isospin dependence of the nuclear equation of state (EOS). While crucial for understanding neutron stars, heavy-ion collisions, and nuclear structure, its density dependence, particularly at subsaturation densities, remains highly uncertain.

• Key Observables: The electric dipole polarizability ($\alpha_D$) and the neutron skin thickness ($\Delta R_{np}$) of $^{208}\text{Pb}$ are considered "clean" probes for constraining $S(\rho)$ at subsaturation densities.
• The Tension: Previous studies using different theoretical frameworks (e.g., Random Phase Approximation with Energy Density Functionals vs. Transport models) have yielded conflicting constraints. Specifically, many models struggle to simultaneously reproduce the PREX-II neutron skin data and the RCNP dipole polarizability data.
• Model Limitations: Semi-classical transport models (like BUU or standard QMD) often fail to accurately describe the neutron skin of $^{208}\text{Pb}$ due to insufficient Pauli blocking. Conversely, mean-field approaches (RPA) often lack the dynamical treatment of fragment formation and cluster correlations found in transport models.

2. Methodology

The authors employ the Antisymmetrized Molecular Dynamics (AMD) model, a microscopic transport model that treats nucleons as Gaussian wave packets and preserves the fermionic character of the system via a Slater determinant wave function.

• Theoretical Framework:

  • Wave Function: The system is described by a Slater determinant of localized Gaussian wave packets.
  • Equations of Motion: Derived from the time-dependent variational principle.
  • Effective Interaction: The authors use a Skyrme-type effective interaction with an extended density-dependent term. This allows them to systematically vary the symmetry energy coefficient at saturation ($S_0$) and the slope parameter ($L$) while keeping the incompressibility ($K_0$) and effective masses ($m^*_s, m^*_v$) fixed.
  • Parameter Sets: 30 distinct parameter sets were generated, varying $S_0$ (30, 32, 34 MeV), $L$ (46–108 MeV), and the neutron-proton effective mass splitting ($f_I = m/m^*_s - m/m^*_v$), covering cases where $m^*_n < m^*_p$ and $m^*_n > m^*_p$.

• Calculation Procedures:

  1. Ground State: The ground state of $^{208}\text{Pb}$ is prepared using a frictional cooling method. The study accounts for "initial state fluctuations" by averaging over 200 events.
  2. IVGDR Simulation: An external electric dipole perturbation is applied to the ground state. The system's time evolution is tracked to calculate the dipole moment $D(t)$.
  3. Strength Function: The transition strength distribution $S(\omega)$ is obtained via Fourier transform of the dipole moment.
  4. Observables:
     • $\alpha_D$: Calculated by integrating the strength function ($\alpha_D = \frac{8\pi}{9} \int \frac{S(\omega)}{\omega} d\omega$).
     • $\Delta R_{np}$: Calculated from the root-mean-square radii of neutrons and protons.

3. Key Contributions

• Unified Framework: This work is the first to simultaneously describe the IVGDR and neutron skin of $^{208}\text{Pb}$ within a single microscopic transport model (AMD) that fully respects the Pauli principle.
• Landau Damping Mechanism: The authors demonstrate that the AMD model naturally reproduces the width of the Giant Dipole Resonance (GDR) without introducing artificial smoothing parameters or two-body collisions. This width arises from Landau damping (phase mixing) caused by the inhomogeneity of local mean fields and effective masses experienced by individual wave packets.
• Sensitivity Analysis: The study quantifies the sensitivity of $\alpha_D$ and $\Delta R_{np}$ to $S_0$, $L$, and the effective mass splitting ($f_I$), revealing distinct dependencies that help break degeneracies in parameter space.

4. Key Results

A. Ground State and Neutron Skin

• The AMD model successfully reproduces experimental binding energies and charge radii.
• $\Delta R_{np}$ Dependence: The neutron skin thickness is strongly correlated with the slope parameter $L$ but weakly dependent on $S_0$ and effective mass splitting.
• Constraint: To reproduce the PREX-II experimental value ($\Delta R_{np} \approx 0.283$ fm), the model requires $L > 67$ MeV.

B. IVGDR and Polarizability

• Centroid Energy: The centroid energy of the IVGDR increases with $S_0$ and the magnitude of effective mass splitting ($|f_I|$), but decreases with increasing $L$. This counter-intuitive behavior (larger $L$ $\to$ lower energy) is attributed to the fact that a stiffer symmetry energy (larger $L$) implies a softer symmetry energy at subsaturation densities, reducing the restoring force.
• $\alpha_D$ Dependence: $\alpha_D$ shows a robust linear dependence on $L$ but also a pronounced dependence on $S_0$ and $f_I$.
• $\chi^2$ Analysis: A parameter set with $S_0 = 30$ MeV and $L = 61$ MeV provided the best fit to the experimental strength function $S(E)$ ($\chi^2_r = 14.17$). However, this set underestimates the neutron skin.

C. Combined Constraints on Symmetry Energy

By simultaneously fitting both $\alpha_D$ (RCNP data) and $\Delta R_{np}$ (PREX-II data), the authors identified a favored region in the parameter space:

• Favored Parameters: $S_0 \approx 34$ MeV and $L \approx 66\text{--}75$ MeV.
• Density Sensitivity: Correlation analysis revealed that both observables are most sensitive to the symmetry energy in the density range $0.20\rho_0 \leq \rho \leq 0.57\rho_0$.
• Extracted Values:
  • $S(\rho = 0.20\rho_0) = 10.5 \pm 0.63$ MeV
  • $S(\rho = 0.57\rho_0) = 23.1 \pm 0.40$ MeV

These constraints are consistent with results from chiral effective field theory and other experimental probes (isospin diffusion, pion ratios) within uncertainties.

5. Significance

• Resolution of Tension: The study provides a unified microscopic explanation for the correlation between $\alpha_D$ and $\Delta R_{np}$, resolving tensions found in previous density functional theory studies.
• Validation of AMD: It validates the AMD model as a powerful tool for nuclear structure and reaction studies, specifically highlighting its ability to handle fermionic statistics and collective excitations without ad-hoc damping mechanisms.
• Precision Constraints: The derived constraints on $S(\rho)$ in the subsaturation region provide critical input for nuclear astrophysics, particularly for modeling the Equation of State of neutron-rich matter in neutron stars.
• Future Directions: The authors note a slight discrepancy between the parameters favoring the strength function shape ($S_0=30$) and those favoring the combined skin/polarizability data ($S_0=34$), suggesting that further exploration of interaction parameter spaces and dynamical ingredients is necessary for a fully quantitative description.