Results for “model order reduction” — Page 24

DOE OSTI2023

Gas-phase Radiolysis in Plutonium Dioxide Powder

A sample of PuO 2 powder undergoes self radiolysis, causing water present in the solid to generate hydrogen, oxygen, and minor species (e.g. H 2 O 2 ). This is a matter of concern for stored PuO 2 , since the gases could react exothermically and rupture the container. In experiments at LANL, most PuO 2 /salt samples generate mostly H 2 , but a few samples generate a mixture of H 2 and O 2 together. In these special cases, the mixture of H 2 and O 2 in the pore volumes of the powder will then be subject to radiolysis, causing the reverse reaction H 2 + O 2 $\Rightarrow$ H 2 O. In order to understand the steady-state quantities of these gases, it is necessary to estimate the amount of energy from the alpha irradiation that drives this reaction. This report outlines one possible approach for this process, yielding a value for the fractional amount of radiation deposited in the gas phase of 294 ppm. This corresponds to the case of a powder with porosity of 78%. We then use a comprehensive chemical kinetics model based on elementary gas phase reactions to calculate the rate of conversion back to water. We estimate a reduction in the H 2 pressure of about 12 kPa/month, which applies only when O 2 is simultaneously present in the gas phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Determining the Electronic Structure and Thermoelectric Properties of MoS 2 /MoSe 2 Type‐I Heterojunction by DFT and the Landauer Approach

Abstract The electronic structure and thermoelectric properties of MoX 2 (X = S, Se) Van der Waals heterojunctions are reported, with the intention of motivating the design of electronic devices using such materials. Calculations indicate the proposed heterojunctions are thermodynamically stable and present a band gap reduction from 1.8 eV to 0.8 eV. The latter effect is highly related to interactions between metallic d‐character orbitals and chalcogen p‐character orbitals. The theoretical approach allows to predict a transition from semiconducting to semi‐metallic behavior. The band alignment indicates a type‐I heterojunction and band offsets of 0.2 eV. Transport properties show clear n‐type nature and a high Seebeck coefficient at 300 K, along with conductivity values (σ/τ) in the order of 10 20 . Lastly, using the Landauer approach and ballistic transport, the proposed heterojunctions can be modeled as a channel material for a typical one‐gate transistor configuration predicting subthreshold values of ≈60 mV dec −1 and field–effect mobilities of ≈160 cm −2 V −1 s −1 .

López‐Galán, Oscar A.↗

DOE OSTI2022

Deep learning to estimate permeability using geophysical data

Time-lapse electrical resistivity tomography (ERT) is a popular geophysical method to estimate three-dimensional (3D) permeability fields from electrical potential difference measurements. Traditional inversion and data assimilation methods are used to ingest this ERT data into hydrogeophysical models to estimate permeability. Due to ill-posedness and the curse of dimensionality, existing inversion strategies provide poor estimates and low resolution of the 3D permeability field. Recent advances in deep learning provide us with powerful algorithms to overcome this challenge. This paper presents a deep learning (DL) framework to estimate the 3D subsurface permeability from time-lapse ERT data. To test the feasibility of the proposed framework, we train DL-enabled inverse models on simulation data. Each measurement in both synthetic and field data is standardized by removing the mean and scaling the time-series to unit variance. This pre-processing step is necessary to bring simulation data closer to field observations. Subsurface process models based on hydrogeophysics are used to generate this synthetic data. Training performed on limited simulation data resulted in the DL model over-fitting. An advanced data augmentation based on mixup is implemented to generate additional training samples to overcome this issue. This mixup technique creates weakly labeled (low-fidelity) samples from strongly labeled (high-fidelity) data. The weakly labeled training data is then used to develop DL-enabled inverse models and reduce over-fitting. As both time-lapse ERT (1133048 features/realization) and 3D permeability (585453 features/realization) data samples are from a high-dimensional space, principal component analysis (PCA) is employed to reduce dimensionality. Encoded ERT and encoded permeability are generated using the trained PCA estimators. A deep neural network is then trained to map the encoded ERT to encoded permeability. This mixup training and unsupervised learning allowed us to build a fast and reasonably accurate DL-based inverse model under limited simulation data. Results show that proposed weak supervised learning can capture salient spatial features in the 3D permeability field. Quantitatively, the average mean squared error (in terms of the natural log) on the strongly labeled training, validation, and test datasets is less than 0.5. The R 2 -score (global metric) is greater than 0.75, and the percent error in each cell (local metric) is less than 10%. Finally, an added benefit in terms of computational cost is that the proposed DL-based inverse model is at least O(10 4 ) times faster than running a forward model once it is trained. Data generation, DL model training, and hyperparameter tuning to identify optimal neural network architectures utilized high-performance computing resources while the DL inference is performed on a standard laptop. Approximately, O(10 5 ) processor hours are used for generating data and DL tuning and training. We acknowledge that the data generation and DL model development are expensive. But once a DL model is trained, it can be re-used for inversion rapidly for the given system, with set physics and domain. Note that traditional inversion may require multiple forward model simulations (e.g., in the order of 10 to 1000), which are very expensive. This computational savings ≈ O(10 5 ) – O(10 7 )) makes the proposed DL-based inverse model attractive for subsurface imaging and real-time ERT monitoring applications due to fast and yet reasonably accurate estimations of permeability field.

