Source of Processable Vitrimer Viscosities: Swap Frequencies and Steric Factors

SJ Rukmani and S Kim and MA Rahman and X Zhao and AP Sokolov and T Saito and L Petridis and JM Carrillo and A Savara, MACROMOLECULES, 57, 11020-11029 (2024).

DOI: 10.1021/acs.macromol.4c01943

Vitrimers exhibit high, processable viscosities, where other polymers do not, and are among the most promising polymers for closed-loop material circularity. We sought to investigate the underlying chemical kinetic factors that result in high viscosities for vitrimers, which are crucial to designing vitrimers with tunable viscosity. To interrogate these factors, we achieved the first simulated predictions of real vitrimer viscosities, using a novel kinetic Monte Carlo molecular dynamics method, overcoming the time and length scale gaps to predict experimental bulk viscosities. The vitrimer architecture investigated is based on poly(dimethylsiloxane) chains and vinylogous urethane bond swaps. We probed the effects of the extent of free swapping groups, %F, the activation energy, E-A, and the steric factor, rho. The steric factor is related to the intrinsic reaction probability for molecules with sufficient energy. All three factors were found to be significant, but the role of rho was found to be the biggest and also the most underappreciated. The results show that the inclusion of accurate rho is of critical importance for viscosity predictions, with the evidence suggesting that the typical assumption of rho = 1 is not valid for vitrimers and that, indeed, very low steric factors are present in bond- swap vitrimers such that values of rho < 10(-10) may be typical. This greatly influences the bond exchange rates and, ultimately, the viscosities. Recognition of this result is necessary for the prediction of vitrimer viscosities from molecular simulations and to make vitrimers by design from molecular dynamics. We also investigated the effects that E-A, rho, and the number of free swapping groups have upon vitreous range temperatures, T-V, with respect to achieving a specific viscosity (eta(V) = 1 x 10(8) Pas), as well as for a commonly reported higher viscosity extrapolation (eta(V) = 1 x 10(12) Pas). The evidence suggests that vitrimers may follow universal curves for E-A vs T-V, as a function of rho. This study achieves the first of these comparisons of molecular simulations to experiments and reveals critical insights toward creating vitrimers by design, while providing a route for the prediction of T-V from kinetic Monte Carlo molecular dynamics simulations.

Return to Publications page