Search

In this work, we developed a compositionally transferable machine learning interatomic potential using atomic cluster expansion potential and PBE-D3 method for (NaCl)1-x(MgCl2)x molten salt and we showed that it is possible to fit a robust potential for this pseudo-binary system by only including data from x={0, 1/3, 2/3, 1}. We also assessed the performance of several DFT methods including PBE-D3, PBE-D4, R2SCAN-D4, and R2SCAN-rVV10 on unary NaCl and MgCl2 salts. Our results show that the R2SCAN-D4 method calculates the thermophysical properties of NaCl and MgCl2 with an overall modestly better accuracy compared to the other three methods., Comment: 40 Pages, 9 Figures

Materials design often relies on human-generated hypotheses, a process inherently limited by cognitive constraints such as knowledge gaps and limited ability to integrate and extract knowledge implications, particularly when multidisciplinary expertise is required. This work demonstrates that large language models (LLMs), coupled with prompt engineering, can effectively generate non-trivial materials hypotheses by integrating scientific principles from diverse sources without explicit design guidance by human experts. These include design ideas for high-entropy alloys with superior cryogenic properties and halide solid electrolytes with enhanced ionic conductivity and formability. These design ideas have been experimentally validated in high-impact publications in 2023 not available in the LLM training data, demonstrating the LLM's ability to generate highly valuable and realizable innovative ideas not established in the literature. Our approach primarily leverages materials system charts encoding processing-structure-property relationships, enabling more effective data integration by condensing key information from numerous papers, and evaluation and categorization of numerous hypotheses for human cognition, both through the LLM. This LLM-driven approach opens the door to new avenues of artificial intelligence-driven materials discovery by accelerating design, democratizing innovation, and expanding capabilities beyond the designer's direct knowledge.

We demonstrate the ability of large language models (LLMs) to perform material and molecular property regression tasks, a significant deviation from the conventional LLM use case. We benchmark the Large Language Model Meta AI (LLaMA) 3 on several molecular properties in the QM9 dataset and 24 materials properties. Only composition-based input strings are used as the model input and we fine tune on only the generative loss. We broadly find that LLaMA 3, when fine-tuned using the SMILES representation of molecules, provides useful regression results which can rival standard materials property prediction models like random forest or fully connected neural networks on the QM9 dataset. Not surprisingly, LLaMA 3 errors are 5-10x higher than those of the state-of-the-art models that were trained using far more granular representation of molecules (e.g., atom types and their coordinates) for the same task. Interestingly, LLaMA 3 provides improved predictions compared to GPT-3.5 and GPT-4o. This work highlights the versatility of LLMs, suggesting that LLM-like generative models can potentially transcend their traditional applications to tackle complex physical phenomena, thus paving the way for future research and applications in chemistry, materials science and other scientific domains.

The integration of machine learning (ML) models enhances the efficiency, affordability, and reliability of feature detection in microscopy, yet their development and applicability are hindered by the dependency on scarce and often flawed manually labeled datasets and a lack of domain awareness. We addressed these challenges by creating a physics-based synthetic image and data generator, resulting in a machine learning model that achieves comparable precision (0.86), recall (0.63), F1 scores (0.71), and engineering property predictions (R2=0.82) to a model trained on human-labeled data. We enhanced both models by using feature prediction confidence scores to derive an image-wide confidence metric, enabling simple thresholding to eliminate ambiguous and out-of-domain images resulting in performance boosts of 5-30% with a filtering-out rate of 25%. Our study demonstrates that synthetic data can eliminate human reliance in ML and provides a means for domain awareness in cases where many feature detections per image are needed.

One compelling vision of the future of materials discovery and design involves the use of machine learning (ML) models to predict materials properties and then rapidly find materials tailored for specific applications. However, realizing this vision requires both providing detailed uncertainty quantification (model prediction errors and domain of applicability) and making models readily usable. At present, it is common practice in the community to assess ML model performance only in terms of prediction accuracy (e.g., mean absolute error), while neglecting detailed uncertainty quantification and robust model accessibility and usability. Here, we demonstrate a practical method for realizing both uncertainty and accessibility features with a large set of models. We develop random forest ML models for 33 materials properties spanning an array of data sources (computational and experimental) and property types (electrical, mechanical, thermodynamic, etc.). All models have calibrated ensemble error bars to quantify prediction uncertainty and domain of applicability guidance enabled by kernel-density-estimate-based feature distance measures. All data and models are publicly hosted on the Garden-AI infrastructure, which provides an easy-to-use, persistent interface for model dissemination that permits models to be invoked with only a few lines of Python code. We demonstrate the power of this approach by using our models to conduct a fully ML-based materials discovery exercise to search for new stable, highly active perovskite oxide catalyst materials.

Oxygen ion conductors are crucial for enhancing the efficiency of various clean energy technologies, including fuel cells, batteries, electrolyzers, membranes, sensors, and more. In this study, LaBi2O4Cl is identified as an ultra-fast oxygen conductor from the MBi2O4X (M=rare-earth element, X=halogen element) family, discovered by a structure-similarity analysis of >60k oxygen-containing compounds. Ab initio studies reveal that LaBi2O4Cl has an ultra-low migration barrier of 0.1 eV for oxygen vacancy, significantly lower than 0.6-0.8 eV for interstitial oxygen. Frenkel pairs are the dominant defects in intrinsic LaBi2O4Cl, facilitating notable oxygen diffusion primarily through vacancies at higher temperatures. LaBi2O4Cl with extrinsic oxygen vacancies (2.8%) exhibits a conductivity of 0.3 S/cm at 25{\deg}C, maintains a 0.1 eV diffusion barrier up to 1100{\deg}C, and transitions from extrinsic to mixed extrinsic and intrinsic behavior as the Frenkel pair concentration increases at higher temperatures. Experimental results on synthesized LaBi2O4Cl and Sr-doped LaBi2O4Cl demonstrate comparable or higher oxygen conductivity than YSZ and LSGM below 400 {\deg}C, with lower activation energies. Further experimental optimization of LaBi2O4Cl, including aliovalent doping and microstructure refinement, could significantly enhance its performance and efficiency, facilitating fast oxygen conduction approaching room temperature.

