Isotropic active colloids: explicit vs. implicit descriptions of propulsion mechanisms

Modeling the couplings between active particles often neglects the possible many-body effects that control the propulsion mechanism. Accounting for such effects requires the explicit modeling of the molecular details at the origin of activity. Here, we take advantage of a recent two-dimensional model of isotropic active particles whose propulsion originates from the interactions between solute particles in the bath. The colloid catalyzes a chemical reaction in its vicinity, which results in a local phase separation of solute particles, and the density fluctuations of solute particles cause the enhanced diffusion of the colloid. In this paper, we investigate an assembly of such active particles, using (i) an explicit model, where the microscopic dynamics of the solute particles is accounted for; and (ii) an implicit model, whose parameters are inferred from the explicit model at infinite dilution. In the explicit solute model, the long-time diffusion coefficient of the active colloids strongly decreases with density, an effect which is not captured by the derived implicit model. This suggests that classical models, which usually decouple pair interactions from activity, fail to describe collective dynamics in active colloidal systems driven by solute-solute interactions.

Enhanced diffusion of tracer particles in nonreciprocal mixtures

We study the diffusivity of a tagged particle in a binary mixture of Brownian particles with nonreciprocal interactions. Numerical simulations reveal that, for a broad class of interaction potentials, nonreciprocity can significantly increase the long-time diffusion coefficient of tracer particles and that this diffusion enhancement is associated with a breakdown of the Einstein relation. These observations are quantified and confirmed via two different and complementary analytical approaches: (i) a linearized stochastic density field theory, which is particularly accurate in the limit of soft interactions, and (ii) a reduced two-body description, which is exact at leading order in the density of particles. The latter reveals that diffusion enhancement can be attributed to the formation of transiently propelled dimers of particles, whose cohesion and speed are controlled by the nonreciprocal interactions.

On analytical theories for conductivity and self-diffusion in concentrated electrolytes

Describing analytically the transport properties of electrolytes, such as their conductivity or the self-diffusion of the ions, has been a central challenge of chemical physics for almost a century. In recent years, this question has regained some interest in light of Stochastic Density Field Theory (SDFT) - an analytical framework that allows the approximate determination of density correlations in fluctuating systems. In spite of the success of this theory to describe dilute electrolytes, its extension to concentrated solutions raises a number of technical difficulties, and requires simplified descriptions of the short-range repulsion between the ions. In this article, we discuss recent approximations that were proposed to compute the conductivity of electrolytes, in particular truncations of Coulomb interactions at short distances. We extend them to another observable (the self-diffusion coefficient of the ions) and compare them to earlier analytical approaches, such as the mean spherical approximation and mode-coupling theory. We show how the treatment of hydrodynamic effects in SDFT can be improved, that the choice of the modified Coulomb interactions significantly affects the determination of the properties of the electrolytes, and that comparison with other theories provides a guide to extend SDFT approaches in this context.

Electrical surface properties of nanoporous alumina membranes: influence of nanochannels' curvature, roughness and composition studied via electrokinetic experiments

Among classical nanoporous oxide membranes, anodic aluminum oxide (AAO) membranes, made of non-connected, parallel and ordered nanochannels, are very interesting nanoporous model systems widely used for multiple applications. Since most of these applications involve local phenomena at the nanochannel surface, the fine description of the electrical surface behavior in aqueous solution is thus of primordial interest. Here, we use an original experimental approach combining several electrokinetic techniques (tangential and transverse streaming potential as well as electrophoretic mobility experiments) to measure the ζ-potential and determine the surface isoelectric points (IEPs) of several AAOs having different characteristic sizes and compositions. Using such an approach, all the different surfaces available in AAOs can be probed: outer surfaces (top and bottom planes), pore wall surfaces (i.e., inner surfaces) and surfaces created by the grinding of the AAOs. We find clear IEP differences between the outer, pore wall and ground surfaces and discuss these in terms of nanochannel and surface morphology (curvature and roughness) and of modifications of the chemical environment of the surface hydroxyl groups. These results highlight the heterogeneities between the different surfaces of these AAO membranes and emphasize the necessity to combine complementary electrokinetic techniques to properly understand the material, an approach which can be extended to many nanoporous systems.

