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.

Analytical solutions for viscoelectric effects in electrokinetic nanochannels

Understanding electrokinetic transport in nanochannels and nanopores is essential for emerging biological and electrochemical applications. The viscoelectric effect is an important mechanism implicated in the increase of local viscosity due to the polarization of a solvent under a strong electric field. However, most analyses of the viscoelectric effect have been limited to numerical analyses. In this work, we present a set of analytical solutions applicable to the physical description of viscoelectric effects in nanochannel electrokinetic systems. To achieve such closed-form solutions, we employ the Debye-Hückel approximation of small diffuse charge layer potentials compared to the thermal potential. We analyze critical parameters, including electroosmotic flow profiles, electroosmotic mobility, flow rate, and channel conductance. We compare and benchmark our analytical solutions with published predictions from numerical models. Importantly, we leverage these analytical solutions to identify essential thermophysical and nondimensional parameters that govern the behavior of these systems. We identify scaling parameters and relations among surface charge density, ionic strength, and nanochannel height.

A critical review of microfluidic systems for CRISPR assays

Reviewed are nucleic acid detection assays that incorporate clustered regularly interspaced short palindromic repeats (CRISPR)-based diagnostics and microfluidic devices and techniques. The review serves as a reference for researchers who wish to use CRISPR-Cas systems for diagnostics in microfluidic devices. The review is organized in sections reflecting a basic five-step workflow common to most CRISPR-based assays. These steps are analyte extraction, pre-amplification, target recognition, transduction, and detection. The systems described include custom microfluidic chips and custom (benchtop) chip control devices for automated assays steps. Also included are partition formats for digital assays and lateral flow biosensors as a readout modality. CRISPR-based, microfluidics-driven assays offer highly specific detection and are compatible with parallel, combinatorial implementation. They are highly reconfigurable, and assays are compatible with isothermal and even room temperature operation. A major drawback of these assays is the fact that reports of kinetic rates of these enzymes have been highly inconsistent (many demonstrably erroneous), and the low kinetic rate activity of these enzymes limits achievable sensitivity without pre-amplification. Further, the current state-of-the-art of CRISPR assays is such that nearly all systems rely on off-chip assays steps, particularly off-chip sample preparation.

Liquid Heterostructures: Generation of Liquid-Liquid Interfaces in Free-Flowing Liquid Sheets

Chemical reactions and biological processes are frequently governed by the structure and dynamics of the interface between two liquid phases, but these interfaces are often difficult to study due to the relative abundance of the bulk liquids. Here, we demonstrate a method for generating multilayer thin film stacks of liquids, which we call liquid heterostructures. These free-flowing layered liquid sheets are produced with a microfluidic nozzle that impinges two converging jets of one liquid onto opposite sides of a third jet of another liquid. The resulting sheet consists of two layers of the first liquid enveloping an inner layer of the second liquid. Infrared microscopy, white light reflectivity, and imaging ellipsometry measurements demonstrate that the buried liquid layer has a tunable thickness and displays well-defined liquid-liquid interfaces and that this inner layer can be only tens of nanometers thick. The demonstrated multilayer liquid sheets minimize the amount of bulk liquid relative to their buried interfaces, which makes them ideal targets for spectroscopy and scattering experiments.

Detection and Discrimination of Single Nucleotide Polymorphisms by Quantification of CRISPR-Cas Catalytic Efficiency

The specificity of CRISPR-Cas12 assays is attractive for the detection of single nucleotide polymorphisms (SNPs) implicated in, e.g., cancer and SARS-CoV-2 variants. Such assays often employ endpoint measurements of SNP or wild type (WT) activated Cas12 trans-cleavage activity; however, the fundamental kinetic effects of SNP versus WT activation remain unknown. We here show that endpoint-based assays are limited by arbitrary experimental choices (like used reporter concentration and assay duration) and work best for known target concentrations. More importantly, we show that SNP (versus WT) activation results in measurable kinetic shifts in the Cas12 trans-cleavage substrate affinity (K_M) and apparent catalytic efficiency (k_cat^*/K_M). To address endpoint-based assay limitations, we then develop an assay based on the quantification of Michaelis-Menten parameters and apply this assay to a 20 base pair WT target of the SARS-CoV-2 E gene. We find that the k_cat^*/K_M measured for WT is 130-fold greater than the lowest k_cat^*/K_M among all 60 measured SNPs (compared to a 4.8-fold for endpoint fluorescence of the same SNP). K_M also offers a strong ability to distinguish SNPs, varies 27-fold over all the cases, and, importantly, is insensitive to the target concentration. Last, we point out trends among kinetic rates and SNP base and location within the CRISPR-Cas12 targeted region.

Electrochemical Methods for Water Purification, Ion Separations, and Energy Conversion

Agricultural development, extensive industrialization, and rapid growth of the global population have inadvertently been accompanied by environmental pollution. Water pollution is exacerbated by the decreasing ability of traditional treatment methods to comply with tightening environmental standards. This review provides a comprehensive description of the principles and applications of electrochemical methods for water purification, ion separations, and energy conversion. Electrochemical methods have attractive features such as compact size, chemical selectivity, broad applicability, and reduced generation of secondary waste. Perhaps the greatest advantage of electrochemical methods, however, is that they remove contaminants directly from the water, while other technologies extract the water from the contaminants, which enables efficient removal of trace pollutants. The review begins with an overview of conventional electrochemical methods, which drive chemical or physical transformations via Faradaic reactions at electrodes, and proceeds to a detailed examination of the two primary mechanisms by which contaminants are separated in nondestructive electrochemical processes, namely electrokinetics and electrosorption. In these sections, special attention is given to emerging methods, such as shock electrodialysis and Faradaic electrosorption. Given the importance of generating clean, renewable energy, which may sometimes be combined with water purification, the review also discusses inverse methods of electrochemical energy conversion based on reverse electrosorption, electrowetting, and electrokinetic phenomena. The review concludes with a discussion of technology comparisons, remaining challenges, and potential innovations for the field such as process intensification and technoeconomic optimization.

