Molecular simulation of thin liquid films: Thermal fluctuations and instability

Threshold current density for diffusion-controlled stability of electrolytic surface nanobubbles

Understanding the stability mechanism of surface micro/nanobubbles adhered to gas-evolving electrodes is essential for improving the efficiency of water electrolysis, which is known to be hindered by the bubble coverage on electrodes. Using molecular simulations, the diffusion-controlled evolution of single electrolytic nanobubbles on wettability-patterned nanoelectrodes is investigated. These nanoelectrodes feature hydrophobic islands as preferential nucleation sites and allow the growth of nanobubbles in the pinning mode. In these simulations, a threshold current density distinguishing stable nanobubbles from unstable nanobubbles is found. When the current density remains below the threshold value, nucleated nanobubbles grow to their equilibrium states, maintaining their nanoscopic size. However, for the current density above the threshold value, nanobubbles undergo unlimited growth and can eventually detach due to buoyancy. Increasing the pinning length of nanobubbles increases the degree of nanobubble instability. By connecting the current density with the local gas oversaturation, an extension of the stability theory for surface nanobubbles [Lohse and Zhang, Phys. Rev. E 91, 031003(R) (2015)] accurately predicts the nanobubble behavior found in molecular simulations, including equilibrium contact angles and the threshold current density. For larger systems that are not accessible to molecular simulations, continuum numerical simulations with the finite difference method combined with the immersed boundary method are performed, again demonstrating good agreement between numerics and theories.

Non-equilibrium molecular simulations of thin film rupture

The retraction of thin films, as described by the Taylor-Culick (TC) theory, is subject to widespread debate, particularly for films at the nanoscale. We use non-equilibrium molecular dynamics simulations to explore the validity of the assumptions used in continuum models by tracking the evolution of holes in a film. By deriving a new mathematical form for the surface shape and considering a locally varying surface tension at the front of the retracting film, we reconcile the original theory with our simulation to recover a corrected TC speed valid at the nanoscale.

Thermal capillary waves on bounded nanoscale thin films

The effect of confining walls on the fluctuation of a nanoscale thin film's free surface is studied using stochastic thin-film equations (STFEs). Two canonical boundary conditions are employed to reveal the influence of the confinement: (1) an imposed contact angle and (2) a pinned contact line. A linear stability analysis provides the wave eigenmodes, after which thermal-capillary-wave theory predicts the wave fluctuation amplitudes. Molecular dynamics (MD) simulations are performed to test the predictions, and a Langevin diffusion model is proposed to capture oscillations of the contact lines observed in MD simulations. Good agreement between the theoretical predictions and the MD simulation results is recovered, and it is discovered that confinement can influence the entire film. Notably, a constraint on the length scale of wave modes is found to affect fluctuation amplitudes from our theoretical model, especially for 3D films. This opens up challenges and future lines of inquiry.

Viscous peeling of a nanosheet

Combining molecular dynamics (MD) and continuum simulations, we study the dynamics of propagation of a peeling front in a system composed of multilayered graphene nanosheets completely immersed in water. Peeling is induced by lifting one of the nanosheet edges with an assigned pulling velocity normal to the flat substrate. Using MD, we compute the pulling force as a function of the pulling velocity, and quantify the viscous resistance to the advancement of the peeling front. We compare the MD results to a 1D continuum model of a sheet loaded with modelled hydrodynamic loads. Our results show that the viscous dependence of the force on the velocity is negligible below a threshold velocity. Above this threshold, the hydrodynamics is mainly controlled by the viscous resistance associated to the flow near the crack opening, while lubrication forces are negligible owing to the large hydrodynamic slip at the liquid-solid boundary. Two dissipative mechanisms are identified: a drag resistance to the upward motion of the edge, and a resistance to the gap opening associated to the curvature of the flow streamlines near the entrance. Surprisingly, the shape of the sheet was found to be approximately independent of the pulling velocity even for the largest velocities considered.

Lattice Boltzmann simulations of stochastic thin film dewetting

