Search NASA⌕ Search

SEARCH · Search NASA

Results for “mixture of Gaussians”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Redshift evolution of the underlying type Ia supernova stretch distribution

The detailed nature of type Ia supernovae (SNe Ia) remains uncertain, and as survey statistics increase, the question of astrophysical systematic uncertainties arises, notably that of the evolution of SN Ia populations. We study the dependence on redshift of the SN Ia SALT2.4 light-curve stretch, which is a purely intrinsic SN property, to probe its potential redshift drift. The SN stretch has been shown to be strongly correlated with the SN environment, notably with stellar age tracers. We modeled the underlying stretch distribution as a function of redshift, using the evolution of the fraction of young and old SNe Ia as predicted using the SNfactory dataset, and assuming a constant underlying stretch distribution for each age population consisting of Gaussian mixtures. We tested our prediction against published samples that were cut to have marginal magnitude selection effects, so that any observed change is indeed astrophysical and not observational in origin. In this first study, there are indications that the underlying SN Ia stretch distribution evolves as a function of redshift, and that the age drifting model is a better description of the data than any time-constant model, including the sample-based asymmetric distributions that are often used to correct Malmquist bias at a significance higher than 5σ. The favored underlying stretch model is a bimodal one, composed of a high-stretch mode shared by both young and old environments, and a low-stretch mode that is exclusive to old environments. The precise effect of the redshift evolution of the intrinsic properties of a SN Ia population on cosmology remains to be studied. The astrophysical drift of the SN stretch distribution does affect current Malmquist bias corrections, however, and thereby the distances that are derived based on SN that are affected by observational selection effects. We highlight that this bias will increase with surveys covering increasingly larger redshift ranges, which is particularly important for the Large Synoptic Survey Telescope.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Blanco DECam Bulge Survey (BDBS)

The Blanco DECam Bulge Survey (BDBS) has imaged more than 200 square degrees of the southern Galactic bulge, providing photometry in the ugrizy filters for ~250 million unique stars. The presence of a strong foreground disk population, along with complex reddening and extreme image crowding, has made it difficult to constrain the presence of young and intermediate age stars in the bulge population. We employed an accurate cross-match of BDBS with the latest data release (EDR3) from the Gaia mission, matching more than 140 million sources with BDBS photometry and Gaia EDR3 photometry and astrometry. We relied on Gaia EDR3 astrometry, without any photometric selection, to produce clean BDBS bulge colour-magnitude diagrams (CMDs). Gaia parallaxes were used to filter out bright foreground sources, and a Gaussian mixture model fit to Galactic proper motions could identify stars kinematically consistent with bulge membership. We applied this method to 127 different bulge fields of 1 deg 2 each, with |$\ell$| ≤ 9.5 ° and –9.5 ° ≤ b ≤ –2.5 °. The astrometric cleaning procedure removes the majority of blue stars in each field, especially near the Galactic plane, where the ratio of blue to red stars is .10%, increasing to values ~20% at higher Galactic latitudes. We rule out the presence of a widespread population of stars younger than 2 Gyr. The vast majority of blue stars brighter than the turnoff belong to the foreground population, according to their measured astrometry. We introduce the distance between the observed red giant branch bump and the red clump as a simple age proxy for the dominant population in the field, and we confirm the picture of a predominantly old bulge. Further work is needed to apply the method to estimate ages to fields at higher latitudes, and to model the complex morphology of the Galactic bulge. We also produce transverse kinematic maps, recovering expected patterns related to the presence of the bar and of the X-shaped nature of the bulge.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Time-dependent saturation and physics-based nonlinear model of cross-beam energy transfer

