A Density Functional Theory Study of Secondary Reactions in n -Butyl Acrylate Free Radical Polymerization

Unveiling Solvent Effects on β-Scissions through Metadynamics and Mean Force Integration

This study introduces a methodology that combines accelerated molecular dynamics and mean force integration to investigate solvent effects on chemical reaction kinetics. The newly developed methodology is applied to the β-scission of butyl acrylate (BA) dimer in polar (water) and nonpolar (xylene and BA monomer) solvents. The results show that solvation in both polar and nonpolar environments reduces the free energy barrier of activation by ∼4 kcal/mol and decreases the pre-exponential factor 2-fold. Employing a hybrid quantum mechanics/molecular mechanics approach with explicit solvent modeling, we compute kinetic rate constants that better match experimental measurements compared to previous gas-phase calculations. This methodology presents promising potential for accurately predicting kinetic rate constants in liquid-phase polymerization and depolymerization processes.

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.

Theoretical Insights into Midchain Radicals and Branching Characteristics in Solution Copolymerization of Ethylene and Vinyl Acetate

The kinetics of intermolecular chain transfer to polymer (CTP) and chain transfer to monomer (CTM), as well as backbiting in solution copolymerization of ethylene and vinyl acetate (VAc), are calculated by density functional theory. An intrinsic rate coefficient for each reaction type, whose Arrhenius parameters are expressed as functions of unreacted monomer composition, is extracted from calculated elementary reactions, and hence it applies to a wide range of fragment compositions of EVA. Backbiting controls the generation of midchain radicals (MCRs) on the backbone, and its rate maximizes at a medium fraction of VAc in the unreacted monomers. CTP plays an insignificant role in MCR generation even at high conversion rates due to its low pre-exponential factor. The concentration of MCRs is quantified and is close to that of ECRs at high conversions and high fractions of VAc in monomers. Additionally, branching characteristics are explored; concentrations of short-chain branching (SCB) and long-chain branching (LCB) are about 0.7-2.5 and 0.1-0.4%, respectively, and drop with VAc fraction rapidly. The role of the migration of MCRs is highlighted, which transforms about 10% of SCB points into LCB points. Disagreeing with insights from laboratory experiments, it is the migration, rather than CTP, which increases the LCB concentration at high conversions.

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.

DFT Modeling of the Alternating Radical Copolymerization and Alder-Ene Reaction between Maleic Anhydride and Olefins

The free radical copolymerization of electron-acceptor and electron-donor vinyl monomers represents a particular case of sequence-controlled polymerization. The reactions of maleic anhydride (MA) or related compounds (acceptor comonomers) with α-olefins (donor comonomers) result in the formation of the alternating copolymers that have clear prospects for petrochemical and biomedical applications. However, in contrast to the well-established polymerization of acrylate monomers, these processes have not been studied theoretically using the density functional theory (DFT) calculations. In our research, we performed a comprehensive theoretical analysis of the free radical copolymerization of MA and closely related maleimide with different structural types of olefins at mpw1pw91/6-311g(d) level of the DFT. The results of our calculations clearly indicated the preference of the alternating reaction mode for the copolymerization of MA with α-olefins, isobutylene and prospective unsaturated monomers, as well as methylenealkanes. The DFT modeling of the thermally induced Alder-ene reaction between MA and olefins allowed to exclude this reaction from the scope of possible side processes at moderately high temperatures. Comparative analysis of MA and N-methylmaleimide (MMI) reactivity shown that the use of MMI instead of MA makes no sense in terms of the reaction rate and selectivity.

Method of Moments Applied to Most-Likely High-Temperature Free-Radical Polymerization Reactions

Many widely-used polymers are made via free-radical polymerization. Mathematical models of polymerization reactors have many applications such as reactor design, operation, and intensification. The method of moments has been utilized extensively for many decades to derive rate equations needed to predict polymer bulk properties. In this article, for a comprehensive list consisting of more than 40 different reactions that are most likely to occur in high-temperature free-radical homopolymerization, moment rate equations are derived methodically. Three types of radicals—secondary radicals, tertiary radicals formed through backbiting reactions, and tertiary radicals produced by intermolecular chain transfer to polymer reactions—are accounted for. The former tertiary radicals generate short-chain branches, while the latter ones produce long-chain branches. In addition, two types of dead polymer chains, saturated and unsaturated, are considered. Using a step-by-step approach based on the method of moments, this article guides the reader to determine the contributions of each reaction to the production or consumption of each species as well as to the zeroth, first and second moments of chain-length distributions of live and dead polymer chains, in order to derive the overall rate equation for each species, and to derive the rate equations for the leading moments of different chain-length distributions. The closure problems that arise are addressed by assuming chain-length distribution models. As a case study, β-scission and backbiting rate coefficients of methyl acrylate are estimated using the model, and the model is then applied to batch spontaneous thermal polymerization to predict polymer average molecular weights and monomer conversion. These predictions are compared with experimental measurements.