58 GEOSCIENCES↗

DOE OSTI2020

Effects of Low-Carbon Energy Adoption on Airborne Particulate Matter Concentrations With Feedbacks to Future Climate Over California

California plans to reduce emissions of long-lived greenhouse gases (GHGs) through adoption of new energy systems that will also lower concentrations of short-lived absorbing soot contained in airborne particulate matter (PM). Here we examine the direct and indirect effects of reduced PM concentrations under a low-carbon energy (GHG-Step) scenario on radiative forcing in California. Simulations were carried out using the source-oriented WRF/Chem (SOWC) model with 12 km spatial resolution for the year 2054. Furthermore, the avoided aerosol emissions due to technology advances in the GHG-Step scenario reduce ground level PM concentrations by ~8.85% over land compared to the Business-as-Usual (BAU) scenario, but changes to meteorological parameters are more modest. Top-of-atmospheric forcing predicted by the SOWC model increased by 0.15 W m -2 , surface temperature warmed by 0.001 K, and planetary boundary layer height (PBLH) increased by 2.20 cm in the GHG-Step scenario compared to the BAU scenario. PM climate feedbacks are small because the significant changes in ground level PM concentrations associated with the GHG-Step scenario are limited to the first few hundred meters of the atmosphere, with little change for the majority of the vertical column above that level. As an order-of-magnitude comparison, the long-term effects of global reductions in GHG emissions (RCP8.5–RCP4.5) lowered average surface temperature over the California study domain by approximately 0.76 K. The effects of long-lived climate pollutants such as CO 2 are much stronger than the effects of short-lived climate pollutants such as PM soot over California in the year 2054.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

FAIR Data and Interpretable AI Framework for Architectured Metamaterials (Final Report)

This research program established a transformative framework for the discovery and design of mechanical metamaterials, which are architected structures engineered to control physical phenomena like sound and vibration in ways natural materials cannot. To overcome the traditional reliance on trial-and-error, the project developed an interpretable Artificial Intelligence (AI) framework that moves beyond "black box" models to reveal the specific geometric patterns—such as "unit-cell templates"—that govern a material’s performance. A major breakthrough was the development of a hierarchical design method, which allows a single material to block vibrations across multiple frequency ranges simultaneously by layering patterns at different scales without them interfering with one another. This was further expanded to include irregular, graph-based designs that use spanning tree algorithms to ensure structural connectivity while allowing for customized, direction-dependent properties like stiffness and acoustic impedance. Beyond design, the project addressed the practicalities of real-world production by developing uncertainty quantification techniques that account for manufacturing defects and material variability, reducing the need for expensive physical testing by orders of magnitude. To speed up the discovery process, the team implemented Gaussian Process Regression and other surrogate models that provide accurate performance predictions at a fraction of the traditional computational cost. The AI-generated designs were successfully validated through fabrication of physical samples and wave propagation experiments, confirming their ability to accurately guide or reflect waves as predicted. By contributing these tools and high-quality FAIR benchmark datasets to the wider scientific community, this work provides a scalable foundation for advancing technologies in aerospace vibration control, medical imaging, and noise reduction.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Continuum Modeling of Porous Electrodes for Electrochemical Synthesis

Electrochemical synthesis possesses substantial promise to utilize renewable energy sources to power the conversion of abundant feedstocks to value-added commodity chemicals and fuels. Of the potential system architectures for these processes, only systems employing 3-D structured porous electrodes have the capacity to achieve the high rates of conversion necessary for industrial scale. However, the phenomena and environments in these systems are not well understood and are challenging to probe experimentally. Fortunately, continuum modeling is well-suited to rationalize the observed behavior in electrochemical synthesis, as well as to ultimately provide recommendations for guiding the design of next-generation devices and components. In this review, we begin by presenting an historical review of modeling of porous electrode systems, with the aim of showing how past knowledge of macroscale modeling can contribute to the rising challenge of electrochemical synthesis. We then present a detailed overview of the governing physics and assumptions required to simulate porous electrode systems for electrochemical synthesis. Here, leveraging the developed understanding of porous-electrode theory, we survey and discuss the present literature reports on simulating multiscale phenomena in porous electrodes in order to demonstrate their relevance to understanding and improving the performance of devices for electrochemical synthesis. Lastly, we provide our perspectives regarding future directions in the development of models that can most accurately describe and predict the performance of such devices and discuss the best potential applications of future models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Proton internal pressure from deeply virtual Compton scattering on collider kinematics

