Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Particle data analysis”

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 973 records · Page 54

Machine Learning-Based Energy Estimation for Sterile Neutrino Searches in the NOνA Experiment

This dissertation presents a search for sterile neutrinos using Monte Carlo datasets and experimental data from the NOνA experiment. This work introduces novel machine learning techniques for energy reconstruction in neutral current events. A deep learning-based energy estimator was developed and integrated into the analysis framework. Using this new energy estimator, the sensitivity to sterile neutrino-induced oscillations is evaluated and presented. Furthermore, the dissertation explores potential methods for further improvement of the analysis.

43 PARTICLE ACCELERATORS↗

Sub-pixel high-resolution imaging of high-energy x-rays inspired by sub-wavelength optical imaging

We have developed and demonstrated an image super-resolution method—XR-UNLOC: X-Ray UNsupervised particle LOCalization—for hard x-rays measured with fast-frame-rate detectors that is an adaptation of the principle of photo-activated localization microscopy (PALM) and stochastic optical reconstruction microscopy (STORM), which enabled biological fluorescence imaging at sub-optical-wavelength scales. We demonstrate the approach on experimental coherent Bragg diffraction data measured with 52 keV x-rays from a nanocrystalline sample. From this sample, we resolve the fine fringe detail of a high-energy x-ray Bragg coherent diffraction pattern to an upsampling factor of 16 of the native pixel pitch of 30 μ m of a charge-integrating fastCCD detector. This was accomplished by analysis of individual photon locations in a series of “nearly-dark” instances of the diffraction pattern that each contain only a handful of photons. Central to our approach was the adaptation of the UNLOC photon fitting routine for PALM/STORM to the hard x-ray regime to handle much smaller point spread functions, which required a different statistical test for photon detection and for sub-pixel localization. A comparison to a photon-localization strategy used in the x-ray community (“droplet analysis”) showed that XR-UNLOC provides significant improvement in super-resolution. We also developed a metric by which to estimate the limit of reliable upsampling with XR-UNLOC under a given set of experimental conditions in terms of the signal-to-noise ratio of a photon detection event and the size of the point spread function for guiding future x-ray experiments in many disciplines where detector pixelation limits must be overcome.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Approximate entropy analysis for nonlinear beam dynamics

In this paper, we applied an approximate entropy (ApEn) analysis to nonlinear beam dynamics. Due to the presence of strong nonlinear magnets, the degree of chaos in beam motion gradually increases proportionally with its amplitude. Such chaos can be quantitatively characterized with ApEn of beam turn-by-turn readings. ApEn analysis is a technique used to quantify the amount of regularity and the unpredictability of fluctuations with respect to time-series data. The ApEn, when applied as a chaos indicator, can then be used for nonlinear lattice optimization and analysis. The National Synchrotron Light Source II (NSLS-II) electron storage ring lattice was used as an example of a real-world application of our technique. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Theory and applications of parton pseudodistributions

We review the basic theory of the parton pseudodistributions approach and its applications to lattice extractions of parton distribution functions. The crucial idea of the approach is the realization that the correlator $M(z,p)$ of the parton fields is a function $M(ν,-z^2)$ of Lorentz invariants $ν=-(zp)$, the Ioffe time, and the invariant interval $z^2$. This observation allows to extract the Ioffe-time distribution $M(ν,-z^2)$ from Euclidean separations z accessible on the lattice. Another basic feature is the use of the ratio $M(ν,-z^2)≡M(ν,-z^2)/M(0,-z^2)$, that allows to eliminate artificial ultraviolet divergence generated by the gauge link for spacelike intervals. The remaining $z^2$-dependence of the reduced Ioffe-time distribution $M(ν,-z^2)$ corresponds to perturbative evolution, and can be converted into the scale-dependence of parton distributions $f(x,μ^2)$ using matching relations. The ν-dependence of $M(ν,-z^2)$ governs the x-dependence of parton densities $f(x,μ^2)$. The perturbative evolution was successfully observed in exploratory quenched lattice calculation. The analysis of its precise data provides a framework for extraction of parton densities using the pseudodistributions approach. It was used in the recently performed calculations of the nucleon and pion valence quark distributions. We also discuss matching conditions for the pion distribution amplitude and generalized parton distributions, the lattice studies of which are now in progress.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Towards a data-driven model of hadronization using normalizing flows

