Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “imaging simulators”

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

Amplitude- and frequency-dependent activation of layer II/III neurons by intracortical microstimulation

Intracortical microstimulation (ICMS) has been used for the development of brain machine interfaces. However, further understanding about the spatiotemporal responses of neurons to different electrical stimulation parameters is necessary to inform the design of optimal therapies. In this study, we employed in vivo electrophysiological recording, two-photon calcium imaging, and electric field simulation to evaluate the acute effect of ICMS on layer II/III neurons. Our results show that stimulation frequency non-linearly modulates neuronal responses, whereas the magnitude of responses is linearly correlated to the electric field strength and stimulation amplitude before reaching a steady state. Temporal dynamics of neurons’ responses depends more on stimulation frequency and their distance to the stimulation electrode. In addition, amplitude-dependent post-stimulation suppression was observed within ~500 μm of the stimulation electrode, as evidenced by both calcium imaging and local field potentials. These findings provide insights for selecting stimulation parameters to achieve desirable spatiotemporal specificity of ICMS.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

NaroNet: Discovery of tumor microenvironment elements from highly multiplexed images

Understanding the spatial interactions between the elements of the tumor microenvironment -i.e. tumor cells. fibroblasts, immune cells- and how these interactions relate to the diagnosis or prognosis of a tumor is one of the goals of computational pathology. We present NaroNet, a deep learning framework that models the multi-scale tumor microenvironment from multiplex-stained cancer tissue images and provides patient-level interpretable predictions using a seamless end-to-end learning pipeline. Trained only with multiplex-stained tissue images and their corresponding patient-level clinical labels, NaroNet unsupervisedly learns which cell phenotypes, cell neighborhoods, and neighborhood interactions have the highest influence to predict the correct label. To this end, NaroNet incorporates several novel and state-of-the-art deep learning techniques, such as patch-level contrastive learning, multi-level graph embeddings, a novel max-sum pooling operation, or a metric that quantifies the relevance that each microenvironment element has in the individual predictions. We validate NaroNet using synthetic data simulating multiplex-immunostained images where a patient label is artificially associated to the -adjustable- probabilistic incidence of different microenvironment elements. We then apply our model to two sets of images of human cancer tissues: 336 seven-color multiplex-immunostained images from 12 high-grade endometrial cancer patients; and 382 35-plex mass cytometry images from 215 breast cancer patients. In both synthetic and real datasets, NaroNet provides outstanding predictions of relevant clinical information while associating those predictions to the presence of specific microenvironment elements.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Multimodal Nanoscale Tomographic Imaging for Battery Electrodes

Accurate representations of the 3D structure within a lithium-ion battery are key to understanding performance limitations. However, obtaining exact reconstructions of electrodes, where the active particles, the carbon black and polymeric binder domain, and the pore space are visualized is challenging. In this work, it is shown that multimodal imaging can be used to overcome this challenge. High-resolution ptychographic X-ray computed tomography are combined with lower resolution but higher contrast transmission X-ray tomographic microscopy to obtain 3D reconstructions of pristine and cycled graphite-silicon composite electrodes. This cross-correlation enables quantitative analysis of the surface of active particles, including the heterogeneity of carbon-black and binder domain and solid-electrolyte interphase coverage. Capturing the active particles as well as the carbon black-binder domain allows using these segmented structures for electrochemical simulations to highlight the influence of the particle embedding on local state of charge heterogeneities.

25 ENERGY STORAGE↗
DOE OSTI2021

Experimental and theoretical insights into the effects of pH on catalysis of bond-cleavage by the lignin peroxidase isozyme H8 from Phanerochaete chrysosporium

