Simulating Retention in Gas–Liquid Chromatography: Benzene, Toluene, and Xylene Solutes

Accurate predictions of retention times, retention indices, and partition con- stants are a long sought-after goal for theoretical studies in chromatography. Although advances in computational chemistry have improved our understanding of molecular iInteractions, little attention has been focused on chromatography, let alone calculations of retention properties. Configurational-bias MoInte Carlo simulations in the isobaric isothermal Gibbs ensemble were used to investigate the partitioning of benzene, toluene, and the three xylene isomers between a squalane liquid phase and a helium vapor phase. The united-atom representa- tion of the TraPPE (transferable potentials for phase equilibria) force field was used for all solutes and squalane. The Gibbs free energies of transfer and Kovats retention indices of the solutes were calculated directly from the partition con- stants (which were averaged over several independent simulations). While the calculated Kovats indices of benzene and toluene at T=403 K are significantly higher than their experimental couInterparts, much better agreement is found for the xylene isomers at T=365 K.