Entanglement effects in model polymer networks

Many fundamental questions for the understanding of polymer networks are more suitably addressed by current computer simulations than by experiments. Details of the microscopic topology, such as the elastically active cluster or loop entanglements, can be identified as well as controlled. In particular, it is possible to isolate and quantify their effects on macroscopic observables such as the elastic modulus. The constraints due to connectivity and conserved topology are more clearly present for networks than for melts. Already for strand lengths between crosslinks which are relatively short, the effect of the conserved topology is important. The mode relaxation in a network is significantly different from that of a melt. For weakly crosslinked systems the melt entanglement length is the relevant scaling parameter. The elastic modulus of a long chain network under ideal conditions reaches an asymptotic value which is about 2.2 times smaller than the prediction of an affine model for a network made of strands of the melt entanglement length. An analysis of the stress reveals that in the linear regime the contribution from the excluded volume is dominant compared to that from the connectivity along the strands. For larger elongations, however, the non-linear elastic response is dominated by the chemically and topologically shortest paths through the system, where chemical crosslinks and topological entanglements between meshes of the network play a similarly important role.     

Dynamic mechanical properties of cross-linked rubbers. V. An approximate analysis of the modulus contributions of trapped and untrapped entanglements

Using the criterion that a coupling entanglement in a polymer network is trapped if all four strands radiating from it terminate in chemical cross-links (an approximation to a more rigorous treatment of Langley), an equation is derived relating the equilibrium modulus to the magnitude of the compliance in the frequency region where all coupling entanglements, whether trapped or not, contribute to the elasticity. The latter value is estimated from the storage compliance at the frequency where the storage compliance of the uncross-linked polymer corresponds to the entanglement compliance derived from integration over the loss compliance. The theory agrees rather well with data on vulcanizates of natural rubber. For 1,4-polybutadiene and styrene–butadiene rubbers, the agreement is somewhat less satisfactory, but the results support the hypothesis that the low-frequency losses observed in lightly cross-linked rubbers are due to relaxation of untrapped entanglements.     

Structure rheological investigations of segment interaction in polymer solutions of the network type

The solution structure of flexible chain polymers is described in terms of an equivalent entanglement network characterized by its apparent chain element A'. Such networks can exhibit free penetration of polymer segments, hindered penetration or no penetration at all. The segment interaction effects a hindrance to penetration. It can be characterized by the concentration dependence of the mass of network strand M_e resp. the entanglement density ν. Another parameter is the entanglement fraction ß, which describes the ratio of the entanglements actually formed to that possible in case of free penetration. Experimental results with different polymers illustrate these ideas.     

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.     

Effect of chain entanglements on plastic deformation behavior of ultra‐high molecular weight polyethylene

Samples of ultra‐high molecular weight polyethylene, in which the chain topology within the amorphous component was altered using two‐stage processing, including crystallization at high pressure in the first step, were produced and their deformation behavior in the plane‐strain compression was studied. Deformation and recovery experiments demonstrated that the state of the molecular network governed by entanglement density is one of the primary parameters controlling the response of the material on the imposed strain, especially at moderate and high strains. Any change in the concentration of entanglements markedly influences the shape of the true stress–true strain curve. The strain hardening modulus decreases while the onset of strain hardening increases with a decrease of the entanglement density within the amorphous component. Density of entanglements also influences the amount of rubber‐like recoverable deformation and permanent plastic flow. In material of the reduced concentration of entanglements permanent flow appears easier and sets in earlier than in the material with a higher entanglement density, becoming a favorable deformation mechanism at moderate strains. As a result, strong strain hardening is postponed to higher strain when compared with the samples of equilibrium entanglement density. In the samples of an increased entanglement density the molecular network becomes stiffer, with a reduced ability of strain induced disentangling of chains. Consequently, there is a less permanent flow and strain hardening begins earlier than in the reference material of an unaltered chain topology. © 2009 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 48: 276–285, 2010     

Relation of elastic modulus to crosslink and entanglement concentrations in rubber networks

