Complex Surface Tension in Surfactant-Stabilized Drops

• The paper investigates the thermodynamic framework governing drop shape transitions, revealing how effective surface tension shifts drive the transformation from spherical to elongated forms.
• It details the nonlinear interfacial equations of state and Marangoni stresses that modulate surfactant adsorption and dictate dynamic drop behavior.
• Experimental and computational studies, including tensiometry and phase-field modeling, validate the complex interplay between surfactant concentration, osmotic pressure, and flow.

Complex surface tension in surfactant-stabilized drops refers to the spatially and temporally varying interfacial tension that emerges when surfactant molecules adsorb at the interface between immiscible fluids (such as oil and water), causing both thermodynamic and hydrodynamic consequences far beyond a simple, constant reduction in surface tension. This phenomenon governs a broad spectrum of nonlinear drop morphologies, interfacial instabilities, flow transitions, and phase behavior in emulsions, foams, thin films, and microfluidic systems. The interplay between the intrinsic (bare) interfacial tension of the fluid interface and the osmotic, rheological, and Marangoni stresses arising from inhomogeneous surfactant coverage lies at the core of this complexity, producing effects ranging from kinetic stabilization against coalescence to dramatic shape transitions and dynamic rigidification or mobilization of the drop interface.

1. Thermodynamic Framework for Surfactant-Modulated Interfacial Tension

The total surface free energy of a droplet with surfactant is governed by both the bare oil–water (or liquid–air) interfacial tension, γ, and the two-dimensional osmotic pressure Πₛurf of the bound surfactant, yielding

$$F_\mathrm{surf}(A,N) = \gamma\, A - \Pi_\mathrm{surf}(A,N)\,A$$

where A is surface area and N is the number of adsorbed surfactant molecules. The osmotic pressure term derives from the free energy Fₛₐ via Πₛurf = -(\partial F_\mathrm{surfads}/\partial A)_N. In the Gibbs grand canonical ensemble, this physics is captured via a partition function parameterized by a conjugate surface tension σ and surfactant chemical potential μ, leading to a sharp phase transition when the “effective” surface tension

$\gamma_\mathrm{eff} = \gamma - \Pi_\mathrm{surf}$

crosses zero. At the critical bulk surfactant density $n_c(T)$, this transition manifests as a discontinuous change in drop morphology from a sphere (surface minimization, γₑff > 0) to an elongated “worm-like” shape (surface maximization, γₑff < 0) (Abacousnac et al., 12 Apr 2025).

2. Interfacial Equations of State and Surfactant Adsorption

The relationship between local surface tension and surfactant surface coverage Γ is typically nonlinear. For many systems, a Langmuir-type equation of state applies: $$\sigma(\Gamma) = \sigma_0 + RT\,\Gamma_\infty \ln(1 - \Gamma/\Gamma_\infty)$$ where σ₀ is the clean tension, Γ∞ the maximal packing, and RT the thermal energy scale (Constante-Amores et al., 2021, McDougall et al., 3 Feb 2026). In the dilute Γ≪Γ∞ regime, this reduces to σ ≈ σ₀ – E_s Γ, with E_s the surface elasticity or Marangoni modulus. The equilibrium surface coverage is controlled by bulk concentration c via an adsorption isotherm such as

$\Gamma(c) = \Gamma_\mathrm{max} \frac{Kc}{1 + Kc}$

with K the adsorption coefficient. When accounting for finite drop size and bulk solubility, mass-balance and corrected Gibbs relations become essential for accurate predictions of Γ, γ, and system stability (Urbina-Villalba et al., 2014).

3. Marangoni Stress, Hydrodynamics, and Drop Dynamics

Gradients in local surfactant concentration generate Marangoni (tangential) stresses, entering the interfacial boundary conditions: $$\mathbf{t}\cdot\boldsymbol{\sigma}\cdot\mathbf{n} = \nabla_s \sigma(\Gamma)$$ where ∇ₛ denotes the surface gradient (Constante-Amores et al., 2021, Marin et al., 2015). These stresses oppose or augment bulk capillarity and modify the flow field around drops, giving rise to phenomena such as:

• Delayed or arrested coalescence due to interfacial rigidification (Bruning et al., 2018, Constante-Amores et al., 2021).
• Enhanced or inverted Marangoni circulations in evaporating droplets, with surface tension gradients sustained far below chemical equilibrium (Marin et al., 2015).
• Tip-streaming jet ejection from drop tips at moderate capillary numbers and high surfactant elasticity (Adami et al., 2010).
• Rigidification or mobilization of the interface depending on the surfactant's rheological properties (Bhamla et al., 2016).

4. Phase Transitions and Morphological Instabilities

The competition between intrinsic surface tension and surfactant osmotic pressure is sufficient to drive first-order morphological transitions in drop behavior. The critical surfactant concentration $n_c(T)$, typically following the form

$n_c(T)/n_0 = (T_0/T)^{1/2} \exp(-T_0/T)$

