12.3: Optimizing Chromatographic Separations

Now that we have defined the solute retention factor, selectivity, and column efficiency we are able to consider how they affect the resolution of two closely eluting peaks. Because the two peaks have similar retention times, it is reasonable to assume that their peak widths are nearly identical. Equation 12.1, therefore, becomes

\[R_\ce{AB} = \dfrac{t_\textrm{r,B} − t_\textrm{r,A}}{0.5(w_\ce{B} +w_\ce{A})} ≈ \dfrac{t_\textrm{r,B} − t_\textrm{r,A}}{0.5(2w_\ce{B})} = \dfrac{t_\textrm{r,B} −t_\textrm{r,A}}{w_\ce{B}}\tag{12.17}\]

where B is the later eluting of the two solutes.

Note

If the number of theoretical plates is the same for all solutes—not strictly true, but not a bad assumption—then from equation 12.15, the ratio t_r/w is a constant. If two solutes have similar retention times, then their peak widths must be similar.

Solving equation 12.15 for w_B and substituting into equation 12.17 leaves us with the following result.

\[R_\ce{AB}= \dfrac{\sqrt{N}}{4} × \dfrac{t_\textrm{r,B} − t_\textrm{r,A}}{t_\textrm{r,B}}\tag{12.18}\]

Rearranging equation 12.8 provides us with the following equations for the retention times of solutes A and B.

\[t_\textrm{r,A} = k_\ce{A}t_\ce{m} + t_\ce{m} \hspace{40px} t_\textrm{r,B} = k_\ce{B}t_\ce{m} + t_\ce{m}\]

After substituting these equations into equation 12.18 and simplifying, we have

\[R_\ce{AB} = \dfrac{\sqrt{N}}{4} × \dfrac{k_\ce{B} - k_\ce{A}}{1 + k_\ce{B}}\]

Finally, we can eliminate solute A’s retention factor by substituting in equation 12.9. After rearranging, we end up with the following equation for the resolution between the chromatographic peaks for solutes A and B.

\[R_\ce{AB}= \dfrac{\sqrt{N}}{4} × \dfrac{α-1}{α} × \dfrac{k_\ce{B}}{1 + k_\ce{B}}\tag{12.19}\]

In addition to resolution, another important factor in chromatography is the amount of time needed to elute a pair of solutes, which we can approximate using the retention time for solute B.

\[t_\textrm{r,B}= \dfrac{16R_\ce{AB}^2H}{u} × \left(\dfrac{α}{α – 1}\right)^2 × \dfrac{(1 + k_\ce{B})^3}{k_\ce{B}^2}\tag{12.20}\]

where u is the mobile phase’s velocity.

Equation 12.19 and equation 12.20 contain terms corresponding to column efficiency, selectivity, and the solute retention factor. We can vary these terms, more or less independently, to improve resolution and analysis time. The first term, which is a function of the number of theoretical plates (for equation 12.19) or the height of a theoretical plate (for equation 12.20), accounts for the effect of column efficiency. The second term is a function of α and accounts for the influence of column selectivity. Finally, the third term in both equations is a function of k_B and accounts for the effect of solute B’s retention factor. A discussion of how we can use these parameters to improve resolution is the subject of the remainder of this section.

12.3.1 Using the Retention factor to Optimize Resolution

One of the simplest ways to improve resolution is to adjust solute B’s retention factor. If all other terms in equation 12.19 remain constant, increasing k_B improves resolution. As shown by the green curve in Figure 12.15, however, the improvement is greatest if the initial value of k_B is small. Once k_B exceeds a value of approximately 10, a further increase produces only a marginal improvement in resolution. For example, if the original value of k_B is 1, increasing its value to 10 gives an 82% improvement in resolution; a further increase to 15 provides a net improvement in resolution of only 87.5%.

Figure 12.15: Effect of k_B on the resolution for a pair of solutes, R_AB, and the retention time for the later eluting solute, t_r,B. The y-axes display the resolution and retention time relative to their respective values when k_B is 1.00.

Any improvement in resolution by increasing the value of k_B generally comes at the cost of a longer analysis time. The red curve in Figure 12.15 shows the relative change in solute B’s retention time as a function of its retention factor. Note that the minimum retention time is for k_B = 2. Increasing k_B from 2 to 10, for example, approximately doubles solute B’s retention time.

Note

The relationship between retention factor and analysis time in Figure 12.15 works to our advantage if a separation produces an acceptable resolution with a large k_B. In this case we may be able to decrease k_B with little loss in resolution and with a significantly shorter analysis time.

