Beyond the Instant: How Nuclei ‘Breathe’ During Reactions

Author: Denis Avetisyan

New theoretical work reveals a surprisingly stable, long-lived intermediate state in nuclear reactions, challenging existing models of how atomic nuclei exchange particles.

The energies associated with the extended Gaussian-Core Model (eGCM) align with the diagonal matrix elements representing the overlap between basis states <span class="katex-eq" data-katex-display="false">\langle\overline{\Phi}\_{k}|\hat{\rm H}|\overline{\Phi}\_{k}\rangle</span>, as indicated by the correspondence between the blue line and the red dashed line; further analysis of the eigenvalues <span class="katex-eq" data-katex-display="false">\nu\_{k}</span> of these norm overlaps confirms a non-zero contribution beyond the trivial case <span class="katex-eq" data-katex-display="false">\nu\_{k} \neq 0</span>. — The energies associated with the extended Gaussian-Core Model (eGCM) align with the diagonal matrix elements representing the overlap between basis states $\langle\overline{\Phi}\_{k}|\hat{\rm H}|\overline{\Phi}\_{k}\rangle$ , as indicated by the correspondence between the blue line and the red dashed line; further analysis of the eigenvalues $\nu\_{k}$ of these norm overlaps confirms a non-zero contribution beyond the trivial case $\nu\_{k} \neq 0$ .

The Enhanced Generator Coordinate Method demonstrates the formation of a compound nucleus in multi-nucleon transfer reactions, aligning with Random Matrix Theory and offering a new perspective on nuclear dynamics.

Despite the longstanding challenge of microscopically describing compound nucleus formation, this work-Multi-Nucleon Transfer Reactions and the Creation and the Evolution of the Compound Nucleus-introduces the enhanced Generator Coordinate Method (eGCM) to model multi-nucleon transfer reactions near the Coulomb barrier. eGCM demonstrates the formation of a long-lived ‘compound nucleus’ exhibiting qualitative differences from both time-dependent Hartree-Fock and standard Generator Coordinate Method calculations. This novel approach successfully predicts compound nucleus formation cross sections, suggesting a nuclear analogue to optical molasses and aligning with predictions from Random Matrix Theory. Will this framework provide a more complete understanding of nuclear reaction dynamics and the evolution of complex nuclear systems?

The Emergence of Complexity in Nuclear Collisions

The simulation of multi-nucleon transfer reactions-where several protons and neutrons are exchanged between colliding nuclei-is fundamentally hampered by the sheer complexity of the “many-body problem.” Unlike simpler reactions involving just a few particles, these transfers involve the simultaneous interactions of dozens, or even hundreds, of nucleons. Each nucleon experiences a complex potential created by all the others, and accurately describing these interwoven forces requires solving the Schrödinger equation for a system with an astronomically large number of degrees of freedom. This quickly exceeds the capabilities of even the most powerful supercomputers, forcing physicists to rely on approximations that, while necessary, inevitably introduce uncertainties into the predictions. Consequently, understanding the detailed mechanisms of these reactions – crucial for fields like nuclear astrophysics and the production of rare isotopes – remains a significant computational challenge, demanding ongoing development of novel theoretical approaches and algorithms.

The intricacies of multi-nucleon transfer reactions are often obscured by the limitations of conventional modeling techniques. These methods typically treat the compound nucleus – the temporary amalgamation of colliding nuclei – as a system where nucleons behave largely independently. However, this simplification fails to capture the complex, collective behavior arising from the strong nuclear force and the quantum mechanical correlations between individual nucleons. Consequently, predictions regarding reaction probabilities and the distribution of transferred nucleons diverge significantly from experimental observations. This inability to accurately describe the dynamic interplay within the compound nucleus hinders a comprehensive understanding of nuclear structure and the processes governing the creation of new isotopes, necessitating the development of more sophisticated theoretical frameworks that account for these crucial correlations.

