M(styrene)=104.15
// styrene bulk polymerization
[styrene] = 8.725 mol/L
// ATRP
[I] = 8.725e-2
[C] = 8.725e-2

kt=1.341e10 * exp(-14.34/RT)
kp=4.2658e+07* exp(-32.51/RT)
km=2.31e+6 * exp(-53.018/RT)
kth=2.2e+05*exp(-114.816/RT)

// special ATRP
ka=4.5e-1
kd=1.1e7

temperature = 110.0 °C

// styrene polymerization
styrene + styrene + styrene --kth--> P(1) + P(1)  
P(n) + styrene --kp--> P(n+1) 
P(n) + P(m) --kt--> D(n+m)
P(n) + styrene --km--> D(n) + P(1)

// ATRP
I + C --ka--> R + CX
R + CX --kd--> I + C
R + styrene --kp--> P(1)
P(n) + CX --kd--> Dr(n) + C
Dr(n) + C --ka--> P(n) + CX

Plugin TerminationHamielec { speciesStyrene = styrene, coefficient = kt }