Ionic solvation at the solid/electrolyte interface: A statistical-mechanical approach

We have studied studied the influence of the size of ions on their adsorbability at a solid surface in the presence of a molecular solvent. Ions and molecules are represented respectively by charged hard spheres and dipolar hard spheres and the surface is just a neutral hard wall. We have found that the electrostatic interaction between ions and molecules can induce the exclusion of small ions from the surface. A pure MSA (mean spherical approximation) calculation would not give any effect of the solvation on the ionic density profile. The present calculation is limited to the case of infinite ionic dilution.     

Effect of attractions on the structure of polymer solutions confined between surfaces: a density functional approach

A density functional theory is presented to study the effect of attractions on the structure of polymer solutions confined between surfaces. The polymer molecules have been modeled as a pearl necklace of freely jointed hard spheres and the solvent as hard spheres, both having Yukawa-type attractions and the mixture being confined between attractive Yukawa-type surfaces. The present theory treats the ideal gas free energy functional exactly and uses weighted density approximation for the hard chain and hard sphere contributions to the excess free energy functional. The attractive interactions are calculated using the direct correlation function obtained from the polymer reference interaction site model theory along with the mean spherical approximation closure. The theoretical predictions on the density profiles of the polymer and the solvent molecules are found to agree quite well with the Monte Carlo simulation results for varying densities, chain lengths, wall separations, and different sets of interaction potentials.     

Forces between surfaces in electrolyte solutions

Density functional theories of solvation forces in charged fluids are extended to treat electrolytes consisting of finitesized ions and neutral solvent particles. The resulting forces display pronounced oscillations whose magnitude is a strong function of the bulk density of the neutral species.     

Solvation forces in asymmetrically charged electrolytes

Density functional theories of solvation forces in charged fluids which self-consistently include the effects of finitesized ions are extended to treat asymmetrically charged electrolytes. For a given electrolyte concentration both the electric potential and the solvation force between charged surfaces are seen to deviate from the results of the classical DLVO theory more markedly for asymmetrically charged electrolytes than for a 1 : 1 electrolyte.     

Charged versus Neutral Hydrogen‐Bonded Complexes: Is There a Difference in the Nature of the Hydrogen Bonds?

A theoretical study on some carboxylic acid dimers formed by positively or negatively charged molecules has been carried out by using DFT methods. The resulting dimers possess either a charge of +2 or ?2. In addition, the corresponding neutral complexes have also been considered. The electron density distribution described by the atoms in molecules and the natural bond orbital methods, as well as the electric field maps of the systems, have been analyzed and compared without finding significant differences between the neutral and ionic complexes. The interaction energy along the dissociation path of the charged dimers shows both a local minimum and a local maximum, defining a stability region between them. When this energetic profile is recalculated by removing the repulsion between the charged groups, it resembles to those of the neutral molecules. Hence, the characteristics of the charged dimers are similar to those of the neutral ones: the addition of a repulsion term for the charged groups permits to retrieve the energetic profiles dependence with the distance in the charged system. The interacting quantum atom (IQA) method has been used to calculate the interaction energy terms, including the classic Coulombic term between the whole molecules and the corresponding of the carboxylic acid groups. The IQA results show repulsive electrostatic interactions when the whole molecules are considered in the ionic complexes, but attractive ones between the carboxylic groups in both neutral and ionic complexes.     

DLVO interactions of tungsten oxide and cobalt oxide surfaces measured with the colloidal probe technique

We have investigated the DLVO surface forces of oxidized tungsten and cobalt surfaces using the atomic force microscope (AFM) colloidal probe technique. It was shown by X-ray photoelectron spectroscopy (XPS) and electrokinetic measurements that this model system is representative of industrial tungsten carbide (WC) and cobalt powders used in the production of hard metals. We found that the attractive van der Waals forces are well described by Hamaker constants, calculated from optical data for WO(3) and CoOOH. The repulsive electrostatic double layer forces between WO(3) surfaces increase with increasing pH due to an increasingly negative surface potential. This surface potential decreases with increasing ionic strength at pH 7.5. The electrostatic interaction between WO(3) and CoOOH is attractive at pH 10, suggesting a positively charged CoOOH surface.     

