Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Cross-slip and easy-glide CRSS in titanium: Theoretical predictions and in-situ TEM measurements

This study investigates the mechanics of prismatic and first-order pyramidal $\langle$a$\rangle$ slip in titanium (Ti), elucidating the physics of easy-glide and cross-slip through a combination of theory and experiments. Screw-character prismatic (Pr) dislocations in Ti are of particular interest because their complex cores can be stable or unstable, leading to activation by either cross-slip or planar glide. To investigate these mechanisms, site-specific micro tensile samples were prepared using focused ion beam (FIB) milling and mounted on a push-to-pull (PTP) device for in-situ transmission electron microscopy (TEM) tensile testing. The in-situ experiments provide direct observations of the onset of dislocation motion and the precise determination of the critical resolved shear stress (CRSS) for the activated mechanisms, and its evolution with load cycling. A comprehensive theory has been developed to predict the CRSS values for easy glide, cross-slip, and multiplication of dislocations. Predicted critical stresses for pyramidal (π)-to-Pr and reverse cross-slip agree closely with the experimental measurements. The latter cross-slip stress is a factor of two higher than that of unobstructed planar slip. The model accounts for overlapping dislocation cores and employs a Wigner-Seitz based cell to evaluate misfit energies. By combining ab initio density functional theory (DFT) with anisotropic elasticity, the framework identifies minimum energy pathways for dislocation glide, which can be intermittent and zig-zag. A simplified expression utilizing (π) and Pr Schmid factor ratios is proposed for critical stress corresponding to (π)-to-Pr cross-slip transition. The results are strongly dependent on crystal orientation, underscoring non-Schmid behavior. Overall, this study explores key critical stress parameters essential for informing higher-scale simulations of plasticity in Ti.

Critical stress↗
DOE OSTI2023

Can a deep-learning model make fast predictions of vacancy formation in diverse materials?

The presence of point defects, such as vacancies, plays an important role in materials design. Here, we explore the extrapolative power of a graph neural network (GNN) to predict vacancy formation energies. We show that a model trained only on perfect materials can also be used to predict vacancy formation energies (E vac ) of defect structures without the need for additional training data. Such GNN-based predictions are considerably faster than density functional theory (DFT) calculations and show potential as a quick pre-screening tool for defect systems. To test this strategy, we developed a DFT dataset of 530 E vac consisting of 3D elemental solids, alloys, oxides, semiconductors, and 2D monolayer materials. We analyzed and discussed the applicability of such direct and fast predictions. We applied the model to predict 192 494 E vac for 55 723 materials in the JARVIS-DFT database. Our work demonstrates how a GNN-model performs on unseen data.

2D materials↗
DOE OSTI2023

Single-Step Selective Oxidation of Methane by Iron-Oxo Species in the Metal–Organic Framework MFU-4l

The direct and selective conversion of methane to methanol can be considered a holy grail for catalysis research. In this work, we study a metal-organic framework known as MFU-4l, modified by design to include highly reactive iron-oxo species for the catalytic C-H bond activation of methane. We investigate the oxidation of methane and the further potential oxidation of the product methanol using N 2 O as an oxygen source and map the potential energy landscape of these reactions using density functional theory calculations. We show that the highest energy barrier encountered during the methane oxidation process is not the C-H bond breaking, but the activation of the iron center by N 2 O. Furthermore, the potential energy landscape for the C-H bond activation exhibits a large, high-energy plateau region instead of a sharp transition state, thus differing from the traditional radical rebound mechanism. Furthermore, this insight offers interesting potential routes to enhance the catalytic activity of the catalyst, to hinder unwanted deactivation pathways, and to reduce the activity towards the over-oxidation of the product.

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

Extending quantum-mechanical benchmark accuracy to biological ligand-pocket interactions

