Molecular dynamics simulations of electrolyte solutions at the (100) goethite surface

Molecular dynamics simulations of electrolyte solutions in contact with a neutral (100) goethite (alpha-FeOOH) surface were used to probe the structure of the mineral-water interface and gain insight into the adsorption properties of monovalent ions. Three electrolyte solutions were considered: NaCl, CsCl, and CsF. The electrolyte ions were chosen to cover a range of ionic sizes and affinities for the aqueous phase. The molecular dynamics simulations indicate the presence of a structured interfacial region resulting from the strong interaction of water with the mineral surface. The specific arrangement and preferred orientation of water that arise from this interaction create adsorption sites in the interfacial region, i.e., as far as 15 A away from the surface, and hence give rise to a strong correlation between the water and ion distributions. The structure of the hydrated ion, its effect on the water arrangement at the interface, and the strength of the ion-water bond are found to be key factors that determine the location and extent of ion adsorption at the interface. Additionally, in all simulations, we find a build up of positive charges near the surface due to cation adsorption, which is compensated by an accumulation of anions in the next few angstr?ms. This creates an excess of negative charges, which is in turn compensated by an excess of positive charges, and so on. As we modeled a neutral surface, the structure of the electrolyte distribution arises from the complex interplay of the interactions between the surface, water, and the electrolyte ions rather than from the need to neutralize a surface charge. In addition, our simulations indicate that the electrolyte distribution does not resemble that of a classical electrical double layer. Indeed, our calculations predict the presence of several condensed layers and oscillations in the net charge away from the 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.     

Comparison of cation adsorption by isostructural rutile and cassiterite

Macroscopic net proton charging curves for powdered rutile and cassiterite specimens with the (110) crystal face predominant, as a function of pH in RbCl and NaCl solutions, trace SrCl(2) in NaCl, and trace ZnCl(2) in NaCl and Na Triflate solutions, are compared to corresponding molecular-level information obtained from static DFT optimizations and classical MD simulations, as well as synchrotron X-ray methods. The similarities and differences in the macroscopic charging behavior of rutile and cassiterite largely reflect the cation binding modes observed at the molecular level. Cation adsorption is primarily inner-sphere on both isostructural (110) surfaces, despite predictions that outer-sphere binding should predominate on low bulk dielectric constant oxides such as cassiterite (ε(bulk) ≈ 11). Inner-sphere adsorption is also significant for Rb(+) and Na(+) on neutral surfaces, whereas Cl(-) binding is predominately outer-sphere. As negative surface charge increases, relatively more Rb(+), Na(+), and especially Sr(2+) are bound in highly desolvated tetradentate fashion on the rutile (110) surface, largely accounting for enhanced negative charge development relative to cassiterite. Charging curves in the presence of Zn(2+) are very steep but similar for both oxides, reflective of Zn(2+) hydrolysis (and accompanying proton release) during the adsorption process, and the similar binding modes for ZnOH(+) on both surfaces. These results suggest that differences in cation adsorption between high and low bulk dielectric constant oxides are more subtly related to the relative degree of cation desolvation accompanying inner-sphere binding (i.e., more tetradentate binding on rutile), rather than distinct inner- and outer-sphere adsorption modes. Cation desolvation may be favored at the rutile (110) surface in part because inner-sphere water molecules are bound further from and less tightly than on the cassiterite (110) surface. Hence, their removal upon inner-sphere cation binding is relatively more favorable.     

Effects of ionic strength and surface charge on protein adsorption at PEGylated surfaces

PEGylated Nb2O5 surfaces were obtained by the adsorption of poly(L-lysine)-g-poly(ethylene glycol) (PLL-g-PEG) copolymers, allowing control of the PEG surface density, as well as the surface charge. PEG (MW 2 kDa) surface densities between 0 and 0.5 nm(-2) were obtained by changing the PEG to lysine-mer ratio in the PLL-g-PEG polymer, resulting in net positive, negative and neutral surfaces. Colloid probe atomic force microscopy (AFM) was used to characterize the interfacial forces associated with the different surfaces. The AFM force analysis revealed interplay between electrical double layer and steric interactions, thus providing information on the surface charge and on the PEG layer thickness as a function of copolymer architecture. Adsorption of the model proteins lysozyme, alpha-lactalbumin, and myoglobin onto the various PEGylated surfaces was performed to investigate the effect of protein charge. In addition, adsorption experiments were performed over a range of ionic strengths, to study the role of electrostatic forces between surface charges and proteins acting through the PEG layer. The adsorbed mass of protein, measured by optical waveguide lightmode spectroscopy (OWLS), was shown to depend on a combination of surface charge, protein charge, PEG thickness, and grafting density. At high grafting density and high ionic strength, the steric barrier properties of PEG determine the net interfacial force. At low ionic strength, however, the electrical double layer thickness exceeds the thickness of the PEG layer, and surface charges "shining through" the PEG layer contribute to protein interactions with PLL-g-PEG coated surfaces. The combination of AFM surface force measurements and protein adsorption experiments provides insights into the interfacial forces associated with various PEGylated surfaces and the mechanisms of protein resistance.     