Knowledge of the domain of applicability of a machine learning model is essential to ensuring accurate and reliable model predictions. In this work, we develop a new approach of assessing model domain and demonstrate that our approach provides accurate and meaningful designation of in-domain versus out-of-domain when applied across multiple model types and material property data sets. Our approach assesses the distance between a test and training data point in feature space by using kernel density estimation and shows that this distance provides an effective tool for domain determination. We show that chemical groups considered unrelated based on established chemical knowledge exhibit significant dissimilarities by our measure. We also show that high measures of dissimilarity are associated with poor model performance (i.e., high residual magnitudes) and poor estimates of model uncertainty (i.e., unreliable uncertainty estimation). Automated tools are provided to enable researchers to establish acceptable dissimilarity thresholds to identify whether new predictions of their own machine learning models are in-domain versus out-of-domain., Comment: 43 pages, 34 figures, journal submission

In nuclear and/or solar applications that involve molten salts, impurities frequently enter the salt as either fission products or via corrosion. Impurities can interact and make complexes, but the impact of such complexation on the properties of the salts and corrosion rates has not been understood. Common impurities in molten salts, such as FLiBe, include Cr, Ni, and Fe. Here, we investigate the complexation of Cr and Ni in FLiBe using molecular dynamics based on a machine learning interatomic potential (MLIP) fitted using the atomic cluster expansion (ACE) method. The MLIP allows us to overcome the challenges of simultaneously needing accurate energetics and long time scale to study complexation. We demonstrate that impurity behavior is more difficult to capture than that of concentrated elements with MLIPs due to less sampling in training data, but that this can be overcome by using active learning strategies to obtain a robust fit. Our findings suggest that there is a weak but potentially significant binding free energy between CrF2 and NiF2 in eutectic FLiBe of -0.112 eV. Under certain conditions this binding creates a significant short-range order between the species and lowers the redox potential of NiF2 in the presence of CrF2 in FLiBe, making Ni dissolution more favorable in the presence of Cr as compared to its dissolution in pure FLiBe. However, we find little impact of this complexation on the diffusivity of Ni and Cr. Overall, the methodology presented here suggests an approach to modeling complexation with MLIPs and suggests that interactions between dissolved cations could be playing a significant role in some salt thermophysical properties.

Ensemble models can be used to estimate prediction uncertainties in machine learning models. However, an ensemble of N models is approximately N times more computationally demanding compared to a single model when it is used for inference. In this work, we explore fitting a single model to predicted ensemble error bar data, which allows us to estimate uncertainties without the need for a full ensemble. Our approach is based on three models: Model A for predictive accuracy, Model $A_{E}$ for traditional ensemble-based error bar prediction, and Model B, fit to data from Model $A_{E}$, to be used for predicting the values of $A_{E}$ but with only one model evaluation. Model B leverages synthetic data augmentation to estimate error bars efficiently. This approach offers a highly flexible method of uncertainty quantification that can approximate that of ensemble methods but only requires a single extra model evaluation over Model A during inference. We assess this approach on a set of problems in materials science., Comment: 14 pages, 4 figures, 1 table

Perovskite SrVO$_3$ has recently been proposed as a novel electron emission cathode material. Density functional theory (DFT) calculations suggest multiple low work function surfaces and recent experimental efforts have consistently demonstrated effective work functions of ~2.7 eV for polycrystalline samples, both results suggesting, but not directly confirming, some fraction of even lower work function surface is present. In this work, thermionic electron emission microscopy (ThEEM) and high-field ultraviolet photoemission spectroscopy are used to study the local work function distribution and measure the work function of a partially-oriented-(110)-SrVO$_3$ perovskite oxide cathode surface. Our results show direct evidence of low work function patches of about 2.1 eV on the cathode surface, with corresponding onset of observable thermionic emission at 750 $^o$C. We hypothesize that, in our ThEEM experiments, the high applied electric field suppresses the patch field effect, enabling the direct measurement of local work functions. This measured work function of 2.1 eV is comparable to the previous DFT-calculated work function value of the SrVO-terminated (110) SrVO$_3$ surface (2.3 eV) and SrO terminated (100) surface (1.9 eV). The measured 2.1 eV value is also much lower than the work function for the (001) LaB$_6$ single crystal cathode (~2.7 eV) and comparable to the effective work function of B-type dispenser cathodes (~2.1 eV). If SrVO$_3$ thermionic emitters can be engineered to access domains of this low 2.1 eV work function, they have potential to significantly improve thermionic emitter-based technologies.

