Simulation of the effects of negative thermal gradients on separation performance of gas chromatography

End Column Reverse Chromatography as a Novel Approach for Enhanced Separation: A Pilot Study

BACKGROUND

Currently, the most popular technique in gas chromatography (GC) is "temperature programming," where the temperature increases from the start of the injection. This leads to faster elution of analytes compared to isothermal methods. However, isothermal methods are considered optimal for separating compounds with similar retention times. Another interesting technique that provides higher resolution is dynamic thermal gradient gas chromatography (TGGC), where separations are achieved as a decreasing thermal gradient. This gradually decreases the positive gas velocity. Nevertheless, it was proven that GC techniques with negative velocity gradients do not improve the resolution of compounds with nearly identical retention times.

OBJECTIVE

Optimizing a new GC approach to combine both the short time from positive temperature ramps programming, and the enhanced separation of the negative ramps of the TGGC, a model under the name of "end column reverse chromatography" (ECRC).

METHODS

The process simply consists of two steps: the first is a normal positive ramp from the start of the injection, and the second step is a negative thermal ramp at a time that is around the retention time of the first eluting peak. This will decrease the solute velocity almost solely for the second compound, leading to relatively enhanced separation.

RESULTS

The optimized ECRC method increased the resolution of two isomers (trans- and cis-chlordane) from 1 (slightly overlapping) in the case of temperature programming to 2.78 as shown in this study. This comes at the expense of the width and intensity of the peaks, where the intensity decreased about 17 and 12% for cis- and trans-chlordane, and the peak width increased with 37 and 77% for the same compounds, respectively.

CONCLUSIONS

ECRC is a novel model for enhanced separation that comes with some drawbacks.

HIGHLIGHTS

It can be an alternative approach to get a fast GC method with enhanced separation for isomers.

Relation between characteristic temperature and elution temperature in temperature programmed gas chromatography - Part II: Influence of column properties

The method development process in gas chromatography can be accelerated by suitable computer simulation tools using knowledge about the solute-column interactions described by thermodynamic retention parameters. Since retention parameters usually are determined under isothermal conditions, the presented work offers a step to estimate one of the most important retention parameters, the characteristic temperature Tchar by less laborious temperature programmed measurements. In the first part an empirical multivariate model was introduced describing the correlation between the elution temperature Telu of a solute and its characteristic temperature Tchar. Now in the second part a simulation model of GC and available retention data from a retention database was used to investigate the correlation between Telu and Tchar for an expanded range of heating rates and initial temperatures. In addition to part I, the simulation is used to investigate the influences of different properties of the separation column such as different phase ratios and column geometries like length and diameter or various stationary phases including SLB-5 ms, SPB-50, Stabilwax, Rtx-Dioxin2, Rxi-17Sil MS, Rxi-5Sil MS, ZB-PAH-CT, DB-5 ms, Rxi-5 ms, Rtx5 and FS5ms. The fit model is valid for all investigated stationary phases. The influence of the phase ratio to the correlation could be determined. Therefore, the model was expanded to this parameter. The expanded range of heating rates and the normalization for the system independent dimensionless heating rate required a further modification of the previously presented correlation model. The model now fits also under isothermal conditions. The results were used for estimation of the Tchar of an analyte from the elution temperature in the temperature program. The prediction performance was investigated and evaluated for 20 different temperature program conditions and at two phase ratios (β=125 and β=250). Under best conditions the estimated and the measured Tchar values show relative differences <0.5 %. With this novel model estimations for Tchar are possible at 20 °C above the initial temperature, which expands the prediction range even for low and medium retained analytes compared to earlier approaches.

Simulation of gas chromatographic separations and estimation of distribution-centric retention parameters using linear solvation energy relationships

For method development in gas chromatography, suitable computer simulations can be very helpful during the optimization process. For such computer simulations retention parameters are needed, that describe the interaction of the analytes with the stationary phase during the separation process. There are different approaches to describe such an interaction, e.g. thermodynamic models like Blumberg's distribution-centric 3-parameter model (K-centric model) or models using chemical properties like the Linear Solvation Energy Relationships (LSER). In this work LSER models for a Rxi-17Sil MS and a Rxi-5Sil MS GC column are developed for different temperatures. The influences of the temperature to the LSER system coefficients are shown in a range between 40 and 200 °C and can be described with Clark and Glew's ABC model as fit function. A thermodynamic interpretation of the system constants is given and its contribution to enthalpy and entropy is calculated. An estimation method for the retention parameters of the K-centric model via LSER models were presented. The predicted retention parameters for a selection of 172 various compounds, such as FAMEs, PCBs and PAHs are compared to isothermal determined values. 40 measurements of temperature programmed GC separations are compared to computer simulations using the differently determined or estimated K-centric retention parameters. The mean difference (RSME) between the measured and predicted retention time is less than 8 s for both stationary phases using the isothermal retention parameters. With the LSER predicted parameters the difference is 20 s for the Rxi-5Sil MS and 38 s for the Rxi-17Sil MS. Therefore, the presented estimation method can be recommended for first method development in gas chromatography.

Theory of linear focusing in chromatographic columns with exponential retention. Part 2: Analysis of special cases

