Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Library”

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 1,783 records · Page 99

Asymptotic consistency of the WSINDy algorithm in the limit of continuum data

In this work we study the asymptotic consistency of the weak-form sparse identification of nonlinear dynamics algorithm (WSINDy) in the identification of differential equations from noisy samples of solutions. We prove that the WSINDy estimator is unconditionally asymptotically consistent for a wide class of models that includes the Navier–Stokes, Kuramoto–Sivashinsky and Sine–Gordon equations. We thus provide a mathematically rigorous explanation for the observed robustness to noise of weak-form equation learning. Conversely, we also show that, in general, the WSINDy estimator is only conditionally asymptotically consistent, yielding discovery of spurious terms with probability one if the noise level exceeds a critical threshold σ c . We provide explicit bounds on σ c in the case of Gaussian white noise and we explicitly characterize the spurious terms that arise in the case of trigonometric and/or polynomial libraries. Furthermore, we show that, if the data is suitably denoised (a simple moving average filter is sufficient), then asymptotic consistency is recovered for models with locally-Lipschitz, polynomial-growth nonlinearities. Our results reveal important aspects of weak-form equation learning, which may be used to improve future algorithms. We demonstrate our findings numerically using the Lorenz system, the cubic oscillator, a viscous Burgers-growth model and a Kuramoto–Sivashinsky-type high-order PDE.

asymptotic consistency↗

Fitness factors impacting survival of a subsurface bacterium in contaminated groundwater

Many factors contribute to the ability of a microbial species to persist when encountering complexly contaminated environments, including time of exposure, the nature and concentration of contaminants, availability of nutritional resources, and possession of a combination of appropriate molecular mechanisms needed for survival. Herein we sought to identify genes that are most important for survival of Gram-negative Enterobacteriaceae in contaminated groundwater environments containing high concentrations of nitrate and metals using the metal-tolerant Oak Ridge Reservation isolate, Pantoea sp. MT58 (MT58). Survival fitness experiments in which a randomly barcoded transposon insertion (RB-TnSeq) library of MT58 was exposed directly to contaminated Oak Ridge Reservation groundwater samples from across a nitrate and mixed metal contamination plume were used to identify genes important for survival with increasing exposure times and concentrations of contaminants, and availability of a carbon source. Genes involved in controlling and using carbon, encoding transcriptional regulators, and related to Gram-negative outer membrane processes were among those found to be important for survival in contaminated Oak Ridge Reservation groundwater. A comparative genomics analysis of 75 Pantoea genus strains allowed us to further separate the survival determinants into core and non-core genes in the Pantoea pangenome, revealing insights into the survival of subsurface microorganisms during contaminant plume intrusion.

59 BASIC BIOLOGICAL SCIENCES↗

Expanding the synthetic biology toolbox of Cupriavidus necator for establishing fatty acid production

Abstract The Gram-negative betaproteobacterium Cupriavidus necator is a chemolithotroph that can convert carbon dioxide into biomass. Cupriavidus necator has been engineered to produce a variety of high-value chemicals in the past. However, there is still a lack of a well-characterized toolbox for gene expression and genome engineering. Development and optimization of biosynthetic pathways in metabolically engineered microorganisms necessitates control of gene expression via functional genetic elements such as promoters, ribosome binding sites (RBSs), and codon optimization. In this work, a set of inducible and constitutive promoters were validated and characterized in C. necator, and a library of RBSs was designed and tested to show a 50-fold range of expression for green fluorescent protein (gfp). The effect of codon optimization on gene expression in C. necator was studied by expressing gfp and mCherry genes with varied codon-adaptation indices and was validated by expressing codon-optimized variants of a C12-specific fatty acid thioesterase to produce dodecanoic acid. We discuss further hurdles that will need to be overcome for C. necator to be widely used for biosynthetic processes.

59 BASIC BIOLOGICAL SCIENCES↗

Application of functional genomics for domestication of novel non-model microbes