Discovering new materials that efficiently catalyze the oxygen reduction and evolution reactions is critical for facilitating the widespread adoption of solid oxide fuel cell and electrolyzer (SOFC/SOEC) technologies. Here, we develop machine learning (ML) models to predict perovskite catalytic properties critical for SOFC/SOEC applications, including oxygen surface exchange, oxygen diffusivity, and area specific resistance (ASR). The models are based on trivial-to-calculate elemental features and are more accurate and dramatically faster than the best models based on ab initio-derived features, potentially eliminating the need for ab initio calculations in descriptor-based screening. Our model of ASR enables temperature-dependent predictions, has well calibrated uncertainty estimates and online accessibility. Use of temporal cross-validation reveals our model to be effective at discovering new promising materials prior to their initial discovery, demonstrating our model can make meaningful predictions. Using the SHapley Additive ExPlanations (SHAP) approach, we provide detailed discussion of different approaches of model featurization for ML property prediction. Finally, we use our model to screen more than 19 million perovskites to develop a list of promising cheap, earth-abundant, stable, and high performing materials, and find some top materials contain mixtures of less-explored elements (e.g., K, Bi, Y, Ni, Cu) worth exploring in more detail.

The discovery and design of new materials which can efficiently catalyze the oxygen reduction and evolution reactions at reduced temperatures is important for facilitating the widespread adoption of fuel cell and electrolyzer technologies. Numerous studies have produced correlations between catalytic properties, such as oxygen surface exchange or electrode area specific resistance (ASR), and properties of the catalyst material. However, correlations have historically been limited in scope (e.g., using only a few materials or at a single temperature) and it has been difficult to provide detailed assessments of their robustness. Here, we assess the ability of the O p-band center electronic structure descriptor, obtained from density functional theory (DFT) calculations, to correlate with oxygen surface exchange rates, diffusivities, and area specific resistances for a large database of perovskite oxide catalytic properties. By data mining the literature, we obtain 747 catalytic property value data points spanning 299 unique perovskite compositions from 313 studies. We assess linear correlations of each property with the O p-band center and find generally modest correlations that are qualitatively useful (prediction mean absolute errors of about 0.5 log units are typical), where the correlations are improved at higher temperatures (e.g., 800 {\deg}C vs. 500 {\deg}C) and significantly improve when considering fits to the subset of materials which have multiple independent measurements. These findings suggest that the spread of property data is significantly influenced by experimental uncertainty, and subsequent measurements of additional materials will likely improve the O p-band center correlations.

Thermionic electron emission cathodes are critical components of various high power and high frequency vacuum electronic devices, electron microscopes, e-beam lithographic devices, and thermionic energy converters, which all demand an efficient and long-lasting low work function cathode. Single phase, polycrystalline perovskite oxide SrVO$_3$, with its intrinsic low effective work function and facile synthesis process, is a promising cathode candidate, where previous works have shown evidence of an effective work function as low as 2.3 eV. However, assessment of the stability over time under conditions relevant for operation and the related interplay of evolving surface chemistry with emission performance are still missing, and necessary for understanding how to best prepare, process and operate SrVO$_3$ cathodes. In this work, we study the vacuum activation process of SrVO$_3$ and find it has promising emission stability over 15 days of continuous high temperature operation. We find that SrVO$_3$ shows surface Sr and O segregation during operation, which we hypothesize is needed to create a positive surface dipole, leading to low effective work function. Emission repeatability from cyclic heating and cooling suggests the promising stability of the low effective work function surface, and additional observations of drift-free emission during one hour of continuous emission testing at high temperature further demonstrates its excellent performance stability.

Materials discovery and design typically proceeds through iterative evaluation (both experimental and computational) to obtain data, generally targeting improvement of one or more properties under one or more constraints (e.g., time or budget). However, there can be great variation in the quality and cost of different data, and when they are mixed together in what we here call multifidelity data the optimal approaches to their utilization are not established. It is therefore important to develop strategies to acquire and use multifidelity data to realize the most efficient iterative materials exploration. In this work, we assess the impact of using multifidelity data through mock demonstration of designing solar cell materials, using the electronic bandgap as the target property. We propose a new approach of using multifidelity data through leveraging machine learning models of both low- and high-fidelity data, where using predicted low-fidelity data as an input feature in the high-fidelity model can improve the impact of a multifidelity data approach. We show how tradeoffs of low- versus high-fidelity measurement cost and acquisition can impact the materials discovery process, and find that the use of multifidelity data has maximal impact on the materials discovery campaign when approximately five low-fidelity measurements per high-fidelity measurement are performed, and when the cost of low-fidelity measurements is approximately 5% or less than that of high-fidelity measurements. This work provides practical guidance and useful qualitative measures for improving materials discovery campaigns that involve multifidelity data.

In this work, we propose a linear machine learning force matching approach that can directly extract pair atomic interactions from ab initio calculations in amorphous structures. The local feature representation is specifically chosen to make the linear weights a force field as a force/potential function of the atom pair distance. Consequently, this set of functions is the closest representation of the ab initio forces given the two-body approximation and finite scanning in the configurational space. We validate this approach in amorphous silica. Potentials in the new force field (consisting of tabulated Si-Si, Si-O, and O-O potentials) are significantly softer than existing potentials that are commonly used for silica, even though all of them produce the tetrahedral network structure and roughly similar glass properties. This suggests that those commonly used classical force fields do not offer fundamentally accurate representations of the atomic interaction in silica. The new force field furthermore produces a lower glass transition temperature ($T_g\sim$1800 K) and a positive liquid thermal expansion coefficient, suggesting the extraordinarily high $T_g$ and negative liquid thermal expansion of simulated silica could be artifacts of previously developed classical potentials. Overall, the proposed approach provides a fundamental yet intuitive way to evaluate two-body potentials against ab initio calculations, thereby offering an efficient way to guide the development of classical force fields., Comment: 11 pages, 9 figures

