Search NASA⌕ Search

SEARCH · Search NASA

Results for “Thermodynamic systems”

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

High-Precision Calculation of Nanoparticle (Nanocrystal) Potentials of Mean Force and Internal Energies

Thermodynamic stability assessment of nanocrystal systems requires precise free energy calculations. This study highlights the importance of meticulous control over various factors, including the thermostat, time step, potential cutoff, initial configuration, sampling method, and overall simulation duration. Free energy computations in dry (solvent-free) systems are on the order of several hundred $k_BT$ but can be obtained with consistent accuracy. However, calculation of internal energies becomes challenging, as they are typically much larger in magnitude than free energies and exhibit significant noise and reduced reliability. To address this limitation, we propose a new internal energy estimate that drastically reduces the noise. Here, we also present formulas that enable the optimization of the parameters of the harmonic bias potential for optimal convergence. Finally, we discuss the implications of these findings for the computation of free energies in nanocrystal clusters and superlattices.

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

Resolving the Coverage Dependence of Surface Reaction Kinetics with Machine Learning and Automated Quantum Chemistry Workflows

Microkinetic models for catalytic systems require estimation of many thermodynamic and kinetic parameters that can be calculated for isolated species and transition states using ab initio methods. However, the presence of nearby coadsorbates on the surface can dramatically alter these thermodynamic and kinetic parameters causing them to be dependent on species coverage fractions. As there are combinatorially many coadsorbed configurations on the surface, computing the coverage dependence of these parameters is far less straightforward. We present a framework for generating and applying machine learning models to predict coverage-dependent parameters for microkinetic models. Our toolkit enables automatic calculation and evaluation of coadsorbed configurations allowing us to sample 2,000 coadsorbed adsorbates and transition states (TSs) for a diverse set of 9 reactions on Cu(111), a challenging surface, with four possible coadsorbates. This dataset was then used to train subgraph isomorphic decision trees (SIDTs) to predict the stability and association energy of configurations. We were able to achieve mean absolute errors (MAEs) of 0.106 eV on adsorbates, 0.172 eV on TSs, and due to natural error cancellation in SIDTs for relative properties, 0.130 eV on reaction energies and 0.180 eV on activation barriers. In conclusion, we describe how to use these models to predict coverage-dependent corrections for adsorbates and TSs and demonstrate on H*, HO*, and O* comparing the generated SIDT model with an iteratively refined version.

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

The Martini 3 Lipidome: Expanded and Refined Parameters Improve Lipid Phase Behavior

Lipid membranes are central to cellular life. Complementing experiments, computational modeling has been essential in unraveling complex lipid-biomolecule interactions, crucial in both academia and industry. The Martini model, a coarse-grained force field for efficient molecular dynamics simulations, is widely used to study membrane phenomena but has faced limitations, particularly in capturing realistic lipid phase behavior. Here, we present refined Martini 3 lipid models with a mapping scheme that distinguishes lipid tails that differ by just two carbon atoms, enhancing the structural resolution and thermodynamic accuracy of model membrane systems including ternary mixtures. The expanded Martini lipid library includes thousands of models, enabling simulations of complex and biologically relevant systems. These advancements establish Martini as a robust platform for lipid-based simulations across diverse fields.

Lipids↗
DOE OSTI2025

Pseudo-equilibrium theory for extrinsic doping control of the topological semimetal Cd 3 As 2

The standard approach for predicting defect equilibria from first principles assumes that the solid-state system is initially in a thermodynamic equilibrium with the external atomic reservoirs. This “growth step” is then often followed by a temperature quench in a “pseudo-equilibrium” in which some or all defect concentrations are frozen in until only the Fermi level E F remains to be equilibrated. However, this protocol does not account for the possibility of site exchanges which can create important defect redistributions as long as short-range defect migration is kinetically permissible. To model this redistribution, we developed an approach to solve for the non-equilibrium chemical potentials as a function of temperature while maintaining the overall defect stoichiometry. We then apply this approach to the Dirac semimetal Cd 3 As 2 to model extrinsic doping with group 1/11 and 14 elements. Undoped Cd 3 As 2 exhibits an undesirable mismatch between E F and the Dirac point. This unintentional electron doping originates from intrinsic defects and is difficult to overcome through adjustment of synthesis conditions alone. Employing our pseudo-equilibrium modeling, we identify extrinsic doping strategies for realizing doping-balanced Cd 3 As 2 at the relatively low temperatures accessible in thin-film growth of this material.

