Mean-field Poisson-Boltzmann Theory Research Articles

Interrelationship of Electric Double Layer Theory and Microfluidic Microbial Fuel Cells: A Review of Theoretical Foundations and Implications for Performance

Microbial fuel cells and their related microfluidic systems have emerged as promising greener energy alternatives for the exploitation of avenues related to combined power and wastewater treatment operations. Moreover, the potential for their application in biosensing technology is large. However, while the fundamental principles of science that govern the design and operation of microbial fuel cells (MFCs) and microfluidic microbial fuel cells (MMFCs) are similar to those found in colloid science, the literature shows that current research lacks sufficient reference to the electrostatic and electrokinetic aspects, focusing mostly on aspects related to the architecture, design, anodes, microbial growth and metabolism, and electron transfer mechanisms. In this regard, research is yet to consider MFCs and MMFCs in the context of electrostatic and electrokinetic aspects. In this extensive review, we show, for the first time, the interrelationship of MFCs and MMFCs with electric double layer theory. Consequently, we show how the analytical solution to the mean field Poisson–Boltzmann theory relates to these systems. Moreover, we show the interrelationship between MFC and MMFCs’ performance and the electric double layer and the associated electrostatic and electrokinetic phenomena. This extensive review will likely motivate research in this direction.

Lipid shape and packing are key for optimal design of pH-sensitive mRNA lipid nanoparticles

The ionizable-lipid component of RNA-containing nanoparticles controls the pH-dependent behavior necessary for an efficient delivery of the cargo-the so-called endosomal escape. However, it is still an empirical exercise to identify optimally performing lipids. Here, we study two well-known ionizable lipids, DLin-MC3-DMA and DLin-DMA using a combination of experiments, multiscale computer simulations, and electrostatic theory. All-atom molecular dynamics simulations, and experimentally measured polar headgroup pKa values, are used to develop a coarse-grained representation of the lipids, which enables the investigation of the pH-dependent behavior of lipid nanoparticles (LNPs) through Monte Carlo simulations, in the absence and presence of RNA molecules. Our results show that the charge state of the lipids is determined by the interplay between lipid shape and headgroup chemistry, providing an explanation for the similar pH-dependent ionization state observed for lipids with headgroup pKa values about one-pH-unit apart. The pH dependence of lipid ionization is significantly influenced by the presence of RNA, whereby charge neutrality is achieved by imparting a finite and constant charge per lipid at intermediate pH values. The simulation results are experimentally supported by measurements of α-carbon 13C-NMR chemical shifts for eGFP mRNA LNPs of both DLin-MC3-DMA and DLin-DMA at various pH conditions. Further, we evaluate the applicability of a mean-field Poisson-Boltzmann theory to capture these phenomena.

Conditional existence of Donnan potential in soft particles and surfaces: Dependence on steric effects mediated by electrolyte ions and structural charges

When a charged layer decorating a particle or a macroscopic surface is equilibrated with an electrolyte solution, a constant Donnan potential is established through that layer due to charge-driven accumulation of counterions and companion exclusion of coions. This situation arises when the thickness of the surface layer well exceeds the screening Debye length, a condition derived from mean-field Poisson-Boltzmann theory within point-like charge approximation. Herein, we revisit this condition underlying the applicability of Donnan electrostatic representation with the account of steric effects mediated by the sizes of the electrolyte ions and structural layer charges. A transcendental equation is derived for the Donnan potential as a function of sizes and valences of anions and cations, electrolyte concentration and size of the layer charges, and a closed-form expression is provided for symmetrical electrolytes. Therefrom we evidence that the existence of a Donnan potential is conditioned not only to large values of the layer thickness compared to a here-defined Debye length operative within the shell, but to additional verification of a criterion that involves space charge density of the layer, solution ionic strength and electrolyte nondiluteness parameter. Illustrative computational examples show how the existence and magnitude of the Donnan potential depend on the key molecular descriptors of the electrolyte and soft interface, and they further quantify the deviations from predictions based on classical Donnan potential expression valid for dilute electrolytes.

Electrostatics of soft (bio)interfaces: Corrections of mean-field Poisson-Boltzmann theory for ion size, dielectric decrement and ion-ion correlations

