Liquid-Phase Effects on Adsorption Processes in Heterogeneous Catalysis

Mechanical Reversal in the Catalytic Capability of Monolayer Transition Metal Dichalcogenides for Hydrogen Evolution Reaction: An Explicit First-Principles Study

Pristine transition metal dichalcogenide (TMD) monolayers are generally regarded as exhibiting low chemical reactivity due to their inert surfaces. Our extensive first-principles calculations, which incorporate an explicit solvation model, reveal that the catalytic performance of pristine TMD MX2 (where M = Mo or W, and X = S, Se or Te) monolayers for hydrogen evolution reaction can be significantly altered and enhanced through mechanically bending deformation. For a WTe2 monolayer, its hydrogen adsorption Gibbs free energy decreases to 0.004 eV under a bending curvature of 0.15 Å-1. The notable reversal in the catalytic capability of curved TMD monolayers can be primarily ascribed to the interplay between elastic energy stimuli and hydrogen adsorption energy barrier, alongside charge transfer to metal atoms facilitated by the weakening of M-X bonds and the exposure of metal atoms. A theoretical model has been established to elucidate the relationship among hydrogen adsorption Gibbs free energy, bending elastic energy, and adsorption energy barrier.

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.

Selective Adsorbent Design with Multifunctional Surfaces: Innovating Solutions for Heterogeneous Catalysis in Water

Heterogeneous catalytic systems with water as the solvent often have the disadvantage of cross-contamination, while concerns about the purification and workup of the aqueous phase after reactions are rare in the lab or industry. In this context, designing and developing the functional selective solid adsorbent and revealing the adsorption mechanism can provide a new strategy and guidelines for constructing supported heterogeneous catalysts to address these issues. Herein, we report the stable composite adsorbent (Fe/ATP@PPy: magnetic Fe3O4/attapulgite with the polypyrrole shell) that features an integrated multifunctional surface, which can effectively tune the selective adsorption processes for Cu2+, Co2+, and Ni2+ ions and nitrobenzene via the cooperative chemisorption/physisorption in an aqueous system. The adsorption experiments showed that Fe/ATP@PPy displayed significantly higher adsorption selectivity for Ni2+ than Cu2+ and Co2+ ions, especially which exhibited an approximate 100.00% removal for both Ni2+ ions and nitrobenzene in the mixture system with a low concentration. Furthermore, combined tracking adsorption of Ni2+ ions and X-ray photoelectron spectroscopy characterization confirmed that the effective adsorption occurs via ion transfer coordination; the pathway was further validated at the molecular level through theoretical modeling. In addition, the selective adsorption mechanism was proposed based on the adsorption experiment, characterization, and the corresponding density functional theory calculation.

Molecularly imprinted solid-phase extraction combined with non-ionic hydrophobic deep eutectic solvents dispersed liquid-liquid microextraction for efficient enrichment and determination of the estrogens in serum samples

Hormonal drugs in biological samples are usually in low concentration and highly intrusive. It is of great significance to enhance the sensitivity and specificity of the detection process of hormone drugs in biological samples by utilizing appropriate sample pretreatment methods for the detection of hormone drugs. In this study, a sample pretreatment method was developed to effectively enrich estrogens in serum samples by combining molecularly imprinted solid-phase extraction, which has high specificity, and non-ionic hydrophobic deep eutectic solvent-dispersive liquid-liquid microextraction, which has a high enrichment ability. The theoretical basis for the effective enrichment of estrogens by non-ionic hydrophobic deep eutectic solvent was also computed by simulation. The results showed that the combination of molecularly imprinted solid-phase extraction and deep eutectic solvent-dispersive liquid-liquid microextraction could improve the sensitivity of HPLC by 33∼125 folds, and at the same time effectively reduce the interference. In addition, the non-ionic hydrophobic deep eutectic solvent has a relatively low solvation energy for estrogen and possesses a surface charge similar to that of estrogen, and thus can effectively enrich estrogen. The study provides ideas and methods for the extraction and determination of low-concentration drugs in biological samples and also provides a theoretical basis for the application of non-ionic hydrophobic deep eutectic solvent extraction.

