Effective interactions between electric double layers

A Gaussian field approach to the solvation of spherical ions in electrolyte solutions

In this work, the electrostatic response of an electrolyte solution to a spherical ion is studied with a Gaussian field theory. In order to capture the ionic correlation effect in concentrated solutions, the bulk dielectric response function is described by a two-Yukawa response function. The modified response function of the solution is solved analytically in the spherical geometry, from which the induced charge density and the electrostatic energy are also derived analytically. Comparisons with results for small ions in electrolyte solutions from the hyper-netted chain theory demonstrate the validity of the Gaussian field theory.

Ion correlation-driven like-charge attraction in multivalent salt solutions

The electrostatic double layer force is key to determining the stability and self-assembly of charged colloids and many other soft matter systems. Fully understanding the attractive force between two like-charged surfaces remains a great challenge. Here, we apply the modified Gaussian renormalized fluctuation theory to study ion correlation-driven like-charge attraction in multivalent salt solutions. The effects of spatially varying ion correlations on the structure of overlapping double layers and their free energy are self-consistently accounted for. In the presence of multivalent salts, increasing surface charge or counterion valency leads to a short-range attraction. We demonstrate that although both overcharging and like-charge attraction are outcomes of ion correlation, there is no causal relationship between them. Our theory also captures the non-monotonic dependence of like-charge attraction on multivalent salt concentration. The reduction of attraction at high salt concentrations could be a contributing factor toward the reentrant stability of charged colloidal suspensions. Our theoretical predictions are consistent with the observations reported in experiments and simulations.

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.

Effective interactions, structure, and pressure in charge-stabilized colloidal suspensions: Critical assessment of charge renormalization methods

Charge-stabilized colloidal suspensions display a rich variety of microstructural and thermodynamic properties, which are determined by electro-steric interactions between all ionic species. The large size asymmetry between molecular-scale microions and colloidal macroions allows the microion degrees of freedom to be integrated out, leading to an effective one-component model of microion-dressed colloidal quasi-particles. For highly charged colloids with strong macroion-microion correlations, nonlinear effects can be incorporated into effective interactions by means of charge renormalization methods. Here, we compare and partially extend several practical mean-field methods of calculating renormalized colloidal interaction parameters, including effective charges and screening constants, as functions of concentration and ionic strength. Within the one-component description, we compute structural and thermodynamic properties from the effective interactions and assess the accuracy of the different methods by comparing predictions with elaborate primitive-model simulations [P. Linse, J. Chem. Phys. 113, 4359 (2000)]. We also compare various prescriptions for the osmotic pressure of suspensions in Donnan equilibrium with a salt ion reservoir and analyze instances where the macroion effective charge becomes larger than the bare one. The methods assessed include single-center cell, jellium, and multi-center mean-field theories. The strengths and weaknesses of the various methods are critically assessed, with the aim of guiding optimal and accurate implementations.

Cation Adsorption in TiO2 Nanotubes: Implication for Water Decontamination

TiO2 nanotubes constitute very promising nanomaterials for water decontamination by the removal of cations. We combined a range of experimental techniques from structural analyses to measurements of the properties of aqueous suspensions of nanotubes, with (i) continuous solvent modeling and (ii) quantum DFT-based simulations to assess the adsorption of Cs+ on TiO2 nanotubes and to predict the separation of metal ions. The methodology is set to be operable under realistic conditions, which, in this case, include the presence of CO2 that needs to be treated as a substantial contaminant, both in experiments and in models. The mesoscopic model, based on the Poisson-Boltzmann equation and surface adsorption equilibrium, predicts that H+ ions are the charge-determining species, while Cs+ ions are in the diffuse layer of the outer surface with a significant contribution only at high concentrations and high pH. The effect of the size of nanotubes in terms of the polydispersity and the distribution of the inner and outer radii is shown to be a third-order effect that is very small when the nanotube layer is not very thick (ranging from 1 to 2 nm). Besides, DFT-based molecular dynamics simulations demonstrate that, for protonation, the one-site and successive association assumption is correct, while, for Cs+ adsorption, the size of the cation is important and the adsorption sites should be carefully defined.

A Gaussian field approach to the planar electric double layer structures in electrolyte solutions

In this work, the planar, electric, double-layer structures of non-polarizable electrodes in electrolyte solutions are studied with Gaussian field theory. A response function with two Yukawa functions is used to capture the electrostatic response of the electrolyte solution, from which the modified response function in the planar symmetry is derived analytically. The modified response function is further used to evaluate the induced charge density and the electrostatic potential near an electrode. The Gaussian field theory, combined with a two-Yukawa response function, can reproduce the oscillatory decay behavior of the electric potentials in concentrated electrolyte solutions. When the exact sum rules for the bulk electrolyte solutions and the electric double layers are used as constraints to determine the parameters of the response function, the Gaussian field theory could at least partly capture the nonlinear response effect of the surface charge density. Comparison with results for a planar electrode with fixed surface charge densities from molecular simulations demonstrates the validity of Gaussian field theory.

Fluids and Electrolytes under Confinement in Single-Digit Nanopores

