Molecular dynamics simulations of end‐grafted polymers with charged side chains

We report molecular dynamics simulations on bottle‐brush polyelectrolytes end‐grafted to a planar surface. For each bottle‐brush polyelectrolyte, flexible charged side chains are anchored to one neutral main chain. The effects of the counterion valence and the grafting density on the density profiles and the structural characteristics of the brush were studied in this work. It is found that the electrostatic repulsion between charged monomers in the side chains leads an extended conformation of the brush in a solution containing monovalent counterions, while strong electrostatic binding of multivalent counterions to the side chains has a significant contribution to the collapse of the brush. For the trivalent case, the distribution of end monomers in the main chains becomes broader upon decreasing the grafting density, as compared with the monovalent case. However, the position of the distribution for the monovalent case is relatively insensitive to the change of the grafting density. Additionally, with increased counterion valence, enhanced electrostatic correlation between counterions and charged side chains also weakens the diffusive ability of counterions. © 2011 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys, 2011     

Effective charges of polyelectrolytes in a salt-free solution based on counterion chemical potential

The phenomenon of counterion condensation around a flexible polyelectrolyte chain with N monomers is investigated by Monte Carlo simulations in terms of the degree of ionization alpha, which is proportional to the effective charge. It is operationally defined as the ratio of observed to intrinsic counterion concentration, alpha = co/ci. The observed counterion concentration in the dilute polyelectrolyte solution is equivalent to an electrolyte solution of concentration co with the same counterion chemical potential. It can be determined directly by thermodynamic experiments such as ion-selective electrode. With the polyelectrolyte fixed at the center of the spherical Wigner-Seitz cell, the polymer conformation, counterion distribution, and chemical potential can be obtained. Our simulation shows that the degree of ionization rises as the polymer concentration decreases. This behavior is opposite to that calculated from the infinitely long charged rod model, which is often used to study counterion condensation. Moreover, we find that, for a specified line charge density, alpha decreases with an increment in chain length and chain flexibility. In fact, the degree of ionization is found to decline with increasing polymer fractal dimension, which can be tuned by varying bending modulus and solvent quality. Those results can be qualitatively explained by a simple model of two-phase approximation.     

Monte Carlo Simulation of Coil-to-Globule Transition of Compact Polymer Chains: Role of Monomer Interacting

Coil-to-globule transitions are fundamental problems existing in polymer science for several decades; however, some features are still unclear, such as the effect of chain monomer interaction. Herein, we use Monte Carlo simulation to study the coil-to-globule transition of simple compact polymer chains. We first consider the finite-size effects for a given monomer interaction, where the short chain exhibits a one-step collapse while long chains demonstrate a two-step collapse, indicated by the specific heat. More interestingly, with the decrease of chain monomer interaction, the critical temperatures marked by the peaks of heat capacity shift to low values. A closer examination from the energy, mean-squared radius of gyration and shape factor also suggests the lower temperature of coil-to-globule transition.     

Langevin dynamics of semiflexible polyelectrolytes: rod-toroid-globule-coil structures and counterion distribution

We have investigated the nature of counterion condensation on uniformly charged semiflexible polyelectrolyte chains and the concomitant configurations by monitoring the role of chain stiffness, chain length, counterion valency, and the strength of electrostatic interaction. The counterion condensation is seen to follow the adsorption process and the effective polymer charge increases with chain stiffness. Size and shape, as calculated through the radius of gyration, effective persistence length, and hydrodynamic radius, are studied. Stable coil-like, globular, folded-chain, toroidal, and rodlike configurations are possible at suitable combinations of values of chain stiffness, chain length, electrostatic interaction strength, and the valency of counterion. For high strengths of electrostatic interactions, sufficiently stiff polyelectrolytes form toroids in the presence of multivalent counterions, whereas flexible polyelectrolytes form disordered globules. The kinetic features of the nucleation and growth of toroids are monitored. Several metastable structures are found to frustrate the formation of toroids. The generic pathway involves the nucleation of one primary loop somewhere along the chain contour, followed by a growth process where the rest of the chain is folded continuously on top of the primary loop. The dependence of the average radii of toroids on the chain length is found to be roughly linear, in disagreement with existing scaling arguments.     