The unified approach to studies of different separation techniques (GC, LC, etc.) adopted in Part 1 continued herein. As before, the column temperature in GC, the solvent composition in LC, etc., are represented by the concept of solute mobilization (y). General equations for dynamic (moving) linear y-gradients (gradients in y) in non-uniform columns developed in Part 1 are reduced herein to special cases of uniform columns. Only the uniform columns, the ideal sample introduction, the positive y-rates, the negative y-gradients, and the eluting solutes are considered. Equations for solute band compression, peak width, peak focusing, peak separation and others derived. All equations are expressed via dimensionless parameters eliminating unessential factors like the mobile phase type, column type and dimensions, etc. The conventional gradient elution LC is treated as a special case of the separations with y-gradients. Effects of the y-gradients on solute retention time, on band compression, on peak focusing and on peak separation are analyzed. The equations developed herein, while describing simple and familiar concepts (retention time, peak width, etc.), are, unfortunately, cumbersome in many cases. The good news is that the equations are exact and the conditions of existence of the equations are explicitly formulated so that there is no question where the equations can and where thy cannot be used.

Estimation of retention parameters from temperature programmed gas chromatography

A fast and reliable method is presented to evaluate retention parameters of the distribution-centric 3-parameter model from temperature programed gas chromatographic measurements. Based on a fully differentiable model of the migration of solutes in a gas chromatographic (GC) system, Newton's method with a trust region is used to determine the three parameters, respectively the three parameters and the column diameter, of several solutes as the minima of the difference between measured and calculated retention times. The determined retention parameters can then be used in method development, using the simulation of GC separation. The results of the retention parameters are compared to the parameters determined using isothermal GC measurements and show good agreement, with deviations of less than 0.5% (1.8 K) for the most important parameter of characteristic temperature T_char. Using the estimated retention parameters, additional GC separations are simulated and compared with measurements. Retention times in additional temperature programmed measurements could be predicted with less than 0.7% deviation. Four to five different temperature programs are enough to determine reliable retention parameters. Unless the column diameter and the column length are exactly known, it is preferable to also estimate the diameter (more precisely the L/d-ratio) together with the retention parameters.

Theory of linear focusing in chromatographic columns with exponential retention. Part 1: Basic solutions

This report is the first of 2-part study of the effect of gradients in column parameters on the column performance. If t, x and p are, respectively, time since sample introduction, distance from column inlet and some parameter of solute migration along the column then ∂p/∂t and ∂p/∂x are, respectively, the rate of changing p and the gradient of p. Unified approach to study of gradients and rates in different chromatographic techniques (LC, GC, etc.) has been developed. To facilitate a unified approach, the umbrella term mobilization (y) representing column temperature (T) in GC, solvent composition (ϕ) in LC, etc. is introduced. Differential equations for migration of a solute band (collection of solute molecules) under the following conditions are formulated and solved:The key solutions describe the time of migration of a solute band and the band width - both as functions of the distance traveled by the band. The solutions are used in Part 2 for the study of the effects of the negative gradients in y on column performance in several practically important cases. A reduction of the key general solutions to much simpler equations for gradient LC has been demonstrated herein as an example.

Comparison of the Dynamic Thermal Gradient to Temperature-Programmed Conditions in Gas Chromatography Using a Stochastic Transport Model

This paper compares dynamic (i.e., temporally changing) thermal gradient gas chromatography (GC) to temperature-programmed GC using a previously published stochastic transport model to simulate peak characteristics for the separation of C12-C40 hydrocarbons. All comparisons are made using chromatographic conditions that give approximately equal analyte retention times (t_R). As shown previously, a static thermal gradient does not improve resolution (R_s) equally for all analytes, which highlights the need for a dynamic thermal gradient. An optimal dynamic thermal gradient should result in constant analyte velocities at any instant in time for those analytes that are actively being separated (i.e., analytes that have low retention factors). The average separation temperature for each analyte is used to determine the thermal gradient profile at different times in the temperature ramp. Because many of the analytes require a similar thermal gradient profile when actively being separated, the thermal gradient profile in this study was held fixed; however, the temperature of the entire thermal gradient was raised over time. From the simulations performed in this study, optimized dynamic thermal gradient conditions are shown to improve R_s by up to 13% over comparative temperature-programmed conditions, even with a perfect injection (i.e., zero injection bandwidth). In the dynamic thermal gradient simulations, all analytes showed improvements in R_s along with slightly shorter t_R values compared to simulations for traditional temperature-programmed conditions.

Residual solvent analysis with hyper-fast gas chromatography-mass spectrometry and a liquid carbon dioxide cryofocusing in less than 90 s

A new hyper-fast gas chromatography method with less than 90 s runtime including the column cool down was developed for the analysis of four gases and 16 residual solvents, combining a CO₂ cryofocusing with a flow-field thermal gradient gas chromatograph (FF-TG-GC) and ToF-MS. The extremely low analysis time can be achieved by combining the new FF-TG-GC and a very short Rxi-624 Sil MS separation column with a small inner diameter and small film thickness (2.05 m × 0.1 mm × 1.0 µm). The column is inserted into a low thermal mass, resistively heated stainless steel capillary. This enables fast temperature programs with heating rates up to 3000 °C/min and a column cool down within a few seconds. In addition to temporal temperature gradients, the FF-TG-GC can generate a spatial temperature gradient that leads to an improved peak shape. Further, an external liquid CO₂ cryo-trap was designed in order to reduce the injection bandwidths of analytes and to take full advantage of the resolving power of the separation column. No modifications are required to the FF-TG-GC for the use of the cryogenic trap, as the cooled spot is heated by the resistively heated stainless steel capillary during the temperature program. With cryofocusing, analyzed residual solvents are baseline separated. R² values over 0.99 for calibration curves and low relative standard deviations (mainly < 3%) for repeatability tests were obtained.