Search NASA⌕ Search

SEARCH · Search NASA

Results for “Complete linearization method”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

Wildfire management decisions outweigh mechanical treatment as the keystone to forest landscape adaptation

Modern land management faces unprecedented uncertainty regarding future climates, novel disturbance regimes, and unanticipated ecological feedbacks. Mitigating this uncertainty requires a cohesive landscape management strategy that utilizes multiple methods to optimize benefits while hedging risks amidst uncertain futures. We used a process-based landscape simulation model (LANDIS-II) to forecast forest management, growth, climate effects, and future wildfire dynamics, and we distilled results using a decision support tool allowing us to examine tradeoffs between alternative management strategies. We developed plausible future management scenarios based on factorial combinations of restoration-oriented thinning prescriptions, prescribed fire, and wildland fire use. Results were assessed continuously for a 100-year simulation period, which provided a unique assessment of tradeoffs and benefits among seven primary topics representing social, ecological, and economic aspects of resilience. Projected climatic changes had a substantial impact on modeled wildfire activity. In the Wildfire Only scenario (no treatments, but including active wildfire and climate change), we observed an upwards inflection point in area burned around mid-century (2060) that had detrimental impacts on total landscape carbon storage. While simulated mechanical treatments (~ 3% area per year) reduced the incidence of high-severity fire, it did not eliminate this inflection completely. Scenarios involving wildland fire use resulted in greater reductions in high-severity fire and a more linear trend in cumulative area burned. Mechanical treatments were beneficial for subtopics under the economic topic given their positive financial return on investment, while wildland fire use scenarios were better for ecological subtopics, primarily due to a greater reduction in high-severity fire. Benefits among the social subtopics were mixed, reflecting the inevitability of tradeoffs in landscapes that we rely on for diverse and countervailing ecosystem services. This study provides evidence that optimal future scenarios will involve a mix of active and passive management strategies, allowing different management tactics to coexist within and among ownerships classes. Our results also emphasize the importance of wildfire management decisions as central to building more robust and resilient future landscapes.

54 ENVIRONMENTAL SCIENCES↗

Elucidating the Mechanical Energy for Cyclization of a DNA Origami Tile

DNA origami has emerged as a versatile method to synthesize nanostructures with high precision. This bottom-up self-assembly approach can produce not only complex static architectures, but also dynamic reconfigurable structures with tunable properties. While DNA origami has been explored increasingly for diverse applications, such as biomedical and biophysical tools, related mechanics are also under active investigation. Here we studied the structural properties of DNA origami and investigated the energy needed to deform the DNA structures. We used a single-layer rectangular DNA origami tile as a model system and studied its cyclization process. This origami tile was designed with an inherent twist by placing crossovers every 16 base-pairs (bp), corresponding to a helical pitch of 10.67 bp/turn, which is slightly different from that of native B-form DNA (~10.5 bp/turn). We used molecular dynamics (MD) simulations based on a coarse-grained model on an open-source computational platform, oxDNA. We calculated the energies needed to overcome the initial curvature and induce mechanical deformation by applying linear spring forces. We found that the initial curvature may be overcome gradually during cyclization and a total of ~33.1 kcal/mol is required to complete the deformation. These results provide insights into the DNA origami mechanics and should be useful for diverse applications such as adaptive reconfiguration and energy absorption.

59 BASIC BIOLOGICAL SCIENCES↗

Dual-energy CT in the obese: a preliminary retrospective review to evaluate quality and feasibility of the single-source dual-detector implementation

Purpose: To determine the feasibility of performing dual-energy CT with a single-source spectral detector system in obese patients. Materials and methods: Retrospective, IRB-approved review of 28 patients weighing ≥ 270 lbs (122 kg) who underwent CT of the abdomen on a single-source spectral detector system was performed. Two blinded, independent radiologists rated relative preference between conventional CT images taken at 120 kVp (CCT120) and monoenergetic 70 keV equivalent (MonoE70) as well as iodine map image quality in the spleen, pancreas, kidneys, and liver. Signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR) were compared between conventional CT and MonoE70 images and correlated with body habitus markers of weight, height, and abdominal diameter. Results: MonoE70 images were preferred by radiologists 100% of the time (1-sample t test, p < 0.0001) over conventional CCT120 images. Noise was significantly lower; SNR and CNR were significantly higher in MonoE70 images than in CCT120 images (paired t tests, p < 0.0001). Mean iodine map rating (scale 1–5) was 4.54 ± 0.58, denoting near homogenous and complete iodine mapping through the spleen, pancreas, kidneys, and liver for the majority of patients. Body habitus markers were not significantly correlated with image preference score; noise; MonoE70 SNR; MonoE70 CNR; change in noise, SNR, or CNR from CCT120 to MonoE70, or iodine map quality; ordinal and linear regression, p = 0.2547, p = 0.6837, p = 0.1888, p = 0.5489, p = 0.9830, p = 0.8849, p = 0.8741, p = 0.1522, respectively. Conclusion: The single-source spectral detector implementation of dual-energy CT provides viable, high-quality imaging for obese patients.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Enhancement of EM Pump Performance through Modeling and Testing

