Molecular dynamics simulations of multilayer polyelectrolyte films: effect of electrostatic and short-range interactions

The effect of the strength of electrostatic and short-range interactions on the multilayer assembly of oppositely charged polyelectrolytes at a charged substrate was studied by molecular dynamics simulations. The multilayer buildup was achieved through sequential adsorption of charged polymers in a layer-by-layer fashion from dilute polyelectrolyte solutions. The strong electrostatic attraction between oppositely charged polyelectrolytes at each deposition step is a driving force behind the multilayer growth. Our simulations have shown that a charge reversal after each deposition step is critical for steady multilayer growth and that there is a linear increase in polymer surface coverage after the first few deposition steps. Furthermore, there is substantial intermixing between chains adsorbed during different deposition steps. We show that the polymer surface coverage and multilayer structure are each strongly influenced by the strength of electrostatic and short-range interactions.     

Forces between nanorods with end-adsorbed chains in a homopolymer melt

Adsorbed or grafted polymers are often used to provide steric stabilization of colloidal particles. When the particle size approaches the nanoscale, the curvature of the particles becomes relevant. To investigate this effect for the case of cylindrical symmetry, I use a classical fluids density functional theory applied to a coarse-grained model to study the polymer-mediated interactions between two nanorods. The rods are coated with end-adsorbing chains and immersed in a polymer melt of chemically identical, nonadsorbing chains. The force between the nanorods is found to be nonmonotonic, with an attractive well when the two brushes come into contact with each other, followed by a steep repulsion at shorter distances. The attraction is due to the entropic phenomenon of autophobic dewetting, in which there is a surface tension between the brush and the matrix chains. These results are similar to previous results for planar and spherical polymer brushes in melts of the same polymer. The depth of the attractive well increases with matrix chain molecular weight and with the surface coverage. The attraction is very weak when the matrix chain molecular weight is similar to or smaller than the brush molecular weight, but for longer matrix chains the magnitude of the attraction can become large enough to cause aggregation of the nanorods.     

Surface forces between telechelic brushes revisited: the origin of a weak attraction

Telechelic polymers are useful for surface protection and stabilization of colloidal dispersions by the formation of polymer brushes. A number of theoretical investigations have been reported on a weak attraction between two telechelic brushes when they are at the classical contact, i.e., when the surface separation is approximately equal to the summation of the brush thicknesses. While recent experiments have confirmed the weak attraction between telechelic brushes, its origin remains elusive because of conflicting approximations used in the previous theoretical calculations. In this paper, we have investigated the telechelic polymer-mediated surface forces by using a polymer density functional theory (PDFT) that accounts for both the surface-adhesive energy and segment-level interactions specifically. Within a single theoretical framework, the PDFT is able to capture both the depletion-induced attraction in the presence of weakly adhesive polymers and the steric repulsion between compressed polymer brushes. In comparison of the solvation forces between telechelic brushes with those between brushes formed by surfactant-like polymers and with those between two asymmetric surfaces mediated by telechelic polymers, we conclude that the weak attraction between telechelic brushes is primarily caused by the bridging effect. Although both the surfactant-like and telechelic polymers exhibit a similar scaling behavior for the brush thickness, a significant difference has been observed in terms of the brush microstructures, in particular, the segment densities near the edges of the polymer brushes.     

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.     

Molecular dynamics simulations of polyelectrolyte adsorption

We have performed molecular dynamics simulations of polyelectrolyte adsorption at oppositely charged surfaces from dilute polyelectrolyte solutions. In our simulations, polyelectrolytes were modeled by chains of charged Lennard-Jones particles with explicit counterions. We have studied the effects of the surface charge density, surface charge distribution, solvent quality for the polymer backbone, strength of the short-range interactions between polymers and substrates on the polymer surface coverage, and the thickness of the adsorbed layer. The polymer surface coverage monotonically increases with increasing surface charge density for almost all studied systems except for the system of hydrophilic polyelectrolytes adsorbing at hydrophilic surfaces. In this case the polymer surface coverage saturates at high surface charge densities. This is due to additional monomer-monomer repulsion between adsorbed polymer chains, which becomes important in dense polymeric layers. These interactions also preclude surface overcharging by hydrophilic polyelectrolytes at high surface charge densities. The thickness of the adsorbed layer shows monotonic dependence on the surface charge density for the systems of hydrophobic polyelectrolytes for both hydrophobic and hydrophilic surfaces. Thickness is a decreasing function of the surface charge density in the case of hydrophilic surfaces while it increases with the surface charge density for hydrophobic substrates. Qualitatively different behavior is observed for the thickness of the adsorbed layer of hydrophilic polyelectrolytes at hydrophilic surfaces. In this case, thickness first decreases with increasing surface charge density, then it begins to increase.     