An essential aspect of extending safe operation of the active nuclear reactors is understanding and predicting the embrittlement that occurs in the steels that make up the Reactor pressure vessel (RPV). In this work we integrate state of the art machine learning methods using ensembles of neural networks with unprecedented data collection and integration to develop a new model for RPV steel embrittlement. The new model has multiple improvements over previous machine learning and hand-tuned efforts, including greater accuracy (e.g., at high-fluence relevant for extending the life of present reactors), wider domain of applicability (e.g., including a wide-range of compositions), uncertainty quantification, and online accessibility for easy use by the community. These improvements provide a model with significant new capabilities, including the ability to easily and accurately explore compositions, flux, and fluence effects on RPV steel embrittlement for the first time. Furthermore, our detailed comparisons show our approach improves on the leading American Society for Testing and Materials (ASTM) E900-15 standard model for RPV embrittlement on every metric we assessed, demonstrating the efficacy of machine learning approaches for this type of highly demanding materials property prediction.

New highly oxygen-active materials may enhance many energy-related technologies by enabling efficient oxygen-ion transport at lower temperatures, e.g., below 400 Celsius. Interstitial oxygen conductors have the potential to realize such performance but have received far less attention than vacancy-mediated conductors. Here, we combine physically-motivated structure and property descriptors, ab initio simulations, and experiments to demonstrate an approach to discover new fast interstitial oxygen conductors. Multiple new families were found which adopt completely different structures from known oxygen conductors. From these families, we synthesized and studied oxygen kinetics in La4Mn5Si4O22+d (LMS), a representative member of perrierite/chevkinite family. We found LMS has higher oxygen ionic conductivity than the widely used yttria-stabilized ZrO2, and among the highest surface oxygen exchange rates at intermediate temperature of known materials. The fast oxygen kinetics is the result of simultaneously active interstitial and interstitialcy diffusion pathways. This work developed and demonstrated a powerful approach for discovering new families of interstitial oxygen conductors and suggests many more such materials remain to be discovered.

We use a random forest model to predict the critical cooling rate (RC) for glass formation of various alloys from features of their constituent elements. The random forest model was trained on a database that integrates multiple sources of direct and indirect RC data for metallic glasses to expand the directly measured RC database of less than 100 values to a training set of over 2,000 values. The model error on 5-fold cross validation is 0.66 orders of magnitude in K/s. The error on leave out one group cross validation on alloy system groups is 0.59 log units in K/s when the target alloy constituents appear more than 500 times in training data. Using this model, we make predictions for the set of compositions with melt-spun glasses in the database, and for the full set of quaternary alloys that have constituents which appear more than 500 times in training data. These predictions identify a number of potential new bulk metallic glass (BMG) systems for future study, but the model is most useful for identification of alloy systems likely to contain good glass formers, rather than detailed discovery of bulk glass composition regions within known glassy systems.

There has been a growing effort to replace manual extraction of data from research papers with automated data extraction based on natural language processing, language models, and recently, large language models (LLMs). Although these methods enable efficient extraction of data from large sets of research papers, they require a significant amount of up-front effort, expertise, and coding. In this work we propose the ChatExtract method that can fully automate very accurate data extraction with minimal initial effort and background, using an advanced conversational LLM. ChatExtract consists of a set of engineered prompts applied to a conversational LLM that both identify sentences with data, extract that data, and assure the data's correctness through a series of follow-up questions. These follow-up questions largely overcome known issues with LLMs providing factually inaccurate responses. ChatExtract can be applied with any conversational LLMs and yields very high quality data extraction. In tests on materials data we find precision and recall both close to 90% from the best conversational LLMs, like ChatGPT-4. We demonstrate that the exceptional performance is enabled by the information retention in a conversational model combined with purposeful redundancy and introducing uncertainty through follow-up prompts. These results suggest that approaches similar to ChatExtract, due to their simplicity, transferability, and accuracy are likely to become powerful tools for data extraction in the near future. Finally, databases for critical cooling rates of metallic glasses and yield strengths of high entropy alloys are developed using ChatExtract., Comment: 7 pages, 3 figures, 1 table

Accurate and comprehensive material databases extracted from research papers are crucial for materials science and engineering, but their development requires significant human effort. With large language models (LLMs) transforming the way humans interact with text, LLMs provide an opportunity to revolutionize data extraction. In this study, we demonstrate a simple and efficient method for extracting materials data from full-text research papers leveraging the capabilities of LLMs combined with human supervision. This approach is particularly suitable for mid-sized databases and requires minimal to no coding or prior knowledge about the extracted property. It offers high recall and nearly perfect precision in the resulting database. The method is easily adaptable to new and superior language models, ensuring continued utility. We show this by evaluating and comparing its performance on GPT-3 and GPT-3.5/4 (which underlie ChatGPT), as well as free alternatives such as BART and DeBERTaV3. We provide a detailed analysis of the method's performance in extracting sentences containing bulk modulus data, achieving up to 90% precision at 96% recall, depending on the amount of human effort involved. We further demonstrate the method's broader effectiveness by developing a database of critical cooling rates for metallic glasses over twice the size of previous human curated databases., Comment: 13 pages, 4 figures