Abstract The unique experimental connection to the QCD energy–momentum tensor offered by generalised parton distributions has been strongly highlighted in the past few years with attempts to extract the pressure and shear forces distributions within the nucleon. If, in principle, this can be performed in a model independent way from experimental data, in practice, the current limited precision and kinematic coverage make such an extraction very challenging. Moreover, the limitation to a leading-order description in the strong coupling of the data has provided only an indirect and weakly sensitive access to gluon degrees of freedom, solely through their mixing to quarks via evolution. In this paper we address this issue by providing a next-to-leading order formalism allowing a reanalysis of global fits with genuine gluonic degrees of freedom. In addition, we provide an estimate of the reduction in uncertainty that could stem from the extended kinematic range relevant for the future Electron Ion Collider. Finally, we stress the connection between the analysis of the dispersion relation in terms of generalised parton distributions and the deconvolution problem.

Dutrieux, H. (ORCID:0000000183344885)↗

DOE OSTI2025

Build-To-Replace Strategy to Reduce O&M Costs of Advanced Reactors

This paper targets the goal of reducing fixed operations and maintenance (O&M) costs to $2/MWh for advanced reactor (AR) designs, an order of magnitude reduction from current nuclear fleets. Such goal will be severely constrained if the underlying assumptions of O&M approaches and practices are not questioned and reexamined, especially today when changes can be implemented effectively and efficiently. Achieving a 90% reduction in AR O&M costs requires a completely new way of thinking – a paradigm shift – not through incremental, technology-focused approaches alone. Here, we address this challenge by evaluating the impact of moving to shorter design life for major structures, systems and components (SSCs) and shorter, more predictable refurbishment cycles as modeled by the commercial airline industry. This paper provides a brief overview of this different mindset to O&M applied to ARs: the build-to-replace approach. We provide a direct example of analysis of a build-to-replace scenario by identifying and evaluating scenarios for reduced system and component lifetimes and associated replacement and refurbishment schedules to evaluate impacts on O&M costs and other lifecycle elements such as SSC reliability.

97 - MATHEMATICS AND COMPUTING↗

DOE OSTI2021

A Transient-Response methodology based on experiments and modeling for Cu-Redox Half-Cycle kinetic analysis on a Cu-SSZ-13 SCR catalyst

In this report the Standard SCR reaction catalyzed by Cu-SSZ-13 is a redox process consisting of a reduction half cycle (RHC) and an oxidation half cycle (OHC) that cycle the active Cu sites between the Cu(II) and Cu(I) states. In the current work, a transient-response methodology consisting of experimental transient response Cu-redox (TRCR) measurements and kinetic modeling was developed for detailed study of individual SCR Cu-redox half cycles. The TRCR protocol allows quantification of the reducible Cu density, Cu(II)-Cu(I) partitioning, relative native RHC and OHC rates, and limiting half cycle during SCR. The half-cycle kinetics are studied over a wide (170–450 °C) temperature, and in differential segments along the catalyst length using spatially resolved capillary inlet mass-spectrometry (SpaciMS). The protocol alone provides two independent measures indicating that OHC increases faster than RHC with temperature, and that at all but the lowest temperature SCR is RHC limited. Introducing the transient kinetic models allows half-cycle reaction pathways, orders and activation energies to be studied and determined; e.g., RHC involves surface and gas-phase NH 3 routes, and OHC with O 2 involves dimer formation. Application of the full methodology provides further confirmation and quantification of the protocol conclusions, and specifically the first measurements of RHC and OHC activation energies; specifically for parallel RHC surface and gas-phase routes, and experimentally obtained OHC activation energy vs. first-principle calculations in the literature. Beyond insights related to the degreened catalyst studied here, the methodology provides a broadly available approach for quantifying how aging processes selectively impact RHC and OHC.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Polyakov model in ’t Hooft flux background: a quantum mechanical reduction with memory

