Inflation of the screening length induced by Bjerrum pairs

Ions and dipoles in electric field: nonlinear polarization and field-dependent chemical reaction

We investigate electric-field effects in dilute electrolytes with nonlinear polarization. As a first example of such systems, we add a dipolar component with a relatively large dipole moment [Formula: see text] to an aqueous electrolyte. As a second example, the solvent itself exhibits nonlinear polarization near charged objects. For such systems, we present a Ginzburg-Landau free energy and introduce field-dependent chemical potentials, entropy density, and stress tensor, which satisfy general thermodynamic relations. In the first example, the dipoles accumulate in high-field regions, as predicted by Abrashikin et al.[Phys.Rev.Lett. 99, 077801 (2007)]. Finally, we consider the case, where Bjerrum ion pairs form a dipolar component with nonlinear polarization. The Bjerrum dipoles accumulate in high-field regions, while field-induced dissociation was predicted by Onsager [J. Chem. Phys.2, 599 (1934)]. We present an expression for the field-dependent association constant K(E), which depends on the field strength nonmonotonically.

Anomalous Underscreening in the Restricted Primitive Model

Underscreening is a collective term for charge correlations in electrolytes decaying slower than the Debye length. Anomalous underscreening refers to phenomenology that cannot be attributed alone to steric interactions. Experiments with concentrated electrolytes and ionic fluids report anomalous underscreening, which so far has not been observed in simulation. We present Molecular Dynamics simulation results exhibiting anomalous underscreening that can be connected to cluster formation. A theory that accounts for ion pairing confirms the trend. Our results challenge the classic understanding of dense electrolytes impacting the design of technologies for energy storage and conversion.

Theory of electrolytes including steric, attractive, and hydration interactions

We present a continuum theory of electrolytes composed of a waterlike solvent and univalent ions. First, we start with a density functional F for the coarse-grained solvent, cation, and anion densities, including the Debye-Hückel free energy, the Coulombic interaction, and the direct interactions among these three components. These densities fluctuate obeying the distribution ∝exp(-F/k_BT). Eliminating the solvent density deviation in F, we obtain the effective non-Coulombic interactions among the ions, which consist of the direct ones and the solvent-mediated ones. We then derive general expressions for the ion correlation, the apparent partial volume, and the activity and osmotic coefficients up to linear order in the average salt density n_s. Second, we perform numerical analysis using the Mansoori-Carnahan-Starling-Leland model [J. Chem. Phys. 54, 1523 (1971)] for three-component hardspheres. The effective interactions sensitively depend on the cation and anion sizes due to competition between the steric and hydration effects, which are repulsive between small-large ion pairs and attractive between symmetric pairs. These agree with previous experiments and Collins' rule [Biophys. J. 72, 65 (1997)]. We also give simple approximate expressions for the ionic interaction coefficients valid for any ion sizes.

Charge oscillations in ionic liquids: A microscopic cluster model

In spite of their enormous applications as alternative energy storage devices and lubricants, room-temperature ionic liquids (ILs) still pose many challenges from a pure scientific viewpoint. We develop an IL microscopic theory in terms of ionic clusters, which describes the IL behavior close to charged interfaces. The full structure factor of finite-size clusters is considered and allows us to retain fine and essential details of the system as a whole. Beside the reduction in the screening, it is shown that ionic clusters cause the charge density to oscillate near charged boundaries, with alternating ion-size thick layers, in agreement with experiments. We distinguish between short-range oscillations that persist for a few ionic layers close to the boundary, as opposed to long-range damped oscillations that hold throughout the bulk. The former can be captured by finite-size ion pairs, while the latter is associated with larger clusters with a pronounced quadrupole (or higher) moment. The long-wavelength limit of our theory recovers the well-known Bazant-Storey-Kornyshev (BSK) equation in the linear regime, and elucidates the microscopic origin of the BSK phenomenological parameters.

Screening length for finite-size ions in concentrated electrolytes

