A neural network model for rapid prediction of analyte focusing in isotachophoresis

We present the development and demonstration of a neural network (NN) model for fast and accurate prediction of whether or not a chosen analyte is focused by an isotachophoresis (ITP) buffer system. The NN model is useful in the rapid evaluation of possible ITP chemistries applicable to analytes of interest. We trained and tested the NN model for univalent species based on extensive data sets of over 10,000 anionic and 10,000 cationic ITP simulations. The NN model uses as inputs the mobilities and the acid dissociation constants of leading electrolyte ion, trailing electrolyte ion, counterion, and a single analyte as well as the leading-to-counterion concentration ratio of the leading zone. The output then indicates whether the chosen electrolyte system yields stable ITP focusing of the analyte. The prediction accuracy of the NN model is over 97.7%. We demonstrate the applicability of the NN by validating its predictions with reported experimental data for anionic and cationic ITP. We have packaged the NN model in a free, web-based application named IONN (isotachophoresis on neural network), which can be used to rapidly screen ITP electrolyte systems.

A review of the zone broadening contributions in free-flow electrophoresis

The broadening of analyte streams, as they migrate through a free-flow electrophoresis (FFE) channel, often limits the resolving power of FFE separations. Under laminar flow conditions, such zonal spreading occurs due to analyte diffusion perpendicular to the direction of streamflow and variations in the lateral distance electrokinetically migrated by the analyte molecules. Although some of the factors that give rise to these contributions are inherent to the FFE method, others originate from non-idealities in the system, such as Joule heating, pressure-driven crossflows, and a difference between the electrical conductivities of the sample stream and background electrolyte. The injection process can further increase the stream width in FFE separations but normally influencing all analyte zones to an equal extent. Recently, several experimental and theoretical works have been reported that thoroughly investigate the various contributions to stream variance in an FFE device for better understanding, and potentially minimizing their magnitudes. In this review article, we carefully examine the findings from these studies and discuss areas in which more work is needed to advance our comprehension of the zone broadening contributions in FFE assays.

Evolution of the theoretical description of the isoelectric focusing experiment: II. An open system isoelectric focusing experiment is a transient, bidirectional isotachophoretic experiment

The carrier ampholytes-based (CA-based) isoelectric focusing (IEF) experiment evolved from Svensson's closed system IEF (constant spatial current density, absence of convective mixing, counter-balancing electrophoretic and diffusive fluxes yielding a steady state pH gradient) to the contemporary open system IEF (absence of convective mixing, large cross-sectional area electrode vessels, lack of counter-balancing electrophoretic- and diffusive fluxes leading to transient pH gradients). Open system IEF currently is described by a two-stage model: In the first stage, a rapid IEF process forms the pH gradient which, in the second stage, is slowly degraded by isotachophoretic processes that move the most acidic and most basic CAs into the electrode vessels. An analysis of the effective mobilities and the effective mobility to conductivity ratios of the anolyte, catholyte, and the CAs indicates that in open system IEF experiments a single process, transient bidirectional isotachophoresis (tbdITP) operates from the moment current is turned on until it is turned off. In tbdITP, the anolyte and catholyte provide the leading ions and the pI 7 CA or the reactive boundary of the counter-migrating H₃ O⁺ and OH^- ions serves as the shared terminator. The outcome of the tbdITP process is determined by the ionic mobilities, pK_a values, and loaded amounts of all ionic and ionizable components: It is constrained by both the transmitted amount of charge and the migration space available for the leading ions. tbdITP and the resulting pH gradient can never reach steady state with respect to the spatial coordinate of the separation channel.

Isotachophoresis: Theory and Microfluidic Applications

Isotachophoresis (ITP) is a versatile electrophoretic technique that can be used for sample preconcentration, separation, purification, and mixing, and to control and accelerate chemical reactions. Although the basic technique is nearly a century old and widely used, there is a persistent need for an easily approachable, succinct, and rigorous review of ITP theory and analysis. This is important because the interest and adoption of the technique has grown over the last two decades, especially with its implementation in microfluidics and integration with on-chip chemical and biochemical assays. We here provide a review of ITP theory starting from physicochemical first-principles, including conservation of species, conservation of current, approximation of charge neutrality, pH equilibrium of weak electrolytes, and so-called regulating functions that govern transport dynamics, with a strong emphasis on steady and unsteady transport. We combine these generally applicable (to all types of ITP) theoretical discussions with applications of ITP in the field of microfluidic systems, particularly on-chip biochemical analyses. Our discussion includes principles that govern the ITP focusing of weak and strong electrolytes; ITP dynamics in peak and plateau modes; a review of simulation tools, experimental tools, and detection methods; applications of ITP for on-chip separations and trace analyte manipulation; and design considerations and challenges for microfluidic ITP systems. We conclude with remarks on possible future research directions. The intent of this review is to help make ITP analysis and design principles more accessible to the scientific and engineering communities and to provide a rigorous basis for the increased adoption of ITP in microfluidics.

Electrophoresis simulations using Chebyshev pseudo-spectral method on a moving mesh

