Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ensemble simulation”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Tuning a RANS k-e model for jet-in-crossflow simulations.

We develop a novel calibration approach to address the problem of predictive k e RANS simulations of jet-incrossflow. Our approach is based on the hypothesis that predictive k e parameters can be obtained by estimating them from a strongly vortical flow, specifically, flow over a square cylinder. In this study, we estimate three k e parameters, C%CE%BC, Ce2 and Ce1 by fitting 2D RANS simulations to experimental data. We use polynomial surrogates of 2D RANS for this purpose. We conduct an ensemble of 2D RANS runs using samples of (C%CE%BC;Ce2;Ce1) and regress Reynolds stresses to the samples using a simple polynomial. We then use this surrogate of the 2D RANS model to infer a joint distribution for the k e parameters by solving a Bayesian inverse problem, conditioned on the experimental data. The calibrated (C%CE%BC;Ce2;Ce1) distribution is used to seed an ensemble of 3D jet-in-crossflow simulations. We compare the ensemble's predictions of the flowfield, at two planes, to PIV measurements and estimate the predictive skill of the calibrated 3D RANS model. We also compare it against 3D RANS predictions using the nominal (uncalibrated) values of (C%CE%BC;Ce2;Ce1), and find that calibration delivers a significant improvement to the predictive skill of the 3D RANS model. We repeat the calibration using surrogate models based on kriging and find that the calibration, based on these more accurate models, is not much better that those obtained with simple polynomial surrogates. We discuss the reasons for this rather surprising outcome.

Lefantzi, Sophia↗
DOE OSTI2022

Nonane and Hexanol Adsorption in the Lamellar Phase of a Nonionic Surfactant: Molecular Simulations and Comparison to Ideal Adsorbed Solution Theory

Adsorption of n-nonane/1-hexanol (C9/C6OH) mixtures into the lamellar phase formed by a 50/50 w/w triethylene glycol mono-n-decyl ether (C10E3)/water system was studied using configurational-bias Monte Carlo simulations in the osmotic Gibbs ensemble. The interactions were described by the Shinoda–Devane–Klein coarse-grained force field. Prior simulations probing single-component adsorption indicated that C9 molecules preferentially load near the center of the bilayer increasing the bilayer thickness, whereas C6OH molecules are more likely to be found near the interface of the polar and non-polar moieties swelling the bilayer in the lateral dimension. Here, we extend this work to binary C9/C6OH adsorption to probe whether the difference in the spatial preferences may lead to a synergistic effect and enhanced loadings for the mixture. Comparing loading trends and the thermodynamics of binary adsorption to unary adsorption reveals that C9–C9 interactions lead to the largest enhancement, whereas C9–C6OH and C6OH–C6OH interactions are less favorable for this bilayer system. As a result, ideal adsorbed solution theory yields satisfactory predictions of the binary loading.

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

A Deep Learning-Driven Sampling Technique to Explore the Phase Space of an RNA Stem-Loop

The folding and unfolding of RNA stem-loops are critical biological processes; however, their computational studies are often hampered by the ruggedness of their folding landscape, necessitating long simulation times at the atomistic scale. Here, we adapted DeepDriveMD (DDMD), an advanced deep learning-driven sampling technique originally developed for protein folding, to address the challenges of RNA stem-loop folding. Although tempering- and order parameter-based techniques are commonly used for similar rare-event problems, the computational costs or the need for a priori knowledge about the system often present a challenge in their effective use. DDMD overcomes these challenges by adaptively learning from an ensemble of running MD simulations using generic contact maps as the raw input. DeepDriveMD enables on-the-fly learning of a low-dimensional latent representation and guides the simulation toward the undersampled regions while optimizing the resources to explore the relevant parts of the phase space. We showed that DDMD estimates the free energy landscape of the RNA stem-loop reasonably well at room temperature. Our simulation framework runs at a constant temperature without external biasing potential, hence preserving the information on transition rates, with a computational cost much lower than that of the simulations performed with external biasing potentials. Here, we also introduced a reweighting strategy for obtaining unbiased free energy surfaces and presented a qualitative analysis of the latent space. This analysis showed that the latent space captures the relevant slow degrees of freedom for the RNA folding problem of interest. Finally, throughout the manuscript, we outlined how different parameters are selected and optimized to adapt DDMD for this system. We believe this compendium of decision-making processes will help new users adapt this technique for the rare-event sampling problems of their interest.