Abstract With the expansion of domesticated microbes producing biomaterials and chemicals to support a growing circular bioeconomy, the variety of waste and sustainable substrates that can support microbial growth and production will also continue to expand. The diversity of these microbes also requires a range of compatible genetic tools to engineer improved robustness and economic viability. As we still do not fully understand the function of many genes in even highly studied model microbes, engineering improved microbial performance requires introducing genome-scale genetic modifications followed by screening or selecting mutants that enhance growth under prohibitive conditions encountered during production. These approaches include adaptive laboratory evolution, random or directed mutagenesis, transposon-mediated gene disruption, or CRISPR interference (CRISPRi). Although any of these approaches may be applicable for identifying engineering targets, here we focus on using CRISPRi to reduce the time required to engineer more robust microbes for industrial applications. One-Sentence Summary The development of genome scale CRISPR-based libraries in new microbes enables discovery of genetic factors linked to desired traits for engineering more robust microbial systems.

59 BASIC BIOLOGICAL SCIENCES↗

Thermophilic site-specific recombination system for rapid insertion of heterologous DNA into the Clostridium thermocellum chromosome

Clostridium thermocellum is an anaerobic thermophile capable of producing ethanol and other commodity chemicals from lignocellulosic biomass. The insertion of heterologous DNA into the C. thermocellum chromosome is currently achieved via a time-consuming homologous recombination process, where a single stable insertion can take 2–4 weeks or more to construct. In this work, we developed a thermostable version of the Serine recombinase Assisted Genome Engineering (tSAGE) approach for gene insertion in C. thermocellum utilizing a site-specific recombinase from Geobacillus sp. Y412MC61, enabling quick and easy insertion of DNA into the chromosome for accelerated genetic tool screening and heterologous gene expression. Using tSAGE, chromosomal insertion of plasmid DNA occurred at a maximum transformation efficiency of 5 × 10 3 CFU/µg, which is comparable to the transformation efficiency of a replicating control plasmid in C. thermocellum. Using tSAGE, we chromosomally integrated and characterized 17 reporter genes, 15 homologous and 31 heterologous constitutive promoters of varying strengths, 4 inducible promoters, and 5 riboswitches in C. thermocellum. We also determined that a 6–7 nucleotide gap between the ribosome binding site (RBS) and the start codon is optimal for high expression by employing a library of superfolder green fluorescent protein expression constructs driven by our strongest tested promoter (P clo1313_1194 ) with different distances between the RBS and start codon. The tools developed here will aid in accelerating C. thermocellum strain engineering for producing sustainable fuels and chemicals directly from plant biomass.

Biofuels↗

Single-cell and spatial omics in plants: from cellular atlases to regulatory mechanisms

Single-cell RNA sequencing (scRNA-seq) has transformed transcriptomic studies by enabling gene expression profiling at the resolution of individual cells within and across a broad range of tissue types, revealing cellular heterogeneity that is obscured in bulk tissue transcriptomes. Over the past decade, improvements in microfluidics and library preparation have drastically increased throughput, allowing tens of thousands of cells to be assayed in a single experiment. Although initially developed in animal systems, scRNA-seq has rapidly emerged as a powerful and widely adopted approach in plant biology. Beyond transcriptomics, the integration of single-cell data with chromatin accessibility, proteomics, metabolomics, and spatial omics is enabling a system-level understanding of plant gene regulation and cellular organization. Network-based analytical frameworks further support the reconstruction of gene regulatory networks and the interpretation of complex single-cell data. In this review, we summarize the current technological landscape of plant single-cell studies, discuss key experimental and analytical challenges, and review emerging strategies for validating single-cell discoveries. We also discuss future directions in applying single-cell technologies to woody perennials plants and bioenergy-relevant crops, emphasizing their potential to accelerate the discovery of cell type-specific regulatory mechanisms underlying growth, stress resilience, and biomass production.

Li, Miaomiao [ORNL] (ORCID:0000000321326168)↗

Random Forest Prediction of Crystal Structure from Electron Diffraction Patterns

Transmission electron microscopy (TEM) diffraction patterns are regularly used to determine the structure of crystalline materials. Electron diffraction is the most common method to solve for unknown or partially known crystal structures, as it provides direct and interpretable feedback on the orientation of crystal grains under the beam [1]. However, it remains a challenge to determine the crystal structure of a new material or even a new phase of an existing material. Analysis of such materials commonly requires manual exploration and comparison with simulated diffraction patterns. This is often a time consuming process with no obvious start point when many similar structures are possible, and this method cannot be used to determine crystal structure or orientation from structures not included in the diffraction libraries. Therefore, we have developed a machine learning model to determine the crystal structure of a material from its electron diffraction pattern.

