Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Confidence interval”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A Standard Siren Measurement of the Hubble Constant Using Gravitational-wave Events from the First Three LIGO/Virgo Observing Runs and the DESI Legacy Survey

We present a new constraint on the Hubble constant H 0 using a sample of well-localized gravitational-wave (GW) events detected during the first three LIGO/Virgo observing runs as dark standard sirens. In the case of dark standard sirens, a unique host galaxy is not identified, and the redshift information comes from the distribution of potential host galaxies. From the third LIGO/Virgo observing run detections, we add the asymmetric-mass binary black hole GW190412 and the high-confidence GW candidates S191204r, S200129m, and S200311bg to the sample of dark standard sirens analyzed in Palmese et al. Our sample contains the top 20% (based on localization) GW events and candidates to date with significant coverage by the Dark Energy Spectroscopic Instrument Legacy Survey. We combine the H 0 posterior for eight dark siren events, finding H 0 = $79.8^{+19.1}_{-12.8}$ km s -1 Mpc -1 (68% highest density interval) for a prior in H 0 uniform between [20, 140] km s -1 Mpc -1 . This result shows that a combination of eight well-localized dark sirens combined with an appropriate galaxy catalog is able to provide an H 0 constraint that is competitive (~20% versus 18% precision) with a single bright standard siren analysis (i.e., assuming the electromagnetic counterpart) using GW170817. When combining the posterior with that from GW170817, we obtain H 0 = $72.77^{+11.0}_{-7.55}$ km s -1 Mpc -1 . This result is broadly consistent with recent H 0 estimates from both the cosmic microwave background and supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Rare Earth Elements (REE) and Critical Minerals (CM) in Middle Pennsylvanian-Age Coals and Associated Sediments in the Central Appalachian Basin, Eastern U.S.A.

The Central Appalachian Basin (CAB) of Kentucky, Tennessee, West Virginia and Virginia has a long history of coal mining and oil and gas extraction that has empowered the regional and national economies, the development of infrastructure, and a highly trained energy resources work force. As our societal demands for advanced technologies have rapidly increased in recent years, coal-related materials are viewed as an important new unconventional domestic source of critical minerals (CM) that are required for telecommunications, aerospace and transportation industries, electronics, the transition to low-carbon emissions energy production, and many consumer products. Coal-related materials encompass coal, associated sediments, coal mining waste materials, produced waters, and ash residues from coal-fired power plants. An important objective of the Evolve Central Appalachia Project (Evolve CAPP), sponsored by the U.S. Department of Energy (DOE) National Energy Technology Laboratory (NETL) is to assess the quantity and distribution of CM resources in the CAB region. The rare earth elements (REE) are considered highest priority, although other important CM such as niobium, gallium, and zirconium are known to occur in the Middle Pennsylvanian-age coals and sediments. Working with coal industry partners who provided access to drill cores, coal-related sediments and waste materials, over 600 samples have been collected for laboratory analysis, and over 730 materials have been scanned using portable x-ray fluorescence (pXRF) equipment. The application of pXRF provides the means for real-time semiquantitative analysis of CM content at very close spacings (typically 2-3 inch intervals) along drill core and in-situ channel samples that span the roof, coal seam, and floor rock. The comparison of pXRF geochemical data with laboratory results, geologic data, and downhole spectral gamma logs can provide high resolution input to lithologic and depositional models for future CM resource evaluations. The preliminary findings show that pXRF is capable of accurately measuring low concentrations of many of the CM with a high level of confidence (Ba, Cr, Ga, K, Nb, Rb, Sr, Th, Y), whereas for others (La, Ce, Co, Mn, Nd, Ni, Sc, Ti, V) the detection limits are very high or spectral interferences increase the uncertainty. Notably, the mean abundances of Y (34 ppm), La (86 ppm), Ga (40 ppm), and V (146 ppm) in coal underclays in the CAB region are up to 7X enriched compared with the overlying coal. These values also exceed the reported concentrations in published reference materials for upper continental crustal rocks (Rudnick and Gao, 2003), North American Shale Composite (Gromet et al., 1984), and North American coal (Finkelman, 1993). The pXRF data are in part verified by laboratory results that indicate the mean Y abundance (36 ppm) is highly correlated (R2 = 0.807) with total REE (ΣREE). The correlation is even higher (R2 = 0.957) with heavy REE (ΣHREE). Applying these correlations to the pXRF data for the coal underclays, the mean estimated values for ΣREE+Y and ΣHREE+Y are 270 ppm and 58 ppm, respectively. Although these average values are not considered high, the range of Y measured by pXRF in the coal underclays extended as high as 114 ppm, which would suggest ΣREE+Y equal to 847 ppm. The mean abundances of Zr (192 ppm) and Th (21 ppm) in coal underclays are also enriched compared with the overlying coal and these results likely reflect the presence of resistant detrital heavy minerals such as monazite, xenotime, and zircon in the underclay matrix. Several of the profiled coal seams and associated wall rocks contained thin volcanic ash layers up to 4-5 inches in thickness. The extent to which these ash fall layers provided a source for CM under the paleoenvironmental conditions that resulted in coal deposits in the CAB remains to be fully studied. Continuing investigations in the Evolve CAPP study area will include laboratory determinations of mineralogic and clay compositions, and evaluations of CM geochemical mobility in the coal and coal underclays.