Liquid sodium has seen development as a nuclear reactor coolant in the United States for over 50 years. While sodium has many desirable characteristics as a heat transfer fluid, its reactivity poses numerous engineering challenges for even the most basic plant components like pumps. Despite the challenges, mechanical pumps have seen success in sodium service albeit at the expense of simplicity and flexibility. In fact, significant failures of these pumps have been recorded which can be attributed to their complexity. Notwithstanding the complexity, mechanical pumps are still widely proposed for the main circuit pump of sodium cooled reactors. However, an alternative to mechanical pumps are Electromagnetic Induction Pumps (EMIPs) which use the principle of the Lorentz Force. EMIPs avoid mechanically moving parts in contact with sodium which simplifies the design by avoiding mechanical wear and seals. Therefore, these pumps are generally safer and more reliable than their mechanical counterparts. That being said, EMIPs suffer from relatively low efficiencies that range from 5% to 45%. Some of the inefficiencies are inherent to the EMIPs pumps; namely, inefficiencies arise in the form of resistive heating losses in the fluid itself. However, other inefficiencies arise from Edge Effects which include Finite Length and Finite Width effects. Due to the significant safety advantages EMIPs have compared to mechanical pumps, there is a strong motivation to integrate these technologies into the design of Generation IV liquid-metal-cooled nuclear power plants. However, EMIPs have low efficiencies partly attributed to Edge Effects. Much of the literature focusing on Edge Effects is dated or simply does not exist in the public domain. More recently, numerical models of EMIPs have shown promising results. However, these models are often custom solvers or use complex software suits which limit their application. This report describes the work completed to understand Edge Effects in two EMIPs designs. These designs are the Annular Linear Induction Pump (ALIP) and Permanent Magnet Induction Pump (PMIP) also known as a Moving Magnet Pump (MMP). The methodology developed during this project aim at measuring key performance parameters (pressure-flowrate and efficiency-slip performance) using high temperature liquid-sodium. Baseline experimental measurements were then used to verify the accuracy of numerical models developed in the software suit ANSYS/FLUENT. On the experimental side, instrumentation was developed to characterize the pressure-flowrate and efficiency-slip performance of both the ALIP and PMIP in high-temperature liquid-sodium. Additionally, instrumentation was developed to measure potential output pressure pulsations and the magnetic field distributions. Moreover, a specialized ALIP was constructed to study the impacts magnetic field tapering, also known as coil grading, on the Finite Length Effect. Lastly, experimental performance curves of the PMIP were collected at 300 oC and 600 oC and compared to analytical theory which accounted for Finite Width Effects. These results showed that the Finite Width Effect and non-uniform magnetic field in a PMIP impact the performance. On the numerical side, models were developed using the software suit FLUENT/ANSYS. These methods were initially compared against a simple MHD flow in a Permanent Magnet Electromagnet Flowmeter (PMEMFM). The results were significant in that they showed agreeance between the numerical model and demonstrated that PMEMFM’s could be calibrated numerically. These modeling methods were then applied to a prototypic PMIP. Results showed that this simple model provided reasonably accurate results of the PMIP p-Q performance.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Total metals, carbon, nitrogen & anion concentration data; Slate River & East River floodplains, Crested Butte, CO; May 2022-October 2022

