# A quantum chemistry study of the free-radical copolymerization propagation kinetics of styrene and 2-hydroxyethyl acrylate.

INTRODUCTIONSeveral works focused on the study of acrylic copolymers, widely used as resins for coatings, appeared recently in the literature. In particular, the free radical copolymerization kinetics was investigated to get a good understanding of the mechanism of copolymer chain-growth (1-7). Copolymers based on 2-hydroxyethyl acry-late (HEA) were reported by copolymerization with ethers, acrylates, or methacrylates (8-13). Notably, HEA is a functional monomer whose copolymers are also useful as biomaterials (14-17). In fact, even if acrylates are generally more toxic than methacrylates (18), (19), poly-HEA (PHEA) based materials are used as cell-culturing media and are applicable for in vitro experiments (20).

Despite of the PHEA potentials, a few studies on the kinetics of copolymerization involving HEA have been reported. The direct synthesis of PHEA is made through atom transfer radical polymerization as controlled free radical polymerization method (21), whereas Bian and Cunningham proposed the nitroxide-mediated living radical polymerization of HEA (22). Jansen et al. investigated the copolymerization of HEA and 2-hydroxyethyl methac-rylate using the real-time Fourier transform infrared spectroscopy (10). With the same technique, it was found that the copolymerization of HEA and n-butyl acrylate, by the use of a nitroxide mediator, leads to an acceleration of HEA conversion and polymerization rate, thanks to hydrogen bonding interactions between the hydroxyl group of HEA monomer and the nitroxide (9). Moreover, the influence of the solvent on the apparent reactivity ratios in free radical copolymerization between HEA and itaconic acid was also studied (23). With particular reference to the HEA/styrene (ST) system, solution copolymerizations were carried out in benzene at 60[degrees]C (24), whereas proton nuclear magnetic resonance and gel permeation chromatography were used to study HEA/ST free radical copolymerization mediated by 2,2,6,6-tetra-methyl-1-piperidinyloxy (TEMPO) at 125 [degrees] C (25). Finally, Monte Carlo simulations were used to predict molecular weight distribution and composition of HEA/ST copolymers (26).

The limited number of investigations on HEA reactivity in free radical polymerization processes is partly due to the fact that HEA polymerization leads to high molecular weight products through crosslinking and transfer to polymer reactions. As a matter of fact in recent years, it was well established that the main issue in acrylate kinetics is the occurrence of secondary reactions, such as chain transfer and the formation of midchain radicals (MCRs) (27-29). The experimental evaluation of kinetic parameters for these reactions is very complicated to achieve due to the reactivity differences between secondary and tertiary carbon radicals (30). However, improvements have been recently made for several polymerization systems applying the pulsed laser polymerization (PLP) technique (28), (30-32). Although previous mentioned studies on this topic did not refer directly to HEA, the same reaction patterns might occur even for this monomer. Therefore, the resulting complexity of the operative kinetic scheme is the main cause of difficulties in the experimental analysis (33), (34).

To overcome such problems, and with the aim of supporting and enriching the experimental efforts, theoretical predictions of kinetics based on quantum mechanics methods could represent an effective tool (35-45). This is especially promising since the potential of density functional theory (DFT) to investigate free radical copolymer-ization kinetics was already shown in the literature (6), (7), (40), (46). In this work, a DFT based computational approach is applied to study the reactivity of HEA in free radical homopolymerization and copolymerization with ST. In particular, the monomer reactivity ratios, [r.sub.UEA] and [r.sub.ST], are evaluated at 50 [degrees] C and 125 [degrees] C and compared with those proposed in literature (8), (25). Finally, the evolution of the copolymer composition is investigated according to the terminal model (47).

COMPUTATIONAL DETAILS

