Brownian Sieving Effect for Boosting the Performance of Microcapillary Hydrodynamic Chromatography. Proof of Concept

An alternative general model for the effective longitudinal diffusion in chromatographic beds filled with ordered porous particles

The two-zone moment-analysis method for the determination of the dispersion tensor in hierarchical retentive porous media has been adopted to compute and model the effective longitudinal diffusion Deff, or equivalently the B-term band broadening, in chromatographic beds filled with ordered porous particles. On the one hand, this approach offers accurate numerical results for Deff while keeping computational expenses low. On the other hand, it also gives direct insight for the analytical modelling, readily revealings the two main essential quantities (resp. referred to as the mobile-zone and stationary-zone effective diffusion factors γm and γs) that contribute to Deff. Modelling these two main parameters provided us with two new analytical models for Deff: a general one, valid for diluted and concentrated packings and accurate in the whole range of relevant intra-particle diffusion coefficient Dpz, and an approximate one, reliable for diluted packings and accurate also for concentrated packings with low to intermediate values of Dpz. The large advantage of both models is that they do not need any fitting parameter because all the required information is incorporated into the experimentally accessible geometric obstruction factor in the mobile phase originating from the tortuosity of the through-pore space (limiting case of fully solid particles without any retention). These models hence serve as an alternative to the Effective Medium Theory (EMT) models used so far in the literature. To validate the theory, five ordered geometries have been investigated. The accuracy of the general model proposed has been quantified and found to be comparable with that of the 3rd order approximate Torquato model for four geometries, even for macro-porosities close to the close-packing limit. The case of a 2-d triangular array of ellipsoidal particles with different elongations is also investigated to show the general validity and applicability of the models.

Moment analysis for predicting effective transport properties in hierarchical retentive porous media

We report on a new homogenization approach to solve, with drastically improved speed and accuracy, the general advection-diffusion equation in hierarchical porous media with localized diffusion and adsorption/desorption processes, thus opening the way to a much deeper understanding of the band broadening process in chromatographic systems. The proposed robust and efficient moment-based approach allows us to compute the exact local and integral concentration moments and hence provides exact solutions for the effective velocity and dispersion coefficients of migrating solute particles. Innovative to the proposed method is also that it not only produces the exact effective transport parameters of the long-time asymptotic solution, but also their entire transient. The analysis of the transient behaviour can be used, for example, to properly identify the time and length scales needed to achieve the macro-transport conditions. If the hierarchical porous media can be represented as the periodic repetition of a unit lattice cell, the method only requires the solution of the time-dependent advection-diffusion equations for the zeroth order and first-order exact local moments, exclusively on the unit cell. This implies an enormous reduction of the computational efforts and a significant improvement of the accuracy of the results when compared to the direct numerical simulation (DNS) approaches which require flow domains that are long enough to achieve steady-state conditions, and hence often cover tens to hundreds of unit cells. The reliability of the proposed method is verified by comparing its predictions with DNS results, in one, two and three dimensions, in both transient and asymptotic conditions. The influence of top and bottom no-slip walls on the separation performance of chromatographic columns with micromachined porous and nonporous pillars is discussed in detail.

Parallel computing for mobilities in periodic geometries

We examine methods for calculating the effective mobilities of molecules driven through periodic geometries in the context of particle-based simulation. The standard formulation of the mobility, based on the long-time limit of the mean drift velocity, is compared to a formulation based on the mean first-passage time of molecules crossing a single period of the system geometry. The equivalence of the two definitions is derived under weaker assumptions than similar conclusions obtained previously, requiring only that the state of the system at subsequent period crossings satisfy the Markov property. Approximate theoretical analyses of the computational costs of estimating these two mobility formulations via particle simulations suggest that the definition based on first-passage times may be substantially better suited to exploiting parallel computation hardware. This claim is investigated numerically on an example system modeling the passage of nanoparticles through the slit-well device. In this case, the traditional mobility formulation is found to perform best when the Péclet number is small, whereas the mean first-passage time formulation is found to converge much more quickly when the Péclet number is moderate or large. The results suggest that, given relatively modest access to modern GPU hardware, this alternative mobility formulation may be an order of magnitude faster than the standard technique for computing effective mobilities of biomolecules through periodic geometries.

50-Fold Reduction of Separation Time in Open-Channel Hydrodynamic Chromatography via Lateral Vortices

Despite its relatively long history, open-channel hydrodynamic chromatography (OC-HDC) still represents a niche technique for determining the size distribution of particle suspensions. Practical limitations of this separation method ultimately arise from the low eluent velocity that is necessary to contain the adverse increase of analyte bandwidth caused by Taylor-Aris dispersion. Because of the micrometric size of the channel cross section, the low eluent velocity translates into order of pL-per-minute flow rates, which introduce a challenge for both the injection and the detection systems. In this article, we provide theoretical/numerical evidence illustrating how a sizable reduction of the Taylor-Aris effect can be obtained by triggering cross-sectional vortices alongside the main pressure-driven axial flow. As a case study, we consider a square channel geometry where the lateral vortices are created by DC-induced electroosmosis. The analysis of particle separation is based on the classical excluded-volume macrotransport approach, which allows to derive the average particle velocity and the axial dispersion coefficient from the solution of two stationary advection-diffusion problems defined onto the channel cross section. We find that lateral vortices can enhance the separation efficiency quantitatively, e.g., by reducing the separation time of a two-species mixture by a 50-fold factor compared to standard OC-HDC.

Brownian sieving enhancement of microcapillary hydrodynamic chromatography. Analysis of the separation performance based on Brenner's macro-transport theory

In a recent article [Analytical Chemistry, 93(17), 6808-6816 (2021)], an unconventional device configuration enforcing a Brownian sieving mechanism was proposed as proof of concept for the efficient implementation of microcapillary hydrodynamic chromatography (MHDC). In this article, we perform a thorough analysis of the device geometry and of operating conditions, in order to single out the optimal configuration maximizing separation resolution. Brenner's macro-transport theory provides the technical picklock to perform the search for the optimum over a wide choice of device geometries and for a range of values of the particle Péclet number covering most conditions encountered in practical implementations of MHDC. Specifically, effective transport coefficients defining the dynamics of the suspended phase are obtained by the solution of a two-dimensional steady-state advection-diffusion equation defined onto the channel cross-section. The eigenvalue/eigenfunction structure of the associated transient problem is exploited in order to quantify the timescale for reaching the macro-transport regime conditions. Based on this timescale and on the effective transport parameters, an estimate of the column length necessary to achieve a prescribed level of separation resolution is obtained. We identify device geometry and operating conditions where the capillary length is shrunk down by a factor above ten compared to the standard MHDC configuration. Lagrangian stochastic statistics of particle ensembles are used to validate the results obtained through Brenner's macro-transport approach. The method proposed can be readily generalized to other classes of device geometries enforcing the same Brownian sieving mechanism.