How Solvent Affects C-H Activation and Hydrogen Production Pathways in Homogeneous Ru-Catalyzed Methanol Dehydrogenation Reactions

Production of hydrogen from alcohols via homogeneous catalytic transformations mediated by molecular transition-metal complexes

Hydrogen obtained from renewable sources such as water and alcohols is regarded as an efficient clean-burning alternative to non-renewable fuels. The use of the so-called bio-H2 regardless of its colour will be a significant step towards achieving global net-zero carbon goals. Challenges still persist however with conventional H2 storage, which include low-storage density and high cost of transportation apart from safety concerns. Global efforts have thus focussed on liquid organic hydrogen carriers (LOHCs), which have shown excellent potential for H2 storage while allowing safer large-scale transformation and easy on-site H2 generation. While water could be considered as the most convenient liquid inorganic hydrogen carrier (LIHC) on a long-term basis, the utilization of alcohols as LOHCs to generate on-demand H2 has tasted instant success. This has helped to draw a road-map of futuristic H2 storage and transportation. The current review brings to the fore the state-of-the-art developments in hydrogen generation from readily available, feed-agnostic bio-alcohols as LOHCs using molecular transition-metal catalysts.

Merging electrocatalytic alcohol oxidation with C-N bond formation by electrifying metal-ligand cooperative catalysts

Electrification of thermal chemical processes could play an important role in creating a more energy efficient chemical sector. Here we demonstrate that a range of MLC catalysts can be successfully electrified and used for imine formation from alcohol precursors, thus demonstrating the first example of molecular electrocatalytic C-N bond formation.This novel concept allowed energy efficiency to be increased by an order of magnitude compared to thermal catalysis. Molecular EAO and the electrification of homogeneous catalysts can thus contribute to current efforts for the electrocatalytic generation of C-N bonds from simple building blocks.

A switchable route for selective transformation of ethylene glycol to hydrogen and glycolic acid using a bifunctional ruthenium catalyst

The developed bifunctional NNN-Ru complex features a high catalytic efficiency for the selective production of hydrogen and glycolic acid from ethylene glycol under mild reaction conditions, where a TON of 6395 was achieved. Tuning the reaction conditions afforded further dehydrogenation of the organic substrate with higher hydrogen production, and a higher TON of 25 225 was attained. The scale-up reaction yielded 1230 mL of pure hydrogen gas under the optimized reaction conditions. The role of the bifunctional catalyst was studied and mechanistic investigations were performed.

Theoretical investigation of solvent and oxidation/deprotonation effects on the electronic structure of a mononuclear Ru-aqua-polypyridine complex in aqueous solution

Mononuclear polypyridine ruthenium (Ru) complexes can catalyze various reactions, including water splitting, and can also serve as photosensitizers in solar cells. Despite recent progress in their synthesis, accurately modeling their physicochemical properties, particularly in solution, remains challenging. Herein, we conduct a theoretical investigation of the structural and electronic properties of a mononuclear Ru-aqua polypyridine complex in aqueous solution, considering five of its possible oxidation/protonation states species: [RuII(H2O)(py)(bpy)2]2+, [RuII(OH)(py)(bpy)2]+, [RuIII(H2O)(py)(bpy)2]3+, [RuIII(OH)(py)(bpy)2]2+ and [RuIV(O)(py)(bpy)2]2+, where py = pyridine and bpy = 2,2'-bipyridine. At first, we investigate the impact of proton-coupled and non-coupled electron transfer reactions on the geometry and electronic structure of the complexes in vacuum and in solution, using an implicit solvent model. Then, using a sequential multiscale approach that combines quantum mechanics and molecular mechanics (S-QM/MM), we examine the explicit solvent effects on the electronic excitations of the complexes, and compare them with the experimental results. The complexes were synthesized, and their absorption spectra measured in aqueous solution. To accurately describe the QM interactions between the metal center and the aqueous ligand in the MM simulations, we developed new force field parameters for the Ru atom. We analyze the solvent structure around the complexes and account for its explicit influence on the polarization and electronic excitations of the complexes. Notably, accounting for the explicit solvent polarization effects of the first solvation shells is essential to correctly describe the energy of the electronic transitions, and the explicit treatment of the hydrogen bonds at the QM level in the excitation calculations improves the accuracy of the description of the metal-to-ligand charge-transfer bands. Transition density matrix analysis is used to characterize all electronic transitions in the visible and ultraviolet ranges according to their charge-transfer (CT) character. This study elucidates the electronic structure of those ruthenium polypyridyl complexes in aqueous solution and underscores the importance of precisely describing solvent effects, which can be achieved employing the S-QM/MM method.