Relaxation of the electrical double layer after an electron transfer approached by Brownian dynamics simulation

In this paper, the dynamical properties of the electrochemical double layer following an electron transfer are investigated by using Brownian dynamics simulations. This work is motivated by recent developments in ultrafast cyclic voltammetry which allow nanosecond time scales to be reached. A simple model of an electrochemical cell is developed by considering a 1:1 supporting electrolyte between two parallel walls carrying opposite surface charges, representing the electrodes; the solution also contains two neutral solutes representing the electroactive species. Equilibrium Brownian dynamics simulations of this system are performed. To mimic electron transfer processes at the electrode, the charge of the electroactive species are suddenly changed, and the subsequent relaxation of the surrounding ionic atmosphere are followed, using nonequilibrium Brownian dynamics. The electrostatic potential created in the center of the electroactive species by other ions is found to have an exponential decay which allows the evaluation of a characteristic relaxation time. The influence of the surface charge and of the electrolyte concentration on this time is discussed, for several conditions that mirror the ones of classical electrochemical experiments. The computed relaxation time of the double layer in aqueous solutions is found in the range 0.1 to 0.4 ns for electrolyte concentrations between 0.1 and 1 mol L(-1) and surface charges between 0.032 and 0.128 C m(-2).     

SANS and SAXS analysis of charged nanoparticle adsorption at oil-water interfaces

A systematic study of the adsorption of charged nanoparticles at dispersed oil-in-water emulsion interfaces is presented. The interaction potentials for negatively charged hexadecane droplets with anionic polystyrene latex particles or cationic gold particles are calculated using DLVO theory. Calculations demonstrate that increased ionic strength decreases the decay length of the electrostatic repulsion leading to enhanced particle adsorption. For the case of anionic PS latex particles, the energy barrier for particle adsorption is also reduced when the surface charge is neutralized through changes in pH. Complementary small-angle scattering experiments show that the highest particle adsorption for PS latex occurs at moderate ionic strength and low pH. For cationic gold particles, simple DLVO calculations also explain scattering results showing that the highest particle adsorption occurs at neutral pH due to the electrostatic attraction between oppositely charged surfaces. This work demonstrates that surface charges of particles and oil droplets are critical parameters to consider when engineering particle-stabilized emulsions.     

Hierarchical multiscale simulation of electrokinetic transport in silica nanochannels at the point of zero charge

Effects of nanoscale confinement and partial charges that stem from quantum calculations are investigated in silica slit channels filled with 1 M KCl at the point of zero charge by using a hierarchical multiscale simulation methodology. Partial charges of both bulk and surface atoms from ab initio quantum calculations that take into account bond polarization and electronegativity are used in molecular dynamics (MD) simulations to obtain ion and water concentration profiles for channel widths of 1.1, 2.1, 2.75, and 4.1 nm. The interfacial electron density profiles of simulations matched well with that of recent X-ray reflectivity experiments. By simulating corresponding channels with no partial charges, it was observed that the partial charges affect the concentration profiles and transport properties such as diffusion coefficients and mobilities up to a distance of about 3 sigma(O)(-)(O) from the surface. Both in uncharged and partially charged cases, oscillations in concentration profiles of K(+) and Cl(-) ions give rise to an electro-osmotic flow in the presence of an external electric field, indicating the presence of an electric double layer at net zero surface charge, contrary to the expectations from classical continuum theory. I-V curves in a channel-bath system using ionic mobilities from MD simulations were significantly different for channels with and without partial charges for channel widths less than 4.1 nm.     

Hydrophobic and ionic interactions in nanosized water droplets

