Search NASA⌕ Search

SEARCH · Search NASA

Results for “Level set method”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Disentangling Sources of Gene Tree Discordance in Phylogenomic Data Sets: Testing Ancient Hybridizations in Amaranthaceae s.l

Gene tree discordance in large genomic data sets can be caused by evolutionary processes such as incomplete lineage sorting and hybridization, as well as model violation, and errors in data processing, orthology inference, and gene tree estimation. Species tree methods that identify and accommodate all sources of conflict are not available, but a combination of multiple approaches can help tease apart alternative sources of conflict. Here, using a phylotranscriptomic analysis in combination with reference genomes, we test a hypothesis of ancient hybridization events within the plant family Amaranthaceae s.l. that was previously supported by morphological, ecological, and Sanger-based molecular data. The data set included seven genomes and 88 transcriptomes, 17 generated for this study. We examined gene-tree discordance using coalescent-based species trees and network inference, gene tree discordance analyses, site pattern tests of introgression, topology tests, synteny analyses, and simulations. We found that a combination of processes might have generated the high levels of gene tree discordance in the backbone of Amaranthaceae s.l. Furthermore, we found evidence that three consecutive short internal branches produce anomalous trees contributing to the discordance. Overall, our results suggest that Amaranthaceae s.l. might be a product of an ancient and rapid lineage diversification, and remains, and probably will remain, unresolved. This work highlights the potential problems of identifiability associated with the sources of gene tree discordance including, in particular, phylogenetic network methods. Our results also demonstrate the importance of thoroughly testing for multiple sources of conflict in phylogenomic analyses, especially in the context of ancient, rapid radiations. We provide several recommendations for exploring conflicting signals in such situations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Sizing ramping reserve using probabilistic solar forecasts: A data-driven method

Ramping products have been introduced or proposed in several U.S. power markets to mitigate the impact of load and renewable uncertainties on market efficiency and reliability. Current methods often rely on historical data to estimate the requirements of ramping products and fail to take into account the effects of the latest weather conditions and their uncertainties, which could lead to overly conservative or insufficient requirements. This study proposes a k-nearest-neighbor-based method to give weather-informed estimates of ramping needs based on short-term probabilistic solar irradiance forecasts. Forecasts from multiple sites are employed in conjunction with principal component analysis to derive numerical classifiers to characterize system-level weather conditions. In addition, we develop a data-driven method to optimize the model parameters in a rolling-forward manner. By using real-world data from the California Independent System Operator, we design two metrics to evaluate method performance: 1) frequency of shortage and 2) oversupply of ramping product. Our proposed method presents advantages in comparison with the baseline and a set of benchmark methods: without compromising system reliability, it reduces system ramping requirements by up to 25%, therefore improving both system reliability and economics.

14 SOLAR ENERGY↗
DOE OSTI2018

On the geometric dependence of the molecular dipole polarizability in water: A benchmark study of higher-order electron correlation, basis set incompleteness error, core electron effects, and zero-point vibrational contributions

