Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “crystal structure prediction”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Integrated Computational Materials and Mechanical Modeling for Additive Manufacturing of Alloys with Graded Structure Used in Fossil Fuel Power Plants

Wire-arc additive manufacturing (WAAM) has demonstrated its unique capability of producing large-size alloy components with a significantly reduced fabrication time and enhanced geometry design freedom. In this project, the team has developed an ICME (Integrated Computational Materials Engineering) modeling framework, which supports the WAAM of the AUSC (Advanced Ultra-Supercritical) power plant components. The manufacturing design has been applied to Inconel 740H, steel P91, as well as the dissimilar alloy components between steel P91 and Inconel 740H. The ICME model framework is developed by considering two types of modeling. First, mechanistic modeling has been applied to control the printing quality and understand the sequence of the dissimilar printing of the wall structure. The following models have been included in the developed ICME framework: finite element thermal model, grain structure model, residual stress simulation, crystal plasticity model, CALPHAD-based precipitation kinetic model, phase stability prediction, thermal expansion predictive model, and heuristic creep model. Secondary, a physics-based machine learning model has also been developed based on the ICME model structure. The machine learning model development is based on the ICME model prediction with calibration of the experiments. In addition, the WAAM has been utilized as a high-throughput experimental tool rapidly generating a gradient of alloy composition to facilitate experimental database generation for process-structure-property relationships. Such a database directly supported the ICME-enhanced machine learning, which further assisted in intermediate composition block design between P91 and 740H. A high-throughput screening study of the oxidation resistance has been performed based on such high-throughput experimentation. Based on the computational design, several dissimilar alloy manufacturing with post-heat treatment have been performed with a comprehensive evaluation of mechanical performance, including hardness mapping, yield strength, creep resistance. In this project, the single component of P91 and 740H processed by WAAM after heat treatment designed by ICME has demonstrated higher performance in yield strength and creep resistance than the wrought materials. The P91 sample prepared by WAAM with ICME-designed heat treatment performs better than P92 in creep resistance. The designed graded alloy printing with intermediate block shows a promising performance that exceeds the traditional welding. Moreover, the current research indicates the high need for location-specific design analysis with uncertainty quantification, an important topic that deserves more dedicated research. The achievement of this project demonstrated the promising future of WAAM in structural alloy manufacturing for energy power plant development. Successful printing requires synergetic efforts made by manufacturing, mechanical, and materials sciences.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2024

THD-C Sheet: A Novel Nonbenzenoid Carbon Allotrope with Tetra-, Hexa-, and Dodeca-Membered Rings

Here, we propose a novel two-dimensional carbon-based structure with tetra-, hexa-, and dodeca-membered rings, which we refer to by the abbreviated name, THD-C. The structure presents a mixture of sp–sp 2 hybridization and can potentially be synthesized by the topological assembly of 4-ethynyldiphenylacetylene molecules. By employing first-principles calculations, the stability and ease of synthesis of the sheet are investigated and compared with various C-allotropes. We predict its metallic behavior and excellent kinetic and dynamic stability. Due to the crystal structure of the sheet, a strong mechanical anisotropy is observed. The effects of functionalization on the electronic properties of the material are also studied, and different semiconducting systems are obtained. The potential of THD-C for energy storage in metal-based batteries, hydrogen storage, and catalysis is also investigated, and we find a superior performance in comparison to graphite and other allotropes. The quantum confinement effect is investigated by constructing nanoribbons and nanotubes of various sizes. For ribbons, we find that tailor-made electronic and magnetic properties can be obtained and explored in potential spintronic devices. Additionally, we observe that nanotubes are conducting irrespective of their chirality and can potentially be used for capture, storage, and separation of industrially relevant small gas molecules.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Analytical nonadiabatic coupling and state-specific energy gradient for the crystal field Hamiltonian describing lanthanide single-ion magnets

