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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

MISPR : an open-source package for high-throughput multiscale molecular simulations

Computational tools provide a unique opportunity to study and design optimal materials by enhancing our ability to comprehend the connections between their atomistic structure and functional properties. However, designing materials with tailored functionalities is complicated due to the necessity to integrate various computational-chemistry software (not necessarily compatible with one another), the heterogeneous nature of the generated data, and the need to explore vast chemical and parameter spaces. The latter is especially important to avoid bias in scattered data points-based models and derive statistical trends only accessible by systematic datasets. Here, we introduce a robust high-throughput multi-scale computational infrastructure coined MISPR (Materials Informatics for Structure–Property Relationships) that seamlessly integrates classical molecular dynamics (MD) simulations with density functional theory (DFT). By enabling high-performance data analytics and coupling between different methods and scales, MISPR addresses critical challenges arising from the needs of automated workflow management and data provenance recording. The major features of MISPR include automated DFT and MD simulations, error handling, derivation of molecular and ensemble properties, and creation of output databases that organize results from individual calculations to enable reproducibility and transparency. In this work, we describe fully automated DFT workflows implemented in MISPR to compute various properties such as nuclear magnetic resonance chemical shift, binding energy, bond dissociation energy, and redox potential with support for multiple methods such as electron transfer and proton-coupled electron transfer reactions. The infrastructure also enables the characterization of large-scale ensemble properties by providing MD workflows that calculate a wide range of structural and dynamical properties in liquid solutions. MISPR employs the methodologies of materials informatics to facilitate understanding and prediction of phenomenological structure–property relationships, which are crucial to designing novel optimal materials for numerous scientific applications and engineering technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Structural switching dynamically controls the doubly pseudoknotted Rous sarcoma virus–programmed ribosomal frameshifting element

A hallmark of retrovirus replication is the translation of two different polyproteins from one RNA through programmed –1 frameshifting. This is a mechanism in which the actively translating ribosome is induced to slip in the 5′ direction at a defined codon and then continues translating in the new reading frame. Programmed frameshifting controls the stoichiometry of viral proteins and is therefore under stringent evolutionary selection. Forty years ago, the first frameshifting stimulatory element was discovered in the Rous sarcoma virus. The ~120 nt RNA segment was predicted to contain a pseudoknot, but its 3D structure has remained elusive. Now, we have determined cryoEM and X-ray crystallographic structures of this classic retroviral element, finding that it adopts a butterfly-like double-pseudoknot fold. One “wing” contains a dynamic pyrimidine-rich helix, observed crystallographically in two conformations and in a third conformation via cryoEM. The other wing encompasses the predicted pseudoknot, which interacts with a second unexpected pseudoknot through a toggle residue, A2546. This key purine switches conformations between structural states and tunes the stability of interacting residues in the two wings. We find that its mutation can modulate frameshifting by as much as 50-fold, likely by altering the relative abundance of different structural states in the conformational ensemble of the RNA. Taken together, our structure–function analyses reveal how a dynamic double pseudoknot junction stimulates frameshifting by taking advantage of conformational heterogeneity, supporting a multistate model in which high Shannon entropy enhances frameshifting efficiency.

Science & Technology - Other Topics↗
DOE OSTI2020

Enhanced tensor polarization in solid-state targets

In this work, we report measurements of enhanced tensor polarization on solid-state targets. The results here represent an increase in tensor polarization over that previously achieved in high energy and nuclear scattering experiments that focused on the measurement of tensor polarized observables. Enhancement techniques are used which require RF produced close to the Larmor frequency of the target spins and use selective semi-saturation resulting from two sources of irradiation, microwave for the DNP process and the additional RF used to manipulate the population of the energy levels in the target material. The spin dynamics of the solid target are used to align the spins enhancing the ensemble average to improve the figure of merit of the scattering experiment. Target rotation at an optimized rate can lead to additional enhancement by applying selective semi-saturation in polycrystalline materials that possess a Pake doublet in their NMR signal.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Visual HPC Workflows for the Analysis of System Dynamics Models

