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 469 records · Page 26

Prediction of Structures and Atomization Energies of Coinage Metals, (M)n, n < 20: Extrapolation of Normalized Clustering Energies to Predict the Cohesive Energy

The geometries of the group 11 coinage metals (n = 2-20) were optimized to determine the lowest energy isomers for each cluster size, singlets for even numbers and doublets for odd numbers. For copper and silver, 2-D (planar) geometries were favored up to n = 6. For gold, 2D (planar) geometries were favored up to n = 13. Normalized clustering energies were plotted as a function of cluster size (n-1/3, for n = 4-20) with various DFT functionals and the CCSD(T)-F12b method and were extrapolated to predict the bulk cohesive energy. In the case of copper and silver, there is excellent agreement between the cohesive energies predicted at the CCSD(T)-F12b level of theory and the experimental values. For gold, the CCSD(T)-F12b values needed to be corrected for spin-orbit relativistic effects to obtain good agreement with experiment. Electronic properties including the HOMO-LUMO gaps for the even clusters and the spin densities for the odd clusters were calculated. The lowest gap is predicted to occur for n = 16 where the HOMO and LUMO are very similar in shape.

Persuad, Rudradatt R.↗

Exploiting electricity market dynamics using flexible electrolysis units for retrofitting methanol synthesis

Here we investigate the economic viability of integrating flexible electrolysis units to produce hydrogen in methanol synthesis processes. Specifically, we investigate whether this approach can help reduce methanol production costs by strategically exploiting dynamics of electricity markets. Our study integrates high-fidelity process simulations, optimization tools, and microkinetic modeling (informed by density functional theory) to conduct detailed techno-economic analyses and to compare performance against traditional processes that use hydrogen produced via steam-methane reforming (SMR). We also use this approach to estimate the levelized cost of hydrogen (LCOH) as a function of time-varying electricity prices (from day-ahead and real-time prices) and of key techno-economic parameters. Our results show that the proposed electrification framework is cost-competitive under certain electricity market conditions. Specifically, we find that, when the electrolysis system is operated in flexible mode (and can respond to dynamics of electricity markets), the associated electricity cost nearly collapses to zero. Conversely, when the unit is not flexible (and cannot respond to markets), the electricity cost comprises 60% of the total cost. Our results also reveal that the LCOH of the flexible electrolysis system participating in real-time electricity markets is 31% lower than the LCOH obtained from SMR. Overall, this indicates that exploiting the dynamics of electricity markets can make hydrogen production cost-competitive and this can lead to viable alternatives to electrify methanol production and other hydrogen-based processes.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Localizing Tungsten Single Atoms around Tungsten Nitride Nanoparticles for Efficient Oxygen Reduction Electrocatalysis in Metal–Air Batteries

Combining isolated atomic active sites with those in nanoparticles for synergizing complex multistep catalysis is being actively pursued in the design of new electrocatalyst systems. However, these novel systems have been rarely studied due to the challenges with synthesis and analysis. Herein, a synergistically catalytic performance is demonstrated with a 0.89 V (vs reversible hydrogen electrode) onset potential in the four-step oxygen reduction reaction (ORR) by localizing tungsten single atoms around tungsten nitride nanoparticles confined into nitrogen-doped carbon (W SAs/WNNC). Through density functional theory calculations, it is shown that each of the active centers in the synergistic entity feature a specific potential-determining step in their respective reaction pathway that can be merged to optimize the intermediate steps involving scaling relations on individual active centers. Impressively, the W SAs/WNNC as the air cathode in all-solid-state Zn-air and Al-air batteries demonstrate competitive durability and reversibility, despite the acknowledged low activity of W-based catalyst toward the ORR.

25 ENERGY STORAGE↗

Background optimization of powder electron diffraction for implementation of the e-PDF technique and study of the local structure of iron oxide nanocrystals

