Influence of molecular architecture on the entanglement network: topological analysis of linear, long- and short-chain branched polyethylene melts via Monte Carlo simulations

Water Dissolved in a Variety of Polymers Studied by Molecular Dynamics Simulation and a Theory of Solutions

The performance of a polymer medium as a separation membrane is determined by the dissolution free energy ΔG and diffusion coefficient D of the permeant. In this work, ΔG and D of water are investigated with all-atom molecular dynamics simulation in a wide variety of polymer species in the amorphous state. The computed ΔG is shown to agree well with the experimental value for linear homopolymers, and the degrees of polymerization of the homopolymers do not affect ΔG when they are beyond ∼10. The copolymers of ethylene-vinylidene difluoride, ethylene-vinyl acetate, and ethylene-acrylamide are then examined by changing the repeating patterns of the constituent monomers in both the periodic and graft forms. It is found that ΔG is determined primarily by the overall compositions of the monomers and is not affected by the copolymerization topology (periodic or graft). The hydrophobicity of the copolymer is enhanced, furthermore, when the hydrophobicity and hydrophilicity of the ethylene and non-ethylene parts are well contrasted and those parts are fragmented along the polymer chain. According to the computed D, the diffusivity of water tends to be larger when the (co)polymer is more hydrophobic and ΔG is more positive. D is actually seen to vary by orders of magnitude with the polymer structures, while the effect of the polymer species on the water permeation is stronger for ΔG than for D.

A generalized tube model of rubber elasticity

In the present paper, a new type of micro-mechanically motivated chain network model for rubber-like materials is proposed. The model captures topological constraints of polymer network chains, in particular, entanglements. The model demonstrates how the local molecular packing constraints modify under deformation and shows the impact of these changes on the macroscopic elasticity of the material. To this end, we combine concepts of a confining tube and a slip-link (reptation) model. In these models, entanglements of polymer chains play an important role. The nature of entanglements is discussed, and relationships governing entanglements are formulated in terms of molecular physics. In the context of nonlinear elasticity, we apply a non-affine concept which captures the liquid-like behavior of polymer networks at smaller scales in a more realistic way. Model predictions show good agreement with experimental results from uniaxial and biaxial tension tests.

Nonequilibrium Monte Carlo simulations of entangled polymer melts under steady shear flow

We present a nonequilibrium Monte Carlo (MC) methodology based on expanded nonequilibrium thermodynamic formalism to simulate entangled polymeric materials undergoing steady shear flow. Motivated by the standard kinetic theory for entangled polymers, a second-rank symmetric conformation tensor based on the entanglement segment vector was adopted in the expanded statistical-ensemble based MC method as the nonequilibrium structural variable that properly represents the deformed structure of the system by the flow. The corresponding (second-rank symmetric) conjugate thermodynamic force variable was introduced to simulate an external flow field. As a test case, we applied the GENERIC MC to C₄₀₀H₈₀₂ entangled linear polyethylene melts in a wide range of shear rates. Detailed analysis of the GENERIC MC results for various structural and rheological properties (such as chain size, chain orientation, mesoscale chain configuration, topological properties, and material functions) was carried out via direct comparison with the corresponding NEMD results. Overall, the GENERIC MC is shown to predict the general trends of the nonequilibrium properties of polymer systems reasonably for a wide range of flow strengths (i.e., 0.5 ≤ De ≤ 540). In conjunction with NEMD, the present MC method can thus be used to extract fundamental thermodynamic information (nonequilibrium entropy and free energy functions) of entangled polymeric systems under various flow types. Despite the general consistency between the GENERIC MC and NEMD results, some quantitative discrepancies appear in the intermediate-to-strong flow regime. This behavior stems from the inherent mean-field nature of the MC force field which is applied to the individual entanglement segments independently and uniformly along the chain, without accounting for any flow-induced structural or dynamical correlations between segments.

Controlling crystal polymorphism of isotactic poly(1-butene) by incorporating long chain branches

Isotactic poly(1-butene) (iPB-1) is a high performance plastic with outstanding properties, such as flexibility, superior creep, environmental stress cracking and abrasive resistance. However, it exhibits a complex crystal polymorphism and polymorphic transformation behavior, which has limited its commercial development. In this paper, the incorporation of long chain branches (LCBs) causes coil contraction in the melt, which favors the direct melt-crystallization of form III that was generally crystallized from solutions and made of unconventional highly twined lamellae. Consequently, low-to-moderately branched iPB-1 samples as-crystallize from the melt into mixtures of form II and form III by compression-molding and fast cooling of the melt to room temperature, and the fraction of crystals of form III (fIII) increases with increasing concentration of LCBs, whereas highly branched samples can as-crystallize into pure form III with uniform crystal size distribution. The corresponding thermomechanical properties can be modified by controlling fIII.

Molecular characteristics of stress overshoot for polymer melts under start-up shear flow

Stress overshoot is one of the most important nonlinear rheological phenomena exhibited by polymeric liquids undergoing start-up shear at sufficient flow strengths. Despite considerable previous research, the fundamental molecular characteristics underlying stress overshoot remain unknown. Here, we analyze the intrinsic molecular mechanisms behind the overshoot phenomenon using atomistic nonequilibrium molecular dynamics simulations of entangled linear polyethylene melts under shear flow. Through a detailed analysis of the transient rotational chain dynamics, we identify an intermolecular collision angular regime in the vicinity of the chain orientation angle θ ≈ 20° with respect to the flow direction. The shear stress overshoot occurs via strong intermolecular collisions between chains in the collision regime at θ = 15°-25°, corresponding to a peak strain of 2-4, which is an experimentally well-known value. The normal stress overshoot appears at approximately θ = 10°, at a corresponding peak strain roughly equivalent to twice that for the shear stress. We provide plausible answers to several basic questions regarding the stress overshoot, which may further help understand other nonlinear phenomena of polymeric systems.

Molecular dynamics for linear polymer melts in bulk and confined systems under shear flow

In this work, we analyzed the individual chain dynamics for linear polymer melts under shear flow for bulk and confined systems using atomistic nonequilibrium molecular dynamics simulations of unentangled (C₅₀H₁₀₂) and slightly entangled (C₁₇₈H₃₅₈) polyethylene melts. While a certain similarity appears for the bulk and confined systems for the dynamic mechanisms of polymer chains in response to the imposed flow field, the interfacial chain dynamics near the boundary solid walls in the confined system are significantly different from the corresponding bulk chain dynamics. Detailed molecular-level analysis of the individual chain motions in a wide range of flow strengths are carried out to characterize the intrinsic molecular mechanisms of the bulk and interfacial chains in three flow regimes (weak, intermediate, and strong). These mechanisms essentially underlie various macroscopic structural and rheological properties of polymer systems, such as the mean-square chain end-to-end distance, probability distribution of the chain end-to-end distance, viscosity, and the first normal stress coefficient. Further analysis based on the mesoscopic Brightness method provides additional structural information about the polymer chains in association with their molecular mechanisms.