Architecture and solution properties of amphiphilic polymer brushes with peripheral charged ions

Functionalized poly(ethylene oxide) (PEO) macromonomers (alpha -tertiary amino and omega -methacryloyl groups) were prepared by ring-opening polymerization of live PEO anions end-capped with styrene oxide using 2-[2-(N ,N -dimethylamino)-ethoxy]ethanol potassium alkoxides as an initiator with methacryloyl chloride. PEO brushes were synthesized by free-radical homopolymerization of such PEO macromonomers. These brushes were converted into peripherally charged brushes by quaternization. We studied the solution properties of both types of brushes from the viewpoint of charge effect. It was found from dynamic light scattering (DLS) that the polymer brushes formed a single macromolecule in solution due to crowding of side chains. It was speculated from angular dependence measurements that the polymer brushes with large aspect ratios took a geometrical anisotropic conformation such as a cylinder. In methanol with a low dielectric constant, radius of gyration (R(G)), and cross-sectional radius of gyration (R(G,C)) of the polymer brushes with charged side chains were smaller than those of the polymer brushes without charges. In a solvent with a low dielectric constant such as methanol, ionic groups do not dissociate and condense. On the other hand, these physical values in an aqueous solution were somewhat larger than those of the polymer brushes without charges. In water with a high dielectric constant, peripherally charged brushes were strongly stabilized by electric double layers.     

Polymer chain conformation in the phase separation process of a binary liquid mixture

A polymer chain conformation change near the critical point of liquid-liquid phase separation was investigated. Poly(N-isopropylacrylamide) labeled with a small amount of carbazolyl group for a fluorophore (P(NIPA-Cz)) was prepared. A ternary system of P(NIPA-Cz)+cyclohexane+methanol was investigated by the fluorescence spectroscopic technique. A mixed solvent of cyclohexane+methanol (CH/MeOH) shows phase separation at the upper critical solution temperature. Light scattering intensity, fluorescence emission intensity and fluorescence anisotropy ratio, as a function of temperature, were measured with quasi statically approaching to the critical demixing point. The fluorescence intensity of the carbazolyl groups attached to the polymer chain decreases with approaching to the critical temperature. This result suggests that the radius of gyration of the polymer decreases upon approaching to the critical demixing point of the solvent. We discuss the collapse and aggregation processes of the polymer based on the fluorescence quenching method. The rotational diffusion coefficient of carbazolyl groups attached to the polymer chain was estimated by the fluorescence depolarization technique. The rotational motion of carbazolyl groups is slowed down upon approaching the critical point.     

端基附壁高分子模型链的塌缩转变

借MonteCarlo和模拟退火方法研究了端基附壁高分子尾形链在不同温度下的形态变化 ,链的模型采用有最近邻相互吸引作用的自避行走 .计算机实验结果表明 ,对于端基附壁的高分子尾形链 ,与自由链一样 ,当温度逐渐降低时高分子链会发生从松散的无规线团到紧密球体的塌缩转变 .计算机模拟得到了端基附壁链的均方末端距及其分量 ,均方回转半径及分布随温度的变化 .由于刚性壁的影响 ,使得有限链长的高分子尾形链与自由链相比 ,其表示链尺寸 温度关系的曲线要稍低 .模拟还发现 ,在高温时壁对链形状的影响比较大 ,壁垂直方向上尺寸明显大于平行方向的尺寸 ,后者接近于自由链的尺寸分量 .然而 ,低于θ温度时 ,尤其是完全塌缩之后 ,壁对链形状的影响已经很小 .     

Encapsulation of a polyelectrolyte chain by an oppositely charged spherical surface

