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 181 records · Page 10

Biaxial Strains Mediated Oxygen Reduction Electrocatalysis on Fenton Reaction Resistant L1 0 ‐PtZn Fuel Cell Cathode

Abstract PtM alloy catalysts (e.g., PtFe, PtCo), especially in an intermetallic L1 0 structure, have attracted considerable interest due to their respectable activity and stability for the oxygen reduction reaction (ORR) in proton exchange membrane fuel cells (PEMFCs). However, metal‐catalyzed formation of ·OH from H 2 O 2 (i.e., Fenton reaction) by Fe‐ or Co‐containing catalysts causes severe degradation of PEM/catalyst layers, hindering the prospects of commercial applications. Zinc is known as an antioxidant in Fenton reaction, but is rarely alloyed with Pt owing to its relatively negative redox potential. Here, sub‐4 nm intermetallic L1 0 ‐PtZn nanoparticles (NPs) are synthesized as high‐performance PEMFC cathode catalysts. In PEMFC tests, the L1 0 ‐PtZn cathode achieves outstanding activity (0.52 A mg Pt −1 at 0.9 V iR ‐free , and peak power density of 2.00 W cm −2 ) and stability (only 16.6% loss in mass activity after 30 000 voltage cycles), exceeding the U.S. DOE 2020 targets and most of the reported ORR catalysts. Density function theory calculations reveal that biaxial strains developed upon the disorder‐order (A1L1 0 ) transition of PtZn NPs would modulate the surface PtPt distances and optimize PtO binding for ORR activity enhancement, while the increased vacancy formation energy of Zn atoms in an ordered structure accounts for the improved stability.

Liang, Jiashun↗

GW with hybrid functionals for large molecular systems

A low-cost approach for stochastically sampling static exchange during time-dependent Hartree–Fock-type propagation is presented. This enables the use of an excellent hybrid density functional theory (DFT) starting point for stochastic GW quasiparticle energy calculations. Generalized Kohn–Sham molecular orbitals and energies, rather than those of a local-DFT calculation, are used for building the Green function and effective Coulomb interaction. The use of an optimally tuned hybrid diminishes the starting point dependency in one-shot stochastic GW, effectively avoiding the need for self-consistent GW iterations.

Chemistry↗

Ab initio electronic density in solids by many-body plane-wave auxiliary-field quantum Monte Carlo calculations

We present accurate many-body results of the electronic densities in several solid materials, including Si, NaCl, and Cu. These results are obtained using the ab initio auxiliary-field quantum Monte Carlo (AFQMC) method working in a plane-wave basis with norm-conserving, multiple-projector pseudopotentials. AFQMC has been shown to be an excellent many-body total energy method. Computation of observables and correlation functions other than the ground-state energy requires back-propagation, whose adaption and implementation in the plane-wave basis AFQMC framework are discussed in the present paper. This development allows us to compute correlation functions, electronic densities, and interatomic forces, paving the way for geometry optimizations and calculations of thermodynamic properties in solids. Finite supercell size effects are considerably more subtle in the many-body framework than in independent-electron calculations. We analyze the convergence of the electronic density, and obtain best estimates for the thermodynamic limit. The densities from several typical density functionals are benchmarked against our near-exact results. The electronic densities we obtained can also be used to help construct improved density functionals.

36 MATERIALS SCIENCE↗

Kohn-Sham Solver (KSSOLV) v2.0

KSSOLV is a MATLAB toolbox for solving Kohn-Sham density functional theory based electronic structure eigenvalue problems. It uses an object oriented features of MATLAB to represent atom, molecules, wavefunctions and Hamiltonians and their operations. It is designed to make it easier for users to prototype and test new algorithms for solving the Kohn-Sham problem. KSSOLV2.0 contains significant improvement over the original KSSOLV described in a paper published in ACM Transaction on Mathematical Software (attached). In addition to performing ground state calculation for small molecules, it can also perform geometry optimization for both molecules and solids. It uses standard pseudopotentials and implements local density approximation, generalized gradient approximation and hybrid functionals. Future releases will also include time-dependent DFT and post DFT calculations such as the GW quasi-particle energy calculation and Bethe-Salpeter equation solver for optical absorption.