Confined fluids and electrolyte solutions in nanopores exhibit rich and surprising physics and chemistry that impact the mass transport and energy efficiency in many important natural systems and industrial applications. Existing theories often fail to predict the exotic effects observed in the narrowest of such pores, called single-digit nanopores (SDNs), which have diameters or conduit widths of less than 10 nm, and have only recently become accessible for experimental measurements. What SDNs reveal has been surprising, including a rapidly increasing number of examples such as extraordinarily fast water transport, distorted fluid-phase boundaries, strong ion-correlation and quantum effects, and dielectric anomalies that are not observed in larger pores. Exploiting these effects presents myriad opportunities in both basic and applied research that stand to impact a host of new technologies at the water-energy nexus, from new membranes for precise separations and water purification to new gas permeable materials for water electrolyzers and energy-storage devices. SDNs also present unique opportunities to achieve ultrasensitive and selective chemical sensing at the single-ion and single-molecule limit. In this review article, we summarize the progress on nanofluidics of SDNs, with a focus on the confinement effects that arise in these extremely narrow nanopores. The recent development of precision model systems, transformative experimental tools, and multiscale theories that have played enabling roles in advancing this frontier are reviewed. We also identify new knowledge gaps in our understanding of nanofluidic transport and provide an outlook for the future challenges and opportunities at this rapidly advancing frontier.

The Role of Oxidation Pattern and Water Content in the Spatial Arrangement and Dynamics of Oxidized Graphene-Based Aqueous Dispersions

In this work, we employ fully atomistic molecular dynamics simulations to elucidate the effects of the oxidation pattern and of the water content on the organization of graphene sheets in aqueous dispersions and on the dynamic properties of the different moieties at neutral pH conditions. Analysis of the results reveals the role of the oxidation motif (peripherally or fully oxidized flakes) in the tendency of the flakes to self-assemble and in the control of key structural characteristics, such as the interlayer distance between the sheets and the average size and the distribution of the formed aggregates. In certain cases, the results are compared to a pertinent experimental system, validating further the relevant computational models. Examination of the diffusional motion of the oxidized flakes shows that different degrees of spatial restriction are imposed upon the decrease in the water content and elucidates the conditions under which a motional arrest of the flakes takes place. At constant water content, the structural differences between the formed aggregates appear to additionally impart distinct diffusional characteristics of a water molecule. A detailed examination of the counterion dynamics describes their interaction with the oxidized flakes and their dependence on the water content and on the oxidation pattern, offering new insight into the expected electrical properties of the dispersions. The detailed information provided by this work will be particularly useful in applications such as molecular sieving, nanofiltration, and in cases where conductive membranes based on oxidized forms of graphene are used.

Strongly enhanced dynamics of a charged Rouse dimer by an external magnetic field

While the dynamics of dimers and polymer chains in a viscous solvent is well understood within the celebrated Rouse model, the effect of an external magnetic field on the dynamics of a charged chain is much less understood. Here, we generalize the Rouse model for a charged dimer to include the effect of an external magnetic field. Our analytically solvable model allows a fundamental insight into the magneto-generated dynamics of the dimer in the overdamped limit as induced by the Lorentz force. Surprisingly, for a dimer of oppositely charged particles, we find an enormous enhancement of the dynamics of the dimer center, which exhibits even a transient superballistic behavior. This is highly unusual in an overdamped system for there is neither inertia nor any internal or external driving. We attribute this to a significant translation and rotation coupling due to the Lorentz force. We also find that magnetic field reduces the mobility of a dimer along its orientation and its effective rotational diffusion coefficient. In principle, our predictions can be tested by experiments with colloidal particles and complex plasmas.

Structural correlations in highly asymmetric binary charged colloidal mixtures

We explore structural correlations of strongly asymmetric mixtures of binary charged colloids within the primitive model of electrolytes considering large charge and size ratios of 10 and higher. Using computer simulations with explicit microions, we obtain the partial pair correlation functions between the like-charged colloidal macroions. Interestingly the big-small correlation peak amplitude is smaller than that of the big-big and small-small macroion correlation peaks, which is unfamiliar for additive repulsive interactions. Extracting optimal effective microion-averaged pair interactions between the macroions, we find that on top of non-additive Yukawa-like repulsions an additional shifted Gaussian attractive potential between the small macroions is needed to accurately reproduce their correct pair correlations. For small Coulomb couplings, the behavior is reproduced in a coarse-grained theory with microion-averaged effective interactions between the macroions. However, the accuracy of the theory deteriorates with increasing Coulomb coupling. We emphasize the relevance of entropic interactions exerted by the microions on the macroions. Our results are experimentally verifiable in binary mixtures of micron-sized colloids and like-charge nanoparticles.

Like-Charge Attraction at the Nanoscale: Ground-State Correlations and Water Destructuring

Like-charge attraction, driven by ionic correlations, challenges our understanding of electrostatics both in soft and hard matter. For two charged planar surfaces confining counterions and water, we prove that, even at relatively low correlation strength, the relevant physics is the ground-state one, oblivious of fluctuations. Based on this, we derive a simple and accurate interaction pressure that fulfills known exact requirements and can be used as an effective potential. We test this equation against implicit-solvent Monte Carlo simulations and against explicit-solvent simulations of cement and several types of clays. We argue that water destructuring under nanometric confinement drastically reduces dielectric screening, enhancing ionic correlations. Our equation of state at reduced permittivity therefore explains the exotic attractive regime reported for these materials, even in the absence of multivalent counterions.

