Search NASA⌕ Search

SEARCH · Search NASA

Results for “nuclear energy density functional”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Stabilization of Charge-Transfer Excited States in Biological Systems: A Computational Focus on the Special Pair in Photosystem II Reaction Centers

Charge-transfer (CT) excited states play an important role in many biological processes. However, many computational approaches often inadequately address the equilibration effects of nuclear and environmental degrees of freedom on these states. One prominent example of systems in which CT states are of utmost importance is reaction centers (RC) in photosystems. Here we use a multiscale approach combined with time-dependent density functional theory to explore the lowest CT excited state of the special pair P D1 –P D2 in the Photosystem II-RC of a cyanobacterium. We find that the nonequilibrium CT excited state resides near the Soret band, making an exciton the lowest-energy excited state. However, accounting for nuclear and state-specific dielectric equilibration along the CT potential energy surface (PES), the CT state P D1 – –P D2 + stabilizes energetically below the excitonic state. Importantly, this underscores the crucial role of state-specific solvation in mapping the PES of CT states, as demonstrated in a simplified dimer model.

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

Bayesian inference of fine features of the nuclear equation of state from future neutron star radius measurements to 0.1 km accuracy

To more precisely constrain the equation of state (EOS) of supradense neutron-rich nuclear matter, future high-precision x-ray and gravitational wave observatories are proposed to measure the radii of neutron stars (NSs) with an accuracy better than about 0.1 km. However, it remains unclear what particular aspects (other than the stiffness generally spoken of in the literature) of the EOS and to what precision they will be better constrained. In this work, within a Bayesian framework using a metamodel EOS for NSs, we infer the posterior probability distribution functions (PDFs) of incompressibility K 0 and skewness J 0 of symmetric nuclear matter (SNM) as well as the slope L, curvature K sym , and skewness J sym characterizing the density dependence of nuclear symmetry energy E sym ⁡(ρ), respectively, from mean values of NS radii consistent with existing observations and an expected accuracy Δ⁢R ranging from about 1.0 to 0.1 km. Here, we found that (1) the Δ⁢R has little effect on inferring the stiffness of SNM at suprasaturation densities, (2) smaller Δ⁢R reveals more accurately not only the PDFs but also pairwise correlations among parameters characterizing high-density E sym ⁡(ρ), (3) a double-peak feature of the PDF(K sym ) corresponding to the strong K sym – J sym and K sym – L anticorrelations is revealed when Δ⁢R is less than about 0.2 km, and the locations of the two peaks are sensitive to the maximum value of J sym reflecting the stiffness of E sym ⁡(ρ) above about 3 times the saturation density ρ 0 of SNM, and (4) the high-precision radius measurement for canonical NSs is more useful than that for massive ones for constraining the EOS of nucleonic matter around (2–3)⁢ρ 0 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Nuclear Computational Low Energy Initiative (NUCLEI)

The NUCLEI project, as defined by the scope of work, developed, implemented and run codes for large-scale computations of many topics in low-energy nuclear physics. Physics studied include the properties of nuclei and nuclear decays, nuclear structure and reactions, and the properties of nuclear matter. The computational techniques used include Quantum Monte Carlo, Configuration Interaction, Coupled Cluster, and Density Functional methods. The research program emphasized areas of high interest to current and possible future DOE nuclear physics facilities, including ATLAS and FRIB (nuclear structure and reactions, and nuclear astrophysics), TJNAF (neutron distributions in nuclei, few body systems, and electroweak processes), NIF (thermonuclear reactions), MAJORANA and FNPB (neutrinoless double-beta decay and physics beyond the Standard Model), and LANSCE (fission studies).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Diffusion Behavior of Oversized Fission Products in bcc Fe Cladding: A First-Principles Study

Fuel-Cladding Chemical Interaction (FCCI) poses significant challenges in nuclear reactors, where fission products from nuclear fuel interact with Fe-based cladding materials, potentially compromising their structural integrity. This study investigates the diffusion behavior of oversized fission products, Pr, Nd, Ce, and La, within bcc Fe cladding using density functional theory (DFT), nudged elastic band (NEB) method, and self-consistent mean field (SCMF) theory. Our results reveal significant long-range vacancy binding energies, particularly up to the 6th nearest neighbor, with La exhibiting the strongest binding affinity, followed by Nd, Ce, and Pr. The NEB calculations indicate significant high barriers for the dissociation of 1nn vacancy-solute pairs for all fission products. The tracer diffusion coefficients of these fission products was derived in Arrhenius form. The significant trapping effect of vacancies by a very dilute amount of fission products reduces vacancy mobility, leading to an oversaturation of point defects, void nucleation, and swelling. These are critical issues for irradiated cladding materials. The tracer diffusion coefficients indicate that Nd diffuses the fastest, followed by La, Ce, and Pr. This study provides essential insights for developing advanced cladding materials and design strategies to mitigate FCCI, ultimately enhancing nuclear reactor safety and performance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Weighted active space protocol for multireference machine-learned potentials

