The Fundamental Noninnocent Role of Water for the Hydrogenation of

Nov 17, 2017 - The latter catalyst displayed the highest TON for the reaction that converts N2O and H2 into innocuous N2 and H2O. This catalyst bears ...
0 downloads 10 Views 1MB Size
Communication Cite This: Inorg. Chem. XXXX, XXX, XXX-XXX

pubs.acs.org/IC

The Fundamental Noninnocent Role of Water for the Hydrogenation of Nitrous Oxide by PNP Pincer Ru-based Catalysts Jesús A. Luque-Urrutia‡ and Albert Poater*,‡ ‡

Institut de Química Computacional i Catàlisi and Departament de Química, Universitat de Girona, Campus Montilivi, 17003 Girona, Catalonia, Spain S Supporting Information *

Scheme 1. Hydrogenation of N2O by the PNP-pincer Ru based catalyst 1 (P = P(iPr)2)

ABSTRACT: The hydrogenation of nitrous oxide by PNP pincer ruthenium complexes supposes a promising way to functionalize a hazardous gas and reduce the greenhouse effect, generating dinitrogen and water. Here, by DFT calculations we describe not only the whole mechanism for such a green transformation but we unravel the fundamental role of water, without which the reaction could not go forward. Water assists mandatorily in the H transfer to generate the hydroxyl group together with the release of dinitrogen.

T

he greenhouse effect (GE) has become a worldwide problem that keeps on affecting the world climate and the beings that live in it. The main perpetrators for this phenomenon are the greenhouse gases: carbon dioxide (CO2), methane (CH4), ozone-depleting substances (ODS), hydrofluorocarbons (HFCs), sulfur hexafluoride (SF6), perfluorocarbons (PFCs), and nitrous oxide (N2O).1,2 The gas with the most impact right now is CO2, being over 60% of the GE cause; however, the other ones should not be overlooked. N2O has been making ground over the years due to the increase in nitrogenized substances for our daily life,3 such as fertilizers and pesticides or burning fossil fuels.4 These processes add N2O to the atmosphere, and due to its steadystate lifetime (120 years), it has been increasing for decades.5,6 The reductive process that consists of the hydrogenation of N2O by means of H2 to generate innocuous N2 and H2O has turned out to be feasible by heterogeneous7 and homogeneous8,9 catalytic systems. For the latter systems, the reaction was indeed not catalytic, but simply stoichiometric,8 and further, computationally, this reaction was scarcely studied. Thus, key mechanistic insights are missing. In order to address these drawbacks, and following the recent Milstein et al. catalytic work,9 we have focused our efforts on unveiling the full mechanism for catalyst 1 in Scheme 1. The latter catalyst displayed the highest TON for the reaction that converts N2O and H2 into innocuous N2 and H2O. This catalyst bears a PNP pincer ligand,10,11,12 which consists of an equatorial tridentate ligand that binds to the ruthenium center by two phosphorus atoms and one nitrogen atom.13−15 Recently, the computational studies by Gonçalves and Huang turned out to be fundamental for the rationalization of the reactivity of pincer based complexes from an aromatic point of view.16 These pincer type ligands have shown many experimental applications such as the formation of NH3 starting from N217 and what © XXXX American Chemical Society

