Accurate simulation models for water interactions with graphene and graphite are important for nanofluidic applications, but existing force fields produce widely varying contact angles. Our extensive review of the experimental literature reveals extreme variation among reported values of graphene–water contact angles and a clustering of graphite–water contact angles into groups of freshly exfoliated (60° ± 13°) and not-freshly exfoliated graphite surfaces. The carbon–oxygen dispersion energy for a classical force field is optimized with respect to this 60° graphite–water contact angle in the infinite-force-cutoff limit, which in turn yields a contact angle for unsupported graphene of 80°, in agreement with the mean of the experimental results. Interaction force fields for finite cutoffs are also derived. A method for calculating contact angles from pressure tensors of planar equilibrium simulations that is ideally suited to graphite and graphene surfaces is introduced. Our methodology is widely applicable to any liquid-surface combination.
