Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Functional Data Analysis”

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

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

At least 631 records · Page 35

Chemoproteogenomic stratification of the missense variant cysteinome

Abstract Cancer genomes are rife with genetic variants; one key outcome of this variation is widespread gain-of-cysteine mutations. These acquired cysteines can be both driver mutations and sites targeted by precision therapies. However, despite their ubiquity, nearly all acquired cysteines remain unidentified via chemoproteomics; identification is a critical step to enable functional analysis, including assessment of potential druggability and susceptibility to oxidation. Here, we pair cysteine chemoproteomics—a technique that enables proteome-wide pinpointing of functional, redox sensitive, and potentially druggable residues—with genomics to reveal the hidden landscape of cysteine genetic variation. Our chemoproteogenomics platform integrates chemoproteomic, whole exome, and RNA-seq data, with a customized two-stage false discovery rate (FDR) error controlled proteomic search, which is further enhanced with a user-friendly FragPipe interface. Chemoproteogenomics analysis reveals that cysteine acquisition is a ubiquitous feature of both healthy and cancer genomes that is further elevated in the context of decreased DNA repair. Reference cysteines proximal to missense variants are also found to be pervasive, supporting heretofore untapped opportunities for variant-specific chemical probe development campaigns. As chemoproteogenomics is further distinguished by sample-matched combinatorial variant databases and is compatible with redox proteomics and small molecule screening, we expect widespread utility in guiding proteoform-specific biology and therapeutic discovery.

Desai, Heta (ORCID:0000000343621707)↗

Population‐level gene expression can repeatedly link genes to functions in maize

SUMMARY Transcriptome‐wide association studies (TWAS) can provide single gene resolution for candidate genes in plants, complementing genome‐wide association studies (GWAS) but efforts in plants have been met with, at best, mixed success. We generated expression data from 693 maize genotypes, measured in a common field experiment, sampled over a 2‐h period to minimize diurnal and environmental effects, using full‐length RNA‐seq to maximize the accurate estimation of transcript abundance. TWAS could identify roughly 10 times as many genes likely to play a role in flowering time regulation as GWAS conducted data from the same experiment. TWAS using mature leaf tissue identified known true‐positive flowering time genes known to act in the shoot apical meristem, and trait data from a new environment enabled the identification of additional flowering time genes without the need for new expression data. eQTL analysis of TWAS‐tagged genes identified at least one additional known maize flowering time gene through trans ‐eQTL interactions. Collectively these results suggest the gene expression resource described here can link genes to functions across different plant phenotypes expressed in a range of tissues and scored in different experiments.

Torres‐Rodríguez, J. Vladimir↗

Exploring physics of ferroelectric domain walls via Bayesian analysis of atomically resolved STEM data

The physics of ferroelectric domain walls is explored using the Bayesian inference analysis of atomically resolved STEM data. We demonstrate that domain wall profile shapes are ultimately sensitive to the nature of the order parameter in the material, including the functional form of Ginzburg-Landau-Devonshire expansion, and numerical value of the corresponding parameters. The preexisting materials knowledge naturally folds in the Bayesian framework in the form of prior distributions, with the different order parameters forming competing (or hierarchical) models. Here, we explore the physics of the ferroelectric domain walls in BiFeO 3 using this method, and derive the posterior estimates of relevant parameters. More generally, this inference approach both allows learning materials physics from experimental data with associated uncertainty quantification, and establishing guidelines for instrumental development answering questions on what resolution and information limits are necessary for reliable observation of specific physical mechanisms of interest.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Enhanced Signal Processing of Distributed Brillouin Fiber Sensors using a Decoupled Radial Basis Function Network

A novel decoupled radial basis function network (D-RBFN) is proposed to accelerate signal processing and address the big data challenges associated with ultra-long distance Brillouin optical time-domain analysis (BOTDA) systems. The proposed frame- work is demonstrated on a dataset measured over a 100 km distance using a bi-directional Raman assisted BOTDA system.