HypothesisElectrostatics of soft (ion-permeable) (bio)particles (e.g. microorganisms, core/shell colloids) in aqueous electrolytes is commonly formulated by the mean-field Poisson-Boltzmann theory and integration of the charge contributions from electrolyte ions and soft material. However, the effects connected to the size of the electrolyte ions and that of the structural charges carried by the particle, to dielectric decrement and ion-ion correlations on soft interface electrostatics have been so far considered at the margin, despite the limits of the Gouy theory for condensed and/or multivalent electrolytes. ExperimentsAccordingly, we modify herein the Poisson-Boltzmann theory for core/shell (bio)interfaces to include the aforementioned molecular effects considered separately or concomitantly. The formalism is applicable for poorly to highly charged particles in the thin electric double layer regime and to unsymmetrical multivalent electrolytes. FindingsComputational examples of practical interests are discussed with emphasis on how each considered molecular effect or combination thereof affects the interfacial potential distribution depending on size and valence of cations and anions, size of particle charges, length scale of ionic correlations and shell-to-Debye layer thickness ratio. The origins of here-evidenced pseudo-harmonic potential profile and ion size-dependent screening of core/shell particle charges are detailed. In addition, the existence and magnitude of the Donnan potential when reached in the shell layer are shown to depend on the excluded volumes of the electrolyte ions.

Complex electric double layers in charged topological colloids

Charged surfaces in contact with liquids containing ions are accompanied in equilibrium by an electric double layer consisting of a layer of electric charge on the surface that is screened by a diffuse ion cloud in the bulk fluid. This screening cloud determines not only the interactions between charged colloidal particles or polyelectrolytes and their self-assembly into ordered structures, but it is also pivotal in understanding energy storage devices, such as electrochemical cells and supercapacitors. However, little is known to what spatial complexity the electric double layers can be designed. Here, we show that electric double layers of non-trivial topology and geometry -including tori, multi-tori and knots- can be realised in charged topological colloidal particles, using numerical modelling within a mean-field Poisson-Boltzmann theory. We show that the complexity of double layers -including geometry and topology- can be tuned by changing the Debye screening length of the medium, or by changing the shape and topology of the (colloidal) particle. More generally, this work is an attempt to introduce concepts of topology in the field of charged colloids, which could lead to novel exciting material design paradigms.

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.

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.

Understanding Depletion Induced Like-Charge Attraction from Self-Consistent Field Model

AbstractThe interaction force between likely charged particles/surfaces is usually repulsive due to the Coulomb interaction. However, the counterintuitive like-charge attraction in electrolytes has been frequently observed in experiments, which has been theoretically debated for a long time. It is widely known that the mean field Poisson-Boltzmann theory cannot explain and predict this anomalous feature since it ignores many-body properties. In this paper, we develop efficient algorithm and perform the force calculation between two interfaces using a set of self-consistent equations which properly takes into account the electrostatic correlation and the dielectric-boundary effects. By solving the equations and calculating the pressure with the Debye-charging process, we show that the self-consistent equations could be used to study the attraction between like-charge surfaces from weak-coupling to mediate-coupling regimes, and that the attraction is due to the electrostatics-driven entropic force which is significantly enhanced by the dielectric depletion of mobile ions. A systematic investigation shows that the interaction forces can be tuned by material permittivity, ionic size and valence, and salt concentration, and that the like-charge attraction exists only for specific regime of these parameters.

Interplay between Depletion and Double-Layer Forces Acting between Charged Particles in Solutions of Like-Charged Polyelectrolytes.

Direct force measurements between negatively charged silica particles in the presence of a like-charged strong polyelectrolyte were carried out with an atomic force microscope. The force profiles can be quantitatively interpreted as a superposition of depletion and double-layer forces. The depletion forces are modeled with a damped oscillatory profile, while the double-layer forces with the mean-field Poisson-Boltzmann theory for a strongly asymmetric electrolyte, whereby an effective valence must be assigned to the polyelectrolyte. This effective valence is substantially smaller than the bare valence due to ion condensation effects. The unusual aspect of the electrical double layer in these systems is the exclusion of the like-charged polyelectrolyte from the vicinity of the surface, leading to a strongly nonexponential diffuse ionic layer that is dominated by counterions and has a well-defined thickness. As the oscillatory depletion force sets in right after this layer, this condition can be used to predict the phase of the oscillatory depletion force.

Determining the Locations of Ions and Water around DNA from X-Ray Scattering Measurements