The fidelity of nuclear reaction predictions hinges on a complete description of nucleon correlations – the intricate ways in which protons and neutrons move in concert within the nucleus. These correlations arise from the strong nuclear force and significantly impact reaction pathways, influencing everything from the production of new isotopes to the energy distribution of emitted particles. However, modeling these many-body effects presents a formidable computational challenge. Unlike simpler, independent-particle models, accurately accounting for nucleon correlations requires considering the simultaneous interactions of numerous particles, leading to exponential scaling in computational cost. Consequently, even with powerful supercomputers, approximations are often necessary, limiting the precision of reaction predictions and underscoring the need for innovative theoretical approaches and efficient algorithms to tame the complexity of the nuclear many-body problem.

Time-dependent Hartree-Fock simulations considered impact parameters <span class="katex-eq" data-katex-display="false"> \textbf{b} </span> within a range of 5-6 fm, as illustrated by representative density profiles at <span class="katex-eq" data-katex-display="false"> t=420 </span> fm/c and at the simulations' final time, with particle number projection performed within the white circles to analyze late-time behavior. — Time-dependent Hartree-Fock simulations considered impact parameters $\textbf{b}$ within a range of 5-6 fm, as illustrated by representative density profiles at $t=420$ fm/c and at the simulations’ final time, with particle number projection performed within the white circles to analyze late-time behavior.

A First Step Towards Accuracy: The Generator Coordinate Method

The Generator Coordinate Method (GCM) addresses the many-body Schrödinger equation by expressing the total wavefunction, Ψ, as a superposition of functions generated by a set of collective coordinates. These coordinates, denoted as α, parameterize the degrees of freedom relevant to the system’s structure. The wavefunction is thus expanded as $\Psi = \in t F(\alpha) | \Phi(\alpha) \rangle d\alpha$ , where $| \Phi(\alpha) \rangle$ represents a generator wavefunction dependent on the collective coordinates and $F(\alpha)$ is a function determining the superposition weights. This approach effectively reduces the complexity of solving the full Schrödinger equation by focusing on collective motion and incorporating correlations through the generator wavefunctions, allowing for calculations of nuclear structure and reaction properties.

The Generator Coordinate Method (GCM) is utilized to describe the structure of compound nuclei formed in nuclear reactions by expanding the many-body wavefunction as a superposition of generator states. While conceptually valuable, standard GCM implementations face limitations stemming from the truncation of the generator space and the Gaussian overlap approximation used in wavefunction normalization. These approximations can lead to inaccuracies in calculated energies and transition probabilities, particularly for systems with significant correlations or complex configurations. Furthermore, the computational cost of GCM scales rapidly with the number of generator states included, restricting the size and complexity of systems that can be realistically modeled with standard implementations.

Advancing the accuracy and computational efficiency of the Generator Coordinate Method (GCM) is essential for modeling complex nuclear reactions, particularly those involving numerous particles or extended nuclear structures. Current limitations in GCM calculations-stemming from the need for extensive basis sets and the computational cost of propagating many-body wavefunctions-restrict its application to relatively simple systems or limited regions of the nuclear chart. Improvements in areas such as optimized basis generation, efficient many-body force calculations, and the development of algorithms for truncating the generator space are critical for extending the reach of GCM to heavier nuclei and more realistic reaction scenarios, including those relevant to astrophysical environments and nuclear technology. Addressing these challenges will enable more reliable predictions of reaction cross-sections, decay properties, and the formation of exotic nuclear species.

The overlap of the Information Propagation Rate (IPR) between the Norm and Hamiltonian effective generalized coordinate mapping reveals a maximum possible IPR, indicated by the dashed line.

Refining Accuracy and Efficiency: The Enhanced Generator Coordinate Method