Gupta, Ayush↗
DOE OSTI2013

Role of collisional broadening in Monte Carlo simulations of terahertz quantum cascade lasers

Using a generalized version of Fermi's golden rule, collisional broadening is self-consistently implemented into ensemble Monte Carlo carrier transport simulations, and its effect on the transport and optical properties of terahertz quantum cascade lasers is investigated. The inclusion of broadening yields improved agreement with the experiment, without a significant increase of the numerical load. Specifically, this effect is crucial for a correct modeling at low biases. In the lasing regime, broadening can lead to significantly reduced optical gain and output power, affecting the obtained current-voltage characteristics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2003

Simulated Radar Characteristics of LBA Convective Systems: Easterly and Westerly Regimes

The 3D Goddard Cumulus Ensemble (GCE) model was used to simulate convection that occurred during the TRMM LBA field experiment in Brazil. Convection in this region can be categorized into two different regimes. Low-level easterly flow results in moderate to high CAPE and a drier environment. Convection is more intense like that seen over continents. Low-level westerly flow results in low CAPE and a moist environment. Convection is weaker and more widespread characteristic of oceanic or monsoon-like systems. The GCE model has been used to study both regimes n order to provide cloud datasets that are representative of both environments in support of TRMM rainfall and heating algorithm development. Two different cases are analyzed: Jan 26, 1999, an eastely regime case, and Feb 23, 1999, a westerly regime case. The Jan 26 case is an organized squall line, while the Feb 23 case is less organized with only transient lines. Radar signatures, including CFADs, from the two simulated cases are compared to each other and with observations. The microphysical processes simulated in the model are also compared between the two cases.

Lang, Stephen E.↗
DOE OSTI2014

Energy dissipation by whistler turbulence: Three-dimensional particle-in-cell simulations

Three-dimensional particle-in-cell simulations of whistler turbulence are carried out on a collisionless, homogeneous, magnetized plasma model. The simulations use an initial ensemble of relatively long wavelength whistler modes and follow the temporal evolution of the fluctuations as they cascade into a broadband, anisotropic, turbulent spectrum at shorter wavelengths. For relatively small levels of the initial fluctuation energy ϵ{sub e}, linear collisionless damping provides most of the dissipation of the turbulence. But as ϵ{sub e} and the total dissipation increase, linear damping becomes less important and, especially at β{sub e} ≪ 1, nonlinear processes become stronger. The PDFs and kurtoses of the magnetic field increments in the simulations suggest that intermittency in whistler turbulence generally increases with increasing ϵ{sub e} and β{sub e}. Correlation coefficient calculations imply that the current structure dissipation also increases with increasing ϵ{sub e} and β{sub e}, and that the nonlinear dissipation processes in these simulations are primarily associated with regions of localized current structures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Unraveling the implications of finite specimen size on the interpretation of dynamic experiments for polycrystalline aluminum through direct numerical simulations