The classical Debye-Hückel (DH) theory clearly accounts for the origin of screening in electrolyte solutions and works rather well for dilute electrolyte solutions. While the Debye screening length decreases with the ion concentration and is independent of ion size, recent surface-force measurements imply that for concentrated solutions, the screening length exhibits an opposite trend; it increases with ion concentration and depends on the ionic size. The screening length is usually defined by the response of the electrolyte solution to a test charge but can equivalently be derived from the charge-charge correlation function. By going beyond DH theory, we predict the effects of ion size on the charge-charge correlation function. A simple modification of the Coulomb interaction kernel to account for the excluded volume of neighboring ions yields a nonmonotonic dependence of the screening length (correlation length) on the ionic concentration, as well as damped charge oscillations for high concentrations.

Dielectric constant of ionic solutions: Combined effects of correlations and excluded volume

The dielectric constant of ionic solutions is known to reduce with increasing ionic concentrations. However, the origin of this effect has not been thoroughly explored. In this paper, we study two such possible sources: long-range Coulombic correlations and solvent excluded-volume. Correlations originate from fluctuations of the electrostatic potential beyond the mean-field Poisson-Boltzmann theory, evaluated by employing a field-theoretical loop expansion of the free energy. The solvent excluded-volume, on the other hand, stems from the finite ion size, accounted for via a lattice-gas model. We show that both correlations and excluded volume are required in order to capture the important features of the dielectric behavior. For highly polar solvents, such as water, the dielectric constant is given by the product of the solvent volume fraction and a concentration-dependent susceptibility per volume fraction. The available solvent volume decreases as a function of ionic strength due the increasing volume fraction of ions. A similar decrease occurs for the susceptibility due to the correlations between the ions and solvent, reducing the dielectric response even further. Our predictions for the dielectric constant fit well with experiments for a wide range of concentrations for different salts in different temperatures, using a single fit parameter related to the ion size.

Bjerrum pairs in ionic solutions: A Poisson-Boltzmann approach

Ionic solutions are often regarded as fully dissociated ions dispersed in a polar solvent. While this picture holds for dilute solutions, at higher ionic concentrations, oppositely charged ions can associate into dimers, referred to as Bjerrum pairs. We consider the formation of such pairs within the nonlinear Poisson-Boltzmann framework and investigate their effects on bulk and interfacial properties of electrolytes. Our findings show that pairs can reduce the magnitude of the dielectric decrement of ionic solutions as the ionic concentration increases. We describe the effect of pairs on the Debye screening length and relate our results to recent surface-force experiments. Furthermore, we show that Bjerrum pairs reduce the ionic concentration in bulk electrolyte and at the proximity of charged surfaces, while they enhance the attraction between oppositely charged surfaces.

Cartilage-like electrostatic stiffening of responsive cryogel scaffolds

Cartilage is a structural tissue with unique mechanical properties deriving from its electrically-charged porous structure. Traditional three-dimensional environments for the culture of cells fail to display the complex physical response displayed by the natural tissue. In this work, the reproduction of the charged environment found in cartilage is achieved using polyelectrolyte hydrogels based on polyvinyl alcohol and polyacrylic acid. The mechanical response and morphology of microporous physically-crosslinked cryogels are compared to those of heat-treated chemical gels made from the same polymers, as a result of pH-dependent swelling. In contrast to the heat-treated chemically-crosslinked gels, the elastic modulus of the physical cryogels was found to increase with charge activation and swelling, explained by the occurrence of electrostatic stiffening of the polymer chains at large charge densities. At the same time, the permeability of both materials to fluid flow was impaired by the presence of electric charges. This cartilage-like mechanical behavior displayed by responsive cryogels can be reproduced in other polyelectrolyte hydrogel systems to fabricate biomimetic cellular scaffolds for the repair of the tissue.

Measured electrical charge of SiO2 in polar and nonpolar media

