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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Including Physics-Informed Atomization Constraints in Neural Networks for Reactive Chemistry

Machine learning interatomic potentials (MLIPs) have emerged as powerful tools for investigating atomistic systems with high accuracy and a relatively low computational cost. However, a common and unaddressed challenge with many current neural network (NN) MLIP models is their limited ability to accurately predict the relative energies of systems containing isolated or nearly isolated atoms, which appear in various reactive processes. To address this limitation, we present a mathematical technique for modifying any existing atom-centered NN architecture to account for the energies of isolated atoms. The result produces a consistent prediction of the atomization energy (AE) of a system using minimal constraints on the model. Using this technique, we build a model architecture that we call hierarchically interacting particle neural network (HIP-NN)-AE, an AE-constrained version of the HIP-NN, as well as ANI-AE, the AE-constrained version of the accurate NN engine for molecular energies (ANI). Our results demonstrate AE consistency of AE-constrained models, which drastically improves the AE predictions for the models. We compare the AE-constrained approach to unconstrained models as well as models from the literature in other scenarios, such as bond dissociation energies, bond dissociation pathways, and extensibility tests. These results show that the constraints improve the model performance in some of these tasks and do not negatively affect the performance on any tasks. The AE constraint approach thus offers a robust solution to the challenges posed by isolated atoms in energy prediction tasks.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Rational optimization of substituted α–MnO 2 cathode for aqueous zinc–ion battery

Density functional theory (DFT) is utilized to explore the effects of increasing concentrations of vanadium (V) and chromium (Cr) substitution on the discharge of α-MnO 2 cathode in hydrated zinc-ion batteries (ZIBs). During H + -intercalation and Zn 2+ -intercalation, Cr-substitution proves to be more effective than V-substitution in promoting discharge behaviors. Transitioning from Mn 0.875 Cr 0.125 O 2 , Mn 0.75 Cr 0.25 O 2 to Mn 0.625 Cr 0.375 O 2 is found to consistently enhance the discharge voltage, along with improved tunnel structure retention and volume expansion suppression. In comparison, the promoting effect of increasing V-substitution is relatively small at initial discharge stages and leads to degradation at later stages, primarily due to an increased concentration of unstable Mn 2+ ion. The superior effect of Cr-substitution is attributed to the unique atomic and electronic structures of substituted Cr 4+ and reduced Cr 3+ ions during discharge. These ions serve as active electron acceptors to limit the formation of Mn 3+ and Mn 2+ ions, and as anchors to stabilize the α-MnO 2 framework and intercalated H + /Zn 2+ ions, respectively. Finally, our study highlights fine-tuning through substitution to enhance the performance of α-MnO 2 -based cathode materials in ZIBs.

25 ENERGY STORAGE↗
DOE OSTI2021

Defect Chemistry and Hydrogen Transport in La/Sr-Based Oxyhydrides

Oxyhydrides in the series La 2 – y Sr y LiH 1 + y O 3 – y are attractive materials for solid-state hydride electrolytes. However, their chemical stability is poor, and thus far only La 2 LiHO 3 (y = 0) has been demonstrated in a working fuel cell. Using a first-principles approach to study defect chemistry and stability, we find that La 2 LiHO 3 can be stabilized over a range of H-rich, O-poor chemical conditions, while Sr 2 LiH 3 O has low stability. Here, defects are critical for ionic transport and for stability, particularly in Sr 2 LiH 3 O, in which anion antisites and H vacancies readily form. In fact, we show that the ground-state structure of Sr 2 LiH 3 O is inherently disordered. La 2 LiHO 3 has low conductivity, which we connect to the high formation energy of defects such as VO and OH that are necessary for hydride conduction. When Sr 2 LiH 3 O is grown under carefully selected conditions, V H + defects can be made prevalent; these defects give rise to much higher ionic conductivity than can be obtained in La 2 LiHO 3 . In general, we show that the choice of synthesis conditions is vital when seeking to optimize the stability and ionic conductivity of these oxyhydrides, thereby tailoring them for use in solid-state hydrogen fuel cells.

08 HYDROGEN↗
DOE OSTI2022

Multi-objective goal-directed optimization of de novo stable organic radicals for aqueous redox flow batteries

