Reverse mapping method for complex polymer systems

We present a comprehensive approach for reverse mapping of complex polymer systems in which the connectivity is created by the simulation of chemical reactions at the coarse‐grained scale. Within the work, we use a recently developed generic adaptive reverse mapping procedure that we adapt to handle the varying connectivity structure resulting from the chemical reactions. The method is independent of the coarse‐grained and fine‐grained force‐fields by design and relies only on a single control parameter. We employ the approach to reverse map four different systems: a three‐component epoxy network, a trimethylol melamine network, a hyperbranched polymer, and a polyethylene terephthalate. In the case of the epoxy network, we use two fine‐scale representations: fully atomistic and united‐atom. Whereas the trimethylol melamine network the hyperbranched polymer and the polyethylene terephthalate are reverse mapped to the fully atomistic description. After the reverse mapping, we examine the fine‐grained structure by comparing the radial distribution functions with respect to the control parameter. Moreover, in the case of the epoxy we perform tensile‐test experiments and examine the resulting Young's modulus. In all cases, we show how the properties of the reverse mapped systems depend on the control parameters. In general, we see that the results are relatively insensitive to the control parameter and the resulting atomistic systems are stable. Only for the trimethylol melamine network, we notice chemically incorrect conformations when the reverse mapping is performed too fast. We provide a remedy for this issue. © 2017 Wiley Periodicals, Inc.

[1]  John D. Hunter,et al.  Matplotlib: A 2D Graphics Environment , 2007, Computing in Science & Engineering.

[2]  P. Khalatur,et al.  Highly Cross-Linked Epoxy Resins: An Atomistic Molecular Dynamics Simulation Combined with a Mapping/Reverse Mapping Procedure , 2007 .

[3]  J. Clarke,et al.  Molecular dynamics simulation of an amorphous polymer under tension. 1. Phenomenology , 1991 .

[4]  Janet E. Jones On the determination of molecular fields. —II. From the equation of state of a gas , 1924 .

[5]  Alexander Lukyanov,et al.  Versatile Object-Oriented Toolkit for Coarse-Graining Applications. , 2009, Journal of chemical theory and computation.

[6]  S. Godtfredsen,et al.  Ullmann ' s Encyclopedia of Industrial Chemistry , 2017 .

[7]  Dirk Reith,et al.  Deriving effective mesoscale potentials from atomistic simulations , 2002, J. Comput. Chem..

[8]  F. Müller-Plathe,et al.  A coarse-grained molecular dynamics – reactive Monte Carlo approach to simulate hyperbranched polycondensation , 2014 .

[9]  S. Crescitelli,et al.  A multiscale approach to triglycerides simulations: from atomistic to coarse-grained models and back. , 2012, Faraday discussions.

[10]  Formation of the Interphase of a Cured Epoxy Resin Near a Metal Surface: Reactive Coarse-Grained Molecular Dynamics Simulations , 2014 .

[11]  G. Odegard,et al.  Molecular modeling of crosslink distribution in epoxy polymers , 2012 .

[12]  V. Varshney,et al.  A Molecular Dynamics Study of Epoxy-Based Networks : Cross-Linking Procedure and Prediction of Molecular and Material Properties , 2008 .

[13]  M. Cho,et al.  The glass transition and thermoelastic behavior of epoxy-based nanocomposites: A molecular dynamics study , 2011 .

[14]  Brian E. Granger,et al.  IPython: A System for Interactive Scientific Computing , 2007, Computing in Science & Engineering.

[15]  H. Berendsen,et al.  Efficient Algorithms for Langevin and DPD Dynamics. , 2012, Journal of chemical theory and computation.

[16]  Pep Español,et al.  Statistical Mechanics of Coarse-Graining , 2004 .

[17]  José Mario Martínez,et al.  PACKMOL: A package for building initial configurations for molecular dynamics simulations , 2009, J. Comput. Chem..

[18]  Z. Zhang,et al.  Multiscale Simulation Study on the Curing Reaction and the Network Structure in a Typical Epoxy System , 2011 .

[19]  Multi-scale modeling of polyimides , 2004 .

[20]  Weijian Xu,et al.  Atomistic molecular modelling of crosslinked epoxy resin , 2006 .

[21]  D. Yan,et al.  Hyperbranched polymers: from synthesis to applications , 2004 .

[22]  Pavel G. Khalatur,et al.  Thermal Properties and Topology of Epoxy Networks: A Multiscale Simulation Methodology , 2015 .

