Influence of bonded-phase coverage in reversed-phase liquid chromatography via molecular simulation

Organic-solvent ditch overlap in reversed-phase liquid chromatography: A molecular dynamics simulation study in cylindrical 6-12 nm-diameter pores

Mass transport through the mesopore space of a reversed-phase liquid chromatography (RPLC) column depends on the properties of the chromatographic interface, particularly on the extent of the organic-solvent ditch that favors the analyte surface diffusivity. Through molecular dynamics simulations in cylindrical RPLC mesopore models with pore diameters between 6 and 12 nm we systematically trace the evolution of organic-solvent ditch overlap due to spatial confinement in the mesopore space of RPLC columns for small-molecule separations. Each pore model of a silica-based, endcapped, C18-stationary phase is equilibrated with two mobile phases of comparable elution strength, namely 70/30 (v/v) water/acetonitrile and 60/40 (v/v) water/methanol, to consider the influence of the mobile-phase composition on the onset of organic-solvent ditch overlap. The simulations show that, as the pore diameter decreases from 9 to 6 nm, the bonded-phase density extends and compacts towards the pore center, which leads to increased accumulation of organic-solvent excess and thus enhanced organic-solvent diffusivity in the ditch. Because the acetonitrile ditch is more pronounced than the methanol ditch, acetonitrile ditch overlap sets in at less severe spatial confinement than methanol ditch overlap. The pore-averaged methanol and acetonitrile diffusivities are considerably raised by ditch overlap in the 6 nm-diameter pore, but also benefit from the ditch (without overlap) in the 7 to 12 nm-diameter pores, whereby local and pore-averaged effects are generally larger for acetonitrile than methanol.

Mobile-Phase Contributions to Organic-Solvent Excess Adsorption and Surface Diffusion in Reversed-Phase Liquid Chromatography

Fast transport of retained analytes in reversed-phase liquid chromatography occurs through surface diffusion in the organic-solvent (OS)-enriched interfacial "ditch" region between the hydrophobic stationary phase and the water (W)-OS mobile phase. Through molecular dynamics simulations that recover the OS excess adsorption isotherms of a typical C₁₈-stationary phase for methanol and acetonitrile, we explore the relation between OS properties, OS excess adsorption, and surface diffusion. The emerging molecular-level picture attributes the mobile-phase contribution to surface diffusion to the hydrogen-bond capability and the eluting power of the OS. The higher affinity of methanol for the formation of W-OS hydrogen bonds at the soft, hydrophobic surface presented by the bonded-phase (C₁₈) chains reduces the OS excess and the related viscosity drop in the ditch. The lower eluting power of methanol, however, translates to increased bonded-phase contacts for analytes, which can increase their mobility gain from surface diffusion above the gain observed with acetonitrile.

Ethanol and Water Adsorption in Conventional and Hierarchical All-Silica MFI Zeolites

Hierarchical zeolites containing both micro- (<2 nm) and mesopores (2-50 nm) have gained increasing attention in recent years because they combine the intrinsic properties of conventional zeolites with enhanced mass transport rates due to the presence of mesopores. The structure of the hierarchical self-pillared pentasil (SPP) zeolite is of interest because all-silica SPP consists of orthogonally intergrown single-unit-cell MFI nanosheets and contains hydrophilic surface silanol groups on the mesopore surface while its micropores are nominally hydrophobic. Therefore, the distribution of adsorbed polar molecules, like water and ethanol, in the meso- and micropores is of fundamental interest. Here, molecular simulation and experiment are used to investigate the adsorption of water and ethanol on SPP. Vapor-phase single-component adsorption shows that water occupies preferentially the mesopore corner and surface regions of the SPP material at lower pressures (P/P ₀ < 0.5) while loading in the mesopore interior dominates adsorption at higher pressures. In contrast, ethanol does not exhibit a marked preference for micro- or mesopores at low pressures. Liquid-phase adsorption from binary water-ethanol mixtures demonstrates a 2 orders of magnitude lower ethanol/water selectivity for the SPP material compared to bulk MFI. For very dilute aqueous solutions of ethanol, the ethanol molecules are mostly adsorbed inside the SPP micropore region due to stronger dispersion interactions and the competition from water for the surface silanols. At high ethanol concentrations (C _EtOH > 700 g L^-1), the SPP material becomes selective for water over ethanol.

Consequences of Cylindrical Pore Geometry for Interfacial Phenomena in Reversed-Phase Liquid Chromatography