could be its counterpart, the cleavage of the N−H bond in NH318 as well as many others. In terms of transition metal centers, Fe,19,20 Ru,21,22 Ir,23,24 Mo,25,26 Ni,27,28 W,29,30 and Pd31,32 have been studied, showing a growing interest in this kind of pincer ligands and a high applicability to many areas.15 From the initial catalyst 1, i.e., a nearly symmetrical species, a N2O molecule is added by the attack of its oxygen to the Ru−H moiety, the bond of which breaks at the same time, and consequently the former hydride protonates the terminal nitrogen atom (see Figure 1).33 Overall, the transition state that defines this concerted attack of the entering N2O turns out to be the highest energy barrier, with 26.7 kcal/mol. The next step implies the formation of the hydroxo ligand on ruthenium with the concerted release of a N2 molecule. However, this step is again not trivial. It requires an energy amount of 33.5 kcal/ mol (see Figure 2a), whereas when assisted by a water molecule, the energy barrier decreases to 25.5 kcal/mol (see Figure 2b). This means that in order to go forward, traces of water in the solvent used, THF, are mandatory.34,35 THF was found to not assist this step, even not facilitating it. Anyway, the barrier corresponding to the I → I−II + H2O step represents the rate-determining step (rds) since it describes the upper energy point of the whole reaction pathway, even though the previous energy barrier corresponding to the step 1 → 1−I is slightly higher in energy by 1.2 kcal/mol.36 Once the hydroxo species II is formed, any of both methylidene moieties of the pincer ligand provides a proton, generating the aqua complex III, overcoming a facile transition state. Then, the H2O liberation means the generation of catalyst 2. From 2, the coordination of H2 to the vacant position of the ruthenium center, followed by a H-transfer to the close Received: October 13, 2017

A

DOI: 10.1021/acs.inorgchem.7b02630 Inorg. Chem. XXXX, XXX, XXX−XXX

Communication

Inorganic Chemistry

Figure 1. Full mechanism for the hydrogenation of N2O. Isopropyl ligand represented. tert-Butyl and phenyl ligands’ energy shown as well (green and pink energies respectively). (Gibbs free energies in kcal/mol, and referred to catalyst 1.)

