Search NASA⌕ Search

SEARCH · Search NASA

Results for “Probability estimates”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Estimating the influence of field inventory sampling intensity on forest landscape model performance for determining high-severity wildfire risk

Abstract Historically, fire has been essential in Southwestern US forests. However, a century of fire-exclusion and changing climate created forests which are more susceptible to uncharacteristically severe wildfires. Forest managers use a combination of thinning and prescribed burning to reduce forest density to help mitigate the risk of high-severity fires. These treatments are laborious and expensive, therefore optimizing their impact is crucial. Landscape simulation models can be useful in identifying high risk areas and assessing treatment effects, but uncertainties in these models can limit their utility in decision making. In this study we examined underlying uncertainties in the initial vegetation layer by leveraging a previous study from the Santa Fe fireshed and using new inventory plots from 111 stands to interpolate the initial forest conditions. We found that more inventory plots resulted in a different geographic distribution and wider range of the modelled biomass. This changed the location of areas with high probability of high-severity fires, shifting the optimal location for management. The increased range of biomass variability from using a larger number of plots to interpolate the initial vegetation layer also influenced ecosystem carbon dynamics, resulting in simulated forest conditions that had higher rates of carbon uptake. We conclude that the initial forest layer significantly affects fire and carbon dynamics and is dependent on both number of plots, and sufficient representation of the range of forest types and biomass density.

Science & Technology - Other Topics↗
DOE OSTI2020

Future projections and uncertainty assessment of precipitation extremes in the Korean peninsula from the CMIP5 ensemble

Abstract Projections of changes in extreme climate are sometimes predicted by multimodel ensemble methods that combine forecasts from individual simulation models using weighted averaging. One method to assign weight to each model is the Bayesian model averaging (BMA) in which posterior probability is used. For the cases of extreme climate, the generalized extreme value distribution (GEVD) is typically used. We applied the approach of GEV‐embedded BMA to a series of 35 years of the annual maximum daily precipitation data (both historical data and data gathered from simulation experiments for future periods) over the Korean peninsula as simulated by the models in the Coupled Model Intercomparison Project Phase 5 (CMIP5). Simulation data under two Representative Concentration Pathway (RCP) scenarios, namely RCP4.5 and RCP8.5, were used. Observed data and 17 CMIP5 models for 12 gird cells in Korea have been examined to predict future changes in precipitation extremes. A simple regional frequency analysis of pooling observations from three stations in each cell was employed to reduce the estimation variance and local fluctuations. A bias correction technique using the regression‐type transfer function was applied to these simulation data. Return levels spanning over 20 and 50 years, as well as the return periods relative to the reference years (1971–2005), were estimated for two future periods, namely Period 1 (2021–2050) and Period 2 (2066–2095). From these analyses, relative increase observed in the spatially averaged 20‐year (50‐year) return level was approximately 23% (16%) and 45% (36%) in the RCP4.5 and RCP8.5 experiments, respectively, by the end of the 21st century. We concluded that extreme rainfalls will likely occur two times and four times more frequently in the RCP4.5 and RCP8.5 scenarios, respectively, as compared to in the reference years by the end of the 21st century.

Lee, Youngsaeng↗
DOE OSTI2021

Quantifying the Solar Energy Resource for Puerto Rico

