Modeling and simulation of chemo-electro-mechanical behavior of pH-electric-sensitive hydrogel

A chemo-electro-mechanical multi-field model, termed the multi-effect-coupling pH-electric-stimuli (MECpHe) model, has been developed to simulate the response behavior of smart hydrogels subject to pH and electric voltage coupled stimuli when the hydrogels are immersed in a pH buffer solution subject to an externally applied electric field. The MECpHe model developed considers multiphysics effects and formulates the fixed charge density with the coupled buffer solution pH and electric voltage effects, expressed by a set of nonlinear partial differential governing equations. The model can be used to predict the hydrogel displacement and the distributive profiles of the concentrations of diffusive ionic species and the electric potential and the fixed charge density in both the hydrogels and surrounding solution. After validation of the model by comparison of current numerical results with experiment data extracted from the literature, one-dimensional steady-state simulations were carried out for equilibrium of the smart hydrogels subject to pH and electric coupled stimuli. The effects of several important physical conditions, including the externally applied electric voltage, on the distributions of the concentrations of diffusive ionic species, the electric potential, the fixed charge density, and the displacement of the hydrogel strip were studied in detail. The effects of the ionic strength on the bending deformation of the hydrogels under the solution pH and electric voltage coupled stimuli are also discussed.     

Electrokinetic transport through the nanopores in cell membrane during electroporation

In electroporation, applied electric field creates hydrophilic nanopores in a cell membrane that can serve as a pathway for inserting biological samples to the cell. It is highly desirable to understand the ionic transfer and fluid flow through the nanopores in order to control and improve the cell transfection. Because of submicron dimensions, conventional theories of electrokinetics may lose their applicability in such nanopores. In the current study, the Poisson-Nernst-Planck equations along with modified Navier-Stokes equations and the continuity equation are solved in order to find electric potential, fluid flow, and ionic concentration through the nanopores. The results show that the electric potential, velocity field, and ionic concentration vary with the size of the nanopores and are different through the nanopores located at the front and backside of the cell membrane. However, on a given side of the cell membrane, angular position of nanopores has fewer influences on liquid flow and ionic transfer. By increasing the radius of the nanopores, the averaged velocity and ionic concentration through the nanopores are increased. It is also shown that, in the presence of electric pulse, electrokinetic effects (electroosmosis and electrophoresis) have significant influences on ionic mass transfer through the nanopores, while the effect of diffusion on ionic mass flux is negligible in comparison with electrokinetics. Increasing the radius of the nanopores intensifies the effect of convection (electroosmosis) in comparison with electrophoresis on ionic flux.     

Electrochemical response and separation in cyclic electric field-flow fractionation

Electric field-flow fractionation (EFFF) is a separation technique that couples a lateral electric field with axial Poiseuille flow to separate particles on the basis of size and/or mobility. In unidirectional EFFF, the field rapidly decreases in time due to charging of the double layer. The field strength could be increased by performing EFFF with cyclic electric fields. In cyclic electric field-flow fractionation (CEFFF), a periodic voltage, which can be either sinusoidal or square-wave, is applied in the lateral direction. In this paper, we measure the electrochemical response of CEFFF, i.e., the current-time response for a given time-dependent voltage and then utilize this electrochemical response in a transport model to predict separation. The CEFFF device studied here comprises two gold-coated glass plates separated by a spacer. The transient current profiles are measured for a step change and cyclic square-shaped voltage. The current profile is compared with the equivalent circuit model, and is fitted to a sum of two decaying exponentials. The dependence of the electrochemical response on voltage, frequency, channel thickness, and salt concentration is studied. Next, the electrochemical data are utilized in the convection-diffusion equation to develop a model for separation by CEFFF. The equations are solved by using a combination of analytical and numerical techniques to determine the mean velocity and the dispersion coefficient of molecules, and to determine the effect of various parameters on the separation efficiency of the EFFF device. Also, the model predictions are compared with experimental data available in the literature.     

Electrokinetic particle translocation through a nanopore containing a floating electrode