Using the ground state dominance approximation and a variational theory, we study the encapsulation of a polyelectrolyte chain by an oppositely charged spherical surface. The electrostatic attraction between the polyelectrolyte and the surface and the entropy loss of the encapsulated polyelectrolyte chain dictate the optimum conditions for encapsulation. Two scenarios of encapsulation are identified: entropy-dominated and adsorption-dominated encapsulation. In the entropy-dominated encapsulation regime, the polyelectrolyte chain is delocalized, and the optimum radius of the encapsulating sphere decreases with increasing the attraction. In the adsorption-dominated encapsulation regime, the polyelectrolyte chain is strongly localized near the surface, and the optimum radius increases with increasing the attraction. After identifying a universal encapsulation parameter, the dependencies of the optimum radius on the salt concentration, surface charge density, polymer charge density, and polymer length are explored.     

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.     

Counterion condensation on charged spheres, cylinders, and planes

We use the framework of counterion condensation theory, in which deviations from linear electrostatics are ascribed to charge renormalization caused by collapse of counterions from the ion atmosphere, to explore the possibility of condensation on charged spheres, cylinders, and planes immersed in dilute solutions of simple salt. In the limit of zero concentration of salt, we obtain Zimm-Le Bret behavior: a sphere condenses none of its counterions regardless of surface charge density, a cylinder with charge density above a threshold value condenses a fraction of its counterions, and a plane of any charge density condenses all of its counterions. The response in dilute but nonzero salt concentrations is different. Spheres, cylinders, and planes all exhibit critical surface charge densities separating a regime of counterion condensation from states with no condensed counterions. The critical charge densities depend on salt concentration, except for the case of a thin cylinder, which exhibits the invariant criticality familiar from polyelectrolyte theory.     

Behavior of a polymer chain immersed in a binary mixture of solvents

The behavior of a polymer chain immersed in a binary solvent mixture is investigated via a single-polymer simulation using an effective Hamiltonian, where the solvent effects are taken into account through a density-functional theory for polymer-solvent admixtures. The liquid-liquid phase separation of the binary solvent mixture is modeled as that of a Lennard-Jones binary fluid mixture with weakly attractive interactions between the different components. Two types of energetic preferences of the polymer chain for the better solvent-(A) no preferential solvophilicity and (B) strong preferential solvophilicity-are employed as polymer-solvent interaction models. The radius of gyration and the polymer-solvent radial distribution functions are determined from the simulations of various molar fractions along an isotherm slightly above the critical temperature of the liquid-liquid phase separation. These quantities near the critical point conspicuously depend on the strength of the preferential solvophilicity. In the case where the polymer exhibits a strong preferential solvophilicity, a remarkable expansion of the polymer chain is observed near the critical point. On the other hand, in the case where the polymer has no preferential solvophilicity, no characteristic variation of the polymer conformation is observed even near the critical point. These results indicate that the expansion of a polymer chain enhances the local phase separation around it, acting as a nucleus of demixing in the vicinity of the critical point. This phenomenon in binary solvents near the liquid-liquid critical point is similar to the expansion of the polymer chain in one-component supercritical solvents near the liquid-vapor critical point, which we have reported [T. Sumi and H. Sekino J. Chem. Phys. 122, 194910 (2005)].     

Collapse kinetics of vibrated granular chains

The kinetics of the collapse of the coil state into condensed states is studied with vibrated granular chain composed of N metal beads partially immersed in water. The radius of gyration of the chain, R(g) is measured. For short chains (N < 140), disk-like condensed state is formed and R(g) decreases with time such that the function ΔR(g)(2) (≡ R(g)(2) - R(g)(2)(∞)) = A e(-t/τ), where the relaxation time τ follows a power-law dependence on the chain length N with an exponent γ = 1.9 ± 0.2. For the chains with length N ≥ 300, rod-like clusters are observed during the initial stage of collapse and R(g)(2) = R(g)(2)(0) - Bt(β), with β = 0.6 ± 0.1. In the coarsening stage, the exponential dependence of ΔR(g)(2) on time still holds, however, the relaxation time τ fluctuates and has no simple dependence on N. Furthermore, the time dependence of the averaged radius of gyration of the individual clusters, R(g,cl) can be described by the theory of Lifshitz and Slyozov. A peak in the structure function of long chains is observed in the initial stage of the collapse transition. The collapse transition in the bead chains is a first order phase transition. However, features of the spinodal decomposition are also observed.     

