First-Principles Molecular Dynamics Analysis of Ligand-Free Suzuki

Dec 19, 2016 - Phone: +81(0)6 68797288. ... water using first-principles molecular dynamics simulations, with solvent H2O molecules explicitly include...
0 downloads 0 Views 4MB Size
Subscriber access provided by BOSTON UNIV

Article

First-Principles Molecular Dynamics Analysis of Ligand-Free SuzukiMiyaura Cross Coupling in Water Solvent: Oxidative Addition Step Teruo Hirakawa, Yuta Uramoto, Daisuke Mimura, Atsuya Takeda, Susumu Yanagisawa, Takashi Ikeda, Kouji Inagaki, and Yoshitada Morikawa J. Phys. Chem. B, Just Accepted Manuscript • DOI: 10.1021/acs.jpcb.6b08644 • Publication Date (Web): 19 Dec 2016 Downloaded from http://pubs.acs.org on December 22, 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 B 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 31

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

First-Principles Molecular Dynamics Analysis of Ligand-Free Suzuki-Miyaura Cross Coupling in Water Solvent: Oxidative Addition Step Teruo Hirakawa,∗,† Yuta Uramoto,† Daisuke Mimura,† Atsuya Takeda,† Susumu Yanagisawa,‡ Takashi Ikeda,¶ Kouji Inagaki,†,§ and Yoshitada Morikawa∗,†,§,∥ Department of Precision Science and Technology, Graduate School of Engineering, Osaka University, 2-1, Yamada-oka, Suita, Osaka 565-0871, Japan , Department of Physics and Earth Sciences, Faculty of Science, University of the Ryukyus, 1 Senbaru, Nishihara, Okinawa 903-0213, Japan, Synchrotron Radiation Research Center, Quantum Beam Science Research Directorate, National Institutes for Quantum and Radiological Science and Technology (QST), 1-1-1 Kouto, Sayo, Hyogo 679-5148, Japan, Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Katsura, Kyoto 615-8520, Japan, and Research Center for Ultra-Precision Science and Technology, Graduate School of Engineering, Osaka University, 2-1, Yamada-oka, Suita, Osaka 565-0871, Japan E-mail: [email protected]; [email protected] Phone: +81(0)6 68797288. Fax: +81(0)6 68797290



To whom correspondence should be addressed Osaka University ‡ University of the Ryukyus ¶ QST § Kyoto University ∥ Research Center, Osaka University †

1 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

Abstract We investigated the oxidative addition step of PhX (X=Cl, Br) to a single Pd(0) atom or a PdX− complex in water by using first-principles molecular dynamics simulations with solvent H2 O molecules explicitly included in calculation models to clarify the origin for the extremely high reactivity of a ligand-free Pd catalyst in an aqueous solution for the Suzuki-Miyaura reaction. The free energy profiles are estimated by using blue moon ensemble (BME) sampling to include the entropy effect in chemical reactions in a water solvent. The free energy barrier of the oxidative addition step is quite low especially for PhBr while the barrier for PhCl is sizable, indicating that the reaction can proceed at room temperature with a high rate for PhBr, but the rate is rather slow for PhCl. We also investigated the effect of the additional halogen anion to the Pd catalyst as a "supporting ligand". The activation barrier of the oxidative addition step is not affected by the supporting halogen ligands, but the stability of the final state is significantly destabilized, which should be important for the following transmetalation step. The solvent effect is also investigated and discussed.

Introduction The Suzuki-Miyaura reaction (SMR) is one of the most efficient palladium-catalyzed crosscoupling reactions for the formation of carbon-carbon bonds, in particular for the formation sp2 -sp2 carbon-carbon of biaryls in organic synthesis. 1–8 The SMR has been developed mainly with ligands, especially phosphine ligands to make the palladium catalyst stable and optimize its reactivity in an organic solvent. The catalytic cycle for the cross-coupling reaction consists of three main steps, i.e., oxidative addition, transmetalation, and reductive elimination (see Figure 1). The SMR without any ligands, that is "ligand free" in a water solvent, was suggested, and it has been attracting enormous interest recently. 9–12 In 2003, Smith and co-workers showed that palladium-containing perovskites are also quite reactive catalysts for the SMR in an aqueous solvent. 9 In 2005, Andrews and co-workers demonstrated that its 2 ACS Paragon Plus Environment

Page 2 of 31

Page 3 of 31

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

turn-over number (TON) is over 399,000/Pd when PhBr is used as one of the reactants, 10 and Arisawa and co-workers reported higher TON up to 2,760,000 using sulfur-modified Ausupported Pd (SAPd), 13,14 which implies that the “ ligand-free ” SMR in water or ethanol solvents is much more efficient than that in organic solvents. Aryl bromide is still suitable as the reactant of the SMR in a water solvent, while aryl chloride cannot be used as aryl halide because of its low rate, which is similar to the ordinary SMR in an organic solvent. Although the ligand-free palladium is very promising and has been studied intensively, the active form of the catalyst is still in debate. Perovskite materials containing palladium were originally prepared as a heterogeneous catalyst, but it was shown that catalytically active Pd species are not on the perovskite surface but dissolved from the perovskite into the water solvent. 10 Reetz and co-workers and de Vries and co-workers proposed the leaching mechanism of the active Pd(0) catalyst from palladium nanoparticles by oxidative addition for the ligand-free Heck reaction, 15,16 and the leaching mechanism is also considered to be applicable for the SMR. 17–19 Leyva-Perez and co-workers observed the increase of the amount of small Pd clusters with a majority size of three or four atoms either by water or aryl bromide and concluded that they are active species as the homogeneous Pd catalyst for the SMR without any ligands in a water solvent. 20 Dhital and co-workers concluded that Pd20 nanoclusters are not so active for the oxidative addition of chlorobenzene using quantum chemical calculations. 21 de Vries and co-workers reported that a dipalladium intermediate has been observed in their experiment by using EXAFS. 22 The dipalladium intermediate can be obtained from the oxidative addition of aryl halide to the Pd catalyst or from the by-products, with some halogen anions (Cl− or I− ). In ligand-free systems, water or halogen anions in a solution may bond to the Pd catalyst and play some roles partly similar to phosphine ligands in conventional catalysts. 23 Since the seminal work by Amatore and Jutand, 24 a Pdhalogen anion complex has been investigated as an active species for cross-coupling reactions. The halogen anion can act as a "ligand" without any phosphine ligands in a water solvent. Moreover, the formation of the anion complex may inhibit the aggregation of Pd to form

3 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