Enzyme Kinetics and Detector Sensitivity Determine Limits of Detection of Amplification-Free CRISPR-Cas12 and CRISPR-Cas13 Diagnostics

Interest in CRISPR-Cas12 and CRISPR-Cas13 detection continues to increase as these detection schemes enable the specific recognition of nucleic acids. The fundamental sensitivity limits of these schemes (and their applicability in amplification-free assays) are governed by kinetic rates. However, these kinetic rates remain poorly understood, and their reporting has been inconsistent. We quantify kinetic parameters for several enzymes (LbCas12a, AsCas12a, AapCas12b, LwaCas13a, and LbuCas13a) and their corresponding limits of detection (LoD). Collectively, we present quantification of enzyme kinetics for 14 guide RNAs (gRNAs) and nucleic acid targets for a total of 50 sets of kinetic rate parameters and 25 LoDs. We validate the self-consistency of our measurements by comparing trends and limiting behaviors with a Michaelis-Menten trans-cleavage reaction kinetics model. For our assay conditions, activated Cas12 and Cas13 enzymes exhibit trans-cleavage catalytic efficiencies between order 105 and 106 M-1 s-1. For assays that use fluorescent reporter molecules (ssDNA and ssRNA) for target detection, the kinetic rates at the current assay conditions result in an amplification-free LoD in the picomolar range. The results suggest that successful detection of target requires cleavage (by an activated CRISPR enzyme) of the order of at least 0.1% of the fluorescent reporter molecules. This fraction of reporters cleaved is required to differentiate the signal from the background, and we hypothesize that this required fraction is largely independent of the detection method (e.g., endpoint vs reaction velocity) and detector sensitivity. Our results demonstrate the fundamental nature by which kinetic rates and background signal limit LoDs and thus highlight areas of improvement for the emerging field of CRISPR diagnostics.

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.

Millisecond timescale reactions observed via X-ray spectroscopy in a 3D microfabricated fused silica mixer. Corrigendum

A figure in the article by Huyke et al. [(2021), J. Synchrotron Rad. 28, 1100-1113] is corrected.

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.

Species Abundance and Reaction Off-Rate Regulate Product Formation in Reactions Accelerated Using Isotachophoresis

We present a model for second-order and pseudo-first-order reversible chemical reactions accelerated using peak-mode isotachophoresis (ITP). In such systems, ITP preconcentrates and co-locates the reactants between the leading and trailing electrolyte zones, and this significantly accelerates chemical reactions. Our model quantifies the effects of reaction rate constants and species abundance on product formation rate. We identify two key non-dimensional parameters, which are specific groupings of reaction rate constants, species concentrations, and influx rates. We then use a regular perturbation to study the effects of reverse reaction rate and relative species abundance (and relative rates of species accumulation) on production rate. We also use this perturbation method to derive an analytical expression for the quasi-steady-state production rate achievable by ITP. Our analytical models and numerical solutions are generally applicable to a wide range of systems, which use ITP to enhance reactions. The model is also an interesting case study of the complex coupling of electric field-driven species transport and reaction kinetics.

Millisecond timescale reactions observed via X-ray spectroscopy in a 3D microfabricated fused silica mixer

Determination of electronic structures during chemical reactions remains challenging in studies which involve reactions in the millisecond timescale, toxic chemicals, and/or anaerobic conditions. In this study, a three-dimensionally (3D) microfabricated microfluidic mixer platform that is compatible with time-resolved X-ray absorption and emission spectroscopy (XAS and XES, respectively) is presented. This platform, to initiate reactions and study their progression, mixes a high flow rate (0.50-1.5 ml min^-1) sheath stream with a low-flow-rate (5-90 µl min^-1) sample stream within a monolithic fused silica chip. The chip geometry enables hydrodynamic focusing of the sample stream in 3D and sample widths as small as 5 µm. The chip is also connected to a polyimide capillary downstream to enable sample stream deceleration, expansion, and X-ray detection. In this capillary, sample widths of 50 µm are demonstrated. Further, convection-diffusion-reaction models of the mixer are presented. The models are experimentally validated using confocal epifluorescence microscopy and XAS/XES measurements of a ferricyanide and ascorbic acid reaction. The models additionally enable prediction of the residence time and residence time uncertainty of reactive species as well as mixing times. Residence times (from initiation of mixing to the point of X-ray detection) during sample stream expansion as small as 2.1 ± 0.3 ms are also demonstrated. Importantly, an exploration of the mixer operational space reveals a theoretical minimum mixing time of 0.91 ms. The proposed platform is applicable to the determination of the electronic structure of conventionally inaccessible reaction intermediates.

CRISPR Enzyme Kinetics for Molecular Diagnostics

CRISPR-diagnostic assays have gained significant interest in the last few years. This interest has grown rapidly during the current COVID-19 pandemic, where CRISPR-diagnostics have been frontline contenders for rapid testing solutions. This surge in CRISPR-diagnostic research prompts the following question: what exactly are the achievable limits of detection and associated assay times enabled by the kinetics of enzymes such as Cas12 and Cas13? To explore this question, we here present a model based on Michaelis-Menten enzyme kinetics theory applied to CRISPR enzymes. We use the model to develop analytical solutions for reaction kinetics and develop back-of-the-envelope criteria to validate and check for consistency in reported enzyme kinetic parameters. We applied our analyses to all studies known to us, which report Michaelis-Menten-type kinetic data for CRISPR-associated enzymes. These studies include all subtypes of Cas12 and Cas13 and orthologs. We found all but one study clearly violate at least two of our three rules and therefore present data that violate basic physical limits. We performed an experimental study of reaction kinetics of LbCas12a with both ssDNA and dsDNA activators and use these data to validate our model and its predicted scaling. The validated model is used to explore CRISPR reaction time scales and the degree of reaction completion for practically relevant target concentrations applicable to CRISPR-diagnostic assays. The results have broad implications for achievable limits of detection and assay times of emerging, amplification-free CRISPR-detection methods.

