Resolving Discrepancies in Disjoining Pressure Predictions for Liquid Nanofilms from Molecular Simulations

This paper resolves significant discrepancies in molecular simulation predictions of disjoining pressure for liquid nanofilms by demonstrating that neglecting long-range dispersion interactions and using inconsistent film thickness definitions in the original Peng method leads to errors, which are corrected by a revised approach that accounts for thickness-dependent surface tension effects and aligns with the Bhatt method.

The Gist

Imagine you have a very thin layer of water, so thin it's only a few molecules wide, sandwiched between two layers of air. This is called a nanofilm. You might find these in soap bubbles, oil droplets, or even in the tiny pores of rocks underground.

Scientists want to know how "sticky" or "stable" this thin layer is. They measure something called Disjoining Pressure. Think of this as the internal "tension" or "push-pull" force trying to either keep the film together or make it collapse.

The Problem: Two Scientists, Two Different Answers

For a long time, two different groups of scientists (let's call them Team Bhatt and Team Peng) were trying to calculate this pressure for water and argon (a gas). They used the same computer models and the same basic setup, but they got wildly different results.

• Team Bhatt said: "The pressure is moderate and stable."
• Team Peng said: "The pressure is huge! And it disappears completely at a certain thickness!"

It was like two chefs using the same recipe but ending up with one delicious cake and one burnt brick. The scientific community was confused. Why the huge difference?

The Investigation: Finding the Missing Ingredients

The authors of this paper (Yang, Zuo, et al.) decided to act as detectives. They realized Team Peng was missing two crucial "ingredients" in their recipe:

1. The Invisible "Long-Range" Glue (Dispersion Interactions)

Imagine molecules are like people at a party.

• Short-range: They can feel the person standing right next to them (a handshake).
• Long-range: They can also feel the general vibe of the whole room, even from across the dance floor (a shout or a song).

Team Peng only counted the handshakes (short-range forces) and ignored the room's vibe (long-range forces).

• The Fix: The authors added the "room vibe" back into the calculation.
• The Surprise: They found that these long-range forces act differently depending on how thick the film is.
  • In thick films: The long-range forces act like extra glue, making the surface tension stronger.
  • In very thin films: It gets weird. Because the film is so squeezed, these long-range forces actually push the molecules to spread out sideways. This lowers the surface tension.
  • The Analogy: Imagine a crowd of people in a hallway. If the hallway is wide, they hold hands tightly (high tension). If the hallway is super narrow, they get squished, and to avoid bumping heads, they spread their arms out sideways, loosening their grip (lower tension).

This "crossover" behavior was the key. Because Team Peng ignored this, their math was off.

2. The Ruler Problem (Defining Thickness)

Team Peng was measuring the thickness of the film with a ruler that was slightly broken.

• The Old Way: They measured the whole box size, including the empty air and some water that had evaporated into the air. It was like measuring a sandwich but including the crumbs that fell off the plate. This made the film look thicker than it really was.
• The New Way: The authors used a "thermodynamic ruler." They calculated the thickness based strictly on how many water molecules were actually in the liquid part, ignoring the air and the crumbs.

The Solution: Putting It All Together

When the authors fixed the recipe:

1. They added the "long-range vibe" (dispersion forces).
2. They used the "correct ruler" for thickness.

The Result?
Team Peng's new numbers suddenly matched Team Bhatt's numbers perfectly! The "burnt brick" turned into a "delicious cake."

They also calculated the Hamaker Constant, which is basically a score for how strong the attraction is between the molecules.

• Old Method (Peng): Gave a score that was too low or too high, depending on the system.
• New Method: Gave a score that matched the "gold standard" (Team Bhatt) and made physical sense.

Why Does This Matter?

This isn't just about math; it's about real-world applications.

• Oil & Gas: Understanding how liquids behave in tiny rock pores helps us extract oil more efficiently.
• Foams & Emulsions: It helps us make better soaps, creams, and food products that don't separate or collapse.
• Nanotechnology: As we build smaller and smaller machines, these tiny forces become the biggest players.

In a nutshell: The paper solved a scientific mystery by realizing that when you squeeze a liquid film very thin, the invisible forces between molecules change their behavior in a surprising way. By fixing how we measure the film and how we count those invisible forces, we finally got the right answer.

Technical Summary

1. Problem Statement

Liquid nanofilms (a few nanometers thick) are critical in systems like foams, emulsions, and georesources. Their stability is governed by disjoining pressure ($\Pi$). While molecular simulations are the primary tool for studying these systems, a significant discrepancy exists between two established methods for calculating $\Pi$:

• The Bhatt Method: Based on an integrated Gibbs-Duhem equation using chemical potential extrapolation.
• The Peng Method: Based on film thermodynamics, relating $\Pi$ to the derivative of surface tension ($\sigma_f$) with respect to film thickness ($h$).

Despite using identical system setups (liquid slabs between vapor phases), previous studies showed poor agreement. For water, the Peng method predicted disjoining pressures nearly double those of the Bhatt method at small thicknesses, with the Peng results dropping to zero much earlier than the Bhatt results. The origin of this large quantitative disagreement was previously unclear.

2. Methodology

The authors conducted Molecular Dynamics (MD) simulations using the LAMMPS code to resolve these discrepancies.

• Systems Studied: Pure water (SPC/E model) and pure argon.
• Simulation Conditions: NVT ensemble with periodic boundary conditions.
  • Water: Simulated at 479 K.
  • Argon: Simulated at 100.05 K and 105.93 K.
• Key Variables Investigated:
  1. Long-Range Dispersion Interactions: Two cases were compared:
     • Cutoff: Truncating Lennard-Jones (LJ) interactions at a finite distance (ignoring long-range tails).
     • PPPM: Using Particle-Particle Particle-Mesh (PPPM) to account for long-range dispersion interactions accurately.
  2. Film Thickness Definition:
     • Original Peng Definition: Based on system size from bulk simulations, susceptible to pressure uncertainty and evaporation artifacts.
     • Thermodynamically Consistent Definition (Ivanov & Toshev): Defined by $h = (N/A - \rho_v L_z) / (\rho_l - \rho_v)$, utilizing bulk liquid and vapor densities to define the dividing surfaces.

The disjoining pressure was calculated using the revised Peng equation: $\Pi = -2 (\partial \sigma_f / \partial h)$.

3. Key Contributions

The paper identifies and corrects two specific flaws in the original application of the Peng method:

1. Necessity of Long-Range Dispersion: The authors demonstrate that neglecting long-range LJ interactions leads to systematic errors. Crucially, they reveal a thickness-dependent crossover effect on surface tension:

   • At large thickness: Including long-range interactions increases surface tension (standard behavior due to stronger in-plane attraction).
   • At small thickness: Including long-range interactions decreases surface tension. This is caused by the enhanced disjoining pressure compressing the film normally, which induces lateral expansion. This expansion counteracts the in-plane contraction, netting a lower surface tension.
   • Impact: Since $\Pi$ is the derivative of $\sigma_f$ with respect to $h$, this non-trivial thickness dependence drastically alters the calculated $\Pi$.

2. Consistent Thickness Definition: The authors show that the original Peng definition of film thickness introduces systematic errors due to unknown pressure conditions and evaporation into vacuum regions. Adopting the thermodynamically consistent definition (based on bulk densities) aligns the Peng method results with the Bhatt method.

4. Results

• Water System (479 K):
  • When using long-range dispersion and the consistent thickness definition, the Peng method results align closely with the Bhatt method.
  • The crossover behavior in surface tension was observed at $\approx 15.5$ Å.
  • Hamaker Constants: The revised method yielded a Hamaker constant ($A_H$) of $3.2 \times 10^{-18}$ J, which is consistent with the Bhatt method ($3.0 \times 10^{-18}$ J). In contrast, the cutoff-based Peng method yielded a significantly lower value ($1.5 \times 10^{-18}$ J), indicating an underestimation of van der Waals attraction.
• Argon System:
  • Argon showed a systematic underestimation of surface tension when long-range interactions were ignored, but the crossover behavior seen in water was absent due to weaker disjoining pressure effects.
  • The revised Peng method (with long-range interactions and consistent thickness) agreed well with the Bhatt method.
  • Hamaker constants derived from the revised method ($3.3 \times 10^{-19}$ J at 100 K) matched the Bhatt method ($2.6 \times 10^{-19}$ J), whereas the cutoff/Peng method yielded values an order of magnitude higher ($1.2 \times 10^{-18}$ J).
• Scaling Laws: The revised simulations showed that $\Pi$ follows the expected $h^{-3}$ scaling predicted by Hamaker theory, whereas cutoff-based treatments deviated significantly from this scaling.

5. Significance

This study resolves a long-standing controversy in the simulation of nanofilms by establishing that methodological consistency is paramount.

• Physical Insight: It reveals a subtle coupling between long-range intermolecular forces, film thickness, and mechanical balance (normal compression vs. lateral expansion) that dictates surface tension behavior in nanofilms.
• Methodological Framework: It provides a corrected, robust framework for the Peng method, requiring:
  1. Accurate treatment of long-range dispersion (e.g., via PPPM).
  2. A thermodynamically consistent definition of film thickness.
• Implications: The findings ensure that predicted disjoining pressures and Hamaker constants are physically meaningful. This is critical for accurately modeling nanoscale interfacial phenomena in applications such as flotation, coating technologies, and the exploitation of georesources, where film stability is a governing factor.