General H2 Activation Modes for Lewis Acid ... - ACS Publications

Jan 27, 2016 - ABSTRACT: A general mechanism for H2 activation by Lewis acid− transition metal ... with an electron deficient Lewis acid site to cat...
0 downloads 0 Views 2MB Size
Subscriber access provided by UNIV OF CALIFORNIA SAN DIEGO LIBRARIES

Article

General H2 Activation Modes for Lewis AcidTransition Metal Bifunctional Catalysts Yinwu Li, Cheng Hou, Jingxing Jiang, Zhihan Zhang, Cunyuan Zhao, Alister J. Page, and Zhuofeng Ke ACS Catal., Just Accepted Manuscript • DOI: 10.1021/acscatal.5b02395 • Publication Date (Web): 27 Jan 2016 Downloaded from http://pubs.acs.org on January 28, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

ACS Catalysis is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 9

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

General H2 Activation Modes for Lewis Acid‐Transition Metal Bifunc‐ tional Catalysts  Yinwu li,1 Cheng Hou,1 Jingxing Jiang,1 Zhihan Zhang,1 Cunyuan Zhao,1 Alister J. Page2 and Zhuofeng Ke1,* 1

MOE Key Laboratory of Bioinorganic and Synthetic Chemistry, School of Chemistry and Chemical Engineering, Sun Yatsen University, Guangzhou 510275, P. R. China 2 Newcastle Institute for Energy and Resources, The University of Newcastle, Callaghan 2308, NSW, Australia ABSTRACT: A general mechanism for H2 activation by Lewis acidtransition metal (LA-TM) bifunctional catalysts has been presented via DFT studies on a representative nickel borane system, (PhDPBPh)Ni. There are four typical H2 activation modes for LA-TM bifunctional catalysts: 1) the cis homolytic mode, 2) the trans homolytic mode, 3) the synergetic heterolytic mode, and 4) the dissociative heterolytic mode. The feature of each activa-

tion mode has been characterized by key transition state structures and natural bond orbital analysis. Among these four typical modes, (PhDPBPh)Ni catalyst most prefers the synergetic heterolytic mode (G‡ = 29.7 kcal/mol), however the cis homolytic mode cannot be totally disregarded (G‡ = 33.7 kcal/mol). In contrast, the trans homolytic mode and dissociative heterolytic mode are less feasible (G‡ = ~42 kcal/mol). The general mechanistic picture presented here is fundamentally important for the development and rational design of LA-TM catalysts in the future. KEYWORDS: H2 activation, Lewis acid, Lewis base, transition metal, nickel, mechanism, density functional theory, homolytic, synergetic, heterolytic, hydrogenation

1. INTRODUCTION  Hydrogen is one of the most important reproducible green energy resources.1 Catalytic H2 activation is widely used in the energy storage,1,2 hydrogenation,3 dehydrogenation,4 hydrodesulfurization,5 hydrodenitrogenation6 and many other transformations. The most challenging aspect of H2 activation is H-H bond cleavage, which arises because of the strength of the H-H  bond.7‐11

Figure 1. The development of TM catalysts, Lewis base-TM catalysts, and Lewis acid-TM catalysts for H2 activation. Transition metals (TMs) are particularly important in the development of homogenous catalysts for H2 activation, due to their ability for d-* back-donation to the H-H bond. The early well-known examplar TM catalyst is Wilkinson's catalyst (RhCl(PPh3)3), discovered in 1966, which is used for the acti-

vation of H2 and the hydrogenation of alkenes.12 Since then, a number of other TM catalysts with various ligands have been developed for the activation of hydrogen.7‐17 “Traditional” TM catalysts (since 1960s, Figure 1), generally subject H2 to oxidative addition onto their relatively electron rich TM center, via a homolytic cleavage mechanism (Figure 2).17,18 Later, another type of catalysts, Lewis base-transition metal (LB-TM) bifunctional catalysts were developed in 1980s (Figure 1).19‐21 LB-TM feature a ligand bearing a Lewis base site, which assists H2 activation.22 The archetypal LB-TM catalysts are the Noyori catalysts.23-25 The LB-TM bifunctional catalysts generally cleave H2 via a heterolytic mechanism, in which the metal center accepts the hydride and the assisting Lewis base attracts the proton (Figure 2). Main group elements with lone pair(s), such as nitrogen, oxygen, or sulfur, are usually utilized as the assisting Lewis base. LB-TMs based on many transition metals, such as Fe,26 Co,27 Rh,28,29 Ir,28,30,31 and Ru28,31-34 etc. have been developed in the past few decades. Very recently, a new family of Lewis acid TM (LA-TM) bifunctional systems has emerged. LA-TMs comprise an electron-rich late transition metal and a Lewis acid ligand35,36 with an electron deficient Lewis acid site to catalyze the activation

of H2. As far as we know, the TM and LB-TM systems have been widely studied for H2 activation, while the LATM catalytic systems did not emerge until 2010s (Figure 1).

ACS Paragon Plus Environment

