SEARCH · Engineering Papers

Results for “Density Functionals”

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 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Toward a better understanding of Ni coarsening in solid oxide cells: NiH on Ni (111) examined at the level of density-functional theory (and KMC)

Ni coarsening in the fuel electrode of solid oxide cells (SOCs) is known to be significantly faster under a humid atmosphere. In this talk, density-functional theory and kinetic Monte Carlo methods are used to explore the hypothesis that the surface diffusion of NiH on Ni may promote Ni coarsening in the SOC operated in electrolysis cell mode. Defining the surface diffusivity as the product of the fractional surface coverage and single-molecule surface diffusivity, the surface diffusivity of NiH on Ni (111) is found to be large enough under a significant overpotential to support the above hypothesis. However, as NiH could dissociate on Ni (111), more work is needed to show that NiH may promote Ni coarsening.

Mantz, Yves↗

DOE OSTI2025

Irreversible oxygen poisoning: Modeling the initial water dissociation kinetics on δ-Pu(111) and δ-Pu(100) through density functional theory

In this work, the initial kinetics of water dissociation on two facets of δ-plutonium, δ-Pu(111) and δ-Pu(100), are explored through density functional theory in order to understand how water dissociation occurs on these facets. We explored the dissociation of water via the formation of hydroxyls, atomic hydrogen, and atomic oxygen species on each facet. We calculate low energetic barriers for water to split to adsorbed hydrogen and hydroxyl species at 0.19 eV for δ-Pu(111) and 0.07 eV for δ-Pu(100). The hydroxyl has a barrier of 0.64 and 0.37 eV to cleave the hydrogen–oxygen bond on δ-Pu(111) and δ-Pu(100), respectively. Due to the highly exergonic adsorption free energy of atomic oxygen of −2.10 eV, the metallic surfaces are found to be fully covered in oxygen, even with the inclusion of oxygen lateral interactions. When combined with the reaction thermodynamics, this free energy of adsorption forms a molecular oxygen desorption barrier greater than 9.51 eV (918 kJ/mol). These results, combined with simulated temperature programed desorption spectra, indicate that oxygen formed via water dissociation on δ-Pu (111) and (100) facets induces an irreversible poison that prevents further reaction of water directly on metallic plutonium surfaces, most likely due to the strong hybridization of the Pu and O valence states. Therefore, these results imply that another mechanism is responsible for the continuously experimentally measured water dissociation, which is most likely due to the Pu oxide.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Activation, Dehydrogenation, and Carbon–Carbon Coupling of Methane by Iridium Cations Studied by Infrared Multiple Photon Dissociation Spectroscopy and Density Functional Theory

Products resulting from the sequential activation of one, two, three, and four methane molecules by atomic iridium cations were characterized by gas-phase infrared multiple photon dissociation spectroscopy and density functional theory (DFT) calculations. Iridium cations were generated using a laser ablation source and reacted with methane in a linear radiofrequency ion trap before mass analysis and spectroscopic interrogation in a Fourier transform ion cyclotron resonance mass spectrometer coupled to the free-electron laser for intracavity experiments (FELICE) beamline. Product ions were irradiated using infrared light over the 250–1500 cm –1 range. Comparisons between the experimental and DFT-calculated spectra enabled structural determination of the products formed. The observed products include HIrCH + , Ir(CH 2 ) 2 + , H s Ir(C 3 H 5 ) + , and Ir(CH 3 ) s (C 3 H 5 ) + , where the subscript s denotes a syn orientation of the two ligands. Furthermore, formation of the latter two products provides evidence for efficient C–H bond activation and subsequent C–C coupling on the atomic iridium cation.

Chemical reactions↗

DOE OSTI2025

A noncollinear density functional theory ansatz for the phononic and thermodynamic properties of α -Pu

Plutonium's phase diagram is host to complex structures and interactions that make the description of its ground state properties elusive. Using all-electron density functional theory calculations, we study the thermodynamic properties of α-Pu. To do this, we build on recent work in the literature by introducing a novel noncollinear magnetic ansatz for α-Pu's ground state. We show that the noncollinear ansatz improves the description of the experimental phonon density of states, heat capacity, and thermal expansion over previous results obtained with a collinear ansatz. Here, despite lacking a bond-equivalent picture of the bonding, these new results for α-Pu, along with recent results on δ-Pu, demonstrate the efficacy of noncollinear ansatzes for the description of plutonium.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Voltage Probability Density Function Shaping Control Strategy Considering Grid Operational Uncertainties

