Computational investigation of longitudinal diffusion, eddy dispersion, and trans-particle mass transfer in bulk, random packings of core–shell particles with varied shell thickness and shell diffusion coefficient

Harnessing an elastic flow instability to improve the kinetic performance of chromatographic columns

Despite decades of research and development, the optimal efficiency of slurry-packed HPLC columns is still hindered by inherent long-range flow heterogeneity from the wall to the central bulk region of these columns. Here, we show an example of how this issue can be addressed through the straightforward addition of a semidilute amount (500 ppm) of a large, flexible, synthetic polymer (18 MDa partially hydrolyzed polyacrylamide, HPAM) to the mobile phase (1% NaCl aqueous solution, hereafter referred to as "brine") during operation of a 4.6 mm × 300 mm column packed with 10μm BEHTM 125 Å particles. Addition of the polymer imparts elasticity to the mobile phase, causing the flow in the interparticle pore space to become unstable above a threshold flow rate. We verify the development of this elastic flow instability using pressure drop measurements of the friction factor versus Reynolds number. In prior work, we showed that this flow instability is characterized by large spatiotemporal fluctuations in the pore-scale flow velocities that may promote analyte dispersion across the column. Axial dispersion measurements of the quasi non-retained tracer thiourea confirm this possibility: they reveal that operating above the onset of the instability improves column efficiency by greater than 100%. These experiments thereby suggest that elastic flow instabilities can be harnessed to mitigate the negative impact of trans-column flow heterogeneities on the efficiency of slurry-packed HPLC columns. While this approach has its own inherent limitations and constraints, our results lay the groundwork for future targeted development of polymers that can impart elasticity when dissolved in commonly used liquid chromatography mobile phases, and can thereby generate elastic flow instabilities to help improve the resolution of HPLC columns.

Multiscale simulation of liquid chromatography: Effective diffusion in macro-mesoporous beds and the B-term of the plate height equation

We performed multiscale simulations of analyte sorption and diffusion in hierarchical porosity models of monolithic silica columns for reversed-phase liquid chromatography to investigate how the mean mesopore size of the chromatographic bed and the analyte-specific interaction with the chromatographic interface influence the analyte diffusivity at various length scales. The reproduced experimental conditions comprised the retention of six analyte compounds of low to moderate solute polarity on a silica-based, endcapped, C18 stationary phase with water‒acetonitrile and water-methanol mobile phases whose elution strength was varied via the volumetric solvent ratio. Detailed information about the analyte-specific interfacial dynamics received from molecular dynamics simulations was incorporated through appropriate linker schemes into Brownian dynamics diffusion simulations in three hierarchical porosity models received from physical reconstructions of silica monoliths with a mean macropore size of 1.23 µm and mean mesopore sizes of 12.3, 21.3, or 25.7 nm. The mean mesopore size was found to have a similar influence on the effective mesopore diffusivity as the analyte polarity and the mobile-phase elution strength, which together determine the analyte residence time on a column. A smaller mesopore size attenuated the increase of the effective mesopore diffusivity with increasing mobile-phase elution strength significantly. The effective bed diffusivity was limited by the analyte residence time rather than by morphological details of the mesopore space. The stronger an analyte was retained by the chromatographic interface inside the mesopores, the slower became the mass transfer between the pore space hierarchies and the lower was the effective bed diffusivity. The B-terms of the plate height equation were finally generated with the bed diffusivities and phase-based retention factors derived from the hierarchical porosity models using additional information about the stationary-phase limit obtained from the analysis of analyte-bonded phase contacts. The B-terms highlight analyte- and mobile phase-specific behavior relevant to isocratic and gradient elution conditions in chromatographic practice. In particular, U-shaped B-term curves are observed due to the dominating contribution of the retention factor and the bed diffusivity to the B-term at low and high elution strength of the mobile phase, respectively.

Multiple-open-tubular column enabling transverse diffusion. Part 3: Simulation of solute dispersion along a real three dimensional-printed column with quadratic channels

