Results for “MIGRATION” — Page 6

DOE OSTI2025

The effect of composition on helium bubble nucleation in complex concentrated alloys

Molecular dynamics simulations of helium implantation have been performed in the MoNbTaTi complex concentrated alloy (CCA) to determine the impact of composition on helium clustering and bubble nucleation. It was found that compositions with lower interstitial atomic volume resulted in smaller cluster sizes, higher helium-to-vacancy ratios, and larger migration barriers. As the atomic volume decreases, it is harder to accommodate additional helium in the lattice and for helium to diffuse between interstitial sites, resulting in smaller cluster sizes. In particular, increasing molybdenum directly correlated with a decreased interstitial atomic volume and increased migration barrier. Niobium exhibited the opposite trend, where more niobium in the material resulted in larger interstitial atomic volumes and cluster sizes but lower migration barriers. The atomic volume was determined to be an indicator of the susceptibility of the material to helium bubble growth and could potentially be used as a metric to screen compositions for resistance to helium damage. In this way, properties such as available interstitial volumes could be used to screen materials more rapidly given the vast compositional space of CCAs.

McCarthy, Megan Jeanne [Sandia National Laboratori↗

DOE OSTI2025

Prediction of vacancy defect diffusion paths in high entropy alloys via machine learning on molecular dynamics data

Identifying the diffusion path of point defects is a critical step in understanding their evolution and the mechanisms of related phenomena. Defect diffusion occurs at small length and time scales, with impacts on material properties that may continue to evolve over ns to μs, ms, and the continuum scale (s, min, etc., and cm, m, etc.). The time scale accessible to molecular dynamics (MD) simulations is limited by small step sizes, typically in the fs range. Thus, surrogate models of MD simulations through machine learning (ML)-based algorithms are of great interest, especially for complex systems such as high entropy alloys (HEAs). In this work, dynamics governing vacancy migration in HEA were approximated with graph convolutional network (GCN) models as ansatzes for kinetic Monte Carlo (KMC) rate catalogs. Network design considered that diffusion in crystalline solids generally depends on interactions between defects and their immediate neighbor atoms. Graphs represented the vacancy surroundings, MD-generated trajectories provided training and comparison datasets, and unsupervised GCN models approximated interatomic dynamics governing vacancy migration in HEAs as ansatzes for KMC. A proof-of-concept model trained on MD data for the Fe, Ni, Cr, Co, and Cu HEA environment was used with two different neighbor interactions to assess the feasibility of training a GCN to predict vacancy defect transition rates in the HEA environment. The resulting setup rapidly generated MD-formatted synthetic trajectories based on dynamics learned from the MD training set, with a time acceleration of roughly two orders of magnitude and a similar diffusion coefficient to MD observations. Additionally, Nudged Elastic Band (NEB) calculations were performed on randomly generated FeNiCrCoCu HEA structures to determine vacancy migration barriers across nearest-neighbor sites. Transition probabilities for each jump, categorized by atomic type, were extracted from these calculations. NEB-based and GCN-based approaches led to similar outcomes.

Reimer, C↗

DOE OSTI2024

Increasing Multilayer Ceramic Capacitor Lifetime With Bipolar Voltage Cycling

Enhancing the lifetime of multilayer ceramic capacitors (MLCCs) is critical in many aerospace, naval, or electrical grid applications, where device failure could lead to catastrophic consequences. The migration of oxygen vacancies from the ceramic to the electrode interface under constant bias is known to reduce the lifetime of oxide-based MLCCs. Bias cycling presents an opportunity to enhance MLCC lifetime by reducing oxygen vacancy migration. The ideal frequency range is expected to lie between frequencies low enough to avoid self-heating but high enough to avoid interfacial defect formation. However, the impact of low-frequency bipolar voltage cycling (BVC) on MLCC degradation mechanisms has not been well studied. This work investigates the impact of periodic BVC on the degradation of MLCCs through highly accelerated lifetime testing (HALT) on X7R capacitors. HALT tests were conducted at 255 °C and 60 V using different switching frequencies: 0 (dc), 0.1, 2.5, and 10 Hz. BVC was found to improve the lifetime of MLCCs compared to dc test conditions. MLCCs tested at 10-Hz BVC showed a 311% increase in average time to failure compared to the dc case. Impedance spectroscopy shows that BVC decreases the rate of resistance degradation within MLCCs, indicating that oxygen vacancy migration to the electrodes is mitigated. The impedance spectra taken on BVC samples highlight how grain boundaries play a vital role in trapping oxygen vacancies. Periodic cycling causes oxygen vacancies to become trapped at grain boundaries, resulting in oxygen vacancies taking longer to reach the electrode interface and thus increasing MLCC lifetime. This work highlights not only how BVC can be used to increase MLCC lifetime but also how periodically cycling MLCCs could increase lifetime in extreme environments, such as at elevated temperatures and electric fields.

Chuong, Kayla Y.↗

DOE OSTI2024

Geologic Characterization of the South Georgia Rift Basin for Source Proximal CO2 Storage

The project Geologic Characterization of the South Georgia Rift Basin for Source Proximal CO2 Storage is one of 9 site characterization projects that were implemented as part of ARRA (American Recovery and Reinvestment Act). Data from this project was used to improve resolution of data in NATCARB in the area of study. Data related to this study has already been incorporated in NATCARB Atlas. The South Carolina Research Foundation and partners evaluated the feasibility of CCS in the Jurassic/ Triassic (J / TR) saline formations of the buried Mesozoic South Georgia Rift (SGR) Basin that extends from South Carolina into Georgia. The J / TR sequence, based on preliminary assessment of limited geologic and geophysical data, appears to have both the appropriate areal extent and multiple horizons to permanently and safely store CO2 The presence of several igneous rock layers within the sequence may potentially provide adequate seals to prevent upward CO2 migration into the Coastal Plain aquifer systems. Approximately 81 kilometers of 2-D seismic reflection data were collected by Bay Geophysical, Inc. to explore a portion of the SGR located in southern Georgia. The 81 kilometers were divided into two lines approximately 40.5 kilometers each, with Line 1 intersecting Georgia well GGS 3457. Line 2 intersects Line 1 at the southern portion of Line 1 to maximize the extent of coverage away from GGS-3457 (a deep well drilled in the 1980s for oil and gas exploration). This well had a set of usable logs, including gamma and neutron logs that provided promising results related to CO2 storage. Results showed sandstone with porosity values greater than 10 percent and a thickness of 120 meters. The design of the seismic shot was to extrapolate information away from the well and to better define the extent of the SGR and the necessary reservoir and caprock for a successful CO2 injection. A numerical simulation model of CO2 Injection and migration was developed based on the geology log for the GGS-3457 well. The simulation model was used to investigate the feasibility of injecting 30 million metric tons of CO2 into SGR J / TA sediments and integrity of the diabase layers as seals to prevent CO2 migration.

2-D seismic↗

DOE OSTI2024

Investigations of Relative Stability and Distribution of Mo Species in Pore Space of Mo/ZSM-5 Systems

Reduced forms of molybdenum (Mo) carbides anchoring in pore space of zeolites are considered as activated catalysts for the methane dehydroaromatization reaction. However, the atomic structure of the active sites and their stability under reactive conditions are still not clearly understood. Herein, we employ a synergistic theoretical and experimental approach to investigate this research gap. The Raman data shows the existence of the distributed binuclear and mononuclear Mo oxides in the pores that are corroborated by the hydrogen temperature-programmed reduction experiments. The binding free energies of mono- and bi-nuclear Mo carbides in straight and sinusoidal channels and their intersections as a function of temperature and anchoring sites have been determined. The catalysts with Mo coordinated at two anchoring sites are more stable compared to those located at the rings with one Al substitution. The activation energies of the catalysts migration to adjacent rings were estimated as a function of temperature and kinetics parameters of such migrations were determined. The migration is facilitated if there is one anchoring site shared between two rings. This process can contribute to an agglomeration mechanism leading to loss of catalytic activity.

Myshakin, Evgeniy↗

DOE OSTI2024

Characterization of Nitrate, Nitrite, Ammonia, and Tritium in D0220 Cores

This report describes characterization of nitrogen species (nitrate, nitrite, ammonia) in aqueous and solid phases in cores taken in borehole D0220 (well 299-E25-245) under crib 216-A-37-1 at the Hanford Site in 2022 to evaluate (a) the types of nitrogen species currently in the vadose zone and (b) the migration of nitrogen species in the vadose zone. Ammonia and tritium from PUREX decladding condensate were sporadically discharged to the crib from 1977 to 1989. During discharge operations, the estimated travel time through the vadose zone to groundwater was 2.5 to 9 months. By 2003, a characterization borehole (C4106) showed residual pore water with elevated tritium and nitrate in the first 100 ft of the vadose zone, which is likely from the crib. However, there was significantly more nitrate present at the shallowest depth (15 to 22 ft), which may indicate a different source for nitrate, such as adsorbed ammonia slowly being oxidized or nitrogen species precipitates slowly dissolving. In this study, more extensive nitrogen species characterization was done on D0220 cores (drilled in 2022) at 40- and 262-ft depths, which included (a) aqueous nitrate, nitrite, and tritium; (b) adsorbed ammonia; (c) nitrogen in carbonates (or other minerals dissolved in acidic acid); (d) nitrogen in iron oxides (or other minerals dissolved in oxalic acid); and (e) nitrogen in hard-to-extract minerals (minerals dissolved in nitric acid). High pore water nitrate (226 to 331 mg/L) at 40-ft depth measured in D0220 (2022) compared to 60 mg/L at 40-ft depth in C4106 (2003) may indicate nitrate is migrating deeper. Tritium concentrations (pore water 132,000 to 148,000 pCi/L) measured in D0220 at 40-ft depth in 2022 were considerably higher than in C4106 at 40-ft depth (160 pCi/L). Additional nitrogen species mass was present in adsorbed and precipitated phases in D0220 cores. Low adsorbed ammonia was measured at 40.1- and 261.7-ft depths. Low concentrations of carbonate associated nitrogen and iron oxide-associated nitrogen were present at 40.1 and 40.6 ft depths. Nitrogen species in solid phase extractions indicate nitrogen precipitates or aqueous nitrate or ammonia trapped in sediment microfractures that are coated by precipitates. Overall, the nitrogen species and tritium characterization at two depths in D0220 showed that additional nitrogen species were present in the vadose zone. In addition, elevated pore water nitrate in D0220 at 40 ft depth from 2022 may indicate vertical migration compared to C4106 nitrate profile from 2003. Analysis of additional depths in D0220 and spatial variability of the nitrate plume along the length of the crib (from surface electrical resistivity) would be extremely useful for this comparison.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2025

Initial Uncertainty Analysis of Carbon Tetrachloride Contamination and Remediation in the Ringold A and Lower Mud Units at the Central Plateau

The long-term effectiveness of groundwater cleanup at the Hanford Site Central Plateau depends on predictive models that can capture key uncertainties in contaminant fate and transport. Carbon tetrachloride (CCl 4 ), a persistent and toxic compound, presents particular challenges due to variability in degradation rates, uncertainty in initial plume distribution, and subsurface heterogeneity. These uncertainties directly influence plume persistence, migration pathways, and remedy performance, and thus must be systematically evaluated to support long-term remediation planning. To address these gaps, a large-scale Monte Carlo analysis was conducted using the Plateau to River (P2R) model framework. The modeling approach parameterized three primary uncertainty factors: (1) degradation rate, (2) initial plume distribution, and (3) hydraulic conductivity. Degradation was represented as a first-order process, with half-lives ranging from 70 to 700 years. Initial plume distributions were created using a geostatistical simulation method (sgsim), which generates many equally plausible versions of how contaminants might be distributed underground. From this, 100 different scenarios were mapped onto the P2R grid. Variability in hydraulic conductivity was represented in a similar way, with 100 scenarios each for the Ringold Lower Mud and Ringold A units (layers 6 and 7), based on fitted exponential variograms and conditioned to well data. In total, more than 1000 realizations were simulated to assess plume behavior under uncertainty. Results demonstrate that degradation kinetics exert the strongest control over plume persistence: Shorter half-lives produced rapid mass reduction, while longer half-lives yielded persistent plumes with limited attenuation. A nonlinear response was observed, with steep mass reductions at half-lives greater than 200 years and near-linear declines beyond this threshold, reflecting interactions between degradation and pumping. The initial plume distribution strongly influenced early transport patterns, with broader sources generating larger plume footprints, although pump-and-treat operations constrained plume migration to managed areas. By comparison, hydraulic conductivity variability in the Ringold units had only a secondary influence, modifying spreading behavior without altering the dominant migration pathways governed by source configuration and hydraulic controls. Overall, the analysis highlights that uncertainty in degradation rate and initial plume configuration are the primary drivers of variability in plume predictions, while conductivity heterogeneity plays a limited role. These findings underscore the need for improved site-specific data on degradation processes and source characterization to enhance the reliability of long-term performance assessments and to better inform remedial decision-making at the Central Plateau.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2024

Deep Learning for Fish Identification from Sonar Data (CRADA 481 Final Report)

In eastern regions of the United States, the American eel is a species of management and regulatory concern because of significant population declines, despite the species’ previous abundance in all tributaries of rivers flowing into the Atlantic Ocean. The American eel is also a candidate for listing under the U.S. Endangered Species Act. While hydropower construction and operation are only one of several factors contributing to this population decline, such a listing could impose additional regulatory challenges for a large number of hydropower projects. In this CRADA project, we improved technologies for identifying migrating eels with the goal of reducing the cost and time required for future American eel hydropower impact assessment and mitigation studies, while maintaining accuracy. We built on results from a previous FOA project (FOA# DE-FOA-0001662), led by the Electric Power Research Institute (EPRI), which developed a highly accurate, deep-learning method for identifying migrating eels from imaging sonar data. The current study aimed to further optimize this deep-learning model, originally designed for image classification, and to develop an object detection software capable of identifying fish from sonar videos in real time, enabling the detection of events like fish migrations and specific species, such as the American eel, at hydropower dams. The data conversion algorithms were packaged as software with a graphical user interface, and the software is evaluated by external collaborators. We focused on the American eel in this project and explored the transferability of the developed deep learning models to the sea lamprey, given the similar body shape and swimming behavior between the two species.

13 HYDRO ENERGY↗

DOE OSTI2025

Exponential Backoff and Its Security Implications for Safety-Critical OT Protocols over TCP/IP Networks

The convergence of Operational Technology (OT) and Information Technology (IT) networks has become increasingly prevalent with the growth of Industrial Internet of Things (IIoT) applications. This shift, while enabling enhanced automation, remote monitoring, and data sharing, also introduces new challenges related to communication latency and cybersecurity. Oftentimes, legacy OT protocols were adapted to the TCP/IP stack without an extensive review of the ramifications to their robustness, performance, or safety objectives. To further accommodate the IT/OT convergence, protocol gateways were introduced to facilitate the migration from serial protocols to TCP/IP protocol stacks within modern IT/OT infrastructure. However, they often introduce additional vulnerabilities by exposing traditionally isolated protocols to external threats. This study investigates the security and reliability implications of migrating serial protocols to TCP/IP stacks and the impact of protocol gateways, utilizing two widely used OT protocols: Modbus TCP and DNP3. Our protocol analysis finds a significant safety-critical vulnerability resulting from this migration, and our subsequent tests clearly demonstrate its presence and impact. A multi-tiered testbed, consisting of both physical and emulated components, is used to evaluate protocol performance and the effects of device-specific implementation flaws. Through this analysis of specifications and behaviors during communication interruptions, we identify critical differences in fault handling and the impact on time-sensitive data delivery. The findings highlight how reliance on lower-level IT protocols can undermine OT system resilience, and they inform the development of mitigation strategies to enhance the robustness of industrial communication networks.

DNP3↗

DOE OSTI2026

[C/N] Ages for Red Giants and Their Implications for Galactic Archaeology

Red giants undergo the first dredge-up, a mixing event that creates a connection between their surface [C/N] and their mass and age. We derive a [C/N]–age relationship for red giants calibrated on Apache Point Observatory Galactic Evolution Experiment (APOGEE) DR17 abundances and APOKASC-3 asteroseismic ages. We find that we can use [C/N] to reliably recover asteroseismic ages between 1 and 10 Gyr with average uncertainties of 1.64 Gyr. We find that [C/N] yields concordant ages, with modest offsets, for stars in different evolutionary states. We also find that the [C/N]–birth mass relationship is robust for luminous giants, and argue that this is an advantage over direct asteroseismology for these stars. We use our ages to infer Galactic birth abundance trends in [Fe/H] and [Mg/H] as a function of position in the Galactic disk. We filter out stars with kinematic or chemical properties consistent with migrators and find the number of migrators to be much lower than expected by standard radial migration prescriptions. The remaining population shows weak chemical evolution trends, on the order of 0.01 dex Gyr −1 , over the last 10 Gyr across a wide range of radii.

Roberts, John D. [The Ohio State Univ., Columbus, ↗

DOE OSTI2025

Full Device Evaluation of Metal Halide Perovskite Solar Cells in Low Earth Orbit Aboard the International Space Station

This study investigates how sustained low Earth orbit (LEO) exposure affects metal halide perovskite (MHP) thin films and perovskite solar cells (PSCs). It examines samples deployed on the Materials International Space Station Experiment (MISSE) 15 and 16 missions. Five methylammonium lead iodide (MAPI) thin films are deployed on MISSE-16, each with distinct UV filters to selectively attenuate AM0 spectral bands. While post-flight optical analysis reveals that the least UV-exposed film exhibits the highest emission and lowest non-radiative recombination rate, no clear correlation is observed among the rest, and all MAPI films maintained excellent integrity throughout the mission. MISSE-15 deployed eight PSCs with diverse structures, MHP compositions, and contact materials. Post-flight analysis reveals stable, highly emissive MHPs, but damaged contacts due to ion migration, which caused loss of electrical response. The MISSE missions demonstrate MHPs' suitability for space applications, while highlighting the need for improved interfacial layers and contact materials to enhance charge carrier mobility, prevent ion migration, and improve charge carrier extraction efficiency.

14 SOLAR ENERGY↗

DOE OSTI2026

A comprehensive first-principles study of the effects of the exchange-correlation functional and magnetism on defect and diffusion properties of the CoCrNi medium-entropy alloy

The present work is a novel, systematic study of the effect of density functional theory input parameters on the vacancy formation energy (VFE), migration barrier for diffusion, and electronic structure for each element in the CoCrNi medium-entropy alloy (MEA). In particular, the novelties include: (1) calculating the aforementioned properties of Co, Cr, or Ni, in the CoCrNi MEA using magnetic and non-magnetic states, and two versions of the generalized gradient approximation: Perdew, Burke, and Ernzerhof (PBE) and the PBE version for solids (PBEsol), and (2) a detailed comparison of 0 K activation energy to experimental creep activation energies. First-principles calculations at 0 K are performed using the Vienna ab-initio simulation package. Special quasirandom structures (SQS) and Widom-type substitution are employed. For each element, Co, Cr, or Ni, non-magnetic calculations result in a higher VFE and larger range of calculated values for the configurations studied. The averaged migration barrier is the highest for Co in the CoCrNi for three of four sets of calculation parameters in the configurations studied. Finally, the results indicate that the average 0 K activation energy for diffusion makes up 70–80% of the experimental creep activation energy, depending on the exchange-correlation functional employed.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Validation of the ERO2.0 code using W7-X and JET experiments and predictions for ITER operation

Abstract The paper provides an overview of recent modelling of global material erosion and deposition in the fusion devices Wendelstein 7-X (W7-X), JET and ITER using the Monte-Carlo code ERO2.0. For validating the modelling tool in a three-dimensional environment, W7-X simulations are performed to describe carbon erosion from the graphite test divertor units, which were equipped in operational phase OP 1.2 and analysed post-mortem. Synthetic spectroscopy of carbon line emission is compared with experimental results from the divertor spectrometer measurement system, showing a good agreement in the e-folding lengths in the radial intensity profiles of carbon. In the case of metallic wall materials, earlier modelling of the Be/W environment in JET and ITER is revisited and extended with an updated set of sputtering and reflection data, as well as including the mixing model for describing the Be/W dynamics in the divertor. Motivated by recent H/D/T isotope experiments in JET, limited and diverted configuration pulses are modelled, showing the expected trend of both Be and W erosion increasing with isotope mass. For the JET diverted configuration pulses, it is shown that Be migrates predominantly to the upper part of the inner divertor where it initially leads to strong W erosion. With longer exposure time, the growth of a Be deposited layer leads to a reduction of W erosion in that region. A similar trend is observed in simulations of the ITER baseline Q = 10 scenario, however with a more symmetric Be migration pattern leading to deposition also on the outer divertor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Microstructure-Based Degradation Modeling in Solid Oxide Cells

The presentation is intended for a graduate class in Georgia Southern University taught by Dr. Hayri Sezer on electrochemistry, SOFC and Li ion battery. It summaries NETL's recent works on microstructure-based modeling of degradation in solid oxide cells. It covers the degradation mechanism of Ni coarsening, Ni migration and the build-up of oxygen partial pressure inside the electrolyte. The phase-field models developed for Ni coarsening and migration are reviewed and the results are compared to experiments in literature to examine different diffusion mechanisms (Ni self-diffusion, Ni(OH)2 diffusion in pores and Ni(OH) surface diffusion) and driving forces (Ni(OH)x concentration and Ni-YSZ contact angle change). It also demonstrates the developed model for oxygen partial pressure and shows that the predicted oxygen partial pressure can explain the Ni oxidation in the hydrogen electrode and crack growth on the electrolyte-oxygen electrode interface under certain operating conditions.

nickel coarsening↗

DOE OSTI2025

Planetesimal Impact Vapor Plumes and Nebular Shocks Form Chondritic Mixtures

The origin of chondrules and the chondritic sedimentary rocks that dominate the meteoritic record is a long-standing problem in planetary science. Here, we develop a physical model for the formation of chondritic mixtures as an outcome of vaporizing collisions between planetesimals that were dynamically excited by the growth and migration of planets. We present calculations of nebular shock waves generated by impact vapor plumes and focus on aspects of the plume interaction with the nebular gas and dust that have been neglected in previous studies of impact ejecta. We find that when water dominates the vapor, the plumes are relatively cool. However, the plume expansion is supersonic and can drive strong shock waves in the dusty nebular gas. Portions of these nebular shock fronts initially melt nebular dust, forming chondrules that are coupled to the moving front. As the shock front expands and cools, the chondrules solidify, while the shock front entrains additional dust. Eventually, the plume expansion stalls and then hydrodynamically collapses, turbulently mixing variably processed dust and size-sorted chondrules. For probable impact parameters and nebular conditions during giant planet growth and migration, the impact-generated mixtures have characteristics that span the range observed in chondritic meteorites, providing an environment for rapid formation of chondritic assemblages after chondrule formation. Our impact vapor and nebular shocks model, named the IVANS model, links chondrule formation to the overall context of planet formation and provides a framework for interpreting the detailed chronological and geochemical record contained in chondritic meteorites.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

DOE OSTI2025

Shape‐Evolving Structured Liquids

Migration, division, and reconfiguration - functions essential to living systems - are driven by active processes. Developing synthetic mimics is an outstanding challenge. Lipid bilayers that bound natural systems are locally deformed by active species, e.g., microtubules, but the resulting non-equilibrium shapes relax when active species motion ceases, and the shape changes lack immediate control. A fully synthetic system is described, driven by active particles encapsulated by a reconfigurable nanoparticle-surfactant membrane that undergoes shape fluctuations reminiscent of living cells. These shape changes are preserved after particle activity stops. Surfactant concentration tunes the interfacial tension over three orders of magnitude, making on-demand shape evolution possible. Directional migration, division, and reconfiguration across multiple scales are possible, leading to a new class of biomimetic, reconfigurable, and responsive materials, paving the way for autonomous synthetic machines.

Kim, Paul Y↗