Search NASA⌕ Search

SEARCH · Search NASA

Results for “Solvated Molecular Catalysts”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Characterize and control the multi-range structure of solution-phase systems with resonant x-ray scattering

This work aims at demonstrating the feasibility and impact of Anomalous X-ray Scattering (AXS) in characterizing the multi-range structures of solution-state systems. We will show preliminary investigation of long-range correlations in concentrated aqueous electrolytes with a combination of AXS and molecular dynamics (MD) simulations. We will also start exploring the capability of AXS to capture intra- and inter-molecular structure of dilute molecular systems, and specifically its sensitivity to the chemical environment surrounding active metal sites. We anticipate that the development of this method will help us to understand ion solvation and transport, which affect the performances of, for instance, electrocatalytic cells, as well as to identify/control the intrinsic structural factors that lead to selectivity, efficiency, and stability control during catalysis.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effect of Framework Composition and NH 3 on the Diffusion of Cu + in Cu-CHA Catalysts Predicted by Machine-Learning Accelerated Molecular Dynamics

Cu-exchanged zeolites rely on mobile solvated Cu + cations for their catalytic activity, but the role of the framework composition in transport is not fully understood. Ab initio molecular dynamics simulations can provide quantitative atomistic insight but are too computationally expensive to explore large length and time scales or diverse compositions. We report a machine-learning interatomic potential that accurately reproduces ab initio results and effectively generalizes to allow multinanosecond simulations of large supercells and diverse chemical compositions. Biased and unbiased simulations of [Cu(NH 3 ) 2 ] + mobility show that aluminum pairing in eight-membered rings accelerates local hopping and demonstrate that increased NH 3 concentration enhances long-range diffusion. The probability of finding two [Cu(NH 3 ) 2 ] + complexes in the same cage, which is key for SCR-NOx reaction, increases with Cu content and Al content but does not correlate with the long-range mobility of Cu + . Supporting experimental evidence was obtained from reactivity tests of Cu-CHA catalysts with a controlled chemical composition.

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

Solvent structure and dynamics over Brønsted acid MWW zeolite nanosheets

In the liquid phase of heterogeneous catalysis, solvent plays an important role and governs the kinetics and thermodynamics of a reaction. Although it is often difficult to quantify the role of the solvent, it becomes particularly challenging when a zeolite is used as the catalyst. This difficulty arises from the complex nature of the liquid/zeolite interface and the different solvation environments around catalytically active sites. Here, we use ab initio molecular dynamics simulations to probe the local solvation structure and dynamics of methanol and water over MWW zeolite nanosheets with varying Brønsted acidity. We find that the zeolite framework and the number and location of the acid sites in the zeolite influence the structure and dynamics of the solvent. In particular, methanol is more likely to be in the vicinity of the aluminum (Al 3+ ) at the T4 site than at T1 due to easy accessibility. The methanol oxygen binds strongly to the Al at the T4 site, weakening the Al–O for the bridging acid site, which results in the formation of the silanol group, significantly reducing the acidity of the site. The behavior of methanol is in direct contrast to that of water, where protons can easily propagate from the zeolite to the solvent molecules regardless of the acid site location. Furthermore, our work provides molecular-level insights into how solvent interacts with zeolite surfaces, leading to an improved understanding of the catalytic site in the MWW zeolite nanosheet.

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

Prediction of hydration energies of adsorbates at Pt(111) and liquid water interfaces using machine learning

Aqueous phase heterogeneous catalysis is important to various industrial processes, including biomass conversion, Fischer–Tropsch synthesis, and electrocatalysis. Accurate calculation of solvation thermodynamic properties is essential for modeling the performance of catalysts for these processes. Explicit solvation methods employing multiscale modeling, e.g., involving density functional theory and molecular dynamics have emerged for this purpose. Although accurate, these methods are computationally intensive. This study introduces machine learning (ML) models to predict solvation thermodynamics for adsorbates on a Pt(111) surface, aiming to enhance computational efficiency without compromising accuracy. In particular, ML models are developed using a combination of molecular descriptors and fingerprints and trained on previously published water–adsorbate interaction energies, energies of solvation, and free energies of solvation of adsorbates bound to Pt(111). These models achieve root mean square error values of 0.09 eV for interaction energies, 0.04 eV for energies of solvation, and 0.06 eV for free energies of solvation, demonstrating accuracy within the standard error of multiscale modeling. Feature importance analysis reveals that hydrogen bonding, van der Waals interactions, and solvent density, together with the properties of the adsorbate, are critical factors influencing solvation thermodynamics. Furthermore, these findings suggest that ML models can provide rapid and reliable predictions of solvation properties. This approach not only reduces computational costs but also offers insights into the solvation characteristics of adsorbates at Pt(111)–water interfaces.