Confined Electrolytes Show Bulk Dynamics Modulated by Hydrodynamic Couplings with the Walls

We use numerical simulations at the mesoscopic scale, namely, multiparticle collision dynamics (MPCD), to investigate the properties of electrolyte solutions in a charged slit pore. The solution is described within the primitive model of electrolytes, where ions are charged hard spheres embedded in a dielectric medium. Hydrodynamic couplings between ions and with the charged walls are precisely accounted for by the MPCD algorithm. We show that the dynamic properties of ions in this situation strongly differ from the behavior at infinite dilution (ideal case), contrary to what is usually assumed in the usual Poisson-Nernst-Planck description of this kind of systems. As a consequence of confinement, the diffusion coefficients of ions unexpectedly increase with the average ionic density in the systems. This is due to a decrease of the proportion of ions that are slowed down by the wall. Moreover, nonequilibrium simulations are used to estimate the electrical conductivity of these confined electrolytes. We show that the simulation results can be accounted for quantitatively by combining bulk descriptions of the electrical conductivity of electrolytes with a simple description of the hydrodynamics of ions in a slit pore.

Conditions for the propulsion of a colloid surrounded by a mesoscale phase separation

We study a two-dimensional model of an active isotropic colloid whose propulsion is linked to the interactions between solute particles of the bath. The colloid catalyzes a chemical reaction in its vicinity, that yields a local phase separation of solute particles. The density fluctuations of solute particles result in the enhanced diffusion of the colloid. Using numerical simulations, we thoroughly investigate the conditions under which activity occurs, and we establish a state diagram for the activity of the colloid as a function of the parameters of the model. We use the generated data to unravel a key observable that controls the existence and the intensity of activity: The filling fraction of the reaction area. Remarkably, we finally show that propulsion also occurs in three-dimensional geometries, which confirms the interest of this mechanism for experimental applications.

Diffusion of a tracer in a dense mixture of soft particles connected to different thermostats

We study the dynamics of a tracer in a dense mixture of particles connected to different thermostats. Starting from the overdamped Langevin equations that describe the evolution of the system, we derive the expression of the self-diffusion coefficient of a tagged particle in the suspension, in the limit of soft interactions between the particles. Our derivation, which relies on the linearization of the Dean-Kawasaki equations obeyed by the density fields and on a path-integral representation of the dynamics of the tracer, extends previous derivations that held for tracers in contact with a single bath. Our analytical result is confronted to results from Brownian dynamics simulations. The agreement with numerical simulations is very good even for high densities. We show how the diffusivity of tracers can be affected by the activity of a dense environment of soft particles that may represent polymer coils-a result that could be of relevance in the interpretation of measurements of diffusivity in biological media. Finally, our analytical result is general and can be applied to the diffusion of tracers coupled to different types of fluctuating environments, provided that their evolution equations are linear and that the coupling between the tracer and the bath is weak.

Electroosmotic Flow Induced Lift Forces on Polymer Chains in Nanochannels

A major objective of research in nanofluidics is to achieve better selectivity in manipulating the fluxes of nano-objects and in particular of biopolymers. Numerical simulations allow one to better understand the physical mechanisms at play in such situations. We performed hybrid mesoscale simulations to investigate the properties of polymers under flows in slit pores at the nanoscale. We use multiparticle collision dynamics, an algorithm that includes hydrodynamics and thermal fluctuations, to investigate the properties of fully flexible and stiff polymers under several types of flow, showing that Poiseuille flows and electroosmotic flows can lead to quantitatively and qualitatively different behaviors of the chain. In particular, a counterintuitive phenomenon occurs in the presence of an electroosmotic flow: When the monomers are attracted by the solid surfaces through van der Waals forces, shear-induced forces lead to a stronger repulsion of the polymers from these surfaces. Such focusing of the chain in the middle of the channel increases its flowing velocity, a phenomenon that may be exploited to separate different types of polymers.

