Molecular dynamics simulations show that water nanodroplets on hierarchically structured copper surfaces (gecko-like pillars bearing nanoscale protuberances) can be more hydrophobic than flat-topped pillars, but the equilibrium contact angles deviate unexpectedly from the classical Wenzel roughness prediction, most strongly at short groove widths. The open problem is twofold and linked: (1) establish how the full hierarchy parameter space — primary groove width, pillar height, protuberance size and density, number of hierarchy levels, and surface characteristic energy, singly and jointly — maps onto nanodroplet contact angle, so that general design principles for nanostructured hydrophobic metal surfaces emerge; and (2) identify the physical mechanism responsible for the nanoscale failure of the Wenzel model at short groove widths and construct a modified, nanoscale-valid wetting model that reduces to the Wenzel prediction in the macroscopic limit. No prior work validates a modified Wenzel-type model against groove-width-dependent nanodroplet simulations or completes the systematic hierarchy-parameter scan.

Background

When a liquid droplet rests on an ideal flat solid, its equilibrium contact angle θ\theta is set by Young's equation, balancing solid-vapor, solid-liquid, and liquid-vapor interfacial tensions. Real surfaces are rough, and the classical correction is Wenzel's 1936 model, which predicts that roughness amplifies the intrinsic wetting behavior: the apparent contact angle on a rough surface obeys cos⁡θW=rcos⁡θY\cos\theta_W = r\cos\theta_Y, where rr is the ratio of true solid surface area to its projected area and θY\theta_Y is Young's angle. A complementary description, the Cassie state, applies when air remains trapped under the droplet so it partially sits on vapor. Many natural water-repellent surfaces (lotus leaves, gecko skin) are hierarchical — they carry roughness on two or more length scales, e.g. microscale pillars decorated with nanoscale protuberances — and such dual-scale textures are central to designing robust superhydrophobic and low-adhesion surfaces. At the nanoscale these continuum pictures become questionable, and the tool of choice is molecular dynamics (MD), in which the trajectories of individual water molecules interacting with explicitly modeled solid atoms are integrated numerically; the apparent contact angle is measured from the equilibrium droplet shape. In MD the strength of the solid-liquid attraction is a tunable parameter (often called the surface characteristic energy, calibrated so that a flat surface reproduces a chosen intrinsic contact angle). Zhang, Ma, Wang, Hao and Zheng (Colloids and Surfaces A 604, 2020, 125291) used MD to study water nanodroplets on copper surfaces textured with gecko-like hierarchical pillars (pillars topped by a single size of smaller protuberances), comparing them with nonhierarchical pillars. They scanned four primary groove widths and three surface characteristic energies and found that the hierarchical pillars tend to be more hydrophobic, but that increasing surface energy or groove width weakens this advantage; when the simulated contact angles were compared with theoretical (Wenzel-type) predictions, an unexpected deviation appeared at the nanoscale, especially for short groove widths. This echoes a broader literature on the limited applicability of the Wenzel and Cassie equations: Nosonovsky (2007) analyzed their range of validity; Yen and Soong (2016) proposed a hybrid Cassie-Wenzel model attributing the breakdown of macroscopic predictions on nanoscale roughness to the finite width of the wall-fluid depletion region (the near-wall layer depleted of fluid) and to coexistence of Cassie and Wenzel states; Kanduc (2017) showed that contact angles of water nanodroplets depend on droplet size through the Tolman correction (the curvature dependence of surface tension) and a contact-angle-dependent generalized line tension (the free-energy cost per unit length of the three-phase contact line); and Klauser, von Kleist-Retzow and Fatikow (2022) reviewed line-tension and drop-size effects at the nanoscale. None of these, however, explains or repairs the groove-width-dependent failure of the Wenzel prediction for nanodroplets on hierarchically structured metal surfaces, and the systematic exploration of the hierarchy parameter space (pillar height, protuberance size and density, multiple hierarchy levels, and their combined effects) called for by the source study has not been carried out.

Problem Statement