Visual analytics supported by high performance computing (HPC) accelerates and enhances the discovery, exploration, and analysis of causal patterns in complex system dynamics (SD) models. We present a suite of visualization-assisted ensemble-based techniques for hypothesis generation and testing, and for sensitivity analysis. By employing HPC to provide parallel, on-demand simulation of SD models, one can “steer” an ensemble of simulated scenarios in real time as one first formulates and then informally tests those hypotheses: this provides rapid feedback for analysts to refine their understanding of the causal relationships emergent from a model. Such understandings can be followed and augmented by rigorous application of statistical methods, namely global variance-based sensitivity analysis, Monte-Carlo filtering, adaptive regional sensitivity analysis, and self-organized maps: here timely computation relies on HPC, while effective presentation emphasizes high-dimensional multivariate data visualization. Immersive visualization in virtual 3D environments provides an excellent adjunct to the traditional 2D graphics typically used for SD models, as it generates an embodied understanding of model behavior and facilitates an active, collaborative critique of model structure and output. Finally, we summarize prospects for HPC-enabled visual analytics applied to SD modeling.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Reduced Models for Chemical Kinetics derived from Parallel Ensemble Simulations of Stirred Reactors

Detailed chemical kinetic models are essential to improve the predictive capability of computational fluid dynamics (CFD) simulations and to impact the development of next generation combustion technologies. However, such detailed kinetic models are beyond the reach of CFD codes due to the computational cost of transporting additional chemical degrees of freedom. Reduced order models for combustion kinetics exploit the presence of lower dimensional manifolds in the state space and require far fewer number of scalars to be transported as part of CFD simulations. In this paper, we describe the software infrastructure developed to allow large ensembles of stirred reactor configurations with detailed chemical kinetics. The canonical reactor simulations are incorporated in a workflow with sparse grid and data reduction methodologies aimed at development of reduced order models. We present results from analyzing the results of the stirred reactor calculations for syngas mixture using principal component analysis (PCA) and neural network methodologies. The results show that nonlinear autoencoder model can surpass PCA in representing the original dataset with far fewer degrees of freedom.

Zhang, Pei↗
DOE OSTI2025

Structure-Based Design of Small-Molecule Inhibitors of Human Interleukin-6

Human Interleukin-6 (hIL-6) is a pro inflammatory cytokine that binds to its receptor, IL-6Rα followed by binding to gp130 and subsequent dimerization to form a hexamer signaling complex. As a critical inflammation mediator, hIL-6 is associated with a diverse range of diseases and monoclonal antibodies in clinical use that either target IL-6Rα or hIL-6 to inhibit signaling. Here, we perform high-throughput structure-based computational screening using ensemble docking for small-molecule antagonists for which the target conformations were taken from 600 ns long molecular dynamics simulations of the apo protein. Prior knowledge of the contact sites from binary complex studies and experimental work was incorporated into the docking studies. The top 20 scoring ligands from the in silico studies after post analysis were subjected to in vitro functional assays. Among these compounds, the ligand with the second-highest calculated binding affinity experimentally showed an ~84% inhibitory effect on IL6-induced STAT3 reporter activity at 10 μM concentration. This finding may pave the way for designing small-molecule inhibitors of hIL-6 of therapeutic significance.

Human Interleukin-6↗
DOE OSTI2022

Influence of Long-Range Forces on the Transition States and Dynamics of NaCl Ion-Pair Dissociation in Water

We study NaCl ion-pair dissociation in a dilute aqueous solution using computer simulations both for the full system with long-range Coulomb interactions and for a well-chosen reference system with short-range intermolecular interactions. Analyzing results using concepts from Local Molecular Field (LMF) theory and the recently proposed AI-based analysis tool “State predictive information bottleneck” (SPIB), we show that the system with short-range interactions can accurately reproduce the transition rate for the dissociation process, the dynamics for moving between the underlying metastable states, and the transition state ensemble. Contributions from long-range interactions can be largely neglected for these processes because long-range forces from the direct interionic Coulomb interactions are almost completely canceled (>90%) by those from solvent interactions over the length scale where the transition takes place. Thus, for this important monovalent ion-pair system, short-range forces alone are able to capture detailed consequences of the collective solvent motion, allowing the use of physically suggestive and computationally efficient short-range models for the dissociation event. Here, we believe that the framework here should be applicable to disentangling mechanisms for more complex processes such as multivalent ion disassociation, where previous work has suggested that long-range contributions may be more important.

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

Nanometer-Scale Acoustic Wave Packets Generated by Stochastic Core-Level Photoionization Events

