DFT Calculations on the Mechanism of Transition-Metal-Catalyzed

Jul 29, 2016 - ... and §NYU−ECNU Center for Computational Chemistry at New York University Shanghai, East China Normal University, Shanghai 200062,...
0 downloads 0 Views 1MB Size
Subscriber access provided by SIMON FRASER UNIV BURNABY

Article

DFT Calculations on the Mechanism of Transition Metal Catalyzed Reaction of Diazo Compounds with Phenols: O-H Insertion versus C-H Insertion Yuan Liu, Zhoujie Luo, John Zenghui Zhang, and Fei Xia J. Phys. Chem. A, Just Accepted Manuscript • DOI: 10.1021/acs.jpca.6b05735 • Publication Date (Web): 29 Jul 2016 Downloaded from http://pubs.acs.org on August 2, 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.

The Journal of Physical Chemistry A 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 28

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

The Journal of Physical Chemistry

DFT Calculations on the Mechanism of Transition Metal Catalyzed Reaction of Diazo Compounds with Phenols: O-H Insertion versus C-H Insertion Yuan Liu,† Zhoujie Luo,‡ John Zenghui Zhang,†,§ Fei Xia*,‡,§ †

School of Physics and Materials Science, East China Normal University, ‡School of Chemistry and Molecular Engineering, East China Normal University, and §NYU-ECNU Center for Computational Chemistry at New York University Shanghai, East China Normal University, Shanghai 200062, China

ABSTRACT The reaction of diazo compounds with transition metal carbenes is an efficient way to achieve the functionalization of chemical bonds in organic molecules, especially for the C-H and O-H bonds. However, the selective mechanisms of C-H and O-H bond insertions by various metal carbenes such as Rh and Cu complexes are not quite clear. In this work, we performed a comprehensively theoretical investigation of the phenol C-H and O-H bonds inserted by Rh and Cu carbenes by using DFT calculations. The calculated results reveal that the nucleophilic additions of phenols to the Rh and Cu carbenes in the C-H bond insertions are the rate-determining steps of whole reactions, which are higher than the barriers in the O-H insertions. In the process of intramolecular [1,3]-H transfer, the Rh and Cu ligands in their carbenes tend to dissociate into solution rather than the intramolecular migration due to their weak metal-carbon bonds. A deeply theoretical analysis on the electronic structures of Rh, Cu and Au carbenes as well as their complexes elucidated their differences in the chemo-selectivity of C-H and O-H insertion products, which accords with the experimental observations well.

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

1. INTRODUCTION Small organic molecules such as phenols are the important natural resources in organic synthesis due to their wide availability and low prices. The specific functionalization of O-H or C-H bonds in organic molecules1-8 is a severe challenge for organic chemists, especially for the selective functionalization of inert C-H bonds. One way to functionalize the C-H and O-H bonds is to employ the metal carbenes generated from the mixture of diazo compounds and transtion metal catalysts.9-16 The generated metal carbenes possess high reactity toward the C-H and O-H bonds so that the methology of employing metal carbenes was rapidly developed these years. In most cases, the widely used transition metal catalysts include the copper,17-23 rhodium24-29 and gold30-34 complexes. The generated transition metal carbenes exhibit remarkable catalysis in activating the C-H and O-H bonds, which arouse the interests of theoretical researchers to explore their mechanisms. The previous theoretical studies35-39 mainly focus on the mechanisms of insertion reactions of C-H bonds in alkanes and O-H bonds in water or alcohol. One hand, for the mechanism of alkyl C-H bond insertion, Nakamura et al.35 firstly carried out a comprehensive study of alkyl C-H bonds catalyzed by dirhodium tetracarboxylate complexes and proposed that the C-H bond activation follows a concerted mechanism via a three-centered transition state. Afterwards, Braga et al.36 performed a further study of alkane C-H bond insertion by copper and silver complexes. Recently, we reported a detailed comparison of the mechanisms of aromatic C-H bonds of anisole and alkyl C-H bonds inserted by gold carbenes,40 which reveals that the aromatic C(sp2)-H bond activation follows a stepwise mechanism actually. Xie et al.41 studied the mechanism of Rhodium-catalyzed C-H functionalization of indole, they concluded the enol pathway is always favored for the phenyl-substituted carbene. On the other hand, Liang et al.37 explored the mechanism of O-H bond insertion of water molecules catalyzed by copper and rhodium complexes. They concluded that the O-H bond insertion catalyzed by copper complex follow the ligand-associated pathway, while the rhodium complex prefers the ligand-free pathway in the reaction process.

ACS Paragon Plus Environment

Page 2 of 28

Page 3 of 28

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

The Journal of Physical Chemistry

