Retention mechanism for polycyclic aromatic hydrocarbons in reversed-phase liquid chromatography with monomeric stationary phases

The Solvation Shell of Small Solutes in Aqueous-Organic Solvent Mixtures and Its Implications for Reversed-Phase Liquid Chromatography

Reversed-phase liquid chromatography (RPLC) operates with water-organic solvent (W-OS) mobile phases where preferential solvation (PS) of solutes is likely. To investigate the relevance of the solute solvation shell in the mobile phase for RPLC retention, we combine data from molecular dynamics simulations of small, neutral solutes (six analytes and two dead time markers) in W-methanol (MeOH) and W-acetonitrile (ACN) mixtures with corresponding retention data obtained on an RPLC column over a wide range of W/OS ratios. Data derived from Kirkwood-Buff integrals show PS by the OS for analytes vs low or negative PS for dead time markers. W-ACN mixtures generate a higher amount of PS than W-MeOH mixtures, which contributes to the higher eluent strength of ACN in RPLC. Difference spatial distribution functions reveal anisotropic solvation shells with OS excess at hydrocarbon elements and W excess at functional groups, predicting that retention by the hydrophobic stationary phase is favored by hydrocarbon elements and limited by functional groups. Analysis of solute-solvent hydrogen bonds pinpoints the hydrogen-bond requirements toward W as the retention-limiting factor. The relation between the solute solvation shell and retention confirms the importance of W-OS and solute-W hydrogen bonding for RPLC retention.

Separations by Shape: Molecular Shape Recognition in Liquid Chromatography

Molecular shape can provide a basis for chromatographic separations that is distinct from other interaction mechanisms often considered in liquid chromatography. Molecular shape recognition, or shape selectivity, is most evident for the separation of isomeric compounds that have constrained molecular structures, such as polycyclic aromatic hydrocarbons. A central feature of shape-selective columns is conformational order within the stationary phase; this aspect of stationary phase morphology is revealed through spectroscopic studies and molecular dynamics simulations, and is correlated with chromatographic performance.

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.

Molecularly Informed Field Theories from Bottom-up Coarse-Graining

Polymer formulations possessing mesostructures or phase coexistence are challenging to simulate using atomistic particle-explicit approaches due to the disparate time and length scales, while the predictive capability of field-based simulations is hampered by the need to specify interactions at a coarser scale (e.g., χ-parameters). To overcome the weaknesses of both, we introduce a bottom-up coarse-graining methodology that leverages all-atom molecular dynamics to molecularly inform coarser field-theoretic models. Specifically, we use relative-entropy coarse-graining to parametrize particle models that are directly and analytically transformable into statistical field theories. We demonstrate the predictive capability of this approach by reproducing experimental aqueous poly(ethylene oxide) (PEO) cloud-point curves with no parameters fit to experimental data. This synergistic approach to multiscale polymer simulations opens the door to de novo exploration of phase behavior across a wide variety of polymer solutions and melt formulations.

Collective Solvation and Transport at Tetrahydrofuran-Silica Interfaces for Separation of Aromatic Compounds: Insight from Molecular Dynamics Simulations

We have performed umbrella sampling molecular dynamics simulations to study the separation mechanism of aromatic compounds at the tetrahydrofuran (THF)-methanol-silica interface by liquid chromatography. Solute molecules with different polarities (naphthol and naphthalene) are selected as representative aromatic compounds. For the polar solute (naphthol), the free energy profile shows a deep minimum near the THF-silica interface, suggesting strong interactions with the polar surface. When methanol is added to the interface, there is a sharp increase in naphthol's free energy minimum, and the corresponding diffusion dynamics also undergoes a dramatic change. These findings explain the fast separation mechanism in recent experiments of separating fused ring compounds in asphaltenes with liquid chromatography. Further solvation structure and orientation analysis suggest that apolar and polar solutes may find their own comfort zones several angstroms away from the interface, and their phenyl ring's orientations would undergo a parallel-to-perpendicular transition as the solute molecule moves away from the surface. Extending our simulation studies to systems with different solute concentrations reveals that there is a decrease in the adsorption free energy accompanied by enhanced surface diffusion as the solute concentration increases, which is related to the crowding in the interfacial layers. Our simulation analysis gives a detailed microscopic description of solute solvation and transport at the THF-silica chromatography interface and will be helpful for improving separation protocols in future applications.

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.

The role of solute polarity on methanol-silica interfacial solvation: a molecular dynamics study

