Search NASA⌕ Search

SEARCH · Search NASA

Results for “ensemble dynamics”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Situation awareness and dynamic ensemble forecasting of abnormal behavior in cyber-physical system

A plurality of monitoring nodes may each generate a time-series of current monitoring node values representing current operation of a cyber-physical system. A feature-based forecasting framework may receive the time-series of and generate a set of current feature vectors using feature discovery techniques. The feature behavior for each monitoring node may be characterized in the form of decision boundaries that separate normal and abnormal space based on operating data of the system. A set of ensemble state-space models may be constructed to represent feature evolution in the time-domain, wherein the forecasted outputs from the set of ensemble state-space models comprise anticipated time evolution of features. The framework may then obtain an overall features forecast through dynamic ensemble averaging and compare the overall features forecast to a threshold to generate an estimate associated with at least one feature vector crossing an associated decision boundary.

Abbaszadeh, Masoud↗
DOE OSTI2021

libEnsemble: A Library to Coordinate the Concurrent Evaluation of Dynamic Ensembles of Calculations

Almost all applications stop scaling at some point; those that don't are seldom performant when considering time to solution on anything but aspirational/unicorn resources. Recognizing these tradeoffs as well as greater user functionality in a near-term exascale computing era, we present libEnsemble, a library aimed at particular scalability- and capability-stretching uses. libEnsemble enables running concurrent instances of an application in dynamically allocated ensembles through an extensible Python library. Here, we highlight the structure, execution, and capabilities of the library on leading pre-exascale environments as well as advanced capabilities for exascale environments and beyond.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

libEnsemble: A complete Python toolkit for dynamic ensembles of calculations

Almost all science and engineering applications eventually stop scaling: their runtime no longer decreases as available computational resources increase. Therefore, many applications will struggle to efficiently use emerging extreme-scale high-performance, parallel, and distributed systems. libEnsemble is a complete Python toolkit and workflow system for intelligently driving ensembles of experiments or simulations at massive scales. It enables and encourages multidisciplinary design, decision, and inference studies portably running on laptops, clusters, and supercomputers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Generative Ensemble Regression: Learning Particle Dynamics from Observations of Ensembles with Physics-informed Deep Generative Models

Here, we propose a new method for inferring the governing stochastic ordinary differential equations (SODEs) by observing particle ensembles at discrete and sparse time instants, i.e., multiple “snapshots.” Particle coordinates at a single time instant, possibly noisy or truncated, are recorded in each snapshot but are unpaired across the snapshots. By training a physics-informed generative model that generates “fake” sample paths, we aim to fit the observed particle ensemble distributions with a curve in the probability measure space, which is induced from the inferred particle dynamics. We employ different metrics to quantify the differences between distributions, e.g., the sliced Wasserstein distances and the adversarial losses in generative adversarial networks. We refer to this method as generative “ensemble-regression” (GER), in analogy to the classic “point-regression,” where we infer the dynamics by performing regression in the Euclidean space. We illustrate the GER by learning the drift and diffusion terms of particle ensembles governed by SODEs with Brownian motions and Lévy processes up to 100 dimensions. We also discuss how to treat cases with noisy or truncated observations. Apart from systems consisting of independent particles, we also tackle nonlocal interacting particle systems with unknown interaction potential parameters by constructing a physics-informed loss function. Finally, we investigate scenarios of paired observations and discuss how to reduce the dimensionality in such cases by proving a convergence theorem that provides theoretical support.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Trajectory Ensemble Methods Provide Single-Molecule Statistics for Quantum Dynamical Systems

The emergence of experiments capable of probing quantum dynamics at the single-molecule level requires the development of new theoretical tools capable of simulating and analyzing these dynamics beyond an ensemble-averaged description. In this article, we present an efficient method for sampling and simulating the dynamics of the individual quantum systems that make up an ensemble and apply it to study the nonequilibrium dynamics of the ubiquitous spin-boson model. We generate an ensemble of single-system trajectories, and we analyze this trajectory ensemble using tools from classical statistical mechanics. Our results demonstrate that the dynamics of quantum coherence is highly heterogeneous at the single-system level due to variations in the initial bath configuration, which significantly affects the transient exchange of coherence between the system and its bath. Here, we observe that single systems tend to retain coherence over time scales longer than that of the ensemble. We also compute a novel thermodynamic entanglement entropy that quantifies a thermodynamic driving force favoring system–bath entanglement.

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

Direct observation of tRNA-chaperoned folding of a dynamic mRNA ensemble