Abstract Advances in the field of goal-directed molecular optimization offer the promise of finding feasible candidates for even the most challenging molecular design applications. One example of a fundamental design challenge is the search for novel stable radical scaffolds for an aqueous redox flow battery that simultaneously satisfy redox requirements at the anode and cathode, as relatively few stable organic radicals are known to exist. To meet this challenge, we develop a new open-source molecular optimization framework based on AlphaZero coupled with a fast, machine-learning-derived surrogate objective trained with nearly 100,000 quantum chemistry simulations. The objective function comprises two graph neural networks: one that predicts adiabatic oxidation and reduction potentials and a second that predicts electron density and local three-dimensional environment, previously shown to be correlated with radical persistence and stability. With no hard-coded knowledge of organic chemistry, the reinforcement learning agent finds molecule candidates that satisfy a precise combination of redox, stability and synthesizability requirements defined at the quantum chemistry level, many of which have reasonable predicted retrosynthetic pathways. The optimized molecules show that alternative stable radical scaffolds may offer a unique profile of stability and redox potentials to enable low-cost symmetric aqueous redox flow batteries.

25 ENERGY STORAGE↗
DOE OSTI2020

Transition Metal Nitrides as Promising Catalyst Supports for Tuning CO/H 2 Syngas Production from Electrochemical CO 2 Reduction

The electrochemical carbon dioxide reduction reaction (CO 2 RR) to produce synthesis gas (syngas) with tunable CO/H 2 ratios has been studied by supporting Pd catalysts on transition metal nitride (TMN) substrates. Combining experimental measurements and density functional theory (DFT) calculations, Pd-modified niobium nitride (Pd/NbN) is found to generate much higher CO and H 2 partial current densities and greater CO Faradaic efficiency than Pd-modified vanadium nitride (Pd/VN) and commercial Pd/C catalysts. In-situ X-ray diffraction identifies the formation of PdH in Pd/NbN and Pd/C under CO 2 RR conditions, whereas the Pd in Pd/VN is not fully transformed into the active PdH phase. DFT calculations show that the stabilized *HOCO and weakened *CO intermediates on PdH/NbN are critical to achieving higher CO 2 RR activity. This work suggests that NbN is a promising substrate to modify Pd, resulting in an enhanced electrochemical conversion of CO 2 to syngas with a potential reduction in precious metal loading.

25 ENERGY STORAGE↗
DOE OSTI2022

Revisiting the Bonding Model for Gold(I) Species: The Importance of Pauli Repulsion Revealed in a Gold(I)–Cyclobutadiene Complex

Understanding the bonding of gold(I) species has been central to the development of gold(I) catalysis. Herein, we present the synthesis and characterization of the first gold(I)-cyclobutadiene complex, accompanied with bonding analysis by state-of-the-art energy decomposition analysis methods. Analysis of possible coordination modes for the new species not only confirms established characteristics of gold(I) bonding, but also suggests that Pauli repulsion is a key yet hitherto overlooked element. Additionally, we obtain a new perspective on gold(I)-bonding by comparison of the gold(I)-cyclobutadiene to congeners stabilized by p-, d-, and f-block metals. Consequently, we refine the gold(I) bonding model, with a delicate interplay of Pauli repulsion and charge transfer as the key driving force for various coordination motifs. Finally, Pauli repulsion is similarly determined as a significant interaction in Au I -alkyne species, corroborating this revised understanding of Au I bonding.

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

Mitigation of the internal p-n junction in CoS 2 -contacted FeS 2 single crystals: Accessing bulk semiconducting transport

Pyrite FeS 2 is an outstanding candidate for a low-cost, nontoxic, sustainable photovoltaic material, but efficient pyrite-based solar cells are yet to materialize. Recent studies of single crystals have shed much light on this by uncovering a p-type surface inversion layer on n-type (S-vacancy doped) crystals, and the resulting internal p-n junction. This leaky internal junction likely plays a key role in limiting efficiency in pyrite-based photovoltaic devices, also obscuring the true bulk semiconducting transport properties of pyrite crystals. Here, we demonstrate complete mitigation of the internal p-n junction in FeS 2 crystals by fabricating metallic CoS 2 contacts via a process that simultaneously diffuses Co (a shallow donor) into the crystal, the resulting heavy n doping yielding direct Ohmic contact to the interior. Low-temperature bulk transport studies of controllably Co- and S-vacancy doped semiconducting crystals then enable a host of previously inaccessible observations and measurements, including determination of donor activation energies (which are as low as 5 meV for Co), observation of an unexpected second activated transport regime, realization of electron mobility up to 2100 cm 2 V –1 s –1 , elucidation of very different mobilities in Co- and S-vacancy-doped cases, and observation of an abrupt temperaturedependent crossover to bulk Efros-Shklovskii variable-range hopping, accompanied by an unusual form of nonlinear Hall effect. Aspects of the results are interpreted with the aid of first-principles electronic structure calculations on both Co- and S-vacancy-doped FeS 2 . Furthermore, this work thus demonstrates unequivocal mitigation of the internal p-n junction in pyrite single crystals, with important implications for both future fundamental studies and photovoltaic devices