Chitosan based surfactant polymers designed to improve blood compatibility on biomaterials

We developed chitosan based surfactant polymers that could be used to modify the surface of existing biomaterials in order to improve their blood compatibility. These polymers consist of a chitosan backbone, PEG side chains to repel non-specific protein adsorption, and hexanal side chains to facilitate adsorption and proper orientation onto a hydrophobic substrate via hydrophobic interactions. Since chitosan is a polycationic polymer, and it is thrombogenic, the surface charge was altered to determine the role of this charge in the hemocompatibility of chitosan. Charge had a notable effect on platelet adhesion. The platelet adhesion was greatest on the positively charged surface, and decreased by almost 50% with the neutralization of this charge. A chitosan surface containing the negatively charged SO(3)(-) exhibited the fewest number of adherent platelets of all surfaces tested. Coagulation activation was not altered by the neutralization of the positive charge, but a marked increase of approximately 5-6 min in the plasma recalcification time (PRT) was displayed with the addition of the negatively charged species. Polyethylene (PE) surfaces were modified with the chitosan surfactant resulting in a significant improvement in blood compatibility, which correlated to the increasing PEG content within the polymer. Adsorption of the chitosan surfactants onto PE resulted in approximately an 85-96% decrease in the number of adherent platelets. The surfactant polymers also reduced surface induced coagulation activation, which was indicated by the PEG density dependent increase in PRTs. These results indicate that surface modification with our chitosan based surfactant polymers successfully improves blood compatibility. Moreover, the inclusion of either negatively charged SO(3)(-) groups or a high density of large water-soluble PEG side chains produces a surface that may be suitable for cardiovascular applications.     

Chemistry of aqueous silica nanoparticle surfaces and the mechanism of selective peptide adsorption

Control over selective recognition of biomolecules on inorganic nanoparticles is a major challenge for the synthesis of new catalysts, functional carriers for therapeutics, and assembly of renewable biobased materials. We found low sequence similarity among sequences of peptides strongly attracted to amorphous silica nanoparticles of various size (15-450 nm) using combinatorial phage display methods. Characterization of the surface by acid base titrations and zeta potential measurements revealed that the acidity of the silica particles increased with larger particle size, corresponding to between 5% and 20% ionization of silanol groups at pH 7. The wide range of surface ionization results in the attraction of increasingly basic peptides to increasingly acidic nanoparticles, along with major changes in the aqueous interfacial layer as seen in molecular dynamics simulation. We identified the mechanism of peptide adsorption using binding assays, zeta potential measurements, IR spectra, and molecular simulations of the purified peptides (without phage) in contact with uniformly sized silica particles. Positively charged peptides are strongly attracted to anionic silica surfaces by ion pairing of protonated N-termini, Lys side chains, and Arg side chains with negatively charged siloxide groups. Further, attraction of the peptides to the surface involves hydrogen bonds between polar groups in the peptide with silanol and siloxide groups on the silica surface, as well as ion-dipole, dipole-dipole, and van-der-Waals interactions. Electrostatic attraction between peptides and particle surfaces is supported by neutralization of zeta potentials, an inverse correlation between the required peptide concentration for measurable adsorption and the peptide pI, and proximity of cationic groups to the surface in the computation. The importance of hydrogen bonds and polar interactions is supported by adsorption of noncationic peptides containing Ser, His, and Asp residues, including the formation of multilayers. We also demonstrate tuning of interfacial interactions using mutant peptides with an excellent correlation between adsorption measurements, zeta potentials, computed adsorption energies, and the proposed binding mechanism. Follow-on questions about the relation between peptide adsorption on silica nanoparticles and mineralization of silica from peptide-stabilized precursors are raised.     

Combined depletion and electrostatic forces in polymer-induced membrane adhesion: a theoretical model

