331: Soft Materials

Kenneth R. ShullDepartment of Materials Science and EngineeringNorthwestern University

Table of Contents

Catalog Description

Different kinds of polymeric materials. Relationships between structure and physical properties; rubber elasticity, the glassy state, crystallinity in polymers. Lectures, laboratory. Prerequisites: 301 or equivalent; 314 or CHEM 342 1.

Course Outcomes

331: Soft Materials

At the conclusion of the course students will be able to:
  1. Given the chemical structure of a common polymer, draw the chemical structures of the monomer(s) from which it was made.
  2. Given the chemical structure of a monomer (or monomers), draw the chemical structure polymers that can be synthesized from it (or them).
  3. Understand, describe and calculate the structural parameters of polymeric materials including monomer units, molecular weight, tacticity, coil dimension, crystallinity, and morphology.
  4. Describe the relationship between the above structural parameters and the mechanical and thermal properties of polymeric materials.
  5. Describe how the structure and mechanical properties of polymeric materials change at the glass transition temperature and at the melting temperature.
  6. Describe how the molecular structure of a surfactant determines its micelle structure.
  7. Describe how surfactant molecules and aggregates affect the optical properties, viscosity and surface tension of polymer solutions.
  8. Calculate the intermolecular and surface forces for molecules and colloids of different geometries.
  9. Design strategies to stabilize or destabilize colloidal systems.
  10. Identify a soft material application in daily life, such as in arts, music, sports or food, and explain how material advancement has imporved the application.
test text

4 Introduction

Since the title of this book contains the term 'Soft Materials', it makes sense to define what we really mean by 'soft'. Here are two ways to think about it:
  1. Soft Materials have Low Elastic Moduli.
    By 'low' we mean significantly lower than the moduli of crystalline metals and ceramics. The jellyfish shown in Fig. 1 is obviously 'soft' in this sense. Metals and ceramics typically have moduli in the range of 100 MPa (see Fig. 2). While the strength of metals can be adjusted by a variety of mechanisms that affect the nature of dislocation motion in these systems, the modulus is set by the nature of the interatomic potentials and there nothing that can really be done to significantly affect the modulus of a given material. Polymers are different, however, and have a much broader range of elastic moduli. The stiffest of these (Kevlar™for example) have elastic moduli in at least one direction that are comparable to the modulus of steel.
  2. Thermal Fluctuations Matter in Soft Materials.
image: 23_home_ken_Mydocs_MSEcore_331_figures_jellyfish.png Figure 1: An example of a 'soft' material.
image: 24_home_ken_Mydocs_MSEcore_331_figures_Young_s_Modulus.png Figure 2: Range of Young's moduli (Pa) for different materials classes.
Exercise: How high above the earths surface must a single oxygen molecule be lifted in order for its gravitational potential energy to be increased by k B T ?
Solution: The gravitational potential energy is mgh , where m is the mass of the object, g is the gravitational acceleration (9.8 m/s 2 ) and h is height. The mass of a single O 2 molecule is obtained by dividing the molecular weight in g/mole by Avogadro's number:
m=32( g mole ) ( mole 6.02x1 0 23 ) =5.3x1 0 -23 g =5.3x1 0 -26 kg
The conversion to kg illustrates the first point that we want to make with this simple calculation: Don't mix units and convert everything to SI units (m-kg-s) to keep yourself sane and avoid unit errors. Boltmann's constant, k B needs to be in SI units as well:
k B =1.38x1 0 -23 J/K
Now we just need to equate k B T to mgh , using something reasonable for the absolute temperature, T . I'll use T =300 K to get the following: h= k B T mg = ( 1.38x1 0 -23 J/K ) ( 300 K ) ( 5.3x1 0 -26 kg ) ( 9.8 m/ s 2 ) =8000 m =8 km

Synthetic Polymers

Most of this book is devoted to synthetic polymers because of their widespread use and importance, and because they illustrate many of the key concepts that are relevant to natural polymers and other soft materials. We begin with a basic introduction to the synthesis of polymeric materials, and the z component is between properties of these materials. We conclude with some examples of soft materials made up of building blocks that are held together by relative weak forces.

5.1 What is a Polymer

A polymer is a large molecule made from many small
repeat units
or 'mers'. There is an inherent anisotropy at the molecular level because both strong and weak bonding interactions are important:
image: 25_home_ken_Mydocs_MSEcore_331_figures_Polymer.png Figure 3: Schematic representation of a polymer

5.2 Classification Scheme

Crystallization and glass formation are the two most important concepts underlying the physical properties of polymers. Polymers crystallize at temperatures below T m (melting temperature) and form glasses at temperatures below T g (glass transition temperature). All polymers will form glasses under the appropriate conditions, but not all polymers are able to crystallize. The
classification scheme
shown in Figure 4 divides polymeric materials based on the locations of T g and T m (relative to the use temperature, T ) and is a good place to start when understanding different types of polymers.
image: 26_home_ken_Mydocs_MSEcore_331_figures_classification_scheme.png
Figure 4: Classification scheme for polymeric materials.

Elastomers

Traditional elastomers are amorphous materials with a glass transition temperature less than the use temperature so that they remain flexible. The are generally crosslinked so that they do not flow over long periods of time. Common examples are shown in Figure 5.
image: 27_home_ken_Mydocs_MSEcore_331_figures_polychloroprene.png
Polychloroprene (Neoprene)
image: 28_home_ken_Mydocs_MSEcore_331_figures_polyisobutylene.png
Polyisobutylene
image: 29_home_ken_Mydocs_MSEcore_331_figures_Silicone_2.png
Silicone
Figure 5: Examples of elastomeric materials.

Glassy Polymers

Glassy polymers are amorphous like elastomers, but their glass transition temperature is above the use temperature. Because of this they behave as rigid solids, with elastic moduli in the range of 1 0 9 Pa. Glassy polymers do not need to be crosslinked, because below T g the molecules and flow of the material is suppressed. Also, because the materials are homogenous over length scales comparable to the wavelength of light, they are transparent. Common examples are illustrated in Figure 6.
image: 30_home_ken_Mydocs_MSEcore_331_figures_polymethylmethacrylate.png
Poly(methyl methacrylate) (PlexiGlas)
image: 31_home_ken_Mydocs_MSEcore_331_figures_polycarbonate.png
Polycarbonate
image: 32_home_ken_Mydocs_MSEcore_331_figures_polyphenyleneoxide.png
Poly(phenylene oxide)
Figure 6: Examples of glassy polymers.

Semicrystalline Polymers

Semicrystalline polymers must have molecular structures that are compatible with the formation of an ordered lattice. Most atactic polymers are amorphous (non-crystalline) for this reason (with the exception being examples like poly(vinyl alcohol where the side group is very small). Another requirement is that T g be less than T m . If the glass transition temperature is higher than T m the material will form a glass before crystallization can occur. In the glassy polymer the material is kinetically trapped in the glassy state, even though the crystalline state has a lower free energy below T m . Examples of some semicrystalline polymers are shown in Figure 7.
image: 33_home_ken_Mydocs_MSEcore_331_figures_polyethyleneterephthalate_bottle.png
Polyethylene terephthalate.
image: 34_home_ken_Mydocs_MSEcore_331_figures_polytetrafluoroethylene.png
Poly(tetrafluoroethylene) (Teflon)
image: 35_home_ken_Mydocs_MSEcore_331_figures_spiderweb.png
Spider Web
Figure 7: Examples of semicrystalline polymers.

5.3 Understanding Polymer Chemistry

Crystallization and glass formation processes are central to our understanding of polymeric materials, we must eventually address the following question:
Before we answer this question, however, we must answer the following question:
In order to answer this question properly, we need to study the processes by which polymeric materials are made. Polymer synthesis involves organic chemistry. After familiarizing ourselves with some of the relevant polymer chemistry, we will be in a position to study the physical properties of polymers. For this reason, our discussion of molecular structure in polymers will include some chemistry.

5.3.1 Covalent Bonding

Polymer molecules consist of atoms (primarily carbon, nitrogen, oxygen and hydrogen) which are covalently bound to one another. The fraction of the periodic table that can form strong covalent bonds is relatively small, corresponding to the ten atoms shown in yellow in Figure 8 (H, C, N, O, F, Si, P, S, Cl, Br). It is useful at this point to recall some of the basic principles governing the bonding between these atoms:
image: 36_home_ken_Mydocs_MSEcore_331_figures_Simple_Periodic_Table_Chart-en.png
Figure 8: Periodic table, illustrating the ten atoms that make up the covalently-bonded portion of polymeric materials.

5.3.2 Lewis Diagrams

Lewis diagrams (Wikipedia link) provide a convenient way of keeping track of the valence electrons in covalently bonded compounds. Several examples are given here. Note how the rules given on the previous page are followed in each case.
image: 37_home_ken_Mydocs_MSEcore_331_figures_Lewis_Diagrams.png Figure 9: Examples of Lewis Diagrams

5.3.3 Bonding

The following principles of covalent bonding in organic materials are very helpful:
(The situation for P and S is actually a bit more complicated when either of these atoms are bonded to oxygen, but these general rules serve our purpose for now.) The chemical structures throughout this book can be seen to obey these rules.

5.3.4 Shorthand Chemical Notation

Most of the chemical structures illustrated in this text are relatively simple, consisting of single and double bonds between atoms. We generally don't bother to write all of the carbons and hydrogens into to the structure. We use the following common conventions.
An example of this convention is shown below in the structure for
benzene
: :
image: 38_home_ken_Mydocs_MSEcore_331_figures_benzene.png Figure 10: Examples of Shorthand Notation when drawing chemical structures
Note the resonance between the two possible ways of drawing the double bonds in the drawing on the left. Molecules with these types of alternating
aromatic compounds
.

Polymerization Reactions

Polymerization is the process by which small molecules react with one another to form large polymer molecules. Polymerization reactions can be broken up into the following two general categories:
Step Growth Polymerizations:
Collections of A and B species react with one another. In linear step growth polymerizations, the ends of molecules react with one another to form longer molecules. A variety of reactions are possible, so you need to know at least a little organic chemistry. We'll focus on just a few of the most common cases.
Chain Growth Polymerizations:
Each polymer chain has one reactive site to which additional monomers are added.
Additional Resource:
The Macrogalleria web site has some excellent, simple descriptions of polymerization reactions. Specific examples are referenced at different points throughout this book.
Click here for Macrogalleria polymerization overview.

6.1 Step-Growth Polymerizations

image: 39_home_ken_Mydocs_MSEcore_331_figures_Step_Growth_Animation_2.png Figure 11: Polymers produced by step growth polymerization. The red and green circles correspond to 'A' and 'B' monomers that react with one another to form the polymer.
In this example green and red monomers can only react with each other. Because there are 5 more red monomers than green monomers, there are 5 molecules remaining at the end of the reaction, with each of these molecules possessing two red end groups. The extent of reaction, p is defined as the fraction of available reactive groups which have actually undergone a reaction. Values close to one are needed in order to obtain useful, high molecular weight polymer. A delicate stoichiometric balance generally needs to be maintained (same amount of red and green monomers) in order to obtain high molecular weight.
The following pages illustrate some of the specific reactions which take place during the polymerization process. To illustrate the concepts involved, we consider the following polymer types, all of which are produced by step-growth polymerization:
  1. Section 6.1.1
  2. Section 6.2.1
  3. Section 6.2.2
  4. Section 6.2.3

6.1.1 Polyamides

image: 40_home_ken_Mydocs_MSEcore_331_figures_amide_formation.png Figure 12: Formation of an amide from an amine and a carboxylic acid
In this example, a primary amine reacts with a carboxylic acid to form an amide linkage. Water is liberated during the condensation reaction to form the amide. Primary amines and acid chlorides undergo a similar reaction:
image: 41_home_ken_Mydocs_MSEcore_331_figures_amide_from_acid_chloride.png Figure 13: Formation of an amide from a primary amine and an acid chloride
Acid chlorides react very rapidly with amines at room temperature, which is very useful for demonstration purposes. Acid chlorides can also react with water to form carboxylic acids, however, and commercial polyamides are generally produced by reaction with carboxylic acids.

6.2 Interfacial Polymerizations

image: 42_home_ken_Mydocs_MSEcore_331_figures_Nylon_Rope_Trick.png Figure 14: Schematic representation of the interfacial polymerization off nylon.

6.2.1 Polyesters

Polyesters can be formed by condensation reactions of alcohols with carboxylic acids:
image: 43_home_ken_Mydocs_MSEcore_331_figures_ester_formation.png Figure 15: Formation of an ester from an alcohol and a carboxylic acid.
Polyethylene terephthalate (Mylar, Dacron, 2L soda bottles) is a common example. Note that under appropriate conditions, the reverse reaction (hydrolysis) reaction can also take place, where the addition of water to an ester bond forms the acid and the alcohol. This reaction is important in a variety of polymers used in biomedical applications, which degrade in the body via hydrolysis of the polymer. Polycaprolactone is one example.
image: 44_home_ken_Mydocs_MSEcore_331_figures_ester_hydrolysis.png Figure 16: Ester hydrolysis: The reverse of the esterification reaction.

6.2.2 Polyurethanes

Polyurethanes are formed by the reaction of isocyanates with alcohols, as shown below:
image: 45_home_ken_Mydocs_MSEcore_331_figures_urethane_formation.png Figure 17: Urethane Formation
Note that this is NOT a condensation reaction, since no byproducts are formed during the reaction. Also note that the urethane linkage contains an oxygen atom in the backbone, whereas the amide linkage does not. (The R 1 and R 2 substituents can have different structures, but are always attached to the illustrated linkages by carbon atoms.)

6.2.3 Epoxies

All epoxies involve reactions of epoxide groups (3-membered rings containing an oxygen atom) with curing agents. Amine curing agents are very common, as illustrated Figure , which shows a primary amine reacint with an epoxide group).
image: 46_home_ken_Mydocs_MSEcore_331_figures_epoxy1.png Figure 18: Reaction of an epoxide with a primary amine.
The secondary amines which remain can react with additional epoxide groups to form a branched structure as shown in Figure 19. The reactive functionality of a primary diamine is 4 when the reaction is with an epoxide (as opposed to its functionality of 2 in the case where the diamine reacts with an acid or acid chloride).
image: 47_home_ken_Mydocs_MSEcore_331_figures_epoxy2.png Figure 19: Reaction of an epoxide with a secondary amine.

Molecular Weight Distributions

