Search NASA⌕ Search

SEARCH · Search NASA

Results for “High-throughput calculations”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Conditional diffusion machine-learning framework for mapping valence electron distribution from convergent beam electron diffraction

Quantitative convergent beam electron diffraction (CBED) enables determination of aspherical valence electron distributions through refinement of low-order structure factors, which are highly sensitive to chemical bonding and charge density variations. However, conventional quantitative CBED (QCBED) requires solving a highly nonlinear inverse problem with many coupled parameters, and computationally intensive dynamical diffraction calculations, making it time-consuming and difficult to apply to complex systems. More broadly, reconstructing charge density and orbital electron distribution from diffraction data has long been a central challenge in both x-ray and electron crystallography. Here, in this study, we introduce an artificial-intelligence (AI)-based framework that replaces traditional refinement with a data-driven inverse solver. Using a large synthetic CBED dataset generated by Bloch-wave simulations, we train a conditional diffusion model to directly infer crystal structural parameters and multipole density formalism parameters, and hence valence electron distributions, from CBED patterns alone. By learning from forward simulations across realistic parameter space, the model effectively solves the inverse problem. Compared with direct regression approaches, the diffusion-based framework provides posterior parameter distributions for rigorous uncertainty quantification while preserving quantitative fidelity and reducing analysis time by orders of magnitude. By eliminating the need for external single-crystal x-ray diffraction data and complex nonlinear refinement, this approach enables practical, high-throughput, and in situ quantitative CBED, enabling real-time mapping of valence electron distributions and their correlation with functional responses in quantum and energy materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Extending SLUSCHI for Automated Diffusion Calculations

We present an extension of the SLUSCHI package (Solid and Liquid in Ultra Small Coexistence with Hovering Interfaces) to enable automated diffusion calculations from first-principles molecular dynamics. While the original SLUSCHI workflow was designed for melting temperature estimation via solid-liquid coexistence, we adapt its input and output handling to isolate the volume search stage and generate one production trajectory suitable for diffusion analysis. Post-processing tools parse VASP outputs, compute mean-square displacements (MSD), and extract tracer diffusivities using the Einstein relation with robust error estimates through block averaging. Diagnostic plots, including MSD curves, running slopes, and velocity autocorrelations, are produced automatically to help identify diffusive regimes. The method has been validated through representative case studies: self-diffusion in Al-Cu liquid alloys, sublattice melting in Li7La3Zr2O12 and Er2O3, interstitial oxygen transport in bcc and fcc Fe, and oxygen diffusivity in Fe-O liquids with variable Si and Al contents. Viscosity and diffusivity are linked through the Stokes-Einstein relation, with composition dependence assessed via simple linear mixing. This capability broadens SLUSCHI from melting-point predictions to transport property evaluation, enabling high-throughput, fully first-principles datasets of diffusion coefficients and viscosities across metals and oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Particle signal considerations for isotope ratio analysis with single particle multi-collector inductively coupled plasma mass spectrometry

Particle analysis has benefitted from the advent of single particle inductively coupled plasma-mass spectrometry (spICP-MS) due to its robustness, sensitivity, and high-throughput nature. Previous methods of spICP-MS have typically utilized quadrupole or time-of-flight mass analyzers and therefore employ electron multiplier-based detectors (such as secondary electron multipliers or microchannel plates). However, to obtain precise measurements on elemental or isotopic ratios within individual particles, multi-collector ICP-MS (MC-ICP-MS) can be used. Here, we investigate Ce isotope ratios, specifically 142 Ce/ 140 Ce, by spMC-ICP-MS using an all-Faraday cup collector array. Using 1 μm (diameter) cerium dioxide particles, integration times of the Faraday cup detectors were varied from 50–500 ms. The signal from the cerium isotopes in the particles was used to determine isotope ratios, which closely matched the expected natural isotopic abundances. Due to the signal decay response from the Faraday cups, the signal from particles lasts much longer than the expected 1–2 ms (up to 100 s of ms). To explore this effect on isotope ratio analysis, multiple ratio analysis methods were used to determine how to obtain optimal precision and accuracy. Relative differences were around 2% for methods that calculated isotope ratios from summing the total signal of an individual particle before calculating the ratio (rather than using every data point individually). It was found that summing all data points per particle, or integrating under the signal peak, yielded both accurate and precise isotope ratios within the particle population. Particles were also sampled off a solid substrate via microextraction, and isotope ratios were determined with relative differences of 0.13% to 9%. This demonstrates the ability to use spMC-ICP-MS to obtain isotope ratios on particles, with little to no relative difference in comparison to the expected ratio, even when operating Faraday detectors at fast 50 ms integration times.

