Understanding Dynamics in Binary Mixtures of Entangled cis-1,4-Polybutadiene Melts at the Level of Primitive Path Segments by Mapping Atomistic Simulation Data onto the Tube Model

Predicting High-Density Polyethylene Melt Rheology Using a Multimode Tube Model Derived Using Non-Equilibrium Thermodynamics

Based on the Generalized bracket, or Beris-Edwards, formalism of non-equilibrium thermodynamics, we recently proposed a new differential constitutive model for the rheological study of entangled polymer melts and solutions. It amended the shortcomings of a previous model that was too strict regarding the values of the convective constraint release parameter for the model not to violate the second law of thermodynamics, and it has been shown capable of predicting a transient stress undershoot (following the overshoot) at high shear rates. In this study, we wish to further examine this model's capability to predict the rheological response of industrial polymer systems by extending it to its multiple-mode version. The comparison with industrial rheological data (High-Density Polyethylene resins), which was based on comparison with experimental data available in (a) Small Amplitude Oscillatory shear, (b) start-up shear, and (c) start-up uniaxial elongation, was noted to be good.

Structure and dynamics of a 1,4-polybutadiene melt in an alumina nanopore: A molecular dynamics simulation

We present results of Molecular Dynamics simulations of a chemically realistic model of 1,4-polybutadiene (PBD)confined in a cylindrical alumina nanopore of diameter 10 nm. The simulations are done at three different temperaturesabove the glass transition temperature T g . We investigate the density layering across the nanopore as well as theorientational ordering in the polymer melt, brought about by the confinement, on both the segmental and chain scales.For the chain scale ordering, the magnitude and orientation of the axes of the gyration tensor ellipsoid of single chainsare studied and are found to prefer to align parallel to the pore axis. Even though double bonds near the wall arepreferentially oriented along the pore walls, studying the nematic order parameter indicates that there is no nematicordering at the melt-wall interface. As for the dynamics in the melt, we focus here on the mean-square-displacement ofthe monomers for several layers across the nanopore as well as the movement of the chain center of mass which bothdisplay a slowing down of the dynamics in the layer at the wall. We also show the strong adsorption of the monomersto the pore wall at lower temperatures.

Thermodynamically consistent incorporation of entanglement spatial fluctuations in the slip-link model

We evaluate the thermodynamic consistency of the anisotropic mobile slip-link model for entangled flexible polymers. The level of description is that of a single chain, whose interactions with other chains are coarse grained to discrete entanglements. The dynamics of the model consist of the motion of entanglements through space and of the chain through the entanglements, as well as the creation and destruction of entanglements, which are implemented in a mean-field way. Entanglements are modeled as discrete slip links, whose spatial positions are confined by quadratic potentials. The confinement potentials move with the macroscopic velocity field, hence the entanglements fluctuate around purely affine motion. We allow for anisotropy of these fluctuations, described by a set of shape tensors. By casting the model in the form of the general equation for the nonequilibrium reversible-irreversible coupling from nonequilibrium thermodynamics, we show that (i) since the confinement potentials contribute to the chain free energy, they must also contribute to the stress tensor, (ii) these stress contributions are of two kinds: one related to the virtual springs connecting the slip links to the centers of the confinement potentials and the other related to the shape tensors, and (iii) these two kinds of stress contributions cancel each other if the confinement potentials become anisotropic in flow, according to a lower-convected evolution of the confinement strength or, equivalently, an upper-convected evolution of the shape tensors of the entanglement spatial fluctuations. In previous publications, we have shown that this cancellation is necessary for the model to obey the stress-optical rule and the Green-Kubo relation, and simultaneously to agree with plateau modulus predictions of multichain models and simulations.

Primitive chain network simulations for the interrupted shear response of entangled polymeric liquids

The non-linear viscoelastic response under interrupted shear flows is one of the interesting characteristics of entangled polymers. In particular, the stress overshoot in the resumed shear has been discussed concerning the recovery of the entanglement network in some studies. In this study, we performed multichain slip-link simulations to observe the molecular structure of an entangled polymer melt. After confirming the reasonable reproducibility of our simulation with the literature data, we analyzed the molecular characteristics following the decoupling approximation. We reasonably found that the segment orientation dominates the stress overshoot even under the resumed shear with minor contributions from the segment stretch and entanglement density. We defined the mitigation function for the recovery of the stress overshoot as a function of the rest time and compared it with the relaxation of the molecular quantities after the initial shear. As a result, we have found that the mitigation of the stress overshoot coincides with the relaxation of entanglement density.

Effect of Bidispersity on Dynamics of Confined Polymer Films

Using Monte Carlo simulations, we studied the effect of bidispersity on the dynamics of polymer films capped between two neutral walls, where we chose three representative compositions for bidispersed polymer films. Our results demonstrate that the characteristic entanglement length is an important parameter to clarify the effect of the bidispersity on the dynamics of polymer films. For the short chains, shorter than the characteristic entanglement length, the average number of near-neighboring particles increases with the decrease of the film thickness and limits the diffusivity of the short chains, which is independent of the film compositions. However, the dynamics of the long chains, of which is above the characteristic entanglement length, is determined by the film's composition. In our previous paper, we inferred from the structures and entanglements of the bidisperse system with short and long chains that the constraint release contributes significantly to the relaxation mechanism of long chains. By calculating the self-diffusion coefficient of long chains, we confirmed this prediction that, with a lower weight fraction of long chains, the self-diffusion coefficient of long chains decreases slowly with the decrease of the film thickness, which is similar to that of short chains. With a higher weight fraction of long chains, the competition between the disentanglement and the increased in the local degree of confinement which resulted in the self-diffusion coefficient of long chains varying non-monotonically with the film thickness. Furthermore, for the bidisperse system with long and long chains, the diffusivity of long chains was not affected by the constraint release, which varied nonmonotonically with the decrease of the film thickness due to the competition between the disentanglement and the enhanced confinement. Herein, compared with the previous work, we completely clarified the relationship between the structures and dynamics for three representative compositions of bidisperse polymer films, which contains all possible cases for bidisperse systems. Our work not only establishes a unified understanding of the dependency of dynamics on the bidispersity of polymer films, but also helps to understand the case of polydispersity, which can provide computational supports for various applications for polymer films.

Effect of Bidispersity on Structure and Entanglement of Confined Polymer Films

Using Monte Carlo simulations combined with a geometric primitive path analysis method (Z1 algorithm), we investigate the effect of bidispersity on the structure and entanglement of polymer films which consist of short (the molecular length is below the characteristic entanglement molecular length) and long (the molecular length is above the characteristic entanglement molecular length) chains between two neutral walls. Our results demonstrate the length-based migrations of chains in bidisperse films (the longer chains reside away from the walls and the shorter chains are close to the walls), which becomes more obvious with the decrease in the weight fraction of long chains. With decreasing the weight fraction of long chains, the number of short-long entanglements exhibits a dramatic increase, whereas the number of long-long entanglements exhibits a slight decrease, which indicates that short chains can significantly affect the local situations of entanglements of bidisperse polymer films. On the basis of the constraint release mechanism, our simulations imply that for the lower weight fraction of long chains, the local degree of confinement instead of the long-long entanglements has a marked effect on the relaxation of long chains, due to the fast relaxation of short chains dilating the tube diameter of long chains. However, for the higher weight fraction of long chains, after the relaxation of short chains, the long-long entanglements are in sufficient quantities to restrict long chains within a tube, which implies that the relaxation of long chains is hardly affected by the number of short-long entanglements. Our work can be helpful for understanding the microscopic structure and entanglement of bidisperse polymer films, which can provide computational support for their various technological applications.