Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Composite models”

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 1,333 records · Page 74

A physics-based model of cladding wastage layer formation rate

In this report, lower length scale simulations to inform an engineering-scale model of fuel-cladding chemical interaction (FCCI) conducted under the auspices of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY22 are described. An overall strategy for the implementation of the BISON model is described, and the past and current lower length scale work on FCCI formation is put into the context of this overall strategy. Density functional theory and kinetic Monte Carlo simulations are used to determine the diffusion coefficient of neodymium in iron. Density functional theory calculations are also used to parameterize a parabolic approximation for the dependence of free energy on composition for the intermetallic compound Fe 17 Nd 2 . A phase-field model of the Fe-Nd system was developed that includes the random solid solution body-centered cubic phase and the intermetallic Fe 17 Nd 2 . The model was used to simulate growth of the intermetallic layer in a diffusion couple and the results were compared to experiment. The model’s prediction of the parabolic growth constant is within 40 % of the experimental value. The model is also used to simulate the formation of a grain boundary attack layer that is observed in experiment. The simulations show enhanced Nd concentration along the grain boundaries, but do not clearly show formation of a Fe 17 Nd 2 layer

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A separations and purification process for improving yields and meeting fuel contaminant specifications for high-octane gasoline produced from dimethyl-ether over a Cu/BEA catalyst

In this work, we have been developing a three-step conversion of biomass-derived syngas to methanol to dimethyl-ether (DME) to non-aromatic hydrocarbons for use as high-octane gasoline and sustainable aviation fuel. This process produces branched alkanes from DME using a Cu/BEA catalyst and is a promising alternative to other syngas conversion processes such as Fischer-Tropsch to linear alkanes and traditional ZSM-5 catalyzed methanol to aromatic gasoline. In this short article we describe some advances in our understanding related to separations and purification via the use of more detailed experimental speciation in an updated process model involving multiple phase equilibrium-based separation steps. Primary modeled reactor outlet constituents (and weight %) are: C3 and lighter hydrocarbon gases (11.1%), C4s (54.5%), H 2 (1.2%), CO 2 (2.9%), water (5.0%), unreacted DME (16.5%), methanol (2.3%), and C5+ hydrocarbons (6.4%). DME (the primary reactant) and H 2 recycle and reuse are important for the overall process efficiency, and the recycle of C4s is important to increase the C5+ yield via reactivation and homologation. Thus H 2 , C4s, and DME are targeted for recycle, while methanol and water need to be removed from the product to conform with fuel specifications. Model predictions from Aspen Plus using the NRTL-RK property method indicate a fuel composition with C5+ content of 97.1 wt%, with minor constituents: 2.4 wt% C4s, 0.3 wt% methanol, 0.1 wt% DME, 0.03 wt% water, and 0.01 wt% C3s. These ranges of minor components conform with fuel quality requirements, and the modeled product is amenable for unconstrained blending to boost gasoline octane ratings.

09 BIOMASS FUELS↗

Systems Engineering of Rhodococcus opacus to Enable Production of Drop-in Fuels from Lignocellulose