After Hurricane Maria, multiple U.S. Department of Energy laboratories studied the state of the electric grid in Puerto Rico and analyzed grid resilience and grid integration of renewable energy. As part of the work done at the National Renewable Energy Laboratory, researchers created new solar resource data, conducted a technical potential and supply curve analysis, and studied the interannual variability of the solar resource. A new methodology was developed to downscale solar resource data from the National Solar Radiation Data Base (NSRDB) from a 4-km x 4-km spatial and 30-minute temporal resolution to a 2 km x 2 km and 5-minute resolution. This methodology primarily used simple physical principles to develop high-resolution cloud properties which were then used to compute solar radiation. The high-resolution datasets were validated against ground measurements and the error metrics were found to be similar to the original lower resolution dataset. Using 20 years of downscaled data from the NSRDB multi-year capacity factors for photovoltaics (PV) were developed for both single-axis tracking and fixed latitude-tilt configurations. Use of the multi-year data provides the ability to understand variability in capacity factors due to variability in weather over a long period of time. For Puerto Rico the coastal regions were found to have significant higher capacity factors than inland. Using land-use and terrain information a technical potential analysis was conducted for Puerto Rico. This analysis restricted single PV plant development to a maximum of 100 MW nameplate capacity. The nameplate capacity for each municipality were then determined. Based on our assumptions, 56 of the 78 total municipalities of Puerto Rico contain some level of solar capacity. Most of the interior municipalities did not have any capacity because of the geographic exclusions used in this study. The lowest capacity for a PV plant observed in a municipality was 10 MW. The maximum capacity within a county was 2,000 MW. Further a supply curve analysis was conducted by taking the results of the technical potential and quantifying system and transmission costs. The levelized cost of energy (LCOE) was calculated for each theoretical PV plant site, and the levelized cost of transmission was added to the LCOE to produce a total cost estimate for each site. The results of the supply curve analysis allow for a relative comparison of the cost for integrating new PV capacity into the grid. This analysis indicates that cheaper total LCOE sites tend to be larger in capacity. The total capacity in this study was found to be far beyond the maximum peak load for the island. However, this study does not consider the economic and market potential for development. The cumulative capacity presented in this study assumes that the best locations are developed first and ignores the complex decision paths for new power plant development. Therefore, this analysis can only be treated as illustrative. Finally this study investigates the impact of inter-annual variability of resource using a variety of metrices including probability of exceedance and variation in capacity factor and LCOE. This study demonstrates that the capacity factor or LCOE could vary by over 10% year to year. This clearly indicates the risks involved in using any particular year of data and clearly points to the use of multi-year data to reduce some of the risks related to variability in weather.

14 SOLAR ENERGY↗
DOE OSTI2022

An Enhanced Computational Framework for Energy Storage Participation in Transmission Planning with Electricity Market Participation: Dynamic Market Participation Restriction Estimation

In the previous study and report, we developed a computational framework to explore the technical and economic viability of the dual use of pumped storage hydropower (PSH) for both transmission assets and market participation. The rationale for the dual-use application is to improve the resource utilization rate and reduce the cost recovery burden. A PSH-based transmission solution generally has a large capacity, which corresponds to a higher investment cost compared to a traditional line solution. Market participation may allow the PSH asset to offset part of its high investment cost. A critical question is whether it is economically and technically viable to allow dual use of a PSH-based transmission project without jeopardizing transmission system reliability. And a more specific follow-up question is how to improve market participation performance while maintaining transmission reliability. In this study, we extend the computational framework by updating the method to determine transmission service obligation requirements dynamically based on the possibility of the need for transmission service in operation. Specifically, we correlate the need for transmission service with the probability of thermal overloading on a monitored transmission line, which is essentially based on the probability distribution of power flow and its probability of exceeding the line rating based on system dispatch and operational uncertainties. In addition, the need depends on the level of risk that a system operator is willing to take. The higher the risk tolerance, the lower the requirement that the system operator would apply to the PSH project as a transmission asset. The requirements are updated every day before the PSH project determines its market participation strategy. With the updated transmission service requirements, it is expected that a PSH project can serve as a transmission asset when the transmission system is at risk and really needs it, while it can participate in a market with sufficient information and certain restrictions. In this way, the transmission service quality and market participation performance can both be improved.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Mapping variations of redshift distributions with probability integral transforms