Adsorption↗
DOE OSTI2023

Unexpected Kinetic Solvent Effects Enhance Activity and Selectivity in Biphasic Systems

Biphasic dehydration of fructose to 5-hydroxymethylfurfural (HMF) has shown unprecedented increases in productivity, but a mechanistic understanding is lacking. Herein, we couple fast experimental reaction kinetics, multiscale modeling (phase behavior, classical molecular dynamics(MD), and quantum mechanics/molecular mechanics MD), in situ sampling, and IR and 13 C-NMR spectroscopy to elucidate the complex effects of nonpolar extracting organic solvents on the kinetics of fructose dehydration. We show that these organic solvents can reach significant mutual solubility with water at reaction temperatures, enabling the partition of the sugar and catalyst into the extracting phase. In the organic-rich environment, the dehydration of fructose proceeds faster and more selectively than in water due to increased relative abundance of the reactive furanose isomer, enhanced water–catalyst–substrate interactions driven by nanophase separation, and higher product stability stemming from preferential solvation. Furthermore, we demonstrate that these solvent effects impact other critical biphasic reactions in biomass upgrading and provide qualitative principles for solvent selection.

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

Stimulus-Responsive Modulation of Solvation Environments in Solid Catalysts

Liquid environments play a crucial role in the biological processes occurring in living organisms as well as in many human-made processes involving electrochemistry, photo-, and thermocatalysis. In the majority of these systems, aqueous phases are ubiquitous due to water’s natural abundance. Water molecules, however, can exert large changes in the chemical environment of catalytically active sites, altering the reaction rates, selectivity, and catalyst stability. These solvation effects induced by water molecules near catalytic sites can drastically change the energy landscape and unlock unique reaction pathways with far more favorable kinetics. In nature, living organisms couple these complex interactions with detection, communication, and actuation mechanisms to induce self-regulatory behavior, ensuring stability of the system and thus long-term durability. Extrapolating this behavior to heterogeneous catalysis is desirable because the resulting “smart materials” can potentially unlock new chemical conversion processes with higher atom efficiency, rates, and stability. The combination of polymer chemistry and heterogeneous catalysis has introduced versatile approaches for creating materials that can respond to cues in the reaction medium that alter the accessibility, intrinsic activity, and selectivity of the catalyst. To achieve this, one could combine stimulus-responsive polymers, which undergo a large volumetric phase transition in response to an external stimulus, with a solid catalyst. This chemo-mechanical response has been employed to create a variety of nanoreactor vessels with stimulus-responsive character that turn on- and off- depending on the reaction conditions. In this Account, we focus on the impact of these polymer coatings on the solvation environment around the active site and the implications of these effects on the reaction energy landscape, molecular arrangement of the solvent, electric fields at the catalyst–liquid interface, binding energy, and mobility of surface reaction intermediates. These seemingly subtle changes in solvent molecules induced by the presence of polymers can have a tremendous impact on the development of bioinspired heterogeneous catalysts, reliable chemical clocks, micro/nanoreactors, and robots. The large library of polymer chemistries offers a plethora of combinations of stimulus-responsive mechanisms (e.g., temperature, pH, light, magnetic field, solvent composition), providing the possibility of creating homeostatic catalysts à la carte.

catalysts↗
DOE OSTI2024

Dynamic Copper Site Redispersion through Atom Trapping in Zeolite Defects

Single-site copper-based catalysts have shown remarkable activity and selectivity for a variety of reactions. However, deactivation by sintering in high-temperature reducing environments remains a challenge and often limits their use due to irreversible structural changes to the catalyst. Here, we report zeolite-based copper catalysts in which copper oxide agglomerates formed after reaction can be repeatedly redispersed back to single sites using an oxidative treatment in air at 550 °C. Under different environments, single-site copper in Cu–Zn–Y/deAlBeta undergoes dynamic changes in structure and oxidation state that can be tuned to promote the formation of key active sites while minimizing deactivation through Cu sintering. For example, single-site Cu 2+ reduces to Cu 1+ after catalyst pretreatment (270 °C, 101 kPa H 2 ) and further to Cu0 nanoparticles under reaction conditions (270–350 °C, 7 kPa EtOH, 94 kPa H 2 ) or accelerated aging (400–450 °C, 101 kPa H 2 ). After regeneration at 550 °C in air, agglomerated CuO was dispersed back to single sites in the presence and absence of Zn and Y, which was verified by imaging, in situ spectroscopy, and catalytic rate measurements. Ab initio molecular dynamics simulations show that solvation of CuO monomers by water facilitates their transport through the zeolite pore, and condensation of the CuO monomer with a fully protonated silanol nest entraps copper and reforms the single-site structure. Importantly, the capability of silanol nests to trap and stabilize copper single sites under oxidizing conditions could extend the use of single-site copper catalysts to a wider variety of reactions and allows for a simple regeneration strategy for copper single-site catalysts.

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

