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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Inelastic Neutron Scattering from Thin Film Biaxially Oriented Polyethylene Terephthalate

Recent interest in emerging processes for polymer manufacturing and bio-based chemistries for direct chemical recycling/upcycling has motivated new research focused on a deeper understanding of atomic-scale polymer properties and how they influence macroscopic phenomena. Uncovering the fundamental properties of polymers that give rise to macroscopic behavior could enable new pathways for improved recyclability or utilization of alternative “greener” polymer analogues. In this study, the neutron vibrational spectrum was measured for a film of biaxially oriented polyethylene terephthalate (BoPET) using inelastic neutron scattering (INS), to investigate the relationship between the structure and dynamics of a widely used polymer. Compared to conventional spectroscopic techniques, the use of INS is advantageous for polymeric materials due to the absence of selection rules (i.e., all transitions are allowed), broad-band energy range, and considerable sensitivity to hydrogen modes. In order to distinguish the vibrational modes caused by trans and gauche rotational isomerism, the normal modes of vibration were calculated from a density functional theory-optimized structure of crystalline PET (cPET), representative of the all-trans state, and compared with INS from “highly crystalline” PET powder. Although in- and out-of-plane wagging of hydrogens on the ring structure exhibit significant contribution to both BoPET and cPET spectra, the wagging, rocking, and twisting modes of hydrogen on the ethylene glycol group are, in most cases, conformation-specific. These results were further rationalized by investigating the role of hyperconjugation in stabilizing both conformations using the natural bond order method. Here, through comparison of experimental and calculated INS results, this work provides the fundamental basis for discovering the role of structure and dynamics in shaping the macroscopic properties of PET and polymer analogues.

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

Data-driven approach to parameterize SCAN + U for an accurate description of 3 d transition metal oxide thermochemistry

Semilocal density-functional theory (DFT) methods exhibit significant errors for the phase diagrams of transition-metal oxides that are caused by an incorrect description of molecular oxygen and the large self-interaction error in materials with strongly localized electronic orbitals. Empirical and semiempirical corrections based on the DFT + U method can reduce these errors, but the parameterization and validation of the correction terms remains an on-going challenge. Here we develop a systematic methodology to determine the parameters and to statistically assess the results by considering interlinked thermochemical data across a set of transition metal compounds. We consider three interconnected levels of correction terms: (1) a constant oxygen binding correction, (2) Hubbard-U correction, and (3) DFT/DFT + U compatibility correction. The parameterization is expressed as a unified optimization problem. We demonstrate this approach for 3d transition metal oxides, considering a target set of binary and ternary oxides. With a total of 37 measured formation enthalpies taken from the literature, the dataset is augmented by the reaction energies of 1710 unique reactions that were derived from the formation energies by systematic enumeration. To ensure a balanced dataset across the available data, the reactions were grouped by their similarity using clustering and suitably weighted. The parameterization is validated using leave-one-out cross validation (CV), a standard technique for the validation of statistical models. We apply the methodology to the strongly constrained and appropriately normed (SCAN) density functional. Based on the CV score, the error of binary (ternary) oxide formation energies is reduced by 40% (75%) to 0.10 (0.03) eV/atom. A simplified correction scheme that does not involve SCAN/SCAN + U compatibility terms still achieves an error reduction of 30% (25%). The method and tools demonstrated here can be applied to other classes of materials or to parameterize the corrections to optimize DFT + U performance for other target physical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Towards atom-level understanding of metal oxide catalysts for the oxygen evolution reaction with machine learning

Abstract Green hydrogen production is crucial for a sustainable future, but current catalysts for the oxygen evolution reaction (OER) suffer from slow kinetics, despite many efforts to produce optimal designs, particularly through the calculation of descriptors for activity. In this study, we develop a dataset of density functional theory calculations of bulk and surface perovskite oxides, and adsorption energies of OER intermediates, which includes compositions up to quaternary and facets up to (555). We demonstrate that per-site properties of perovskite oxides such as Bader charge or band center can be tuned through element substitution and faceting, and develop a machine learning model that accurately predicts these properties directly from the local chemical environment. We leverage these per-site properties to identify promising perovskites with high theoretical OER activity. The identified design principles and promising materials provide a roadmap for closing the gap between current artificial catalysts and biological enzymes such as photosystem II.

Chemistry↗
DOE OSTI2020

Quantum Chemistry-Informed Active Learning to Accelerate the Design and Discovery of Sustainable Energy Storage Materials

