Catalytic Ring-Opening Polymerisation of Cyclic Ethylene Carbonate: Importance of Elementary Steps for Determining Polymer Properties Revealed via DFT-MTD Simulations Validated Using Kinetic Measurements

The production of CO2-containing polymers is still very demanding in terms of controlling the synthesis of products with pre-defined CO2 content and molecular weight. An elegant way of synthesising these polymers is via CO2-containing building blocks, such as cyclic ethylene carbonate (cEC), via catalytic ring-opening polymerisation. However, to date, the mechanism of this reaction and control parameters have not been elucidated. In this work, using DFT-metadynamics simulations for exploiting the potential of the polymerisation process, we aim to shed more light on the mechanisms of the interaction between catalysts (in particular, the catalysts K3VO4, K3PO4, and Na2SnO3) and the cEC monomer in the propagation step of the polymeric chain and the occurring CO2 release. Confirming the simulation results via subsequent kinetics measurements indicates that, depending on the catalyst’s characteristics, it can be attached reversibly to the polymeric chain during polymerisation, resulting in a defined lifetime of the activated polymer chain. The second anionic oxygen of the catalyst can promote the catalyst’s transfer to another electrophilic cEC monomer, terminating the growth of the first chain and initiating the propagation of the new polymer chain. This transfer reaction is an essential step in controlling the molecular weight of the products.


Introduction
Increasing concern about climate change, connected to the use of (diminishing) fossil resources, necessitates the development of new, sustainable technologies from renewable sources for the generation of chemicals, energy, and materials [1][2][3][4][5].To this end, cyclic ethylene carbonate (cEC in short), due to ecological and economic reasons, and its ready availability from ethylene oxide and carbon dioxide, seems to be a very attractive monomer and reagent for polyether carbonate production [6,7].The obtained aliphatic polycarbonates are of interest due to their new and valuable applications in biodegradable and biocompatible materials.Hence, the presence of carbonate functional groups in the polymer chain facilitates the depolymerisation of aliphatic polycarbonates.Interestingly, in ring-opening polymerisation (ROP) of cEC, the polymers formed may have lower densities than the monomers (volume expansion may accompany polymerisation), which can be beneficial for industrial applications [8][9][10][11].The catalysed ROP of cEC is a low-pressure technology, representing a technical advantage compared to the direct copolymerisation of ethylene oxide (EO) and CO 2 [12][13][14][15][16][17][18][19].Both thermodynamic and kinetic factors are crucial for the ROP of cEC.It is worth noting that the maximum temperature for the homopolymerisation Polymers 2024, 16, 136 2 of 24 of CEC is below 25 • C, but polymerisations have also occurred at over 100 • C [17].At high polymerisation temperatures, the resultant polymers' repeat units are a mixture of monomeric units (carbonate units) and the corresponding oxide units, meaning that CO 2 release occurs during the polymerisation.The release of the CO 2 makes the overall reaction ∆S positive, so this polymerisation becomes thermodynamically possible only at high temperatures.In cEC polymerisation, the purpose is not to completely avoid CO 2 release, but keeping CO 2 to some extent in the final polymer product is advantageous [17].To overcome the kinetic stability of cEC and polymerise the cEC, modern catalysts and catalytic processes have been developed in the last few decades to achieve the ROP quickly.Notably, an efficient catalyst must conduct the reaction in such a way that CO 2 release is minimised during the polymerisation of cEC.A series of recently investigated catalysts approved for cEC polymerisation is of particular industrial interest, i.e., they are structurally simple, commercially available, and efficient catalysts.The relatively high conversion values and CO 2 content percentages that were experimentally obtained and approved by the Covestro company [4,5,20] in the ROP of cEC confirm that these catalysts are exciting candidates with the potential for obtaining a more fundamental understanding with kinetic studies and computational methods.Among the examined series of 30 different catalysts, we selected K 3 VO 4 , K 3 PO 4 , and Na 2 SnO 3 as representative examples for the current study.The cyclic carbonate conversion reactions cover carbonate formation and ether linkages in the final product, along with some CO 2 release (Scheme 1).
technology, representing a technical advantage compared to the direct copolymerisation of ethylene oxide (EO) and CO2 [12][13][14][15][16][17][18][19].Both thermodynamic and kinetic factors are crucial for the ROP of cEC.It is worth noting that the maximum temperature for the homopolymerisation of CEC is below 25 °C, but polymerisations have also occurred at over 100 °C [17].At high polymerisation temperatures, the resultant polymers' repeat units are a mixture of monomeric units (carbonate units) and the corresponding oxide units, meaning that CO2 release occurs during the polymerisation.The release of the CO2 makes the overall reaction ΔS positive, so this polymerisation becomes thermodynamically possible only at high temperatures.In cEC polymerisation, the purpose is not to completely avoid CO2 release, but keeping CO2 to some extent in the final polymer product is advantageous [17].To overcome the kinetic stability of cEC and polymerise the cEC, modern catalysts and catalytic processes have been developed in the last few decades to achieve the ROP quickly.Notably, an efficient catalyst must conduct the reaction in such a way that CO2 release is minimised during the polymerisation of cEC.A series of recently investigated catalysts approved for cEC polymerisation is of particular industrial interest, i.e., they are structurally simple, commercially available, and efficient catalysts.The relatively high conversion values and CO2 content percentages that were experimentally obtained and approved by the Covestro company [4,5,20] in the ROP of cEC confirm that these catalysts are exciting candidates with the potential for obtaining a more fundamental understanding with kinetic studies and computational methods.Among the examined series of 30 different catalysts, we selected K3VO4, K3PO4, and Na2SnO3 as representative examples for the current study.The cyclic carbonate conversion reactions cover carbonate formation and ether linkages in the final product, along with some CO2 release (Scheme 1).Scheme 1. Catalytic ROP of cEC and exemplary molecular structure of the respective polyether carbonate polyol formed.
For the selected catalysts, the cations are either equal (K3PO4 and K3VO4) or, if they are different, no significant effect from the cationic counterpart is observed in simulations and experiments (Na2SnO3 vs. K2SnO3).Hence, the main difference in catalytic performance for the selected catalysts is caused by the anionic part, and we have considered the anionic pathway (Scheme 2) for the current study.Anionic ROP (Scheme 2) begins with a nucleophilic attack by an anionic initiator, such as the anionic counterpart of the catalyst, or can be generated, e.g., via proton abstraction from starter alcohol.The nucleophile can either attack the carbonyl carbon or the methylene carbon of the cyclic carbonate [21][22][23].
Scheme 1. Catalytic ROP of cEC and exemplary molecular structure of the respective polyether carbonate polyol formed.
For the selected catalysts, the cations are either equal (K 3 PO 4 and K 3 VO 4 ) or, if they are different, no significant effect from the cationic counterpart is observed in simulations and experiments (Na 2 SnO 3 vs.K 2 SnO 3 ).Hence, the main difference in catalytic performance for the selected catalysts is caused by the anionic part, and we have considered the anionic pathway (Scheme 2) for the current study.Anionic ROP (Scheme 2) begins with a nucleophilic attack by an anionic initiator, such as the anionic counterpart of the catalyst, or can be generated, e.g., via proton abstraction from starter alcohol.The nucleophile can either attack the carbonyl carbon or the methylene carbon of the cyclic carbonate [21][22][23].
Polymers 2024, 16, x FOR PEER REVIEW 2 of 27 technology, representing a technical advantage compared to the direct copolymerisation of ethylene oxide (EO) and CO2 [12][13][14][15][16][17][18][19].Both thermodynamic and kinetic factors are crucial for the ROP of cEC.It is worth noting that the maximum temperature for the homopolymerisation of CEC is below 25 °C, but polymerisations have also occurred at over 100 °C [17].At high polymerisation temperatures, the resultant polymers' repeat units are a mixture of monomeric units (carbonate units) and the corresponding oxide units, meaning that CO2 release occurs during the polymerisation.The release of the CO2 makes the overall reaction ΔS positive, so this polymerisation becomes thermodynamically possible only at high temperatures.In cEC polymerisation, the purpose is not to completely avoid CO2 release, but keeping CO2 to some extent in the final polymer product is advantageous [17].
To overcome the kinetic stability of cEC and polymerise the cEC, modern catalysts and catalytic processes have been developed in the last few decades to achieve the ROP quickly.Notably, an efficient catalyst must conduct the reaction in such a way that CO2 release is minimised during the polymerisation of cEC.A series of recently investigated catalysts approved for cEC polymerisation is of particular industrial interest, i.e., they are structurally simple, commercially available, and efficient catalysts.The relatively high conversion values and CO2 content percentages that were experimentally obtained and approved by the Covestro company [4,5,20] in the ROP of cEC confirm that these catalysts are exciting candidates with the potential for obtaining a more fundamental understanding with kinetic studies and computational methods.Among the examined series of 30 different catalysts, we selected K3VO4, K3PO4, and Na2SnO3 as representative examples for the current study.The cyclic carbonate conversion reactions cover carbonate formation and ether linkages in the final product, along with some CO2 release (Scheme 1).
Scheme 1. Catalytic ROP of cEC and exemplary molecular structure of the respective polyether carbonate polyol formed.
For the selected catalysts, the cations are either equal (K3PO4 and K3VO4) or, if they are different, no significant effect from the cationic counterpart is observed in simulations and experiments (Na2SnO3 vs. K2SnO3).Hence, the main difference in catalytic performance for the selected catalysts is caused by the anionic part, and we have considered the anionic pathway (Scheme 2) for the current study.Anionic ROP (Scheme 2) begins with a nucleophilic attack by an anionic initiator, such as the anionic counterpart of the catalyst, or can be generated, e.g., via proton abstraction from starter alcohol.The nucleophile can either attack the carbonyl carbon or the methylene carbon of the cyclic carbonate [21][22][23].

