Joshua T. W.
Bray
a,
Mark J.
Ford
b,
Peter B.
Karadakov
a,
Adrian C.
Whitwood
a and
Ian J. S.
Fairlamb
*a
aDepartment of Chemistry, University of York, Heslington, York, YO10 5DD, UK. E-mail: ian.fairlamb@york.ac.uk
bBayer Aktiengesellschaft, Crop Science Division, Industriepark Höchst G836, 65926 Frankfurt am Main, Germany. E-mail: mark.ford@bayer.com
First published on 20th November 2018
Pd-Catalyzed arylcyanation reactions are widely applied in academic and industrial problems. We demonstrate for the first time that water acts as a powerful trigger for heterogeneous catalysis in arylcyanation reactions. A selection of trans-Pd(OAc)2(HNR2)2 precatalysts have been synthesized, characterized and kinetically scrutinized as precatalysts for the Pd-catalyzed arylcyanation of aryl bromides. Kinetic and mechanistic studies using in operando infrared spectroscopic analysis reveal the key role played by water, specifically deriving from K4[Fe(CN)6]·3H2O, in facilitating a remarkable switch in mechanism from a homogeneous to heterogeneous catalytic manifold. Evidence for heterogeneity includes sigmoidal kinetic traces, an order in Pd of 0.26 and catalyst turnover sequestration upon addition of excess Hg. The TOF and TON values recorded using K4[Fe(CN)6] (<220 ppm H2O) were higher than when using K4[Fe(CN)6]·3H2O (>4000 ppm H2O). Catalytic cycles are configured according to the experimental evidence gained from this study. We conclude that water plays a critical role in Pd catalyst speciation in arylcyanation chemistry.
Excess cyanide has been shown to give PdII cyanide complexes, e.g. Pd(CN)42−, as a dominant off-cycle catalytically inactive species6 and previous attempts to tackle this issue include the addition of diamines,7 Zn reductant8 or cyanide-shuttling agents.9 Further mechanistic work, using stoichiometric Pd and a bidentate phosphine ligand, has revealed the productive stage of a homogeneous catalytic cycle to be a migratory reductive elimination.10
In these reported systems, all are presumed to operate by a homogeneous catalytic cycle, i.e. involving single Pd atoms. The role of higher order Pd species (e.g. Pd clusters and nanoparticles) as catalytically competent species has received comparatively less attention, despite their confirmed role within the field of reductive Pd catalysis.11 The work of Frech et al.12 is an exception; they presented sigmoidal kinetic profiles characteristic of the presence of catalytically active Pd nanoparticles formed from a molecular PdII precursor by slow-nucleation and an autocatalytic growth mechanism.13–16
Water has previously been shown to play an integral role in Pd-catalyzed Suzuki–Miyaura reactions by aiding the leaching of Pd atoms from Pd nanoparticles (PdNPs) into the reaction mixture.17,18 It also accelerates the autocatalytic formation of Irn0 nanoparticles from polyoxoanion-supported IrI precatalysts in the hydrogenation of cyclohexene.19 Trace air plays a similar role20 and is also reported to assist in the activation of Pd(PtBu)2 to form catalytically-active PdNPs.21
In the context of arylcyanation, water aids the solubility of the cyanating agent, Zn(CN)2, changing the reaction from zeroth order in bromobenzene (anhydrous conditions) to first order in the presence of 1.7 vol% water (in DMF).22 However, despite there being a role for water, the use of either K4[Fe(CN)6] and K4[Fe(CN)6]·3H2O appears to be interchangeable throughout the literature, with little rationale for selection being made or a role for water being indicated for effective arylcyanations.
A kinetic study on the catalytic behaviour of Pd species, deriving from exemplar PdII precatalysts, in arylcyanation reactions is detailed within this paper. The use of in situ infrared spectroscopy in operando facilitated quantitative analysis of the rates of appropriate reactions, and the unprecedented, qualitative study of intermediate, solubilized K4[Fe(CN)6] concentrations throughout the course of arylcyanation reactions at Pd. The study led to the surprising finding that the presence of water triggers the reaction from being homogeneous to heterogeneous in terms of catalyst behaviour.
We determined that trans-Pd(OAc)2(pip)2 (pip = piperidine, 1-pip) is a competent PdII precatalyst for the direct C–H bond functionalization of 2′-deoxyadenosine using aryl halides.26 PdNPs readily derive from 1-pip under reducing reaction conditions (spherical PdNPs were found to be 1.5–2 nm in size26). Both Boykin and co-workers27 and Süss-Fink and co-workers28 showed that Suzuki–Miyaura cross-couplings mediated by similar PdII precatalysts, containing amine ligands, exhibited remarkable catalyst efficacy. The kinetic studies conducted by Süss-Fink and co-workers28 supported heterogeneous Pd catalysis, also implicating PdNPs. Conversely, trans-Pd(OAc)2(HNR2)2 complexes were found to be viable catalytic intermediates in selective oxidative C–H bond activation processes, e.g. functionalization of 2,2,6,6-tetramethylpiperidine (TMP).29
For this study, five well-defined mononuclear trans-Pd(OAc)2(L)2 complexes (L = HNR2), containing different sized secondary amine ligands, were synthesized by reaction of high purity Pd3(OAc)6 material (note: Pd(OAc)2 is usually a cluster of this type)30,31 with the prerequisite secondary amine ligands (Fig. 2). It was necessary to employ two equivalents of secondary amine; any less results in the formation of a dinuclear PdII complex containing only one amine ligand at each PdII centre. The full experimental characterization for these PdII complexes can be found in the experimental section (see ESI†), but several interesting structural features are worthy of mention. Inspection of the single crystal X-ray diffraction structures of 1-pyr, 1-dmp and 1-aze (Fig. 2) indicate that there is intramolecular H-bonding between the acetate ligand and the amino proton in both 1-pyr and 1-dmp, which is inkeeping the solid-state structure of both 1-pip26 and 1-tmp.29 Boykin conjectured that this type of interaction adds to the stability of similar complexes.27 In light of recent work on concerted metalation deprotonation processes,32,33 we propose that acetate-assisted deprotonation facilitates reductive elimination of the oxidized amino ligand, with concomitant formation of Pd0 (Fig. 3). The mechanism is consistent with what has previously been reported for 1-pip.24,26 Interestingly, there is no apparent H-bond in 1-aze, presumably due to increased bulk of the azepane ring. However, upon expanding the unit cell of the X-ray data, extended intermolecular H-bonding and π–π stacking between acetate carbonyl groups is apparent (see ESI†).
Fig. 3 Postulated mechanism for activation of 1-pip; we propose that the related trans-Pd(OAc)2(L)2 are activated in forming Pd0 species by a similar mechanism. |
N-Benzyl-4-bromo-6-methyl-2-pyridone, 2, was chosen as a benchmark substrate (Fig. 4). Confirmation of its aromatic character, by comparison with N-methyl-4-bromo-6-methyl-2-pyridone, benzene, pyridine and bromobenzene, was gained by computational (DFT) studies – see ESI.†
Fig. 4 Precatalyst activity screen. Comparison of precatalyst activity in the cyanation of N-benzyl-4-bromo-6-methyl-2-pyridinone, 2 (see chemical scheme). Reaction conditions: [2]0: 0.72 M, catalyst (0.01 mol%, 72 μM Pd), K4[Fe(CN)6] (0.22 eq.), Na2CO3 (1 eq.), DMAc, N2, 140 °C. Data cut off at 10 h. The Pd precatalyst structures are given in Fig. 2. |
The initial catalyst screening studies involved reaction of 2 mediated by anhydrous K4[Fe(CN)6]35 and confirmed the reaction as a binary system – no homocoupling or dehalogenation was observed under any conditions tested. Displacement of bromide by cyanide requires a Pd catalyst species. The catalytic activity of the Pd precatalyst series was examined using low loadings of Pd (0.01 mol%, 72 μM Pd in the system).
The kinetic profiles showing the formation of product 3 over time are presented in Fig. 4. The best precatalyst under the reaction conditions examined was found to be 1-aze, containing azepane ligands, or 1-tmp containing 2,2,6,6-tetramethylpiperidine ligands. The piperidine-derived precatalyst, 1-pip, exhibited comparable activity up to 50% conversion to product 3 while the pyrrolidine-derived precatalyst 1-pyr was considerably less active than either 1-aze or 1-pip. In all reactions, a first order consumption with respect to 2 was determined. By comparison with 1-pip, the cis-2,6-dimethylpiperidine-derived precatalyst 1-dmp, was poor, suffering from significant catalyst deactivation after 2 h (ca. 15% conversion).
In continuation of the kinetic studies, efforts were focused on 1-pip. As a precatalyst 1-pip: (a) exhibits very good activity for the reaction of 2 → 3; (b) is the most straight-forward to prepare within the series; (c) the post-catalysis PdNPs derived from it after the direct arylation of 2′-deoxyadenosine have been characterized by TEM as spherical PdNPs.24,26
The catalyst efficiency of 1-pip was optimum at increased concentration of Pd while maintaining constant catalyst loadings (at 0.01 and 0.005 mol%). At a Pd concentration of 72 μM (0.01 mol%) the reaction proceeds to 77% conversion to 3, after 6 hours at 140 °C (Fig. 5). This compares with 40% conversion at [Pd] = 36 μM (0.01 mol%). Similar catalyst turnover numbers (TONs)‡ were recorded for 0.01 mol% (72 μM) and 0.005 mol% (72 μM) catalyst loadings (7700 and 7600 respectively). Comparison of the catalyst turnover frequencies (TOFs)‡ revealed that catalyst efficiency was highest when the reaction was charged with 0.01 mol% 1-pip at [Pd] = 72 μM (Fig. 5) and a substrate (2) concentration of 0.72 M. This provided a TOF of 0.527 s−1 and TON of 7700.36 It was these optimum conditions under which the precatalyst series were screened (see Fig. 4).
The order with respect to Pd under pseudo-anhydrous conditions was determined using the Burés normalized time scale method (Fig. 6, top).37 The rate of reaction shows a second-order dependence upon Pd catalyst concentration, indicating a homogeneous reaction mechanism and rate-determining state involving two discrete catalyst molecules38,39 or one involving a dominant dinuclear [Pd2] complex.40 Given this, the mechanism must go via a rate determining step involving either bi-metallic Pd-to-Pd transmetallation (i.e. independent of [2]) or oxidative addition involving two atoms of Pd (i.e. dependent on [2]). The reaction was determined to be first order with respect to the concentration of 2 (Fig. 6, bottom), again using the Burés time normalized scale method.41 Therefore, due to the limited solubility of K4[Fe(CN)6] in the reaction solvent a rate-determining oxidative addition step involving two atoms of Pd is thus confirmed.
Fig. 6 Normalized time scale reaction profiles. (top) Analysis of two different concentrations of Pd using the Burés37 normalized time scale method, shows that under ‘low water’ conditions the arylcyanation reaction is second order in 1-pip. Equation used: [3] = t[Pd]n where n is the order with respect to [Pd] when the traces overlay. (bottom) By performing the cyanation reaction at two different concentrations of 2 under ‘low water’ conditions, the Burés normalized time scale method shows that the order of reaction is first order in 2.41 Equation used: , where a is the order with respect to [A] when the traces overlay. |
The non-identical curves for the reactions conducted in the presence of either K4[Fe(CN)6]·3H2O (black squares) and K4[Fe(CN)6] with additional water (3 equivalents, light green circles) is to be expected as added water is fully soluble (free to influence catalysis), whereas the water derived from K4[Fe(CN)6]·3H2O is not.
The observation of a sigmoidal kinetic curve for a metal-catalyzed reaction under reducing conditions is powerful evidence for the in situ formation of a heterogeneous catalyst.14 Extended induction periods indicate that mononuclear Pd (as 1-pip) is catalytically incompetent under the conditions used, requiring activation before productive catalysis can occur. The Finke–Watzky 2-step mechanism13 was developed precisely to elucidate the rate of the pseudoelementary steps involved in the formation of catalytically active metal nanoclusters from a catalytically inactive precursor by slow continuous nucleation (k1) and fast autocatalytic surface growth (k2).
Using the Finke–Watzky 2-step mechanism, it was possible to use a non-linear least squares fit to support our hypothesis of heterogeneous catalyst behaviour when using 1-pip together with K4[Fe(CN)6]·3H2O. The values of k1 and k2 for each catalyst loading could be measured by this approach (Fig. 8, left). Analysis of all the 1-pip catalyst loadings showed that a good fit was only obtained with a 1-pip catalyst loading below and including 0.5 mol% ([Pd] = 1.8 mM), indicating that precatalyst activation was occurring via the Finke–Watzky 2-step mechanism at lower catalyst loadings. Further information regarding the goodness of fit and subsequent error analysis using Fisher's F distribution is included in the ESI.†42 Furthermore when a large excess of metallic Hg was added to a reaction performed with 1-pip (0.5 mol%, 1.8 mM Pd) at 93 min, no further product 3 was formed (Fig. 8, right). The Hg drop test provides a means of showing that a metallic surface is playing some role in terms of productive catalysis;14 importantly, this test is only supportive evidence for the presence of heterogeneous catalyst species, and cannot be interpreted as conclusive proof for catalyst heterogeneity. Other powerful poisoning tests, such as the quantitative addition of 1,10-phenanthroline,43,44 were considered but deemed to be chemically incompatible with the precatalyst (1,10-phenantroline would outcompete piperidine as a ligand at PdII, resulting in a chelated N,N-Pd(OAc)2 species), inhibiting the catalyst activation mechanism proposed vide supra.
Plotting TOF vs. catalyst45 loading indicates a counter-intuitive, negative effect upon catalyst efficiency caused by increased catalyst loadings (Fig. 9). At low Pd catalyst loading and concentration (0.05 mol%, 180 μM), the catalyst turned-over steadily for ca. 12 hours until reaction completion, however when the Pd catalyst loading was increased to 0.1 mol% a higher TOF was achieved (Fig. 9, left). This indicates that in the presence of water the catalyst turns over most efficiently at a [Pd] of ca. 360 μM and Pd:substrate ratio of 1:1000. However, upon increasing the Pd catalyst loading and concentration further to 0.5 mol% (1.8 mM Pd) and above, the catalyst TOF suffers a significant decrease. The lower TOF resulting from the increased [Pd] can be attributed to fewer catalytically active sites per mole of Pd, indicating the presence and role for larger Pd clusters at higher [Pd]. Increased TON and TOF upon decreasing catalyst loading has also been attributed to catalyst degradation by aggregation in Sonogashira reactions: Fairlamb et al. report higher TOFs at 0.001 mol% than at 0.01 mol% (ref. 46) catalyst loadings and a significant decrease in catalyst activity above that catalyst concentration. A similar outcome noted by De Vries et al. where higher yields were obtained in a Heck47,48 reaction at 0.02 and 0.08 mol% Pd catalyst loadings than at 0.00125 and 1.28 mol%. When stabilized PdNPs of a range of sizes (1.8–4 nm) were utilized in a Suzuki reaction it was reported that TOF decreased upon increased NP size.49 An inverse TOF/[Pd] relationship therefore indicates that heterogeneous species are likely involved in the reaction. Extraction of the observed rate constants from peak rates of the sigmoidal reaction profiles allows for the calculation of the order with respect to Pd via the differential method. This revealed a reaction order of 0.257 (±0.02) (Fig. 9, right) and shows that 1 in every 3.89 Pd atoms are catalytically competent; providing further evidence for the involvement of small Pd clusters when K4[Fe(CN)6]·3H2O is used as the cyanating agent, i.e. under ‘excess water’ conditions.
It is worthy of note that the hexacyanoferrate and additive counterions were also found to have a significant effect on the efficiency of catalysis. Using kinetic methods, this has been extensively studied by following the characteristic [Fe(CN)6]4− infrared absorbance (seen at 2045 cm−1) in correlation with reaction progression. The additive, Na2CO3, has been shown to assist salt metathesis with K4[Fe(CN)6] which tensions dissociative ligand exchange reactivity and complex solubility. The use of alternative carbonate additives (K2CO3 and Cs2CO3) results in no catalytic turnover; this has been attributed to a differing hexacyanoferrate solubilities following partial counterion exchange. The direct addition of Na4[Fe(CN)6] as a cyanating agent resulted in inhibited catalytic turnover caused by its increased solubility in hot DMAc (shown by in situ IR spectroscopy). These findings are summarized in the context of a complete mechanism in Fig. 10, however further discussion, including kinetic traces, is included in the ESI.†
For example, when water is introduced into the arylcyanation reaction in the form of K4[Fe(CN)6]·3H2O, it triggers the reaction to proceed via heterogeneous catalytic manifold (H2OTOTAL > 4000 ppm, as measured using Karl-Fisher titration) while inhibiting homogeneous catalysis. The proposal is supported by an inverse TOF/catalysis loading relationship, a partial order in Pd (0.26), sigmoidal kinetics fitted to the Finke–Watzky mechanism, and a positive Hg drop test (testing for surface catalysis). The order in Pd suggests that 1 in 3.8 Pd atoms are catalytically active, which taken alongside the positive Hg drop test and sigmoidal kinetic behaviour, confirms the involvement of higher order Pd species in the catalysis. Water facilitates the leaching of atomic Pd from PdNPs.17 In our study, ‘excess water’ is the trigger for a catalytic cycle involving heterogeneous Pd species.
Under ‘low water’ conditions, using anhydrous K4[Fe(CN)6] as a cyanide source, kinetic behaviour representative of a homogeneous catalytic cycle has been evidenced (H2OTOTAL < 220 ppm). Under these pseudo-anhydrous conditions the rate of reaction shows a second order dependence upon 1-pip concentration, indicating that two Pd atoms, i.e. a dinuclear Pd2 species, is involved in the rate-determining transition state.40 Under pseudo-anhydrous conditions it is possible to achieve lower catalyst loadings (down to 0.005 mol%) and optimum conditions were found to be 0.01 mol% with a Pd concentration of 72 μM. Under such conditions 1-pip achieved a TOF of 0.527 s−1. The relative activity of Pd(OAc)2(HNR2)2 precatalysts was found to be in order: 1-tmp > 1-aze > 1-pip > 1-pyr > 1-dmp.
Lastly, further experimental evidence confirmed that sodiated additives assist the solubility of K4[Fe(CN)6] by partial counterion swapping to NaxKy[Fe(CN)6], where x + y = 4. This process facilitates the formation of NaBr (lattice enthalpy: 747 kJ mol−1) upon dissociative ligand exchange with Pd, rather than KBr (lattice enthalpy: 682 kJ mol−1). Other carbonate salts do not facilitate productive catalytic turnover, but significantly influence the concentration of [Fe(CN)6]4− observed in solution (see ESI†). Using Na4[Fe(CN)6] results in a significantly reduced reaction efficiency, i.e. low product conversions and slow reaction, due to increased solubility.
A proposed mechanistic rationale from the observations deduced from this study is depicted in Fig. 10 for 1-pip. Reduction of 1-pip occurs almost instantaneously (typically within seconds), delivering a catalytically active, bimetallic [LnPd0]2 species to the reaction system. In the presence of ‘low water’ (<220 ppm), [LnPd0]2 is then able to enter homogeneous catalytic cycle 1. In both cycles, the oxidative addition product then undergoes dissociative ligand exchange with NaK3[Fe(CN)6], which is formed through counterion exchange with Na2CO3 (not displayed in Fig. 10). Reductive elimination from LnPdII(Ar)(CN) yields the ArCN product and regenerates the active catalyst species. Under any of the reaction conditions tested, LnPd0 is in equilibrium with aggregated Pd particles, either as a reservoir, moribund form or catalytically-competent form. The ‘excess H2O’ shifts the position of the equilibrium from LnPd0 to Pdn0 (representative of the total active aggregated Pd), which allows cycle 2 to become operative, with concomitant shutdown of cycle 1. The dissociative ligand exchange step is undoubtedly complicated, as there are several species that can act as the cyanide source. The steps proposed within the catalytic manifold involve salt metathesis and eventually lead to NaK3[Fe(CN)5X], the fate of which is likely halogenated Fe species.
After backfilling the round bottomed flask with dry N2, a background spectrum was collected, and then dry, degassed DMAc added (6 mL). At this point a reference spectrum was recorded with stirring. After sequence initiation, the reaction was stirred at room temperature for 5 minutes, before being heated to reaction temperature (as monitored by the internal thermocouple) by being placed in a silicon oil bath set to an appropriate level above the desired reaction temperature. After the reaction mixture reached the correct temperature (determined using an external, digital thermometer), the solvent was stirred for a further five minutes and a spectrum used as a reference. To the hot solvent was then added the Na2CO3 (0.25–1 eq., freshly ground in a pestle and mortar) and was stirred for a further 5 minutes. After the system reached equilibrium, the cyanating agent (K4[Fe(CN)6]·3H2O or K4[Fe(CN)6], 0.167–0.5 eq., freshly ground in a pestle and mortar) was added and stirred for 5 minutes. A spectrum was then selected as a second reference. To the mixture was then added N-benzyl-4-bromo-6-methylpyridin-2-one (2, 1 eq.), and the sides of the flask rinsed with dry, degassed DMAc (0.5–1 mL, adding the appropriate amount to make the final total reaction volume up to 7 mL) via a septum. The reaction mixture was then stirred for 5 minutes, to allow for temperature equilibration back to the desired reaction temperature and dissolution of the substrate. During this time the appropriate amount of dry, degassed DMAc was added to the small Schlenk containing the precatalyst to make a stock solution. All catalyst stock solutions were made to a concentration of 25 mmol dm−3 prior to their addition to reaction mixtures. After the reaction mixture had achieved a constant temperature and immediately after complete dissolution of the catalyst, the reaction was initiated by addition of an appropriate quantity of catalyst stock solution using a syringe to the round bottomed flask via a septum.
The reaction progress was monitored by measuring the increase in absorbance at 2238 cm−1, corresponding to the product's nitrile stretch (in 3). After the absorbance at 2238 cm−1 had reached a constant value (indicating reaction completion) the experiment was stopped, and the reaction mixture allowed to cool to room temperature. Once the mixture had cooled, the final conversion of each reaction was measured by direct 1H NMR spectroscopic analysis of the crude reaction mixture. An aliquot of the reaction mixture was filtered through a small Celite/cotton wool plug to remove any organic-insoluble salts and diluted with CDCl3 before running the NMR experiment. The final conversion was calculated by measuring the ratio of the diagnostic peaks at 6.80 ppm for N-benzyl-4-bromo-6-methylpyridin-2-one (2) and 6.85 ppm for N-benzyl-4-cyano-6-methylpyridin-2-one (3).
All peak absorption data for 2238 cm−1, 2045 cm−1 and any other peaks of interest were exported into Excel in three forms: i) with no reference subtraction; ii) after subtraction of a reference spectrum of DMAc at reaction temperature; iii) after subtraction of a reference spectrum of the reaction mixture before substrate addition, at reaction temperature °C. All subsequent analysis was carried out on dataset iii). All kinetic traces (conversions) were normalised using the final reaction conversion (as calculated vide supra) and the absorption value of an appropriate data point at 2238 cm−1 (mid data scatter) after reaction completion. All data was smoothed with a moving four-point average. An annotated figure detailing the experimental set up is included in the ESI.†
Footnotes |
† Electronic supplementary information (ESI) available: Experimental and computational details, X-ray details, results of all kinetic studies and characterization data accompanies this paper. See DOI: 10.1039/c8re00178b |
‡ The methods used for determining turnover numbers (TONs)/turnover frequencies (TOFs) are detailed in the viewpoint paper reported by Kozuch and Martin.36 The TOF errors were calculated using the gradient of TON vs. time during the period where the reactions were at their optimum rates. |
This journal is © The Royal Society of Chemistry 2019 |