Microscopic Simulations of Electrochemical Double-Layer Capacitors

Electrochemical double-layer capacitors (EDLCs) are devices allowing the storage or production of electricity. They function through the adsorption of ions from an electrolyte on high-surface-area electrodes and are characterized by short charging/discharging times and long cycle-life compared to batteries. Microscopic simulations are now widely used to characterize the structural, dynamical, and adsorption properties of these devices, complementing electrochemical experiments and in situ spectroscopic analyses. In this review, we discuss the main families of simulation methods that have been developed and their application to the main family of EDLCs, which include nanoporous carbon electrodes. We focus on the adsorption of organic ions for electricity storage applications as well as aqueous systems in the context of blue energy harvesting and desalination. We finally provide perspectives for further improvement of the predictive power of simulations, in particular for future devices with complex electrode compositions.

The Role of Energy Scales for the Structure of Ionic Liquids at Electrified Interfaces: A Theory-Based Approach

Ionic liquids offer unique bulk and interfacial characteristics as battery electrolytes. Our continuum approach naturally describes the electrolyte on a macroscale. An integral formulation for the molecular repulsion, which can be quantitatively determined by both experimental and theoretical methods, models the electrolyte on the nanoscale. In this article, we perform a systematic series expansion of this integral formulation, derive a description of chemical potentials in terms of higher-order concentration gradients, and rationalize the appearance of fourth-order derivative operators in modified Poisson equations, as recently proposed in this context. In this way, we formulate a rigorous multiscale methodology from atomistic quantum chemistry calculations to phenomenological continuum models. We apply our generalized framework to ionic liquids near electrified interfaces and perform analytical asymptotic analysis. Three energy scales describing electrostatic forces between ions, molecular repulsion, and thermal motion determine the shape and width of the long-ranging charged double layer. We classify the charge screening mechanisms dependent on the system parameters as dielectricity, ion size, interaction strength, and temperature. We find that the charge density of electrochemical double layers in ionic liquids either decays exponentially, for negligible molecular repulsion, or oscillates continuously. Charge ordering across several ion diameters occurs if the repulsion between molecules is comparable with thermal energy and Coulomb interactions. Eventually, phase separation of the bulk electrolyte into ionic layers emerges once the molecular repulsion becomes dominant. Our framework predicts the exact phase boundaries among these three phases as a function of temperature, dielectricity, and ion size.

Self-consistent relaxation theory of collective ion dynamics in Yukawa one-component plasmas under intermediate screening regimes

The self-consistent relaxation theory is employed to describe the collective ion dynamics in strongly coupled Yukawa classical one-component plasmas. The theory is applied to equilibrium states corresponding to intermediate screening regimes with appropriate values of the structure and coupling parameters. The information about the structure (the radial distribution function and the static structure factor) and the thermodynamics of the system are sufficient to describe collective dynamics over a wide range of spatial scales, namely, from the extended hydrodynamic to the microscopic dynamics scale. The main experimentally measurable characteristics of the equilibrium collective dynamics of ions-the spectrum of the dynamic structure factor, the dispersion parameters, the speed of sound, and the sound attenuation-are determined within the framework of the theory without using any adjustable parameters. The results demonstrate agreement with molecular dynamics simulations. Thus a direct realization is presented of the key idea of statistical mechanics: for the theoretical description of the collective particle dynamics in equilibrium fluids it is sufficient to know the interparticle interaction potential and the structural characteristics. Comparison with alternative or complementary theoretical approaches is provided.

Application of nanoarchitectonics in moist-electric generation

The consumption of energy is an important resource that cannot be ignored in modern society. Non-renewable forms of energy, such as coal, natural gas, and oil, have always been important strategic resources and are always facing a crisis of shortage. Therefore, there is an urgent need for green renewable forms of energy. As an emerging green energy source, the moist-electric generator (MEG) has been studied in recent years and may become an energy source that can be utilized in daily life. Along with the advancement of technological means, nanoarchitectonics play an important role in MEG devices. This review aims to provide a comprehensive summary of the fundamentals of the MEG from the perspective of different material classifications and to provide guidance for future work in the field of MEGs. The effects of various parameters and structural designs on the output power, recent important literature and works, the mechanism of liquid-solid interactions at the nanoscale, and the application status and further potential of MEG devices are discussed in this review. It is expected that this review may provide valuable knowledge for future MEG research.

Pimples reduce and dimples enhance flat dielectric surface image repulsion

In solid-liquid, or liquid-liquid, interfaces with dielectric contrast, charged particles interact with the induced polarization charge of the interface. These interactions contribute to an effective self-energy of the bulk ions and mediate ion-ion interactions. For flat interfaces, the self-energy and the mediated interactions are neatly constructed by the image charge method. For other geometries, explicit results are scarce and the problem must be treated via approximations or direct computation. The case of interfaces with roughness is of great practical importance. This article provides analytical results, valid to first-order in perturbation theory, for the self-energy of particles near rough substrates. Explicit formulas are provided for the case of a sinusoidal deformation of a flat surface. Generic deformations can be treated by superposition. In addition to results for the self-energy, the surface polarization charge is presented as a quadrature. The interaction between an ion and the deformed surface is modified by the change in relative distance as well as by the local curvature of the surface. Solid walls, with a lower dielectric constant than the liquid, repel all ions. We show that the repulsion is reduced by local convexity and enhanced by concavity; dimples are more repulsive than pimples.

