Multiscale modeling of polyisoprene on graphite.

The local dynamics and the conformational properties of polyisoprene next to a smooth graphite surface constructed by graphene layers are studied by a multiscale methodology. First, fully atomistic molecular dynamics simulations of oligomers next to the surface are performed. Subsequently, Monte Carlo simulations of a systematically derived coarse-grained model generate numerous uncorrelated structures for polymer systems. A new reverse backmapping strategy is presented that reintroduces atomistic detail. Finally, multiple extensive fully atomistic simulations with large systems of long macromolecules are employed to examine local dynamics in proximity to graphite. Polyisoprene repeat units arrange close to a parallel configuration with chains exhibiting a distribution of contact lengths. Efficient Monte Carlo algorithms with the coarse-grain model are capable of sampling these distributions for any molecular weight in quantitative agreement with predictions from atomistic models. Furthermore, molecular dynamics simulations with well-equilibrated systems at all length-scales support an increased dynamic heterogeneity that is emerging from both intermolecular interactions with the flat surface and intramolecular cooperativity. This study provides a detailed comprehensive picture of polyisoprene on a flat surface and consists of an effort to characterize such systems in atomistic detail.

Detailed atomistic Monte Carlo simulations of a polymer melt on a solid surface and around a nanoparticle.

The molecular factors that govern interfacial interactions between a polymer melt and a solid surface remain largely unclear despite significant progress made in the last years. Simulations are increasingly employed to elucidate these features, however, equilibration and sampling with models of long macromolecules in such heterogeneous systems present significant challenges. In this study, we couple the application of preferential sampling techniques with connectivity-altering Monte Carlo algorithms to explore the configurational characteristics of a polyethylene melt in proximity to a surface and a highly curved nanoparticle. Designed algorithms allow efficient sampling at all length scales of large systems required to avoid finite-size effects. Using detailed atomistic models for the polymer and realistic structures for a silica surface and a fullerene, we find that at the extreme limit where particles are comparable to the polymer Kuhn segment length, curvature penalizes the formation of long train segments. As a result, an increased number of shorter contacts belonging to different chains are made competing with the anticipated decrease of the bound layer thickness with particle size if polymer adsorbed per unit area remained constant. For very small nanoparticles, formation of new train segments cannot compete with the overall reduction of adsorbance which is present irrespective of the enthalpic interactions; a result that demonstrates the need for an accurate description of polymer rigidity at these length scales.