Production of drop-in fuels from lignocellulose using Rhodococcus opacus PD630 (hereafter R. opacus) is a challenging goal. During the grant period we have pushed the field forward significantly in several areas of research. Towards the end goal of accelerating the adoption of R. opacus in biofuel production, during the grant period we have expanded the phenotypic characterization of R. opacus grown in single aromatic (model lignocellulosic) compounds or their mixtures, modeling the growth conditions in lignocellulosic biomass. Harnessing the power of adaptive evolution, we produced evolved R. opacus isolates with superior lignin valorization capabilities and identified differentially expressed genes and pathways after adaptation. We used next generation multi-omic techniques such as genomic, transcriptomic, and metabolomic analyses, to identify the catabolic pathways used by R. opacus to degrade aromatic compounds and funnel these degradation products into central metabolism, as well as the aromatic transport genes required for increased tolerance and utilization. Taking this information one step further, we identified endogenous transcription factors and regulatory mechanisms important for degradation of five model aromatic compounds. To accurately estimate R. opacus growth and consumption on model lignin compounds we pioneered the use of novel extraction procedures prior to GC-MS analysis. Alongside 13 C-metabolic flux analysis, we have elucidated the metabolic routes preferred by Rhodococcus opacus during aromatic compound degradation. Finally, we used in tandem lipidomics and high-resolution mass spectrometry to identify the modulation of mycolic acids and phospholipid membrane composition modification as a strategy for aromatic tolerance in R. opacus. Being a non-model organism, R. opacus lacks the breadth of tools and technical foundation which drive biofuel research in more well-understood microbes such as Escherichia coli. To reduce this burden for use, we designed and produced new tools for genomic manipulation and engineering in R. opacus. These engineering breakthroughs support efficient genomic editing, enabling gene overexpression, repression, and genetic alteration. Using these tools, we have generated synthetically engineered strains with increased lipogenesis and growth, both positive traits required for increased lignin valorization. Optimizing engineered strains for biofuel production from lignocellulose requires extremely sophisticated synthetic rewiring of metabolism. To facilitate systems-level reorganization of metabolism in R. opacus, we created a genome-scale model that accurately predicts metabolic flux and growth rates on the aromatic compound phenol. Lignin requires extensive pre-treatment before biological degradation by R. opacus. Towards an eventual goal of degrading real-world lignin, we developed new depolymerization processes to generate lignin breakdown products (LBP). We optimized LBP storage and composition analysis techniques, enabling accurate prediction of specific LBP compound integration into cell wall components. Overall, through the work funded by this grant we generated 20 manuscripts (17 published, 3 in review/preparation), methods for increased accuracy in metabolomics of aromatic compounds, multiple genetic tools for altering the R. opacus genome, genome scale models for predicting flux through metabolic pathways, as well as multi-omic data for community use. The work funded by this grant has increased the knowledge of aromatic degradation in bacteria and advanced our efforts to optimize R. opacus for lignin valorization.

09 BIOMASS FUELS↗

The effects of introducing elasticity using different interpolation schemes to the grand potential phase field model

Introducing elastic energy in the phase field method has been shown to influence interfacial energy, depending on the elastic interpolation scheme. This study investigates the impact of the elastic energy when using a grand potential-based phase field method, comparing the result of Khachaturyan’s strain interpolation scheme (KHS) and Voight-Taylor’s elastic energy interpolation scheme (VTS). The KHS model leads to a decrease in the interfacial energy, while the VTS model leads to an increase. The change in interfacial energy is greater with the VTS model than the KHS model, which suggests that the KHS model is more appropriate to limit the artificial impact of the elastic energy on the interfacial energy. When the contribution at the interface is not negligible, it is shown that both the microstructure evolution kinetics and the equilibrium microstructure can be influenced by the choice of the elastic scheme being used. In addition, this paper shows that the grand potential model might not be appropriate when the system requires the introduction of a composition-dependent term in the elastic energy contribution. This limitation is due to the need for an explicit and invertible relation between the total potential and the composition.

36 MATERIALS SCIENCE↗

Two-neutron knockout as a probe of the composition of states in Mg 22 , Al 23 , and Si 24

Simpson and Tostevin proposed that the width and shape of exclusive parallel momentum distributions of the A - 2 residue in direct two-nucleon knockout reactions carry a measurable sensitivity to the nucleon single-particle configurations and their couplings within the wave functions of exotic nuclei. We further report here on the first benchmarks and use of this new spectroscopic tool. Exclusive parallel momentum distributions for states in the neutron-deficient nuclei Mg 22 , Al 23 , and Si 24 populated in such direct two-neutron removal reactions were extracted and compared to predictions combining eikonal reaction theory and shell-model calculations. For the well-known Mg 22 and Al 23 nuclei, measurements and calculations were found to agree, supporting the dependence of the parallel momentum distribution width on the angular momentum composition of the shell-model two-neutron amplitudes. In Si 24 , a level at 3439(9) keV, of relevance for the important Al 23 ( p , γ ) Si 24 astrophysical reaction rate, was confirmed to be the 2 2 + state, whereas the 4 1 + state, expected to be strongly populated in two-neutron knockout, was not observed. This puzzle is resolved by theoretical considerations of the Thomas-Ehrman shift, which also indicates that a previously reported 3471-keV state in Si 24 is, in fact, the ( 0 2 + ) level with one of the largest experimental mirror-energy shifts ever observed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A FFT-based mesoscale continuum dislocation mechanics with defect energy: Applications to composites and polycrystals

