Computed Ligand Effects on the Oxidative Addition of Phenyl Halides to Phosphine Supported Palladium ( 0 ) Catalysts

Table of


Introduction
Palladium-catalyzed cross-coupling reactions have been developed to accommodate a variety of substrates, 1 permitting the formation of new C-C or C-Y bonds, where Y is a heteroatom.The importance of this group of reactions was internationally recognized in 2010, when the Nobel prize for Chemistry was jointly awarded to Heck, 2 Negishi 3 and Suzuki 4 for their research on C-C bond formation.Other researchers who have also contributed to this field include Kumada, 5 Stille, 6 and Hiyama, 7 as well as Buchwald 8 and Hartwig, 9 who have made significant contributions to a range of carbon-heteroatom bond forming reactions.Computational studies of C-C cross-coupling reactions have recently been reviewed, 10 as have studies combining experimental and computational data to gain insight into palladium catalysis. 11me 1 Oxidative addition step for an aryl halide (ArX) with a palladium(0) catalyst A key step of these reactions, in common with some palladium-catalyzed C-H activation reactions, 12 is the oxidative addition of an aryl or alkyl halide to the palladium centre (Scheme 1).The activation of the Pd(0) catalyst, followed by insertion of the metal atom into the C ipso and halide (X) (or triflate, SO 3 CF 3 ) bond, has been shown to be rate limiting in certain conditions. 13roadly, catalyst design to support oxidative addition can be formulated as requiring an electron rich palladium centre, which supports the increase in oxidation state from Pd(0) to Pd(II). 14he catalyst also needs to tolerate an increase in coordination number during the oxidative addition by starting out as lowcoordinate or through the facile loss of one/several ligands. 15hermore, the system must remain active through multiple cycles, either by avoiding resting states and inactive reservoir species, 16 or by re-entering the catalytic cycle relatively easily.Ligands are often key to the fine-tuning of the activity and selectivity of organometallic catalysts.In general, electron-50 donating spectator ligands are favoured for oxidative additions, 17 and steric bulk can be used to support lower coordination numbers for the initial [PdL n ] species; the ligands shown in Scheme 2 fulfill these criteria and patented systems for crosscoupling reactions have been reviewed. 17 With a view to exploring and enhancing the role computational chemistry can play in the process of ligand-driven catalyst discovery and design we have been pursuing a computational methodology for the analysis and prediction of ligand effects in homogenous organometallic catalysis, combining DFT-calculated 65 ligand property parameters 22 with the analysis and prediction of both experimental 23 and calculated 24 data capturing catalyst performance, such as yield, rates and barriers to reaction.Accurate computational studies of catalytic cycles, especially where a manifold of competing reaction pathways exists, are an 70 important component of this methodology and we have evaluated the impact of computational method effects to establish and P t Bu 3 ligand because its large steric bulk restricts the number of accessible reaction pathways in the mechanistic manifold and limits the number of complexes and conformers which need to be considered.
For other, slightly smaller ligands such as PPh 3 and PCy 3 , a 10 number of additional steps/pathways may become energetically accessible and hence need to be taken into account (Scheme 3), making the study of these ligands computationally more demanding.In order to permit energetic comparison of competing pathways and comparison with experiment, the computational approach used must be sufficiently accurate, which requires it to take account of solvation, 26 dispersion 27 and vibrational corrections. 28he exact energetic balance of such competing pathways will be influenced most profoundly by both the ligand and the 20 aryl/alkyl halide used.Here we report calculated results for the oxidative addition of two aryl halides (PhBr and PhCl) to a palladium centre supported by a number of synthetically useful ligands (PCy 3 , PPh 3 , P t Bu 3 and the SPhos ligand as a representative of the family of biarly ligands developed in 25 Buchwald's group 21 ).Results for the PCy 3 ligand and a phenyl bromide substrate will be used to map out competing dissociative and associative pathways; this ligand has been used experimentally in a range of cross-coupling reactions with less activated aryl chloride and aryl triflate substrates 29 and has been the subject of several experimental kinetic studies, 13a,30 providing experimental data to validate calculated results.Ligand effects on the mechanism have also been considered for additional, synthetically relevant ligands (PPh 3 , SPhos, P t Bu 3 ) and the interplay between ligand and substrate effects has been evaluated by comparing data for PhBr and PhCl substrates in the presence of these ligands.These results have been used to illustrate what can be achieved with the methodology used, but also to discuss where computational improvements might be made in the future.

40
As indicated in Scheme 3, a number of different pathways and species can be envisaged as part of the mechanistic manifold for palladium-catalyzed oxidative addition.Depending on ligand size and donor strength, as well as the catalyst precursor used, the initial coordination number of the palladium(0) species [PdL n ] 45 can take values of n = 1-4 (1, 2, 11 or 12), although the monoligated complex 1 has not been detected experimentally for any ligand in solution; 31 solvent coordination might occur in this case.13b,31a,32 For large ligands it has often been assumed that the active  7); this will be denoted as the bisphosphine pathway (Path C, Scheme 3,   80   shown in blue).A further alternative, where solvent is weakly coordinated by [PdL] 1, 13b,31a,32 has also been considered for L = P t Bu 3 (section SI2 †).
Depending on the pathway followed, a range of different complexes can result from the oxidative addition.Monoligation 85 during the oxidative addition transition state, accessed by pathways A and B, will initially produce three possible isomers of a T-shaped complex, [Pd(L)(Ar)(X)] 5.The isomer with the phosphine ligand trans to the halide has been observed crystallographically for P t Bu 3 , whose steric bulk and ability to 90 form γ-agostic interactions with the metal can protect the adjacent vacant coordination site on the metal. 15,37For smaller ligands, the monoligated product 5 can undergo ligand addition to form square-planar cis or trans isomers of [PdL 2 (Ar)(X)], c-9 or t-9; both isomers have been observed crystallographically (see, for ).The cis isomer c-9 can also be reached directly by oxidative addition to a bisligated palladium center (Path C), which may be followed by isomerisation to give t-9.
Alternatively, 5 can coordinate a second T-shaped complex to form a halide bridged dinuclear complex [(µ-X) 2 Pd 2 (L) 2 (Ph) 2 ] 10, which again adopts a square-planar geometry around each metal center for either an anti or a syn dimer (a-10 and s-10 respectively), observed crystallographically for P t Bu 3 and P(1-Ad) t Bu 2 , 13a as well as P(o-tolyl) 3 19a,b and CataCXiumA. 39Based 5 on a computational study by Lledós, Espinet and co-workers, albeit without consideration of dispersion corrections, the observed higher activity of bulky ligands may be related to avoiding the formation of side products and reservoir species such as 9 and 10. 40

