Hierarchical modeling of polystyrene: From atomistic to coarse-grained simulations

We present a hierarchical approach that combines atomistic and mesoscopic simulations that can generally be applied to vinyl polymers. As a test case, the approach is applied to atactic polystyrene (PS). First, a specific model for atactic PS is chosen. The bonded parameters in the coarse-grained force field, based on data obtained from atomistic simulations of isolated PS dimers, are chosen in a way which allows to differentiate between meso and racemic dyads. This approach in principle allows to study isotactic and syndiotactic melts as well. Nonbonded interactions between coarse-grained beads were chosen as purely repulsive. The proposed mesoscopic model reproduces both the local structure and the chain dimensions properly. An explicit time mapping is performed, based on the atomistic and CG mean-square displacements of short chains, demonstrating an effective speed up of about 3 orders of magnitude compared to brute force atomistic simulations. Finally the equilibrated coarse-grained chains are back mapped onto the atomistic systems. This opens new routes for obtaining well equilibrated high molecular weight polymeric systems and also providing very long dynamic trajectories at the atomistic level for these polymers.