We construct a compactification of Polyakov model on T 2 × $\mathbb{R}$ down to quantum mechanics which remembers non-perturbative aspects of field theory even at an arbitrarily small area. Standard compactification on small T 2 × $\mathbb{R}$ possesses a unique perturbative vacuum (zero magnetic flux state), separated parametrically from higher flux states, and the instanton effects do not survive in the Born-Oppenheimer approximation. By turning on a background magnetic GNO flux in co-weight lattice corresponding to a non-zero ’t Hooft flux, we show that N-degenerate vacua appear at small torus, and there are N - 1 types of flux changing instantons between them. We construct QM instantons starting with QFT instantons using the method of replicas. For example, SU(2) gauge theory with flux reduces to the double-well potential where each well is a fractional flux state. Despite the absence of a mixed anomaly, the vacuum structure of QFT and the one of QM are continuously connected. We also compare the quantum mechanical reduction of the Polyakov model with the deformed Yang-Mills, by coupling both theories to TQFTs. In particular, we compare the mass spectrum for dual photons and energy spectrum in the QM limit. We give a detailed description of critical points at infinity in the semi-classical expansion, and their role in resurgence structure.

nonperturbative effects↗

DOE OSTI2023

Modeling simultaneous particle shrinkage, dissolution and breakage using the modified moving grid technique

Simultaneous shrinkage, dissolution and breakage are important particle size reduction phenomena that characterize processes like the reactive degradation of solid chemicals. The dynamics of the particle size distribution (PSD) for such processes are non-trivial to model due to the number expending processes brought about by the eventual dissolution of particles. To this end, Population Balance Model (PBM) resolved through the sectional techniques is the natural approach. Here, we introduce a modified Moving Grid technique (m-MGT) to accurately resolve the particle size reduction phenomena. Our technique mimics the perpetual particle shrinkage through a continuously left-moving size grid and incorporates a strategic grid removal routine to capture the disappearance of particles. Coupled with the Fixed Pivot (FP) discretization for breakage, our m-MGT not only preserves the moment-related properties, but also benchmarked very well against the analytical number densities and exhibited a minimum of first-order convergence in all assessed case studies.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

GHG policy impacts for Seattle’s buildings: targets, timing, and scope

Many US cities are addressing climate change by setting goals to reduce their greenhouse gas (GHG) emissions by a specified amount within a specified time period. In order to achieve these goals, reducing emissions from existing buildings is crucial. Many cities are passing legislation to target existing buildings through benchmarking, auditing, retuning, or energy or emissions performance standards. As cities design legislation, they must consider the timing of the policies, how to prioritize building types and sizes, and how these design decisions will impact the city’s emissions. This paper addresses these questions for one particular US city: Seattle, Washington. A model of Seattle’s building stock was created with benchmarking and tax assessor data. It was then used to predict GHG emissions reductions due to different policy implementations for existing commercial and multifamily buildings. Key findings are: (1) the proposed emissions policy is expected to reduce cumulative emissions from buildings by 19% between 2020 and 2050; (2) delaying the implementation of the policy by five years could limit savings to 12%; and (3) including smaller buildings in the policy could increase savings to 34%. The lessons learned and how this can be used by other cities are discussed.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2021

Analysis of Electro-static Assisted Air Dehumidification

The removal of water vapor from the air to reduce relative humidity is a well-known indoor environmental comfort requirement. Common dehumidification approaches require a substantial amount of energy and usually involve the cooling of atmospheric humid air below its dew point or the use of absorbent/adsorbent materials to extract water vapor out of the air. More recently, researchers investigated the effect of electrostatic forces for enhancing water vapor condensation. However, the studies are limited, and there is a lack of correlations that can predict the dehumidification rate. Also, the findings from the literature focused mainly on small flow rates on the order of one to two cfm. Electrically-enhanced condensation consists of the use of highly charged particles, preferably highly charged water droplets, that attract polar water vapor molecules to their surfaces and promote condensation, a phenomenon known as dielectrophoresis. An effect of the electric charge is the reduction of the vapor pressure on the droplets' surface with respect to the saturated pressure predicted by the Kelvin equation. Consequently, the equilibrium between evaporation and condensation is shifted towards condensation. Following the application of the modified Kelvin-Thomson theory, we developed a preliminary physics-based model to predict an effective size range of the charged droplets for optimal dehumidification. The range resulted in about 2 to 4 μm in diameter, under few simplifying assumptions. The effect of the size and the charge of the electrosprayed droplets on the overall dehumidification rate was briefly discussed. The use of electrosprays to produce small but highly charged droplets was broadly discussed in this paper. The cone-jet mode was identified as the most suitable electrospray operational mode, and it generated droplets of small size and high electrical charge. The cone-jet stability was also analyzed in detail. The preliminary data of the present work and the model results indicated that several electrospray heads were required to achieve a 5% dehumidification rate for airflow rates of about 5 cfm.

