 Open Access Article
 Open Access Article
      
        
          
            Maria 
            Schlüter
          
        
      ab, 
      
        
          
            Ryota 
            Enomoto
          
        
      a, 
      
        
          
            Shin 
            Makino
          
        
      a, 
      
        
          
            Lisa 
            Weihs
          
        
      c, 
      
        
          
            Cyra Lina 
            Stamm
          
        
      a, 
      
        
          
            Kerstin 
            Wohlgemuth
          
        
       b and 
      
        
          
            Christoph 
            Held
b and 
      
        
          
            Christoph 
            Held
          
        
       *ac
*ac
      
aTU Dortmund University, Department of Biochemical and Chemical Engineering, Laboratory of Fluid Separations, 44227 Dortmund, Germany. E-mail: christoph.held@tu-dortmund.de;   Tel: +49 231 755 2086
      
bTU Dortmund University, Department of Biochemical and Chemical Engineering, Laboratory of Plant and Process Design, 44227 Dortmund, Germany
      
cTU Dortmund University, Department of Biochemical and Chemical Engineering, Laboratory of Thermodynamics, 44227 Dortmund, Germany
    
First published on 5th September 2024
Glycolysis is the most promising chemical recycling method to depolymerize poly(ethylene terephthalate) (PET) with ethylene glycol (EG) into the monomer bis(2-hydroxyethyl) terephthalate (BHET). Boosting the depolymerization kinetics while staying under comparatively mild and green reaction conditions is required to bring glycolysis to industrial scale utilization. This work suggests achieving this goal by a combined pressure, temperature and co-solvent addition approach. By using the environmentally friendly γ-valerolactone (GVL) as a suitable co-solvent in the traditional PET glycolysis system, and slight temperature and pressure elevation, the kinetics was boosted by almost two orders of magnitude compared to the standard literature process. A kinetic model was employed to describe the kinetics as a function of temperature and GVL concentration. The optimized condition allowed nearly full conversion after 2 minutes only.
Nearly full conversion is achievable in PET glycolysis above 170 °C,20,21 but reaction kinetics is one of the bottlenecks to bring PET glycolysis towards a commercial scale. Reaction kinetics is influenced by the type and concentration of catalyst, PET particle size, concentration of reactants and temperature as well as using co-solvents that serve as kinetic activators. The homogeneous metal acetate catalyst zinc acetate (ZnAc2) has proven to be a reliable depolymerization catalyst, providing good selectivity towards BHET.21–24 The most common studies investigate PET glycolysis at 190 °C, since higher temperatures will cause significant EG evaporation. Nevertheless, there are some studies at temperatures up to 275 °C, which was achieved by pressurizing the reaction mixtures.25,26 Another strategy to increase kinetics is the addition of co-solvents to the initial reaction mixture, such as aniline, nitrobenzene, 1-methyl-2-pyrrolidinone, and dimethyl sulfoxide (DMSO) or a 1,3-dimethylurea (1,3-DMU)/ZnAc2 deep eutectic solvent system.27,28 Recently, Le et al.29 reported beneficial effects of the green solvent anisole on the kinetics.
Based on a reference kinetic dataset, which we determined for the ZnAc2-catalysed PET glycolysis at 190 °C and 1 bar, the objectives of this study are to: first, improve kinetics by using the biomass-derived green platform molecule γ-valerolactone (GVL)30 as additional co-solvent in the traditional reaction mixture. Moreover, individual experiments were performed to understand the kinetic boost of GVL in detail. Second, increase kinetics by elevating the temperature without solvent loss, which we will achieve by increasing the pressure. Third, make use of the combined effect of GVL as co-solvent and a temperature increase to ultimately boost the reaction kinetics. The influence of these individual factors on the reaction kinetics was evaluated by kinetic constants, obtained by the pseudo-first order reversible reaction model, as utilized in other studies.20,31,32 This combined temperature and co-solvent optimization approach will achieve a kinetic boost of PET glycolysis of about two orders of magnitude over the conventional process, thus contributing to a greener chemical recycling strategy.
| Chemical | Abbreviation | CAS | M/g mol−1 | Supplier | Type | Purity/wt% | 
|---|---|---|---|---|---|---|
| a Mass average molecular weight determined by the ASTM D 4603 standard test method in a previous work.33 | ||||||
| Poly(ethylene terephthalate) | PET | 25038-59-9 | 58 ![[thin space (1/6-em)]](https://www.rsc.org/images/entities/char_2009.gif) 000a | Colorless single-use post-consumer bottles (“JA”) | ||
| Ethylene glycol | EG | 107-21-1 | 62.07 | Sigma-Aldrich | ReagentPlus | 99 | 
| Zinc acetate | ZnAc2 | 557-34-6 | 183.48 | Sigma-Aldrich | 99.99 | |
| γ-Valerolactone | GVL | 108-29-2 | 100.12 | Sigma-Aldrich | ReagentPlus | 99 | 
| Water | H2O | 7789-18-5 | 18.015 | Deionized | ||
| Acetone | — | 67-64-1 | 58.08 | VWR international | TECHNICAL | 99 | 
| Bis(2-hydroxyethyl) terephthalate | BHET | 959-26-2 | 254.238 | Sigma-Aldrich | 98.1 | |
| Methanol | — | 67-56-1 | 32.04 | VWR international | HiPerSolv | 99.8 | 
The PET used in this study was sourced from colorless single-use post-consumer bottles. Caps and labels were taken off and surface impurities were removed using small amounts of acetone. The remaining material was manually cut into 5 × 5 × 0.2 mm pieces. Following this, the PET particles were shred into small particles using a grinder. These crushed PET particles underwent fractionation in a sieving tower. Unless otherwise stated, the experiments were carried out with a particle size fraction of 0.2 ≤ dp < 1 mm.
 was set to 50, and the initial mass ratio of solvent to PET
 was set to 50, and the initial mass ratio of solvent to PET  was maintained at 6. Finally, a small magnetic stirrer was added.
 was maintained at 6. Finally, a small magnetic stirrer was added.
      
      
        
        |  | (1) | 
|  | (2) | 
|  | (3) | 
|  | (4) | 
 = 50, 0.2 ≤ dp < 1 mm and a stirrer speed of 1400 rpm
 = 50, 0.2 ≤ dp < 1 mm and a stirrer speed of 1400 rpm
		Table S2 in the ESI† summarises the experimental and modeling results of all kinetic experiments. The experimental results include the approximate reaction time until reaching the reaction equilibrium teq, the determined equilibrium conversion XeqPET, the equilibrium composition of the reaction mixture Kx, the BHET process yield at the reaction equilibrium Yprocess,eqBHET and the monomer mole fraction xproductMonomer in the BHET product. Regarding the modeling, Table S2† contains the initial composition of the reaction mixture as well as the retrieved values for the kinetic constants k. In addition to the kinetic investigations in Table S2,† individual experiments were conducted to gain a detailed understanding of the kinetic enhancement provided by GVL.
|  | ||
| Fig. 2  (a) PET conversion XPET over reaction time t for different GVL ratios. Experimental data (symbols) vs. modeling (lines) using eqn (1), (3) and (4); (b) kinetic constants k obtained from the pseudo-first order reversible reaction model (eqn (3)) applied to the kinetic experiments with respect to the reference kinetic series kref with  of 0 (grey = reference), 0.118 (red), 0.236 (blue) and 0.473 (green). Other reaction conditions are: T = 190 °C, p = 1 bar,  = 6,  = 50, 0.2 ≤ dp < 1 mm and a stirrer speed of 1400 rpm. A detailed summary of all results can be found in Table S2.† | ||
It becomes evident from Fig. 2 (a) that the reaction mixture without GVL required approximately 21 minutes to achieve an equilibrium conversion of around 93%. In comparison to other PET glycolysis reactions catalysed by ZnAc2, this reaction kinetics is already fast. This is attributed to the small-sized PET particles (0.2 ≤ dp < 1 mm) used in our approximately 3.5 ml reaction mixture. Notably, the same equilibrium conversion was reached much faster by using GVL co-solvent without influencing the reaction equilibrium, i.e. still allowing nearly full conversion. Further, boiling of the liquid phase during the experiments was not observed probably due to the high boiling temperature of GVL (205 °C at 1 bar (ref. 35)) and the absence of any influencing temperature minimum azeotrope between EG and GVL under the investigated conditions. Besides, Fig. 2 (a) illustrates the results of the kinetic modeling, wherein the kinetic constant k was retrieved from the experimental data using Kx from experimental data obtained at the reaction equilibrium (see Table S2 in the ESI†). For a more detailed understanding of the kinetic model, Fig. S3 in the ESI† displays the changing mole fractions of the components during depolymerization over reaction time for experimental series 4  with the highest amount of GVL, where EG was present still in excess. Consequently, the concentration of EG barely changes along the progressing reaction time (according to reviews from paper 3) (−2.12%). That is why the reaction can be categorized as a pseudo-first order reversible reaction even at highest studied GVL concentrations.
 with the highest amount of GVL, where EG was present still in excess. Consequently, the concentration of EG barely changes along the progressing reaction time (according to reviews from paper 3) (−2.12%). That is why the reaction can be categorized as a pseudo-first order reversible reaction even at highest studied GVL concentrations.
Additionally, Fig. 2 (b) provides an overview of the obtained kinetic constants in relation to the reference kinetic series. The kinetic constant increases linearly with the higher GVL/EG ratio in the mixture. A ratio  of 0.118 already doubles the kinetic constant, and a ratio of 0.473 increases the value by nearly five times compared to the reference k. These values underline the accelerating effect of GVL on the reaction kinetics at isothermal and isobaric conditions (190 °C, 1 bar).
 of 0.118 already doubles the kinetic constant, and a ratio of 0.473 increases the value by nearly five times compared to the reference k. These values underline the accelerating effect of GVL on the reaction kinetics at isothermal and isobaric conditions (190 °C, 1 bar).
|  | ||
| Fig. 3  (a) PET conversion XPET over reaction time t for different pressures and temperatures. Experimental data (symbols) vs. modeling (lines) using eqn (1), (3) and (4); (b) kinetic constants k obtained from the pseudo-first order reversible reaction model (eqn (3)) applied to the kinetic experiments with respect to the reference kinetic series kref with T = 190 °C & p = 1 bar (dark grey = reference), T = 215 °C & p = 2 bar (mid-grey) and T = 232 °C & p = 3 bar (light grey). Other reaction conditions are:  = 0,  = 6,  = 50, 0.2 ≤ dp < 1 mm and a stirrer speed of 1400 rpm. A detailed summary of all results can be found in Table S2.† | ||
As expected, the pressure-assisted temperature increase accelerated the reaction, whereby the elevated reaction temperatures were reached after the dead time of 3 minutes, as can be seen in Fig. S2 in the ESI.† The temperature increase had also a slight effect on the reaction equilibrium of the endothermic reaction because the mean equilibrium conversion increased from 93 over 93.5 to 96%. Additionally, Fig. 3 (a) shows the results of the kinetic modeling using eqn (1), (3) and (4), wherein the kinetic constant k was retrieved from the experimental data using Kx from experimental data obtained at the reaction equilibrium (see Table S2 in the ESI†). The modeled and experimental data exhibit good agreement, reinforcing the characterization of the reaction as a pseudo-first order reversible reaction, as described in the previous section.
Fig. 3 (b) also provides an overview of all kinetic constants in relation to the reference kinetic series. A 25-fold increase with respect to the reference kinetic series for a temperature of 232 °C at a pressure of 3 bar underlines the substantial temperature effect on the reaction kinetics as known from the Arrhenius approach.
 of 0.118 was selected, and a temperature of 232 °C at 3 bar was applied as can be seen in Table 2, reaction series no. 7. Fig. 4 (a) presents the results of the PET conversion over time for this kinetic series, along with the reference kinetic series 1 and the optimal outcomes derived from the individual effects of GVL (series 4) and temperature (series 6) on the kinetics according to Table 2. Fig. 4 (b) provides an overview of the specific kinetic constants in relation to the reference kinetic series 1.
 of 0.118 was selected, and a temperature of 232 °C at 3 bar was applied as can be seen in Table 2, reaction series no. 7. Fig. 4 (a) presents the results of the PET conversion over time for this kinetic series, along with the reference kinetic series 1 and the optimal outcomes derived from the individual effects of GVL (series 4) and temperature (series 6) on the kinetics according to Table 2. Fig. 4 (b) provides an overview of the specific kinetic constants in relation to the reference kinetic series 1.
        |  | ||
| Fig. 4  (a) PET conversion XPET over reaction time t for different temperatures, pressures and GVL ratios. Experimental data (symbols) vs. modeling (lines) using eqn (1), (3) and (4); (b) kinetic constants k obtained from the pseudo-first order reversible reaction model (eqn (3)) applied to the kinetic experiments with respect to the reference kinetic series kref with T = 190 °C & p = 1 bar &  = 0 (dark grey = reference), T = 190 °C & p = 1 bar &  = 0.473 (green), T = 232 °C & p = 3 bar &  = 0 (light grey) and T = 232 °C & p = 3 bar &  = 0.118 (orange). Other reaction conditions are:  = 6,  = 50, 0.2 ≤ dp < 1 mm and a stirrer speed of 1400 rpm. A detailed summary of all results can be found in Table S2.† | ||
Evidently, from Fig. 4, it can be concluded that introducing a small quantity of GVL into the mixture, and simultaneously elevating the reaction conditions to 232 °C at 3 bar, whereby the reaction temperature was reached after the dead time of 3 minutes, (see Fig. S2 in the ESI†) allowed for a 70-fold increase in reaction kinetics over the reference conditions. Just after 15 seconds of reaction time, the first sample of kinetic series 7 reached a near-equilibrium state. Because experimental data was not available between the reaction times of 0 and 15 seconds, fitting the kinetic constant k using Kx from experimental data obtained at the reaction equilibrium (see Table S2 in the ESI†) might not be that accurate. Nonetheless, the primary objective of the model in this study is to quantify the acceleration of the reaction kinetics through the kinetic constant. However, all obtained experimental data points align closely with the respective modeled kinetic curve in Fig. 4 (a). The significant amplification of almost two orders of magnitude relative to the reference kinetic series shown in Fig. 4 (b) underlines the substantial impact of temperature increase on the GVL-assisted PET glycolysis reaction.
In order to evaluate the overall PET glycolysis process, the BHET process yield (see eqn (2)) at the reaction equilibrium was determined for all kinetic experimental series (see Table 2) conducted in this work. Values of about 80% were determined for the reference series no. 1 (190 °C, 1 bar), as well as for series no. 5 and 6 due to identical liquid composition of EG and H2O in the downstream process, as anticipated. Comparable literature processes also yielded values ranging between 75 and 82%.20,27,28 Apart from the beneficial effect of GVL on the kinetics, the overall process yield obtained from maximum conversion samples decreased from about 80% without GVL to 55% with the highest amount of GVL investigated in this study. This decline is attributed to the increasing fraction of GVL in the liquid phase of the downstream process and the enhanced solubility of BHET in GVL.33
|  | ||
| Fig. 5  PET conversion XPET over reaction time t for different experimental procedures. Conventional procedure of putting complete prepared samples in the heating block and setting the dead time due to heating-up to 3 minutes (cyan); pre-heating the solvents and ZnAc2 at the reaction temperature prior to PET dosage and using no dead time (yellow); pre-dissolving PET in GVL36 prior to adding pre-heated ZnAc2/EG catalyst solution and using no dead time (purple). Reaction conditions were  = 0.946,  = 6, T = 190 °C, p = 1 bar,  = 50, PET flakes = 5 × 5 × 0.2 mm and a stirrer speed of 1400 rpm adapted from Chen et al.36 and modified to enable complete dissolution of PET in GVL. | ||
It can be observed from Fig. 5 that there were no significant differences in the results concerning XPET among the three procedures. Please note that all samples had the same composition of reactants and solvents as well as the same reaction time. This is a surprising effect contrary to the expectations and experience from the literature.27 These experiments disproved the first hypothesis. Continuing with the remaining hypothesis, we know from other works (“Activity-Based Models to Predict Kinetics of Levulinic Acid Esterification”37) that co-solvents manipulate the thermodynamic activity of the reactants and the catalysts. Hence, we postulate based on the results from Fig. 5 that GVL has a beneficial effect on the activity of ZnAc2, EG or PET. Following hypothesis three, we conducted experiments without the ZnAc2 catalyst with and without GVL in the reaction system. To ensure a measurable conversion within a reasonable timeframe, samples were measured at 232 °C and 3 bar as described in sections 2.2 and 2.3. The conversion of the samples was compared after 7 minutes of reaction time and the results are depicted in Fig. 6.
Under the same reaction conditions, a 20% conversion was achieved without GVL, whereas only a 13% conversion was attained for the samples with GVL. The results obtained without the catalyst suggest that the enhanced kinetics due to the presence of GVL in the reaction mixture is not attributed to the interaction of GVL with the reactants. Instead, they are more likely attributed to the interaction between GVL and the ZnAc2 catalyst, possibly involving acid–base synergistic catalysis.34 Detailed descriptions of reaction mechanisms involving systems containing ZnAc2/Lewis base catalysts in PET glycolysis have already been published in the literature.28,34,38,39 Applied to our Lewis acid (ZnAc2) catalysed reaction system, it is plausible that GVL functions as a Lewis base. This is because the oxygen atom in the lactone ring has a lone pair of electrons, which can be donated to form a bond with the electron-deficient species (Lewis acid). Consequently, GVL can coordinate with ZnAc2 to promote bond cleavage and catalytic activity.
 of 0.473 accelerated the reaction by a factor of about 5 with respect to the reference kinetic series. Furthermore, it was concluded that the enhanced kinetics due to the presence of GVL in the reaction mixture is attributed to the interaction between ZnAc2 and GVL possibly due to acid–base synergistic catalysis. Apart from that, the influence of slightly elevated reaction conditions in terms of temperature and pressure was explored. The temperature increase assisted by pressure elevation markedly accelerated the reaction, with 5 times the acceleration at 215 °C and 25 times the acceleration at 232 °C, underscoring the substantial temperature effect on the reaction kinetics. Lastly, the combined temperature increase on the GVL-assisted PET glycolysis reaction was assessed. At a temperature of 232 °C, pressure of 3 bar and a relatively small GVL to EG ratio
 of 0.473 accelerated the reaction by a factor of about 5 with respect to the reference kinetic series. Furthermore, it was concluded that the enhanced kinetics due to the presence of GVL in the reaction mixture is attributed to the interaction between ZnAc2 and GVL possibly due to acid–base synergistic catalysis. Apart from that, the influence of slightly elevated reaction conditions in terms of temperature and pressure was explored. The temperature increase assisted by pressure elevation markedly accelerated the reaction, with 5 times the acceleration at 215 °C and 25 times the acceleration at 232 °C, underscoring the substantial temperature effect on the reaction kinetics. Lastly, the combined temperature increase on the GVL-assisted PET glycolysis reaction was assessed. At a temperature of 232 °C, pressure of 3 bar and a relatively small GVL to EG ratio  of 0.118, a 70-fold amplification relative to the reference kinetics was achieved. This emphasizes the positive combined effect on the reaction kinetics of all the influencing factors investigated in this study. All investigated approaches significantly decreased the reaction time of the PET glycolysis while staying under comparatively mild reaction conditions, thus contributing to a greener PET glycolysis reaction. Nevertheless, further investigations are required regarding the downstream of the PET glycolysis process. While only elevating the temperature and pressure does not affect the downstream process, the introduction of GVL presents a trade-off between accelerating reaction kinetics and reducing process yield. Finding an optimal balance becomes crucial. Additionally, recycling of ZnAc2, EG and GVL should be implemented to minimize the need for downstream purification and to reduce costs. Further, repolymerisation of the retrieved BHET product, which contains both monomers and dimers, into PET is necessary. This step will help evaluate whether additional purification steps to separate monomers and dimers are required.
 of 0.118, a 70-fold amplification relative to the reference kinetics was achieved. This emphasizes the positive combined effect on the reaction kinetics of all the influencing factors investigated in this study. All investigated approaches significantly decreased the reaction time of the PET glycolysis while staying under comparatively mild reaction conditions, thus contributing to a greener PET glycolysis reaction. Nevertheless, further investigations are required regarding the downstream of the PET glycolysis process. While only elevating the temperature and pressure does not affect the downstream process, the introduction of GVL presents a trade-off between accelerating reaction kinetics and reducing process yield. Finding an optimal balance becomes crucial. Additionally, recycling of ZnAc2, EG and GVL should be implemented to minimize the need for downstream purification and to reduce costs. Further, repolymerisation of the retrieved BHET product, which contains both monomers and dimers, into PET is necessary. This step will help evaluate whether additional purification steps to separate monomers and dimers are required.
    
    
      | Footnote | 
| † Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4re00235k | 
| This journal is © The Royal Society of Chemistry 2024 |