Correction: A system for the high-throughput measurement of the shear modulus distribution of human red blood cells

Correction for 'A system for the high-throughput measurement of the shear modulus distribution of human red blood cells' by Amir Saadat et al., Lab Chip, 2020, 20, 2927-2936, DOI: 10.1039/D0LC00283F.

Effects of Weak Electrolytes on Electric Double Layer Ion Distributions

Many common experimental systems have electric double layers containing weak electrolytes, including systems with buffers. The pH at the boundary of the diffuse layer is an important parameter for determining the physicochemical state of the system, including surface charge density. We show that the Boltzmann equilibrium relation can be used as an exact solution for weak electrolyte electric double layers. Using these results, we provide a closed-form relation for the maximum pH change in a buffered electric double layer, in terms of the boundary potential. Importantly, our results suggest that equilibrium electric double layer concepts developed for strong electrolytes can be expanded to include weak electrolytes.

Simultaneous optical and infrared thermal imaging of isotachophoresis

Joule heating in isotachophoresis (ITP) can limit minimum assay times and efforts to scale up processed sample volumes. Despite its significance, the dynamics of Joule heating on spatiotemporal temperature fields in ITP systems have not been investigated. We here present novel measurements of spatiotemporal temperature and electromigration fields in ITP. To achieve this, we obtain simultaneous and registered optical and infrared thermal images of the ITP process. We conduct a series of experiments at constant current operation and vary the leading electrolyte concentration to study and highlight the importance of buffer-dependent ionic conductivity on the resulted temperature rise. The measurements demonstrate a substantial increase of temperature in the adjusted trailing electrolyte region, and the propagation of a thermal wave in the ITP channel with a velocity equal to that of the electromigration front. We present scaling of the experimental data that indicates the dependence of front velocity and temperature rise on current density and ionic conductivity. The current study has direct application to the design and optimization of scaled-up ITP systems and the validation of numerical models of Joule heating.

Process design tools and techno-economic analysis for capacitive deionization

Capacitive deionization (CDI) devices use cyclical electrosorption on porous electrode surfaces to achieve water desalination. Process modeling and design of CDI systems requires accurate treatment of the coupling among input electrical forcing, input flow rates, and system responses including salt removal dynamics, water recovery, energy storage, and dissipation. Techno-economic analyses of CDI further require a method to calculate and compare between a produced commodity (e.g. desalted water) versus capital and operational costs of the system. We here demonstrate a new modeling and analysis tool for CDI developed as an installable Matlab program that allows direct numerical simulation of CDI dynamics and calculation of key performance and cost parameters. The program is provided for free and is used to run open-source Simulink models. The Simulink environment sends information to the program and allows for a drag and drop design space where users can connect CDI cells to relevant periphery blocks such as grid energy, battery, solar panel, waste disposal, and maintenance/labor cost streams. The program allows for simulation of arbitrary current forcing and arbitrary flow rate forcing of one or more CDI cells. We employ validated well-mixed reactor formulations together with a non-linear circuit model formulation that can accommodate a variety of electric double layer sub-models (e.g. for charge efficiency). The program includes a graphical user interface (GUI) to specify CDI plant parameters, specify operating conditions, run individual tests or parameter batch-mode simulations, and plot relevant results. The techno-economic models convert among dimensional streams of species (e.g. feed, desalted water, and brine), energy, and cost and enable a variety of economic estimates including levelized water costs.

A system for the high-throughput measurement of the shear modulus distribution of human red blood cells

Reduced deformability of red blood cells (RBCs) can affect the hemodynamics of the microcirculation and reduce oxygen transport efficiency. It is also well known that reduced RBC deformability is a signature of various physical disorders, including sepsis, and that the primary determinant of RBC deformability is the membrane shear modulus. To measure the distribution of an individual's RBC shear modulus with high throughput, we a) developed a high-fidelity computational model of RBCs in confined microchannels to inform design decisions; b) created a novel experimental system combining microfluidic flow, imaging, and image analysis; and c) performed automated comparisons between measured quantities and numerical predictions to extract quantitative measures of the RBC shear modulus for each of thousands of cells. We applied our computational simulation platform to construct the appropriate deformability figure(s) of merit to quantify RBC stiffness based on an experimentally measured, steady-state cell shape in flow through a microchannel. In particular, we determined a shape parameter based on the second moment of the cell shape that is sensitive to the changes in the membrane stiffness and cell size. We then conducted microfluidic experiments and developed custom automated image processing codes to identify and track the position and shape of individual RBCs within micro-constrictions. The fabricated microchannels include a square cross-section imaging region (7 by 7 μm) and we applied order 10 kPa pressure differences to induce order 10 mm s^-1 cell velocities. The combination of modeling, microfluidics, and imaging enables, for the first time, quantitative measurement of the shear moduli of thousands of RBCs in human blood samples. We demonstrate the high-throughput features by sensitive quantification of the changes in the distribution of RBC stiffness with aging. This combined measurement and computational platform is ultimately intended to diagnose blood cell disorders in patients.

On the competition between mixing rate and uniformity in a coaxial hydrodynamic focusing mixer