Venketeswaran, Abhishek↗

19-LW-045 Full Length Final Report. Molecular Mechanisms of Bacterial Pathogenesis: Waging the Arms Race with Superbugs

As the current global pandemic makes abundantly clear, we need a better understanding of infectious disease to safeguard human health, the economy and global security. Modern omics techniques hold the promise of providing a comprehensive understanding of the molecular mechanisms of life, including causes of pathogenesis from infectious disease at the molecular level, but we there is a serious gap in annotation of gene function. For as much as half of the genes and gene products encoded in genomes the molecular and/or cellular function is unknown or only partially understood. Recent innovations in fluorescence microscopy for live cell imaging and genetic engineering make it possible to determine the temporal correlation between molecular events, such as a gene being expressed due to host-pathogen interaction, and cellular events, such as bacterial invasion of immune cells. This is turn allows us to gain new insight as to the molecular and cellular role of individual genes and will enable the discovery and validation of new molecular mechanisms essential for infectious disease. Knowing the molecular mechanisms of disease processes will provide new therapeutic targets or novel countermeasure strategies. We aimed to develop a lattice light sheet fluorescence microscope as a unique resource at LLNL for long time course live cell imaging experiments; to develop the reagents and cell lines needed to monitor molecular events during the course pathogenic bacteria infecting mammalian immune cells; and to demonstrate that we could capture molecular events during an infection. We fully commissioned the LLNL lattice light sheet microscope and conducted initial proof of principle imaging experiments on mammalian immune cells and pathogenic bacteria. It is clear from the experience gained that long time course live cell imaging has tremendous potential to help elucidate molecular mechanisms of host-pathogen interactions and to help annotate gene function, which would establish a basis for new countermeasures. It is also clear that if live cell imaging is to realize its full potential new data processing and analysis tools will need to be developed to facilitate analysis of molecular events within cells; new sample chambers and stages could facilitate studies with a wider range of cell and tissue types; and alternative molecular tagging methods need to be explored to enable more facile engineering of cells labeled with molecular specificity.

59 BASIC BIOLOGICAL SCIENCES↗

Coupled Bunch Instabilities Growth in the Fermilab Booster during Acceleration Cycle

Currently, Fermilab Booster accelerates ~4.5E12 protons per pulse (ppp) in 81 bunches from 400 MeV to 8 GeV at 15 Hz to provide beam to multiple HEP experiments and is being upgraded to handle higher beam intensity >6.7E12 ppp at a repetition rate of 20 Hz. In the current mode of operation, we control longitudinal coupled bunch instabilities (LCBI) using combination of passive and active dampers. The issues with LCBI are expected to worsen at higher beam intensities. Here, we would like to investigate at what time in the beam cycle a particular mode is going to originate and how much it contributes at a different time of the acceleration cycle using the collected wall current monitor data from injection to extraction. Also, measure growth rates as a function of beam intensity for some dominant modes. Finally, this paper presents results from the observed LCBI analysis for the highest available beam intensities, unexpected findings, and future plans to mitigate the LCBI at higher beam intensities

43 PARTICLE ACCELERATORS↗

A Science Gateway for the Repeatable Analysis of Machine Learning Predicted Gravity Anomalies

In recent years, deep learning has become an increasingly popular alternative for modeling in geoscience applications due to its scalability and efficiency. However, the interpretability, compute, data volume, and hyperparameter tuning requirements of deep learning models make development and monitoring difficult. Furthermore, model explainability and communicating results obtained by these models to users or domain experts is a challenge, as domain experts in geoscience also need to have a deep understanding of how those models function in order to support their scientific works. Here, we describe a science gateway and machine learning pipeline for predicting gravity anomalies from geophysical data. The gateway, built on open-source technologies, provides a holistic view of the pipeline through interactive visualizations aimed at enabling efficient exploratory data analysis. The repeatability, reproducibility, and monitoring capabilities of this overall system allow us to iterate and analyze at scale. Using this pipeline and gateway, we can repeatedly produce accurate high-resolution gravity anomaly datasets. By describing the underlying technologies, implementation, and results, here we provide a foundation for the broader adoption of science gateways into cross-cutting geoscience and machine learning research projects as a means to improve the scientific discovery and collaboration in the geophysics and computational sciences community.