It is well-known that power systems operation always affected by various uncertainties which make the bus voltage a random process that can be characterized by its probability density function (PDF) at any time instant. In this context, this paper presents a novel PDF-based voltage control framework for power systems. By modeling voltage as a stochastic process, we formulate a stochastic differential equationthat captures grid uncertainties. The associated Fokker-Planck-Kolmogorov equation is derived to describe the evolution of the voltage PDF, which enables the formulation of a PDF-shaping control strategy. To simplify the PDF control formulation, a B-spline neural network is introduced for real-time estimation and regulation of the voltage distribution. The proposed PDF control law updates voltage references for energy storage systems and synchronous generators using real-time PDF measurements and feedback signals. The proposed method is validated on a modified Kundur’s two-area system. Simulation results demonstrate that the controller can significantly improve the voltage stability under stochastic conditions, highlighting its effectiveness in modern inverter-rich grids.

Gui, Yonghao [ORNL] (ORCID:0000000250435534)↗

DOE OSTI2026

Insights on carbon dioxide adsorption in a flexible 1-D coordination polymer from in situ X-ray scattering and density functional theory

A combination of in situ small-angle X-ray scattering (SAXS) microstructure characterization over scales ranging from 1 nm to 10 µm and powder X-ray diffraction (XRD) structure characterization under various gas/pressure/temperature conditions with density functional theory (DFT) calculations provides new insights for the CO2 sorption behavior of a one-dimensional porous coordination polymer: catena-bis(dibenzoylmethanato)(4,4′-bipyridyl)nickel(II), denoted NiDBM-Bpy. The NiDBM-Bpy chains are held together by van der Waals forces, but the structure of guest-free NiDBM-Bpy is unsolved due to a lack of suitable crystals and high-quality powder XRD patterns. Nevertheless, SAXS and powder XRD can follow microstructural and structural changes as a function of gas pressure, composition and temperature. Both mixed-gas flow and static supercritical CO2 regimes are explored experimentally. DFT calculations are used to model the structural variation associated with XRD changes and hysteresis in the sorption isotherms. XRD and DFT calculations suggest that an orthorhombic Fddd structure with two CO2 per Ni emerges following a transition from the structure with lower molar volume and symmetry that exists without CO2 present.

Allen, Andrew J.↗

DOE OSTI2025

Extending Orbital-Optimized Density Functional Theory to L-Edge XPS and Beyond: Spin–Orbit Coupling via Nonorthogonal Quasi-Degenerate Perturbation Theory

Quantum mechanical calculations of core electron binding energies (CEBEs) are relevant to interpreting X-ray photoelectron spectroscopy (XPS). Orbital-optimized density functional theory (OO-DFT) accurately predicts K-edge CEBEs but is challenged by the presence of significant spin–orbit coupling (SOC) at L- and higher edges involving inner-shell orbitals with nonzero angular momentum. Here, to extend OO-DFT to L-edges and higher, our method utilizes scalar-relativistic, spin-restricted open-shell OO-DFT to construct a minimal, quasi-degenerate basis of core-hole states corresponding to a chosen inner-shell (e.g., ionizing all six possible 2p spin orbitals). Nonorthogonal configuration interaction (NOCI) is then used to obtain the matrix elements of the full Hamiltonian including SOC in this quasi-degenerate model space of determinants. Using a screened 1-electron SOC operator parametrized with the Dirac-Coulomb-Breit (DCB) Hamiltonian results in doublet splitting (DS) values for third row elements that are nearly in quantitative agreement with experiment. The resulting NOCI eigenvalues are shifted by the average of the (scalar) OO-DFT CEBEs to yield CEBEs (split by SOC) corrected for dynamic correlation. Comparing calculations on gas phase molecules with experimental results establishes that NO-QDPT with the SCAN functional (NO-QDPT/SCAN), using the DCB screened 1-electron SOC operator is accurate to about 0.2 eV for L-edge CEBEs of molecules containing third row atoms. However, this NO-QDPT approach becomes less accurate for fourth-row elements starting in the middle of the 3d transition metal series, with errors increasing as atomic number increases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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