Lassetter, Billy↗
DOE OSTI2024

Combination of Searches for Higgs Boson Pair Production in pp Collisions at $\sqrt{s}$ = 13 TeV with the ATLAS Detector

This Letter presents results from a combination of searches for Higgs boson pair production using 126–140 fb−1 of proton-proton collision data at √𝑠 =13 TeV recorded with the ATLAS detector. At 95% confidence level (CL), the upper limit on the production rate is 2.9 times the standard model (SM) prediction, with an expected limit of 2.4 assuming no Higgs boson pair production. Constraints on the Higgs boson self-coupling modifier 𝜅 𝜆 = 𝜆 𝐻⁡𝐻⁡𝐻 /𝜆$^{SM}_{𝐻⁡𝐻⁡𝐻}$, and the quartic 𝐻⁡𝐻⁡𝑉⁢𝑉 coupling modifier 𝜅 2⁢𝑉 = 𝑔 𝐻⁡𝐻⁡𝑉⁢𝑉⁡ /𝑔$^{SM}_{𝐻⁡𝐻⁡𝑉⁢𝑉}$, are derived individually, fixing the other parameter to its SM value. The observed 95% CL intervals are −1.2<𝜅 𝜆 < 7.2 and 0.6<𝜅 2⁢𝑉 < 1.5, respectively, while the expected intervals are −1.6<𝜅 𝜆 < 7.2 and 0.4<𝜅 2⁢𝑉 < 1.6 in the SM case. Constraints obtained for several interaction parameters within Higgs effective field theory are the strongest to date, offering insights into potential deviations from SM predictions.

Hadron colliders↗
DOE OSTI2025

Performance Evaluation of Intelligent Solar Control Software Through Hardware-in-the-Loop (CRADA Final Report)

Recent research has highlighted the potential for solar to act as a zero-marginal-cost and zero-emission flexibility resource on the bulk power system when operated with advanced control systems. To increase the performance of these systems, leading technologies, including machine learning (ML) and hierarchical inverter set point allocation, have been developed by Latimer Controls, Inc. to estimate the headroom of large PV plants for grid operation and control; however, these technologies lack comprehensive validation under real-world application scenarios. Latimer Controls, Inc. received two voucher awards for research at a national laboratory from the Department of Energy American Made Solar Prize Round 6. The National Renewable Energy Laboratory (NREL) was selected to collaborate with Latimer staff to conduct a performance evaluation of Latimer PV control software. The NREL team will develop a hardware-in-the-loop (HIL) testbed to perform testing and validation of the Latimer PV control technology in a de-risked yet realistic testbed environment. Latimer and NREL worked together to analyze the test data, draw conclusions from the results, and disseminate the resulting scientific findings. In this CRADA work, we propose to test and validate the real-world application of the Latimer Control solution in an HIL environment. We evaluate the performance of different flexible solar technologies in responding to automatic generation control signals in a closed-loop fashion. In particular, a data-driven potential high limit (PHL) estimation is developed for large solar plants to accurately estimate their headroom so that they have fast and short-time regulation and control capability to participate in grid services and respond to grid signals in real time (e.g., AGC). This PHL estimation algorithm is embedded in a hardware power plant controller (PPC) and tested with an IEEE-39 bus system model developed in RTDS. To account for the varying cloud conditions and diverse inverter dispatches, we developed a 135-MW PV plant with detailed modeling of 27 individual PV modules and inverters using RTDS. The real-world communications used in such big plants, such as ModBus TCP/IP for inverter level and DNP3 for plant level, were developed to emulate the real-world applications in big PV plants. The ML-based PHL estimation method is tested under nine separate weather scenarios against the ‘reference-control’ solution, hereafter referred to as the baseline solution. The baseline method reserves a subset of inverters (reference group) to operate at their PHL at all times and dispatches only the remaining inverters (control group) at curtailed levels to fulfill the flexibility need. Despite being successfully piloted by NREL in California in 2017 and Chile in 2020, there exist two gaps in the state of the art to fully unlock the flexibility of PV plants: a. There is a trade-off between the PHL estimation accuracy and the flexibility range. b. There lacks granularity in the PHL estimation to capture the variation across inverters. The Latimer solution seeks to address these gaps by applying machine learning methods to improve PHL estimation accuracy while accounting for variability at every inverter. Performance metrics were taken from the 2023 Georgia Power CARES utility-scale RFP. The results demonstrate that the ML-based approach outperforms the traditional baseline method in PHL estimation accuracy for 7 of 9 scenarios. The average PHL error across the nine scenarios was 7.40% for the ML-based method, 2.06% less than the 9.46% PHL error average across scenarios that was exhibited by the baseline method. Additionally, the PHL error was below 5% for at least 95% of the testing interval for 3 of 9 tested intervals with the ML approach, whereas it did not achieve this metric for any of the baseline tests. Overall, simulation results indicate the superior performance of an ML-based approach compared to the conventional baseline reference-control approach, showcasing its potential to support grid stability and operational efficiency. This laboratory HIL testing using real PPC, representative power system simulation models in real-time with detailed PV plant and inverter models, and real-world communication protocols gives us confidence that this machine learning based PHL estimation algorithm works well in the hardware PPC and therefore de-risks future field commissioning. The end goal of this project is to advance grid technology to address the grid operation challenges brought by solar plant’s variability and uncertainties in power generation.