36 MATERIALS SCIENCE↗

Evaluation of probabilistic photometric redshift estimation approaches for The Rubin Observatory Legacy Survey of Space and Time (LSST)

Many scientific investigations of photometric galaxy surveys require redshift estimates, whose uncertainty properties are best encapsulated by photometric redshift (photo-z) posterior probability density functions (PDFs). A plethora of photo-z PDF estimation methodologies abound, producing discrepant results with no consensus on a preferred approach. We present the results of a comprehensive experiment comparing twelve photo-z algorithms applied to mock data produced forLarge Synoptic Survey Telescope The Rubin Observatory Legacy Survey of Space and Time (lsst) Dark Energy Science Collaboration (desc). By supplying perfect prior information, in the form of the complete template library and a representative training set as inputs to each code, we demonstrate the impact of the assumptions underlying each technique on the output photo-z PDFs. In the absence of a notion of true, unbiased photo-z PDFs, we evaluate and interpret multiple metrics of the ensemble properties of the derived photo-z PDFs as well as traditional reductions to photo-z point estimates. Additionally, we report systematic biases and overall over/under-breadth of the photo-z PDFs of many popular codes, which may indicate avenues for improvement in the algorithms or implementations. Furthermore, we raise attention to the limitations of established metrics for assessing photo-z PDF accuracy; though we identify the conditional density estimate (CDE) loss as a promising metric of photo-z PDF performance in the case where true redshifts are available but true photo-z PDFs are not, we emphasize the need for science-specific performance metrics.

79 ASTRONOMY AND ASTROPHYSICS↗

iMaNGA: mock MaNGA galaxies based on IllustrisTNG and MaStar SSPs – I. Construction and analysis of the mock data cubes

ABSTRACT Galaxy formation and evolution simulations are essential tools to probe poorly known astrophysics processes, but particular care is needed to compare simulations with galaxy observations, as observed data need to be modelled as well. We present a method to generate mock galaxies from the hydro-dynamical IllustrisTNG simulations which are suited to compare with integral field spectroscopic observation of galaxies from the SDSS-IV/Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. First, we include the same instrumental effects and procedures as adopted in the acquisition and analysis of real data. Furthermore, we generate the galaxy spectra from the simulations using new stellar population models based on the MaNGA stellar library (MaStar). In this way, our mock data cubes have the same spatial sampling, cover the same wavelength range (3600–10 300 Å), and share the same spectral resolution (R ≈ 1800) and flux calibration of real MaNGA galaxy spectra. In this first paper, we demonstrate the method over an early- and a late-type simulated galaxy from TNG50. We analyse the correspondent mock MaNGA-like data cubes with the same full spectral fitting code, firefly, which was used for the observed spectra. We find that the intrinsic and recovered age and metallicity gradients are consistent within 1σ, with residuals over all tassels consistent with 0 at the 68 ${{\ \rm per\ cent}}$ confidence level. We also perform the challenging test at comparing intrinsic and recovered star formation histories, finding a close resemblance between input and output. In follow-up papers, we will present a full simulated MaNGA-like catalogue (≈10 000 galaxies) with a comprehensive comparison of TNG50 simulations to MaNGA observational results.

79 ASTRONOMY AND ASTROPHYSICS↗

Efficient analysis routines for single and double peaked Type 2 AGN spectra