The nonlinear physics of cross-beam energy transfer (CBET) for multi-speckled laser beams is examined using large-scale particle-in-cell simulations for a range of laser and plasma conditions relevant to indirect-drive inertial confinement fusion (ICF) experiments. The time-dependent growth and saturation of CBET involve complex, nonlinear ion and electron dynamics, including ion trapping-induced enhancement and detuning, ion acoustic wave (IAW) nonlinearity, oblique forward stimulated Raman scattering (FSRS), and backward stimulated Brillouin scattering (BSBS) in a CBET-amplified seed beam. Ion-trapping-induced detuning of CBET is captured in the kinetic linear response by a new δf-Gaussian-mixture algorithm, enabling an accurate characterization of trapping-induced non-Maxwellian distributions. Ion trapping induces nonlinear processes, such as changes to the IAW dispersion and nonlinearities (e.g., bowing and self-focusing), which, together with pump depletion, FSRS, and BSBS, determine the time-dependent nature and level of CBET gain as the system approaches a steady state. Using VPIC simulations at intensities at and above the onset threshold for ion trapping and the insight from the time-dependent saturation analyses, we construct a nonlinear CBET model from local laser and plasma conditions that predicts the CBET gain and the energy deposition into the plasma. This model is intended to provide a more accurate, physics-based description of CBET saturation over a wide range of conditions encountered in ICF hohlraums compared with linear CBET gain models with ad hoc saturation clamps often used in laser ray-based methods in multi-physics codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Nonlinear models for coupling the effects of stimulated Raman scattering to inertial confinement fusion codes

Laser plasma instabilities (LPI) reduce driver-target coupling, alter implosion symmetry, and therefore can fundamentally limit fusion performance in inertial confinement fusion (ICF). Developing a predictive modeling capability for LPI effects can critically advance the success of the field. We perform vector particle-in-cell simulations of multi-speckled laser beams undergoing stimulated Raman scattering (SRS) at various densities and intensities relevant to mainly indirectly driven and a subset of parameter space for directly driven ICF systems, focusing on the regimes with intensities above the onset of electron trapping. Based on the wavenumber of the SRS daughter electron plasma wave, we identify several regions with underpinning SRS saturation physics: the electron-trapping dominated region with intermediate kλD values, the strong Landau damping region at larger kλD values, and the region where the Langmuir decay instability arises at lower kλD values. We develop a nonlinear SRS reflectivity model that features the base trapping-dominated scaling of (kλD)−4 and its modifications. Electron trapping deforms the initialized electron distribution functions, and we have developed a new δf-Gaussian-mixture algorithm for an accurate characterization of the trapped hot electron population. With this SRS hot electron description, we construct a nonlinear energy deposition model and a hot electron source model—based on a modified Manley–Rowe relation—suitable for including SRS effects as a sub-grid module in a high-fidelity ICF design code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A comparison of probabilistic generative frameworks for molecular simulations

Generative artificial intelligence is now a widely used tool in molecular science. Despite the popularity of probabilistic generative models, numerical experiments benchmarking their performance on molecular data are lacking. Here, in this work, we introduce and explain several classes of generative models, broadly sorted into two categories: flow-based models and diffusion models. We select three representative models: neural spline flows, conditional flow matching, and denoising diffusion probabilistic models, and examine their accuracy, computational cost, and generation speed across datasets with tunable dimensionality, complexity, and modal asymmetry. Our findings are varied, with no one framework being the best for all purposes. In a nutshell, (i) neural spline flows do best at capturing mode asymmetry present in low-dimensional data, (ii) conditional flow matching outperforms other models for high-dimensional data with low complexity, and (iii) denoising diffusion probabilistic models appear the best for low-dimensional data with high complexity. Our datasets include a Gaussian mixture model and the dihedral torsion angle distribution of the Aib9 peptide, generated via a molecular dynamics simulation. We hope our taxonomy of probabilistic generative frameworks and numerical results may guide model selection for a wide range of molecular tasks.

Artificial intelligence↗
DOE OSTI2024

Secondary structure determines electron transport in peptides