ACS Catalysis

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

LB-TM Catalysts  

TM Catalysts  

M

D

homolytic













H H



H

H A

M M cis homolytic

A

M tr ans homolytic

A

M

mechanism? 

H

H H

A

M

heterolytic

H H





H H

H H



LA-TM Catalysts

M

synergetic heterolytic

A



H



M dissociative heterolytic

Proposed Four H 2 Activation Modes for LA-TM Catalysts

Figure 2. H2 activation mechanisms for TM catalysts, LB-TM catalysts, and LA-TM catalysts. In former studies, reversible binding of H2 in LA-TM system has been found.37,38 In 2011, Owen et al. reported that H2 could be activated by the rhodium borane catalytic system at 358.15 K.39 In 2012, Peter’s group synthetized a nickel borane complex, which catalyzed H2 activation at room temperature and low pressure, demonstrating the potential of LA-TM catalysts in H2 activation and hydrogenation.40 Erker and coworkers reported the hydrogen activation by an intermolecular zirconocene borane pair under mild conditions,41 and suggested a mechanism of heterolytic splitting of H2 similar to frustrated Lewis pairs.42 To date, a number of LA-TM catalytic systems have been developed, including Fe,43 Co,44 Ni,40 Ru,45 and Pt46 complexes. In these catalytic system, earth-abundant transition metals in first row47 such as Fe, Co, and Ni has attracted more interests, and boron is usually used as the Lewis acid site in ligands.48,49

Figure 3. H2 activation catalyzed by the (PhDPBPh)Ni system. In contrast to traditional TM and LB-TM catalysts, the H2 activation mechanism for LA-TM systems still remains to be fully understood. In a former study, Peter’s group found that H2 coordinates with (PhDPBPh)Ni forming a tetrahedral complex.50 Then the H2 cleavage leads to a dihydride product (Figure 3). Two potential mechanisms were proposed: H2 homolytic cleavage on the nickel center, and H2 heterolytic cleavage on the Ni-B bond. The heterolytic cleavage was suggested to be preferred due to its lower activation free energy.

Page 2 of 9

Sakaki’s group studied the different electronic effect51 between homolytic and heterolytic, and suggested that boron play an important role in the H2 cleavage in this LA-TM system.52 Herein, we investigate the mechanism of LA-TM catalyzed H2 activation using the representative (PhDPBPh)Ni nickel borane catalyst reported by Peter et al.,40 and present the first general mechanistic picture for LA-TM catalyzed H2 activation. We are particularly interested in the following pertinent questions concerning H2 activation with LA-TM catalysts: 1) Are there other types of H2 activation mechanisms for LA-TM catalysts? 2) Does the assistance of Lewis acid occur during H2 cleavage, or after cleavage to stabilize the hydride? 3) Does the LATM interaction remain during H2 cleavage? With these questions in mind, a comprehensive theoretical study was carried out on the (PhDPBPh)Ni system to unveil the general mechanism of H2 Activation by LA-TM catalysts. In this research, we presented four typical H2 activation modes for LA-TM catalysts (Figure 2): 1) The cis homolytic mode to a cis dihydride intermediate. 2) The trans homolytic mode to a trans dihydride species. 3) The synergetic heterolytic mode, where the H2 heterolytically cleaves onto the LA-TM bond. 4) The dissociative heterolytic mode, where the separate LA and TM activate H2 in a dissociative manner. These results are expected to provide a general mechanistic framework for future development and rational design of LATM catalysts.

2. COMPUTATIONAL DETAILS  All structures were optimized in gas-phase by density functional theory (DFT),53 using the M06-L54 functional with basis sets I (BSI, lanl2dz55 for metal atom and 6-31G (d, p) for nonmetal atoms) in the gas phase. It should be noted that the optimizations were performed in the gas-phase, with the ideal gas conditions at 0 K. The optimized structure is taken as a single molecular cluster, which may be different from the reaction conditions in a solution. Frequency analysis calculations for optimized structures were performed to characterize the structures to be minima (no imaginary frequency) or transition states (one imaginary frequency). IRC calculations were carried out to confirm the connection between two minima for each transition state. Based on M06-L/BSI optimized geometries, the energy results were further refined by calculating the single point energy at the M06-L/BSII level of theory (BSII designates SDD56 for metal atom and 6-311++G**57 for nonmetal atoms). The bulky solvation effect of benzene (= 2.3) were simulated by SMD56 continuum solvent mode at the M06-L/BSII level of theory. The standard Gibbs free energy in solution is calculate by the equation (1), °



°





°

ln 24.5

(1)

Eo(gas)

Where the is the gas-phase electronic energy, ZPE is the zero-point energy, Gsolv is the solvation free energy in benzene, and Go298K is the additional thermal contribution to Gibbs free energy at 298.15 K and 1 atm. The final term accounts for the free energy change from an ideal gas of 1 atm (24.5 L,

2

ACS Paragon Plus Environment

Page 3 of 9

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

298.15 K) to 1M solution. All natural bond orbital (NBO) analyses were performed using the NBO 3.0 package,58 and all electronic structure calculations were performed using Gaussian 09 program.59 The 3D optimized structures were displayed by CYLview visualization program.60   

