Search NASA⌕ Search

Engineering topics

Yuk, Simuck F.

Publications and source records attributed to Yuk, Simuck F..

At least 19 records
All sourcesNASA NTRSDOE OSTI
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↗
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 OSTI2023

Dynamically Formed Active Sites on Liquid Boron Oxide for Selective Oxidative Dehydrogenation of Propane

Boron-based catalysts have been shown to be both active and selective for driving the oxidative dehydrogenation of propane (ODHP) without the use of precious metals. This reaction occurs at temperatures that melt the oxide catalyst which challenges our ability to identify the liquid structures of the boron oxide phase under reaction conditions, hindering the understanding of its active sites and reaction mechanism. By combining ab initio molecular dynamics simulation, in-situ Raman characterization, and microkinetic modeling, we propose that the di-coordinated boron sites (BO2) in liquid boron oxide are the active species for O2 activation under reaction conditions. The formed peroxy-like species (>B-O-O-B<) can be viewed as a moderate oxidant for ODHP. The dynamical >B-O* dangling bond originated from >B-O-O-B< site as well as the liquid B2O3 structure itself, plays a critical role in the abstraction of H atoms from propane (C3H7 radical formation). Microkinetic modeling reveals C3H7 radical formation to be the main rate controlling step (~75% degree of rate control) with the dehydration of boron hydroxyls (B-OHs) to recover the di-coordinated boron active sites controlling the remainder of the rate (~25% degree of rate control). Moreover, the activation barriers are found to strongly depend upon the surface B-OH concentration. These findings provide significant insights into the active site and reaction mechanisms on boron-based catalysts for ODHP and underlie the importance of understanding the liquid nature of the catalyst to account for the catalytic activity.

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

Assessing entropy for catalytic processes at complex reactive interfaces

When chemical reactions are accelerated by a catalyst, entropy differences between reactants and their transient intermediates can be the driving force behind the promotion or inhibition of desired and parasitic chemical pathways. Understanding and controlling catalytic processes therefore requires both a fundamental and practicable understanding of entropy in addition to enthalpy. In unstructured media such as the vapor phase equilibrated with sparsely covered surfaces, entropy can be adequately accounted for by well-established approaches based on translational, rotational, and harmonic vibrational partition functions. However, these approximations become inadequate in more complex condensed phase environments, e.g., solid liquid interfaces of confined reaction spaces. In this chapter, we provide an overview of the state-of-art in the computational quantification of entropy and its known ramifications on catalysis. The fundamental roles of thermodynamics and kinetics in catalysis are covered in enough detail to appreciate and contextualize the computational methods employed to compute chemically accurate estimates of entropy. These methods are discussed in appropriate detail and range from the ubiquitous harmonic oscillator approximation where entropy unrelated to high frequency oscillations is typically underestimated, to enhanced free energy sampling with molecular dynamics where the desired accuracy must be weighed against the associated computational cost of obtaining it. The rising importance of machine learning and artificial intelligence in accelerating methodological progress in this field is touched upon, as well. Finally, applications, successes, and pitfalls of using these methods are provided to showcase past and present accomplishments while clarifying where improvements in both understanding and methodology are still needed.

Kollias, Loukas↗
DOE OSTI2021

Impact of functional groups on the electrocatalytic hydrogenation of aromatic carbonyls to alcohols

We report that electrocatalytic hydrogenation (ECH) of biomass-derived feedstocks has a critical dependence on the molecular structure of the organic and its adsorption on the electrode surface. In this study, we investigated the role of functional groups in the adsorption of the organic molecule on the charged Pd (111) surface and its subsequent effect on organic reduction in electrochemical hydrogenation of organic molecules. With three aromatic carbonyls of benzaldehyde (BZD), acetophenone (ACE), and vanillin (VAN), we rationalize molecular-scale adsorption and interfacial charge transfer processes by employing density-functional-theory based ab initio molecular dynamics simulations. We observe that the functional group and electrode charge strongly affect the proximity of organic molecule to the Pd (111) surface, where distances of aromatic ring and carbonyl group of the organic on the electrode change distinctively with functional groups and charge state of electrode, which strongly impact reduction of organics on the surface. Calculations of differential electron density show the strongest reduction with benzaldehyde via interfacial electron transfer from the charged Pd surface. We also observe that the interaction between the functional groups and solvent (VAN > BZD > ACE) significantly influence the organic interaction with the charged electrode (BZD > VAN > ACE), resulting in the net interaction energy between the organic and the electrode in the order of BZD > ACE > VAN. Experimental measurement of ECH rate also show the same trend of the net interaction energy. These results demonstrate the significance of solvent effect on the reducibility of organic molecules on electrodes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ab initio molecular dynamics with enhanced sampling in heterogeneous catalysis