Fluoride salts are prospective materials for applications in some next generation nuclear reactors and their thermophysical properties at various conditions are of interest. Experimental measurement of the properties of these salts is often difficult and, in some cases, unfeasible due to challenges from high temperatures, impurity control, and corrosivity. Therefore, accurate theoretical methods are needed for fluoride salt property prediction. In this work, we used moment tensor potentials (MTP) to approximate the potential energy surface of eutectic FLiBe (0.66 LiF 0.33 BeF2) predicted by the ab initio (DFT D3) method. We then used the developed potential and molecular dynamics to obtain several thermophysical properties of FLiBe, including radial distribution functions, density, self-diffusion coefficients, thermal expansion, specific heat capacity, bulk modulus, viscosity, and thermal conductivity. Our results show that the MTP potential approximates the potential energy surface accurately and the overall approach yields very good agreement with experimental values. The converged fitting can be obtained with less than 600 configurations generated from DFT calculations, which data can be generated in just 1200 core hours on today's typical processors. The MTP potential is faster than many machine learning potentials and about one order of magnitude slower than widely used empirical molten salt potentials such as Tosi Fumi., Comment: 46 Pages, 21 Figures

Short-timescale atomic rearrangements are fundamental to the kinetics of glasses and frequently dominated by one atom moving significantly (a rearrangement), while others relax only modestly. The rates and directions of such rearrangements (or hops) are dominated by the distributions of activation barriers (Eact) for rearrangement for a single atom and how those distributions vary across the atoms in the system. We have used molecular dynamics simulations of Cu50Zr50 metallic glass below Tg in an isoconfigurational ensemble to catalog the ensemble of rearrangements from thousands of sites. The majority of atoms are strongly caged by their neighbors, but a tiny fraction has a very high propensity for rearrangement, which leads to a power-law variation in the cage-breaking probability for the atoms in the model. In addition, atoms generally have multiple accessible rearrangement vectors, each with its own Eact. However, atoms with lower Eact (or higher rearrangement rates) generally explored fewer possible rearrangement vectors, as the low Eact path is explored far more than others. We discuss how our results influence future modeling efforts to predict the rearrangement vector of a hopping atom.

Nuclear materials are often demanded to function for extended time in extreme environments, including high radiation fluxes and transmutation, high temperature and temperature gradients, stresses, and corrosive coolants. They also have a wide range of microstructural and chemical makeup, with multifaceted and often out-of-equilibrium interactions. Machine learning (ML) is increasingly being used to tackle these complex time-dependent interactions and aid researchers in developing models and making predictions, sometimes with better accuracy than traditional modeling that focuses on one or two parameters at a time. Conventional practices of acquiring new experimental data in nuclear materials research are often slow and expensive, limiting the opportunity for data-centric ML, but new methods are changing that paradigm. Here we review high-throughput computational and experimental data approaches, especially robotic experimentation and active learning that based on Gaussian process and Bayesian optimization. We show ML examples in structural materials ( e.g., reactor pressure vessel (RPV) alloys and radiation detecting scintillating materials) and highlight new techniques of high-throughput sample preparation and characterizations, and automated radiation/environmental exposures and real-time online diagnostics. This review suggests that ML models of material constitutive relations in plasticity, damage, and even electronic and optical responses to radiation are likely to become powerful tools as they develop. Finally, we speculate on how the recent trends in artificial intelligence (AI) and machine learning will soon make the utilization of ML techniques as commonplace as the spreadsheet curve-fitting practices of today.

Speckle patterns produced by coherent X-ray have a close relationship with the internal structure of materials but quantitative inversion of the relationship to determine structure from speckle patterns is challenging. Here, we investigate the link between coherent X-ray speckle patterns and sample structures using a model 2D disk system and explore the ability of machine learning to learn aspects of the relationship. Specifically, we train a deep neural network to classify the coherent X-ray speckle patterns according to the disk number density in the corresponding structure. It is demonstrated that the classification system is accurate for both non-disperse and disperse size distributions.

Accurately quantifying swelling of alloys that have undergone irradiation is essential for understanding alloy performance in a nuclear reactor and critical for the safe and reliable operation of reactor facilities. However, typical practice is for radiation-induced defects in electron microscopy images of alloys to be manually quantified by domain-expert researchers. Here, we employ an end-to-end deep learning approach using the Mask Regional Convolutional Neural Network (Mask R-CNN) model to detect and quantify nanoscale cavities in irradiated alloys. We have assembled the largest database of labeled cavity images to date, which includes 400 images, >34k discrete cavities, and numerous alloy compositions and irradiation conditions. We have evaluated both statistical (precision, recall, and F1 scores) and materials property-centric (cavity size, density, and swelling) metrics of model performance, and performed in-depth analysis of materials swelling assessments. We find our model gives assessments of material swelling with an average (standard deviation) swelling mean absolute error based on random leave-out cross-validation of 0.30 (0.03) percent swelling. This result demonstrates our approach can accurately provide swelling metrics on a per-image and per-condition basis, which can provide helpful insight into material design (e.g., alloy refinement) and impact of service conditions (e.g., temperature, irradiation dose) on swelling. Finally, we find there are cases of test images with poor statistical metrics, but small errors in swelling, pointing to the need for moving beyond traditional classification-based metrics to evaluate object detection models in the context of materials domain applications.

The information content of atomic resolution scanning transmission electron microscopy (STEM) images can often be reduced to a handful of parameters describing each atomic column, chief amongst which is the column position. Neural networks (NNs) are a high performance, computationally efficient method to automatically locate atomic columns in images, which has led to a profusion of NN models and associated training datasets. We have developed a benchmark dataset of simulated and experimental STEM images and used it to evaluate the performance of two sets of recent NN models for atom location in STEM images. Both models exhibit high performance for images of varying quality from several different crystal lattices. However, there are important differences in performance as a function of image quality, and both models perform poorly for images outside the training data, such as interfaces with large difference in background intensity. Both the benchmark dataset and the models are available using the Foundry service for dissemination, discovery, and reuse of machine learning models.

