Chuntian
Hu
a,
Kevin H.
Shaughnessy
b and
Ryan L.
Hartman
*a
aDepartment of Chemical and Biomolecular Engineering, New York University, Brooklyn, NY 11201, USA. E-mail: ryan.hartman@nyu.edu
bDepartment of Chemistry, The University of Alabama, Tuscaloosa, AL 35487, USA
First published on 4th November 2015
The influence of water on deprotonation and ionic mechanisms of a Heck alkynylation and its resultant E-factors were investigated. Estimation of the Hatta modulus, MH < 0.02, in cationic deprotonation, anionic deprotonation, and the ionic mechanism each separately confirmed an infinitely slow rate of reaction with respect to the diffusive flux within the thin film of the immiscible aqueous–organic interface. As a consequence, intrinsic kinetic expressions for far-equilibrium conditions were derived from first principles for each mechanism. Analyses of Gibbs free energies revealed that water potentially switched the rate-determining steps of cationic and anionic deprotonation to any of oxidative addition of organohalide to form Pd-complex (ΔG++ = 97.6 kJ mol−1), coordination of the alkyne with the oxidative addition adduct (ΔG++ = 97.6 kJ mol−1), or ligand substitution to form the cationic Pd-complex (ΔG++ = 94.9 kJ mol−1). Hydrogen-bonding in the transfer mechanism might account for the switch. Water, in general, was found to influence which step governs each catalytic cycle and the magnitude of its Gibbs free energy. Transformation of the synthesis from batch to continuous-flow was also studied by analyses of E-factors within the thin film. The amount of waste generated, as indicted by estimations of E-factors, was less in continuous-flow operation than in batch when the fastest step of deprotonation (ligand substitution) was infinitely fast with respect to the diffusive flux. The concentration of hydrophilic phosphine ligand was observed to influence mass transport limitations and the E-factor. Increasing ligand concentrations beyond (10.5), (13.3), and (23.2) × 10−3 mol L−1 for reaction temperatures of 353, 343, and 323 K increased the E-factor above its minimum value of 4.7, and it also induced mass-transfer-limitations. The switch from intrinsic to mass-transport-limited kinetics by finite changes in the ligand concentration explains ambiguity when performing aqueous-phase catalyzed Heck alkynylations and possibly multiphase Pd-catalyzed C–C cross-couplings in general. The potential exists to inadvertently mask the reactivity of useful ligands during discovery and to force mass transport limitations during manufacture. Understanding why the E-factor can be minimized is vital to the sustainable discovery and manufacture of fine chemicals, materials, natural products, and pharmaceuticals.
The field of homogeneous, organometallic catalysis has been widely investigated with much emphasis on the development of new chemical methods using transition metals (e.g., Pt, Pd, Rh, Ru) and for good reason. Albeit Pd, as an example, is scarce, its versatility in C–C and C–N bond formations (e.g., Buchwald–Hartwig amination3,4 and Heck,5–11 Hiyama,12,13 Negishi,14–19 Sonogashira,20–25 Stille,26–28 and Suzuki couplings10,29–33) continues to broadly impact society. For instance, the Pd-catalyzed Cu-free Sonogashira reaction (a.k.a., the Heck alkynylation) is the C–C coupling of an alkyne's terminal sp hybridized C with a sp2 C of an aryl or vinyl halide, or triflate (see Scheme 1). Heck alkynylation is widely used in the preparation of pharmaceuticals,34 natural products,35 and materials with unique optical properties36 (see Fig. 1).
![]() | ||
Fig. 1 Materials (e.g., carbonyl (2-ethyl-hexyloxy)-PPE 4),36 natural products (e.g., (+)-(S)-laudanosine 5, a kind of benzylisoquinoline alkaloids),35 and pharmaceuticals (e.g., Tazarotene (AGN-190168) 6, for treatment of acne and psoriasis)34 prepared by Heck alkynylations. |
Continuous production of such compounds introduces engineering challenges. Solids handling of inorganic halide salt by-products, and the separation and recycling of homogeneous catalysts from organic product streams are roadblocks.37 Catalyst recovery is of significant importance mainly because these transition metals (e.g., Pt, Pd, Rh, Ru) are rare in the planet. The design of water-soluble hydrophilic phosphines in biphasic systems introduces the opportunity for catalyst recovery, and it provides a reaction medium that solubilizes salt by-products. Water, however, has the potential to influence the reaction chemistry.
Two primary mechanisms that have been proposed for the Heck alkynylation are the carbopalladation38 and the deprotonation24 mechanisms. According to Ujaque, et al.'s density functional theory (DFT) calculations on a model system, the carbopalladation route for Heck alkynylations (in the absence of water) can be discarded.20 Mårtensson, et al. also demonstrated experimentally that the carbopalladation cycle is not operative in organic systems.23 However, the presence of water potentially switches the cross-coupling's mechanism from deprotonation to carbopalladation by reducing the Gibbs free energy barrier of the reductive elimination step in the carbopalladation catalytic cycle.39 In the present work, the influence of water on cationic and anionic derivatives of the deprotonation mechanism of a Heck alkynylation are investigated to further examine whether or not water actually switches the governing mechanism to carbopalladation.
For the deprotonation pathway, two alternative routes have been proposed: the cationic and anionic mechanisms23 shown in Fig. 2. The difference between the two variants of the deprotonation mechanism is the sequential order of the deprotonation steps. Both catalytic cycles, shown in Fig. 2a and b, undergo the reversible oxidative addition of the organohalide R1X to the Pd0Ln complex 7 (or 7′). The resulting oxidative addition adduct, species 8 (or 8′), then undergoes reversible complexation to the alkyne, forming an alkyne-Pd complex 9 (or 9′). At this point, the two mechanisms proceed differently to form complex 11 (or 11′). In the cationic mechanism (Fig. 2a), the ligand L is generally thought to substitute X in species 9 leading to the formation of cationic Pd complex 10, which undergoes the deprotonation of alkyne by the addition of base. By a comparison, with the anionic mechanism (Fig. 2b), the deprotonation of alkyne occurs first giving rise to the anionic complex 10′. The L–X ligand substitution then takes place to form complex 11′.
![]() | ||
Fig. 2 Proposed (a) cationic and (b) anionic mechanisms for the Pd-catalyzed Cu-free cross-coupling between a terminal alkyne and an aryl halide. |
Both the cationic and the anionic mechanisms finish by reductive elimination of the alkyne product to regenerate the Pd(0) catalyst 7 (or 7′). General understanding that water influences each step is ascertained in the present study through a comparison of DFT calculations20 with intrinsic reaction kinetics and the influence of mass transport limitations.
Solvent waste is a sustainability problem of synthetic chemistry that the engineering of continuous-flow syntheses could help address. The mass of the wastes generated relative to the mass of the desired products formed, the “E-factor”, typically ranges from 5 to 100 for fine chemicals and pharmaceuticals syntheses.40 Values of the E-factor are not independent of the reactor design where the potential for transport limitations exists. Calculations must be performed up front to avoid misinterpretations of mechanistic theory, overlooking novel chemical methods, and the development of falsified kinetics that increase E-factor values. This is especially important for multiphase reactions where interfacial contact area and film thickness of biphasic systems can differ in batch from flow. Much opportunity remains to minimize the mass of chemical waste generated during synthesis and production.
The present study reports the influence of water on deprotonation (both cationic and anionic) and ionic mechanisms of a Heck alkynylation. The potential to transform synthetic chemistry from batch processing to continuous-flow is also examined by estimations of E-factors. The optimum hydrophilic ligand concentration for different temperatures is presented. Analyses of Hatta moduli undergird that higher ligand concentrations could force mass transport limited conditions, which mask the discovery and the development of novel chemical methods. The discoveries reported in the present work could significantly reduce the mass of chemical wastes generated during fine chemicals and pharmaceuticals syntheses, both in the laboratory and production.
![]() | ||
Fig. 3 Gibbs free energy (ΔG++) profiles of the Heck alkynylation of an aryl iodide predicted by DFT calculations20 (dashed-line) and calculated using eqn (S1)–(S4) and (S7)–(S13)† applied to wet chemistry experiments (solid-line) for (a) cationic and (b) anionic deprotonation. |
The alkyne's substituent (R2 from Scheme 1) is expected to influence DFT calculated ΔG++ values by as much as 11.2 kJ mol−1 (for R2 = H, CF3, OMe, NMe2),20 which is reported by the error bars of Fig. 3. The nature of the phosphine also introduces conformational diversity from 1.3 to 47.7 kJ mol−1.44 Albeit the phosphine and the substitution of the aryl iodide may influence Gibbs free energy to some extent, our comparison of the DFT calculations20 with our data calculated (using eqn (S1)–(S13)†) from wet chemistry experiments using an aryl bromide provides fundamental insight on the participation of water, from step-to-step in the catalysis.
Water influences the reaction steps of cationic deprotonation. The Gibbs free energies of each step of cationic deprotonation are illustrated in Fig. 3a. As can be seen, the DFT calculations20 support that coordination of the alkyne 2 with adduct 8 exhibits the largest ΔG++ (at 96.2 kJ mol−1) of any reaction step, and thus the formation of adduct 9 in organic solvent is considered rate-determining. In the presence of water, however, any step of, oxidative addition of organohalide 1 to complex 7 (ΔG++ = 97.6 kJ mol−1), coordination of the alkyne 2 with adduct 8 (ΔG++ = 97.6 kJ mol−1), or ligand substitution to form the cationic Pd complex 10 (ΔG++ = 94.9 kJ mol−1) could be the rate-determining step. Water clearly influences which step governs the catalytic cycle.
An evaluation of anionic deprotonation yields similar observations. It is generally understood that coordination of the alkyne 2 with adduct 8′ is the rate-determining step with ΔG++ = 96.2 kJ mol−1, as shown in Fig. 3b. Hydrogen-bonding in the transfer mechanism might account for the switch to any step of, oxidative addition of organohalide 1 to complex 7′ (ΔG++ = 97.6 kJ mol−1), coordination of the alkyne 2 with adduct 8′ (ΔG++ = 97.6 kJ mol−1), or ligand substitution to form the Pd complex 11′ (ΔG++ = 94.9 kJ mol−1), as the rate-determining step. It is possible that H-bonding impedes either the substitution of the hydrophilic ligand or the oxidative addition of the alkyne 2, resulting in higher free energy barriers. The finite difference in free energy is an important consideration as water, an unconventional solvent, also solubilizes halide salt by-products that would otherwise impede performing the reaction in flow. The ionic strength of the aqueous phase, a thermodynamic consequence of introducing halide salt-product, could further influence the steps of either cationic or anionic deprotonation.
The first step, similar to the cationic and anionic mechanisms (Fig. 2), is the oxidative addition of the organohalide 1 to the complex 7′′. The cycle proceeds by the addition of base with the oxidative addition adduct 8′′. Here, it is possible that the corresponding phenylacetylide (i.e., R2––Pd) exists in solution as a result of excess base.20 Addition of this phenylacetylide to complex 9′′ forms Pd complex 10′′, which then undergoes reductive elimination to generate complex 11′′ and the alkyne product 3 (see Fig. 4). The ionic mechanism has also been studied by DFT20 in order to ascertain its theoretical likelihood as the dominant catalytic pathway of Heck alkynylation.
The Gibbs free energy profiles for the ionic mechanism are reported in Fig. 5 with and without any water. In either scenario, the first step, oxidative addition of the organohalide 1 to the complex 7′′, is rate-determining where in the presence of water ΔG++ = 97.6 kJ mol−1 and without aqueous solvent ΔG++ = 71.2 kJ mol−1 (as predicted by DFT calculations20). The obvious distinction between the two scenarios is shown by the ΔG++ values of the phenylacetylide addition to complex 9′′. In the purely organic system, no free energy barrier is reported. However, water appears to influence the free energy by 11.8 kJ mol−1. The observation undergirds that water molecules stabilize the ionic complex 9′′, potentially by H-bonding of the ligand. As we will next learn, the ligand plays a key role, not only in the cross-coupling's kinetics, but also in the generation of chemical waste.
![]() | ||
Fig. 5 Gibbs free energy (ΔG++) profiles of the Heck alkynylation of an aryl iodide predicted by DFT calculations20 (dashed-line) and calculated using eqn (S5)–(S13)† applied to wet chemistry experiments (solid-line) for the ionic mechanism. |
As expected in Fig. 6a, the E-factor of the Heck alkynylation decreases with increasing reaction temperature because the reaction accelerates (i.e., Arrhenius temperature dependence) and thermodynamic solubility tends to increase.
By its definition, water is excluded from the parameter. However, for the sake of a comparison, E-factor values calculated on the basis of water and organic solvent are also included in Fig. 6a. The larger E-factor values (where water is included) should not be interpreted as a system that will generate more waste. The opposite is true, case-by-case, depending on the chemistry, the number of total synthesis steps, and when the entire process is considered. Water, in this case of using hydrophilic ligand for Heck alkynylation, potentially minimizes the solvent mass required in downstream separation steps by two orders of magnitude.39 Merit exists to redefine the definition of the E-factor where water is an unconventional solvent in organic synthesis.
A general consideration in any chemical manufacturing strategy is the implementation of batch or continuous-flow operation. The E-factor can be applied to simplify the decision making process. Let us consider the thin film of thickness x0 at the aqueous–organic interface. The E-factor can be calculated based on the maximum attainable conversion within the thin film. As shown in Fig. 6b, E-factors were estimated to decrease with increasing kca values in the range of 1 to 100 min−1. Here, kc is the mass transfer coefficient and a is the liquid–liquid interfacial area. Values of kca are typically < 10 min−1 in batch (in the range of 2.5–5.2 min−1 in the present work) and from 50–100 min−1 for 0.5–1.0 mm I.D. capillary microchannels where segmented flow is achieved.45 In this example, the amount of waste generated is less in continuous-flow operation than in batch when the reaction rate is assumed to be infinitely fast with respect to the mass transfer rate.41
Let us now consider the case where the reaction rate within the thin film of thickness x0 (at the aqueous–organic interface) may or may not be infinitely fast with respect to the mass transfer rate.46,47 The Hatta modulus (MH) is a useful parameter that distinguishes reaction rates from diffusive fluxes of substrates, catalytic intermediates or products (see eqn (S14)†). As shown in Table 1, estimation of the specific reaction rates of the formation of complex 10 and 11′ (k3 and k4′ from Fig. 2) and of the molecular diffusivities41,48 enable estimation of the maximum MH values (see eqn (S14)†) where α = 1, β = 0, γ = 0, and δ = 0.
Entry | T (K) | x 0 (μm) | k i (M−1 s−1) | ![]() |
Diaryl alkyne yield (%) | M H (×103) |
---|---|---|---|---|---|---|
1 | 353 | 0.1 | 1373 | 7.2 | 37 ± 7 | 4.10 |
2 | 343 | 0.1 | 998 | 5.8 | 70 ± 10 | 3.86 |
3 | 323 | 0.1 | 497 | 3.8 | 81 ± 10 | 3.37 |
From Table 1, the maximum MH values of (4.10), (3.86), and (3.37) × 10−3 were calculated at 353, 343, and 323 K each separately for a film thickness on the order of 0.10 μm, which is the upper limit of the typical film thickness range for immiscible emulsions (in the range of 0.001 to 0.10 μm).49 Values of MH < 0.02 indicate an infinitely slow reaction,47 and consequently diffusion did not control the kinetics for all temperatures reported herein.
A trade-off exists between the chemical waste generated and the amount of ligand and precious metal consumed. On one hand, a greater mass of desired product is possible for a given residence time and with less waste generated by stoichiometrically increasing the ligand concentration with respect to the Pd loading (see eqn (S1–S6)†). On the other hand, accelerating the kinetics of this biphasic system any further by increasing the ligand loading has the potential to switch from intrinsic to mass-transport limited conditions. Fig. 7a illustrates the important point when considering either the cationic or the anionic deprotonation mechanism. Either deprotonation mechanism exhibits an overall lower free energy of product compared to the ionic mechanism (see Fig. 3 and 5). The fastest reaction step, whether the mechanism is cationic deprotonation or anionic deprotonation, is the ligand substitution. As a result, the E-factor is reported in Fig. 7a as a function of the maximum MH values for each temperature studied. One observes that the E-factor decreases with increasing MH values until the theoretical minimum is achieved (at 4.7) while varying the ligand concentration. Any excess ligand in the reaction mixture would, in principle, be considered waste because it would only increase the E-factor without any additional product generated. In this example, the maximum ligand concentrations are shown in Fig. 7a to be (10.5), (13.3), and (23.2) × 10−3 mol L−1 for 353, 343, and 323 K, respectively. It is also important to note in Fig. 7a that the theoretical minimum E-factor is attainable where the reaction is infinitely slow with respect to mass transfer (i.e., MH < 0.02).
Increasing the ligand concentration beyond the critical point where the E-factor is minimized has the potential to force mixed mass-transfer-reaction-rate and mass-transfer limitations. The existence of transport limitations can mask the discovery of useful ligands and actually contribute to the amount of waste generated at both the laboratory and production scales. Fig. 7b demonstrates that E-factor values climb (after its minimum has been reached) into the intermediate region (i.e., 0.02 < MH < 2 and mixed mass-transfer-reaction-rate limited)47 and beyond into the region where the reaction only occurs within the thin film (i.e., MH > 2 and mass-transfer limited).47 Note that the maximum conversion within the film (i.e., 100%) is achieved for values of MH < 0.02. The transition point depends mainly on the film thickness x0. In micro-scale segmented flow (e.g., Re < 1), for instance of two immiscible phases, the flow rate controls x0.50 Typical x0 values in micro-scale segmented flows range from 1 to 10 μm.45Fig. 8 illustrates the resultant E-factors at 353 K for MH values from 0.01 to 10.
![]() | ||
Fig. 8 E-factors at 353 K for MH values from 0.01 to 10 and aqueous–organic film thicknesses (x0) from 1 to 10 μm. |
One observes that values of x0 influence both the amount of waste generated and the range of mass transport limitations. Where the reaction is infinitely slow, the E-factor decreases with decreasing values of x0 whereas the opposite trend is observed for values of MH > 2. The same trend is shown in Fig. S2 and S3† for 343 and 323 K, respectively.
The existence of the ionic mechanism, in theory, could result in the diffusive flux being less than the reaction rate for one step of the mechanism. The fastest reaction step, the addition of the phenylacetylide to form Pd complex 10′′ of Fig. 4, exhibits a much lower ΔG++ value compared to any other step (see Fig. 5). The corresponding MH values, in the range of 24.6–39.6, are due to ki values on the order of 1012 min−1. However, the ionic mechanism depends on the existence of phenylacetylide in solution, which has yet to be experimentally validated. Furthermore, oxidative addition of the organohalide 1 to the complex 7′′, the rate-determining step of the ionic mechanism, predicts values of MH < 0.02. It is therefore reasonable to conclude that intrinsic kinetics govern the Heck alkynylation, water influences each of the three catalytic cycles, and that the ligand concentration plays a key role in the existence of mass transport limitations and in the generation of unnecessary chemical waste.
Analyses of Gibbs free energies revealed that water potentially switched the rate-determining steps of cationic and anionic deprotonation. According to DFT calculations, the coordination of alkyne substrate to form alkyne-Pd complex governs the overall reaction in organic solvent. In the presence of water, however, any step of, oxidative addition of organohalide to form Pd-complex, coordination of the alkyne with the oxidative addition adduct, or ligand substitution to form the cationic Pd-complex could be the rate-determining step. Hydrogen-bonding in the transfer mechanism might account for the switch. It is possible that H-bonding impedes either the substitution of the hydrophilic ligand or the oxidative addition of the alkyne in either mechanism, resulting in higher free energy barriers. Water also introduced a Gibbs free energy barrier for the addition of the phenylacetylide to form Pd-complex in the ionic mechanism. Water, in general, was found to influence which step governs each catalytic cycle and the magnitude of its Gibbs free energy.
Estimations of the E-factor of the thin film at aqueous–organic interfaces made progress toward understanding why the chemical waste generated during Heck alkynylations could be minimized. The amount of waste generated, in the present work, was less in continuous-flow operation than in batch when the fastest step of deprotonation in the Heck alkynylation (ligand substitution) was infinitely fast with respect to the diffusive flux. Fundamental understandings that minimize E-factor values could reduce operating costs and capital expenditures with the positive outcome of an overall greener planet.
The concentration of hydrophilic phosphine ligand was observed to play a key role in the existence of mass transport limitations and in the generation of unnecessary chemical waste (i.e., by approximation of the E-factor). Increasing ligand concentrations beyond minimum E-factor values induced mass-transfer-limitations, which could impact both laboratory- and production-scale syntheses. The switch from intrinsic to mass-transport-limited kinetics by finite changes in the ligand concentration explains ambiguity when performing Heck alkynylations and possibly multiphase Pd-catalyzed C–C cross-couplings in general. The potential exists to mask the reactivity of useful ligands during discovery.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c5re00034c |
This journal is © The Royal Society of Chemistry 2016 |