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 775 records · Page 43

The Milky Way’s bulge star formation history as constrained from its bimodal chemical abundance distribution

ABSTRACT We conduct a quantitative analysis of the star formation history (SFH) of the Milky Way’s (MW) bulge by exploiting the constraining power of its stellar [Fe/H] and [Mg/Fe] distribution functions. Using Apache Point Observatory Galactic Evolution Experiment survey data, we confirm the previously established bimodal [Mg/Fe]–[Fe/H] distribution within 3 kpc of the inner Galaxy. To fit the chemical bimodal distribution, we use a simple but flexible star formation framework, which assumes two distinct stages of gas accretion and star formation, and systematically evaluate a wide multidimensional parameter space. We find that the data favour a three-phase SFH that consists of an initial starburst, followed by a rapid star formation quenching episode, and a lengthy, quiescent secular evolution phase. The metal-poor, high-α bulge stars ([Fe/H] < 0.0 and [Mg/Fe] > 0.15) are formed rapidly (<2 Gyr) during the early starburst. The density gap between the high- and low-α sequences is due to the quenching process. The metal-rich, low-α population ([Fe/H] > 0.0 and [Mg/Fe] < 0.15) then accumulates gradually through inefficient star formation during the secular phase. This is qualitatively consistent with the early SFH of the inner disc. Given this scenario, a notable fraction of young stars (age <5 Gyr) is expected to persist in the bulge. Combined with extragalactic observations, these results suggest that a rapid star formation quenching process is responsible for bimodal distributions in both the MW’s stellar populations and in the general galaxy population and thus plays a critical role in galaxy evolution.

Lian, Jianhui↗

Dark Energy Survey Year 3 results: cosmology with moments of weak lensing mass maps – validation on simulations

ABSTRACT We present a simulated cosmology analysis using the second and third moments of the weak lensing mass (convergence) maps. The second moment, or variances, of the convergence as a function of smoothing scale contains information similar to standard shear two-point statistics. The third moment, or the skewness, contains additional non-Gaussian information. The analysis is geared towards the third year (Y3) data from the Dark Energy Survey (DES), but the methodology can be applied to other weak lensing data sets. We present the formalism for obtaining the convergence maps from the measured shear and for obtaining the second and third moments of these maps given partial sky coverage. We estimate the covariance matrix from a large suite of numerical simulations. We test our pipeline through a simulated likelihood analyses varying 5 cosmological parameters and 10 nuisance parameters and identify the scales where systematic or modelling uncertainties are not expected to affect the cosmological analysis. Our simulated likelihood analysis shows that the combination of second and third moments provides a 1.5 per cent constraint on S8 ≡ σ8(Ωm/0.3)0.5 for DES Year 3 data. This is 20 per cent better than an analysis using a simulated DES Y3 shear two-point statistics, owing to the non-Gaussian information captured by the inclusion of higher order statistics. This paper validates our methodology for constraining cosmology with DES Year 3 data, which will be presented in a subsequent paper.

79 ASTRONOMY AND ASTROPHYSICS↗

Dark energy survey year 3 results: cosmological constraints from the analysis of cosmic shear in harmonic space

ABSTRACT We present cosmological constraints from the analysis of angular power spectra of cosmic shear maps based on data from the first three years of observations by the Dark Energy Survey (DES Y3). Our measurements are based on the pseudo-Cℓ method and complement the analysis of the two-point correlation functions in real space, as the two estimators are known to compress and select Gaussian information in different ways, due to scale cuts. They may also be differently affected by systematic effects and theoretical uncertainties, making this analysis an important cross-check. Using the same fiducial Lambda cold dark matter model as in the DES Y3 real-space analysis, we find ${S_8 \equiv \sigma _8 \sqrt{\Omega _{\rm m}/0.3} = 0.793^{+0.038}_{-0.025}}$, which further improves to S8 = 0.784 ± 0.026 when including shear ratios. This result is within expected statistical fluctuations from the real-space constraint, and in agreement with DES Y3 analyses of non-Gaussian statistics, but favours a slightly higher value of S8, which reduces the tension with the Planck 2018 constraints from 2.3σ in the real space analysis to 1.5σ here. We explore less conservative intrinsic alignments models than the one adopted in our fiducial analysis, finding no clear preference for a more complex model. We also include small scales, using an increased Fourier mode cut-off up to $k_{\rm max}={5}\, {h}\, {\rm Mpc}^{-1}$, which allows to constrain baryonic feedback while leaving cosmological constraints essentially unchanged. Finally, we present an approximate reconstruction of the linear matter power spectrum at present time, found to be about 20 per cent lower than predicted by Planck 2018, as reflected by the lower S8 value.

