SEARCH · Search NASA

Results for “shape classification”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Whole-Genome Comparisons of Ergot Fungi Reveals the Divergence and Evolution of Species within the Genus Claviceps Are the Result of Varying Mechanisms Driving Genome Evolution and Host Range Expansion

The genus Claviceps has been known for centuries as an economically important fungal genus for pharmacology and agricultural research. Only recently have researchers begun to unravel the evolutionary history of the genus, with origins in South America and classification of four distinct sections through ecological, morphological, and metabolic features (Claviceps sects. Citrinae, Paspalorum, Pusillae, and Claviceps). The first three sections are additionally characterized by narrow host range, whereas section Claviceps is considered evolutionarily more successful and adaptable as it has the largest host range and biogeographical distribution. However, the reasons for this success and adaptability remain unclear. Our study elucidates factors influencing adaptability by sequencing and annotating 50 Claviceps genomes, representing 21 species, for a comprehensive comparison of genome architecture and plasticity in relation to host range potential. Our results show the trajectory from specialized genomes (sects. Citrinae and Paspalorum) toward adaptive genomes (sects. Pusillae and Claviceps) through colocalization of transposable elements around predicted effectors and a putative loss of repeat-induced point mutation resulting in unconstrained tandem gene duplication coinciding with increased host range potential and speciation. Alterations of genomic architecture and plasticity can substantially influence and shape the evolutionary trajectory of fungal pathogens and their adaptability. Furthermore, our study provides a large increase in available genomic resources to propel future studies of Claviceps in pharmacology and agricultural research, as well as, research into deeper understanding of the evolution of adaptable plant pathogens.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Deep Generative Modeling of Periodic Variable Stars Using Physical Parameters

The ability to generate physically plausible ensembles of variable sources is critical to the optimization of time domain survey cadences and the training of classification models on data sets with few to no labels. Traditional data augmentation techniques expand training sets by reenvisioning observed exemplars, seeking to simulate observations of specific training sources under different (exogenous) conditions. Unlike fully theory-driven models, these approaches do not typically allow principled interpolation nor extrapolation. Moreover, the principal drawback of theory-driven models lies in the prohibitive computational cost of simulating source observables from ab initio parameters. In this work, we propose a computationally tractable machine learning approach to generate realistic light curves of periodic variables capable of integrating physical parameters and variability classes as inputs. Our deep generative model, inspired by the transparent latent space generative adversarial networks, uses a variational autoencoder (VAE) architecture with temporal convolutional network layers, trained using the OGLE-III optical light curves and physical characteristics (e.g., effective temperature and absolute magnitude) from Gaia DR2. A test using the temperature–shape relationship of RR Lyrae demonstrates the efficacy of our generative “physics-enhanced latent space VAE” (PELS-VAE) model. Such deep generative models, serving as nonlinear nonparametric emulators, present a novel tool for astronomers to create synthetic time series over arbitrary cadences.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Mycorrhiza better predict soil fungal community composition and function than aboveground traits in temperate forest ecosystems

We sampled soil in the organic and mineral horizons beneath two AM-associated (Fraxinus americana, Thuja occidentalis) and two ECM-associated tree species (Betula alleghaniensis, and Tsuga canadensis), with an evergreen and deciduous species in each mycorrhizal group. To characterize fungal communities and organic matter decomposition beneath each tree species, we sequenced the ITS1 region of fungal DNA and measured the potential activity of carbon and nitrogen-targeting extracellular enzymes. Each tree species harbored distinct fungal communities, supporting the need to consider both mycorrhizal type and leaf habit. However, between tree characteristics, mycorrhizal type better predicted fungal communities. Across fungal guilds, saprotrophic fungi were the most important group in shaping fungal community differences in soils beneath all tree species. The effect of leaf habit on carbon and nitrogen-targeting hydrolytic enzymes depended on tree mycorrhizal association in the organic horizon, while oxidative enzyme activities were higher beneath EcM-associated trees across both soil horizons and leaf habits. These data include extracellular hydrolytic and oxidative enzyme activities, ITS sequencing fungal community data, soil carbon to nitrogen, and soil pH. Site level data include climate (mean annual temperature, precipitation), elevation, and soil series and order information.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Topography and functional traits shape the distribution of key shrub plant functional types in low-Arctic tundra