Paramagnetic molecules with a metal ion as an electron spin center are promising building blocks for molecular qubits and high-density memory arrays. However, fast spin relaxation and decoherence in these molecules lead to a rapid loss of magnetization and quantum information. Nonadiabatic coupling (NAC), closely related to spin-vibrational coupling, is the main source of spin relaxation and decoherence in paramagnetic molecules at higher temperatures. Predicting these couplings using numerical differentiation requires a large number of computationally intensive ab initio or crystal field electronic structure calculations. To reduce computational cost and improve accuracy, we derive and implement analytical NAC and state-specific energy gradient for the ab initio parametrized crystal field Hamiltonian describing single-ion molecular magnets. Our implementation requires only a single crystal field calculation. In addition, the accurate NACs and state-specific energy gradients can be used to model spin relaxation using sophisticated nonadiabatic molecular dynamics, which avoids the harmonic approximation for molecular vibrations. To test our implementation, we calculate the NAC values for three lanthanide complexes. Finally, the predicted values support the relaxation mechanisms reported in previous studies.

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

Evidence of Ba-substitution induced spin-canting in the magnetic Weyl semimetal EuCd 2 As 2

Recently EuCd 2 As 2 was predicted to be a magnetic Weyl semimetal with a lone pair of Weyl nodes generated by A-type antiferromagnetism and protected by a rotational symmetry. However, it was soon discovered that the actual magnetic structure broke the rotational symmetry and internal pressure was later suggested as a route to stabilize the desired magnetic state. In this work we test this prediction by synthesizing a series of Eu 1-x BaxCd 2 As 2 single crystals and studying their structural, magnetic, and transport properties via both experimental techniques and first-principles calculations. We find that small concentrations of Ba (~3%–10%) lead to a small out-of-plane canting of the Eu moment. However, for higher concentrations this effect is suppressed and a nearly in-plane model is recovered. Studying the transport properties we find that all compositions show evidence of an anomalous Hall effect dominated by the intrinsic mechanism as well as large negative magnetoresistances in the longitudinal channel. A nonmonotonic evolution of the transport properties is seen across the series which correlates to the proposed canting suggesting canting may enhance the topological effects. Careful density functional theory calculations using an all-electron approach revise prior predictions finding a purely ferromagnetic ground state with in-plane moments for both the EuCd 2 As 2 and Eu 0.5 Ba 0.5 Cd 2 As 2 compounds, corroborating our experimental findings. This work suggests that Ba substitution can tune the magnetic properties in unexpected ways which correlate to changes in measures of topological properties, encouraging future work to locate the ideal Ba concentration for Eu moment canting.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Simultaneous prediction of structural properties in epitaxially–grown GaN with quantum and conventional multi–output learning algorithms

Hundreds of GaN thin film crystal plasma–assisted molecular beam epitaxy synthesis experiment records spanning two decades were organized into a dataset correlating the growth experiment design parameters with discrete, binary determinations of crystallinity and surface morphology. Conventional data science techniques as well as both quantum and classical multi–output supervised machine learning algorithms were implemented to investigate the relationships between the operating parameter data and the structural figures of merit. Correlation coefficients, decision tree nodes, p–values, and SHAP values all support substrate temperature and gallium effusion cell conditions as being statistically significant for simultaneously influencing GaN crystallinity and surface morphology. Here, a conventional deep neural network learned best from the data, followed by a quantum–classical hybrid gradient boosting algorithm. When combined with calculations of uncertainty intervals based on VennAbers predictors, machine learning predictions of both structural properties show good agreement with results reported in published experimental literature.

36 MATERIALS SCIENCE↗
DOE OSTI2025

High-Entropy Alloys for Accelerator Beam Window Applications

Development of novel high-entropy alloys (HEAs) is currently underway for potential use as beam windows in future multi-megawatt target systems at Fermilab. HEAs encompass a new class of materials with a vast design space allowing for material properties to be tailored for particular applications and to potentially offer improved resistance to beam-induced radiation damage and thermal shock effects. The alloy systems being studied consist of several compositions of AlCoCrMnTiV with 4 6 component elements. These alloys are all predicted by CALPHAD simulation to have a single-phase BCC crystal structure and low density, with some compositions displaying ordered, nanoscale precipitates. This presentation will briefly discuss alloy design and synthesis before giving a detailed description of the characterization studies of these HEAs in both the pristine state and post-irradiation by low-energy heavy ions to high damage levels. Electron microscopy techniques to quantify elemental homogeneity and composition, determine grain size, shape, and orientation, and quantify lattice parameters, defect structures and precipitate phases are all being used to study alloy microstructures. Mechanical properties of the alloys at the microscale will be reported. The evolution of these properties as a function of radiation damage will also be described. Bulk thermal characteristics of these HEAs have been tested to measure specific heat capacity and coefficient of thermal expansion as a function of temperature. To determine bulk tensile properties a miniature tensile testing apparatus is under development; it s commissioning will be covered briefly. The talk will conclude with our plans for alloy down-selection.

