SEARCH · Engineering Papers

Results for “applied statistics”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,135 records · Page 63

All sources NASA NTRS DOE OSTI

DOE OSTI2023

The Atacama Cosmology Telescope: map-based noise simulations for DR6

Abstract The increasing statistical power of cosmic microwave background (CMB) datasets requires a commensurate effort in understanding their noise properties. The noise in maps from ground-based instruments is dominated by large-scale correlations, which poses a modeling challenge. This paper develops novel models of the complex noise covariance structure in the Atacama Cosmology Telescope Data Release 6 (ACT DR6) maps. We first enumerate the noise properties that arise from the combination of the atmosphere and the ACT scan strategy. We then prescribe a class of Gaussian, map-based noise models, including a new wavelet-based approach that uses directional wavelet kernels for modeling correlated instrumental noise. The models are empirical, whose only inputs are a small number of independent realizations of the same region of sky. We evaluate the performance of these models against the ACT DR6 data by drawing ensembles of noise realizations. Applying these simulations to the ACT DR6 power spectrum pipeline reveals a ∼ 20% excess in the covariance matrix diagonal when compared to an analytic expression that assumes noise properties are uniquely described by their power spectrum. Along with our public code,mnms, this work establishes a necessary element in the science pipelines of both ACT DR6 and future ground-based CMB experiments such as the Simons Observatory (SO).

Astronomy & Astrophysics↗

DOE OSTI2022

Variable wildfire impacts on the seasonal water temperatures of western US streams: A retrospective study

Recent increases in the burn area and severity of wildfires in the western US have raised concerns about the impact on stream water temperature–a key determinant of cold-water fish habitats. However, the effect on seasonal water temperatures of concern, including winter and summer, are not fully understood. In this study, we assessed the impact of wildfire burns at Boulder Creek (Oregon), Elk Creek (Oregon), and Gibbon River (Wyoming) watersheds on the downstream winter and summer water temperatures for the first three post-fire years. To obtain results independent of the choice of the analytical method, we evaluated the consequence of each burn using three different statistical approaches that utilize local water temperature data. Our results from the three approaches indicated that the response of water temperatures to wildfire burns varied across seasons and sites. Wildfire burns were associated with a median increase of up to 0.56°C (Standard Error; S.E. < 0.23°C) in the summer mean water temperatures (MWT) and 62 degree-day Celsius (DDC; S.E. < 20.7 DDC) in the summer accumulated degree days (ADD) for the three subsequent years across studied stream sites. Interestingly, these burns also corresponded to a median decrease of up to 0.49°C (S.E. < 0.45°C) in the winter MWT and 39 DDC (S.E. < 40.5 DDC) in the winter ADD for the same period across sites. Wildfire effects on the downstream water temperatures diminished with increasing site distance from the burn perimeter. Our analyses demonstrated that analytical methods that utilize local watershed data could be applied to evaluate fire effects on downstream water temperatures.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Total Measurement Uncertainty in Neutron Coincidence Multiplicity Analysis

Neutron multiplicity counting is the most commonly used nondestructive assay technique for determining the plutonium mass within containers of scrap PuO 2 or mixed oxide (MOX). In multiplicity analysis, the 240 Pu eff mass, leakage multiplication, and alpha ratio (the ratio of [α, n]-to-spontaneous fission neutron production) are the three primary unknown sample properties. They must be determined simultaneously. To solve for these three unknowns in a multiplicity assay, three measured values are needed: the singles, doubles, and triples neutron count rates. While the analysis is limited to solving for three unknowns, there are many additional factors that impact the observed count rates and contribute to the measurement uncertainty. In this study we investigate the various uncertainty contributors for the multiplicity analysis through a combination of traditional uncertainty propagation techniques supplemented by Monte Carlo simulations to address the dependences not explicitly expressed by the point source model. Uncertainties arising from counting statistics, calibration parameters, calibration method, nuclear data, and various material characteristics (isotopic abundances, chemical form, density, and impurities) are considered. A Total Measurement Uncertainty (TMU) estimate is then developed from these uncertainty contributors. This study is confined to multiplicity analysis of items commonly encountered in international safeguards applications. That is, the study focused on Pu oxides and MOX materials for the masses ranging up to 4000 grams total Pu. Multiplicity measurements were simulated using MCNP V6 based on the Plutonium Scrap Multiplicity Counter (PSMC), Epithermal Multiplicity Counter (ENMC), Pyrochemical Multiplicity Counter, and Large Epithermal Multiplicity Counter (LEMC) for this study; however, this report focuses on the parameterization of the uncertainties for the PSMC. The performance differences between the PSMC and the other multiplicity counting systems are relatively small, primarily manifesting in the impact on measurement precision so that the evaluation developed for the PSMC can be applied to the other multiplicity counting systems. Finally an analysis tool, the Multiplicity TMU Estimator, was developed from this study to serve as an aid for evaluation of the total measurement uncertainty of multiplicity assay results obtained from the commonly used INCC acquisition and analysis software.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