Yang, Chao↗

Machine-Learned Force Field for Molecular Dynamics Simulations of Nonequilibrium Ammonia Synthesis on Iron Catalysts

Ammonia (NH 3 ) is one of the most important industrial chemicals. The conventional NH 3 synthesis method-the Haber–Bosch process-converts atmospheric nitrogen (N 2 ) into NH 3 using H 2 with an iron (Fe) catalyst. However, this process requires high pressures (100–200 atm) and temperatures (700–800 K) near thermal equilibrium. Recently, Fe-based nanocatalysts have been reported to produce promising NH 3 yields under atmospheric pressures and temperature-modulated nonequilibrium conditions. Understanding the mechanism of nonequilibrium catalysis with programmed temperature variation could help to optimize this fully electrified and less energy-intensive process. Although reactive molecular dynamics (RMD) simulations can be a useful tool to model nonequilibrium catalytic processes, they require the development of accurate force fields (i.e., interatomic potentials). Here, we present a machine-learned (ML) force field within the Deep Potential MD (DPMD) framework, trained using periodic density functional theory (DFT) calculations, to model NH 3 synthesis on Fe catalysts with various surface adsorbates such as *N, *H, *N 2 , *H 2 , *NH, *NH 2 , and *NH 3 . Here, we generated the DFT data from static models of elementary reactions on the most stable (110) surface of body-centered cubic Fe, which then were augmented by data from constant number of particles–volume–temperature (NVT) DFT-MD trajectories at various temperatures. Finally, we utilized the fully optimized ML force field to investigate reaction dynamics at an Fe(110) surface at linearly increasing temperatures using NVT-DPMD simulations. Our simulations indicate that pulsed temperature ramping could prove favorable for NH3 synthesis. For example, we conducted ramping under multiple sets of conditions: (i) from 900 to 1200 K over periods of 0.1–0.3 ns for Fe surfaces precovered with N or NH along with H; and (ii) from 300 to 600 K over 0.1–0.3 ns for Fe surfaces precovered with NH 3 . While our simulations so far are limited to short time scales (very rapid heating), these observations shed light on the mechanism of the high NH 3 synthesis rate achieved in a novel temperature-modulated nonequilibrium catalytic reactor using pulsed heating and cooling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Binding of uranyl cations to a Zr-based metal-organic framework by density functional theory

Zr-based metal–organic frameworks (Zr-MOFs) have been widely used as ion adsorbents for the removal or extraction of toxic and/or radionuclide species from aqueous solutions. However, the mechanisms by which uranyl cations (UO 2 2+ ) interact with Zr-MOFs have not been established. In this work, the nature of the bonding of uranyl cations with a Zr-MOF was determined using density functional theory for nineteen structurally distinct candidate complexes. Further, the results showed that in all cases the binding energy was of the order of 1.5 eV, but depended on the specific bonding site. The most stable structure involved coordination of the uranyl cation and two structurally distinct oxygens in the Zr-MOF metal node. It was also found that higher degree of deprotonation in Zr-MOF correlated with higher binding energy between the Zr-MOF and uranyl cations. These insights can aid in the design of Zr-MOFs with optimized features for efficient capture of uranyl cations.

36 MATERIALS SCIENCE↗

Taming the fixed-node error in diffusion Monte Carlo via range separation