Excluded volume effects in macromolecular forces and ion-interface interactions

A charged Yukawa liquid confined in a slit nanopore is studied in order to understand excluded volume effects in the interaction force between the pore walls. A previously developed self-consistent scheme [S. Buyukdagli, C. V. Achim, and T. Ala-Nissila, J. Stat. Mech. 2011, P05033] and a new simpler variational procedure that self-consistently couple image forces, surface charge induced electric field, and pore modified core interactions are used to this aim. For neutral pores, it is shown that with increasing pore size, the theory predicts a transition of the interplate pressure from an attractive to a strongly repulsive regime associated with an ionic packing state, an effect observed in previous Monte Carlo simulations for hard core charges. We also establish the mean-field theory of the model and show that for dielectrically homogeneous pores, the mean-field regime of the interaction between the walls corresponds to large pores of size d > 4 ?. The role of the range of core interactions in the ionic rejection and interplate pressure is thoroughly analyzed. We show that the physics of the system can be split into two screening regimes. The ionic packing effect takes place in the regime of moderately screened core interactions characterized with the bare screening parameter of the Yukawa potential b ? 3/l(B), where l(B) is the Bjerrum length. In the second regime of strongly screened core interactions b ? 3/l(B), solvation forces associated with these interactions positively contribute to the ionic rejection driven by electrostatic forces and enhance the magnitude of the attractive pressure. For weakly charged pores without a dielectric discontinuity, core interactions make a net repulsive contribution to the interplate force and also result in oscillatory pressure curves, whereas for intermediate surface charges, these interactions exclusively strengthen the external pressure, thereby reducing the magnitude of the net repulsive interplate force. The pronounced dependence of the interplate pressure and ionic partition coefficients on the magnitude and the range of core interactions indicates excluded volume effects as an important ion specificity and a non-negligible ingredient for the stability of macromolecules in electrolyte solutions.     

Theory and simulations of charged polymers: From solution properties to polymeric nanomaterials

Charged polymers are macromolecules with ionizable groups. These polymeric systems demonstrate unique properties that are qualitatively different from their neutral counterparts. In this review I survey the recent progress made in understanding properties of the solutions of charged polymers, swelling of polyelectrolyte gels, conformational transformations of charged dendrimers, complexation between charged macromolecules, adsorption of charged polymers at surfaces and interfaces, and multilayer assembly in ionic systems.     

Charge-dependent cavity radii for an accurate dielectric continuum model of solvation with emphasis on ions: aqueous solutes with oxo, hydroxo, amino, methyl, chloro, bromo, and fluoro functionalities

Dielectric continuum solvation models are widely used because they are a computationally efficacious way to simulate equilibrium properties of solutes. With advances that allow for molecular-shaped cavities, they have reached a high level of accuracy, in particular for neutral solutes. However, benchmark tests show that existing schemes for defining cavities are unable to consistently predict accurately the effects of solvation on ions, especially anions. This work involves the further development of a protocol put forth earlier for defining the cavities of aqueous solutes, with resulting advances that are most striking for anions. Molecular cavities are defined as interlocked spheres around atoms or groups of atoms in the solute, but the sphere radii are determined by simple empirically based expressions involving the effective atomic charges of the solute atoms (derived from molecular electrostatic potential) and base radii. Both of these terms are optimized for the different types of atoms or functional groups in a training set of neutral and charged solutes. Parameters in these expressions for radii were fitted by minimizing residuals between calculated and measured standard free energies of solvation (DeltaG(s)*), weighted by the uncertainty in the measured value. The calculations were performed using density functional theory with the B3LYP functional and the 6-311+G** basis set and the COnductor-like Screening MOdel (COSMO). The optimized radii definitions reproduce DeltaG(s)* of neutral solutes and singly charged ions in the training set to within experimental uncertainty and, more importantly, accurately predict DeltaG(s)* of compounds outside the training set, in particular anions (J. Phys. Chem. A 2003, 107, 5778). Inherent to this approach, the cavity definitions reflect the strength of specific solute-water interactions. We surmise that this feature underlies the success of the model, referred to as the CD-COSMO model for Charge-Dependent (also Camaioni-Dupuis) COSMO model. These findings offer encouragement that we can keep extending this scheme to other functional groups and obtain better accuracy in using continuum solvation models to predict equilibrium properties of aqueous ionic solutes. The approach is illustrated for a number of test cases, including the determination of acidities of an amine base, a study of the tautomerization equilibrium of a zwitterionic molecule (glycine), and calculating solvation energies of transition states toward a full characterization of reaction pathways in aqueous phase, here in S(N)2 exchange reactions. The calculated reaction barriers in aqueous solution are in excellent agreement with experimental values.     

