Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TDHF”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

The TDHF code Sky3D version 1.1

The nuclear mean-field model based on Skyrme forces or related density functionals has found widespread application to the description of nuclear ground states, collective vibrational excitations, and heavy-ion collisions. The code Sky3D solves the static or dynamic equations on a three-dimensional Cartesian mesh with isolated or periodic boundary conditions and no further symmetry assumptions. Pairing can be included in the BCS approximation for the static case. The code is implemented with a view to allow easy modifications for including additional physics or special analysis of the results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Density constrained TDHF [Book Chapter]

In this manuscript we provide an outline of the numerical methods used in implementing the density constrained time-dependent Hartree-Fock (DC-TDHF) method and provide a few examples of its application to nuclear fusion. In this approach, dynamic microscopic calculations are carried out on a three-dimensional lattice and there are no adjustable parameters, the only input is the Skyrme effective NN interaction. After a review of the DC-TDHF theory and the numerical methods, we present results for heavy-ion potentials V (R), coordinate-dependent mass parameters M(R), and precompound excitation energies E*(R) for a variety of heavy-ion reactions. Using fusion barrier penetrabilities, we calculate total fusion cross sections σ(E c.m. ) for reactions between both stable and neutron-rich nuclei. We also determine capture cross sections for hot fusion reactions leading to the formation of superheavy elements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Using TDHF to Study Quasifission Dynamics [Book Chapter]

We show that the microscopic TDHF approach provides an important tool to shed some light on the nuclear dynamics leading to the formation of superheavy elements. In particular, we discuss studying quasifission dynamics and calculating ingredients for compound nucleus formation probability calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

TDHF Investigations of the U+U Quasifission Process

The use of actinide collisions have been suggested as a way to produce neutron-rich isotopes of high-Z nuclei. The collision dynamics of these reactions can be studied using unrestricted time-dependent Hartree–Fock (TDHF) calculations. Here, we report on the recent studies of quasifission for the 238 U+ 238 U system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

The TDHF code Sky3D version 1.2

The Sky3D code has been widely used to describe nuclear ground states, collective vibrational excitations, and heavy-ion collisions. The approach is based on Skyrme forces or related energy density functionals. The static and dynamic equations are solved on a three-dimensional grid, and pairing is been implemented in the BCS approximation. This updated version of the code aims to facilitate the calculation of nuclear strength functions in the regime of linear response theory, while retaining all existing functionality and use cases. The strength functions are benchmarked against available RPA codes, and the user has the freedom of choice when selecting the nature of external excitation (from monopole to hexadecapole and more). Some utility programs are also provided that calculate the strength function from the time-dependent output of the dynamic calculations of the Sky3D code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Absence of hindrance in a microscopic 12 C + 12 C fusion study

Studies of low-energy fusion of light nuclei are important in astrophysical modeling, with small variations in reaction rates having a large impact on nucleosynthesis yields. Due to the lack of experimental data at astrophysical energies, extrapolation and microscopic methods are needed to model fusion probabilities. To investigate deep sub-barrier 12 C+ 12 C fusion cross sections and establish trends for the S factor. Microscopic methods based on static Hartree-Fock and time-dependent Hartree-Fock (TDHF) mean-field theory are used to obtain 12 C+ 12 C ion-ion fusion potentials. Fusion cross sections and astrophysical S factors are then calculated using the incoming wave boundary condition method. Both density-constrained frozen Hartree-Fock (DCFHF) and density-constrained TDHF (DC-TDHF) predict a rising S factor at low energies, with DC-TDHF predicting a slight damping in the deep sub-barrier region (≈1 MeV). Comparison between DC-TDHF calculations and maximum experimental cross sections in the resonance peaks are good. However, the discrepancy in experimental low-energy results inhibits interpretation of the trend. Using the fully microscopic DCFHF and DC-TDHF methods, no S factor maximum is observed in the 12 C+ 12 C fusion reaction. In addition, no extreme sub-barrier hindrance is predicted at low energies. Furthermore, the development of a microscopic theory of fusion including resonance effects, as well as further experiments at lower energies must be done before the deep sub-barrier behavior of the reaction can be established.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

The Effective Fragment Molecular Orbital Method: Achieving High Scalability and Accuracy for Large Systems

