Water wettability of graphene and graphite, optimization of solid-liquid interaction force fields, and insights from mean-field modeling

Body Forces Drive the Apparent Line Tension of Sessile Droplets

The line tension of a three-phase contact line is implicated in a wide variety of interfacial phenomena, but there is ongoing controversy, with existing measurements spanning six orders of magnitude in both signs. Here, we show that computationally obtained magnitudes, sign changes, and nontrivial variations of apparent line tension can be faithfully reproduced in a parsimonious model that incorporates only liquid-substrate interactions. Our results suggest that the origin for the remarkable variation lies in the failure of a widely used estimation method to eliminate body forces, leading measured line tensions to behave like an extensive quantity.

Comparative first principles-based molecular dynamics study of catalytic mechanism and reaction energetics of water oxidation reaction on 2D-surface

The study of the water-splitting process, which can proceed in 2e^- as well as 4e^- pathway, reveals that the process is entirely an uphill process, and the third step, that is, the oxooxo bond formation is the rate-determining step. The kinetic barrier of the oxygen evolution reaction (OER) on the 2D material catalysts in the presence of explicit solvents is scarcely studied. Here, we investigate the dynamics of the OER on the undoped graphene and the activation energy barrier of each step using first principles molecular dynamics simulations. Here we provide a detailed analysis of the kinetics of all the 4e^- transfer steps of OER on the graphene surface. We also compare the accuracy of one of the density functional theory (DFT) functionals and density functional based tight binding (DFTB) method in explaining the OER steps. The comparative study reveals that DFTB can be used for performing metadynamics simulations quipped with much less computational cost than DFT functionals. By both Perdew-Burke-Ernzerhof and DFTB methods, the third step is revealed to be the rate-determining step with an energy barrier of 21.19 ± 0.51 and 20.23 ± 0.20 kcal mol^-1 , respectively. DFTB gives an impression of being successful in predicting the energy barriers of OER in 4e- transfer pathway and comparable to the DFT method, and we would like to extend the use of DFTB for further studies with a sizable and complex system.

Hydrocarbon contamination in angström-scale channels

Nonspecific molecular adsorption such as airborne contamination occurs on most surfaces including those of 2D materials and alters their properties. While surface contamination is studied using a plethora of techniques, the effect of contamination on confined systems such as nanochannels/pores leading to their clogging is still lacking. We report a systematic investigation of hydrocarbon adsorption in angstrom (Å) slit channels of varying heights. Hexane is chosen to mimic the hydrocarbon contamination and the clogging of the Å-channels is evaluated via a helium gas flow measurement. The level of hexane adsorption, in other words, the degree of clogging depends on the size difference between the channels and hexane. A dynamic transition of the clogging and revival process is shown in sub-2 nm thin channels. Long-term storage and stability of our Å-channels are demonstrated here for up to three years, alleviating the contamination and unclogging the channels using thermal treatment. This study highlights the importance of the nanochannels' stability and demonstrates the self-cleansing nature of sub-2 nm thin channels enabling a robust platform for molecular transport and separation studies. We provide a method to assess the cleanliness of nanoporous membranes, which is vital for the practical applications of nanofluidics in various fields such as molecular sensing, separation and power generation.

Effects of the Interfacial Modeling Approach on Equilibrium Calculations of Slip Length for Nanoconfined Water in Carbon Slits

In this investigation, equilibrium molecular dynamics simulations were conducted to assess the influence of the interface modeling approach on the calculation of hydrodynamic slip in carbon nanochannels. A Green-Kubo formalism was implemented for the calculation of the slip length in water confined by graphite layers. The nonbonded interactions between solid and liquid atoms (interface models) were modeled using parameters optimized to represent the wetting behavior and adsorption energy curves from electronic structure calculations. Conventional carbon-oxygen-only interaction models were compared against comprehensive models able to represent the molecular-orientation-dependent energy of interaction. Quasi-universal relationships built under the premise of the slip length dependence on the water-graphite affinity and characterized by macroscopic wettability were critically assessed. It was found that the wetting behavior cannot fully characterize the hydrodynamic slip because interface models that produced the same surface wettability yielded different values of the friction coefficient. Alternatively, the density depletion length, used to characterize the interfacial liquid structuring and the availability of momentum carriers (interfacial water molecules), was able to accurately represent the slip length trends independently of the interface model. These findings reassert the importance of physically sound interface models to study interfacial transport properties and the need of reliable parameters and characterization procedures to support theoretical models that seek to unveil the inconsistencies in hydrodynamic slip calculations.

Dynamic wetting of solid-liquid-liquid system by molecular kinetic theory

HYPOTHESIS

The mechanisms of dynamic wetting of solid-liquid-liquid (SLL) system, especially the viscosity effects of two liquids, can be investigated by the molecular kinetic theory (MKT).

METHODS

The molecular kinetic theory combined with published data was used to study the roles of a fluid viscosity and a solid surface in dynamic wetting.

FINDINGS

First, the MKT on dynamic wetting was introduced and its limitation was analyzed. Second, a viscosity effect and a solid surface effect were considered. The viscosity effect was divided into three parts for the first time, including two pure liquid zones and a mixing zone. Third, a coefficient activation free energy model was proposed, considering the effects of mixing liquids and a solid surface. Finally, the key parameters in the MKT and the application and validation of the coefficient activation free energy model were discussed in detail. This model can explain the energy dissipation in a vicinity of a three-phase contact-line successfully in a SLL wetting system. This work sheds light on the physical mechanisms of fluid and solid surface properties on the dynamic wetting in a SLL system.