Globules of annealed amphiphilic copolymers: Surface structure and interactions

A mean-field theory of globules of random amphiphilic copolymers in selective solvents is developed for the case of an annealed copolymer sequence: each unit can be in one of two states, H (insoluble) or P (soluble or less insoluble). The study is focussed on the regime when H and P units tend to form long blocks, and when P units dominate in the dilute phase, but are rare in the globule core. A first-order coil-to-globule transition is predicted at some T = T _cg. The globule core density at the transition point increases as the affinity of P units to the solvent, ˜, is increased. Two collapse transitions, coil → “loose” globule and “loose” globule → “dense” globule, are predicted if ˜ is high enough and P units are marginally soluble or weakly insoluble. H and P concentration profiles near the globule surface are obtained and analyzed in detail. It is shown that the surface excess of P units rises as ˜ is increased. The surface tension decreases in parallel. Considering the interaction between close enough surfaces of two globules, we show that they always attract each other at a complete equilibrium. It is pointed out, however, that such equilibrium may be difficult to reach, so that partially equilibrium structures (defined by the condition that a chain forming one globule does not penetrate into the core of the other globule) are relevant. It is shown that at such partial equilibrium the interaction is repulsive, so the globules may be stabilized from aggregation. The strongest repulsion is predicted at the coil-to-globule transition point T _cg: the repulsion force decreases with the distance between the surfaces according to a power law. In the general case (apart from T _cg) the force vs. distance decay becomes exponential; the decay length ξ diverges as T → T _cg. The developed theory explains certain anomalous properties observed for globules of amphiphilic homopolymers.     

Charge relaxation in polyampholytes of various statistics

We discuss theoretically the relaxation of charge fluctuations in polyampholyte solutions. It has been shown previously by some of us (J. Wittmer et al. Europhys. Lett. 24, 263 (1993)) that the charge distribution along the polyampholyte backbone has a dramatic influence on the polarization energy and hence on the solubility. Here it is demonstrated that a similar effect exists for the charge relaxation. The charge relaxation mechanism qualitatively depends on the statistics: for alternating polyampholytes the relaxation is mainly due to local dipole inversion and is not primarily driven by electrostatic interactions, whereas for random polyampholytes it is driven by electrostatic interactions. Intermediate statistics (with short-ranged (exponential) correlations) appear as a combination of these two limiting cases: short-wavelength modes are insensitive to the loss of correlations along the backbone, whereas long-wavelength modes correspond to a random statistics with renormalized charges. The relaxation of the dielectric constant is also calculated. Received: 20 December 2002 / Accepted: 13 March 2003 / Published online: 24 April 2003 RID="a" ID="a"e-mail: johner@ics.u-strasbg.fr     

Phase behavior of a suspension of hard spherocylinders plus ideal polymer chains