Xie et al.39 calculated the mechanism of O-H bonds in allyl alcohol inserted by rhodium carbenes and suggested that the O-H insertion products could be obtained through a [1,3]-proton shift in the enol intermediate or a [1,2]-proton shift in the oxonium ylide, which were regulated by the orientation of the ester groups in carbenes. The most of current experiments attempt to improve the efficiency of C-H and O-H bond insertions catalyzed by metal carbenes. In particular, how to enhance the chemo-selectivity toward C-H bond insertion remains a severe challenge for organic chemists. Very recently, Zhang42 and Shi’s group43 independently found that the triphenylphosphine oxide gold (Au-TPPO) carbenes have an unprecedented selectivity toward the C-H bond insertion in phenols. The prominent catalysis of gold catalysts stimulated us to perform a detailed theoretical and experimental study toward its selective mechanisms.44 Our recently published results reveal that the stability of enol intermediates play an important role in the selectivity. It is the first case in which the selective mechanisms of C-H and O-H bond insertions were detailedly elucidated by theoretical calculations.44 However, the mechanisms of selective C-H and O-H bond insertions catalyzed by other important transition metal Rh and Cu carbenes remain unclear, especially for the pivotal question how the effects of electronic structures of the coordinated transition metals affect the selectivity.20, 42, 44 In this work, we carried out a comprehensively theoretical investigation on the C-H and O-H bonds inserted by Rh and Cu carbenes. The reactions of phenols and methyl α-diazophenylacetate catalyzed by the catalysts with the ligands copper bisoxazoline (Cu-box), copper pyridine bisoxazoline (Cu-pybox) and dirhodium tetraacetate (Rh-TC) were chosen as the research targets, since the three complexes above were widely used in lab to generate the related carbenes for reacting with the diazo compounds.20, 27 By using the density functional theory(DFT) methods, we calculated the detailed pathways and free energy profiles of C-H and O-H bonds inserted by Rh and Cu carbenes. Based on a deeply theoretical analysis and comparison of electronic structures of key intermediates of Rh, Cu and

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Au carbenes, we elucidate the reason why the Rh, Cu and Au carbenes give rise to different C-H and O-H bond insertion products. 2. COMPUTATIONAL METHODS All the DFT calculations throughout this work were carried out using the Gaussian 09 software package.45 The geometric structures of intermediates and transition states were optimized using the M06 functional.46,47 The M06 functional has been demonstrated that it could give accurate descriptions for the intramolecular interactions in organometallic systems, especially for the noncovalent interactions between aromatic groups. The Lanl2dz basis set48,49 combined with the effective core potential was used to describe the heavy elements Rh, Cu and Au, and the 6-31G* basis set50 was utilized to describe the nonmetallic elements C, N, O, P and H, as the same as the computational levels that we used in the previous studies.40,44 The frequency analyses were performed on the optimized structures to verify that the intermediates are local minima and transition states are stationary points that possess only one imaginary frequency. The solvent effect of CH2Cl2 in the reactions was evaluated using the integral equation formalism model (IEFPCM)51 with a dielectric constant ε = 8.93, based on the structures obtained in gas phase. All the calculated energies of intermediates and transition states discussed here refer to the standard Gibbs free energies ∆Gsol in kcal/mol calculated at the temperature 298 Kelvin, including the solvation corrections evaluated from the IEFPCM51 model. The natural bond orbital (NBO)52 analysis and Wiberg bond order53 were performed at the same computational level. The contour images of molecular orbitals of key structures were plotted using the visualized package GaussView.54 The Cartesian coordinates of the intermediates and transition states along the calculated lowest energy pathways are provided in the section 1 of Supporting Information. 3. RESULTS AND DISCUSSION

ACS Paragon Plus Environment

Page 4 of 28

Page 5 of 28

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

The Journal of Physical Chemistry

3.1 Pathways of O-H Bond Insertion Catalyzed by Rh Carbenes The DFT calculations performed by Liang37 and Liu et al.40 demonstrate that the metal catalysts could react with diazo compounds to release dinitrogen molecules and form active metal carbenes as the precursors for the insertion reactions. The optimized structures of Rh-TC, Cu-pybox, Cu-box and Au-TPPO carbenes as well as the phenol are displayed in Figure 1. The transition metal carbenes and phenols are regarded as the reactants of the whole O-H insertion reactions, with their sum of free energies being 0.0 kcal/mol. Figure 2 shows the calculated reaction pathways and free energy profiles of the O-H bond insertion by the Rh-TC carbenes. The experiments performed by Zhang et al.55 have demonstrated that the attacking of the hydroxyl oxygen of benzyl alcohol to the carbene carbon leads to the species of oxonium ylide. By means of DFT calculations, we found that the oxonium ylide Int-o1 forms through the transition state TS-o1 with a modest barrier 14.4 kcal/mol. However, the energy of ylide Int-o1 is higher than the reactants by 13.8 kcal/mol and quite unstable in thermodynamics. There exist three possible reaction pathways37,44 connecting the ylide Int-o1 on the energy surfaces. One of pathways is that Int-o1 transforms to the ligand-free ylide Int-o2* by releasing the Rh-TC ligand into solution. To reach the product Pro-o1, the hydroxyl hydrogen in Int-o2* needs to transfer to the carbene carbon, through a water mediated TS-o2* with an endothermicity of 23.3 kcal/mol. The second pathway involves the formation of the metal-associated ylide Int-o2** with an endothermicity of 20.4 kcal/mol, where the Rh-TC ligand migrates to the aromatic ring in carbenes. The third pathway is the direct [1,2]-proton transfer in TS-o2**. However, the barriers of three pathways relative to Int-o1 above are too high and should be ruled out from the possible mechanisms. In our recent study of the O-H bonds inserted by Au-TPPO carbenes, we proposed that the hydroxyl hydrogen in ylide tends to transfer to the adjacent carboxyl oxygen and form enol complexes.44 Here, we optimized the enol complex Int-o2 as well as TS-o2 and found that Int-o2 is more stable than Int-o1 by 16.3 kcal/mol. Subsequently, the Int-o2 transforms to the enol Int-o3 by

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

