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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Consequences of the failure of equipartition for the p – V behavior of liquid water and the hydration free energy components of a small protein

Earlier we showed that in the molecular dynamics simulation of a rigid model of water it is necessary to use an integration time-step δ t ≤ 0.5 fs to ensure equipartition between translational and rotational modes. Here we extend that study in the NVT ensemble to NpT conditions and to an aqueous protein. We study neat liquid water with the rigid, SPC/E model and the protein BBA (PDB ID: 1FME) solvated in the rigid, TIP3P model. We examine integration time-steps ranging from 0.5 fs to 4.0 fs for various thermostat plus barostat combinations. We find that a small δ t is necessary to ensure consistent prediction of the simulation volume. Hydrogen mass repartitioning alleviates the problem somewhat, but is ineffective for the typical time-step used with this approach. The compressibility, a measure of volume fluctuations, and the dielectric constant, a measure of dipole moment fluctuations, are also seen to be sensitive to δ t . Using the mean volume estimated from the NpT simulation, we examine the electrostatic and van der Waals contribution to the hydration free energy of the protein in the NVT ensemble. These contributions are also sensitive to δ t . In going from δ t = 2 fs to δ t = 0.5 fs, the change in the net electrostatic plus van der Waals contribution to the hydration of BBA is already in excess of the folding free energy reported for this protein.

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

A Polynomial-Time Classical Algorithm for Noisy Quantum Circuits

We provide a polynomial-time classical algorithm for noisy quantum circuits. The algorithm computes the expectation value of any observable for any circuit, with a small average error over input states drawn from an ensemble (e.g., the computational basis). Our approach is based upon the intuition that noise exponentially damps nonlocal correlations relative to local correlations. This enables one to classically simulate a noisy quantum circuit by keeping track of only the dynamics of local quantum information. Our algorithm also enables sampling from the output distribution of a circuit in quasipolynomial time, so long as the distribution anticoncentrates. A number of implications are discussed, including a fundamental limit on the efficacy of noise mitigation strategies: For constant noise rates, any quantum circuit for which error mitigation succeeds in polynomial-time on most input states can also be classically simulated in polynomial-time on most input states. Our algorithms scale exponentially in the inverse noise rate, which is fundamental and makes them impractical for current quantum devices.

decoherence↗
DOE OSTI2026

The Sparks tracking code

This note introduces a new tracking code, Sparks, built on the foundations laid by hemod4 and designed for fast tracking of large ensembles of single particles over thousands of turns. The code is intended for the rapid evaluation of nonlinear lattice behavior through observables such as detuning, Poincaré sections, and smear, providing a computationally efficient alternative to large-scale dynamic aperture studies during lattice design and optimization.

43 PARTICLE ACCELERATORS↗
NASA NTRS1973

Statistical mechanics of negative temperature states

The dynamics of two-dimensional interacting line vortices is identical to that of the two-dimensional electrostatic guiding center plasma. Both are Hamiltonian systems and are therefore susceptible to statistical mechanical treatments. The predictions of the microcanonical ensemble are explored for this system. Interest focuses primarily on the regime of total positive interaction energy, which should be above the Onsager negative temperature threshold. Calculations of the probability distribution for a component by means of the central limit theorem are carried out in the manner of Khinchin. The probability distribution of a component reduced to the usual Gibbs distribution in the regime of positive temperatures, and is still explicitly calculable for negative temperatures. The negative temperature states are neither quiescent nor spatially uniform. Expressions for the temperature are explicitly provided in terms of the total particle energy and particle number. A BBGKY hierarchy can be derived for both temperature regimes. Numerical simulations involving solutions of the equations of motion of 4008 particles are presented.

Montgomery, D. C.↗
DOE OSTI2024

Elucidating and contrasting the mechanisms for Mg and Ca sulfate ion-pair formation with multi-level embedded quantum mechanics/molecular dynamics simulations

