Hydrogen-Bond Heterogeneity Boosts Hydrophobicity of Solid Interfaces

Icing and Adhesion Behaviors on Surfaces with Varied Lattice Constants

Investigating droplet wetting and icing behavior is crucial for comprehending the principles of surface icing and the design of anti-icing surfaces. In this study, we present the evidence from molecular dynamics (MD) simulations that reveal a hitherto unreported behavior of droplet wetting and icing adhesion on surfaces with lattice constants from 2.7 to 4.5 Å. Here, we observe that the contact angles (CA) of droplets on a face-centered cubic (FCC) lattice surface consistently correlate positively with the lattice constant. Further examination of droplet behavior on an idealized crystal surface reveals that hydrophilic surfaces (e.g., CA = 85°) inhibit freezing more effectively than hydrophobic surfaces (e.g., CA = 97°). This finding contradicts the conventional explanation that hydrophobic surfaces reduce heterogeneous nucleation, thereby delaying icing. This study introduces a mechanistic explanation for the promotion of water icing by hydrophobic surfaces and offers a novel design concept for the development of anti-ice surfaces in future applications.

Surprising Hydrophobic Polymer Surface with a High Content of Hydrophilic Polar Groups

The wetting property of a solid surface has been a hotspot for centuries, and many studies suggest that the hydrophobicity is highly related to the polar components. However, the underlying mechanism of polar moieties on the hydrophobicity remains unclear. Here, we tailor the surface polar moieties of epoxy resin (EP) by ozone modification and assess their wetting properties. Our results show that, for the modified EP with more (60.54%) polar moieties, the polar effect on hydrophobicity cannot be empirically observed. To reveal the underlying mechanism, the absorption parameters, including equilibrium distance, adsorption radius, and effective adsorption sites for water on EP before and after ozone treatment, are calculated on the basis of molecular simulations. After ozone modification, the equilibrium distance (from 1.95 to 1.70 Å), adsorption radius (from 3.80 to 4.50 Å), and effective adsorption sites (from 1 to 2) change slightly and the EP surface remains hydrophobic, although the polar groups significantly increase. Therefore, it is concluded that the wetting properties of solid surfaces are dominated by the equilibrium distance, adsorption radius, and effective adsorption sites for water on solids, and the nonlinear relationship between polar groups and hydrophilicity is clarified.

Partition Constant of Binary Mixtures for the Equilibrium between a Bulk and a Confined Phase

It is well-known that the thermodynamic, kinetic and structural properties of fluids, and in particular of water and its solutions, can be drastically affected in nanospaces. A possible consequence of nanoscale confinement of a solution is the partial segregation of its components. Thereby, confinement in nanoporous materials (NPM) has been proposed as a means for the separation of mixtures. In fact, separation science can take great advantage of NPM due to the tunability of their properties as a function of nanostructure, morphology, pore size, and surface chemistry. Alcohol-water mixtures are in this context among the most relevant systems. However, a quantitative thermodynamic description allowing for the prediction of the segregation capabilities as a function of the material-solution characteristics is missing. In the present study we attempt to fill this vacancy, by contributing a thermodynamic treatment for the calculation of the partition coefficient in confinement. Combining the multilayer adsorption model for binary mixtures with the Young equation, we conclude that the liquid-vapor surface tension and the contact angle of the pure substances can be used to predict the separation ability of a particular material for a given mixture to a semiquantitative extent. Moreover, we develop a Kelvin-type equation that relates the partition coefficient to the radius of the pore, the contact angle, and the liquid-vapor surface tensions of the constituents. To assess the validity of our thermodynamic formulation, coarse grained molecular dynamics simulations were performed on models of alcohol-water mixtures confined in cylindrical pores. To this end, a coarse-grained amphiphilic molecule was parametrized to be used in conjunction with the mW potential for water. This amphiphilic model reproduces some of the properties of methanol such as enthalpy of vaporization and liquid-vapor surface tension, and the minimum of the excess enthalpy for the aqueous solution. The partition coefficient turns out to be highly dependent on the molar fraction, on the interaction between the components and the confining matrix, and on the radius of the pore. A remarkable agreement between the theory and the simulations is found for pores of radius larger than 15 Å.

Impact of surface nanostructure and wettability on interfacial ice physics