We present the implementation and demonstration of the Chebyshev pseudo-spectral method coupled with an adaptive mesh method for performing fast and highly accurate electrophoresis simulations. The Chebyshev pseudo-spectral method offers higher numerical accuracy than all other finite difference methods and is applicable for simulating all electrophoresis techniques in channels with open or closed boundaries. To improve the computational efficiency, we use a novel moving mesh scheme that clusters the grid points in the regions with poor numerical resolution. We demonstrate the application of the Chebyshev pseudo-spectral method on a moving mesh for simulating nonlinear electrophoretic processes through examples of isotachophoresis (ITP), isoelectric focusing (IEF), and electromigration-dispersion in capillary zone electrophoresis (CZE) at current densities as high as 1000 A/m 2 . We also show the efficacy of our moving mesh method over existing methods that cluster the grid points in the regions with large concentration gradients. We have integrated the adaptive Chebyshev pseudo-spectral method in the open-source SPYCE simulator and verified its implementation with other electrophoresis simulators.

Web-Based Open-Source Tool for Isotachophoresis

We present the development of a client-side web-based simulator for complex electrophoresis phenomena, including isotachophoresis. The simulation tool is called Client-based Application for Fast Electrophoresis Simulation (CAFES). CAFES uses the broad cross-browser compatibility of JavaScript to provide a rapid and easy-to-use tool for coupled unsteady electromigration, diffusion, and equilibrium electrolyte reactions among multiple weak electrolytes. The code uses a stationary grid (for simplicity) and an adaptive time step to provide reliable estimates of ion concentration dynamics (including pH profile evolution), requiring no prior installation nor compilation. CAFES also offers a large database of commonly used species and their relevant physicochemical properties. We present a validation of predictions from CAFES by comparing them to experimental data of peak- and plateau-mode isotachophoresis experiments. The code yields accurate estimates of interface velocity, plateau length and relative intensity, and pH variations while significantly reducing the computation time compared to existing codes. The tool is open-source and available for free at https://microfluidics.stanford.edu/cafes.

Microscale electrokinetic-based analysis of intact cells and viruses

Miniaturized electrokinetic methods have proven to be robust platforms for the analysis and assessment of intact microorganisms, offering short response times and higher integration than their bench-scale counterparts. The present review article discusses three types of electrokinetic-based methodologies: electromigration or motion-based techniques, electrode-based electrokinetics, and insulator-based electrokinetics. The fundamentals of each type of methodology are discussed and relevant examples from recent reports are examined, to provide the reader with an overview of the state-of-the-art on the latest advancements on the analysis of intact cells and viruses with microscale electrokinetic techniques. The concluding remarks discuss the potential applications and future directions.

Dynamic computer simulations of electrophoresis: 2010-2020

The transport of components in liquid media under the influence of an applied electric field can be described with the continuity equation. It represents a nonlinear conservation law that is based upon the balance laws of continuous transport processes and can be solved in time and space numerically. This procedure is referred to as dynamic computer simulation. Since its inception four decades ago, the state of dynamic computer simulation software and its use has progressed significantly. Dynamic models are the most versatile tools to explore the fundamentals of electrokinetic separations and provide insights into the behavior of buffer systems and sample components of all electrophoretic separation methods, including moving boundary electrophoresis, CZE, CGE, ITP, IEF, EKC, ACE, and CEC. This article is a continuation of previous reviews (Electrophoresis 2009, 30, S16–S26 and Electrophoresis 2010, 31, 726–754) and summarizes the progress and achievements made during the 2010 to 2020 time period in which some of the existing dynamic simulators were extended and new simulation packages were developed. This review presents the basics and extensions of the three most used one‐dimensional simulators, provides a survey of new one‐dimensional simulators, outlines an overview of multi‐dimensional models, and mentions models that were briefly reported in the literature. A comprehensive discussion of simulation applications and achievements of the 2010 to 2020 time period is also included.

Simul 6: A fast dynamic simulator of electromigration

Simul 6 is a 1D dynamic simulator of electromigration based on the mathematical model of electromigration in free solutions. The model consists of continuity equations for the movement of electrolytes in a separation channel, acid-base equilibria of weak electrolytes, and the electroneutrality condition. It accounts for any number of multivalent electrolytes or ampholytes and provides a complete picture about dynamics of electromigration and diffusion in the separation channel. The equations are solved numerically using software means which allow for parallelization and multithreaded computation. Simul 6 has a user-friendly graphical interface. It is typically used for inspection of system peaks (zones) in electrophoresis, stacking and preconcentrating analytes, optimization of separation conditions, method development in either capillary zone electrophoresis, isotachophoresis, and isoelectric focusing. Simul 6 is the successor of Simul 5, and has been launched as a free software available for download at https://simul6.app/.

High-resolution numerical simulations of electrophoresis using the Fourier pseudo-spectral method

