Multiscale Simulations of Graphite-Capped Polyethylene Melts: Brownian Dynamics/Kinetic Monte Carlo Compared to Atomistic Calculations and Experiment

AP Sgouros and GG Vogiatzis and G Megariotis and C Tzoumanekas and DN Theodorou, MACROMOLECULES, 52, 7503-7523 (2019).

DOI: 10.1021/acs.macromol.9b01379

We propose a multiscale simulation strategy capable of investigating and predicting the interfacial structure and dynamics of polymer/solid interfaces over a wide range of time and length scales. The lowest level of description involves atomistic molecular dynamics simulations providing useful information for the parametrization of the higher level of description, a mesoscopic hybrid particle/field slip-spring model Macromolecules 2017, 50, 3004, 4524. The mesoscopic model is extended so as to capture the interfacial properties of polymer/ solid interfaces observed in atomistic simulations and to make predictions for long polymer chains exhibiting relaxation phenomena over times inaccessible by conventional atomistic simulations. As a case study, our multiscale approach is applied to polyethylene/graphite interfaces; our approach is nevertheless generic and can be implemented for various other combinations of polymers and solid surfaces. The conformational properties of the mesoscopic chains in the vicinity of the interface are resolved through the incorporation of springs with variable, walldistance-dependent stiffness. The interfacial kinetics is described in terms of adsorption and desorption processes of the mesoscopic beads, which are classified as adsorbed and free in relation to solid substrates; these processes are tracked by the kinetic Monte Carlo in the course of the mesoscopic simulation, with the relevant rate constants derived by hazard plot analysis of the atomistic molecular dynamics. For the classification of the polymer segments we develop a new mesoscopic definition for the adsorbed state, taking into account the position, energy, and previous history of the polymer beads, which yields results that are insensitive to the time-step of the simulation. The residence time of the macromolecules in the adsorbed state is quantified through various approaches (rate constants from cumulative hazard slopes and monitoring of detachment times); overall, the characteristic detachment time of long polymer chains adsorbed from the melt is found to increase with molar mass, M, as similar to M-2.5 while the relaxation time from the autocorrelation of the end-to-end vector increases as similar to M-3.5; the latter scaling has been observed for bulk samples also. The impact of the time-step and of diffusion-driven recrossings between adsorbed and free states of beads is assessed through an analytical solution to a simple model for "adsorption" near repulsive walls.

Return to Publications page