A crystal plasticity elastoviscoplastic FFT (fast Fourier transform) formulation with a mesoscale continuum field dislocation mechanics model is presented, which incorporates a defect energy density that depends on GND densities and an associated material length scale. This allows to thermodynamically derive internal length scale dependent intra-crystalline backstress and Peach–Koehler force acting on GND densities. The model considers GND density evolution through a filtered numerical spectral approach, which is coupled with stress equilibrium through the elastoviscoplastic FFT algorithm. The discrete Fourier transform (DFT) method together with finite difference (FD) schemes is applied to solve both the backstress tensor and the Fourier–Green operator. Numerical results are first reported for two-phase laminate composites with plastic single crystal channels and elastic precipitates for shear loadings. Channel size effects are simulated and analyzed on the overall and local hardening behaviors during monotonous loadings. In addition, the evolutions of GND densities and the role of their associated backstress on size effects are examined during reversible shear loading. In a second part, the role of the defect energy internal length scale on polycrystal’s hardening during tension–compression is discussed. The results are compared to those obtained using FFT-based continuum field dislocation mechanics without defect energy.

36 MATERIALS SCIENCE↗

Counter Data Paucity through Adversarial Invariance Encoding: A Case Study on Modeling Battery Thermal Runaway

Lithium-ion batteries, widely used for their durability and high energy storage, face the risk of internal short circuits leading to catastrophic thermal runaway events. These events, triggered by external stimuli like mechanical loads, pose safety concerns in applications such as electric vehicles. Detecting and understanding thermal runaway events is crucial, but physics-driven models struggle to explain the non-linear evolution of battery temperature during these events, considering factors like material composition and state-of-charge. Due to the rarity of these events and the cost of data collection, we propose a deep learning (DL) model to predict battery temperature responses during thermal runaway. The challenge lies in the scarcity of data, making traditional DL models prone to overfitting and learning low-quality representations of the complex process.Our approach introduces a novel few-shot architecture that incorporates an adversarially governed invariant encoding process. This architecture aims to distill "invariant" relationships by addressing distributional shifts in data across various battery properties, facilitating the detection of thermal runaway events. Specifically, our results demonstrate that deep learning models conditioned on these "invariant" representations outperform state-of-the-art baselines, achieving a remarkable 96.8% performance improvement in terms of the popular metric MAPE. This framework presents a promising direction for enhancing battery safety modeling, particularly in the context of rare and complex events like thermal runaway. Our code and code and dataset used for the paper are public1.

Tabassum, Anika [ORNL] (ORCID:0000000254600955)↗

Convergence of Weak-SINDy Surrogate Models

In this paper, we give an in-depth error analysis for surrogate models generated by a variant of the Sparse Identification of Nonlinear Dynamics (SINDy) method. We start with an overview of a variety of nonlinear system identification techniques, namely SINDy, weak-SINDy, and the occupation kernel method. Under the assumption that the dynamics are a finite linear combination of a set of basis functions, these methods establish a linear system to recover coefficients. We illuminate the structural similarities between these techniques and establish a projection property for the weak-SINDy technique. Following the overview, we analyze the error of surrogate models generated by a simplified version of weak-SINDy. In particular, under the assumption of boundedness of a composition operator given by the solution, we show that (i) the surrogate dynamics converges towards the true dynamics and (ii) the solution of the surrogate model is reasonably close to the true solution. Finally, as an application, we discuss the use of a combination of weak-SINDy surrogate modeling and proper orthogonal decomposition (POD) to build a surrogate model for partial differential equations (PDEs).

97 MATHEMATICS AND COMPUTING↗

A general mechanism for gel layer formation on borosilicate glass under aqueous corrosion

Mineral and glass dissolution is a scientific topic deeply investigated but incompletely understood and of a great interest for the geochemical and materials science communities. If the interfacial dissolution/reprecipitation mechanism seems to be applicable to most of silicate minerals, the debate remains open concerning glass. Here we studied two model glasses, a ternary borosilicate (CJ1) and the same glass doped with 4.1 mol% of Al2O3 (CJ2). The two glasses were altered at 90°C, pH 9, and in conditions far and close to saturation with respect to amorphous silica, to determine the initial and residual rates. Moreover, a specific experiment was conducted for a short duration with a solution highly enriched with 18O and 29Si isotopes to understand how passivating gels form. SEM, TEM and ToF-SIMS characterization, along with Monte Carlo simulations were used to understand the rate limiting reactions at play and infer the role of Al. We show that Al yields a slower matrix dissolution in dilute conditions. However, it slows down he formation and the maturation of the passivating gel and favors alteration by partial hydrolysis of Si and Al entities followed by in-situ reorganization/relaxation into a porous network. Unexpectedly, CJ1 experienced both interfacial dissolution/reprecipitation and partial hydrolysis followed by in situ reorganization of the silicate network during the course of a single experiment. This study offers a unified concept that can pave the way for the future development of a predictive kinetic model based on a detailed description of bond breaking and bond forming as a function of glass composition and alteration conditions.