We present a method for mapping variations between probability distribution functions and apply this method within the context of measuring galaxy redshift distributions from imaging survey data. This method, which we name PITPZ for the probability integral transformations it relies on, uses a difference in curves between distribution functions in an ensemble as a transformation to apply to another distribution function, thus transferring the variation in the ensemble to the latter distribution function. This procedure is broadly applicable to the problem of uncertainty propagation. In the context of redshift distributions, for example, the uncertainty contribution due to certain effects can be studied effectively only in simulations, thus necessitating a transfer of variation measured in simulations to the redshift distributions measured from data. We illustrate the use of PITPZ by using the method to propagate photometric calibration uncertainty to redshift distributions of the Dark Energy Survey Year 3 weak lensing source galaxies. For this test case, we find that PITPZ yields a lensing amplitude uncertainty estimate due to photometric calibration error within 1 per cent of the truth, compared to as much as a 30 per cent underestimate when using traditional methods.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Kinetic theory of particle-in-cell simulation plasma and the ensemble averaging technique

Abstract We derive the kinetic theory of fluctuations in physically and numerically stable particle-in-cell (PIC) simulations of electrostatic plasmas. The starting point is the single-time correlation at the start of the simulation between the statistical fluctuations of the weighted densities of macroparticle centers in the plasma particle phase-space. The fluctuations are associated with different initial conditions, typically due to the random initial conditions (in velocity space) of the macroparticles/simulation plasma, assigned according to their initial distribution of probability. The single-time correlations at all time steps and in each spatial grid cell are then determined from the Laplace–Fourier transforms of the discretized Klimontovich-like equation for the macroparticles and Maxwell’s equations for the fields, as computed by modern PIC codes. We recover the expressions for the electrostatic field and the plasma particle density fluctuation autocorrelation spectra as well as the kinetic equations describing the average evolution of PIC-simulated plasma particles, first derived by Langdon (1970b Proc. 4th Conf. Numerical Simulation of Plasmas ) using a test macroparticle approach perturbing a discretized Vlasovian plasma and then averaging the obtained physical quantity over the initial macroparticle velocity distribution. We generalize and extend these results to the modern algorithms in PIC codes using arbitrary macroparticle weights. Analytical estimates of statistical fluctuation amplitudes are derived as a function of the plasma simulation parameters, using the central limit theorem in the limit of a large number of macroparticles per cell. The theory is then used to analyze the ensemble averaging technique of PIC simulations where statistical averages are performed over ensembles of PIC simulations, modeling the same plasma physics problem but using different statistical realizations of the initial distribution functions of the macroparticles. This method is illustrated by linear Landau damping uncovering (from noise, which is usually considered numerical) the physical fluctuations driven by a single small amplitude electrostatic wave perturbing a PIC simulation plasma in equilibrium.

fluctuations correlations↗
DOE OSTI2023

Identification of Galaxy–Galaxy Strong Lens Candidates in the DECam Local Volume Exploration Survey Using Machine Learning

We perform a search for galaxy–galaxy strong lens systems using a convolutional neural network (CNN) applied to imaging data from the first public data release of the DECam Local Volume Exploration Survey, which contains ~520 million astronomical sources covering ~4000 deg 2 of the southern sky to a 5σ point–source depth of g = 24.3, r= 23.9, i = 23.3, and z = 22.8 mag. Following the methodology of similar searches using Dark Energy Camera data, we apply color and magnitude cuts to select a catalog of ~11 million extended astronomical sources. After scoring with our CNN, the highest-scoring 50,000 images were visually inspected and assigned a score on a scale from 0 (not a lens) to 3 (very probable lens). We present a list of 581 strong lens candidates, 562 of which are previously unreported. We categorize our candidates using their human-assigned scores, resulting in 55 Grade A candidates, 149 Grade B candidates, and 377 Grade C candidates. We additionally highlight eight potential quadruply lensed quasars from this sample. Due to the location of our search footprint in the northern Galactic cap (b > 10 deg) and southern celestial hemisphere (decl. < 0 deg), our candidate list has little overlap with other existing ground-based searches. Where our search footprint does overlap with other searches, we find a significant number of high-quality candidates that were previously unidentified, indicating a degree of orthogonality in our methodology. We report properties of our candidates including apparent magnitude and Einstein radius estimated from the image separation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

