The interpretation of results from chromatographic measurements can often be augmented with an appropriate mathematical treatment of the solute retention that is observed. The goal of the treatment is to make the resulting metric as independent of the instrument as possible. A typical situation that arises from analysis by gas chromatography with mass spectrometry is that the library search routine produces “hits” that are ambiguous if not nonsensical (Ref. 1). The correct interpretation of the mass spectrum must then be done manually (Refs. 2, 3), and the mass spectral data should be augmented by additional analytical techniques. The specific techniques that should be used must be determined on a case-by-case basis by a qualified person. One additional datum that is typically already present in gas chromatography with mass spectrometry detection is chromatographic retention. The raw datum from a chromatographic measurement is a retention time, tr, of each eluted peak, and a corresponding intensity. Here, we will not treat other aspects of the output, such as the width and shape of the chromatographic signal. The retention time (for a given stationary phase) is dependent on the column temperature, column pressure, column geometry (length and inside diameter; phase ratio), and ambient (atmospheric) pressure.
If the volumetric carrier gas flow rate (at the column exit) is measured and multiplied by the retention time, the retention volume, VR, is obtained. The adjusted retention volume, VR′, is the retention volume corrected for the void volume (or mobile phase holdup) of the column. It is obtained by simply subtracting the retention volume of an unretained solute (VM):
VR′ = VR − VM (1)
While it is possible to calculate the corresponding adjusted retention time, tr′, by subtracting the retention time of the unretained peak, tm, it is better to work with volumes because average flow-rate variations between individual analyses are then accounted for. Note that with each level of refinement beyond the raw retention time, a facet of instrument dependence from the resulting parameter is removed. VR is independent of the flow rate; VR′ is, further, independent of column geometry. Continuing with this approach, the net retention volume, VN, is defined, by applying a factor, j, to account for the pressure drop across the column:
VN = jVR′ (2)
where j is usually the Martin-James compressibility factor3:
(3)
In Eq. 3, Pi is the inlet pressure (absolute) and Po is the outlet pressure (usually atmospheric pressure). The net retention volume is important because it is independent of the inlet and outlet pressures, as well as being independent of the carrier gas flow rate and column geometry.
The specific retention volume, Vg, corrects the net retention volume for the amount of stationary phase actually on the column, and the column temperature is adjusted or corrected to 0 °C:
(4)
where Tcol is the column temperature (in K), and Ws is the mass of the stationary phase in the column. The Vg value is a characteristic for a particular solute on a particular stationary phase and is instrument independent. This is a quantity that may be compared from instrument to instrument and laboratory to laboratory with a high level of confidence provided the stationary phase used is a single, pure compound or a well characterized mixture. If the mass of the stationary phase is not known, or is not meaningful, one may use the net retention volume directly, or one may correct the net retention volume to a column temperature of 0 °C (represented by VN0) by simply not including the term for Ws (that is, setting it equal to unity).
It is also extremely valuable to calculate a relative retention, ra/b:
(5)
where the numerical superscripts refer to the retention volumes of solutes “a” and “b.” In this case, solute “a” is a reference compound. The relative retention is dependent only on the column temperature and the type of stationary phase. For reasons of operational simplicity, this parameter is usually one of the best to use for qualitative analysis. It can account for small differences in the column temperature, stationary phase considerations, column history, and minor disturbances in the carrier gas flow rate. It is possible to account for the column temperature by plotting the logarithm of the retention parameters against 1/T, where T is the thermodynamic temperature. The column pressure is accounted for by variations in the volume measurement; therefore, there is no pressure dependence to these parameters.
We can go beyond the simple retention parameters discussed earlier to incorporate a logarithmic interpolation on a uniform scale by use of the Kovats retention index (Ref. 4). The isothermal Kovats retention index is calculated by use of the following defining equation:
(6)
Here, Isample is the dimensionless Kovats retention index that is a function of both temperature and the stationary phase employed. The terms represented by X are retention parameters of the sample and standards. Following this convention, XS is the retention parameter of the sample under consideration. Any retention parameter, such as the adjusted retention time, t′, the net retention volume, VN, the adjusted net retention volume, VN0, and the relative retentions, ra/b, can be used. XS is the retention parameter of the sample under consideration, XL is the retention parameter of a normal alkane (that is, straight chain or unbranched) of carbon number nL that elutes earlier than the sample, and XH is the retention parameter of a normal alkane having a carbon number greater than nL+1 that elutes after the sample. The retention index of a sample is, therefore, 100 multiplied by the carbon number of a hypothetical normal alkane that shows the same retention parameter on the stationary phase at that temperature. Thus, a sample that has a retention index of 785, for example, would co-elute with a hypothetical normal alkane that has 7.85 carbon atoms. By definition, the retention indices of the normal alkanes (on any stationary phase) are equal to 100 multiplied by the carbon number. Thus, for n-hexane, I = 600, and similarly for the other n-alkanes in this homologous series. The zero point in the scale is defined for hydrogen, for which I = 0. Kovats retention indices have been determined for selected compounds on the most common stationary phases (Refs. 5-7).
The temperature dependence of I for a given sample is known to follow a hyperbolic form similar to the familiar Antoine equation used to represent vapor pressure:
(7)
In this equation, A, B, and C are empirically determined constants, and T is the temperature (in °C). A nonlinear fitting routine should be used to determine these constants. When retention indices are available for at least three temperatures, initial values for A, B, and C can be determined by use of the following equations:
(8)
(9)
B = (I2 − A)(T2 + C) (10)
In these equations, I1, I2, and I3 are retention indices of the sample measured at temperatures T1, T2, and T3. When additional retention indices are available at other temperatures, we advocate the use of minimum deviation estimates from these three equations to furnish the starting values for the nonlinear fit. When retention indices at four temperatures are available, the best starting values are obtained from the I1, I2, and I3 triplet that minimizes the deviation with the experimental value with that produced by Eq. 7. This approach provides the fastest convergence, and also helps avoid converging to local minima. Predictions made by use of Eq. 7 can be used for retention indices within the measured temperature range as well as extrapolation somewhat beyond that range on a case-by-case basis.
It is also of value to report and use the temperature dependence as a slope coefficient, δIsample/10 °C, the variation of Isample for a particular stationary phase over a particular temperature range. While not as reliable as the Antoine-type fit, this coefficient is useful for predictions within the range of the measured results.
In all of the above discussion, one must understand that the column temperature is fixed. The Kovats retention indices can be made applicable to temperature programmed analysis by use of
(11)
where I is the retention index, n is the number of carbon atoms in the smaller n-alkane, N is the number of carbon atoms in the larger n-alkane, and tr is the retention time of the indicated peak.
A useful alternative to the Kovats system is the Lee retention index (isothermal and temperature dependent), based on the polynuclear aromatic hydrocarbon (PAH) standard compounds: naphthalene (I = 200), phenanthrene (I = 300), chrysene (I = 400), and picene or benzo(g,h,i)perylene (I = 500). Isothermal (and temperature-dependent) Kovats and Lee retention indices for many compounds are tabulated in the NIST Chemistry Web Book (Ref. 8).