liberating the Rh-TC ligand into solution instead of the migration to form Int-o3*. Finally, the [1,3]-H transfer occurs via TS-o3, which leads to the final stable product Pro-o1. In addition, in order to compare to the pathway from the carbonyl side shown in Figure 2, we have calculated the pathway of nucleophilic addition from the methoxy group side and the corresponding energy profiles are shown in Figure S1 of Supporting Information. Obviously, the lowest energy reaction pathway of that from the methoxy group side is less favorable in energetics than that from the carbonyl side. 3.2 Pathways of C-H Bond Insertion Catalyzed by Rh Carbenes The energy profiles of C-H insertion catalyzed Rh-TC carbene are shown in Figure 3. Since the para site of phenol is the most favorable site for the phenol C-H bond insertion, we only calculated the reaction pathway occurred at the para-position of phenol. Our recent theoretical study40 indicates that the first step of the aromatic C(sp2)-H insertion is the direct addition of the C(sp2)-C(sp2) atoms between the phenol and metal carbene. The barrier of C(sp2)-C(sp2) addition via the transition state TS-c1 is calculatedly 19.2 kcal/mol, which is obviously higher than 14.4 kcal/mol of the O-H bond nucleophilic attack in Figure 2. The addition intermediate Int-c1 is unstable in thermodynamics, with the energy 12.2 kcal/mol higher than the reactants. The [1,3]-migration intermediate Int-c2* can be hardly obtained because of its high endothermicity. Similar to the process of O-H insertion, the carbonyl oxygen atom is ready to abstract the proton to form the enol complex Int-c2, with a small barrier of 3.8 kcal/mol via TS-c2. The enol complex Int-c2 could be further stabilized by releasing the Rh-TC ligand into solution, leading to a free enol Int-c3. Through the TS-c3, the hydrogen transfer occurs in Int-c3 and the final product Pro-c1 is produced. The rate-determining step for the C-H bond insertion is the C(sp2)-C(sp2) addition with the barrier of 19.2 kcal/mol. It is higher than that of [1,3]-hydrogen transfer in the O-H insertion by 2 kcal/mol. Based on the kinetic viewpoint, the difference of 2 kcal/mol in free energy will leads to a large discrepancy in the reaction

ACS Paragon Plus Environment

Page 6 of 28

Page 7 of 28

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

The Journal of Physical Chemistry

rates. Therefore, the O-H insertion in phenol is more facile than the C-H insertion catalyzed by the Rh-TC carbene, which accords with the experimental results. 3.3 Pathways of O-H and C-H Bond Insertion Catalyzed by Cu-Carbenes Figure 4 shows the calculated lowest energy profiles of phenol O-H and C-H bond insertions catalyzed by Cu-pybox carbenes, denoted by the red and black curves, respectively. For the O-H insertion, the O-H nucleophilic attack of phenol to the Cu-pybox carbene is still facile, requiring an activation free energy of 16.1 kcal/mol via transition state TS-o4. Through TS-o5, the enol complex Int-o5 is generated and further transforms to a more stable enol Int-o3 by liberating the Cu-pybox ligand into solution. Finally, the product Pro-o1 is formed through a water-assisted transition state TS-o3. The crucial step in the O-H bond insertion is the [1,3]-hydrogen transfer and the total energy barrier is 17.2 kcal/mol. For the C-H insertion, the reaction process is similar to that in Figure 3. The rate-determining step is the first step of nucleophilic attack via TS-c4, with the barrier 19.3 kcal/mol. Figure 5 shows the lowest energy pathways of C-H and O-H bond insertions catalyzed by Cu-box carbenes. As for the first steps of nucleophilic attack of C-H and O-H bond insertions, the energy difference between the transition states TS-c6 and TS-o6 is 3.5 kcal/mol. For the C-H insertion, it is still the rate-determining step in the whole reaction. The crucial barriers of [1,3]-H transfer for TS-o3 and TS-c3 in the O-H and C-H bond insertions are almost the same, 17.2 and 17.5 kcal/mol, respectively. Nevertheless, the [1,3]-H transfer is still the rate-determining step for the O-H insertion, lower than the crucial step of C-H insertion by 3 kcal/mol. Thus, the O-H bond insertion by Cu-box carbenes is also more favorable than C-H insertion in kinetics. The calculated results explain the experimental selectivity of C-H and O-H bonds of phenol well. 3.4 Nucleophilic Addition of C-H Bond Insertion

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Based on the pathways calculated for the Rh-TC, Cu-pybox and Cu-box carbenes shown above, we found that the first addition steps in the C-H bond insertions always possess higher barriers, with 19.2, 19.3 and 20.2 kcal/mol for the three carbenes respectively. Why the first steps in the C-H insertions have such high barriers? The answer could be obtained by analyzing the electronic structures of transition metal carbenes. Figure 1(b) shows that the Wiberg bond order of Rh-C bond in Rh-TC carbene is 0.83, which is larger than the values of other three metal-carbon bonds in Figure 1(c)-(e). The molecular orbital of Rh-TC carbene in Figure.6(a) shows that a typical π orbital of Rh-TC carbene is formed by the dxy orbital of rhodium and the px orbital of carbon, and meanwhile Figure 6(b) shows the sdx2 hybrid orbital of rhodium overlaps with the sp2 orbital of carbon and form the σ bond between them. Therefore, we conclude that the Rh-C bond in the Rh-TC carbene is actually a double bond. Figure 6(c) and Figure 6(d) show the occupied orbitals of the Cu-box and Au-TPPO carbenes. The Cu-C and Au-C bonds exhibit the similar orbital composition, which consist of sd hybrid orbital of gold and the sp2 hybrid orbital of carbon. Due to the lack of the π conjugation, the Cu-C and Au-C bonds are actually single bonds that are weaker than the Rh-C double bond. Even though the bond orders of Cu-C bonds in Cu-pybox carbene and Cu-box carbene are lower than that in Rh-TC carbene, we noticed that the electronic effects of the coordinated ligands have an important influence on the barriers of addition. The Cu-pybox and Cu-box ligands are typical electron-donating groups in experiments, while the Au-TPPO ligand is an electron-withdrawing one. The electron-donating effect makes the carbene carbons more electronegative so that the nucleophilic addition becomes more difficult and vice versa. It explains the reason why the addition in the C-H insertion catalyzed by Au-TPPO carbenes has a lower barrier of 15.6 kcal/mol,44 lower than that of Cu-pybox and Cu-box carbenes by 3.7 and 4.6 kcal/mol, respectively. 3.5 Analysis of Key Enol Intermediates