For the sake of a direct comparison, after replicating their calculations (M05-2X/6-311++G(2d,p) ∼ sdd //M05-2X/631G(d,p) ∼ sdd) with the computational scheme used here (M06/aug-cc-PVTZ ∼ sdd//M06/TZVP ∼ sdd), we obtained an energy difference of 4.7 kcal/mol. We attribute the increase of the preference for 1 with respect to 3 to the increased flexibility of the pincer ligand that bears 1, thanks to the second P atom, over a N atom of 3. This is explained due to the increased preference for the sp 3 hybridization of the phosphorus over the nitrogen atom, together with the smaller size of the latter atom (98 pm for P vs 65 pm for N in Bohr’s atomic radius model).38 Alternative pathways were also tested. First, instead of coordination to the ruthenium center of N2O via the O atom, it was done via the terminal N atom, which turns out to be 12.9 kcal/mol more expensive. Then, without the presence of traces of water, the formation of the isomer I′ is favored by 1.4 kcal/ mol with respect to the direct formation of II (see Figure 1). This step has to be considered as a possible catalyst death since it leads to a peculiar Ru−P bond cleavage, where the oxygen atom is inserted between both atoms (see Figure 3). Once we achieve this undesired species I′, the energy required to reorder the atoms in order to generate II is 58.2 kcal/mol, which is an energy barrier too high bearing the experimental conditions. Thus, this confirms that at an initial stage, traces of water must assist the mechanism development before water molecules formed as a byproduct of the hydrogenation of N2O are released. Furthermore, the assistance of either more water or THF molecules was tested. Among the computed data, for the hydrogenation step IV → 1, a water molecule increased the energy barrier by 25.9 kcal/mol, while a second water molecule also increased the energy barrier of steps I → II and II → III by 9.1 and 2.0 kcal/mol, respectively. Actually, going into detail,

Figure 2. Transition states I → II (a) not assisted or (b) assisted by a water molecule (most of hydrogens omitted for clarity, selected distances given in Å).

methylylidene of the pincer ligand, generates the dihydride species 1, thus closing the catalytic cycle, in a strongly exergonic manner (71.7 kcal/mol). Here, the first thing that should be remarked is that the thermodynamic equilibrium between catalysts 1 and 2 with intermediate IV in between is strongly displaced toward 1, by 14.3 kcal/mol. This was unexpected since the same equilibrium for complex 3 of Wang et al. (see Scheme 2),37 which resembles our catalyst, was found to be only slightly displaced to the hydrogenated species, by only 2.2 kcal/mol. Scheme 2. Catalyst Screening for Hydrogenation of N2O

B

DOI: 10.1021/acs.inorgchem.7b02630 Inorg. Chem. XXXX, XXX, XXX−XXX

Communication

Inorganic Chemistry

therefore, we can correlate this increase in bulkiness to its poorer catalytic performance. However, we noticed that system 2-Ph is indeed less hindered than 2, yet it is slightly worse. This poorer performance bearing phenyls may be due to the stronger C−P double bond that exists between the phosphorus moiety and the pyridine of the PNP pincer in 2-Ph that later must be hydrogenated to regenerate the catalytic hydride species 1. According to the Mayer bond order data for the bond between the P atom and the methylylidene fragment that gets hydrogenated, we observed a clear difference: for 2 being 1.079, while for 2-Ph it is 1.131, whereas upon hydrogenation of the CH moiety, MBOs for species 1 are 0.894 for the Ph and 0.881 for iPr. We have determined the whole reaction mechanism for the pincer based catalyst 1. The possibility of a water autocatalysis at the initial stages of the reaction was determined, while it would be hindered as time progresses and water is formed as a product of the hydrogenation of N2O. In complete absence of water, a peculiar bond cleavage of any of the Ru−P bonds in order to form a Ru−O−P moiety has been computationally highly stable, and we consider it a plausible catalyst poisoning due to the high energy required in order to continue the mechanism and close the catalytic cycle. Between the three studied ligands on P atoms, iPr is indeed the best one, with tBu having increased bulkiness and Ph gaining some conjugation between the phenyls and the pyridine, thus, requiring more energy to overcome the upper barrier of the whole mechanism. Further studies are planned to in silico predict a modified PNP pincer Ru-based catalyst that does not need the presence of any water molecule, together with the control of the potential Ru− P bond cleavage that here has been demonstrated to lead toward the decomposition of the catalyst.

Figure 3. Alternative pathway through intermediate I′: (a) transition state I → I′, (b) intermediate I′ and (c) transition state I′ → II (most of hydrogens omitted for clarity, selected distances given in Å).

the six-membered cycle of the water assisted transition state I → II displayed in Figure 2b is fundamental to explaining its relative high stability. To understand the nature of the pincer ligand and its role in the catalytic studied reaction here, we have also performed an analysis of catalyst 1 replacing the iPr ligands bonded to the phosphorus atoms of the pincer ligand by tert-butyl (tBu) groups (1-tBu) and phenyl (Ph) ligands (1-Ph). The results for both latter catalysts, with poorer catalytic performance with respect to 1, are compared in Figure 1, as well. According to our data, we have found that, just as Milstein et al. found, the isopropyl ligands are better than the tert-butyl and phenyl ligands. Looking into the energy data in Figure 1, we can see that 1-tBu and 1-Ph are both higher in energy in all cases (except III for 1-Ph, vide infra) going up to a maximum energy difference of 12.9 kcal (I−II) for 1-tBu and 7.7 kcal/ mol (I) for 1-Ph. Even though the iPr and the Ph ligands quite resemble one another, for the tBu there are increased differences with respect to the other ligands, which is caused by the higher steric hindrance of the latter ligand. As we can see in Figure 4, the steric hindrance of the tBu groups is remarkable compared to the Ph and the iPr ones;



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.inorgchem.7b02630. Computational details and all XYZ coordinates, energies and 3D structures of all species (PDF)



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

Albert Poater: 0000-0002-8997-2599 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS J.A.L.-U. thanks Universitat de Girona for a IFUdG2017 predoctoral fellowship. A.P. thanks the Spanish MINECO for a project CTQ2014-59832-JIN, and the EU for a FEDER fund (UNGI08-4E-003)). We thank Profs. David Milstein and Rong Zheng for helpful comments.



