Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermodynamic equilibrium”

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

Uncertainty quantification for equations of state: copper as an example

Equations of state are essential for providing a fundamental description of materials properties in thermodynamic equilibrium and are used to provide closure relations for hydrodynamics simulations. Generally, equations of state rely on simple physics-based parameterized materials models to inform on the free energy of a material through out a given thermodynamic state space. Historically the parameters of these models have been tuned by hand to fit various experimental data. However, modern optimization and uncertainty quantification techniques allow us to quickly test thousands of parameter combinations and obtain meaningful uncertainty estimates on the parameters, opening opportunities for assessing systematic uncertainties in experiments, assessing model adequacy, and more. In this report, we use Bayesian inference to fit the solid (fcc) equation of state of copper. We focus on fitting five different experimental datasets, including the isobaric density, isobaric heat capacity, room temperature isotherm, principal isentrope, and principal Hugoniot. We fit all five data types simultaneously, and then explore the extent to which combinations of 2 subsets of the 5 datasets can constrain the EOS parameters, as compared to the fit to all 5. This information is useful for investigating the extent to which different datasets can con strain EOS models and thereby help guide experimental investigations in order to best constrain the EOS. We also discuss ways that the methodologies can be used to investigate systematic discrepancies between experiments, as well as how the methods can be used to assess model uncertainty. The framework we develop is general, in that it can be used with a variety of optimization or uncertainty quantification techniques and with a variety of data sources, including both experimental and ab-inito data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Molecular beam epitaxy growth of IrO 2 using plasma-only oxidation

Growth of IrO 2 films via molecular beam epitaxy (MBE) using an oxygen plasma as the sole oxidant is demonstrated for the first time. We investigate the oxidizing conditions required for IrO 2 by characterizing the species present in the plasma using optical emission spectroscopy and comparing the thermodynamic equilibrium between Ir/IrO 2 under atomic and molecular oxygen. Across the pressure range accessible, monatomic oxygen was the primary reactive species present in the plasma. IrO 2 films were grown under varying Ir flux, oxygen pressure, and substrate temperature to understand growth thermodynamics and kinetics. Structural characterization assessed film phase and composition, crystallinity, strain, and surface morphology. The optimized films from this study validate plasma-only reactive growth by MBE as a successful method for growing IrO 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2021

Elastic forces drive nonequilibrium pattern formation in a model of nanocrystal ion exchange

The widely used process of nanocrystal ion exchange operates out of thermodynamic equilibrium and can require mixing components of varying sizes. Here we use theory and computer simulation to study a simple model which captures these two basic features of ion exchange reactions. We show that a strong driving force for exchange among different-sized species creates nonequilibrium patterns within model nanocrystals. We further demonstrate that such patterns can be thermodynamically stable within core/shell nanocrystals. These results help us understand the heterostructures formed in ion-exchanged nanocrystals and suggest strategies for leveraging elasticity to design patterns in nanoscale materials. Chemical transformations, such as ion exchange, are commonly employed to modify nanocrystal compositions. Yet the mechanisms of these transformations, which often operate far from equilibrium and entail mixing diverse chemical species, remain poorly understood. Here we explore an idealized model for ion exchange in which a chemical potential drives compositional defects to accumulate at a crystal’s surface. These impurities subsequently diffuse inward. We find that the nature of interactions between sites in a compositionally impure crystal strongly impacts exchange trajectories. In particular, elastic deformations which accompany lattice-mismatched species promote spatially modulated patterns in the composition. These same patterns can be produced at equilibrium in core/shell nanocrystals, whose structure mimics transient motifs observed in nonequilibrium trajectories. Moreover, the core of such nanocrystals undergoes a phase transition—from modulated to unstructured—as the thickness or stiffness of the shell is decreased. Our results help explain the varied patterns observed in heterostructured nanocrystals produced by ion exchange and suggest principles for the rational design of compositionally patterned nanomaterials.

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

Exploring the Structural, Dynamic, and Functional Properties of Metal–Organic Frameworks through Molecular Modeling