Szakas, Sarah E. [Oak Ridge National Laboratory (O↗
DOE OSTI2025

Design of lightweight BCC multi-principal element alloys with enhanced hydrogen storage using a machine learning-driven genetic algorithm

Body-centered cubic (BCC) based multi-principal element alloy (MPEA) hydrides have demonstrated significant potential for compact and efficient hydrogen storage. In this work, we first leverage machine learning (ML) models to predict the hydrogen affinity, storage capacity and phase stability of BCC MPEAs, creating a unique hydrogen-to-metal (H/M) predictor for materials with unprecedented performance. We developed a metaheuristic optimizer high-throughput framework by interfacing ML models with a genetic algorithm for the accelerated search of {Mg, Al, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Nb, Mo} based lightweight BCC MPEAs with improved hydrogen storage characteristics. We report five new MPEAs with a predicted gravimetric hydrogen storage capacity of around 3.5 wt% or more, including Cr 0.09 Mg 0.73 Ti 0.18 (4.25 wt% H) and Cr 0.21 Nb 0.11 Ti 0.35 V 0.33 (3.5 wt% H). The electronic structure of the top-performing composition, Cr 0.09 Mg 0.73 Ti 0.18 , was analyzed using density functional theory (DFT) to understand the reasons for its improved hydrogen storage properties compared to TiFe (1.90 wt% H), LaNi 5 (1.37 wt% H) or BCC MPEAs like TiVNbCr (3.70 wt% H). Temperature-dependent molecular dynamics (MD) studies were further performed on optimized BCC MPEAs to qualitatively study hydrogen mobility and analyze the effect of different elemental composition on bulk hydrogen diffusion. Our findings demonstrate how a ML assisted genetic algorithm framework can be used for efficient search of stable, lightweight and cost-effective MPEAs while minimizing the need for expensive ab initio calculations.

DFT↗
DOE OSTI2026

Accurate Electron-Phonon Interactions from Advanced Density Functional Theory

Electron-phonon coupling (EPC) is key to understanding charge transport, band renormalization, and superconductivity in energy materials, including correlated transition-metal oxides, ferroelectric perovskites, optoelectronic semiconductors, and phonon-mediated superconductors. Although first-principles density functional theory (DFT)-based EPC calculations are used widely, their predictive power is limited by the accuracy, transferability, and efficiency of the underlying exchange-correlation functionals. These limitations become exacerbated in complex 𝑑- and 𝑓-electron materials, where beyond-DFT approaches and additional corrections, such as the Hubbard 𝑈, are commonly invoked. Here, using the examples of Co⁢O and Ni⁢O, we show how the 𝑟 2 SCAN density functional correctly captures strong EPC effects in transition-metal oxides without requiring the introduction of Hubbard 𝑈 parameters. We also find that 𝑟 2 SCAN successfully describes the subtle interplay between ionic and covalent bonding, and strong EPC effects in the low-temperature rhombohedral phase of the prototypical ferroelectric Ba⁢Ti⁢O 3 , without requiring Hubbard 𝑈 and intersite 𝑉 corrections. We further demonstrate the ability of 𝑟 2 SCAN to accurately model the EPC of the main-group semiconductor Ga⁢As and the phonon-mediated superconducting properties of Mg⁢B 2 , with reliable electronic bands and phonons. Our study establishes 𝑟 2 SCAN-based EPC as a transferable, parameter-free, and computationally efficient framework for predictive material-specific modeling of EPC in energy materials, and opens a practical route toward high-throughput screening of superconductors, thermoelectrics, optoelectronic semiconductors, and oxide electrodes and catalysts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Discovery of tunable and soluble organic emitters for solid-state lasers with a self-driving laboratory

We have recently demonstrated the ability of using self-driving laboratories for AI-driven searches of organic emitters for solid-state lasing devices. Our past workflow featured solubility challenges for such large molecular moieties. In this next-generation study, we return to the drawing board to explore a family of com-pounds that are much solution processable and composed of a set of electronic cores that provide a broader color response. Out of 252 potential candidates,and with guidance from DFT calculations, we selectively perform a compre-hensive study exploring 51 fluorene-based A-B-A type organic laser oligomers, armed with our self-driving lab. The candidates range from simple hydrocarbon molecules to complex hetero atom-mixed molecules. As a result of this study, we highlight diketopyrrolopyrrole and benzodiazole derivatives for their largely red-shifted emissions. Furthermore, we investigate the effect of color change aris-ing from hetero atom permutation, fluorine addition, thiophene coupling, and a combination of fluorine addition and thiophene coupling. Amplified spontaneous emission (ASE) measurements in the solid state further corroborate the lasing potential of selected candidates, reinforcing their suitability for future device applications. The computational study with density functional theory confirms the experimental results.

fluorescence↗
DOE OSTI2025

Rapid Bayesian High Entropy Alloy Designs Fabricated via Wire Arc Additive Manufacturing

Purpose: This project seeks to demonstrate a new high-throughput (rapid) alloy design technique applied to creating new high entropy alloys (HEAs) for extreme environments. High entropy alloys shift the design paradigm from being focused on a single principal element (e.g. nickel-based alloys) to target alloys that include high atomic fractions (X >10%) of multiple elements. These HEA materials can exhibit sluggish diffusion and enhanced corrosion resistance, ideal for potential applications in advanced ultra supercritical (A-USC) steam cycles for power generation. Scope: The addition of multiple elements in high atomic fractions creates an enormous design space that cannot easily be investigated by traditional material design strategies such as designed of experiments (DOE). This project utilizes a Bayesian machine learning algorithm that has been modified to work with calculation of phase diagrams (CALPHAD) software. This Bayesian algorithm reduces manual inputs and increase the likelihood of achieving an optimal solution. Compositional inputs to this algorithm will be assessed using existing material property models for high temperature strength and corrosion resistance. The target for alloy performance will be a 15% (~100 ⁰C) increase in allowable service temperature beyond heat-resistant stainless steels while maintaining or improving alloy cost and corrosion resistance. Haynes 230 was selected as a baseline, which is 57 wt% Ni with 22 wt% Cr 14 wt% W, and 2 wt% Mo as solid solution strengtheners. In addition to rapid design via Bayesian machine learning, the alloys were rapidly fabricated using a multi-wire arc additive manufacturing (mWAAM) technique which allows for precise control of alloy composition and assessing of alloy design “windows” to study composition effects. Build speeds for wire-arc additive processes are among the highest for additive technologies enabling rapid and reliable sample fabrication when compared to conventional methods such as arc button melting. The mWAAM samples will be rapidly characterized via instrumented indentation for room temperature modulus and strength and for elevated temperature strength via hot hardness tests. After being screened with hardness testing, potential alloys will be further evaluated with conventional microscopy techniques including scanning electron microscopy (SEM) and transmission electron microscopy (TEM) to assess agreement with modeling results. The most promising compositions will also be evaluated by printing full sized tensile specimens for mechanical behavior tests at elevated temperatures. Results: Bayesian machine learning of a single performance function was initially used to optimize five performance metrics: 1) single phase stability, 2) yield strength, 3) creep resistance (low diffusion coefficient), 4) freezing range (weldability), and 5) material cost. The single performance function was suboptimal as assumptions had to be made about the results while formulating the optimization. A goal-oriented Bayesian optimization strategy (Hanaoka, 2021) was implemented with CALPHAD for use with the five metrics above. This multi-objective Bayesian optimization (MOBO) enabled the design of NiCrCoFe alloys with V and W additions. A base composition of NiCoCr was selected as Ni provides a stable FCC matrix, Cr aids corrosion/oxidation resistance, and Co is a solid-solutions strengthener that also improves creep by increasing the activation energy. Fe helps reduce diffusion coefficients and cost. Finally, V and W were selected for their reasonable solubility and high atomic misfit to aid in solid solution strengthening. Cracking of the mWAAM specimens was an early issue, and the Easton solidification cracking model (Easton et al., 2014a) was selected for addition to the MOBO function. High performing alloys fabricated by mWAAM included Ni 28 Cr 25 Co 26 Fe 15 V 8 and Ni 62 Cr 18 Co 1 Fe 3 W 15 . It was observed that even after adapting the mWAAM process for W, the W did not fully dissolve. To fully evaluate the Ni 62 Cr 18 Co 1 Fe 3 W 15 composition, a cored wire (80-20 NiCr sheath/powder core) was manufactured and printed via WAAM, and HIP’ing was utilized to homogenize and densify the printed alloy. The V and W alloys produced met metrics 1 (solid solution), 4 (solidification cracking), and 5 (cost). However, an unmodeled mechanism of thermal stress cracking was identified in the WAAM produced materials, perhaps exacerbated by the lack of grain boundary strengthening elements (B, C). Conclusions & Recommendations: A high-throughput (rapid) alloy design technique was applied to designing and manufacturing new high entropy alloys (HEAs) for extreme environments utilizing MOBO and mWAAM. The developed process was rapid and effective in addressing the mechanisms included in the model. The lack of grain boundary strengthening element additions (e.g., B, C) was a simplification that likely produced thermal stress cracking that turned into a large part of the investigation. Additions on the order of 0.005 wt% B and 0.05 wt% C likely would have minimized thermal stress grain boundary cracking. Overall, the high throughput design strategy is promising for rapid design of metrics-driven alloys for advanced ultra supercritical (A-USC) steam cycles for power generation. The MOBO and mWAAM process could be commercialized to accelerate metrics-driven alloy design. In addition, the cored-wire process utilized for scale-up is a promising high-volume process for WAAM alloy development and scale-up.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ab initio property predictions of quinary solid solutions using small binary cells