The interfacial phenomena behind analyte separation in a reversed-phase liquid chromatography column take place nearly exclusively inside the silica mesopores. Their cylindrical geometry can be expected to shape the properties of the chromatographic interface with consequences for the analyte density distribution and diffusivity. To investigate this topic through molecular dynamics simulations, we introduce a cylindrical pore inside a slit pore configuration, where the inner curved and outer planar silica surface bear the same bonded phase. The present model replicates an average-sized (9 nm) mesopore in an endcapped C₁₈ column equilibrated with a mobile phase of 70/30 (v/v) water/acetonitrile. Simulations performed for ethylbenzene and acetophenone show that the surface curvature shifts the bonded phase and analyte density toward the pore center, decreases the solvent density in the bonded-phase region, increases the acetonitrile excess in the interfacial region, and considerably enhances the surface diffusivity of both analytes. Overall, the cylindrical pore provides a more hydrophobic environment than the slit pore. Ethylbenzene density is decidedly increased in the cylindrical pore, whereas acetophenone density is nearly equally distributed between the cylindrical and slit pore. The cylindrical pore geometry thus sharpens the discrimination between the apolar and moderately polar analytes while enhancing the mass transport of both.

Insights from molecular simulations about dead time markers in reversed-phase liquid chromatography

Among the most popular compounds to estimate the hold-up time in reversed-phase liquid chromatography (RPLC) are acetone and uracil, which are considered as too small and too polar, respectively, for retention by the hydrophobic stationary phase, although their observed elution behavior does not fully support this assumption. We investigate how acetone and uracil as solutes interact with the chromatographic interface through molecular dynamics simulations in an RPLC mesopore model of a silica-supported, endcapped, C₁₈ phase equilibrated with a water (W)‒acetonitrile (ACN) mobile phase. The simulation results provide a molecular-level explanation for the observed elution behavior of acetone and uracil, but also question whether true dead time markers for RPLC exist. Both solutes have a density maximum in the interfacial region in addition to a low presence in the bonded-phase region, but these density peaks clearly differ from the adsorption and partitioning peaks of true analytes. Acetone partially behaves like a co-solvent of ACN and partially like the analyte acetophenone. Like ACN, acetone can be found in the first and second layer of solvent molecules at the silica surface; like acetophenone, acetone adsorbs to the bonded-phase chains by orienting its polar group to the bulk region to sustain hydrogen bonds with W molecules. Uracil behavior is governed by a need for extensive hydrogen-bond coordination by W molecules. Uracil adsorbs to the very edge of the bonded-phase chains, on the bulk-region side of the ACN density maximum in the interfacial region. Further penetration into the chains is prevented by the absence of W molecules, which are not found deeper in the bonded phase, except at the silica surface. Contrary to true analytes, accumulation of uracil and acetone in the interfacial region ceases at an equimolar presence of W and ACN in the mobile phase (at 70‒80% ACN volume fraction). Uracil achieves a closer approximation of the stationary-phase limit than acetone, but carries the risk of HILIC retention at high ACN fraction in the mobile phase.

Column Characterization and Selection Systems in Reversed-Phase High-Performance Liquid Chromatography

Reversed-phase high-performance liquid chromatography (RP-HPLC) is the most popular chromatographic mode, accounting for more than 90% of all separations. HPLC itself owes its immense popularity to it being relatively simple and inexpensive, with the equipment being reliable and easy to operate. Due to extensive automation, it can be run virtually unattended with multiple samples at various separation conditions, even by relatively low-skilled personnel. Currently, there are >600 RP-HPLC columns available to end users for purchase, some of which exhibit very large differences in selectivity and production quality. Often, two similar RP-HPLC columns are not equally suitable for the requisite separation, and to date, there is no universal RP-HPLC column covering a variety of analytes. This forces analytical laboratories to keep a multitude of diverse columns. Therefore, column selection is a crucial segment of RP-HPLC method development, especially since sample complexity is constantly increasing. Rationally choosing an appropriate column is complicated. In addition to the differences in the primary intermolecular interactions with analytes of the dispersive (London) type, individual columns can also exhibit a unique character owing to specific polar, hydrogen bond, and electron pair donor-acceptor interactions. They can also vary depending on the type of packing, amount and type of residual silanols, "end-capping", bonding density of ligands, and pore size, among others. Consequently, the chromatographic performance of RP-HPLC systems is often considerably altered depending on the selected column. Although a wide spectrum of knowledge is available on this important subject, there is still a lack of a comprehensive review for an objective comparison and/or selection of chromatographic columns. We aim for this review to be a comprehensive, authoritative, critical, and easily readable monograph of the most relevant publications regarding column selection and characterization in RP-HPLC covering the past four decades. Future perspectives, which involve the integration of state-of-the-art molecular simulations (molecular dynamics or Monte Carlo) with minimal experiments, aimed at nearly "experiment-free" column selection methodology, are proposed.