Counterion adsorption theory of dilute polyelectrolyte solutions: apparent molecular weight, second virial coefficient, and intermolecular structure factor

Polyelectrolyte chains are well known to be strongly correlated even in extremely dilute solutions in the absence of additional strong electrolytes. Such correlations result in severe difficulties in interpreting light scattering measurements in the determination of the molecular weight, radius of gyration, and the second virial coefficient of charged macromolecules at lower ionic strengths from added strong electrolytes. By accounting for charge-regularization of the polyelectrolyte by the counterions, we present a theory of the apparent molecular weight, second virial coefficient, and the intermolecular structure factor in dilute polyelectrolyte solutions in terms of concentrations of the polymer and the added strong electrolyte. The counterion adsorption of the polyelectrolyte chains to differing levels at different concentrations of the strong electrolyte can lead to even an order of magnitude discrepancy in the molecular weight inferred from light scattering measurements. Based on counterion-mediated charge regularization, the second virial coefficient of the polyelectrolyte and the interchain structure factor are derived self-consistently. The effect of the interchain correlations, dominating at lower salt concentrations, on the inference of the radius of gyration and on molecular weight is derived. Conditions for the onset of nonmonotonic scattering wave vector dependence of scattered intensity upon lowering the electrolyte concentration and interpretation of the apparent radius of gyration are derived in terms of the counterion adsorption mechanism.     

Development of knotting during the collapse transition of polymers

A dynamic Monte Carlo simulation of the collapse transition of polymer chains is presented. The chains are represented as self-avoiding walks on the simple cubic lattice with a nearest-neighbor contact potential to model the effect of solvent quality. The knot state of the chains is determined using the knot group procedure presented in the accompanying paper. The equilibrium knot spectrum and the equilibrium rms radius of gyration as functions of the chain length and the contact potential are reported. The collapse transition was studied following quenches from good-to poor-solvent conditions. Our results confirm the prediction that the newly formed globule is not yet at equilibrium, since it has not yet achieved its equilibrium knot spectrum. For our model system, the relaxation of the knot spectrum is about an order of magnitude slower than that of the radius of gyration. The collapse transition is also studied for a model in which both ends of the chain remain in good-solvent conditions. Over the time scale of these simulations, knot formation is frustrated in this inhomogeneous model, verifying that the mechanism of knotting is the tunneling of chain ends in and out of the globule.     

Dynamics of liquid–liquid demixing: mesoscopic simulations of polymer solutions

In this article, the mesoscopic simulation method dissipative particle dynamics (DPD) is applied to study the dynamics of polymer–solvent liquid–liquid phase separation. It will be shown that the degree of branching has a pronounced effect on the radius of gyration and the centre of mass diffusion of the polymer. Based on the simulation results it can be concluded that the difference in chemical potential between the mixed and the demixed state is the main driving force behind the centre of mass diffusion (and thus phase separation), rather than the reduced radius of gyration due to to polymer chain collapse.     

Collapse of a two-state sticky hard sphere chain

The collapse of a homopolymer gaussian chain into a globule is represented as a transition between two states, viz., extended and collapsed. Appropriately, this model has been labeled as the all-or-none view of chain collapse. In the collapsed state, the single polymer partition function is expressed by a single Mayer diagram with the maximum number of f-bonds arising from nonbonded square well interactions. Our target is the dependence of the transition temperature on chain length and the interaction range of the square well, as indicated through the behavior of the radius of gyration and the constant volume heat capacity. Properties of the collapse transition are calculated exactly for chains with three to six backbone atoms and heuristically for long chains using arguments derived from the small chains and from conditions of integrability. Comparison with simulation studies is made.     

Adsorption of weakly charged polyelectrolytes onto oppositely charged spherical colloids

