Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical thermodynamics”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Projection algorithm for state preparation on quantum computers

We present an efficient method to prepare states of a many-body system on quantum hardware, first isolating individual quantum numbers and then using time evolution to isolate the energy. Our method in its simplest form requires only one additional auxiliary qubit. The total time evolved for an accurate solution is proportional to the ratio of the spectrum range of the trial state to the gap to the lowest excited state, a substantial improvement over other projection algorithms, and the accuracy increases exponentially with the time evolved. Isolating the quantum numbers is efficient because of the known eigenvalues and increases the gap thus shortening the propagation time required. The success rate of the algorithm, or the probability of producing the desired state, is a simple function of measurement times and phases and is dominated by the square overlap of the original state to the desired state. We present examples from the nuclear shell model and the Heisenberg model. In conclusion, we compare this algorithm to previous algorithms for short evolution times and discuss potential further improvements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Application of the variational autoencoder to detect the critical points of the anisotropic Ising model

We generalize the previous study on the application of variational autoencoders to the two-dimensional Ising model to a system with anisotropy. Due to the self-duality property of the system, the critical points can be located exactly for the entire range of anisotropic coupling. This presents an excellent test bed for the validity of using a variational autoencoder to characterize an anisotropic classical model. Furthermore, we reproduce the phase diagram for a wide range of anisotropic couplings and temperatures via a variational autoencoder without the explicit construction of an order parameter. Considering that the partition function of ($d$ + 1)-dimensional anisotropic models can be mapped to that of the $d$-dimensional quantum spin models, the present study provides numerical evidence that a variational autoencoder can be applied to analyze quantum systems via the quantum Monte Carlo method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Dataset for manuscript "Equipartition and the temperature of maximum density of TIP4P/2005 water"

We simulate TIP4P/2005 water in the temperature range of 257 K to 318 K with time-steps 0.25, 0.50, 1.00, 2.00, and 4.00 fs. The density-temperature behavior obtained using 0.25 or 0.50 fs are in excellent agreement with each other but differ from those obtained using time-steps that have been shown earlier to lead to a breakdown of equipartition. The temperature of maximum density (TMD) is 277.15 K with time-step 0.25 or 0.50 fs, but is shifted to progressively lower values for longer time-steps, a trend that holds for different thermostat/barostat combinations. Enhancing the water-water dispersion interaction, as has been recommended for simulating disordered proteins in TIP4P/2005, degrades the description of the liquid-vapor phase envelope. We present a simple physically transparent reasoning to highlight the separation of the time-scales between translational and rotational motion. We also develop a metric, Chi, that we term the equipartition anomaly, to detect equipartition violations in simulations that include molecules that are treated as rigid objects. Calculating Chi is shown to be straightforward and sensitive to equipartition violations. A key takeaway from this study is that using sufficiently short time-steps (less than or equal to 0.5 fs) to preserve equipartition is essential for obtaining meaningful liquid water properties and for producing reliable simulation data, as correct-ensemble sampling is fundamental to ensure reproducibility across codes and simulation alogrithms. The included dataset provides the raw data used in the preparation of the graphs noted in the manuscript.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Single T gate in a Clifford circuit drives transition to universal entanglement spectrum statistics

Clifford circuits are insufficient for universal quantum computation or creating t t -designs with t\ge 4 t ≥ 4 . While the entanglement entropy is not a telltale of this insufficiency, the entanglement spectrum of a time evolved random product state is: the entanglement levels are Poisson-distributed for circuits restricted to the Clifford gate-set, while the levels follow Wigner-Dyson statistics when universal gates are used. In this paper we show, using finite-size scaling analysis of different measures of level spacing statistics, that in the thermodynamic limit, inserting a single T (\pi/8) ( π / 8 ) gate in the middle of a random Clifford circuit is sufficient to alter the entanglement spectrum from a Poisson to a Wigner-Dyson distribution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Bringing discrete-time Langevin splitting methods into agreement with thermodynamics

In light of the recently published complete set of statistically correct Grønbech–Jensen (GJ) methods for discrete-time thermodynamics, we revise a differential operator splitting method for the Langevin equation in order to comply with the basic GJ thermodynamic sampling features, namely, the Boltzmann distribution and Einstein diffusion, in linear systems. This revision, which is based on the introduction of time scaling along with flexibility of a discrete-time velocity attenuation parameter, provides a direct link between the ABO splitting formalism and the GJ methods. This link brings about the conclusion that any GJ method has at least weak second order accuracy in the applied time step. It further helps identify a novel half-step velocity, which simultaneously produces both correct kinetic statistics and correct transport measures for any of the statistically sound GJ methods. Explicit algorithmic expressions are given for the integration of the new half-step velocity into the GJ set of methods. Finally, numerical simulations, including quantum-based molecular dynamics (QMD) using the QMD suite Los Alamos Transferable Tight-Binding for Energetics, highlight the discussed properties of the algorithms as well as exhibit the direct application of robust, time-step-independent stochastic integrators to QMD.

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