Three dimensional cluster analysis for atom probe tomography using Ripley’s K-function and machine learning

The size and structure of spatial molecular and atomic clustering can significantly impact material properties and is therefore important to accurately quantify. Ripley’s K-function (K(r)), a measure of spatial correlation, can be used to perform such quantification when the material system of interest can be represented as a marked point pattern. This work demonstrates how machine learning models based on K (r)-derived metrics can accurately estimate cluster size and intra-cluster density in simulated three dimensional (3D) point patterns containing spherical clusters of varying size; over 90% of model estimates for cluster size and intra-cluster density fall within 11% and 18% error of the true values, respectively. These K (r)-based size and density estimates are then applied to an experimental APT reconstruction to characterize MgZn clusters in a 7000 series aluminum alloy. Here we find that the estimates are more accurate, consistent, and robust to user interaction than estimates from the popular maximum separation algorithm. Using K (r) and machine learning to measure clustering is an accurate and repeatable way to quantify this important material attribute.

36 MATERIALS SCIENCE↗

DOE OSTI2023

A forward modeling approach to analyzing galaxy clustering with S IM BIG

We present cosmological constraints from a simulation-based inference (SBI) analysis of galaxy clustering from the SimBIG forward modeling framework. SimBIG leverages the predictive power of high-fidelity simulations and provides an inference framework that can extract cosmological information on small nonlinear scales. In this work, we apply SimBIG to the Baryon Oscillation Spectroscopic Survey (BOSS) CMASS galaxy sample and analyze the power spectrum, P ℓ (k), to k max = 0.5 h/Mpc. We construct 20,000 simulated galaxy samples using our forward model, which is based on 2,000 high-resolution Quijote N -body simulations and includes detailed survey realism for a more complete treatment of observational systematics. We then conduct SBI by training normalizing flows using the simulated samples and infer the posterior distribution of ΛCDM cosmological parameters: Ω m , Ω b , h, n s , σ 8 . We derive significant constraints on Ω m and σ 8 , which are consistent with previous works. Our constraint on σ 8 is 27% more precise than standard P ℓ analyses because we exploit additional cosmological information on nonlinear scales beyond the limit of current analytic models, k > 0.25 h/Mpc. This improvement is equivalent to the statistical gain expected from a standard P ℓ analysis of galaxy sample ~ 60% larger than CMASS. While we focus on P ℓ in this work for validation and comparison to the literature, SimBIG provides a framework for analyzing galaxy clustering using any summary statistic. We expect further improvements on cosmological constraints from subsequent SimBIG analyses of summary statistics beyond P ℓ .

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Detection of significant antiviral drug effects on COVID-19 with reasonable sample sizes in randomized controlled trials: A modeling study