Normal and Pressure-shear plate impact (NPI and PSPI) experiments are popular experimental techniques for studying the mean-field macroscopic behavior of polycrystalline metals under high-rate dynamic loading. However, since both configurations rely upon geometry for subjecting the specimen to high strain rates, these experiments often involve a limited specimen size. Moreover, because of the inherent heterogeneities present within polycrystalline metals, it is difficult to ascertain if the size of the specimen and/or regions where measurements are made are sufficiently large for making representative inferences about the mean-field macroscopic properties from single-point velocity measurements. In the present study, we quantify the expected measurement variability on observable point measurements in NPI and PSPI experiments by carrying out direct numerical simulations (DNS) of statistically representative polycrystalline microstructures subjected to dynamic compression and compression-shear loading. In particular, we consider the role of specific material heterogeneities (e.g. the grain-to-grain difference in size, crystallographic orientation) on dispersion in the normal and transverse particle velocity records and on local fluctuations in key state variables (e.g. velocity, accumulated plastic strain) by incorporating these effects directly into a representative synthetic microstructure geometry and crystalline description of pure polycrystalline aluminum. The form of the present study is a large parametric investigation, consisting of ten ensembles of one hundred simulations. Each of the thousand simulations reflects a randomly realized synthetic microstructure in one of five cases of decreasing average grain size for the two loading configurations. Our analysis of the DNS results demonstrates that for both of these experimental configurations, the grain size directly correlates with the coefficient of variation (CV) in simulated point measurements, showing a convergent decrease in CV to zero (i.e. particle velocity record approaches the mean-field value) with decreasing grain size. Remarkably, the magnitude of variations in the particle velocity record is shown to be largest where the deviatoric stresses are most significant. In the case of NPI, this occurs at the elastic and plastic wavefront, whereas, in the case of PSPI, the magnitude of fluctuations are approximately constant throughout the experimental window time. The reasoning for the scatter in particle velocity due to the heterogeneous microstructure is demonstrated to be dependent on the mechanisms for accommodating deformation and on the interaction of reflection waves generated at sites of heterogeneities occurring at the scale of grains. Lastly, we develop a power-law description for the magnitude of scattering versus characteristic length, which provides a statistical framework for assessing the required number of grains per characteristic specimen dimension for minimizing scatter within these two experimental configurations (NPI, PSPI).

36 MATERIALS SCIENCE↗
DOE OSTI2015

Structure and dynamics of POPC bilayers in water solutions of room temperature ionic liquids

Molecular dynamics simulations in the NPT ensemble have been carried out to investigate the effect of two room temperature ionic liquids (RTILs), on stacks of phospholipid bilayers in water. We consider RTIL compounds consisting of chloride ([bmim][Cl]) and hexafluorophosphate ([bmim][PF{sub 6}]) salts of the 1-buthyl-3-methylimidazolium ([bmim]{sup +}) cation, while the phospholipid bilayer is made of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC). Our investigations focus on structural and dynamical properties of phospholipid and water molecules that could be probed by inelastic and quasi-elastic neutron scattering measurements. The results confirm the fast incorporation of [bmim]{sup +} into the lipid phase already observed in previous simulations, driven by the Coulomb attraction of the cation for the most electronegative oxygens in the POPC head group and by sizeable dispersion forces binding the neutral hydrocarbon tails of [bmim]{sup +} and of POPC. The [bmim]{sup +} absorption into the bilayer favours the penetration of water into POPC, causes a slight but systematic thinning of the bilayer, and further stabilises hydrogen bonds at the lipid/water interface that already in pure samples (no RTIL) display a lifetime much longer than in bulk water. On the other hand, the effect of RTILs on the diffusion constant of POPC (D{sub POPC}) does not reveal a clearly identifiable trend, since D{sub POPC} increases upon addition of [bmim][Cl] and decreases in the [bmim][PF{sub 6}] case. Moreover, because of screening, the electrostatic signature of each bilayer is only moderately affected by the addition of RTIL ions in solution. The analysis of long wavelength fluctuations of the bilayers shows that RTIL sorption causes a general decrease of the lipid/water interfacial tension and bending rigidity, pointing to the destabilizing effect of RTILs on lipid bilayers.

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

Clay Swelling in Dry Supercritical Carbon Dioxide: Effects of Interlayer Cations on the Structure, Dynamics, and Energetics of CO 2 Intercalation Probed by XRD, NMR, and GCMD Simulations