Mesoscale simulations of the behavior of charged polymer brushes under normal compression and lateral shear forces

Dissipative particle dynamics (DPD) was used to investigate the behavior of two opposing end-grafted charged polymer brushes in aqueous media under normal compression and lateral shear. The effect of polymer molecular weight, degree of ionization, grafting density, ionic strength, and compression on the polymer conformation and the resulting shear force between the opposing polymer layers were investigated. The simulations were carried out for the poly(tert-butyl methacrylate)-block-poly(sodium sulfonate glycidyl methacrylate) copolymer, referred as PtBMA-b-PGMAS, end-attached to a hydrophobic surface for comparison with previous experimental data. Mutual interpenetration of the opposing end-grafted chains upon compression is negligible for highly charged polymer brushes for compression ratios ranging from 2.5 to 0.25. Under electrostatic screening effects or for weakly charged polymer brushes, a significant mutual interpenetration was measured. The variation of interpenetration thickness with separation distance, grafting density, and polymer size follows the same scaling law as the one observed for two opposing grafted neutral brushes in good solvent. However, compression between two opposing charged brushes results in less interpenetration relative to neutral brushes when considering equivalent grafting density and molecular weight. The friction coefficient between two opposing polymer-coated surfaces sliding past each other is shown to be directly correlated with the interpenetration thickness and more specifically to the number of polymer segments within the interpenetration layer.     

分子动力学模拟多巴在自组装膜上的黏附性

为了更好地理解贻贝在表面的黏附机理,实现水下胶黏,采用分子动力学方法研究了多巴在自组装膜上的黏附性:采用伞形取样和加权柱状图分析方法计算了多巴在不同自组装膜表面的黏附自由能,使用拉伸分子动力学模拟研究了多巴在不同自组装膜表面上黏附后的脱附力.结果表明,多巴在带负电的羧基自组装膜上的黏附能比在带正电的氨基自组装膜上的大,多巴更容易黏附到带负电表面;多巴在带电表面的黏附能比未带电表面的黏附能更强,表明在带电表面黏附更稳定.进一步分析了多巴在不同表面的取向分布,发现多巴与不同表面相互作用的方式不同:与疏水表面主要通过苯环相互作用;与亲水表面主要通过羟基相互作用;与负电表面主要通过氨基相互作用;与正电表面主要通过羧基相互作用.通过模拟比较了多巴在不同自组装膜上的脱附力,发现多巴在带电表面的脱附力比在未带电表面的大,与黏附能的趋势一致.对比4种非带电表面的脱附力,发现多巴在疏水性甲基自组装膜表面的脱附力最大,黏附更稳定,随着表面疏水性的增加,脱附力增大,黏附稳定性增强.本工作可为研发新型水下胶黏剂提供理论指导.     

Streptococcus mutans and Streptococcus intermedius adhesion to fibronectin films are oppositely influenced by ionic strength