The selected computational approach adopted is based on DFT. In particular, the Becke 3 parameter and Lee-Yang-Parr functional (B3LYP) were used in the DFT calculations to evaluate exchange and correlation energies (48), (49). All quantum chemical calculations of reactants and products were performed with a spin multiplicity of 2 using an unrestricted wave function (UB3LYP). B3LYP methods provide excellent low-cost performance, as demonstrated in previous works reported in the literature (6), (7), (39), (45), (50-52). The all electron 6-31 basis set with added polarization functions (6-31G(d,p)) was used as basis set. All geometries were fully optimized with the Berny algorithm and were followed by frequency calculations. The geometry of each molecular structure was considered stable only after calculating vibrational frequencies and force constants and if no imaginary vibrational frequency was found. Transition state structures are located adopting the synchronous transit-guided quasi Newton method and are characterized by a single imaginary vibrational frequency [53]. The corresponding kinetic constants were determined through the classical transition state theory as:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

where [k.sub.b] and h are Boltzmann and Plank constant, respectively, T is the temperature, [E.sub.a] is the activation energy of the process, and Q represents the product of the partition functions ([q.sup.vib] [q.sup.rot] [q.sup.(rot, int)]) for transition states ([not equal to]) and reactants (R). In particular, [q.sup.vib] and [q.sup.rot] are the vibrational and rotational partition functions, respectively, calculated according to Eqs. 2 and 3.

where [v.sub.i] are the vibrational frequencies, [I.sub.x], [I.sub.y], [I.sub.z] are the rotational constants, and [sigma] is the rotational symmetry number. Moreover, [q.sup.(rot, int)] is the internal rotation partition function. In fact, previous works showed that the computational accuracy is improved by treating the low-vibralional frequencies, indicative of a small energy barrier for the relative motion of a partion of the molecule with respect to the other one, as internal rotations [35, 36, 54, 55]. More in detail, for internal motions corresponding to the lower vibrational frequencies, the potential energy was calculated, as a function of the rotational angle, at the B3LYP/6-31G(d) level of theory, and it was then interpolated using cubic splines. As reported in the literature [51], the calculated potential energy, V ([empty set]), was used to solve the rotational ID Schrodinger equation:

- [h.sup.2] / 2 [I.sub.m] x [a.sub.2] [PHI] / a [[empty set].sup.2] + V ([empty set]) x [PHI] ([empty set]) = [epsilon] x [PHI] ([epsilon])

where h is the ratio between the Planck constant and 2 [pi], [PHI] is the energy levels associated with each vibrational frequency, [PHI] is the rotational angle, and [I.sub.m] is the reduced moment of inertia of a part of the molecule with respect to the other one calculated either as (56):

[I.sub.m] = [I.sub.1] x [I.sub.2] / [I.sub.1] + [I.sub.2]

where [I.sub.1] and [I.sub.2] are the moments of inertia of the two rotating moieties. The internal rotation partition function was then calculated according to Eq. 6.

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (6)

where zk is the calculated eigenvalue, [g.sub.k] is the degeneracy of rotational energy level [[epsilon].sub.k], and [[sigma].sub.int] is the symmetry number of the internal rotation.

All quantum chemistry calculations are performed with the Gaussian 03 suite of programs and all pictures are drawn with PyMol 1.3 [57, 58].

RESULTS AND DISCUSSION

As mentioned in the Introduction section, there are few reported studies on HEA/ST free radical copolymerization; moreover, data reported in the literature are quite scattered. According to the Q-e scheme proposed by Alfrey and Price [59], the monomer reactivity ratios are estimated as [r.sub.HEA] = 0.31 and [r.sub.ST] = 0.36 [26]. Chow measured the same parameters by solution polymerization of HEA with ST in benzene at 60 [degrees] C obtaining values of [r.sub.HEA] = 0-34 and [r.sub.ST] = 0.38 [24]. Later, Penlidis and coworkers proposed a series of observations about the HEA/ST polymerization, examining reactivity ratios and full conversion kinetics as a function of experimental variables such as feed composition, temperature, and initiation concentration. Copolymer composition data at 50 [degrees] C were proposed and reactivity ratios of [r.sub.HEA] = 0.20 and [r.sub.ST] = 0.46 were determined. The work emphasized all the difficulties related to the experimental investigation of HEA/ST system, showing in particular the variability in the reactivity ratio values probably due to the polarity change in the copolymer and monomer mixture, depending on the initial comonomer feed composition [8]. Recently, the living free-radical copolymerizations of ST with a series of polar monomers, including HEA, mediated by TEMPO and initiated by benzoyl peroxide at 125 [degrees] C were studied [25]. Monomer reactivity ratios of [r.sub.HEA] = 0.43 and [r.sub.ST] = 0.28 were found, but limits in controlling the process during the polymerization are highlighted even in this work.