By combining density-functional theory (DFT) and wave function theory via the range separation (RS) of the interelectronic Coulomb operator, we obtain accurate fixed-node diffusion Monte Carlo (FN-DMC) energies with compact multi-determinant trial wave functions. In particular, we combine here short-range exchange-correlation functionals with a flavor of selected configuration interaction known as configuration interaction using a perturbative selection made iteratively (CIPSI), a scheme that we label RS-DFT-CIPSI. One of the take-home messages of the present study is that RS-DFT-CIPSI trial wave functions yield lower fixed-node energies with more compact multi-determinant expansions than CIPSI, especially for small basis sets. Indeed, as the CIPSI component of RS-DFT-CIPSI is relieved from describing the short-range part of the correlation hole around the electron-electron coalescence points, the number of determinants in the trial wave function required to reach a given accuracy is significantly reduced as compared to a conventional CIPSI calculation. Importantly, by performing various numerical experiments, we evidence that the RS-DFT scheme essentially plays the role of a simple Jastrow factor by mimicking short-range correlation effects, hence avoiding the burden of performing a stochastic optimization. Considering the 55 atomization energies of the Gaussian-1 benchmark set of molecules, we show that using a fixed value of mu = 0.5 bohr(-1) provides effective error cancellations as well as compact trial wave functions, making the present method a good candidate for the accurate description of large chemical systems.

Scemama, Anthony↗

Accelerating the Structure Exploration of Diverse Bi–Pt Nanoclusters via Physics‐Informed Machine Learning Potential and Particle Swarm Optimization

Bimetallic Bi–Pt nanoclusters exhibit diverse structural motifs, including core-shell, Janus, and mixed alloy configurations, due to the unique bonding characteristics between Bi and Pt atoms. Using density functional theory refinements from ChIMES physically machine-learned potential and CALYPSO particle swarm optimization global searches, 34 Bi20-Pt20 nanoclusters are systematically classified. The results reveal that Bi atoms predominantly occupy surface sites, driven by charge transfer effects. Cohesive energy trends alone prove insufficient for structure differentiation, necessitating a data-driven approach employing principal component analysis and K-means clustering. Furthermore, vibrational, electronic, and infrared spectral analyses provide additional insights into structure-property relationships. The findings offer an original framework for the automated classification and analysis of bimetallic nanoclusters, enhancing the understanding of their stability and functional properties.

bimetallic nanoparticles↗

Superconductivity, charge density waves, and bipolarons in the Holstein model

The electron-phonon ( e-ph ) interaction remains of great interest in condensed matter physics and plays a vital role in realizing superconductors, charge density waves (CDW), and polarons. We study the two-dimensional Holstein model for e-ph coupling using determinant quantum Monte Carlo across a wide range of its phase diagram as a function of temperature, electron density, dimensionless e-ph coupling strength, and the adiabatic ratio of the phonon frequency to the Fermi energy. We describe the behavior of the CDW correlations, the competition between superconducting and CDW orders and polaron formation, the optimal conditions for superconductivity, and the transition from the weak-coupling regime to the strong-coupling regime. Superconductivity is optimized at intermediate e-ph coupling strength and intermediate electron density, and the superconducting correlations increase monotonically with phonon frequency. The global maximum for superconductivity in the Holstein model occurs at large phonon frequency, the limit where an attractive Hubbard model effectively describes the physics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Crystal Structure Prediction of Binary Alloys via Deep Potential

Predicting crystal structure has been a challenging problem in physics and materials science for a long time. A reliable energy calculation engine combined with an efficient global search algorithm, such as particle swarm optimization algorithm or genetic algorithm, is needed to conduct crystal structure prediction. In recent years, machine learning-based interatomic potential energy surface models have been proposed, potentially allowing us to perform crystal structure prediction for systems with the accuracy of density functional theory (DFT) and the speed of empirical force fields. In this paper, we employ a previously developed Deep Potential model to predict the intermetallic compound of the aluminum–magnesium system, and find six meta-stable phases with negative or nearly zero formation energy. In particular, Mg 12 Al 8 shows excellent ductility and Mg 5 Al 27 has a high Young's modulus. Based on our benchmark results, we propose a relatively robust structure screening criterion that selects potentially stable structures from the Deep Potential-based convex hull and performs DFT refinement. By using this criterion, the computational cost needed to construct the convex hull with ab initio accuracy can be dramatically reduced.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Linking Perfluorosulfonic Acid Ionomer Chemistry and High-Current Density Performance in Fuel-Cell Electrodes

