Molecular dynamics simulation of high slip flow of water confined between graphene nanochannels at experimentally accessible shear rates

C Carmelo Riccardo Civello (Department of Mathematics, Faculty of Science, Computing and Emerging Technologies, Swinburne University of Technology 1 , P.O. Box 218, Hawthorn, Victoria 3122,) L Luca Maffioli (Department of Mathematics, Faculty of Science, Computing and Emerging Technologies, Swinburne University of Technology 1 , P.O. Box 218, Hawthorn, Victoria 3122,) E Edward R. Smith (Mechanical and Aerospace Engineering, Brunel University London 2 , Kingston Lane, Uxbridge, UB8 3PH London,) J James P. Ewen (Department of Mechanical Engineering, Imperial College London 3 , South Kensington Campus, London, SW7 2AZ London,) P Peter J. Daivis (Department of Physics, RMIT University 5 , GPO Box 2476, Melbourne 3001, Victoria,) D Daniele Dini (Department of Mechanical Engineering, Imperial College London 3 , South Kensington Campus, London, SW7 2AZ London,) B B. D. Todd (Department of Mathematics, Faculty of Science, Computing and Emerging Technologies, Swinburne University of Technology 1 , P.O. Box 218, Hawthorn, Victoria 3122,)

Abstract

The transient time correlation function (TTCF) method has emerged as a powerful methodology for accurately probing systems at low shear rates. In the present study, TTCF was used to evaluate the shear rate dependence of the slip length in a high-slip system consisting of water confined between graphene walls at experimentally accessible shear rates, for which classical nonequilibrium molecular dynamics (NEMD) is unfeasible. The corresponding Navier friction coefficient was computed for all shear rates spanning six orders of magnitude and compared with the equilibrium limit. We report for the first time NEMD results obtained at experimentally accessible shear rates using the TTCF approach for a system that has attracted significant interest over the past decades. The slip length calculated with TTCF is in good agreement with previous equilibrium molecular dynamics simulations and experiments. Our aim here is to highlight the extraordinary power of TTCF, particularly for high-slip (low effective shear rate) systems, and to verify that equilibrium methods directly match NEMD measurements at experimentally accessible shear rates.

Article Details

Volume / Issue Vol. 164, Issue 23
Published June 21, 2026
ISSN 0021-9606
Publisher American Institute of Physics

Journal Info

The Journal of Chemical Physics

American Institute of Physics

ISSN: 0021-9606 Physical Sciences

Authors (7)

C

Carmelo Riccardo Civello

Department of Mathematics, Faculty of Science, Computing and Emerging Technologies, Swinburne University of Technology 1 , P.O. Box 218, Hawthorn, Victoria 3122,

L

Luca Maffioli

Department of Mathematics, Faculty of Science, Computing and Emerging Technologies, Swinburne University of Technology 1 , P.O. Box 218, Hawthorn, Victoria 3122,

E

Edward R. Smith

Mechanical and Aerospace Engineering, Brunel University London 2 , Kingston Lane, Uxbridge, UB8 3PH London,

J

James P. Ewen

Department of Mechanical Engineering, Imperial College London 3 , South Kensington Campus, London, SW7 2AZ London,

P

Peter J. Daivis

Department of Physics, RMIT University 5 , GPO Box 2476, Melbourne 3001, Victoria,

D

Daniele Dini

Department of Mechanical Engineering, Imperial College London 3 , South Kensington Campus, London, SW7 2AZ London,

B

B. D. Todd

Department of Mathematics, Faculty of Science, Computing and Emerging Technologies, Swinburne University of Technology 1 , P.O. Box 218, Hawthorn, Victoria 3122,