Abstract Background Lignin peroxidases catalyze a variety of reactions, resulting in cleavage of both β-O-4′ ether bonds and C–C bonds in lignin, both of which are essential for depolymerizing lignin into fragments amendable to biological or chemical upgrading to valuable products. Studies of the specificity of lignin peroxidases to catalyze these various reactions and the role reaction conditions such as pH play have been limited by the lack of assays that allow quantification of specific bond-breaking events. The subsequent theoretical understanding of the underlying mechanisms by which pH modulates the activity of lignin peroxidases remains nascent. Here, we report on combined experimental and theoretical studies of the effect of pH on the enzyme-catalyzed cleavage of β-O-4′ ether bonds and of C–C bonds by a lignin peroxidase isozyme H8 from Phanerochaete chrysosporium and an acid stabilized variant of the same enzyme. Results Using a nanostructure initiator mass spectrometry assay that provides quantification of bond breaking in a phenolic model lignin dimer we found that catalysis of degradation of the dimer to products by an acid-stabilized variant of lignin peroxidase isozyme H8 increased from 38.4% at pH 5 to 92.5% at pH 2.6. At pH 2.6, the observed product distribution resulted from 65.5% β-O-4′ ether bond cleavage, 27.0% C α -C 1 carbon bond cleavage, and 3.6% C α -oxidation as by-product. Using ab initio molecular dynamic simulations and climbing-image Nudge Elastic Band based transition state searches, we suggest the effect of lower pH is via protonation of aliphatic hydroxyl groups under which extremely acidic conditions resulted in lower energetic barriers for bond-cleavages, particularly β-O-4′ bonds. Conclusion These coupled experimental results and theoretical explanations suggest pH is a key driving force for selective and efficient lignin peroxidase isozyme H8 catalyzed depolymerization of the phenolic lignin dimer and further suggest that engineering of lignin peroxidase isozyme H8 and other enzymes involved in lignin depolymerization should include targeting stability at low pH.

09 BIOMASS FUELS↗
DOE OSTI2021

Dynamics and energy coupling of gas puff Z-pinches on a fast linear transformer driver

Gas puff Z-pinch experiments with annular Ar and Ne gas shells have been conducted on the Compact Experimental System for Z-pinch and Ablation Research (CESZAR) linear transformer driver (LTD) with 500 kA current and 160 ns rise time. Here, we present results from the first systematic gas puff Z-pinch experiments using a fast (≤ 200 ns) LTD as a driver, in which we show that 7% of the stored energy in the capacitors is coupled to plasma kinetic energy as estimated via self-emission and laser schlieren images. 0D and 1D simulations—which do not allow instability growth and thus reach greater maximum average velocities—using initial conditions inferred from experimental implosion trajectories predict coupling in excess of 10% of the stored energy. Additionally, the Ar and Ne implosions were comparably massed and thus achieved similar maximum kinetic energies, though the Ne pinches were more stable and the x-ray pulses were longer and produced higher yield: 2–5 ns and 0.21–0.52 J (0.15–0.37 J/cm) of Ar K-shell and 12–25 ns and 2.2–3.9 J (1.6–2.6 J/cm) of Ne K-shell, respectively. The difference in stability is most likely attributed to variations in initial conditions such as density distribution and gas breakdown initiation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Evaluation of Techniques for Vegetation Removal from UAV-Based Photogrammetric DSM near WIPP Land Withdrawal Act Boundary, NM

The Waste Isolation Pilot Plant is located within New Mexico's karst landscape and is the United States' only deep geologic radioactive waste repository which isolates transuranic waste from defense activities underground in a bedded salt formation. Characteristic karst topography in this region is made up of material (such as limestone) that can be dissolved by water over a long period of time, and includes features such as springs, surface streams, sinking streams, caves, and sinkholes. Karst features can potentially impact the integrity of the waste repository in the future and compromise WIPP's performance, particularly when coupled with incompatible land-use activities within the Land Withdrawal Act (LWA) boundary as well as increased water withdrawals outside the LWA. In order to better capture the ground surface topography to support hydrology modeling efforts, the U.S. DOE has enlisted FIU to develop a high-resolution digital elevation model (DEM). In a previous field study, a high-resolution digital surface model (DSM) of an adjacent representative site, Basin 6, was developed by photogrammetry using aerial images captured with unmanned aerial vehicles. The DSM generated was based on point clouds representing the entire landscapes, including points of terrain, vegetation and infrastructure. These above-ground features need to be removed to create a DSM/DEM with good representation of the bare ground. In recent years, many techniques including both pre- and post-processing of DSMs, were developed for vegetation removal with varied success. This research focuses on applying and comparing various vegetation removal methods using point clouds generated of the Basin 6 pilot study area. The 3D RGB-based point cloud classification using Python scripts was a modification of the procedure by Themistocleous (2019). The method using Liblas classifies the point cloud three-dimensionally, whereas the method by Themistocleous (2019) is two-dimensional. The preliminary result of the method using Liblas, however, still had a large amount of vegetation remaining compared to the Pix4D Machine Learning method. The vegetation layer remains highly visible in the DEMs processed using the vegetation indices as opposed to the Pix4D Machine Learning method, which produces a DEM where mostly bare ground is seen. Further modifications to the input parameters are required to improve the vegetation vs. bare ground classification using this approach. The next step will be to improve the method using Liblas and to also test other vegetation removal methods found in the literature review to compare them against the two described here.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Imaging surface topography with coherent x-ray reflectivity: Theory, kinematics, and simulations