Pd black, extending the TON of catalysts. The oxidative addition with ligand systems has been studied by Kozuch and co-workers, and they reported that a tri-coordinated anionic Pd catalyst can still be reactive. 25 Surawatanawong and Hall suggested that the halogens can play roles similar to ligands to enhance the oxidative addition step under phosphine-free conditions for the Heck reaction. 26 The SMR without any ligands in a water solvent has the higher TON; therefore, the energy barriers for the reaction processes are expected to be lower than those of the conventional SMR with phosphine-ligands. However, as far as we know, the reaction mechanisms of the SMR in a water solvent have not been studied theoretically. In many cases, the oxidative addition is the rate determining step of the catalytic cycle. 27,28 Thus, in the present study we focus on the oxidative addition step of the SMR in a water solvent explicitly including solvent molecules in our calculation models and perform full quantum mechanical molecular dynamics (MD) simulations. Most theoretical studies for the SMR reported so far, used a simple solvent model to simulate the cross coupling mechanisms in an organic solvent. 21,25,26,29–64 In these studies, the solvent effect is approximately treated using a polarizable continuum model (PCM). Sikk and co-workers showed the importance of the solvent effect for oxidative addition steps, especially to the Pd-anion complex. 51 Only a few works used explicit solvent models. 57,61,62,65 However, solvent molecules may play important roles in ligand-free Pd cross coupling reactions. Vidossich and co-workers reported solvent coordination in the Pd-phosphine system by performing first-principles MD simulations explicitly including solvent molecules. 62 Furthermore, Besora and co-workers pointed out that a conformation search is important in cross-coupling reactions with bulky ligands. 48

4 ACS Paragon Plus Environment

Page 4 of 31

Page 5 of 31

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

Methods All of the calculations are based on the density functional theory (DFT) within the generalized gradient approximation(GGA-PBE) as implemented in the STATE-Senri code. 66 We used ultrasoft pseudo potentials and a plane-wave basis set with cutoff energies of 25 Ry and 225 Ry for wave functions and charge densities, respectively. We performed first-principles MD simulations for the oxidative addition step of PhBr/PhCl over the ligand-free palladium catalyst with water solvent molecules explicitly included (see Figure 2). The supercell size used in our calculations is 12.38 Å × 12.38 Å × 12.38 Å and a single k-point was sampled for the Brillouin zone sampling. In the cell, 64 H2 O molecules are contained in the case of normal water (1.0 g/cm3 ). To simulate the oxidative addition step, one PhBr and one Pd atom or PdBr− anion were introduced into the cell and five H2 O molecules were instead removed from the cell to keep the water density constant. The catalyst model of Pd is still unknown, and here, we assumed Pd single atom or Pd-X− (X= Br or Cl) anion complexes. Constant temperature MD simulations were performed using the velocity Verlet integration scheme with a time-step of 0.75 fs per MD step. This time step is slightly longer than those typically used in classical MD simulations for water. Therefore, we actually used a deuteron mass rather a hydrogen mass for the H atoms in the water solvent. A Nosé-Hoover chain (NHC) with a chain length of eight was used, and the temperature target was set to 400K. Because the GGA-PBE tends to overestimate the hydrogen bond strength and accordingly the diffusion constant of water molecules is underestimated, the target temperature was set about 60 K higher than that in the actual experiments to mimic the real dynamics of the water solvent at experimental conditions. In order to accurately estimate the free energy barriers for the oxidative addition step in the water solvent, we used a blue moon ensemble (BME) method and a constrained MD to sample rare events, to calculate the mean force as a function of an arbitrary reaction coordinate ξ = ξ(r), and to construct the free energy profile. 67,68 In calculations of free energy profiles with the constrained MD by using the mean force, the error was estimated 5 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

by the method of Jacucci and Rahman, using the concept of the correlation length between those samples. 69,70 We imposed a constraint on the distance between the Pd atom and the bond center between the C1 atom of the benzene ring and the halogen atom X (X = Br or Cl), and it is hereafter denoted by RPd−BC (see Figure 3). The accuracy of free-energy calculations with BME depends on the reaction coordinates, and RPd−BC seems to be reasonably good for a wide range of values from 0 Å to 7 Å. We noticed, however, a small jump of atomic geometries close to the transition states. Therefore, we used an angle between Pd-Br and Pd-C1 bonds (θBrPdC ), and the results are shown and discussed in the Results and Discussion section. We performed over 5000 MD simulations (time step: 0.75 fs, total time: 3.75 ps) for each point of constraint value along the reaction coordinate.

Results and Discussion Oxidative Addition Step In the oxidative addition reaction for the SMR in an organic solvent, two competitive pathways have been suggested. Senn and Ziegler found an SN 2 pathway, 29,30 while several other research groups supported a concerted pathway. 34,42,44,49 Other competitive pathways are proposed by Hartwig, namely oxidative addition and phosphine dissociation. 71–74 Recent studies concluded that favored pathways depend on both ligands and solvents, for example, Besora and co-workers reported that in the case of palladium bis-phosphine in polar solvent, the SN 2 mechanism is favored, while in the case of palladium mono-phosphine with a polar solvent, both of the two mechanisms become competitive. 49 Mollar and co-workers pointed out that the SN 2 pathway is preferred for the activation of C sp3 -Br bond, while the concerted pathway is preferred for the activation of C sp2 -Br bond. 45 Therefore, in the present study, we investigated the oxidative addition step of X-Ph in the concerted pathway. We use a single Pd(0) atom and PdX− anion complex in a water solvent as the simplest 6 ACS Paragon Plus Environment

Page 6 of 31

Page 7 of 31

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

and most probable models for the ligand-free Pd catalyst in water and investigate the free energy profile for the concerted pathway of the oxidative addition to discuss the effects of the ligand and solvents. To observe the process of the oxidative addition step of the ligand-free SMR in the water solvent, we performed constrained MD with the BME sampling, starting with a Pd(0) atom and a PhX (X=Cl, Br) molecule separated and solvated in water, and ending in a Pd(II) XPh complex including a process of the isomerization from cis to trans. We used a single Pd atom as a basic model of a catalyst. Figure 4 shows the relative free energy profile as a function of RPd−BC (upper graph) along with the number of H2 O molecules connected to the Pd atom (coordination number, NC lower graph) and the atomic configurations of a PhX-Pd complex (top and bottom panels). The number of coordinated H2 O molecules to Pd is defined by the number of O-Pd bonds shorter than 2.4 Å. Before going into the details of the oxidative addition process of PhX to Pd in the water solvent, we check the accuracy of our calculations by comparing the present results with previously reported results. We calculated the free energy barrier for the oxidative addition of PhBr to mono-ligated PdPPh3 , and the free energy profile is shown with a green line in Figure 4. The free energy barrier for the oxidative addition of PhBr to monoligated PdPPh3 is calculated to be 2.4 kcal/mol, which is in good agreement with the value calculated for the same process by DFT with B3LYP hybrid functional. 41 Recently, the Maseras group reported a higher value of 9.6 kcal/mol as the activation energy for a similar oxidative addition by using a more sophisticated methodology, employing DFT with the M06 functional. 75 M06 includes dispersion correction, which increases the activation energies of the oxidative addition step systematically. However, the tendencies of the oxidative addition reaction are in good agreement between both functionals and we thus think that the main conclusion of our work is not affected by the inclusion of the van der Waals correction. In the initial state (RPd−BC = 7.0 Å), the Pd atom is solvated by two H2 O molecules, forming a bisligated Pd(H2 O)2 , as seen in Figures 2 and 4 (right upper panel). As the

7 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

