Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “peak finding”

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 19 records

Gaussian Process Classification for Galaxy Blend Identification in LSST

Abstract A significant fraction of observed galaxies in the Rubin Observatory Legacy Survey of Space and Time (LSST) will overlap at least one other galaxy along the same line of sight, in a so-called “blend.” The current standard method of assessing blend likelihood in LSST images relies on counting up the number of intensity peaks in the smoothed image of a blend candidate, but the reliability of this procedure has not yet been comprehensively studied. Here we construct a realistic distribution of blended and unblended galaxies through high-fidelity simulations of LSST-like images, and from this we examine the blend classification accuracy of the standard peak-finding method. Furthermore, we develop a novel Gaussian process blend classifier model, and show that this classifier is competitive with both the peak finding method as well as with a convolutional neural network model. Finally, whereas the peak-finding method does not naturally assign probabilities to its classification estimates, the Gaussian process model does, and we show that the Gaussian process classification probabilities are generally reliable.

79 ASTRONOMY AND ASTROPHYSICS↗

Beam Loss Assessment Through Use of Photomultiplier Tubes

The first machine in the Fermilab Accelerator chain, the Linac, delivers a 400 MeV proton beam. The first portion of the Fermilab Linac, the Drift Tube Linac, lacks the degree of instrumentation necessary for beam tuning. To compensate for this, photomultiplier tube (PMT s) based loss monitors were installed on either side of the first two drift tube tanks, but are not yet operational. One of the main goals in this is to make PMT's operational beam loss monitors for tuning. Noise reduction and peak finding on the PMT data is a requirement for this. A method for noise reduction and peak finding has been developed and implemented to produce a consistent and stable output. Future work includes integration with ACNET to automate input and output of data for analysis of beam loss.

Waggoner, Alexander↗

End-to-end deep learning pipeline for real-time Bragg peak segmentation: from training to large-scale deployment

X-ray crystallography reconstruction, which transforms discrete X-ray diffraction patterns into three-dimensional molecular structures, relies critically on accurate Bragg peak finding for structure determination. As X-ray free electron laser (XFEL) facilities advance toward MHz data rates (1 million images per second), traditional peak finding algorithms that require manual parameter tuning or exhaustive grid searches across multiple experiments become increasingly impractical. While deep learning approaches offer promising solutions, their deployment in high-throughput environments presents significant challenges in automated dataset labeling, model scalability, edge deployment efficiency, and distributed inference capabilities. We present an end-to-end deep learning pipeline with three key components: (1) a data engine that combines traditional algorithms with our peak matching algorithm to generate high-quality training data at scale, (2) a modular architecture that scales from a few million to hundreds of million parameters, enabling us to train large expert-level models offline while deploying smaller, distilled models at the edge, and (3) a decoupled producer-consumer architecture that separates specialized data source layer from model inference, enabling flexible deployment across diverse computing environments. Using this integrated approach, our pipeline achieves accuracy comparable to traditional methods tuned by human experts while eliminating the need for experiment-specific parameter tuning. Although current throughput requires optimization for MHz facilities, our system's scalable architecture and demonstrated model compression capabilities provide a foundation for future high-throughput XFEL deployments.

Wang, Cong↗

Detecting Dwarf Galaxies Around NGC 55

In order to improve our understanding of the satellite populations of low-mass galaxies, the DECam Local Volume Exploration (DELVE) – DEEP Survey is performing a search for satellites around the isolated, low-mass galaxy, NGC 55. As part of this search, we characterize our dwarf galaxy detection sensitivity by injecting artificial dwarf galaxies into our data set and measuring the fraction that we can recover. Additionally, we inject and recover artificial dwarfs to test how the blending of resolved stars affects the detection of dwarf galaxies. We inject 650 artificial dwarfs across 26 DELVE-DEEP coadded images of the NGC 55 halo and use the software tool, Source Extractor, to perform source recovery. A dwarf-search pipeline is then run on our dwarf-injected catalogs, utilizing the photutils function, find peaks, and detections are recorded as a function of dwarf size and luminosity. Our results demonstrate that we are most sensitive to large, bright dwarf galaxies, and our sensitivity begins to decrease as we approach DELVE-DEEP’s detection limit or blending causes significant reductions in recovered magnitude. We also find that our results are hindered by the limitations of star-galaxy separation applied to distant, faint sources, which is possibly a consequence of the blending of small star clusters within dwarf galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