Here, solutions and minerals containing sulfate (SO 4 2- ), and Ca 2+ and Mg 2+ cations, are ubiquitous throughout the lithosphere and are significant components of seawater, thus presenting a prototypical system for the study of strong electrolytes and crystal nucleation mechanisms. However, despite their relative abundance, key questions remain unanswered about the most fundamental atomic-level steps of their mineralization pathways and aqueous dynamics. Here, we carry out enhanced sampling multi-level molecular dynamics (MD) embedded correlated wavefunction theory simulations to elucidate ion-pairing mechanisms for Mg–SO 4 and Ca–SO 4 in concentrated aqueous solution, accurately capturing effects arising from both structural dynamics and electron exchange–correlation. We predict contact-ion-pair formation to be barrierless and highly exoergic for Ca–SO 4 , in agreement with its minimal solubility, whereas for Mg–SO 4 , solvent-shared and contact ion pairs have similar free energies, qualitatively consistent with its higher solubility. Finally, we demonstrate that brief high-temperature pre-equilibration may be utilized to accelerate convergence of free energies in blue-moon-ensemble enhanced-sampling MD.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2017

An Intercomparison of Changes Associated with Earth's Lower Tropospheric Temperature Using Traditional and AMIP-Style Reanalyses

Reanalyses have become an integral tool for evaluating regional and global climate variations, and an important component of this is modifications to the energy budget. Reductions in Arctic Sea ice extent has induced an albedo feedback, causing the Arctic to warm more rapidly than anywhere else in the world, referred to as "Arctic Amplification." This has been demonstrated by observations and numerous reanalyses, including the Modern Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2). However, the Arctic Amplification signal is non-existent in a ten member ensemble of the MERRA-2 Atmospheric Model Intercomparison Project (M2AMIP) simulations, using the same prescribed climate forcing, including Sea Surface Temperature (SST) and ice. An evaluation of the temperature tendency within the lower troposphere due to radiation, moisture, and dynamics as well as the surface energy budget in MERRA-2 and M2AMIP will demonstrate that despite identical prescribed SSTs and sea ice in both versions, enhanced warming in the Arctic in MERRA-2 is in response to the analysis increment tendency due to temperature observations. Furthermore, the role of boundary conditions, model biases and changes in observing systems on the Arctic Amplification signal will be assessed. Literature on the topic of Arctic Amplification demonstrates that the enhanced warming begins in the mid-1990s. Anomalously warm Arctic SST in the early time period of MERRA-2 can mute the trend in Arctic lower troposphere temperature without the constraint of observations in M2AMIP. Additionally, MERRA-2 uses three distinct datasets of SST and sea ice concentration, which also plays a role.

Arcti↗
NASA NTRS2020

Dynamical Forecasts of Tropical Terrestrial Carbon Fluxes with the NASA S2S Retrospective Forecast System

Recent advances in the ability to predict climate anomalies at sub-seasonal to seasonal (S2S) timescales allow us to explore the possibility of forecasting carbon flux anomalies. Although carbon flux forecasting is a relatively new concept, it is potentially beneficial as it can help us better understand global and regional land-atmosphere carbon feedbacks associated with climate variations and can provide guidance for future field mission design. Here we evaluate the skill of forecasted terrestrial carbon anomalies generated from meteorological anomalies produced with the NASA Global Modeling and Assimilation Office (GMAO) S2S forecast system. We focus here on three representative time periods (the most recent 2015-2016 El Nino, 2011 La Nina, and 2014 as a neutral year), with each corresponding 9-month forecast comprising four ensemble members initialized in the preceding December. The meteorological variables produced by the GMAO forecast system were bias-corrected using a climatology derived from the Modern Era Retrospective analysis for Research and Applications, version 2 (MERRA-2) before being used to drive a suite of offline simulations with the NASA Catchment-CN terrestrial biosphere model, a model that computes water-energy-carbon dynamics. Forecasts are evaluated by comparing against satellite-driven estimates of gross primary production (GPP) and inverse model estimates of net carbon flux that incorporate satellite carbon dioxide measurements. We find that the restrospectively predicted carbon fluxes in the tropics reasonably reproduce the signs and magnitudes of the observed anomalies between the 2015-2016 El Nino and the 2011 La Nina for both net flux and GPP. For instance, for the El Nino period, the magnitude of the forecasted negative GPP anomaly in the South American tropics (which undergoes anomalously warm and dry conditions) agrees with the observed GPP anomaly at leads of up to three or four months. Overall, this study demonstrates potential skill in the forecast of biospheric carbon fluxes a few months in advance, a capability that could contribute to attribution studies focusing on carbon flux variations and support innovative observation strategies in the future.