The local structural characterization of iron oxide nanoparticles is explored using a total scattering analysis method known as pair distribution function (PDF) (also known as reduced density function) analysis. The PDF profiles are derived from background-corrected powder electron diffraction patterns (the e-PDF technique). Due to the strong Coulombic interaction between the electron beam and the sample, electron diffraction generally leads to multiple scattering, causing redistribution of intensities towards higher scattering angles and an increased background in the diffraction profile. In addition to this, the electron–specimen interaction gives rise to an undesirable inelastic scattering signal that contributes primarily to the background. The present work demonstrates the efficacy of a pre-treatment of the underlying complex background function, which is a combination of both incoherent multiple and inelastic scatterings that cannot be identical for different electron beam energies. Therefore, two different background subtraction approaches are proposed for the electron diffraction patterns acquired at 80 kV and 300 kV beam energies. From the least-square refinement (small-box modelling), both approaches are found to be very promising, leading to a successful implementation of the e-PDF technique to study the local structure of the considered nanomaterial.

36 MATERIALS SCIENCE↗

Achieving High‐Performance and 2000 h Stability in Anion Exchange Membrane Fuel Cells by Manipulating Ionomer Properties and Electrode Optimization

The primary function of the ionomers that are incorporated into fuel cell electrode catalyst layers is to provide pathways for ion transport between the catalyst active sites and the electrolyte. This is influenced by many variables, including the ion-exchange capacity, water uptake, and molecular weight. In anion exchange membrane fuel cells (AEMFCs), controlling ionomer water uptake is particularly important and tailoring this property in each electrode is an important consideration when looking to maximize cell performance. In this study, three poly(norbornene) tetrablock copolymer ionomers with a range of physical properties are synthesized and incorporated into AEMFC anode and cathode electrodes. Systematic electrode engineering with these ionomers allows the peak power density to be increased by 100% (1.6 W cm -2 → 3.2 W cm -2 ) and the current density at 0.2 V to be increased by 59% (5.9 A cm -2 → 9.4 A cm -2 ). Moreover, the top-performing electrode configuration is tested in an operating AEMFC at the US Department of Energy defined current density of 600 mA cm -2 for 2000 h, showing a record-low voltage decay rate of 15.36 µV h -1 – only 3.65% –a over 2000 h. Here, this work sets a new bar for AEMFCs, reporting the best combination of performance and durability of any AEMFC to date.

08 HYDROGEN↗

Computer-aided design of stability enhanced nicotinamide cofactor biomimetics for cell-free biocatalysis

Cell-free biocatalysis (CFB) is an efficient and environmentally friendly method to synthesize molecules such as pharmaceuticals, biochemicals, and biofuels through the in vitro use of enzyme cascades. These enzymes often require redox cofactors to drive chemical reactions. Natural redox cofactors (NAD(P)H) are expensive to isolate, motivating synthetic nicotinamide cofactor biomimetics (NCBs) as a cost-effective solution. A select handful of NCBs have been identified as potential NAD(P)H alternatives with comparable or improved redox capabilities, however, they display a tendency to degrade in common buffers. In this study, a library of 132 NCB candidates is systematically generated, over 85% of which have not been characterized in the literature, to expand the diversity of currently explored NCBs. The decomposition mechanism of NCBs in phosphate is evaluated using density functional theory (DFT), revealing protonation at the nicotinamide C5 position as a reporter of cofactor stability. Based on this result, we trained a linear regression model on DFT calculated descriptors to predict NCB stability in phosphate buffer, achieving mean absolute error (MAE) and root mean squared error (RMSE) values within computational accuracy. Analysis of key atomic descriptors and qualitative trends in our dataset informed the design of novel NCB candidates we propose with optimized stability. This work enables researchers to predict the relative stability of NCBs before synthesis, thereby streamlining the process to make CFB more affordable and viable at industry scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evaluation of WRF simulation of deep convection in the US Southern Great Plains