T-box riboswitches are multi-domain noncoding RNAs that surveil individual amino acid availabilities in most Gram-positive bacteria. T-boxes directly bind specific tRNAs, query their aminoacylation status to detect starvation, and feedback control the transcription or translation of downstream amino-acid metabolic genes. Most T-boxes rapidly recruit their cognate tRNA ligands through an intricate three-way stem I-stem II-tRNA interaction, whose establishment is not understood. Using single-molecule FRET, SAXS, and time-resolved fluorescence, we find that the free T-box RNA assumes a broad distribution of open, semi-open, and closed conformations that only slowly interconvert. tRNA directly binds all three conformers with distinct kinetics, triggers nearly instantaneous collapses of the open conformations, and returns the T-box RNA to their pre-binding conformations upon dissociation. This scissors-like dynamic behavior is enabled by a hinge-like pseudoknot domain which poises the T-box for rapid tRNA-induced domain closure. This study reveals tRNA-chaperoned folding of flexible, multi-domain mRNAs through a Venus flytrap-like mechanism.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Portable, heterogeneous ensemble workflows at scale using libEnsemble

libEnsemble is a Python-based toolkit for running dynamic ensembles, developed as part of the DOE Exascale Computing Project. The toolkit utilizes a unique generator–simulator–allocator paradigm, where generators produce input for simulators, simulators evaluate those inputs, and allocators decide whether and when a simulator or generator should be called. The generator steers the ensemble based on simulation results. Generators may, for example, apply methods for numerical optimization, machine learning, or statistical calibration. libEnsemble communicates between a manager and workers. Flexibility is provided through multiple manager–worker communication substrates each of which has different benefits. These include Python’s multiprocessing, mpi4py, and TCP. Multisite ensembles are supported using Balsam or Globus Compute. We overview the unique characteristics of libEnsemble as well as current and potential interoperability with other packages in the workflow ecosystem. We highlight libEnsemble’s dynamic resource features: libEnsemble can detect system resources, such as available nodes, cores, and GPUs, and assign these in a portable way. These features allow users to specify the number of processors and GPUs required for each simulation; and resources will be automatically assigned on a wide range of systems, including Frontier, Aurora, and Perlmutter. Such ensembles can include multiple simulation types, some using GPUs and others using only CPUs, sharing nodes for maximum efficiency. We also describe the benefits of libEnsemble’s generator–simulator coupling, which easily exposes to the user the ability to cancel, and portably kill, running simulations based on models that are updated with intermediate simulation output. We demonstrate libEnsemble’s capabilities, scalability, and scientific impact via a Gaussian process surrogate training problem for the longitudinal density profile at the exit of a plasma accelerator stage. In conclusion, the study uses gpCAM for the surrogate model and employs either Wake-T or WarpX simulations, highlighting efficient use of resources that can easily extend to exascale.

Dynamic ensembles↗
DOE OSTI2024

Unlocking the unfolded structure of ubiquitin: Combining time-resolved x-ray solution scattering and molecular dynamics to generate unfolded ensembles

The unfolding dynamics of ubiquitin were studied using a combination of x-ray solution scattering (XSS) and molecular dynamics (MD) simulations. The kinetic analysis of the XSS ubiquitin signals showed that the protein unfolds through a two-state process, independent of the presence of destabilizing salts. In order to characterize the ensemble of unfolded states in atomic detail, the experimental XSS results were used as a constraint in the MD simulations through the incorporation of x-ray scattering derived potential to drive the folded ubiquitin structure toward sampling unfolded states consistent with the XSS signals. We detail how biased MD simulations provide insight into unfolded states that are otherwise difficult to resolve and underscore how experimental XSS data can be combined with MD to efficiently sample structures away from the native state. Our results indicate that ubiquitin samples unfolded in states with a high degree of loss in secondary structure yet without a collapse to a molten globule or fully solvated extended chain. Finally, we propose how using biased-MD can significantly decrease the computational time and resources required to sample experimentally relevant nonequilibrium states.

Chemistry↗
DOE OSTI2025

Early events in G-quadruplex folding captured by time-resolved small-angle X-ray scattering