Method Effects on Prediction
Before describing our results on the different mechanisms and ligands, we first consider how the results depend on the computational protocol, for the case of speciation in solution between the different [Pd(PCy 3 ) n ] complexes (n = 1 -3).As in our previous work, 25b the geometry of all species has been optimized using the B3LYP functional in vacuum, with a medium basis set (denoted as BS1), followed by single-point energy calculations of the dispersion correction and with larger basis sets (BS2).(See computational details below and the ESI †.)The effect of geometry optimisation with dispersion-corrected functionals was found to be generally modest, but might make a more important contribution for crowded complexes involved in Path C, see discussion below and in section SI3 †.We have then computed vibrational frequencies and used ideal gas statistical mechanics to compute relative free energies in vacuum.These have been corrected by solvation free energies computed using continuum methods. 41e note that some authors have argued that this approach exaggerates entropic effects, which are smaller in solution than in the gas phase. 42However, such effects should in principle be treated adequately by the continuum model.In our experience, discrepancies with experiment observed when using this approach, instead of being due to such entropic factors, can be the result of inaccurate electronic structure theory.For relative energies, this can arise, for example, from the neglect of dispersion interactions in many density functional theory methods, underestimating the stabilization of higher coordination numbers.We note also that for complex systems such as those considered here, it is essential to identify correctly the lowestenergy isomer and conformer of each of the minima and 50 transition states involved. 43Finally, while anionic metal complexes have been observed to play a role in Heck and crosscoupling reactions performed in polar solvents and with some palladium precursors (Pd(OAc) 2 , Pd(dba) 2 ), 44 here we focused on reactions in toluene, where these species are less likely to play a 55 dominant role.30b The same applies to mechanisms for oxidative addition of alkyl halides, RX, in which the metal centre attacks the carbon atom to displace X -, which only subsequently adds to the metal.While there is evidence for such mechanisms with RX and in more polar conditions, they should play little role here.60   Table 1a and Fig. 1 show calculated relative free energies of different coordination numbers for the three [Pd(PCy 3 ) n ] species (n = 1-3), at different levels of theory.‡ Experimentally, solutions are found to contain a mixture of the n = 2 and n = 3 species.Measurement of the equilibrium constant at temperatures 65 between -68 and -85 °C yields an experimental enthalpy for binding of L to PdL 2 in solution of about -5 kcal mol -1 .30a This is not readily comparable to computation as it includes solvent enthalpic effects, which are not easily computed using continuum solvent methods.The standard Gibbs energy for binding is 70 negative at lower temperatures, but increases at higher temperatures.Extrapolation of these experimental data suggests a ∆G° of binding of -0.7 kcal mol -1 at 25 °C, and of +0.3 kcal mol - 1 at 100 °C.
It can be noted that at the B3LYP level of theory, PdL 3 lies 75 higher in energy than PdL 2 + L, and much higher in Gibbs energy as ligand dissociation is entropically favourable.This clearly disagrees with the experimental values.This has also been observed for L = PPh 3 by Ahlquist and Norrby. 35 Fig. 1 Method effects on the coordination number of [PdLn] for L = PCy3 at 90 °C (all energies are given in kcal mol -1 and are calculated relative to [PdL2], see also Table 1a).
On the other hand, upon including an approximate treatment of dispersion (∆E-D(BS2)), the [PdL 3 ] species is found to lie much 85 lower in potential energy compared to [PdL 2 ] and L in vacuum.The computed Gibbs energy of binding is now close to zero, and, in excellent agreement with experiment, is calculated to be small and negative near room temperature.The agreement with experiment is less good at higher and lower temperatures, perhaps because the temperature dependence of the solvation Gibbs energy has been neglected here.The continuum model is parameterized to give accurate results at room temperature and may be less reliable at higher and lower temperatures.Nevertheless, the present tests suggest that our chosen "best" computational protocol yields results within a few kcal mol -1 of experiment.In the rest of the paper, only these '"best"' computed standard free energies will be discussed, with all results corrected for a reaction temperature of 90 °C.A breakdown of energy contributions for each ligand may be found Tables S2-5 †.

