Partial pair distribution functions, XRD patterns, segmental dynamics, elastic constants and glass temperature in EPON862/DETDA epoxy predicted through molecular dynamics.

We examine the mechanical response of single layer graphene nanoribbons (GNR) under constant compressive loads through molecular dynamics simulations. Compressive stress-strain curves are presented for GNRs of various lengths and widths. The dependence of GNR’s buckling resistance on its size, aspect ratio, and chiral angle is discussed and approximate corresponding relations are provided. A single master curve describing the dependence of the critical buckling stress of GNRs on their aspect ratio is presented. Our findings were compared to the continuum elasticity theories for wide plates and wide columns. In the large width limit, the response of the GNRs agrees with the predictions of the wide plates theory and thus, with that of wide graphenes. In the small width limit, the behavior of graphene nanoribbons deviates from that of periodic graphenes due to various edge related effects which govern the stiffness and the stability of the graphene membranes, but it qualitatively agrees with the theory of wide columns. In order to assess the effect of thermal fluctuations on the critical buckling stress a wide range of temperatures is examined. The findings of the current study could provide important insights regarding the feasibility and the evaluation of the performance of graphene-based devices.

In previous work by the authors, a new methodology was developed for Brownian dynamics/kinetic Monte Carlo (BD/kMC) simulations of polymer melts. In this study, this methodology is extended for dynamical simulations of crosslinked polymer networks in a coarse-grained representation, wherein chains are modeled as sequences of beads, each bead encompassing a few Kuhn segments. In addition, the C++ code embodying these simulations, entitled Engine for Mesoscopic Simulations for Polymer Networks (EMSIPON) is described in detail. A crosslinked network of cis-1,4-polyisoprene is chosen as a test system. From the thermodynamic point of view, the system is fully described by a Helmholtz energy consisting of three explicit contributions: entropic springs, slip springs and non-bonded interactions. Entanglements between subchains in the network are represented by slip springs. The ends of the slip springs undergo thermally activated hops between adjacent beads along the chain backbones, which are tracked by kinetic Monte Carlo simulation. In addition, creation/destruction processes are included for the slip springs at dangling subchain ends. The Helmholtz energy of non-bonded interactions is derived from the Sanchez–Lacombe equation of state. The isothermal compressibility of the polymer network is predicted from equilibrium density fluctuations in very good agreement with the underlying equation of state and with experiment. Moreover, the methodology and the corresponding C++ code are applied to simulate elongational deformations of polymer rubbers. The shear stress relaxation modulus is predicted from equilibrium simulations of several microseconds of physical time in the undeformed state, as well as from stress-strain curves of the crosslinked polymer networks under deformation.

The mechanical response of single and multiple graphene sheets under uniaxial compressive loads was studied with molecular dynamics (MD) simulations, using different semi-empirical force fields at different boundary conditions or constrains. Compressive stress–strain curves were obtained and the critical stress/strain values were derived. The MD results are compared to the linear elasticity continuum theory for loaded slabs. Concerning the length dependence of critical values, qualitatively similar behavior is observed between the theory and numerical simulations for single layer graphenes, as the critical stress/strain for buckling was found to scale to the inverse squared length. However discrepancies were noted for multilayer graphenes, where the critical buckling stress also decreased with increasing length, though at a slower rate than expected from elastic buckling analysis.