The adsorption of a long weakly charged flexible polyelectrolyte in a salt solution onto an oppositely charged spherical surface is investigated. An analytical solution for Green's function is derived, which is valid for any sphere radius and consistently recovers the result of a planar surface in the limit of large sphere radii, by substituting the Debye-Hückel potential via the Hulthén potential. Expressions for critical quantities like the critical radius and the critical surface charge density are provided. In particular, we find a universal critical line for the sphere radius as a function of the screening length separating adsorbed from desorbed states. Moreover, results for the monomer density distribution, adsorbed layer thickness, and the radius of gyration are presented. A comparison of our theoretical results with experiments and computer simulations yields remarkably good agreement.     

Structure of polyelectrolyte complexes by Brownian dynamics simulation: effects of the bond length asymmetry of the polyelectrolytes

Brownian dynamics simulations were performed to study the structure of polyelectrolyte complexes formed by two flexible, oppositely charged polyelectrolyte chains. The distribution of monomers in the complex as well as the radius of gyration and structure factor of complexes and individual polyelectrolytes are reported. These structural properties were calculated for polyelectrolyte chains with equal number of monomers, keeping constant the bond length of the negative chain and increasing the bond length of the positive chain. This introduces an asymmetry in the length of the chains that modulates the final structure of the complexes. In the symmetric case the distribution of positive and negative monomers in the complex are identical, producing clusters that are locally and globally neutral. Deviations from the symmetric case lead to nonuniform, asymmetric monomer distributions, producing net charge oscillations inside the complex and large changes in the radius of gyration of the complex and individual chains. From the radius of gyration of the polyelectrolyte chains it is shown that the positive chain is much more folded than the negative chain when the chains are asymmetric, which is also confirmed through the scaling behavior of the structure factors.     

Monte Carlo simulations of charged dendrimer-linear polyelectrolyte complexes and explicit counterions

We study complexes composed of one dendrimer of generation G = 4 (G4 dendrimer) with N(t) = 32 charged terminal groups and an oppositely charged linear polyelectrolyte accompanied by neutralizing counterions in an athermal solvent using Monte Carlo simulations based on the bond fluctuation model. In our study both the full Coulomb potential and the excluded volume interactions are taken into account explicitly with the reduced temperature τ and the chain length N(ch) as the main simulation parameters. Our calculations indicate that there exist three temperature ranges that determine the behavior of such complexes. At τ(complex) stable charged dendrimer-linear polyelectrolyte complexes are formed first, which are subsequently accompanied by selective counterion localization within the complex interior at τ(loc) ≤ τ(complex), and counterion condensation as temperature is further decreased below τ(cond) < τ(loc). In particular, we observe that condensation takes place exclusively on the excess charges in the complex and thus no condensation is observed at the compensation point (N(ch) = N(t)), irrespective of τ. For N(ch) ≠ N(t) the complex is overally charged. Furthermore, we discuss the size and structure of the dendrimer and the linear polyelectrolyte within the complex, as well as spatial distributions of monomers and counterions. Conformations of the chain in the bound state are analysed in terms of loops, trains, and tails.     

Simulation Study of the Conformational Properties of Diblock Polyelectrolytes in Salt Solutions

Coarse-grained molecular dynamics simulations are performed to understand the behavior of diblock polyelectrolytes in solutions of divalent salt by studying the conformations of chains over a wide range of salt concentrations. The polymer molecules are modeled as bead spring chains with different charged fractions and the counterions and salt ions are incorporated explicitly. Upon addition of a divalent salt, the salt cations replace the monovalent counterions, and the condensation of divalent salt cations onto the polyelectrolyte increases, and the chains favor to collapse. The condensation of ions changes with the salt concentration and depends on the charged fraction. Also, the degree of collapse at a given salt concentration changes with the increasing valency of the counterion due to the bridging effect. As a quantitative measure of the distribution of counterions around the polyelectrolyte chain, we study the radial distribution function between monomers on different polyelectrolytes and the counterions inside the counterion worm surrounding a polymer chain at different concentrations of the divalent salt. Our simulation results show a strong dependence of salt concentration on the conformational properties of diblock copolymers and indicate that it can tune the self-assembly behaviors of such charged polyelectrolyte block copolymers.