To increase k_B without significantly changing the selectivity, α, any change to the chromatographic conditions must result in a general, nonselective increase in the retention factor for both solutes. In gas chromatography, we can accomplish this by decreasing the column’s temperature. Because a solute’s vapor pressure is smaller at lower temperatures, it spends more time in the stationary phase and takes longer to elute. In liquid chromatography, the easiest way to increase a solute’s retention factor is to use a mobile phase that is a weaker solvent. When the mobile phase has a lower solvent strength, solutes spend proportionally more time in the stationary phase and take longer to elute.

Adjusting the retention factor to improve the resolution between one pair of solutes may lead to unacceptably long retention times for other solutes. For example, suppose we need to analyze a four-component mixture with baseline resolution and with a run-time of less than 20 min. Our initial choice of conditions gives the chromatogram in Figure 12.16a. Although we successfully separate components 3 and 4 within 15 min, we fail to separate components 1 and 2. Adjusting the conditions to improve the resolution for the first two components by increasing k_B,2 provides a good separation of all four components, but the run-time is now too long (Figure 12.16b). This problem of finding a single set of acceptable operating conditions is known as the general elution problem.

One solution to the general elution problem is to make incremental adjustments to the retention factor over time. Thus, we choose our initial chromatographic conditions to optimize the resolution for early eluting solutes. As the separation progresses, we adjust the chromatographic conditions to decrease k_B—and, therefore, decrease the retention time—for each of the later eluting solutes (Figure 12.16c). In gas chromatography this is accomplished by temperature programming. The column’s initial temperature is selected such that the first solutes to elute are fully resolved. The temperature is then increased, either continuously or in steps, to bring off later eluting components with both an acceptable resolution and a reasonable analysis time. In liquid chromatography the same effect is obtained by increasing the solvent’s eluting strength. This is known as a gradient elution. We will have more to say about each of these in later sections of this chapter.

Figure 12.16: Example showing the general elution problem in chromatography. See text for details.

12.3.2 Using Selectivity to Optimize Resolution

A second approach to improving resolution is to adjust the selectivity, α. In fact, for α ≈ 1 it usually is not possible to improve resolution by adjusting the solute retention factor, k_B, or column efficiency, N. A change in α often has a more dramatic effect on resolution than a change in k_B. For example, changing α from 1.1 to 1.5, while holding all other terms constant, improves resolution by 267%. In gas chromatography, we usually adjust α by changing the stationary phase, and we usually change the composition of the mobile phase in liquid chromatography.

To change α we need to selectively adjust individual solute retention factors. Figure 12.17 shows one possible approach for the liquid chromatographic separation of a mixture of substituted benzoic acids. Because the components are weak acids, their retention times vary with the pH of the mobile phase, as shown in Figure 12.17a. The intersections of these curves show pH values where two solutes co-elute. For example, at a pH of 3.8 terephthalic acid and p-hydroxybenzoic acid elute as a single chromatographic peak.

Note

Let’s use benzoic acid, C₆H₅COOH, to explain why pH can affect a solute’s retention time. The separation uses an aqueous mobile phase and a nonpolar stationary phase. At lower pHs, benzoic acid is predominately in its weak acid form, C₆H₅COOH, and easily partitions into the nonpolar stationary phase. At more basic pHs, however, benzoic acid is in its weak base form, C₆H₅COO^–. Because it now carries a charge, its solubility in the mobile phase increases and its solubility in the nonpolar stationary phase decreases. As a result, it spends more time in the mobile phase and has a shorter retention time.

Figure 12.17a shows that there are many pH values where some separation is possible. To find the optimum separation, we plot α for each pair of solutes. The red, green, and orange curves in Figure 12.17b show the variation in a with pH for the three pairs of solutes that are hardest to separate (for all other pairs of solutes, α > 2 at all pH levels). The blue shading shows windows of pH values in which at least a partial separation is possible—this figure is sometimes called a window diagram—and the highest point in each window gives the optimum pH within that range. The best overall separation is the highest point in any window, which, for this example, is a pH of 3.5. Because the analysis time at this pH is more than 40 min (see Figure 12.17a), choosing a pH between 4.1–4.4 might produce an acceptable separation with a much shorter analysis time.

Figure 12.17: Example showing how the mobile phase pH in liquid chromatography affects selectivity: (a) retention times for four substituted benzoic acids as a function of the mobile phase’s pH; (b) alpha values for three pairs of solutes that are difficult to separate. See text for details. The mobile phase is an acetic acid/sodium acetate buffer and the stationary phase is a nonpolar hydrocarbon. Data from Harvey, D. T.; Byerly, S.; Bowman, A.; Tomlin, J. “Optimization of HPLC and GC Separations Using Response Surfaces,” J. Chem. Educ. 1991, 68, 162–168.