36 MATERIALS SCIENCE↗
DOE OSTI2025

pH-dependent reactivity of water at MgO(100) and MgO(111) surfaces

Facet-dependent surface charging of metal oxides in water dominates the ion transport behavior across the interface, in turn impacting many natural and industrial processes such as adsorption, the formation and stabilization of nanoparticle suspensions, corrosion, and heterogeneous catalysis. Here we investigated the pH-dependent surface chemistry of two low-index MgO single crystal surfaces, namely MgO(100) and MgO(111), using vibrational sum frequency generation (vSFG) spectroscopy. This allowed us to evaluate facet-dependent pH effects on the hydration and hydroxylation at the solid/aqueous interface and point-of-zero charge (PZC) values. The MgO system is complicated by its thermodynamic instability with respect to Mg(OH) 2 in water at ambient conditions. Here, for both hydroxylated MgO(100) and MgO(111) surfaces, the PZC is found to be around pH ∼ 12, which compares well with reported values for MgO single crystal and nanoparticle surfaces. However, structure specific differences in the molecular water hydrogen bonding network near the surface are evident at mildly acidic pH. To our knowledge, this is the first account of the PZC values for the MgO(111) single-crystal surface, an electrostatically unstable MgO termination that is prone to reconstruction.

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

Electron-scale temperature gradients in kinetic equilibrium: MMS observations and Vlasov–Maxwell solutions

For an ideal gas, the notion of thermodynamic equilibrium requires that the system’s temperature be spatially uniform. We demonstrate that for a collisionless plasma, kinetic equilibrium can be achieved even in the presence of a temperature gradient. Motivated by recent MMS magnetopause observations of thin current layers exhibiting electron-scale temperature gradients, here we present an exact solution to the Vlasov-Maxwell system which accounts for the observed electron temperature gradient and reproduces much of the plasma’s spatial variation along one dimension. Furthermore, the current layer in the Vlasov-Maxwell model is self-consistently supported by a prominent crescent shaped electron velocity distribution that naturally arises due to the temperature transition, which agrees with the electron distribution structures measured by MMS during two asymmetric current layers encountered at the magnetopause.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Experimental electronic structures of the Fe IV =O bond in S=1 heme vs. nonheme sites: Effect of the porphyrin ligand

High-valent Fe IV =O species are common intermediates in biological and artificial catalysts. Heme and nonheme S=1 Fe IV =O sites have been synthesized and studied for decades but little quantitative experimental comparison of their electronic structures has been available, due to the lack of direct methods focused on the iron. This study allows a rigorous determination of the electronic structure of a nonheme Fe IV =O center and its comparison to an Fe IV =O heme site using 1s2p resonant inelastic X-ray scattering (RIXS) and Fe L-edge X-ray absorption spectroscopy (XAS). Further, variable temperature magnetic circular dichroism (VT-MCD) of the ligand field transitions, combined with nuclear resonance vibrational spectroscopy of the two S=1 Fe IV =O systems show that the equatorial ligand field decreases from a nonheme to a heme Fe IV =O site. Alternatively, RIXS and Fe L-edge XAS combined with MCD show that the Fe dπ orbitals are unperturbed in the Fe IV =O heme relative to the nonheme site because the strong axial Fe-O bond uncouples the Fe dπ orbitals from the porphyrin π-system. As a consequence, the thermodynamics and kinetics of the H-atom abstraction reactions are actually very similar for heme compound II and nonheme Fe IV =O active sites.

bioinorganic chemistry↗
DOE OSTI2026

Primary biomolecular adsorption energetics of core–shell nanocomplexes: Implications for biological interactions

Molecular organization at the nano-bio interface governing the colloidal stability, reactivity, immune recognition, and drug delivery performance of nanoparticles remains difficult to predict. Quantifying the primary hydration energetics of biomolecule-coated nanomaterials can determine those interactions and provide a basis for engineered nanocarriers with tailored behavior in biological systems. Here, we measured the thermodynamics of water adsorption on patchy dry magnetite (Fe 3 O 4 ) nanoparticles coated with three model biomolecules, bovine serum albumin, potato starch, and lauric acid and compared these properties to the hydration energetics of the corresponding free dry biomolecules. The results demonstrate how the surface functionalization alters the hydrophilicity, the accessible hydrophilic surface, and the interaction potential of the nanocomplex surface with biological media. The protein coating increases the interaction potential of the surface of the nanocomplex. The weaker interaction potential of the polysaccharide coating and the relatively large hydrophilic surface area allow dynamic and reversible binding, while the fatty acid rearranges into a partial bilayer with very strong hydrophilicity. The findings establish the hydration enthalpy as a quantitative basis to determine and interpret nanoparticle interactions with proteins, membranes, and biological fluids, and provide a thermodynamic foundation for designing nanocarriers with predictable biological reactivity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Thermalization from quantum entanglement: Jet simulations in the massive Schwinger model