Here we employed Density Functional Theory (DFT) to compute oxidation potentials of 1,400 homobenzylic ether molecules to search for the ideal sustainable redoxmer design. The generated data were used to construct an active learning model based on Bayesian optimization (BO) that targets candidates with desired oxidation potentials utilizing only a minimal number of DFT calculations. The active learning model demonstrated not only significant efficiency improvement over the random selection approach but also robust capability in identifying desired candidates in an untested set of 112,000 homobenzylic ether molecules. Our findings highlight the efficacy of quantum chemistry-informed active learning to accelerate the discovery of materials with desired properties from a vast chemical space.

25 ENERGY STORAGE↗
DOE OSTI2021

Dark Energy Survey Year 3 results: Optimizing the lens sample in a combined galaxy clustering and galaxy-galaxy lensing analysis

We investigate potential gains in cosmological constraints from the combination of galaxy clustering and galaxy-galaxy lensing by optimizing the lens galaxy sample selection using information from Dark Energy Survey (DES) Year 3 data and assuming the DES Year 1 metacalibration sample for the sources. We explore easily reproducible selections based on magnitude cuts in i-band as a function of (photometric) redshift, zphot, and benchmark the potential gains against those using the well-established redMaGiC [E. Rozo et al., Mon. Not. R. Astron. Soc. 461, 1431 (2016)MNRAA40035-871110.1093/mnras/stw1281] sample. We focus on the balance between density and photometric redshift accuracy, while marginalizing over a realistic set of cosmological and systematic parameters. Our optimal selection, the MagLim sample, satisfies i<4zphot+18 and has ∼30% wider redshift distributions but ∼3.5 times more galaxies than redMaGiC. Assuming a wCDM model (i.e. with a free parameter for the dark energy equation of state) and equivalent scale cuts to mitigate nonlinear effects, this leads to 40% increase in the figure of merit for the pair combinations of Ωm, w, and σ8, and gains of 16% in σ8, 10% in Ωm, and 12% in w. Similarly, in ΛCDM, we find an improvement of 19% and 27% on σ8 and Ωm, respectively. We also explore flux-limited samples with a flat magnitude cut finding that the optimal selection, i<22.2, has ∼7 times more galaxies and ∼20% wider redshift distributions compared to MagLim, but slightly worse constraints. We show that our results are robust with respect to the assumed galaxy bias and photometric redshift uncertainties with only moderate further gains from increased number of tomographic bins or the inclusion of bin cross-correlations, except in the case of the flux-limited sample, for which these gains are more significant.

79 ASTRONOMY AND ASTROPHYSICS↗
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

Net Electronic Charge as an Effective Electronic Descriptor for Oxygen Release and Transport Properties of SrFeO 3 -based Oxygen Sorbents

Perovskite oxides, as oxygen sorbents, exhibit excellent potential in thermochemical redox applications such as chemical looping air separation (CLAS), resulting from their excellent redox properties and high tunability. The structural and compositional flexibility of perovskites, while making these mixed oxides highly versatile, also poses challenges for their design and optimization. Moreover, the elevated operating temperature and dynamic oxygen partial pressure swings during redox processes add additional complications for oxygen sorbent design. From a fundamental standpoint, finding simple yet effective descriptors for the effect of cation substitutions on oxygen sorbents’ redox and oxygen transport properties are of significant importance for rational optimization of perovskite-based oxygen sorbents. In the present study, a series of A/B-site-doped SrFeO 3-δ perovskites are investigated using density functional theory (DFT) + U simulations. From these simulations, an effective electronic structural descriptor, the net electronic charge (Δe) of oxygen anions, is extracted based on its excellent correlations with vacancy formation energy, experimentally derived oxygen capacity, and oxygen vacancy migration barrier. Overall, these findings provide an effective tool to correlate the oxygen sorbents’ redox performances with a simple yet physically meaningful descriptor for rational design and optimization of perovskites in the context of CLAS.

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

Exploring the Effects of Node Topology, Connectivity, and Metal Identity on the Binding of Nerve Agents and Their Hydrolysis Products in Metal–Organic Frameworks