Spontaneous propulsion of an isotropic colloid in a phase-separating environment

The motion of active colloids is generally achieved through their anisotropy, as exemplified by Janus colloids. Recently, there was a growing interest in the propulsion of isotropic colloids, which requires some local symmetry breaking. Although several mechanisms for such propulsion were proposed, little is known about the role played by the interactions within the environment of the colloid, which can have a dramatic effect on its propulsion. Here, we propose a minimal model of an isotropic colloid in a bath of solute particles that interact with each other. These interactions lead to a spontaneous phase transition close to the colloid, to directed motion of the colloid over very long timescales and to significantly enhanced diffusion, in spite of the crowding induced by solute particles. We determine the range of parameters where this effect is observable in the model, and we propose an effective Langevin equation that accounts for it and allows one to determine the different contributions at stake in self-propulsion and enhanced diffusion.

Comparison of different coupling schemes between counterions and charged nanoparticles in multiparticle collision dynamics

We applied the multiparticle collision dynamics (MPC) simulation technique to highly asymmetric electrolytes in solution, i.e., charged nanoparticles and their counterions in a solvent. These systems belong to a domain of solute size which ranges between the electrolyte and the colloidal domains, where most analytical theories are expected to fail, and efficient simulation techniques are still missing. MPC is a mesoscopic simulation method which mimics hydrodynamics properties of a fluid, includes thermal fluctuations, and can be coupled to a molecular dynamics of solutes. We took advantage of the size asymmetry between nanoparticles and counterions to treat the coupling between solutes and the solvent bath within the MPC method. Counterions were coupled to the solvent bath during the collision step and nanoparticles either through a direct interaction force or with stochastic rotation rules which mimic stick boundary conditions. Moreover, we adapted the simulation procedure to address the issue of the strong electrostatic interactions between solutes of opposite charges. We show that the short-ranged repulsion between counterions and nanoparticles can be modeled by stochastic reflection rules. This simulation scheme is very efficient from a computational point of view. We have also computed the transport coefficients for various densities. The diffusion of counterions was found in one case to increase slightly with the volume fraction of nanoparticles. The deviation of the electric conductivity from the ideal behavior (solutes at infinite dilution without any direct interactions) is found to be strong.

Critical effect of pore characteristics on capillary infiltration in mesoporous films

Capillary phenomena governing the mass-transport (capillary filling, condensation/evaporation) has been experimentally investigated in around 20 different silica thin films exhibiting various porosities with pores dimension ranging from 2 to 200 nm. Films have been prepared by sol-gel chemistry combined with soft-templating approaches and controlled dip coating process. Environmental ellipsometric porosimetry combined with electronic microscopy were used to assess the porosity characteristics. Investigation of lateral capillary filling was performed by following the natural infiltration of water and ionic liquids at the edge of a sessile drop in open air or underneath a PDMS cover. The Washburn model was applied to the displacement of the liquid front within the films to deduce the kinetic constants. The role of the different capillary phenomena were discussed with respect to the porosity characteristics (porosity vol%, pore dimensions and constrictions). We show that correlation between capillary filling rate and pore dimensions is not straightforward. Generally, with a minimum of constrictions, faster filling is observed for larger pores. In the case of mesopores (<50 nm in diameter), the presence of bottle necks considerably slows down the infiltration rate. At such a small dimension, evaporation/capillary condensation dynamics, taking place at the meniscus inside the porosity, has to be considered to explain the transport mode. This fundamental study is of interest for applications involving liquids at the interface of mesoporous networks such as nanofluidics, purification, separation, water harvesting or heat transfer.

Nonideal effects in electroacoustics of solutions of charged particles: combined experimental and theoretical analysis from simple electrolytes to small nanoparticles