Quantifying charge-state transition energy levels of impurities in semiconductors is critical to understanding and engineering their optoelectronic properties for applications ranging from solar photovoltaics to infrared lasers. While these transition levels can be measured and calculated accurately, such efforts are time-consuming and more rapid prediction methods would be beneficial. Here, we significantly reduce the time typically required to predict impurity transition levels using multi-fidelity datasets and a machine learning approach employing features based on elemental properties and impurity positions. We use transition levels obtained from low-fidelity (i.e., local-density approximation or generalized gradient approximation) density functional theory (DFT) calculations, corrected using a recently proposed modified band alignment scheme, which well-approximates transition levels from high-fidelity DFT (i.e., hybrid HSE06). The model fit to the large multi-fidelity database shows improved accuracy compared to the models trained on the more limited high-fidelity values. Crucially, in our approach, when using the multi-fidelity data, high-fidelity values are not required for model training, significantly reducing the computational cost required for training the model. Our machine learning model of transition levels has a root mean squared (mean absolute) error of 0.36 (0.27) eV vs high-fidelity hybrid functional values when averaged over 14 semiconductor systems from the II-VI and III-V families. As a guide for use on other systems, we assessed the model on simulated data to show the expected accuracy level as a function of bandgap for new materials of interest. Finally, we use the model to predict a complete space of impurity charge-state transition levels in all zinc blende III-V and II-VI systems.

In this work, we use molecular dynamics simulations to study the enhancement of surface over bulk diffusion (surface enhanced diffusion) in (PbO)x(SiO2)1-x glasses. This work is motivated to better understand surface diffusion in glasses and its connection to fragility, and to enhance surface diffusion in silica and related glasses for greater thermodynamic stability during vapor-deposition. By adding PbO to silica, the fragility of glass increases continuously for 10% <= x <= 70% during experiments. The increase in fragility may correspond to an increase in surface enhanced diffusion, as fragility and surface diffusion are correlated. We observe that for the silicates investigated, while surface enhanced diffusion increases with fragility, the enhancement is quite small. The slower diffusing Si and O atoms have higher enhancements, which could allow for some surface stabilization effects. We demonstrate that there are only small changes in atomic arrangements, consistent with the similar diffusion rates, at the surface as compared to bulk. Finally, we examine the trend of bulk versus surface diffusion in view of previous observations in organic and metallic glasses and found that in oxides, fragility increase may not be strongly linked to enhanced surface diffusion.

In this work, we revisit the fragile-to-strong transition (FTS) in the simulated BKS silica from the perspective of microscopic dynamics in an effort to elucidate the dynamical behaviors of fragile and strong glass-forming liquids. Softness, which is a machine-learned feature from local atomic structures, is used to predict the microscopic activation energetics and long-term dynamics. The FTS is found to originate from a change in the temperature dependence of the microscopic activation energetics. Furthermore, results suggest there are two diffusion channels with different energy barriers in BKS silica. The fast dynamics at high temperatures is dominated by the channel with small energy barriers ($<\sim$1 eV), which is controlled by the short-range order. The rapid closing of this diffusion channel when lowering temperature leads to the fragile behavior. On the other hand, the slow dynamics at low temperatures is dominated by the channel with large energy barriers controlled by the medium-range order. This slow diffusion channel changes only subtly with temperature, leading to the strong behavior. The distributions of barriers in the two channels show different temperature dependences, causing a crossover at $\sim$3100 K. This transition temperature in microscopic dynamics is consistent with the inflection point in the configurational entropy, suggesting there is a fundamental correlation between microscopic dynamics and thermodynamics., Comment: 9 pages, 5 figures, additional 6 pages of supplementary materials

In a current density versus temperature (J-T) (Miram) curve in thermionic electron emission, experimental measurements demonstrate there is a smooth transition between the exponential region and the saturated emission regions, which is sometimes referred to as the "roll-off" or "Miram curve knee". The shape of the Miram curve knee is an important figure of merit for thermionic vacuum cathodes. Specifically, cathodes with a sharp Miram curve knee at low temperature with a flat saturated emission current are typically preferred. Our previous work on modeling nonuniform thermionic emission revealed that the space charge effect and patch field effect are key pieces of physics which impact the shape of the Miram curve knee. This work provides a more complete understanding of the physical factors connecting these physical effects and their relative impact on the shape of the knee, including the smoothness, the temperature, and the flatness of the saturated emission current density. For our analyses, we use a periodic, equal-width striped ("zebra crossing") work function distribution as a model system and illustrate how the space charge and patch field effects restrict the emission current density near the Miram curve knee. The results indicate there are three main physical parameters which significantly impact the shape of the Miram curve. Such physical knowledge directly connects the patch size, work function values, anode-cathode voltage, and anode-cathode gap distance to the shape of the Miram curve, providing new understanding and a guide to the design of thermionic cathodes used as electron sources in vacuum electronic devices (VEDs).