This review spotlights the role of atomic-level modeling in research on metal-organic frameworks (MOFs), especially the key methodologies of density functional theory (DFT), Monte Carlo (MC) simulations, and molecular dynamics (MD) simulations. The discussion focuses on how periodic and cluster-based DFT calculations can provide novel insights into MOF properties, with a focus on predicting structural transformations, understanding thermodynamic properties and catalysis, and providing information or properties that are fed into classical simulations such as force field parameters or partial charges. Classical simulation methods, highlighting force field selection, databases of MOFs for high-throughput screening, and the synergistic nature of MC and MD simulations, are described. By predicting equilibrium thermodynamic and dynamic properties, these methods offer a wide perspective on MOF behavior and mechanisms. Additionally, the incorporation of machine learning (ML) techniques into quantum and classical simulations is discussed. These methods can enhance accuracy, expedite simulation setup, reduce computational costs, as well as predict key parameters, optimize geometries, and estimate MOF stability. By charting the growth and promise of computational research in the MOF field, the aim is to provide insights and recommendations to facilitate the incorporation of computational modeling more broadly into MOF research.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Digital quantum simulation of Floquet symmetry-protected topological phases

Quantum many-body systems away from equilibrium host a rich variety of exotic phenomena that are forbidden by equilibrium thermodynamics. A prominent example is that of discrete time crystals, in which time-translational symmetry is spontaneously broken in periodically driven systems. Pioneering experiments have observed signatures of time crystalline phases with trapped ions, solid-state spin systems, ultracold atoms and superconducting qubits. Here we report the observation of a distinct type of non-equilibrium state of matter, Floquet symmetry-protected topological phases, which are implemented through digital quantum simulation with an array of programmable superconducting qubits. We observe robust long-lived temporal correlations and subharmonic temporal response for the edge spins over up to 40 driving cycles using a circuit of depth exceeding 240 and acting on 26 qubits. We demonstrate that the subharmonic response is independent of the initial state, and experimentally map out a phase boundary between the Floquet symmetry-protected topological and thermal phases. Our results establish a versatile digital simulation approach to exploring exotic non-equilibrium phases of matter with current noisy intermediate-scale quantum processors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Magnetism in the two-dimensional dipolar XY model

Here, motivated by a recent experiment on a square-lattice Rydberg atom array realizing a long-range dipolar XY model [C. Chen et al., Nature (London) 616, 691 (2023)], we numerically study the model's equilibrium properties. We obtain the phase diagram, critical properties, entropies, variance of the magnetization, and site-resolved correlation functions. We consider both ferromagnetic and antiferromagnetic interactions and apply quantum Monte Carlo and pseudo-Majorana functional renormalization group techniques, generalizing the latter to a U⁡(1) symmetric setting. Our simulations perform extensive thermometry in dipolar Rydberg atom arrays and establish conditions for adiabaticity and thermodynamic equilibrium. On the ferromagnetic side of the experiment, we determine the entropy per particle S/N≈0.5, close to the one at the critical temperature, S c /N=0.585⁢(15). The simulations suggest the presence of an out-of-equilibrium plateau at large distances in the correlation function, thus motivating future studies on the nonequilibrium dynamics of the system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Quantifying and visualizing the microscopic degrees of freedom of grain boundaries in the Wigner–Seitz cell of the displacement-shift-complete lattice

Here we introduce a grain boundary (GB) translation vector, t ws , to describe and quantify the domain of the microscopic degrees of freedom of GBs. It has long been recognized that for fixed macroscopic degrees of freedom of a GB there exists a large multiplicity of states characterized by different relative grain translations. More recently another degree of freedom, [n], the number of GB atoms, has emerged and is now recognized as an equally important component of GB structural multiplicity. In this work, we show that all GB microstates can be uniquely characterized by their value of t ws , which is located within the Wigner–Seitz (WS) cell of the Displacement-Shift-Complete lattice (DSCL) of the GB. The GB translation vector captures information about both the translation state and the number of GB atoms. We show that the density of GB microstates inside the WS cell of the DSCL is not uniform and can form clusters that correspond to different GB phases. Further, the vectors connecting the centers of the clusters correspond to the Burgers vectors of GB phase junctions, which can be predicted without building the junctions. Using t ws , we quantify GB excess shear and argue that it is defined up to a DSCL vector, which has implications for thermodynamic equilibrium conditions. Additionally, this work generalizes the definition of the number of GB atoms [n] to asymmetric boundaries.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Comparison of Equilibrium and Nonequilibrium Approaches for Relative Binding Free Energy Predictions