Morcelli, Stefano↗

DOE OSTI2021

Analysis of New Data of Electro-static Assisted Air Dehumidification Processes

The removal of water vapor from the air to reduce relative humidity is a well-known indoor environmental comfort requirement. Common dehumidification approaches require a substantial amount of energy and usually involve the cooling of atmospheric humid air below its dew point or the use of absorbent/adsorbent materials to extract water vapor out of the air. More recently, researchers investigated the effect of electrostatic forces for enhancing water vapor condensation. However, the studies are limited, and there is a lack of correlations that can predict the dehumidification rate. Also, the findings from the literature focused mainly on small flow rates on the order of one to two cfm. Electrically-enhanced condensation consists of the use of highly charged particles, preferably highly charged water droplets, that attract polar water vapor molecules to their surfaces and promote condensation, a phenomenon known as dielectrophoresis. An effect of the electric charge is the reduction of the vapor pressure on the droplets' surface with respect to the saturated pressure predicted by the Kelvin equation. Consequently, the equilibrium between evaporation and condensation is shifted towards condensation. Following the application of the modified Kelvin-Thomson theory, we developed a preliminary physics-based model to predict an effective size range of the charged droplets for optimal dehumidification. The range resulted in about 2 to 4 m in diameter, under few simplifying assumptions. The effect of the size and the charge of the electrosprayed droplets on the overall dehumidification rate was briefly discussed. The use of electrosprays to produce small but highly charged droplets was broadly discussed in this paper. The cone-jet mode was identified as the most suitable electrospray operational mode, and it generated droplets of small size and high electrical charge. The cone-jet stability was also analyzed in detail. The preliminary data of the present work and the model results indicated that several electrospray heads were required to achieve a 5% dehumidification rate for airflow rates of about 5 cfm.

Morcelli, Stefano↗

DOE OSTI2026

Enhancing Sterile Neutrino Oscillation Sensitivities using SBND-PRISM at the Short-Baseline Neutrino Programme

The Short-Baseline Neutrino (SBN) Programme at Fermilab is comprised of two detectors, SBND and ICARUS, placed at $110$\,m and $600$\,m along the Booster Neutrino Beam at Fermilab. The key physics aim of the programme is to definitively test the sterile neutrino hypothesis, a proposed fourth flavour of neutrino that may explain certain experimental anomalies seen regarding standard model neutrino oscillations. To be able to detect the existence of sterile neutrinos, the uncertainties of the programme must be well constrained. To enable this, a robust analysis must be constructed that can consistently identify the correct values of systematic parameters and generate accurate predictions of what the neutrino energy spectrum at ICARUS should look like. SBND's design allows for the introduction of a technique called PRISM. In PRISM, the detector is divided into regions of different off axis angles from the beam, forming different samples where systematics impact each one in a distinct way. This allows any fits performed to obtain a better understanding of the correct value of the systematic parameters. The first analysis included in this thesis focuses on the improvements seen to the sensitivity of SBN to sterile oscillation parameters when using a PRISM configuration instead of treating SBND as a single whole. Focusing on the $5\sigma$ exclusion contour from the $\numu$ disappearance channel using three PRISM samples, an improvement on the order of $30$\,\% is seen, extending the parameter space for which the null hypothesis can be excluded. This thesis also presents a series of mock data studies comparing the abilities of SBND and PRISM analyses, concluding that in the case of simple changes between Monte Carlo (MC) and data, PRISM produces predictions of the ICARUS event rate spectrum that are more accurate and have smaller uncertainties. When moving to more realistic mock data, using different models to create the mock data than were used for the MC, the postfit predictions at ICARUS had a smaller difference between the postfit and mock data reconstructed energy spectra across all mock data samples tested. Finally, a covariance matrix defined by the maximum discrepancy between the postfit and mock data spectra at ICARUS from each of the SBND and PRISM fits across all the samples was constructed. The resultant $1\sigma$ fractional error induced by this bias systematic on the postfit spectrum has a smaller magnitude for PRISM than SBND, with the improvements ranging from $2.21$\,\% to $6.26$\,\% depending on the mock data samples. This summarises the reduction in systematic error when using PRISM instead of treating SBND as a single detector. *************************************** AUTHOR = Slater, Bethany University of Liverpool b.slater2@liverpool.ac.uk TITLE = Enhancing Sterile Neutrino Oscillation Sensitivities using SBND-PRISM at the Short-Baseline Neutrino Programme PAGES = 218 NOTE = Ph.D. University of Liverpool March 2026 ABSTRACT = The Short-Baseline Neutrino (SBN) Programme at Fermilab is comprised of two detectors, SBND and ICARUS, placed at $110$\,m and $600$\,m along the Booster Neutrino Beam at Fermilab. The key physics aim of the programme is to definitively test the sterile neutrino hypothesis, a proposed fourth flavour of neutrino that may explain certain experimental anomalies seen regarding standard model neutrino oscillations. To be able to detect the existence of sterile neutrinos, the uncertainties of the programme must be well constrained. To enable this, a robust analysis must be constructed that can consistently identify the correct values of systematic parameters and generate accurate predictions of what the neutrino energy spectrum at ICARUS should look like. SBND's design allows for the introduction of a technique called PRISM. In PRISM, the detector is divided into regions of different off axis angles from the beam, forming different samples where systematics impact each one in a distinct way. This allows any fits performed to obtain a better understanding of the correct value of the systematic parameters. The first analysis included in this thesis focuses on the improvements seen to the sensitivity of SBN to sterile oscillation parameters when using a PRISM configuration instead of treating SBND as a single whole. Focusing on the $5\sigma$ exclusion contour from the $\numu$ disappearance channel using three PRISM samples, an improvement on the order of $30$\,\% is seen, extending the parameter space for which the null hypothesis can be excluded. This thesis also presents a series of mock data studies comparing the abilities of SBND and PRISM analyses, concluding that in the case of simple changes between Monte Carlo (MC) and data, PRISM produces predictions of the ICARUS event rate spectrum that are more accurate and have smaller uncertainties. When moving to more realistic mock data, using different models to create the mock data than were used for the MC, the postfit predictions at ICARUS had a smaller difference between the postfit and mock data reconstructed energy spectra across all mock data samples tested. Finally, a covariance matrix defined by the maximum discrepancy between the postfit and mock data spectra at ICARUS from each of the SBND and PRISM fits across all the samples was constructed. The resultant $1\sigma$ fractional error induced by this bias systematic on the postfit spectrum has a smaller magnitude for PRISM than SBND, with the improvements ranging from $2.21$\,\% to $6.26$\,\% depending on the mock data samples. This summarises the reduction in systematic error when using PRISM instead of treating SBND as a single detector.