The Relationship between Precipitation and Precipitable Water in CMIP6 Simulations and Implications for Tropical Climatology and Change

It is well documented that over the tropical oceans, column-integrated precipitable water (pw) and precipitation (P) have a nonlinear relationship. In this study moisture budget analysis is used to examine this P–pw relationship in a normalized precipitable water framework. It is shown that the parameters of the nonlinear relationship depend on the vertical structure of moisture convergence. Specifically, the precipitable water values at which precipitation is balanced independently by evaporation versus by moisture convergence define a critical normalized precipitable water, pw nc . This is a measure of convective inhibition that separates tropical precipitation into two regimes: a local evaporation-controlled regime with widespread drizzle and a precipitable water–controlled regime. Most of the 17 CMIP6 historical simulations examined here have higher pw nc compared to ERA5, and more frequently they operate in the drizzle regime. When compared to observations, they overestimate precipitation over the high-evaporation oceanic regions off the equator, thereby producing a ‘‘double ITCZ’’ feature, while underestimating precipitation over the large tropical landmasses and over the climatologically moist oceanic regions near the equator. The responses to warming under the SSP585 scenario are also examined using the normalized precipitable water framework. It is shown that the critical normalized precipitable water value at which evaporation versus moisture convergence balance precipitation decreases as a result of the competing dynamic and thermodynamic responses to warming, resulting in an increase in drizzle and total precipitation. Statistically significant historical trends corresponding to the thermodynamic and dynamic changes are detected in ERA5 and in lowintensity drizzle precipitation in the PERSIANN precipitation dataset.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Acceleration of Solvation Free Energy Calculation via Thermodynamic Integration Coupled with Gaussian Process Regression and Improved Gelman–Rubin Convergence Diagnostics

The determination of the solvation free energy of ions and molecules holds profound importance across a spectrum of applications spanning chemistry, biology, energy storage, and the environment. Molecular dynamics simulations are powerful tools for computing this critical parameter. Nevertheless, the accurate and efficient calculation of the solvation free energy becomes a formidable endeavor when dealing with complex systems characterized by potent Coulombic interactions and sluggish ion dynamics and, consequently, slow transition across various metastable states. Here, in the present study, we expose limitations stemming from the conventional calculation of the statistical inefficiency g in the thermodynamic integration method, a factor that can hinder the determination of convergence of the solvation free energy and its associated uncertainty. Instead, we propose a robust scheme based on Gelman–Rubin convergence diagnostics. We leverage this improved estimation of uncertainties to introduce an innovative accelerated thermodynamic integration method based on the Gaussian Process regression. This methodology is applied to the calculation of the solvation free energy of trivalent rare-earth elements immersed in ionic liquids, a scenario in which the aforementioned challenges render standard approaches ineffective. The proposed method proves to be effective in computing solvation free energy in situations where traditional thermodynamic integration methods fall short.

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

Elucidation of the physical factors that control activated transport of penetrants in chemically complex glass-forming liquids

Understanding the activated transport of penetrant or tracer atoms and molecules in condensed phases is a challenging problem in chemistry, materials science, physics, and biophysics. Many angstrom- and nanometer-scale features enter due to the highly variable shape, size, interaction, and conformational flexibility of the penetrant and matrix species, leading to a dramatic diversity of penetrant dynamics. Based on a minimalist model of a spherical penetrant in equilibrated dense matrices of hard spheres, a recent microscopic theory that relates hopping transport to local structure has predicted a novel correlation between penetrant diffusivity and the matrix thermodynamic dimensionless compressibility, S 0 ( T ) (which also quantifies the amplitude of long wavelength density fluctuations), as a consequence of a fundamental statistical mechanical relationship between structure and thermodynamics. Moreover, the penetrant activation barrier is predicted to have a factorized/multiplicative form, scaling as the product of an inverse power law of S 0 ( T ) and a linear/logarithmic function of the penetrant-to-matrix size ratio. This implies an enormous reduction in chemical complexity that is verified based solely on experimental data for diverse classes of chemically complex penetrants dissolved in molecular and polymeric liquids over a wide range of temperatures down to the kinetic glass transition. The predicted corollary that the penetrant diffusion constant decreases exponentially with inverse temperature raised to an exponent determined solely by how S 0 ( T ) decreases with cooling is also verified experimentally. Our findings are relevant to fundamental questions in glassy dynamics, self-averaging of angstrom-scale chemical features, and applications such as membrane separations, barrier coatings, drug delivery, and self-healing.