3. RESULTS AND DISCUSSION  H2 cleavage initiates from intermediate (PhDPBPh)Ni(H2) (Figure 4), where the H2 is activated by nickel center through -coordination. This (PhDPBPh)Ni(H2) complex is located to be a tetrahedral structure, in agreement with previous study.50 The boron in (PhDPBPh)Ni(H2) is sp3 hybridized, indicating that the d10 Ni0 center contributes a pair of d electrons to boron’s empty p orbital, as expected for a typical LA-TM system. As shown in Figure 4, the H2  dative interaction and the d(Ni)σ*(H2) back donation both play an important role in H2 activation. The H–H bond length has been elongated to 0.824 Å in the (PhDPBPh)Ni(H2) complex, which is the evidence for d-σ* back donation.

Ph

Ph

Figure 4. Optimized ( DPB )Ni(H2) structure (left);  dative interaction and nickel d-σ* back donation (right).

Figure 5. Mechanism of the out-of-plane (cis) homolytic mode for the LA-TM catalyzed H2 activation (PPh2 is simplified as P for clarity; M = nickel).

Cis Homolytic Mode. The mechanism of the cis homolytic mode of LA-TM catalyzed H2 activation is depicted in Figure 5. This mechanism involves two steps: 1)the H2 cleavage and 2) the hydride isomerized to final product PC. In the first step, as shown in Figure 5, H2 is cleaved via a pyramidal transition

state 1-TS, which exhibits out-of-plane cis cleavage character, i.e. the boron atom, is out of the H2 cleavage plane (M-H-H). The Hb-Ni-Ha-B and P2-Ni-Ha-B dihedral angles in 1-TS1 are 93° and 95°, respectively (Figure 6). Thus, boron severs as a Lewis acid to stabilize 1-TS1 by interacting with the axial filled dz2 orbital of nickel center. The calculated activation free energy of this step is 27.1 kcal/mol. IRC analysis found that this H2 cleavage via 1-TS1 leads to the cis dihydride intermediate 1-IM, which has a pseudo pyramidal structure (Figure 6). In 1-IM, the boron atom bridges to Ni-Hb Bond, where the NiHa distance is 1.473 Å, and the B-Ha and B-Ni distances are 1.363 Å, 2.397 Å, respectively. Such pyramidal intermediates, in which the boron Lewis acid binds to the M-H bond and the M-(μ2-H)-B orients in the metal dz2 orientation, have been observed previously in related systems.61-63 This cis dihydride NiII intermediate is less stable than RC by 27.4 kcal/mol. In the second step, hydride Ha transfers trans to Hb, via the isomerization transition state 1-TS2, resulting in the square planar product PC. The activation free energy for this process is 33.7 kcal/mol relative to RC. In PC, Ha is trans to Hb with the Hb-Ni-Ha = 166°. The Ni-B distance is 2.402 Å, the NiHa is 1.468 Å, and the Ni-Hb is 1.637 Å. These values indicate that Ha moves closer to boron atom and bridges between B and Ni, forming a Ni-(μ2-H)-B structure. Compared with the Ni-Ha bond, the Ni-Hb bond in Ni-(μ2-H)-B is shorter, which is in accord with classical Ni-(μ2-H)-B structures.40 The out-of-plane H2 cleavage observed in the cis homolytic mode can be well illustrated by NBO analysis. In 1-TS1, the Ha-Hb distance is 1.821 Å, and so this bond is essentially broken. Furthermore, the two localized Ni-H σ bonds are nearly totally formed, as can be seen in the orbital occupancies in Figure 7. NBO analysis along the reaction coordinates provides more details for the changes in orbitals in this cis homolytic H2 activation process, as depicted in Figure 7. Prior to H2 cleavage, the dx2-y2(Ni) and dxy(Ni) orbitals contribute prominently to H2 activation. The dxy(Ni) orbital occupancy drops from 1.87 to 1.83 (see SI), before increasing to 1.94, due to back-donation to the σ*(H2) orbital. More importantly, the dx2y2(Ni) occupancy drops continuously from 1.67 to 1.58, until it is eventually oxidized in generating the σ(Ni-H) bonds. Meanwhile the σ(H2) bond vanishes at the transition state region. The H2 cleavage driven by the dx2-y2(Ni) and dxy(Ni) orbitals is indicative of a typical homolytic process, where the Ni center is oxidized from Ni0 to NiII and the boron stays out of the xy plane. In the later stage of H2 cleavage, the Ni-H orbital interacts with boron’s empty p* orbital, generating a Ni-(μ2-H)-B structure. In the subsequent isomerization step, Ha moves to the trans position of Hb with the assistance of phosphorus. As shown in Figure 7, the σ(P-H) bond generates temporarily during isomerization, simultaneously with an increase in the dx2-y2(Ni) occupancy. In the late stage of isomerization, the Ni-Hb bond forms again as the σ(P-H) orbital becomes unoccupied. The cis homolytic mode induces significant changes in coordination geometry along the H2 activation pathway, from tetrahedral RC to pyramidal transition state 1-TS1, and finally to a square planar PC. Therefore, a cis homolytic mode may be prefered for LA-TM catalysts with structural flexibility.