Alchemical relative binding free energy calculations have recently found important applications in drug optimization. A series of congeneric compounds are generated from a preidentified lead compound, and their relative binding affinities to a protein are assessed in order to optimize candidate drugs. While methods based on equilibrium thermodynamics have been extensively studied, an approach based on nonequilibrium methods has recently been reported together with claims of its superiority. However, these claims pay insufficient attention to the basis and reliability of both methods. Here we report a comparative study of the two approaches across a large data set, comprising more than 500 ligand transformations spanning in excess of 300 ligands binding to a set of 14 diverse protein targets. Ensemble methods are essential to quantify the uncertainty in these calculations, not only for the reasons already established in the equilibrium approach but also to ensure that the nonequilibrium calculations reside within their domain of validity. If and only if ensemble methods are applied, we find that the nonequilibrium method can achieve accuracy and precision comparable to those of the equilibrium approach. Compared to the equilibrium method, the nonequilibrium approach can reduce computational costs but introduces higher computational complexity and longer wall clock times. There are, however, cases where the standard length of a nonequilibrium transition is not sufficient, necessitating a complete rerun of the entire set of transitions. This significantly increases the computational cost and proves to be highly inconvenient during large-scale applications. Our findings provide a key set of recommendations that should be adopted for the reliable implementation of nonequilibrium approaches to relative binding free energy calculations in ligand-protein systems.

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

Time-dependent water vapor desorption isotherm model of hardened cement paste

Highlights: • Short-term and long-term desorption isotherms of hardened cement paste are reviewed. • Water content at RH depends on both thermodynamic and microstructural equilibria. • An extended GAB equation is proposed to fit the cavitation drop during desorption. • A dynamic time-dependent sorption isotherm model is shown and used for simulations. • Numerical results capture anomalous drying in the same RH as seen experimentally. This study aims to enhance moisture transport modeling by elucidating the so-called anomalous water diffusion in cementitious materials. Water desorption isotherms are measured from samples at various drying stages using techniques with different durations to quantify the impact of the drying duration on the water sorption ability. A single water sorption isotherm does not solely give a relationship between the water content and relative humidity but also represents the state of microstructure. The continuous evolution of a desorption isotherm due to drying-induced microstructural rearrangement is demonstrated. For numerical modeling, the microstructural alteration can be explicitly considered through a dynamic desorption isotherm model, which governs the local thermodynamic equilibrium at the capillary meniscus. This approach is implemented into a multiphase transport code, whose ability to predict drying is validated using literature data. Finally, the effect of prolonged drying on the colloidal nature of the calcium-silicate-hydrate gel is discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Soft X-ray and EUV emission spectra of beryllium plasma produced by neodymium-glass laser radiation with broad frequency and angular spectra

We present the results of an experimental study of soft X-ray (SXR) and extreme ultraviolet emission spectra of the plasma produced by exposing a plane solid beryllium target to laser radiation with broad frequency and angular spectra. SXR lines up to 1s – 9p of Be IV as well as the plasma continuum are recorded for a laser focal-spot intensity of 5.3 × 10{sup 13} W cm{sup −2}. To model the SXR beryllium plasma spectra, simulations are carried out using the INDHAUS programme and the FLYCHK code in the framework of local thermodynamic equilibrium model, which agree nicely with experimentally obtained results. (laser plasma)

36 MATERIALS SCIENCE↗
DOE OSTI2020

Impact of thermodynamic fluctuations and pattern size on the nucleation behavior during area selective deposition

In this work, the authors explore the impact that thermodynamic fluctuations have on the spontaneous appearance of defects during area selective deposition on patterned surfaces that are fully covered by adsorbates under local thermodynamic equilibrium, such as self-assembled monolayers. By using a simple lattice gas model for the adsorbed monolayer, the authors were able to track the spontaneous formation of defects in the monolayer as a function of the pattern width. The results indicate that, for pattern widths of the order of tens of nanometers, roughening effects at the pattern edge can be the leading source for the spontaneous appearance of nucleation defects. This leads to an enhancement of the density of defects that can be up to three orders of magnitude higher than those expected in uniform (not patterned) surfaces. The model also predicts a density of defects that is inversely proportional to the pattern width. Finally, if the dynamic nucleation of defects during area selective deposition is driven by thermodynamic fluctuations, the model predicts that the nucleation rate should be proportional to the total precursor fluency and independent of purge times. Moreover, a tight confinement of the monolayer through a high quality smooth interface in the patterned substrate and strong cohesive interactions between adsorbates should each contribute to a reduction of the overall defect density.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Collisional Excitation of HCN by CO to Refine the Modeling of Cometary Comae