A theoretical formalism is described for understanding coherent x-ray reflectivity (CXR) from the surface of a semi-infinite crystal having a variable surface topography, described by the height profile ℎ(𝑥,𝑦). The surface topography is imaged as a complex “effective density,” obtained from the phasing and inversion of the coherent x-ray reflectivity data, measured through a rocking scan centered at a vertical momentum transfer 𝑄$^{0}_{𝑧}$ and a vertical range Δ⁢𝑄 𝑧 . The formalism predicts that the effective density has an amplitude with a maximum located at the surface height for each position within the surface plane. The phase of the effective density has a lateral variation that is controlled by the surface height and a vertical variation that reflects a combination of the interfacial structure and specific choice of measurement conditions. This understanding enables direct observation of nanometer-scale interfacial topography, i.e., ℎ⁡(𝑥,𝑦)⁢𝑐 𝑠 (where 𝑐 𝑠 is the vertical substrate lattice parameter) with Å-scale sensitivity to surface height. Numerical simulations illustrate and confirm the theoretical results. These results show how the interpretation of the interfacial density phase obtained by CXR data inversion (i.e., surface topography with respect to a flat surface) is conceptually similar to that previously known for Bragg coherent diffraction imaging (BCDI) measurements of isolated nanoparticles (i.e., lattice displacements with respect to an ideal crystal lattice). This suggests that CXR can be thought of as a form of dark field imaging with respect to the bright field BCDI approach. An implication of these results is that interfacial imaging may bypass some of the significant challenges associated with BCDI imaging of multiple particles having different orientations.

X-ray imaging↗
DOE OSTI2020

The impact of mineral reactive surface area variation on simulated mineral reactions and reaction rates

Reactive transport modeling is an essential tool to simulate complex geochemical reactions in porous media that can impact formation properties including porosity and permeability. However, simulating these reactions is challenging due to uncertainties in model parameters, particularly mineral surface areas. Imaging has emerged as a powerful means of estimating model parameters including porosity and mineral abundance, accessibility and accessible surface area. However, these parameters, particularly mineral accessible surface area, vary with image resolution. This work aims to enhance understanding of the impact of image resolution and other means of estimating mineral reactive surface area on simulated mineral reactions and reaction rates. Mineral surface areas calculated from images with resolutions of 0.34 μm and 5.71 μm were used to simulate mineral reactions in the context of geologic CO 2 sequestration in the Paluxy formation at the continuum scale. Additional simulations were carried out using BET surface areas collected from the literature and geometric surface areas. Simulations were run for 7300 days and mineral volume fractions and effluent ion concentrations tracked and compared. Variations in mineral surface areas measured from images are within 1 order of magnitude and yield similar simulation results, indicating the impact of image resolution on simulated reactions and reaction rates is minimum for the resolutions and sample considered. In comparison, surface areas obtained from BET and geometric approaches are 1–4 orders of magnitude higher than image-obtained surface areas and result in greater simulated reaction rates and extents. Minerals with high reaction rates (calcite and siderite) are most impacted by surface area values at short times where simulated mineral volume fractions at longer times agree relatively well, even for simulations with several orders of magnitude variation in surface area. Phases with lower reaction rates, such as K-feldspar and muscovite, are predominantly impacted over longer times where variations in surface areas impact reaction extents and porosity evolution. Overall, variations in surface areas due to image resolution are small and result in little variation in simulated reactions and reaction rates while there are significant variations in simulation results when other surface area estimates are used. These variations, however, depend on the reactivity of the mineral phase where surface areas of fast-reacting phases largely impact simulations at short (hours to days) timescales and surface areas of slower-reactive phases impact longer term simulations (years).

58 GEOSCIENCES↗
DOE OSTI2024

Lattice Symmetry-Guided Charge Transport in 2D Supramolecular Polymers Promotes Triplet Formation

