Analytic Theory for Dilute Colloids in a Charged Slit

Electrostatics of Nanoparticle-Wall Interactions within Nanochannels: Role of Double-Layer Structure and Ion-Ion Correlations

We perform computational investigations of the electrolyte-mediated interactions of charged nanoparticles with the walls of nanochannels. We investigate the role of discrete ion effects, valence, and electrolyte strength on nanoparticle-wall interactions. We find for some of the multivalent charge regimes that the like-charged nanoparticles and walls can have attractive interactions. We study in detail these interactions and the free-energy profile for the nanoparticle-wall separation. We find there are energy barriers and energy minima giving preferred nanoparticle locations in the channel near the center and at a distance near to but separated from the channel walls. We characterize contributions from surface overcharging, condensed layers, and overlap of ion double layers. We perform our investigations using coarse-grained particle-level simulations with Brownian dynamics, classical density functional theory, and the mean-field Poisson-Boltzmann theory. We discuss the implications of our results for phenomena in nanoscale devices.

Ion correlations in nanofluidic channels: effects of ion size, valence, and concentration on voltage- and pressure-driven currents

The effects of ion-ion and ion-wall correlations in nanochannels are explored, specifically how they influence voltage- and pressure-driven currents and pressure-to-voltage energy conversion. Cations of different diameters (0.15, 0.3, and 0.9 nm) and different valences (+1, +2, and +3) at concentrations ranging from 10(-6) M to 1 M are considered in 50-nm- and 100-nm-wide nanoslits with wall surface charges ranging from 0 C/m(2) to -0.3 C/m(2). These parameters are typical of nanofluidic devices. Ion correlations have significant effects on device properties over large parts of this parameter space. These effects are the result of ion layering (oscillatory concentration profiles) for large monovalent cations and charge inversion (more cations in the first layer near the wall than necessary to neutralize the surface charge) for the multivalent cations. The ions were modeled as charged, hard spheres using density functional theory of fluids, and current was computed with the Navier-Stokes equations with two different no-slip conditions.

Energy variational analysis of ions in water and channels: Field theory for primitive models of complex ionic fluids

Ionic solutions are mixtures of interacting anions and cations. They hardly resemble dilute gases of uncharged noninteracting point particles described in elementary textbooks. Biological and electrochemical solutions have many components that interact strongly as they flow in concentrated environments near electrodes, ion channels, or active sites of enzymes. Interactions in concentrated environments help determine the characteristic properties of electrodes, enzymes, and ion channels. Flows are driven by a combination of electrical and chemical potentials that depend on the charges, concentrations, and sizes of all ions, not just the same type of ion. We use a variational method EnVarA (energy variational analysis) that combines Hamilton's least action and Rayleigh's dissipation principles to create a variational field theory that includes flow, friction, and complex structure with physical boundary conditions. EnVarA optimizes both the action integral functional of classical mechanics and the dissipation functional. These functionals can include entropy and dissipation as well as potential energy. The stationary point of the action is determined with respect to the trajectory of particles. The stationary point of the dissipation is determined with respect to rate functions (such as velocity). Both variations are written in one Eulerian (laboratory) framework. In variational analysis, an "extra layer" of mathematics is used to derive partial differential equations. Energies and dissipations of different components are combined in EnVarA and Euler-Lagrange equations are then derived. These partial differential equations are the unique consequence of the contributions of individual components. The form and parameters of the partial differential equations are determined by algebra without additional physical content or assumptions. The partial differential equations of mixtures automatically combine physical properties of individual (unmixed) components. If a new component is added to the energy or dissipation, the Euler-Lagrange equations change form and interaction terms appear without additional adjustable parameters. EnVarA has previously been used to compute properties of liquid crystals, polymer fluids, and electrorheological fluids containing solid balls and charged oil droplets that fission and fuse. Here we apply EnVarA to the primitive model of electrolytes in which ions are spheres in a frictional dielectric. The resulting Euler-Lagrange equations include electrostatics and diffusion and friction. They are a time dependent generalization of the Poisson-Nernst-Planck equations of semiconductors, electrochemistry, and molecular biophysics. They include the finite diameter of ions. The EnVarA treatment is applied to ions next to a charged wall, where layering is observed. Applied to an ion channel, EnVarA calculates a quick transient pile-up of electric charge, transient and steady flow through the channel, stationary "binding" in the channel, and the eventual accumulation of salts in "unstirred layers" near channels. EnVarA treats electrolytes in a unified way as complex rather than simple fluids. Ad hoc descriptions of interactions and flow have been used in many areas of science to deal with the nonideal properties of electrolytes. It seems likely that the variational treatment can simplify, unify, and perhaps derive and improve those descriptions.