In this work, we investigate how geometric changes influence the static dipole polarizability ($α$) of a water molecule by explicitly computing the corresponding dipole polarizability surface (DPS) across 3125 total (1625 symmetry-unique) geometries using linear response coupled cluster theory including single, double, and triple excitations (LR-CCSDT) and the doubly augmented triple-$ζ$ basis set (d-aug-cc-pVTZ). Analytical formulae based on power series expansions of this $ab initio$ surface are generated using linear least-squares analysis and provide highly accurate estimates of this quantity as a function of molecular geometry (i.e., bond and angle variations) in a computationally tractable manner. An additional database, which consists of 25 representative molecular geometries and incorporates a more thorough treatment of both basis sets and core electron effects, is provided as a current benchmark for this quantity and the corresponding leading-order $C_6$ dispersion coefficient. This database has been utilized to assess the importance of these effects as well as the relative accuracy that can be obtained using several quantum chemical methods and a library of density functional approximations. In addition to high-level electron correlation methods (like CCSD) and our analytical least-squares formulae, we find that the SCAN0, PBE0, MN15, and B97-2 hybrid functionals yield the most accurate descriptions of the molecular polarizability tensor in H 2 O. Using first-order perturbation theory, we compute the zero-point vibrational correction to $α$ at the CCSDT/d-aug-cc-pVTZ level and find that this correction contributes approximately 3% to the isotropic ($α$ iso ) and nearly 50% to the anisotropic ($α$ aniso ) polarizability values. In doing so, we find that $α$ iso = 9.8307 bohr 3 , which is in excellent agreement with the experimental value of 9.83 ± 0.02 bohr 3 provided by Russell and Spackman. The DPS reported herein provides a benchmark-quality quantum mechanical estimate of this fundamental quantity of interest and should find extensive use in the development (and assessment) of next-generation force fields and machine-learning based approaches for modeling water in complex condensed-phase environments.

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

Discriminative versus generative approaches to simulation-based inference

Most of the fundamental, emergent, and phenomenological parameters of particle and nuclear physics are determined through parametric template fits. Simulations are used to populate histograms which are then matched to data. This approach is inherently lossy, since histograms are binned and low-dimensional. Deep learning has enabled unbinned and high-dimensional parameter estimation through neural likelihood(-ratio) estimation. We compare two approaches for neural simulation-based inference (NSBI): one based on discriminative learning (classification) and one based on generative modeling. These two approaches are directly evaluated on the same datasets, with a similar level of hyperparameter optimization in both cases. In addition to a Gaussian dataset, we study NSBI using a Higgs boson dataset from the FAIR Universe Challenge. We find that both the direct likelihood and likelihood ratio estimation are able to effectively extract parameters with reasonable uncertainties. For the numerical examples and within the set of hyperparameters studied, we found that the likelihood ratio method is more accurate and/or precise. Both methods have a significant spread from the network training and would require ensembling or other mitigation strategies in practice.

high energy physics↗
DOE OSTI2025

CSGL: chemical synthesis graph learning for molecule representation

Abstract Motivation Molecule representation learning (MRL) translates molecules into a real vector space, serving as input to downstream tasks in biology, chemistry, and computer science. This article introduces a chemical synthesis graph learning (CSGL) framework, which enhances MRL by considering both the atomic structures of molecules and their roles in chemical reactions through a hierarchical graph representation. Specifically, molecules are first modeled based on their molecular graphs, which capture atomic-level structural information. They are then further refined using a chemical synthesis graph, where nodes represent reactant and product molecule sets, and edges encode chemical transformations between reactants and products (e.g. changes in molecular structures). CSGL optimizes molecular embeddings of reactant and product nodes in a fashion that ensures the embeddings conform to a chemical balance constraint. Results Experimental results show that our method CSGL achieves strong performance on a variety of tasks, including product prediction, reaction classification, and molecular property prediction. Availability and implementation https://github.com/li-2023/CSGL.

Biochemistry & Molecular Biology↗
DOE OSTI2020

The Architecture of Metabolism Maximizes Biosynthetic Diversity in the Largest Class of Fungi