Science & Technology - Other Topics↗
DOE OSTI2021

Large-Scale Synoptic Systems and Fog During the C-FOG Field Experiment

Abstract The goal of this work is to summarize synoptic meteorological conditions during the Coastal Fog (C-FOG) field project that took place onshore and offshore of the Avalon Peninsula, Newfoundland, from 25 August until 8 October 2018. Visibility was measured at three locations at the Ferryland supersite that are about 1 km from each other, and at two additional sites 66 and 76 km to the north. Supporting meteorological measurements included surface winds, air temperature, humidity, pressure, radiation, cloud-base height, and atmospheric thermodynamic profiles from radiosonde soundings. Statistics are presented for surface measurements during fog events including turbulence kinetic energy, net longwave radiation, visibility, and precipitation. Eleven fog events are observed at Ferryland. Each significant fog event is related to a large-scale cyclonic system. The longest fog event is due to interaction of a northern deep low and a tropical cyclone. Fog occurrence is also examined across Atlantic Canada by including Sable Island, Yarmouth, Halifax, and Sydney. It is concluded that at Ferryland, all significant fog events occur under a cyclonic system while at Sable Island all significant fog events occur under both cyclonic and anticyclonic systems. The fog-formation mechanism involves cloud lowering and stratus broadening or only stratus broadening for the cyclonic systems while for the anticyclonic systems it is stratus broadening or radiation. Although widely cited as the main cause of fog in Atlantic Canada, advection fog is not found to be the primary or sole fog type in the events examined.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Effect of Collective Dynamics and Anharmonicity on Entropy in Heterogenous Catalysis: Building the Case for Advanced Molecular Simulations

In this paper, we present a perspective on the computational determination of entropy and its effects and consequences on heterogeneous catalysis. Special attention is paid to the role of anharmonicity and the deviations from the standard harmonic oscillator approximations which can fail to provide a reliable assessment of entropy. To address these challenges, advanced methodologies are needed that can explicitly account for these thermodynamic drivers through the appropriate statistical sampling of reactive free energy surfaces. We discuss where anharmonicity should be expected, where it has been observed from a theoretical perspective, and the methods currently employed to address it. We concentrate on three types of systems where we have observed major, nonnegligible anharmonic effects: (1) supported nanoparticles, where the migration of metal atoms, complexes, and entire clusters exhibit anharmonic behavior in their dynamic motion; (2) porous solids where confinement effects distort potential energy surfaces and hinder molecular motions resulting in large entropic terms; and (3) solid/liquid interfaces, where interactions between solvent molecules and adsorbed species can result in large solvent organization free energy and unique reactivity.

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

Thermodynamically consistent incorporation of the Langmuir adsorption model into compressible fluctuating hydrodynamics

For a gas–solid interfacial system where chemical species undergo reversible adsorption, we develop a mesoscopic stochastic modeling method that simulates both gas-phase hydrodynamics and surface coverage dynamics by coupling the Langmuir adsorption model with compressible fluctuating hydrodynamics. To this end, we derive a thermodynamically consistent mass–energy update scheme that accounts for how the mass and energy variables in the gas and surface subsystems should be updated according to the changes in the number of molecules of each species in each subsystem due to adsorption and desorption events. By performing a stochastic analysis for the ideal Langmuir model and the full hydrodynamic system, we analytically confirm that our mass–energy update scheme captures thermodynamic equilibrium predicted by equilibrium statistical mechanics. We find that an internal energy correction term is needed, which is attributed to the difference in the mean kinetic energy of gas molecules colliding with the surface from that computed from the Maxwell–Boltzmann distribution. By performing an equilibrium simulation study for an ideal gas mixture of CO and Ar, with CO undergoing reversible adsorption, we validate our overall simulation method and implementation.

Adsorption↗
DOE OSTI2021

Radiosonde Observations of Environments Supporting Deep Moist Convection Initiation during RELAMPAGO-CACTI

The Remote Sensing of Electrification, Lightning, and Mesoscale/Microscale Processes with Adaptive Ground Observations (RELAMPAGO) and Cloud, Aerosol, and Complex Terrain Interactions (CACTI) projects deployed a high-spatiotemporal-resolution radiosonde network to examine environments supporting deep convection in the complex terrain of central Argentina. This study aims to characterize atmospheric profiles most representative of the near-cloud environment (in time and space) to identify the mesoscale ingredients affecting storm initiation and growth. Spatiotemporal autocorrelation analysis of the soundings reveals that there is considerable environmental heterogeneity, with boundary layer thermodynamic and kinematic fields becoming statistically uncorrelated on scales of 1–2 h and 30 km. Using this as guidance, we examine a variety of environmental parameters derived from soundings collected within close proximity (30 km in space and 30 min in time) of 44 events over 9 days where the atmosphere either: 1) supported the initiation of sustained precipitating convection, 2) yielded weak and short-lived precipitating convection, or 3) produced no precipitating convection in disagreement with numerical forecasts from convection-allowing models (i.e., Null events). There are large statistical differences between the Null event environments and those supporting any convective precipitation. Null event profiles contained larger convective available potential energy, but had low free-tropospheric relative humidity, higher freezing levels, and evidence of limited horizontal convergence near the terrain at low levels that likely suppressed deep convective growth. Here, we also present evidence from the radiosonde and satellite measurements that flow–terrain interactions may yield gravity wave activity that affects CI outcome.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Statistical perspective on embrittling potency for intergranular fracture

