# The Chain Length Distribution of an Ideal Reversible Deactivation Radical Polymerization

Simon Harrisson

► To cite this version:

Simon Harrisson. The Chain Length Distribution of an Ideal Reversible Deactivation Radical Polymerization. Polymers, 2018, Polymerization Kinetics, 10 (8), pp.887. 10.3390/polym10080887. hal-02194890

HAL Id: hal-02194890

https://hal.science/hal-02194890v1

Submitted on 1 Jul 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# The Chain Length Distribution of an Ideal Reversible Deactivation Radical Polymerization

Simon Harrisson ID

Laboratoire des IMRCP, Université de Toulouse, CNRS UMR 5623, Université Paul Sabatier, 118 route de Narbonne, 31062 Toulouse CEDEX 9, France; harrisson@chimie.ups-tlse.fr; Tel.: +33-5-6155-6968

Received: 18 July 2018; Accepted: 4 August 2018; Published: 8 August 2018

<figure><img src="image_2.png" /></figure>

check for
updates

**Abstract:** The chain length distribution (CLD) of a reversible deactivation radical polymerization at full conversion is shown to be a negative binomial distribution with parameters that are simple functions of the number-average degree of polymerization and either the chain transfer constant (in the case of polymerizations that incorporate a reversible chain transfer step) or the concentrations of dormant polymer chains and deactivating agent and the rate constants of propagation and deactivation (other types of RDRP). Expressions for the CLD at intermediate conversions are also derived, and shown to be consistent with known expressions for the number-average degree of polymerization and dispersity. It is further demonstrated that these CLDs are well-approximated by negative binomial distributions with appropriate choice of parameters. The negative binomial distribution is thus a useful model for CLDs of reversible deactivation radical polymerizations.

**Keywords:** polymer; reversible deactivation radical polymerization; chain length distribution; negative binomial distribution

## 1. Introduction

Reversible deactivation radical polymerization (RDRP) encompasses a range of polymerization techniques developed over the last 25 years that permit the synthesis of well-defined, relatively narrowly-dispersed polymers via a radical process. The fundamental feature of RDRP is that dormant polymer chains are reversibly activated to form propagating radicals which can grow by adding one or more monomer units before returning to the dormant state (Scheme 1) [1].

<figure><img src="image_3.png" /><figcaption>**Scheme 1.** General mechanism of the reversible deactivation step in reversible deactivation radical polymerization.</figcaption></figure>
RDRP techniques are widely used to generate highly functional materials with complex architectures (e.g., block, star, comb, etc.). In particular, several recent studies have applied RDRP to the preparation of highly complex structures containing large numbers of functional blocks [2–9].

The true precision of such complex structures is difficult to quantify. One approach is to assume near-perfect control of the polymerization, with the independent addition of monomers to each chain [10]. This is analogous to an anionic polymerization, and results in a Poisson distribution of chain lengths, with dispersity, D, given by Equation (1) [11].

$$D = 1 + 1/DP_n \qquad (1)$$

In reality, monomers do not add independently of one another, as multiple monomers may be added during a single activation-deactivation cycle (with rate constant $k_p$). Assumption of a Poisson chain length distribution (CLD) will thus lead to an underestimation of the true degree of structural variation. It would be useful to have a simple model that more accurately reflects the CLD of polymers produced by RDRP techniques.

Early work on RDRP kinetics led to the derivation of a formula for the dispersity of polymers produced by RDRP processes as a function of conversion, c (Equation (2)) [12–14].

$$D = 1 + 1/DP_n + \left(\frac{2-c}{c}\right) \left(\frac{k_p[PS]}{k_d[S^*]}\right) \qquad (2)$$

In Equation (2), [S*] is the concentration of deactivating agent (assumed to be constant with rate constant of deactivation $k_d$). This could represent nitroxide radicals in nitroxide-mediated radical polymerization (NMRP, S* represented as S* in Scheme 1) [15], copper(II) species in atom transfer radical polymerization (ATRP, S* represented as MS in Scheme 1) [16], or reversible addition-fragmentation chain transfer agents (RAFT agents) in RAFT polymerization (in this case S* and PS are identical) [17].