The Set of Small Ordered Structures (SSOS) approach is an ab initio technique for modelling random solid solutions in which many small structures are averaged so that their correlation functions match those of a desired composition. SSOS has been shown to be effective in reducing the cost of density functional theory calculations relative to other well-known techniques such as cluster expansions and special quasirandom structures for modelling solid solutions. Here in this work, we demonstrate that SSOS’s can be constructed using cells with only a subset of elements while still accurately modelling multi-component systems. Specifically, we show that small binary cells can effectively model two quinary high entropy alloys – NbTaTiHfZr and MoNbTaVW – accurately capturing properties such as formation energy, lattice parameters, elastic constants, and root-mean-square atomic displacements. Overall, this insight is useful for those looking to construct databases of such small structures for predicting the properties of multi-component solid solutions, as it greatly decreases the number of structures that needs to be considered.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Interpretable, extensible linear and symbolic regression models for charge density prediction using a hierarchy of many-body correlation descriptors

Here, density functional theory (DFT) is routinely used to make electronic structure predictions for high-throughput screening of materials and molecules for technologically relevant areas, like the identification of better catalysts, electronic materials, and drug discovery. However, the DFT formalism is limited by (a) its poor (quadratic-to-quartic) scaling, and (b) the need to perform repeated eigenvalue computations of the electronic Hamiltonian as part of its self-consistent field (SCF) iteration procedure to obtain the converged ground state electron density, ρ (r). Approaches that directly predict ρ (r) of a structure with high accuracy can accelerate conventional SCF calculations and can also be used in linearly scaling methods such as orbital-free DFT. To this end, we present a procedure to predict the ground state electron density of molecular and periodic three-dimensional systems directly from the atomic structure with a particular emphasis on physical interpretability. In our framework, ρ (r) is modeled using many-body correlation descriptors that accurately capture the effects of local atomic arrangements in the neighborhood of a grid point. Our use of a linear regression scheme to fit to charge density data enables transparent analysis of the relative contributions of various types of local atomic correlations. By systematically including increasingly complex correlations, our model is shown to accurately predict ρ (r) for a variety of chemically and electronically diverse systems — amorphous Ge, Al(001) slab, crystalline Ga 2 O 3 , molecular benzene, and polyethylene. We then demonstrate a symbolic regression-based protocol to construct easily computable, interpretable features from lower-order correlations that significantly improves our electron density predictions with effectively no increase in the computational cost.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