Numerical Studies on Electrostatic Interaction Forces and the Free Energy between Parallel Colloidal Rods of Finite Size in Skewed Configurations

Formulas for interaction forces F(s) and the free energy G(s) between two parallel charged prismatic rods of various scaled values of d, ψ_s, and L in skewed configurations are obtained, where s is the lengthwise positional difference between the front-end faces of the respective rods, and d is the minimal distance between the opposing faces of the rods, ψ_s is the electric surface potential, L is the length of the rods. To obtain the free-energy function G(s), (i) 3D spatial distributions of the electric potential ψ around two rods were determined by numerically solving the nonlinear Poisson-Boltzmann equation with a finite element method, (ii) with the ψ distributions so determined, the lengthwise interaction electrostatic Maxwell stress tangential to the midplane between the rods was calculated to obtain the (discrete) s dependence of the stress, and (iii) by introducing two different fitting functions, the discrete s dependence was transformed into a continuous force function, F(s), which was then lengthwise integrated to derive G(s). It was found that the curves of G(s) linearly decreased with increasing s between 1 and L + 1 due to a localization of the stress. Although natural, it is of interest that the values of G(0) calculated for rods of various values of d, ψ_s, and L were in good agreement with those of the interaction free energy obtained in our preceding work by the widthwise integration of repulsive electrostatic forces normal to the midplane between the parallel rods in nonskewed configurations.

Tunability of Interactions between the Core and Shell in Rattle-Type Particles Studied with Liquid-Cell Electron Microscopy

Yolk-shell or rattle-type particles consist of a core particle that is free to move inside a thin shell. A stable core with a fully accessible surface is of interest in fields such as catalysis and sensing. However, the stability of a charged nanoparticle core within the cavity of a charged thin shell remains largely unexplored. Liquid-cell (scanning) transmission electron microscopy is an ideal technique to probe the core-shell interactions at nanometer spatial resolution. Here, we show by means of calculations and experiments that these interactions are highly tunable. We found that in dilute solutions adding a monovalent salt led to stronger confinement of the core to the middle of the geometry. In deionized water, the Debye length κ-1 becomes comparable to the shell radius Rshell, leading to a less steep electric potential gradient and a reduced core-shell interaction, which can be detrimental to the stability of nanorattles. For a salt concentration range of 0.5-250 mM, the repulsion was relatively long-ranged due to the concave geometry of the shell. At salt concentrations of 100 and 250 mM, the core was found to move almost exclusively near the shell wall, which can be due to hydrodynamics, a secondary minimum in the interaction potential, or a combination of both. The possibility of imaging nanoparticles inside shells at high spatial resolution with liquid-cell electron microscopy makes rattle particles a powerful experimental model system to learn about nanoparticle interactions. Additionally, our results highlight the possibilities for manipulating the interactions between core and shell that could be used in future applications.

A temperature-dependent critical Casimir patchy particle model benchmarked onto experiment

Synthetic colloidal patchy particles immersed in a binary liquid mixture can self-assemble via critical Casimir interactions into various superstructures, such as chains and networks. Up to now, there are no quantitatively accurate potential models that can simulate and predict this experimentally observed behavior precisely. Here, we develop a protocol to establish such a model based on a combination of theoretical Casimir potentials and angular switching functions. Using Monte Carlo simulations, we optimize several material-specific parameters in the model to match the experimental chain length distribution and persistence length. Our approach gives a systematic way to obtain accurate potentials for critical Casimir induced patchy particle interactions and can be used in large-scale simulations.

Phase Transitions and Electrochemical Properties of Ionic Liquids and Ionic Liquid-Solvent Mixtures

Recent advances in studies of ionic liquids (IL) and ionic liquid-solvent mixtures are reviewed. Selected experimental, simulation, and theoretical results for electrochemical, thermodynamical, and structural properties of IL and IL-solvent mixtures are described. Special attention is paid to phenomena that are not predicted by the classical theories of the electrical double layer or disagree strongly with these theories. We focus on structural properties, especially on distribution of ions near electrodes, on electrical double layer capacitance, on effects of confinement, including decay length of a dissjoining pressure between confinig plates, and on demixing phase transition. In particular, effects of the demixing phase transition on electrochemical properties of ionic liquid-solvent mixtures for different degrees of confinement are presented.

Smoothed profile method for direct numerical simulations of hydrodynamically interacting particles

