Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reconstruction algorithms”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Predicting the single-site and multi-site event discrimination power of dual-phase time projection chambers

Dual-phase xenon time projection chambers (TPCs) are widely used in searches for rare dark matter and neutrino interactions, in part because of their excellent position reconstruction capability in 3D. Despite their millimeter-scale resolution along the charge drift axis, xenon TPCs face challenges in resolving single-site (SS) and multi-site (MS) interactions in the transverse plane. In this paper, we build a generic TPC model with an idealized signal readout, and use Fisher Information (FI) to predict its theoretical capability of differentiating SS and MS events using the electroluminescence signal. We also demonstrate via simulation that, when only statistical photon noise is present, the theoretical limits can be approached with conventional reconstruction algorithms like maximum likelihood estimation, and with a convolutional neural network classifier. The implications of this study on future TPC experiments will be discussed.

Physics↗
DOE OSTI2021

PRISM3D: a 3-D reference seismic model for Iberia and adjacent areas

SUMMARY We present PRISM3D, a 3-D reference seismic model of P- and S-wave velocities for Iberia and adjacent areas. PRISM3D results from the combination of the most up-to-date earth models available for the region. It extends horizontally from 15°W to 5°E in longitude, 34°N to 46°N in latitude and vertically from 3.5 km above to 200 km below sea level, and is modelled on a regular grid with 10 and 0.5 km of grid node spacing in the horizontal and vertical directions, respectively. It was designed using models inferred from local and teleseismic body-wave tomography, earthquake and ambient noise surface wave tomography, receiver function analysis and active source experiments. It includes two interfaces, namely the topography/bathymetry and the Mohorovičić (Moho) discontinuity. The Moho was modelled from previously published receiver function analysis and deep seismic sounding results. To that end we used a probabilistic surface reconstruction algorithm that allowed to extract the mean of the Moho depth surface along with its associated standard deviation, which provides a depth uncertainty estimate. The Moho depth model is in good agreement with previously published models, although it presents slightly sharper gardients in orogenic areas such as the Pyrenees or the Betic-Rif system. Crustal and mantle P- and S-wave wave speed grids were built separately on each side of the Moho depth surface by weighted average of existing models, thus allowing to realistically render the speed gradients across that interface. The associated weighted standard deviation was also calculated, which provides an uncertainty estimation on the average wave speed values at any point of the grid. At shallow depths (<10 km), low P and S wave speeds and high VP/VS are observed in offshore basins, while the Iberian Massif, which covers a large part of western Iberia, appears characterized by a rather flat Moho, higher than average VP and VS and low VP/VS. Conversely, the Betic-Rif system seems to be associated with low VP and VS, combined with high VP/VS in comparison to the rest of the study area. The most prominent feature of the mantle is the well known high wave speed anomaly related to the Alboran slab imaged in various mantle tomography studies. The consistency of PRISM3D with previous work is verified by comparing it with two recent studies, with which it shows a good general agreement.The impact of the new 3-D model is illustrated through a simple synthetic experiment, which shows that the lateral variations of the wave speed can produce traveltime differences ranging from –1.5 and 1.5 s for P waves and from –2.5 and 2.5 s for S waves at local to regional distances. Such values are far larger than phase picking uncertainties and would likely affect earthquake hypocentral parameter estimations. The new 3-D model thus provides a basis for regional studies including earthquake source studies, Earth structure investigations and geodynamic modelling of Iberia and its surroundings.

Arroucau, P.↗
DOE OSTI2023

Evaluation of Electron Tomography Capabilities for Shale Imaging

Despite the advantageous resolution of electron tomography (ET), reconstruction of three-dimensional (3D) images from multiple two-dimensional (2D) projections presents several challenges, including small signal-to-noise ratios, and a limited projection range. This study evaluates the capabilities of ET for thin sections of shale, a complex nanoporous medium. A numerical phantom with 1.24 nm pixel size is constructed based on the tomographic reconstruction of a Barnett shale. A dataset of 2D projection images is numerically generated from the 3D phantom and studied over a range of conditions. First, common reconstruction techniques are used to reconstruct the shale structure. The reconstruction uncertainty is quantified by comparing overall values of storage and transport metrics, as well as the misclassification of pore voxels compared to the phantom. We then select the most robust reconstruction technique and we vary the acquisition conditions to quantify the effect of artifacts. We find a strong agreement for large pores over the different acquisition workflows, while a wider variability exists for nanometer-scale features. Furthermore, the limited projection range and reconstruction are identified as the main experimental bottlenecks, thereby suggesting that sample thinning, advanced holders, and advanced reconstruction algorithms offer opportunities for improvement.

