🔬 condensed matter

Solid adsorption: the missing mechanism for surfactant contact lines -- a phase-field approach

This paper presents a thermodynamically consistent phase-field model for soluble surfactants that incorporates solid surface adsorption, revealing this mechanism as the key factor enabling predictive modeling of contact line dynamics by correctly capturing the experimentally observed trend of increased hydrophilicity across all contact angles.

Original authors: Parvathy K. Kannan, Kazi T. Iqbal, Diego Díaz, Ilse Mateman, Shahab Mirjalili, Gustav Amberg, Shervin Bagheri, Outi Tammisola

Published 2026-06-26

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Parvathy K. Kannan, Kazi T. Iqbal, Diego Díaz, Ilse Mateman, Shahab Mirjalili, Gustav Amberg, Shervin Bagheri, Outi Tammisola

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Picture: The Missing Piece of the Puzzle

Imagine you are trying to predict how a drop of soapy water spreads out on a table. Scientists have been building computer models to do this for years. These models are like sophisticated video games that simulate physics. However, the paper argues that all previous versions of this "game" were missing a crucial rule.

Because of this missing rule, the computer models kept making a specific mistake: they predicted that soap would make wet surfaces wetter and dry surfaces drier. But in the real world, experiments show that soap almost always makes surfaces wetter, regardless of whether they started out wet or dry.

This paper introduces a new "rule" called solid adsorption. By adding this rule to their computer model, the authors finally made the simulation match reality. They proved that for soap to behave the way it does in real life, it must stick to the solid surface (the table), not just float in the water or sit on the water's skin.

The Characters in Our Story

To understand how this works, let's break down the three main players:

The Drop (The Fluid): Think of a water droplet as a bouncy ball.
The Soap (The Surfactant): Imagine the soap molecules as tiny, two-faced magnets. One side loves water, the other side hates it. They want to sit right on the edge where water meets air.
The Table (The Solid Wall): This is the surface the drop sits on.

The Old Way vs. The New Way

The Old Model (The "Just Soap" Theory):
Previous scientists thought the soap only cared about the surface of the water (the skin of the drop).

The Analogy: Imagine the soap is a crowd of people only interested in standing on the roof of a house (the water surface). They pull the roof down, making the house flatten out.
The Problem: If the house was already leaning one way (hydrophobic/dry), pulling the roof down just made it lean more that way. If it was leaning the other way (hydrophilic/wet), it leaned even more. The model predicted that soap would exaggerate whatever the surface was already doing.
The Reality Check: In real life, soap acts like a universal "wetting agent." It makes dry surfaces wet and wet surfaces even wetter. The old model couldn't explain this.

The New Model (The "Solid Adsorption" Theory):
The authors realized the soap molecules also care about the table.

The Analogy: Now, imagine the soap molecules are like sticky notes. They don't just stick to the roof (water surface); they also stick to the walls of the house (the solid table).
The Mechanism: When the soap sticks to the table, it changes the "stickiness" of the table itself. It makes the table more inviting to the water.
The Result: Because the soap is actively changing the table's personality to be more "water-friendly," the drop spreads out. This happens even if the table was originally very dry and repelled water. This matches what we see in experiments.

The "Autophobing" Surprise

The paper also explains a weird phenomenon called autophobing. This is when you add soap, and instead of spreading out, the drop suddenly shrinks back up and becomes a tight ball.

The Analogy: Imagine the soap molecules are a group of people.
- In the normal case, they stand on the table under the drop, making the table friendly, so the drop spreads.
- In the "autophobing" case, the soap molecules are picky. They refuse to stand under the drop. Instead, they run to the dry table outside the drop and stick there.
The Result: By sticking to the outside, they make the outside of the table very slippery and "dry" for the water. The water feels like it's being pushed away from the edges, so it recoils and pulls itself into a tight ball. The paper's new model is the first to successfully simulate this "recoiling" behavior.

Why This Matters (According to the Paper)

The authors built a new mathematical framework (a "phase-field model") that includes these sticky interactions between the soap and the solid wall.

It Fixes the Math: It resolves the contradiction between what computers used to say and what scientists see in the lab.
It Explains the "Why": It shows that the reason soap makes things wetter is because it physically coats the solid surface, changing its energy.
It Captures the Weird Stuff: It can now predict both the spreading of drops and the shrinking (autophobing) of drops, depending on where the soap molecules decide to stick.

Summary

Think of the old models as trying to understand a dance by only watching the dancers' feet (the water surface). This paper says, "Wait, you have to watch how they hold onto the floor (the solid wall) too." Once you include the grip on the floor, the dance makes perfect sense, and the computer simulation finally matches the real-world performance.

Technical Summary: Solid Adsorption: The Missing Mechanism for Surfactant Contact Lines – A Phase-Field Approach

