Deformation of Surface Nanobubbles Induced by Substrate Hydrophobicity

Recent experimental measurements have shown that there exists a population of nanobubbles with different curvature radii, whereas both computer simulations and theoretical analysis indicated that the curvature radii of different nanobubbles should be the same at a given supersaturation. To resolve such inconsistency, we perform molecular dynamics simulations on surface nanobubbles that are stabilized by heterogeneous substrates either in the geometrical heterogeneity model (GHM) or in the chemical heterogeneity model (CHM) and propose that the inconsistency could be ascribed to the substrate-induced nanobubble deformation. We find that, as expected from theory and computer simulation, for either the GHM or the CHM, there exists a universal upper limit of contact angle for the nanobubbles, which is determined by the degree of supersaturation alone. By analyzing the evolution of the shape of nanobubbles as a function of substrate hydrophobicity that is controlled here by the liquid-solid interaction, two different origins of nanobubble deformation are identified. For substrates in the GHM, where the contact line is pinned by surface roughness, variation in the liquid-solid interaction changes only the location of the contact line and the measured contact angle, without causing a change in the nanobubble curvature. For substrates in the CHM, however, the liquid-solid interaction exerted by the bottom substrate can deform the vapor-liquid interface, resulting in variations in both the curvature of the vapor-liquid interface and the contact angle.