3

ACS Paragon Plus Environment

ACS Catalysis

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 9

Figure 6. Key structures for the mechanism of cis homolytic mode (each phenyl group in phosphines is visualized as one carbon atom for clarity).

ligand field may lead to large deformation energy due to the existence of filled anti-bonding orbital (dz2) in the transition state. We cannot totally exclude that in some LA-TM cases, without a filled anti-bonding dz2 orbital, the trans homolytic mode may be possible, and the functionality of the Lewis acid is primarily to help stabilize the hydride in the product.

Figure 7. NBO analysis along the cis homolytic reaction pathway of H2 activation.

Trans Homolytic Mode. The mechanism of the trans homolytic mode for H2 activation is depicted in Figure 8. In this mode, H2 is cleaved via a pyramidal transition state 2-TS, in which the boron atom lies in the Ni-H-H cleavage plane. The dihedral angle Ha-Ni-Hb-B is 174°, as shown in Figure 8. The dihedral angle P1-Ha-P2-B is -8.6°, indicating that Ha lies in the Ni-P1-P2 plane, while Hb is located axial to the NiP1-P2 plane during the H2 cleavage. In 2-TS, the Ha-Hb distance is 1.964 Å, whereas the Hb-B distance is 3.122 Å, suggesting a homolytic H2 cleavage without the direct cooperation of boron. However, the boron atom remains interacting with metal center datively (dNi-B = 2.272 Å). This mode that directly leads to the trans dihydride product PC, has a high activation free energy of 41.5 kcal/mol. Therefore, this mode is unlikely to be feasible for the (PhDPBPh)Ni(H2) system, which is in good agreement with a previous study of the (PhDPBPh)Ni system.50 NBO analysis provides further insight into the in-plane character of the trans homolytic mode (Figure 9). From RC to 2-TS, dxz(Ni) orbital occupancy drops from 1.87 to 1.67, while the dz2(Ni) orbital occupancy decreases from 1.67 to 1.54. This indicates that dxz(Ni) and dz2(Ni) orbitals are involved in the H2 cleavage to generate new σ(Ni-H) bonds. Figure 9 also shows that the σ(H2) occupancy drops from 1.77 to 0 and σ(Ni-Hb) appears before 2-TS. dz2(Ni) orbital occupancy disappears after 2-TS, and this is accompanied with the formation of the σ(Ni-Ha) bond. The interchange of σ(Ni-Ha) and σ(B-H) bonds, and the disappearance of p*(B) orbital occupancy indicate the formation of the Ni-(μ2-H)-B bridge structure observed in the product PC. There is a significant change in the coordination geometry during the trans homolytic mode. The change of

Figure 8. Mechanism and 2-TS structure of the trans homolytic mode for the LA-TM catalyzed H2 activation (P denotes PPh2 for clarity; M = nickel).

Figure 9. NBO analysis along the trans homolytic reaction pathway of H2 activation.

Synergetic Heterolytic Mode: The transition state (3-TS) for synergetic heterolytic H2 cleavage is shown in Figure 10. This transition state exhibits a quasi-tetrahedral structure, with a B-Ni-Ha-Hb dihedral angle of 163°. This is similar to the in-plane character of homolytic transition state 2-TS. However,

4

ACS Paragon Plus Environment

Page 5 of 9

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

in the heterolytic 3-TS, the HaHb distance is 1.228 Å, indicating a relatively early transition state compared to 1-TS1 and 2TS. The HbB distance in 3-TS is 1.505 Å, shorter than that in 2-TS (3.122 Å). This means that the boron Lewis acid interacts synergetically with the Hb atom during H2 cleavage. The activation free energy of this mode is 29.7 kcal/mol, which is comparable to previous results (∆G‡ = 28.9 kcal/mol) for the (PhDPBiPr)Ni model.50

Figure 10. Mechanism and key structure for the synergetic heterolytic mode of LA-TM catalyzed H2 activation (P denotes PPh2 for clarity; M = nickel).