Burleigh, A. [Fermilab]↗
DOE OSTI2020

Phase Diagram and Structure Map of Binary Nanoparticle Superlattices from a Lennard-Jones Model

A first-principles prediction of the binary nanoparticle phase diagram assembled by solvent evaporation has eluded theoretical approaches. In this paper, we show that a binary system interacting through the Lennard-Jones (LJ) potential contains all experimental phases in which nanoparticles are effectively described as quasi hard spheres. We report a phase diagram consisting of 53 equilibrium phases, whose stability is quite insensitive to the microscopic details of the potentials, thus giving rise to some type of universality. Furthermore, we show that binary lattices may be understood as consisting of certain particle clusters, i.e., motifs, that provide a generalization of the four conventional Frank–Kasper polyhedral units. Our results show that metastable phases share the very same motifs as equilibrium phases. Finally, we discuss the connection with packing models, phase diagrams with repulsive potentials, and the prediction of likely experimental superlattices.

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

Predicting Elastic Properties of Materials from Electronic Charge Density Using 3D Deep Convolutional Neural Networks

Materials representation plays a key role in machine learning-based prediction of materials properties and new materials discovery. Currently both graph and three-dimensional (3D) voxel representation methods are based on the heterogeneous elements of the crystal structures. Here, we propose to use electronic charge density (ECD) as a generic unified 3D descriptor for materials property prediction with the advantage of possessing close relation with the physical and chemical properties of materials. We developed an ECD-based 3D convolutional neural networks (CNNs) for predicting the elastic properties of materials, in which CNNs can learn effective hierarchical features with multiple convolving and pooling operations. Extensive benchmark experiments over 2170 $Fm\bar3m$ face-centered-cubic materials show that our ECD-based CNNs can achieve good performance for elasticity prediction. Especially, our CNN models based on the fusion of elemental Materials-Agnostic Platform for Informatics and Exploration features and ECD descriptors achieved the best fivefold cross-validation performance. More importantly, we showed that our ECD-based CNN models can achieve significantly better extrapolation performance when evaluated over nonredundant data sets, where there are few neighbor-training samples around test samples. As an additional validation, we evaluated the predictive performance of our models on 329 materials of space group $Fm\bar3m$ by comparing to density functional theory calculated values, which shows a better prediction power of our model for bulk modulus than shear modulus. Because of the unified representation power of ECD, it is expected that our ECD-based CNN approach can also be applied to predict other physical and chemical properties of crystalline materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Data-Driven Kinetic Reaction Networks for Separation Chemistry

Understanding complex, multistep chemical reactions at the molecular level is a major challenge whose solution would greatly benefit the design and optimization of numerous chemical processes. The separation of rare-earth (4f) and actinide (5f) elements is an example where improving our chemical understanding is important for designing and optimizing new chemistries, even with a limited number of observations. Here, in this work, we leverage data-driven artificial intelligence and machine-learning approaches to develop kinetic reaction networks that describe the liquid–liquid extraction mechanism of uranium using N,N-di-2-ethylhexyl-isobutyramide (DEHiBA). Specifically, we compare and contrast the properties of two classes of models: (1) purely data-driven models that are regularized using chemistry-agnostic, L1 regression and (2) chemistry-informed models that are regularized using relative reaction energies provided by quantum mechanical calculations. We observe that purely data-driven models are unbiased, simple, and accurate in their predictions of experimental measurements when provided with sufficient data but are difficult to fully constrain and interpret. In contrast, chemistry-informed models exhibit significantly improved chemical interpretability and consistency, providing a detailed description of the separation process while achieving high accuracy through ensemble averaging. Overall, the dominant species predicted to be extracted into the organic phase is UO 2 (NO 3 ) 2 (DEHiBA) 2 , agreeing with experimental slope analysis, thermodynamic modeling, EXAFS, and crystal structures. This work demonstrates that leveraging the fundamental structure of the problem can lead to efficient learning schemes that provide both accurate predictions and chemical insights at a low computational cost.

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