Proteins play a key role in biological electron transport, but the structure–function relationships governing the electronic properties of peptides are not fully understood. Despite recent progress, understanding the link between peptide conformational flexibility, hierarchical structures, and electron transport pathways has been challenging. Here, we use single-molecule experiments, molecular dynamics (MD) simulations, nonequilibrium Green’s function-density functional theory (NEGF-DFT), and unsupervised machine learning to understand the role of secondary structure on electron transport in peptides. Our results reveal a two-state molecular conductance behavior for peptides across several different amino acid sequences. MD simulations and Gaussian mixture modeling are used to show that this two-state molecular conductance behavior arises due to the conformational flexibility of peptide backbones, with a high-conductance state arising due to a more defined secondary structure (beta turn or 3 10 helices) and a low-conductance state occurring for extended peptide structures. These results highlight the importance of helical conformations on electron transport in peptides. Conformer selection for the peptide structures is rationalized using principal component analysis of intramolecular hydrogen bonding distances along peptide backbones. Molecular conformations from MD simulations are used to model charge transport in NEGF-DFT calculations, and the results are in reasonable qualitative agreement with experiments. Projected density of states calculations and molecular orbital visualizations are further used to understand the role of amino acid side chains on transport. Overall, our results show that secondary structure plays a key role in electron transport in peptides, which provides broad avenues for understanding the electronic properties of proteins.

Science & Technology - Other Topics↗
DOE OSTI2023

Automated characterization of spatial and dynamical heterogeneity in supercooled liquids via implementation of machine learning

Abstract A computational approach by an implementation of the principle component analysis (PCA) with K -means and Gaussian mixture (GM) clustering methods from machine learning algorithms to identify structural and dynamical heterogeneities of supercooled liquids is developed. In this method, a collection of the average weighted coordination numbers ( W C N s ‾ ) of particles calculated from particles’ positions are used as an order parameter to build a low-dimensional representation of feature (structural) space for K -means clustering to sort the particles in the system into few meso-states using PCA. Nano-domains or aggregated clusters are also formed in configurational (real) space from a direct mapping using associated meso-states’ particle identities with some misclassified interfacial particles. These classification uncertainties can be improved by a co-learning strategy which utilizes the probabilistic GM clustering and the information transfer between the structural space and configurational space iteratively until convergence. A final classification of meso-states in structural space and domains in configurational space are stable over long times and measured to have dynamical heterogeneities. Armed with such a classification protocol, various studies over the thermodynamic and dynamical properties of these domains indicate that the observed heterogeneity is the result of liquid–liquid phase separation after quenching to a supercooled state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Tracking atomic structure evolution during directed electron beam induced Si-atom motion in graphene via deep machine learning

Using electron beam manipulation, we enable deterministic motion of individual Si atoms in graphene along predefined trajectories. Structural evolution during the dopant motion was explored, providing information on changes of the Si atom neighborhood during atomic motion and providing statistical information of possible defect configurations. The combination of a Gaussian mixture model and principal component analysis applied to the deep learning-processed experimental data allowed disentangling of the atomic distortions for two different graphene sublattices. This approach demonstrates the potential of e-beam manipulation to create defect libraries of multiple realizations of the same defect and explore the potential of symmetry breaking physics. The rapid image analytics enabled via a deep learning network further empowers instrumentation for e-beam controlled atom-by-atom fabrication. Here, the analysis described in the paper can be reproduced via an interactive Jupyter notebook at https://git.io/JJ3Bx.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A robust estimator of mutual information for deep learning interpretability

Abstract We develop the use of mutual information (MI), a well-established metric in information theory, to interpret the inner workings of deep learning (DL) models. To accurately estimate MI from a finite number of samples, we present GMM-MI (pronounced ‘Jimmie’), an algorithm based on Gaussian mixture models that can be applied to both discrete and continuous settings. GMM-MI is computationally efficient, robust to the choice of hyperparameters and provides the uncertainty on the MI estimate due to the finite sample size. We extensively validate GMM-MI on toy data for which the ground truth MI is known, comparing its performance against established MI estimators. We then demonstrate the use of our MI estimator in the context of representation learning, working with synthetic data and physical datasets describing highly non-linear processes. We train DL models to encode high-dimensional data within a meaningful compressed (latent) representation, and use GMM-MI to quantify both the level of disentanglement between the latent variables, and their association with relevant physical quantities, thus unlocking the interpretability of the latent representation. We make GMM-MI publicly available in this GitHub repository.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

