Brownian dynamics simulations of polyelectrolyte adsorption in shear flow: effects of solvent quality and charge patterning

We probe the effects of solvent quality and charge patterning on polyelectrolyte adsorption in shear flow using Brownian dynamics simulations with hydrodynamic interaction (HI). The polyelectrolyte is modeled as a freely jointed bead-rod chain, and electrostatic and non-electrostatic interactions are accounted for by using screened Coulombic and Lennard-Jones potentials, respectively. In the absence of flow, the conformation of a polyelectrolyte molecule adsorbed onto a uniformly charged surface changes from flat to globular with an increase in bead-bead attraction (hydrophobicity), consistent with prior experimental observations. In the presence of flow, migration due to bead-wall HI and, as a consequence, desorption decrease with an increase in bead-bead attraction, implying that flow-induced desorption is more difficult under poor-solvent conditions. When bead-bead non-electrostatic attraction is strong, desorption can be enhanced by increasing bead-bead electrostatic repulsion. Analogous to the effect of bead-surface electrostatic attraction, an increase in the strength of bead-surface non-electrostatic attraction reduces desorption. We also study the effect of shear flow on the adsorption of a polyelectrolyte molecule onto surfaces decorated with periodic arrays of charged patches. An increase in patch periodicity increases desorption even when the effective surface charge density is kept the same. The results of this work suggest mechanisms for controlling the desorption of polyelectrolyte molecules in shear flows.     

Brownian dynamics simulations of polyelectrolyte adsorption in shear flow with hydrodynamic interaction

The adsorption of single polyelectrolyte molecules in shear flow is studied using Brownian dynamics simulations with hydrodynamic interaction (HI). Simulations are performed with bead-rod and bead-spring chains, and electrostatic interactions are incorporated through a screened Coulombic potential with excluded volume accounted for by the repulsive part of a Lennard-Jones potential. A correction to the Rotne-Prager-Yamakawa tensor is derived that accounts for the presence of a planar wall. The simulations show that migration away from an uncharged wall, which is due to bead-wall HI, is enhanced by increases in the strength of flow and intrachain electrostatic repulsion, consistent with kinetic theory predictions. When the wall and polyelectrolyte are oppositely charged, chain behavior depends on the strength of electrostatic screening. For strong screening, chains get depleted from a region close to the wall and the thickness of this depletion layer scales as N(1/3)Wi(2/3) at high Wi, where N is the chain length and Wi is the Weissenberg number. At intermediate screening, bead-wall electrostatic attraction competes with bead-wall HI, and it is found that there is a critical Weissenberg number for desorption which scales as N(-1/2)kappa(-3)(l(B)|sigmaq|)(3/2), where kappa is the inverse screening length, l(B) is the Bjerrum length, sigma is the surface charge density, and q is the bead charge. When the screening is weak, adsorbed chains are observed to align in the vorticity direction at low shear rates due to the effects of repulsive intramolecular interactions. At higher shear rates, the chains align in the flow direction. The simulation method and results of this work are expected to be useful for a number of applications in biophysics and materials science in which polyelectrolyte adsorption plays a key role.     

Adsorption of charged dendrimers: a Brownian dynamics study

The adsorption of isolated charged dendrimers onto oppositely charged flat surfaces is studied in this work using Brownian dynamics simulations. The dendrimer is modeled as a freely jointed bead-rod chain in which excluded-volume interactions are modeled by a repulsive Lennard-Jones potential and bead-bead and bead-surface electrostatic interactions are described by screened Coulombic potentials. Adsorption behavior is studied as a function of inverse screening length, dendrimer generation, and dendrimer charge distribution. Adsorbed dendrimers adopt a disclike conformation in which they flatten in the direction normal to the surface and expand in the direction parallel to the surface. As the inverse screening length increases, the dendrimer expands in the normal direction and contracts in the parallel direction, adopting a conformation that is more stretched in the normal direction. When the inverse screening length becomes sufficiently large, the dendrimer desorbs and adopts a spherelike conformation. Bead density profiles show that adsorbed dendrimers form a two-layer structure, with one layer corresponding to adsorbed beads and a second, less dense layer corresponding to beads one rod length away from the surface. They also reveal how the distribution of monomers within the dendrimer and near the surface can be tailored by changing various problem parameters. The results presented here are expected to be helpful in providing qualitative guidance for dendrimer design in various applications.     

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.     

Brownian dynamics simulations of polyelectrolyte adsorption onto charged patterned surfaces

