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.

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.

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.

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.

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.

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.

Coarse-Grained Models of Aqueous Solutions of Polyelectrolytes: Significance of Explicit Charges.

The structure of polyelectrolytes is highly sensitive to small changes in interactions between their monomers. In particular, interactions mediated by counterions play a significant role and are affected by both specific molecular effects and generic concentration effects. The ability of coarse-grained models to reproduce the structural properties of an atomic model is thus a challenging task. Our present study compares the ability of different kinds of coarse-grained models: (i) to reproduce the structure of an atomistic model of a polyelectrolyte (the sodium polyacrylate) and (ii) to reproduce the variations of this structure with the number of monomers and with the concentration of different species. We show that adequate scalings of the gyration radius of the polymer Rg with the number of monomers N and with the box size Lbox are only obtained, first, if the monomer charges and the counterions are explicitly described and, second, if an attractive Lennard-Jones contribution is added to the interaction between distant monomers. Also, we show that implicit ion models are relevant only to the high electrostatic screening regime.

Computation of the Hydrodynamic Radius of Charged Nanoparticles from Nonequilibrium Molecular Dynamics.

We have used nonequilibrium molecular dynamics to simulate the flow of water molecules around a charged nanoparticle described at the atomic scale. These nonequilibrium simulations allowed us to compute the friction coefficient of the nanoparticle and then to deduce its hydrodynamic radius. We have compared two different strategies to thermostat the simulation box, since the low symmetry of the flow field renders the control of temperature non trivial. We show that both lead to an adequate control of the temperature of the system. To deduce the hydrodynamic radius of the nanoparticle we have employed a partial thermostat, which exploits the cylindrical symmetry of the flow field. Thereby, only a part of the simulation box far from the nanoparticle is thermostated. We have taken into account the finite concentration of the nanoparticle when calculating the friction force acting on it. We have focused on the case of polyoxometalate ions, which are inorganic charged nanoparticles. It appears that, for a given structure of the nanoparticle at the atomic level, the hydrodynamic radius significantly increases with the nanoparticles charge, a phenomenon that had not been quantified so far using molecular dynamics. The presence of an added salt only slightly modifies the hydrodynamic radius.

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.

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.

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.

Self-diffusion and activity coefficients of ions in charged disordered media.

Self-diffusion and single ion activity coefficients of ions of size symmetric electrolytes were studied in the presence of a collection of charged obstacles (called matrix) within a "soft" version of the primitive model of electrolyte solutions. The matrix subsystem possesses a net charge, depending on the concentration and charge of obstacles. The brownian dynamics method was used to calculate the self-diffusion coefficients of mobile species. The replica Ornstein-Zernike theory for the partly quenched systems was used to calculate the individual activity coefficients of mobile ionic species. The results reflect the competition between attractive (obstacle-counterion, co-ion-counterion), and repulsive (obstacle-co-ion) interactions in these model systems. For the simplest possible system of symmetric monovalent ions the latter effect wins: Co-ions are excluded from the area around obstacles, and this slows down their diffusion compared to that of counterions. Introduction of divalent charges into the system changes this result when the concentration of obstacles is low. We compare these results to those obtained for the corresponding fully annealed systems, i.e., where all the species are mobile. In most cases the self-diffusion and activity coefficients of counterions and co-ions in the presence of charged obstacles follow the trends of the fully annealed solution, which are dictated by the composition of the mixture. In few situations, however, the presence of charged obstacles modifies these trends. Our study allows us to clearly identify the effects due to obstacles, and to separate them from those arising from the composition of the solution. In the case of charge and size symmetric systems, the results for the individual activity coefficients fully support the hypothesis of the "electrostatic excluded volume". Thermodynamic and dynamic results are consistent in explaining the behavior of the systems studied.

Effective interaction between charged nanoparticles and DNA.