The repeat unit for this polymer consists of a red and green unit together. The degree of polymerization ( N ) is the number of repeating units for a single molecule. Because of the statistical nature of the reaction process, not all molecules at any given point in the polymerization will have the same size. This point is illustrated in Figure 20, which shows schematic representations of polymer with and N=10.5 . The molecular weight, M , of a polymer is obtained from N by multiplying by M 0 , the molecular weight of a repeating unit of a polymer chain:
M=N M 0 (7.1)
image: 48_home_ken_Mydocs_MSEcore_331_figures_Molecular_Weight_1.png
image: 49_home_ken_Mydocs_MSEcore_331_figures_Molecular_Weight_2.png
Figure 20: Linear, step-growth polymers with N=3 (left) and N=10.5 (right).
Some key definitions:
N : Degree of polymerization - the number of repeat units in a single molecule
M : Molecular weight - molecular weight of a polymer molecule (g/mole)
M 0 : Molecular weight per repeat unit (g/mole)
The complete distribution of molecular sizes can be described by a histogram. There are two different quantities that we can keep track of: n N and w N , which are defined as follows:
The total number of molecules ( n ) or total weight ( w ) is obtained by summing over all possible values of N :
n= n 1 + n 2 + n 3 +... n = N n N (7.2)
w= N w N = M 0 N N n N (7.3)
Note: We've assumed here for simplicity that N takes on integer values. This is true for chain growth polymerizations, but for step growth polymers N can formally have non-integer values (like the value of 10.5 in Figure 20). In this case you can just imagine that distribution has divided into 'bins' of width '1'.
image: 50_home_ken_Mydocs_MSEcore_331_figures_Polymerization_Histogram.png Figure 21: Histogram of molecular sizes. The same distribution is shown, with the numbers of molecules of each size ( n N 1*2/1tted on the left, and the total weight corresponding to each of the different molecular sizes ( w N 53plotted on the right.

7.1 Number Average Molecular Weight

The
number average molecular weight
, M n , is simply the total weight of all the polymer, divided by the number of polymer molecules (typically in moles):
M n =w/n (7.4)
In more formal terms, we can replace w and by the respective sums from Eqs. 7.2 and 7.3:
M n = M 0 N N n N N n N (7.5)
By dividing by the repeat unit molecular weight, M 0 , we obtain the
number average degree of polymerization
, N n :
N n = N N n N N n N (7.6)

7.2 Weight Average Molecular Weight

The
weight average molecular weight
, M w , is similar to number average, but we use the distribution of weights instead of the distribution of molecules:

M w = M 0 N w N N w = M 0 N w N N N w N = M 0 N n N N 2 N n N N (7.7)
The
weight average degree of polymerization
, N w , is obtained simply by dividing M w by M 0 , i.e. N w = M w / M 0 . The
polydispersity index
, given by M w / M N must be greater than or equal to one. Polymers with a polydispersity index close to 1 are said to be monodisperse, with a narrow distribution of molecular weights. While Eq. 7.7 is useful in a formal sense, it's a bit cumbersome to use. Suppose we have a series of molecular weight fractions, with w i between the weight of each fraction and M i being the weight average molecular weight of each fraction. The weight average molecular weight for this collection of molecules is given by the following summation:
M w = i w i M i (7.8)
Exercise:
Suppose 3g of a monodisperse polymer with M=50,000 g/mol is blended with 5g of a monodisperse polymer with M=85,000 g/mol. What are the values of M n , and the polydispersity index for the resultant mixture?
Solution:
To get M n , we need to divide the total weight of polymer (3g + 5g = 8g) by n, the total number of moles of polymer molecules:
n= 3g 50,000 g/mol + 5g 85,000 g/mol =1.9x1 0 -4 mol
M n = 8g 1.19x1 0 -4 mol =67,200 g/mol
The average molecular weight is obtained from the definition of this quantity:
M w = M 0 N=1 w N N w = 1 w N=1 w N M N (7.9)
In our example, we obtain the following:
M w = 1 8g ( 3g*50,000 g/mol +5g*85,000g/mol ) =71,875 g/mol
Polydispersity Index = M w M n = 71,875 g/mol 67,200 g/mol =1.07
Exercise:
What if the polymers in the previous example are not monodisperse? Suppose the molecular weights quoted are weight averages, the low molecular weight polymer has a polydispersity of 1.6, and the high molecular weight polymer has a polydispersity of 2.5.
Solution:
We divide by the polydispersities to get number average molecular weights, and use these values to get the total number of molecules:
M n ( poly.a ) = 50kg/mol 1.6 =31kg/mol
M n ( poly.b ) = 85kg/mol 2.5 =34kg/mol
n= 3g 31,000g/mol + 5g 34,000g/mol =2.44x1 0 -4 mol
M n ( blend ) = 8g 2.44x1 0 -4 mol =32,800g/mol
Due to the form of the expression for the average molecular weight, we can modify the standard expression for the Mw by summing over all different components of the blend, and using the weight average molecular weight for each component:
M w ( blend ) = 1 w i w i M wi (7.10)
In our case M w is the same as it was in the previous example:
M w = 1 8g ( 3g*50,000g/mol+5g*85,000g/mol ) =71,875g/mol
Polydispersity index = M w M n = 71,875g/mol 32,800g/mol =2.19

7.3 Molecular Weight Evolution during Step Growth Polymerization

Suppose n a moles of difunctional A molecules with molecular weight M a react with n b difunctional B molecules with molecular weight M b in a step growth polymerization. If n a > n b , then all of the B molecules can react, and n a - n b molecules will remain at the end of the polymerization:
n min = n a - n b
For condensation reactions, a mole of condensation product with molecular weight M c material as a wholeecule of B that reacts. The total weight, w , of polymer produced is:
w= n a M a + n b M b -2 n b M c (7.11)
The maximum number average molecular weight, M n max , is obtained by using this value for w , an setting n equal to n min :
M n max = w n min = n a M a + n b M b -2 n b M c n a - n b (7.12)
For high molecular weight polymers, we can make some additional approximations, since know that n a - n b << n b in this case. We can therefore simplify M n to the following:
M n max n b ( M a + M b -2 M c ) n a - n b (7.13)
The repeat unit molecular weight, M 0 is obtained by adding the molecular weights for the A and B monomers, and subtracting out the molecular weight of the 2 condensation products that are formed (since a repeat unit contains 2 A-B linkages, and therefore involves 2 condensation reactions):
M 0 = M a + M b -2 M c (7.14)
From this we obtain the following expression for M n max :
M n max n a M 0 n a - n b (7.15)
The next thing we would like to sort out is how n , the total number of molecules, varies with the extent of reaction, p . The trick is to realize that every time a reaction takes place, the number of molecules in the system decreases by 1. Each molecule has two groups (one at each end) that can react, and two groups take part in each reaction. For n we obtain the following:
n= n a + n b -p( n a + n b ) =( n a + n b ) ( 1-p ) (7.16)
In this equation n a + n b is the initial number of molecules that we start off with (when p =0), and p( n a + n b ) gives the decrease in the number of molecules resulting from all of the reactions that have taken place for some finite value of p . The total mass of the polymer that has been produced is ( n a + n b ) M 0 /2 (ignoring the fact that the ends of all the molecules some extra mass that may potentially leave as part of a condensation product). The weight average molecular weight is simply w/n , from which we obtain the following:
M n = w n = ( n a + n b ) M 0 2( n a + n b ) ( 1-p ) = M 0 2( 1-p ) (7.17)
As mentioned above, the degree of polymerization that can be obtained from a homogeneous mixture of A and B monomers in a linear stop growth polymerization is limited by stoichiometry. Reaction of all A monomers requires the existence of unreacted B monomers, and vice versa. If n a n b , then some of the functional groups on either the A or B monomers must remain unreacted, so that p has a maximum value somewhat less than 1. Let's consider the case where n a - n b , in which case we know there are n a - n b molecules left when all the B groups have reacted. Each of these molecules has an unreacted A group at each end, so there are a total of 2( n a - n b ) unreacted groups at the end of the reaction. At the beginning of the reaction (when p =0) we had n a + n b reactive groups, so the total number that have reacted is 2( n a + n b ) -2( n a - n b ) =4 n b . This corresponds to a fraction of reacted groups which we refer to as p max :
p max = 4 n b 2( n a + n b ) = 2 n b n a + n b (7.18)
Another way to obtain M n max is to use Eq. 7.17 for M n , with p= p max :

M n max = M 0 2( 1- p max ) = M 0 ( n a + n b 2( ( n a - n b ) (7.19)
which with n a n b is equivalent to Eq. 7.15.
From the information given so far we know how the number average degree of polymerization evolves with p , but in many cases we would like to know what the full molecular weight distribution might look like. To do this we need to account for the statistics, and use the more detailed definitions for M n (Eq. ) and M w (Eq. ) A molecule formed from a linear step growth polymerization that has a degree of polymerization of N was made from 2N molecules ( N of type A and N of type B). Because these 2N original molecules are all bifunctional, there are a total of 4N reactive groups in the monomers that reacted together to form a polymer of length N . Only two of these groups remain unreacted (the end groups). Because each reaction involves two groups, this means that 2N-1 reactions have taken place ( 4N-2 groups have reacted).
The key assumption is that all groups react with equal probability. This probability must therefore be equal to the fraction of functional groups which has reacted, which has already been defined as the extent of reaction,p. The probability that any single one of these 2N-1 reactions has taken place is therefore equal to p . We have to multiply all of these probabilities together to get the probability of having a polymer with a degree of polymerization of N , so we end up with the following:
n N p 2N-1 (7.20)
The weight, w N , corresponding to each fraction is proportional to n N N :
w N α N p 2N-1 (7.21)
The shapes of these distributions for n N and w N are shown in the plots below, and can be replotted with different values of p in the MATLAB modulue referred to below.
image: 51_home_ken_Mydocs_MSEcore_331_figures_Molecular_Sizes_Ideal_Distribution_Interactive_2.png Figure 22: Distribution of molecular sizes during an ideal step growth polymerization.
  1. What happens to the polydispersity index as the degree of polymerization increases?
  2. Where are the number and weight averaged degrees of polymerization in comparison to the peak in w N ?
  3. How does the shape of the molecular weight distribution change as N gets very large?
  4. What sizes of molecules are present in the highest numbers?
  5. For what values of p is Nn equal to 0.5, 1, and infinity? Can you develop a simple explanation for these results? What are the values of the polydispersity for these values of p?
Analytic expressions can be derived for number-average and weight average degrees of polymerization for ideal step growth polymerizations that follow the distribution described by Eq. 7.20. Substitution of this expression into Eq. 7.5 gives the following for M n :
M n = M 0 NN p ( 2N-1 ) N p ( 2N-1 ) (7.22)
We assume that our material is isotropic N ,i.e., N =1/2, 1, 3/2, 2, etc. Half integer values of N are allowed because N is the number of repeat units in the molecule. (A value of 3/2 for N would correspond to a molecule with one A monomer between two monomers, for example). Alternatively, we can replace N by L/2 , where L takes on all positive integral values:
M n = M 0 2 L p ( L-1 ) p ( L-1 ) (7.23)
We can follow a similar procedure to get an expression for M w . Combination of the ideal, linear step growth size distribution (Eq. 7.20) with the definition of M w (Eq. 7.7), gives
M w = M 0 N N 2 p ( 2N-1 ) NN p ( 2N-1 ) (7.24)
we can again let L=2N, so that L takes on all integer values, so that the expression for M w becomes:
M w = 1 2 M 0 N L 2 p ( L-1 ) NL p ( L-1 ) (7.25)
These expressions for M n and M w can be simplified by using the following mathematical formulas:
L=1 p L-1 = 1 1-p L=1 L p L-1 = 1 ( 1-p ) 2 L=1 L 2 p L-1 = ( 1+p ) ( 1-p ) 3 (7.26)
Substituting the results for these sums into Eqs.7.23 and7.25 results in the following:
M n = M 0 2( 1-p ) ; N n = 1 2( 1-p ) ; (7.27)
M w = M 0 2 { 1+p 1-p } (7.28)
The expression for M n is the same as the result that we got from the previous section, just by keeping track of the total number of molecules in the system. Calculation of M w is a bit more complicated though, and requires the statistical approach used in this section. For the polydispersity index we get:
M w M n =1+p (7.29)
Useful, high molecular weight polymers produced by step-growth polymerization are therefore have a polydispersity index that is close to 2.

Non-Linear Step Growth Polymerizations

image: 52_home_ken_Mydocs_MSEcore_331_figures_Nonlinear_Step_Growth_Polymerization.png Figure 23: Nonlinear Step Growth Polymerization
Nonlinear step-growth polymerizations involve the use of monomers that have a reactive functionality, f , which is greater than 2.

8.1 Gelation

A network (or gel) is formed when the molecular weight of the largest molecule in the system reaches infinity. Everything which is covalently attached to this single very large molecule is part of the gel fraction.
image: 53_home_ken_Mydocs_MSEcore_331_figures_Gelation.png Figure 24: Molecular schematic of gel network.
Unattached molecules which are not part of the gel fraction make up thesol fraction . Thegel point refers to the point in the polymerization reaction at which point a gel fraction first appears, and the overall process is referred to asgelation .

8.2 Prepolymers

In many cases, the length of the polymer between branch points is determined by the molecular weight of a prepolymer . Prepolymers generally have a reactive functionality of 2.
image: 54_home_ken_Mydocs_MSEcore_331_figures_Prepolymer.png Figure 25: Example of a Prepolymer
The following factors determine the property of the network:
  1. Molecular weight of prepolymer
  2. Composition of prepolymer
  3. Functionality of curing agent (4 for diamine curing of epoxies)
  4. Extent of reaction (number of dangling ends)
  5. Network structure (loops, entanglements, etc.)

8.3 Corothers Theory of Gelation

We begin with n a molecules with a reactive functionality of f a that are able to react with n b molecules with a reactive functionality of f b . For illustrative purposes, consider the perfectly stoichiometric case where the total functionality of all the A monomers is equal to the total functionality of all the B monomers:
n a f a = n b f b (8.1)
We define an average reactive functionality, f av , so that ( n a + n b ) f av is equal to the total number of functional groups:
( n a + n b ) f av = n a f a + n b f b (8.2)
The number of initial molecules corresponding to a given number of functional groups is inversely proportional to f av . For large values of f av , p does not have to reach unity in order for the total number of molecules to be reduced to one. An expression for N n can be obtained substituting p f av /2 for p in Eq. , describing the evolution of M n with p for linear step growth polymerization:
N n = 1 2-p f av (8.3)
The quantity p f av describes the average number of times that any given molecule has reacted. The number average molecular weight is predicted to diverge to infinity when this average number of reactions per molecule is 2. This occurs when p= p gel , where p gel Standardextent of reaction at the endpoint. Setting p gel f av =2 gives:
p gel = 2 f av (8.4)
This simple equation for the gel point is very useful, although it is important to remember that this derivation suffers from the following deficiencies:
  1. Gelation actually occurs for finite N n , since N= only for the largest molecule.
  2. The theory neglects loop formation (reactions between two portions of the same molecule). These reactions change the shape of a given molecule, but do not decrease the overall number of molecules.
The first assumption tends to overestimate p gel , and the second assumption tends to underestimate p gel . A more accurate description of gelation requires a much more detailed theory than the one presented here.
When one molecule which spans the entire sample, it can no longer flow like a liquid. One the characteristics of non-linear step growth polymerizations is that the viscosity (resistance to flow) of the reaction mixture increases as the extent of reaction increases, and eventually diverges at the gel point. Also, the reaction rate generally decreases as the reaction proceeds and the reactive molecules become larger and larger. These features are illustrated conceptually in Figure 26. Here we show schematic representations of the time-dependent extent of reaction, along with the time dependence of the number average degree of polymerization, N n , and of two physical properties of the material: η 0 and E . Here η 0 is the limiting viscosity measured at very low shear rates, and E is the elastic modulus measured at very long times. The viscosity is a characteristic of a liquid material, and characterizes the material for values of the extent of reaction, p , that are less than the extent of reaction corresponding to the gel point, p gel . For values of p that exceed p gel , the material behaves as a solid, and has a finite elastic modulus, E . At the gel point, the viscosity diverges to infinity, since it is no longer possible for the material to flow. Instead, the material becomes an elastic solid, with a value of E that increases from 0 at the gel point.
Note that that materials near the gel point are viscoelastic, and can no longer be described by a single value of the viscosity or modulus. Instead, these properties depend on the timescale of the measurement, as described in more detail on the section on viscoelasticity.
image: 55_home_ken_Mydocs_MSEcore_331_figures_Gelation_properties.png Figure 26: Time-dependence of property development in a nonlinear step growth polymerization.
Exercise:
A very simple epoxy formulation has bisphenol A diglycidyl ether as the epoxide component and hexamethylene diamine as the crosslinker:
H2N( CH2 ) 6NH2
A) Hexamethylenediamine
image: 56_home_ken_Mydocs_MSEcore_331_figures_Bisphenol_A_diglycidyl_ether.png
B) Bisphenol A diglycidyl ether
  1. How many grams of hexamethylene diamine should be added to 1g of Bisphenol A diglycidyl ether in order to optimize the polymerization reaction in a homogeneous solution?
  2. Estimate the fraction of amine groups that need to react in order to reach the gel point.
Solution:
n a = n b f b n a = n b 2
We have 1g of B, which we divide by it's molecular weight to get the number of moles:
So we know we need about 1.45x1 0 -3 moles of the hexamethylenediamine. Multiplying by its molecular weight of 116 g/mole gives a total mass of 0.17g. Note that we need a lot less of A than B in this case. If we want our two component epoxy to consist of two parts that we combine in nearly equal volumes, we'll have to add some extra stuff to part A (the diamine part) that doesn't take part in the curing reaction.
To figure out the average functionality, we arrange Eq. 8.2 to give the following for f av :
f av = n a f a + n b f b n a + n b
Combining this with Eq. 8.4 for the gel point gives:
p gel = 2( n a + n b ) n a f a + n b f b
Substituting our stochiometric condition that n a = n b /2 gives p gel =0.75 .

9 Chain-Growth Polymerizations

In chain growth polymerization monomers are added one at a time to a reactive site at that typically remains at the end of the monomer that has been most recently added. There are three phases of the reaction that need to be considered in general:
  1. Initiation: In this step the a reactive species is formed that is able to form covalent bond with a monomer, reforming the active species at the end of the growing polymer chain.
  2. Propagation: Monomers are added one at a time to the growing chain.
  3. Termination: The reactive undergoes a reaction of some sort that causes it to lose its reactivity toward other monomers. The polymer chain is 'dead' at this point, and no longer grows in length.
image: 57_home_ken_Mydocs_MSEcore_331_figures_Chain_Growth_Polymerization_Animation_2.png
Figure 27: Chain growth polymerization.

9.1 Chain Growth Polymerization Mechanisms

Chain growth polymerizations generally occur either by addition to a double bond, or by opening a ring. In both cases the overall number of bonds is conserved, and there are no condensation products. If the cases where chain growth occurs by addition to double bond between two carbon atoms, the following takes place:
image: 58_home_ken_Mydocs_MSEcore_331_figures_chain_growth_types.png
Figure 28: The two basic types of chain growth polymerizations.
In a ring opening polymerization, the following takes place:

9.2 Reactive Species for Chain Growth Polymerization

Three of the most common active sites are a free radical, negatively charged anion, and a positively charged cation, illustrated in Figure 29 for the addition to the double bond of vinyl polymers. Vinyl polymers have the general structure C H 2 CHR , where R is something other than hydrogen. In each of these cases, a single monomer repeat unit is added to the end of the chain, and the reactive site moves to the end of the chain, on the repeat unit that has just been added. Reactions in organic chemistry are all about keeping track of what the bonding electrons are doing, and Lewis diagrams are very helpful in this sense. In the figures designed to illustrate different propagation and termination reactions, we just show the Lewis structures for some of the bonds, to make it simpler to keep track of the situation before and after the reaction has taken place. For our purpose we are not as interested in the detailed reaction mechanism, which would require that we provide a bit more information about the structure of some of the short-lived reaction intermediate.
Free Radical Polymerization
image: 59_home_ken_Mydocs_MSEcore_331_figures_radical_addition.png
Anionic Polymerization
image: 60_home_ken_Mydocs_MSEcore_331_figures_anionic_addition.png
Cationic Polymerization
image: 61_home_ken_Mydocs_MSEcore_331_figures_cationic_addition.png
Figure 29: Propagation mechanisms for three common chain growth polymerization mechanisms.

9.3 Initiation

The previous section outlines the chemistry of some different propagation reactions, where monomers are added to an active site. Initiation is the step by which the active site is produced at the beginning. Here we give two simple examples, one for initiation of a radical polymerization, and the second being the initiation of a an anionic polymerization.

9.3.1 Initiation of a free radical polymerization

In this example, illustrated in Figure 30, a peroxide bond between two oxygen atoms splits to form two free radicals. These radicals are able to add to the double bond in a vinyl monomer like styrene to initiate the polymerization reaction.
image: 62_home_ken_Mydocs_MSEcore_331_figures_dicumyl_peroxide.png
image: 63_home_ken_Mydocs_MSEcore_331_figures_dicumyl_peroxide_initiation.png
Figure 30: Initiation off a radical polymerization reaction with an epoxide that decomposes into two radical species at elevated temperature.

9.3.2 Initiation of an anionic polymerization

In the example shown in Figure 31 is very simple conceptually. The reaction is done in solution, and a small amount of secondary butyl lithium is simply added to a solution of styrene molecule in an appropriate organic solvent. The carbon-lithium bond is very reactive, and has a lot of ionic character. It can really be viewed as an existing anion that is ready to react directly with the styrene monomer.
image: 64_home_ken_Mydocs_MSEcore_331_figures_anionic_initiation.png
Figure 31: Initiation off a radical polymerization reaction with an epoxide that decomposes into two radical species at elevated temperature.

9.4 Termination Mechanisms

Termination is the process by which reactive chain ends become unreactive. Termination reactions can be avoided in 'living polymerizations' like anionic and cationic polymerizations.

9.4.1 Radical Recombination

Termination reactions cannot be eliminated in free radical polymerizations, since two free radicals can always combine directly to form a single bond, as illustrated below in Figure 32.
image: 65_home_ken_Mydocs_MSEcore_331_figures_radical_recombination.png
Figure 32: Termination by radical recombination.
As unpaired electrons, free radicals are very reactive towards one another. As illustrated in this example, two free radicals can readily combine with one another to form a covalent bond. Combination is therefore one type of termination reaction which is very prominent with free radical polymerizations. Note that the number of molecules decreases by one during the combination reaction.

9.4.2 Disproportionation

image: 66_home_ken_Mydocs_MSEcore_331_figures_disproportionation.png
Figure 33: Termination by disproportionation. The circled proton, along with an electron, moves from the molecule on the right, terminating the free radical on the right. The remaining electrons on the molecule on the left rearrange themselves into a double bond.
Plain Laytionation reactions can be viewed as the transfer of a proton andEntanglementsfrom one active molecule to another. The animation here illustrates how the transfer results in the termination of both molecules, with the formation of a double bond. In the final state, the carbon atoms all have eight electrons (including shared electrons) in the valence shell. The molecular weights of the two polymer molecules remain essentially unchanged (with the exception of the transfer of a single proton).

9.4.3 Chain Transfer

image: 67_home_ken_Mydocs_MSEcore_331_figures_chain_transfer.png Figure 34: Chain transfer to toluene.
Chain transfer refers to the migration of the active free radical from one molecule to the other. In this example, the active radical from the growing polymer chain is terminated by the addition of a proton and an electron from the toluene molecule on the right. The net result is that the free radical is transferred from the polymer molecule to the solvent molecule (toluene), which can then initiate the polymerization of additional monomer. Note that the process is very similar to disproportionation, except that the species which donates the proton and electron does not already have a free radical.

9.5 Intramolecular Chain Transfer

image: 68_home_ken_Mydocs_MSEcore_331_figures_chain_transfer_-_intramolecular.png
Figure 35: An Animation of Intramolecular Chain Transfer
In this example, the active radical moves from the end of the polymer chain to a different portion on the same polymer chain. Polymerization continues from this radical, resulting in the formation of a short branch, in this case consisting of four carbons. Because the branches are randomly placed along the polymer backbone, they interfere with the polymer's ability to crystallize.
Exercise:
  1. Chain transfer agents, like toluene in the previous example, can be intentionally added to give some control over the molecular weights of polymers synthesized by free radical polymerizations. In a qualitative sense, how will the addition of chain transfer agents change the number average molecular weight of the polymer?
  2. How does your answer to this question change if the chain transfer is always to other sites within the same polymer molecule?
  3. Termination reactions involving free radical polymerizations involve the elimination of two radicals to form a new bond. Where does this new bond appear for combination and disproportionation reactions?

9.6 Polymerization of Dienes

Dienes are monomers with two double bonds, with a single bond between them. Examples include butadiene, isoprene and chloroprene. The situation for isoprene is illustrated in Figure 36. The monomer has 4 carbon atoms in a line, which we number 1 through 4. Double bonds connect carbons 1 to 2 and 3 to 4, and a single bond connects carbon 2 to 3. Carbons 2 and 3 are distinguishable from one another because carbon 2 is bonded to a methyl group that is not present on carbon 2. When a propagating polymer chain interacts with the isoprene monomer, addition can occur in any of 4 different ways, which are illustrated in Figure 36:
Natural rubber is a naturally occurring version of polyisoprene that is harvested from certain tropical trees. The excellent elastomeric properties of this material arise from the pure cis 1-4 microstructure that is produced during the natural polymerization of this polymer. Because of the zigzag nature of the polymer backbone, it does not readily crystallize. A segment of the backbone of the polymer is shown in Figure 37. This schematic shows a portion of the molecule in a fully extended form. Because of the easy rotation about single bonds, the molecules actually exist in a collection of random configurations.
image: 69_home_ken_Mydocs_MSEcore_331_figures_isoprene_polymerizations.png
Figure 36: Polymerization of Isoprene
image: 70_home_ken_Mydocs_MSEcore_331_figures_polyisoprene_cis_1-4.png
Figure 37: The chemical structure of natural rubber (cis 1-4 polyisoprene). .

9.7 Living Polymerizations

Living polymerizations are chain growth polymerizations that proceed without termination or chain transfer reactions. Relatively monodisperse polymers ( M w / M N < 1.1) can typically be obtained when the initiation rate is faster than than the propagation rate. Block copolymers are formed by the sequential addition of two (or more) different monomer types in a living polymerization.

9.8 Vinyl Polymers and Tacticity

Vinyl polymers
are an important class of polymers produced by chain growth polymerizations. They are produced from vinyl monomers with the following general structure:
image: 71_home_ken_Mydocs_MSEcore_331_figures_vinyl_monomer.png Figure 38: Chemical structure of a vinyl monomer.
When a vinyl monomer is polymerized and we stretch the resulting polymer chain out so that the backbone C-C bonds all lie on the same plane, the 'R' groups end up on different sides of the molecule. The distribution of these R groups determines the
tacticity
, according to the following definitions:
These different tacticities are illustrated in Figure 39.
image: 72_home_ken_Mydocs_MSEcore_331_figures_tacticity.png Figure 39: Illustration of different polymer tacticities.
Tacticity is significant because it determines the ability of a polymer to crystallize. The disordered structure of an atactic polymer is inconsistent with the ordered structure of a crystalline polymer. As a result, atactic polymers generally cannot crystallize.Standardons to this rule include polymers where the 'R' group is very small, so that this group can be incorporated into an ordered crystalline array, even if it is randomly placed along the polymer chain. For this reason polyvinyl chloride can be partially crystalline even if the polymer is atactic. Atactic polystyrene and atactic polypropylene, however, are always amorphous.

9.9 Crosslinking

Crosslinking is the process by which separate molecules in an amorphous polymer are chemically attached to one another, as shown schematically in Figure 40. The black lines indicate individual molecules, and the red dots represent crosslink points where these molecules are 'tied' together.
image: 73_home_ken_Mydocs_MSEcore_331_figures_Crosslinking.png
Figure 40: Diagram representation of crosslinking.
A representative and commercially important crosslinking system, known by the commercial name of Sylgard 184, consists of the two different parts shown in Figure 41. Part 1 consists of poly(dimethyl silioxane) polymer with vinyl groups (double bonds) at each end of the polymer molecules. Part 2 contains a silicone crosslinker with several silane groups (Si-H) that are able to react with the vinyl groups as shown in Figure 42. The reaction can be viewed as a step growth reaction involving a component with a reactive functionality of 2 (part 1) and a second component with a reactive functionality that is much larger than 2 (part 2).
image: 74_home_ken_Mydocs_MSEcore_331_figures_silicone_divinyl.svg
Figure 41: A diagram of silicone crosslinking.
image: 75_home_ken_Mydocs_MSEcore_331_figures_Silicone_Crosslinking.png
Figure 42: Reaction diagram of silicone crosslinking.

Common Polymers

Here we list some common polymer produced by the different synthesis methods introduced in the previous sections.

10.1 Chain Growth: Addition to a Double Bond

10.1.1 Polyethylene

The simplest polymer from a structural standpoint is polyethylene, with the structure shown below in Figure 43.Dynamicity is not relevant in this case, since there are no substituents other than hydrogen on the carbon backbone.
image: 76_home_ken_Mydocs_MSEcore_331_figures_polyethylene.png
Figure 43: Chemical structure of polyethylene.
High Density Polyethylene:
High density polyethylene refers to version with very little chain branching, thus resulting in a high degree of crystallinity. Completely linear polyethylene has a melting point of 138 °C, and a glass transition temperature near -100 °C.
T g =-120 C, T m =138 C (perfectly linear)
Low Density Polyethylene:
Low density polyethylenes (LDPE's) and high density polyethylenes (HDPE's) are identical in their chemical structure at the atomic level. They are actually structural isomers of one another. Chain branching within low density polyethylene inhibits crystallization, resulting in a material with a melting point lower than 138 degrees C. The decreased crystallinity of LDPE results in a material which is more flexible (lower elastic modulus) than HDPE.
image: 77_home_ken_Mydocs_MSEcore_331_figures_Low_Density_Polyethylene_Zoomed.png
Figure 44: Structure of low density polyethylene
The chain branches responsible for inhibiting crystallization in low density polyethylene are typically short. This illustration shows a 3-carbon (propyl) branch, potentially resulting from intramolecular chain transfer during the polymerization reaction.
T g =-120 C, T m <138 C (depending on branching)
The Importance of Molecular Weight:
Polymeric materials generally have favorable mechanical properties only when the molecular weight is very large - typically hundreds of thousands of g/mol. The point is illustrated with polyethylene:

10.1.2 Polypropylene

image: 78_home_ken_Mydocs_MSEcore_331_figures_polypropylene.png
Figure 45: Structure of polypropylene
The most widely used form of polypropylene is isostatic, with a melting point of 171 °C (for the perfectly isotactic version - a few degrees lower for the actual commercial versions), and a glass transition temperature which is well below room temperature ( 5 C). Single crystals of polypropylene have lower moduli than single crystals of polyethylene along the chain direction, because of the helical structure of propylene. The modulus of isotropic semicrystalline polypropylene is often larger than that of high density polyethylene, however, because of the details of the semicrystalline structure that is formed. The uses of polypropylene and high density polypropylene are similar.

10.1.3 Polybutene-1

This polymer included to illustrate the evolution o the polymer properties when we continue to make the side chain longer.
image: 79_home_ken_Mydocs_MSEcore_331_figures_polybutene-1.png
Figure 46: Polybutene-1

10.1.4 Poly(methyl methacrylate):

Poly(methyl methacrylate) (PMMA) is one of the most common materials used to make polymer glass. It is commonly known by the DuPont tradename Plexiglas,™ and has a glass transition temperature between 100 °C and 125 °C, depending on the tacticity. It is also forms the basis for many biomaterials, including dental adhesives.
T g =100-125 C, no T m (atactic)
image: 80_home_ken_Mydocs_MSEcore_331_figures_polymethylmethacrylate_out.png image: 30_home_ken_Mydocs_MSEcore_331_figures_polymethylmethacrylate.png
Figure 47: Poly(methyl methacrylate)

10.1.5 Poly(methyl acrylate)

Poly(methyl acrylate) is not a widely used polymer, primarily because it's glass transition temperature is too low ( 5 C for the atactic polymer) to be useful as a rigid polymer glass, and too high to be useful as an elastomer. It is included here to illustrate the effect that removing the extra methyl group from the polymer backbone has on the glass transition of the polymer.
T g =5 C, no T m (atactic)
image: 81_home_ken_Mydocs_MSEcore_331_figures_polymethylacrylate.png
Figure 48: Poly(methyl acrylate).

10.1.6 Neoprene

Neoprene, also called polychloroprene, is a material commonly used in wetsuits. Like polyisoprene, it can be polymerized in different forms, corresponding to 1-2, 3-4, cis 1-4 and trans 1-4 addition of the monomer (trans 1-4 addition shown below). This pictured wetsuit has a 0.5mm layer of neoprene sandwiched between layers of nylon and another synthetic material.
image: 82_home_ken_Mydocs_MSEcore_331_figures_neoprene.png image: 27_home_ken_Mydocs_MSEcore_331_figures_polychloroprene.png
Figure 49: Neoprene (polychloroprene).

10.1.7 Polyisobutylene

Polyisobutylene is a common material used to make elastomers, referred to more simply as 'butyl' rubber. It is generally copolymerized by with a small amount of isoprene, so that the resulting double bonds can be used to form a crosslinked material. It is more resistant to solvent penetration than most elastomers, and is often used in applications (like the gloves above) where barrier resistance is needed.
T g =-75 C, T m =2 C
image: 28_home_ken_Mydocs_MSEcore_331_figures_polyisobutylene.png image: 83_home_ken_Mydocs_MSEcore_331_figures_polyisobutylene_gloves.JPG
Figure 50: Polyisobutylene.

10.1.8 Polystyrene

Polystyrene is almost always atactic, and therefore amorphous. It has a glass transition temperature of 100 °C, and is therefore a glassy polymer at room temperature. Its uses are typically in packing material "Styrofoam", and for making cheap plastic objects, like the vials shown above. When suitably modified by the addition of other types of polymers, it is the basis for relatively high performance plastics such as high impact polystyrene (HIPS) and acrylonitrile-butadiene-styrene (ABS).
T g =100 C, no T m (atactic)
image: 84_home_ken_Mydocs_MSEcore_331_figures_polystyrene.png
Figure 51: Structure of polystyrene.

10.1.9 Poly(tetrafluroethylene) (PTFE)Poly(tetrafluoroethylene) (PTFE)

Poly(tetrafluoroethylene) is more commonly known by its DuPont trade name,
Teflon
. Its most outstanding properties are its low surface energy and its low friction against a variety of other materials. It has very poor mechanical properties, however, and is difficult to process by its melting temperature exceeds the temperature at which it begins to degrade.
T g =130 C (by one report[7]); T m 330 C.
image: 85_home_ken_Mydocs_MSEcore_331_figures_polytetrafluorethylene.png image: 34_home_ken_Mydocs_MSEcore_331_figures_polytetrafluoroethylene.png
Figure 52: Monomer unit of poly(tetrafluroethylene) (Image from[2]).

10.1.10 Poly(vinyl acetate)

Poly(vinyl acetate) is often used as base for chewing gum. It is glassy at room temperature but becomes softer at body temperature, which is just above T g .
T g =30 C, no T m (atactic)
image: 86_home_ken_Mydocs_MSEcore_331_figures_polyvinylacetate.png
Figure 53: Monomer unit of poly(vinyl acetate)

10.1.11 Poly(vinyl chloride) PVC

Poly(vinyl chloride) can be partially crystalline even if the material is atactic, because the "R" group in this case is a chlorine atom, which is relatively small. The glass transition temperature of the material is 85 °C, although the addition of small molecules as "plasticizers" can reduce T g to below room temperature. When a material is referred to as "vinyl", it is probably PVC. Record albums (before the age of compact disks) and water pipes are commonly made out of poly(vinyl chloride).
T g =85 C, no T m (atactic)
image: 87_home_ken_Mydocs_MSEcore_331_figures_polyvinylchloride.png Figure 54: Structure of poly(vinyl chloride)

10.1.12 Poly(vinyl pyridine)

Poly(vinyl pyridine) is very similar to polystyrene, and its physical properties (entanglement molecular weight, etc.) are quite similar to the properties of polystyrene. It exists in one of two forms, poly(2-vinyl pyridine) (P2VP and poly(4-vinyl pyridine) (P4VP), based on the location of the nitrogen in the phenyl ring. Both types interact strongly with metals. The polymers are not used in wide quantities, but have been useful in a range of model studies of polymer behavior, often when incorporated with another material as part of a block copolymer.
T g =100 C, no T m (atactic)
image: 88_home_ken_Mydocs_MSEcore_331_figures_polyvinylpyridine.png
Figure 55: Structure of two types of poly(vinyl pyridine)

10.2 Chain Growth: Ring Opening

10.2.1 Poly(ethylene oxide)

Poly(ethylene oxide) (PEO) is generally formed by the ring opening polymerization of ethylene oxide. It is also referred to as polyethylene glycol, although this generally refers to lower molecular weight versions with hydroxyl end groups. PEO is water soluble, and is used in a wide range of biomedical applications, often in a gel form. Lithium salts are also soluble in PEO, and PEO/Li complexes are often used as an electrolyte in battery and fuel cell applications.
T g =-65 C; T m 65 C
image: 89_home_ken_Mydocs_MSEcore_331_figures_peo.png
Figure 56: Monomer unit of poly(ethylene oxide)

10.2.2 Polycaprolactam :

Polycaprolactam is the polyamide equivalent of polycaprolactone, and is synthesized by the ring opening polymerization of the corresponding cyclic amide. It is often referred to as Nylon 6 , since there are 6 carbon atoms in the repeating unit of the polymer. Note that this is different than Nylon 6,6 produced by condensation polymerization, where the repeat unit has 6 carbons originating from each of the two monomers used in the polymerization reaction.
T g 50 C; T m 220 C
image: 90_home_ken_Mydocs_MSEcore_331_figures_polycaprolactam.png
Figure 57: Polycaprolactam.

10.2.3 Polycaprolactone :

Polycaprolactone somewhat unique in that is a polyester that is synthesized by ring opening polymerization of a cyclic ester. It can be viewed as a polyester version of the polyamide, polycaprolactam. Contrary to step growth polymerization of polyesters, the ester linkage is not formed during the polymerization reaction, but is already present in the monomer. Polycaprolactone is biodegradable because the polymer slowly degrades by ester hydrolysis over time.
T g -50 C; T m =60 C
image: 91_home_ken_Mydocs_MSEcore_331_figures_polycaprolactone.png
Figure 58: Polycaprolactone.

10.3 Step Growth Polymers

A variety of common polymers are discussed briefly in the pages below.

10.3.1 Kevlar

Kevlar™ is a trademark of DuPont, Inc. The name actually is used to refer to a variety of aromatic polyamides, or aramids. As the name suggests, the polymers have phenyl groups in the backbone of the chain, and the repeat units are joined by amide linkages. The simplest possible aramid has the structure shown in Figure 59.
T m : above degradation temperature for the polymer.
image: 92_home_ken_Mydocs_MSEcore_331_figures_kevlar.png
Figure 59: Monomer unit of one version of Kevlar

10.3.2 Polycarbonate

A variety of polycarbonates exist. The most common one, (with the GE trademark of Lexan) has a glass transition temperature of 150 °C. It is used for compact disks, eyeglass lenses, and shatterproof glass. See http://www.pslc.ws/macrog/pcsyn.htm for a good description of the synthesis of polycarbonate via a step growth, condensation reaction involving a phenolic di-alcohol and phosgene ( COC l 2 ) .
T g 150 C
image: 93_home_ken_Mydocs_MSEcore_331_figures_polycarbonate_out.png image: 31_home_ken_Mydocs_MSEcore_331_figures_polycarbonate.png
Figure 60: Polycarbonate.

10.3.3 Polyethylene Terephthalate (PET)

Polyethylene terephthalate (trade names include Mylar and Dacron) is produced in fiber form for textiles, and in film form for recyclable bottles, etc. Its degree of crystallinity is highly dependent on the processing conditions, since it can easily be quenched to a glassy state before crystallization is able to occur.
=80 C
=260 C
image: 94_home_ken_Mydocs_MSEcore_331_figures_PET.png
image: 33_home_ken_Mydocs_MSEcore_331_figures_polyethyleneterephthalate_bottle.png image: 95_home_ken_Mydocs_MSEcore_331_figures_polyethyleneterephthalate_dacron.png
Figure 61: Polyethylene Terephthalate

10.3.4 Poly(phenylene oxide)

Polyphenylene oxide is a high performance polymer has many varied uses, largely because of its excellent performance at high temperatures. 190 C
image: 96_home_ken_Mydocs_MSEcore_331_figures_polyphenyleneoxide_out.png image: 32_home_ken_Mydocs_MSEcore_331_figures_polyphenyleneoxide.png
Figure 62: Poly(phenylene oxide).http://plastiquarian.com/ppo.htm

10.3.5 Ultem Polyetherimide

Ultem (a trademark of GE) is a form of polyetherimide. It is a high performance polymer that combines high strength and rigidity at elevated temperatures with long term heat resistance ( T g =215 C). The repeat unit is illustrative of the complex chemical structure of many modern, high performance polymers.
image: 97_home_ken_Mydocs_MSEcore_331_figures_Ultem_Polyetherimide_1.png image: 98_home_ken_Mydocs_MSEcore_331_figures_Ultem_Polyetherimide_2.png
Figure 63: Structure of Ultem polyetherimide, with some materials that have been made from it.http://www.alcanairex.com/products/e/100/120p01_e.htm

10.3.6 Silicones

Silicones are an important class of synthetic polymers which do not have carbon in the backbone. Instead, the backbone consists of alternating silicon and oxygen atoms. Different classes of silicones are specified by the substituents on the silicone atoms. Poly(dimethyl siloxane) (PDMS), with methyl substituents, is the most important silicone. Its glass transition and melting temperatures are very low, so that it remains flexible at very low temperatures. It also has a very low surface energy, and forms a hydrophobic surface that is very water resistant.
=-130 ºC ; =-45ºC
image: 99_home_ken_Mydocs_MSEcore_331_figures_pdms_structure.png image: 29_home_ken_Mydocs_MSEcore_331_figures_Silicone_2.png
Figure 64: Poly(dimethyl siloxane).
We list PDMS here as a step growth polymer because it can be produced from a self-condensation of silanol (SiOH) groups. The starting point is actually dimethyl-dichlorosilane. In th presence of water the SiCl bonds hydrolyze to SiOH:
image: 100_home_ken_Mydocs_MSEcore_331_figures_chlorosilane_hydrolysis.png
Figure 65: Hydrolysis of chlorosilane bonds.
The resultant silanol groups than can then condense by the elimination of water:
image: 101_home_ken_Mydocs_MSEcore_331_figures_silanol_condensation.png
Figure 66: Condensation of silanol groups.
One of the interesting features of silicones is that they can also be synthesized by anionic, ring opening polymerization of cyclic, oligomeric forms of PDMS. Here's one example:
image: 102_home_ken_Mydocs_MSEcore_331_figures_silione_anionic_polym.png
Figure 67: Anionic ring opening polymerization of a cyclic silicone oligomer.

10.4 Copolymers

The polymers considered so far are homopolymers, which means they consist of a single chemical repeat unit. Copolymers consist of more than one repeat unit. Examples include
statistical copolymers
, where the different repeating units are distributed throughout a molecule, as shown schematically in Figure 68.
image: 103_home_ken_Mydocs_MSEcore_331_figures_Copolymers_1.png
Figure 68: Schematic representation of a statistical copolymer.
In
block copolymers
, the different repeating units appear as distinct 'blocks'. A diblock copolymer for example, can be viewed as two different homopolymer molecules that have been covalently joined to one another as shown in Figure 69. A wide variety of interesting block copolymer morphologies have been observed. Diblock copolymers are most often synthesized by living, chain growth polymerizations.
image: 104_home_ken_Mydocs_MSEcore_331_figures_Copolymers_2.png
Figure 69: Schematic representation of a block copolymer.
Exercise:
Describe a procedure for making a diblock copolymer of polystyrene and poly(methyl methacrylate).
Solution:
The details of how to do this go a bit beyond the scope of this text, but conceptually we know we need to do the following in sequence: Add and appropriate initiator for a living chain growth polymerization, add monomer 'A' until all of it reacts, add monomer B until all of it reacts, terminate the reaction. Here are the details that can work for this system.
  1. Add an appropriate initiator. We don't want to use a radical initiator, because we need to avoid the inevitable termination reactions. From the information given in this text, either anionic or cationic polymerization. It turns out that the only anionic polymerization works in this case. We'll choose tetrahydrofuran as a solvent, clean it up very well to get rid of any impurities, and add some secondary butyllithium (see Fig. 31) as the initiator. The molar quantity of initiator will be equal to the number of moles of polymer we get at the end of the reaction.
  2. Add the first monomer to the initiator solution. The order matters in this case. We need to add the styrene first, and choose reaction conditions so that the initiation is much faster than the subsequent polymerization. This ensures that the polymers all have roughly the same length, since they all started polymerizing a the same time.
  3. Once all the styrene reacts, we add the methyl methacyrylate monomer to the solution. Because the polymerization is 'living', the active chain ends from the polystyrene part of the polymerization are still available to react with the methyl methacrylate.
  4. Terminate the reaction. This is very easy to do, and can be done by adding a small amount of methanol to the reaction, to produce lithium methoxide and a 'dead' chain after the proton from alcohol trades places with the lithium that is associated with the negatively charged carbon at the end of the polymer chain.

11 The Glass Transition

The information presented so far about the structure of amorphous polymers applies both to elastomers like natural rubber, in addition to engineering thermoplastics like polystyrene. The properties of rigid plastics and rubber are obviously quite different, in spite of these similarities. The difference is that below a certain temperature, molecular motion is no longer possible, and the material becomes a relatively hard solid. This transition between “soft” and “hard” behavior occurs at the
glass transition temperature
, T g .

11.1 Free volume

A very useful way to think about the glass transition involves the concept of unoccupied space, or 'free volume' in a polymer. This concept is illustrated conceptually in Figure 70. The red dots connected together by lines represent polymer molecules. In this lattice model of polymers, the red dots are constrained to fit on a lattice. (The lattice is completely artificial, and does not exist in reality). These dots represent the '
occupied volume
' . Empty lattice sites represent
free volume
.
The lattice on the left has no free volume, so there is no way for the molecules to move by the 'hopping' of segments into a small region of unoccupied space. The lattice on the right does have free volume, however, so that molecular motion is relatively easy.
image: 105_home_ken_Mydocs_MSEcore_331_figures_Glass_Transition.png
Figure 70: Conceptual representations of polymers with and without free volume.
The conceptual free volume model of the glass transition makes the assumption that
volume thermal expansion coefficient
,
α is defined in the following manner:
α = 1 V dV dT (11.1)
Note that the linear thermal expansion coefficient, defined in terms of the linear dimensions of the sample rather than the volume, is equal to α/3. If the fractional change in free volume is small ( Δ V/V1 ) then α describes the linear relationship between the change in volume in the change in temperature:
Δ V V = α Δ T (11.2)
In the glassy phase ( T< T g ), the free volume is assumed to remain constant, so that increases in the volume are governed entirely by increases in the occupied volume:
Δ V occ V = α Δ T (11.3)
In the liquid phase ( T> T g ), the free volume and occupied volumes both increase with temperature, and we have:
Δ V V = Δ V f V + Δ V occ V = α Δ T (11.4)
Combining11.3 and gives the following expression for the free volume in the liquid state:
Δ V f V =( α - α g ) Δ T (11.5)
Now we define a new temperature, T which is the theoretical temperature at which the free volume would go to zero, if the material were to continue to follow the liquid behavior below T g . It is illustrated schematically in Figure 71. To determine the fractional free volume at the glass transition, we use Eq. 11.5 with Δ T= T g - T , recognizing that V f /V=0 at T= T :
Typically, T is about 50K below the measured glass transition temperature, and the fractional free volume at T g is in the range of a couple percent.
image: 106_home_ken_Mydocs_MSEcore_331_figures_Glass_Transition_Drawing.png
Figure 71: The glass transition and the free volume concept.
The concept of free volume can be helpful in sorting out how the glass transition depends on the structure of the polymer, based on the following two guidelines:
  1. Changes to the polymer structure that increase the free volume needed in and r for the backbone of the polymer to move will increase T g . Making the backbone very still (by the incorporation of phenyl groups, as in increase T g . Also, bulky, 'fat' substituents (like methyl or phenyl groups) tend to increase T g when added to the backbone.
  2. Changes to the polymer structure that introduce extra free volume will decrease T g . This is the case for long, thin additions to the backbone, such as alkyl chains ( - ( C H 2 ) n - ) .
Without experience it's difficult to know which of these effects is greater. If I add something to the backbone for example, does it create more free volume that it needs to move? The best way to get a handle on this is to take a look at the glass transition temperatures of atactic polymers with the general chemical structure shown in Figure 72. The glass transition temperatures for these different polymers are listed in Table 1.
image: 107_home_ken_Mydocs_MSEcore_331_figures_repeat_with_R1R2.png
Figure 72: Structure of the repeat units for the materials listed in Table1.
R 1 =H
R 1 =CH 3
R 2 : H
-125
-13
R 2 : CH 3
-13
-75
R 2 :image: 108_home_ken_Mydocs_MSEcore_331_figures_phenyl_side_group.png
100
180
R 2 :image: 109_home_ken_Mydocs_MSEcore_331_figures_methyl_ester.png
10
105
R 2 :image: 110_home_ken_Mydocs_MSEcore_331_figures_butyl_ester.png
-54
20
Table 1: Glass transition temperatures (in C) for polymers with the general structure shown in Figure 72 (From ref.[4]).
Exercise:
Polystyrene has a volume thermal expansion coefficient of 6 × 1 0 -4 K -1 for temperatures above T g , and a thermal expansion coefficient of 2 × 1 0 -4 K -1 for temperatures below T g . Estimate the ratio of the free volume to the total volume of sample in the glassy state. Assume that T is 50K below T g .
Solution:
Since the fractional free volume is assumed to be fixed below T g ,we are looking for the for the fractional free volume at the glass transition, as given by Eq. 11.5. With α =6 × 1 0 -4 K -1 , a g =2 × 1 0 -4 K -1 and T g - T =50 K, we obtain:
V f V =( 4 × 1 0 -4 K -1 ) ( 50 K ) =0.02
This value of 0.02 is a typical value for the fractional free volume in a glassy polymer.

11.2 Enthalpy and Heat Capacity

image: 111_home_ken_Mydocs_MSEcore_331_figures_Enthalpy_and_Heat_Capacity_Diagram.png
Figure 73: A graph of enthalpy vs temperature
The heat capacity ( c p ) of the liquid is larger than the heat capacity of the glass. Note that the enthalpy itself is continuous at T= T g . The glass transition is therefore a second order transition - thermodynamic quantities like volume, enthalpy and entropy are continuous at the transition but the derivatives of these quantities with respect to temperature are discontinuous.

11.3 Differential Scanning Calorimetry (DSC)

Glass transition temperatures are most commonly measured by differential scanning calorimetry . The technique can be used to measure the glass transition temperature of a polymer, in addition to the temperature at which a polymer melts or crystallizes. The heat of fusion (heat required to melt a semicrystalline polymer) can also be obtained. Heat flowing into the sample is endothermic, and heat flowing out of the sample is exothermic. We have plotted endothermic heat flow in the positive direction, and exothermic heat flow in the negative direction. Unfortunately, there is no universally followed sign convention for displaying DSC data. Sometimes endothermic heats are plotted in the negative direction. In order to avoid confusion, the endothermic direction is commonly indicated on the DSC plot, as we have done in the curve shown in Figure 75. Students are encouraged to look at the Macrogalleria web sites on the glass transition[3] and differential scanning calorimetry[1] for discussions at the appropriate level for this class.
image: 112_home_ken_Mydocs_MSEcore_331_figures_Differential_Scanning_Calorimetry_1.png
Figure 74: Simple schematic representation of a calorimeter.
image: 113_home_ken_Mydocs_MSEcore_331_figures_Differential_Scanning_Calorimetry_2.png
Figure 75: Graph of heat flow vs temperature.

Structure of Amorphous Polymers

Materials science is the study and application of the relationships between the properties of a material, the structure required to obtain these properties, and the processing methods which can be used to obtain these properties. Our discussion of polymer synthesis was motivated primarily by a need to understand processing methods for polymer materials. We also introduced the structure of different polymers at the atomic level, corresponding to the arrangements of individual atoms in the polymeric repeat units. Our discussion now moves to the structure of polymers on a molecular scale. We begin with amorphous polymers, and follow with a discussion of semicrystalline polymers. Ultimately, we will find that our understanding of many of the important properties of polymers can be related to these structural features by simple, yet remarkably accurate theories. In particular, the theory of rubber elasticity relates the elastic properties of an elastomer to the molecular structure described here.
The squiggly lines in Figure 76 represent the backbones of individual molecules in an amorphous polymer. (One molecule is highlighted in red) Our goal is to relate the properties of an amorphous polymer to the distribution of shapes of these molecules. Ultimately, we will show that many of the mechanical properties of amorphous polymers at temperatures above the glass transition are related to this distribution of polymer shapes.
image: 114_home_ken_Mydocs_MSEcore_331_figures_Amorphous_Polymers.png
Figure 76: A diagram of an amorphous polymer.

12.1 End-to-end Vector for Polymer Molecules

image: 115_home_ken_Mydocs_MSEcore_331_figures_end-to-end-vector.png
Figure 77: The end-to-end vector, R , for a polymer molecule.
The shape of this particular molecule can only be completely specified by describing the path taken by the molecule from one end to the other. Fortunately, the most important quantity is much simpler. This quantity
end-to-end vector
, R which is simply the vector spanning the two ends of a given polymer molecule. This vector will in general be different for different polymer molecules in a sample, but the distribution of vectors can be accurately predicted. The shapes of amorphous polymer molecules are random, and random walk statistics can be used to describe the distributions of end-to-end vectors that are obtained.

12.2 Random Walk in One Dimension

We can understand much of what we need to understand by considering the simplest case of a random walk in 1 dimension. In this case our 'walker' takes N x steps, each of which is a randomly chosen move to either the left or right. At the end of these N x steps we write down how many steps removed from the starting point our walker is. We refer to this number as i . For example, if the walker ends up 4 steps to the right of where he ended up, then i=4. Where our walker ended up after one particular trip is not very important though. We need to have him take lots of trips, and collect statistics to see where he ended up after taking a large number of trips. We'll let him take a total of M trips, and write down the number of times, n, that the walker has ended up i steps from where we started. If the walk is really random, then the relationship between n and i will be approximated by Eq. 12.1
image: 116_home_ken_Mydocs_MSEcore_331_figures_Random_Walk_Animation_2.png
Figure 78: Random walk animation
n( i ) =2M ( 1 2 π N x ) 1/2 exp ( - i 2 2 N x ) (12.1)
We can check that the summation is appropriately normalized. In other words, is the following expression valid?
i n i =M (12.2)
The summation here is over all possible values of i . Note that for a random walk with an even number of steps, only even values of i are possible. Similarly, for a walk with an odd number of steps, only odd values of i are possible. We can get the summation we want by summing over all integer values of i and then dividing by two:
in( i ) = M 2 π N x i=- exp ( - i 2 2 N x ) M 2 π N x - exp ( - i 2 2 N x ) di (12.3)
Functions of the form f( x ) exp - ( x/A ) 2 (where A is a constant) are called Gaussian functions, and generally describe random processes. The following integrals of a Gaussian function will be very useful for us:
- exp - ( x/A ) 2 dx=A π ( a ) - x exp - ( x/A ) 2 dx=0 ( b ) - x 2 exp - ( x/A ) 2 dx= A 3 π /2 ( c ) (12.4)
In our case, A= 2 N x , so i n( i ) =M and our expression for the Gaussian distribution of random walks is indeed normalized.

12.3 Average of a Function

The entire distribution represented by the Gaussian function will be useful to us, but it is still useful to have some averages. We encountered averages already in our discussion of the number and weight average distribution. In general, if P x ( x ) is the probability that a discrete variable has a value of x , than the average value of x (referred to as x ) is obtained from the following expression:
x =xx P x ( x ) (12.5)
where the sum is over all possible values of x . Similarly, the average of x 2 is given by: x 2 =x x 2 P x ( x ) dx (12.6)
The procedure can be generalized to calculate the average of any function of x :
f( x ) =xf( x ) P x ( x ) (12.7)
Exercise:
Suppose I put the following 10 numbers into a drawer: 1, 4, 8, 12, 19, 25, 28, 33, 37, 45. I randomly pick a number from the drawer and then return it. What is the average value of all the numbers that I pick if I continue with this exercise?
Solution:
We just use Eq. 12.5. Since there are 10 numbers and I am equally likely to be picked, they each have a value of 0.1 for P r . So in this case we have:
=0.1( 1+4+8+12+19+25+28+33+37+45 ) =21.2
In this case the answer is just the average of all the numbers, which is about as simple it gets.
In general for a function of i, where i is only able to take on discrete values, we have:
<f( i ) >=if( i ) P r ( i )
What if we have a function of a continuous variable, x ? In this case, we must replace P(i), with P r ( x ) dx , where P r ( x ) dx is the probability that the continuous variable has a value between x and x+dx . We must also replace the summation with the appropriate integral to obtain the following result:
f( x ) = f( x ) P r ( x ) x (12.8)
Exercise:
Suppose I put the following 10 numbers into a drawer: 1, 4, 8, 12, 19, 25, 28, 33, 37, 45. I randomly pick a number from the drawer and then return it. What is the average value of all the numbers that I pick if I continue with this exercise?
Solution:
Calculate the average value of x 2 , assuming that x is equally likely to take on all values between 1 and 10, and that no values outside this range are possible. We know that P r has some constant value between 1 and 10. In order for our procedure to work, we need to make sure that the probability distribution is normalized, with x min x max P r ( x ) x =1. In our case this gives P r =1/9 and we get:
x 2 = 1 9 1 10 x 2 x = 1 27 . x 3 1 9 = 1 0 3 -1 27 = 999 27 =37

12.4 Averages for Random Walks

If the length of each step taken by the random walker is a , then the distance between the beginning and the end of the walk is simply i times this step length:
R x =ia
We can substitute R x /a for i and d R x /a for di in Eq. 12.3 to get an expression for the total number of walks M. We'll generalize a bit further and replace the limits of - and for i to limits of R min and R max for R x . We obtain the following for M( R min , R max ) , the total number of walks with values of R x between R min and R max :
M( R min , R max ) =M R min R max P x d R x (12.9)
with the probability density function, P x given as follows:
P x = 1 2 π N x a 2 exp ( - R x 2 2 N x a 2 ) (12.10)
To calculate the average values of R x and R x 2 , referred to as R x and R x 2 , respectively, we use Eq. 12.8 to calculate the averages:
R x = - R x P x d R x = 1 2 π N x a 2 - R x exp ( - R x 2 2 N x a 2 ) =0 (12.11)
R x 2 = - R x 2 P x d R x = 1 2 π N x a 2 - R x 2 exp ( - R x 2 2 N x a 2 ) = N x a 2 (12.12)
We know intuitively that Eq. 12.11 must be true because the Gaussian distribution function is symmetric: positive and negative values of R x with equal magnitudes are equally likely. This means that any contribution to R x will be exactly offset by a contribution with the opposite sign. For this reason, the root mean square end-to-end distance, R x 2 = N x a , is typically used to describe the spatial extent of a random walk.

12.5 Random Walks in 3 Dimensions

Now we can extend some of these results to describe the shapes of real polymer molecules. We begin with a discussion of two dimensional random walks. The simplest extension is to give our random walker 4 choices for his movements. In addition to moving to the left and right, he can now move up and down, as illustrated in Figure 79a. If we move the constraint that there are only 2 directions in which the walker can move, we get the example shown in Figure 79b. These trajectories are now beginning to look like the shapes of actual polymer molecules, which is why this whole exercise is relevant.
image: 117_home_ken_Mydocs_MSEcore_331_figures_randwomwalks2d.png
Figure 79: 2-D Random Walk of 1000 steps on a square lattice (a) and with random bond angles (b). The end-to-end vector, R , is shown in each case.
Much of our understanding of amorphous polymers is based on our treatment of individual molecules as random walks in three dimensions. We assume that we can treat a polymer of degree of polymerization N as a random walk of N steps, each with an effective length of a, which we refer to as the
Statistical segment length
.
The statistical segment length does not correspond to a real polymer dimension that we can get directly by looking at the structure of the polymer molecule. We need more detailed models of chain dimensions for that. For now, you just need to know that it is a length that gives the right overall molecular size. To simplify our treatment and make the connection back to the statistics that we already developed for one-dimensional random walks, we make the following simplifying assumptions:
The value of R x that we studied in the one dimensional random walk case is just the x component of the end-to-end vector, which also has y and z components:
R = R x x ˆ + R y y ˆ + R z z ˆ (12.13)
where x ˆ , y ˆ and z ˆ are unit vectors in the x, y and z directions, respectively. The probability, P xyz that a given random walk of N steps in three dimensions has an end-to end vector of with x, y and z components of R x , R y and R z is given by multiplying the probabilities for the individual components:
P xyz = P x P y P z = R 0 3 ( 3 2 π ) 3/2 exp ( -3 R 2 2 R 0 2 ) (12.14)
Here we have assumed that P y and P z have the same form as P x from Eq. 12.10, with y or z substituting for x as appropriate. We have also used the following relationships:
R 2 = R x 2 + R y 2 + R z 2 , (12.15)
N x = N y = N z =N/3, (12.16)
along with the following definition of R 0 :
R 0 = N a (12.17)
Recall that P xyz ( R x , R y , R z ) d R x d R y d R z is the probability that x component of the end-to-end vector is between R x and R x +dx , the y component is between R y and R y +dy , and the z component is between R z and R z +dz . We can also define a probability density P R , where P R dR is the probability that the magnitude of the end-to-end vector lies between R and R+dR. Because the probability only depends on R for an isotropic system (as must be the case for a truly random walk), we can substitute 4 π R 2 dR for d R x d R y d R z to obtain:
P R =4 π R 2 P xyz = 4.15 R 2 R 0 3 exp ( -3 R 2 2 R 0 2 ) (12.18)
Note that R 0 2 is the average of R 2 , obtained from the following version of Eq. 12.8:
R 2 = 0 R 2 P R dR= 4.15 R 0 3 0 R 4 exp ( -3 R 2 2 R 0 2 ) = R 0 2 (12.19)
image: 118_home_ken_Mydocs_MSEcore_331_figures_Probability_Density.png
Figure 80: Probability density.
Note that the distribution is normalized: 0 P r d R =1 .
Exercise:
What fraction of molecules in an amorphous polymer have end to end vectors that are within 10% of R 0 ?
Solution:
P R =4.15 R 2 exp ( -1.5 R 2 )
Note that this expression for P R is normalized, with 0 P R R =1. We need to integrate the probability density from R =0.9 to R = 1.1:
P r tot =4.150.91.1 R 2 exp ( -1.5 R 2 ) d R
This corresponds to the shaded area in the following figure:
image: 119_home_ken_Mydocs_MSEcore_331_figures_Probability_Density_Practice.png
Numerical integration of the equation above gives P r tot = 0.18, which corresponds to the red shaded area. This is the total fraction of molecules with R between 0.9 R 0 and 1.1 R 0 .

13 Models of Chain Dimensions

Our discussion of random walk statistics tells us what the shapes of the molecules look like when we know the statistical segment length, but it doesn't tell us what the statistical segment length actually is for a particular polymer. Clearly the statistical segment length must be related to the detailed chemical structure of the repeat unit itself. For example, we wouldn't expect that polyvinyl chloride and polyethylene terephthalate would have the same statistical segment length, based on the very different chemical structures of these polymers. To obtain a prediction for the statistical segment length we need to zoom in and consider the actual backbone bonds that are connecting the polymer together. In general every bond is vector, r with a length and an orientation. As illustrated in Figure 81, the end-to-end vector is obtained by adding up all N b bond vectors ( N b =12 in this example):
R = i=1 N b r i = r 1 + r 2 + r 3 +...+ r N b (13.1)
image: 120_home_ken_Mydocs_MSEcore_331_figures_Quantitative_Chain_Dimension_Calculations.png
Figure 81: Vector diagram of a polymer.

13.1 General Considerations

We are interested in the average magnitude of the end-to-end vector. As discussed in the previous section, the average of R itself is not useful, since it just averages to zero. We are interested in the average value of the square of the magnitude of the end-to end vector. This quantity is equal to the dot product of the end-to-end vector with itself:
R 2 = R R =( r 1 + r 2 + r 3 +...+ r N b ) ( r 1 + r 2 + r 3 +...+ r N b ) (13.2)
The dot product has N b 2 terms, which can be represented in the following matrix form:
R 2 ={ r 1 r 1 + r 1 r 2 + r 1 r 3 +... + r 1 r N b + r 2 r 1 + r 2 r 2 + r 2 r 3 +... + r 2 r N b + r 3 r 1 + r 3 r 2 + r 3 r 3 +... + r 3 r N b +... +... +... +... +... + r N b r 1 + r N b r 2 + r N b r 3 +... + r N b r N b } (13.3)
The dot product between two vectors, A and B is given in general by | A | | B | cos θ where θ is the angle between the two vectors, illustrated below in Figure 82. At this point we will make our first simplifying assumption, which is that the lengths of all bonds along the backbone are identical. This assumption is certainly valid for vinyl polymers, and other polymers which have only C-C single bonds along the backbone.
image: 121_home_ken_Mydocs_MSEcore_331_figures_Dot_Product.png
Figure 82: Definition of the angle between two vectors.
If this bond length is , then R 2 can be represented as follows, where θ ij is the angle between bondi and bondj:
R 2 = 2 { cos θ 11 + cos θ 12 + cos θ 13 +... + cos θ 1 N b + cos θ 21 + cos θ 22 + cos θ 23 +... + cos θ 2 N b + cos θ 31 + cos θ 32 + cos θ 33 +... + cos θ 3 N b +... +... +... +... +... + cos θ N b 1 + cos θ N b 2 + cos θ N b 3 +... + cos θ N b N b } (13.4)
The average value of R, will be determined by the average values of θ ij In mathematical terms, we have the following expression:
R 2 = 2 { cos θ 11 + cos θ 12 + cos θ 13 +... + cos θ 1 N b + cos θ 21 + cos θ 22 + cos θ 23 +... + cos θ 2 N b + cos θ 31 + cos θ 32 + cos θ 33 +... + cos θ 3 N b +... +... +... +... +... + cos θ N b 1 + cos θ N b 2 + cos θ N b 3 +... + cos θ N b N b } (13.5)
In more compact notation we can write:
R 2 = 2 i,j=1 N b cos θ ij (13.6)
So we have reduced the problem to figuring out what the values of cos θ ij are in the system. This requires a specific model of how the bonds are joined to one another. In the following sections we consider two specific models (the freely jointed and freely rotating models), and a general model where the bond rotation angles are constrained.

13.2 Freely Jointed Chain Model

The simplest (and most unrealistic) model of chain dimensions is the freely jointed chain model, where all bond orientations are equally likely. In the freely jointed chain model cos θ ij =0 for ij , because orientations giving positive and negative values of cos θ ij are equally likely. The only exception is for the N b diagonal terms for which i=j . These represent the projection of a vector onto itself, and all have θ =0 and cos θ =1 .
R 2 = N b 2 (13.7)

13.3 Freely Rotating Chain Model

In this model, bond angles are fixed, but rotation about bonds is possible. Rotation about a bond sweeps out a cone as shown below. The bond rotation angle is defined as φ , and in the freely rotating model all values of φ are assumed to be equally likely.
image: 122_home_ken_Mydocs_MSEcore_331_figures_Free_Rotating_Chain_Model.png
Figure 83: Freely Rotating Chain Models with views perpendicular to the axis of the rotating bond and along the axis of the rotating bond
In the freely rotating chain model, the angle between adjacent bonds is fixed at certain angle, referred to here simply as θ (without the subscripts). Adjacent bonds,i.e. those with i-j=1 or i-j=-1 , have cos θ ij = cos θ . It can be shown that terms with i-j =2 or -2 have cos θ ij = cos 2 θ , terms with i-j =3 or -3 have cos θ ij = cos θ , etc. The matrix of terms making up the summation to give R 2 therefore has the following form:
R 2 = 2 { 1 + cos θ + cos 2 θ +... + cos ( N b -1 ) θ + cos θ +1 + cos θ +... + cos ( N b -2 ) θ + cos 2 θ + cos θ +1 +... + cos ( N b -3 ) θ +... +... +... +... +... + cos ( N b -1 ) θ + cos ( N b -2 ) θ + cos ( N b -3 ) θ +... +1 } (13.8)
Of the N b 2 terms in the summation, N b terms (the diagonal terms with i=j ) are equal to 1. There are N b -1 terms with i-j =1, and an additional N b -1 terms with i-j =-1. There are therefore 2( N b -1) terms with a magnitude equal to cos θ . Similarly, there are 2( N b -2) terms equal to cos 2 θ , 2( N b -3) terms equal to cos 3 θ , etc. The expression for R 2 can therefore be written as follows:
R 2 = 2 { N b +2( N b -1 ) cos θ +2( N b -2 ) ( cos θ ) 2 +...+2 ( cos θ ) N b -1 }
In the more compact notation, we have:
R 2 = N b 2 + 2 i=1 N b -1 2( N b -i ) ( cos θ ) i (13.9)
Equation 13.9 is an exact expression that is valid for any value of N b . Unfortunately, it's not a very useful equation. Fortunately, we can get a greatly simplified expression that is very nearly exact for any reasonably large value of N b . We are able to do this because cos θ is less than one, so that ( cos θ ) i rapidly decreases asi increases. For example, for θ =7 1 (the value corresponding to C-C single bonds), ( cos θ ) i = 1.3x1 0 -5 for i=10 . For large values of N b , we only need to consider contributions from values of i which are much smaller than N b . A very good approximation is obtained by substituting N b for N b -i and extending the sum to i= , which leads to the following:
R 2 = N b 2 { 1+2 i=1 ( cos θ ) i } (13.10)
At this point we are in a position to use the following standard sum, valid for any value of x with an absolute magnitude less than 1:
i=1 x i = x 1-x (13.11)
Use of this expression with x= cos θ gives (after a little algebraic rearrangement):
R 2 = N b 2 { 1+ cos θ 1- cos θ } (13.12)
For singly bonded carbon along the backbone, θ =7 1 and R 2 =2 N b 2 .

13.4 Bond Angle Restrictions

The freely rotating model is overly simplistic, because all bond rotation angles cannot be equally likely. At certain angles, functional groups attached to adjacent carbon atoms are too close to one another, as illustrated here.
image: 123_home_ken_Mydocs_MSEcore_331_figures_Restricted_Chain_Model.png
Figure 84: Bond angle restrictions shown perpendicular and along the axis of rotation (Newman projection ).
Interferences between functional groups cause the free energy to vary as a function of the bond rotation angle, φ , with an energy vs. φ curve that looks like the one shown in Figure 85. The
trans configuration
, where φ =0 , has the lowest free energy, because the backbone carbon bonds are as far from each other as possible in this configuration. Secondary minima in the energy exist for φ = ± 12 0 , where there is minimal overlap between the additional groups bonded to the carbon carbon molecules. These configurations are referred to as the
gauche configurations
.
image: 124_home_ken_Mydocs_MSEcore_331_figures_Karplus_Curve.png
Figure 85: Energy as a function of bond rotation angle.
We know that the freely rotating chain model must also be overly simplistic, because some bond rotation angles (those corresponding to the gauche and trans configurations, for example) are more likely than others. In this case R 2 is given by the following expression, which we present here without proof.
R 2 = N b 2 { 1+ cos θ 1- cos θ } { 1+ cos φ 1- cos φ } (13.13)
Here cos φ is the average of cos φ , taking into account the fact that not all angles are equally likely. Note that cos φ =0 for the freely rotating chain model, in which case Eq. 13.13 reduces to Eq. 13.12.

13.5 Characteristic Ratio

The
characteristic ratio
, C , is defined as the ratio of the actual value of R 2 to the value which would be obtained from the freely jointed chain model: C R 2 N b 2 (13.14)
Each model of chain dimensions has it's own prediction for the value of C . For a chain with fixed bond angles, comparison to Eq. 13.13 gives the following:
C ={ 1+ cos θ 1- cos θ } { 1+ cos φ 1- cos φ } (13.15)
Exercise:
What would the characteristic ratio be for polystyrene if 70% of the bonds were in trans configurations and the remaining 30% were in gauche configurations?
Solution:
We use Eq. 13.15 for C , with θ =7 1 , the correct bond angle for singly bonded carbon, which we have along the polystyrene backbone. Gauche bonds have φ = ± 12 0 ( cos φ =-0.5) trans bonds have φ =0 ( cos φ =1) . With our assumed distribution of gauche and trans bonds we obtain the following for cos φ :
cos φ =0.7( 1 ) +0.3( -0.5 ) =0.55
This gives C =6.9 , which is actually a pretty reasonable value. Most polymers with singly bonded carbons along the backbone have a value for C that is between 6 and 10.
Exercise:
How is the statistical segment length of a vinyl polymer related to C ?
Solution:
The characteristic ratio and statistical segment length are both defined in terms of an expression for R2. All we need to do is equate these two expressions (Eqs.12.17 and13.14):
C N b 2 =N a 2 (13.16)
The repeat unit for a vinyl polymer contains two C-C single bonds along the backbone, so N b =2N . Also, single bonds between carbon atoms have a bond length of 1.54 Å. The statistical segment length is therefore related to C as follows:
a=2.18 Å C (13.17)

13.6 Self-Avoiding Random Walks

image: 125_home_ken_Mydocs_MSEcore_331_figures_Random_Walk_Self-Avoid.png
Figure 86: Random walk diagram
An isolated polymer molecule in a solvent no longer obeys random walk statistics. If the polymer is surrounded by solvent molecules, these surroundings no longer look the same as other parts of the same molecule. The chain will adopt shapes that avoid direct polymer/polymer contacts, and the end-to-end distance will be larger than the value predicted by random walk models. The rms end-to-end distance increases with the 0.6 power of the chain length, instead of the 0.5 power that was obtained for random walk statistics:
R 0 = a s N 0.6 (13.18)
Note: The value of a s is not exactly the same as the statistical segment length describing chain dimensions in pure polymer melts, but it will be close to it.

Elasticity Basics

In this section, we discuss the mechanical properties of polymers, defined as the relationships which are observed between an applied stress state and the resulting strain. Are focus in this course is on elasticity and viscoelasticity. We begin with a refresher on the fundamental definitions of stress and strain, and then discuss the three most important experimental geometries for mechanical measurements of soft materials.

14.1 Specification of the State of Stress

The two dimensional stress state is illustrated below in Figure 87. In two dimensions there are 3 independent components of the stress tensor. Two of these are normal forces where the stress is perpendicular to the relevant surface. These are σ xx and σ yy , with forces acting in the x and y directions, respectively. There is also a shear stress, σ xy , which act parallel to each of the 4 surfaces of a square chunk of material. Note that the magnitudes of the shear forces acting on the 4 different surfaces are all equal to one another, a fact imposed by the requirement that there must be no net force or torques acting on this overall square of material (otherwise the material would be accelerating in response to this net force). The values of these three stress components depend on the way we define the orientation of the coordinate system. It is always possible to define a coordinate system so that the shear components disappear, and we are only left with normal stresses. These are the
principal stresses
, which in two dimensions are σ 1 p and σ 2 p . The same holds for a 3-dimensional stress state: we can always find a coordinate system for any given stress state so that only normal stresses. In three dimensions we have three principal stresses, σ 1 p , σ 2 p and σ 3 p . These act in the three, mutually perpendicular principal stress directions as illustrated in Figure 88. In summary, any state of stress is completely specified by the three principal stresses ( σ 1 p , σ 2 p , σ 3 p ), and the orientation of the coordinate system used to define these principal stresses.
image: 126_home_ken_Mydocs_MSEcore_331_figures_stress_state_2d.png
Figure 87: Representation of a two-dimensional stress state.
image: 127_home_ken_Mydocs_MSEcore_331_figures_principal_stresses.png
Figure 88: 3-D representation of stress, expressed in terms of the principal stresses, σ 1 p , σ 2 p and σ 3 p .

14.2 Extension Ratios and the Strain Ellipsoid

A similar concept of principal axes applies to the strain state as well. For polymeric systems (and elastomers in particular), we need a description that can be used to describe both small and large strains. We the concept of the
strain ellipsoid
, illustrated in Figure 89. An undeformed sphere of radius r 0 is deformed into an ellipsoid with principal radii of r 0 λ 1 , r 0 λ 2 and r 0 λ 3 . Here λ 1 , λ 2 and λ 3 are the
principal extension ratios
. The values of these principal extension ratios, along with the orientation of the principal axes of the strain ellipsoid, completely specify the state of strain.
image: 128_home_ken_Mydocs_MSEcore_331_figures_strain_ellipsoid.png
Figure 89: Deformation of a sphere of radius r 0 into an ellipsoid with principal radii of r 0 λ 1 , r 0 λ 2 and r 0 λ 3 .

14.3 Uniaxial Extension or Compression

A common mechanical testing geometry corresponds to the application of a single normal force to the sample, as illustrated in Figure 90. For our purposes we define the coordinate system so that the normal force is applied in the z direction, and we refer to this force as P z . We also assume that our sample has isotropic mechanical properties and is in the shape of a cylinder. The result of the applied force is to increase the length of of the sample from 0 to , and to decrease it's diameter from d 0 to d .
image: 129_home_ken_Mydocs_MSEcore_331_figures_uniaxial_extension.png
Figure 90: Uniaxial tensile deformation.
According to the first law of thermodynamics , the change in free energy of a system is equal to the work done on the system, provided that the process is reversible. In this example, an applied deformation force, P z , produces a small displacement, d . The work done on the system is P z d , which is equal to the change in the deformation free energy change, Δ F d :
P z d=d( Δ F d ) (14.1)
We can write this force in terms of λ z , the extension ratio in this direction:
P z = d( Δ F d ) d = d( Δ F d ) d λ z d λ z d (14.2)
With λ z =/ 0 , we have d λ z d =1/ 0 . We also have V 0 = A 0 0 , so and we obtain the following for the engineering stress, σ eng :
σ eng = P z A 0 = d d λ z ( Δ F d V 0 ) (14.3)
Here A 0 , 0 and V 0 are the cross sectional area, length and volume of the undeformed material.
The deformation free energy normalized by the undeformed volume is the deformation free energy density, Δ f d :
Δ f d = Δ F d V 0 (14.4)
With these definitions the following simple expression gives the engineering stress in terms of the deformation free energy density, for uniaxial deformation along the z axis:
σ eng = d d λ z ( Δ f d ) (14.5)

14.4 Simple Shear

In a simple shear deformation, two parallel surfaces within a material remain at a fixed separation as they slide past one another, as illustrated in Figure 91. In this example an applied shear force, P xy is applied to the two surfaces, which each have an area, A. The separation between the two surfaces is h , and the relative displacement of the two surfaces is u . For these quantities we obtain the shear sress, σ xy , the engineering shear strain, e xy or γ , and the shear modulus, G :
Note that e xy and γ are both commonly used symbols to represent the shear strain.
image: 130_home_ken_Mydocs_MSEcore_331_figures_Shear_Deformation.png
Figure 91: Diagrammatic representation of shear deformation.
The force is again given by the derivative of the free energy with respect to the displacement in the direction of the applied force:
P xy = d( Δ F d ) du = d( Δ F d ) d γ d γ du = 1 h d d γ ( Δ F d ) (14.6)
The shear stress is then obtained by dividing by the cross-sectional area (note: hA=V , ν el = n el /V, σ xy = P xy /A :
σ xy = d d γ ( Δ f d ) (14.7)
One complicating factor of simple shear in comparison to uniaxial extension or compression is that the strain ellipsoid is rotated continuously during the shear experiment. As a result we cannot connect the principal axes of the coordinate system in a straightforward way. The situation for shear in the x-y plane is illustrated in Figure 92. Shear is a two-dimensional stress state, so the material is undeformed in the third direction, in this case the z direction, so λ z =1. The principal strains in the x-y plane are λ 1 and λ 2 . While the relationship between the orientation of the strain ellipsoid is a complicated one, the shear strain itself is given by the difference between λ 1 and λ 2 :
γ = e xy = λ 1 - λ 2 (14.8)
Eq. 14.8 will enable us to determine the deformation free energy of a material undergoing a shear deformation, and is used below in Section15.4.
image: 131_home_ken_Mydocs_MSEcore_331_figures_Shear_Deformation_princiapl_extension_ratios.png
Figure 92: Extension ratios for simple shear.

14.5 Torsion

An important geometry for characterizing soft materials is the torsional geometry shown in Figure 93. In this material a cylindrical or disk-shaped material is twisted about an axis of symmetry. The material could be a long, thin fiber (Figure 93a) or a flat disk sandwiched between two plates (Figure 93b). We obtain the shear modulus by looking at the torsional stiffness of mateiral, i.e. , the Torque, T required to rotate the top and bottom of the fiber by an angle θ 0 .
image: 132_home_ken_Mydocs_MSEcore_331_figures_shear_in_torsion.png Figure 93: Fiber torsion.
We define a cylindrical system with a z axis along the fiber axis. The other axes in this coordinate system are the distance r from this axis of symmetry, and the angle θ around the z axis. The shear strain in the θ -z plane depends only on r , and is given by:
e θ z =r d θ dz =r θ 0 (14.9)
The corresponding shear stress is obtained by multiplying by the shear modulus, G :
σ θ z =Gr θ 0 / (14.10)
We integrate the shear stress to give the torque, T :
T= 0 d/2 r σ θ z 2 π r r = π G θ 0 d 4 32 (14.11)
This geometry is commonly used in an oscillatory mode, where θ is oscillated at a specified frequency. In this case the torque response is obtained by using the dynamic shear moduli defined in Section16.6 in place of G .

15 Rubber Elasticity

Crosslinked rubbers are unique in that thermodynamic arguments can be used to predict their elastic moduli with remarkable accuracy. Our starting point will be a description of the free energy of an elastomer as a function of its deformation. This free energy is dominated by entropic contributions arising from restrictions on the number of different conformations (or shapes) that polymer strands are able to adopt. The detailed descriptions of polymer chain statistics given earlier were developed so that we would be in a position to describe the mechanical properties of rubbery materials. The description of rubber elasticity given here is proof that statistics and thermodynamics are actually useful!!

15.1 Molecular Deformation

Different types of molecular deformations, and the characteristic lengths and forces, are illustrated in Figure 94. In rubber elasticity we are interested in deformations due to changes in the overall shapes of the molecules, with typical molecular forces in the range of 10 pN and deformations in the range of 10 nm. In this force regime we don't need to worry about the specific chemistry of the polymer backbone. All that matters is the distance between the two ends of the molecule. Figure 95 shows the situation in this regime of 'entropic elasticity' in more detail.
image: 133_home_ken_Mydocs_MSEcore_331_figures_Molecular_Deformation_1.png
Figure 94: Molecular deformation occurring at different force and length scales (from From Osterhelt et al.: http://www.iop.org/EJ/abstract/1367-2630/1/1/006/).
image: 134_home_ken_Mydocs_MSEcore_331_figures_single_molecule_deformation.svg
Figure 95: Illustration of force, P d , that must be applied to a single molecule to maintain and end-to-end separation of R .
The origins of entropic elasticity can be understood in terms of the probability density, P xyz , that describes the relative probability that the end-to-end vector has a specific set of x, y and z components. This probability depends on the magnitude of the end-to-end vector normalized by R 0 , the root-mean-squared value of the end-to-end vector . For our purposes here we'll see that we don't actually need the normalization for P xyz - we just need to understand how it depends on | R | . For simplicity we refer to | R | simply as R in our discussion, with P xyz given as follows:
P xyz exp ( -3 R 2 2 R 0 2 ) (15.1)
The value of P xyz is proportional to the number of molecular configurations, Ω , that are consistent with this particular value of R , so we can write:
Ω ( R ) =C exp ( -3 R 2 2 R 0 2 ) (15.2)
where C is some constant that we don't know.
The entropy is obtained from the number of possible molecular shapes:
S d ( R ) = k B ln Ω ( R ) = k B ln C- 3 k B R 2 2N a 2 (15.3)
The free energy is assumed to be entirely dominated by this entropic contribution:
F d ( R ) =-T S d ( R ) =- k B T ln C+ 3 k B T R 2 2N a 2 (15.4)
In reality, there is a small enthalpic contribution to the stretching free energy as well, because stretching the chain increases the relative proportion of higher energy gauche bonds, but we're not going to worry about that correction here. Because we are only every interested in changes in the free energy, constant terms that don't depend on R are going to drop out when we compare free energies, so the detailed value of this constant is not going to matter.
This expression for the elastic free energy of an individual molecule in terms of its end-to-end distance is remarkably simple. It is one of the most important and widely used results in polymer science. In the following section, we will see how it is used to obtain estimates for the elastic modulus of a crosslinked elastomer.

15.2 Free energy of a stretched rubber

When a rubber is stretched, the free energy increases. In general the free energy increase has enthalpic and entropic contributions ( F d = H d -T S d ). The basic assumption of rubber elasticity theory is that the free energy increase due to deformation is dominated by the decrease in the entropy. In other words | T Δ S d | >>| Δ H d | . The second assumption we will make is that the deformation at a microscopic level mimics the deformation at a macroscopic level. In other words, relative changes in the spacings between crosslink points are identical to relative changes in the overall sample demonstrations. This assumption is referred to as the affine deformation assumption.

15.3 Free energy change due to deformation

image: 135_home_ken_Mydocs_MSEcore_331_figures_Free_Energy_Change_from_Deformation.png
Figure 96: Contrasting free energies before and after deformation
What really matters when determining the mechanical properties of a material is the change in free energy resulting from a deformation, which we express in terms of the molecular extension ratios in the x, y and z directions. Prior to deformation, we have R = R x x ˆ + R y y ˆ + R z z ˆ , so that R 2 = R x 2 + R y 2 + R z 2 . The undeformed free energy, F 0 , is given by using this expression for R 2 in Eq. 15.4:
F 0 =- k B T ln C+ 3 k B T 2N a 2 { R x 2 + R y 2 + R z 2 } (15.5)
After deformation, the values of R x , R y and R z are each multiplied by the respective extension ratio ( R = λ x R x x ˆ + λ y R y y ˆ + λ z R z z ˆ ), so that R 2 = λ x 2 R x 2 + λ y 2 R y 2 + λ z 2 R z 2 . The deformed free energy, F def is:
F def =- k B T ln C+ 3 k B T 2N a 2 { λ x 2 R x 2 + λ y 2 R y 2 + λ z 2 R z 2 } (15.6)
The free energy change due to deformation of the molecule is given as follows:
Δ F d = F def - F 0 = 3 k B T{ ( λ x 2 -1 ) R x 2 +( λ y 2 -1 ) R y 2 +( λ z 2 -1 ) R z 2 } 2N a 2 (15.7)
This result is for the deformation of a single polymer molecule, which for a crosslinked elastomer corresponds to a segment that connects crosslink points. There are a huge number of these segments in a macroscopic chunk of rubber. To get the free energy change for the material as a whole, we need to replace R x 2 by n el R x 2 , where n el is the total number of network strands and R x 2 is the average value of R x 2 for these segments. We need to make similar substitutions for R y 2 and R z 2 to obtain the following:
Δ F d = 3 k B T n el { ( λ x 2 -1 ) R x 2 +( λ y 2 -1 ) R y 2 +( λ z 2 -1 ) R z 2 } 2N a 2 (15.8)
Now we assume that the material was isotropic when it was crosslinked, so that R x 2 = R y 2 = R z 2 = R 2 /3 . With this assumption, and with R 0 2 -N a 2 , we obtain the following result:
Δ F d = k B T n el R 2 { λ x 2 + λ y 2 + λ z 2 -3 } 2 R 0 2 (15.9)
Note that R 2 is the mean square end-to-end distance of the polymer strands that span the crosslinks, and that R 0 2 is the value of R 2 when the polymer strands obey random walk statistics. The relationship between these two quantities depends on the conditions of the crosslinking reaction.
It is often useful to work in terms of intensive free energy changes (free energy per unit volume). The free energy of deformation per unit volume Δ f d is obtained very simply from Δ F d by dividing by V , the volume of a sample. We retain the same expression as shown above, but with the strand concentration, ν el substituted for the number of strands, n el :
Δ f d = Δ F d V = k B T v s β { λ x 2 + λ y 2 + λ z 2 -3 } 2 (15.10)
where ν el and β are defined as follows:
ν el = n el V ; β = R 2 R 0 2 (15.11)
Like the intensive free energy, the strand concentration, v s is useful because it does not depend on the overall size and shape of the elastomer. The assumptions used to develop the free energy expression given above correspond to the simplest model of rubber elasticity, which is often referred to as the Neohookean model . In our discussion of the specific cases of shear and uniaxial deformation that we discuss below, we make the additional assumption that the material is incompressible, in which case the product of the three principal extension ratios is one:
λ 1 λ 2 λ 3 =1 (15.12)

15.4 Shear deformation of an elastomer

For pure shear we need to relate the shear strain, γ , to the extension ratios, λ 1 , λ 2 and λ 3 . Pure shear is a two dimensional stress state, with λ 3 = 1 so the Neohookean strain energy function (Eq. 15.10) reduces to the following:
Δ f d = k B T v s β { λ 1 2 + λ 2 2 -2 } 2 = k B T v s β ( λ 1 - λ 2 ) 2 2 = k B T v s β γ 2 2 (15.13)
The shear stress is obtained by differentiation of the strain energy function with respect to γ (Eq. 14.7):
σ xy = d d γ ( Δ f d ) =G γ (15.14)
where G is the shear modulus, given as follows:
G= k B T v s β (15.15)
Note that the shear stress is proportional to the shear strain, even for large values of the strain.

15.5 Uniaxial deformation of an elastomer

For uniaxial extension or compression the deformation is applied along one axis, which we define as the z axis. We assume that our material is isotropic, so that the extensions in the x and z directions are identical to one another, i.e. λ x = λ y . In addition, the material is assumed to be incompressible, so λ x λ y λ z = 1. We therefore have λ x = λ y = λ z -1/2 . We can therefore write the free energy of a deformed elastomer as a function of the single extension ratio, λ z :
Δ f d = G 2 { λ z 2 + 2 λ z -3 } (15.16)
The engineering stress is obtained by differentiating with respect to λ z (Eq. 14.5):
σ eng = d d λ ( Δ f d ) =G{ λ z - 1 λ z 2 } (15.17)
Young's modulus (E) is defined as the derivative of the stress with respect to strain, evaluated at low strain ( λ z =1 ):
E σ eng λ z | λ z =1 =G . { 1+ 2 λ z 3 } | λ z =1 =3G=3 v s β k B T, (15.18)
The concentration of network strands ( v s ) is the inverse of the volume per strand. This can be calculated from the molecular weight and density. Defining M s as the number average molecular weight of a network strand (molecular weight between crosslinks) gives:
v s ( strands volume ) = ρ ( mass/volume ) M s ( mass/mole ) N av ( strands/mole ) (15.19)
where N av is Avogadro's number ( 6.02x1 0 23 ). With R= k B N av =8.314 J/moleK the expression
E= 3 ρ β RT M s
Finally, note that 1 J/m 3 = 1 Pa. Stresses and elastic moduli have units of force/area or energy/volume.
Exercise:
Calculate the expected value of Young's modulus for a crosslinked polyisoprene that has a number average molecular between crosslinks of 4000 g/mole. The density of polyisoprene is 0.9 g/cm 3 . Assume that the polymer was crosslinked under equilibrium conditions in the melt state.
Solution:
Because the polymer was crosslinked under equilibrium conditions, the network strands obey random walk statistics, with = 1. To keep everything in SI units, we need the density in units of kg/m 3 (1 g/ c m 3 = 1000 kg/ m 3 ), and M s in kg/mol. In our example, ρ = 900 kg/ m 3 and M s = 4 kg/mol. We'll also assume that we are interested in the elastic modulus near room temperature (T 300 K).
E= 3( 900kg/ m 3 ) ( 8.314 J molK ) ( 300K ) 4kg/mol =1.6x1 0 6 J/ m 3 =1.6x1 0 6 Pa .

Linear Viscoelasticity

16.1 Intro to Time-Dependent Behavior

Many of the concepts of rubber elasticity still apply under certain situations to materials that are not actually crosslinked. Consider, for example, the behavior of Silly-Putty. Silly-putty is based on silicones that are uncrosslinked. Nevertheless, this material still bounces like an ordinary crosslinked rubber ball. This is because over very short times it behaves elastically. The deformation energy stored in the material as it comes into contact with the surface and deforms is available to propel the material back into the air as the deformation relaxes and this strain energy is converted to kinetic energy. If we let the sample sit on a surface for a long period of time, however, it eventually flows and behaves as a liquid.
(a)
image: 136_home_ken_Mydocs_MSEcore_331_figures_Short_Term_Viscoelastic1.png
(b)
image: 137_home_ken_Mydocs_MSEcore_331_figures_Short_Term_Viscoelastic2.png
(c)
image: 138_home_ken_Mydocs_MSEcore_331_figures_Short_Term_Viscoelastic3.png
(d)
image: 139_home_ken_Mydocs_MSEcore_331_figures_Long_Term_Viscoelastic1.png
(e)
image: 140_home_ken_Mydocs_MSEcore_331_figures_Long_Term_Viscoelastic2.png
(f)
image: 141_home_ken_Mydocs_MSEcore_331_figures_Long_Term_Viscoelastic3.png
Figure 97: Illustration of viscoelastic behavior. Images of a viscoelastic material behaving elastically as it bounces from a surface (a-c) and the same material spreading over a surface at long times (d-f).
Before we discuss the way that the properties of viscoelastic materials are characterized, it is useful to discuss the limiting behaviors of purely elastic solids and perfectly viscous liquids. We illustrate this by showing how different materials respond to a shear strain, γ that begins to increase at some constant rate to some final value, and is then fixed at this final value as illustrated in the top panel of Figure 99. The response for solids, liquids and viscoelastic materials is as described below.
Solids:
In this case the stress is proportional to the strain, and the rate at which the strain is applied does not matter. As a result the time dependence of the stress looks just like the time dependence of the strain. The slopes of the time-dependent stress and strain curves (the top two curves in Figure 99) are related to one another by the shear modulus of the material.
Liquids:
In liquids the stress is proportional to the rate at which the strain is applied, and is independent of the current strain. Liquids do not store any strain energy, and as soon as the strain stops changing, the stress drops back to zero. The stress for the time dependent strain shown in Figure 99 is constant while the stress is increasing and is zero otherwise. The stress is given by the shear viscosity , η , as follows:
σ xy = η d γ dt = η γ ˙ (16.1)
As a simple illustration of the significance of the viscosity, we can calculate the time it takes for a small sphere of metal to drop to the bottom of pool of water. Suppose we use an iron particle with a radius of 1 μ m. The situation is as shown in Figure 98. The gravitational force, F, causing the ball to sink is proportional to the volume of the sphere, and the difference in densities between the solid and liquid:
image: 142_home_ken_Mydocs_MSEcore_332_figures_falling_ball.svg
Figure 98: Terminal velocity of a sphere descending in a viscous fluid.
F=g 4 3 π R 3 ( ρ s - ρ ) (16.2)
Here ρ s is the density of the solid sphere, ρ is the density of the surrounding liquid, R is the radius of the sphere and g is the gravitational acceleration (9.8 m/s 2 ). This force is balanced by the viscous force exerted by the water on the sphere as the sphere moves the liquid with a velocity, v . For a viscous liquid this force is proportional to the velocity. It is also proportional to the viscosity of the liquid and to the radius of the sphere. The specific relationship is as follows:
F=6 π V η R (16.3)
At steady state, the velocity reaches a value that is obtained by equating the forces in Eqs.16.2 and16.3. In this way we obtain:
V= 2 9 g R 2 η ( ρ s - ρ ) (16.4)
Iron has a density of 7.87 g/cm 3 (7870 kg/m 3 ) and water has a density of 1000 kg/m 3 and a viscosity of 10 -3 Pa-sec. From this we get a velocity of 16 μ m/s.
Viscoelastic Materials:
Viscoelastic materials have characteristics of both solids and liquids. The shear stress depends on both the shear stress and the strain rate, and is not a simple function of either one. The shear stress actually depends on the details of the previous strain history, as described in more detail below.
image: 143_home_ken_Mydocs_MSEcore_331_figures_generic_viscoelasticity.png
Figure 99: Interactive contrast of elastic, viscous, and viscoelastic response.
While the shear geometry is most commonly used to characterize liquid-like materials, solid-like materials are more commonly investigated in a uniaxial tensile geometry, as illustrated in Figure 100. Viscoelastic solids will have a value for Young's moduls, E , that depends on the time scale of the measurement. Imagine a tensile experiment where a strain of e 0 is instantaneously applied to the sample. In a viscoelastic material, the resulting stress will decay with time while we maintain the strain at this fixed value. The time dependence of the resulting tensile stress, σ , enables us to define a time-dependent relaxation modulus , E( t ) :
E( t ) = σ ( t ) e 0 (16.5)
image: 144_home_ken_Mydocs_MSEcore_332_figures_uniaxial_viscoelastic.svg
Figure 100: Experimental geometry and time dependent strain for for the determination of the time-dependent relaxation modulus and the frequency-dependent dynamic moduli.
We are often interested in the application of an oscillatory strain to a material. Examples include the propagation of sound waves, where wave propagation is determined by the response of the material at the relevant frequency of the acoustic wave that is propagating through the material. In an oscillatory experiment, referred to as a dynamic mechanical experiment in Figure 100, the applied shear strain is an oscillatory function with an angular frequency, ω , and an amplitude, e 0 :
e( t ) = e 0 sin ( ω t ) (16.6)
Note that the strain rate, de dt , is also an oscillatory function, with the same angular frequency, but shifted with respect to the strain by 90 :
de dt = e 0 ω cos ( ω t ) = γ 0 ω sin ( ω t+ π /2 ) (16.7)
The resulting stress is also an oscillatory function with an angular frequency of ω , and is described by its amplitude and by the phase shift of the relative to the applied strain:
σ ( t ) = σ 0 sin ( ω t+ φ ) (16.8)
Now we can define a complex modulus with real and imaginary components as follows:
E * = E +i E =| E * | e i φ (16.9)
There are a couple different ways to think about the complex modulus, E * . As a complex number we can express it either in terms of its real and imaginary components ( E and E , respectively), or in terms of its magnitude, | E * | and phase, φ . The magnitude of the complex modulus is simply the stress amplitude normalized by the strain amplitude:
| E * | = σ 0 / e 0 (16.10)
The phase angle, φ , describes the lag between the stress and strain in the sample. Examples for φ =9 0 (the maximum value, characteristic of a liquid) and φ =4 5 are shown in Figure 101.
image: 145_home_ken_Mydocs_MSEcore_332_figures_sinewave.png
Figure 101: Time dependent stress and strain in a dynamic mechanical experiment.
In order to understand the significance of the real and imaginary components of E * we begin with the Euler formula for the exponential of an imaginary number:
e i φ = cos φ +i sin φ (16.11)
use this expression in Eq. 16.9, we see that the storage modulus , E gives the stress that is in phase with the strain (the solid-like part), and is given by the following expression:
E =| E * | cos φ (16.12)
Similarly, the loss modulus , E , gives the response of the material that is in phase with the strain rate (the liquid-like part):
E =| E * | sin φ (16.13)
We can combine Eqs. 16.12 and 16.13 to get the following expression for tan φ , commonly referred to simply as the loss tangent:
tan φ = E E (16.14)
The loss tangent gives the ratio of the energy dissipated in one cycle of an oscillation to the maximum stored elastic energy during this cycle.

16.2 Shear Relaxation Modulus

In a shear geometry (Fig. 102) we can also conduct relaxation experiments or oscillatory experiments. The shear relaxation modulus is defined as the ratio of the resulting time dependent stress to this initial strain:
G( t ) σ xy ( t ) γ 0 (16.15) As with E( t ) , the shear modulus is independent of γ 0 , provided that γ 0 is sufficiently small. This low strain regime defines the regime of linear viscoelasticity. In the remainder of this section we confine ourselves to this linear regime, assuming a shear geometry so that the relevant time and frequency-dependent property is the shear modulus, G .
image: 146_home_ken_Mydocs_MSEcore_331_figures_shear_strain.png
Figure 102: Application of an instantaneous shear strain.

16.3 Boltzmann Superposition Principle

In the linear viscoelastic regime the effects of strains applied at different times are additive, a concept known as the
Boltzmann superposition principle
. A simple example is illustrated in Figure 103, where the following step strains are applied to the sample:
image: 147_home_ken_Mydocs_MSEcore_331_figures_boltzmann_superposition.png
Figure 103: Material response to a series of step strains.
The stress contribution from these different strains is obtained by multiplying by the shear relaxation modulus, evaluated at the time that has passed since the strain was applied. If strain was applied at t 1 , for example, the stress at some later time, t , is obtained by multiplying by G( t- t 1 ) . Adding up the contributions from the three different step strains in our example leads to the following:
σ xy = Δ γ 1 G( t- t 1 ) + Δ γ 2 G( t- t 2 ) + Δ γ 3 G( t- t 3 ) (16.16)
This expression can easily be extended to include an arbitrary number of step strains:
σ xy =i Δ γ i G( t- t i ) (16.17)
By taking a sum of very small strain increments, we can generalize this expression to account for a continuously changing strain. We begin by writing Δ γ i in the following way:
Δ γ i = d γ dt ( t i ) Δ t i = γ ˙ ( t i ) Δ t i (16.18)
This substitution leads to the following expression for the stress:
σ xy =i γ ˙ ( t i ) G( t- t i ) Δ t i (16.19)
We can write this in integral form by taking Δ t i 0 and replacing the summation by an integral over all times less than the current time:
σ xy ( t ) = - t γ ˙ ( t i ) G( t- t i ) d t i (16.20)
This expression gives the current stress that remains as a result of all of the strains introduced at different times in the past. It is often convenient to change variables so that s=t- t i , d t i =-ds . Note that s = when t i =- , and s=0 when t i =t . The integral can therefore be rewritten as follows:
σ xy ( t ) =0 γ ˙ ( t-s ) G( s ) ds (16.21)
Note that because s represents time, G( s ) and G( t ) represent the same time-dependent relaxation modulus. We are just using different variables ( s and t ) to represent time.
The simplest application of Eq. is 'steady shear', where a constant strain rate is applied. Because γ ˙ is independent of time in this case, it can be brought outside the integral, giving the following:
σ xy = γ ˙ 0G( t ) dt (16.22)
The following expression is obtained for the viscosity:
η 0 = σ xy γ ˙ =0G( t ) dt (16.23)
We refer to this viscosity as the
zero shear viscosity
, because it assumes that the material is in the linear viscoelastic regime, where the applied stress is low enough so that the stress is proportional to the magnitude of the applied strain. The subscript '0' is a reminder to us that the zero shear viscosity corresponds to a limiting value of the viscosity that is obtained at a very low strain rate.

16.4 Idealized Relaxation Curves

One thing to keep in mind when looking at the viscoelastic properties of materials is that processes occur over a very large range of time scales. To capture all of these time scales we typically plot the time on a logarithmic scale. The same is true for the frequency-domain experiments discussed below. The relaxation modulus can often vary of several orders of magnitude, so that we also plot the relaxation modulus itself on a logarithmic scale.
Three different idealized forms of the relaxation modulus are shown in Figure 104. The simplest behavior corresponds to a
Maxwell model
, where a single relaxation time, τ g , describes the decay in the modulus from the glassy value of G g :
G( t ) = G g exp ( -t/ τ g ) (16.24) This relaxation time could, for example, describe the decay in the modulus from a glassy value of 1 0 9 Pa. This relaxation time can be viewed as a 'glass transition time', and its existence is a clue that some equivalence between time and temperature must exist. The relaxation time characterizing the glass transition depends very strongly on the temperature, but is independent of the polymer molecular weight. If the polymer molecular weight is very high, however, the polymers become entangled with one another, and behave elastically for times that are too short for these entanglements (shown schematically in Figure 105) to relax by molecular diffusion. In addition to τ g , there is a second transition time, τ e , determined by the lifetime of these molecular entanglements. Addition of this second relaxation results in the following expression for the relaxation modulus:
G( t ) = G g exp ( -t/ τ g ) + G e exp ( -t/ τ e ) (16.25)
Here G e is the plateau modulus (often referred to as G N 0 in the literature). Entanglements behave like crosslinks, but they have finite lifetimes as the molecules diffuse and the entanglements are eliminated and reformed elsewhere. The plateau modulus is given by the concentration of entanglements, v e , or equivalently, by the number average molecular weight between entanglements, M e (referred to as the entanglement molecular weight):
G e = v e k B T= ρ RT M e (16.26)
image: 148_home_ken_Mydocs_MSEcore_331_figures_viscoelasticity_relaxation_curves.png
Figure 104: Characteristic relaxation curves.
image: 149_home_ken_Mydocs_MSEcore_331_figures_Molecular_Entanglement.png
Figure 105: Schematic representation of a molecular entanglement.
The actual relaxation behavior of polymeric materials is complex, and can generally not be described in detail by including only one or two relaxation times. The actual relaxation behavior can always be described with sufficient accuracy by a
generalized Maxwell model
, where we include a large number of individual exponential relaxation processes:
G( t ) = i=1 N r G i exp ( -t/ τ i ) (16.27)
Inclusion of a sufficiently large number of relaxation processes (large N r ) enables very complicated relaxation behavior to be modeled accurately.

16.5 Temperature Dependence

A characteristic of single-component polymeric materials is that the temperature dependence of all of the different relaxation times is the same. The concept is illustrated in Figure 106, which shows the relaxation modulus at two different temperatures for an idealized material with two different relaxation times. The relaxation times scale with temperature according to the same temperature dependent shift factor, a T . The shift factor is equal to one at a reference temperature, T ref , and increases as the temperature is decreased. Because both relaxation times are multiplied by the same factor, the curve is shifted linearly along the logarithmic x axis.
image: 150_home_ken_Mydocs_MSEcore_331_figures_Temp_Shift.png
Figure 106: Temperature dependence of a material with two relaxation times.
Because all of the relaxation times are multiplied by a T when the temperature is changed, data obtained at different temperatures superpose when plotted as a function of t/ a T . As a result the relaxation behavior over a large range of times can be obtained by measuring the relaxation spectrum at different times, and shifting the data to the reference temperature.
The
Vogel equation
is an empirical equation that is often used to describe the temperature dependence of the shift factor, a T :
log ( a T ) =A+ B T- T (16.28)
The factor A is determined from the requirement that a T =1 at T= T ref :
A=- B T ref - T (16.29)
The relaxation times diverge to at T= T , which in free volume theory is the temperature at which the free volume of the equilibrium liquid goes to zero.

16.6 Dynamic Moduli

Experimentally, a wide variation in time scales is accessed by oscillating or vibrating the sample and measuring the frequency response of the material. Consider, for example, an oscillatory shear strain:
γ = γ 0 sin ( ω t ) (16.30)
The strain rate is also sinusoidal:
γ ˙ ( t ) = d γ ( t ) dt = ω γ 0 cos ( ω t ) = ω γ 0 sin ( ω t+ π /2 ) (16.31)
Note that the strain and the strain rate are out of phase by π /2 (90 ). This concept of a phase difference is very important in understanding the frequency dependent dynamic moduli. In this case we use Boltzmann superposition to obtain an expression for the stress:
σ xy ( t ) =0 γ ˙ ( t-s ) G( s ) ds=- ω γ 0 0 cos { ω ( t-s ) } G( s ) ds (16.32)
Now we make use of the following trigonometric identity:
cos ( a-b ) = sin ( a ) sin ( b ) + cos ( a ) cos ( b ) (16.33)
We can therefore write the time dependent stress in the following way.
σ xy ( t ) = γ 0 ω [ 0G( s ) sin ( ω s ) ds ] sin ( ω t ) + γ 0 ω [ 0G( s ) cos ( ω s ) ds ] cos ( ω t ) (16.34)
We define G (storage modulus) and G (loss modulus) such that G describes the component of the stress that is in phase with the strain (the elastic component) and G describes the component of the stress that is in phase with the strain rate (the viscous component).
σ xy ( t ) γ 0 = G ( ω ) sin ( ω t ) + G ( ω ) cos ( ω t ) (16.35)
By comparing Equations 16.34 and 16.35, we obtain the following for G and G :
G ( ω ) = ω 0G( t ) sin ( ω t ) dt G ( ω ) = ω 0G( t ) cos ( ω t ) dt (16.36)
In is useful to consider the behavior of G and G in limiting cases where the system is a perfectly elastic solid with no viscous character, and where the material is a Newtonian liquid with no elastic character:
Perfectly Elastic System:
In this case the shear modulus is independent of both time and frequency. For γ ( t ) = γ 0 sin ( ω t ) we have:
σ xy ( t ) =G γ ( t ) =G γ 0 sin ( ω t ) (16.37)
Comparing to Eq. 16.35 gives:
G =G G =0 (16.38)
Perfectly Viscous System:
We again have γ ( t ) = γ 0 sin ( ω t ) , but this time the shear stress depends only on the strain rate:
σ xy ( t ) = η d γ ( t ) dt = η ω γ 0 cos ( ω t ) (16.39)
Comparing to Eq. 16.35 gives:
G =0 G ( ω ) = ω η (16.40)

16.6.1 Phase Angle and Loss Tangent

We define different components of the dynamic shear modulus in a way that is exactly analogous to the dynamic tensile moduli discussed earlier. The quantities G and G can be viewed as the real and imaginary components of a a complex modulus, G * :
G * ( ω ) = G ( ω ) +i G ( ω ) (16.41)
As with E * we can represent G * by its magnitude, which we refer to as | G * | and its phase angle, which we refer to as φ :
| G * ( ω ) | = { ( G ( ω ) ) 2 + ( G ( ω ) ) 2 } 1/2 φ = tan -1 ( G ( ω ) G ( ω ) ) (16.42)
The significance of | G * | and δ are illustrated by writing the sinusoidal displacement and the resulting stress in the following way:
γ ( t ) γ 0 = sin ( ω t ) (16.43)
σ xy ( t ) γ 0 =| G * | sin ( ω t+ φ ) (16.44)
| G * | gives the ratio of the stress amplitude to the strain amplitude, and φ is the phase lag between the stress and strain. The following expressions relate | G * | , G , G and φ :
G =| G * | cos ( φ ) G =| G * | sin ( φ ) tan ( φ ) = G G (16.45)
The
loss tangent
, tan( φ ), gives the ratio of dissipated and stored elastic energy during a given cycle of deformation. This quantity is maximized when the inverse of the frequency corresponds to a characteristic relaxation time of the material. For a material with two characteristic relaxation times, τ 1 and τ 2 , tan( δ ) will be maximized at ω =2 π / τ 1 and ω =2 π / τ 2 .

16.7 Torsional Resonator

We can now consider the specific example of a torsional resonator, shown below in Figure 107. The equation of motion for the spring as it is being twisted is[8]:
T= K t θ +I d 2 θ d t 2 (16.46)
where θ is the rotation angle, T is the torque on the resonator, K t is the torsional stiffness and I is the moment of inertia of the system, which depends on the details of the inertial bar used in the experiment. For a cylindrical fiber K t is obtained from Eq. 14.11:
K t T θ = π G d 4 32 (16.47)
image: 151_home_ken_Mydocs_MSEcore_332_figures_torsional_resonator_no_numbers.svg
Figure 107: Torsional resonator.
The resonant angular frequency of the oscillator, ω n is given by the following expression:
ω n = K t I (16.48)
The solution to Eq. 16.46 for the case where T=0 and we just let the fiber move (by twisting to a certain angle and letting it go, for example) is as follows:
θ = θ 0 cos ( ω n t ) (16.49)
This analysis assumes that the system is entirely elastic. What if the spring has some viscoelastic character to it? We begin by using Euler's formula (Eq. eq:euler formula ) to rewrite Eq. 16.49 in the following way:
θ = θ 0 Real ( exp ( i ω n t ) ) (16.50)
Now all we need to do is to replace the modulus, G, with the complex modulus, G * , and everything works out fine. Here's what happens:
To understand this last point, suppose write ω n * in the following way:
ω n * = ω n +i ω n (16.51)
Putting this back into Eq. 16.50 gives:
θ = θ 0 Real ( exp ( i ω n t ) ) exp ( - ω n t ) (16.52)
Using the Euler formula (Eq. eq:euler formula) to expand exp ( i ω n t ) and taking the real part gives:
θ = θ 0 cos ( ω n t ) exp ( - ω n t ) (16.53)
So the imaginary part of the complex frequency describes the decay of oscillation with time.

16.8 Viscoelastic Models

Models of viscoelasticity can be represented visually by connecting solid-like and liquid-like elements together. In materials science we're used to the solid-like elements. They are simply elastic springs where the force is proportional to the extension of the spring. In terms of stress and strain, the stress, σ is proportional to the strain, e . There is not time-dependence to the behavior of an ideal spring. The stress is proportional to the strain, regardless of how fast or slow the strain is applied. In many real materials the rate at which the strain is applied matters as well, and this is where the dashpots come in. A dashpot is a liquid-like element where the stress is proportional to the rate at which the strain is changing, with the viscosity being the ratio between the stress and the strain rate. Schematic representations of springs and dashpots are shown in Figure 108. We assume a tensile geometry in our discussion of the viscoelastic models, so the relevant modulus is E as opposed to G , although this same formalism is generally applied to the shear moduli as well.
image: 152_home_ken_Mydocs_MSEcore_332_figures_springs_and_dashpots.svg
Figure 108: Schematic representations of springs and dashpots used to represent the time-dependent response of viscoelastic materials.

16.8.1 Maxwell model

The simplest viscoelastic model that contains both liquid-like and solid-like elements is the Maxwell model consisting of a linear dashpot with viscosity η in series with a linear spring with a modulus, E (see Figure 109). Because the elements are in series the strain is the same in each one of them. This stress can be related to the strain e 1 in the dashpot and the strain e 2 in the spring through the following expressions:
σ = η d e 1 dt (16.54)
σ =E e 2 (16.55)
image: 153_home_ken_Mydocs_MSEcore_332_figures_maxwell_model.svg
Figure 109: Maxwell model.
The total strain, e is e 1 + e 2 , so the time derivative of the total strain is:
de dt = d e 1 dt + d e 2 dt = σ η + 1 E d σ dt (16.56)
In a stress relaxation experiment the strain instantaneous increased to an initial value of e 0 and we follow the relaxation of the stress at this fixed strain. Because the dashpot cannot respond instantaneously, all of the initial strain is in the spring, i.e., e 0 = e 2 and the initial stress, σ 0 is equal to E e 0 . The solution to Eq. 16.56 for this initial condition and for de/dt=0 is:
σ ( t ) =E e 0 exp ( -t/ τ ) (16.57)
with τ = η /E . We divide by e 0 to obtain the relaxation modulus: E( t ) σ ( t ) / e 0 =E exp ( -t/ τ ) (16.58)
We can also obtain the solution for the stress in the case where we apply an oscillatory strain. The math is a bit more complicated, but we end up with the following complex modulus:
E * ( ω ) = Ei ω τ 1+i ω τ (16.59)
We can obtain the real and imaginary parts of E * by multiplying the top and bottom of Eq. 16.59 by 1+i/( ω τ ) , remembering that i 2 =-1 . Taking the real and imaginary components of E * gives E and E :
E ( ω ) = E ( ω τ ) 2 1+ ( ω τ ) 2 (16.60)
E ( ω ) = E ω τ 1+ ( ω τ ) 2 (16.61)

16.8.2 Standard Linear Solid

For a single maxwell element the relaxation modulus decays to zero at long times ( t τ ). In many real systems we are interested in describing the day of the relaxation modulus from a large value at short times to a much smaller value at larger times. The standard linear solid accounts for this by putting an elastic element with a modulus of E r (the 'relaxed' modulus). In parallel with a single Maxwell element as shown in Figure 110. Extra spring adds directly to stress response, so all we have to do is add E r to the expressions for E( t ) and E * ( w ) that were obtained from the Maxwell model:
E( t ) = E r + E 1 exp ( -t/ τ 1 ) (16.62)
E * ( ω ) = E r + E 1 i ω τ 1 1+i ω τ 1 (16.63)
image: 154_home_ken_Mydocs_MSEcore_332_figures_standard_linear_solid.svg Figure 110: Standard linear solid.
So we see that the standard linear solid describes the relaxation of the modulus from an initial value of E 1 + E r at very short times to a relaxed value of E r at long times. The standard linear solid is the simplest model for describing the transition of an amorphous, crosslinked polymer from glassy polymer behavior to rubbery behavior with typical values of E 1 and E r being 1 0 9 Pa and 1 0 6 Pa, respectively.

16.8.3 Generalized Maxwell Model

For real systems, the behaviors of E( t ) and E * ( ω ) are usually much more complicated than given by the predictions of the Maxwell model or standard linear solid. To describe the behavior for real materials, we can add an arbitrary number of Maxwell elements in parallel, resulting in the generalized Maxwell model shown in Figure 111. The stresses for each of the parallel elements are additive, so E( t ) and E * ( w ) are given by the following expressions: E( t ) = E r +j E j exp ( -t/ τ j ) (16.64)
E * ( ω ) = E r +j E j i ω τ j 1+i ω τ j (16.65)
with τ j = η j / E j .
image: 155_home_ken_Mydocs_MSEcore_332_figures_generalized_maxwell.svg Figure 111: Generalized Maxwell model.

16.8.4 Kelvin-Voigt Model

The Kelvin-Voigt model consists of a spring and dashpot in parallel, as shown in Figure 112. In this case the stresses in two elements are additive and the strains in the two elements are the same: σ = η de dt +Ee (16.66)
The Kelvin-Voigt model is the simplest model for the description of a creep experiment, where the stress jumps instantaneously from 0 to σ = σ 0 and we track the strain as a function of time. We have:
e( t ) = σ 0 E ( 1- exp ( t/ τ ) ) ; τ = η /E (16.67) The dynamic modulus for the Kelvin-Voigt element is given by the following expression:
E * =E+i ω η (16.68)
So that E is simply equal to E, and E is equal to ω η .
image: 156_home_ken_Mydocs_MSEcore_332_figures_Kelvin_Voight_model.svg Figure 112: Kelvin-Voigt Model.

16.8.5 Time-Temperature Superposition of Dynamic Mechanical Data

image: 157_home_ken_Mydocs_MSEcore_331_figures_Dynamic_Mechanical_Data.png
Figure 113: An interactive showing dynamic mechanical data
Measurement of G N 0 and M e
image: 158_home_ken_Mydocs_MSEcore_331_figures_Measurements_of_G_and_M.png
Figure 114: Plot showing temperature shifted curve
In this plot we show the temperature shifted master curve forG', from the PTBA data. The specific value of G N 0 is generally taken as the inflection point in G', or as the value of G' at the frequency at which δ is minimized. Once G' is known, the entanglement molecular weight, Me can be determined from the relevant expression.
Viscocity from Rheological Data
The zero shear viscosity, η 0 , at the temperature where a T = 1 is determined from the low frequency data, in the regime where G' is proportional to ω . The full temperature dependence of n is then given by the measured temperature dependence a dependence of a T .
image: 159_home_ken_Mydocs_MSEcore_331_figures_Rheological_Data.png
Figure 115: Rheological data

Semicrystalline Polymers

Many polymers of commercial importance are able to crystallize to a certain extent. The most widespread crystallizable polymers are polyethylene and isotactic polypropylene. As we will see, crystallization is never complete, and there is always a certain amorphous fraction in any polymer. For this reason, many of the concepts which apply to amorphous polymers are still valid for semicrystalline polymers. For example, semicrystalline polymers are often quite brittle at temperatures below the glass transition temperature of the amorphous fraction.
In the following pages, we introduce some common semicrystalline polymers, and discuss their basic structural feature polymers, in addition to kinetic issues that are important in their processing.

17.1 Structural Hierarchy in Semicrystalline Polymers

The structure of semicrystalline polymers is much more complex than the structure of amorphous polymers, and our discussion in the following pages is not as quantitative as our discussion of the random walk configurations of polymer chains in amorphous polymers. Instead, we illustrate the features of semicrystalline polymers by considering the following basic structural features:
  1. Helices formed by individual polymer molecules
  2. Perfect crystals formed by the lateral packing of these helices
  3. Lamellar crystallites formed by folded helices
  4. Semicrystalline regions formed by stacking of lamellar crystallites
  5. Spherulitic morphologies formed by radial growth of lamellar crystallites

17.2 The Structural Repeat Unit

The molecular structure is an important factor in the ability of a polymeric material to crystallize. One relevant parameter related to the structure vinyl polymers is the number of backbone carbons
structural repeat unit
, which is distinct from the
monomeric repeat unit
that is defined in terms of the structure of the monomer itself. Polyethylene, example, can be viewed as a repeating string of C H 2 units, so there one carbon atom per structural repeat unit, although there are two carbon atoms per monomeric repeat unit. Tacticity needs to be accounted for as well. Isotactic vinyl polymers have two backbone carbon atoms per structural repeat, because all non-hydrogen substituent groups are on the same side of the molecule. Syndiotactic vinyl polymers have four backbone carbon atoms per structural repeat, because the substituent groups alternate between different sides of the molecule.
Exercise:
What are the structural repeats for each of the following polymers:
  1. Poly(tetrafluoroethylene)
  2. Syndiotactic Polypropylene
  3. Isotactic Polypropylene
  4. Atactic Polystyrene
Solution:
  1. The repeating unit here is simply C F 2 since I can reconstruct the polymer simply by connecting these units to one another.
  2. The repeating unit here contains 4 backbone carbons and two methyl groups, since there is a methyl group on every other carbon, and these methyl groups alternate from one side of the backbone to the other.
  3. In this case there are two backbone carbons and a single methyl group (a single propylene unit).
  4. Since the polymer is atactic, it has an inherently disordered structure, with no structural repeat. This is why it is always amorphous.

17.3 Helix Formation

The first step in understanding the structure of crystalline polymer is to understand the shapes that the polymer molecules take when they pack within the crystalline structure. The simplest situation is for a molecule like polyethylene, for which the molecules in the crystalline structure are all fully extended, with all of their bonds in the trans configuration as shown in Figure 116.
image: 160_home_ken_Mydocs_MSEcore_331_figures_Planar_Zigzag.png image: 161_home_ken_Mydocs_MSEcore_331_figures_Helix_Problem_Diagram.png
Figure 116: Diagram of a planar zigzag (left), with the specific geometry for tetrahedrally bonded carbon (right). Note that the bond angles appear distorted in the leftmost portion of the figure because the horizontal and vertical scales are not necessarily equivalent.)
The planar zigzag configuration is one special case of the shapes of polymers along the chain axis. In the more general case, the backbone carbon atoms all lie along the a helix that is formed about some central axis. Helices of vinyl polymers form by rotation about the C-C single bonds which make up the backbones of the molecules. The curved line in the side view of the planar zigzag in Figure 116 represents the helix that passes through each of the carbon atoms in a single molecule in the polyethylene crystal structure.
The general notation describing the helical structure is x*y/z where x , y and z are defined as follows:
The planar zigzag can be formally viewed as a 1*2/1 helix, with 3 turns of the helix shown in Figure 116. From the top view, we see that all of the backbone atoms of this helix are in the same plane. More complicated helices form because of steric interactions between the other substituents that are placed on the carbon backbones. It's easy for polyethylene molecules to adopt an all-trans configuration because the hydrogen atoms are small and don't interfere with each other when the polymer chain adopts this particular shape. Fluorine atoms are larger than hydrogens, however, so poly(tetrafluorethylene) (Teflon TM ) adopts a more complicated, 13/7 helix, which is is shown below in Figure 117.
Exercise:
What is length of a polyethylene molecule with M = 10,000 g/mole, if the molecule is in all-trans (planar zigzag) configuration.
Solution:
For bond angle, θ , for singly bonded carbon atoms are 71 , and the C-C bond length is 1.54 Å. From the drawing below, we see that the projected length of this bond along the direction of the zigzag is equal to 1.25 Å. We just need to multiply this length by the number of carbons in the polyethylene chain, which we get by dividing the molecular weight of 10,000 g/mole by the molecular weight of a single C H 2 unit (14 g/mole). This gives the following for the total length, L :
L=1.25 Å ( 10,000 g/mole 14 g/mole ) =960 Å
image: 162_home_ken_Mydocs_MSEcore_331_figures_Helix_Formation.png
Figure 117: Representation of a 13/7 helix.

17.4 Crystalline unit cells

The structure of a perfect, three dimensional crystal, is defined by the lateral packing of polymer helices into repeating structure. As with any crystalline material, the structure is defined by the unit cell dimensions, and by the contents of the unit cell. The dimensions of the unit cell are defined by three lengths ( a , b and c ) and by three angles ( α , β and γ ), as illustrated in Figure 118.
image: 163_home_ken_Mydocs_MSEcore_331_figures_Crystalline_Unit_Cells.png
Figure 118: Diagram of a crystalline unit cell.
The perfect crystal structure of a polymer is defined by the helix type, unit cell dimensions and by the packing of the helices in the unit cell. Polyethylene can exist in multiple crystal structures. One of these, referred to as polyethylene I, is illustrated in Figure 119. This crystal structure has two 1*2/1 helices packed into an othorhombic unit cell ( α = β = γ = 90 ). As is often the case, the helix axis corresponds to the c direction of the unit cell.
image: 164_home_ken_Mydocs_MSEcore_331_figures_Polyethylene_Crystal_Structure.png
α = β = γ =9 0
a = 4.946 A
b = 7.418 A
c = 2.546 A
2 helices per unit cell
Figure 119: Structure of polyethylene I.
Polymer
Helix Type
a, b, c (Å)
α , β , γ
structural units/cell
polyethylene I
1*2/1
7.41, 4.95, 2.55
90, 90, 90
4
polyethylene II
1*2/1
8.09, 4.79, 2.53
90, 107.9, 90
4
poly(tetrafluoroethylene ) I
1*13/6
5.59, 5.59, 16.88
90, 90, 113.3
13
polypropylene (syndiotactic)
4*2/1
14.50, 5.60, 7.40
90, 90, 90
8
polypropylene (isotactic)
2*3/1
6.66, 20.78, 6.50
90, 99.6, 90
12
PET
12*1/1
4.56, 5.96, 10.75
98.5, 118, 112
Polyisoprene (cis 1,4)
8*1/1
12.46, 8.86, 8.1
90, 90, 90
8
Exercise:
Calculate the density of perfectly crystalline polyethylene, given assuming the crystal structure from the previous page.
Solution:
Density is mass/volume, so we just need to calculate the volume of a unit cell, and the total mass of material included in the unit cell:
1. Unit cell volume: abc = 9.341 × 1 0 -23 c m 3
2. Mass: 14 g/mole = 2.326 × 1 0 -23 g per C H 2 unit. Since there are 4 C H 2 units in a unit cell, mass per unit cell = 4( 2.326 × 1 0 -23 g ) =9.302 × 1 0 -23 g .
3. Density = mass/volume = 0.996 g/cm 3 .
Note that the density is close to 1 g/cm 3 , which is generally true for most polymers.

17.5 Stiffness of Polymer Single Crystals

image: 165_home_ken_Mydocs_MSEcore_331_figures_Polymer_Stiffness.png
Figure 120: Macromol. 23, 2365-2370 (1990).
Single crystal polyethylene has the following calculated values for the modulus along the 3 crystal axes:
Amorphous fraction: E=0.001 GPa
The actual modulus is almost always much lower than this because of imperfect crsytal alignment, chain folding and the presence of an amorphous fraction. See the following pages for details.

17.6 Chain Folding

image: 166_home_ken_Mydocs_MSEcore_331_figures_Chain_Folding_1.png
Figure 121: Diagram of polymer chain folding
The solid line in this animation represents the axis of a helix. As the polymer crystallizes, the helix folds back and forth on itself, producing a flat, plate-like lamellar crystallite. The chain folds are energetically unfavorable - the configuration with the lowest free energy would be a for the helices to pack in a completely extended manner, with no chain folds at all. However, the chain-folded structure is favored kinetically. Thin crystallites (low values of thecrystal thickness , λ ), are able to form much more quickly than crystallites with high values of λ .

17.6.1 Models of chain folding

image: 167_home_ken_Mydocs_MSEcore_331_figures_Chain_Folding_Models.png
Figure 122: Extreme views of chain folded crystals
These drawings illustrate two extreme views of the structure of chain folded crystals. In the adjacent reentry model, the molecules are assumed to fold back on themselves in a very well-organized way. In the random switchboard model, chains which leave the crystal at the tops and bottom surfaces reenter the crystal at random positons. Reality lies somewhere between these two extremes.

17.6.2 Amorphous Fraction and Tie Molecules

image: 168_home_ken_Mydocs_MSEcore_331_figures_Amorphous_Regions_and_Ties.png
Figure 123: Lamellar crystallites
Isolatedlamellar crystallites can be obtained from polymers crystallized from dilute solution, but melt-crystallized polymers often have structures similar to what is illustrated here. The lamellar, platelike crystallites are separated by amorphous regions. Individual crystallites are bridged bytie molecules . These tie molecules span the amorphous region, and are incorporated into adjacent crystallites.

17.6.3 Percent Crystallinity for Lamellar Crystallites

image: 169_home_ken_Mydocs_MSEcore_331_figures_Percent_Crystallinity_Diagram.png
Figure 124: Diagram of percent crystallinity of lamellar crystallites
The amorphous fraction in a semicrystalline polymer resides between the lamellar crystallites as depicted here. If the crystallites are assumed to be perfectly crystalline, the percent crystallinity is given by the ratio of thecrystal thickness ( λ c ) to thecrystal period ( d c ). The crystal period is determined from x-ray scattering measurements. The overall percent crystallintity is determined by measuring the density of the sample.

17.7 Spherulitic morphologies

image: 170_home_ken_Mydocs_MSEcore_331_figures_Spherulites.png
Figure 125: A representation of a spherulite
A spherulite originates from the radial growth of lamellar crystallites from a central nucleation point. Not all semicrystalline polymers have spherulitic morphologies, but this morphology is quite common.
image: 171_home_ken_Mydocs_MSEcore_331_figures_Spherulites_1.png
Figure 126: Formation of a spherulite
During crystallization, the spherulite grows radially outward, with each branch extending in the growth direction while maintinaing a constant crystal thickness, λ c .
image: 172_home_ken_Mydocs_MSEcore_331_figures_Spherulites_2.png
Figure 127: Geometric representation of a spherulite
This is the geometric requirement for addition of new lamellae; the additional space taken up by the spherulite as a whole is filled by new lamellae.
image: 173_home_ken_Mydocs_MSEcore_331_figures_Crystal_Growth_Vid_2.png
Figure 128: Video of the formation of a spherulite

17.8 Birefringence

image: 174_home_ken_Mydocs_MSEcore_331_figures_Birefringence_Interactive.png
Figure 129: Interactive model of birefringence
  1. Light emerging from the sample for ψ =45 º and Δ φ =90 º is said to be circularly polarized. Why is this a sensible description of the polarization state you observe under these circumstances?
  2. For what values of ψ is no light at all transmitted through the analyzer?
  3. The maximum light transmission through the analyzer is obtained for ϕ =45 º and Δ φ =180 º . What fraction of the light emerging from the sample is transmitted under these circumstances?

17.8.1 Birefringence and Radial Symmetry

image: 175_home_ken_Mydocs_MSEcore_331_figures_Birefringence_Diagram.png
Figure 130: Diagram showing birefringence
http://bly.colorado.edu/lcphysics/textures/
If the object between the polarizer and analyzer has spherical symmetry, as is the case with these liquid crystalline droplets, the'Maltesecross ' pattern shown here will be obtained. Dark patches correspond to Ψ = 0 and 90°, and the bright patches are obtained at Ψ = ± 45°.

17.9 Growth of a Lamellar Crystallite

The overall free energy change ( Δ F ) associated with the formation of a small crystallite has two components:
  1. A surface energy component, which is always positive.
  2. A bulk free energy component, which is negative for temperatures below the equilibrium melting temperature, T m 0 .
The bulk component can be written in terms of the undercooling ( T m 0 -T) by writing the bulk free energy change per unit volume ( Δ f c ) in terms of its entropic and enthalpic components:
Δ f c = Δ h c -T Δ s c (17.1)
Here Δ h c and Δ s c are the enthalpy and entropy of melting for a perfect crystal, normalized by the volume of the crystal. We can use the fact that Δ f c =0 when T= T m 0 to eliminate Δ s c from the equation and write Δ f c in terms of the enthalpy of melting and the undercooling , Δ T, defined as T m 0 -T :
Δ f c =- Δ h c Δ T T m (17.2) Note that we have used the following sign conventions (as is generally the case):
The growth of a lamellar crystallite occurs by the addition of molecular strands (typically corresponding to the c axis of the cyrstalline structure) to the growing crystallite as illustrated in Figure 131. In addition to the
image: 176_home_ken_Mydocs_MSEcore_331_figures_Secondary_Nucleation.png
Figure 131: Schematic representation of a growing lamellar crystallite.
If we account for the surface energies associated with the new crystal/amorphous interface that is created, the increase in free energy for adding n "strands" to the face of a lamellar crystallite is given by the following expression:
Δ F=2a' λ c γ s +2na'b' γ e +na'b' λ c Δ f c (17.3)
If we assume that n is large we can neglect the term involving γ s , so that Δ F<0 when Δ f c < 2 γ e λ c . With Δ f c given by Eq. 17.2, we can rewrite the condition that Δ F<0 in the following form:
λ c > 2 γ e T m 0 Δ h c Δ T (17.4)
A crystal will grow by the addition of additional strands when this criterion is met. Thick crystals (high λ c ) are thermodynamically favorable, but thin crystals (low λ c ) are kinetically favorable. In this model, the crystal thickness which is actually obtained is the lowest, thermodynamically possible value. The inequality from the previous page therefore becomes an approximate equality:
λ c 2 γ e T m 0 Δ h Δ T (17.5)
image: 177_home_ken_Mydocs_MSEcore_331_figures_Nucleation_Graph.png
Figure 132: Graphical representation of the crystal thickness prediction from Eq. 17.5.
Because thick crystals are thermodynamically more stable than thin crystals, the thick crystals will melt at a slightly higher temperature than thin crystals. One therefore expects a relationship between the temperature at which the crystals will form, and the temperature at which the crystals will melt. Crystals crystallized at higher temperatures will also melt at higher temperatures. At an undercooling of zero the actual melting temperature, T m , will be equal to the equilibrium melting temperature. For very low undercoolings, however, the crystallization rate becomes very slow, so crystallization will never take place at the exact equilibrium melting temperature. Values of the equilibrium melting temperature are obtained from measurements of T m at higher undercoolings, and extrapoloating to Δ T=0 as illustrated in the following graph:
image: 178_home_ken_Mydocs_MSEcore_331_figures_Equilibrium_Melting_Temp.png
Figure 133: Melting temp ( T m ) as a function of the crystallization temperature, T c .

17.10 Density of Semicrystalline Polymers

The percent crystallinity in a polymer sample is generally determined by measuring the density of the material. A linear relationship between density and percent crystallinity is assumed, so that the following relationship holds:
%Crystallinity= ρ - ρ amorphous ρ crystalline - ρ amorphous x100 (17.6)
image: 179_home_ken_Mydocs_MSEcore_331_figures_Semicrystalline_Density.png
Figure 134: Density vs percent crystallinity

18 Liquid Crystals

Liquid crystals are materials that that are intermediate between true traditional crystalline materials, and liquids. The molecules in a liquid crystalline phase are not organized into a full, 3- dimensionally periodic structure, but they do have some degree of orientational or positional order. Examples of molecular organization in some liquid crystalline phases are shown below.
image: 180_home_ken_Mydocs_MSEcore_331_figures_Liquid_Crystals_1.png
Figure 135: Structure of liquid crystals
http://commons.wikimedia.org/wiki/Liquid_crystal

18.1 Nematic Liquid Crystals

image: 181_home_ken_Mydocs_MSEcore_331_figures_Nematic_Liquid_Crystals.png
Figure 136: Structure of nematic liquid crystals
In anematic liquid crystal phase , the centers of the molecules have liquid-like order, but the molecules are aligned in the same direction, defined by thedirector , represented for these needleshaped and disc-shaped molecules by the red arrow.

18.2 Nematic Liquid Crystals - Texture

image: 182_home_ken_Mydocs_MSEcore_331_figures_Texture_of_Nematic_Liquid_Crystals.png
Figure 137: Image of a nematic liquid crystal
http://commons.wikimedia.org/wiki/Liquid_crystal
This photograph was taken between crossed polarizers, uing a microscope geometry like that described in the description of Section 17.8.

18.3 Liquid Crystal Displays

image: 183_home_ken_Mydocs_MSEcore_331_figures_Liquid_Crystal_Displays.png
Figure 138: Structure of an LCD
http://en.wikipedia.org/wiki/Liquid_crystal_display

19 Solutions and Blends

We now move from our discussion of the shapes of individual polymer molecules in a homogeneous, one component system, to the behavior of solutions and blends consisting of different types of molecules that have been mixed together. We are interested in studying these systems for many reasons, including the following:
  1. Characterization: many of the properties of polymer molecules are determined by the properties of polymer solutions.
  2. Processing: many applications of polymers involve require that the polymer molecules be dissolved in an appropriate solvent.
  3. Improved materials: new materials with improved properties can often be formed by blending or "alloying" different polymers.
  4. Recycling: polymers need to be separated before they can be recycled because the different polymers do not mix favorably.
The key quantity when discussing thermodynamics of polymer solutions (polymer + small-molecule solvent) and polymer blends (polymer A + polymer B) is the free energy of mixing. We therefore begin with a discussion of polymer solution thermodynamics. In general, we consider the case where we mix n a A molecules with a degree of polymerization N a with n b B molecules with a degree of polymerization N b . Polymer solutions will correspond to the case where one of the degrees of polymerization is very small (typically 1). We will make the assumption that the sizes of the repeat units for the components of the mixture are identical. This assumption is not as restrictive as it would initially appear, since one can always define "effective" repeat units with the desired volume. In this sense, N a and N b are determined by the relative molecular volumes of the components in the mixture, and are defined in terms of some reference volume.
Exercise:
The room temperature densities of polystyrene and toluene are 1.05 g/c m 3 and 0.87 g/c m 3 , respectively. What reference volume, V 0 , should we use to characterize the thermodynamics of polystyrene solutions in toluene? What degrees of polymerization would characterize a solution where the polystyrene molecular weight is 100,000 g/mole?
Solution:
The obvious reference volume in this case is the volume of a toluene molecule, which we take as having N=1. This volume is obtained from the density and molecular weight of toluene. Toluene ( C 6 H 8 ) has a molecular weight of 92 g/mole.
V 0 = ( 92g mole ) ( c m 3 0.87g ) =106c m 3 /mole
Note that we have used the molar volume for V 0 . This convention will be followed throughout our discussion.
The value of N for the polymer molecules will be given by the molecular volume, divided by the reference volume:
N polystyrene = ( 100,000g mole ) ( c m 3 1.05g ) ( mole 106c m 3 ) =898
In our discussion of the thermodynamics of mixtures, we will refer to N as the number of "segments", defined in terms of a reference volume as described here. In some cases this reference volume will correspond to a repeat unit of the polymer, but this does not need to be the case in general.

19.1 Chemical Potentials and Free Energy of Mixing

The chemical potential of A molecules, μ a , is defined as the free energy required to add an additional 'A' molecule to the system. Because all A molecules are equivalent, it is also equal to the free energy per A molecule in the system. The following expressions are helpful
F= n a μ a + n b μ b (19.1)
. F n a | n b = μ a . F n b | n a = μ b (19.2)
The derivative with respect to n a is taken with n b held constant, and vice versa. Note that F is the total extensive free energy of the system, since it increases as the numbers of A and B molecules ( n a and n b ) increases.
image: 184_home_ken_Mydocs_MSEcore_331_figures_Free_Energy_Mixture_2.png
Figure 139: Animation of the chemical potential of a molecule
The chemical potentials of A and B molecules of the pure components are defined as u a 0 and u b 0 . The free energy of the two components before mixing is therefore given by:
F unmixed = n a μ a 0 + n b μ b 0 (19.3)
The free energy of mixing is equal to the difference in the free energies of the mixed and unmixed systems:
Δ F mix =F- F unmixed = n a μ a + n b μ b -( n a μ a 0 + n b μ b 0 ) = n a ( μ a - μ a 0 ) + n b ( μ b - μ b 0 ) (19.4)
Because we are only interested in changes in chemical potentials due to mixing, we can set the chemical potentials of the pure, unmixed components to zero ( μ a 0 = μ b 0 =0 ).
By defining the chemical potentials of the pure components as zero, we obtain:
Δ F mix = n a μ a + n b μ b (19.5)
Equation 19.5 is completely general, and valid for mixtures of molecules of any size. It's not that useful, however, because we generally don't write down values of n a and n b when we want to specify the composition of a multicomponent mixture. For mixtures of small molecules (metal alloys, for example), we commonly use mole fractions, X a and X b , to specify the composition:
X a = n a n a + n b X b = n b n a + n b (19.6)
This makes sense for metals because the molar volumes of the different components are not that different from one another. In polymer systems, however, we are often dealing with systems where the components have molar volumes that can differ by orders of magnitude. This is the situation if we are interested in the solubility of a high molecular weight polymer (which can have a molecular weight of several hundred thousand g/mole) in a solvent (which will typically have a molecular weight of 100 g/mole). Organic solutions are typically specified by the weight fractions of the different components. In discussions of polymer solution thermodynamics it is more common to work in terms of volume fractions of the different components in the systems.
The concepts in this section can be easily extended to systems with more than two components, but we confine ourselves here to binary solutions of just two components: A molecules of length N a and B molecules of length N b . In our discussion of polymer thermodynamics it is important to keep in mind N a N b are not necessarily true degrees of polymerization as defined earlier in this text. Instead they are obtained by dividing the molar volumes of the different components by some reference volume, V 0 . For polymer/solvent mixtures this reference volume is typically taken as the molar volume of the solvent, and for polymer/polymer mixtures it is typically taken as the volume per repeat unit for one of the polymers. The total volume, V , and volume fractions, φ a and φ b are given by the following expressions:
V= n a N a V 0 + n b N b V 0 (19.7)
φ a = n a N a V 0 V φ b = n b N b V 0 V (19.8)
The free energy of mixing per unit volume, Δ f mix , is obtained from Δ F mix by dividing by the volume of the system:
Δ f mix = n a μ a + n b μ b n a N a V 0 + n b N b V 0 = μ a φ a N a V 0 + μ b φ b N b V 0 (19.9)
By rearranging the equation at the bottom of the previous page we obtain:
Δ f mix V 0 = μ a φ a N a + μ b φ b N b (19.10)
The quantities Δ f mix , μ a / N a and μ b / N b all depend on the composition of the material. The relationship between these quantities is illustrated by the tangent construction shown in Figure 140, where the composition dependence of Δ f mix is plotted as a function of φ b for a hypothetical material. For a given composition (denoted φ b in Figure 140), the chemcial potentials are obtained by first drawing a tangent to the free energy curve at φ b = φ b . The value of μ a / N a is determined by extrapolating this tangent line to φ b =0 and μ a / N a is obtained by extrapolating the tangent line to φ b =1. This result is consistent with Eq. 19.10, which shows that normalized free energy of mixing at a certain composition is determined by the weighted average of the normalized chemical potentials. As a result a straight line drawn between the normalized chemical potentials on the left and right axes must hit the free energy curve at φ b = φ b .
To fully validate the tangent construction illustrated in Figure 140 we still need to show that the straight line between the chemical potentials must be equal to the slope of the free energy curve. This can done in formal terms by combining Eqs. 19.2 and 19.8, which after simplification leads to the following relationship:
( Δ f mix V 0 ) φ b = μ b N b - μ a N a =- ( Δ f mix V 0 ) φ a (19.11)
The fact that the tangent must give the difference between μ b / N b and μ a / N b makes sense conceptually, because the only way changing the composition by some amount requires that an A segment be replaced by a B segment.
image: 185_home_ken_Mydocs_MSEcore_331_figures_Solution_Thermodynamics.png
Figure 140: Relationship of the free energy of mixing to the chemical potentials.

19.2 Ideal Entropy of Mixing for Polymers

For small molecules the ideal free energy of mixing is given by the following expression:
Δ s ideal V 0 R =- X a ln X b - X b ln X b (19.12)
The quantity s ideal V 0 is the entropy of mixing per mole of molecules. If the molar volumes of the A and B molecules are the same, than we can simply replace the mole fractions, X a and X b with the corresponding volume fractions, φ a and φ b . The entropy of mixing per molecule is independent of the size of the molecule. Therefore, the entropy of mixing per repeat unit is smaller than the entropy of mixing per molecule by a factor of N , where N is the degree of polymerization. In general, we have the following expression for the ideal free energy of mixing per volume for polymers with molar volumes of N a V 0 and N b V 0 :
s ideal V 0 R =- φ a ln φ a N a - φ b ln φ b N b (19.13)

19.3 Idealized Enthalpy of Mixing

The key result is that the entropy of mixing for a given volume is much smaller for large molecules than it is for small molecules. This result explains why different types of polymers almost never mix with one another in the liquid state. The driving force for two liquids to mix is normally the entropy of mixing, but this driving force is very low for mixtures of very large polymer molecules. In fact, the entropy of mixing for polymer mixtures is sometimes negative. While the ideal entropy of mixing must be positive, it is very small, and can be overwhelmeed by non-ideal contributions to the entropy of mixing which can be positive or negative. Before discussing these effects, however, we will continue in our derivation of the Flory Huggins equation for the free energy of mixing. The Flory-Huggins theory combines the ideal free energy of mixing with a very simple form for the enthalpy of mixing. Unlike the entropy of mixing, the enthalpy of mixing is assumed to be independent of the size of the molecules. We can therefore use small molecules to illustrate the origins of the equation for the enthalpy of mixing.
image: 186_home_ken_Mydocs_MSEcore_331_figures_Mixing_Simulation_2.png
Figure 141: An animated simulation of mixing
image: 187_home_ken_Mydocs_MSEcore_331_figures_Enthalpy_of_Mixing.png
Figure 142: Pictoral representations of mixing
In this example, the blue square is placed into an environment where one of its four nearest neighbors is blue, and three are red. Red-blue contacts have an energy (or enthalpy) of E ab , red-red contacts have an energy of E aa , and blue-blue contacts have an energy of E bb . The enthalpy change, Δ H, associated with the removal of the square from a pure blue phase on the left (4 blue neighbors) to the mixture on the right (1 blue neighbor and 3 red neighbors) is therefore given by the following formula:
Δ H=3 E ab + E bb -4 E bb =3( E ab - E bb ) (19.14)
image: 188_home_ken_Mydocs_MSEcore_331_figures_Enthalpy_of_Mixing_II.png
Figure 143: Diagram showing contacts
How do we calculate the energy of a randomly mixed collection of A and B molecules? We begin by assuming that the energy of an A-A contact is E aa , the energy of an B-B contact is E bb , and the energy of an A-B contact is E ab . Also assume that each molecule z nearest neighbors ( z=4 for the 2-dimensional array of squares shown above.) The number of A-B contacts in a mixture of n a A molecules and n b B molecules is given by:
n aa = n a z p a /2 n bb = n b z p b /2 n ab = n a z p b (19.15)
where p a is the probability that a nearest neighbor is an A molecule. The factor of 2 in the expressions for n aa and n bb is needed to avoid double counting. (For example, if n a =2, there is only one interaction between the two molecules, not 2.) Now we invoke the mean field (or random mixing) approximation, which is that p a = φ a , and p b = φ b . We obtain the follwing.
n aa = n a z φ a /2 n bb = n b z φ b /2 n ab = n a z φ b (19.16)
The energy of the mixed state is obtained by multiplying the numbers of different contacts by their energy:
E mixed = E aa n aa + E bb n bb + E ab n ab (19.17)
Substitution for n aa , n bb and n ab gives the following:
E mixed = E aa n a z φ a /2+ E bb n b z φ b /2+ E ab n a z φ b (19.18)
The energy of the unmixed state is obtained in a similar fashion, remembering that in this state all A's are surrounded by other A's, and all B's are surrounded by other B's:
E unmixed = E aa n a z/2+ E bb n b z/2 (19.19)
The energy change due to mixing is obtained by subtracting the energy of the mixed state from the energy of the mixed state:
Δ E mix = E mixed - E unmixed = E aa n a z φ a 2 + E bb n b z φ b 2 + E ab n a z φ b - E aa n a z 2 - E bb n b z 2 (19.20)
This expression can be further simplified by recalling that n a = φ a n and n b = φ b n , where n is the total number of molecules ( n= n a + n b ):
Δ E mix =nz[ E aa φ a 2 2 + E bb φ b 2 2 + E ab φ a φ b - E aa φ a 2 - E bb φ b 2 ] (19.21)
After some rearrangement we obtain the following expression:
Δ E mix =nz φ a φ b [ E ab - E aa + E bb 2 ] (19.22)
Note that the mixing energy is positive if the energy of a-b contacts is larger than the average energy of a-a and b-b contacts. Because z, E aa , E bb , and E ab are actually not measurable quantities, it is common to lump them together into a parameter called the
Chi parameter
represented by the Greek letter 'chi':
χ z k B T [ E ab - E aa + E bb 2 ] (19.23)
With this definition of χ , the energy of mixing for A and B molecules (with the same size) is:
Δ H mix = Δ E mix =RT χ n φ a φ b (19.24)
Strictly speaking, the 'enthalpy' and the 'energy' of a system are not the same thing. The Gibbs free energy involves the enthalpy, and is the appropriate thermodynamic quantity at a fixed temperature and pressure. The Helmholtz free energy involves the energy, and is the appropriate thermodynamic quantity at fixed temperature and volume. For compressible systems the two quantities differ slightly. We use the expression from the previous page for the extensive enthalpy of mixing, and divide by the volume of the system ( n V 0 ) to obtain the intensive enthalpy of mixing:
Δ h mix = Δ H mix n V 0 = RT χ φ a φ b v 0 (19.25)
Normalizing as we did for the ideal entropy of mixing gives:
Δ h mix V 0 RT = χ φ a φ b (19.26)

19.4 Flory-Huggins Free Energy of Mixing

The most commonly used expression used to describe the free energy of mixing for polymer/polymer and polymer/solvent mixtures combines mean-field enthalpy of mixing with the ideal combinatorial entropy of mixing:
Δ f mix = Δ h mix -T Δ s ideal (19.27)
Combining these forms for the enthalpy and entropy gives the following expression, referred to as the Flory-Huggins free energy of mixing:
Δ f mix V 0 RT = φ a ln φ N a + φ b ln φ b N b + χ φ a φ b (19.28)
It is important to recognize that many assumptions went into the derivation of this equation, and that they can't possibly all be true. Nevertheless, this equation is very commonly used in polymer science, and needs to be understood. In fact, this equation is often used as a definition for χ . Based on our earlier definition of χ in terms of energetic parameters, one would expect that it's magnitude would decrease as the temperature increases. In fact, the actual free energy of mixing can have very different behavior, and the free energy of mixing is generally accurately described by the Flory- Huggins expression only when χ is allowed to vary with temperature and composition. We begin on the next page by constructing phase diagrams, with the assumption that the free energy of mixing is adequately described by a value of χ that is independent of composition. Free energy curves for N a = N b =1 are shown in Figure 144 for χ ranging from 0 to 4. For N a = N b =1, the critical value of χ is equal to 2. When χ is larger than this critical value there will be a region of compositions where the free energy curve has a negative curvature in the free energy curve.
image: 189_home_ken_Mydocs_MSEcore_331_figures_floryplotsna1nb1.png
Figure 144: Regular solution form or for the free energy of mixing of small molecules.

19.5 The Coexistence (binodal) and Spinodal Curves

The critical value of χ decreases as N a and N b increase. In addition, the free energy curves are no longer symmetric about φ b =0.5 when N a N b . Both of these points are illustrated in Figure 145, where we plot the Flory-Huggins free energy of mixing (Eq. 19.28) for N a =200 , N b =500 and χ = 0.012. A two-phase region will exist for average compositions of the alloy that are between the compositions between φ 1 and φ 2 , the two compositions defined by the common tangent construction in Figure 145. Because the tangent to the curve at these two compositions are identical, they have the same values of μ a and μ b, which is a requirement for equilibrium. In addition, the average free energy for a mixture of these two phases (dashed line in Figure 145) is less than the free energy for a homogenous phase of this average composition (solid line in Figure 145). By plotting these two composition as a function of χ , we obtain the binodal curve shown as the solid line in Figure 146. The binodal curve is also referred to as the coexistence curve because it shows the compositions of the phases that are in equilibrium with one another in the two-phase region of the phase diagram.
image: 190_home_ken_Mydocs_MSEcore_331_figures_floryplotsna200nb500.svg
Figure 145: Free energy of mixing for a system with asymmetric molecular sizes in the two phase regime.
image: 191_home_ken_Mydocs_MSEcore_331_figures_Na200Nb300chidiagram.png
Figure 146: Binodal and spinodal curves for N a =200 and N b =300 .
The spinodal curve (or stability limit ) represents the boundary in the phase diagram where the free energy of mixing changes from a positive curvature to a negative curvature. It occurs when χ = χ s , where χ s defines the spinodal curve. For N a =200 and N b =300 the spinodal curve is shown as the dashed line in Figure 146. The change in curvature from positive to negative values occurs when the second derivative of the free energy with respect to composition is equal to zero. We obtain the following expression for χ s by setting φ b =1- φ a in the free energy expression:
2 φ a 2 ( φ a ln φ a N a + ( 1- φ a ) ( ln 1- φ a ) N b + χ s φ a ( 1- φ a ) ) =0 (19.29)
Evaluating the second derivative leads to the following:
1 φ a N a + 1 ( 1- φ a ) N b -2 χ s =0 (19.30)
which can be rearranged to give the following for χ s :
χ s = 1 2 φ a N a + 1 2 φ b N b (19.31)

19.6 Critical Point

The phase diagrams obtained from the free energy of mixing expression show under what conditions polymeric liquids are completely miscible, and under what conditions phase separation will occur. The critical point is the point on the phase diagram where a liquid is just beginning to undergo phase separation. Mathematically, the critical point corresponds to the point on the spinodal curve where χ s has the lowest possible value. We obtain this point by setting the derivative of χ s to zero:
φ a ( 1 2 φ a N a + 1 2( 1- φ a ) N b ) =0 (19.32)
The solution to this equation is given by φ a = φ a,crit , at which point χ s is equal to χ crit :
φ a,crit = N b N b + N a (19.33)
χ crit = N a + N b +2 N a N b 2 N a N b (19.34)
Exercise:
The value of for blends of polystyrene (PS) and poly(methyl methacrylate) (PMMA) is 0.037 at 170 °C, based on a reference volume equal to the reference volume of a polystyrene repeat unit. Suppose monodisperse samples of PS and PMMA of equal molecular weights are mixed with one another. What range of molecular weights will be completely miscible with one another at 170 °C? (Assume the densities of the two polymers are identical).
Solution:
First we calculate χ crit :
χ crit = N a + N b +2 N a N b 2 N a N b = 2 N (19.35) for N a = N b =N
The polymers are completely miscible for χ < χ crit , which occurs when N < 2/ χ . With χ = 0.037, we have complete miscibility for N < 54. The repeat unit molecular weight for polystyrene is 104 g/mol, so the critical molecular weight below which the polymers are completely miscible is 54*104 g/mol = 5616 g/mol.

19.7 Spinodal Decomposition

"Spinodal Decomposition" is the mechanism of phase separation where small perturbations in the composition of the single phase grow because this decreases the overall free energy. This happens inside the spinodal curve defined previously. The microstructure is generally characterized by a well defined wavelenth, as illustrated below.
image: 192_home_ken_Mydocs_MSEcore_331_figures_Spinodal_Decomposition.png
Figure 147: Characteristic microstructures for phase separation by spinodal decompositon (left) and nucleation and growth (right).

19.8 Phase Diagrams in Temperature Space

Actual values of χ for real systems are generally obtained by fitting the Flory Huggins expression to the free energy of mixing to actual experimental data. This procedure can give actual values of χ which either increase or decreaes with increasing temperature. Smaller values of χ (where mixing is favored) can therefore be obtained at either high or low temperatures, giving the following two types of phase diagrams:
image: 193_home_ken_Mydocs_MSEcore_331_figures_Phase_Diagrams.png
Figure 148: Phase diagrams in temperature space.
Exercise:
The following equation describes the temperature dependence of χ for the PMMA/butanol polymer/solvent system:
χ =1.45-0.0115T (19.36)
Here T is expressed in C. For what range of temperatures will high molecular weight PMMA completely dissolve in butanol?
Solution:
From the phase diagram, we see that the two components are completely miscible for χ < χ crit . If the A phase is the solvent ( N a =1 ) we obtain the following for :
χ crit = 1+ N b +2 N b 2 N b (19.37)
A cautionary note: Note that the polymer volume fraction at the critical point is very low (i.e., is very close to 1 for N b >> N a ). The detailed predictions of the Flory Huggins theory will not be accurate in this regime, for reasons that are discussed more fully in the section on dilute solution thermodynamics.

19.9 Chemical Potentials

Recall that the Section 19.1 of species A is determined by differentiating the extensive free energy with respect to the number of A molecules in the system. The extensive free energy is obtained by multiplying the free energy density, f, by the volume of the system, V. For a binary mixture of A and B polymers we have:
Δ f mix = RT v 0 { φ a ln φ a N a + φ b ln φ b N b + χ φ a φ b } V (19.38)
Using the expressions for φ a and φ b (Eq. 19.8) and V (Eq. 19.7) we have:
Δ F mix RT ={ n a ln φ a + n b ln φ b + χ n a N a φ b } (19.39)
We are now in a position to obtain expressions for the chemical potentials by differentiating the extensive free energy of mixing. We make the usual assumption that the chemical potentials of the pure components are equal to zero (activities of the pure components equal to one):
μ a = Δ F mix n a =RT n a ( n a ln φ a + n b ln φ b + χ n a N a φ b ) (19.40)
In this way we obtain the following result:
μ a RT = ln φ a + φ b ( 1- N a N b ) + χ N a φ b 2 (19.41)
When discussing polymer/solvent mixtures we replace the 'a' and 'b' subscripts with 'p' (for polymer) and 's' (for solvent). In addition, we take N s = 1. Formally, N p is the ratio of the polymer volume to the solvent volume. With these substitutions, we obtain the following expression for μ s , the solvent chemical potential, which is valid for all concentrations:
μ s RT = ln ( 1- φ p ) + φ p ( 1- 1 N p ) + χ φ p 2 (19.42)

19.9.1 Limiting case for high molecular weight polymer

For N p the solvent chemical potential reduces to μ s :
μ s RT = ln ( 1- φ p ) + φ p + χ φ p 2 (19.43)
In addition the critical point from Eqs. 19.33 and 19.34 occurs at φ p =0 and χ =0.5. This means that for χ >0.5 pure solvent will be in equilibrium with a phase that contains both polymer and solvent. The polymer volume fraction in this phase is determined from the requirement that the solvent chemical potential must be zero. In other words, if a chunk of high molecular weight polymer is immersed in a solvent, solvent diffuses into the polymer until φ p decreases to the point where μ s =0. In general, this solution must be solved numerically. In Figure 149 we plot the concentration dependence of μ s for χ =0 , χ =0.5 and χ =1 .
image: 194_home_ken_Mydocs_MSEcore_331_figures_musinf.svg
Figure 149: Concentration dependence of the solvent chemical potential from Eq.19.43.

19.9.2 Limiting case for dilute solutions

For small values of φ p the log term can be expanded as follows:
ln ( 1- φ p ) =- φ p - φ p 2 2 -... (19.44)
Retaining only the first two terms in the expansion gives the following for the solvent chemical potential:
μ s RT =- φ p N +( χ -1/2 ) φ p 2 (19.45)
Exercise:
What is the polymer volume fraction for a polymer solution with a concentration of 0.05 g/cm 3 ? Assume the solvent has a density, ρ s , of 0.85 g/cm 3 and the polymer has a density, ρ p , of 1.1 g/cm 3 and there is no volume change on mixing.
Solution:
We just need to calculate the total weight of solution in a volume of 1 c m 3 :
c= φ p ρ p =0.05(1.1 g/c m 3 )=0.055g/c m 3 (19.46)
Note: Actual solution concentrations are generally reported as weight of polymer per total volume of solution - the quantity that we refer to as c. From a theoretical standpoint, it is generally more convenient to use the polymer volume fraction, f c . These quantities are linearly related to one another through the polymer density, as illustrated here. The solvent density is not needed for this calculation, although it is needed to determine the ratio of polymer to solvent.

19.10 Osmotic Pressure

The concept of osmotic pressure can be illustrated by considering the operation of the membrane osmometer illustrated conceptually in Figure 150. The device consists of a polymer/solvent mixture that is separated from a bath of pure solvent by a membrane that is permeable to solvent but impermeable to the polymer. The chemical potential of the pure solvent is zero, by definition. Because the solvent chemical potential in the polymer/solvent mixture is less than zero, there is a driving force for solvent to move into this side of the membrane from the pure solvent side. Solvent will continue to move into this side of the device until the hydrostatic pressure generated by the increased depth of the liquid layer provides counteracts the driving force associated with the chemical potential balance. This pressure is the osmotic pressure, and is obtained by dividing the chemical potential difference by the solvent volume:
Π =- μ s / V s (19.47)
image: 195_home_ken_Mydocs_MSEcore_331_figures_Osmotic_Pressure.png
Figure 150: Diagram demonstrating the principles of osmotic pressure.
One important aspect of an osmotic pressure measurement is that it can be used to provide a direct measure of the number average molecular weight of the polymer. Osmotic pressure effects are important in many other situations as well. For example, the osmotic pressure is the minimum pressure required to purify a solution by forcing it through a membrane. The osmotic pressure of aqueous solutions is also important in a wide variety of situations in biology. The height of many trees is limited by the maximum osmotic pressure that can be sustained in the leaves, since this is the driving force that draws water to the upper reaches of the tree. Here we explore these phenomena in more detail, and develop some quantitative descriptions of the osmotic pressure.
By combining the Flory-Huggins expression for the solvent chemical potential (Eq. 19.42) with the definition of the osmotic pressure (Eq. 19.47), we obtain the following expression for Π :
Π = RT V s ( - ln ( 1- φ p ) - φ p ( 1- 1 N p ) - χ φ p 2 ) (19.48)
If N p is large the osmotic pressure depends only on the polymer concentration, and not on the degree of polymerization. The following result is obtained for N p >>1:
Π = RT V s ( - ln ( 1- φ p ) - φ p - χ φ p 2 ) (19.49)

19.10.1 Osmometry for Molecular Weight Determination

The number average molecular weight can be obtained by measuring the osmotic pressure of a dilute polymer solution. Using the dilute solution form of the Flory-Huggins expression for the solvent chemical potential (Eq. 19.45) gives the following for Π :
Π RT = φ p N p V s +( 1/2- χ ) φ p 2 V s (19.50)
The polymer volume fraction and overall solution volume (V) are given by the following expressions:
φ p = n p N p V s V (19.51)
Use of this expression for φ p gives the following for π :
Π RT = n p V +( 0.5- χ ) N p 2 V s ( n p V ) 2 +... (19.52)
Here it is helpful to rewrite n p /V in the following way:
n p V =( w V ) ( n p w ) = C M n (19.53)
where c is simply the polymer concentration (kg/m 3 in SI units, g/cm 3 in the units that are more commonly used) and M n is the number average molecular weight (kg/mole in SI units). The osmotic pressure is therefore related to the number average molecular weight and concentration as follows:
Π RT = C M n +( 0.5- χ ) N p 2 V s ( 1 M n ) 2 C 2 +... (19.54)
This is expression is often written as a virial expansion in powers of the concentration:
Π RT = C M n + A 2 C 2 + A 3 C 3 +... (19.55)
A 2 is the second virial coefficient , which describes effective two-body interactions between polymer segments within the solvent. A 2 = 0 for χ =0.5 , which means that the polymer chain does not 'see' itself under these conditions. For χ =0.5 , the polymer chain behaves as a random walk, with chain dimensions very similar to what is obtained for the pure, amorphous polymer in the absence of solvent. For χ >0.5 , polymer/polymer contacts are more favorable than polymer/solvent contacts. Under these conditions, high molecular weight polymers will not dissolve.

19.10.2 Scaling Theory of Osmotic Pressure

The previous derivation of the expressions for the osmotic pressure are useful, but they are not packed with physical insight regarding the actual meaning of osmotic pressure. Some simple physical arguments are more useful in this sense. We can start by making an analogy to the ideal gas law, which must be valid in the dilute regime. We start with the familiar ideal gas law:
PV=nRT (19.56)
If we replace the pressure, P , with the osmotic pressure, Π , and replace the number of moles of gas molecules, n, with the number of polymer molecules n p . By equating n p /V with C/ M n (Eq. 19.53), we recover the leading term in virial expansion for the osmotic pressure (Eq. 19.55):
Π = CRT M n (19.57)
In the dilute regime, molecular collisions are rare, and the pressure and free energy are dominated by the entropic penalty associated with confinement of the molecules into a fixed volume. Virial coefficients quantify the deviations from this ideal, dilute limit.
image: 196_home_ken_Mydocs_MSEcore_331_figures_Overlap_Concentration.png
Figure 151: Polymer solution at a concentration less than the overlap concentration, c * . At c * the spheres overlap and completely fill the space.
As the solution gets more and more concentrated, we eventually get to the point where the individual polymer molecules overlap with one another. If we treat the polymers as spheres with a volume equal to R 0 3 , the concentration at which this happens, defined as the
overlap concentration
, c * , is given by the average concentration in one of these spheres:
C * ( g/c m 3 ) = M ( g/mole ) R 0 3 ( c m 3 ) N av ( mol e -1 ) (19.58)
Below the overlap concentration, we can specify the polymer concentration by ξ , the average distance between polymer molecules in solution, as illustrated in Figure 152. The concentration of polymer molecules (molecules per volume) is 1/ ξ 3 and the osmotic pressure is given by the following expression:
Π = RT ξ 3 (19.59)
image: 197_home_ken_Mydocs_MSEcore_331_figures_Correlation_Length_for_Dilute_Solutions.png
Figure 152: Correlation length and dilute solutions
Exercise:
What is the overlap concentration poly(methyl methacrylate) with M = 100,000 g/mole in toluene?
Solution:
We can combine Eq. 19.58 for the overlap concentration with Eq. 13.18 for R 0 to obtain the following:
C * = M R 0 3 N av = M a s 3 N 1.8 N av = M 0 a s 3 N 0.8 N av (19.60)
Note that we have used the self-avoiding walk expression for R 0 , which is the relevant expression for good solvent conditions ( χ <0.5). PMMA has M 0 = 100 g/mole, so N=1000 for the polymer in our example. With these values of N and M 0 , and an assumed value of 7 Å ( 7 × 1 0 -8 cm) for a s , we obtain c * =1.9 × 1 0 -3 g/c m -3 . The density of toluene is slightly less than 1 g/cm 3 , so this corresponds to a polymer weight fraction of a couple tenths of a percent.
The basic assumption of scaling theory of osmotic pressure is that Π is still given by Eq. 19.59 at concentrations above the overlap concentration, but that ξ has a different meaning in this regime. The approach works in the 'semidilute' solution regime where the overall concentration is still low ( φ p < 0.25), but where there is substantial overlap between polymer molecules. Each overlap point represents a constraint that contributes an entropic contribution of ~ k B T to the free energy of the system. The osmotic pressure is therefore ~ k B T times the concentration of overlap points, or k B T/ ξ 3 , where ξ is the correlation length , or average distance between overlap points. As the solution gets more concentrated, the average distance between polymer contacts decreases., as shown schematically in Figure 153.
image: 198_home_ken_Mydocs_MSEcore_331_figures_Correlation_Length_and_Concentration.png
Figure 153: Change in correlation length with solution concentration.
To get the concentration dependence of the correlation length, we can start with the postulate that the corrections to the low-conentration form of the osmotic pressure must depend only on c/ c * , so that the osmotic pressure has the following form:
Π M CRT =1+f( C/ C * ) (19.61)
Here f( C/ C * ) is some function that we still need to figure out. It's easy to sort out what this function must look like, because for C C * we don't expect that the molecular weight of the polymer to matter any more. At these concentrations the osmotic pressure is dominated by polymer/polymer contacts, and the solution at this local level does not depend at all on the molecular weight. The only way for the osmotic pressure to be independent of the concentration for high concentrations is for f( C/ C * ) /M to be independent of the M . Because C * M -4/5 , this is only possible if f( C/ C * ) ( C/ C * ) 5/4 . This means that for C> C * we have Π C 9/4 , or equivalently, Π φ p 9/4 . The correlation length, which is related to the osmotic pressure through Eq. 19.59, is then given by the following expression for C> C * :
ξ = a 0 φ p -3/4 (19.62)
Here a 0 is a typical monomer size, comparable in magnitude to the cube root of the molecular volume, or to the statistical segment length. This value for ξ can then be combined with Eq. 19.59 to give the following expression for the osmotic pressure in the semidilute concentration regime:
Π = k B T a 0 3 φ p 9/4 (19.63)
The full expression for the osmotic pressure for dilute and semidilute concentrations is obtained by using this value of ξ for C> C * and adding the osmotic pressure that dominates in the dilute limit:
Π = CRT M + k B T a 0 3 φ p 9/4 (19.64)
While derived in completely different ways, Eqs. 19.49 and 19.63 actually behave similarly in the concentration regime where they both apply. This is illustrated in Figure 154., where these two expressions are compared for a specific combination of a 0 and χ .
image: 199_home_ken_Mydocs_MSEcore_331_figures_osmoticpressureplots.svg
Figure 154: Concentration dependence of the Flory-Huggins form for the osmotic pressure (Eq. 19.49) with the semidilute scaling form (Eq. 19.63), with χ =0.4 and a 0 3 =4 v s .

19.10.3 Typical Magnitude of the Osmotic Pressure

The fundamental pressure scale for water is the thermal energy, RT divided by the solvent molar volume, V s . For water, we have (at 25 ºC):
V s = 18g mole 1 0 -6 m 3 1g =1.8x1 0 -5 m 3 =18c m 3 (19.65)
RT V s = ( 8.314 J/mole-K ) ( 298 K ) 3.0 × 1 0 -29 m 3 =1.4 × 1 0 8 Pa (19.66)
This pressure is relatively large. The actual osmotic pressure is obtained by multiplying by a prefactor involving the Flory-Huggins interaction parameter and the polymer volume fraction. Osmotic pressures of many atmospheres (1 atm = 1 0 5 Pa) can easily be obtained.
Exercise:
What is the osmotic pressure of seawater?
Solution:
The osmotic pressure arises from all of the dissolved salts in seawater. Typical concentrations are tabulated below.[5] The total concentration of all dissolved salts is 1.12 moles/kg, which after multiplying by the density of seawater (1025 kg/m 3 ) gives an overall salt concentration of 1.15 moles/m 3 . The osmotic pressure is obtained by multiplying this concentration by RT :
Π =( 8.314 J/molK ) ( 300 K ) ( 1025 kg/ m 3 ) =2.6x1 0 6 Pa
Species
moles/kg
H 2 O
53.6
C l -
0.546
N a +
0.469
M g 2+
0.0528
S O 4 2-
0.0282
C a 2+
0.0103
K +
0.0102
B r -
0.000844
S r 2+
0.000091
F -
0.000068
Total Ions
1.12

19.11 Equilibrium Swelling of a Neohookean Material

Suppose a crosslinked material is immersed in a solvent (or exposed to solvent vapor) so that it swells isotropically, increasing it's volume from an initial volume of V dry to a final volume of V wet . The amount of swelling that we get is determined by the concentration of network strands, which also determines the elastic modulus. So by measuring the swelling we can also obtain a reasonably accurate estimation of the shear modulus of the dry, unswollen polymer, which we refer to here as G dry . The modulus of solvent-swollen material, which we'll call G wet is also of interest in many applications. We'll start by assuming that the material was crosslinked at equilibrium, with the network strands obeying random walk statistics, so that the strain energy density is given by Eq. 15.10, with β =1 and G dry = k B T ν el :
Δ f d =( V dry V wet ) G dry { λ x 2 + λ y 2 + λ z 2 -3 } 2 (19.67)
The The factor of V dry / V wet accounts for the fact that as the volume of the gel increases due to solvent swelling, the concentration of network strands decreases accordingly. The extension ratios appearing in Eq. 19.67 are referenced to the undeformed, dry elastomer. To calculate the shear modulus of the solvent-swollen elastomer we need to consider a deformation process that occurs in two steps: istropic expansion of the gel to an equilibrium solvent-swollen state and shear deformation of this isotropically deformed gel.

19.11.1 Equilibrium Swelling of the Gel

The geometry of an experiment used to measure the solvent swelling is shown in Figure 155. The polymer sample of interest is placed in a sealed container containing wither pure solvent (solvent activity, a s , equal to 1), or a solvent diluted with soluble polymer to give a s <1. When water is the solvent the solvent activity is simply the relative humidity. Elasticity of the gel introduces an elastic pressure, p el that limits the amount of solvent that is able to diffuse into the gel. To calculate this pressure we begin with the deformation free energy given in Eq. 19.67, with the following values for the extension ratios and relative volumes of the dry and wet states:
λ x = λ y = λ z λ s (19.68)
V dry / V wet = φ p = λ s -3 (19.69)
Here φ p is the volume fraction of polymer in the swollen gel.
image: 200_home_ken_Mydocs_MSEcore_331_figures_Neohookean_swelling.svg
Figure 155: Experimental geometry for obtaining the equilibrium swelling.
A general expression for the equilibrium swelling is obtained by writing the deformation free energy for the swelling (step 1 in the previous subsection) in terms of the volume:
Δ F d = V wet Δ f d = 3 V dry 2 G dry { ( V wet V dry ) 2/3 -1 } (19.70)
The elastic pressure obtained directly from this expression by differentiating with respect to the sample volume:
p el = Δ F d V wet = G dry ( V dry V wet ) 1/3 = G dry φ p 1/3 (19.71)
The equilibrium polymer volume fraction in the gel is obtained by equating the solvent chemical potential in the solvent bath with the solvent chemical potential inside the gel, accounting for the elastic pressure, p el . We illustrate the point here by using the Flory-Huggins expression for the solvent chemical potential, with N p =, since the polymer is crosslinked (see Eq. 19.43):
μ s =RT ln a s =RT{ ln ( 1- φ p ) + φ p + χ φ p 2 } + p el V s (19.72)
The last term accounts for the fact that the pressure is larger inside the gel than outside the gel. This term involves the shear modulus of the dry material, equal to RT multiplied by the concentration of elastic strands in the network (Eq. 15.15, where we take β =1 ). This concentration can be expressed in terms of the average molar volume of a network strand, V el . We further define N el , an effective degree of polymerization of an elastic strand, as the ratio of V el to V s :
G dry = ν el RT= RT V el = RT N el V s (19.73) Now we can combine Eqs. 19.71, 19.72 and 155 to obtain an expression relating the solvent activity to φ p :
ln a s = ln ( 1- φ p ) + φ p + χ φ p 2 + φ p 1/3 N el (19.74)
This equation can be solved numerically to obtain φ p (or its inverse, which gives the swelling ratio V wet / V dry ). Example curves of the swelling ratio obtained in this way are shown in Figure 156.
image: 201_home_ken_Mydocs_MSEcore_331_figures_swelling_curves.svg
Figure 156: Example swelling curves for a fixed value of χ and different values of N el (part a) and for a fixed value of N el and different values of χ (part b).

19.11.2 Shear Deformation of the Swollen Elastomer

Suppose that we now apply a shear strain to the swollen gel. The , γ = λ 1 - λ 2 , with λ 1 λ 2 =1 and λ 3 = λ s (see Section 15.1). For the swollen and sheared state we have λ x = λ s λ 1 , λ y = λ s λ 2 and λ 3 = λ s .
With these extension ratios we have the following expressions for volume ratio the free energy change from state 1 (the isotropically swollen state) to state 2 (the swollen and subsequently sheared state):
Δ f 12 = G dry 2 λ s { λ 1 2 - λ 2 2 +2 } = G dry 2 λ s { λ 1 2 - λ 2 2 } = G dry 2 λ s γ 2 (19.75)
We differentiate twice with respect to γ (see Section 15.4) to obtain the shear modulus of the swollen elastomer:
G wet = d 2 d γ 2 ( Δ f 12 ) = G dry λ s = G dry φ p 1/3 (19.76)
Note that the elastic swelling pressure calculated from the previous section is equal to the shear modulus of the solvent-swollen polymer.

20 Surfaces and Interfaces

20.1 Equilibrium Contact Angle

The shape of a liquid droplet on a solid surface is determined by a balance of forces where the solid, liquid and ambient environment (typically air) come into contact. These three phases meet at the contact line. The surface energy of the liquid ( γ ), the surface energy of the solid ( γ s ) and the interfacial energy of the solid/liquid interface ( γ sl ) are related to one anther by the following expression, obtained by balancing horizontal forces at the contact line.
γ s = γ sl + γ l cos θ e
image: 202_home_ken_Mydocs_MSEcore_331_figures_Contact_angle.png
Figure 157: Diagram of the equilibrium contact angle

20.2 Polymer/polymer Interfacial Free Energy

image: 203_home_ken_Mydocs_MSEcore_331_figures_Interfacial_structure.png
Figure 158: Schematic representation of the interface between immiscible polymers.
The interface between two immiscible polymers is characterized by a finite interfacial width, w . The actual value of w is determined by a balance of enthalpic considerations, which favor a small value of w , and entropic considerations, which favor a large value of w . Below we obtain an expression for the interfacial energy as a function of w , and minimize this with respect to w to get expressions for both w and the interfacial free energy between immiscible A and B polymers.
The enthalpic contribution to the interfacial free energy is obtained by integrating the enthalpy per unit volume, Δ h , associated with contacts between A and B repeat units in the interfacial zone. Here we use the same regular solution form that forms the basis of the Flory-Huggins expression for the free energy of mixing:
Δ h( z ) = RT V 0 χ φ a ( z ) φ b ( z ) (20.1)
If we assume as simple linear expression for the concentration profile in the interfacial zone, we have the following for 0<z<w (as illustrated in Figure 159).
φ a ( z ) =1-z/w; (20.2) φ b ( z ) =z/w (20.3)
image: 204_home_ken_Mydocs_MSEcore_331_figures_Enthalpy.png
Figure 159: Linear approximation to the concentration profile across an interface between immiscible polymers .
These expressions result in the following for γ χ , the enthalpic contribution to the interfacial free energy. We refer to this contribution as γ χ because it includes all of the thermodynamic interactions that are included in χ . In addition to all the enthalpy of mixing, this parameter also includes the non-ideal entropy of mixing:
γ χ =0w Δ h( z ) dz= RT χ V 0 0w( z w ) ( 1- z w ) dz= RTw χ 6 V 0 (20.4)
A simple argument for the entropic contribution every time a random walk goes from one side of the interface to the other, it must turn around, so that we use one half of the available molecular configurations. If the number of molecular configurations is Ω , the configurational entropy of associated with these configuration is k B ln Ω . If the number of configurations decreases by a factor or 2, then the entropy decreases by k B ln ( 2 ) and the free energy increases by k B T ln ( 2 ) . Because our theory is already approximate, we'll assume that this free energy increase per forced turnaround of a polymer molecule is simply k B T .
Our calculation of the entropic contribution to the free energy proceeds in the four steps outlined below:
  1. Determine N w , number of repeat units in random walk with R 0 =w .
    a N w =w; N w = ( w/a ) 2
  2. Calculate v w , the volume of each of these random walks:
    v w = N w v 0 = v 0 ( w/a ) 2
  3. Calculate Γ w , the number of these random walks per unit area in the interfacial region:
    Γ w = w v w = a 2 w v 0
  4. multiply by Γ i by k B T to get γ s , the entropic contribution to the interfacial free energy:
    γ s = k B T a 2 w
The overall interfacial free energy, γ ab , is obtained by adding γ χ and γ s :
γ ab = γ χ + γ s = k B T v 0 ( χ w 6 + a 2 w )
The actual value of w is obtained by minimizing γ with respect to w . We obtain the following:
w=a 6 χ
γ ab = 2a k B T v 0 χ /6 = 2aRT V 0 χ /6
This approximate approach gives the correct scaling ( wa χ -1/2 ; γ ab a k B T v 0 χ 1/2 ), but the exact numerical prefactors are not obtained with quantitative accuracy with this approach. A more complete treatment results in the following for w and γ :
w= 2a 6 χ (20.5)
γ ab = aRT V 0 χ /6 (20.6)

20.3 Block Copolymer Morphologies

Block copolymers consisting of two different polymers characterized by a sufficiently large value of χ form a wide variety of interesting morphologies. In these cases two different polymer types that prefer to separate into different phases are covalently linked to one another along a single molecule. The morphologies observed for simple diblock copolymer molecules (a sequence of A repeat units linked to a sequence of B repeat units) is shown in Figure 160. If the length of the A block is much longer than the length of the B block, then the B blocks segregate into isolated spheres in a continuous matrix of A polymer formed from the A blocks. If the relative length of the B block is increased, the spheres transition to a hexagonal array of cylinders, and then to a morphology where the B cylinders form an interconnected, 3-dimensional network within the A matrix. If the A and B block lengths are similar in length, then the morphology is a simple alternating array of A and B lamellae, as shown schematically in Figure 161. An electron micrograph of a lamellar diblock copolymer morphology is shown in Figure 162.
image: 205_home_ken_Mydocs_MSEcore_331_figures_Diblock_Copolymers.png
Figure 160: Diblock copolymer morphologies.
image: 206_home_ken_Mydocs_MSEcore_331_figures_lamellar_diblock_morphology.png
Figure 161: Molecular orientation in a lamellar diblock copolymer morphology.
image: 207_home_ken_Mydocs_MSEcore_331_figures_block_copolymer_lamellae.png Figure 162: Electron micrograph of a lamellar diblock copolymer morphology.
The free energy per molecule of a lamellar-forming diblock copolymer is the sum of a chain stretching free energy, F s , and a free energy associated with the interface between the A and B domains. We refer to this interfacial free energy per chain as F γ . For our purposes we assume that the A and B blocks have the same statistical segment length, a . The overall volume fractions of the two different blocks, f a and f b , are given by the relative degrees of polymerization of the two blocks:
f a = N a N a + N b = N a N f a = N b N a + N b = N b N (20.7)
Here N is the overall degree of polymerization of the molecule,i.e., N= N a + N b . The overall width of the A region of the block copolymer is f a , and the width of the B region is f b . The stretching free energy is obtained from the requirement that each block of the molecule must have an rms end-to-end distance equal to half the width of the respective region:
R 0 a = f a /2 R 0 b = f b /2 (20.8)
We define f a and If the degrees of polymerization of the two blocks are N a and N b , and the statistical segment length for each of the two different repeating units is a , we have:
F s = 3 2 RT ( R 0 a ) 2 N a a 2 + 3 2 RT ( R 0 b ) 2 N b a 2 = 3 8 2 N a 2 (20.9)
The interfacial free energy is obtained by multiplying the interfacial free energy by the area per molecule along the A/B interface, which we refer to as Σ :
F γ = γ ab Σ (20.10)
This area per molecule is given by dividing the molecular volume, N V 0 , by the length per molecule, /2 :
Σ = 2N V 0 (20.11)
The overall free energy per chain, F tot , is given by the sum of F s and F γ . Combination of Eq. 20.6 for γ ab with Eqs. 20.9, 20.10 and 20.11 gives:
F tot = 3 8 RT 2 N a 2 + 2aNRT χ /6 (20.12)
The equilibrium value of , which we refer to as eq is obtained as the value of for which F tot is minimized. Setting d F tot / d to 0 gives:
eq = ( 8 3 6 ) 1/3 a N 2/3 χ 1/6 =1.03a N 2/3 χ 1/6 (20.13)
We can use the fact that R 0 = N 1/2 a to write the expression for eq in the following way:
eq R 0 =1.03 ( χ N ) 1/6 (20.14)
We can now substitute eq for in Eq. 20.12 to obtain the following expression for the equilibrium free energy per molecule in the lamellar microstructure:
F eq =1.19RT ( χ N ) 1/3 (20.15)
Equation 20.15 gives us the molar free energy for a lamellar phase, but it still does not tell us if this lamellar phase has a lower free energy for the disordered phase, where the A and B blocks remain completely mixed. In order to do this we need to compare to the molar free energy in the disordered phase, where the A and B blocks remain completely mixed with one another. To do this we use the Flory-Huggins free energy of mixing (Eq. 19.28). Because the block copolymer is a one component system, there is no entropy of mixing, so we don't include the ideal entropy of mixing terms. We just have the term involving χ , which accounts for the interactions between the A and B repeat units in the mixed, disordered phase. The free energy per molecule is obtained by multiplying the free energy per unit volume by the molecular volume, N V 0 . We get the following simple expression for F dis , the molar free energy of the disordered phase:
F dis =RT χ φ a φ b =RT χ f a f b (20.16)
For a symmetric diblock copolymer with f a = f b =0.5 , F dis = F eq when χ N is equal to 10.4. When χ N is larger than this value, the ordered phase free energy given by Eq. 20.15 is lower than the disordered phase free energy given by Eq. 20.16. As a result we expect that an ordered lamellar phase will be formed for f a =0.5 when χ N is greater than 10.4. A disordered phase will be formed for χ N<10.4 . This value of 10.5 is remarkably similar to the value of 10.5 that is obtained by a much more detailed calculation. These more detailed calculations are necessary in order to understand how the preferred geometry of the ordered phase (spheres, cylinders, etc.) depends on the f a . A full calculation of the phase diagram is shown in Figure 163.
image: 208_home_ken_Mydocs_MSEcore_331_figures_bcp_phase_diagram.png Figure 163: Block copolymer phase diagram (adapted from ref.[6]).

Case Study: Thermoreversible Gels from Triblock Copolymer Solutions

21.1 Introduction

Many of the principles of this course are illustrated by the behavior solutions of acrylic triblock copolymers in higher alcohols. The specific molecular structures of the polymer and the solvent are shown in Figure 164 and consist of a midblock of poly(n-butyl acrylate) (PnBA) between two blocks of poly(methyl methacrylate) (PMMA). The value of χ for the PnBA/solvent system is less than 0.5 for all temperatures of interest, but the χ for the PMMA/solvent system has the following temperature dependence:
χ =1.45-0.0115T (21.1)
with T in C.
The following points are illustrated by this example:
  1. Polymers are soluble in solvents when χ <0.5.
  2. For χ >0.5, the solvent content within the polymer can be obtained from the chemical potential expression.
  3. To understand the temperature dependence of the solubility you need to know the temperature dependence of χ .
  4. Polymer liquids have higher heat capacities than the corresponding polymer glasses.
  5. Calorimetry can be used to detect the enthalpy recovery peak from an aged polymer glass.
  6. Gelation occurs when the average functionality is 2.
  7. Time-temperature superposition works when the structure does not change appreciably with Temperature.
  8. The viscosity is obtained by integrating G( t ) or from the response at very low frequencies.
(a)
Triblock Copolymer
image: 209_home_ken_Mydocs_MSEcore_331_figures_triblock_comp.png
(b)
Solvents
image: 210_home_ken_Mydocs_MSEcore_331_figures_butanol.png
(butanol)
image: 211_home_ken_Mydocs_MSEcore_331_figures_ethylhexanol.png
(ethyl hexanol)
Figure 164: Chemical structure of gel-forming triblock copolymer gels (a), and the alcohols used as the solvents (b).
From Eq. 21.1 we see that the value of χ characterizing the PMMA/solvent interaction is greater than 0.5 for T<83 C. At temperatures above 83 C the PMMA and PnBA blocks are both in good solvent conditions, and the polymer behaves as a normal polymer solution. At lower temperatures, however, the PnBA midblock remains in good solvent conditions, but the PMMA endblocks are no longer soluble. As a result these endblocks aggregate to form the micellar structures illustrated schematically in Figure 165. As a solution is cooled, χ increases and solvent is expelled from these aggregates. When sufficient solvent has been expelled from the PMMA aggregates, they become glassy, and the system behaves as a solid.
image: 212_home_ken_Mydocs_MSEcore_331_figures_triblock_regimes.png Figure 165: Schematic representation of the temperature and concentration dependence of the structure of the triblock copolymer solutions.

21.2 Thermoreversible Gelcasting of Ceramics

The transition from 'liquid-like' to 'solid-like' behavior as the temperature is changes enable these materials to be used in a variety of applications, including the ceramics processing application illustrated in Figures166 and167. In this application ceramic particles are added to the triblock copolymer at a relatively high volume fraction (around 50% by volume). The thermoreversible nature of the triblock copolymer solution allows the suspension to be poured into a silicone mold at high tempertures (typically 80 C). When the system is cooled to room temperature, the triblock copolymer forms a solid gel and the material can be removed from the mold. A ceramic material with the shape of the mold is obtained by letting the solvent evaporate from the suspension, and firing the object at elevated temperature to burn off the polymer and sinter the ceramic particles together.
image: 213_home_ken_Mydocs_MSEcore_331_figures_alumina_dispersion_schematic.png Figure 166: Schematic representation of the alumina particles dispersed within a thermoreversible triblock copolymer gel.
image: 214_home_ken_Mydocs_MSEcore_331_figures_gelcasting_schematic.png Figure 167: Illustration of the thermoreversible gelcasting process.

21.3 Quantifying the Solid/Liquid Transition

The temperature dependence of the mechanical properties of the triblock copolymer solutions are quantified by the frequency dependence of the storage and loss moduli, G ( ω ) and G ( ω ) . It is conceptually simpler to think in terms of the time-dependence of the relaxation modulus, G( t ) , which can be written as a sum of exponential relaxations:
G( t ) = i G i exp ( -t/ τ i ) (21.2)
In the frequency domain, the values storage and loss moduli are given by the following expressions:
G ( ω ) = i G i τ i 2 ω 2 1+ τ i 2 ω 2 (21.3)
G ( ω ) = i G i τ i ω 1+ τ i 2 ω 2 (21.4)
The most important thing for us is the viscosity, which we can obtain in terms of the values of G i and τ i , or in terms of the loss modulus at very low frequencies:
η = i G i τ i = lim ω 0 G ω (21.5)
We generally obtain data as a master plot (where the frequency is multiplied by the temperature shift factor, a T ). In this case we obtain the temperature dependence of the viscosity by multiplying by a T :
η = a T lim a T ω 0 G a T ω (21.6)
We can also obtain an expression for the limiting modulus at high frequency (or low temperature):
G 0 = i G i (21.7)
Finally, we define an average relaxation time, τ av in the following way:
τ av = η G 0 = G i τ i G i (21.8)
Values of the storage and loss moduli for the triblock copolymer solutions are shown in Figure 168, along with the values of G i and τ i used to fit the data. The temperature dependent viscosity for a solution with φ p =0.15 is shown in Figure 169, and the temperature dependent relaxation times for a variety of solution concentrations are shown in Figure 170.
image: 215_home_ken_Mydocs_MSEcore_331_figures_gel_rheology.png Figure 168: a) Master curves for G and G (symbols) with fits to Eqs.21.3 and21.4. b) Values of G i and τ i used to fit the data from part a (the relaxation times, τ i are referred to by λ i in the figure). The reference temperature (defined as the temperature where a T =1 ) is 65 C for these data.
image: 216_home_ken_Mydocs_MSEcore_331_figures_triblock_viscosity.png Figure 169: Temperature dependence of the viscosity for a triblock copolymer solution with φ p = 0.15. in 2-ethyl hexanol. The x's represent the solvent viscosity.
image: 217_home_ken_Mydocs_MSEcore_331_figures_Relaxation_times_triblock.png Figure 170: Values of the average relaxation time, τ av as a function of temperature. Data are included for φ p =0.05( ) φ p =0.10( ) φ p =0.15( ) φ p =0.20( ) φ p =0.25( ) φ p =0.30( ) .

21.4 Characterizing the Glass Transition in PMMA domains

The glass transition is typically observed as step change in heat flow as a sample is cooled through the glass transition. In our case the PMMA domains go through a glass transition during cooling, but the PMMA concentration is so low that the glass transition is not clearly observed during cooling or when when the sample is immediately reheated after cooling (see Figure 171). However, a large enthalpy peak is observed after the sample has been left to 'age' at room temperature for a long time before reheating the sample. As shown in Figure 172, the enthalpy peak gets larger and moves to higher temperature as the aging time at room temperature is increased. This peak has the appearance of a melting peak associated with crystallization, but in this case it arises from the glassy PMMA domains.
The origins of the enthalpy peak observed during the aging experiments can be understood by realizing that the glass transition separates an equilibrium, liquid regime at temperatures above T g from a non-equilibrium, glassy regime at temperatures below T g . At temperatures below the glass transition, the enthalpy content of the sample is higher than the equilibrium enthalpy content defined by the extrapolation of the liquid behavior (the dashed line in Figure 173). As a result, the enthalpy slowly decreases toward this equilibrium line as the sample is aged below the glass transition. The decrease in enthalpy during this aging process is labeled as Δ H a in Figure 173. When the sample is reheated, the enthalpy increases with temperature according to the glassy heat capacity, C p g . As a result the enthalpy content eventually crosses the equilibrium line and becomes less than the equilibrium enthalpy content. At a temperature somewhere just above T g , the sample is able to equilibrate, and the enthalpy increases by an amount Δ H r in order to catch up to the equilibrium value. This is the enthalpy corresponding to the peaks in Figure 172. As the aging time decreases, Δ H a and Δ H r both decrease, and the temperature at which the enthalpy is recovered moves closer to T g . The actual value of the T g can be estimated by measuring the area of the enthalpy recovery peak and plotting against the location of this peak. Doing this for the data shown in Figure 172 results in an estimate for T g of the PMMA domains of 35 C.
image: 218_home_ken_Mydocs_MSEcore_331_figures_DSC_heat_flow_for_triblocks.png Figure 171: Measured heat flow during cooling, immediate reheating, and reheating after room temperature storage for 130 days.
image: 219_home_ken_Mydocs_MSEcore_331_figures_enthalpy_peak_aging.png
Figure 172: Dependence on the enthalpy recovery peak on the aging time.
image: 220_home_ken_Mydocs_MSEcore_331_figures_enthalpy_recovery.png
Figure 173: Origins of the enthalpy recovery peak.
As a final check on all of this, we can see if it makes sense that the glass transition of the PMMA domains should be around 35 C. The glass transition of PMMA is about 125 C, which is certainly a lot higher than the value of 35 C that we are claiming here. However, we need to remember that the PMMA domains contain a lot of solvent, and this solvent will depress T g significantly. How much solvent do we expect there to be in the PMMA domains? From Eq. 21.1 we obtain χ =1.05 at T=35 C. The relationship between χ and φ s , the solvent volume fraction in the PMMA cores of the micelles is obtained by setting N p to and setting the solvent chemical potential (from Eq19.43) to 0:
ln ( φ s ) +1- φ s + χ ( 1- φ s ) 2 =0 ln ( 1- φ p ) + φ p + χ φ p 2 =0 (21.9)
This equation needs to be solved numerically to obtain φ s (or φ p ) from a specified value of χ . A simple MATLAB script used to do this is given in Section25.2. A value of 1.05 for χ gives φ s =0.289. Independent measures of the glass transition temperature of bulk PMMA samples with solvent loadings comparable to this do in fact show measured glass transition temperatures close to 35 C. The relationship between χ and φ p given by Eq. 21.9 is plotted in Figure 174.
image: 221_home_ken_Mydocs_MSEcore_331_figures_polymerswelling.png
Figure 174: Relationship between φ p and χ obtained from Eq. 21.9.

21.5 Concentration Dependence of the Gel Modulus

Elasticity of the triblock copolymer solutions at low temperatures arises from the fact that the middle, PnBA blocks of the copolymer can span different PMMA aggregates, thereby linking the whole structure together. The PMMA aggregates behave as physical crosslinks with a functionality given by the number of midblocks that that bridge different aggregates. This functionality is obtained by multiplying the aggregation number,i.e. the number of PMMA endblocks in a single aggregate, by the probability f that a PnBA midlbock spans two different PMMA aggregates (we also have to divide by two to account for the fact that there are two PMMA blocks on each triblock copolymer molecule). Aggregation numbers depend on the polymer concentration, and are typically very large, as illustrated in Figure 175. For high molecular weight polymers the percolation threshold, where the average functionality of a micelle is 2, is quite low, corresponding to φ p 0.035. Above this percolation threshold, the shear modulus is given by the following expression:
G 0 = ν k B T D 2 R 0 2 (21.10)
Here D is the average distance between micelle cores, which can be measured experimentally by x-ray scattering. We also have R 0 = N 1/2 a , where N is the midblock degree of polymerization and a is the statistical segment length for the midblock. The quantity ν is the conentration of 'load bearing strands', which in this case is the concentration of triblock copolymer chains with bridging midblocks:
ν = f φ p ρ N av M
where M is the molecular weight of the triblock copolymer molecule and ρ is the polymer density. The modulus is strongly concentration dependent because of the concentration dependence of the bridging fraction, f , which is shown in Figure 176. Measured and calculated values of G 0 are shown in Figure 177.
image: 222_home_ken_Mydocs_MSEcore_331_figures_aggregation_number.png Figure 175: Concentration dependence of aggregation number for three different PMMA-PnBA-PMMA block copolymer solutions. The three data sets correspond to three differnt triblock copolymers, with different molecular weights for the PMMA and PnBA blocks.
image: 223_home_ken_Mydocs_MSEcore_331_figures_structrual_gelation.png Figure 176: Concentration dependence of the network functionality.
image: 224_home_ken_Mydocs_MSEcore_331_figures_conc_dep_of_G_for_gels.png Figure 177: Calculated and experimental concentration dependence of the network functionality.

21.6 Hydrogels: Water as the Solvent

As a final illustration of what can be done with these sorts of triblock copolymer gels, we consider materials where the midblock is replaced with poly(acrylic acid), a polymer that is water soluble at neutral pH. The structure of these polymer is shown in Figure 178. This figure also shows a scheme for forming gels from these materials. Instead of adjusting χ between the solvent and the PMMA endblocks by changing temperature, we do this by adding a small amount of water to the solvent (which is initially dimethyl sulfoxide). Addition of just a small amount of water increases the effective value of χ characterizing the solvent/PMMA interaction. The result is that the relaxation times for the triblock copolymer solution increase dramatically, as illustrated by the rheological data in Figure 179. The effects of solvent composition can be illustrated by introducing a shift factor, a s that depends on the composition of the solvent. It's use is illustrated in Figure 180. The viscosity at a given temperature is proportional to a s , so we see that small increases in the water content of the solvent result in an increase in the solution viscosity by several orders of magnitude. This occurs because water induces the aggregation of PMMA blocks in to discrete domains, just as reducing temperature did for the case where alcohol was used as the solvent.
image: 225_home_ken_Mydocs_MSEcore_331_figures_ionic_gels.png Figure 178: Formation and ionic crosslinking of PMMA-PMAA-PMMA triblock copolymer hydrogels (PMAA = poly(methacrylic acid).
image: 226_home_ken_Mydocs_MSEcore_331_figures_dmso_water_curves.png
Figure 179: Master plots of storage and loss moduli for PMMA-PMAA-PMMA triblock copolymer gels in DMSO/water mixtures, for water contents within the solvent of 4, 6 and 8 wt. %.
image: 227_home_ken_Mydocs_MSEcore_331_figures_dmso_water_master.png Figure 180: 'Super Master' plot of the data from Figure.
A useful feature of the hydrogels formed from the triblock copolymers with the poly(acrylic acid midblocks) is that the properties can be changed dramatically by immersing the gels in solutions containing divalent cations like Ca 2+ or Zn 2+ . These ions add additional crosslinks to the material, increasing the modulus by a factor of 100 or more (Figure 181). The materials also have high toughness, and can be extended to several times their original length prior to failure (Figure 182).
image: 228_home_ken_Mydocs_MSEcore_331_figures_ionic_crosslinking.png Figure 181: Elastic modulus for triblock copolymer gels that have been ionically crosslinked to varying degrees.
image: 229_home_ken_Mydocs_MSEcore_331_figures_ionic_gels_in_tension.png
Figure 182: Tensile stress/strain curves for ionically crosslinked triblock copolymer gels.

Summary

22.1 Classification Scheme

image: 26_home_ken_Mydocs_MSEcore_331_figures_classification_scheme.png
Figure 183: Summary of polymer properties and classes

22.2 Molecular Weight Averages

N: Degree of polymerization
M: Molecular weight
M 0 Molecular weight per repeat unit
M n =w/n
M w = 1 w i w i M i = Weight Average Molecular
w = total weight of polymer
n = total number of polymer molecules
w i = weight of fraction i (grams)
M i = molecular weight of fraction i

22.3 Polymerizations

22.3.1 Linear step growth polymerization

M n = M 0 2( 1-p )
M w M n =1+p
n= n a - n b
(p = extent of reaction)

22.3.2 Nonlinear step growth polymerization

M n = M 0 2-p f av
( n a + n b ) f av = n a f a + n b f b
f a , f b = functionality of A and B monomers
n a , n b = numbers of A and B monomers
p f av = 2 at the gel point

22.3.3 Chain Growth Polymerization

Active site is typically a free radical, anion, or cation.
Propagation corresponds to addition to a double bond, or opening of a ring.
Termination can occur by Disproportionation or combination.
Intermolecular chain transfer reduces molecular weight.
Intramolecular chain transfer results in chain branching.

22.4 Chain Dimensions

22.4.1 Structure of Amorphous Polymers

The probability density distribution for the end to end vector, R (with components R x , R y and R z ) is:
Here a is the statistical segment length, defined by the following expression:

22.4.2 Structure of amorphous vinyl polymers

R 0 2 = C N b l 2 =N a 2
C ={ 1+ cos θ 1- cos θ } { 1+< cos φ > 1-< cos φ > }
θ = bond angle
φ = rotation angle

22.4.3 Glass Transition

V f V =( α l - α g ) ( T g - T ) (Free volume theory)

22.4.4 Differential Scanning Calorimetry

Thermal signatures of melting, crystallization and glass formation

22.5 Crystallization

Polymer Crystallization
image: 230_home_ken_Mydocs_MSEcore_331_figures_Summary_2.png
Figure 184: Polymer crystallization
λ 2 γ e T m 0 Δ h Δ T
% crystallinity = /d
x*y/z notation for helix formation
spherulitic morphologies

22.6 Rubber Elasticity

Single Chain:
Δ F d = 3 k B T R 2 2N a 2
Bulk Elastomer:
Δ f d = v s k B T 2 { λ x 2 + λ y 2 + λ z 2 -3 }
Shear: σ xy =G γ ; G= v s k B T= ρ RT M
Uniaxial Extension or Compression: σ eng =G{ λ y - 1 λ y 2 }
Δ f d = deformation free energy (per volume)
σ eng = engineering stress
σ xy = shear stress
γ = shear strain
v s = strand concentration
λ x , λ y , λ z = extension ratios

22.7 Viscoelasticity

Response to step strain:
G( t ) σ xy ( t ) γ ( 0 )
Boltzmann superposition principle:
σ xy ( t ) = 0G( s ) γ ( t-s ) ds
image: 231_home_ken_Mydocs_MSEcore_331_figures_Summary_10.png
Figure 185: Bolzmann superposition principle on a graph
γ ˙ = strain rate = d γ dt
Dynamic Mechanical Spectroscopy
γ = γ 0 sin ( ω t )
γ ( t ) = ω γ 0 cos ( ω t )
σ 12 ( t ) γ 0 =G'( w ) sin ( ω t ) +G"( w ) cos ( ω t )
tan δ = G"( ω ) G'( ω )
G = storage modulus (elastic response)
G = loss modulus (viscous response)
δ = phase angle
η 0 = lim ω 0 G ( ω ) ω =0G( t ) dt=A exp ( B T- T )
Time-temperature superposition: Changes in temperature shift the G(t), G'( ω ) and G"( ω ) curves horizontally (when plotted against log(t) or log( ω )), but the shapes of the curves do not

22.8 Solutions and Blends

Free energy of mixing
Δ f mix = RT V 0 { φ a ln φ a N a + φ b ln φ N b + χ φ a φ b }
π =- μ s v 0
Chemical potential:
μ a RT = ln φ a + φ b ( 1- N a N b ) + χ N a φ b 2
π = RTc M n in dilute limit
Spinodal curve: χ s = 1 2 φ a N a + 1 2( 1- φ a ) N b
Critical point: φ a,crit = N b N b + N a
χ crit = N a + N b +2 N a N b 2 N a N b ={ 0.5 for N a =1, N b 2 N for N a = N b =N .

22.9 Interfaces

Polymer/Polymer Interfaces
w= 2a 6 χ
γ = a V 0 RT χ /6

Review Questions

23.1 Polymer Synthesis

23.2 Molecular Weight Distributions

23.3 Gels and Elastomers

23.4 Glass transition

23.5 Chain dimensions of amorphous polymers

23.6 Rubber Elasticity

23.7 Viscoelasticity

23.8 Semicrystalline Polymers

23.9 Solution Thermodynamics

23.10 Interfaces and Self-Assembly

331 Problems

Basic chemistry review

Some basic problems having to do with the chemical structures of polymers.
Calculate the repeat unit molecular weights for each of the 22 polymers shown in the 'common polymers' section of the text.
Draw Lewis diagrams illustrating the valence shell configurations for polystyrene, poly(methyl methacrylate) and poly(dimethylsiloxane).
Draw Lewis diagrams illustrating the valence shell configurations for amide, ester and urethane linkages.

Polymerization Reactions

Suppose a polyester is formed by combining 1.00 mol of a diacid (2 acid groups per molecule) with 1.01 mol of a dialcohol in a homogeneous solution. The molecular weight of the diacid is 160g/mol and the molecular weight of the dialcohol is 180 g/mol.
  1. What is the highest value of M n that can be expected from this combination?
  1. How many grams of water will be produced during the polymerization reaction?
Consider the following 5 monomers:
image: 232_home_ken_Mydocs_MSEcore_331_figures_2-1-problem.png
  1. Draw the repeat units for three linear polymers that can be produced by reactions between the monomers in this list. Identify these polymers according to their type (polyamides, polyesters, etc.), and indicate whether each polymerization reaction is condensation reaction or not.
  2. What combination of monomers from this list would you choose in order to produce a three-dimensional network?
  3. Draw the repeat unit structure of a linear polyamide formed by a condensation polymerization involving a combination of the molecules from this list, and identify the condensation product.
Suppose a diamine (primary amines at each of the molecule) is reacted with a di-epoxide (epoxide groups at each end of the molecule).
  1. What ratio of diamine molecules to di-epoxide molecules would you use in order to optimize the ability of the system to gel (i.e., form a network)?
  2. Estimate the fraction of amine groups will need to react in order for a gel to form? Identify any assumptions that you have made in arriving at this number.
Why is tacticity relevant for polystyrene or polypropylene, but not polyethylene?
Describe methods that could be used to make the following two block copolymers:
  1. A polystyrene/polyisoprene diblock polymer
  2. A mutiblock copolymer consisting of alternating Nylon 6-6 and Nylon 6-10 blocks.
  3. Comment on the polydispersity of block lengths, and of the number of blocks within a molecule, which can obtained by the methods described in parts a and b.
A step-growth reactions takes place between 1.2 moles of dicarboxylic acid, 0.4 moles of glycerol (a trifunctional alcohol) and 0.6 moles of ethylene glycol (a difunctional alcohol). Calculate the critical extent of reaction for gelation using the Carothers theory. Describe why you would expect gelation to occur for values of p that are different from this predicted value. (Hint: you will need to generalize the expression for the average functionality, f av , to include the possibility of more than two different monomers.)
In 1963 Karl Ziegler and Giulio Natta received the Nobel prize in chemistry for the development of processes which could be used to produce (among other things) linear, unbranched isotactic polypropylene. Why was this an important development? What useful properties does isotactic polypropylene have when compared to atactic polypropylene? Describe some applications where isotactic polypropylene could be used, but where it would be impractical to use atactic polypropylene. (Note: the glass transition temperature of both forms of polypropylene is well below room temperature.)
“Living” polymerizations are chain growth polymerizations that proceed without termination or chain transfer reactions.
  1. Discuss the conditions that must be met in order for polymers with a low polydispersity to be formed by a living polymerization.
  2. What type of reactions must be avoided in order to produce block copolymers by radical chain growth polymerization.
  3. Proteins are perfectly disperse copolymers, with a polydispersity index of exactly 1. While relatively monodisperse polymers ( M w / M n <1.05 ) can be made by living polymerization methods like anionic polymerization, it is not possible to make perfectly monodisperse polymer by these methods, even when termination and chain transfer events are completely eliminated. Why is that?
Suppose 100 moles of styrene monomer are polymerized by anionic polymerization, using 0.1 mole of initiator.
  1. Draw the repeat unit structure of the polymer which is formed in this reaction.
  2. Below what temperature will this polymer behave as a rigid material?
  3. Suppose that there are no chain transfer or termination events during the resulting polymerization. Calculate the number average molecular weight ( M n ) of the polymer that is produced.
Suppose nylon is formed in two separate batches by reacting the same sets of monomers with one another. In the first reaction, an extent of reaction of 0.999 is obtained. In the second reaction, an extent of reaction of 0.9999 is obtained. Assume that an ideal distribution of molecular weights is obtained for each reaction.
  1. What is the ratio of M w for the second reaction to M w for the first reaction?
  2. Suppose equal weights of the polymers formed by these two reactions are mixed with one another. What is the polydispersity index of the resultant mixture?
Draw the chemical structures of the monomers from which the following polymers are formed.
image: 233_home_ken_Mydocs_MSEcore_331_figures_9-1-problem.png
Suppose you want to sell a set of cheap plastic mugs which are suitable for drinking coffee. Briefly discuss the potential applicability of the following materials for this purpose:
  1. atactic polystyrene
  2. atactic poly(vinyl chloride)
  3. high density polyethylene
Consider the the piece of the highly crosslinked chemical structure shown below.
  1. Draw the chemical structures of the two different monomers that were used to form this material. (Note that the R groups just represent connections to the rest of the material, they are not part of the structures you need to draw).
    image: 234_home_ken_Mydocs_MSEcore_331_figures_rolayer.png
  1. What type of polymer is this (polyester, polyamide, polyurethane, epoxy, etc.)
  2. Next to each structure that you have drawn in part a, indicate its chemical functionality.
  3. What is the relative amount of each of the two monomers from part a (in moles) that you should be added to one another to optimize the polymerization?
  4. Assuming that the two monomers are added together in optimum proportions you mentioned in part d, what fraction of reactive groups need to react to form an insoluble gel.

The Glass Transition

An amorphous polymer has T =5 0 C. The thermal expansion coefficient in the glassy state is 3 × 1 0 -4 / C , and the thermal expansion coefficient in the liquid state is 4.5 × 1 0 -4 / C.
  1. Calculate the glass transition temperature if the fractional free volume at the glass transition is 0.025.
  1. The glass transition temperature measured at a different cooling rate is found to correspond to a fractional free volume of 0.021. Is this cooling rate higher of lower than the cooling rate from part a?
  2. Rank these three polymers in the order of high T g to low T g , and provide an explanation for your prediction:
    Polystyrene
    image: 84_home_ken_Mydocs_MSEcore_331_figures_polystyrene.png
    Poly(para butylstyrene)
    image: 235_home_ken_Mydocs_MSEcore_331_figures_polybutylstyrene.png
    Poly( α -methylstyrene)
    image: 236_home_ken_Mydocs_MSEcore_331_figures_polymethylstyrene.png
The DSC curves below correspond to the 3 polymers listed below. Identify which is which, and identify on the curves any crystallization (T c ), glass transition (T g ) or melting (T m ) temperatures, and give your reasoning in each case.
image: 237_home_ken_Mydocs_MSEcore_331_figures_dsc_slide1.png
  1. Atactic polystyrene being heated from room temperature to 200 o C.
  2. High density polyethylene being cooled from 200 o C to room temperature.
  3. Poly(ethylene terephthalate) that had originally been cooled rapidly from 300 o C to room temperature, and is now being slowly heated from room temperature back up to 300 o C.
Short, bulky side groups (like methyl of phenyl groups) increase the fractional free volume required for the motion of polymer chain segments, whereas long, thin groups (like the n-butyl group - C H 2 C H 2 C H 2 C H 3 ) increase the available free volume.
  1. Describe how each of these two types of groups is expected to affect the glass transition temperature.
  2. Polymers with very little flexibility along the backbone of the polymer chain have very high glass transition temperatures. Give a brief explanation for this in terms of free volume theory.

Chain Dimensions

The value of R 0 , the rms (root mean squared) end-to-end distance for a polystyrene molecule in the bulk material (no solvent) is 200Å. The molecular weight, M, of this molecule is 100, 000 g/mol.
  1. What is the rms end-to-end distance for a polystyrene molecule with M=250,000 g/mol?
  2. What is the characteristic ratio ( C ) for polystyrene?
  3. What is the statistical segment length (equivalent length of a repeat unit in a random walk model of chain dimensions) for polystyrene?
  4. The volume pervaded by a single molecule is given approximately by R 0 3 . For a polystyrene molecule with M=250, 000 g/mol, compare this volume to the volume actually occupied by the polystyrene molecule. How many different molecules overlap with one another in polystyrene with this molecular weight?
  5. Consider a simple model of polystyrene for which all of the bonds are in either gauche or trans configurations. What fraction of bonds are in the trans configuration?
Consider a simple model for the chain dimensions of amorphous poly(vinyl chloride) where only trans and gauche bonds are possible.
Poly(vinyl chloride)
Trans Configuration
Gauche Configuration
image: 238_home_ken_Mydocs_MSEcore_331_figures_17-1-PVC.png
image: 239_home_ken_Mydocs_MSEcore_331_figures_17-1-trans.png
image: 240_home_ken_Mydocs_MSEcore_331_figures_17-1-gauche.png
  1. What fraction of the backbone bonds are in the gauche configuration if C for the polymer is 8?
  2. Calculate the statistical segment length of a repeat unit, assuming the value of C from part a.
What fraction of molecules in an equilibrated amorphous polymer melt have end-to-end vectors with magnitudes larger than the root-mean-square end-to-end distance? (note – this will require a numerical or graphical solution.)

Mechanical Properties

A thermoset is cooled to a very low temperature. As the temperature is increased, the elastic modulus remains roughly constant until T g , at which point the elastic modulus rapidly decreases. Above the glass transition temperature the elastic modulus slowly increases. Describe in qualitative terms what is going on at the glass transition, and why the modulus actually increases slightly above T g .
Values of Young’s moduli ( E ) for elastomers are typically near 1 0 6 Pa (1 MPa).
  1. What is the average molecular weight between crosslinks for a polymer with this value of E ? Assume the polymer has a specific gravity of 1. What fraction of cis 1-4 polyisoprene repeat units contain crosslinks if the modulus is 1 0 6 Pa?
  2. Plot the predicted tensile stress/strain curve for this material for engineering strains from -0.5 to 2 (extension ratios from 0.5 to 3).
The yield point of a material is often defined as the onset of nonlinear behavior in the stress-strain curve. Can this definition be used for an elastomeric material that is deformed in tension? Why or why not?
Suppose a sample of Sylgard 184 is cured, and the degree of polymerization of the prepolymer is 60. Estimate the value for both the shear modulus and Young's modulus of this polymer. Describe any assumptions that you've made.
Consider the deformation behavior of a uniaxially stretched rubber.
  1. Plot the free energy of a uniaxially stretched rubber as a function of the extension ratio from λ = 0.1 (compressive strain to λ =4 (tensile strain). (The extension ratio of interest here corresponds to the direction in which the rubber is stretched.)
  2. At what value of the strain is the free energy minimized?
  3. For what compressive value of the extension ratio is the free energy of the rubber equal to the free energy it has for an extension ratio equal to 4?
Calculate the force required to fix the ends of a 1,4 polyisoprene molecule with M = 100,000 g/mol at the separation of 1000Å. Assume C = 8, and that the double and single bonds have equal lengths of 1.54 Å.
Consider a material which is defined between two circular plates that are parallel to one another. Each plate has a diameter of 2.5 cm, and the plates are separated by 1mm.
  1. Calculate the maximum shear strain at the outer periphery of the disks, if the top disk is rotated by 1 with respect to the bottom disk.
  2. Calculate the torque required to rotate the top plate by 1 , in terms of the magnitude of the complex shear modulus, G* . (You will need to integrate contributions from all portions of the disk, keeping in mind that the local stress is equal to the shear modulus multiplied by the local strain.)
  3. Suppose the material between the two plates is a purely viscous fluid ( G'=0 at the appropriate frequency), and that the viscosity of the fluid is 1000 Pa-sec. Calculate the magnitude of the torque which must be applied to give a sinusoidal rotation with a magnitude of 1 . Assume that the angular frequency of the applied torque is 100 radians/second.
Consider a polymer with a density of 1 g/c m 3 that has the following relaxation modulus at room temperature (the molecular weight is much higher than the average molecular weight between entanglements).
G( t ) = G 1 exp ( -t τ 1 ) + G 2 exp ( -t τ 2 )
Suppose G 1 =1 0 9 Pa , G 2 = 3x1 0 5 Pa , τ 1 = 1 0 -6 s and τ 2 = 1 s.
  1. Calculate the zero shear viscosity for this polymer.
  2. Is this polymer crosslinked? How do you know?
  3. Calculate the average molecular weight between entanglements.
  4. Describe what happens to τ 1 and τ 2 when the temperature is increased (do they remain constant, increase or decrease)?
  5. What vibrational frequencies do you expect to be most efficiently damped by this polymer?
  6. Describe what happens to τ 1 and τ 2 when the molecular weight of the polymer is increased (do they remain constant, increase or decrease)?
Consider a polymer with the behavior of the relaxation modulus ( G( t ) ) and temperature shift factor ( a T ) plotted below. Assume the polymer has a density of 1 g/cm 3 .
image: 241_home_ken_Mydocs_MSEcore_331_figures_gtplot.png
image: 242_home_ken_Mydocs_MSEcore_331_figures_atplot.png
  1. Describe in qualitative terms what is going on at each of the two characteristic values of t/ a T where G( t ) decreases substantially.
  2. Is this polymer crosslinked? Briefly explain your answer.
  3. Estimate the entanglement molecular weight for this polymer.
  4. Estimate the viscosity of this polymer that is obtained at very low shear rates, at a temperature of 40 C .
  5. What is the ratio of viscosity of the polymer at -40 C to the viscosity at 40 C?
  6. Suppose that G ( ω ) is measured at a very high frequency of 1 0 9 s -1 . What value of G do you expect to obtain from this experiment?
Use the PtBA data in the appendix to address the following:
  1. Determine the shift factors as a function of temperature for poly(t-butyl acrylate) (PtBA), using 75 C as the reference temperature. Plot these shift factors as a function of temperature, and determine a value for T by fitting to the Vogel equation. Compare this value of T to T g for PtBA. Values of G and G at different temperatures are included in the data tables at the end of the text.
  2. Determine the zero shear viscosity of PtBA with M=350,000 g/mol at 7 5 C .
  3. Determine the entanglement molecular weight, M e , for PTBA. Assume a density of 1 g/cm 3 for PtBA.
Using the appendix data for poly(2-vinylpyridine), repeat the previous 3 sub-problems for poly(2-vinylpyridine) with a molecular weight of 82,000 g/mol, using 16 0 C as the reference temperature (where a T =1) , and 13 0 C for the viscosity calculation.
An elastomer is made from a poly(n-butyl acrylate) precursor molecule with double bonds at either end that are able to react with a tetrafunctional crosslinker. The chemical structure is shown at the left below, and a schematic representation of the crosslinked structure is shown at the right.
image: 243_home_ken_Mydocs_MSEcore_331_figures_pnba.png image: 244_home_ken_Mydocs_MSEcore_331_figures_crosslinked_structure.png
  1. Calculate the shear modulus, G, for this material at room temperature.
  2. Give Young's modulus (E) and the Poisson's ratio ( ν ) for this material.
  3. What is the root mean square end-to-end distance R 2 1/2 for crosslink points at either end of the one of the precursor chains when the material is in its undeformed state?
  4. What is the value of C for poly(n-butyl acrylate)?
  5. Express the deformation free energies of the elastomer, divided by the shear modulus, G , for the following strain states:
    1. Simple shear at a shear strain, γ , of 1.
    2. Uniaxial extension to an engineering strain of 100 % ( λ =2) .
    3. Uniaxial compression to a engineering strain of -50% ( λ =0.5) .
For each if the following pairs of quantities, circle the largest one. If you expect them to have the same value, write 'same'.
quantity 1
quantity 2
The thermal expansion coefficient of poly(n-butyl acrylate) at temperatures above -45 C.
The thermal expansion coefficient of poly(n-butyl acrylate) at temperatures below -45 C.
the free volume for a poly(n-butyl acrylate) sample above -45 C.
the free volume for a poly(n-butyl acrylate) sample below -45 C.
The free energy of an undeformed poly(n-butyl acrylate) sample.
The free energy of poly(n-butyl acrylate) sample at a shear strain of 0.3
The entropy of an undeformed poly(n-butyl acrylate) sample.
The entropy of poly(n-butyl acrylate) sample at a shear strain of 0.3
T g for poly(n-butyl acrylate)
T g for poly(methyl acrylate) (where the 4-carbon butyl group in poly(n-butyl acrylate) is replaced by a methyl group.

Crystallinity

Consider a polyethylene molecule with a molecular weight of 150,000 g/mol.
  1. Estimate the number of times the molecule folds back on itself within a lamellar crystallite with a thickness of 100 Å. Assume for the purposes of this calculation that the molecule remains in the crystallite, and does not move into the amorphous regions of the polymer.
  2. From your answer to part a, what can you say about the ability to stretch a polymer molecule in a semicrystalline polymer without breaking any bonds?
Describe how axis of the unit cell in polyethylene is oriented with respect to a spherulite, with respect to an individual lamellar crystallite, and with respect to the molecular chain axis.
Consider a polymer that crystallizes with an orthorhombic unit cell (all angle =9 0 ). The polymer has a carbon backbone, and crystallizes with a 4*1/1 helical structure. The unit cell dimension are a=10.4 , b=5.3Å with 1.27Å per backbone carbon along the c axis. The molecular weight of a polymeric repeat unit is 125 g/mol, and there are two of these units per unit cell.
  1. What is the theoretical density of a completely crystalline version of this polymer (100% crystalline)?
  2. The percent crystallinity of a sample with a density of 1.35 g/cc is 60%. What is the density of the amorphous regions of this polymer?
Poly(ethylene terephthalate) has a relatively high equilibrium melting point (265 °C), but the polymer at room temperature is almost completely amorphous. Based on this result, do you expect that T g for PET is above room temperature or below room temperature? Explain your reasoning.
Consider two perfect single crystals of polyethylene and polypropylene. Polypropylene forms a 3/1 helix, and polyethylene forms a 2/1 helix. Which of these single crystals would you expect to have a higher modulus along the C axis? Include in your discussion the relative importance of bond bending and bond rotation during elastic deformation of the two materials. Also mention the importance of the lateral packing (i.e., the number of helices per unit area in the AB plane).
High density polyethylene has an equilibrium melting temperature of 138 °C. Two identical batches of this polymer are crystallized at different temperatures: one at 100 °C and one at 125 °C. Describe one important difference in the resulting structures of these two semicrystalline polymers, and explain the origin of the difference.
A linear homopolymer was crystallized from the melt at crystallization temperatures ( T c ) within the range 270K to 330K. Following complete crystallization the following melting temperatures ( T m ) were measured by differential scanning calorimetry (DSC).
Given:
T c (K)
T m (K)
270
300
280
306.5
290
312.5
300
319.0
310
325.0
320
331.0
330
337.5
Graphically determine the equilibrium melting temperature, T m 0 .
Small angle x-ray scattering experiments using x-rays with a wavelength of 1.54 Å were used to measure the distance between crystalline regions. Peaks corresponding to the lamellar period were obtained at the following angles:
T c ( K )
θ ( degrees )
270
0.44
280
0.39
290
0.33
300
0.27
310
0.22
320
0.17
330
0.11
  1. Use the Bragg equation to calculate the values of the lamellar period for each crystallization temperature.
  2. The degree of crystallinity was measured to be 45% for all samples. Calculate the lamellar thickness in each case, and use this information to graphically obtain the fold surface energy, γ e . Assume the enthalpy of fusion of the polymer is 150 J/c m 3 , and the value of the equilibrium melting temperature from the previous problem.
Poly(ethylene oxide) (PEO) is a semicrystalline polymer with melting point of 65 C and a glass transition temperature of -60 C. Suppose a series of crystalline samples is made by performing the crystallization at temperatures of 40 C and 60 C.
  1. Describe any differences you expect rate of crystallization, and the crystal thickness c , at these two temperatures. Be as quantitative as possible.
  2. Suppose the two samples are slowly heated. Describe any differences you expert to observe in the melting temperatures for these two samples.
  3. PEO is always crystalline at room temperature, but PET is often amorphous at room temperature, even though PET has a much higher melting temperature than PEO. Why is this?

Liquid Crystals

Consider the liquid crystalline texture from the nematic liquid crystal texture page of a nematic liquid crystalline material confined between two crossed polarizers.
  1. Give an example of the type of defect that can lie at the places where 4 dark bands come together. (show how the molecules are oriented in these regions)
  2. Describe why alternating light and dark bands emanate from these defect points.
  3. Is it possible to obtain a completely black image from a sample that has good nematic order? If so, describe a situation where this would be the case.

Thermodynamics of Solutions and Blends

Which contribution to the free energy (enthalpy of entropy) dominates at low temperature?
What are the signs of h mix and s mix for a system which has an LCST phase diagram?
For what values of χ will a polymer be able to completely dissolve in a small-molecular solvent?
Measured values of χ for a blend of polymer A with polymer B have the following temperature dependence:
χ = 0.326 T -2.3 × 1 0 -4
where T is the absolute temperature.
Consider a blend with N a =2500, N b =3500
  1. What is the predicted critical composition, and the predicted critical temperature for this blend?
  2. Plot free energy curves (free energy as a function of composition) for χ / χ crit = 1.1, 1.2, 1.3, 1.4, and 1.5. Using the common tangent construction from these five curves, construct a phase diagram with χ / χ crit on the vertical axis and the blend composition on the horizontal axis.
Consider the following phase diagram:
image: 245_home_ken_Mydocs_MSEcore_331_figures_28-1-problem.png
  1. Is the enthalpy of mixing positive or negative for this system?
  2. Is the entropy of mixing positive or negative for this system?
  3. Draw free energy of mixing vs. composition curves at 10 0 ,15 0 ,20 0 that are consistent with this phase diagram. Show any common tangents points and points of inflection.
Consider a polymer/solvent mixture for which the ratio of polymer to solvent volume is 10,000.
  1. Determine the critical value of the chi parameter ( χ crit )
  2. What is polymer volume fraction at the critical point?
  3. Plot the solvent chemical potential for this system as a function of polymer volume fraction for χ / χ crit =0.9, χ / χ crit =1.0, χ / χ crit =1.2, χ / χ crit =1.5 .
  4. Suppose that the polymer is dipped into a pure solvent (solvent chemical potential equals zero). Solvent will diffuse into the polymer until the solvent chemical potential in the polymer reaches zero. From the plots made for part c, determine the equilibrium solvent volume fractions for χ / χ crit =1.2 , and χ / χ crit =1.5.
Describe in your own words why most pairs of different types of polymers do not mix with each other in the liquid state.
Most liquids mix better at high temperatures than they do at low temperatures. However, polymeric liquids often mix better at low temperatures than high temperatures. Explain in as much detail as possible why this is the case.
Consider the following osmotic pressure data obtained for polyethylene oxide solutions in water (1MPa = 1 0 6 Pa) at 20 °C.
φ P
Π ( MPa )
0.0457
0.0411
0.0917
0.145
0.138
0.346
0.185
0.685
0.232
1.212
0.280
2.010
0.329
3.161
0.377
4.868
0.427
7.323
0.476
10.77
0.516
14.54
  1. Plot the osmotic pressure as a function of the polymer volume fraction, using logarithmic axes for both the pressure and the polymer volume fraction. On the same plot, include the Flory Huggins expression for Π , using the value of χ that gives the best fit to the data.
  2. Repeat part a, but now compare the data to the scaling expression for φ in the semidilute regime (Eq.), where a o 3 / v s is the adjustable parameter. (Note that v s is the solvent molecular volume, obtained from the molar volume, V s , by dividing by N av ).
Use the PEO solution data from the previous problem to calculate the following:
  1. The concentration of PEO at 2 0 C that gives an osmotic pressure large enough to support a column of water that is 100 meters tall.
  2. What is the correlation length ξ , for a PEO solution in water with φ P =0.1?
For charged molecules, the osmotic pressure is often dominated by the counterions. Suppose a solution of sodium polyacrylate at a concentration of 0.05 g/c m 3 is made in water. Estimate the osmotic pressure, assuming that the counterions can be treated as an ideal gas.
image: 246_home_ken_Mydocs_MSEcore_331_figures_polysodiumacrylate.png
The interfacial free energy characterizing the interface high molecular weight versions of polystyrene (PS) and poly(2-vinyl pyridine) (P2VP) is 3x1 0 -3 J/m 2 at 180 °C. Assume the volume per repeat unit is the same for both polymers (165 3 , and use this as your reference volume. Also assume that the statistical segment length for both polymers is 6.7 Å.
PS
image: 84_home_ken_Mydocs_MSEcore_331_figures_polystyrene.png
P2VP
image: 247_home_ken_Mydocs_MSEcore_331_figures_poly2vinylpyridine.png
  1. What is χ for the PS/P2VP system at 180 °C?
  2. Suppose PS and P2VP polymers with equal degrees of polymerization are mixed with one another. What range of molecular weights will form mixtures that are completely miscible with one another at all compositions at 180 °C?
  3. Suppose monodisperse PS and P2VP samples with degrees of polymerization of 1000 are mixed with one another. What is the solubility of the PS in the P2VP at 180 °C? (hint: for a symmetric system like this, where both degrees of polymerization are equal to one another, the common tangent construction is satisfied at the point where the first derivative of the free energy of mixing with respect to volume fraction is equal to zero).
Use the data from the PS/P2VP system from the previous problem to answer the following questions:
  1. A diblock copolymer is synthesized that has a PS block with a degree of polymerization of 500 attached to a 2nd block with a degree of polymerization of 500. Calculate the periodicity of the lamellar structure that is formed by this material at equilibrium.
  2. Describe in qualitative terms the equilibrium geometries of the structures formed when the P2VP block length is fixed, but the PS block degree of polymerization is steadily increased.
Poly(n-isopropylacrylamide) (chemical structure shown below) exhibits a lower critical solution temperature in water. It is fully soluble in water at temperatures below 32 C, but the phase diagram with water has a two phase region for temperatures above 32 C.
image: 248_home_ken_Mydocs_MSEcore_331_figures_polyisopropylacrylamide.png
  1. The temperature dependence of χ for the Poly(n-isopropylacrylamide)/water system can be written in the following form: χ =A+BT , where A and B are constants and T is the temperature in C. What is the sign (positive or negative) of B in this equation?
  2. Is the entropy of mixing for this system positive or negative? Briefly explain your answer.
  3. Calculate the osmotic pressure of an aqueous solution of high molecular weight poly(n-isopropylacrylamide) with a polymer volume fraction of 0.2.

Capillarity

A hollow glass tube is immersed in water as shown in the illustration below:
image: 249_home_ken_Mydocs_MSEcore_331_figures_44-1-problem.png
  1. Calculate the internal tube diameter that gives a capillary rise of 10 cm, assuming that water completely wets the glass (contact angle of zero).
  2. Calculate the capillary rise if the same glass tube from part a is coated with a thin layer of polystyrene, changing the contact angle to 90 .
Answer the following:
  1. Calculate the force required to push a cylindrical fiber with a radius of 1 μ m into water, if the contact angle the fiber makes with the water is 13 0 .
  2. What stress can be supported if you create an array of these fibers on a square lattice, with a spacing between fiber centers of 4 μ m ? (Assume the density of the fibers is equal to the density of water, so you don't need to worry about forces due to buoyancy).
  3. How much weight could be supported by an object like this with a total area of 1 c m 2 ?

331 Computational Exercises

In binary blends where the two components have very different molar volumes (often the case with polymers), the following form of the free energy expression is commonly used:
v 0 G v k B T = φ b ln φ b N b + ( 1- φ b ) ln ( 1- φ b ) N a + χ φ b ( 1- φ b )
The expressions for the chemical potentials for the A and B molecules are given by the following expression:
μ a k B T = ln ( 1- φ b ) + φ b ( 1- N a / N b ) + χ N a φ b 2
μ b k B T = ln ( φ b ) +( 1- φ b ) ( 1- N b / N a ) + χ N b ( 1- φ b ) 2 Here G v is the free energy per volume, φ b is the volume fraction of B in the system, and N a and N b are proportional to the molecular volumes of the two components. The quantity v 0 is a reference volume, so the molecular volume of the ' A ' molecule is v 0 N a and the molecular volume of the ' B ' molecule is v 0 N b . Adapt the MATLAB script below to do the following:
global na nb chi % these values also get used in function definitions
set(0,'defaultaxesfontsize',16)
set(0,'defaultlinelinewidth',2)
phi=linspace(0.001,0.999,1000);  
na=150;
nb=100;
chi=0.02;
% fv is the expression for the free energy of mixing
fv=@(phi) phi.*log(phi)/nb+(1-phi).*log(1-phi)./na+chi.*phi.*(1-phi);

%mua and  mub are the chemical potentials of A and B
mua=@(phi) log(1-phi)+phi*(1-na/nb)+chi*na*phi.^2; 
mub=@(phi) log(phi)+(1-phi)*(1-nb/na)+chi*nb*(1-phi)^2;

% now write the function that is equal to zero when the A and B chemical
% potentials are equal to one another for phi=phi(1) and phi=phi(2)

ftosolve=@(phi) [mua(phi(1))-mua(phi(2));
    mub(phi(1))-mub(phi(2))];

plot(phi,fv(phi));
xlabel('\phi_{b}')
ylabel('G_{v}V_{0}/RT')

% start with a guess for the equilibrium volume fractions
phiguess = [0.2; 0.8];     % Make a starting guess at the equilibrium compositions
[phicalc,fval] = fsolve(ftosolve,phiguess);  % Call solver 

% now we add the tangent line
slope=(fv(phicalc(2))-fv(phicalc(1)))/(phicalc(2)-phicalc(1));
intercept=fv(phicalc(1))-slope*phicalc(1);
tangentline=intercept+slope.*phi;
hold on
plot(phi,tangentline,'r')
hold off
title(['\chi=' num2str(chi) ', N_a=' num2str(na), ', N_b=' num2str(nb)], 'fontsize', 11)

% now save the plot as a jpg file
saveas(gcf,'commontangent.jpg');

% this saves the file a .eps file, used to embed the figure into the
% solution set - students can comment out this next line if they don't wan
% the .svg file
print(gcf,'-dsvg','commontangent.svg')



  1. Plot the free energy expression χ =0.02 , N a =200 , and N b =100 , along with its common tangent construction, and report the values of the equilibrium compositions of the two phases. Include the plot with your solution.
  1. Repeat the previous calculation for χ =0.025 and χ =0.05 (you may need to adjust the initial guess for the volume fractions).
  2. Decrease χ to the critical value, χ crit , which can be calculated analytically. Include a plot of G v as a function of φ b for value of χ .
Extend work from the previous problem to obtain a MATLAB code that generates the full phase diagram (with χ on the vertical axis and φ b on the horizontal axis) for a system with N a =200 and N b =100 . Show both the bimodal and spinodal curves, and turn in a printout of your actual MATLAB code. Hint: You'll want to solve the equations for a variety of χ values above the critical value of χ , which you can calculate analytically. You can use a 'for' loop in MATLAB to do this. Calculation of the spinodal curve is pretty easy, because there's an analytic expression for that. When you're developing your initial guesses for the bimodal compositions to pass to the 'fsolve' command, keep in mind that the bimodal compositions must lie outside the spinodal compositions. This can help you develop sensible guesses to use for this quantity.

Appendix 1: MATLAB scripts

25.1 Symbolic Math

The following script uses MATLAB's symbolic math capability to evaluate the integrals in Eq. eq:maxwell model.
% symbolic_visco.m
syms tau w go t real positive % specify symbolic variables and insist that they be real and positive
g=go*exp(-t/tau); % define exponential  relaxation function
storage=w*int(g*sin(w*t),t,0,inf); % evaluate the ingetral for G'
loss=w*int(g*cos(w*t),t,0,inf); % evaluate the integral for G"
pretty(storage) % visualize the expressions in a simpler form
pretty(loss)
The 'pretty' command returns the solutions in a format that is relatively easy to interpret, and is shown in Figure 186. Note that the first expression (for G ) needs to be manipulated a little bit to give the form shown in Eq. eq:maxwell model.
image: 250_home_ken_Mydocs_MSEcore_331_figures_matlab_g_integration.png
Figure 186: Output generated by symbolicvisco.m.

25.2 Solution to a Nonlinear Equation

The following script returns φ s =0.2885 as the solution to Eq. 21.9 for χ =1.05 .
chi=1.05;  % this is the value of chi we want to use
mus=@(phis)log(phis)+1-phis+chi*(1-phis)^2; % solvent chemical potential expression
phiguess=0.2;  % initial guess for phiguess
phisval=fsolve(mus,phiguess)  % returns the value of phis for which mus=0, with an initial guess for phis of phiguess

26 Appendix 2: Data Tables

26.1 Dynamic Moduli

26.1.1 Poly(2-vinyl pyridine)

130 C
log ( ω ) (s -1 )
log ( G ) (kPa)
log ( G ) (kPa)
-0.8
4.48
4.25
-0.6
4.55
4.25
-0.4
4.6
4.3
-0.2
4.66
4.28
0
4.71
4.35
0.2
4.76
4.39
0.4
4.81
4.43
0.6
4.87
4.46
0.8
4.92
4.49
1
4.97
4.5
1.2
5.02
4.49
1.4
5.05
4.44
1.6
5.08
4.36
1.8
5.12
4.23
2
5.13
4.3
140 C
log ( ω ) (s -1 )
log ( G ) (kPa)
log ( G ) (kPa)
-1
3.6
4.01
-0.8
3.88
4.13
-0.6
4.12
4.2
-0.4
4.28
4.24
-0.2
4.4
4.25
0
4.48
4.26
0.2
4.54
4.25
0.4
4.6
4.27
0.6
4.65
4.3
0.8
4.71
4.33
1
4.76
4.36
1.2
4.81
4.4
1.4
4.87
4.44
1.6
4.91
4.44
1.8
4.97
4.45
2
5.02
4.44
160 C
log ( ω ) (s -1 )
log ( G ) (kPa)
log ( G ) (kPa)
-1
2.1
2.83
-0.8
2.21
3.03
-0.6
2.36
3.22
-0.4
2.58
3.42
-0.2
2.86
3.6
0
3.19
3.79
0.2
3.5
3.96
0.4
3.8
4.1
0.6
4.05
4.19
0.8
4.24
4.24
1
4.37
4.25
1.2
4.46
4.25
1.4
4.53
4.25
1.6
4.6
4.26
1.8
4.65
4.26
2
4.7
4.29
180 C
log ( ω ) (s -1 )
log ( G ) (kPa)
log ( G ) (kPa)
-0.4
1.72
2.53
-0.2
1.79
2.73
0
1.9
2.92
0.2
2.06
3.12
0.4
2.31
3.32
0.6
2.59
3.51
0.8
2.92
3.7
1
3.26
3.88
1.2
3.59
4.05
1.4
3.9
4.19
1.6
4.16
4.27
1.8
4.35
4.31
2
4.47
4.3
200 C
log ( ω ) (s -1 )
log ( G ) (kPa)
log ( G ) (kPa)
0.4
1.71
2.63
0.6
1.81
2.82
0.8
1.97
3.02
1
2.18
3.21
1.2
2.45
3.41
1.4
2.77
3.6
1.6
3.11
3.79
1.8
3.46
3.97
2
3.79
4.11

26.1.2 Poly(t-butyl acrylate)

75 C
log ( ω ) (s -1 )
log ( G ) (kPa)
log ( G ) (kPa)
-1
5.01
4.37
-0.83
5.03
4.39
-0.67
5.05
4.42
-0.5
5.07
4.46
-0.33
5.09
4.51
-0.17
5.11
4.58
0
5.13
4.66
0.17
5.16
4.74
0.33
5.19
4.83
0.5
5.23
4.92
0.67
5.26
5.02
0.83
5.31
5.11
1
5.36
5.21
1.17
5.41
5.31
1.33
5.47
5.41
1.5
5.53
5.51
1.67
5.6
5.61
1.83
5.68
5.72
2
5.76
5.82
85 C
log ( ω ) (s -1 )
log ( G ) (kPa)
log ( G ) (kPa)
-1
4.88
4.39
-0.83
4.91
4.38
-0.67
4.94
4.37
-0.5
4.97
4.37
-0.33
4.99
4.37
-0.17
5.01
4.38
0
5.03
4.4
0.17
5.06
4.44
0.33
5.08
4.48
0.5
5.1
4.54
0.67
5.12
4.61
0.83
5.14
4.69
1
5.17
4.77
1.17
5.2
4.86
1.33
5.24
4.95
1.5
5.28
5.05
1.67
5.33
5.14
1.83
5.38
5.24
2
5.43
5.34
100 C
log ( ω ) (s -1 )
log ( G ) (kPa)
log ( G ) (kPa)
-1
4.6
4.41
-0.83
4.66
4.42
-0.67
4.72
4.42
-0.5
4.78
4.41
-0.33
4.82
4.41
-0.17
4.86
4.4
0
4.9
4.39
0.17
4.93
4.38
0.33
4.95
4.38
0.5
4.98
4.38
0.67
5
4.39
0.83
5.03
4.4
1
5.05
4.43
1.17
5.06
4.47
1.33
5.09
4.52
1.5
5.11
4.59
1.67
5.13
4.66
1.83
5.16
4.74
2
5.19
4.83
120 C
log ( ω ) (s -1 )
log ( G ) (kPa)
log ( G ) (kPa)
-1
3.79
4.15
-0.83
4.01
4.25
-0.67
4.20
4.32
-0.5
4.35
4.37
-0.33
4.46
4.40
-0.17
4.56
4.41
0
4.63
4.42
0.17
4.70
4.42
0.33
4.75
4.42
0.5
4.80
4.42
0.67
4.84
4.41
0.83
4.88
4.41
1
4.91
4.40
1.17
4.94
4.39
1.33
4.97
4.39
1.5
5.00
4.40
1.67
5.02
4.41
1.83
5.04
4.44
2
5.06
4.47
140 C
log ( ω ) (s -1 )
log ( G ) (kPa)
log ( G ) (kPa)
-1
2.72
3.57
-0.83
2.98
3.73
-0.67
3.24
3.88
-0.5
3.50
4.02
-0.33
3.75
4.14
-0.17
3.98
4.25
0
4.17
4.32
0.17
4.33
4.38
0.33
4.45
4.41
0.5
4.55
4.42
0.67
4.63
4.43
0.83
4.70
4.44
1
4.75
4.44
1.17
4.80
4.44
1.33
4.85
4.43
1.5
4.89
4.42
1.67
4.92
4.42
1.83
4.95
4.41
2
4.98
4.41

27 Other Resources

27.1 Polymer Synthesis

The Macrogalleria web site is an excellent online source describing the reaction mechanisms for most of the polymerization reactions that we have considered:
http://www.pslc.ws/macrog/level4.htm
This site includes reaction mechanisms for free radical, anionic anionic, cationic and step growth polymerizations.

References

1"Differential Scanning Calorimetry".
2, "Teflon", Wikipedia, the free encyclopedia ().
3"The Glass Transition".
4"Thermal Transitions of Homopolymers: Glass Transition & Melting Point".
5, "Seawater", Wikipedia, the free encyclopedia (2015).
6Frank S. Bates and Glenn H. Fredrickson, "Block Copolymers—Designer Soft Materials", Physics Today 52, 2 (2008), pp. 32--38.
7Hamdy F.M. Mohamed, E.E. Abdel-Hady, and Salwa S. Mohamed, "Temperature Dependence of the Free Volume in Polytetrafluoroethylene Studied by Positron Annihilation Spectroscopy", Radiation Physics and Chemistry 76, 2 (2007), pp. 160--164.
8 Wikipedia, "Torsion Spring", Wikipedia, the free encyclopedia ().

Index