Open Access Article
Matt J. Pricea,
Anna Bachs-Herrera
bc,
Arianna Brandolese
a,
Yuya Watanabe
a,
Matthew N. Grayson
*b,
Antoine Buchard
*c and
Andrew P. Dove
*a
aSchool of Chemistry, University of Birmingham, Edgbaston, Birmingham, B15 2TT, UK
bDepartment of Chemistry, University of Bath, Claverton Down, Bath, BA2 7AY, UK
cDepartment of Chemistry, Green Chemistry Centre of Excellence, University of York, York, YO10 5DD, UK
First published on 15th April 2026
Polyethylene terephthalate (PET) is one of the most prominent single-use plastics. It can be readily circularised via chemical depolymerisation to monomer and subsequent repolymerisation to deliver virgin quality PET from waste. Chemical recycling-to-monomer (CRM) approaches rely on catalysts to enhance the activity of the depolymerisation process. Recently, dual-catalytic systems that combine inexpensive and readily available Lewis acids with simple organobases have displayed cooperative catalytic activity, among the highest reported for these processes; however, their mechanism of action, alongside the large number of possible combinations, makes identifying ideal catalyst species challenging. In this manuscript, we systematically explore the relationship between metal Lewis acids and organobases and their impact on PET glycolysis activity. These studies reveal a relationship between metal acetate electronegativity and the pKa of the organobase that can be used to predict the most active combination for a given metal salt or organobase. Furthermore, through preliminary DFT calculations, we demonstrate that the glycolysis likely follows a mechanism in which the base coordinates to the Lewis acid, acting as a single catalytic entity, and in which the acetate ligand assists with proton transfer.
Where PET is collected, the waste processing industry applies mechanical recycling and reprocessing methods that can result in a product that can be used for the same initial purpose. However, over time, the properties of the polymer degrade, and it becomes contaminated with dirt and other debris that impact both its performance and appearance, leading to “downcycling” into lower-value products. Chemical recycling-to-monomer (CRM) offers a valuable alternative to mechanical recycling, enabling the depolymerisation of waste PET back to monomers from which recycled polymers that possess the same properties as virgin material can be made.
To date, PET chemical recycling has been exploited through hydrosilylation,8–10 hydrogenolysis,11 and recently by both glycolysis12–16 and alcoholysis17,18 processes; the latter approaches present the most viable solutions for depolymerisation/repolymerisation cycles at scale.12,19–26 Most focus has been directed towards PET glycolysis to bis(2-hydroxyethylene)terephthalate (BHET) using both metal-based and organic catalysts.22,27,28 The requirement to operate at high temperatures in the presence of protic sources makes catalyst design challenging, hence simple metal salts, such as Zn(OAc)2 are among the most widely studied, as well as being used in industrial PET glycolysis processes. Ionic liquids have also been widely studied and employed as recyclable catalysts, including in commercial settings, offering good depolymerisation activity and separation potential.29–31 More recently, organic salt catalysts, made from organic acids and bases, have also been applied for PET glycolysis.28 These salts have been almost all found to be more thermally stable than their component parts, and hence remain more highly active over longer time periods (or catalyst recycle loops) compared to the base catalysts alone.22,32,33 Several other studies, including our own, have shown a significant rate enhancement when metal acetates (e.g. zinc(II), manganese(II), copper(II) and nickel(II)) are combined with organic bases (e.g. tropine, 4-dimethylaminopyridine, DMAP, etc.), with high catalyst stability.14,22,34–36 Notably, the combination of species results in a synergistic effect where the activity of the combination is greater than the sum of the parts, proposed to be a result of a dual catalytic mechanism in which the Lewis acid and the base act as separate entities, but in tandem.
The large number of potential combinations of inexpensive, commercially available and earth-abundant metal salts and organic bases that could act as PET glycolysis catalysts presents significant potential to enhance catalyst activity in a practical manner. In turn, this could enable depolymerisation to occur in a more energy-efficient manner, for example, at a lower temperature, or to enable access to a wider range of alcoholysis reactions that might, at present, be excluded as a consequence of the high reaction temperatures applied. To this end, to advance catalyst design towards these goals, we sought to develop a fundamental understanding of the proposed dual catalytic mechanism through experimental and computational methods and hence elucidate the key features of the catalyst components that allowed for predictable selection of highly active catalytic combinations of metal salts and organic bases.
:
1 combination of Lewis acids and organobases for PET depolymerisation at 180 °C, with 1 mol% catalyst loading with respect to each component (Fig. 1A). As expected, lower catalyst loadings resulted in reduced activity, with 1 mol% representing the optimal compromise between catalyst loading and reaction time (Fig. S17). Initially, organobases with pKaH ranging from 5 to 11 were screened using Zn(OAc)2 as the metal salt catalyst. Depolymerisation was assessed by monitoring the appearance of BHET in solution using 1H NMR spectroscopy of samples taken from the reactions, which contain N-methyl-2-pyrrolidone (NMP) as an internal standard. The initial rate constants, k′, were derived from the conversion-time plots (Fig. S8–S11) and analysed using a ‘shrinking core’ model (eqn (S1)) on account of the heterogeneous nature of the depolymerisation reaction. Analysis of the initial rate constants with respect to pKaH of the organobase revealed no clear trend, with a peak of activity appearing for the combination of Zn(OAc)2 with imidazole, Im (pKaH = 6.9), a weak base, and lowest activity when 1,4-diazabicyclo[2.2.2]octane, DABCO (pKaH = 8.8) was used as the cocatalyst.
In search of more active catalyst combinations, we sought to develop an understanding of how the activity trend differed when the bases were combined with other metal acetates for catalysing PET glycolysis. To survey the impact of Lewis acid, three further metal acetates, La(OAc)3, Mg(OAc)2, Mn(OAc)2, were selected on account of their wide range of electronegativity (as determined by Pauling electronegativity (χ) values) and studied for depolymerisation of PET with each of the organic bases applied previously with Zn(OAc)2. Determination of the initial depolymerisation rate for each metal salt/base combination revealed that a different base provided peak activity for each different metal salt (Fig. 1B, Fig. S8–S11). La(OAc)3/base dual catalysts display a relatively clear trend of increased base pKaH correlating to increased activity (Fig. S10), likely a consequence of base dominated activity on account of the relative weakness of La(OAc)3 as a Lewis acid and poor activity as a glycolysis catalyst, with peak activity achieved with 1,8-diazabicyclo(5.4.0)undec-7-ene, DBU (pKaH = 13.5). Similar analysis for Mg(OAc)2 and Mn(OAc)2 as the Lewis acid revealed peak activities with 1-azabicyclo[2.2.2]octane, ABCO (pKaH = 11.0) and N-methyl imidazole, NMI (pKaH = 7.4), respectively. An analogous correlation is observed when considering the bases pKaHs values in MeCN and the metal acetate electronegativity (Fig. S12).
To explain the combinations that led to peak activity for all the catalyst combinations, we examined possible relationships between characteristics of the Lewis acid (Pauling electronegativity, ionic radius, Pearson hardness and Gagne–Hawthorne length) and organic base pKaH. Interestingly, plotting the Pauling electronegativity (χ) of the metal and the pKaH of the organic base for the combinations that displayed the highest activity for depolymerisation revealed a clear linear trend in which the metal acetate with the highest Pauling electronegativity was most active with the base with the lowest pKaH and conversely, the metal acetate with the lowest Pauling electronegativity was most active with the base with the highest pKaH (Fig. 1C and S12). To verify this observation, scandium acetate (Sc(OAc)3) was selected as it has a Pauling electronegativity value of χ = 1.36, between that of Mg(OAc)2 and Mn(OAc)2. The correlation predicted that a base with pKaH around 10 would provide the combination of peak activity, which was confirmed with depolymerisation tests using Sc(OAc)3/base combinations with DABCO, DMAP and ABCO in which the peak activity was achieved with Sc(OAc)3 and DMAP (base pKaH = 9.6), which falls closest to the linear relationship between peak k′, pKaH and χ (red sphere, Fig. 1C). To elucidate the catalyst structure, additional studies were conducted to investigate the potential coordination of ethylene glycol and NMP to the metal acetate (Fig. S13–S16). Zn(OAc)2 was examined in the presence of varying ratios of ethylene glycol and NMP, and the mixtures were monitored by NMR spectroscopy. Only negligible spectral shifts were observed for Zn(OAc)2, indicating that ethylene glycol and NMP are unlikely to coordinate to the zinc centre under the glycolysis conditions tested. Additionally, all tested catalysts were insoluble in ethylene glycol at 180 °C (Fig. S18).
To gain a deeper insight into the role of each single catalyst component for PET glycolysis and confirm that the peak dual catalyst combinations retained the synergistic effect in which the combined activity was greater than the sum of the catalytic activity of the individual components, we screened PET depolymerisation activity of the Lewis acids and organic bases as standalone catalysts. Simply surveying PET depolymerisation activity when catalysed by the range of organic bases studies reveals a simple trend that relates initial depolymerisation rate to pKaH, in which the activity generally increases as the pKaH increases, in line with stronger bases being better at activating the ethylene glycol towards nucleophilic attack on the PET (Fig. S3). Screening metal acetate catalysts found that the activity was in the order of Mn(OAc)2 > Mg(OAc)2 > Zn(OAc)2 ≈ La(OAc)3 (Table S2). While Zn(OAc)2 is most commonly used as a catalyst,37 these results agree with previously reported works, which place Mn(OAc)2 as more efficient than Zn(OAc)2.38 Correlation between activity and metal acetate physical parameters was conducted through the aid of a linear regression analysis of the proxy variables with activity. Regression analysis showed an inversely proportional relationship with both electronegativity and ionic radius, implying that smaller, less electronegative metal ions have higher activity. This indicates that what defines ultimate catalyst activity may be more complicated than simply evaluation of the metal acetate electronegativity, in line with the observation that the most active dual catalyst system was Mn(OAc)2/NMI. Notably, however, in each case, the combination of metal acetate and organic base displayed between 1.5 and 3× higher k′ than the sum of the metal acetate and organic base k′ values (Table S2 and Fig. S2–S5). This strongly indicates that the combination of metal acetate and organic base results in a more active catalyst species than could be obtained from both species acting independently.
The property-activity relationships suggest that the origins of the synergistic effect between metal acetate and organic base may be more complex than initially assumed in the dual catalyst mechanism (Fig. 2). Therefore, Density Functional Theory (DFT) calculations were undertaken using Gaussian 16, Revisions C.01 (ref. 39) and C.02 (ref. 40) to elucidate the more intricate relationships between mechanism, properties and activity. As outlined previously, we had anticipated that interactions between the metal acetate and organic base would lead to the base being bound to the metal centre at room temperature, but to unbind at the reaction temperature of 453 K where it could undertake a dual catalyst mechanism. The possibility of a complex formation between the metal acetate and base was therefore explored by DFT calculations, focusing on Zn(OAc)2 as a representative system. Calculation of the binding Gibbs free energies (ΔGbinding) of Zn(OAc)2 and each of the organic bases employed experimentally revealed a reversible process with ΔGbinding ≈ 0 for bases with pKaH between 6.9 and 11. At the extremes of pKaH a clear thermodynamic preference for being bound or unbound was observed. Binding of DMA (pKaH = 5.0) was unfavourable (ΔGbinding = +10.8 kcal mol−1), whereas, in contrast, interaction with the strong base, DBU, (pKaH = 13.5) revealed that the binding of the base to metal was thermodynamically favourable, ΔGbinding = −4.2 kcal mol−1. ΔGbinding for bases with pKaH between these extremes were closer to zero (Fig. 3). This finding suggests that base strength may influence the mechanism by which depolymerisation occurs, with stronger bases interacting more strongly with Zn2+, although steric factors may also play an important role in the binding of the two species.
![]() | ||
| Fig. 2 Schematics of the two proposed mechanisms of PET glycolysis: bound (left) and unbound (right) pathways. | ||
![]() | ||
| Fig. 3 Zn(OAc)2/base complex geometries and binding Gibbs free energies. The systems were optimised at ωB97X-D/6-31G(d) level of theory and with single point energy calculation at the ωB97X-D/def2-TZVP level of theory, both with the CPCM implicit solvent model for 1,2-ethanediol at 453 K (ε = 19.1).45 Images were created using CYLview 1.0 beta.46 | ||
The observation of binding Gibbs free energies being close to zero with most bases suggests that two possible mechanistic scenarios are possible: one in which the base coordinates to the Lewis acid, acting as a single entity (from now on named bound pathway), and another in which they act separately (unbound pathway) with the base acting as a Brønsted base (Fig. 2). To gain further insight into the catalyst's mode of action, both pathways were explored using small molecule model compounds. Methyl benzoate was chosen to mimic the polymer substrate, along with ethanol, as a mimic for ethylene glycol. Zn(OAc)2/Im was chosen as the benchmark dual catalyst. For both envisaged pathways, a conformational search for intermediates and transition states (TSs) was carried out using CREST, version 2.1241,42 with GFN2-xTB,43 followed by DFT optimisation and single point energy calculations, and final computation of the quasi-harmonic Gibbs free energies with GoodVibes.44 The large difference in the activation barriers showed in the energy profiles of the glycolysis of PET catalysed by Zn(OAc)2/Im (ΔΔG‡ = +7.9 kcal mol−1, Fig. 4) suggested that the reaction is likely to follow the bound pathway, where the proton transfer occurs with one of the acetates (TSIII–IVa).
![]() | ||
| Fig. 4 PET glycolysis catalysed by Zn(OAc)2/Im in the bound (green) and unbound (orange) pathways: calculated Gibbs free energy profile at 453.15 K. The structures were optimised at ωB97X-D/6-31G(d) level of theory and with single point energy calculation at the ωB97X-D/def2-TZVP level of theory, both with the CPCM implicit solvent model for 1,2-ethanediol at 453 K (ε = 19.1).45 | ||
The influence of the base on the catalytic mechanism was next probed by investigating the TSs of the bound and unbound pathways for catalytic systems Zn(OAc)2/DMA, Zn(OAc)2/ABCO and Zn(OAc)2/DBU, which feature bases with a lower (5.0) and greater (11.0, 13.5) pKaH, respectively, than that of imidazole (6.9). These TSs were modelled from the lowest energy conformers of the Zn(OAc)2/Im system, by replacing Im with DMA, ABCO and DBU (see SI). After DFT optimisations and single point energy calculations, the Gibbs free energy barriers for the two pathways were plotted for each dual catalyst (Fig. 5).
![]() | ||
| Fig. 5 Gibbs free energy barriers (bound and unbound pathways) as a function of the Zn(OAc)2/base systems (base = DMA, Im, ABCO, DBU). The systems were optimised at ωB97X-D/6-31G(d) level of theory and with single point energy calculation at the ωB97X-D/def2-TZVP level of theory, both with the CPCM implicit solvent model for 1,2-ethanediol at 453 K (ε = 19.1).45 | ||
Although no correlation was found between the energy barriers and the pKaH of the base or the ΔGbinding of the dual catalysts, the bound pathway is favoured over the unbound pathway for the four bases. Furthermore, Zn(OAc)2/Im presents a comparable barrier to that of Zn(OAc)2/ABCO (34.1 and 33.9 kcal mol−1, respectively) despite the higher activity of the former, while the bound pathway for Zn(OAc)2/DBU presents a higher barrier than Zn(OAc)2/Im. At the experimental temperature of 453.15 K, the differences in the kinetics observed result in differences of Gibbs free energies of the order of 0.1 kcal mol−1, which DFT is unable to reproduce. For this reason, an explanation for the faster PET glycolysis rate promoted by Zn(OAc)2/Im cannot be provided at this stage.
The DFT data collected subsequently allowed us to formulate a plausible reaction mechanism. In the first step, the zinc ion adopts a square pyramidal geometry, with imidazole coordinated in axial position, and both acetate ligands coordinated in a η3 fashion and bidentate (Fig. 4, II). This is followed by the zinc ion returning to a tetrahedral geometry upon coordination of the ethanol, with one of the two now monodentate η1 acetate ligands, activating ethanol by hydrogen-bonding, while some π–π stacking happens between the imidazole and the PET model (IIIa, Fig. S23). In the transition state (TSIII–IVa), the carbonyl oxygen of the PET coordinates to the zinc ion, which maintains the tetrahedral geometry, replacing the ethanol. The ethanol is directed towards the carbonyl's electrophilic centre, attacking it as it is deprotonated by the acetate. This results in high-energy intermediate IVa, in which the tetrahedral carbon readily rearranges into a new carbonyl and an ethyl ester (V and VI). However, the higher-than-expected energy barriers calculated for Zn(OAc)2/DMA (compared to experimental ones) suggest that this dual catalyst may follow a different depolymerisation mechanism. Future work will aim at exploring other possible mechanisms, including the involvement of ethylene glycol as a ligand in the catalyst.
Supplementary information: detailed experimental procedures and additional characterisation data; DFT calculation details, data and associated digital repository. See DOI: https://doi.org/10.1039/d5cy01499a.
| This journal is © The Royal Society of Chemistry 2026 |