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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Stochastic symplectic reduced-order modeling for model-form uncertainty quantification in molecular dynamics simulations in various statistical ensembles

Here, this work focuses on the representation of model-form uncertainties in molecular dynamics simulations in various statistical ensembles. In prior contributions, the modeling of such uncertainties was formalized and applied to quantify the impact of, and the error generated by, pair-potential selection in the microcanonical ensemble (NVE). In this work, we extend this formulation and present a linear-subspace reduced-order model for the canonical (NVT) and isobaric (NPT) ensembles. The symplectic reduced-order basis is randomized on the tangent space of the Stiefel manifold to provide topological relationships and capture model-form uncertainty. Using the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS), we assess the relevance of these stochastic reduced-order atomistic models on canonical problems involving a Lennard-Jones fluid and an argon crystal melt.

42 ENGINEERING↗
DOE OSTI2025

Entropy is an important design principle in the photosystem II supercomplex

Photosystem II (PSII) can achieve near-unity quantum efficiency of light harvesting in ideal conditions and can dissipate excess light energy as heat to prevent the formation of reactive oxygen species (ROS) under light stress. Understanding how this pigment–protein complex accomplishes these opposing goals is a topic of great interest that has so far been explored primarily through the lens of the system energetics. Despite PSII’s known flat energy landscape, a thorough consideration of the entropic effects on energy transfer in PSII is lacking. In this work, we aim to discern the free energetic design principles underlying the PSII energy transfer network. To accomplish this goal, we employ a structure-based rate matrix and compute the free energy terms in time following a specific initial excitation to discern how entropy and enthalpy drive ensemble system dynamics. We find that the interplay between the entropy and enthalpy components differ among each protein subunit, which allows each subunit to fulfill a unique role in the energy transfer network. This individuality ensures that PSII can accomplish efficient energy trapping in the reaction center (RC), effective nonphotochemical quenching (NPQ) in the periphery, and robust energy trapping in the other-monomer RC if the same-monomer RC is closed. We also show that entropy, in particular, is a dynamically tunable feature of the PSII free energy landscape accomplished through regulation of LHCII binding. These findings help rationalize natural photosynthesis and provide design principles for more efficient solar energy harvesting technologies.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1976

Dispersion of ensembles of non-interacting particles

The dynamics of an ensemble of noninteracting particles dispersing from a common origin and moving in a common force field with an initial distribution of momenta is analyzed using an approach where the particles are considered as a continuum described by a phase-space distribution function. General solutions are obtained for both the distribution function and the associated spatial density function. The linear case of small departures from circular orbits in an axisymmetric gravitational field is treated along with the specific case of particle dispersion from an object in a circular orbit in the same type of field. Numerical results are presented for the latter case, and consideration is given to the inverse problem of determining the initial time and velocity distribution from knowledge of the ensemble structure at a later time. Explicit results are provided for the case of an ellipsoidal distribution of initial momenta, and a numerical procedure is indicated for treating more general cases.

Heard, W. B.↗
DOE OSTI2025

Uncertainty of 21st Century western U.S. snowfall loss derived from regional climate model large ensemble

Abstract The western United States is dependent on winter snowfall over its major mountain ranges, which gradually melts each year, serving as a natural reservoir for water resources. In a future warmer climate, much of this snowfall could be replaced by rain, making it more challenging to capture and store water. In this study, we utilize an ensemble of dynamically downscaled simulations forced by 14 global climate models (GCMs). These GCMs project wildly different futures, in terms of both temperature and precipitation change, producing significant uncertainty in snowfall projections. Here we exploit the robust statistics of the downscaled ensemble, and diagose the sensitivity of end-of-century snowfall loss across the region to both warming and regional wetting/drying in the driving GCM. The windward slopes of the Sierra Nevada and Cascades are particularly sensitive to warming (losing ~ 15% annual snowfall per degree warming), with little influence of precipitation. By contrast, snowfall loss in the inter-mountain west is less sensitive to warming (~ 5% K −1 ), but is significantly offset/exacerbated by precipitation changes (~ 0.5% snow per 1% precipitation). Combining such sensitivities with the warming and regional precipitation signals in the full CMIP6 ensemble, we can fully quantify likely snowfall loss and its uncertainty at any location, for any emissions scenario. We find that the western U.S. as a whole will lose 34 ± 8% of its total volumetric snowfall by end-of-century under the high-emissions SSP3-7.0 scenario, but 25 ± 6% and 17 ± 6% under the lower-emissions SSP2-4.5 and SSP1-2.6 scenarios.