Fast microfluidic mixers for use with line-of-sight integrating detection schemes pose unique challenges. Such detectors typically cannot discriminate signal from slow moving (e.g. near internal walls) and fast-moving portions of the fluid stream. This convolves reaction rate dynamics with fluid flow residence time dynamics. Further, the small cross sections of typical three-dimensional hydrodynamic focusing devices lead to lower detection signals. The current study focuses on achieving both small time scales of mixing and homogenous residence times. This is achieved by injecting sample through a center capillary and hydrodynamically focusing using a sheath flow within a tapered second capillary. The current design also features a third, larger coaxial capillary. The mixed stream flows into the large cross-section of this third capillary to decelerate and expand the stream by up to 14-fold to improve line-of-sight signal strength of reaction products. Hydrodynamic focusing, mixing, and expansion are studied using analytical and numerical models and also studied experimentally using a fluorescein-iodide quenching reaction. The experimentally validated models are used to explore trade-offs between mixing rate and uniformity. For the first time, this work presents detailed analysis of the Lagrangian time history of species transport during mixing inside coaxial capillaries to measure mixing nonuniformity. The mixing region enables order 100 μs mixing times and residence time widths of the same order (140 μs).

Using Ultramicroporous Carbon for the Selective Removal of Nitrate with Capacitive Deionization

The contamination of water resources with nitrate is a growing and significant problem. Here we report the use of ultramicroporous carbon as a capacitive deionization (CDI) electrode for selectively removing nitrate from an anion mixture. Through moderate activation, we achieve a micropore-size distribution consisting almost exclusively of narrow (<1 nm) pores that are well suited for adsorbing the planar, weakly hydrated nitrate molecule. Cyclic voltammetry measurements reveal an enhanced capacitance for nitrate when compared to chloride as well as significant ion sieving effects when sulfate is the only anion present. We measure high selectivities (S) of both nitrate over sulfate (S_NO3/SO4 = 17.8 ± 3.6 at 0.6 V) and nitrate over chloride (S_NO3/Cl = 6.1 ± 0.4 at 0.6 V) when performing a constant voltage CDI separation on 3.33 mM/3.33 mM/1.67 mM Cl/NO₃/SO₄ feedwater. These results are particularly encouraging considering that a divalent interferant was present in the feed. Using molecular dynamics simulations, we examine the solvation characteristics of these ions to better understand why nitrate is preferentially electrosorbed over sulfate and chloride.

Simultaneous RNA purification and size selection using on-chip isotachophoresis with an ionic spacer

We present an on-chip method for the extraction of RNA within a specific size range from low-abundance samples. We use isotachophoresis (ITP) with an ionic spacer and a sieving matrix to enable size-selection with a high yield of RNA in the target size range. The spacer zone separates two concentrated ITP peaks, the first containing unwanted single nucleotides and the second focusing RNA of the target size range (2-35 nt). Our ITP method excludes >90% of single nucleotides and >65% of longer RNAs (>35 nt). Compared to size selection using gel electrophoresis, ITP-based size-selection yields a 2.2-fold increase in the amount of extracted RNAs within the target size range. We also demonstrate compatibility of the ITP-based size-selection with downstream next generation sequencing. On-chip ITP-prepared samples reveal higher reproducibility of transcript-specific measurements compared to samples size-selected by gel electrophoresis. Our method offers an attractive alternative to conventional sample preparation for sequencing with shorter assay time, higher extraction efficiency and reproducibility. Potential applications of ITP-based size-selection include sequencing-based analyses of small RNAs from low-abundance samples such as rare cell types, samples from fluorescence activated cell sorting (FACS), or limited clinical samples.

High water recovery and improved thermodynamic efficiency for capacitive deionization using variable flowrate operation

Water recovery is a measure of the amount of treated water produced relative to the total amount of water processed through the system, and is an important performance metric for any desalination method. Conventional operating methods for desalination using capacitive deionization (CDI) have so far limited water recovery to be about 50%. To improve water recovery for CDI, we here introduce a new operating scheme based on a variable (in time) flow rate wherein a low flow rate during discharge is used to produce a brine volume which is significantly less than the volume of diluent produced. We demonstrate experimentally and study systematically this novel variable flowrate operating scheme in the framework of both constant current and constant voltage charge-discharge modes. We show that the variable flowrate operation can increase water recovery for CDI to very high values of ∼90% and can improve thermodynamic efficiency by about 2- to 3-fold compared to conventional constant flowrate operation. Importantly, this is achieved with minimal performance reductions in salt removal, energy consumption, and volume throughput. Our work highlights that water recovery can be readily improved for CDI at very minimal additional cost using simple flow control schemes.

Performance metrics for the objective assessment of capacitive deionization systems

In the growing field of capacitive deionization (CDI), a number of performance metrics have emerged to describe the desalination process. Unfortunately, the separation conditions under which these metrics are measured are often not specified, resulting in optimal performance at minimal removal. Here we outline a system of performance metrics and reporting conditions that resolves this issue. Our proposed system is based on volumetric energy consumption (Wh/m³) and throughput productivity (L/h/m²) reported for a specific average concentration reduction, water recovery, and feed salinity. To facilitate and rationalize comparisons between devices, materials, and operation modes, we propose a nominal standard separation of removing 5 mM from a 20 mM NaCl feed solution at 50% water recovery. We propose this particular separation as a standard, but emphasize that the rationale presented here applies irrespective of separation details. Using our proposed separation, we compare the desalination performance of a flow-through electrode (fte-CDI) cell and a flow between membrane (fb-MCDI) device, showing how significantly different systems can be compared in terms of generally desirable desalination characteristics. In general, we find that performance analysis must be considered carefully so to not allow for ambiguous separation conditions or the maximization of one metric at the expense of another. Additionally, for context and clarity, we discuss a number of important underlying performance indicators and cell characteristics that are not performance measures in and of themselves but can be examined to better understand differences in performance.

