Coarse grained models of graphene and graphene oxide for use in aqueous solution

Obtaining stable aqueous dispersions of graphene-based materials is a major obstacle in the development and widespread use of graphene in nanotechnology. The efficacy of atomistic simulations in obtaining a molecular-level insight into aggregation and exfoliation of graphene/graphene oxide (GO) is hindered by length and time scale limitations. In this work, we developed coarse-grained (CG) models of graphene/GO sheets, compatible with the polarizable Martini water model, using molecular dynamics, iterative Boltzmann inversion and umbrella sampling simulations. The new CG models accurately reproduce graphene/GO–water radial distribution functions and sheet–sheet aggregation free energies for small graphene (−316 kJ mol −1 ) and GO (−108 kJ mol −1 ) reference sheets. Deprotonation of carboxylic acid functionalities stabilize the exfoliated state by electrostatic repulsion, providing they are present at sufficiently high surface concentration. The simulations also highlight the pivotal role played by entropy in controlling the propensity for aggregation or exfoliation. The CG models improve the computational efficiency of simulations by an order of magnitude and the framework presented is transferrable to sheets of different sizes and oxygen contents. They can now be used to provide fundamental physical insights into the stability of dispersions and controlled self-assembly, underpinning the computational design of graphene-containing nanomaterials.