Transport phenomena are key in controlling the performance of electrochemical energy-conversion technologies and can be highly complex, involving multiple length scales and materials/phases. Material designs optimized for one reactant species transport however may inhibit other transport processes. We explore such trade-offs in the context of polymer-electrolyte fuel-cell electrodes, where ionomer thin films provide the necessary proton conductivity but retard oxygen transport to the Pt reaction site and cause interfacial resistance due to sulfonate/Pt interactions. We examine the electrode overall gas-transport resistance and its components as a function of ionomer content and chemistry. Low-equivalent-weight ionomers allow better dissolved-gas and proton transport due to greater water uptake and low crystallinity but also cause significant interfacial resistance due to the high density of sulfonic acid groups. These effects of equivalent weight are also observed via in situ ionic conductivity and CO displacement measurements. Of critical importance, the results are supported by ex situ ellipsometry and X-ray scattering of model thin-film systems, thereby providing direct linkages and applicability of model studies to probe complex heterogeneous structures. Structural and resultant performance changes in the electrode are shown to occur above a threshold sulfonic-group loading, highlighting the significance of ink-based interactions. Furthermore, our findings and methodologies are applicable to a variety of solid-state energy-conversion devices and material designs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Screened Range-Separated Hybrid Functional with Polarizable Continuum Model Overcomes Challenges in Describing Triplet Excitations in the Condensed Phase Using TDDFT

Long range-corrected (LRC) or range-separated hybrid (RSH) functionals where the long-range (LR) limit of electronic interactions is set to the exact exchange have been shown to correct the tendency of traditional density functional theory (DFT) to underestimate the frontier orbital gap. Consequently, the use of such functionals in calculating electronic excited states using linear response based time-dependent DFT (TDDFT) has been successful in correcting the tendency for underestimating the energies of charge transfer states by DFT-based calculations. More recently formulations of functionals that attenuate the LR limit to address condensed-phase effects to polarize the electronic density have been reported. In particular screened RSH (SRSH) combined with polarizable continuum model (PCM) was benchmarked successfully in reproducing the fundamental gap and charge transfer state energies of molecular systems in the condensed phase. Here we use SRSH-PCM to address triplet excited states, and show its success in obtaining correspondence of the low-lying triplet states to the singlet–triplet gap in a similar way that the fundamental orbital gap corresponds to electron removal and addition energies. Importantly, the accuracy of the SRSH-PCM in calculating triplet excitations stands on the polarization consistent framework in addressing the scalar dielectric constant and without affecting the optimal tuning by triplet energies. Finally, the prospect of even further improving the SRSH-PCM accuracy in calculating triplet states can be achieved by optimal tuning on the basis of the spin multiplicity gap.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Can classical DEM simultaneously capture compressibility and flowability of milled biomass?

Accurate prediction of the rheological behavior of biomass is essential for the design and operation of hoppers, feeders, and storage systems in biorefineries. This study examines whether the classical, coarse-grained discrete element method (DEM) formulation can simultaneously reproduce the compressibility and flowability of milled herbaceous biomass, using Miscanthus × giganteus as a representative material. The model represents particles as rigid spheres interacting through Hertz-Mindlin elastic-frictional contacts augmented with an area-dependent cohesion term. Laboratory cyclic compression and wedge-shaped hopper discharge experiments were used as calibration benchmarks. Although the model can independently reproduce each behavior by appropriately tuning particle Young's modulus E and cohesion energy density k, an extensive parametric investigation comprising more than 600 simulations reveals that the optimal parameter regions for compression and hopper flow are distinct and non-overlapping in (E, k) space. Surrogate surface analysis further shows that the corresponding objective-function valleys exhibit similar trends but are approximately parallel and spatially offset, precluding a unified calibration within the explored domain. Sensitivity analysis indicates that compressibility is governed predominantly by stiffness and cohesion, whereas the slope of the mass flow rate-opening relation in hopper discharge is primarily controlled by tangential friction. Extensions incorporating particle size distribution and clumped-sphere representations do not eliminate the incompatibility. These results systematically reveal, for the first time, the structural limitation of simplified DEM formulations in representing biomass rheological behavior, underscoring the necessity for models incorporating additional physical mechanisms, such as particle deformability or enhanced interlocking, to achieve unified predictive capability for biomass handling behavior.