Effects of surface charge and cluster size on the electrochemical dissolution of platinum nanoparticles using COMB3 and continuum electrolyte models

We study the site-dependent dissolution of platinum nanoparticles under electrochemical conditions to assess their thermodynamic stability as a function of shape and size using empirical molecular dynamics and electronic-structure models. The third-generation charge optimized many-body potential is employed to determine the validity of uniform spherical representations of the nanoparticles in predicting dissolution potentials (the Kelvin model). To understand the early stages of catalyst dissolution, implicit solvation techniques based on the self-consistent continuum solvation method are applied. It is demonstrated that interfacial charge and polarization can shift the dissolution energies by amounts on the order of 0.74 eV depending on the surface site and nanoparticle shape, leading to the unexpected preferential removal of platinum cations from highly coordinated sites in some cases.

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

Effect of Solvent on the Local Structure, Dynamics, and Vibrational Density of States in Sn-BEA Zeolite

Lewis acid zeolites are attractive catalysts for epoxidation and biomass valorization, as they are highly active and selective in the liquid phase and can operate at or near ambient conditions. While a rich experimental literature exists on liquid-phase Lewis acid zeolite catalysis, our understanding of the molecular organization and solvent dynamics in the vicinity of Lewis acid sites with differing metal site speciation remains limited. In this work, we investigate the molecular coordination and diffusion of two common solvents (methanol and water) around the closed and open Sn-BEA zeolite active sites using molecular dynamics simulations with a machine-learned interatomic potential trained on ab initio molecular dynamics trajectories. Molecular dynamics simulations reveal that introducing active sites significantly enhances local order in the first and second solvation shells compared to the pure silica case. For methanol, both closed and open active sites are singly coordinated, while more than two water molecules coordinate the open site. In contrast to methanol, we observed that water molecules dissociate, leading to the formation of additional Sn-OH and silanol groups away from the active site. The diffusion coefficients of water and methanol are functions of the solvent population in the pore. Here, our work provides insights into how active site speciation in Lewis acid zeolites affects solvent coordination, diffusion, and vibrational signature. This information is foundational for catalyst design and optimization of liquid-phase catalytic processes in zeolites. It also demonstrates the suitability of machine-learned interatomic potentials for modeling reactive systems, enabling sufficiently long trajectories for appropriate statistical averaging.

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

A first principles analysis of potential-dependent structural evolution of active sites in Fe-N-C catalysts

Fe-N-C (iron–nitrogen–carbon) electrocatalysts have emerged as potential alternatives to precious metal-based materials for the oxygen reduction reaction (ORR). However, the structure of these materials under electrochemical conditions is not well understood, and their poor stability in acidic environments poses a formidable challenge for successful adoption in commercial fuel cells. To provide molecular-level insights into these complex phenomena, we combine periodic density functional theory (DFT) calculations, exhaustive treatment of coadsorption effects for ORR reaction intermediates, including O and OH, and comprehensive analysis of solvation stabilization effects to construct voltage-dependent ab initio thermodynamic phase diagrams that describe the in situ structure of the active sites. These structures are further linked to activity and stability descriptors that can be compared with experimental parameters such as the half-wave potential for ORR and the onset potential for carbon corrosion and CO 2 evolution. The results indicate that pyridinic Fe sites at zigzag carbon edges, as well as other edge sites, exhibit high activity for ORR compared to sites in the bulk. However, edges neighboring the active sites are prone to instability via overoxidation and consequent site loss. The results suggest that it could be beneficial to synthesize Fe-N-C catalysts with small sizes and large perimeter edge lengths to enhance ORR activity, while voltage fluctuations should be limited during fuel cell operation to prevent carbon corrosion of overoxidized edges.

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

Discerning Influences from Enthalpy and Entropy at Aqueous Interfaces Involved in Biomass Conversions in Porous Catalysts

Project Summary: The goal in this proposal is to learn how solvent influences the enthalpies and entropies of catalytic species in zeolite pores. Specifically, the physical, chemical, and structural features of solvent, catalytic species, and zeolite pores that determine enthalpies and entropies of solvation in solution phase biomass conversions will be interrogated using multiscale simulations and machine learning. Multiscale simulations are based off of existing strategies and employ quantum mechanics and classical molecular dynamics, providing an excellent balance between chemical accuracy and computational expense. They are capable of calculating enthalpies and entropies of solvation separately and have been validated in prior work to achieve high accuracy compared to their parent methods.

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