We develop a semi-quantitative analytical theory to describe adhesion between two identical planar charged surfaces embedded in a polymer-containing electrolyte solution. Polymer chains are uncharged and differ from the solvent by their lower dielectric permittivity. The solution mimics physiological fluids: It contains 0.1 M of monovalent ions and a small number of divalent cations that form tight bonds with the headgroups of charged lipids. The components have heterogeneous spatial distributions. The model was derived self-consistently by combining: (a) a Poisson-Boltzmann like equation for the charge densities, (b) a continuum mean-field theory for the polymer profile, (c) a solvation energy forcing the ions toward the polymer-poor regions, and (d) surface interactions of polymers and electrolytes. We validated the theory via extensive coarse-grained Molecular Dynamics (MD) simulations. The results confirm our analytical model and reveal interesting details not detected by the theory. At high surface charges, polymer chains are mainly excluded from the gap region, while the concentration of ions increases. The model shows a strong coupling between osmotic forces, surface potential and salting-out effects of the slightly polar polymer chains. It highlights some of the key differences in the behaviour of monomeric and polymeric mixed solvents and their responses to Coulomb interactions. Our main findings are: (a) the onset of long-ranged ion-induced polymer depletion force that increases with surface charge density and (b) a polymer-modified repulsive Coulomb force that increases with surface charge density. Overall, the system exhibits homeostatic behaviour, resulting in robustness against variations in the amount of charges. Applications and extensions of the model are briefly discussed.     

Water Structure and Phase Transition Near a Surface

This article concerns the density and orientational structure of water near uncharged (hydrophobic) and charged (hydrophilic) surfaces and the surface-induced phase transition accompanying a powerful, attractive surface force. Previously reported studies on the water structure using the integral equation theories and computer simulations are critically reviewed, and some controversial aspects are reexamined by additional calculations. Effects of the surface curvature and surface charge density on the structure are analyzed in detail. A new analysis is performed on the phase transition for water containing a trace amount of a hydrophobic component at surfaces. It is shown that the transition occurs even at charged surfaces and the attraction induced between surfaces predominates over the strong Coulombic repulsion at sufficiently small separations. At highly charged surfaces, however, the transition does not occur.     

Electrostatic surface interactions in mixtures of symmetric and asymmetric electrolytes: a Monte Carlo study

Canonical Monte Carlo simulations of the interaction between a uniformly charged spherical particle and a discretely charged planar surface in solutions of symmetric and asymmetric electrolytes were performed. To assess the nature of the interactions, the force exerted on the colloidal particle perpendicular to the planar surface was calculated. Attractive minima in the interaction force between the similarly charged surfaces reveal the occurrence of two phenomena: long-range attraction of electrostatic origin and short-range attraction due to depletion effects. The degree of electrostatic coupling determines the magnitude and range of like-charge attraction between the two surfaces.     

Long-range attraction between surfaces: existence and amplitude?

The possible occurrence of a long-range attraction between hydrophobic surfaces is fundamental for understanding the kinetics of protein folding or self-assembling structures, such as biological membranes, the stability of emulsions and inorganic dispersions. Direct force measurements have revealed the presence in water of a long-range attraction between macroscopic and hydrophobic surfaces. Nevertheless, the existence of this so-called ‘hydrophobic force’ for smaller objects is still under discussion. For macroscopic surfaces, it appears that electrostatic contributions due to surface heterogeneities and gas effects and/or bubble bridging have to be taken into account, but do they define an intrinsic ‘hydrophobic force’? For charged colloidal particles, theoretical predictions of electrostatic attraction and a phase separation when the counterions are multivalent, are partially confirmed by experiments, and recent experimental evidences of an attraction in a confined geometry put an exciting challenge for theoreticians.     

Simulating equilibrium surface forces in polymer solutions using a canonical grid method

A new simulation method for nonuniform polymer solutions between planar surfaces at full chemical equilibrium is described. The technique uses a grid of points in a two-dimensional thermodynamic space, labeled by surface area and surface separations. Free energy differences between these points are determined via Bennett's optimized rates method in the canonical ensemble. Subsequently, loci of constant chemical potential are determined within the grid via simple numerical interpolation. In this way, a series of free energy versus separation curves are determined for a number of different chemical potentials. The method is applied to the case of hard sphere polymers between attractive surfaces, and its veracity is confirmed via comparisons with established alternative simulation techniques, namely, the grand canonical ensemble and isotension ensemble methods. The former method is shown to fail when the degree of polymerization is too large. An interesting interplay between repulsive steric interactions and attractive bridging forces occurs as the surface attraction and bulk monomer density are varied. This behavior is further explored using polymer density functional theory, which is shown to be in good agreement with the simulations. Our results are also discussed in light of recent self-consistent field calculations which correct the original deGennes results for infinitely long polymers. In particular, we look at the role of chain ends by investigating the behavior of ring polymers.     