Lee, Eunjee↗
DOE OSTI2025

Stochastic Ensemble Generation for Improved Characterization of Representing Geologic Variability in a Reservoir: IBDP Case Study for SMART Initiative

This document is a poster covering the findings from activities on training data generation, specifically geologic ensemble generation. The generated geologic realizations captured the range of possible permeability distributions of the subsurface at the Illinois Basin - Decatur Project (IBDP) site, based on available well log variabilities. The percentages of reservoirs and baffles in the injection zone and a truncation of baffle permeability led to more variance in the simulations. This will be used to build forward modeling, history matching, and optimization workflows. The geologic realizations were also ranked according to dynamic measures of hydraulic diffusivity, and simulations confirm a greater contrast between the reservoir and the baffles during injection.

stochastic ensemble generation↗
DOE OSTI2024

Unravelling the orbits of cluster galaxy populations according to their dominant gas ionization source

ABSTRACT We investigate the kinematical and dynamical properties of cluster galaxy populations classified according to their dominant source of gas ionization, namely: star-forming (SF) galaxies, optical active galactic nuclei (AGNs), mixed SF plus AGN ionization (transition objects, T), and quiescent (Q) galaxies. We stack 8892 member galaxies from 336 relaxed galaxy clusters to build an ensemble cluster and estimate the observed projected profiles of numerical density and velocity dispersion, $\sigma _P(R)$, of each galaxy population. The MAMPOSSt code and the Jeans equations inversion technique are used to constrain the velocity anisotropy profiles of the galaxy populations in both parametric and non-parametric ways. We find that Q (SF) galaxies display the lowest (highest) typical cluster-centric distances and velocity dispersion values. Transition galaxies are more concentrated and tend to exhibit lower velocity dispersion values than SF galaxies. Galaxies that host an optical AGN are as concentrated as Q galaxies but display velocity dispersion values similar to those of the SF population. MAMPOSSt is able to find equilibrium solutions that successfully recover the observed $\sigma _P(R)$ profile only for the Q, T, and AGN populations. We find that the orbits of all populations are consistent with isotropy in the inner regions, becoming increasingly radial with the distance from the cluster centre. These results suggest that Q galaxies are in equilibrium within their clusters, while SF galaxies have more recently arrived in the cluster environment. Finally, the T and AGN populations appear to be in an intermediate dynamical state between those of the SF and Q populations.

Valk, Greique A. (ORCID:0009000827731299)↗
DOE OSTI2024

Earth's record-high greenness and its attributions in 2020

Terrestrial vegetation is a crucial component of Earth's biosphere, regulating global carbon and water cycles and contributing to human welfare. Despite an overall greening trend, terrestrial vegetation exhibits a significant inter-annual variability. The mechanisms driving this variability, particularly those related to climatic and anthropogenic factors, remain poorly understood, which hampers our ability to project the long-term sustainability of ecosystem services. Here, in this work, by leveraging diverse remote sensing measurements, we pinpointed 2020 as a historic landmark, registering as the greenest year in modern satellite records from 2001 to 2020. Using ensemble machine learning and Earth system models, we found this exceptional greening primarily stemmed from consistent growth in boreal and temperate vegetation, attributed to rising CO 2 levels, climate warming, and reforestation efforts, alongside a transient tropical green-up linked to the enhanced rainfall. Contrary to expectations, the COVID-19 pandemic lockdowns had a limited impact on this global greening anomaly. Our findings highlight the resilience and dynamic nature of global vegetation in response to diverse climatic and anthropogenic influences, offering valuable insights for optimizing ecosystem management and informing climate mitigation strategies.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1998

The Dynamics of Coupled Oscillator Phase Control

