skip to main content


Title: Elucidating the role of network topology dynamics on the coil-stretch transition hysteresis in extensional flow of entangled polymer melts

Dissipative particle dynamics (DPD) simulations are performed on coarse-grained replicas of linear, monodisperse entangled polyethylene melts [Formula: see text] and [Formula: see text] undergoing both steady-state and transient planar elongational flow (PEF). The fidelity of the DPD simulations is verified by direct comparison of flow topological and rheological properties of a 334-particle chain liquid against the united-atom [Formula: see text] liquid, simulated using nonequilibrium molecular dynamics (NEMD). These DPD simulations demonstrate that a flow-induced coil-stretch transition (CST) and its associated hysteresis caused by configurational microphase separation, as observed in previous NEMD simulations of PEF, can be replicated using a more computationally efficient coarse-grained system. Results indicate that the breadth of the CST hysteresis loop is enlarged for the longer molecule liquid relative to the shorter one. Furthermore, relaxation simulations reveal that reducing the applied flow Deborah number ([Formula: see text]) from a high value corresponding to a homogeneous phase of highly stretched molecules to a [Formula: see text] within the biphasic region results in a two-stage relaxation process. There is a fast initial stratification of the kinetically trapped highly stretched chains into regions of highly extended and less extended chains, which displays similar behavior to a system undergoing a spinodal decomposition caused by spatial configurational free energy fluctuations. After a short induction period of apparently random duration, the less extended chain regions experience a stochastic nucleation event that induces configurational relaxation to domains composed of coiled molecules over a much longer time scale, leaving the more highly extended chains in surrounding sheetlike domains. The time scales of these two relaxation processes are of the same order of magnitude as the Rouse and disengagement times of the equilibrium liquids.

 
more » « less
NSF-PAR ID:
10364636
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
Society of Rheology
Date Published:
Journal Name:
Journal of Rheology
Volume:
66
Issue:
3
ISSN:
0148-6055
Page Range / eLocation ID:
p. 551-569
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    Ternary block copolymer (BCP)‐homopolymer (HP) blends offer a simple method for tuning nanostructure sizes to meet application‐specific demands. Comprehensive dissipative particle dynamic (DPD) simulations were performed to study the impact of polymer interactions, molecular weight, and HP volume fraction (φHP) on symmetric ternary blend morphological stability and domain spacing. DPD reproduces key features of the experimental phase diagram, including lamellar domain swelling with increasingφHP, the formation of an asymmetric bicontinuous microemulsion at a critical HP concentration , and macrophase separation with further HP addition. Simulation results matched experimental values for and lamellar swelling as a function of HP to BCP chain length ratio,α = NHP/NBCP. Structural analysis of blends with fixedφHPbut varyingαconfirmed that ternary blends follow the wet/dry brush model of domain swelling with the miscibility of HPs and BCPs depending onα. Longer HPs concentrate in the center of domains, boosting their swelling efficiencies compared to shorter chains. These results advance our understanding of BCP‐HP blend phase behavior and demonstrate the value of DPD for studying polymeric blends. © 2019 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys.2019,57, 794–803

     
    more » « less
  2. Herein, we have employed a supervised learning approach combined with Core-Modified Dissipative Particle Dynamics Simulations (CM-DPD) in order to develop and design a reliable physics-based computational model that will be used in studying confined flow of suspensions. CM-DPD was recently developed and has shown promising performance in capturing rheological behavior of colloidal suspensions; however, the model becomes problematic when the flow of the material is confined between two walls. Wall-penetration by the particles is an unphysical phenomenon that occurs in coarse-grained simulations such as Dissipative Particle Dynamics (DPD) that mostly rely on soft inter-particle interactions. Different solutions to this problem have been proposed in the literature; however, no reports have been given on how to deal with walls using CM-DPD. Due to complexity of interactions and system parameters, designing a realistic simulation model is not a trivial task. Therefore, in this work we have trained a Random Forest (RF) for predicting wall penetration as we vary input parameters such as interaction potentials, flow rate, volume fraction of colloidal particles, and confinement ratio. The RF predictions were compared against simulation tests, and a sufficiently high accuracy and low errors were obtained. This study shows the viability and potentiality of ML combined with DPD to perform parametric studies in complex fluids. 
    more » « less
  3. Nonequilibrium work relations have fundamentally advanced our understanding of molecular processes. In recent years, fluctuation theorems have been extensively applied to understand transitions between equilibrium steady-states, commonly described by simple control parameters such as molecular extension of a protein or polymer chain stretched by an external force in a quiescent fluid. Despite recent progress, far less is understood regarding the application of fluctuation theorems to processes involving nonequilibrium steady-states such as those described by polymer stretching dynamics in nonequilibrium fluid flows. In this work, we apply the Crooks fluctuation theorem to understand the nonequilibrium thermodynamics of dilute polymer solutions in flow. We directly determine the nonequilibrium free energy for single polymer molecules in flow using a combination of single molecule experiments and Brownian dynamics simulations. We further develop a time-dependent extensional flow protocol that allows for probing viscoelastic hysteresis over a wide range of flow strengths. Using this framework, we define quantities that uniquely characterize the coil-stretch transition for polymer chains in flow. Overall, generalized fluctuation theorems provide a powerful framework to understand polymer dynamics under far-from-equilibrium conditions. 
    more » « less
  4. Abstract

    Conjugated polymers (CPs), characterized by rigid conjugation backbones and flexible peripheral side chains, hold significant promise in various organic optoelectronic applications. In this study, we employ coarse‐grained molecular dynamics (CG‐MD) simulations to investigate the intricate interplay of solvent quality, temperature, and chain architecture (e.g., side‐chain length and molecular mass) on the conformational behaviors of CPs in dilute solutions. Our research uncovers distinctive conformational behaviors under varying solvent conditions, highlighting the versatile nature of polymer chains, which can adopt extended configurations in good solvents and transition to aggregated states in poor solvents. Additionally, the mass scaling exponent , a robust structural descriptor, consistently described CPs behavior across diverse architectures and solvent conditions. Furthermore, our study shows that a CP with longer side‐chain exhibits improved solubility, which is further confirmed by experimental observations. Moreover, our analysis of the shape descriptor provided valuable insights into the symmetry and dimensionality of CPs under varying solvent conditions. These findings offer a comprehensive understanding of conformational behaviors of CPs in dilute solution, which are helpful in guiding the conformational design of polymer for specific applications.

     
    more » « less
  5. Monodisperse suspensions of Brownian colloidal spheres crystallize at high densities, and ordering under shear has been observed at densities below the crystallization threshold. We perform large-scale simulations of a model suspension containing over [Formula: see text] particles to quantitatively study the ordering under shear and to investigate its link to the rheological properties of the suspension. We find that at high rates, for [Formula: see text], the shear flow induces an ordering transition that significantly decreases the measured viscosity. This ordering is analyzed in terms of the development of layering and planar order, and we determine that particles are packed into hexagonal crystal layers (with numerous defects) that slide past each other. By computing local [Formula: see text] and [Formula: see text] order parameters, we determine that the defects correspond to chains of particles in a squarelike lattice. We compute the individual particle contributions to the stress tensor and discover that the largest contributors to the shear stress are primarily located in these lower density, defect regions. The defect structure enables the formation of compressed chains of particles to resist the shear, but these chains are transient and short-lived. The inclusion of a contact friction force allows the stress-bearing structures to grow into a system-spanning network, thereby disrupting the order and drastically increasing the suspension viscosity. 
    more » « less