No Arabic abstract
Drawing an analogy to the paradigm of quasi-elastic neutron scattering, we present a general approach for quantitatively investigating the spatiotemporal dependence of structural anisotropy relaxation in deformed polymers by using small-angle neutron scattering. Experiments and non-equilibrium molecular dynamics simulations on polymer melts over a wide range of molecular weights reveal that their conformational relaxation at relatively high momentum transfer $Q$ and short time can be described by a simple scaling law, with the relaxation rate proportional to $Q$. This peculiar scaling behavior, which cannot be derived from the classical Rouse and tube models, is indicative of a surprisingly weak direct influence of entanglement on the microscopic mechanism of single-chain anisotropy relaxation.
The flow and deformation of macromolecules is ubiquitous in nature and industry, and an understanding of this phenomenon at both macroscopic and microscopic length scales is of fundamental and practical importance. Here we present the formulation of a general mathematical framework, which could be used to extract, from scattering experiments, the molecular relaxation of deformed polymers. By combining and modestly extending several key conceptual ingredients in the literature, we show how the anisotropic single-chain structure factor can be decomposed by spherical harmonics and experimentally reconstructed from its cross sections on the scattering planes. The resulting wavenumber-dependent expansion coefficients constitute a characteristic fingerprint of the macromolecular deformation, permitting detailed examinations of polymer dynamics at the microscopic level. We apply this approach to survey a long-standing problem in polymer physics regarding the molecular relaxation in entangled polymers after a large step deformation. The classical tube theory of Doi and Edwards predicts a fast chain retraction process immediately after the deformation, followed by a slow orientation relaxation through the reptation mechanism. This chain retraction hypothesis, which is the keystone of the tube theory for macromolecular flow and deformation, was critically examined by analyzing the fine features of the two-dimensional anisotropic spectra from small-angle neutron scattering by entangled polystyrenes. It is shown that the unique scattering patterns associated with the chain retraction mechanism were not experimentally observed. This result calls for a fundamental revision of the current theoretical picture for nonlinear rheological behavior of entangled polymeric liquids.
Evolving structure and rheology across Kuhn scale interfaces in entangled polymer fluids under flow play a prominent role in processing of manufactured plastics, and have numerous other applications. Quantitative tracking of chain conformation statistics on the Kuhn scale is essential for developing computational models of such phenomena. For this purpose, we formulate here a two-scale/two-mode model of entangled polymer chains under flow. Each chain is partitioned by successive entanglements into strands that are in one of two modes: entangled or dangling. On the strand scale, conformation statistics of ideal (non-interacting) strands follows a differential evolution equation for the second moment of its end-to-end distance. The latter regulates persistent random walks sampling conformation statistics of ideal entangled strands on the Kuhn scale, as follows from a generalized Green-Kubo relation and the Maximum Entropy Principle. We test it numerically for a range of deformation rates at the start-up of simple elongational and shear flows. A self-consistent potential, representing segmental interactions, modifies strand conformation statistics on the Kuhn scale, as it renormalizes the parameters controlling the persistent random walk. The generalized Green-Kubo relation is then inverted to determine how the second moment of the strand end-to-end distance is changed by the self-consistent potential. This allows us to devise a two-scale propagation scheme for the statistical weights of subchains of the entangled chain. The latter is used to calculate local volume fractions for each chemical type of Kuhn segments in entangled chains, thus determining the self-consistent potential.
We present a simple reaction kinetics model to describe the polymer synthesis used by Lusignan et al. (PRE, 60, 5657, 1999) to produce randomly branched polymers in the vulcanization class. Numerical solution of the rate equations gives probabilities for different connections in the final product, which we use to generate a numerical ensemble of representative molecules. All structural quantities probed by Lusignan et al. are in quantitative agreement with our results for the entire range of molecular weights considered. However, with detailed topological information available in our calculations, our estimate of the `rheologically relevant linear segment length is smaller than that estimated by them. We use a numerical method based on tube model of polymer melts to calculate the rheological properties of such molecules. Results are in good agreement with experiment, except that in the case of the largest molecular weight samples our estimate of the zero-shear viscosity is significantly lower than the experimental findings. Using acid concentration as an indicator for closeness to the gelation transition, we show that the high-molecular-weight polymers considered are at the limit of mean-field behavior - which possibly is the reason for this disagreement. For a truly mean-field gelation class of model polymers, we numerically calculate the rheological properties for a range of segment lengths. Our calculations show that the tube theory with dynamical dilation predicts that, very close to the gelation limit, contribution to viscosity for this class of polymers is dominated by the contribution from constraint-release Rouse motion and the final viscosity exponent approaches Rouse-like value.
The interplay of nematic order and phase separation in solutions of semiflexible polymers in solvents of variable quality is investigated by density functional theory (DFT) and molecular dynamics (MD) simulations. We studied coarse-grained models, with a bond-angle potential to control chain stiffness, for chain lengths comparable to the persistence length of the chains. We varied both the density of the monomeric units and the effective temperature that controls the quality of the implicit solvent. For very stiff chains only a single transition from an isotropic fluid to a nematic is found, with a phase diagram of swan-neck topology. For less stiff chains, however, also unmixing between isotropic fluids of different concentration, ending in a critical point, occurs for temperatures above a triple point. The associated critical behavior is examined in the MD simulations and found compatible with Ising universality. Apart from this critical behavior, DFT calculations agree qualitatively with the MD simulations.
Using a coarse-grained bead-spring model for semi-flexible macromolecules forming a polymer brush, structure and dynamics of the polymers is investigated, varying chain stiffness and grafting density. The anchoring condition for the grafted chains is chosen such that their first bonds are oriented along the normal to the substrate plane. Compression of such a semi-flexible brush by a planar piston is observed to be a two-stage process: for small compressions the chains contract by buckling deformation whereas for larger compression the chains exhibit a collective (almost uniform) bending deformation. Thus, the stiff polymer brush undergoes a 2-nd order phase transition of collective bond reorientation. The pressure, required to keep the stiff brush at a given degree of compression, is thereby significantly smaller than for an otherwise identical brush made of entirely flexible polymer chains! While both the brush height and the chain linear dimension in the z-direction perpendicular to the substrate increase monotonically with increasing chain stiffness, lateral (xy) chain linear dimensions exhibit a maximum at intermediate chain stiffness. Increasing the grafting density leads to a strong decrease of these lateral dimensions, compatible with an exponential decay. Also the recovery kinetics after removal of the compressing piston is studied, and found to follow a power-law / exponential decay with time. A simple mean-field theoretical consideration, accounting for the buckling/bending behavior of semi-flexible polymer brushes under compression, is suggested.