Equilibrium molecular dynamics evaluation of the solid-liquid friction coefficient: role of timescales

Solid-liquid friction plays a key role in nanofluidic systems. Yet, despite decades of method development to quantify solid-liquid friction using molecular dynamics (MD) simulations, an accurate and widely applicable method is still missing. Here, we propose a method to quantify the solid-liquid friction coefficient (FC) from equilibrium MD simulations of a liquid confined between parallel solid walls. In this method, the FC is evaluated by fitting the Green-Kubo (GK) integral of the S-L shear force autocorrelation for the range of time scales where the GK integral slowly decays with time. The fitting function was derived based on the analytical solution considering the hydrodynamic equations in our previous work [H. Oga et al., Phys. Rev. Research 3, L032019 (2021)], assuming that the timescales related to the friction kernel and to the bulk viscous dissipation can be separated. By comparing the results with those of other equilibrium MD-based methods and those of non-equilibrium MD for a Lennard-Jones liquid between flat crystalline walls with different wettability, we show that the FC is extracted with excellent accuracy by the present method, even in wettability regimes where other methods become innacurate. We then show that the method is also applicable to grooved solid walls, for which the GK integral displays a complex behavior at short times. Overall, the present method extracts efficiently the FC for various systems, with easy implementation and low computational cost.Comment: 22 pages, 6 Figure

Equilibrium molecular dynamics evaluation of the solid-liquid friction coefficient: Role of timescales

Solid-liquid friction plays a key role in nanofluidic systems. Following the pioneering work of Bocquet and Barrat, who proposed to extract the friction coefficient (FC) from the plateau of the Green-Kubo (GK) integral of the solid-liquid shear force autocorrelation, the so-called plateau problem has been identified when applying the method to finite-sized molecular dynamics simulations, e.g., with a liquid confined between parallel solid walls. A variety of approaches have been developed to overcome this problem. Here, we propose another method that is easy to implement, makes no assumptions about the time dependence of the friction kernel, does not require the hydrodynamic system width as an input, and is applicable to a wide range of interfaces. In this method, the FC is evaluated by fitting the GK integral for the timescale range where it slowly decays with time. The fitting function was derived based on an analytical solution of the hydrodynamics equations [Oga et al., Phys. Rev. Res. 3, L032019 (2021)], assuming that the timescales related to the friction kernel and the bulk viscous dissipation can be separated. By comparing the results with those of other GK-based methods and non-equilibrium molecular dynamics, we show that the FC is extracted with excellent accuracy by the present method, even in wettability regimes where other GK-based methods suffer from the plateau problem. Finally, the method is also applicable to grooved solid walls, where the GK integral displays complex behavior at short times.Oga H., Omori T., Joly L., et al. Equilibrium molecular dynamics evaluation of the solid-liquid friction coefficient: Role of timescales, Journal of Chemical Physics, 159(2), 024701, 14 July 2023, © 2023 American Chemical Society. https://doi.org/10.1063/5.0155628

Shear force measurement of the hydrodynamic wall position in molecular dynamics

Flows in nanofluidic systems are strongly affected by liquid-solid slip, which is quantified by the slip length and by the position where the slip boundary condition applies. Here, we show that the viscosity, slip length, and hydrodynamic wall position (HWP) can be accurately determined from a single molecular dynamics (MD) simulation of a Poiseuille flow, after identifying a relation between the HWP and the wall shear stress in that configuration. From this relation, we deduce that in gravity-driven flows, the HWP identifies with the Gibbs dividing plane of the liquid-vacuum density profile. Simulations of a generic Lennard-Jones liquid confined between parallel frozen walls show that the HWP for a pressure-driven flow is also close to the Gibbs dividing plane (measured at equilibrium), which therefore provides an inexpensive estimate of the HWP, going beyond the common practice of assuming a given position for the hydrodynamic wall. For instance, we show that the HWP depends on the wettability of the surface, an effect usually neglected in MD studies of liquid-solid slip. Overall, the method introduced in this article is simple, fast, and accurate and could be applied to a large variety of systems of interest for nanofluidic applications.Cecilia Herrero, Takeshi Omori, Yasutaka Yamaguchi, and Laurent Joly, "Shear force measurement of the hydrodynamic wall position in molecular dynamics", The Journal of Chemical Physics 151, 041103 (2019) https://doi.org/10.1063/1.5111966

Full characterization of the hydrodynamic boundary condition at the atomic scale using an oscillating channel: Identification of the viscoelastic interfacial friction and the hydrodynamic boundary position

Omori T., Inoue N., Joly L., et al. Full characterization of the hydrodynamic boundary condition at the atomic scale using an oscillating channel: Identification of the viscoelastic interfacial friction and the hydrodynamic boundary position. Physical Review Fluids, 4(11), 114201. https://doi.org/10.1103/physrevfluids.4.114201