Manufacturers

Thermodynamics of Size Exclusion Chromatography

The thermodynamics and the separation mechanism of SEC was first treated by Casassa, Gidding, Doi and Yau et al.1-9 Their treatments are based on the assumption that the solute molecules continuously move between the stationary phase (pore volume) and mobile phase (interstitial volume) and assume a dynamic equilibrium throughout the separation column. A thermodynamic equilibrium is reached when the chemical potential of the macromolecules in both phases become equal. At low concentrations, the polymer solutionis ideal. Then the chemical potential of the macromolecules in the mobile phase of concentration (cM) is given by

μM = μM,0 + kT · ln(cM / cM,0)

where μM,0 is the chemical potential of macromolecules in the reference state of concentration cM,0.

CROW logo

When a molecule diffuses into a pore (stationary phase) its entropy changes by ΔS and its enthalpy by ΔH. The change in entropy is caused by conformational and translational restrictions in the pores. In other words, a polymer squeezed into a pore can adopt a smaller number of conformations and the movement of the entire molecule is restricted by the pore walls. When a polymer enters a pore, surface–solvent contacts will be replaced by monomer–surface contacts, resulting in an enthalpy change. The change in enthalpy is due to interactions of the polymers with the pore surface. ΔH can be positive or negative depending on whether attractive or repulsive forces dominate. The chemical potential in the stationary phase (μS) of concentration (cS) is then given by

μS = μS,0 + kT · ln(cS / cS,0) - TΔS + ΔH

The concentration equilibrium is reached when μS = μM. Then

kT · ln(cS / cS,0) - TΔS + ΔH = kT · ln(cM / cM,0)

Rearranging the terms gives the partition coefficient K = cS/cM11,

K = exp(ΔS / k + ΔH / kT)

The enthalpic part ΔH is difficult to predict because it depends on the type of solvent, polymer and porous material. The entropic part ΔS, on the other hand, is universal because it is determined solely by the geometry and size of the pores and the hydrodynamic size of the polymer molecules. In an ideal SEC, the stationary phase has been designed to interact purely entropically with the eluent which requires a theta solvent and an inert porous packing. Then

K = exp(ΔS / k)

The separation effect of SEC also depends on the shape of the pores which are often of a highly convoluted complex shape. To understand the effect of pore size and pore geometry, several theoretical partition models have been developed that employ simple geometric equivalents of the polymers and the pore systems including (Gaussian) chains in slits, square tubes, cylinders as well as pores formed by randomly arranged spheres and randomly oriented planes.3,9,13,14 One of the oldest models was developed by Cassasa who calculated the partition coefficient K of ideal polymers in pores from the Green function G(r, r') (also called segment distribution function). The simplest pore geometry is a slit-shaped pore formed by two parallel plates of width d extending in x and y directions. In the case of a Gaussian chain, the segment distribution function G(r, r') can be factored into three independent components Gx, Gy and Gz. The segment distribution parallel to the plates (Gx, Gy) will not be affected by the slit. In other words, only Gz experiences a change. As has been shown by Cassada and others, the Gaussian segment distribution for this type of confinement is given by1-3, a

whereas the Gaussian probability distribution for an unperturbed chain is given by

The partition coefficient of a polymer in a slit Kslit can be calculated as the average of the Green function with respect to z and z':10

For narrow slits, Rg >> d, the first term dominates the series. Then

The figure below compares the partitioning coefficient of Gaussian chains with hard sphere molecules in slit-like pores of width d. The plot of K is given as a function of 2Rg/d. At 2Rg = d, the concentration in the pores is approximately 1/3 of the concentration in the mobile phase, whereas the concentration of hard spheres asymptotically approaches zero. For very large pores, (2Rg/d → 0), the partition coefficient K approaches unity for both hard spheres and flexible macromolecules. Then the molecule concentration inside the pores approaches asymptotically the concentration outside the pores (total permeation limit). Thus, the separation column becomes insensitive for polymers much smaller than the pores. If, on the other hand, the molecules are much larger than the pores, all macromolecules are excluded from the pores and the concentration inside the pores approaches zero (total exclusion limit). In other words, all molecules that pass through the column are eluted at approximately the same time tM and no separation occurs.

A real polymer chain in a pore will behave very differently than an ideal chain (chain in a theta solvent) since it is more extended in a confined space. We therefore expect real chains to have very different partition coefficients than ideal chains. Unfortunately, an exact treatment of real polymers in confinement is not possible, that is, their partition coefficients even in simple confining geometries are unknown. As has been shown by Teraoka11, it is possible, however, to obtain asymptotic expressions between ΔS and N for sufficiently long chains in simple geometrical confinement using scaling arguments.

Let us consider a real chain confined in a cylindrical pore of diameter d. If its radius of gyration (RgN3/5) is much smaller than the pore radius, the chain conformations are not much affected by the tube walls. However, when Rg exceeds the pore radius, the chain in the tube is forced to adopt conformations that are extended in the direction of the cylinder. As has been shown by Teraoka11, Rubinstein18 and others, a chain confined in a tube can be envisioned as a train of blobs of diameter ξ = d along the pore channel (see figure below).

