Hybrid Quantum Mechanical, Molecular Mechanical, and Machine Learning Potential for Computing Aqueous-Phase Adsorption Free Energies on Metal Surfaces

Performing reliable computer simulations of elementary processes occurring at metal-water interfaces is pivotal for novel catalyst design in sustainable energy applications. Computational catalyst design hinges on the ability to reliably and efficiently compute the potential energy surface (PES) of the system. Due to the large system sizes needed for studying processes at liquid water-metal interfaces, these systems can currently not be described using density functional theory (DFT). In this work, we used a hybrid quantum mechanical, molecular mechanical, and machine learning potential for studying the adsorption behavior of phenol, atomic hydrogen, 2-butanol, and 2-butanone on the (0001) facet of Ru under reducing conditions when Ru is not oxidized. Specifically, we describe the adsorbate and the surrounding metal atoms at the DFT level of theory. Here, we also considered the electrostatic field effect of the water molecules on adsorbate-metal interactions. Next, for the water-water and water-adsorbate interactions, we used established classical force fields. Finally, for the water-Ru surface interaction, for which no reliable force fields have been published, we used Behler-Parrinello high-dimensional neural network potentials (HDNNPs). Employing this setup, we used our explicit solvation for metal surface (eSMS) approach to compute the aqueous-phase effect on the low-coverage adsorption of selected molecules and atoms on the (0001) facet of Ru. In agreement with previous experimental and computational studies of oxygenated molecules over transition metal facets, we found that liquid water destabilizes the tested adsorbates on Ru(0001). Interestingly, our findings indicate that adsorbates on Ru are less affected by the presence of an aqueous phase than on other transition metals (e.g., Pt), highlighting the necessity of experimental investigations of Ru-based catalytic systems in liquid water.

Liquid-Phase Effects on Adsorption Processes in Heterogeneous Catalysis

Aqueous solvation free energies of adsorption have recently been measured for phenol adsorption on Pt(111). Endergonic solvent effects of ∼1 eV suggest solvents dramatically influence a metal catalyst's activity with significant implications for the catalyst design. However, measurements are indirect and involve adsorption isotherm models, which potentially reduces the reliability of the extracted energy values. Computational, implicit solvation models predict exergonic solvation effects for phenol adsorption, failing to agree with measurements even qualitatively. In this study, an explicit, hybrid quantum mechanical/molecular mechanical approach for computing solvation free energies of adsorption is developed, solvation free energies of phenol adsorption are computed, and experimental data for solvation free energies of phenol adsorption are reanalyzed using multiple adsorption isotherm models. Explicit solvation calculations predict an endergonic solvation free energy for phenol adsorption that agrees well with measurements to within the experimental and force field uncertainties. Computed adsorption free energies of solvation of carbon monoxide, ethylene glycol, benzene, and phenol over the (111) facet of Pt and Cu suggest that liquid water destabilizes all adsorbed species, with the largest impact on the largest adsorbates.

Direct Observation of the Dynamics of Ylide Solvation by Hydrogen-bond Donors Using Time-Resolved Infrared Spectroscopy

The photoexcitation of α-diazocarbonyl compounds produces singlet carbene intermediates that react with nucleophilic solvent molecules to form ylides. The zwitterionic nature of these newly formed ylides induces rapid changes in their interactions with the surrounding solvent. Here, ultrafast time-resolved infrared absorption spectroscopy is used to study the ylide-forming reactions of singlet carbene intermediates from the 270 nm photoexcitation of ethyl diazoacetate in various solvents and the changes in the subsequent ylide-solvent interactions. The results provide direct spectroscopic observation of the competition between ylide formation and C-H insertion in reactions of the singlet carbene with nucleophilic solvent molecules. We further report the specific solvation dynamics of the tetrahydrofuran (THF)-derived ylide (with a characteristic IR absorption band at 1636 cm-1) by various hydrogen-bond donors and the coordination by lithium cations. Hydrogen-bonded ylide bands shift to a lower wavenumber by -19 cm-1 for interactions with ethanol, -14 cm-1 for chloroform, -10 cm-1 for dichloromethane, -9 cm-1 for acetonitrile or cyclohexane, and -16 cm-1 for Li+ coordination, allowing the time evolution of the ylide-solvent interactions to be tracked. The hydrogen-bonded ylide bands grow with rate coefficients that are close to the diffusional limit. We further characterize the specific interactions of ethanol with the THF-derived ylide using quantum chemical (MP2) calculations and DFT-based atom-centered density matrix propagation trajectories, which show preferential coordination to the α-carbonyl group. This coordination alters the hybridization character of the ylidic carbon atom, with the greatest change toward sp2 character found for lithium-ion coordination.