Ab initio prediction of an order-disorder transition in Mg 2 Ge O 4 : Implication for the nature of super-Earth's mantles

Here we present an ab initio prediction of an order-disorder transition (ODT) from a I¯42d-type to Th 3 P 4 -type phase in the cation sublattices of Mg 2 GeO 4 , a post-post-perovskite phase. This uncommon type of prediction is achieved by carrying out a high-throughput sampling of atomic configurations in a 56-atom supercell followed by a Boltzmann ensemble statistics calculation. Mg 2 GeO 4 is a low-pressure analog of I¯42d-type Mg 2 SiO 4 , a predicted major planet-forming phase of super-Earths' mantles. Therefore, a similar ODT is anticipated in I¯42d-type Mg 2 SiO 4 as well, which should impact the internal structure and dynamics of these planets. Furthermore, the prediction of this Th 3 P 4 -type phase in Mg 2 GeO 4 further enhances the relationship between the crystal structures of Earth/planet-forming silicates and oxides at extreme pressures and those of rare-earth sesquisulfides at low pressures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Rapid and robust antibody Fab fragment crystallization utilizing edge-to-edge beta-sheet packing

Antibody therapeutics are one of the most important classes of drugs. Antibody structures have become an integral part of predicting the behavior of potential therapeutics, either directly or as the basis of modeling. Structures of Fab:antigen complexes have even greater value. While the crystallization and structure determination of Fabs is easy relative to many other protein classes, especially membrane proteins, broad screening and optimization of crystalline hits is still necessary. Through a comprehensive review of rabbit Fab crystal contacts and their incompatibility with human Fabs, we identified a small secondary structural element from the rabbit light chain constant domain potentially responsible for hindering the crystallization of human Fabs. Upon replacing the human kappa constant domain FG loop (HQGLSSP) with the two residue shorter rabbit loop (QGTTS), we dramatically improved the crystallization of human Fabs and Fab:antigen complexes. Our design, which we call “Crystal Kappa”, enables rapid crystallization of human fabs and fab complexes in a broad range of conditions, with less material in smaller screens or from dilute solutions.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2026

Beyond Melting: Amorphous Bonding for Joining and Consolidation

Crystallization may be the hidden constraint in thermoplastic composite manufacturing. It requires tightly controlled cooling, induces residual stresses through shrinkage, and introduces path-dependent behavior that complicates predictive modeling yet remains essential for structural performance. This work asks: can bonding be achieved without relying on melt-driven crystallization? To address this, thin (5–20 μm) polyetherimide (PEI) interlayers are pre-healed to slow-cooled polyaryletherketone (PAEK) in two contexts. The first, Thermabond®, is sub-melt joining of low melt-PAEK laminates. Results show that bond quality is governed primarily by processing (i.e., adequate healing and film handling) rather than modest changes in interlayer thickness. This concept is then extended to laminate-scale manufacturing through an architecture known as OATMEAL (Out-of-autoclave Amorphous/semicrystalline Thermoplastic Material for Energy-efficient Aerospace-grade Laminates). PEI is healed to carbon fiber reinforced polyetheretherketone (PEEK) at the prepreg and excess PEI is then ablated from the surface. Crystallinity is developed off-line during prepreg fabrication, while subsequent consolidation occurs below the melt temperature to preserve it. Cross-ply warpage experiments show that, contrary to intuition, repeated amorphous interfaces reduce global curvature by lowering the effective stress lock-in temperature and eliminating crystallization shrinkage from the lamina response. Correspondingly, laminate behavior is accurately predicted using classical laminate theory (CLT) with a single effective stress-free temperature, whereas conventional CF/PEEK requires accounting for crystallization-driven effects. By decoupling interfacial healing from crystallization, OATMEAL enables sub-melt consolidation, reduces energy consumption by up to 75%, and increases manufacturing throughput by fivefold. These results demonstrate that amorphous bonding is not only a joining strategy, but a pathway to more predictable and scalable thermoplastic composite manufacturing.

solidification↗
DOE OSTI2026

Mechanistic Insights into Defect-Mediated Crystallization Revealed by Lattice Strain Evolution