A Six-Membered Concerted Mechanism for CO 2 Capture by Amines Studied under Charged Microdroplet Reaction Conditions

Carbon dioxide (CO 2 ) capture and storage represents an important technological challenge. A mechanistic understanding of interactions involved in the capture process is necessary not only for technological development but also for efficient conversion of captured CO 2 into value-added materials. Herein, we present a novel contained secondary electrospray ionization platform for studying the interactions of gaseous amines and CO 2 gas under microdroplet reaction conditions, which enables mass spectrometry (MS) characterization of CO 2 capture products and intermediates in real time. We detected [2 M + CO 2 + H] + species, which corresponds to a six-membered intermediate. DFT calculations confirmed the high stability of the protonated six-membered ring intermediate. This finding provides a plausible concerted mechanism in the microdroplet environment that excludes the involvement of thermodynamically disfavored ionic species. The carbamic acid counterpart of the final product/salt was readily characterized by tandem MS. The carbamic acid/amine salt was also isolated and characterized by Fourier transform infrared spectroscopy. By virtue of the fact that headspace vapors of amines are sampled, we were able to establish a high-throughput platform that enabled the CO 2 capture capacity of five different amines to be studied in under 2 min. The same device also enabled the absolute quantification of capture capacity.

Amines↗
DOE OSTI2026