In situ X-ray diffraction (XRD) and nuclear magnetic resonance (NMR) experiments combined with molecular dynamics simulations using the grand canonical ensemble [grand canonical molecular dynamics (GCMD)] show that the cation size, charge, and solvation energy play critical roles in determining the interlayer expansion of smectite clay minerals when exposed to dry supercritical scCO2 (scCO2) under conditions relevant to petroleum reservoirs and geological CO2 sequestration conditions. The GCMD results show that the smectite mineral, hectorite, containing interlayer alkali and alkaline earth cations with small ionic radii and high solvation energies (e.g., Na+ and Ca2+) does not intercalate CO2 and that the fully collapsed interlayer structure is the most energetically stable configuration. With Cs+ and Ba2+, the monolayer structure is the stable configuration, and CO2 should spontaneously enter the interlayer. With Cs+, there is not even an energy barrier for CO2 intercalation, in agreement with the XRD and NMR results. 13C NMR and simulations show that the average orientation of the intercalated CO2 is with their O-C-O axes parallel to the basal clay surface and that they undergo a rapid rotation about an axis perpendicular to the main molecular axis. The simulations show that the strength of the interaction between the exchangeable cation and the clay structure dominates the intercalation energetics in dry scCO2.

Loganathan, Narasimhan (ORCID:0000000302515430)↗
NASA NTRSDate not supplied

Improving the Representation of Land Surface Processes using the Data Assimilation Research Testbed (DART)

The land surface is a critical part of the earth system as processes related to water, carbon, energy and nitrogen cycling have important implications for climate forcing, air quality, water availability and seasonal atmospheric forecasting. Despite advances in land surface modeling, land surface model performance is often limited because of errors related to initial and boundary conditions, model structure, and parameters. Data assimilation (DA) techniques combined with an expanding network of earth system observations present an opportunity to reduce these errors and improve simulations. Here we apply an Ensemble Kalman Filter DA system as part of the Data Assimilation Research Testbed to a variety of land surface simulations. First, we describe the use of remotely sensed biomass observations to provide improved simulations of plant phenology, carbon and water cycling for regions highly sensitive to climate change (Western US, China, and Arctic). We discuss approaches to account for systemic biases between models and observations, including the use of spatially-varying adaptive ensemble inflation as an alternative approach to re-scaling soil moisture observations. Finally, we discuss a strategy to incorporate complementary observations (snow water equivalent, solar-induced fluorescence) to better constrain the representation of carbon and water cycling across complex terrain.

DART↗
DOE OSTI2016

Enhancing pairwise state-transition weights: A new weighting scheme in simulated tempering that can minimize transition time between a pair of conformational states

Simulated tempering (ST) is a widely used enhancing sampling method for Molecular Dynamics simulations. As one expanded ensemble method, ST is a combination of canonical ensembles at different temperatures and the acceptance probability of cross-temperature transitions is determined by both the temperature difference and the weights of each temperature. One popular way to obtain the weights is to adopt the free energy of each canonical ensemble, which achieves uniform sampling among temperature space. However, this uniform distribution in temperature space may not be optimal since high temperatures do not always speed up the conformational transitions of interest, as anti-Arrhenius kinetics are prevalent in protein and RNA folding. Here, we propose a new method: Enhancing Pairwise State-transition Weights (EPSW), to obtain the optimal weights by minimizing the round-trip time for transitions among different metastable states at the temperature of interest in ST. The novelty of the EPSW algorithm lies in explicitly considering the kinetics of conformation transitions when optimizing the weights of different temperatures. We further demonstrate the power of EPSW in three different systems: a simple two-temperature model, a two-dimensional model for protein folding with anti-Arrhenius kinetics, and the alanine dipeptide. The results from these three systems showed that the new algorithm can substantially accelerate the transitions between conformational states of interest in the ST expanded ensemble and further facilitate the convergence of thermodynamics compared to the widely used free energy weights. We anticipate that this algorithm is particularly useful for studying functional conformational changes of biological systems where the initial and final states are often known from structural biology experiments.

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

