Results for “time-dependent” — Page 13

DOE OSTI2025

Composite Behavior of Nanopore Array Large Memristors

Synthetic nanopores were recently demonstrated with memristive and nonlinear voltage-current behaviors, akin to ion channels in a cell membrane. Such ionic devices are considered a promising candidate for the development of brain-inspired neuromorphic computing techniques. In this work, we show the composite behavior of nanopore-array large memristors, formed with different membrane materials, pore sizes, electrolytes, and device arrangements. Anodic aluminum oxide (AAO) membranes with 5 nm and 20 nm diameter pores and track-etched polycarbonate (PCTE) membranes with 10 nm diameter pores are tested and shown to demonstrate memristive and nonlinear behaviors with approximately 107–1010 pores in parallel when electrolyte concentration across the membranes is asymmetric. Ion diffusion through the large number of channels induces time-dependent electrolyte asymmetry that drives the system through different memristive states. The behaviors of series composite memristors with different configurations are also presented. In addition to helping understand fluidic devices and circuits for neuromorphic computing, the results also shed light on the development of field-assisted ion-selection-membrane filtration techniques as well as the investigations of large neurons and giant synapses. Further work is needed to de-embed parasitic components of the measurement setup to obtain intrinsic large memristor properties.

Chemistry↗

DOE OSTI2026

NASA’s Pandora SmallSat Mission: Simulating the Impact of Stellar Photospheric Heterogeneity and Its Correction

Stellar photospheric heterogeneity is a dominant astrophysical systematic impacting exoplanet transmission spectroscopy. NASA’s Pandora SmallSat Mission is designed to address this challenge through contemporaneous visible-band photometry and near-infrared spectroscopy of exoplanet host stars. Here, we present an end-to-end simulation study quantifying Pandora’s ability to infer stellar photospheric properties and correct stellar contamination using out-of-transit observations. We construct eight representative stellar activity scenarios and generate 160 simulated Pandora datasets, incorporating time-dependent stellar spectra, instrument response, and noise. Given accurate models, Bayesian retrievals of joint visible photometry (0.4–0.7 μm) and near-infrared spectroscopy (0.9–1.6 μm, R ≈ 120) recover photospheric temperatures with typical uncertainties of ≈30 K, with no significant bias. Models with two spectral components (i.e., a quiescent photosphere and spots) are strongly favored in 95% of cases; one-component models are preferred when true spot filling factors fall below a detection threshold of ≈0.3%. We propagate the true and inferred stellar parameters to compute true, inferred, and residual contamination signals under physically motivated spot geometries. For simple spot distributions, contamination signals of 10 2 –10 3 ppm are reduced to ≲10 ppm—well below Pandora’s expected transmission spectroscopy precision (30–100 ppm). For more complex spot distributions, geometric degeneracies limit deterministic corrections, leaving residual contamination at the 10 3 ppm level that must be mitigated using additional constraints, such as spot-crossing events and joint stellar–planetary retrievals of transmission spectra. These results define regimes in which stellar contamination can be corrected from stellar observations alone and show how Pandora stellar observations can identify cases where additional information is required.

Astronomy and AstroPhysics↗

DOE OSTI2025

Thermal Radiation Transport with Tensor Trains

We present a novel tensor network algorithm to solve the time-dependent, gray thermal radiation transport equation. The method invokes a tensor train (TT) decomposition for the specific intensity. The efficiency of this approach is dictated by the rank of the decomposition. When the solution is “low rank,” the memory footprint of the specific intensity solution vector may be significantly compressed. The algorithm, following a step-then-truncate approach of a traditional discrete ordinates method, operates directly on the compressed state vector, thereby enabling large speedups for low-rank solutions. To achieve these speedups, we rely on a recently developed rounding approach based on the Gram-SVD. We detail how familiar S N algorithms for (gray) thermal transport can be mapped to this TT framework and present several numerical examples testing both the optically thick and thin regimes. The TT framework finds low-rank structure and supplies up to ≃60× speedups and ≃1000× compressions for problems demanding large angle counts, thereby enabling previously intractable SN calculations and supplying a promising avenue to mitigate ray effects.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

An Integral-based Technique to Accelerate the Monte Carlo Radiative Transfer Computation for Supernovae