boron atom plays a very important role in the H2 cleavage. On the other hand, the dz2(Ni) occupancy drops from 1.67 to zero at 3-TS. As the dz2(Ni) orbital becomes unoccupied, another σ(Ni-H) bond is formed. Boron and nickel thus act synergeticallyin this H2 activation mode, leading to a tetrahedral transition state in which H2 cleavage is heterolytic, with the boron associated with both the nickel center and the Hb atom. It is reasonable to anticipate that this synergetic heterolytic mode is more likely to be preferred, if the LA-TM catalyst has an appropriate LATM distance and the penalty in the weakening of LA-TM dative bond is insignificant. Dissociative Heterolytic Mode. This dissociative heterolytic H2 activation involves two steps (Figure 12): 1) the insertion of H2 into the B-Ni bond; and 2) the H2 heterolytic cleavage by the dissociative boron atom and nickel center. In the first step, H2 inserts into the B-Ni bond via transition state 4-TS1 to generate tetrahedral intermediate 4-IM, which is saturated with a solvent molecule, benzene. As show in Figure 13, H2 is activated in 4-IM by the dissociative boron and nickel centers (dHa-Hb = 0.839 Å), which facilitates the subsequent H2 cleavage. As compared to RC, the Ni-B distance increases from 2.210 to 3.557 Å in 4-IM. The activation free energy of the first step is 34.7 kcal/mol. This mechanism is similar to that observed with frustrated Lewis pairs catalysts.2,11,42,64 In the second step, H2 is cleaved by dissociative Ni and B centers via transition state 4-TS2. In 4-TS2, dH-H is 1.189 Å, and the NiH and BH distances are 1.483 and 1.276 Å, respectively. The Ha-Ni-B-Hb dihedral angle is -8°, indicating that H2 cleavage occurs in the Ni-H-H-B plane. The activation free energy for the dissociative heterolytic H2 cleavage is 42.2 kcal/mol.

Figure 11. NBO analysis along the synergetic heterolytic reaction pathway of H2 activation. The distinctive synergetic character for this heterolytic H2 activation mode can be revealed by NBO analysis (Figure 11). From RC to 3-TS, σ(Ha-Hb) occupancy drops slightly from 1.77 to 1.64 before disappearing at 3-TS. At the same time, the occupation of the LP(Hb) orbital increases to ~ 1, which clearly indicates heterolytic H2 cleavage. After 3-TS, the LP(Hb) orbital becomes unoccupied, because the LP(Hb) interacts with the LP*(B) to generate the Ni-(μ2-H)-B bridge bond in PC. Thus the assistance of boron atom plays a very important role in the H2 cleavage. After 3-TS, the LP(Hb) orbital becomes unoccupied, and this is accompanied by the LP(Hb) interacting with the LP*(B) orbital. This interaction generates the Ni-(μ2-H)-B bridge bond in PC, wherein the assistance of

Figure 12. Mechanism of the dissociative heterolytic mode for the LA-TM catalyzed H2 activation (P denotes PPh2 for clarity;

M = nickel).

5

ACS Paragon Plus Environment

ACS Catalysis

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 9

Figure 13. Key structures for the mechanism of dissociative heterolytic mode (each phenyl group in phosphines is visualized as one carbon atom for clarity).

Figure 14. NBO analysis along the dissociative heterolytic reaction pathway of H2 activation.

This dissociative heterolytic mode provides an alternative avenue for LA-TM catalyzed H2 activation. NBO analysis was performed to gain further insight into this interesting mode (Figure 14). In the insert step, σ(Ha-Hb) drops from 1.76 to 1.62 to generate the intermediate 4-IM. NBO analysis shows that σ(Ha-Hb) is activated by LP*(B) in the 4-IM,while σ*(HaHb) is activated by dz2(Ni). In the H2 cleavage step, σ(Ha-Hb) drops from 1.63 to 1.58 before ultimately disappearing to generate a σ(B-Hb) bond. This is accompanied by the LP*(B) orbital becoming unoccupies. Concomitantly, dz2(Ni) drops from 1.73 to 1.58 before being oxidized to generate a σ(Ni-Ha) bond. H2 activation by this mode proceeds via a dissociative LA and TM structure, and therefore requires the LA-TM dative interaction being overcome. This could potentially be achieved by using a frustrated LA-TM catalyst, or by weakening the LATM directly. Overview of the Mechanism for (PhDPBPh)Ni Catalyzed H2 Activation. Having presented detailed mechanistic aspects of these four general H2 activation modes for LA-TM catalysts, we now discuss their relative feasibility for the [(PhDPBPh)Ni] system. The potential free energy surfaces for all H2 activation modes are summarized in Figure 15.

Figure 15. Free energy profiles of the four H2 activation modes catalyzed by the [(PhDPBPh)Ni] complex (P denotes PPh2 for clarity).

6

ACS Paragon Plus Environment

Page 7 of 9

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