Low temperature catalytic hydrodeoxygenation of lignin-derived phenols to cyclohexanols over the Ru/SBA-15 catalyst

Cyclohexanol and its derivatives are widely used as chemical intermediates and fuel additives. Herein, Ru/SBA-15 catalysts were prepared via impregnation, and used for the production of cyclohexanols from lignin-derived phenols. The catalyst samples were characterized by XRD, XPS, TEM, etc., where the Ru0 species was speculated as the active phase. 5 wt% Ru/SBA-15 with small Ru particle size (4.99 nm) and high Ru dispersion (27.05%) exhibited an excellent hydrogenation activity. A high cyclohexanol yield of >99.9% was achieved at 20 °C for 5 h in an aqueous phase, and the catalyst indicated stable activity and selectivity after five runs. Crucially, Ru/SBA-15 exhibited a zero-order reaction rate with an apparent activation energy (Ea) as low as 10.88 kJ mol-1 and a TON of 172.84 at 80 °C. Simultaneously, demethoxylation activity was also observed in the hydrodeoxygenation (HDO) of G- and S-type monophenols, and a high yield of 37.4% of cyclohexanol was obtained at 80 °C and 4 h when using eugenol as substrate.

Hydrodeoxygenation of guaiacol over orthorhombic molybdenum carbide: a DFT and microkinetic study

The hydrodeoxygenation of guaiacol is modelled over a (100) β-Mo2C surface using density functional theory and microkinetic simulations.

The Effect of Gold Nanoparticles on the Catalytic Activity of NiTiO3 for Hydrodeoxygenation of Guaiacol

Guaiacol is a typical model compound used to investigate and understand the hydrodeoxygenation behaviour of bio-oils, which is critical to their application as an alternative to fossil resources. While extensive research has been carried out on developing catalysts for guaiacol hydrodeoxygenation, the true active sites in these catalysts are often illusive. This study investigated the effect of Au-loading on the catalytic activity of NiTiO3 for the hydrodeoxygenation of guaiacol. It showed that metallic Ni formed by the partial reduction in NiTiO3 was responsible for its catalytic activity. Au-loading in NiTiO3 effectively reduces the temperature required for the NiTiO3 reduction from 400 °C to 300 °C. Consequently, at an Au-loading of 0.86 wt%, the 0.86 Au/NiTiO3-300 °C catalyst was found to deliver a guaiacol conversion of ~32%, more than 6 times higher than that of the pure NiTiO3-300 °C catalyst.

Aqueous-phase effects on ethanol decomposition over Ru-based catalysts

Liquid water decelerates ethanol reforming over Ru(0001) but increases the H2 selectivity due to accelerated WGS and suppressed methanation.

Dependency of solvation effects on metal identity in surface reactions

Solvent interactions with adsorbed moieties involved in surface reactions are often believed to be similar for different metal surfaces. However, solvents alter the electronic structures of surface atoms, which in turn affects their interaction with adsorbed moieties. To reveal the importance of metal identity on aqueous solvent effects in heterogeneous catalysis, we studied solvent effects on the activation free energies of the O-H and C-H bond cleavages of ethylene glycol over the (111) facet of six transition metals (Ni, Pd, Pt, Cu, Ag, Au) using an explicit solvation approach based on a hybrid quantum mechanical/molecular mechanical (QM/MM) description of the potential energy surface. A significant metal dependence on aqueous solvation effects was observed that suggests solvation effects must be studied in detail for every reaction system. The main reason for this dependence could be traced back to a different amount of charge-transfer between the adsorbed moieties and metals in the reactant and transition states for the different metal surfaces.