This data package includes processed and undiluted measurements for metal, total carbon, total nitrogen, and anion concentrations from pore water (groundwater) and surface water samples from the Slate River and East River floodplains of Crested Butte, CO, focus field sites for the SLAC Floodplain Hydro-Biogeochemistry SFA. The data was generated as part of the work targeting the overarching research question for the SLAC SFA: How do ubiquitous subsurface interfaces mediate molecular-scale biogeochemical processes and groundwater quality in floodplains and watersheds? Samples were collected between May and October of 2022. These measurements were all recorded at the Arizona Laboratory for Emerging Contaminants (ALEC) at the University of Arizona located in Tucson, AZ. Groundwater samples were extracted from a network of installed rhizon (Rhizosphere Research Products, part no. 19.60.21F, 0.6 micrometer mesh size) and piezometer wells within the river floodplain. All water samples were shaded from sun exposure during extraction from the subsurface and preserved at 4C until measured at ALEC.Analysis by ICP-MS (metals):Measurements for total metals were made on the Agilent 7700x ICP-MS (for total metals) – Agilent Technologies, Santa Clara, CA. The analytical QA/QC protocol was adapted from US EPA Method 200.8 for analysis by ICP-MS. Calibration standards were prepared from multi-element stock solutions (SPEX Certiprep, Metuchen, NJ). Calibration curves include at least 7 points with correlation coefficients > 0.995. The QC protocol includes a continuing calibration blank (CCB), a continuing calibration verification (CCV) solution and at least one quality control sample (QCS) to be analyzed just after calibration and again after every 12 samples and at the completion of the run. The QCS solutions are from an independent source, such as NIST SRM 1643e - Trace elements in water, or QCS solutions from High Purity Standards (Charleston, SC). Acceptable QC responses must be between 90 and 110% of the certified value. Lastly, a suitable internal standard (usually Rh, In, Ga or Ge) is added using on-line addition into the sample line and mixing tee.Analysis by Shimadzu TOC-L (TOC/TN):The TOC-L system is a combustion technique where liquid samples are injected and combusted into CO2 for carbon detection by non-dispersive infrared (NDIR) and NO for detection by chemiluminescence. A calibration curve using five standard solutions between 0.1 and 7 ppm for carbon and 0.05 and 3.5 ppm for nitrogen is made for each type of measurement with a linearity >0.99. All samples, standards, and QC’s are prepared in 24mL scintillation vials that have been baked for 4hrs at 475 Cº and made using RO water (18.2mΩ). QC’s include a calibration blank check (CCB), continuing calibration check (CCC), and a certified reference material check (CRM). All QC’s are within ±10% error and are run before and after each batch of samples. Samples are diluted and rerun if any measurement concentrations are above the highest standard.Analysis by Ion Chromatography (Anions):The instrument used is the Thermo Scientific Dionex ICS-6000 using AS+AG22 column set for anion analysis with sodium carbonate eluent. A calibration curve using five standard solutions between 5 and 250 umol/L is made with a linearity >0.99. Standards and QC’s are prepared in 15mL polypropylene conical tubes, pipetted along with the samples into 1.5mL polypropylene vials. Dilutions are made using RO water (18.2mΩ). QC’s include a calibration blank check (CCB), continuing calibration check (CCC), and a certified reference material check (CRM). All QC’s are within ±10% error and are run before and after each batch of samples. Samples are diluted and rerun if any measurement concentrations are above the highest standard.All files are in csv format.

54 ENVIRONMENTAL SCIENCES↗

Non-linear evolution in QCD at high-energy beyond leading order

The standard formulation of the high-energy evolution in perturbative QCD, based on the Balitsky-Kovchegov equation, is known to suffer from severe instabilities associated with radiative corrections enhanced by double transverse logarithms, which occur in all orders starting with the next-to-leading one. Over the last years, several methods have been devised to resum such corrections by enforcing the time-ordering of the successive gluon emissions. We observe that the instability problem is not fully cured by these methods: various prescriptions for performing the resummation lead to very different physical results and thus lack of predictive power. We argue that this problem can be avoided by using the rapidity of the dense target (which corresponds to Bjorken x) instead of that of the dilute projectile as an evolution time. This automatically ensures the proper time-ordering and also allows for a direct physical interpretation of the results. We explicitly perform this change of variables at NLO. We observe the emergence of a new class of double logarithmic corrections, potentially leading to instabilities, which are however less severe, since disfavoured by the typical BK evolution for “dilute-dense” scattering. We propose several prescriptions for resumming these new double-logarithms to all orders and find only little scheme dependence: different prescriptions lead to results which are consistent to each other to the accuracy of interest. We restore full NLO accuracy by completing one of the resummed equations (non-local in rapidity) with the remaining NLO corrections.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