14 SOLAR ENERGY↗
DOE OSTI2024

Challenges of COVID-19 Case Forecasting in the US, 2020–2021

During the COVID-19 pandemic, forecasting COVID-19 trends to support planning and response was a priority for scientists and decision makers alike. In the United States, COVID-19 forecasting was coordinated by a large group of universities, companies, and government entities led by the Centers for Disease Control and Prevention and the US COVID-19 Forecast Hub ( https://covid19forecasthub.org ). We evaluated approximately 9.7 million forecasts of weekly state-level COVID-19 cases for predictions 1–4 weeks into the future submitted by 24 teams from August 2020 to December 2021. We assessed coverage of central prediction intervals and weighted interval scores (WIS), adjusting for missing forecasts relative to a baseline forecast, and used a Gaussian generalized estimating equation (GEE) model to evaluate differences in skill across epidemic phases that were defined by the effective reproduction number. Overall, we found high variation in skill across individual models, with ensemble-based forecasts outperforming other approaches. Forecast skill relative to the baseline was generally higher for larger jurisdictions (e.g., states compared to counties). Over time, forecasts generally performed worst in periods of rapid changes in reported cases (either in increasing or decreasing epidemic phases) with 95% prediction interval coverage dropping below 50% during the growth phases of the winter 2020, Delta, and Omicron waves. Ideally, case forecasts could serve as a leading indicator of changes in transmission dynamics. However, while most COVID-19 case forecasts outperformed a naïve baseline model, even the most accurate case forecasts were unreliable in key phases. Further research could improve forecasts of leading indicators, like COVID-19 cases, by leveraging additional real-time data, addressing performance across phases, improving the characterization of forecast confidence, and ensuring that forecasts were coherent across spatial scales. In the meantime, it is critical for forecast users to appreciate current limitations and use a broad set of indicators to inform pandemic-related decision making.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Searches for exclusive Higgs and Z boson decays into a vector quarkonium state and a photon using 139 fb$^{-1}$ of ATLAS $\sqrt{s}=13$ TeV proton–proton collision data

Searches for the exclusive decays of Higgs and Z bosons into a vector quarkonium state and a photon are performed in the μ + μ - γ final state with a proton– proton collision data sample corresponding to an integrated luminosity of 139 fb -1 collected at $\sqrt{s}=13$ TeV with the ATLAS detector at the CERN Large Hadron Collider. The observed data are compatible with the expected back grounds. The 95% confidence-level upper limits on the branching fractions of the Higgs boson decays into J/ψ γ , ψ(2S) γ , and $Υ$(1S, 2S, 3S) γ are found to be 2.0 × 10 -4 , 10.5×10 -4 , and (2.5, 4.2, 3.4)×10 -4 , respectively, assuming Standard Model production of the Higgs boson. The corresponding 95% CL upper limits on the branching fractions of the Z boson decays are 1.2 × 10 -6 , 2.4 × 10 -6 , and (1.1, 1.3, 2.4) × 10 -6 . An observed 95% CL interval of (-133, 175) is obtained for the κ c /κ γ ratio of Higgs boson coupling modifiers, and a 95% CL interval of (-37, 40) is obtained for κ b /κ γ .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Autonomous monitoring of algal biomass: Success stories and lessons learned from long-term field deployment

Autonomous, high-frequency monitoring of outdoor algal ponds is needed to quantify biomass productivity and detect culture decline in environments prone to contamination, grazers, and variable operating conditions. We report successes and lessons learned in translating a laboratory spectroradiometric monitoring approach to a multi-year autonomous field deployment at the Arizona Center for Algae Technology and Innovation (AzCATI). The system measures spectrally resolved pond reflectance by ratioing upwelling radiance from each raceway to simultaneous downwelling sky irradiance using fiber-coupled spectrometers. A physics-based reflectance model (ASHARP) is fit to each spectrum pair to estimate optical parameters, including a biomass-proxy coefficient (C a ) which enables near-real-time tracking of biomass accumulation and culture state at 2–5 min intervals. From May 2022 through September 2025 the platform operated continuously while scaling from two to six raceway ponds. Several strains of algae were monitored successfully, including the high productivity Tetraselmis striata and Picochlorum celeri. Transitioning data acquisition from a Windows laptop to a Raspberry Pi improved uptime from 57% (2022) to ~89% (2024–2025) and enabled routine real-time analysis. Further, we converted relative biomass estimates to absolute ash-free dry weight (AFDW) using experimentally-derived calibrations, providing field-relevant biomass predictions with conservative confidence bounds. These results demonstrate the feasibility of long-term, autonomous optical monitoring for well-mixed open-raceway algal cultivation and provide practical guidance for reliable field operation and scaling.

Katinas, Christopher Michael [Sandia National Labo↗
DOE OSTI2023

Uncertainty quantification of machine learning models to improve streamflow prediction under changing climate and environmental conditions

Machine learning (ML) models, and Long Short-Term Memory (LSTM) networks in particular, have demonstrated remarkable performance in streamflow prediction and are increasingly being used by the hydrological research community. However, most of these applications do not include uncertainty quantification (UQ). ML models are data driven and can suffer from large extrapolation errors when applied to changing climate/environmental conditions. UQ is required to quantify the influence of data noises on model predictions and avoid overconfident projections in extrapolation. In this work, we integrate a novel UQ method, called PI3NN, with LSTM networks for streamflow prediction. PI3NN calculates Prediction Intervals by training 3 Neural Networks. It can precisely quantify the predictive uncertainty caused by the data noise and identify out-of-distribution (OOD) data in a non-stationary condition to avoid overconfident predictions. We apply the PI3NN-LSTM method in the snow-dominant East River Watershed in the western US and in the rain-driven Walker Branch Watershed in the southeastern US. Results indicate that for the prediction data which have similar features as the training data, PI3NN precisely quantifies the predictive uncertainty with the desired confidence level; and for the OOD data where the LSTM network fails to make accurate predictions, PI3NN produces a reasonably large uncertainty indicating that the results are not trustworthy and should avoid overconfidence. PI3NN is computationally efficient, robust in performance, and generalizable to various network structures and data with no distributional assumptions. It can be broadly applied in ML-based hydrological simulations for credible prediction.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Feature review of photovoltaic modeling software utilizing blind performance assessment

While confidence in photovoltaic (PV) modeling software has always been essential, the rapid pace of new PV plant developments makes accuracy and credibility more critical than ever. Independent assessments, particularly through blind modeling comparisons, are therefore necessary to ensure unbiased benchmarking across PV modeling software. Previous studies have been limited by a narrow range of models compared, anonymized results, or system size. This study presents results from the first-ever onymous blind modeling comparison, evaluated using both lab- and utility-scale fixed-tilt, monofacial, south-facing systems at sub-hourly time intervals. Seven commercially used PV software tools were compared: 3E SynaptiQ, PlantPredict, PVsyst, RatedPower, SAM, SolarFarmer, and Solargis Evaluate. Predictions were submitted directly by software representatives, providing unique insights into each software’s implementation and resulting prediction behavior. Notable features, including plane-of-array (POA) transposition model, module temperature model, shading model, and performance model were analyzed and compared. Four summary tables compile these features of the software, serving as a resource to help users understand the methodological differences and select the most suitable software for their applications. The software tools show deviations from mean error in annual yield up to 2.5 % in the lab-scale system, increasing to 6.0 % for the utility-scale system. These differences arise from a combination of user decisions and the inherent behavior of the software, indicating the need for continuous and rigorous validation of modeling methods using these software tools against complex, real-world systems.

14 SOLAR ENERGY↗
DOE OSTI2023

Population of Merging Compact Binaries Inferred Using Gravitational Waves through GWTC-3

We report on the population properties of compact binary mergers inferred from gravitational-wave observations of these systems during the first three LIGO-Virgo observing runs. The Gravitational-Wave Transient Catalog 3 (GWTC-3) contains signals consistent with three classes of binary mergers: binary black hole, binary neutron star, and neutron star–black hole mergers. We infer the binary neutron star merger rate to be between 10 and 1700 Gpc − 3 yr − 1 and the neutron star–black hole merger rate to be between 7.8 and 140 Gpc − 3 yr − 1 , assuming a constant rate density in the comoving frame and taking the union of 90% credible intervals for methods used in this work. We infer the binary black hole merger rate, allowing for evolution with redshift, to be between 17.9 and 44 Gpc − 3 yr − 1 at a fiducial redshift ( z = 0.2 ). The rate of binary black hole mergers is observed to increase with redshift at a rate proportional to ( 1 + z ) κ with κ = 2. 9 − 1.8 + 1.7 for z ≲ 1 . Using both binary neutron star and neutron star–black hole binaries, we obtain a broad, relatively flat neutron star mass distribution extending from 1.2 − 0.2 + 0.1 to 2.0 − 0.3 + 0.3 M ⊙ . We confidently determine that the merger rate as a function of mass sharply declines after the expected maximum neutron star mass, but cannot yet confirm or rule out the existence of a lower mass gap between neutron stars and black holes. We also find the binary black hole mass distribution has localized over- and underdensities relative to a power-law distribution, with peaks emerging at chirp masses of 8.3 − 0.5 + 0.3 and 27.9 − 1.8 + 1.9 M ⊙ . While we continue to find that the mass distribution of a binary’s more massive component strongly decreases as a function of primary mass, we observe no evidence of a strongly suppressed merger rate above approximately 60 M ⊙ , which would indicate the presence of a upper mass gap. Observed black hole spins are small, with half of spin magnitudes below χ i ≈ 0.25 . While the majority of spins are preferentially aligned with the orbital angular momentum, we infer evidence of antialigned spins among the binary population. We observe an increase in spin magnitude for systems with more unequal-mass ratio. We also observe evidence of misalignment of spins relative to the orbital angular momentum. Published by the American Physical Society 2023

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Utah FORGE: Well 16B(78)-32 Drill Core Fracture Analysis Images and Data

This dataset contains drilling core data from well 16B(78)-32, including PDF documents with flattened core images annotated by feature type and core interval, as well as spreadsheets detailing feature morphologies by depth, planar feature measurements, and planar feature orientations rotated to in situ conditions. Core was recovered from three intervals, one per stimulation stage, in the crystalline rocks affected by the stimulation of well 16A(78)-32. Seven core runs were conducted, yielding 135.8 feet of recovered core. Features in the core were categorized into planar fractures, semi-planar fractures, unbroken mineralized fractures, rough fractures, curviplanar fractures, concave-convex surfaces, and planar compositional features such as mylonite or dike-like structures. Planar features were measured while the core was positioned horizontally, with the core axis aligned to a downhole azimuth of 42 degrees. Planar core measurements from stimulations 2 and 3 that could be confidently correlated with FMI data were rotated to in situ orientations. This was done by rotating the planes along vertical and horizontal axes to match the azimuth and inclination data recorded in the directional survey of well 16B(78)-32, as well as applying an axial rotation to resemble the fracture orientations observed in the FMI log at corresponding depths. Coherent sets of planar fracture measurements were made by aligning the core within each 3-foot section of the dissected core barrel, and between adjacent 3-foot sections within a core run by matching rock fabrics, saw cuts and/or tool marks. Where coherent fracture measurements could not be made within a core run, data sets are denoted by a subscript (i.e. 2-Ta and 2-Tb both come from tangent core run number 2).

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Improving deep learning-based protein distance prediction in CASP14

Abstract Motivation Accurate prediction of residue–residue distances is important for protein structure prediction. We developed several protein distance predictors based on a deep learning distance prediction method and blindly tested them in the 14th Critical Assessment of Protein Structure Prediction (CASP14). The prediction method uses deep residual neural networks with the channel-wise attention mechanism to classify the distance between every two residues into multiple distance intervals. The input features for the deep learning method include co-evolutionary features as well as other sequence-based features derived from multiple sequence alignments (MSAs). Three alignment methods are used with multiple protein sequence/profile databases to generate MSAs for input feature generation. Based on different configurations and training strategies of the deep learning method, five MULTICOM distance predictors were created to participate in the CASP14 experiment. Results Benchmarked on 37 hard CASP14 domains, the best performing MULTICOM predictor is ranked 5th out of 30 automated CASP14 distance prediction servers in terms of precision of top L/5 long-range contact predictions [i.e. classifying distances between two residues into two categories: in contact (<8 Angstrom) and not in contact otherwise] and performs better than the best CASP13 distance prediction method. The best performing MULTICOM predictor is also ranked 6th among automated server predictors in classifying inter-residue distances into 10 distance intervals defined by CASP14 according to the precision of distance classification. The results show that the quality and depth of MSAs depend on alignment methods and sequence databases and have a significant impact on the accuracy of distance prediction. Using larger training datasets and multiple complementary features improves prediction accuracy. However, the number of effective sequences in MSAs is only a weak indicator of the quality of MSAs and the accuracy of predicted distance maps. In contrast, there is a strong correlation between the accuracy of contact/distance predictions and the average probability of the predicted contacts, which can therefore be more effectively used to estimate the confidence of distance predictions and select predicted distance maps. Availability and implementation The software package, source code and data of DeepDist2 are freely available at https://github.com/multicom-toolbox/deepdist and https://zenodo.org/record/4712084#.YIIM13VKhQM. Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Search for photons above 10 18 eV by simultaneously measuring the atmospheric depth and the muon content of air showers at the Pierre Auger Observatory

The Pierre Auger Observatory is the most sensitive instrument to detect photons with energies above 1 0 17 eV . It measures extensive air showers generated by ultrahigh energy cosmic rays using a hybrid technique that exploits the combination of a fluorescence detector with a ground array of particle detectors. The signatures of a photon-induced air shower are a larger atmospheric depth of the shower maximum ( X max ) and a steeper lateral distribution function, along with a lower number of muons with respect to the bulk of hadron-induced cascades. In this work, a new analysis technique in the energy interval between 1 and 30 EeV ( 1 EeV = 1 0 18 eV ) has been developed by combining the fluorescence detector-based measurement of X max with the specific features of the surface detector signal through a parameter related to the air shower muon content, derived from the universality of the air shower development. No evidence of a statistically significant signal due to photon primaries was found using data collected in about 12 years of operation. Thus, upper bounds to the integral photon flux have been set using a detailed calculation of the detector exposure, in combination with a data-driven background estimation. The derived 95% confidence level upper limits are 0.0403, 0.01113, 0.0035, 0.0023, and 0.0021 km − 2 sr − 1 yr − 1 above 1, 2, 3, 5, and 10 EeV, respectively, leading to the most stringent upper limits on the photon flux in the EeV range. Compared with past results, the upper limits were improved by about 40% for the lowest energy threshold and by a factor 3 above 3 EeV, where no candidates were found and the expected background is negligible. The presented limits can be used to probe the assumptions on chemical composition of ultrahigh energy cosmic rays and allow for the constraint of the mass and lifetime phase space of super-heavy dark matter particles. Published by the American Physical Society 2024

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Search for resonant and nonresonant production of pairs of dijet resonances in proton-proton collisions at $ \sqrt{s} $ = 13 TeV

A search for pairs of dijet resonances with the same mass is conducted in final states with at least four jets. Results are presented separately for the case where the four jet production proceeds via an intermediate resonant state and for nonresonant production. The search uses a data sample corresponding to an integrated luminosity of 138 fb$^{−1}$ collected by the CMS detector in proton-proton collisions at $ \sqrt{s} $ = 13 TeV. Model-independent limits, at 95% confidence level, are reported on the production cross section of four-jet and dijet resonances. These first LHC limits on resonant pair production of dijet resonances via high mass intermediate states are applied to a signal model of diquarks that decay into pairs of vector-like quarks, excluding diquark masses below 7.6 TeV for a particular model scenario. There are two events in the tails of the distributions, each with a four-jet mass of 8 TeV and an average dijet mass of 2 TeV, resulting in local and global significances of 3.9 and 1.6 standard deviations, respectively, if interpreted as a signal. The nonresonant search excludes pair production of top squarks with masses between 0.50 TeV to 0.77 TeV, with the exception of a small interval between 0.52 and 0.58 TeV, for supersymmetric R-parity-violating decays to quark pairs, significantly extending previous limits. Here, the most significant excess above the predicted background occurs at an average dijet mass of 0.95 TeV, for which the local and global significances are 3.6 and 2.5 standard deviations, respectively.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Identification of Faults Susceptible to Induced Seismicity (Final Report)

Central to the work documented in this report is the capability of geocellular models to represent the geologic conceptual model updated with fault identification from machine learning and joint inversion modeling of microseismic data measured and recorded as a consequence of CO 2 injection at a field demonstration site: the Illinois Basin - Decatur Project (IBDP). This work required seven unique geocellular models with 100s of simulated variations to gain a very high degree of confidence in the identification of geologic features present that contributed to induced microseismicity at IBDP. All forward modeling: pressure modeling, stress modeling, and seismic modeling used the same geologic conceptual model and representations of that model at different scales. The pressure modeling and poroelastic modeling created “snapshots” of pore pressure and stress field changes at different times during CO 2 injection, in which microseismic events were clustered (in time). These pressure and stress snapshots, within the framework and architecture of the geologic conceptual model via the geocellular model, informed the single fault and fault network models to ascertain the likelihood of fault movement (seismic or aseismic). The outcomes of the pressure, stress, and fault/fault network (seismic) modeling confirmed that the faults in the geologic conceptual model in Task 2 were likely the source of microseismic events measured at IBDP and acted as conduits for pressure to be transmitted from the injection interval into the Precambrian crystalline basement rock. This closely coordinated and integrated unique modeling approach was conducted to prove the viability of our proposed workflow 1) to better resolve crystalline basement faults, 2) detect subseismic faults that could be activated by injection, 3) increase the certainty in fault detection and their susceptibility to release seismic energy, and 4) understand transmission of pressure vertically from the well to the underlying fractured crystalline basement. The proposed methodology was effective in guiding an iterative process of calibrating forward modeling results based on similar geocellular models while honoring the geologic conceptual model (i.e., characterization data and knowledge of regional geology); this led to higher level of certainty in the identification of fault/faults zones to control seismicity and transmission of pressure to the regions of recorded and located injection induced seismicity.