Ab Initio Metadynamics Simulations of Hexafluoroisopropanol Solvent Effects: Synergistic Role of Solvent H-Bonding Networks and Solvent-Solute C-H/π Interactions

The solvent effects in Friedel-Crafts cycloalkylation of epoxides and Cope rearrangement of aldimines were investigated by using ab initio molecular dynamics simulations. Explicit molecular treatments were applied for both reactants and solvents. The reaction mechanisms were elucidated via free energy calculations based on metadynamics simulations. The results reveal that both reactions proceed in a concerted fashion. Key solvent-substrate interactions are identified from the structures of transition states with explicit solvent molecules. The remarkable promotion effect of hexafluoroisopropanol solvent is ascribed to the synergistic effect of H-bonding networks and C-H/π interactions with substrates.

The Mechanism of Markovnikov-Selective Epoxide Hydrogenolysis Catalyzed by Ruthenium PNN and PNP Pincer Complexes

The homogeneous catalysis of epoxide hydrogenolysis to give alcohols has recently received significant attention. Catalyst systems have been developed for the selective formation of either the Markovnikov (branched) or anti-Markovnikov (linear) alcohol product. Thus far, the reported catalysts exhibiting Markovnikov selectivity all feature the potential for Noyori/Shvo-type bifunctional catalysis, with either a RuH/NH or FeH/OH core structure. The proposed mechanisms of epoxide ring-opening have involved cooperative C-O bond hydrogenolysis involving the metal hydride and the acidic pendant group on the ligand, in analogy to the well-documented mechanism of polar double-bond hydrogenation exhibited by catalysts of this type. In this work, we present a combined computational/experimental study of the mechanism of epoxide hydrogenolysis catalyzed by Noyori-type PNP and PNN complexes of ruthenium. We find that, at least for these ruthenium systems, the previously proposed bifunctional pathway for epoxide ring-opening is energetically inaccessible; instead, the ring-opening proceeds through opposite-side nucleophilic attack of the ruthenium hydride on the epoxide carbon, without the involvement of the ligand N-H group. For both catalyst systems, the rate law and overall barrier predicted by density functional theory (DFT) are consistent with the results from kinetic studies.

Electrification of a Milstein-type catalyst for alcohol reformation

Novel energy and atom efficiency processes will be keys to develop the sustainable chemical industry of the future. Electrification could play an important role, by allowing to fine-tune energy input and using the ideal redox agent: the electron. Here we demonstrate that a commercially available Milstein ruthenium catalyst (1) can be used to promote the electrochemical oxidation of ethanol to ethyl acetate and acetate, thus demonstrating the four electron oxidation under preparative conditions. Cyclic voltammetry and DFT-calculations are used to devise a possible catalytic cycle based on a thermal chemical step generating the key hydride intermediate. Successful electrification of Milstein-type catalysts opens a pathway to use alcohols as a renewable feedstock for the generation of esters and other key building blocks in organic chemistry, thus contributing to increase energy efficiency in organic redox chemistry.

Solvent-Assisted Ketone Reduction by a Homogeneous Mn Catalyst

The choice of a solvent and the reaction conditions often defines the overall behavior of a homogeneous catalytic system by affecting the preferred reaction mechanism and thus the activity and selectivity of the catalytic process. Here, we explore the role of solvation in the mechanism of ketone reduction using a model representative of a bifunctional Mn-diamine catalyst through density functional theory calculations in a microsolvated environment by considering explicit solvent and fully solvated ab initio molecular dynamics simulations for the key elementary steps. Our computational analysis reveals the possibility of a Meerwein–Ponndorf–Verley (MPV) type mechanism in this system, which does not involve the participation of the N–H moiety and the formation of a transition-metal hydride species in ketone conversion. This path was not previously considered for Mn-based metal–ligand cooperative transfer hydrogenation homogeneous catalysis. The MPV mechanism is strongly facilitated by the solvent molecules present in the reaction environment and can potentially contribute to the catalytic performance of other related catalyst systems. Calculations indicate that, despite proceeding effectively in the second coordination sphere of the transition-metal center, the MPV reaction path retains the enantioselectivity preference induced by the presence of the small chiral N,N′-dimethyl-1,2-cyclohexanediamine ligand within the catalytic Mn(I) complex.