In many forms of RDRP, the deactivating agent is derived from the dormant polymer (e.g., an alkyl halide in ATRP or alkoxyamine in NMRP), designated as PS in Equation (2) and Scheme 1. In reversible chain transfer polymerizations such as RAFT polymerization, the deactivating agent and the dormant polymer are identical. If the ratio of [S*] to [PS] is assumed to be constant (equal to 1 for reversible chain transfer polymerizations), an apparent chain transfer constant, $C_S$, can be defined (Equation (3)).

$$C_S = \frac{k_d[S^*]}{k_p[PS]} \qquad (3)$$

Substitution into Equation (2) gives Equation (4) [12,13].

$$D = 1 + 1/DP_n + \left(\frac{2-c}{c}\right) \cdot 1/C_S \qquad (4)$$

This simple formula has been extended by Zhu et al. to account for side reactions such as irreversible termination, which cause broadening of the CLD as the reaction proceeds to high conversion [18].

While numerous expressions exist for the dispersity under more or less realistic assumptions, there are relatively few examples of expressions for the entire CLD [19]. Gold proposed a modified Poisson distribution which takes into account unequal initiation and propagation rates, but assumes the independent addition of monomers [20]. Muller et al. proposed an expression for the CLD of an RDRP process with slow equilibration between active and dormant species, but assumed a constant monomer concentration [21]. Tobita [22,23] has derived an expression for the CLD of an RDRP continuous stirred tank reactor (constant monomer concentration) process as a hypergeometric function combining the Poisson and geometric (most probable) distributions. Finally, Konkolewicz et al. have developed a recursive equation for the CLD of a batch RDRP by convoluting Poisson and geometric distributions
while taking into account the changing monomer concentration in the batch process [24–26]. The resulting CLD is given by Equation (5) [24].

$$ \Pr(DP = k) = \sum_{j=0}^{\infty} \frac{\mu_{\text{decay}}^{j} e^{-\mu_{\text{decay}}} j!}{j!} \Pr(DP = k | N = j) \quad (5) $$

In Equation (5), $\mu_{\text{decay}}$ is the mean number of activation cycles that the chains have undergone, while $\Pr(DP = k | N = j)$ represents the probability that a chain that has been activated $j$ times has a degree of polymerization of $k$. This probability is given by the recursive formula of Equation (6).

$$ \Pr(DP = k | N = j) = \sum_{i=0}^{m} \Pr(DP = k - i | N = j - 1) \frac{[M]_{0}}{C_{S}[S^{*}]} \left(1 - \frac{1}{C_{S}}\right)^{j-1} \Pr(DP = k | N = 0) = \begin{cases} 1, & k = 0 \\ 0, & k > 0 \end{cases} \quad (6) $$

In this article, a similar approach to that of Tobita and Konkolewicz et al. is used to show that the CLD of a polymer produced by an ideal RDRP batch process at full conversion is given by the negative binomial distribution with parameters $C_S$ and $DP_n/(DP_n + C_S)$ (Equation (7)).

$$ \Pr(DP = k) = \frac{\Gamma(k + C_S)}{k! \Gamma(C_S)} \left( \frac{C_S}{DP_n + C_S} \right)^{C_S} \left( \frac{DP_n}{DP_n + C_S} \right)^k \quad (7) $$

The distribution at intermediate conversions is also derived, as well as parameters for a negative binomial distribution that closely approximates the true distribution. First derived in 1714 [27], the negative binomial distribution is one of the fundamental distributions of statistics, and is frequently used to model discrete data whose variance is greater than their mean [28]. The considerable body of statistical literature on this distribution could thus be applied to the study of polymer CLDs.

## 2. Materials and Methods

All simulations were carried out using Excel 2013 software on a desktop computer. The code used to generate simulated distributions of 10,000 chains is given in the Supporting Information.

## 3. Results

We define an ideal RDRP as one in which all rate constants are independent of chain length and conversion, and no termination or other side reactions take place.

The rate of propagation is given by Equation (8).

$$ R_p = - \frac{d[M]}{dt} = k_p [M] [P^*] \quad (8) $$

In this equation, $[M]$ represents the monomer concentration, and $[P^*]$, the concentration of active chains (assumed to be constant).

In terms of conversion, $c$:

$$ R_p = [M]_0 \frac{dc}{dt} = k_p [M]_0 (1-c) [P^*] \quad (9) $$

$$ c = 1 - e^{-k_p [P^*] t} \quad (10) $$

Chain deactivation occurs at a constant rate, $R_d$, given by Equation (11).