The 1996 Mid-Atlantic Winter Flood: Exploring Climate Risk through a Storyline Approach

Abstract This article explores the application of thermodynamic perturbations to a historical midlatitude, wintertime, rain-on-snow flood event to evaluate how similar events may evolve under different climate forcings. In particular, we generate a hindcast of the 1996 Mid-Atlantic flood using an ensemble of 14-km variable-resolution simulations completed with the U.S. Department of Energy’s global Energy Exascale Earth System Model (E3SM). We show that the event is skillfully reproduced over the Susquehanna River Basin (SRB) by E3SM when benchmarked against in situ observational data and high-resolution reanalyses. In addition, we perform five counterfactual experiments to simulate the flood under preindustrial conditions and four different levels of warming as projected by the Community Earth System Model Large Ensemble. We find a nonlinear response in simulated surface runoff and streamflow as a function of atmospheric warming. This is attributed to changing contributions of liquid water input from a shallower initial snowpack (decreased snowmelt), increased surface temperatures and rainfall rates, and increased soil water storage. Flooding associated with this event peaks from around +1 to +2 K of global average surface warming and decreases with additional warming beyond this. There are noticeable timing shifts in peak runoff and streamflow associated with changes in the flashiness of the event. This work highlights the utility of using storyline approaches for communicating climate risk and demonstrates the potential nonlinearities associated with hydrologic extremes in areas that experience ephemeral snowpack, such as the SRB.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

Simulation of ENSO Teleconnections to Precipitation Extremes over the United States in the High-Resolution Version of E3SM

Abstract We evaluate the simulated teleconnection of El Niño–Southern Oscillation (ENSO) to winter season precipitation extremes over the United States in a long (98 years) 1950 control high-resolution version (HR; 25-km nominal atmosphere model horizontal resolution) of the U.S. Department of Energy’s (DOE) Energy Exascale Earth System Model version 1 (E3SMv1). The model bias and spatial pattern of ENSO teleconnections to mean and extreme precipitation in HR overall are similar to the low-resolution model’s (LR; 110 km) historical simulation (four-member ensemble, 1925–59). However, over the southeastern United States (SE-U.S.), HR produces stronger El Niño–associated extremes, reducing LR’s model bias. Both LR and HR produce weaker than observed increase in storm track activity during El Niño events there, but HR improves the ENSO-associated variability of moisture transport over SE-U.S. During El Niño, stronger vertical velocities in HR produce stronger large-scale precipitation, causing larger latent heating of the troposphere that pulls in more moisture from the Gulf of Mexico into the SE-U.S. This positive feedback also contributes to the stronger mean and extreme precipitation response in HR. Over the Pacific Northwest, LR’s bias of stronger than observed La Niña associated extremes is amplified in HR. Both models simulate stronger than observed moisture transport from the Pacific Ocean into the region during La Niña years. The amplified HR bias there is due to stronger orographically driven vertical updrafts that create stronger large-scale precipitation, despite weaker La Niña–induced storm track activity. Significance Statement New high-resolution Earth system models (ESMs) solve mathematical equations of fluid flow at much smaller spatial scales than prevalent ESMs, and thus are prohibitively expensive to compute. However, they can be useful for simulating accurate details of regional climate extremes that are driven by naturally occurring climate oscillations like El Niño–Southern Oscillation (ENSO). Here, we evaluate the simulation of ENSO-driven precipitation extremes over the United States in the high-resolution version of the U.S. Department of Energy’s new Energy Exascale Earth System Model version 1. We find that the high-resolution model improves upon its low-resolution counterpart over the southeastern United States by producing a better transport of moisture into the region from the Gulf of Mexico during El Niño. Over the U.S. Pacific Northwest, the high-resolution model simulates the atmospheric flow in more detail over the complex mountainous terrain. However, it also brings in more moisture from the Pacific Ocean just like the low-resolution model. This causes it to produce precipitation extremes during La Niña years there that are stronger than that observed in the real world.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2016