The theory of rubber elasticity relates the elastic modulus of unfilled amorphous rubber to the concentration of elastically effective strands. A theoretical relation between this concentration and the concentrations of potential entanglements, random tetrafunctional crosslinks, and chain ends was proposed recently. In the present work, the new relation was combined with the theory of rubber elasticity and verified experimentally. Polydimethylsiloxane samples were cured by ⁶⁰Co irradiation and were extensively extracted to determine gel fraction, which was used to calculate concentrations of crosslinking and scission due to irradiation. Equilibrium modulus values determined from creep tests were in excellent agreement with those calculated using the new relation if the average spacing between potential entanglements is 116 (CH₃)₂SiO units. Thus, in typical commercial silicone rubbers, the contribution to the modulus from trapped entanglements is greater than the direct contribution from crosslinks. The new relation allows the calculation of crosslink concentrations from modulus measurements on other unfilled rubbers once the potential entanglement spacing of the polymer is determined.     

Pre-Averaged Sampling On the Entanglement Kinetics for Polymer Dynamics

A new model for entangled polymer dynamics based on pre-averaged sampling of the entanglement structure is proposed. Although it has been reported that sliplink simulations are powerful and promising to predict entangled polymer dynamics, it is still unpractical to calculate polymers with many entanglements. In the present study, a possible approach to achieve fast calculation is proposed by pre-averaged sampling of entanglement structure with skipping detail kinetics of entanglements dominated by chain ends in conventional sliplink models. To achieve time development of the chain conformation and entanglement structure, i) number of entanglement per chain and number of monomers for each segment are randomly obtained from the equilibrium distribution proposed by Schieber [J. Chem. Phys. 2003 , 118, 5162] and ii) the renewed entanglement structure is mechanically equilibrated. The established power-laws on molecular weight dependence of chain dimension, the longest relaxation time and self-diffusion coefficient were reasonably reproduced. Comparison on linear viscoelastic response is also discussed.     

Experimental study of molecular entanglement in polymer networks

Elastomeric molecular networks have been prepared by endlinking polydimethyl siloxane molecules having functional chain ends, both in the presence of an unreactive polymeric diluent and in the undiluted state. Values of tensile (Young) modulus were found to be in good agreement with the simple molecular theory of rubberlike elasticity for networks prepared in a highly diluted state. For concentrated systems the modulus was anomalously high, however. The discrepancy can be attributed to chain entanglements. A second interpenetrating network was introduced into networks formed in the diluted state by replacing the diluent polymer by reactive polymer, which was then gelled in situ. The modulus of these combined networks was much higher than the sum of the moduli of the constituent networks, implying a large contribution from molecular entanglements. © 1994 John Wiley & Sons, Inc.     

Models for adhesion at weak polymer interfaces

The weak interfaces between immiscible polymer pairs typically fail through chain scission. The critical facture toughness for such interfaces is closely related to the density of intermolecular entanglements at the interface. From scaling analysis, a simple correlation between facture toughness and chain entanglement was developed. It predicts well the interfacial adhesion for many immiscible polymer pairs found in the literature. For an interface with block copolymer reinforcement, its critical fracture toughness comes from both intermolecular entanglements of homopolymers and copolymer bridges. In the chain scission regime (low copolymer coverage), the block copolymer contribution is found proportional to copolymer interfacial coverage, with the coefficient being the energy to stretch and break a copolymer chain. The chain‐breaking energy for different copolymers was evaluated and compared to literature data. © 2009 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 47: 2313–2319, 2009     

Effect of chemical structure of polycarbonates on entanglement spacing

The master curves of a series of aliphatic polycarbonates(APCs) with different lengths of methylene segments in the repeat unit were obtained by dynamic rheological measurements.The plateau modulus and entanglement molecular weight were determined and cross-checked by different methods.Though having distinct difference in chemical structure of repeat units,both APCs and bisphenol-A polycarbonates have the similar entanglement weight and entanglement spacing.On the other side,the plateau modulus decreases with increasing the length of the side group of aliphatic polycarbonates with different side-chain lengths in the literature.The packing length model can explain the relationship between chain structure and entanglements.     

Strain hardening of polymer glasses: Effect of entanglement density,temperature, and rate

The strain hardening behavior of model polymer glasses is studied with simulations over a wide range of entanglement densities, temperatures, strain rates, and chain lengths. Entangled polymers deform affinely at scales larger than the entanglement length as assumed in entropic network models of strain hardening. The dependence of strain hardening on strain and entanglement density is also consistent with these models, but the temperature dependence has the opposite trend. The dependence on temperature, rate, and interaction strength can instead be understood as reflecting changes in the flow stress. Microscopic analysis of local rearrangements and the primitive paths between entanglements is used to test models of strain hardening. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 3487–3500, 2006     

