文摘
Likely candidates for the global potential energy minima of (H2O)n clusters with n 21 on the (0001) surfaceof graphite are found using basin-hopping global optimization. The potential energy surfaces are constructedusing the TIP4P intermolecular potentials for the water molecules (the TIP3P is also explored as a secondarychoice), a Lennard-Jones water-graphite potential, and a water-graphite polarization potential that is builtfrom classical electrostatic image methods and takes into account both the perpendicular and the parallelelectric polarizations of graphite. This potential energy surface produces a rather hydrophobic water-graphiteinteraction. As a consequence, the water component of the lowest graphite-(H2O)n minima is quite closelyrelated to low-lying minima of the corresponding TIP4P (H2O)n clusters. In about half of the cases, thegeometrical substructure of the water molecules in the graphite-(H2O)n global minimum coincides with thatof the corresponding free water cluster. Exceptions occur when the interaction with graphite induces a changein geometry. A comparison of our results with available theoretical and experimental data is performed.