Improved Approach for ab Initio Calculations of Rate Coefficients for Secondary Reactions in Acrylate Free-Radical Polymerization

Secondary reactions in radical polymerization pose a challenge when creating kinetic models for predicting polymer structures. Despite the high impact of these reactions in the polymer structure, their effects are difficult to isolate and measure to produce kinetic data. To this end, we used solvation-corrected M06-2X/6-311+G(d,p) ab initio calculations to predict a complete and consistent data set of intrinsic rate coefficients of the secondary reactions in acrylate radical polymerization, including backbiting, β-scission, radical migration, macromonomer propagation, mid-chain radical propagation, chain transfer to monomer and chain transfer to polymer. Two new approaches towards computationally predicting rate coefficients for secondary reactions are proposed: (i) explicit accounting for all possible enantiomers for reactions involving optically active centers; (ii) imposing reduced flexibility if the reaction center is in the middle of the polymer chain. The accuracy and reliability of the ab initio predictions were benchmarked against experimental data via kinetic Monte Carlo simulations under three sufficiently different experimental conditions: a high-frequency modulated polymerization process in the transient regime, a low-frequency modulated process in the sliding regime at both low and high temperatures and a degradation process in the absence of free monomers. The complete and consistent ab initio data set compiled in this work predicts a good agreement when benchmarked via kMC simulations against experimental data, which is a technique never used before for computational chemistry. The simulation results show that these two newly proposed approaches are promising for bridging the gap between experimental and computational chemistry methods in polymer reaction engineering.

Beauty of Explicit Dispersity (Đ) Equations in Controlled Polymerizations

Dispersity (Đ) as a critical parameter indicates the level of uniformity of the polymer molar mass or chain length. In the past several decades, the development of explicit equations for calculating Đ experiences a continual revolution. This viewpoint tracks the historical evolution of the explicit equations from living to reversible-deactivation polymerization systems. Emphasis is laid on displaying the charm of explicit Đ equations in batch reversible-deactivation radical polymerization (RDRP), with highlights of the relevant elegant mathematical manipulations. Some representative emerging applications enabled by the existing explicit equations are shown, involving nitroxide-mediated polymerization (NMP), atom transfer radical polymerization (ATRP), and reversible addition-fragmentation chain transfer (RAFT) polymerization systems. Stemming from the several outlined challenges and outlooks, sustained concerns about the explicit Đ equations are still highly deserved. It is expected that these equations will continue to play an important role not only in traditional polymerization kinetic simulation and design of experiments but also in modern intelligent manufacturing of precision polymers and classroom education.

Analysis of UV Curing Strategy on Reaction Heat Control and Part Accuracy for Additive Manufacturing

In this research, the relationship between the curing strategies and geometrical accuracy of parts under UV light was investigated. An IR camera was utilized to monitor the process using different combinations of photosensitive resin and curing strategies. The influences of curing strategies on different material compositions were studied with single-factor analysis. With the different exposure frequencies of the UV light, the peak temperature was adjusted to avoid overheating. The three-dimensional geometry of casting tensile bars was measured to investigate the shrinkage and warpage during the curing process. Different material compositions were also selected to study the effects of the maximum temperature on the shrinkage of the parts. The findings of this work show that, with the same amount of energy input, a more fragmented exposure allows for a more controllable max temperature, while one-time exposure leads to a high temperature during the process. With the decrease of the released heat from the reaction, the shrinkage of the casting part has a slightly increasing tendency. Moreover, the warpage of the parts decreased drastically with the decrease of temperature. The addition of fillers enhances the control over temperature and increases the geometrical accuracy.

Nanoconfinement in miniemulsion increases reaction rates of thiol–ene photopolymerization and yields high molecular weight polymers

Photoinitiated thiol–ene polymerization was performed in bulk and miniemulsion. We show that the compartmentalization of the reaction inside nanodroplets led to faster reaction kinetics and yielded polymers with higher molecular weight.

State-Of-The-Art Quantification of Polymer Solution Viscosity for Plastic Waste Recycling

Solvent-based recycling is a promising approach for closed-loop recovery of plastic-containing waste. It avoids the energy cost to depolymerize the plastic but still allows to clean the polymer of contaminants and additives. However, viscosity plays an important role in handling the polymer solutions at high concentrations and in the cleaning steps. This Review addresses the viscosity behavior of polymer solutions, available data, and (mostly algebraic) models developed. The non-Newtonian viscosity models, such as the Carreau and Yasuda-Cohen-Armstrong models, pragmatically describe the viscosity of polymer solutions at different concentrations and shear rate ranges. This Review also describes how viscosity influences filtration and centrifugation processes, which are crucial steps in the cleaning of the polymer and includes a polystyrene/styrene case study.

Connecting Gas-Phase Computational Chemistry to Condensed Phase Kinetic Modeling: The State-of-the-Art