The electric signal induced by an ultrasonic wave in aqueous solutions of charged species is measured and modeled without any adjustable parameter.

Pore network model of electrokinetic transport through charged porous media

We introduce a method for the numerical determination of the steady-state response of complex charged porous media to pressure, salt concentration, and electric potential gradients. The macroscopic fluxes of solvent, salt, and charge are computed within the framework of the Pore Network Model (PNM), which describes the pore structure of the samples as networks of pores connected to each other by channels. The PNM approach is used to capture the couplings between solvent and ionic flows which arise from the charge of the solid surfaces. For the microscopic transport coefficients on the channel scale, we take a simple analytical form obtained previously by solving the Poisson-Nernst-Planck and Stokes equations in a cylindrical channel. These transport coefficients are upscaled for a given network by imposing conservation laws for each pores, in the presence of macroscopic gradients across the sample. The complex pore structure of the material is captured by the distribution of channel diameters. We investigate the combined effects of this complex geometry, the surface charge, and the salt concentration on the macroscopic transport coefficients. The upscaled numerical model preserves the Onsager relations between the latter, as expected. The calculated macroscopic coefficients behave qualitatively as their microscopic counterparts, except for the permeability and the electro-osmotic coupling coefficient when the electrokinetic effects are strong. Quantitatively, the electrokinetic couplings increase the difference between the macroscopic coefficients and the corresponding ones for a single channel of average diameter.

Dynamics of solutes with hydrodynamic interactions: comparison between Brownian dynamics and stochastic rotation dynamics simulations

The dynamics of particles in solution or suspension is influenced by thermal fluctuations and hydrodynamic interactions. Several mesoscale methods exist to account for these solvent-induced effects such as Brownian dynamics with hydrodynamic interactions and hybrid molecular dynamics-stochastic rotation dynamics methods. Here we compare two ways of coupling solutes to the solvent with stochastic rotation dynamics (SRD) to Brownian dynamics with and without explicit hydrodynamic interactions. In the first SRD scheme [SRD with collisional coupling (CC)] the solutes participate in the collisional step with the solvent and in the second scheme [SRD with central force coupling (CFC)] the solutes interact through direct forces with the solvent, generating slip boundary conditions. We compare the transport coefficients of neutral and charged solutes in a model system obtained by these simulation schemes. Brownian dynamics without hydrodynamic interactions is used as a reference to quantify the influence of hydrodynamics on the transport coefficients as modeled by the different methods. We show that, in the dilute range, the SRD CFC method provides results similar to those of Brownian dynamics with hydrodynamic interactions for the diffusion coefficients and for the electrical conductivity. The SRD CC scheme predicts diffusion coefficients close to those obtained by Brownian dynamic simulations without hydrodynamic interactions, but accounts for part of the influence of hydrodynamics on the electrical conductivity.

Numerical homogenization of electrokinetic equations in porous media using lattice-Boltzmann simulations

We report the calculation of all the transfer coefficients which couple the solvent and ionic fluxes through a charged pore under the effect of pressure, electrostatic potential, and concentration gradients. We use a combination of analytical calculations at the Poisson-Nernst-Planck and Navier-Stokes levels of description and mesoscopic lattice simulations based on kinetic theory. In the absence of added salt, i.e., when the only ions present in the fluid are the counterions compensating the charge of the surface, exact analytical expressions for the fluxes in cylindrical pores allow us to validate a new lattice-Boltzmann electrokinetics (LBE) scheme which accounts for the osmotic contribution to the transport of all species. The influence of simulation parameters on the numerical accuracy is thoroughly investigated. In the presence of an added salt, we assess the range of validity of approximate expressions of the fluxes computed from the linearized Poisson-Boltzmann equation by a systematic comparison with LBE simulations.

Models of electrolyte solutions from molecular descriptions: the example of NaCl solutions

