Flory-Huggins Theory

3.012 Fundamentals of Materials Science Fall 2005 Lecture 24: 12.07.05 Flory-Huggins Theory Today: LAST TIME .........................................................................................................................................................................................2 Lattice Models of Solutions .........................................................................................................................................................2 ENTROPY OF MIXING IN THE FLORY-HUGGINS MODEL ...................................................................................................................3 CONFIGURATIONS OF A SINGLE CHAIN.............................................................................................................................3 COUNTING CONFIGURATIONS FOR A COLLECTION OF CHAINS ................................................................................4 THE ENERGY OF MIXING IN THE FLORY-HUGGINS MODEL ..............................................................................................................7 The Bragg-Williams approximation ...........................................................................................................................................9 PREDICTIONS FROM THE FLORY-H UGGINS MODEL .......................................................................................................................11 REFERENCES ...................................................................................................................................................................................14 Reading: Dill and Bromberg Ch. 31 ‘Polymer Solutions,’ pp. 593-605. Supplementary Reading: Announcements: Final quiz: Tuesday, Dec. 20 1:30pm-4:30pm DuPont Same format and length as first 2 quizzes; you will have 3 hrs but quiz should only take 2 hrs Equation sheet is allowed, as in prior quizzes Review sessions (tentatively): Thermo/Stat Mech review: Friday, Dec. 16, 5-7pm (normal recitation room) Structure review: Sunday, Dec. 18, time TBD (normal recitation room) Thermo coverage on Quiz 3: lectures 17-24 3.012 Fundamentals of Materials Science Fall 2005 Last time Lattice Models of Solutions • Lattice models are a form of COARSE GRAINING, where only the most important molecular details of a system are retained in a statistical mechanical model. This is a powerful approach to create stat mech models where meaningful predictions can be made for complex materials. We began a derivation of the free energy of mixing for a polymer solution based on Flory-Huggins theory: np = number polymer chains N = number of segments per chain ns = number of solvent molecules Figure by MIT OCW. ! M = total lattice sites = ns + Nnp 3.012 Fundamentals of Materials Science Fall 2005 Entropy of mixing in the Flory-Huggins model CONFIGURATIONS OF A SINGLE CHAIN  How do we calculate W, the number of uniques states in the polymer system? We start by looking at the number of ways to place a single polymer chain on the lattice, ν 1: o What are the number of conformations for first bead? o With the first bead placed on the lattice, what is the number of possible locations for the second segment of the chain? o Moving on to placement of the third segment of the chain: 3.012 Fundamentals of Materials Science o Fall 2005 We repeat this process to place all N segments of the chain on the lattice, and arrive at ν 1, the total number of configurations for a single chain: COUNTING CONFIGURATIONS FOR A COLLECTION OF CHAINS o We can follow the same procedure used for a single chain to obtain the number of configurations possible for an entire set of np chains. We start by placing the FIRST SEGMENT OF ALL np CHAINS. The number of configurations for the first segment of all np chains is ν first: " first = M(M #1)(M # 2)(M # 3) $ $ $ (M # (n p #1)) = o M! (M # n p )! The number of configurations for the (N - 1) remaining segments of all np chains is ν subsequent: ! $ M # N ( n #1) # ( N #1) '* $ M # n p '-* $ M # n p #1'- * $ M # 2n p '- * p , )/ " subsequent = ,( z #1)& z #1 0 0 0 z #1 0 0 0 z #1 )& )& ( )&& )/,( )/ ,( )/ )/ M M M M % (.+ % (. + % (. ,+ + % (. ! " subsequent o $ z #1' n p (N #1) ( M # n p )! =& ) % M ( (M # Nn p )! Putting these two configuration counts together, we have the total number of configurations for the collection of np chains of N segments each: ! W = " first" subsequent n p!  ! The factor of np! Corrects for the over-counting since the polymer chains are indistinguishable, and we can’t tell the difference between two configurations with the same polymer distributions but different chain identities: Figure by MIT OCW. Pure Polymer Figure by MIT OCW. Figure by MIT OCW. 3.012 Fundamentals of Materials Science o Applying Stirling’s approximation: Fall 2005 ln x!" x ln x # x , …arriving at a final result: "#S mix = $k b [ n s ln % s + n p ln % p ] ! ! 3.012 Fundamentals of Materials Science Fall 2005 The energy of mixing in the Flory-Huggins model  Molecules form bonds with one another, attract, and repel one another by electrostatic, polar, or Van der Waals interactions- changing the internal energy of the system. Figure by MIT OCW. Figure by MIT OCW.  To account for these interactions in a lattice model, we need to first calculate the enthalpy of mixing. At constant pressure, the enthalpy of mixing will be related to the internal energy of mixing in the following way:  We will assume in the lattice model that the change of volume on mixing is zero- thus the second term mix mix drops out and ΔH = ΔU .  We can directly formulate an expression for the energy of mixing in the model:  The internal energy of the solution is obtained by considering the contacts between molecules on the lattice: mij =# i - j contacts " ij = energy per i - j contact # energy & U solution = " ( # contacts)% ( $ contact ' !  ! The total number of contacts made by polymer segments or solvent molecules on the lattice can be related to the coordination number and the contact numbers mij: 3.012 Fundamentals of Materials Science total solvent contacts: total polymer contacts:  Substituting into the expression for the internal energy of the solution, we have: U solution = m pp" pp + mss" ss + m ps" ps $ zNn p # m ps ' $ zn s # m ps ' U solution = & " + ) pp & )" BB + " psm ps 2 2 % ( % ( $ $ " + " ss '' zNn p zn U solution = " pp + s " BB + &" ps # & pp )) m ps 2 2 2 % (( %  ! The internal energy of the umixed state has a simpler form: U unmixed = ( # s - s contacts)" ss + ( # p - p contacts)" pp !  This leads us finally to an expression for the energy of mixing: "U mix = U solution # U unmixed ! Fall 2005 3.012 Fundamentals of Materials Science Fall 2005 The Bragg-Williams approximation  Our next step is to simplify the last equation- how can we determine mps from our model? o One way to achieve this is to make the approximation that the polymer and solvent molecules mix randomly on the lattice. In that case, all we need to do to find mps is calculate the probability of PS contacts for random mixing. This is called the mean field approximation or the BraggWilliams approximation. Note that we have already implicitly invoked this approximation in our excluded volume correction in the calculation of the entropy of mixing.  The probability that a lattice site contains a polymer segment is:  For a single solvent molecule, the average number of contacts with polymer segments is then:  Thus, the total number of P-S contacts is: 3.012 Fundamentals of Materials Science Fall 2005  This gives us a completely known expression for the energy of mixing in terms of the interaction energies:  From this expression, the molar enthalpy of mixing is:  We define the Flory-Huggins interaction parameter as: o  χ is a unitless equivalent of Ω, the interaction parameter from the regular solution model. Note that Ω is thus directly related to the molecular energy of interaction between the two components of a binary system! The change in internal energy now reduces to: 3.012 Fundamentals of Materials Science Fall 2005 Predictions from the Flory-Huggins model  Now, this expression, derived directly from our molecule-scale model, can be combined with the lattice model entropy of mixing derived above to provide the free energy of mixing for the regular solution: 0.05 "Smix 0.4 0 0 0.3 0.2 0.2 0.4 0.6 -0.05 N=100 "Hmix "Gmix free energy, enthalpy, entropy of mixing 0.5 enthalpy of mixing -0.1 free energy of mixing 0.1 -0.15 0 0 0.2 0.4 0.6 -0.1 -0.2 !p 0.8 "Gmix 1 -0.2 -0.25 !p 0.8 1 3.012 Fundamentals of Materials Science  Fall 2005 The Flory-Huggins model predicts major trends in the behavior of real polymer solutions: T (0C) 30 a One phase b 25 c 20 d 15 10 5 0.0 Experimental phase diagrams of polystyrene in cyclohexane. Molecular weights ( a ) 1.27 x 106, ( b ) 2.5 x 105, ( c ) 8.9 x 104, ( d ) 4.3 x 104. The dashed lines indicate the Flory-Huggins theory predictions for the first and third curves from the top. Two phases 0.1 0.2 � 0.3 0.4 0.5 Figure by MIT OCW. 3.5 I/TC x 105 3.4 3.3 3.2 3.1 0 0 2 4 6 (1/x1/2 . 1/2x) x 102 8 10 A plot of the reciprocal of the critical temperature against the molecular size function for polystyrene fractions in cyclohexane ( ) and for polyisobutylene fractions in diisobutyl ketone ( ). Figure by MIT OCW. 3.012 Fundamentals of Materials Science  Fall 2005 …And this theory is used to predict new behavior in polymers in current research: Image removed for copyright reasons. Scanned image of article: Gonzalez-Leon, Juan, Metin Acar, Ryu, Sang-Woog Ryu, and Ruzette, Anne-Valerie Ruzette, and Anne M. Mayes. "Low-temperature processing of 'baroplastics' by pressure-induced flow." Nature 426 (2003): 424-428. 3.012 Fundamentals of Materials Science References 1. 2. Dill, K., and S. Bromberg. 2003. Molecular Driving Forces, New York. Flory, P. J. 1953. Principles of Polymer Chemistry. Cornell University Press, Ithaca. Fall 2005