The solvation structure and dynamics of small organic molecules at the methanol-silica interface are important for understanding the reaction dynamics in heterogeneous catalysis as well as the transport mechanisms in liquid chromatography. The role of solute polarity in interfacial solvation at the methanol-silica interface has been investigated via umbrella sampling molecular dynamics (MD) simulations and 1,3-propanediol and n-pentane were selected as representative species of polar and apolar solutes. Free energy calculations reveal that it took a similar free energy cost to transfer both solute molecules from the interface to the bulk, despite the huge difference in their polarities. The 1,3-propendiol molecule can penetrate the adsorbed methanol layer and form hydrogen bonds with the silica surface with its backbone perpendicular to the silica surface, resulting in a significant slowdown of translational and rotational dynamics. Further analysis of solvent density and solute orientations suggest that at the minimum of the adsorption free energy curve, the 1,3-propanediol molecule is in a desolvated state, while n-pentane is in a solvated state. The collective effect of solute concentration has also been studied by unbiased MD simulations, and the free energy barriers of transferring the solute molecule from the interface to bulk, as well as the parallel diffusion coefficients at the interface, show a non-monotonic dependence on solute concentration, which can be related to the crowded environment in the interfacial layers.

Characterization of stationary phases in supercritical fluid chromatography including exploration of shape selectivity

Achiral packed column supercritical fluid chromatography (SFC) has shown an important regain of interest in academic and industrial laboratories in the recent years. In relation to this increased concern, major instrument manufacturers have designed some stationary phases specifically for SFC use. SFC stationary phases have been widely examined over the last two decades, based on the use of linear solvation energy relationships (LSER), which relate analyte retention to its properties and to the interaction capabilities of the chromatographic system. The method provides some understanding on retention mechanisms (normal phase, reversed phase or mixed-mode) and the possibility to compare stationary phases on a rational basis, especially through a spider diagram providing a visual classification. The latter can be used as a primary tool to select complementary stationary phases to be screened for any separation at early stages of method development, before optimization steps. In this context, the characterization of the 14 columns from the Shim-pack UC series (Shimadzu Corporation, Kyoto, Japan), which are dedicated to SFC and more broadly to unified chromatography (UC), was performed, using the LSER methodology. As in previous works, seven descriptors, including five Abraham descriptors (E, S, A, B, V) and two descriptors describing positive and negative charges (D^- and D⁺) were first employed to describe interactions with neutral and charged analytes. Secondly, two more descriptors were introduced, which were previously employed solely for the characterization of enantioselective systems and expressing shape features of the analytes (flexibility F and globularity G). They brought additional insight into the retention mechanisms, showing how spatial insertion of the analytes in some stationary phases is contributing to shape separation capabilities and how folding possibilities in flexible molecules is unfavorable to retention in other stationary phases.

Molecular dynamics simulations of amino acid adsorption and transport at the acetonitrile–water–silica interface: the role of side chains

The solvation and transport of amino acid residues at liquid–solid interfaces have great importance for understanding the mechanism of separation of biomolecules in liquid chromatography. This study uses umbrella sampling molecular dynamics simulations to study the adsorption and transport of three amino acid molecules with different side chains (phenylalanine (Phe), leucine (Leu) and glutamine (Gln)) at the silica–water–acetonitrile interface in liquid chromatography. Free energy analysis shows that the Gln molecule has stronger binding affinity than the other two molecules, indicating the side chain polarity may play a primary role in adsorption at the liquid–solid interface. The Phe molecule with a phenyl side chain exhibits stronger adsorption free energy than Leu with a non-polar side chain, which can be ascribed to the better solvated configuration of Phe. Further analysis of molecular orientations found that the amino acid molecules with apolar side chains (Phe and Leu) have ‘standing up’ configurations at their stable adsorption state, where the polar functional groups are close to the interface and the side chain is far from the interface, whereas the amino acid molecule with a polar side chain (Gln) chooses the ‘lying’ configuration, and undergoes a sharp orientation transition when the molecule moves away from the silica surface. Extending our simulation studies to systems with different solute concentrations reveals that there is a decrease in the adsorption free energy as well as surface diffusion as the solute concentration increases, which is related to the crowding in the interfacial layers. This simulation study gives a detailed microscopic description of amino acid molecule solvation and transport at the acetonitrile–water–silica interface in liquid chromatography and will be helpful for understanding the retention mechanism for amino acid separation.

The solvation and transport of amino acid residues at liquid–solid interfaces have great importance for understanding the mechanism of separation of biomolecules in liquid chromatography.

The combination of partition, size exclusion, and hydrodynamic models in chromatography, and application to bonded phases on porous supports