We introduce a model of hadronization based on invertible neural networks that faithfully reproduces a simplified version of the Lund string model for meson hadronization. Additionally, we introduce a new training method for normalizing flows, termed MAGIC, that improves the agreement between simulated and experimental distributions of high-level (macroscopic) observables by adjusting single-emission (microscopic) dynamics. Our results constitute an important step toward realizing a machine-learning based model of hadronization that utilizes experimental data during training. Finally, we demonstrate how a Bayesian extension to this normalizing-flow architecture can be used to provide analysis of statistical and modeling uncertainties on the generated observable distributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Testing the Jump Finding Code for NEXUS Qubit Analysis

At Fermilab, the NEXUS dilution refrigerator is repeating the same experiment from Wilen et al. 1 underground to further investigate the effect of ionizing radiation and cosmic ray muons on qubit decoherence. This refrigerator is 100m underground in the MINOS cavern to significantly lower the incident rate of cosmic ray muons, meaning the majority of incidents is due to ionizing radiation. To gain a better understanding of how the rate of ionizing radiation affects the qubits, we use 4 radiation configurations to collect data from the chip. In this work, we examine the efficiency of our jump detection code to find jumps of various sizes. We do so using simulated data with a known number of injected jumps of known sizes. These jumps are charge offsets of the qubit, thought to be due to quasiparticle poisoning. To test this, we run the jump finding code on an array of simulated data in which two main parameters, the smoothing factor and the threshold at which a charge offset is considered a jump, are changed to test an array of these parameters. We then perform a χ2 analysis to find the best combination of these two parameters.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Aerosol impacts on isolated deep convection: findings from TRACER

This study focuses on quantifying the conditional relationship between aerosol and convective precipitation properties of isolated deep convective clouds (DCCs) in the Houston–Galveston region, after adjusting for confounding effects. We leverage comprehensive ground-based observations from the TRacking Aerosol Convection interactions ExpeRiment (TRACER) to estimate aerosol effects on convective echo top height (ETH), intensity, and area separately. Our results show that greater aerosol number concentrations generally have a limited impact on these convective properties, showing relationships consistent with the possibility of both invigoration and suppression effects. Under certain conditions, where ultrafine particles are abundant, aerosols exhibit a positive effect on ETH, increasing it by about 1 km. However, it is inevitable to consider measurement uncertainties and the limitations of temporal and spatial resolution in the data, as these factors can further contribute to uncertainties in our estimates. In DCCs associated with sea breezes, the estimated aerosol effects on DCCs are found to be more pronounced. However, this heightened effect could be attributed to the exclusion of key confounders such as boundary layer updrafts in the analysis.

54 ENVIRONMENTAL SCIENCES↗

Measurement of the low-energy antitriton inelastic cross section

In this Letter, the first measurement of the inelastic cross section for antitriton–nucleus interactions is reported, covering the momentum range of 0.8 ≤ p < 2.4 GeV/c. The measurement is carried out using data recorded with the ALICE detector in pp and Pb–Pb collisions at a centre-of-mass energy per nucleon of 13 TeV and 5.02 TeV, respectively. The detector material serves as an absorber for antitriton nuclei. The raw yield of (anti)triton nuclei measured with the ALICE apparatus is compared to the results from detailed ALICE simulations based on the GEANT4 toolkit for the propagation of (anti)particles through matter, allowing one to quantify the inelastic interaction probability in the detector material. This analysis complements the measurement of the inelastic cross section of antinuclei up to A = 3 carried out by the ALICE Collaboration, and demonstrates the feasibility of the study of the isospin dependence of inelastic interaction cross section with the analysis techniques presented in this Letter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Uncertainty Quantification of Monolithic Tungsten Target Block Fatigue