Development of an effective antiviral drug for Coronavirus Disease 2019 (COVID-19) is a global health priority. Although several candidate drugs have been identified through in vitro and in vivo models, consistent and compelling evidence from clinical studies is limited. The lack of evidence from clinical trials may stem in part from the imperfect design of the trials. We investigated how clinical trials for antivirals need to be designed, especially focusing on the sample size in randomized controlled trials. A modeling study was conducted to help understand the reasons behind inconsistent clinical trial findings and to design better clinical trials. We first analyzed longitudinal viral load data for Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) without antiviral treatment by use of a within-host virus dynamics model. The fitted viral load was categorized into 3 different groups by a clustering approach. Comparison of the estimated parameters showed that the 3 distinct groups were characterized by different virus decay rates (p-value < 0.001). The mean decay rates were 1.17 d -1 (95% CI: 1.06 to 1.27 d -1 ), 0.777 d -1 (0.716 to 0.838 d -1 ), and 0.450 d -1 (0.378 to 0.522 d -1 ) for the 3 groups, respectively. Such heterogeneity in virus dynamics could be a confounding variable if it is associated with treatment allocation in compassionate use programs (i.e., observational studies). Subsequently, we mimicked randomized controlled trials of antivirals by simulation. An antiviral effect causing a 95% to 99% reduction in viral replication was added to the model. To be realistic, we assumed that randomization and treatment are initiated with some time lag after symptom onset. Using the duration of virus shedding as an outcome, the sample size to detect a statistically significant mean difference between the treatment and placebo groups (1:1 allocation) was 13,603 and 11,670 (when the antiviral effect was 95% and 99%, respectively) per group if all patients are enrolled regardless of timing of randomization. The sample size was reduced to 584 and 458 (when the antiviral effect was 95% and 99%, respectively) if only patients who are treated within 1 day of symptom onset are enrolled. We confirmed the sample size was similarly reduced when using cumulative viral load in log scale as an outcome. We used a conventional virus dynamics model, which may not fully reflect the detailed mechanisms of viral dynamics of SARS-CoV-2. The model needs to be calibrated in terms of both parameter settings and model structure, which would yield more reliable sample size calculation. In this study, we found that estimated association in observational studies can be biased due to large heterogeneity in viral dynamics among infected individuals, and statistically significant effect in randomized controlled trials may be difficult to be detected due to small sample size. The sample size can be dramatically reduced by recruiting patients immediately after developing symptoms. We believe this is the first study investigated the study design of clinical trials for antiviral treatment using the viral dynamics model.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

Reassessing the Origins and Contemporary Relevance of ck Acceptability Parameters: Evolving Perspectives on Similarity

“Sensitivity and Uncertainty Analyses Applied to Criticality Safety Validation,” introduces sensitivity and uncertainty methods to address challenges in defining and extending areas of applicability for criticality safety validation. These areas are traditionally defined by the bounds or limits on key parameters, but establishing valid ranges and managing complex parameter variations remain challenging. NUREG/CR-6655 introduces ck and other integral indices, as well as concepts such as the completeness of benchmark coverage, to better quantify system similarities. The work proposed herein seeks to evaluate these foundational concepts to ensure that the bounds remain effective in guiding the assessment of similarity and applicability in modern applications. The concept of completeness, along with other parameters envisioned within the framework, serves as an example of the foundational ideas that have been established, though their effectiveness in practice may not be fully understood. Advancements in scripting tools, coupled with the speed and efficiency of modern computing and statistical models, now allow for faster and more thorough assessments than previously possible. These advancements also enable the identification of trends within the data, which could provide additional insight into system behavior and further broaden the scope of previously performed benchmarks. By leveraging these capabilities, we will revisit and expand the scope of these foundational methods to determine whether the necessary elements for robust similarity evaluation are already embedded, partially realized, or remain untapped.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Advanced EXAFS analysis techniques applied to the L -edges of the lanthanide oxides

The unique properties of the lanthanide (Ln) elements make them critical components of modern technologies, such as lasers, anti-corrosive films and catalysts. Thus, there is significant interest in establishing structure–property relationships for Ln-containing materials to advance these technologies. Extended X-ray absorption fine structure (EXAFS) is an excellent technique for this task considering its ability to determine the average local structure around the Ln atoms for both crystalline and amorphous materials. However, the limited availability of EXAFS reference spectra of the Ln oxides and challenges in the EXAFS analysis have hindered the application of this technique to these elements. The challenges include the limited k-range available for the analysis due to the superposition of L-edges on the EXAFS, multielectron excitations (MEEs) creating erroneous peaks in the EXAFS and the presence of inequivalent absorption sites. Herein, we removed MEEs to model the local atomic environment more accurately for light Ln oxides. Further, we investigated the use of cubic and non-cubic lattice expansion to minimize the fitting parameters needed and connect the fitting parameters to physically meaningful crystal parameters. The cubic expansion reduced the number of fitting parameters but resulted in a statistically worse fit. The non-cubic expansion resulted in a similar quality fit and showed non-isotropic expansion in the crystal lattice of Nd 2 O 3 . In total, the EXAFS spectra and the fits for the entire set of Ln oxides (excluding promethium) are included. The knowledge developed here can assist in the structural determination of a wide variety of Ln compounds and can further studies on their structure–property relationships.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Detecting Thermally Induced Spinodal Decomposition with Picosecond Ultrasonics in Cast Austenitic Stainless Steels