$$ R_d = k_d [S^*] [P^*] \quad (11) $$
The probability, *p*, that an active chain will propagate is given by

$$p = \frac{R_p}{R_p + R_d} = \frac{k_p[M]}{k_p[M] + k_d[S^*]} = \frac{[M]/[PS]}{[M]/[PS] + k_d[S^*]/k_p[PS]} \quad (12)$$

Defining *T* as the target degree of polymerization at full conversion, equal to [M]₀/[PS], we have

$$p = \frac{T(1-c)}{T(1-c) + C_S} \quad (13)$$

In a single activation-deactivation cycle, an active chain will add *n* monomers, where *n* is a geometrically distributed random variable [22–26]. The probability that *k* monomers are added in a single activation/deactivation cycle at time *t* is given by

$$\mathrm{Pr}(n = k)_t = p^k(1-p) \quad (14)$$

As the concentration of monomer changes during the polymerization, the instantaneous distribution of *n* also varies. The total distribution of *n*, representing all activation/deactivation cycles from the beginning of the polymerization (*t* = 0) to a time *t* corresponding to conversion *c*, is obtained by integrating with respect to time:

$$\mathrm{Pr}(n = k) = \frac{\int_{0}^{t} \left( \frac{T(1-c)}{T(1-c) + C_S} \right)^k \left( \frac{C_S}{T(1-c) + C_S} \right) dt}{\int_{0}^{t} dt} \quad (15)$$

In terms of conversion, this is

$$\mathrm{Pr}(n = k) = \frac{\int_{0}^{c} \left( \frac{T(1-c)}{T(1-c) + C_S} \right)^k \left( \frac{C_S}{T(1-c) + C_S} \right) \frac{dc}{1-c}}{-\ln(1-c)} \quad (16)$$

Evaluation of the integral gives

$$\mathrm{Pr}(n = k) = \begin{cases} 1 + \frac{\ln\left(\frac{1-p_f}{1-p_0}\right)}{\ln(1-c)}, & k = 0 \\ \frac{(p_0)^k - (p_f)^k}{-k \ln(1-c)}, & k \ge 1 \end{cases} \quad (17)$$

where

$$p_0 = \frac{T}{T + C_S}, \quad p_f = \frac{T(1-c)}{T(1-c) + C_S}$$

This distribution has a probability generating function [29], $G_n(z)$, given by Equation (18):

$$G_n(z) = \sum_{k=0}^{\infty} z^k \mathrm{Pr}(n = k) = 1 + \frac{\ln\left[\frac{(1-p_0 z)(1-p_f)}{(1-p_0)(1-p_f z)}\right]}{\ln(1-c)} \quad (18)$$

Throughout the polymerization, activation/deactivation cycles occur at a constant rate. The number, *N*, of activation/deactivation cycles experienced by each chain is Poisson distributed, with an expected number of cycles per chain equal to the total number of cycles divided by the total number of chains [22–26].

$$E(N) = \frac{k_d[S^*][P^*]t}{[PS]} = -C_S \ln(1-c) \quad (19)$$

The corresponding probability generating function $G_N(z)$ is given by Equation (20):

$$G_N(z) = e^{-C_S \ln(1-c)(z-1)} \quad (20)$$
Each chain consists of N segments of length $n_1, n_2, \dots, n_N$, where $n_i$ are independent, identically distributed random variables with probability generating function $G_n(z)$ and N is Poisson-distributed with expected value of $-C_S \ln(1-c)$.

The distribution of total chain lengths, $Y = n_1 + n_2 + \dots + n_N$, is thus a compound Poisson distribution with probability generating function $G_Y(z)$ given by Equation (21):

$$G_Y(z) = G_N(G_n(z)) = \left[ \frac{(1-p_0)(1-p_f z)}{(1-p_0 z)(1-p_f)} \right]^{C_S} \quad (21)$$

The probability mass function of the CLD can be obtained from the derivatives of $G_Y(z)$ evaluated at $z = 0$, leading to the recursive formula of Equation (22):

$$\mathrm{Pr}(Y = k) = \frac{G_Y^{(k)}(0)}{k!} = \begin{cases} \left(\frac{1-p_0}{1-p_f}\right)^{C_S} = \left[\frac{T(1-c)+C_S}{T+C_S}\right]^{C_S}, & k=0 \\ \frac{C_S}{k} \sum_{i=1}^{k} ((p_0)^i - (p_f)^i) \mathrm{Pr}(Y = k-i), & k \ge 1 \end{cases} \quad (22)$$