Solvent-mediated outer-sphere CO2 electro-reduction mechanism over the Ag111 surface

The electrocatalytic CO₂ reduction reaction (CO₂RR) is one of the key technologies of the clean energy economy. Molecular-level understanding of the CO₂RR process is instrumental for the better design of electrodes operable at low overpotentials with high current density. The catalytic mechanism underlying the turnover and selectivity of the CO₂RR is modulated by the nature of the electrocatalyst, as well as the electrolyte liquid, and its ionic components that form the electrical double layer (EDL). Herein we demonstrate the critical non-innocent role of the EDL for the activation and conversion of CO₂ at a high cathodic bias for electrocatalytic conversion over a silver surface as a representative low-cost model cathode. By using a multiscale modeling approach we demonstrate that under such conditions a dense EDL is formed, which hinders the diffusion of CO₂ towards the Ag111 electrocatalyst surface. By combining DFT calculations and ab initio molecular dynamics simulations we identify favorable pathways for CO₂ reduction directly over the EDL without the need for adsorption to the catalyst surface. The dense EDL promotes homogeneous phase reduction of CO₂via electron transfer from the surface to the electrolyte. Such an outer-sphere mechanism favors the formation of formate as the CO₂RR product. The formate can undergo dehydration to CO via a transition state stabilized by solvated alkali cations in the EDL.

In addition to the commonly accepted inner-sphere mechanism for e⁻ transfer, we show that an outer-sphere electron transfer from the cathode to CO₂ is operable at high overpotentials.

Computational mechanistic studies of ruthenium catalysed methanol dehydrogenation

Homogeneous ruthenium catalysed methanol dehydrogenation could become a key reaction for hydrogen production in liquid fuel cells. In order to improve existing catalytic systems, mechanistic insight is paramount in directing future studies. Herein, we describe what computational mechanistic research has taught us so far about ruthenium catalysed dehydrogenation reactions. In general, two mechanistic pathways can be operative in these reactions: a metal-centered or a metal-ligand cooperative (Noyori-Morris type) minimum energy reaction pathway (MERP). Discerning between these mechanisms on the basis of computational studies has proven to be highly input dependent, and to circumvent pitfalls it is important to consider several factors, such as solvent effects, metal-ligand cooperativity, alternative geometries, and complex electronic structures of metal centres. This Frontiers article summarizes the reported computational research performed on ruthenium catalyzed dehydrogenation reactions performed in the past decade, and serves as a guide for future research.

Small-Pore Zeolite Membranes: A Review of Gas Separation Applications and Membrane Preparation

There have been significant advancements in small-pore zeolite membranes in recent years. With pore size closely related to many energy- or environment-related gas molecules, small-pore zeolite membranes have demonstrated great potential for the separation of some interested gas pairs, such as CO2/CH4, CO2/N2 and N2/CH4. Small-pore zeolite membranes share some characteristics but also have distinctive differences depending on their framework, structure and zeolite chemistry. Through this mini review, the separation performance of different types of zeolite membranes with respect to interested gas pairs will be compared. We aim to give readers an idea of membrane separation status. A few representative synthesis conditions are arbitrarily chosen and summarized, along with the corresponding separation performance. This review can be used as a quick reference with respect to the influence of synthesis conditions on membrane quality. At the end, some general findings and perspectives will be discussed.

Beyond Continuum Solvent Models in Computational Homogeneous Catalysis

In homogeneous catalysis solvent is an inherent part of the catalytic system. As such, it must be considered in the computational modeling. The most common approach to include solvent effects in quantum mechanical calculations is by means of continuum solvent models. When they are properly used, average solvent effects are efficiently captured, mainly those related with solvent polarity. However, neglecting atomistic description of solvent molecules has its limitations, and continuum solvent models all alone cannot be applied to whatever situation. In many cases, inclusion of explicit solvent molecules in the quantum mechanical description of the system is mandatory. The purpose of this article is to highlight through selected examples what are the reasons that urge to go beyond the continuum models to the employment of micro-solvated (cluster-continuum) of fully explicit solvent models, in this way setting the limits of continuum solvent models in computational homogeneous catalysis. These examples showcase that inclusion of solvent molecules in the calculation not only can improve the description of already known mechanisms but can yield new mechanistic views of a reaction. With the aim of systematizing the use of explicit solvent models, after discussing the success and limitations of continuum solvent models, issues related with solvent coordination and solvent dynamics, solvent effects in reactions involving small, charged species, as well as reactions in protic solvents and the role of solvent as reagent itself are successively considered.