From in silica to in silico: retention thermodynamics at solid–liquid interfaces

The dynamics of solvated molecules at the solid/liquid interface is essential for a molecular-level understanding for the solution thermodynamics in reversed phase liquid chromatography (RPLC).

Molecular Mechanisms Underlying Solute Retention at Heterogeneous Interfaces

Despite considerable effort, a molecular-level understanding of the mechanisms governing adsorption/desorption in reversed-phase liquid chromatography is still lacking. This impedes rational design of columns and the development of reliable, computationally more efficient approaches to predict the selectivity of a particular column design. Using state-of-the art, validated force fields and free-energy simulations, the adsorption thermodynamics of benzene derivatives is investigated in atomistic detail and provides a quantitative microscopic understanding of retention when compared with experimental data. It is found that pure partitioning or pure adsorption is rather the exception than the rule. Typically, a pronounced ∼1 kcal/mol stabilization on the surface is accompanied by a broad trough indicative of partitioning before the probe molecule incorporates into the mobile phase. The present findings provide a quantitative and rational basis to develop improved effective, coarse-grained computational models and to design columns for specific applications.

Comparative Study of the Effect of Defects on Selective Adsorption of Butanol from Butanol/Water Binary Vapor Mixtures in Silicalite-1 Films

A promising route for sustainable 1-butanol (butanol) production is ABE (acetone, butanol, ethanol) fermentation. However, recovery of the products is challenging because of the low concentrations obtained in the aqueous solution, thus hampering large-scale production of biobutanol. Membrane and adsorbent-based technologies using hydrophobic zeolites are interesting alternatives to traditional separation techniques (e.g., distillation) for energy-efficient separation of butanol from aqueous mixtures. To maximize the butanol over water selectivity of the material, it is important to reduce the number of hydrophilic adsorption sites. This can, for instance, be achieved by reducing the density of lattice defect sites where polar silanol groups are found. The density of silanol defects can be reduced by preparing the zeolite at neutral pH instead of using traditional synthesis solutions with high pH. In this work, binary adsorption of butanol and water in two silicalite-1 films was studied using in situ attenuated total reflectance-Fourier transform infrared (ATR-FTIR) spectroscopy under equal experimental conditions. One of the films was prepared in fluoride medium, whereas the other one was prepared at high pH using traditional synthesis conditions. The amounts of water and butanol adsorbed from binary vapor mixtures of varying composition were determined at 35 and 50 °C, and the corresponding adsorption selectivities were also obtained. Both samples showed very high selectivities (100-23 000) toward butanol under the conditions studied. The sample having low density of defects, in general, showed ca. a factor 10 times higher butanol selectivity than the sample having a higher density of defects at the same experimental conditions. This difference was due to a much lower adsorption of water in the sample with low density of internal defects. Analysis of molecular simulation trajectories provides insights on the local selectivities in the zeolite channel network and at the film surface.

Hydrophobicity of carbohydrates and related hydroxy compounds

The hydrophobic interaction of carbohydrates and other hydroxy compounds with a C18-modified silica gel column was measured with pure water as eluent, thereby expanding the range of measurements already published. The interaction is augmented by structure strengthening salts and decreasing temperature. Although the interaction of the solute with the hydrophobic interface is expected to only imperfectly reflect its state in aqueous bulk solution, the retention can be correlated to hydration numbers calculated from molecular mechanics studies given in the literature. No correlation can be established towards published hydration numbers obtained by physical methods (isentropic compressibility, O-17 NMR relaxation, terahertz spectroscopy, and viscosity). The hydrophobicity is discussed with respect to the chemical structure. It increases with the fraction and size of hydrophobic molecular surface regions.

Tuning the Flight Length of Molecules Diffusing on a Hydrophobic Surface

Transport at solid-liquid interfaces is critical to self-assembly, biosensing, and heterogeneous catalysis, but surface diffusion remains difficult to characterize and rationally manipulate, due to the inherent heterogeneity of adsorption on solid surfaces. Using single-molecule tracking, we characterized the diffusion of a fluorescent long-chain surfactant on a hydrophobic surface, which involved periods of confinement alternating with bulk-mediated "flights". The concentration of methanol in solution was varied to tune the strength of the hydrophobic surface-molecule interaction. The frequency of confinement had a nonmonotonic dependence on methanol concentration that reflected the relative influence of anomalously strong adsorption sites. By carefully accounting for the effect of this surface heterogeneity, we demonstrated that flight lengths increased monotonically as the hydrophobic attraction decreased, in agreement with theoretical predictions for bulk-mediated surface diffusion. The theory provided an accurate description of surface diffusion, despite the system being heterogeneous, and can be leveraged to optimize molecular search and assembly processes.