HEA* + HEA* [right arrow] [k.sub.HEA - HEA] HEA - HEA*

HEA* + ST [right arrow] [k.sub.HEA - ST] HEA - ST*

ST* + ST [right arrow] [k.sub.pST - ST] SY - ST*

ST* + HEA [right arrow] [k.sub.ST - HEA] ST - HEA

SCHEME 1. Elementary reactions involved in the determination of the monomer reactivity ratios.

With the aim to compare our results with the copolymer composition data available in literature [8, 25], the HEA/ST copolymerization was studied by the computational method at 50 [degrees] C and 125 [degrees] C. To obtain rHEA and rST values, the elementary reactions in Scheme 1 were simulated. Monomer reactivity ratios were determined as defined in Eqs. 7 and 8.

[r.sub.HEA] = [k.sub.HEA] HEA / [k.sub.HEA] - ST [r.sub.ST] = [k.sub.ST] - ST / [k.sub.ST] - HEA The Arrhenius kinetic parameters obtained by the simulations at the two temperatures are collected in Table 1 and the corresponding values of monomer reactivity ratio are summarized in Table 2. Finally, the transition state geometries for the reactions reported in Scheme 1 are shown in Fig. 1.

TABLE 1. Arrhenius kinetic parameters for HEA/ST theoretically determined at 50 [degrees] C and 125 [degrees] C. T = 50 degrees C [Log.sub.10] [E.sub.a] [k.sub.i-j] [Log.sub.10] (A) (A) HEA* + HEA 8.04 6.03 9,221 8.35 [right arrow] [k.sub.HEA-HEA] HEA-HEA HEA* + ST [right 8.27 5.10 65,491 8.62 arrow] [k.sub.HEA-ST] HEA-ST* ST* + ST [right 8.43 8.28 672 8.47 arrow] [k.sub.ST-ST] ST-ST* ST* + HEA [right 7.67 6.40 2,167 7.97 arrow] [k.sub.ST-HEA] ST-HEA* T = 125 degrees C [E.sub.a] [k.sub.i-j] HEA* + HEA 6.03 110,326 [right arrow] [k.sub.HEA-HEA] HEA-HEA HEA* + ST [right 5.10 664,618 arrow] [k.sub.HEA-ST] HEA-ST* ST* + ST [right 8.28 8,314 arrow] [k.sub.ST-ST] ST-ST* ST* + HEA [right 6.40 28,560 arrow] [k.sub.ST-HEA] ST-HEA* Activation energies ([E.sub.a]) in kcal [mol.sup.-1], A and [k.sub.i -j] in L [mol.sup.-1] [s.sup.-1]. [k.sub.i -j] = A * exp (-[E.sub.a]/R/T). TABLE 2. Monomer reactivity ratios for copolymerization of 2-hydroxyethyl acrylate (HEA) and styrene (ST). T = 50 [degrees] C f = 125 [degrees] C [r.sub.HEA] 0.14 0.17 [r.sub.ST] 0.31 0.29