09 BIOMASS FUELS↗

Achieving industrial ammonia synthesis rates at near-ambient conditions through modified scaling relations on a confined dual site

The production of ammonia through the Haber–Bosch process is regarded as one of the most important inventions of the 20th century. Despite significant efforts in optimizing the process, it still consumes 1 to 2% of the worldwide annual energy for the high working temperatures and pressures. The design of a catalyst with a high activity at milder conditions represents another challenge for this reaction. Herein, we combine density functional theory and microkinetic modeling to illustrate a strategy to facilitate low-temperature and -pressure ammonia synthesis through modified energy-scaling relationships using a confined dual site. Our results suggest that an ammonia synthesis rate two to three orders of magnitude higher than the commercial Ru catalyst can be achieved under the same reaction conditions with the introduction of confinement. Such strategies will open pathways for the development of catalysts for the Haber–Bosch process that can operate at milder conditions and present more economically viable alternatives to current industrial solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analysis and prediction of reaction kinetics using the degree of rate control

“Degree of rate control” (DRC) analysis provides a quantitative approach for analysing the kinetics of multi-step reaction mechanisms that has been widely applied to both heterogeneous and homogeneous catalysis research, as well as electrocatalysis. The DRC of any given transition state or intermediate is defined as a partial derivative such that it approximately equals the fractional increase in net rate to the product of interest per differential decrease in its standard-state free energy for that species (÷RT), holding constant the standard-state free energies of all other transition states and intermediates. Even very complex mechanisms usually have only a few species with non-zero DRCs and are thus the species whose interactions with the catalyst most strongly affect the net rate. These key DRC values thus offer a simple and intuitive route to optimize catalyst materials, especially with the assistance of computational methods like density functional theory (DFT). These high-DRC species are also the species whose energetics must be most accurately measured or calculated to achieve an accurate kinetic model for any reaction mechanism. In simple cases with a single “rate-determining step”, the DRC for its transition state (TS) is + 1. Catalyst-bound intermediates, on the other hand, often have negative DRCs equal to a small integer times their fractional occupancy of catalyst sites. The apparent activation energy equals a weighted average of the standard-state enthalpies (relative to reactants) of all the species (intermediates, transition states and products) in the reaction mechanism, each weighted by its DRC (+RT). It has been shown that the apparent transfer coefficient in electrocatalysis, an inverted form of the Tafel slope, is a weighted average of the number of electrons transferred to generate each intermediate or product species in the mechanism, weighted again by the DRC. Quantitative analysis of kinetic isotope effects (KIEs, or the ratio of net rates for different reactant isotopes) in complex mechanisms has shown that the logarithm of the KIE equals the weighted average over all species in the mechanism of the difference between the two isotopes in their standard-state free energies (÷RT), again weighted by the DRC. The reaction orders with respect to fluid-phase concentrations of reactants, products and intermediates have also been proven to be directly related to DRCs. Thus, there are numerous experimental observables which equate to short linear combinations of DRCs, so that combinations of experimental measurements might provide access to DRC values. Since its invention for transition states in 1994 and its generalization to include intermediates in 2009, the DRC has thus far mainly been calculated for microkinetic models based either entirely on DFT or on DFT with the key energies (i.e., those for high-DRC species) being fine-tuned to match experiments. The relationships summarized in this work provide new opportunities for using experiments earlier in the development and optimization of microkinetic models that require input from computational catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Phase Diagrams of Alloys and Their Hydrides via On-Lattice Graph Neural Networks and Limited Training Data