Predicting the binding affinity of ligands to protein pockets is key in the drug design pipeline. The flexibility of ligand-pocket motifs arises from a range of attractive and repulsive electronic interactions during binding. Accurately accounting for all interactions requires robust quantum-mechanical (QM) benchmarks, which are scarce for ligand-pocket systems. Additionally, disagreement between “gold standard” Coupled Cluster (CC) and Quantum Monte Carlo (QMC) methods casts doubt on many benchmarks for larger non-covalent systems. We introduce the “QUantum Interacting Dimer” (QUID) benchmark framework containing 170 non-covalent (non-)equilibrium systems modeling chemically and structurally diverse ligand-pocket motifs. Symmetry-adapted perturbation theory shows that QUID broadly covers non-covalent binding motifs and energetic contributions. Robust binding energies are obtained using complementary CC and QMC methods, achieving agreement of 0.5 kcal/mol. The benchmark data analysis reveals that several dispersion-inclusive density functional approximations provide accurate energy predictions, though their atomic van der Waals forces differ in magnitude and orientation. Contrarily, semiempirical methods and empirical force fields require improvements in capturing non-covalent interactions (NCIs) for out-of-equilibrium geometries. The wide span of NCIs, highly accurate interaction energies, and analysis of molecular properties take QUID beyond the “gold standard” for QM benchmarks of ligand-protein systems.