Electrokinetic particle translocation through a nanopore containing a floating electrode is investigated by solving a continuum model, composed of the coupled Poisson-Nernst-Planck (PNP) equations for the ionic mass transport and the modified Stokes equations for the flow field. Two effects due to the presence of the floating electrode, the induced-charge electroosmosis (ICEO) and the particle-floating electrode electrostatic interaction, could significantly affect the electrokinetic mobility of DNA nanoparticles. When the electrical double layers (EDLs) of the DNA nanoparticle and the floating electrode are not overlapped, the particle-floating electrode electrostatic interaction becomes negligible. As a result, the DNA nanoparticle could be trapped near the floating electrode arising from the induced-charge electroosmosis when the applied electric field is relatively high. The presence of the floating electrode attracts more ions inside the nanopore resulting in an increase in the ionic current flowing through the nanopore; however, it has a limited effect on the deviation of the current from its base current when the particle is far from the pore.     

Electrokinetic particle translocation through a nanopore

Nanoparticle electrophoretic translocation through a single nanopore induces a detectable change in the ionic current, which enables the nanopore-based sensing for various bio-analytical applications. In this study, a transient continuum-based model is developed for the first time to investigate the electrokinetic particle translocation through a nanopore by solving the Nernst-Planck equations for the ionic concentrations, the Poisson equation for the electric potential and the Navier-Stokes equations for the flow field using an arbitrary Lagrangian-Eulerian (ALE) method. When the applied electric field is relatively low, a current blockade is expected. In addition, the particle could be trapped at the entrance of the nanopore when the electrical double layer (EDL) adjacent to the charged particle is relatively thick. When the electric field imposed is relatively high, the particle can always pass through the nanopore by electrophoresis. However, a current enhancement is predicted if the EDL of the particle is relatively thick. The obtained numerical results qualitatively agree with the existing experimental results. It is also found that the initial orientation of the particle could significantly affect the particle translocation and the ionic current through a nanopore. Furthermore, a relatively high electric field tends to align the particle with its longest axis parallel to the local electric field. However, the particle's initial lateral offset from the centerline of the nanopore acts as a minor effect.     

Electrokinetics in nanochannels: part I. Electric double layer overlap and channel-to-well equilibrium

In this paper a new model is described for calculating the electric potential field in a long, thin nanochannel with overlapped electric double layers. Electrolyte concentration in the nanochannel is predicted self-consistently via equilibrium between ionic solution in the wells and within the nanochannel. Differently than published models that require detailed iterative numerical solutions of coupled differential equations, the framework presented here is self-consistent and predictions are obtained solving a simple one-dimensional integral. The derivation clearly shows that the electric potential field depends on three new parameters: the ratio of ion density in the channel to ion density in the wells; the ratio of free-charge density to bulk ion density within the channel; and a modified Debye-Hückel thickness, which is the relevant scale for shielding of surface net charge. For completeness, three wall-surface boundary conditions are analyzed: specified zeta-potential; specified surface net charge density; and charge regulation. Predictions of experimentally observable quantities based on the model proposed here, such as depth-averaged electroosmotic flow and net ionic current, are significantly different than results from previous overlapped electric double layer models. In this first paper of a series of two, predictions are presented where channel depth is varied at constant well concentration. Results show that under conditions of electric double layer overlap, electroosmosis contributes only a small fraction of the net ionic current, and that most of the measurable current is due to ionic conduction in conditions of increased counterion density in the nanochannel. In the second of this two-paper series, predictions are presented where well-concentration is varied and the channel depth is held constant, and the model described here is employed to study the dependence of ion mobility on ionic strength, and compare predictions to measurements of ionic current as a function of channel depth and ion density.     

A coarse‐grained MARTINI‐like force field for DNA unzipping in nanopores

In nanopore force spectroscopy (NFS) a charged polymer is threaded through a channel of molecular dimensions. When an electric field is applied across the insulating membrane, the ionic current through the nanopore reports on polymer translocation, unzipping, dissociation, and so forth. We present a new model that can be applied in molecular dynamics simulations of NFS. Although simplified, it does reproduce experimental trends and all‐atom simulations. The scaled conductivities in bulk solution are consistent with experimental results for NaCl for a wide range of electrolyte concentrations and temperatures. The dependence of the ionic current through a nanopore on the applied voltage is symmetric and, in the voltage range used in experiments (up to 2 V), linear and in good agreement with experimental data. The thermal stability and geometry of DNA is well represented. The model was applied to simulations of DNA hairpin unzipping in nanopores. The results are in good agreement with all‐atom simulations: the scaled translocation times and unzipping sequence are similar. © 2015 Wiley Periodicals, Inc.     