The Enhanced Generator Coordinate Method (eGCM) builds upon the Generator Coordinate Method by utilizing specifically designed basis sets to represent the many-body wavefunctions of complex nuclear systems. These basis sets are constructed from Bloch wavefunctions, which are particularly suited for describing collective motion, and Wannier wavefunctions, which provide a localized description of single-particle states. The combination of these wavefunction types allows eGCM to more effectively capture the correlations between nucleons and accurately represent the complex shapes and configurations observed in nuclear reactions and structures. This approach contrasts with traditional GCM methods that may employ simpler, less flexible basis sets, potentially limiting their ability to accurately model these systems.

The Enhanced Generator Coordinate Method (eGCM) demonstrates improved accuracy in modeling complex nuclear systems through the utilization of advanced basis sets. Simulations utilizing eGCM have yielded a compound nucleus formation probability of 0.34, representing a quantifiable measure of its effectiveness in predicting reaction outcomes. This probability is determined through calculations based on the generated basis sets and reflects the likelihood of successful compound nucleus formation within the modeled system, providing a benchmark for evaluating the method’s predictive power and validating its application to increasingly complex nuclear scenarios.

The Enhanced Generator Coordinate Method (eGCM) incorporates the impact parameter to refine the modeling of reaction dynamics. Simulations utilize an impact parameter range of 5-6 fm, allowing for a more detailed analysis of the relative angular momentum between colliding nuclei. This parameter is crucial because it directly influences the fusion probability and the excitation function of the resulting compound nucleus. By varying the impact parameter within this range, the eGCM framework can account for different collision geometries and their corresponding contributions to the overall reaction cross-section, resulting in a more accurate representation of the nuclear reaction process.

The Enhanced Generator Coordinate Method (eGCM) necessitates the construction of accurate basis sets to represent the many-body wavefunctions of complex nuclear systems. These basis sets are frequently generated using Time-Dependent Hartree-Fock (TDHF) trajectories, which provide initial configurations for subsequent refinement. A defining characteristic of the eGCM approach is the scale of the GCM basis employed; current implementations utilize a basis set size of 39,630, enabling a detailed and comprehensive description of the nuclear many-body problem, despite the associated computational demands.

The significant divergence between Time-Dependent Hartree-Fock (TDHF) and extended Generator Coordinate Method (eGCM) probabilities for heavy fragments is demonstrated by comparing their respective probability distributions <span class="katex-eq" data-katex-display="false">P_H(N)</span>, <span class="katex-eq" data-katex-display="false">P_H(Z)</span> (red lines) with those of lighter fragments <span class="katex-eq" data-katex-display="false">P_L(N)</span>, <span class="katex-eq" data-katex-display="false">P_L(Z)</span> (solid lines), and the overall probabilities are represented by solid and dotted black and green lines. — The significant divergence between Time-Dependent Hartree-Fock (TDHF) and extended Generator Coordinate Method (eGCM) probabilities for heavy fragments is demonstrated by comparing their respective probability distributions $P_H(N)$ , $P_H(Z)$ (red lines) with those of lighter fragments $P_L(N)$ , $P_L(Z)$ (solid lines), and the overall probabilities are represented by solid and dotted black and green lines.

Validating with Random Matrix Theory: A Statistical Benchmark

Random Matrix Theory (RMT) offers a statistical approach to analyzing the properties of systems with a vast number of interacting quantum states, such as the compound nucleus formed in heavy-ion collisions. Unlike traditional quantum mechanical methods that attempt to solve for individual energy levels, RMT focuses on the statistical distribution of these levels, assuming that the detailed microscopic interactions are not fully known or are too complex to model directly. This approach relies on the universality observed in the energy level statistics of diverse systems – from nuclear physics to quantum chaos and even certain financial models – suggesting that the underlying statistical behavior is independent of the specific details of the system. By characterizing the energy level distribution using ensembles like the Gaussian Orthogonal Ensemble (GOE), RMT provides a benchmark against which the predictions of more detailed models, such as the extended Gaussian Core Model (eGCM), can be rigorously tested.