The Southern Great Plains (SGP) exhibits a relatively high frequency of periods with extremely high rainfall rates (RR) and hail. Here, seven months of 2017 are simulated using the Weather Research and Forecasting (WRF) model applied at convection permitting resolution with the Mibrandt-Yau microphysics scheme. Simulation fidelity is evaluated, particularly during intense convective events, using data from ASOS stations, dual-polarization RADAR, gridded data sets and observations at the DoE Atmospheric Radiation Measurement site. The spatial gradients and temporal variability of precipitation and the cumulative density functions for both RR and wind speeds exhibit fidelity. Odds ratios >1 indicate WRF is also skillful in simulating high composite reflectivity (cREF, used as a measure of widespread convection) and RR > 5 mmhr –1 over the domain. Detailed analyses of the ten days with highest spatial coverage of cREF >30 dBZ show spatially similar reflectivity fields and high RR in both RADAR data and WRF simulations. However, during periods of high reflectivity, WRF exhibits a positive bias in terms of very high RR (> 25 mmhr –1 ) and hail occurrence, and during the summer and transition months, maximum hail size is underestimated. For some renewable energy applications fidelity is required with respect to the joint probabilities of wind speed and RR and/or hail. While partial fidelity is achieved for the marginal probabilities, performance during events of critical importance to these energy applications is currently not sufficient. Further research into optimal WRF configurations in support of potential damage quantification for these applications is warranted.

54 ENVIRONMENTAL SCIENCES↗

Structure of chalcogen overlayers on Au(111): Density functional theory and lattice-gas modeling

Ordering of different chalcogens, S, Se, and Te, on Au(111) exhibit broad similarities but also some distinct features, which must reflect subtle differences in relative values of the long-range pair and many-body lateral interactions between adatoms. We develop lattice-gas (LG) models within a cluster expansion framework, which includes about 50 interaction parameters. These LG models are developed based on density functional theory (DFT) analysis of the energetics of key adlayer configurations in combination with the Monte Carlo (MC) simulation of the LG models to identify statistically relevant adlayer motifs, i.e., model development is based entirely on theoretical considerations. The MC simulation guides additional DFT analysis and iterative model refinement. Given their complexity, development of optimal models is also aided by strategies from supervised machine learning. The model for S successfully captures ordering motifs over a broader range of coverage than achieved by previous models, and models for Se and Te capture the features of ordering, which are distinct from those for S. More specifically, the modeling for all three chalcogens successfully explains the linear adatom rows (also subtle differences between them) observed at low coverages of ~0.1 monolayer. The model for S also leads to a new possible explanation for the experimentally observed phase with a (5 × 5)-type low energy electron diffraction (LEED) pattern at 0.28 ML and to predictions for LEED patterns that would be observed with Se and Te at this coverage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Etched Silicon Planar CRL Optics for the High-Energy X-ray Diffraction Beamlines 11-ID-B and 11-ID-C at the APS

We report the focusing of high-energy X-rays is most commonly achieved by commercially available compound refractive lenses (CRL), which are added together in larger structures. There is a need for high-energy X-ray beamlines to enable focusing in only one dimension down to a few micrometers to study layered materials and thin films. High-energy X-rays above 50 keV require often 100 lenses or more, therefore it is more efficient to etch the complete arrays of lenses into a monolithic substrate like a silicon wafer. Custom lens arrays were fabricated for the high-energy x-ray beamlines 11-ID-B, and 11-ID-C at the Argonne National Laboratory. The fabrication required tuning parameters of the Bosch process to optimize the etching of deep vertical sidewalls at a reasonable speed without overheating. The beamline 11-ID-B is equipped with lens arrays of 85-110 lenses for 59 keV and 190-250 lenses for 87 keV photon energy, focusing the beam vertically below 5.0 µm at the sample positions between 1600 mm to 2100 mm. 11-ID-C has received arrays of 435 lenses that focus 106 keV X-rays to a vertical height of 2.0 µm at the sample position. In both cases, the flux density gain in the focal spot is larger than a factor of 20, relative to using slits to cut the beam size. The lens setup is permanently installed and can be moved in and out of the beam during the experiment without further alignment. The lenses have enabled operando diffraction studies on interfaces in batteries and the development of pair distribution function (PDF) analysis using grazing incidence geometry.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Manufacturing of Fabric Electrodes using a High-Throughput Screening Platform for Redox Flow Batteries