79 ASTRONOMY AND ASTROPHYSICS↗

Why Dissolving Salt in Water Decreases Its Dielectric Permittivity

The dielectric permittivity of salt water decreases on dissolving more salt. For nearly a century, this phenomenon has been explained by invoking saturation in the dielectric response of the solvent water molecules. Herein, we employ an advanced deep neural network (DNN), built using data from density functional theory, to study the dielectric permittivity of sodium chloride solutions. Notably, the decrease in the dielectric permittivity as a function of concentration, computed using the DNN approach, agrees well with experiments. Detailed analysis of the computations reveals that the dominant effect, caused by the intrusion of ionic hydration shells into the solvent hydrogen-bond network, is the disruption of dipolar correlations among water molecules. Accordingly, the observed decrease in the dielectric permittivity is mostly due to increasing suppression of the collective response of solvent waters.

74 ATOMIC AND MOLECULAR PHYSICS↗

Information-Theoretic Exploration of Multivariate Time-Varying Image Databases

Modern scientific simulations produce very large datasets, making interactive exploration of such data computationally prohibitive. An increasingly common data reduction technique is to store visualizations and other data extracts in a database. The Cinema project is one such approach, storing visualizations in an image database for post hoc exploration and interactive image-based analysis. This work focuses on developing efficient algorithms that can quantify various types of multivariate dependencies existing within multi-variable datasets. It applies specific mutual information measures for the quantification of salient regions from multivariate image data. Here, using such information measures, the opacity of the images is modulated so that the salient regions are automatically highlighted and the domain scientists can interactively explore the most relevant regions for scientific discovery.

97 MATHEMATICS AND COMPUTING↗

A shift between mineral and nonmineral sources of iron and sulfur causes proteome-wide changes in Methanosarcina barkeri

ABSTRACT Iron (Fe) and sulfur (S) are required elements for life, and changes in their availability can limit the ecological distribution and function of microorganisms. In anoxic environments, soluble Fe typically exists as ferrous iron [Fe(II)] and S as sulfide (HS − ). These species exhibit a strong affinity that ultimately drives the formation of sedimentary pyrite (FeS 2 ). Recently, paradigm-shifting studies indicate that Fe and S in FeS 2 can be made bioavailable by methanogens through a reductive dissolution process. However, the impact of the utilization of FeS 2 , as opposed to canonical Fe and S sources, on the phenotype of cells is not fully understood. Here, shotgun proteomics was utilized to measure changes in the phenotype of Methanosarcina barkeri MS grown with FeS 2 , Fe(II)/HS − , or Fe(II)/cysteine. Shotgun proteomics tracked 1,019 proteins overall, with 307 observed to change between growth conditions. Functional characterization and pathway analyses revealed these changes to be systemic and largely tangential to Fe/S metabolism. As a final step, the proteomics data were viewed with respect to previously collected transcriptomics data to deepen the analysis. Presented here is evidence that M. barkeri adopts distinct phenotypes to exploit specific sources of Fe and S in its environment. This is supported by observed protein abundance changes across broad categories of cellular biology. DNA adjacent metabolism, central carbon metabolism methanogenesis, metal trafficking, quorum sensing, and porphyrin biosynthesis pathways are all features in the phenotypic differentiation. Differences in trace metal availability attributed to complexation with HS − , either as a component of the growth medium [Fe(II)/HS − ] or generated through reduction of FeS 2 , were likely a major factor underpinning these phenotypic differences. IMPORTANCE The methanogenic archaeon Methanosarcina barkeri holds great potential for industrial bio-mining and energy generation technologies. Much of the biochemistry of this microbe is poorly understood, and its characterization will provide a glimpse into biological processes that evolved close to life’s origin. The discovery of its ability to extract iron and sulfur from bulk, solid-phase minerals shifted a longstanding paradigm that these elements were inaccessible to biological systems. The full elucidation of this process has the potential to help scientists and engineers extract valuable metals from low-grade ore and mine waste generating energy in the form of methane while doing so.

59 BASIC BIOLOGICAL SCIENCES↗

ALGEBRA v.1.27

The ALGEBRA program allows the user to manipulate data from a finite element analysis before it is plotted. The finite element output data is in the form of variable values (e.g., stress, strain, and velocity components) in an EXODUS II database. The ALGEBRA program evaluates user-supplied functions of the data and writes the results to an output EXODUS II database that can be read by plot programs.