Bacterial adhesion to protein-coated surfaces is mediated by an interplay of specific and nonspecific interactions. Although nonspecific interactions are ubiquitously present, little is known about the physicochemical mechanisms of specific interactions. The aim of this paper is to determine the influence of ionic strength on the adhesion of two streptococcal strains to fibronectin films. Streptococcus mutans LT11 and Streptococcus intermedius NCTC11324 both possess antigen I/II with the ability to bind fibronectin from solution, but S. intermedius binds approximately 20x less fibronectin than does the S. mutans strain under identical conditions. Both strains as well as fibronectin films are negatively charged in low ionic strength phosphate buffered saline (PBS, 10x diluted), but bacteria appear uncharged in high ionic strength PBS. Physicochemical modeling on the basis of overall cell surface properties (cell surface hydrophobicity and zeta potentials) demonstrates that both strains should favor adhesion to fibronectin films in a high ionic strength environment as compared to in a low ionic strength environment, where electrostatic repulsion between equally charged surfaces is dominant. Adhesion of S. intermedius to fibronectin films in a parallel plate flow chamber was completely in line with this modeling, while in addition atomic force microscopy (AFM) indicated stronger adhesion forces upon retraction between fibronectin-coated tips and the cell surfaces in high ionic strength PBS than in low ionic strength PBS. Thus, the dependence of the interaction on ionic strength is dominated by the overall negative charge on the interacting surfaces. Adhesion of S. mutans to fibronectin films, however, was completely at odds with theoretical modeling, and the strain adhered best in low ionic strength PBS. Moreover, AFM indicated weaker repulsive forces upon approach between fibronectin-coated tips and the cell surfaces in low ionic strength PBS than in high ionic strength PBS. This indicated that the dependence of the interaction on ionic strength is dominated by electrostatic attraction between oppositely charged, localized domains on the interacting surfaces, despite their overall negative charge. In summary, this study shows that physicochemical modeling of bacterial adhesion to protein-coated surfaces is only valid provided the number of specific interaction sites on the cell surfaces is low, such as on S. intermedius NCTC11324. Nonspecific interactions are dominated by specific interactions if the number of specific interaction sites is large, such as on S. mutans LT11. Its ionic strength dependence indicates that the specific interaction is electrostatic in nature and operative between oppositely charged domains on the interacting surfaces, despite the generally overall negatively charged character of the surfaces.     

Friction between brush layers of charged and neutral bottle-brush macromolecules. molecular dynamics simulations

Using molecular dynamics simulations, we study the lubricating properties of neutral and charged bottle-brush coatings as a function of the compression and shear stresses and brush grafting density. Our simulations have shown that in charged bottle-brush systems under shear there is a layer with excess counterions located in the middle between brush-bearing surfaces. The main deformation mode of the charged bottle-brush layers is associated with the backbone deformation, resulting in the backbone deformation ratio, α, and shear viscosity, η, being universal functions of the Weissenberg number. In the case of neutral bottle-brush systems, in addition to the backbone deformation there is also side chain deformation. The coupling between backbone and side chain deformation violates universality in the deformation ratio, α, dependence on the Weissenberg number and results in scaling exponents varying with the compression stress and brush grafting density. The existence of different length scales controlling deformation of neutral bottle brushes manifests itself in the shear viscosity, η, dependence on the shear rate, ?γ. Shear viscosity, η, as a function of the shear rate, ?γ, has two plateaus and two shear thinning regimes. The low shear rate plateau and shear thinning regime correspond to the backbone deformation, while the second plateau and shear thinning regime at moderate shear rates are due to side chain deformation. For both systems the value of the friction coefficient increases with increasing shear rate. The values of the friction coefficient for charged bottle-brush systems are about ten times smaller than corresponding values for neutral systems at the same shear rate.     

Single ion diffusive transport within a Poly(styrene sulfonate) polymer brush matrix probed by fluorescence correlation spectroscopy

Diffusive transport within complex environments is a critical piece of the chemistry occurring in such diverse membrane systems as proton exchange and bilayer lipid membranes. In the present study, fluorescence correlation spectroscopy was used to evaluate diffusive charge transport within a strong polyelectrolyte polymer brush. The fluorescent cation rhodamine-6G was used as a counterion probe molecule, and the strong polyelectrolyte poly(styrene sulfonate) was the polymer brush. Such strong polyelectrolyte brushes show promise for charge storage applications, and thus it is important to understand and tune their transport efficiencies. The polymer brush demonstrated preferential solvation of the probe counterion as compared to solvation by the aqueous solvent phase. Additionally, diffusion within the polymer brush was strongly inhibited, as evidenced by a decrease in diffusion constant of 4 orders of magnitude. It also proved possible to tune the transport characteristics by controlling the solvent pH, and thus the ionic strength of the solvent. The diffusion characteristics within the charged brush system depend on the brush density as well as the effective interaction potential between the probe ions and the brush. In response to changes in ionic strength of the solution, it was found that these two properties act in opposition to each other within this strong polyelectrolyte polymer brush environment. A stochastic random walk model was developed to simulate interaction of a diffusing charged particle with a periodic potential, to show the response of characteristic diffusion times to electrostatic field strengths. The combined results of the experiments and simulations demonstrate that responsive diffusion characteristics in this brush system are dominated by changes in Coulombic interactions rather than changes in brush density. More generally, these results support the use of FCS to evaluate local charge transport properties within polyelectrolyte brush systems, and demonstrate that the technique shows promise in the development of novel polyelectrolyte films for charge storage/transport materials.     

