Propylene Polymerization using Ziegler-Natta Catalysts Mathematical Modeling, Polymerization Kinetics and Polymer Characterization Study Institute for Polymer Research Symposium May 10, 2011 University it of Waterloo Waterloo, Ontario Ahmad Alshaiban and João B. P. Soares Department of Chemical Engineering g University of Waterloo Waterloo, Ontario, Canada
Outline Background Polypropylene Ziegler-Natta Catalysts and Electron Donors Results Modeling Polymerization : Kinetic Reactor Setup Estimation of Kinetic Constants Polymerization : Characterization Molecular Weight / GPC Tacticity / 13 C NMR Crystallinity / CEF Conclusions 2
Polypropylene Background Propylene CH 3 CH CH2 Isotactic Atactic Syndiotactic CH CH 3 CH 3 CH 3 CH 3 CH 3 CH CH CH CH CH 2 CH 2 CH 2 CH 2 CH 2 CH 2 CH 3 CH 3 CH 3 CH 3 CH 3 CH CH CH CH CH CH 2 CH 2 CH 2 CH 2 CH 2 CH 2 CH 3 CH 3 CH 3 CH 3 CH 3 CH CH CH CH CH CH 2 CH 2 CH 2 CH 2 CH 2 CH 2 3
Ziegler Natta Catalyst Background Electron Donor Functionality Low isotactic Ti Mg Cl Donor :D Donor :D Vacancy :D Inactive Polymer :D Atactic Highly isotactic Ti Mg Cl Vacancy Polymer Active site models for MgCl 2 TiCl 4 (Kakugo et al., 1988) 4
Isotacticity and MWD Background Polypropylene tacticity and MWD have a significant effect on its properties. Electron donors control polypropylene tacticity. (Chadwick et al, 2001) 5
Modeling Results Mathematical Modeling 6
Modeling Modeling Results Site Transformation: Electron Donor D D Atactic isotactic Ti Mg Cl Vacancy Polymer Propagation: Ti + Ti Termination: Ti + H H Ti H + P p R p R p R tr R tf P tr R p R tr R tr R tf P tf R p R tf R tr R tf 7
Dynamic Solution Simulation Modeling Results eight (g/mol) Molecular W 25 20 15 10 5 0-5 x 10 M w M n 4 Overall : Molecular l Weights vs. Time Normal Conditions 2 [Do] 2 [H 2 ] 2 [M] 1 2 3 4 5 6 7 8 Time "Sec" " x 10 4 A. Alshaiban and J. B. P. Soares, Macromol. Symp., 285, 8 (2009) 8 A. Alshaiban and J.B.P. Soares, Macromol. React. Eng., 5, 96 (2011)
Monte Carlo Simulations Modeling Results dw / d log(mw) 1.50 I Isotactic i (799%) (79.9 Stereoblock Isotactic (91.5 %) Atactic (13.0 %) 1.50 Isotactic Atactic (2.1 %) Stereoblock Isotactic Stereoblocks (7.1 %) Stereoblocks (6.4 %) 1.00 Atactic 050 0.50 0.00 1 2 3 4 5 6 log(mw) dw / d log(mw) 1.00 Atactic 0.50 0.00 1 2 3 4 5 6 log(mw) d 0.5 x Do Do 2 x Do A. Alshaiban and J.B.P. Soares, Macromol. React. Eng., 5, 96 (2011) 23.50 23.00 22.50 22.00.50.00 20.50 20.00 19.50 19.00 9 ppm downfield TMS
Polymerization Polymerization 10
Reactor Set up Polymerization Kinetic Propylene PI FI Cooling Water rpm control TIC PC Cooling Water Elect. Heat Level-1 Level-2 DCPDMS (D-Donor) Di cyclo pentyl di methoxy silane Do/Ti (mol/mol) 1.5 0 H 2 (psi) 16 0 T ( o C) 60 70 11
Estimation of Kinetic Constants Polymerization Kinetic Propylene Uptake R p, Exp V dm dtt 2 R p k p Mathematical Formulation MC 1 e min k K k p a k d [ R p, Exp R p 0 ] 2 K a k t(1 k 1 d d K K a a ) e k d t V 12 *
Polymerization Polymerization Kinetic 60 o C Factor I II III IV Al/Ti (mol/mol) 900 (± 7.1%) 900 (± 16%) 900 (± 7.7%) 900 (± 7.2%) Do/Ti (mol/mol) 1.5 (± 7.1%) 0 1.5 (± 9.6%) 0 H 2 (psi) 0 0 16 (± 10%) 16 (± 10%) Yield (g PP g cat 1 min 1 ) 7.6 (±15%) 6.1 (±15%) 17.3 (±2.5%) 8.4 (±6.5%) l/min) Rp (mo 0.04 0.03 002 0.02 0.01 60_I_36 60_I_37 60_I_40 0.00 0 5 10 15 20 25 30 35 Time (min) 13 *
Polymerization Polymerization Kinetic mol/min) Rp (m R p Experimental vs. Calculated at 60ºC 0.08 60_I_exp 0.06 0.04 0.02 0 (D, ) 60_II_exp (, )) 60_III_exp 60_IV_exp 60_I_cal 60_II_cal 60_III_cal 60_IV_cal 0 5 10 15 20 25 30 35 Time (min) (D, H) (, H) 14
Polymerization Results mol/min) Rp (m R p Experimental vs. Calculated at 70ºC 0.08 70_I_exp 70_II_exp 0.06 0.04 0.02 0 0 5 10 15 20 25 30 Time (min) 70_III_exp 70_IV_exp 70_I_cal 70_II_cal 70_III_cal 70_IV_cal (D, ) (, )) (D, H) (, H) 15
Polymerization 2 Group II (D, ) min[( Rp Exp Rp ) ( Rp Exp Rp ) min) Rp (mol/ 0.040 0.030 0.020 0.010 0.000 k K k p a k d, 60 C, 60_II_exp 60_II_cal 70_II_exp 70_II_cal 0 5 10 15 20 25 30 35 Time (min) 60 70 E kcal/mol o C K a (min) 1 0.33 0.89 22.8 k 1 1 p (L mol 1 min ) 1.7x10 3 2.5x10 3 86 8.6 o C 70 C ] Polymerization 2 Kinetic k d (min) 1 0.010 0.041 31.8 16 *
Polymerization Polymerization Kinetic Activation Ene ergy (kcal.mo l 1) 35 30 25 20 15 10 5 0 EA Ep Ed I (D, ) II (, ) III (D, H) IV (, H) 17 *
Summary of the Estimated Kinetic Constants Polymerization Kinetic K a increases as T increases. K a increases by adding H 2 in absence of Do. K a get s the highest value at the presence of H 2 only (IV). K p increases as T increases. k p get s the highest value at the presence of Do & H 2 (III). k p increases by adding H 2 (III) & (IV). k d increases as T increases. k d increases by adding H 2 in absence of Do (IV). 18 *
Conclusion Conclusion Kinetic Developed a detailed mathematical model that describes propylene polymerization kinetics and polypropylene microstructure, taking in consideration the effect of external electron donors. ( * ) Estimated the activation energies of activation, propagation, and deactivation of a commercial heterogeneous Ziegler Natta catalyst for propylene polymerization in order to be integrated with the commercial process simulator. ( * ) A. Alshaiban and J. B. P. Soares, Macromol. Symp., 285, 8 (2009) A. Alshaiban and J.B.P. Soares, Macromol. React. Eng., 5, 96 (2011) 19
Polymerization Characterization GPC Analysis
GPC Analysis Polymerization Characterization M n decreases with H 2 at the same T. M n increases with T in the presence of Do and vice versa. PDI decreases with H 2. 2 PDI decreases with increasing T.
GPC Analysis Polymerization Characterization dw W/d(log MW) 0.8 0.7 0.6 0.5 0.4 03 0.3 0.2 2 0.35 0.3 0.25 0.2 015 0.15 Exp.# 039 0.1 5 Site Types 0.05 2 = 0.00496 0 I (donor, no H 2 ) II (no donor, no H 2 ) III (donor, H 2 ) IV (no donor, H 2 ) 2 3 4 5 6 7 Number of Site Type 0.1 0 2 3 4 5 6 7 8 9 log MW
GPC Analysis Polymerization Characterization Parameter I (D, ) II (, ) 60 o C III (D, H) IV (, H) M n (g mol 1 ) 30.1 x 10 3 22.5 x 10 3 18 x 10 3 19 x 10 3 M w (g mol 1 ) 151.6 x 10 3 140.3 x 10 3 82.9 x 10 3 77.3 x 10 3 PDI 5.0 6.2 4.6 4.1 Number of site types (n) 5 5 5 5 Parameter I (D, ) II (, ) 70 o C III (D, H) IV (, H) M 1 94x10 n (g mol ) 84.1 x 10 3 13.7 x 10 3 27.5 x 10 3 9.4 3 M w (g mol 1 ) 499.4 x 10 3 72.5 x 10 3 95.5 x 10 3 32.6 x 10 3 PDI 5.9 5.3 3.5 3.5 Number of site types (n) 5 5 4 4
Polymerization Characterization 13 C NMR Analysis
13 C NMR Analysis Polymerization 19.868 19.943 20.090 20.16 20.29 20.329 20.605 20.858.061.181.357.594.675.69.714.730.738.746.754.762.70.78.786.794.851.84.892.90.908.916.924.932 Characterization 19.867 19.946 20.090 20.169 20.298 20.342 20.691 20.857 20.898 20.936.060.182.371.593.632.674.698.729.737.745.753.761.769.777.785.793.850.883.891.899.907.915 19 19 20 20 20 20 20 20 17.5 18.0 18.5 19.0 19.5 20.0 20.5.0.5 22.0 22.5 ppm 1.052 0.873 0.904 1.258 1.841 2.252 56.43 1. 0 III (D, H) IV ( H) 17.5 18.0 18.5 19.0 19.5 20.0 20.5.0.5 22.0 22.5 ppm 849 805 051 622 948 207 101 170 000 IV (, H) 1.8 1.8 2.0 3.6 3.9 5.2 0.1 59.1 1.0
13 C NMR Analysis Seq.# Range ( d )* I (D, ) II (, ) III (D, H) Polymerization Characterization IV (, H) 1 mmmm 22.0.7 83.78 61.25 93.77 76.46 2 mmmr.7.4 3.73 9.51 3.06 6.73 3 rmmr.4.2 0.53 1.87 0.00 0.00 4 mmrr.2.0 4.17 7.76 2.48 4.96 5 mmrm + rmrr.0 20-7 2.01 6.38 0.00 4.60 6 rmrm 20.7 20.5 0.00 1.67 0.00 0.00 7 rrrr 20.5 20.25 2.90 3.55 0.00 2.61 8 rrrm 20.25 20.0 1.67 5.18 0.00 2.29 9 mrrm 20.0 19.7 1. 2.82 0.69 2.35 Range of d reported by Busico et al.(2001)
13 C NMR Analysis: Simulated H 2 effect on polypropylene tacticity Polymerization Characterization 1/2 x H2 4.28% H 2 2 x H2 12.70% 6.04% 7.19% 7.45% 3.90% 83.01% 86.77% 88.65% Isotactic chains in wt% Atactic chains in wt% Stereoblock chains in wt% A. Alshaiban, 2008, MASc. Thesis, University of Waterloo A. Alshaiban and J. B. P. Soares, Macromol. Symp., 285, 8 (2009) A. Alshaiban and J.B.P. Soares, Macromol. React. Eng., 5, 96 (2011)
Polymerization Characterization TREF CEF Analysis Injection Crystallization Elution F E T i T f T i T f CEF F c F E
CEF Analysis Polymerization Characterization dw/d dt 35 30 25 20 15 10 5 (D, ) (, ) A_70_I_158 A_70_I_159 A_70_II_157 A_70_II_147 ( D, H) (, H) A_70_III_154 A_70_III_155 A_70_IV_151 A_70_IV_153 0 18 38 58 78 98 118 Temperature ( o C)
CEF Analysis Polymerization Characterization dw W/dT 20 15 10 5 (D, ) (, ) ( D, H) (, H) A_60_I A_60_II A_60_III A_60_IV 0 18 38 58 78 98 118 Temperature ( o C)
CEF Analysis 20 116.7 C Polymerization Characterization dw/dt dw/dt 15 A60I A_60_I A70I A_70_I 10 5 0 112.1 C 18 38 58 78 98 118 30 Temperature ( o C) 25 20 15 A_60_II 10 5 0 A_70_II 109.9 C 110.8 C 18 38 58 78 98 118 Temperature ( o C) M n g/mol Group I (D, ) PDI mmmm % 30.1 x 10 3 5.0 83.78 84.1 x 10 3 5.9 95.5 M n g/mol Group II (, ) PDI mmmm % 22.5 x 10 3 6.2 61.25 13.7 x 10 3 53 5.3 63.00
dw/dt dw/dt CEF Analysis 116.3 C 15 10 A_60_III A_70_III 114.9 C 5 0 18 38 58 78 98 118 20 Temperature ( o C) 15 10 A_60_IV A_70_IV 111.5 C 111.9 C 5 M n g/mol Group III (D, H) PDI Polymerization Characterization mmmm % 18 x 10 3 4.6 93.77 27.5 x 10 3 3.5 96.90 Group IV (, H) M n PDI mmmm g/mol % 19 x 10 3 4.1 76.46 9.4 x 10 3 3.5 74.85 0 18 38 58 78 98 118 Temperature ( o C)
Conclusion Conclusion Conducted a systematic polymerization kinetics and microstructural studies for polypropylene produced using 4 th generation ZN catalyst and compared the hydrogen effect with the simulation results obtained from the developed mathematical model. Adding Do increases M n ; and at Do presence, M n increases with T as in groups I (D, ) & III (D,H). M n decreases with H 2 at the same T as we go from II (, ) toiv(, H) and from I (D, ) toiii (D,H) Number of site types decreased by one at high T in the presence of H 2. No significant change in pentad assignments with T except for group I (Do, ) Introducing H 2 tends to increase the tacticity [ I (D, ) III (D,H) ] and [ II (, ) IV (, H) ] which h is in agreement with our simulation of the developed d mathematical ti model. Group III (D, H) shows the highest crystallization peak temperature. Crystallinity increases with T in the presence of Do, I (D, ) &III (D,H). A. Alshaiban and J. B. P. Soares, Macromol. Symp., 285, 8 (2009) A. Alshaiban and J. B. P. Soares, Macromol. Symp., Accepted A. Alshaiban and J. B. P. Soares, Macromol. React. Eng., 5, 96 (2011)
The End