We study by molecular dynamics simulations the driving force for the hydrophobic interaction between graphene sheets of different sizes down to the atomic scale. Similar to the prediction by Lum, Chandler, and Weeks for hard-sphere solvation [J. Phys. Chem. B 1999, 103, 4570-4577], we find the driving force to be length-scale dependent, despite the fact that our model systems do not exhibit dewetting. For small hydrophobic solutes, the association is purely entropic, while enthalpy favors dissociation. The latter is demonstrated to arise from the enhancement of hydrogen bonding between the water molecules around small hydrophobes. On the other hand, the attraction between large graphene sheets is dominated by enthalpy which mainly originates from direct solute-solute interactions. The crossover length is found to be inside the range of 0.3-1.5 nm(2) of the surface area of the hydrophobe that is eliminated in the association process. In the large-scale regime, different thermodynamic properties are scalable with this change of surface area. In particular, upon dimerization, a total and a water-induced stabilization of approximately 65 and 12 kJ/mol/nm(2) are obtained, respectively, and on average around one hydrogen bond is gained per 1 nm(2) of graphene sheet association. Furthermore, the potential of mean force between the sheets is also scalable except for interplate distances smaller than 0.64 nm which corresponds to the region around the barrier for removing the last layer of water. It turns out that, as the surface area increases, the relative height of the barrier for association decreases and the range of attraction increases. It is also shown that, around small hydrophobic solutes, the lifetime of the hydrogen bonds is longer than in the bulk, while around large hydrophobes it is the same. Nevertheless, the rearrangement of the hydrogen-bond network for both length-scale regimes is slower than in bulk water.
Read full abstract