Ion selectivity in capacitive deionization with functionalized electrode: Theory and experimental validation

Capacitive deionization (CDI) is a promising technique for salt removal and may have potential for highly selective removal of ion species. In this work, we take advantage of functional groups usually used with ionic exchange resins and apply these to CDI. To this end, we functionalize activated carbon with a quaternary amines surfactant and use this surface to selectively and passively (no applied field) trap nitrate ions. We then set the cell voltage to a constant value to regenerate these electrodes, resulting in an inverted capacitive deionization (i-CDI) operation. Unlike resins, we avoid use of concentrated chemicals for regeneration. We measure the selectivity of nitrate versus chloride ions as a function of regeneration voltage and initial chloride concentration. We experimentally demonstrate up to about 6.5-fold (observable) selectivity in a cycle with a regeneration voltage of 0.4 V. We also demonstrate a novel multi-pass, air-flush i-CDI operation to selectively enrich nitrate with high water recovery. We further present a dynamic, multi-species electrosorption and equilibrium solution-to-surface chemical reaction model and validate the model with detailed measurements. Our i-CDI system exhibits higher nitrate selectivity at lower voltages; making it possible to reduce NaNO3 concentrations from ∼170 ppm to below the limit of maximum allowed values for nitrate in drinking water of about 50 ppm NaNO3.

A method for quantifying in plane permeability of porous thin films

The in-plane permeability of porous thin films is an important fluid mechanical property that determines wicking and pressure-driven flow behavior in such materials. This property has so far been challenging to measure directly due to the small sidewall cross-sectional area of thin films available for flow. In this work, we propose and experimentally demonstrate a novel technique for directly measuring in-plane permeability of porous thin films of arbitrary thicknesses, in situ, using a manifold pressed to the top surface of the film. We both measure and simulate the influence of the two dimensional flow field produced in a film by the manifold and extract the permeability from measurements of pressure drop at fixed flow rates. Permeability values measured using the technique for a periodic array of channels are comparable to theoretical predictions. We also determine in-plane permeability of arrays of pillars and electrodeposited porous copper films. This technique is a robust tool to characterize permeability of thin films of arbitrary thicknesses on a variety of substrates. In Supplementary material, we provide a solid model, which is useful in three-dimensional printer reproductions of our device.

Frequency analysis and resonant operation for efficient capacitive deionization

Capacitive deionization (CDI) performance metrics can vary widely with operating methods. Conventional CDI operating methods such as constant current and constant voltage show advantages in either energy or salt removal performance, but not both. We here develop a theory around and experimentally demonstrate a new operation for CDI that uses sinusoidal forcing voltage (or sinusoidal current). We use a dynamic system modeling approach, and quantify the frequency response (amplitude and phase) of CDI effluent concentration. Using a wide range of operating conditions, we demonstrate that CDI can be modeled as a linear time invariant system. We validate this model with experiments, and show that a sinusoid voltage operation can simultaneously achieve high salt removal and strong energy performance, thus very likely making it superior to other conventional operating methods. Based on the underlying coupled phenomena of electrical charge (and ionic) transfer with bulk advection in CDI, we derive and validate experimentally the concept of using sinusoidal voltage forcing functions to achieve resonance-type operation for CDI. Despite the complexities of the system, we find a simple relation for the resonant time scale: the resonant time period (frequency) is proportional (inversely proportional) to the geometric mean of the flow residence time and the electrical (RC) charging time. Operation at resonance implies the optimal balance between absolute amount of salt removed (in moles) and dilution (depending on the feed volume processed), thus resulting in the maximum average concentration reduction for the desalinated water. We further develop our model to generalize the resonant time-scale operation, and provide responses for square and triangular voltage waveforms as two examples. To this end, we develop a general tool that uses Fourier analysis to construct CDI effluent dynamics for arbitrary input waveforms. Using this tool, we show that most of the salt removal (∼95%) for square and triangular voltage forcing waveforms is achieved by the fundamental Fourier (sinusoidal) mode. The frequency of higher Fourier modes precludes high flow efficiency for these modes, so these modes consume additional energy for minimal additional salt removed. This deficiency of higher frequency modes further highlights the advantage of DC-offset sinusoidal forcing for CDI operation.

Efficient Production of On-Target Reads for Small RNA Sequencing of Single Cells Using Modified Adapters

Although single-cell mRNA sequencing has been a powerful tool to explore cellular heterogeneity, the sequencing of small RNA at the single-cell level (sc-sRNA-seq) remains a challenge, as these have no consensus sequence, are relatively low abundant, and are difficult to amplify in a bias-free fashion. We present two methods of single-cell-lysis that enable sc-sRNA-seq. The first method is a chemical-based technique with overnight freezing while the second method leverages on-chip electrical lysis of plasma membrane and physical extraction and separation of cytoplasmic RNA via isotachophoresis. We coupled these two methods with off-chip small RNA library preparation using CleanTag modified adapters to prevent the formation of adapter dimers. We then demonstrated sc-sRNA-seq with single K562 human leukemic cells. Our approaches offer a relatively short hands-on time of 6 h and efficient generation of on-target reads. The sc-sRNA-seq with our approaches showed detection of miRNA with various abundances ranging from 16 000 copies/cell to about 10 copies/cell. We anticipate this approach will create a new opportunity to explore cellular heterogeneity through small RNA expression.

Tailoring Permeability of Microporous Copper Structures through Template Sintering