Ligand Effects
Table 2 shows results for ligands discussed here; PCy 3 , PPh 3 , and SPhos, with previously published 25b and additional data for P t Bu 3 also included for comparison.Both PCy

a) PCy 3
As discussed in the previous section, at our "best" level of theory, the [Pd(PCy 3 ) 2 ] complex 2 Cy is slightly favoured over the trisligated complex 11 Cy at 90 °C, while trisligation becomes more favourable at lower temperatures (see also Table 1b).Attempts to optimize a tetrakis [PdL 4 ] complex for this ligand were unsuccessful and led to ligand dissociation.

Dissociative Pathway (A)
These results suggest that ligand dissociation to form the lowcoordinate [PdL] complex 1 Cy is energetically accessible (Figure 2, labeled in grey, Tables 2 and S2 †).This might exist as a transient species or as a solvated intermediate with a solvent molecule coordinated to the palladium centre.13b,30c,31a Note that the inclusion of continuum-based solvation free energies in our "best" computed energies should implicitly describe such solvent binding, at least roughly.We have argued previously 25b that the continuum solvent method may slightly overestimate the solvent stabilization of this highly unsaturated metal centre.In contrast, 50 several computational studies of other ligands have assumed that the active catalyst is the monophosphine complex 1, formed after ligand dissociation from bisligated [PdL 2 ], 2. 10a,33a,b,d,46 In addition to the dissociated complex 1 Cy shown, there will be a point of maximum Gibbs energy along the reaction path for addition of a ligand, or of the aryl halide, to [PdL].There is no potential energy barrier to such additions, with the Gibbs energy 65 barrier being due to loss of entropy upon approach.As in our previous work, 25b we estimate the Gibbs energy barrier to addition to be ca.4.5 kcal/mol at 90 °C, based on the rate constants for reactions known to be diffusion-controlled.With this estimate, the transition state (TS) for oxidative addition to the 70 monoligated metal centre [4 Cy ] ‡ is predicted to be very slightly lower in Gibbs energy (at 22.7 kcal mol -1 ) than the TSs for ligand loss from PdL 2 or for aryl halide addition to PdL (both at ca. 23.0 kcal mol -1 ).However, these values are very close, preventing calculations from distinguishing reliably between these options.75

Associative Displacement Pathway (B)
The highest Gibbs energy transition state on this pathway (B, Fig. 2) is for associative displacement of one PCy 3 ligand by the PhBr substrate ([6 Cy ] ‡ , 28.6 kcal mol -1 ), which has the highest barrier of 80 all the pathways considered here, presumably due to considerable ligand steric hindrance.This is one of the two mechanisms that are consistent with the experimentally observed order of reaction in Hartwig's kinetic study.13a It is also the mechanism predicted to be favoured with L = P t Bu 3 .25b However, for PCy 3 it can be 85 discounted based on the present calculations, since the barrier is so much higher than for Path C; this is also consistent with experiment as discussed below.

Bisphosphine Pathway (C)
The bisligated oxidative addition, [8 Cy ] ‡ (shown as the blue pathway in Fig. 2) has the lowest calculated barrier (20.6 kcal mol -1 ) of the competing transition states, [4 Cy ] ‡ , [6 Cy ] ‡ and [8 Cy ] ‡ studied here.This is in agreement with Hartwig's interpretation of the experimental kinetics, 13a whereby the irreversible step during the oxidative addition involves a bisligated complex.It is to be noted, though, that this observation is not enough to exclude the associative displacement mechanism discussed above.In experiments carried out in the presence of excess PCy 3 , the results suggest that the contribution from the monoligated Path A is suppressed, but the analysis is complicated due to the fact that some of the Pd(0) is now present as PdL 3 which is unlikely to participate in the reaction.
Hartwig et al. have also studied the kinetics for this reaction, at the slightly lower temperature of 10 °C, and with an excess of ligand present.They observe a slight decrease in reactivity as more L is added, due to slight changes in the position of the equilibrium with unreactive PdL 3 .Reactivity is first-order with respect to [PhBr], and, neglecting the complications arising from equilibrium with PdL 3 , their observed reactivity corresponds to a rate constant of 3 × 10 -4 M -1 s -1 , or to a ∆G ‡ of 21 kcal mol -1 .$ This is consistent with the other experiments, and with our computational results predicting the bisphosphine pathway (C) to be most favourable.The considerable sensitivity to method effects discussed earlier makes it difficult to distinguish confidently between paths A and C based on computation alone, and both pathways could indeed in principle be operating to some extent, though the bisligated path C is more consistent with experiment.
Computational work by Schoenebeck and Houk 29b has shown that, while calculations suggest the monoligated pathway (A) to be favoured in terms of Gibbs energies, experimentally observed selectivities can only be explained by considering a bisligated oxidative addition pathway.(Table S6 † gives an overview of the computational approaches used in computational studies referred to.)An earlier study involving the sterically and electronically similar 47 P i Pr 3 ligand that included solvation corrections, but not dispersion effects, had also ruled out the bisligated oxidative 60 addition (C), as well as associative displacement (path B), as too high in terms of Gibbs energy.33a The discussion of method effects above helps to explain these discrepancies, with the neglect of dispersion effects artificially favouring the dissociative pathway (A).65