Recent studies have shown that metal–organic frameworks (MOFs) built from hexanuclear M(IV) oxide cluster nodes are effective catalysts for nerve agent hydrolysis, where the properties of the active sites on the nodes can strongly influence the reaction energetics. The connectivity and metal identity of these M6 nodes can be easily tuned, offering extensive opportunities for computational screening to predict promising new materials. Thus, we used density functional theory (DFT) to examine the effects of node topology, connectivity, and metal identity on the binding energies of multiple nerve agents and their corresponding hydrolysis products. By computing an optimization metric based on the relative binding strengths of key hydrolysis reaction species (water, agent, and bidentate-bound products), we predicted optimal M6 nodes for hydrolyzing specific nerve agent and simulant molecules, where our results are in qualitative agreement with observed experimental trends. This analysis highlighted the notion that no single metal or node topology is optimal for all possible organophosphates, suggesting that MOFs should be selected based on the agent of interest. Using the large amount of data generated from our DFT calculations, we then derived quantitative structure–activity relationship (QSAR) models to help explain the complex trends observed in the binding energies. Through linear regression, we identified the most important descriptors for describing the binding of nerve agents and their hydrolysis products to M6 nodes. These results suggested that both molecular and node properties, including both structural and chemical features, collectively contribute to the binding energetics. By performing a thorough statistical analysis, we showed that our QSAR models are capable of making quantitatively accurate binding energy predictions for nerve agents and their hydrolysis products in a wide variety of M(IV)-MOFs. In conclusion, the insights gained herein can be used to guide future experiments for the synthesis of MOFs with enhanced catalytic activity for organophosphate hydrolysis.

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

2D in-Plane Ordered MXene Nanosheets Derived from (Mo 2/3 Er 1/3 ) 2 AlC Rare-Earth i-MAX for Energy Storage Applications

MXenes have become one of the most versatile families of two-dimensional (2D) materials due to their high conductivity, hydrophilicity, and remarkable electrochemical performance. This has stimulated intense efforts to design and synthesize MXenes, including structurally unique in-plane ordered 2D MXenes called i-MXenes. Here, we have synthesized the quaternary rare earth (RE)-based i-MAX phase (Mo 2/3 Er 1/3 ) 2 AlC using an arc melting method, and the corresponding 2D i-MXene was then obtained through a LiF/HCl soft etching process. Literature studies have shown that Al and the RE element are etched out during the etching process, leading to the formation of pure vacancy-ordered Mo1.33C 2D i-MXene. However, our investigation reveals that upon exposure to a fluorine solution, the i-MAX phase forms RE fluoride impurities, which are challenging to remove through HCl−DI water washing and persist in the final product, resulting in impure Mo 1.33 C@Er i-MXene. These results were confirmed by various characterizations such as X-ray diffraction, Raman spectroscopy, X-ray photoelectron spectroscopy, and scanning transmission electron microscopy. Although the Mo 1.33 C@Er electrode showed a 24-fold increase in specific capacitance compared to its parent i-MAX phase, it still exhibited a high charge-transfer resistance arising from the insulating nature of RE fluoride byproducts, which adversely influence the overall capacitance behavior of the synthesized 2D Mo 1.33 C@Er i-MXenes. This study contributes to identifying pathways for the preparation of pure 2D i-MXenes from RE-based i-MAX phases and developing improved synthesis methods. With additional process optimization, the 2D i-MXene holds a strong potential for electrochemical energy storage applications. Additionally, the electronic structures of Mo 1.33 C were theoretically studied using first-principles density functional theory calculations, which revealed that pristine Mo 1.33 C is metallic, and this metallic nature is preserved even with −O, −F, and mixed functionalization.

chemical structure↗
DOE OSTI2021

Importance of Van der Waals Interactions in Hydrogen Adsorption on a Silicon-carbide Nanotube Revisited with vdW-DFT and Quantum Monte Carlo

Density functional theory (DFT) is a valuable tool for calculating adsorption energies toward designing materials for hydrogen storage. However, dispersion forces being absent from the local/semi-local theory, it remains unclear as to how the consideration of van der Waals (vdW) interactions affects such calculations. For the first time, we applied diffusion Monte Carlo (DMC) to evaluate the adsorption characteristics of a hydrogen molecule on a (5,5) armchair silicon-carbide nanotube (H 2 -SiCNT). Within the DFT framework, we benchmarked various exchange-correlation functionals, including those recently developed for treating dispersion or vdW interactions. We found that the vdW-corrected DFT methods agree well with DMC, whereas the local (semilocal) functional significantly over (under)-binds. Furthermore, we fully optimized the H 2 -SiCNT geometry within the DFT framework and investigated the correlation between the structure and charge density. The vdW contribution to the adsorption was found to be non-negligible at similar to 1 kcal/mol per hydrogen molecule, which amounts to 9-29% of the ideal adsorption energy required for hydrogen storage applications.

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

Computational catalyst discovery: Active classification through myopic multiscale sampling

