Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic ensembles”

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 163 records · Page 9

Dynamic response of a pulsed Burke-Schumann diffusion flame

Turbulent flames are often envisioned as an ensemble of random vortices interacting with the combustion process. A better understanding of the vortex-flame interactions therefore would be useful in improving the modeling of turbulent diffusion flames. Substantial simplification may be made by investigating controlled interactions in a laminar flame, as opposed to random interactions in a turbulent flame. The general goals of the research project are to improve our understanding of (1) the influence of buoyancy on co-flow diffusion flames and (2) the effects of buoyancy on vortex-flame interactions in co-flow diffusion flames. As a first step toward objective (2), we conducted a joint experimental and numerical investigation of the vortex-flame interaction. Vortices were produced by mechanically pulsing the fuel flow at a low frequency, e.g., 10 Hz. Experiments were conducted using a nonflickering Burke-Schumann flame in both microgravity (mu-g) and normal gravity (1g) as a means of varying the buoyant force without modification of the pressure (i.e., density). The effects of buoyant convection may then be determined by a comparison of the mu-g and 1g results. The mu-g results may also reveal the important mechanisms which are masked or overwhelmed by buoyant convection in 1g. A numerical investigation was conducted using a validated, time-accurate numerical code to study the underlying physics during the flame interaction and to assist the interpretation of the experimental results.

Sheu, Jyh-Cherng↗

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↗

A preliminary study of the tropical water cycle and its sensitivity to surface warming

The Goddard Cumulus Ensemble Model (GCEM) has been used to demonstrate that cumulus-scale dynamics and microphysics play a major role in determining the vertical distribution of water vapor and clouds in the tropical atmosphere. The GCEM is described and is the basic structure of cumulus convection. The long-term equilibrium response to tropical convection to surface warming is examined. A picture of the water cycle within tropical cumulus clusters is developed.

Lau, K. M.↗

On the gravitational stability of the interstellar medium in spiral arms

The dynamics of self-gravitating gaseous density perturbations in the tightly wound spiral-arm model of the forcing potential (Brown 1969) is investigated. The cloud ensemble of the cold phase of the interstellar medium is regarded as a collisional gas whose 'molecules' are the clouds themselves, and the dynamical evolution of density perturbations is calculated assuming this gas obeys a polytropic pressure-density law. It is suggested that the formation of extreme Population I objects found predominantly along the arms of Sb/Sc systems is regulated by an overall energy balance in the cloud fluid. A single parameter is found which appears to characterize the flow's gravitational response to perturbations; this spiral parameter, Q-sub-sp, differs from the Toomre Q parameter by a factor of the square root of the density enhancement at the shock. It is necessary but not sufficient that Q-sub-sp be less than unity for significant growth to occur.

Balbus, S. A.↗

Deconvolution of dynamic heterogeneity in protein structure

Heterogeneity is intrinsic to the dynamic process of a chemical reaction. As reactants are converted to products via intermediates, the nature and extent of heterogeneity vary temporally throughout the duration of the reaction and spatially across the molecular ensemble. The goal of many biophysical techniques, including crystallography and spectroscopy, is to establish a reaction trajectory that follows an experimentally provoked dynamic process. It is essential to properly analyze and resolve heterogeneity inevitably embedded in experimental datasets. We have developed a deconvolution technique based on singular value decomposition (SVD), which we have rigorously practiced in diverse research projects. In this review, we recapitulate the motivation and challenges in addressing the heterogeneity problem and lay out the mathematical foundation of our methodology that enables isolation of chemically sensible structural signals. We also present a few case studies to demonstrate the concept and outcome of the SVD-based deconvolution. Finally, we highlight a few recent studies with mechanistic insights made possible by heterogeneity deconvolution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computing Nonequilibrium Responses with Score-Shifted Stochastic Differential Equations

Using equilibrium fluctuations to understand the response of a physical system to an externally imposed perturbation is the basis for linear response theory, which is widely used to interpret experiments and shed light on microscopic dynamics. For nonequilibrium systems, perturbations cannot be interpreted simply by monitoring fluctuations in a conjugate observable and general response results rely on path ensemble averaging. Furthermore, these techniques do not apply to perturbations that affect the diffusion tensor in a stochastic system. Here, we introduce an “effective” physical process that represents the diffusion perturbed dynamics and enables accurate calculations of responses to a change in the diffusion. Interestingly, the effective dynamics contain an additional drift involving the instantaneous “score” of the system, and we leverage score matching algorithms to carry out nonequilibrium response calculations on systems for which the exact stationary distribution is unknown.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-Performance Semiempirical Excited-State Molecular Dynamics Powered by Graphics Processing Units