The current STS spallation target design consists of 21 water cooled, tantalum clad, monolithic tungsten blocks. The target wheel assembly rotates such that each target block receives a direct proton pulse every 1.4 seconds (21 blocks / 15 Hz pulses). The quasi-steady heating induced by the repetitive beam pulses creates a significant mean stress distribution in the target. On top of this, the rapid energy deposition from a proton pulse induces dynamic stress waves that are of the same order of magnitude as the mean stresses. Over the design goal of 5,000 operational hours per year for a decade, each target must withstand over 130 million repetitions of these dynamic stress events. The fatigue life prediction of the target blocks is a complex sequence of simulations that depend upon a set of input parameters. For the tungsten block, the end goal is the fatigue factor-of-safety (FOS). This study documents a new structural analysis routine using software available to U.S. national laboratories (e.g. CUBIT, Sierra Multiphysics, Dakota) that enables more rapid target fatigue analysis. The increased throughput is amenable to parameter sensitivity and uncertainty quantification (UQ) studies. A set of input parameters (ultimate strength, fatigue strength, density, elastic modulus, Poisson’s ratio, thermal conductivity, coefficient of thermal expansion, specific heat, flow rate) with expected ranges are documented based on literature data, recent experimental testing, and engineering judgement. UQ bounds the resulting fatigue FOS at [0.5, 3.4] using epistemic interval analysis. These bounds are expected to be wider due to model effects (discretization) and beam effects (focus, offset) that were not included. Parameter sensitivity underscores the importance of ultimate strength and fatigue strength of irradiated tungsten on the fatigue calculations.

43 PARTICLE ACCELERATORS↗

Even Higher-Level Synthesis: An Exploration of AI Hardware Accelerators using HLS4ML

With the rise of artificial intelligence, the popularization of deep learning, and a constantly evolving industry, the demand for flexible and efficient tools has never been greater. As algorithms grow more complex, their runtime and energy consumption increase exponentially. Customized hardware accelerators, long used for specific mathematical operations, remain essential for managing modern applications' computational and power demands. Hardware accelerators can speed up complex computations by orders of magnitude, but their manual design and verification processes are often challenging and time-consuming. High-Level Synthesis (HLS) provides a solution by transforming high-level algorithm descriptions, typically written in C++ or SystemC, into synthesizable RTL suitable for hardware implementation. This approach reduces development time for RTL engineers while offering flexibility beyond what traditional handwritten RTL can provide. We extended this capability to the machine-learning domain with the open-source framework hls4ml, which allows neural networks trained in Python frameworks like Tensorflow or PyTorch to be synthesized into efficient hardware representations for the traditional FPGA and ASIC flows. This breakthrough addresses the growing need for reduced design turnaround and easy verification of ML hardware accelerators with low latency and power efficiency constraints. During this tutorial, we will demonstrate how Python complements HLS by simplifying the ML design process, bridging the gap between software and hardware development. Attendees will explore how we translate neural networks modeled in Python into fixed-point C++ models suitable for HLS workflows. We will dive into strategies like Value-Range Analysis and Quantization-Aware Training, which optimize these designs for deployment and evaluate their accuracy, power consumption, and energy efficiency. To exemplify these concepts, experts from Fermilab will share their experiences applying this technology to high-energy physics experiments, where real-time, low-latency processing is critical. Over the years, Fermilab engineers have demonstrated how deep neural networks, optimized for hardware using hls4ml, can meet the stringent requirements of trigger systems at the CERN Large Hadron Collider. These systems rely on rapid decision-making to process immense data volumes while retaining only the most relevant events for further analysis. The application of hls4ml has also been extended to innovative technologies like smart pixel arrays. These smart pixels integrate ML inference capabilities directly into sensor devices, enabling localized data processing at the pixel level. This approach drastically reduces the need to transmit raw data to external processing units, significantly decreasing power consumption and latency. By embedding neural networks within the pixel architecture, the smart pixels can identify and prioritize relevant data in real time, providing a highly efficient solution for edge computing in scenarios such as particle detectors and imaging systems. Fermilab's work highlights the potential of hardware-accelerated ML in scenarios where both speed and power efficiency are mission-critical. Through this tutorial, attendees will gain valuable insights into the challenges and solutions of deploying ML in hardware. Understanding how HLS and hls4ml streamline the development of neural network-based hardware accelerators is fundamental for the industry's future. Participants will learn how these technologies are shaping the future of AI and scientific computing.