The effective fragment molecular orbital (EFMO) method has been developed to predict the total energy of a very large molecular system accurately (with respect to the underlying quantum mechanical method) and efficiently by taking advantage of the locality of strong chemical interactions and employing a two-level hierarchical parallelism. The accuracy of the EFMO method is partly attributed to the accurate and robust intermolecular interaction prediction between distant fragments, in particular, the many-body polarization and dispersion effects, which require the generation of static and dynamic polarizability tensors by solving the coupled perturbed Hartree–Fock (CPHF) and time-dependent HF (TDHF) equations, respectively. Solving the CPHF and TDHF equations is the main EFMO computational bottleneck due to the inefficient (serial) and I/O-intensive implementation of the CPHF and TDHF solvers. In this work, the efficiency and scalability of the EFMO method are significantly improved with a new CPU memory-based implementation for solving the CPHF and TDHF equations that are parallelized by either message passing interface (MPI) or hybrid MPI/OpenMP. Here, the accuracy of the EFMO method is demonstrated for both covalently bonded systems and noncovalently bound molecular clusters by systematically examining the effects of basis sets and a key distance-related cutoff parameter, R cut . R cut determines whether a fragment pair (dimer) is treated by the chosen ab initio method or calculated using the effective fragment potential (EFP) method (separated dimers). Decreasing the value of Rcut increases the number of separated (EFP) dimers, thereby decreasing the computational effort. It is demonstrated that excellent accuracy (<1 kcal/mol error per fragment) can be achieved when using a sufficiently large basis set with diffuse functions coupled with a small R cut value. With the new parallel implementation, the total EFMO wall time is substantially reduced, especially with a high number of MPI ranks. Given a sufficient workload, nearly ideal strong scaling is achieved for the CPHF and TDHF parts of the calculation. For the first time, EFMO calculations with the inclusion of long-range polarization and dispersion interactions on a hydrated mesoporous silica nanoparticle with explicit water solvent molecules (more than 15k atoms) are achieved on a massively parallel supercomputer using nearly 1000 physical nodes. In addition, EFMO calculations on the carbinolamine formation step of an amine-catalyzed aldol reaction at the nanoscale with explicit solvent effects are presented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Heavy-ion collisions and fission dynamics with the time-dependent Hartree–Fock theory and its extensions

Microscopic methods and tools to describe nuclear dynamics have considerably been improved in the past few years. They are based on the time-dependent Hartree–Fock (TDHF) theory and its extensions to include pairing correlations and quantum fluctuations. The TDHF theory is the lowest level of approximation of a range of methods to solve the quantum many-body problem, showing its universality to describe many-fermion dynamics at the mean-field level. The range of applications of TDHF to describe realistic systems allowing for detailed comparisons with experiment has considerably increased. For instance, TDHF is now commonly used to investigate fusion, multi-nucleon transfer and quasi-fission reactions. Thanks to the inclusion of pairing correlations, it has also recently led to breakthroughs in our description of the saddle to scission evolution, and, in particular, the non-adiabatic effects near scission. Beyond mean-field approaches such as the time-dependent random-phase approximation (TDRPA) and stochastic mean-field methods have reached the point where they can be used for realistic applications. We review recent progresses in both techniques and applications to heavy-ion collision and fission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Hybrid algorithm for the time-dependent Hartree–Fock method using the Yang–Baxter equation on quantum computers *

Abstract The time-dependent Hartree–Fock (TDHF) method is an approach to simulate the mean field dynamics of electrons within the assumption that the electrons move independently in their self-consistent average field and within the space of single Slater determinants. One of the major advantages of performing time dynamics within Hartree–Fock theory is the free fermionic nature of the problem, which makes TDHF classically simulatable in polynomial time. Here, we present a hybrid TDHF implementation for quantum computers. This quantum circuit grows with time; but with our recent work on circuit compression via the Yang–Baxter equation (YBE), the resulting circuit is constant depth. This study provides a new way to simulate TDHF with the aid of a quantum device as well as provides a new direction for the application of YBE symmetry in quantum chemistry simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Quantal diffusion approach for multinucleon transfer processes in the 58,64 Ni + 208 Pb reactions: Toward the production of unknown neutron-rich nuclei