Driven by the imminent need to rapidly process and classify millions of AGN spectra drawn from next generation astronomical facilities, we present a spectral fitting routine for Type 2 AGN spectra optimized for high volume processing, using the quasar spectral fitting library (qsfit). We analyse an optically selected sample of 813 luminous Type 2 AGN spectra at z < 0.83 from the Sloan Digital Sky Survey (SDSS) to qualify its performance. We report a median narrow line H α/H β Balmer decrement of 4.5 ± 0.8, alluding to the presence of dust in the narrow line region (NLR). We publish a specialized QSFIT fitting routine for high signal-to-noise ratio spectra and general fitting routine for double peaked Type 2 AGN spectra applied on a subsample of 45 spectra from our parent sample. We report a median red and blue peak velocity separation of 390 ± 60 kms −1 . No trend is found for red or blue peaks to exhibit systematically different luminosity or ionization properties. Emission line diagnostics show that the double peaks in all sources are illuminated by an AGN-powered ionizing continuum. Finally, we examine the morphology of host galaxies of our double peaked sample. We find double peaked Type 2 AGN reside in merging systems at a comparable frequency to single peaked AGN. This suggests that the double peaked AGN phenomenon is likely to have a bi-conical outflow origin in the majority of cases. We publicly release the code used for spectral analysis and produced catalogues used in this work.

79 ASTRONOMY AND ASTROPHYSICS↗

Colour gradients of low-redshift galaxies in the DESI Legacy Imaging Survey

ABSTRACT Radial colour gradients within galaxies arise from gradients of stellar age, metallicity, and dust reddening. Large samples of colour gradients from wide-area imaging surveys can complement smaller integral-field spectroscopy data sets and can be used to constrain galaxy formation models. Here, we measure colour gradients for low-redshift galaxies (z < 0.1) using photometry from the DESI Legacy Imaging Survey DR9. Our sample comprises ∼93 000 galaxies with spectroscopic redshifts and ∼574 000 galaxies with photometric redshifts. We focus on gradients across a radial range 0.5Reff to Reff, which corresponds to the inner disc of typical late-type systems at low redshift. This region has been the focus of previous statistical studies of colour gradients and has recently been explored by spectroscopic surveys such as MaNGA. We find that the colour gradients of most galaxies in our sample are negative (redder towards the centre), consistent with the literature. We investigate empirical relationships between colour gradient, average g − r and r − z colour, Mr, M⋆, and sSFR. Trends of gradient strength with Mr (M⋆) show an inflection around Mr ∼ −21 ($\log _{10} \, M_\star /\mathrm{M_\odot }\sim 10.5$). Below this mass, colour gradients become steeper with increasing M⋆, whereas colour gradients in more massive galaxies become shallower. We find that positive gradients (bluer stars at smaller radii) are typical for galaxies of $M_{\star }\sim 10^{8}\, \mathrm{M_\odot }$. We compare our results to age and metallicity gradients in two data sets derived from fits of different stellar population libraries to MaNGA spectra, but find no clear consensus explanation for the trends we observe. Both MaNGA data sets seem to imply a significant contribution from dust reddening, in particular, to explain the flatness of colour gradients along the red sequence.

Astronomy & Astrophysics↗

A differentiable model of the evolution of dark matter halo concentration

ABSTRACT We introduce a new model of the evolution of the concentration of dark matter haloes, c(t). For individual haloes, our model approximates c(t) as a power law with a time-dependent index, such that at early times, concentration has a nearly constant value of c ≈ 3–4, and as cosmic time progresses, c(t) smoothly increases. Using large samples of halo merger trees taken from the Bolshoi–Planck and MultiDark Planck 2 cosmological simulations, we demonstrate that our three-parameter model can approximate the evolution of the concentration of individual haloes with a typical accuracy of 0.1 dex for $t\gtrsim 2\, {\rm Gyr}$ for all Bolshoi–Planck and MultiDark Planck 2 haloes of present-day peak mass $M_{0}\gtrsim 10^{11.5}\, {\rm M}_{\odot }$. We additionally present a new model of the evolution of the concentration of halo populations, which we show faithfully reproduces both average concentration growth and the diversity of smooth trajectories of c(t), including capturing correlations with halo mass and halo assembly history. Our publicly available source code, diffprof, can be used to generate Monte Carlo realizations of the concentration histories of cosmologically representative halo populations. diffprof is differentiable due to its implementation in the jax autodiff library, which facilitates the incorporation of our model into existing analytical halo model frameworks.

Stevanovich, Dash↗

The simulated catalogue of optical transients and correlated hosts (SCOTCH)