How Useful are Conventional I–Vs for Performance Calibration of Single- and Two-Junction Perovskite Solar Cells? A Statistical Analysis of Performance Data on ≈200 Cells from 30 Global Sources

As perovskite photovoltaics (PV) advance from the laboratory to commercial prototypes, their accurate and reliable performance testing is becoming increasingly important. The well-documented dynamic response of perovskite solar cells to an external applied voltage has led to the development of steady-state performance measurement methods; however, these methods have not been widely adopted by the perovskite PV community. A key reason for this is that steady-state measurement methods take tens of minutes to complete, as opposed to conventional "fast" current-voltage (I-V) measurements usually lasting a few seconds. Fast I-Vs arise from a snapshot, almost always not a steady-state condition of the device; however, given their widespread use, the question arises: how do performance parameters of perovskite PV compare when measured with fast I-V and with a steady-state method? Results compiled from approximately 200 perovskite PV cells, including single junction, and two-terminal perovskite-perovskite and perovskite-Si tandems, show that fast I-Vs can provide a useful measure of the open-circuit voltage of the devices, while the short-circuit current and the overall efficiency can be widely misestimated. Here, the implications of these findings on performance testing protocols are discussed and possible options for fast and accurate testing of perovskite PV are proposed.

14 SOLAR ENERGY↗

The International Linear Collider. A Global Project

A large, world-wide community of physicists is working to realise an exceptional physics program of energy-frontier, electron-positron collisions with the International Linear Collider (ILC). This program will begin with a central focus on high-precision and model-independent measurements of the Higgs boson couplings. This method of searching for new physics beyond the Standard Model is orthogonal to and complements the LHC physics program. The ILC at 250 GeV will also search for direct new physics in exotic Higgs decays and in pair-production of weakly interacting particles. Polarised electron and positron beams add unique opportunities to the physics reach. The ILC can be upgraded to higher energy, enabling precision studies of the top quark and measurement of the top Yukawa coupling and the Higgs self-coupling. The key accelerator technology, superconducting radio-frequency cavities, has matured. Optimised collider and detector designs, and associated physics analyses, were presented in the ILC Technical Design Report, signed by 2400 scientists. There is a strong interest in Japan to host this international effort. A detailed review of the many aspects of the project is nearing a conclusion in Japan. Now the Japanese government is preparing for a decision on the next phase of international negotiations, that could lead to a project start within a few years. The potential timeline of the ILC project includes an initial phase of about 4 years to obtain international agreements, complete engineering design and prepare construction, and form the requisite international collaboration, followed by a construction phase of 9 years.

43 PARTICLE ACCELERATORS↗

Discontinuous Galerkin discretization for quantum simulation of chemistry

All-electron electronic structure methods based on the linear combination of atomic orbitals method with Gaussian basis set discretization offer a well established, compact representation that forms much of the foundation of modern correlated quantum chemistry calculations—on both classical and quantum computers. Despite their ability to describe essential physics with relatively few basis functions, these representations can suffer from a quartic growth of the number of integrals. Recent results have shown that, for some quantum and classical algorithms, moving to representations with diagonal two-body operators can result in dramatically lower asymptotic costs, even if the number of functions required increases significantly. We introduce a way to interpolate between the two regimes in a systematic and controllable manner, such that the number of functions is minimized while maintaining a block-diagonal structure of the two-body operator and desirable properties of an original, primitive basis. Techniques are analyzed for leveraging the structure of this new representation on quantum computers. Empirical results for hydrogen chains suggest a scaling improvement from O(N 4.5 ) in molecular orbital representations to O(N 2.6 ) in our representation for quantum evolution in a fault-tolerant setting, and exhibit a constant factor crossover at 15 to 20 atoms. Moreover, we test these methods using modern density matrix renormalization group methods classically, and achieve excellent accuracy with respect to the complete basis set limit with a speedup of 1–2 orders of magnitude with respect to using the primitive or Gaussian basis sets alone. These results suggest our representation provides significant cost reductions while maintaining accuracy relative to molecular orbital or strictly diagonal approaches for modest-sized systems in both classical and quantum computation for correlated systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The PAU Survey: an improved photo- z sample in the COSMOS field