A general method is presented for computing the motions of hydrodynamically interacting particles in various kinds of host fluids for arbitrary Reynolds numbers. The method follows the standard procedure for performing direct numerical simulations (DNS) of particulate systems, where the Navier-Stokes equation must be solved consistently with the motion of the rigid particles, which defines the temporal boundary conditions to be satisfied by the Navier-Stokes equation. The smoothed profile (SP) method provides an efficient numerical scheme for coupling the continuum fluid mechanics with the dispersed moving particles, which are allowed to have arbitrary shapes. In this method, the sharp boundaries between solid particles and the host fluid are replaced with a smeared out thin shell (interfacial) region, which can be accurately resolved on a fixed Cartesian grid utilizing a SP function with a finite thickness. The accuracy of the SP method is illustrated by comparison with known exact results. In the present paper, the high degree of versatility of the SP method is demonstrated by considering several types of active and passive particle suspensions.

Bayesian unsupervised learning reveals hidden structure in concentrated electrolytes

Electrolytes play an important role in a plethora of applications ranging from energy storage to biomaterials. Notwithstanding this, the structure of concentrated electrolytes remains enigmatic. Many theoretical approaches attempt to model the concentrated electrolyte by introducing the idea of ion pairs, with ions either being tightly "paired" with a counter-ion or "free" to screen charge. In this study, we reframe the problem into the language of computational statistics and test the null hypothesis that all ions share the same local environment. Applying the framework to molecular dynamics simulations, we find that this null hypothesis is not supported by data. Our statistical technique suggests the presence of two distinct local ionic environments at intermediate concentrations, whose differences surprisingly originate in like charge correlations rather than unlike charge attraction. Through considering the effect of these "aggregated" and "non-aggregated" states on bulk properties including effective ion concentration and dielectric constant, we identify a scaling relation between the effective screening length and theoretical Debye length, which applies across different dielectric constants and ion concentrations.

Electric-field-based Poisson-Boltzmann theory: Treating mobile charge as polarization

Mobile charge in an electrolytic solution can in principle be represented as the divergence of ionic polarization. After adding explicit solvent polarization a finite volume of an electrolyte can then be treated as a composite nonuniform dielectric body. Writing the electrostatic interactions as an integral over electric-field energy density we show that the Poisson-Boltzmann functional in this formulation is convex and can be used to derive the equilibrium equations for electric potential and ion concentration by a variational procedure developed by Ericksen for dielectric continua [J. L. Ericksen, Arch. Rational Mech. Anal. 183, 299 (2007)AVRMAW0003-952710.1007/s00205-006-0042-4]. The Maxwell field equations are enforced by extending the set of variational parameters by a vector potential representing the dielectric displacement which is fully transverse in a dielectric system without embedded external charge. The electric-field energy density in this representation is a function of the vector potential and the sum of ionic and solvent polarization making the mutual screening explicit. Transverse polarization is accounted for by construction, lifting the restriction to longitudinal polarization inherent in the electrostatic potential based formulation of Poisson-Boltzmann mean field theory.

Reversible heat production during electric double layer buildup depends sensitively on the electrolyte and its reservoir

Several modern technologies for energy storage and conversion are based on the screening of electric charge on the surface of porous electrodes by ions in an adjacent electrolyte. This so-called electric double layer (EDL) exhibits an intricate interplay with the electrolyte's temperature that was the focus of several recent studies. In one of them, Janssen et al. [Phys. Rev. Lett. 119, 166002 (2017)] experimentally determined the ratio Q_rev/W_el of reversible heat flowing into a supercapacitor during an isothermal charging process and the electric work applied therein. To rationalize that data, here, we determine Q_rev/W_el within different models of the EDL using theoretical approaches such as density functional theory (DFT) as well as molecular dynamics simulations. Applying mainly the restricted primitive model, we find quantitative support for a speculation of Janssen et al. that steric ion interactions are key to the ratio Q_rev/W_el. Here, we identified the entropic contribution of certain DFT functionals, which grants direct access to the reversible heat. We further demonstrate how Q_rev/W_el changes when calculated in different thermodynamic ensembles and processes. We show that the experiments of Janssen et al. are explained best by a charging process at fixed bulk density or in a "semi-canonical" system. Finally, we find that Q_rev/W_el significantly depends on parameters such as pore and ion size, salt concentration, and valencies of the cations and anions of the electrolyte. Our findings can guide further heat production measurements and can be applied in studies on, for instance, nervous conduction, where reversible heat is a key element.

Size and charge correlations in spherical electric double layers: a case study with fully asymmetric mixed electrolytes within the solvent primitive model

Size and charge correlations in spherical electric double layers are investigated through Monte Carlo simulations and density functional theory, through a solvent primitive model representation. A fully asymmetric mixed electrolyte is used for the small ions, whereas the solvent, apart from being a continuum dielectric, is also treated as an individual component. A partially perturbative density functional theory is adopted here, and for comparison, a standard canonical ensemble Monte Carlo simulation is used. The hard-sphere free energy is treated within a weighted density approach and the residual ionic contribution is estimated through perturbation around the uniform density. The results from both methods corroborate each other quantitatively over a wide range of physical parameters. The importance of structural correlations is envisaged through the size and charge asymmetry of the supporting electrolytes that includes the solvent as a component.

Size and charge correlations in spherical electric double layers are investigated through Monte Carlo simulations and density functional theory, through a solvent primitive model representation.

Molecular dynamics simulations of the colloidal interaction between smectite clay nanoparticles in liquid water