The objective of this project is to establish a new manufacturing methodology with machine learning- based high-throughput screening for the design and development of hierarchical structured, high-performance fabric electrodes for redox flow batteries (RFBs). The end goal of the project is to design and manufacture fabric electrodes for RFB applications that can provide 250 mA/cm2 current density operation for 100-cycles with 80% average energy efficiency. This was accomplished by first examining the structure-performance-property linkages of the electrodes provided by our partner, AvCarb. The electrodes’ microstructure was characterized by determining their pore size distribution, tortuosity, specific surface area, and porosity. The ohmic, charge transfer and mass transfer resistances were then calculated using electrochemical impedance spectroscopy. Carbon cloth electrodes showed the greatest resistance, which was dominated by charge transfer resistance, which we believe is related to the surface functionalization. Full cell cycling was used in order to determine the area specific resistance and energy efficiency of the cells. All of this experimental data and the results of the mathematical model (to increase the amount of inputs with parametric sweeping) were used to develop a machine learning-based model for the design of high-performance fabric electrodes. Using the results from the machine learning tool, optimized electrodes were fabricated by AvCarb. The ohmic, charge transfer and mass transfer resistances for these new electrodes were measured, and both performed better than any of the initial samples which had been provided by AvCarb.

25 ENERGY STORAGE↗

Resolving the Solvation Structure and Transport Properties of Aqueous Zinc Electrolytes from Salt-in-Water to Water-in-Salt Using Neural Network Potential

Zn Cl 2 solutions are promising electrolytes for aqueous zinc-ion batteries. Here, we report a joint computational and experimental study of the structural and dynamic properties of aqueous Zn Cl 2 electrolytes with concentrations ranging from salt-in-water to water-in-salt (WIS). By developing a neural network potential (NNP) model, we perform molecular dynamics (MD) simulations with accuracy but at much larger lengths and longer timescales. The NNP predicted structures are validated by the structure factors measured by X-ray total scattering experiments. The MD trajectories provide a comprehensive and quantitative picture of the Zn 2 + solvation shell structures. Additionally, we find that the O − H covalent bonds in water are strengthened with increasing salt concentration, thus expanding the electrochemical stability window of aqueous electrolytes. In terms of dynamic properties, the calculated and experimentally measured conductivities are in good agreement. Through the analysis of the calculated cation transference number, we propose a three-stage charge carrier transport mechanism with increasing concentration: independent ion transport, strongly correlated ion transport, and small positive charge carrier diffusion through negatively charged polymeric clusters. Our study provides fundamental atomic scale insights into the structure and transport properties of the Zn Cl 2 electrolyte that can aid the optimization and development of WIS electrolytes. Published by the American Physical Society 2025

25 ENERGY STORAGE↗

Chemical Modifications of Ag Catalyst Surfaces with Imidazolium Ionomers Modulate H 2 Evolution Rates during Electrochemical CO 2 Reduction

Bridging polymer design with catalyst surface science is a promising direction for tuning and optimizing electrochemical reactors that could impact long-term goals in energy and sustainability. Particularly, the interaction between inorganic catalyst surfaces and organic-based ionomers provides an avenue to both steer reaction selectivity and promote activity. In this work, we studied the role of imidazolium-based ionomers for electrocatalytic CO 2 reduction to CO (CO 2 R) on Ag surfaces and found that they produce no effect on CO 2 R activity yet strongly promote the competing hydrogen evolution reaction (HER). By examining the dependence of HER and CO 2 R rates on concentrations of CO 2 and HCO 3 –, we developed a kinetic model that attributes HER promotion to intrinsic promotion of HCO 3 – reduction by imidazolium ionomers. We also show that varying the ionomer structure by changing substituents on the imidazolium ring modulates the HER promotion. This ionomer-structure dependence was analyzed via Taft steric parameters and density functional theory calculations, which suggest that steric bulk from functionalities on the imidazolium ring reduces access of the ionomer to both HCO 3 – and the Ag surface, thus limiting the promotional effect. Our results help develop design rules for ionomer–catalyst interactions in CO 2 R and motivate further work into precisely uncovering the interplay between primary and secondary coordination in determining electrocatalytic behavior.

36 MATERIALS SCIENCE↗

Stability and Performance of 3d Transition Metal Carbo‐Sulfides: A Density Functional Theory Exploration for Li‐Ion Battery Anodes

