Atomistic Insights into Aqueous Phase Carbonyl Hydrogenation over Ru Catalysts from a Hybrid QM/MM/ML Potential

02 December 2025, Version 1
This content is an early or alternative research output and has not been peer-reviewed by Cambridge University Press at the time of posting.

Abstract

The aqueous-phase hydrogenation of methyl ethyl ketone (2-butanone, MEK) to 2-butanol over Ru(0001) was investigated as a model system to elucidate condensed-phase effects in the catalytic hydrogenation of carbonyl-containing feedstocks. A hybrid quantum mechanical/molecular mechanical/machine learning free-energy perturbation (QM/MM/ML-FEP) framework was employed to quantify solvent effects on the relative stabilization of surface intermediates and transition states. Microkinetic modeling reveals that the Alkoxy pathway is favored over the Hydroxy pathway across all environments—dry gas, co-fed water vapor, and liquid water. Water plays a dual mechanistic role: it directly promotes hydrogenation of the carbonyl oxygen via proton shuttling/hydrogen-atom transfer from adsorbed H2O/OH species, and it acts as a solvent that modulates the adsorption strength of surface intermediates and transition states. In the co-fed water vapor system, the second hydrogenation step proceeds through H-mediated (63%), H₂O-mediated (31%), and OH-mediated (5%) pathways, whereas in liquid water, the reaction occurs exclusively through the H₂O-mediated pathway, enhancing the overall turnover frequency by an order of magnitude and shifting the rate controlling step from the second to the first hydrogenation of the carbonyl carbon. These findings demonstrate that the enhanced activity of aqueous-phase Ru-catalyzed carbonyl hydrogenation arises from both the direct participation of water in the mechanism and its solvent-induced stabilization and destabilization of key surface states.

Keywords

Heterogeneous Catalysis
Solid-Liquid Interface
Quantum Mechanical/ Molecular Mechanical Model
QM/MM
eSMS
Machine Learning (ML)
RuNNer

Supplementary materials

Title
Description
Actions
Title
Supporting Information for: Atomistic Insights into Aqueous Phase Carbonyl Hydrogenation over Ru Catalysts from a Hybrid QM/MM/ML Potential
Description
Lennard-Jones parameters used in MD simulations, average rotational correlation times, charges on the QM clusters, values for the parameter δ in the removal of LJ interactions, optimized geometry of clusters, lateral interaction parameters employed in the microkinetic model, average adsorption energy as a function of H coverage on the surface, average adsorption energy as a function of H2O coverage on the surface, linear and quadratic lateral interaction model for H-H and H2O-H2O, average free-energy profiles for aqueous-phase effects on the binding strength of transition states, height distribution function of the oxygen atoms of the aqueous phase as a function of height (perpendicular distance) from the metal surface.
Actions
Title
Adsorbate and TS coordinates
Description
TXT file includes all fractional coordinates in VASP POSCAR format and XYZ coordinate of the QM clusters format of all DFT-optimized ground states and transition states utilized in this work.
Actions

Comments

Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting and Discussion Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.