Products of Oxidative Addition
The free energies of a range of monoligated and bisligated oxidative addition products are given in Table 2 (complexes 5, 9, 10).Isomer interconversion has been shown by Maseras to be 70 facile for three-and five-coordinate PH 3 complexes, which can be accessed by square-planar complexes if a ligand is dissociated in the former case or ligand/solvent is coordinated in the latter.42a Therefore, we have assumed here that interconversion can take place easily to reach the isomer of lowest Gibbs energy.In line 75 with experimental results, where the trans complex of 9 Cy has been observed as the only product of the oxidative addition, 15,30b this complex is the only one which could be optimised for the PCy 3 ligand.for some cross-coupling reactions. 4833a,d,35,40,43a,49 Fig. 3 summarises the calculation results for oxidative addition (see also Table 2 and Table S3 †).95   Complex 12 Ph , the tetrakis triphenylphosphine palladium(0) complex, has been crystallographically characterized 34b and is a commercially available precursor for PPh 3 -ligated catalysts.In solution, ligand dissociation to the [PdL 3 ] species 11 Ph is known to be near quantitative except at very low temperatures and in the 100 presence of excess phosphine. 36,50Hence, the equilibrium with 12 Ph , while it can be observed by 31 P-NMR at low temperatures, 50 is not relevant here.Dissociation of 11 Ph to form bisligated 2 Ph is found experimentally to be unfavourable, with a dissociation equilibrium constant of the order of 10 -4 M. 34a Our calculations underestimate the stability of 12 Ph significantly: this species is predicted to be much higher in Gibbs energy than 11 Ph + L. This is probably an artifact associated with carrying out geometry optimization at a level of theory that does not account for dispersion.Indeed, test optimisations with B3LYP-D2 (section SI4 †) show a significant change in geometry.This is also in agreement with previous work, 27,35 showing reasonable stability for 12 Ph when using methods that account for dispersion in optimization as well as energy calculations; for a more general overview, see e.g. a recent study by Jensen et al.. 51 It appears that this shortcoming in our calculations mainly affects the very crowded complexes, such as 12 Ph ; see the supporting information for a more complete discussion.The calculated ligand binding Gibbs energy required to transform 11 Ph into 2 Ph is 2.9 kcal mol -1 at 25 °C, corresponding to a dissociation equilibrium constant of 8 × 10 -3 .The experimental estimate for this dissociation constant places it as being significantly smaller than 10 -4 M. 34a Considering the three possible oxidative pathways starting from the [PdL 2 ], 2 Ph , complex, as with PCy 3 our calculations predict that the bisphosphine pathway C is favoured, proceeding via [8 Ph ] ‡ , with a barrier of 21.8 kcal mol -1 vs. 2 Ph , and a barrier of 23.7 kcal mol -1 vs. 11 Ph (90 °C, blue pathway in Fig. 3) to be favoured.The competing associative displacement route (B, red in Fig. 3) has a highest barrier ([6 Ph ] ‡ ) of 24.9 kcal mol -1 above 2 Ph , and the highest point of the monoligated oxidative addition pathway ([4 Ph ] ‡ ) lies 24.7 kcal mol -1 above reactants.The transbisligated product complex (t-9) is predicted to be most stable for this ligand and has indeed been observed crystallographically. 38b The bisphosphine pathway (C) was also found to be favoured by Kozuch and Martin. 16An earlier study by Fu, Liu et al. compared the three different pathways, but, in the absence of dispersion corrections, found both paths B and C to involve much higher Gibbs energy barriers than path A, in line with the analysis of method effects detailed above.33a Experimentally, kinetic data for oxidative addition of various aryl halides to solutions of 12 Ph is available for comparison.30c,34a,52 For the case of PhBr, an apparent rate constant k app = 9 × 10 -4 M -1 s -1 was measured at 25 °C for reaction with 12 Ph .30c In fact, under the conditions used, 12 Ph will have fully dissociated to 11 Ph , which will itself be in equilibrium with a small amount of 2 Ph .Assuming that the latter reacts with PhBr with a rate constant k, then 11 Ph will decay with k app = kK/[PPh 3 ].
Given that [PPh 3 ] was equal to 0.002 M in the experiments, this means that kK equals 1.8 × 10 -6 s -1 , equivalent to an activation Gibbs energy of 25.2 kcal mol -1 , in reasonable agreement with the value of 24.0 kcal mol -1 at 25 °C calculated for pathway C here.