Slater, Bethany [Liverpool U.]↗

DOE OSTI2026

Enhancing Sterile Neutrino Oscillation Sensitivities using SBND-PRISM at the Short-Baseline Neutrino Programme

The Short-Baseline Neutrino (SBN) Programme at Fermilab is comprised of two detectors, SBND and ICARUS, placed at $110$\,m and $600$\,m along the Booster Neutrino Beam at Fermilab. The key physics aim of the programme is to definitively test the sterile neutrino hypothesis, a proposed fourth flavour of neutrino that may explain certain experimental anomalies seen regarding standard model neutrino oscillations. To be able to detect the existence of sterile neutrinos, the uncertainties of the programme must be well constrained. To enable this, a robust analysis must be constructed that can consistently identify the correct values of systematic parameters and generate accurate predictions of what the neutrino energy spectrum at ICARUS should look like. SBND's design allows for the introduction of a technique called PRISM. In PRISM, the detector is divided into regions of different off axis angles from the beam, forming different samples where systematics impact each one in a distinct way. This allows any fits performed to obtain a better understanding of the correct value of the systematic parameters. The first analysis included in this thesis focuses on the improvements seen to the sensitivity of SBN to sterile oscillation parameters when using a PRISM configuration instead of treating SBND as a single whole. Focusing on the $5\sigma$ exclusion contour from the $\numu$ disappearance channel using three PRISM samples, an improvement on the order of $30$\,\% is seen, extending the parameter space for which the null hypothesis can be excluded. This thesis also presents a series of mock data studies comparing the abilities of SBND and PRISM analyses, concluding that in the case of simple changes between Monte Carlo (MC) and data, PRISM produces predictions of the ICARUS event rate spectrum that are more accurate and have smaller uncertainties. When moving to more realistic mock data, using different models to create the mock data than were used for the MC, the postfit predictions at ICARUS had a smaller difference between the postfit and mock data reconstructed energy spectra across all mock data samples tested. Finally, a covariance matrix defined by the maximum discrepancy between the postfit and mock data spectra at ICARUS from each of the SBND and PRISM fits across all the samples was constructed. The resultant $1\sigma$ fractional error induced by this bias systematic on the postfit spectrum has a smaller magnitude for PRISM than SBND, with the improvements ranging from $2.21$\,\% to $6.26$\,\% depending on the mock data samples. This summarises the reduction in systematic error when using PRISM instead of treating SBND as a single detector. *************************************** AUTHOR = Slater, Bethany University of Liverpool b.slater2@liverpool.ac.uk TITLE = Enhancing Sterile Neutrino Oscillation Sensitivities using SBND-PRISM at the Short-Baseline Neutrino Programme PAGES = 218 NOTE = Ph.D. University of Liverpool March 2026 ABSTRACT = The Short-Baseline Neutrino (SBN) Programme at Fermilab is comprised of two detectors, SBND and ICARUS, placed at $110$\,m and $600$\,m along the Booster Neutrino Beam at Fermilab. The key physics aim of the programme is to definitively test the sterile neutrino hypothesis, a proposed fourth flavour of neutrino that may explain certain experimental anomalies seen regarding standard model neutrino oscillations. To be able to detect the existence of sterile neutrinos, the uncertainties of the programme must be well constrained. To enable this, a robust analysis must be constructed that can consistently identify the correct values of systematic parameters and generate accurate predictions of what the neutrino energy spectrum at ICARUS should look like. SBND's design allows for the introduction of a technique called PRISM. In PRISM, the detector is divided into regions of different off axis angles from the beam, forming different samples where systematics impact each one in a distinct way. This allows any fits performed to obtain a better understanding of the correct value of the systematic parameters. The first analysis included in this thesis focuses on the improvements seen to the sensitivity of SBN to sterile oscillation parameters when using a PRISM configuration instead of treating SBND as a single whole. Focusing on the $5\sigma$ exclusion contour from the $\numu$ disappearance channel using three PRISM samples, an improvement on the order of $30$\,\% is seen, extending the parameter space for which the null hypothesis can be excluded. This thesis also presents a series of mock data studies comparing the abilities of SBND and PRISM analyses, concluding that in the case of simple changes between Monte Carlo (MC) and data, PRISM produces predictions of the ICARUS event rate spectrum that are more accurate and have smaller uncertainties. When moving to more realistic mock data, using different models to create the mock data than were used for the MC, the postfit predictions at ICARUS had a smaller difference between the postfit and mock data reconstructed energy spectra across all mock data samples tested. Finally, a covariance matrix defined by the maximum discrepancy between the postfit and mock data spectra at ICARUS from each of the SBND and PRISM fits across all the samples was constructed. The resultant $1\sigma$ fractional error induced by this bias systematic on the postfit spectrum has a smaller magnitude for PRISM than SBND, with the improvements ranging from $2.21$\,\% to $6.26$\,\% depending on the mock data samples. This summarises the reduction in systematic error when using PRISM instead of treating SBND as a single detector.