Colloidal interactions between clay nanoparticles have been studied extensively because of their strong influence on the hydrology and mechanics of many soils and sedimentary media. The predominant theory used to describe these interactions is the Derjaguin-Landau-Verwey-Overbeek (DLVO) model, a framework widely applied in colloidal and interfacial science that accurately predicts the interactions between charged surfaces across water films at distances greater than ~ 3 nm (i.e., ten water monolayers). Unfortunately, the DLVO model is inaccurate at the shorter interparticle distances that predominate in most subsurface environments. For example, it inherently cannot predict the existence of equilibrium states wherein clay particles adopt interparticle distances equal to the thickness of one, two, or three water monolayers. Molecular dynamics (MD) simulations have the potential to provide detailed information on the free energy of interaction between clay nanoparticles; however, they have only been used to examine clay swelling and aggregation at interparticle distances below 1 nm. We present the first MD simulation predictions of the free energy of interaction of smectite clay nanoparticles in the entire range of interparticle distances from the large interparticle distances where the DLVO model is accurate (>3 nm) to the short-range swelling states where non-DLVO interactions predominate (<1 nm). Our simulations examine a range of salinities (0.0 to 1.0 M NaCl) and counterion types (Na, K, Ca) and establish a detailed picture of the breakdown of the DLVO model. In particular, they confirm previous theoretical suggestions of the existence of a strong non-DLVO attraction with a range of ~ 3 nm arising from specific ion-clay Coulomb interactions in the electrical double layer.

Strong-coupling theory of counterions with hard cores between symmetrically charged walls

By a combination of Monte Carlo simulations and analytical calculations, we investigate the effective interactions between highly charged planar interfaces, neutralized by mobile counterions (salt-free system). While most previous analysis have focused on pointlike counterions, we treat them as charged hard spheres. We thus work out the fate of like-charge attraction when steric effects are at work. The analytical approach partitions counterions in two subpopulations, one for each plate, and integrates out one subpopulation to derive an effective Hamiltonian for the remaining one. The effective Hamiltonian features plaquette four-particle interactions, and it is worked out by computing a Gibbs-Bogoliubov inequality for the free energy. At the root of the treatment is the fact that under strong electrostatic coupling, the system of charges forms an ordered arrangement, that can be affected by steric interactions. Fluctuations around the reference positions are accounted for. To dominant order at high coupling, it is found that steric effects do not significantly affect the interplate effective pressure, apart at small distances where hard-sphere overlap are unavoidable, and thus rule out configurations.

Disappearance of the polyelectrolyte peak in salt-free solutions

We investigate the nature of the polyelectrolyte peak in salt-free solutions by molecular dynamics simulations using a minimal model of polyelectrolyte solutions that includes an explicit solvent and counterions and small angle scattering experiments. It is found that the polyelectrolyte peak progressively disappears as the strength of solvation for the charged species is increased and the scattering profiles start to resemble those of neutral polymer solutions. The disappearance of the polyelectrolyte peak coincides with the emergence of attractive interchain interactions over a wide range of length scales. These findings provide insights into the microscopic origin of the polyelectrolyte peak.

Modeling hydration-mediated ion-ion interactions in electrolytes through oscillating Yukawa potentials

Classical Poisson-Boltzmann theory represents a mean-field description of the electric double layer in the presence of only Coulomb interactions. However, aqueous solvents hydrate ions, which gives rise to additional hydration-mediated ion-ion interactions. Experimental and computational studies suggest damped oscillations to be a characteristic feature of these hydration-mediated interactions. We have therefore incorporated oscillating Yukawa potentials into the mean-field description of the electric double layer. This is accomplished by allowing the decay length of the Yukawa potential to be complex valued. Ion specificity emerges from assigning individual strengths and phases to the Yukawa potential for anion-anion, anion-cation, and cation-cation pairs as well as for anions and cations interacting with an electrode or macroion. Excluded volume interactions between ions are approximated by replacing the ideal gas entropy by that of a lattice gas. We derive mean-field equations for the Coulomb and Yukawa potentials and use their solutions to compute the differential capacitance for an isolated planar electrode and the pressure that acts between two planar, like-charged macroion surfaces. Attractive interactions appear if the surface charge density of the macroions is sufficiently small.

Numerical Studies on Electrical Interaction Forces and Free Energy between Colloidal Plates of Finite Size

By solving the nonlinear Poisson-Boltzmann (PB) equation with a finite element method (FEM), three-dimensional (3D) spatial distributions of the electric potential (ψ, scaled) in electrolyte solutions having two charged parallel finite plates (including cubes and prismatic rods) are determined for various separations (d, scaled by the Debye length, κ^-1), surface potentials (ψ_s), and plate dimensions (length × width × thickness, each scaled by κ^-1). The total interaction force between two plates, F, is the sum of the electrostatic double-layer (EDL) repulsion (the osmotic pressure, F_osm) and the Maxwell electrostatic stress (F_es). The EDL repulsion is estimated using the distribution of ψ not only between the facing surfaces of two parallel plates but also around the other extremities of the plates. The Maxwell stress (F_es) is localized near the extremities to act as a repulsive force on the midplane between the two plates. The ratio F_es/F is 0.07-0.5, depending on d, ψ_s, and dimensions. It is found that, with increasing dimensions, the total F values per unit area calculated for finite plates, F̃, decreasingly approach the exact ones for parallel infinite plates, F̃^inf; for example, at d = 1 and ψ_s = 5, the ratio F̃/F̃^inf is 2.83 for plates with dimensions of 1 × 1 × 1 and 1.18 for plates of 10 × 10 × 1. The repulsions arising from the extremities cannot be neglected for plates with dimensions <10 × 10 × 1. Furthermore, the total interaction forces (F) are calculated at a series of discrete d values, respectively, for parallel plates. We introduce a force fitting function, F_f(d), with parameters that can be determined so that F_f(d) fits well to the calculated serial F values. By integrating the F_f(d), we obtain the interaction free energy, G(d), for finite parallel plates that consists of two Γ functions.