The Effect of Future Ambient Air Pollution on Human Premature Mortality to 2100 Using Output from the ACCMIP Model Ensemble

Ambient air pollution from ground-level ozone and fine particulate matter (PM(sub 2.5)) is associated with premature mortality. Future concentrations of these air pollutants will be driven by natural and anthropogenic emissions and by climate change. Using anthropogenic and biomass burning emissions projected in the four Representative Concentration Pathway scenarios (RCPs), the ACCMIP ensemble of chemistry climate models simulated future concentrations of ozone and PM(sub 2.5) at selected decades between 2000 and 2100. We use output from the ACCMIP ensemble, together with projections of future population and baseline mortality rates, to quantify the human premature mortality impacts of future ambient air pollution. Future air-pollution-related premature mortality in 2030, 2050 and 2100 is estimated for each scenario and for each model using a health impact function based on changes in concentrations of ozone and PM(sub 2.5) relative to 2000 and projected future population and baseline mortality rates. Additionally, the global mortality burden of ozone and PM(sub 2.5) in 2000 and each future period is estimated relative to 1850 concentrations, using present-day and future population and baseline mortality rates. The change in future ozone concentrations relative to 2000 is associated with excess global premature mortality in some scenarios/periods, particularly in RCP8.5 in 2100 (316 thousand deaths per year), likely driven by the large increase in methane emissions and by the net effect of climate change projected in this scenario, but it leads to considerable avoided premature mortality for the three other RCPs. However, the global mortality burden of ozone markedly increases from 382000 (121000 to 728000) deaths per year in 2000 to between 1.09 and 2.36 million deaths per year in 2100, across RCPs, mostly due to the effect of increases in population and baseline mortality rates. PM(sub 2.5) concentrations decrease relative to 2000 in all scenarios, due to projected reductions in emissions, and are associated with avoided premature mortality, particularly in 2100: between 2.39 and 1.31 million deaths per year for the four RCPs. The global mortality burden of PM(sub 2.5) is estimated to decrease from 1.70 (1.30 to 2.10) million deaths per year in 2000 to between 0.95 and 1.55 million deaths per year in 2100 for the four RCPs due to the combined effect of decreases in PM(sub 2.5) concentrations and changes in population and baseline mortality rates. Trends in future air-pollution-related mortality vary regionally across scenarios, reflecting assumptions for economic growth and air pollution control specific to each RCP and region. Mortality estimates differ among chemistry climate models due to differences in simulated pollutant concentrations, which is the greatest contributor to overall mortality uncertainty for most cases assessed here, supporting the use of model ensembles to characterize uncertainty. Increases in exposed population and baseline mortality rates of respiratory diseases magnify the impact on premature mortality of changes in future air pollutant concentrations and explain why the future global mortality burden of air pollution can exceed the current burden, even where air pollutant concentrations decrease.

Silva, Raquel A.↗
NASA NTRS2022

Model Evaluation of Short-Lived Climate Forcers for the Arctic Monitoring and Assessment Programme: A Multi-Species, Multi-Model Study