Direct numerical simulation of electrokinetic translocation of a cylindrical particle through a nanopore using a Poisson-Boltzmann approach

Nanopore-based sensing of single molecules is based on a detectable change in the ionic current arising from the electrokinetic translocation of individual nanoparticles through a nanopore. In this study, we propose a continuum-based model to investigate the dynamic electrokinetic translocation of a cylindrical nanoparticle through a nanopore and the corresponding ionic current response. It is the first time to simultaneously solve the Poisson-Boltzmann equation for the ionic concentrations and the electric field contributed by the surface charges of the nanopore and the nanoparticle, the Laplace equation for the externally applied electric field, and the modified Stokes equations for the flow field using an arbitrary Lagrangian-Eulerian method. Current blockade due to the particle translocation is predicted when the electric double layers (EDLs) of the particle and the nanopore are not overlapped, which is in qualitative agreement with existing experimental observations. Effects due to the electric field intensity imposed, the EDL thickness, the nanopore's surface charge, the particle's initial orientation and lateral offset from the nanopore's centerline on the particle translocation including both translation and rotation, and the ionic current response are comprehensively investigated. Under a relatively low electric field imposed, the particle experiences a significant rotation and a lateral movement. However, the particle is aligned with its longest axis parallel to the local electric field very quickly due to the dielectrophoretic effect when the external electric field is relatively high.     

Flow of an Anisotropic Dielectric Liquid in a Diverging Channel in a Strong Transverse Electric Field

A steady plane flow of an anisotropically polarizable liquid in a channel with nonparallel walls was considered. One of the walls was grounded, and the other was under a high electric potential. The polarization anisotropy was described in terms of a unit vector whose direction was determined by a relaxation equation. The dependence of the polarization of the liquid on the strength of the electric field and the anisotropy vector was specified using an equilibrium relation. Such a model can describe, for example, a suspension of anisotropically polarizable particles in a highly insulating liquid. The velocity, pressure, polarization, anisotropy vector, and electric field distributions in the liquid were determined and investigated. It was shown that, at some critical Reynolds number, backflows are initiated near the channel walls. The dependence of the critical Reynolds number on the diverging angle of the channel and on the properties of a liquid in a strong electric field was determined. The applied electric field increases the critical Reynolds number, which provides a means of controlling the regime of the considered flow using electrical methods.     

Effective screening length of isotropic liquid samples submitted to an applied voltage

A cell of isotropic liquid in the shape of a slab of thickness d and containing ionic impurities is considered. It is shown that the screening effect produced by the ionic charges on the external field is characterized by an effective surface length, lambda(S)(U), depending on the applied voltage U. The analysis indicates that lambda(S)(U)) < lambda(D) when the applied voltage is very large, and lambda(S)(U) --> lambda(D) for very small values of the applied voltage, where lambda(D) is the Debye screening length. The presence of the ions is responsible also for a counterpotential, v, that for small U is such to cancel the effective electric field in the sample, whereas in the opposite limit it is inversely proportional to the applied difference of potential.     

Numerical modeling of Joule heating-induced temperature gradient focusing in microfluidic channels

Temperature gradient focusing (TGF) is a recently developed technique for spatially focusing and separating ionic analytes in microchannels. The temperature gradient required for TGF can be generated either by an imposed temperature gradient or by Joule heating resulting from an applied electric field that also drives the flow. In this study, a comprehensive numerical model describing the Joule heating induced temperature development and TGF is developed. The model consists of a set of governing equations including the Poisson-Boltzmann equation, the Laplace equation, the Navier-Stokes equations, the energy equations and the mass transport equation. As the thermophysical and electrical properties including the liquid dielectric constant, viscosity, and electric conductivity are temperature-dependent, these governing equations are coupled, and therefore the coupled governing equations are solved numerically by using a CFD-based numerical method. The numerical simulations agree well with the experimental results, suggesting the valid mathematical model presented in this study.     

The Electrophoretic Mobility and Electric Conductivity of a Concentrated Suspension of Colloidal Spheres with Arbitrary Double-Layer Thickness