58 GEOSCIENCES↗

On the interaction of a wind turbine wake with a conventionally neutral atmospheric boundary layer

In this work, we investigate the dynamics of wind turbine tip-vortex breakdown in a conventionally neutral atmospheric boundary layer (ABL). To this end, high-resolution data are collected from large-eddy simulations of a wind turbine operating within a neutral ABL and studied by means of proper orthogonal decomposition (POD) and Fourier analysis. The high resolution of the generated data in both space and time allows us to gain insight into the tip-vortex breakdown mechanisms by (i) capturing the energy modes of the coherent structures, (ii) studying their contribution to the tip-vortex breakdown through their power spectra functions and mean kinetic energy (MKE) flux, and (iii) analysing the growth rate of each contributing perturbation frequency along tip vortices. Our analysis shows that under a fully turbulent scenario, the growth rate of perturbations along the tip vortices is largest for low wave numbers, i.e. long-wave perturbations. Additionally, the MKE flux reaches its highest value at two diameters downstream of the rotor plane, a behaviour that can be attributed to the coexistence of multiple interacting POD modes, with the streamwise vortex roller mode being the primary contributor to the total MKE flux budget, contributing approximately . Finally, comparisons with a laminar, uniform flow scenario subject to a single-frequency perturbation highlight the differences between the two ambient flow conditions. In the non-turbulent, uniform flow scenario, the growth rate attains its maximum value at a wave number corresponding to the out-of-phase mutual-inductance mechanism, whereas the MKE flux exhibits local minima and maxima along the wake and at different downstream locations depending on the perturbation frequency. Our analyses suggest that the breakdown of the wind turbine tip vortices under a fully turbulent neutral ABL inflow is due to complex interactions across a range of excitation frequencies, in which the mutual-inductance instability may not be the dominant one.

17 WIND ENERGY↗

BioXTAS RAW 2 : new developments for a free open-source program for small-angle scattering data reduction and analysis

BioXTAS RAW is a free open-source program for reduction, analysis and modelling of biological small-angle scattering data. Here, the new developments in RAW version 2 are described. These include improved data reduction using pyFAI ; updated automated Guinier fitting and D max finding algorithms; automated series ( e.g. size-exclusion chromatography coupled small-angle X-ray scattering or SEC-SAXS) buffer- and sample-region finding algorithms; linear and integral baseline correction for series; deconvolution of series data using regularized alternating least squares ( REGALS ); creation of electron-density reconstructions using electron density via solution scattering ( DENSS ); a comparison window showing residuals, ratios and statistical comparisons between profiles; and generation of PDF reports with summary plots and tables for all analysis. Furthermore, there is now a RAW API, which can be used without the graphical user interface (GUI), providing full access to all of the functionality found in the GUI. In addition to these new capabilities, RAW has undergone significant technical updates, such as adding Python 3 compatibility, and has entirely new documentation available both online and in the program.

97 MATHEMATICS AND COMPUTING↗

RWRtoolkit: multi-omic network analysis using random walks on multiplex networks in any species

Abstract We introduce RWRtoolkit, a multiplex generation, exploration, and statistical package built for R and command-line users. RWRtoolkit enables the efficient exploration of large and highly complex biological networks generated from custom experimental data and/or from publicly available datasets, and is species agnostic. A range of functions can be used to find topological distances between biological entities, determine relationships within sets of interest, search for topological context around sets of interest, and statistically evaluate the strength of relationships within and between sets. The command-line interface is designed for parallelization on high-performance cluster systems, which enables high-throughput analysis such as permutation testing. Several tools in the package have also been made available for use in reproducible workflows via the KBase web application.