Problem Statement
Surfactants play a critical role in multiphase systems by modifying interfacial tension and wetting behavior. A robust experimental finding across diverse systems is that surfactants systematically reduce equilibrium contact angles, shifting both hydrophilic and hydrophobic surfaces toward increased wettability. However, existing numerical frameworks for soluble surfactants in moving contact line problems share a fundamental limitation: they neglect surfactant adsorption at solid walls. Current models, which only account for fluid-fluid interfacial tension reduction, predict that surfactants merely rescale Young's equation. This leads to a qualitative discrepancy with experiments: these models predict that hydrophilic substrates become more hydrophilic and hydrophobic substrates become more hydrophobic. Furthermore, they fail to explain phenomena such as "autophobing" (spontaneous droplet recoil), which requires asymmetric modification of solid-liquid and solid-gas interfacial tensions. The authors identify the absence of solid surface adsorption as the missing mechanism required to resolve this mismatch.

Methodology
The authors develop a thermodynamically consistent phase-field model for soluble surfactants in two-phase flows, incorporating both interfacial and solid surface adsorption. The framework is derived via variational principles consistent with the second law of thermodynamics.

Governing Equations: The model couples the incompressible Navier-Stokes equations with Cahn-Hilliard equations for the phase-field variable ( $\phi$ ) and surfactant concentration ( $\psi$ ). The bulk free energy includes mixing energy, surfactant entropy of mixing, interfacial adsorption energy, and a bulk penalty for surfactant solubility.
Solid Adsorption Mechanism: A key innovation is the formulation of wall free energy. The authors incorporate the Langmuir-Szyszkowski equation of state to model the reduction of solid-fluid surface tension as a function of local surfactant concentration. This allows for asymmetric adsorption strengths ( $\beta_{s1}$ and $\beta_{s2}$ ) on the solid-liquid and solid-gas interfaces.
Boundary Conditions: The variational derivation yields kinetic boundary conditions for both the phase field and surfactant concentration at the walls. Crucially, the inclusion of a gradient penalty term ( $|\nabla \psi|^2$ ) in the surfactant free energy elevates the surfactant equation to fourth order in space. This mathematical structure is essential for prescribing non-trivial normal gradients of surfactant concentration at the boundary, thereby communicating solid adsorption effects into the bulk fluid distribution.
Numerical Implementation: The system is solved on a staggered Cartesian grid using a second-order Adams-Bashforth scheme for temporal integration and a fractional-step method for the Navier-Stokes equations. A directionally dependent forward-backward splitting strategy is employed to handle the variable-coefficient Laplacian arising from degenerate surfactant mobility.

Key Results
The model was validated against analytical benchmarks for clean droplet spreading and compared against experimental data and prior numerical studies.

Resolution of the Wettability Discrepancy: When solid adsorption is included, the model predicts a universal shift toward increased hydrophilicity (reduced contact angles) across all initial contact angles. This aligns with experimental observations where both hydrophilic and hydrophobic surfaces become more wettable upon surfactant addition. In contrast, simulations without wall adsorption reproduced the erroneous trend of amplifying the substrate's existing wetting character.
Interfacial Curvature and Marangoni Stresses: The inclusion of solid adsorption reveals a subtle but significant bending of the fluid-fluid interface near the contact line. This curvature arises from a coupling mechanism: solid adsorption creates a surfactant concentration gradient near the wall, generating Marangoni stresses that drive bulk flow. Incompressibility constraints convert this flow into a local pressure field, which modifies the normal stress jump and forces the interface to adjust its curvature. This feature is absent in models lacking solid adsorption.
Prediction of Autophobing: By reversing the adsorption asymmetry (strong adsorption on the solid-ambient interface and weak adsorption on the solid-liquid interface), the model successfully predicts autophobing. In this regime, droplets spontaneously retract to larger equilibrium contact angles. The simulations showed contact angle increases of approximately $15^\circ$ to $19^\circ$ , consistent with experimental ranges reported in literature. This phenomenon is driven by the preferential accumulation of surfactants on the "dry" substrate ahead of the contact line, drastically reducing the solid-ambient surface tension relative to the solid-liquid tension.
Scaling Laws: The authors derived scaling laws for the thickness of the surfactant adsorption layer at the wall. They identified two regimes: one where the layer thickness depends on surfactant diffusion parameters (low gradient regime) and another where it depends on the surfactant strength at the wall (high gradient regime).

Significance and Claims
The paper claims that the inclusion of solid surface adsorption provides the "missing mechanism" necessary for the predictive modeling of surfactant-laden contact line dynamics. The authors assert that without this mechanism, numerical models cannot qualitatively reproduce the universal experimental trend of increased wettability or explain complex phenomena like autophobing.

The study establishes that:

Solid adsorption is not merely a secondary effect but a central driver of wetting dynamics that alters the fundamental balance of forces at the contact line.
The thermodynamic consistency of the proposed model, particularly the derivation of boundary conditions from a free energy functional including gradient penalties, is essential for capturing the coupling between surfactant distribution, Marangoni stresses, and interface curvature.
The model successfully bridges the gap between prior numerical predictions and experimental reality, offering a unified framework to explain both enhanced spreading and spontaneous dewetting based solely on the asymmetry of wall adsorption.

The authors conclude that their work resolves the qualitative discrepancy between prior formulations and experiments, demonstrating that solid adsorption is the critical factor required to accurately model surfactant-laden wetting phenomena.