Di Guglielmo, Giuseppe [Fermilab]↗

Cross Section Evaluation for Exclusive Channels of K+Λ and K+Σ0 Electroproduction off Protons Using CLAS Detector Data

In this work, a method for evaluating the cross sections of electroproduction of K+Λ0 and K+Σ0 off protons in the region of invariant masses of final hadrons MK + MY <W < 2.65 GeV (MK and MY being the masses of the kaon and hyperon, respectively) and squares of four-momentum transfers of virtual photons, i.e., photon virtualities 0 < Q2 < 5 GeV2, is developed based on experimental data of these exclusive channels’ cross sections measured by the CLAS detector in Hall B at Jefferson Lab. A set of algorithms has been implemented to evaluate the differential cross sections of these channels, along with their statistical and systematic uncertainties. A program was developed for the evaluation of differential cross sections and structure functions using C++ and Python libraries. An interactive website was created for working with the program, enabling the analysis of one-dimensional and two-dimensional dependences of structure functions and differential cross sections. The evaluation of differential cross sections for the K+Λ and K+Σ0 electroproduction channels is necessary for extracting the structure function σLT from the data on the polarization asymmetry of electroproduction reactions of these final states with longitudinally polarized electrons. The obtained results are also important for the development of realistic Monte Carlo event generators in planning future experiments and for evaluating the efficiency of detecting final particles when extracting reaction cross sections from experimental data.

Golda, A. V.↗

Observation and Modeling of the Thermal Runaway of High-Capacity Pouch Cells Due to an Internal Short Induced by an Indenter

Owing to the complexity of coupling mechanical and electrical solvers for finite element analysis (FEA) of mechanical abuse of a battery, the simulation should be performed by incorporating a local short circuit model activated based on the status of battery elements. However, the approach should also consider various other conditions of the battery simulated, such as the initial state of charge (SOC), for precise simulations of the internal short. This study proposes an approach developed with a mechanical model and single-particle battery model, along with a reaction kinetics model for thermal abuse using LS-DYNA. Several key parameters were selected, such as the activation criteria, kinetic abuse and element integration points during the internal short circuit induced by the indenter. The impact of these key parameters in the FEA on thermal runaway response for various initial SOCs from 25 to 100 % were analyzed. The FEA results obtained were compared against experimental results with various parameters such as activation criterion of internal short, kinetic abuse based on accelerating rate calorimetry (ARC) data. Finally, the root causes of the discrepancy between the FEA and experimental results are discussed regarding integration point what solver can make difference during calculation.

ENERGY STORAGE↗

Halo effective field theory analysis of one-neutron knockout reactions of Be 11 and C 15

Background: One-nucleon knockout reactions provide insightful information on the single-particle structure of nuclei. When applied to one-neutron halo nuclei, they are purely peripheral, suggesting that they could be properly modeled by describing the projectile within a halo effective field theory (halo-EFT). Purpose: We reanalyze the one-neutron knockout measurements of 11 Be and 15 C —both one-neutron halo nuclei—on beryllium at about 60 MeV/nucleon. We consider halo-EFT descriptions of these nuclei which already provide excellent agreement with breakup and transfer data. Method: Here we include a halo-EFT description of the projectile within an eikonal-based model of the reaction and compare its outcome to existing data. Results: Excellent agreement with experiment is found for both nuclei. The asymptotic normalization coefficients inferred from this comparison confirm predictions from ab initio nuclear-structure calculations and values deduced from transfer data. Conclusions: Halo-EFT can be reliably used to analyze one-neutron knockout reactions measured for halo nuclei and test predictions from state-of-the-art nuclear structure models on these experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Development of techniques for gantryless associated-particle imaging

To move fast-neutron radiography using the associated-particle imaging technique from the laboratory to the field, the development of new analysis techniques is required. In particular, the relative positions of the source and detectors need to be determined when they have been placed by hand, the normalization for a particular source–detector geometry needs to be determined without a measurement in the same geometry with no object present, and accurate image stitching is required when multiple detector positions are necessary to image an object. The present work describes methods that employ transmission neutron data to localize a fast-neutron imaging panel with respect to the neutron source, calculate a normalization for a given source–detector geometry, stitch images together, and describes the required system calibrations. The reported techniques enable in-field neutron radiography for cases in which the source–detector geometry is not well known a priori and where operational constraints preclude a normalization measurement.