PhX molecule approaches the Pd atom, one of the two H2 O molecules bonded to the Pd atom is replaced with the PhX molecule, forming an η 2 coordinated PhX to the Pd atom (RPd−BC ∼ 4.5 Å). During the replacement, the free energy is downhill and the system is stabilized by ∼ 20 kcal/mol, indicating that the process is significantly exothermic. After forming the η 2 coordination of PhX to the Pd atom, the Pd atom moves along the benzene ring and comes to the C1 atom which is directly bonded to X (RPd−BC ∼ 2.5 Å). Then Pd is inserted in between C1 and X atoms, and at the same time, the number of coordinated H2 O molecules increases from one to two, forming a square planar cis-Pd(II) XPh(H2 O)2 complex. The activation free energy barriers in this process are 3.4 kcal/mol for PhBr and 8.5 kcal/mol for PhCl. The activation free energies of the oxidative addition of PhBr to H2 O ligated Pd (3.4 kcal/mol) is slightly larger than that to the mono-ligated PdPPh3 (2.4 kcal/mol) mentioned above. According to the report by Jover et al., the reactivity of the oxidative addition mainly depends on the σ-donation property. 43 It is known that the water ligand bonded to the palladium atom has a weaker σ-donation property than the phosphine ligand. 76 Therefore, the trend in the activation energies of the oxidative addition to H2 O ligated Pd and phosphine Pd follows the general rule. These activation energies are in reasonable agreement with previous results for oxidative addition of PhX (X=Cl, Br) in mono-phosphine pathway, and it was pointed out that the difference in the activation barriers between PhBr and PhCl can be ascribed to the difference in the bond strength of the C-X bonds in the PhX. 39 The cis-Pd(II) XPh(H2 O)2 complex is more stable than the η 2 complex by −15.2 kcal/mol for PhBr and −11.3 kcal/mol for PhCl (RPd−BC ∼ 1.5 Å). Finally, the cis-trans isomerization process is simulated, and the free energy barriers are estimated as about 26.0 kcal/mol for PhBr and 17.9 kcal/mol for PhCl (RPd−BC ∼ 0.6 Å). As discussed in the Introduction, the oxidative addition step is postulated to be the rate limiting step of the whole catalytic cycle. However, in the present study we show that the oxidative addition step of PhBr to ligand-free Pd in the water solvent has a small barrier of 3.4 kcal/mol, indicating that the process proceeds quite smoothly at room temperature or above.

8 ACS Paragon Plus Environment

Page 8 of 31

Page 9 of 31

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

The low activation barrier should be the origin for the very high reactivity of the ligand-free Pd catalyst for the SMR. Furthermore, the large exothermic energy of the oxidative addition step suggests that if a ligand-free Pd atom exists in a solution, it immediately reacts with PhBr in the solution and forms a cis-Pd(II) XPh(H2 O)2 complex, inhibiting the aggregation of the Pd catalyst. This causes the long lifetime and the high TON of the catalyst. The activation barrier of the oxidative addition step for PhCl, however, is 8.5 kcal/mol and much higher than that for PhBr, slowing down the process. The higher barrier for the oxidative addition step leads to a longer lifetime of the η 2 complex, and this may cause the aggregation of the Pd catalyst as discussed by de Vries and co-workers. 16 In the isomerization from the cis-form to the trans-form of the Pd(II) XPh(H2 O)2 complex, there is a high energy barrier, namely, 26.0 kcal/mol for PhBr and 17.9 kcal/mol for PhCl; furthermore, the process is endothermic by 24.9 kcal/mol for PhBr and 15.1 kcal/mol for PhCl. Therefore, it is difficult to isomerize from the cis-form to the trans-form and even if it happens, the reverse reaction is easy to move back to the cis-form. Therefore, the product in the oxidative addition step of the ligand-free SMR in the water solvent is the cis-form of the Pd complex. This result is significantly different from that of Braga because Braga reported that the cis-trans isomerization happens easily due to the trans effect of PPh3 and H2 O. Thus, we think the transmetalation mechanism with bisligated Pd(H2 O)2 is rather different from that in the conventional SMR. This process is now under investigation and will be reported in our future work. Next, let us investigate the ligand effect of a halogen anion X− to Pd(0) on the oxidative addition reaction. Here, we investigated a single-atom Pd catalyst ligated by Br− or Cl− . In the ligand-free SMR, the Pd catalyst is solvated in a solution in which halogen anions exist and there is a possibility of forming PdX− anion complexes. We performed the BME sampling for the oxidative addition of PhX to PdX− (X=Cl, Br), and the results are shown in Figure 5 (right panel) along with the results without halogen anions. In the initial state, PhX is bonded to PdX− in η 2 coordination (see the lower right panel in Figure 5), while in the

9 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

final state, Pd forms the square planar (PdPhX2 H2 O)− complex. (see the lower left panel in Figure 5). The free energy barriers for the oxidative addition step are 5.1 kcal/mol for PhBr to PdBr− and 8.4 kcal/mol for PhCl to PdCl− , and exothermic energies are 10.1 kcal/mol for (PdPhBr2 H2 O)− and 8.1 kcal/mol for (PdPhCl2 H2 O)− . The coordination number of H2 O changed from zero to one during the oxidative addition to form the square planar (PdPhX2 H2 O)− complex. Compared with the oxidative addition to halogen-free Pd(0) , the oxidative addition to PdBr− /PdCl− has almost the same free energy barriers. Conversely, the exothermic energies to form the (PdPhX2 H2 O)− complexes are significantly reduced from those to form the PdPhX(H2 O)2 complex. The destabilization of the products of the oxidative addition process may accelerate the following reaction step, i.e, the transmetalation step. After the oxidative addition to PdBr− /PdCl− , two halogen anions are connected to Pd in trans-position to each other and Ph is connected to Pd in cis-position to the halogen anions. We also checked the stability of the Ph connected to Pd in trans-position to one of the two halogen anions, but we could not identify the stable final state, indicating that the trans-form of the Pd(II) X2 PhH2 O complex is much less stable than the cis-form of the complex, which is similar to the Pd(II) XPh(H2 O)2 complex as discussed above. Therefore, the following transmetalation step may be different from the conventional SMR with phosphine ligated Pd catalysts. To clarify the solvent effect, we investigated the oxidative addition steps of PhX to a Pd(0) H2 O complex and to a PdX− anion complex in a vacuum. In the case of a Pd(0) H2 O complex, we included one H2 O molecule because during the oxidation step one H2 O molecule is always bonded to the central Pd metal, as seen in Figure 4. In these calculations, we did not perform finite-temperature MD simulations, but we carried out simpler static reaction path calculations using a climbing image nudged elastic band method proposed by Henkelman and co-workers 77 because without a solvent the entropy effect is expected to be small. The calculated energy profiles are shown in the upper left panel of Figure 5. As seen in Figure 5, the effect of X− on the activation barriers of the oxidative addition is small, namely, the

10 ACS Paragon Plus Environment

Page 10 of 31

Page 11 of 31

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