ACS Paragon Plus Environment

Page 8 of 28

Page 9 of 28

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

The Journal of Physical Chemistry

Apart from the nucleophilic additions in the C-H insertions, the steps of [1,3]-H transfers also play an important role in the pathways of O-H bond insertions. The [1,3]-H transfers in the O-H insertions involve the crucial precursors, namely, the free enols which were generated from the corresponding enol complexes. Whether it is easy for the enol complexes to liberate the attached metal ligands has an important influence on the subsequent barriers of [1,3]-H transfers. Moreover, it also determines whether the enol complexes prefer a ligand-associated mechanism or ligand-free along the reaction pathways. Thus, a deeply theoretical analysis upon the stability of enol complexes is necessary. We evaluated the NBO charges and the Wiberg bond orders of four enol complexes, Int-o2, Int-o5, Int-o7 and Au-enol, as shown in Figure 7. The bond orders of the metal-carbon bonds in Int-o2, Int-o5, Int-o7 and Au-enol are 0.26, 0.15, 0.21 and 0.38 respectively, which reveals that Au-enol has the strongest Au-C bond among the enol complexes. The NBO analysis shows that the charge of the rhodium atom in the Rh-TC ligand is 0.68, while the charge of the rhodium atom in the complex reduced to 0.50. The reduction of charge with the 0.18 unit indicates that the rhodium atom accepts the electron donation from the carbon atom through the orbital interactions between them. Similarly, the NBO charges in Cu-pybox, Cu-box and Au-TPPO decrease by 0.04, 0.15 and 0.23, respectively. Obviously, the Au-C bond has the largest bond order value 0.38, which means that the stronger interaction exits between the Au and C atoms. This explains the reason why the Au-TPPO ligand prefers the migration to the liberation into the solution in the C-H and O-H bond insertions.44 In contrast, the Rh and Cu carbenes tend to liberate the ligands into solution due to the weak interactions between the metal and carbon atoms in the enol complexes. 3.6 Composition of Metal-Carbon Bonds In a bonding analysis further, Figure 8 shows a schematic comparison of the orbital composition of metal-carbon bonds in Int-o2, Int-o7 and Au-enol. Figure 8(a)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

indicates that there exists a weak σ interaction between the dxy orbital of rhodium and the px orbital of the carbene carbon. Compared to the Rh-C double bond in Figure 6(a), the Rh-C interaction in Int-o2 has decreased dramatically since a NBO analysis reveals that the π conjugation between them was broken. It means that the dissociation of the Rh-TC ligand from Int-o2 is a quite facile, which explains the reason why Int-o2 has a low barrier for dissociation. Figure 8(b) and Figure 8(c) display the HOMO orbitals of Int-o7 and Au-enol, respectively. Although the Cu and Au atoms in the carbenes adopt a pattern of sd hybridization to overlap with the sp2 orbitals of carbon atoms, it is evident that the overlap in Figure 8(c) is more effective than that in Figure 8(b). The reason lies in the fact that the 6s orbital in Au atom has larger orbital radius so that it improves the overlaps of atomic orbitals between the gold an carbon atoms. In our previous study,44 the dissociation of Au-TPPO ligand of Au-enol in Figure.7(g) into solution needs an endothermic energy of 7.1kcal/mol. Thus, the Au-TPPO ligand prefers to migrate to the aromatic ring in Au-enol and form new interactions to compensate for the energy loss in the dissociation. 4. CONCLUTIONS In this work, we performed a comprehensively theoretical investigation on the mechanisms of selective functionalization of phenol O-H and C-H bond insertions catalyzed by the Rh-TC, Cu-pybox and Cu-box carbenes by using the DFT methods. By comparison of the mechanisms of Rh-TC, Cu-pybox and Cu-box carbenes to that of Au-TPPO carbene, we found that two remarkable differences exist in the mechanisms of selectivity between the Au-TPPO carbene and others. For the C-H and O-H bond insertions catalyzed by Rh-TC, Cu-pybox and Cu-box carbenes, the first steps of nucleophilic additions of phenols to carbenes in the C-H insertions have considerably high barriers with the values of 19.2, 19.3 and 20.2 kcal/mol, respectively. They are higher than the barriers of the nucleophilic additions in the O-H insertions as well as that of [1,3]-H transfers in the C-H and O-H insertions. Therefore, the nucleophilic additions in the C-H bond insertions are