Kainer, David (ORCID:0000000172714676)↗

Audacity of huge: overcoming challenges of data scarcity and data quality for machine learning in computational materials discovery

Machine learning (ML)-accelerated discovery requires large amounts of high-fidelity data to reveal predictive structure–property relationships. For many properties of interest in materials discovery, the challenging nature and high cost of data generation has resulted in a data landscape that is both scarcely populated and of dubious quality. Data-driven techniques starting to overcome these limitations include the use of consensus across functionals in density functional theory, the development of new functionals or accelerated electronic structure theories, and the detection of where computationally demanding methods are most necessary. When properties cannot be reliably simulated, large experimental data sets can be used to train ML models. In the absence of manual curation, increasingly sophisticated natural language processing and automated image analysis are making it possible to learn structure–property relationships from the literature. Finally, models trained on these data sets will improve as they incorporate community feedback.

36 MATERIALS SCIENCE↗

The landscape of regulatory element evolution in a C4 perennial grass

Gene regulatory evolution is a well-known source of phenotypic diversity and adaptive evolution. Although cis-regulatory elements (CREs) play a vital role in gene expression evolution, the molecular evolution of CREs remains mostly unknown due to the difficulty in identifying and characterizing these functional elements. Comparative genomic analyses of noncoding DNA can be leveraged to identify conserved noncoding sequences (CNS), many of which may harbor functional CREs conserved by purifying selection. However, purely computational inference of CREs from putative CNS can be erroneous due to the complex genomic architecture in plants. One promising experimental approach to identify CREs is by profiling accessible chromatin regions (ACRs) that are often associated with the location of CREs. In this study, we use comparative genomics along with the profiling of ACRs to study the molecular evolution of putative functional noncoding regulatory regions in Panicoid grasses. We identified sets of CNS that varied in relationship to the degree of evolutionary divergence among the studied taxa, including identifying core-Panicoid-CNS. We augmented this analysis by profiling ACRs in Panicum hallii ecotypes using ATAC-seq. ACRs had low SNP density at the summit, harbored a high frequency of core-Panicoid-CNS, and were enriched with expression QTL. These data help to annotate the P. hallii genome for putative functional elements and suggest that a large proportion of these ACRs are evolving under purifying selection. Turnover in CNS and ACR between ecotypes of P. hallii identifies a small set of putatively divergent CREs that may underlie differences in gene regulation between genotypes from inland and coastal habitats. In summary, we profiled ACRs in Panicoid grasses and integrated this data with our putative CNS prediction framework, which provides unique insight into patterns of polymorphism and divergence in CREs in C4 perennial grasses.

59 BASIC BIOLOGICAL SCIENCES↗

Revisiting the definition of field capacity as a functional parameter in a layered agronomic soil profile beneath irrigated maize

The soil water content at the condition of field capacity (θ FC ) is a key parameter in irrigation scheduling and has been suggested to be determined by running a synthetic drainage experiment until the flux rate (q) at the bottom of the soil profile achieves a predefined negligible value (q FC ). We question the impact of q FC on the assessment of field capacity. Moreover, calculating θ FC as the integral mean of the water content profile when q is equal to q FC is strictly valid only for uniform soil profiles. By contrast, this practice is ambiguous and biased for stratified soil profiles due to the soil water content discontinuity at the layer interfaces. In this study, the concept of field capacity was revisited and adapted to practical agronomic heuristics. By resorting to the assessment of root-zone water storage capacity (W), we envision field capacity as a functional hydraulic parameter derived from synthetic irrigation scheduling scenarios to minimize drought stress, drainage, and nitrate leachate below the root zone. A functional analysis was carried out on a 135-cm-thick layered soil profile beneath maize in eastern Nebraska. On-farm irrigation scheduling applications and agricultural practices were recorded for 20 years (2001–2020) at a daily time step. Hydrus-1D was calibrated and validated with direct measurements of the soil water retention curve and soil water content data, respectively, in each soil layer. A set of functional field capacity values was derived from 24 irrigation scheduling scenarios, and the optimal water storage capacity at field capacity (W FC ) was approximately 50 cm (corresponding to about 80% saturation in the soil profile). An average irrigation amount of 217.5 mm distributed over 21 events was obtained by using optimal irrigation scheduling, which was initiated when the matric pressure head took on a value of –700 cm and the irrigation rate was set at 1.0 cm d –1 . This irrigation practice ensured water storage at approximately the same level (ideally at W FC ) by sustaining only evapotranspiration fluxes in the uppermost portion of the root zone and by limiting excessive drainage. This protocol can be transferred to other agricultural fields.

