1. Introduction
Free radical polymerization has been widely used to produce a diverse range of commercial products for decades [
1]. With advancements in technology, there is a continuous push to create modified products with consistent characteristics. Acrylamide-acrylate copolymer solutions are replacing more complex mixtures used in consumer personal care items, including in hair gels, creams, and waxes [
2]. Solvent is an essential part of commercial recipes, as it reduces viscosity and improves heat removal from the reaction systems designed to produce copolymers with uniform properties. While many commercial systems use non-polar organic solvents for synthesis, these are being replaced by eco-friendly and biocompatible solvents, such as water and alcohols, particularly for biomedical and personal care product applications [
3,
4]. However, polar solvents can form hydrogen bonds (H-bonds) with monomers and/or disrupt monomer–monomer H-bonding, affecting polymerization kinetics and polymer properties compared to the bulk system or polymerization in a non-polar solvent. With few literature studies of solution radical polymerization kinetics of partially water-miscible monomers available, this work studies the polymerization of methyl acrylate (MA) and
N-
tert-butylacrylamide (t-BuAAm) in a solution mixture of ethanol (EtOH) and H
2O.
Agboluaje et al. [
5] applied the pulsed-laser polymerization–size-exclusion chromatography (PLP-SEC) technique to study the effect of an alcohol/water solvent on the propagation rate coefficient (
) of both MA and a second sparingly soluble monomer, 2-methoxyethyl acrylate (MEA). The
values increased both with a decreasing monomer concentration in the solution and with an increasing fraction of water in the alcohol/water mixtures: at 30 °C, the
for 5 wt% MA (
wmon = 0.05) in an EtOH/H
2O mixture containing a weight fraction of ethanol (
) of 0.75 was 3 times greater than the bulk value, with the value increasing to a factor of 16 higher in pure water [
5]. The experimental data for MA were well-represented by the following:
where
T is the temperature in kelvin, the reference temperature
Tref is 303 K,
R is the universal gas constant in J/mol·K, the weight fraction of the solvent
is calculated on a polymer-free basis, and the Arrhenius equation for
is taken from an IUPAC benchmark study [
6].
Previous studies have shown that the polymerization of acrylamides is also affected by the reaction medium. Ganachaud et al. reported that the
of
N-isopropyl acrylamide increases by a factor of 2 at infinite dilution compared to
= 0.95 in water [
4]. Solvent composition also plays a role, with the overall rate of acrylamide polymerization the highest in water and decreasing with the addition of alcohol to the solvent mixture [
7]. Recently, the PLP-SEC technique has been applied to determine that the
value for t-BuAAm, a monomer with reduced water solubility compared to MA, increases by 15% when
is reduced from 1 to 0.75 for
wmon = 0.10 at 30 °C [
8].
In addition to its influence on chain-growth kinetics in polar solutions, the monomer concentration also impacts the relative importance of backbiting on the polymerization rate in acrylate and acrylamide systems. The intramolecular H-atom abstraction process (i.e., backbiting, with rate coefficient
transfers the radical functionality from the chain-end secondary propagation radical (SPR) to a tertiary carbon-centered mid-chain radical (MCR). As the configuration of the MCR is more stable than that of the SPR due to its tertiary nature, the propagation rate coefficient for the monomer addition to the MCR (
) is much less than that for the addition to the SPR
, resulting in a reduction in the rate of overall monomer consumption. The influence can be captured by defining an effective propagation rate coefficient (
that is lowered significantly compared to
[
9,
10]:
While this expression is written for homopolymerization, backbiting reduces the averaged propagation rate coefficient for copolymerization in a similar fashion. Due to its increased stability, a significant fraction of the total radicals can be MCRs, even at lower temperatures when the total number of short-chain branches (SCBs) formed by the reaction is low (<1% of total polymer units). Both the MCR fraction (which influences rate) and the extent of branching are functions of the monomer concentration and temperature [
11,
12].
In our previous work, we quantified the influence of solvent choice on MA homopolymerization by using an in situ NMR technique to follow monomer conversions for batch reactions conducted in toluene and in an ethanol/water mixture. In addition, the product molar mass distributions (MMDs) and weight-average molar masses (
Mw) were characterized using SEC, and the percentage of short-chain branches (%SCB, the number of branch points per 100 monomer repeat units) was measured using
13C NMR [
13]. The faster rate of polymerization observed in the ethanol/water mixture was attributed to the increased
value (due to hydrogen bonding) in the ethanol/water system relative to toluene [
13]. A kinetic model that included the formation and consumption of MCRs was developed to capture the influence of solvent composition (
between 0.75 and 0.95) and initial MA content (
between 0.10 and 0.40) on monomer conversion profiles, polymer MMDs, and branching levels over a range of operating temperatures and initiator concentrations. A second publication extended the model to MA/t-BuAAm copolymerization, accounting for the backbiting of both monomers [
14]. The influence of the solvent on chain-growth kinetics was determined in a PLP-SEC study of the copolymerization system [
8], with reactivity ratios estimated by fitting the drift in the comonomer composition with overall monomer conversion during batch operation [
14]. The detailed model was verified by comparing predictions to the MA/t-BuAAm copolymer MMDs, comonomer compositions, and monomer conversions measured for an extensive set of batch copolymerizations [
14].
While batch reactors are simple to operate, it is difficult to control the monomer composition and, hence, the rate of polymerization due to the unequal consumption rates of the two monomers, resulting in a broadened copolymer composition distribution. For polymerization in water, the decrease in the monomer concentration with time also leads to an increase in
with conversion [
15,
16]. As composition heterogeneity is not desired for most copolymer applications, semi-batch operation is a common industrial practice for controlling copolymer molecular weight, dispersity, and composition. An added benefit is that controlling the monomer inventory through feeding policy provides a means to control the rate of heat generation in the system and to minimize the potential hazard of a thermal runaway [
17,
18]. Typical operation involves adding monomers continuously to a mixed reactor at a fixed rate such that the reactor volume and polymer content in the vessel increase. The relative amounts of monomer in the reactor are affected by both the polymerization kinetics and the rates of monomer addition, with “starved-feed” operation (i.e., maintaining the instantaneous conversion of monomer to polymer high throughout the feeding period) often used in industry to control the polymer molar mass and copolymer composition [
17].
The ability of a polymerization model to represent both batch and semi-batch operations is an important test of its generality. In a batch system, there is a higher monomer concentration throughout most of the polymerization that favors chain growth (a second-order reaction with a rate proportional to the monomer concentration) over backbiting (a first-order reaction). A direct comparison of branching levels in batch and semi-batch systems has demonstrated this for the 2-hydroxyethyl acrylate (HEA)/butyl acrylate (BA) copolymerization system in an organic solvent, with the study also finding that the influence of H-bonding on backbiting differs due to the lower amount of HEA (the H-bonding monomer) in the semi-batch system [
19]. The influence of relative monomer and solvent levels should be even more important for MA/t-BuAAm copolymerized in ethanol/water, as both the
(Equation (1)) and
(Equation (2)) values in the semi-batch system will differ significantly from those in the batch system. Thus, in this study, we compare the evolution of copolymer properties (MMDs and branching levels) produced by batch and semi-batch operations, as well as verifying that the kinetic model developed from the batch study provides a good representation of semi-batch operation.
2. Materials and Methods
Deuterated solvents CDCl
3 (99% purity) and methanol-D4 (MeOD-D4, 99% purity) were obtained from Sigma-Aldrich (Oakville, ON, Canada), and acetone-D6 (99.9% purity) and acetonitrile-D3 (99.8% purity) were obtained from Cambridge Isotope (Tewksbury, MA, USA); all deuterated solvents were used as received. Fully protonated solvents tetrahydrofuran (THF, ≥99.0% purity, 250 ppm BHT inhibitor, Sigma-Aldrich), ethanol-H6 (EtOH-H6, anhydrous, Commercial Alcohols, Brampton, ON, Canada), methanol-H4 (MeOH-H4, ≥99.8% purity, Sigma-Aldrich), diethyl ether (VWR Chemicals, Mississauga, ON, Canada), hexanes (≥98.5% purity, Fisher Scientific, Ottawa, ON, Canada), and deionized water (H
2O) were also used without purification, as well as the inhibitor hydroquinone (HQ, ≥99.9% HPLC-Grade, Fisher Scientific). The thermal initiator 2,2′-azobis(2-methylpropionitrile) (AIBN, 98% purity, Sigma-Aldrich) was recrystallized in MeOH-H4 before use for lab batch experiments. It has been demonstrated that removing the inhibitor from MA (99% purity, ≤100 ppm monomethyl ether hydroquinone (MeHQ), Sigma-Aldrich) does not influence homopolymerization conversion profiles; thus, MA was used as received throughout this work, without further purification [
13]. t-BuAAm (97% purity, Sigma-Aldrich), a powder at room temperature, was used as received.
Semi-batch experiments were performed at two temperatures, 60 and 70 °C, under a nitrogen atmosphere in a 3-neck 250 mL round-bottom flask (RBF) equipped with a magnetic stir bar. The reactor was maintained at the operating temperature in a thermostatted silicone oil bath heated with an IKA C-MAG HS7 hot plate for a typical reaction time of 150 min. The experimental setup was equipped with a Graham condenser; an oil bubbler for monitoring the nitrogen gas flow; and a needle valve to control the nitrogen flowrate at 30–40 bubbles/min, as observed flowing through the oil bubbler.
A typical semi-batch reaction had a total mass of 120 g, with a portion of the solvent precharged along with the initiator to the vessel, and the remainder of the material (monomer with some solvent required to dissolve t-BuAAm) was fed at a constant rate using a REGLO Ismatec ISM4208 peristaltic pump delivered through Pharmed BPT Tubing (1.52 ID, Cole- Parmer, Montreal PQ). All components were purged for 10 min prior to the reactor and feed content preparation. After precharging with the solvent (EtOH and water) and initiator, the reactor flask was set in a silicone oil bath to heat the contents to the desired temperature under agitation, as verified by inserting a thermometer into the vessel. The reaction was initiated by the start of the monomer feed stream to the reactor.
The pump was calibrated before each reaction to deliver the desired feed at a constant rate over 150 min, although the actual delivery time sometimes varied slightly from this target. The calibration was performed by comparing the amount of fluid delivered over a 10 min period with the calculated amount and then adjusting the pump settings. The densities of MA, EtOH, and H
2O were taken from the literature [
13], and the density of t-BuAAm dissolved in solvent was assumed to be 1 g/cm
3.
For the batch operation, the desired monomer and solvent components were added to the reactor, along with a magnetic stir bar, and purged with nitrogen for 10 min, with the pre-weighed initiator mass added once the desired reaction temperature was reached, as verified by inserting a thermometer in the vessel.
Samples were collected at specified time intervals throughout the reaction into vials containing the HQ inhibitor, and they were immediately chilled in an ice-cold bath. A portion of the aliquots (containing the monomer, polymer, and solvent) was diluted for an NMR analysis of the monomer conversion and instantaneous comonomer composition. The remaining aliquots were dried under an air stream at room temperature and then precipitated with a 50:50 diethyl ether:hexane solution to remove residual unreacted t-BuAAm. The precipitated polymer was then dried again before preparation for an SEC analysis of the polymer MMDs and an
1H NMR analysis of the copolymer composition. The analytical equipment and procedures used, the proton NMR used to determine the overall monomer conversion and instantaneous comonomer composition of the reactor aliquots diluted in CDCl
3, and the SEC used to measure the molecular mass distributions (MMD) are explained in previous studies, including the details of SEC calibration [
13,
14].
As previously documented [
13,
14], the monomer conversion profiles and polymer MMDs obtained from the batch polymerizations conducted in the stirred reactor setup were nearly identical to those obtained using the in situ NMR technique under identical conditions. For copolymerization under semi-batch operation, the overall fraction of MA in a sample can be calculated according to Equation (3), based upon the NMR determinations of the overall monomer conversion (
and the comonomer (
) and copolymer (
) compositions:
The resulting value, , reflects the molar fraction of MA contained in both the unreacted monomer and the copolymer relative to the total amount of monomer (MA and t-BuAAm), a value that should equal the fraction of MA contained in the reactor feed, , if the analytical methods are reliable. For all the copolymerization semi-batch experiments, the values remained within ±0.02 of , indicating that the NMR analyses for the copolymer system are quite precise.
4. Conclusions
In this study, a small-scale mixed lab reactor was used, and associated analytical techniques were developed to study the batch and semi-batch solution radical homo- and co-polymerizations of MA and t-BuAAm in EtOH and EtOH/H2O solvent mixtures. Polymer MMDs broadened and shifted to lower values over the course of the batch reactions as monomer concentrations decreased. In contrast, for the semi-batch operation with the monomer fed at a constant flowrate, the monomer concentration and composition remained relatively constant over the course of the reaction such that the polymer properties ( and Ð, as well as copolymer composition) did not drift with the reaction time. Thus, this study demonstrates the advantages of the semi-batch reaction mode used widely in industry.
For the semi-batch copolymerization, the fraction of MA in the formed copolymer was higher than that in the unreacted comonomer, indicating that MA was preferentially incorporated. The polymerization rates and polymer molar masses both increased as the t-BuAAm fraction increased, due to the lowered influence of backbiting. All trends were well-captured by the model developed to represent the copolymerization system, including the influence of the EtOH/H2O solvent level and composition on the monomer propagation rate coefficients. Previously only tested against batch experimental data, the ability of the model to also represent semi-batch (co)polymerizations conducted at low free-monomer levels that increase values but reduce validates the assumptions made and the set of rate coefficients developed to describe the MA/t-BuAAm system.