Structural defects and lattice strain are intrinsic to many crystalline materials, yet their roles in controlling chemical reaction mechanisms and directing crystallization pathways remain poorly understood. Here, in this study, we revealed the three-dimensional evolution of strain and dislocation defects at the nanoscale during the growth of heterogeneously nucleated barite (BaSO 4 ) and calcite (CaCO 3 ) crystals by using coherent X-ray scattering, electron microscopy, and molecular simulations. Unlike barite, which formed with minimal internal strain, calcite developed dislocation defects and exhibited spatially varying strain that increased during growth. During growth in Sr-rich solutions, calcite likely incorporates Sr 2+ into the defects, which further modulates the local lattice structure and increases both the compressive and tensile strain. These findings suggest that calcite crystallization was likely dominated by attachment of precursor phases, which gave rise to defect-enriched domain structures not predicted by classical growth models. By linking defect formation to ion incorporation and growth dynamics, this work provides fundamental insight into how lattice-level strain heterogeneity governs the chemical reactivity of ionic crystals.

Bragg coherent diffractive imaging↗
DOE OSTI2023

Defect graph neural networks for materials discovery in high-temperature clean-energy applications

We present a graph neural network approach that fully automates the prediction of defect formation enthalpies for any crystallographic site from the ideal crystal structure, without the need to create defected atomic structure models as input. Here we used density functional theory reference data for vacancy defects in oxides, to train a defect graph neural network (dGNN) model that replaces the density functional theory supercell relaxations otherwise required for each symmetrically unique crystal site. Interfaced with thermodynamic calculations of reduction entropies and associated free energies, the dGNN model is applied to the screening of oxides in the Materials Project database, connecting the zero-kelvin defect enthalpies to high-temperature process conditions relevant for solar thermochemical hydrogen production and other energy applications. The dGNN approach is applicable to arbitrary structures with an accuracy limited principally by the amount and diversity of the training data, and it is generalizable to other defect types and advanced graph convolution architectures. In conclusion, it will help to tackle future materials discovery problems in clean energy and beyond.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1995

Finding the global minimum: a fuzzy end elimination implementation

The 'fuzzy end elimination theorem' (FEE) is a mathematically proven theorem that identifies rotameric states in proteins which are incompatible with the global minimum energy conformation. While implementing the FEE we noticed two different aspects that directly affected the final results at convergence. First, the identification of a single dead-ending rotameric state can trigger a 'domino effect' that initiates the identification of additional rotameric states which become dead-ending. A recursive check for dead-ending rotameric states is therefore necessary every time a dead-ending rotameric state is identified. It is shown that, if the recursive check is omitted, it is possible to miss the identification of some dead-ending rotameric states causing a premature termination of the elimination process. Second, we examined the effects of removing dead-ending rotameric states from further considerations at different moments of time. Two different methods of rotameric state removal were examined for an order dependence. In one case, each rotamer found to be incompatible with the global minimum energy conformation was removed immediately following its identification. In the other, dead-ending rotamers were marked for deletion but retained during the search, so that they influenced the evaluation of other rotameric states. When the search was completed, all marked rotamers were removed simultaneously. In addition, to expand further the usefulness of the FEE, a novel method is presented that allows for further reduction in the remaining set of conformations at the FEE convergence. In this method, called a tree-based search, each dead-ending pair of rotamers which does not lead to the direct removal of either rotameric state is used to reduce significantly the number of remaining conformations. In the future this method can also be expanded to triplet and quadruplet sets of rotameric states. We tested our implementation of the FEE by exhaustively searching ten protein segments and found that the FEE identified the global minimum every time. For each segment, the global minimum was exhaustively searched in two different environments: (i) the segments were extracted from the protein and exhaustively searched in the absence of the surrounding residues; (ii) the segments were exhaustively searched in the presence of the remaining residues fixed at crystal structure conformations. We also evaluated the performance of the method for accurately predicting side chain conformations. We examined the influence of factors such as type and accuracy of backbone template used, and the restrictions imposed by the choice of potential function, parameterization and rotamer database. Conclusions are drawn on these results and future prospects are given.

NASA Program Exobiology↗
DOE OSTI2021

Predicting Energetics Materials’ Crystalline Density from Chemical Structure by Machine Learning

