Amjad Alia,
Muhammad Nadeemb,
Jinwei Lua,
Jamile Mohammadi Moradiana,
Tahir Rasheedc,
Tariq Azizc,
Chanez Maouchea,
Yintian Guoc,
Muhammad Awaise,
Fan Zhiqiangd and
Li Quo*a
aResearch School of Polymeric Materials Science & Engineering, Jiangsu University, Zhenjiang, 212013, PR China. E-mail: liguo@ujs.edu.cn
bDepartment of Environmental Engineering, Wuhang University of Technology, Wuhan, 430223, PR China
cInterdisciplinary Research Center for Advanced Materials, King Fahd University of Petroleum and Minerals, Dhahran 31261, Saudi Arabia
dMOE Key Laboratory of Macromolecular Synthesis and Functionalization, Department of Polymer Science and Engineering, Zhejiang University, Hangzhou, 310027, PR China
eResearch Center of Fluid Machinery Engineering and Technology, Jiangsu University, Zhenjiang, 212013, PR China
First published on 27th September 2021
The kinetics and mechanism of ethylene and cyclic diene 5-ethylidene-2-norbornene (ENB) copolymerization catalyzed by rac-Et(Ind)2ZrCl2/[Ph3C][B(C6F5)4]/triisobutylaluminium (TIBA) were investigated using a quench-labeling procedure using 2-thiophenecarbonyl chloride (TPCC). The E/ENB copolymers were characterized by gel permeation chromatography (GPC), differential scanning calorimetry (DSC), scanning electron microscopy (SEM), and 1H nuclear magnetic resonance (NMR) spectroscopy and sulfur analysis. To reduce the errors of the ethylene–diene copolymerization for the kinetics study, we selected E/ENB with steric and electronic features that permit us to elucidate the metallocene catalyst behavior against dienes. A quantitative approach of catalyst speciation, stereodynamics, and micro-kinetics assisted the resolution of mechanistic problems, such as the elastomeric synthesis of ethylene propylene diene monomer rubber (EPDM), the catalyst resting state nature, and how much ion-pairing occurs during polymerization. We report here the precise observation of metal–polymer species, explanation of the dynamics of their initiation, propagation, and termination, and ethylene ENB copolymer development. An approach based on acyl chloride was used to selectively quenched transition metal–polymer bonds to evaluate the polymeric catalyst in terms of its reaction rate, Rp, propagation rate content, kp, and mole fraction of active centers. It is noted that the decline in catalytic activity in the range of 1800 s, and the active center [Zr]/[*C] fraction significantly increased during the initial 1000 s and then tended towards a steady figure of 86%. It is suggested that nearly complete initiation of all olefins catalysts can be obtained after a sufficiently extended reaction. The quick increase in active sites in the first stage can be described by the immediate initiation of active sites positioned on the surfaces of catalyst particles. The initial polymerization rate, Rp, is high and the crystalline properties of the E/ENB copolymer are low due to the greater incorporation of ENB in the polymer backbone, and later the polymerization reaction rates remained stable with a lower mol% of ENB. The melting temperature (Tm) ranges from 108 to 127 °C, whereas the crystalline temperature ranges from 63 to 108 (J g−1). In the E–ENB copolymers, the value of kpE is much greater than that of kpENB; at 120 s, the kpE and kpENB values are 9115 and 431 L mol−1 s−1, respectively, implying smaller diffusion barriers in the early stages, which are close to the actual propagation rate constant.
Since the discovery of olefin polymerization catalysts in the 1950s, they have been a source of fascination for scientists. Catalysts are essential in manufacturing ethylene/viny (ENB,VNB) derivatives (ENB, VNB) copolymers. The majority of catalysts in the literature are group 4 single-sight metallocenes, particularly ansa-zirconocenes, which are suitable polymerization catalysts for ethylene, and higher α-olefins.18 The potential of these catalysts to yield ethylene-based copolymers with homogeneous composition and insert higher olefins into the PE backbone is one of their most attractive features.7 Regardless of their renowned industrial accomplishments over the last 60 years, catalytic olefin polymerizations are still far from being completely understood.19,20 Surprisingly, there are no particularly well-defined rate laws for the basic initiation, propagation, and termination mechanisms for olefin polymerizations with beneficial productivity.20–22 Without understanding the concentration and speciation of catalytically active species, comparing empirical rate laws with molecular mechanisms is complicated. Since their discovery, it has been challenging to determine the concentration of [*C]/[Zr] active centers and the corresponding molecular kinetic constants in polyolefin catalysts.20,23,24 The first investigation, conducted in the late 1950s, involved adding a labeling agent to the end of propagating polyolefin chains. Several methods have been reported in the literature using different compounds, including CO, O2, I2, CS2, H, CH3COCl, and MeOH; however, they have disadvantages. In the case of carbon monooxide (CO), more than five insertions of the ligand exceed the capacity of the active site. Also, electrophilic and nucleophilic labeling compounds such as CS2, MeOH, and the oxidative labeling agents I2 and O2 are not precisely selective toward carbon–metal bonds and consequently also label “dead” polymer chains bound to (alkylaluminum MAO, TIBA, TEA) cocatalysts.21,25–28
In the 2020s, Bochmann and Landis et al. used the quenched-flow (often called “stopped-flow”) polymerization technique to develop a new way of measuring active centers and propagation rate constant (kp).23,28,29,32 The time correlation of polymer molecular weight (MW) and polymer yield have been well defined under controlled conditions, i.e., before chain transfer and termination processes become noticeable (characteristically shorter than 1 s).21,28–30 On the other hand, Landis et al. presented a method based on the employment of chromophore labeling agents.31 The growing polymer labelled chains formulated by MW, using this approach in combination with ultraviolet gel permeation chromatography (UV-GPC), might give insight into metallocene-catalyzed olefin polymerization mechanism theory. Due to a lack of conclusive evidence and the fact that the quench-labeling agent can re-join with the “stabilized alkylaluminum”, it is still unclear if the quenching agent is compatible with alkylaluminum.32,33
Herein, the efficiency of a symmetrical bridge zirconocenes catalyst in incorporating ethylene and propylene homo- and copolymerization is determined using a kinetics model developed in previous studies. By experimenting with different alkylaluminum materials, borates, and metallocene skeleton modification and substitutions, we explored the impact that symmetrical Cp ligands, co-catalyst, and polymerization time have on catalyst activation. To reduce the errors in the ethylene–diene copolymerization for kinetics study, we select E/ENB that has steric and electronic features that permit us to elucidate the metallocene catalyst behavior against the dienes. A quantitative approach of catalyst speciation, stereodynamics, and micro-kinetics could assist the resolution of mechanistic problems, such as the elastomeric synthesis of EPDM, the catalyst resting state nature, and how much ion-pairing occurs during polymerization.34 We report herein the precise observation of metal–polymer species; an explanation of the dynamics of their initiation, propagation, and termination; and ethylene ENB copolymer development and the clear distinction between irregular and continuous behavior. In addition, these results strongly confirm early concepts on metallocene deactivation in solution and why complete initiation of all olefins catalysts can be obtained after a sufficiently extended reaction.
At 50 °C, we performed a series of E–ENB copolymerization reactions using a fixed amount of ENB 0.06 mol L−1 in toluene, pre-contact time, polymerization time, and quenching time. The sulfur (S) content of the E/ENB copolymer was measured to calculate the concentration of the active centers using the equation [S] = [C], and the change in [C*]/[Zr] fraction with reaction time was calculated, with the results shown in Tables S1 and S2.† The rate equation Rp = kp[C*][M], which has been well established for most homo- and heterogeneous catalyzed olefins reactions, was used to determine the time-dependent change in the chain propagation rate constant (kp).40 The olefin polymerization rate (Rp) was calculated from differentiations of the curve of the obtained polymers against reaction time (see the ESI, Fig. S4 and S5†). According to the literature, the starting monomer concentration was taken as [M] for this system, which represents an equilibrium of the ethylene concentration in toluene under normal pressure (1 atm and 50 °C 0.085 mol L−1).40 Table 1 summarizes the copolymerization results conclusions, including the time-dependent [C*]/[Zr] fraction, ENB mol% in the polymer, molecular weight, molecular weight distribution, and thermal characteristics data.
Run | Time (s) | Yield (g) | ENB in polb (mol%) | Activity (107 g poly per Mt per h) | *Cc (%) | MWd (g mol−1) | PDd | Tme (°C) | ΔHme (J g−1) |
---|---|---|---|---|---|---|---|---|---|
a Reaction conditions: catalyst (Et(Ind)2ZrCl2) = 1.25 μmol, borate (Ph3C)B(C6F5)4 = 2.50 μmol, ethylene pressure = 0.1 MPa, ENB = 0.06 mol L−1, solvent toluene = 50 mL, activator TIBA 1000 μmol, and TPCC 2000 μmol.b Determined by 1H NMR.c *C calculated from the quenching reaction of the active center [C*]/[Mt] by TPCC (see the ESI).d Determined by GPC.e Determined by DSC. | |||||||||
1.1 | 120 | 0.47 | 5.10 | 1.13 | 16.43 | 58690 | 3.06 | 108.56 | 108.56 |
1.2 | 240 | 0.68 | 4.87 | 0.82 | 29.91 | 59543 | 3.87 | 109.03 | 65.60 |
1.3 | 360 | 0.87 | 4.61 | 0.70 | 45.38 | 84068 | 4.07 | 117.01 | 78.39 |
1.4 | 480 | 1.09 | 3.97 | 0.65 | 54.68 | 85067 | 3.88 | 111.00 | 63.38 |
1.5 | 840 | 1.44 | 3.43 | 0.49 | 75.42 | 86099 | 3.74 | 126.11 | 87.06 |
1.6 | 1200 | 1.59 | 2.98 | 0.38 | 79.2 | 89300 | 3.82 | 127.24 | 109.9 |
1.7 | 1500 | 1.7 | 2.87 | 0.33 | 82.67 | 97692 | 4.07 | 126.72 | 117.9 |
1.8 | 1800 | 1.76 | 2.83 | 0.28 | 86.06 | 100774 | 4.58 | 127.24 | 103.5 |
The influence that the metallocene/methylaluminoxane (MAO) catalyst system has on cyclic dienes is not fully understood, including the Et(Ind)2ZrCl2/TIBA/borate catalyst system. It is well known that for cyclic and linear dienes, increasing the diene content in the feed decreases the polymerization catalyst activity, as expected. The E/ENB copolymer activity decline in the reaction period of 120 to 1800 s, and active center fraction [Zr]/[*C] significantly increased in the initial 1000 s and subsequently tended towards a steady figure of 86%, which is higher than those of ethylene and propylene homopolymerization, but lower than their copolymerization, see Fig. 2. It is suggested that complete initiation of all olefin catalysts can be obtained after a sufficiently extended reaction. The fraction of active sites increased in two stages. In the range of 120–800 s, the active sites quickly reached 75%, and then little change was observed in the period of 800 s to 1800 s. The quick increase in the active sites in the first stage can be described by the immediate initiation of active sites positioned on the surfaces of catalyst particles. Subsequently, these active sites are easily reachable by the monomer and cocatalyst, they can be activated during the pre-contact process, and their initiation can be accomplished fast. The first steady level in the range of 120–800 s can be attributed to the initial induction-steady stage of active sites. In our previous study, we defined the development of the dormant or deactivated active site by the interaction of TMA with the metallocene to the incomplete activation of the metallocene complex. In addition, chain transfer with TIBA or another alkylaluminum will lead to saturated PE chain ends. The β-H transfer of a propagating chain results in a PE chain with a vinyl termination.
Fig. 2 In ethylene/ENB copolymerization, the catalytic activity decreases, and the active center fraction increases with the reaction time. |
Fig. 3 illustrates the polymerization rate curve obtained using a metallocene/TIBA/borate catalyst system at 50 °C. We estimated the changes in ethylene and ENB consumption rate (RpE and RpENB, respectively) over time using the copolymer composition data and yield shown in Table 1. The following equations were used to calculate the rate constants for ethylene and ENB insertion separately:
RpE = [C*][E] | (a) |
RpENB = [C*][ENB] | (b) |
Fig. 3 The active centers and rate of chain propagation changes with reaction time in ethylene/ENB copolymerization. |
After introducing the borate cocatalyst, the polymerization rate slightly decreased in a reaction time of 1800 s (see Fig. 3). During this polymerization time, a significant change in the solvent viscosity was noted, and the construction of fibrous polymer entangled on the stirrer was also observed (see Fig. 4). Scanning electron microscopy (SEM) was used to examine the morphology using a Hitachi S-4800 system connected to an energy-dispersive X-ray spectroscopy (EDS) system. The growing active sites are evidently encapsulated in precipitated polymer throughout this phase, preventing monomers from diffusing into the active sites, resulting in declining rates.
Fig. 4 In ethylene/ENB copolymerization, the construction of polymer fibres in (a) and (b) at 120 s, (c) and (d) at 480 and (e) and (f) in the 1800 s range. |
In the situation that the initial polymerization rate is high and crystalline properties of the E/ENB copolymer are low due to the greater incorporation of ENB in the polymer backbone, later the polymerization reaction rates stay stable with a lower ENB mol%. As the polymerization time increases, the Rp and ENB mol% of the E–ENB copolymers decreases.
The molecular weight of the E–ENB copolymers was lower in the initial period and significantly increased with reaction time. When compared to the molecular weights of ethylene–propylene homopolymers and their copolymers produced with Mt-II/TIBA and Mt/MMAO under similar conditions, EP copolymers produced with lower molecular weight but PE molecular weights are comparable. This means that the chain transfer reaction of the active centers with TIBA in E/ENB copolymerization is less dominant than in the EP system. The molecular weight distributions (MWDs) of E/ENB copolymers were broad and became broader as the reaction progressed, MWD evidently larger than 2, the theoretical polydispersity index of a true single-site metallocene catalyst system, suggesting the presence of multiple active sites in the catalytic system (see Fig. S1†). Furthermore, 1H NMR spectroscopy was used to determine the integrated comonomer content. Fig. 5 shows the curve of ENB mol% in the copolymer against the active centers in the catalyst system. Surprisingly, the E/ENB copolymer system displays a higher reactivity of ENB than E/P/ENB with the same catalyst system in the 1200 s range.11 In addition, ENB is a non-conjugated cyclic diene and exhibits two different double bonds, one cis-CHCH– and another vinyl bond.11
Fig. 5 In ethylene/ENB copolymerization, the incorporation rate and active centers change with reaction time. |
However, the vinyl bond is nearly nonreactive with most of the Ziegler–Natta catalyst systems. ENB polymerizes as an α-olefin with a bulky alkyl group using an endocyclic double (cis-CHCH–) bond, and has more significant ring strain, leaving the vinyl group unreacted. These results look to be closely related to the comonomer's bulkiness.10,43–45
The above-discussed reaction has been confirmed from the peaks of the vinyl groups in the assembled polymer chain (see the ESI, Fig. S2†). In addition, the mol% of ENB in the E/ENB copolymers is higher at the initial time than later on. The active center also quickly increases over the first 120–480 s, meaning that the cyclic diene ENB quickly activates the active site in the initial time that was inactive for ethylene.
The DSC heating curve of the ethylene–ENB copolymers with time is shown in Fig. S3.† Usually, with an increase in ENB from 5 to 10 mol%, the ENB concentration and polymerization time substantially influence the crystallinity and melting temperature (Tm) because ENB is bulkier than ethylene, and propylene prevents the rotational movement of the polymer chains. The melting temperature (Tm) ranges from 108 to 127 °C, whereas the crystallinity temperature ranges from 63 to 108 (J g−1). These results appear to be closely related to the copolymer composition, implying that composition drift might have an unpredictably negative impact on thermal characteristics. In principle, the existence of ENB-rich segments and chains at the initial time of the reaction could indicate the presence of crystallizable segments with varying levels of ethylene sequences, such as chain segments with higher ENB content that melt at lower temperatures and blocks with a lower comonomer content that melt at higher temperatures.
Table 2 summarizes the copolymerization results, including the time-dependent [C*]/[Zr] fraction, with kinetics data such as RpE, RpENB, kpE, and kpENB. As seen in Fig. 6 and 7, both the kpE and kpENB values moderately decreased with reaction time. The effects are comparable to those observed in E and P homo- and copolymerization systems. In the E–ENB copolymerization system, the value of kpE is much greater than that of kpENB. At 120 s, the kpE and kpENB values are 9115 and 431 L mol−1 s−1, respectively, implying more minor diffusion barriers in the early stages, which are close to the actual propagation rate constant.
Run | Time/s | *Cb/% | RpEc/molcat molENB molcat−1 s−1 | RpENBc/molcat molENB molcat−1 s−1 | kpEd/L mol−1 s−1 | kpENBd/L mol−1 s−1) |
---|---|---|---|---|---|---|
a Reaction conditions: catalyst (Et(Ind)2ZrCl2) = 1.25 μmol, borate (Ph3C)B(C6F5)4 = 2.50 μmol, ethylene pressure = 0.1 MPa, ENB = 0.06 mol L−1, solvent toluene = 50 mL, activator TIBA 1000 μmol, and TPCC 2000 μmol.b *C calculated from the quenching reaction of the active center [C*]/[Mt] by TPCC (see the ESI).c RpE and RpENB in units of mmolpoly molMt−1 s−1.d The propagation constant of ethylene and ENB. | ||||||
1.1 | 120 | 16.43 | 62.9 | 3.21 | 9115 | 431 |
1.2 | 240 | 29.91 | 51.2 | 1.96 | 4076 | 165 |
1.3 | 360 | 45.38 | 41.8 | 1.14 | 2193 | 72 |
1.4 | 480 | 54.68 | 34 | 0.689 | 1480 | 38.6 |
1.5 | 840 | 75.42 | 18.5 | 0.145 | 584 | 6.79 |
1.6 | 1200 | 79.2 | 10 | 0.0228 | 301 | 0.987 |
1.7 | 1500 | 82.67 | 6.2 | 0.0075 | 179 | 0.322 |
1.8 | 1800 | 86.06 | 3.875 | 0.00228 | 107 | 0.0965 |
Fig. 6 Variation of the propagation rate constant (kpE) and active centers [C*]/[Zr] with reaction time in ethylene/ENB copolymerization. |
Fig. 7 Variation of the propagation rate constant (kpENB) and active centers [C*]/[Zr] with reaction time in ethylene/ENB copolymerization. |
The previously reported initial kpE values of the ethylene homopolymerization and E/ENB copolymerization kpE are different from each other (4110 vs. 9115). This is because the active center [C*]/[Zr] ratio of E homopolymerization at the initial stage was significantly lower than that of the E/ENB copolymerization, meaning that the [C*]/[Zr] in these two systems must be quite different.12,37,40 The [C*]/[Zr] in the E/ENB copolymerization system that produces copolymers may be made up of loosely associated ion-pairs with large kp values, whereas [C*]/[Zr] in homopolymerization could be made up of a large number of contact ion-pairs with reduced kp values. To make comparisons among the kp values of the E, E/P, and E/ENB homo- and copolymerization systems that are more obvious, the curves of kp vs. the polymerization time are shown in the same graph in Fig. 8.37,42
Fig. 8 Comparison between the propagation rate constant (kpE) of ethylene homo-, ethylene/propylene co-, and ethylene–ENB copolymerization under the same catalyst and conditions. |
The initial value of kp for E/ENB copolymerization is higher than that for E/P copolymerization because the E/P copolymers are almost amorphous and may be dissolved in toluene, while E/ENB exhibits lower crystallinity, and this difference could be due to a higher diffusion barrier in the E/ENB copolymerization system.
The moderate decline in the kp values with polymerization time could also be attributed to the lower diffusion barrier in the E/ENB polymerization system (see Fig. 9). Although, the effect of the active center ratio on the kp should also be considered. The active center fraction in homopolymerization is lower than in E/ENB copolymerization, particularly in the later stages. The active centers generated later in the copolymerization process appear to have a lower kp value than those formed earlier. The active center fraction of the previously reported E and P homopolymerization and their copolymerization with E/ENB make it more interesting. It is well known that P is bulkier than E. The metallocene catalyst in the form of contact ion-pairs may be incapable of accepting P for further coordination, with the result of it becoming dormant or inactive in homopolymerization. In conclusion, in propylene polymerization, the [C*]/[Zr] ratio is lower than that in ethylene polymerization. Similarly, the active center fraction in E/ENB copolymerization is higher than that in ethylene propylene homopolymerization but lower than the copolymerization, meaning that ENB can also activate the active sites in the metallocene catalyst system that were inactive in the homopolymerization of E. According to Tsai et al., zirconocene complexes generally work as precatalysts that need to be activated by alkylaluminum and a borate cocatalyst. The active catalyst species were produced through the reaction of cocatalyst and precatalyst. Scheme 1 shows the ethylene insertion into the Zr–C bond and insertion of ENB into Zr–PE. Even though ENB is sterically more challenging than propylene, the π complexes disappear or are weaker, and the insertion barriers are more significant than for propylene insertions into Zr–C. The addition of ENB to the system resulted in a substantial decline in inactivity. In addition, the steric impact in ENB insertion is mainly caused by the repulsion between ENB and the precatalyst ligand.
Fig. 9 Comparison between the ethylene and ENB propagation rate constants in ethylene/ENB copolymerization. |
Scheme 1 Schematic diagram showing ethylene insertion into the Zr–C bond and the insertion of ENB into Zr–PE. |
In the proposed mechanistic model (Scheme 2), the catalyst active sites bearing a propagation chain (A) could undergo β-H migration, resulting in the generation of L2Zr–H species (B), which may be converted to species (C) via the insertion of ethylene into the L2Zr–H bond. As suggested by Kissin and co-workers, in their model, they clarified the influence of co-monomer activation in the ethylene with propylene and α-olefin copolymerizations using an industrialized heterogeneous polymeric Ziegler–Natta catalyst system.46,47 However, (C) is active for ethylene and higher α-olefin polymerization; it is also capable of reverting to a dormant state (D) due to robust β-agostic contacts between hydrogen on the methyl and Zr metal. The lower percentage of these inactive or dormant sites in the E/ENB copolymerization is the main reason for a higher fraction of active centers than in the previously reported ethylene homopolymerization. According to Fan et al., in E/P copolymerization, the insertion of P in the hydrogen zirconium bond (Zr–H) can bypass the chemical reaction, leading to an increased active center fraction, which is higher than their homopolymerization.20,47
Scheme 2 Schematic diagram showing a mechanistic model of the copolymerization of ethylene and ENB, how the inactive sites are formed, and how ENB can reactivate them. |
ENB polymerizes as an α-olefin with a bulky alkyl group using an endocyclic double (cis-CHCH–) bond and has more significant ring strain. Endocyclic double bond insertion of ENB in the hydrogen zirconium bond can sidestep the reactions, increasing the [C*]/[Zr] fraction ratio in the system through (E). The above-discussed information has been confirmed from the peaks of the vinyl groups in the assembled polymer chain (Fig. S1†). In addition, the formation of (D) from (C) may need suitable adjustments of zirconium and the alkyl (Zr–iPr) moiety, which require enough space among the cation and anion. It is logical to assume that [L2Zr–iPr] Cl–borate species with associated contact ion-pairs are hardly transferred in the dormant state than the loosely associated ion-pairs. This assumption was supported by a noticeable increase in the ENB incorporation rate during the initial 120–360 s of copolymerization (Table 1). As illustrated in Fig. 4, the incorporation rate of ENB is higher at the initial polymerization time, and the active center fraction reached 45%. Upon a further increase in the polymerization time, this fraction reached 85%. This indicates that the catalytic species activated during the initial 120–360 s can insert more ENB than those later activated catalytic species.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1ra06243c |
This journal is © The Royal Society of Chemistry 2021 |