For water nanodroplets on copper surfaces textured with hierarchical (gecko-like) nanostructures, the source study poses a compound open problem whose single overarching objective is a quantitative, nanoscale-valid account of how surface geometry determines wettability, usable as a design principle for nanostructured hydrophobic metal surfaces. It comprises two linked parts: (1) Determine how the equilibrium contact angle of water nanodroplets depends on the hierarchical geometry — primary groove width, pillar height, size and density of the protuberances on the pillar tops, and the number of hierarchy levels — singly and jointly with the surface characteristic energy, extending beyond the four groove widths and three energies already simulated, so that general design principles for hierarchical hydrophobic surfaces can be stated. (2) Explain the observed deviation of the simulated contact angles from the Wenzel prediction at the nanoscale, most pronounced at short groove widths: identify the responsible mechanism (candidates from adjacent literature include line tension and its generalizations, Tolman-type curvature corrections to surface tension, the finite width of the wall-fluid depletion region, partial penetration or coexistence of Cassie and Wenzel states, and finite-range solid-fluid interactions), and construct a modified Wenzel-type model, reducing to the classical prediction in the macroscopic limit, that predicts nanodroplet contact angles on such structured surfaces to within simulation accuracy. Part (1) supplies the empirical geometry-to-contact-angle map; part (2) supplies the predictive model that makes that map trustworthy at the nanoscale, and neither alone answers the source's question of how to design hierarchical nanostructures for controlled wettability.

References

  1. Zhang M., Ma L., Wang Q., Hao P., Zheng X., 'Wettability behavior of nanodroplets on copper surfaces with hierarchical nanostructures', Colloids and Surfaces A: Physicochemical and Engineering Aspects 604:125291 (2020). DOI 10.1016/j.colsurfa.2020.125291, https://doi.org/10.1016/j.colsurfa.2020.125291
  2. Wenzel R. N., 'Resistance of solid surfaces to wetting by water', Industrial & Engineering Chemistry 28:988-994 (1936). DOI 10.1021/ie50320a024, https://doi.org/10.1021/ie50320a024
  3. Nosonovsky M., 'On the Range of Applicability of the Wenzel and Cassie Equations', Langmuir 23:9919-9920 (2007). DOI 10.1021/la701324m, https://doi.org/10.1021/la701324m
  4. Yen T.-H., Soong C.-Y., 'Hybrid Cassie-Wenzel model for droplets on surfaces with nanoscale roughness', Physical Review E 93:022805 (2016). DOI 10.1103/physreve.93.022805, https://doi.org/10.1103/PhysRevE.93.022805
  5. Kanduč M., 'Going beyond the standard line tension: Size-dependent contact angles of water nanodroplets', The Journal of Chemical Physics 147:174701 (2017). DOI 10.1063/1.4990741, https://doi.org/10.1063/1.4990741
  6. Klauser W., von Kleist-Retzow F. T., Fatikow S., 'Line Tension and Drop Size Dependence of Contact Angle at the Nanoscale', Nanomaterials 12:369 (2022). DOI 10.3390/nano12030369, https://doi.org/10.3390/nano12030369
  7. Qiang W., Lan Z., Liu Y., Liu Y., Xu W., Wen R., Ma X., 'Excellent hydrophobicity and high mobility of condensate droplets on hierarchical nanostructured surfaces: Insights from MD simulations', Journal of Molecular Liquids 407:125139 (2024). DOI 10.1016/j.molliq.2024.125139, https://doi.org/10.1016/j.molliq.2024.125139
  8. He X., Fan J., Cui K., Tian Y., Wang T., Wang S., 'Wetting of nanoscale water films on hierarchically structured surfaces', Physics of Fluids 36:052016 (2024). DOI 10.1063/5.0205294, https://doi.org/10.1063/5.0205294
  9. Wang Y., You M., Sun Q., 'Effects of Two-Level Surface Roughness on Superhydrophobicity', Coatings 15:1269 (2025). DOI 10.3390/coatings15111269, https://doi.org/10.3390/coatings15111269
  10. Wang X., Fu C., Zhang C., Qiu Z., Wang B., 'A Comprehensive Review of Wetting Transition Mechanism on the Surfaces of Microstructures from Theory and Testing Methods', Materials 15:4747 (2022). DOI 10.3390/ma15144747, https://doi.org/10.3390/ma15144747