barriers are slightly increased by about 3 kcal/mol in the present study, and the barriers are slightly decreased by ∼ 5 kcal/mol in previous studies on a CH3 X oxidative addition to Pd/PdCl− . 38,78 Conversely, the effect of the X− anion on the reaction energies is large, namely, the product of the oxidative addition step is stabilized by 8∼10 kcal/mol , which agrees with previous works. 38,78 We notice that the comparison with those previous studies on the CH3 X oxidative addition is valid for checking of the difference of solvent effects between the explicit solvent model and the continuum models. As for the solvent effect, the effect is small for the case of the oxidative addition step of PhX to a single Pd(0) H2 O complex. Conversely, for the oxidative addition step of PhX to the PdX− anion complex, the activation barriers do not depend on the inclusion of the solvent, but the exothermic energies do depend significantly on the inclusion of the solvent. This effect was observed previously using a continuous solvent model 38 but the effect seems to be larger in the present study using explicit solvent molecules. As seen in the right hand panel of Figure 5, by adding the halogen anion (X− ) to the Pd catalyst along with the solvent effect, the exothermic energy is reduced by about 5 kcal/mol for Br and by about 3 kcal/mol for Cl. The reduction of the exothermic energy will be important to enhance the following transmetalation step. Finally, we compare the solvent effects resulting from our explicit solvent model with that of the continuum model reported by Jong and Bickelhaupt. 38 In the case of the PhCl oxidative addition to Pd, the activation energy we calculated does not depend on the existence of the water solvent, and the final state energy is slightly more stabilized than that of the gas phase by about 3 kcal/mol, which is in good agreement with the case of the CH3 Cl oxidative addition to Pd using the continuum model. Whereas, in the case of the PhCl oxidative addition to PdCl− , the activation energy is 4.8 kcal/mol lower in solution than that in the gas phase; however, the activation barrier for the case of the CH3 Cl oxidative addition to PdCl− by using the continuum model is 4.9 kcal/mol higher in solution than that in the gas phase, having an opposite effect on the activation energy. The reaction energy is destabilized by 9.9 kcal/mol in solution more than that in the gas phase, which is in good agreement

11 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

with the destabilization of 10.3 kcal/mol in solution calculated by the continuum model. Therefore, especially for the anionic catalyst, the explicit solvent model is important for the calculation of the activation energy.

Tendency of the Pd Aggregation To validate the single Pd catalyst model, we investigated the possibility of Pd aggregation, namely, the formation of a Pd2 dimer from two single Pd catalysts or from two single PdBr− catalysts. We performed first-principles MD simulations of two single Pd (or PdBr− ) catalysts located close to each other and monitored whether two Pd atoms form a Pd2 dimer or not. If a Pd2 dimer is formed, this suggests that the Pd single catalyst tends to aggregate to form large Pd clusters. For each single Pd catalyst, namely, Pd or PdBr− , we considered two possible cases, i.e., a dimer formation before and after the oxidative addition step. Before the oxidative addition step, one of the two reactants, PhBr, forms an η 2 complex with the Pd catalyst (see Figure 6). We assumed that two η 2 complexes would come close to each other in the water solvent and one of the two Pd (or PdBr− ) moved from one η 2 complex to the other one, forming two Pd (or PdBr− ) catalysts bonded to a single PhBr with a Pd-Pd distance of 2.68 Å. After the oxidative addition step, the Pd catalyst forms a square-shaped Pd complex, namely, PhPdBr(H2 O)2 in the case of the single Pd catalyst (Figure 4) and PhPdBr2 (H2 O)− in the case of the single PdBr− catalyst (Figure 5). We assumed that the two square-shaped Pd complexes would come close to each another in the water solvent. Starting from these four initial configurations, we performed first-principles MD simulations and investigated the tendency of the Pd catalyst aggregation. For these simulations, we used slightly different unit cell models from those used in the simulations of the oxidative addition step. Prior to the oxidation step for the Pd catalyst without Br− , two Pd atoms bonded to one PhBr are introduced together with 54 H2 O molecules in a cubic box of 12.38 Å× 12.38 Å× 12.38 Å, while for the PdBr− catalyst, two PdBr− anions are introduced together with 86 H2 O molecules in a cubic box of 14.21 Å× 14.21 Å× 14.21 Å. Following the oxidation 12 ACS Paragon Plus Environment

Page 12 of 31

Page 13 of 31

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

step for the Pd catalyst without Br− (PdBr− catalyst), two PhPdBr (PhPdBr− 2 ) complexes are introduced together with 86 H2 O molecules in a cubic box of 14.21 Å× 14.21 Å× 14.21 Å. Figure 6 shows the time evolution of the Pd-Pd distances for the four cases. On the one hand, if two single Pd catalysts come close to each other before the oxidative addition, for both catalyst cases, namely, the single Pd catalyst (red dashed line) and the single PdBr− catalyst (black solid line), the Pd-Pd bond distances fluctuate between 2.5 and 3.0 Å, which is close to the metallic Pd-Pd distance of 2.75 Å, indicating the initial stage of the Pd aggregation. On the other hand, if two single Pd catalysts come close to each other after the oxidative addition, the Pd-Pd bond exceeds 3.0 Å for the case of the single Pd catalyst (blue dashed line), and it exceeds 3.5 Å for the case of the single PdBr− catalyst (green solid line); both bond lengths are much larger than the Pd-Pd metallic bond, suggesting that the Pd aggregation is unlikely. The results indicate that if the single Pd (or PdBr− ) catalysts come close to each other before the oxidative addition step, they easily form Pd-Pd metallic bonds and tend to aggregate, while if they meet each other after the oxidative addition step, they never form Pd-Pd metallic bonds and do not aggregate. Therefore, the lifetime of the η 2 complex is very important; namely, if the barrier of the oxidative addition step is large, and the lifetime of the η 2 complex is long, then the single Pd catalysts have a large probability to aggregate. Conversely, if the barrier of the oxidative addition step is small, and the lifetime of the η 2 complex is short, then the single Pd catalysts have a small probability to aggregate. The discussion is in accordance with the de Vries hypothesis. 16 Even after the oxidative addition, the Pd catalysts may form a Pd(Br2 )Pd complex as seen in Figure 6. In the case of the single PdBr− catalyst, however, the electrostatic repulsive interaction may significantly reduce the possibility of two PdBr− catalysts coming close to each other, making the aggregation more unlikely. Thus, we conclude that the Pd aggregation is unlikely especially for the case of the single PdBr− catalyst, provided that the free energy barrier for the oxidative addition step is reasonably low.

13 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

Validity of the Constraint for the Free Energy Calculation The quality of the free energy difference calculated using the BME method depends on the choice of reaction coordinates. The reaction coordinate RPd−BC we chose in our simulations is reasonably good for a wide range of values from 0 Å to ∼7 Å. We noticed, however, a small gap in geometries near the transition steps of the oxidative addition steps. The right panel of Figure 7 shows the values for RPd−BC and the lengths of the C1 -Cl along the paths obtained using constraints on RPd−C or on θClPdC . As seen in the figure, the length of the C1 -Cl bond for BME with the constraint on RPd−BC changes suddenly at around RPd−BC ∼ 2.0 Å, indicating that the constraint on the RBd−BC is not adequate to describe the present reaction path. Conversely, both RPd−BC and RC1 −Cl for BME with the constraint on θClPdC1 are smooth, and there is no gap. Therefore, θClPdC seems to be a better reaction coordinate to describe the present path. Still, the right hand panel of Figure 7 shows that the relative free energies calculated by the two BME simulations agree rather well with each other, indicating that the relative free energies can be calculated by both of the constraints. Although θClPdC1 seems to be better for describing the present reaction path, especially near the transition state, the error in the free energies calculated by the constraint on θClPdC1 becomes large; therefore, θClPdC1 cannot describe the whole range of the reaction path we have investigated. Therefore, we used RPd−BC as the constraint in the present simulations.