We present measurements of the net electrical surface charge of silicon dioxide (SiO₂) in contact with solvents of dielectric constants between 5 and 80. Our experimental approach relies on observing the thermal motion of single silica particles confined in an electrostatic fluidic trap created by SiO₂ surfaces. We compare the experimentally measured functional form of the trapping potential with that from free energy calculations and thereby determine the net surface charge in the system. Our findings clearly demonstrate that contrary to popular perception, even in the absence of surfactants, the net electrical charge of ionizable surfaces in contact with apolar solvents can be large enough to lead to significant repulsive forces. A charge regulation model for SiO₂ surfaces with a single tunable parameter explains our measurements. This model may find general applicability in estimating the net charge of ionizable surfaces, given system parameters such as the dissociation or association constants of the ionizable groups and the pH, ionic strength, and dielectric constant of the solvent phase.

Electrodynamics of lipid membrane interactions in the presence of zwitterionic buffers

Due to thermal motion and molecular polarizability, electrical interactions in biological systems have a dynamic character. Zwitterions are dipolar molecules that typically are highly polarizable and exhibit both a positive and a negative charge depending on the pH of the solution. We use multilamellar structures of common lipids to identify and quantify the effects of zwitterionic buffers that go beyond the control of pH. We use the fact that the repeat spacing of multilamellar lipid bilayers is a sensitive and accurate indicator of the force balance between membranes. We show that common buffers can in fact charge up neutral membranes. However, this electrostatic effect is not immediately recognized because of the concomitant modification of dispersion (van der Waals) forces. We show that although surface charging can be weak, electrostatic forces are significant even at large distances because of reduced ionic screening and reduced van der Waals attraction. The zwitterionic interactions that we identify are expected to be relevant for interfacial biological processes involving lipid bilayers, and for a wide range of biomaterials, including amino acids, detergents, and pharmaceutical drugs. An appreciation of zwitterionic electrodynamic character can lead to a better understanding of molecular interactions in biological systems and in soft materials in general.

Interaction regimes for oppositely charged plates with multivalent counterions

Within a mean-field treatment of the interaction between two oppositely charged plates in a salt-free solution, the distance at which a transition from an attractive to a repulsive regime appears can be computed analytically. The mean-field description, however, breaks down under strong Coulombic couplings, which can be achieved at room temperature with multivalent counterions and highly charged surfaces. Making use of the contact theorem and simple physical arguments, we propose explicit expressions for the equation of state in several situations at short distances. The possibility of Bjerrum pair formation is addressed and is shown to have profound consequences on the interactions. To complete the picture, we consider the large-distance limit, from which schematic phase diagram discriminating attractive from repulsive regions can be proposed.

Dressed counterions: polyvalent and monovalent ions at charged dielectric interfaces

We investigate the ion distribution and overcharging at charged interfaces with dielectric inhomogeneities in the presence of asymmetric electrolytes containing polyvalent and monovalent ions. We formulate an effective "dressed counterion" approach by integrating out the monovalent salt degrees of freedom and show that it agrees with results of explicit Monte Carlo simulations. We then apply the dressed counterion approach within the framework of the generalized strong-coupling theory, valid for polyvalent ions at low concentrations, which enables an analytical description for salt effects as well as dielectric inhomogeneities in the limit of strong Coulomb interactions. Limitations and applicability of this theory are examined by comparing the results with simulations.

Dressed counterions: strong electrostatic coupling in the presence of salt

We reformulate the theory of strong electrostatic coupling in order to describe an asymmetric electrolyte solution of monovalent salt ions and polyvalent counterions using field-theoretical techniques and Monte Carlo simulations. The theory is based on an asymmetric treatment of the different components of the electrolyte solution. The weak coupling Debye-Hückel approach is used in order to describe the monovalent salt ions while a strong coupling approach is used to tackle the polyvalent counterions. This combined weak-strong coupling approach effectively leads to dressed interactions between polyvalent counterions and thus directly affects the correlation attraction mediated by polyvalent counterions between like-charged objects. The general theory is specifically applied to a system composed of two uniformly charged plane-parallel surfaces in the presence of salt and polyvalent counterions. In the strong coupling limit for polyvalent counterions, the comparison with Monte Carlo simulations shows good agreement for large enough values of the electrostatic coupling parameter. We delineate two limiting laws that in fact encompass all the Monte Carlo data.