Given the existential climate crisis faced by mankind and the world, the lifetime and sustainability of nuclear reactors as a carbon-free source of renewable energy depend on the susceptibility of their structural components to environmental degradation. In particular, critical components for light water reactors (LWRs) evolve over decades in service, losing ductility and toughness due to thermal and irradiation aging. Techniques to monitor their health cannot be easily applied in the field due to their destructive, expensive, or immobile nature. Thus, non-destructive evaluation (NDE) methods are sought to monitor and evaluate the health of major LWR components such as core barrels, steam generator tubes, or primary coolant pipes and are often required by policy, such as NRC policy #10-CFR-50.65. In this work we demonstrate the use of gigahertz, non-contact ultrasonics to gauge the state of cast austenitic stainless steels (CASS), used in some of the largest components in LWR primary systems. We do so by linking changes in their surface acoustic wave (SAW) characteristics using transient grating spectroscopy (TGS) to transmission electron microscopy (TEM)-verified evidence of spinodal decomposition and G-phase precipitation. In this thesis, thermal aging is shown to induce SAW peak splitting in spinodally decomposed CASS alloys, correlated strongly with lowered toughness and decreased ductility. Furthermore, statistical testing on the number of SAW peak splits observed show that the second SAW peak significantly appears more frequently and is significantly different in frequency in comparison to counts and frequencies measured in unaged specimens. The ability of this technique to non-destructively detect microstructural degradation at a distance in a predictive manner in the case of CASS motivates extending gigahertz ultrasonics to detect other LWR material degradation modes as an in-vessel inspection technique, such as reactor pressure vessel (RPV) embrittlement. This allows for the greater use of NDE techniques for confident monitoring of LWR structural material health to 80 years and beyond, saving costs by minimizing structural replacements until needed and maximizing energy production by preventing early decommission until necessary.

36 MATERIALS SCIENCE↗

DOE OSTI2025

An analytically tractable marked power spectrum

The increasing precision of cosmology data in the modern era is calling for methods to allow the extraction of non-Gaussian information using tools beyond two-point statistics. The marked power spectrum has the potential to extract beyond two-point information in a computationally efficient way while using much of the infrastructure already available for the power spectrum. In this work we explore the marked power spectrum from an analytical perspective. In particular, we explore a low-order polynomial for the mark that allows us to better control the theoretical uncertainties and we show that with minimal new degrees of freedom the analytical results match measurements from N-body simulations for both the matter field and biased tracers in redshift space. Finally, we show that even within the limited forms of mark that we consider, there are degeneracies that can be broken by inclusion of the marked auto-spectrum or the cross-spectrum with the unmarked field. I n conclusion, we discuss future theoretical developments that would enable us to apply this approach to survey data.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

An extended radar relative calibration adjustment (eRCA) technique for higher-frequency radars and range–height indicator (RHI) scans

