Development of non-bonded interaction parameters between graphene and water using particle swarm optimization.

New Lennard-Jones parameters have been developed to describe the interactions between atomistic model of graphene, represented by REBO potential, and five commonly used all-atom water models, namely SPC, SPC/E, SPC/Fw, SPC/Fd, and TIP3P/Fs by employing particle swarm optimization (PSO) method. These new parameters were optimized to reproduce the macroscopic contact angle of water on a graphene sheet. The calculated line tension was in the order of 10-11 J/m for the droplets of all water models. Our molecular dynamics simulations indicate the preferential orientation of water molecules near graphene-water interface with one OH bond pointing toward the graphene surface. Detailed analysis of simulation trajectories reveals the presence of water molecules with ≤∼1, ∼2, and ∼4 hydrogen bonds at the surface of air-water interface, graphene-water interface, and bulk region of the water droplet, respectively. Presence of water molecules with ≤∼1 and ∼2 hydrogen bonds suggest the existence of water clusters of different sizes at these interfaces. The trends observed in the libration, bending, and stretching bands of the vibrational spectra are closely associated with these structural features of water. The inhomogeneity in hydrogen bond network of water at the air-water and graphene-water interface is manifested by broadening of the peaks in the libration band for water present at these interfaces. The stretching band for the molecules in water droplet shows a blue shift as compared to the pure bulk water, which conjecture the presence of weaker hydrogen bond network in a droplet.