Bradley J.
Butler
and
Jason B.
Harper
*
School of Chemistry, University of New South Wales, UNSW Sydney 2052, Australia. E-mail: j.harper@unsw.edu.au
First published on 16th October 2014
The effect of the ionic liquid 1-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide on the ethanolysis of diethyl chlorophosphate was investigated. The addition of the salt increases the rate constant, though the extent of this increase was dependent on the mole fraction used and the organic co-solvent. Temperature-dependent kinetic analyses were used to demonstrate that the rate enhancement was as a result of a decrease in the activation enthalpy offsetting a decrease in the activation entropy. Altering the salt used showed the importance of the balance of these enthalpic and entropic effects.
There are a growing number of studies on how ionic liquids affect reaction outcomes, such as reaction rate and selectivity, for examples see the work of Welton,10–20 Chiappe15,21–31 and D'Anna32–40 along with that from our own group.2,5,7,41–54 These investigations have typically been limited, particularly for kinetic studies, to well-defined organic processes focusing on reactions at carbon centres such as the reaction of benzyl halides with pyridine44,49,51,52 and the methanolysis of 3-chloro-3,7-dimethyloctane.41,44,49,52 Thus, while there is a growing understanding as to the effect of ionic liquids on reactions at carbon centres, these is little in the way of investigation at alternate centres such as phosphorus. Reactions at phosphorus centres are important industrially in the preparation of plasticisers55 and pesticides,56 along with on the smaller scale such as the preparation of oligonucleotides.57
While there have been some studies into phosphorus chemistry in ionic liquids, they have been primarily qualitative in nature (i.e. no kinetic analyses undertaken), focusing on the stability of the phosphorus compounds in these solvents and the ability of ionic liquids to prevent unwanted side reactions such as hydrolysis. Initial work58 demonstrated the different activity of water in a range of ionic liquids and highlighted the importance of the nature of the water present. In a subsequent study,59 the same group considered the reaction of phosphorus trichloride with isopropyl amine in either one of a variety of hydrophobic ionic liquids or the molecular solvent hexane; when carried out in the ionic liquids there was no hydrolysis of phosphorous reagent (cf. a significant extent in the hexane case) and a greater extent of conversion to the desired product. Subsequently, the synthetic utility of this type of reaction control has been demonstrated.60,61
A study where kinetic data was obtained involved reaction of Paroxon with piperidine in a series of ionic liquids.62 It was noted that the ionic liquids affected both the rate constant of reaction of the phosphorus species and the product distribution; in the example used nucleophilic substitution at the phosphorus centre (SN2@P), an aliphatic carbon centre (SN2@C) and at an aromatic carbon centre (SNAr) were all possible. Whilst of interest in noting the different reactivity in ionic liquids relative to molecular solvents (and between ionic liquids), no activation parameter data was determined so the origins of the rate changes in terms of entropy and enthalpy of activation were not discussed. Further, no dependence of reaction outcome on the proportion of ionic liquid in the reaction mixture was considered; this has been shown to be significant in previous reports.5,41,42,51,54
The work described herein aims to consider the effects of an ionic liquid on the rate of a reaction at a phosphorus centre and to correlate this with what is known about how ionic liquids affect reactions at carbon centres.
The reaction chosen was the ethanolysis of diethyl chlorophosphate (Scheme 1), which proceeds through a bimolecular substitution mechanism.63 Diethyl chlorophosphate is a commonly used analogue for both pesticides and nerve agents due to it having comparable reactivity while being of a much low volatility64 and therefore safer to handle.65–67 The ionic liquid 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([bmim][(CF3SO2)2N]) was chosen as it is well studied and allows direct comparisons with previous work where this ionic liquid was also used.7,41,42,44,51,53
![]() | ||
Scheme 1 Ethanolysis of diethyl chlorophosphate by ethanol-d6 to afford the partially deuterated triethylphosphate. |
Kinetic analyses were carried out in solutions containing diethyl chlorophosphate (ca. 0.3 mol L−1) at a given temperature and at a given mole fraction of salt. For the cases where the co-solvent was ethanol-d6 it also served as reagent. In cases where the co-solvent was acetonitrile, ethanol-d6 was also present at a concentration of ca. 1.0 mol L−1. The reaction was followed using 1H NMR spectroscopy‡ with spectra recorded at regular intervals with a minimum of 100 spectra recorded to a minimum of one half life. For cases where the observed rate constant was <8 × 10−6 s−1in situ monitoring was impractical and the sample was held at the desired temperature in a water bath. Spectra were acquired periodically and initial rates determined up to ca. 10% conversion.
For cases where a large excess of nucleophile was used, the extent of conversion was determined through monitoring the change in the integration either of the signal at δ ca. 4.3 due to the CH2 protons of the diethyl chlorophosphate or of the signal at δ ca. 4.1 corresponding to the CH2 protons of the resultant triethylphosphate. From these integrations, pseudo-first order rate constants and further, the corresponding second order rate constants could be calculated. In the case of second order reactions, the relative concentrations of the diethyl chlorophosphate and the triethylphosphate were determined through comparison of signal at δ ca. 4.3 allowing the absolute concentrations to be calculated from the initial concentration of the reagent; from these and the concentration of the nucleophile, the second order rate constant was determined. The bimolecular form of the Eyring equation70 was utilised to calculate the activation parameters.
There are potential competing reactions to that shown in Scheme 1; both hydrolysis due to adventitious water and nucleophilic substitution at a methylene group (SN2@C) are possible. The former process would generate the diester diethyl hydrogen phosphate, while the latter would generate a phosphochloridate that may react further with ethanol to generate the same diester. Under all of the reaction conditions described here only a single byproduct (identified through spiking experiments to be the diester) was observed, and in negligible amounts (with the exception of the chloride salt, discussed further below), using either 1H or 31P NMR spectroscopy (for representative spectra, see ESI†).
As represented in Fig. 1 there is a clear trend in the rate constant for the reaction shown in Scheme 1 as the mole fraction of the ionic liquid in the reaction mixture changes. Initially, the rate constant for the process increases up to a maximum value at ca. χIL 0.3, corresponding to a ca. four-fold (Table 1) increase from the case in ethanol (χIL 0). As the mole fraction increases further, there is a decrease in the rate constant up until the maximum mole fraction used (χIL 0.72); at this point there is still an increase in rate constant compared to the molecular solvent case (of ca. two-fold).
This trend does not exactly parallel those trends observed previously for reactions at carbon centres. Arguably the closer comparison mechanistically – the reaction of benzyl halides with pyridine in acetonitrile proceeding through a bimolecular mechanism – showed an increase in the rate constant with proportion of ionic liquid concentration to a maximum at χIL 0.77 followed by a slight decrease to χIL 0.85.51 The extent of rate enhancement was similar to that described here, with the maximum enhancement being ca. 4.5 times.
The alternate comparison is with the SN1 methanolysis of 3-chloro-3,7-dimethyloctane in mixtures of methanol and [bmim][(CF3SO2)2N].41 In that case, an initial increase in rate constant with proportion of ionic liquid, followed by a decrease as more ionic liquid was added, was also observed. In that case, however, the maximum effect occurred at a much lower proportion of ionic liquid – a relative rate enhancement of ca. two-fold at χIL 0.02.
It is interesting to note that this trend of rate constant enhancement has similar elements to those previous studies, particularly an increase in rate constant to a given mole fraction of ionic liquid followed by a decrease. In this case, the concentration at which the maximum effect occurs is different to those described previously, lying approximately halfway between the two cases. This suggests that the mechanism of reaction is not the only variable that must be considered when analysing the effect of ionic liquid concentration on the rate of a particular reaction.
Unlike the previous case, which was carried out under pseudo-first order conditions, the low concentration of ethanol used meant that kinetic analysis had to take into account the change in concentration of this reagent. As ethanol-d6 is NMR silent the change in concentration of this reagent had to be calculated utilising the known starting concentration and using the assumption that any absolute change in concentration of the diethyl chlorophosphate was equivalent to the absolute change in concentration of the ethanol-d6. This, and the need to determine the concentrations of multiple species, complicated analysis in this case.
As seen in Fig. 2, the general trend of an increase in rate constant to a maximum at χILca. 0.3 followed by a decrease, mirrors that seen in the previous case of [bmim][(CF3SO2)2N] in ethanol-d6 (Fig. 1). However the actual extent of rate constant enhancement is greater (Table 2), with a ca. sixteen-fold increase. This is primarily due to the rate constant in the χIL 0 case being markedly lower in acetonitrile ((1.4 ± 0.2) × 10−6 mol L−1 s−1) than in ethanol-d6 ((7.5 ± 0.3) × 10−6 mol L−1 s−1). This is most likely due to the relative natures of ethanol and acetonitrile with a larger rate constant in a protic solvent likely due to more effective stabilisation of the charge-separated transition state. The most interesting observation here is that the effect from the presence of the [bmim][(CF3SO2)2N] is much greater than any affect caused by changing the molecular co-solvent.
At χIL 0, the activation enthalpy was positive and the activation entropy large and negative, as would be expected for a non-spontaneous reaction that proceeds through an associative mechanism. Upon moving to χIL 0.3, there is a marked decrease in both the activation parameters. A decrease in the activation enthalpy corresponds to a decrease in the activation energy (an enthalpic benefit) whilst a decrease in the activation enthalpy results in an increase in the activation energy (an entropic cost). Increasing the proportion of ionic liquid further does not significantly affect the activation parameters significantly.
These data are consistent with the observed dependencies of the rate constant of the reaction on the proportion of the ionic liquid [bmim][(CF3SO2)2N] in the reaction mixture. At low proportions of [bmim][(CF3SO2)2N], the decrease in the activation enthalpy outweighs the decrease in the activation entropy; the result is a decrease in the activation energy and an increase in the rate constant. Whilst the uncertainties in the activation parameters negate detailed interpretation, it is likely that the balance of enthalpic and entropic effects changes slightly with the latter dominating and the rate constant decreasing at high mole fractions of ionic liquid.
It is worth considering the microscopic origins of the observed changes in the activation parameters. The decrease in both parameters suggests either a stabilisation of the transition state through interaction with the ionic liquid (which also involves some ordering) or decreased organisation about and stabilisation of the starting materials in the presence of the salt. Given the charge development in the transition state, the former is considered more reasonable.
The outcomes of this reaction can be compared to the previously discussed SN2 Menschutkin reactions of a benzyl halide with pyridine.7,44,49 In those cases, rather than an enthalpic benefit causing the rate constant enhancement, an increase in the activation entropy was observed on moving to ionic liquid solvent; there is also corresponding increase in the activation enthalpy. This is consistent with organisation of the solvent about the starting material being more significant in the ionic liquid case; this interaction (which was shown to be between the cation and the lone pair of the nucleophile49 is not present in the transition state. Therefore, while mechanistically similar, the cause of the rate constant increase in the ethanolysis of diethyl phosphate on addition of ionic liquid is not the same.
In comparing the effects of an ionic liquid on the activation parameters of the reaction shown in Scheme 1 with other systems, the SN1 methanolysis of 3-chloro-3,7-dimethyloctane is affected by [bmim][(CF3SO2)2N] in a similar fashion with decreases in both activation parameters observed.42 The origin of the similarities in outcomes is unclear though stabilisation of developing charge in the transition state likely dominates the interaction with the comparatively less interacting nucleophile in the phosphate ester case.
While previously there was an increase in rate constant on increasing proportion of the ionic liquid in the reaction mixture, when lithium bis(trifluoromethylsulfonyl)imide was added, there is a significant decrease in the rate constant of the reaction of ca. thirty-fold at χsalt 0.24 (Table 4). Therefore any positive effects of the parent [bmim][(CF3SO2)2N] are negated. In contrast, when the chloride salt was added, a large increase in the rate constant for ethanolysis of diethyl chlorophosphate was observed, the effect was greater than in the [bmim][(CF3SO2)2N] case (Table 4).
χ salt |
k
2/10−6 mol L−1 s−1![]() |
k
2/k2(ethanol)![]() |
---|---|---|
a Uncertainties quoted are the standard error of the mean derived from at least three replicate measurements at each mole fraction. b Uncertainties are compounded from rate data. | ||
0 | 7.5 ± 0.3 | 1 |
0.24 [Li][(CF3SO2)2N] | 0.26 ± 0.06 | 0.035 ± 0.008 |
0.22 [bmim][Cl] | 83 ± 5 | 11.0 ± 0.8 |
The first thing to note directly is that the activation parameters are markedly different in the presence of the lithium salt. The increase in the activation enthalpy (the large uncertainty in the activation entropy makes further discussion on that point difficult) suggests either increased stabilisation of the starting materials or decreased interaction of the solvent mixture with the incipient charges in the transition state. The former seems unlikely (given the charge development in the transition state) and hence this change is likely the result of coordination of the oxophilic lithium71 with the nucleophile.
The [bmim][Cl] case however is more enlightening. Here, the change in activation parameters is again an enthalpic benefit with an entropic cost much like the [bmim][(CF3SO2)2N] case. The extent of these changes however is markedly different with a marginally lower decrease in enthalpy of activation and a much lower decrease in the entropy of activation. This change in enthalpic benefit/entropic cost balance results in a much greater rate constant enhancement. Of most interest here is the fact that any interactions between the salt and the transition state are reduced in the [bmim][Cl] case when compared to the [bmim][(CF3SO2)2N] case. This may suggest preferential ion pairing of the [bmim] cation with the chloride anion, reducing interactions with the transition state. It is important to note that rate constant enhancement is not proportional to the level of interaction, rather the balance of enthalpic and entropic effects is most important; this has been seen previously for reaction at carbon centres.54 Altering the makeup of both the cation and anion in an effort to alter potential interactions in solution is clearly of great interest.
Footnotes |
† Electronic supplementary information (ESI) available: Representative NMR spectra demonstrating the extent of byproduct formation. Kinetic data on which Fig. 1–4 and Tables 1, 2 and 4 are based. Kinetic data and Eyring plots from which the data in Tables 3 and 5 are derived. Kinetic data, Eyring plot and activation parameters for the ethanolysis of diethyl chlorophosphate in acetonitrile. See DOI: 10.1039/c4nj01224k |
‡ 31P NMR spectroscopy was considered, particularly in order to identify any formation of byproducts. The overlap of signals in such spectra for many of the solvents systems used meant that 31P NMR spectroscopy was subsequently not used for kinetic analysis. For more details, see ESI.† |
This journal is © The Royal Society of Chemistry and the Centre National de la Recherche Scientifique 2015 |