ACS Paragon Plus Environment

Page 10 of 28

Page 11 of 28

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

The Journal of Physical Chemistry

rate-determining steps and unfavorable in kinetics. The calculated pathways account for the experimental observations why the major O-H bond insertion products catalyzed by Rh-TC, Cu-pybox and Cu-box carbenes have been obtained. Secondly, the Rh-TC, Cu-pybox and Cu-box ligands tend to dissociate from the enol complexes to from free enols which are followed by the [1,3]-H transfers. The Au-TPPO ligand in Au-carbene prefers the migration to the dissociation so that the subsequent [1,3]-H transfer occurs in ligand-associated complex. Based on a detailed electronic structure analysis on the metal-carbon bonds in complexes, we elucidated the reason why the differences in the ligand dissociation and migration exist between the Rh, Cu and Au carbenes. The mechanistic insights presented above give an unambiguous explanation in the chemo-selectivity of O-H and C-H bonds of phenols inserted by the Rh-TC, Cu-pybox, Cu-box and Au-TPPO carbenes, especially for the electronic structure effects of various transition metals on the reaction mechanisms. This will help us to take deep insight into the chemo-selectivity for different metal catalysts and find the most efficient catalyst for the functionalization of chemical bonds in organic molecules. Supporting Information The Cartesian coordinates of intermediates and transition states are listed in section 1 and the calculated pathways of O-H insertion from the methoxy group side are shown in Figure S1 of section 2. Corresponding Authors *E-mail: [email protected]

Tel: +86-(0)21-20596009.

ACKNOWLEDGEMENTS This work was supported by the National Natural Science Foundation of China (Grants No. 21433004 and 21473056) and Natural Science Foundation of Shanghai

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 12 of 28

(14ZR1411800). We acknowledge the support of the NYU-ECNU Center for Computational Chemistry at NYU Shanghai. We also thank the supercomputer center of ECNU for providing computer time.

Reference: (1)Kuninobu, Y.; Matsuki, T.; Takai, K., Rhenium-Catalyzed Regioselective Alkylation of Phenols. J. Am. Chem. Soc. 2009, 131, 9914-9915. (2)Ciana, C. L.; Phipps, R. J.; Brandt, J. R.; Meyer, F. M.; Gaunt, M. J., A Highly Para-Selective Copper(II)-Catalyzed Direct Arylation of Aniline and Phenol Derivatives. Angew. Chem. Int. Ed. 2011, 50, 458-462. (3)Dai, H. X.; Li, G.; Zhang, X. G.; Stepan, A. F.; Yu, J. Q., Pd(II)-Catalyzed ortho- or meta-C-H Olefination of Phenol Derivatives. J. Am. Chem. Soc. 2013, 135, 7567-7571. (4)Huang, Z. L.; Jin, L. Q.; Feng, Y.; Peng, P.; Yi, H.; Lei, A. W., Iron-Catalyzed Oxidative Radical Cross-Coupling/Cyclization between Phenols and Olefins. Angew. Chem. Int. Ed. 2013, 52, 7151-7155. (5)Sevov,

C.

S.;

Hartwig,

J.

F.,

Iridium-Catalyzed

Intermolecular

Hydroetherification of Unactivated Aliphatic Alkenes with Phenols. J. Am. Chem. Soc. 2013, 135, 9303-9306. (6)Wu, Z. Q.; Luo, F. H.; Chen, S.; Li, Z. K.; Xiang, H. F.; Zhou, X. G., Palladium-catalyzed Para-selective Arylation of Phenols with Aryl Iodides in Water. Chem. Commun. 2013, 49, 7653-7655. (7)Zhu, R. Y.; Wei, J. B.; Shi, Z. J., Benzofuran Synthesis via Copper-mediated Oxidative Annulation of Phenols and Unactivated Internal Alkynes. Chem. Sci. 2013, 4, 3706-3711.

ACS Paragon Plus Environment

Page 13 of 28

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

The Journal of Physical Chemistry