mzapy : An Open-Source Python Library Enabling Efficient Extraction and Processing of Ion Mobility Spectrometry-Mass Spectrometry Data in the MZA File Format

We have recently reported MZA, a new and simple mass spectrometry data structure based on the broadly supported HDF5 format and created to facilitate software development. While this format is inherently supportive of application development, the availability of a core library with standard mass spectrometry utilities greatly facilitates fast software development. Here, we present a Python library, mzapy, for efficient extraction and processing of mass spectrometry data in the MZA format. In addition to raw data extraction, mzapy contains supporting utilities enabling tasks including calibration, signal processing, peak finding, and generating plots. Being implemented in pure Python with minimal and largely standardized dependencies makes mzapy uniquely suited to application development in the multi-omics domain. The free and open source mzapy is built with extensibility in mind, and future development will support cloud computing and artificial intelligence/machine learning applications. The software source code is freely available at https://github.com/PNNL-m-q/mzapy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Prototype acoustic positioning system for the Pacific Ocean Neutrino Experiment

We present the design and initial performance characterization of the prototype acoustic positioning system intended for the Pacific Ocean Neutrino Experiment. It comprises novel piezo-acoustic receivers with dedicated filtering- and amplification electronics installed in P-ONE instruments and is complemented by a commercial system comprised of cabled and autonomous acoustic pingers for sub-sea installation manufactured by Sonardyne Ltd. We performed an in-depth characterization of the acoustic receiver electronics and their acoustic sensitivity when integrated into P-ONE pressure housings. These show absolute sensitivities of up to -125 dB re V2/μPa2 in a frequency range of 10–40 kHz. We furthermore conducted a positioning measurement campaign in the ocean by deploying three autonomous acoustic pingers on the seafloor, as well as a cabled acoustic interrogator and a P-ONE prototype module deployed from a ship. Using a simple peak-finding detection algorithm, we observe high accuracy in the tracking of relative ranging times at approximately 230–280 μs at distances of up to 1600 m, which is sufficient for positioning detectors in a cubic-kilometer detector and which can be further improved with more involved detection algorithms. The tracking accuracy is further confirmed by independent ranging of the Sonardyne system and closely follows the ship's drift in the wind measured by GPS. The absolute positioning shows the same tracking accuracy with its absolute precision only limited by the large uncertainties of the deployed pinger positions on the seafloor.

Data analysis↗

Pure spin current injection of single-layer monochalcogenides

We compute the spectrum of pure spin current injection in ferroelectric single-layer SnS, SnSe, GeS, and GeSe. The formalism takes into account the coherent spin dynamics of optically excited conduction states split in energy by spin–orbit coupling. The velocity of the electron's spins is calculated as a function of incoming photon energy and angle of linearly polarized light within a full electronic band structure scheme using density functional theory. We find peak speeds of 520, 360, 270 and 370 Km s -1 for SnS, SnSe, GeS and GeSe, respectively which are an order of magnitude larger than those found in bulk semiconductors, e.g., GaAs and CdSe. Interestingly, the spin velocity is almost independent of the direction of polarization of light in a range of photon energies. Our results demonstrate that single-layer SnS, SnSe, GeS and GeSe are candidates to produce on demand spin-current in spintronics applications.

2D-monochalcogenides↗

Deep learning-based spatiotemporal multi-event reconstruction for delay line detectors

Abstract Accurate observation of two or more particles within a very narrow time window has always been a challenge in modern physics. It creates the possibility of correlation experiments, such as the ground-breaking Hanbury Brown–Twiss experiment, leading to new physical insights. For low-energy electrons, one possibility is to use a Microchannel plate with subsequent delay lines for the readout of the incident particle hits, a setup called a Delay Line Detector. The spatial and temporal coordinates of more than one particle can be fully reconstructed outside a region called the dead radius. For interesting events, where two electrons are close in space and time, the determination of the individual positions of the electrons requires elaborate peak finding algorithms. While classical methods work well with single particle hits, they fail to identify and reconstruct events caused by multiple nearby particles. To address this challenge, we present a new spatiotemporal machine learning model to identify and reconstruct the position and time of such multi-hit particle signals. This model achieves a much better resolution for nearby particle hits compared to the classical approach, removing some of the artifacts and reducing the dead radius a factor of eight. We show that machine learning models can be effective in improving the spatiotemporal performance of delay line detectors.