Here, this Letter introduces excited-state molecular dynamics in PYSEQM, a GPU-accelerated semiempirical quantum chemistry engine implemented in PyTorch. The new module enables Born–Oppenheimer molecular dynamics (BOMD) using configuration-interaction singles and random phase approximation for excited states, allowing long trajectories and large statistical ensembles to be simulated efficiently on a single GPU. We also implement an extended Lagrangian excited-state BOMD (XL-ESMD) scheme that propagates auxiliary electronic variables, enabling relaxed ground and excited-state convergence thresholds without compromising energy conservation. The excited-state BOMD implementation scales smoothly from small chromophores to a nearly 900-atom dendrimer (taking 6.5 s per MD step). PYSEQM also supports batched execution, allowing many geometries or trajectories to be evaluated in a single GPU launch, substantially increasing throughput and making ensemble-based protocols routine. As a demonstration, we compute absorption, emission, and infrared spectra from trajectories propagated on the ground and first excited states. The XL-ESMD scheme yields identical spectra at significantly lower computational cost, establishing the role of extended Lagrangian based dynamics for efficient excited-state BOMD simulations. Beyond raw performance, PYSEQM’s PyTorch foundation provides automatic differentiation for forces, efficient GPU batching, and seamless interfacing with machine learning models. These capabilities position PYSEQM as a practical platform for machine learning-augmented excited-state dynamics and lay the foundation for future data-driven nonadiabatic excited-state dynamics modeling of ultrafast spectroscopic probes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A collision operator for describing dissipation in noncanonical phase space

The phase space of a noncanonical Hamiltonian system is partially inaccessible due to dynamical constraints (Casimir invariants) arising from the kernel of the Poisson tensor. When an ensemble of noncanonical Hamiltonian systems is allowed to interact, dissipative processes eventually break the phase space constraints, resulting in a thermodynamic equilibrium described by a Maxwell–Boltzmann distribution. However, the time scale required to reach Maxwell–Boltzmann statistics is often much longer than the time scale over which a given system achieves a state of thermal equilibrium. Examples include diffusion in rigid mechanical systems, as well as collisionless relaxation in magnetized plasmas and stellar systems, where the interval between binary Coulomb or gravitational collisions can be longer than the time scale over which stable structures are self-organized. Here, we focus on self-organizing phenomena over spacetime scales such that particle interactions respect the noncanonical Hamiltonian structure, but yet act to create a state of thermodynamic equilibrium. We derive a collision operator for general noncanonical Hamiltonian systems, applicable to fast, localized interactions. This collision operator depends on the interaction exchanged by colliding particles and on the Poisson tensor encoding the noncanonical phase space structure, is consistent with entropy growth and conservation of particle number and energy, preserves the interior Casimir invariants, reduces to the Landau collision operator in the limit of grazing binary Coulomb collisions in canonical phase space, and exhibits a metriplectic structure. We further show how thermodynamic equilibria depart from Maxwell–Boltzmann statistics due to the noncanonical phase space structure, and how self-organization and collisionless relaxation in magnetized plasmas and stellar systems can be described through the derived collision operator.

Boltzmann equation↗

Projected evolution of droughts and human exposure in the contiguous United States under SSP5-8.5: a regional downscaling perspective

The increasingly unpredictable and extreme weather patterns under a warming climate underscore the urgency of accurate regional assessments of future drought risk. This study evaluates the projected drought evolution in the contiguous United States under the high-emission shared socioeconomic pathway 5–8.5 climate scenario for the coming decades. Using a multi-model ensemble of six Coupled Model Intercomparison Project Phase 6 global climate models combined with dynamical downscaling techniques, we analyzed near-term (2020–2039) and mid-term (2040–2059) drought patterns using the self-calibrating palmer drought severity index (ScPDSI), the standardized precipitation index (SPI-12), and the Standardized Precipitation-Evapotranspiration Index (SPEI-12). Results reveal a widespread increase in abnormally dry (D0) and moderate drought (D1) conditions, particularly in urban areas, while severe (D2), extreme (D3), and exceptional (D4) droughts are expected to become less common in many regions. Meanwhile, persistent and intensifying droughts are projected in the western and southwestern U.S., driven by long-term soil moisture deficits. The ScPDSI projects that 1.1 million urban residents will be affected by D0 conditions in 2050, while SPI-12 suggests a decrease in the total affected populations after 2040. ScPDSI indicates prolonged droughts in the West, and SPI-12 captures transient variability. Although the total drought-exposed population is expected to decrease, urban areas will continue to bear a greater burden, particularly for mild droughts (D0, D1). These findings highlight a shift toward more frequent mild droughts, fewer severe droughts, and persistent drying in the Southwest, emphasizing the need for region-specific adaptation strategies.