Although these four types of H2 activation modes are, in general, possible LA-TM catalysts, the trans homolytic mode and dissociative heterolytic mode are less feasible for the [(PhDPBPh)Ni] system. As shown in Figure 15, Both modes have an activation free energy of around 42 kcal/mol relative to RC (2-TS, G‡ = 41.5 kcal/mol; 4-TS2, G‡ = 42.2 kcal/mol). In the case of the trans homolytic mode, the transition state has high activation free energy arises for two reasons. 1)The electronic configuration of metal center changes from d10 to d8 during H2 cleavage. The anti-bonding orbital dz2(Ni) is filled in pyramidal transition state 2-TS (Figure 8), and this leads to a large deformation energy for structure changing from tetrahedron to pyramid.65 2) The boron atom doesn’t assist the H2 cleavage (dB-H distance in 2-TS is too long). As for the dissociative heterolytic mode, it is less feasible because it requires overcoming the boron-nickel dative interaction, which is expected to be very strong in the [(PhDPBPh)Ni] system. This is evidenced by the fact thatthe intermediate 4-IM is already higher than RC by 34.5 kcal/mol after the dissociation of boron-nickel bond. The synergetic heterolytic mode and the cis homolytic mode are more favored for (PhDPBPh)Ni catalyzed H2 activation. The respective rate-determining activation free energies, relative to RC, are 33.7 and 29.7 kcal/mol (Figure 15). The synergetic heterolytic mode is relative more preferred by 4.0 kcal/mol. However, the cis homolytic mode is still kinetically competitive, with a lower H2 cleavage transition state (1-TS1, G‡ = 27.1 kcal/mol) than that of the synergetic heterolytic mode (3TS, G‡ = 29.7 kcal/mol). Due to the thermodynamic instablility of the cis dihydride intermediate 1-IM, its subsequent isomerization to the trans dihydride product PC requires overcoming slightly higher activation free energy (1-TS2, G‡ = 33.7 kcal/mol). Therefore, the synergetic heterolytic mode is the most plausible mechanism for (PhDPBPh)Ni catalyzed H2 activation, while the cis homolytic mode cannot be totally excluded. The synergetic heterolytic mode is favored because the H2 cleavage can be cooperatively assisted by Lewis acid (boron) and transition metal (nickel). Furthermore, the (PhDPBPh)Ni complex has an insignificant change in coordination geometry during the synergetic heterolytic H2 cleavage. Of these two homolytic activation modes, the cis homolytic mode is more favorable than the trans homolytic mode on account of electronic factors. Firstly, there is no filled anti-bonding d(Ni) orbitals in 1-TS1 (Figure 6), which is a very late transition state that can be regarded as a NiII species. Evidence can be seen in Figure 7 that the dx2 -y2 (Ni) is totally oxidized to be an empty anti-bonding orbital. Secondly, the filled dz2(Ni) orbital can be stabilized by the empty boron p orbital.

4. CONCLUSION  In summary, a general H2 activation mechanism for the newly emerging LA-TM bifunctional catalysts has been elucidated. There are four typical H2 activation modes: (1) the cis homolytic mode through a cis dihydride intermediate; (2) the trans homolytic mode where H2 cleaves in the axial position to trans dihydride product; (3) the synergetic heterolytic mode, in which H2 cleaves cooperatively assisted by both the Lewis

acid and transition metal; and (4) the dissociative heterolytic mode, in which the Lewis acid and transition metal assist H2 cleavage in dissociative manner. The feature of each activation mode has been generalized by key transition state and natural bond orbital analysis. For the (PhDPBPh)Ni system, the calculated activation free energy is 33.7 kcal/mol for the cis homolytic mode, 41.5 kcal/mol for the trans homolytic mode, 29.7 kcal/mol for the synergetic heterolytic mode, and 42.2 kcal/mol for the dissociative heterolytic mode, respectively. Thus, the synergetic heterolytic mode is found to be the most plausible mechanism for (PhDPBPh)Ni catalyst, although the cis homolytic mode cannot be totally excluded. In contrast, the trans homolytic mode and dissociative heterolytic mode are less feasible for the [(PhDPBPh)Ni] system. Considering the structural diversity of LA-TM catalysts, it can be expected that other LA-TM systems may differ from [(PhDPBPh)Ni] and activate H2 via different scenarios. These results should provide a general mechanistic framework for the development and rational design of LA-TM catalysts in the future.

ASSOCIATIVE CONTENT  Supporting Information. Complete ref. 59; details of orbital occupancy calculations, and Cartesian coordinates. This material is available free of charge via the Internet at http://pubs.acs.org.

AUTHOR INFORMATION  Corresponding Author  *E-mail: [email protected] Notes  The authors declare no competing financial interest.

ACKNOWLEDGMENTS  This work was supported by the NSFC Fundation (Grants 21203256 and 21473261), the Guangdong Natural Science Funds for Distinguished Young Scholar (No. 2015A030306027). Computing facilities were supported in part by the high performance grid computing platform of Sun Yat-sen University, the Guangdong Province Key Laboratory of Computational Science and the Guangdong Province Computational Science Innovative Research Team. We thank Prof. Jonas C. Peters for very helpful discussions.

REFERENCES  (1) Satyapal, S.; Petrovic, J.; Read, C.; Thomas, G.; Ordaz, G. Catal. Today 2007, 120, 246256. (2) Hamilton, C. W.; Baker, R. T.; Staubitz, A.; Manners, I. Chem. Soc. Rev. 2009, 38, 279293. (3) Clapham, S. E.; Hadzovic, A.; Morris, R. H. Coord. Chem. Rev. 2004, 248, 22012237. (4) Keaton, R. J.; Blacquiere, J. M.; Baker, R. T. J. Am. Chem. Soc. 2007, 129, 18441845. (5) Sattler, A.; Parkin, G. J. Am. Chem. Soc. 2011, 133, 37483751. (6) Suresh, C.; Santhanaraj, D.; Gurulakshmi, M.; Deepa, G.; Seivaraj, M.; Rekha, N. R. S.; Shanthi, K. Acs Catal 2012, 2, 127134.