The adsorption of single polyelectrolyte molecules onto surfaces decorated with periodic arrays of charged patches was studied using Brownian dynamics simulations. A free-draining, freely jointed bead-rod chain was used to model the polyelectrolyte, and electrostatic interactions were incorporated using a screened Coulombic potential with the excluded volume accounted for by a hard-sphere potential. The simulations predicted that the polyelectrolyte lies close to the adsorbing surface if the patch length, surface charge density, and screening length are sufficiently large. Chain conformations were found to be very sensitive to patch length, patch spacing, and the nature of the charge on adjacent patches. This is due both to the size of the polymer relative to patch length and spacing and to the structure of the electric field near the surface. In some cases, the component of the radius of gyration parallel to the surface can be made smaller than its free-solution value, which is contrary to what is observed for a uniformly charged surface. Isolated charged patches were also considered, and significant adsorption was observed above a critical surface charge density. The results demonstrate how polyelectrolyte conformations can be controlled by the design of the charged patches and may be useful for applications in which adsorbed polyelectrolyte films play a key role.     

Strong and weak adsorptions of polyelectrolyte chains onto oppositely charged spheres

We investigate the complexation of long thin polyelectrolyte (PE) chains with oppositely charged spheres. In the limit of strong adsorption, when strongly charged PE chains adapt a definite wrapped conformation on the sphere surface, we analytically solve the linear Poisson-Boltzmann equation and calculate the electrostatic potential and the energy of the complex. We discuss some biological applications of the obtained results. For weak adsorption, when a flexible weakly charged PE chain is localized next to the sphere in solution, we solve the Edwards equation for PE conformations in the Hulthen potential, which is used as an approximation for the screened Debye-Huckel potential of the sphere. We predict the critical conditions for PE adsorption. We find that the critical sphere charge density exhibits a distinctively different dependence on the Debye screening length than for PE adsorption onto a flat surface. We compare our findings with experimental measurements on complexation of various PEs with oppositely charged colloidal particles. We also present some numerical results of the coupled Poisson-Boltzmann and self-consistent field equation for PE adsorption in an assembly of oppositely charged spheres.     

Brownian dynamics simulations of polyelectrolyte adsorption onto topographically patterned surfaces

The effect of patterned surface topography on the adsorption of single polyelectrolyte molecules is explored using Brownian dynamics simulations. The polyelectrolyte is modeled as a free-draining, freely jointed bead-rod chain, and electrostatic interactions are incorporated using a screened Coulombic potential with excluded volume interactions accounted for by the repulsive part of a Lennard-Jones potential. Topography consisting of periodically spaced valleys of square cross section separated by flat hills is considered. Chain conformations are characterized for a wide range of valley widths, depths, and spacings as well as for several different types of surface charge distributions. Depending on the parameter values describing the topography, the chains are found to adopt conformations ranging from flat and extended to those associated with bridge-, brush-, or semi-bridge-like structures. The formation of these structures is rationalized on the basis of a free-energy model that takes into account the increase in free energy due to entropic confinement, excluded volume interactions, and chain stretching as well as the decrease in free energy due to bead-surface electrostatic attraction. The results of this work are expected to be useful in designing patterned surface topography to control the conformations of adsorbed polyelectrolyte molecules.     

Universal consequences of the presence of excluded volume interactions in dilute polymer solutions undergoing shear flow

The role of solvent quality in determining the universal material properties of dilute polymer solutions undergoing steady simple shear flow is examined. A bead-spring chain representation of the polymer molecule is used, and the influence of solvent molecules on polymer conformations is modelled by a narrow Gaussian excluded volume potential that acts pairwise between the beads of the chain. Brownian dynamics simulations data, acquired for chains of finite length, and extrapolated to the limit of infinite chain length, are shown to be model independent. This feature of the narrow Gaussian potential, which leads to results identical to a delta-function repulsive potential, enables the prediction of both universal crossover scaling functions and asymptotic behavior in the excluded volume limit. Universal viscometric functions, obtained by this procedure, are found to exhibit increased shear thinning with increasing solvent quality. In the excluded volume limit, they are found to obey power law scaling with the characteristic shear rate beta, in close agreement with previously obtained renormalization group results. The presence of excluded volume interactions is also shown to lead to a weakening of the alignment of the polymer chain with the flow direction.     

Polyelectrolyte–particle complex formation. Polyelectrolyte linear charge density and ionic concentration effects. Monte Carlo simulations

