Layering Theory of Liquids at Solid Interfaces: Interfacial Layering Oscillator Model

The structural organization of liquids near solid interfaces profoundly influences phenomena such as wettability, nanofluidic transport, and interfacial heat transfer. This study introduces the Interfacial Layering Oscillator Model (ILOM), a concise, semi-phenomenological framework that accurately captures the oscillatory density profiles of liquids adjacent to planar solid surfaces. By deriving a second-order differential equation rooted in classical statistical mechanics and calibrated with molecular dynamics simulations, ILOM predicts the amplitude, decay rate, and wavelength of interfacial density layering with exceptional computational efficiency. This versatile model applies to both hydrophilic and hydrophobic surfaces and extends to liquids beyond water, including methanol, providing valuable insights into critical interfacial properties that advance nanoscale fluid mechanics and material design.

💡 Research Summary

This paper presents a novel theoretical framework called the Interfacial Layering Oscillator Model (ILOM) to describe and predict the oscillatory density profiles of liquids near solid interfaces. Recognizing the profound impact of nanoscale liquid structuring on phenomena like wettability and nanofluidic transport, the authors address the need for a model that balances physical accuracy with computational efficiency, bridging the gap between rigorous but complex statistical mechanical theories and computationally expensive molecular dynamics (MD) simulations.

The core of the work lies in deriving a simplified, second-order ordinary differential equation (ODE) that captures the essential physics of interfacial layering. Starting from the foundational Yvon-Born-Green (YBG) equation of classical statistical mechanics, the authors apply linearization and short-range approximations to transform the integro-differential equation into a more tractable form. The resulting governing equation (Eq. 1.2) resembles a damped, driven harmonic oscillator for the normalized density deviation, ℎ(𝑧). This formulation intuitively incorporates terms representing the curvature of the density profile (inertia), energy dissipation due to interactions (damping), the periodic packing of molecules (oscillation), and fluid-fluid correlation effects (driving force).

The study introduces three progressive variants of the model: the simple harmonic oscillator (SHO) model with constant parameters yielding an analytical solution; the SHO1 model with distance-dependent damping and frequency parameters for greater accuracy; and the SHO2 model which includes an external forcing term. The position-dependent parameters in SHO1, such as 𝛾(𝑧) and 𝜔²(𝑧), are designed to decay exponentially from the wall into the bulk, capturing the localized nature of interfacial effects.

To calibrate and validate the ILOM, the authors performed equilibrium MD simulations using the LAMMPS package. Systems included water reservoirs over flexible graphene (hydrophobic) and rigid hexagonal boron nitride (hBN, hydrophilic) surfaces, as well as methanol over graphene. Simulations employed the SPC/E water model, Lennard-Jones plus Coulomb potentials, and were conducted in the NVT ensemble.

The key outcome is the successful parameterization of the ILOM against MD data. Optimized parameters for the SHO model and more complex functions for the SHO1 model are provided for each liquid-surface combination. The results demonstrate that the ILOM can accurately reproduce the amplitude, decay length, and wavelength of the interfacial density oscillations observed in detailed MD simulations, but at a fraction of the computational cost. The model’s versatility across different surface chemistries (hydrophilic/hydrophobic) and liquids (water/methanol) highlights its potential as a general and efficient predictive tool. This work provides a valuable theoretical bridge, offering deep physical insight and a practical method for analyzing and designing nanoscale interfacial systems in fluid mechanics and materials science.