Classical nucleation of vapor between hydrophobic plates

The surface tension of Martini 3 water mixtures

The Martini model, a coarse-grained forcefield for biomolecular simulations, has experienced a vast increase in popularity in the past decade. Its building-block approach balances computational efficiency with high chemical specificity, enabling the simulation of organic and inorganic molecules. The modeling of coarse-grained beads as Lennard-Jones particles poses challenges for the accurate reproduction of liquid-vapor interfacial properties, which are crucial in various applications, especially in the case of water. The latest version of the forcefield introduces refined interaction parameters for water beads, tackling the well-known artifact of Martini water freezing at room temperature. In addition, multiple sizes of water beads are available for simulating the solvation of small cavities, including the smallest pockets of proteins. This work focuses on studying the interfacial properties of Martini water, including surface tension and surface thickness. Employing the test-area method, we systematically compute the liquid-vapor surface tension across various combinations of water bead sizes and for temperatures from 300 to 350 K. These findings are of interest to the Martini community as they allow users to account for the low interfacial tension of Martini water by properly adjusting observables computed via coarse-grained simulations to allow for accurate matching against all-atom or experimental results. Surface tension data are also interpreted in terms of local enrichment of the various mixture components at the liquid-vapor interface by means of Gibbs' adsorption formalism. Finally, the critical scaling of the Martini surface tension with temperature is reported to be consistent with the critical exponent of the 3D Ising universality class.

Contact angle hysteresis on nonwetting microstructured surfaces: Effect of randomly distributed pillars or holes

We present a numerical study of the advancing and receding apparent contact angles for a liquid meniscus in contact with an ultrahydrophobic surface with randomly distributed microsized pillars or holes in the Cassie's wetting regime. We study the Wilhelmy plate system in the framework of the full capillary model to obtain these angles using the heterogeneous surface approximation model for a broad interval of values of pillar or hole concentration and for both square and circular shapes of the pillars or holes cross-section. Three types of random placing of defects on the plate are investigated, i.e., two with restrictions: (1) with maximum and (2) with minimum distance between the defects (in these cases the defects are isolated), and (3) without restrictions (the defects can overlap). The results show that the type of defect distribution and also the type of the defects shape (circular or square) does not affect the magnitude of the two angles. The results of the numerical simulations showed that the retention force for a plate with randomly located defects is not greater, and for larger concentrations of pillars or holes, it is smaller than that for periodically spaced ones. Comparisons with experimental results for the receding contact angle on surfaces with pillars and with the advancing contact angle on surfaces with periodically arranged holes is carried out.

Reversible Surface Energy Storage in Molecular-Scale Porous Materials

Forcible wetting of hydrophobic pores represents a viable method for energy storage in the form of interfacial energy. The energy used to fill the pores can be recovered as pressure-volume work upon decompression. For efficient recovery, the expulsion pressure should not be significantly lower than the pressure required for infiltration. Hysteresis of the wetting/drying cycle associated with the kinetic barrier to liquid expulsion results in energy dissipation and reduced storage efficiency. In the present work, we use open ensemble (Grand Canonical) Monte Carlo simulations to study the improvement of energy recovery with decreasing diameters of planar pores. Near-complete reversibility is achieved at pore widths barely accommodating a monolayer of the liquid, thus minimizing the area of the liquid/gas interface during the cavitation process. At the same time, these conditions lead to a steep increase in the infiltration pressure required to overcome steric wall/water repulsion in a tight confinement and a considerable reduction in the translational entropy of confined molecules. In principle, similar effects can be expected when increasing the size of the liquid particles without altering the absorbent porosity. While the latter approach is easier to follow in laboratory work, we discuss the advantages of reducing the pore diameter, which reduces the cycling hysteresis while simultaneously improving the stored-energy density in the material.

Hydrophobicity of molecular-scale textured surfaces: The case of zeolitic imidazolate frameworks, an atomistic perspective

Hydrophobicity has proven fundamental in an inexhaustible amount of everyday applications. Material hydrophobicity is determined by chemical composition and geometrical characteristics of its macroscopic surface. Surface roughness or texturing enhances intrinsic hydrophilic or hydrophobic characteristics of a material. Here we consider crystalline surfaces presenting molecular-scale texturing typical of crystalline porous materials, e.g., metal-organic frameworks. In particular, we investigate one such material with remarkable hydrophobic qualities, ZIF-8. We show that ZIF-8 hydrophobicity is driven not only by its chemical composition but also its sub-nanoscale surface corrugations, a physical enhancement rare amongst hydrophobes. Studying ZIF-8's hydrophobic properties is challenging as experimentally it is difficult to distinguish between the materials' and the macroscopic corrugations' contributions to the hydrophobicity. The computational contact angle determination is also difficult as the standard "geometric" technique of liquid nanodroplet deposition is prone to many artifacts. Here, we characterise ZIF-8 hydrophobicity via: (i) the "geometric" approach and (ii) the "energetic" method, utilising the Young-Dupré formula and computationally determining the liquid-solid adhesion energy. Both approaches reveal nanoscale Wenzel-like bathing of the corrugated surface. Moreover, we illustrate the importance of surface linker termination in ZIF-8 hydrophobicity, which reduces when varied from sp3 N to sp2 N termination. We also consider halogenated analogues of the methyl-imidazole linker, which promote the transition from nanoWenzel-like to nanoCassie-Baxter-like states, further enhancing surface hydrophobicity. Present results reveal the complex interface physics and chemistry between water and complex porous, molecular crystalline surfaces, providing a hint to tune their hydrophobicity.

What keeps nanopores boiling

The liquid-to-vapor transition can occur under unexpected conditions in nanopores, opening the door to fundamental questions and new technologies. The physics of boiling in confinement is progressively introduced, starting from classical nucleation theory, passing through nanoscale effects, and terminating with the material and external parameters that affect the boiling conditions. The relevance of boiling in specific nanoconfined systems is discussed, focusing on heterogeneous lyophobic systems, chromatographic columns, and ion channels. The current level of control of boiling in nanopores enabled by microporous materials such as metal organic frameworks and biological nanopores paves the way to thrilling theoretical challenges and to new technological opportunities in the fields of energy, neuromorphic computing, and sensing.