We demonstrate that the absorption of femtosecond hard x-ray pulses excites quasispherical, high-amplitude, and high-wave-vector coherent acoustic phonon wave packets using an all hard-x-ray pump-probe scattering experiment. The time- and momentum-resolved diffuse scattering signal is consistent with an ensemble of 3D strain wave packets induced by the rapid electron cascade dynamics following photoionization at uncorrelated excitation centers. We quantify key parameters of this process, including the localization size of the stress field and the photon energy conversion efficiency into elastic energy. The parameters are determined by the photoelectron and Auger electron cascade dynamics, as well as the electron-phonon interaction. In particular, we obtain the localization size of the observed strain wave packet to be 1.5 and 2.5 nm for bulk SrTiO 3 and KTaO 3 single crystals, respectively. The results provide crucial information on the mechanism of x-ray energy deposition into matter and shed light on the shortest collective length scales accessible to coherent acoustic phonon generation using x-ray excitation. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2020

Power System Event Classification and Localization Using a Convolutional Neural Network

Detection and timely identification of power system disturbances are essential for situation awareness and reliable electricity grid operation. Because records of actual events in the system are limited, ensemble simulation-based events are needed to provide adequate data for building event-detection models through deep learning; e.g., a convolutional neural network (CNN). An ensemble numerical simulation-based training data set have been generated through dynamic simulations performed on the Polish system with various types of faults in different locations. Such data augmentation is proven to be able to provide adequate data for deep learning. The synchronous generators’ frequency signals are used and encoded into images for developing and evaluating CNN models for classification of fault types and locations. With a time-domain stacked image set as the benchmark, two different time-series encoding approaches, i.e., wavelet decomposition-based frequency-domain stacking and polar coordinate system-based Gramian Angular Field (GAF) stacking, are also adopted to evaluate and compare the CNN model performance and applicability. The various encoding approaches are suitable for different fault types and spatial zonation. With optimized settings of the developed CNN models, the classification and localization accuracies can go beyond 84 and 91%, respectively.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Nonadiabatic dynamics with classical trajectories: The problem of an initial coherent superposition of electronic states

Advances in coherent light sources and development of pump–probe techniques in recent decades have opened the way to study electronic motion in its natural time scale. When an ultrashort laser pulse interacts with a molecular target, a coherent superposition of electronic states is created and the triggered electron dynamics is coupled to the nuclear motion. A natural and computationally efficient choice to simulate this correlated dynamics is a trajectory-based method where the quantum-mechanical electronic evolution is coupled to a classical-like nuclear dynamics. These methods must approximate the initial correlated electron–nuclear state by associating an initial electronic wavefunction to each classical trajectory in the ensemble. Different possibilities exist that reproduce the initial populations of the exact molecular wavefunction when represented in a basis. We show that different choices yield different dynamics and explore the effect of this choice in Ehrenfest, surface hopping, and exact-factorization-based coupled-trajectory schemes in a one-dimensional two-electronic-state model system that can be solved numerically exactly. Furthermore, this work aims to clarify the problems that standard trajectory-based techniques might have when a coherent superposition of electronic states is created to initialize the dynamics, to discuss what properties and observables are affected by different choices of electronic initial conditions and to point out the importance of quantum-momentum-induced electronic transitions in coupled-trajectory schemes.

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

How May the Choice of Downscaling Techniques and Meteorological Reference Observations Affect Future Hydroclimate Projections?

Abstract We present an intercomparison of a suite of high‐resolution downscaled climate projections based on a six‐member General Circulation Model (GCM) ensemble from Coupled Models Intercomparison Project (CMIP6). The CMIP6 GCMs have been downscaled using dynamical and statistical downscaling techniques based on two meteorological reference observations over the conterminous United States. We use the regional climate model, RegCM4, for dynamical downscaling, double bias correction constructed analogs method for statistical downscaling, and Daymet and Livneh datasets as the reference observations for statistical training and bias‐correction. We evaluate the performances of downscaled data in both historical and future periods under the SSP585 scenario. While dynamical downscaling improves the simulation of some performance evaluation indices, it adds an extra bias in others, highlighting the need for statistical correction before its use in impact assessments. Downscaled datasets after bias‐correction compare exceptionally well with observations. However, the choice of downscaling techniques and the underlying reference observations influence the hydroclimate characteristics of downscaled data. For instance, the statistical downscaling generally preserves the GCMs climate change signal but overestimates the frequency of hot extremes. Similarly, simulated future changes are sensitive to the choice of reference observations, particularly for precipitation extremes that exhibit a higher projected increase in the ensembles trained and/or corrected by Daymet than Livneh. Overall, these results demonstrate that multiple factors, including downscaling techniques and reference observations, can substantially influence the outcome of downscaled climate projections and stress the need for a comprehensive understanding of such method‐based uncertainties.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Arrested-motility states in populations of shape-anisotropic active Janus particles

