Numerical study of density functional theory with mean spherical approximation for ionic condensation in highly charged confined electrolytes.

We investigate numerically a density functional theory (DFT) for strongly confined ionic solutions in the canonical ensemble by comparing predictions of ionic concentration profiles and pressure for the double-layer configuration to those obtained with Monte Carlo (MC) simulations and the simpler Poisson-Boltzmann (PB) approach. The DFT consists of a bulk (ion-ion) and an ion-solid part. The bulk part includes nonideal terms accounting for long-range electrostatic and short-range steric correlations between ions and is evaluated with the mean spherical approximation and the local density approximation. The ion-solid part treats the ion-solid interactions at the mean-field level through the solution of a Poisson problem. The main findings are that ionic concentration profiles are generally better described by PB than by DFT, although DFT captures the nonmonotone co-ion profile missed by PB. Instead, DFT yields more accurate pressure predictions than PB, showing in particular that nonideal effects are important to describe highly confined ionic solutions. Finally, we present a numerical methodology capable of handling nonconvex minimization problems so as to explore DFT predictions when the reduced temperature falls below the critical temperature.

Molecular simulation of aqueous solutions at clay surfaces.

We report a molecular simulation study of aqueous solutions at montmorillonite clay surfaces. Unlike most previous studies, ours does not focus on the interlayer nanopores, but looks at both kinds of external surfaces of clay particles: basal surfaces along the clay layers, and lateral surfaces through which interlayer and larger interparticle pores are linked. We present results on structural, dynamic and thermodynamic properties and phenomena, including hydration complexes of ions, H bonding networks, modification of the water dynamics with respect to the bulk, and the role of water in the cation exchange between interlayer and interparticle pores.

Two-scale Brownian dynamics of suspensions of charged nanoparticles including electrostatic and hydrodynamic interactions.

We propose here a multiscale strategy based on continuous solvent Brownian dynamics (BD) simulations to study the dynamical properties of aqueous suspensions of charged nanoparticles. We extend our previous coarse-graining strategy [V. Dahirel et al., J. Chem. Phys. 126, 114108 (2007)] to account for hydrodynamic interactions between solute particles. Within this new procedure, two BD simulations are performed: (1) The first one investigates the time scales of the counterions and coions (the microions) with only one nanoparticle in the simulation box but explicit microions, (ii) the second one investigates the larger time scale of the nanoparticles with numerous nanoparticles in the simulation box but implicit microions. We show how individual and collective transport coefficients can be computed from this two-scale procedure. To ensure the validity of our procedure, we compute the transport coefficients of a 10-1 model electrolyte in aqueous solution with a 1-1 added salt. We do a systematic comparison between the results obtained within the new procedure and those obtained with explicit BD simulations of the complete system containing several nanoparticles and explicit microions. The agreement between the two methods is found to be excellent: Even if the new procedure is much faster than explicit simulations, it allows us to compute transport coefficients with a good precision. Moreover, one step of our procedure also allows us to compute the individual transport coefficients relative to the microions (self-diffusion coefficients and electrophoretic mobility).

Signature of low-dimensional diffusion in complex systems.

We present a clear signature of the dimensionality of water diffusion in a powder sample of a synthetic hectorite (a model clay), by analyzing the corresponding neutron scattering functions. The data follow the theoretical predictions for a powder-averaged two-dimensional diffusion, with a two-dimensional diffusion coefficient of 0.75 x 10(-9) m2 s(-1). Neutron scattering data of bulk water are used as a reference, representing motion in three dimensions. The approach is based on analyzing the scattered intensity at zero energy transfers, along with the broadening of the scattering functions, collected at a wide range of energy resolutions. The mathematical relationship between these two quantities follows, for a given shape of the resolution function, a universal master curve, independent of the diffusion coefficient, but strongly dependent on the dimensionality of the motion, which can thus be determined with clarity.

Ion-mediated interactions between charged and neutral nanoparticles.

Monte-Carlo simulations are used to study the ion-mediated effective interaction between weakly charged and highly charged nanoparticles in an implicit solvent. Three models of nanoparticles are successively studied, from crude charged hard spheres to dipolar and non-spherical nanoparticles. The analysis of the effective potential revealed that in an electrolyte solution, even a neutral nanoparticle feels an important repulsive force in the presence of a charged nanoparticle, with a typical range similar to the Debye length. When the two nanoparticles carry charges of opposite sign, we have shown that this repulsion can reverse the effect of the direct attractive electrostatic potential at short distances. This also yields the change of sign of the effective potential as a function of the relative orientations of two anisotropic nanoparticles. Moreover, we found that the 3-body terms of the effective potentials can overcome the 2-body terms, which is not observed in the case of symmetrically charged nanoparticles.