3D scanning transmission electron microscopy↗
DOE OSTI2023

Validation of semi-analytical, semi-empirical covariance matrices for two-point correlation function for early DESI data

ABSTRACT We present an extended validation of semi-analytical, semi-empirical covariance matrices for the two-point correlation function (2PCF) on simulated catalogs representative of luminous red galaxies (LRGs) data collected during the initial 2 months of operations of the Stage-IV ground-based Dark Energy Spectroscopic Instrument (DESI). We run the pipeline on multiple effective Zel’dovich (EZ) mock galaxy catalogs with the corresponding cuts applied and compare the results with the mock sample covariance to assess the accuracy and its fluctuations. We propose an extension of the previously developed formalism for catalogs processed with standard reconstruction algorithms. We consider methods for comparing covariance matrices in detail, highlighting their interpretation and statistical properties caused by sample variance, in particular, non-trivial expectation values of certain metrics even when the external covariance estimate is perfect. With improved mocks and validation techniques, we confirm a good agreement between our predictions and sample covariance. This allows one to generate covariance matrices for comparable data sets without the need to create numerous mock galaxy catalogs with matching clustering, only requiring 2PCF measurements from the data itself. The code used in this paper is publicly available at https://github.com/oliverphilcox/RascalC.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The Calar Alto Legacy Integral Field Area survey: extended and remastered data release

ABSTRACT This paper describes the extended data release (eDR) of the Calar Alto Legacy Integral Field Area (CALIFA) survey. It comprises science-grade quality data for 895 galaxies obtained with the Potsdam Multi Aperture Spectograph/PPak instrument at the 3.5-m telescope at the Calar Alto Observatory along the last 12 yr, using the V500 setup [3700–7500 Å, 6 Å/full-width at half-maximum (FWHM)] and the CALIFA observing strategy. It includes galaxies of any morphological type, star formation stage, a wide range of stellar masses (∼107–1012 M⊙), at an average redshift of ∼0.015 (90 per cent within 0.005 < z < 0.05). Primarily selected based on the projected size and apparent magnitude, we demonstrate that it can be volume corrected resulting in a statistically limited but representative sample of the population of galaxies in the nearby Universe. All the data were homogeneous re-reduced, introducing a set of modifications to the previous reduction. The most relevant is the development and implementation of a new cube-reconstruction algorithm that provides with an (almost) seeing-limited spatial resolution (FWHMPSF ∼ 1.0 arcsec). To illustrate the usability and quality of the data, we extracted two aperture spectra for each galaxy (central 1.5 arcsec and fully integrated), and analyse them using pyFIT3D. We obtain a set of observational and physical properties of both the stellar populations and the ionized gas, that have been compared for the two apertures, exploring their distributions as a function of the stellar masses and morphologies of the galaxies, comparing with recent results in the literature.

Sánchez, S. F. (ORCID:0000000164449307)↗
DOE OSTI2021

Cosmic Ray Background Rejection with Wire-Cell LArTPC Event Reconstruction in the MicroBooNE Detector

For a large liquid-argon time-projection chamber (LArTPC) operating on or near the Earth’s surface to detect neutrino interactions, the rejection of cosmogenic background is a critical and challenging task because of the large cosmic-ray flux and the long drift time of the time-projection chamber. We introduce a superior cosmic background rejection procedure based on the Wire-Cell three-dimensional (3D) event reconstruction for LArTPCs. From an initial 1:20 000 neutrino to cosmic-ray background ratio, we demonstrate these tools on data from the MicroBooNE experiment and create a high-performance generic neutrino event selection with a cosmic contamination of 14.9% (9.7%) for a visible energy region greater than O(200) MeV. The neutrino interaction selection efficiency is 80.4% and 87.6% for inclusive ν μ charged-current and ν e charged-current interactions, respectively. Here, this significantly improved performance compared with existing reconstruction algorithms marks a major milestone toward reaching the scientific goals of LArTPC neutrino oscillation experiments operating near the Earth’s surface.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Adversarial methods to reduce simulation bias in neutrino interaction event filtering at liquid argon time projection chambers