Solvation of furfural at metal–water interfaces: Implications for aqueous phase hydrogenation reactions

Metal–water interfaces are central to understanding aqueous-phase heterogeneous catalytic processes. However, the explicit modeling of the interface is still challenging as it necessitates extensive sampling of the interfaces’ degrees of freedom. Herein, we use ab initio molecular dynamics (AIMD) simulations to study the adsorption of furfural, a platform biomass chemical on several catalytically relevant metal–water interfaces (Pt, Rh, Pd, Cu, and Au) at low coverages. We find that furfural adsorption is destabilized on all the metal–water interfaces compared to the metal–gas interfaces considered in this work. This destabilization is a result of the energetic penalty associated with the displacement of water molecules near the surface upon adsorption of furfural, further evidenced by a linear correlation between solvation energy and the change in surface water coverage. To predict solvation energies without the need for computationally expensive AIMD simulations, we demonstrate OH binding energy as a good descriptor to estimate the solvation energies of furfural. Using microkinetic modeling, we further explain the origin of the activity for furfural hydrogenation on intrinsically strong-binding metals under aqueous conditions, i.e., the endothermic solvation energies for furfural adsorption prevent surface poisoning. In conclusion, our work sheds light on the development of active aqueous-phase catalytic systems via rationally tuning the solvation energies of reaction intermediates.

aqueous-phase hydrogenation↗
DOE OSTI2025

Mechanistic Insights into the Demethylation of Lignin‐Derived Structures Using Protic Ionic Liquids: A Density Functional Theory Study

Lignin valorization is restricted by the stability of its methoxy groups, creating a critical need for efficient demethylation strategies. Here, in this study, density functional theory (DFT) is employed to dissect the mechanistic pathways of demethylation in lignin model compounds, guaiacol and syringol, using protic ionic liquids (PILs) that act as both solvent and catalyst. Conductor‐like Screening Model for Real Solvents (COSMO‐RS) analysis identifies monoethanolammonium acetate ([MEOA][Ace]) as the most promising medium, attributed to its strong hydrogen bonding network and solvation ability. By integrating implicit and explicit solvation models, it is revealed that an acid‐catalyzed hydrolytic mechanism governs demethylation, with PILs stabilizing crucial transition states and intermediates. Complementary electronic structure evaluations, including highest occupied molecular orbital–lowest unoccupied molecular orbital (HOMO‐LUMO) gap analysis, charge distribution, and electrostatic potential mapping, demonstrate how PILs lower energetic barriers and enhance reactivity. To simulate realistic environments, this study is extended to lignin dimer complexes with varying water content, uncovering how water‐bridged solvation reverses demethylation preference from guaiacyl to syringyl units. This mechanistic shift aligns with experimental observations showing faster S‐unit reactivity in hydrated systems. Together, these findings provide atomic‐level insight into lignin demethylation dynamics and highlight how tuning acid concentration in PILs can accelerate kinetics, enabling a rational pathway toward next‐generation biomass conversion technologies.

COSMO-RS↗
DOE OSTI2025

Probing Condensed-Phase Structure and Dynamics in Hierarchical Zeolites and Nanosheets for Catalytic Upgradation of Biomass (Final Report)

Understanding complex reaction pathways in systems governed by multi-scale collective interactions across time and length scales remains a central scientific challenge. This project was guided by the hypothesis that the interplay among oligomers, solvents, and active sites can be tuned by a suitable choice of solvation environment and pore architecture in solid-acid catalysts to direct chemical transformations relevant to biomass conversion. Zeolites and zeolite nanosheets were used as model platforms, allowing for the interaction of macromolecules with the surface of the zeolite nanosheets and with smaller pores that host catalytically active sites. To investigate these coupled phenomena, we employ a multi-scale computational framework that integrates molecular-level descriptions with advanced sampling approaches to capture key physical and chemical interactions. Our work through this project improved fundamental understanding of how reactants and solid-acid catalysts interact in solvent-rich environments, thereby enabling the rational design of catalytic systems that upgrade biomass with enhanced selectivity and energy efficiency. In addition, the project developed advanced sampling methodologies critical for disentangling complex, reactive processes in multi-component catalytic environments.

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

Interactive multiscale modeling to bridge atomic properties and electrochemical performance in Li-CO 2 battery design

