Multi-Nucleon Transfer Reactions and the Creation and the Evolution of the Compound Nucleus

This paper introduces the enhanced Generator Coordinate Method (eGCM) as the first microscopic quantum approach capable of describing compound nucleus formation in multi-nucleon transfer reactions by incorporating quantum fluctuations neglected in traditional mean-field theories like TDHF.

The Gist

The Big Picture: Why Do Atoms Stick Together?

Imagine you are trying to build a new, super-heavy Lego castle. You have two smaller castles (nuclei) and you want to smash them together to make one giant, stable one. In the world of nuclear physics, this is called a Compound Nucleus.

For 90 years, scientists have known this happens. They have a "rulebook" (phenomenological models) that predicts how often it happens, but they don't actually understand how the pieces lock together at the deepest, quantum level. It's like knowing a car engine works because you've seen it run for decades, but never having opened the hood to see the pistons moving.

The problem is that the current best tools for looking inside the engine are too "classical." They treat the nucleus like a smooth, predictable fluid. But in reality, the nucleus is a chaotic, quantum mess of particles jumping around. The old tools miss the "quantum fluctuations"—the tiny, random jitters that are actually crucial for the reaction to work.

The Old Way vs. The New Way

The Old Way (TDHF):
Think of the old method (Time-Dependent Hartree-Fock) as watching a single, perfect movie of two cars crashing.

• In this movie, the cars hit, bounce off, or stick together based on a single, smooth path.
• The problem? Real life isn't one smooth path. It's a million different possibilities happening at once. The old method ignores the fact that the cars could vibrate, spin, or wobble in a thousand different ways during the crash. It misses the "quantum noise."

The New Way (eGCM):
The authors, Matthew Kafker and Aurel Bulgac, have invented a new tool called eGCM (enhanced Generator Coordinate Method).

• Instead of watching one movie, imagine you are watching 48,000 different movies of the same crash simultaneously.
• In some movies, the cars hit head-on. In others, they graze each other. In some, they spin wildly; in others, they wobble.
• The eGCM takes all these different "what-if" scenarios (called Slater determinants) and mixes them together like a giant, quantum smoothie. It calculates how all these different possibilities interfere with each other.

The "Nuclear Molasses" Discovery

When the scientists ran their massive simulation (using a supercomputer with 48,000 graphics cards—enough to power a small city!), they found something shocking.

The Old Prediction:
The old "single movie" method said: "When these two nuclei (Calcium-48 and Lead-208) crash, they will touch, swap a few particles, and then bounce apart. They will never stick together to form a heavy, stable compound."

The New Reality:
The new eGCM method said: "Actually, because of all those quantum jitters and the mixing of all those different crash scenarios, there is a 34% chance they will get stuck together!"

They formed a super-heavy nucleus (No-256) that acts like nuclear molasses.

• The Analogy: Imagine throwing two sticky balls of clay at each other. The old method thought they would bounce off like rubber balls. The new method shows that because of the "quantum glue" (interference between all the different ways they could collide), they get stuck in a gooey, long-lasting state.
• This "compound nucleus" is incredibly stable in the simulation. It doesn't break apart even after a very long time. This is the first time a theory has successfully predicted this "sticking" from first principles.

Why Does This Matter?

1. It's a "Quantum Interference" Party: The reason the nuclei stick is due to destructive interference. In quantum mechanics, waves can cancel each other out. Here, the "waves" of the different collision paths cancel out the possibility of them bouncing apart, leaving only the path where they stick together. It's like a choir where everyone sings a different note, but somehow, the noise cancels out to create a perfect, silent pause that traps the particles.
2. Creating New Elements: This is crucial for understanding how new elements are made in labs and in space (like in colliding neutron stars). If we can't predict how nuclei stick together, we can't predict how to make the next element on the periodic table.
3. The "Compound Nucleus" is Real: For decades, the "compound nucleus" was just a theoretical guess (a conjecture). This paper provides the first microscopic proof that it actually exists and explains why it forms.