Multiple-open-tubular columns enabling transverse diffusion (MOTTD) consist of straight and parallel flow-through channels separated by a mesoporous stationary phase. In Part 1, a stochastic model of band broadening along MOTTD columns accounting for longitudinal diffusion, trans-channel velocity bias, and mass transfer resistance in the stationary phase was derived to demonstrate the intrinsic advantage of MOTTD columns over classical particulate columns. In Part 2, the model was refined for the critical contribution of the channel-to-channel polydispersity and applied to address the best trade-off between analysis speed and performance. In this Part 3, a MOTTD column with a square array of quadratic channels is fabricated by 3D-printing (combining polymer stereolithography with photolithography using photomasks) to deliver unprecedently small apparent channel diameters of 117.6 ± 5.0 μm. The colors in the microscopy photographs of the actual 3D-printed channels are binarized to delimitate the mobile phase volume from the stationary phase volume. The same numerical simulations as those in Part 2 are then performed for two MOTTD columns (external porosity ϵ_e=31.7%, same apparent channel diameter 117.6 μm): one containing 16 virtual perfect quadratic channels and the other 16 real 3D-printed channels. The reduced velocities (or Peclet numbers) are varied over a wide range from 0.2 to 5000 and the zone retention factors were fixed at k₁=1.04, 5, and 25. The results demonstrate that smoothing the edges of the targeted quadratic channels by the 3D-printed technique is advantageous in terms of solute dispersion. It outperforms the negative effect of the channel-to-channel polydispersity which is mitigated by transverse diffusion of the analyte in the stationary phase. For Peclet numbers larger than 50, the HETP of the 3D-printed MOTTD column is found 7%, 15%, and 16% smaller than that of the MOTTD column consisting of a square array of perfect quadratic channels. This confirms the known effect of channel geometry on solute dispersion in microfluidic systems. Flow channels in fabricated MOTTD columns are preferred to be circular so that the distribution of transverse diffusion lengths across the open channels remains as tight as possible. Finally, the general theory of nonuniform columns of Giddings reveals that the polydispersity of the cross-sectional area (RSD 8.4%) along a single 3D-printed channel has no negative impact on solute dispersion in MOTTD columns. Overall, MOTTD columns could become a serious alternative technology to conventional particulate columns. This implies a novel fabrication process that delivers circular channel diameters smaller than 10 μm, cross-sectional area polydispersity no larger than 25%, external porosities in a range from 15% (high speed separations) to 75% (high performance separations), and conventional mesoporous silica as the stationary phase. It adresses new synthesis routes based on either organic fibers or tubular micelle templating agents in suspension with silica gel solutions.

Quantitating Diffusion Enhancement in Pore Hierarchies

A microporous continuum traversed by a set of mutually perpendicular channels is considered to be a model for a hierarchically porous system of the mesoporous zeolite type. Transient profiles of molecular uptake as determined by kinetic Monte Carlo (kMC) simulation are found to be in excellent agreement with the result attained by the application of the two-region model (the Kärger model) of molecular diffusion. In particular, it is found that, in the two limiting cases referred to as fast exchange and slow exchange, there exist two simple analytical expressions for the rate of molecular uptake and hence for the quantification of transport enhancement in comparison with the purely microporous adsorbent. In the general case, transport enhancement is simply recognized by the reciprocal addition of the expressions in the two limiting cases.

Single-Molecule Dynamics Reflect IgG Conformational Changes Associated with Ion-Exchange Chromatography

Conformational changes of antibodies and other biologics can decrease the effectiveness of pharmaceutical separations. Hence, a detailed mechanistic picture of antibody-stationary phase interactions that occur during ion-exchange chromatography (IEX) can provide critical insights. This work examines antibody conformational changes and how they perturb antibody motion and affect ensemble elution profiles. We combine IEX, three-dimensional single-protein tracking, and circular dichroism spectroscopy to investigate conformational changes of a model antibody, immunoglobulin G (IgG), as it interacts with the stationary phase as a function of salt conditions. The results indicate that the absence of salt enhances electrostatic attraction between IgG and the stationary phase, promotes surface-induced unfolding, slows IgG motion, and decreases elution from the column. Our results reveal previously unreported details of antibody structural changes and their influence on macroscale elution profiles.

Multiple-open-tubular column enabling transverse diffusion. Part 2: Channel size distribution and structure optimization

Multiple-open-tubular columns enabling transverse diffusion (MOTTD) are made of straight, parallel, and cylindrical flow channels separated by a mesoporous stationary phase. In Part 1, a model of band broadening along MOTTD columns accounting for longitudinal diffusion, the trans-channel velocity bias, and mass transfer resistance in the stationary phase was proposed and validated. In this Part 2, the model is completed by considering the impact of short-range inter-channel velocity biases on the MOTTD plate number. These velocity biases are caused by the wide distribution of the channel diameters. Different ratios, ρ, of the average inner diameter, 2<r_c>, of the flow channels to their closest center-to-center distance d (d= 5 μm, ρ= 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, and 0.9) with a relative standard deviation (RSD) increasing from 0 to 50% are considered. The zone retention factor k₁ was increased from 1 to 25. The complete model of band broadening is validated after adjustment to dispersion data obtained by 1) the lattice-Boltzmann method for modeling fluid flow, 2) a random-walk particle-tracking (RWPT) technique to address advective-diffusive transport, and 3) by considering two distinct populations of flow channels (inner radii r_c,1=<r_c>(1-RSD) and r_c,2=<r_c>(1+RSD)) arranged at the nodes of a hexagonal compact array. The completed model of band broadening in MOTTD columns reveals that the RSD of the channel diameters has only a moderate impact on the optimum plate number of MOTTD columns: the relative increase of the minimum plate height do not exceed 30% even for the largest RSDs. However, when the mass transfer of the analyte is governed by its slow rate of transverse diffusion across the MOTTD column, the plate height can be increased by up to 100% at high average velocities. Regarding the best trade-off between analysis speed and column performance at a fixed pressure drop of 400 bar, irrespective of the zone retention factor and RSD of the distribution of the channel diameters, the fastest analyses are recommended for MOTTD columns having a small structural parameter ρ. In contrast, for the longest analysis times, the largest values of ρ are required to maximize the performance of MOTTD columns.

