Understanding the asymmetry between advancing and receding microscopic contact angles

Measuring line tension: Thermodynamic integration during detachment of a molecular dynamics droplet

The contact line (CL) is where solid, liquid, and vapor phases meet, and Young's equation describes the macroscopic force balance of the interfacial tensions between these three phases. These interfacial tensions are related to the nanoscale stress inhomogeneity appearing around the interface, and for curved CLs, e.g., a three-dimensional droplet, another force known as the line tension must be included in Young's equation. The line tension has units of force, acting parallel to the CL, and is required to incorporate the extra stress inhomogeneity around the CL into the force balance. Considering this feature, Bey et al. [J. Chem. Phys. 152, 094707 (2020)] reported a mechanical approach to extract the value of line tension τℓ from molecular dynamics (MD) simulations. In this study, we show a novel thermodynamics interpretation of the line tension as the free energy per CL length, and based on this interpretation, through MD simulations of a quasi-static detachment process of a quasi-two-dimensional droplet from a solid surface, we obtained the value τℓ as a function of the contact angle. The simulation scheme is considered to be an extension of a thermodynamic integration method, previously used to calculate the solid-liquid and solid-vapor interfacial tensions through a detachment process, extended here to the three-phase system. The obtained value agreed well with the result by Bey et al. and showed the validity of thermodynamic integration at the three-phase interface.

The receding contact line cools down during dynamic wetting

When a contact line (CL)-where a liquid-vapor interface meets a substrate-is put into motion, it is well known that the contact angle differs between advancing and receding CLs. Using non-equilibrium molecular dynamics simulations, we reveal another intriguing distinction between advancing and receding CLs: while temperature increases at an advancing CL-as expected from viscous dissipation, we show that temperature can drop at a receding CL. Detailed quantitative analysis based on the macroscopic energy balance around the dynamic CL showed that the internal energy change of the fluid due to the change of the potential field along the pathline out of the solid-liquid interface induced a remarkable temperature drop around the receding CL, in a manner similar to latent heat upon phase changes. This result provides new insights for modeling the dynamic CL, and the framework for heat transport analysis introduced here can be applied to a wide range of nanofluidic systems.

Local stress tensor calculation by the method-of-plane in microscopic systems with macroscopic flow: A formulation based on the velocity distribution function

In this work, we developed a calculation method of local stress tensor applicable to non-equilibrium molecular dynamics (NEMD) systems, which evaluates the macroscopic momentum advection and the kinetic term of the stress in the framework of the Method-of-Plane (MoP), in a consistent way to guarantee the mass and momentum conservation. From the relation between the macroscopic velocity distribution function and the microscopic molecular passage across a fixed control plane, we derived a method to calculate the basic properties of the macroscopic momentum conservation law including the density, the velocity, the momentum flux, and the two terms of the stress tensor, i.e., the interaction and the kinetic terms, defined on a surface with a finite area. Any component of the streaming velocity can be obtained on a control surface, which enables the separation of the kinetic momentum flux into the advection and stress terms in the framework of MoP, and this enables strict satisfaction of the mass and momentum conservation for an arbitrary closed control volume (CV) set in NEMD systems. We validated the present method through the extraction of the density, velocity, and stress distributions in a quasi-one-dimensional steady-state Couette flow system and in a quasi-2D steady-state NEMD system with a moving contact line. We showed that with the present MoP, in contrast to the volume average method, the conservation law was satisfied even for a CV set around the moving contact line, which was located in a strongly inhomogeneous region.

Thermally Stimulated Wettability Transformations on One-Cycle Atomic Layer Deposition-Coated Cellulosic Paper: Applications for Droplet Manipulation and Heat Patterned Paper Fluidics

Cellulosic materials are widely used in daily life for paper products and clothing as well as for emerging applications in sustainable packaging and inexpensive medical diagnostics. Cellulose has a high density of hydroxyl groups that create strong intra- and interfiber hydrogen bonding. These abundant hydroxyl groups also make cellulose superhydrophilic. Schemes for hydrophobization and spatially selective hydrophobization of cellulosic materials can expand the application space for cellulose. Cellulose is often hydrophobized through wet chemistry surface modification methods. This work reports a new modification method using a combination of atomic layer deposition (ALD) and atmospheric heating to alter the wettability of purely cellulosic chromatography paper. We find that once the cellulosic paper is coated with a single ALD cycle (1cy-ALD) of Al₂O₃, it can be made sticky superhydrophobic after a 150 °C ambient post-ALD heating step. An X-ray photoelectron spectroscopy investigation reveals that the ALD-modified cellulosic surface becomes more susceptible to adsorption of adventitious carbon upon heating than an untreated cellulosic surface. This conclusion is further supported by the ability to use alternating air plasma and heat treatments to reversibly transition between the hydrophilic and hydrophobic states. We attribute the apparent abruptness of this wetting transition to a Cassie-Wenzel-like phenomenon, which is also consistent with the sticky hydrophobic wetting behavior. Using scanning probe methods, we show that the surfaces have roughness at multiple length scales. Using a Cassie-Wenzel model, we show how a small change in the surface's Young's contact angle-upon adsorption of adventitious carbon-can lead to an abrupt increase in hydrophobicity for surfaces with such roughnesses. Finally, we demonstrate the ability to spatially pattern the wettability on these 1cy-ALD-treated cellulosic papers via selective heating. This ALD-treated hydrophobic paper also shows promise for microliter droplet manipulation and patterned lab-on-paper devices.