Arrays of coupled oscillators have been proposed as means of realizing high power rf sources via coherent spatial power combining. In such applications, a uniform phase distribution over the aperture is usually desired. However, it has been shown that by detuning some of the oscillators away from the oscillation frequency of the ensemble of oscillators, one may achieve other useful aperture phase distributions. Of particular interest among those achievable are linear phase distributions because these result in steering of the output rf beam away from the broadside direction. The theory describing the behavior of such arrays of coupled oscillators is quite complicated since the phenomena involved are inherently nonlinear. However, a simplified theory has been developed which facilitates intuitive understanding. This simplified theory is based on a "continuum model" in which the aperture phase is represented by a continuous function of the aperture coordinates. A challenging aspect of the development of this theory is the derivation of appropriate boundary conditions at the edges or ends of the array.

Pogorzelski, R. J.↗
NASA NTRS1981

Inertial measurement unit pre-processors and post-flight STS-1 comparisons

The flight results show that the relative tri-redundant Inertial Measurement Unit IMU performance throughout the entire entry flight was within the expected accuracy. Comparisons are presented which show differences in the accumulated sensed velocity changes as measured by the tri-redundant IMUs (in Mean Equator and Equinox of 1950.0), differences in the equivalent inertial Euler angles as measured with respect to the M50 system, and finally, preliminary instrument calibrations determined relative to the ensemble average measurement set. Also, differences in the derived body axes rates and accelerations are presented. Because of the excellent performance of the IMUs during the STS-1 entry, the selection as to which particular IMU would best serve as the dynamic data source for entry reconstruction is arbitrary.

Findlay, J. T.↗
NASA NTRS2018

Multidecadal Changes in Lower Stratospheric Ozone: Variability Vs. Trends

As upper stratospheric ozone appears to be recovering as a result of decreasing chlorine loading following the implementation of the Montreal Protocol and its amendments and in agreement with model projections, several recent studies report an apparent decline of ozone concentrations in the lower stratosphere in the last two decades, particularly in the extratropics. Our previous work as well as at least two other studies provide evidence that this decline results from transport changes rather than an intensification of chemical depletion. It remains unclear whether these changes represent long-term internal variability or are a consequence of a climate forcing. Here we perform free-running ensembles of the recent past (1980-2016) using the Goddard Earth Observing System Model (GEOS) at the cubed sphere C180 (approximately half degree) resolution. Two suites of 10-member ensembles are performed, one in which observed sea surface temperature (SSTs) are fully prescribed, and the other in which the linear SST trend over the recent past is removed so as to only retain internal variability. We evaluate the trends in both ozone as well as two idealized tracers with prescribed uniform loss that are used to isolate the role of transport from chemistry and emissions. Probability-distribution-functions of the trends in both ozone and idealized tracers are compared among ensemble members and with observed trends in order to evaluate the likelihood of recent observed declines in lower stratospheric ozone, relative to large internal variability. Moreover, comparisons among simulations with and without imposed SST trends indicate the extent to which dynamically-driven ozone trends reflect forced trends or internal variability in lower stratospheric dynamics.

Wargan, Krzysztof↗
NASA NTRS2021

Assimilation of DAWN Doppler Wind Lidar Data During the 2017 Convective Processes Experiment (CPEX): Impact on Precipitation and Flow Structure

An improved representation of 3-D air motion and precipitation structure through forecast models and assimilation of observations is vital for improvements in weather forecasting capabilities. However, there are few independent data to properly validate a model forecast of precipitation structure when the underlying dynamics are evolving on short convective timescales. Using data from the JPL Ku/Ka-band Airborne Precipitation Radar (APR-2) and the 2 μmDoppler Aerosol Wind (DAWN) lidar collected during the2017 Convective Processes Experiment (CPEX), the NASA Unified Weather Research and Forecasting (WRF) Ensemble Data Assimilation System (EDAS) modeling system was used to quantify the impact of high-resolution sparsely sampled DAWN measurements on the analyzed variables and on the forecast when the DAWN winds were assimilated. Over-all, the assimilation of the DAWN wind profiles had a discernible impact on the wind field as well as the evolution and timing of the 3-D precipitation structure. Analysis of individual variables revealed that the assimilation of the DAWN winds resulted in important and coherent modifications of the environment. It led to an increase in the near-surface convergence, temperature, and water vapor, creating more favorable conditions for the development of convection exactly where it was observed (but not present in the control run). Comparison to APR-2 and observations by the Global Precipitation Measurement (GPM) satellite shows a much-improved forecast after the assimilation of the DAWN winds – development of precipitation where there was none, more organized precipitation where there was some, and a much more intense and organized cold pool, similar to the analysis of the dropsonde data. The onset of the vertical evolution of the precipitation showed similar radar-derived cloud-top heights, but delayed in time. While this investigation was limited to a single CPEX flight date, the investigation design is appropriate for further investigation of the impact of airborne Doppler wind lidar observations upon short-term convective precipitation forecasts

