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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Safe Electrolytes for Batteries

The demand for reliable and safe lithium-ion batteries (LIBs) for electric vehicles (EVs) and energy storage systems (ESS) necessitates the exploration of nonflammable phosphorus-based electrolytes as alternatives to the traditional flammable carbonate-based ones. However, integrating phosphorus-based electrolyte poses challenges, including capacity fading and compatibility issues with graphite material. Here, a high-throughput (HTP) electrochemical characterization method, similar to PH test paper, is introduced to fast screen compatible phosphorus-based electrolyte for use with graphite. Among 1,740 combinations of phosphorus-based electrolytes and graphite materials, 101 promising combinations are identified for further evaluation. These identified phosphorus-based electrolytes and graphite materials are evaluated and optimized in Li/Graphite half-cells and Graphite/LiFePO 4 (LFP) full-cells using commercial-level electrodes. The desolvation energy of a complex with one Li+ and four solvent molecules of phosphate or carbonate is calculated by density functional theory (DFT). Key parameters such as viscosity, ionic conductivity, and flammability of the electrolytes are thoroughly tested and optimized. The modified phosphate-dominant electrolyte (2 M LiFSI TEP/DME/EC (6/2/2, volume ratio) + 5 wt.% VEC + 5 wt.% FEC) demonstrates excellent thermal stability with lithiated graphite and showcases superior cycling performance in the Graphite/LFP full cell, surpassing prior research findings in phosphate-dominant electrolytes. By the rapid HTP screening and identification of compatible electrolyte-graphite combinations, this approach contributes to expediting the development of safer LIBs for EVs and ESS.

25 ENERGY STORAGE↗
DOE OSTI2021

Optimization and supervised machine learning methods for fitting numerical physics models without derivatives

Here, we address the calibration of a computationally expensive nuclear physics model for which derivative information with respect to the fit parameters is not readily available. Of particular interest is the performance of optimization-based training algorithms when dozens, rather than millions or more, of training data are available and when the expense of the model places limitations on the number of concurrent model evaluations that can be performed. As a case study, we consider the Fayans energy density functional model, which has characteristics similar to many model fitting and calibration problems in nuclear physics. We analyze hyperparameter tuning considerations and variability associated with stochastic optimization algorithms and illustrate considerations for tuning in different computational settings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Surrogate Hessian accelerated structural optimization for stochastic electronic structure theories

In this work, we present an efficient energy-based method for structural optimization with stochastic electronic structure theories, such as diffusion quantum Monte Carlo (DMC). This method is based on robust line-search energy minimization in reduced parameter space, exploiting approximate but accurate Hessian information from a surrogate theory, such as density functional theory. The surrogate theory is also used to characterize the potential energy surface, allowing for simple but reliable ways to maximize statistical efficiency while retaining controllable accuracy. We demonstrate the method by finding the minimum DMC energy structures of the selected flake-like aromatic molecules, such as benzene, coronene, and ovalene, represented by 2, 6, and 19 structural parameters, respectively. In each case, the energy minimum is found within two parallel line-search iterations. The method is near-optimal for a line-search technique and suitable for a broad range of applications. It is easily generalized to any electronic structure method where forces and stresses are still under active development and implementation, such as diffusion Monte Carlo, auxiliary-field Monte Carlo, and stochastic configuration interaction, as well as deterministic approaches such as the random-phase approximation. Accurate and efficient means of geometry optimization could shed light on a broad class of materials and molecules, showing high sensitivity of induced properties to structural variables.

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

Elucidation of Bottom-up Growth of CaCO3 Involving Prenucleation Clusters from Structure Predictions and Decomposition of Globally Optimized (CaCO3)n Nanoclusters

