Twenty-fifth anniversary of the retention index system in gas—liquid chromatography

Comprehensive study of retention influencing gas chromatographic parameters affecting linear retention indices

Retention in gas chromatographic systems has a central role in the identification of compounds even if detectors providing spectral information are used. But linear retention indices (LRI) of a single compound originating from multiple sources tend to vary greatly, probably due to differences in the experimental settings of the determinations. The effect of gas chromatographic parameters on LRI has been investigated using 41 compounds - previously identified from food contact plastics - and n-alkanes (n-C7-n-C40) used as reference series. As the reproducibility of LRIs under the same conditions is generally very good, the smallest changes in the settings often caused statistically significant, though irrelevant changes in the LRI values. Therefore, a multicriterial scoring-ranking system has been worked out to highlight the LRI value differences. Our results highlight that column length, heating rate, and film thickness can all be the reasons of the varying published LRI values. We also demonstrated that for the reproduction of LRI data, the chemistry (and not simply the polarity) of the stationary phase is crucial.

AIRI: Predicting Retention Indices and Their Uncertainties Using Artificial Intelligence

The Kováts retention index (RI) is a quantity measured using gas chromatography and is commonly used in the identification of chemical structures. Creating libraries of observed RI values is a laborious task, so we explore the use of a deep neural network for predicting RI values from structure for standard semipolar columns. This network generated predictions with a mean absolute error of 15.1 and, in a quantification of the tail of the error distribution, a 95th percentile absolute error of 46.5. Because of the Artificial Intelligence Retention Indices (AIRI) network's accuracy, it was used to predict RI values for the NIST EI-MS spectral libraries. These RI values are used to improve chemical identification methods and the quality of the library. Estimating uncertainty is an important practical need when using prediction models. To quantify the uncertainty of our network for each individual prediction, we used the outputs of an ensemble of 8 networks to calculate a predicted standard deviation for each RI value prediction. This predicted standard deviation was corrected to follow the error between the observed and predicted RI values. The Z scores using these predicted standard deviations had a standard deviation of 1.52 and a 95th percentile absolute Z score corresponding to a mean RI value of 42.6.

Nontargeted Screening for the Verification of Allergenic Ingredients and Perfume Authenticity by GC-ecTOF-MS

Fraudulent products are ubiquitous in all markets. Besides the financial aspect, a major issue regarding products adulteration for society and the environment is the missing regulation and control of these goods. Therefore, harmful and toxic compounds in fraudulent products may become a risk for human health and the environment. Even small amounts of toxic substances can still have damaging effects. Thus, a sensitive and reliable identification of compounds is needed. Novel technologies are necessary to use identifying and preventing fraud and related risks. Goods from the food, flavor, and fragrance markets often contain volatile organic compounds (VOCs), which include most allergenic fragrances. For the detection and identification of these substances, gas chromatographic separation hyphenated with high resolution mass spectrometry (GC–HRMS) is an ideal instrumental technique. GC–simultaneous electron and chemical ionization (ec) time-of-flight (TOF)-MS generates various types of information via simultaneous ec–HRMS. Advantages are given for target, known unknown, and unknown unknown data analysis by generating various types of ions within one single experimental GC–MS run. In this study, the experimental nontargeted screening approach and corresponding data analysis workflows—simultaneously using molecular ion information and structural information—are presented for the molecular identification and authenticity verification process from a brand perfume using GC–ecTOF-MS.

Development and validation of analytical method for identification of new psychoactive substances using linear retention indexes and gas chromatography-mass spectrometry

New Psychoactive Substances (NPS) are quickly developing to evade legislation, posing unprecedented challenges to public health and law enforcement authorities around the world. The aim of this work was to develop and validate a simple and reliable non-target gas chromatography/mass spectrometry (GC/MS) analytical method based on linear retention indexes for the expeditious identification of NPS without the need of analytical standards. The method was optimized and validated for 22 different drugs covering ten categories: phenethylamines (amphetamine, MDMA, methamphetamine, 25CNBOMe, 2-FA, 5-MAPB), "classic" drugs (cocaine, ephedrine, THC, heroine), synthetic cannabinoids (JWH-081, AM-2201, JWH-210, MAM-2201), piperazines (o-CPP, p-CPP), tryptamines (5-MeO-MiPT), synthetic cathinones (N-ethylpentylone), synthetic opioids (U-47700), aminoindanes (5-IAI), plant-based substances (Salvinorin-A) and "other" (methiopropamine). Three figures of merit (Selectivity, Precision and Robustness) were evaluated with retention index confidence intervals ranging from 0.5 to 20.6 i.u. and relative standard deviations in the range of 0.003% to 0.027% (repeatability) and 0.02% to 0.29% (intermediate precision). A general equation for estimating linear retention index variation as a function of retention time tolerance has been derived. This result in combination with a 2_III^6-3 fractional factorial design allowed to conclude column polarity to be only statistically relevant factor as compared to gas flow, split ratio, injection temperature, temperature program offset and column brand.