Microporous metals are used extensively for applications that combine convective and conductive transport and benefit from low resistance to both modes of transport. Conventional fabrication methods, such as direct sintering of metallic particles, however, often produce structures with limited fluid transport properties due to the lack of control over pore morphologies such as the pore size and porosity. Here, we demonstrate control and improvement of hydraulic permeability of microporous copper structures fabricated using template-assisted electrodeposition. Template sintering is shown to modify the fluid transport network in a manner that increases permeability by nearly an order of magnitude with a less significant decrease (∼38%) in thermal conductivity. The measured permeabilities range from 4.8 × 10^-14 to 1.3 × 10^-12 m² with 5 μm pores, with the peak value being roughly 5 times larger than the published values for sintered copper particles of comparable feature sizes. Analysis indicates that the enhancement of permeability is limited by constrictions, i.e., bottlenecks between connecting pores, whose dimensions are highly sensitive to the sintering conditions. We further show contrasting trends in permeability versus conductivity of the electrodeposited microporous copper and conventional sintered copper particles and suggest these differing trends to be the result of their inverse structural relationship.

Thermodynamics of Ion Separation by Electrosorption

We present a simple, top-down approach for the calculation of minimum energy consumption of electrosorptive ion separation using variational form of the (Gibbs) free energy. We focus and expand on the case of electrostatic capacitive deionization (CDI). The theoretical framework is independent of details of the double-layer charge distribution and is applicable to any thermodynamically consistent model, such as the Gouy-Chapman-Stern and modified Donnan models. We demonstrate that, under certain assumptions, the minimum required electric work energy is indeed equivalent to the free energy of separation. Using the theory, we define the thermodynamic efficiency of CDI. We show that the thermodynamic efficiency of current experimental CDI systems is currently very low, around 1% for most existing systems. We applied this knowledge and constructed and operated a CDI cell to show that judicious selection of the materials, geometry, and process parameters can lead to a 9% thermodynamic efficiency and 4.6 kT per removed ion energy cost. This relatively high thermodynamic efficiency is, to our knowledge, by far the highest thermodynamic efficiency ever demonstrated for traditional CDI. We hypothesize that efficiency can be further improved by further reduction of CDI cell series resistances and optimization of operational parameters.

Self similarities in desalination dynamics and performance using capacitive deionization

Charge transfer and mass transport are two underlying mechanisms which are coupled in desalination dynamics using capacitive deionization (CDI). We developed simple reduced-order models based on a mixed reactor volume principle which capture the coupled dynamics of CDI operation using closed-form semi-analytical and analytical solutions. We use the models to identify and explore self-similarities in the dynamics among flow rate, current, and voltage for CDI cell operation including both charging and discharging cycles. The similarity approach identifies the specific combination of cell (e.g. capacitance, resistance) and operational parameters (e.g. flow rate, current) which determine a unique effluent dynamic response. We here demonstrate self-similarity using a conventional flow between CDI (fbCDI) architecture, and we hypothesize that our similarity approach has potential application to a wide range of designs. We performed an experimental study of these dynamics and used well-controlled experiments of CDI cell operation to validate and explore limits of the model. For experiments, we used a CDI cell with five electrode pairs and a standard flow between (electrodes) architecture. Guided by the model, we performed a series of experiments that demonstrate natural response of the CDI system. We also identify cell parameters and operation conditions which lead to self-similar dynamics under a constant current forcing function and perform a series of experiments by varying flowrate, currents, and voltage thresholds to demonstrate self-similarity. Based on this study, we hypothesize that the average differential electric double layer (EDL) efficiency (a measure of ion adsorption rate to EDL charging rate) is mainly dependent on user-defined voltage thresholds, whereas flow efficiency (measure of how well desalinated water is recovered from inside the cell) depends on cell volumes flowed during charging, which is determined by flowrate, current and voltage thresholds. Results of experiments strongly support this hypothesis. Results show that cycle efficiency and salt removal for a given flowrate and current are maximum when average EDL and flow efficiencies are approximately equal. We further explored a range of CC operations with varying flowrates, currents, and voltage thresholds using our similarity variables to highlight trade-offs among salt removal, energy, and throughput performance.

Self-Cleaning Porous Surfaces for Dry Condensation

Efficient water removal from a cool surface during condensation is critically important to the enhancement of a variety of heat transfer applications. Previous work has focused on the fabrication of superhydrophobic surfaces which promote water droplets and removal via droplet shedding or jumping. Here, we report a novel strategy with a droplet self-cleaning surface which spontaneously transports all of the droplets from the condensation surface to the back side. We fabricate the self-cleaning surface by simply tailoring the wettability of the two sides of a porous membrane and demonstrate that the hydrophobic side is effective in clearing off droplets of a wide range of diameters. Even during rapid impingement of droplets smaller than 10 μm, this surface remains dry. We further demonstrate a "dry condensation" process wherein a surface undergoing rapid condensation is maintained free of droplets. This minimizes the essential thermal resistance of the process, and we estimate a twofold increase in condensation rate compared with a simple copper surface under the same conditions. Our method is tailorable, extendable to a wide range of materials and geometries, and shows great potential for a broad range of condensation processes.

SINC-seq: correlation of transient gene expressions between nucleus and cytoplasm reflects single-cell physiology

We report a microfluidic system that physically separates nuclear RNA (nucRNA) and cytoplasmic RNA (cytRNA) from a single cell and enables single-cell integrated nucRNA and cytRNA-sequencing (SINC-seq). SINC-seq constructs two individual RNA-seq libraries, nucRNA and cytRNA, per cell, quantifies gene expression in the subcellular compartments, and combines them to create novel single-cell RNA-seq data. Leveraging SINC-seq, we discover distinct natures of correlation among cytRNA and nucRNA that reflect the transient physiological state of single cells. These data provide unique insights into the regulatory network of messenger RNA from the nucleus toward the cytoplasm at the single-cell level.

