Computing nuclear response functions with time-dependent coupled-cluster theory

This paper validates a time-dependent coupled-cluster theory approach for computing nuclear response functions by demonstrating its accuracy against static frameworks and revealing collective dipole resonances and chaotic behavior in light nuclei under strong electric fields.

The Gist

The Big Picture: Listening to the Nucleus Sing

Imagine an atomic nucleus not as a static rock, but as a tiny, vibrating drum made of protons and neutrons. When you hit this drum, it vibrates in specific patterns. In physics, these vibrations are called resonances.

For decades, scientists have been trying to predict exactly how these "nuclear drums" sound. They usually do this by taking a "snapshot" of the drum while it's vibrating (a static approach). But this paper introduces a new way of listening: Time-Dependent Coupled-Cluster (TDCC) theory.

Instead of taking a snapshot, the authors built a super-advanced movie camera. They simulate the nucleus being hit by a tiny electric "tap," and then they watch, frame-by-frame, how the protons and neutrons dance, wiggle, and oscillate over time. By recording this movie and running it backward through a mathematical filter (a Fourier transform), they can figure out exactly what notes the nucleus is playing.

The Cast of Characters

• The Nucleus: Think of it as a crowded dance floor. The protons and neutrons are the dancers. They are constantly bumping into each other and holding hands in complex ways (these are called "correlations").
• The "Tap" (Perturbation): The scientists give the dance floor a tiny, quick shove using an electric field. It's like flicking a drum skin.
• The Movie Camera (TDCC): This is the new method. It tracks every single dancer's movement over time, accounting for how they influence each other, rather than just assuming they move in a simple, average pattern.

What They Did (The Experiment)

The team tested their new "movie camera" on three specific nuclei: Helium-4 (very small), Oxygen-16 (medium), and Oxygen-24 (neutron-rich).

1. The Check-Up (Validation): First, they wanted to make sure their movie camera was accurate. They compared the "movie" results against the old "snapshot" results (which are already known to be very good for simple cases).

   • The Result: The movie and the snapshot matched almost perfectly. This proved their new method works and is reliable.

2. Watching the Dance (Density Fluctuations): This is where the movie camera shines. In the old snapshot method, you can't really see how the dancers move relative to each other.

   • The Giant Dipole Resonance (GDR): In Oxygen-16, they saw the protons and neutrons moving in perfect opposition, like a seesaw. The protons go left while the neutrons go right, then they switch. This is a classic, collective dance move.
   • The Pygmy Dipole Resonance (PDR): In Oxygen-24 (which has extra neutrons), they saw something cooler. The "core" of the nucleus (protons and neutrons) stayed mostly still, while the "skin" of extra neutrons wiggled back and forth like a loose jacket. This is a specific dance move that only happens in heavy, neutron-rich nuclei.

3. The Chaos Test (Non-Linear Regime): Finally, they decided to turn up the volume. Instead of a tiny tap, they hit the nucleus with a massive electric field.

   • The Result: When the hit was too hard, the orderly dance turned into chaos. The neat, repeating patterns broke down, and the system started behaving unpredictably (chaotically). It's like hitting a drum so hard it starts rattling and making weird, random noises instead of a clear tone. Interestingly, this chaotic behavior matched what other, simpler theories predicted, giving scientists confidence that even in chaos, the laws of physics hold up.

Why This Matters

Why should you care about watching nuclear movies?

• Stellar Alchemy: The universe creates heavy elements (like gold and uranium) in the hearts of dying stars and exploding supernovas. To understand how these elements are made, we need to know exactly how nuclei react to energy. This new method helps us predict those reactions more accurately.
• The Future of Energy: Understanding how nuclei vibrate and break apart is crucial for nuclear fusion (clean energy) and fission (current nuclear power).
• Beyond the Average: Old methods often treated the nucleus like a smooth, average blob. This new method sees the individual "dancers" and how they interact, giving us a much sharper, more realistic picture of the quantum world.

The Bottom Line

The authors successfully built a high-definition "time machine" for nuclear physics. They proved it works by showing it matches old methods, but they also showed it can see things the old methods can't: the specific choreography of protons and neutrons, and how the nucleus behaves when pushed to its absolute limits. It's a step forward from taking a blurry photo to watching a high-definition documentary of the atomic world.

Technical Summary

1. Problem Statement

Accurate prediction of nuclear dynamical processes (e.g., fusion rates, neutron capture, fission) is crucial for astrophysics, particularly for nuclei far from stability where experimental data is scarce.

• Limitations of Current Methods: Time-Dependent Density Functional Theory (TDDFT) is the standard for nuclear dynamics but relies on mean-field approximations, failing to capture quantum many-body correlations (e.g., tunneling below barriers).
• Limitations of Static Ab Initio Methods: While static ab initio approaches (like Coupled-Cluster, CC) successfully describe ground states and some excited states in medium-mass nuclei (up to $^{208}$Pb), they struggle with continuum states and real-time dynamical evolution.
• The Gap: There is a need for a method that combines the high accuracy of ab initio many-body correlations with the ability to simulate real-time nuclear dynamics and response functions.

2. Methodology

