Direct Numerical Simulations of Electrophoresis of Charged Colloids

We propose a numerical method to simulate electrohydrodynamic phenomena in charged colloidal dispersions. This method enables us to compute the time evolutions of colloidal particles, ions, and host fluids simultaneously by solving Newton, advection-diffusion, and Navier--Stokes equations so that the electrohydrodynamic couplings can be fully taken into account. The electrophoretic mobilities of charged spherical particles are calculated in several situations. The comparisons with approximation theories show quantitative agreements for dilute dispersions without any empirical parameters, however, our simulation predicts notable deviations in the case of dense dispersions.Comment: 4pages, 3figures, to appear in Phys. Rev. Let

Electrophoretic mobility of a charged colloidal particle: A computer simulation study

We study the mobility of a charged colloidal particle in a constant homogeneous electric field by means of computer simulations. The simulation method combines a lattice Boltzmann scheme for the fluid with standard Langevin dynamics for the colloidal particle, which is built up from a net of bonded particles forming the surface of the colloid. The coupling between the two subsystems is introduced via friction forces. In addition explicit counterions, also coupled to the fluid, are present. We observe a non-monotonous dependence of the electrophoretic mobility on the bare colloidal charge. At low surface charge density we observe a linear increase of the mobility with bare charge, whereas at higher charges, where more than half of the ions are co-moving with the colloid, the mobility decreases with increasing bare charge.Comment: 15 pages, 8 figure

Colloidal electrophoresis: Scaling analysis, Green-Kubo relation, and numerical results

We consider electrophoresis of a single charged colloidal particle in a finite box with periodic boundary conditions, where added counterions and salt ions ensure charge neutrality. A systematic rescaling of the electrokinetic equations allows us to identify a minimum set of suitable dimensionless parameters, which, within this theoretical framework, determine the reduced electrophoretic mobility. It turns out that the salt-free case can, on the Mean Field level, be described in terms of just three parameters. A fourth parameter, which had previously been identified on the basis of straightforward dimensional analysis, can only be important beyond Mean Field. More complicated behavior is expected to arise when further ionic species are added. However, for a certain parameter regime, we can demonstrate that the salt-free case can be mapped onto a corresponding system containing additional salt. The Green-Kubo formula for the electrophoretic mobility is derived, and its usefulness demonstrated by simulation data. Finally, we report on finite-element solutions of the electrokinetic equations, using the commercial software package COMSOL.Comment: To appear in Journal of Physics: Condensed Matter - special issue on occasion of the CODEF 2008 conferenc

Numerical electrokinetics

A new lattice method is presented in order to efficiently solve the electrokinetic equations, which describe the structure and dynamics of the charge cloud and the flow field surrounding a single charged colloidal sphere, or a fixed array of such objects. We focus on calculating the electrophoretic mobility in the limit of small driving field, and systematically linearise the equations with respect to the latter. This gives rise to several subproblems, each of which is solved by a specialised numerical algorithm. For the total problem we combine these solvers in an iterative procedure. Applying this method, we study the effect of the screening mechanism (salt screening vs. counterion screening) on the electrophoretic mobility, and find a weak non-trivial dependence, as expected from scaling theory. Furthermore, we find that the orientation of the charge cloud (i. e. its dipole moment) depends on the value of the colloid charge, as a result of a competition between electrostatic and hydrodynamic effects.Comment: accepted for publication in Journal of Physics Condensed Matter (proceedings of the 2012 CODEF conference

Electronic interactions in fullerene spheres

The electron-phonon and Coulomb interactions inC$_{60}$, and larger fullerene spheres are analyzed. The coupling between electrons and intramolecular vibrations give corrections $\sim 1 - 10$ meV to the electronic energies for C$_{60}$, and scales as $R^{-4}$ in larger molecules. The energies associated with electrostatic interactions are of order $\sim 1 - 4$ eV, in C$_{60}$ and scale as $R^{-1}$. Charged fullerenes show enhanced electron-phonon coupling, $\sim 10$ meV, which scales as $R^{-2}$. Finally, it is argued that non only C$_{60}^{-}$, but also C$_{60}^{--}$ are highly polarizable molecules. The polarizabilities scale as $R^3$ and $R^4$, respectively. The role of this large polarizability in mediating intermolecular interactions is also discussed.Comment: 12 pages. No figure

Structural and Electronic Instabilities in Polyacenes: Density Matrix Renormalization Group Study of a Long--Range Interacting Model

We have carried out Density Matrix Renormalization Group (DMRG) calculations on the ground state of long polyacene oligomers within a Pariser-Parr-Pople (PPP) Hamiltonian. The PPP model includes long-range electron correlations which are required for physically realistic modeling of conjugated polymers. We have obtained the ground state energy as a function of the dimerization $\delta$ and various correlation functions and structure factors for $\delta=0$. From energetics, we find that while the nature of the Peierls' instabilityin polyacene is conditional and strong electron correlations enhance the dimerization. The {\it cis} form of the distortion is favoured over the {\it trans} form. However, from the analysis of correlation functions and associated structure factors, we find that polyacene is not susceptible to the formation of a bond order wave (BOW), spin density wave (SDW) or a charge density wave (CDW) in the ground state.Comment: 31 pages, latex, 13 figure

Selective trapping of DNA using glass microcapillaries

We show experimentally that a cheap glass microcapillary can accumulate {\lambda}-phage DNA at its tip and deliver the DNA into the capillary using a combination of electro-osmotic flow, pressure-driven flow, and electrophoresis. We develop an efficient simulation model for this phenomenon based on the electrokinetic equations and the finite-element method. Using our model, we explore the large parameter space of the trapping mechanism by varying the salt concentration, the capillary surface charge, the applied voltage, the pressure difference, and the mobility of the analyte molecules. Our simulation results show that this system can be tuned to capture a wide range of analyte molecules, such as DNA or proteins, based on their electrophoretic mobility. Our method for separation and pre-concentration of analytes has implications for the development of low-cost lab-on-a-chip devices.Comment: 9 pages, 4 figure

A Dissipative-Particle-Dynamics Model for Simulating Dynamics of Charged Colloid

A mesoscopic colloid model is developed in which a spherical colloid is represented by many interacting sites on its surface. The hydrodynamic interactions with thermal fluctuations are taken accounts in full using Dissipative Particle Dynamics, and the electrostatic interactions are simulated using Particle-Particle-Particle Mesh method. This new model is applied to investigate the electrophoretic mobility of a charged colloid under an external electric field, and the influence of salt concentration and colloid charge are systematically studied. The simulation results show good agreement with predictions from the electrokinetic theory.Comment: 17 pages, 8 figures, submitted to the proceedings of High Performance Computing in Science & Engineering '1

Reducing spurious flow in simulations of electrokinetic phenomena

Electrokinetic transport phenomena can strongly influence the behaviour of macromolecules and colloidal particles in solution, with applications in, e.g., DNA translocation through nanopores, electro-osmotic flow in nanocapillaries, and electrophoresis of charged macromolecules. Numerical simulations are an important tool to investigate these electrokinetic phenomena, but are often plagued by spurious fluxes and spurious flows that can easily exceed physical fluxes and flows. Here, we present a method that reduces one of these spurious currents, spurious flow, by several orders of magnitude. We demonstrate the effectiveness and generality of our method for both electrokinetic lattice-Boltzmann and finite-element-method based algorithms by simulating a charged sphere in an electrolyte solution, and flow through a nanopore. We also show that previous attempts to suppress these spurious currents introduce other sources of error.Comment: 13 pages, 7 figure