The influence of the linear charge density (LCD) of a polyelectrolyte on its adsorption on an oppositely charged colloidal particle is investigated by Monte Carlo simulations. Adsorption characteristics are studied at different linear charge densities and ionic concentrations and for a given polyelectrolyte/particle size ratio so that particle curvature has full effect. The isolated polyelectrolyte goes through a smooth transition from a collapsed structure to an extended rod-like conformation with increasing the linear charge density in the low ionic concentration regime. In the high ionic concentration regime, the polyelectrolyte is less sensitive to the increase in the linear charge density and adopts a coil conformation. We found that complex formation is promoted by decreasing the ionic concentration and increasing the linear charge density and that large changes in the polymer dimensions are observed at the adsorption-desorption limit. By adjusting the linear charge density and ionic strength, we demonstrate that the adsorption-desorption limit corresponds to a sharp transition from non-adsorbed to adsorbed conformations and that the mean adsorption energy per monomer has to be less than -0.4 kT to achieve adsorption. We calculated that the linear charge density at the adsorption-desorption limit is related to the Debye-Hückel length according to LCD^crit ~3². At small values of the linear charge density and low ionic strength (no adsorption is observed at high ionic strength), a large amount of monomers are present in loops and tails. By increasing LCD, the amount of monomers in trains reaches a maximum value and the polyelectrolyte adopt flat conformation at the surface of the particle.     

Rheology of dilute polyelectrolyte solutions in narrow channels

The shear flow of dilute polyelectrolyte solutions bounded by either neutral or repulsive walls is modeled using a nonlinear dumbbell with conformation-dependent friction. Assuming that the configurational probability density function depends on the internal coordinates (r) and the distance of the center of mass of the molecule to the walls, coupled differential equations for the tensor moments <rr> are obtained. Coulombic repulsion between beads is considered to simulate the charge repulsion between ionized sites distributed along the backbone of a real polyelectrolyte. The repulsive interaction between the polyelectrolyte molecule and the charged walls is that of the DLVO model and the molecule is considered to be a charged sphere. Numerical solutions for the components of the tensor <rr> are worked out with the preaverage approach, and only when neutral walls considered are exact solutions obtained. Viscosity results show that in the limit of very wide channels, the corresponding viscosity in the bulk is obtained. The wall repulsion on the charged molecules produces migration of molecules towards the center of the channel resulting in a depleted layer with lower viscosity next to the walls. The calculated slip phenomenon using the method employed by Grisafi and Brunn is dependent on the beads repulsion and the shear rate. The slip velocity obtained with the Mooney method shows similarities with available experimental results for polyelectrolyte solutions. Birefringence calculations are performed in narrow and wide channels for different bead repulsions, with interesting results for both flexible and rigid molecules. Received: 26 September 1998 Accepted in revised form: 11 March 1999     

Molecular dynamics simulations of polyelectrolyte multilayering on a charged particle

Molecular dynamics simulations of polyelectrolyte multilayering on a charged spherical particle revealed that the sequential adsorption of oppositely charged flexible polyelectrolytes proceeds with surface charge reversal and highlighted electrostatic interactions as the major driving force of layer deposition. Far from being completely immobilized, multilayers feature a constant surge of chain intermixing during the deposition process, consistent with experimental observations of extensive interlayer mixing in these films. The formation of multilayers as well as the extent of layer intermixing depends on the degree of polymerization of the polyelectrolyte chains and the fraction of charge on its backbone. The presence of ionic pairs between oppositely charged macromolecules forming layers seems to play an important role in stabilizing the multilayer film.     

Polyelectrolyte adsorption: electrostatic mechanisms and nonmonotonic responses to salt addition

The main question addressed in this work is as follows: Under pure electrosorption conditions, that is, disregarding nonelectrostatic effects, how does the net adsorption of a polyelectrolyte at an oppositely charged surface respond to the addition of simple salt? Previous simulations and mean-field calculations have suggested that the polymers will desorb. However, we will demonstrate that an increased adsorption also is possible, even for pure electrosorption, at low and intermediate levels of salt. As this is a correlation-driven effect, mean field approaches will fail to capture it. Using simulations, one will in general need to simulate large systems and relatively long polymers. Also important is the presence of a proper bulk solution, with a finite and well-defined polyelectrolyte concentration. We have performed a theoretical study of polyelectrolyte adsorption, assuming screened Coulomb interactions between monomers; that is, the salt is implicit. This work focuses on the effects from ionic screening and polymer length. Specifically, the adsorption at a weakly charged colloidal particle, with a diameter of 200 nm, is monitored for various salt concentrations, in the presence of highly charged chains. Using simulations, we investigate polymers with two different degrees of polymerization: 40 and 160, respectively. These simulations are complemented by predictions from classical polymer density functional theory, utilizing a recently developed correlation-correction (Forsman, J.; Nordholm, S. Langmuir, in press). The agreement with corresponding simulations is semiquantitative, and because the calculations run many orders of magnitude faster than the simulations, longer and more realistic polymers could be studied with this approach. However, switching off the correlation-correction leads to a mean-field theory, which fails to even qualitatively reproduce the simulated response to screening.     