CMIP6↗

Implementation of variable time step stochastic dynamics for electronically inelastic gas-surface collisions

A variable time step algorithm has been implemented for solving the stochastic equations of motion for gas-surface collisions. It has been tested for a simple model of electronically inelastic collisions with an insulator surface in which the phonon manifold acts as a heat bath and electronic states are localized. In addition to reproducing the accurate nuclear dynamics of the surface atoms, numerical calculations have shown the algorithm to yield accurate ensemble averages of physical observables such as electronic transition probabilities and total energy loss of the gas atom to the surface. This new algorithm offers a gain in efficieny of up to an order of magnitude compared to fixed time step integration.

Garrett, Bruce C.↗

ATLAS: Airborne Tunable Laser Absorption Spectrometer for stratospheric trace gas measurements

The ATLAS instrument is an advanced technology diode laser based absorption spectrometer designed specifically for stratospheric tracer studies. This technique was used in the acquisition of N2O tracer data sets on the Airborne Antarctic Ozone Experiment and the Airborne Arctic Stratospheric Expedition. These data sets have proved valuable for comparison with atmospheric models, as well as in assisting in the interpretation of the entire ensemble of chemical and meteorological data acquired on these two field studies. The N2O dynamical tracer data set analysis revealed several ramifications concerning the polar atmosphere: the N2O/NO(y) correlation, which is used as a tool to study denitrification in the polar vertex; the N2O Southern Hemisphere morphology, showing subsidence in the winter polar vortex; and the value of the N2O measurements in the interpretation of ClO, O3, and NO(y) measurements and of the derived dynamical tracer, potential vorticity. Field studies also led to improved characterization of the instrument and to improved accuracy.

Loewenstein, Max↗

Spin squeezing of macroscopic nuclear spin ensembles

Spin squeezing has been explored in atomic systems as a tool for quantum sensing, improving experimental sensitivity beyond the spin standard quantum limit for certain measurements. To optimize absolute metrological sensitivity, it is beneficial to consider macroscopic spin ensembles, such as nuclear spins in solids and liquids. Coupling a macroscopic spin ensemble to a parametrically-modulated resonant circuit can create collective spin squeezing by generating spin correlations mediated by the circuit. We analyze the squeezing dynamics in the presence of decoherence and finite spin polarization, showing that achieving 7 dB spin squeezing is feasible in several nuclear spin systems. The metrological benefit of squeezing a macroscopic spin ensemble lies in the suppression of technical noise sources in the spin detection system relative to the spin projection noise. This expands the experimental sensitivity bandwidth when searching for signals of unknown frequency and can improve the resonant signal-to-noise ratio. Squeezing macroscopic spin ensembles may prove to be a useful technique for fundamental physics experiments aimed at detecting spin interactions with oscillating background fields, such as ultralight dark matter. Published by the American Physical Society 2025

Boyers, Eric↗

Integrated machine learning-molecular dynamics framework for electrolyte property prediction

Electrochemical stability windows determine the operating range of battery electrolytes, yet accurate prediction remains challenging because stability emerges from statistical ensembles of local solvation environments rather than single ground-state molecular structures. Traditional density functional theory calculations on energy-minimized clusters cannot capture the thermal variations in local coordination environments and geometries that govern decomposition, while SMILES-based machine learning methods lack explicit representation of three-dimensional solvation structure and ion pairing. Here, we introduce a structure-aware machine learning framework that predicts frontier orbital energies (HOMO and LUMO) directly from molecular dynamics-sampled solvation configurations, achieving sub-0.6 eV accuracy at computational costs 3–4 orders of magnitude lower than first-principles methods. Across twelve representative battery electrolytes, we demonstrate that solvent-separated and contact ion pairs exhibit strong size- and local chemistry dependent electronic stability, with variations in coordination shifts of HOMO or LUMO level by 2–3 eV, and that extended solvation structure and partially desolvated environment further modulate stability by up to 3 eV. By encoding the statistical nature of electrochemical failure through ensemble sampling of explicit solvation geometries, our approach enables high-throughput screening and rational design of next-generation battery electrolytes with mechanistic understanding of structure–property relationships.