Singlet-to-triplet intersystem crossing (ISC) in organic molecules is intimately connected with their geometries: by modifying the molecular shape, symmetry selection rules pertaining to spin-orbit coupling can be partially relieved, leading to extra matrix elements for increased ISC. As an analog to this molecular design concept, the study finds that the lattice symmetry of supramolecular polymers also defines their triplet formation efficiencies. A supramolecular polymer self-assembled from weakly interacting molecules is considered. Its 2D oblique unit cell effectively renders it as a coplanar array of 1D molecular columns weakly bound to each other. Using momentum-resolved photoluminescence imaging in combination with Monte Carlo simulations, the study found that photogenerated charge carriers in the supramolecular polymer predominantly recombine as spin-uncorrelated carrier pairs through inter-column charge transfer states. This lattice-defined recombination pathway leads to a substantial triplet formation efficiency (≈60%) in the supramolecular polymer. These findings suggest that lattice symmetry of micro-/macroscopic structures relying on intermolecular interactions can be strategized for controlled triplet formation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Atomic Fingerprinting of Heteroatoms Using Noncontact Atomic Force Microscopy

We report immense strides have been made in increasing the resolution of scanning probe microscopy. Noncontact atomic force microscopy (n-AFM) now offers one the ability to characterize and visualize single molecules with subatomic resolution. Specifically, nc-AFM with a carbon monoxide (CO) functionalized tip has the ability to discriminate functional groups (–C≡C–, –CH 2 , –C$=$O, …), although discriminating atomic species often remains as an ongoing challenge. Here, real-space pseudopotentials constructed within density functional theory are employed to accurately simulate nc-AFM images of molecules containing heteroatoms (S, I, and N) within dibenzothiophene (DBT), 2-iodotriphenylene (ITP), acridine (ACR) and ferrous phthalocyanine (FePc). It is found that S and I atoms can be easily identified from C based on their unique features. For N atoms, a use of tip functionalization is proposed to effectively discriminate them from C atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Importance of viewing angle: Hotspot effect improves the ability of satellites to track terrestrial photosynthesis

The product of near-infrared reflectance of vegetation and photosynthetic active radiation (NIRvP) is a new tool for monitoring gross primary productivity (GPP) dynamics in terrestrial ecosystems, due to the discovered linear correlation between NIRvP and GPP. While remote sensing-based NIRvP is considerably influenced by sensor geometry, such geometry impacts on the NIRvP-GPP relationship remain underexplored. In this study, we calculate NIRvP using observations from the Deep Space Climate Observatory (DSCOVR) that provide unique hotspot observation geometry in which the sensor viewing angle coincides with the sun direction. We evaluated the linear correlation between NIRvP and GPP in both the common nadir direction and the special hotspot direction. The results indicate that NIRvP in the hotspot direction significantly outperforms that in the nadir direction for tracking GPP variations across different ecosystems from diurnal to daily scales. This conclusion is further supported by data from the MODerate resolution Imaging Spectroradiometer (MODIS) and simulations using the Soil Canopy Observation Photosynthesis Energy (SCOPE) model. Finally, our research highlights the value of using the unconventional hotspot-based sun-tracking satellite observations for a more accurate characterization of GPP dynamics in terrestrial ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Clarification of creep deformation mechanism in heat-affected zone of 9Cr steels with In Situ experiments

This work quantified nonuniform creep deformation across the heterogeneous heat-affected zone (HAZ) of Grade 91 steel with sophisticated experiments, including an electric-thermal finite element model–assisted Gleeble thermomechanical simulation and a high-temperature creep testing with in situ digital image correlation (DIC). High temperature creep properties of HAZ sub-zones were quantitatively measured by the DIC. Furthermore, by utilizing peak temperature, hardness, local creep strain, and underlying microstructures, creep deformation mechanisms in HAZ were further understood. DIC measurements reveal a creep-vulnerable zone (CVZ) exposed to a peak temperature of 932°C (close to A C3 ) in the intercritical HAZ experienced the fastest creep strength degradation instead of the soft zone with the lowest hardness prior to creep. The significantly reduced precipitation strengthening from misplacement of undissolved and coarsened M 23 C 6 carbides led to a faster recrystallization of tempered martensite in the CVZ. Weak untransformed tempered martensite (ferrite grains) stabilized by local Cr enrichment from dissolved M 23 C 6 also harmed the CVZ's creep resistance.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Competing Intermolecular and Molecule–Surface Interactions: Dipole–Dipole-Driven Patterns in Mixed Carborane Self-Assembled Monolayers

