Polymer Solutions, Solvent Quality & Flory-Huggins Lattice Thermodynamics

The dissolution of a high molecular weight polymer in a low molecular weight solvent is governed by the fundamental laws of chemical thermodynamics. Spontaneous mixing at constant temperature $T$ and pressure $P$ requires that the Gibbs free energy of mixing be negative:

However, a negative $\Delta G_m$ is merely a necessary condition for dissolution; it is not sufficient to guarantee complete, single-phase thermodynamic stability across all concentrations.

Thermodynamic Criteria for Phase Stability

For a binary mixture to remain permanently homogeneous without spontaneous phase separation into two liquid phases (liquid-liquid demixing), the free energy of mixing curve must be strictly concave upward (positive second derivative with respect to composition) across the entire composition range:

where $\phi_2$ is the volume fraction of the polymer. If there exists an intermediate composition interval where the curvature is negative:

the homogeneous solution is intrinsically thermodynamically unstable against spontaneous infinitesimal composition fluctuations, triggering spontaneous phase separation via spinodal decomposition.

The Fundamental Challenge of Macromolecular Mixing

In simple low molecular weight binary liquids (e.g., benzene + toluene), the combinatorial entropy of mixing $\Delta S_m$ is large and positive, driving spontaneous dissolution even when mixing is moderately endothermic ($\Delta H_m > 0$). In polymers, however, connecting $N$ monomer segments together into an unbroken covalent chain severely restricts their translational freedom. When $N$ independent solvent molecules are replaced by one $N$-mer chain, the translational combinatorial entropy of mixing per unit volume is reduced by a factor of nearly $1/N \sim 10^{-3} - 10^{-5}$. Because $T \Delta S_m$ is exceptionally small, polymer dissolution is dominated by the enthalpy of mixing $\Delta H_m$. Even a very modest endothermic interaction energy per segment ($k_B T \chi$) is sufficient to overcome the tiny combinatorial entropy, driving macroscopic phase separation.

Formulated independently by Paul Flory and Maurice Huggins in 1941–1942, the Flory-Huggins lattice theory revolutionized polymer thermodynamics by calculating the combinatorial entropy of placing flexible polymer chains onto a regular lattice.

The Lattice Geometry

Consider a lattice with total volume $V$ divided into $N_0$ identical spatial cells each of volume $v_0$ (equal to the volume of a single solvent molecule). The lattice has coordination number $z$ (number of nearest-neighbor cells adjacent to any site). The binary solution consists of:

• $n_1$ solvent molecules, each occupying $1$ cell ($r_1 = 1$).
• $n_2$ polymer chains, each consisting of $x$ connected segments each occupying $1$ cell ($r_2 = x$).

The total number of lattice cells is:

The volume fractions of solvent ($\phi_1$) and polymer ($\phi_2$) are:

Sequential Placement & Number of Microstates $\Omega$

To calculate the total number of distinguishable configurations $\Omega$, Flory calculated the number of ways $\nu_{i+1}$ to insert the $(i+1)$-th polymer chain into the lattice when $i$ chains already occupy $i x$ cells:

For the $(i+1)$-th chain:

• The first segment has $N_0 - i x$ empty cells available.
• The second segment must occupy one of the $z$ adjacent cells, with probability of finding an empty cell given by the mean vacancy fraction $f_i = (N_0 - i x) / N_0$.
• The third and subsequent segments have $(z - 1)$ directional options, each with probability $f_i$.

Multiplying these factors and applying Stirling's approximation ($\ln n! \approx n \ln n - n$) yields the celebrated Flory-Huggins combinatorial entropy of mixing:

In molar units (where $N_A k_B = R$):

For monomeric liquid mixtures ($x = 1$), $\phi_1 = X_1$ and $\phi_2 = X_2$, recovering ideal Raoult's law mixing. When $x \gg 1$, the $n_2 \ln \phi_2$ term is heavily suppressed relative to $n_1 \ln \phi_1$, formalizing the dramatic loss of combinatorial entropy.

The energetic change accompanying mixing is modeled through a mean-field Bragg-Williams nearest-neighbor contact exchange mechanism.

Nearest-Neighbor Contact Exchange Energetics

Prior to mixing:

• Pure solvent contains only $1-1$ molecular contacts with contact energy $w_{11}$.
• Pure polymer contains only $2-2$ segment contacts with contact energy $w_{22}$.

Upon mixing, new $1-2$ contacts are created at the expense of $1-1$ and $2-2$ contacts:

The energy change associated with forming a single solvent-segment contact pair is:

Definition of the Interaction Parameter $\chi$

The dimensionless Flory-Huggins interaction parameter $\chi$ (or $\chi_{12}$) is defined as the exchange energy $\Delta w_{12}$ per solvent molecule scaled by thermal energy $k_B T$:

where $z$ is the lattice coordination number. In mean-field approximation, the average number of solvent-polymer contact pairs in the mixture is $z n_1 \phi_2$. Therefore, the total enthalpy of mixing is:

In molar units:

Physical Meaning and Temperature Dependence of $\chi$

From regular solution theory and Hildebrand solubility parameters:

where $V_1$ is the molar volume of the solvent. In real polymer solutions, $\chi$ is not purely enthalpic; it also contains an excess entropic contribution $\chi_s$ arising from non-random packing, specific orientation constraints, and solvent local ordering:

Typical empirical values of $\chi_s$ range from $0.2$ to $0.4$.

• $\chi < 0$: Highly favorable exothermic mixing (strong specific interactions such as hydrogen bonding).
• $0 \le \chi < 0.5$: Athermal to moderately endothermic mixing in a good solvent.
• $\chi = 0.5$: The exact theta ($\Theta$) state where repulsive excluded volume is canceled by attractive interactions.
• $\chi > 0.5$: Poor solvent; phase separation occurs above critical polymer concentrations.

Combining the combinatorial entropy and the contact exchange enthalpy yields the complete Flory-Huggins Free Energy of Mixing:

Dividing by the total number of moles of lattice sites $N_0 = n_1 + x n_2$:

Solvent Chemical Potential $\mu_1 - \mu_1^\circ$

The chemical potential of the solvent relative to pure solvent is obtained by differentiating $\Delta G_m$ with respect to $n_1$ at constant $T, P, n_2$:

Applying the chain rule:

Dilute Solution Series Expansion & Osmotic Pressure

In dilute polymer solutions ($\phi_2 \ll 1$), we expand $\ln(1 - \phi_2)$ in a Taylor series:

Substituting into the solvent chemical potential:

From classical thermodynamics, the osmotic pressure $\Pi$ is related to solvent chemical potential by $\Pi V_1 = -(\mu_1 - \mu_1^\circ)$, where $V_1$ is the molar volume of the solvent:

Substituting the mass concentration $c = \phi_2 / \bar{v}$ (where $\bar{v}$ is the polymer specific volume, so $V_1 / (x \bar{v}) = M_2$):

where the second virial coefficient $A_2$ is:

• In a good solvent ($\chi < 0.5$): $A_2 > 0$ (osmotic pressure exceeds ideal van 't Hoff value due to excluded volume repulsion).
• At the theta temperature ($\chi = 0.5$): $A_2 = 0$ (ideal thermodynamic behavior; $\Pi/c = R T / M_2$ across finite concentrations).
• In a poor solvent ($\chi > 0.5$): $A_2 < 0$ (attractive interactions dominate, leading toward phase separation).

The phase behavior of polymer solutions as a function of temperature and composition is characterized by liquid-liquid phase separation envelopes.

The Spinodal Condition

The boundary between metastable and intrinsically unstable thermodynamic states is the spinodal curve, defined by the vanishing second derivative of the free energy of mixing:

Differentiating the Flory-Huggins free energy density $\frac{\Delta G_m}{N_0 R T} = (1 - \phi_2)\ln(1 - \phi_2) + \frac{\phi_2}{x}\ln \phi_2 + \chi(1 - \phi_2)\phi_2$: First derivative:

Second derivative:

Solving for the interaction parameter $\chi$ along the spinodal:

The Critical Point

The critical point occurs where the third derivative of the free energy also vanishes (the inflection point of the spinodal curve):

Rearranging:

Solving for the critical polymer volume fraction $\phi_{2,c}$:

For high polymers where $x \gg 1$:

Substituting $\phi_{2,c}$ back into the spinodal equation yields the critical interaction parameter $\chi_c$:

In the limit of infinite molecular weight ($x \to \infty$):

This reveals a crucial hallmark of polymer thermodynamics: the critical point is strongly asymmetric, shifting to extremely dilute polymer concentrations ($\phi_{2,c} \to 0$).

The Binodal (Coexistence) Curve

The true phase boundary is the binodal curve, defined by the equality of chemical potentials of each component in both coexisting phases ($'$ and $''$):

Between the binodal and spinodal lies the metastable nucleation-and-growth regime. Inside the spinodal curve lies the spontaneous spinodal decomposition regime, characterized by interconnected bicontinuous morphologies.

The physical conformation and dimensions of a macromolecule in solution depend intimately on the thermodynamic interaction between monomer segments and the solvent.

The Three Thermodynamic Regimes of Solvent Quality

1. Good Solvent ($\chi < 0.5, A_2 > 0$):

• Monomer-solvent contacts are energetically favored over monomer-monomer contacts.
• The chain expands to maximize solvent contact, exhibiting repulsive excluded volume.
• The coil dimension swells beyond unperturbed Gaussian dimensions, following Flory's self-avoiding walk scaling:

(Renormalization group theory yields $\nu \approx 0.588$).

2. Theta ($\Theta$) Solvent ($\chi = 0.5, A_2 = 0$):

• The thermodynamic temperature where the attractive segment-segment van der Waals interactions exactly balance the steric repulsive excluded volume.
• The effective excluded volume parameter $v_{\text{ex}}$ vanishes:

• The polymer behaves as an ideal, non-interacting random walk (phantom chain):

3. Poor Solvent ($\chi > 0.5, A_2 < 0$):

• Segment-segment attractive interactions dominate over solvent interactions.
• The coil contracts below unperturbed dimensions into a dense, collapsed globule:

• At temperatures below $\Theta$, macroscopic phase separation into polymer-rich and polymer-poor phases occurs.

The Linear Expansion Factor $\alpha$

The swelling of a real polymer coil relative to its unperturbed theta state is quantified by the linear expansion factor $\alpha$:

Flory derived the fundamental closed-form equation for chain expansion:

where $C_T$ is a molecular constant dependent on chain flexibility and molar volume.

• At $T = \Theta$: $\alpha = 1$ (unperturbed state).
• For $T > \Theta$: $\alpha > 1$ (swollen coil in good solvent).
• For $T < \Theta$: $\alpha < 1$ (contracted coil).

Polymer-solvent systems exhibit two distinct types of liquid-liquid phase separation envelopes governed by temperature.

Upper Critical Solution Temperature (UCST)

• The classical Flory-Huggins model with enthalpic interaction $\chi \propto 1/T$ predicts that raising the temperature decreases $\chi$, increasing solvent quality and restoring single-phase homogeneity.
• Phase separation occurs upon cooling below a critical temperature, termed the Upper Critical Solution Temperature (UCST):

• As molecular weight approaches infinity ($x \to \infty$), $T_c^{\text{UCST}} \to \Theta$.
• Classical UCST systems: Polystyrene in cyclohexane ($\Theta = 34.5^\circ\text{C}$), PMMA in 1-chlorobutane.

Lower Critical Solution Temperature (LCST)

• Many polymer solutions, including water-soluble polymers (e.g., poly(N-isopropylacrylamide), PNIPAM; poly(ethylene oxide), PEO) and hydrocarbon solutions at elevated temperatures near the solvent vapor-liquid critical point, exhibit phase separation upon heating.
• The boundary temperature above which two phases appear is the Lower Critical Solution Temperature (LCST).
• Physical Mechanisms of LCST:

1. Aqueous Systems (Hydrophobic Effect): Below the LCST, dissolution is driven by exothermic hydrogen bonding between water and polar groups, accompanied by rigid, highly ordered iceberg water cages surrounding non-polar moieties (negative excess entropy $\Delta S^E < 0$). Upon heating, hydrogen bonds dissociate and the release of structured water cages provides a large entropic gain for phase separation ($T \Delta S^E < 0$). For PNIPAM, the LCST occurs sharply at $32^\circ\text{C}$.

2. Non-Aqueous Systems (Free Volume Disparity): Near the boiling point of the solvent, the thermal expansion coefficient of the solvent liquid far exceeds that of the tethered polymer chain. This macroscopic free volume mismatch produces an unfavorable equation-of-state compressibility entropy penalty upon mixing, driving phase separation at high temperatures.

Because synthetic polymers are polydisperse (possessing a broad distribution of molecular weights), the strong dependence of the critical interaction parameter $\chi_c$ and binodal coexistence curve on chain length $x$ provides the fundamental mechanism for polymer fractionation by solubility.

Principle of Fractional Precipitation

From the Flory-Huggins partition coefficient:

where $\phi_{2,x}''$ is the volume fraction of species $x$ in the precipitated concentrated phase, $\phi_{2,x}'$ is its fraction in the dilute supernatant phase, and $\sigma$ is a positive thermodynamic separation parameter proportional to the degree of supersaturation:

Because the partition coefficient scales exponentially with degree of polymerization $x$, the highest molecular weight species preferentially partition into the concentrated precipitated phase:

Standard Fractionation Procedures

1. Non-Solvent Addition (Titration):

• A dilute ($< 1\text{ wt}\%$) polymer solution in a good solvent is stirred at constant temperature.
• A miscible non-solvent is added dropwise until faint turbidity appears (coacervation point).
• The mixture is warmed until clear, then cooled slowly to equilibrium.
• The dense, polymer-rich coacervate phase is isolated by decantation or centrifugation, recovering the highest molecular weight fraction.
• Repeating this process stepwise yields narrow molecular weight cuts.

2. Temperature-Induced Precipitation:

• For a UCST system, a dilute polymer solution is cooled in precise temperature increments $\Delta T = 0.5 - 2^\circ\text{C}$.
• The highest molecular weight chains precipitate first because their critical temperature $T_c(x)$ is highest.

3. Coacervation & Microencapsulation:

• Liquid-liquid phase separation where the polymer-rich phase separates as liquid droplets (coacervates) rather than solid precipitates; used widely in pharmaceutical drug microencapsulation and food colloids.