Quantifying the flow efficiency in constant-current capacitive deionization

Here we detail a previously unappreciated loss mechanism inherent to capacitive deionization (CDI) cycling operation that has a substantial role determining performance. This mechanism reflects the fact that desalinated water inside a cell is partially lost to re-salination if desorption is carried out immediately after adsorption. We describe such effects by a parameter called the flow efficiency, and show that this efficiency is distinct from and yet multiplicative with other highly-studied adsorption efficiencies. Flow losses can be minimized by flowing more feed solution through the cell during desalination; however, this also results in less effluent concentration reduction. While the rationale outlined here is applicable to all CDI cell architectures that rely on cycling, we validate our model with a flow-through electrode CDI device operated in constant-current mode. We find excellent agreement between flow efficiency model predictions and experimental results, thus giving researchers simple equations by which they can estimate this distinct loss process for their operation.

Assay for Listeria monocytogenes cells in whole blood using isotachophoresis and recombinase polymerase amplification

We present a new approach which enables lysis, extraction, and detection of inactivated Listeria monocytogenes cells from blood using isotachophoresis (ITP) and recombinase polymerase amplification (RPA). We use an ITP-compatible alkaline and proteinase K approach for rapid and effective lysis. We then perform ITP purification to separate bacterial DNA from whole blood contaminants using a microfluidic device that processes 25 μL sample volume. Lysis, mixing, dispensing, and on-chip ITP purification are completed in a total of less than 50 min. We transfer extracted DNA directly into RPA master mix for isothermal incubation and detection, an additional 25 min. We first validate our assay in the detection of purified genomic DNA spiked into whole blood, and demonstrate a limit of detection of 16.7 fg μL^-1 genomic DNA, the equivalent of 5 × 10³ cells per mL. We then show detection of chemically-inactivated L. monocytogenes cells spiked into whole blood, and demonstrate a limit of detection of 2 × 10⁴ cells per mL. Lastly, we show preliminary experimental data demonstrating the feasibility of the integration of ITP purification with RPA detection on a microfluidic chip. Our results suggest that ITP purification is compatible with RPA detection, and has potential to extend the applicability of RPA to whole blood.

Isotachophoresis applied to biomolecular reactions

This review discusses research developments and applications of isotachophoresis (ITP) to the initiation, control, and acceleration of chemical reactions, emphasizing reactions involving biomolecular reactants such as nucleic acids, proteins, and live cells. ITP is a versatile technique which requires no specific geometric design or material, and is compatible with a wide range of microfluidic and automated platforms. Though ITP has traditionally been used as a purification and separation technique, recent years have seen its emergence as a method to automate and speed up chemical reactions. ITP has been used to demonstrate up to 14 000-fold acceleration of nucleic acid assays, and has been used to enhance lateral flow and other immunoassays, and even whole bacterial cell detection assays. We here classify these studies into two categories: homogeneous (all reactants in solution) and heterogeneous (at least one reactant immobilized on a solid surface) assay configurations. For each category, we review and describe physical modeling and scaling of ITP-aided reaction assays, and elucidate key principles in ITP assay design. We summarize experimental advances, and identify common threads and approaches which researchers have used to optimize assay performance. Lastly, we propose unaddressed challenges and opportunities that could further improve these applications of ITP.

Equilibria model for pH variations and ion adsorption in capacitive deionization electrodes