Gin, Stephane↗

A hybrid calorimetry-simulation model of mixing enthalpy for molten salt

Calorimetric determination of enthalpies of mixing (ΔH mix ) in multicomponent molten salts is often interpreted using empirical models that lack physically meaningful parameters. However, for improving pyrochemical separation of spent nuclear fuel, where lanthanides are major fission products and critical elements, a deeper thermodynamic understanding of the link between excess thermodynamic properties and solvation structure is critically needed. In this work, we implement a hybrid and physics-informed framework, MIVM+Calorimetry+AIMD, which integrates experimentally measured ΔH mix (via high temperature drop calorimetry) with solvation structures from ab initio molecular dynamics (AIMD). This approach is demonstrated using LaCl 3 mixed with eutectic LiCl-KCl (58 mol% – 42 mol%) at 873 K and 1133 K. MIVM-derived parameters enable extrapolation of excess Gibbs energy and La 3+ activity across compositions. In contrast, direct ΔH mix predictions from AIMD and polarizable ion model simulations deviate significantly. By incorporating experimentally benchmarked solvation structures into an interpretable thermodynamic model, the MIVM+Calorimetry+AIMD formalism achieves higher accuracy and generalizable method for studying molten salts, offering a robust path for understanding and optimizing molten salt chemistry relevant to nuclear fuel cycles and separation science.

Goncharov, Vitaliy G. [Washington State Univ., Pul↗

Integrating very-high-resolution UAS data and airborne imaging spectroscopy to map the fractional composition of Arctic plant functional types in Western Alaska

Widespread changes in vegetation cover and composition are driving strong impacts on Arctic ecosystem functioning and global climate feedbacks. An accurate characterization of tundra vegetation composition is required to understand how the Arctic will respond to future climate change. However, quantifying tundra vegetation composition over large areas is challenging as commonly-used satellite observations are too coarse, spatially and spectrally, to differentiate low-lying tundra vegetation types. Recent airborne and spaceborne imaging spectroscopy platforms provide better data to characterize vegetation composition. Yet, our ability to characterize vegetation composition with imaging spectroscopy remains largely unexplored in the Arctic, particularly due to a lack of ground observations needed to train and test classification models. To address this problem, we collected very-high-resolution (VHR, ~5 cm) unoccupied aerial system (UAS) imagery at three low-Arctic tundra sites located on the Seward Peninsula, western Alaska. In this paper, we examine the feasibility of integrating imagery from the UAS and the hyperspectral Airborne Visible/Infrared Imaging Spectrometer, Next Generation (AVIRIS-NG) airborne instrument to map the fractional composition of 12 key Arctic plant functional types (PFTs). To this end, we first mapped the 12 PFTs from our VHR UAS imagery using random forest classification. We then used these UAS-derived PFT maps as ground truth to develop partial least squares regression (PLSR) models to predict the fractional cover (FCover) of each PFT from AVIRIS-NG imagery. Further, we evaluated the performance of our PLSR models using reserved UAS samples, as well as by mapping PFT FCover and dominant PFT for large tundra landscapes. Our results show that 1) Arctic PFTs can be effectively mapped using VHR UAS imagery, with overall accuracy between 86% and 92%, 2) when the UAS mapped PFTs were used to inform PLSR scaling models, the FCover of the 12 PFTs could be effectively estimated from AVIRIS-NG imagery with a mean absolute error (MAE) <0.13, and 3) our PLSR models outperformed traditional, fully constrained least-squares (FCLS) linear mixture analysis and produced high-quality, spatially contiguous PFT FCover and PFT maps that captured vegetation spatial patterns with similar accuracy to those developed from UAS imagery. The developed PLSR models have the potential to be broadly applied for quantifying vegetation composition with AVIRIS-NG images to help monitor tundra vegetation dynamics and improve process-based modeling of tundra ecosystems.