Computer Science↗

Neural architecture codesign for fast physics applications

We develop a pipeline to streamline neural architecture codesign for physics applications to reduce the need for ML expertise when designing models for novel tasks. Our method employs neural architecture search and network compression in a two-stage approach to discover hardware efficient models. This approach consists of a global search stage that explores a wide range of architectures while considering hardware constraints, followed by a local search stage that fine-tunes and compresses the most promising candidates. We exceed performance on various tasks and show further speedup through model compression techniques such as quantization-aware-training and neural network pruning. We synthesize the optimal models to high level synthesis code for FPGA deployment with the hls4ml library. Additionally, our hierarchical search space provides greater flexibility in optimization, which can easily extend to other tasks and domains. We demonstrate this with two case studies: Bragg peak finding in materials science and jet classification in high energy physics, achieving models with improved accuracy, smaller latencies, or reduced resource utilization relative to the baseline models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Preliminary Thermoluminescent Dosimeter Glow Curve Analysis with Automated Glow Peak Identification for LiF Mg,Ti

When appropriately analyzed, thermoluminescent dosimeter glow curve analysis allows for improved quantification of thermoluminescent material behavior while flagging abnormalities. The mathematical separation of a glow curve into contributions from energetically unique trap states, or glow curve analysis, may be used to remove undesired effects of signal fading for complex materials. A generalized glow curve analysis software for the separation of glow curves is presented in this paper. Written in C ++ , the software uses the first-order kinetics model with automatic peak identification. The automatic identification of peaks is achieved through a unique peak-finding algorithm. Here, the program was performance tested using experimental glow curve data from LiF:Mg,Ti, and comparative results are presented.

47 OTHER INSTRUMENTATION↗

Inferring the shape of data: a probabilistic framework for analysing experiments in the natural sciences

A critical step in data analysis for many different types of experiments is the identification of features with theoretically defined shapes in N -dimensional datasets; examples of this process include finding peaks in multi-dimensional molecular spectra or emitters in fluorescence microscopy images. Identifying such features involves determining if the overall shape of the data is consistent with an expected shape; however, it is generally unclear how to quantitatively make this determination. In practice, many analysis methods employ subjective, heuristic approaches, which complicates the validation of any ensuing results—especially as the amount and dimensionality of the data increase. Here, we present a probabilistic solution to this problem by using Bayes’ rule to calculate the probability that the data have any one of several potential shapes. This probabilistic approach may be used to objectively compare how well different theories describe a dataset, identify changes between datasets and detect features within data using a corollary method called Bayesian Inference-based Template Search; several proof-of-principle examples are provided. Altogether, this mathematical framework serves as an automated ‘engine’ capable of computationally executing analysis decisions currently made by visual inspection across the sciences.

Science & Technology - Other Topics↗

Mu2e - Extinction Monitor Research & Development

Current efforts are being conducted at Fermi National Laboratory to study potential violations in accepted theory that would otherwise suggest a restructuring of our fundamental understanding of the universe. Mu2e is one of these frontier projects that studies charged lepton flavor violation (CLFV) which if observed, would suggest physics beyond the Standard Model. Therefore, this note encompasses several projects that contribute to the fruition of Mu2e investigations. Due to the broad range of disciplinary inconsistencies that each project requires, all the work is being presented as a means of justifying contribution to Mu2e. The projects are comprised of a G4Beamline simulation analyzing 8GeV proton beam interaction with a titanium window of Recycler ring extinction rates using three Cherenkov radiation-based detectors and supplemental work for the implementation of a micro–Telecommunications Computing Architecture (MicroTCA) crate to establish a peak finding algorithm to ensure that the out-of-time beam is less than 10^-10 fractional level along with inefficiency analysis on scintillation counters for the Cosmic-Ray Veto (CRV) analysis. Preliminary results have been achieved for extinction rate simulation by achieving coincidence rates for 2/3-fold and 3/3-fold on the detectors in the order of 10^-9 and 10^-10, respectively. Only preliminary results of a triangular counter and four rectangular di-counters for the CRV have been realized but other non-experimental contributions were made to the development of the MicroTCA crate.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Mu2e - Research & Development