Electrostatic relaxation and hydrodynamic interactions for self-diffusion of ions in electrolyte solutions.

The concentration dependence of self-diffusion of ions in solutions at large concentrations has remained an interesting yet unsolved problem. Here we develop a self-consistent microscopic approach based on the ideas of mode-coupling theory. It allows us to calculate both contributions which influence the friction of a moving ion: the ion atmosphere relaxation and hydrodynamic interactions. The resulting theory provides an excellent agreement with known experimental results over a wide concentration range. Interestingly, the mode-coupling self-consistent calculation of friction reveal a nonlinear coupling between the hydrodynamic interactions and the ion atmosphere relaxation which enhances ion diffusion by reducing friction, particularly at intermediate ion concentrations. This rather striking result has its origin in the similar time scales of the relaxation of the ion atmosphere relaxation and the hydrodynamic term, which are essentially given by the Debye relaxation time. The results are also in agreement with computer simulations, with and without hydrodynamic interactions.

Interpretation of conductivity results from 5 to 45 degrees C on three micellar systems below and above the CMC.

Electrical conductivity has been used at different temperatures to study three micellar systems: tetradecyltrimethylammonium chloride (TTACl), dodecyltrimethylammonium chloride (DTACl), and decyltrimethylammonium chloride (DeTACl). A phenomenon of premicellization is observed for DeTACl and DTACl below the critical micellar concentration (CMC). Association constants are introduced in the MSA-transport theory to correctly reproduce experimental conductivity and also calculate the effective charge of the micelles and their degree of dissociation. Various mechanisms are considered to explain premicellization. The formation of a neutral pair followed by an association involving two monomers and a counterion appears to be the most probable first step in the premicellization process.

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.

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.

Ion dynamics in compacted clays: derivation of a two-state diffusion-reaction scheme from the lattice Fokker-Planck equation.

We show how a two-state diffusion-reaction description of the mobility of ions confined within compacted clays can be constructed from the microscopic dynamics of ions in an external field. The diffusion-reaction picture provides the usual interpretation of the reduced ionic mobility in clays, but the required partitioning coefficient K(d) between trapped and mobile ions is generally an empirical parameter. We demonstrate that it is possible to obtain K(d) from the microscopic dynamics of ions interacting with the clay surfaces by evaluating the ionic mobility using a novel lattice implementation of the Fokker-Planck equation. The resulting K(d) allows a clear-cut characterization of the trapping sites on the clay surfaces and determines the adsorption/desorption rates. The results highlight the limitations of standard approximation schemes and pinpoint the crossover from jump to Brownian diffusion regimes.

Diffusion of water in clays on the microscopic scale: modeling and experiment.

Diffusion of water in montmorillonite clays at low hydration has been studied on the microscopic scale by two quasi-elastic neutron scattering techniques, neutron spin-echo (NSE) and time-of-flight (TOF), and by classical microscopic simulation. Experiment and simulation are compared both directly on the level of intermediate scattering functions, I(Q, t), and indirectly on the level of relaxation times after a model of atomic motion is applied. Regarding the dynamics of water in Na- and Cs-monohydrated montmorillonite samples, the simulation and NSE results show a very good agreement, both indicating diffusion coefficients of the order of (1-3) x 10(-10) m(2) s(-1). The TOF technique significantly underestimates water relaxation times (therefore overestimates water dynamics), by a factor of up to 3 and 7 in the two systems, respectively, primarily due to insufficiently long correlation times being probed. In the case of the Na-bihydrated system, the TOF results are in closer agreement with the other two techniques (the techniques differ by a factor of 2-3 at most), giving diffusion coefficients of (5-10) x 10(-10) m(2) s(-1). Attention has been also paid to the elastic incoherent structure factor, EISF(Q). Simulation has played a key role in understanding the various contributions to EISF(Q) in clay systems and in clearly distinguishing the signatures of "apparent" and true confinement. Indirectly, simulation highlights the difficulty in interpreting the EISF(Q) signal from powder clay samples used in experiments.

Molecular dynamics simulation of hydrogen fluoride mixtures with 1-ethyl-3-methylimidazolium fluoride: a simple model for the study of structural features.

