Search NASA⌕ Search

SEARCH · Search NASA

Results for “optimized energy 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.

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

At least 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Optimization of a Secondary Air Injector for a Rich-Quench-Lean (RQL) Ammonia Combustor Using Computational Fluid Dynamics

Ammonia has emerged as a promising medium for moving hydrogen around the globe due to its energy density, existing infrastructure for production, transportation, and storage, and multiple applications from fertilizer to power generation. However, one of the most significant challenges with ammonia combustion for power generation is the formation of nitrogen oxides (NOx) during combustion. Several combustion strategies have been developed to minimize the formation of NOx. One such strategy, known as rich-quench-lean (RQL), is a method that combusts NH3 in a fuel-rich environment, followed by a quick mix section and a lean burnout section. Rapid mixing of secondary air before lean burnout is thought to be important to minimize the formation of NOx. A genetic algorithm (GA) is used to parametrically vary the secondary air injector design, including the diameter, count, and angles over a constrained design space. The designs are evaluated using a non-reacting OpenFOAM model, with the objective function being the uniformity index of the secondary air (modeled as a scalar). Optimization campaigns show that traditional correlation-based designs might not be adequate. After running 414 OpenFOAM models, the optimal design is 32, 1 mm diameter, angled air injectors, increasing the uniformity index at the outlet of the quick mix zone by 36.8%. The most promising designs will be manufactured and tested in a small RQL burner setup, which is expected to lead to validation and insight into the practical usage of NH3 combustion for power generation.

ammonia combustion↗
DOE OSTI2024

Crystal structure of anthraquinone-2-carboxylic acid, C 15 H 8 O 4