Puleva, Mirela [University of Luxembourg, Luxembou↗
DOE OSTI2023

2023 Roadmap on molecular modelling of electrochemical energy materials

New materials for electrochemical energy storage and conversion are the key to the electrification and sustainable development of our modern societies. Molecular modelling based on the principles of quantum mechanics and statistical mechanics as well as empowered by machine learning techniques can help us to understand, control and design electrochemical energy materials at atomistic precision. Therefore, this roadmap, which is a collection of authoritative opinions, serves as a gateway for both the experts and the beginners to have a quick overview of the current status and corresponding challenges in molecular modelling of electrochemical energy materials for batteries, supercapacitors, CO 2 reduction reaction, and fuel cell applications.

25 ENERGY STORAGE↗
DOE OSTI2023

Decagonal Sn clathrate on d-Al-Ni-Co

Decagonal quasiperiodic ordering of Sn thin films on d-Al-Ni-Co is shown based on scanning tunneling microscopy (STM), low-energy electron diffraction, and density functional theory (DFT). Interestingly, the decagonal structural correlations are partially retained even up to a large film thickness of 10 nm grown at 165±10K. The nucleation centers called “Sn white flowers” identified by STM at submonolayer thickness are recognized as valid patches of the decagonal clathrate structure with low adsorption energies of these motifs. Due to the excellent lattice matching (to within 1%) between columns of Sn dodecahedra in the clathrate structure and pentagonal motifs at the d-Al-Ni-Co surface, the interfacial energy favors clathrate over the competing Sn crystalline forms. DFT study of the Sn/Al-Ni-Co composite model shows good mechanical stability, as shown by the work of separation of Sn from Al-Ni-Co slab that is comparable to the clathrate self-separation energy. Finally, the relaxed surface terminations of the R 2 T 4 clathrate approximant are in self-similarity correspondence with the motifs observed in the STM images from monolayer to the thickest Sn film.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Thermodynamics of the Metal Carbonates and Bicarbonates of Mn, Co, Ni, Cu, and Zn Relevant to Mineral Energetics

The gas phase heats of formation of ground-state MCO 3 , M(HCO 3 ) 2 , and M(HCO 3 )(OH), where M = Mn, Co, Ni, Cu, and Zn, have been predicted using the correlated molecular orbital theory at the CCSD(T) level extrapolated to the complete basis set limit using the Feller–Peterson–Dixon (FPD) approach. Cohesive energies of the carbonates were predicted based on the calculated gas phase and experimental solid heats of formation. Coulombic dissociation energies (CDEs) between metal cations and anions show a near-linear correlation with Shannon metal cation atomic radii, yet no correlation is found with the hardness of these cations. The total reaction dissociation energies (TRDEs) of transition metals are higher than their CDEs for the di-bicarbonates, in contrast to those for Mg and Ca based on our prior work. In addition to differences in the energies needed to prepare the transition metal dications, electron donation from the ligands to the 3d orbitals of open-shell transition metal dications from lone pairs of adjacent O atoms also plays a role. No electron donation from the ligands to the fully occupied 3d orbitals of Zn and Cd was found. Decomposition energies for generating MO, CO 2 , and/or H 2 O were calculated. Gas phase metal exchange energies only partially correlate with the electrochemical series for M(s) → M 2+( aq). The FPD heats of formation were used to benchmark a range of density functional theory exchange–correlation functionals, including those commonly used in solid-state mineral calculations. None of the functionals provided chemical accuracy agreement (±1 kcal/mol) with the FPD results. So, the best agreement with the FPD results is predicted for the $\tau$-HCTH functional with an average unsigned error of 8.3 kcal/mol.

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

Continuum radiated power density $(P^{\textrm{C}}_{\textrm{rad}})$ and effective charge $(Z_{\textrm{eff}})$ estimates from multi-energy photon-counting measurements

Multi-energy soft x-ray pinhole cameras have been designed, built, calibrated, and operated at Madison symmetric Torus, Alcator C-Mod, and more recently at Tokamak a configuration variable and Tungsten Environment in Steady-State Tokamak (WEST), to measure plasma emission across multiple energy ranges. Here we describe a new methodology to estimate the local continuum radiated power density and the plasma effective charge (Z eff ) directly from photon-counting measurements of the line-free continuum emission (Bremsstrahlung and Radiative Recombination) in several energy bands between 11 and 18 keV. This capability is particularly valuable for confinement systems using metal plasma-facing components, where x-ray losses from interactions with the sputtered wall can represent a significant fraction of the total radiated power (P rad ). The approach leverages a well-characterized detector responsivity, modeled by a complementary error function, and interprets the differential multi-energy measurements between adjacent energy levels through the probability density function of a Gaussian distribution. The implementation of this diagnostic technique is currently under development on the WEST tokamak, aiming at the goal of providing real-time P rad and Z eff measurements during long-pulse operation (up to 1000 s) in the 2026 campaign.

WEST↗
DOE OSTI2022

Theoretical Modeling of Electrochemical Proton-Coupled Electron Transfer

Proton-coupled electron transfer (PCET) plays an essential role in a wide range of electrocatalytic processes. Here, a vast array of theoretical and computational methods have been developed to study electrochemical PCET. These methods can be used to calculate redox potentials and pKa values for molecular electrocatalysts, proton-coupled redox potentials and bond dissociation free energies for PCET at metal and semiconductor interfaces, and reorganization energies associated with electrochemical PCET. Periodic density functional theory can also be used to compute PCET activation energies and perform molecular dynamics simulations of electrochemical interfaces. Various approaches for maintaining a constant electrode potential in electronic structure calculations and modeling complex interactions in the electric double layer (EDL) have been developed. Theoretical formulations for both homogeneous and heterogeneous electrochemical PCET spanning the adiabatic, nonadiabatic, and solvent-controlled regimes have been developed and provide analytical expressions for the rate constants and current densities as functions of applied potential. The quantum mechanical treatment of the proton and inclusion of excited vibronic states have been shown to be critical for describing experimental data, such as Tafel slopes and potential dependent kinetic isotope effects. The calculated rate constants can be used as input to microkinetic models and voltammogram simulations to elucidate complex electrocatalytic processes.

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

Enhancing Cycling Stability and Capacity Retention of NMC811 Cathodes by Reengineering Interfaces via Electrochemical Fluorination

High-capacity cathodes (LiNi 0.8 Mn 0.1 Co 0.1 O 2 ) that can boost the energy density of lithium-ion batteries are promising candidates for vehicle electrification. However, several factors specific to high energy density materials entailing electrode reactions inhibit their application. Fluorination has shown a promising ability to combat the detrimental electrochemical performances of cathode materials, however, it remains difficult to achieve the desired functionality. In this work, a novel electrochemical fluorination (ECF) that demonstrates a promising electrochemical performance enhancement via stabilization of the cathode–electrolyte-interphase (CEI) by forming conformal LiF is proposed. Besides LiF surface layer formation, ECF reduces the degree of fluorination-induced Ni/Li disordering and enhances the layered structural stability as probed by X-ray diffraction. Because of the robust CEI, ECF-NMC811 cathodes deliver 203.0 mAh g -1 first discharge capacity at the current rate of C/10, with ≈98% capacity retention up to 100 cycles. Similarly, it delivers ≈180 mAh g -1 capacity at a 1 C rate with 86.4% capacity retention up to 200 cycles with average coulombic efficiency of > 99.5%. Comprehensive characterization with a multitude of probes reveals that ECF enhances the cycling stability of the electrode without altering bulk structure and morphology.

25 ENERGY STORAGE↗
DOE OSTI2022

Systematic Evaluation of Counterpoise Correction in Density Functional Theory

A widespread belief persists that the Boys–Bernardi function counterpoise (CP) procedure “overcorrects” supramolecular interaction energies for the effects of basis-set superposition error. To the extent that this is true for correlated wave function methods, it is usually an artifact of low-quality basis sets. The question has not been considered systematically in the context of density functional theory, however, where basis-set convergence is generally less problematic. We present a systematic assessment of the CP procedure for a representative set of functionals and basis sets, considering both benchmark data sets of small dimers and larger supramolecular complexes. The latter include layered composite polymers with ~150 atoms and ligand–protein models with ~300 atoms. Provided that CP correction is used, we find that intermolecular interaction energies of nearly complete-basis quality can be obtained using only double-ζ basis sets. Furthermore, this is less expensive as compared to triple-ζ basis sets without CP correction. CP-corrected interaction energies are less sensitive to the presence of diffuse basis functions as compared to uncorrected energies, which is important because diffuse functions are expensive and often numerically problematic for large systems. Our results upend the conventional wisdom that CP “overcorrects” for basis-set incompleteness. In small basis sets, CP correction is mandatory in order to demonstrate that the results do not rest on error cancellation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Comparing ab initio and quantum-kinetic approaches to electron transport in warm dense matter

Accurate knowledge of the electronic transport properties of warm dense matter is one of the main concerns of research in high-energy-density physics. Three modern approaches with vastly different levels of fidelity are reviewed and compared: the Kubo–Greenwood (KG) approach based on density-functional-theory molecular dynamics simulations (QMD), quantum kinetic theory based on average-atom models, and time-dependent density functional theory. Throughout, emphasis is placed on the connection between static properties of the electrons (e.g., density of states) and transport properties. Overall, it is found that whenever the conduction electrons can be modeled as being nearly free, fair to excellent agreement is found between QMD and kinetic theory approaches. Such a circumstance is required for modeling warm dense matter as a plasma of ions and free electrons, which is assumed in most kinetic theory approaches. The sensitivity of transport properties to the electronic structure is further highlighted by comparing different exchange–correlation approximations in QMD and KG calculations. It is found that the inclusion of exact exchange via thermal hybrid functionals can make a pronounced impact on electrical and thermal conduction in warm dense matter. We also investigate dynamic screening physics via kinetic theory and time-dependent density functional theory calculations of the mean free path of an electron in a hot dense plasma. In sum, we identify three axes along which to make progress in predicting electron transport in warm dense matter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Importance of finite-size corrections for accurate ab initio modeling of carrier capture at semiconductor defects: A case study of substitutional C N in GaN

In ab initio studies of carrier-capture processes in defective semiconductor materials, the single-effective-mode formalism and the static-coupling approximation have become the predominant theoretical approaches for determining carrier-capture coefficients. The single-mode formalism relies on accurate nonequilibrium defect energies obtained from density-functional theory (DFT), where required inputs are a series of configurationally displaced, defect-containing supercells obtained using an interpolative ansatz, and where the DFT outputs are corresponding total energies that have traditionally been postprocessed using a long-established ground-state formulation of finite-size corrections and defect-formation energies. This formulation remains commonly used even though the defects that form a configuration-coordinate (CC) diagram typically exist as structures that are displaced from the ground state. To remedy this inconsistency, Kumagai has recently proposed novel methods for implementing finite-size corrections specifically intended for DFT calculations of the defect energies used to construct CC diagrams and implement the single-mode formalism [Y. Kumagai, Phys. Rev. B 107, L220101 (2023)]. Kumagai's approach builds on the latest finite-size-correction methods introduced to describe vertical charge-state transitions for charge-localizing point defects in semiconductors and insulators [T. Gake et al., Phys. Rev. B 101, 020102 (2020); S. Falletta et al., Phys. Rev. B 102, 041115 (2020)]. The newly identified finite-size artifact treated in these studies is the polarization charge induced on a configurationally frozen defect and its subsequent interaction with a vertical transition in charge state. In this work, we evaluate Kumagai's proposed methodology by applying it in a high-precision DFT study of carrier capture by substitutional C N in GaN, a well-characterized and technologically relevant defect and material. We have rigorously calculated C N defect energies across various supercell sizes for each defect configuration and charge state on the hole-capture CC diagram of C N (𝑞=−1), enabling a direct comparison of the slopes of the defect energies versus inverse cell size with those predicted by Kumagai. The most consequential prediction of Kumagai's method is that these slopes distinctly vary as the square of the linear-interpolation parameter used to construct the nonequilibrium defect configurations. Our results quantitatively support this prediction. Moreover, with these new finite-size corrections and multiple-cell-size DFT calculations in place, we find that the classical energy barrier for hole capture by C N (𝑞=−1) in GaN decreases to 0.092–0.127 eV. This finding confirms the recent ≈ 0.1 eV prediction of Reshchikov based on the weak temperature dependence for hole capture observed in photoluminescence experiments [M. A. Reshchikov, J. Appl. Phys. 129, 121101 (2021)]. These results stand in stark contrast to previously calculated barriers of 0.486 and 0.73 eV, which also used the single-mode formalism but were obtained by instead using ground-state-based finite-size corrections. Our reduced classical barrier for capture increases the temperature-dependent hole-capture coefficient of a C N (𝑞=−1) defect by more than two to four orders of magnitude for temperatures of 100–600 K, compared to the previous 0.486 eV results. While other defects may not be as dramatically affected as here, we suggest that incorporating proper finite-size corrections for the vertical-transition-like states embedded within CC diagrams is an essential, yet previously unrecognized, component of accurate modeling of carrier-capture when using the single-effective-mode formalism.

dielectric properties↗
DOE OSTI2026

Characterizing Defect Dynamics in Silicon Carbide Using Symmetry-Adapted Collective Variables and Machine Learning Interatomic Potentials

Silicon carbide (SiC) divacancies are attractive candidates for spin-defect qubits possessing long coherence times and optical addressability. The high activation barriers associated with SiC defect formation and motion pose challenges for their study by first-principles molecular dynamics. In this work, we develop and deploy machine learning interatomic potentials (MLIPs) to accelerate defect dynamics simulations while retaining ab initio accuracy. We employ an active learning strategy comprising symmetry-adapted collective variable discovery and enhanced sampling to compile configurationally diverse training data, calculation of energies and forces using density functional theory (DFT), and training of an E(3)-equivariant MLIP based on the Allegro model. Here, the trained MLIP reproduces DFT-level accuracy in defect transition activation free energy barriers, enables the efficient and stable simulation of multidefect 216-atom supercells, and permits an analysis of the temperature dependence of defect thermodynamic stability and formation/annihilation kinetics to propose an optimal annealing temperature to maximally stabilize VV divacancies.

Computer simulations↗
DOE OSTI2025

Structural Evidence of Interanionic Hydrogen Bonding in Phosphoric Acid Solutions

Interanionic hydrogen bonding (IAHB) is a noncovalent interaction between like-charged ions that challenges conventional electrostatic understanding. This study provides direct structural evidence of IAHB in concentrated aqueous phosphoric acid (PA) solutions, which exhibit >60% dissociation under these conditions. Oxygen K-edge X-ray absorption fine structure spectroscopy, combined with electron affinity time-dependent density functional theory calculations, reveals the formation of stable, cyclic phosphate-phosphate IAHB dimers at PA concentrations ≥7 M. Extended X-ray absorption fine structure data show distinct long-range ordering consistent with these dimers, and near-edge X-ray absorption fine structure spectra confirm a concentration-dependent transition from monomeric to dimeric species. Energy decomposition analysis through density functional theory shows that the formation of solution-phase IAHB is energetically favored and is attributed to polarization of and the charge transfer between the two fragments driven by the surrounding solvent molecules, in addition to the permanent electrostatics. These findings offers crucial structural insights into the H-bonded networks in concentrated PA, highlighting the critical role of solvent in facilitating anion–anion association.

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

Spin-state gaps and self-interaction-corrected density functional approximations: Octahedral Fe(II) complexes as case study

Accurate prediction of a spin-state energy difference is crucial for understanding the spin crossover phenomena and is very challenging for density functional approximations, especially for local and semi-local approximations due to delocalization errors. Here, we investigate the effect of the self-interaction error removal from the local spin density approximation (LSDA) and Perdew–Burke–Ernzerhof generalized gradient approximation on the spin-state gaps of Fe(II) complexes with various ligands using recently developed locally scaled self-interaction correction (LSIC) by Zope et al. [J. Chem. Phys. 151, 214108 (2019)]. The LSIC method is exact for one-electron density, recovers the uniform electron gas limit of the underlying functional, and approaches the well-known Perdew–Zunger self-interaction correction (PZSIC) as a particular case when the scaling factor is set to unity. Our results, when compared with reference diffusion Monte Carlo results, show that the PZSIC method significantly overestimates spin-state gaps favoring low spin states for all ligands and does not improve upon density functional approximations. The perturbative LSIC-LSDA using PZSIC densities significantly improves the gaps with a mean absolute error of 0.51 eV but slightly overcorrects for the stronger CO ligands. Finally, the quasi-self-consistent LSIC-LSDA, such as coupled-cluster single double and perturbative triple [CCSD(T)], gives a correct sign of spin-state gaps for all ligands with a mean absolute error of 0.56 eV, comparable to that of CCSD(T) (0.49 eV).

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

From Intermolecular Interaction Energies and Observable Shifts to Component Contributions and Back Again: A Tale of Variational Energy Decomposition Analysis

Quantum chemistry in the form of density functional theory (DFT) calculations is a powerful numerical experiment for predicting intermolecular interaction energies. However, no chemical insight is gained in this way beyond predictions of observables. Energy decomposition analysis (EDA) can quantitatively bridge this gap by providing values for the chemical drivers of the interactions, such as permanent electrostatics, Pauli repulsion, dispersion, and charge transfer. These energetic contributions are identified by performing DFT calculations with constraints that disable components of the interaction. This review describes the second-generation version of the absolutely localized molecular orbital EDA (ALMO-EDA-II). The effects of different physical contributions on changes in observables such as structure and vibrational frequencies upon complex formation are characterized via the adiabatic EDA. Example applications include red- versus blue-shifting hydrogen bonds; the bonding and frequency shifts of CO, N 2 , and BF bound to a [Ru(II)(NH 3 ) 5 ] 2 + moiety; and the nature of the strongly bound complexes between pyridine and the benzene and naphthalene radical cations. Additionally, the use of ALMO-EDA-II to benchmark and guide the development of advanced force fields for molecular simulation is illustrated with the recent, very promising, MB-UCB potential.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