A novel method to describe the interaction pressure between charged plates with application of the weighted correlation approach

Based on the Euler-Lagrange equation for ion density distribution in an inhomogeneous, charged, and hard-sphere fluid, a novel method is proposed to determine the interaction pressure between charged plates. The resulting expression is a sum of distinct physical contributions to the pressure, which involves different contributions to the single-particle direct correlation function. It can, therefore, be conveniently used in any density functional approach to facilitate analysis of the pressure components. In this study, the so-called fundamental measure theory (FMT)∕weighted correlation approach (WCA) approach is applied to estimate both the hard-sphere and the electric residual contributions to the single-particle direct correlation function, upon the calculation of the ionic density profiles between charged plates. The results, against the Monte Carlo simulations, show that the FMT∕WCA approach is superior to the typical FMT∕mean spherical approximation approach of the density functional theory in predicting the interaction pressure between charged plates immersed in an electrolyte solution upon various conditions in the primitive model. The FMT∕WCA approach can well capture the fine features of the pressure-separation dependence, to reproduce not only the shoulder shape and the weak attractions in monovalent electrolytes but also the strongly oscillatory behavior of pressure in divalent electrolytes where pronounced attractions are observed. In addition, it is found that the FMT∕WCA approach even has an advantage over the anisotropic, hyper-netted chain approach in that it agrees with the Monte Carlo results to a very good extent with, however, much less computational effort.     

Friction and normal interaction forces between irreversibly attached weakly charged polymer brushes

Polyelectrolyte brushes were built on mica by anchoring polystyrene-poly(acrylic acid) (PS-b-PAA) diblock copolymers at a controlled surface density in a polystyrene monolayer covalently attached to OH-activated mica surfaces. Compared to physisorbed polymer brushes, these irreversibly attached charged brushes allow the polymer grafting density to remain constant upon changes in environmental conditions (e.g., pH, salt concentration, compression, and shear). The normal interaction and friction forces as a function of surface separation distance and at different concentrations of added salt (NaCl) were investigated using a surface forces apparatus. The interaction force profiles were completely reversible both on loading and receding and were purely repulsive. For a constant polymer grafting density, the influence of the polyelectrolyte charges and the Debye screening effect on the overall interaction forces was investigated. The experimental interaction force profiles agree very well with scaling models developed for neutral and charged polymer brushes. The variation of the friction force between two PAA brushes in motion with respect to each other as a function of surface separation distance appeared to be similar to that observed with neutral brushes. This similarity suggests that the increase in friction is associated with an increase in mutual interpenetration upon compression as observed with neutral polymers. The effect of the PAA charges and added ions was more significant on the repulsive normal forces than on the friction forces. The reversible characteristics of the normal force profiles and friction measurements confirmed the strong attachment of the PAA brushes to the mica substrate. High friction coefficients (ca 0.3) were measured at relatively high pressures (40 atm) with no surface damage or polymer removal.     

Electrostatic interactions and ionization effect in immobilized artificial membrane retention. A comparative study with octanol-water partitioning