58 GEOSCIENCES↗
DOE OSTI2025

Restoring Historic Forest Disturbance Frequency Would Partially Mitigate Droughts in the Central Sierra Nevada Mountains

Forest thinning and prescribed fire are expected to improve the climate resilience and water security of forests in the western U.S., but few studies have directly modeled the hydrological effects of multi-decadal landscape-scale forest disturbance. By updating a distributed process-based hydrological model (DHSVM) with vegetation maps from a distributed forest ecosystem model (LANDIS-II), we simulate the water resource impacts of forest management scenarios targeting partial or full restoration of the pre-colonial disturbance return interval in the central Sierra Nevada mountains. In a fully restored disturbance regime that includes fire, thinning, and insect mortality, reservoir inflow increases by 4%–9% total and 8%–14% in dry years. At sub-watershed scales (10–100 km2), thinning dense forests can increase streamflow by >20% in dry years. In a thinner forest, increased understory transpiration compensates for decreased overstory transpiration. Consequentially, 73% of streamflow gains are attributable to decreased overstory rain and snow interception loss. Thinner forests can increase headwater peak flows, but reservoir-scale peak flows are almost exclusively influenced by climate. Uncertainty in future precipitation causes high uncertainty in future water yield, but the additional water yield attributable to forest disturbance is about five times less sensitive to annual precipitation uncertainty. This partial decoupling of the streamflow disturbance response from annual precipitation makes disturbance especially valuable for water supply during dry years. Our study can increase confidence in the water resource benefits of restoring historic forest disturbance frequencies in the central Sierra Nevada mountains, and our modeling framework is widely applicable to other forested mountain landscapes.