Embrittling potency is a thermodynamic metric that assesses the influence of solute segregation to a grain boundary (GB) on intergranular fracture. Historically, authors of studies have reported embrittling potency as a single scalar value, assuming a single segregation site of importance at a GB and a particular cleavage plane. However, the topography of intergranular fracture surfaces is not generally known a priori. Accordingly, we, in this study, present a statistical ensemble approach to compute embrittling potency, where many free surface (FS) permutations are systematically considered to model fracture of a GB. The result is a statistical description of the thermodynamics of GB embrittlement. As a specific example, embrittling potency distributions are presented for Cr segregation to sites at two Ni $\langle 111 \rangle$ symmetric tilt GBs using atomistic simulations. We show that the average embrittling potency for a particular GB site, considering an ensemble of FS permutations, is not equal to the embrittling potency computed using the lowest energy pair of FSs. A mean GB embrittlement is proposed, considering both the likelihood of formation of a particular FS and the probability of solute occupancy at each GB site, to compare the relative embrittling behavior of two distinct GBs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Evaluating Tropical Precipitation Relations in CMIP6 Models with ARM Data

Abstract A set of diagnostics based on simple, statistical relationships between precipitation and the thermodynamic environment in observations is implemented to assess phase 6 of the Coupled Model Intercomparison Project (CMIP6) model behavior with respect to precipitation. Observational data from the Atmospheric Radiation Measurement (ARM) permanent field observational sites are augmented with satellite observations of precipitation and temperature as an observational baseline. A robust relationship across observational datasets between column water vapor (CWV) and precipitation, in which conditionally averaged precipitation exhibits a sharp pickup at some critical CWV value, provides a useful convective onset diagnostic for climate model comparison. While a few models reproduce an appropriate precipitation pickup, most models begin their pickup at too low CWV and the increase in precipitation with increasing CWV is too weak. Convective transition statistics compiled in column relative humidity (CRH) partially compensate for model temperature biases—although imperfectly since the temperature dependence is more complex than that of column saturation. Significant errors remain in individual models and weak pickups are generally not improved. The conditional-average precipitation as a function of CRH can be decomposed into the product of the probability of raining and mean precipitation during raining times (conditional intensity). The pickup behavior is primarily dependent on the probability of raining near the transition and on the conditional intensity at higher CRH. Most models roughly capture the CRH dependence of these two factors. However, compensating biases often occur: model conditional intensity that is too low at a given CRH is compensated in part by excessive probability of precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Synchrotron X-ray diffraction analysis of constituent phases in transition joint between nickel alloy 738LC and a MnFeCoNiCu alloy

Synchrotron x-ray diffraction (XRD) analysis was performed on transition joints between a single-phase MnFeCoNiCu alloy and Ni-base Alloy 738LC to efficiently identify the constituent phases across the interface, with different levels of material intermixing generated by laser-welding with variable power. Here, the intermixing extent was quantified by postmortem energy dispersive x-ray spectroscopy mapping. Diffraction-based analyses on complex composition spaces with limited prior data present significant challenges because elemental substitution in both disordered and ordered phases is more extensive than in conventional alloy spaces, which may effectuate relatively large shifts in the observed lattice parameters that convolute the analysis. Therefore, thermodynamic simulations and crystallographic literature data were employed to construct a system-specific diffraction library of twelve prospective phases for the composition space investigated. Subsequently, for predicted disordered cubic phases, statistical hard-sphere models were established to estimate the lattice parameters and predict diffraction peak positions for inclusion in the library. The library was then employed to analyze diffraction profiles measured from the variably intermixed transition joints, with focus on accounting for both high and low-intensity peaks. 99.0 % of diffraction peaks with relative intensity greater than 0.001 were assigned to phases from the system-specific library, exemplifying rigorous peak accounting and indicating that no unexpected phases were present. Up to six of the twelve library phases were experimentally found in the transition joints. The lattice parameters predicted by the statistical hard-sphere model based upon thermodynamic simulations agree reasonably well with the measured values for the disordered FCC matrix phase.

36 MATERIALS SCIENCE↗