We report the recent boom in computational chemistry has enabled several projects aimed at discovering useful materials or catalysts. We acknowledge and address two recurring issues in the field of computational catalyst discovery. First, calculating macro-scale catalyst properties is not straightforward when using ensembles of atomic-scale calculations [e.g., density functional theory (DFT)]. We attempt to address this issue by creating a multi-scale model that estimates bulk catalyst activity using adsorption energy predictions from both DFT and machine learning models. The second issue is that many catalyst discovery efforts seek to optimize catalyst properties, but optimization is an inherently exploitative objective that is in tension with the explorative nature of early-stage discovery projects. In other words, why invest so much time finding a “best” catalyst when it is likely to fail for some other, unforeseen problem? We address this issue by relaxing the catalyst discovery goal into a classification problem: “What is the set of catalysts that is worth testing experimentally?” Here, we present a catalyst discovery method called myopic multiscale sampling, which combines multiscale modeling with automated selection of DFT calculations. It is an active classification strategy that seeks to classify catalysts as “worth investigating” or “not worth investigating” experimentally. Our results show an ~7–16 times speedup in catalyst classification relative to random sampling. These results were based on offline simulations of our algorithm on two different datasets: a larger, synthesized dataset and a smaller, real dataset.

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

Exact-two-component block-localized wave function: A simple scheme for the automatic computation of relativistic ΔSCF

Block-localized wave function is a useful method for optimizing constrained determinants. In this article, we extend the generalized block-localized wave function technique to relativistic two-component framework. Optimization of excited state determinants for two-component wave functions presents a unique challenge because the excited state manifold is often quite dense with degenerate states. Furthermore, we test the degree to which certain symmetries result naturally from the ΔSCF optimization such as time reversal symmetry and symmetry with respect to the total angular momentum operator on a series of atomic systems. Variational optimizations may often break symmetry in order to lower the overall energy, just as unrestricted Hartree-Fock breaks spin symmetry. Overall, we demonstrate that time reversal symmetry is roughly maintained when using Hartree-Fock, but less so when using Kohn-Sham density functional theory. Additionally, maintaining total angular momentum symmetry appears to be system dependent and not guaranteed. Finally, we were able to trace the breaking of total angular momentum symmetry to the relaxation of core electrons.

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

Information Content of the Parity-Violating Asymmetry in Pb 208

The parity-violating asymmetry A PV in 208 Pb, recently measured by the PREX-2 Collaboration, is studied using modern relativistic (covariant) and nonrelativistic energy density functionals. We first assess the theoretical uncertainty on A PV which is intrinsic to the adopted approach. To this end, we use quantified functionals that are able to accommodate our previous knowledge on nuclear observables such as binding energies, charge radii, and the dipole polarizability α D of 208 Pb. We then add the quantified value of A PV together with α D to our calibration dataset to optimize new functionals. Based on these results, we predict a neutron skin thickness in 208 Pb r skin = 0.19 ± 0.02 fm and the symmetry-energy slope L = 54 ± 8 MeV. These values are consistent with other estimates based on astrophysical data and are significantly lower than those recently reported using a particular set of relativistic energy density functionals. Here, we also make a prediction for the A PV value in 48 Ca that will be soon available from the CREX measurement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Mechanistic Insights and Rational Design of Ca-Doped CeO 2 Catalyst for Acetic Acid Ketonization

Carboxylic acid ketonization has recently gained significant attention to produce biomass-derived hydrocarbon fuels as it not only removes the highly reactive carboxylic functional group but also increases the size of the carbon chain. In this work, Ca-doped CeO 2 -based catalysts were investigated for acetic acid ketonization using a combined experimental and computational approach. Acetic acid conversion was performed across a range of temperatures including higher temperatures relevant to catalytic hot gas filtration (450 °C). Ca addition slightly decreases overall acetic acid ketonization reactivity yet stabilizes the catalyst at the higher temperatures necessary for catalytic hot gas filtration. From density functional theory calculations of the ketonization reaction mechanism, the C–C coupling and water formation steps are identified as two of the most energy-consuming steps on a CeO 2 surface with a proximal oxygen vacancy and the presence of a Ca dopant stabilizes the key intermediates. Calculations predict an optimal structure comprising three Ca ensembles to minimize the reaction free energies for C–C coupling and water formation steps. These findings provide a priori information to guide future experiments for ketonization catalyst design and development.

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

Redox Potentials for Hydrogen Species in Acetonitrile and Tetrahydrofuran