In recent years, substantial efforts have been made for the study of multinucleon transfer reactions at energies around the Coulomb barrier both experimentally and theoretically, aiming at the production of unknown neutron-rich heavy nuclei. It is crucial to provide reliable theoretical predictions based on microscopic theories with sufficient predictive power. Purpose: This article aims to clarify the applicability of the quantal diffusion approach based on the stochastic mean-field (SMF) theory for multinucleon transfer processes. Isotope production cross sections are evaluated for the reactions of 64 Ni + 208 Pb at E c.m. = 268 MeV and 58 Ni + 208 Pb at E c.m. = 270 MeV and are compared with available experimental data. Methods: Three-dimensional time-dependent Hartree-Fock (TDHF) calculations are carried out for a range of initial orbital angular momenta with Skyrme SLy4d functional. Quantal diffusion equations, derived based on the SMF theory, for variances and covariance of neutron and proton numbers of reaction products are solved, with microscopic drift and diffusion coefficients obtained from time evolution of occupied single-particle orbitals in TDHF. Secondary de-excitation processes, both particle evaporation and fission, are simulated by a statistical compound-nucleus de-excitation model, gemini++. Dynamics of a fast isospin equilibration process followed by a slow drift toward the mass symmetry are commonly observed, as expected. Various reaction outcomes are evaluated, including average mass and charge numbers of reaction products, total kinetic energy loss (TKEL), scattering angle, contact time, and production cross sections for primary and secondary products. By comparing with the experimental data, we find that SMF and TDHF quantitatively reproduce experimental data for few-nucleon-transfer channels around the average values. In contrast, for many-nucleon-transfer channels, we find that the SMF approach provides much better description of the experimentally measured isotopic distributions. The results underline the importance of beyond-mean-field effects, especially one-body (mean-field) fluctuations and correlations, in describing multinucleon transfer processes. Moreover, through a combined analysis of SMF with a statistical model, gemini++, we find a significant contribution of transfer-induced fission, which is consistent with the experimental observation. In some cases, the SMF approach overestimates the isotopic width, requiring further improvements of the theoretical description. Possible ways to improve the description are discussed. The SMF approach is designed to describe the quantum many-body problem according to an ensemble of mean-field trajectories, taking into account part of many-body correlations in the description. As it requires feasible computational costs comparable to the ordinary TDHF approach, together with further model improvements, it will be a promising tool in the search for optimal reaction conditions to produce yet-unknown neutron-rich heavy nuclei through the multinucleon transfer reaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Size-dependent errors in real-time electron density propagation

Real-time (RT) electron density propagation with time-dependent density functional theory (TDDFT) or Hartree–Fock (TDHF) is one of the most popular methods to model the charge transfer in molecules and materials. However, both RT-TDHF and RT-TDDFT within the adiabatic approximation are known to produce inaccurate evolution of the electron density away from the ground state in model systems, leading to large errors in charge transfer and erroneous shifting of peaks in absorption spectra. Given the poor performance of these methods with small model systems and the widespread use of the methods with larger molecular and material systems, here we bridge the gap in our understanding of these methods and examine the size-dependence of errors in RT density propagation. We analyze the performance of RT density propagation for systems of increasing size during the application of a continuous resonant field to induce Rabi-like oscillations, during charge-transfer dynamics, and for peak shifting in simulated absorption spectra. We find that the errors in the electron dynamics are indeed size dependent for these phenomena, with the largest system producing the results most aligned with those expected from linear response theory. The results suggest that although the RT-TDHF and RT-TDDFT methods may produce severe errors for model systems, the errors in charge transfer and resonantly driven electron dynamics may be much less significant for more realistic, large-scale molecules and materials.

Chemistry↗
DOE OSTI2024

Scalable learning of potentials to predict time-dependent Hartree–Fock dynamics

We propose a framework to learn the time-dependent Hartree–Fock (TDHF) inter-electronic potential of a molecule from its electron density dynamics. Although the entire TDHF Hamiltonian, including the inter-electronic potential, can be computed from first principles, we use this problem as a testbed to develop strategies that can be applied to learn a priori unknown terms that arise in other methods/approaches to quantum dynamics, e.g., emerging problems such as learning exchange–correlation potentials for time-dependent density functional theory. We develop, train, and test three models of the TDHF inter-electronic potential, each parameterized by a four-index tensor of size up to 60 × 60 × 60 × 60. Two of the models preserve Hermitian symmetry, while one model preserves an eight-fold permutation symmetry that implies Hermitian symmetry. Across seven different molecular systems, we find that accounting for the deeper eight-fold symmetry leads to the best-performing model across three metrics: training efficiency, test set predictive power, and direct comparison of true and learned inter-electronic potentials. All three models, when trained on ensembles of field-free trajectories, generate accurate electron dynamics predictions even in a field-on regime that lies outside the training set. To enable our models to scale to large molecular systems, we derive expressions for Jacobian-vector products that enable iterative, matrix-free training.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Microscopic predictions for the production of neutron-rich nuclei in the reaction Yb 176 + Yb 176