On the Recovery of PLP-Molar Mass Distribution at High Laser Frequencies: A Simulation Study

Due to the inherent difficulties in determination of the degree of branching for polymers produced in pulsed laser polymerization (PLP) experiments, the behavior of the degree of branching and backbiting reaction in high laser frequency and relatively high reaction temperatures have not been well-established. Herein, through a simulation study, the validity of different explanations on the recovery of PLP-molar mass distribution at high laser frequencies is discussed. It is shown that the reduction of the backbiting reaction rate at high laser frequency, and consequent decrease in the degree of branching, is not a necessary condition for recovering the PLP-molar mass distribution. The findings of this work provide simulation support to a previous explanation about the possibility of using high laser frequency for reliable determination of the propagation rate coefficient for acrylic monomers.

Identification of β scission products from free radical polymerizations of butyl acrylate at high temperature

Unsaturated low molar mass species were identified via ESI-MS after fractionation of poly(butyl acrylate) from high temperature radical polymerization.

Characterization of oligo(acrylic acid)s and their block co-oligomers

Oligo(acrylic acid), oligoAA are important species currently used industrially in the stabilization of paints and also for the production of self-assembled polymer structures which have been shown to have useful applications in analytical separation methods and potentially in drug delivery systems. To properly tailor the synthesis of oligoAA, and its block co-oligomers synthesized by Reversible-Addition Fragmentation chain Transfer (RAFT) polymerization to applications, detailed knowledge about the chemical structure is needed. Commonly used techniques such as Size Exclusion Chromatography (SEC) and Electrospray Ionization-Mass Spectrometry (ESI-MS) suffer from poor resolution and non-quantitative distributions, respectively. In this work free solution Capillary Electrophoresis (CE) has been thoroughly investigated as an alternative, allowing for the separation of oligoAA by molar mass and the RAFT agent end group. The method was then extended to block co-oligomers of acrylic acid and styrene. Peak capacities up to 426 were observed for these 1D CE separations, 10 times greater than what has been achieved for Liquid Chromatography (LC) of oligostyrenes. To provide a comprehensive insight into the chemical structure of these materials ¹H and ¹³C Nuclear Magnetic Resonance (NMR) spectroscopy was used to provide an accurate average chain length and reveal the presence of branching. The chain length at which branching is detected was investigated with the results showing a degree of branching of 1% of the monomer units in oligoAA with an average chain length of 9 monomer units, which was the shortest chain length at which branching could be detected. This branching is suspected to be a result of both intermolecular and intramolecular transfer reactions. The combination of free solution CE and NMR spectroscopy is shown to provide a near complete elucidation of the chemical structure of oligoAA including the average chain length and branching as well as the chain length and RAFT agent end group distribution. Furthermore, the purity in terms of the dead chains and unreacted RAFT agent was quantified. The use of free solution CE and ¹H NMR spectroscopy demonstrated in this work can be routinely applied to oligoelectrolytes and their block co-oligomers to provide an accurate characterization which allows for better design of the materials produced from these oligomers.

Mid-Chain Radical Migration in the Radical Polymerization of n-Butyl Acrylate

The occurrence of intramolecular transfer to polymer in the radical polymerization of acrylic monomers has been extensively documented in the literature. Whilst it has been largely assumed that intramolecular transfer to polymer leads to short chain branches, there has been some speculation over whether the mid-chain radical can migrate. Herein, by the matrix-assisted laser desorption/ionization time of flight (MALDI-TOF) mass spectrometry (MS) of poly(n-butyl acrylate) synthesized by solution polymerization under a range of conditions, it is shown that this mid-chain radical migration does occur in the radical polymerization of acrylates conducted at high temperatures, as is evident from the shape of the molecular weight distribution. Using a mathematical model, an initial approximation of the rate at which migration occurs is made and the distribution of branching lengths formed in this scenario is explored. It is shown that the polymerizations carried out under a low monomer concentration and at high temperatures are particularly prone to radical migration reactions, which may affect the rheological properties of the polymer.

Access to the β-scission rate coefficient in acrylate radical polymerization by careful scanning of pulse laser frequencies at elevated temperature

A novel method to estimate the β-scission rate coefficient in acrylate radical polymerization is presented.

Analyzing the discrepancies in the activation energies of the backbiting and β-scission reactions in the radical polymerization of n-butyl acrylate

Activaton energies for backbiting and β-scission reactions for the polymerization of n-BA were determined by fitting experimental data to a mathematical model. The activation energy for backbiting was higher and that for β-scission lower than those accepted in the literature.

Pulsed-laser and quantum mechanics study of n-butyl cyanoacrylate and methyl methacrylate free-radical copolymerization

Improved control over polymer microstructure is achieved by radical copolymerization of n-butyl cyanoacrylate with methacrylates.

Critically evaluated rate coefficients in radical polymerization – 7. Secondary-radical propagation rate coefficients for methyl acrylate in the bulk

Benchmark propagation rate coefficient (k_p) data for the radical polymerization of methyl acrylate are provided.