Free radicals computational approaches

The transient behavior of n was investigated by a number of studies (e.g., see references 13-16). However, the pseudo-steady-state assumption is generally accepted for the mole balance of free radicals in the polymerization system due to the very high reactivity and low concentration of free radicals. This approach can simplify the analysis of the emulsion polymerization kinetics and save computation time significantly. 

As stated above, the thermochemistry of free radicals can also be estimated by the group additivity method, if group values are available. With the exception of a few cases reported in Benson (1976), however, such information presently does not exist. Therefore, we rely on the model compound approach (for S and Cp) and bond dissociation energy (BDE) considerations and computational quantum mechanics for the determination of the heats of formation of radicals. 

EPR Spectra of Organic Free Radicals in Solution from an Integrated Computational Approach... 

Recognition of systems approach generation of databases of kinetics and energetics advances in computation fundamental understanding of chemistry (free radicals) instrumentation (remote sensing, space based, advances in sensor technology)... 

The use of quantum chemistry to obtain the individual rate coefficients of a free-radical polymerization process frees them from errors due to kinetic model-based assumptions. However, this approach introduces a new source of error in the model predictions the quantum chemical calculations themselves. As is well known, as there are no simple analytical solutions to a many-electron Schrodinger equation, numerical approximations are required. While accurate methods exist, they are generally very computationally intensive and their computational cost typically scales exponentially with the size of the system under study. The apphcation of quantum chemical methods to radical polymerization processes necessarily involves a compromise in which small model systems are used to mimic the reactions of their polymeric counterparts so that high levels of theory may be used. This is then balanced by the need to make these models as reahstic as possible hence, lower cost theoretical procedures are frequently adopted, often to the detriment of the accuracy of the calculations. Nonetheless, aided by rapid and continuing increases to computer power, chemically accurate predictions are now possible, even for solvent-sensitive systems [8]. In this section we examine the best-practice methodology required to generate accurate gas- and solution-phase predictions of rate coefficients in free-radical polymerization. 

Applying the mechanistic modeling approach to the petroleum refining scale is a big challenge because of the size and stifftiess (wide range of time and length scales) of the problem. It is not obvious how this approach can be used to describe coke formation which is an important yet incompletely understood process. A computationally less intensive approach is to sidestep transient intermediates (carboeations or free radicals) and focus on dominant reaction pathways. 

Using the results for the moments from this approach, the PDI is computed in Equation 1.51. Because q is the probability of propagation compared to chain inactivation events, the value for q must be very close to 1 for a polymer of any appreciable length to be produced. This finding shows that the PDI for a steady-state free radical polymerization terminated exclusively by disproportionation should be 2. 

FIGURE 13.13 The computed approach index A (Equation 13.64) and instantaneous conversion versus t. Reprinted with permission from Kreft T, Reed WF. Predictive control and verification of conversion kinetics and polymer molecular weight in semibatch free radical homopolymer reactions. Eur Polym J 2009 45 2288-2303. 

The ab initio computation of nuclear hyperhne tensors of small free-radical systems has a long history [43, 47, 68-83]. Our group recently validated a general computational approach rooted in DFT to the analysis of spin-probing and spinlabeling experiments by providing accurate description of thermodynamic and spectroscopic properties of several aliphatic nitroxides as proxy 1 and tempo [56]. The performances of the model for a typical problem were tuned not only by the choice of the right density functional and basis set but also by a proper account of stereoelectronic, vibrational, and environmental effects [56]. 

Chern [42] developed a mechanistic model based on diffusion-controlled reaction mechanisms to predict the kinetics of the semibatch emulsion polymerization of styrene. Reasonable agreement between the model predictions and experimental data available in the literature was achieved. Computer simulation results showed that the polymerization system approaches Smith-Ewart Case 2 kinetics (n = 0.5) when the concentration of monomer in the latex particles is close to the saturation value. By contrast, the polymerization system under the monomer-starved condition is characterized by the diffusion-con-trolled reaction mechanisms (n > 0.5). The author also developed a model to predict the effect of desorption of free radicals out of the latex particles on the kinetics of the semibatch emulsion polymerization of methyl acrylate [43]. The validity of the kinetic model was confirmed by the experimental data for a wide range of monomer feed rates. The desorption rate constant for methyl acrylate at 50°C was determined to be 4 x 10 cm s ... 

There are two possible methods of approach. Either the space integral rates Siu and the associated derivative terms a are evaluated at the beginning of the time interval, and the dependent variables are updated together at the end of the interval or the appropriate 5, and Sjy a are evaluated immediately prior to the solution of each equation and each dependent variable is updated immediately after the solution of its equation. When using the latter method, the order of solution of the equations becomes important. The most satisfactory order treats the most reactive free radicals first, then less reactive intermediates, initial reactants, reaction products and inert species in that order, and finally the enthalpy. The latter method is potentially the more economical computationally, but it forfeits precise chemical conservation during single time steps on the approach to the steady state. This has been found to cause instability in certain diffusion flame problems to be outlined in Section 7.2(b). The former method does not suffer from this disadvantage. 

The partial equilibrium assumptions by themselves in conjunction with Eqs. (9.1), (9.2), and (9.3), and a reaction mechanism as outlined above, do not permit construction of a complete model from which an eigenvalue burning velocity and full profiles may be computed ab initio. On the other hand the assumptions are extremely useful when dealing with H-N-C-O ffame systems, since their application to reactions (i), (ii), (iii), and (xviii) above allows us to calculate many of the species profiles, and particularly the free radical profiles, on close approach to full equilibrium. The time-dependent computation does not economically do this directly. The computations require an input mass flux or burning velocity which must be either a measured or a separately calculated value. For composite flux calculation purposes the overall radical pool is chosen so as to represent a total flux of free electron spins, that is, spins belonging to H, O, OH, and O2 (Dixon-Lewis et al, 1975). 

In the present work we attempt to gain a better understanding of the destruction mechanism of dibenzofuranyl + O2 system as well as the reaction pathways important in its oxygen-free decomposition. The difficulties encountered now are related to the large size of aromatic species contained in the dibenzofuranyl + O2 system making high level calculations very costly or not possible. The approach we have taken to circumvent the problem is to reduce the system to the smallest representative unit and then to proceed with the computation of the thermochemical properties (Figure 1.1). Based on this approach dibenzofuranyl peroxy radical (A) can be concentrated around the phenyl peroxy radical reaction system (B), which itself can be reduced to vinyl peroxy radical system (C). 

---