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.

SourceResearch connectionFrom yearThrough year

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

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Hyperfine interactions for small systems including transition-metal elements using self-interaction corrected density-functional theory

The interactions between the electronic magnetic moment and the nuclear spin moment, i.e., magnetic hyperfine (HF) interactions, play an important role in understanding electronic properties of magnetic systems and in realizing platforms for quantum information science applications. We investigate the HF interactions for atomic systems and small molecules, including Ti or Mn, by using Fermi–Löwdin orbital (FLO) based self-interaction corrected (SIC) density-functional theory. We calculate the Fermi contact (FC) and spin-dipole terms for the systems within the local density approximation (LDA) in the FLO-SIC method and compare them with the corresponding values without SIC within the LDA and generalized-gradient approximation (GGA), as well as experimental data. For the moderately heavy atomic systems (atomic number Z ≤ 25), we find that the mean absolute error of the FLO-SIC FC term is about 27 MHz (percentage error is 6.4%), while that of the LDA and GGA results is almost double that. Therefore, in this case, the FLO-SIC results are in better agreement with the experimental data. For the non-transition-metal molecules, the FLO-SIC FC term has the mean absolute error of 68 MHz, which is comparable to both the LDA and GGA results without SIC. For the seven transition-metal-based molecules, the FLO-SIC mean absolute error is 59 MHz, whereas the corresponding LDA and GGA errors are 101 and 82 MHz, respectively. Therefore, for the transition-metal-based molecules, the FLO-SIC FC term agrees better with experiment than the LDA and GGA results. We observe that the FC term from the FLO-SIC calculation is not necessarily larger than that from the LDA or GGA for all the considered systems due to the core spin polarization, in contrast to the expectation that SIC would increase the spin density near atomic nuclei, leading to larger FC terms.

Chemistry↗

DOE OSTI2025

The Electron Thermal Conductivity of Pu and Zr Substituted $\mathcal{γ}$-U

Uranium alloys are attractive recycled nuclear fuels because of their high thermal conductivity (𝑘) and fissile density. Limited experimental studies of the 𝑘 of U-Pu-Zr alloys in the range of 15 to 20 wt% Pu and 6 to 15 wt% Zr indicate that increasing the content of either Zr or Pu tends to lower 𝑘. However, which element has the greater effect on 𝑘, and the associated mechanisms, remains unclear. Here, in this study, the electron thermal conductivity (𝑘 𝑒 ) of U-Pu-Zr compositions are calculated using density functional theory. The electronic structure is evaluated to understand the effects of plutonium (Pu) and zirconium (Zr) substitution on the 𝑘 𝑒 of 𝛾-U. Alloys of up to 37.5 at. % Pu and 37.5 at. % Zr are examined. Two methods are applied to calculate 𝑘 𝑒 ; we find that the accuracy of each method depends on the electronic and mass similarities between the solute and solvent atoms. Specifically, when the solute atom is similar in electronic structure and mass, the more accurate method is that which employs the electron relaxation time of 𝛾-U, while if the elements are dissimilar, a mixed method that mixes several parameters associated with JNW_S⁢3033426825100132 from each element in the alloy is best. The introduction of all alloying elements decreases 𝑘 𝑒 ; however, in binary compounds, Pu and Zr have different effects. Pu flattens the electronic bands but compensates for this deleterious effect by increasing electron density near the Fermi level. Zr flattens the electronic bands more severely without adding electron density near the Fermi level. Therefore, Zr decreases 𝑘 𝑒 more than Pu in binary compounds. In ternary compounds, the difference between Pu and Zr is minimal due to the phononic change from the large mass change of Zr substitution, even at 12.5 at. %. Thus, we predict that higher loadings of Pu, and potentially other actinides, can be added to U-Pu-Zr compositions for faster recycling of spent fuel without sacrificing 𝑘. We also note that these 𝑘 𝑒 calculation methods can be applied to non-fuel alloys that require 𝑘 𝑒 predictions, such as cladding, heat exchanger, and structural materials.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Zirconium‐Based Metal–Organic Frameworks with Free Hydroxy Groups for Enhanced Perfluorooctanoic Acid Uptake in Water