c) SPhos
Developed in the late 1990s, 53 palladium ligated complexes of SPhos, PCy 2 (C 6 H 4 -2,6-(OMe) 2 -C 6 H 3 ) (Fig. 1), enhance reactivity for very hindered substrates in Suzuki cross-coupling reactions, 21a and in particular for aryl chlorides, 54 at low catalyst loadings.5] However, it became apparent that 60 palladium complexes can exhibit interactions between Pd and the second ring of the biaryl group which extends the catalyst's lifetime, 21a preventing cyclometallation and the formation of palladacycles.The steric bulk of alkyl substituents and additional groups on the biaryl unit also serves to increase the stability and 65 concentration of the monoligated Pd(L) species (1), considered key in the oxidative addition of aryl chlorides.21a Computational modelling of this ligand is challenging, as both the cyclohexyl and biaryl substituents can adopt a range of conformations and their motion is likely to be correlated.We 70 have found that the nature of the preferred conformer changes for the different species involved in oxidative addition, as shown in Table 3.As discussed in SI3 †, both biaryl and cyclohexyl group orientations respond to the coordination environment and these have been sampled extensively.While the effect of biaryl rotation 75 has been explored previously, 56 consideration of cyclohexyl conformational preferences was less complete and our present work has used a more extensive approach, combining database mining of the CSD 45 and DFT calculations to identify the lowest energy conformer for each complex (section SI3 †).Table 2 80 shows results for the conformer compatible with all steps, assuming that barriers to conformational change are lower in Gibbs energy than 'reactive' barriers lying along the oxidative addition mechanistic route.Low coordination numbers for this ligand are stabilized by the biaryl group coming into close proximity to the palladium centre, 56 either providing sites for secondary interactions with the 95 π-system (S a ) or the oxygens of the methoxy groups (S b ), as well as hampering the coordination of solvent/ligand/substrate through steric bulk.This is discussed in section SI3 †.
For SPhos the dissociative pathway (A, [4 SP ] ‡ , barrier = 19.6 kcal mol -1 ) presents the lowest barrier, while the associative 100 displacement pathway (B, [6 SP ] ‡ ) is significantly higher in Gibbs energy.We have not considered the bisligated pathway (path C) due to the steric hindrance exerted by two SPhos ligands in cis coordination sites.The overall calculated Gibbs energy of activation with this ligand is lower than with the other ligands 105 covered in this study, accounting in part for the success of this ligand in catalysis.To the best of our knowledge, no detailed kinetic studies have been reported for this ligand, preventing further validation against experimental data.
The bromide bridged dimer is quite similar in energy to the monoligated T-shaped product, and indeed a chloride-bridged dimer of SPhos has been isolated and characterized crystallographically. 57 This presumably also contributes to the synthetic utility of the ligand, by minimising the amount of metal complexes present as unreactive reservoir species under catalytic conditions.