We present a method to derive implicit solvent models of electrolyte solutions from all-atom descriptions; providing analytical expressions of the thermodynamic and structural properties of the ions consistent with the underlying explicit solvent representation. Effective potentials between ions in solution are calculated to perform perturbation theory calculations, in order to derive the best possible description in terms of charged hard spheres. Applying this method to NaCl solutions yields excellent agreement with the all-atom model, provided ion association is taken into account.

Nonspecific DNA-protein interaction: why proteins can diffuse along DNA

Recent single molecule experiments have reported that DNA binding proteins (DNA-BPs) can diffuse along DNA. This suggests that interactions between proteins and DNA play a role during the target search even far from their specific site on DNA. Here we show by means of Monte Carlo simulations and analytical calculations that there is a counterintuitive repulsion between the two oppositely charged macromolecules at a nanometer range. For the concave shape of DNA-BPs, and for realistic protein charge densities, we find that the DNA-protein interaction free energy has a minimum at a finite surface-to-surface separation, in which proteins can easily slide. When a protein encounters its target, the free energy barrier is completely counterbalanced by the H-bond interaction, thus enabling the sequence recognition.

Brownian Simulations Contribution to the Study of Ionic Dynamics in Aqueous Solutions

We review some results concerning the dynamics of ionic solutions that we have recently obtained by Brownian dynamics simulations; we also present new results concerning the dynamical behavior of asymmetrical electrolytes (1–10 and 2–20 electrolytes). The Brownian dynamics treats the solution in the framework of the continuous solvent model, with a ‘soft-core’ version of the primitive model. Both direct interactions and hydrodynamic interactions between solutes are taken into account in the calculations. The method allows one to obtain the self-diffusion coefficient of each ion and the electrical conductivity of the solution. The computed transport coefficients are in good agreement with experimental determinations in various cases (in aqueous solutions of 1–1 electrolytes as well as in aqueous solutions of micelles). We show that hydrodynamic interactions must be taken into account to obtain electrical conductivities in agreement with experiments. The effect of hydrodynamic interactions on the self-diffusion is also striking, especially in solutions of asymmetrical electrolytes. If this effect remains weak in simple electrolyte solutions and for small ions in 1–20 and 2–20 electrolyte solutions, it is great for macroions in latter solutions (increase of about 15 to 40%).

Coarse-graining in suspensions of charged nanoparticles

A coarse-grain description of nanocolloidal suspensions in the presence of an added salt is presented here. It enables us to simulate trajectories of the nanoparticles from effective functions that depend on average densities of salt ions. In practice, the ion-averaged effective potential is used as input of a Brownian dynamics (BD) simulation. This potential may be derived by various methods, ranging from purely analytical to fully numerical ones. For the description of dynamical properties, this simulation also requires an effective diffusion coefficient that must be calculated or experimentally determined, and that accounts for the effects of microions on the mobility of the nanoparticles. The different versions of our coarse-graining procedure are applied to the case of a maghemite suspension, for which an explicit description of all ions would be very time-consuming.

How the excluded volume architecture influences ion-mediated forces between proteins

The effective interactions between model proteins of various shapes are computed by means of Monte Carlo simulations. In particular, we determine how the modification of the excluded volume architecture influences both entropic and purely electrostatic ion-mediated forces between proteins. We find that interprotein interactions are strongly affected by protein shape, which results in a high decrease of electrostatic screening for typical active site geometries. Effective interactions are then closer to the direct Coulombic interactions, and both affinity and selectivity are enhanced by several orders of magnitude.

Toward the description of electrostatic interactions between globular proteins: Potential of mean force in the primitive model

Monte Carlo simulations are used to calculate the exact potential of mean force between charged globular proteins in aqueous solution. The aim of the present paper is to study the influence of the ions of the added salt on the effective interaction between these nanoparticles. The charges of the model proteins, either identical or opposite, are either central or distributed on a discrete pattern. Contrarily to Poisson-Boltzmann predictions, attractive, and repulsive direct forces between proteins are not screened similarly. Moreover, it has been shown that the relative orientations of the charge patterns strongly influence salt-mediated interactions. More precisely, for short distances between the proteins, ions enhance the difference of the effective forces between (i) like-charged and oppositely charged proteins, (ii) attractive and repulsive relative orientations of the proteins, which may affect the selectivity of protein/protein recognition. Finally, such results observed with the simplest models are applied to a more elaborate one to demonstrate their generality.