The emergence of large-scale collective phenomena from simple interactions between individual units is a hallmark of active matter systems. Active colloids with alignment-dominated interparticle interactions tend to develop orientational order and form motile coherent states, such as flocks and swarms. Alternatively, a combination of self-propulsion and excluded-volume interactions results in self-trapping and active phase separation into dense clusters. Here, we reveal unconventional arrested-motility states in ensembles of active discoidal particles powered by induced-charge electrophoresis. Combining experiments and computational modeling, we demonstrate that the shape asymmetry of the particles promotes the hydrodynamically assisted formation of active particles’ bound states in a certain range of excitation parameters, ultimately leading to a spontaneous collective state with arrested motility. Unlike the jammed clusters obtained through self-trapping, the arrested-motility phase remains sparse, dynamic, and reconfigurable. The demonstrated mechanism of phase separation seeded by bound state formation in ensembles of oblate active particles is generic and should be applicable to other active colloidal systems.

42 ENGINEERING↗
DOE OSTI2021

Assessing the representation of the Australian carbon cycle in global vegetation models

Australia plays an important role in the global terrestrial carbon cycle on inter-annual timescales. While the Australian continent is included in global assessments of the carbon cycle such as the global carbon budget, the performance of dynamic global vegetation models (DGVMs) over Australia has rarely been evaluated. We assessed simulations of net biome production (NBP) and the carbon stored in vegetation between 1901 to 2018 from 13 DGVMs (TRENDY v8 ensemble). We focused our analysis on Australia's short-term (inter-annual) and long-term (decadal to centennial) terrestrial carbon dynamics. The TRENDY models simulated differing magnitudes of NBP on inter-annual timescales, and these differences resulted in significant differences in long-term vegetation carbon accumulation (–4.7 to 9.5 Pg C). We compared the TRENDY ensemble to several satellite-derived datasets and showed that the spread in the models' simulated carbon storage resulted from varying changes in carbon residence time rather than differences in net carbon uptake. Differences in simulated long-term accumulated NBP between models were mostly due to model responses to land-use change. The DGVMs also simulated different sensitivities to atmospheric carbon dioxide (CO 2 ) concentration, although notably, the models with nutrient cycles did not simulate the smallest NBP response to CO 2 . Our results suggest that a change in the climate forcing did not have a large impact on the carbon cycle on long timescales. However, the inter-annual variability in precipitation drives the year-to-year variability in NBP. We analysed the impact of key modes of climate variability, including the El Niño–Southern Oscillation (ENSO) and the Indian Ocean Dipole (IOD), on NBP. While the DGVMs agreed on sign of the response of NBP to El Niño and La Niña and to positive and negative IOD events, the magnitude of inter-annual variability in NBP differed strongly between models. In addition, we find that differences in the timing of simulated phenology and fire dynamics are associated with differences in simulated or prescribed vegetation cover and process representation. We further find model disagreement in simulated vegetation carbon, phenology, and apparent carbon residence time, indicating that the models have different types and coverage of vegetation across Australia (whether prescribed or emergent). Our study highlights the need to evaluate parameter assumptions and the key processes that drive vegetation dynamics, such as phenology, mortality, and fire, in an Australian context to reduce uncertainty across models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Large Ensemble Particle Filter for Spatial Climate Reconstructions Using a Linear Inverse Model

Proxy records that document the last 2000 years of climate provide evidence for the wide range of the natural climate variability from inter-annual to secular timescales not captured by the short window of recent direct observations. Assessing climate models ability to reproduce such natural variations is crucial to understand climate sensitivity and impacts of future climate change. Paleoclimate data assimilation (PDA) offers a powerful way to extend the short instrumental period by optimally combining the physics described by General Circulation Climate Models (GCMs) with information from available proxy records while taking into account their uncertainties. Here we present a new PDA approach based on a sequential importance resampling (SIR) Particle filter (PF) that uses Linear Inverse Modeling (LIM) as an emulator of several CMIP-class GCMs. We examine in a perfect-model framework the skill of the various LIMs to forecast the dynamics of the surface temperatures and provide spatial field reconstructions over the last millennium in a SIR PF. Our results show that the LIMs allow for skillful ensemble forecasts at 1-year lead-time based on GCMs dynamical knowledge with best prediction in the tropics and the North Atlantic. The PDA further provides a set of physically consistent spatial fields allowing robust uncertainty quantification related to climate models biases and proxy spatial sampling. Our results indicate that the LIM yields dynamical memory improving climate variability reconstructions and support the use of the LIM as a GCM-emulator in real reconstruction to propagate large ensembles of particles at low cost in SIR PF.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Sea Level Rise in the CESM Large Ensemble: The Role of Individual Climate Forcings and Consequences for the Coming Decades