Phenomenon of dual- and single-retention behaviors of solutes and its validation by computational simulation in linear programmed temperature gas chromatography

The current theory of programmed temperature gas chromatography considers that solutes are focused by the stationary phase at the column head completely and does not explicitly recognize the different effects of initial temperature (To ) and heating rate (rT ) on the retention time or temperature of a homologue series. In the present study, n-alkanes, 1-alkenes, 1-alkyl alcohols, alkyl benzenes, and fatty acid methyl esters standards were used as model chemicals and were separated on two nonpolar columns, one moderately polar column and one polar column. Effects of To and rT on the retention of nonstationary phase focusing solutes can be explicitly described with isothermal and cubic equation models, respectively. When the solutes were in the stationary phase focusing status, the single-retention behavior of solutes was observed. It is simple, dependent upon rT only and can be well described by the cubic equation model that was visualized through four sequential slope analyses. These observed dual- and single-retention behaviors of solutes were validated by various experimental data, physical properties, and computational simulation.

A new accurate quadratic equation model for isothermal gas chromatography and its comparison with the linear model

The gas holdup time (tM) is a dominant parameter in gas chromatographic retention models. The difference equation (DE) model proposed by Wu et al. (J. Chromatogr. A 2012, http://dx.doi.org/10.1016/j.chroma.2012.07.077) excluded t(M). In the present paper, we propose that the relationship between the adjusted retention time t'RZ and carbon number z of n-alkanes follows a quadratic equation (QE) when an accurate tM is obtained. This QE model is the same as or better than the DE model for an accurate expression of the retention behavior of n-alkanes and model applications. The QE model covers a larger range of n-alkanes with better curve fittings than the linear model. The accuracy of the QE model was approximately 2-6 times better than the DE model and 18-540 times better than the LE model. Standard deviations of the QE model were approximately 2-3 times smaller than those of the DE model.

Difference equation model for isothermal gas chromatography expresses retention behavior of homologues of n-alkanes excluding the influence of holdup time

A difference equation (DE) model is developed using the methylene retention increment (Δtz) of n-alkanes to avoid the influence of gas holdup time (tM). The effects of the equation orders (1st-5th) on the accuracy of a curve fitting show that a linear equation (LE) is less satisfactory and it is not necessary to use a complicated cubic or higher order equation. The relationship between the logarithm of Δtz and the carbon number (z) of the n-alkanes under isothermal conditions closely follows the quadratic equation for C3-C30n-alkanes at column temperatures of 24-260 °C. The first and second order forward differences of the expression (Δlog Δtz and Δ2log Δtz, respectively) are linear and constant, respectively, which validates the DE model. This DE model lays a necessary foundation for further developing a retention model to accurately describe the relationship between the adjusted retention time and z of n-alkanes.

Isothermal retention indices on poly(3,3,3-trifluoropropylmethylsiloxane) stationary phases

Isothermal retention indices (I) of 56 solutes of varied chemical nature and polarity were determined by a well-proven accurate method on poly(3,3,3-trifluoropropylmethylsiloxanes) with trifluoropropyl group percentages between 0 and 50% within the 353-413 K temperature range. The column temperature effect on the retention index for a representative set of solutes was studied. The linear model and the recent expression I = a + bT(-1) + c ln T were discussed. The retention index dependence on the polymer polarity was also investigated. Principal component analysis and multiple linear regression (MLR) were applied to obtain retention models relating the I values to physicochemical and molecular properties of the solutes chromatographed. The best model of prediction of the retention in these fluorinated stationary phases was obtained by using the Abraham's descriptors.

Development of a database of gas chromatographic retention properties of organic compounds

A comprehensive database of gas chromatographic retention properties of chemical compounds has been developed using multiple literature sources. The National Institute of Standards and Technology (NIST) database of retention data for non-polar and polar stationary phases currently contains 292,924 data records for 42,888 compounds. The database includes data for Kováts indices, linear indices, Lee indices, retention times and retention volumes. The first release of this database for non-polar stationary phases is available with NIST/US Environmental Protection Agency (EPA)/National Institutes of Health (NIH) Mass Spectral Database (June 2005) and through the internet (NIST Chemistry WebBook). The paper describes the database and the process by which it has been compiled. The format of data presentation and the quality control procedures are described. Data sources of gas chromatographic retention data are also discussed.

Estimation of Kováts retention indices using group contributions

We have constructed a group contribution method for estimating Kováts retention indices by using observed data from a set of diverse organic compounds. Our database contains observed retention indices for over 35,000 different molecules. These were measured on capillary or packed columns with polar and nonpolar (or slightly polar) stationary phases under isothermal or nonisothermal conditions. We neglected any dependence of index values on these factors by averaging observations. Using 84 groups, we determined two sets of increment values, one for nonpolar and the other for polar column data. For nonpolar column data, the median absolute prediction error was 46 (3.2%). For data on polar columns, the median absolute error was 65 (3.9%). While accuracy is insufficient for identification based solely on retention, it is suitable for the rejection of certain classes of false identifications made by gas chromatography/mass spectrometry.

Conversion of programmed-temperature retention indices from one set of conditions to another

In order to make programmed-temperature retention index (PTRI) data be shared by other chromatographers and laboratories, conversion of PTRI from one set of experimental conditions to another is investigated in detail in this work. It was found that the differences between the PTRIs at different heating rates are structurally dependent, especially the number of ring in molecules. Thus, with the help of molecule constitutional descriptors, equations of PTRI conversion to certain initial temperature, heating rate, and stationary phase were obtained with high correlation coefficients and low standard deviations. Calculation errors of PTRI conversion between different heating rates and between different initial temperatures were from 1.1 to 2.9 retention index units (i.u.), which is in the same order with experiment errors. It is well known that reproducibility of PTRI on a polar column is not as good as that on an apolar column because of the apolarity of the n-alkane homologues. Thus, topological descriptors were used for PTRI conversion between two columns with different polar stationary phases, giving better results than those obtained by constitutional descriptors. This shows that topological descriptors could provide more molecular structural information than constitutional descriptors. However, as constitutional descriptor has the advantages of clear physical meaning and very simple calculation, it is our first selection when the PTRI calculation accuracy is satisfied. The method developed is simple in calculation, easy to be performed with high accuracy.

Retention in gas–liquid chromatography with a polyethylene oxide stationary phase: Molecular simulation and experiment

Configurational-bias Monte Carlo simulations in the isobaric-isothermal Gibbs ensemble were carried out to investigate the partitioning of normal alkanes, primary and secondary alcohols, symmetric alkyl ethers and arenes between a helium vapor phase and a polyethylene oxide stationary phase (M(W)=382 g mol(-1)). The united-atom version of the transferable potentials for phase equilibria force field was used to model all solutes, polyethylene oxide and helium. The Gibbs free energies of transfer and Kovats retention indices of the solutes were calculated directly from the partition constants at two different temperatures, 353 and 393 K. Chromatographic experiments on a Carbowax 20M retentive phase were performed for the same set of solutes and temperatures ranging from 333 to 413 K. The predicted retention indices for alcohols, ethers and arenes are overestimated by about 120, 70 and 20 retention index units, respectively, pointing to an overestimation of the first-order electrostatic interactions in the model system. Molecular-level analysis shows that hydrogen-bonding and dipole-dipole interactions lead to orientational ordering for the alcohol and ether analytes, whereas the weaker dipole-quadrupole interactions for the arene solutes are not sufficient to induce orientational ordering. The retention indices of alcohols and ethers decrease with increasing temperature because of the large entropic cost of hydrogen-bonding and orientational ordering. In contrast, the retention indices for arenes increase with increasing temperature because the entropic cost of cavity formation is smaller for arenes than for comparable alkanes.

Classification and comparison of capillary columns by determination of the solution enthalpy of polar and non polar probes

The behaviour and separation characteristics of capillary columns containing different stationary phases (bonded methyl- and methylphenylsiloxanes and polyglycols, and carbon layer and porous polymer) were classified and compared by measuring the values of the solution enthalpy DeltaHs of n-alkanes and 1-alcohols at various temperatures. The difference in DeltaHs values between straight-chain 1-alcohols and n-alkanes with the same number of carbon atom (DeltaDeltaHs) does not change with changing temperature when gas-liquid partition stationary phases are used, but depends on the column composition and polarity and therefore permits to evaluate the influence of polarity on the solute-solvent interaction. The trend of the DeltaDeltaHs values is correlated with that of the DeltaC parameter and with the McReynolds' polarity constants.

Temperature Dependence of Hydrogen Bonding: An Investigation of the Retention of Primary and Secondary Alcohols in Gas−Liquid Chromatography

Configurational-bias Monte Carlo simulations in the Gibbs Ensemble were carried out to investigate the analyte partitioning of n-pentane, n-hexane, n-heptane, 1-propanol, and 2-propanol into a dioctyl ether retentive (stationary) phase used in gas-liquid chromatography. The united-atom version of the TraPPE (transferable potentials for phase equilibria) force field was used to model all analytes and the solvent. The analyte partition coefficients, Gibbs free energies of transfer, and Kovats retention indexes were calculated at four different temperatures ranging from 303.15 to 348.15 K. Although hydrogen bonding is a major contributor to the retention of the alcohol analytes over the entire temperature range, its importance for the separation factor between the primary and secondary alcohol decreases substantially with increasing temperature. The enthalpies and entropies for hydrogen bond formation were also estimated from the temperature dependence of the corresponding equilibrium constants. In agreement with experimental measurements, it is observed that the hydrogen bond involving 1-propanol is enthalpically favored, but entropically disfavored compared to 2-propanol.

Sensitivity of the methylbenzenes and chlorobenzenes retention index to column temperature, stationary phase polarity, and number and chemical nature of substituents

Retention indices of methylbenzenes and chlorobenzenes on two fused silica capillary columns, HP-5 (diphenylsiloxane 5% diphenyldimethylsiloxane) and ZB-WAX (polyethylene glycol), have been calculated at various isothermal temperatures and compared with literature data. The retention index temperature effect was studied for each solute, finding greater retention index the higher the column temperature. A comparison between the straight line fit and the fit to the recently proposed equation I = A + B/T +C ln T was carried out. The effect of the stationary phase polarity on the retention index was checked. In general, a greater retention index was found for the more polar stationary phase. The retention indices of the chlorobenzenes are greater than the retention indices of the methylbenzenes, irrespective of the stationary phase and the column temperature. In addition, the influence of the methyl/chlorine substitution on the benzene molecule was investigated at each temperature. The retention indices increased as the number of substituents (methyl/chlorine) increased. The retention index increments of methyl and chloro derivatives are also discussed, which permits to compare the effect of both, methyl or chlorine, chemical functions, for a fixed substituent number in the benzene molecule.

Temperature dependence of the Kovats retention indices for alkyl 1, 3-diketones on a DB-5 capillary column

A series of alkyl 1,3-diketones were used to study the temperature dependence of the Kovats retention indices in the temperature range 130?190 ?C (403?463 K). The temperature dependence is described by the empirical equation I = B B/T + ClnT. On the basis of this equation, the activation enthalpy, ?H?, and the chemical potential of the partitioning of one methylene group between the two phases of the chromatographic system, ??p(CH2), were calculated. Also the Kovats retention indices ? boiling point correlations (linear and reciprocal) for alkyl 1,3-diketones were studied and ??p(CH2) was calculated.

Isothermal Kováts retention indices of sulfur compounds on a poly(5% diphenyl-95% dimethylsiloxane) stationary phase

Isothermal Kováts retention indices of 21 sulfur compounds relevant to the fuel gas and food industries are reported on a poly(5% diphenyl-95% dimethylsiloxane) capillary column stationary phase. Measurements were performed at four temperatures and the temperature dependence of the values modeled with Antoine-type equations. Indices were calculated using a non-linear technique, and the predicted values were found to agree with values obtained using traditional logarithmic predictions. We demonstrate that there is sufficient separation between retention indices to predict the identity of a compound by its retention index.

Retention index calculation without n-alkanes--the virtual carbon number

For the fast gas chromatographic identification of separated components the retention index is still one of the most often used tools, although mass spectrometry is available in almost all analytical laboratories. For the calculation of the retention indices it is not necessary to use n-alkanes or any other homologous series. If the analyte contains some compounds, not necessarily belonging to a homologous series, with well-known retention indices those compounds can be used as index references and the index of the other compounds can be calculated as is done using n-alkanes. The only difference is that instead of the carbon number of the n-alkanes, virtual carbon numbers of the reference compounds should be used. The method of calculation, and the effect of this method of calculation on the reproducibility are discussed in this paper, and analyses of some halogenated compounds and marjoram oil are used as experimental examples.

Minimum in the temperature dependence of the Kováts retention indices of nitroalkanes and alkanenitriles on an apolar phase

The Kováts retention indices (I) of 1-nitroalkanes and alkanenitriles were determined on polydimethylsiloxane and Innowax (polyethylene glycol) columns in a wide temperature range. The temperature dependence of the retention indices exhibits a definite minimum for the early members of the homologous series. The position of the minimum shifts to lower temperatures with increasing carbon atom number of the solute. The thermodynamic explanation of an extreme in the I vs. T function is the higher solvation heat capacities of nitroalkanes and alkanenitriles relative to those of the reference n-alkanes, owing to the deviation from the ideal state in the solution. A novel equation was derived which describes the minimum in the I vs. T function, too.

Temperature dependence of Kováts indices in gas chromatography revisited

Temperature dependence of the Kováts retention index (I) was measured for some aliphatic ketones and aldehydes on a poly(dimethyl siloxane) (HP-1) stationary phase. An interesting minimum (non-linearity) was observed for the I versus isothermal column temperature (T) relationships. A novel empirical model is proposed: I=A+B/T+C ln T, where A, B and Care equation constants and B/C = T(min). A detailed statistical analysis clearly shows superiority of the extended model (i.e., of this containing the logarithm of the temperature (ln T) term) over the earlier established Antoine-type reciprocal equation. The minimum temperature (and the energy like quantity=RT(min), where R is the gas constant) changes in a systematic manner. The factors effecting the (RT(min)) term are as follows: (i) this term decreases with the increase of the molecular mass of the respective oxo compounds; (ii) ketones have higher absolute values of (RT(min) than aldehydes; (iii) branching of the carbon chain lowers the mentioned (RT(min)). This enthalpy term is unambiguously bound to the polarity of solutes.

Temperature effects on the retention of n-alkanes and arenes in helium-squalane gas-liquid chromatography. Experiment and molecular simulation

Experiments and molecular simulations were carried out to study temperature effects (in the range of 323 to 383 K) on the absolute and relative retention of n-hexane, n-heptane, n-octane, benzene, toluene and the three xylene isomers in gas-liquid chromatography. Helium and squalane were used as the carrier gas and retentive phase, respectively. Both the experiments and the simulations show a markedly different temperature dependence of the retention for the n-alkanes compared to the arenes. For example, over the 60 K temperature range studied, the Kovats retention index of benzene is found to increase by about 16 or 18+/-10 retention index units determined from the experiments or simulations, respectively. For toluene and the xylenes, the experimentally measured increases are similar in magnitude and range from 14 to 17 retention index units for m-xylene to o-xylene. The molecular simulation data provide an independent method of obtaining the transfer enthalpies and entropies. The change in retention indices is shown to be the result of the larger entropic penalty and the larger heat capacity for the transfer of the alkane molecules.

Improving the accuracy of Kováts' retention indices in isothermal gas chromatography

Isothermal Kováts' retention indices are currently reported as whole numbers, and are frequently deduced from a linear least mean squares fitting of the logarithms of adjusted retention times of a number of n-alkanes versus carbon number, following an iterative method that minimises errors. The currently accepted accuracy is about one retention index unit for apolar stationary phases, and lower for polar stationary phases. This paper presents results that show how the accuracy of the retention index may be safely reported to one-tenth of a retention index unit by the use of a non-linear equation, with present day gas chromatographs without electronic flow controllers. Results are presented that prove the correctness of the retention indices found for several substances on one particular capillary column. Hints on the minimum retention times needed to achieve the 0.1 retention index accuracy are mentioned, for retention times recorded in minutes and in seconds. According to results of this paper, two chromatograms, run under the appropriate conditions, are sufficient to obtain the desired accuracy. The method proposed in this paper does not require knowledge of the hold-up time of the chromatogram.

Analysis of the equations for the temperature dependence of the retention index. I. Relation between equations

The relation between the parameters of several equations for the retention index temperature dependence was established, taking the hyperbola deduced from the retention theory as starting point. The transformation factors depend only on methylene contributions to the thermodynamic functions of solution and temperature. Their evolution with the mean temperature of the range was illustrated for SE-30 and Carbowax-20M glass capillary columns. On this basis the post-run standardisation of dI/dT values at a reference mean temperature is possible. Examples and statistical correlation between series of parameters from different equations for perfumery solutes were shown.