The authors develop and apply Time-Dependent Coupled-Cluster (TDCC) theory to solve the time-dependent $A$-body Schrödinger equation.

• Theoretical Framework:

  • Ansatz: The nuclear wavefunction is evolved as $|\Psi(t)\rangle = e^{T(t)} |\Phi_0\rangle$, where $|\Phi_0\rangle$ is a static Hartree-Fock reference state and $T(t)$ is the cluster operator truncated at the CCSD level (Singles and Doubles).
  • Bivariational Approach: To calculate expectation values (like transition moments) with non-Hermitian similarity-transformed Hamiltonians, a left state $\langle \tilde{\Psi}(t)| = \langle \Phi_0| L(t) e^{-T(t)}$ is introduced, where $L(t)$ is the de-excitation operator (also truncated at CCSD).
  • Perturbation: A time-dependent perturbation $W(t) = \Theta \delta(t)$ (or a Gaussian pulse) is applied to the ground state. The response function $R(E)$ is extracted via the Fourier transform of the time-dependent transition moment $\Theta(t) = \langle \tilde{\Psi}(t)| \Theta |\Psi(t)\rangle$.
  • Equations of Motion: The time evolution of the cluster amplitudes ($T$) and de-excitation amplitudes ($L$) yields a system of coupled ordinary differential equations (ODEs).

• Computational Implementation:

  • Hamiltonian: Uses the chiral nucleon-nucleon interaction NNLOopt. Three-nucleon forces are currently excluded.
  • Solver: The ODE system is stiff, especially for larger nuclei. The authors utilize the CVODE solver from the SUNDIALS library, employing adaptive time-stepping (Adams-Moulton for light nuclei, BDF methods for heavier ones) and complex-valued vector data structures.
  • Post-Processing: The time signal is windowed (using a cosine function) before applying a Discrete Fourier Transform (DFT) to minimize spectral artifacts and extract energy spectra.

3. Key Contributions

• Validation of TDCC: The paper provides a rigorous benchmark by comparing TDCC results against equivalent static calculations using the Lorentz Integral Transform (LIT) method combined with CC.
• Extension to Medium-Mass Nuclei: While previous TDCC applications were limited to spherical symmetry or very light nuclei, this work pushes the method to medium-mass nuclei ($^{16}$O, $^{24}$O) while maintaining axial symmetry.
• Real-Time Density Evolution: The method allows for the direct visualization of proton and neutron density fluctuations in real-time, offering insights into collective modes that are difficult to extract from static frameworks.
• Exploration of Non-Linear Regimes: The framework is used to investigate nuclear behavior under strong electric fields, transitioning from linear response to chaotic dynamics.

4. Key Results

• Benchmarking ($^4$He and $^{16}$O):

  • TDCC results for electric dipole response functions show negligible discrepancies (average deviation $\approx 2\%$, max $\approx 4\%$) compared to static LIT-CC results for sum rules ($\alpha_D$, $m_0$, $m_1$).
  • The Giant Dipole Resonance (GDR) peak positions are stable and consistent with experimental data and other ab initio methods (e.g., NCSM).
  • Convergence with respect to the model space size ($N_{max}$) and simulation time ($T_{max}$) was established. A simulation time of $T_{max} = 2000$ fm/c provides sufficient resolution ($\sim 0.6$ MeV).

• Density Fluctuations and Collective Modes:

  • GDR in $^{16}$O: The time evolution clearly shows protons and neutrons oscillating in counter-phase, visually confirming the GDR as a collective motion.
  • Pygmy Dipole Resonance (PDR) in $^{24}$O: In neutron-rich $^{24}$O, the method isolates low-energy modes where excess neutrons oscillate against a proton-neutron core, visualizing the PDR and the development of a neutron skin.

• Non-Linear and Chaotic Regimes:

  • By increasing the perturbation strength $\epsilon$ to $\sim 100$ MeV/fm, the system enters a non-linear regime.
  • Chaotic Behavior: Phase-space trajectories ($D(t)$ vs. $D(t+\tau)$) reveal the emergence of chaotic dynamics at high field strengths.
  • Spectral Changes: In the non-linear regime, the GDR peak intensity decreases by $\sim 40\%$, and the spectrum becomes more fragmented with enhanced low-energy strength. These qualitative findings align with previous TDDFT results.

5. Significance and Outlook

• Validation of Ab Initio Dynamics: The work proves that TDCC is a viable, accurate alternative to TDDFT for calculating nuclear response functions, capturing correlations beyond the mean field.
• Access to New Physics: The method provides unique access to real-time density evolution and non-linear phenomena (chaos) that are inaccessible to static perturbation theory.
• Future Directions:
  • The current implementation uses a static basis, limiting the description of continuum effects (damping widths, particle emission). Future work will incorporate lattice formulations of CC and time-dependent single-particle orbitals to better handle open quantum systems and nuclear collisions.
  • The study of strong-field regimes suggests potential experimental relevance for future facilities like the Gamma Factory at CERN, where high-energy photons could induce non-linear nuclear responses.

In summary, this paper establishes TDCC as a powerful tool for nuclear structure and reaction theory, successfully bridging the gap between high-precision static ab initio methods and the dynamic evolution of nuclear systems.