DAWN↗
NASA NTRS2018

Impact of Satellite Sea Surface Salinity Observations on ENSO Predictions from the GMAO S2S Forecast System

Sea surface salinity (SSS) observations from space allow us to investigate if improved estimates of near-surface density stratification and associated mixing can positively impact seasonal to interannual variability of tropical Pacific Ocean dynamics as well as dynamical ENSO forecasts. For the first part of the presentation, we utilize our intermediate-complexity coupled model. Baseline experiments assimilate satellite sea level (multi-satellite gridded AVISO, 2013), SST (Reynolds et al., 2004), and in situ subsurface temperature and salinity observations (GTSPP NODC, 2006). These baseline experiments are then compared with experiments that additionally assimilate Aquarius (V5.0 Lilly and Lagerloef, 2008) and SMAP (V4.0 Fore et al., 2016) SSS. Twelve-month forecasts are initialized for each month from September 2011 to September 2017. For initialization of the coupled forecast, the positive impact of SSS assimilation is brought about by surface freshening near the eastern edge of the western Pacific warm pool and density changes that lead to shallower mixed layer between 10S-5N. This pattern enhances air/sea interaction and amplifies the equatorial Kelvin wave signal. We find that including satellite SSS significantly improves NINO3.4 sea surface temperature anomaly validation over most forecast lead times. We next assess how different satellite SSS products impact the validation of ENSO forecasts. SMAP V4 reduces the salty bias in the western Pacific and so is an improvement upon SMAP V2 and SMOS V2 (Boutin et al., 2017) has similar validation characteristics as a combination of Aquarius and SMAP V4. Next we shift to present results from the NASA GMAO Sub-seasonal to seasonal (S2S_v2.1) production coupled model (i.e. the same model that contributes ENSO forecasts to the North American Multi-Model Ensemble Experiment). From March to June 2015, the availability of two overlapping satellite SSS instruments, Aquarius and SMAP, allows a unique opportunity to compare and contrast forecasts initialized with the benefit of these two satellite SSS observation types. We assess the impact of satellite sea surface salinity (SSS) observations on dynamical ENSO forecasts for the big 2015 El Nino event. We will present distinct experiments for the overlap period that include 1) freely evolving SSS (i.e. no satellite SSS as the production system), 2) Aquarius, and 3) SMAP initialization. Our results show that using Aquarius slightly improves validation of the reanalysis (including sea level and temperature statistics). Our production system without SSS assimilation generated too warm forecasts for the 2015 El Nino from March initial conditions. Incorporating Aquarius into initialization of the coupled system leads to a deeper, more realistic MLD that acts to damp the downwelling Kelvin signal and slightly cool NINO3.4 SST. With Aquarius the forecasts better match the observed amplitude of the 2015 event. On the other hand, SMAP V2 relaxation generally degrades validation statistics. At forecast initialization, SMAP is much too salty within 10o of the equator, leading to deeper MLD east of 165W. This deeper MLD leads to over-damping of the downwelling signal (i.e. relative upwelling), in turn leading to relatively too cool ENSO forecasts.

Hackert, E.↗
NASA NTRS2002

Self-Consistent Magnetosphere-Ionosphere Coupling and Associated Plasma Energization Processes