Abstract Ecological diversity in fungi is largely defined by metabolic traits, including the ability to produce secondary or “specialized” metabolites (SMs) that mediate interactions with other organisms. Fungal SM pathways are frequently encoded in biosynthetic gene clusters (BGCs), which facilitate the identification and characterization of metabolic pathways. Variation in BGC composition reflects the diversity of their SM products. Recent studies have documented surprising diversity of BGC repertoires among isolates of the same fungal species, yet little is known about how this population-level variation is inherited across macroevolutionary timescales. Here, we applied a novel linkage-based algorithm to reveal previously unexplored dimensions of diversity in BGC composition, distribution, and repertoire across 101 species of Dothideomycetes, which are considered the most phylogenetically diverse class of fungi and known to produce many SMs. We predicted both complementary and overlapping sets of clustered genes compared with existing methods and identified novel gene pairs that associate with known secondary metabolite genes. We found that variation among sets of BGCs in individual genomes is due to nonoverlapping BGC combinations and that several BGCs have biased ecological distributions, consistent with niche-specific selection. We observed that total BGC diversity scales linearly with increasing repertoire size, suggesting that secondary metabolites have little structural redundancy in individual fungi. We project that there is substantial unsampled BGC diversity across specific families of Dothideomycetes, which will provide a roadmap for future sampling efforts. Our approach and findings lend new insight into how BGC diversity is generated and maintained across an entire fungal taxonomic class.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Limits on Stellar-Mass Compact Objects as Dark Matter from Gravitational Lensing of Type Ia Supernovae

We present that the nature of dark matter (DM) remains unknown despite very precise knowledge of its abundance in the Universe. An alternative to new elementary particles postulates DM as made of macroscopic compact halo objects (MACHO) such as black holes formed in the very early Universe. Stellar-mass primordial black holes (PBHs) are subject to less robust constraints than other mass ranges and might be connected to gravitational-wave signals detected by the Laser Interferometer Gravitational-Wave Observatory (LIGO). New methods are therefore necessary to constrain the viability of compact objects as a DM candidate. Here we report bounds on the abundance of compact objects from gravitational lensing of type Ia supernovae (SNe). Current SNe data sets constrain compact objects to represent less than 35.2% (Joint Lightcurve Analysis) and 37.2% (Union 2.1) of the total matter content in the Universe, at 95% confidence level. The results are valid for masses larger than ~ 0.01 M ⊙ (solar masses), limited by the size SNe relative to the lens Einstein radius. We demonstrate the mass range of the constraints by computing magnification probabilities for realistic SNe sizes and different values of the PBH mass. Our bounds are sensitive to the total abundance of compact objects with M ≳ 0.01 M ⊙ and complementary to other observational tests. These results are robust against cosmological parameters, outlier rejection, correlated noise, and selection bias. PBHs and other MACHOs are therefore ruled out as the dominant form of DM for objects associated to LIGO gravitational wave detections. Lastly, these bounds constrain early-Universe models that predict stellar-mass PBH production and strengthen the case for lighter forms of DM, including new elementary particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Model-ready benchmarks for NPP, ANPP, litter fluxes, and recruitment into the 1 cm dbh size class

The intended use of this dataset is to serve as an observational benchmark to evaluate model predictions of NPP, ANPP, litter fluxes, and recruitment at Barro Colorado Island, Panama. This dataset contains four CSV files and one text file. “Benchmarks-NPP-ANPP-R-L.csv” provides estimates of annual ecosystem-level reproductive litter flux (R), leaf litter flux (L), aboveground net primary productivity (ANPP), and net primary productivity (NPP) for 61 field plots throughout tropical, temperate, and boreal forest biomes. An additional 499 plots (n = 550) include estimates of just R, L, and R/L. Each row reports a distinct set of estimates for one sampling interval at one plot. “Metadata-Benchmarks-NPP-ANPP-R-L.csv” contains field descriptions for all data fields in “Benchmarks-NPP-ANPP-R-L.csv”. “References-Benchmarks-NPP-ANPP-R-L.txt” contains full references to the original studies used to produce the observations at each plot included in the data. “Benchmarks-Recruitment.csv” provides estimates of species-level recruitment rates into the 1 cm dbh size class at four CTFS-ForestGeo sites using methods that account for unobserved mortality of new recruits between census intervals (Kohyama et al., 2018). “Metadata-Benchmarks-Recruitment.csv” contains field descriptions for all data fields in “Benchmarks-Recruitment.csv”.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Ultrasonic Measurements of Temperature Profile and Heat Fluxes in Coal-Fired Power Plants (Final Report)