Abstract Time-resolved small-angle X-ray experiments are reported here that capture and quantify a previously unknown rapid collapse of the unfolded oligonucleotide as an early step in the folding of hybrid 1 and hybrid 2 telomeric G-quadruplex structures. The rapid collapse, initiated by a pH jump, is characterized by an exponential decrease in the radius of gyration from 24.3 to 12.6 Å. The collapse is monophasic and is complete in <600 ms. Additional hand-mixing pH-jump kinetic studies show that slower kinetic steps follow the collapse. The folded and unfolded states at equilibrium were further characterized by SAXS studies and other biophysical tools, showing that G4 unfolding was complete at alkaline pH, but not in LiCl solution as is often claimed. The SAXS Ensemble Optimization Method analysis reveals models of the unfolded state as a dynamic ensemble of flexible oligonucleotide chains with a variety of transient hairpin structures. These results suggest a G4 folding pathway in which a rapid collapse, analogous to molten globule formation seen in proteins, is followed by a confined conformational search within the collapsed particle to form the native contacts ultimately found in the stable folded form.

Biochemistry & Molecular Biology↗
DOE OSTI2023

Causes and evolution of winter polynyas north of Greenland

During the 42-year period (1979–2020) of satellite measurements, four major winter (December–March) polynyas have been observed north of Greenland: one in December 1986 and three in the last decade, i.e., February of 2011, 2017, and 2018. The 2018 polynya was unparalleled in its magnitude and duration compared to the three previous events. Given the apparent recent increase in the occurrence of these extreme events, this study aims to examine their evolution and causality, in terms of forced versus natural variability. The limited weather station and remotely sensed sea ice data are analyzed combining with output from the fully coupled Regional Arctic System Model (RASM), including one hindcast and two ensemble simulations. We found that neither the accompanying anomalous warm surface air intrusion nor the ocean below had an impact (i.e., no significant ice melting) on the evolution of the observed winter open-water episodes in the region. Instead, the extreme atmospheric wind forcing resulted in greater sea ice deformation and transport offshore, accounting for the majority of sea ice loss in all four polynyas. Our analysis suggests that strong southerly winds (i.e., northward wind with speeds greater than 10 m s –1 ) blowing persistently over the study region for at least 2 d or more were required over the study region to mechanically redistribute some of the thickest Arctic sea ice out of the region and thus to create open-water areas (i.e., a latent heat polynya). To assess the role of internal variability versus external forcing of such events, we carried out and examined results from the two RASM ensembles dynamically downscaled with output from the Community Earth System Model (CESM) Decadal Prediction Large Ensemble (DPLE) simulations. Out of 100 winters in each of the two ensembles (initialized 30 years apart: one in December 1985 and another in December 2015), 17 and 16 winter polynyas were produced north of Greenland, respectively. The frequency of polynya occurrence had no apparent sensitivity to the initial sea ice thickness in the study area pointing to internal variability of atmospheric forcing as a dominant cause of winter polynyas north of Greenland. We assert that dynamical downscaling using a high-resolution regional climate model offers a robust tool for process-level examination in space and time, synthesis with limited observations, and probabilistic forecasts of Arctic events, such as the ones being investigated here and elsewhere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Reconfigurable structure and tunable transport in synchronized active spinner materials

Ensembles of actuated colloids are excellent model systems to explore emergent out-of-equilibrium structures, complex collective dynamics, and design rules for the next generation materials. Here, we demonstrate that ferromagnetic microparticles suspended at an air-water interface and energized by an external rotating magnetic field spontaneously form dynamic ensembles of synchronized spinners in a certain range of the excitation field parameters. Each spinner generates strong hydrodynamic flows, and collective interactions of the multiple spinners promote a formation of dynamic lattices. On the basis of experiments and simulations, we reveal structural transitions from liquid to nearly crystalline states in this novel active spinner material and demonstrate that dynamic spinner lattices are reconfigurable, capable of self-healing behavior and that the transport of embedded inert cargo particles can be remotely tuned by the parameters of the external excitation field. Our findings provide insights into the behavior of active spinner materials with reconfigurable structural order and tunable functionalities.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Characterization of Heparin’s Conformational Ensemble by Molecular Dynamics Simulations and Nuclear Magnetic Resonance Spectroscopy

Heparin is a highly charged, polysulfated polysaccharide and serves as an anticoagulant. Heparin binds to multiple proteins throughout the body, suggesting a large range of potential therapeutic applications. Although its function has been characterized in multiple physiological contexts, heparin’s solution conformational dynamics and structure–function relationships are not fully understood. Molecular dynamics (MD) simulations facilitate the analysis of a molecule’s underlying conformational ensemble, which then provides important information necessary for understanding structure–function relationships. However, for MD simulations to afford meaningful results, they must both provide adequate sampling and accurately represent the energy properties of a molecule. The aim of this study is to compare heparin’s conformational ensemble using two well-developed force fields for carbohydrates, known as GLYCAM06 and CHARMM36, using replica exchange molecular dynamics (REMD) simulations, and to validate these results with NMR experiments. In this work, the anticoagulant sequence, an ultra-low-molecular-weight heparin, known as Arixtra (fondaparinux, sodium), was simulated with both parameter sets. The results suggest that GLYCAM06 matches experimental nuclear magnetic resonance three-bond J-coupling values measured for Arixtra better than CHARMM36. In addition, NOESY and ROESY experiments suggest that Arixtra is very flexible in the sub-millisecond time scale and does not adopt a unique structure at 25 C. Moreover, GLYCAM06 affords a much more dynamic conformational ensemble for Arixtra than CHARMM36.

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