Here we present – and make publicly available – accurate and precise photometric redshifts in the ACS footprint from the COSMOS field for objects with i AB ≤ 23. The redshifts are computed using a combination of narrow-band photometry from PAUS, a survey with 40 narrow bands spaced at $100\,\mathring{\rm A}$ intervals covering the range from 4500 to $8500\,\mathring{\rm A}$, and 26 broad, intermediate, and narrow bands covering the UV, visible and near-infrared spectrum from the COSMOS2015 catalogue. We introduce a new method that models the spectral energy distributions as a linear combination of continuum and emission-line templates and computes its Bayes evidence, integrating over the linear combinations. The correlation between the UV luminosity and the $\mathrm{O\,{\small II}}$ line is measured using the 66 available bands with the zCOSMOS spectroscopic sample, and used as a prior which constrains the relative flux between continuum and emission-line templates. The flux ratios between the $\mathrm{O\,{\small II}}$ line and H α , H β and $\mathrm{O\,{\small III}}$ are similarly measured and used to generate the emission-line templates. Comparing to public spectroscopic surveys via the quantity Δ z ≡ (z photo – z spec )/(1 + z spec ), we find the photometric redshifts to be more precise than previous estimates, with σ 68 (Δ z ) ≈ (0.003, 0.009) for galaxies at magnitude i AB ~ 18 and i AB ~ 23, respectively, which is three times and 1.66 times tighter than COSMOS2015. Additionally, we find the redshifts to be very accurate on average, yielding a median of the Δ z distribution compatible with |median(Δ z )| ≤ 0.001 at all redshifts and magnitudes considered. Both the added PAUS data and new methodology contribute significantly to the improved results. The catalogue produced with the technique presented here is expected to provide a robust redshift calibration for current and future lensing surveys, and allows one to probe galaxy formation physics in an unexplored luminosity-redshift regime, thanks to its combination of depth, completeness, and excellent redshift precision and accuracy.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantifying Rates of Erosion Using Historical Aerial Photos - 19478

Erosion can play an important role in the fate and transport of waste at sites where disposal of long-lived waste is anticipated. The statistical characterization of key processes related to erosion is essential to the understanding of site stability through time. In this work multi-temporal historical aerial images are analyzed in conjunction with orthophotography and Lidar data to develop probability distributions for variables representing important erosion processes. Scans of photographic negatives were obtained for multiple photographs and several negatives were scanned with a calibrated photogrammetric scanner by Geomni (formerly Blue Skies Consulting) in Belen, New Mexico. Other negatives were provided by USDA or scanned on a high-resolution, large format media scanner in the United States National Archives Cartographic Research Room, College Park, Maryland (1939). An additional set of aerial photographs was obtained from the USGS EROS Data Center web site. For other years, photographic prints were available for aerial photography. Several prints were scanned with a large-format media scanner at 1200 dpi. Ground control points for selected digital images used in the stereoscopic analyses were identified in ArcMap. Ground control points were collected by matching features which could be identified in the aerial photography as well as the orthophotography or shaded relief imagery from the Lidar dataset. Such features included fence line intersections, road intersections, railroad bridges, and structures, depending on the age of the photography. Many old fence lines and abandoned roads remained visible in the shaded relief imagery while they were completely obscured by tree canopy in modern photography. Ground control elevations were sampled from the Lidar-based DEM or directly from point clouds. The image coordinates and matching three-dimensional projected coordinates of the ground control points were exported to IDL where parameters for a three-dimensional direct linear transformation (DLT) were obtained. Built-in singular value decomposition routines were used to solve an overdetermined set of equations for the eleven parameters of the DLT for a particular historical photograph. A solution to the DLT parameters was obtained for each historical aerial photograph in which features of interest were identified using the stereoscope. Ultimately, this process resulted in the identification of locations for features like gully heads using both the Lidar dataset and projection of the estimated location from the historical aerial photo under consideration back to the Lidar dataset. The distance between the estimate of the location from the historical image and that of the Lidar dataset was the estimated distance the feature has moved. This distance was then divided by the number of years between the Lidar dataset and the year of the historical aerial image of interest to get a rate of movement through time. This was done for several dozen points on each historical aerial image. There was inevitable uncertainty associated with the estimation of points from the historical aerial imagery. This includes estimation of changes in locations of features like gully heads through time. In this work a method to estimate the uncertainty associated with measurements of locations of features using the stereoscopic analysis is developed and applied. For features with locations in space that do not change (e.g., road intersections, etc.) the difference in measurements of locations from the Lidar dataset and that of the estimate from the historical aerial image provides an observation of the measurement error. These points used to estimate the measurement error are referred to as control points. For each of the historical aerial images, a different measurement error distribution is estimated using control points. These measurement errors are combined with the estimated distances from the stereoscopic analyses to produce distributions of gully retreat rates and hillslope erosion rates used along with output from Erosion Working Group (EWG) Landscape Evolution Models (LEMs) to provide statistical probability distributions for input into the PPA Model. Specifically, gully head retreat rate and gully widening rate were characterized using statistical probability distributions. The central limit theorem was then applied to develop rates that represent long-term averages that are used to inform a PPA Model built in GoldSim. The importance of these rates with respect to characterization of future risk can then assessed using sensitivity analyses. This provides an approach to determine if subsequent reduction of the uncertainty associated with the probability distributions of erosion rates will correspond to meaningful reduction in the estimation of future risk from the site. Ultimately this is a general method that utilizes historical aerial imagery to estimate rates of erosion with uncertainty. This information can be used to inform PPA models to project future risks from a given site. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Benchmark Solutions for Radiation Transport in Stochastic Media with Inhomogeneous Material Statistics

