A change, denoted by Δ, is the value of a variable for a solution or mixture minus the values for the pure components considered separately. The objective is to find explicit formulas for ΔHm and ΔSm, the enthalpy and entropy increments associated with the mixing process.
The right-hand side is a function of the number of molesn1 and volume fraction φ1 of solvent (component1), the number of moles n2 and volume fraction φ2 of polymer (component 2), with the introduction of a parameter chiχ to take account of the energy of interdispersing polymer and solvent molecules. R is the gas constant and T is the absolute temperature. The volume fraction is analogous to the mole fraction, but is weighted to take account of the relative sizes of the molecules. For a small solute, the mole fractions would appear instead, and this modification is the innovation due to Flory and Huggins.
Derivation
We first calculate the entropy of mixing, the increase in the uncertainty about the locations of the molecules when they are interspersed. In the pure condensed phases — solvent and polymer — everywhere we look we find a molecule.[3] Of course, any notion of "finding" a molecule in a given location is a thought experiment since we can't actually examine spatial locations the size of molecules. The expression for the entropy of mixing of small molecules in terms of mole fractions is no longer reasonable when the solute is a macromolecularchain. We take account of this dissymmetry in molecular sizes by assuming that individual polymer segments and individual solvent molecules occupy sites on a lattice. Each site is occupied by exactly one molecule of the solvent or by one monomer of the polymer chain, so the total number of sites is
N1 is the number of solvent molecules and N2 is the number of polymer molecules, each of which has x segments.[4]
These are also the probabilities that a given lattice site, chosen at random, is occupied by a solvent molecule or a polymer segment, respectively. Thus
For a small solute whose molecules occupy just one lattice site, x equals one, the volume fractions reduce to molecular or mole fractions, and we recover the usual equation from ideal mixing theory.
In addition to the entropic effect, we can expect an enthalpy change.[5] There are three molecular interactions to consider: solvent-solvent w11, monomer-monomer w22 (not the covalent bonding, but between different chain sections), and monomer-solvent w12. Each of the last occurs at the expense of the average of the other two, so the energy increment per monomer-solvent contact is
The total number of such contacts is
where z is the coordination number, the number of nearest neighbors for a lattice site, each one occupied either by one chain segment or a solvent molecule. That is, xN2 is the total number of polymer segments (monomers) in the solution, so xN2z is the number of nearest-neighbor sites to all the polymer segments. Multiplying by the probability φ1 that any such site is occupied by a solvent molecule,[6] we obtain the total number of polymer-solvent molecular interactions. An approximation following mean field theory is made by following this procedure, thereby reducing the complex problem of many interactions to a simpler problem of one interaction.
The enthalpy change is equal to the energy change per polymer monomer-solvent interaction multiplied by the number of such interactions
The polymer-solvent interaction parameter chi is defined as
It depends on the nature of both the solvent and the solute, and is the only material-specific parameter in the model. The enthalpy change becomes
Assembling terms, the total free energy change is
where we have converted the expression from molecules N1 and N2 to moles n1 and n2 by transferring Avogadro's numberNA to the gas constantR = kNA.
where Vseg is the actual volume of a polymer segment.
This treatment does not attempt to calculate the conformational entropy of folding for polymer chains. (See the random coil discussion.) The conformations of even amorphous polymers will change when they go into solution, and most thermoplastic polymers also have lamellar crystalline regions which do not persist in solution as the chains separate. These events are accompanied by additional entropy and energy changes.
^ "Thermodynamics of High PolymerSolutions," Paul J. FloryJournal of Chemical Physics, August 1941, Volume 9, Issue 8, p. 660 Abstract. Flory suggested that Huggins' name ought to be first since he had published several months earlier: Flory, P.J., "Thermodynamics of high polymer solutions," J. Chem. Phys.10:51-61 (1942) Citation Classic No. 18, May 6, 1985
^ We are ignoring the free volume due to molecular disorder in liquids and amorphous solids as compared to crystals. This, and the assumption that monomers and solute molecules are really the same size, are the main geometric approximations in this model.
^ The enthalpy is the internal energy corrected for any pressure-volumework at constant (external) P. We are not making any distinction here. This allows the approximation of Helmholtz free energy, which is the natural form of free energy from the Flory-Huggins lattice theory, to Gibbs free energy.
^ In fact, two of the sites adjacent to a polymer segment are occupied by other polymer segments since it is part of a chain; and one more, making three, for branching sites, but only one for terminals.