ABSTRACT As we observe a rapidly growing number of astrophysical transients, we learn more about the diverse host galaxy environments in which they occur. Host galaxy information can be used to purify samples of cosmological Type Ia supernovae, uncover the progenitor systems of individual classes, and facilitate low-latency follow-up of rare and peculiar explosions. In this work, we develop a novel data-driven methodology to simulate the time-domain sky that includes detailed modelling of the probability density function for multiple transient classes conditioned on host galaxy magnitudes, colours, star formation rates, and masses. We have designed these simulations to optimize photometric classification and analysis in upcoming large synoptic surveys. We integrate host galaxy information into the snana simulation framework to construct the simulated catalogue of optical transients and correlated hosts (SCOTCH, a publicly available catalogue of 5-million idealized transient light curves in LSST passbands and their host galaxy properties over the redshift range 0 < z < 3. This catalogue includes supernovae, tidal disruption events, kilonovae, and active galactic nuclei. Each light curve consists of true top-of-the-galaxy magnitudes sampled with high (≲2 d) cadence. In conjunction with SCOTCH, we also release an associated set of tutorials and transient-specific libraries to enable simulations of arbitrary space- and ground-based surveys. Our methodology is being used to test critical science infrastructure in advance of surveys by the Vera C. Rubin Observatory and the Nancy G. Roman Space Telescope.

Astronomy & Astrophysics↗

A resolved study of the inner regions of nearby galaxies with an excess of young massive stars: missing link in the AGN–starburst connection?

ABSTRACT We have selected galaxies with very high levels of H α emission (EQW(H α) >700 Å.) in their central regions from the final data release of the MaNGA survey. Our study focuses on 14 very well-resolved nearby galaxies with stellar masses in the range 9.5 < log M*/(M⊙) < 11.5. We investigate a variety of procedures for selecting galaxy regions that are likely to harbour excess populations of young massive stars, finding that selection in the 2D space of extinction-corrected H α EQW and [S iii]/[S ii] line ratio produces the best results. By comparing stacked spectra covering these regions with stacked spectra covering normal starburst regions with 100 Å<EQW(H α) <200 Å, we obtain the following main results: (1) Clear signatures of excess Wolf–Rayet stars are found in half of the H α excess regions, (2) galaxy regions containing excess Wolf–Rayet stars are more often associated with the presence of high-ionization emission lines characteristic of accreting black holes. Excess [Ne iii] is detected in 4 out of 8 of the WR regions and there are tentative [Fe x] detections in 2 galaxies. (3) Regions of the galaxy with excess Wolf–Rayet stars are located where the interstellar medium has larger ionized gas turbulent velocities and higher neutral gas overdensities. We make a first attempt to constrain changes in the high-mass end of the stellar initial mass function (IMF) using the HR-pyPopStar evolutionary population synthesis models that include high-wavelength-resolution theoretical atmosphere libraries for Wolf–Rayet stars.

Kauffmann, Guinevere↗

DSPS: Differentiable stellar population synthesis

ABSTRACT Models of stellar population synthesis (SPS) are the fundamental tool that relates the physical properties of a galaxy to its spectral energy distribution (SED). In this paper, we present DSPS: a python package for SPS. All of the functionality in DSPS is implemented natively in the JAX library for automatic differentiation, and so our predictions for galaxy photometry are fully differentiable, and directly inherit the performance benefits of JAX, including portability onto GPUs. DSPS also implements several novel features, such as i) a flexible empirical model for stellar metallicity that incorporates correlations with stellar age, ii) support for the Diffstar model that provides a physically-motivated connection between the star formation history of a galaxy (SFH) and the mass assembly of its underlying dark matter halo. We detail a set of theoretical techniques for using autodiff to calculate gradients of predictions for galaxy SEDs with respect to SPS parameters that control a range of physical effects, including SFH, stellar metallicity, nebular emission, and dust attenuation. When forward modelling the colours of a synthetic galaxy population, we find that DSPS can provide a factor of 5 speed-up over standard SPS codes on a CPU, and a factor of 300-400 on a modern GPU. When coupled with gradient-based techniques for optimization and inference, DSPS makes it practical to conduct expansive likelihood analyses of simulation-based models of the galaxy–halo connection that fully forward model galaxy spectra and photometry.