36 MATERIALS SCIENCE↗
DOE OSTI2022

Computing x-ray absorption spectra from linear-response particles atop optimized holes

State specific orbital optimized density functional theory (OO-DFT) methods, such as restricted open-shell Kohn–Sham (ROKS), can attain semiquantitative accuracy for predicting x-ray absorption spectra of closed-shell molecules. OO-DFT methods, however, require that each state be individually optimized. In this Communication, we present an approach to generate an approximate core-excited state density for use with the ROKS energy ansatz, which is capable of giving reasonable accuracy without requiring state-specific optimization. Herein, this is achieved by fully optimizing the core-hole through the core-ionized state, followed by the use of electron-addition configuration interaction singles to obtain the particle level. This hybrid approach can be viewed as a DFT generalization of the static-exchange (STEX) method and can attain ~0.6 eV rms error for the K-edges of C–F through the use of local functionals, such as PBE and OLYP. This ROKS(STEX) approach can also be used to identify important transitions for full OO ROKS treatment and can thus help reduce the computational cost of obtaining OO-DFT quality spectra. ROKS(STEX), therefore, appears to be a useful technique for the efficient prediction of x-ray absorption spectra.

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

Deep potential molecular dynamics simulations of low-temperature plasma-surface interactions

Machine learning approaches to potential generation for molecular dynamics (MD) simulations of low-temperature plasma-surface interactions could greatly extend the range of chemical systems that can be modeled. Empirical potentials are difficult to generalize to complex combinations of multiple elements with interactions that might include covalent, ionic, and metallic bonds. This work demonstrates that a specific machine learning approach, Deep Potential Molecular Dynamics (DeepMD), can generate potentials that provide a good model of plasma etching in the Si-Cl-Ar system. Comparisons are made between MD results using DeepMD models and empirical potentials, as well as experimental measurements. Pure Si properties predicted by the DeepMD model are in reasonable agreement with experimental results. Simulations of Si bombardment by Ar + ions demonstrate the ability of the DeepMD method to predict sputtering yields as well as the depth of the amorphous-crystalline interface. Etch yields as a function of flux ratio and ion energy for simultaneous Cl 2 and Ar + impacts are in good agreement with previous simulation results and experiment. Predictions of etch yields and etch products during plasma-assisted atomic layer etching of Si-Cl 2 -Ar are shown to be in good agreement with MD predictions using empirical potentials and with experiment. Finally, good agreement was also seen with measurements for the spontaneous etching of Si by Cl atoms at 300 K. Further, the demonstration that DeepMD can reproduce results from MD simulations using empirical potentials is a necessary condition to future efforts to extend the method to a much wider range of systems for which empirical potentials may be difficult or impossible to obtain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Photoionization of C 60 at high energies

Calculations of the photoionization cross section of C60 including all valence and core subshells have been carried out for photon energies up to 1 keV using density functional theory (DFT) and time-dependent DFT (TDDFT) within the framework of a fully molecular model. The high-energy valence subshell cross sections behave rather differently than the results from model potentials where the cross section falls far too rapidly with energy. This unphysical behavior of the model calculations is traced to smearing out of the carbon nuclei in the model potentials which makes it difficult for momentum to be conserved at the higher energies. Comparison with 60 times carbon atomic cross section sheds light on similarities and differences. The high energy behavior of individual valence cross sections reflects the amount of C 2s/2p mixing in the molecular orbitals.

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

Metal Independent Correlations for Site-specific Binding Energies of Relevant Catalytic Intermediates