Current efforts are being conducted at Fermi National Laboratory to study potential violations in accepted theory that would otherwise suggest a restructuring of our fundamental understanding of the universe. Mu2e is one of these frontier projects that studies Charged Lepton Flavor Violation (CLFV) which if observed, would suggest physics beyond the Standard Model. Therefore, this note encompasses several projects that contribute to the fruition of Mu2e investigations. Due to the broad range of disciplinary inconsistencies that each project requires, all the work is being presented as a means of justifying contribution to Mu2e. The projects are comprised of a simulation exploring the extinction level of proton pulses after Recycler ring re-bunching by using G4beamline to simulate an 8GeV proton beam interaction with a titanium target, three Cherenkov radiation-based detectors and 2/3-fold and 3/3-fold coincidence rate analysis. Additionally, supplemental work for the implementation of a Micro Telecommunications Computing Architecture (TCA) crate to establish a peak finding algorithm to ensure that the out-of-time beam is less than 10$^{−10}$ fractional level along with single-layer inefficiency analysis on scintillation counters for the Cosmic-Ray Veto (CRV) analysis to ensure the overall inefficiency is 10$^{−4}$. Preliminary results have been achieved for the G4beamline simulation 2/3-fold and 3/3-fold coincidences which are in the order of 10$^{−9}$ and 10$^{−10}$, respectively. Only preliminary results of a triangular counter and four rectangular di-counters for the CRV have been realized. The microTCA crate development is still ongoing.

43 PARTICLE ACCELERATORS↗

Efficient analysis routines for single and double peaked Type 2 AGN spectra

Driven by the imminent need to rapidly process and classify millions of AGN spectra drawn from next generation astronomical facilities, we present a spectral fitting routine for Type 2 AGN spectra optimized for high volume processing, using the quasar spectral fitting library (qsfit). We analyse an optically selected sample of 813 luminous Type 2 AGN spectra at z < 0.83 from the Sloan Digital Sky Survey (SDSS) to qualify its performance. We report a median narrow line H α/H β Balmer decrement of 4.5 ± 0.8, alluding to the presence of dust in the narrow line region (NLR). We publish a specialized QSFIT fitting routine for high signal-to-noise ratio spectra and general fitting routine for double peaked Type 2 AGN spectra applied on a subsample of 45 spectra from our parent sample. We report a median red and blue peak velocity separation of 390 ± 60 kms −1 . No trend is found for red or blue peaks to exhibit systematically different luminosity or ionization properties. Emission line diagnostics show that the double peaks in all sources are illuminated by an AGN-powered ionizing continuum. Finally, we examine the morphology of host galaxies of our double peaked sample. We find double peaked Type 2 AGN reside in merging systems at a comparable frequency to single peaked AGN. This suggests that the double peaked AGN phenomenon is likely to have a bi-conical outflow origin in the majority of cases. We publicly release the code used for spectral analysis and produced catalogues used in this work.

79 ASTRONOMY AND ASTROPHYSICS↗

Peak Reduction Using Mode Adjustment of Heat Pump Water Heaters in a Residential Neighborhood

Building electrification is putting pressure on distribution grid worldwide. Peak reduction is an important concern that can help reduce the growing stress and allow to defer investments in new capacity. Water heaters represent a convenient way of reducing peak because they are less dependent on weather, and their storage volume allows for asynchronous water heating and hot water use. Previous empirical studies investigated the ability of water heaters to reduce peak through the adjustment of the temperature setpoint. However, not all equipment vendors offer this option. This study aims at understanding the feasibility of peak reduction with an alternative configuration available in the market - by adjusting the device mode rather than temperature setpoint. The peak reduction methodology is tested in an occupied 46-townhome neighborhood located in Atlanta, GA. We find that peak shifting is possible with the adjustable mode approach, with the change in the peak load by 30-60%.

demand response↗

Origin and Suppression of Beam Damage-Induced Oxygen-K Edge Artifact from γ-Al 2 O 3 using Cryo-EELS