A number of situations such as protein folding in confined spaces, lubrication in tight spaces, and chemical reactions in confined spaces require an understanding of water-mediated interactions. As an illustration of the profound effects of confinement on hydrophobic and ionic interactions, we investigate the solvation of methane and methane decorated with charges in spherically confined water droplets. Free energy profiles for a single methane molecule in droplets, ranging in diameter (D) from 1 to 4 nm, show that the droplet surfaces are strongly favorable as compared to the interior. From the temperature dependence of the free energy in D = 3 nm, we show that this effect is entropically driven. The potentials of mean force (PMFs) between two methane molecules show that the solvent separated minimum in the bulk is completely absent in confined water, independent of the droplet size since the solute particles are primarily associated with the droplet surface. The tendency of methanes with charges (M(q+) and M(q-) with q(+) = |q(-)| = 0.4e, where e is the electronic charge) to be pinned at the surface depends dramatically on the size of the water droplet. When D = 4 nm, the ions prefer the interior whereas for D < 4 nm the ions are localized at the surface, but with much less tendency than for methanes. Increasing the ion charge to e makes the surface strongly unfavorable. Reflecting the charge asymmetry of the water molecule, negative ions have a stronger preference for the surface compared to positive ions of the same charge magnitude. With increasing droplet size, the PMFs between M(q+) and M(q-) show decreasing influence of the boundary owing to the reduced tendency for surface solvation. We also show that as the solute charge density decreases the surface becomes less unfavorable. The implications of our results for the folding of proteins in confined spaces are outlined.     

Discrete solvation layering in confined binary liquids

The solvation force profiles of squalane/octamethylcyclotetrasiloxane (OMCTS) mixtures confined between Si3N4 tips and highly oriented pyrolytic graphite (HOPG) and hexadecane/OMCTS confined between alkanethiol-functionalized tips and freshly cleaved mica have been measured by atomic force microscopy. Measurements on HOPG reveal oscillatory behavior where discrete solvation layers of both squalane and OMCTS are observed in a single force curve. The large repulsive force of the first solvation layer (squalane) on HOPG indicates that it is strongly bound. Oscillatory behavior is also observed for hexadecane/OMCTS on mica excepting that the oscillations are found in the attractive regime. The OMCTS layers in this case are less ordered with slightly larger (approximately 1 A) periodicities. These results are in agreement with computer simulations for binary liquid mixtures but differ qualitatively from surface force apparatus experiments.     

Oscillatory forces of nanoparticle suspensions confined between rough surfaces modified with polyelectrolytes via the layer-by-layer technique

This paper addresses the systematic study of surface roughness effects on the internal structuring of silica nanoparticle suspensions under confinement. The confining surfaces are modified by physisorption of layers of oppositely charged polyelectrolytes with the so-called layer-by-layer technique. The layer-by-layer technique modifies the surface roughness without changing the surface potential of a multilayer with the same outermost layer, by increasing the number of constituent layers and ionic strength of the polyelectrolyte solutions and by selecting an appropriate pair of polyelectrolytes. The oscillatory forces of nanoparticle suspensions with a particle diameter of 26 nm are measured by a colloidal-probe atomic force microscope (CP-AFM). The characteristic lengths of the oscillatory force, i.e., wavelength, which indicates interparticle distance, and decay length, or particle correlation length, are not affected by the surface roughness. The corresponding reduction in the oscillatory amplitude and the shift in the phase correlate with an increase in surface roughness. Increasing surface roughness further induces a disappearance of the oscillations, and both confining surfaces contribute to the effect of surface roughness on the force reduction. In order to show an oscillatory force, the particles have to show positional correlation over a reasonably long range perpendicular to the surface, and the correlation function should be the same over a larger lateral area. This requires that both the particles and the surfaces have a high degree of order or symmetry; otherwise, the oscillation does not occur. A roughness of a few nanometers on a single surface, which corresponds to about 10% of the nanoparticle diameter, is sufficient to eliminate the oscillatory force.     

The role of variable surface charge and surface complexation in the adsorption of humic acid on magnetite

The ionic strength dependence of humic acid (HA) adsorption on magnetite (Fe₃O₄) was investigated at pH 5, 8 and 9, where variable charged magnetite is positive, neutral and negative, respectively. The adsorption studies revealed that HA has high affinity to magnetite surface especially at lower pH, where interacting partners have opposite charges. However, in spite of electrostatic repulsion at pH 9 notable amounts of humate are adsorbed. Increasing ionic strength enhances HA adsorption at each pH due to charge screening. The dominant interaction is probably a ligand-exchange reaction, nevertheless the Coulombic contribution to the organic matter accumulation on oxide surface is also significant under acidic condition. The results from size exclusion chromatography demonstrate that the smaller size HA fractions enriched with functional groups are adsorbed preferentially on the surface of magnetite at pH 8 in dilute NaCl solution.     