End-to-End Automated Segmentation Framework for Four-Dimensional Scanning Transmission Electron Microscopy Data

Four-dimensional scanning transmission electron microscopy (4D-STEM) is powerful for rapidly characterizing arrays of nanoparticles produced via high-throughput synthesis. However, such 4D-STEM datasets typically contain thousands of nanoparticles, each characterized by thousands of diffraction patterns spatially distributed across the nanoparticle, necessitating efficient and comprehensive analysis. We propose an end-to-end segmentation framework to automatically segment each nanoparticle into regions with distinct composition/orientation of crystal grains, using only the 4D-STEM data. Bragg disk information is extracted in a physics-informed manner from the diffraction patterns at each spatial location and combined with the real space coordinates to form feature vectors. These feature vectors are then used as inputs to a Gaussian mixture model (GMM) to segment the nanoparticle into distinct regions. We also develop two visualization tools based on the GMM outputs to infer the interface transition and the degree of superposition. Our framework comprehensively integrates machine learning tools and physics knowledge, and provides a basis for substantially compressing enormous 4D-STEM datasets, e.g., by replacing the full 4D-STEM dataset for each nanoparticle with only a single set of Bragg disk features for each distinct crystal grain identified in the nanoparticle. In this article, we demonstrate the power of our framework by presenting results for real, complex datasets.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Blanco DECam Bulge Survey (BDBS) II: project performance, data analysis, and early science results

ABSTRACT The Blanco DECam Bulge Survey (BDBS) imaged more than 200 sq deg of the Southern Galactic bulge using the ugrizY filters of the Dark Energy Camera, and produced point spread function photometry of approximately 250 million unique sources. In this paper, we present details regarding the construction and collation of survey catalogues, and also discuss the adopted calibration and dereddening procedures. Early science results are presented with a particular emphasis on the bulge metallicity distribution function and globular clusters. A key result is the strong correlation (σ ∼ 0.2 dex) between (u − i)o and [Fe/H] for bulge red clump giants. We utilized this relation to find that interior bulge fields may be well described by simple closed box enrichment models, but fields exterior to b ∼ −6° seem to require a secondary metal-poor component. Applying scaled versions of the closed box model to the outer bulge fields is shown to significantly reduce the strengths of any additional metal-poor components when compared to Gaussian mixture models. Additional results include: a confirmation that the u band splits the subgiant branch in M22 as a function of metallicity, the detection of possible extratidal stars along the orbits of M 22 and FSR 1758, and additional evidence that NGC 6569 may have a small but discrete He spread, as evidenced by red clump luminosity variations in the reddest bands. We do not confirm previous claims that FSR 1758 is part of a larger extended structure.

Johnson, Christian I.↗
DOE OSTI2022

Connecting galaxy evolution in clusters with their radial profiles and phase space distribution: results from the IllustrisTNG hydrodynamical simulations

ABSTRACT We study the population of galaxies around galaxy clusters in the hydrodynamic simulation suite IllustrisTNG 300-1 to study the signatures of their evolutionary history on observable properties. We measure the radial number density profile, phase space distribution, and splashback radius for galaxies of different masses and colours over the redshift range z = 0−1. The three primary physical effects which shape the galaxy distribution within clusters are the galaxy quenching, angular momentum distribution, and dynamical friction. We find three distinct populations of galaxies by applying a Gaussian mixture model to their distribution in colour and mass. They have distinct evolutionary histories and leave distinct signatures on their distribution around cluster haloes. We find that low-mass red galaxies show the most concentrated distribution in clusters and the largest splashback radius, while high-mass red galaxies show a less concentrated distribution and a smaller splashback radius. Blue galaxies, which mostly quench into the low-mass red population, have the shallowest distribution within the clusters, with those on radial orbits quenched rapidly before reaching pericentre. Comparison with the distribution of galaxies from the Dark Energy Survey survey around Sunyaev–Zeldovich clusters from the Atacama Cosmology Telescope and South Pole Telescope surveys shows evidence for differences in galaxy evolution between simulations and data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