d) P t Bu 3
29b,32b,33a,40,46,49,59 Our earlier study of oxidative addition 25b considered only this ligand and showed that computed activation parameters for the two monoligated mechanisms (pathways A and B) could account for the experimentally 20 observed kinetics in a near-quantitative way.25b Here we report additional results 25 relating to dimer formation to give 10 tBu 40 and the likely stability of a trans square planar product complex t-9 tBu , as shown in Table 2.We have also considered associative displacement of one ligand by solvent for this system and these results are included, along with barriers to isomerisation of the T-shaped product 5 tBu and the intermediate adduct 3 tBu , in section SI4 †.Calculation results have been summarized in Fig. 5 to allow comparison with data for other ligands.
For this ligand, the complex 2 tBu [Pd(P t Bu 3 ) 2 ] has been characterized crystallographically 60 and no higher coordination numbers have been observed by NMR.30a In addition, both the Tshaped oxidative addition product 5 tBu 15 and a chloride bridgedanalogue of 10 with an ortho-substituted aryl group 13a have been isolated and characterized crystallographically.The threecoordinate product 5 tBu is stabilized by agostic interactions 40 between the ligand and the palladium center, 15,61 which reduce the likelihood of dimerisation to form 10 tBu , 21c although the results here suggest that this may be finely balanced energetically, with 5 tBu and both isomers of 10 tBu lying within less than 2.5 kcal mol -1 of each other.Slight modifications of the 45 ligand or the substrate might thus lead to dimerisation and could trap some of the palladium in an unproductive pathway.The trans square-planar complex t-9 tBu lies higher in energy (10.0 kcal mol -1 ) than the other possible product species.
As we have reported previously, 25b mechanistically, this ligand 50 presents fewer pathways for evaluation, as the bisphosphine pathway (C) is not accessible due to the considerable steric hindrance of two PtBu 3 ligands, which would hamper their cis coordination.The preferred mechanism with PhBr involves passing over the associative displaement TS (path B) followed by 55 oxidative addition.The key TS lies higher in Gibbs energy ([6 tBu ] ‡ , barrier = 29.2kcal mol -1 ) than the monoligated oxidative addition (path A, [4 tBu ] ‡ , barrier = 26.9kcal mol -1 ).Our best estimate places formation of PdL at a Gibbs energy of around 30 kcal mol -1 , 25b which, with the variational TSs for addition of L or 60 PhBr, means that path B, i.e. the associative displacement route, should be favoured over a dissociative pathway (A) for this ligand.The calculated barrier to associative displacement is in reasonable agreement with available experimental data (see reference 25b for a more detailed discussion).
Before considering substrate effects in greater detail, we can take stock of ligand effects on the likely pathway and energetics of oxidative addition.The calculations reported here, validated where possible by comparison with available experimental data, suggest that very bulky ligands (P t Bu 3 and SPhos) favour low-75 coordinate pathways (A and B).Also, with these ligands, the Tshaped three-coordinate oxidative addition products are reasonably stable towards formation of dimer or bisligated complexes.This is presumably favourable for the reactions following on from oxidative addition in typical catalytic cycles.80   Smaller ligands, which are less "privileged" in experimental usage, can support higher coordination numbers around the metal centre.This means that prior to oxidative addition, ligand dissociation must occur.Nevertheless, contrary to the more bulky ligands where the C-X bond cleavage event occurs with only a 85 single phosphine coordinated to the palladium centre, for PCy 3 and PPh 3 , the bisphosphine pathway (C) is predicted to dominate.
Though the barriers to oxidative addition are similar to or indeed lower than those calculated for bulkier ligands, unproductive product complexes are energetically much more accessible after this step, suggesting that the rate-limiting step for cross-coupling reactions might occur later in the catalytic cycle.Our calculated barriers are generally in good agreement with experimental data where this is available, although we note that detailed kinetic studies of oxidative addition are rare.

Halide Effects
The cost of reagents for a chemical reaction can be an important factor, especially when considering the adoption of a new catalytic route.For oxidative addition, this can impact the choice of halide substrate used.Ideally, ligand design can enable the use of cheaper and structurally more versatile aryl chlorides as substrates. 62The geometries optimized with phenyl bromide provided a convenient starting point for expanding this study of ligand effects to also consider phenyl chloride as the substrate."Best method" Gibbs energies for 90 °C to allow comparison with the bromide results are shown in Table 4, with a more 20 detailed breakdown of method effects on energies for all ligands included in Tables S7-S10 †.
Experimental and computational studies suggest that the rate limiting step for phenyl chloride oxidative addition is the monophosphine transition state [4] ‡ (accessed by paths A and B), 25 with limited ligand effects on the likely mechanism for the systems considered.10a,13a In Hartwig's experimental kinetic study exploring both ligand and substrate effects, 13a different rate limiting steps were assigned for the different phenyl halides, with each rate law first order with respect to the concentration of 30 [PdL 2 ] species (2).Chloride substrates showed a dependence on the concentration of the ligand as well as the aryl chloride, thereby indicating that dissociation of a ligand must occur prior to rate-limiting ArCl oxidative addition ([4] ‡ ); this transition state contributes to both the dissociative pathway (A) and the associative displacement pathway (B) considered here.As discussed above, the trisligated complex 11 is unlikely to contribute to experimental observations at higher temperature.Comparison of the chloro species (Table 4) with the equivalent bromo complexes (Table 2) show that the former have higher relative free energies for transition states, and to some extent also for intermediates, relative to the bisligated reference complex, 2. The most pronounced substrate effects can be seen for the With additional results in hand, we can now consider the interplay between ligands and substrates more fully: In the case of PPh 3 no change in favoured mechanism or rate limiting step is 70 predicted by the calculations, although higher barriers suggest a slower rate of reaction, as expected for the chlorides.Similarly, for SPhos, there is no change for the different halide substrates, with the dissociative pathway (A) most likely.
Experimental data is available for the PCy For PCy 3 , there is the complication that the PdL 2 /PdL 3 equilibrium affects observed kinetics (as discussed in the section on method effects above).However, the available experimental data were obtained at fairly low concentration of free L, so that mostly PdL 2 should be present, and we will not consider this equilibrium here.In the study by Hartwig et al.,13a the measured rates again vary in proportion to substrate concentration, and inversely to concentration of free ligand.This again supports 10 mechanism A or B, with rate limiting oxidative addition through [4] ‡ ; the measured rate constant at 70 °C, k = 1.07 × 10 -6 s -1 , yields a ∆G ‡ = 29.4kcal/mol.In another experimental study, 30c the data appear to have been analysed using a rate law that assumes no dependence on concentration of free ligand, and without added ligand.The proposed rate constant of 0.015 M -1 s -1 at room temperature is at first sight much larger than that reported by Hartwig et al., 13a but considering the likely very low concentration of free ligand in this study, there is probably no inconsistency.We thus focus our comparison with computation 20 on the results of Hartwig et al. 13a As for PtBu 3 , our calculated transition state ([ Cl 4 Cy ] ‡ ) is found to lie at a relative Gibbs energy (31.1 kcal mol -1 ) that is consistent with experiment.However, our calculations show a significantly lower relative Gibbs energy for the bisligated transition state ([ Cl 8 Cy ] ‡ ), suggesting that path C should instead be favoured.This mechanism is not compatible with the observed kinetics, so it appears that the computational protocol used is not sufficiently accurate for this particular ligand/substrate combination, for reasons that are not yet entirely clear.There are of course many possible sources of inaccuracybasis set, functional, treatment of dispersion, statistical mechanics for entropy correction, and solvent treatment; some of these effects have been explored in greater detail (see SI4b-d †).As stated earlier, these effects seem to combine to yield errors of a few kcal mol -1 in most cases, but the present case is more 35 sensitive.