Establishing energy correlations among different metals can accelerate the discovery of efficient and cost-effective catalysts for complex reactions. Using a recently introduced coordination-based model, we can predict site specific metal binding energies (∆EM) that can be used as a descriptor for chemical reactions. In this study, we have examined a range of metals including Ag, Au, Co, Cu, Ir, Ni, Os, Pd, Pt, Rh, and Ru and found linear correlations between predicted ∆EM and adsorption energies of CH and O (∆ECH and ∆EO) at various coordination environment for all the considered metals. Interestingly, all the metals correlate with one another under specific surface site coordination, indicating that different metals are interrelated in a particular coordination environment. Furthermore, we have tested and verified for PtPd and PtIr-based alloys that they follow a similar behavior. Moreover, we have expanded the metal space by taking some early transition metals along with few s-block metals and shown a cyclic behavior of the adsorbate binding energy (∆EA) versus ∆EM. Therefore, ∆ECH and ∆EO can be efficiently interpolated between metals, alloys and intermetallics based on information related to one metal only. This simplifies the process of screening new metal catalyst formulations and their reaction energies.

Density function theory↗
NASA NTRS1993

Two-integral distribution functions for axisymmetric galaxies

The new method presented for finding distribution functions, which depend only on the classical integrals of energy and angular momentum for stellar systems with known axisymmetric densities, is the analog for the axisymmetric case of Eddington's classical solution for the isotropic distribution function, depending only on energy, of a known spherical density. It is required that density be expressed as a function of the potential and of a radial coordinate. Our solution is also an integral which is derived directly from the density, and hence can be used with complicated densities. A numerical quadrature is generally required to evaluate this solution, but contour integrals can be computed accurately by numerical quadrature; this is preferable to an explicit evaluation if the latter is an infinite series, such as is obtained using Fricke's method. We give several examples, including some for which our distribution functions are new. Our method can be extended simply to the related problems of finding anisotropic distribution functions for spherical or disk systems.

Hunter, C.↗
DOE OSTI2024

A comparison of QTP functionals against coupled-cluster methods for EAs of small organic molecules

EA-EOM-CCSD electron affinities and LUMO energies of various Kohn–Sham density functional theory (DFT) methods are calculated for an a priori IP benchmark set of 64 small, closed-shell molecules. The purpose of these calculations was to investigate whether the QTP KS-DFT functionals can emulate EA-EOM-CC with only a mean-field approximation. We show that the accuracy of DFT—relative to CCSD—improves significantly when elements of correlated orbital theory are introduced into the parameterization to define the QTP family of functionals. In particular, QTP(02), which has only a single range separation parameter, provides results accurate to a MAD of <0.15 eV for the whole set of 64 molecules compared to EA-EOM-CCSD, far exceeding the results from the non-QTP family of density functionals.

Chemistry↗
DOE OSTI2023

A High‐Rate Li–CO 2 Battery Enabled by 2D Medium‐Entropy Catalyst

Abstract Lithium‐air batteries based on CO 2 reactant (Li–CO 2 ) have recently been of interest because it has been found that reversible Li/CO 2 electrochemistry is feasible. In this study, a new medium‐entropy cathode catalyst, (NbTa) 0.5 BiS 3 , that enables the reversible electrochemistry to operate at high rates is presented. This medium entropy cathode catalyst is combined with an ionic liquid‐based electrolyte blend to give a Li–CO 2 battery that operates at high current density of 5000 mA g −1 and capacity of 5000 mAh g −1 for up to 125 cycles, far exceeding reported values in the literature for this type of battery. The higher rate performance is believed to be due to the greater stability of the multi‐element (NbTa) 0.5 BiS 3 catalyst because of its higher entropy compared to previously used catalysts with a smaller number of elements with lower entropies. Evidence for this comes from computational studies giving very low surface energies (high surface stability) for (NbTa) 0.5 BiS 3 and transmission electron microscopystudies showing the structure being retained after cycling. In addition, the calculations indicate that Nb‐terminated surface promotes Li–CO 2 electrochemistry resulting in Li 2 CO 3 and carbon formation, consistent with the products found in the cell. These results open new direction to design and develop high‐performance Li–CO 2 batteries.

25 ENERGY STORAGE↗
DOE OSTI2021

A General Strategy to Immobilize Single‐Atom Catalysts in Metal–Organic Frameworks for Enhanced Photocatalysis

