Integral equation theory for athermal solutions of linear polymers

An integral equation model is developed for athermal solutions of flexible linear polymers with particular reference to good solvent conditions. Results from scaling theory are used in formulating form factors for describing the single chain structure, and the impact of solvent quality on the chain fractal dimension is accounted for. Calculations are performed within the stringlike implementation of the polymer reference interaction site model with blobs (as opposed to complete chains) treated as the constituent structural units for semidilute solutions. Results are presented for the second virial coefficient between polymer coils and the osmotic compressibility as functions of the chain length and polymer volume fraction, respectively. Findings from this model agree with results from scaling theory and experimental measurements, as well as with an earlier investigation in which self-avoiding chains were described using Gaussian form factors with a chain length and concentration-dependent effective statistical segment length. The volume fractions at the threshold for connectedness percolation are evaluated within a coarse-grained closure relation for the connectedness Ornstein-Zernike equation. Results from these calculations are consistent with the usual interpretation of the semidilute crossover concentration for model solutions of both ideal and swollen polymer coils.     

A highly coarse-grained model to simulate entangled polymer melts

We introduce a highly coarse-grained model to simulate the entangled polymer melts. In this model, a polymer chain is taken as a single coarse-grained particle, and the creation and annihilation of entanglements are regarded as stochastic events in proper time intervals according to certain rules and possibilities. We build the relationship between the probability of appearance of an entanglement between any pair of neighboring chains at a given time interval and the rate of variation of entanglements which describes the concurrence of birth and death of entanglements. The probability of disappearance of entanglements is tuned to keep the total entanglement number around the target value. This useful model can reflect many characteristics of entanglements and macroscopic properties of polymer melts. As an illustration, we apply this model to simulate the polyethylene melt of C(1000)H(2002) at 450 K and further validate this model by comparing to experimental data and other simulation results.     

Monte Carlo simulations of polymer dynamics: Recent advances

A brief review is given of applications of Monte Carlo simulations to study the dynamical properties of coarse-grained models of polymer melts, emphasizing the crossover from the Rouse model toward reptation, and the glass transition. The extent to which Monte Carlo algorithms can mimic the actual chain dynamics is critically examined, and the need for the use of coarse-grained rather than fully atomistic models for such simulations is explained. It is shown that various lattice and continuum models yield qualitatively similar results, and the behavior agrees with the findings of corresponding molecular dynamics simulations and experiments, where available. It is argued that these simulations significantly enhance our understanding of the theoretical concepts on the dynamics of dense macromolecular systems. © 1997 John Wiley & Sons, Inc.     

Effect of molecular weight on the capillary absorption of polymer droplets

We study capillary absorption of small polymer droplets into nonwettable capillaries using coarse-grained molecular dynamics simulations and a simple analytical model. Studies of droplets of simple fluids have revealed that the capillary process depends on the ratio of tube-to-droplet radii [Willmott Faraday Discuss., 2010, 146, 233; Marmur J. Colloid Interface Sci. 1988, 122, 209]. Here we consider the absorption of droplets of polymers and study the effect of polymer chain length on the capillary absorption process. Our simulations reveal that for droplets of the same size (radius), the critical tube radius, below which there is no absorption, increases with the length of the polymer chains that constitute the droplets. We propose a model to explain this effect, which incorporates an entropic penalty for polymer confinement and find that this model agrees quantitatively with the simulations. We also find that the absorption dynamics is sensitive to the polymer chain length. In some cases during the capillary uptake transient partial absorption states, where the droplet is partially in and partially out of the tube, were observed. Such dynamics cannot be explained by a generalized Lucas-Washburn approach.     

Osmotic pressure and polymer size in semidilute polymer solutions under good-solvent conditions