We study isotropic-isotropic and isotropic-nematic phase transitions of fluid mixtures containing hard spherocylinders (HSC) and added non-adsorbing ideal polymer chains using scaled particle theory (SPT). First, we investigate isotropic-nematic (I -N phase coexistence using SPT in the absence of polymer. We compare the results obtained using a Gaussian form of the orientational distribution function (ODF) to minimize the free energy versus minimizing numerically. We find that formal numerical minimization gives results that are much closer to computer simulation results. In order to describe mixtures of HSC plus ideal chains we studied the depletion of ideal chains around a HSC. We analyze the density profiles of ideal chains near a hard cylinder and find the depletion thickness δ is a function of the ratio of the polymer's radius of gyration R_g and the cylinder radius R_c. Our results are compared with a common approximation in which the depletion thickness is taken equal to the radius of gyration of the polymer chain. We incorporate the correct depletion thickness into SPT and find that for R _g/R _c < 1.56 using ideal chains gives phase transitions at smaller polymer concentrations, whereas for R _g/R _c > 1.56 , which is a common experimental situation, the phase transitions are found at larger polymer concentrations with respect to δ = R _g . The differences are significant, especially for R _g ≫ R _c , so we can conclude it is essential to take into account the properties of ideal polymer chains and the resulting depletion near a cylinder. Finally, we present phase diagrams for rod-polymer mixtures which could be realized under experimental conditions.     

Spherical polyelectrolyte block copolymer micelles: Structural change in presence of monovalent salt

Spherical polyelectrolyte block copolymer micelles were investigated as a function of added NaCl salt concentration using Small-Angle Neutron Scattering (SANS) and Light Scattering (LS). The micelles are formed by the self-association of charged-neutral copolymers made of a long deuterated polyelectrolyte moiety (NaPSS_d)₂₅₁ and a short hydrophobic moiety (PEP)₅₂. In presence of salt, the core shape and the aggregation number of the micelles are not affected. The hydrodynamic radius of the micelle is found to be identical to the radius of the whole micelle deduced from neutron scattering and thus the hydrodynamic radius is a valid measure of the corona thickness. At the lowest salt concentrations investigated the thickness of the corona, R_s, remains essentially constant and a contraction is observed above an added-salt concentration c_s of 2×10^-2 M where this crossover concentration corresponds to the average ionic strength of the free counterions in the corona. The contraction takes place while maintaining a rod-like behavior of the chains at short scale and obeys to: R_s c_s^-0.18. The exponent 0.18 suggests an electrostatic persistence length proportional to the Debye screening length.     

Conformation of poly(L-lysine)-graft-poly(ethylene glycol) molecular brushes in aqueous solution studied by small-angle neutron scattering

Small-angle neutron scattering (SANS) has been employed for the analysis of conformations of poly(L-lysine)-graft-poly(ethylene glycol) (PLL-g -PEG) molecular bottle brushes in aqueous solutions. The degree of polymerisation of the PEG chains was systematically varied in order to unravel dependence of the conformational properties of the bottle brushes on the molecular weight of the grafted chains. The grafting density was kept constant and high enough to ensure strong overlap of the PEG chains. The scattering spectra were fitted on the basis of the model of an effective worm-like chain with the account of average radial distribution and local fluctuations of the PEG density in the bottle brush. The results of the fits indicate that molecular brushes retain weakly bent configuration on the length scale of the order of (or larger than) the brush thickness. This finding is in agreement with earlier simulation and recent theoretical results.     

Curved polymer and polyelectrolyte brushes beyond the Daoud-Cotton model

We revise the classical Daoud-Cotton (DC) model to describe conformations of polymer and polyelectrolyte chains end-grafted to convex spherical and cylindrical surfaces. In the framework of the DC model, local stretching of chains in the brush does not depend on the degree of polymerization of grafted chains, and the polymer density profile follows a single-exponent power law. This model, however, does not correspond to a minimum in free energy of the curved brush. The nonlocal (NL) approximation exploited in the present paper implies the minimization of the overall free energy of the brush and predicts that the polymer density profile does not follow a single-exponent power law. In the limit of large surface curvature the NL approximation provides the same scaling laws for brush thickness and free energy as the local DC model. Numerical prefactors are however different. Extra extension of chains in the brush interior region leads to larger equilibrium brush thickness and lower free energy per chain. A significant difference between outcomes of the two models is found for brushes formed by ionic polymers, particularly for weakly dissociating (p H-sensitive) polyelectrolytes at low solution salinity.     

Nematic ordering in anisotropic elastomers: Effect of frozen anisotropy

Nematic ordering in anisotropic non-Gaussian elastomers is considered theoretically using mean field approximation. We focus on the effect of anisotropy during network cross-linking on the system elasticity and, in particular, on the so-called soft deformation mode. As the main result, we calculate the dependence of the elastomer free energy on the angle between the axis of “frozen” anisotropy and the nematic director. The dependence of the isotropic-nematic transition point on the orientational field acting on the monomers during the cross-linking process is also calculated. Received: 5 November 1997 / Revised and Accepted: 29 June 1998     

Melt fracture of entangled polymers*

We construct a model for a slippage plane in a sheared melt, based on a balance between reptation bridging and shear debonding. The resulting state could show up at rather low shear rates and be locally stable. But it is not easy to nucleate: the conventional entangled state is also locally stable. We propose that slippage occurs on solid walls: either at the container surface, or on dust particles floating in the melt. Slippage at solid/melt interfaces was studied (experimentally and theoretically) long ago. There is a critical stress for slippage: our estimate (for strong adsorption of melt chains on the solid) gives (plateau modulus) for typical cases. Thus, melt fracture is expected at moderate stresses, in agreement with observations by S.Q. Wang and coworkers.     

Micelles of poly(isoprene-b-2-vinylpyridine-b-ethylene oxide) terpolymers in aqueous media and their interaction with surfactants

Well-defined poly(isoprene-b-2-vinylpyridine-b-ethylene oxide) (PI-P2VP-PEO) triblock terpolymers were synthesized by anionic polymerization high-vacuum techniques. The terpolymers formed spherical three-layer (onion-type) micelles in neutral and acidic pH aqueous media as evidenced by static and dynamic light scattering. In pure water, kinetically frozen micelles with a core composed of a soft PI inner part and a hard P2VP outer shell and protected by a neutral PEO corona were formed. In acidic media the core was formed by the soft PI hydrophobic segment, whereas the corona consisted of an inner cationic polyelectrolyte P2VPH⁺ part and an outer PEO shell. The aggregation numbers were found to be high in all cases, due to the high hydrophobicity of the core-forming blocks. In the latter case an increase in size was observed due to the electrostatic repulsions between the P2VPH⁺ chains in the inner part of the corona, which is also responsible for the lower aggregation numbers observed in the acidic solutions. The interaction of these onion-type micelles with cationic (DTMAB) and anionic (SDS) surfactants led to the formation of mixed polymer/surfactant aggregates. Their structural characteristics could be varied by combining changes in surfactant type and concentration, solution pH and type of electrostatic interaction, leading to interesting, block-copolymer-based, environmentally responsive colloidal systems.     

Effects of ionic strength and charge annealing in star-branched polyelectrolytes

We have developed a scaling theory that describes the conformations of weak star-branched polyelectrolytes in dilute solutions. The dependences of the overall star size on the number of branches and on the ionic strength of the solution (tuned by the addition of low molecular weight salt) are analyzed. The intrinsic structure of the polyelectrolyte stars in salt-free and salt-added solutions is discussed in terms of concentration and elastic blobs. In contrast to neutral stars, the swollen corona of the polyelectrolyte star consists of blobs which are not closely packed. We have shown that the size of star polyelectrolytes is less sensitive to the variation in the ionic strength than the size of linear polyelectrolytes. The effects of the ionization-recombination balance in the star polyelectrolyte were taken into account. For polyelectrolytes with small ionization constant, the size of the star depends nonmonotonically on the number of branches and on the ionic strength of the solution due to recombination of counterions with charged monomers. Received: 10 November 1997 / Revised: 16 February 1998 / Accepted: 1st April 1998     

Excluded-volume polymer-induced depletion interaction between parallel flat plates

The interaction between two parallel plates due to non-adsorbing polymer chains with excluded volume is calculated using the adsorption method. The adsorption is calculated from the profile of the polymer segment concentration between the plates, which is obtained from the product function of the concentration profile near a single wall, involving the correlation length. The renormalization group theory provides expressions for the osmotic pressure and consequently for the osmotic compressibility, chemical potential and correlation length of a polymer solution. Both the local polymer concentration profiles as well as the minimum of the interaction potential between the plates agree with recently published self-avoiding random walk computer simulations. Received 9 August 2001     

Dielectrophoresis of nanocolloids: A molecular dynamics study

Dielectrophoresis (DEP), the motion of polarizable particles in non-uniform electric fields, has become an important tool for the transport, separation, and characterization of microparticles in biomedical and nanoelectronics research. In this article we present, to our knowledge, the first molecular dynamics simulations of DEP of nanometer-sized colloidal particles. We introduce a simplified model for a polarizable nanoparticle, consisting of a large charged macroion and oppositely charged microions, in an explicit solvent. The model is then used to study DEP motion of the particle at different combinations of temperature and electric field strength. In accord with linear response theory, the particle drift velocities are shown to be proportional to the DEP force. Analysis of the colloid DEP mobility shows a clear time dependence, demonstrating the variation of friction under non-equilibrium. The time dependence of the mobility further results in an apparent weak variation of the DEP displacements with temperature.     

Stratification of foam films made from polyelectrolyte solutions

Studies of thin liquid films, made from semidilute polyelectrolyte solutions, are presented. The disjoining pressure variation with film thickness exhibits oscillations, corresponding to film stratification. The oscillations become sharper as the polymer concentration c increases, and disappear when salt is added. The period of the oscillations scales as c ^-1/2. The observed stratification is related to the polymer network and the size of the steps to the mesh size ξ. Received 25 April 2000 and Received in final form 3 October 2000     

Kinetics of phase separation in ramified polymer blends of arbitrary topology

We present here a theoretical study of kinetics of phase separation within a mixture made of two chemically incompatible ramified polymers. For simplicity, we assume that they have the same topology. We are interested in the variation of the relaxation rate, _q, versus the wave number q, in the vicinity of the spinodal temperature. The kinetics is governed by local (Rouse) and reptation motions (faster and slower modes). For qR_G 1 (R_G being the gyration radius), kinetics is entirely controlled by local motions where each chain moves inside its own tube, and we show that the corresponding characteristic frequency, ^{-1}_q, scales as ^{-1}_{q G}q⁶, where _G is a known topological factor. For qR_G 1, however, kinetics is rather dominated by long-wavelength (reptation) motions where unlike ramified polymers creep inside a long tube. For this case, we find that ^{-1}_q ( 0 )q² (_c - ), where ( 0 ) is another known topological factor that represents the total mobility of free monomers belonging to connected chains and reticulation points, and _c accounts for the critical value of the segregation parameter. Finally, the derived relaxation rate must be compared to that relative to a linear polymer mixture.     

Dynamics of strongly entangled polymer systems: activated reptation

The effect of excluded-volume interactions on the reptation dynamics of long polymer chains is considered theoretically. It is shown that interactions give rise to an exponential increase of the reptation time, , if polymer chains are long enough: , where is the number of monomers per entanglement. We propose a novel dynamical mechanism of activated reptation implying that neighboring chains exchange conformations of their terminal fragments. It is shown that the exchange mechanism is compatible with the equilibrium polymer chain statistics and that it provides a bridge between the previous theories. Received: 25 July 1997 / Accepted: 8 October 1997     

Langevin dynamics of the glass forming polymer melt: Fluctuations around the random phase approximation

In this paper the Martin-Siggia-Rose (MSR) functional integral representation is used for the study of the Langevin dynamics of a polymer melt in terms of collective variables: mass density and response field density. The resulting generating functional (GF) takes into account fluctuations around the random phase approximation (RPA) up to an arbitrary order. The set of equations for the correlation and response functions is derived. It is generally shown that for cases whenever the fluctuation-dissipation theorem (FDT) holds we arrive at equations similar to those derived by Mori-Zwanzig. The case when FDT in the glassy phase is violated is also qualitatively considered and it is shown that this results in a smearing out of the ideal glass transition. The memory kernel is specified for the ideal glass transition as a sum of all “water-melon” diagrams. For the Gaussian chain model the explicit expression for the memory kernel was obtained and discussed in a qualitative link to the mode-coupling equation. Received: 9 January 1998 / Revised: 24 April 1998 / Accepted: 2 July 1998     

Variational theory for a single polyelectrolyte chain

Variational methods are applied to a single polyelectrolyte chain. The polymer is modeled as a Gaussian chain with screened electrostatic repulsion between all monomers. As a variational Hamiltonian, the most general Gaussian kernel, including the possibility of a classical or mean polymer path, is employed. The resulting self-consistent equations are systematically solved both for large and small monomer-monomer separations along the chain. In the absence of screening, the polymer is stretched on average. It is described by a straight classical path with Gaussian fluctuations around it. If the electrostatic repulsion is screened, the polymer is isotropically swollen for large separations, and for small separations the polymer correlation function is calculated as an analytic expansion in terms of the monomer-monomer separation along the chain. The electrostatic persistence length and the electrostatic blobsize are inferred from the crossover between distinct scaling ranges. We perform a global analysis of the scaling behavior as a function of the screening length and electrostatic interaction strength , where is the Bjerrum length and A is the distance of charges along the polymer chain. We find three different scaling regimes. i) A Gaussian-persistent regime with Gaussian behavior at small, persistent behavior at intermediate, and isotropically swollen behavior at large length scales. This regime occurs for weakly charged polymers and only for intermediate values of the screening length. The electrostatic persistence length is defined as the crossover length between the persistent and the asymptotically swollen behavior and is given by and thus disagrees with previous (restricted) variational treatments which predict a linear dependence on the screening length .ii) A Gaussian regime with Gaussian behavior at small and isotropically swollen behavior at large length scales. This regime occurs for weakly charged polymers and/or strong screening, and the electrostatic repulsion between monomers only leads to subfluent corrections to Gaussian scaling at small separations. The concept of a persistence length is without meaning in this regime. iii) A persistent regime , where the chain resembles a stretched rod on intermediate and small scales. Here the persistence length is given by the original Odijk prediction, , if the overstretching of the chain is avoided. We also investigate the effects of a finite polymer length and of an additional excluded-volume interaction, which modify the resultant scaling behavior. Applications to experiments and computer simulations are discussed. Received 24 December 1997     