All-purpose, all-sky photometric redshifts for the Legacy Imaging Surveys Data Release 8

ABSTRACT In this paper, we present photometric redshift (photo-z) estimates for the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Surveys, currently the most sensitive optical survey covering the majority of the extragalactic sky. Our photo-z methodology is based on a machine-learning approach, using sparse Gaussian processes augmented with Gaussian mixture models (GMMs) that allow regions of parameter space to be identified and trained separately in a purely data-driven way. The same GMMs are also used to calculate cost-sensitive learning weights that mitigate biases in the spectroscopic training sample. By design, this approach aims to produce reliable and unbiased predictions for all parts of the parameter space present in wide area surveys. Compared to previous literature estimates using the same underlying photometry, our photo-zs are significantly less biased and more accurate at z > 1, with negligible loss in precision or reliability for resolved galaxies at z < 1. Our photo-z estimates offer accurate predictions for rare high-value populations within the parent sample, including optically selected quasars at the highest redshifts (z > 6), as well as X-ray or radio continuum selected populations across a broad range of flux (densities) and redshift. Deriving photo-z estimates for the full Legacy Imaging Surveys Data Release 8, the catalogues provided in this work offer photo-z estimates predicted to be of high quality for ≳9 × 108 galaxies over ${\sim}19\, 400\, \text{deg}^{2}$ and spanning 0 < z ≲ 7, offering one of the most extensive samples of redshift estimates ever produced.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A unified exploration of the chronology of the Galaxy

The Milky Way has distinct structural stellar components linked to its formation and subsequent evolution, but disentangling them is non-trivial. With the recent availability of high-quality data for a large numbers of stars in the Milky Way, it is a natural next step for research in the evolution of the Galaxy to perform automated explorations with unsupervised methods of the structures hidden in the combination of large-scale spectroscopic, astrometric, and asteroseismic data sets. We determine precise stellar properties for 21 076 red giants, mainly spanning 2–15 kpc in Galactocentric radii, making it the largest sample of red giants with measured asteroseismic ages available to date. We explore the nature of different stellar structures in the Galactic disc by using Gaussian mixture models as an unsupervised clustering method to find substructure in the combined chemical, kinematic, and age subspace. The best-fitting mixture model yields four distinct physical Galactic components in the stellar disc: the thin disc, the kinematically heated thin disc, the thick disc, and the stellar halo. We find hints of an age asymmetry between the Northern and Southern hemisphere, and we measure the vertical and radial age gradient of the Galactic disc using the asteroseismic ages extended to further distances than previous studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Artificial-intelligence-assisted analysis of 28 Si * → 7⁢𝛼 breakup data

Mid-weight 𝛼-conjugate nuclei are predicted to possess exotic toroid like resonances with high angular momenta. The search for these states in 28 Si* is the main point of two published experimental investigations of the peripheral 28 Si + 12 C reaction by Cao and collaborators and by Hannaman and collaborators. In this work, we develop a novel artificial intelligence (AI) based machine learning method utilizing the Gaussian Mixture Model (GMM) to analyze available experimental and theoretical data. Here, we additionally study the reaction with the Hybrid 𝛼-Cluster (H⁡𝛼⁢C) model. In all the examined data, our results suggest the presence of underlying structure which is close to that predicted for toroidal states.

Breakup reactions↗
DOE OSTI2023

