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

Redox Poise during Rhodospirillum rubrum Phototrophic Growth Drives Large-scale Changes in Macromolecular Synthesis Pathways

During photoheterotrophic growth on organic substrates, purple nonsulfur photosynthetic bacteria like Rhodospirillum rubrum can acquire electrons by multiple means, including oxidation of organic substrates, oxidation of inorganic electron donors (e.g., H2), and by reverse electron flow from the photosynthetic electron transport chain. These electrons are stored as reduced electron-carrying cofactors (e.g., NAD(P)H and ferredoxin). The overall ratio of oxidized to reduced cofactors (e.g., NAD(P)+:NAD(P)H), or ’redox poise’, is difficult to understand or predict, as are the cellular processes for dissipating these reducing equivalents. Using physics-based models that capture mass action kinetics consistent with the thermodynamics of reactions and pathways, a range of redox conditions for heterophototrophic growth are evaluated, from conditions in which the NADP+/NADPH levels approach thermodynamic equilibrium to conditions in which the NADP+/NADPH ratio is far above the typical physiological values. Modeling predictions together with experimental measurements indicate that the redox poise of the cell results in large-scale changes in the activity of biosynthetic pathways and, thus, changes in cell macromolecule levels (DNA, RNA, proteins, and fatty acids). Furthermore, modeling predictions indicate that during phototrophic growth, reverse electron flow from the quinone pool is a minor contributor to the production of reduced cofactors (e.g., NAD(P)H) compared to other oxidative processes (H2 and carbon substrate oxidation). Instead, the quinone pool primarily operates to aid ATP production. The high level of ATP, in turn, drives reduction processes even when NADPH levels are relatively low compared to NADP+ by coupling ATP hydrolysis to the reductive processes. The model, in agreement with experimental measurements of macromolecule ratios of cells growing on different carbon substrates, indicates that the dynamics of nucleotide versus lipid and protein production is likely a significant mechanism of balancing oxidation and reduction in the cell.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Technical note: A modified formulation of dynamic energy budget theory for faster computation of biological growth

Abstract. The mass conservation equation in the presence of boundary fluxes and chemical reactions from non-equilibrium thermodynamics is used to derive a modified dynamic energy budget (mDEB) model. Compared to the standard dynamic energy budget (sDEB) model (Kooijman, 2009), this modified formulation does not place the dilution effect in the mobilization kinetics of reserve biomass, and it maintains the partition principle for reserve mobilization dynamics for both linear and non-linear kinetics. Overall, the mDEB model shares most features with the sDEB model. However, for biological growth that requires multiple nutrients, the mDEB model is computationally much more efficient by not requiring numerical iterations for obtaining the specific growth rate. In an example of modeling the growth of Thalassiosira weissflogii in a nitrogen-limiting chemostat, the mDEB model was found to have almost the same accuracy as the sDEB model while requiring almost half of the computing time of the sDEB model. Since the sDEB model has been successfully applied in numerous studies, we believe that the mDEB model can help improve the modeling of biological growth and the associated ecosystem processes in various contexts.

Tang, Jinyun↗
DOE OSTI2024

Spectroscopic Investigations of Molecular Symmetry Breakdown

This project addressed precise measurements of the spectra of molecules to address fundamental questions in chemical physics in the realm of nuclear spin symmetries and their breakdown in molecules, and parity non-conservation caused by the weak nuclear force. Nuclear spin symmetry is intimately connected to parity in symmetrical molecules and plays a central role in the physical and chemical properties of symmetric molecules, their collisional relaxation, and the establishment of thermodynamic equilibrium in an ensemble. Such effects are present in every molecular system or process, but their small size means they are often only directly observable under special conditions. However, their measurement is of fundamental, and often practical importance, and strongly relates to the themes of coherence in light and matter in the recent Basic Energy Sciences Advisory Committee (BESAC) report “Challenges at the Frontiers of Matter and Energy: Transformative Opportunities for Discovery Science”.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Method for electrochemical gas separation

A method for separating a Lewis acid gas from a fluid mixture, comprising contacting the fluid mixture with a reduced electroactive species; a non-aqueous electrolyte; and a stabilizing additive to form an anion adduct between the Lewis acid gas and the reduced electroactive species, wherein the electroactive species comprises an oxidized state, and at least one reduced state that bonds with the Lewis acid gas to form the anion adduct, wherein the stabilizing additive comprises a cationic Lewis acid, a hydrogen-bond donor, or a combination thereof, and the stabilizing additive is present in an effective amount to kinetically favor the forming of the anion adduct from the reduced electroactive species and thermodynamically favor the forming of the anion adduct in the thermodynamic equilibrium between the anion adduct and the reduced electroactive species.