Gamma-alumina (γ-Al 2 O 3 ), like other low-Z oxides, is readily damaged when exposed to an electron beam. This typically results in the formation of a characteristic pre-edge peak in the oxygen-K edge of electron energy-loss spectra (EELS) acquired during or after the damage process. This artifact can mask the presence of intrinsic O-K edge fine structure that would reveal chemical properties of the material; therefore, its suppression is key. In this work, we systematically investigate the conditions that give rise to the damage-induced O-K pre-edge peak and show that it can be effectively suppressed by performing EELS experiments at cryogenic (cryo) temperatures. Prolonged exposure of γ-Al 2 O 3 to a focused electron beam results in a hole bored through the sample; this was used as a reproducible beam damage condition. O-K edge EELS spectra were collected from a single-crystal γ-Al 2 O 3 sample both during and after focused 2 electron beam hole drilling, and at room and cryo temperatures, using a monochromated scanning transmission electron microscope (STEM). Furthermore, the characteristic 531 eV pre-edge peak visible in the room temperature EELS spectra was completely suppressed in the cryo-EELS spectra, even in the presence of a visible drilled hole. We then correlated these experimental observations with multiple-scattering EELS simulations to determine the likely atomistic origin of the damage-induced O-K pre-edge peak. The findings indicate that the pre-edge peak is caused primarily by the presence of surface O dimer (O-O) bonds formed during beam damage, and that operating at cryo temperature suppresses the formation of surface O-O bonds, thus preventing formation of the O-K pre-edge peak. Additionally, Al-L 2,3 edge EELS spectra revealed Al loss primarily from tetrahedral sites during hole drilling.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Magic Gap Ratio for Optimally Robust Fermionic Condensation and Its Implications for High- $T_c$ Superconductivity

Bardeen-Schrieffer-Cooper (BCS) and Bose-Einstein condensation (BEC) occur at opposite limits of a continuum of pairing interaction strength between fermions. A crossover between these limits is readily observed in a cold atomic Fermi gas. Whether it occurs in other systems such as the high temperature superconducting cuprates has remained an open question. We uncover here unambiguous evidence for a BCS-BEC crossover in the cuprates by identifying a universal magic gap ratio 2Δ/k B T c ≈ 6.5 (where Δ is the pairing gap and T c is the transition temperature) at which paired fermion condensates become optimally robust. At this gap ratio, corresponding to the unitary point in a cold atomic Fermi gas, the measured condensate fraction N0 and the height of the jump δγ(T c ) in the coefficient γ of the fermionic specific heat at T c are strongly peaked. In the cuprates, δγ(T c ) is peaked at this gap ratio when Δ corresponds to the antinodal spectroscopic gap, thus reinforcing its interpretation as the pairing gap. We find the peak in δγ(T c ) also to coincide with a normal state maximum in γ, which is indicative of a pairing fluctuation pseudogap above T c .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

TARDIS. II. Synergistic Density Reconstruction from Lyα Forest and Spectroscopic Galaxy Surveys with Applications to Protoclusters and the Cosmic Web

Looking at this work, we expand upon the Tomographic Absorption Reconstruction and Density Inference Scheme (TARDIS) in order to include multiple tracers while reconstructing matter density fields at Cosmic Noon (z ~ 2–3). In particular, we jointly reconstruct the underlying density field from simulated Lyα forest observations at z ~ 2.5 and an overlapping galaxy survey. We find that these data are synergistic, with the Lyα forest providing reconstruction of low-density regions and galaxy surveys tracing the density peaks. We find a more accurate power spectra reconstruction going to higher scales when fitting these two data sets simultaneously than when using either one individually. When applied to cosmic web analysis, we find that performing the joint analysis is equivalent to an Lγα survey with significantly increased sight-line spacing. Because we reconstruct the velocity field and matter field jointly, we demonstrate the ability to evolve the mock observed volume further to z = 0, allowing us to create a rigorous definition of a “protocluster” as regions that will evolve into clusters. We apply our reconstructions to study protocluster structure and evolution, finding for realistic survey parameters that we can provide accurate mass estimates of the z ≈ 2 structures and their z = 0 fate.

79 ASTRONOMY AND ASTROPHYSICS↗