The properties of the fundamental hydrogen species, proton (H+), hydrogen atom (H•), and hydride anion (H-) are critical to a vast range of chemical processes, yet their thermodynamic properties in nonaqueous solvents are not well established. A hybrid supermolecule-continuum approach is used to predict the Gibbs free energies of solvation (?G°solv) and standard redox potentials (E°) for the 2H+/H2 and H•/H- couples in acetonitrile (MeCN) and tetrahydrofuran (THF) following the approach previously used for water. Gas phase geometry optimizations were done using density functional theory (DFT) with the ?B97XD functional, Møller-Plesset second-order perturbation theory (MP2) with augmented correlation-consistent basis sets and composite-correlated molecular orbital theory (G3MP2 and G3MP2B3). Additional single point calculations in gas phase were performed for H- affinity for MeCN and THF at the coupled cluster CCSD(T) level using MP2 geometries. Solvation was included using the self-consistent reaction field SMD model. For H+, the ?G°solv values are predicted to be -252.2 kcal/mol in MeCN and -261.0 kcal/mol in THF. The redox potentials for the 2H+/H2 couple are +0.49 V in MeCN and +0.11 V in THF relative to the aqueous SHE, consistent with the available experimental data. The ?G°solv for H- is predicted to be -80 kcal/mol in MeCN and -68 in THF, and for H•, ?G°solv is predicted to be -1.3 kcal/mol in THF and -2.0 kcal/mol in MeCN. These solvation energies yield respective calculated redox potentials for the H•/H- couple of -0.14 V in MeCN and -0.68 V in THF relative to the aqueous SHE. H2 solvation is slightly exergonic in MeCN (-3.3 kcal/mol) and near thermoneutral in THF (+1.3 kcal/mol). The pKa(H2) is predicted to be 45.4 in MeCN and 47.8 in THF. The proton-coupled electron transfer (PCET) accounting term CG was evaluated giving values on the ferrocene scale within approximately 2 kcal/mol of literature in both solvents. These results provide benchmarks for hydrogen redox chemistry in MeCN and THF and improve consistency across E°, ?G°solv, pKa, and CG. These results provide benchmarks for hydrogen redox chemistry in nonaqueous media challenging previous estimates based on indirect extrapolations and assumptions.

Duda, Damian P.↗
DOE OSTI2024

Tailoring the Electronic Structures and Spectral Properties of ZnO with Irradiation Defects Generated Under Intense Electronic Excitation: A Combined Experimental and DFT Approach

Here, the relationship between the composition of the internal defect states, spectral properties, and correlated electronic structures of wurtzite zinc oxide (ZnO) crystals under 645 MeV Xe 35+ irradiation is systematically investigated, employing experimental characterizations combined with first-principle calculations. Based on the ion irradiation-induced thermal expansion and relaxation processes, the high concentration of vacancy/interstitial defects produced from the transient disordered phase in molten track states trigger photoelectric changes, as follows: i) the generation of internal defect states effectively reduces the intrinsic bandgap (3.25 eV → 2.66 eV); ii) a large number of defective active sites inhibits the recombination between electron–hole pairs, causing dark conductance and photoconductance to increase with increasing damage levels until optimal fluence is achieved. Based on the density functional theory (DFT) with the GGA + U (GGA = generalized gradient approximation) method, the defective models associated with the different electronic structures, density of states, formation energy, and the nature of the chemical bonding are established. The narrowing of the bandgap observed experimentally and the enhancement of carrier concentration originating from the internal electron defect states are qualitatively verified, therefore laying the foundation for designing future nanoscale photoelectronic devices and microelectronics applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Architecting the Third Dimension of Electrochemical Energy Storage

Three-dimensional (3D) architectural design has emerged as a powerful strategy to push electrochemical energy storage (EES) devices beyond the intrinsic limitations of conventional two-dimensional (2D) electrodes. While planar architectures enable high packing density and mature manufacturing, they suffer from limited ion transport and low active-material loading. In contrast, 3D architectures introduce low-tortuosity networks and high surface area that enhance charge and mass transport while supporting thick, high mass-loading electrodes. However, their practicality remains hindered by challenges in volumetric density, mechanical stability, and large-scale manufacturability. Here, this Perspective examines the key evaluation and design principles that govern 3D device performance. We discuss the fundamental trade-offs between porosity, volumetric density, and mechanical stability that shape 3D design and highlight emerging strategies for integrating materials engineering, structural optimization, device integration, computational modeling, and scalable manufacturing. By aligning structural functionality with manufacturability, 3D architectures can evolve from laboratory prototypes to commercially viable energy storage systems.

25 ENERGY STORAGE↗