54 ENVIRONMENTAL SCIENCES↗

A detailed genome-scale metabolic model of Clostridium thermocellum investigates sources of pyrophosphate for driving glycolysis

Lignocellulosic biomass is an abundant and renewable source of carbon for chemical manufacturing, yet it is cumbersome in conventional processes. A promising, and increasingly studied, candidate for lignocellulose bioprocessing is the thermophilic anaerobe Clostridium thermocellum given its potential to produce ethanol, organic acids, and hydrogen gas from lignocellulosic biomass under high substrate loading. Possessing an atypical glycolytic pathway which substitutes GTP or pyrophosphate (PP i ) for ATP in some steps, including in the energy-investment phase, identification, and manipulation of PP i sources are key to engineering its metabolism. Previous efforts to identify the primary pyrophosphate have been unsuccessful. Here, we explore pyrophosphate metabolism through reconstructing, updating, and analyzing a new genome-scale stoichiometric model for C. thermocellum, iCTH669. Hundreds of changes to the former GEM, iCBI655, including correcting cofactor usages, addressing charge and elemental balance, standardizing biomass composition, and incorporating the latest experimental evidence led to a MEMOTE score improvement to 94%. We found agreement of iCTH669 model predictions across all available fermentation and biomass yield datasets. The feasibility of hundreds of PP i synthesis routes, newly identified and previously proposed, were assessed through the lens of the iCTH669 model including biomass synthesis, tRNA synthesis, newly identified sources, and previously proposed PP i -generating cycles. In all cases, the metabolic cost of PP i synthesis is at best equivalent to investment of one ATP suggesting no direct energetic advantage for the cofactor substitution in C. thermocellum. Even though no unique source of PP i could be gleaned by the model, by combining with gene expression data two most likely scenarios emerge. First, previously investigated PP i sources likely account for most PP i production in wild-type strains. Second, alternate metabolic routes as encoded by iCTH669 can collectively maintain PP i levels even when previously investigated synthesis cycles are disrupted. Model iCTH669 is available at github.com/maranasgroup/iCTH669.

59 BASIC BIOLOGICAL SCIENCES↗

Chromium versus Aluminum: Impact of Nickel Alloy Composition and Interfacial Kinetics on High-Temperature Passivating Oxide Formation

High-temperature corrosion resistance depends critically on the formation of a passivating surface oxide, which is highly sensitive to alloy composition and structure. Such details often elude experimental investigation, and simplified analytical models fail to provide a truly chemical view of passivating oxide evolution. Here, we explicitly compare the fundamental chemistry of Cr and Al as prototypical passivating elements in Ni alloys by directly simulating competing reaction and diffusion processes within the oxide film using kinetic Monte Carlo and density functional theory. We find that the origin and expression of passivating behavior during early-stage thermal oxidation are qualitatively different between the two alloy systems. Ni–Cr alloys feature a sudden onset of passivation associated with a sharp phase transition upon Cr enrichment that directly couples oxidation kinetics to phase transformation behavior. In contrast, Ni–Al alloys display more continuous oxide phase variation with Al enrichment, ultimately resulting in a lower composition threshold for passivation and a thinner passivating layer. In addition, we elucidate the nonobvious role of metal exchange within the alloy near the oxide boundary, which fundamentally alters film composition and passivating behavior. Furthermore, our results have key implications for engineering improved corrosion-resistant alloys, both in terms of compositional variation and processing.

Alloys↗

How do soil microbes shape ecosystem biogeochemistry in the context of global change?

Earth system models disagree on the future of global soil carbon C, with projections ranging from losses of 70 Pg to gains of 250 Pg by the year 2100 (Todd-Brown et al., 2014). This disagreement is due, in part, to the high degree of uncertainty in the mechanisms that control soil C, the largest dynamic pool of organic C on Earth. It remains unclear how climate change, rising atmospheric CO2 levels, and other environmental changes will combine to affect global soil C. As key engineers of biogeochemical cycling, microbes play a crucial role in soil C responses to these environmental changes. Yet, the path toward representing microbes in soil models remains uncertain. Furthermore, we address four burning questions relevant for modelling microbial ecological, evolutionary, and biogeochemical processes in soils: (1) What is the role of soil microbes in biogeochemical cycling? (2) Does microbial community composition matter for soil carbon cycling? (3) Which ecological and evolutionary processes contribute to functional changes? (4) How should we model microbial eco-evolutionary mechanisms in soil C models?