The emergence of a spatial pattern in the externally forced response (FR) of dynamic sea level (DSL) during the altimeter era has recently been demonstrated using climate models but our understanding of its initial emergence, drivers, and implications for the future is poor. Here the anthropogenic forcings of the DSL pattern are explored using the Community Earth System Model Large Ensemble (CESM-LE) and Single-Forcing Large Ensemble, a newly available set of simulations where values of individual forcing agents remain fixed at 1920 levels, allowing for an estimation of their effects. Statistically significant contributions to the DSL FR are identified for greenhouse gases (GHGs) and industrial aerosols (AERs), with particularly strong contributions resulting from AERs in the mid-twentieth century and GHGs in the late twentieth and twenty-first century. Secondary, but important, contributions are identified for biomass burning aerosols in the equatorial Atlantic Ocean in the mid-twentieth century, and for stratospheric ozone in the Southern Ocean during the late twentieth century. Key to understanding regional DSL patterns are ocean heat content and salinity anomalies, which are driven by surface heat and freshwater fluxes, ocean dynamics, and the spatial structure of seawater thermal expansivity. Finally, potential implications for the interpretation of DSL during the satellite era and the longer records from tide gauges are suggested as a topic for future research.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Steering Molecular Dynamics Simulations of Membrane-Associated Proteins with Neutron Reflection Results

Herein, we present a novel method to incorporate structural results from surface-sensitive scattering, such as X-ray or neutron reflectometry, into molecular dynamics simulations. While reflectometry techniques generally provide a means to determine the molecular-scale structures of organized interfacial films, they were recently shown to offer the capability to characterize the structures of protein-membrane complexes supported by a solid substrate. One-dimensional information inherent in the experimental results is used in the form of component volume occupancy (CVO) profiles, which describe the distribution of molecular components within an interfacial architecture, to construct real-space constraints in the form of a biasing potential for the simulation that vanishes when the simulated and experimental profiles agree. This approach improves the correspondence between simulation and experiment, as shown in the re-evaluation of an neutron-reflection-derived structure which was approximated by an independent molecular dynamics simulation in earlier work, and it also leads to faster equilibration of ensemble structures. We further show that time averaging the CVO profile that develops in the simulation while biasing with this approach permits fluctuations about the average that are necessary for conformational exploration of the system. This method is particularly valuable for studies of proteins at interfaces that contain disordered regions since the conformation of such regions is difficult to judge from the analysis of one-dimensional experimental profiles and may take prohibitively long to equilibrate in simulations.

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

Operando studies reveal active Cu nanograins for CO 2 electroreduction

Carbon dioxide electroreduction facilitates the sustainable synthesis of fuels and chemicals. Although Cu enables CO 2 -to-multicarbon product (C 2+ ) conversion, the nature of the active sites under operating conditions remains elusive. Importantly, identifying active sites of high-performance Cu nanocatalysts necessitates nanoscale, time-resolved operando techniques. Here, we present a comprehensive investigation of the structural dynamics during the life cycle of Cu nanocatalysts. A 7 nm Cu nanoparticle ensemble evolves into metallic Cu nanograins during electrolysis before complete oxidation to single-crystal Cu 2 O nanocubes following post-electrolysis air exposure. Operando analytical and four-dimensional electrochemical liquid-cell scanning transmission electron microscopy shows the presence of metallic Cu nanograins under CO 2 reduction conditions. Correlated high-energy-resolution time-resolved X-ray spectroscopy suggests that metallic Cu, rich in nanograin boundaries, supports undercoordinated active sites for C–C coupling. Quantitative structure–activity correlation shows that a higher fraction of metallic Cu nanograins leads to higher C 2+ selectivity. A 7 nm Cu nanoparticle ensemble, with a unity fraction of active Cu nanograins, exhibits sixfold higher C 2+ selectivity than the 18 nm counterpart with one-third of active Cu nanograins. Importantly, the correlation of multimodal operando techniques serves as a powerful platform to advance our fundamental understanding of the complex structural evolution of nanocatalysts under electrochemical conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