Search NASA⌕ Search

SEARCH · Search NASA

Results for “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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Classifying thermodynamic cloud phase using machine learning models

Vertically resolved thermodynamic cloud-phase classifications are essential for studies of atmospheric cloud and precipitation processes. The Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Thermodynamic Cloud Phase (THERMOCLDPHASE) value-added product (VAP) uses a multi-sensor approach to classify the thermodynamic cloud phase by combining lidar backscatter and depolarization, radar reflectivity, Doppler velocity, spectral width, microwave-radiometer-derived liquid water path, and radiosonde temperature measurements. The measured pixels are classified as ice, snow, mixed phase, liquid (cloud water), drizzle, rain, and liq_driz (liquid+drizzle). We use this product as the ground truth to train three machine learning (ML) models to predict the thermodynamic cloud phase from multi-sensor remote sensing measurements taken at the ARM North Slope of Alaska (NSA) observatory: a random forest (RF), a multi-layer perceptron (MLP), and a convolutional neural network (CNN) with a U-Net architecture. Evaluations against the outputs of the THERMOCLDPHASE VAP with 1 year of data show that the CNN outperforms the other two models, achieving the highest test accuracy, F1 score, and mean intersection over union (IOU). Analysis of ML confidence scores shows that ice, rain, and snow have higher confidence scores, followed by liquid, while mixed, drizzle, and liq_driz have lower scores. Feature importance analysis reveals that the mean Doppler velocity and vertically resolved temperature are the most influential data streams for ML thermodynamic cloud-phase predictions. Lidar measurements exhibit lower feature importance due to rapid signal attenuation caused by the frequent presence of persistent low-level clouds at the NSA site. The ML models' generalization capacity is further evaluated by applying them at another Arctic ARM site in Norway using data taken during the ARM Cold-Air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) field campaign. The models demonstrated similar performance to that observed at the NSA site. Finally, we evaluate the ML models' response to simulated instrument outages and signal degradation and show that a CNN U-Net model trained with input channel dropouts performs better when input fields are missing.

ARM Aerial Facility↗
DOE OSTI2025

Inverse Thermodynamics: Designing Interactions for Targeted Phase Behavior

The traditional goal of inverse self-assembly is to design interactions that drive particles toward a desired target structure. However, achieving successful self-assembly also requires tuning the thermodynamic conditions under which the structure is stable. In this work, we extend the inverse design paradigm to explicitly address this challenge by developing a framework for inverse thermodynamics, i.e., the design of interaction potentials that realize specific thermodynamic behavior. As a step in this direction, using patchy particle mixtures as a model system, we demonstrate how precise control over both bonding topology and bond energetics enables the programming of targeted phase behavior. In particular, we establish design principles for azeotropic demixing and show how to create mixtures that exhibit azeotropy at any prescribed composition. Our predictions are validated through Gibbs-ensemble simulations [Panagiotopoulos, Mol. Phys. 1987, 61, 813−826]. These results highlight the necessity of coupling structural design with thermodynamic engineering, and provide a blueprint for controlling complex phase behavior in multicomponent systems.

Azeotropes↗
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 OSTI2025

A thermodynamically consistent discretization of 1D thermal-fluid models using their metriplectic 4-bracket structure

Thermodynamically consistent models in continuum physics, i.e. models which satisfy the first and second laws of thermodynamics, may be expressed using the metriplectic formalism. In this work, we leverage the structures underlying this modeling formalism to preserve thermodynamic consistency in discretizations of a fluid model. The procedure relies (1) on ensuring that the spatial semi-discretization retains certain symmetries and degeneracies of the Poisson and metriplectic 4-brackets, and (2) on the use of an appropriate energy conserving time-stepping method. Here, the minimally simple yet nontrivial example of a one-dimensional thermal-fluid model is treated. It is found that preservation of the requisite symmetries and degeneracies of the 4-bracket is relatively simple to ensure in Galerkin spatial discretizations, suggesting a path forward for thermodynamically consistent discretizations of more complex fluid models using more specialized Galerkin methods.

Hamiltonian structure↗
DOE OSTI2024

Thermodynamically consistent Cahn–Hilliard–Navier–Stokes equations using the metriplectic dynamics formalism