Density Functional Tight-Binding Enables Tractable Studies of Quantum Plasmonics

Routine investigations of plasmonic phenomena at the quantum level present a formidable computational challenge due to the large system sizes and ultrafast time scales involved. This Feature Article highlights the use of density functional tight-binding (DFTB), particularly its real-time time-dependent formulation (RT-TDDFTB), as a tractable approach to study plasmonic nanostructures from a quantum mechanical purview. We begin by outlining the theoretical framework and limitations of DFTB, emphasizing its efficiency in modeling systems with thousands of atoms over picosecond time scales. Applications of RT-TDDFTB are then explored in the context of optical absorption, nonlinear harmonic generation, and plasmon-mediated photocatalysis. We demonstrate how DFTB can reconcile classical and quantum descriptions of plasmonic behavior, capturing key phenomena such as size-dependent plasmon shifts and plasmon coupling in nanoparticle assemblies. Lastly, we showcase DFTB’s ability to model hot carrier generation and reaction dynamics in plasmon-driven H2 dissociation, underscoring its potential to model photocatalytic processes. Collectively, these studies establish DFTB as a powerful, yet computationally efficient tool to probe the emergent physics of materials at the limits of space and time.

Electrical energy↗

DOE OSTI2025

Revealing symmetry-broken superconducting configurations by density functional theory

A coherent theory for the superconductivity of both conventional and unconventional superconductors is currently lacking. Here we show that superconductivity arises from the formation of a symmetry-broken superconducting configuration (SCC) due to atomic perturbation of the normal conducting configuration (NCC). This electron–phonon interaction creates straight one-dimensional tunnels (SODTs) for charge density of electrons and/or holes as revealed by the calculations based on density functional theory (DFT). The SODTs act as resistance-free superhighways and are correlated to the Cooper pairs in the Bardeen–Cooper–Schrieffer (BCS) theory. The formation of SODTs implies that the electron–phonon interaction in the BCS theory can be represented by the difference in charge densities between SCC and NCC predicted by DFT. The present work highlights that in conventional superconductors, SODTs are embedded within the bulk materials and are easily destroyed by phonon vibrations, resulting in a low critical superconducting temperature (T C ). Conversely, in unconventional superconductors such as YBa 2 Cu 3 O 7 (YBCO 7 ), SODTs are protected by a layered pontoon structure with very weak bonding to the bulk materials, maintaining SODTs’ stability at higher temperatures and leading to a much higher T C . The present approach is validated for 14 conventional superconductors of 18 pure elements and MgB 2 examined in this work, including the presently predicted superconductivity in Cu, Ag, Au, Sb, Bi, and MgB 2 at 0 K and 0 GPa, and one unconventional superconductor of YBCO 7 . Our discovery indicates that DFT can be a practical tool for predicting superconductors, enabling a systematic search for new superconducting materials in the future.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Strain-Induced Shifts in Defective Graphite Phonon Modes Predicted by Density Functional Theory

Carbon fiber composites have gained attention as a structural material because of their high strength-to-weight ratio, and understanding the effect of defects on reactivity and mechanical properties is important for the longevity and safety of the composite. Although it is known that strain causes the underlying graphitic vibrational modes to redshift, it is not clear how strain may alter reactivity and defect-induced vibrational changes. To investigate the strain-induced phonon changes of defective carbon fiber composites, density functional theory calculations of graphite are used, including intercalated hydrogen and fluorine defects. By comparing changes in the bond lengths, formation energies, and phonon density of states for uniaxially and biaxially strained graphite, strain was found to generally make defect formation more favorable and the specific behavior changes are dependent on the strain direction and defect identity. Specifically, intercalated fluorine phonons are more sensitive to strain than hydrogen intercalation phonons, and strain applied along the zigzag direction alters the calculated properties more than strain along the armchair direction. Finally, these results highlight the importance of understanding the microstructural effect of deviations from the ideal material because small changes in strain or defect type can significantly alter the behavior of the carbon fiber composite core.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Real-Space Constrained Density Functional Theory Investigation of Site-Specific, Interfacial Charge Recombination Dynamics Across the Au Nanoparticle/TiO 2 Heterojunction