Molecular Insights on Solvent Effects in Organic Reactions as Obtained through Computational Chemistry Tools

Molecular understanding of the role of protic solvents in a gamut of organic transformations can be developed using density functional and ab initio computational studies focused on the reaction mechanism. Inclusion of explicit solvent molecules in the vital TSs has been proven to be valuable toward improving the energetic estimates of organocatalytic as well as transition-metal-catalyzed organic reactions. Herein, we provide an overview of the importance of an explicit-implicit solvation model using a number of interesting examples.

Mechanism of iron complexes catalyzed in the N-formylation of amines with CO2 and H2: the superior performance of N-H ligand methylated complexes

CO2 hydrogenation into value-added chemicals not only offer an economically beneficial outlet but also help reduce the emission of greenhouse gases. Herein, the density functional theory (DFT) studies have been carried out on CO2 hydrogenation reaction for formamide production catalyzed by two different N-H ligand types of PNP iron catalysts. The results suggest that the whole mechanistic pathway has three parts: (i) precatalyst activation, (ii) hydrogenation of CO2 to generate formic acid (HCOOH), and (iii) amine thermal condensation to formamide with HCOOH. The lower turnover number (TON) of a bifunctional catalyst system in hydrogenating CO2 may attribute to the facile side-reaction between CO2 and bifunctional catalyst, which inhibits the generation of active species. Regarding the bifunctional catalyst system addressed in this work, we proposed a ligand participated mechanism due to the low pKa of the ligand N-H functional in the associated stage in the catalytic cycle. Remarkably, catalysts without the N-H ligand exhibit the significant transfer hydrogenation through the metal centered mechanism. Due to the excellent catalytic nature of the N-H ligand methylated catalyst, the N-H bond was not necessary for stabilizing the intermediate. Therefore, we confirmed that N-H ligand methylated catalysts allow for an efficient CO2 hydrogenation reaction compared to the bifunctional catalysts. Furthermore, the influence of Lewis acid and strong base on catalytic N-formylation were considered. Both significantly impact the catalytic performance. Moreover, the catalytic activity of PNMeP-based Mn, Fe and Ru complexes for CO2 hydrogenation to formamides was explored as well. The energetic span of Fe and Mn catalysts are much closer to the precious metal Ru, which indicates that such non-precious metal catalysts have potentially valuable applications.

Trends in computational molecular catalyst design

Computational methods have emerged as a powerful tool to augment traditional experimental molecular catalyst design by providing useful predictions of catalyst performance and decreasing the time needed for catalyst screening. In this perspective, we discuss three approaches for computational molecular catalyst design: (i) the reaction mechanism-based approach that calculates all relevant elementary steps, finds the rate and selectivity determining steps, and ultimately makes predictions on catalyst performance based on kinetic analysis, (ii) the descriptor-based approach where physical/chemical considerations are used to find molecular properties as predictors of catalyst performance, and (iii) the data-driven approach where statistical analysis as well as machine learning (ML) methods are used to obtain relationships between available data/features and catalyst performance. Following an introduction to these approaches, we cover their strengths and weaknesses and highlight some recent key applications. Furthermore, we present an outlook on how the currently applied approaches may evolve in the near future by addressing how recent developments in building automated computational workflows and implementing advanced ML models hold promise for reducing human workload, eliminating human bias, and speeding up computational catalyst design at the same time. Finally, we provide our viewpoint on how some of the challenges associated with the up-and-coming approaches driven by automation and ML may be resolved.

Alcohol-Assisted Hydrogenation of Carbon Monoxide to Methanol Using Molecular Manganese Catalysts