Carboranedithiol isomers adsorbing with opposite orientations of their dipoles on surfaces are self-assembled together to form mixed monolayers where both lateral dipole–dipole and lateral thiol–thiolate (S–H···S) interactions provide enhanced stability over single-component monolayers. Here, we demonstrate the first instance of the ability to map individual isomers in a mixed monolayer using the model system carboranedithiols on Au{111}. The addition of methyl groups to one isomer provides both an enhanced dipole moment and extra apparent height for differentiation via scanning tunneling microscopy (STM). Associated computational investigations rationalize favorable interactions of mixed pairs and the associated stability changes that arise from these interactions. Both STM images and Monte Carlo simulations yield similarly structured mixed monolayers, where approximately 10% of the molecules have reversed dipole moment orientations but no direct chemical attachment to the surface, leading to homogeneous monolayers with no apparent phase separation. Deprotonating the thiols by depositing the molecules under basic conditions eliminates the lateral S–H···S interactions while accentuating the dipole–dipole forces. The molecular system investigated is composed of isomeric molecules with opposite orientations of dipoles and identical surface packing, which enables the mapping of individual molecules within the mixed monolayers and enables analyses of the contributions of the relatively weak lateral interactions to the overall stability of the assemblies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Low Temperature Oxygen Activation on the NiAg(100) Single-Atom Alloy Surface

Silver-catalyzed ethylene epoxidation remains the only industrially viable route for ethylene oxide (EO) production. However, this process requires chlorine and other promoters to achieve a high EO selectivity while still generating substantial CO 2 emissions. A recent theory-guided approach identified Ni, in single-atom alloy (SAA) form, as a new promoter of this reaction. Specifically, the addition of Ni to Ag nanoparticles supported on α-Al 2 O 3 at a highly diluted ratio (1 Ni per 200 Ag atoms) increased catalyst selectivity to EO by ∼25%, the same increase afforded by the ubiquitous industrial promoter chlorine. To better understand the effect of Ni, we investigated the interaction of O 2 with NiAg(100) SAA surfaces by using scanning tunneling microscopy (STM) and density functional theory (DFT). While only molecular O 2 was present when pure Ag(100) was exposed to O 2 at 78 K, a distinct NiO 2 species, indicative of O 2 dissociation at Ni atom sites, was identified on the NiAg(100) SAA under the same conditions. High-resolution STM imaging backed by DFT simulations elucidated the formation of an O−Ni−O species with the oxygen atoms in 4-fold hollow sites. These findings provide direct experimental evidence that Ni atoms are very effective at O 2 activation, even at cryogenic temperatures. This suggests that, in addition to the known role of Ni in stabilizing the unselective nucleophilic oxygen on Ag, it could also accelerate O 2 dissociation, which can be rate limiting.

Catalysts↗
DOE OSTI2024

Magnetic Solitons in Hierarchical 3D Magnetic Nanoarchitectures of Nanoflower Shape

Curvilinear magnetism emerged as a new route to tailor properties of magnetic solitons by the choice of geometry and topology of a magnetic architecture. Here, we develop an anodized aluminum oxide template-based approach to realize hierarchical 3D magnetic nanoarchitectures of nanoflower shape. The technique provides defect-free regular arrays of magnetic nanoflowers of tunable shape with a period of 400 nm over cm 2 areas. We combined advanced magnetic imaging methods with micromagnetic simulations to study complex magnetic states in nanoflowers originating due to magnetostatics-driven symmetry break in curvilinear nanomembranes. An interaction between surface and volume magnetostatic charges in 3D curved nanoflowers leads to the stabilization of asymmetric and shifted vortices as well as states with two Bloch lines. Ordered large area arrays of complex-shaped magnetic nanoarchitectures developed in this work are relevant for prospective research on 3D magnonics and spintronics.

3D nanoarchitectures↗
DOE OSTI2025

Unraveling the Mechanism of Interfacial Charge Transfer and Photoresponsivity of WS 2 Quantum Dots/MoS 2 (0D–2D) Heterostructure-Based Transistors

