Importance of sample intraparticle diffusivity in investigations of the mass transfer mechanism in liquid chromatography

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.

Preparation and investigation of two-component silica-modified hydrophobic layer for minimizing retention loss of reversed-phase chromatographic columns using fully aqueous mobile phases

Chromatographers often employ fully aqueous mobile phases to retain highly polar compounds in reversed-phase liquid chromatography (RPLC). However, when the flow rate is interrupted, either accidentally or intentionally, a substantial loss in retention occurs due to the spontaneous dewetting of water from the hydrophobic surface of conventional RPLC-C18 particles. Previous studies have shown that maintaining a low C18 surface coverage (approximately 1.5 μmol/m2) can mitigate water dewetting by increasing chain disorder, facilitating the intercalation of water clusters between the C18-bonded chains, and keeping the mesopores wetted. In this research, we explore the potential and additional benefits of using two-component surface bonding materials (C8/C18 and PhenylHexyl (PhHx)/C18) at a constant and low total surface coverage of 1.51 ± 0.15 μmol/m2. We synthesized seven one- and two-component modified silica particles with a volume average particle size of 5.22 μm and an average mesopore size of 104 Å. The surface coverage was increased from 0 to 0.54, 1.00, and to 1.66 μmol2 for C8 chains and from 0 to 0.52, 0.70, and to 1.65 μmol2 for PhHx ligands. To prevent interactions between water and any unreacted silanols, all seven derivatized particles were heavily endcapped with trimethylsilane (TMS) reagent. The fraction of the surface area remaining in contact with water was determined by measuring the retention times of weakly (thiourea) and strongly (thymine) retained compounds at intervals of 1, 2, 4, 8, 16, 32, and 64 minutes following the cessation of flow. Two distinct column temperatures, 24°C and 60°C, were employed in the experiments. Retention losses were found to be minimized in the presence of a small quantity of C8 chains (less than 40% of the total surface coverage). Additionally, it is essential to consider substantial fractions of PhHx chains, as long as the presence of the PhHx ligand does not significantly impact retention and selectivity. Combining mixed RPLC bondings with a low total surface coverage of approximately 1.5 μmol/m2 emerges as a viable solution for further minimizing retention loss in standard C18-bonded RPLC columns, particularly within the surface coverage range of 2.5-3.0 μmol/m2.

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.

Retention and mass transfer properties of the series of unbonded, amide-bonded, and alkylsulfobetaine-bonded ethylene bridged hybrid hydrophilic interaction liquid chromatography columns

This work investigates the link between the retentivity and the stationary phase to mobile phase mass transfer resistance of hydrophilic interaction liquid chromatography (HILIC) columns packed with the same base ethylene-bridged hybrid particles (BEH). The retention volumes, the plate heights, and the volume of the adsorbed water layer were measured for the ACQUITY^TM UPLC^TM BEH^TM 130 Å HILIC Column (unbonded BEH), ACQUITY UPLC BEH 130 Å Amide Column (amide group attached), and Atlantis^TM Premier BEH 95 Å Z-HILIC (zwitterionic group attached) Column. The method of Guo (toluene retention volumes in pure acetonitrile and in the HILIC eluent) was validated from the UNIFAC group-contribution method and applied to measure accurately the water layer volumes in these columns. A strong correlation was found between the retention volumes of most neutral polar analytes and the volume of the water layer adsorbed in the HILIC column. The fraction of the pore volume occupied by the water layer increases significantly from the BEH HILIC Column to the BEH Amide Column, and to the BEH Z-HILIC Column. This is explained by the water solvation of the attached ligands in the pore volume of the BEH Particles and to the smaller average mesopore size of the BEH Z-HILIC Particles. A second and strong correlation is also observed between the water content in the HILIC particle and the stationary phase to mobile phase mass transfer resistance of the HILIC columns at high mobile phase linear velocities. The measured intra-particle diffusivity normalized to the bulk diffusion coefficient decreased from 0.33 (BEH HILIC Column) to 0.10 (BEH Amide Column) and to only 0.03 (BEH Z-HILIC Column) for comparable retention of cytosine. These results are fully consistent with the higher viscosity of the internal eluent (higher water content) and higher internal obstruction for diffusion (smaller mesopores and internal porosity) in the BEH Z-HILIC Particles. Still, in gradient elution mode, the peak capacity was found to be 18% higher for the BEH Z-HILIC Column than that on the BEH Amide Column because the retention factors at elution were smaller when maintaining the same analysis time and starting eluent composition.

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.