We present an integral-based technique (IBT) algorithm to accelerate supernova (SN) radiative transfer calculations. The algorithm utilizes “integral packets,” which are calculated by the path integral of the Monte Carlo (MC) energy packets, to synthesize the observed spectropolarimetric signal at a given viewing direction in a 3D time-dependent radiative transfer program. Compared to the event-based technique (EBT) proposed by M. Bulla et al., our algorithm significantly reduces the computation time and increases the MC signal-to-noise ratio (S/N). Using a 1D spherical symmetric Type Ia SN ejecta model DDC10 and its derived 3D model, the IBT algorithm has successfully passed the verification of spherical symmetry and cross comparison on a 3D SN model with the direct-counting technique and EBT. Notably, with our algorithm implemented in the 3D MC radiative transfer code SEDONA, the computation time is faster than EBT by a factor of 10−30, and the S/N is better by a factor of 1.5−3, with the same number of MC quanta.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

A Search for Millimeter-bright Blazars as Astrophysical Neutrino Sources

The powerful jets of blazars have been historically considered as likely sites of high-energy cosmic-ray acceleration. However, the particulars of the launched jet and the locations of leptonic and hadronic jet loading remain unclear. In the case when leptonic and hadronic particle injection occur jointly, a temporal correlation between synchrotron radiation and neutrino production is expected. We use a first catalog of millimeter wavelength (95–225 GHz) blazar light curves from the Atacama Cosmology Telescope for a time-dependent correlation with 12 yr of muon neutrino events from the IceCube South Pole Neutrino Observatory. Such millimeter emission traces activity of the bright jet base, which is often self-absorbed at lower frequencies and potentially gamma-ray opaque. We perform an analysis of the population, as well as analyses of individual, selected sources. We do not observe a significant signal from the stacked population. TXS 0506+056 is found as the most significant, individual source, though this detection is not globally significant in our analysis of selected active galactic nuclei. Our results suggest that the majority of millimeter-bright blazars are neutrino dim. In general, it is possible that many blazars have lighter, leptonic jets, or that only selected blazars provide exceptional conditions for neutrino production.