Li-CO 2 batteries are promising energy storage systems due to their high theoretical energy density and CO 2 fixation capability, relying on reversible Li 2 CO 3 /C formation during discharge/charge cycles. Here, we present a multiscale modeling framework integrating Density Functional Theory (DFT), Ab-Initio Molecular Dynamics (AIMD), classical Molecular Dynamics (MD), and Finite Element Analysis (FEA) to investigate atomic and cell-level properties. The considered Li-CO 2 battery consists of a lithium metal anode, an ionic liquid electrolyte, and a carbon cloth cathode with Sb 0.67 Bi 1.33 Te 3 catalyst. DFT and AIMD determined the electrical conductivities of Sb 0.67 Bi 1.33 Te 3 and Li 2 CO 3 using the Kubo–Greenwood formalism and studied the CO 2 reduction mechanism on the cathode catalyst. MD simulations calculated the CO 2 diffusion coefficient, Li + transference number, ionic conductivity, and Li + solvation structure. The FEA model, parameterized with atomistic simulation data, reproduced the available experimental voltage–capacity profile at 1 mA/cm 2 and revealed spatio-temporal variations in Li 2 CO 3 /C deposition, porosity, and CO 2 concentration dependence on discharge rates in the cathode. Accordingly, Li 2 CO 3 can form large and thin film deposits, leading to dispersed and local porosity changes at 0.1 mA/cm 2 and 1 mA/cm 2 , respectively. The capacity decreases exponentially from 81,570 mAh/g at 0.1 mA/cm 2 to 6200 mAh/g at 1 mA/cm 2 , due to pore clogging from excessive discharge product deposition that limits CO 2 transport to the cathode interior. Therefore, the performance of Li-CO 2 batteries can be improved by enhancing CO 2 transport, regulating Li 2 CO 3 deposition, and optimizing cathode architecture.

Battery performance↗
DOE OSTI2021

Solvent Effects on Catalytic Reactions and Related Phenomena at Liquid-Solid Interfaces

Catalytic reactions involve the direct interaction of reactants, intermediates and products with the catalyst surface. We not only need to control the atomic structure and electronic properties of the active site, but also explore the multiple molecular interactions that occur beyond the active site; they play an essential role in altering the binding and reactivity of surface species. In liquid-phase catalysis, solvents provide additional degrees of freedom in the design of the catalytic process for desirable activity and selectivity. Here, the multi-faceted effects of solvents have a profound impact on the catalyst performance by restricting the mass transfer to the site, tuning the chemical potential of the surface species, competing for active sites, stabilizing the initial and transition states, and causing mechanistic changes by participating in the kinetically relevant elementary steps. This review addresses the different aspects of solvent effects, using a few prototype solid-liquid interfaces to illustrate these fundamental features. Recent experimental and computational studies that provide new insight at the molecular level are examined. Solvent structures in the proximity of the catalyst surface are discussed along with their influence in molecular binding and reaction at the solid-liquid interfaces. Furthermore, opportunities to alter such a solid-liquid interaction by tuning the wettability of the catalyst surfaces are explored.

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

Protonation Dynamics of Confined Ethanol–Water Mixtures in H-ZSM-5 from Machine Learning-Driven Metadynamics

Zeolites are indispensable heterogeneous catalysts in industrial chemical processes, valued for their strong Brønsted acidity, well-defined microporous frameworks, and tunable pore structures. Their catalytic activity arises primarily from Brønsted acid sites (BAS), typically present as bridging hydroxyl groups (Si–OH–Al). Under aqueous reaction conditions, these protons interact dynamically with water and alcohol molecules, leading to complex solvation and protonation behavior within confined pores. In this study, we investigate the protonation equilibrium occurring between ethanol and water at the BAS of acidic zeolites under varying hydration levels, i.e., C2H5OH–(H2O)n, n=1–4. Local structure was analyzed through an adaptive-learning global optimization algorithm, while enhanced sampling molecular dynamics simulations with Well-Tempered Metadynamics (WMetaD) and machine learning interatomic potentials (MLPs) provide free-energy surfaces (FES) at variable hydration levels. The results reveal a strong dependence of proton localization on the degree of hydration. At low hydration (1 water molecule), the proton resides predominantly on ethanol; with 2 water molecules, it shifts toward water, and at higher hydration (3 or more water molecules), it becomes extensively delocalized over the water cluster. These findings underscore the critical role of solvation in modulating acid site behavior and suggest that a minimum of three water molecules is necessary to fully stabilize the proton on water within the zeolite framework. This solvation threshold has significant implications for catalytic processes, particularly in biomass conversion reactions where alcohol protonation is a key step in dehydration mechanisms.

machine learning↗