Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear density functional theory”

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

Prediction of Redox Potentials for U, Np, Pu, and Am in Aqueous Solution

The redox properties of the actinides in aqueous solution are important for fuel production/reprocessing and understanding the environmental impact of nuclear waste. The redox potentials for U, Np, Pu, and Am in oxidation states from 0 up to VII (as appropriate) in aqueous solutions have been predicted at the density functional theory level with the B3LYP functional, Stuttgart small core pseudopotential basis sets for the actinides, and explicit (30H 2 O molecules)/implicit treatment of the aqueous solvent using the self-consistent reaction field COSMO and SMD approaches for the implicit solvation. The predictions of the structural parameters of clusters incorporating first and second solvation shells are consistent with the available experimental data., Our results are typically within 0.2 V of the available experimental data using two explicit solvation shells with an implicit solvent model. The use of the PW91 functional substantially improved the prediction of the Pu(VI/V) redox couple. The redox couples for An(VI/IV) and An(V/IV) which involve the addition of protons and removal of the actinyl oxygens led to slightly larger differences from experiment. Here, the An(IV/0) and An(III/0) couples were reliably predicted with our approach. Predictions of the unknown An(II/I) redox potentials were negative, consistent with expectations, and predictions for unknown An(VII/VI), An(III/II), and An(II/0) redox couples improve prior estimates.

Actinides↗
DOE OSTI2025

Prediction of Solute Segregation at Metal/Oxide Interfaces Using Machine Learning Approaches

The atomic structure and chemistry at metal/oxide interfaces play a crucial role in determining their properties. However, studying semi-coherent metal/oxide interfaces that include misfit dislocations through density functional theory (DFT) is often computationally expensive due to the large number of atoms involved, ranging from hundreds to thousands. In this study, we explore solute segregation behavior at the Fe/Y 2 O 3 interface—an important model interface for cladding applications in nuclear fission reactors—by combining DFT calculations with a machine learning (ML) approach. ML models are trained using DFT-calculated segregation energies (𝐸 𝑆𝑒𝑔 ) to identify the key chemical and geometric factors influencing solute segregation at metal/oxide interfaces, revealing the competition between these features in determining 𝐸 𝑆𝑒𝑔 . Moreover, the segregation behavior at a specific Fe/Y 2 O 3 interface is predicted with high accuracy using ML models trained on data from this interface. Furthermore, it is found that the ML models could also predict solute segregation at a different Fe/Y 2 O 3 interface with a new orientation relationship (OR), at a computational cost of less than 1/45 of that required for similar DFT calculations.

36 - MATERIALS SCIENCE↗
DOE OSTI2024

Effect of impurities on hydrogen defect stability and migration barrier in yttrium dihydride crystal

The impurity or alloying atoms in YH 2 can alter the local electronic structure and so the hydrogen defect stability, as well as the H migration barrier energy. Thus, DFT calculations were employed to determine the effect of foreign elements from alkali and alkaline earth metals to transition metals and one critical impurity element, O, on H vacancy stability and retention characteristics in YH2. Results revealed that alloying elements act as hydrogen vacancy sinks by reducing the vacancy formation energy at neighboring sites. The implantation of non-magnetic foreign elements (s1, s2, and d10 valence electrons) in hydrogen energy landscape was calculated to be minor; while the hydrogen vacancy formation energy was reduced from 1.37 eV to 1.00 eV, the migration energy barrier of hydrogen was increased from 0.87 eV to 1.15 eV for non-magnetic foreign elements. The migration energy barrier monotonically decreased with increasing d-shell occupancy, reaching as low as 0.4 eV for Cr, Mo(d4), and Fe (d4). Alloying with late transition metals (d8 and d9) moderately impacted the hydrogen vacancy formation. Finally, it was found to be O addition into the YH 2- lattice did not alter the energy landscape of hydrogen vacancies. Since alloyed YH 2 has not been studied extensively, this study provides an atomistic understanding how alloying elements and impurities trap vacancies and affects hydrogen mobility YH 2 . Meanwhile, the main findings of this study may serve as guidelines for introducing alloying elements in ZrH 2 as well.

08 HYDROGEN↗
DOE OSTI2025