Scheme 2. The anionic ROP of cyclic carbonates (considered in this work).
Polymers 2024, 16, 136 3 of 24 In this work, we aim to understand the details of the interaction between the catalysts and cEC monomer, employing DFT-MTD simulations and kinetic measurements to shed more light on the elementary steps for the CO 2 release mechanism and differences between various considered catalysts.The main questions that we address are the consumption of the entire starter molecule in the initiation step, which pathways are dominant in the propagation step throughout the interaction between catalyst and monomer, and the possible effects on the CO 2 content and molecular weight of the polymer are.
Recent findings obtained using advanced ab initio molecular dynamics (AIMD) simulations showed that the possibility of identification of alternative pathways of various complex reactions could be challenged in a dynamic picture [24][25][26][27][28][29].This motivated us to investigate the ROP of cEC in the presence of the catalysts, as mentioned above, using DFT-metadynamics simulations as the computational method of choice.We employ metadynamics simulations to calculate the free energy surfaces (FESs) for the ROP and CO 2 release mechanisms at a defined temperature.Simulation of a free energy landscape rather than a zero Kelvin potential energy surface allows the incorporation of entropic contributions due to the flexibility of motion of free molecules in the reaction environment.According to previous studies, the entropic contributions affect the formation of molecular complexes from small molecules and their respective transition state barriers of formation [30][31][32][33].
Understanding the underlying mechanisms and identifying the critical steps of catalysed ring opening and chain growth of cyclic carbonates and CO 2 release pathways is the basis for computational chemistry-assisted research, leading to improved catalyst structures and performance.Using the results and insights obtained in this work, novel catalyst candidates can be rationally designed for better performance and process optimisation.
To make this possible, comparing the computer-calculated values with experimental ones is an elegant approach.By determining the kinetics, it is possible to understand the reaction better.Understanding the kinetics also opens up the control possibilities of ring-opening polymerisation.The mechanism of the ROP is not fully clarified, so it is still unclear why the basicity has an influence and why specific catalysts no longer react during the reaction [6].These questions are, therefore, the central focus of this paper, and with them, the clarification of the mechanism of ring-opening polymerisation.
These polymers give access to a wide-ranging product portfolio.For example, using a hydrophobic initiator, functional non-ionic surfactants can be produced that are both environmentally friendly due to their excellent biodegradability and can be recycled.For example, these surfactants have already been proven to purify the microplastic contamination of water [34].The production of low-molecular polymers (700-1200 g/mol) for the production of surfactants was the aim of the funded project on which this work is based.However, the findings of this work can also be applied directly to higher degrees of polymerisation.

Computational Details
All DFT calculations, including geometry optimisations and ab initio molecular dynamics (AIMD) simulations, were performed using the CP2K program [35] with the Gaussian and plane-wave (GPW) method.The valence orbitals were expanded in the DZVP-MOLOPT Gaussian basis set in combination with Goedecker, Teter, and Hutter pseudopotentials and were used with a plane wave cutoff energy of 280 Ry.We used the PBE density functional [36] augmented with Grimme D3 dispersion correction [37].The criterion for self-consistent field convergence was set to 5.0 × 10 −7 .The AIMD simulations were done in the NVT ensemble, with the temperature controlled by a CSVR (canonical sampling through velocity rescaling) thermostat, set at various temperatures (413, 423, 433, and 443 K), and a period of 500 fs.The MD time step was 0.5 fs.To investigate the reaction mechanism at finite temperature, characterise the reaction pathway, and identify the transition state region between reactant and product states, we performed metady-namics simulations using the PLUMED 2.8 plugin in combination with CP2K [38][39][40].The simulations were initiated from the optimised molecular structures and conducted until several transition state re-crossing events could be observed.To prevent sampling from unnecessary regions of the FES, harmonic walls with force constant K = 250 kcal/molÅ 2 were used for CVs.The Gaussian bias potentials were added every 100 steps.Multiple test simulations were run to set the computational parameters, including the number and type of CVs, the height and width of Gaussian bias potentials, and quadratic walls for each part of the reaction mechanism.
We performed our simulations on a molecular cluster schematically shown in Scheme 3. As depicted in Scheme 3, the molecular cluster consists of six cEC monomers, one alcohol molecule (F = 1 for simplification), and a catalyst molecule.The typical catalyst molecule includes the anionic (MO 4 3− or MO 3 2− ) and cationic (A + ) counterions.
Polymers 2024, 16, x FOR PEER REVIEW 4 of 27 reaction mechanism at finite temperature, characterise the reaction pathway, and identify the transition state region between reactant and product states, we performed metadynamics simulations using the PLUMED 2.8 plugin in combination with CP2K [38][39][40].The simulations were initiated from the optimised molecular structures and conducted until several transition state re-crossing events could be observed.To prevent sampling from unnecessary regions of the FES, harmonic walls with force constant K = 250 kcal/molÅ 2 were used for CVs.The Gaussian bias potentials were added every 100 steps.Multiple test simulations were run to set the computational parameters, including the number and type of CVs, the height and width of Gaussian bias potentials, and quadratic walls for each part of the reaction mechanism.We performed our simulations on a molecular cluster schematically shown in Scheme 3. As depicted in Scheme 3, the molecular cluster consists of six cEC monomers, one alcohol molecule (F = 1 for simplification), and a catalyst molecule.The typical catalyst molecule includes the anionic (MO4 3− or MO3 2− ) and cationic (A + ) counterions.We considered six cEC monomers in our simulations since six cECs can make the first shell interact with the catalyst with a reasonable computational cost for AIMD calculations.In fact, we investigated the interaction between the catalyst anion and cEC monomer, and a significant number of molecules participating in reaction mechanisms is one of each species.To estimate the interactions between cEC monomers, ROH starter, and catalyst-anion, we have probed the variation of average H-bond distances between ROH and six surrounding cEC monomers (ROH…O(cEC)) in the first solvation shell through an unbiased equilibrium MD simulation in Figure S1A in SI.The ROH-cEC distances are compared with the ROH-O(anion-CAT).The main interactions involved are those between ROH and the anion (shorter distances) that result in proton transfer in the initiation step of polymerisation.The interaction between ROH and O-cEC (according to distances) is weaker and O-cEC mainly interacts with the surrounding cations.Hence, ROH does not exist that long in the reaction medium and is deprotonated by the anion of the catalyst.After the proton abstraction, the main H bond is between the protonated catalyst anion and the surrounding molecules (cEC monomers and alkoxide), which can cover the protonated anion sufficiently.The reaction temperature in the experiment is above 150 °C, and the catalysts are expected to be solubilised entirely under these conditions.Additionally, lab results indicate that the catalyst cannot be separated adequately via filtration after the reaction.Both facts conclude that considering homogenous reaction conditions is a valid starting point.The collective variables (CVs) were the relevant distances in each reaction step.We used VMD 1.9.4 software to visualise the MD trajectories and snapshots [41].We considered six cEC monomers in our simulations since six cECs can make the first shell interact with the catalyst with a reasonable computational cost for AIMD calculations.In fact, we investigated the interaction between the catalyst anion and cEC monomer, and a significant number of molecules participating in reaction mechanisms is one of each species.To estimate the interactions between cEC monomers, ROH starter, and catalyst-anion, we have probed the variation of average H-bond distances between ROH and six surrounding cEC monomers (ROH• • • O(cEC)) in the first solvation shell through an unbiased equilibrium MD simulation in Figure S1A in SI.The ROH-cEC distances are compared with the ROH-O(anion-CAT).The main interactions involved are those between ROH and the anion (shorter distances) that result in proton transfer in the initiation step of polymerisation.The interaction between ROH and O-cEC (according to distances) is weaker and O-cEC mainly interacts with the surrounding cations.Hence, ROH does not exist that long in the reaction medium and is deprotonated by the anion of the catalyst.After the proton abstraction, the main H bond is between the protonated catalyst anion and the surrounding molecules (cEC monomers and alkoxide), which can cover the protonated anion sufficiently.The reaction temperature in the experiment is above 150 • C, and the catalysts are expected to be solubilised entirely under these conditions.Additionally, lab results indicate that the catalyst cannot be separated adequately via filtration after the reaction.Both facts conclude that considering homogenous reaction conditions is a valid starting point.The collective variables (CVs) were the relevant distances in each reaction step.We used VMD 1.9.4 software to visualise the MD trajectories and snapshots [41].

Experimental Method
All reactions proceeded at the same conditions using a 50 mL reactor with a heatable jacket.The reactor has a GL25 opening for measurement with an in situ IR probe and two feed valves for operation under inert gas.The reactor was heated with a Haake F6 thermostat.Silicone oil AP200 was chosen as the operating fluid for the thermostat.The IKA RCT Basic magnetic stirrer was used for mixing.A Mettler Toledo Reakt IR 15 was used for the in situ measurements.The spectrometer has a Si probe connected to the spectrometer via AgX 9.5 mm × 1.5 mm fibre optics.The measured wavenumber range is from 4000 cm −1 to 800 cm −1 with a resolution of 4 cm −1 .The reactor was brought to experimental temperature according to Tables S2 and S6 (see Supporting Information), and the in situ IR probe was supplied with fresh nitrogen (l) for cooling the detector.The blank value of the IR probe was measured with the empty reactor, in which only the magnetic stirring bar was present.After measuring the blank spectrum, cyclic ethylene carbonate (20 g, 1 eq.), which had been preheated in a 60 • C water bath, was added to the reactor, the stirrer was set to 600 RPM, and the in situ measurement was started At least 10 min waiting time was used to obtain a baseline measurement for cEC.Ethylene glycol (10 mol%, 0.1 eq.) and catalyst (1 mol%, 0.01 eq.) were weighed beforehand and added quickly to the reactor.As the reaction progressed, the changes in the characteristic bands were monitored in the Software IC IR 7 (see Figure 1): 1800-1760 cm −1 [42] for cEC (C-O-C plugin vibrations, strained system) [42], 1750-1700 cm −1 [42] for CO 2 (linear unstrained system) [42] and the ether band at 1260 cm −1 [42], (these are the expected signals for a nascent poly(ethylene ether carbonate)) [42][43][44][45].The reaction was allowed to proceed until the band typical for cEC at 1800-1750 cm −1 was diminished entirely.