Visualizing the strongly reshaped skyrmion Hall effect in multilayer wire devices

Magnetic skyrmions are nanoscale spin textures touted as next-generation computing elements. When subjected to lateral currents, skyrmions move at considerable speeds. Their topological charge results in an additional transverse deflection known as the skyrmion Hall effect (SkHE). While promising, their dynamic phenomenology with current, skyrmion size, geometric effects and disorder remain to be established. Here we report on the ensemble dynamics of individual skyrmions forming dense arrays in Pt/Co/MgO wires by examining over 20,000 instances of motion across currents and fields. The skyrmion speed reaches 24 m/s in the plastic flow regime and is surprisingly robust to positional and size variations. Meanwhile, the SkHE saturates at ~22 ° , is substantially reshaped by the wire edge, and crucially increases weakly with skyrmion size. Particle model simulations suggest that the SkHE size dependence - contrary to analytical predictions - arises from the interplay of intrinsic and pinning-driven effects. These results establish a robust framework to harness SkHE and achieve high-throughput skyrmion motion in wire devices.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Effects of pH on an IDP conformational ensemble explored by molecular dynamics simulation

The conformational ensemble of intrinsically disordered proteins, such as α-synuclein, are responsible for their function and malfunction. Misfolding of α-synuclein can lead to neurodegenerative diseases, and the ability to study their conformations and those of other intrinsically disordered proteins under varying physiological conditions can be crucial to understanding and preventing pathologies. In contrast to well-folded peptides, a consensus feature of IDPs is their low hydropathy and high charge, which makes their conformations sensitive to pH perturbation. We examine a prominent member of this subset of IDPs, α-synuclein, using a divide-and-conquer scheme that provides enhanced sampling of IDP structural ensembles. We constructed conformational ensembles of α-synuclein under neutral (pH ~ 7) and low (pH ~ 3) pH conditions and compared our results with available information obtained from smFRET, SAXS, and NMR studies. Specifically, α-synuclein has been found to in a more compact state at low pH conditions and the structural changes observed are consistent with those from experiments. We also characterize the conformational and dynamic differences between these ensembles and discussed the implication on promoting pathogenic fibril formation. We find that under low pH conditions, neutralization of negatively charged residues leads to compaction of the C-terminal portion of α-synuclein while internal reorganization allows α-synuclein to maintain its overall end-to-end distance. Here, we also observe different levels of intra-protein interaction between three regions of α-synuclein at varying pH and a shift towards more hydrophilic interactions with decreasing pH.

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

A stochastic covariance shrinkage approach to particle rejuvenation in the ensemble transform particle filter

Rejuvenation in particle filters is necessary to prevent the collapse of the weights when the number of particles is insufficient to properly sample the high-probability regions of the state space. Rejuvenation is often implemented in a heuristic manner by the addition of random noise that widens the support of the ensemble. This work aims at improving canonical rejuvenation methodology by the introduction of additional prior information obtained from climatological samples; the dynamical particles used for importance sampling are augmented with samples obtained from stochastic covariance shrinkage. A localized variant of the proposed method is developed. Numerical experiments with the Lorenz '63 model show that modified filters significantly improve the analyses for low dynamical ensemble sizes. Furthermore, localization experiments with the Lorenz '96 model show that the proposed methodology is extendable to larger systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Parton distribution functions from scalar light-front parton gas model

Here, we propose an application of a microcanonical ensemble with light-front kinematics to model the phase-space distribution of relativistic constituents of a bound state. These constituents denoted by partons are treated as classical spin-zero particles confined inside the bound state with inter-parton collisions as their only interaction. The microcanonical molecular dynamics ensemble is applied to obtain the phase-space distribution of such a thermodynamic system. We sample this phase-space distribution using Monte Carlo algorithms to obtain the parton distribution functions (PDFs) in scenarios with 3, 4, and 5 identical partons. In addition PDFs when a selected number of massless partons are mixed with 3 massive partons are also presented.

Microcanonical ensemble↗