REFERENCES

(1) Freedman, B.; Frost, R. Gale Encyclopedia of Science; Gale Group, 1995; pp 1875−1880. (2) Rodhe, H. A Comparison of the Contribution of Various Greenhouse Gases to the Greenhouse Effect. Science 1990, 248, 1217− 1219.

Figure 4. Steric maps for the pincer ligands of the catalysts 1, 1-tBu, 1Ph. Ru metal center and CO ligand omitted for clarity (tBu center is out of scale). C

DOI: 10.1021/acs.inorgchem.7b02630 Inorg. Chem. XXXX, XXX, XXX−XXX

Communication

Inorganic Chemistry (3) (a) Konsolakis, M. Recent Advances on Nitrous Oxide (N2O) Decomposition over Non-Noble-Metal Oxide Catalysts: Catalytic Performance, Mechanistic Considerations, and Surface Chemistry Aspects. ACS Catal. 2015, 5, 6397−6421. (b) Severin, K. Synthetic Chemistry with Nitrous Oxide. Chem. Soc. Rev. 2015, 44, 6375−6386. (c) Parmon, V. N.; Panov, G. I.; Uriarte, A.; Noskov, A. S. Nitrous Oxide in Oxidation Chemistry and Catalysis: Application and Production. Catal. Today 2005, 100, 115−131. (d) Tolman, W. B. Angew. Chem., Int. Ed. 2010, 49, 1018−1024. (4) López, J. C.; Quijano, G.; Souza, T. S. O.; Estrada, J. M.; Lebrero, R.; Muñoz, R. Biotechnologies for Greenhouse Gases (CH4, N2O, and CO2) Abatement: State of the Art and Challenges. Appl. Microbiol. Biotechnol. 2013, 97, 2277−2303. (5) Hsu, J.; Prather, M. J. Global Long-Lived Chemical Modes Excited in a 3-D Chemistry Transport Model: Stratospheric N2O, NOy, O3 and CH4 Chemistry. Geophys. Res. Lett. 2010, 37, L07805. (6) Montzka, S. A.; Dlugokencky, E. J.; Butler, J. H. Non-CO2 Greenhouse Gases and Climate Change. Nature 2011, 476, 43−50. (7) (a) Benton, A. F.; Thacker, C. M. The Kinetics of Reaction between Nitrous Oxide and Hydrogen at a Silver Surface. J. Am. Chem. Soc. 1934, 56, 1300−1304. (b) Dixon, J. K.; Vance, J. E. The Reaction between Nitrous Oxide and Hydrogen on Platinum. J. Am. Chem. Soc. 1935, 57, 818−821. (c) Miyamoto, A.; Baba, S.; Mori, M.; Murakami, Y. Ion-Solvent Molecule Interactions in the Gas Phase. The Potassium Ion and Benzene. J. Phys. Chem. 1981, 85, 3117−3122. (d) Delahay, G.; Mauvezin, M.; Guzmán-Vargas, A.; Coq, B. Effect of the Reductant Nature on the Catalytic Removal of N2O on Fe-Zeolite-β Catalysts. Catal. Commun. 2002, 3, 385−389. (e) Nobukawa, T.; Yoshida, M.; Okumura, K.; Tomishige, K.; Kunimori, K. Effect of Reductants in N2O Reduction Over Fe-MFI Catalysts. J. Catal. 2005, 229, 374−388. (8) (a) Kaplan, A. W.; Bergman, R. G. Nitrous Oxide Mediated Oxygen Atom Insertion into a Ruthenium−Hydride Bond. Synthesis and Reactivity of the Monomeric Hydroxoruthenium Complex (DMPE) 2Ru(H)(OH). Organometallics 1997, 16, 1106−1108. (b) Lee, J.-H.; Pink, M.; Tomaszewski, J.; Fan, H.; Caulton, K. G. Facile Hydrogenation of N2O by an Operationally Unsaturated Osmium Polyhydride. J. Am. Chem. Soc. 2007, 129, 8706−8707. (c) Doyle, L. E.; Piers, W. E.; Borau-Garcia, J. Ligand Cooperation in the Formal Hydrogenation of N2O Using a PCsp2P Iridium Pincer Complex. J. Am. Chem. Soc. 2015, 137, 2187−2190. (d) Gianetti, T. L.; Annen, S. P.; Santiso-Quinones, G.; Reiher, M.; Driess, M.; Grützmacher, H. Nitrous Oxide as a Hydrogen Acceptor for the Dehydrogenative Coupling of Alcohols. Angew. Chem., Int. Ed. 2016, 55, 1854−1858. (e) Gianetti, T. L.; Rodríguez-Lugo, R. E.; Harmer, J. F.; Trincado, M.; Vogt, M.; Santiso-Quinones, G.; Grützmacher, H. Zero-Valent Amino-Olefin Cobalt Complexes as Catalysts for Oxygen Atom Transfer Reactions from Nitrous Oxide. Angew. Chem., Int. Ed. 2016, 55, 15323−15328. (9) Zeng, R.; Feller, M.; Ben-David, Y.; Milstein, D. Hydrogenation and Hydrosilylation of Nitrous Oxide Homogeneously Catalyzed by a Metal Complex. J. Am. Chem. Soc. 2017, 139, 5720−5723. (10) Yu, H.; Jia, G.; Lin, Z. Theoretical Studies on O-Insertion Reactions of Nitrous Oxide with Ruthenium Hydride Complexes. Organometallics 2008, 27, 3825−3833. (11) Zhang, J.; Balaraman, E.; Leitus, G.; Milstein, D. Electron-Rich PNP- and PNN-Type Ruthenium(II) Hydrido Borohydride Pincer Complexes. Synthesis, Structure, and Catalytic Dehydrogenation of Alcohols and Hydrogenation of Esters. Organometallics 2011, 30, 5716−5724. (12) Zell, T.; Ben-David, Y.; Milstein, D. Highly Efficient, General Hydrogenation of Aldehydes Catalyzed by PNP Iron Pincer Complexes. Catal. Sci. Technol. 2015, 5, 822−826. (13) Schneider, S.; Meiners, J.; Askevold, B. Cooperative Aliphatic PNP Amido Pincer Ligands − Versatile Building Blocks for Coordination Chemistry and Catalysis. Eur. J. Inorg. Chem. 2012, 2012, 412−429. (14) van der Vlugt, J. I. Cooperative Catalysis with First-Row Late Transition Metals. Eur. J. Inorg. Chem. 2012, 2012, 363−375.