Polymer blend/organoclay nanocomposite with poly(ethylene oxide) and poly(methyl methacrylate)

A series of polymer blend/organoclay nanocomposite at a fixed blending ratio was prepared using equal ratio of poly(ethylene oxide) (PEO) and poly(methyl methacrylate) (PMMA) via solvent casting method. With respect to nanoscale internal structure, we found that PMMA chains have better affinity with organoclay than PEO, based on the results from X-ray diffraction. Direct visualization via transmission electron microscopy (TEM) also supported the better affinity of PMMA with organoclay by indicating the existence of hybrid structures of mainly intercalated but with some exfoliated forms. The miscible nature of the blend and homogeneous dispersion state of layered silicate in the blend system were investigated via the microscopic fractured surface morphologies. From rheological measurements (storage and loss modulus), we discovered the role of the physical network structure between polymer and organoclay to be a main factor for the enhancement of elastic properties.     

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     

Viscoelastic properties of entangled polymers. II. The interchain–intrachain entanglement model

Our earlier model of entangled chain dynamics represented the elastic effects at entanglement points as a coupling with the surrounding medium, which was propagated weakly to every other entanglement site on the parent chain. In this way, a great variety of linear viscoelastic phenomena could be successfully predicted, with the only significant deficiency being a quantitative failure in curve fitting the local minimum in dynamic loss modulus G″(ω). Here we introduce the “intrachain entanglement,” or the “internal entanglement,” in which the parent chain entangles directly with itself. Elastic forces between pairs of internal entanglements are assumed to be stronger than those with the medium (“external entanglement”). Predictions of the new model are compared with rheological data on monodisperse polystyrene. Good agreement is obtained between theory and experiment, including G″(ω).     

Shear and extensional rheology of entangled polymer melts: Similarities and differences

This work extends our previous understanding concerning the nonlinear responses of entangled polymer solutions and melts to large external deformation in both simple shear and uniaxial extension. Many similarities have recently been identified for both step strain and startup continuous deformation, including elastic yielding, i.e., chain disentanglement after cessation of shear or extension, and emergence of a yield point during startup deformation that involves a deformation rate in excess of the dominant molecular relaxation rate. At a sufficiently high constant Hencky rate, uniaxial extension of an entangled melt is known to produce window-glass-like rupture. The present study provides evidence against the speculation that chain entanglements tie up into "dead knots" in constant-rate extension because of the exponentially growing chain stretching with time. In particular, it is shown that even Instron-style tensile stretching, i.e., extending a specimen by applying a constant velocity on both ends, results in rupture. Yet, in the same rate range, the same entangled melt only yields in simple shear, and the resulting shear banding is clearly not a characteristic of rupture. Thus, we conclude that chain entanglements respond to simple shear in the manner of yielding whereas uniaxial extension is rather effective in causing some entanglements to lock up, making it impossible for the entanglement network to yield at high rates.     

Effect of the functionality of junctions on the elasticity of polymacromonomer networks: Computer simulation

The elastic properties of polymer networks formed via the radical polymerization of macromonomers with two polymerizable end groups are studied via computer simulation. It is shown that variation in the average functionality of network junctions, f _avg, in a wide range (∼5–55) leads to a significant change in the shear modulus of the network. According to experiments with real networks (gels of poly(ethylene oxide) macromonomers), the shear modulus increases as f _avg increases. This effect is not due only to a decrease in the fluctuations of positions of network junctions. The main cause of the increase in the modulus is that the modulus component due to interaction between polymer chains (entanglements) increases as the functionality of junctions in the investigated networks increases. The conclusion is made that these networks gain entanglements during the formation of network junctions with high functionality rather than inherit them from the solution of macromonomer chains.     

Local and chain dynamics in miscible polymer blends: A Monte Carlo simulation study