For current and future neutrino oscillation experiments using large liquid argon time projection chambers (LAr-TPCs), a key challenge is identifying neutrino interactions from the pervading cosmic-ray background. Rejection of such background is often possible using traditional cut-based selections, but this typically requires the prior use of computationally expensive reconstruction algorithms. This work demonstrates an alternative approach of using a 3D submanifold sparse convolutional network trained on low-level information from the scintillation light signal of interactions inside LAr-TPCs. This technique is applied to example simulations from ICARUS, the far detector of the short baseline neutrino program at Fermilab. The results of the network, show that cosmic background is reduced by up to 76.3% whilst neutrino interaction selection efficiency remains over 98.9%. We further present a way to mitigate potential biases from imperfect input simulations by applying domain adversarial neural networks (DANNs), for which modified simulated samples are introduced to imitate real data and a small portion of them are used for adversarial training. A series of mock-data studies are performed and demonstrate the effectiveness of using DANNs to mitigate biases, showing neutrino interaction selection efficiency performances significantly better than that achieved without the adversarial training.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Measurement of the depth of maximum of air-shower profiles with energies between 10 18.5 and 10 20 eV using the surface detector of the Pierre Auger Observatory and deep learning

We report an investigation of the mass composition of cosmic rays with energies from 3 to 100 EeV ( 1 EeV = 10 18 eV ) using the distributions of the depth of shower maximum X max . The analysis relies on ∼ 50 , 000 events recorded by the surface detector of the Pierre Auger Observatory and a deep-learning-based reconstruction algorithm. Above energies of 5 EeV, the dataset offers a 10-fold increase in statistics with respect to fluorescence measurements at the Observatory. After cross-calibration using the fluorescence detector, this enables the first measurement of the evolution of the mean and the standard deviation of the X max distributions up to 100 EeV. Our findings are threefold: (i) The evolution of the mean logarithmic mass toward a heavier composition with increasing energy can be confirmed and is extended to 100 EeV. (ii) The evolution of the fluctuations of X max toward a heavier and purer composition with increasing energy can be confirmed with high statistics. We report a rather heavy composition and small fluctuations in X max at the highest energies. (iii) We find indications for a characteristic structure beyond a constant change in the mean logarithmic mass, featuring three breaks that are observed in proximity to the ankle, instep, and suppression features in the energy spectrum. Published by the American Physical Society 2025

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Solving key challenges in collider physics with foundation models

Foundation models are neural networks that are capable of simultaneously solving many problems. Large language foundation models like ChatGPT have revolutionized many aspects of daily life, but their impact for science is not yet clear. In this paper, we use a new foundation model for hadronic jets to solve three key challenges in collider physics. In particular, we show how experiments can (1) save significant computing power when developing reconstruction algorithms, (2) perform a complete uncertainty quantification for high-dimensional measurements, and (3) search for new physics with model agnostic methods using low-level inputs. In each case, there are significant computational or methodological challenges with current methods that limit the science potential of deep learning algorithms. By solving each problem, we take jet foundation models beyond proof-of-principle studies and into the toolkit of practitioners.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Proton radiography inversions with source extraction and comparison to mesh methods

In this article, proton radiography is a central diagnostic technique for measuring electromagnetic (EM) fields in high-energy-density, laser-produced plasmas. In this technique, protons traverse the plasma where they accumulate small EM deflections which lead to variations in the proton fluence pattern on a detector. Path-integrated EM fields can then be extracted from the fluence image through an inversion process. In this work, experiments of laser-driven foils were conducted on the OMEGA laser and magnetic field reconstructions were performed using both “fluence-based” techniques and high-fidelity “mesh-based” methods. We implement nonzero boundary conditions into the inversion and show their importance by comparing against mesh measurements. Good agreement between the methods is found only when nonzero boundary conditions are used. We also introduce an approach to determine the unperturbed proton source profile, which is a required input in fluence reconstruction algorithms. In this approach, a fluence inversion is embedded inside of a mesh region, which provides overconstrained magnetic boundary conditions. A source profile is then iteratively optimized to satisfy the boundary information. This method substantially enhances the accuracy in recovering EM fields. Lastly, we propose a scheme to quantify uncertainty in the final inversion that is introduced through errors in the source retrieval.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Analyzer-free linear dichroic ptychography