We consider the lattice Domb-Joyce model at a value of the coupling for which scaling corrections approximately vanish and determine the universal scaling functions associated with the osmotic pressure and the polymer size for semidilute polymer solutions (c/c( *)相似文献   

Universal properties of a single polymer chain in slit: Scaling versus molecular dynamics simulations

We revisit the classical problem of a polymer confined in a slit in both of its static and dynamic aspects. We confirm a number of well known scaling predictions and analyze their range of validity by means of comprehensive molecular dynamics simulations using a coarse-grained bead-spring model of a flexible polymer chain. The normal and parallel components of the average end-to-end distance, mean radius of gyration and their distributions, the density profile, the force exerted on the slit walls, and the local bond orientation characteristics are obtained in slits of width D=4/10 (in units of the bead diameter) and for chain lengths N=50/300. We demonstrate that a wide range of static chain properties in normal direction can be described quantitatively by analytic model-independent expressions in perfect agreement with computer experiment. In particular, the observed profile of confinement-induced bond orientation is shown to closely match theory predictions. The anisotropy of confinement is found to be manifested most dramatically in the dynamic behavior of the polymer chain. We examine the relation between characteristic times for translational diffusion and lateral relaxation. It is demonstrated that the scaling predictions for lateral and normal relaxation times are in good agreement with our observations. A novel feature is the observed coupling of normal and lateral modes with two vastly different relaxation times. We show that the impact of grafting on lateral relaxation is equivalent to doubling the chain length.     

Identifying the primitive path mesh in entangled polymer liquids

Similar to entangled ropes, polymer chains cannot slide through each other. These topological constraints, the so-called entanglements, dominate the viscoelastic behavior of high-molecular-weight polymeric liquids. Tube models of polymer dynamics and rheology are based on the idea that entanglements confine a chain to small fluctuations around a primitive path which follows the coarse-grained chain contour. To establish the microscopic foundation for these highly successful phenomenological models, we have recently introduced a method for identifying the primitive path mesh that characterizes the microscopic topological state of computer-generated conformations of long-chain polymer melts and solutions. Here we give a more detailed account of the algorithm and discuss several key aspects of the analysis that are pertinent for its successful use in analyzing the topology of the polymer configurations. We also present a slight modification of the algorithm that preserves the previously neglected self-entanglements and allows us to distinguish between local self-knots and entanglements between distant sections of the same chain. Our results indicate that the latter make a negligible contribution to the tube and that the contour length between local self-knots, N_lk is significantly larger than the entanglement length N_e. © 2005 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 43: 917–933, 2005     

Chain overlap and entanglements in dilute polymer solutions: Brownian dynamics simulation

We have analyzed chain conformations and the existence — or otherwise — of chain overlaps and entanglements in dilute polymer solutions (at concentrations c < C^*, c^* = critical concentration). The fundamental problem of existence of chain overlaps in dilute solutions is also related to the drag reduction phenomenon (DR). Some experimental results pertinent to DR are explained in terms of entanglements even for solutions at concentrations defined in ppm. We report results of Brownian dynamics simulations of polymer solutions in which the equations of motion of the chains are solved by using the Langevin equation. Chains move according to actions of a systematic frictional force and a randomly fluctuating force w(t), where t is time. In addition, a shear flow field can be introduced into the model. To evaluate the structure of polymer chains in solution we have devised a measure of interchain contacts and two different measures of entanglements. The results for c = 0.3 c^* demonstrate that both chain entanglements and overlaps take place even in dilute solution. They also confirm predictions from an earlier combinatorial model.     

Long-range correlations, geometrical structure, and transport properties of macromolecular solutions. The equivalence of configurational statistics and geometrodynamics of large molecules

A special theory of Brownian relativity was previously proposed to describe the universal picture arising in ideal polymer solutions. In brief, it redefines a Gaussian macromolecule in a 4-dimensional diffusive spacetime, establishing a (weak) Lorentz-Poincaré invariance between liquid and polymer Einstein's laws for Brownian movement. Here, aimed at inquiring into the effect of correlations, we deepen the extension of the special theory to a general formulation. The previous statistical equivalence, for dynamic trajectories of liquid molecules and static configurations of macromolecules, and rather obvious in uncorrelated systems, is enlarged by a more general principle of equivalence, for configurational statistics and geometrodynamics. Accordingly, the three geodesic motion, continuity, and field equations could be rewritten, and a number of scaling behaviors were recovered in a spacetime endowed with general static isotropic metric (i.e., for equilibrium polymer solutions). We also dealt with universality in the volume fraction and, unexpectedly, found that a hyperscaling relation of the form, (average size) x (diffusivity) x (viscosity)1/2 ~f(N0, phi0) is fulfilled in several regimes, both in the chain monomer number (N) and polymer volume fraction (phi). Entangled macromolecular dynamics was treated as a geodesic light deflection, entaglements acting in close analogy to the field generated by a spherically symmetric mass source, where length fluctuations of the chain primitive path behave as azimuth fluctuations of its shape. Finally, the general transformation rule for translational and diffusive frames gives a coordinate gauge invariance, suggesting a widened Lorentz-Poincaré symmetry for Brownian statistics. We expect this approach to find effective applications to solutions of arbitrarily large molecules displaying a variety of structures, where the effect of geometry is more explicit and significant in itself (e.g., surfactants, lipids, proteins).     

Dissipative particle dynamics simulation of depletion layer and polymer migration in micro- and nanochannels for dilute polymer solutions

The flows of dilute polymer solutions in micro- and nanoscale channels are of both fundamental and practical importance in variety of applications in which the channel gap is of the same order as the size of the suspended particles or macromolecules. In such systems depletion layers are observed near solid-fluid interfaces, even in equilibrium, and the imposition of flow results in further cross-stream migration of the particles. In this work we employ dissipative particle dynamics to study depletion and migration in dilute polymer solutions in channels several times larger than the radius of gyration (Rg) of bead-spring chains. We compare depletion layers for different chain models and levels of chain representation, solvent quality, and relative wall-solvent-polymer interactions. By suitable scaling the simulated depletion layers compare well with the asymptotic lattice theory solution of depletion near a repulsive wall. In Poiseuille flow, polymer migration across the streamlines increases with the Peclet and the Reynolds number until the center-of-mass distribution develops two symmetric off-center peaks which identify the preferred chain positions across the channel. These appear to be governed by the balance of wall-chain repulsive interactions and an off-center driving force of the type known as the Segre-Silberberg effect.     

Equilibrium conformational dynamics of a polymer in a solvent

Molecular dynamics simulations were used to study the conformational dynamics of a bead-spring model polymer in an explicit solvent under good solvent conditions. The dynamics of the polymer chain were investigated using an analysis of the time autocorrelation functions of the Rouse coordinates of the polymer chain. We have investigated the variation of the correlation functions with polymer chain length N, solvent density rho, and system size. The measured initial decay rates gamma(p) of the correlation functions were compared with the predictions from a theory of polymer dynamics which uses the Oseen tensor to describe hydrodynamic interactions between monomers. Over the range of chain lengths considered (N = 30-60 monomers), the predicted scaling of gamma(p) proportional to N(-3nu) was observed at high rho, where nu is the polymer scaling exponent. The predicted gamma(p) are generally higher than the measured values. This discrepancy increases with decreasing rho, as a result in the breakdown in the conditions required for the Oseen approximation. The agreement between theory and simulation at high rho improves considerably if the theoretical expression for gamma(p) is modified to avoid sum-to-integral approximations, and if the values of (R(p)2), which are used in the theory, are taken directly from the simulation rather than being calculated using approximate scaling relations. The observed finite-size scaling of gamma(p) is not quantitatively consistent with the theoretical predictions.     

Microscopic theory of gelation and elasticity in polymer-particle suspensions

A simplified mode-coupling theory (MCT) of ergodic-nonergodic transitions, in conjunction with an accurate two-component polymer reference interaction site model (PRISM) theory for equilibrium structural correlations, has been systematically applied to investigate gelation, localization, and elasticity of flexible polymer-hard particle suspensions. The particle volume fraction at the fluid-gel transition is predicted to depend exponentially on reduced polymer concentration and size asymmetry ratio at relatively high colloid concentrations. In contrast, at lower particle volume fractions, a power-law dependence on polymer concentration is found with effective exponents and prefactors that depend systematically on the polymer/particle size ratio. Remarkable power-law and near universal scaling behavior is found for the localization length and elastic shear modulus. Multiple experiments for gel boundaries and shear moduli are in good agreement with the no adjustable parameter theory. The one exception is the absolute magnitude of the shear modulus which is strongly overpredicted, apparently due to nonequilibrium dense cluster formation. The simplified MCT-PRISM theory also captures the qualitative aspects of the weak depletion-driven "glass melting" phenomenon at high particle volume fractions. Calculations based on an effective one-component model of structure within a low particle volume fraction framework yield qualitatively different features than the two-component approach and are apparently all in disagreement with experiments. This suggests that volume fraction and size asymmetry dependent many-body screening of polymer-mediated depletion attractions at finite particle concentrations are important.     

Microscopic theory of orientational order, structure and thermodynamics in strained polymer liquids and networks

A microscopic integral equation theory of the segmental orientational order parameter, structural correlations and thermodynamics of strained polymer solutions, melts and networks has been developed. The nonclassical problem of the consequences of intermolecular excluded volume repulsions and chain connectivity is addressed. The theory makes several novel predictions, including effective power law dependences of the orientational order parameter on monomer concentration and chain degree of polymerization, and strain hardening of the bulk modulus. The predictions of a nearly classical strain dependence, and supralinear scaling with segment concentration, of the strain-induced nematic order parameter is in agreement with nuclear magnetic resonance experiments. The absolute magnitudes of the a priori calculated orientational order parameter agree with simulations and experiments to within a factor of 2. The possible complicating influence of "trapped entanglements" in crosslinked networks is discussed. Extensions of the theory are possible to treat the mechanical response of flexible polymer liquids and rubbers, and the structure, thermodynamics, and mechanical properties of strained liquid crystal forming polymers.     

Computer simulation of linkage of two ring chains

We performed off-lattice Monte Carlo simulations of links of two model ring chains with chain length N up to 32,768 in the theta solution or amorphous bulk state by using a random walk model (Model I), and molecular dynamics simulations of two model ring chains in solution with excluded volume interaction (Model II) to investigate topological effects on the geometry of link and ring conformation. In the case of Model I, the mean squared linking number, its distribution, and the size of two chains with fixed linking number are investigated. Our simulation results confirm the previous theoretical prediction that the mean squared linking number decays as pe(-qs(2)) with the distance of centers of chain mass s, where p and q are found to be chain length dependent and q asymptotically approaches to 0.75 as chain length increases. The linking number distribution of two chains has a universal form for long chains, but our simulation results clearly show that the distribution function deviates from the Gaussian distribution, a fact not predicted by any previous theoretical work. A scaling prediction is proposed to predict the link size, and is checked for our simulations for the Model II. The simulation results confirmed the scaling prediction of the blob picture that the link with linking number m occupies a compact volume of m blobs, and the size of the link is asymptotic to R(L) ≈ bN(ν)m(1/3-ν), where N is the chain length, and v is the Flory exponent of polymer in solutions.     

Molecular dynamics simulations of supramolecular polymer rheology

Using equilibrium and nonequilibrium molecular dynamics simulations, we studied the equilibrium and rheological properties of dilute and semidilute solutions of head-to-tail associating polymers. In our simulation model, a spontaneous complementary reversible association between the donor and the acceptor groups at the ends of oligomers was achieved by introducing a combination of truncated pseudo-Coulombic attractive potential and Lennard Jones repulsive potential between donor, acceptor, and neighboring groups. We have calculated the equilibrium properties of supramolecular polymers, such as the ring/chain equilibrium, average molecular weight, and molecular weight distribution of self-assembled chains and rings, which all agree well with previous analytical and computer modeling results. We have investigated shear thinning of solutions of 8- and 20-bead associating oligomers with different association energies at different temperatures and oligomer volume fractions. All reduced viscosity data for a given oligomer length can be collapsed into one master curve, exhibiting two power-law regions of shear-thinning behavior with an exponent of -0.55 at intermediate ranges of the reduced shear rate β and -0.8 (or -0.9) at larger shear rates. The equilibrium viscosity of supramolecular solutions with different oligomer lengths and associating energies is found to obey a power-law scaling dependence on oligomer volume fraction with an exponent of 1.5, in agreement with the experimental observations for several dilute or semidilute solutions of supramolecular polymers. This implies that dilute and semidilute supramolecular polymer solutions exhibit high polydispersity but may not be sufficiently entangled to follow the reptation mechanism of relaxation.     

Anomalous nanoparticle diffusion in polymer solutions and melts: a mode-coupling theory study

Mode-coupling theory is employed to study diffusion of nanoparticles in polymer melts and solutions. Theoretical results are directly compared with molecular dynamics simulation data for a similar model. The theory correctly reproduces the effects of the nanoparticle size, mass, particle-polymer interaction strength, and polymer chain length on the nanoparticle diffusion coefficient. In accord with earlier experimental, simulation, and theoretical work, it is found that when the polymer radius of gyration exceeds the nanoparticle radius, the Stokes-Einstein relation underestimates the particle diffusion coefficient by as much as an order of magnitude. Within the mode-coupling theory framework, a microscopic interpretation of this phenomenon is given, whereby the total diffusion coefficient is decomposed into microscopic and hydrodynamic contributions, with the former dominant in the small particle limit, and the latter dominant in the large particle limit. This interpretation is in agreement with previous mode-coupling theory studies of anomalous diffusion of solutes in simple dense fluids.     

高分子相分离的分子动力学模拟

用粗粒化分子动力学(MD)模拟方法从分子层次研究两组分聚合物共混体系相分离过程中的动力学. 在相分离初期, 相区尺寸不随时间增加而变化; 在相分离中期, 相区尺寸与时间有很好的标度关系, 标度指数(α=1/3)符合Lifshiz-Slyozov提出的以扩散为主导的蒸发-凝聚机理的标度预测; 在相分离后期, 体系实现宏观相分离, 相区尺寸不再随时间改变而变化. 体积分数小的高分子链尺寸在相分离过程中先收缩再扩张, 在实现宏观相分离后, 高分子链尺寸又回到本体状态尺寸.     

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.