In recent decades, quantum chemical calculations (QCC) have increased in accuracy, not only providing the ranking of chemical reactivities and energy barriers (e.g., for optimal selectivities) but also delivering more reliable equilibrium and (intrinsic/chemical) rate coefficients. This increased reliability of kinetic parameters is relevant to support the predictive character of kinetic modeling studies that are addressing actual concentration changes during chemical processes, taking into account competitive reactions and mixing heterogeneities. In the present contribution, guidelines are formulated on how to bridge the fields of computational chemistry and chemical kinetics. It is explained how condensed phase systems can be described based on conventional gas phase computational chemistry calculations. Case studies are included on polymerization kinetics, considering free and controlled radical polymerization, ionic polymerization, and polymer degradation. It is also illustrated how QCC can be directly linked to material properties.

Going Beyond the Carothers, Flory and Stockmayer Equation by Including Cyclization Reactions and Mobility Constraints

A challenge in the field of polymer network synthesis by a step-growth mechanism is the quantification of the relative importance of inter- vs. intramolecular reactions. Here we use a matrix-based kinetic Monte Carlo (kMC) framework to demonstrate that the variation of the chain length distribution and its averages (e.g., number average chain length x_n), are largely affected by intramolecular reactions, as mostly ignored in theoretical studies. We showcase that a conventional approach based on equations derived by Carothers, Flory and Stockmayer, assuming constant reactivities and ignoring intramolecular reactions, is very approximate, and the use of asymptotic limits is biased. Intramolecular reactions stretch the functional group (FG) conversion range and reduce the average chain lengths. In the likely case of restricted mobilities due to diffusional limitations because of a viscosity increase during polymerization, a complex x_n profile with possible plateau formation may arise. The joint consideration of stoichiometric and non-stoichiometric conditions allows the validation of hypotheses for both the intrinsic and apparent reactivities of inter- and intramolecular reactions. The kMC framework is also utilized for reverse engineering purposes, aiming at the identification of advanced (pseudo-)analytical equations, dimensionless numbers and mechanistic insights. We highlight that assuming average molecules by equally distributing A and B FGs is unsuited, and the number of AB intramolecular combinations is affected by the number of monomer units in the molecules, specifically at high FG conversions. In the absence of mobility constraints, dimensionless numbers can be considered to map the time variation of the fraction of intramolecular reactions, but still, a complex solution results, making a kMC approach overall most elegant.

Jacket temperature regulation allowing well-defined non-adiabatic lab-scale solution free radical polymerization of acrylates

Conventional batch solution free radical polymerization of n-butyl acrylate with thermal initiators such as AIBN is known to be strongly exothermic and influenced by highly activated side reactions such as backbiting and β-scission.

Impact of side reactions on molar mass distribution, unsaturation level and branching density in solution free radical polymerization of n-butyl acrylate under well-defined lab-scale reactor conditions

The paper describes the influence of side reactions in isothermal solution free-radical polymerization of n-butyl acrylate accounting for chain-length dependent diffusional limitations on termination.

Coupled matrix kinetic Monte Carlo simulations applied for advanced understanding of polymer grafting kinetics

An innovative coupled matrix-based Monte Carlo (CMMC) concept has been applied to successfully assess the detailed description of the molecular build-up of linear and non-linear chains in the free-radical induced grafting of linear precursors chains.

Progress in Reaction Mechanisms and Reactor Technologies for Thermochemical Recycling of Poly(methyl methacrylate)

Chemical or feedstock recycling of poly(methyl methacrylate) (PMMA) by thermal degradation is an important societal challenge to enable polymer circularity. The annual PMMA world production capacity is over 2.4 × 10⁶ tons, but currently only 3.0 × 10⁴ tons are collected and recycled in Europe each year. Despite the rather simple chemical structure of MMA, a debate still exists on the possible PMMA degradation mechanisms and only basic batch and continuous reactor technologies have been developed, without significant knowledge of the decomposition chemistry or the multiphase nature of the reaction mixture. It is demonstrated in this review that it is essential to link PMMA thermochemical recycling with the PMMA synthesis as certain structural defects from the synthesis step are affecting the nature and relevance of the subsequent degradation reaction mechanisms. Here, random fission plays a key role, specifically for PMMA made by anionic polymerization. It is further highlighted that kinetic modeling tools are useful to further unravel the dominant PMMA degradation mechanisms. A novel distinction is made between global conversion or average chain length models, on the one hand, and elementary reaction step-based models on the other hand. It is put forward that only by the dedicated development of the latter models, the temporal evolution of degradation product spectra under specific chemical recycling conditions will become possible, making reactor design no longer an art but a science.

Connecting polymer synthesis and chemical recycling on a chain-by-chain basis: a unified matrix-based kinetic Monte Carlo strategy

Polymer synthesis and subsequent depolymerisation/degradation are linked at the molecular level.

Visualization and design of the functional group distribution during statistical copolymerization