Nonrandom adsorption of polyelectrolyte chains on finite regularly charged surfaces

Adsorption phenomena are relevant in a wide variety of subjects, from biophysics to technological applications. Different aspects, such as molecular recognition, multilayer deposition, and dynamics of polymer adsorption have been addressed. The methodologies used range from analytical and numerical methods to molecular dynamics or Monte Carlo simulations. In this work, a coarse‐grained model is used to explore the adsorption of charged backbones to oppositely charged regions of a surface. These regions encompass those small enough to prevent complete adsorption, but extend to surfaces sufficiently large to promote adsorption with minimal effect on the three‐dimensional conformation in bulk. Apart from the different surface areas explored, variations on the surface charge density, polyelectrolyte chain length, and chain stiffness were also considered. The degree of compaction of the polyelectrolyte, on adsorption, is different from that found in the bulk. Also, results indicate an nonuniform adsorption pattern on regularly charged surfaces. © 2013 Wiley Periodicals, Inc.     

Packaging of a flexible polyelectrolyte inside a viral capsid: effect of salt concentration and salt valence

The effect of salt on the location and structure of a flexible polyelectrolyte confined inside a viral capsid and the Donnan equilibrium of the salt across the capsid have been examined using a coarse-grained model solved by Monte Carlo simulations. The polyelectrolyte was represented by a linear jointed chain of charged beads, and the capsid was represented by a spherical shell with embedded charges. At low salt concentration, the polyelectrolyte was strongly adsorbed onto the inner capsid surface, whereas at high salt concentration it was located preferentially in the central part of the capsid. Under the condition of equal Debye screening length, the electrostatic screening increased as the valence of the polyelectrolyte counterion was increased. The distribution of the small cations and anions was unequal across the capsid. An excess of polyelectrolyte counterions occurred inside the capsid, and the excess increased with the salt concentration. A simplified representation of the small ions through the use of the screened Coulomb potential provided only a qualitatively correct picture; the electrostatic screening originating from the small ions was exaggerated.     

Shear‐Induced Diffusion in Dilute Suspensions of Charged Colloids

We propose a model for the nonequilibrium enhancement of colloidal self‐diffusion in an externally imposed shear flow in charged systems. The diffusion enhancement is calculated in terms of electrostatic, two‐body interactions between the particles in shear flow. In the high‐shear rate, low‐volume fraction limit in which the model is valid, we compare these calculations to the experiments of Qiu et al. [PRL 61, 2554 (1988)] and simulations of Chakrabarti et al. [PRE 50, R3326 (1994)] and find good agreement on scaling and magnitude to within experimental uncertainty of the electrostatic parameters.     

Molecular dynamics simulations of polyelectrolyte-polyampholyte complexes. Effect of solvent quality and salt concentration

Complexation between polyelectrolyte and polyampholyte chains in poor solvent conditions for the polyelectrolyte backbone has been studied by molecular dynamics simulations. In a poor solvent a polyelectrolyte forms a necklace-like structure consisting of polymeric globules (beads) connected by strings of monomers. The simulation results can be explained by assuming the existence of two different mechanisms leading to the necklace formation. In the case of weak electrostatic interactions, the necklace formation is driven by optimization of short-range monomer-monomer attraction and electrostatic repulsion between charged monomers on the polymer backbone. In the case of strong electrostatic interactions, the necklace structure appears as a result of counterion condensation. While the short-range attractions between monomers are still important, the correlation-induced attraction between condensed counterions and charged monomers and electrostatic repulsion between uncompensated charges provide significant contribution to optimization of the necklace structure. Upon forming a complex with both random and diblock polyampholytes, a polyelectrolyte chain changes its necklace conformation by forming one huge bead. The collapse of the polyelectrolyte chain occurs due to the neutralization of the polyelectrolyte charge by polyampholytes. In the case of the random polyampholyte, the more positively charged sections of the chain mix with negatively charged polyelectrolyte forming the globular bead while more negatively charged chain sections form loops surrounding the collapsed core of the aggregate. In the case of diblock polyampholyte, the positively charged block, a part of the negatively charged block, and a polyelectrolyte chain form a core of the aggregate with a substantial section of the negatively charged block sticking out from the collapsed core of the aggregate. In both cases the core of the aggregate has a layered structure that is characterized by the variations in the excess of concentration of monomers belonging to polyampholyte and polyelectrolyte chains throughout the core radius. These structures appear as a result of optimization of the net electrostatic energy of the complex and short-range attractive interactions between monomers of the polyelectrolyte chain.     