Here, we present the first dataset of collisional (de)-excitation rate coefficients of HCN induced by CO, one of the main perturbing gases in cometary atmospheres. The dataset spans the temperature range of 5–50 K. It includes both state-to-state rate coefficients involving the lowest ten and nine rotational levels of HCN and CO, respectively, and the so-called “thermalized” rate coefficients over the rotational population of CO at each kinetic temperature. The derivation of these coefficients exploited the good performance of the statistical adiabatic channel model (SACM) on top of an accurate interaction potential computed at the CCSD(T)-F12b/CBS level of theory. The reliability of the SACM approach was validated by comparison with full quantum calculations restricted at the lowest total angular momentum of the system. These results provide essential input to accurately model the distribution among the rotational energy levels and the abundance of HCN in cometary atmospheres, accounting for deviations from local thermodynamic equilibrium that typically occurs in such environments.

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

Weak scaling of the contact distance between two fluctuating interfaces with system size

A pair of flat parallel surfaces, each freely diffusing along the direction of their separation, will eventually come into contact. If the shapes of these surfaces also fluctuate, then contact will occur when their centers-of-mass remain separated by a nonzero distance ℓ. An example of such a situation is the motion of interfaces between two phases at conditions of thermodynamic coexistence, and in particular the annihilation of domain wall pairs under periodic boundary conditions. Here we present a general approach to calculate the probability distribution of the contact distance ℓ and determine how its most likely value ℓ* depends on the surfaces' lateral size L. Using the Edward-Wilkinson equation as a model for interfaces, we demonstrate that ℓ* scales weakly with system size, i.e., the dependence of ℓ* on L for both (1+1)- and (2+1)-dimensional interfaces is such that lim L→∞ (ℓ*/L) = 0. In particular, for (2+1)-dimensional interfaces ℓ* is an algebraic function of logL, a result that is confirmed by computer simulations of slab-shaped domains formed under periodic boundary conditions. Overall, this weak scaling implies that such domains remain topologically intact until ℓ becomes very small compared to the lateral size of the interface, contradicting expectations from equilibrium thermodynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Equilibrium distribution of point defects in Fe-Y-O as a typical representative of nanocluster-strengthened alloys

The statistical mechanic approach has been used to estimate the equilibrium defect distribution in a multiphase system by minimizing its free energy using microstructural characteristics and a finite set of defects obtained by first-principles calculations. Recently the approach was extended to take in consideration defects at the interface of the precipitate and matrix. Herein we apply the developed approach to the investigation of the stability of bcc Fe containing yttria, Y 2 O 3 , nanoclusters as a prototype of nanostructured ferritic alloys. It has been obtained that (100)Fe-O interface is unstable with respect to vacancies production. In a contrast to (100)Fe-O interface, the so called Klim interface is stable, i.e. local vacancy concentration at this interface at 600 K is below 10 -12 . It has been demonstrated that due to large defect formation energies the ODS particles are extremely stable and the main defect corresponds to Fe atoms substitute Y in Y 2 O 3 precipitate. Moreover, under thermodynamic equilibrium condition, the preexist vacancies in bulk Fe do not accumulate oxygen atom. The later observation not necessarily forbid the existence of a large amount of preexisting Fe vacancy - oxygen atom clusters at the initial stages of alloy formation far from equilibrium.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ionic Peltier effect in Li-ion electrolytes