The Technical "Secret Sauce"

How did they do it?

• They simulated the collision of Calcium and Lead.
• They didn't just look at the start and end; they looked at the collision at thousands of tiny time steps.
• They used a massive "basis set" (a library of possible states) that was 39,000 times larger than anything used before.
• They found that the energy levels of the resulting nucleus matched a famous mathematical pattern called Random Matrix Theory (the Wigner-Dyson surmise). This is the "fingerprint" of a true, chaotic compound nucleus.

The Bottom Line

This paper is a breakthrough because it moves nuclear physics from "guessing based on patterns" to "calculating from the ground up."

They showed that if you account for the chaotic, quantum nature of the nucleus (by mixing millions of possibilities), you discover that heavy nuclei can get stuck together in a "molasses-like" state much more often than we thought. It's the first time we've successfully described the birth of a compound nucleus using the fundamental laws of quantum mechanics, rather than just guessing.

In short: They finally opened the hood of the nuclear engine, turned on the lights, and saw the pistons moving in a way that explains why the engine actually runs.

Technical Summary

1. Problem Statement

The formation of a compound nucleus in heavy-ion collisions has long been a theoretical conjecture based on Niels Bohr's qualitative arguments (1936). Despite its central role in nuclear physics, no microscopic quantum approach based on the time-dependent Schrödinger equation has successfully described this process.

• Limitations of Current Methods: The most advanced existing microscopic method, Time-Dependent Hartree-Fock (TDHF), relies on a mean-field approximation. This approach treats the mean field as a classical expectation value, thereby neglecting crucial quantum fluctuations.
• Limitations of Previous GCM: The Generator Coordinate Method (GCM) and its time-dependent extensions (TDGCM) attempt to include correlations but have historically been limited by assumptions of adiabatic evolution or Gaussian overlap approximations, failing to capture the full complexity of non-equilibrium reactions like Multi-Nucleon Transfer (MNT).
• The Gap: There is a lack of a framework capable of predicting compound nucleus formation cross-sections and explaining the "nuclear molasses" phenomenon (analogous to optical molasses) where nuclei stick together for extended periods.

2. Methodology: Enhanced GCM (eGCM)

The authors introduce and implement the Enhanced Generator Coordinate Method (eGCM), a novel extension of the GCM designed to operate within the continuum for nuclear reactions.

• Core Concept: eGCM is formulated as a Configuration Interaction (CI) framework. It constructs the many-body wavefunction as a superposition of a vast number of linearly independent, non-orthogonal, time-dependent Slater determinants.
• Generator Coordinates: Unlike traditional GCM which uses static shapes, eGCM uses a "bouillabaisse" of trajectories generated by Time-Dependent Density Functional Theory (TDDFT). The generator coordinates include:
  • Impact parameters ($b$): Defining the collision geometry.
  • Time ($\tau$): A crucial coordinate distinguishing eGCM from previous approaches, allowing the mixing of wavefunctions at different stages of the collision trajectory.
• Computational Scale:
  • System: $^{48}\text{Ca} + ^{208}\text{Pb}$ collision at a center-of-mass energy of 235 MeV (slightly above the Coulomb barrier).
  • Data Volume: The study utilized TDDFT simulations storing single-particle wavefunctions for 152 impact parameters and 385 time steps.
  • Resources: Required ~80 TB of memory and utilized 48,000 GPUs on the Frontier supercomputer (Oak Ridge National Lab) to construct and diagonalize a GCM Hamiltonian matrix of size 39,630. This is arguably the largest basis set ever used in a nuclear GCM calculation.