The values of the [k.sub.ST] - ST rate coefficient in Table 1 are in reasonable agreement with those proposed in the literature for the ST homopolymerization and measured by PLP (237 and 2318 [Lmol.sub.-1] [s.sub.-1] at 50 [degrees] C and 125 [degrees] C, respectively) [60]. In particular, the Arrhenius pre-expo-nential for this elementary reaction differs slightly from that published in an our previous work [6] due to the modification introduced in its evaluation, as illustrated in the Computational Details section. No check is possible for the other kinetic constants due to the lack of experimental data. On the other hand, another interesting remark can be done comparing the rate coefficient values to the transition state geometries [7]. The calculated values exhibit the following ranking: [k.sub.HEA-ST] > [k.sub.HEA-HEA] > [k.sub.ST-HEA] > [k.sub.ST-ST] This trend in reactivity is reflected by the values of the distance (d) between the two carbon atoms involved in each radical-monomer reaction, shown in Fig. 1. In fact, d decreases as the corresponding rate coefficients ([d.sub.HEA-ST] < [d.sub.HEA-HEA] < [d.sub.ST-HEA] < [d.sub.ST-ST]): larger distances indicate higher reactivities, as the radical and the monomer react dissipating less energy to approach each other.

Moving to the analysis of the monomer reactivity ratios, they exhibit a weak dependence on temperature. Moreover, to understand the reliability of these parameters, a plot of HEA mole fraction in the copolymer (F1) as a function of HEA mole fraction in the monomer phase ([f.sub.1]) at the two investigated temperatures is reported in Figs. 2 (50 [degrees] C) and 3 (125 [degrees] C), respectively. In each figure, experimental values of copolymer composition from the literature (Ref. 8 at 50 [degrees] C and Ref. 25 at 125 [degrees] C) are compared with the prediction of the Mayo-Lewis terminal copolymerization model (Eq. 9) calculated using the monomer reactivity ratios evaluated by our computational method [47],

The results in Figs. 2 and 3 confirm the reliability of the calculated values of reactivity ratios: the terminal model reproduces the experimental data at both temperatures. Even though the prediction is less accurate at 125 [degrees] C with HEA mole fraction in the monomer phase larger than 0.5, such result is quite satisfactory. Moreover, considering that chain transfer to polymer and formation of MCRs cause the main problem in the acrylate polymerization as mentioned in the Introduction section, some deviation between experimental and calculated composition data were predictable. In fact, the calculations performed are not affected by the side reactions and yield data for the hypothetical case in which all the side reactions could be switched off. Despite of this simplification, the proposed approach is able to reproduce with reasonable agreement the real system.

[F.sub.1] = [r.sub.1] [[f.sup.2].sub.1] / [r.sub.1] [[f.sup.2].sub.1] + 2 [f.sup.1] [f.sup.2] + [r.sub.2] [[f.sup.2].sub.2]

CONCLUSIONS

A computational study on the free-radical copolymerization propagation kinetics of HEA and ST was proposed. HEA/ST copolymerization was investigated at 50 [degrees] C and 125 [degrees] C, and relative kinetic parameters were determined adopting DFT. In particular, simulations at 50 [degrees] C led to monomer reactivity ratios of [r.sub.HEA] = 0.14 and [r.sub.ST] = 0.31, whereas calculation performed at 125 [degrees] C gave values of [r.sub.HEA] = 0.17 and [r.sub.ST] = 0.29. Moreover, taking as reference experimental data proposed in previous works, the copolymer composition as a function of monomer composition was investigated, at the same temperatures indicated above. In both cases, the prediction by the terminal copolymerization model, obtained using computational parameters, well described the composition trend showed in literature, confirming the potentiality of the computational approach proposed. It can also provide useful information for the lack of experimental data relative to composition-averaged copolymerization propagation rate coefficient, even extending the set of computational data here proposed with the determination of kinetic parameters necessary to study the presence of a possible penultimate unit effect [61].

[FIGURE 1 OMITTED]

[FIGURE 2 OMITTED]

[FIGURE 3 OMITTED]

REFERENCES

(1.) H.A.S. Schoonbrood, A.M. Aerdts, A.L. German, and G.P.M. van der Velden, Macromolecules, 28, 5518 (1995).

