Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “RICHLAND”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Interfacial Hydrophilicity Controls Mineral Transformation Outcomes for Enstatite and Amorphous MgSiO 3

Subsurface injection of carbon dioxide (CO 2 ) into mafic-ultramafic rocks for permanent storage via mineralization is being studied to reduce emissions. Here, we investigated the carbonation products of enstatite (MgSiO 3 ) to assess its efficiency in sequestering CO 2 for safe and permanent storage as carbonate minerals. This was accomplished by conducting variable temperature carbonation reactions with samples of differing crystallinities and surface chemistries. Reaction progress was monitored utilizing in situ X-ray diffraction, and the presence of carbonate products was confirmed using additional techniques, such as thermogravimetric analysis coupled with mass spectrometry and scanning electron microscopy with energy dispersive spectrometry. Our results show that crystalline enstatite produces small amounts of the anhydrous form of MgCO 3 (magnesite), while amorphous MgSiO 3 , which was used to simulate mafic glass, more readily converts to the hydrated/hydroxylated hydromagnesite [Mg 5 (CO 3 ) 4 (OH) 2 ·4H 2 O]. These results, supplemented with dynamic vapor sorption experiments, suggest that surface properties play a significant role in the pathway and degree of carbonation. These developments concerning the reactivity of CO 2 with reactive mafic phases will help further our understanding of the reactivity of these mafic-ultramafic minerals with implications for permanent carbon storage and other subsurface engineering scenarios involving reactive reservoirs.

basalt↗
DOE OSTI2025

Bench-Scale Apparatus for Rapid, Simultaneous, Comprehensive Vapor–Liquid Equilibrium and Kinetic Property Measurements to Advance CO 2 Capture Solvent Technology

We present here, the design, fabrication, and testing of a new "pressure, volume, temperature (PVT)" apparatus that enables rapid standardized testing for viscosity, vapor–liquid equilibria, and kinetics for water-lean carbon capture solvents. This unit is the first of its kind where equilibrium data are collected during operation as a PTx (pressure–temperature–composition) cell while kinetic data are collected simultaneously with an internal mini wetted-wall contactor (WWC) using controlled adjustments of CO 2 injections to allow for measurements of gas flux (in and out of the liquid). Additionally, in situ measurements of viscosity data are also continuously collected while a solvent is in circulation during gas absorption. This cell empowers comprehensive testing of critical CO 2 capture solvent properties in a single measurement, ensuring all data are collected at the same temperature, pressure, and CO 2 loading. This apparatus also expedites screening of materials since less than 50 mL of sample is needed as compared to 2–3 L needed to get similar data from a conventional WWC. In conclusion, we describe here the methodology of data collected on this new PVT cell, nicknamed “Gary” in honor of Professor Gary Rochelle, for multiple water-lean amine solvents, which we compare to data collected from conventional instrumentation for amine testing.

Zheng, Richard [STARS Technology Corporation, Rich↗
DOE OSTI2025

Nitrilotriacetamides: Nitric Acid and Water Extraction and Physicochemical Properties

Nitrilotriacetamides are an emerging class of extractants that have been identified as promising candidates for a variety of separation schemes. Systems incorporating these ligands are prone to third-phase formation at high loading, but this can be mitigated by increasing the alkyl chain length or adding a phase modifier (e.g., 1-octanol). In this work, the impacts of these modifications on HNO 3 extraction, H 2 O extraction, and physiochemical properties (interfacial tension, viscosity, and density) are examined. The studies show that ligand alkyl chain length has very little impact on HNO 3 and H 2 O extraction, whereas the addition of 10% octanol significantly increases extraction of both HNO 3 and H 2 O. Viscosity and interfacial ten-sion are impacted by HNO 3 concentration, alkyl chain length, and the presence of phase modifier but remain in a range that should be compatible with commonly used solvent contacting equipment. Additional studies of physiochemical properties at variable ligand concentration suggest that a change in aggregation occurs above approximately 10-50 mM ligand. Thermodynamic modeling of distribution ratios and spectroscopic studies aimed at understanding speciation in the presence of HNO 3 are also reported. Further, these experiments suggest that the first equivalent of extracted HNO 3 interacts with both the central amine and carbonyl groups of the ligand. The spectroscopic changes at higher HNO 3 concentrations are relatively minor, suggesting additional HNO 3 is associated in an outer-sphere manner

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

Solvation Structure of 237 Np 4+ in a Noncomplexing Environment

Here, the solvation structure of an Np 4+ ion in an aqueous, noncomplexing and nonoxidizing environment of trifluoromethanesulfonic (triflic) acid was investigated with X-ray absorption spectroscopy (XAS) combined with ab initio molecular dynamics (AIMD) and time-dependent density functional theory (TDDFT) calculations. Np L III -edge X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) data were collected for Np 4+ in 1, 3, and 7 M triflic acid using a laboratory-scale spectrometer and separately at a synchrotron facility, producing data sets in excellent agreement. TDDFT calculations revealed a weak pre-edge feature not previously reported for Np L III -edge XANES. AIMD modeling results showed differences in the hydration shell of the Np 4+ ion at different concentrations of triflic acid; these results are supported by the experiment. EXAFS fit models to the experiment resulted in similar coordination of Np 4+ in noncomplexing aqueous media as reported in the literature for 1 M perchloric acid but, together with calculations, revealed more than one distance between Np and O atoms in 7 M triflic acid. These results imply monodentate coordination with sulfonate groups in 7 M triflic acid and suggest the possibility of proto-neptunyl species in relatively low-concentration Np 4+ acid solutions.