We investigate the emergence of thermalization in a quantum-field-theoretic model mimicking the production of jets in QCD: the massive lattice Schwinger model coupled to external sources. Specifically, we compute the expectation values of local operators as functions of time and compare them to their thermal counterparts, quantify the overlap between the evolving density matrix and the thermal one, and compare the dynamics of the energy-momentum tensor to predictions from relativistic hydrodynamics. Through these studies, we find that the system approaches thermalization at late times and elucidate the mechanisms by which quantum entanglement drives thermalization in closed field-theoretic systems. Our results show how thermodynamic behavior emerges in real time from unitary quantum dynamics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Novel Remediation of Mercury-Contaminated Sediments Using Manganese Oxide and Activated Carbon Amendments - 20394

Mercury contamination in sediments poses a challenging remediation problem because of its ability to biogeochemically cycle and biomagnify in the food chain as methylmercury from low resident sediment and porewater concentrations. One mitigation approach is to lower the overall rate of net methylmercury production by creating conditions unfavorable for sulfate and iron reduction or methanogenesis, the primary pathways for methylation of inorganic mercury. Manganese oxides and other redox-active amendments are a promising and novel in situ amendment technology aimed at reducing net methylation of mercury though control of redox conditions in mercury-contaminated sediment-water systems. This paper reviews the thermodynamic and kinetic framework for this remediation technology, presents experimental and spectroscopic data illustrating the approach using solid manganese oxides, and discusses potential advantages, disadvantages, and applications. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Sampling lattices in semi-grand canonical ensemble with autoregressive machine learning

Calculating thermodynamic potentials and observables efficiently and accurately is key for the application of statistical mechanics simulations to materials science. However, naive Monte Carlo approaches, on which such calculations are often dependent, struggle to scale to complex materials in many state-of-the-art disciplines such as the design of high entropy alloys or multi-component catalysts. To address this issue, we adapt sampling tools built upon machine learning-based generative modeling to the materials space by transforming them into the semi-grand canonical ensemble. Furthermore, we show that the resulting models are transferable across wide ranges of thermodynamic conditions and can be implemented with any internal energy model U, allowing integration into many existing materials workflows. We demonstrate the applicability of this approach to the simulation of benchmark systems (AgPd, CuAu) that exhibit diverse thermodynamic behavior in their phase diagrams. Finally, we discuss remaining challenges in model development and promising research directions for future improvements.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Recycling of Lead Pastes from Spent Lead–Acid Batteries: Thermodynamic Constraints for Desulphurization

Lead–acid batteries are important to modern society because of their wide usage and low cost. The primary source for production of new lead–acid batteries is from recycling spent lead–acid batteries. In spent lead–acid batteries, lead is primarily present as lead pastes. In lead pastes, the dominant component is lead sulfate (PbSO 4 , mineral name anglesite) and lead oxide sulfate (PbO•PbSO 4 , mineral name lanarkite), which accounts for more than 60% of lead pastes. In the recycling process for lead–acid batteries, the desulphurization of lead sulfate is the key part to the overall process. In this work, the thermodynamic constraints for desulphurization via the hydrometallurgical route for recycling lead pastes are presented. The thermodynamic constraints are established according to the thermodynamic model that is applicable and important to recycling of lead pastes via hydrometallurgical routes in high ionic strength solutions that are expected to be in industrial processes. The thermodynamic database is based on the Pitzer equations for calculations of activity coefficients of aqueous species. The desulphurization of lead sulfates represented by PbSO 4 can be achieved through the following routes. (1) conversion to lead oxalate in oxalate-bearing solutions; (2) conversion to lead monoxide in alkaline solutions; and (3) conversion to lead carbonate in carbonate solutions. Among the above three routes, the conversion to lead oxalate is environmentally friendly and has a strong thermodynamic driving force. Oxalate-bearing solutions such as oxalic acid and potassium oxalate solutions will provide high activities of oxalate that are many orders of magnitude higher than those required for conversion of anglesite or lanarkite to lead oxalate, in accordance with the thermodynamic model established for the oxalate system. An additional advantage of the oxalate conversion route is that no additional reductant is needed to reduce lead dioxide to lead oxide or lead sulfate, as there is a strong thermodynamic force to convert lead dioxide directly to lead oxalate. As lanarkite is an important sulfate-bearing phase in lead pastes, this study evaluates the solubility constant for lanarkite regarding the following reaction, based on the solubility data, PbO•PbSO 4 + 2H + ⇌ 2Pb 2+ + SO 4 2– + H 2 O(l).