Multireference methods such as multiconfiguration pair-density functional theory accurately capture electronic correlation in systems with strong multiconfigurational character, but their cost precludes direct use in molecular dynamics. Combining these methods with machine-learned interatomic potentials (MLPs) can extend their reach. However, the sensitivity of multireference calculations to the choice of the active space complicates the consistent evaluation of energies and gradients across structurally diverse nuclear configurations. To overcome this limitation, we introduce the weighted active space protocol (WASP), a systematic approach to assign a consistent active space for a given system across uncorrelated configurations. By integrating WASP with MLPs and enhanced sampling techniques, we propose a data-efficient active learning cycle that enables the training of an MLP on multireference data. We demonstrated the approach on the TiC + -catalyzed C–H activation of methane, a reaction that poses challenges for Kohn–Sham density functional theory due to its significant multireference character. This framework enables accurate and efficient modeling of catalytic dynamics, establishing a paradigm for simulating complex reactive processes beyond the limits of conventional electronic-structure methods.

enhanced sampling↗
DOE OSTI2023

Nuclear–Electronic Orbital QM/MM Approach: Geometry Optimizations and Molecular Dynamics

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods allow simulations of chemical reactions in atomistic solvent and heterogeneous environments such as proteins. Herein, the nuclear–electronic orbital (NEO) QM/MM approach is introduced to enable the quantization of specified nuclei, typically protons, in the QM region using a method such as NEO-density functional theory (NEO-DFT). This approach includes proton delocalization, polarization, anharmonicity, and zero-point energy in geometry optimizations and dynamics. Expressions for the energies and analytical gradients associated with the NEO-QM/MM method, as well as the previously developed polarizable continuum model (NEO-PCM), are provided. Geometry optimizations of small organic molecules hydrogen bonded to water in either dielectric continuum solvent or explicit atomistic solvent illustrate that aqueous solvation can strengthen hydrogen-bonding interactions for the systems studied, as indicated by shorter intermolecular distances at the hydrogen-bond interface. We then performed a real-time direct dynamics simulation of a phenol molecule in explicit water using the NEO-QM/MM method. Furthermore, these developments and initial examples provide the foundation for future studies of nuclear–electronic quantum dynamics in complex chemical and biological environments.

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

Second Report of the Nuclear Data Subcommittee of the Nuclear Science Advisory Committee