Many industrial processes are inaccessible or inhospitable to characterization by traditional temperature measurement methods, such as thermocouples, especially over prolonged exposure to harsh environments. Ultrasound is an established characterization technology with diverse applications ranging from medical imaging to therapies to flaw detection to nondestructive evaluation. Ultrasound may characterize solid materials and components noninvasively as a nondestructive evaluation modality and obtain internal measurements of material properties. For example, the speed of ultrasound propagation changes with Young’s modulus and Poisson’s ratio, which can be found from its measurements. Traditional ultrasonic characterization assumes all material properties remain constant with the position. When this assumption holds, a property of interest may be measured by relating it to the speed of ultrasound propagation (or a speed of sound, SOS) and measuring the SOS by timing the ultrasound propagation through a known distance. However, when a property of interest is spatially distributed, the propagation time depends on the SOS changing with the position along the ultrasound propagation path. The multiple temperature distributions may lead to an identical time of flight (TOF). Temperature is one property that impacts the speed of ultrasound and often cannot be assumed to remain constant with the position. Previously, in the context of temperature, we addressed the challenge of ultrasonic characterization of spatially distributed properties by developing a method for measuring segmental temperature distributions (MSTD). This method divides the ultrasonic propagation into segments bound by echogenic features. These features provide ultrasonic interfaces where some energy is reflected toward the receiving transducer, and the rest continues through the medium. The time-of-flight between the echoes reflected from echogenic features characterizes the spatial distribution in the properties of interest in the corresponding segment of the ultrasonic propagation path. This project demonstrated the application of the MSTD method in industrial conditions of the coal-fired power plant. We implemented the MSTD using metals and alloys waveguides, which may be the existing structure for which the temperature distribution is characterized or purposefully designed waveguides added to the structure by welding or other means specifically to quantify thermal properties using the MSTD method. Previous iterations of the MSTD method used ceramic and cementitious waveguides, which significantly attenuate ultrasound. On the other hand, low attenuation in metallic waveguides creates interactions between echogenic features which compilates the signal analysis in the segmental TOF measurements. We have established the WG design principles that minimize the interferences between trailing and primary echoes and, in some cases, eliminate them. The waveguides in which echoes do not interfere improve the timing accuracy and the robustness of ultrasonic measurements of the spatial distributions in material properties. Our emphasis remained on the estimation of the temperature distributions. We have developed general recommendations for designing ultrasonically segmented waveguides with the reduced influence of trailing echoes. Two of our waveguide designs were tested in the industry. The first waveguide was designed for insertion into a combustion zone of the utility-scale coal-fired power plant boiler. The second design allows the characterization of temperature distribution in the direction normal to the boiler’s water wall, a large heat exchanger converting the chemical energy released during combustion to the steam driving the electrical power generation turbines. These waveguides were designed to operate within a restrictive space of thermally insulated water wall and incorporate densely located echogenic features while combatting the influence of trailing echoes. The project has successfully demonstrated the feasibility of using the developed method for accurate, continuous, and robust temperature measurements in extreme environments of power generation and other industrial processes. It, therefore, has achieved its overarching goal of advancing the technology readiness level of the novel Ultrasound Measurements of Segmental Temperature Distribution (US-MSTD) method for real-time measurements of the temperature distribution and heat fluxes closer to commercial availability, developing a prototype multipoint measurement system, and validating its performance on coal-fired utility boilers. The success of this project was achieved in collaboration with the power generator, Rocky Mountain Power, and set the stage for the transfer of this technology from the laboratory to the industry.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2026

Estimating Flexibility Envelopes for Residential Customers From Utility Smart Meter Data: Preprint