25 ENERGY STORAGE↗
DOE OSTI2021

Thermodynamic non-ideality and disorder heterogeneity in actinide silicate solid solutions

Non-ideal thermodynamics of solid solutions can greatly impact materials degradation behavior. We have investigated an actinide silicate solid solution system (USiO 4 –ThSiO 4 ), demonstrating that thermodynamic non-ideality follows a distinctive, atomic-scale disordering process, which is usually considered as a random distribution. Neutron total scattering implemented by pair distribution function analysis confirmed a random distribution model for U and Th in first three coordination shells; however, a machine-learning algorithm suggested heterogeneous U and Th clusters at nanoscale (~2 nm). The local disorder and nanosized heterogeneous is an example of the non-ideality of mixing that has an electronic origin. Partial covalency from the U/Th 5f–O 2p hybridization promotes electron transfer during mixing and leads to local polyhedral distortions. The electronic origin accounts for the strong non-ideality in thermodynamic parameters that extends the stability field of the actinide silicates in nature and under typical nuclear waste repository conditions.

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

Walking with the Atoms in a Chemical Bond: A Perspective Using Quantum Phase Transition

Phase transitions happen at critical values of the controlling parameters, such as the critical temperature in classical phase transitions, and system critical parameters in the quantum case. However, true criticality happens only at the thermodynamic limit, when the number of particles goes to infinity with constant density. To perform the calculations for the critical parameters, a finite-size scaling approach was developed to extrapolate information from a finite system to the thermodynamic limit. With the advancement in the experimental and theoretical work in the field of ultra-cold systems, particularly trapping and controlling single atomic and molecular systems, one can ask: do finite systems exhibit quantum phase transition? To address this question, finite-size scaling for finite systems was developed to calculate the quantum critical parameters. The recent observation of a quantum phase transition in a single trapped 171 Yb + ion indicates the possibility of quantum phase transitions in finite systems. This perspective focuses on examining chemical processes at ultra-cold temperatures, as quantum phase transitions—particularly the formation and dissociation of chemical bonds—are the basic processes for understanding the whole of chemistry.

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

Organic Rankine Cycle Integration and Optimization for High Efficiency CHP Genset Systems (Final Technical Report)

This project successfully advanced the integration of Organic Rankine Cycle (ORC) technology with reciprocating engine–based combined heat and power (CHP) systems to improve electrical efficiency, total CHP efficiency, and grid-responsive operation. Over three budget periods, the work progressed from high-temperature ORC component development and thermodynamic model validation to next-generation system design, working fluid transition, and techno-economic analysis. Key technical accomplishments include development and validation of a thermodynamic model capable of accurately predicting ORC performance across an expanded temperature and pressure envelope; successful identification and validation of low-global-warming-potential (GWP) working fluids—most notably R1233zd(E)—as viable replacements for R245fa; and demonstration of scalable ORC architectures suitable for integration with 1–20 MW class reciprocating engines. These advances enable flexible CHP configurations that can increase electrical output while maintaining high overall utilization of available thermal energy. The project also produced a clean-sheet design for a next-generation ORC system targeting substantially higher power output per unit, supported by detailed component selection, heat exchanger evaluation, and system-level modeling. Techno-economic analyses indicate that ORC-enabled flexible CHP systems can meet or exceed Department of Energy (DOE) efficiency targets while providing value to both facility operators and the electric grid.. Late-stage testing of the largest next-generation ORC prototype identified limitations related to pump net positive suction head (NPSH) requirements and condenser flooding under certain operating conditions. Although these issues constrained full validation of that configuration within the project timeframe, they provided clear and actionable design guidance for future system refinements. Importantly, validated modeling, smaller-scale testing, and working fluid evaluations confirmed the technical viability of the overall approach. In aggregate, this project met its core objectives by establishing validated design tools, de-risking key ORC technologies for CHP applications, and defining a credible pathway toward commercialization of flexible, high-efficiency CHP systems. The results form a strong foundation for continued development and deployment beyond the conclusion of the DOE-funded effort.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2026