CROW logo

Each of the partial chains within each blob assumes a conformation of a real chain in the absence of confinement. Then the size of the subchain follows the same scaling laws as those in unconfined space:

ξb g3/5

where b is the length of a statistical segment and g is the number of segments in a blob. The polymer chain consists of nd = N/g blobs each having the same size d. Assuming a simple cubic lattice and placing the spherical blobs on the grid points, each neighboring blob of the sequence of self-avoiding spheres can be placed on any of the five adjacent grid points, whereas in the tube the blobs have only one choice. Then the nd blobs in free space can be arranged in approximately 5nd different ways, whereas in the pore only one linear arrangement of blobs is possible. As has been shown above, the partition coefficient can be calculated as the exponent of the entropy, K = exp(S/kB), which is identical to the ratio of all possible chain conformations in the mobile phase to those in the pore. Since the monomers inside a blob are unperturbed, the loss in entropy is solely due to restriction in blob mobility. Then K can be calculated as the number of arrangements of blobs outside the pore (mobile phase) to those inside the pore:11

K ≅ 1/5nd

The change in entropy ΔS is given as

- ΔS/k = - lnKporend = N/gN (b/d)5/3 ≅ (Rg/d)5/3

or

Kpore ≅ exp[-(Rg/d)5/3]

In the case of a slit shaped pore, the blobs are arranged in a two-dimensional way. Then a sequence of blob in a slit can be placed on an imaginary 2-dimensional cubic grid in 3nd different ways. Using the same relationships we find that

K ≅ (3/5)nd

The arrangement of blobs in a slit is equivalent to a two-dimensional self-avoiding walk. The length of such a walk scales as

R||d (N/g)3/4d (Rg/d)5/4

where we made use of the scaling relationships g ≅ (d/b)5/3/b and Rgb N3/5. Similar scaling relationships can be derived for polymers in a θ-solvent; when ideal chains are squeezed into a cylindrical pore, the polymer blobs can assume more than one arrangement since neighboring blobs are allowed to overlap. Assuming a one-dimensional lattice, the number of arrangement of blobs in the pore is 2nd whereas in the free solution (mobile phase), the same blobs can be arranged in 6nd different ways. Then K = 1/3nd and The change in entropy ΔS is11,13

- ΔS/k = - lnKporeN/gN (b/d)2 ≅ (Rg/d)2

References, Notes & Further Readings
  1. E.F. Casassa, J. Polym. Sci. Part B: Polym. Lett. 5, 773 (1967)
  2. E.F. Casassa and Y. Tagami, Macromolecules, 2, 14 (1969)
  3. E.F. Casassa and Y. Tagami, J. Polym. Sci. Part A-2: Polym. Phys. 6, 63 (1968)
  4. E.F. Casassa,d J. Polym. Chem. 75, 3929 (1971)
  5. E.F. Casassa, Sep. Sci. 6, 305 (1971)
  6. W.W. Yau, C.P. Malon, and S.W. Fleming, J. Polym. Sci. B, 6,803 (1968)
  7. W.W. Yau., J. Polym. Sci. A27, 483 (1969)
  8. M. Doi, Trans. Faraday Soc. I1,71, 1720 (1975)
  9. J.C. Giddings et al., J. Phys. Chem. 72, 4397 (1968)
  10. M. Doi, S.F. Edwards, The Theory of Polymer Dynamics, Oxford University Press, New York 198
  11. I. Teraoka, Polymer Solutions: An Introduction to Physical Properties, John Wiley (2002)
  12. H.G. Barth, Early Development of Size-Exclusion Chromatography, LCGC 31, 7, 550-558 (2013)
  13. J.C. Day, B. Alince and A.A. Robertson, Can. J. Chem., 56, 2951-2958 (1978)
  14. K.H. Altgelt and L. Segal, Gel Permeation Chromatography, Marcel Dekker, New York. NY. 1971
  15. J.V. Dawkins, Pure & Appi. Chem. 54, 281 (1982)
  16. N. Donkai, A. Nakazawa, H. Inagaki, Bull. Inst. Chem. Res. 48, 79, 1970
  17. D. Berek, J. Sep. Sci. 33, 315 (2010)
  18. M. Rubinstein and R. Colby, Polymer Physics, 1st Ed., Oxford University Press (2003)

December 12, 2021

  • Summary

    Size EXclusion Chromatography

    also known as gel filtration chromatography, is a high throughput analytical method that is a special mode of high-pressure liquid chromatography (HPLC).

  • SEC employs porous particles (stationary phase) and solvent (mobile phase) in a long column to separate molecules (solute) of different size in solution. The separation process is based on the continuous size-depending molecule partitioning between the stationary and mobile phase.

  • The solute continually and recurrently permeates in and out of the pores of the stationary phase as the solute bands travel along with the solvent down the separation column.

  • The time a macromolecule requires to travel through the column depends on its hydrodynamic size and shape relative to the size and shape of the pores of the stationary phase.

  • Smaller macromolecules are retained in the pores longer than larger macromolecules. As a result, molecules with a larger hydrodynamic volume elute from the column faster than smaller molecules.

Read More