Molecular Dynamics Simulation of Wetting Behavior: Contact Angle Dependency on Water Potential Models

The underlying principle of surface wettability has obtained great attentions for the development of novel functional surfaces. Molecular dynamics simulations has been widely utilized to obtain molecular-level details of surface wettability that is commonly quantified in term of contact angle of a liquid droplet on the surface. In this work, the sensitivity of contact angle calculation at various degrees of surface hydrophilicity to the adopted potential models of water: SPC/E, TIP4P, and TIP5P, is investigated. The simulation cell consists of a water droplet on a structureless surface whose hydrophilicity is modified by introducing a scaling factor to the water-surface interaction parameter. The simulation shows that the differences in contact angle described by the potential models are systematic and become more visible with the increase of the surface hydrophilicity. An alternative method to compute a contact angle based on the height of center-of-mass of the droplet is also evaluated, and the resulting contact angles are generally larger than those determined from the liquid-gas interfacial line.