Abstract Single‐atom catalysts (SACs) are witnessing rapid development due to their high activity and selectivity toward diverse reactions. However, it remains a grand challenge in the general synthesis of SACs, particularly featuring an identical chemical microenvironment and on the same support. Herein, a universal synthetic protocol is developed to immobilize SACs in metal–organic frameworks (MOFs). Significantly, by means of SnO 2 as a mediator or adaptor, not only different single‐atom metal sites, such as Pt, Cu, and Ni, etc., can be installed, but also the MOF supports can be changed (for example, UiO‐66‐NH 2 , PCN‐222, and DUT‐67) to afford M 1 /SnO 2 /MOF architecture. Taking UiO‐66‐NH 2 as a representative, the Pt 1 /SnO 2 /MOF exhibits approximately five times higher activity toward photocatalytic H 2 production than the corresponding Pt nanoparticles (≈2.5 nm) stabilized by SnO 2 /UiO‐66‐NH 2 . Remarkably, despite featuring identical parameters in the chemical microenvironment and support in M 1 /SnO 2 /UiO‐66‐NH 2 , the Pt 1 catalyst possesses a hydrogen evolution rate of 2167 µmol g –1 h –1 , superior to the Cu 1 and Ni 1 counterparts, which is attributed to the differentiated hydrogen binding free energies, as supported by density‐functional theory (DFT) calculations. This is thought to be the first report on a universal approach toward the stabilization of SACs with identical chemical microenvironment on an identical support.

Sui, Jianfei↗
DOE OSTI2020

A First–Principles–Based Sub–Lattice Formalism for Predicting Off–Stoichiometry in Materials for Solar Thermochemical Applications: The Example of Ceria

Theoretical models that reliably can predict off-stoichiometry in materials via accurate descriptions of underlying thermodynamics are crucial for energy applications. For example, transition-metal and rare-earth oxides that can tolerate a large number of oxygen vacancies, such as CeO 2 and doped CeO 2 , can split water and carbon dioxide via a two-step, oxide-based solar thermochemical (STC) cycle. The search for new STC materials with a performance superior to that of state-of-the-art CeO 2 can benefit from predictions accurately describing the thermodynamics of oxygen vacancies. The sub-lattice formalism, a common tool used to fit experimental data and build temperature-composition phase diagrams, can be useful in this context. Here, sub-lattice models are derived solely from zero-temperature quantum mechanics calculations to estimate fairly accurate temperature- and oxygen-partial-pressure-dependent off-stoichiometries in CeO 2 and Zr-doped CeO 2 . Physical motivations for deriving some of the “excess” sub-lattice model parameters directly from quantum mechanical calculations, instead of fitting to minimize deviations from experimental and/or theoretical data, are identified. As a result, important limitations and approximations of the approach used are specified and extensions to multi-cation oxides are also suggested to help identify novel candidates for water and carbon dioxide splitting and related applications.

08 HYDROGEN↗
DOE OSTI2022

In‐Silico Device Performance Prediction of Cosensitizer Dye Pairs for Dye‐Sensitized Solar Cells

Abstract Endeavors in the field of dye‐sensitized solar cells (DSCs) have shown great promise when adopting a data‐driven approach to materials discovery, such as successful molecular‐scale predictions of light‐harvesting chromophores. However, predictions of DSC dyes would become much more sophisticated if a molecular‐to‐macroscopic DSC device prediction methodology existed. Thereby, a fully computational pipeline is presented that predicts device‐performance parameters of DSCs which contain varying dye combinations. Optimal pairing of complementary dyes is identified via a data‐driven workflow that affords cosensitized DSCs with maximum power‐conversion efficiencies. Six high‐performing DSC dyes are paired with partner dyes that are screened from a database of 8488 compounds using sequential heuristic filters. Existing models that predict short‐circuit‐current density ( J SC ) and open‐circuit voltage ( V OC ) parameters are adapted to predict singly sensitized and cosensitized DSC performance. The predictions for J sc values of singly sensitized devices match experimental literature values with comparable accuracy to more computationally costly methods. Five out of six dye pairings are predicted to have greater J SC values when cosensitized compared to their corresponding singly sensitized devices, including two pairs that show strong J sc boosts of +13% and +12% when cosensitized. Thus, the prospect of an entirely in‐silico prediction pipeline for DSC performance that can be used to realize the fully automated design of optimized cosensitized DSCs is demonstrated.

14 SOLAR ENERGY↗