Boardman, Eli N. [University of Nevada, Reno, NV (↗
DOE OSTI2025

Temporal relationships among lunar crustal rocks

Temporal relationships among the three most common suites of lunar crustal rocks have been investigated by obtaining new high precision ages on Felsic/Alkali-suite Quartz monzodiorite Clast B from breccia 15405 and Magnesian-suite norite 78235/6/8/55/56 and comparing them to previously dated ferroan anorthosite sample 60025. The weighted average age of 4337.19 ± 0.49 Ma of 15405 Clast B is defined by zircon U-Pb and Pb-Pb ages as well as mineral isochron Sm-Nd and Nd-Nd ages. It is identical to the weighted average age for Apollo 17 norite 78235/6/8/55/56 of 4334.1 ± 3.5 Ma which is defined by Pb-Pb ages measured on baddeleyites in this investigation and less precise Pb-Pb and Sm-Nd ages reported in the literature. Both ages are ∼ 25 Ma younger than the weighted average of Sm-Nd and Pb-Pb ages reported in the literature on ferroan anorthosite 60025 of 4359.3 ± 2.3 Ma. The fact that ages of all three samples are defined by multiple U-Pb, Pb-Pb, Sm-Nd, and 142 Nd- 143 Nd chronometers provide confidence that they record the igneous crystallization history of the samples and do not represent disturbances or mixing lines with no temporal significance. Here, the extent to which these three ages represent broader scale magmatism is difficult to evaluate. Nevertheless, the age defined for 15405 Clast B, 78235/6/8/55/56, and 60025 are contemporaneous with the peak of ages observed in detrital zircons from the Apollo 12, 14, 15, and 17 landing sites (4340 ± 20 Ma), a Mg-suite Sm-Nd whole rock isochron defined by samples from Apollo 14, 15, 16, and 17 landing sites (4348 ± 25 Ma), and a Ferroan Anorthosite-suite Sm-Nd whole rock isochron defined by samples from the Apollo 15 and 16 landing sites (4354 ± 29 Ma). This implies that Ferroan Anorthosite-suite magmatism is temporally distinct and earlier than magmatism associated with the Mg-suite and the Felsic/Alkali-suite, as predicted by the lunar magma ocean model of lunar differentiation. The short 35 ± 10 Ma interval between primary ferroan anorthosite magmatism and secondary magmatism suggests that the lunar crust formed over a limited period of time. Although heat from decay of long-lived isotopes, large impacts, tidal heating associated with interactions between the Earth and Moon, and density driven overturn of the magma ocean have all been invoked to explain production of ancient secondary crustal magmatism, only tidal heating and cumulate overturn are consistent with the apparent short duration of secondary crustal magmatism and the great depth of crystallization implied for some Mg-suite samples. The initial ε 143 Nd values derived from the 15405 Clast B and 78238 Mg-suite norite isochrons, as well as a Mg-suite whole rock isochron are −0.23 ± 0.11, −0.27 ± 0.74, and −0.25 ± 0.09, respectively. They are identical within uncertainty indicating that Mg-suite and Felsic/Alkali-suite magmas were derived from materials that had the same time averaged Sm/Nd ratios since the formation of the solar system. This, combined with the contemporaneous nature of 15405 Clast B and 78235/6/8/55/56 Mg-suite norite, is consistent with evolution of both samples, and likely both magma suites, from a common source through closed system fractional crystallization or partial melting processes.

Chronology↗
DOE OSTI2024

Search for the Chiral Magnetic Effect with charge-dependent azimuthal correlations in Xe–Xe collisions at s NN = 5.44 TeV

Charge-dependent two- and three-particle correlations measured in Xe–Xe collisions at $\sqrt{s_{NN}}$ = TeV are presented. Results are obtained for charged particles in the pseudorapidity range |η| and transverse momentum interval 02. ≤ p T < GeV/c for different collision centralities. The three-particle correlator γ αβ ≡ $\langle$cos(φ α + φ β - 2Ψ a )$\rangle$, calculated for different combinations of charge sign α and β, is expected to be sensitive to the presence of the Chiral Magnetic Effect (CME). Its magnitude is similar to the one observed in Pb–Pb collisions in contrast to a smaller CME signal in Xe–Xe collisions than in Pb–Pb collisions predicted by Monte Carlo (MC) calculations including a magnetic field induced by the spectator protons. These observations point to a large non-CME contribution to the correlator. Furthermore, the charge dependence of γ αβ can be described by a blast wave model calculation that incorporates background effects and by the Anomalous Viscous Fluid Dynamics model with values of the CME signal consistent with zero. The Xe–Xe and Pb–Pb results are combined with the expected CME signal dependence on the system size from the MC calculations including a magnetic field to obtain the fraction of CME contribution in γ αβ , $f$ CME . The CME fraction is compatible with zero for the 30% most central events in both systems and then becomes positive. This yields an upper limit of 2% (3%) and 25% (32%) at 95% (99.7%) confidence level for the CME signal contribution to γ αβ in the 0–70% Xe–Xe and Pb–Pb collisions, respectively.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