Cahn–Hilliard–Navier–Stokes (CHNS) systems describe flows with two-phases, e.g., a liquid with bubbles. Obtaining constitutive relations for general dissipative processes for such systems, which are thermodynamically consistent, can be a challenge. We show how the metriplectic 4-bracket formalism (Morrison and Updike, 2024) achieves this in a straightforward, in fact algorithmic, manner. First, from the noncanonical Hamiltonian formulation for the ideal part of a CHNS system we obtain an appropriate Casimir to serve as the entropy in the metriplectic formalism that describes the dissipation (e.g. viscosity, heat conductivity and diffusion effects). General thermodynamics with the concentration variable and its thermodynamics conjugate, the chemical potential, are included. Having expressions for the Hamiltonian (energy), entropy, and Poisson bracket, we describe a procedure for obtaining a metriplectic 4-bracket that describes thermodynamically consistent dissipative effects. The 4-bracket formalism leads naturally to a general CHNS system that allows for anisotropic surface energy effects. Furthermore, this general CHNS system reduces to cases in the literature, to which we can compare.

Cahn–Hilliard↗
DOE OSTI2025

Proton, Electron, and Hydrogen-Atom Transfer Thermodynamics of the Metal–Organic Framework, Ti-MIL-125, Are Intrinsically Correlated to the Structural Disorder

Interfacial charge transfer reactions involving protons and/or electrons are fundamental to heterogeneous catalysis and many other reactions relevant to energy, chemical, and biological sectors. Metal–organic frameworks (MOFs) with redox-active metal-oxo nodes have emerged as candidate materials to examine these reactions with near-atomic-level precision, given their crystalline nature. Here, we employed a colloidally stable, Ti-based MOF, Ti-MIL-125, with different crystal sizes to examine catalytically relevant charge transfer thermodynamics. The Ti 8 (μ 2 -O) 8 (μ 2 -OH) 4 nodes structurally mimic TiO 2 , which has shown some PCET reactivity toward reactions of H 2 , O 2 , and others. In this report, we have demonstrated that a change in crystal size induces different amounts of structural disorder to the Ti-oxo node, further changing the thermodynamics of proton/electron/hydrogen-atom transfer reactions. Using electrochemical open-circuit potential (E OCP ) measurements, we have determined that all crystallites undergo a 1H + /1e – redox reaction, which, given the stoichiometry, can be considered as a net H atom transfer (HAT) reaction. The thermodynamics of this HAT reaction, the Ti 3+ O–H bond dissociation free energy (BDFE), was dependent on the crystal size of the MOF, as the decrease in crystal size induced more structural disorder. Our computational calculations have indicated that this difference in BDFE is due to a local change in the geometry of Ti cations, rather than the commonly invoked defects, such as the “missing-linker” defect sites. Individual proton/electron transfer (PT/ET) thermodynamics were also highly dependent on the crystal sizes. These were probed using pK a or band gaps (E g ), respectively. These findings suggest that, particularly when MOFs are nanosized with a large amount of structural disorder, they should no longer be considered “true” single-site catalysts; this is an implicit, but widespread assumption within the MOF-based catalysis field. Implications of these findings will be contrasted with structurally similar metal oxides like TiO 2 and other redox-active MOFs.

Bond dissociation free energy↗
DOE OSTI2024

Pyrophosphate-free glycolysis in Clostridium thermocellum increases both thermodynamic driving force and ethanol titers