Conclusion We performed first-principles MD simulations of the oxidative addition step in a SMR with a ligand-free Pd catalyst in a water solvent to clarify the origin for the extremely high catalytic activity. We used a single Pd atom solvated in water and solvent H2 O molecules are explicitly included in our calculation model. We precisely calculated the relative free energy profiles including the solvent entropy by using the BME method. The results show that the energy barrier on the oxidative addition step is quite low especially for PhBr, which

14 ACS Paragon Plus Environment

Page 14 of 31

Page 15 of 31

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

is significantly different from the conventional SMR with ligands in organic solvents. This may be the reason why the SMR in water without any ligands can achieve a high TON and low-leaching of the Pd catalyst. Conversely, the barrier of the oxidative addition for PhCl is sizable, indicating that the reaction can proceed at room temperature or above, but the rate is rather slow compared with PhBr. Therefore, the lifetime of the η 2 complex of PhCl-Pd is elongated, causing the sintering of Pd to degrade the catalyst, which agrees with the de Vries idea for a ligand-free Heck reaction. 16 We also analyzed the by-product anion effect to the oxidative addition reaction of the ligand-free SMR in water. We observed that the activation barriers of the oxidative addition of PhX to the PdX− anion complex are similar to those of the oxidative addition of PhX to a single Pd(0) atom in a water solvent. The final state of the oxidative addition step, however, becomes significantly destabilized by the additional coordination of a halogen anion to the central Pd atom, and this should have an important effect for the following transmetalation step. The coordination of anions to the Pd atom should also prevent the aggregation of the Pd atom to form the Pd black. Therefore, the halogen anion can act as a "ligand" in a ligand-free system.

Acknowledgements We would like to thank Prof. Nobuaki Kambe, Prof. Takanori Iwasaki, Drs. Daiju Matsumura and Yasuo Nishihata for valuable discussions. The present study was partly supported by Grants- in Aid for Scientific Research on Innovative Areas 3D Active-Site Science (No. 26105010 and No. 26105011) and Scientific Research (C) (No. 26410014) from the Japan Society for the Promotion of Science (JSPS), the Elements Strategy Initiative for Catalysts and Batteries (ESICB) supported by the Ministry of Education, Culture, Sports, Science, and Technology, Japan (MEXT), and the JSPS Core-to-Core Program (Type A) "Advanced Research Networks: Computational Materials Design on Green Energy." The numerical calculations have been done with the facilities of the Supercomputer Center, In-

15 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

stitute for Solid State Physics, the University of Tokyo, Nagoya University, and Tohoku University.

References (1) Miyaura, N.; Suzuki, A. Stereoselective Synthesis of Arylated (E)-Alkenes by the Reaction of Alk-1-Enylboranes with Aryl Halides in the Presence of Palladium Catalyst. J. Chem. Soc., Chem. Commun. 1979, 19, 866–867. (2) Miyaura, N.; Yamada, K.; Suzuki, A. A New Stereospecific Cross-Coupling by the Palladium-Catalyzed Reaction of 1-Alkenylboranes with 1-Alkenyl or 1-Alkynyl Halides. Tetrahedron Lett. 1979, 20, 3437–3440. (3) Miyaura, N.; Yanagi, T.; Suzuki, A. The Palladium-Catalyzed Cross-Coupling Reaction of Phenylboronic Acid with Haloarenes in the Presence of Bases. Synth. Commun. 1981, 11, 513–519. (4) Miyaura, N.;

Suzuki, A. Palladium-Catalyzed Cross-Coupling Reactions of

Organoboron Compounds. Chem. Rev. 1995, 95, 2457–2483. (5) Suzuki, A. Recent Advances in the Cross-Coupling Reactions of Organoboron Derivatives with Organic Electrophiles, 1995–1998. J. Organomet. Chem. 1999, 576, 147–168. (6) Suzuki, A. Cross-Coupling Reactions via Organoboranes. J. Organomet. Chem. 2002, 653, 83–90. (7) Nicolaou, K.; Bulger, P. G.; Sarlah, D. Palladium-Catalyzed Cross-Coupling Reactions in Total Synthesis. Angew. Chem., Int. Ed. 2005, 44, 4442–4489. (8) Suzuki, A.; Yamamoto, Y. Cross-Coupling Reactions of Organoboranes: An Easy Method for C-C Bonding. Chem. Lett. 2011, 40, 894–901.

16 ACS Paragon Plus Environment

Page 16 of 31

Page 17 of 31

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

(9) Smith, M. D.; Stepan, A. F.; Ramarao, C.; Brennan, P. E.; Ley, S. V. PalladiumContaining Perovskites: Recoverable and Reuseable Catalysts for Suzuki Couplings. Chem. Commun. 2003, 2652–2653. (10) Andrews, S. P.; Stepan, A. F.; Tanaka, H.; Ley, S. V.; Smith, M. D. Heterogeneous or Homogeneous? A Case Study Involving Palladium-Containing Perovskites in the Suzuki Reaction. Adv. Synth. Catal. 2005, 347, 647–654. (11) Leadbeater, N. E. Fast, Easy, Clean Chemistry by Using Water as a Solvent and Microwave Heating: The Suzuki Coupling as an Illustration. Chem. Commun. 2005, 2881–2902. (12) Hill, N. J.; Bowman, M. D.; Esselman, B. J.; Byron, S. D.; Kreitinger, J.; Leadbeater, N. E. Ligand-Free Suzuki–Miyaura Coupling Reactions Using an Inexpensive Aqueous Palladium Source: A Synthetic and Computational Exercise for the Undergraduate Organic Chemistry Laboratory. J. Chem. Educ. 2014, 91, 1054–1057. (13) Hoshiya, N.; Shimoda, M.; Yoshikawa, H.; Yamashita, Y.; Shuto, S.; Arisawa, M. Sulfur Modification of Au via Treatment with Piranha Solution Provides Low-Pd Releasing and Recyclable Pd Material, SAPd. J. Am. Chem. Soc. 2010, 132, 7270–7272. (14) Arisawa, M.; Sato, T.; Hoshiya, N.; Al-Amin, M.; Kogami, Y.; Shuto, S. Ligand-Free Suzuki–Miyaura Coupling with Sulfur-Modified Gold-Supported Palladium in the Synthesis of a Conformationally-Restricted Cyclopropane Compound Library with ThreeDimensional Diversity. ACS Comb. Sci. 2014, 16, 215–220. (15) Reetz, M. T.; Westermann, E. Phosphane-Free Palladium-Catalyzed Coupling Reactions: the Decisive Role of Pd Nanoparticles. Angew. Chem., Int. Ed. 2000, 39, 165– 168. (16) de Vries, A. H. M.; Mulders, J. M. C. A.; Mommers, J. H. M.; Henderickx, H. J. W.;