Demand response from residential customers has significant potential to support power system operations, but accurate flexibility estimation is challenging due to the limited resolution of advanced metering infrastructure (AMI) data. Most utility AMI measurements are recorded at hourly intervals, with only a small portion at higher resolutions, and even fewer households have appliance-level energy usage data. To address this issue, this paper proposes a two-stage long short-term memory (LSTM) framework for estimating household flexibility envelopes from low-resolution AMI data. In the first stage, the heating, ventilating, and air-conditioning (HVAC) load and non-HVAC loads are estimated by using a model trained on a small set of households with appliance-level profiles. These estimated data are then used to compute the upper- and lower-flexibility bounds, which are subsequently down-sampled to lower-resolution data. In the second stage, these flexibility bounds serve as training inputs for another LSTM model, enabling direct prediction of flexibility envelopes for households with only hourly AMI data. This method is validated using Pecan Street data from two different areas, and the results demonstrate its applicability and effectiveness.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2018

An overview of methods to identify and manage uncertainty for modelling problems in the water-environment-agriculture cross-sector

Uncertainty pervades the representation of systems in the water–environment–agriculture cross-sector. Successful methods to address uncertainties have largely focused on standard mathematical formulations of biophysical processes in a single sector, such as partial or ordinary differential equations. More attention to integrated models of such systems is warranted. Model components representing the different sectors of an integrated model can have less standard, and different, formulations to one another, as well as different levels of epistemic knowledge and data informativeness. Thus, uncertainty is not only pervasive but also crosses boundaries and propagates between system components. Uncertainty assessment (UA) cries out for more eclectic treatment in these circumstances, some of it being more qualitative and empirical. Here in this paper, we discuss the various sources of uncertainty in such a cross-sectoral setting and ways to assess and manage them. We have outlined a fast-growing set of methodologies, particularly in the computational mathematics literature on uncertainty quantification (UQ), that seem highly pertinent for uncertainty assessment. There appears to be considerable scope for advancing UA by integrating relevant UQ techniques into cross-sectoral problem applications. Of course this will entail considerable collaboration between domain specialists who often take first ownership of the problem and computational methods experts.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Robotic Solutions for Remote Access and Retrieval of High-Level Radioactive Waste - 19544

This paper describes the ongoing development of robotic solutions for accessing and retrieving high-level radioactive waste from underground storage bins as part of the Calcine Retrieval Project at the INL Site. The waste is in the form of a solid powder, referred to as calcine, stored in a set of stainless steel bins within an underground concrete vault. Due to the radiation field produced by the calcine, all operations must be performed remotely. The Calcine Retrieval Project is currently carrying out proof-of-concept testing of methods for the various stages of calcine retrieval. The bins were designed with no means for access, so retrieval operations will require welding access riser pipes to the bins, followed by cutting through the bin walls from within the risers. These tasks are complicated by the fact that four of the risers must be placed off center on a domed surface, resulting in a weld joint that is curved and tilted relative to the riser axis. The Access Riser Connection robot is being developed to perform these operations using interchangeable welding and plasma cutting attachments. The majority of the calcine will be removed using a vacuum system and pneumatically transferred to storage. This system is not described in this paper. However, the bins contain several internal features that are expected to retain calcine, most notably circumferential stiffening rings on the walls and radial ribs on the floors. One of the solutions being developed for the task of retrieving this residual calcine is a vacuum crawler, a wheeled robot using vacuum to adhere to the interior walls of the bins. The vacuum crawler would use a separate vacuum system to remove residual calcine or would use compressed air to move calcine to a location from which it can be retrieved more easily. While comparable welding and wall-climbing robots have been developed before, this particular application has several constraints that are not satisfied by existing systems. For this reason, both the access riser welder and the vacuum crawler, although still under development, may be of interest to other projects faced with similar constraints. (author)

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

Opinion: Making the Most of Colour Rendition Research