Itinerant Magnetism in Hydride-Synthesized CaCo 12 B 6

A new compound in the underexplored Ca–Co–B phase space has been discovered, validating high-throughput computations from the Open Quantum Materials Database, which predicted thermodynamic stability for CaCo 12 B 6 in the SrNi 12 B 6 structure type. The synthetic effects of different boron precursors and the advantages of using CaH 2 instead of Ca metal were demonstrated by the short synthesis duration and high purity of CaCo 12 B 6 , in contrast with traditional synthesis routes. Powder X-ray diffraction (PXRD) confirmed that CaCo 12 B 6 shares the SrNi 12 B 6 structure ( R $\bar{3}$m (#166), a = 9.469(4) Å, c = 7.468(2) Å, Z = 3) and is water- and air-stable. High-temperature in situ PXRD indicates that CaCo 12 B 6 is stable below 1050 K under vacuum in a silica capillary. CaCo 12 B 6 decomposes between 693 and 773 K during spark-plasma sintering. Density functional theory calculations indicate that CaCo 12 B 6 is metallic with a ferromagnetic ground state. X-ray absorption near-edge spectroscopy and Bader charge analysis indicate that Co atoms in CaCo 12 B 6 lack ionic character. Magnetometry reveals room-temperature paramagnetism with μ eff = 1.7(1)μ B per Co atom and a Weiss constant of +190(10)K. Ferromagnetic ordering occurs below 172(1)K, resulting in a saturation moment of 0.46 μ B per Co atom. Our findings demonstrate that the hydride route is a viable strategy for discovery of new ternary alkaline-earth-transition metal borides analogous to rare-earth-containing counterparts.

diffraction↗
DOE OSTI2024

Millimeter-Scale Exfoliation of hBN with Tunable Flake Thickness for Scalable Encapsulation

As a two-dimensional (2D) dielectric material, hexagonal boron nitride (hBN) is in high demand for applications in photonics, nonlinear optics, and nanoelectronics. Unfortunately, the high-throughput preparation of macroscopic-scale, high-quality hBN flakes with controlled thickness is an ongoing challenge, limiting device fabrication and technological integration. In this report we present a metal thin-film exfoliation method to prepare hBN flakes with millimeter-scale dimension, near-unity yields, and tunable flake thickness distribution from 1-7 layers, a substantial improvement over scotch tape exfoliation. The single crystallinity and high quality of the exfoliated hBN are demonstrated with optical microscopy, atomic force microscopy, Raman spectroscopy, and second harmonic generation. We further explore a possible mechanism for the effectiveness and selectivity based on thin-film residual stress measurements, density functional theory calculations, and transmission electron microscopy imaging of the deposited metal films. We find that the magnitude of the residual tensile stress induced by thin film deposition plays a key role in determining exfoliated flake thickness in a manner which closely resembles 3D semiconductor spalling. Lastly, we demonstrate that our exfoliated, large-area hBN flakes can be readily incorporated as encapsulating layers for other 2D monolayers. Altogether, this method brings us one step closer to the high throughput, mass production of hBN-based 2D photonic, optoelectronic, and quantum devices.

2D materials↗
DOE OSTI2025

Data-Driven Discovery of Bimetallic Nanoparticles Catalysts for the Hydrogenolysis of Polyethylene

Supported platinum nanoparticles are known to convert polyolefins to high-quality liquid hydrocarbons with hydrogen under relatively mild conditions. However, no systematic study has been undertaken using bimetallic catalysts for polyethylene upcycling. Specifically, a total of 98 monometallic and bimetallic combinations (Ag, Cr, Co, Cu, Fe, Ga, In, Mn, Ni, Pd, Pt, Rh, Ru, Zr) on alumina were synthesized utilizing surface organometallic chemistry (SOMC) technique via robotic platform. These were investigated at a small scale (10 mg of catalyst and 50 mg of polyethylene) for their activity for the hydrogenolysis of polyethylene in a high-throughput batch reactor. Combinations of Ni and Co were selected as candidates with high activity toward conversion into paraffin oils. Reaction conditions were optimized with Ni/Co/Al 2 O 3 catalyst at a larger scale (300 mg catalyst and 3 g polyethylene) to obtain a high yield (93.1%) of paraffin wax with desired properties (M n = 380 Da) and low polydispersity (Đ = 1.2). Ni/Co/Al 2 O 3 was compared against Co/Ni/Al 2 O 3 to understand the role of the deposition sequence. When Co is deposited before Ni, a layer of cobalt aluminate is formed upon reduction, stabilizing the deposition of 5 nm metallic Ni particles. When nickel is deposited before Co, particles are larger (average >20 nm) and more oxidized (Ni δ+ in NiAl 2 O 4 ), decreasing the availability of the catalytically active metallic Ni. In conclusion, the difference in electronic environments was also described by DFT calculations, which revealed that smaller 3D clusters of Ni are preferred on CoAl2O4 over the 3D clusters on NiAl 2 O 4 and that these smaller clusters are more reducible, as confirmed experimentally.

Polymer↗
DOE OSTI2024

High-Throughput Screening and Accurate Prediction of Ionic Liquid Viscosities Using Interpretable Machine Learning

Ionic liquids (ILs) are a novel group of green solvents with great promise for various industrial applications, including carbon capture and lignocellulosic biomass deconstruction. However, the use of ILs at the industrial scale remains challenging due to their high viscosities at ambient temperatures. To develop ILs with lower viscosities, a systematic study of their quantitative structure–property relationship (QSPR) is desirable. Here, we developed four machine learning (ML) models to predict viscosity at various temperature and pressure ranges, trained over a wide range of ILs consisting of various cationic and anionic families. ML methods including two-factor polynomial regression (two-factor PR), support vector regression (SVR), feed-forward neural networks (FFNN), and categorical boosting (CATBoost) were developed based on features that have proven useful in previous ML studies: COSMO-RS (conductor-like screening model for real solvents)-derived surface screening charge densities (sigma profiles). FFNN and CATBoost were the most accurate in predicting IL viscosities with lower average absolute relative deviation and higher R2 values on the test set. Tanimoto similarity scores were calculated to characterize the chemical space and structural similarity of the investigated ions. Furthermore, SHapley Additive exPlanation (SHAP) analysis was employed to interpret the ML results. Temperature, the polar area of ILs, and the nonpolar regions of ions are key features that influence the viscosity predictions. Importantly, the IL viscosity prediction here is the most accurate reported to date.

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

Cryptate binding energies towards high throughput chelator design: metadynamics ensembles with cluster–continuum solvation

A tiered forcefield/semiempirical/meta-GGA pipeline together with a thermodynamic scheme designed with error cancellation in mind was developed to calculate binding energies of [2.2.2] cryptate complexes of mono- and divalent cations. Stable complexes of Na, K, Rb, Ca, Zn and Pb were generated, revealing consistent cation–N lengths but highly variable cation–O lengths and an amine stacking mechanism potentially augmenting the cation size selectivity. Metadynamics, used for searching the high-dimensional potential energy surface, together with a cluster–continuum model for affordable – yet accurate – solvation modeling, enabled the discovery of more stable geometries than those previously reported. Similar solvation energy curve shapes for lone vs. coordinated ions enabled rapid solvation convergence via the cancellation of errors stemming from finite cluster sizes. In conclusion, an R 2 of 0.850 vs. experimental aqueous binding energies was obtained, validating this scheme as the backbone of a high-throughput workflow for chelator design.

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

PyLRO: A Python calculator for analyzing long-range structural order

We present PyLRO, an open-source Python calculator designed to detect, quantify, and display long-range order in periodic structures. The program’s design methodology, workflow, and approach to order quantification are described and demonstrated using a simple toy model. Additionally, we apply PyLRO to a series of metastable AlPO 4 structural intermediates from a prior high-pressure study, demonstrating how to compute and visualize structural order in all directions on a Miller sphere. We further highlight the program’s capabilities through a high-throughput analysis of structural patterns in the pressure-induced amorphization of AlPO 4 , revealing atomistic insights into specific energy regions of massive amorphous structures. These results suggest that PyLRO can be a valuable tool for investigating crystal–amorphous transition in materials research.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Wide-ranging predictions of new stable compounds powered by recommendation engines

The computational search for new stable inorganic compounds is faster than ever, thanks to high-throughput density functional theory (DFT). However, stable compound searches remain highly expensive because of the enormous search space and the cost of DFT calculations. To aid these searches, recommendation engines have been developed. We conduct a systematic comparison of the performance of previously developed recommendation engines, specifically ones based on elemental substitution, data mining, and neural network prediction of formation enthalpy. After identifying ways to improve the recommendation engines, we find the neural network to be superior at recommending stable Heusler compounds. Armed with improved recommendation engines, we identify tens of thousands of compounds that are stable at zero temperature and pressure, now available in the Open Quantum Materials Database. We summarize this diverse pool of compounds, including the elusive mixed anion compounds, and two of their many applications: thermoelectricity and solar thermochemical fuel production.

Science & Technology - Other Topics↗
DOE OSTI2024

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↗