Polyelectrolyte character of rigid rod peptide bundlemer chains constructed via hierarchical self-assembly

Short α-helical peptides were computationally designed to self-assemble into robust coiled coils that are antiparallel, homotetrameric bundles. These peptide bundle units, or 'bundlemers', have been utilized as anisotropic building blocks to construct bundlemer-based polymers via a hierarchical, hybrid physical-covalent assembly pathway. The bundlemer chains were constructed using short linker connections via 'click' chemistry reactions between the N-termini of bundlemer constituent peptides. The resulting bundlemer chains appear as extremely rigid, cylindrical rods in transmission electron microscopy (TEM) images. Small angle neutron scattering (SANS) shows that these bundlemer chains exist as individual rods in solution with a cross-section that is equal to that of a single coiled coil bundlemer building block of ≈20 Å. SANS further confirms that the interparticle solution structure of the rigid rod bundlemer chains is heterogeneous and responsive to solution conditions, such as ionic-strength and pH. Due to their peptidic constitution, the bundlemer assemblies behave like polyelectrolytes that carry an average charge density of approximately 3 charges per bundlemer as determined from SANS structure factor data fitting, which describes the repulsion between charged rods in solution. This repulsion manifests as a correlation hole in the scattering profile that is suppressed by dilution or addition of salt. Presence of rod cluster aggregates with a mass fractal dimension of ≈2.5 is also confirmed across all samples. The formation of such dense, fractal-like cluster aggregates in a solution of net repulsive rods is a unique example of the subtle balance between short-range attraction and long-rage repulsion interactions in proteins and other biomaterials. With computational control of constituent peptide sequences, it is further possible to deconvolute the underlying sequence driven structure-property relationships in the modular bundlemer chains.

Electric double layers with surface charge modulations: Exact Poisson-Boltzmann solutions

Poisson-Boltzmann theory is the cornerstone for soft matter electrostatics. We provide exact analytical solutions to this nonlinear mean-field approach for the diffuse layer of ions in the vicinity of a planar or a cylindrical macroion. While previously known solutions are for homogeneously charged objects, the cases worked out exhibit a modulated surface charge-or equivalently, surface potential-on the macroion (wall) surface. In addition to asymptotic features at large distances from the wall, attention is paid to the fate of the contact theorem, relating the contact density of ions to the local wall charge density. For salt-free systems (counterions only), we make use of results pertaining to the two-dimensional Liouville equation, supplemented by an inverse approach. When salt is present, we invoke the exact two-soliton solution to the 2D sinh-Gordon equation. This leads to inhomogeneous charge patterns, that are either localized or periodic in space. Without salt, the electrostatic signature of a charge pattern on the macroion fades exponentially with distance for a planar macroion, while it decays as an inverse power law for a cylindrical macroion. With salt, our study is limited to the planar geometry and reveals that pattern screening is exponential.

Long-range forces and charge inversions in model charged colloidal dispersions at finite concentration

In charged colloidal dispersion systems the interest is in finding their stability conditions, phase transitions, and transport properties, either in bulk or confinement, among other physicochemical quantities, for which the knowledge of the dispersions' molecular structure and the associated macroion-macroion forces is crucial. To investigate these phenomena simple models have been proposed. Most of the theoretical and simulation studies on charged particles suspensions are at infinite dilution conditions. Hence, these studies have been focused on the electrolyte structure around one or two isolated central particle(s), where phenomena as charge reversal, charge inversion and surface charge amplification have been shown to be relevant. However, experimental studies at finite volume fraction exhibit interesting phenomenology which imply very long-range correlations. A simple, yet useful, model is the Colloidal Primitive Model, in which the colloidal dispersion is modeled as a mixture of size (and charge) asymmetrical hard spheres, at finite volume fraction. In this paper we review recent integral equations solutions for this model, where very long-range attractive-repulsive forces, as well as new long-range, giant charge inversions are reported. The calculated macroions radial distribution functions, charge distributions, and macroion-macroion forces are qualitatively consistent with existing experimental results, and Monte Carlo and molecular dynamics simulations.

Onset of transverse (shear) waves in strongly-coupled Yukawa fluids