Experimental Method
All reactions proceeded at the same conditions using a 50 mL reactor with a heatable jacket.The reactor has a GL25 opening for measurement with an in situ IR probe and two feed valves for operation under inert gas.The reactor was heated with a Haake F6 thermostat.Silicone oil AP200 was chosen as the operating fluid for the thermostat.The IKA RCT Basic magnetic stirrer was used for mixing.A Mettler Toledo Reakt IR 15 was used for the in situ measurements.The spectrometer has a Si probe connected to the spectrometer via AgX 9.5 mm × 1.5 mm fibre optics.The measured wavenumber range is from 4000 cm -1 to 800 cm -1 with a resolution of 4 cm -1 .The reactor was brought to experimental temperature according to Tables S2 and S6 (see Supporting Information), and the in situ IR probe was supplied with fresh nitrogen (l) for cooling the detector.The blank value of the IR probe was measured with the empty reactor, in which only the magnetic stirring bar was present.After measuring the blank spectrum, cyclic ethylene carbonate (20 g, 1 eq.), which had been preheated in a 60 °C water bath, was added to the reactor, the stirrer was set to 600 RPM, and the in situ measurement was started At least 10 min waiting time was used to obtain a baseline measurement for cEC.Ethylene glycol (10 mol%, 0.1 eq.) and catalyst (1 mol%, 0.01 eq.) were weighed beforehand and added quickly to the reactor.As the reaction progressed, the changes in the characteristic bands were monitored in the Software IC IR 7 (see Figure 1): 1800-1760 cm -1 [42] for cEC (C-O-C plugin vibrations, strained system) [42], 1750-1700 cm -1 [42] for CO2 (linear unstrained system) [42] and the ether band at 1260 cm -1 [42], (these are the expected signals for a nascent poly(ethylene ether carbonate)) [42][43][44][45][46].The reaction was allowed to proceed until the band typical for cEC at 1800-1750 cm -1 was diminished entirely.

Results and Discussion
Our previous study considered the ring-opening polymerisation of cEC in the presence of the starter ROH, which begins with the proton abstraction from the starter ROH molecule by the catalyst anion (step I [46]).The generated RO − attacks the C=O group of the first neighbouring cEC.Then, the second generated alkoxide attacks the next cEC (step II).This mechanism is shown in Figure 2, considering our model cluster.The rate-determining step for this mechanism is the cEC ring opening (step II).The values of the ∆G ̸ = for the two steps are reported in Table 1.molecule by the catalyst anion (step I [47]).The generated RO attacks the C=O group of the first neighbouring cEC.Then, the second generated alkoxide attacks the next cEC (step II).This mechanism is shown in Figure 2, considering our model cluster.The rate-determining step for this mechanism is the cEC ring opening (step II).The values of the ΔG ≠ for the two steps are reported in Table 1.Step II, the ring opening of the second cEC via the second nucleophilic attack by alkoxide and ring opening.
Table 1.The calculated ΔG ≠ for steps I and II (depicted in Figure 1).

Catalyst
ΔG ≠ [kcal/mol] I/II Na2SnO3 5.8/11.2K3VO4 5.5/7.7 K3PO4 4.5/6.9 As shown in Table 1, the barriers are low in the presence of ROH (initiation).The barriers, shown in Table 1 for step II, are correlated to the second cEC opening through a nucleophilic attack by the first cleaved cEC.This step is the rate-determining step since the proton abstraction (step I) has the lower barrier.However, considering that the concentration of the starter molecule is less than the cEC monomer at a particular stage in the reaction, the starter is consumed, and only cEC and catalyst are present in the reaction environment.Therefore, the reaction proceeds throughout the catalyst-cEC interaction (propagation), and investigating this pathway is the primary purpose of the current study.
In this work, we have considered Na2SnO3, K3VO4, and K3PO4 as the catalysts for the cEC ring-opening polymerisation.From a previous study, we found that the M-O (central atom-oxygen) bond length has a crucial impact on the catalyst basicity, i.e., a longer M-O bond leads to more basicity (easier proton abstraction) but less nucleophilicity for a nucleophilic attack to the CH2 of the cEC ring.By the basicity, we refer to the basicity of the Step II, the ring opening of the second cEC via the second nucleophilic attack by alkoxide and ring opening.
Table 1.The calculated ∆G ̸ = for steps I and II (depicted in Figure 1).As shown in Table 1, the barriers are low in the presence of ROH (initiation).The barriers, shown in Table 1 for step II, are correlated to the second cEC opening through a nucleophilic attack by the first cleaved cEC.This step is the rate-determining step since the proton abstraction (step I) has the lower barrier.However, considering that the concentration of the starter molecule is less than the cEC monomer at a particular stage in the reaction, the starter is consumed, and only cEC and catalyst are present in the reaction environment.Therefore, the reaction proceeds throughout the catalyst-cEC interaction (propagation), and investigating this pathway is the primary purpose of the current study.

Catalyst
In this work, we have considered Na 2 SnO 3 , K 3 VO 4 , and K 3 PO 4 as the catalysts for the cEC ring-opening polymerisation.From a previous study, we found that the M-O (central atom-oxygen) bond length has a crucial impact on the catalyst basicity, i.e., a longer M-O bond leads to more basicity (easier proton abstraction) but less nucleophilicity for a nucleophilic attack to the CH 2 of the cEC ring.By the basicity, we refer to the basicity of the transition metal/main element oxides (K 3 PO 4 , K 3 VO 4 , Na 2 SnO 3 ), inherently weaker bases than KOH and NaOH.The order of decreasing the M-O bond distance of the considered catalysts is Na 2 SnO 3 > K 3 VO 4 > K 3 PO 4 .These catalysts' average M-O bond lengths are 1.939, 1.734, and 1.571 (in Å).The average O• • • cations distances are 2.271, 2.617, and 2.604 (in Å) for Na 2 SnO 3 , K 3 VO 4 , and K 3 PO 4 , respectively.The differences in the molecular structures of the catalysts severely affect their catalytic performance.Hence, a longer M-O bond distance performs better in keeping a higher CO 2 content in the product.The planar structure of the SnO 3 2− catalyst anion causes complexation of the alkoxide anion (RO − generated from ROH) to the Sn centre for a while during the initiation step.This complexation delays the reaction in case of SnO 3 2− vs. VO 4 3− , 9.5 h vs. 5.0 h, respectively, in experiment, which is in accordance with the higher barrier of the cEC ring opening for SnO 3 2− vs. VO 4 3− (Table 2).For a detailed analysis we refer to the [46].In addition, the cation's lability also plays a role in more efficient catalytic performance.Especially for the smaller anions such as PO 4 3− , bigger cations perform better than smaller ones.The ∆G ̸ = values of the ring opening decrease from H + to K + for phosphate anion, and in the experiment, the CO 2 contribution increases with larger cations.This effect is due to the longer anion. ..ation distance and less-covalent character of the anion-cation interaction, leading to more accessibility of the catalyst anion in the reaction.For bigger anions such as SnO 3 2− and VO 4 3− , switching between Na + and K + does not affect the catalyst performance, and due to the stronger ionic character, the ∆G ̸ = values and CO 2 -contribution do not change with changing Na + to K + .The complete results of the cation influence on the ∆G ̸ = and the experimental CO 2 contribution are reported in Table 2.We note that for H 3 PO 4 and Li 3 PO 4 catalysts in Table 2 only minor conversions (due to the reduced activity) were observed.However, the conversion was higher than 80 percent for the rest of the catalysts.from ROH) to the Sn centre for a while during the initiation step.This complexation delays the reaction in case of SnO3 2− vs. VO4 3− , 9.5 h vs. 5.0 h, respectively, in experiment, which is in accordance with the higher barrier of the cEC ring opening for SnO3 2− vs. VO4 3− (Table 2).For a detailed analysis we refer to the reference 47.In addition, the cation's lability also plays a role in more efficient catalytic performance.Especially for the smaller anions such as PO4 3− , bigger cations perform better than smaller ones.The ΔG ≠ values of the ring opening decrease from H + to K + for phosphate anion, and in the experiment, the CO2 contribution increases with larger cations.This effect is due to the longer anion…ation distance and less-covalent character of the anion-cation interaction, leading to more accessibility of the catalyst anion in the reaction.For bigger anions such as SnO3 2− and VO4 3− , switching between Na + and K + does not affect the catalyst performance, and due to the stronger ionic character, the ΔG ≠ values and CO2-contribution do not change with changing Na + to K + .The complete results of the cation influence on the ΔG ≠ and the experimental CO2 contribution are reported in Table 2.We note that for H3PO4 and Li3PO4 catalysts in Table 2 only minor conversions (due to the reduced activity) were observed.However, the conversion was higher than 80 percent for the rest of the catalysts.The examination of the CO 2 release path of Figure 3 with MTD simulations confirmed that the higher nucleophilicity of the phosphate anion is the leading cause for easier CO 2 release during the polymerisation of cEC.Considering that, the mechanism shown in Figure 3, which results in the attachment of the catalyst to the cEC ring, can proceed in parallel with the pathway shown in Figure 2.This work analyses this pathway and the elementary steps to conduct this molecular complex (between the catalyst anion and the cleaved cEC) to further chain growth and polymerisation.The experimental observations show that higher catalyst concentration increases the CO 2 release.Additionally, at the beginning of polymerisation, the primary cause of the CO 2 release is due to the cEC decomposition.Hence, the hidden elementary steps in the pathway shown in Figure 3 can affect the molecular weight of the polymer as well as the CO 2 content and the overall reaction rate.We investigated this pathway in the current study to better understand the catalyst-cEC interaction and to address the role of catalyst characteristics.