As the demand for high-performance and reliable energy storage devices continues to rise, identifying new anode materials is crucial for advancing Li-ion battery (LIB) technology. Inspired by recent experimental breakthroughs in synthesizing two-dimensional transition metal carbo-chalcogenides (2D-TMCCs), density functional theory calculations are performed to systematically explore their sulfide variants (TM 2 S 2 C) spanning all 3d transition metals in three possible phases. Through comprehensive evaluations of thermodynamic, dynamic, mechanical, and thermal stabilities, seven stable 2D-TMCC candidates are identified, four of which exhibit superior battery performance. Notably, V-based 2D-TMCCs across all three phases deliver moderate open-circuit voltages (OCV), efficient Li diffusion, and substantial capacities, making them promising candidates for industrial applications without requiring specific phase controls. A Cr-based 2D-TMCC (with sulfur atoms above carbon atoms) offers the highest capacity of 515.40 mAh g −1 , the lowest Li diffusion barrier, and an optimal OCV, highlighting its appealing potential as an anode material for LIBs. Furthermore, significant Li–Li spacing and pronounced electron delocalization in these four 2D-TMCCs suggest a reduced risk of dendrite formation. This work expands the 2D-TMCC family and identifies up-and-coming candidates for next-generation LIB anodes.

anode materials↗

Embedded random phase approximation for magnetic systems: H 2 dissociative adsorption on Fe(110)

The random phase approximation (RPA), a method for treating electron correlation, has been shown to be superior to standard density functional theory (DFT) approximations in numerous cases. However, the RPA’s computational cost is substantially higher than that of DFT, particularly restricting its application to extended surfaces. The recently introduced embedded RPA (emb-RPA) approach [Wei et al., J. Chem. Phys. 159(19), 194108 (2023)] reduces this computational cost by approximately two orders of magnitude. While previous applications of emb-RPA focused on non-spin-polarized systems, here we extend the approach to ferromagnetic ones. Unlike other embedded correlated wavefunction methods, such as embedded complete active space self-consistent field theory, emb-RPA is advantageous for spin-polarized systems because the RPA is compatible with unrestricted DFT solutions, which are eigenfunctions of the spin angular momentum operator S z but not the total spin-squared operator S 2 . By applying emb-RPA with specific magnetization constraints, we achieved a speedup of two to three orders of magnitude (one order when accounting for the one-time embedding potential optimization cost) with only small errors (∼50 meV) compared to full periodic RPA. Moreover, emb-RPA significantly reduces the over-binding errors of DFT approximations. In conclusion, we anticipate that the acceleration enabled by the spin-polarized emb-RPA approach will broaden the applicability of RPA to magnetic materials.

Density functional theory↗

Fermi-Löwdin-orbital self-interaction correction using the optimized-effective-potential method within the Krieger-Li-Iafrate approximation

Perdew-Zunger self-interaction correction (PZ-SIC) offers a route to remove self interaction errors on an orbital-by-orbital basis. A recent formulation of PZ-SIC by Pederson, Ruzsinszky and Perdew proposes restricting the unitary transformation to localized orbitals called Fermi-Löwdin orbitals. This formulation, called the FLOSIC method, simplifies PZ-SIC calculations and was implemented self-consistently using a Jacobi-like (FLOSIC-Jacobi) iteration scheme. In this work we implement the FLOSIC approach using the Krieger-Li-Iafrate (KLI) approximation to the optimized effective potential (OEP). We compare the results of present FLOSIC-KLI approach with FLOSIC-Jacobi scheme for atomic energies, atomization energies, ionization energies, barrier heights, polarizability of chains of hydrogen molecules etc. to validate the FLOSIC-KLI approach. The FLOSICKLI approach, which is within the realm of Kohn-Sham theory, predicts smaller energy gaps between frontier orbitals due to the lowering of eigenvalues of the lowest unoccupied orbitals. Results show that atomic energies, atomization energies, ionization energy as an absolute of highest occupied orbital eigenvalue, and polarizability of chains of hydrogen molecules between the two methods agree within 2%. Lastly the FLOSIC-KLI approach is used to determine the vertical ionization energies of water clusters.

74 ATOMIC AND MOLECULAR PHYSICS↗

Co-optimization of fuel properties, combustion system geometry, and injection strategy for conventional diesel fuel