The central importance of the nuclear data curated by the US Nuclear Data Program (USNDP) for clean energy generation, national security, nonproliferation, medical applications, and space exploration as well as basic science was described in a prior report issued by the DOE/NSF Nuclear Science Advisory Committee subcommittee on Nuclear Data (NSAC-ND) in September 2022. In this report, we present a set of fourteen (14) recommendations that will enhance and advance DOE-NP's stewardship of nuclear data. The first three recommendations focus on the existing core USNDP capabilities, namely: 1) Support the nuclear structure evaluation workforce to improve the currency, consistency, and accessibility of the Evaluated Nuclear Structure Data File (ENSDF); 2) Enhance nuclear reaction evaluation within the USNDP in support of the Evaluated Nuclear Data File (ENDF) through expansion of the workforce and integration of high-performance computing, automation, and machine learning and; 3) Continue atomic mass evaluation in support AME and NUBASE databases. This is followed by eight (8) recommendations representing new cross-cutting initiatives involving both measurement and evaluation to address outstanding nuclear data needs. These new initiatives require a highly trained, diverse workforce that includes personnel with expertise from both inside and outside the nuclear physics community from which evaluators have traditionally been recruited. As such, many of these initiatives are accomplished via a Topical Nuclear Data Collaborations (TNDC). A TNDC is made up of domestic and international stakeholders, subject matter and nuclear data experts, and nuclear data evaluators and features a workforce development plan to ensure that nuclear data evaluators maintain currency in the relevant applications and are seen as equity partners in the endeavor. These include: 1) Establish a coordinated effort to improve evaluation and modeling in nuclear astrophysics for stellar dynamics, multi-messenger astronomy and nucleosynthesis; 2) Initiate a TNDC to develop and maintain nuclear structure evaluation beyond discrete states, including nuclear level densities, photon strength functions and photonuclear data for improved reaction modeling, and exploring nuclear structure at finite temperature; 3) Create a TNDC to perform correlated fission data evaluation, including cross sections, fragment yields, v(A), v(E n ) for nuclear energy, national security, nonproliferation and basic science; 4) From a panel of subject matter experts to establish and annually update a roster of key decay data to nurture its accelerated dissemination including both measurement and evaluation for targeted high-value nuclides for national security, nonproliferation and medical applications; 5) Comprehensive, consistent neutron-induced structure and reaction data for nuclear energy, national security, nonproliferation and planetary nuclear spectroscopy; 6) Charged-particle stopping powers for detector design, space effects and ion beam therapy; 7) High-energy reactions for space exploration and medical nuclide production, and; 8) The creation of an infrastructure for open data and data preservation for use by the entire nuclear physics community. All told, these initiatives require approximately $6.5M increase in NP support of the USNDP in fiscal year 2023 dollars and would require at least 3-5 years to carry out due to the length of time needed to recruit and train new nuclear data researchers. This relatively modest investment would help ensure that the fruits of the nuclear data research carried out by DOE-NP and its collaborators would be brought to bear to address some of the most important needs of our nation and the world. To ensure effective execution of this plan, we present an overview of recruitment, training, and retention goals for the USNDP, the centerpiece of which is a mutually agreed upon code of conduct. Finally, we identify the facility and instrumentation needed to perform the recommended experimental activities. This includes a short review of target fabrication capabilities, reactors, neutron beam, light- and heavy-stable ion, gamma-ray, high-energy and radioactive ion beam facilities. Lastly, a more complete appendix of experimental facilities previously compiled is included with new input provided for 6 facilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Nuclei with Up to $$\varvec{A=6}$$ Nucleons with Artificial Neural Network Wave Functions

The ground-breaking works of Weinberg have opened the way to calculations of atomic nuclei that are based on systematically improvable Hamiltonians. Solving the associated many-body Schrödinger equation involves non-trivial difficulties, due to the non-perturbative nature and strong spin-isospin dependence of nuclear interactions. Artificial neural networks have proven to be able to compactly represent the wave functions of nuclei with up to $A=4$ nucleons. In this work, we extend this approach to $^6$Li and $^6$He nuclei, using as input a leading-order pionless effective field theory Hamiltonian. We successfully benchmark their binding energies, point-nucleon densities, and radii with the highly-accurate hyperspherical harmonics method.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Comparative Critical Mass Calculations for NNL and ENDF/B-VIII.0 Zirconium Hydride Thermal Neutron Scattering Laws

Zirconium hydride (ZrH x ) is a moderator material for TRIGA reactors and historical space reactor systems, such as SNAP-10A. Thermal neutron scattering laws (TSL) for two phases of this material, δ and ε, have been previously evaluated by Naval Nuclear Laboratory (NNL) and submitted to the National Nuclear Data Center (NNDC) for inclusion in the US national ENDF/B-VIII.1 nuclear data library. In contrast to the current ENDF/B-VIII.0 TSL evaluations, which consider only a single phase, the new evaluations are derived from separate ab initio calculations for both phases and include coherent elastic effects of the zirconium sublattice. To estimate the impact of these changes to the TSL evaluation of this material, comparative critical mass calculations were performed with MC21 for homogenous mixtures of high- enriched uranium (HEU) and ZrH x in bare and water reflected sphere configurations. These calculations yield an impact on the estimated critical mass as a function of 235 U loading density with maximum differences as large as 1% - 5% for over-moderated thermal spectrum systems. Consequently, the NNL TSL evaluations are anticipated to have a small impact on criticality calculations of thermal reactor systems regardless of the material phase. Nevertheless, characteristic differences exist in the predicted thermal spectra as function of energy for the two sets of TSL evaluations, which are attributed to difference in the underlying phonon density of states of hydrogen bound in ZrH x .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Insights into the possible existence of a soft dipole mode in He 8