(2.) M. Sanchez-Chaves, G. Martinez, and E.L. Madruga, J. Sci. Part A. Polyin. Chem., 37, 2941 (1997).

(3.) D. Li, N. Li, and R.A. Hutchinson, Macromolecules, 39, 4366 (2006).

(4.) W. Wang and R.A. Hutchinson, Macromol. React. Eng., 2, 199 (2008).

(5.) W. Wang and R.A. Hutchinson, Macromolecules, 41, 9011 (2008).

(6.) K. Liang, M. Dossi, D. Moscatelli, and R.A. Hutchinson, Macromolecules, 42, 7736 (2009).

(7.) M. Dossi, K. Liang, R.A. Hutchinson, and D. Moscatelli,.J. Phys. Chem. B, 114, 4213 (2010).

(8.) N.T. McManus, J.D. Kim, and A. Penlidis, Polym. Bull., 41, 661 (1998).

(9.) J.R. Lizotte and T.E. Long, Macromol. Chem. Phys., 205, 692 (2004).

(10.) J.F.G.A. Jansen, E.E.J.E. Houben, P.H.G. Tummers, D. Wienke, and J. Hoffmann, Macromolecules, 37, 2275 (2004).

(11.) G.A. Mun, Z.S. Nurkeeva, G.T. Akhmetkalieva, S.N. Shmakov, V.V. Khutoyanskiy, S.C. Lee, and K. Park, J. Polym. Sci. Part B: Polym. Phys., 44, 195 (2006)

(12.) G.A. Mun, Z.S. Nurkeeva, A.B. Beissegul, A.V. Dubolazov, P.I. Urkimbaeva, K. Park, and V.V. Khutoyanskiy, Macromol. Chem. Phys., 208, 979 (2007).

(13.) A.V. Khutoryanskaya, Z.A. Mayeva, G.A. Mun, and V.V. Khutoryanskiy, Biomacrotnolecules, 9, 3353 (2008).

(14.) M. Luck, B.R. Paulke, W. Schroder, T. Blunk, and R.H. Muller, J. Monied. Res., 39, 478 (1998).

(15.) Z. Lu, G. Liu, and S. Duncan, J. Membr. Sci., 221, 113 (2003).

(16.) A. Arun and B.S.R. Reddy, Bionzaterials, 26, 1185 (2005).

(17.) Y. Chan, T. Wong, F. Byrne, M. Kavallaris, and V. Buimus, Biomacromolecules, 9, 1826 (2008).

(18.) C.L. Russom, R.A. Drummond, and A.D. Hoffman. Bull. Environ. Contain. Toxicol., 41, 589 (1988).

(19.) K.L. Deartield, C.S. Millis, K. Harrington-Brock, C.L. Doerr, and M.M. Moore, Mutagenesis, 5, 381 (1989).

(20.) C.M. Ramos, S. Lainez, F. Sancho, M.A.G. Esparza, R. PlaNells-Cases, J.M.G. Verdugo, J.L.G. Ribelles, M.S. Sanchez, M.M. Pradas, J.A. Barcia, and J.M. Soria, Tissue Eng. A, 14, 1365 (2008).

(21.) S. Coca, C.B. Jasieczek, K.L. Beers, and K. Matyjaszewski, J. Polym. Sci. Part A: Polym. Chem., 36, 1417 (1998).

(22.) K. Bian and M.F. Cunningham, Macromolecules, 38, 695 (2005).

(23.) J.M.G. Cowie, I.J. McEwen, and D.J. Yule, Eur. Polym. J., 36, 1795 (2000).

(24.) C.D. Chow, J. Polym. Sci. Part A: Polym. Chem., 13, 309 (1975).

(25.) H. Jianying, C. Jiayan, Z. Jiaming, C. Yihong, D. Lizong, and Z. Yousi, J. Appl. Polym. Sci., 100, 3531 (2006).