Accurately solving implicit Monte Carlo (IMC) thermal photon transport problems with mixed material cells is important in realistic applications. The production IMC package at LLNL treats mixed material cells arising from ALE remap and hydrodynamics using the same approximate model. The new Imp IMC thermal photon transport package currently under development has both a material interface reconstruction (MIR) algorithm and a Levermore-Pomraning (LP) stochastic medium algorithm for treating mixed material cells. Existing stochastic medium algorithms for treating mixed material cells in IMC lack a complete theoretical basis. The IMC LP algorithm implementation has been demonstrated to reproduce published deterministic LP solutions for the particular case of spatially homogeneous material statistics. Realistic simulations will include spatially inhomogeneous material statistics (material mean chord lengths). In a previous investigation, the LP-model for transport in binary stochastic media in rod geometry was generalized to accommodate spatially varying material chord lengths, i.e., the mixing statistics were allowed to be nonhomogeneous. Analytical solutions were obtained and used to produce a verifi cation suite for the Imp IMC Levermore-Pomraning implementation for different spatial variations of the chord lengths. However, the accuracy of the LP model when the mixing statistics are nonhomogeneous has not been assessed and leaves open the question of whether local accuracy is improved or further degraded when chord lengths are not uniform. This shortcoming is rectifi ed here by developing benchmark analytic solutions for transport in binary Markovian stochastic mixtures in rod geometry with nonhomogeneous mixing statistics, using spatially varying chord lengths considered in the previous investigation based on the LP model. Methods for sampling a nonhomogeneous Poisson process (NHPP) are first described and used to construct individual realizations of the binary mixtures in rod geometry. Analytic solutions are then obtained for the forward and backward directed fluxes on a given realization, now viewed as a deterministic medium with alternating layers of the two materials with known interface locations. Finally, material averaged scalar fluxes are obtained using these sampling schemes with spatially linear and quadratic chord lengths and used to assess the accuracy of the previously obtained LP-model results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Emergence of local scaling relations in adsorption energies on high-entropy alloys

Alloying has been proposed to circumvent scaling relations between the adsorption energies thus allowing for the complete optimization of multistep reactions. Herein the fidelity of scaling rules on high-entropy alloy (HEA) surfaces is assessed focusing on hydrogen-containing molecules, *AH x for A = C and N (x = 0, 1, 2, 3), A = S (x = 0, 1, 2) and A = O (x = 0, 1). Using an adsorbate- and site-specific deep learning model to rapidly compute the adsorption energies on CoMoFeNiCu HEA surfaces, the energies of *AH x and *A are shown to be linearly correlated if *A and *AH x have identical adsorption site symmetry. However, a local linear dependence emerges between the configuration-averaged adsorption energies irrespective of the site symmetry. Although these correlations represent a weaker form of the scaling relationships, they are sufficient to prohibit the optimization of multistep reactions. The underpinning of this behavior is twofold (1) the nearsightedness principle and (2) the narrow distribution of the adsorption energies around the mean-field value. While the nearsightedness is general for all electronic systems, the second criterion applies in HEAs with relatively strong reactive elements. The present findings strongly suggest that alloys may not generally enable the breaking of scaling relationships.