Slater, Bethany [Liverpool U.]↗

DOE OSTI2026

Enhancing Sterile Neutrino Oscillation Sensitivities using SBND-PRISM at the Short-Baseline Neutrino Programme

The Short-Baseline Neutrino (SBN) Programme at Fermilab is comprised of two detectors, SBND and ICARUS, placed at $110$\,m and $600$\,m along the Booster Neutrino Beam at Fermilab. The key physics aim of the programme is to definitively test the sterile neutrino hypothesis, a proposed fourth flavour of neutrino that may explain certain experimental anomalies seen regarding standard model neutrino oscillations. To be able to detect the existence of sterile neutrinos, the uncertainties of the programme must be well constrained. To enable this, a robust analysis must be constructed that can consistently identify the correct values of systematic parameters and generate accurate predictions of what the neutrino energy spectrum at ICARUS should look like. SBND's design allows for the introduction of a technique called PRISM. In PRISM, the detector is divided into regions of different off axis angles from the beam, forming different samples where systematics impact each one in a distinct way. This allows any fits performed to obtain a better understanding of the correct value of the systematic parameters. The first analysis included in this thesis focuses on the improvements seen to the sensitivity of SBN to sterile oscillation parameters when using a PRISM configuration instead of treating SBND as a single whole. Focusing on the $5\sigma$ exclusion contour from the $\numu$ disappearance channel using three PRISM samples, an improvement on the order of $30$\,\% is seen, extending the parameter space for which the null hypothesis can be excluded. This thesis also presents a series of mock data studies comparing the abilities of SBND and PRISM analyses, concluding that in the case of simple changes between Monte Carlo (MC) and data, PRISM produces predictions of the ICARUS event rate spectrum that are more accurate and have smaller uncertainties. When moving to more realistic mock data, using different models to create the mock data than were used for the MC, the postfit predictions at ICARUS had a smaller difference between the postfit and mock data reconstructed energy spectra across all mock data samples tested. Finally, a covariance matrix defined by the maximum discrepancy between the postfit and mock data spectra at ICARUS from each of the SBND and PRISM fits across all the samples was constructed. The resultant $1\sigma$ fractional error induced by this bias systematic on the postfit spectrum has a smaller magnitude for PRISM than SBND, with the improvements ranging from $2.21$\,\% to $6.26$\,\% depending on the mock data samples. This summarises the reduction in systematic error when using PRISM instead of treating SBND as a single detector. *************************************** AUTHOR = Slater, Bethany University of Liverpool b.slater2@liverpool.ac.uk TITLE = Enhancing Sterile Neutrino Oscillation Sensitivities using SBND-PRISM at the Short-Baseline Neutrino Programme PAGES = 218 NOTE = Ph.D. University of Liverpool March 2026 ABSTRACT = The Short-Baseline Neutrino (SBN) Programme at Fermilab is comprised of two detectors, SBND and ICARUS, placed at $110$\,m and $600$\,m along the Booster Neutrino Beam at Fermilab. The key physics aim of the programme is to definitively test the sterile neutrino hypothesis, a proposed fourth flavour of neutrino that may explain certain experimental anomalies seen regarding standard model neutrino oscillations. To be able to detect the existence of sterile neutrinos, the uncertainties of the programme must be well constrained. To enable this, a robust analysis must be constructed that can consistently identify the correct values of systematic parameters and generate accurate predictions of what the neutrino energy spectrum at ICARUS should look like. SBND's design allows for the introduction of a technique called PRISM. In PRISM, the detector is divided into regions of different off axis angles from the beam, forming different samples where systematics impact each one in a distinct way. This allows any fits performed to obtain a better understanding of the correct value of the systematic parameters. The first analysis included in this thesis focuses on the improvements seen to the sensitivity of SBN to sterile oscillation parameters when using a PRISM configuration instead of treating SBND as a single whole. Focusing on the $5\sigma$ exclusion contour from the $\numu$ disappearance channel using three PRISM samples, an improvement on the order of $30$\,\% is seen, extending the parameter space for which the null hypothesis can be excluded. This thesis also presents a series of mock data studies comparing the abilities of SBND and PRISM analyses, concluding that in the case of simple changes between Monte Carlo (MC) and data, PRISM produces predictions of the ICARUS event rate spectrum that are more accurate and have smaller uncertainties. When moving to more realistic mock data, using different models to create the mock data than were used for the MC, the postfit predictions at ICARUS had a smaller difference between the postfit and mock data reconstructed energy spectra across all mock data samples tested. Finally, a covariance matrix defined by the maximum discrepancy between the postfit and mock data spectra at ICARUS from each of the SBND and PRISM fits across all the samples was constructed. The resultant $1\sigma$ fractional error induced by this bias systematic on the postfit spectrum has a smaller magnitude for PRISM than SBND, with the improvements ranging from $2.21$\,\% to $6.26$\,\% depending on the mock data samples. This summarises the reduction in systematic error when using PRISM instead of treating SBND as a single detector.