With an extreme neutron-to-proton ratio of N/Z = 3, 8 He provides an ideal laboratory for the study of a variety of exotic phenomena, such as the emergence of a soft dipole mode that is dominated by transitions into the continuum. In this contribution, a covariant density-functional theory (DFT) framework is used to compute ground-state properties and the dipole response of 8 He. Although 8 He is admittedly too light for DFT to be applicable, the great merit of the approach is that the spurious contamination associated with the center-of-mass motion is guaranteed to decouple from the physical response. Given that a strong mixing between the isoscalar and isovector dipole modes is expected for a system with such a large neutron-proton asymmetry as 8 He, the narrow structures that emerged at low energies in the isovector dipole response are attributed to the shift of the spurious strength to zero (or near zero) excitation energy. Furthermore, the theoretical framework implemented here disfavors the emergence of a soft dipole mode in 8 He.

6 ≤ A ≤ 19↗
DOE OSTI2020

Symmetric Nuclear Matter from the Strong Interaction

We study the equation of state of symmetric nuclear matter at zero temperature over a wide range of densities using two complementary theoretical approaches. At low densities, up to twice nuclear saturation density, we compute the energy per particle based on modern nucleon-nucleon and three-nucleon interactions derived within chiral effective field theory. For higher densities, we derive for the first time constraints in a Fierz-complete setting directly based on quantum chromodynamics using functional renormalization group techniques. Here, we find remarkable consistency of the results obtained from both approaches as they come together in density and the natural emergence of a maximum in the speed of sound 𝑐𝑆 at supranuclear densities. The presence of this maximum appears tightly connected to the formation of a diquark gap. Notably, this maximum is observed to exceed the asymptotic value 𝑐$^{2}_{𝑆}$ =1/3 while its exact position in terms of the density cannot yet be determined conclusively.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

A physics-based model of cladding wastage layer formation rate

In this report, lower length scale simulations to inform an engineering-scale model of fuel-cladding chemical interaction (FCCI) conducted under the auspices of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY22 are described. An overall strategy for the implementation of the BISON model is described, and the past and current lower length scale work on FCCI formation is put into the context of this overall strategy. Density functional theory and kinetic Monte Carlo simulations are used to determine the diffusion coefficient of neodymium in iron. Density functional theory calculations are also used to parameterize a parabolic approximation for the dependence of free energy on composition for the intermetallic compound Fe 17 Nd 2 . A phase-field model of the Fe-Nd system was developed that includes the random solid solution body-centered cubic phase and the intermetallic Fe 17 Nd 2 . The model was used to simulate growth of the intermetallic layer in a diffusion couple and the results were compared to experiment. The model’s prediction of the parabolic growth constant is within 40 % of the experimental value. The model is also used to simulate the formation of a grain boundary attack layer that is observed in experiment. The simulations show enhanced Nd concentration along the grain boundaries, but do not clearly show formation of a Fe 17 Nd 2 layer

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Gamow-Teller Strength in Ca 48 and Ni 78 with the Charge-Exchange Subtracted Second Random-Phase Approximation

Here, we develop a fully self-consistent subtracted second random-phase approximation for charge-exchange processes with Skyrme energy-density functionals. As a first application, we study Gamow-Teller excitations in the doubly magic nucleus 48 Ca, the lightest double-β emitter that could be used in an experiment, and in 78 Ni, the single-beta-decay rate of which is known. The amount of Gamow-Teller strength below 20 or 30 MeV is considerably smaller than in other energy-density-functional calculations and agrees better with experiment in 48 Ca, as does the beta-decay rate in 78 Ni. These important results, obtained without ad hoc quenching factors, are due to the presence of two-particle–two-hole configurations. Their density progressively increases with excitation energy, leading to a long high-energy tail in the spectrum, a fact that may have implications for the computation of nuclear matrix elements for neutrinoless double-β decay in the same framework.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1972

Deformed Brueckner-Hartree-Fock calculations

The renormalized Brueckner-Hartree-Fock (RBHF) theory for many-body nuclear systems is generalized to permit calculations for intrinsic states having permanent deformation. Both Hartree-Fock and Brueckner self-consistencies are satisfied, and details of the numerical techniques are discussed. The Hamada-Johnston interaction is used in a study of deformations, binding, size, and separation energies for several nuclei. Electromagnetic transition rates, moments, and electron scattering form factors are calculated using nuclear wave functions obtained by angular momentum projection. Comparison is made to experiment as well as to predictions of ordinary and density-dependent Hartree-Fock Theory.