Conclusions
The manifold of reaction pathways for the oxidative addition of phenyl halide substrates to phosphine-modified palladium(II) complexes has been investigated with dispersion-corrected 40 density functional theory for a range of synthetically relevant ligands.Three mechanistic possibilities were considered (Scheme 3: A) a dissociative pathway accessing a low-coordinate [PdL] complex 1, which undergoes oxidative addition; B) concerted associative displacement to [PdL 2 ] 2, followed by monoligated oxidative addition; C) a bisligated pathway, where oxidative addition occurs directly to the [PdL 2 ] complex 2. Depending on ligand concentration, aryl halide substrate, solvent and temperature, different pathways can be accessed and in some cases, equilibria and competing pathways may need to be considered.Detailed experimental kinetic analysis can supply mechanistic insights and rate constants/barriers in these cases, and here we have been able to use such data to test and validate our computational methodology.
On the whole, the calculated barriers and favoured pathways 55 agree well with the available experimental data, allowing the computational prediction of likely reaction pathway (and hence a rate law), as well as the quantitative analysis of intermediates and transition states.√ In line with ligand design criteria derived from experimental studies, the bulky and electron-rich ligands P t Bu 3 60 and SPhos can access low-coordinate complexes, the most catalytically active species, easily throughout the cycle.While the bisphosphine oxidative addition step is reasonably facile for the smaller PCy 3 and PPh 3 ligands, ligand dissociation to access reactive palladium complexes becomes more important and the 65 catalyst is more likely to become trapped in unreactive intermediates.This work has demonstrated that a detailed evaluation of ligand effects is feasible and can support the interpretation of experimental data, allowing some pathways to be ruled out with 70 certainty.In addition, the energetic balance of competing reaction pathways has been shown to be quite subtle, illustrating and illuminating experimentally observed sensitivity to both ligand and substrate effects.Finally, this work has highlighted that multiple competing mechanisms may need to be considered for a 75 full evaluation of ligand effects, and that both Gibbs energy and dispersion corrections are necessary to achieve reasonable agreement with available experimental data.

Computational Details
Structures were fully optimized in Gaussian (G03, see ESI for full citation †) with the standard B3LYP density functional 63 and a flexible double-ξ (triple-ξ with ECP on Pd and Br) polarized basis set, denoted as BS1 (full details are given in the ESI †).B3LYP/BS1 harmonic frequencies were used to identify 85 stationary points and compute zero-point energy, enthalpic and entropic corrections.A polarizable continuum model was used to obtain single point solvation free energies with toluene solvent.Single point energy calculations were also carried out with a larger (augmented triple-ξ) basis set, again with ECP on Pd and 90 Br (denoted as BS2), and with the B3LYP-D2 functional 64 as implemented in ORCA. 65Full details of the computational methodology used, as well as further details about transition state scans and conformational searches, are given in the ESI †.
£ We have used Grimme's "D2" correction here to facilitate comparison with our PtBu3 results published previously (reference 25b).We note that the "D3" dispersion correction has since been described (described in S. Grimme, J.Antony, S. Ehrlich, H. Krieg, J. Chem.Phys.2010, 132, 154104), which Grimme recommends as superior.Our own test calculations (some of which are summarised in Table SI4f †) suggest that the dispersion correction does not substantially alter the trends and conclusions described here.§ Different steric measures can be considered, e.g.PPh3 has a Tolman cone angle 66 of 145 º and a He8_steric parameter 47 of 8.0 kcal mol -1 , whereas the corresponding data for PCy3 are 170 º, and 15.5 kcal mol -1 respectively.Both sets of data indicate that PCy3 is larger.$ They also studied the reaction with PhI, in the presence of excess L and at the low temperature of -80 °C.Under these conditions, the Pd will be present predominantly as PdL3, and accordingly inverse first-order 20 kinetics with respect to L are observed.The free energy of dissociation of L at this temperature is 2.3 kcal mol -1 .30a The observed rate-constant obtained by Hartwig et al, 13a is 8.5 × 10 -4 M -1 s -1 , corresponding to a ∆G ‡ of 13.8 kcal mol -1 .This would correspond to a ∆G ‡ with respect to PdL2 of 11.5 kcal mol -1 , which is lower than calculated here for PhBr -as expected for the more reactive iodide.√ With the exception of the PhCl/PCy3 system, where the calculated barrier for monoligated oxidative addition (path A) matches experimental data well, but calculations suggest that the bisligated oxidative addition of path C provides an alternative route with lower barriers.