Attention-Augmented Parametric Kernel Graph Neural Network (APKGNN) for Node Classification

We present a new graph neural network, the Attention-based Parametric-Kernel augmented Graph Neural Network (APKGNN), developed for node classification tasks. Despite extensive work on modeling multi-faceted relationships between connected nodes of a graph, the effect of attention on edge features mapped to relationships has not yet been analyzed through learning representation. This study derives such an attention vector by first calculating node features corresponding to endpoints of an edge and then aggregating these with extracted local intrinsic patches of a given graph to generate augmented local patch vectors. This process uses a parametric kernel based on Gaussian mixture models (GMMs) to embed local neighborhoods of the graph in local patches. The patch vectors then convolve with the above node features to produce an updated node representation. We show that this new learning representation (APKGNN) achieves higher node classification accuracy on tasks - both standard benchmarks (Cora, PubMed, Citeseer) and new experimental short text corpora where nodes correspond to text documents and words. This implementation of the GNN convolution layer outperforms state-of-the-art (SOTA) algorithms, achieving higher training, validation, and test accuracy by a significant margin on three standard benchmark data sets under both SOTA experimental settings and those for new testbeds.

Bose, Avishek↗
DOE OSTI2020

Thermal Overloading Risk Mitigation With a Semi-Analytical Probabilistic Model on Branch Current

A semi-analytical formulation is presented in this paper for the probability computation of branch current in multiphase systems. The developed formula is derived based on the linear power flow model in rectangular coordinates. The system uncertainty injections can be renewable energy resources or loads and are modeled using a Gaussian mixture model (GMM). The developed formula can be used to compute the line current violation probability as well as integrate into optimal power flow problem as chance-constraint relaxation. Here, the proposed formula is first compared with the Matlab embedded numerical integration function to show its performance. Besides, the semi-analytical formula is validated and compared with the Monte Carlo simulation method using the IEEE 123-bus system, EPRI Ckt5, and Ckt7 systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Evaluating proxies for the drivers of natural gas productivity using machine-learning models

We report the extensive development of unconventional reservoirs using horizontal drilling and multistage hydraulic fracturing has generated large volumes of reservoir characterization and production data. The analysis of this abundant data using statistical methods and advanced machine-learning (ML) techniques can provide data-driven insights into well performance. Most predictive modeling studies have focused on the impact that different well completion and stimulation strategies have on well production but have not fully exploited the available in situ rock property data to determine its role in reservoir productivity. We have used machine-learning techniques to rank rock mechanical properties, microseismic attributes, and stimulation parameters in the order of their significance for predicting natural gas production from an unconventional reservoir. The data for this study came from a hydraulically fractured well in the Marcellus Shale in Monongalia County, West Virginia. The data classes included measurements aggregated by well completion stage that included (1) gas production, (2) well-log-derived measurements including bulk density, elastic moduli, shear impedance, compressional impedance, brittleness, and gamma measurements, (3) microseismic attributes, (4) long-period long-duration (LPLD) event counts, (5) fracture counts, and (6) stimulation parameters that included the fluid injection volume and average pumping pressure. To identify observable proxies for the drivers of gas production, we evaluated five commonly used ML approaches including multivariate adaptive regression spline, Gaussian mixture model, random forest, gradient boosting, and neural network. We selected five variables including LPLD event count, seismogenic b-value, hydraulic diffusivity, cumulative moment, and fluid volume as the features most likely to impact gas productivity at the stage level in the study area. The data-driven selection of these parameters for their importance in determining gas production can help reservoir engineers design more effective hydraulic-fracture treatments in the Marcellus Shale and other similar unconventional reservoirs. Plain language summary: We use machine-learning methods and data-driven selection of reservoir parameters to rank and better understand their importance in determining gas production, which can help reservoir engineers design more effective hydraulic-fracture treatments in the Marcellus Shale and other similar unconventional reservoirs.

58 GEOSCIENCES↗