(8)Seoane, A.; Casanova, N.; Quinones, N.; Mascarenas, J. L.; Gulias, M., Straightforward Assembly of Benzoxepines by Means of a Rhodium(III)-Catalyzed C-H Functionalization of o-Vinylphenols. J. Am. Chem. Soc. 2014, 136, 834-837. (9)Davies, H. M. L.; Beckwith, R. E. J., Catalytic Enantioselective C-H Activation by means of Metal-carbenoid-induced C-H Insertion. Chem. Rev. 2003, 103, 2861-2903. (10)Davies, H. M. L.; Manning, J. R., Catalytic C-H Functionalization by Metal Carbenoid and Nitrenoid Insertion. Nature 2008, 451, 417-424. (11)Díaz-Requejo, M. M.; Pérez, P. J., Coinage Metal Catalyzed C-H Bond Functionalization of Hydrocarbons. Chem. Rev. 2008, 108, 3379-3394. (12)Doyle, M. P.; Duffy, R.; Ratnikov, M.; Zhou, L., Catalytic Carbene Insertion into C-H Bonds. Chem. Rev. 2010, 110, 704-724. (13)Zhu, S. F.; Zhou, Q. L., Transition-Metal-Catalyzed Enantioselective Heteroatom-Hydrogen Bond Insertion Reactions. Acc. Chem. Res. 2012, 45, 1365-1377. (14)Gillingham, D.; Fei, N., Catalytic X-H Insertion Reactions based on Carbenoids. Chem. Soc. Rev. 2013, 42, 4918-4931. (15)Xie, J.; Pan, C. D.; Abdukader, A.; Zhu, C. J., Gold-catalyzed C(sp3)-H Bond Functionalization. Chem. Soc. Rev. 2014, 43, 5245-5256. (16)Liu, L.; Zhang, J. L., Gold-catalyzed Transformations of Alpha-diazocarbonyl Compounds: Selectivity and Diversity. Chem. Soc. Rev. 2016, 45, 506-516. (17)Morilla, M. E.; Molina, M. J.; Díaz-Requejo, M. M.; Belderrain, T. R.; Nicasio, M. C.; Trofimenko, S.; Pérez, P. J., Copper-catalyzed Carbene Insertion into O-H Bonds: High Selective Conversion of Alcohols into Ethers. Organometallics 2003, 22, 2914-2918. (18)Fructos, M. R.; de Fremont, P.; Nolan, S. P.; Díaz-Requejo, M. M.; Pérez, P. J., Alkane carbon-hydrogen bond functionalization with (NHC)MCl precatalysts (M = Cu, Au ; NHC=N-heterocyclic carbene). Organometallics 2006, 25, 2237-2241. (19)Maier, T. C.; Fu, G. C., Catalytic Enantioselective O-H Insertion Reactions. J. Am. Chem. Soc. 2006, 128, 4594-4595.

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 14 of 28

(20)Chen, C.; Zhu, S. F.; Liu, B.; Wang, L. X.; Zhou, Q. L., Highly Enantioselective Insertion of Carbenoids into O-H Bonds of Phenols: An Efficient Approach to Chiral Alpha-aryloxycarboxylic Esters. J. Am. Chem. Soc. 2007, 129, 12616-12617. (21)Zhu, S. F.; Chen, C.; Cai, Y.; Zhou, Q. L., Catalytic Asymmetric Reaction with Water: Enantioselective Synthesis of Alpha-Hydroxyesters by a Copper-carbenoid O-H Insertion Reaction. Angew. Chem. Int. Ed. 2008, 47, 932-934. (22)Yates, P., The Copper-Catalyzed Decomposition of Diazoketones. J. Am. Chem. Soc. 1952, 74, 5376-5381. (23)Osako, T.; Panichakul, D.; Uozumi, Y., Enantioselective Carbenoid Insertion into Phenolic O-H Bonds with a Chiral Copper(I) Imidazoindolephosphine Complex. Org. Lett. 2012, 14, 194-197. (24)Demonceau,

A.;

Noels,

A.

F.;

Hubert,

A.

J.;

Teyssie,

P.,

Transition-Metal-Catalyzed Reactions of Diazoesters - Insertion into C-H Bonds of Paraffins by Carbenoids. J. Chem. Soc. Chem. Commun. 1981, 14, 688-689. (25)Demonceau,

A.;

Noels,

A.

F.;

Hubert,

A.

J.;

Teyssie,

P.,

Transition-Metal-Catalyzed Reactions of Diazoesters - Insertion into C-H Bonds of Paraffins Catalyzed by Bulky Rhodium(Ii) Carboxylates - Enhanced Attack on Primary C-H Bonds. Bull. Soc. Chim. Belg. 1984, 93, 945-948. (26)Padwa, A.; Austin, D. J.; Price, A. T.; Semones, M. A.; Doyle, M. P.; Protopopova, M. N.; Winchester, W. R.; Tran, A., Ligand Effects on Dirhodium(II) Carbene Reactivities - Highly Effective Switching between Competitive Carbenoid Transformations. J. Am. Chem. Soc. 1993, 115, 8669-8680. (27)Qu, Z. H.; Shi, W. F.; Wang, J. B., Linear Free Energy Correlation Analysis on the Electronic Effects of Rh(II) carbene O-H Insertion. J. Org. Chem. 2004, 69, 217-219. (28)Hu, F. D.; Xia, Y.; Ye, F.; Liu, Z. X.; Ma, C.; Zhang, Y.; Wang, J. B., Rhodium(III)-Catalyzed

ortho

Alkenylation

of

N-Phenoxyacetamides

with

N-Tosylhydrazones or Diazoesters through C-H Activation. Angew. Chem. Int. Ed. 2014, 53, 1364-1367.

ACS Paragon Plus Environment

Page 15 of 28

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

The Journal of Physical Chemistry