a
3 and SPhos ligands can adopt a number of different conformers and these have been explored by mining of the Cambridge Structural Database (CSD) 45 and DFT calculations in different coordination environments, see section SI3 † for details.All energies shown are given relative to [PdL 2 ] + PhX and corrected for a reaction temperature of 90 °C, chosen to allow direct comparison with previous work.25b Table 2 B3LYP-D2/BS2 Gibbs energies in kcal mol -1 for the oxidative addition of PhBr to [PdL n ], see Scheme 3 for details and Tables S2-S5 for detailed energy contributions.† ‡ ∆G° (90 °C) + ∆Gsolv PCy3 PPh3 SPhos P t ½[(µ-Br)2Pd2L2(Ph)2] -14.1 -13.Optimisation unsuccessful; b not attempted for this ligand.

Fig. 2 60 (
Fig. 2 Potential Gibbs energy surface (∆G° (90 °C) + ∆Gsolv) for the 55 Mitchell and Baird carried out careful kinetics studies on the reaction of Pd(PCy 3 ) 2 with PhBr, both in the absence and presence of additional PCy 3 , at room temperature.30b Added PCy 3 inhibits oxidative addition, decreasing the relative amounts of [Pd(PCy 3 )] (1 Cy ) and [Pd(PCy 3 ) 2 ] (2 Cy ) due to formation of [Pd(PCy 3 ) 3 ] (11 Cy ).30bIn the absence of additional ligand, the Pd(0) is mostly present in the bisligated form.This is observed to disappear, with a pseudo-first order rate law, with an apparent rate constant given by k[PhBr] + k'.The first term here can be accounted for by the bisphosphine pathway (C).The value of k in their experiments is 1.3 × 10 -3 M -1 s -1 , corresponding to a Gibbs energy of activation of 21.3 kcal mol -1 .This compares very well to our calculated Gibbs energy of activation (22.8 kcal mol -1 at 10 °C, see TableS2b†).The second term can be accounted for by the dissociative pathway (A), with phosphine loss assumed to be rate-limiting, given the low concentration of free L and high concentration of PhBr in these experiments.The value of k' measured, 6.1 × 10 -4 s -1 , would correspond to a Gibbs energy of activation of 21.7 kcal mol -1 .Again, this matches our calculated barriers reasonably well (ca.24.5 kcal mol -1 at 10 °C from energy for 1 Cy + Gibbs energy barrier to ligand addition estimated as detailed above).

2
for labelling conventions used.
; b not attempted for this ligand; c looser convergence criteria used for BS2 single point calculation.45

50 monophosphine
and bisphosphine oxidative addition transition states, [4] ‡ and [8] ‡ respectively, as for these barriers different Chalide bonds are broken.For the monoligated oxidative addition, the increase in the barrier ranges from 2.2 -8.4 kcal mol -1 , with the largest difference observed for [4 Cy ] ‡ , whilst for the bisligated 55 pathway (C) both PPh 3 and PCy 3 see barrier increases of 2.9 and 3.0 kcal mol -1 respectively.The lower reactivity of aryl chlorides has been attributed to the strength of the C-Cl bond, 29a compared to weaker C-Br and C-I bonds.Our previous work 25b showed that for a P t Bu 3 catalyst the 60 change in rate limiting step is observed for the associative displacement step (path B), with the bromo system limited by the associative displacement of PhBr and a ligand ([ Br 6 tBu ] ‡ ), whilst the chloro system was restricted by the monophosphine oxidative addition of PhCl ([ Cl 4 tBu ] ‡ ), part of paths A and B. For this ligand, 65 dimerisation of the product becomes slightly less likely with the smaller chloride.
Entry for: Computed Ligand Effects on the Oxidative Addition of Phenyl Halides to Phosphine Supported Palladium(0) Catalysts Claire L. McMullin, a Natalie Fey* b and Jeremy N. Harvey* b Suggested graphic (Scheme 3 in manuscript): Suggested text: Reaction pathways for oxidative addition of phenyl halides to palladium complexes have been investigated computationally for L = PCy 3 , PPh 3 , PtBu 3 and SPhos.

Table 4
B3LYP-D2/BS2 Gibbs energies in kcal mol -1 for the oxidative 3and PCytBu 2 cases, We have not considered PCytBu 2 computationally, but, as in our previous work, 25b suggest that the calculated properties for P t Bu 3 should be rather similar.As seen in Table4, our calculations agree well with experiment, showing [4] ‡ as rate-limiting, and with a relative Gibbs energy of 30.7 kcal mol -1 .
This suggests that path A or B is followed, with insertion into the C-Cl bond rate-limiting.The measured rate constant suggests an 80 activation Gibbs energy of 29.3 kcal mol -1 .