Cascade of coil-globule conformational transitions of single flexible polyelectrolyte molecules in poor solvent

We show that hydrophobic flexible polyelectrolyte molecules of poly(2-vinylpyridine) and poly(methacryloyloxyethyl dimethylbenzylammonium chloride) are trapped and frozen due to adsorption on the mica surface, and the observed AFM single molecule structures reflect the molecular conformation in solution. An increase of the ionic strength of the solution induces the cascade of abrupt conformational transitions due to the intrachain segregation from elongated coil to compact globule conformations through intermediate pearl necklace-globule conformations with different amounts of beads per chain. The length of the necklaces and the number of beads decrease, while the diameter of beads increases with the increase of ionic strength. Coexistence at the same time of extended coils, necklaces with different amounts of beads, and compact globules indicates the cascade of the first-order-type phase transitions.     

Polyion adsorption onto catanionic surfaces. A Monte Carlo study

The adsorption of a single and negatively charged polyion with varying flexibility onto a surface carrying both negative and positive charges representing a charged membrane surface has been investigated by using a simple model employing Monte Carlo simulations. The polyion was represented by a sequence of negatively charged hard spheres connected with harmonic bonds. The charged surface groups were also represented by charged hard spheres, and they were positioned on a hard surface slightly protruding into the solution. The surface charges were either frozen in a liquidlike structure or laterally mobile. With a large excess of positive surface charges, the classical picture of a strongly adsorbed polyion with an extended and flat configuration emerged. However, adsorption also appeared at a net neutral surface or at a weakly negatively charged surface, and at these conditions the adsorption was stronger with a flexible polyion as compared to a semiflexible one, two features not appearing in simpler models containing homogeneously charged surfaces. The presence of charged surface patches (frozen surface charges) and the ability of polarization of the surface charges (mobile surface charges) are the main reasons for the enhanced adsorption. The stronger adsorption with the flexible chain is caused by its greater ability to spatially correlate with the surface charges.     

Structure of flexible and semiflexible polyelectrolyte chains in confined spaces of slit micro/nanochannels

Understanding the behavior of a polyelectrolyte in confined spaces has direct relevance in design and manipulation of microfluidic devices, as well as transport in living organisms. In this paper, a coarse-grained model of anionic semiflexible polyelectrolyte is applied, and its structure and dynamics are fully examined with Brownian dynamics (BD) simulations both in bulk solution and under confinement between two negatively charged parallel plates. The modeling is based on the nonlinear bead-spring discretization of a continuous chain with additional long-range electrostatic, Lennard-Jones, and hydrodynamic interactions between pairs of beads. The authors also consider the steric and electrostatic interactions between the bead and the confining wall. Relevant model parameters are determined from experimental rheology data on the anionic polysaccharide xanthan reported previously. For comparison, both flexible and semiflexible models are developed accompanying zero and finite intrinsic persistence lengths, respectively. The conformational changes of the polyelectrolyte chain induced by confinements and their dependence on the screening effect of the electrolyte solution are faithfully characterized with BD simulations. Depending on the intrinsic rigidity and the medium ionic strength, the polyelectrolyte can be classified as flexible, semiflexible, or rigid. Confined flexible and semiflexible chains exhibit a nonmonotonic variation in size, as measured by the radius of gyration and end-to-end distance, with changing slit width. For the semiflexible chain, this is coupled to the variations in long-range bond vector correlation. The rigid chain, realized at low ionic strength, does not have minima in size but exhibits a sigmoidal transition. The size of confined semiflexible and rigid polyelectrolytes can be well described by the wormlike chain model once the electrostatic effects are taken into account by the persistence length measured at long length scale.