New coarse-graining procedure for the dynamics of charged spherical nanoparticles in solution

A multiscale strategy based on the Brownian dynamics (BD) simulation method is presented here. It leads to an approximate but realistic reproduction of the dynamics of charged nanoparticles in suspension. This method is particularly suited to systems containing highly dissymmetric electrolytes with added salts, such as micellar suspensions or protein solutions. The coarse-graining procedure leads to a description where only the translational degrees of freedom of the nanoparticles are left, all the degrees of freedom related to the smallest solutes being rigorously averaged out. The authors' contribution aims at quantitatively evaluating the influence of the eliminated forces on the dynamics of the nanoparticles. For this purpose, an effective diffusion coefficient has to be calculated. In practice, this effective diffusion coefficient is taken as an input of a coarse-grained simulation that uses the potential of mean force between nanoparticles. The procedure has been validated by the quantitative comparison between the coarse-grained calculations and BD simulations at the "microscopic" level of description (which explicitly include microions). For a model of aqueous solutions of 10-1 electrolyte with a 1-1 added salt, the agreement is found to be excellent. This new method allows us to compute the diffusion coefficients of nanoparticles with a computation time at least one order of magnitude lower than with explicit BD.

Determining the Radius and the Apparent Charge of a Micelle from Electrical Conductivity Measurements by Using a Transport Theory: Explicit Equations for Practical Use

We propose here a procedure which combines experiments and simple analytical formulas that allows us to determine good estimations of the size and charge of ionic micelles above the critical micellar concentration (cmc). First, the conductivity of n-tetradecyltrimethylammonium bromide and chloride (TTABr and TTACl, respectively) aqueous solutions was measured at 25 degrees C, before and above their cmc. Then, an analytical expression for the concentration dependence of the conductance of an ionic mixture with three species (monomers, micelles, and counterions) was developed and applied to the analysis of the experiments. The theoretical calculations use the mean spherical approximation (MSA) to describe equilibrium properties. Here, we propose new expressions for the electrical conductivity, adapted to the case of electrolytes that are dissymmetric in size, and applicable up to a total surfactant concentration of 0.1 mol L(-1). Moreover, we show that they are good approximations of the corresponding numerical results obtained from Brownian dynamics simulations. Since the analytical formulas given in the present paper involve a small number of unknown parameters, they allow one to derive the size and charge of macroions in solution from conductivity measurements.

Brownian dynamics investigation of magnetization and birefringence relaxations in ferrofluids

Brownian dynamics simulations are used to investigate the dynamics of orientational properties of real charge-stabilized ferrofluids, i.e. stable colloidal dispersions of magnetic nanoparticles. The relaxation times of the magnetization and of the birefringence, data accessible by experimental techniques, have been computed at several volume fractions. Besides, the effect of hydrodynamic interactions has been investigated. Equilibrium simulations without field are found to be inadequate to determine the aforementioned relaxation times for the systems under study, the dipolar interaction being too weak. Thus a nonequilibrium simulation procedure that mimics the experimental operating mode has been developed. After equilibrium simulations under a magnetic field, both birefringence and magnetization decays are recorded once the field is suppressed. Birefringence and magnetization decays are markedly impeded as the volume fraction increases, whereas they are barely enhanced when the intensity of the initial magnetic field is increased at a fixed volume fraction. Eventually, hydrodynamic interactions exhibit a slight but systematic lengthening of the relaxation times.

Structural properties of charge-stabilized ferrofluids under a magnetic field: A Brownian dynamics study