Figure 1 shows the comparison of calculated CLDs obtained from Equation (22) and simulated CLDs for $C_S$ of 0.5 and 5 and conversions of 25–100%. Additional comparisons of calculated and simulated CLDs for a wider range of chain transfer constants may be found in the Supporting Information (Figures S1–S4).

The expected value and dispersity of the CLD can also be obtained from $G_Y(z)$, and give

$$DP_n = \lim_{z \to 1^{-}} G'_{Y}(z) = cT \quad (23)$$

$$\mathcal{D} = \lim_{z \to 1^{-}} \left( \frac{G''_{Y}(z) + G'_{Y}(z)}{[G'_{Y}(z)]^2} \right) = 1 + \frac{1}{DP_n} + \left( \frac{2-c}{c} \right) \cdot 1/C_S \quad (24)$$

Equation (24) is identical to the expression derived by Müller et al. [13] (Equation (4)) using the method of moments [30].

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Comparison of calculated CLDs (solid lines) with simulated distributions (open circles, 10,000 chains) for polymerizations with $T = 100$, $C_S = 0.5$ (a) or 5 (b), conversion = 0.25, 0.5, 0.75 and 1.</figcaption></figure>

The above expressions include chains of length 0, which would not normally be included in the determination of molecular weight. Corrected values for $DP_n$ and $\mathcal{D}$ can be obtained by excluding chains of zero length [13] to give:

$$DP_n^{k \neq 0} = \frac{cT}{\left(1 - \left[\frac{T(1-c)+C_S}{T+C_S}\right]^{C_S}\right)} \quad (25)$$
$$ \mathcal{D}^{k \neq 0} = \left(1 + \frac{1}{cT} + \frac{2-c}{c \cdot C_S}\right) \left(1 - \left[\frac{T(1-c) + C_S}{T + C_S}\right]^{C_S}\right) \quad (26) $$

As the conversion approaches 1, $p_f$ approaches 0, and $G_Y(z)$ approaches:

$$ \lim_{X \to 1} G_Y(z) = \left[ \frac{(1 - p_0)}{(1 - p_0 z)} \right]^{C_S} \quad (27) $$

This is the probability generating function of a negative binomial distribution with parameters $C_S$ and $p_0$. Thus, the CLD of a polymer prepared under ideal RDRP batch conditions at full conversion is a negative binomial distribution, with probability mass function given by Equation (7), number average chain length equal to $T$, and dispersity given by Equation (28).

$$ \mathcal{D} = 1 + 1/DP_n + 1/C_S \quad (28) $$

## 4. Discussion

While the negative binomial distribution strictly describes the CLD only at full conversion, it also provides a good approximation to CLDs at intermediate conversions. The full CLD (including zero-length chains) is approximated by the random variable $Z$ with probability mass function given by Equation (29), where $Z$ is a mixture of a negative binomial distribution with parameters $\rho$ and $\pi$ and a constant random variable representing unreacted chains:

$$ \begin{aligned} \Pr(Z = k) &= \begin{cases} (1-c)^{C_S} + (1-(1-c)^{C_S})(1-\pi)^{\rho}, & k = 0 \\ (1-(1-c)^{C_S}) \left( \frac{k+\rho-1}{k} \right) (1-\pi)^{\rho}\pi^k, & k \ge 1 \end{cases} \\ \rho &= \left( \left( 1 + \frac{2-c}{c \cdot C_S} \right) (1-(1-c)^{C_S}) - 1 \right)^{-1} \\ \pi &= 1 - \left( 1 + cT + \frac{T(2-c)}{C_S} - \frac{cT}{1-(1-c)^{C_S}} \right)^{-1} \end{aligned} \quad (29) $$

As shown in Figure 2, the approximation is best for low conversions (when the CLD is approximately geometric), high conversions (when the CLD is approximately negative binomial) and $C_S$ close to 1 (when the CLD is approximately geometric) or $\gtrsim 50$ (when the CLD is approximately Poisson). A comparison between approximate and calculated distributions for a wide range of $C_S$ and $T$ is given in the Supporting Information (Figures S5–S8). The quality of the approximation is comparable to that of the commonly-used [31–34] Schulz-Zimm distribution [35,36] (Equation (30)).

