Open Access Article
Peter
Deglmann
*a,
Sara
Machleit
b,
Cesare
Gallizioli
b,
Susanne M.
Rupf
b and
Alex J.
Plajer
*b
aBASF SE, Carl-Bosch-Strasse 38, 67056 Ludwigshafen am Rhein, Germany. E-mail: peter.deglmann@basf.com
bInstitut für Chemie und Biochemie, Freie Universität Berlin, Fabeckstraße 34-36, 14195 Berlin, Germany. E-mail: plajer@zedat.fu-berlin.de
First published on 10th April 2023
The catalytic construction of well-defined materials from mixtures of building blocks is an important challenge in sustainable catalysis. In this regard, we have recently reported a new type of selective ring-opening terpolymerisation (ROTERP), in which three monomers (A, B, C) are selectively enchained into a (ABA′C)n sequence, but the reasons behind this unusual selectivity remained unanswered. In this study, we present a detailed investigation into the full ROTERP mechanism based on the reactivity of model intermediates, computational studies investigating >100 possible intermediates and transition states and reaction kinetics. Experimental verification of the intermediate speciation, the primary insertion steps and the side-reactions lets us show that although most insertions and side-reactions are thermodynamically viable, kinetic selection processes at the propagating chain end determine the sequence selectivity. Computational studies elucidate the special role and speciation of the lithium catalyst which during the catalytic cycle involves mono-metallic, bi-metallic and charge separated transition states comprising both coordinative activation of incoming monomers and functional groups of the polymer backbone adjacent to the propagating chain. Our study not only deciphers the mechanism of a rare selective terpolymerisation but also helps answering open questions relevant to ring-opening copolymerisation (ROCOP) and alkali-metal catalysis in general, thus guiding the design of future polymerisation catalysis for degradable materials.
O)SLi] = 216.0 ppm, δ [R(C
O)OR] = 169.2 ppm,
(C
O) = 1701.5 cm−1) and with CS2 generating lithium dithiocarbonate DT (δ [RO(C
S)SLi] = 244.6 ppm,
(C
S) = 1037.1 cm−1) as confirmed by multinuclear NMR, IR and a yellow discolouration due to formation of the (C
S) chromophore. In contrast no reaction is observed between A and BO (which we employed as the model epoxide throughout the study due to its higher boiling point compared to PO) showing that epoxide insertion by the lithium alkoxide is associated with significantly higher barriers explaining the absence of ether linkages during ROTERP. However, this leaves the question unanswered why PTA is selected over CS2 by A in ROTERP and this could either be due to insertion kinetics or intermediate thermodynamics. In terms of the latter, we found that CS2 insertion into A to form DT is at least partially reversible. Heating of DT under dynamic vacuum at 80 °C for 1 min results in formation of A as visually indicated by a vanishing of the yellow colour and confirmed by 1H and 7Li NMR. However, some insoluble decomposition products also form. Analogous attempts to revert PTA insertion into A starting from TC at high temperatures failed. If CS2 into A is reversible but PTA insertion isn't, we hypothesized that transformation of DT into TC should be feasible. Addition of PTA to DT in d8-THF initially results in no reaction at room temperature but heating to 80 °C overnight results in complete transformation to TC. This clearly indicates that PTA insertion from ATT is the thermodynamically favoured propagation event. However, as ROTERP also occurs at room temperature with identical sequence selectivity and DT cannot be converted to TC at room temperature intermediate thermodynamics cannot be the only reason behind the insertion selectivity of ATT. To determine the kinetically favoured propagation step, we conducted a competition experiment in which we reacted A with a 10
:
1 mixture of CS2 and PTA at room temperature in d8-THF (note that a 10 CS2
:
1 PTA ratio is also employed in the required monomer mixtures of ROTERP). Interestingly we see 96% selectivity for PTA insertion yielding TC and only 4% CS2 insertion yielding DT although an excess of CS2 was provided. We hence propose that lithium alkoxide ATT selects PTA on kinetic grounds. Nevertheless, the same insertion is likewise favoured thermodynamically, and CS2-PTA exchange could play a role at increased reaction temperatures.
O)SLi] = 216.0 ppm,
= 1689.1 cm−1 and δ [R(C
O)SR] = 193.7 ppm,
= 1556.3 cm−1) and with CS2 generating the lithium trithiocarbonate TT (δ [RS(C
S)SLi] = 244.6 ppm,
= 996.5 cm−1) as confirmed by multinuclear NMR, IR and a yellow discolouration due to the C
S chromophore. Since both CS2 and PTA insertion are again possible from T this raises the question whether monomer selection in this step during ROTERP occurs on a thermodynamic or kinetic basis. Again, we found CS2 insertion to be reversible. Heating TT at 80 °C under vacuum cleanly regenerates T as observed by 1H and 7Li NMR which is also visible by the disappearance of the yellow (C
S) chromophore. Reversible CS2 insertion could again provide a pathway for CS2 to PTA exchange. However, in contrast to DT, this reaction does not occur cleanly for TT and results in complex product mixtures. Although some PTA ring opening is observed indicating this insertion to be thermodynamically favoured over CS2 insertion, monomer exchange pathways are unlikely to play a role during ROTERP and we thus investigated the potential for kinetic selection. Reacting T with a 10
:
1 mixture of CS2 to PTA, which is the monomer ratio present during ROTERP, results in 99% selective CS2 insertion yielding TT. We hence suggest the monomer selection in step (5) of Fig. 1 to occur on a kinetic basis; note that this is the opposite chemoselectivity that A shows. Surprisingly we observe ring opening of BO by T to form thioethers in d8-THF at room temperature. However only trace amounts of thioether links are formed during ROTERP and this is likely because of faster CS2 than BO insertion. Accordingly reacting T with a 2
:
1 mix of CS2 and BO resulted in no thioether and only heterocarbonate formation. Another aspect concerning TE is the origin of erroneous thioester links (Fig. 1) which either form from TE if it inserts PTA or if ATE propagates without O/S exchange. However, we previously found O/S exchange from TC to always occur quantitatively after BO insertion, making this thioester formation pathway unlikely.49 On the other hand, we observed that lower CS2 loadings in the starting monomer feed result in more thioester errors making PTA over CS2 insertion from TE the more likely mechanistic origin of errors. Seeking to confirm this hypothesis we conducted are series of competition experiments like those depicted in Fig. 3b in which T was combined with mixtures of 1 eq. PTA and 5 eq., 2 eq. or 1 eq. CS2 in d8-THF. In contrast to the competition experiment featuring a 1
:
10 PTA
:
CS2 ratio, we observe PTA ring opening to an increasing extent as the amount of CS2 is lowered. Hence kinetic CS2 over PTA selection by the thiolate T is not achieved at lower relative CS2 concentrations making erroneous PTA insertion by the thiolate intermediate T the likely origin of thioester errors.
S)SR] = 224.1 ppm,
= 1054.0 cm−1) and an ester link (δ [RC(
O)OR] = 168.5 ppm,
= 1700.0 cm−1) as well as a lithium thiocarboxylate chain end (δ [RC(
O)SLi] = 217.1 ppm) and tail-selectively ring-opened BO. Going one step further, all selection processes occurring in steps (1) to (5) (Fig. 1) can be coupled to each other as shown in Fig. 5. Here, we reacted T with a mixture of 10 eq. CS2, 1 eq. PTA and 1 eq. BO to likewise form the thiocarboxylate TC* while preserving the unscrambled trithiocarbonate link.
:
1 mixture of TC and TT with 1 eq. BO in d8-THF/CS2 at room temperature. Under these conditions negligible reaction of TC occurs whereas TT reacts to form organic heterocarbonates (Fig. 6(b)) by BO insertion and initiation of CS2/BO coupling. We then compared the reaction kinetics of PTA/CS2/BO ROTERP with that of the parent PTA/BO and CS2/BO ROCOP reactions which complemented our findings (Fig. 6(c)). Here we tracked the reaction progress of Li catalysed PTA/CS2/BO ROTERP (at 1 eq. LiOBn: 100 eq. PTA: 800 eq. BO: 1600 eq. CS2), CS2/BO ROCOP (at 1 eq. LiOBn: 800 eq. BO: 1600 eq. CS2) and PTA/BO ROCOP (at 1 eq. LiOBn: 100 eq. PTA: 800 eq. BO) employing an equal volume of toluene in place of CS2 to reach comparable reactant concentrations. The reaction kinetics at room temperature assessed by 1H NMR aliquot analysis as shown in Fig. 6(c) reveal that ROTERP (TOF ca. 59 h−1) is slower than CS2/BO ROCOP (ca. 234 h−1) but faster than PTA/BO ROCOP (ca. 13 h−1) under comparable reaction conditions with regards to BO turnover. This suggests that reaction steps during ROTERP like PTA/BO ROCOP (i.e. PTA insertion and BO insertion from TC) are slower than the ones like CS2/BO ROCOP (i.e. CS2 insertion and BO insertion from TT). As ROTERP is 0th order in PTA this leaves BO insertion from TC as the slowest reaction step as indicated by our prior model intermediate studies. We therefore propose that the intermediate TC is the likely resting state of ROTERP. Furthermore, the experiment outlined in Fig. 5 can be viewed as a ROTERP at high Li-initiator loading (i.e. 100 mol%) as all monomer are present which are enchained to the TC* oligomer. Hence this provides direct spectroscopic evidence for the thiocarboxylate chain-end in the resting state of the ROTERP which confirms the conclusions made in this paragraph.
![]() | ||
Fig. 6 (a) Linear fit of the PTA conversion vs. time at 1 eq. LiOBn: 100 eq. PTA: 800 eq. BO: 1600 eq. CS2 at RT. (b) Competition experiment reacting a 1 : 1 mixture of TC and TT with 1 eq. BO at room temperature in d8-THF (overlaid 400 MHz 1H NMR shown); note that the previously reported tBu derivative of TC was employed to avoid signal overlap.49 (c) Comparison of the BO turnover versus time for ROTERP and the two parent ROCOP reactions as outlined in the main text; TOF as discussed in the main text obtained by linear fitting over the entire displayed range. | ||
:
1 mixture (by molarity) of CS2 and PO; a temperature of 80 °C was assumed for this solvation treatment as well as for computing the statistic thermodynamics contributions to allow for an optimal comparison with our initial work.49,59 For all transition states, it was checked which combination of reactants and products they actually link together by slight distortions (corresponding to 30 K) along both directions of the imaginary vibrational mode and subsequent structural relaxation. Before entering into the discussion of our computational results, we would like to point out that identification of reaction intermediates and transition states is non-trivial due to the manifold of possible ligands for coordination to the Li centres, which can be chain-ends, functional groups of the polymer chain adjacent to the chain-ends and the monomers. With respect to monomer coordination, it could be shown that PO (as a model alkylene oxide) binds more strongly than CS2 as well as PTA, so that in the results shown in the following always PO was used to saturate available coordination sites at Li. It became also obvious that phthalic diester groups of the product are in spite of their chelating nature not superior as ligands compared to two PO, so that coordination to remote functional groups of the propagating chain or groups of other polymer chains is not expected (see ESI†). Another complication represents the tendency of lithium salts to aggregate into dimers and even more aggregated structures in solution. The most favoured intermediates and transition states from a Gibbs free energy perspective were identified assessing over 100 potential species for all of which searches for the lowest energy conformer at the level of structure optimization were performed; the reader is referred to the ESI† for a summary of the structures and their associated Gibbs free energies. Despite of these efforts, our investigations cannot deliver an exhaustive list of all possible lithium complexes involved in the ROTERP cycle, but rather deliver one possible reactive pathway which nevertheless allows to draw some general conclusions about the transition state geometries, intermediate speciation, and the role of the lithium catalyst. A summary of the computed catalytic cycle is given in Fig. 7. The starting point represents the suspected lithium thiocarboxylate resting state in its monomeric form Im for which coordination of three additional epoxides results in the most stable monomeric structure. We used thiobenzoate as a model polymer chain end and propylene oxide as model alkylene oxide throughout this study for reasons of computational simplicity. Aggregation to Id results in a slightly more stable dimeric structure with a Li2O2 core bridged by the thiocarboxylate oxygens. It should be mentioned here, that for all aggregation equilibria, lithium thiocarboxylates were inferred as reagent, as these represent the suspected chain ends in the resting state of the catalytic cycle and must therefore be primarily available in the reaction mixture. From Id, ring opening occurs through a bimetallic transition state TSI (ΔG‡ = 125.4 kJ mol−1) in which the reacting thiocarboxylate is coordinated by one Li while the inserted PO ring is coordinatively activated by the second Li centre. It should be noted that tail-selective PO ring-opening at the CH2 position is energetically most favourable as it represents the SN2 reaction at the less substituted carbon atom. Monometallic PO ring-opening transition states starting from Im result in significantly higher ΔG‡ as a consequence of the angular constraints around the attacked C-atom at which in SN2 reactions nucleophile and leaving group optimally assume an angle of 180°. Further bimetallic transition states in which two monomeric units are not linked by an anionic bridging ligand are likewise less favoured as a consequence of the fact that either reactant or product have to be of zwitterionic nature. TSI leads to a dimeric intermediate IId with the thiocarboxylate bridging through the O and S centres (generally, μ2- or κ2-bridging by either O only or O and S, respectively, was found to be energetically quite close for thiocarboxylate links between two Li) as well as a bridging alkoxide. Notably the carbonyl oxygen of the adjacent thioester link coordinates in this most stable form of IId to one of the Li centres and this already hints the subsequent activation of the thioester group for nucleophilic attack in the O/S exchange step. Intramolecular attack (ΔG‡ = 48.4 kJ mol−1) of the thioester by the alkoxide viaTSII transforms the alkoxide into a thiolate and this exchange step is orders of magnitude-faster than the PO ring-opening before (and expected to be also faster than any reaction with CS2 or PTA, when comparing the analogous differences in Gibbs free energies between Vd and transition states TSV or TSV*). TSII actually represents the initial alkoxide attack to form as primary product a tetrahedral intermediate in which the two O-atoms as in the transition state are linked to Li; however, any attempt of optimizing analogous species with the negatively charged O and S coordinated to Li – as a prerequisite for thiolate formation – led to spontaneous C–S bond cleavage, so that it can be concluded that this carbonyl addition–elimination reaction is limited by the addition step and that elimination is fast and does not require more than dissociation of a neutral donor atom from Li and coordination of another. Notably the thiolate IIId is 24.4 kJ mol−1 more stable than its alkoxide counterpart IId demonstrating that O/S exchange at the chain-end is a thermodynamically favourable process. Monometallic O/S exchange transition states starting from monomeric lithium alkoxides lead to significantly higher barriers, which fits to the general observation that the strongly coordinating alkoxide anion has a high preference for dinuclear species which is much less pronounced for less coordinating anions like thiolate or anions with negative charge delocalization like thiocarboxylates or trithiocarbonates. Moreover, transition states similar to TSII for which the thiocarboxylate spectator ligand acts as a μ2-bridge via O only (instead of κ2-bridging via O and S) also lie energetically higher, which results from the fact that the corresponding structures contain too many small rings annulated with each other, showing that a modulation of the Li⋯Li separation through different bridging modes of the thiocarboxylates plays a role in the ROTERP mechanism. From IIId CS2 insertion has to occur in order to continue propagation, however in this case bimetallic transition states are less favoured than monometallic ones. We infer that this is partially due to steric encumbrance of the thiolates by the two flanking Li centres which also withdraw electron density from the nucleophilic S. In order to achieve monometallic CS2 insertion, fragmentation of the dimers IIId into their monomeric form IIIm under coordination of an additional equivalent of PO per Li occurs, overall requiring 21.7 kJ mol−1. The four-membered transition state TSIII (ΔG‡ = 45.1 kJ mol−1 with respect to IIIm) that leads to CS2 insertion also lies from a free energy perspective lower than the prior PO ring-opening. Importantly CS2 insertion is only associated with a slightly lower barrier than PTA insertion (TSIII*, ΔG‡ = 46.5 kJ mol−1) confirming that although CS2 insertion is kinetically favoured, monomer selection by the lithium
thiolate chain end is also concentration driven and requires a large excess of CS2 as found in our experimental studies in order to achieve highly regular product structures. These findings also explain that erroneous PTA insertion from lithium thiolates viaTSIII* cause the aforementioned thioester errors. Progressing in the catalytic cycle, TSIII leads to the lithium trithiocarbonate IVm which then undergoes the thermodynamically slightly favourable dimerization to IVd with a thiocarboxylate Im. Surprisingly from here the energetically most favourable PO ring-opening pathway involves initial dissociation of the trithiocarbonate ligand under additional PO coordination to form the zwitterionic intermediate IVz (the given free energy refers to then solvent separated ion pair) which lies energetically uphill by 87.7 kJ mol−1 compared to IVd. This ligand exchange at Li effectively swaps a μ2-S for a μ2-O donor and could be facilitated by the high oxophilicity of Li. PO ring-opening in TSIV from here occurs with a very small additional activation barrier (ΔG‡ = 24.4 kJ mol−1) so that the overall barrier for PO ring-opening starting from IVd that needs to be overcome amounts to 112.1 kJ mol−1. PO ring-opening by the lithium trithiocarbonate from IVd to Vd is therefore predicted to occur more easily by 13.3 kJ mol−1 than PO ring-opening by lithium thiocarboxylates, confirming TSI to be the rate-determining step of the reaction. The zwitterionic nature of the most favourable TSIV can be understood as a consequence of the generally weaker binding between Li and trithiocarbonate anions compared to thiocarboxylates (as can be seen in the ESI by comparing the speciation of I and IV). From here, several possible propagation pathways with PTA are possible; in case of the monometallic ones some of the most favourable ones are actually barrierless, i.e. are spontaneous reactions once the corresponding monometallic precursor with a free coordination site has been generated by the energetically uphill (ΔG = 40.0 kJ mol−1) dissociation of dimers Vd into monomers Vm (see ESI† for more details) and this is in line with the experimentally determined 0th reaction order in PTA concentration; in cases with more steric crowding around Li, PTA inserts in a four-membered transition state TSV, which however does not seem the most favoured option. In contrast, CS2 insertion occurs most easily starting directly from Vd with an insertion barrier to TSV* of ΔG‡ = 69.2 kJ mol−1 confirming that PTA insertion is in fact kinetically favoured from this alkoxide intermediate that does – because of the more different barriers involved – less depend on monomer concentration than the competition of TSIII and TSIII*. One final question concerns the circumstance why no O/S exchange occurs from Vm (corresponding to ATT in Fig. 1). In this regard a cyclised variant Vc observed without explicitly searching for it upon a conformational scan of Vm suggests that cyclization initiating the O/S exchange process is a possible (and obviously reversible) process that occurs with a negligible barrier. In two cases it was found, that releasing a dialkyl dithiocarbonate and a thiolate chain end from this is actually slightly exergonic with respect to dialkyl trithiocarbonate and alkoxide chain end. However, from the activation barriers computed for TSIII, one would expect the corresponding transition states of CS2 addition around 20 kJ mol−1 higher (see ESI† for V) than the most favourable TSV. As PTA insertion is irreversible (vide supra), any Vm is thus in the presence of PTA rapidly scavenged which suppresses any O/S scrambling as predicted by our model intermediate studies.
![]() | ||
| Fig. 7 Energy profile for the ROTERP propagation cycle computed on the M06-2X/def2-TZVPD‖TPSSH/def-TZVP/COSMO(ε = ∞) level of theory. | ||
Footnote |
| † Electronic supplementary information (ESI) available. CCDC 2236905 and 2236906. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d3cy00301a |
| This journal is © The Royal Society of Chemistry 2023 |