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.
Supplementary materials
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



![Author ORCID: We display the ORCID iD icon alongside authors names on our website to acknowledge that the ORCiD has been authenticated when entered by the user. To view the users ORCiD record click the icon. [opens in a new tab]](https://www.cambridge.org/engage/assets/public/coe/logo/orcid.png)