High-fidelity scaling relationships for determining dissipative particle dynamics parameters from atomistic molecular dynamics simulations of polymeric liquids

Atomistic Simulation of Flow-Induced Microphase Separation and Crystallization of an Entangled Polyethylene Melt Undergoing Uniaxial Elongational Flow and the Role of Kuhn Segment Extension

Atomistic simulations of the linear, entangled polyethylene C₁₀₀₀H₂₀₀₂ melt undergoing steady-state and startup conditions of uniaxial elongational flow (UEF) over a wide range of flow strength were performed using a united-atom model for the atomic interactions between the methylene groups constituting the polymer macromolecules. Rheological, topological, and microstructural properties of these nonequilibrium viscoelastic materials were computed as functions of strain rate, focusing on regions of flow strength where flow-induced phase separation and flow-induced crystallization were evident. Results of the UEF simulations were compared with those of prior simulations of planar elongational flow, which revealed that uniaxial and planar flows exhibited essentially a universal behavior, although over strain rate ranges that were not completely equivalent. At intermediate flow strength, a purely configurational microphase separation was evident that manifested as a bicontinuous phase composed of regions of highly stretched molecules that enmeshed spheroidal domains of relatively coiled chains. At high flow strength, a flow-induced crystallization (FIC) occurred, producing a semicrystalline material possessing a high degree of crystallinity and primarily a monoclinic lattice structure. This FIC phase formed at a temperature (450 K) high above the quiescent melting point (≈400 K) and remained stable after cessation of flow for temperature at or below 435 K. Careful examination of the Kuhn segments constituting the polymer chains revealed that the FIC phase only formed once the Kuhn segments had become essentially fully extended under the UEF flow field. Thermodynamic properties such as the heat of fusion and heat capacity were estimated from the simulations and found to compare favorably with experimental values.

Microstructural evolution and reverse flow in shear-banding of entangled polymer melts

The temporal and spatial evolution of shear banding under startup of shear flow was simulated for highly entangled, linear, monodisperse polyethylene melts of differing molecular weight, C₇₅₀H₁₅₀₂, C₁₂₀₀H₂₄₀₂, and C₃₀₀₀H₆₀₀₂, using a high-fidelity coarse-grained dissipative particle dynamics method. It was determined that shear stress was dominated by segmental orientation of entangled strands at low shear rates, but at a critical shear rate below the reciprocal of the Rouse time, flow-induced disentanglement resulted in the onset of chain tumbling that reduced the average degree of orientation, leading to a regime of decreasing shear stress, with a commensurate onset of increasing average chain extension imposed by the strong flow kinematics that ultimately drove the steady-state shear stress higher. During startup of shear flow, shear band development began immediately after the maximum in the first normal stress difference, where distinct fast and slow bands formed. The slow bands consisted of relatively entangled and coiled molecules, whereas the fast bands consisted of more disentangled and extended chains that experienced quasiperiodic rotation/retraction cycles. The simulation results often exhibited a generation of temporary reverse flow, in which the local fluid velocity was temporarily opposite to that of the bulk flow direction, at the onset of the shear-banding phenomena; this effect was consistent with earlier experiments and theoretical results. The physical mechanism for the generation of reverse flow during shear-band formation was investigated and found to be related to the recoil of the molecules comprising the slow band. Overall, the phenomenon of shear banding appeared to arise due to flow-induced disentanglement from orientational ordering and segmental stretching that affected individual chains to different degrees, ultimately resulting in regions of relatively coiled and entangled chains that evolved into a slow band, whereas the locally disentangled chains, experiencing quasiperiodic stretch-rotation cycles, formed a fast band. The transitional period resulted in a kinematic instability that generated the temporary reverse-flow phenomenon.

Multiscale simulations of polyzwitterions in aqueous bulk solutions and brush array configurations