Magnetosphere-Ionosphere (MI) coupling and associated with this process electron and ion energization processes have interested scientists for decades and, in spite of experimental and theoretical research efforts, are still ones of the least well known dynamic processes in space plasma physics. The reason for this is that the numerous physical processes associated with MI coupling occur over multiple spatial lengths and temporal scales. One typical example of MI coupling is large scale ring current (RC) electrodynamic coupling that includes calculation of the magnetospheric electric field that is consistent with the ring current (RC) distribution. A general scheme for numerical simulation of such large-scale magnetosphere-ionosphere coupling processes has been presented earlier in many works. The mathematical formulation of these models are based on "modified frozen-in flux theorem" for an ensemble of adiabatically drifting particles in the magnetosphere. By tracking the flow of particles through the inner magnetosphere, the bounce-averaged phase space density of the hot ions and electrons can be reconstructed and the magnetospheric electric field can be calculated such that it is consistent with the particle distribution in the magnetosphere. The new a self-consistent ring current model has been developed that couples electron and ion magnetospheric dynamics with calculation of electric field. Two new features were taken into account in addition to the RC ions, we solve an electron kinetic equation in our model, self-consistently including these results in the solution. Second, using different analytical relationships, we calculate the height integrated ionospheric conductances as the function of precipitated high energy magnetospheric electrons and ions as produced by our model. This results in fundamental changes to the electric potential pattern in the inner magnetosphere, with a smaller Alfven boundary than previous potential formulations would predict but one consistent with recent satellite observations. This leads to deeper penetration of the plasma sheet ions and electrons into the inner magnetosphere and more effective ring current ions and electron energization.

Khazanov, G. V.↗
DOE OSTI2024

On learning what to learn: Heterogeneous observations of dynamics and establishing possibly causal relations among them

Abstract Before we attempt to (approximately) learn a function between two sets of observables of a physical process, we must first decide what the inputs and outputs of the desired function are going to be. Here we demonstrate two distinct, data-driven ways of first deciding “the right quantities” to relate through such a function, and then proceeding to learn it. This is accomplished by first processing simultaneous heterogeneous data streams (ensembles of time series) from observations of a physical system: records of multiple observation processes of the system. We determine (i) what subsets of observables are common between the observation processes (and therefore observable from each other, relatable through a function); and (ii) what information is unrelated to these common observables, therefore particular to each observation process, and not contributing to the desired function. Any data-driven technique can subsequently be used to learn the input–output relation—from k-nearest neighbors and Geometric Harmonics to Gaussian Processes and Neural Networks. Two particular “twists” of the approach are discussed. The first has to do with the identifiability of particular quantities of interest from the measurements. We now construct mappings from a single set of observations from one process to entire level sets of measurements of the second process, consistent with this single set. The second attempts to relate our framework to a form of causality: if one of the observation processes measures “now,” while the second observation process measures “in the future,” the function to be learned among what is common across observation processes constitutes a dynamical model for the system evolution.

Sroczynski, David W.↗
DOE OSTI2024

Dataset for manuscript "Consequences of the failure of equipartition for the p-V behavior of liquid water and the hydration free energy components of a small protein"

Previously, we showed that in the molecular dynamics simulation of a rigid model of water it is necessary to use an integration time-step dt that is less than or equal to 0.5 fs to ensure equipartition between translational and rotational modes. We extended that study in the NVT ensemble to NpT conditions and to an aqueous protein. We study neat liquid water with the rigid, SPC/E model and the protein BBA (PDB ID: 1FME) solvated in the rigid, TIP3P model. We examined integration time-steps ranging from 0.5 fs to 4.0 fs for various thermostat plus barostat combinations. We find that a small time-step, dt, is necessary to ensure consistent prediction of the simulation volume. Hydrogen mass repartitioning alleviates the problem somewhat, but is ineffective for the typical time-step used with this approach. The compressibility, a measure of volume fluctuations, is seen to be sensitive to dt. Using the mean volume estimated from the NpT simulation, we examined the electrostatic and van der Waals contribution to the hydration free energy of the protein in the NVT ensemble. These contributions are also sensitive to dt. In going from a time-step of 2 fs to a time-step of 0.5 fs, the change in the net electrostatic plus van der Waals contribution to the hydration of BBA is already in excess of the folding free energy reported for this protein. The data-set contains the simulation metadata and log files that support the claims noted above.

59 BASIC BIOLOGICAL SCIENCES↗