Table 12.1: Minimum Number of Theoretical Plates to Achieve Desired Resolution for Selected Values of k_B and α

	RAB = 1.00		RAB = 1.25			RAB = 1.50
kB	α = 1.05	α = 1.10	α = 1.05		α = 1.10	α = 1.05		α = 1.10
0.5	63500	17400	99200		27200	143000		39200
1.0	28200	7740	44100		12100	63500		17400
1.5	19600	5380	30600		8400	44100		12100
2.0	15900	4360	24800		6810	35700		9800
3.0	12500	3440	19600		5380	28200		7740
5.0	10200	2790	15900		4360	22900		6270
10.0	8540	2340	13300		3660	19200		5270

12.3.3 Using Column Efficiency to Optimize Resolution

A third approach to improving resolution is to adjust the column’s efficiency by increasing the number of theoretical plates, N. If we have values for k_B and α, then we can use equation 12.19 to calculate the number of theoretical plates for any resolution. Table 12.1 provides some representative values. For example, if α = 1.05 and k_B = 2.0, a resolution of 1.25 requires approximately 24 800 theoretical plates. If our column provides only 12 400 plates, half of what is needed, then a separation is not possible. How can we double the number of theoretical plates? The easiest way is to double the length of the column, although this also doubles the analysis time. A better approach is to cut the height of a theoretical plate, H, in half, providing the desired resolution without changing the analysis time. Even better, if we can decrease H by more than 50%, it may be possible to achieve the desired resolution with an even shorter analysis time by decreasing k_B or α.

To decrease the height of a theoretical plate we need to understand the experimental factors that affect band broadening. There are several theoretical treatments of band broadening. We will consider one approach that considers four contributions: variations in paths length, longitudinal diffusion, mass transfer in the stationary phase, and mass transfer in the mobile phase.

Multiple Paths: Variations in Path Length

As solute molecules pass through the column they travel paths that differ in length. Because of this difference in path length, solute molecules entering the column at the same time, exit the column at different times. The result, as shown in Figure 12.18, is band broadening. The contribution of multiple paths to the height of a theoretical plate, H_p, is

\[H_\ce{p} = 2λd_\ce{p}\tag{12.21}\]

where d_p is the average diameter of the particulate packing material, and λ is a constant that accounts for the consistency of the packing. A smaller range of particle sizes and a more consistent packing produce a smaller value for λ. For a column without packing material, H_p is zero and there is no contribution to band broadening from multiple paths.

Note

Inconsistent packing can produces channels through the column that allow some solute molecules to travel quickly through the column. It also can lead to pockets that can temporarily trap some solute molecules, slowing their progress through the column. A more uniform packing minimizes these two problems.

Figure 12.18: The effect of multiple paths on a solute’s band broadening. The solute’s initial band profile is rectangular. As this band travels through the column, individual solute molecules travel different paths, three of which are shown by the meandering colored paths (the actual lengths of these paths are shown by the straight arrows at the bottom of the figure). Most solute molecules travel paths with lengths similar to that shown in blue, with a few traveling much shorter paths (green) or much longer paths (red). As a result, the solute’s band profile at the end of the column is broader and Gaussian in shape.

Longitudinal Diffusion

The second contribution to band broadening is the result of the solute’s longitudinal diffusion in the mobile phase. Solute molecules are constantly in motion, diffusing from regions of higher solute concentration to regions where the concentration of solute is smaller. The result is an increase in the solute’s band width (Figure 12.19). The contribution of longitudinal diffusion to the height of a theoretical plate, H_d, is

\[H_\ce{d} = \dfrac{2γD_\ce{m}}{u}\tag{12.22}\]

where D_m is the solute’s diffusion coefficient in the mobile phase, u is the mobile phase velocity, and γ is a constant related to the efficiency of column packing. Note that the effect of H_d on band broadening is inversely proportional to the mobile phase velocity—a higher velocity provides less time for longitudinal diffusion. Because a solute’s diffusion coefficient is larger in the gas phase than in a liquid phase, longitudinal diffusion is a more serious problem in gas chromatography.

Figure 12.19: The effect of longitudinal diffusion on a solute’s band broadening. Two horizontal crosssections through the column and the corresponding concentration versus distance profiles are shown, with (a) being earlier in time. The red arrow shows the direction in which the mobile phase is moving.

Mass Transfer

As the solute passes through the column it moves between the mobile phase and the stationary phase. We call this movement between phases mass transfer. As shown in Figure 12.20, band broadening occurs if the solute’s movement within the mobile phase or within the stationary phase is not fast enough to maintain an equilibrium partitioning of solute between the two phases. On average, solute molecules in the mobile phase move further down the column before passing into the stationary phase. Solute molecules in the stationary phase, on the other hand, take longer than expected to move back into the mobile phase. The contributions of mass transfer in the stationary phase, H_s, and mass transfer in the mobile phase, H_m, are given by the following equations