Phase diagram for the adsorption of weak polyelectrolytes at a soft charged surface

We have experimentally studied the adsorption of polyelectrolytes at oppositely charged surfaces. A weak flexible polyelectrolyte, poly(acrylic acid), was adsorbed from dilute solutions on a Langmuir film of a cationic amphiphile, dimethyldioctadecylammonium bromide. The polymer surface coverage, Gamma, at equilibrium was measured by two reflectivity techniques-ellipsometry and polarization modulated infrared reflection absorption spectroscopy (PM-IRRAS)-as a function of the surface charge density, sigma, and of the polymer ionization degree, alpha. Different adsorption regimes were evidenced. For weakly charged surfaces, sigma < sigma sat, Gamma increases with sigma and with 1/alpha, as expected for a neutralization of the surface by the adsorbed polymers. For highly charged surfaces, sigma > sigma sat, the adsorption of polyelectrolytes saturates. The mean orientation of the adsorbed chains also depends on the value of sigma: it is parallel to the surface for sigma < sigma (< sigma sat) and orthogonal to the surface for sigma > sigma. We have measured the values of sigma sat and sigma as a function of alpha and compared the results with existing theories.     

Modeling of confinement-induced phase transitions for surfactant layers on amphiphilic surfaces

A self-consistent field model is used to consider a solution of positively charged surfactants up to its critical micellization concentration adsorbing onto two surfaces in close proximity. Each surface mimics a polystyrene sulfonate interface; that is, hydrophobic properties are combined with a (fixed) negative charge. We observe large and sudden changes in adsorption as a function of separation, which are not normally considered when interpreting surface force measurements. The parameters are chosen such that the adsorbed surfactant layer is of a monolayer type when the surfaces are far apart. A typical interaction curve is presented for a fixed surfactant chemical potential, which is extracted from the set of adsorption isotherms each with a fixed slit width. When the slit width approaches the thickness of the two surfactant layers, a first-order phase transition takes place, which is driven by the unfavorable hydrophobic-water contacts. At the transition, the average orientation of the surfactants switches from a high concentration of tails at the surface to a bilayer configuration where tail profiles from both sides merge in the center. The headgroups are pulled slightly away from the surface. The interaction force jumps from a weak electrostatic repulsion at large distances (two effectively positively charged surface layers repel each other) to a strong electrostatic attraction at short distances (the central surfactant bilayer is attracted to the oppositely charged surfaces). The amount of adsorbed surfactants tend to decrease with decreasing distance between the surfaces but suddenly increases at the transition. Because of this, we anticipate that in surface force experiments, for example, there is a hysteresis associated with this transition: the forces and also the adsorbed amounts depend not only on the distance between the surfaces but also on the history if nonsufficient equilibration times are implemented.     

Adsorption of comb copolymers on weakly attractive solid surfaces

In this work continuum and lattice Monte Carlo simulation methods are used to study the adsorption of linear and comb polymers on flat surfaces. Selected polymer segments, located at the tips of the side chains in comb polymers or equally spaced along the linear polymers, are attracted to each other and to the surface via square-well potentials. The rest of the polymer segments are modeled as tangent hard spheres in the continuum model and as self-avoiding random walks in the lattice model. Results are presented in terms of segment-density profiles, distribution functions, and radii of gyration of the adsorbed polymers. At infinite dilution the presence of short side chains promotes the adsorption of polymers favoring both a decrease in the depletion-layer thickness and a spreading of the polymer molecule on the surface. The presence of long side chains favors the adsorption of polymers on the surface, but does not permit the spreading of the polymers. At finite concentration linear polymers and comb polymers with long side chains readily adsorb on the solid surface, while comb polymers with short side chains are unlikely to adsorb. The simple models of comb copolymers with short side chains used here show properties similar to those of associating polymers and of globular proteins in aqueous solutions, and can be used as a first approximation to investigate the mechanism of adsorption of proteins onto hydrophobic surfaces.     

Adsorption of Polyelectrolyte-Nanoparticle Systems on Silica: Influence on Interaction Forces