The expansion of shrubs in the Arctic tundra fundamentally modifies land-atmosphere interactions. However, it remains unclear how shrub distribution and expansion differ across key species due to challenges with discriminating tundra plant species at regional scales. Here, we combined multi-scale, multi-platform remote sensing and in situ trait measurements to elucidate the distribution patterns and primary controls of two representative deciduous-tall-shrub (DTS) genera, Alnus and Salix, in low-Arctic tundra. We show that topographic features were a key control on DTSs, creating heterogeneous, but predictable distributions of Alnus and Salix fractional cover (fCover). Alnus was more tolerant of elevation and slope and was found on hilly uplands (slope >10°) within a specific elevational band (200–400 m above sea level [MSL]). In contrast, Salix occurred at lower elevations (50–300 m MSL) on gentler slopes (3-10°) and required adequate soil moisture associated with its profligate water use. We also show that niche differentiation between Alnus and Salix changed with patch size, where larger patches were more specialized in resource requirements than individual plants of Alnus and Salix. To understand what constrains the growth of DTSs at locations with low fCover, we developed environmental limiting factor models, which showed that topography limits the upper bound of Alnus and Salix fCover in 69.2% and 48.7% of the landscape, respectively. These findings highlight a critical need to better understand and represent topography-controlled processes and functional traits in regulating shrub distribution, as well as a need for more detailed species classification to predict shrubification in the Arctic.

alder↗

DOE OSTI2019

Late Quaternary Tectonics, Incision, and Landscape Evolution of the Calchaquí River Catchment, Eastern Cordillera, NW Argentina

Unraveling the relative impacts of climate, tectonics, and lithology on landscape evolution is complicated by the temporal and spatial scale over which observations are made. We use soil and desert pavement classification, longitudinal river profiles, 10 Be–derived catchment mean modern and paleo–erosion rates, and vertical incision rates to test whether, if we restrict our analyses to a spatial scale over which climate is relatively invariant, tectonic and lithologic factors will dominate the late Quaternary landscape evolution of the Calchaquí River Catchment, NW Argentina. We find that the spatial distribution of erosion rates, normalized channel steepness indices, and concavity indices reflect active tectonics and lithologic resistance. Knickpoints are spatially coincident with tectonic and/or lithologic discontinuities, indicating local base–level control by faulting. Catchment mean erosion rates, ranging from 22.5 ± 2.6 to 121.9 ± 13.7 mm/kyr, and paleo–erosion rates, ranging from 56 +43 / –19 to 105 +60 / –33 mm/kyr, are similar, possibly suggesting that Quaternary climate changes have not had a strong enough influence on erosion rates to be detected using cosmogenic 10 Be. However, punctuated abandonment of pediment and strath terraces at 43.6 +15.0 / –11.6 , 91.2 +54.2 / –22.2 , and 151 +92.7 / –34.1 ka and disparities between vertical incision rates and catchment mean erosion rates could suggest periods of landscape transience, possibly reflecting climate cyclicity. Furthermore, our results emphasize the role of tectonic uplift and lithologic contrasts in shaping long–term erosion rates and channel morphology at the relatively local scale of the Calchaqui River Catchment, in contrast to regional–scale studies which find precipitation to exert the dominant control.

10Be catchment mean erosion rates↗

DOE OSTI2025

Novel CHI3L1 ‐Associated Angiogenic Phenotypes Define Glioma Microenvironments: Insights From Multi‐Omics Integration