Confinement Effects on Distribution and Transport of Neutral Solutes in a Small Hydrophobic Nanopore

Molecular dynamics simulations are used to study confinement effects in small cylindrical silica pores with extended hydrophobic surface functionalization as realized, for example, in reversed-phase liquid chromatography (RPLC) columns. In particular, we use a 6 nm cylindrical and a 10 nm slit pore bearing the same C₁₈ stationary phase to compare the conditions inside the smaller-than-average pores within an RPLC column to column-averaged properties. Two small, neutral, apolar to moderately polar solutes are used to assess the consequences of spatial confinement for typical RPLC analytes with water (W)-acetonitrile (ACN) mobile phases at W/ACN ratios between 70/30 and 10/90 (v/v). The simulated data show that true bulk liquid behavior, as observed over an extended center region in the 10 nm slit pore, is not recovered within the 6 nm cylindrical pore. Instead, the ACN-enriched solvent layer around the C₁₈ chain ends (the ACN ditch), a general feature of hydrophobic interfaces equilibrated with aqueous-organic liquids, extends over the entire pore lumen of the small cylindrical pore. This renders the entire pore a highly hydrophobic environment, where, contrary to column-averaged behavior, neither the local nor the pore-averaged sorption and diffusion of analytes scales directly with the W/ACN ratio of the mobile phase. Additionally, the solute polarity-related discrimination between analytes is enhanced. The consequences of local ACN ditch overlap in RPLC columns are reminiscent of ion transport in porous media with charged surfaces, where electrical double-layer overlap occurring locally in smaller pores leads to discrimination between co- and counterionic species.

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.

Enantioselective adsorption dynamics of leucyl-leucine in a Chirobiotic R column

The dynamics of adsorption of the Leu-Leu stereoisomers in a chromatographic column packed with the Chirobiotic R chiral stationary phase bearing grafted antibiotic ristocetin A was studied by means of measurement and analysis of van Deemter plots. Similar measurements were carried out with weakly retained Gly-Gly for the sake of comparison. The bulk diffusion coefficients of the investigated dipeptides were also determined. It is found that the van Deemter plots of both the Leu-Leu stereoisomers and Gly-Gly have an uncommon convex-upward shape. Besides, the van Deemter B coefficients for the Leu-Leu stereoisomers, but not for Gly-Gly, have unusually high values. It is suggested that a high transcolumn contribution to eddy dispersion, which turned out to be enantioselective, accounts for these findings. Adsorption kinetics of all the dipeptides considered is relatively slow, the adsorption rate constant (k_ads) being of order of magnitude 20-60 s^-1. k_ads does not depend on the configuration of Leu-Leu stereoisomers, although their affinity toward the chiral selector depends on this factor. This supports the above hypothesis that eddy dispersion is mainly responsible for the observed peculiarities in the dynamic behavior of dipeptides, and adsorption kinetics has secondary importance in this phenomenon.

Estimating and leveraging protein diffusion on ion-exchange resin surfaces

Protein mobility at solid-liquid interfaces can affect the performance of applications such as bioseparations and biosensors by facilitating reorganization of adsorbed protein, accelerating molecular recognition, and informing the fundamentals of adsorption. In the case of ion-exchange chromatographic beads with small, tortuous pores, where the existence of surface diffusion is often not recognized, slow mass transfer can result in lower resin capacity utilization. We demonstrate that accounting for and exploiting protein surface diffusion can alleviate the mass-transfer limitations on multiple significant length scales. Although the surface diffusivity has previously been shown to correlate with ionic strength (IS) and binding affinity, we show that the dependence is solely on the binding affinity, irrespective of pH, IS, and resin ligand density. Different surface diffusivities give rise to different protein distributions within the resin, as characterized using confocal microscopy and small-angle neutron scattering (length scales of micrometer and nanometer, respectively). The binding dependence of surface diffusion inspired a protein-loading approach in which the binding affinity, and hence the surface diffusivity, is modulated by varying IS. Such gradient loading increased the protein uptake efficiency by up to 43%, corroborating the importance of protein surface diffusion in protein transport in ion-exchange chromatography.

Unusual van Deemter plots of optical isomers on a chiral brush-type liquid chromatography column