(26.) M.N. Galbraith, G. Moad, D.H. Solomon, and T.H. Spurting, Macromolecules, 20, 675 (1987).

(27.) R.S. Lehrle and C.S. Pattenden, Polym. Degrad. Stab., 63, 153 (1999).

(28.) A.N. Nikitin, R.A. Hutchinson, M. Buback, and P. Hesse, Macromolecules, 40, 8631 (2007).

(29.) M. Buback, P. Hesse, T. Junkers, T. Sergeeva, and T. Theis, Macromolecules, 41, 288 (2008).

(30.) J.M. Asua, S. Beuermann, M. Buback, P. Castignolles, B. Charleux, R.G. Gilbert, R.A. Hutchinson, J.R. Leiza, A.N. Nikitin, J.P. Vairon, and A.M. van Herk, Macromol. Chem. Phys., 205, 2151 (2004).

(31.) M. Buback, Macromol. Symp., 90, 275 (2008).

(32.) J. Barth, M. Buback, P. Hesse, and T. Sergeeva, Macromolecules, 43, 4023 (2010).

(33.) W.T.K. Stevenson, R.A. Evangelista, R.L. Broughton, and M.V. Sefton, J. Appl. Polym. Sci., 34, 65 (1987).

(34.) R.H. Yocum and E.B. Nyquist, Functional Monomers: Their Preparation, Polymerization and Application, Marcel Dekker, New York, 299 (1973).

(35.) J.P.A. Heuts, R.G. Gilbert, and L. Radom, Macromolecules, 28, 8771 (1995).

(36.) J.P.A. Heuts, R.G. Gilbert, and L. Radom, J. Phys. Chem., 100, 18997 (1996).

(37.) H. Fischer and L. Radom, Angew. Chem. Int. Ed. Engl., 40, 1340 (2001).

(38.) E.I. Izgorodina and M.L. Coote, Chem. Phys., 324, 96 (2006).

(39.) V. Van Speybroeck, K. Van Cauter, B. Coussens, and M. Waroquier, ChemPhysChem, 6, 180 (2005).

(40.) D. Moscatelli, C. Cavallotti, and M. Morbidelli, Macromolecules, 39, 9641 (2006).

(41.) 1. Degirmeci, D. Avci, V. Aviyente, K. Van Cauter, V. Van Speybroeck, and M. Waroquier, Macromolecules, 40, 9599 (2007).

(42.) D. Moscatelli, M. Dossi, C. Cavallotti, and G. Storti, Macromol. Symp., 259, 337 (2007).

(43.) X. Yu, J. Pfaendtner, and L.J. Broadbelt,.1. Phys. Chem. A, 112, 6772 (2008).

(44.) 1. Degirmeci, V. Van Speybroeck, and M. Waroquier, Macromolecules, 42, 3033 (2009).

(45.) M. Dossi, G. Storti, and D. Moscatelli, Macromol. Theory Simul., 19, 170 (2010).

(46.) X. Yu, S.E. Levine, and L.J. Broadbelt, Macromolecules, 41, 8242 (2008).

(47.) F.R. Mayo and F.M. Lewis, J. Am. Chem. Soc., 66, 1594 (1944).

(48.) A.D. Becke, J. Chem. Phys., 98, 5648 (1993).

(49.) C. Lee, W. Yang, and R.G. Parr, Phys. Rev. B,37, 785 (1988).

(50.) M.W. Wong and L. Radom, J. Phys. Chem. A, 102, 2237 (1998).

(51.) V. Van Speybroeck, D. Van Neck, M. Waroquier, S. Wauters, M. Saeys, and G.B. Marin,.J. Phys. Chem. A, 104, 10939 (2000).

(52.) (a) R. Gomez-Balderas, M.L. Coote, D.J. Henry, and L. Radom, J. Phys. Chem. A, 108, 2874 (2004); (b) M.L. Coote, J. Phys. Chem. A, 108, 3865 (2004).