ABSTRACT The CHI3L1 signaling pathway significantly influences glioma angiogenesis, but its role in the tumor microenvironment (TME) remains elusive. We propose a novelCHI3L1‐associated vascular phenotype classification for glioma through integrative analyses of multiple datasets with bulk and single‐cell transcriptome, genomics, digital pathology, and clinical data. We investigated the biological characteristics, genomic alterations, therapeutic vulnerabilities, and immune profiles within these phenotypes through a comprehensive multi‐omics approach. We constructed the vascular‐related risk (VR) score based onCHI3L1‐associated vascular signatures (CAVS) identified by machine learning algorithms. Utilizing unsupervised consensus clustering, gliomas were stratified into three distinct vascular phenotypes: Cluster A, marked by high vascularization and stromal activation with a relatively low levels of tumor‐infiltrating lymphocytes (TILs); Cluster B, characterized by moderate vascularization and stromal activity, coupled with a high density of TILs; and Cluster C, defined by low vascularization and sparse immune cell infiltration. We observed that the CAVS effectively indicated glioma‐associated angiogenesis and immune suppression by single‐cell RNA‐seq analysis. Moreover, the high‐VR‐score group exhibited enhanced angiogenic activity, reduced immune response, resistance to immunotherapy, and poorer clinical outcomes. The VR score independently predicted glioma prognosis and, combined with a nomogram, provided a robust clinical decision‐making tool. Potential drug prediction based on transcription factors for high‐risk patients was also performed. Our study reveals thatCHI3L1‐associated vascular phenotypes shape distinct immune landscapes in gliomas, offering insights for optimizing therapeutic strategies to improve patient outcomes.

Oncology↗

DOE OSTI2026

Investigating the Effects of Bars on Star Formation and Nuclear Activity of Galaxies Using DESI Survey Data

We present a statistical analysis of the connections between galactic bars, star formation, and active galactic nucleus (AGN) activity using 33,201 disk galaxies (0.01 < z < 0.05) from Dark Energy Spectroscopic Instrument Data Release 1 (DESI DR1) cross-matched with Galaxy Zoo DESI. Based on morphological classifications, we identify 3508 strongly barred and 8335 weakly barred systems. We find that barred galaxies exhibit a clear bimodal distribution in color–mass space: weak bars are preferentially found in bluer, lower-mass disks, whereas strong bars are more common in massive, redder systems. Strongly barred galaxies are on average more massive and metal rich than unbarred systems. In addition, strong bars enhance central star formation rates (SFRs) in low-mass galaxies but reduce specific SFRs in massive systems, reflecting a dual role where bars initially trigger central star formation but eventually promote quenching by accelerating gas consumption. In terms of nuclear activity, barred galaxies display a higher incidence of AGN activity. The presence of a bar is also associated with an increased fraction of powerful AGN, with the highest proportions found in strongly barred systems. However, the correlations between AGN activity and detailed bar structural parameters are weak, suggesting that the link between bars and nuclear activity is indirect and regulated by multiple factors. Overall, our results support a scenario in which bars facilitate angular momentum transport and gas inflow, thereby driving central star formation and fueling supermassive black hole accretion while operating alongside other processes that shape galaxy evolution.