Efficient prediction of sampling-intensive thermodynamic properties is needed to evaluate material performance and permit high-throughput materials modeling for a diverse array of technology applications. To alleviate the prohibitive computational expense of high-throughput configurational sampling with density functional theory (DFT), surrogate modeling strategies like cluster expansion are many orders of magnitude more efficient but can be difficult to construct in systems with high compositional complexity. We therefore employ minimal-complexity graph neural network models that accurately predict and can even extrapolate to out-of-train distribution formation energies of DFT-relaxed structures from an ideal (unrelaxed) crystallographic representation. This enables the large-scale sampling necessary for various thermodynamic property predictions that may otherwise be intractable and can be achieved with small training data sets. Two exemplars, optimizing the thermodynamic stability of low-density high-entropy alloys and modulating the plateau pressure of hydrogen in metal alloys, demonstrate the power of this approach, which can be extended to a variety of materials discovery and modeling problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exploring Structure-Sensitive Relations for Small Species Adsorption Using Machine Learning

Accurate prediction of adsorption energies on heterogeneous catalyst surfaces is crucial to predicting reactivity and screening materials. Adsorption linear scaling relations have been developed extensively but often lack accuracy and apply to one adsorbate and a single binding site type at a time. These facts undermine their ability to predict structure sensitivity and optimal catalyst structure. Using machine learning on nearly 300 density functional theory calculations, we demonstrate that generalized coordination number scaling relations hold well for oxygen- and high-valency carbon-binding species but fail for others. Here we reveal that the valency and the electronic coupling of a species with the surface, along with the site type and its coordination environment, are critical for small species adsorption. The model simultaneously predicts the adsorption energy and preferred site and significantly outperforms linear scalings in accuracy. It can expose the structure sensitivity of chemical reactions and enable enhanced catalyst activity via engineering particle shape and facet defects. The generality of our methodology is validated by training the model with transition metal data and transferring it to predict adsorption energies on single-atom alloys.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanistic insight into the effects of A-site doping in perovskite nanorods on interfacial Li-ion and Na-ion transport in polymer–perovskite composite electrolytes

Composite polymer electrolytes (CPEs) offer a promising pathway to safer and higher energy density electrochemical energy storage yet optimizing interfacial ion transport remains a critical challenge due to high resistances and a lack of understanding of the complex interfacial electrostatic interactions. This work investigates the mechanistic origins of the interfacial ion transport enhancement in single-ion conducting (SIC) polymer electrolytes filled with A-site disordered perovskite LaxMyTiOz (LMTO) (M = Li, Na, K) nanorods. By integrating Density Functional Theory (DFT) with broadband dielectric spectroscopy (BDS) and pulsed-field-gradient NMR (PFG-NMR), we analyze the ion transport behavior in Li-based and Na-based SIC systems containing compositionally distinct LMTO fillers (LMTO800 and LMTO900). Specifically, we introduce the site-to-site energy difference (ΔE) descriptor, which extends beyond single-site descriptors such as adsorption energy and ion binding energies. We demonstrate that ΔE together with ion binding energies provide a more comprehensive picture of the energy landscape at the polymer–ceramic interface. The combined DFT and experimental results show that A-site doping in the ceramic filler can effectively modulate ΔE and ion binding energies, creating a favorable potential energy landscape that promotes enhanced ion transport along the interface. This comprehensive study illustrates the impact of ceramic compositional tuning on interfacial ion transport and provides materials design rules towards real-world applications.

B shepard, Lauren [Pennsylvania State University]↗