Ab initio molecular dynamics simulations combined with enhanced sampling techniques are becoming widespread methods to investigate chemical phenomena in catalytic systems. These techniques automatically include finite temperature effects, anharmonicity, and collective dynamics in their robust description of enthalpic and entropic contributions, which can have significant impact on reaction free energy landscapes. This contrasts with standard ab initio static approaches that are based on assessing reaction free energies from various coarse-grained descriptions of the reaction potential energy surface. Enhanced sampling ab initio molecular dynamics opens the way to first principles simulations of systems of increasing complexity like solid/liquid catalytic interfaces. Here, we aim at guiding the reader through the basis of these techniques, summarizing their fundamental theoretical and practical aspects, and reviewing the relevant literature in the field. After a brief introduction to the problem, we will illustrate the advantage of using molecular simulations to include finite temperature effects, examine the most common ab initio techniques currently in use, describe their application to solid state heterogeneous catalysts, and finally critically review the most popular enhanced sampling techniques used in computational catalysis.

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

Environment of Metal–O–Fe Bonds Enabling High Activity in CO 2 Reduction on Single Metal Atoms and on Supported Nanoparticles

Single-atom catalysts are often reported to have catalytic properties that surpass those of nanoparticles, while a direct comparison of sites common and different for both is lacking. In this work, we show that single atoms of the Pt-group embedded into the surface of Fe 3 O 4 have a greatly enhanced interaction strength with CO 2 compared with Fe 3 O 4 surface. The strong CO 2 adsorption on single Rh atoms and corresponding low activation energies lead to two-orders-of-magnitude higher conversion rates of CO 2 compared to Rh nanoparticles. This high activity of single atoms stems from the partially oxidic state imposed by their coordination to the support. Fe 3 O 4 -supported Rh nanoparticles follow the behavior of single atoms for CO 2 interaction and reduction, which is attributed to the dominating role of partially oxidic sites at the Fe 3 O 4 -Rh interface. Thus, we show a likely common catalytic chemistry for two kinds of materials thought to be different, and we show that single atoms of Pt-group metals on Fe 3 O 4 are an especially successful material for catalyzed reactions that depend primarily upon sites with the metal-O-Fe environment.

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

Mechanistic understanding of catalytic conversion of ethanol to 1-butene over 2D-pillared MFI zeolite

Ethanol is an important C2 platform molecule for producing value-added chemicals and distillate hydrocarbon fuels (e.g., jet and diesel). Among these, catalytic upgrading of ethanol to butenes can generate valuable commodity chemicals (e.g., 1-butene) and provide C4 olefin intermediates that can be further upgraded to jet/diesel fuels. 2D zeolites offer hierarchical mesoporous structures, leading to improved mass transport and reduced diffusion length, which can help to address the coking challenges faced by ethanol conversion to hydrocarbons over 3D zeolites. In this study, we investigate the acid-catalyzed conversion of ethanol to 1-butene over the Brønsted acid sites (BAS) in 2D-pillared MFI zeolite (2D-PMFI) using ab initio molecular dynamics (AIMD) simulations, in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS), and calorimetric measurements. A detailed thermodynamic analysis, using quasi-harmonic approximation (QHA), on the Gibbs free energy pathway of ethanol conversion shows that the consideration of entropy is critical to accurately capture the detailed thermodynamic profiles. The formation of bulky ethoxy-ethene complex is found to be a potential rate-determining step (RDS), proceeding via stepwise mechanism. The reactivity of 2D-PMFI can be further tuned by manipulating RDS through a careful control of number of BAS and operating temperatures. The accuracy of computational model is confirmed through the comparison of calculated vibrational density of states (VDOS) against the experimental DRIFTS measurements. Overall, our study provides mechanistic insights into ethanol upgrading over the 2D-PMFI and shows the importance of evaluating entropic effects in such confined system.

Ethanol conversion, 2D-pillared zeolite, MFI, Anha↗
DOE OSTI2020

Perspectives on van der Waals Density Functionals: The Case of TiS 2