7

ACS Paragon Plus Environment

ACS Catalysis

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(7) Torrent, M.; Sola, M.; Frenking, G. Chem. Rev. 2000, 100, 439493. (8) Niu, S. Q.; Hall, M. B. Chem. Rev. (Washington, DC, U. S.) 2000, 100, 353405. (9) He, T.; Tsvetkov, N. P.; Andino, J. G.; Gao, X.; Fullmer, B. C.; Caulton, K. G. J. Am. Chem. Soc. 2010, 132, 910911. (10) Kimura, T.; Koiso, N.; Ishiwata, K.; Kuwata, S.; Ikariya, T. J. Am. Chem. Soc. 2011, 133, 88808883. (11) Camaioni, D. M.; Ginovska-Pangovska, B.; Schenter, G. K.; Kathmann, S. M.; Autrey, T. J. Phys. Chem. A 2012, 116, 72287237. (12) Osborn, J. A.; Jardine, F. H.; Young, J. F.; Wilkinson, G. J. Chem. Soc. A 1966, 17111732. (13) Maseras, F.; Lledos, A.; Clot, E.; Eisenstein, O. Chem. Rev. (Washington, DC, U. S.) 2000, 100, 601636. (14) Heinekey, D. M.; Lledos, A.; Lluch, J. M. Chem. Soc. Rev. 2004, 33, 175182. (15) McGrady, G. S.; Guilera, G. Chem. Soc. Rev. 2003, 32, 383392. (16) Andrews, L. Chem. Soc. Rev. 2004, 33, 123132. (17) Kubas, G. J. Chem. Rev. 2007, 107, 41524205. (18) Menzel, M.; Winkler, H. J.; Ablelom, T.; Steiner, D.; Fau, S.; Frenking, G.; Massa, W.; Berndt, A. Angew. Chem., Int. Ed. 1995, 34, 13401343. (19) Karvembu, R.; Prabhakaran, R.; Natarajan, K. Coord. Chem. Rev. 2005, 249, 911918. (20) Blum, Y.; Czarkie, D.; Rahamim, Y.; Shvo, Y. Organometallics 1985, 4, 14591461. (21) Shvo, Y.; Czarkie, D. J. Organomet. Chem. 1986, 315, C25C28. (22) Kuzu, I.; Krummenacher, I.; Meyer, J.; Armbruster, F.; Breher, F. Dalton Trans. 2008, 58365865. (23) Haack, K.-J.; Hashiguchi, S.; Fujii, A.; Ikariya, T.; Noyori, R. Angew. Chem., Int. Ed. 1997, 36, 285288. (24) Hashiguchi, S.; Fujii, A.; Haack, K.-J.; Matsumura, K.; Ikariya, T.; Noyori, R. Angew. Chem., Int. Ed. 1997, 36, 288290. (25) Doucet, H.; Ohkuma, T.; Murata, K.; Yokozawa, T.; Kozawa, M.; Katayama, E.; England, A. F.; Ikariya, T.; Noyori, R. Angew. Chem., Int. Ed. 1998, 37, 17031707. (26) Song, H.; Kang, B.; Hong, S. H. Acs Catal. 2014, 4, 28892895. (27) Zhang, G. Q.; Hanson, S. K. Org. Lett. 2013, 15, 650653. (28) Wu, X. F.; Li, X. H.; Zanotti-Gerosa, A.; Pettman, A.; Liu, J. K.; Mills, A. J.; Xiao, J. L. Chem.-Eur. J. 2008, 14, 22092222. (29) Matharu, D. S.; Morris, D. J.; Clarkson, G. J.; Wills, M. Chem. Commun. 2006, 32323234. (30) Heiden, Z. M.; Rauchfuss, T. B. J. Am. Chem. Soc. 2009, 131, 35933600. (31) Lei, M.; Zhang, W. C.; Chen, Y.; Tang, Y. H. Organometallics 2010, 29, 543548. (32) Sandoval, C. A.; Ohkuma, T.; Muniz, K.; Noyori, R. J. Am. Chem. Soc. 2003, 125, 1349013503. (33) Abdur-Rashid, K.; Clapham, S. E.; Hadzovic, A.; Harvey, J. N.; Lough, A. J.; Morris, R. H. J. Am. Chem. Soc. 2002, 124, 1510415118. (34) Ikariya, T.; Blacker, A. J. Acc. Chem. Res. 2007, 40, 13001308. (35) Amgoune, A.; Bourissou, D. Chem. Commun. 2011, 47, 859871. (36) Braunschweig, H.; Dewhurst, R. D. Dalton Trans. 2011, 40, 549558. (37) Alcaraz, G.; Helmstedt, U.; Clot, E.; Vendier, L.; SaboEtienne, S. J. Am. Chem. Soc. 2008, 130, 1287812879. (38) Alcaraz, G.; Grellier, M.; Sabo-Etienne, S. Acc. Chem. Res. 2009, 42, 16401649. (39) Tsoureas, N.; Kuo, Y. Y.; Haddow, M. F.; Owen, G. R. Chem. Commun. 2011, 47, 484486.

