The modeling of overloaded elution band profiles in supercritical fluid chromatography

A new metric for relating macroscopic chromatograms to microscopic surface dynamics: the distribution function ratio (DFR)

Heterogeneous stationary phase chemistry causes chromatographic tailing that lowers separation efficiency and complicates optimizing mobile phase conditions. Model-free metrics are attractive for assessing optimal separation conditions due to the low quantity of information required, but often do not reveal underlying mechanisms that cause tailing, for example, heterogeneous retention modes. We report a new metric, which we call the Distribution Function Ratio (DFR), based on a graphical comparison between the chromatogram and Gaussian cumulative distribution functions, achieving correspondence to ground truth surface dynamics with a single chromatogram. Using a Monte Carlo framework, we show that the DFR can predict the prevalence of heterogeneous retention modes with high precision when the relative desorption rate between modes is known, as in during surface dynamics experiments. Ground truth comparisons reveal that the DFR outperforms both the asymmetry factor and skewness by yielding a one-to-one correspondence with heterogeneous retention mode prevalence over a broad range of experimentally realistic values. Perhaps of more value, we illustrate that the DFR, when combined with the asymmetry factor and skewness, can estimate microscopic surface dynamics, providing valuable insights into surface chemistry using existing chromatographic instrumentation. Connecting ensemble results to microscopic quantities through the lens of simulation establishes a new chemistry-driven route to measuring and advancing separations.

Predictions of overloaded concentration profiles in supercritical fluid chromatography

Here, overloaded concentration profiles were predicted in supercritical fluid chromatography using a combined two-dimensional heat and mass transfer model. The heat balance equation provided the temperature and pressure profiles inside the column. From this the density, viscosity, and mobile phase velocity profiles in the column were calculated. The adsorption model is here expressed as a function of the density and temperature of the mobile phase. The model system consisted of a Kromasil Diol column packed with 2.2-µm particles (i.e., a UHPSFC column) and the solute was phenol eluted with neat carbon dioxide at three different outlet pressures and five different mobile phase flow rates. The proposed model successfully predicted the eluted concentration profiles in all experimental runs with good agreement even with high-density drops along the column. It could be concluded that the radial temperature and density gradients did not significantly influence the overloaded concentration elution profiles.