Liu, Jianfei [Chinese Academy of Sciences (CAS), B↗

DOE OSTI2025

Knowledge Oriented Graph Unified Transformer (KOGUT) v0.1

KOGUT — Knowledge Oriented Graph Unified Transformer KOGUT implements the Relational Graph Transformer (RelGT) architecture for knowledge graph link prediction in biological domains, with a primary focus on microbial growth media prediction. While the original RelGT (arXiv:2505.10960) targets relational tables, time series, and multi-table databases, KOGUT adapts this architecture for heterogeneous biological knowledge graphs, providing first-in-class AI predictive models for microbial cultivation. Key Adaptations Beyond Original RelGT: - Knowledge Graph Focus: Applied to biological KGs with semantic node types (taxa, chemicals, media, phenotypes, environments) versus generic relational database tables, trained on the KG-Microbe knowledge graph (1.3M entities, 2.9M edges, 24 relation types). - Multimodal Node Encoding: Integrates node labels, categories, descriptions, and synonyms from KG metadata through learned embedding layers—adapting relational column features to graph node attributes with textual semantics. - Extended K-Hop Subgraph Strategy: Optimized neighborhood sampling (3-hop default, configurable up to 200 nodes) tuned for sparse biological networks, building on the original local-global attention framework with biological relation preservation. - Biolink Predicate Preservation: Type-specific transformations for 24 biological edge semantics (occurs_in, consumes, produces, has_phenotype, subclass_of) beyond standard relational foreign keys, enabling multi-relation link prediction. - Inductive Learning Support: Enables zero-shot predictions for novel taxa through feature-based embeddings (temperature, oxygen requirements, gram stain, cell shape), extending the original transductive relational benchmark scope to uncultured microorganisms. CheapSOTA Performance Optimizations (This Distribution): - VQ-EMA Centroid Attention: Vector quantization with exponential moving average for improved global context modeling (+5-10% MRR improvement). - HDF5 Precomputed Data Loading: One-time preprocessing of k-hop subgraphs to eliminate redundant graph traversals (2-5× training speedup). - Distributed Data Parallel Training: Multi-GPU support for scaling to larger knowledge graphs (tested on 4× NVIDIA A100 GPUs at NERSC Perlmutter). - Mixed Precision Training: Automatic mixed precision (AMP) for memory efficiency and faster training. Advantages Over Standard Knowledge Graph Embedding Models: Combines RelGT's proven multi-element tokenization (features, type, hop, structure) with graph-native biological representations, enabling interpretable link prediction across heterogeneous entities that standard embedding models (TransE, RotatE, ComplEx) and table-based transformers cannot directly model. Achieves near-perfect performance on microbial growth media prediction (MRR: 0.9966, Precision@1: 0.9932, Hit@10: 1.0000) while maintaining explainability through attention-based reasoning over biological pathways. Training Data: - KG-Microbe merged knowledge graph: 1,379,337 nodes, 2,960,472 edges - 24 biological relation types including taxonomic hierarchies, metabolic interactions, phenotype associations, and environmental relationships - Primary prediction task: Growth media suitability for microbial taxa (biolink:occurs_in, 50K edges) - Multi-relation capability: Predicts links for any of the 24 relation types, including chemical consumption/production, phenotype associations, and taxonomic classification Citation: Original RelGT Architecture: Dwivedi et al., "Relational Graph Transformer", arXiv:2505.10960, 2025 KOGUT Implementation: Knowledge Oriented Graph Unified Transformer for Microbial Growth Media Prediction Developed at Lawrence Berkeley National Laboratory (LBNL) Trained on NERSC Perlmutter supercomputer

Joachimiak, Marcin [Lawrence Berkeley National Lab↗

DOE OSTI2023

Evolution of neuronal cell classes and types in the vertebrate retina

The basic plan of the retina is conserved across vertebrates, yet species differ profoundly in their visual needs1. Retinal cell types may have evolved to accommodate these varied needs, but this has not been systematically studied. Here we generated and integrated single-cell transcriptomic atlases of the retina from 17 species: humans, two non-human primates, four rodents, three ungulates, opossum, ferret, tree shrew, a bird, a reptile, a teleost fish and a lamprey. We found high molecular conservation of the six retinal cell classes (photoreceptors, horizontal cells, bipolar cells, amacrine cells, retinal ganglion cells (RGCs) and Müller glia), with transcriptomic variation across species related to evolutionary distance. Major subclasses were also conserved, whereas variation among cell types within classes or subclasses was more pronounced. However, an integrative analysis revealed that numerous cell types are shared across species, based on conserved gene expression programmes that are likely to trace back to an early ancestral vertebrate. The degree of variation among cell types increased from the outer retina (photoreceptors) to the inner retina (RGCs), suggesting that evolution acts preferentially to shape the retinal output. Finally, we identified rodent orthologues of midget RGCs, which comprise more than 80% of RGCs in the human retina, subserve high-acuity vision, and were previously believed to be restricted to primates. By contrast, the mouse orthologues have large receptive fields and comprise around 2% of mouse RGCs. Projections of both primate and mouse orthologous types are overrepresented in the thalamus, which supplies the primary visual cortex. We suggest that midget RGCs are not primate innovations, but are descendants of evolutionarily ancient types that decreased in size and increased in number as primates evolved, thereby facilitating high visual acuity and increased cortical processing of visual information.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Modeling Nondestructive Defect Detection in Additively Manufactured Metallic Structures for Nuclear Applications

The future of quickly, economically produced metallic nuclear reactor parts with minimal supply-chain dependence lies in Laser Powder Bed Fusion (LPBF) Additive Manufacturing (AM): a 3D printing method involving laser melting and net shaping stainless steel and Inconel metallic powder into a solid structure. However, intrinsic features in LPBF frequently leads to the formation of materials defects, such as pores, within 3D printed structures. As safe long-term use in energy applications requires knowledge of all relevant defects before deployment in a reactor, we must develop methods for nondestructive detection of these defects. We are investigating Pulsed Thermal Tomography (PTT), which is a non-contact nondestructive imaging method scalable to arbitrary structure size. Thermal tomography (TT) is a computational method for 3D spatial reconstruction of material thermal effusivity from flash or pulsed thermography temperature data cube. Thermography data cube consists of 2D surface temperature measurements at different times. The objective of the present work is to investigate limits on defect detection in AM metallic structures with PTT. To this effect, we modeled PTT with COMSOL heat transfer computer simulations. We developed a layered media COMSOL simulation consisting of a Stainless Steel 316 (SS316) plate with an internal layer of un-sintered SS316 powder. Thermophysical properties of the powder layer were modeled with equivalent volume mixing model. To account for partial sintering at the boundary of the defect, the transition between solid and powder layers was modeled as a Gaussian. Using data from COMSOL simulations, we reconstructed depth-dependent thermal effusivity, which allowed defect visibility estimation. A series of parametric studies determined that at 1mm depth, 50µm is the smallest detectable defect. In addition, classification of the defects which can lead to early fatigue of the metallic structure in a reactor is briefly discussed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Performance of Pulsed Thermal Tomography Imaging with Machine Learning-Based Classification of Defects in Additively Manufactured Structures

Additive manufacturing (AM) is an emerging method for cost-efficient fabrication of complex topology nuclear reactor parts from high-strength corrosion resistance alloys, such as stainless steel and Inconel. AM of metallic structures for nuclear energy applications is currently based on laser powder bed fusion (LPBF) process. Some of the challenges with using LPBF method for nuclear manufacturing include the possibility of introducing pores into metallic structures. Integrity of AM structures needs to be evaluated nondestructively because material flaws could lead to premature failures in high temperature nuclear reactor environment. Currently, there exist limited capabilities to evaluate actual AM structures non-destructively. Pulsed Thermal Tomography Imaging (PTT) provides a capability for non-destructive evaluation (NDE) of subsurface defects in arbitrary size structures. The PTT method is based on recording material surface temperature transients with infrared (IR) camera following thermal pulse delivered on material surface with flash light. The PTT method has advantages for NDE of actual AM structures because the method involves one-sided non-contact measurements and fast processing of large sample areas captured in one image. Following initial qualification of an AM component for deployment in a nuclear reactor, a PTT system can also be used for in-service nondestructive evaluation (NDE) applications. In this report, we describe recent progress in enhancing PTT capabilities in detecting and visualizing microscopic defects in metallic specimens. The thermal tomography (TT) algorithm obtains depth reconstructions of spatial effusivity from the data cube of sequentially recorded surface temperatures. However, interpretation of TT images is non-trivial because of blurring of images with increasing depth. To address this challenge, we have developed a deep learning convolutional neural network (CNN) to classify size and orientation subsurface defects in simulated TT images. CNN is trained on a database of TT images created for a set of simulated metallic structures with elliptical subsurface voids. Test of CNN performance demonstrate the ability to classify radii and angular orientation of subsurface defects in TT images. In addition, we have shown that CNN trained on elliptical defects is capable of classifying irregular-shaped defects obtained from scanning electron microscopy (SEM) of stainless steel sections printed with LPBF.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2017

Recent publications from the Alzheimer's Disease Neuroimaging Initiative: Reviewing progress toward improved AD clinical trials

Introduction The Alzheimer's Disease Neuroimaging Initiative (ADNI) has continued development and standardization of methodologies for biomarkers and has provided an increased depth and breadth of data available to qualified researchers. This review summarizes the over 400 publications using ADNI data during 2014 and 2015. Methods We used standard searches to find publications using ADNI data. Results (1) Structural and functional changes, including subtle changes to hippocampal shape and texture, atrophy in areas outside of hippocampus, and disruption to functional networks, are detectable in presymptomatic subjects before hippocampal atrophy; (2) In subjects with abnormal β‐amyloid deposition (Aβ+), biomarkers become abnormal in the order predicted by the amyloid cascade hypothesis; (3) Cognitive decline is more closely linked to tau than Aβ deposition; (4) Cerebrovascular risk factors may interact with Aβ to increase white‐matter (WM) abnormalities which may accelerate Alzheimer's disease (AD) progression in conjunction with tau abnormalities; (5) Different patterns of atrophy are associated with impairment of memory and executive function and may underlie psychiatric symptoms; (6) Structural, functional, and metabolic network connectivities are disrupted as AD progresses. Models of prion‐like spreading of Aβ pathology along WM tracts predict known patterns of cortical Aβ deposition and declines in glucose metabolism; (7) New AD risk and protective gene loci have been identified using biologically informed approaches; (8) Cognitively normal and mild cognitive impairment (MCI) subjects are heterogeneous and include groups typified not only by “classic” AD pathology but also by normal biomarkers, accelerated decline, and suspected non‐Alzheimer's pathology; (9) Selection of subjects at risk of imminent decline on the basis of one or more pathologies improves the power of clinical trials; (10) Sensitivity of cognitive outcome measures to early changes in cognition has been improved and surrogate outcome measures using longitudinal structural magnetic resonance imaging may further reduce clinical trial cost and duration; (11) Advances in machine learning techniques such as neural networks have improved diagnostic and prognostic accuracy especially in challenges involving MCI subjects; and (12) Network connectivity measures and genetic variants show promise in multimodal classification and some classifiers using single modalities are rivaling multimodal classifiers. Discussion Taken together, these studies fundamentally deepen our understanding of AD progression and its underlying genetic basis, which in turn informs and improves clinical trial design.

Weiner, Michael W.↗

DOE OSTI2024

Data and scripts associated with a manuscript on a meta-analysis synthesizing stream biogeochemical response to wildfires across space and time (v2)

This data package is associated with the publication “Catchment characteristics modulate the influence of wildfires on nitrate and dissolved organic carbon in lotic systems across space and time: A meta-analysis” submitted to Global Biogeochemical Cycles (Cavaiani et al. 2025). This study uses meta-analytical techniques to evaluate the effect of wildfire on in-stream responses in burned and unburned watersheds. The study aims to provide additional insight into the range of responses and net influences that wildfires have on hydro-biogeochemistry across broad spatial scales, burn extents, and the persistence of water-quality change. This study compiles data and metadata from 18 total publications that includes 1) surface water geochemistry data (dissolved organic carbon; nitrate), 2) climate classifications, 3) year of the wildfire, 4) the time lag between when the fire occurred and when the sampling occurred, and 5) study design of the publication. In total, this meta-analysis draws data that spans 8 climate guilds, 3 biomes, 62 watersheds, and 20 unique wildfires. See Sites_meta_data.csv for citations of the papers used in this meta-analysis. All R scripts and the associated data can also be found on GitHub at This data package was originally published in March 2024. It was updated in April 2025 (v2; new and modified files). See the change history section in the readme for more details. This data package contains five primary folders that include the following: (1) inputs; (2) output for analysis; (3) initial plots; (4) R scripts; and (5) GIS data. The data package also contains a data dictionary (dd) that provides column header definitions and a file-level metadata (flmd) file that describes every file. The “inputs” folder contains a list of all publications identified during the formal web search and an indication of whether each publication was included in the final analysis. Additionally, it includes site-level metadata, catchment characteristics, and GIS data for all publications included in the final analysis. The “Output_for_analysis” folder contains all data frames and figures generated from each R script used for additional data analysis. The “initial_plots” folder includes all exploratory figures that will be included in a supplemental and figures that will be submitted with the manuscript for publication. The “R_scripts” folder contains the scripts that perform all the data manipulations, statistical analyses, and plots. The “gis_data” folder includes shape files for each fire included in this meta-analysis. This data package contains the following file types: csv, pdf, jpeg, cpg, dbf, prj, shp, shp.ea.iso.xml, shp.iso.xml, shx.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Reduction of detection limit and quantification uncertainty due to interferent by neural classification with abstention