Hydrolysis Reaction Pathways of Thorium Oxide Nanoclusters

Density functional theory benchmarked by correlated molecular orbital theory is used to develop a fundamental and predictive understanding of the interaction of thorium oxide nanoclusters with gas phase water to provide insight into nuclear-waste storage, production of thorium nuclear reactor fuels, and reprocessing of spent fuel. The structures of Th n O 2n (n = 3 – 6) clusters and their interactions with water have been studied at the B3LYP, MP2, and CCSD(T) levels. Hydrolysis is initiated by the formation of Lewis acid-base adducts, with relative H 2 O binding energies (physisorption) ranging from −15 kcal/mol to −22 kcal/mol. The initial H 2 O physisorption energy is ca. −21 kcal/mol regardless of the cluster size and is consistent with the experimentally obtained initial adsorption energy on a thorium dioxide surface. The physisorption enthalpies for additional water molecules can be affected by the presence of terminal groups OH groups generated by proton transfer to a Th-O near the site of adsorption. The hydrolysis products (chemisorption) form either bridging or terminal hydroxides. More exothermic hydrolysis steps were predicted for the formation of terminal hydroxides as compared to the formation of bridging hydroxides. Here, the calculated transition state barriers for transfer of protons from bound water complexes to form the chemisorption products are very low. Overall, water readily reacts with thorium oxide clusters preferring hydroxide products over hydrated complexes. First and second order fits were predicted for the combined physisorption and chemisorption energies for the hydrolysis of thorium oxide clusters. Finally, ionization energies and electron affinities were calculated as were HOMO-LUMO gaps to provide additional insights into the properties of the thorium oxide and hydroxide clusters.

Adsorption↗
DOE OSTI2026

An Ab Initio Molecular Dynamics Study of Key Thermodynamic Input Parameters for Computer Simulation of U-6Nb Solidification

The key to metallic fuel development is the fabrication of uranium metal and alloys into fuel forms. U-Nb alloys are one of the best candidates for a metallic fuel alloy with high-temperature strength sufficient to support the core, acceptable nuclear properties, good fabricability, and compatibility with usable coolant media. Melt processing has been a key component of the metallic fuel cycle, and process models require thermophysical parameters at elevated temperatures, particularly above the melting temperatures, regarding which experimental data are scarce, for accurate simulations and process development. By means of ab initio density-functional theory (DFT) quantum molecular dynamics (QMD), we have calculated the main thermophysical parameters—the density, thermal expansion coefficient, specific heat, thermal conductivity, melting temperature, latent heat of fusion, and viscosity—used in the modeling of the U-6 wt.% Nb alloy casting. The melting temperature of the U-6 wt.% Nb alloy at ambient pressure is obtained by means of QMD simulations using the Z-method. The ambient volume change and latent heat of melting of U-6 wt.% Nb are also derived from QMD simulations in conjunction with analytical fitting for the energy and pressure. The thermal conductivity for the solid U-Nb alloy is calculated from the semi-classical Boltzmann transport equation combined with an estimate of the electron relaxation time obtained from DFT simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2019

A NICER view of PSR J0030+0451: Implications for the Dense Matter Equation of State

Both the mass and radius of the millisecond pulsar PSRJ0030+0451 have been inferred via pulse-profile modeling of X-ray data obtained by NASA’s Neutron Star Interior Composition Explorer (NICER) mission. In this Letter we study the implications of the mass–radius inference reported for this source by Riley et al. for the dense matter equation of state (EoS), in the context of prior information from nuclear physics at low densities. Using a Bayesian framework we infer central densities and EoS properties for two choices of high-density extensions: a piecewise-polytropic model and a model based on assumptions of the speed of sound in dense matter. Around nuclear saturation density these extensions are matched to an EoS uncertainty band obtained from calculations based on chiral effective field theory interactions, which provide a realistic description of atomic nuclei as well as empirical nuclear matter properties within uncertainties. We further constrain EoS expectations with input from the current highest measured pulsar mass; together, these constraints offer a narrow Bayesian prior informed by theory as well as laboratory and astrophysical measurements. The NICER mass–radius likelihood function derived by Riley et al. using pulse-profile modeling is consistent with the highest-density region of this prior. The present relatively large uncertainties on mass and radius for PSR J0030+0451 offer, however, only a weak posterior information gain over the prior. We explore the sensitivity to the inferred geometry of the heated regions that give rise to the pulsed emission, and find a small increase in posterior gain for an alternative (but less preferred) model. Lastly, we investigate the hypothetical scenario of increasing the NICER exposure time for PSRJ0030+0451.