Page 8 of 9

(40) Harman, W. H.; Peters, J. C. J. Am. Chem. Soc. 2012, 134, 50805082. (41) Podiyanachari, S. K.; Frohlich, R.; Daniliuc, C. G.; Petersen, J. L.; Muck-Lichtenfeld, C.; Kehr, G.; Erker, G. Angew. Chem., Int. Ed. 2012, 51, 88308833. (42) Stephan, D. W.; Erker, G. Angew. Chem., Int. Ed. 2010, 49, 4676. (43) Fong, H.; Moret, M. E.; Lee, Y.; Peters, J. C. Organometallics 2013, 32, 30533062. (44) Lin, T. P.; Peters, J. C. J. Am. Chem. Soc. 2013, 135, 1531015313. (45) Boone, M. P.; Stephan, D. W. J. Am. Chem. Soc. 2013, 135, 85088511. (46) Barnett, B. R.; Moore, C. E.; Rheingold, A. L.; Figueroa, J. S. J. Am. Chem. Soc. 2014, 136, 1026210265. (47) Kaye, S. S.; Long, J. R. J. Am. Chem. Soc. 2005, 127, 6506-6507. (48) Emslie, D. J. H.; Cowie, B. E.; Kolpin, K. B. Dalton Trans. 2012, 41, 1101-1117. (49) Dang, L.; Lin, Z.; Marder, T. B. Chem. Commun. 2009, 39873995. (50) Harman, W. H.; Lin, T. P.; Peters, J. C. Angew. Chem., Int. Ed. 2014, 53, 10811086. (51) Ohnishi, Y. Y.; Nakao, Y.; Sato, H.; Sakaki, S. J. Phys. Chem. A 2007, 111, 79157924. (52) Zeng, G.; Sakaki, S. Inorg. Chem. 2013, 52, 28442853. (53) Hohenberg, P.; Kohn, W. Phys. Rev. 1964, 136, B864B871. (54) Zhao, Y.; Truhlar, D. G. Theor. Chem. Acc. 2008, 120, 215241. (55) Hay, P. J.; Wadt, W. R. J. Chem. Phys. 1985, 82, 299310. (56) Marenich, A. V.; Cramer, C. J.; Truhlar, D. G. J. Phys. Chem. B 2009, 113, 45384543. (57) Rassolov, V. A.; Pople, J. A.; Ratner, M. A.; Windus, T. L. J. Chem. Phys. 1998, 109, 12231229. (58) NBO Version 3.1, Glendening, E. D.; Reed, A. E.; Carpenter, J. E.; Weinhold, F. (59) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; Nakatsuji, H.; Caricato, M.; Li, X.; Hratchian, H. P.; Izmaylov, A. F.; Bloino, J.; Zheng, G.; Sonnenberg, J. L.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Montgomery, J. A.; Peralta, J. E.; Ogliaro, F.; Bearpark, M.; Heyd, J. J.; Brothers, E.; Kudin, K. N.; Staroverov, V. N.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Rega, N.; Millam, J. M.; Klene, M.; Knox, J. E.; Cross, J. B.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Zakrzewski, V. G.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Dapprich, S.; Daniels, A. D.; Farkas, Ö.; Foresman, J. B.; Ortiz, J. V.; Cioslowski, J.; Fox, D. J. Gaussian 2009, Revision D.01; Gaussian, Inc. Wallingford CT, 2009. (60) Legault, C. Y., CYLview, 1.0b; Université de Sherbrooke 2009, http://www.cylview.org. (61) Tsoureas, N.; Hope, R. F.; Haddow, M. F.; Owen, G. R. Eur. J. Inorg. Chem. 2011, 52335241. (62) Owen, G. R.; Tsoureas, N.; Hope, R. F.; Kuo, Y. Y.; Haddow, M. F. Dalton Trans. 2011, 40, 59065915. (63) Foreman, M. R. S.; Hill, A. F.; Tshabang, N.; White, A. J. P.; Williams, D. J. Organometallics 2003, 22, 55935596. (64) Stephan, D. W.; Erker, G. Angew. Chem., Int. Ed. 2015, 54, 64006441. (65) Hou, C.; Jiang, J.; Li, Y.; Zhang, Z.; Zhao, C.; Ke, Z. Dalton Trans. 2015, 44, 1657316585.

8

ACS Paragon Plus Environment

Page 9 of 9

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

TOC

9

ACS Paragon Plus Environment