Cavity formation at metal-water interfaces

The free energy cost of forming a cavity in a solvent is a fundamental concept in rationalizing the solvation of molecules and ions. A detailed understanding of the factors governing cavity formation in bulk solutions has inter alia enabled the formulation of models that account for this contribution in coarse-grained implicit solvation methods. Here, we employ classical molecular dynamics simulations and multistate Bennett acceptance ratio free energy sampling to systematically study cavity formation at a wide range of metal-water interfaces. We demonstrate that the obtained size- and position-dependence of cavitation energies can be fully rationalized by a geometric Gibbs model, which considers that the creation of the metal-cavity interface necessarily involves the removal of interfacial solvent. This so-called competitive adsorption effect introduces a substrate dependence to the interfacial cavity formation energy that is missed in existing bulk cavitation models. Using expressions from scaled particle theory, this substrate dependence is quantitatively reproduced by the Gibbs model through simple linear relations with the adsorption energy of a single water molecule. Besides providing a better general understanding of interfacial solvation, this paves the way for the derivation and efficient parametrization of more accurate interface-aware implicit solvation models needed for reliable high-throughput calculations toward improved electrocatalysts.

A Combined Reaction Path Search and Hybrid Solvation Method for the Systematic Exploration of Elementary Reactions at the Solid-Liquid Interface

We present a combined simulation method of single-component artificial force induced reaction (SC-AFIR) and effective screening medium combined with the reference interaction site model (ESM-RISM), termed SC-AFIR+ESM-RISM. SC-AFIR automatically and systematically explores the chemical reaction pathway, and ESM-RISM directly simulates the precise electronic structure at the solid-liquid interface. Hence, SC-AFIR+ESM-RISM enables us to explore reliable reaction pathways at the solid-liquid interface. We applied it to explore the dissociation pathway of an H2O molecule at the Cu(111)/water interface. The reaction path networks of the whole reaction and the minimum energy paths from H2O to H2 + O depend on the interfacial environment. The qualitative difference in the energy diagrams and the resulting change in the kinematically favored dissociation pathway upon changing the solvation environments are discussed. We believe that SC-AFIR+ESM-RISM will be a powerful tool to reveal the details of chemical reactions in surface catalysis and electrochemistry.

Unraveling the reaction mechanisms for furfural electroreduction on copper

Electrochemical routes for the valorization of biomass-derived feedstock molecules offer sustainable pathways to produce chemicals and fuels. However, the underlying reaction mechanisms for their electrochemical conversion remain elusive. In particular, the exact role of proton-electron coupled transfer and electrocatalytic hydrogenation in the reaction mechanisms for biomass electroreduction are disputed. In this work, we study the reaction mechanism underlying the electroreduction of furfural, an important biomass-derived platform chemical, combining grand-canonical (constant-potential) density functional theory-based microkinetic simulations and pH dependent experiments on Cu under acidic conditions. Our simulations indicate the second PCET step in the reaction pathway to be the rate- and selectivity-determining step for the production of the two main products of furfural electroreduction on Cu, i.e., furfuryl alcohol and 2-methyl furan, at moderate overpotentials. We further identify the source of Cu's ability to produce both products with comparable activity in their nearly equal activation energies. Furthermore, our microkinetic simulations suggest that surface hydrogenation steps play a minor role in determining the overall activity of furfural electroreduction compared to PCET steps due to the low steady-state hydrogen coverage predicted under reaction conditions, the high activation barriers for surface hydrogenation and the observed pH dependence of the reaction. As a theoretical guideline, low pH (<1.5) and moderate potential (ca. -0.5 V vs. SHE) conditions are suggested for selective 2-MF production.