DEMONSTRATION OF MULTI-PHYSICS EVALUATION OF FUEL ROD BURST PROBABILITY WITH BURNUP EXTENSION UNDER LBLOCA CONDITIONS

Extending the fuel discharge burnup level, e.g., from the current limit of rod averaged discharge burnup limit of 62 GWD/MT to a proposed new limit of 75 GWD/MT, can provide significant economic benefits to the current fleet of operating light water reactors (LWRs). It allows for longer operating cycles and improved resource utilization. The major economic gain of longer operating cycles is attributable to the increased capacity factor resulting from decreased refueling time as a fraction of total operating time, as well as fewer assemblies to be discharged for a given amount of energy produced. The main licensing challenges for higher burnup fuel are to ensure fuel rod safety under design basis accident conditions, especially under large-break loss-of-coolant accident (LBLOCA) and reactivity insertion accident (RIA). In this work, two-year cycle core design for a typical 4-loop pressurized water reactor (PWR) is performed with enrichment increased up to 6% and burnup extended to 75 GWD/MT. The fuel rod burst potential evaluations under large-break loss-of-coolant accident (LBLOCA) conditions are subsequently performed using the multi-physics best estimate plus uncertainty analysis framework LOTUS (LOCA Toolkit for the U.S. LWRs) and the preliminary results are presented.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

An Integrated Framework for Risk Assessment of High Safety Significant Safety-related Digital Instrumentation and Control Systems in Nuclear Power Plants: Methodology and Demonstration

This report documents the activities performed by Idaho National Laboratory (INL) during Fiscal Year (FY) 2022 for the U.S. Department of Energy (DOE) Light Water Reactor Sustainability (LWRS) Program, Risk Informed Systems Analysis (RISA) Pathway, digital instrumentation and control (DI&C) risk assessment project. In FY 2019, the RISA Pathway initiated a project to develop a risk assessment strategy for delivering a technical basis to support effective and secure DI&C technologies for digital upgrades/designs. A framework was proposed for this strategy, which aims to (1) provide a best-estimate, risk-informed capability to quantitatively and accurately estimate the risk impact of plant modernization, considering the introduction of high safety-significant safety-related (HSSSR) DI&C systems, (2) support and supplement existing risk-informed DI&C design guides by providing quantitative risk information and evidence, (3) offer a capability of design architecture evaluation of various DI&C systems, (4) assure the long-term safety and reliability of HSSSR DI&C systems, and (5) reduce uncertainty in costs and support integration of DI&C systems in the plant. To achieve these technical goals, the framework provides a means to address relevant technical issues by: (1) defining a risk-informed analysis process for DI&C upgrade, that integrates hazard analysis, reliability analysis, and consequence analysis, (2) applying risk-informed tools to address common cause failures (CCFs) and quantify corresponding failure probabilities for DI&C technologies, particularly software CCFs, (3) evaluating the impact of digital failures at the component level, system level, and plant level, and (4) providing insights and suggestions on designs to manage the risks, thus to support the development and deployment of advanced DI&C technologies on nuclear power plant (NPPs).

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Bayesian Cross-matching of High Proper-motion Stars in Gaia DR2 and Photometric Metallicities for ∼1.7 million K and M Dwarfs