Ice accumulation on solid surfaces is a severe problem for safety and functioning of a large variety of engineering systems, and its control is an enormous challenge that influences the safety and reliability of many technological applications. The use of molecular dynamics (MD) simulations is popular, but as ice nucleation is a rare event when compared to simulation timescales, the simulations need to be accelerated to force ice to form on a surface, which affects the accuracy and/or applicability of the results obtained. Here, we present an alternative seeded MD simulation approach, which reduces the computational cost while still ensuring accurate simulations of ice growth on surfaces. In addition, this approach enables, for the first time, brute-force all-atom water simulations of ice growth on surfaces unfavorable for nucleation within MD timescales. Using this approach, we investigate the effect of surface wettability and structure on ice growth in the crucial surface-ice interfacial region. Our main findings are that the surface structure can induce a flat or buckled overlayer to form within the liquid, and this transition is mediated by surface wettability. The first overlayer and the bulk ice compete to structure the intermediate water layers between them, the relative influence of which is traced using density heat maps and diffusivity measurements. This work provides new understanding on the role of the surface properties on the structure and dynamics of ice growth, and we also present a useful framework for future research on surface icing simulations.

Dynamic behaviors of interfacial water on the self-assembly monolayer (SAM) heterogeneous surface

Dynamic behaviors of water molecules near the surface with mixed hydrophobic and hydrophilic areas are studied by molecular dynamics simulation. More specifically, the diffusion coefficient and hydrogen bond lifetime of interfacial water on the self-assembly monolayer composed of hydrophobic and hydrophilic groups and their dependence on the mixing ratio are studied. The diffusion dramatically slows down, and the hydrogen bond lifetime considerably increases when a few hydrophilic groups are added to the hydrophobic surface. When the percentage of hydrophilic groups increases to 25%, the behavior of interfacial water is similar to the case of the pure hydrophilic surface. The sensitivity to the hydrophilic group can be attributed to the fact that the grafted hydrophilic groups can not only retard the directly bound water molecules but also affect indirectly bound water by stabilizing hydrogen bonds among interfacial water molecules.

Dynamical insights into the mechanism of a droplet detachment from a fiber

Quantifying the detachment behavior of a droplet from a fiber is important in many applications such as fog harvesting, oil-water separation, or water management in fuel cells. When the droplets are forcibly removed from hydrophilic fibers, the ease of detachment strongly depends on droplet volume and the rate of the process controlled by the applied force. Experiments, conducted on a ferrofluid under magnetic force, as well as continuum level calculations from fluid mechanics have so far been unable to resolve the time-dependent dynamics of droplet detachment and, most importantly, to assess the role of the applied force as the key determinant of the volume of the droplet residue remaining on the fiber after detachment. In the present work, we study the mechanism of water droplet detachment and retention of residual water on smooth hydrophilic fibers using nonequilibrium molecular dynamics simulations. We investigate how the applied force affects the breakup of a droplet and how the minimal detaching force per unit mass decreases with droplet size. We extract scaling relations that allow extrapolation of our findings to larger length scales that are not directly accessible by molecular models. We find that the volume of the residue on a fiber varies nonmonotonically with the detaching force, reaching the maximal size at an intermediate force and associated detachment time. The strength of this force decreases with the size of the drop, while the maximal residue increases with the droplet volume, V, sub-linearly, in proportion to the V2/3.

What Controls the Limit of Supercooling and Superheating of Pinned Ice Surfaces?

Cold-adapted organisms produce antifreeze proteins and glycoproteins to control the growth, melting and recrystallization of ice. It has been proposed that these molecules pin the crystal surface, creating a curvature that arrests the growth and melting of the crystal. Here we use thermodynamic modeling and molecular simulations to demonstrate that the curvature of the superheated or supercooled surface depends on the temperature and distances between ice-binding molecules, but not the details of their interactions with ice. We perform simulations of ice pinned with the antifreeze protein TmAFP, polyvinyl alcohol with different degrees of polymerization, and model ice-binding molecules to determine the thermal hystereses on melting and freezing, i.e. the maximum curvature that can be attained before, respectively, ice melts or grows irreversibly over the ice-binding molecules. We find that the thermal hysteresis is controlled by the bulkiness of the ice-binding molecules and their footprint at the ice surface. We elucidate the origin of the asymmetry between freezing and melting hysteresis found in experiments and propose guidelines to design synthetic antifreeze molecules with potent thermal hysteresis activity.

Hydrophobicity of proteins and nanostructured solutes is governed by topographical and chemical context