The extended Gaussian Orthogonal Ensemble (GOE) provides a statistically derived benchmark against which to validate the accuracy of energy level distributions predicted by the extended Green’s function coupled channel method (eGCM). Specifically, the GOE predicts the distribution of energy level spacings in complex quantum systems exhibiting chaotic behavior; by comparing the statistical properties of eGCM-calculated energy levels to those predicted by the GOE, researchers can assess whether eGCM accurately captures the underlying quantum dynamics of the compound nucleus. This rigorous test involves analyzing parameters such as the average level spacing and the variance of spacing distribution, allowing for quantitative confirmation of the method’s reliability in reproducing expected energy level characteristics.

The Wigner-Dyson surmise, a statistical prediction derived from Random Matrix Theory (RMT), details the expected distribution of energy levels in complex quantum systems. Specifically, the Gaussian Orthogonal Ensemble (GOE) predicts a characteristic form for the spacing between adjacent energy levels, described by a formula involving $\frac{2\sqrt{\pi}}{\Gamma(1/2)}$ . By comparing the observed energy level spacing from extended Gaussian Core Model (eGCM) calculations with this surmise, the accuracy of eGCM in representing the quantum chaotic behavior of the compound nucleus can be rigorously assessed. Significant deviation from the predicted Wigner-Dyson distribution would indicate inadequacies in the eGCM’s treatment of the system’s complexity.

Analysis of energy level spacing derived from the extended Gaussian Core Model (eGCM) demonstrates a strong correlation with predictions from Random Matrix Theory (RMT). Specifically, eGCM simulations yield an average level spacing of 4 keV with a standard deviation of 0.5273. This value is remarkably close to the theoretical Gaussian Orthogonal Ensemble (GOE) standard deviation of 0.5227, indicating that the statistical properties of the energy levels generated by eGCM accurately reflect the characteristics expected of complex quantum systems governed by RMT. This agreement serves as a key validation point, bolstering confidence in the reliability and predictive power of the eGCM method for describing the compound nucleus.

Looking Ahead: Expanding the Horizon of Nuclear Simulations

The advent of the extended Gaussian-Core Model (eGCM) represents a significant leap forward in the simulation of complex nuclear reactions, particularly when coupled with the power of Time-Dependent Density Functional Theory (TDDFT). This combination allows researchers to probe the dynamic behavior of nuclei with unprecedented detail, moving beyond static approximations. The inclusion of advanced functionals, such as SeaLL1, within the TDDFT framework further refines these simulations by providing a more accurate description of the strong nuclear force and its influence on nuclear processes. This methodology isn’t limited to simply modeling known reactions; it facilitates the investigation of previously inaccessible phenomena and opens doors to predicting the outcomes of nuclear encounters under extreme conditions, offering insights into the creation of elements in stars and the fundamental nature of nuclear matter itself.

The convergence of enhanced methods like the energy-constrained Gaussian-constrained mean field (eGCM) and time-dependent density functional theory (TDDFT) presents an unprecedented opportunity to investigate the full spectrum of nuclear behavior. From detailing the intricate arrangements of protons and neutrons within individual nuclei – revealing the fundamental architecture of nuclear structure – these simulations extend to the extreme conditions found in astrophysical environments. Researchers can now model the creation of elements within stars, the explosive events of supernovae, and the dynamics of neutron star mergers, bridging the gap between laboratory experiments and cosmic phenomena. This capability allows for a deeper understanding of nucleosynthesis, the process by which heavier elements are formed, and provides crucial insights into the origin of elements found throughout the universe, effectively linking the microcosm of the nucleus to the macrocosm of astrophysical events.

Continued development of the extended Gaussian-based Configuration Mixing (eGCM) framework promises increasingly precise insights into the intricacies of the nuclear realm. Current research focuses on constructing more sophisticated basis sets, capable of representing the complex wavefunctions of atomic nuclei with greater fidelity, and integrating relativistic effects – crucial for accurately describing the behavior of nucleons at high energies and in extreme environments. These enhancements will not only refine existing simulations but also unlock the ability to explore previously inaccessible nuclear phenomena, pushing the boundaries of nuclear physics and astrophysics by providing a more complete and nuanced understanding of nuclear structure and reactions. The pursuit of these advancements will ultimately lead to more reliable predictions of nuclear properties and behaviors, fostering breakthroughs in fields ranging from fundamental nuclear theory to the modeling of stellar evolution and nucleosynthesis.