The van der Waals interaction is of foundational importance for a wide variety of physical systems. In particular, van der Waals forces lie at the heart of potential device technologies that may be realized from the functional organization of layered two-dimensional (2D) nanomaterials. Furthermore, for intermediate to large-scale applications modeling, van der Waals density functionals have become the de facto choice for first-principles calculations. In particular, the vdW-DF family of functionals have provided a systematic approach to this theoretically challenging problem. While much progress has been made, there remains room for improvement in the microscopic description of vdW forces from these density functionals. In this work, we compute benchmark results for the binding energy and the electronic density response to binding in TiS2 via accurate diffusion quantum Monte Carlo calculations. We compare these benchmark data to results obtained from local, semilocal, and van der Waals functionals. In particular, we gauge the quality of the original vdW-DF/vdW-DF2 functionals, as well as updated variants such as vdW-DF-C09, vdW-DF-optB88, vdW-DF-optB86b, and vdW-DF2-B86R. We find a close relationship between the accuracy of predicted interlayer separation distances and binding energies for TiS 2 , with the vdW-DF-optB88 functional performing very well in terms of both quantities. In general, the more recently developed functionals are systematic improvements over older ones. However, when considering the response of the electron density to binding, we find that local-density approximation (LDA) and PBEsol generally outperform the vdW-DF functionals in describing the interlayer charge accumulation with vdW-DF-C09 variants performing the best overall.

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

Hydrogen Bonding Enhances the Electrochemical Hydrogenation of Benzaldehyde in the Aqueous Phase

Abstract The hydrogenation of benzaldehyde to benzyl alcohol on carbon‐supported metals in water, enabled by an external potential, is markedly promoted by polarization of the functional groups. The presence of polar co‐adsorbates, such as substituted phenols, enhances the hydrogenation rate of the aldehyde by two effects, that is, polarizing the carbonyl group and increasing the probability of forming a transition state for H addition. These two effects enable a hydrogenation route, in which phenol acts as a conduit for proton addition, with a higher rate than the direct proton transfer from hydronium ions. The fast hydrogenation enabled by the presence of phenol and applied potential overcompensates for the decrease in coverage of benzaldehyde caused by competitive adsorption. A higher acid strength of the co‐adsorbate increases the intensity of interactions and the rates of selective carbonyl reduction.

Sanyal, Udishnu↗
DOE OSTI2020

First-principle investigation on catalytic hydrogenation of benzaldehyde over Pt-group metals

We report that understanding the hydrogenation of organic compounds in the aqueous phase has always been fundamentally important for improving carbon neutral pathways to fuels and value-added chemicals. In this study, we investigated both thermodynamic and kinetic profiles of benzaldehyde hydrogenation over the Pd(111) and Pt(111) metal surfaces using density functional theory (DFT) and ab initio molecular dynamic (AIMD) simulations. The adsorption of H 2 shows the mixed preference of H adsorption sites on the Pt(111), while the fcc adsorption site is dominant for H on the Pd(111). When benzaldehyde is added to the systems, we observe a strong reduction of benzaldehyde on charged Pd (111) surface compared with that on neutral surface. In contrast, charged state of the Pt(111) surface does not change their interaction. Subsequent hydrogenation reaction of benzaldehyde over Pd(111), proceeding via Langmuir-Hinshelwood mechanism, is affected by two major factors: the presence of H 2 O solvent and surface charge. The presence of H 2 O solvent greatly reduces the activation energy of C-H and O-H bond formation during the hydrogenation process. Furthermore, the hydrogenation step via C-H bond formation is preferred thermodynamically and kinetically over O-H bond formation during thermocatalytic hydrogenation, while the opposite trend holds true during electrocatalytic hydrogenation.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Effect of Collective Dynamics and Anharmonicity on Entropy in Heterogenous Catalysis: Building the Case for Advanced Molecular Simulations

In this paper, we present a perspective on the computational determination of entropy and its effects and consequences on heterogeneous catalysis. Special attention is paid to the role of anharmonicity and the deviations from the standard harmonic oscillator approximations which can fail to provide a reliable assessment of entropy. To address these challenges, advanced methodologies are needed that can explicitly account for these thermodynamic drivers through the appropriate statistical sampling of reactive free energy surfaces. We discuss where anharmonicity should be expected, where it has been observed from a theoretical perspective, and the methods currently employed to address it. We concentrate on three types of systems where we have observed major, nonnegligible anharmonic effects: (1) supported nanoparticles, where the migration of metal atoms, complexes, and entire clusters exhibit anharmonic behavior in their dynamic motion; (2) porous solids where confinement effects distort potential energy surfaces and hinder molecular motions resulting in large entropic terms; and (3) solid/liquid interfaces, where interactions between solvent molecules and adsorbed species can result in large solvent organization free energy and unique reactivity.

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

