Descreening of Field Effect in Electrically Gated Nanopores

This modeling work investigates the electrical modulation characteristics of field-effect gated nanopores. Highly nonlinear current modulations are observed in nanopores with non-overlapping electric double layers, including those with pore diameters 100 times the Debye screening length. We attribute this extended field-effect gating to a descreening effect, i.e. the counter-ions do not fully relax to screen the gating potential due to the presence of strong ionic transport

EH-DPD: a dissipative particle dynamics approach to electrohydrodynamics

Abstract: Electrohydrodynamics is crucial in many nanofluidic and biotechnological applications. In such small scales, the complexity due to the coupling of fluid dynamics with the dynamics of ions is increased by the relevance of thermal fluctuations. Here, we present a mesoscale method based on the Dissipative Particle Dynamics (DPD) model of the fluid. Two scalar quantities, corresponding to the number of positive and negative ions carried by each DPD particle, are added to the standard DPD formulation. We introduced a general framework that, given the definition of the free-energy of the DPD particle, allows to derive a fluctuation-dissipation relation and the expression for ionic fluxes between the DPD particles. This provides a link between the dynamics of the system and its equilibrium properties. The model is then validated simulating a planar electroosmotic flow for the cases of overlapping and non overlapping electric double layers. It is shown that using a Van der Waals equation of state the effect of ionic finite size can be accounted, leading to significant effects on the concentration and velocity profiles with respect to the ideal solution case. Graphic abstract: [Figure not available: see fulltext.]

Charge relaxation dynamics of an electrolytic nanocapacitor

Understanding ion relaxation dynamics in overlapping electric double layers (EDLs) is critical for the development of efficient nanotechnology based electrochemical energy storage, electrochemomechanical energy conversion and bioelectrochemical sensing devices as well as controlled synthesis of nanostructured materials. Here, a Lattice Boltzmann (LB) method is employed to simulate an electrolytic nanocapacitor subjected to a step potential at t = 0 for various degrees of EDL overlap, solvent viscosities, ratios of cation to anion diffusivity and electrode separations. The use of a novel, continuously varying and Galilean invariant, molecular speed dependent relaxation time (MSDRT) with the LB equation recovers a correct microscopic description of the molecular collision phenomena and enhances the stability of the LB algorithm. Results for large EDL overlaps indicated oscillatory behavior for the ionic current density in contrast to monotonic relaxation to equilibrium for low EDL overlaps. Further, at low solvent viscosities and large EDL overlaps, anomalous plasma-like spatial oscillations of the electric field were observed that appeared to be purely an effect of nanoscale confinement. Employing MSDRT in our simulations enabled a modeling of the fundamental physics of the transient charge relaxation dynamics in electrochemical systems operating away from equilibrium wherein Nernst-Einstein relation is known to be violated.Comment: Accepted for publication in the Journal of Physical Chemistry C on October 30 2014. Supplementary info available free of charge via the Internet at http://pubs.acs.org. Revised version includes more details on the computation of the molecular speed dependent relaxation time (MSDRT) and emphasizes the Galilean invariance of the computed MSDR

Nonlinear Dynamic Modeling, Simulation And Characterization Of The Mesoscale Neuron-electrode Interface

Extracellular neuroelectronic interfacing has important applications in the fields of neural prosthetics, biological computation and whole-cell biosensing for drug screening and toxin detection. While the field of neuroelectronic interfacing holds great promise, the recording of high-fidelity signals from extracellular devices has long suffered from the problem of low signal-to-noise ratios and changes in signal shapes due to the presence of highly dispersive dielectric medium in the neuron-microelectrode cleft. This has made it difficult to correlate the extracellularly recorded signals with the intracellular signals recorded using conventional patch-clamp electrophysiology. For bringing about an improvement in the signalto-noise ratio of the signals recorded on the extracellular microelectrodes and to explore strategies for engineering the neuron-electrode interface there exists a need to model, simulate and characterize the cell-sensor interface to better understand the mechanism of signal transduction across the interface. Efforts to date for modeling the neuron-electrode interface have primarily focused on the use of point or area contact linear equivalent circuit models for a description of the interface with an assumption of passive linearity for the dynamics of the interfacial medium in the cell-electrode cleft. In this dissertation, results are presented from a nonlinear dynamic characterization of the neuroelectronic junction based on Volterra-Wiener modeling which showed that the process of signal transduction at the interface may have nonlinear contributions from the interfacial medium. An optimization based study of linear equivalent circuit models for representing signals recorded at the neuron-electrode interface subsequently iv proved conclusively that the process of signal transduction across the interface is indeed nonlinear. Following this a theoretical framework for the extraction of the complex nonlinear material parameters of the interfacial medium like the dielectric permittivity, conductivity and diffusivity tensors based on dynamic nonlinear Volterra-Wiener modeling was developed. Within this framework, the use of Gaussian bandlimited white noise for nonlinear impedance spectroscopy was shown to offer considerable advantages over the use of sinusoidal inputs for nonlinear harmonic analysis currently employed in impedance characterization of nonlinear electrochemical systems. Signal transduction at the neuron-microelectrode interface is mediated by the interfacial medium confined to a thin cleft with thickness on the scale of 20-110 nm giving rise to Knudsen numbers (ratio of mean free path to characteristic system length) in the range of 0.015 and 0.003 for ionic electrodiffusion. At these Knudsen numbers, the continuum assumptions made in the use of Poisson-Nernst-Planck system of equations for modeling ionic electrodiffusion are not valid. Therefore, a lattice Boltzmann method (LBM) based multiphysics solver suitable for modeling ionic electrodiffusion at the mesoscale neuron-microelectrode interface was developed. Additionally, a molecular speed dependent relaxation time was proposed for use in the lattice Boltzmann equation. Such a relaxation time holds promise for enhancing the numerical stability of lattice Boltzmann algorithms as it helped recover a physically correct description of microscopic phenomena related to particle collisions governed by their local density on the lattice. Next, using this multiphysics solver simulations were carried out for the charge relaxation dynamics of an electrolytic nanocapacitor with the intention of ultimately employing it for a simulation of the capacitive coupling between the neuron and the v planar microelectrode on a microelectrode array (MEA). Simulations of the charge relaxation dynamics for a step potential applied at t = 0 to the capacitor electrodes were carried out for varying conditions of electric double layer (EDL) overlap, solvent viscosity, electrode spacing and ratio of cation to anion diffusivity. For a large EDL overlap, an anomalous plasma-like collective behavior of oscillating ions at a frequency much lower than the plasma frequency of the electrolyte was observed and as such it appears to be purely an effect of nanoscale confinement. Results from these simulations are then discussed in the context of the dynamics of the interfacial medium in the neuron-microelectrode cleft. In conclusion, a synergistic approach to engineering the neuron-microelectrode interface is outlined through a use of the nonlinear dynamic modeling, simulation and characterization tools developed as part of this dissertation research