• Mathematical Formulation:
  • The method solves a static eigenvalue problem (Eq. 8) to find eigenstates $|\Psi_n\rangle$ and eigenvalues $E_n$.
  • Time evolution is handled analytically via $|\Psi(t)\rangle = \sum c_n(0) e^{-iE_n t/\hbar} |\Psi_n\rangle$, eliminating the need for a separate time-dependent GCM code.
  • Particle Number Projection: Used to determine the probability of finding specific neutron ($N$) and proton ($Z$) numbers in the reaction fragments.

3. Key Contributions

1. First Microscopic Prediction of Compound Nucleus Formation: The paper presents the first theoretical approach capable of predicting the cross-section for compound nucleus formation from first principles, moving beyond phenomenological models.
2. Validation of the Compound Nucleus Hypothesis: The authors demonstrate that the eGCM spectrum satisfies the Wigner-Dyson surmise (Random Matrix Theory), showing level spacing statistics consistent with the Gaussian Orthogonal Ensemble (GOE). This provides a microscopic derivation of Bohr's compound nucleus concept.
3. Quantum Interference and "Nuclear Molasses": The study reveals that destructive quantum interference among the massive number of configurations leads to the formation of a long-lived "compound-like" state, explaining why certain nuclei do not immediately separate (fission) after contact.
4. Superiority over TDHF and GCMR: The paper establishes that previous methods (TDHF and the Reinhard et al. GCM extension, GCMR) ignore a vast number of relevant many-body configurations. The Inverse Participation Ratio (IPR) analysis shows eGCM wavefunctions are orders of magnitude more complex and mixed than GCMR.

4. Key Results

• Compound Nucleus Formation: In the $^{48}\text{Ca} + ^{208}\text{Pb}$ reaction, eGCM predicts the formation of a compound nucleus $^{256}_{102}\text{No}_{154}$ (Nobelium) with an excitation energy of $\approx 78$ MeV.
  • Probability: The formation probability is calculated at $P \approx 0.34$ (34%).
  • Stability: Unlike TDHF trajectories where nuclei separate within $\approx 2,000$ fm/c, the eGCM compound state is extremely long-lived and does not break apart even in long-time simulations.
• Discrepancy with TDHF:
  • TDHF: Predicts that all trajectories separate before the simulation ends ($t_f = 2,482$ fm/c). It fails to predict the formation of a compound nucleus.
  • eGCM: Predicts significant nucleon transfer in both directions, with a dominant probability ($\approx 0.85$) for the heavy fragment to retain $N \ge 127$ and $Z \ge 83$.
• Level Statistics: The eGCM energy spectrum shows an average level spacing of $\approx 4$ keV. The normalized distribution of level spacings has a standard deviation of 0.5273, closely matching the GOE value of 0.5227.
• IPR Analysis: The IPR values for eGCM are more than two orders of magnitude larger than the maximum possible value for GCMR, confirming that eGCM captures the "phase transition" of nucleon pairs spreading over degenerate orbitals.

5. Significance

• Theoretical Breakthrough: This work bridges the gap between phenomenological models and microscopic quantum mechanics, providing the first rigorous derivation of the compound nucleus from the many-body Schrödinger equation.
• Implications for Nuclear Synthesis: Understanding the mechanism of compound nucleus formation is critical for predicting the creation of superheavy elements and understanding nucleosynthesis in astrophysical events (neutron star mergers, supernovae).
• Methodological Advancement: The successful implementation of eGCM on a massive scale (using 48,000 GPUs) demonstrates that high-performance computing can now tackle the full complexity of non-equilibrium nuclear many-body problems, including quantum fluctuations previously thought intractable.
• New Physics: The identification of "nuclear molasses" suggests a new regime of nuclear dynamics where quantum interference stabilizes transient nuclear systems, challenging the classical view of heavy-ion collisions as purely dissipative or fission-prone processes.

In conclusion, Kafker and Bulgac demonstrate that by fully incorporating quantum fluctuations and time-dependent correlations through the eGCM framework, it is possible to describe the formation and evolution of a compound nucleus, a phenomenon that has remained a theoretical conjecture for 90 years.