Local chain structure and local environment play an important role in the dynamics of polymer chains in miscible blends. In general, the friction coefficients that describe the segmental dynamics of the two components in a blend differ from each other and from those of the pure melts. In this work, we investigate polymer blend dynamics with Monte Carlo simulations of a generalized bond fluctuation model, where differences in the interaction energies between nonbonded nearest neighbors distinguish the two components of a blend. Simulations employing only local moves and respecting a no bond crossing condition were carried out for blends with a range of compositions, densities, and chain lengths. The blends investigated here have long time dynamics in the crossover region between Rouse and entangled behavior. In order to investigate the scaling of the self-diffusion coefficients, characteristic chain lengths N(c) are calculated from the packing length of the chains. These are combined with a local mobility mu determined from the acceptance rate and the effective bond length to yield characteristic self-diffusion coefficients D(c)=muN(c). We find that the data for both melts and blends collapse onto a common line in a graph of reduced diffusion coefficients DD(c) as a function of reduced chain length NN(c). The composition dependence of dynamic properties is investigated in detail for melts and blends with chains of length N=20 at three different densities. For these blends, we calculate friction coefficients from the local mobilities and consider their composition and pressure dependence. The friction coefficients determined in this way show many of the characteristics observed in experiments on miscible blends.     

Chain orientation and intrinsic anisotropy in birefringence-free polymer blends

We found 10 miscible pairs of dissimilar polymers that rendered birefringence-free polymer blends. Each pair consists of a positive birefringence polymer and a negative one. Zero birefringence in the blend is attained as the result of compensation of positive and negative contributions to the overall birefringence. This concept was confirmed by characterizing the chain orientation of individual chains in the blend using polarized Fourier-transform infrared spectroscopy and by estimating theoretically the intrinsic optical anisotropy of the chain molecules.     

Phase behavior and viscoelastic properties of entangled block copolymer gels

Triblock copolymers in midblock‐selective solvents can form physical gels. However, at low triblock contents (near the percolation threshold), the bridging of chains between micelles can lead to macrophase separation. Adding a styrene–isoprene diblock to a styrene–isoprene–styrene triblock copolymer in squalane can eliminate macrophase separation, yielding a wide range of stable, single‐phase gels with a disordered arrangement of micelles. The plateau modulus of these triblock gels scales with the 2.2 power of polymer content, indicating the importance of entanglements in dictating the modulus. Comparing gels made from the midblock‐saturated derivative of the same polymer [styrene‐(ethylene‐alt‐propylene)‐styrene] in squalane reveals that the modulus differences in the gels are a direct consequence of the difference in the entanglement molecular weight of the midblock homopolymer in bulk. Finally, the broad relaxation spectrum of these triblocks is well‐described by a recent theory for the dynamics of entangled star polymers, with the breadth of the relaxation spectrum dictated by the number of entanglements per midblock in the gel. © 2001 John Wiley & Sons, Inc. J Polym Sci Part B: Polym Phys 39: 2183–2197, 2001     

Flow birefringence of concentrated polymer solutions in two-dimensional flows

The results of flow birefringence measurements are reported for polymer solutions of moderate concentration subjected to a wide range of two-dimensional flows. These flows were generated in a four-roll mill which enables one to systematically vary the ratio of the vorticity to the rate of strain in the flow while holding the velocity gradient constant. It is shown that steady-state birefringence data collected over a wide range of flow types can be correlated against the eigenvalue of the velocity gradient tensor, in agreement with criterion for strong and weak flows from model calculations. Transient birefringence measurements in which purely extensional flows were started from rest are also reported. It was observed that the birefringence went through a pronounced overshoot in time for two different polymer/solvent systems. Flow induced increases in the solution turbidity were also observed and the increased turbidity remained constant over a period of many hours after extensional flows were arrested. The birefringence, on the other hand, decayed to zero almost immediately after the flows were stopped. These changes in the turbidity suggest that crystallization of the polymer was occurring. The qualitative results of experiments are compared to recent network model calculations using the theory of Yamamoto for concentrated polymer systems. It is found that this model can predict qualitatively many of the experimental observations if the function describing the breakage of polymer chain entanglements is allowed to depend on the conformation of the polymer segments bridging the entanglements. In particular, this dependency of the entanglement breakage on the conformation of the network segments leads to a predicted overshoot of birefringence when purely extensional flows are started from rest. It is also demonstrated through this model that birefringence data taken over a wide range of flow types can be used to estimate the degree to which the network deforms affinely with the flow field.