We present a Bayesian method to cross-match 5,827,988 high proper-motion Gaia sources (μ > 40 mas yr{sup −1}) to various photometric surveys: Two Micron All Sky Survey, AllWISE data release from the Wide-field Infrared Explorer (WISE) mission, Galaxy Evolution Explorer, Radial Velocity Experiment, Sloan Digital Sky Survey, and Panoramic Survey Telescope and Rapid Response System (Pan-STARRS). To efficiently associate these objects across catalogs, we develop a technique that compares the multidimensional distribution of all sources in the vicinity of each Gaia star to a reference distribution of random field stars obtained by extracting all sources in a region on the sky displaced 2′. This offset preserves the local field stellar density and magnitude distribution, allowing us to characterize the frequency of chance alignments. The resulting catalog with Bayesian probabilities >95% has a marginally higher match rate than current internal Gaia data release 2 (DR2) matches for most catalogs. However, a significant improvement is found with Pan-STARRS, where ∼99.8% of the sample within the Pan-STARRS footprint is recovered, as compared to a low ∼20.8% in Gaia DR2. Using these results, we train a Gaussian process regressor to calibrate two photometric metallicity relationships. For dwarfs of 3500 < T{sub eff} < 5280 K, we use metallicity values of 4378 stars from the Apache Point Observatory Galactic Evolution Experiment and Hejazi et al. to calibrate the relationship, producing results with a 1σ precision of 0.12 dex and few systematic errors. We then indirectly infer the metallicity of 4018 stars with 2850 < T{sub eff} < 3500 K, which are wide companions of primaries whose metallicities are estimated with our first regressor, to produce a relationship with a 1σ precision of 0.21 dex and significant systematic errors. Additional work is needed to better remove unresolved binaries from this second sample to reduce these systematic errors.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Weldability of irradiated Stainless Steel 304 materials harvested from the National Research Universal (NRU) reactor

In this work, the weldability of the neutron-irradiated Stainless Steel 304 (SS 304) materials containing estimated helium levels of 12 to 44 appm has been studied through a Collaborative Research and Development Agreement between the Canadian Nuclear Laboratories and Oak Ridge National Laboratory. The work utilizes a SS 304 component harvested from the National Research Universal reactor. Laser beam welding was used to explore suitable welding conditions and associated parameters, to determine helium concentration limits for successful welding of irradiated stainless-steels. The experimental results show that maximum fusion zone void sizes, fusion zone void quantity and total heat affected zone helium-induced crack length generally increased with increasing He levels and effective weld heat input. A proper combination of welding parameters with low effective heat input, such as, laser power, weld travel speed and wire feed speed, were shown to improve the welding quality, and potentially reduce the length and probability of the formation of grain boundary helium-induced cracking and voids of irradiated SS 304 even with He concentrations as high as about 44 appm.

36 MATERIALS SCIENCE↗
DOE OSTI2023

AutoSourceID-Classifier: Star-galaxy classification using a convolutional neural network with spatial information

Aims.Traditional star-galaxy classification techniques often rely on feature estimation from catalogs, a process susceptible to introducing inaccuracies, thereby potentially jeopardizing the classification’s reliability. Certain galaxies, especially those not manifesting as extended sources, can be misclassified when their shape parameters and flux solely drive the inference. We aim to create a robust and accurate classification network for identifying stars and galaxies directly from astronomical images. Methods.The AutoSourceID-Classifier (ASID-C) algorithm developed for this work uses 32x32 pixel single filter band source cutouts generated by the previously developed AutoSourceID-Light (ASID-L) code. By leveraging convolutional neural networks (CNN) and additional information about the source position within the full-field image, ASID-C aims to accurately classify all stars and galaxies within a survey. Subsequently, we employed a modified Platt scaling calibration for the output of the CNN, ensuring that the derived probabilities were effectively calibrated, delivering precise and reliable results. Results.We show that ASID-C, trained on MeerLICHT telescope images and using the Dark Energy Camera Legacy Survey (DECaLS) morphological classification, is a robust classifier and outperforms similar codes such as SourceExtractor. To facilitate a rigorous comparison, we also trained an eXtreme Gradient Boosting (XGBoost) model on tabular features extracted by SourceExtractor. While this XGBoost model approaches ASID-C in performance metrics, it does not offer the computational efficiency and reduced error propagation inherent in ASID-C’s direct image-based classification approach. ASID-C excels in low signal-to-noise ratio and crowded scenarios, potentially aiding in transient host identification and advancing deep-sky astronomy.

Astronomy & Astrophysics↗
DOE OSTI2024

Aging of Limestones and Silane–Siloxane-Based Protective Hydrophobics: The Impact of Heating–Cooling and Freeze–Thaw Cycles