$$ \begin{aligned} n(k) &= \frac{k^{z-1}y^z e^{-yk}}{\Gamma(z)} \\ z &= \frac{1}{\mathcal{D}-1}, \quad y = \frac{z}{DP_n} \end{aligned} \quad (30) $$

Comparisons of the negative binomial and Schulz-Zimm distributions are presented in the Supporting Information, Table S1, Figures S9–S12.

The results obtained in this work assume an initial degree of polymerization of 0, and are thus valid for chains grown from a small molecule initiator or chain transfer agent, or individual blocks of a block copolymer. The total CLD of a block copolymer is given by the sum of the CLDs of its individual blocks. This sum is not a negative binomial distribution unless all blocks were polymerized to full conversion and $DP_n/(DP_n + C_S)$ is the same for each block. In real polymerizations, the presence of side reactions such as termination and irreversible chain transfer may result in significant deviations from the ideal case considered here [18]. However, when relatively short chains are targeted, such reactions may be disregarded [6,7]. Thus the negative binomial distribution may usefully serve as a model for the CLD of short blocks within multiblock copolymers prepared by RDRP [4–9]. Additionally, the simplicity of the negative binomial distribution and its implementation in readily-available software
such as Excel™ facilitates the visualization of CLDs as an aid to teaching or comprehension. This allows features such as the number of unreacted chains or the relative proportions of high and low molecular weight polymer to be estimated, and provides a more nuanced picture of the overall CLD than the dispersity alone.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Comparison of negative binomial approximation (thick red lines) with calculated CLDs (thin black lines) for polymerizations with *T* = 100, *C*<sub>S</sub> = 0.5 (a) or 5 (b), conversion = 0.25, 0.5, 0.75 and 1.</figcaption></figure>

**5. Conclusions**

The chain length distribution of an ideal RDRP at full conversion is given by a negative binomial distribution with parameters C<sub>S</sub> and DP<sub>n</sub>/(DP<sub>n</sub> + C<sub>S</sub>). Unlike previously derived expressions for the CLD of an RDRP process, this distribution can be expressed as a simple function of two fundamental parameters of the polymerization. This simplicity is achieved, however, at the cost of assumptions such as the complete absence of termination, which will not hold in real polymerizations. Nevertheless, such a simple model will be useful in situations where side-reactions may be assumed to be negligible, and complements more complex models such as those of Tobita [22,23] and Konkolewicz et al. [24–26] that deal with more realistic systems that incorporate irreversible termination and transfer reactions.

**Supplementary Materials:** The following are available online at http://www.mdpi.com/2073-4360/10/8/887/s1. Simulated CLDs for a wide range of C<sub>S</sub> and targeted degrees of polymerization, accompanied by calculated CLDs and negative binomial approximations. Comparison of goodness of fit of negative binomial approximation and Schulz-Zimm distribution to calculated CLDs. Details of procedure used to simulate the CLD.

**Funding:** This work was funded by the Agence Nationale de la Recherche, grant number ANR-15-CE08-0039.

**Acknowledgments:** S.H. acknowledges funding from the Agence Nationale de la Recherche (Project ASYMCOPO, ANR-15-CE08-0039).

**Conflicts of Interest:** The authors declare no conflict of interest.

**References**

1. Jenkins, A.D.; Jones, R.G.; Moad, G. Terminology for reversible-deactivation radical polymerization previously called “controlled” radical or “living” radical polymerization (IUPAC Recommendations 2010). *Pure Appl. Chem.* **2010**, *82*, 483–491. [CrossRef]

2. Zamfir, M.; Lutz, J.-F. Ultra-precise insertion of functional monomers in chain-growth polymerizations. *Nat. Commun.* **2012**, *3*, 1138. [CrossRef]

3. Lutz, J.-F.; Schmidt, B.V.K.J.; Pfeifer, S. Tailored Polymer Microstructures Prepared by Atom Transfer Radical Copolymerization of Styrene and N-substituted Maleimides. *Macromol. Rapid Commun.* **2011**, *32*, 127–135. [CrossRef] [PubMed]