17 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

de Vries, J. G. Homeopathic Ligand-Free Palladium as a Catalyst in the Heck Reaction. A Comparison with a Palladacycle. Org. Lett. 2003, 5, 3285–3288. (17) Pagliaro, M.; Pandarus, V.; Ciriminna, R.; Beland, F.; Demma Carà, P. Heterogeneous versus Homogeneous Palladium Catalysts for Cross-Coupling Reactions. ChemCatChem 2012, 4, 432–445. (18) Misch, L. M.; Birkel, A.; Figg, C. A.; Fors, B. P.; Hawker, C. J.; Stucky, G. D.; Seshadri, R. Rapid Microwave-Assisted Sol–Gel Preparation of Pd-Substituted LnFeO3 (Ln= Y, La): Phase Formation and Catalytic Activity. Dalton Trans. 2014, 43, 2079– 2087. (19) Kashin, A. N.; Ganina, O. G.; Cheprakov, A. V.; Beletskaya, I. P. The Direct NonPerturbing Leaching Test in the Phosphine-Free Suzuki–Miyaura Reaction Catalyzed by Palladium Nanoparticles. ChemCatChem 2015, 7, 2113–2121. (20) Leyva-Pérez, A.; Oliver-Meseguer, J.; Rubio-Marqués, P.; Corma, A. Water-Stabilized Three-and Four-Atom Palladium Clusters as Highly Active Catalytic Species in LigandFree C-C Cross-Coupling Reactions. Angew. Chem., Int. Ed. 2013, 52, 11554–11559. (21) Dhital, R. N.; Kamonsatikul, C.; Somsook, E.; Bobuatong, K.; Ehara, M.; Karanjit, S.; Sakurai, H. Low-Temperature Carbon–Chlorine Bond Activation by Bimetallic Gold/Palladium Alloy Nanoclusters: An Application to Ullmann Coupling. J. Am. Chem. Soc. 2012, 134, 20250–20253. (22) de Vries, J. G. A Unifying Mechanism for All High-Temperature Heck Reactions. The Role of Palladium Colloids and Anionic Species. Dalton Trans. 2006, 421–429. (23) Setiawan, D.; Kalescky, R.; Kraka, E.; Cremer, D. Direct Measure of Metal–Ligand Bonding Replacing the Tolman Electronic Parameter. Inorg. Chem. 2016, 55, 2332– 2344.

18 ACS Paragon Plus Environment

Page 18 of 31

Page 19 of 31

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

(24) Amatore, C.; Azzabi, M.; Jutand, A. Role and Effects of Halide Ions on the Rates and Mechanisms of Oxidative Addition of Iodobenzene to Low-Ligated Zerovalent Palladium Complexes Pd0 (PPh3 )2 . J. Am. Chem. Soc. 1991, 113, 8375–8384. (25) Kozuch, S.; Amatore, C.; Jutand, A.; Shaik, S. What makes for a Good Catalytic Cycle? A Theoretical Study of the Role of an Anionic Palladium(0) Complex in the Cross-Coupling of an Aryl Halide with an Anionic Nucleophile. Organometallics 2005, 24, 2319–2330. (26) Surawatanawong, P.; Hall, M. B. Theoretical Study of Alternative Pathways for the Heck Reaction through Dipalladium and “ Ligand-Free ” Palladium Intermediates. Organometallics 2008, 27, 6222–6232. (27) Smith, G. B.; Dezeny, G. C.; Hughes, D. L.; King, A. O.; Verhoeven, T. R. Mechanistic Studies of the Suzuki Cross-Coupling Reaction. J. Org. Chem. 1994, 59, 8151–8156. (28) Christmann, U.; Vilar, R. Monoligated Palladium Species as Catalysts in CrossCoupling Reactions. Angew. Chem., Int. Ed. 2005, 44, 366–374. (29) Bickelhaupt, F. M.; Ziegler, T. Oxidative Insertion as Frontside SN 2 Substitution: A Theoretical Study of the Model Reaction System Pd + CH3 Cl. Organometallics 1995, 14, 2288–2296. (30) Senn, H. M.; Ziegler, T. Oxidative Addition of Aryl Halides to Palladium(0) Complexes: A Density-Functional Study Including Solvation. Organometallics 2004, 23, 2980–2988. (31) Goossen, L. J.; Koley, D.; Hermann, H. L.; Thiel, W. Mechanistic Pathways for Oxidative Addition of Aryl Halides to Palladium(0) Complexes: A DFT Study. Organometallics 2005, 24, 2398–2410. (32) Braga, A. A.; Morgon, N. H.; Ujaque, G.; Maseras, F. Computational Characterization

19 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

of the Role of the Base in the Suzuki-Miyaura Cross-Coupling Reaction. J. Am. Chem. Soc. 2005, 127, 9298–9307. (33) Diefenbach, A.; de Jong, G. T.; Bickelhaupt, F. M. Activation of H–H, C–H, C–C and C–Cl Bonds by Pd and PdCl− . Understanding Anion Assistance in C–X Bond Activation. J. Chem. Theory Comput. 2005, 1, 286–298. (34) Braga, A. A.; Morgon, N. H.; Ujaque, G.; Lledós, A.; Maseras, F. Computational Study of the Transmetalation Process in the Suzuki–Miyaura Cross-Coupling of Aryls. J. Organomet. Chem. 2006, 691, 4459–4466. (35) Goossen, L. J.; Koley, D.; Hermann, H. L.; Thiel, W. Palladium Monophosphine Intermediates in Catalytic Cross-Coupling Reactions: A DFT Study. Organometallics 2006, 25, 54–67. (36) Ahlquist, M.; Fristrup, P.; Tanner, D.; Norrby, P.-O. Theoretical Evidence for LowLigated Palladium(0):[Pd-L] as the Active Species in Oxidative Addition Reactions. Organometallics 2006, 25, 2066–2073. (37) Ahlquist, M.; Norrby, P.-O. Oxidative Addition of Aryl Chlorides to Monoligated Palladium(0): A DFT-SCRF Study. Organometallics 2007, 26, 550–553. (38) de Jong, G. T.; Bickelhaupt, F. M. Catalytic Carbon-Halogen Bond Activation: Trends in Reactivity, Selectivity, and Solvation. J. Chem. Theory Comput. 2007, 3, 514–529. (39) Lam, K. C.; Marder, T. B.; Lin, Z. DFT Studies on the Effect of the Nature of the Aryl Halide Y-C6 H4 -X on the Mechanism of its Oxidative Addition to Pd0 L versus Pd0 L2 . Organometallics 2007, 26, 758–760. (40) Huang, Y.-L.; Weng, C.-M.; Hong, F.-E. Density Functional Studies on PalladiumCatalyzed Suzuki–Miyaura Cross-Coupling Reactions Assisted by N- or P-Chelating Ligands. Chem. Eur. J. 2008, 14, 4426–4434. 20 ACS Paragon Plus Environment