Sjaardema, Gregory↗

Topsoil bulk geochemical compositions - An updated harmonized global dataset

Mineral weathering is a key biogeochemical process because of the capacity of minerals to stabilize organic matter. However, predicting soil weathering status across large spatial areas still isn’t possible due to a lack of global data and theoretical frameworks. To address this knowledge gap, multiple global datasets of bulk topsoil geochemical compositions have been harmonized using R. These datasets document topsoil bulk geochemical compositions across five continents (n = ~16,000 observations). Source data for these observations include the EuroGEOSurveys Geochemical Baseline Database (FOREGS), the US Geological Survey National Geochemical Database (NASGLP), the Geochemical Atlas of Australia (GAA), the US Geological Survey Alaska Geochemical Database (AGD84), the National Cooperative Soil Survey (NCSS), the European Geochemical Mapping of Agricultural Soil (GEMAS), Ecorespira-Amazon (ERA), the New Zealand Geochemical Baseline Survey (NZ_GBS), and the African Soil Information Service (AFSIS). Major elements observed include Aluminum (Al), Calcium (Ca), Iron (Fe), Potassium (K), Magnesium (Mg), Sodium (Na), Titanium (Ti), Manganese (Mn), Phosphorus (P), Carbon (C), and Sulfur (S). This data package includes the harmonized dataset itself, and the R scripts necessary to harmonize these datasets, in addition to metadata that describes all columns, files, and databases used in this project. Methods & Sampling Step 1 – Databases of geochemical data identified This study aimed to leverage existing measurements of topsoil geochemical data. Databases were first identified and deemed appropriate for inclusion if they were measuring soils and performed these measurements on the <2mm soil fraction. Databases such as NCSS and AGD84 needed more post processing to include in the database and this was done using the NCSS_datamerge_031626 R file and Alaska_USGSmerge_031626 R file, respectively. Step 2 – Database harmonization Once appropriate databases were identified, they were harmonized for ease of analysis using the R script Database_Harmonization_031826. This included removing columns from original datasets that would not be used in analysis (removed columns are noted in the code). Then, data cleaning procedures specific to each dataset were undertaken. This includes standardizing columns to include units and adding metadata columns regarding procedures for analyzing specific elements. Functions for standardizing measurements and units are outline in R files: calculate element_mg_kg_031626, calculate_oxide_wt_perc_031626, change_oxide_caps_031626, and conv_2_numeric_031626. This also included adding a unique identifier for each sample to identify it with its respective database (see CD_ID in data dictionary). Geographic information: Data reflect a compilation of datasets collected globally. Geographic areas covered by each of the datasets include: - EuroGEOSurveys Geochemical Baseline Database (FOREGS) - European continent - North American Soil Geochemical Landscapes (NASGLP) - continental United States and limited parts of Canada (see database key for more details) - National Geochemical Survey of Australia (GAA) - Australia - Alaska geochemical database (AGDB4) - Alaska - National Cooperative Soil Survey (NCSS) - Global measurements, but concentrated in the continental United States - Geochemical data for arable land and land under permanent grass cover in continental Europe (GEMAS) - continental Europe - Ecorespira-Amazon (ERA) - Geochemical data from the Amazon basin - Geochemical baseline data for New Zealand (NZGBS) - New Zealand - Geochemical data collected across continental Africa (AfSIS) - Measurements across Africa

EARTH SCIENCE > LAND SURFACE > SOILS↗

Charged point defect benchmark of Hybrid and GGA-PBE

Data used by the publication "High-throughput calculations of charged point defect properties with semi-local density functional theory - performance benchmarks for materials screening applications." This work presented an in-depth benchmark analysis of automated, semi-local point defect calculations with a-posteriori corrections, compared to 245 “gold standard” hybrid calculations previously published. We considered three different a-posteriori correction sets for semi-local calculations, implemented in a fully automated workflow, and consider the qualitative and quantitative differences for four different categories of defect information: thermodynamic transition levels, formation energies, fermi levels, and dopability limits. We highlighted the type of qualitative information about point defect properties that can be extracted from high-throughput calculations based on semi-local DFT methods, while also demonstrating the limits of quantitative accuracy that can be achieved by these approaches.

Broberg, Danny↗

KBase Narrative - Genomic and environmental controls on Castellaniella biogeography in an anthropogenically disturbed site