[23]  Xin-She Yang,et al.  Introduction to Algorithms , 2021, Nature-Inspired Optimization Algorithms.

[24]  J. A. Barker,et al.  Monte Carlo studies of the dielectric properties of water-like models , 1973 .

[25]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[26]  You-Liang Zhu,et al.  A kinetic chain growth algorithm in coarse‐grained simulations , 2016, J. Comput. Chem..

[27]  Siewert J Marrink,et al.  Going Backward: A Flexible Geometric Approach to Reverse Transformation from Coarse Grained to Atomistic Models. , 2014, Journal of chemical theory and computation.

[28]  Denis J. Evans,et al.  The Nose–Hoover thermostat , 1985 .

[29]  Pierre de Buyl,et al.  H5MD: A structured, efficient, and portable file format for molecular data , 2013, Comput. Phys. Commun..

[30]  A. Khokhlov,et al.  Influence of cross-linking rate on the structure of hypercrosslinked networks: Multiscale computer simulation , 2016 .

[31]  Matthew Ming Fai Yuen,et al.  Material properties of the cross-linked epoxy resin compound predicted by molecular dynamics simulation , 2007 .

[32]  Giovanni Samaey,et al.  Coarse‐grained molecular dynamics simulations of polymerization with forward and backward reactions , 2018, J. Comput. Chem..

[33]  P. B. Warren,et al.  DISSIPATIVE PARTICLE DYNAMICS : BRIDGING THE GAP BETWEEN ATOMISTIC AND MESOSCOPIC SIMULATION , 1997 .

[34]  Kurt Kremer,et al.  Multiscale simulation of soft matter systems – from the atomistic to the coarse-grained level and back , 2009 .

[35]  Ilpo Vattulainen,et al.  Novel Methods in Soft Matter Simulations , 2013 .

[36]  Andrzej J. Rzepiela,et al.  Reconstruction of atomistic details from coarse‐grained structures , 2010, J. Comput. Chem..

[37]  F. Müller-Plathe,et al.  A Simple Reverse Mapping Procedure for Coarse-Grained Polymer Models with Rigid Side Groups , 2011 .

[38]  D. Roose,et al.  A molecular dynamics study of water transport inside an epoxy polymer matrix , 2015 .

[39]  H. Komber,et al.  Conversion dependence of the structural units and the degree of branching of a hyperbranched polyester based on 4,4-bis-(4'-hydroxyphenyl)pentanoic acid determined by NMR spectroscopy , 1999 .

[40]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

[41]  E. Peters,et al.  Multi-scale simulations for predicting material properties of a cross-linked polymer , 2015 .

[42]  Jakub Krajniak lab-tools: Set of Python scripts for analysing molecular dynamics simulations of polymers: release v1.3 , 2016 .

[43]  F. E. Karasz,et al.  Nature of the water-epoxy interaction , 1985 .

[44]  F. Müller-Plathe,et al.  Reactive molecular dynamics with material-specific coarse-grained potentials: growth of polystyrene chains from styrene monomers. , 2010, The journal of physical chemistry. B.

[45]  G. Samaey,et al.  Generic Adaptive Resolution Method for Reverse Mapping of Polymers from Coarse-Grained to Atomistic Descriptions. , 2016, Journal of chemical theory and computation.

[46]  William L. Jorgensen,et al.  Optimized intermolecular potential functions for liquid hydrocarbons , 1984 .

[47]  Bernd Ensing,et al.  Coarse Grained-to-Atomistic Mapping Algorithm: A Tool for Multiscale Simulations , 2007 .

[48]  Berk Hess,et al.  GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers , 2015 .

[49]  Thomas Brandes,et al.  ESPResSo++: A modern multiscale simulation package for soft matter systems , 2013, Comput. Phys. Commun..

[50]  E. Nies,et al.  A parallel algorithm for step- and chain-growth polymerization in molecular dynamics. , 2014, The Journal of chemical physics.

[51]  S. Neyertz,et al.  Preparation of bulk melt chain configurations of polycyclic polymers , 2001 .

[52]  P. Khalatur,et al.  Fully atomistic molecular dynamics simulation of nanosilica-filled crosslinked polybutadiene , 2016 .

[53]  A. Khokhlov,et al.  Hypercrosslinked polystyrene networks: an atomistic molecular dynamics simulation combined with a mapping/reverse mapping procedure. , 2014, The Journal of chemical physics.