Energy - Storage↗

HydraGNN_Predictive_GFM_2024 - Ensemble of predictive graph foundation models for ground state atomistic materials modeling

We provide the ensemble of fifteen pre-trained graph foundation models (GFMs) for atomistic materials modeling applications. Each one of the fifteen GFMs has been trained on five open-source datasets that (once aggregated) amount to over 154 million atomistic structures, which cover over two-thirds of the natural elements of the periodic table and that comprises a broad set of organic and inorganic compounds. This vast set of atomistic structures comprises ground state configurations that are dynamically stable (i.e., equilibrated structures with atomic forces approximately close to zero values) as well as dynamically unstable structures (i.e., non-equilibrium structures with non-negligible non-zero values of atomic forces). The ensemble of datasets aggregated does NOT include excited states. The datasets have been curated to remove atomistic structures with spectral norm of the force tensor above 100 eV/angstrom. Moreover, a linear term of the energy was computed for each dataset using a linear regression model that uses the chemical concentration of each natural element as regressor. The linear term predicted by the linear regression model has been subtracted from each original energy value to perform a re-alignment of the energy values across different electronic structures approximation theories performed to generate the diverse multi-source, multi-fidelity datasets. The folder "ADIOS_files" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "ADIOS_files" directory contains 6 sub-directories named as follows: - ANI1x-v3.bp - MPTrj-v3.bp - OC2020-20M-v3.bp - OC2020-v3.bp - OC2022-v3.bp - qm7x-v3.bp Each sub-directory contains the pre-processed datasets converted in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used to the development, training, and performance testing of the ensemble go predictive graph foundation models. Each GFM was developed using HydraGNN (https://github.com/ORNL/HydraGNN) as underlying graph neural network (GNN) architecture. The multi-task learning (MTL) capability of HydraGNN was used to simultaneously train the GFMs on labeled values for direct predictions of energy (a total system property of an atomistic structure that measures the chemical stability) and atomic forces (an atomic level property of an atomistic structure that measures the dynamical stability). The hyper parameters of the GFM have been tuned using scalable hyperparameter optimization (HPO) algorithms implemented in the software DeepHyper (https://github.com/deephyper/deephyper). The pre-training of each HPO trial was performed using distributed data parallelism (DDP) to scale the training across 128 compute nodes of the exascale OLCF supercomputer Frontier. Each HPO trial was trained only for 10 epochs and an early stopping was performed to avoid wasting significant computational resources on GNN architectures that were clearly underperforming. For each HPO trial, the 'omnistat' tool developed by (AMD Research - Advanced Micro Device) was used to measure the total energy consumption in kWh. The ensemble of GFMs was obtained by selecting the fifteen best performing HPO trials. Four models have been selected for their clear advantage in accuracy, and these are the GFMs with IDs 229, 156, 147, 260. Additional eleven models have been selected based on judicious balance between accuracy and energy consumption needed for training, and these are the GFMs with IDs 165, 78, 137, 1, 175, 171, 181, 67, 179, 167, 351. Each selected GFM of the ensemble was continued to cumulate a total of at most 30 epochs. In some cases, the total number of epochs actually performed was les than 30 due to two combined factors: (1) the size of the GFM (i.e., the number of model parameters to train) and (2) the total wall-clock time for which the computational resources could be allocated on OLCF-Frontier. The "Ensemble_of_models" directory contains 15 sub-directories named as follows: - gfm_0.229 - gfm_0.156 - gfm_0.147 - gfm_0.260 - gfm_0.165 - gfm_0.78 - gfm_0.137 - gfm_0.1 - gfm_0.175 - gfm_0.171 - gfm_0.181 - gfm_0.67 - gfm_0.179 - gfm_0.167 - gfm_0.351 Each one of these sub-directories refers to one of the fifteen HPO trials that have been selected to continue the pre-training with at most 30 epochs. With each sub-directory associated with a specific HPO trial, the following files can be found: - config.json: file for argument parsing to develop and train an HydraGNN architecture - gfm_0.ID_epoch_N.pk: file with model parameters for HPO ID trial after N epochs of training The ensemble of fifteen GFM architectures was used for (1) ensemble averaging to stabilize the predictions of energy and atomic forces after pre-training for post-processing analysis and (2) ensemble uncertainty quantification (UQ). The code used to develop, pre-train, and load the pre-trained models for post-processing analysis is available on the ORNL-GitHub at the following link: https://github.com/ORNL/HydraGNN/tree/Predictive_GFM_2024

36 MATERIALS SCIENCE↗

Laplace-transform technique for deriving thermodynamic equations from the classical microcanonical ensemble

A direct and convenient method is presented for deriving expressions which equate any thermodynamic state function to averages of specific dynamical functions and their fluctuations over the classical microcanonical distribution. Specific expressions are obtained for a variety of thermodynamic quantities. The effect of various entropy definitions on the results are assessed, and the latter are compared to previous work in the literature. The derived formulas are applied to the analysis of molecular-dynamics computer simulations.

Pearson, E. M.↗

Application of an Ensemble Smoother to Precipitation Assimilation

Assimilation of precipitation in a global modeling system poses a special challenge in that the observation operators for precipitation processes are highly nonlinear. In the variational approach, substantial development work and model simplifications are required to include precipitation-related physical processes in the tangent linear model and its adjoint. An ensemble based data assimilation algorithm "Maximum Likelihood Ensemble Smoother (MLES)" has been developed to explore the ensemble representation of the precipitation observation operator with nonlinear convection and large-scale moist physics. An ensemble assimilation system based on the NASA GEOS-5 GCM has been constructed to assimilate satellite precipitation data within the MLES framework. The configuration of the smoother takes the time dimension into account for the relationship between state variables and observable rainfall. The full nonlinear forward model ensembles are used to represent components involving the observation operator and its transpose. Several assimilation experiments using satellite precipitation observations have been carried out to investigate the effectiveness of the ensemble representation of the nonlinear observation operator and the data impact of assimilating rain retrievals from the TMI and SSM/I sensors. Preliminary results show that this ensemble assimilation approach is capable of extracting information from nonlinear observations to improve the analysis and forecast if ensemble size is adequate, and a suitable localization scheme is applied. In addition to a dynamically consistent precipitation analysis, the assimilation system produces a statistical estimate of the analysis uncertainty.

Zhang, Sara↗

Summertime Transport Pathways from Different Northern Hemisphere Regions into the Arctic

Trace gases and aerosols play an important role in Arctic chemistry and climate. As most Arctic tracers and aerosols are transported from midlatitude source regions, long-range transport into the Arctic is one of the key factors to understand the current and future states of Arctic climate. While previous studies have investigated the airmass fraction and transit time distribution in the Arctic, the actual transport pathways and their underlying dynamics and efficiencies are yet to be understood. In this study, we implement a large ensemble of idealized tagged pulse passive tracers in the Whole Atmosphere Community Climate Model version 5 to identify and analyze summertime transport pathways from different Northern Hemisphere surface regions into the Arctic. Three different transport pathways are identified as those associated with fast, intermediate and slow time scales. Midlatitude tracers can be transported into the Arctic in the troposphere via the fast transport pathway (~8 days), which moves tracers northward from the source region mainly through transient eddies. For the intermediate transport pathway, which happens on 1~3 weeks’ time scales, midlatitude tracers are first zonally transported by the jet stream, and then advected northward into the Arctic over Alaska and northern North Atlantic. Tropical and subtropical tracers are transported into the Arctic lower stratosphere via the slow transport pathway (1~3 months), as the tracers are lifted upward into the tropical and subtropical lower stratosphere, and then transported into the Arctic following the isentropic surfaces.

Trace gases↗

Quantifying the role of ocean coupling in Arctic amplification and sea-ice loss over the 21st century

The enhanced warming of the Arctic, relative to other parts of the Earth, a phenomenon known as Arctic amplification, is one of the most striking features of climate change, and has important climatic impacts for the entire Northern Hemisphere. Several mechanisms are believed to be responsible for Arctic amplification; however, a quantitative understanding of their relative importance is still missing. Here, using ensembles of model integrations, we quantify the contribution of ocean coupling, both its thermodynamic and dynamic components, to Arctic amplification over the 20th and 21st centuries. We show that ocean coupling accounts for ~80% of the amplification by 2100. In particular, we show that thermodynamic coupling is responsible for future amplification and sea-ice loss as it overcomes the effect of dynamic coupling which reduces the amplification and sea-ice loss by ~35%. Our results demonstrate the utility of targeted numerical experiments to quantify the role of specific mechanisms in Arctic amplification, for better constraining climate projections.

Ocean heat transport↗