G Raaijmakers↗
DOE OSTI2025

NNL.Fe.qSNAP-ZBL.2024.1: A Fe Spectral Neighbor Analysis Potential for Radiation Damage Simulations

The NNL.Fe.qSNAP-ZBL.2024.1 machine-learned potential (MLP) has been generated to support the development of an elemental body-centered cubic (BCC) Fe athermal recombination corrected neutron damage model and simulations of primary recoil atom (PRA) cascades in BCC Fe. This MLP is a quadratic spectral neighbor analysis potential (qSNAP) hybridized with the universal Ziegler-Beirsack-Littmark (ZBL) potential at short-range and is named according to Naval Nuclear Laboratory MLP naming conventions (NNL.material-system.MLP-type.year.version). Training set calculations for Fe are presented along with the subsequent MLP fitting procedure. A key criterion of the fitting procedure is that ZBL describes the short-range interaction with minimal impact on the MLP. The MLP is compared to density functional theory (DFT) predicted properties relevant to radiation damage simulation, including threshold displacement energies, for validation. The NNL.Fe.qSNAP-ZBL.2024.1 potential is considered suitable for molecular dynamics (MD) simulations of radiation defects up to 800 K and PRA cascades in BCC Fe up to around 10 keV. The potential can additionally be used on a limited basis for recoils of 10–20 keV, within which range the emergence of structures outside the training set in cascade simulations may cause system instabilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1983

Numerical models for cometary dust environments

The development of a numerical model for the dust environment near a comet is described. The major components of the model are a simple treatment of the acceleration of the dust by the gas outflow from an isothermal icy nucleus, a detailed treatment of the dust particle trajectories controlled by solar gravity and radiation pressure, including the formation of the dust envelope and tail, and a simple treatment of the radiation from the dust in the visible and thermal infrared. The results of recent comet observations are used to estimate parameter values for the model, including dust-to-gas ratios by light and by mass, the distribution of dust and gas emission on the nuclear surface, the distributions of dust particle size, density, mass, and radiative properties, and other quantities. Gas and dust production rates as functions of time are derived from observed light curves using the simple photometric theory. Sample results are presented for the dust concentration, flux and fluence along the trajectory of a spacecraft which flies by Comet P/Halley in 1986.

Divine, N.↗
DOE OSTI2025

Electronic and Geometric Contributors to Hydrogen Binding in Uranium Oxide Grain Boundaries

Hydrogen induced corrosion of uranium, which leads to the formation of toxic and pyrophoric UH 3 , raises significant safety concerns for long-term storage of nuclear materials. Previous work suggests hydrogen diffuses through the grain boundaries (GBs) of the passivating oxide layer to initiate hydriding reactions. However, the atomistic mechanisms underlying this phenomenon and the structural factors that control its initiation are not well understood. To address this knowledge gap, here we use a high-throughput density functional theory (DFT) workflow to investigate the adsorption of H and H 2 in the defective bulk UO 2 . Specifically, we have exhaustively investigated the adsorption of H (107 sites) and H 2 (26 sites) in three different coincident site lattice (CSL) GBs: Σ3, Σ5, and Σ9. Compared to the binding energies in pristine UO 2 , we observe significantly stronger hydrogen adsorption at these GB sites. Interestingly, we find that the trends in H and H 2 adsorption vary considerably across the three GB models. In particular, while a small number of sites in Σ5 and Σ9 show exothermic adsorption of H and H 2 , respectively, no such sites are found in Σ3. These results provide fundamental atomistic insights that could guide the development of future corrosion mitigation strategies for the storage of nuclear materials.

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

Electron-Ion Dynamics with Time-Dependent Density Functional Theory: Towards Predictive Solar Cell Modeling