(15) Qu, S.; Dang, Y.; Song, C.; Wen, M.; Huang, K.; Wang, Z. Catalytic Mechanisms of Direct Pyrrole Synthesis via Dehydrogenative Coupling Mediated by PNP-Ir or PNN-Ru Pincer Complexes: Crucial Role of Proton-Transfer Shuttles in the PNP-Ir System. J. Am. Chem. Soc. 2014, 136, 4974−4991. (16) Gonçalves, T. P.; Huang, K. Metal-Ligand Cooperative Reactivity in the (Pseudo)-Dearomatized PNx(P) Systems: The Influence of the Zwitterionic Form in Dearomatized Pincer Complexes. J. Am. Chem. Soc. 2017, 139, 13442−13449. (17) Kuriyama, S.; Arashiba, K.; Nakajima, K.; Tanaka, H.; Kamaru, N.; Yoshizawa, K.; Nishibayashi, Y. Catalytic Reduction of Dinitrogen to Ammonia by Use of Molybdenum−Nitride Complexes Bearing a Tridentate Triphosphine as Catalysts. J. Am. Chem. Soc. 2014, 136, 9719−9731. (18) Chang, Y.; Nakajima, Y.; Tanaka, H.; Yoshizawa, K.; Ozawa, F. Facile N−H Bond Cleavage of Ammonia by an Iridium Complex Bearing a Noninnocent PNP-Pincer Type Phosphaalkene Ligand. J. Am. Chem. Soc. 2013, 135, 11791−11794. (19) Fillman, K. L.; Bielinski, E. A.; Schmeier, T. J.; Nesvet, J. C.; Woodruff, T. M.; Pan, C. J.; Takase, M. K.; Hazari, N.; Neidig, M. L. Flexible Binding of PNP Pincer Ligands to Monomeric Iron Complexes. Inorg. Chem. 2014, 53, 6066−6072. (20) Mellone, I.; Gorgas, N.; Bertini, F.; Peruzzini, M.; Kirchner, K.; Gonsalvi, L. Selective Formic Acid Dehydrogenation Catalyzed by FePNP Pincer Complexes Based on the 2,6-Diaminopyridine Scaffold. Organometallics 2016, 35, 3344−3349. (21) Nakajima, Y.; Okamoto, Y.; Chang, Y.; Ozawa, F. Synthesis, Structures, and Reactivity of Ruthenium Complexes with PNP-pincer Type Phosphaalkene Ligands. Organometallics 2013, 32, 2918−2925. (22) Alberico, E.; Lennox, A. J. J.; Vogt, L. K.; Jiao, H.; Baumann, W.; Drexler, H.; Nielsen, M.; Spannenberg, A.; Checinski, M. P.; Junge, H.; Beller, M. Unravelling the Mechanism of Basic Aqueous Methanol Dehydrogenation Catalyzed by Ru−PNP Pincer Complexes. J. Am. Chem. Soc. 2016, 138, 14890−14904. (23) Cheng, C.; Kim, B. G.; Guironnet, D.; Brookhart, M.; Guan, C.; Wang, D. Y.; Krogh-Jespersen, K.; Goldman, A. S. Synthesis and Characterization of Carbazolide-Based Iridium PNP Pincer Complexes. Mechanistic and Computational Investigation of Alkene Hydrogenation: Evidence for an Ir(III)/Ir(V)/Ir(III) Catalytic Cycle. J. Am. Chem. Soc. 2014, 136, 6672−6683. (24) Feller, M.; Ben-Ari, E.; Diskin-Posner, Y.; Carmieli, R.; Weiner, L.; Milstein, D. O2 Activation by Metal−Ligand Cooperation with IrI PNP Pincer Complexes. J. Am. Chem. Soc. 2015, 137, 4634−4637. (25) Ö ztopcu, Ö .; Holzhacker, C.; Puchberger, M.; Weil, M.; Mereiter, K.; Veiros, L. F.; Kirchner, K. Synthesis and Characterization of Hydrido Carbonyl Molybdenum and Tungsten PNP Pincer Complexes. Organometallics 2013, 32, 3042−3052. (26) Wei, Z.; Junge, K.; Beller, M.; Jiao, H. Hydrogenation of PhenylSubstituted C ≡ N, C = N, C ≡ C, C = C and C = O Functional Groups by Cr, Mo and W PNP Pincer Complexes − a DFT Study. Catal. Sci. Technol. 2017, 7, 2298−2307. (27) Vogt, M.; Rivada-Wheelaghan, O.; Iron, M. A.; Leitus, G.; Diskin-Posner, Y.; Shimon, L. J. W.; Ben-David, Y.; Milstein, D. Anionic Nickel(II) Complexes with Doubly Deprotonated PNP Pincer-type Ligands and Their Reactivity Toward CO2. Organometallics 2013, 32, 300−308. (28) Kreye, M.; Freytag, M.; Jones, P. G.; Williard, P. G.; Bernskoetter, W. H.; Walter, M. D. Homolytic H2 Cleavage by a Mercury-Bridged Ni(i) Pincer Complex [{(PNP)Ni}2{μ-Hg}]. Chem. Commun. 2015, 51, 2946−2949. (29) Moha, V.; Leitner, W.; Hölscher, M. Hölscher M. H3 Synthesis in the N2/H2 Reaction System using Cooperative Molecular Tungsten/Rhodium Catalysis in Ionic Hydrogenation: A DFT Study. Chem. - Eur. J. 2016, 22, 2624−2628. (30) Arashiba, K.; Sasaki, K.; Kuriyama, S.; Miyake, Y.; Nakanishi, H.; Nishibayashi, Y. Synthesis and Protonation of Molybdenum− and Tungsten−Dinitrogen Complexes Bearing PNP-Type Pincer Ligands. Organometallics 2012, 31, 2035−2041. D