(29)Jia, S. K.; Xing, D.; Zhang, D.; Hu, W. H., Catalytic Asymmetric Functionalization

of

Aromatic

C-H

Bonds

by

Electrophilic

Trapping

of

Metal-Carbene-Induced Zwitterionic Intermediates. Angew. Chem. Int. Ed. 2014, 53, 13098-13101. (30)Fructos, M. R.; Belderrain, T. R.; de Fremont, P.; Scott, N. M.; Nolan, S. P.; Díaz-Requejo, M. M.; Pérez, P. J., A Gold Catalyst for Carbene-transfer Reactions from Ethyl Diazoacetate. Angew. Chem. Int. Ed. 2005, 44, 5284-5288. (31)Rivilla, I.; Gomez-Emeterio, B. P.; Fructos, M. R.; Díaz-Requejo, M. M.; Pérez, P. J., Exclusive Aromatic vs Aliphatic C-H Bond Functionalization by Carbene Insertion with Gold-Based Catalysts. Organometallics 2011, 30, 2855-2860. (32)Kilpin, K. J.; Paul, U. S. D.; Lee, A. L.; Crowley, J. D., Gold(I) "click'' 1,2,3-triazolylidenes: Synthesis, Self-assembly and Catalysis. Chem. Commun. 2011, 47, 328-330. (33)Pérez, P. J.; Díaz-Requejo, M. M.; Rivilla, I., Gold-catalyzed Naphthalene Functionalization. Beilstein J. Org. Chem. 2011, 7, 653-657. (34)Delgado-Rebollo, M.; Beltran, A.; Prieto, A.; Díaz-Requejo, M. M.; Echavarren, A. M.; Pérez, P. J., Catalytic Hydrocarbon Functionalization with Gold Complexes Containing N-Heterocyclic Carbene Ligands with Pendant Donor Groups. Eur. J. Inorg. Chem. 2012, 2012, 1380-1386. (35)Nakamura, E.; Yoshikai, N.; Yamanaka, M., Mechanism of C-H Bond Activation/C-C Bond Formation Reaction between Diazo Compound and Alkane Catalyzed by Dirhodium Tetracarboxylate. J. Am. Chem. Soc. 2002, 124, 7181-7192. (36)Braga, A. A. C.; Maseras, F.; Urbano, J.; Caballero, A.; Díaz-Requejo, M. M.; Pérez, P. J, Mechanism of Alkane C-H Bond Activation by Copper and Silver Homoscorpionate Complexes. Organometallics 2006, 25, 5292-5300. (37)Liang, Y.; Zhou, H. L.; Yu, Z. X., Why Is Copper(I) Complex More Competent Than Dirhodium(II) Complex in Catalytic Asymmetric O-H Insertion Reactions? A Computational Study of the Metal Carbenoid O-H Insertion into Water. J. Am. Chem. Soc. 2009, 131, 17783-17785.

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 16 of 28

(38)Bonge, H. T.; Hansen, T., Computational Study of C-H Insertion Reactions with Ethyl Bromodiazoacetate. Eur. J. Org. Chem. 2010, 2010, 4355-4359. (39)Xie, Z. Z.; Liao, W. J.; Cao, J.; Guo, L. P.; Verpoort, F.; Fang, W. H., Mechanistic Insight into the Rhodium-Catalyzed O-H Insertion Reaction: A DFT Study. Organometallics 2014, 33, 2448-2456. (40)Liu, Y.; Y, Z. Z.; Luo, Z. J.; Zhang, Z. H.; Liu, L.; Xia, F., Mechanistic Investigation of Aromatic C(sp2-H) and Alkyl C(sp3-H) Bond Insertion by Gold Carbenes. J. Phys. Chem. A 2016, 120, 1925-1932. (41)Xie, Q.; Song, X. S.; Qu, D. Y.; Guo, L. P.; Xie, Z. Z., DFT Study on the Rhodium(II)-Catalyzed C-H Functionalization of Indoles: Enol versus Oxocarbenium Ylide. Organometallics 2015, 34, 3112-3119. (42)Yu, Z. Z.; Ma, B.; Chen, M. J.; Wu, H. H.; Liu, L.; Zhang, J. L., Highly Site-Selective

Direct

C-H

Bond

Functionalization

of

Phenols

with

alpha-Aryl-alpha-diazoacetates and Diazooxindoles via Gold Catalysis. J. Am. Chem. Soc. 2014, 136, 6904-6907. (43)Xi, Y. M.; Su, Y. J.; Yu, Z. Y.; Dong, B. L.; McClain, E. J.; Lan, Y.; Shi, X. D., Chemoselective

Carbophilic

Addition

of

alpha-Diazoesters

through

Ligand-Controlled Gold Catalysis. Angew. Chem. Int. Ed. 2014, 53, 9817-9821. (44)Liu, Y.; Yu, Z. Z.; Zhang, J. Z.; Liu, L.; Xia, F.; Zhang, J. L., Origins of Unique Gold-catalysed chemo- and site-selective C-H Functionalization of Phenols with Diazo Compounds. Chem. Sci. 2016, 7, 1988-1995. (45)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.; et al. Gaussian 09, Gaussian, Inc.: Wallingford, CT, USA, 2009. (46)Zhao, Y.; Truhlar, D. G., The M06 Suite of Density Functionals for Main Group Thermochemistry, Thermochemical Kinetics, Noncovalent Interactions, Excited States, and Transition Elements: Two New Functionals and Systematic Testing of Four M06-Class Functionals and 12 Other Functionals. Theor. Chem. Acc. 2008, 120, 215-241.