The electrophoresis in a monodisperse suspension of dielectric spheres with an arbitrary thickness of the electric double layers is analytically studied. The effects of particle interactions are taken into account by employing a unit cell model, and the overlap of the double layers of adjacent particles is allowed. The electrokinetic equations, which govern the ionic concentration distributions, the electric potential profile, and the fluid flow field in the electrolyte solution surrounding the charged sphere in a unit cell, are linearized assuming that the system is only slightly distorted from equilibrium. Using a perturbation method, these linearized equations are solved with the surface charge density (or zeta potential) of the particle as the small perturbation parameter. Analytical expressions for the electrophoretic mobility of the colloidal sphere in closed form correct to O(zeta) are obtained. Based on the solution of the electrokinetic equations in a cell, a closed-form formula for the electric conductivity of the suspension up to O(zeta(2)) is derived from the average electric current density. Comparisons of the results of the cell model with different conditions at the outer boundary of the cell are made for both the electrophoretic mobility and the electric conductivity. Copyright 2001 Academic Press.     

Diffusioosmotic flows in slit nanochannels

Diffusioosmotic flows of electrolyte solutions in slit nanochannels with homogeneous surface charges induced by electrolyte concentration gradients in the absence of externally applied pressure gradients and potential differences are investigated theoretically. A continuum mathematical model consisting of the strongly coupled Nernst-Planck equations for the ionic species' concentrations, the Poisson equation for the electric potential in the electrolyte solution, and the Navier-Stokes equations for the flow field is numerically solved simultaneously. The induced diffusioosmotic flow through the nanochannel is computed as functions of the externally imposed concentration gradient, the concentration of the electrolyte solution, and the surface charge density along the walls of the nanochannel. With the externally applied electrolyte concentration gradient, a strongly spatially dependent electric field and pressure gradient are induced within the nanochannel that, in turn, generate a spatially dependent diffusioosmotic flow. The diffusioosmotic flow is opposite to the applied concentration gradient for a relatively low bulk electrolyte concentration. However, the electrolyte solution flows from one end of the nanochannel with a higher electrolyte concentration to the other end with a lower electrolyte concentration when the bulk electrolyte concentration is relatively high. There is an optimal concentration gradient under which the flow rate attains the maximum. The induced flow is enhanced with the increase in the fixed surface charge along the wall of the nanochannel for a relatively low bulk electrolyte concentration.     

Threading synthetic polyelectrolytes through protein pores

We have measured the ionic current signatures of sodium poly(styrene sulfonate) as its single molecules translocate through an alpha-hemolysin pore embedded into a bilayer in a salty aqueous medium under an externally applied electric field. As in the previous experiments involving DNA and RNA, the pore current, which is a measure of the ionic conductivity of the low molar mass electrolyte ions, is significantly reduced when the polymer molecule translocates through the pore. The magnitude and the duration of the reduction in the pore current are measured for each of the translocation events. By studying thousands of events of reduction in the ionic current, we have constructed distribution functions for the extent of the reduced current and for the translocation time. The details of these distribution functions are significantly different from those for DNA and RNA. By investigating over two orders of magnitude in the molecular weight of the polymer, the average translocation time is found to be proportional to the molecular weight and inversely proportional to the applied voltage. This demonstration of threading a synthetic polyelectrolyte through a protein pore opens up many opportunities to systematically explore the fundamental physical principles behind translocation of single macromolecules, by resorting to the wide variety of synthetically available polymers without the complexities arising from the sequences of biological polymers. In addition, the present experiments suggest yet another experimental protocol for separation of polymer molecules directly in aqueous media.     

Model‐based analysis of a dielectrophoretic microfluidic device for field‐flow fractionation

We present the development of a dynamic model for predicting the trajectory of microparticles in microfluidic devices, employing dielectrophoresis, for Hyperlayer field‐flow fractionation. The electrode configuration is such that multiple finite‐sized electrodes are located on the top and bottom walls of the microchannel; the electrodes on the walls are aligned with each other. The electric potential inside the microchannel is described using the Laplace equation while the microparticles' trajectory is described using equations based on Newton's second law. All equations are solved using finite difference method. The equations of motion account for forces including inertia, buoyancy, drag, gravity, virtual mass, and dielectrophoresis. The model is used for parametric study; the geometric parameters analyzed include microparticle radius, microchannel depth, and electrode/spacing lengths while volumetric flow rate and actuation voltage are the two operating parameters considered in the study. The trajectory of microparticles is composed of transient and steady state phases; the trajectory is influenced by all parameters. Microparticle radius and volumetric flow rate, above the threshold, do not influence the steady state levitation height; microparticle levitation is not possible below the threshold of the volumetric flow rate. Microchannel depth, electrode/spacing lengths, and actuation voltage influence the steady‐state levitation height.     