In this work, we have studied the interfacial properties of cationic polyelectrolyte (PE) and silica nanoparticle (NP) systems at macroscopic silica surfaces by means of ellipsometry. The influence of adsorbed layers on the interactions between silica surfaces was also investigated using the bimorph surface force apparatus. Added nanoparticles were observed to strongly swell the interfacial polyelectrolyte layers, an effect partly related to neutralization of charged polyelectrolyte groups. The effect was more pronounced for low charged than for highly charged polyelectrolytes. Overall, the presence of nanoparticles seemed to increase the repulsive interaction measured between silica surfaces. The force measured on approach was long range and quite strongly repulsive. On separation, an attractive bridging interaction was measured for polyelectrolyte-covered surfaces. For the low charged polyelectrolyte used in the study, the force turned repulsive on addition of nanoparticles. For the highly charged polyelectrolyte used, a change from a very strong attraction (involving a jump of the surfaces out of contact) to a very long-range elastic attractive force was observed on adding nanoparticles. The long-range elastic force indicates that polymer chains and nanoparticles form a transient network in the gap between the surfaces. The observed difference in the outward force curves may explain why the addition of nanoparticles appears to improve, e.g., shear-resistance and reflocculation characteristics of polymeric flocculants. Copyright 2000 Academic Press.     

Theory of polyelectrolyte adsorption onto surfaces patterned with charge and topography

Mean-field theory is used to derive criteria for the adsorption of a weakly charged polyelectrolyte molecule from salt solution onto surfaces patterned with charge and topography. For flat surfaces patterned with periodic arrays of charged patches, the adsorbed layer thickness predicted using mean-field theory and that found by Brownian dynamics simulations are in quantitative agreement in the strong-adsorption regime, which corresponds to sufficiently small kappa or sufficiently large |sigma(eff)q|, where kappa is the inverse Debye screening length, sigma(eff) is an effective surface charge density, and q is the charge on each segment of the polyelectrolyte. Qualitative agreement is obtained in the weak-adsorption regime, and for the case where surfaces are patterned with both charge and topography. For uniformly charged, sinusoidally corrugated surfaces, the theory predicts that the critical temperature required for adsorption can be greater than or less than the corresponding value for a flat surface depending on the relative values of kappa and the corrugation wave number. If the surface charge is also allowed to vary sinusoidally, then adsorption is predicted to occur only when the topography crests have a surface charge opposite to that of the polyelectrolyte. Surfaces patterned with rectangular indentations having charged bottoms which are separated by flat charged plateaus are investigated as well. Adsorption is predicted to occur even when the net surface charge is zero, provided that the plateaus have a charge opposite to that of the polyelectrolyte. If the charge on the plateaus and polyelectrolyte is the same, adsorption may still occur if electrostatic attraction from the indentation bottoms is sufficiently strong.     

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.     

Nanoscale actuation of electrokinetic flows on thermoreversible surfaces

We report on a novel approach for controlling nanohydrodynamic properties at the solid-liquid interfaces through the use of stimuli-responding polymer coatings. The end-tethered polymers undergo a phase separation upon external activation. The reversible change in the thickness and polarity of the grafted polymers yields in a dynamic control of the surface-generated, electrokinetic phenomena. Nonactivated, swollen polymers are thicker than the electrical double layer (EDL) and prohibit the development of an EOF even on charged surfaces. On the other hand, activated polymer chains shrink and become thinner than the EDL and allow for the EOF to build up unimpeded. We show here that, for given experimental conditions, the EOF velocity on the shrunken surface is 35 times greater than the one on the nonactivated surface. Furthermore, we reveal that coupling of such surfaces with dense arrays of thermal actuators developed in our laboratory can lead to novel micro- and nanofluidic devices.     

Adsorption of different amphiphilic molecules onto polystyrene latices

In order to know the influence of the surface characteristics and the chain properties on the adsorption of amphiphilic molecules onto polystyrene latex, a set of experiments to study the adsorption of ionic surfactants, nonionic surfactants and an amphiphilic synthetic peptide on different latex dispersions was performed. The adsorbed amount versus the equilibrium surfactant concentration was determined. The main adsorption mechanism was the hydrophobic attraction between the nonpolar tail of the molecule and the hydrophobic regions of the latex surface. This attraction overcame the electrostatic repulsion between chains and latex surface with identical charge sign. However, the electrostatic interactions chain-surface and chain-chain also played a role. General patterns for the adsorption of ionic chains on charged latex surfaces could be established. Regarding the shape, the isotherms presented different plateaus corresponding to electrostatic effects and conformational changes. The surfactant size also affects the adsorption results: the higher the hydrophilic moiety in the surfactant molecule the lower the adsorbed amount.