Low-energy minima structures for (CaCO3)n, n = 28, are predicted using bottom-up genetic algorithms in conjunction with density functional theory electronic structure calculations, in comparison with the frozen and relaxed top-down clusters generated by cuts from the calcite, vaterite, and aragonite crystal structures. Similarities in structural motifs for the bottom-up and relaxed top-down are revealed using a fragment recognition technique. Fragment energy decomposition analysis shows that the bottom-up and relaxed top-down clusters belong to two classes of amorphous clusters with distinct intracluster energy distributions, despite their structural similarity. The bottom-up clusters with >20 formula units are surface stabilized with negative surface energy densities. In contrast, the top-down clusters are interior stabilized with positive surface energy densities. We prove that the sign of the surface energy density determines whether the nucleation reaction energy as a function of nuclear size has a maximum or a minimum. The surface-stabilized bottom-up clusters are proposed to be a type of prenucleation cluster at the minimum of the nucleation reaction energy. A mechanism for mineralization of CaCO3 involving prenucleation clusters and nonclassical growth pathway is proposed on the basis of our theoretical findings, which is consistent with previous titration experiments.

calcium carbonate, prenucleation cluster, biominer↗
DOE OSTI2026

Maximizing machine learning interatomic potential transferability for the discovery of the novel stellated octadecagon Bi18-Pt24 cage structure

Achieving true transferability remains the central challenge for Machine Learning Interatomic Potentials (ML-IAPs) in modeling complex bimetallic nanoclusters across their vast potential energy surfaces. We systematically investigate data selection strategies to optimize the Chebyshev Interaction Model for Efficient Simulation (ChIMES) potential for the Bi-Pt nanoclusters by comparing three innovative sampling methods: Principal Component Analysis (PCA)/k-means (structural diversity), t-distributedStochasticNeighborEmbedding (t-SNE)/k-means (force-space diversity), and hierarchical clustering. Quantitatively, the PCA/k-means strategy proved most effective for global accuracy, yielding the lowest force errors and achieving energy root mean square errors (RMSE) values competitive with Density Functional Theory (DFT), demonstrating excellent accuracy (19.16meV/atom). Structural validation on 34 unique DFT-optimized isomers further confirmed the potential’s high fidelity, with the best model PCA/k-means reproducing structures with an average root mean square deviation (RMSD) of 0.10 Å. However, the t-SNE methods, by maximizing diversity in the force space, demonstrated superior extrapolative power, leading to the more precise prediction of a novel stellated octadecagon Bi18⁢Pt24 cage structure, demonstrating the potential for exploring previously unseen morphologies. Our results establish a clear methodology for strategic data sampling that successfully maximizes ML-IAP transferability, providing an accurate and computationally efficient tool that accelerates the theoretical discovery of complex bimetallic architectures.