While carbon dioxide is the main cause for global warming, modeling short-lived climate forcers (SLCFs) such as methane, ozone, and particles in the Arctic allows us to simulate near-term climate and health impacts for a sensitive, pristine region that is warming at 3 times the global rate. Atmospheric modeling is critical for understanding the long-range transport of pollutants to the Arctic, as well as the abundance and distribution of SLCFs throughout the Arctic atmosphere. Modeling is also used as a tool to determine SLCF impacts on climate and health in the present and in future emissions scenarios. In this study, we evaluate 18 state-of-the-art atmospheric and Earth system models by assessing their representation of Arctic and Northern Hemisphere atmospheric SLCF distributions, considering a wide range of different chemical species (methane, tropospheric ozone and its precursors, black carbon, sulfate, organic aerosol, and particulate matter) and multiple observational datasets. Model simulations over 4 years (2008-2009 and 2014-2015) conducted for the 2022 Arctic Monitoring and Assessment Programme (AMAP) SLCF assessment report are thoroughly evaluated against satellite, ground, ship, and aircraft-based observations. The annual means, seasonal cycles, and 3-D distributions of SLCFs were evaluated using several metrics, such as absolute and percent model biases and correlation coefficients. The results show a large range in model performance, with no one particular model or model type performing well for all regions and all SLCF species. The multi-model mean (mmm) was able to represent the general features of SLCFs in the Arctic and had the best overall performance. For the SLCFs with the greatest radiative impact (CH4, 03, BC, and SO(sup 2-)(sub 4)), the mmm was within ±25 % of the measurements across the Northern Hemisphere. Therefore, we recommend a multi-model ensemble be used for simulating climate and health impacts of SLCFs. Of the SLCFs in our study, model biases were smallest for C"4 and greatest for OA. For most SLCFs, model biases skewed from positive to negative with increasing latitude. Our analysis suggests that vertical mixing, long-range transport, deposition, and wildfires remain highly uncertain processes. These processes need better representation within atmospheric models to improve their simulation of SLCFs in the Arctic environment. As model development proceeds in these areas, we highly recommend that the vertical and 3-D distribution of SLCFs be evaluated, as that information is critical to improving the uncertain processes in models.

Arctic Monitoring and Assessment Programme↗
DOE OSTI2017

First-Principles Molecular Dynamics Study of a Deep Eutectic Solvent: Choline Chloride/Urea and Its Mixture with Water

First-principles molecular dynamics simulations in the canonical ensemble at temperatures of 333 and 363 K and at the corresponding experimental densities are carried out to investigate the behavior of the 1:2 choline chloride/urea (reline) deep eutectic solvent and its equimolar mixture with water. Analysis of atom–atom radial and spatial distribution functions and of the H-bond network reveals the microheterogeneous structure of these complex liquid mixtures. In neat reline, the structure is governed by strong H-bonds of the trans- and cis-H atoms of urea to the chloride ion. In hydrous reline, water competes for the anions, and the H atoms of urea have similar propensities to bond to the chloride ions and the O atoms of urea and water. Finally, the vibrational spectra exhibit relatively broad peaks reflecting the heterogeneity of the environment. Although the 100 ps trajectories allow only for a qualitative assessment of transport properties, the simulations indicate that water is more mobile than the other species and its addition also fosters faster motion of urea.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Conditional channel simulation

In this work we design a specific simulation tool for quantum channels which is based on the use of a control system. This allows us to simulate an average quantum channel which is expressed in terms of an ensemble of channels, even when these channel-components are not jointly teleportation-covariant. This design is also extended to asymptotic simulations, continuous ensembles, and memory channels. As an application, we derive relative-entropy-of-entanglement upper bounds for private communication over various channels, including non-Gaussian mixtures of bosonic lossy channels. Among other results, we also establish the two-way quantum and private capacity of the so-called “dephrasure” channel.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

The Climate Response to Emissions Reductions Due to COVID-19: Initial Results From CovidMIP

Many nations responded to the COVID-19 pandemic by restricting travel and other activities during 2020, resulting in temporarily reduced emissions of CO 2 , other greenhouse gases and ozone and aerosol precursors. We perform a coordinated Intercomparison, CovidMIP, of Earth System model simulations to assess the impact on climate of these emissions reductions. Eleven models performed multiple initial-condition ensembles to produce over 280 simulations spanning both initial condition and model structural uncertainty. We find model consensus on reduced aerosol amounts (particularly over East Asia) and associated increases in surface shortwave radiation levels. However, any impact on near-surface temperature or rainfall during 2020-2024 is extremely small and is not detectable in this initial analysis. Regional analyses on a finer scale, and closer attention to extremes (especially linked to changes in atmospheric composition and air quality) are required to test the impact of COVID-19-related emission reductions on near-term climate.

54 ENVIRONMENTAL SCIENCES↗