separates spherical from worm-like morphologies, with the transition line tunable by temperature or surfactant chemical potential. This equilibrium phase behavior is validated experimentally in lipid-stabilized silicone oil droplets, where small changes in [SDS] or temperature toggle droplet shape, and signatures such as pearling instabilities and “worm-star” formation confirm the underlying osmotic mechanism (Abacousnac et al., 12 Apr 2025).

5. Experimental and Computational Methodologies

A range of experimental and computational approaches elucidate the dynamics of complex surface tension:

• Pendant-drop tensiometry and Young–Laplace profile analysis provide direct measurement of γ and associated adsorption isotherms (Harikrishnan et al., 2016).
• High-speed 3D particle-tracking velocimetry resolves spatiotemporal surface-tension gradients and interfacial velocities, e.g., in evaporating drops (Marin et al., 2015).
• Hybrid interface-tracking/level-set and phase field (Cahn–Hilliard) numerical methods fully couple surfactant transport, nonlinear equations of state, interfacial forces, and flow, capturing coalition, tipstreaming, Marangoni-induced flow patterns, and phase transitions in the drop population (Soligo et al., 2018, Huet et al., 2023, Constante-Amores et al., 2021, Soligo et al., 2023).
• Direct quantification of apparent line tension versus true line tension in wetting measurements establishes that even trace background surfactant fundamentally alters observed drop contact angles through size-dependent adsorption and attendant γ(R) corrections (Staniscia et al., 2022).

6. Rheological and Flow-Topology Implications

Surfactant-induced surface-tension gradients modify not only individual drop shape but also the collective rheological and topological properties of emulsions:

• Bulk viscosity and normal-stress differences in electrohydrodynamic emulsions are markedly affected by the balance between charge convection and Marangoni stresses, with critical shear rates and electric field thresholds set by interfacial elasticity and surfactant diffusivity (Poddar et al., 2019, McDougall et al., 3 Feb 2026).
• In turbulent suspensions, surfactant-laden drops introduce local changes in the interface-vorticity alignment and shift the flow-topology statistics toward more elongational events at the interface—effects scaling with elasticity number and Weber number and modulating mixing, coalescence, and phase inversion (Soligo et al., 2023).
• Apparent interfacial rigidity induced by surfactant can prevent the formation of satellite drops in bubble bursting events, dramatically reducing aerosol production at critical concentrations well below the CMC, with this suppression traceable to Marangoni blockade rather than simple tension reduction (Pierre et al., 2022).

7. Complementary Effects: Line Tension, Nanoparticles, and Rheological Complexity

Complex surface tension effects are further complicated by additional interfacial phenomena:

• Apparent line tension measurements are susceptible to substantial artifacts due to size-dependent surfactant adsorption, masking the molecular-scale intrinsic τ_true with much larger τ_app contributions (Staniscia et al., 2022).
• The inclusion of nanoparticles in surfactant-laden systems enables both higher and lower γ than achievable with either constituent alone via synergistic surfactant transport, competitive adsorption, and nanoparticle-induced interfacial “inoculation” (Harikrishnan et al., 2016).
• Surface rheology (dilatational and shear elasticity/viscosity) couples with surfactant dynamics to produce rich time-dependent interface responses, including frequency-dependent complex surface tension and thin-film stabilization or destabilization (Bhamla et al., 2016).

---

References

• "Worm-like emulsion droplets" (Abacousnac et al., 12 Apr 2025)
• "Apparent line tension induced by surface-active impurities" (Staniscia et al., 2022)
• "Effects of interplay of nanoparticles, surfactants and base fluid on the interfacial tension of nanocolloids" (Harikrishnan et al., 2016)
• "Delayed coalescence of surfactant containing sessile droplets" (Bruning et al., 2018)
• "Surfactant-dependent contact line dynamics and droplet spreading on textured substrates: derivations and computations" (Gao et al., 2021)
• "Electrorheology of a dilute emulsion of surfactant-covered drops" (Poddar et al., 2019)
• "Nonlinear electrohydrodynamics of a surfactant-laden leaky dielectric drop" (McDougall et al., 3 Feb 2026)
• "Tipstreaming of a drop in simple shear flow in the presence of surfactant" (Adami et al., 2010)
• "Role of surfactant-induced Marangoni stresses in drop-interface coalescence" (Constante-Amores et al., 2021)
• "Surfactant-driven flow transitions in evaporating droplets" (Marin et al., 2015)
• "A Novel Algorithm for the Estimation of the Surfactant Surface Excess at Emulsion Interfaces" (Urbina-Villalba et al., 2014)
• "Effect of surfactant-laden droplets on turbulent flow topology" (Soligo et al., 2023)
• "Coalescence of surfactant-laden drops by Phase Field Method" (Soligo et al., 2018)
• "Simulation of high-speed impact of surfactant-laden drops" (Huet et al., 2023)
• "Interfacial mechanisms for stability of surfactant-laden films" (Bhamla et al., 2016)
• "Influence of surfactant concentration on drop production by bubble bursting" (Pierre et al., 2022)