Source record↗
DOE OSTI2026

Determining the Solubility Behavior of Kogarkoite in Simulated Nuclear Waste

Kogarkoite (Na 3 FSO 4 ) is a sparingly soluble fluoride–sulfate double salt that has been identified in high level nuclear waste sludge at the Hanford Site and, more recently, in sludge batch compilation samples at the Savannah River Site (SRS). Due to its complex dissolution behavior, which exhibits an inverse dependence on sodium ion activity, the presence of this mineral poses significant challenges to waste retrieval and processing. Incomplete dissolution during sludge washing can lead to the retention of fluoride and sulfate in the high-level waste feed, potentially causing the formation of corrosive, immiscible molten salt layers, known as "glass gall,” in vitrification melters. Current efforts to optimize flowsheet parameters and wash-water volumes are hindered by the absence of a commercially available, certified reference material, which prevents the accurate calibration of analytical methods and the verification of dissolution kinetics. To address this critical gap, this research focuses on the laboratory synthesis of pure Kogarkoite to serve as a standard for comprehensive solubility and washing performance testing. A coupled synthesis and simulant campaign was executed using an evaporative crystallization protocol designed to replicate the dynamic concentration effects observed in tank farm operations. Thirteen simulant matrices were prepared by dissolving systematically varied ratios of sodium fluoride (NaF) and sodium sulfate (Na 2 SO 4 ) in deionized water under three distinct caustic regimes: 0.0 g (control), 4.0 g (~1 M), and 12.0 g (~3 M) sodium hydroxide (NaOH). While thermodynamic equilibrium models suggest that high-caustic environments should favor the stability of the double salt7, results from this evaporative study at 25 0 C revealed a distinct kinetic divergence. Simulants with high hydroxide loading predominantly yielded large, blocky crystals of sodium sulfate decahydrate (Na 2 SO 4 .10H 2 O). Successful synthesis of pure Kogarkoite was achieved exclusively in specific NaOH-free compositional windows, where the precipitate manifested as fine, opaque granular aggregates. Ion chromatography (IC) analysis confirmed phase purity through the simultaneous stoichiometric depletion of both fluoride and sulfate from the supernatant. This successful synthesis establishes a reproducible route to generate bulk Kogarkoite, enabling the subsequent phase of quantitative dissolution testing using inhibited water to optimize sludge-batch assembly.

