Constrained Rouse model of rubber viscoelasticity

In this work we use a new approach to investigate the equilibrium and linear dynamic-mechanical response of a polymer network. The classical Rouse model is extended to incorporate quenched constraints on its end-boundary conditions; a microscopic stress tensor for the network system is then derived in the affine deformation limit. To test the model we calculate the macroscopic stress in equilibrium, corresponding to the long-time limit of relaxation. Particular attention is paid to the treatment of compressibility and hydrostatic pressure in a sample with open boundaries. Although quite different in general, for small strains the model compares well with the classic equilibrium rubber-elasticity models. The dynamic shear modulus is obtained for a network relaxing after an instantaneous step strain by keeping track of relaxation of consecutive Rouse modes of constrained network strands. The results naturally cover the whole time range--from the dynamic glassy state down to the equilibrium incompressible rubber plateau.     

Theory of relaxation properties of two‐dimensional polymer networks, 2. Local dynamic characteristics

Using normal mode transformation obtained in Part 1 of this series^[1], the exact analytical expressions for the mean‐square displacements of junctions and non‐junction beads, the autocorrelation functions of the end‐to‐end chain vectors between neighboring junctions, and those of subchain vectors of a two‐dimensional regular network consisting of "bead and spring" Rouse chains are obtained. Contributions of intra‐ and interchain relaxation processes to the local dynamic characteristics considered are compared. The time behavior of dynamic quantities obtained is estimated for different scales of motions. The possibility of describing long‐time relaxation of a two‐dimensional network by a simplified coarse‐grained network model is demonstrated. It is shown that the local relaxation properties of a two‐dimensional polymer network (as well as a three‐dimensional network) on scales smaller than the average distance between cross‐links are very close to those of a single Rouse chain. The large‐scale collective relaxation of the polymer networks having a two‐dimensional connectivity differs considerably from that of the three‐dimensional networks.     

Theory of relaxation properties of two‐dimensional polymer networks, 1. Normal modes and relaxation times

The local relaxation properties of polymer networks with a two‐dimensional connectivity are considered. We use the mesh‐like network model in which the average positions of junctions form the regular spatial structure consisting of square repeating units (network cells). The two‐dimensional polymer network consisting of “bead and spring” Rouse chains and the simplified coarse‐grained network model describing only the large‐scale collective relaxation of a network are studied. For both dynamic network models the set of relaxation times and the transformation from Cartesian coordinates of network elements to normal modes are obtained. Using the normal mode transformation obtained, in Part 2 of this series the exact analytical expressions for various local dynamic characteristics of the polymer network having a two‐dimensional connectivity will be calculated.     

Linear viscoelastic properties of transient networks formed by associating polymers with multiple stickers

We have developed a single-chain theory that describes dynamics of associating polymer chains carrying multiple associative groups (or stickers) in the transient network formed by themselves and studied linear viscoelastic properties of this network. It is shown that if the average number N of stickers associated with the network junction per chain is large, the terminal relaxation time τ(A) that is proportional to τ(X)N(2) appears. The time τ(X) is the interval during which an associated sticker goes back to its equilibrium position by one or more dissociation steps. In this lower frequency regime ω<1/τ(X), the moduli are well described in terms of the Rouse model with the longest relaxation time τ(A). The large value of N is realized for chains carrying many stickers whose rate of association with the network junction is much larger than the dissociation rate. This associative Rouse behavior stems from the association/dissociation processes of stickers and is different from the ordinary Rouse behavior in the higher frequency regime, which is originated from the thermal segmental motion between stickers. If N is not large, the dynamic shear moduli are well described in terms of the Maxwell model characterized by a single relaxation time τ(X) in the moderate and lower frequency regimes. Thus, the transition occurs in the viscoelastic relaxation behavior from the Maxwell-type to the Rouse-type in ω<1/τ(X) as N increases. All these results are obtained under the affine deformation assumption for junction points. We also studied the effect of the junction fluctuations from the affine motion on the plateau modulus by introducing the virtual spring for bound stickers. It is shown that the plateau modulus is not affected by the junction fluctuations.     

Monte Carlo simulations of stress relaxation of entanglement-free Fraenkel chains. I. Linear polymer viscoelasticity