54 ENVIRONMENTAL SCIENCES↗

Angular clustering properties of the DESI QSO target selection using DR9 Legacy Imaging Surveys

ABSTRACT The quasar target selection for the upcoming survey of the Dark Energy Spectroscopic Instrument (DESI) will be fixed for the next 5 yr. The aim of this work is to validate the quasar selection by studying the impact of imaging systematics as well as stellar and galactic contaminants, and to develop a procedure to mitigate them. Density fluctuations of quasar targets are found to be related to photometric properties such as seeing and depth of the Data Release 9 of the DESI Legacy Imaging Surveys. To model this complex relation, we explore machine learning algorithms (random forest and multilayer perceptron) as an alternative to the standard linear regression. Splitting the footprint of the Legacy Imaging Surveys into three regions according to photometric properties, we perform an independent analysis in each region, validating our method using extended Baryon Oscillation Spectroscopic Survey (eBOSS) EZ-mocks. The mitigation procedure is tested by comparing the angular correlation of the corrected target selection on each photometric region to the angular correlation function obtained using quasars from the Sloan Digital Sky Survey (SDSS) Data Release 16. With our procedure, we recover a similar level of correlation between DESI quasar targets and SDSS quasars in two-thirds of the total footprint and we show that the excess of correlation in the remaining area is due to a stellar contamination that should be removed with DESI spectroscopic data. We derive the Limber parameters in our three imaging regions and compare them to previous measurements from SDSS and the 2dF QSO Redshift Survey.

79 ASTRONOMY AND ASTROPHYSICS↗

Luminosity determination using Z boson production at the CMS experiment

The measurement of Z boson production is presented as a method to determine the integrated luminosity of CMS data sets. The analysis uses proton–proton collision data, recorded by the CMS experiment at the CERN LHC in 2017 at a center-of-mass energy of 13 TeV. Events with Z bosons decaying into a pair of muons are selected. The total number of Z bosons produced in a fiducial volume is determined, together with the identification efficiencies and correlations from the same data set, in small intervals of 20 pb –1 of integrated luminosity, thus facilitating the efficiency and rate measurement as a function of time and instantaneous luminosity. Using the ratio of the efficiency-corrected numbers of Z bosons, the precisely measured integrated luminosity of one data set is used to determine the luminosity of another. For the first time, a full quantitative uncertainty analysis of the use of Z bosons for the integrated luminosity measurement is performed. The uncertainty in the extrapolation between two data sets, recorded in 2017 at low and high instantaneous luminosity, is less than 0.5%. We show that the Z boson rate measurement constitutes a precise method, complementary to traditional methods, with the potential to improve the measurement of the integrated luminosity.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Progress Towards A Measurement of Neutrino Induced Charged Current Neutral Pion Production in the MicroBooNE Experiment