The present study is a continuation of our efforts to investigate the effect of electrostatic interactions and ionization on immobilized artificial membrane (IAM) retention. The previous set of neutral and basic drugs was extended to include acids and ampholytes and analogous buffer conditions in the mobile phase were used, namely morpholinepropanesulfonic acid and phosphate buffer saline, adjusted at pH 7.4. The important contribution of electrostatic forces in IAM retention of positively charged species was further justified by the results of the present study, while analogous electrostatic interactions for ionized acidic drugs were not found to affect significantly the affinity for the IAM stationary phase. The critical role of shielding or exposure of the charged centers on the IAM surface, as a result of the effect of the aqueous component of the mobile phase, was evaluated by the use of water instead of buffer for a number of drugs. Measurements at pH 5.0 demonstrated the effect of ionization in IAM retention despite the partial compensation by electrostatic interactions in the case of protonated basic drugs. Silanophilic interactions were also found to play a potential role as secondary interactions in IAM retention. IAM chromatographic indices were compared to octanol-water distribution coefficients and the corresponding relationships established. Finally, solvation analysis was applied in the aim to gain insight in the balance of forces between IAM retention and octanol-water partitioning. The results showed that apart from electrostatic interactions, there is no significant differentiation between the two systems.     

Surface excitations involving headgroup conformations and molecular protrusions. a monte carlo study of surface forces.

The force between two parallel zwitterionic surfaces has been calculated using Monte Carlo computer simulations. The zwitterions are modelled as two oppositely charged hard spheres joined by a string of length L with parameters chosen to mimic a phospholipid system. All centers interact by a homogeneous Coulomb interaction and by a hard sphere exclusion. The anchoring of the negative centers to the surfaces has been treated within two different models. They were either anchored by a parabolic potential or by a protrusion potential, i.e., a potential proportional to the distance between the center and the surface. The latter model gives a more realistic picture of the interaction between amphiphilic surfaces in aqueous solution but here only repulsive forces could be calculated. The first model also allowed the calculation of attractive forces.For distances D between the surfaces, as defined by the location of the negative centers, that are larger than 2L there is an attractive force of the classical van der Waals type. When, on the other hand, D < 2L a strong repulsive force appears, which in the limit D ⪡ 2L is analogous to a double layer force.Recently it was suggested (Israelachvili and Wennerström, Langmuir, 6 (1990) 873 ) that the repulsive so called hydration force observed for biological lipid systems has its origin in confinements of surface excitations induced by a second surface. Here we demonstrate how this mechanism works in a particular microscopic model of the surface and we show that it gives an important contribution to the total force. Although still simplistic several qualitative features of the force in the phospholipid systems are reproduced in the calculations. For example, a reduction of the size of the charged centers leads to a decrease in the repulsive force. This mimics the observed difference between phosphatidyl choline and phosphatidyl ethanolamine.     

Direct surface force measurements of polyelectrolyte multilayer films containing nanocrystalline cellulose

Polyelectrolyte multilayer films containing nanocrystalline cellulose (NCC) and poly(allylamine hydrochloride) (PAH) make up a new class of nanostructured composite with applications ranging from coatings to biomedical devices. Moreover, these materials are amenable to surface force studies using colloid-probe atomic force microscopy (CP-AFM). For electrostatically assembled films with either NCC or PAH as the outermost layer, surface morphology was investigated by AFM and wettability was examined by contact angle measurements. By varying the surrounding ionic strength and pH, the relative contributions from electrostatic, van der Waals, steric, and polymer bridging interactions were evaluated. The ionic cross-linking in these films rendered them stable under all solution conditions studied although swelling at low pH and high ionic strength was inferred. The underlying polymer layer in the multilayered film was found to dictate the dominant surface forces when polymer migration and chain extension were facilitated. The precontact normal forces between a silica probe and an NCC-capped multilayer film were monotonically repulsive at pH values where the material surfaces were similarly and fully charged. In contrast, at pH 3.5, the anionic surfaces were weakly charged but the underlying layer of cationic PAH was fully charged and attractive forces dominated due to polymer bridging from extended PAH chains. The interaction with an anionic carboxylic acid probe showed similar behavior to the silica probe; however, for a cationic amine probe with an anionic NCC-capped film, electrostatic double-layer attraction at low pH, and electrostatic double-layer repulsion at high pH, were observed. Finally, the effect of the capping layer was studied with an anionic probe, which indicated that NCC-capped films exhibited purely repulsive forces which were larger in magnitude than the combination of electrostatic double-layer attraction and steric repulsion, measured for PAH-capped films. Wherever possible, DLVO theory was used to fit the measured surface forces and apparent surface potentials and surface charge densities were calculated.