4. Anastasaki, A.; Waldron, C.; Wilson, P.; Boyer, C.; Zetterlund, P.B.; Whittaker, M.R.; Haddleton, D. High Molecular Weight Block Copolymers by Sequential Monomer Addition via Cu(0)-Mediated Living Radical Polymerization (SET-LRP): An Optimized Approach. *ACS Macro Lett.* **2013**, *2*, 896–900. [CrossRef]
5. Gody, G.; Maschmeyer, T.; Zetterlund, P.B.; Perrier, S. Rapid and quantitative one-pot synthesis of sequence-controlled polymers by radical polymerization. *Nat. Commun.* **2013**, *4*, 2505. [CrossRef] [PubMed]

6. Gody, G.; Maschmeyer, T.; Zetterlund, P.B.; Perrier, S. Pushing the Limit of the RAFT Process: Multiblock Copolymers by One-Pot Rapid Multiple Chain Extensions at Full Monomer Conversion. *Macromolecules* **2014**, *47*, 3451–3460. [CrossRef]

7. Gody, G.; Maschmeyer, T.; Zetterlund, P.B.; Perrier, S. Exploitation of the Degenerative Transfer Mechanism in RAFT Polymerization for Synthesis of Polymer of High Livingness at Full Monomer Conversion. *Macromolecules* **2014**, *47*, 639–649. [CrossRef]

8. Engelis, N.G.; Anastasaki, A.; Nurumbetov, G.; Truong, N.P.; Nikolaou, V.; Shegiwal, A.; Whittaker, M.R.; Davis, T.P.; Haddleton, D.M. Sequence-controlled methacrylic multiblock copolymers via sulfur-free RAFT emulsion polymerization. *Nat. Chem.* **2017**, *9*, 171–178. [CrossRef] [PubMed]

9. Wenn, B.; Martens, A.C.; Chuang, Y.M.; Gruber, J.; Junkers, T. Efficient multiblock star polymer synthesis from photo-induced copper-mediated polymerization with up to 21 arms. *Polym. Chem.* **2016**, *7*, 2720–2727. [CrossRef]

10. Gody, G.; Zetterlund, P.B.; Perrier, S.; Harrisson, S. The limits of precision monomer placement in chain growth polymerization. *Nat. Commun.* **2016**, *7*. [CrossRef] [PubMed]

11. Flory, P.J. Molecular Size Distribution in Ethylene Oxide Polymers. *J. Am. Chem. Soc.* **1940**, *62*, 1561–1565. [CrossRef]

12. Goto, A.; Fukuda, T. Kinetics of living radical polymerization. *Prog. Polym. Sci.* **2004**, *29*, 329–385. [CrossRef]

13. Mueller, A.H.E.; Zhuang, R.; Yan, D.; Litvinenko, G. Kinetic Analysis of “Living” Polymerization Processes Exhibiting Slow Equilibria. 1. Degenerative Transfer (Direct Activity Exchange between Active and “Dormant” Species). Application to Group Transfer Polymerization. *Macromolecules* **1995**, *28*, 4326–4333. [CrossRef]

14. Müller, A.H.E.; Litvinenko, G.; Yan, D. Kinetic Analysis of “Living” Polymerization Systems Exhibiting Slow Equilibria. 4. “Dissociative” Mechanism of Group Transfer Polymerization and Generation of Free Ions in Cationic Polymerization. *Macromolecules* **1996**, *29*, 2346–2353. [CrossRef]

15. Nicolas, J.; Guillaneuf, Y.; Lefay, C.; Bertin, D.; Gigmes, D.; Charleux, B. Nitroxide-mediated polymerization. *Prog. Polym. Sci.* **2013**, *38*, 63–235. [CrossRef]

16. Matyjaszewski, K.; Xia, J. Atom Transfer Radical Polymerization. *Chem. Rev.* **2001**, *101*, 2921–2990. [CrossRef] [PubMed]

17. Moad, G.; Rizzardo, E.; Thang, S.H. Living Radical Polymerization by the RAFT Process. *Aust. J. Chem.* **2005**, *58*, 379–410. [CrossRef]

18. Mastan, E.; Zhu, S. A Molecular Weight Distribution Polydispersity Equation for the ATRP System: Quantifying the Effect of Radical Termination. *Macromolecules* **2015**, *48*, 6440–6449. [CrossRef]

19. Mastan, E.; Li, X.; Zhu, S. Modeling and theoretical development in controlled radical polymerization. *Prog. Polym. Sci.* **2015**, *45*, 71–101. [CrossRef]