Abbasi, R. [Loyola University, Chicago, IL (United↗

DOE OSTI2026

Exploring the Origins of Anti-Ambipolarity in BBL Polymer: Links to Redox Chemistry, Electronic Structure, and Structural Dynamics

We examine the intrinsic physical-chemical properties of the conjugated ladder-type polymer poly(benzimidazobenzophenanthroline) (BBL) in response to electron transfer. We aim at explaining the origin of the anti-ambipolar behavior behind the observed BBL nonlinear response associated with specific device architectures. To elucidate this point, we use theory and computation based on first principles, including density functional theory optimizations, ab initio molecular dynamics, time-dependent DFT, and Marcus-theory analysis. Our results reveal that this redox response is not simply monotonic but follows an alternating odd/even pattern in which gap narrowing and reopening occur sequentially before near-gapless behavior emerges at high charging. Converging theoretical evidence in this work demonstrates that bell shaped conductivity in BBL originates in its fundamental electronic structure and supramolecular organization.

FOS: Physical sciences↗

DOE OSTI2025

MOSAiC studies of long-lasting mixed-phase cloud events and analysis of the liquid-phase properties of Arctic clouds

Vertically resolved observations of the temporal evolution of mixed-phase clouds (MPCs) were performed over the central Arctic during the MOSAiC (Multidisciplinary drifting Observatory for the Study of Arctic Climate) expedition, which lasted from October 2019 to September 2020. The research icebreaker Polarstern , drifting with the pack ice for more than 7 months, mostly at latitudes > 85° N, served as a platform for state-of-the-art remote sensing of aerosols and clouds. The use of the recently introduced dual field-of-view (FOV) polarization lidar technique in combination with the well-established lidar-radar retrieval technique provided, for the first time, a robust instrumental basis to monitor the evolution of the liquid and the ice phase of MPCs and the interplay between the two phases. Two long-lasting Arctic MPC events observed close to the North Pole in mid-winter (December 2019) and late summer (September 2020) are discussed to provide new insight into Arctic MPC evolution processes. In the second part of the article, cloud statistics, covering all seasons of a year, are presented. The focus is on the optical and microphysical properties of the liquid phase. These results are solely derived from the dual-FOV lidar observations. The key findings of the study can be summarized as follows: persistent activation of aerosol particles to form water droplets is of great importance for the longevity of MPCs. The observations confirm that ice formation occurs predominantly via immersion freezing. The field studies suggest that the free tropospheric reservoirs of cloud condensation nuclei (CCN) and of ice-nucleating particles (INPs) were always well filled, i.e., the clouds did not exhaust their supply of activatable and activated particles. The observation of long-lasting MPC events, low ice production rates, and a sufficiently large INP reservoir leads to the recommendation to use a time-dependent immersion freezing parameterization in MPC modeling efforts.

Jimenez, Cristofer [Leibniz Inst. for Tropospheric↗

DOE OSTI2025

QCD–Gravity Double Copy in Regge Asymptotics: From $2\rightarrow n$ Amplitudes to Radiation in Shockwave Collisions

This paper discusses multi-particle production in QCD and in gravity at ultrarelativistic energies, their double-copy relations, and strong parallels in emergent shockwave dynamics. Dispersive techniques are applied to derive the BFKL equation for multi-gluon production in Regge asymptotics. Identical methods apply in gravity and are captured by a gravitational Lipatov equation. The building blocks in both cases are Lipatov vertices and reggeized propagators satisfying double-copy relations; in gravity, Weinberg’s soft theorem is recovered as a limit of the Lipatov framework. BFKL evolution in QCD generates wee parton states of maximal occupancy characterized by an emergent semi-hard saturation scale. Renormalization group equations in the Color Glass Condensate (CGC) EFT describe wee parton correlations and their rapidity evolution. A shockwave picture of deeply inelastic scattering and hadron–hadron collisions follows, with multi-particle production described by Cutkosky’s rules in strong time-dependent fields. Gluon radiation in the CGC EFT has a double copy in gravitational shockwave collisions, with a similar correspondence applicable between gluon and graviton shockwave propagators. Possible extensions of this semi-classical double copy are outlined for computing multi-particle production in gravitational shockwave collisions, self-force and tidal contributions, and classical and quantum noise in the focusing of geodesics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Development of an In Situ Fission Gas Release Instrument for Fuel Sample Irradiations in the High Flux Isotope Reactor

Experimental measurement of gaseous fission product release with respect to temperature and burnup is a critical aspect of understanding nuclear fuel performance, validating predictive models, and qualifying new fuels. To measure this phenomenon in real-time, Oak Ridge National Laboratory has developed an instrument for measuring in situ fission gas release from small-scale fuel samples irradiated in the High Flux Isotope Reactor (HFIR). The instrument uses a continuous flow of Heover the fuel samples to sweep gaseous fission products from a sealed capsule in the HFIR Be reflector to an instrument station adjacent to the reactor. The instrument station houses two high-purity germanium (HPGe) detectors that measure decay gamma rays from fission products passing through a room temperature dwell chamber placed over the detector crystal. The sealed capsules in the reactor are designed to modulate fuel sample temperatures between 700 and 1,100°C by changing the Ar/He gas mixture surrounding the capsules during irradiation. N-type thermocouples are incorporated into the capsule housing to record real-time fuel temperatures. The capsules are heated primarily by prompt gamma rays emitted from the HFIR core with minimal heat contributions from fission in the fuel samples to minimize temperature gradients in the specimens for separate-effects characterization of the material. This paper describes modeling of time-dependent nuclear heating and fission product formation in fuel samples, thermal characteristics of the in-core capsules, and expected gaseous fission product gamma spectra at the HPGe instrument station.

Mulligan, Padhraic L [ORNL] (ORCID:000000025826540↗

DOE OSTI2026

MARVEL Reactor Fuel Performance Report

The Microreactor Applications Research Validation and EvaLuation (MARVEL) project is producing a high temperature liquid metal-cooled nuclear test bed at Idaho National Laboratory (INL) to ultimately improve the integration of microreactors to end-user applications. This ambitious effort seeks to design, authorize, construct, test, and operate the reactor within five years. In order to construct and operate the MARVEL reactor in a timely manner, the system will utilize materials and component designs which have already been used, qualified, or licensed from previous reactors. The MARVEL reactor will be located at the INL Transient Reactor Test (TREAT) facility in the north high-bay equipment pit and will use the existing 304 stainless steel-clad U-ZrH1.6 pin-type fuel system developed by General Atomics and purchased from TRIGA International. This fuel has been previously qualified under the United States Department of Energy’s (US DOE) Reduced Enrichment for Research and Test Reactors (RERTR) Program. Even though the regulator of the MARVEL reactor is the US DOE, the standards and overall approach recommended by the Nuclear Regulatory Commission is well-defined and utilized here. Following NUREG-1537 regulatory guidance, this report documents the authorization case for the MARVEL fuel system’s application to MARVEL and establishes stable and predictable fuel performance during the most thermophysically unfavorable conditions achievable in the MARVEL reactor. To that end, this report provides a comprehensive survey of the known thermophysical properties, performance, and quantitative relationships associated with the MARVEL reactor fuel element and uses this information to determine its mechanical integrity and risk of reaching unacceptable conditions during the most extreme accident scenarios predicted for the reactor using the most conservative assumptions available. The information contained herein is compiled from a combination of historical reports and peer reviewed scientific publication manuscripts. Known mechanisms under which the fuel is susceptible to failure are highlighted and compared to conditions that could exist in the MARVEL reactor during an unanticipated transient or accident scenario. The two scenarios considered for analysis in this report are (1) an unprotected loss of flow accident and (2) a hypothetical unprotected loss of coolant accident during the loss of flow accident. Preliminary 2D steady-state analyses herein indicate that both fuel-cladding chemical interactions and fuel-cladding mechanical interactions are negligible throughout the fuel’s operational cycle under both normal and high temperature accident scenario conditions. Although higher fidelity 3D time-dependent modeling and simulations are planned, the following may be concluded presently. The MARVEL fuel element maintains its geometric stability and structural integrity during the most extreme accident scenarios predicted for the MARVEL reactor. The hoop stress during the unprotected loss of flow accident reaches about -1.27 MPa; this negative stress indicates that it is compressive rather than tensile. The compressive stress is a result of the NaK pressure on the outside of the fuel element, caused by the restricted thermal expansion of the NaK coolant, exceeding the internal pressures generated inside of the fuel element. The hoop stress generated in the cladding during the unprotected loss of coolant accident reaches a maximum of about approximately 10 MPa, which is nearly an order of magnitude less than the predicted yield strength of the cladding under high-temperature accident scenario conditions. Calculations were compared with results from high performance computational simulations using BISON and are in very close agreement. A conservative MARVEL fuel meat peak temperature limit of 900 °C is recommended presently, which is about 180 °C higher than the peak fuel temperature predicted to occur during the most extreme accident. Based on the known properties and behavior of the MARVEL fuel element, the fuel successfully meets its design and safety requirements under normal and most extreme accident conditions with a large safety margin.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

PSA 2025 DPRA for Cyber Optimization

Cyberattacks can have many different attack paths, durations, and goals. There are also many different mitigation options involving hardware, software, and/or humans. Evaluating defense options should include quantitative evaluation of overall effectiveness to make cost and risk-informed decisions. Typical cyberattack modeling methods only provide a qualitative evaluation and have difficulty with time dependent scenarios. The main areas of cybersecurity are confidentiality, integrity, and availability. For companies with cyber-physical systems such as advanced nuclear reactors, cyber-related integrity is a requirement set by the U.S. Nuclear Regulatory Commission. But companies are also concerned about availability or reliability as a business case. As cyber threats are evolving to a business-for-hire structure, more attacks focus on disrupting business success and reliability, causing financial and economic stability risk. Companies want reliability analysis while optimizing cost, which requires more than safety modeling methods. Dynamic-state-based and Markov-based modeling provides a method for better cyber scenario modeling with timing and conditional features not found in other numerical evaluation methods. EMRALD (Event Modeling Risk Assessment using Lined Diagrams) is a dynamic risk analysis modeling and simulation tool and has features that reduce modeling issues such as state-base explosion found in Markov-based tools. It has been used to model different time-dependent events including plant behavior and operator procedures. As a general modeling tool, EMRALD can also be used to model cyberattack scenarios with varying mitigation options and quantify effectiveness, producing numerical data for risk-informed decisions. This paper uses EMRALD to demonstrate that dynamic risk analysis can be used for cyber threat modeling to provide insights for design decision-making and optimize defense strategies.

97 - MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Constraints on Dynamical Dark Energy from Multiple Probes in the Full Dark Energy Survey

We present results on dark energy evolution, assuming a time-dependent equation of state $w(a)=w_0+w_a(1-a)$, from growth and geometric probes using the full six-year Dark Energy Survey dataset: type Ia supernovae, baryon acoustic oscillations, and weak gravitational lensing and galaxy clustering (3$\times$2pt). The combination yields $w_0=-0.84^{+0.10}_{-0.10}$ and $w_a=-0.44^{+0.60}_{-0.55}$, the tightest constraints ever obtained from a single survey, with $2.2σ$ deviation from a cosmological constant. Adding the DESI DR2 BAO data yields $w_0=-0.84^{+0.06}_{-0.07}$ and $w_a=-0.53^{+0.33}_{-0.28}$, representing the most stringent low-redshift-only test of dynamical dark energy to date, with a $2.3σ$ deviation. In this combination, adding 3$\times$2pt doubles the constraining power. Finally, when combined with primary CMB information, we obtain $w_0=-0.82^{+0.05}_{-0.05}$, $w_a=-0.63^{+0.21}_{-0.18}$, with a $3.0σ$ deviation. We find that including 3$\times$2pt in the previously studied SN + DESI BAO + CMB combination leaves the significance essentially unchanged ($3.2 σ$ to $3.0σ$) while improving the figure of merit by $\sim$10%. We systematically investigate the impact of leaving out each one of the probes and find that the significance of the deviation from a cosmological constant ranges from 2.3 to 3.2$σ$, with best-fit parameters consistently in the region $w_0 >-1$ and $w_a <0$. Excluding SN from the all data combination yields a $2.6σ$ departure from $Λ$CDM, providing a cross-check independent of supernova photometric calibration. These results support the weak preference for evolving dark energy reported by several recent cosmological analyses. By combining growth and geometric probes from a single survey, this work realizes the multi-probe dark energy program envisioned at the inception of DES.

Abbott, T. M.C. [Cerro-Tololo InterAmerican Obs.] ↗

DOE OSTI2026

Constraining the Galactic bar using the M92 stellar stream

Stellar streams are excellent probes of the gravitational potential in which they evolve. In the Milky Way (MW), globular cluster (GC) streams are routinely used to infer properties about time-dependent perturbations of the underlying potential. This implies that streams with Galactocentric radii small enough to be perturbed by the MW bar should offer constraints on it, such as its pattern speed, which currently has a wide range of values reported in the literature and is important when studying stellar kinematics. The GC M92 has a small pericentre and should be affected by the bar. It has a diffuse stellar stream, but confirming stream members has previously been hindered by a lack of spectroscopic data. In this paper, we use Dark Energy Spectroscopic Instrument (DESI) observations together with photometric and astrometric data to obtain spectroscopic members of the M92 stream for the first time. We identify a clear spatial distribution and gradients in distance moduli, proper motions, and radial velocities that confirm the stream's existence. We compare the observed stream to mock streams generated in different barred potentials and estimate the MW bar's pattern speed $Î©= 29.1^{+0.7}_{-0.4}$ km s$^{-1}$ kpc$^{-1}$ and $\dot Î©= 0.7^{+3.5}_{-2.3}$ km s$^{-1}$ kpc$^{-1}$ Gyr$^{-1}$. This is the first time a stellar stream is used to probabilistically infer these bar properties, and it opens up an exciting realm of inner Galactic potential characterisation using stellar streams.

Byström, Amanda [Edinburgh U., Inst. Astron.] (ORC↗

DOE OSTI2026

Integrating Immersive Visualization in Molten-Salt Reactor Waste Management for Experimental Design and Planning

Molten-salt reactors (MSRs) represent a promising solution for next-generation nuclear energy, offering advantages in safety, fuel efficiency, and waste minimization. However, their liquid-fueled design presents unique challenges for spent fuel management, making post-shutdown waste characterization essential for developing effective strategies. Despite this need, there is a notable absence of visualization platforms specifically tailored to the unique characteristics and analytical requirements of MSR waste management. Existing tools in the nuclear industry are primarily designed for reactor operations or generic data exploration and lack both integration with MSR-specific multiphysics frameworks and the ability to simultaneously visualize time-dependent thermal fields, chemical composition evolution, and radiation distribution patterns. To address these limitations, this paper presents an immersive virtual reality (VR) visualization platform that processes and displays high-fidelity multiphysics simulation output from the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework in real-time, using Unity. The platform visualizes MSR waste characteristics such as nuclide decay, salt cooling, and corrosion by using Exodus II output data and running on a VR headset. It includes a user-friendly interface with features such as visibility toggling, cross-sectional slicing, and time-series animation for exploring simulation data. These capabilities support experimental design, stakeholder engagement, and public communication by making complex reactor behavior more accessible and understandable. By enhancing spatial reasoning and reducing cognitive load, this immersive environment fosters more effective communication and decision-making in MSR waste management.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

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

Zinc(II) Monomeric, Dimeric, and Trimeric Photosensitizers with Microsecond‐Lived Intra‐ligand Charge Transfer Excited States Investigated through Time‐Resolved Optical and X‐Ray Spectroscopy

Zn II photosensitizers relative to Cu I complexes have received less attention due to their energetically higher metal‐to‐ligand charge transfer states. Three Zn II complexes, namely a Monomer, a bimetallic helicate, and a trimetallic helicate, bearing phenanthroline ligands are hereby studied through time‐resolved X‐ray absorption (tr‐XAS) and femto‐microsecond optical transient absorption spectroscopy (OTA). The formation of intraligand singlet charge transfer ( I ILCT) excited states is achieved within femtoseconds, followed by intersystem crossing (ISC) in nanoseconds to generate microsecond‐lived triplet ( 3 ILCT) states. Femtosecond OTA shows that the 1 ILCT states in the Monomer, Dimer, and Trimer occur within 235 fs, 683 fs, and 730 fs, respectively, while nano‐microsecond OTA and tr‐XAS show their 3 ILCT states to decay within 1.00 µs, 1.48 µs, and 1.51 µs. The ISC from the 1 ILCT to the 3 ILCT state for the Trimer is 42.8 ns compared to the Monomer and Dimer with ISC rates of less than 13 ns. These differences arise due to the stabilization by π‐π and CH‐π noncovalent interactions of the phenanthroline ligands. The dihedral and torsional angles indicate stronger ligand strains in the excited states of the Dimer and Trimer versus the Monomer. DFT calculations for the electrochemical oxidation potentials further highlight their capability in inducing photoredox processes.

earth-abundant Zn-based photosensitizers↗

DOE OSTI2025

srlife : A software tool for estimating the life of high temperature concentrating solar receivers. Part II – Ceramic receivers

As Concentrating Solar Power (CSP) technologies aim for higher operating temperatures to enhance efficiency and meet industrial process heat demands, high-temperature metallic materials, including nickel-based superalloys, face challenges in maintaining structural integrity. Advanced ceramics offer a promising alternative due to their superior high-temperature strength. However, accurately assessing the performance of ceramic components requires a fundamentally different approach from that used for metallic components. This Part II of a two-part paper describes the integration of ceramic statistical failure models within srlife – an open-source tool for predicting the life of high-temperature CSP receivers. These models account for the inherent variability in ceramic strength, as well as the effects of subcritical crack growth (SCG) under high temperature cyclic loads. Here, the paper includes an example problem that demonstrate the process of evaluating ceramic receivers using srlife. Part I details the life estimation process for metallic receivers (i.e. creep-fatigue life) along with input and output data structure, thermohydraulic analysis, and structural analysis. The complete tool is available as open-source software at https://github.com/srlife-project/srlife and can be installed via the PyPi package manager (https://pypi.org). By supporting both ceramic and metallic receiver analyses, srlife facilitates fair comparisons between competing metallic and ceramic designs, enabling accurate evaluations of plant efficiency and the economic benefits of ceramic solar receivers and other components.

High temperature ceramic receivers↗