The role of the amorphous fraction for the equilibrium shape of polymer single crystals

The equilibrium state of polymer single crystals is considered by explicitly taking into account the amorphous fraction formed by loops and tails of the chains using a statistical model introduced by Muthukumar (Philos. Trans. R. Soc. London, Ser. A 361, 539 (2003)). We show that under realistic conditions below the equilibrium melting temperature, tight loops and close re-entries are favored, and that the amorphous fraction can be mapped into an excess surface free energy. The model is extended to many-chain crystals where it is shown that the lamellar thickness increases with the number of chains in the crystal and extended-chain conformations are thermodynamically favored if the number of chains in the crystal is sufficiently large. The number of chains necessary to form an extended-chain crystal in thermodynamic equilibrium scales with the square of the degree of polymerization of the chains. We discuss the temperature behavior of the equilibrium crystal thickness in the under-cooled state.     

Monte carlo studies of phase transitions in polymer blends and block copolymer melts

Summary The unmixing transition of both symmetrical polymer blends AB (i.e. chain lengthsN _A=N _B=N) and asymmetrical ones (N _B/N _A=2,3) is studied by large-scale Monte Carlo simulations of the bond fluctuation model. Combination of semi-grand-canonical simulation techniques, ?histogram reweighting? and finitesize scaling allows an accurate location of the coexistence curve in the critical region. The variation of the critical temperature with chain length (N) is studied and compared to theoretical predictions. For the symmetrical case, use of chain lengths up toN=512 allows a rough estimation of crossover scaling functions for the crossover from Ising to mean-field exponents. The order-disorder transitions in melts of both symmetric (compositionf=N _A/(N _A+N _B)=1/2) and asymmetric (f=3/4) block copolymers is studied for very short chains (16≤N≤60). The interplay between structure and chain configuration is emphasized. Qualitative evidence for ?dumbell formation? of chains and vacancy enrichment in A-B-interfaces and near hard walls is presented. Paper presented at the I International Conference on Scaling Concepts and Complex Fluids, Copanello, Italy, July 4–8, 1994.