Analysis of electrolyte transport through charged nanopores

We revisit the classical problem of flow of electrolyte solutions through charged capillary nanopores or nanotubes as described by the capillary pore model (also called “space charge” theory). This theory assumes very long and thin pores and uses a one-dimensional flux-force formalism which relates fluxes (electrical current, salt flux, and fluid velocity) and driving forces (difference in electric potential, salt concentration, and pressure). We analyze the general case with overlapping electric double layers in the pore and a nonzero axial salt concentration gradient. The 3 × 3 matrix relating these quantities exhibits Onsager symmetry and we report a significant new simplification for the diagonal element relating axial salt flux to the gradient in chemical potential. We prove that Onsager symmetry is preserved under changes of variables, which we illustrate by transformation to a different flux-force matrix given by Gross and Osterle [J. Chem. Phys. 49, 228 (1968)JCPSA60021-960610.1063/1.1669814]. The capillary pore model is well suited to describe the nonlinear response of charged membranes or nanofluidic devices for electrokinetic energy conversion and water desalination, as long as the transverse ion profiles remain in local quasiequilibrium. As an example, we evaluate electrical power production from a salt concentration difference by reverse electrodialysis, using an efficiency versus power diagram. We show that since the capillary pore model allows for axial gradients in salt concentration, partial loops in current, salt flux, or fluid flow can develop in the pore. Predictions for macroscopic transport properties using a reduced model, where the potential and concentration are assumed to be invariant with radial coordinate (“uniform potential” or “fine capillary pore” model), are close to results of the full model.Stanford University. Global Climate and Energy ProjectUnited States. Dept. of Energy. Office of Basic Energy Sciences. SUNCAT Center for Interface Science and Catalysi

Accurate Multi-physics Numerical Analysis of Particle Preconcentration Based on Ion Concentration Polarization

This paper studies mechanism of preconcentration of charged particles in a straight micro-channel embedded with permselective membranes, by numerically solving coupled transport equations of ions, charged particles and solvent fluid without any simplifying assumptions. It is demonstrated that trapping and preconcentration of charged particles are determined by the interplay between drag force from the electroosmotic fluid flow and the electrophoretic force applied trough the electric field. Several insightful characteristics are revealed, including the diverse dynamics of co-ions and counter ions, replacement of co-ions by focused particles, lowered ion concentrations in particle enriched zone, and enhanced electroosmotic pumping effect etc. Conditions for particles that may be concentrated are identified in terms of charges, sizes and electrophoretic mobilities of particles and co-ions. Dependences of enrichment factor on cross-membrane voltage, initial particle concentration and buffer ion concentrations are analyzed and the underlying reasons are elaborated. Finally, post priori a condition for validity of decoupled simulation model is given based on charges carried by focused charge particles and that by buffer co-ions. These results provide important guidance in the design and optimization of nanofluidic preconcentration and other related devices.Comment: 18 pages, 11 firgure

Observation of the induced pressure in a hybrid micro/nano-channel

This paper studies the flow characteristics in micro/nano-channels subjected to an applied electric field. The nano-channel flow was observed by means of the fluorescence Calcein. A Fluorescence Concentration Gradient Interface (FCGI) was observed across the nano-channel array. The front of the FCGI was shown to have an analogous parabolic shape. The propagation of this interface reflects indirectly the induced pressure at the micro/nano-channel junction, where the enrichment-depletion processes are known to take place. This induced pressure was predicted by numerical simulations, and this paper gives the first experimental evidence

Hydronium-dominated ion transport in carbon-dioxide-saturated electrolytes at low salt concentrations in nanochannels

Nanochannel ion transport is known to be governed by surface charge at low ionic concentrations. In this paper, we show that this surface charge is typically dominated by hydronium ions arising from dissolution of ambient atmospheric carbon dioxide. Taking the hydronium ions into account, we model the nanochannel conductance at low salt concentrations and identify a conductance minimum before saturation at a value independent of salt concentration in the dilute limit. Via the Poisson-Boltzmann equation, our model self-consistently couples chemical-equilibrium dissociation models of the silica wall and of the electrolyte bulk, parametrized by the dissociation reaction constants. Experimental data with aqueous KCl solutions in 165-nm-high silica nanochannels are described well by our model, both with and without extra hydronium from added HCl