Shear stress relaxation modulus GS(t) curves of entanglement-free Fraenkel chains have been calculated using Monte Carlo simulations based on the Langevin equation, carrying out both in the equilibrium state and following the application of a step shear deformation. While the fluctuation-dissipation theorem is perfectly demonstrated in the Rouse-chain model, a quasiversion of the fluctuation-dissipation theorem is observed in the Fraenkel-chain model. In both types of simulations on the Fraenkel-chain model, two distinct modes of dynamics emerge in GS(t), giving a line shape similar to that typically observed experimentally. Analyses show that the fast mode arises from the segment-tension fluctuations or reflects the relaxation of the segment tension created by segments being stretched by the applied step strain-an energetic-interactions-driven process-while the slow mode arises from the fluctuations in segmental orientation or represents the randomization of the segmental-orientation anisotropy induced by the step deformation-an entropy-driven process. Furthermore, it is demonstrated that the slow mode is well described by the Rouse theory in all aspects: the magnitude of modulus, the line shape of the relaxation curve, and the number-of-beads (N) dependence of the relaxation times. In other words, one Fraenkel segment substituting for one Rouse segment, it has been shown that the entropic-force constant on each segment is not a required element to give rise to the Rouse modes of motion, which describe the relaxation modulus of an entanglement-free polymer over the long-time region very well. This conclusion provides an explanation resolving a long-standing fundamental paradox in the success of Rouse-segment-based molecular theories for polymer viscoelasticity-namely, the paradox between the Rouse segment size being of the same order of magnitude as that of the Kuhn segment (each Fraenkel segment with a large force constant HF can be regarded as basically equivalent to a Kuhn segment) and the meaning of the Rouse segment as defined in the Rouse-chain model. The general agreement observed in the comparison of the simulation and experimental results indicates that the Fraenkel-chain model, while being still relatively simple, has captured the key element in energetic interactions--the rigidity on the segment--in a polymer system.     

Analysis of Polymer Dynamics by NMR Modulated Gradient Spin Echo

NMR modulated gradient spin echo method, which allows the quantification of polymer segmental displacement via the measurement of the velocity autocorrelation, requires the formulation of theoretical predictions in the frequency domain in order to test their validity. We worked out the velocity autocorrelation spectrum of the segmental motion according to the Rouse and the tube/reptation model to compare it to the experimental results obtained by the new NMR technique. The analysis of polybutadiene shows the segmental velocity autocorrelation spectrum typical for the reptation-like motion of polymer in a “tube”. The measurement of bulk water indicates a kind of Rouse motion in a network of hydrogen bonds.     

二维三角格子上高分子动态 Monte Carlo模拟的协同运动算法

在动态Monte Carlo模拟的协同运动算法中,几个相邻的链节可以同时运动,这可以理解为高分子链中张力的作用引起的协同运动。将这一算法用于二维三角格子模型上RW链和SAW链的模拟。结果表明RW链的动力学行为符合Rouse理论,说明说明该算法可以用于高分子动力学研究,其优点是不需要使用键长涨落模型。     

Long time dynamics of a polymer with rigid body monomer units relating to a protein model: Comparison with the rouse model

Simulation data from an off-lattice polymer model are compared with data from the Rouse model. The polymer model is built of sequentially connected rigid monomer units that represent the amide planes of a protein backbone. The time propagation of the dynamics of the polymer model is performed by a Monte Carlo method. The elementary Monte Carlo steps correspond to local confomational changes in a window of three consecutive monomer units. The time autocorrelation functions of end-to-end distances from segments within the linear chain molecule are considered in detail. Both models exhibit a stretched exponential decay pattern. A comparison of the data from the Rouse model and the computer simulation provide an estimate of the time unit of 15 ps for a full scan of the algorithm for local conformational changes along the chain. With a conservative estimate of the parameters governing the Rouse model this time unit is four orders of magnitude longer than the elementary time step of a conventional computer simulation of polymer dynamics based upon the classical equations of motion for all atoms. © 1992 by John Wiley & Sons, Inc.     

高分子动态Monte Carlo模拟的协同运动算法

在高分子动力学的研究中 ,动态 Monte Carlo模拟发挥了重要的作用 [1] .动态 Monte Carlo模拟的关键是选择具有物理真实性的高分子运动算法 .目前广为采用的算法是经 Hilhorst和 Deutch[2 ] 修正的 Verdier- Stockmayer算法 [3] 以及 Carmesin和 Kremer等的键长涨落算法 [4 ] ,陆建明和杨玉良曾提出一种高分子动态算法 [5] ,在他们的算法中很长一段链节可能作蛇行运动 ,但是冯捷等 [6 ] 指出这种运动模式不满足微观可逆性条件 .本文对该运动模式进行修正 ,得到一种协同运动算法 ,并对其动力学行为进行检验 .1　算　法在平面正方形格子…     

Two tube models of rubber elasticity

Polymer entanglements lead to complicated topological constraints and interactions between neighboring chains in a dense solution or melt. Entanglements can be treated in a mean field approach, within the famous reptation model, since they effectively confine each individual chain in a tube-like geometry. In polymer networks, due to crosslinks preventing the global reptation and constraint release, entanglements acquire a different topological meaning and have a much stronger effect on the resulting mechanical response. In this article we discuss two different models of rubber elasticity, both utilizing the reptation ideas. First, we apply the classical ideas of reptation statistics to calculate the effective rubber-elastic free energy of an entangled rubbery network. In the second approach, we examine the classical Rouse dynamics of chains with quenched constraints at their ends by crosslinks, and along the primitive path by entanglements. We then proceed to average a microscopic stress tensor for the network system and present it in a manageable form in the equilibrium t → ∞ limit. Particular attention is paid to the treatment of compressibility and hydrostatic pressure in a sample with open boundaries. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 2679–2697, 2006     