54 ENVIRONMENTAL SCIENCES↗

Using Complex Probability Amplitudes to Simulate Solute Transport in Composite Porous Media

Probability amplitudes are fundamental to quantum mechanics and offer robust descriptions of complicated systems, which have allowed physicists to explain behaviors inaccessible to classical physics. This article ponders how some of the same conceptual underpinnings of the mathematics used for modeling quantum systems might be applied to subsurface water resources problems and speculates how these tools could facilitate applications on quantum computers. A probability amplitude-based model for describing advective-dispersive transport in porous media using linear operators is investigated. The proposed complex valued model decomposes spreading into two “sub-continuum partial dispersion” coefficients, and this recovers classical spreading when the sum of these coefficients is the Fickian dispersion coefficient. However, the probability amplitudes have a manyto-one relationship to a probability distribution, so it embeds a level of heterogeneity into seemingly equivalent functions. Two propagators with different sub-continuum coefficients may have the same macroscopic behavior when either is considered in isolation, but when they act on the other the system’s behavior changes. Additionally, differences in the amplitudes cause a reduction in spreading as velocity correlations are disrupted, despite both propagators having identical dispersion coefficients, and this cannot be achieved using classical methods without changing the dispersion coefficient. The main point is that these amplitude-based models offer a way to embed information about the system into the propagators, instead of just “averaging it out” when making an upscaled model.

54 ENVIRONMENTAL SCIENCES↗

Regeneration of Active Surface Alloys during Cyclic Oxidation and Reduction: Oxidation of H 2 on Pd/Ag(111)

The surface morphology and composition of a catalyst during excursions between oxidizing and reducing conditions can change substantially, especially in bimetallic alloys. Both thermodynamic and kinetic factors play a role in determining the properties of alloy surfaces where the active phase may be a metastable state. Previously, Ag oxide reduction was shown to be dramatically enhanced when Pd is on the surface; however, Pd is more stable when dissolved in Ag, raising the question as to whether a highly active Pd surface state will persist over multiple reaction cycles—a requirement for catalytic function. Experiments herein demonstrate that the enhanced chemical functionality due to the presence of Pd on the surface is retained, based on the enhanced rate of silver oxide reduction over multiple oxidation/reduction cycles for a Pd/Ag(111) model. Repeated oxidation and reduction promote PdAg alloying and reversible structural and compositional changes are detected using X-ray photoelectron spectroscopy. Furthermore, this study establishes that metastable phases can persist in reactive processes on surfaces, indicating their potential in heterogeneous catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Applying Machine Learning to the Classification of DC-DC Converters (Milestone 2 Deliverable Report)

Since extending the Autodetector to a convolutional neural network (CNN) machine learning classifier model, an effort has been executed to demonstrate its ability to distinguish not only a switching DC-DC converter as high voltage, but identify the make and model of a converter on which it was trained. This was achieved by collecting data in a noisy environment, pre-processing the time domain data to obtain composite images using a method that improves upon that of the prior research, then validating a trained CNN model to an accuracy of 100% on a selected candidate converter.

42 ENGINEERING↗

The northeast materials database for magnetic materials

The discovery of magnetic materials with high operating temperature ranges and optimized performance is essential for advanced applications. Current data-driven approaches are limited by the lack of accurate, comprehensive, and feature-rich databases. This study aims to address this challenge by using Large Language Models (LLMs) to create a comprehensive, experiment-based, magnetic materials database named the Northeast Materials Database (NEMAD), which consists of 67,573 magnetic materials entries (www.nemad.org). The database incorporates chemical composition, magnetic phase transition temperatures, structural details, and magnetic properties. Enabled by NEMAD, we trained machine learning models to classify materials and predict transition temperatures. Our classification model achieved an accuracy of 90% in categorizing materials as ferromagnetic (FM), antiferromagnetic (AFM), and non-magnetic (NM). The regression models predict Curie (Néel) temperature with a coefficient of determination (R 2 ) of 0.87 (0.83) and a mean absolute error (MAE) of 56K (38K). These models identified 25 (13) FM (AFM) candidates with a predicted Curie (Néel) temperature above 500K (100K) from the Materials Project. This work shows the feasibility of combining LLMs for automated data extraction and machine learning models to accelerate the discovery of magnetic materials.

Ferromagnetism↗