Perfluorooctanoic acid (PFOA) is a highly recalcitrant organic pollutant, and its bioaccumulation severely endangers human health. While various methods are developed for PFOA removal, the targeted design of adsorbents with high efficiency and reusability remains largely unexplored. Here, in this work, the rational design and synthesis of two novel zirconium‐based metal‒organic frameworks (MOFs) bearing free ortho ‐hydroxy sites, namely noninterpenetrated PCN‐1001 and twofold interpenetrated PCN‐1002, are presented. Single crystal analysis of the pure ligand reveals that intramolecular hydrogen bonding plays a pivotal role in directing the formation of MOFs with free hydroxy groups. Furthermore, the transformation from PCN‐1001 to PCN‐1002 is realized. Compared to PCN‐1001, PCN‐1002 displays higher chemical stability due to interpenetration, thereby demonstrating an exceptional PFOA adsorption capacity of up to 632 mg g −1 (1.53 mmol g −1 ), which is comparable to the reported record values. Moreover, PCN‐1002 shows rapid kinetics, high selectivity, and long‐life cycles in PFOA removal tests. Solid‐state nuclear magnetic resonance results and density functional theory calculations reveal that multiple hydrogen bonds between the free ortho ‐hydroxy sites and PFOA, along with Lewis acid‐base interaction, work collaboratively to enhance PFOA adsorption.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Density functional theory calculations of the mixing enthalpy of ternary uranium carbide compounds

The high melting point of uranium-zirconium carbides (U,Zr)C makes them an ideal fuel for nuclear thermal propulsion (NTP) reactors. Gaps remain in the current understanding of the U-Zr-C system due to the difficulty of conducting thermodynamic experiments at NTP operation conditions. Density functional theory calculations using the Hubbard U model (DFT+U) were performed using orbital matrix occupation (OMC) to obtain the mixing enthalpy for UC and ZrC for (U,Zr)C ternary compounds. Similarly, DFT+U calculations were also carried out for the (U,Nb)C and (U,Ta)C systems. In conclusion, the DFT results are envisioned to be used in thermodynamic assessments of the uranium carbide systems based on the CALPHAD approach to supplement the lack of experimental data for the mixing thermodynamics.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

Structure and Dynamics of Imidazolium in an Ionic Liquid-PEGDA Iongel via IR, 2D-IR, and NMR Spectroscopy

Fourier transform infrared (FTIR) absorption spectroscopy, ultrafast two-dimensional infrared (2D-IR) spectroscopy, nuclear magnetic resonance (NMR) spectroscopy, and density functional theory (DFT) calculations show that, in iongels composed of 1-ethyl-3-methylimidazolium bis(trifluoromethyl\-sulfonyl)imide ([C 2 C 1 Im][Tf 2 N]}) and cross-linked poly(ethylene glycol)diacrylate (cl-PEGDA), imidazolium cations solvate PEGDA chains in tight and loose complexes. H/D isotope substitution at the 2-position of the imidazolium ring simplifies the CH-stretching vibrational bands. Tight complexes are characterized by stronger, directional hydrogen bonding interactions between the 2-position of the imidazolium ring and the ethereal oxygens. The loose complexes resemble bulk ionic liquid (IL) non-specifically perturbed by the PEGDA chains.

ionic liquids, carbon dioxide, cross-linked PEGDA,↗

DOE OSTI2025

Tuning the Properties of Lignin‐Derived Deep Eutectic Solvents for Biomass Processing