To assess the efficacy of a mixed-dimensional van der Waals (vdW) heterostructure in modulating the optoelectronic responses of nanodevices, the charge transport properties of the transition-metal dichalcogenide (TMD)-based heterostructure comprising zero-dimensional (0D) WS 2 quantum dots (QDs) and two-dimensional (2D) MoS 2 flakes are critically analyzed. Herein, a facile strategy was materialized in developing an atomically thin phototransistor assembled from mechanically exfoliated MoS 2 and WS 2 QDs synthesized using a one-pot hydrothermal route. The amalgamated photodetectors exhibited a high responsivity of ∼8000 A/W at an incident power of 0.05 nW of white light, surpassing that of the pristine MoS 2 devices. Furthermore, the detectivity of pristine MoS 2 , which was on the order of 10 10 , increased to 10 12 Jones for the WS 2 QDs/MoS 2 heterostructure photodetector, outperforming other WS 2 -based materials. The quasiparticle band gap and density of states (DOS) are further analyzed to elucidate the photophysics of the WS 2 QD/MoS 2 hybrid assembly. The difference in the work function between MoS 2 and WS 2 QDs gives rise to an electric field across the 0D–2D interface, facilitating effective charge separation and migration and contributing to the enhancement of photoresponsivity. The analysis of optical responses using density functional theory (DFT) revealed stronger absorption and less reflection over a broader spectrum of wavelengths for the heterostructure compared to the pristine materials. The estimated optical conductivity aligns well with the experimentally predicted maximum photoresponsivity under visible light, which is attributed to the high absorbance of 2D MoS 2 . In conclusion, combining diverse spectroscopic and imaging techniques with quantum simulation provides insights that clarify the pertinence of 0D–2D TMDs in designing phototransistors.

WS2 QDs/MoS2 heterostructure↗
DOE OSTI2020

Nanoscale Mass Spectrometry Multimodal Imaging via Tip-Enhanced Photothermal Desorption

Materials ranging from adhesives, pharmaceuticals, lubricants, and personal care products are traditionally studied using macroscopic characterization techniques. However, their functionality is in reality defined by details of chemical organization on often noncrystalline matter with characteristic length scales on the order of microns to nanometers. Additionally, these materials are traditionally difficult to analyze using standard vacuum-based approaches that provide nanoscale chemical characterization due to their volatile and beam-sensitive nature. Therefore, approaches that operate under ambient conditions need to be developed that allow probing of nanoscale chemical phenomena and correlated functionality. Here, we demonstrate a tool for probing and visualizing local chemical environments and correlating them to material structure and functionality using advanced multimodal chemical imaging on a combined atomic force microscopy (AFM) and mass spectrometry (MS) system using tip-enhanced photothermal desorption with atmospheric pressure chemical ionization (APCI). We demonstrate enhanced performance metrics of the technique for correlated imaging and point sampling and illustrate the applicability for the analysis of trace chemicals on a human hair, additives in adhesives on paper, and pharmaceuticals samples notoriously difficult to analyze in a vacuum environment. Overall, this approach of correlating local chemical environments to structure and functionality is key to advancing research in many fields ranging from biology, to medicine, to material science.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Radiation image reconstruction and uncertainty quantification using a Gaussian process prior

We propose a complete framework for Bayesian image reconstruction and uncertainty quantification based on a Gaussian process prior (GPP) to overcome limitations of maximum likelihood expectation maximization (ML-EM) image reconstruction algorithm. The prior distribution is constructed with a zero-mean Gaussian process (GP) with a choice of a covariance function, and a link function is used to map the Gaussian process to an image. Unlike many other maximum a posteriori approaches, our method offers highly interpretable hyperparamters that are selected automatically with the empirical Bayes method. Furthermore, the GP covariance function can be modified to incorporate a priori structural priors, enabling multi-modality imaging or contextual data fusion. Lastly, we illustrate that our approach lends itself to Bayesian uncertainty quantification techniques, such as the preconditioned Crank–Nicolson method and the Laplace approximation. The proposed framework is general and can be employed in most radiation image reconstruction problems, and we demonstrate it with simulated free-moving single detector radiation source imaging scenarios. We compare the reconstruction results from GPP and ML-EM, and show that the proposed method can significantly improve the image quality over ML-EM, all the while providing greater understanding of the source distribution via the uncertainty quantification capability. Furthermore, significant improvement of the image quality by incorporating a structural prior is illustrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