Liquid-liquid partition chromatography has been used for many years as a model and teaching introduction to column chromatography. However, the partition model does not describe separations on bonded phases with porous supports particularly well, especially regarding the thermodynamics controlling solute distribution. Further difficulties arise when more than one mechanism is involved in solute retention. Nomenclature is not perfectly aligned with the underlying thermodynamic descriptors and is inconsistently applied over various chromatographic techniques. Presented here is a general description of retention that spans partition, size exclusion, and hydrodynamic separation processes, and is then applied to bonded-phase separations on porous supports. The model provides a general description applicable to adsorption, reversed-phase, hydrophilic interaction, size-exclusion, hydrodynamic chromatography, and any combination of these techniques including liquid chromatography at the critical condition. Further expansion to include retention by ion-exchange and field-flow fractionation appears to be possible. Recommendations on retention factor definition and evaluation are given.

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.

Ultrathin Molecular-Layer-by-Layer Polyamide Membranes: Insights from Atomistic Molecular Simulations

In this study, we present an atomistic simulation study of several physicochemical properties of polyamide (PA) membranes formed from interfacial polymerization or from a molecular-layer-by-layer (mLbL) on a silicon wafer. These membranes are composed of meta-phenylenediamine (MPD) and benzene-1,3,5-tricarboxylic acid chloride (TMC) for potential reverse osmosis (RO) applications. The mLbL membrane generation procedure and the force field models were validated, by comparison with available experimental data, for hydrated density, membrane swelling, and pore size distributions of PA membranes formed by interfacial polymerization. Physicochemical properties such as density, free volume, thickness, the degree of cross-linking, atomic compositions, and average molecular orientation (which is relevant for the mLbL membranes) are compared for these different processes. The mLbL membranes are investigated systematically with respect to TMC monomer growth rate per substrate surface area, MPD/TMC ratio, and the number of mLbL deposition cycles. Atomistic simulations show that the mLbL deposition generates membranes with a constant film growth if both the TMC monomer growth rate and MPD/TMC monomer ratio are kept constant. The film growth rate increases with TMC monomer growth rate or MPD/TMC ratio. Furthermore, it was found on one hand that the mLbL membrane density and free volume varies significantly with respect to the TMC monomer growth rate, while on the other hand the degree of cross-linking and the atomic composition varies considerably with the MPD/TMC ratio. Additionally, it was found that both TMC and MPD orient at a tilted angle with respect to the substrate surface, where their angular distribution and average angle orientation depend on both the TMC growth rate and the number of deposition cycles. This study illustrates that molecular simulations can play a crucial role in the understanding of structural properties that can empower the design of the next generation RO membranes created from molecular-layer-by-layer (mLbL) on a silicon wafer.

Influence of phase type and solute structure on changes in retention with pressure in reversed-phase high performance liquid chromatography

The influence of pressure on the retention of several types of solute, including acids, bases and neutrals, was studied by the use of restriction capillaries added to the end of various monomeric and polymeric octadecylsilyl-modified 5μm particle size columns. Although it appeared that certain polymeric columns could give somewhat greater increases in retention with pressure, differences in behaviour between these different C18 columns were rather small. Differences in solute molecular size were most important in determining increases in retention with pressure. However, solute structure such as polarity and planarity were also influential. A prototype C30 column gave interesting selectivity changes between planar and non-planar solutes as a function of pressure. Considerable selectivity differences with pressure were shown when diverse mixtures of solutes were analysed. For the solutes studied, only minor effects of increased pressure on column efficiency and peak shape were noted.

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.

Multicomponent adsorption of alcohols onto silicalite-1 from aqueous solution: isotherms, structural analysis, and assessment of ideal adsorbed solution theory

Configurational-bias Monte Carlo (CBMC) simulations in the isobaric-isothermal version of the Gibbs ensemble (GE) were carried out to probe the adsorption from aqueous solutions of methanol and/or ethanol onto silicalite-1. This methodology does require neither specification of the chemical potential nor any reference to activity models based on experimental data. The CBMC-GE methodology can be applied to the complete range of mixture compositions from pure water to pure alcohols and can also be used when multiple solute types are present at high concentration. The simulations demonstrate high selectivities for the alcohols (α(ethanol) > α(methanol)) almost over the entire composition range. The ideal adsorbed solution theory is found to substantially underpredict the amount of sorbed water and leads to very large errors for low alcohol solution concentrations. The simulations indicate that, at lower loadings, the adsorbed alcohol molecules can serve as seeds for water adsorption but, at higher loadings, alcohols displace water molecules from their preferred region.