Linear dichroism is an important tool to characterize the transmission matrix and determine the crystal or orbital orientation in a material. To achieve high-resolution mapping of transmission properties, we introduce the linear-dichroism scattering model in ptychographic imaging and develop an efficient two-stage reconstruction algorithm. Using the proposed algorithm on a uniaxial material, the dichroic transmission matrix can be recovered without an analyzer by using ptychography measurements with as few as three different polarization angles, with the help of an empty region to remove phase ambiguities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Upscaling X-ray nanoimaging to macroscopic specimens

Upscaling X-ray nanoimaging to macroscopic specimens has the potential for providing insights across multiple length scales, but its feasibility has long been an open question. By combining the imaging requirements and existing proof-of-principle examples in large-specimen preparation, data acquisition and reconstruction algorithms, the authors provide imaging time estimates for howX-ray nanoimaging can be scaled to macroscopic specimens. To arrive at this estimate, a phase contrast imaging model that includes plural scattering effects is used to calculate the required exposure and corresponding radiation dose. The coherent X-ray flux anticipated from upcoming diffraction-limited light sources is then considered. This imaging time estimation is in particular applied to the case of the connectomes of whole mouse brains. To image the connectome of the whole mouse brain, electron microscopy connectomics might require years, whereas optimized X-ray microscopy connectomics could reduce this to one week. Furthermore, this analysis points to challenges that need to be overcome (such as increased X-ray detector frame rate) and opportunities that advances in artificial-intelligence-based `smart' scanning might provide. While the technical advances required are daunting, it is shown that X-ray microscopy is indeed potentially applicable to nanoimaging of millimetre- or even centimetre-size specimens.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Depth-resolved Laue microdiffraction with coded apertures

A rapid data acquisition and reconstruction method is introduced to image the crystalline structure of materials and the associated strain and orientations at micrometre resolution using Laue diffraction. The method relies on scanning a coded aperture across the diffracted X-ray beam from broadband illumination and a reconstruction algorithm to resolve Laue microdiffraction patterns as a function of depth along the incident illumination path. It provides rapid access to full diffraction information for sub-micrometre volume elements in bulk materials. Both the theory and the experimental validation of this imaging approach are presented in this report.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Optimizing Error-Bounded Lossy Compression for Scientific Data on GPUs

Error-bounded lossy compression is a critical technique for significantly reducing scientific data volumes. With ever-emerging heterogeneous high-performance computing (HPC) architecture, GPU-accelerated error-bounded compressors (such as CUSZ and cuZFP) have been developed. However, they suffer from either low performance or low compression ratios. To this end, we propose CUSZ+ to target both high compression ratios and throughputs. We identify that data sparsity and data smoothness are key factors for high compression throughputs. Our key contributions in this work are fourfold: (1) We propose an efficient compression workflow to adaptively perform run-length encoding and/or variable-length encoding. (2) We derive Lorenzo reconstruction in decompression as multidimensional partial-sum computation and propose a fine-grained Lorenzo reconstruction algorithm for GPU architectures. (3) We carefully optimize each of CUSZ kernels by leveraging state-of-the-art CUDA parallel primitives. (4) We evaluate CUSZ+ using seven real-world HPC application datasets on V100 and A100 GPUs. Experiments show CUSZ+ improves the compression throughputs and ratios by up to 18.4x and 5.3x, respectively, over CUSZ on the tested datasets.

Tian, Jiannan↗
DOE OSTI2025

Individual and Simultaneous Imaging of ⁹⁹mTc and ¹⁷⁷Lu With a Preclinical Broad Energy-Spectrum CZT-Based SPECT

Radiopharmaceutical therapy has demonstrated a high efficacy in the treatment of various tumor types. One of the radionuclides already used in the clinic is 177Lu, a beta emitter that also emits several photons imageable with SPECT. Quantitative imaging of 177Lu is critical for developing new radiopharmaceuticals. Energy resolution is an important factor when imaging multiple photon emissions. Solid-state detectors offer a superior performance over scintillators, that are commonly used in commercially-available preclinical SPECT scanners. This study demonstrates the feasibility of 99m Tc and 177Lu quantitative imaging in mouse phantoms, individually and simultaneously, with a SPECT prototype built with four CdZnTe (CZT) detector heads and a custom-designed and energy-optimized parallel-hole tungsten collimator. With a custom implementation of the one-step late (OSL) image reconstruction algorithm, the system is capable of imaging energies from ~70 keV to 250 keV. Above 250 keV, images were significantly affected by septal penetration, consistent with the collimator design. A recovery coefficient within 25% was obtained for activities as low as 2 kBq/mL for 99m Tc and 45% for 177Lu. Compared to a commercial NaI-based preclinical SPECT (VECTor4/CT), our prototype showed a superior energy resolution (< 5% at 140 keV), a similar uniformity with a high-compact design.

Encarnação, Pedro M C C↗
DOE OSTI2020

The Single-Volume Scatter Camera

The multi-institution Single-Volume Scatter Camera (SVSC) collaboration led by Sandia National Laboratories (SNL) is developing a compact, high-efficiency double-scatter neutron imaging system. Kinematic emission imaging of fission-energy neutrons can be used to detect, locate, and spatially characterize special nuclear material. Neutron-scatter cameras, analogous to Compton imagers for gamma ray detection, have a wide field of view, good event-by-event angular resolution, and spectral sensitivity. Existing systems, however, suffer from large size and/or poor efficiency. We are developing high-efficiency scatter cameras with small form factors by detecting both neutron scatters in a compact active volume. This effort requires development and characterization of individual system components, namely fast organic scintillators, photodetectors, electronics, and reconstruction algorithms. In this presentation, we will focus on characterization measurements of several SVSC candidate scintillators. The SVSC collaboration is investigating two system concepts: the monolithic design in which isotropically emitted photons are detected on the sides of the volume, and the optically segmented design in which scintillation light is channeled along scintillator bars to segmented photodetector readout. For each of these approaches, we will describe the construction and performance of prototype systems. We will conclude by summarizing lessons learned, comparing and contrasting the two system designs, and outlining plans for the next iteration of prototype design and construction.

Manfredi, Juan↗
DOE OSTI2020

Fast widefield imaging of neuronal structure and function with optical sectioning in vivo

Optical microscopy, owing to its noninvasiveness and subcellular resolution, enables in vivo visualization of neuronal structure and function in the physiological context. Optical-sectioning structured illumination microscopy (OS-SIM) is a widefield fluorescence imaging technique that uses structured illumination patterns to encode in-focus structures and optically sections 3D samples. However, its application to in vivo imaging has been limited. In this study, we optimized OS-SIM for in vivo neural imaging. We modified OS-SIM reconstruction algorithms to improve signal-to-noise ratio and correct motion-induced artifacts in live samples. Incorporating an adaptive optics (AO) module to OS-SIM, we found that correcting sample-induced optical aberrations was essential for achieving accurate structural and functional characterizations in vivo. With AO OS-SIM, we demonstrated fast, high-resolution in vivo imaging with optical sectioning for structural imaging of mouse cortical neurons and zebrafish larval motor neurons, and functional imaging of quantal synaptic transmission at Drosophila larval neuromuscular junctions.

42 ENGINEERING↗
DOE OSTI2023

Highly multicolored light-emitting arrays for compressive spectroscopy

Miniaturized, multicolored light-emitting device arrays are promising for applications in sensing, imaging, computing, and more, but the range of emission colors achievable by a conventional light-emitting diode is limited by material or device constraints. In this work, we demonstrate a highly multicolored light-emitting array with 49 different, individually addressable colors on a single chip. The array consists of pulsed-driven metal-oxide-semiconductor capacitors, which generate electroluminescence from microdispensed materials spanning a diverse range of colors and spectral shapes, enabling facile generation of arbitrary light spectra across a broad wavelength range (400 to 1400 nm). When combined with compressive reconstruction algorithms, these arrays can be used to perform spectroscopic measurements in a compact manner without diffractive optics. As an example, we demonstrate microscale spectral imaging of samples using a multiplexed electroluminescent array in conjunction with a monochrome camera.

42 ENGINEERING↗