Mohd J. Kamaruddina,
Nam T. Nguyenb,
Georgios A. Dimitrakisa,
Jaouad El harfiab,
Eleanor R. Binnera,
Sam W. Kingmana,
Edward Lestera,
John P. Robinsona and
Derek J. Irvine*ab
aNational Centre for Industrial Microwave Processing, Faculty of Engineering, Department of Chemical and Environmental Engineering, Energy and Sustainability Research Division, University of Nottingham, Nottingham, NG7 2RD, UK. E-mail: derek.irvine@nottingham.ac.uk
bSchool of Chemistry, University of Nottingham, Nottingham, NG7 2RD, UK
First published on 18th December 2013
This paper demonstrates that direct, “in situ” measurement of a reaction mixture using a coaxial probe technique can be used to accurately follow the progress of a chemical reaction. Such a system was shown to clearly indicate the presence/onset of and define the magnitude of key reaction parameters such as induction periods and end-points over a broad range of temperatures and viscosities. Thus it allowed the reaction to be conducted for the ideal time period, so maximising reactor through-put, energy efficiency and end product quality. Furthermore, by relating these ‘in situ’ measurements to a pre-prepared calibration curve, key experimentally achieved reaction rates could be determined. Finally, these continuously acquired, non-intrusive ‘in situ’ measurements were validated by comparison to conventional and industry accepted off-line measurement techniques.
These difficulties in following the true reaction progress are evident in both laboratory and industrially applied ROP of cyclic monomers. Many of these are obtained from natural or agricultural sources, so represent a real opportunity to replace petrochemically derived polymers with more environmentally acceptable alternatives. ROP is a chain growth polymerisation, and many of the synthetic methods used to successfully achieve this exhibit pseudo living/controlled characteristics (i.e. linear growth of molecular weight with both time and conversion).15 A variety of cyclic monomers have been evaluated for ROP such as cyclic ethers,15b,16 amides,16b,17 lactams,16a and esters, which include lactones,18 lactides,19 glycoglides.15,20 Typically, these are CP systems because their mechanisms contain fast, reversible or temporary deactivation reactions as detailed by Duda et al.15a These influence the kinetics by limiting/removing the overall potential for side reactions or termination processes in the similar way to RAFT or ATRP control in FRP.15b,21
The dielectric properties of a particular molecular entity can be expressed in a number of ways to best rationalise/predict how it will respond to an incident electromagnetic (EM) field.22 The real part of its complex permittivity (dielectric constant, ε′) defines the extent to which it will store energy through polarisation.23 Meanwhile, the imaginary part of complex permittivity (dielectric loss factor, ε′′) expresses its ability to dissipate this stored energy into heat.23 A third parameter loss tangent (tan
δ = ε′′/ε′) can then be defined, which indicates the potential for the material to heat under the influence of an applied electric field. Clearly, if an alternating electric field of high frequency is applied to a chemical reaction mixture, then the overall measured dielectric response will be a function of the respective dielectric properties of all the molecular species present and the interactions that occur between them.24 Thus, by definition it should be possible to follow the progress of a specific chemical reaction by monitoring the dielectric changes in the medium if they are altered as a result of the target chemical reaction. There have been only been a small number of papers in which the progress of polymerisation reactions have been followed using dielectric properties across a significantly broad conversion range. Most of these studies have all involved the use of either (a) an exclusive off-line measurement of dielectric properties25 and/or (b) specialised detectors which are inappropriate for use with polymerisations because they have μm gaps which can easily be blocked as the viscosity of the reaction mixture increases.26 In general, direct measurement techniques have found limited application in industrial polymerisation processes due to the inherent difficulty in making measurements in systems that are typically either heterogeneous or exhibit large changes in viscosity. Some use of near infra-red techniques has been made, but many processes still rely on a strategy of sampling followed by an off-line assessment such as acid/hydroxyl value, refractive index measurement, GPC/HPLC data, etc. Furthermore, as these reported dielectric monitoring strategies typically have a requirement to introduce a sampling system and/or specific analysis loop or mesh detector, they are not suited to the constant, direct monitoring regimes required to follow a controlled polymerisation in real time. This can also be said to be true of AC impedance measurements, which would need to be made at much lower frequency to prevent radiation of EM energy. However, even at these lower frequencies, this measurement would involve a capacitor containing a small inter-plate gap being placed into the polymerisation mixture, which would result in a high potential for blockage. Typically, the small number of dielectric property studies that have involved ROP of lactones: have only focused on the monomer without considering the other precursors; have not dealt with the concentration effects of the actual reactant mixture,26b,c and have been collected over a limited temperature range that often did not include the actual reaction temperature used. Thus, because the dielectric property data has been shown to exhibit significant, non-linear changes with temperature, this data is of little real use for trying to monitor the reaction.22 A recent study by Nakamura et al.27 did investigate the change in dielectric properties during the polycondensation of lactic acid by following the dielectric properties of the individual components as a function of temperature. However, this study was limited to the individual components and it did not consider how the dielectric properties of the reaction mixture behaved with reaction time or conversion.
The aim of this study is to rigorously quantify the dielectric properties of all precursors, reactant mixtures and product polymers involved in the tin catalysed ROP of ε-caprolactone (CL) across a broad temperature range (ambient to 180 °C) at a frequency of 2.45 GHz.18h,28 This data is used to follow the progress of this ROP using a coaxial probe.
:
[Sn(Oct)2] ratio was kept constant at 1
:
0.012. For example to achieve a DP of 20, Sn(Oct)2 (17.50 mL of a 0.06 M Sn(Oct)2/toluene solution, 1.05 mmol) and BzOH (9.07 mL, 87.57 mmol) were added via syringe and this solution was mechanically stirred until the mixture was homogeneous, and then heated to the target reaction temperature (150 °C) by being immersed in a pre-heated oil bath. The mixture was rigorously stirred throughout the reaction and small samples were removed via syringe at set times and immediately quenched by immersing the sample vial in liquid nitrogen to halt the reaction so that it remained representative of the reaction system at the time it was taken. Additionally, at this point a direct, ‘in situ’ coaxial probe measurement of the systems dielectric properties was collected. These measurements were compared to off-line analysis of the sampled material which included cavity perturbation (below), gel permeation chromatography (GPC) and nuclear magnetic resonance (NMR).
At the outset of each measurement the properties of the empty quartz sample holders were measured. The samples then were placed into the holders and introduced into the cavity with the use of the linear stage and during the measurements they were repeatedly moved from the furnace to the cavity and vice versa to enable measurements at high temperatures. The different values of frequency shift and quality factor between the empty and the loaded quartz sample holders were recorded, and the complex permittivity of the samples was calculated with the use of the simple perturbation equations. Further details of the experimental technique can be found in the literature.29 The test compound or mixture was then introduced into the tube using a scaled syringe (10–500 mL). Typically, the sample was prepared such that the ratio between its diameter and height is ∼0.1. It was then positioned in the furnace at ambient temperature at the outset of the experiment. The dielectric properties of the sample material introduced into the tube were studied in two ways. (i) Single set temperature measurement: The prepared sample was heated directly to the reaction temperature chosen for the subsequent ROP (150 °C). Once the sample reached the target temperature, a period of 10 minutes was allowed for the material of equilibrate. It was then introduced into the resonant cavity and measurements of the frequency shift and alteration in the quality factor in the cavity upon the introduction of the sample were conducted at 2470 and 912 MHz respectively. (ii) Temperature range measurement: the temperature of the furnace was increased at a rate of 1 °C per minute. When a pre-set measurement temperature was reached, the furnace temperature was kept constant for 10 minutes to stabilise sample temperature before the dielectric measurements were taken. The dielectric properties of the samples were assessed across a range from ambient −180 °C at 10 °C intervals. In both cases, each data point represents the average of 5 measurements, and the standard deviation of these was used to estimate error. Any volume change in the sample during the temperature sweep was quantified and the data corrected. The method and equations used to calculate the permittivity, dielectric loss and dissipation factor (tan
δ) are those detailed in previous publications.22
400 Da calibrated with poly(styrene) narrow standards. All equipment and standards were supplied by Polymer Laboratories (Varian Inc). GPC data was analysed using the Cirrus GPC Off-line software package.![]() | ||
| Fig. 2 Variation in ε′ and ε′′ of CL at 2.45 GHz over temperature range 20–180 °C via coaxial probe technique. | ||
![]() | ||
| Fig. 3 Variation of ε′ and ε′′ of PCL DP56 at 2.45 GHz over temperature range 20–180 °C via cavity perturbation technique. | ||
The dielectric properties of both materials were found to depend strongly upon and non-linearly with temperature. The monomer is a liquid at ambient temperature and, as expected from previous work, its dielectric loss (ε′′) was found to decrease across the temperature range measured.6,7 With the polymer, ε′′ was observed to be very low up to 35 °C. However, at approximately 40–60 °C, it underwent a sharp increase to reach a maximum value at about 55 °C, from which it too gradually decreased until 180 °C. Meanwhile, the behaviour of ε′ for the monomer and polymer were both observed to increase steadily from ambient to 60 °C. From 60 °C to 180 °C, the ε′ value of the polymer remained relatively constant, whereas that of the monomer decreased very slightly. These trends have been attributed to the influence of a number of competing physical and chemical effects.6,7 Below 35 °C, the polymer has a definite solid physical form and thus its dielectric relaxation characteristics are influenced by the presence of intermolecular/inter-chain interactions. Between 45 and 60 °C the material is progressing through its melting point transition (PCL m.p. = ∼55 °C and Tg = ∼−60 °C).6,7 This results in a change in the amount/strength of the intermolecular interactions, a process that reflects a change (increase) in the interaction of the material's electronic structure with the applied electric field. Finally, between 60 °C and 180 °C, the behaviour of ε′′ has been related to the rise in sample temperature resulting in a reduction in sample viscosity/density and thus intermolecular interaction within the sample. Meanwhile, the differing trends in ε′ exhibited by the liquid and solid phases are proposed to be related to the different relaxation characteristics of the materials, where the respective relaxation times are temperature dependent, and a study of this effect will be the subject of a subsequent publication. The key observation that is central to this study is that the measured dielectric responses for the monomer and its related polymer are an order of magnitude different (i.e. a ε′ value of 25 versus 2 respectively). This difference is thought to be partially due to the reduction in the degrees of freedom that each bond/functional group within the polymer has to rotate as the chain is built, which in turn affects the overall dipolar moment of the molecule. A more detailed discussion of the molecular effects is not relevant to this publication, which is concerned with exploiting this difference, but will be included in the future publication mentioned above.
The dielectric responses of the other precursors (ESI Fig. 2†) and two polymerisation mixtures containing different Mw polymer (Fig. 4 and 5) were then evaluated via the cavity perturbation technique. This shows that (a) the dielectric properties did indeed predominantly follow the properties of the monomer and (b) the data was essentially independent of the Mw of the polymer, as the data following a PCL DP 20–monomer mixture exhibits the same trend to that of a PCL DP 200–monomer combination.
![]() | ||
| Fig. 4 Variation in ε′ of monomer and ROP polymerisation mixture at 2.47 GHz vs. temperature via coaxial probe technique. | ||
![]() | ||
| Fig. 5 Variation in ε′′ of monomer and ROP polymerisation mixture at 2.47 GHz vs. temperature via coaxial probe technique. | ||
As shown in a previous paper by the authors,39 the major conclusion drawn from this data is that the key component that will dominate the measured dielectric properties of this specific ROP polymerisation mixture at the chosen reaction temperature will be the monomer.
| PCL (DP) | Molar ratios CL : BzOH : Sn(Oct)2 |
Theoretical Mn at >90% conv. (g mol−1) | Observed Mn at >90% conv. (g mol−1) | Ð | Time to detect polymer (min) | No. of samples >90% conv. |
|---|---|---|---|---|---|---|
| 20 | 20 : 1 : 0.012 |
2054 | 2071 | 1.48 | <30 | 3 |
| 87 | 87 : 1 : 0.012 |
8937 | 8077 | 1.43 | 60 | 2 |
| 200 | 200 : 1 : 0.012 |
20545 | 18751 | 1.62 | 120 | 1 |
Fig. 6 plots a comparison of both ε′ change and Mn build of the polymer (measure off-line by GPC) against reaction time. This clearly shows that the dielectric measurement trends are consistent with loss of monomer as it is converted into polymer.
![]() | ||
| Fig. 6 Comparison plot of ε′ and measured Mn versus time for the same ROP via coaxial probe technique. | ||
:
[CL] ratio and to show that this measurement can be successfully used to follow the course of such a controlled ROP reaction, a calibration curve was constructed using the method detailed in the literature.39 To further evaluate the influence of the measurement temperature on the ε′ value, another identical series of pre-prepared molar mixtures of precipitated PCL (DP = 20) and Cl were prepared via the procedure described above and the dielectric properties of the mixtures were measured using the coaxial probe technique at 80, 100, 120 and 150 °C. Fig. 7 and 8 show the plots of the ε′ and ε′′ values versus molar ratio for these mixtures at different temperatures.
![]() | ||
| Fig. 7 Plot of ε′ vs. molar ratio for pre-prepared PCL/CL mixtures at different temperatures via coaxial probe technique. | ||
![]() | ||
| Fig. 8 Plot of ε′′ vs. molar ratio for pre-prepared PCL/CL mixtures at different temperatures via coaxial probe technique. | ||
Fig. 7 and 8 show that both the recorded ε′ and ε′′ data decrease linearly with increasing molar ratio at a given temperature. Thus, this experiment demonstrates that these plots can be used directly as calibration curves for any PCL ROP polymerisation at a single specific reaction temperature as the data has been shown to be Mn independent. As expected from the survey of the monomer's dielectric properties (Fig. 2), the ε′′ data was observed to be more temperature dependant than ε′. In fact Fig. 8 shows that ε′′ is so dependent on temperature that a new calibration curve would have to be drawn for every change in reaction temperature. However, at a specific temperature the ε′′ relationship is clearly noted to be linear across the whole temperature range, suggesting that it may in fact be possible to use this relationship to indicate the actual reaction temperature achieved for a specific reaction using a specific reactor geometry/conditions. Meanwhile, the data in Fig. 7 shows that the ε′ data is far less dependent on temperature. Furthermore, it demonstrates that above a molar ratio of 47
:
53,which is equivalent to ∼50% conversion, the ε′ values can be regarded as temperature independent across all the temperatures investigated in this study (provided that the polymer is above its melting point). This is because the variance in the gradient of the calibration curves is less that the measurement errors of the technique. The data in Fig. 2 and 3 show that whilst the ε′ value for the polymer is essentially constant from ∼75 to 180 °C, that of the monomer exhibits a significant decrease as the temperature increases from 20 °C to 180 °C. Therefore, this difference in the monomer response to the applied field results in the calibration curves for the lower temperature mixtures presenting a slightly different gradient, which in turn produces a spreading of the curves at low conversion, i.e. high monomer concentration. Closer inspection of the data reveals that at temperatures above 100 °C the gradient difference in the calibration curves had reduced and this meant that the temperature independent window of the ε′ data had been extended from the 50% conversion and above for the 80–100 °C range to approximately 25% conversion and above. This is because the monomer's ε′ reduction is at its most non-linear between 80 and 100 °C, as it is in the transition from cresting the maximum to beginning its linear decrease.
Thus, it was concluded for mixtures that are above the polymer's melting point that; (a) it has been confirmed that ε′ rather than ε′′ is the better property to use for ‘in situ’ monitoring of a broad range of temperatures on the grounds of temperature independence, (b) the size of the operational ε′ window can be related back to the original dielectric survey of the pre-cursors and (c) ε′′ also defines a linear relationship over the whole conversion range of an individual set temperature thus is as conversion sensitive as ε′, but that its greater sensitivity may make it possible to relate the ε′′ data to the actual reaction temperature achieved in a particular reaction vessel/chemical process. However, it should be noted that the measurement errors are greater in the determination of ε′′ because of the lower magnitude of the overall response. Furthermore, in this specific case where the typical reaction temperature used for a ROP reaction of CL is between 120 and 180 °C, the ‘in situ’ data from the polymerisation experiments can be regarded as being temperature independent from 25 to 100% conversion. Therefore, for reaction systems that were operating within this conversion range, it would be necessary to only develop a single correlation between the monomer conversion and ε′ value. Consequently, work is now underway to ascertain if these calibration curves and temperature window estimations can be predicted directly by predictive mathematical calculation from the original dielectric property survey of the reaction precursors, preventing the need to experimentally generate the calibration curve.
![]() | ||
| Fig. 9 Plot of ε′ vs. reaction time for ROP to synthesize PCL of three different DPs at 150 °C via coaxial probe technique. | ||
However, the gradient of this linear relationship increases with the pre-catalyst concentration and could, therefore, be used to rank the polymerisation rates of these comparative reactions, defining the empirically achieved rate of reaction in the polymerisation.
In addition, information on three other features can be extracted from this data.
![]() | ||
| Fig. 10 Plot of ln(1/1 − c) against time for the DP200 and DP87 ROP polymerisations using different Sn catalyst levels. | ||
Fig. 10 defined that the kapp predicted from the coaxial dielectric monitoring of the polymerisation to synthesise the DP 87 polymer is 0.0089 min−1, which is over twice that of the DP 200 reaction (0.0042 min−1). The predicted propagation constant (kp) can then be calculated by using the relationship kp = kapp/[I]o, where [I]o is the initial concentration of initiator (DP 87 [BzOH]o = 0.103 mol L−1 and DP 200 [BzOH]o = 0.043 mol L−1). This gives a propagation rate for DP 87 of 0.086 mol−1 L min−1 and for DP 200 of 0.098 mol−1 L min−1. This data was compared and contrasted with results from an identical DP 87 polymerisation which had its conversion assessed in the conventional manner, i.e. by precipitating the polymer and taking the gravimetric yield. By this method, the predicted kapp and kp were 0.009 mol−1 and 0.092 mol−1 L min−1 respectively. The fact that the actual initiator corrected propagation rates of the dielectric monitoring (both DP 87 and DP 200) and the gravimetric monitoring are essentially the same supports the conclusion that this dielectric monitoring is an accurate method to determine the rate of propagation for a polymerisation.
It had been hoped to use this method of kapp assessment to determine the activation energy (Ea) of the process by plotting a graph of kapp versus 1/Temperature for a series of experiments run at different reaction temperatures. In this case the gradient of the resultant linear relationship is equivalent to Ea/R (where R is the gas constant). However, due to the fact that the ε′ and ε′′ values exhibited some degree of temperature dependence which may influence this calculation, it was concluded that fundamental reasons behind this temperature dependence would need to be further investigated before this analysis was conducted.
| Time of sample (min) | ε′a | Conv.b (%) | Mnc (g mol−1) | Ð | Corrected theoretical Mn (g mol−1) | Conversion observed Mn (g mol−1) | |
|---|---|---|---|---|---|---|---|
a Measured using online coaxial probe at 2.45 GHz.b Determined using 1H NMR (300 MHz, CDCl3).c Determined using GPC calibrated with polystyrene standards (580–370 000 g mol−1). |
|||||||
| 1 | 0 | 18.59 | 0 | 200 | 1.1 | 180 | 200 |
| 2 | 30 | 6.74 | 82 | 1720 | 1.27 | 1880 | 1960 |
| 3 | 45 | 4.12 | 97 | 1880 | 1.47 | 2200 | 2210 |
| 4 | 60 | 3.66 | 98 | 1900 | 1.58 | 2240 | 2250 |
| 5 | 75 | 3.29 | 100 | 1940 | 1.61 | 2280 | 2260 |
The end point of these reactions (i.e. targeted as 90% conversion in this case) was defined as an ε′ value of ∼7 by experimentally measuring the ε′ value of a pre-prepared solution of 10% monomer–90% polymer. A series of polymerisations targeting a product of DP 20 were then conducted, which were terminated as the dielectric response reached a value of 7 to 8, which should approximate to just under 90% conversion. For comparison, the data from a typical DP 20 polymerisation which has its end-point above 90% has been included in Table 1 as entry 1. This data was found to map well on to the experimental data from following the polymerisation in real time (Fig. 9), where the number of samples that were: (a) below a ε′ value of 7, (b) defined by NMR to have greater than 90% conversion (Table 1) and (c) noted to exhibit a deviation away from the linear propagation relationship via GPC (Fig. 9) were found to be identical. This confirms that that the dielectric data can be accurately used to achieve a particular desired molecular structure of product.
The most appropriate dielectric property to use for this ‘in situ’ monitoring process was concluded to be the dielectric constant, ε′, because it exhibits; (a) the largest scale measurement change thus minimising error and (b) is applicable across the broadest temperature range. For these polymerisations, which are conducted above 100 °C due to the physical characteristics of the reagents, the measurements were shown to be temperature independent above 25% conversion to polymer. Thus the generation of a single calibration curve will be sufficient for all reaction conditions applied.
Additionally, by using this monitoring process it has been found possible to calculate the true experimentally achieved reaction rates, i.e. kapp and kp, hence, application of this technique significantly reduces the barriers to the industrial operation of these controlled catalytic polymerisations mechanisms at scale.
Furthermore, this combination of capabilities also has the potential to allow successful product range extension/development of novel differentiated materials, without the need to generate novel chemistry. For example, by indicating the exact point that a second monomer should be added, block copolymers can be accurately manufactured at scale. Hence, the coaxial probe technique exhibits significant potential importance in the field of improving the efficiency of not only polymer but general chemical manufacture. It should have broad applicability to any chemical synthesis/manufacturing processes where detectable changes in dielectric properties result from conducting the desired molecular transformation.
Footnote |
| † Electronic supplementary information (ESI) available. See DOI: 10.1039/c3ra46941g |
| This journal is © The Royal Society of Chemistry 2014 |