Lignin-derived deep eutectic solvents (DESs) have been investigated as sustainable green media for biomass processing. However, the properties and processability of DESs have not been fully understood with the chemical structures of their constituents for biomass fractionation. Here, in this article, the properties of the phenolic DESs are discussed with different numbers of functional groups, such as –OCH 3 and –CHO in their hydrogen bond donor (HBD) structures. The formation of DES is significantly related to the hydrogen bond between its constituents, identified by nuclear magnetic resonance (NMR) analysis and density functional theory calculation (DFT). Lower viscosity and net basicity of DES are achieved with fewer –OCH 3 on HBD structures, resulting in enhanced processability and fractionation efficiency. The thermal stability of the DES is also influenced by the –OCH 3 and –CHO of HBD, as indicated by its onset temperature. The recyclability of the phenolic DES is confirmed by the fractionation performance of the recycled DES. Understanding the structural impacts of DES constituents on the properties and performance is crucial for designing solvents in biorefinery applications.

Ryu, Jiae [State University of New York (SUNY), Sy↗

DOE OSTI2026

Uncovering the linear boron environment in Na 3 BP 2 through solid-state 11 B NMR spectroscopy

Boron-based compounds exhibit a wide range of structural diversity, with potential applications spanning organic and inorganic chemistry. Herein, we focus on the characterization of the linear boron-phosphorus unit P═B═P in Na 3 BP 2 using solid-state nuclear magnetic resonance (ssNMR) spectroscopy and density functional theory (DFT) calculations. High-resolution 11 B ssNMR spectra were recorded at two fields, and key parameters such as chemical shift anisotropy (CSA), quadrupolar coupling constants (C Q ), and electric field gradient (EFG) tensors were extracted. The 11 B NMR results revealed a distinct chemical environment for the two-coordinate boron atom, with a CSA span (Ω) of 280 ppm and a C Q of 3.0 MHz. These values were further validated through periodic plane-wave DFT calculations, which showed good agreement with experimental results. The obtained spectral parameters are compared to other linear boron units, such as the BO 2 motif, providing a broader context for understanding boron coordination in inorganic compounds. This work expands the body of NMR knowledge on boron-containing materials, particularly for linear boron motifs. The findings contribute to the growing field of boron chemistry and its potential applications in advanced materials.

Porter, Andrew P. [Ames Laboratory (AMES), Ames, I↗

DOE OSTI2024

Molecular NMR shieldings, J -couplings, and magnetizabilities from numeric atom-centered orbital based density-functional calculations

This paper reports and benchmarks a new implementation of nuclear magnetic resonance shieldings, magnetizabilities, and J-couplings for molecules within semilocal density functional theory, based on numeric atom-centered orbital (NAO) basis sets. NAO basis sets are attractive for the calculation of these nuclear magnetic resonance (NMR) parameters because NAOs provide accurate atomic orbital representations especially near the nucleus, enabling high-quality results at modest computational cost. Moreover, NAOs are readily adaptable for linear scaling methods, enabling efficient calculations of large systems. Here, the paper has five main parts: (1) It reviews the formalism of density functional calculations of NMR parameters in one comprehensive text to make the mathematical background available in a self-contained way. (2) The paper quantifies the attainable precision of NAO basis sets for shieldings in comparison to specialized Gaussian basis sets, showing similar performance for similar basis set size. (3) The paper quantifies the precision of calculated magnetizabilities, where the NAO basis sets appear to outperform several established Gaussian basis sets of similar size. (4) The paper quantifies the precision of computed J-couplings, for which a group of customized NAO basis sets achieves precision of ~Hz for smaller basis set sizes than some established Gaussian basis sets. (5) The paper demonstrates that the implementation is applicable to systems beyond 1000 atoms in size.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

Structural Studies of Alloyed and Nanoparticle Transition Metal Dichalcogenides by Selenium-77 Solid-State Nuclear Magnetic Resonance Spectroscopy

Layered transition metal dichalcogenides (TMDCs) such as MoS 2 , MoSe 2 , and WSe 2 are under intense investigation because they are atomically thin semiconductors with photophysical properties that can be tuned by changing their composition or morphology. Mechanochemical processing has been proposed as a method to obtain alloyed TMDCs in the series Mo 1–x W x S y Se 2–y (x = 0–1; y = 0–2). However, elucidating the chemical transformations occurring at the atomic scale following mechanochemical processing can be challenging because the products are often amorphous or microcrystalline. To address this challenge, we probe TMDC mixing and alloying by using a combination of powder X-ray diffraction, Raman spectroscopy, diffuse reflectance spectroscopy, 77 Se solid-state nuclear magnetic resonance (SSNMR) spectroscopy, and planewave density functional theory (DFT) calculations. The nature of the milling material and reaction atmosphere are shown to be essential factors in limiting the formation of undesired oxide byproducts. We demonstrate acquisition of 77 Se SSNMR spectra using different combinations of Carr-Purcell Meiboom-Gill acquisition (CPMG) pulse sequences, magic angle spinning (MAS), and MAS dynamic nuclear polarization. Further, the combination of SSNMR with the other characterization methods clearly demonstrates that high energy impact ball milling induces molecular level alloying of Mo, W and chalcogen atoms in the family Mo 1–x W x S y Se 2–y . Gauge including projector augmented wave DFT calculations yield accurate 77 Se chemical shift (CS) tensor components. 77 Se SSNMR spectroscopy was also applied to study the structure of WSe 2 nanocrystals intercalated with ethylenediamine. The intercalated WSe 2 nanocrystals exhibit a more positive isotropic 77 Se CS as compared to bulk WSe 2 , however, the 77 Se CS anisotropy is the same, confirming the WSe 2 layers have a similar structure as in their bulk counterparts.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Localizing tetrahedral aluminum in nitrate-bearing gibbsite to constrain defect-impurity coupling

The enhanced radiolytic stability of gibbsite (α-Al(OH) 3 ) containing trace nitrate (NO 3 − ) is a phenomenon in nuclear waste management, but its structural origins remain unresolved. Motivated by the detection of minority tetrahedral aluminum (T d ) defects in synthetic gibbsite, we hypothesized that these sites may participate in NO 3 − retention or mediate H 2 suppression. To evaluate this, we combined orthogonal techniques comprised of spatially selective solid-state 27 Al MAS NMR, comparative spectroscopy, and density functional theory (DFT) modeling. Paramagnetic editing and dynamic nuclear polarization (DNP) MAS NMR confirm that T d defects are confined to the particle interior. DFT calculations reveal no energetic stabilization of NO 3 − near T d sites. Comparative NMR analysis shows that T d is also present in chloride-bearing gibbsite, which exhibits high radiolytic hydrogen yields. These three independent disqualifications rule out T d as a structural contributor to nitrate-mediated suppression and narrow the scope of defect-driven explanations. The findings redirect mechanistic attention away from coordination defects and toward redox-active impurity pathways, providing a refined foundation for understanding radiation tolerance in Al(OH) 3 .

Graham, Trent R. [Pacific Northwest National Labor↗

DOE OSTI2025

Strong Proton‐Phonon Coupling Drives Fast Ion Transport in Perovskites

Conduction of protons in solids is a cooperative process propelled by phonons, with molecular details obscured by the irregular movements in the thermal bath. It is shown that substitution with Y forms an imaginary phonon mode, instrumental for the function as proton conductor and effectively lowering the activation barrier for proton transport. To untangle the interplay in the exemplary proton conductor BaSn 0.9 Y 0.1 O 3 , its crystallographic structure is determined with high resolution neutron diffractometry and its phonon density of states with density functional theory calculations, experimentally validated by element specific nuclear resonant vibration spectroscopy. Based on phonon analysis, a quantitative transport model is present, which predicts the activation energy and performance by the ratio of ionic radii. Rather than individual vibrational modes, it is the oxygen sub-lattice which exerts its momentum on the protons. The extent of this momentum transfer is governed by the ratio of ionic radii. This model extends the transition state theory by the phonon-phonon interaction and complements the previously proposed idea that lattice dynamics is decisive for proton transport and specifies which properties of the material exactly define the vibration properties.

activation barriers↗

DOE OSTI2026

Interstitialcy-mediated transport in Hematite

Defect transport mechanisms drive the performance of materials in both energy storage and energy generating systems. Iron oxide systems are of particular interest with use cases in both realms. Previous studies of a-Fe2O3 (hematite) determined that the primary point defects at elevated temperature are O vacancies and Fe interstitials, and these defects determine transport rates at temperatures >850°C. However, in advanced nuclear reactor designs, operating temperatures will range between 400-800°C. In this study, we utilize isotopic tracers and the nanoscale resolution of atom probe tomography to determine diffusivities at lower temperatures relevant to nuclear energy generating systems. Supported by a density functional theory informed rate theory model, we find a change from vacancy to interstitialcy-mediated transport < 800°C in the anion lattice.

Yano, Kayla H.↗

DOE OSTI2025

Nuclear quantum effects of metal surface-mediated C–H activation

The nuclear quantum effects of surface-mediated C–H activation of surface CH 3 are considered for the pristine Pt(111) and Au(111) surfaces at 300 K. The kinetic barriers without nuclear quantum effects are calculated using both static density functional theory calculations and ab initio molecular dynamics. Static calculations are performed using the harmonic approximation while the free energy pathway is calculated using enhanced sampling molecular dynamics. Machine learning potentials are trained using generated datasets and validated against the ab initio molecular dynamics generated free energy pathways. The machine learning potentials are used to perform centroid molecular dynamics to consider the nuclear quantum effects of C–H activation. Nuclear quantum effects are found to have a very significant effect on the free energy pathway, with reduced importance at higher temperatures and in the CD 3 case.

Bunting, Rhys J. [Lawrence Livermore National Labo↗

DOE OSTI2024

Crystal structure of perfluorononanoic acid, C 9 HF 17 O 2

The crystal structure of perfluorononanoic acid (PFNA) was solved via parallel tempering using synchrotron powder diffraction data obtained from the Brockhouse X-ray Diffraction and Scattering (BXDS) Wiggler Lower Energy (WLE) beamline at the Canadian Light Source. PFNA crystallizes in monoclinic space group P2 1 /c (#14) with lattice parameters a = 26.172(1) Å, b = 5.6345(2) Å, c = 10.9501(4) Å, and β = 98.752(2)°. The crystal structure is composed of dimers, with pairs of PFNA molecules connected by hydrogen bonds via the carboxylic acid functional groups. The Rietveld-refined structure was compared to a density functional theory-optimized structure, and the root-mean-square Cartesian difference was larger than normally observed for correct powder structures. The powder data likely exhibited evidence of disorder which was not successfully modeled.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2024

Accelerating multiscale electronic stopping power predictions with time-dependent density functional theory and machine learning

Knowing the rate at which particle radiation releases energy in a material, the “stopping power,” is key to designing nuclear reactors, medical treatments, semiconductor and quantum materials, and many other technologies. While the nuclear contribution to stopping power, i.e., elastic scattering between atoms, is well understood in the literature, the route for gathering data on the electronic contribution has for decades remained costly and reliant on many simplifying assumptions, including that materials are isotropic. We establish a method that combines time-dependent density functional theory (TDDFT) and machine learning to reduce the time to assess new materials to hours on a supercomputer and provide valuable data on how atomic details influence electronic stopping. Our approach uses TDDFT to compute the electronic stopping from first principles in several directions and then machine learning to interpolate to other directions at a cost of 10 million times fewer core-hours. We demonstrate the combined approach in a study of proton irradiation in aluminum and employ it to predict how the depth of maximum energy deposition, the “Bragg Peak,” varies depending on the incident angle—a quantity otherwise inaccessible to modelers and far outside the scales of quantum mechanical simulations. The lack of any experimental information requirement makes our method applicable to most materials, and its speed makes it a prime candidate for enabling quantum-to-continuum models of radiation damage. The prospect of reusing valuable TDDFT data for training the model makes our approach appealing for applications in the age of materials data science.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Nature of the Reactive Biferric Peroxy Intermediate P′ in the Arylamine Oxygenases and Related Binuclear Fe Enzymes

Binuclear nonheme iron enzymes activate O 2 to perform a wide range of chemical transformations. The process of O 2 activation typically involves a biferric peroxy-level intermediate P. It has been previously found that this intermediate undergoes further activation, either protonation or rearrangement to form P′ or further oxidation to form high-valent intermediates Q or X. Here, this study defines the structure of the P′ intermediate in the N-oxygenases CmlI (and AurF based on previous data) using nuclear resonance vibrational spectroscopy (NRVS) in conjugation with density functional theory (DFT) calculations. These results, combined with variable temperature variable field (VTVH) magnetic circular dichroism (MCD) spectroscopy on the 1-electron cryoreduced P′, define the structure of the P′ intermediate as a μ-1,2-hydroxoperoxo biferric site with a second hydroxide bridge. Reaction coordinate calculations demonstrate that single electron transfer (SET) is facilitated by protonation of the peroxo, activating its reductive cleavage, and that the additional hydroxide bridge does not impact this reaction. VTVH MCD studies further reveal that the hydroxide bridge is absent in the biferrous site, suggesting that during the O 2 reaction with the biferrous site, a water molecule forms the hydroxide bridge in providing the proton that activates the peroxide in P′ for reactivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Computational insights into the structure of anhydrous Pu(III) oxalate

Despite the plutonium oxalate method's wide use in plutonium reprocessing, the method's mechanistic details remain unclear, particularly the identity of the plutonium oxidation state during conversion from an oxalate hydrate to the oxide. Recently, the optical vibrational spectra of Pu(III) oxalate during calcination were measured, providing an experimental reference for computational studies aimed at elucidating the oxalate structure. Here, in this work, we compare the vibrational and electronic properties of two candidate anhydrous Pu(III) oxalate structures calculated using density functional theory with recent experiments. We find that both structures are plausible and may coexist at experimental temperatures, providing insights into the broad features measured in the Raman and infrared spectra.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Machine learning-accelerated path integral molecular dynamics simulations of reactive organic electrolytes

Hydrogen bonded electrolytes that exhibit accelerated proton transport via sequential reactive hops have drawn interest for their promise in clean energy applications. Molecular dynamics simulations of these electrolytes offer the opportunity to uncover microscopic mechanistic details that could be used to design and tune the properties of candidate electrolyte technologies. However, accurately modeling the proton transfer reactions and transport properties that give rise to high charge conductivites in these electrolytes proves computationally challenging because of the need to perform lengthy condensed phase simulations, treating both the electronic and nuclear degrees of freedom quantum mechanically. In this paper, we demonstrate that such a modeling task can be efficiently achieved with the use of density functional theory (DFT)-trained machine learning potentials (MLP) to accelerate path integral molecular dynamics (PIMD) simulations. We highlight the practical utility of this approach by using it to benchmark how closely PIMD simulations employing different DFT exchange–correlation functionals reproduce the composition-dependent densities, diffusion coefficients, and electrical conductivities of mixtures consisting of imidazole and levulinic acid. Even with the speedup afforded by our MLPs, PIMD simulations remain quite expensive. Furthermore, in order to render PIMD more computationally tractable, we introduce and benchmark the accuracy of a ring polymer contraction approach that leverages a computationally efficient short-range MLP to accelerate our PIMD simulations by an additional factor of four.

Chemical bonding↗