This paper provides an accurate theoretical defect energy database for pure and Bi-containing III-V (III-V:Bi) materials and investigates efficient methods for high-throughput defect calculations based on corrections of results obtained with local and semi-local functionals. Point defects as well as nearest-neighbor and second-nearest-neighbor pair defects were investigated in charge states ranging from -5 to 5. Ga-V:Bi systems (GaP:Bi, GaAs:Bi, and GaSb:Bi) were thoroughly investigated with significantly slower, higher fidelity hybrid Heyd-Scuseria-Ernzerhof (HSE) and significantly faster, lower fidelity local density approximation (LDA) calculations. In both approaches spurious electrostatic interactions were corrected with the Freysoldt correction. The results were verified against available experimental results and used to assess the accuracy of a previous band alignment correction. Here, a modified band alignment method is proposed in order to better predict the HSE values from the LDA ones. The proposed method allows prediction of defect energies with values that approximate those from the HSE functional at the computational cost of LDA (about 20x faster for the systems studied here). Tests of selected point defects in In-V:Bi materials resulted in corrected LDA values having a mean absolute error (MAE)=0.175 eV for defect levels vs. HSE. The method was further verified on an external database of defects and impurities in CdX (X=S, Se, Te) systems, yielding a MAE=0.194 eV. These tests demonstrate the correction to be sufficient for qualitative and semi-quantitative predictions, and may suggest transferability to many semiconductor systems without significant loss in accuracy. Properties of the remaining In-V:Bi defects and all Al-V:Bi defects were predicted with the use of the modified band alignment method.

Experimental observations of thermionic electron emission demonstrate a smooth transition between TL and FSCL regions of the emitted-current-density-versus-temperature (J-T) (Miram) curve and the emitted-current-density-versus-voltage (J-V) curve. Knowledge of the temperature and shape of the TL-FSCL transition is important in evaluating the thermionic electron emission performance of cathodes, including predicting the lifetime. However, there have been no first-principles physics-based models that can predict the smooth TL-FSCL transition region for real thermionic cathodes without applying physically difficult to justify a priori assumptions or empirical phenomenological equations. Previous work detailing the nonuniform thermionic emission found that the effects of 3-D space charge, patch fields, and Schottky barrier lowering can lead to a smooth TL-FSCL transition region from a model thermionic cathode surface with a checkerboard spatial distribution of work function values. In this work, we construct a physics-based nonuniform emission model for commercial dispenser cathodes for the first time. This emission model is obtained by incorporating the cathode surface grain orientation via electron backscatter diffraction (EBSD) and the facet-orientation-specific work function values from density functional theory (DFT) calculations. The model enables construction of two-dimensional emitted current density maps of the cathode surface and corresponding J-T and J-V curves. The predicted emission curves show excellent agreement with experiment, not only in TL and FSCL regions but, crucially, also in the TL-FSCL transition region. This model improves the understanding on the relationship between thermionic emission and cathode microstructure, which is beneficial to the design of vacuum electronic devices.

In this work, we perform semantic segmentation of multiple defect types in electron microscopy images of irradiated FeCrAl alloys using a deep learning Mask Regional Convolutional Neural Network (Mask R-CNN) model. We conduct an in-depth analysis of key model performance statistics, with a focus on quantities such as predicted distributions of defect shapes, defect sizes, and defect areal densities relevant to informing modeling and understanding of irradiated Fe-based materials properties. To better understand the performance and present limitations of the model, we provide examples of useful evaluation tests which include a suite of random splits, and dataset size-dependent and domain-targeted cross validation tests. Overall, we find that the current model is a fast, effective tool for automatically characterizing and quantifying multiple defect types in microscopy images, with a level of accuracy on par with human domain expert labelers. More specifically, the model can achieve average defect identification F1 scores as high as 0.8, and, based on random cross validation, have low overall average (+/- standard deviation) defect size and density percentage errors of 7.3 (+/- 3.8)% and 12.7 (+/- 5.3)%, respectively. Further, our model predicts the expected material hardening to within 10-20 MPa (about 10% of total hardening), which is about the same error level as experiments. Our targeted evaluation tests also suggest the best path toward improving future models is not expanding existing databases with more labeled images but instead data additions that target weak points of the model domain, such as images from different microscopes, imaging conditions, irradiation environments, and alloy types. Finally, we discuss the first phase of an effort to provide an easy-to-use, open-source object detection tool to the broader community for identifying defects in new images.

We explore the use of characteristic temperatures derived from molecular dynamics to predict aspects of metallic Glass Forming Ability (GFA). Temperatures derived from cooling curves of self-diffusion, viscosity, and energy were used as features for machine learning models of GFA. Multiple target and model combinations with these features were explored. First, we use the logarithm of critical casting thickness, $log_{10}(D_{max})$, as the target and trained regression models on 21 compositions. Application of 3-fold cross-validation on the 21 $log_{10}(D_{max})$ alloys showed only weak correlation between the model predictions and the target values. Second, the GFA of alloys were quantified by melt-spinning or suction casting amorphization behavior, with alloys that showed crystalline phases after synthesis classified as Poor GFA and those with pure amorphous phases as Good GFA. Binary GFA classification was then modeled using decision tree-based methods (random forest and gradient boosting models) and were assessed with nested-cross validation. The maximum F1 score for the precision-recall with Good Glass Forming Ability as the positive class was $0.82 \pm 0.01$ for the best model type. We also compared using simple functions of characteristic temperatures as features in place of the temperatures themselves and found no statistically significant difference in predictive abilities. Although the predictive ability of the models developed here are modest, this work demonstrates clearly that one can use molecular dynamics simulations and machine learning to predict metal glass forming ability.