A simple practical approach to describe transverse (shear) waves in strongly-coupled Yukawa fluids is presented. Theoretical dispersion curves, based on hydrodynamic consideration, are shown to compare favorably with existing numerical results for plasma-related systems in the long-wavelength regime. The existence of a minimum wave number below which shear waves cannot propagate and its magnitude are properly accounted in the approach. The relevance of the approach beyond plasma-related Yukawa fluids is demonstrated by using experimental data on transverse excitations in liquid metals Fe, Cu, and Zn, obtained from inelastic x-ray scattering. Some potentially important relations, scalings, and quasi-universalities are discussed. The results should be interesting for a broad community in chemical physics, materials physics, physics of fluids and glassy state, complex (dusty) plasmas, and soft matter.

Asymmetric electrolytes near structured dielectric interfaces

The ion distribution of electrolytes near interfaces with dielectric contrast has important consequences for electrochemical processes and many other applications. To date, most studies of such systems have focused on geometrically simple interfaces, for which dielectric effects are analytically solvable or computationally tractable. However, all real surfaces display nontrivial structure at the nanoscale and have, in particular, a nonuniform local curvature. Using a recently developed, highly efficient computational method, we investigate the effect of surface geometry on ion distribution and interface polarization. We consider an asymmetric 2:1 electrolyte bounded by a sinusoidally deformed solid surface. We demonstrate that even when the surface is neutral, the electrolyte acquires a nonuniform ion density profile near the surface. This profile is asymmetric and leads to an effective charging of the surface. We furthermore show that the induced charge is modulated by the local curvature. The effective charge is opposite in sign to the multivalent ions and is larger in concave regions of the surface.

Electrostatics of Nanoparticle-Wall Interactions within Nanochannels: Role of Double-Layer Structure and Ion-Ion Correlations

We perform computational investigations of the electrolyte-mediated interactions of charged nanoparticles with the walls of nanochannels. We investigate the role of discrete ion effects, valence, and electrolyte strength on nanoparticle-wall interactions. We find for some of the multivalent charge regimes that the like-charged nanoparticles and walls can have attractive interactions. We study in detail these interactions and the free-energy profile for the nanoparticle-wall separation. We find there are energy barriers and energy minima giving preferred nanoparticle locations in the channel near the center and at a distance near to but separated from the channel walls. We characterize contributions from surface overcharging, condensed layers, and overlap of ion double layers. We perform our investigations using coarse-grained particle-level simulations with Brownian dynamics, classical density functional theory, and the mean-field Poisson-Boltzmann theory. We discuss the implications of our results for phenomena in nanoscale devices.

Visualizing charges accumulated in an electric double layer by three-dimensional open-loop electric potential microscopy

Charges accumulated in an electric double layer (EDL) play key roles in various interfacial phenomena and electronic devices. However, direct imaging of their spatial distribution has been a great challenge, which has hindered our nano-level understanding of the mechanisms of such interfacial phenomena and functions. In this study, we present direct imaging of charges accumulated at an electrode-electrolyte interface using three-dimensional open-loop electric potential microscopy (3D-OL-EPM). Conventional OL-EPM allows us to visualize two-dimensional potential distributions in liquid yet the zero of the measured potential is not well defined due to the influence of the long-range (LR) interaction between the cantilever and the sample. Here, we present practical ways to reduce such an influence by improving the equation for the potential calculation and subtracting the LR contribution estimated from a Z potential profile. These improvements enabled the calibration of the measured potential values with respect to the bulk solution potential. With these improvements, we visualized opposite charge accumulation behaviors on a polarizable and non-polarizable electrode with a varying electrode potential. Combining OL-EPM with a 3D tip scanning method, we also performed a 3D-OL-EPM measurement on a Cu fine wire and visualized the nanoscale distribution of the charges accumulated at the interface. Such real-space information on the charge distributions in an EDL should provide valuable insights into the mechanisms of interfacial phenomena and functions that are important in various academic and industrial research on electronic devices, electrochemistry, tribology and life sciences.

Ionic size effects on the Poisson-Boltzmann theory

In this paper, we develop a simple theory to study the effects of ionic size on ionic distributions around a charged spherical particle. We include a correction to the regular Poisson-Boltzmann equation in order to take into account the size of ions in a mean-field regime. The results are compared with Monte Carlo simulations and a density functional theory based on the fundamental measure approach and a second-order bulk expansion which accounts for electrostatic correlations. The agreement is very good even for multivalent ions. Our results show that the theory can be applied with very good accuracy in the description of ions with highly effective ionic radii and low concentration, interacting with a colloid or a nanoparticle in an electrolyte solution.

Clustering and assembly dynamics of a one-dimensional microphase former

Both ordered and disordered microphases ubiquitously form in suspensions of particles that interact through competing short-range attraction and long-range repulsion (SALR). While ordered microphases are more appealing materials targets, understanding the rich structural and dynamical properties of their disordered counterparts is essential to controlling their mesoscale assembly. Here, we study the disordered regime of a one-dimensional (1D) SALR model, whose simplicity enables detailed analysis by transfer matrices and Monte Carlo simulations. We first characterize the signature of the clustering process on macroscopic observables, and then assess the equilibration dynamics of various simulation algorithms. We notably find that cluster moves markedly accelerate the mixing time, but that event chains are of limited help in the clustering regime. These insights will inspire further study of three-dimensional microphase formers.