Unusual dynamic behavior of the enantiomers of 9-bromo-11b-(tert-butyl)-2,3,6,11b-tetrahydrooxazolo[3',2':1,5]pyrrolo[3,4-b]quinoline-5,11-dione (I) was observed on a Nautilus-R column packed with silica grafted with antibiotic ristocetin. It consisted in (i) antibatic behavior of the van Deemter plots of the enantiomers and (ii) high and strongly enantiomer dependent values of the A- and B-terms of the van Deemter equation. Although rare, such a pattern has been found earlier in chiral chromatography, with all reported cases limited to brush-type chiral stationary phases. Adsorption dynamics in this system was studied by means of the moment method and the peak parking technique; hydrodynamic properties of the column were explored by using unretained tracers. It was shown that the peculiar shape of the van Deemter curves for the enantiomers of I is conditioned by imperfect packing of the stationary phase, which result in high transcolumn eddy dispersion, and by slow adsorption/desorption kinetics. It was proven that the whole void volume of the column available to an eluent is not accessible to the studied analyte because it cannot penetrate the space between neighboring grafted ligands. Its mass transfer in pores is also affected by the fact that the stagnant layer of the binary water-acetonitrile mobile phase differs in composition from the bulk liquid due to preferential adsorption of water that influences the apparent molecular diffusivity of solutes. An effect of the structures of analyte and chiral selector on the adsorption kinetics is also briefly discussed.

Retention and effective diffusion of model metabolites on porous graphitic carbon

The study of metabolites in biological samples is of high interest for a wide range of biological and pharmaceutical applications. Reversed phase liquid chromatography is a common technique used for the separation of metabolites, but it provides little retention for polar metabolites. An alternative to C18 bonded phases, porous graphitic carbon has the ability to provide significant retention for both non-polar and polar analytes. The goal of this work is to study the retention and effective diffusion properties of porous graphitic carbon, to see if it is suitable for the wide injection bands and long run times associated with long, packed capillary-scale separations. The retention of a set of standard metabolites was studied for both stationary phases over a wide range of mobile phase conditions. This data showed that porous graphitic carbon benefits from significantly increased retention (often >100 fold) under initial gradient conditions for these metabolites, suggesting much improved ability to focus a wide injection band at the column inlet. The effective diffusion properties of these columns were studied using peak-parking experiments with the standard metabolites under a wide range of retention conditions. Under the high retention conditions, which can be associated with retention after injection loading for gradient separations, Deff/Dm∼0.1 for both the C18-bonded and porous graphitic carbon columns. As C18 bonded particles are widely, and successfully utilized for long gradient separations without issue of increasing peak width from longitudinal diffusion, this suggests that porous graphitic carbon should be amenable for long runtime gradient separations as well.

The rationale for the optimum efficiency of columns packed with new 1.9μm fully porous Titan-C18 particles-a detailed investigation of the intra-particle diffusivity

In a previous report, it was reported that columns packed with fully porous 1.9μm Titan-C18 particles provided a minimum reduced plate height as small as 1.7 for the most retained compound (n-octanophenone) under RPLC conditions. These particles are characterized by a relatively narrow size distribution with a relative standard deviation (RSD) of only 10%. A column packed with classical 5μm Symmetry-C18 particles, used as a reference RPLC column, generated a minimum reduced plate height of 2.1 for the same retained compound. This work demonstrates that this was due to an unusually low intra-particle diffusivity across these particles, which leads to a small longitudinal diffusion coefficient along the column. The demonstration is based on the combination of accurate measurements of the height equivalent to a theoretical plate (HETP), inverse size exclusion chromatography (ISEC), peak parking (PP), and minor disturbance method (MDM) experiments. The experimental results show that the reduced eddy dispersion HETP term (A=0.8 for a reduced velocity of 5), the internal particle porosity (ϵp=0.35), and the enrichment of acetonitrile in the pore volume (75% acetonitrile in the bulk, 85% inside the mesoporous volume) are identical on both the Titan-C18 and Symmetry-C18 columns. The difference between the internal structures of these two brands of RPLC-C18 fully porous particles lies in the values of the internal obstruction factor γp, which is 0.42 for the Symmetry-C18 but only 0.26 for the Titan-C18 particles. This is in part related to the diffusion hindrance due to the small average pore size of the Titan-C18 particles, around 59Å versus 77Å for Symmetry-C18 particles. A simple model of constriction along diffusion paths having the shape of a truncated cone suggests that the width of the pore size distribution (RSD of 30% and 20% for Titan-C18 and Symmetry-C18 particles) is mostly responsible for the difference in their obstruction factors.