Sarker, Md Sharif [Florida International Univ. (FI↗
DOE OSTI2024

Quantitative Predictive Theories through Integrating Quantum, Statistical, Equilibrium, and Nonequilibrium Thermodynamics

Today's thermodynamics is largely based on the combined law for equilibrium systems and statistical mechanics derived by Gibbs in 1873 and 1901, respectively, while irreversible thermodynamics for nonequilibrium systems resides essentially on the Onsager Theorem as a separate branch of thermodynamics developed in 1930s. Between them, quantum mechanics was invented and was quantitatively solved in terms of density functional theory (DFT) in 1960s. Furthermore, these three scientific domains operate based on different principles and are very much separated from each other. In analogy to the parable of the blind men and the elephant articulated by Perdew, they individually represent different portions of a complex system and thus are incomplete by themselves alone, resulting in the lack of quantitative agreement between their predictions and experimental observations. Over the last two decades, the author's group has developed a multiscale entropy approach (recently termed as zentropy theory) that integrates DFT-based quantum mechanics and Gibbs statistical mechanics and is capable of accurately predicting entropy and free energy of complex systems. Furthermore, in combination with the combined law for nonequilibrium systems developed by Hillert, the author developed the theory of cross phenomena beyond the phenomenological Onsager Theorem. The zentropy theory and theory of cross phenomena jointly provide quantitative predictive theories for systems from electronic to any observable scales as reviewed in the present work.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

On Gibbs Equilibrium and Hillert Nonequilibrium Thermodynamics

During his time at Royal Institute of Technology (Kungliga Tekniska högskolan) in Sweden, the present author learned nonequilibrium thermodynamics from Mats Hillert. The key concepts are the separation of internal and external variables of a system and the definitions of potentials and molar quantities. In equilibrium thermodynamics derived by Gibbs, the internal variables are not independent and can be fully evaluated from given external variables. While irreversible thermodynamics led by Onsager focuses on internal variables though often mixed with external variables. Hillert integrated them together by first emphasizing their differences and then examining their connections. His philosophy was reflected by the title of his book “Phase Equilibria, Phase Diagrams and Phase Transformations” that puts equilibrium, nonequilibrium, and internal processes on equal footing. Here, in the present paper honoring Hillert, the present author reflects his experiences with Hillert and his work in last 40 years and expresses his gratitude for all the wisdom and support from him in terms of “Hillert nonequilibrium thermodynamics” and discusses some recent topics that the present author has been working on.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Corrosion of 316 SS in direct contact with Zircaloy-4 in repository relevant conditions

The effects of interfacial contact between 316 stainless steel (316 SS) waste canister and spent fuel cladding material (Zircaloy-4) in case of water intrusion into a deep geological repository are investigated. Hydrothermal corrosion experiments of 316 SS and Zircaloy-4 in direct contact were performed at elevated temperatures in MilliQ water and simulated Mont Terri groundwater. Hematite, magnetite, chromite, and baddeleyite form as secondary phases. pH of the Mont Terri solution increases with corrosion duration and magnetite forms preferentially instead of hematite. Thermodynamic modeling shows that iron oxidation in 316 SS controls the pH and redox potential of the aqueous environment and is in good agreement with experimental Mont Terri results, showing that experimental Mont Terri results trend to equilibrium. All oxides except hematite are predicted to exist at thermodynamic equilibrium based on calculations, implying that experimentally formed hematite will be reduced or transformed to magnetite at equilibrium. The thermodynamically predicted chromite was observed to form only at the interface where 316 SS was in contact with Zircaloy-4 at 250°C, showing that contact with Zircaloy-4 seems to accelerate the 316 SS corrosion to thermodynamic equilibrium. Finally, Mont Terri groundwater readily trends towards thermodynamic equilibrium, leading to a pH and redox buffered system that is dominated by iron in solution coming from the accelerated corrosion of 316 SS caused by Zircaloy-4 contact.

316 stainless steel↗
DOE OSTI2022

Temperature screening and cross-field impurity accumulation from a thermodynamic perspective

In a variety of different systems, high-Z ion species show a marked tendency to accumulate in regions of high plasma density. It has previously been suggested that the apparent universality of this behavior could be explained thermodynamically, in terms of the maximum-entropy state attainable when the system must obey an ambipolarity condition. Furthermore, the previous treatment did not allow for the possibility of temperature gradients. Here, tools from non-equilibrium (Onsager) thermodynamics are used to show that ambipolarity continues to play a key role in producing this behavior in the presence of temperature gradients, and to recover well-known temperature screening effects that appear in these cases.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Quantum Foundations of Classical Reversible Computing

The reversible computation paradigm aims to provide a new foundation for general classical digital computing that is capable of circumventing the thermodynamic limits to the energy efficiency of the conventional, non-reversible digital paradigm. However, to date, the essential rationale for, and analysis of, classical reversible computing (RC) has not yet been expressed in terms that leverage the modern formal methods of non-equilibrium quantum thermodynamics (NEQT). In this paper, we begin developing an NEQT-based foundation for the physics of reversible computing. We use the framework of Gorini-Kossakowski-Sudarshan-Lindblad dynamics (a.k.a. Lindbladians) with multiple asymptotic states, incorporating recent results from resource theory, full counting statistics and stochastic thermodynamics. Important conclusions include that, as expected: (1) Landauer’s Principle indeed sets a strict lower bound on entropy generation in traditional non-reversible architectures for deterministic computing machines when we account for the loss of correlations; and (2) implementations of the alternative reversible computation paradigm can potentially avoid such losses, and thereby circumvent the Landauer limit, potentially allowing the efficiency of future digital computing technologies to continue improving indefinitely. We also outline a research plan for identifying the fundamental minimum energy dissipation of reversible computing machines as a function of speed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Adsorption of terbium (III) on DGA and LN resins: Thermodynamics, isotherms, and kinetics

Two commercially available extraction chromatography (EXC) resins containing N,N,N’,N’-tetra-n-octyldiglycolamide (DGA Resin, Normal, 50 – 100 μm) and Bis(2-ethylhexyl) phosphate (LN Resin, 100 – 150 μm) were used as adsorbents to study fundamental adsorption properties such as thermodynamic values, equilibrium isotherms, and kinetic uptake models for terbium(III) adsorption. Weight distribution ratios (D w ) for terbium on DGA and LN resins were measured using a [ 160 Tb]Tb 3+ radiometric tracer in nitric acid as a function of acidity, temperature, initial analyte concentration, and equilibrium time. The D w values showed increasing binding affinity for DGA resin at high nitric acid concentrations and decreasing binding affinity for LN resins. Thermodynamic studies for DGA and LN resins revealed that the Gibbs free energy (ΔG) increased consistently with temperature. To model equilibrium data, increasingly higher parameter equilibrium isotherm models (Henry (1) < Langmuir, Freundlich (2) < Redlich-Peterson (3) < Fritz-Schluender (4)) were compared on their root mean squared errors (RMSE) and adjusted determination coefficients to determine the most applicable model. In all cases, the empirical four-parameter Fritz-Schluender isotherm demonstrated a superior fit. Similar comparisons for reaction-based kinetic models (Pseudo-first-order < Pseudo-second-order < Pseudo-n-order) revealed that the higher-order PNO model yielded a superior fit of kinetic data for both resins. Furthermore, in some cases, adsorption isotherms and kinetic models could also be modeled by a lower-order model with minimal change in error parameters. Weber-Morris plots revealed that two linear sections are observed for each resin, where the first linear segment is attributed to fast (film diffusion) adsorption of terbium, followed by slower intraparticle diffusion of terbium through the pores as the rate-limiting step. Based on the Weber-Morris plot, both film and intraparticle diffusion are involved in controlling the kinetic rate of adsorption for DGA and LN resins.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Nonadiabatic Force Matching for Alchemical Free-Energy Estimation

We propose a method to compute free-energy differences from nonadiabatic alchemical transformations by using flow-based generative models. The method, nonadiabatic force matching, hinges on estimating the dissipation along an alchemical switching process in terms of a nonadiabatic force field that can be learned through stochastic flow matching. The learned field can be used in conjunction with short-time trajectory data to evaluate upper and lower bounds on the alchemical free energy that variationally converge to the exact value if the field is optimal. Applying the method to evaluate the alchemical free energy of atomistic models shows that it can substantially reduce the simulation cost of a free-energy estimate at a negligible loss of accuracy when compared with thermodynamic integration.

Computational chemistry↗
DOE OSTI2026

Chemically Recyclable Analogs of Styrene–Butadiene Copolymers Enabling Perfectly Linear Ethylene–Styrene Materials with Random Phenyl Distribution

Copolymerization of cyclopentene (CP) and 4-phenylcyclopentene (4PCP) at a full range of comonomer feed ratios is reported using Ru-based ring-opening metathesis polymerization (ROMP) yielding homogeneous copolymers analogous to poly(styrene-ran-1,4-butadiene) and poly(ethylene-ran-styrene) copolymers following hydrogenation under mild conditions. In all cases, total monomer conversions of 86%–92% yielded copolymers with compositions within 4% of monomer feeds. Analysis of equilibrium copolymerization thermodynamics, rarely performed on two cycloolefin monomers with low ring strain energies, provides rational design strategies for negotiating two monomers with different equilibrium monomer concentrations. Inverse-gated decoupled 13 C NMR analysis of dyad sequences on the resulting copolymer microstructures concludes a near-random distribution of comonomer units. The copolymers produced from ROMP have number-average molar masses up to 60 kg mol –1 , moderate dispersities (1.5 ± 0.1), and high trans olefin content (86% ± 2%) while glass-transitions temperatures follow the Fox equation and span the full range between homopolymer extremities of PCP (−96 °C) and P4PCP (17 °C). Unlike most prevulcanized elastomers, these materials undergo facile chemical recycling to monomer, producing complete ring-closing metathesis depolymerization (RCMD) of the polymer back to the CP comonomers. Quantitative olefin hydrogenation produced perfectly linear polyethylene with 4%–16% of the backbone carbons containing a phenyl pendant, analogous to ES copolymers with up to 71.5% w/w styrene units but with random distribution of the aromatic pendants. Thermal properties of these materials are discussed, which span from semicrystalline to amorphous, and with T g values notably less than the reported ES copolymer analogs at similar compositions.

animal feed↗
DOE OSTI2020

Thermodynamics and kinetics of solution combustion synthesis: Ni(NO 3 ) 2 + fuels systems

Solution combustion synthesis (SCS) utilizes exothermic self-propagating reactions to prepare nanoscale materials that can be used widely in energy, electronics, and biomedical technologies and other applications. SCS is a specific variety of a more general combustion synthesis (CS) method. Investigations of the thermodynamics, kinetics, and the mechanisms of SCS reactions, are not as well studied as the other CS processes. This work reports on a systematic study of the thermodynamics and kinetics of SCS reactions involving Ni(NO 3 ) 2 , an oxidizer, and either glycine (C 2 H 5 NO 2 ) or hexamethylenetetramine (HMT, C 6 H 12 N 4 ) as fuels. A thermodynamic modeling approach, based on the Gibbs free energy minimization principle, is applied to the simultaneous calculations of the adiabatic temperatures and compositions of the equilibrium products. Our calculations reveal the influence of fuel-to-oxidizer ratio, amount of water, and the oxygen in air on the combustion temperature under adiabatic conditions and the composition of the resulting products. We have, in turn, measured the combustion temperature and phase composition of products and compared them with the calculations. Variations of drying times for the solutions yield precursor gels with varying water contents. This approach enables the manipulation of combustion parameters and confirms the use of calculated activation energies for reactions using the Merzhanov-Khaikin method. Here, the results show that SCS reactions in fuel-lean solutions producing NiO have higher activation energy in contrast to reactions with fuel-rich solutions that form Ni. Reduction of activation energies due to the increase in the fuel-to-oxidizer ratio could be related to the observed change of the rate-limiting stages of the endothermic decomposition of the individual reactants to the exothermic decomposition of coordinate compounds formed between the reactants.

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

Predicting Morphological Evolution during Coprecipitation of MnCO 3 Battery Cathode Precursors Using Multiscale Simulations Aided by Targeted Synthesis

The performance of lithium-ion batteries is intimately linked to both the structure and the morphology of the cathode material, which in turn is critically linked to the synthesis conditions. However, few studies focus on understanding synthesis, especially during the coprecipitation of metal oxide precursors, a process that largely determines the final morphology of the material. In this paper, we go beyond the typical equilibrium particle shape analysis conducted in the literature and incorporate kinetic aspects of morphology evolution. We perform these studies using controlled synthesis on a well-defined metal salt system (MnCO 3 ) combined with multiscale simulations and high-resolution microscopy. Results show that with increasing metal concentration, the particles transition from rhombohedral to cubic to spherical shapes. Computational analysis using density functional theory (DFT) reveals that rhombohedral shaped particles evolve under equilibrium conditions. Phase field techniques indicate that at higher metal concentrations, fast growth kinetics of the precipitates result in the transition to cubic and, subsequently, spherical shapes, accompanied by a decrease in particle size. This study, while limited to the one metal salt system, provides an approach to shed light on the synthesis process of mixed transition metal salts, gradient materials, and other cathode materials of interest to the battery community.

25 ENERGY STORAGE↗
DOE OSTI2022

Q-POP-Thermo: A general-purpose thermodynamics solver for ferroelectric materials

We report that Q-POP-Thermo is a program designed to compute thermodynamic monodomain equilibrium states and their properties for ferroelectric single crystals and thin films based on the Landau-Ginzburg-Devonshire (LGD) Theory. Utilizing symbolic manipulation with the SymPy Library, the governing equations along with appropriate boundary conditions are solved for speedy minimization of the free energy of a crystal. Utilizing the popular Differential Evolution algorithm, with appropriate hybridization, multiple phase diagrams, such as the pressure-temperature phase diagram for bulk single crystals and the common strain-temperature phase diagram for monodomain thin-film systems can be readily generated. Furthermore, a variety of material properties of stable ferroelectric phases, including dielectric, piezoelectric, and electrocaloric properties, can simultaneously be calculated. Validation studies are presented for both thin-film and single crystal systems to test the effectiveness and capability of the open-source program.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Thermal exchange-correlation functionals: Capturing quantum electron behavior in warm, dense plasmas

We summarize and give perspective upon recent progress in developing non-empirical constraint-based thermal (i.e., free energy) exchange-correlation (XC) density functionals essential for accurate description of the quantum behavior of electrons in warm, dense plasmas. After delineating the critical role of ground-state functionals for zero-temperature, time-dependent DFT, we outline the underpinnings of local density approximation, generalized gradient approximation (GGA), and meta-GGA XC free-energy functionals. Two basic thermalization principles for upgrading ground-state XC functionals to successful thermal ones are emphasized. Then, we turn to a long-standing challenge, assessment of the accuracy of well-founded functionals. Unlike the ground state, there are a few exact results for large T and P. An exception is path integral Monte Carlo (PIMC) data for dense H/D and He plasmas. For those, we did ab initio molecular dynamics simulations under selected thermodynamic conditions employing five thermal XC functionals: two approximate thermal GGAs, fully thermal GGA, an approximate meta-GGA, and fully thermal meta-GGA. Comparisons with the PIMC data show that functionals thermalized by augmenting a non-thermal functional with a lower-level thermal contribution are inferior to functionals with thermal XC and spatial inhomogeneity effects taken into account at the same level of refinement. We believe this and similar evidence should be convincing to the high-energy density physics community of the necessity of use of proper thermal XC functionals in simulation studies of finite-temperature quantum effects in warm, dense plasmas.

Ab-initio molecular dynamics↗