Simulation of the Catalyst-cEC Attachment Path
To run metadynamics simulations on the path shown in Figure 3, we have considered the distance between the nucleophilic oxygen atom of the catalyst anion and the methylene carbon in the cEC monomer, i.e., O − • • • C(CH 2 -O), as the first collective variable (CV1) and the C(H 2 )• • • O(ethereal) inside the cEC as the second collective variable (CV2).The MTD simulations have been performed at 413, 423, 433, and 443 K.For brevity, we only show the results of 423 K; for the other temperatures, we include the data in SI (Figures S1-S3).We have shown the variation of the CV1 versus CV2 in Figure 4 for the trajectory of cEC cleavage and catalyst attachment at 423 K, in which the two reaction coordinates (two distances) are plotted versus each other (x, y axis for CV1 and CV2, respectively).Investigation of the patterns created by the CVs in each plot indicates some similarities and differences between the path of catalyst attachment for different catalysts that can lead us to different elementary steps.We note that the three catalysts' computational parameters and temperatures are equal.As shown in Figure 4, for Na 2 SnO 3 and K 3 VO 4 catalysts, the catalyst attachment path goes through an intermediate region of CV1 ≥ 2.5 Å and CV2 ≥ 2.5 Å.On the other hand, in the case of K 3 PO 4 , a diagonal transition path (that directly connects the reactant state to the catalyst-cEC complex) with a more significant density of the transient states, which is invisible for the SnO 3 2− and VO 4 3− , appears in the plot.Variation of CV1 versus CV2 in Figure 4 identifies different regions based on the density of the dots (transient structures).These regions are singled out in Figure 4 based on the values of the CVs.As can be seen in Figure 5 (left side), three stable states are observed: the reactant state (CV1 ≥ 2.5 Å, CV2 ≤ 1.6 Å), the intermediate state (CV1 ≥ 2.5 Å, CV2 ≥ 2.5 Å), and the state in which the CAT-cEC bond (complex) is formed (CV1 ≤ 1.6 Å, CV2 ≥ 2.5).The regions between these states indicate the transition state areas, i.e., TS vert., TS hor., and TS Nu−attack .On the right side of Figure 5, the schematic molecular structures corresponding to the highlighted regions are depicted.).The intermediate either returns to the closed cEC or releases CO 2 .In the case of K 3 PO 4 , an additional path can be identified: the cEC cleavage and nucleophilic attack by the catalyst anion happen simultaneously, and a diagonal transition state region is observed (TS Nu−attack ).These two distinguishable patterns can be identified in the calculated FESs according to the relative energies in kcal/mol.Figure 6 shows the calculated FESs at 423 K for the three catalysts.As seen in Figure 6, the stepwise path, which is the transition The FES of the K 3 PO 4 identifies a new concerted path in addition to the stepwise, which corresponds to the diagonal transition in that the two CVs are involved simultaneously.The numbers at the FESs show the energies in kcal/mol.The more stable intermediate, in the case of Na 2 SnO 3 and K 3 VO 4 , can conduct the reaction in both directions; that is, turning into the catalyst attachment and CO 2 -release or into the backward direction and generating the cEC, which involves an alkoxide attack and ring opening (Figure 2 Step II).In the case of K 3 PO 4 , a pattern with direct catalyst attachment (diagonal with simultaneous involvement of CVs) can be observed.).The intermediate either returns to the closed cEC or releases CO2.In the case of K3PO4, an additional path can be identified: the cEC cleavage and nucleophilic attack by the catalyst anion happen simultaneously, and a diagonal transition state region is observed (TS Nu−attack ).These two distinguishable patterns can be identified in the calculated FESs according to the relative energies in kcal/mol.Figure 6 shows the calculated FESs at 423K for the three catalysts.As seen in Figure 6, the stepwise path, which is the transition through the intermediate, can be observed in the pattern of the FES of Na2SnO3 and K3VO4.The FES of the K3PO4 identifies a new concerted path in addition to the stepwise, which corresponds to the diagonal transition in that the two CVs are involved simultaneously.The numbers at the FESs show the energies in kcal/mol.The more stable intermediate, in the case of Na2SnO3 and K3VO4, can conduct the reaction in both directions; that is, turning into the catalyst attachment and CO2-release or into the backward direction and generating the cEC, which involves an alkoxide attack and ring opening (Figure 2 Step II).In the case of K3PO4, a pattern with direct catalyst attachment (diagonal with simultaneous involvement of CVs) can be observed.In summary, the different patterns observed for stannate/vanadate vs. phosphate reflect the different natures of the chemical interactions between the catalyst and cEC (nonsimilar catalyst characteristics).The more-nucleophilic character of the catalyst anion results in the sampling of a concerted pathway of the catalyst attachment.The following section compares the calculated results with model parameters determined from kinetic measurements to address the correlations/differences between simulations and experiments.

Experimental Analysis of the Kinetics of Polymerisation of cEC, including CO2 Release
The experimental investigations on the kinetics of the polymerisation of cEC indicate that higher CO2 content can be obtained by the stannate and vanadate anions, which confirms the better performance of these catalysts for the ROP of cEC.Notably, the kinetic measurements confirmed that, at the beginning of the reaction, the primary source of the CO2 release is the cEC cleavage (decomposition).The schematic presentation of the experimentally observed reaction network for SnO3  In summary, the different patterns observed for stannate/vanadate vs. phosphate reflect the different natures of the chemical interactions between the catalyst and cEC (non-similar catalyst characteristics).The more-nucleophilic character of the catalyst anion results in the sampling of a concerted pathway of the catalyst attachment.The following section compares the calculated results with model parameters determined from kinetic measurements to address the correlations/differences between simulations and experiments.