Zwitterionic polymers are very promising candidates for antifouling materials that exhibit high chemical stability as compared to polyethylene glycol-based systems. A number of simulation and experimental studies have emerged over recent years for the investigation of sulfobetaine-based zwitterionic polymers. Investigating the structural and thermodynamic properties of such polymers requires access to broad time and length regimes, thus necessitating the development of multiscale simulation strategies. The present article advocates a mesoscopic dissipative particle dynamics (DPD) model capable of addressing a wide range of time and length scales. The mesoscopic force field was developed hand-in-hand with atomistic simulations based on the OPLS force field through a bottom-up parameterization procedure that matches the atomistically calculated strand-length, strand-angle and pair distribution functions. The DPD model is validated against atomistic simulations conducted in this work, and against relevant atomistic simulation studies, theoretical predictions and experimental correlations from the literature. Properties examined include the conformations of SPE polymers in dilute bulk aqueous solution, the density profile and thickness of brush arrays as functions of the grafting density and chain length. In addition, we compute the potential of mean force of an approaching hydrophilic or hydrophobic foulant via umbrella sampling as a function of its position relative to the poly-zwitterion-covered surface. The aforementioned observables lead to important insights regarding the conformational tendencies of grafted polyzwitterions and their antifouling properties.

Dissipative particle dynamics simulations in colloid and Interface science: a review

Dissipative particle dynamics (DPD) is one of the most efficient mesoscale coarse-grained methodologies for modeling soft matter systems. Here, we comprehensively review the progress in theoretical formulations, parametrization strategies, and applications of DPD over the last two decades. DPD bridges the gap between the microscopic atomistic and macroscopic continuum length and time scales. Numerous efforts have been performed to improve the computational efficiency and to develop advanced versions and modifications of the original DPD framework. The progress in the parametrization techniques that can reproduce the engineering properties of experimental systems attracted a lot of interest from the industrial community longing to use DPD to characterize, help design and optimize the practical products. While there are still areas for improvements, DPD has been efficiently applied to numerous colloidal and interfacial phenomena involving phase separations, self-assembly, and transport in polymeric, surfactant, nanoparticle, and biomolecules systems.

Effects of chain length and polydispersity on shear banding in simple shear flow of polymeric melts

The characteristics of shear banding were investigated in entangled, polydisperse, linear polymer melts under steady-state and startup conditions of simple shear flow. This virtual experimentation was conducted using course-grained nonequilibrium dissipative particle dynamics simulations expressed in terms of a force-field representation that faithfully models the atomistic system dynamics. We examined melts with two mean molecular bead numbers of Nn = 2 50 and 400 and polydispersity indexes of 1.0, 1.025, and 1.05. The wide range of relaxation timescales in the polydisperse melts decreased the nonmonotonic character of the steady-state shear stress vs. shear rate profile compared to a monodisperse linear melt. The polydispersity level required to observe a stress plateau in the shear stress profile at intermediate shear rates was correlated with the nominal entanglement density. Startup of shear flow simulations revealed the development of spatial inhomogeneities and dynamic instabilities in polydisperse fluids containing both monotonic and nonmonotonic shear stress flow curves. Although the shape and duration of instabilities were found to be correlated with the monotonicity of the shear stress profile, the onset and underlying mechanism leading to the formation of shear bands were generally universal. The simulations revealed that perturbations arose soon after the occurrence of a large stress overshoot under startup conditions, and that banded structures stemmed from local reorientation and subsequent deconstruction of the entanglement network. Furthermore, data indicated that the inception of strain localization occurred at shear rates near the reciprocal of the Rouse characteristic timescale, [small gamma, Greek, dot above] > τR-1. Transient shear banding was observed in shorter chain melts undergoing startup of shear flow in which instabilities arose after the appearance of a stress overshoot. These instabilities eventually decayed, but only long after the stresses had attained their steady-state values. The longer chain melt exhibited a shear band structure that remained indefinitely after the stresses had attained steady state.