Here, studies have shown that fuel properties can impact an engine’s operation in several ways, including ignition delay, sooting tendency, mixture formation, and combustion temperature. In mixing-controlled compression ignition (MCCI) engines, the fuel system design and piston bowl geometry significantly affect combustion performance and emissions. Based on current information, it is difficult to draw conclusions about fuel property effects and sensitivities. The central fuel hypothesis approach used in the US Department of Energy Co-Optima program has worked well for spark ignition fuels: identifying critical fuel property ranges is sufficient to screen fuel blends that are expected to maximize efficiency and reduce pollutant emissions. However, for MCCI-relevant fuels, the information gained from past studies is not sufficient to build such a merit function or to allow for performing a similar screening of fuel blends. It is hypothesized that a co-optimization of a fuel’s physical and chemical properties, combustion system geometry, and injection strategy could leverage synergies between the effects of the fuel properties and geometries, resulting in improved performance over state-of-the-art. A machine learning–assisted unconstrained global optimization algorithm was used to explore a design space comprising 23 independent variables. The results show that physical property effects were minimal even for large variations in fuel properties, and the only interaction effect that was observed was the effect of varied fuel density parameters on fuel/air mixture formation. Nevertheless, these interactions were not sufficient in magnitude to significantly affect optimization results. Therefore, analysis of the results suggests that fuel physical properties cannot be leveraged in a co-optimization context to increase engine efficiency.

33 ADVANCED PROPULSION SYSTEMS↗

3D printing of architected sulfur cathodes with dual-site atomic catalysts for accelerated polysulfide kinetics and Li-ion transport in high areal-loading lithium–sulfur batteries

The practical deployment of lithium–sulfur batteries (LSBs) is hindered by fundamental limitations in conventional slurry-cast cathodes, including poor sulfur utilization, sluggish ion transport, and low areal capacity, particularly in thick electrodes required for high energy density. To address these challenges, we present direct ink writing (DIW) as an additive manufacturing strategy to fabricate advanced current-collector-free, 3D-printed sulfur cathodes (3DP S@CoNi-DSACs/NC) with hierarchically porous architectures that enhance lithium-ion diffusion, promote electrolyte penetration, and reduce interfacial resistance. The synergistic effects of Co/Ni dual-atom sites accelerate redox kinetics and mitigate polysulfide shuttling. As a result, the optimized 3DP cathode with a sulfur loading of 5.4 mg cm −2 demonstrated excellent rate capability, delivering a high reversible capacity of 1041.4 mAh g −1 at 1C with 85.5% capacity retention after 1000 cycles, significantly outperforming its cast counterpart. Remarkably, even at a higher sulfur loading of 8.1 mg cm −2 , the 3DP cathode maintains outstanding performance, achieving a discharge capacity of 1538.4 mAh g −1 and an areal capacity of 12.5 mAh cm −2 at 0.1C. This study not only demonstrates the functional integration of catalytically active materials into 3D printable sulfur cathode architectures but also offers a scalable and transformative platform for building high-performance LSBs beyond conventional electrode manufacturing methods.

3D electrode↗

Global study of separable pairing interaction in covariant density functional theory

A systematic global investigation of pairing properties based on all available experimental data on pairing indicators has been performed for the first time in the framework of covariant density functional theory. It is based on the separable pairing interaction. The optimization of the scaling factors of this interaction to experimental data clearly reveals its isospin dependence in the neutron subsystem. However, the situation is less certain in the proton subsystem since similar accuracy of the description of pairing indicators can be achieved both with isospin-dependent and mass-dependent scaling factors. The differences in the functional dependencies of scaling factors lead to uncertainties in the prediction of proton and neutron pairing properties which are especially pronounced at high isospin and could have a significant impact on some physical observables. For a given part of the nuclear chart the scaling factors for spherical nuclei are smaller than those for deformed ones; this feature exists also in nonrelativistic density functional theories. Furthermore, its origin is traced back to particle-vibration coupling in odd-A nuclei which is missing in all existing global studies of pairing. Although the present investigation is based on the NL5(E) covariant energy density functional (CEDF), its general conclusions are expected to be valid also for other CEDFs built at the Hartree level.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