The coupled transport of charge and heat provide fundamental insights into the microscopic thermodynamics and kinetics of materials. We describe a sensitive ac differential resistance bridge that enables measurements of the temperature difference on two sides of a coin cell with a resolution of better than 10 μK. We use this temperature difference metrology to determine the ionic Peltier coefficients of symmetric Li-ion electrochemical cells as a function of Li salt concentration, solvent composition, electrode material, and temperature. The Peltier coefficients Π are negative, i.e., heat flows in the direction opposite to the drift of Li ions in the applied electric field, large, –Π > 30 kJ mol –1 , and increase with increasing temperature at T > 300 K. The Peltier coefficient is approximately constant on time scales that span the characteristic time for mass diffusion across the thickness of the electrolyte, suggesting that heat of transport plays a minor role in comparison to the changes in partial molar entropy of Li at the interface between the electrode and electrolyte. Here, our work demonstrates a new platform for studying the non-equilibrium thermodynamics of electrochemical cells and provides a window into the transport properties of electrochemical materials through measurements of temperature differences and heat currents that complement traditional measurements of voltages and charge currents.

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

Model comparisons for two-temperature plasma equations of state

When a plasma is generated in the laboratory, energy is often deposited preferentially into either the electrons or the ions, giving rise to a quasiequilibrium state in which the two species, electrons and ions, are well described by two effective temperatures, T e and T i . Accurate hydrodynamic modeling of such a two-temperature plasma requires an equation of state that captures the relevant many-body physics without assuming a strict local thermodynamic equilibrium. Several models have been proposed within the literature, which extend conventional statistical approaches, each employing a different combination of assumptions for modifying the equilibrium equations. Here, in this work, we compare the predictions for several models, presenting derivations of the internal energy and pressure for each microscopic model within a unified framework so that the assumptions of each model may be more easily compared to one another. We find that for sufficiently weak coupling, all models agree with one another. However, as the coupling strength is increased, the disagreement between the models becomes more pronounced. Moreover, the relative sizes of the corrections predicted by each model depend on which species has the higher temperature, T e >T i vs T e < T i .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Quantifying Energy Conversion in Higher-Order Phase Space Density Moments in Plasmas

Weakly collisional and collisionless plasmas are typically far from local thermodynamic equilibrium (LTE), and understanding energy conversion in such systems is a forefront research problem. The standard approach is to investigate changes in internal (thermal) energy and density, but this omits energy conversion that changes any higher-order moments of the phase space density. In this Letter, we calculate from first principles the energy conversion associated with all higher moments of the phase space density for systems not in LTE. Particle-in-cell simulations of collisionless magnetic reconnection reveal that energy conversion associated with higher-order moments can be locally significant. Furthermore, the results may be useful in numerous plasma settings, such as reconnection, turbulence, shocks, and wave-particle interactions in heliospheric, planetary, and astrophysical plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Boron nitride nanotube precursor formation during high-temperature synthesis: kinetic and thermodynamic modelling

In this work, we performed integrated modelling of the chemical pathways of formation for boron nitride nanotube (BNNT) precursors during high-temperature synthesis in a B/N 2 mixture. Integrated modelling includes quantum chemistry, Quantum–classical molecular dynamics, thermodynamic modelling, and kinetic approaches. We demonstrate that BN compounds are formed via the interaction of molecular nitrogen with small boron clusters, rather than through interactions with less reactive liquid boron. (This process can also be described as N 2 molecule fixation.) Liquid boron evaporates to produce these boron clusters (B m with m ≤ 5), which are subsequently converted into B m N n chains. The production of such chains is crucial to the growth of BNNTs because these chains form the building blocks of bigger and longer BN chains and rings, which are in turn the building blocks of fullborenes and BNNTs. Additionally, kinetic modelling revealed that B 4 N 4 and B 5 N 4 species in particular play a major role in the N 2 molecule fixation process. Furthermore, the formation of these species via reactions with B 4 and B 5 clusters is not adequately described under the assumption of thermodynamic equilibrium, as is demonstrated in our kinetic modelling. Thus, the accumulation of both B 4 N 4 and B 5 N 4 depends on the background gas pressure and the gas cooling rate. Long BN chains and rings, which are precursors of the fullborene and BNNT growth, form via self-assembly of components B 4 N 4 and B 5 N 4 . Our modelling results—particularly the increased densities of B 4 N 4 and B 5 N 4 species at higher gas pressures—explain the experimentally observed effect of gas pressure on the yield of high-quality BNNTs. The catalytic role of hydrogen was also studied; it is shown that HBNH molecules can be the main precursor of BNNT synthesis in the presence of hydrogen.

77 NANOSCIENCE AND NANOTECHNOLOGY↗