Norris, Jesse (ORCID:0000000289883326)↗
DOE OSTI2025

Dynamically downscaled seasonal heat wave projections in the CONUS

Heat waves are a well-documented hazard that are projected to increase in intensity, duration, and frequency with climate change. Regions of the US experience widely varying temperatures; for example, 35 °C is extremely hot for spring in the Northeast but not for summer in the Southeast. It is important to evaluate projections within a regional context and at a high enough resolution to understand the risks to populations. We identify heat waves across the Conterminous US (CONUS) under SSP5–8.5 from 2020 to 2059 with an ensemble of dynamically downscaled Coupled Model Intercomparison Project Phase 6 (CMIP6) model outputs. We demonstrate that there are regional differences caused by seasonal and local drivers of persistent hot temperatures. Summer heat waves are increasing in intensity and duration faster than winter heat waves because of the atmospheric conditions that promote these events. Our analysis emphasizes the value of fine-resolution modeling for projecting future climate risks.

Rubin, Hannah [University of Tennessee, Knoxville ↗
DOE OSTI2024

Adjusting the Energy Profile for CH–O Interactions Leads to Improved Stability of RNA Stem-Loop Structures in MD Simulations

The role of ribonucleic acid (RNA) in biology continues to grow, but insight into important aspects of RNA behavior is lacking, such as dynamic structural ensembles in different environments, how flexibility is coupled to function, and how function might be modulated by small molecule binding. In the case of proteins, much progress in these areas has been made by complementing experiments with atomistic simulations, but RNA simulation methods and force fields are less mature. It remains challenging to generate stable RNA simulations, even for small systems where well-defined, thermostable structures have been established by experiments. Further many different aspects of RNA energetics have been adjusted in force fields, seeking improvements that are transferable across a variety of RNA structural motifs. In this work, the role of weak CH···O interactions is explored, which are ubiquitous in RNA structure but have received less attention in RNA force field development. By comparing data extracted from high-resolution RNA crystal structures to energy profiles from quantum mechanics and force field calculations, it is shown that CH···O interactions are overly repulsive in the widely used Amber RNA force fields. A simple, targeted adjustment of CH···O repulsion that leaves the remainder of the force field unchanged was developed. Then, the standard and modified force fields were tested using molecular dynamics (MD) simulations with explicit water and salt, amassing over 300 μs of data for multiple RNA systems containing important features such as the presence of loops, base stacking interactions as well as canonical and noncanonical base pairing. In this work and others, standard force fields lead to reproducible unfolding of the NMR-based structures. Including a targeted CH···O adjustment in an otherwise identical protocol dramatically improves the outcome, leading to stable simulations for all RNA systems tested.

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

Examining the Robustness of Weakened Orographic Influence on Precipitation in Downscaled Climate Projections Over the Western US

Assessing local climate change impacts often requires downscaling coarse global climate model (GCM) output to finer resolution. Two main approaches exist: dynamical downscaling using high-resolution regional climate models, and statistical downscaling based on historical relationships between large-scale and local variables. In a recent analysis of five dynamically downscaled simulations over the western United States, Koszuta et al. (2024, https://doi.org/10.1029/2023gl107298) found that warming weakens orographic influence on winter precipitation, damping increases on windward slopes and amplifying them in rain-shadowed regions. Here we show that this effect is robust across seasons and multiple dynamically downscaled ensembles, and is more pronounced at higher model resolutions. However, it is absent in projections from a widely used statistical model (LOCA2), even when trained on high-resolution future simulations (LOCA2-Hybrid). This highlights a key limitation of many statistical downscaling methods: their preservation of parent GCM trends, which usually fail to capture emergent changes in orographic precipitation patterns.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Trends and Drivers of Terrestrial Sources and Sinks of Carbon Dioxide: An Overview of the TRENDY Project

The terrestrial biosphere plays a major role in the global carbon cycle, and there is a recognized need for regularly updated estimates of land-atmosphere exchange at regional and global scales. An international ensemble of Dynamic Global Vegetation Models (DGVMs), known as the “Trends and drivers of the regional scale terrestrial sources and sinks of carbon dioxide” (TRENDY) project, quantifies land biophysical exchange processes and biogeochemistry cycles in support of the annual Global Carbon Budget assessments and the REgional Carbon Cycle Assessment and Processes, phase 2 project. DGVMs use a common protocol and set of driving data sets. A set of factorial simulations allows attribution of spatio-temporal changes in land surface processes to three primary global change drivers: changes in atmospheric CO 2 , climate change and variability, and Land Use and Land Cover Changes (LULCC). Here, we describe the TRENDY project, benchmark DGVM performance using remote-sensing and other observational data, and present results for the contemporary period. Simulation results show a large global carbon sink in natural vegetation over 2012–2021, attributed to the CO 2 fertilization effect (3.8 ± 0.8 PgC/yr) and climate (–0.58 ± 0.54 PgC/yr). Forests and semi-arid ecosystems contribute approximately equally to the mean and trend in the natural land sink, and semi-arid ecosystems continue to dominate interannual variability. The natural sink is offset by net emissions from LULCC (–1.6 ± 0.5 PgC/yr), with a net land sink of 1.7 ± 0.6 PgC/yr. Despite the largest gross fluxes being in the tropics, the largest net land-atmosphere exchange is simulated in the extratropical regions.

58 GEOSCIENCES↗
DOE OSTI2025

Vegetation biogeography is a main source of uncertainty in modelling the land carbon cycle

The terrestrial biosphere exchanges a large amount of CO 2 with the atmosphere through photosynthesis and respiration, determining the magnitude of land carbon sink and consequently influencing the rate of global warming. The magnitudes of global photosynthesis and respiration, however, vary widely across models (100-200 PgC/year), constituting a key and persistent source of uncertainty in carbon cycle and climate modelling. Here, we argue that the uncertainty in the land carbon cycle modelling is largely attributable to the uncertainty in biogeography – the distribution of plant functional types (PFTs). Using an ensemble of dynamic global vegetation models (DGVMs), we find a strong dependence of total photosynthesis on total area for each PFT. The dependence allows us to reduce the spread of land carbon cycle estimates by ~75% using remote sensing-based PFT maps. We further find that 56 ± 21% of climate-driven changes in global photosynthesis modelled by DGVMs are caused by changes in PFT distribution in the last two decades. Our study identifies vegetation biogeography as a main controlling factor of uncertainty in land carbon cycle modelling and highlights the importance of biogeography-climate interactions in carbon cycle and climate studies.

Zhao, Ruiying [National Univ. of Singapore (Singap↗
DOE OSTI2025

Upper bounds for 21st-century surface air temperatures in the Western United States

The last decade has seen a large number of severe heatwaves that were unprecedented in the observational record, highlighting challenges associated with observationally-based statistical quantification of the likelihood and magnitude of future extreme temperatures. An alternative to such probabilistic assessments is identification of upper bounds that quantify the hottest surface air temperatures that can possibly be achieved by the end of the 21st century. Theory, simulations, and observational analyses support the existence of a finite upper bound for surface air temperature; however, estimates for future upper-bound values that are realistic and usable for planning remain unavailable. Here, we combine atmospheric theory with large ensembles of dynamically downscaled projections to estimate historical and end-of-century upper bounds for surface air temperatures. A number of physical mechanisms can influence upper bounds, and at the end of the 21st century, estimates based on mechanisms that yield more moderate upper-bounds produce values around 60∘C for much of the Western United States and in excess of 80∘C for the hottest parts of the domain. Even cooler high-altitude locations have end-of-century upper bounds over 50∘C. Although these upper-bound estimates might seem implausibly large, increases in the upper bounds over the 21st century are similar to increases in dynamically downscaled peak surface temperatures after adjusting those downscaled temperatures to eliminate the possibly biased model trends in surface specific humidity. While upper bound estimates are high relative to historical observations, they nonetheless suggest that heatwave intensity risk is bounded, with uncertainty dominated by projections of surface and upper-level humidity.

Risser, Mark D↗
NASA NTRS1996

Studying Turbulence Using Numerical Simulation Databases: Proceedings of the 1996 Summer Program - Part 6

Topics considered include: New approach to turbulence modeling; Second moment closure analysis of the backstep flow database; Prediction of the backflow and recovery regions in the backward facing step at various Reynolds numbers; Turbulent flame propagation in partially premixed flames; Ensemble averaged dynamic modeling. Also included a study of the turbulence structures of wall-bounded shear flows; Simulation and modeling of the elliptic streamline flow.

Source record↗
NASA NTRS2020

Causes of slowing‐down seasonal CO2 amplitude at Mauna Loa

Changing amplitude of the seasonal cycle of atmospheric CO2 (SCA) in the northern hemisphere is an emerging carbon cycle property. Mauna Loa (MLO) station (20°N, 156°W), which has the longest continuous northern hemisphere CO2 record, shows an increasing SCA before the 1980s (p < .01), followed by no significant change thereafter. We analyzed the potential driving factors of SCA slowing‐down, with an ensemble of dynamic global vegetation models (DGVMs) coupled with an atmospheric transport model. We found that slowing‐down of SCA at MLO is primarily explained by response of net biome productivity (NBP) to climate change, and by changes in atmospheric circulations. Through NBP, climate change increases SCA at MLO before the 1980s and decreases it afterwards. The effect of climate change on the slowing‐down of SCA at MLO is mainly exerted by intensified drought stress acting to offset the acceleration driven by CO2 fertilization. This challenges the view that CO2 fertilization is the dominant cause of emergent SCA trends at northern sites south of 40°N. The contribution of agricultural intensification on the deceleration of SCA at MLO was elusive according to land–atmosphere CO2 flux estimated by DGVMs and atmospheric inversions. Our results also show the necessity to adequately account for changing circulation patterns in understanding carbon cycle dynamics observed from atmospheric observations and in using these observations to benchmark DGVMs.

Mauna Loa↗
DOE OSTI2025

Solid-state platform for cooperative quantum dynamics driven by correlated emission

While traditionally regarded as an obstacle to quantum coherence, recent breakthroughs in quantum optics have shown that the dissipative interaction of a qubit with its environment can be leveraged to protect quantum states and synthesize many-body entanglement. Inspired by this progress, here we set the stage for the—as yet uncharted—exploration of analogous cooperative phenomena in hybrid solid-state platforms. We develop a comprehensive formalism for the quantum many-body dynamics of an ensemble of solid-state spin defects interacting with the magnetic field fluctuations of a common solid-state reservoir. Our framework applies to any solid-state reservoir whose fluctuating spin, pseudospin, or charge degrees of freedom generate magnetic fields. To understand whether correlations induced by dissipative processes can play a relevant role in a realistic experimental setup, we apply our model to a qubit array interacting via the spin fluctuations of a ferromagnetic bath. Our results show that the low-temperature collective relaxation rates of the qubit ensemble can display clear signatures of super- and subradiance, i.e., forms of cooperative dynamics traditionally achieved in atomic ensembles. We find that the solid-state analog of these cooperative phenomena is robust against spatial disorder in the qubit ensemble and thermal fluctuations of the magnetic reservoir, providing a route for their feasibility in near-term experiments. Finally, our work lays the foundation for a multiqubit approach to quantum sensing of solid-state systems and the direct generation of many-body entanglement in spin-defect ensembles. Furthermore, we discuss how the tunability of solid-state reservoirs opens up novel pathways for exploring cooperative phenomena in regimes beyond the reach of conventional quantum optics setups.

NV centers↗
NASA NTRS2009

Dynamics and Predictability of Hurricane Humberto (2007) Revealed from Ensemble Analysis and Forecasting

This study uses short-range ensemble forecasts initialized with an Ensemble-Kalman filter to study the dynamics and predictability of Hurricane Humberto, which made landfall along the Texas coast in 2007. Statistical correlation is used to determine why some ensemble members strengthen the incipient low into a hurricane and others do not. It is found that deep moisture and high convective available potential energy (CAPE) are two of the most important factors for the genesis of Humberto. Variations in CAPE result in as much difference (ensemble spread) in the final hurricane intensity as do variations in deep moisture. CAPE differences here are related to the interaction between the cyclone and a nearby front, which tends to stabilize the lower troposphere in the vicinity of the circulation center. This subsequently weakens convection and slows genesis. Eventually the wind-induced surface heat exchange mechanism and differences in landfall time result in even larger ensemble spread. 1

Sippel, Jason A.↗
DOE OSTI2024

RASPA3

RASPA3, a molecular simulation code for computing adsorption and diffusion in nanoporous materials and thermodynamic and transport properties of fluids. It implements force field based classical Monte Carlo/molecular dynamics in various ensembles. RASPA3 is rewritten from the ground up in C++23 with speed and code readability in mind. Transition-matrix Monte Carlo is added to compute the density of states and free energies. The Monte Carlo code for rigid molecules is based on quaternions, and the atomic positions needed in the energy evaluation are recreated from the center of mass position and quaternion orientation. The expanded ensemble methodology for fractional molecules, with a scaling parameter λ between 0 and 1, now also keeps track of analytic expressions of dU/dλ, allowing independent verification of the chemical potential using thermodynamic integration. The source code is freely available under the MIT license on GitHub.

Dubbeldam, David↗
NASA NTRS2001

Statistical Ensemble of Large Eddy Simulations

A statistical ensemble of large eddy simulations (LES) is run simultaneously for the same flow. The information provided by the different large scale velocity fields is used to propose an ensemble averaged version of the dynamic model. This produces local model parameters that only depend on the statistical properties of the flow. An important property of the ensemble averaged dynamic procedure is that it does not require any spatial averaging and can thus be used in fully inhomogeneous flows. Also, the ensemble of LES's provides statistics of the large scale velocity that can be used for building new models for the subgrid-scale stress tensor. The ensemble averaged dynamic procedure has been implemented with various models for three flows: decaying isotropic turbulence, forced isotropic turbulence, and the time developing plane wake. It is found that the results are almost independent of the number of LES's in the statistical ensemble provided that the ensemble contains at least 16 realizations.

Carati, Daniele↗
DOE OSTI2024

Coherent Phonon Dynamics in Plasmonic Gold Tetrahedral Nanoparticle Ensembles

Coherent phonon modes supported by plasmonic nanoparticles offer prospective applications in chemical and biological sensing. Whereas characterization of these phonon modes often requires single particle measurements, synthetic routes to narrow size distributions of nanoparticles permit ensemble investigations. Recently, the synthesis of highly monodisperse gold tetrahedral nanoparticles with tunable edge length and corner sharpness has been developed. Herein, we characterize a size series of these nanoparticles in colloidal dispersion via transient absorption spectroscopy to examine their mechanical and plasmonic responses upon photoexcitation. Oscillations of transient absorption signals are observed in the plasmon resonance and correspond to the lowest order radial breathing modes of the nanoparticles, the frequencies of which are affected by edge length and truncation of the corners. Furthermore, quality factor values ranging from 10 to 15 are observed for the oscillations that convey potential utility in mass-sensing and plasmon-exciton coupling photonics schemes. Finite-difference time domain and finite element analysis calculations establish the specific optically-relevant phonon modes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