Page 20 of 31

Page 21 of 31

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

(41) Li, Z.; Fu, Y.; Guo, Q.-X.; Liu, L. Theoretical Study on Monoligated Pd-Catalyzed Cross-Coupling Reactions of Aryl Chlorides and Bromides. Organometallics 2008, 27, 4043–4049. (42) Gourlaouen, C.; Ujaque, G.; Lledós, A.; Medio-Simon, M.; Asensio, G.; Maseras, F. Why Is the Suzuki–Miyaura Cross-Coupling of sp3 Carbons in α-Bromo Sulfoxide Systems Fast and Stereoselective? A DFT Study on the Mechanism. J. Org. Chem. 2009, 74, 4049–4054. (43) Jover, J.; Fey, N.; Purdie, M.; Lloyd-Jones, G. C.; Harvey, J. N. A Computational Study of Phosphine Ligand Effects in Suzuki–Miyaura Coupling. J. Mol. Catal. A: Chem. 2010, 324, 39–47. (44) McMullin, C. L.; Jover, J.; Harvey, J. N.; Fey, N. Accurate Modelling of Pd(0)+ PhX Oxidative Addition Kinetics. Dalton Trans. 2010, 39, 10833–10836. (45) Mollar, C.; Besora, M.; Maseras, F.; Asensio, G.; Medio-Simón, M. Competitive and Selective Csp3 -Br versus Csp2 -Br Bond Activation in Palladium-Catalysed Suzuki CrossCoupling: An Experimental and Theoretical Study of the Role of Phosphine Ligands. Chem. Eur. J. 2010, 16, 13390–13397. (46) Proutiere, F.; Schoenebeck, F. Solvent Effect on Palladium-Catalyzed Cross-Coupling Reactions and Implications on the Active Catalytic Species. Angew. Chem., Int. Ed. 2011, 50, 8192–8195. (47) Kozuch, S.; Martin, J. M. L. What Makes for a Bad Catalytic Cycle? A Theoretical Study on the Suzuki–Miyaura Reaction within the Energetic Span Model. ACS Catal. 2011, 1, 246–253. (48) Besora, M.; Braga, A. A.; Ujaque, G.; Maseras, F.; Lledós, A. The Importance of Conformational Search: A Test Case on the Catalytic Cycle of the Suzuki–Miyaura Cross-Coupling. Theor. Chem. Acc. 2011, 128, 639–646. 21 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

(49) Besora, M.; Gourlaouen, C.; Yates, B.; Maseras, F. Phosphine and Solvent Effects on Oxidative Addition of CH3 Br to Pd(PR3 ) and Pd(PR3 )2 Complexes. Dalton Trans. 2011, 40, 11089–11094. (50) Mesganaw, T.; Silberstein, A. L.; Ramgren, S. D.; Nathel, N. F. F.; Hong, X.; Liu, P.; Garg, N. K. Nickel-Catalyzed Amination of Aryl Carbamates and Sequential SiteSelective Cross-Couplings. Chem. Sci. 2011, 2, 1766–1771. (51) Sikk, L.; Tammiku-Taul, J.; Burk, P.; Kotschy, A. Computational Study of the Sonogashira Cross-Coupling Reaction in the Gas Phase and in Dichloromethane Solution. J. Mol. Model. 2012, 18, 3025–3033. (52) Quasdorf, K. W.; Antoft-Finch, A.; Liu, P.; Silberstein, A. L.; Komaromi, A.; Blackburn, T.; Ramgren, S. D.; Houk, K.; Snieckus, V.; Garg, N. K. Suzuki- Miyaura CrossCoupling of Aryl Carbamates and Sulfamates: Experimental and Computational Studies. J. Am. Chem. Soc. 2011, 133, 6352–6363. (53) Li, Z.; Jiang, Y.-Y.; Yeagley, A. A.; Bour, J. P.; Liu, L.; Chruma, J. J.; Fu, Y. Mechanism of the Pd-Catalyzed Decarboxylative Allylation of α-Imino Esters: Decarboxylation via Free Carboxylate Ion. Chem. Eur. J. 2012, 18, 14527–14538. (54) Dong, L.; Wen, J.; Qin, S.; Yang, N.; Yang, H.; Su, Z.; Yu, X.; Hu, C. Iron-Catalyzed Direct Suzuki–Miyaura Reaction: Theoretical and Experimental Studies on the Mechanism and the Regioselectivity. ACS Catal. 2012, 2, 1829–1837. (55) García-Melchor, M.; Braga, A. A.; Lledós, A.; Ujaque, G.; Maseras, F. Computational Perspective on Pd-catalyzed C–C Cross-Coupling Reaction Mechanisms. Acc. Chem. Res. 2013, 46, 2626–2634. (56) Zhang, T.-X.; Li, Z. A DFT Study on Pd-Catalyzed Suzuki Cross-Coupling Polycondensation of Aryl Bromide Monomers. Comput. Theor. Chem. 2013, 1016, 28–35.

22 ACS Paragon Plus Environment

Page 22 of 31

Page 23 of 31

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

(57) Lyngvi, E.; Schoenebeck, F. Oxidative Addition Transition States of Pd(0) Complexes in Polar Solvent—A DFT Study Involving Implicit and Explicit Solvation. Tetrahedron 2013, 69, 5715–5718. (58) Lennox, A. J. J.; Lloyd-Jones, G. C. Transmetalation in the Suzuki–Miyaura Coupling: The Fork in the Trail. Angew. Chem., Int. Ed. 2013, 52, 7362–7370. (59) McMullin, C. L.; Fey, N.; Harvey, J. N. Computed Ligand Effects on the Oxidative Addition of Phenyl Halides to Phosphine Supported Palladium(0) Catalysts. Dalton Trans. 2014, 43, 13545–13556. (60) Boekfa, B.; Pahl, E.; Gaston, N.; Sakurai, H.; Limtrakul, J.; Ehara, M. C–Cl Bond Activation on Au/Pd Bimetallic Nanocatalysts Studied by Density Functional Theory and Genetic Algorithm Calculations. J. Phys. Chem. C 2014, 118, 22188–22196. (61) Ortuño, M. A.; Lledós, A.; Maseras, F.; Ujaque, G. The Transmetalation Process in Suzuki–Miyaura Reactions: Calculations Indicate Lower Barrier via Boronate Intermediate. ChemCatChem 2014, 6, 3132–3138. (62) Vidossich, P.; Ujaque, G.; Lledós, A. Palladium Monophosphine Pd(PPh3 ): Is It Really Accessible in Solution? Chem. Commun. 2014, 50, 661–663. (63) Sperger, T.; Sanhueza, I. A.; Kalvet, I.; Schoenebeck, F. Computational Studies of Synthetically Relevant Homogeneous Organometallic Catalysis Involving Ni, Pd, Ir, and Rh: An Overview of Commonly Employed DFT Methods and Mechanistic Insights. Chem. Rev. 2015, 115, 9532–9586. (64) Aufiero, M.; Scattolin, T.; Proutiere, F.; Schoenebeck, F. Air-Stable Dinuclear IodineBridged Pd(I) Complex-Catalyst, Precursor, or Parasite? The Additive Decides. Systematic Nucleophile-Activity Study and Application as Precatalyst in Cross-Coupling. Organometallics 2015, 34, 5191–5195.