Ion adsorption and equilibrium between electrolyte and microstructure of porous electrodes are at the heart of capacitive deionization (CDI) research. Surface functional groups are among the factors which fundamentally affect adsorption characteristics of the material and hence CDI system performance in general. Current CDI-based models for surface charge are mainly based on a fixed (constant) charge density, and do not treat acid-base equilibria of electrode microstructure including so-called micropores. We here expand current models by coupling the modified Donnan (mD) model with weak electrolyte acid-base equilibria theory. In our model, surface charge density can vary based on equilibrium constants (pK's) of individual surface groups as well as micropore and electrolyte pH environments. In this initial paper, we consider this equilibrium in the absence of Faradaic reactions. The model shows the preferential adsorption of cations versus anions to surfaces with respectively acidic or basic surface functional groups. We introduce a new parameter we term "chemical charge efficiency" to quantify efficiency of salt removal due to surface functional groups. We validate our model using well controlled titration experiments for an activated carbon cloth (ACC), and quantify initial and final pH of solution after adding the ACC sample. We also leverage inductively coupled plasma mass spectrometry (ICP-MS) and ion chromatography (IC) to quantify the final background concentrations of individual ionic species. Our results show a very good agreement between experiments and model. The model is extendable to a wide variety of porous electrode systems and CDI systems with applied potential.

Rapid Hydrogen-Deuterium Exchange in Liquid Droplets

The rate of hydrogen-deuterium exchange (HDX) in aqueous droplets of phenethylamine has been determined with submillisecond temporal resolution by mass spectrometry using nanoelectrospray ionization with a theta-capillary. The average speed of the microdroplets is measured using microparticle image velocimetry. The droplet travel time is varied from 20 to 320 μs by changing the distance between the emitter and the heated inlet to the mass spectrometer and the voltage applied to the emitter source. The droplets were found to accelerate by ∼30% during their observable travel time. Our droplet imaging shows that the theta-capillary produces two Taylor cone-jets (one per channel), causing mixing to take place from droplet fusion in the Taylor spray zone. Phenethylamine (ϕCH₂CH₂NH₂) was chosen to study because it has only one functional group (-NH₂) that undergoes rapid HDX. We model the HDX with a system of ordinary differential equations. The rate constant for the formation of -NH₂D⁺ from -NH₃⁺ is 3660 ± 290 s^-1, and the rate constant for the formation of -NHD₂⁺ from -NH₂D⁺ is 3330 ± 270 s^-1. The observed rates are about 3 times faster than what has been reported for rapidly exchangeable peptide side-chain groups in bulk measurements using stopped-flow kinetics and NMR spectroscopy. We also applied this technique to determine the HDX rates for a small 10-residue peptide, angiotensin I, in aqueous droplets, from which we found a 7-fold acceleration of HDX in the droplet compared to that in bulk solution.

Energy breakdown in capacitive deionization

We explored the energy loss mechanisms in capacitive deionization (CDI). We hypothesize that resistive and parasitic losses are two main sources of energy losses. We measured contribution from each loss mechanism in water desalination with constant current (CC) charge/discharge cycling. Resistive energy loss is expected to dominate in high current charging cases, as it increases approximately linearly with current for fixed charge transfer (resistive power loss scales as square of current and charging time scales as inverse of current). On the other hand, parasitic loss is dominant in low current cases, as the electrodes spend more time at higher voltages. We built a CDI cell with five electrode pairs and standard flow between architecture. We performed a series of experiments with various cycling currents and cut-off voltages (voltage at which current is reversed) and studied these energy losses. To this end, we measured series resistance of the cell (contact resistances, resistance of wires, and resistance of solution in spacers) during charging and discharging from voltage response of a small amplitude AC current signal added to the underlying cycling current. We performed a separate set of experiments to quantify parasitic (or leakage) current of the cell versus cell voltage. We then used these data to estimate parasitic losses under the assumption that leakage current is primarily voltage (and not current) dependent. Our results confirmed that resistive and parasitic losses respectively dominate in the limit of high and low currents. We also measured salt adsorption and report energy-normalized adsorbed salt (ENAS, energy loss per ion removed) and average salt adsorption rate (ASAR). We show a clear tradeoff between ASAR and ENAS and show that balancing these losses leads to optimal energy efficiency.

Rapid Slow Off-Rate Modified Aptamer (SOMAmer)-Based Detection of C-Reactive Protein Using Isotachophoresis and an Ionic Spacer

We present an on-chip electrophoretic assay for rapid protein detection with a SOMAmer (Slow Off-Rate Modified Aptamer) reagent. We used isotachophoresis (ITP) coupled with an ionic spacer to both react and separate SOMAmer-protein complex from free SOMAmer reagent. ITP accelerates the reaction kinetics as the ionic spacer concurrently separates the reaction products. We developed a numerical and analytical model to describe ITP spacer assays, which involve low-mobility, nonfocusing targets that are recruited into the ITP zone by higher-mobility, ITP-focused probes. We demonstrated a proof-of-concept of this assay using C-reactive protein (CRP) in buffer, and achieved a 2 nM limit of detection (LOD) with a combined 20 min assay time (10 min off-chip preparation of reagents and 10 min on-chip run). Our findings suggest that this approach has potential as a simple and rapid alternative to other homogeneous immunoassays. We also explore the extension of this assay to a diluted serum sample spiked with CRP, where we observe decreased sensitivity (an LOD of 25 nM in 20-fold diluted serum). We describe the challenges in extending this assay to complex samples and achieving higher sensitivity and specificity for clinical applications.

Isotachophoresis for fractionation and recovery of cytoplasmic RNA and nucleus from single cells

There is a substantial need for simultaneous analyses of RNA and DNA from individual single cells. Such analysis provides unique evidence of cell-to-cell differences and the correlation between gene expression and genomic mutation in highly heterogeneous cell populations. We present a novel microfluidic system that leverages isotachophoresis to fractionate and isolate cytoplasmic RNA and genomic DNA (gDNA) from single cells. The system uniquely enables independent, sequence-specific analyses of these critical markers. Our system uses a microfluidic chip with a simple geometry and four end-channel electrodes, and completes the entire process in <5 min, including lysis, purification, fractionation, and delivery to DNA and RNA output reservoirs, each containing high quality and purity aliquots with no measurable cross-contamination of cytoplasmic RNA versus gDNA. We demonstrate our system with simultaneous, sequence-specific quantitation using off-chip RT-qPCR and qPCR for simultaneous cytoplasmic RNA and gDNA analyses, respectively.

Transient delivery of modified mRNA encoding TERT rapidly extends telomeres in human cells

Telomere extension has been proposed as a means to improve cell culture and tissue engineering and to treat disease. However, telomere extension by nonviral, nonintegrating methods remains inefficient. Here we report that delivery of modified mRNA encoding TERT to human fibroblasts and myoblasts increases telomerase activity transiently (24-48 h) and rapidly extends telomeres, after which telomeres resume shortening. Three successive transfections over a 4 d period extended telomeres up to 0.9 kb in a cell type-specific manner in fibroblasts and myoblasts and conferred an additional 28 ± 1.5 and 3.4 ± 0.4 population doublings (PDs), respectively. Proliferative capacity increased in a dose-dependent manner. The second and third transfections had less effect on proliferative capacity than the first, revealing a refractory period. However, the refractory period was transient as a later fourth transfection increased fibroblast proliferative capacity by an additional 15.2 ± 1.1 PDs, similar to the first transfection. Overall, these treatments led to an increase in absolute cell number of more than 10(12)-fold. Notably, unlike immortalized cells, all treated cell populations eventually stopped increasing in number and expressed senescence markers to the same extent as untreated cells. This rapid method of extending telomeres and increasing cell proliferative capacity without risk of insertional mutagenesis should have broad utility in disease modeling, drug screening, and regenerative medicine.