SEARCH · Engineering Papers
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.
Choline based eutectic solvents are developed for CO 2 capture. The physical properties as a function of the H-bond donor and acceptor components are examined. Accordingly, molecular descriptors are identified for future property predictions.
QSPR analyses can be used to identify useful descriptors leading to statistical models for membrane crossover. This data-driven approach can be used to evaluate ROMs for asymmetric non-aqueous redox flow batteries.
Although cerium oxide (CeO 2 ) is widely used as a catalyst support, its limited defect sites and surface oxygen vacancy/mobility should be improved. The incorporation of zirconium (Zr) in the cerium (Ce) lattice is shown to increase the number of oxygen vacancies and improve catalytic activity. Using a fixed surface density (SD) of copper (∼2.3 Cu atoms per nm 2 ) as a surface species, the role of the support (Ce y Zr 1−y O 2 (y = 1.0, 0.9, 0.6, 0.5, and 0.0)) and defect site effects in the CO oxidation reaction was investigated. Spectroscopic (e.g., Raman, XRD, XPS) and microscopic (e.g., SEM-EDX, HR-TEM) characterization techniques were applied to evaluate the defect sites, crystallite size, lattice parameters, chemical composition, oxidation states of elements and microstructure of the catalysts. Here, the CO oxidation reaction with varied CO : O 2 ratios (1 : 5, 1 : 1, and 1 : 0.5 (stoichiometric)) was used as a model reaction to describe the relationship between the structure and the catalytic performance of each catalyst. Based on the characterization results of Ce y Zr 1−y O 2 materials, the addition of Zr causes physical and chemical changes to the overall material. The inclusion of Zr into the structure of CeO 2 decreased the overall lattice parameter of the catalyst and increased the number of defect sites. The prepared catalysts were able to reach complete CO conversion (∼100%) at low temperature conditions (<200 °C), each showing varied reaction activity. The difference in CO oxidation activity was then analyzed and related to the structure, wherein Cu loading, surface oxygen vacancies, reduction–oxidation ability, CuO x –support interaction and oxygen mobility in the catalyst were the crucial descriptors.
High entropy oxides (HEOs) have emerged as promising catalysts for several important chemical transformations including alkane activation. Hydrogen adsorption energy (HAE) has been used as a key descriptor for many reactions including methane C–H activation and hydrogen evolution reactions. Hence, understanding the relationship between HAEs and the surface chemistry of HEO surfaces could lay the foundation for meaningful correlations among methane C–H activation, HAE, and the complex, local environment of HEO surfaces. Here, we used a medium-entropy oxide as a prototypical system – Mg 0.25 Ni 0.25 Cu 0.25 Zn 0.25 O with a rock-salt structure – to interrogate these relationships. We sampled 2000 different surfaces of its (100) plane and calculated the HAEs at randomly chosen surface O sites using density functional theory (DFT). Our analysis of the 2000 data points reveals that the HAEs at the surface O sites are significantly influenced by the local environment around the adsorption sites, particularly the nature of the metal atom directly below the surface O site where H adsorbs. After comparing several popular graph-neural-network-based machine learning models, we found that the DimeNet++ model performed best achieving satisfactory accuracy in predicting HAEs for both Mg 0.25 Ni 0.25 Cu 0.25 Zn 0.25 O and slightly varied compositions. Our work underscores the promise of such models and the need for further refinement to address the complexity of HEOs.
Exploring the structural, chemical, and physical properties of matter on the nano- and atomic scales has become possible with the recent advances in aberration-corrected electron energy-loss spectroscopy (EELS) in scanning transmission electron microscopy (STEM). However, the current paradigm of STEM-EELS relies on the classical rectangular grid sampling, in which all surface regions are assumed to be of equal a priori interest. However, this is typically not the case for real-world scenarios, where phenomena of interest are concentrated in a small number of spatial locations, such as interfaces, structural and topological defects, and multi-phase inclusions. One of the foundational problems is the discovery of nanometer- or atomic-scale structures having specific signatures in EELS spectra. Herein, we systematically explore the hyperparameters controlling deep kernel learning (DKL) discovery workflows for STEM-EELS and identify the role of the local structural descriptors and acquisition functions in experiment progression. In agreement with the actual experiment, we observe that for certain parameter combinations the experiment path can be trapped in the local minima. We demonstrate the approaches for monitoring the automated experiment in the real and feature space of the system and knowledge acquisition of the DKL model. Based on these, we construct intervention strategies defining the human-in-the-loop automated experiment (hAE). This approach can be further extended to other techniques including 4D STEM and other forms of spectroscopic imaging. The hAE library is available on Github at https://github.com/utkarshp1161/hAE/tree/main/hAE.
The interfacial hydrostatic stress of the working electrode, not its surface-normal stress, is a more general and widely applicable stress descriptor to use for the mechanics-modified equilibrium potential Δ U .
Binding energies of adsorbed CO and H are key descriptors governing the activity and selectivity of the co-electrolysis of CO 2 and H 2 O to produce syngas with desired CO/H 2 ratios. Palladium hydride (PdH), which forms in situ at negative overpotentials, has been identified as the active Pd phase for CO 2 reduction to syngas. Herein, binding energies of CO and H are determined using temperature programmed desorption (TPD) of CO and H 2 from Pd(111), PdH/Pd(111), and Cu/PdH/Pd(111) under ultra-high vacuum (UHV) conditions. TPD results reveal that desorption of H 2 from subsurface PdH occurs at 460 K, while desorption from surface PdH is more facile at 320 K. CO desorption temperatures shift 20 K lower on PdH/Pd(111) compared to on Pd(111). The presence of 0.7 ML Cu further increases the desorption temperature of H 2 by 30 K while simultaneously reducing CO desorption temperatures by 70 K. Density functional theory (DFT) calculations show that CO adsorption onto Pd sites is hindered on the 0.7 ML Cu/PdH/Pd(111) surface while the kinetic barrier for H 2 desorption is increased. The trends in the binding energies of CO and H on model surfaces are consistent with electrochemical measurements of CuPd powder catalysts in a membrane electrode assembly (MEA), where H 2 evolution is reduced while CO production is enhanced compared to unmodified Pd catalysts. Overall, the results from model surface studies (TPD and DFT) provide a prediction and explanation for the activity and CO/H 2 ratios observed in electrochemical experiments. This study also demonstrates that CuPd is a promising catalyst with reduced Pd-loading to produce CO-rich syngas.
The growing significance of lignin-first biorefineries, which focus on upgrading the aromatics resulting from lignin depolymerization, presents opportunities for bioproduct synthesis using microbial strains capable of funneling a diverse array of phenolics into a single commodity chemical. In this study, we evaluated a biomass-to-bioproduct pipeline involving the reductive catalytic fractionation (RCF) of poplar biomass followed by biological funneling with a Novosphingobium aromaticivorans strain that produces 2-pyrone-4,6-dicarboxylic acid (PDC), a potential bioplastic precursor. Considering the impact of solvent on RCF reactor operating pressure, and the potential inhibitory effects of solvent on downstream microbial funneling, we performed an analysis of six pure solvents, namely methanol, ethanol, isopropanol, isobutanol, 1,4-dioxane and ethylene glycol, and different variations of their aqueous mixtures comprising 5 to 50 vol% water. For each pure solvent and solvent/water system, we measured phenolic monomer yields in the RCF process and PDC yields from the phenolic monomers. We then developed correlation models that relate phenolic monomer yields from RCF-derived samples to Hansen solubility parameters to determine solvent descriptors that contribute to high yields. Furthermore, we developed an integrated biorefinery system to estimate the minimum selling price (MSP) of PDC and the associated carbon footprint to identify solvent systems with better costs and sustainability metrics. These analyses resulted in the 50 vol% methanol/water system being identified as optimal because it reduces RCF reactor pressure and is compatible with microbial funneling with N. aromaticivorans. This solvent system produced 63 g PDC per kg biomass (264 g PDC per kg lignin) from 85 g phenolic monomers per kg biomass at a reduced reactor pressure of 48 bar (reduced by 26% compared to our previous poplar-to-PDC pipeline). The MSP for this system is $\$$13.98 per kg of purified PDC (carbon footprint of 1.47 kg CO 2 e per kg), which is about 24% lower than a previously described poplar-to-PDC pipeline and 46% lower than a lignin-to-PDC pipeline that used pure methanol as the solvent. The results from this study illustrate improvements that can be made in lignocellulosic biorefineries that are compatible with the hybrid chemical and biological processes needed to gain value from lignin.
Due to their high volumetric hydrogen storage capacity under moderate storage conditions, TiFe alloys have been widely investigated as candidates for practical solid-state hydrogen storage. Partially substituting Ti or Fe sites can improve the key characteristics of TiFe alloys, such as the first hydrogen absorption step (activation) and the equilibrium hydrogen pressure (thermodynamic properties). However, the selection of substitution elements has heavily relied on intuition and trial-and-error. Also, conventional substitution strategies have mainly focused on single-element substitution within the TiFe alloy, limiting the design space and tunability for target applications. Here, to address this limitation, we report a multi-element substitution strategy motivated by an efficient, data-driven machine learning (ML) approach combined with corroborating density functional theory (DFT) calculations. Our models successfully predict experimentally measured hydride stability in five selected alloys using only compositional descriptors. Most importantly, the multi-element substitution leads to enhanced activation properties compared to pure TiFe, achieving near room-temperature activation behaviour. This work provides a method for on-demand tuning of hydrogen storage and activation properties, which may have broad implications for data-driven discovery of energy storage materials.
The lithium argyrodite sulfide solid electrolyte Li 6 PS 5 Cl has attracted considerable interest for all-solid-state batteries owing to its high ionic conductivity, which can be further enhanced through ionic substitution. Although a variety of substitutions have been investigated, thioarsenate argyrodites remain comparatively underexplored. Here, we systematically investigate the phase stability and Li-ion conduction mechanisms in superionic Br-incorporated thioarsenate argyrodites using first-principles calculations and molecular dynamics simulations based on machine learning interatomic potentials (MLIPs). Systematic variation of S/Br site inversion reveals that an optimal degree of anion disorder significantly enhances inter-cage connectivity and facilitates long-range Li-ion diffusion. Configurational entropy serves as an effective quantitative descriptor of anion disorder, exhibiting a strong correlation with ionic conductivity. While greater anion disorder induced by site inversion and higher Br content enhances ionic conductivity up to 50 mS cm −1 , it simultaneously reduces structural stability. This trade-off results in an optimal window in which a moderate level of disorder yields conductivities exceeding 20 mS cm −1 while maintaining synthetic feasibility. In conclusion, this work highlights the reliability and efficiency of MLIPs for elucidating ion-transport mechanisms and accelerating the design of novel superionic argyrodites.
Proton exchange membrane water electrolyzers (PEMWEs) are promising zero-emission technologies. However, their high cost remains a barrier to widespread adoption. Iridium oxide is commonly used as an oxygen evolution reaction (OER) catalyst, and its cost and scarcity make it essential to reduce its loading while increasing its activity. Evaluation of iridium oxide activity should be carried out in the membrane electrode assembly (MEA) configuration to replicate realistic operating conditions. Herein, we present a comprehensive benchmarking framework to accurately evaluate the amorphous and crystalline iridium oxides at the MEA level. By systematically varying the catalyst loading, this study confirmed that each MEA was utilized uniformly, presenting intrinsic electrochemical properties independent of the loading. Through intrinsic charge density determined by voltammetry, we established two electrochemical descriptors to evaluate catalyst redox reactions. The mass activity was evaluated by correlating current vs. loading, and the slope provides loading-independent mass activity. The effect of the porous transport layer on OER activity was discussed, identifying a ‘background’ current at zero-loading. In conclusion, this study highlights potential pitfalls in MEA-level catalyst screening and underscores the importance of the loading study for reliable results.
The universal tendency in scanning probe microscopy (SPM) over the last two decades is to transition from simple 2D imaging to complex detection and spectroscopic imaging modes. The emergence of complex SPM engines brings forth the challenge of reliable data interpretation, i.e., conversion from detected signals to descriptors specific to tip–surface interactions and subsequently to material’s properties. In this work, we implemented a Bayesian inference approach for the analysis of the image formation mechanisms in band excitation SPM. Compared to the point estimates in classical functional fit approaches, Bayesian inference allows for the incorporation of extant knowledge of materials and probe behavior in the form of corresponding prior distribution and return the information on the material functionality in the form of readily interpretable posterior distributions. We explore the nonlinear mechanical behaviors spatially in a classical ferroelectric material, PbTiO 3 . We observe the non-trivial evolution of the Duffing stiffness term and the nonlinearity of the sample surface, determine spatial clustering of the nonlinear response, and perform a Landau analysis on predicting the nonlinear coefficient, which indicates that ferroelectric behavior can be a cause of the observed results. These observations suggest that the spectrum of anomalous behaviors at the ferroelectric domain walls may be broader than previously believed and can extend to non-conventional mechanical properties in addition to static and microwave conductance.
We present a machine learning approach for accurately predicting formation energies of binary compounds in the context of crystal structure predictions. The success of any machine learning model depends significantly on the choice of representation used to encode the relevant physical information into machine-learnable data. We test different representation schemes based on partial radial and angular distribution functions (RDF+ADF) on Al–Ni and Cd–Te structures generated using our genetic algorithm for structure prediction. We observe a remarkable improvement in predictive accuracy upon transitioning from global to atom-centered representations, resulting in a threefold decrease in prediction errors. We show that a support vector regression model using a combination of atomic radial and angular distribution functions performs best at the formation energy prediction task, providing small root mean squared errors of 3.9 meV/atom and 10.9 meV/atom for Al–Ni and Cd–Te, respectively. We test the performance of our models against common traditional descriptors and find that RDF- and ADF-based representations significantly outperform many of those in the prediction of formation energies. The high accuracy of predictions makes our machine learning models great candidates for the exploration of energy landscapes.
We present pyflosic, an open-source, general-purpose python implementation of the Fermi–Löwdin orbital self-interaction correction (FLO-SIC), which is based on the python simulation of chemistry framework (pyscf) electronic structure and quantum chemistry code. Thanks to pyscf, pyflosic can be used with any kind of Gaussian-type basis set, various kinds of radial and angular quadrature grids, and all exchange-correlation functionals within the local density approximation, generalized-gradient approximation (GGA), and meta-GGA provided in the libxc and xcfun libraries. A central aspect of FLO-SIC is the Fermi-orbital descriptors, which are used to estimate the self-interaction correction. Importantly, they can be initialized automatically within pyflosic; they can also be optimized within pyflosic with an interface to the atomic simulation environment, a python library that provides a variety of powerful gradient-based algorithms for geometry optimization. Although pyflosic has already facilitated applications of FLO-SIC to chemical studies, it offers an excellent starting point for further developments in FLO-SIC approaches, thanks to its use of a high-level programming language and pronounced modularity.
The plasticity of amorphous solids undergoing shear is characterized by quasi-localized rearrangements of particles. While many models of plasticity exist, the precise relationship between the plastic dynamics and the structure of a particle’s local environment remains an open question. Previously, machine learning was used to identify a structural predictor of rearrangements called “softness.” Although softness has been shown to predict which particles will rearrange with high accuracy, the method can be difficult to implement in experiments where data are limited and the combinations of descriptors it identifies are often difficult to interpret physically. Here, we address both of these weaknesses, presenting two major improvements to the standard softness method. First, we present a natural representation of each particle’s observed mobility, allowing for the use of statistical models that are both simpler and provide greater accuracy in limited datasets. Second, we employ persistent homology as a systematic means of identifying simple, topologically informed, structural quantities that are easy to interpret and measure experimentally. We test our methods on two-dimensional athermal packings of soft spheres under quasi-static shear. We find that the same structural information that predicts small variations in the response is also predictive of where plastic events will localize. We also find that an excellent accuracy is achieved in athermal sheared packings using simply a particle’s species and the number of nearest neighbor contacts.
The Perdew–Zunger self-interaction correction (PZ-SIC) improves the performance of density functional approximations for the properties that involve significant self-interaction error (SIE), as in stretched bond situations, but overcorrects for equilibrium properties where SIE is insignificant. This overcorrection is often reduced by local scaling self-interaction correction (LSIC) of the PZ-SIC to the local spin density approximation (LSDA). Here, we propose a new scaling factor to use in an LSIC-like approach that satisfies an additional important constraint: the correct coefficient of the atomic number Z in the asymptotic expansion of the exchange–correlation (xc) energy for atoms. LSIC and LSIC+ are scaled by functions of the iso-orbital indicator z σ , which distinguishes one-electron regions from many-electron regions. LSIC+ applied to the LSDA works better for many equilibrium properties than LSDA-LSIC and the Perdew, Burke, and Ernzerhof generalized gradient approximation (GGA), and almost close to the strongly constrained and appropriately normed (SCAN) meta-GGA. LSDA-LSIC and LSDA-LSIC+, however, fail to predict interaction energies involving weaker bonds, in sharp contrast to their earlier successes. It is found that more than one set of localized SIC orbitals can yield a nearly degenerate energetic description of the same multiple covalent bond, suggesting that a consistent chemical interpretation of the localized orbitals requires a new way to choose their Fermi orbital descriptors. To make a locally scaled down SIC to functionals beyond the LSDA requires a gauge transformation of the functional’s energy density. Here, the resulting SCAN-sdSIC, evaluated on SCAN-SIC total and localized orbital densities, leads to an acceptable description of many equilibrium properties including the dissociation energies of weak bonds.
A machine learning (ML) methodology that uses a histogram of interaction energies has been applied to predict gas adsorption in metal–organic frameworks (MOFs) using results from atomistic grand canonical Monte Carlo (GCMC) simulations as training and test data. In this work, the method is first extended to binary mixtures of spherical species, in particular, Xe and Kr. In addition, it is shown that single-component adsorption of ethane and propane can be predicted in good agreement with GCMC simulation using a histogram of the adsorption energies felt by a methyl probe in conjunction with the random forest ML method. Here, the results for propane can be improved by including a small number of MOF textural properties as descriptors. We also discuss the most significant features, which provides physical insight into the most beneficial adsorption energy sites for a given application.
In this paper, we investigate the electronic structure of a planar mononuclear Cu-based molecule [Cu(C 6 H 4 S 2 ) 2 ] z in two oxidation states ( z = –2, –1) using density-functional theory (DFT) with Fermi–Löwdin orbital (FLO) self-interaction correction (SIC). The dianionic Cu-based molecule was proposed to be a promising qubit candidate. Self-interaction error within approximate DFT functionals renders severe delocalization of electron and spin densities arising from 3 d orbitals. The FLO-SIC method relies on optimization of Fermi–Löwdin orbital descriptors (FODs) with which localized occupied orbitals are constructed to create SIC potentials. Starting with many initial sets of FODs, we employ a frozen-density loop algorithm within the FLO-SIC method to study the Cu-based molecule. We find that the electronic structure of the molecule remains unchanged despite somewhat different final FOD configurations. In the dianionic state (spin S = 1/2), FLO-SIC spin density originates from the Cu d and S p orbitals with an approximate ratio of 2:1, in quantitative agreement with multireference calculations, while in the case of SIC-free DFT, the orbital ratio is reversed. Overall, FLO-SIC lowers the energies of the occupied orbitals and, in particular, the 3 d orbitals unhybridized with the ligands significantly, which substantially increases the energy gap between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) compared to SIC-free DFT results. The FLO-SIC HOMO–LUMO gap of the dianionic state is larger than that of the monoanionic state, which is consistent with experiment. Our results suggest a positive outlook of the FLO-SIC method in the description of magnetic exchange coupling within 3 d -element-based systems.