We study numerically the effect of thermal fluctuations and of variable fluid-substrate interactions on the spontaneous dewetting of thin liquid films. To this aim, we use a recently developed lattice Boltzmann method for thin liquid film flows, equipped with a properly devised stochastic term. While it is known that thermal fluctuations yield shorter rupture times, we show that this is a general feature of hydrophilic substrates, irrespective of the contact angle θ. The ratio between deterministic and stochastic rupture times, though, decreases with θ. Finally, we discuss the case of fluctuating thin film dewetting on chemically patterned substrates and its dependence on the form of the wettability gradients.

Relaxation of Thermal Capillary Waves for Nanoscale Liquid Films on Anisotropic-Slip Substrates

The relaxation dynamics of thermal capillary waves for nanoscale liquid films on anisotropic-slip substrates are investigated using both molecular dynamics (MD) simulations and a Langevin model. The anisotropy of slip on substrates is achieved using a specific lattice plane of a face-centered cubic lattice. This surface's anisotropy breaks the simple scalar proportionality between slip velocity and wall shear stress and requires the introduction of a slip-coefficient tensor. The Langevin equation can describe both the growth of capillary wave spectra and the relaxation of capillary wave correlations, with the former providing a time scale for the surface to reach thermal equilibrium. Temporal correlations of interfacial Fourier modes, measured at thermal equilibrium in MD, demonstrate that (i) larger slip lengths lead to a faster decay in wave correlations and (ii) unlike isotropic-slip substrates, the time correlations of waves on anisotropic-slip substrates are wave-direction-dependent. These findings emerge naturally from the proposed Langevin equation, which becomes wave-direction-dependent, agrees well with MD results, and allows us to produce experimentally verifiable predictions.

Modification of the Fluctuation Dynamics of Ultrathin Wetting Films

We report on the effect of intermolecular forces on the fluctuations of supported liquid films. Using an optically induced thermal gradient, we form nanometer-thin films of wetting liquids on glass substrates, where van der Waals forces are balanced by thermocapillary forces. We show that the fluctuation dynamics of the film interface is strongly modified by intermolecular forces at lower frequencies. Data spanning three frequency decades are in excellent agreement with theoretical predictions accounting for van der Waals forces. Our results emphasize the relevance of intermolecular forces on thermal fluctuations when fluids are confined at the nanoscale.

Nanoscale thin-film flows with thermal fluctuations and slip

The combined effects of thermal fluctuations and liquid-solid slip on nanoscale thin-film flows are investigated using stochastic lubrication equations (SLEs). The previous no-slip SLE for films on plates is extended to consider slip effects and a new SLE for films on fibers is derived, using a long-wave approximation to fluctuating hydrodynamics. Analytically derived capillary spectra, which evolve in time, are found from the new SLEs and compared to molecular dynamics simulations. It is shown that thermal fluctuations lead to the generation and growth of surface waves, and slip accelerates this growth. SLEs developed here provide useful tools to study nanoscale film dewetting, nanofiber coating, and liquid transport using nanofibers.

Molecular physics of jumping nanodroplets

Next-generation processor-chip cooling devices and self-cleaning surfaces can be enhanced by a passive process that requires little to no electrical input, through coalescence-induced nanodroplet jumping. Here, we describe the crucial impact thermal capillary waves and ambient gas rarefaction have on enhancing/limiting the jumping speeds of nanodroplets on low adhesion surfaces. By using high-fidelity non-equilibrium molecular dynamics simulations in conjunction with well-resolved volume-of-fluid continuum calculations, we are able to quantify the different dissipation mechanisms that govern nanodroplet jumping at length scales that are currently difficult to access experimentally. We find that interfacial thermal capillary waves contribute to a large statistical spread of nanodroplet jumping speeds that range from 0-30 m s-1, where the typical jumping speeds of micro/millimeter sized droplets are only up to a few m s-1. As the gas surrounding these liquid droplets is no longer in thermodynamic equilibrium, we also show how the reduced external drag leads to increased jumping speeds. This work demonstrates that, in the viscous-dominated regime, the Ohnesorge number and viscosity ratio between the two phases alone are not sufficient, but that the thermal fluctuation number (Th) and the Knudsen number (Kn) are both needed to recover the relevant molecular physics at nanoscales. Our results and analysis suggest that these dimensionless parameters would be relevant for many other free-surface flow processes and applications that operate at the nanoscale.