基于石英纳米孔道的单颗粒尺寸分布分析

发展了一种基于石英纳米孔道的单颗粒电化学动态分析方法, 用于单个CdSe/ZnS量子点纳米颗粒的尺寸分布分析. 其机制是向石英纳米孔道两端施加电压, 表面带有正电荷的单个CdSe/ZnS量子点纳米颗粒在电场力驱动下由管内向管外运动, 当量子点纳米颗粒穿过纳米孔道尖端狭小的限域空间时, 其表面正电荷使石英纳米孔道内电荷密度增加, 孔道内的电化学限域效应进一步将电荷密度增加的信息放大并转变为可读的离子流增强信号. 通过对动态离子流信号解析可实时获取具有2种不同尺寸的量子点纳米颗粒所导致的2类过孔事件信息, 从而对在限域空间内运动的纳米颗粒进行尺寸分布分析.     

Co-ion and ion competition effects: ion distributions close to a hydrophobic solid surface in mixed electrolyte solutions

We consider within a modified Poisson-Boltzmann theory an electrolyte, with different mixtures of NaCl and NaI, near uncharged and charged solid hydrophobic surfaces. The parametrized potentials of mean force acting on Na+, Cl-, and I- near an uncharged self-assembled monolayer were deduced from molecular simulations with polarizable force fields. We study what happens when the surface presents negative charges. At moderately charged surfaces, we observe strong co-ion adsorption and clear specific ion effects at biological concentrations. At high surface charge densities, the co-ions are pushed away from the interface. We predict that Cl- ions can also be excluded from the surface by increasing the concentration of NaI. This ion competition effect (I- versus Cl-) may be relevant for ion-specific partitioning in multiphase systems where polarizable ions accumulate in phases with large surface areas.     

Dressed ion theory of size-asymmetric electrolytes: effective ionic charges and the decay length of screened Coulomb potential and pair correlations

The effects of ionic size asymmetry on long-range electrostatic interactions in electrolyte solutions are investigated within the primitive model. Using the formalism of dressed ion theory we analyze correlation functions from Monte Carlo simulations and the hypernetted chain approximation for size asymmetric 1:1 electrolytes. We obtain decay lengths of the screened Coulomb potential, effective charges of ions, and effective permittivity of the solution. It is found that the variation of these quantities with the degree of size asymmetry depends in a quite intricate manner on the interplay between the electrostatic coupling and excluded volume effects. In most cases the magnitude of the effective charge of the small ion species is larger than that of the large species; the difference increases with increasing size asymmetry. The effective charges of both species are larger (in absolute value) than the bare ionic charge, except for high asymmetry where the effective charge of the large ions can become smaller than the bare charge.     

Electrooxidation of Glucose at Nanoporous Gold Surfaces: Structure Dependent Electrocatalysis and Its Application to Amperometric Detection

Electrooxidation of glucose is investigated at nanoporous gold (NPG) with controlled surface structures by applying different deposition charges during the formation of Ag? Au layers. As the deposition charge increases, the NPG surfaces exhibit smaller ligament/pore structures and the electrocatalytic oxidation of glucose becomes more effective. Voltammetric responses of NPG suggest that the electrocatalytic oxidation arises from the enrichment of (110) or (100) surface orientation of gold with higher deposition charges. The electrooxidation of glucose is retained at NPG surfaces with higher deposition charges in the presence of Cl^?, which suggests possible applications to the amperometric glucose detection in biological samples.     

van der Waals forces in presence of free charges: an exact derivation from equilibrium quantum correlations

We study interatomic forces in a fluid consisting of a mixture of free charges and neutral atoms in the framework of the quantum many-body problem at nonzero temperature and nonzero density. Of central interest is the interplay between van der Waals forces and screening effects due to free charges. The analysis is carried out in a partially recombined hydrogen plasma in the Saha regime. The effective potentials in the medium between two atoms, or an atom and a charge, or two charges, are determined from the large-distance behavior of equilibrium proton-proton correlations. We show, in a proper low-temperature and low-density scaling limit, that those potentials all decay as r(-6) at large distance r, while the corresponding amplitudes are calculated exactly. In particular, the presence of free charges only causes a partial (nonexponential) screening of the atomic potential, and it does not modify its typical r(-6) decay. That potential reduces to the standard van der Waals form for two atoms in vacuum when the temperature is driven to zero. The analysis is based on first principles: it does not assume preformed atoms and takes into account in a coherent way all effects, quantum mechanical binding, ionization, and collective screening, which originate from the Coulomb potential. Our method relies on the path integral representation of the quantum Coulomb gas.     