Mixtures of hydrogen fluoride with ionic liquids show unique physicochemical properties, including their ability to form polyfluoride species (pointed out for the first time in this media by von Rosenvinge et al. J. Chem. Phys. 1997, 107, 8012). Among those systems the acidic 1-ethyl-3-methylimidazolium fluoride (EMIF.2.3HF) has been widely studied experimentally since it is the more promising for electrochemical applications. Recent studies (Hagiwara et al. J. Electrochem. Soc. 2002, 149, D1), while yielding many results, raised some questions about structural features of the liquid: absence of hydrogen bonds between the EMI+ ring hydrogen atoms and the fluoride anions, persistence of stacks and layers of cations similar to those existing in the crystal, and interpretation of the X-ray diffraction spectra. To address these questions, we have developed a simple molecular dynamics model. Our simulations are very consistent with experimental results and complete them, providing an atomic scale interpretation.

Frequency-dependent dielectric permittivity of salt-free charged lamellar systems.

We present a new model to analyze dielectric spectroscopy measurements on charged lamellar systems, with the following improvements with respect to the hitherto available models: (i) it does not rely on the hypothesis of local electro-neutrality, and allows to treat the salt-free case; (ii) the chemical exchange governing the partition between free and bound ions is properly taken into account; (iii) a fully analytical solution is provided. The variation of the frequency-dependent dielectric permittivity with both thermodynamic and kinetic characteristics of the free-bound ion equilibrium is presented. In particular, the relative weights of both relaxation modes (exchange and transport), and their characteristic frequencies are discussed. This study opens the way to the analysis of systems for which the usual models are irrelevant, such as salt-free clay gels or membranes.

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.

Study of individual Na-montmorillonite particles size, morphology, and apparent charge.

Size, morphology, and apparent charge of individual Na-montmorillonite particles of natural MX-80 sodium montmorillonite were investigated in the present study by the use of three coupling methods. In the first part of this work, natural and synthetic montmorillonite clays were studied with atomic force microscopy (AFM) and photo-correlation spectroscopy (PCS). Both techniques exhibit the presence of two clay populations with a high dispersion of the length distribution. Microscopic analysis of the system revealed that clay particles could be reasonably approximated at low concentrations to ellipsoidal tactoids about 1.2 nm high. Average dimensions of the first population were typically 320-400 nm long/250 nm wide and 200-250 nm long/120 nm wide for natural and synthetic clays, respectively. The second population exhibits smaller sizes: 65 and 50 nm long and 35 and 25 nm wide for natural and synthetic clays, respectively. The statistics obtained for natural clay were then verified by PCS experiments on sodium montmorillonite suspensions. Both techniques reveal an important length dispersion. However, the relative proportions of the two kinds of particles could not be established properly because of both lack of statistics and limitations of the employed techniques. In the following part, conductivity measurements were performed on dilute montmorillonite clay suspensions. Raw data were then interpreted with the sizes and morphological information gained in the first part of the present work. The apparent charge of the clay sheets was found to be 8% of the structural charge.

Analytical theories of transport in concentrated electrolyte solutions from the MSA.

Ion transport coefficients in electrolyte solutions (e.g., diffusion coefficients or electric conductivity) have been a subject of extensive studies for a long time. Whereas in the pioneering works of Debye, Hückel, and Onsager the ions were entirely characterized by their charge, recent theories allow specific effects of the ions (such as the ion size dependence or the pair association) to be obtained, both from simulation and from analytical theories. Such an approach, based on a combination of dynamic theories (Smoluchowski equation and mode-coupling theory) and of the mean spherical approximation (MSA) for the equilibrium pair correlation, is presented here. The various predicted equilibrium (osmotic pressure and activity coefficients) and transport coefficients (mutual diffusion, electric conductivity, self-diffusion, and transport numbers) are in good agreement with the experimental values up to high concentrations (1-2 mol L(-1)). Simple analytical expressions are obtained, and for practical use, the formula are given explicitly. We discuss the validity of such an approach which is nothing but a coarse-graining procedure.

An analytical model for probing ion dynamics in clays with broadband dielectric spectroscopy.

A simple two-state model is proposed to explicitly derive the ionic contribution to the frequency-dependent dielectric permittivity of clay. This model is based on a separation of time scales and accounts for two possible solvation modes (inner/outer-sphere complexes) for ions in the interlayer spacing and a possible chemical exchange between both forms. The influence on the permittivity of thermodynamic (distribution constant K(d)) and dynamic (diffusion coefficient, chemical relaxation rate) parameters is discussed. In turn, this model is used to analyze experimental data obtained with Na-montmorillonite for two relative humidities. The values of the parameters extracted from these measurements, and their variation with water content, show that the proposed model is at least reasonable.

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.

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.

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).