Genome assemblies were imported into KBase using the Batch Import Assembly from Staging Area (v1.0.57) function. All assemblies were annotated using the Annotated Multiple Microbial Assemblies with RASTtk - v1.073 tool. Annotated genomes were grouped into sets using the Add Genomes to GenomeSet - v1.7.6 function. Individual annotated genomes can be found both below and in the Data menu to the left. Taxonomy was assigned using the Classify Microbes with GTDB-Tk-v1.7.0 tool. The results of this analysis are shown below. Analysis of the Castellaniella pangenome was performed using the Compute Pangenome (v0.0.7) tool. Using the same method, we also computed the ORR-specific and non-ORR Castellaniella pangenomes. All pangenome results (including the presence/absence matrix) can be found below.

Szink, Elizabeth↗

Biological insights from multi-omics analysis strategies: Complex pleotropic effects associated with autophagy

Research strategies that combine molecular data from multiple levels of genome expression (i.e., multi-omics data), often referred to as a systems biology strategy, has been advocated as a route to discovering gene functions. In this study we conducted an evaluation of this strategy by combining lipidomics, metabolite mass-spectral imaging and transcriptomics data from leaves and roots in response to mutations in two AuTophaGy-related ( ATG ) genes of Arabidopsis . Autophagy is an essential cellular process that degrades and recycles macromolecules and organelles, and this process is blocked in the atg7 and atg9 mutants that were the focus of this study. Specifically, we quantified abundances of ~100 lipids and imaged the cellular locations of ~15 lipid molecular species and the relative abundance of ~26,000 transcripts from leaf and root tissues of WT, atg7 and atg9 mutant plants, grown either in normal (nitrogen-replete) and autophagy-inducing conditions (nitrogen-deficient). The multi-omics data enabled detailed molecular depiction of the effect of each mutation, and a comprehensive physiological model to explain the consequence of these genetic and environmental changes in autophagy is greatly facilitated by the a priori knowledge of the exact biochemical function of the ATG7 and ATG9 proteins.

59 BASIC BIOLOGICAL SCIENCES↗

A Computational Study of the Shear Behavior of Reinforced Concrete Beams Affected from Alkali–Silica Reactivity Damage

In this study, an investigation of the shear behavior of full-scale reinforced concrete (RC) beams affected from alkali–silica reactivity damage is presented. A detailed finite element model (FEM) was developed and validated with data obtained from the experiments using several metrics, including a force–deformation curve, rebar strains, and crack maps and width. The validated FEM was used in a parametric study to investigate the potential impact of alkali–silica reactivity (ASR) degradation on the shear capacity of the beam. Degradations of concrete mechanical properties were correlated with ASR expansion using material test data and implemented in the FEM for different expansions. The finite element (FE) analysis provided a better understanding of the failure mechanism of ASR-affected RC beam and degradation in the capacity as a function of the ASR expansion. The parametric study using the FEM showed 6%, 19%, and 25% reduction in the shear capacity of the beam, respectively, affected from 0.2%, 0.4%, and 0.6% of ASR-induced expansion.

36 MATERIALS SCIENCE↗

Pantropical Tree Leaf Hydraulic Properties, 1986 - 2015

This data package contains leaf pressure-volume parameters (osmotic potential at full turgor, turgor loss point, bulk elastic modulus, and others). Data was obtained by combining new data I gleaned from original publications not already in the following compilations (after removing the extratropical data): Bartlett, M. K., Zhang, Y., Kreidler, N., Sun, S., Ardy, R., Cao, K., and Sack, L.: Global analysis of plasticity in turgor loss point, a key drought tolerance trait, Ecol Lett, 17, 1580-1590, 2014. Bartlett, M. K., Scoffoni, C., and Sack, L.: The determinants of leaf turgor loss point and prediction of drought tolerance of species and biomes: a global meta-analysis, Ecology Letters, 15, 393-405, 2012. Maréchaux, I., Bartlett, M. K., Sack, L., Baraloto, C., Engel, J., Joetzjer, E., Chave, J., and Kitajima, K.: Drought tolerance as predicted by leaf water potential at turgor loss point varies strongly across species within an Amazonian forest, Functional Ecology, 29, 1268-1277, 2015. Data are from natural tropical and subtropical forests or savannas and span rainfall gradients from deciduous dry tropical forests to everwet tropical forests, and come from countries all throughout the pantropics. The attached zip file contains data organized into three CSV files.

54 ENVIRONMENTAL SCIENCES↗

Selection function of clusters in Dark Energy Survey year 3 data from cross-matching with South Pole Telescope detections