Alcohol-assisted hydrogenation of carbon monoxide (CO) to methanol was achieved using homogeneous molecular complexes. The molecular manganese complex [Mn(CO)2Br[HN(C2H4P i Pr2)2]] ([HN(C2H4P i Pr2)2] = MACHO- i Pr) revealed the best performance, reaching up to turnover number = 4023 and turnover frequency 857 h-1 in EtOH/toluene as solvent under optimized conditions (T = 150 °C, p(CO/H2) = 5/50 bar, t = 8-12 h). Control experiments affirmed that the reaction proceeds via formate ester as the intermediate, whereby a catalytic amount of base was found to be sufficient to mediate its formation from CO and the alcohol in situ. Selectivity for methanol formation reached >99% with no accumulation of the formate ester. The reaction was demonstrated to work with methanol as the alcohol component, resulting in a reactive system that allows catalytic "breeding" of methanol without any coreagents.

Homogeneous Reforming of Aqueous Ethylene Glycol to Glycolic Acid and Pure Hydrogen Catalyzed by Pincer-Ruthenium Complexes Capable of Metal-Ligand Cooperation

Glycolic acid is a useful and important α-hydroxy acid that has broad applications. Herein, the homogeneous ruthenium catalyzed reforming of aqueous ethylene glycol to generate glycolic acid as well as pure hydrogen gas, without concomitant CO₂ emission, is reported. This approach provides a clean and sustainable direction to glycolic acid and hydrogen, based on inexpensive, readily available, and renewable ethylene glycol using 0.5 mol % of catalyst. In-depth mechanistic experimental and computational studies highlight key aspects of the PNNH-ligand framework involved in this transformation.

Accurate and rapid prediction of pKa of transition metal complexes: semiempirical quantum chemistry with a data-augmented approach

Data-augmented high-throughput QM approach to compute pK_a of transition metal hydride complexes with hDFT accuracy and low cost.

Homogeneous and heterogeneous catalysts for hydrogenation of CO2 to methanol under mild conditions

This review summarizes the concepts, mechanisms, drawbacks and challenges of the state-of-the-art catalysis for CO₂ to MeOH under mild conditions. Thoughtful guidelines and principles for future research are presented and discussed.

Hydrogenative Depolymerization of Nylons

The widespread crisis of plastic pollution demands discovery of new and sustainable approaches to degrade robust plastics such as nylons. Using a green and sustainable approach based on hydrogenation, in the presence of a ruthenium pincer catalyst at 150 °C and 70 bar H₂, we report here the first example of hydrogenative depolymerization of conventional, widely used nylons and polyamides, in general. Under the same catalytic conditions, we also demonstrate the hydrogenation of a polyurethane to produce diol, diamine, and methanol. Additionally, we demonstrate an example where monomers (and oligomers) obtained from the hydrogenation process can be dehydrogenated back to a poly(oligo)amide of approximately similar molecular weight, thus completing a closed loop cycle for recycling of polyamides. Based on the experimental and density functional theory studies, we propose a catalytic cycle for the process that is facilitated by metal–ligand cooperativity. Overall, this unprecedented transformation, albeit at the proof of concept level, offers a new approach toward a cleaner route to recycling nylons.

Elucidating cation effects in homogeneously catalyzed formic acid dehydrogenation

In this work, we use density functional theory based molecular dynamics with an explicit description of methanol solvent to study the effect of cations on formic acid dehydrogenation catalyzed by a ruthenium PNP pincer complex (RuPNP). Formic acid dehydrogenation is a two step process that involves the reorientation of the formate moiety bound via its oxygen to a H bound intermediate, followed by the hydride transfer step to form CO2 and the hydrogenated catalyst. We find the reorientation step to proceed with a low barrier in methanol solvent and in the presence of a Li+ cation, while the hydride transfer is significantly hindered by the presence of cations (Li+ and K+). The cation seems to strongly stabilize the negatively charged formate moiety, hindering complete hydride transfer and resulting in a high barrier for this step. This study is a first step towards addressing the exact role of cations in formic acid dehydrogenation reactions.

Determination of pKa Values via ab initio Molecular Dynamics and its Application to Transition Metal-Based Water Oxidation Catalysts