Experimental Analysis of the Kinetics of Polymerisation of cEC, including CO 2 Release
The experimental investigations on the kinetics of the polymerisation of cEC indicate that higher CO 2 content can be obtained by the stannate and vanadate anions, which confirms the better performance of these catalysts for the ROP of cEC.Notably, the kinetic measurements confirmed that, at the beginning of the reaction, the primary source of the The experimental investigations on the kinetics of the polymerisation of cE that higher CO2 content can be obtained by the stannate and vanadate anions, firms the better performance of these catalysts for the ROP of cEC.Notably, measurements confirmed that, at the beginning of the reaction, the primary so CO2 release is the cEC cleavage (decomposition).The schematic presentation o imentally observed reaction network for SnO3 2− and VO4 3− anions is shown in and for PO4 3− in Scheme 5. k1 and k2 in Scheme 4 represent the rate constants of ing polymerisation of cEC and the chain propagation with CO2 release, respec k3 and k4 show the polymeric chain (PECn) decay due to decomposition or C from the chain, respectively.The kinetics of these reactions can be described with first-order rate Equ (4).The determination of the conversion and the CO2 content was performed via the 1 H-NM of the products dissolved in CD2Cl2 or CDCl3 and the measurement of 31 P-NMR in CD2C (see Figure 8).Similar to the case of stannate, k3, and k4 are also negligible for phosphat catalysis.As Table 4 shows, the ratio of k2/k1 is very different from what was observed wit stannate.For the first entry, it is close to one, which means the propagation reaction wit CO2 release and the ring-opening polymerisation have similar rates.However, for entr two, the ring-opening polymerisation is faster than the propagation with CO2 release.Th difference can also be seen from the high CO2 content.However, since the reaction did no proceed to full conversion and achieved only 7% yield, this measurement is highly sus ceptible to experimental error.For entries 3 and 4, the rate of the reaction with CO2 releas is faster than cEC ring-opening polymerisation, producing a polymer with much lowe The kinetics of these reactions can be described with first-order rate Equations (1)-(4).
The material balances of the involved components were combined from these rate equations and solved numerically.Using a Runge-Kutta algorithm, the rate constants were determined via fitting the model to the experimental data.For fitting the absorption curves, the concentrations of the corresponding species were multiplied with the extinction coefficients of these species, estimated from single component spectra (Figure 7 and Table 3).(For a detailed description of the parameter fitting, see Supporting Information, Tables S2-S6.For representative results of fitting the kinetic model to experimental data for all experiments at different temperatures, see Supporting Information, Figure S13-S15).
were determined via fitting the model to the experimental data.For fitting the absorption curves, the concentrations of the corresponding species were multiplied with the extinction coefficients of these species, estimated from single component spectra (Figure 7 and Table 3).(For a detailed description of the parameter fitting, see Supporting Information, Tables S2  Table 3 reports the experimentally measured rate constants corresponding to the reaction network in Scheme 4 (for SnO3 2− anion as the catalyst) at various temperatures.As seen in Table 3, k3 and k4 are negligibly low.Hence, the primary source of the CO2 release is a result of the interaction between monomer and catalyst (CO2 release from cEC).This result agrees with the mechanism shown in Figure 3, in which the CO2 release is due to the catalyst attachment.In Table 3, we also show the ratio of k2/k1 to compare the rate of CO2 release versus polymerisation of cEC.The relative values show that CO2 release occurs either slightly slower than polymerisation (entries 1-3) or at a similar rate to polymerisation (entry 4).The k2/k1 ratio is not strongly dependent on the temperature.This experimental observation agrees well with the computationally observed mechanism in Figures  Table 3 reports the experimentally measured rate constants corresponding to the reaction network in Scheme 4 (for SnO 3 2− anion as the catalyst) at various temperatures.As seen in Table 3, k 3 and k 4 are negligibly low.Hence, the primary source of the CO 2 release is a result of the interaction between monomer and catalyst (CO 2 release from cEC).This result agrees with the mechanism shown in Figure 3, in which the CO 2 release is due to the catalyst attachment.In Table 3, we also show the ratio of k 2 /k 1 to compare the rate of CO 2 release versus polymerisation of cEC.The relative values show that CO 2 release occurs either slightly slower than polymerisation (entries 1-3) or at a similar rate to polymerisation (entry 4).The k 2 /k 1 ratio is not strongly dependent on the temperature.This experimental observation agrees well with the computationally observed mechanism in Figures 4 and 5 (for stannate and vanadate anions).Hence, the possibility of proceeding via catalyst attachment followed by the CO 2 release or returning to the closed cEC and further involvement in the main polymerisation path (propagation via alkoxide attack) is identified.
Furthermore, the calculated FESs in Figure 6 show similar values of the ∆G ̸ = for converting the intermediate to the CAT-CH 2 bond and forming the reactant state (18.5 vs. 18.7 kcal/mol, respectively).According to the experimental measurements, a similar trend of the k 1 and k 2 rate constant ratios can be observed for VO 4 3− , which implies similar catalytic characteristics of the VO 4 3− and SnO 3 2− anions.Additionally, as shown in Figure 6, in the middle part, for the backward path, i.e., returning to the reactant state from the intermediate state, the ∆G ̸ = is more favoured (i.e., 9.2 kcal/mol from intermediate to the reactant state).The results of kinetic measurements of VO 4  3− and the rate constants are reported in Table S1 in SI.
For the PO 4 3− anion, the kinetic measurement results differ from those with SnO 3 2− and VO 4 3− as catalysts.The experimentally observed network of the reactions catalysed by PO 4 3− is shown in Scheme 5.In this Scheme, k des indicates the rate constant of catalyst deactivation, which is observed for PO 4  3− .The results of the kinetic measurements and the rate constants are reported in Table 4.The literature (also shown in Table 4) shows that the ROP of cEC catalysed by phosphate has lower CO 2 content in the final polymer (10% lower than K 3 VO 4 and Na 2 SnO 3 ) [46].In addition to that, according to 31 P-NMR, insertion of phosphate is observed during the polymerisation ( 31 P NMR (500 MHz, CD 2 Cl 2 ) δ 17.84 (s, 1P), 1.64 (d, J = 169.7 Hz, 7P), 1.08 (s, 2P)), shown in Figure 8.The peaks at 0-2 ppm can be assigned to potassium phosphate [47].However, there is a signal at 17.84 ppm, which does not correspond to the signals of potassium phosphate.Instead, signals related to phosphonates can be found in this range [47].4 and 5 (for stannate and vanadate anions).Hence, the possibility of proceeding via catalyst attachment followed by the CO2 release or returning to the closed cEC and further involvement in the main polymerisation path (propagation via alkoxide attack) is identified.
Furthermore, the calculated FESs in Figure 6 show similar values of the ΔG ≠ for converting the intermediate to the CAT-CH2 bond and forming the reactant state (18.5 vs. 18.7 kcal/mol, respectively).According to the experimental measurements, a similar trend of the k1 and k2 rate constant ratios can be observed for VO4 3− , which implies similar catalytic characteristics of the VO4 3− and SnO3 2− anions.Additionally, as shown in Figure 6, in the middle part, for the backward path, i.e., returning to the reactant state from the intermediate state, the ΔG ≠ is more favoured (i.e., 9.2 kcal/mol from intermediate to the reactant state).The results of kinetic measurements of VO4 3− and the rate constants are reported in Table S1 in SI.
For the PO4 3− anion, the kinetic measurement results differ from those with SnO3 2− and VO4 3− as catalysts.The experimentally observed network of the reactions catalysed by PO4 3− is shown in Scheme 5.In this Scheme, kdes indicates the rate constant of catalyst deactivation, which is observed for PO4 3− .The results of the kinetic measurements and the rate constants are reported in Table 4.The literature (also shown in Table 4) shows that the ROP of cEC catalysed by phosphate has lower CO2 content in the final polymer (10% lower than K3VO4 and Na2SnO3) [47].In addition to that, according to 31   In all experiments, below 170 °C, no complete reaction could be observed.Even adding fresh potassium phosphate did not lead to any further reaction.Only above 170 °C is it possible to complete the reaction, since the polymerisation reaction is faster than the deactivation of the catalyst.However, this leads to lower CO2 contents of the products.In all experiments, below 170 • C, no complete reaction could be observed.Even adding fresh potassium phosphate did not lead to any further reaction.Only above 170 • C is it possible to complete the reaction, since the polymerisation reaction is faster than the deactivation of the catalyst.However, this leads to lower CO 2 contents of the products.The determination of the conversion and the CO 2 content was performed via the Polymers 2024, 16, 136 13 of 24 1 H-NMR of the products dissolved in CD 2 Cl 2 or CDCl 3 and the measurement of 31 P-NMR in CD 2 Cl 2 (see Figure 8).Similar to the case of stannate, k 3 , and k 4 are also negligible for phosphate catalysis.
As Table 4 shows, the ratio of k 2 /k 1 is very different from what was observed with stannate.For the first entry, it is close to one, which means the propagation reaction with CO 2 release and the ring-opening polymerisation have similar rates.However, for entry two, the ring-opening polymerisation is faster than the propagation with CO 2 release.This difference can also be seen from the high CO 2 content.However, since the reaction did not proceed to full conversion and achieved only 7% yield, this measurement is highly susceptible to experimental error.For entries 3 and 4, the rate of the reaction with CO 2 release is faster than cEC ring-opening polymerisation, producing a polymer with much lower CO 2 contents.Due to the simultaneous deactivation of the phosphate catalyst, the evaluation of the kinetic parameter is greatly hampered, therefore resulting in large error bars of the deduced values.
Furthermore, the values of k des versus k 1 and k 2 are considered.In that case, a closer correlation between the propagation with CO 2 release (k 2 ) and catalyst deactivation can be found, e.g., for entry four at 425 K (k 2 and k des are very close).In comparison to the computational results, this is in line with the concerted mechanism in the case of phosphate where the catalyst attachment to the cleaved cEC appears faster than in the case of stannate, due to a lower barrier (10 vs. 18.5 kcal/mol, respectively, Figure 6).This situation provides conditions for a faster CO 2 release which is correlated to the lower barrier in simulated FES of K 3 PO 4 for the diagonal path.
Considering the theoretical predictions and experimental results we have shown so far, the main question that still needs to be addressed is how the polymerisation proceeds through the generated molecular complex, in which the catalyst is bound to the monomer via the CH 2 moiety (Figures 4 and 5), and what is the probability of CO 2 release alongside the chain propagation path starting from the catalyst-monomer molecular complex.

Alternative Pathways of Chain Growth Revealed by Simulations
One of the possible chain growth paths starts with the catalyst-anion attached to the monomer is shown in Figure 9.We have examined this path using MTD simulations with the considered CVs depicted in Figure 9.The considered CVs correspond to the COO − nucleophilic attack and cEC cleavage (CV1 and CV2, respectively).

Alternative Pathways of Chain Growth Revealed by Simulations
One of the possible chain growth paths starts with the catalyst-anion attached monomer is shown in Figure 9.We have examined this path using MTD simulation the considered CVs depicted in Figure 9.The considered CVs correspond to the nucleophilic attack and cEC cleavage (CV1 and CV2, respectively).The analysis of the corresponding MTD trajectories for this path at four temper for the three catalysts indicates that, since the COO − terminus has a relatively weak ophilic character than the RO − (due to the negative charge delocalisation on the tw gen atoms), the barrier for the nucleophilic attack to the second cEC increases.O other hand, since the ΔG ≠ for the CH2-O(ethereal) bond cleavage is, on average kcal/mol for various molecular systems and temperatures, the possibility of CO2 upon cEC cleavage enhances before COO − attachment.Alongside examination pathway, the possibility of 30% random CO2 release for some of the trajectories w served (four trajectories ended with CO2 release out of 12 trajectories).Thus, to p CO2 liberation and obtain the MTD trajectory for calculating the ΔG ≠ values, we had The analysis of the corresponding MTD trajectories for this path at four temperatures for the three catalysts indicates that, since the COO − terminus has a relatively weak nucleophilic character than the RO − (due to the negative charge delocalisation on the two oxygen atoms), the barrier for the nucleophilic attack to the second cEC increases.On the other hand, since the ∆G ̸ = for the CH 2 -O(ethereal) bond cleavage is, on average, 8-11 kcal/mol for various molecular systems and temperatures, the possibility of CO 2 release upon cEC cleavage enhances before COO − attachment.Alongside examination of this pathway, the possibility of 30% random CO 2 release for some of the trajectories was observed (four trajectories ended with CO 2 release out of 12 trajectories).Thus, to prohibit CO 2 liberation and obtain the MTD trajectory for calculating the ∆G ̸ = values, we had to fix the corresponding bond lengths leading to CO 2 liberation.These bonds, responsible for CO 2 release and kept fixed during simulations, are highlighted in red in Figure 9.The calculated barriers for this reaction (∆G ̸ = ) are shown in Table 5 for the rate-determining step, which is the nucleophilic attack by COO − to the second cEC monomer (CH 2 moiety).
Table 5.The calculated ∆G ̸ = for the COO − attack to the cEC (rate determining step) at the pathway shown in Figure 9 for the three catalysts calculated at four temperatures.As seen from Table 5, generally, the barriers decrease with increasing temperature, and for K 3 PO 4 , the barriers are lowest.However, chain propagation can proceed via this pathway without the starter ROH, with a higher barrier and more probability of CO 2 liberation than the path including the starter ROH.This finding agrees with the experimental measurements that the majority of CO 2 liberation is caused by cEC decomposition at the beginning of the reaction.In all paths, one reaction coordinate (CV2) is responsible for the ring opening of cEC to convert it to a nucleophile (COO − ) that can further proceed with the ring opening of the next cEC (propagation).

Catalyst/∆G
As an alternative possibility, the insertion of the cEC monomer to the O(CAT)-CH 2 (cEC) bond, which is already formed in the attached catalyst-monomer complex, is examined.With the insertion paths, we estimate if a cEC can perform as a nucleophilic species to cleave the O(CAT)-CH 2 (cEC) bond.The insertion paths are shown in Figure 10.In the insertion paths I and II, we try to cleave the cEC ring and convert it to a nucleophile.However, due to the inherent weak nucleophilicity of cEC by itself, CO 2 release alongside the insertion of cEC into the O(CAT)-CH 2 bond is observed with a predominant possibility (in 80% of the trajectories).
This leads us to the conclusion that cEC is inherently a more electrophilic species than nucleophilic, especially when it interacts with the cations in the surrounding environment that pre-activate the cEC ring.These observations introduce an alternative pathway to the previously proposed mechanism, in which cEC attacks as a nucleophile to a second cEC (Scheme 6) [6].This leads us to the conclusion that cEC is inherently a more electrophilic species th nucleophilic, especially when it interacts with the cations in the surrounding environm that pre-activate the cEC ring.These observations introduce an alternative pathway to previously proposed mechanism, in which cEC attacks as a nucleophile to a second c (Scheme 6) [6].Scheme 6.The previously proposed mechanism of cEC polymerisation, in which cEC attacks a nucleophile to a second cEC [6].Scheme 6.The previously proposed mechanism of cEC polymerisation, in which cEC attacks as a nucleophile to a second cEC [6].Scheme 6.The previously proposed mechanism of cEC polymerisation, in which cEC attacks as nucleophile to a second cEC [6].
Figure 10.The computationally examined insertion paths starting from the attached CAT-cEC mo lecular complex, considering the second cEC as a nucleophile.The CO2 release is the predominan possibility when the cEC is considered a nucleophile.
In insertion path III, the SN2 mechanism via the nucleophilic attack by the O(CAT) to the CH2 moiety of a neighbouring cEC is examined.We have run the insertion paths fo In insertion path III, the S N 2 mechanism via the nucleophilic attack by the O(CAT) to the CH 2 moiety of a neighbouring cEC is examined.We have run the insertion paths for all three catalysts at four temperatures and, after analysing all trajectories, besides CO 2 release, two new pathways are revealed based on the MTD trajectories.We note that these new revealed paths were unknown.In one observed path, the second O of the PO 4 3− anion attacks the CH 2 of the second cEC, which means that chain growth can occur through a new channel via the second cEC cleavage.In the second observed path, the catalyst detachment for SnO 3 2− and VO 4 3− occurred.The two new observed paths are shown in Figure 11A,B.Path A was mainly observed for phosphate and path B for stannate and vanadate catalysts.The variation of the CVs versus the time evolution is depicted in Figure 12 for the trajectories with the two alternative paths.In addition to the CVs, the variation of the non-CV distances that lead to the new-unknown paths in Figure 11 is also depicted in Figure 12.As can be seen in Figure 12, alternative bond formation/dissociation can be probed in these trajectories instead of the considered CVs. Figure 12   the possibility of these methods of polymerisation and for the chain transfer (catalyst detachment) for the three catalysts throughout the calculation of the corresponding MD trajectories and ΔG ≠ values.Since the barriers of the catalyst's second oxygen attack are lower than the catalyst detachment (Table 6), the second nucleophilic attack begins faster and can promote the chain transfer/catalyst detachment.This finding means that with the formation of the second cleaved cEC and the start of a new polymerisation channel, the older generated chain leaves the catalyst, leading to polymer chains with specific shorter lengths and constant molecular weights.To calculate the ∆G ̸ = for the two observed paths, we considered the appropriate CVs and simulated the MTD trajectories at 423 K.These CVs are shown in Figure 13, and the calculated ∆G ̸ = values are reported in Table 6.The ∆G ̸ = values correlate to the rate-determining step, corresponding to the nucleophilic attack by the O − (CV1) in each pathway.In the chain transfer/catalyst-detachment path in Figure 13, we performed the transfer of a partially formed chain with an external nucleophilic attack (RO − ).As seen in Table 6, the catalyst's second oxygen attack has a lower barrier in the case of K 3 PO 4 than K 3 VO 4 and Na 2 SnO 3 , in agreement with the stronger nucleophilicity of the phosphate anion.
The chain transfer through the catalyst detachment has a high barrier, which indicates that the generated chain is firmly attached to the catalyst.Nevertheless, chain transfer occurs at a specific time during the growth period of each individual chain.
Since the barriers of the catalyst's second oxygen attack are lower than the catalyst detachment (Table 6), the second nucleophilic attack begins faster and can promote the chain transfer/catalyst detachment.This finding means that with the formation of the second cleaved cEC and the start of a new polymerisation channel, the older generated chain leaves the catalyst, leading to polymer chains with specific shorter lengths and constant molecular weights.6.
These insights obtained using calculations are well in line with the experimen servations and results.The calculated rate constants of the chain transfer and ca detachment (ktrans), according to the kinetic measurements in the experiment, are re in Table 7 compared to the rate constants of the ring-opening polymerisation and th agation with CO2 release (Tables 3 and 4).As seen in Table 7, the ktrans (chain trans constant) is approximately an order of magnitude smaller than the k1 and k2 (the ra stants of polymerisation).This observation agrees with the calculated high ΔG ≠ val catalyst detachment in Table 6.
For the reaction with the catalyst systems K2SnO3 and K3PO4, LC-ESI-MS o mass spectrometry measurements were made from samples withdrawn during th tion at defined time intervals.These mass spectra are shown in Figures 14 and 15 6.
These insights obtained using calculations are well in line with the experimental observations and results.The calculated rate constants of the chain transfer and catalystdetachment (k trans ), according to the kinetic measurements in the experiment, are reported in Table 7 compared to the rate constants of the ring-opening polymerisation and the propagation with CO 2 release (Tables 3 and 4).As seen in Table 7, the k trans (chain transfer rate constant) is approximately an order of magnitude smaller than the k 1 and k 2 (the rate constants of polymerisation).This observation agrees with the calculated high ∆G ̸ = values for catalyst detachment in Table 6.For the reaction with the catalyst systems K 2 SnO 3 and K 3 PO 4 , LC-ESI-MS orbitrap mass spectrometry measurements were made from samples withdrawn during the reaction at defined time intervals.These mass spectra are shown in Figures 14 and 15.It can be observed from Figures 14 and 15 that after a short reaction time, the masses of the polymers primarily correspond to the masses present at the end of the reaction.Substantial changes in the molecular weight during the progress of the reaction were not observed, indicating no influence of the decreasing monomer concentration on the degree of polymerisation.Therefore, the ring-opening polymerisation reaction behaves similarly to a free radical polymerisation with strong transfer to monomer.Such polymerisation reactions are primarily defined by initiation, growth, termination, and transfer reactions.The termination itself occurs through combination and transfer reactions.If the transfer step is the predominant step for the termination of the growth of the individual chain, the average degree of polymerisation is given by the ratio of the propagation rate to the trans-  It can be observed from Figures 14 and 15 that after a short reaction time, the masses of the polymers primarily correspond to the masses present at the end of the reaction.Substantial changes in the molecular weight during the progress of the reaction were not observed, indicating no influence of the decreasing monomer concentration on the degree of polymerisation.Therefore, the ring-opening polymerisation reaction behaves similarly to a free radical polymerisation with strong transfer to monomer.Such polymerisation reactions are primarily defined by initiation, growth, termination, and transfer reactions.The termination itself occurs through combination and transfer reactions.If the transfer step is the predominant step for the termination of the growth of the individual chain, the average degree of polymerisation is given by the ratio of the propagation rate to the trans- It can be observed from Figures 14 and 15 that after a short reaction time, the masses of the polymers primarily correspond to the masses present at the end of the reaction.Substantial changes in the molecular weight during the progress of the reaction were not observed, indicating no influence of the decreasing monomer concentration on the degree of polymerisation.Therefore, the ring-opening polymerisation reaction behaves similarly to a free radical polymerisation with strong transfer to monomer.Such polymerisation reactions are primarily defined by initiation, growth, termination, and transfer reactions.The termination itself occurs through combination and transfer reactions.If the transfer step is the predominant step for the termination of the growth of the individual chain, the average degree of polymerisation is given by the ratio of the propagation rate to the transfer rate.If both reactions only depend on catalyst and monomer concentration, the concentrations cancel out, and Pn is given by the propagation rate constant to the transfer rate constant (see Equation ( 5)).

Pn =
M polymer M Monomer = k prop.k transfer (5) Using LC-ESI-MS orbitrap mass spectrometry, the number average molecular weights of the polymers can be obtained.Since an ethylene oxide unit (EO) is always incorporated into the polymer chain during a propagation step, the mass of CO 2 can be subtracted from the number average molar mass (M n ) after determining the CO 2 content via 1 H-NMR.These masses can be used to calculate the degree of polymerisation via dividing them by 44 g/mol (M ethylene oxide).
Since the propagation rate is defined by the reaction rate constants k 1 and k 2 , which describe the propagation with and without CO 2 incorporation, resulting in the CO 2 selectivity, these constants can be added to the rate constant of chain growth (k w ).Using these LC-ESI-MS data and Equation ( 5), the rate constant for the transfer step can be determined, as shown in Table 7.
Plotting the degree of polymerisation against the catalyst concentration, it is apparent that the degree of polymerisation is constant within the experimental error.Since the catalysts are bases, this agrees with the observations of Rokicki et al. [6].Furthermore, the calculated barriers of nucleophilic attack by the catalyst anions to the partially formed polymer chains are 25 and 38 kcal/mol for PO 4 3− and SnO 3 2− , respectively, supporting chain degradation at high temperatures triggered by an increased rate of back binding and acceleration of the transfer reaction.See Figure 16.
Polymers 2024, 16, x FOR PEER REVIEW 23 of 27 fer rate.If both reactions only depend on catalyst and monomer concentration, the concentrations cancel out, and Pn is given by the propagation rate constant to the transfer rate constant (see Equation ( 5)).
Using LC-ESI-MS orbitrap mass spectrometry, the number average molecular weights of the polymers can be obtained.Since an ethylene oxide unit (EO) is always incorporated into the polymer chain during a propagation step, the mass of CO2 can be subtracted from the number average molar mass (Mn) after determining the CO2 content via 1 H-NMR.These masses can be used to calculate the degree of polymerisation via dividing them by 44 g/mol (M ethylene oxide).
Since the propagation rate is defined by the reaction rate constants k1 and k2, which describe the propagation with and without CO2 incorporation, resulting in the CO2 selectivity, these constants can be added to the rate constant of chain growth (kw).Using these LC-ESI-MS data and Equation ( 5), the rate constant for the transfer step can be determined, as shown in Table 7.
Plotting the degree of polymerisation against the catalyst concentration, it is apparent that the degree of polymerisation is constant within the experimental error.Since the catalysts are bases, this agrees with the observations of Rokicki et al. [6].Furthermore, the calculated barriers of nucleophilic attack by the catalyst anions to the partially formed polymer chains are 25 and 38 kcal/mol for PO4 3− and SnO3 2− , respectively, supporting chain degradation at high temperatures triggered by an increased rate of back binding and acceleration of the transfer reaction.See Figure 16.By plotting the rate constants of the propagation and transfer reaction (Figure 17), it can be seen that the rate constant of the transfer reaction is smaller than the growth reaction by a factor of about 6 to 10.This difference makes the reaction possible in the first By plotting the rate constants of the propagation and transfer reaction (Figure 17), it can be seen that the rate constant of the transfer reaction is smaller than the growth reaction by a factor of about 6 to 10.This difference makes the reaction possible in the first place but inevitably reduces the degree of polymerisation.This finding explains the variations in the molecular weights with different catalysts already observed in the literature [6,13].

Conclusions
In this work, we have investigated alternative paths of the ring-opening polymerisation of cEC monomer with and without CO2 release, originating from the catalyst-monomer interaction in the propagation step.The kinetic measurements, in combination with metadynamics simulations, indicate that the mechanism of chain growth (starting from the catalyst-monomer molecular complex with the catalyst attached to the cleaved cEC monomer) has a higher barrier with a more significant probability for CO2 release than the two-step mechanism involving the ROH starter molecule in the initiation step.We examined the insertion of the cEC monomer as a nucleophile into the CAT-CH2 molecular complex and the detachment of the catalyst from the chain, followed by the chain transfer.The barriers of chain transfer are high; this finding is in accordance with the increase of the possibility of CO2 release throughout the cEC ring cleavage.Our mechanistic investigations confirmed that the cEC monomer can rarely be considered as a nucleophile, but instead predominantly as an electrophilic species, which is in accordance with the anionic ROP of cEC proposed in the literature.These computational findings are supported by the experimental kinetic measurements and the results that the calculated rate constants of the chain transfer and catalyst detachment are almost an order of magnitude higher than the chain propagation and CO2 release.The calculated barriers confirm that the nucleophilic attack throughout the second anionic oxygen of the catalyst anion can promote Figure 17.Left: Arrhenius plots of the growth reactions for the incorporation and non-incorporation of CO 2 and the termination reaction (transfer) when using sodium orthovanadate (Na 3 VO 4 ).The system used: catalyst (0.01 eq.), ethylene glycol (0.10 eq.), cyclic ethylene carbonate (1.00 eq.).Right: Arrhenius plot of the growth reactions for the incorporation and non-incorporation of CO 2 and the termination reaction (transfer) when using potassium stannate (K 2 SnO 3 ).The system used: catalyst (0.01 eq.), ethylene glycol (0.10 eq.), cyclic ethylene carbonate (1.00 eq.).

Conclusions
In this work, we have investigated alternative paths of the ring-opening polymerisation of cEC monomer with and without CO 2 release, originating from the catalyst-monomer interaction in the propagation step.The kinetic measurements, in combination with metadynamics simulations, indicate that the mechanism of chain growth (starting from the catalystmonomer molecular complex with the catalyst attached to the cleaved cEC monomer) has a higher barrier with a more significant probability for CO 2 release than the two-step mechanism involving the ROH starter molecule in the initiation step.We examined the insertion of the cEC monomer as a nucleophile into the CAT-CH 2 molecular complex and the detachment of the catalyst from the chain, followed by the chain transfer.The barriers of chain transfer are high; this finding is in accordance with the increase of the possibility of CO 2 release throughout the cEC ring cleavage.Our mechanistic investigations confirmed that the cEC monomer can rarely be considered as a nucleophile, but instead predominantly as an electrophilic species, which is in accordance with the anionic ROP of cEC proposed in the literature.These computational findings are supported by the experimental kinetic measurements and the results that the calculated rate constants of the chain transfer and catalyst detachment are almost an order of magnitude higher than the chain propagation and CO 2 release.The calculated barriers confirm that the nucleophilic attack throughout the second anionic oxygen of the catalyst anion can promote chain transfer and catalyst detachment.This realisation means that with the formation of the second cleaved cEC and the start of a new polymerisation channel, the former generated chain leaves the catalyst.This transfer can lead to polymer chains with specific shorter lengths and keeps the molecular weights constant.The results obtained in this study can help in designing novel catalysts with improved characteristics to produce polymers with higher CO 2 content.For higher molecular weights of the polymers, catalysts with higher basicity should be applied, with a lower tendency to detach from the growing chain.

Supplementary Materials:
The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/polym16010136/s1.Figures S1-S3 show the variation of the two CVs versus each other for the three catalysts at four temperatures; Figures S4-S6 show the variation of CV1 and CV2 versus time corresponding to Figure 3  Figure S10 shows the variation of the CV1 and CV2 for the trajectories of the upper part of Figure 13 and Figure S11 shows the variation of the CV1 and CV2 for the lower part of Figure 13; Table S1 shows experimental measurements and rate constants regarding vanadate catalyst.From Section S3 onwards, the supporting information for kinetic measurements is presented.The reactor used is shown in Figure S12.The experimental parameters used are shown in Tables S2-S6.Section S4 deals with the results of the simulation, showing both the fit of the calculated adsorption of the kinetic model to the experimental values and the resulting concentration-time curves for the catalysts stannates (Figure S13), phosphates (Figure S14), and orthovanadates (Figure S15).Section S5 shows the NMR spectra of the 1 H-NMR of the products of the kinetic measurement.The NMR spectra for stannates are in Figures S16-S19 for phosphates, Figures S20-S23, and for the orthovanadate catalyst, Figures S24-S27.Funding: We gratefully acknowledge the RWTH-HPC (project number p0020212) for computation time and technical support.This research work has been carried out within the project DreamRe-sourceConti (033R222A), which has received public funding from the BMBF (Bundesminustrerim für Bildung und Forschung).Technical University of Berlin work has been carried out within the project "Dream ResourceConti" (033R222C).The German Federal Ministry of Education and Research (BMBF) funded the project within the funding priority "CO 2 Plus-Stoffliche Nutzung von CO 2 zur Verbesserung der Rohstoffbasis".

Scheme 3 .
Scheme 3. The considered molecular cluster for simulations; in this model cluster, there are six cEC molecules, one starter alcohol, and one typical catalyst molecule.

Scheme 3 .
Scheme 3. The considered molecular cluster for simulations; in this model cluster, there are six cEC molecules, one starter alcohol, and one typical catalyst molecule.

Figure 1 .
Figure 1.Example 3D spectrum of the in situ measurement of the reaction at 600RPM 140 °C.

Figure 1 .
Figure 1.Example 3D spectrum of the in situ measurement of the reaction at 600RPM 140 • C.

Figure 2 .
Figure 2. The considered mechanism of the ROP of cEC in a previous study.Step I, the proton abstraction by the catalyst anion and nucleophilic attack by the RO − to the carbonyl carbon of the cEC.Step II, the ring opening of the second cEC via the second nucleophilic attack by alkoxide and ring opening.

Figure 2 .
Figure 2. The considered mechanism of the ROP of cEC in a previous study.Step I, the proton abstraction by the catalyst anion and nucleophilic attack by the RO − to the carbonyl carbon of the cEC.Step II, the ring opening of the second cEC via the second nucleophilic attack by alkoxide and ring opening.

Figure 3
Figure 3 depicts the considered CO 2 release mechanism triggered by the interaction between the catalyst anion and the CH 2 of the cEC ring.

Figure 3
Figure 3 depicts the considered CO2 release mechanism triggered by the interaction between the catalyst anion and the CH2 of the cEC ring.

Figure 3 .
Figure 3.The nucleophilic attack by the catalyst anion to the CH2 of the cEC ring results in a terminal COO − and leads to the CO2 release.CV1 and CV2 are the distances used for MTD simulation.Figure 3. The nucleophilic attack by the catalyst anion to the CH 2 of the cEC ring results in a terminal COO − and leads to the CO 2 release.CV1 and CV2 are the distances used for MTD simulation.

Figure 3 .
Figure 3.The nucleophilic attack by the catalyst anion to the CH2 of the cEC ring results in a terminal COO − and leads to the CO2 release.CV1 and CV2 are the distances used for MTD simulation.Figure 3. The nucleophilic attack by the catalyst anion to the CH 2 of the cEC ring results in a terminal COO − and leads to the CO 2 release.CV1 and CV2 are the distances used for MTD simulation.

Figure 4 .
Figure 4. Variation of the O −… C(CH2-O) distance (CV1) versus the C(H2) … O(ethereal) distance (CV2).Different regions can be identified based on the density of the dots (transient structures).The solid red bond shows the CH2-O(cEC) that is cleaved via nucleophilic attack by the catalyst anion and the dashed red line indicates the CH2-O(CAT) bond that is formed.

Figure 5 .
Figure 5. Left: various regions according to the density of the transient states in Figure 4. Right: the schematic structures that correlate to the highlighted regions on the left side.

Figure 4 .
Figure 4. Variation of the O − • • • C(CH 2 -O) distance (CV1) versus the C(H 2 )• • • O(ethereal) distance (CV2).Different regions can be identified based on the density of the dots (transient structures).The solid red bond shows the CH 2 -O(cEC) that is cleaved via nucleophilic attack by the catalyst anion and the dashed red line indicates the CH 2 -O(CAT) bond that is formed.In the case of Na 2 SnO 3 and K 3 VO 4 , the mechanism begins with ring cleavage (ethereal CH 2 -O bond cleavage and intermediate formation) by passing a vertical transition state region (TS vert.).The intermediate either returns to the closed cEC or releases CO 2 .In the case of K 3 PO 4 , an additional path can be identified: the cEC cleavage and nucleophilic attack by the catalyst anion happen simultaneously, and a diagonal transition state region is observed (TS Nu−attack ).These two distinguishable patterns can be identified in the calculated FESs according to the relative energies in kcal/mol.Figure6shows the calculated FESs at 423 K for the three catalysts.As seen in Figure6, the stepwise path, which is the transition intermediate, can be observed in the pattern of the FES of Na 2 SnO 3 and K 3 VO 4 .

Figure 4 .
Figure 4. Variation of the O −… C(CH2-O) distance (CV1) versus the C(H2) … O(ethereal) distance (CV2).Different regions can be identified based on the density of the dots (transient structures).The solid red bond shows the CH2-O(cEC) that is cleaved via nucleophilic attack by the catalyst anion and the dashed red line indicates the CH2-O(CAT) bond that is formed.

Figure 5 .
Figure 5. Left: various regions according to the density of the transient states in Figure 4. Right: the schematic structures that correlate to the highlighted regions on the left side.

Figure 5 . 27 Figure 6 .
Figure 5. Left: various regions according to the density of the transient states in Figure 4. Right: the schematic structures that correlate to the highlighted regions on the left side.Polymers 2024, 16, x FOR PEER REVIEW 11 of 27 2− and VO4 3− anions is shown in Scheme 4 and for PO4 3− in Scheme 5. k1 and k2 in Scheme 4 represent the rate constants of ring-opening polymerisation of cEC and the chain propagation with CO2 release, respectively, and

Figure 6 .
Figure 6.The FESs of the CO 2 release mechanism and the energetics of the regions correspond to the areas depicted in Figure 5 for the three catalysts K 3 VO 4 , Na 2 SnO 3 , and K 3 PO 4 at 423 K.The energy values (numbers) are in kcal/mol.Two different patterns can be identified in the case of K 3 VO 4 and Na 2 SnO 3 (stepwise via an intermediate) versus K 3 PO 4 (concerted).
is the cEC cleavage (decomposition).The schematic presentation of the experimentally observed reaction network for SnO 3 2− and VO 4 3− anions is shown in Scheme 4 and for PO 4 3− in Scheme 5. k 1 and k 2 in Scheme 4 represent the rate constants of ringopening polymerisation of cEC and the chain propagation with CO 2 release, respectively, and k 3 and k 4 show the polymeric chain (PECn) decay due to decomposition or CO 2 release from the chain, respectively.

r = k • c Scheme 4 .
Scheme 4. The experimentally observed reaction network of cEC ROP, catalysed by SnO 3 2− and VO 4 3− anions.k 1 , k 2 , k 3 , and k 4 indicate the corresponding rate constants of ring-opening polymerisation (PECn), chain propagation with CO 2 release, decomposition of the chain, and CO 2 release from the chain, respectively.

Scheme 5 .
Scheme 5.The experimentally observed reaction network was catalysed by PO4 3− anions.k1, k2, k and k4 indicate the corresponding rate constants of polymerization, propagation with CO2 releas decomposition of the chain, and CO2 release from the chain, respectively.kdes indicates the rate con stant of catalyst deactivation.

Scheme 5 .
Scheme 5.The experimentally observed reaction network was catalysed by PO 4 3− anions.k 1 , k 2 , k 3 , and k 4 indicate the corresponding rate constants of polymerization, propagation with CO 2 release, decomposition of the chain, and CO 2 release from the chain, respectively.k des indicates the rate constant of catalyst deactivation.

Figure 7 .
Figure 7. Representative results of fitting the kinetic model to the experimental data.Experimental conditions: TR = 170 °C, ccEC = 13.83mol/L, cEthylenglycol = 1.38 mol/L, 0.01 eq.Potassium stannate.(A) Concentrate against the time of the reactants and products.(B) Fit the IR adsorption bands of the experiment to the experimental results.

Figure 7 .
Figure 7. Representative results of fitting the kinetic model to the experimental data.Experimental conditions: T R = 170 • C, c cEC = 13.83mol/L, c Ethylenglycol = 1.38 mol/L, 0.01 eq.Potassium stannate.(A) Concentrate against the time of the reactants and products.(B) Fit the IR adsorption bands of the experiment to the experimental results.

Figure 9 .
Figure 9.The computationally considered pathway for chain growth starts from the catalys omer molecular complex, while the catalyst is attached to the opened cEC.

Figure 9 .
Figure 9.The computationally considered pathway for chain growth starts from the catalyst-monomer molecular complex, while the catalyst is attached to the opened cEC.
mers 2024, 16, x FOR PEER REVIEW 17 of As an alternative possibility, the insertion of the cEC monomer to the O(CA CH2(cEC) bond, which is already formed in the attached catalyst-monomer complex examined.With the insertion paths, we estimate if a cEC can perform as a nucleoph species to cleave the O(CAT)-CH2(cEC) bond.The insertion paths are shown in Figure In the insertion paths I and II, we try to cleave the cEC ring and convert it to a nucleoph However, due to the inherent weak nucleophilicity of cEC by itself, CO2 release alongs the insertion of cEC into the O(CAT)-CH2 bond is observed with a predominant possibil (in 80% of the trajectories).

Figure 10 .
Figure10.The computationally examined insertion paths starting from the attached CAT-cEC molecular complex, considering the second cEC as a nucleophile.The CO 2 release is the predominant possibility when the cEC is considered a nucleophile.
upper part shows the second O(CAT)-CH 2 bond is formed as soon as the cEC ring is cleaved (blue line).As shown for vanadate and stannate, two new bonds between the O(CAT)• • • CH 2 + and O − (cEC)• • • central atom (V or Sn) are formed when the ring is cleaved (blue and orange lines, respectively).Observation of the paths shown in Figure 11 motivated us to examine the possibility of these methods of polymerisation and for the chain transfer (catalyst detachment) for the three catalysts throughout the calculation of the corresponding MD trajectories and ∆G ̸ = values.

Figure 11 .
Figure 11.(A) The observed pathway alongside the MTD simulation (instead of the CVs showed in Figure 10); the second oxygen of the catalyst attacks the second cEC, and the chain propagation can also go via the second monomer.(B) The second observed path shows the catalyst detachment from the partially formed chain.

Figure 11 . 27 Figure 12 .Figure 12 .
Figure 11.(A) The observed pathway alongside the MTD simulation (instead of the CVs showed in Figure 10); the second oxygen of the catalyst attacks the second cEC, and the chain propagation can also go via the second monomer.(B) The second observed path shows the catalyst detachment from the partially formed chain.Polymers 2024, 16, x FOR PEER REVIEW 20 of 27

Figure 12 .
Figure 12.The variation of the CVs and non-CV distances alongside the simulation of the insertion paths (Figure 11).Variation of the non-CV distances triggered us to consider the new pathways of catalyst detachment and chain transfer, as shown in Figure 13.CV2 is the CH 2 -O(ethereal) distance inside the cEC, which indicates the cEC cleavage.

Polymers 2024 ,Figure 13 .
Figure 13.Upper part: is the nucleophilic attack by the second O − of the catalyst anion to the cEC, lower part: is the chain transfer and catalyst detachment path.The corresponding bar these trajectories are shown in Table6. .

Figure 13 .
Figure 13.Upper part: is the nucleophilic attack by the second O − of the catalyst anion to the second cEC, lower part: is the chain transfer and catalyst detachment path.The corresponding barriers for these trajectories are shown in Table6.

Polymers 2024 , 27 Figure 14 .
Figure 14.Plot of LC-ESI-MS spectra over time for monitoring the development of the molecular weight of the polymer using K2SnO3 as the catalyst.

Figure 15 .
Figure 15.Plot of LC-ESI-MS spectra over time for monitoring the development of the molecular weight of the polymer using K3PO4 as the catalyst.

Figure 14 . 27 Figure 14 .
Figure 14.Plot of LC-ESI-MS spectra over time for monitoring the development of the molecular weight of the polymer using K 2 SnO 3 as the catalyst.

Figure 15 .
Figure 15.Plot of LC-ESI-MS spectra over time for monitoring the development of the molecular weight of the polymer using K3PO4 as the catalyst.

Figure 15 .
Figure 15.Plot of LC-ESI-MS spectra over time for monitoring the development of the molecular weight of the polymer using K 3 PO 4 as the catalyst.

Figure 16 .
Figure 16.Degree of polymerisation at different concentrations of catalyst.A significant reduction in the degree of polymerisation can be observed by increasing the catalyst concentration and thereby increasing the basicity.Tested system Starter: ethylene glycol (0.1 eq.) Catalyst: K2SnO3 T: 150 °C.

Figure 16 .
Figure 16.Degree of polymerisation at different concentrations of catalyst.A significant reduction in the degree of polymerisation can be observed by increasing the catalyst concentration and thereby increasing the basicity.Tested system Starter: ethylene glycol (0.1 eq.) Catalyst: K 2 SnO 3 T: 150 • C.
in the main text; FiguresS7-S9show the variations of the CV1 and CV2 for the MD trajectories corresponding to Figure9in the main text;

Table 2 .
The activation free energies (∆G ̸ = ) and CO 2 contents for various cations for the PO 4

Table 4 .
Experimentally measured rate constants corresponding to Scheme 5 (for PO4 3− anion as the catalyst).CO2 contents were determined via 1 H-NMR in CD2Cl2.

Table 3 .
At various temperatures, experimentally measured rate constants corresponding to the reaction network in Scheme 4 (for SnO3 2− anion as the catalyst).CO2 values were determined via 1 H-NMR in CD2Cl2.

Table 3 .
At various temperatures, experimentally measured rate constants corresponding to the reaction network in Scheme 4 (for SnO 3 2− anion as the catalyst).CO 2 values were determined via 1 H-NMR in CD 2 Cl 2 .

Table 4 .
Experimentally measured rate constants corresponding to Scheme 5 (for PO 4 3− anion as the catalyst).CO 2 contents were determined via 1 H-NMR in CD 2 Cl 2 .

Table 6 .
The barriers of the depicted pathways in Figure13at 423 K.

Table 7 .
Experimentally measured rate constants for the chain transfer (k trans ) compared to the polymerisation and CO 2 release rate constant for the stannate and vanadate catalysts.