To expedite new molecular compound development, a long-sought goal within the chemistry community has been to predict molecules’ bulk properties of interest a priori to synthesis from a chemical structure alone. In this work, we demonstrate that machine learning methods can indeed be used to directly learn the relationship between chemical structures and bulk crystalline properties of molecules, even in the absence of any crystal structure information or quantum mechanical calculations. We focus specifically on a class of organic compounds categorized as energetic materials called high explosives (HE) and predicting their crystalline density. An ongoing challenge within the chemistry machine learning community is deciding how best to featurize molecules as inputs into machine learning models—whether expert handcrafted features or learned molecular representations via graph-based neural network models—yield better results and why. We evaluate both types of representations in combination with a number of machine learning models to predict the crystalline densities of HE-like molecules curated from the Cambridge Structural Database, and we report the performance and pros and cons of our methods. Our message passing neural network (MPNN) based models with learned molecular representations generally perform best, outperforming current state-of-the-art methods at predicting crystalline density and performing well even when testing on a data set not representative of the training data. However, these models are traditionally considered black boxes and less easily interpretable. Here, to address this common challenge, we also provide a comparison analysis between our MPNN-based model and models with fixed feature representations that provides insights as to what features are learned by the MPNN to accurately predict density.

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

End-to-end optimization for battery materials and molecules by combining graph neural networks and reinforcement learning

The National Renewable Energy Laboratory (NREL), together with the Colorado School of Mines (CSM) and Colorado State University (CSU), has developed a machine learning-enhanced approach to design new battery materials. Currently, such materials are designed in part via numerous expensive high-fidelity computational simulations that predict the performance of a given composition. Even with computational screening tools, the vast landscape of possible molecular or crystal structures exceeds current and future computational capacity. Improving the efficiency by which new materials can be optimized will therefore disrupt the cost, risk, and time required to bring new energy solutions to the marketplace. Predicting the properties of an organic molecule or periodic crystalline material given its structure has grown increasingly common. These approaches leverage large-scale computational and experimental databases and ML approaches such as graph neural networks. The inverse design problem of finding a material that possesses desired properties is substantially more challenging, since enumerating all valid material structures is not feasible. In this project, we leveraged recent success in reinforcement learning to efficiently navigate this high-dimensional search space. Just as algorithms can find the optimal chess moves from nearly limitless options, we train an approach to evolve a simple starting structure into a complex structure that possess the desired properties. Our solution has been demonstrated by applying it to two related design application tasks for short- and long-term energy storage, respectively: (1) the design of solid-state ion conductors and (2) the design of organic redox-active materials. The project has resulted an open-source software library for material design, documented examples of applying the library to both organic and inorganic material optimization, and peer-reviewed publications detailing the data, computational models, and resulting candidate materials.

25 ENERGY STORAGE↗
DOE OSTI2022

End-to-End Optimization for Battery Materials and Molecules by Combining Graph Neural Networks and Reinforcement Learning

The National Renewable Energy Laboratory (NREL), together with the Colorado School of Mines (CSM) and Colorado State University (CSU), has developed a machine learning-enhanced approach to the design of new battery materials. Currently, such materials are designed in part via numerous expensive high-fidelity computational simulations that predict the performance of a given composition. Even with computational screening tools, the vast landscape of possible molecular or crystal structures exceeds current and future computational capacity. Improving the efficiency by which new materials can be optimized will therefore disrupt the cost, risk, and time required to bring new energy solutions to the marketplace. Predicting the properties of an organic molecule or periodic crystalline material given its structure has grown increasingly common. These approaches leverage large-scale computational and experimental databases and ML approaches such as graph neural networks. The inverse design problem of finding a material that possesses desired properties is substantially more challenging, since enumerating all valid material structures is not feasible. In this project, we leveraged recent success in reinforcement learning to efficiently navigate this high-dimensional search space. Just as algorithms can find the optimal chess moves from nearly limitless options, we train an approach to evolve a simple starting structure into a complex structure that possess the desired properties. Our solution has been demonstrated by applying it to two related design application tasks for short- and long-term energy storage, respectively: (1) the design of solid-state ion conductors and (2) the design of organic redox-active materials. The project has resulted an open-source software library for material design, documented examples of applying the library to both organic and inorganic material optimization, and peer-reviewed publications detailing the data, computational models, and resulting candidate materials.

25 ENERGY STORAGE↗