The culmination of this research delivers a pathway towards substantially more precise calculations of nuclear reaction rates and cross sections, data vital across a spectrum of scientific and technological fields. Accurate predictions in these areas underpin advancements in nuclear astrophysics – allowing for more refined modeling of stellar evolution and nucleosynthesis – and are essential for the development of nuclear technologies, including energy production and medical isotope creation. Remarkably, this improved predictive capability was achieved within a computationally manageable timeframe of 2,482 femtometers per c, demonstrating the efficiency of the implemented methods and suggesting scalability for even more complex systems. This benchmark establishes a robust foundation for future investigations demanding high-fidelity simulations of nuclear processes.

Analysis of the Information Potential Ratio (IPR) using expansion coefficients <span class="katex-eq" data-katex-display="false">a_k(t)</span> (blue/red) and <span class="katex-eq" data-katex-display="false">d(\alpha|t)</span> (green/purple) reveals the timing of nuclear contact (dashed line) and separation (gray band) in Time-Dependent Hartree-Fock (TDHF) simulations, as corroborated by the inset displaying the squared magnitudes of these coefficients alongside <span class="katex-eq" data-katex-display="false">c_n(t)</span>. — Analysis of the Information Potential Ratio (IPR) using expansion coefficients $a_k(t)$ (blue/red) and $d(\alpha|t)$ (green/purple) reveals the timing of nuclear contact (dashed line) and separation (gray band) in Time-Dependent Hartree-Fock (TDHF) simulations, as corroborated by the inset displaying the squared magnitudes of these coefficients alongside $c_n(t)$ .

The study illuminates how complex nuclear landscapes emerge not from imposed design, but from the interplay of local nucleon exchanges. The creation of a long-lived compound nucleus, demonstrated via the Enhanced Generator Coordinate Method, isn’t a feat of engineered stability, but a consequence of the system finding its own robust configuration. This resonates deeply with the notion that order doesn’t require architects; it arises from the rules governing the interactions themselves. As Richard Feynman once observed, “The best way to understand something is to build it.” Here, the ‘building’ isn’t a conscious act, but the system’s natural exploration of configurations, revealing robustness inherent in the dynamics of multi-nucleon transfer, not pre-planned within it.

The Landscape Ahead

The demonstration of prolonged compound nucleus formation via multi-nucleon transfer, as revealed through the Enhanced Generator Coordinate Method, subtly shifts the understanding of nuclear reaction dynamics. It suggests that the transition from direct reaction to full equilibration is not the abrupt event often depicted, but a more gradual emergence – a coastline blurring between land and sea. The alignment with Random Matrix Theory is less a validation of a particular model, and more a recognition that complexity, at a certain threshold, behaves according to statistical principles; the forest evolves without a forester, yet follows rules of light and water.

Remaining are questions not of how this long-lived compound nucleus forms, but of its influence. Does this prolonged existence meaningfully alter reaction pathways, favoring specific decay channels? Current explorations largely treat the compound nucleus as a transient state; future work might examine its potential to act as a catalyst, shaping the distribution of reaction products. Furthermore, the computational demands of eGCM present a practical limitation; simplifying assumptions, while necessary, introduce the risk of obscuring subtle but significant effects.

The field now faces a choice: to refine existing theoretical frameworks, seeking ever-greater precision in describing known phenomena, or to embrace the inherent unpredictability of complex systems. Order is the result of local interactions, not directives. Perhaps the most fruitful path lies in focusing not on controlling the reaction, but on understanding the rules governing its spontaneous organization.

Original article: https://arxiv.org/pdf/2604.21845.pdf

Contact the author: https://www.linkedin.com/in/avetisyan/