Closely adjacent gold nanoparticles linked by chemisorption of neutral rhodamine 123 molecules providing enormous SERS intensity

Addition of neutral R123 molecules (10(-7) M) to an as-prepared gold nanoparticles (AuNPs) suspension generated flocculates that are a small number of closely adjacent particles. Formation of AuNP flocculates was evidenced by the coupled localized plasmon peak at 720-750 nm. The AuNP flocculates provided pronounced SERS spectra of adsorbed neutral R123 molecules (SERS-A) as anticipated by FDTD (Finite Difference Time Domain) simulations. The observed SERS spectra are significantly different from those of cationic R123(+) molecules (SERS-B), which electrostatically adsorbed on Cl(-)-treated AuNPs. The difference is not simply due to deprotonation but reflects a distinct difference in adsorption nature between neutral R123 and cationic R123(+) molecules. Indeed neutral R123 molecules exclusively gave an Au-N stretching band at 202 cm(-1), showing the chemisorption on Au surfaces through lone pair electrons at the amino groups. The different adsorption nature is further evidenced by the observation that cationic R123(+) molecules adsorbed on as-prepared (without NaCl addition) AuNP flocculates gave both SERS-A and SERS-B spectra. Thus, the cationic R123(+) molecules form the flocculates both by chemisorption and electrostatic adsorption owing to modest surface charge on as-prepared AuNPs.     

Electro-osmosis at inhomogeneous charged surfaces: hydrodynamic versus electric friction

Electrokinetic methods are efficient in probing the electrostatic surface properties of charged systems. However, anomalies observed in experiments indicate that the classical electrokinetic theory should be reconsidered. Using Green's function methods and hydrodynamic simulations, we investigate electro-osmosis driven by electric-field-induced ion motion near a charged planar substrate with smooth or rough boundary. First, a reformulation of electro-osmotic theory for planar charged surfaces employing Green's functions shows that the Helmholtz-Smoluchowski (HS) relation between electrostatic potential and solvent velocity is exact for smooth surfaces, even in the presence of ion correlations. Deviations from HS theory are caused by combined hydrodynamic and electric surface friction, as our hydrodynamic simulations of ions at smooth and corrugated charged surfaces in lateral electric fields demonstrate. Within the simulations, hydrodynamic interactions are treated in the continuum limit and the presence of a no-slip boundary condition at the surface is taken into account. While electrofriction is relevant in highly charged system and/or for multivalent ions, hydrodynamic friction is dominant in systems with moderate surface charge density and/or low ionic valency. We also derive the effective electrokinetic surface charge from the electro-osmotic solvent profiles, which is substantially reduced when compared with the bare value and shows qualitative agreement with the experimental tendency.     

Interactions between charged surfaces immersed in a polyelectrolyte solution

With grand canonical simulations invoking a configurationally weighted scheme, we have calculated interactions between charged surfaces immersed in a polyelectrolyte solution. In contrast to previous simulations of such systems, we have imposed full equilibrium conditions (i.e., we have included diffusive equilibrium with a bulk solution). This has a profound impact on the resulting interactions: even at modest surface charge densities, oppositely charged chains will, at sufficiently large separations, adsorb strongly enough to overcompensate for the nominal surface charge. This phenomenon, known as charge inversion, generates a double-layer repulsion and a free-energy barrier. Simpler canonical approaches, where the chains are assumed to neutralize the surfaces perfectly, will not capture this stabilizing barrier. The barrier height increases with the length of the polyions. Interestingly enough, the separation at which the repulsion becomes attractive is independent of chain length. The short chains here are unable to reach across from one surface to the other. We therefore conclude that the transition to an attractive regime is not provided by the formation of such "intersurface" bridges. With long chains and at large separations, charge inversion displays decaying oscillatory behavior (i.e., the apparent surface charge switches sign once again). This is due to polyion packing effects. We have also investigated responses to salt addition and changes in polyelectrolyte concentration. Our results are in qualitative and semiquantitative agreement with experimental findings, although it should be noted that our chains are comparatively short, and the experimental surface charge density is poorly established.