Self‐Limiting Polymerization‐Induced Crystallization‐Driven Self‐Assembly (SL‐PI‐CDSA) Enables Templated Synthesis of Chiral Plasmonic, Hybrid 2D Hexagonal Assemblies

Examples of self-regulating synthetic self-assembly are relatively few, with most known chemical systems relying on kinetic rather than thermodynamic control. Herein, we demonstrate the rapid generation (t ≤ 5 min) of size-tunable ultralow dispersity (Ð ≤ 1.01) 2D hexagonal nanosheets governed by self-limiting self-assembly (SLSA). Self-assembly in natural systems occurs with exquisite control of structure, function, and dimension. We demonstrate that key aspects of biological assembly can be rationally applied toward the development of bottom-up approaches for the construction of chiral nanomaterials. To this end, self-limiting polymerization-induced crystallization-driven self-assembly (SL-PI-CDSA) of modular and templating aryl isocyanide (AIC) monomers yields functional 2D assemblies permitting post-polymerization/assembly modifications. Detailed study of the internal and external structure of the hexagonal nanosheets reveals topological defects which offer mechanistic insights into both their assembly and subsequent utilization. Specifically, these features enable fabrication of chiral hybrid organic–inorganic nanomaterials incorporating chiral plasmonic metal nanoparticles (MNPs). Our results suggest that the synergistic interplay of template-driven confinement and hierarchical chirality induce symmetry breaking of in situ-generated gold MNPs. We anticipate that the platform presented will facilitate fabrication of new hybrid, chiral organic–inorganic nanostructures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Interfacial Effects on Transport Coefficient Measurements in Li-ion Battery Electrolytes

Development of Li + -containing electrolytes with improved transport properties requires reliable, reproducible, and ideally low volume techniques to rigorously understand ion-transport with varying composition. Precisely measuring the complete set of transport coefficients in liquid electrolytes under battery-relevant operating conditions is difficult and the reliability of these methods are sparsely described in electrolyte transport literature. In this work, we apply a potentiostatic polarization-based transport characterization approach typically used for polymer electrolytes to liquid electrolyte systems in an attempt to fully measure all transport coefficients (conductivity, total salt diffusion coefficient, thermodynamic factor and transference number) for the model system of LiPF 6 in an ethylene carbonate—ethyl methyl carbonate (EC:EMC) mixture. Using systematic timescale and statistical analyses, we find that transport coefficients measured using potentiostatic polarization of Li-Li symmetric cells exhibit strong correlation to Li electrode interfacial resistance, indicating that such methods are probing both bulk and interfacial phenomena. This reveals a major roadblock in characterizing electrolyte systems where the interfacial resistance is significantly larger than ohmic electrolyte resistance. As a result, we find that methods that rely on potentiostatic Li metal stripping/plating do not readily result in reliable liquid electrolyte transport coefficients, unlike similar methods for solid polymer electrolytes, where interfacial resistances are typically smaller than electrolyte resistances at the elevated temperatures typically of interest for such electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2020

Machine learning approaches for structural and thermodynamic properties of a Lennard-Jones fluid

Predicting the functional properties of many molecular systems relies on understanding how atomistic interactions give rise to macroscale observables. However, current attempts to develop predictive models for the structural and thermodynamic properties of condensed-phase systems often rely on extensive parameter fitting to empirically selected functional forms whose effectiveness is limited to a narrow range of physical conditions. Here, we illustrate how these traditional fitting paradigms can be superseded using machine learning. Specifically, we use the results of molecular dynamics simulations to train machine learning protocols that are able to produce the radial distribution function, pressure, and internal energy of a Lennard-Jones fluid with increased accuracy in comparison to previous theoretical methods. The radial distribution function is determined using a variant of the segmented linear regression with the multivariate function decomposition approach developed by Craven et al. [J. Phys. Chem. Lett. 11, 4372 (2020)]. The pressure and internal energy are determined using expressions containing the learned radial distribution function and also a kernel ridge regression process that is trained directly on thermodynamic properties measured in simulation. The presented results suggest that the structural and thermodynamic properties of fluids may be determined more accurately through machine learning than through human-guided functional forms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