79 ASTRONOMY AND ASTROPHYSICS↗

Emission line predictions for mock galaxy catalogues: a new differentiable and empirical mapping from DESI

ABSTRACT We present a simple, differentiable method for predicting emission line strengths from rest-frame optical continua using an empirically determined mapping. Extensive work has been done to develop mock galaxy catalogues that include robust predictions for galaxy photometry, but reliably predicting the strengths of emission lines has remained challenging. Our new mapping is a simple neural network implemented using the JAX Python automatic differentiation library. It is trained on Dark Energy Spectroscopic Instrument Early Release data to predict the equivalent widths (EWs) of the eight brightest optical emission lines (including H α, H β, [O ii], and [O iii]) from a galaxy’s rest-frame optical continuum. The predicted EW distributions are consistent with the observed ones when noise is accounted for, and we find Spearman’s rank correlation coefficient ρs > 0.87 between predictions and observations for most lines. Using a non-linear dimensionality reduction technique, we show that this is true for galaxies across the full range of observed spectral energy distributions. In addition, we find that adding measurement uncertainties to the predicted line strengths is essential for reproducing the distribution of observed line-ratios in the BPT diagram. Our trained network can easily be incorporated into a differentiable stellar population synthesis pipeline without hindering differentiability or scalability with GPUs. A synthetic catalogue generated with such a pipeline can be used to characterize and account for biases in the spectroscopic training sets used for training and calibration of photo-z’s, improving the modelling of systematic incompleteness for the Rubin Observatory LSST and other surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

guppy i : a code for reducing the storage requirements of cosmological simulations

ABSTRACT As cosmological simulations have grown in size, the permanent storage requirements of their particle data have also grown. Even modest simulations present a major logistical challenge for the groups which run these boxes and researchers without access to high performance computing facilities often need to restrict their analysis to lower quality data. In this paper, we present guppy, a compression algorithm and code base tailored to reduce the sizes of dark matter-only cosmological simulations by approximately an order of magnitude. guppy is a ‘lossy’ algorithm, meaning that it injects a small amount of controlled and uncorrelated noise into particle properties. We perform extensive tests on the impact that this noise has on the internal structure of dark matter haloes, and identify conservative accuracy limits which ensure that compression has no practical impact on single-snapshot halo properties, profiles, and abundances. We also release functional prototype libraries in C, Python, and Go for reading and creating guppy data.

79 ASTRONOMY AND ASTROPHYSICS↗

Carbon- and Oxygen-rich stars in MaStar: identification and classification

ABSTRACT Carbon- and Oxygen-rich stars populating the Thermally Pulsing Asymptotic Giant Branch (TP-AGB) phase of stellar evolution are relevant contributors to the spectra of ∼1 Gyr old populations. Atmosphere models for these types are uncertain, due to complex molecules and mass-loss effects. Empirical spectra are then crucial, but samples are small due to the short (∼3 Myr) TP-AGB lifetime. Here we exploit the vastness of the MaNGA Stellar library MaStar (∼60 000 spectra) to identify C and O-rich type stars. We define an optical colour selection with cuts of (g − r) > 2 and (g − i) < 1.55(g − r) − 0.07, calibrated with known C- and O-rich spectra. This identifies C- and O-rich stars along clean, separated sequences. An analogue selection is found in V, R, I bands. Our equation identifies C-rich and O-rich spectra with predictive performance metric F1-scores of 0.72 and 0.74 (over 1), respectively. We finally identify 41 C- and 87 O-rich type AGB stars in MaStar, 5 and 49 of which do not have a SIMBAD counterpart. We also detect a sample of non-AGB, dwarf C-stars. We further design a fitting procedure to classify the spectra into broad spectral types, by using as fitting templates empirical C- and O-rich spectra. We find remarkably good fits for the majority of candidates and categorize them into C- and O-rich bins following existing classifications, which correlate to effective temperature. Our selection models can be applied to large photometric surveys (e.g. Euclid, Rubin). The classified spectra will facilitate future evolutionary population synthesis models.

Hill, Lewis (ORCID:0000000180528080)↗