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

Influence of the symmetry energy on the nuclear binding energies and the neutron drip line position

A clear connection can be established between properties of nuclear matter and finite-nuclei observables, such as the correlation between the slope of the symmetry energy and the dipole polarizability, or between compressibility and the isoscalar monopole giant resonance excitation energy. Establishing a connection between realistic atomic nuclei and an idealized infinite nuclear matter leads to a better understanding of underlying physical mechanisms that govern nuclear dynamics. In this work, we aim to study the dependence of the binding energies and related quantities (e.g., location of drip lines, the total number of bound even-even nuclei) on the symmetry energy S 2 (ρ). The properties of finite nuclei are calculated by employing the relativistic Hartree-Bogoliubov model, assuming even-even axial and reflection symmetric nuclei. Calculations are performed by employing two families of relativistic energy density functionals, based on different effective Lagrangians, constrained to a specific symmetry energy at the saturation density J within the interval of 30–36 MeV. Nuclear binding energies and related quantities of bound nuclei are calculated between 8 ≤ Z ≤ 104 from the two-proton to the two-neutron drip line. As the neutron drip line is approached, the interactions with stiffer J tend to predict more bound nuclei, resulting in a systematic shift of the two-neutron drip line towards more neutron-rich nuclei. Consequentially, a correlation between the number of bound nuclei N nucl and S 2 (ρ) is established for a set of functionals constrained using the similar optimization procedures. Furthermore the direction of the relationship between the number of bound nuclei and the symmetry energy highly depends on the density under consideration.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Electronic Structure of Metallophthalocyanines, MPc (M = Fe, Co, Ni, Cu, Zn, Mg) and Fluorinated MPc

In this work, we compute the electronic structure and optical excitation energies of metal-free and transition-metal phthalocyanines (H2Pc and MPc for M = Fe, Co, Ni, Cu, Zn, Mg) using density functional theory with optimally tuned range-separated hybrid functionals (OT-RSH). We show that the OT-RSH approach provides photoemission spectra in quantitative agreement with experiments as well as optical band gaps within 10% of their experimental values, capturing the interplay of localized d-states and delocalized pi-pi* states for these organometallic compounds. We examine the tunability of MPcs and H 2 Pc through fluorination, resulting in quasi-rigid shifts of the molecular orbital energies by up to 0.7 eV. Our comprehensive data set provides a new computational benchmark for gas-phase phthalocyanines, significantly improving upon other density-functional-theory-based approaches.

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

An experimental and computational view of the photoionization of diol–water clusters

In the interstellar medium, diols and other prebiotic molecules adsorb onto icy mantles surrounding dust grains. Water in the ice may affect the reactivity and photoionization of these diols. Ethylene glycol (EG), 1,2-propylene glycol, and 1,3-propylene glycol clusters with water clusters were used as a proxy to study these interactions. The diol–water clusters were generated in a continuous supersonic molecular beam, photoionized by synchrotron-based vacuum ultraviolet light from the Advanced Light Source, and subsequently detected by reflectron time-of-flight mass spectrometry. The appearance energies for the detected clusters were determined from the mass spectra, collected at increasing photon energy. Clusters of both diol fragments and unfragmented diols with water were detected. The lowest energy geometry optimized conformers for the observed EG–water clusters and EG fragment–water clusters have been visualized using density functional theory (DFT), providing insight into hydrogen bonding networks and how these affect fragmentation and appearance energy. As the number of water molecules clustered around EG fragments (m/z 31 and 32) increased, the appearance energy for the cluster decreased, indicating a stabilization by water. This trend was supported by DFT calculations. Fragment clusters from 1,2-propylene glycol exhibited a similar trend, but with a smaller energy decrease, and no trend was observed from 1,3-propylene glycol. In conclusion, we discuss and suggest that the reactivity and photoionization of diols in the presence of water depend on the size of the diol, the location of the hydroxyl group, and the number of waters clustered around the diol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗
NASA NTRS2000

Langmuir Probe Measurements in an Inductively Coupled Ar/CF4 Plasmas

Technological advancement in the microelectronics industry requires an understanding of the physical and chemical processes occurring in plasmas of fluorocarbon gases, such as carbon tetrafluoride (CF4) which is commonly used as an etchant, and their mixtures to optimize various operating parameters. In this paper we report data on electron number density (ne), electron temperature'(Te), electron energy distribution function (EEDF), mean electron energy, ion number density (ni), and plasma potential (Vp) measured by using Langmuir probe in an inductively coupled 13.56 MHz radio frequency plasmas generated in 50%Ar:50%CF4 mixture in the GEC cell. The probe data were recorded at various radial positions providing radial profiles of these plasma parameters at 10-50 mTorr pressures and 200 W and 300 W of RF power. Present measurements indicate that the electron and ion number densities increase with increase in pressure and power. Whereas the plasma potential and electron temperature decrease with increase in pressure, and they weakly depend on RF power. The radial profiles exhibit that the electron and ion number densities and the plasma potential peak at the center of the plasma with an exponential fall away from it, while the electron temperature has a minimum at the center and it increases steadily towards the electrode edge. The EEDFs have a characteristic drop near the low energy end at all pressures and pressures and their shapes represent non-Maxwellian plasma and exhibit more like Druyvesteyn energy distribution.v

Rao, M. V. V. S.↗
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↗