Stones are traditionally used in construction and architectural applications as building elements due to their aesthetic and technical/structural performance. Like other environmental factors (rain, humidity, moisture, salt presence, biological activity, etc.), heating–cooling and freeze–thaw cycles significantly threaten the longevity of stone materials. Hence, considering the socio-economic and cultural value of stones, preventive actions such as hydrophobic coatings are applied to prevent or mitigate damage. The scope of this study is the performance assessment of limestones with different characteristics and the efficiency of various commercial silane/siloxane-based hydrophobic coatings when exposed to thermal variation and freeze–thaw. For that purpose, the standards EN 14066:2013 (determination of resistance to aging by thermal shock) and EN 12371:2010 (determination of frost resistance) were followed. Open porosity and static contact angles were estimated to assess the stone durability and water protection capabilities of the hydrophobics. Additionally, sound speed propagation velocity, quality of building material index, elastic modulus and flexural strength were measured to evaluate the variation of mechanical properties. Static contact angle revealed that the coatings maintained an efficient level of hydrophobicity even after thermal-shock and freeze–thaw weathering tests. The study also revealed a critical interaction between freeze–thaw cycles, hydrophobic coatings and structural integrity of the stones, mostly on more porous ones. When they are subjected to harsh environmental conditions, untreated porous limestones keep structural cohesion, allowing for the natural absorption and release of water during freezing and thawing. On the contrary, when limestones are treated, the hydrophobic coatings can moderately obstruct the water release due to the partial saturation of the porous framework by the products. It also probably resulted from the different mechanical behavior between the inner matrix and layer of stone coated, resulting in a premature breakout and mechanical damage of the stone.

Lisci, Carla (ORCID:000000032179130X)↗
DOE OSTI2022

Parameter inference with non-linear galaxy clustering: accounting for theoretical uncertainties

ABSTRACT We implement euclidemulator (version 1), an emulator for the non-linear correction of the matter power spectrum, into the Markov chain Monte Carlo forecasting code montepython. We compare the performance of halofit, hmcode, and euclidemulator1, both at the level of power spectrum prediction and at the level of posterior probability distributions of the cosmological parameters, for different cosmological models and different galaxy power spectrum wavenumber cut-offs. We confirm that the choice of the power spectrum predictor has a non-negligible effect on the computed sensitivities when doing cosmological parameter forecasting, even for a conservative wavenumber cut-off of $0.2\, h\, {\rm Mpc}^{-1}$. We find that euclidemulator1 is on average up to 17 per cent more sensitive to the cosmological parameters than the other two codes, with the most significant improvements being for the Hubble parameter of up to 42 per cent and the equation of state of dark energy of up to 26 per cent, depending on the case. In addition, we point out that the choice of the power spectrum predictor contributes to the risk of computing a significantly biased mean cosmology when doing parameter estimations. For the four tested scenarios we find biases, averaged over the cosmological parameters, of between 0.5σ and 2σ (from below 1σ up to 6σ for individual parameters). This paper provides a proof of concept that this risk can be mitigated by taking a well-tailored theoretical uncertainty into account as this allows to reduce the bias by a factor of 2 to 5, depending on the case under consideration, while keeping posterior credibility contours small: the standard deviations are amplified by a factor of ≤1.4 in all cases.

Astronomy & Astrophysics↗
DOE OSTI2021

Sparse Bayesian mass mapping with uncertainties: hypothesis testing of structure