高分子动力学的单链模型

高分子单链模型是高分子稀溶液理论研究的基本模型.对其进行深入地分析,不仅有助于解决高分子稀溶液体系中溶液黏度和分子链扩散等基本问题,而且能够增进人们对高分子链结构与溶液性质间关联性的理解.虽然基于经典连续性介质力学的流体动力学理论可以定性,甚至半定量地获得稀溶液的一些重要性质,但是,随着科学技术的发展,人们从分子水平上建立了许多描述高分子稀溶液性质的模型和理论,期望能够定量地描述高分子稀溶液的性质.本文以高分子稀溶液中3个典型的单链模型为例(包括:不含流体力学相互作用的Rouse模型、含二体流体力学相互作用的Zimm模型和含多体流体力学相互作用的部分穿透球模型),综述高分子稀溶液的重要性质,并详细地给出其动力学方程的推导过程及其重要的研究进展.特别是,对于Rouse模型,本文还将其预言结果拓展到了短链高分子流体体系;此外,还介绍了这一领域的关键科学问题、发展前景和研究方向.     

Coherent elastic scattering of a reptating chain in a melt during transient relaxation

We calculate the elastic scattering of a long chain in a polymer melt during the process of relaxation after a sudden deformation, i.e., the return to equilibrium with time, while the deformed shape is maintained. The scattering thus depends on the duration t of the relaxation, as well as on the scattering vector, and is a physical quantity characteristic of the dynamics of long chains. The reptation model of de Gennes is used as developed for deformed melts by Doi and Edwards: the chain is confined by other chains in a tube, from which disengagement by the ends is the only way of renewing its configuration. The tube diameter is taken as a parameter. We give both an analytic form and numerical evaluations. On comparison with calculations based on the Rouse model, in which the chain is assumed to be free in a viscous medium, it is seen that the experimental data (given elsewhere) should enable one to distinguish between the two models.     

Rheological study of transient networks with junctions of limited multiplicity

We theoretically study the viscoelastic and thermodynamic properties of transient gels comprised of telechelic associating polymers. We extend classical theories of transient networks so that correlations among polymer chains through the network junctions are taken into account. This extension enables us to investigate how rheological quantities such as elastic modulus, viscosity, and relaxation time are affected by the association equilibrium, and how these quantities are related to the aggregation number (or multiplicity) of the junctions. In this paper, we assume, in the conventional manner, that chains are elastically effective if both their ends are connected with other chains. It is shown that the dynamic shear moduli are well described in terms of the Maxwell model. As a result of the correlation, the reduced moduli (moduli divided by the polymer concentration) increase with the concentration, but become independent of the concentration in the high-concentration range. The fraction of pairwise junctions is larger at lower concentrations, indicating the presence of concatenated chains in the system, which decreases as the concentration increases. This leads to a network relaxation time that increases with the concentration.     

Combined Molecular Dynamics Simulation and Rouse Model Analysis of Static and Dynamic Properties of Unentangled Polymer Melts with Different Chain Architectures

Chain architecture effect on static and dynamic properties of unentangled polymers is explored by molecular dynamics simulation and Rouse mode analysis based on graph theory. For open chains, although they generally obey ideal scaling in chain dimensions, local structure exhibits nonideal behavior due to the incomplete excluded volume(EV) screening, the reduced mean square internal distance(MSID) can be well described by Wittmer' theory for linear chains and the resulting chain swelling is architecture dependent, i.e., the more branches a bit stronger swelling. For rings, unlike open chains they are compact in term of global sizes. Due to EV effect and nonconcatenated constraints their local structure exhibits a quite different non-Gaussian behavior from open chains, i.e., reduced MSID curves do not collapse to a single master curve and fail to converge to a chain-length-independent constant, which makes the direct application of Wittmer's theory to rings quite questionable.Deviation from ideality is further evidenced by limited applicability of Rouse prediction to mode amplitude and relaxation time at high modes as well as the non-constant and mode-dependent scaled Rouse mode amplitudes, while the latter is architecture-dependent and even molecular weight dependent for rings. The chain relaxation time is architecture-dependent, but the same scaling dependence on chain dimensions does hold for all studied architectures. Despite mode orthogonality at static state, the role of cross-correlation in orientation relaxation increases with time and the time-dependent coupling parameter rises faster for rings than open chains even at short time scales it is lower for rings.     