Abstract. This study extends the relative calibration adjustment technique for calibration of weather radars to higher-frequency radars as well as range–height indicator (RHI) scans. The calibration of weather radars represents one of the most dominant sources of error for their use in a variety of fields including quantitative precipitation estimation and model comparisons. While most weather radars are routinely calibrated, the frequency of calibration is often less than required, resulting in miscalibrated time periods. While full absolute calibration techniques often require the radar to be taken offline for a period of time, there have been online calibration techniques discussed in the literature. The relative calibration adjustment (RCA) technique uses the statistics of the ground clutter surrounding the radar as a monitoring source for the stability of calibration but has only been demonstrated to work at S- and C-band for plan-position indicator (PPI) scans at a constant elevation. In this work the RCA technique is modified to work with higher-frequency radars, including Ka-band cloud radars. At higher frequencies the properties of clutter can be much more variable. This work introduces an extended clutter selection procedure that incorporates the temporal stability of clutter and helps to improve the operational stability of RCA for relatively higher-frequency radars. The technique is also extended to utilize range–height scans from radars where the elevation is varied rather than the azimuth. These types of scans are often utilized with research radars to examine the vertical structure of clouds. The newly extended technique (eRCA) is applied to four Department of Energy Atmospheric Radiation Measurement (DOE ARM) weather radars ranging in frequency from C- to Ka-band. Cross comparisons of three co-located radars with frequencies C, X, and Ka at the ARM Cloud, Aerosol, and Complex Terrain Interactions (CACTI) site show that the technique can determine changes in calibration. Using an X-band radar at the ARM Eastern North Atlantic (ENA) site, we show how the technique can be modified to be more resilient to clutter fields that show increased variability, in this case due to sea clutter. The results show that this technique is promising for a posteriori data calibration and monitoring.

47 OTHER INSTRUMENTATION↗

DOE OSTI2023

Relationship between Grain Size and Sample Thickness on the Creep-Rupture Performance of Thin Metallic Sheets of INCONEL Alloy 740H

A study was conducted on INCONEL® alloy 740H® to examine the role of grain size and sheet thickness on the alloy’s creep-rupture behavior. Three different starting sheet thicknesses were utilized and multiple heat-treatment conditions anticipated for compact heat exchanger (CHX) manufacturing were applied to produce a range of grain sizes. Creep-rupture testing was conducted at 750 °C for times up to about 6000 h and the results were compared to wrought databases. The data show that both creep strength and ductility were important factors in the overall creep performance of the sheets. Reductions in performance were observed due either to accelerated creep when grain size was fine or loss of rupture ductility when grain sizes approached the sheet thickness. Some combinations of heat-treatment and thickness were able to produce typical expected wrought creep properties. Historically a ‘rule-of-thumb’ requirement for creep testing suggests 3-5 grains per sample minimum dimension to ensure homogeneous behavior. This research shows that to ensure representative wrought creep performance (i.e. no effect of sample size), the sample minimum cross-section should be 10 times the average ASTM grain size. Statistical analysis of the microstructures suggests the population of larger grains as a controlling feature in creep failure.

14 SOLAR ENERGY↗

DOE OSTI2024

Precision Computations in Strongly Coupled Conformal Field Theories (Final Technical Report)

Conformal Field Theories (CFTs) are quantum field theories that are invariant under the conformal symmetry group (which includes translations and rotations, but also local rescalings of spacetime). They are building blocks of general quantum field theories, and appear in many areas of physics, including statistical physics, condensed matter physics, particle physics, and quantum gravity. Because of their extra symmetries, the mathematical structure of CFTs is tightly constrained, and this leads to the idea of the ``conformal bootstrap," which is to use these mathematical structures to constrain, and in some cases determine, CFT observables. A new numerical implementation of the conformal bootstrap idea appeared in 2008 with the work of Rattazzi, Rychkov, Tonni, and Vichi. Their observation was that certain bootstrap constraints (conformal symmetry and unitarity) could be combined to yield a convex optimization problem that constraints CFT data. By solving this convex optimization problem on a computer, one could obtain bounds on observables like critical exponents and operator product expansion (OPE) coefficients. Over the course of this award, the PI has improved numerical bootstrap techniques by optimizing known algorithms and finding new ones for performing the required convex optimization computations. The PI has applied these techniques to compute high-precision observables in several important strongly-coupled systems. The PI has also explored both analytical and numerical bootstrap methods for constraining the space of low energy effective field theories of quantum gravity, and developed new analytical techniques for CFT and QFT more broadly.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

An automated integrated web-based smart tool for open stope design