Nucleic acids carry a negative charge, attracting salt ions and water. Interactions with these components of the solvent drive DNA to condense, RNA to fold, and proteins to bind. To understand these biological processes, knowledge of solvent structure around the nucleic acids is critical. Yet, because they are often disordered, ions and water evade detection by x-ray crystallography and other high-resolution methods. Small-angle x-ray scattering (SAXS) is uniquely sensitive to the spatial correlations between solutes and the surrounding solvent. Thus, SAXS provides an experimental constraint to guide or test emerging solvation theories. However, the interpretation of SAXS profiles is nontrivial because of the difficulty in separating the scattering signals of each component: the macromolecule, ions, and hydration water. Here, we demonstrate methods for robustly deconvoluting these signals, facilitating a more straightforward comparison with theory. Using SAXS data collected on an absolute intensity scale for short DNA duplexes in solution with Na+, K+, Rb+, or Cs+ counterions, we mathematically decompose the scattering profiles into components (DNA, water, and ions) and validate the decomposition using anomalous scattering measurements. In addition, we generate a library of physically motivated ion atmosphere models and rank them by agreement with the scattering data. The best-fit models have relatively compact ion atmospheres when compared to predictions from the mean-field Poisson-Boltzmann theory of electrostatics. Thus, the x-ray scattering methods presented here provide a valuable measurement of the global structure of the ion atmosphere that can be used to test electrostatics theories that go beyond the mean-field approximation.

Accurate predictions of forces in the presence of multivalent ions by Poisson-Boltzmann theory.

Forces between positively and negatively charged colloidal particles across aqueous salt solutions containing multivalent ions are measured directly with the atomic force microscope (AFM). The measurements are interpreted quantitatively with Poisson-Boltzmann (PB) theory. Thereby, the surface potentials and regulation properties of the particle surfaces are extracted from symmetric measurements between the same types of particles. This information is used to predict force profiles in the asymmetric situations involving different types of particles without any adjustable parameters. These predictions turn out to be very accurate, which demonstrates that the mean-field PB theory is reliable down to distances of about 5 nm. While various reports in the literature indicate that this theory should fail due to neglect of ion correlations, such effects seem important only at higher concentrations and smaller distances.

Polarization as a field variable from molecular dynamics simulations

A theoretical and computational framework for systematically calculating the macroscopic polarization density as a field variable from molecular dynamics simulations is presented. This is done by extending the celebrated Irving and Kirkwood [J. Chem. Phys. 18, 817 (1950)] procedure, which expresses macroscopic stresses and heat fluxes in terms of the atomic variables, to the case of electrostatics. The resultant macroscopic polarization density contains molecular dipole, quadrupole, and higher-order moments, and can be calculated to a desired accuracy depending on the degree of the coarse-graining function used to connect the molecular and continuum scales. The theoretical and computational framework is verified by recovering the dielectric constant of bulk water. Finally, the theory is applied to calculate the spatial variation of the polarization vector in the electrical double layer of a 1:1 electrolyte solution. Here, an intermediate asymptotic length scale is revealed in a specific region, which validates the application of mean field Poisson-Boltzmann theory to describe this region. Also, using the existence of this asymptotic length scale, the lengths of the diffuse and condensed/Stern layers are identified accurately, demonstrating that this framework may be used to characterize electrical double layers over a wide range of concentrations of solutions and surface charges.

Effective interaction between charged nanoparticles and DNA

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

Langevin Poisson-Boltzmann equation: point-like ions and water dipoles near a charged surface

Water ordering near a charged membrane surface is important for many biological processes such as binding of ligands to a membrane or transport of ions across it. In this work, the mean-field Poisson-Boltzmann theory for point-like ions, describing an electrolyte solution in contact with a planar charged surface, is modified by including the orientational ordering of water. Water molecules are considered as Langevin dipoles, while the number density of water is assumed to be constant everywhere in the electrolyte solution. It is shown that the dielectric permittivity of an electrolyte close to a charged surface is decreased due to the increased orientational ordering of water dipoles. The dielectric permittivity close to the charged surface is additionally decreased due to the finite size of ions and dipoles.

Nonextensive statistical mechanics of ionic solutions

Classical mean-field Poisson–Boltzmann theory of ionic solutions is revisited in the theoretical framework of nonextensive Tsallis statistics. The nonextensive equivalent of Poisson–Boltzmann equation is formulated revisiting the statistical mechanics of liquids and the Debye–Hückel framework is shown to be valid for highly diluted solutions even under circumstances where nonextensive thermostatistics must be applied. The lowest order corrections associated to nonadditive effects are identified for both symmetric and asymmetric electrolytes and the behavior of the average electrostatic potential in a homogeneous system is analytically and numerically analyzed for various values of the complexity measurement nonextensive parameter q .