The theory and practice of ultrafast liquid chromatography: A tutorial

Modern high-throughput experimentation and challenging analytical problems of academic/industrial research have put the responsibility on separation scientists to develop new fast separation approaches. With the availability of high-pressure pumps, small particles with hydrolytically stable surface chemistries, reduced extra-column band broadening, and low volume detectors with fast signal processing, it is now feasible to do sub-minute to sub-second chromatography. Herein, the fundamental theoretical principles of ultrafast chromatography, along with practical solutions, are reviewed. Approaches for rapid separations in packed beds, narrow open tubular columns, and monoliths are demonstrated, along with the challenges that were faced. The instrumentation requirements (pumps, injection systems, detectors, column packing process) for using short columns ranging from 0.5 to 5 cm are examined, followed by real applications. One of the main problems in ultrafast chromatography is partial or complete peak overlap. As per Gidding's statistical overlap theory, peak overlap cannot be avoided for a completely random sample for a column with a given peak capacity. Signal processing techniques based on Fourier transform deconvolution of band broadening, power laws, derivatives, and iterative curve fitting are explained to help improve the chromatographic resolution. An example of ten peaks separated in under a second is shown and discussed. Other ultrafast separations in supercritical fluid chromatography or capillary electrophoresis are briefly mentioned to provide a complete understanding of this emerging field.

Shedding light on mechanisms leading to convex-upward van Deemter curves on a cellulose tris(4-chloro-3-methylphenylcarbamate)-based chiral stationary phase

An unusual convex-upward van Deemter curve was observed for the more retained enantiomer of a chiral sulfoxide (2-(benzylsulfinyl)benzamide) on a cellulose tris(4-chloro-3-methylphenylcarbamate)-based chiral stationary phase (CSP), prepared on silica particles of 1000 Å pore size. In contrast, the firstly eluted enantiomer of the same molecule exhibited the traditional convex-downward van Deemter curve. A detailed kinetic and thermodynamic investigation has revealed that this unusual phenomenon, which however has already been observed in chiral chromatography, originates when the adsorption of the compound is very strong and the solid-phase diffusion negligible. Experimentally, the intraparticle diffusion of the more retained enantiomer of the sulfoxide was found to be one order of magnitude smaller than that of the first eluted one. Overall, this translates into very little longitudinal diffusion (b-term of van Deemter curve) accompanied by high solid-liquid mass transfer resistance (c-term). Finally the comparison with another, differently-substituted chiral sulfoxide (whose enantiomers both exhibit traditional van Deemter curve behavior) has allowed to correlate these findings to the specific characteristics of the molecule.

Morphology-transport relationships for SBA-15 and KIT-6 ordered mesoporous silicas

Quantitative morphology-transport relationships are derived for ordered mesoporous silicas through direct numerical simulation of hindered diffusion in realistic geometrical models of the pore space obtained from physical reconstruction by electron tomography. We monitor accessible porosity and effective diffusion coefficients resulting from steric and hydrodynamic interactions between passive tracers and the pore space confinement as a function of λ = d_tracer/d_meso (ratio of tracer diameter to mean mesopore diameter) in SBA-15 (d_meso = 9.1 nm) and KIT-6 (d_meso = 10.5 nm) silica samples. For λ = 0, the pointlike tracers reproduce the true diffusive tortuosities. For 0 ≤λ < 0.5, the derived hindrance factor quantifies the extent to which diffusion of finite-size tracers through the materials is hindered compared with free diffusion in the bulk liquid. The hindrance factor connects the transport properties of the ordered silicas to their mesopore space morphologies and enables quantitative comparison with random mesoporous silicas. Key feature of the ordered silicas is a narrow, symmetric mesopore size distribution (∼10% relative standard deviation), which engenders a sharper decline of the accessible-porosity window with increasing λ than observed for random silicas with their wide, asymmetric mesopore size distributions. As support structures, ordered mesoporous silicas should offer benefits for applications where spatial confinement effects and molecular size-selectivity are of prime importance. On the other hand, random mesoporous silicas enable higher diffusivities for λ > 0.3, because the larger pores carry most of the diffusive flux and keep pathways open when smaller pores have closed off.