Boglaienko, Daria V. [Pacific Northwest National L↗
DOE OSTI2025

Perspective on Many-Body Methods for Molecular Polaritonic Systems

Recent advances in strong light–matter interactions have revealed a wealth of new physical phenomena in molecules embedded in optical cavities, including modified chemical reactivity, altered excitation spectra, and novel quantum correlations. To describe these effects from first-principles, the field of ab initio quantum electrodynamics (QED) has emerged as a compelling extension of quantum chemistry that treats electronic and photonic degrees of freedom on equal footing. In this Perspective, we review the growing landscape of many-body QED methods, including Hartree–Fock, density functional theory (QEDFT), time-dependent DFT (QED-TDDFT), configuration interaction (QED-CI), complete active space (QED-CASSCF), coupled cluster (QED-CC), quantum Monte Carlo (QED-QMC), and density matrix renormalization group (QED-DMRG), highlighting recent developments and implementations. We further explore real-time methods, gradient and Hessian formalisms, and the integration of nonadiabatic nuclear dynamics. Applications range from benchmark simulations of polaritonic chemistry to quantum simulations on emerging quantum hardware. We conclude by outlining future directions for theory development and interdisciplinary efforts at the interface of quantum chemistry, condensed matter, and quantum optics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Quantum Electrodynamics Coupled-Cluster at Scale: High-Performance Implementation for Complex Systems

Coupled-cluster theory (CC) is a highly accurate and versatile method for simulating complex interactions within quantum systems. The extension of CC theory to model mixed electron-photon processes with quantum electrodynamics (QED) has improved our capability to predict cavity-modified chemistry, a field where photons are used as cost-effective and eco-friendly alternatives to catalyze/inhibit chemical reactions. However, calculations with CC methods, even without incorporating QED effects, are often prohibitively expensive. Simulations of larger systems require scalable infrastructures that exist for traditional CC methods but not for QED-CC methods. As such, we present a GPU-enabled, high-performance, open-source implementation of the quantum electrodynamics coupled-cluster method with single and double excitations (QED-CCSD) within the ExaChem quantum chemistry software package. ExaChem relies on the Tensor Algebra for Many-body Methods (TAMM) infrastructure: a parallel heterogeneous tensor library designed to achieve scalable performance on modern heterogeneous supercomputing platforms. Furthermore, we discuss theoretical foundations, algorithmic details, and numerical benchmarks to showcase the larger systems that ExaChem can simulate and how the integration of photonic degrees-of-freedom alters their ground-state properties.

Basis sets↗
DOE OSTI2026

Strategic Melamine Coating on Lithium Metal for Li 3 N-Rich Solid Electrolyte Interphase and Improved Battery Cycling Stability

The practical applications of lithium (Li) metal batteries (LMBs) are limited by challenges such as dendrite formation and unstable solid electrolyte interphase (SEI), especially at higher C-rates. Here, this study introduces melamine-coated Li metal anodes (LMAs), forming a Li 3 N-rich SEI layer that improves ionic conductivity and mechanical stability. The optimized melamine-coated LMA demonstrated uniform coverage resulting in denser Li deposition, nearly doubled cycle life (~148 cycles at 0.5 C, 1C = 4.1 mA cm -2 ), compared to Bare-Li. These findings emphasize that coating materials-induced beneficial SEI components could lead to improvement of LMB performance.

Batteries↗
DOE OSTI2025

Demonstration of an Automated System for Vertical Profiles of Volatile Organic Compounds

Volatile organic compounds (VOCs) play important roles throughout the atmosphere, many of which are altitude dependent. This highlights the need for easily deployable devices to sample VOCs across different atmospheric layers. To address this, we present the design and initial application of a Time Resolved Automated Volatile organIc compounds Sampling system (TRAVIS). VOCs are collected on sorbent tubes, which are subsequently analyzed by a thermal desorption gas chromatography mass spectrometry pipeline. TRAVIS leverages a piezoelectric pump with an integrated pressure sensor for precise (0.1% flow rate relative standard deviation) and accurate (−3 ± 2% error in VOC quantitation) measurements. Via deployment on a tethered balloon system over an agricultural area, TRAVIS is used to show consistent vertically resolved VOC concentrations in a well-mixed (i.e., turbulent) atmosphere (e.g., 5% relative standard deviation for isoprene) and vertically dependent concentrations for a stratified atmosphere (e.g., prior to boundary layer development). Furthermore, we also show VOC information from an intermittent plume via both targeted and untargeted analysis, highlighting future applications for spurious events in agriculture, air quality monitoring, and environmental impact. Overall, the development of TRAVIS represents a lightweight, accurate, sensitive, and precise VOC sampling module for the scientific community.

Aerosols↗
DOE OSTI2025

Secondary Organic Aerosol from OH Oxidation of Acyclic Terpenes Is More Viscous and Less Volatile than That of Their Cyclic Analogs

Biogenic volatile organic compounds (BVOCs), a dominant source of secondary organic aerosol (SOA) globally, exhibit emission rates and composition that are plant species-specific and vary with environmental stressors. A common outcome of plant stress is increased emissions of acyclic terpenes. The paucity of information about acyclic terpene SOA chemistry contributes to uncertainties in predictions of SOA global loadings and impacts on Earth’s radiative budget, particularly in a changing climate where acyclic terpene emissions could become more prominent. This study compared properties of SOA derived from OH oxidation of acyclic and cyclic monoterpenes (ß-ocimene, a-pinene) and sesquiterpenes (ß-farnesene, ß-caryophyllene). Single particle mass spectrometry was used for assessing shape, density, and evaporation kinetics of size-selected SOA particles, and nanospray desorption electrospray ionization high-resolution mass spectrometry (nano-DESI-HRMS) was used to measure molecular composition of SOA. Acyclic terpene SOA exhibited higher viscosity and lower volatility compared to cyclic terpene SOA, and had a greater volume fraction remaining (VFR) after ~24 hours of evaporation - approximately 1.3-1.6 times higher VFR than that of cyclic terpene SOA. Additionally, HRMS analysis revealed greater chemical diversity and higher fractions of extremely low volatility compounds (56-62% ELVOC/LVOC) in acyclic terpene SOA compared to cyclic counterparts (25-37% ELVOC/LVOC). Our findings highlight the potential importance of accounting for acyclic terpene aerosol chemistry under conditions of plant stress to improve predictions of SOA loadings and impacts.

SOA physical properties↗
DOE OSTI2025

Developing Fluorescence-Based Sensors to Support Rare Earth Element Separation

Rare earth elements (REEs) are essential to most renewable energy technologies. Unfortunately, as we transition to sustainable energy production, the demand for REEs is rapidly growing well beyond current rates of production. As a result, novel means of efficient, scalable, and easily adaptable methods for processing primary and recycle feedstocks are needed. Development and integration of sensors for highly selective in-line monitoring can support more efficient design and testing of such novel separation processes, as well as more cost-effective deployment of those separation flowsheets. Work here will explore the application of fluorescence spectroscopy, a highly sensitive and selective technique, to quantify multiple lanthanides in complex mixtures including known interferents or quenching agents. Results include identification of the optimal excitation wavelength and the limit of detection of various rare earth elements as well as the performance of data-science-based quantification approaches in streams where “unknowns” are present. Overall, the data science tools in conjunction with optical sensor data were able to quantify analytes in the presence of other lanthanides which can be anticipated in the actual industrial stream. Here we include characterization of lanthanides in a microfluidic device similar to those used in new process development. This study demonstrates the capability of utilizing fluorescence spectroscopy to quantify analytes in a complicated solution matrix, suggesting this is a successful approach for in-line monitoring to optimize the separation efficiency in an industrial stream.

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

Tuning Catalytic Reactivity via Wetting Control through Oxygen Vacancies: Ru Clusters on Anatase TiO 2 and CeO 2 Supports

The shape of supported metal particles regulates their catalytic reactivity and is determined by the degree of wetting between the metal particle and the support surface. Flattened particles that wet support surfaces were reported in various catalytic systems, particularly in the subnanometer size regime. Such consequential metal–support wetting phenomena are poorly understood, and methods to study them on powder catalysts under realistic conditions are lacking. Here, we investigate the size-dependent wetting behaviors of Ru particles on two reducible-oxide supports, anatase TiO 2 (TiO 2 -A) and CeO 2 , under reducing catalytic conditions. X-ray absorption spectroscopy (XAS), low-energy ion scattering (LEIS), and density functional theory (DFT) are combined to determine the shape of Ru particles. Ru particles remain three-dimensional without wetting the TiO 2 -A support within the coverage range studied (0.06–0.98 Ru nm –2 ). In contrast, at low coverages (<0.25 Ru nm –2 ), Ru wets the CeO 2 support to form flat, disordered structures. The higher wettability of CeO 2 than TiO 2 -A is attributed to oxygen vacancies in the near-surface region. The shape difference between small Ru particles or clusters on the two supports leads to drastically contrasting catalytic reactivities in polyolefin hydrogenolysis, despite similar diameters. As a result, this work highlights the implications of metal–support wetting, or cluster shape, on catalytic behaviors of small metal clusters, while establishing the foundation for future systematic studies of such a phenomenon in realistic systems, by delivering a multitechnique methodology and revealing governing fundamental principles.

42 ENGINEERING↗
DOE OSTI2025

Decomposing a Compound Flood Event in an Urban Pacific Northwest Estuary: Primary Drivers and Projections for the Future

The Duwamish River Estuary (DRE) of Washington is prone to compound flooding during atmospheric river (AR) events. The processes contributing to such flooding (coastal and fluvial) have remained opaque to municipalities that are increasingly impacted. Here, we conduct a suite of coupled atmosphere-hydrology-ocean model simulations with varying forcing combinations (tide, surge, and/or river discharge) to identify the primary drivers of compound flooding during a recent AR event. We also test year 2,100 climate forcing to project how flooding and drivers may change in the future for the same event. We identify a clear distinction between dynamics in the downstream, engineered portion of the DRE compared to the upstream, “natural” river. Downstream, tides dominate water levels but contributions from storm surge and nonlinear tide-surge interaction elevate tide-only high waters from no flooding to major flooding. Upstream, total water levels during the event are ∼6 cm higher than downstream due to an increasing influence of river discharge over surge and tides. Notably, nonlinear surge-river and tide-river interactions act to reduce upstream water levels up to 50% compared to estimates which linearly sum tides, surge, and river, likely reducing flood vulnerability. Under two future climate scenarios: one with only sea level rise (SLR) and another with SLR plus atmospheric warming, we find little change in mechanism contributions to water levels. Expanded flooding in both cases is largely due to SLR, as a ∼50% increase to river discharge under the warming scenario has no impact downstream and marginally increases (∼3 cm) water level upstream.

atmospheric river↗
DOE OSTI2025

Composition and Plume Gas Interaction Control Iron Fractional Solubility More Than Particle Size in Volcanic Ash: Implications for Fertilization of the North Atlantic

Deposition of volcanic ash is thought to impact marine biogeochemical cycling by adding soluble iron (Fe) to the surface ocean. The magnitude of this input is a function of the amount of ash deposited, the total Fe content in the ash, and ash‐derived Fe's fractional solubility. However, the relative importance of chemical composition, acidic processing by the volcanic plume, and ash particle size in determining solubility is unclear. We paired an aerosol leach meant to provide an upper limit for fractional Fe solubility with chemical analyses of ash from the Cumbre Vieja (CV) and La Soufrière eruptions, which both impacted the North Atlantic in 2021. Fe in the ash samples is <6% soluble, but Fe fractional solubility in CV ash is approximately triple that of La Soufrière ash. Compared to La Soufrière, a larger proportion of the Fe in CV ash is in silicate rather than oxide minerals, which release more soluble Fe. Elevated levels of surficial fluorine (F) also suggest that CV ash was subjected to a more fluorine‐rich eruption plume and underwent more acidic processing. Particle size does not appear to be a primary control on Fe release. We estimate that the CV eruption had a much larger impact on dissolved Fe (DFe) concentration in the surface ocean than the La Soufrière eruption because of differences in soluble Fe content and particle deposition velocity. These differences may help explain why some eruptions elicit a biological response in the ocean while others do not.

Elliott, Hope E. [Univ. of Miami, FL (United State↗
DOE OSTI2025

Implementation and Evaluation of Emission‐Driven Land‐Atmosphere Coupled Simulation in E3SMv2.1

Emissions-driven (prognostic CO 2 ) simulations are essential for representing two-way carbon-climate feedback in Earth System Models. We present an emissions-driven land–atmosphere coupled biogeochemistry (BGC) configuration (BGCLNDATM_progCO2) in version 2.1 of the Energy Exascale Earth System Model (E3SMv2.1). This is the first E3SM configuration that performs land-atmosphere emission-hindcasts. Here, we document its implementation, evaluate the model's performance against observations and other models, and propose a structured evaluation protocol for such emissions-driven simulations. We conducted transient historical simulations (1850–2014) with BGCLNDATM_progCO2 and compare them to reference simulations—a land-atmosphere coupled simulation without BGC and a standalone land simulation with BGC, both using prescribed CO 2 concentrations—and to observations. BGCLNDATM_progCO2 overestimates atmospheric CO 2 concentrations by 11–23 ppm yet stays within the 40-ppm spread CMIP6 emission-driven models and retains physical climate properties comparable to the reference runs. The CO 2 biases are partly attributed to underrepresented oceanic CO 2 uptake and inadequate representations of some terrestrial processes. In general, introducing prognostic CO 2 did not change physical climate metrics at the global scale but had larger regional effects, particularly over land where spatially heterogeneous CO 2 and prognostic leaf area index influenced surface energy balance. Finally, we propose a general evaluation protocol including spin-up assessment, atmospheric CO 2 benchmarking, physical climate evaluation, and land biogeochemical analysis to support scientific rigor and facilitate inter-model comparisons. The new configuration lays the groundwork for future enhancements, including improved terrestrial biogeochemical processes, integrated marine biogeochemistry, and additional human–Earth system interactions. These developments advance E3SM toward fully coupled emissions-driven simulations, enabling more accurate carbon–climate feedback projections and informing mitigation policy by providing physically consistent carbon-budget metrics for mitigation scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Evaluation of E3SM Simulated Aerosols and Aerosol‐Cloud Interactions Across GCM and Convection‐Permitting Scales

This paper introduces an Earth system modeling testbed for predicting aerosols and aerosol‐cloud interactions (ACIs) at convection‐permitting scales. Using the Energy Exascale Earth System Model (E3SM) version 2 with a four‐mode Modal Aerosol Module, we conduct simulations at 3.25 km resolution on a regionally refined mesh (RRM) across four regions with distinct aerosol and cloud regimes. Results are compared with the standard 100 km E3SM configuration and evaluated against satellite, aircraft, and ground‐based observations. We find that increasing model resolution improves heavy precipitation simulation but amplifies positive bias in light drizzle at coarse resolution. These resolution‐induced changes affect cloud and aerosol properties to varying degrees across regions. Generally, cloud cover and liquid water path (LWP) show better agreement with satellite retrievals at 3.25 km, though surface‐based comparisons suggest otherwise. Aerosol composition remains poorly represented at both resolutions. The RRM increases Aitken mode aerosol number concentrations via enhanced new particle formation. However, accumulation mode aerosols are decreased at higher resolution as aerosol removals become more efficient. This partially contributes to fewer cloud condensation nuclei (CCN) and lower cloud droplet number concentrations (N d ), which produces larger model biases in some scenarios. These findings suggest that solely increasing horizontal resolution to kilometer scales is insufficient to broadly improve aerosol and cloud predictions without concurrent advancements in physical and chemical process representations. Nonetheless, the RRM moderately improves key ACI relationships such as CCN‐N d correlation, reflecting enhanced aerosol activation representation. The LWP‐N d relationship is also better captured by RRM, suggesting a better characterization of LWP adjustment.

Huang, Meng [Pacific Northwest National Laboratory↗
DOE OSTI2026

Impact of solvent forces and broken symmetry on the assembly of designed proteins at a liquid-solid interface

The era of protein design has enabled the creation of hybrid protein-inorganic interfaces, leading to both surface-directed self-assembly of de novo protein architectures and protein-directed formation of inorganic materials. However, the resulting patterns of protein assembly are often unexpected, implying that essential interactions are not accounted for in current design platforms. Here, we use high-speed atomic force microscopy (AFM) analyzed through machine learning to follow the assembly of protein nanorods in aqueous electrolytes on two types of mica exhibiting disparate symmetry elements, which are imprinted on the overlying hydration structure. Using Monte Carlo simulations, we reproduce the observed phases and show that an observed smectic phase, previously thought to be unstable for non-interacting rods in two dimensions, emerges when crystal symmetry introduces a directional bias. The findings demonstrate the importance of incorporating solvent forces as modulated by the hydration structure inherent to interfacial systems when designing protein assemblies at liquid-crystal interfaces. Coupling physics-based simulations that can account for these factors to de novo protein design algorithms can lead to improved design platforms for bio-inspired, hybrid materials.

bioinspired materials↗
DOE OSTI2026

Continental-scale integration of soil metagenomes and organic matter chemistry reveals ubiquitous microbial capacity for chemically-recalcitrant carbon decomposition

Soil organic matter (SOM) decomposition by microorganisms is a major uncertainty in predicting terrestrial carbon–atmosphere feedbacks, partly because we lack understanding of the microbial diversity involved in depolymerizing different carbon pools across environmental gradients. We address this gap using a continental-scale dataset pairing shotgun metagenomes with high-resolution SOM chemistry, assembling 0.76 Tbp of prokaryotic MAGs (828 genomes) and identifying 66,727 SOM molecules from 47 standardized U.S. soil cores selected using respiration rates from 106 soils. Integrating these datasets reveals widespread microbial potential for depolymerizing chemically-recalcitrant SOM previously considered stable. We uncover complementary metabolic specialization between genera affiliated with two abundant bacterial orders, Rhizobiales and Chthoniobacterales, and an archaeal order, Nitrososphaerales. This metabolic partitioning is consistent across soil depths and activity levels, suggesting coordinated decomposition of complex SOM through distinct but complementary biochemical strategies. The metabolic potential for depolymerization of chemically-recalcitrant compounds is supported by the abundance of these molecules across the soils, as indicated by Fourier-Transform Ion Cyclotron Resonance Mass Spectrometry (FTICR-MS), and by flux balance analysis of metabolic models. Our results show that a substantial portion of ostensibly stable SOM remains vulnerable to microbial decomposition, a mechanism not captured in current Earth System Models.

Song, Young C. [Pacific Northwest National Laborat↗
DOE OSTI2025

Molecular dynamics study of grain boundaries as defect sinks under irradiation in LiAlO 2 and LiAl 5 O 8

Lithium aluminate ceramics, LiAlO 2 and LiAl 5 O 8 , show promise in nuclear environments due to their excellent radiation tolerance. Molecular dynamics simulations investigate grain boundaries (GB) and their role in defect evolution. Results reveal that GBs act as efficient defect sinks, with Li and Al atoms exhibiting distinct behaviors during displacement cascades. Tritium migration in LiAlO 2 is also studied, showing rapid diffusion and stable configurations with oxygen, corroborated by ab initio simulations from the literature. The calculated tritium diffusion coefficient of 1.33 × 10 - ¹⁴ m²/s aligns with the literature, validating the model. LiAl 5 O 8 demonstrates superior defect healing compared to LiAlO 2 , attributed to enhanced atomic transfer between grains and GBs. These findings reveal key insights into defect dynamics, providing essential insights for their application in tritium-producing burnable absorber rods (TPBARs).

36 MATERIALS SCIENCE↗