Background: Clostridium thermocellum is a promising candidate for production of cellulosic biofuels, however, its final product titer is too low for commercial application, and this may be due to thermodynamic limitations in glycolysis. Previous studies in this organism have revealed a metabolic bottleneck at the phosphofructokinase (PFK) reaction in glycolysis. In the wild-type organism, this reaction uses pyrophosphate (PP i ) as an energy cofactor, which is thermodynamically less favorable compared to reactions that use ATP as a cofactor. Previously we showed that replacing the PP i -linked PFK reaction with an ATP-linked reaction increased the thermodynamic driving force of glycolysis, but only had a local effect on intracellular metabolite concentrations, and did not affect final ethanol titer. Results: In this study, we substituted PP i -pfk with ATP-pfk, deleted the other PPi-requiring glycolytic gene pyruvate:phosphate dikinase (ppdk), and expressed a soluble pyrophosphatase (PPase) and pyruvate kinase (pyk) genes to engineer PP i -free glycolysis in C. thermocellum. We demonstrated a decrease in the reversibility of the PFK reaction, higher levels of lower glycolysis metabolites, and an increase in ethanol titer by an average of 38% (from 15.1 to 21.0 g/L) by using PP i -free glycolysis. Conclusions: By engineering PP i -free glycolysis in C. thermocellum, we achieved an increase in ethanol production. These results demonstrate that optimizing the thermodynamic landscape through metabolic engineering can enhance product titers. While further increases in ethanol titers are necessary for commercial application, this work represents a significant step toward engineering glycolysis in C. thermocellum to increase ethanol titers.

09 BIOMASS FUELS↗
DOE OSTI2025

Thermodynamic Cloud Phase Classifications Using Machine Learning at NSA and ANX

Vertically resolved thermodynamic cloud phase classifications are essential for studies of atmospheric cloud and precipitation processes. The Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) THERMOCLDPHASE Value-Added Product (VAP) uses a multi-sensor approach to classify thermodynamic cloud phase by combining lidar backscatter and depolarization, radar reflectivity, Doppler velocity, spectral width, microwave radiometer-derived liquid water path, and radiosonde temperature measurements. The measured voxels are classified as ice, snow, mixed-phase, liquid (cloud water), drizzle, rain, and liq_driz (liquid+drizzle). We use this product as the ground truth to train three machine learning (ML) models to predict the thermodynamic cloud phase from multi-sensor remote sensing measurements taken at the ARM North Slope of Alaska (NSA) observatory: a random forest (RF), a multilayer perceptron (MLP), and a convolutional neural network (CNN) with a U-Net architecture. Evaluations against the outputs of the THERMOCLDPHASE VAP with one year of data show that the CNN outperforms the other two models, achieving the highest test accuracy, F1-score, and mean Intersection over Union (IOU). Analysis of ML confidence scores shows ice, rain, and snow have higher confidence scores, followed by liquid, while mixed, drizzle, and liq_driz have lower scores. Feature importance analysis reveals that the mean Doppler velocity and vertically resolved temperature are the most influential datastreams for ML thermodynamic cloud phase predictions. The ML models’ generalization capacity is further evaluated by applying them at another Arctic ARM site in Norway using data taken during the ARM Cold-Air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) field campaign. Finally, we evaluate the ML models’ response to simulated instrument outages and signal degradation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Insights into controlling bacterial cellulose nanofiber film properties through balancing thermodynamic interactions and colloidal dynamics

In recent years, nanocellulose has emerged as a sustainable and environmentally friendly alternative to traditional petroleum-derived structural polymers. Sourced either from plants, algae, or bacteria, nanocellulose can be processed into colloid, gel, film and fiber forms. However, the required fundamental understanding of process parameters that govern the morphology and structure–property relationships of nanocellulose systems, from colloidal suspensions to bulk materials, has not been developed and generalized for all forms of cellulose. This further hinders the more widespread adoption of this biopolymer in applications. Our study investigates the dispersion of cellulose nanofibers (CNFs) produced by a bacterial–yeast co-culture, in solvents, highlighting the role of thermodynamic interactions in influencing their colloidal behavior. By adjusting Hansen solubility parameters, we controlled the thermodynamic relationship between CNFs and solvents across various concentrations, studying the dilute to semi-dilute regimes. Rheological measurements revealed that the threshold at which a concentration-based regime transition occurs is distinctly solvent-dependent. Complementing rheological analysis with small angle X-ray scattering and zeta potential measurements, our findings reveal that enhancing CNF–solvent interactions increases excluded volume in the dilute regime, emphasizing the importance of the balance between fiber–fiber and fiber–solvent interactions. Moreover, we investigated the transition from colloidal to solid state by creating films from dispersions with varying interaction parameters in semi-dilute regimes. Through mechanical testing and scanning electron microscopy imaging of the fracture surfaces, we highlight the significance of electrokinetic effects in such transitions, as dispersions with higher electrokinetic stabilization gave rise to stronger and tougher films despite having less favorable thermodynamic interaction parameters. Finally, our work provides insights into the thermodynamic and electrokinetic interplay that governs bacterial CNF dispersion, offering a foundation for future application and a deeper understanding of nanocellulose's colloidal and structure-property relationships.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Computing chemical potentials with machine-learning-accelerated simulations to accurately predict thermodynamic properties of molten salts

