SummaryIn a copolymerization process two (or more) monomer species M1 and M2 chemically react to form longchained polymer molecules. During this process the ratio of the rates of reaction of the monomers M1 and M2 usually does not equal the ratio of the monomer concentrations in the solution. In other words, the composition of the copolymer chains produced in a very small time period usually differs from the monomer feed composition. For that reason the mole fractions of the momomers in the feed, X and (1–X), change as the copolymerization reaction proceeds. Consequently, the mole fractions of the momomers found in the polymer chains that are being produced, Y and (1–Y), change in course of time too. Considering a reactive site at the end of the chain of a growing bipolymer, it could be a radical, or ionic species, or a coordination complex. Moreover, the active site can be either a remainder of monomer M1, or a remainder of M2, depending on which one of both monomers has been the latest to couple with the growing chain. Using a rewritten version of the socalled "copolymer equation", which is applicable to many polymer systems, the program calculates the copolymer composition Y from a given X and reactivity ratios r_{1} and r_{2}. To track the composition drift, it performs reiterative calculations in which the computed instantaneous copolymer composition Y is used to compute the next feed composition X. The composition drift is visualized graphically by plotting the copolymer composition versus the feed composition. As the program runs, the compositions passed through are being represented as a lengthening curve, starting from the given X and ending as soon as the compositions no longer change. To set an example: the program lets you track the composition drift coming up for discussion in problem 7 of the preliminary examination example Polymeren I. Here the initial value of the mole fraction in the feed X equals 0,50 and the reactivity ratios r_{1} and r_{2} are 0,0 and 1,0 respectively (see also the solution help at page 3 of the document). You can download Polymeren I from here. 
Program Listing 
Elucidation and Help 
PROGRAM:COPOLYM  Title 
GridOn  > 2nd, FORMAT > GridOn, ENTER 
AxesOn:LabelOn  
FnOff  > VARS > YVARS > On/Off... > FnOff 
Plotsoff  > 2nd, STAT PLOT > PlotsOff 
0 → Xmin  "→" Symbolizes the STO> key > VARS > Window... > Xmin 
94/62 → Xmax  
.2 → Xscl : 0 → Ymin  
1→Ymax : .2→Yscl  
1 → Xres  
5/62→K  
1.5 → R  Example values r1=1.5, r2=0.1, X=0.8 
.1 → S  
.8 → F  Initial mole fraction X (monomer M1) 
ClrDraw  > 2nd, DRAW > ClrDraw 
Lbl PS  > PRGM > CTL > Lbl (Label) PS: Polymer System 
ClrHome  > PRGM > I/O > ClrHome 
Disp "Enter R1, R2, X",R,S,F  
Input R  Reactivity ratio r1 
Input S  Reactivity ratio r2 
Input F  Initial feed composition X 
max(.001,R) → R  > MATH > NUM > max( 
max(.001,S) → S  
Line (0,1,1,1)  > 2nd, DRAW > Line( 
Line (1,1,1,0)  
Line (0,0,1,1)  
F → X  
X+.0001 → U  
(RX²+X(1–X)) / (RX²+2X(1–X)+ S(1–X)²) + .0001 → V 

Lbl CD  CD: Composition Drift 
Repeat abs(V–Y)<K  > PRGM > CTL > Repeat > MATH > NUM > abs( > 2nd, TEST > "<" 
U+.7(X–U)→X  
(RX²+X(1–X))/(RX²+2X(1–X)+S(1–X)²) → Y  
End  > PRGM > CTL > End 
Line (U,V,X,Y)  
abs(Y–X) → J  
If J<.0005  > PRGM > CTL > If 
Goto CC  > PRGM > CTL > Goto 
(X–Y) / J → H  
Y → V : X → U  
X+Hmin(K,.9J) → X  > MATH > NUM > min( 
min(1,X) → X  
max(0,X) → X  
Goto CD  
Lbl CC  CC: Composition remains Constant 
Circle(X,Y,1/62)  > 2nd, DRAW > Circle( 
Pause  > PRGM > CTL > Pause 
Goto PS 
While running the program, a composition drift curve is drawn. When the program execution is interrupted, you may press the ENTER key to continue and add another curve to the graph; or press the ON key to stop the program. 
Hein ten Horn
Top  Programs Texas Instruments  Home