The crystal structure of anthraquinone-2-carboxylic acid has been solved and refined using synchrotron X-ray powder diffraction data, and optimized using density functional theory techniques. Anthraquinone-2-carboxylic acid crystallizes in space groupP-1 (#2) witha= 3.7942(2),b= 13.266(5),c= 22.835(15) Å,α= 73.355(30),β= 89.486(6),γ= 86.061(1)°,V= 1098.50(7) Å 3 , andZ= 4. The crystal structure contains two independent molecules of anthraquinone-2-carboxylic acid. Although the expected hydrogen-bonded dimers are present, the dimers are not centrosymmetric. The dimer contains one molecule of each planar low-energy conformation. The crystal structure consists of a herringbone array of centrosymmetric pairs of molecules parallel to thebc-plane. The molecules stack along the shorta-axis. The powder pattern has been submitted to ICDD® for inclusion in the Powder Diffraction File™ (PDF®).

Materials Science↗
DOE OSTI2026

Neural network potentials with effective charge separation for non-equilibrium dynamics of ionic solids: a ZnO case study

Developing neural network potentials (NNPs) accurate under non-equilibrium dynamics is challenging, as such systems require extensive sampling beyond equilibrium phases. Here we construct high-fidelity NNPs for zinc oxide (ZnO), a polymorphic ionic solid, using density functional theory (DFT) reference data. To efficiently capture transitional configurations, we combine enhanced-sampling molecular dynamics with empirical potentials, data distillation, and pretraining on short-range atomic energies (A-Train), followed by transfer learning with DFT-relabeled datasets. This hierarchical approach improves transferability across polymorphs and stress states. We further introduce effective charge separation, treating long-range Coulombic terms analytically while short-range residual interactions are learned by the NNP. The optimal effective charges fall in the range 0.5–1.0 q e , consistent with dielectric-screened values derived from formal charges but distinct from Bader estimates. Motivated by this observation, we propose a simple data-driven protocol in which effective charges are optimized by comparing DFT reference energies with explicit Coulomb calculations, without additional NNP training. This strategy improves accuracy and transferability in DFT-level predictions of energies, forces, and stress. Together, these results provide a practical charge-selection framework for robust NNP development in ionic solids, enabling reliable simulation of polymorphic phase transformations and non-equilibrium dynamics.

Chemistry↗
DOE OSTI2024

Crystal structure of ractopamine hydrochloride, C 18 H 24 NO 3 Cl

The crystal structure of ractopamine hydrochloride has been solved and refined using synchrotron X-ray powder diffraction data, and optimized using density functional theory techniques. Ractopamine hydrochloride crystallizes in space groupPbca(#61) witha= 38.5871(49),b= 10.7691(3),c= 8.4003(2) Å,V= 3490.75(41) Å 3 , andZ= 8. The ractopamine cation contains two chiral centers, and the sample consists of a mixture of the S,S/R,R/S,R and R,S forms. Models for the two diastereomers S,S and S,R were refined, and yielded equivalent residuals, but the S,R form is significantly lower in energy. The crystal structure consists of layers of molecules parallel to thebc-plane. In each structure one of the H atoms on the protonated N atom acts as a donor in a strong discrete N–H⋯Cl hydrogen bond. Hydroxyl groups act as donors in O–H⋯Cl and O–H⋯O hydrogen bonds. Both the classical and C–H⋯Cl and C–H⋯O hydrogen bonds differ between the forms, helping to explain the large microstrain observed for the sample. The powder pattern has been submitted to ICDD® for inclusion in the Powder Diffraction File™ (PDF®).

Materials Science↗
DOE OSTI2024

Influence of Alkali Metal Cations on the Oxygen Reduction Activity of Pt 5 Y and Pt 5 Gd Alloys

Electrolyte species can significantly influence the electrocatalytic performance. In this work, we investigate the impact of alkali metal cations on the oxygen reduction reaction (ORR) on active Pt 5 Gd and Pt 5 Y polycrystalline electrodes. Due to the strain effects, Pt alloys exhibit a higher kinetic current density of ORR than pure Pt electrodes in acidic media. In alkaline solutions, the kinetic current density of ORR for Pt alloys decreases linearly with the decreasing hydration energy in the order of Li + > Na + > K + > Rb + > Cs + , whereas Pt shows the opposite trend. To gain further insights into these experimental results, we conduct complementary density functional theory calculations considering the effects of both electrode surface strain and electrolyte chemistry. The computational results reveal that the different trends in the ORR activity in alkaline media can be explained by the change in the adsorption energy of reaction intermediates with applied surface strain in the presence of alkali metal cations. Our findings provide important insights into the effects of the electrolyte and the strain conditions on the electrocatalytic performance and thus offer valuable guidelines for optimizing Pt-based electrocatalysts.

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

Comparing measurement of limiting current in block copolymer electrolytes as a function of salt concentration with theoretical predictions

Optimizing electrolyte performance is crucial for the widespread adoption of electrochemical energy storage. We demonstrate that limiting current provides a robust criterion for determining the optimum electrolyte. Experiments were conducted on rigid block copolymer electrolytes comprising mixtures of polystyrene-block-poly(ethylene oxide) copolymer (SEO) and lithium bis(trifluoromethanesulfonyl) imide salt (LiTFSI) over a salt concentration range from r av = 0.04 to r av = 0.20 (r av is the molar ratio of lithium ions to ethylene oxide). Here we show that the maximum limiting current density is 4.3 mA cm –2 at r av = 0.12. The dependence of limiting current on salt concentration is in good agreement with predictions from Newman's concentrated solution theory with no adjustable parameters.

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

Exciton Transfer Between Extended Electronic States in Conjugated Inter-Polyelectrolyte Complexes

Artificial light harvesting, a process that involves converting sunlight into chemical potential energy, is considered to be a promising part of the overall solution to address urgent global energy challenges. Conjugated polyelectrolyte complexes (CPECs) are particularly attractive for this purpose due to their extended electronic states, tunable assembly thermodynamics, and sensitivity to their local environment. Importantly, ionically assembled complexes of conjugated polyelectrolytes can act as efficient donor-acceptor pairs for electronic energy transfer (EET). However, to be of use in material applications, we must understand how modifying the chemical structure of the CPE backbone alters the EET rate beyond spectral overlap considerations. In this report we investigate the dependence of the EET efficiency and rate on the electronic structure and excitonic wave function of the CPE backbone. To do so, we synthesized a series of alternating copolymers where the electronic states are systematically altered by introducing comonomers with electron withdrawing and electron-rich character while keeping the linear ionic charge density nearly fixed. We find evidence that the excitonic coupling may be significantly affected by the exciton delocalization radius, in accordance with analytical models based on the line-dipole approximation and quantum chemistry calculations. Our results imply that care should be taken when selecting CPE components for optimal CPEC EET. These results have implications for using CPECs as key components in water-based light-harvesting materials, either as standalone assemblies or as adsorbates on nanoparticles and thin films.

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

Pulsed Light Synthesis of High Entropy Nanocatalysts with Enhanced Catalytic Activity and Prolonged Stability for Oxygen Evolution Reaction

The ability to synthesize compositionally complex nanostructures rapidly is a key to high-throughput functional materials discovery. In addition to being time-consuming, a majority of conventional materials synthesis processes closely follow thermodynamics equilibria, which limit the discovery of new classes of metastable phases such as high entropy oxides (HEO). Herein, a photonic flash synthesis of HEO nanoparticles at timescales of milliseconds is demonstrated. By leveraging the abrupt heating and cooling cycles induced by a high-power-density xenon pulsed light, mixed transition metal salt precursors undergo rapid chemical transformations. Hence, nanoparticles form within milliseconds with a strong affinity to bind to the carbon substrate. Oxygen evolution reaction (OER) activity measurements of the synthesized nanoparticles demonstrate two orders of magnitude prolonged stability at high current densities, without noticeable decay in performance, compared to commercial IrO 2 catalyst. This superior catalytic activity originates from the synergistic effect of different alloying elements mixed at a high entropic state. It is found that Cr addition influences surface activity the most by promoting higher oxidation states, favoring optimal interaction with OER intermediates. The proposed high-throughput method opens new pathways toward developing next-generation functional materials for various electronics, sensing, and environmental applications, in addition to renewable energy conversion.

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

Mesoscale Modeling Approach for Quantifying Microstructure-Aware Micromechanical Responses in Metal Hydrides

Metal hydrides can undergo significant volume changes upon hydrogen uptake and release, which induce a mechanical response that depends not only on the evolving hydrogen composition but also on the microstructure. We present a comprehensive mesoscale modeling framework based on microelasticity theory to quantify the micromechanical responses of metal hydrides, specifically focusing on a hydrogenating polycrystalline MgH 2x particle within a host material as a model micromechanical system. Utilizing digitally generated realistic microstructures and density-functional-theory-derived parameters, we analyzed highly nonuniform local stress profiles in the polycrystalline hydrides under the clamping force exerted by the host during hydrogenation. Our framework also allows us to predict the corresponding strain energy accumulation and mechanical hot spots formation in the hydrides, highlighting their roles in thermodynamic destabilization and mechanical failure, respectively. Through extensive parametric simulations, we further quantified the influence of interface type, crystallinity, grain size, loading ratio, and host stiffness, providing practical guidance for optimizing microstructural design and host material selection. This proposed approach is broadly applicable to micromechanical systems with complex microstructural features involving chemical reaction- and/or phase-transformation-induced deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Complexions at the iron-magnetite interface

Synthesizing distinct phases and controlling crystalline defects are key concepts in materials design. These approaches are often decoupled, with the former grounded in equilibrium thermodynamics and the latter in nonequilibrium kinetics. By unifying them through defect phase diagrams, we can apply phase equilibrium models to thermodynamically evaluate defects—including dislocations, grain boundaries, and phase boundaries—establishing a theoretical framework linking material imperfections to properties. Using scanning transmission electron microscopy (STEM) with differential phase contrast (DPC) imaging, we achieve the simultaneous imaging of heavy Fe and light O atoms, precisely mapping the atomic structure and chemical composition at the iron-magnetite (Fe/Fe 3 O 4 ) interface. We identify a well-ordered two-layer interface-stabilized phase state (referred to as complexion) at the Fe[001]/Fe 3 O 4 [001] interface. Using density-functional theory (DFT), we explain the observed complexion and map out various interface-stabilized phases as a function of the O chemical potential. The formation of complexions increases interface adhesion by 20% and alters charge transfer between adjacent materials, impacting transport properties. Our findings highlight the potential of tunable defect-stabilized phase states as a degree of freedom in materials design, enabling optimized corrosion protection, catalysis, and redox-driven phase transitions, with applications in materials sustainability, efficient energy conversion, and green steel production.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Parameterizing empirical interatomic potentials for predicting thermophysical properties via an irreducible derivative approach: the case of ThO 2 and UO 2

The accuracy of classical physical property predictions using molecular dynamics simulations is determined by the quality of the interatomic potentials. Here we introduce a training approach for empirical interatomic potentials (EIPs) which is well suited for capturing phonons and phonon-related properties. Our approach is based on direct comparisons of the second- and third-order irreducible derivatives (IDs) between an EIP and the Born–Oppenheimer potential within density functional theory (DFT) calculations. IDs fully exploit space group symmetry and allow for training without redundant information. We demonstrate the fidelity of our approach in the context of ThO 2 and UO 2 , where we optimize parameters of an embedded-atom method potential in addition to core–shell interactions. Our EIPs provide thermophysical properties in good agreement with DFT and outperform widely utilized EIPs for phonon dispersion and thermal conductivity predictions. Reasonable estimates of thermal expansion and formation energies of Frenkel pairs are also obtained.

empirical interatomic potential↗
DOE OSTI2025

Understanding Redox Organic Behavior in Deep Eutectic Solvents: Considerations for Molecular Design

Electrolytes based on deep eutectic solvents (DESs) coupled with redox active organic molecules have shown potential as a versatile and energy dense electrochemical energy storage system. However, progress in these systems has been held back by a lack of understanding of the irregular behavior displayed when redox active organic molecules are transitioned from other solvent systems. In this work, the hydrogen bonding characteristics of a series of redox organic molecules were investigated through infrared spectroscopy and molecular modeling. New understanding of these interactions was then used to explain their electrochemical behavior in a DES electrolyte. A model was used to predict the behavior of new derivatives towards the design of an optimized redox organic-DES system. Hydrogen bonding between the redox molecules and the solvent was found to significantly shift the potential of a redox reaction more positive when a hydrogen bond forms at the redox active site. It was predicted that functionalizing a molecule with electron withdrawing groups to lower the electron density of the redox active functional group lowers the strength of the hydrogen bond and thus alleviates the undesirable potential shift. This hypothesis was demonstrated by the addition of nitro groups to fluorenones.

25 ENERGY STORAGE↗
DOE OSTI2020

Reaction Mechanism of Vapor-Phase Formic Acid Decomposition over Platinum Catalysts: DFT, Reaction Kinetics Experiments, and Microkinetic Modeling

A combination of periodic density functional theory (DFT, PW91-GGA) calculations, reaction kinetics experiments, and mean-field microkinetic modeling is used to derive insights on the reaction mechanism and determine the nature of the active site under reaction conditions for the vapor-phase decomposition of formic acid (FA, HCOOH) over Pt/C catalysts. Microkinetic models formulated using DFT energetics derived on the clean Pt(100) and Pt(111) required large parameter adjustments to reproduce the experimentally measured apparent activation energies and reaction orders. Further, these models predicted high surface coverage of adsorbed carbon monoxide (CO*), inconsistent with the environment of the active site in the DFT calculations on the clean surfaces. Consequently, we reperformed DFT calculations for the entire reaction network on partially CO*-covered (4/9 monolayer, ML) Pt(111) and Pt(100). The resultant microkinetic models, with thermochemistry and kinetics explicitly dependent on CO* coverage, were able to reproduce the experimentally determined activation energies and reaction orders, in addition to being self-consistent in CO* coverage. Our results suggest that Pt(100) is likely poisoned by CO* under typical reaction conditions and does not contribute significantly to the experimentally observed reactivity. Instead, we find that Pt(111) better represents the active site for FA decomposition reaction on Pt/C catalysts. The optimized model on 4/9 ML CO*-covered Pt(111) suggests that the reaction occurs via the carboxyl (COOH*) intermediate and that the spectator CO*-assisted pathways play a significant role under reaction conditions. Lastly, this study underscores the importance of spectator species on the energetics and the mechanism of a catalytic reaction and their key role in developing a model that better addresses the nature of the active site under realistic catalytic reaction conditions.

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

Facet-Dependent Photodegradation of Methylene Blue by Hematite Nanoplates in Visible Light

The expression of specific crystal facets in different nanostructures is known to play a vital role in determining the sensitivity toward photodegradation of organics, which can generally be ascribed to differences in surface structure and energy. Herein, we report the synthesis of hematite nanoplates with controlled relative exposure of basal (001) and edge (012) facets, enabling us to establish direct correlation between surface structure and the photocatalytic degradation efficiency of methylene blue (MB) in the presence of hydrogen peroxide. MB adsorption experiments showed that the capacity on (001) is about three times larger than on (012), and this is supported by density functional theory calculations suggesting adsorption energies about 1.5 times higher. However, this alone cannot explain the observed MB photodegradation rate, which is around 14.5 times faster on (001) than (012). The enhancement we attribute to a higher availability of photoelectron accepting surface Fe3+ sites on the (001) facet. This facilitates more efficient iron valence cycling and the heterogeneous photo-Fenton reaction yielding MB-oxidizing hydroxyl radicals at the surface. Our findings help establish a rational basis for design and optimization of hematite nanostructures as photocatalysts for environmental remediation.

Zong, Meirong↗
DOE OSTI2022

237 Np Mössbauer Isomer Shifts: A Lesson About the Balance of Static and Dynamic Electron Correlation in Heavy Element Complexes

A large set of neptunium compounds with different oxidation states (III to VII) was assembled to study the Mössbauer isomer shift by wave function calculations and better understand covalency in f-elements complexes. The contact density approach was used to calculate the isomer shift using complete active space self-consistent field (CASSCF) multiconfiguration wave functions, as well as matrix product states [from Density Matrix Renormalization Group (DMRG) algorithms] for large active spaces. Dynamic correlation effects for the isomer shifts were treated via CASPT2 energy derivatives with respect to the nuclear radius. The CASSCF calculations appear to produce different orbital overlocalization errors for low and high Np oxidation states. For compounds with low Np oxidation numbers, the errors can be attributed to the overlocalization of the 5f orbitals. For the compounds with high Np oxidation numbers, the main errors arise from an overlocalization of ligand orbitals and concomitant to weak donation bonding. Attempts to mitigate the overlocalization errors with large active spaces using DMRG were only partially successful, showing that explicit treatment of dynamic correlation is necessary for accurate predictions of Mössbauer isomer shifts. The CASPT2 calculations perform very satisfactorily. For a subset of Np compounds, both static and dynamic correlation effects were substantial. Furthermore, a rational active space selection based on orbital entanglement diagrams proved beneficial for determining the optimal reference wave function.

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

Density Functionals for Hydrogen Storage: Defining the H2Bind275 Test Set with Ab Initio Benchmarks and Assessment of 55 Functionals

Efficient and high-capacity storage materials are indispensable for a hydrogen-based economy. In silico tools can accelerate the process of discovery of new adsorbent materials with optimal hydrogen adsorption enthalpies. Density functional theory is well-poised to become a very useful tool for enabling high-throughput screening of potential materials. In this work, we have identified density functional approximations that provide good performance for hydrogen binding applications following a two-pronged approach. First, we have compiled a data set (H2Bind275) that comprehensively represents the hydrogen binding problem capturing the chemical and mechanistic diversity in the binding sites encountered in hydrogen storage materials. We have also computed reference interaction energies for this data set using coupled-cluster theory. Second, we have assessed the performance of 55 density functional approximations for predicting H 2 interaction energies and have identified two hybrid density functionals (ωB97X-V and ωB97M-V), two double hybrid density functionals (DSD-PBEPBE-D3(BJ) and PBE0-DH), and one semilocal density functional (B97M-V) as the best performing ones. We have recommended the addition of empirical dispersion corrections to systematically underbinding density functionals such as revPBE, BLYP, and B3LYP for improvements in performance at negligible additional cost. We have also recommended the usage of the def2-TZVPP basis set as it represents a good compromise between accuracy and cost, limiting the finite basis set errors to less than 1 kJ/mol.

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

Accelerating defect predictions in semiconductors using graph neural networks

First-principles computations reliably predict the energetics of point defects in semiconductors but are constrained by the expense of using large supercells and advanced levels of theory. Machine learning models trained on computational data, especially ones that sufficiently encode defect coordination environments, can be used to accelerate defect predictions. Here, we develop a framework for the prediction and screening of native defects and functional impurities in a chemical space of group IV, III–V, and II–VI zinc blende semiconductors, powered by crystal Graph-based Neural Networks (GNNs) trained on high-throughput density functional theory (DFT) data. Using an innovative approach of sampling partially optimized defect configurations from DFT calculations, we generate one of the largest computational defect datasets to date, containing many types of vacancies, self-interstitials, anti-site substitutions, impurity interstitials and substitutions, as well as some defect complexes. We applied three types of established GNN techniques, namely crystal graph convolutional neural network, materials graph network, and Atomistic Line Graph Neural Network (ALIGNN), to rigorously train models for predicting defect formation energy (DFE) in multiple charge states and chemical potential conditions. We find that ALIGNN yields the best DFE predictions with root mean square errors around 0.3 eV, which represents a prediction accuracy of 98% given the range of values within the dataset, improving significantly on the state-of-the-art. We further show that GNN-based defective structure optimization can take us close to DFT-optimized geometries at a fraction of the cost of full DFT. The current models are based on the semi-local generalized gradient approximation-Perdew–Burke–Ernzerhof (PBE) functional but are highly promising because of the correlation of computed energetics and defect levels with higher levels of theory and experimental data, the accuracy and necessity of discovering novel metastable and low energy defect structures at the PBE level of theory before advanced methods could be applied, and the ability to train multi-fidelity models in the future with new data from non-local functionals. The DFT-GNN models enable prediction and screening across thousands of hypothetical defects based on both unoptimized and partially optimized defective structures, helping identify electronically active defects in technologically important semiconductors.

Rahman, Md Habibur (ORCID:000000027705984X)↗
DOE OSTI2020

Towards defect engineering in hexagonal MoS 2 nanosheets for tuning hydrogen evolution and nitrogen reduction reactions

Combined computational and experimental approaches were used to evaluate defective 2H-MoS 2 nanosheets for their activity and selectivity for hydrogen evolution reaction (HER) and nitrogen reduction reaction (NRR). Density functional theory calculations were used to understand the relationship between HER and NRR activity on the ideal basal MoS 2 plane, seven grain-boundaries, ten single-/few-atom vacancies and anti-sites, and zigzag and armchair edge sites. The results confirm that 2H-MoS 2 should contain several defects with high activity for HER: armchair and zigzag edges, VS vacancy, MoS 2 anti-site, and S-S and Mo-Mo grain boundaries. Considering Gibbs free energy change for all the steps in the NRR mechanism and kinetic barriers for a key NRR step, we have found that activation of perspective NRR selective sites in 2H-MoS 2 , namely V MoS6 and clusters of S-vacancies, would require large overpotential, conditions at which HER dominates. The DFT conclusions are supported by the electrochemical studies of NRR activity and selectivity under aqueous conditions, which show an increase in NRR activity but a decrease in Faradaic efficiency as applied cell potential becomes more negative. The results of this work therefore highlight the challenges in activating natural 2H-MoS 2 for NRR, which would require additional material engineering or reaction condition optimization as a way to suppress HER, decrease the NRR overpotential or preferentially both.

2H-MoS2↗