Open Access Article
Yifei
He
a and
Christine K.
Luscombe
*b
aDepartment of Materials Science and Engineering, University of Washington, Seattle, USA
bPi-Conjugated Polymers Unit, Okinawa Institute of Science and Technology, Okinawa, Japan. E-mail: christine.luscombe@oist.jp
First published on 10th June 2024
Polythiophenes are one of the most widely studied conjugated polymers. With the discovery of the chain mechanism of Kumada catalyst-transfer polymerisation (KCTP), various polythiophene copolymer structures, such as random, block, and gradient copolymers, have been synthesized via batch or semi-batch (sequential addition) methods. However, the lack of quantitative kinetic data for thiophene monomers brings challenges to experimental design and structure prediction when synthesizing the copolymers. In this study, the reactivity ratios and the polymerisation rate constants of 3-hexylthiophene with 4 thiophene comonomers in KCTP are measured by adapting the Mayo–Lewis equation and the first-order kinetic behaviour of chain polymerisation. The obtained kinetic information highlights the impact of the monomer structure on the reactivity in the copolymerisations. The kinetic data are used to predict the copolymer structure of equimolar batch copolymerisations of the 4 thiophene derivatives with 3-hexylthiophene, with the experimental data agreeing well with the predictions. 3-Dodecylthiophene and 3-(6-bromo)hexylthiophene, which have higher structural similarity to 3-hexylthiophene, show nearly equivalent reactivity to 3-hexylthiophene and give random copolymers in the batch copolymerisation. 3-(2-Ethylhexyl)thiophene with a branched side chain is less reactive compared to 3-hexylthiophene and failed to homopolymerize at room temperature, but produced gradient copolymers with 3-hexylthiophene. Finally, the bulkiest 3-(4-octylphenyl)thiophene, despite its ability to homopolymerize, failed to maintain chain polymerisation in the copolymerisation with 3-hexylthiophene, possibly due to the large steric hindrance caused by the phenyl ring directly attached to the thiophene center. This study highlights the importance of monomer structures in copolymerisations and the need for accurate kinetic data.
Poly(3-hexylthiophene) (P3HT) and its derivatives are the most extensively studied model systems synthesized by KCTP. Various polythiophene copolymers have been developed to reveal the structure–property relationship in CPs. For example, random,15,16 block,15–19 and gradient polythiophene copolymers15,16,20,21 have shown different solid-state morphologies. Compared to that of the random copolymers of the same composition, the morphology of the block copolymers shows further microphase separation, in which the size of the crystalline and the amorphous domains depends on the feed ratio of the comonomers.16,17 In between the random and block distribution patterns, gradient copolymers tend to have more ordered and separated domains than the random copolymers, but with different properties than the pure block copolymers.16 In the design of copolymers with more than two components, the sequence of each segment also plays a role in controlling the morphology.18,19 An important lesson revealed by the preceding studies in copolymer design is the necessity of more precise control over the monomer sequence in tuning the morphology and thus the properties.
Understanding the effect of the monomer structure on the kinetics of polymerisation can lay the foundation for the precise placement of comonomers in KCTP. There is currently no systematic and quantitative study about the kinetic behaviour of thiophene monomers. The lack of kinetic information leads to laborious tuning of experimental parameters for copolymerisation and a misunderstanding of the structure of the final product. A contradiction between the measured monomer reactivity and the predicted copolymer structure due to insufficient kinetic information is observed in the literature. For example, Schmode et al. have qualitatively measured that the reactivity of 2,5-dibromo-3-(6-bromohexyl)thiophene is higher than 2,5-dibromo-3-hexylthiophene, and thus the batch copolymerisation method was adopted to obtain a gradient copolymer.22,23 However, Palermo et al. obtained similar reactivity of 2-bromo-3-(6-bromohexyl)-5-iodothiophene (3BrHT) and 2-bromo-3-hexyl-5-iodothiophene (3HT) and obtained a random copolymer instead of a gradient copolymer in the batch copolymerization.20
With the above inspiration, this study examines the reactivity of the 3HT monomer and a series of thiophene comonomers (Fig. 1(a)–(d)) having different degrees of structural similarity with 3HT in KCTP for one particular concentration. The chosen comonomers vary by the side chains at their 3 positions and have been shown to be compatible with KCTP in either homopolymerisation or copolymerisation.17,24–26
In addition to the reactivity ratios commonly applied to describe the relative kinetics of comonomers,21,27,28 we further calculated the four rate constants, k11, k12, k22, and k21 that define the ratios. The more specific kinetic information will help determine the kinetic difference related to the order of monomer addition, which is helpful in the synthesis of more complex copolymer architectures such as triblock copolymers. We also characterize the monomer consumption rate and the copolymer structure in the equimolar batch copolymerisation to verify the obtained reactivity values.
k appp of each monomer can thus be calculated from the slope of the linear regime of the semi-logarithmic conversion plot in Fig. 2 from triplicated polymerisation of each monomer. The NMR characterization and estimated Mn of each homopolymer are summarized in Fig. S2–S5.† We see that increasing side chain steric hindrance lowers the homopolymerisation rate constants. The three thiophene monomers bearing linear side chains, 3HT, 3DDT and 3BrHT, possess similar homopolymerisation rate constants, while the reactivity of 3OPT homopolymerisation is significantly lower. 3EHT failed to homopolymerize at room temperature using our experimental conditions. The synthesis of P3EHT homopolymers or blocks in copolymers via KCTP using the same Grignard and monomer has been commonly performed at temperatures above 35 °C.17,30,31 The required temperature elevation is due to the increased steric bulk of 3EHT that impedes the coupling of incoming monomers and chain ends.32
![]() | ||
| Fig. 2 The monomer consumption behaviour in homopolymerisation. The homopolymerisation of 3EHT failed at room temperature. | ||
To ensure the strict control of all other factors that may also impact the reaction kinetics, we approximate the reactivity of P3EHT homopolymerisation as infinitesimal in the later quantitative analysis rather than perform the synthesis at an elevated temperature.
Mayo et al. further correlated the four governing parameters in the forms of reactivity ratios as shown below.35 The experiment design and the calculation of the reactivity ratio using Mayo–Lewis equation were based on previous study.21 [M1] and [M2] represent the initial feed of the comonomers measured by t0 aliquots. d[M1] and d[M2] were approximated by the monomer consumption between 1 min and 2.5 min after the initiation of copolymerisation. A series of 6 initial feed ratios covering 2 orders of magnitude of concentration were prepared for each comonomer combination. Each series was then triplicated to generate 15–18 valid datasets for the non-linear fitting. k11 and k22 were approximated by the rate constants measured from the homopolymerisation kinetic studies of each monomer as described in section 3.1.36 Thus, the heteropolymerisation constants k12 and k21 were calculated based on the corresponding reactivity ratio and homopolymerisation constants. The fitted reactivity ratios and calculated rate constants of each pair of comonomer are summarized in Table 2.
| M 1 | M 2 | r 1 | M 2 | R 2 | k 11 [min−1] | k 22 [min−1] | k 12 [min−1] | k 21 [min−1] |
|---|---|---|---|---|---|---|---|---|
| 3HT | 3DDT | 0.89 ± 0.07 | 0.73 ± 0.23 | 0.952 | 0.023 ± 0.002 | 0.018 ± 0.007 | 0.026 ± 0.003 | 0.025 ± 0.01 |
| 3HT | 3EHT | 2.15 ± 0.22 | 0.42 ± 0.23 | 0.939 | 0.023 ± 0.002 | Infinitesimal | 0.011 ± 0.001 | Infinitesimal |
| 3HT | 3OPT | 1.47 ± 0.50 | 1.47 × 10–18 ± 0.46 | 0.217 | 0.023 ± 0.002 | 0.010 ± 0.002 | 0.016 ± 0.006 | 6.80 × 1015 ± 5 × 1015 |
| 3HT | 3BrHT | 1.08 ± 0.32 | 1.24 ± 0.84 | 0.783 | 0.023 ± 0.002 | 0.027 ± 0.003 | 0.022 ± 0.007 | 0.022 ± 0.01 |
The trend of reactivity ratios also correlates with the dissimilarity of the chemical structures of comonomers with respect to 3HT. 3DDT and 3BrHT both have reactivity ratios close to 1 with respect to 3HT, which suggests that these two monomers are kinetically equivalent in the copolymerisation with 3HT. The result is consistent with their chemical structures, where both 3DDT and 3BrHT possess simple and linear side chains. Based on this, the final products of 3HT with either of these two monomers are expected to be random copolymers. In the batch copolymerisation of 3EHT and 3HT, r1 is greater than 1 and r2 is smaller than 1, which suggests that the 3HT monomer is preferably consumed compared to 3EHT. In addition, the product of the reactivity ratios is close to 1, which predicts that the final copolymer will have gradient monomer distribution along the chain. Moreover, by setting kapp22 as infinitesimal based on the failure of homopolymerisation of 3EHT, kapp21 becomes infinitesimal as well. The four rate constants suggest the necessity to build up a 3HT block first before the 3EHT block in the semi-batch copolymerisation of block copolymers. The reactivity ratio of the 3HT and 3OPT monomer pair also indicates the faster consumption of the 3HT monomer in batch copolymerisation, as r1 is greater than 1 and r2 is close to zero. This pair of reactivity ratios predicts that the final product of the batch copolymerisation will have primary 3HT blocks along the chain. Additionally, in the data fitting of 3HT and 3OPT reactivity ratio, the accuracy of the model decreases as shown by the R2 value. The decrease in the accuracy is likely the result of direct attachment of the phenyl ring to the thiophene, which leads to the loss of the catalyst transfer feature under the current experimental conditions thus preventing a chain polymerisation.
:
3DDT predicts that in batch copolymerisation, the final polymer will be a random copolymer with the composition ratio equal to the monomer feed ratio. Given the high similarity of the chemical structures of the two comonomers, the copolymer composition is hard to distinguish solely by NMR signals. Thus, we combined monomer consumption behaviour via GC-MS and polymer weight distribution via MALDI-TOF to characterize the copolymer composition as presented in Fig. 3. 3HT and 3DDT were consumed at a similar rate in the copolymerisation as shown in Fig. 3(a). Fig. 3(b) and (c) presents MALDI-TOF spectra of aliquots from the growing and plateau stages. Both aliquots contain mass peaks that could only represent copolymers as labeled on the spectra. In addition, the ratios of calculated units of 3HT and 3DDT from the signature peaks are close to 1
:
1. This result indicates that the 3DDT composition is present in the copolymer throughout the polymerisation and holds a composition around 50%. The resulting product is a copolymer with random comonomer distribution, which is aligned with the prediction from the reactivity ratio.
![]() | ||
| Fig. 3 (a) Respective monomer consumption behaviour of 3HT and 3DDT, (b) MALDI-TOF spectrum of P(3HT-co-3DDT) at t = 4 min, and (c) MALDI-TOF spectrum of P(3HT-co-3DDT) at t = 60 min. | ||
The reactivity ratio of 3HT
:
3EHT predicts that the batch copolymerisation of this combination will produce copolymers with the monomer composition gradually shifting from 3HT to 3EHT. The prediction was supported by the monomer consumption rate and NMR characterisation over a series of polymer aliquots taken throughout the reaction. 3EHT was less reactive than 3HT, as indicated by the slower monomer consumption rate shown in Fig. 4(a). As shown in Fig. S7,† the 3EHT composition could be represented by the singlet appearing at 6.94 ppm while 3HT signal shows up at 6.98 ppm. By taking the ratio of the two thiophene hydrogens, we obtain the composition change plot shown in Fig. 4(b). The linear increment of the 3EHT composition in the copolymer demonstrates the gradient distribution of monomers in the final product. We hypothesize that the different initiating dimers of 3EHT formed in copolymerisation and homopolymerisation may lead to contrasting results of the two polymerisations. Ni(dppp)Cl2 as the catalyst leads to the in situ formation of tail-to-tail dimers as the intiator.14 In homopolymerisation, only 3EHT-3EHT dimers would be formed as the initiation point of polymer chains. While in the copolymerisation, less sterically hindered dimers such as 3HT–3HT and 3HT–3EHT offer more efficient paths to add incoming monomers, which results in the successful copolymerisation of P(3HT-co-3EHT) at room temperature.
![]() | ||
| Fig. 4 (a) Respective monomer consumption behaviour of 3HT and 3EHT and (b) compositional change of 3EHT in P(3HT-co-3EHT). | ||
The close values of the reactivity ratios of 3HT and 3BrHT predict a random copolymer structure with the comonomer composition equal to the feed ratio. The consumption rate of both comonomers is close to each other as expected by the kinetic values. In Fig. S9,† the characteristic triplet at around 3.42 ppm represents the fraction of 3BrHT in the polymer. Taking the ratio of the triplet over the alpha methylene signal, we obtain the composition plot as shown in Fig. 5(b) with a stable 3BrHT composition throughout the synthesis. The result shows that the final copolymer has random monomer distribution.
![]() | ||
Fig. 5 (a) Respective monomer consumption behaviour of 3HT and 3BrHT and (b) compositional change of 3BrHT in P(3HT-co-3BrHT) (initial feed ratio was 0.65 : 0.48). | ||
:
3OPT
:
3OPT has shown difficulty given the decreased R2 value, we hypothesized that the deviation from the model resulted from the failed chain polymerisation for this pair. We validate the hypothesis by evaluating the dispersity change of 3OPT homopolymers and random copolymers as shown in Fig. 6. The homopolymerisation of P3OPT gives steady dispersity values of around 1.4 throughout the reaction and linear growth of the molecular weight with monomer conversion, as shown in Fig. 6(a). On the other hand, the dispersities of the equimolar random copolymerisation of 3HT and 3OPT have increased significantly with the apparent loss of linear relationship between MW and monomer conversion [Fig. 6(a) and (b)]. The other combinations showed clear linearity between DP and monomer conversion (Fig. S10–S12†). This implies that the loss of the chain polymerisation comes from the difference in the structures of the two comonomers.
Further qualitative investigation was performed to understand the kinetics of hetero-monomer addition in the failed KCTP of 3HT
:
3OPT. The syntheses of diblock copolymers of P(3HT-b-3OPT) and P(3OPT-b-3HT) were attempted and characterized with MALDI-TOF. As shown in Fig. 7(a), m/z values of the most abundant peaks reflect the growth of 3OPT blocks from the 3HT chain ends. However, in the case of P(3OPT-b-3HT), we observed very limited growth of 3HT segments from the 3OPT chain ends (Fig. 7(b)). Moreover, the MALDI spectrum of P(3HT-b-3OPT) also implies the loss of end group control as Br/Br chain ends become predominant, as indicated in Fig. 7(a). These observations possibly result from the interfered association of the 3OPT monomer with the Ni catalyst due to the phenyl group.37
As a final point, we address the differences observed in the polymerization of 3EHT and 3OPT. Both are bulkier than the other monomers. 3EHT fails to homopolymerize, while 3OPT does. 3EHT successfully copolymerizes while 3OPT does not. We speculate that these differences may arise from the differences in reactive species present due to the Schlenk equilibrium, where the alkyl thiophenes may form more dimers exasperating steric issues for 3EHT,38 or because the additional phenyl unit on 3OPT alters the association of the Ni catalyst.37 These issues may also affect the differences in monomer conversion, as shown in Fig. 3–5(a). Teasing out these details will be a topic of future study.
:
3OPT. Our systematic study thus indicates the need for special handling of experimental conditions in the cases of bulky comonomers, such as temperature, catalyst ligands,37 and choice of Grignard reagents.39 Understanding the impact of the monomer structure on the reaction kinetics allows the efficient selection of appropriate comonomers and experimental conditions without tremendous commitment to the trial-and-error process.
Footnote |
| † Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4py00009a |
| This journal is © The Royal Society of Chemistry 2024 |