Au nanoparticle (NP)/TiO 2 heterojunction is a representative system to study interfacial charge transfer in photocatalysis and photovoltaics, where suppressing recombination from TiO 2 to Au can enhance hot carrier extraction. We apply real-space constrained density functional theory (CDFT) with Marcus theory to quantify charge recombination time scales across Au/TiO 2 . This approach enables direct control and visualization of charge-separated states, aligning with site-specific probes like time-resolved X-ray photoelectron spectroscopy (trXPS). We find that the charge-separated state features a bipolaron, with recombination dominated by TiO 2 LUMO to Au HOMO transitions, primarily at interfacial Au sites. Marcus rate predictions are benchmarked with surface hopping methods, quantifying differences in time scales and computational efficiency. Lastly, we examine how the Au cluster size affects the free energy change (ΔG) and reorganization energy (λ), explaining trends in closed-shell systems and highlighting challenges for open-shell extrapolations. Overall, CDFT + Marcus theory provides efficient, mechanistically transparent interfacial charge transfer modeling, and we clearly defined its applicability and limitation.

Glenna, Drew M. [Univ. of Idaho, Idaho Falls, ID (↗

DOE OSTI2026

Excitation energy of fission fragments within nuclear time-dependent density functional theory

The number and properties of the neutrons and photons emitted in nuclear fission are directly related to the excitation energy of the fission fragments when they are formed at scission. Though not observable experimentally because of the extremely short timescales, the excitation energy of fission fragments can be predicted by microscopic theory based on time-dependent density functional theory (TDDFT). Initial results on the value of the total kinetic energy of fission reactions were very promising, but could not probe all possible fragmentations. Here, in this work, we perform large-scale TDDFT calculations in 240 Pu enabled by the development of a new TDDFT solver. We obtain TDDFT trajectories covering nearly all possible fragmentations. We find that the total kinetic energy is close to experimental values only for the most likely fission while it is severely underestimated at both small and large asymmetries. This conclusion seems rather independent of the parametrization of the energy functional, both in its particle-hole and particle-particle channels.

Physics - Nuclear physics and radiation physics↗

DOE OSTI2025

Development of a Transferable Density-Functional Tight-Binding Model for Organic Molecules at the Water/Platinum Interface

A computationally efficient and transferable approach for modeling reactions at metal/water interfaces could significantly accelerate our understanding and ultimately the development of new catalytic transformations, particularly in the context of the emerging field of biomass conversion. Here, we present a parametrization of Pt–X (X = H, O, C) density-functional tight-binding (DFTB) for addressing this need. We first constructed Pt–H, Pt–O, and Pt–C repulsive potential splines. These pairwise parameters were then augmented to include many-body interactions using the Chebyshev Interaction Model for Efficient Simulation (ChIMES). We compare the geometrical and energetic performances of both DFTB and DFTB/ChIMES methods with DFT reference data across a variety of organic molecules at the platinum surface from nanoparticles to single-crystal surfaces. DFTB shows limited transferability between extended crystal surfaces and small nanoparticles. This transferability is significantly improved through the introduction of three-body interactions with Pt in DFTB/ChIMES, which provides consistent results across various systems, with reductions in the RMSD from around 30 kcal/mol in DFTB to around 10 kcal/mol. We demonstrate the stability and reliability of the obtained parameters by performing metadynamic simulations for the adsorption of phenol on Pt(111). We observe that DFTB itself is undersolvating the surface, leading to only one or two chemisorbed water molecules in a c(4 × 6) unit cell. In contrast, DFTB/ChIMES leads to a coverage of about 0.5 ML and successfully captures the chemisorbed mode of phenol at both the solid/liquid and the solid/gas interfaces. Furthermore, in agreement with experimental measurements, the adsorption at the solid/liquid interface is significantly weaker than that at the solid/gas interface. As a result, we highlight that even with DFTB, where we can accumulate dynamics for more than 1 ns for a given system, the simulations are not fully converged.

Adsorption↗

DOE OSTI2026

Development of an ab initio learned model of electron deposition range in deuterium-tritium plasmas through time-dependent density functional theory calculations and machine learning

Accurate hydrodynamic modeling for laser-direct-drive (LDD) inertial-confinement-fusion (ICF) relies on precise calculations of the electron thermal conduction in all target materials. The nonlocal stopping range of electrons in ICF plasmas directly influences thermal conduction; yet, no first principles model exists for the electron mean free path in the conduction-zone regime. This work utilized time-dependent stochastic density-functional theory (TD-sDFT) to calculate the electron stopping power in deuterium-tritium (DT) plasmas at (ρ, T) conditions relevant to the conduction zone and the compressed shell in ICF. Using a combination of our TD-sDFT data and already established analytical models, we developed and trained an artificial neural network to create a global model for the nonlocal electron deposition range, λ E . We compared our machine-learning (ML) based model for λ E to the currently-used modified-Lee-More model in LDD radiation-hydrodynamic codes, such as lilac, and saw an overall decrease in the deposition range. To understand the effects of λ E on LDD ICF implosion dynamics, we implemented the ML-based model into lilac; specifically, we looked at designs consistent with a current experiment on the OMEGA laser and for a newly designed LDD-ICF target for the future OMEGA-Next facility. In both cases, we saw an overall drop in predicted ablation pressure, peak areal density, and neutron yield due to the reduced thermal conduction (smaller λ E ) in DT plasmas. Comparisons with the experiment on OMEGA are also made.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Creation of an Fe 3 P Schreibersite Density Functional Tight Binding Model for Astrobiological Simulations

The mineral schreibersite, e.g., Fe 3 P, is commonly found in iron-rich meteorites and could have served as an abiotic phosphorus source for prebiotic chemistry. However, atomistic calculations of its degradation chemistry generally require quantum simulation approaches, which can be too computationally cumbersome to study sufficient time and length scales for this process. In this regard, we have created a computationally efficient semiempirical quantum density functional tight binding (DFTB) model for iron and phosphorus-containing materials by adopting an existing semiautomated workflow that represents many-body interactions by linear combinations of Chebyshev polynomials. Here, we have utilized a relatively small training set to optimize a DFTB model that is accurate for schreibersite physical and chemical properties, including its bulk properties, surface energies, and water absorption. We then show that our model shows strong transferability to several iron phosphide solids as well as multiple allotropes of iron metal. Our resulting DFTB parametrization will allow us to interrogate schreibersite aqueous decomposition at longer time and length scales than standard quantum approaches, providing for more detailed investigations of its role in prebiotic chemistry on early Earth.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Gold-Standard Chemical Database 137 (GSCDB137): A Diverse Set of Accurate Energy Differences for Assessing and Developing Density Functionals

We present GSCDB137, a rigorously curated benchmark library of 137 data sets (8377 entries) covering main-group and transition-metal reaction energies and barrier heights, (intra- and intermolecular) noncovalent interactions, dipole moments, polarizabilities, electric-field response energies, and vibrational frequencies. Legacy data from GMTKN55 and MGCDB84 have been updated to today's best reference values; redundant or low-quality points were removed, and many new, property-focused sets were added. Testing 29 popular density functional approximations (DFAs) confirms the expected Jacob's-ladder hierarchy overall but also reveals notable exceptions: functional performance for frequencies and electric-field properties correlates poorly with that for other ground-state energetics. ωB97M-V and ωB97X-V are the most balanced hybrid meta-GGA and hybrid GGA, respectively; B97M-V and revPBE-D4 lead the meta-GGA and GGA classes. Double hybrids lower mean errors by about 30% versus their hybrid analogues but demand careful frozen-core, basis set, and spin contamination treatment. GSCDB137 offers a comprehensive, openly documented platform for rigorous validation of DFA and universal machine learning potentials, and training of the next generation of exchange-correlation functionals.

Liang, Jiashu [University of California, Berkeley,↗

DOE OSTI2024

Second-harmonic generation tensors from high-throughput density-functional perturbation theory

Optical materials play a key role in enabling modern optoelectronic technologies in a wide variety of domains such as the medical or the energy sector. Among them, nonlinear optical crystals are of primary importance to achieve a broader range of electromagnetic waves in the devices. However, numerous and contradicting requirements significantly limit the discovery of new potential candidates, which, in turn, hinders the technological development. In the present work, the static nonlinear susceptibility and dielectric tensor are computed via density-functional perturbation theory for a set of 579 inorganic semiconductors. The computational methodology is discussed and the provided database is described with respect to both its data distribution and its format. Several comparisons with both experimental and ab initio results from literature allow to confirm the reliability of our data. The aim of this work is to provide a relevant dataset to foster the identification of promising nonlinear optical crystals in order to motivate their subsequent experimental investigation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