Many measurements in the physical sciences can be cast as counting experiments, where the number of occurrences of a physical phenomenon informs the prevalence of the phenomenon's source. Often, detection of the physical phenomenon (termed signal) is difficult to distinguish from naturally occurring phenomena (termed background). In this case, the discrimination of signal events from background can be performed using classifiers, and they may range from simple, threshold-based classifiers to sophisticated neural networks. These classifiers are often trained and validated to obtain optimal accuracy, however we show that the optimal accuracy classifier does not generally coincide with a classifier that provides the lowest detection limit, nor the lowest quantification uncertainty. Here, we present a derivation of the detection limit and quantification uncertainty in the classifier-based counting experiment case. We also present a novel abstention mechanism to minimize the detection limit or quantification uncertainty a posteriori. We illustrate the method on two data sets from the physical sciences, discriminating Ar-37 and Ar-39 radioactive decay from non-radioactive events in a gas proportional counter, and discriminating neutrons from gammas in an inorganic scintillator and report results therefrom.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2016

Three-dimensional reconstruction of the size and shape of protein microcrystals using Bragg coherent diffractive imaging

Three-dimensional imaging of protein crystals during X-ray diffraction experiments opens up a range of possibilities for optimising crystal quality and gaining new insights into the fundamental processes that drive radiation damage. Obtaining this information at the appropriate lengthscales however is extremely challenging. One approach that has been recently demonstrated as a promising avenue for charactering the size and shape of protein crystals at nanometre lengthscales is Bragg Coherent Diffractive Imaging (BCDI). BCDI is a recently developed technique that is able to recover the phase of the continuous diffraction intensity signal around individual Bragg peaks. When data is collected at multiple points on a rocking curve a Reciprocal Space Map (RSM) can be assembled and then inverted using BCDI to obtain a three-dimensional image of the crystal. The first demonstration of two-dimensional BCDI of protein crystals was reported by Boutet at al., recently this work was extended to the study of radiation damage of micron-sized crystals. Here we present the first three-dimensional reconstructions of a Lysozyme protein crystal using BDI. The results are validated against RSM and TEM data and have implications for both radiation damage studies and for developing new approaches to structure retrieval from micron-sized protein crystals.

