Large effect of lateral box size in molecular dynamics simulations of liquid-solid friction

Abstract

Molecular dynamics simulations are a powerful tool to characterize liquid-solid friction. A slab configuration with periodic boundary conditions in the lateral dimensions is commonly used, where the measured friction coefficient could be affected by the finite lateral size of the simulation box. Here we show that for a very wetting liquid close to its melting temperature, strong finite size effects can persist up to large box sizes along the flow direction, typically ∼30 particle diameters. We relate the observed decrease of friction in small boxes to changes in the structure of the first adsorbed layer, which becomes less commensurable with the wall structure. Although these effects disappear for lower wetting cases or at higher temperatures, we suggest that the possible effect of the finite lateral box size on the friction coefficient should not be automatically set aside when exploring unknown systems.

Slip and Stress From Low Strain-Rate Nonequilibrium Molecular Dynamics: The Transient-Time Correlation Function Technique

“…A recent work 59 pointed out that small inhomogeneous systems, under certain conditions, tend to underestimate the solid-liquid friction. However, in this work we have not accounted for this effect, as we focussed solely on the derivation and application of the TTCF formalism under constant pressure conditions for the study of highly confined fluid systems.…”

Slip and stress from low shear rate nonequilibrium molecular dynamics: The transient-time correlation function technique

“…The high density of the wall is due to the reduced diameter of the wall particles and corresponds to ρ w = 1 if the quantities were normalized, assuming σ w = 1. A recent work 59 pointed out…”

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

“…Figure 1 (a) shows the MD simulation system of a quasi-1D Couette-type flow, where the basic setup is a standard one similar to our previous study. 31,32 The two parallel solid walls were fcc crystals and every pair of the nearest neighbors in the walls was bound through a harmonic potential Φ h (r) = k 2 (r − r eq ) 2 , with r being the interparticle distance, r eq = 0.277 nm, and k = 46.8 N/m. Interactions between fluid particles and between fluid and solid particles were modeled by a 12-6 LJ potential…”