Context. Galaxy clusters selected based on overdensities of galaxies in photometric surveys provide the largest cluster samples. However, modeling the selection function of such samples is complicated by noncluster members projected along the line of sight (projection effects) and the potential detection of unvirialized objects (contamination). Aims. We empirically constrained the magnitude of these effects by cross-matching galaxy clusters selected in the Dark Energy Survey data with the redMaPPer algorithm with significant detections in three South Pole Telescope surveys (SZ, pol-ECS, pol-500d). Methods. For matched clusters, we augmented the redMaPPer catalog with the SPT detection significance. For unmatched objects we used the SPT detection threshold as an upper limit on the SZe signature. Using a Bayesian population model applied to the collected multiwavelength data, we explored various physically motivated models to describe the relationship between observed richness and halo mass. Results. Our analysis reveals a clear preference for models with an additional skewed scatter component associated with projection effects over a purely log-normal scatter model. We rule out significant contamination by unvirialized objects at the high-richness end of the sample. While dedicated simulations offer a well-fitting calibration of projection effects, our findings suggest the presence of redshift-dependent trends that these simulations may not have captured. Our findings highlight that modeling the selection function of optically detected clusters remains a complicated challenge that requires a combination of simulation and data-driven approaches.

79 ASTRONOMY AND ASTROPHYSICS↗

Estimating SHmax azimuth with P sources and vertical geophones: Use P-P reflection amplitudes or use SV-P reflection times?

We compared two methods for extracting the azimuth of maximum horizontal stress (SHmax) from 3D land-based seismic data generated by a P source and recorded with vertical geophones. In the first method, we used the direct-SV mode that is produced by all land-based P sources. P sources generate SV illumination that radiates in all azimuth directions from a source station and creates SV-P reflections that are recorded by vertical geophones. Unless stratigraphy has steep dip, SV-P raypaths recorded by vertical geophones are the reverse of P-SV raypaths recorded by horizontal geophones. Thus, SV-P data provide the same S-wave sensitivity to stress fields as popular P-SV data do. In the second method, we retrieved P-P reflections and then performed an amplitude-variation-with-azimuth (AVA) analysis of the amplitude-gradient behavior of P-P reflection wavelets. We did this analysis in narrow azimuth corridors to determine the gradient of reflection-wavelet amplitudes as a function of azimuth. This P-P AVA amplitude-gradient method has been of great interest in the reflection seismology community since it was introduced in the late 1990s. Each of these methods, AVA analysis of the gradient of P-P reflection amplitudes and azimuth-dependent arrival times of SV-P reflections, can be used to determine the azimuth of SHmax stress. We compare the results of the two methods with ground truth measurements of SHmax azimuth at a CO 2 sequestration site in the Michigan Basin. SHmax azimuths were determined from P-P and SV-P data at three major boundaries at depths of approximately 3500 ft (1067 m), 5500 ft (1676 m), and 7500 ft (2286 m). Two estimates of SHmax azimuth (one using SV-P data and one using P-P data) were made at each stacking bin inside a 24 mi 2 (62 km 2 ) image space. The result was approximately 98,000 estimates of SHmax azimuth across each of these three boundaries for each of these two prediction strategies. Histogram displays of PP AVA gradient estimates had peaks at correct azimuths of SHmax at all three depths, but the spread of the distributions widened with depth and split into two peaks at the deepest boundary. In contrast, each histogram of SHmax azimuth predicted by azimuth-dependent SV-P traveltimes had a single, definitive peak that was positioned at the correct SHmax azimuth at all three boundary depths.

Geochemistry & Geophysics↗

Physics in the Machine: Integrating Physical Knowledge in Autonomous Phase-Mapping

Application of artificial intelligence (AI), and more specifically machine learning, to the physical sciences has expanded significantly over the past decades. In particular, science-informed AI, also known as scientific AI or inductive bias AI, has grown from a focus on data analysis to now controlling experiment design, simulation, execution and analysis in closed-loop autonomous systems. The CAMEO (closed-loop autonomous materials exploration and optimization) algorithm employs scientific AI to address two tasks: learning a material system’s composition-structure relationship and identifying materials compositions with optimal functional properties. By integrating these, accelerated materials screening across compositional phase diagrams was demonstrated, resulting in the discovery of a best-in-class phase change memory material. Key to this success is the ability to guide subsequent measurements to maximize knowledge of the composition-structure relationship, or phase map. In this work we investigate the benefits of incorporating varying levels of prior physical knowledge into CAMEO’s autonomous phase-mapping. This includes the use of ab-initio phase boundary data from the AFLOW repositories, which has been shown to optimize CAMEO’s search when used as a prior.

97 MATHEMATICS AND COMPUTING↗