Charging mechanism in colloidal particles leading to a linear relation between charge and size

Colloidal particles that are dissolved in a liquid can become charged by preferential adsorption or dissociation of surface molecules. The assumption of a constant surface charge density leads to a particle charge that is proportional to the square of the radius. Recent experiments however have shown that in a number of cases the charge is linearly proportional to its size. We present a simple model for the charging of colloidal particles leading to a linear relation between charge and size in the framework of the mean-field Poisson-Boltzmann theory. The model is based on preferential adsorption or dissociation of surface molecules, in combination with a large concentration of chargeable surface sites. The model also predicts a surface potential that is independent of the particle size. Numerical simulations show that this model is in good agreement with analytical approximations.

Ion-Mediated Nucleic Acid Helix-Helix Interactions

Salt ions are essential for the folding of nucleic acids. We use the tightly bound ion (TBI) model, which can account for the correlations and fluctuations for the ions bound to the nucleic acids, to investigate the electrostatic free-energy landscape for two parallel nucleic acid helices in the solution of added salt. The theory is based on realistic atomic structures of the helices. In monovalent salt, the helices are predicted to repel each other. For divalent salt, while the mean-field Poisson-Boltzmann theory predicts only the repulsion, the TBI theory predicts an effective attraction between the helices. The helices are predicted to be stabilized at an interhelix distance ∼26–36 Å, and the strength of the attractive force can reach −0.37 k B T/bp for helix length in the range of 9–12 bp. Both the stable helix-helix distance and the strength of the attraction are strongly dependent on the salt concentration and ion size. With the increase of the salt concentration, the helix-helix attraction becomes stronger and the most stable helix-helix separation distance becomes smaller. For divalent ions, at very high ion concentration, further addition of ions leads to the weakening of the attraction. Smaller ion size causes stronger helix-helix attraction and stabilizes the helices at a shorter distance. In addition, the TBI model shows that a decrease in the solvent dielectric constant would enhance the ion-mediated attraction. The theoretical findings from the TBI theory agree with the experimental measurements on the osmotic pressure of DNA array as well as the results from the computer simulations.

Nonlinear screening and effective electrostatic interactions in charge-stabilized colloidal suspensions

A nonlinear response theory is developed and applied to electrostatic interactions between spherical macroions, screened by surrounding microions, in charge-stabilized colloidal suspensions. The theory describes leading-order nonlinear response of the microions (counterions, salt ions) to the electrostatic potential of the macroions and predicts microion-induced effective many-body interactions between macroions. A linear response approximation [A.R. Denton, Phys. Rev. E 62, 3855 (2000)] yields an effective pair potential of screened-Coulomb (Yukawa) form, as well as a one-body volume energy, which contributes to the free energy. Nonlinear response generates effective many-body interactions and essential corrections to both the effective pair potential and the volume energy. By adopting a random-phase approximation (RPA) for the response functions, and thus neglecting microion correlations, practical expressions are derived for the effective pair and triplet potentials and for the volume energy. Nonlinear screening is found to weaken repulsive pair interactions, induce attractive triplet interactions, and modify the volume energy. Numerical results for monovalent microions are in good agreement with available ab initio simulation data and demonstrate that nonlinear effects grow with increasing macroion charge and concentration and with decreasing salt concentration. In the dilute limit of zero macroion concentration, leading-order nonlinear corrections vanish. Finally, it is shown that nonlinear response theory, when combined with the RPA, is formally equivalent to the mean-field Poisson-Boltzmann theory and that the linear response approximation corresponds, within integral-equation theory, to a linearized hypernetted-chain closure.

A Monte Carlo simulation study of electrostatic forces between hexagonally packed DNA double helices

The electrostatic contribution to the force between hexagonally packed B-DNA double helices has been studied using different statistical mechanical descriptions and the Monte Carlo simulation method. The effects of small ion correlations and sizes are considered, and comparison with experimental results is made. It is found that for monovalent counterions, the mean field Poisson–Boltzmann theory can give a reasonable reproduction of the experimental data, even though the simulations show that its description of the electrostatic interaction can be qualitatively wrong. The chemical equilibrium of the ordered DNA phase with the surrounding bulk salt solution is found to be an important feature. Furthermore, the simulations show that the correlation between the ion clouds of different DNA polyions gives rise to a significant attractive contribution to the interaction when divalent cations are present. It is suggested that this electrostatic attraction is an important driving force for the experimentally observed condensation and collapse of DNA solutions to ordered systems, which can be induced by multivalent ions.