23 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

(65) Vidossich, P.; Lledós, A.; Ujaque, G. First-Principles Molecular Dynamics Studies of Organometallic Complexes and Homogeneous Catalytic Processes. Acc. Chem. Res. 2016, 49, 1271–1278. (66) Morikawa, Y. Further Lowering of Work Function by Oxygen Adsorption on the K/Si (001) Surface. Phys. Rev. B 1995, 51, 14802. (67) Carter, E. A.; Ciccotti, G.; Hynes, J. T.; Kapral, R. Constrained Reaction Coordinate Dynamics for the Simulation of Rare Events. Chem. Phys. Lett. 1989, 156, 472–477. (68) Sprik, M.; Ciccotti, G. Free Energy from Constrained Molecular Dynamics. J. Chem. Phys. 1998, 109, 7737–7744. (69) Jacucci, G.; Rahman, A. Comparing the Efficiency of Metropolis Monte Carlo and Molecular-Dynamics Methods for Configuration Space Sampling. II Nuovo Cimento D 1984, 4, 341–356. (70) Allen, M. P.; Tildesley, D. J. Computer Simulation of Liquids; Oxford University Press: New York, 1989. (71) Hartwig, J. F. Transition Metal Catalyzed Synthesis of Arylamines and Aryl Ethers from Aryl Halides and Triflates: Scope and Mechanism. Angew. Chem., Int. Ed. 1998, 37, 2046–2067. (72) Stambuli, J. P.; Bühl, M.; Hartwig, J. F. Synthesis, Characterization, and Reactivity of Monomeric, Arylpalladium Halide Complexes with a Hindered Phosphine as the Only Dative Ligand. J. Am. Chem. Soc. 2002, 124, 9346–9347. (73) Stambuli, J. P.; Incarvito, C. D.; Bühl, M.; Hartwig, J. F. Synthesis, Structure, Theoretical Studies, and Ligand Exchange Reactions of Monomeric, T-Shaped Arylpalladium(II) Halide Complexes with an Additional, Weak Agostic Interaction. J. Am. Chem. Soc. 2004, 126, 1184–1194. 24 ACS Paragon Plus Environment

Page 24 of 31

Page 25 of 31

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

(74) Yamashita, M.; Hartwig, J. F. Synthesis, Structure, and Reductive Elimination Chemistry of Three-Coordinate Arylpalladium Amido Complexes. J. Am. Chem. Soc. 2004, 126, 5344–5345. (75) Vikse, K.; Naka, T.; McIndoe, J. S.; Besora, M.; Maseras, F. Oxidative Additions of Aryl Halides to Palladium Proceed through the Monoligated Complex. ChemCatChem 2013, 5, 3604–3609. (76) Crabtree, R. H. The Organometallic Chemistry of the Transition Metals, 5th ed.; John Wiley & Sons: New York, 2009. (77) Henkelman, G.; Uberuaga, B. P.; Jónsson, H. A Climbing Image Nudged Elastic Band Method for Finding Saddle Points and Minimum Energy Paths. J. Chem. Phys. 2000, 113, 9901–9904. (78) Steinmetz, M.; Grimme, S. Benchmark Study of the Performance of Density Functional Theory for Bond Activations with (Ni, Pd)-Based Transition-Metal Catalysts. ChemistryOpen 2013, 2, 115–124.

25 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

Figure 1: Overall cycle of the Suzuki-Miyaura reaction. The red-circled PhBr and PhB(OH)2 are the two reactants. The blue-circled PhPh is the product. PhB(OH)2 becomes PhB(OH)− 3 in alkaline solutions.

26 ACS Paragon Plus Environment

Page 26 of 31

Page 27 of 31

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: Model system of the simulation. The periodic boundary conditions are imposed on the system. PhX(X=Br or Cl) is placed in the center; halogen (X=Cl, Br) anions are shown in pink, carbon (C) atoms are shown in green, oxygen (O) atoms are shown in red, and hydrogen (H) atoms are shown in white. PhX is surrounded by H2 O molecules. Pd is represented in gold. The number of water molecules in a unit cell is 59. The cell size is 12.38 Å × 12.38 Å × 12.38 Å.

Figure 3: Constrained coordinate used in the BME simulations. The distance between Pd and the midpoint between the halogen atom and C in the ipso-position of the benzene ring is denoted by RPd−BC .

27 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 4: Relative free energy profile and the number of coordinated H2 O to Pd for the oxidative addition step. The top and bottom panels show snapshots of the oxidative addition step, while the middle panel shows the free energy profile (upper graph) and the number of coordinated H2 O molecules to the Pd atom (lower graph). The solid red and blue lines in the upper graph represent the free energy profiles of the oxidative addition of PhBr and PhCl to Pd in the water solvent, respectively. The solid green line represents the free energy profile of the oxidative addition step of PhBr to mono-ligated PdPPh3 in a vacuum. The vertical bar is the standard deviation of the estimated free energies. The red dotted and blue dotted lines represent the number of coordinated H2 O molecules to the Pd atom during the oxidative addition of PhBr and PhCl, respectively.

28 ACS Paragon Plus Environment

Page 28 of 31

Page 29 of 31

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 5: Energy profiles of the oxidative addition step. Upper left panel: static reaction path calculations using a climbing image nudged elastic band (CI-NEB) method for the oxidative addition of PhX (X=Cl, Br) to a PdX− anion complex or to a Pd(0) H2 O neutral complex in a vacuum. Upper right panel: free-energy profiles calculated on the basis of firstprinciples molecular dynamics simulations with BME sampling for the oxidative addition of PhX (X=Cl, Br) to the PdX− anion complex or to a Pd(0) neutral atom in water. Lower panels: Snap shot of the oxidative addition step of PhCl to the PdX− anion complex in water.

29 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 6: Trajectories of the distance between one Pd and another Pd during finite temperature MD for over 8 ps (upper-right panel) and final configurations. The black (red) line shows the Pd-Pd distances in case of PdBr− (Pd) before the oxidative addition (O.A.). The green (blue) line shows the Pd-Pd distances in case of PdBr− (Pd) after the oxidative addition. The black (red) framed panel represents the final structure in case of PdBr− (Pd) before the oxidative addition. The green (blue) framed panel represents the final structure in case of PdBr− (Pd) after the oxidative addition.

Figure 7: Relative free energy for the oxidative addition step of PhCl to Pd in the water solvent calculated using two reaction coordinates, namely the distance between Pd and the bond center of C1 of the phenyl ring and a Cl atom (RPd−BC ) and the angle between Pd-Cl and Pd-C1 bonds (θClPdC ) (left panel). Profiles of the reaction coordinate RPd−BC and the length of the Cl-C1 bond (right panel)

30 ACS Paragon Plus Environment

Page 30 of 31

Page 31 of 31

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

TOC graphics

31 ACS Paragon Plus Environment