We investigate the effective interaction mediated by salt ions between charged nanoparticles (NPs) and DNA. DNA is modeled as an infinite cylinder with a constant surface charge in an implicit solvent. Monte Carlo simulations are used to compute the free energy of the system described in the framework of the primitive model of electrolytes, which accounts for excluded volumes of salt ions. A mean-field Poisson-Boltzmann theory also allows us to compute the free energy and provides us with explicit formulae for its main characteristics (position and depth of the minimum). We intend here to identify the physical parameters that have a major impact on the NP-DNA interaction, in an attempt to evaluate physico-chemical properties which could play a role in genotoxicity or, which could be exploited for therapeutic use. Thus, we investigate the influence on the effective interaction of: the shape of the nanoparticle, the magnitude of the nanoparticle charge and its distribution, the value of the pH of the solution, the magnitude of Van der Waals interactions depending on the nature of the constitutive material of the NP (metal vs. dielectric). We show that for positively charged concave NPs the effective interaction is repulsive at short distance, so that it presents a minimum at distance from the DNA. This short-range repulsion is specific to indented particles and is a robust property that holds for a large range of materials and charge densities.

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.

Salt exclusion in charged porous media: a coarse-graining strategy in the case of montmorillonite clays.

We study the exclusion of salt from charged porous media (Donnan effect), by using a coarse-grained approach. The porous medium is a lamellar system, namely a Montmorillonite clay, in contact with a reservoir, which contains an electrolyte solution. We develop a specific coarse-graining strategy to investigate the structural properties of this system. Molecular simulations are used to calibrate a mesoscopic model of the clay micropore in equilibrium with a reservoir. Brownian Dynamics simulations are then used to predict the structure of ions in the pore and the amount of NaCl salt entering the pore as a function of the pore size (the distance L between clay surfaces) and of the electrolyte concentration in the reservoir. These results are also compared to the predictions of a Density Functional Theory, which takes into account the excluded volumes of ions. We show that the calibration of the mesoscopic model is a key point and has a strong influence on the result. We observe that the salt exclusion increases when kappaL decreases (where kappa is the inverse of the Debye length) and that this effect is modulated by the correlations between ions. Two different regimes are revealed. At low concentrations in the reservoir, we observe a regime controlled by electrostatics: the Coulomb attraction between ions increases the amount of salt in the interlayer space. On the opposite, at high concentrations in the reservoir, the excluded volume effect dominates.

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.

Self-diffusion coefficients of ions in the presence of charged obstacles.

The self-diffusion coefficient of ions of the charge- and size-symmetric +1:-1 (or +2:-2) electrolyte was studied in the presence of ionic obstacles (matrix) representing disordered media. For this purpose the Brownian dynamics method was used, complemented with the replica Ornstein-Zernike theory for the partly-quenched systems. The matrix was prepared by a rapid quench of the size-symmetric +1:-1 (in few cases also of +2:-2) electrolyte solution being in equilibrium at (temperature, relative permittivity) T0, epsilon. Within the matrix the charge- and size-symmetric (+1:-1 or +2:-2) electrolyte at T1, epsilon1 was distributed. This component was fully mobile (annealed) and in thermodynamic equilibrium with the matrix. In this study a special attention was paid to the self-diffusion of the annealed ions. The ratio D/D degrees, where D degrees is the self-diffusion coefficient of ions at infinite dilution, has been studied for various model parameters varying the concentration of all species in the system. The presence of charged obstacles decreases the self-diffusion of the annealed electrolyte; the D/D degrees values are lower in the partly-quenched mixtures than in the fully annealed electrolyte of the same concentration. In the investigated range of concentrations and solvent dielectric constants, the D/D degrees values first increased with the increased concentration of annealed electrolyte present and then decreased. An increase of the strength of the Coulomb interaction between annealed ions, or between annealed and quenched charges, yielded a decrease of the self-diffusion. In the range of concentrations investigated in this work, the decrease is mainly due to the Coulomb interaction with the matrix, since the presence of neutral obstacles did not modify the diffusion properties with respect to the situation without obstacles.