Background: Production of neutron-rich nuclei is of vital importance to both understanding nuclear structure far from stability and to informing astrophysical models of the rapid neutron capture process (r-process). Multinucleon transfer (MNT) in heavy-ion collisions offers a possibility to produce neutron-rich nuclei far from stability. Purpose: The 176 Yb + 176 Yb reaction has been suggested as a potential candidate to explore the neutron-rich region surrounding the principal fragments. The current study has been conducted with the goal of providing guidance for future experiments wishing to study this (or similar) system. Methods: Time-dependent Hartree-Fock (TDHF) and its time-dependent random-phase approximation (TDRPA) extension are used to examine both scattering and MNT characteristics in 176 Yb + 176 Yb. TDRPA calculations are performed to compute fluctuations and correlations of the neutron and proton numbers, allowing for estimates of primary fragment production probabilities. Results: Both scattering results from TDHF and transfer results from the TDRPA are presented for different energies, orvientations, and impact parameters. In addition to fragment composition, scattering angles and total kinetic energies, as well as correlations between these observables are presented. Conclusions: 176 Yb + 176 Yb appears to be an interesting probe for the midmass neutron-rich region of the chart of nuclides. Lastly, the predictions of both TDHF and TDRPA are speculative, and will benefit from future experimental results to test the validity of this approach to studying MNT in heavy, symmetric collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Production of neutron-rich heavy nuclei in deep-inelastic 208 Pb + 208 Pb collisions within the stochastic mean-field theory

In deep-inelastic collisions of heavy nuclei, reaction products with a wide range of mass and charge are produced. Such collisions been considered as a possible way to produce superheavy nuclei, as an alternative to fusion reactions. To provide reliable theoretical predictions, it is desired to develop microscopic approaches that correctly and accurately describe nucleon transfer processes in dissipative collisions of heavy nuclei. The purpose of the present work is (1) to investigate the mechanism of nucleon transfers in dissipative collisions of two heavy nuclei, and (2) to explore possible pathways to produce neutron-rich heavy nuclei, through detailed theoretical analyses of fluctuations and correlations in nucleon transfers in 208 Pb + 208 Pb reactions. Three-dimensional time-dependent Hartree-Fock (TDHF) calculations are performed for the collisions of 208 Pb + 208 Pb at 𝐸 c.m. = 832, 936, and 1040 MeV, using the Skyrme SLy4d energy density functional. To calculate fluctuations and correlations in nucleon transfers, we employ the stochastic mean-field (SMF) theory, and the results are compared with another theoretical framework currently available, the time-dependent random phase approximation (TDRPA). Primary and secondary production cross sections are calculated with the SMF theory combined with a statistical model, GEMINI ++ . Using information of nucleon flow across a neck of colliding nuclei in TDHF calculations, we solve quantal diffusion equations for fluctuations and correlations in nucleon transfers based on the SMF theory. From the SMF calculations, we obtain the time evolution of diffusion coefficients as well as fluctuations and correlations in nucleon transfers for a range of initial orbital angular momenta. We compare the results of the SMF calculations with those of TDRPA, showing that TDRPA tends to predict substantially larger fluctuations and correlations in strongly damped collisions of heavy nuclei, which exhibit complex initial angular momentum dependence, while the SMF results provide almost constant (stable) values. Using the obtained fluctuations and correlations, we calculate primary and secondary production cross sections for the 208 Pb + 208 Pb collisions. From the results, we find that both lighter and heavier reaction products as compared to 208 Pb are produced for a wide region in the 𝑁−𝑍 plane as primary products, thanks to the quantal diffusion mechanism in the dissipative collisions. However, we show that cross sections for production of heavy nuclei with 𝑍 ≳ 90 or 𝑁 ≳ 135 are washed out due to secondary particle evaporation and/or fission processes. On the other hand, we find that there remain sizable cross sections for production of neutron-rich nuclei along 𝑁 = 126 with 𝑍< 82, even after secondary disintegration processes. We demonstrate that the secondary production cross sections depend weakly on incident energies, but lower (higher) energy is slightly preferred for production of nuclei with smaller (larger) atomic numbers as compared to 𝑍 = 82. Here, based on the microscopic SMF calculations, it has been shown that deep-inelastic collisions of heavy nuclei, such as 208 Pb + 208 Pb examined in this study, can be a promising means to produce neutron-rich heavy nuclei along 𝑁=126. Discrepancies between the SMF and TDRPA approaches are left unsolved for future investigations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Microscopic study of the fusion reactions 40,48 Ca+ 78 Ni and the effect of the tensor force