We present the formulation, implementation, and performance evaluation of the Fourier pseudo-spectral method for performing fast and accurate simulations of electrophoresis. We demonstrate the applicability of this method for simulating a wide variety of electrophoretic processes such as capillary zone electrophoresis, transient-isotachophoresis, field amplified sample stacking, and oscillating electrolytes. Through these simulations, we show that the Fourier pseudo-spectral method yields accurate and stable solutions on coarser computational grids compared with other nondissipative spatial discretization schemes. Moreover, due to the use of coarser grids, the Fourier pseudo-spectral method requires lower computational time to achieve the same degree of accuracy. We have demonstrated the application of the Fourier pseudo-spectral method for simulating realistic electrophoresis problems with current densities as high as 5000 A/m² with over tenfold speed-up compared to the commonly used second-order central difference scheme, to achieve a given degree of accuracy. The Fourier pseudo-spectral method is also suitable for simulating electrophoretic processes involving a large number of concentration gradients, which render the adaptive grid-refinement techniques ineffective. We have integrated the numerical scheme in a new electrophoresis simulator named SPYCE, which we offer to the community as open-source code.

Simulation and experiment of asymmetric electrode placement for electrophoretic exclusion in a microdevice

Electrophoretic exclusion (EE) is a counterflow gradient technique that exploits hydrodynamic flow and electrophoretic forces to exclude, enrich, and separate analytes. Resolution for this technique has been theoretically examined and the smallest difference in electrophoretic mobilities that can be completely separated is estimated to be 10^-13  cm² /Vs. Traditional and mesoscale systems have been used, whereas microfluidics offers a greater range of geometries and configurations towards approaching this theoretical limit. To begin to understand the impact of seemingly subtle changes to the entrance flow and the electric field configurations, three closely related microfluidic interfaces were modeled, fabricated, and tested. These interfaces consisted of systematically varying placement of an asymmetric electrode relative to a channel entrance: leading electrode placed outside the channel entrance, leading electrode aligned with the channel, and leading electrode placed within the channel. A charged fluorescent dye is used as a sensitive and accurate probe for the model and to test the concentration variation at these interfaces. Models and experiments focused on visualizing the concentration profile in areas of high temporal dynamics, thus providing a severe test of the models. Experimental data and simulation results showed strong qualitative agreement. The complexity of the electric and flow fields about this interface and the agreement between models and testing suggests the theoretical assessment capabilities can be used to faithfully design novel and more efficient interfaces.

Electrohydrodynamic instability of ion-concentration shock wave in electrophoresis

Capillary electrophoresis techniques often involve ion-concentration shock waves in an electrolyte solution, propagating under the effect of an external electric field. These shock waves are characterized by self-sharpening gradients in ion concentrations and electrical conductivity that are collinear with the electric field. The coupling of electric field and fluid motion at the shock interface sometimes leads to an undesirable electrohydrodynamic (EHD) instability. Using linear stability analysis, we describe the motion of small-amplitude disturbances of an electrophoretic shock wave. Our analysis shows that the EHD instability results due to the competition between destabilizing electroviscous flow and stabilizing electromigration of the shock wave. The ratio of timescales corresponding to electroviscous flow and electromigration yields a threshold criterion for the onset of instability. We present a validation of this threshold criterion with published experimental data and also describe the physical mechanism underlying the EHD instability of the electrophoretic shock wave.

Stability of electrophoretic transport of ions

We present an investigation of instability during electrophoretic transport of ions in a class of electrolytes called oscillating electrolytes. We analyze the onset of instability in electrophoretic transport in a binary electrolyte by modeling growth of small concentration disturbances over a base state with uniform acid and base concentrations. Our linear stability analysis shows that the growth rate of low wave-number concentration disturbances increases with an increase in wave number. Whereas, the growth rate of high wave-number disturbances decreases with increasing wave number due to the stabilizing effect of molecular diffusion. Our analysis also yields the scaling for growth rates and the wave number of most unstable mode with electric field. In addition, we show that the electrophoretic system exhibits instability only for a certain range of species concentrations. We also discuss the physical mechanism underlying the instability of transport process. We show that the instability is exhibited by those binary electrolytes that consist of a multivalent species with unusually high electrophoretic mobility in higher ionization states. Throughout, we provide verification of our linear stability analysis with full nonlinear simulations.

A method for non-invasive full-field imaging and quantification of chemical species

We present a novel method for full-field scalar visualization and quantification of species concentration fields. We term this method species-altered fluorescence imaging (SAFI). The method employs electrically neutral fluorescent dyes whose quantum yields are selectively quenched or enhanced by species of interest. SAFI enables simultaneous imaging of material interfaces and provides non-invasive, scalar-field quantitation of two-dimensional species concentration fields. We describe criteria for choosing SAFI dyes and tabulate 35 promising SAFI dyes and their relevant properties. Next, we describe species concentration quantification with SAFI via Stern-Volmer quenching and discuss the sensitivity and resolution of our method. We demonstrate this method with two dyes, 6-methoxy-N-(3-sulfopropyl)quinolinium (SPQ) and 10-(3-sulfopropyl)acridinium betaine (SAB). We demonstrate our method in full-field visualization of several challenging electrokinetic flows: isotachophoresis (ITP) in both cationic and anionic modes, and in a convective electrokinetic instability (EKI) flow. Through these experiments we collectively quantify ion concentration shock velocities, simultaneously measure concentrations of five species, and quantify the development of an unsteady, chaotic, 2D flow.