Why does the rouse model fairly describe the dynamic characteristics of polymer melts at molecular masses below critical mass?

Generalization of the Rouse model without any use of the postulates concerning the Gaussian distribution of the vector connecting the ends of segments is advanced. In the initial (in general, nonlinear) Langevin equations, self-averaging over continuous fragments of a macromolecule naturally defines a linear term for the tagged chain, and this term differs from the entropy term of the classical Rouse model only by the numerical coefficient. According to the inertia-free approximation, the initial decay rates of correlation functions for the normal modes are described by the Rouse model independently of the character of fluctuations of the vector connecting the ends of the Kuhn segment. This statement is valid for any moment if the initial Langevin equations are treated in terms of the approximation of dynamic self-consistency. Simulation of the Fraenkel chains by the method of Brownian dynamics shows that decay of autocorrelation functions of shortwave normal modes is fairly described by the linearized equations for a given model of a chain and that the Rouse equation can be used for the long-wave modes. The results of this study make it possible to explain a marked difference between the lengths of the Kuhn and Rouse segments that is estimated from static and dynamic experiments.     

From caging to Rouse dynamics in polymer melts with intramolecular barriers: a critical test of the mode coupling theory

By means of computer simulations and solution of the equations of the mode coupling theory (MCT), we investigate the role of the intramolecular barriers on several dynamic aspects of nonentangled polymers. The investigated dynamic range extends from the caging regime characteristic of glass-formers to the relaxation of the chain Rouse modes. We review our recent work on this question, provide new results, and critically discuss the limitations of the theory. Solutions of the MCT for the structural relaxation reproduce qualitative trends of simulations for weak and moderate barriers. However, a progressive discrepancy is revealed as the limit of stiff chains is approached. This disagreement does not seem related with dynamic heterogeneities, which indeed are not enhanced by increasing barrier strength. It is not connected either with the breakdown of the convolution approximation for three-point static correlations, which retains its validity for stiff chains. These findings suggest the need of an improvement of the MCT equations for polymer melts. Concerning the relaxation of the chain degrees of freedom, MCT provides a microscopic basis for time scales from chain reorientation down to the caging regime. It rationalizes, from first principles, the observed deviations from the Rouse model on increasing the barrier strength. These include anomalous scaling of relaxation times, long-time plateaux, and nonmonotonous wavelength dependence of the mode correlators.     

Stiffness and excluded volume effects on conformation and dynamics of polymers: A simulation study

This work investigates the effects of the excluded volume and especially those of the chain stiffness on the structural and dynamical properties of a model polymer chain.The theoretical framework is the same as in the recent works by Steinhauser et al.,where a Rouse approach is adopted.Our model differs in that our chains have a finite average bending angle.As in the works by Steinhauser et al.,Langevin dynamic simulations were performed without hydrodynamic interactions.Whereas this doesn’t impact the static properties we obtain,it also allows us to compare our results on dynamic properties to those predicted by Rouse theory,where hydrodynamic interactions are also neglected.Our results show that the structural properties are very sensitive to the chain stiffness,whereas the dynamic scaling laws remain the same as those by Rouse theory,with the prefactor depending on the persistence length.     

Relaxation of polymer thin films in isothermal temperature‐jump measurements

The dynamic behavior of thin polymer films is of interest in the fabrication of microelectronics and optoelectronics and in the coatings industry. It is known that polymer relaxation is affected by film thickness and the particular substrate/polymer pair. We previously used a spectroscopic ellipsometer to investigate the glass transition in thin films. In addition to information on the modification of thermal transitions such as the glass‐transition temperature, the speed of data acquisition in an automated, spectroscopic ellipsometer, operated at a single wavelength of 780 nm, allows for the direct observation of the isothermal dimensions of a thin polymer film as a function of time after a rapid temperature change. In this article, we discuss recent results from the observation of the time dependence of film‐normal thickness and normalized, in‐plane, lateral dimension as well as simple fits to this relaxation behavior in terms of a normalized viscosity and relaxation time. The results support a highly asymmetric initial thermal expansion normal to the film followed by close to isotropic relaxation and anisotropic “flow” (the flow in response to the vanishingly small shears of thermal expansion). These features may clarify issues involving the observation of chain confinement in thin polymer films in terms of potential differences between equilibrium and dynamic measurements. © 2000 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 38: 2929–2936, 2000     

Dynamics and conformational properties of Rouse polymers in random layered flows

Here we study the dynamics and the conformational properties of Rouse polymers in random layered flows. We obtain explicit expressions for the mean squared displacement of the center of mass of the polymer as a function of time and number of beads. The averages involve both the thermal noise and the realizations of the random flows. We find also an explicit dependence of the equilibrium value of the end-to-end vector on the number of beads and other dimensional parameters. Our results are compared to that for Rouse polymers immersed in non-random steady shear flows.