This project focussed on two aspects of computational modeling with a view toward application for photovoltaic design: (i) increased reliability of exchange-correlation functionals in time-dependent density functional theory (TDDFT) (a method of choice of the calculation of electronic spectra and dynamics), especially for time-resolved non-perturbative dynamics, and (ii) the development of a practical but rigorously-based method for coupling electronic and nuclear motion via the exact-factorization (EF) approach (a relatively new framework for developing approximations).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

First-principles investigation of cerium and neodymium diffusion in BCC chromium and vanadium via vacancy-mediated transport

Lanthanide transport plays a crucial role in the performance and longevity of metallic nuclear fuels. This study examines the diffusion behavior of Ce and Nd—two major fission products—in body-centered cubic (BCC) Cr and V, which are potential liner or coating materials for mitigating fuel-cladding chemical interactions (FCCI). Using density functional theory (DFT) calculations and self-consistent mean-field (SCMF) analysis, the vacancy-mediated diffusion coefficients are evaluated. Our findings reveal that Ce and Nd act as oversized solutes and are strongly bound to vacancies in BCC Cr and V, with diffusivities in Cr significantly lower than in V and in hexagonal closed-packed (HCP) Zr, as investigated in our previous work. The activation energies for Ce and Nd diffusion are 3.39 and 3.32 eV, respectively, in BCC Cr, and 2.56 and 2.33 eV, respectively, in BCC V. Analysis of vacancy drag and partial diffusion coefficient ratios indicates a strong tendency for lanthanide enrichment at vacancy sinks in BCC Cr, and to a lesser extent in BCC V, with this effect persisting up to the melting point in Cr and remaining substantial for Nd in V at high temperatures. Under irradiation, the increase in vacancy concentration is expected to enhance lanthanide transport, potentially accelerating interactions at liner-cladding interfaces. Although BCC Cr exhibits relatively low lanthanide diffusivities under equilibrium conditions, the expected segregation tendencies under irradiation suggest that Zr liners may be a more favorable option. Further investigations using rate theory, cluster dynamics, and phase-field modeling are required to quantitatively assess the performance of these materials in reactor environments.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

First-Principles insights into group-V impurities and their impact on germanium detector performance

The outstanding properties of high-purity germanium (HPGe) detectors, such as excellent energy resolution, high energy sensitivity, and a low background-to-signal ratio, make them essential and ideal candidates for detecting particle signatures in nuclear processes such as neutrino-less double beta decay (0νββ). However, the presence of defects and impurities in HPGe crystals can lead to charge trapping, which affects carrier mobility and results in significant energy resolution degradation. In this work, we employ density functional theory with a hybrid functional to study the energetics of possible point defects in Ge. Our findings indicate that group-V impurities form more readily in Ge compared to vacancy and interstitial of Ge. Unlike N dopants, which yield deep trap states, P, As, and Sb create shallow traps close to the conduction band edge of Ge. Furthermore, we predict that group-V defects can condense into defect complexes with Ge vacancies. These vacancy-impurity complexes form deep traps in Ge, similar to Ge vacancies, suggesting that both vacancies and vacancy-impurity complexes contribute to charge trapping in these detectors, thereby diminishing their performance.

36 MATERIALS SCIENCE↗
DOE OSTI2024

In situ characterization of metastable Pb3O5 and Pb2O3 phases during thermal decomposition of PbO2 to PbO

Nonstoichiometric lead oxides play a key role in the formation and cycling of the positive electrodes in a lead acid battery. These phases have been linked to the underutilization of the positive active material but also play a key role in the battery’s cycle life, providing inter-particle adhesion and the connection to the underlying lead grid. Similar phases have previously been identified by mass loss or color change during thermal annealing of PbO2 to PbO, suggesting that at least two intermediate PbOx phases exist. Using multiple in situ analysis techniques (powder diffraction, x-ray absorption, x-ray photoelectron spectroscopy) and ex situ nuclear magnetic resonance measurements, the structural conversion and changes in the lead oxidation states were identified during this process. Isolation of the PbOx phases enabled confirmation of Pb3O5 and Pb2O3 by diffraction and the first 207Pb NMR measurement of these intermediates.

Kinnibrugh, Tiffany L.↗
DOE OSTI2025

Coherent Spins in van der Waals Semiconductor GeS 2 at Ambient Conditions