An analysis of MicroBooNE data with a signal of one muon, one neutral pion, and no charged pions is presented. Studying neutral pion production in the MicroBooNE detector provides an opportunity to better understand neutrino-argon interactions, and is crucial for future accelerator-based neutrino oscillation experiments. This analysis presents the progress towards the first measurement of the differential cross section for charged current (CC) $π^0$ production in neutrino-argon interactions. Using a dataset corresponding to about 7 × 10 20 protons on target (POT), we present an analysis which aims to measure the single differential cross sections as a function of the $π^0$ kinematic variables such as momentum and scattering angle. The Wiener-SVD technique for unfolding the measurement is presented and demonstrated using multiple generator predictions. A future iteration of this analysis will compare an unfolded data measurement to these models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Global gross primary production simulated by three terrestrial biosphere models to analyse smoothing of limiting rates

Processed global terrestrial gross primary production (GPP) from three terrestrial biosphere models (TBMs). Two simulations were run with each model, one that selects the rate limiting process of photosynthesis based on the minimum of the potential rates (following Farquhar et al. 1980) and the other that uses quadratic smoothing between the potential rates (following Collatz et al. 1991). The three TBMs are: 1) the Energy Exascale Earth System Model (E3SM) land model (ELM) (release: v1.1.0) (Burrows et al., 2020), a coupled carbon, nitrogen, and phosphorus model with sun/shade big-leaf canopy photosynthesis scaling. 2) The Functionally Assembled Terrestrial Ecosystem Simulator (FATES, tag: sci.1.30.0_api.8.0.0) (Koven et al., 2019), coupled with the Community Land Model (CLM, version 5; Lawrence et al., 2019), a carbon-only vegetation demography model with multi-leaf and multi-canopy layers for scaling photosynthesis, and with a leaf area index (LAI) optimisation scheme. 3) The Sheffield Dynamic Global Vegetation Model (SDGVM, tag: Walkeretal2020_GCB), a carbon-only model with multi-leaf-layer canopy photosynthesis scaling and also with an LAI optimisation scheme (Walker et al., 2017; Woodward & Lomas, 2004). These data were generated to evaluate the effect of Collatz smoothing of limiting rates on simulated GPP across a range of terrestrial biosphere models. Analysis of results can be found in Walker et al. (2020). This dataset is comprised of R scripts to process data and produce figures and .RDS files, a binary format specific to R and read with the 'readRDS' function.

54 ENVIRONMENTAL SCIENCES↗

Constructing Data-Driven Predictions at the Far Detector for NOvA's Neutrino Oscillation Analysis.

NOvA, is a two-detector, long-baseline neutrino oscillation experiment located at Fermilab, Batavia, IL, USA. It is designed primarily to constrain neutrino oscillation parameters using $\nu_\mu \ (\bar{\nu}_\mu)$ disappearance and $\nu_e \ (\bar{\nu}_e)$ appearance data. The Neutrinos at Main Injector (NuMI) beamline at Fermilab provides a high purity 900 KW intense beam of neutrinos and anti-neutrinos to NOvA. The NOvA Near Detector, located 100m underground and 1km away from the beam source, observes the un-oscillated $\nu_\mu \ (\bar{\nu}_\mu)$ and beam $\nu_e \ (\bar{\nu}_e)$ event spectrum. The Far Detector, located in Ash River, MN, USA, is 809 km from the ND and records the oscillated $\nu_e \ (\bar{\nu}_e)$ and the un-oscillated $\nu_\mu \ (\bar{\nu}_\mu)$ event spectrum. NOvA uses a data-driven technique called extrapolation to predict the expected number of $\nu_\mu \ (\bar{\nu}_\mu)$ and $\nu_e \ (\bar{\nu}_e)$ events at the Far Detector using the Near Detector data. The use of data from a functionally equivalent Near Detector provides a powerful constraint on the systematic uncertainties in NOvA neutrino oscillation analyses. As NOvA continues to add data statistics, a robust constraint on systematics becomes more crucial for neutrino oscillation analysis. The details of the NOvA neutrino oscillation analysis framework and how it constrains dominant systematic uncertainties using the Near Detector data will be discussed in this poster.

43 PARTICLE ACCELERATORS↗