After decades of debate, new methods for characterizing color rendition have been formalized, providing an improved system of tools to aid technology developers and lighting specifiers. Solid-state lighting, previously burdened by legacy measures based on outdated science, is positioned to take advantage of these new tools to provide tailored, efficient lighting that better meets the needs of specific architectural lighting applications. This transition requires new guidance, however, presenting an opportunity for researchers. There are several key methodological considerations specific to experiments investigating color rendition that will help build a base of guidance that can be widely applied. The importance of these ideas mostly stems from new knowledge gained in the last five years, but their origins date back much further. 1.The spectral power distribution must be adequately varied to stress the measure(s) of color rendition being evaluated. Often, this means varying the gamut shape—the pattern of shift across different hues—to probe weaknesses of hue-averaged characterizations, avoiding misleading conclusions resulting from limited variation in the stimuli. 2.Evaluating the role of chromaticity on color quality evaluations requires careful accommodation of chromatic adaptation to prevent a centering bias. Asking for delayed responses may be insufficient if chromaticity changes are apparent to participants, just as side-by-side viewing of different chromaticities may not produce results that can be reliably applied in the field. It is also critical that color rendition and chromaticity conditions are varied independently. 3.The stimulus being viewed is the interaction of light and objects’ spectral reflectances. This interaction should be quantified in order to assess how derived guidelines will apply measures using standardized color sample sets. Although light source characterizations are typically used to model responses, the objects and their presentation may strongly influence evaluations. 4.Context is important. Many factors—some investigated, others hypothesized—can influence subjective evaluations of color quality, including illuminance level, type of scene, object familiarity, color composition, viewing duration, and viewer traits. Although their influence on research results is unavoidable, the effects of these contextual factors must be carefully considered in developing guidance. Appropriate experimental design can often mitigate negative consequences arising from these considerations. At a minimum, it is crucial to understand how each of these considerations can influence experimental outcomes, as well as the interpretation and application of the results in practice. Consulting other experts prior to conducting an experiment can help ensure appropriate, up-to-date research methods will be used, so that the results will be novel and useful—and not rejected during peer review. There is an interesting feedback loop: technology development provides new research tools, which allow for refined research methods and new results that lead to more technology development. It is time to make a concerted, collective effort to improve color rendition research so that the full potential of solid-state lighting and new characterization methods can be realized, ensuring that high-efficiency, high-utility, high-appeal lighting systems become standard.

Royer, Michael P.↗
DOE OSTI2022

Photovoltaic-assisted heat pump water heater system and method

A solar power system has a photovoltaic module that supplies electrical current derived from sunlight to a heater system of a hot water tank. The heater system uses the electrical energy from the photovoltaic module, and optionally also from a local power grid, to heat water up to a selected maximum stored hot water temperature. Control circuitry senses when the level of electrical power is below a threshold solar-power level, and in response sets the stored hot water temperature for the heater system to a first lower temperature value. The control circuitry also senses when the level of electrical power is above the threshold solar-power level, an in response sets the stored hot water temperature for the heater system to a second elevated temperature value higher than the lower temperature, making advantageous use of solar power when available. A mixing unit adds cold water to set the dispensed water temperature.

14 SOLAR ENERGY↗
DOE OSTI2018

A new hybrid algorithm combining a new chaos optimization approach with gradient descent for high dimensional optimization problems

Chaos optimization algorithm (COA) is a recently developed method for global optimization based on chaos theory. It has the features of easy implementation, short execution time and robust mechanisms of escaping from local optima compared with the existing stochastic searching algorithm. By integrating the gradient descent method with a new COA approach using new strategies including symmetrization and levelling, a novel hybrid method, called SCGD, is proposed in this paper. The main idea is to combine the good ability of exploitation of the gradient descent with the high ability of exploration of COA. The purpose is to avoid being trapped into local optima and to improve convergence in large space and high-dimension optimization problems. The simulation results based on a set of 50 benchmark functions with diverse properties and different dimensions show the efficiency and the abilities of the proposed method in finding the global optima compared with the existing methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Low-lying states and total internal partition sums of CH