Hydrophobic interactions drive many important biomolecular self-assembly phenomena. However, characterizing hydrophobicity at the nanoscale has remained a challenge due to its nontrivial dependence on the chemistry and topography of biomolecular surfaces. Here we use molecular simulations coupled with enhanced sampling methods to systematically displace water molecules from the hydration shells of nanostructured solutes and calculate the free energetics of interfacial water density fluctuations, which quantify the extent of solute-water adhesion, and therefore solute hydrophobicity. In particular, we characterize the hydrophobicity of curved graphene sheets, self-assembled monolayers (SAMs) with chemical patterns, and mutants of the protein hydrophobin-II. We find that water density fluctuations are enhanced near concave nonpolar surfaces compared with those near flat or convex ones, suggesting that concave surfaces are more hydrophobic. We also find that patterned SAMs and protein mutants, having the same number of nonpolar and polar sites but different geometrical arrangements, can display significantly different strengths of adhesion with water. Specifically, hydroxyl groups reduce the hydrophobicity of methyl-terminated SAMs most effectively not when they are clustered together but when they are separated by one methyl group. Hydrophobin-II mutants show that a charged amino acid reduces the hydrophobicity of a large nonpolar patch when placed at its center, rather than at its edge. Our results highlight the power of water density fluctuations-based measures to characterize the hydrophobicity of nanoscale surfaces and caution against the use of additive approximations, such as the commonly used surface area models or hydropathy scales for characterizing biomolecular hydrophobicity and the associated driving forces of assembly.

Influence of Order within Nonpolar Monolayers on Hydrophobic Interactions

We report an experimental investigation of the influence of molecular-level order (crystallinity) within nonpolar monolayers on hydrophobic interactions. The measurements were performed using gold film-supported monolayers formed from alkanethiols (CH₃(CH₂)_nSH, with n ranging from 3 to 17), which we confirmed by using polarization-modulation infrared reflection-adsorption spectroscopy to exhibit chain-length-dependent order (methylene peak moves from 2926 to 2919 cm^-1, corresponding to a transition from liquid- to quasi-crystalline-like order) in the absence of substantial changes in chain density (constant methyl peak intensity). By using monolayer-covered surfaces immersed in either aqueous triethanolamine (TEA, 10 mM, pH 7.0) or pure methanol, we quantified hydrophobic and van der Waals contributions to adhesive interactions between identical pairs of surfaces (measured using an atomic force microscope) as a function of the length and order of the aliphatic chains within the monolayers. In particular, we measured pull-off forces arising from hydrophobic adhesion to increase in a nonlinear manner with chain length (abrupt increase between n = 5 and 6 from 2.1 ± 0.3 to 14.1 ± 0.7 nN) and to correlate closely with a transition from a liquid-like to crystalline-like monolayer phase. In contrast, adhesion in methanol increased gradually with chain length from 0.3 ± 0.1 to 3.2 ± 0.3 nN for n = 3 to 7 and then did not change further with an increase in chain length. These results lead to the hypothesis that order within nonpolar monolayers influences hydrophobic interactions. Additional support for this hypothesis was obtained from measurements reported in this paper using long-chain alkanethiols (ordered) and alkenethiols (disordered). The results are placed into the context of recent spectroscopic studies of hydrogen bonding of water at nonpolar monolayers. Overall, our study provides new insight into factors that influence hydrophobic interactions at nonpolar monolayers.

Parameterization of a coarse-grained model with short-ranged interactions for modeling fuel cell membranes with controlled water uptake

The coarse-grained model FF_pvap reproduces the experimental activity coefficient of water in tetramethylammonium chloride solutions over a wide range of concentrations, with a hundred-fold gain in computing efficiency with respect to atomistic models.

The Many Faces of Heterogeneous Ice Nucleation: Interplay Between Surface Morphology and Hydrophobicity

What makes a material a good ice nucleating agent? Despite the importance of heterogeneous ice nucleation to a variety of fields, from cloud science to microbiology, major gaps in our understanding of this ubiquitous process still prevent us from answering this question. In this work, we have examined the ability of generic crystalline substrates to promote ice nucleation as a function of the hydrophobicity and the morphology of the surface. Nucleation rates have been obtained by brute-force molecular dynamics simulations of coarse-grained water on top of different surfaces of a model fcc crystal, varying the water-surface interaction and the surface lattice parameter. It turns out that the lattice mismatch of the surface with respect to ice, customarily regarded as the most important requirement for a good ice nucleating agent, is at most desirable but not a requirement. On the other hand, the balance between the morphology of the surface and its hydrophobicity can significantly alter the ice nucleation rate and can also lead to the formation of up to three different faces of ice on the same substrate. We have pinpointed three circumstances where heterogeneous ice nucleation can be promoted by the crystalline surface: (i) the formation of a water overlayer that acts as an in-plane template; (ii) the emergence of a contact layer buckled in an ice-like manner; and (iii) nucleation on compact surfaces with very high interaction strength. We hope that this extensive systematic study will foster future experimental work aimed at testing the physiochemical understanding presented herein.