Slater, Bethany [Liverpool U.]↗

DOE OSTI2025

Decay heat analysis for advanced reactor spent fuel transportation and storage applications

Accurate characterization of nuclide inventories and decay heat in spent nuclear fuel is critical for ensuring its safe handling, storage, transportation, and disposal. Although extensive research has been conducted on light-water reactor fuel, advanced reactors present unique challenges due to their diverse core configurations, fuel characteristics, neutron energy spectra, and burnup levels. Building upon previous efforts that developed representative reactor core models for various advanced reactor types and fuels, this study evaluates reactor-specific decay heat characteristics. The results highlight significant variations across advanced reactor types as well as across reactor designs within the same reactor type, and they provide comparison to typical commercial light-water reactor fuel. For example, thermal-spectrum reactor fuels were observed to have an approximately 100-fold decrease in decay heat over the first decade of cooling, whereas the reduction was 10-fold for fast-spectrum reactor fuels. Mass-specific decay heat at discharge can differ by three orders of magnitude among the fast and thermal reactor systems considered. Overall, for the analyzed advanced reactor fuel, fewer than 17 nuclides account for over 99% of total decay heat at 0.5 years, and that number drops to fewer than 7 nuclides at 100 years of cooling. By quantifying reactor-specific decay heat trends and nuclide contributions, this work provides a technical basis to support the development of spent fuel management strategies for advanced reactor fuels as well as safety evaluations for storage, transportation, and long-term waste disposal.

Advanced reactors↗