In this work, we provide a microscopic description of the fusion reactions between 40,48 Ca and 78 Ni. The internuclear potentials are obtained using the density-constrained (DC) time-dependent Hartree-Fock (TDHF) approach and fusion cross sections are calculated via the incoming wave boundary condition method. By performing DC-TDHF calculations at several selected incident energies, the internuclear potentials for both systems are obtained and the energy-dependence of fusion barrier are revealed. The influence of tensor force on internuclear potentials of 48 Ca+ 78 Ni is more obvious than those of 40 Ca+ 78 Ni. By comparing the calculated fusion cross sections between 40 Ca+ 78 Ni and 48 Ca+ 78 Ni, an interesting enhancement of subbarrier fusion cross sections for the former system is found, which can be explained by the narrow width of internuclear potential for 40 Ca+ 78 Ni while the barrier heights and positions are very close to each other. The tensor force suppresses the subbarrier fusion cross sections of both two systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Description of the multinucleon transfer mechanism for Ca 48 + Pu 244 and Kr 86 + Pt 198 reactions in a quantal transport approach

Multinucleon transfer (MNT) reactions involving heavy projectile and target combinations stand as a promising method for synthesizing new neutron-rich exotic nuclei, which may not be possible using hot or cold fusion reactions or fragmentation. Exploring the mechanisms behind MNT reactions is essential and it requires a comprehensive theoretical framework that can explain the physical observables in these reactions. This work aims to show that the quantal diffusion approach based on the stochastic mean-field (SMF) theory is capable of explaining the reaction dynamics observed in MNT reactions. Primary product mass distributions in 48 Ca + 244 Pu reaction at E c.m. = 203.2 MeV and 86 Kr + 198 Pt reaction at E c.m. = 324.2 MeV are calculated and compared with the available experimental data. In this work, we utilize the time-dependent Hartree-Fock (TDHF) calculations to analyze the mean-field reaction dynamics computationally in the reactions 48 Ca + 244 Pu and 86 Kr + 198 Pt for a broad range of initial angular momenta. Quantal transport description based on the SMF approach is used to calculate quantal diffusion coefficients and mass variances in 48 Ca + 244 Pu and 86 Kr + 198 Pt systems. The primary products arising from quasifission reactions are described by joint probability distribution in the SMF approach and those arising from fusion-fission are estimated by using the statistical deexcitation code gemini + +. Mean values of charge and mass numbers, scattering angles of the primary reaction products, and the total kinetic energies after the collision are calculated within the TDHF framework for a broad range of initial angular momenta. Throughout all the collisions, drift toward the mass symmetry and large mass dispersion associated with this drift are observed. Here, the calculated primary fragment and mass distributions using the SMF approach successfully explain experimental observations for the 48 Ca + 244 Pu and 86 Kr + 198 Pt systems. The primary mass distributions, mean values of binary products, and mass dispersions are determined and results are compared with the available experimental data. The observed agreement between the experimental data and SMF results highlights the effectiveness of the quantal diffusion mechanism based on the SMF approach, which does not include any adjustable parameters other than standard parameters of Skyrme energy density functional.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Beyond the Coulson–Fischer point: characterizing single excitation CI and TDDFT for excited states in single bond dissociations

Linear response time dependent density functional theory (TDDFT), which builds upon configuration interaction singles (CIS) and TD-Hartree–Fock (TDHF), is the most widely used class of excited state quantum chemistry methods and is often employed to study photochemical processes. This paper studies the behavior of the resulting excited state potential energy surfaces beyond the Coulson–Fischer (CF) point in single bond dissociations, when the optimal reference determinant is spin-polarized. Many excited states exhibit sharp kinks at the CF point, and connect to different dissociation limits via a zone of unphysical concave curvature. In particular, the unrestricted M S = 0 lowest triplet T 1 state changes character, and does not dissociate into ground state fragments. The unrestricted M S = ±1 T 1 CIS states better approximate the physical dissociation limit, but their degeneracy is broken beyond the CF point for most single bond dissociations. On the other hand, the M S = ±1 T 1 TDHF states reach the asymptote too soon, by merging with the ground state from the CF point onwards. Use of local exchange–correlation functionals causes M S = ±1 T 1 TDDFT states to resemble their unphysical M S = 0 counterpart. The 2 orbital, 2-electron model system of minimal basis H 2 is analytically treated to understand the origin of these issues, revealing that the lack of double excitations is at the root of these remarkable observations. In conclusion, the behavior of excited state surfaces is also numerically examined for species like H 2 , NH 3 , C 2 H 6 and LiH in extended basis sets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