Recent advances in capillary ultrahigh pressure liquid chromatography

In the twenty years since its initial demonstration, capillary ultrahigh pressure liquid chromatography (UHPLC) has proven to be one of most powerful separation techniques for the analysis of complex mixtures. This review focuses on the most recent advances made since 2010 towards increasing the performance of such separations. Improvements in capillary column preparation techniques that have led to columns with unprecedented performance are described. New stationary phases and phase supports that have been reported over the past decade are detailed, with a focus on their use in capillary formats. A discussion on the instrument developments that have been required to ensure that extra-column effects do not diminish the intrinsic efficiency of these columns during analysis is also included. Finally, the impact of these capillary UHPLC topics on the field of proteomics and ways in which capillary UHPLC may continue to be applied to the separation of complex samples are addressed.

Bed morphological features associated with an optimal slurry concentration for reproducible preparation of efficient capillary ultrahigh pressure liquid chromatography columns

Column wall effects and the formation of larger voids in the bed during column packing are factors limiting the achievement of highly efficient columns. Systematic variation of packing conditions, combined with three-dimensional bed reconstruction and detailed morphological analysis of column beds, provide valuable insights into the packing process. Here, we study a set of sixteen 75μm i.d. fused-silica capillary columns packed with 1.9μm, C18-modified, bridged-ethyl hybrid silica particles slurried in acetone to concentrations ranging from 5 to 200mg/mL. Bed reconstructions for three of these columns (representing low, optimal, and high slurry concentrations), based on confocal laser scanning microscopy, reveal morphological features associated with the implemented slurry concentration, that lead to differences in column efficiency. At a low slurry concentration, the bed microstructure includes systematic radial heterogeneities such as particle size-segregation and local deviations from bulk packing density near the wall. These effects are suppressed (or at least reduced) with higher slurry concentrations. Concomitantly, larger voids (relative to the mean particle diameter) begin to form in the packing and increase in size and number with the slurry concentration. The most efficient columns are packed at slurry concentrations that balance these counteracting effects. Videos are taken at low and high slurry concentration to elucidate the bed formation process. At low slurry concentrations, particles arrive and settle individually, allowing for rearrangements. At high slurry concentrations, they arrive and pack as large patches (reflecting particle aggregation in the slurry). These processes are discussed with respect to column packing, chromatographic performance, and bed microstructure to help reinforce general trends previously described. Conclusions based on this comprehensive analysis guide us towards further improvement of the packing process.

Larger voids in mechanically stable, loose packings of 1.3μm frictional, cohesive particles: Their reconstruction, statistical analysis, and impact on separation efficiency

Lateral transcolumn heterogeneities and the presence of larger voids in a packing (comparable to the particle size) can limit the preparation of efficient chromatographic columns. Optimizing and understanding the packing process provides keys to better packing structures and column performance. Here, we investigate the slurry-packing process for a set of capillary columns packed with C18-modified, 1.3μm bridged-ethyl hybrid porous silica particles. The slurry concentration used for packing 75μm i.d. fused-silica capillaries was increased gradually from 5 to 50mg/mL. An intermediate concentration (20mg/mL) resulted in the best separation efficiency. Three capillaries from the set representing low, intermediate, and high slurry concentrations were further used for three-dimensional bed reconstruction by confocal laser scanning microscopy and morphological analysis of the bed structure. Previous studies suggest increased slurry concentrations will result in higher column efficiency due to the suppression of transcolumn bed heterogeneities, but only up to a critical concentration. Too concentrated slurries favour the formation of larger packing voids (reaching the size of the average particle diameter). Especially large voids, which can accommodate particles from>90% of the particle size distribution, are responsible for a decrease in column efficiency at high slurry concentrations. Our work illuminates the increasing difficulty of achieving high bed densities with small, frictional, cohesive particles. As particle size decreases interparticle forces become increasingly important and hinder the ease of particle sliding during column packing. While an optimal slurry concentration is identified with respect to bed morphology and separation efficiency under conditions in this work, our results suggest adjustments of this concentration are required with regard to particle size, surface roughness, column dimensions, slurry liquid, and external effects utilized during the packing process (pressure protocol, ultrasound, electric fields).