Influence of Ag metal dispersion on the thermal conversion of ethanol to butadiene over Ag-ZrO 2 /SiO 2 catalysts

Atomistic scale models were developed and coupled with experimental investigation to deliver a functional understanding of catalytic activity and selectivity in the conversion of ethanol to 1,3-butadiene over Ag/ZrO 2 /SiO 2 . A detailed evaluation of the structural and electronic properties of the resultant catalyst models led to the identification of critical active sites of the catalyst. More importantly, the extent of Ag dispersion on the SiO 2 support and relative proximity to ZrO 2 were found to vary with the oxidation state of Ag and local coordination environment (Ag-O SiO2 ) allowing for critical control of ethanol conversion towards butadiene or ethylene. Simulations revealed that less dispersed or clustered Ag contain predominantly Ag 0 charge state and promote conversion of ethanol to ethylene. The well-dispersed Ag/ZrO 2 /SiO 2 catalyst instead contain a larger fraction of cationic Ag δ+ , and predominantly promote ethanol dehydrogenation and subsequent production of butadiene. The theoretical insights drawn were validated and confirmed experimentally using TEM, XRD and reactivity measurements demonstrating the effect of Ag dispersion on the selectivity of ethanol conversion.

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

Selective acetylene hydrogenation over single metal atoms supported on Fe 3 O 4 (001): A first-principle study

Supported single-atom catalysts (SACs) have gained increasing attention for improved catalytic activity and selectivity for industrially-relevant reactions. In this study, we explore the hydrogenation of acetylene over single Pt, Ru, Rh, Pd, and Ir atoms supported on the Fe 3 O 4 (001) surface using density functional theory (DFT) calculations. The thermodynamic profile of H diffusion is significantly changed with the type of single metal atoms used and suggests that H spillover from the single atom dopant to the Fe 3 O 4 (001) surface is favored and will likely lead to high H coverages in the functioning catalyst. Correspondingly, as the surface H coverage increases, the important desorption step of ethylene becomes energetically competitive against the detrimental hydrogenation steps of ethylene to ethane. Overall, we show herein that the selectivity of SACs towards ethylene production can be tuned by considering the proper type of metal and controlling the redox state of the support.

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

Single Atom Catalysis: An Analogy between Heterogeneous and Homogeneous Catalysts

In recent years, enormous efforts have been invested to improve the atom-efficiency of metal catalysts by decrease the size of metal particles to the limit of isolated atoms in heterogenous catalysis, termed single atom catalysts (SACs). In addition to being atom-efficient, these SACs are increasingly gaining attention due to the unique reactivity and selectivity they have exhibited in industrially relevant reactions, including both oxidation and hydrogenation reactions. Clear parallels can be drawn between heterogeneous SACs and their analogous homogeneous catalysts. In this chapter, we highlight SACs that demonstrate this relationship and report on relevant, recent developments in the field of single atom catalysis.

Yuk, Simuck F.↗
DOE OSTI2019

Low-Temperature Oxidation of Methanol to Formaldehyde on a Model Single-Atom Catalyst: Pd Atoms on Fe 3 O 4 (001)

Single atom catalysis has been a topic of increasing interest due to the potential for improved selectivity, reactivity, and catalyst cost. However, single atom catalysts are still difficult to characterize under realistic reaction conditions, leading to controversy regarding the capabilities of single atoms and a need for model studies. Herein, we examine the reaction of methanol on single Pd atoms supported on Fe 3 O 4 (001) under ultra-high vacuum conditions. On clean Fe 3 O 4 (001) a small fraction of methanol is converted to formaldehyde through a methoxy intermediate at 516 K. The addition of single Pd atoms lowers the barrier to C-H bond cleavage by a factor of two, resulting in formaldehyde desorption by 290 K. Furthermore, Pd atoms sinter by 300 K in the presence of methanol, and Pd clusters do not exhibit the same chemistry. Single atom sites significantly lower the barrier to the oxidation of methanol, although their stability remains an issue.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