Even though functional copolymers with a low percentage of functional comonomer units (up to 20 mol%) are widely used, for instance for the development of polymer therapeutics and hydrogels, insights in the functional group distribution over the actual chains are lacking and the average composition is conventionally used to describe the functionalization degree. Here we report the visualization of the monomer distribution over the different polymer chains by a synergetic combination of experimental and theoretical analysis aiming at the construction of functionality-chain length distributions (FUNC-CLDs). A successful design of the chemical structure of the comonomer pair, the initial functional comonomer amount (13 mol%), and the temperature (100 °C) is performed to tune the FUNC-CLD of copoly(2-oxazoline)s toward high functionalization degree for both low (100) and high (400) target degrees of polymerization. The proposed research strategy is generic and extendable to a broad range of copolymerization chemistries, including reversible deactivation radical polymerization.

Modeling of the Free Radical Copolymerization Kinetics of n-Butyl Acrylate, Methyl Methacrylate and 2-Ethylhexyl Acrylate Using PREDICI®

Kinetic modeling of the bulk free radical copolymerizations of n-butyl acrylate (BA) and 2-ethylhexyl acrylate (EHA); methyl methacrylate (MMA) and EHA; as well as BA, MMA and EHA was performed using the software PREDICI®. Predicted results of conversion versus time, composition versus conversion, and molecular weight development are compared against experimental data at different feed compositions. Diffusion-controlled effects and backbiting for BA were incorporated into the model as they proved to be significant in these polymerizations. The set of estimated global parameters allows one to assess the performance of these copolymerization systems over a wide range of monomer compositions.

Modeling of Miniemulsion Polymerization of Styrene with Macro-RAFT Agents to Theoretically Compare Slow Fragmentation, Ideal Exchange and Cross-Termination Cases

A 5-dimensional Smith-Ewart based model is developed to understand differences for reversible addition-fragmentation chain transfer (RAFT) miniemulsion polymerization with theoretical agents mimicking cases of slow fragmentation, cross-termination, and ideal exchange while accounting for chain length and monomer conversion dependencies due to diffusional limitations. The focus is on styrene as a monomer, a water soluble initiator, and a macro-RAFT agent to avoid exit/entry of the RAFT leaving group radical. It is shown that with a too low RAFT fragmentation rate coefficient it is generally not afforded to consider zero-one kinetics (for the related intermediate radical type) and that with significant RAFT cross-termination the dead polymer product is dominantly originating from the RAFT intermediate radical. To allow the identification of the nature of the RAFT retardation it is recommended to experimentally investigate in the future the impact of the average particle size (dp) on both the monomer conversion profile and the average polymer properties for a sufficiently broad dp range, ideally including the bulk limit. With decreasing particle size both a slow RAFT fragmentation and a fast RAFT cross-termination result in a stronger segregation and thus rate acceleration. The particle size dependency is different, allowing further differentiation based on the variation of the dispersity and end-group functionality. Significant RAFT cross-termination is specifically associated with a strong dispersity increase at higher average particle sizes. Only with an ideal exchange it is afforded in the modeling to avoid the explicit calculation of the RAFT intermediate concentration evolution.

β-Myrcene/isobornyl methacrylate SG1 nitroxide-mediated controlled radical polymerization: synthesis and characterization of gradient, diblock and triblock copolymers

β-Myrcene (My), a natural 1,3-diene, and isobornyl methacrylate (IBOMA), from partially bio-based raw materials sources, were copolymerized by nitroxide-mediated polymerization (NMP) in bulk using the SG1-based BlocBuilder™ alkoxyamine functionalized with an N-succinimidyl ester group, NHS-BlocBuilder, at T = 100 °C with initial IBOMA molar feed compositions f IBOMA,0 = 0.10-0.90. Copolymer reactivity ratios were r My = 1.90-2.16 and r IBOMA = 0.02-0.07 using Fineman-Ross, Kelen-Tudos and non-linear least-squares fitting to the Mayo-Lewis terminal model and indicated the possibility of gradient My/IBOMA copolymers. A linear increase in molecular weight versus conversion and a low dispersity (Đ ≤ 1.41) were exhibited by My/IBOMA copolymerization with f IBOMA,0 ≤ 0.80. My-rich and IBOMA-rich copolymers were shown to have a high degree of chain-end fidelity by performing subsequent chain-extensions with IBOMA and/or My, and by 31P NMR analysis. The preparation by NMP of My/IBOMA thermoplastic elastomers (TPEs), mostly bio-sourced, was then attempted. IBOMA-My-IBOMA triblock copolymers containing a minor fraction of My or styrene (S) units in the outer hard segments (M n = 51-95 kg mol-1, Đ = 1.91-2.23 and F IBOMA = 0.28-0.36) were synthesized using SG1-terminated poly(ethylene-stat-butylene) dialkoxyamine. The micro-phase separation was suggested by the detection of two distinct T gs at about -60 °C and +180 °C and confirmed by atomic force microscopy (AFM). A plastic stress-strain behavior (stress at break σ B = 3.90 ± 0.22 MPa, elongation at break ε B = 490 ± 31%) associated to an upper service temperature of about 140 °C were also highlighted for these triblock polymers.