Ford, W. F.↗
NASA NTRS1972

Deformed Brueckner-Hartree-Fock calculations.

The renormalized Brueckner-Hartree-Fock (RBHF) theory for many-body nuclear systems has been generalized to permit calculations for intrinsic states having permanent deformation. Both Hartree-Fock and Brueckner self-consistencies are satisfied, and details of the numerical techniques are discussed. The Hamada-Johnston interaction is used in a study of deformations, binding, size, and separation energies for several nuclei. Electromagnetic transition rates, moments, and electron scattering form factors are calculated using nuclear wave functions obtained by angular momentum projection. Comparison is made to experiment as well as to predictions of ordinary and density-dependent Hartree-Fock theory.

Ford, W. F.↗
DOE OSTI2020

Efficient computational search for lanthanide-binding additive dopants for advanced U-Zr based fuels

A rational design strategy combining theory and experiment is highly desirable for the development of U-Zr based fuels that are resistant to lanthanide-induced fuel-cladding chemical interaction (FCCI). Here, we performed rapid computational screening of FCCI-mitigating (via strong lanthanide-binding) elements across the periodic table by utilizing three simple criteria and density functional theory (DFT) calculated total energies of relevant elements and compounds retrieved from high-throughput DFT databases. In addition to successfully identifying previously demonstrated additives, our search leads to the discovery of Bi as a highly effective element among all 68 elements screened, and its effectiveness was experimentally confirmed in this study.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Diffusion behavior of lanthanide fission products in bcc Fe cladding: A first-principles study

Fuel-cladding chemical interaction poses significant challenges in nuclear reactors, where fission products generated from nuclear fuel interact with Fe-based cladding materials, potentially compromising their structural integrity. This study investigates the diffusion behavior of lanthanide fission products, Lanthanum (La), Cerium (Ce), Praseodymium (Pr), and Neodymium (Nd), within body-centered cubic (bcc) Fe cladding using the density functional theory, nudged elastic band method, and self-consistent mean field theory. Our results reveal significant vacancy binding energies, particularly with the 1st and 2nd nearest neighbors, which diminish beyond the 5th nearest neighbor, with La exhibiting the strongest binding affinity, followed by Nd, Ce, and Pr. The nudged elastic band calculations indicate significant high barriers for the dissociation of 1st nearest neighbor vacancy-solute pairs for all fission products. The tracer diffusion coefficients of these fission products were derived in an Arrhenius form, with a magnetic correction that accounts for the high-temperature paramagnetic state. The significant trapping effect of vacancies caused by a very dilute concentration of fission products reduces vacancy mobility, potentially leading to modifications in point defect supersaturation, void nucleation, and swelling under irradiation. These represent critical challenges for irradiated cladding materials. The tracer diffusion coefficients indicate that Nd diffuses the fastest, followed by La, Ce, and Pr. Furthermore, this study provides essential insights for understanding fission product transport in cladding materials and informs future design strategies to mitigate fuel-cladding chemical interaction, ultimately enhancing nuclear reactor safety and performance.

Diffusion↗
DOE OSTI2024

Symmetry Breaking in the Lowest-Lying Excited-State of CCl4: Valence Shell Spectroscopy in the 5.0–10.8 eV Photon Energy Range

We report absolute high-resolution vacuum ultraviolet (VUV) photoabsorption cross-sections of carbon tetrachloride (CCl4) in the photon energy range 5.0–10.8 eV (248–115 nm). The molecular spectrum and electronic structure have been comprehensively investigated together with quantum chemical calculations, providing geometries, bond lengths, vertical excitation energies and oscillator strengths. The major electronic excitations have been assigned to valence and Rydberg transitions which are also accompanied by vibrational excitation assigned to degenerate stretching, v3′t2 and degenerate deformation v4′t2 modes. The rather complex nuclear dynamics along the degenerate deformation mode, v4′t2, have been thoroughly investigated by Time-Dependent Density Functional Theory (TD-DFT) method. The relevant Jahn–Teller distortion operative within the lowest-lying electronic excited-state is shown here for the first time in order to yield a weak absorption feature at 6.156 eV. Further calculations on the potential energy curves for the singlet excited-states along the C–Cl stretching coordinate show the relevance of efficient C–Cl bond excision.

Biochemistry & Molecular Biology↗