Alexandre
Simula
a,
Fernando
Ruipérez
a,
Nicholas
Ballard
a,
José R.
Leiza
a,
Steven
van Es
ab and
José M.
Asua
*a
aPOLYMAT and Kimika Aplikatua Saila, University of the Basque Country UPV/EHU, Joxe Mari Korta Zentroa, Tolosa Hiribidea 72, 20018, Donostia/San Sebastián, Spain. E-mail: jm.asua@ehu.eus
bDispoltec BV, Urmonderbaan 22, 6167 RD Geleen, The Netherlands
First published on 15th November 2018
The alkoxyamine 3-(((2-cyanopropan-2-yl)oxy)(cyclohexyl)amino)-2,2-dimethyl-3-phenylpropanenitrile (Dispolreg 007) has recently emerged as a robust regulator for the controlled (co)polymerization of methacrylates and styrene by nitroxide mediated polymerization. However, whilst the successful use of this alkoxyamine in a number of systems has now been demonstrated, the fundamental reason behind its ability to overcome the long-standing issues that normally arise in the nitroxide mediated polymerization of methacrylates has not been conclusively determined. Herein, a combination of quantum chemical calculations, mathematical modelling and experimental data is utilized to compare the behaviour of Dispolreg 007 with alternative nitroxides such as SG1 and DPAIO. The chemical structure of Dispolreg 007 leads to a nitroxide which has a higher rate of combination and a lower rate of β-H transfer with the propagating radical compared to SG1. This is not the case for DPAIO which shows relatively similar rate constants. The comparison between modelled and experimental data allowed for a better understanding of the system, with small amounts of macromonomers being produced mainly by β-H transfer at high monomer conversions (>80%). Finally, the estimation of the rate constants was key to optimize the polymerization conditions to achieve good control and high livingness at high monomer conversions (90%).
Extensive work has been carried out on the development of new nitroxides/alkoxyamines to shift the dormant/active equilibrium and to avoid undesired disproportionation events. The use of α-hydrogen nitroxides/alkoxyamines which include 2,2,5-trimethyl-4-phenyl-3-azahexane-3-nitroxide (TIPNO) and N-tert-butyl-N-[1-diethylphosphono-(2,2-dimethylpropyl)] (SG1) nitroxide has improved the ability to polymerize styrene and acrylate derivatives with good control, though the polymerization of methacrylates is still problematic.6–8 More specifically for SG1/BlocBuilder, β-hydrogen transfer could be reduced for one monomeric unit, unfortunately not for the polymer, due to a strong penultimate effect.9 The penultimate effect also affected the rate constant for recombination (kC) and decomposition (kD) of the macroalkoxyamine and as a result did not allow good control and “livingness” according to theoretical calculations, which was confirmed by experimental data.10 This issue could be overcome through copolymerization of methacrylates with small amounts of co-monomers, as introduced by Charleux et al.11,12
The development of more complex structures introduced by Grubbs et al. could regulate the polymerization of methacrylates, but their versatility towards other monomers was limited.13–15 The development of imidazolines and imidazolidines gave rise to various trends as intramolecular and/or intermolecular disproportionation could be favoured depending on the nature of the substituents but could not be avoided.5 2,2-Diphenyl-3-phenylimino-2,3-dihydro-indol-1-yloxyl (DPAIO) and its derivatives yielded the first successful polymerization of methacrylates, obtaining good control and minimal undesired disproportionation.16,17 Nevertheless, the limited versatility of the nitroxide towards the polymerization of styrene and acrylates is still problematic, though some modifications on the alkoxyamine structures improved the control in the case of styrene.18 The development of new nitroxide/alkoxyamine structures has also been motivated to reduce the tendency of undesired disproportionation between two nitroxides (self-disproportionation).19,20 The instability of the nitroxide can be favoured by the presence of hydrogen atoms on the carbon α to the N–O bond, which can affect the control of the polymerization thereafter.15,19,21,22
Recently, our research group has introduced Dispolreg 007 (3-(((2-cyanopropan-2-yl)oxy)(cyclohexyl)amino)-2,2-dimethyl-3-phenyl-propanenitrile), which proved to be a robust and readily obtainable regulator for the polymerization of methacrylates and styrene.23,24 In addition, copolymerization (random, gradient and block) could be achieved in both solution and dispersed media.25–27 Nevertheless, the underlying reasons for the success of Dispolreg 007 in the polymerization of methacrylates remain unclear and little is known about the reactivity of the corresponding nitroxide with itself and the propagating radical. Herein, we seek to gain a deeper insight into the stability of the generated nitroxide, the polymerization of methyl methacrylate (MMA) with Dispolreg 007, and the nature and importance of the β-H transfer reaction between the propagating radical and the nitroxide. A combination of quantum chemical calculations, mathematical modelling and experimental data is utilized to estimate the values of the various rate constants and compare the level of control exerted by Dispolreg 007 with that of other known alkoxyamines/nitroxides, such as BlocBuilder/SG1 and DPAIO.
In order to get a deeper insight into the ability of this apparently stable nitroxide to mediate the polymerization of MMA, we first focused on the structural differences between poly(methyl methacrylate) based macroalkoxyamines, with SG1 (PMMA-SG1), DPAIO (PMMA-DPAIO) and D7 (the nitroxide part of Dispolreg 007 – PMMA-D7) as end groups. The presence of various functional groups along the N–O bond greatly influences the stability of the macro-alkoxyamines, which is reflected in the different values of kD and kC. The values of kD estimated from experimental data have been reported for D7 and SG1, with kPMMA-D7D = 1.4 × 10−2 s−1 (at 90 °C), kPMMA-SG1D = 1 × 10−2 s−1 (at 45 °C). In addition, the value for DPAIO has been estimated from the cleavage temperature to be kPMMA-DPAIOD = 2.02 × 10−4 s−1 (at 120 °C).10,16,24 In order to understand the differences in the values of kD and estimate the value of kC for PMMA-D7, quantum chemical calculations were carried out. For this, high level methods are available and have been utilized to determine the dissociation energies of model alkoxyamines.36 However, in spite of their accuracy for relatively small molecules (∼100 electrons),37 these high level methods are extremely time consuming for larger molecules (PMMA-D7 has 292 electrons; the details are given in the ESI†). Therefore, DFT was used in this work, as it gives a relatively good estimation of bond dissociation energies (BDE) and geometrical parameters.38,39 The molecular model used for the poly(methyl methacrylate) radical contains two MMA units (Fig. S1† and n = 1 in Fig. 1). The geometrical parameters corresponding to the macro-alkoxyamines/nitroxide and the respective Cartesian coordinates are detailed in the ESI (Tables S1 and S2†). The reliability of the level of quantum calculation was tested by comparison with the predictions of higher level calculations (CBS-4 M,40 Table S3†) and in the determination of the geometrical parameters of known nitroxides31,41 and the BDE of known alkoxyamines,39 with minimal error observed (Tables S4–S6†). It should be stressed that while the quantum chemical calculations provide valuable details on the molecular and electronic structures, this can only yield a qualitative insight into the order of magnitude of the rate constants.
The analysis of the C–O bond strength and stability of the different macroalkoxyamines can give an estimation of the order of magnitude of kD for PMMA-D7, PMMA-SG1 and PMMA-DPAIO.39 The results are summarized in Table 1 and show that the bond dissociation energy (BDE) of the C–O bond is slightly higher for PMMA-D7 (142.30 kJ mol−1) compared with that of PMMA-DPAIO (140.61 kJ mol−1) and substantially higher than that for PMMA-SG1 (84.87 kJ mol−1), which might be explained by several features.
BDE | R e(C–O) | η(σCO*) | |
---|---|---|---|
Dispolreg 007 | 142.30 | 1.439 | 0.041 |
DPAIO | 140.61 | 1.447 | 0.043 |
SG1 | 84.87 | 1.454 | 0.047 |
First, on comparing the C–O bond distance for the macroalkoxyamines, it is observed that C–O is 0.015 Å shorter in PMMA-D7 than in PMMA-SG1 (1.439 vs. 1.454 Å), which may be ascribed to a lower charge transfer to the antibonding sigma (σCO*) orbital. This is reflected in a smaller occupation number in that orbital for PMMA-D7 (0.041) compared to that for PMMA-SG1 (0.047), thus strengthening the C–O bond. The slightly lower BDE of PMMA-DPAIO (140.61 kJ mol−1) compared to that of PMMA-D7 (142.30 kJ mol−1) might also be due to a small elongation of the C–O bond (0.008 Å) and a slightly larger population of the antibonding σCO* orbital in PMMA-DPAIO (0.043 compared to 0.041 for PMMA-D7). The higher stability of PMMA-D7 compared to that of PMMA-SG1 suggests that the recombination of the nitroxide with the propagating radical is favored; hence kPMMA-D7C/kPMMA-D7D > kPMMA-SG1C/kPMMA-SG1D. In addition, the extensive stabilization provided by the aromatic rings in DPAIO should yield kPMMA-D7C/kPMMA-D7D relatively similar to kPMMA-DPAIOC/kPMMA-DPAIOD.
There is a discrepancy between the values of BDE obtained herein and the activation energies of the macroalkoxyamines PMMA-D7 (109 kJ mol−1),24 PMMA-SG1 (≈100 kJ mol−1)2 and PMMA-DPAIO (126 kJ mol−1)16 estimated using experimental data. However, comparison is not straightforward because the confidence intervals of the previously reported values were not given and analysis of the procedures used to estimate these activation energies indicate that these intervals may be significant (see the ESI† for details).
The second set of calculations focused on the undesired disproportionation reaction between a propagating radical and the nitroxide. Disproportionation can either be intermolecular or intramolecular (Scheme 1) with the associated rate constants kTNC and kTND respectively. When SG1 and DPAIO are employed, it has been demonstrated that the intramolecular path plays a minor role in the formation of macromonomers.7,17 Hence, it is important to check if intramolecular disproportionation also plays a minor role with PMMA-D7. It should be noted that the intramolecular disproportionation reaction cannot be predicted by either the structure of the nitroxide/propagating radical or by the geometry of the transition state.5 The comparison of the transition state energy of the intra- and intermolecular disproportionation highlighted that for PMMA-D7, the intramolecular path (ΔHTSintra = 149.06 kJ mol−1) is kinetically disfavored compared to the intermolecular one (see below, ΔHTSinter = 85.74 kJ mol−1). The same behavior is observed for SG1 where values of ΔHTSintra = 133.85 kJ mol−1 and ΔHTSinter = 76.56 kJ mol−1 were obtained respectively. Overall, the transition states are lower in energy for SG1 than for D7, which suggests that the unlikely intramolecular disproportionation reaction is kinetically easier for SG1.
The predominance of the intermolecular disproportionation process required further calculations. Hence, the potential energy surfaces (PES) were calculated and they are depicted in Fig. 1. The values of the transition states (Fig. 1 and Table S7†) reveal that the intermolecular disproportionation between PMMA and the nitroxide based on Dispolreg 007 (D7) shows a much higher activation energy (ΔHTS = 85.74 kJ mol−1) compared to DPAIO (77.29 kJ mol−1) and SG1 (76.56 kJ mol−1), which presents the lowest value. Consequently, the β-H transfer is kinetically easier for SG1 and DPAIO compared to D7. Thermodynamically, the reaction is favored with SG1 which shows the most exothermic enthalpy of reaction (ΔHreac = −90.68 kJ mol−1). This process is thermodynamically less likely to occur with D7, which presents the least negative enthalpy of reaction (ΔHreac = −81.98 kJ mol−1), whilst an intermediate value is found for DPAIO (ΔHreac = −86.72 kJ mol−1). Consequently, we should expect the rate constant kPMMA-D7TNC to be lower than the recently reported value for PMMA-SG1 of kPMMA-SG1TNC = 7.6 × 103 M−1 s−1 at 90 °C.42
The smaller barrier observed for the intermolecular disproportionation in PMMA-SG1 may be explained in terms of the geometrical structure of the transition state. In the case of PMMA-SG1, the transition state resembles the reactants (early transition state43), thus facilitating the molecular rearrangement to surpass the reaction barrier. For PMMA-D7, the transition state occurs later, with geometrical features of the products more developed, as shown in Table 2 and further displayed in Fig. 2 and S2.†
![]() | ||
Fig. 2 Molecular models used in the quantum chemical simulations shown during the transition states of the intermolecular disproportionation reaction. The general molecular structure of the transition state in the intermolecular disproportionation reaction can be found in Fig. S1 and S2 of the ESI.† |
R e(C–C) | R e(C–H) | R e(O–H) | α(CHO) | |
---|---|---|---|---|
D7 | 1.412 | 1.282 | 1.380 | 92.7 |
DPAIO | 1.410 | 1.250 | 1.431 | 93.0 |
SG1 | 1.421 | 1.232 | 1.453 | 97.1 |
During the formation of the macromonomer, the C–C bond in the transition state is shorter for D7 than for SG1 (RCC = 1.412 Å vs. 1.421 Å), which means that with D7 the CC bond is more developed at this stage. Similarly, the hydrogen is farther from the carbon of the incipient double bond of the PMMA–C
CH2 macromonomer in the case of D7 (RCH = 1.282 Å) compared to that of SG1 (RCH = 1.232 Å). The developing hydroxylamine also presents interesting features, with a shorter NO–H bond in the case of D7-H (ROH = 1.380 Å) than of SG1-H (ROH = 1.453 Å). DPAIO-H shows geometrical features closer to those of D7-H (see Table 2). This suggests that the transition state for the β-H transfer reaction with DPAIO also occurs later, thus hindering the disproportionation compared to SG1. Nevertheless, the obtained geometrical values are rather similar, which suggests that the electronic distribution might also play a key role in the reactivity of Dispolreg 007.
In summary, three characteristics are likely to influence the control over the polymerization of MMA with Dispolreg 007: (1) the stronger C–ON bond of the macroalkoxyamine PMMA-D7 compared to that of PMMA-SG1, due to a higher BDE. This difference might be the primary cause for control over the polymerization of MMA; (2) the higher energy transition state in the intermolecular disproportionation reaction kinetically disfavors the process compared to that of the PMMA radical and SG1; and (3) thermodynamically, this process is the least exothermic compared with SG1 and DPAIO nitroxides.
Nevertheless, it is important to point out that quantum chemical calculation can only give a hint on the behavior of the system. In order to get a deeper insight into the polymerization and to get a better estimation of the kinetic parameters (e.g. kC and kTNC), we relied on experimental data and mathematical modeling. Hence, polymerizations of MMA at 90 °C in the presence of various concentrations of Dispolreg 007 were carried out. The initial monomer concentration was maintained at [M]0 = 4.35 mol L−1 in toluene, and the initial alkoxyamine concentration [RN]0 was modified to target various chain lengths, DPn = [M]0/[RN]0 of 50, 100 and 200. The polymerizations were carefully monitored by 1H NMR (Fig. S3 and S4†) to follow the monomer conversion and the production of the macromonomer (bearing CCH2 functionalities, δ = 6.04 and 5.38 ppm) by either radical–radical termination by disproportionation, or β-H transfer (intermolecular disproportionation) to the nitroxide. The influence of intramolecular disproportionation of the macroalkoxyamine is neglected herein; hence kTND = 0 s−1. This hypothesis is in keeping with our quantum chemical calculations. The evolution of molar mass as a function of monomer conversion was monitored by SEC.
The evolutions of monomer conversion as a function of time, of molar mass as a function of monomer conversion and of the percentage of macromonomers (%MM, calculated as a ratio of the integrals of H2CC of the macromonomers at δ = 6.04 and 5.38 ppm and the integral of the polymeric CH3; ESI†) as a function of monomer conversion were fit using a mathematical model (coded in MATLAB, Scheme S2†), based on the reaction in Scheme 1. The influence of oxygen inhibition is introduced, which affects the radical concentration (and monomer conversion) at the beginning of polymerization. The concentration of the inhibitor selected herein is [O2]0 = 1.94 × 10−4 mol L−1, which is 10 times lower than the highest [O2]0 that can be dissolved in toluene, as calculated from the molar fraction of oxygen that can be dissolved in toluene at 25 °C.44 The unknown rate constants for alkoxyamine recombination (kCA), macroalkoxyamine recombination (kC), disproportionation with nitroxide (kTNC) and inhibition with oxygen (kinh) are estimated by fitting the experimental data. Other rate constants and initial values are taken from the literature and are summarized in Table S8.† It is important to note that a better fit was obtained when using a chain length dependent rate constant for termination defined as ki, it = k1, 1t× i−αS for i ≤ 100 and αS = 0.63, hence where radicals of chain lengths CL ≤ 100 are involved.45–47 The rate constant k1, 1t corresponds to the monomeric rate constant for termination k1,1t = 2.33 × 1010 × exp(−9000/RT) M−1 s−1.47 For radicals with chain lengths CL > 100, the chain length dependent rate constant of termination is expressed as ki, it = k1, 1t × 100(−αS+αL)× i−αL, with αL = 0.16.47,48 Finally, for radicals of chain lengths CL > 1000, the average rate of termination <kt> was approximated to a singular rate constant of termination kt = 1.984 × 108 × exp(−5890/RT) M−1 s−1.48
Fig. 3 shows that the experimental evolutions of monomer conversion as a function of time, of molar mass as a function of monomer conversion and of percentage of macromonomers as a function of monomer conversion were well fitted by the mathematical model using the following values (at 90 °C) of kCA = 3.1 × 106 M−1 s−1 (±4.1 × 103 M−1 s−1), kC = 3.7 × 105 M−1 s−1 (±2.0 × 101 M−1 s−1), kTNC = 8.3 × 102 M−1 s−1 (±4.0 × 10−2 M−1 s−1) and kinh = 34 × 107 M−1 s−1 (±4.7 × 105 M−1 s−1). The experimental data in Fig. 3A and B are obtained from online NMR measurements, while the experimental data in Fig. 3C are obtained from offline SEC measurements. Offline NMR measurements were also performed (Fig. S5†) and relatively good fitting could be achieved with the mathematical model.
At low reaction times (<20 min), the fitting of the monomer conversion for high DPn is not optimum. This might be ascribed to a strong penultimate effect and to side reactions between radicals and oxygen, which might explain the high value obtained for kinh herein. Both effects do not play a role further on, allowing good fitting between the experimental and modeled data at high monomer conversions (>50%) for all targeted DPns. The importance of using chain length dependent rate constants for termination is further highlighted in Fig. S6,† wherein a singular value is selected for kt. Using the previously utilized value of kt = 6.07 × 107 M−1 s−1,16 poorer fitting is obtained for both the evolution of monomer conversion and of macromonomer production. Fig. 3B shows that the %MM increases when DPn decreases due to a higher nitroxide concentration which is likely to increase the probability of the disproportionation reaction between the nitroxide and a propagating radical. The relative importance of the radical–radical termination by disproportionation vs. nitroxide–radical disproportionation was estimated using the mathematical model. In Fig. S7,† the evolution of the percentage of macromonomers generated by radical–radical termination by disproportionation or nitroxide–radical disproportionation is presented. At DP50, it becomes clear that the disproportionation reaction with a nitroxide is predominant at high conversions, above 80%. In Fig. 3C, the mathematical model could fit the experimental values of molar mass relatively well for all targeted DPns. At low conversions the model predicted relatively high molar masses. The reason is that at the beginning of the process the concentration of the nitroxide was very low and the polymer chains grow to high values.
The estimated value of kCA is in the same order of magnitude as that for an SG1-based alkoxyamine with an initiating tertiary radical (from 24 to 120 °C).49,50 This is not the case for the macroalkoxyamine, where kPMMA-D7C = 3.7 × 105 M−1 s−1 is one order of magnitude higher than kPMMA-SG1C = 1.4 × 104 M−1 s−1 (at 45 °C).10 First, the decrease of an order of magnitude for kC compared to that of the parent alkoxyamine is expected, due to a penultimate effect.29 Then, the difference in the two macroalkoxyamines coincides with the quantum calculations above, suggesting a higher stability of the Dispolreg 007 based macroalkoxyamine. The estimated value of kPMMA-D7C is also lower than kPMMA-DPAIOC = 1.45 × 106 M−1 s−1 (at 40 °C).16
Regarding the β-H transfer reaction, the estimated rate constant for kTNC = 8.3 × 102 M−1 s−1 (at 90 °C) is lower for Dispolreg 007 than for SG1 (kPMMA-SG1TNC = 7.6 × 103 M−1 s−1 at 90 °C).42 This observation is in line with the quantum chemical calculations wherein disproportionation is kinetically less favored in the case of Dispolreg 007. Using this estimated value of kTNC, it is possible to calculate the theoretical fraction of disproportionation, introduced by Souaille and Fischer, as fD = kTNC/(kTNC + kC) = 0.22% (at 90 °C).51 This is indeed far lower than what was calculated for PMMA-SG1 (fD = 7.8% at 70 °C) as a result of a higher kC and a lower kTNC.7 The fraction of disproportionation cannot be accurately estimated for DPAIO due to the presence of two breakable bonds C–O and N–O.17
The maximum monomer conversion Cmax one could achieve with Dispolreg 007 can be calculated utilizing the equations derived by Souaille and Fischer (see the ESI† for details).52 More specifically, for an initial alkoxyamine concentration [RN]0 = 2.17 × 10−2 mol L−1 (corresponding to a targeted DP of 200 using the conditions above), the polymerization of MMA with Dispolreg 007 can reach 98% conversion theoretically, whilst only 90% can be achieved with DPAIO and 69% with SG1.
It is important to stress that the various equations derived by Souaille and Fischer on the features of nitroxide mediated polymerization are obtained assuming chain length independent rate constants.51–53 Herein, the same assumption clearly cannot be made, due to the importance of the chain length dependency of the rate constant for termination kt for a more accurate fitting of our experimental data. Therefore, the calculated values of fD and Cmax are given above as approximate indications and should be treated with caution.
The (kC, kD) Fischer's diagram often used to predict the features of polymerization is based on chain length independent rate constants and in the absence of β-H transfer between the radicals and the nitroxide. Therefore, this type of diagram is not valid herein.51,53 Consequently, we relied on mathematical modeling to get an insight into the predictive “livingness” of the polymerization process at high monomer conversions (90%).
Fig. 4 presents the effect of the value of targeted DPn and the polymerization temperature on the fraction of “living” chains at 90% monomer conversion for the batch solution polymerization of MMA ([M]0 = 4.35 mol L−1). The simulations were carried out with the parameters in Table S4,† assuming no inhibition. Some inaccuracies are expected due to the unavailability of the activation energies of kCA, kC, and kTNC, which forces us to assume a weak influence of temperature on these rate constants.
![]() | ||
Fig. 4 Temperature and targeted DPn conditions for high “livingness” at high monomer conversion (90%). |
Fig. 4 shows that at a reaction temperature of 363 K (90 °C), the percentage of “living” chains at 90% monomer conversion can change dramatically, ranging from 77% to 7% in a range of targeted DPn of 5 to 595. This is in part due to a difference of the time needed to reach 90% conversion, which will affect the percentage of “dead” chains.54 As such, “livingness” can be maximized at lower monomer conversions (and reaction times) during the formation of the PMMA based macro-alkoxyamine thus resulting in good chain extension. This was indeed the strategy utilized in our previous report,24 wherein the polymerization of MMA was stopped at 68% monomer conversion (64% theoretical “livingness”) to yield a well-defined PMMA-b-poly(butyl methacrylate) block copolymer (Mn = 19800 g mol−1, Đ = 1.46).
Under the batch conditions herein, a decrease of 10 °C of the reaction temperature could significantly increase the “livingness” above 90% conversion, from 84% to 23% for DPn ranging from 5 to 595. The influence of reaction time is highlighted in Fig. S8,† wherein the percentage of “living” chains is plotted as a function of temperature and time to reach 90% monomer conversion.
In short, high “livingness” can be achieved for the polymerization of MMA in the presence of Dispolreg 007, with minimum macromonomer formation from β-H transfer with the nitroxide. This can be achieved over a range of temperatures and molar masses. However, targeting high “livingness” at high monomer conversion can come at a cost of long reaction time or narrow range of targeted molar masses. It should be emphasized that for shorter reaction times, optimized conditions for “livingness” and control can also be achieved using a slight excess of the nitroxide at the onset of the polymerization or mild reducing agents.55–57
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c8py00900g |
This journal is © The Royal Society of Chemistry 2019 |