20. Gold, L. Statistics of Polymer Molecular Size Distribution for an Invariant Number of Propagating Chains. *J. Chem. Phys.* **1958**, *28*, 91–99. [CrossRef]

21. Mueller, A.H.E.; Yan, D.; Litvinenko, G.; Zhuang, R.; Dong, H. Kinetic Analysis of “Living” Polymerization Processes Exhibiting Slow Equilibria. 2. Molecular Weight Distribution for Degenerative Transfer (Direct Activity Exchange between Active and “Dormant” Species) at Constant Monomer Concentration. *Macromolecules* **1995**, *28*, 7335–7338. [CrossRef]

22. Tobita, H. Molecular Weight Distribution of Living Radical Polymers. *Macromol. Theory Simul.* **2006**, *15*, 12–22. [CrossRef]

23. Tobita, H. Fundamental Molecular Weight Distribution of RAFT Polymers. *Macromol. React. Eng.* **2008**, *2*, 371–381. [CrossRef]

24. Konkolewicz, D.; Hawkett, B.S.; Gray-Weale, A.; Perrier, S. RAFT Polymerization Kinetics: Combination of Apparently Conflicting Models. *Macromolecules* **2008**, *41*, 6400–6412. [CrossRef]

25. Konkolewicz, D.; Hawkett, B.S.; Gray-Weale, A.; Perrier, S. RAFT polymerization kinetics: How long are the cross-terminating oligomers? *J. Polym. Sci. A Polym. Chem.* **2009**, *47*, 3455–3466. [CrossRef]

26. Konkolewicz, D.; Siauw, M.; Gray-Weale, A.; Hawkett, B.S.; Perrier, S. Obtaining Kinetic Information from the Chain-Length Distribution of Polymers Produced by RAFT. *J. Phys. Chem. B* **2009**, *113*, 7086–7094. [CrossRef] [PubMed]
27. Bartko, J.J. A Note on the Negative Binomial Distribution. *Technometrics* **1962**, *4*, 609–610. [CrossRef]

28. Hilbe, J.M. *Negative Binomial Regression*, 2nd ed.; Cambridge University Press: Cambridge, UK, 2011.

29. Grimmett, G.R.; Stirzaker, D.R. *Probability and Random Processes*, 3rd ed.; Oxford University Press: Oxford, UK, 2001; pp. 148–155.

30. Mastan, E.; Zhu, S. Method of moments: A versatile tool for deterministic modeling of polymerization kinetics. *Eur. Polym. J.* **2015**, *68*, 139–160. [CrossRef]

31. Qi, S.; Klushin, L.I.; Skvortsov, A.M.; Schmid, F. Polydisperse Polymer Brushes: Internal Structure, Critical Behavior, and Interaction with Flow. *Macromolecules* **2016**, *49*, 9665–9683. [CrossRef]

32. Patil, R.R.; Turgman-Cohen, S.; Šrogl, J.; Kiserow, D.; Genzer, J. On-Demand Degrafting and the Study of Molecular Weight and Grafting Density of Poly(methyl methacrylate) Brushes on Flat Silica Substrates. *Langmuir* **2015**, *31*, 2372–2381. [CrossRef] [PubMed]

33. Schaeffel, D.; Yordanov, S.; Staff, R.H.; Kreyes, A.; Zhao, Y.; Schmidt, M.; Landfester, K.; Hofkens, J.; Butt, H.-J.; Crespy, D.; Koynov, K. Fluorescence Correlation Spectroscopy in Dilute Polymer Solutions: Effects of Molar Mass Dispersity and the Type of Fluorescent Labeling. *ACS Macro Lett.* **2015**, *4*, 171–176. [CrossRef]

34. Lynd, N.A.; Meuler, A.J.; Hillmyer, M.A. Polydispersity and block copolymer self-assembly. *Prog. Polym. Sci.* **2008**, *33*, 875–893. [CrossRef]

35. Schulz, G.V. Über die Kinetik der Kettenpolymerisationen. *V. Z. Phys. Chem.* **1939**, *B43*, 25–46. [CrossRef]

36. Zimm, B.H. Apparatus and Methods for Measurement and Interpretation of the Angular Variation of Light Scattering; Preliminary Results on Polystyrene Solutions. *J. Chem. Phys.* **1948**, *16*, 1099–1117. [CrossRef]

<figure><img src="image_6.png" /></figure>

© 2018 by the author. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).