DOI: 10.1021/acs.inorgchem.7b02630 Inorg. Chem. XXXX, XXX, XXX−XXX

Communication

Inorganic Chemistry (31) Bailey, W. D.; Kaminsky, W.; Kemp, R. A.; Goldberg, K. I. Synthesis and Characterization of Anionic, Neutral, and Cationic PNP Pincer PdII and PtII Hydrides. Organometallics 2014, 33, 2503−2509. (32) Yang, Z.; Xia, C.; Liu, D.; Liu, Y.; Sugiya, M.; Imamoto, T.; Zhang, W. P-Stereogenic PNP Pincer-Pd Catalyzed Intramolecular Hydroamination of Amino-1,3-dienes. Org. Biomol. Chem. 2015, 13, 2694−2702. (33) DFT static calculations were performed with Gaussian09 using the BP86-D3BJ functional, and the SDD basis set for ruthenium and the TZVP basis set for the other atoms. Single-point calculations on the optimized geometries with M06 functional and the aug-cc-PVTZ basis set. Solvent effects were included by PCM using THF as a solvent. (34) Values at a pressure of 1354 atm using the approach of Martin and co-workers, excluding the potential overestimation of the entropy contribution, confirm the message of the paper since the barrier assisted by a water molecule in step I→II is again 8.0 kcal/mol exactly lower in energy, whereas the undesired step I→I′ is still favored by 1.3 kcal/mol without the assistance of a water molecule. Further, the no reversibility character of the 2→1 step is confirmed, since the energy difference increases 4.6 kcal/mol. (35) (a) Martin, R. L.; Hay, P. J.; Pratt, L. R. Hydrolysis of Ferric Ion in Water and Conformational Equilibrium. J. Phys. Chem. A 1998, 102, 3565−3573. (b) Manzini, S.; Poater, A.; Nelson, D. J.; Cavallo, L.; Slawin, A. M. Z.; Nolan, S. P. Insights into the Decomposition of Olefin Metathesis Pre-Catalysts. Angew. Chem., Int. Ed. 2014, 53, 8995−8999. (c) Poater, A.; Pump, E.; Vummaleti, S. V. C.; Cavallo, L. The Right Computational Recipe for Olefin Metathesis with Ru-Based Catalysts: the Whole Mechanism of Ring-Closing Olefin Metathesis. J. Chem. Theory Comput. 2014, 10, 4442−4448. (36) Kozuch, S.; Shaik, S. How to Conceptualize Catalytic Cycles? The Energetic Span Model. Acc. Chem. Res. 2011, 44, 101−110. (37) Lu, G.; Li, H.; Zhao, L.; Huang, F.; Schleyer, P.; Wang, Z.-X. v. R.; Wang, Z. Designing Metal-Free Catalysts by Mimicking TransitionMetal Pincer Templates. Chem. - Eur. J. 2011, 17, 2038−2043. (38) Bohr, N. On the Constitution of Atoms and Molecules. Phys. Mag. 1913, 26, 1−25.

E

DOI: 10.1021/acs.inorgchem.7b02630 Inorg. Chem. XXXX, XXX, XXX−XXX