36 MATERIALS SCIENCE↗

Uncertainty Quantification of Calculated Temperatures for the AGR 5/6/7 Experiment

This report documents the quantification of uncertainty of the calculated temperature data for the Advanced Gas Reactor (AGR) 5/6/7 fuel irradiation experiment conducted in the Advanced Test Reactor at Idaho National Laboratory in support of the Advanced Reactor Technologies? research and development program. Recognizing uncertainties inherent in physics and thermal simulations of the AGR 5/6/7 capsules, the results of the numerical simulations are used in combination with statistical analysis methods to improve qualification of measured data. The calculated fuel temperatures for AGR tests are also used for validation of the fission product transport and fuel performance simulation models. These crucial roles of the calculated fuel temperatures in ensuring achievement of the AGR experimental program objectives require accurate determination of the model temperature uncertainties. This report covers temperature uncertainty results for each of the five AGR 5/6/7 capsules. To quantify the uncertainty of calculated temperatures determined using the ABAQUS finite element heat transfer code, this study identifies and analyzes model parameters of potential importance to the calculated temperatures of fuel compacts and thermocouples. The selection of input parameters for uncertainty quantification is based on the ranking of their influences upon temperature predictions. Thus, selected input parameters include those with high sensitivity and those with the largest uncertainty. Propagation of model parameter uncertainty and sensitivity is then used to quantify the overall uncertainty of calculated temperatures. Measurement uncertainty, analysis of modeling assumptions, and expert judgment are used as the basis to quantify the uncertainty range for selected input parameters. The input uncertainties are dynamic, accounting for the effect of unplanned events and changes in thermal properties of capsule components over extended exposure to high temperatures and fast neutron irradiation. The sensitivity analysis performed in this work went beyond the traditional local sensitivity. Using experimental design, analysis of pairwise interactions of model parameters was performed to establish sufficiency of the time dependent first order (linear) expansion terms in constructing the temperature response surface. To achieve completeness, uncertainty propagation made use of pairwise noise correlations of model parameters. Furthermore, using an interpolation scheme over the input parameter domain, the analysis obtains time dependent sensitivity over the test campaign duration. This allows computation of uncertainty for the calculated fuel temperatures and the calculated graphite temperatures at thermocouple locations during the entire irradiation period.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A multi-dimensional parametric study of variability in multi-phase flow dynamics during geologic CO 2 sequestration accelerated with machine learning

Successful geologic CO 2 storage projects depend on numerical simulations to predict reservoir performance during site selection, injection verification, and post-injection monitoring phases of the project. These numerical simulations solve non-linear sets of coupled partial differential equations, while accounting for multi-phase fluid dynamics on the basis of constitutive equations that are embedded into the solution scheme. As a consequence, individual simulations often require tens to hundreds of hours to complete on high-performance computing clusters. Moreover, laboratory experiments reveal that parametric functions for capillary pressure and relative permeability exhibit substantial variability, even within the same rock type. This combination of computational expense and wide-ranging parametric variability means that there remains substantial uncertainty in the behavior of multi-phase CO 2 -water systems, particularly in the context of feedbacks between relative permeability and capillary pressure. To bridge this knowledge gap, here we develop a novel workflow that utilizes physics-based numerical simulation to train an artificial neural network (ANN) emulator for interrogating the multivariate parameter space that governs both capillary pressure and relative permeability. With this approach, the ANN is trained to emulate both fluid pressure distribution and CO 2 saturation, which are then interrogated quantitatively to generate parametric response surface mappings with high-fidelity resolution. Results from this study initially show that capillary entry pressure is the dominant control on both CO 2 plume geometry and fluid pressure propagation when considering the combined effects of capillary pressure and relative permeability, particularly when phase interference is low and residual CO 2 saturation is high. Moreover, the ANN emulator provides tremendous computational speed-up by computing 2691 individual simulations in several minutes; whereas, the same simulation ensemble would have required ~3 years of simulation time using only physics-based simulation methods (25,000 times speed up).

58 GEOSCIENCES↗

Mass Spectrometer Transient Analysis