(53.) C. Peng, P.Y. Ayala, H.B. Schlegel, and M.J. Frisch, J. Comput. Chem., 17, 49 (1996).

(54.) M.K. Sabbe, M.F. Reyniers, V. Van Speybroeck, M. Waroquier, and G.B. Marin, ChemPhysChem, 9, 124 (2008).

(55.) S. Fascella, C. Cavallotti, R. Rota, and S. Carra, J. Phys. Chem. A, 108, 3829 (2004).

(56.) R.G. Gilbert and S.C. Smith, Theory of Unimolecular and Recombination Reactions, Blackwell Scientific Publications, Oxford, England (1990).

(57.) M.J. Frisch, G.W. Trucks, H.B. Schlegel, G.E. Scuseria, M.A. Robb, J.R. Cheeseman, J.A. Montgomery Jr., T. Vreven, K.N. Kudin, J.C. Burant, J.M. Millam, S.S. Iyengar, J. Tomasi, V. Barone, B. Mennucci, M. Cossi, G. Scalmani, N. Rega, G.A. Petersson, H. Nakatsuji, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, M. Klene, X. Li, J.E. Knox, H.P. Hratchian, J.B. Cross, V. Bakken, C. Adarno, J. Jaramillo, R. Gomperts, R.E. Stratmann, O. Yazyev, A.J. Austin, R. Cammi, C. Pomelli, J.W. Ochterski, P.Y. Ayala, K. Morokuma, G.A. Voth, P. Salvador, J.J. Dannenbcrg, V.G. Zakrzewski, S. Dapprich, A.D. Daniels, M.C. Strain, O. Farkas, D.K. Malick, A.D. Rabuck, K. Raghavachari, Foresman, J.V. Ortiz, Q. Cui, A.G. Baboul, S. Clifford, J. Cioslowski, B.B. Stefanov, G. Liu, A. Liashenko, P. Piskorz, I. Komaromi, R.L. Martin, D.J. Fox, T. Keith, M.A. Al-Laham, C.Y. Peng, A. Nanayakkara, M. Challacombe, P.M.W. Gill, B. Johnson, W. Chen, M.W. Wong, C. Gonzalez, and J.A. Pople, Gaussian 03, Revision C.02, Gaussian, Inc., Wallingford, CT (2004).

(58.) W.L. DeLano, The PyMOL Molecular Graphics System, DeLano Scientific, San Carlos, CA (2002).

(59.) T. Alfrey and C.C. Price, J. Pol)'jrt. Sci., 2, 101 (1947).

(60.) M. Buback, R.G. Gilbert, R.A. Hutchinson, B. Klumperman, F.D. Kuchta, B.G. Manders, K.F. Odriscoll, G.T. Russell, and J. Schweer, J. Macromol. Chem. Phys., 196, 3267 (1995).

(61.) T. Fukuda, Y.-D. Ma, and H. Inagaki, Macromolecules, 18, 17 (1985).

Correspondence to: Davide Moscatelli; e-mail: davidc.moscateIIi@polimi.it DOI 10.1002/pen.22045

Published online in Wiley Online Library (wilcyonlinelibrary.com). [c] 2011 Society of Plastics Engineers

Marco Dossi, Giuseppe Storti, Davide Moscatelli Dipartimento di Chimica, Materiali e Ingegneria Chimica "Giulio Natta", Politecnico di Milano, 20131 Milano, Italy

Printer friendly Cite/link Email Feedback | |

Author: | Dossi, Marco; Storti, Giuseppe; Moscatelli, Davide |
---|---|

Publication: | Polymer Engineering and Science |

Article Type: | Report |

Geographic Code: | 4EUIT |

Date: | Oct 1, 2011 |

Words: | 3885 |

Previous Article: | Polymerization of ethylene using amido functional half-sandwich complexes of group 4 metals. |

Next Article: | Erratum to "a mathematical model to study the effect of the air jet in the film blowing process". |

Topics: |