The Stability Graph is a widely used tool for the design of open stopes in underground mining. Many users of the Stability Graph still apply this design method manually. Although the manual approach has benefits, using multiple graphs and stability number computation charts for each stope surface is time-consuming, even for the experienced mining engineer. Current practice in the use of the method also limits data sharing. This paper presents a StopeSoft web-based tool for open stope stability prediction that is developed on the basis of the Stability Graph method and is available at openstope.com. StopeSoft incorporates flexibility in terms of Stability Graph options and incorporates additional critical factors often overlooked. As a web-based tool, StopeSoft encourages and makes data sharing possible globally, focused on expanding the database and improving the current limitations of the Stability Graph to provide practical, reliable solutions for mining engineers, consultants, and academics. The StopeSoft automated process facilitates the process of open stope stability prediction, saving time and minimizing potential human errors. Statistical treatment of the data accounts for the variability of input parameters to emphasize the probabilistic nature of the Stability Graph method. The probabilistic interpretation of the stability states of stope surfaces eliminates the false feeling of absolute stope performance based on its location on the Stability Graph , as implied by the deterministic approach.

58 GEOSCIENCES↗

DOE OSTI2020

Simulation scaling studies of reactor core two-phase flow using direct numerical simulation

Tremendous growth in supercomputing power in recent years has resulted in the emergence of high-resolution flow analysis methods as an advanced research tool to evaluate single and two-phase flow behavior. In particular, unstructured mesh-based methods have been applied to analyze flows in complex reactor core geometries, including those of light water reactors (LWR). The finite-element based code, PHASTA, is utilized to perform large-scale simulations of two-phase bubbly flows in LWR geometries. Given the large computational cost of direct numerical simulation (DNS) coupled with interface tracking methods (ITM), typical domains encompass a portion of a single subchannel. In the presented research, the state-of-the-art analysis of turbulent two-phase flows in complex LWR subchannel geometries are demonstrated at both prototypical reactor parameters as well as scaled low pressure conditions. Three different cases are studied, a high-pressure simulation in prototypical reactor subchannel geometry, a low-pressure case in prototypical geometry and a final low-pressure case in a geometry scaled up to conserve the ratio between the bubble size and the domain pitch. Utilizing advanced statistical processing tools, these simulation conditions are compared to shed light on the relevancy of two-phase flow characteristics given the significant differences between LWR and low-pressure conditions. These findings can lead to the generation of useful guiding principles when researchers need to scale the two-phase flow behavior captured at low pressure and temperature conditions to those at reactor operating conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2026

Machine learning method for enforcing variable independence in background estimation with LHC data: ABCDisCoTEC

A novel solution is presented for the problem of estimating the backgrounds of a signal search using observed data while simultaneously maximizing the sensitivity of the search to the signal. The 'ABCD method' provides a reliable framework for background estimation by partitioning events into one signal-enhanced region (A) and three background-enhanced control regions (B, C, and D) via two smoothly varying, statistically independent variables. In practice, even slight correlations between the two variables can significantly undermine the method's performance. Thus, choosing appropriate variables by hand can present a formidable challenge, especially when background and signal differ only subtly. To address this issue, the ABCD with distance correlation (ABCDisCo) method was developed to construct two learned variables via a neural network trained to provide strong signal-background discrimination with small values of the distance correlation (DisCo) measure between the two learned variables. However, relying solely on minimizing the DisCo can result in learned variables that may not have distributions of background events that are smoothly varying and localized at extreme values, as necessary for the validity of the background estimation. The ABCDisCo training enhanced with closure (ABCDisCoTEC) method is introduced to solve this issue by directly minimizing the nonclosure, expressed as a dedicated differentiable loss term. This extended method is applied to a data set of proton-proton collisions at a center-of-mass energy of 13 TeV recorded by the CMS detector at the CERN Large Hadron Collider. Additionally, given the complexity of the minimization problem with constraints on multiple loss terms, the modified differential method of multipliers is applied and shown to greatly improve the stability and robustness of the ABCDisCoTEC method, compared to grid search hyperparameter optimization procedures.

Hayrapetyan, Aram [Yerevan Phys. Inst.]↗

DOE OSTI2024

Visibility and annoyance of the phantom array effect varies with age and history of migraine