Bragg coherent diffractive imaging↗

DOE OSTI2015

GOTTCHA Database, Version 1

One major challenge in the field of shotgun metagenomics is the accurate identification of the organisms present within the community, based on classification of short sequence reads. Though microbial community profiling methods have emerged to attempt to rapidly classify the millions of reads output from contemporary sequencers, the combination of incomplete databases, similarity among otherwise divergent genomes, and the large volumes of sequencing data required for metagenome sequencing has led to unacceptably high false discovery rates (FDR). Here we present the application of a novel, gene-independent and signature-based metagenomic taxonomic profiling tool with significantly smaller FDR, which is also capable of classifying never-before seen genomes into the appropriate parent taxa.The algorithm is based upon three primary computational phases: (I) genomic decomposition into bit vectors, (II) bit vector intersections to identify shared regions, and (III) bit vector subtractions to remove shared regions and reveal unique, signature regions.In the Decomposition phase, genomic data is first masked to highlight only the valid (non-ambiguous) regions and then decomposed into overlapping 24-mers. The k-mers are sorted along with their start positions, de-replicated, and then prefixed, to minimize data duplication. The prefixes are indexed and an identical data structure is created for the start positions to mimic that of the k-mer data structure.During the Intersection phase -- which is the most computationally intensive phase -- as an all-vs-all comparison is made, the number of comparisons is first reduced by four methods: (a) Prefix restriction, (b) Overlap detection, (c) Overlap restriction, and (d) Result recording. In Prefix restriction, only k-mers of the same prefix are compared. Within that group, potential overlap of k-mer suffixes that would result in a non-empty set intersection are screened for. If such an overlap exists, the region which intersects is first reduced by performing a binary search of the boundary suffixes of the smaller set into the larger set, which defines the limits of the zipper-based intersection process. Rather than recording the actual k-mers of the intersection, another data structure of identical "shape" is created which consists of only bit vectors so that only a 1 or 0 will be stored in the location of the k-mer suffix that was found in the intersection. This reduces the amount of data generated and stored considerably.During the Subtraction phase, relevant intersection bitmasks are first unionized together to form a single bitmask which is then applied over the original genome to reveal only those regions of the genome that are unique. These regions are then exported to disk in FASTA format and used in the application of determining the constituents of an unknown metagenomic community.The DATABASE provided is the result of the algorithm described.