Optically active spin defects in van der Waals (vdW) materials have emerged as versatile quantum sensors, enabling applications for a wide range of quantum phenomena in condensed matter systems. Their ease of exfoliation and compatibility with device integration make them promising candidates for future quantum technologies. Here we report the observation and room-temperature coherent control of ensemble spin defects in the high-temperature crystalline phase of germanium disulfide (β-GeS 2 ), a two-dimensional (2D) semiconductor with low nuclear spin density. The defects exhibit spin-1/2 behavior, and their dynamics can be explained by a weakly coupled spin-pair model. We implement dynamical decoupling techniques to extend the coherence time (T 2 ) by a factor of 20. Finally, we use density functional theory (DFT) calculations to estimate the structures and spin densities of two possible spin defect candidates. This work will help to expand the field of quantum sensing with spin defects in 2D materials.

2D materials↗
DOE OSTI2025

Determining the Conformation of Supported Complexes Using an 17 O TEDOR-like NMR Experiment

Dynamic nuclear polarization surface-enhanced nuclear magnetic resonance (NMR) spectroscopy has enabled the determination of the three-dimensional configuration of surface sites, in particular supported metal complexes of relevance to single-site heterogeneous catalysis. These approaches have chiefly leveraged the application of NMR double-resonance experiments that either reveal the complex conformation via point-to-point intramolecular distances between spin-labeled atoms or the complex-surface orientation via distances between the spins and the surface plane. Either method typically requires expensive isotope labeling and each reports on different structural features. The application of an experiment that simultaneously reveals both types of distances with chemical resolution would be ideal. Here, in this article, we describe an 17 O{ 1 H} pseudo-3D correlation experiment that achieves this goal. Specifically, Si–O–Si and Si–O–M oxygens are well-resolved by 17 O NMR; therefore, distances can be simultaneously measured radially, between Si– 17 O–M and the 1 H’s of the ligands, and vertically to the Si– 17 O–Si linkages of the silica support. We demonstrate the experiment using supported yttrium and zirconium complexes. Good agreement is obtained when comparing the experimental results to theoretical predictions from density functional theory calculations, highlighting the reliability of this relatively simple experiment.

alkyls↗
DOE OSTI2025

Reliable p K a Prediction through Efficient Incorporation of Anharmonicity within the Nuclear–Electronic Orbital Framework

Accurate pK a prediction is critical for understanding chemical reactivity and molecular properties across a wide range of applications. Computational approaches usually invoke a harmonic treatment of the vibrational modes for zero-point energies, as well as thermal and entropic contributions. Herein, we present a general protocol for relative pK a prediction that incorporates the significant anharmonic effects using nuclear–electronic orbital (NEO) theory. This protocol is validated against experimental data for a range of molecules in acetonitrile, including protonated nitrogen bases, nitrophenols, anilines, and diamines, as well as cobalt electrocatalysts. For simple acids, the NEO approach offers only a slight improvement over conventional density functional theory with the standard harmonic vibrational treatment, whereas for hydrogen-bonded acids, the NEO approach offers more significantly improved performance at a comparable computational cost. This accessible methodology provides a practical route for accurate pKa prediction in challenging systems and is extendable to related thermodynamic properties such as hydricities and proton-coupled redox potentials.

Density functional theory↗
DOE OSTI2024

Generative deep-learning reveals collective variables of Fermionic systems

Complex processes of fermionic systems ranging from protein folding to nuclear fission often follow a low-dimensional reaction path parametrized in terms of a few collective variables. In nuclear theory, variables related to the shape of the nuclear density in a mean-field picture are key to describing the large amplitude collective motion of the neutrons and protons. Exploring the adiabatic energy landscape spanned by these degrees of freedom reveals the possible reaction channels while simulating the dynamics in this reduced space yields their respective probabilities. Unfortunately, this theoretical framework breaks down whenever the systems encounters a quantum phase transition with respect to the collective variables. Here, in this study, we introduce a novel generative deep-learning algorithm designed to build reaction paths that ensure that the many-fermion wave function stays differentiable with respect to the collective variables. This approach is applicable to any fermionic system described by a coherent state. We use the case of potential energy curves in the 16 O nucleus within the Hartree-Fock theory to illustrate its main features.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