The p K a values are important for the in-depth elucidation of catalytic processes, the computational determination of which has been challenging. The first simulation protocols employing ab initio molecular dynamics simulations to calculate p K a values appeared almost two decades ago. Since then several slightly different methods have been proposed. We compare the performance of various evaluation methods in order to determine the most reliable protocol when it comes to simulate p K a values of transition metal-based complexes, such as the here investigated Ru-based water oxidation catalysts. The latter are of high interest for sustainable solar-light driven water splitting, and understanding of the underlying reaction mechanism is crucial for their further development.

Modeling the Catalyst Activation Step in a Metal–Ligand Radical Mechanism Based Water Oxidation System

Designing catalysts for water oxidation (WOCs) that operate at low overpotentials plays an important role in developing sustainable energy conversion schemes. Recently, a mononuclear ruthenium WOC that operates via metal–ligand radical coupling pathway was reported, with a very low barrier for O–O bond formation, that is usually the rate-determining step in most WOCs. A detailed mechanistic understanding of this mechanism is crucial to design highly active oxygen evolution catalysts. Here, we use density functional theory based molecular dynamics (DFT-MD) with an explicit description of the solvent to investigate the catalyst activation step for the [Ru(bpy) 2 (bpy–NO)] 2 + complex, that is considered to be the rate-limiting step in the metal–ligand radical coupling pathway. We find that a realistic description of the solvent environment, including explicit solvent molecules and thermal motion, is crucial for an accurate description of the catalyst activation step, and for the estimation of the activation barriers.

Chemoenzymatic Hydrogen Production from Methanol through the Interplay of Metal Complexes and Biocatalysts

Microbial methylotrophic organisms can serve as great inspiration in the development of biomimetic strategies for the dehydrogenative conversion of C₁ molecules under ambient conditions. In this Concept article, a concise personal perspective on the recent advancements in the field of biomimetic catalytic models for methanol and formaldehyde conversion, in the presence and absence of enzymes and co-factors, towards the formation of hydrogen under ambient conditions is given. In particular, formaldehyde dehydrogenase mimics have been introduced in stand-alone C₁ -interconversion networks. Recently, coupled systems with alcohol oxidase and dehydrogenase enzymes have been also developed for in situ formation and decomposition of formaldehyde and/or reduced/oxidized nicotinamide adenine dinucleotide (NADH/ NAD⁺ ). Although C₁ molecules are already used in many industries for hydrogen production, these conceptual bioinspired low-temperature energy conversion processes may lead one day to more efficient energy storage systems enabling renewable and sustainable hydrogen generation for hydrogen fuel cells under ambient conditions using C₁ molecules as fuels for mobile and miniaturized energy storage solutions in which harsh conditions like those in industrial plants are not applicable.

DFT Provides Insight into the Additive-Free Conversion of Aqueous Methanol to Dihydrogen Catalyzed by [Ru(trop2dad)]: Importance of the (Electronic) Flexibility of the Diazadiene Moiety

The mechanism for complete dehydrogenation of aqueous methanol to CO₂ and three equivalents of H₂ catalyzed by [Ru(trop₂dad)] was investigated with DFT (trop₂dad = 1,4-bis(5H-dibenzo[a,d]cyclohepten-5-yl)-1,4-diazabuta-1,3-diene). To date, this is the only catalyst that promotes the acceptorless dehydrogenation of aqueous methanol in homogeneous phase under mild conditions without the addition of an additive (base, acid, or a secondary catalyst). A detailed understanding of the mechanism of this transformation may therefore be of significant importance for the conversion of liquid organic fuels. Previous computational studies using simplified models of the catalyst suggested entirely ligand-centered reaction pathways with rather high-energy barriers for complete dehydrogenation of aqueous methanol. These are, however, not consistent with the experimental data. In the present paper, we reveal a different reaction mechanism for aqueous methanol dehydrogenation that involves metal–ligand cooperativity involving the diazadiene (dad) ligand and has substantially lower barriers, in good agreement with the experimental data. The dad moiety of the ligand actively participates in the alcohol activation mechanism. In the first step of the reaction, the dad ligand rearranges from a σ- to a π-bound coordination mode. This adjusts the electronic structure of both the metal and the ligand, leading to an enhanced Brønsted basicity of the nitrogen centers and higher Lewis acidity of the ruthenium center. As a result, concerted proton-hydride transfer to/from metal-hydride and N-protonated dad-ligand moieties becomes possible, leading to low-barrier metal–ligand cooperative elementary steps for alcohol activation and H₂ elimination.