ACS Paragon Plus Environment

Page 17 of 28

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

The Journal of Physical Chemistry

(47)Zhao, Y.; Truhlar, D. G., Density Functionals with Broad Applicability in Chemistry. Acc. Chem. Res. 2008, 41, 157-167. (48)Hay, P. J.; Wadt, W. R., Abinitio Effective Core Potentials for Molecular Calculations - Potentials for the Transition-Metal Atoms Sc to Hg. J. Chem. Phys. 1985, 82, 270-283. (49)Wadt, W. R.; Hay, P. J., Abinitio Effective Core Potentials for Molecular Calculations - Potentials for Main Group Elements Na to Bi. J. Chem. Phys. 1985, 82, 284-298. (50)Rassolov, V. A.; Ratner, M. A.; Pople, J. A., Semiempirical Models for Image Electrostatics. I. Bare External Charge. J. Chem. Phys. 2001, 114, 2062-2066. (51)Scalmani, G.; Frisch, M. J., Continuous Surface Charge Polarizable Continuum Models of Solvation. I. General formalism. J. Chem. Phys. 2010, 132,114110. (52)Glendening, E. D. R., A. E.; Carpenter, J. E.; Weinhold, F., NBO, Version3.1. (53)Wiberg, K. B., Application of Pople-Santry-Segal Cndo Method to Cyclopropylcarbinyl and Cyclobutyl Cation and to Bicyclobutane. Tetrahedron 1968, 24, 1083-1096. (54)Dennington, R.; Keith, T.; Millam, J.; Inc., S.; Shawnee Mission, K. 2009. GaussView, Version 5. (55)Zhang, X.; Huang, H. X.; Guo, X.; Guan, X. Y.; Yang, L. P.; Hu, W. H., Catalytic Enantioselective Trapping of an Alcoholic Oxonium Ylide with Aldehydes: Rh-II/Zr-IV-co-catalyzed Three-component Reactions of Aryl Diazoacetates, Benzyl Alcohol, and Aldehydes. Angew. Chem. Int. Ed. 2008, 47, 6647-6649.

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Figures and Captions: Figure 1. Optimized structures of (a) Phenol, (b) Rh-TC carbene, (c) Cu-pybox carbene, (d) Cu-box carbene and (e) Au-TPPO carbene, where the values of NBO charges and Wiberg bond orders are shown in black and red, respectively. Figure 2. Calculated free energy profiles ∆Gsol of the phenol O-H bond insertion catalyzed by Rh-TC carbene as well as the corresponding structures of intermediates and transition states along different pathways. The lowest energy pathway is shown in red and the other marked in black. All the values of free energies are in the units of kcal/mol.

Figure 3. Calculated free energy profiles ∆Gsol of the phenol C-H bond insertion catalyzed by Rh-TC carbene as well as the corresponding structures of intermediates and transition states along different pathways. The lowest energy pathway is shown in red and the other marked in black. All the values of free energies are in the units of kcal/mol.

Figure 4. Calculated competitive pathways of the phenol O-H and C-H insertions catalyzed by Cu-pybox carbene. The red and black curves denote the calculated lowest-energy pathways of O-H and C-H bond insertions, respectively.

Figure 5. Calculated competitive pathways of the phenol O-H and C-H insertions catalyzed by Cu-box carbene. The red and black curves denote the calculated lowest-energy pathways of O-H and C-H bond insertions, respectively.

Figure 6. Contour maps and schematic representations for (a) the 116th occupied molecular orbital of Rh-TC carbene, (b) the 115th occupied molecular orbital of Rh-TC carbene, (c) the HOMO of Cu-box carbene and (d) the 118th occupied molecular orbital of Au-TPPO carbene.

ACS Paragon Plus Environment

Page 18 of 28

Page 19 of 28

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

The Journal of Physical Chemistry

Figure 7. Optimized structures of intermediates and metal ligands, (a) Int-o2, (b) Rh-TC, (c) Int-o5, (d) Cu-pybox, (e) Int-o7, (f) Cu-box, (g) Au-enol and (h) Au-TPPO, where the values of NBO charges and Wiberg bond orders are shown in black and red. Figure 8. Contour maps of (a) the 142th occupied molecular orbital of Int-o2, (b) the HOMO of Int-o7 and (c) the HOMO of Au-enol.

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Figure 1.

ACS Paragon Plus Environment

Page 20 of 28

Page 21 of 28

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

The Journal of Physical Chemistry

Figure 2.

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Figure 3.

ACS Paragon Plus Environment

Page 22 of 28

Page 23 of 28

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

The Journal of Physical Chemistry

Figure 4.

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Figure 5.

ACS Paragon Plus Environment

Page 24 of 28

Page 25 of 28

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

The Journal of Physical Chemistry

Figure 6.

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Figure 7.

ACS Paragon Plus Environment

Page 26 of 28

Page 27 of 28

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

The Journal of Physical Chemistry

Figure 8.

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

TOC:

ACS Paragon Plus Environment

Page 28 of 28