Freitas, Tracey↗

DOE OSTI2015

GOTTCHA Database, Version 1

Freitas, Tracey↗

DOE OSTI2015

GOTTCHA, Version2

One major challenge in the field of shotgun metagenomics is the accurate identification of the organisms present within the community, based on classification of short sequence reads. Though microbial community profiling methods have emerged to attempt to rapidly classify the millions of reads output from contemporary sequencers, the combination of incomplete databases, similarity among otherwise divergent genomes, and the large volumes of sequencing data required for metagenome sequencing has led to unacceptably high false discovery rates (FDR). Here we present the application of a novel, gene-independent and signature-based metagenomic taxonomic profiling tool with significantly smaller FDR, which is also capable of classifying never-before seen genomes into the appropriate parent taxa.The algorithm is based upon three primary computational phases: (I) genomic decomposition into bit vectors, (II) bit vector intersections to identify shared regions, and (III) bit vector subtractions to remove shared regions and reveal unique, signature regions. In the Decomposition phase, genomic data is first masked to highlight only the valid (non-ambiguous) regions and then decomposed into overlapping 24-mers. The k-mers are sorted along with their start positions, de-replicated, and then prefixed, to minimize data duplication. The prefixes are indexed and an identical data structure is created for the start positions to mimic that of the k-mer data structure. During the Intersection phase -- which is the most computationally intensive phase -- as an all-vs-all comparison is made, the number of comparisons is first reduced by four methods: (a) Prefix restriction, (b) Overlap detection, (c) Overlap restriction, and (d) Result recording. In Prefix restriction, only k-mers of the same prefix are compared. Within that group, potential overlap of k-mer suffixes that would result in a non-empty set intersection are screened for. If such an overlap exists, the region which intersects is first reduced by performing a binary search of the boundary suffixes of the smaller set into the larger set, which defines the limits of the zipper-based intersection process. Rather than recording the actual k-mers of the intersection, another data structure of identical "shape" is created which consists of only bit vectors so that only a 1 or 0 will be stored in the location of the k-mer suffix that was found in the intersection. This reduces the amount of data generated and stored considerably. During the Subtraction phase, relevant intersection bitmasks are first unionized together to form a single bitmask which is then applied over the original genome to reveal only those regions of the genome that are unique. These regions are then exported to disk in FASTA format and used in the application of determining the constituents of an unknown metagenomic community.

Freitas, Tracey↗