The electronic structure and spin-orbit states of the CH radical have been systematically investigated using multi-reference configuration interaction (MRCI) and single-reference coupled-cluster (CC) methods. These calculations were performed in conjunction with large correlation-consistent basis sets of quadruple-, quintuple-, and sextuple-ζ quality. To achieve high accuracy, electronic energies for all states were extrapolated to the complete basis set (CBS) limit, enabling the detailed construction of potential energy curves and determination of reliable spectroscopic constants. Spin-orbit coupling effects were explicitly incorporated, and vibrational energy levels were computed via Numerov analysis. Furthermore, the resulting values exhibit good to excellent agreement with available experimental data. Dipole moment and transition dipole moment curves were evaluated to assess the opacity characteristics of CH, revealing that transitions such as Χ 2 Π (u′′ = 0) → Α 2 Δ (u′ = 0), Χ 2 Π (u′′ = 0) → Β 2 Σ − (u′ = 0), Χ 2 Π (u′′ = 0) → C 2 Σ + (u′ = 0), and Χ 2 Π (u′′ = 0) → D 2 Σ + (u′ = 3) are particularly probable. Finally, the total internal partition function sum (TIPS) of CH was computed over a broad temperature range (10–30,000 K) based on our high-accuracy ab initio results.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

Robust predictions of specialized metabolism genes through machine learning

Plant specialized metabolism (SM) enzymes produce lineage-specific metabolites with important ecological, evolutionary, and biotechnological implications. Using Arabidopsis thaliana as a model, we identified distinguishing characteristics of SM and GM (general metabolism, traditionally referred to as primary metabolism) genes through a detailed study of features including duplication pattern, sequence conservation, transcription, protein domain content, and gene network properties. Analysis of multiple sets of benchmark genes revealed that SM genes tend to be tandemly duplicated, coexpressed with their paralogs, narrowly expressed at lower levels, less conserved, and less well connected in gene networks relative to GM genes. Although the values of each of these features significantly differed between SM and GM genes, any single feature was ineffective at predicting SM from GM genes. Using machine learning methods to integrate all features, a prediction model was established with a true positive rate of 87% and a true negative rate of 71%. In addition, 86% of known SM genes not used to create the machine learning model were predicted. We also demonstrated that the model could be further improved when we distinguished between SM, GM, and junction genes responsible for reactions shared by SM and GM pathways, indicating that topological considerations may further improve the SM prediction model. Application of the prediction model led to the identification of 1,220 A. thaliana genes with previously unknown functions, each assigned a confidence measure called an SM score, providing a global estimate of SM gene content in a plant genome.

data integration↗
DOE OSTI2023

Non-local contribution from small scales in galaxy–galaxy lensing: comparison of mitigation schemes

ABSTRACT Recent cosmological analyses with large-scale structure and weak lensing measurements, usually referred to as 3 × 2pt, had to discard a lot of signal to noise from small scales due to our inability to accurately model non-linearities and baryonic effects. Galaxy–galaxy lensing, or the position–shear correlation between lens and source galaxies, is one of the three two-point correlation functions that are included in such analyses, usually estimated with the mean tangential shear. However, tangential shear measurements at a given angular scale θ or physical scale R carry information from all scales below that, forcing the scale cuts applied in real data to be significantly larger than the scale at which theoretical uncertainties become problematic. Recently, there have been a few independent efforts that aim to mitigate the non-locality of the galaxy–galaxy lensing signal. Here, we perform a comparison of the different methods, including the Y-transformation, the point-mass marginalization methodology, and the annular differential surface density statistic. We do the comparison at the cosmological constraints level in a combined galaxy clustering and galaxy–galaxy lensing analysis. We find that all the estimators yield equivalent cosmological results assuming a simulated Rubin Observatory Legacy Survey of Space and Time (LSST) Year 1 like set-up and also when applied to DES Y3 data. With the LSST Y1 set-up, we find that the mitigation schemes yield ∼1.3 times more constraining S8 results than applying larger scale cuts without using any mitigation scheme.

79 ASTRONOMY AND ASTROPHYSICS↗