We present Brownian dynamics simulations of real charge-stabilized ferrofluids, which are stable colloidal dispersions of magnetic nanoparticles, with and without the presence of an external magnetic field. The colloidal suspensions are treated as collections of monodisperse spherical particles, bearing point dipoles at their centers and undergoing translational and rotational Brownian motions. The overall repulsive isotropic interactions between particles, governed by electrostatic repulsions, are taken into account by a one-component effective pair interaction potential. The potential parameters are fitted in order that computed structure factors are close to the experimental ones. Two samples of ferrofluid differing by the particle diameter and consequently by the intensity of the magnetic interaction are considered here. The magnetization and birefringence curves are computed: a deviation from the ideal Langevin behaviors is observed if the dipolar moment of particles is sufficiently large. Structure factors are also computed from simulations with and without an applied magnetic field H: the microstructure of the repulsive ferrofluid becomes anisotropic under H. Even our simple modeling of the suspension allows us to account for the main experimental features: an increase of the peak intensity is observed in the direction perpendicular to the field whereas the peak intensity decreases in the direction parallel to the field.

Relaxation of the electrical double layer after an electron transfer approached by Brownian dynamics simulation

In this paper, the dynamical properties of the electrochemical double layer following an electron transfer are investigated by using Brownian dynamics simulations. This work is motivated by recent developments in ultrafast cyclic voltammetry which allow nanosecond time scales to be reached. A simple model of an electrochemical cell is developed by considering a 1:1 supporting electrolyte between two parallel walls carrying opposite surface charges, representing the electrodes; the solution also contains two neutral solutes representing the electroactive species. Equilibrium Brownian dynamics simulations of this system are performed. To mimic electron transfer processes at the electrode, the charge of the electroactive species are suddenly changed, and the subsequent relaxation of the surrounding ionic atmosphere are followed, using nonequilibrium Brownian dynamics. The electrostatic potential created in the center of the electroactive species by other ions is found to have an exponential decay which allows the evaluation of a characteristic relaxation time. The influence of the surface charge and of the electrolyte concentration on this time is discussed, for several conditions that mirror the ones of classical electrochemical experiments. The computed relaxation time of the double layer in aqueous solutions is found in the range 0.1 to 0.4 ns for electrolyte concentrations between 0.1 and 1 mol L(-1) and surface charges between 0.032 and 0.128 C m(-2).

Transport coefficients of aqueous dodecyltrimethylammonium bromide solutions: Comparison between experiments, analytical calculations and numerical simulations

We study dynamical properties of ionic species in aqueous solutions of dodecyltrimethylammonium bromide, for several concentrations below and above the critical micellar concentration (cmc). New experimental determinations of the electrical conductivity are given which are compared to results obtained from an analytical transport theory; transport coefficients of ions in these solutions above the cmc are also computed from Brownian dynamics simulations. Analytical calculations as well as the simulation treat the solution within the framework of the continuous solvent model. Above the cmc, three ionic species are considered: the monomer surfactant, the micelle and the counterion. The analytical transport theory describes the structural properties of the electrolyte solution within the mean spherical approximation and assumes that the dominant forces which determine the deviations of transport processes from the ideal behavior (i.e., without any interactions between ions) are hydrodynamic interactions and electrostatic relaxation forces. In the simulations, both direct interactions and hydrodynamic interactions between solutes are taken into account. The interaction potential is modeled by pairwise repulsive 1/r(12) interactions and Coulomb interactions. The input parameters of the simulation (radii and self-diffusion coefficients of ions at infinite dilution) are partially obtained from the analytical transport theory which fits the experimental determinations of the electrical conductivity. Both the electrical conductivity of the solution and the self-diffusion coefficients of each species computed from Brownian dynamics are compared to available experimental data. In every case, the influence of hydrodynamic interactions (HIs) on the transport coefficients is investigated. It is shown that HIs are crucial to obtain agreement with experiments. In particular, the self-diffusion coefficient of the micelle, which is the largest and most charged species in the present system, is enhanced when HIs are included whereas the diffusion coefficients of the monomer and the counterion are roughly not influenced by HIs.