The phantom array effect (PAE) is a series of repeated images that may be perceived when a person moves their eyes in large saccades across a light source (or a specular reflection of that light source) that is modulating in output over time. Fifty-five people, including a group of 25 who experience migraine, evaluated the visibility and annoyingness of phantom arrays produced by 85 unique temporal light modulation waveforms (including sine, rectangular, complex and DC waveforms) generated using an LED placed against a black background. Those with migraine exhibited higher average visibility compared to those without migraine ( p = 0.019) and were relatively more sensitive at higher frequencies ( p < 0.001). Younger participants also found more stimuli to be visible ( p < 0.001). The threshold sensitivity function was similar to that developed for the phantom array visibility measure (PAVM), and PAVM was effective in predicting visibility ( R 2 = 0.87 for the relevant region of PAVM < 3). While those in the migraine group did not report seeing the PAE more often in everyday life at a statistically significant level, they reported being more annoyed by it and having more unwanted physiological responses (headaches, eye fatigue and distraction/disorientation). Members of the migraine group were also more likely to have changed their behaviour in architectural spaces (such as leaving a restaurant with ‘flickering’ lights). In the four hours after completing the experiment, 64% of the migraine group (vs. 19% of the non-migraine group) reported experiencing discomfort or an adverse reaction. In particular, 41% reported experiencing a headache (vs. 8% for those in the non-migraine group).

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

Searches for New Physics With Muon Conversion at Fermilab and Triboson Production at the LHC

We report on several efforts to search for physics beyond the standard model of particle physics at broad energy scales. The Mu2e experiment at Fermilab will search for charged lepton flavor violation via the muon to electron conversion process, which is suppressed in the Standard Model. Mu2e will be operated at a low energy, yet can probe New Physics at very high mass scales (O(1e3 - 1e4 ) TeV). At high energies, the CMS experiment at the CERN LHC continues to deliver an impressive suite of Standard Model measurements and limits on a variety of New Physics signatures. Mu2e is under construction and slated to collect its first physics data in the coming years. This thesis describes work done during the construction phase of Mu2e and focuses on two critical areas: magnetic field modeling and statistical analysis. We describe a novel method for field modeling which we validate using a simulated dataset representing the expected magnetic field in the Detector Solenoid. This method blends a standard least-squares fitting technique that utilizes physically motivated analytical model functions with a novel physics informed network that is constructed to obey Maxwell’s equations. We show the technique can model the field with an accuracy of 10−7 despite the presence of injected noise in the pseudo-measurements at the 10−5 level. We then present preliminary results of the calibration of 3D Hall probes at the sub-10−4 level. These probes will be used to directly measure the Mu2e Detector Solenoid magnetic field on a sparse grid; these measurements serve as the input to the field model fitting. Finally, we describe the first implementation of both an unbinned shape analysis and a Bayesian interpretation applied to Mu2e pseudo-data. Up to 20% tighter limits can be set by the shape analysis compared to a standard cut & count analysis. The AlCap experiment collected data at PSI in 2015 to measure several important quantities related to nuclear muon capture on an aluminum target, which is a significant background process for Mu2e. The neutron emission from muon capture can introduce background hits in the Mu2e detectors and can increase radiation damage in various elements of the apparatus. We present measurements of the neutron group fluence and mean neutron multiplicity for muon capture on aluminum nuclei. Finally, we discuss an analysis of triboson production at CMS using an Effective Field Theory framework. Standard Model triboson production, which was first observed at CMS in 2020, has a relatively small cross section and provides direct access to both anomalous triple gauge couplings and quartic gauge couplings. These couplings, interpreted in the Standard Model Effective Field Theory, are studied in the present work. We target the boosted regime where the background rate is low and yields are enhanced when dimension-6 and dimension-8 Wilson coefficients are non-zero. We do not observe an excess in the data and therefore set bounds on the Wilson coefficients. For dimension-6 coefficients the tightest observed (expected) bounds are set on cW /Λ2 where Λ is the mass scale of new physics; the bounds are [−0.13, 0.12] TeV−2 ([−0.12, 0.12] TeV−2 ) at 95% CL. The tightest bounds in dimension-8 are set on fT,0 / Λ4 ; the observed (expected) bounds at 95% CL are [−0.63, 0.69] TeV−4 ([−0.54, 0.62] TeV−4 ). Additional results are presented which include scenarios where multiple Wilson coefficients are non-zero, the application of signal model clipping to address unitarity violation in Effective Field Theories, and a novel template fit developed for easier reinterpretation of our results.

Kampa, Cole Erik [Northwestern U. (main)] (ORCID:0↗