Heath, Matthew↗

Girth and groomed radius of jets recoiling against isolated photons in lead-lead and proton-proton collisions at $\sqrt{s_{NN}} = 5.02$ TeV

This Letter presents the first measurements of the groomed jet radius R g and the jet girth g in events with an isolated photon recoiling against a jet in lead-lead (PbPb) and proton-proton (pp) collisions at the LHC at a nucleon-nucleon center-of-mass energy of 5.02 TeV. The observables R g and g provide a quantitative measure of how narrow or broad a jet is. The analysis uses PbPb and pp data samples with integrated luminosities of 1.7 nb -1 and 301 pb -1 , respectively, collected with the CMS experiment in 2018 and 2017. Events are required to have a photon with transverse momentum $p^{γ}_{T} > 100$ GeV and at least one jet back-to-back in azimuth with respect to the photon and with transverse momentum $p^{jet}_{T}$ such that $p^{jet}_{T}/p^{γ}_{T} > 0.4$ . The measured R g and g distributions are unfolded to the particle level, which facilitates the comparison between the PbPb and pp results and with theoretical predictions. It is found that jets with $p^{jet}_{T}/p^{γ}_{T} > 0.8$ , i.e., those that closely balance the photon $p^{γ}_{T}$ , are narrower in PbPb than in pp collisions. Relaxing the selection to include jets with $p^{jet}_{T}/p^{γ}_{T} > 0.4$ reduces the narrowing of the angular structure of jets in PbPb relative to the pp reference. This shows that selection bias effects associated with jet energy loss play an important role in the interpretation of jet substructure measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Improved heavy-ion PID using scintillation light detector with neural network analysis: a Monte Carlo simulation study

The photon collection efficiency of gaseous scintillator detectors varies according to the position of the impinging charged particles in the medium that generates scintillation light. Thus, when impinging particles are distributed over a large area, the intrinsic photon-number resolution of the system is affected by a large variation. This work presents and discusses a method for adjusting the total number of detected photons to account for variation in the photon collection efficiency as a function of the position of the light source within the scintillating medium. The method was developed and validated by processing data from systematic simulation studies based on GEANT4 that model the response of the Energy Loss Optical Scintillation System (ELOSS) detector. The position of the charged particle is calculated using a deep neural network algorithm. This is accomplished by analyzing the distribution of scintillation light recorded by the array of photosensors. The estimated particle position is then used to calculate the correction factor and adjust the amount of captured light to account for variations in the photon collection efficiency. The neural network algorithm provides excellent tracking capabilities, achieving sub-millimeter position resolution and an angular resolution of 12 mrad, approaching the performance of traditional tracking detectors (e.g., drift chambers). The present method can be generalized to any optical scintillation system where the photon collection efficiency depends on the position of the impinging particle.

Heavy-ion detectors↗

The ICARUS light detection system (Fermilab Italian summer student program 2023)

The ICARUS T600 detector is the largest Liquid Argon Time Projection Chamber ever used on a neutrino beam, and it acts as the far detector of the Short Baseline Neutrino program at Fermilab. Its purpose is the study of the possibility of the existence of a fourth sterile neutrino in the O(eV2) mass range that could take part in the neutrino oscillations. The light detection system of the ICARUS detector plays the role of localising in space and time the neutrino interactions taking place inside the detector, collecting the argon scintillation photons generated when an event occurs. This light detection system also has a crucial role in the trigger system and the rejection of the huge amount of background cosmic events, working together with the trigger request signals sent by the beam complex. In August 2023 part of the cables that are used to carry the signals from the photomultipliers to the electronics were replaced with a new model, because a deterioration of the quality of the signal inside the cables had previously been observed. An analysis of the performance of the new cables was carried out, comparing laser run data collected before and after the change of the cable model, and it was proved that the introduction of the new cables lead to a general improvement of the quality of the signals reaching the electronics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