The successful design and deployment of next-generation nuclear technologies heavily rely on thermodynamic data for relevant molten salt systems. However, the lack of accurate force fields and efficient methods has limited the quality of thermodynamic predictions from atomistic simulations. Here we propose an efficient free energy framework for computing chemical potentials, which is the central free energy quantity behind many thermodynamic properties. We accelerate our simulations without sacrificing accuracy by using machine learning interatomic potentials trained on density functional theory (DFT) data. Using lithium chloride as our model system, we compute chemical potentials with DFT-accuracy for solid and liquid phases by transmuting ions into noninteracting particles. Notably, in the liquid phase, we demonstrate consistency whether we transmute one ion pair or the entire system into ideal gas particles. By locating the temperature where the chemical potential of solid and liquid phases cross, we predict a melting point of 880 ± 18 K for lithium chloride, which is remarkably close to the experimental value of 883 K. With this successful demonstration, we lay the foundation for high-throughput thermodynamic predictions of many properties that can be derived from the chemical potentials of the minority and majority components in molten salts.

Gibson, Luke D. [Oak Ridge National Laboratory (OR↗
DOE OSTI2024

RASPA3

RASPA3, a molecular simulation code for computing adsorption and diffusion in nanoporous materials and thermodynamic and transport properties of fluids. It implements force field based classical Monte Carlo/molecular dynamics in various ensembles. RASPA3 is rewritten from the ground up in C++23 with speed and code readability in mind. Transition-matrix Monte Carlo is added to compute the density of states and free energies. The Monte Carlo code for rigid molecules is based on quaternions, and the atomic positions needed in the energy evaluation are recreated from the center of mass position and quaternion orientation. The expanded ensemble methodology for fractional molecules, with a scaling parameter λ between 0 and 1, now also keeps track of analytic expressions of dU/dλ, allowing independent verification of the chemical potential using thermodynamic integration. The source code is freely available under the MIT license on GitHub.

Dubbeldam, David↗
DOE OSTI2025

Thermodynamics of calcium binding to heparin: Implications of solvation and water structuring for polysaccharide biofunctions

Heparan sulfates are found in all animal tissues and have essential roles in living systems. This family of biomacromolecules modulates binding to calcium ions (Ca 2+ ) in low free energy reactions that influence biochemical processes from cell signaling and anticoagulant efficacy to biomineralization. Despite their ubiquity, the thermodynamic basis for how heparans and similarly functionalized biomolecules regulate Ca 2+ interactions is not yet established. Using heparosan (Control) and heparins with different positions of sulfate groups, we quantify how SO 3 − and COO − content and SO 3 − position modulate Ca 2+ binding by isothermal titration calorimetry. The free energy of all heparin-Ca2+ interactions (ΔG rxn ) is dominated by entropic contributions due to favorable water release from polar, hydrophilic groups. Heparin with both sulfate esters (O-SO 3 − ) and sulfamides (N-SO 3 − ) has the strongest binding to Ca 2+ compared to heparosan and to heparin with only O-SO 3 − groups (~3X). By linking Ca 2+ binding thermodynamics to measurements of the interfacial energy for calcite (CaCO 3 ) crystallization onto polysaccharides, we show molecule-specific differences in nucleation rate can be explained by differences in water structuring during Ca 2+ interactions. A large entropic term (-TΔS rxn ) upon Ca 2+ –polysaccharide binding correlates with high interfacial energy to CaCO 3 nucleation. Combining our measurements with literature values indicates many Ca 2+ –polysaccharide interactions have a shared thermodynamic signature. The resulting enthalpy–entropy compensation relationship suggests these interactions are generally dominated by water restructuring involving few conformational changes, distinct from Ca 2+ –protein binding. Our findings quantify the thermodynamic origins of heparin-specific interactions with Ca 2+ and demonstrate the contributions of solvation and functional group position during biomacromolecule-mediated ion regulation.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

M3SF-24LL010302062-NEA-TDB Management and International Collaborations in Sorption and Thermodynamic Modeling

This progress report (Level 3 Milestone Number M3SF-24LL010302062) summarizes research conducted at Lawrence Livermore National Laboratory (LLNL) within the Crystalline International Collaborations Work Package Number SF-24LL01030206. The activity is focused on our long-term commitment to engaging our partners in international nuclear waste repository research. This includes participation in the Nuclear Energy Agency Thermochemical Database (NEA-TDB) Project and development of methodologies for integrating US and international thermodynamic databases for use in SFWST Generic Disposal System Assessment (GDSA) efforts. A continuing focus for FY24 efforts is to support the US participation in the NEA-TDB effort. The focus of FY24 activities was the development of an agreement for a Phase 7 activity that will start in Q1 of 2025. Mavrik Zavarin is now the US representative on both the Management Board and the Executive group to the NEA-TDB. He is also the POC for the Cements State of the Art Report that is undergoing peer review in FY24. In FY24, we used our position on the NEA-TDB MB and EG to facilitate the integration of NEA-TDB thermochemical data with LLNL’s SUPCRTNE thermodynamic database that supports the SFWST GDSA activities. This effort is coordinated with the Argillite work package SUPCRTNE database development efforts (Wolery, 2024). The goal is to provide a downloadable database that will be hosted on LLNL’s thermodynamics website which incorporates NEA-TDB data into the LLNL database where appropriate. We also began engagement with the EURAD2 program that was initiated in FY24 by our European collaborators at the Karlsruhe Institute of Technology (KIT), Germany. The primary focus of the engagement is with WP20: DITUSC Thermodynamic database evaluation program. A kickoff meeting for this activity is planned for early FY25. Finally, we have been selected to co-host (with Clemson University) the International Conference on Chemistry and Migration Behaviour of Actinides and Fission Products in the Geosphere in 2025 (Migration2025). The meeting will be held September 21-26, 2025, in New Orleans, Louisiana, and will focus on international efforts to understand the risks of radionuclide releases into the environment. This central focus of this conference is on international efforts to develop safe disposal options for nuclear wastes. As such, we are developing a theme focused on US underground nuclear waste repository science.

58 GEOSCIENCES↗
DOE OSTI2025

CASM Monte Carlo: Calculations of the thermodynamic and kinetic properties of complex multicomponent crystals

Monte Carlo techniques play a central role in statistical mechanics approaches that connect macroscopic thermodynamic and kinetic properties to the electronic structure of a material. This paper describes the implementation of Monte Carlo techniques for the study of multicomponent crystalline materials within the Clusters Approach to Statistical Mechanics (CASM) software suite, and demonstrates their use in model systems to calculate free energies and kinetic coefficients, study phase transitions, and construct phase diagrams from first principles. Many crystal structures are complex, with multiple sublattices occupied by differing sets of chemical species, along with the presence of vacancies or interstitial species. This imposes constraints on concentration variables, the form of thermodynamic potentials, and the values of kinetic transport coefficients. The framework used by CASM to formulate thermodynamic potentials and kinetic transport coefficients accounting for arbitrarily complex crystal structures is presented and demonstrated with examples of increasing complexity. Additionally, an overview of the capabilities of the CASM software specific to Monte Carlo methods is given, and a new CASM software package is introduced, casm-flow, which helps automate the setup, submission, management, and analysis of Monte Carlo simulations.

Cluster expansion↗
DOE OSTI2025

Detecting thermodynamic phase transition via explainable machine learning of photoemission spectroscopy

Identifying thermodynamic signatures of electronic phases, such as superconductivity, is challenging in low-dimensional materials due to strong fluctuations and low probing volume. Spectroscopic methods are often used to identify new bulk phases, but their main measurable quantity—electronic energy gaps—is no longer an effective order parameter in low-dimensional and fluctuating systems. Combining angle-resolved photoemission with a domain-adversarial neural network, we report a data-driven method to identify thermodynamic phase transitions solely based on single-particle spectra. We demonstrate 97.6% accuracy in cuprate superconductor Bi 2 Sr 2 CaCu 2 O 8+δ with strong superconducting fluctuations. This model notably compensates for the scarcity of experimental data by leveraging virtually inexhaustible simulated data. Further, its explainability reveals the crucial role of in-gap spectral weight in detecting phase fluctuations and thermodynamic transitions. Our work pinpoints the spectroscopic signatures of fluctuating orders and enables using spectroscopy for machine-learning-assisted material discovery for low-dimensional and strong coupling systems.

2D materials↗
DOE OSTI2024

Adsorption Thermodynamics for Process Simulation

Adsorption has rapidly evolved in recent decades and is an established separation technology extensively practiced in gas separation industries and others. However, rigorous thermodynamic modeling of multicomponent adsorption equilibrium remains elusive, and industrial practitioners rely heavily on expensive and time-consuming trial-and-error pilot studies to develop adsorption units. Here, this article highlights the need for rigorous adsorption thermodynamic models and the limitations and deficiencies of existing models such as the extended Langmuir isotherm, dual-process Langmuir isotherm, and adsorbed solution theory. It further presents a series of recent advances in the generalization of the classical Langmuir isotherm of single-component adsorption by deriving an activity coefficient model to account for the adsorbed phase adsorbate–adsorbent interactions, substituting adsorbed phase adsorbate and vacant site concentrations with activities, and extending to multicomponent competitive adsorption equilibrium, both monolayer and multilayer. Requiring a minimum set of physically meaningful model parameters, the generalized Langmuir isotherm for monolayer adsorption and the generalized Brunauer–Emmett–Teller isotherm for multilayer adsorption address various thermodynamic modeling challenges including adsorbent surface heterogeneity, isosteric enthalpies of adsorption, BET surface areas, adsorbed phase nonideality, adsorption azeotrope formation, and multilayer adsorption. Also discussed is the importance of quality adsorption data that cover sufficient temperature, pressure, and composition ranges for reliable determination of the model parameters to support adsorption process simulation, design, and optimization.

09 BIOMASS FUELS↗
DOE OSTI2024

Role of Metal–Organic Framework Topology on Thermodynamics of Polyoxometalate Encapsulation

Polyoxometalates (POMs) are discrete anionic clusters whose rich redox properties, strong Brønsted acidity, and high availability of active sites poise them as potent catalysts for oxidation reactions. Here, metal–organic frameworks (MOFs) have emerged as tunable, porous platforms to immobilize POMs, thus increasing their solution stability and catalytic activity. While POM@MOF composite materials have been widely used for a variety of applications, little is known about the thermodynamics of the encapsulation process. Here, we utilize an up-and-coming technique in the field of heterogeneous materials, isothermal titration calorimetry (ITC), to obtain full thermodynamic profiles (ΔH, ΔS, ΔG, K a ) of POM binding. Six different 8-connected hexanuclear Zr-MOFs were investigated to determine the impact of MOF topology (csq, scu, the) on POM encapsulation thermodynamics.

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

First-Principles Insights into the Thermodynamics of Variable-Temperature Ammonia Synthesis on Transition-Metal-Doped Cu (100) and (111)

Ammonia (NH 3 ) is one of the most produced chemicals worldwide. NH 3 synthesis predominantly utilizes the Haber–Bosch (HB) process, requiring high temperatures and pressures. Despite significant process advances, ample opportunity remains for improving the rate, selectivity, catalyst stability, and energy efficiency. Inspired by a recently developed programmable heating and quenching (PHQ) technique, we present in this paper a first-principles screening of candidate single-atom alloy catalysts generated from doping (111) and (100) surfaces of copper (Cu), an ineffective HB catalyst in its pure form. We predict the thermodynamics of two rate-limiting reactions, N 2 dissociative adsorption and the final hydrogenation step leading up to NH 3 release, at 400 and 900 K. Thermodynamically, the former reaction is favored at low temperatures, while the latter is favored at high temperatures. Vanadium-, chromium-, and molybdenum-doped Cu surfaces, due to intermediate M–N covalent bonding character, emerge as appealing candidate catalysts for PHQ NH 3 synthesis, as they balance the thermodynamics of the above-mentioned reaction steps at their respective optimal temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