Vangheluwe, Raphaël [Université Paris-Saclay, CNRS↗
DOE OSTI2026

Structural effects of composition tuning in A-site disordered perovskite La0.5Li0.5−xMxTiO3 (M = Na, K) nanorods for fast interfacial transport for solid composite electrolyte design

Polymer-ceramic composite electrolytes (CPE) have emerged as promising replacements for currently used liquid organic electrolytes, which pose as safety hazards, in Li-ion batteries. Limits in conductivity enhancement in CPEs is often attributed to a high resistance for Li-ion transport along the interface due to the incompatibility of the polymer and ceramic phases. A clear understanding of the interfacial structure and how this impacts interfacial Li-ion transport is needed in order to efficiently design a CPE with optimal ionic conductivity. In this study, density functional theory (DFT) calculations, in conjunction with scanning transmission electron microscopy (STEM) and electron energy loss spectroscopy (EELS), are used to unveil the bulk and surface structure of La0.5Li0.5−xMxTiO3 (M = Na, K) (LMTO) nanorods, a newly-reported ceramic system with promising applications in CPEs, and the LMTO interactions with the poly(lithium sulfonyl (trifluoromethane sulfonyl)imide methacrylate) (p(MTFSILi)) polymerized ionic liquid. Substitution of lithium for sodium and potassium onto the A-site perovskite lattice is shown to increase the stability of the preferred pseudocubic perovskite phase because of their increased cation size which allows them to reduce the TiO6 octahedral rotations in the bulk and at the surface. STEM and EELS results show that LMTO nanorods with differing compositions have Ti-enriched (110)-oriented surfaces. Increased sodium and potassium compositions in LMTO is also shown using the DFT calculations to weaken the binding of the lithium atom to the MTFSI unit when LiMTFSI is adsorbed to the LMTO surface. This weakening of lithium binding to the polymer at the LMTO interface indicates that the lithium mobility is increased, correlating to an increase in interfacial Li-ion transport. Overall, this work provides insights into how to tune the interfacial interactions between the polymer (p(MTFSI)) and ceramic (LMTO) for optimal CPE design.

Shepard, Lauren B [Pennsylvania State University, ↗
DOE OSTI2021

Tuning the electronic structure of Ag-Pd alloys to enhance performance for alkaline oxygen reduction

Alloying is a powerful tool that can improve the electrocatalytic performance and viability of diverse electrochemical renewable energy technologies. Herein, we enhance the activity of Pd-based electrocatalysts via Ag-Pd alloying while simultaneously lowering precious metal content in a broad-range compositional study focusing on highly comparable Ag-Pd thin films synthesized systematically via electron-beam physical vapor co-deposition. Cyclic voltammetry in 0.1 M KOH shows enhancements across a wide range of alloys; even slight alloying with Ag (e.g. Ag 0.1 Pd 0.9 ) leads to intrinsic activity enhancements up to 5-fold at 0.9 V vs. RHE compared to pure Pd. Based on density functional theory and x-ray absorption, we hypothesize that these enhancements arise mainly from ligand effects that optimize adsorbate–metal binding energies with enhanced Ag-Pd hybridization. This work shows the versatility of coupled experimental-theoretical methods in designing materials with specific and tunable properties and aids the development of highly active electrocatalysts with decreased precious-metal content.

25 ENERGY STORAGE↗
DOE OSTI2021

Density Functional Theory Solute-Dislocation Interactions

DFT-optimized structures of dislocations with a solute at different sites. The dataset include the VASP input files (except POTCAR) and output files OSZICAR and OUTCAR. There are also text files containing the DFT-computed interaction energies. The README file contains expressions that can be used to compute the interaction energies and chemical misfit energies using the volumetric size misfit for any solute species in Mg.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Properties of La 0.7 Ca 0.3 MnO 3 under extreme tensile strain

The complex phase diagram of manganites with simultaneously active spin, charge, orbital, and lattice degrees of freedom continues providing surprises. In a recent groundbreaking experiment, membranes of the perovskite manganite La 0.7 Ca 0.3 MnO 3 (LCMO) deposited on a flexible polymer layer were strained up to 8% [S. S. Hong et al., Science 368, 71 (2020)], much more than achieved by regular strain induced by a rigid substrate. By increasing this strain, a metal-insulator transition was reported. Here we reproduce the results of the experiments using Monte Carlo simulations of the two-orbital double-exchange model including Jahn-Teller distortions at hole density x=1/3. The full phase diagram with varying temperature and Jahn-Teller coupling λ is presented. When the bandwidth W of mobile electrons is reduced, thus when the effective Jahn-Teller coupling λ/W is increased, a metal-insulator transition is found in our simulations, between a ferromagnetic metallic state with uniform charge distribution and an insulator with diagonal charge stripes that retains its ferromagnetic character. In between the hole-rich diagonals, staggered orbital order occurs. We also report resistivity and magnetization measurements alongside with spin correlations and charge structure factors. Our overall conclusions are in agreement with the recent experimental and density functional theory results by Hong et al., and we confirm much earlier ground state predictions of striped ferromagnetic order using energy optimization techniques by T. Hotta et al. [Phys. Rev. Lett. 86, 4922 (2001)]. The experimental observation of one of the states predicted by theory suggests that diagonal stripes could be achieved at other hole densities as well, such as x=1/4, if LCMO membranes with that hole doping were subject to similar strains.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Hydrogen Bond Benchmark: Focal‐Point Analysis and Assessment of DFT Functionals

We performed a hierarchical, convergent ab initio benchmark study and systematically analyzed the performance of density functional approximations for describing hydrogen bonds in small neutral, cationic, and anionic complexes, as well as in larger systems involving amide, urea, deltamide, and squaramide moieties. Focal point analyses (FPA), extrapolating to the ab initio limit, were carried out using correlated wave function methods up to CCSDT(Q) for the small complexes and CCSD(T) for the larger systems, together with correlation-consistent Gaussian basis sets up to the complete basis set limit. Optimized geometries and vibrational frequencies were obtained at the CCSD(T) level. The resulting FPA hydrogen-bond energies converge within a few tenths of a kcal mol −1 . These reference data were used to evaluate 60 density functionals (including 12 dispersion-corrected), spanning the local-density approximation (LDA), generalized gradient approximations (GGAs), meta-GGAs, hybrids, meta-hybrids, double-hybrids, and range-separated hybrids. Overall, the meta-hybrid M06-2X provides the best performance for both hydrogen bond energies and geometries, while the dispersion-corrected GGAs BLYP-D3(BJ) and BLYP-D4 also yield accurate hydrogen-bond data and can serve as cost-effective options for studying large and complex systems.

coupled cluster theory↗
DOE OSTI2020

Spectral-density estimation with the Gaussian integral transform

The spectral-density operator $\hat{ρ}(ω) = δ(ω–\hat{H})$ plays a central role in linear response theory as its expectation value, the dynamical response function, can be used to compute scattering cross sections. In this work, we describe a near optimal quantum algorithm providing an approximation to the spectral density with energy resolution $\Delta$ and error $\epsilon$ using $O(\sqrt{\text{log}_2 (1/ε)[\text{log}_2 (1 / Δ) + \text{log}_2 (1/ε)]/ Δ)}$ operations. This is achieved without using expensive approximations to the time-evolution operator, but instead exploiting qubitization to implement an approximate Gaussian integral transform of the spectral density. Finally, we also describe appropriate error metrics to assess the quality of the spectral function approximations more generally.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Increasing Compressed Gas Energy Storage Density Using CO2–N2 Gas Mixture

This paper demonstrates a new method by which the energy storage density of compressed air systems is increased by 56.8% by changing the composition of the compressed gas to include a condensable component. A higher storage density of 7.33 MJ/m3 is possible using a mixture of 88% CO2 and 12% N2 compared to 4.67 MJ/m3 using pure N2. This ratio of gases representing an optimum mixture was determined through computer simulations that considered a variety of different proportions from pure CO2 to pure N2. The computer simulations are based on a thermodynamic equilibrium model that predicts the mixture composition as a function of volume and pressure under progressive compression to ultimately identify the optimal mixture composition (88% CO2 + 12% N2). The model and simulations predict that the optimal gas mixture attains a higher energy storage density than using either of the pure gases.

25 ENERGY STORAGE↗
DOE OSTI2024

Mitigating Green's function Monte Carlo signal-to-noise problems using contour deformations

The Green's function Monte Carlo (GFMC) method provides accurate solutions to the nuclear many-body problem and predicts properties of light nuclei starting from realistic two- and three-body interactions. Controlling the GFMC fermion sign problem is crucial, as the signal-to-noise ratio decreases exponentially with imaginary time, requiring significant computing resources. Inspired by similar scenarios in lattice quantum field theory and spin systems, in this work, we employ integration contour deformations to improve the GFMC signal-to-noise ratio. Machine learning techniques are used to select optimal contours with minimal variance from parametrized families of deformations. As a proof of principle, we consider the deuteron binding energies and Euclidean density response functions. Here, we only observe mild signal-to-noise improvement for the binding energy case. On the other hand, we achieve an order of magnitude reduction of the variance for Euclidean density response functions, paving the way for computing electron- and neutrino-nucleus cross sections of larger nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Mitigating Green's function Monte Carlo signal-to-noise problems using contour deformations

The Green's function Monte Carlo (GFMC) method provides accurate solutions to the nuclear many-body problem and predicts properties of light nuclei starting from realistic two- and three-body interactions. Controlling the GFMC fermion-sign problem is crucial, as the signal-to-noise ratio decreases exponentially with Euclidean time, requiring significant computing resources. Inspired by similar scenarios in lattice quantum field theory and spin systems, in this work, we employ integration contour deformations to improve the GFMC signal-to-noise ratio. Machine learning techniques are used to select optimal contours with minimal variance from parameterized families of deformations. As a proof of principle, we consider the deuteron binding energies and Euclidean density response functions. We only observe mild signal-to-noise improvement for the binding energy case. On the other hand, we achieve an order of magnitude reduction of the variance for Euclidean density response functions, paving the way for computing electron- and neutrino-nucleus cross-sections of larger nuclei.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

An Analytical Tool to Evaluate Defect Thermodynamics of (La,Ba)Fe1-xMxO3-δ Perovskites for Solid-Oxide Cell Applications

A modeling tool of the defect thermodynamics of (La,Ba)Fe1-xMxO3-δ perovskites which includes energetic information about oxygen vacancy formation, hydration, hydride formation, and charge disproportionation reactions has been developed1. This tool incorporates defect energies and entropies expressed as sixth order polynomial functions to allow refinements of the defect reaction equilibrium constants in the thermodynamic analysis. Calculation of (La,Ba)Fe1-xMxO3-δ Brouwer diagrams as a function of pO2/pH2O and pH2/H2O in a range of temperatures of interest is facilitated by this modeling tool. The results obtained can provide direct guidance how the electronic and ionic defect concentrations of the triple conducting perovskite materials can be used to optimize performance of solid oxide cells for energy applications. The impact of magnetic and electronic structures of the perovskites on the defect reaction energies and entropies as obtained from density function theory modeling and the role played by hydride defect species will also be discussed.

Lee, Yueh-Lin↗
DOE OSTI2023

Spectral-partitioned Kohn-Sham density functional theory

Here we introduce a general, variational scheme for systematic approximation of a given Kohn-Sham free-energy functional by partitioning the density matrix into distinct spectral domains, each of which may be spanned by an independent diagonal representation without requirement of mutual orthogonality. It is shown that by generalizing the entropic contribution to the free energy to allow for independent representations in each spectral domain, the free energy becomes an upper bound to the exact (unpartitioned) Kohn-Sham free energy, attaining this limit as the representations approach Kohn-Sham eigenfunctions. A numerical procedure is devised for calculation of the generalized entropy associated with spectral partitioning of the density matrix. The result is a powerful framework for Kohn-Sham calculations of systems whose occupied subspaces span multiple energy regimes. As a case in point, we apply the proposed framework to warm- and hot-dense matter described by finite-temperature density functional theory, where at high energies the density matrix is represented by that of the free-electron gas, while at low energies it is variationally optimized. We derive expressions for the spectral-partitioned Kohn-Sham Hamiltonian, atomic forces, and macroscopic stresses within the projector-augmented wave (PAW) and the norm-conserving pseudopotential methods. It is demonstrated that at high temperatures, spectral partitioning facilitates accurate calculations at dramatically reduced computational cost. Moreover, as temperature is increased, fewer exact Kohn-Sham states are required for a given accuracy, leading to further reductions in computational cost. Finally, it is shown that standard multiprojector expansions of electronic orbitals within atomic spheres in the PAW method lack sufficient completeness at high temperatures. Spectral partitioning provides a systematic solution for this fundamental problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Engineering Thermally Resilient and Kinetically Active Reversible Protonic Ceramic Cells via Interfacial Design

Achieving concurrent fast electrode kinetics and long-term thermo-mechanical durability remains a critical challenge for reversible protonic ceramic electrochemical cells (R-PCECs). Herein, we report an interfacial engineering strategy that integrates a perovs.kite-type PrBaRu0.1Co1.9O5+δ (PBRC) nanoparticle layer onto a PrBa0.5Sr0.5Co1.5Fe0.5O5+δ (PBSCF) substrate (PBRC-PBSCF), together with a modified pellet-assisted sintering approach to fabricate dense BaZr0.4Ce0.4Y0.1Yb0.1O3-δ (BZCYYb4411) electrolytes. The in situ reconstructed heterointerface enhances oxygen reduction/evolution reaction (ORR/OER) kinetics, promotes H2O adsorption/dissociation, and improves steam tolerance, as verified by electrochemical measurements and interfacial microstructural analyses. Density functional theory reveals that Ru-induced electronic modulation at the PBRC-PBSCF interface lowers the energy of oxygen vacancy formation and optimizes the position of the O 2p band center, thereby accelerating oxygen redox kinetics and stabilizing the interface. The resulting R-PCECs deliver an excellent peak power density of 1.112 W cm−2 and an electrolysis current density of −1.257 A cm−2 at 1.3 V in 3% H2O wet air at 600°C, with a reasonable faradaic efficiency. Furthermore, the cells demonstrate excellent stability, sustaining 100 h of thermal cycling (400–600°C, 200°C h−1) in both fuel cell and electrolysis modes, with 600 h of stability in electrolysis mode (600°C, −0.5 to −2 A cm−2).

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2020

Predicting Excitation Energies of Twisted Intramolecular Charge-Transfer States with the Time-Dependent Density Functional Theory: Comparison with Experimental Measurements in the Gas Phase and Solvents Ranging from Hexanes to Acetonitrile

Electronically excited states characterized by intramolecular charge transfer play an essential role in many biological processes and optical devices. The ability to make quantitative ab initio predictions of the relative energetics involved is a challenging yet desirable goal, especially for large molecules in solution. In this work, we present a data set of 61 experimental measurements of absorption and emission processes, both in the gas phase and in solvents representing a broad range of polarities, which involve intramolecular charge transfer mediated by a nonzero, "twisted" dihedral angle between one or more donor and acceptor subunits. Among a variety of density functionals investigated within the framework of linear-response theory, the "optimally tuned" LRC-ωPBE functional, which utilizes a system-specific yet nonempirical procedure to specify the range-separation parameter, emerges as the preferred choice. For the entire set of excitation energies, involving changes in dipole moment ranging from 4 to >20 Debye, the mean signed and absolute errors are 0.02 and 0.18 eV, respectively (compared, e.g., to -0.30 and 0.30 for PBE0, 0.44 and 0.47 for LRC-ωPBEh, 0.83 and 0.83 for ωB97X-V). We analyze the performance of polarizable continuum solvation models available in Q-Chem that partition the solvent response into fast and slow time scales, and clear trends emerge when measurements corresponding to the four small 4-(dimethylamino)benzonitrile (DMABN)-like molecules and a charged species are excluded. We make the case that the large errors found only for small molecules in the gas phase and weak solvents cannot be expected to improve via the optimal tuning procedure, which enforces a condition that is exact only in the well-separated donor-acceptor limit, and present empirical evidence implicating the outsized importance for small donor-acceptor systems of relaxation effects that cannot be accounted for by the linear-response time-dependent density functional theory within the adiabatic approximation. Finally, we demonstrate the utility of the optimally tuned density functional approach by targeting the charge-transfer states of a large biomimetic model system for light-harvesting structures in Photosystem II.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