Molecular dynamics (MD) simulations are an excellent tool for developing both fundamental microscopic understanding and quantitative predictions of many key aspects of heterogeneous nucleation. However, MD simulations of vapor bubble nucleation are challenging, due to very large nucleation barriers and small nucleation rates. Consequently, direct MD simulations of vapor nucleation are mainly limited to rapid heating, e.g., by a hot solid surface that brings the liquid close or to the spinodal line, where the nucleation barrier becomes zero. Such studies were performed for flat surfaces and around nanoparticles . While they can shed light on rapid laser heating induced vapor formation, the kinetics of this highly thermally driven and non-equilibrium process has little relevance to standard nucleation processes. In particular, it was observed that rapid heating of hydrophilic surfaces enables vapor formation, while hydrophobic surfaces are less effective in this respect . This behavior is completely opposite to what has been observed experimentally in "close to equilibrium" nucleation processes that are promoted by lower nucleation barriers present on weakly wetting surfaces . The likely reason for the behavior observed in aforementioned studies is that during rapid heating by hot solid hydrophilic surface, the rate of temperature increase of the liquid is faster due to lower interfacial thermal resistance. This effect is not present near equilibrium, when solid and liquids are at almost the same temperature regardless of the interfacial wetting properties.