Adsorption-induced changes of the structure of the tethered chain layers in a simple fluid

We use density functional theory to study the influence of fluid adsorption on the structure of grafted chain layer. The chains are modeled as freely jointed spheres. The chain segments and spherical molecules of the fluid interact via the Lennard-Jones potential. The fluid molecules are attracted by the substrate. We calculate the excess adsorption isotherms, the average height of tethered chains, and the force acting on selected segments of the chains. The parameters that were varied include the length of grafted chains, the grafting density, the parameters characterizing fluid-chain and fluid-surface interactions, the bulk fluid density, and temperature. We show that depending on the density of the bulk fluid the height of the bonded layer increases, remains constant, or decreases with increasing temperature.

Molecular simulation studies of reversed-phase liquid chromatography

Over the past 20 years, molecular simulation methods have been applied to the modeling of reversed-phase liquid chromatography (RPLC). The purpose of these simulations was to provide a molecular-level understanding of: (i) the structure and dynamics of the bonded phase and its interface with the mobile phase, (ii) the interactions of analytes with the bonded phase, and (iii) the retention mechanism for different analytes. However, the investigation of chromatographic systems poses significant challenges for simulations with respect to the accuracy of the molecular mechanics force fields and the efficiency of the sampling algorithms. This review discusses a number of aspects concerning molecular simulation studies of RPLC systems including the historical development of the subject, the background needed to understand the two prevalent techniques, molecular dynamics (MD) and Monte Carlo (MC) methods, and the wealth of insight provided by these simulations. Examples from the literature employing MD approaches and from the authors' laboratory using MC methods are discussed. The former can provide information on chain dynamics and transport properties, whereas the latter techniques are uniquely suited for the investigation of phase and sorption equilibria that underly RPLC retention, and both can be used to elucidate the bonded-chain conformations and solvent distributions.

Molecular simulations of retention in chromatographic systems: use of biased Monte Carlo techniques to access multiple time and length scales

The use of configurational-bias Monte Carlo simulations in the Gibbs ensemble allows for the sampling of phenomena that occur on vastly different time and length scales. In this review, applications of this simulation approach to probe retention in gas and reversed-phase liquid chromatographic systems are discussed. These simulations provide an unprecedented view of the retention processes at the molecular-level and show excellent agreement with experimental retention data.

Computational study of enantioseparation by amylose tris(3,5-dimethylphenylcarbamate)-based chiral stationary phase

The mechanism of chiral separation on amylose tris(3,5-dimethylphenylcarbamate) is studied with docking simulations of enantiomers by molecular dynamics. All-atom models of amylose tris(3,5-dimethylphenylcarbamate) on the modified silica gel surface were constructed for the docking simulations of metalaxyl and benalaxyl. The elution orders and energetic differences were also predicted based on the intermolecular interactions, which were in agreement with the experimental results. The radial distribution function was employed to analyze the structural features of the enantiomer-chiral stationary phase complex and used to elucidate the mechanism of chiral separation. The separation of metalaxyl and benalaxyl is mainly controlled by the hydrogen bond. And the binding sites had slight differences for the pair of enantiomers, but obvious differences between different chemicals.

Study of solvent adsorption on chemically bonded stationary phases by microcalorimetry and liquid chromatography

A detailed, molecular-level description of the sorption mechanism in reversed-phase liquid chromatography is of great interest to analytical chemists. For this purpose, solvent adsorption in the octadecyl stationary bonded phase was investigated. Preferential adsorption of solvents from an acetonitrile-water and methanol-water mobile phase was measured on a series of non-end-capped octadecyl bonded phases with different coverage densities of bonded ligands using the minor disturbance method. For a comparison, a microcalorimetric study of organic solvent adsorption on the stationary phase was executed. The results from the excess isotherm measurement agree well with the experimental measurement of the heat of immersion of the bonded stationary phases by the test solvents. The microcalorimetric measurement is another method for determination of solvation processes of the stationary phases. Changes of the heat of immersion provide information about the surface accessibility for interaction with solvent molecules. The increase of the stationary phase coverage density reduces the free space between bonded chains and penetration of solvent between organic chains.