ABSTRACT A crucial aspect of mass mapping, via weak lensing, is quantification of the uncertainty introduced during the reconstruction process. Properly accounting for these errors has been largely ignored to date. We present a new method to reconstruct maximum a posteriori (MAP) convergence maps by formulating an unconstrained Bayesian inference problem with Laplace-type l1-norm sparsity-promoting priors, which we solve via convex optimization. Approaching mass mapping in this manner allows us to exploit recent developments in probability concentration theory to infer theoretically conservative uncertainties for our MAP reconstructions, without relying on assumptions of Gaussianity. For the first time, these methods allow us to perform hypothesis testing of structure, from which it is possible to distinguish between physical objects and artefacts of the reconstruction. Here, we present this new formalism, and demonstrate the method on simulations, before applying the developed formalism to two observational data sets of the Abell 520 cluster. Initial reconstructions of the Abell 520 catalogues reported the detection of an anomalous ‘dark core’ – an overdense region with no optical counterpart – which was taken to be evidence for self-interacting dark matter. In our Bayesian framework, it is found that neither Abell 520 data set can conclusively determine the physicality of such dark cores at $99{{\ \rm per\ cent}}$ confidence. However, in both cases the recovered MAP estimators are consistent with both sets of data.

Price, M. A.↗
DOE OSTI2023

The MIGDAL experiment: Measuring a rare atomic process to aid the search for dark matter

We present the Migdal In Galactic Dark mAtter expLoration (MIGDAL) experiment aiming at the unambiguous observation and study of the so-called Migdal effect induced by fast-neutron scattering. It is hoped that this elusive atomic process can be exploited to enhance the reach of direct dark matter search experiments to lower masses, but it is still lacking experimental confirmation. Our goal is to detect the predicted atomic electron emission which is thought to accompany nuclear scattering with low, but calculable, probability, by deploying an Optical Time Projection Chamber filled with a low-pressure gas based on CF 4 . Initially, pure CF 4 will be used, and then in mixtures containing other elements employed by leading dark matter search technologies — including noble species, plus Si and Ge. High resolution track images generated by a Gas Electron Multiplier stack, together with timing information from scintillation and ionisation readout, will be used for 3D reconstruction of the characteristic event topology expected for this process — an arrangement of two tracks sharing a common vertex, with one belonging to a Migdal electron and the other to a nuclear recoil. Different energy-loss rate distributions along both tracks will be used as a powerful discrimination tool against background events. In this article we present the design of the experiment, informed by extensive particle and track simulations and detailed estimations of signal and background rates. In pure CF 4 we expect to observe 8.9 (29.3) Migdal events per calendar day of exposure to an intense D–D (D–T) neutron generator beam at the NILE facility located at the Rutherford Appleton Laboratory (UK). With our nominal assumptions, 5σ median discovery significance can be achieved in under one day with either generator.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Red supergiant candidates for multimessenger monitoring of the next Galactic supernova

ABSTRACT We compile a catalogue of 578 highly probable and 62 likely red supergiants (RSGs) of the Milky Way, which represents the largest list of Galactic RSG candidates designed for continuous follow-up efforts to date. We match distances measured by Gaia DR3, 2MASS photometry, and a 3D Galactic dust map to obtain luminous bright late-type stars. Determining the stars’ bolometric luminosities and effective temperatures, we compare to Geneva stellar evolution tracks to determine likely RSG candidates, and quantify contamination using a catalogue of Galactic AGB in the same luminosity-temperature space. We add details for common or interesting characteristics of RSG, such as multistar system membership, variability, and classification as a runaway. As potential future core-collapse supernova progenitors, we study the ability of the catalogue to inform the Supernova Early Warning System (SNEWS) coincidence network made to automate pointing, and show that for 3D position estimates made possible by neutrinos, the number of progenitor candidates can be significantly reduced, improving our ability to observe the progenitor pre-explosion and the early phases of core-collapse supernovae.

Astronomy & Astrophysics↗
DOE OSTI2022

Integration of Omics into a New Comprehensive Rate Law for Competitive Terminal Electron-Accepting Processes in Reactive Transport Models: Application to N, Fe, S, and Contaminant Transformations in Stream and Wetland Sediments