Materials with low work functions are critical for an array of applications requiring the facile removal or efficient transport of electrons through a device. Perovskite oxides are a promising class of materials for finding low work functions, and here we target applications in thermionic and field electron emission. Perovskites have highly malleable compositions which enable tunable work function values over a wide range, robust stability at high temperatures, and high electronic conductivities. In this work, we screened over 2900 perovskite oxides in search of stable, conductive, low-work-function materials using Density Functional Theory (DFT) methods. Our work provides insight into the materials chemistry governing the work function value of a perovskite, where materials with barely filled d bands possess the lowest work functions. Our screening has resulted in a total of seven promising compounds, such as BaMoO3 and SrNb0.75Co0.25O3 with work functions of 1.1 eV and 1.5 eV, respectively. These promising materials and others presented in this study may find use as low work function electron emitters in high power vacuum electronic and thermionic energy conversion devices. Moreover, the database of calculated work functions and materials chemistry trends governing the value of the work function may aid in the engineering of perovskite heterojunction devices.

LaB6 has been used as a commercial electron emitter for decades. Despite the large number of studies on the work function of LaB6, there is no comprehensive understanding of work function trends in the hexaboride materials family. In this study, we use Density Functional Theory (DFT) calculations to calculated trends of rare earth hexaboride work function and rationalize these trends based on the electronegativity of the metal element. We predict that alloying LaB6 with Ba can further lower the work function by ~0.2 eV. Interestingly, we find that alloyed (La, Ba)B6 can have lower work functions than either LaB6 or BaB6, benefitting from an enhanced surface dipole due to metal element size mismatch. In addition to hexaborides we also investigate work function trends of similar materials families, namely tetraborides and transition metal nitrides, which, like hexaborides, are electrically conductive and refractory and thus may also be promising materials for electron emission applications. We find that tetraborides consistently have higher work functions than their hexaboride analogues as the tetraborides having less ionic bonding and smaller positive surface dipoles. Finally, we find that HfN has a low work function of about 2.2 eV, making HfN a potentially promising new electron emitter material.

MAST-SEY is an open-source Monte Carlo code capable of calculating secondary electron emission using input data generated entirely from first principle (density functional theory) calculations. It utilizes the complex dielectric function and Penn's theory for inelastic scattering processes, and relativistic Schr\"odinger theory by means of a partial-wave expansion method to govern elastic scattering. It allows the user to include explicitly calculated momentum dependence of the dielectric function, as well as to utilize first-principle density of states in secondary electron generation, which provides a more complete description of the underlying physics. In this paper we thoroughly describe the theoretical aspects of the modeling, as used in the code, and present sample results obtained for copper and aluminum.

The bulk and surface dynamics of Cu50Zr50 metallic glass were studied using classical molecular dynamics (MD) simulations. As the alloy undergoes cooling, it passes through liquid, supercooled, and glassy states. While bulk dynamics showed a marked slowing down prior to glass formation, with increasing activation energy, the slowdown in surface dynamics was relatively subtle. The surface exhibited a lower glass transition temperature than the bulk, and the dynamics preceding the transition were accurately described by a temperature-independent activation energy. Surface dynamics were much faster than bulk at a given temperature in the supercooled state, but surface and bulk dynamics were found to be very similar when compared at their respective glass transition temperatures. The manifestation of dynamical heterogeneity, as characterized by the non-Gaussian parameter and breakdown of the Stokes-Einstein equation, was also similar between bulk and surface for temperatures scaled by their respective glass transition temperatures. Individual atom motion was dominated by a cage and jump mechanism in the glassy state for both the bulk and surface. We utilize this cage and jump mechanisms to separate the activation energy for diffusion into two parts: (i) cage-breaking barrier (Q1), associated with the rearrangement of neighboring atoms to free up space and (ii) the subsequent jump barrier (Q2). It was observed that Q1 dominates Q2 for both bulk and surface diffusion, and the difference in activation energies for bulk and surface diffusion mainly arose from the differences in cage-breaking barrier Q1., Comment: Manuscript and Supplementary Information all as a single PDF

Electron microscopy is widely used to explore defects in crystal structures, but human detecting of defects is often time-consuming, error-prone, and unreliable, and is not scalable to large numbers of images or real-time analysis. In this work, we discuss the application of machine learning approaches to find the location and geometry of different defect clusters in irradiated steels. We show that a deep learning based Faster R-CNN analysis system has a performance comparable to human analysis with relatively small training data sets. This study proves the promising ability to apply deep learning to assist the development of automated microscopy data analysis even when multiple features are present and paves the way for fast, scalable, and reliable analysis systems for massive amounts of modern electron microscopy data.

Videos captured using Transmission Electron Microscopy (TEM) can encode details regarding the morphological and temporal evolution of a material by taking snapshots of the microstructure sequentially. However, manual analysis of such video is tedious, error-prone, unreliable, and prohibitively time-consuming if one wishes to analyze a significant fraction of frames for even videos of modest length. In this work, we developed an automated TEM video analysis system for microstructural features based on the advanced object detection model called YOLO and tested the system on an in-situ ion irradiation TEM video of dislocation loops formed in a FeCrAl alloy. The system provides analysis of features observed in TEM including both static and dynamic properties using the YOLO-based defect detection module coupled to a geometry analysis module and a dynamic tracking module. Results show that the system can achieve human comparable performance with an F1 score of 0.89 for fast, consistent, and scalable frame-level defect analysis. This result is obtained on a real but exceptionally clean and stable data set and more challenging data sets may not achieve this performance. The dynamic tracking also enabled evaluation of individual defect evolution like per defect growth rate at a fidelity never before achieved using common human analysis methods. Our work shows that automatically detecting and tracking interesting microstructures and properties contained in TEM videos is viable and opens new doors for evaluating materials dynamics.