\[H_\ce{s} = \dfrac{qkd_\ce{f}^2}{(1 + k)^2D_\ce{s}}u\tag{12.23}\]

\[H_\ce{m} = \dfrac{fn(d_\ce{p}^2,d_\ce{c}^2)}{D_\ce{m}}u\tag{12.24}\]

where d_f is the thickness of the stationary phase, d_c is the column’s diameter, D_s and D_m are the solute’s diffusion coefficient in the stationary phase and the mobile phase, k is the solute’s retention factor, and q is a constant related to the column packing material. Although the exact form of H_m is not known, it is a function of particle size and column diameter. Note that the effect of H_s and H_m on band broadening is directly proportional to the mobile phase velocity—smaller velocities provide more time for mass transfer.

Note

The abbreviation fn in equation 12.24 means “function of.”

Figure 12.20 Effect of mass transfer on band broadening: (a) Ideal equilibrium Gaussian profiles for the solute in the mobile phase and in the stationary phase. (b, c) If we allow the solute’s band to move a small distance down the column, an equilibrium between the two phases no longer exits. The red arrows show the movement of solute—what we call the mass transfer of solute—from the stationary phase to the mobile phase, and from the mobile phase to the stationary phase. (d) Once equilibrium is reestablished, the solute’s band is now broader.

Putting it All Together

The height of a theoretical plate is a summation of the contributions from each of the terms affecting band broadening.

\[H = H_\ce{p} + H_\ce{d} + H_\ce{s} + H_\ce{m}\tag{12.25}\]

An alternative form of this equation is the van Deemter equation

\[H = A + \dfrac{B}{u}+ Cu\tag{12.26}\]

which emphasizes the importance of the mobile phase’s velocity. In the van Deemter equation,

• The \(A\) term accounts for the contribution of multiple paths (H_p),
• The \(B/u\) term accounts for the contribution of longitudinal diffusion (H_d), and
• The \(Cu\) term accounts for the combined contribution of mass transfer in the stationary phase and mobile phase (H_s and H_m).

There is some disagreement on the best equation for describing the relationship between plate height and mobile phase velocity.⁵ In addition to the van Deemter equation, other equations include

\[H = \dfrac{B}{u}+ (C_\ce{s} + C_\ce{m})u\]

where C_s and C_m are the mass transfer terms for the stationary phase and the mobile phase and

\[H = Au^{1/3} + \dfrac{B}{u} + Cu\]

All three equations, and others, have been used to characterize chromatographic systems, with no single equation providing the best explanation in every case.⁶

To increase the number of theoretical plates without increasing the length of the column we need to decrease one or more of the terms in equation 12.25. The easiest way to decrease H is by adjusting the velocity of the mobile phase. For smaller mobile phase velocities, column efficiency is limited by longitudinal diffusion, and at higher velocities efficiency is limited by the two mass transfer terms. As shown in Figure 12.21—which uses the van Deemter equation—the optimum mobile phase velocity is the minimum in a plot of H as a function of u.

Figure 12.21: Plot showing the relationship between the height of a theoretical plate, H, and the mobile phase’s velocity, u, based on the van Deemter equation.

The remaining parameters affecting the terms in equation 12.25 are functions of the column’s properties and suggest other possible approaches to improving column efficiency. For example, both H_p and H_m are a function of the size of the particles used for packing the column. Decreasing particle size, therefore, is another useful method for improving efficiency.

Note

The smaller the particles, the more pressure is needed to push the mobile phase through the column. As a result, for any form of chromatography there is a practical limit to particle size.

Perhaps the most important advancement in chromatography columns is the development of open-tubular, or capillary columns. These columns have very small diameters (d_c ≈ 50–500 μm) and contain no packing material (d_p = 0). Instead, the interior wall of a capillary column is coated with a thin film of the stationary phase. Plate height is reduced because the contribution to H from H_p (equation 12.21) disappears and the contribution from H_m (equation 12.24) becomes smaller. Because the column does not contain any solid packing material, it takes less pressure to move the mobile phase through the column, allowing for much longer columns. The combination of a longer column and a smaller height for a theoretical plate increases the number of theoretical plates by approximately 100~. Capillary columns are not without disadvantages. Because they are much narrower than packed columns, they require a significantly smaller amount of sample, which may be difficult to inject reproducibly. Another approach to improving resolution is to use thin films of stationary phase, which decreases the contribution to H from H_s (equation 12.23).

Contributors

David Harvey (DePauw University)