Surface waters represent important sources of alternate energy and drinking water in the United States, and characterizing the biogeochemical processes that affect surface water quality is relevant to the DOE-BER mission. Sediment biogeochemical processes regulate the release of carbon (C), nutrients, and contaminants to surface waters and thus influence water quality. Sediment biogeochemical processes are dynamic and affected by the deposition and remobilization of solid material and changes in environmental conditions driven by water discharge variations. Wetlands are important natural filters of surface waters which may either trap, metabolize, or mobilize nutrients and contaminants. Despite their importance, biogeochemical processes regulating nutrient and contaminant release and C transformation in stream and wetland sediments cannot be predicted accurately by current mathematical models. These reactive transport models largely rely on detectable changes in geochemical conditions to activate metabolic processes, do not accurately account for the competition between microbial processes, and poorly constrain effects of hydrological perturbations on biogeochemical processes. In this BER-SBR exploratory project, metagenomic and geochemical signatures were combined to identify microbially-mediated redox processes in anaerobic stream and wetland sediments from the Savannah River Site (SRS, ANL SFA) and East Fork Poplar Creek (EFPC, ORNL SFA) that play important roles in C, uranium (U), and mercury (Hg) transformations. In addition, sediment incubations were conducted to examine the competition between anaerobic respiration processes Finally, new rate laws were developed for reactive transport models that rely on complementary metagenomic and geochemical signatures to identify the underlying anaerobic microbial processes in stream and wetland sediments, describe the competition between the dominant metabolic processes involved in nutrient release and U and Hg mobilization, and more accurately quantify carbon transformation and the response of microbial processes to changes in redox conditions associated with hydrological forcing. These rate laws were optimized in batch reactors with SRS wetland sediments, where iron and sulfate reduction dominate. Anaerobic carbon remineralization processes followed the expected thermodynamic sequence of microbial respiration with depth in both sediments, except that geochemical signals indicated that sulfate reduction was inactive in EFPC sediments and moderate in SRS wetland sediments. Estimates indicated that microbial iron reduction contributed to at least half of the production of reduced iron in these sediments. Incubations demonstrated that nitrate reduction, denitrification, and dissimilatory nitrate reduction to ammonium were active in the natural EFPC sediment and activated upon nitrate amendment in these nitrate-rich sediments. In turn, these processes were outcompeted by the addition of either iron oxides or sulfate as alternative terminal electron acceptors. Although geochemical products of sulfate reduction were not detected in the incubations, the abundance of sulfate reduction genes increased with depth in the sediment and was equally more pronounced in treatments amended with either iron oxides or sulfate. Simultaneously, anaerobic sulfide oxidizing bacteria coupling sulfide oxidation to DNRA (and not conventional denitrification) were apparently enriched over time, regardless of the treatments. These findings indicate that sulfate reduction is important in freshwater stream sediments and probably catalyzed by a cryptic sulfur cycle involving nitrogen species, in which the sulfur products from sulfate reduction are immediately removed by side reactions and not detectable by geochemical measurements alone. Similar experiments in SRS sediments, however, demonstrated little interaction between nitrogen and sulfur cycling microorganisms. Sulfate reduction was impacted by the addition of more thermodynamically favorable electron acceptors, suggesting either that iron-reducing microorganisms outcompeted sulfate-reducing microorganisms for organic substrate, depleting the stock of electron donor available for sulfate reduction, or that the cryptic sulfur cycle was shunted by the precipitation of FeS generated as a result of the abiotic reduction of iron oxides by dissolved sulfide. A diagnostic modeling exercise was conducted to further investigate the competition between terminal electron accepting processes. As conventional kinetic models typically do not account for cryptic cycles and used inaccurate formulations to describe competition between microbial communities, new metabolic rate laws were developed that explicitly express the electron acceptor-specific enzyme of each energetically favorable metabolic process based on gene abundance detected in the incubations. The model was tested with the sediment slurry incubation data to determine whether substrate competition could explain the decrease in sulfate reduction observed in the presence of iron oxide competitor. The model was able to reproduce geochemical concentrations really well in each treatment once the model was calibrated with the unamended control, suggesting that microbial competition was indeed driven by thermodynamic considerations. Overall, carbon remineralization processes and rates will be reproduced much more realistically with the new metabolic rate laws.

54 ENVIRONMENTAL SCIENCES↗