This software implements a complete preprocessing pipeline for transient mass spectrometry (MS) data collected during TAP (Temporal Analysis of Products) experiments. It is designed to extract chemically meaningful fluxes from overlapping ion signals by applying a calibrated defragmentation matrix and solving the resulting linear system using non-negative least squares (NNLS) regression. The core script, preprocess_mass_spec.py, performs the following operations: Gain correction: Applies amplifier gain scalars derived from inert-packed calibration pulses to normalize signal intensities across AMUs and acquisition settings. Background subtraction: Removes experiment baselines using user-defined time windows, ensuring compatibility with slow-diffusing species and preventing negative values that would interfere with NNLS. Options to subtract before and after defragmentation. Defragmentation: Constructs a fragmentation matrix A from zeroth moments of calibration pulses (equal molar gas:inert mixtures) and solves Ax=b at each time point, where b is the raw MS signal and x is the estimated species flux. The matrix is normalized to inert signals and accounts for instrument-specific fragmentation behavior. Pulse-mode handling: Supports both averaged and individual pulse modes, enabling statistical treatment of fluxes and calculation of standard deviations. Integration and output: Computes zeroth moments (integrated fluxes) and exports time-resolved and integrated data in CSV format, suitable for downstream kinetic modeling. The software is validated using both virtual TAP simulations (VTAP) and experimental data from propane dehydrogenation (PDH) on CrOx/Al2O3 catalysts. It preserves temporal resolution by applying NNLS point-by-point across the pulse duration (typically 6,000+ time slices per pulse), leveraging the linear superposition principle to reconstruct full flux profiles. The defragmented outputs are compatible with kinetic extraction methods such as the G and Y procedures, which are used to derive rate–concentration relationships from TAP data. The details of these validations are discussed in detail in the supporting manuscript and supporting information. Example data and output files are also included. The methodology is robust to experimental noise and drift, with calibration protocols that account for pulse size effects, MS aging, and inert gas normalization. The software is modular, reproducible, and tailored for high-throughput TAP-MS workflows in catalysis research.

Kristy, Stephen [Idaho National Laboratory (INL), ↗

Quasi-Lindblad pseudomode theory for open quantum systems

Here, we introduce a new framework to study the dynamics of open quantum systems with linearly coupled Gaussian baths. Our approach replaces the continuous bath with an auxiliary discrete set of pseudomodes with dissipative dynamics, but we further relax the complete positivity requirement in the Lindblad master equation and formulate a quasi-Lindblad pseudomode theory. We show that this quasi-Lindblad pseudomode formulation directly leads to a representation of the bath correlation function in terms of a complex weighted sum of complex exponentials, an expansion that is known to be rapidly convergent in practice and thus leads to a compact set of pseudomodes. The pseudomode representation is not unique and can differ by a gauge choice. When the global dynamics can be simulated exactly, the system dynamics is unique and independent of the specific pseudomode representation. However, the gauge choice may affect the stability of the global dynamics, and we provide an analysis of why and when the global dynamics can retain stability despite losing positivity. We showcase the performance of this formulation across various spectral densities in both bosonic and fermionic problems, finding significant improvements over conventional pseudomode formulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A declining major merger fraction with redshift in the local Universe from the largest-yet catalogue of major and minor mergers in SDSS

It is difficult to accurately identify galaxy mergers and it is an even larger challenge to classify them by their mass ratio or merger stage. In previous work we used a suite of simulated mergers to create a classification technique that uses linear discriminant analysis to identify major and minor mergers. Here, we apply this technique to 1.3 million galaxies from the SDSS DR16 photometric catalogue and present the probability that each galaxy is a major or minor merger, splitting the classifications by merger stages (early, late, post-coalescence). We present publicly available imaging predictor values and all of the above classifications for one of the largest-yet samples of galaxies. We measure the major and minor merger fraction (f merg ) and build a mass-complete sample of galaxies, which we bin as a function of stellar mass and redshift. For the major mergers, we find a positive slope of fmerg with stellar mass and negative slope of f merg with redshift between stellar masses of 10.5 < M * (log M ⊙ ) < 11.6 and redshifts of 0.03 < z < 0.19. We are able to reproduce an artificial positive slope of the major merger fraction with redshift when we do not bin for mass or craft a complete sample, demonstrating the importance of mass completeness and mass binning. Further, we determine that the positive trend of the major merger fraction with stellar mass is consistent with a hierarchical assembly scenario. The negative trend with redshift requires that an additional assembly mechanism, such as baryonic feedback, dominates in the local Universe.

79 ASTRONOMY AND ASTROPHYSICS↗