Liquid transport in rectangular microchannels by electroosmotic pumping

The characteristics of electroosmotic flow in rectangular microchannels were investigated in this paper. A 2D Poisson–Boltzmann equation and the 2D momentum equation were used to model the electric double layer field and the flow field. The numerical solutions show significant influences of the channel cross-section geometry (i.e. the aspect ratio) on the velocity field and the volumetric flow rate. Also, the numerical simulation of the electroosmotic flow reveals how the velocity field and the volumetric flow rate depend on the ionic concentration, zeta potential, channel size and the applied electrical field strength.     

DNA capture into a nanopore: interplay of diffusion and electrohydrodynamics

We present a detailed analysis of the process of voltage driven capture of DNA molecules by nanopores. We show that ionic current generates a nonuniform electric field that acts on both the DNA and on its counterions and that the response of DNA to the electric field is affected by its electroosmotic coupling to the mobile counterions. We calculate the voltage and molecular mass dependence of the radius of capture and of the capture rate in the diffusion limited regime. We argue that electroosmotic flow through the DNA coil is suppressed in the vicinity of the pore and present a tentative estimate of the capture rate in the barrier limited regime.     

3-D transient electrophoretic motion of a spherical particle in a T-shaped rectangular microchannel

This paper considers the electrophoretic motion of a spherical particle in an aqueous electrolyte solution in a T-shaped rectangular microchannel, where the size of the channel is close to that of the particle. This is a complicated transient process where the electric field, the flow field, and the particle motion are coupled together. A theoretical model was developed to investigate the influences of the applied electric potentials, the zeta potentials of the channel and the particle, and the size of the particle on the particle motion. A direct numerical simulation method using the finite element method is employed. This method employs a generalized Galerkin finite element formulation that incorporates both equations of the fluid flow and equations of the particle motion into a single variational equation where the hydrodynamic interactions are eliminated. The ALE method is used to track the surface of the particle at each time step. The numerical results show that the electric field in the T-shaped microchannel is influenced by the presence of the particle, and that the particle motion is influenced by the applied electric potentials and the zeta potentials of the channel and the particle. The path of the particle motion is dominated by the local electric field and the ratio of the zeta potential of the channel to that of the particle. The particle's velocity is also dependent on its size in a small channel.     

Electrokinetic transport through nanochannels

This article presents a numerical study of the electrokinetic transport phenomena (electroosmosis and electrophoresis) in a three-dimensional nanochannel with a circular cross-section. Due to the nanometer dimensions, the Boltzmann distribution of the ions is not valid in the nanochannels. Therefore, the conventional theories of electrokinetic flow through the microchannels such as Poisson-Boltzmann equation and Helmholtz-Smoluchowski slip velocity approach are no longer applicable. In the current study, a set of coupled partial differential equations including Poisson-Nernst-Plank equation, Navier-Stokes, and continuity equations is solved to find the electric potential field, ionic concentration field, and the velocity field in the three-dimensional nanochannel. The effects of surface electric charge and the radius of nanochannel on the electric potential, liquid flow, and ionic transport are investigated. Unlike the microchannels, the electric potential field, ionic concentration field, and velocity field are strongly size-dependent in nanochannels. The electric potential gradient along the nanochannel also depends on the surface electric charge of the nanochannel. More counter ions than the coions are transported through the nanochannel. The ionic concentration enrichment at the entrance and the exit of the nanochannel is completely evident from the simulation results. The study also shows that the flow velocity in the nanochannel is higher when the surface electric charge is stronger or the radius of the nanochannel is larger.     

A study of the transport of ions against their concentration gradient across ion-exchange membranes using the network method

The network method has been used to analyze the conditions that favour the uphill transport across ion-exchange membranes. A model for the Nernst-Planck-Poisson equations describing the ionic transport in such system is proposed, including the Donnan equilibrium relations at the membrane/solution interfaces. With this model and the electric circuit simulation program PSPICE, the transient response of the system under open circuit conditions (I=0) and the response of the system subject to an applied potential difference are simulated. The ionic concentrations and electric potential profiles, as well as the electric current density, the ionic fluxes and the charge density, have been obtained as a function of time.