Dynamic Contact Angle Hysteresis

• Dynamic contact angle hysteresis is the rate- and history-dependent difference between advancing and receding contact angles affecting wettability and fluid transport.
• It arises from the interplay of contact-line pinning, molecular friction, and hydrodynamic dissipation, with models like Cox–Voinov quantifying its dynamics.
• This phenomenon is critical for applications in multiphase flows, inkjet printing, and programmable surfaces where precise wetting control is required.

Dynamic contact angle hysteresis refers to the rate-dependent, history-dependent difference between advancing and receding contact angles during the motion of a three-phase contact line on a solid, critically influencing wettability, capillarity, and fluid transport in many applications. In contrast to static hysteresis, which arises from metastable energy barriers due to surface heterogeneity or roughness, dynamic contact angle hysteresis fundamentally couples microscopic dissipation, contact-line friction, and hydrodynamic effects. This phenomenon controls macroscopic processes ranging from wet granular mechanics and multiphase flows in porous media to “sticky-on-demand” surfaces and dynamic wetting under oscillatory or steady motion.

1. Fundamental Mechanisms of Dynamic Contact Angle Hysteresis

Dynamic contact angle hysteresis emerges due to the interplay of contact-line pinning, molecular-scale friction, surface heterogeneity, and hydrodynamic dissipation. When a contact line advances, it remains pinned until the local angle reaches an upper bound (advancing angle $\theta_a$); upon recession, it remains pinned until reaching a lower bound (receding angle $\theta_r$). The interval $[\theta_r,\,\theta_a]$ defines the intrinsic hysteresis window and is set by the energy cost of surmounting microscopic defects or chemical heterogeneities (Akbari et al., 2015).

As the contact line moves with nonzero velocity $U_\text{cl}$:

• Pinning and Depinning Dynamics: The contact line remains immobilized (pinned) so long as $\theta$ lies within $[\theta_r,\theta_a]$. Once a threshold is crossed, the line moves with a finite, velocity-dependent dynamic angle $\theta_d$ (Mani et al., 2015), with the transition manifesting as stick–slip motion (especially on heterogeneous or rough substrates) (Zhang et al., 2021).
• Hydrodynamic Dissipation: Viscous losses in the wedge-shaped region near the moving line (the Huh–Scriven paradox) produce nontrivial, rate-dependent variation in $\theta_d$, leading to separable advancing and receding dynamic branches (Shi et al., 2017, Bao et al., 2018).
• Frictional and Molecular Effects: Atomic- or molecular-scale friction at the contact line, revealed via large-scale MD on “absolutely smooth” versus atomically structured walls, establishes that both static and dynamic hysteresis vanish in the absence of lateral energy barriers, unifying both under a contact-line friction law $f_\text{fric} = \xi U$ (Liu et al., 2019).

2. Mathematical Models and Scaling Laws

Multiple mathematical frameworks describe dynamic contact angle hysteresis, distinguished by system geometry, surface structure (homogeneous/heterogeneous), and dominating physical effects:

• Cox–Voinov Law and Extensions: On smooth substrates and in the hydrodynamic regime, the dynamic angle obeys

$$\theta_d^3 = \theta_e^3 \pm 9\,\mathrm{Ca}\ln \left( \frac{L}{\ell} \right)$$

where $\mathrm{Ca} = \mu U/\gamma$ (capillary number); $+/-$ for advancing/receding contact lines (Akbari et al., 2015, Shen et al., 2024). Empirically, in LBM and MD, $\cos\theta_d = \cos\theta_e - k\,\mathrm{Ca}$ with $k\approx1$ universally across wetting regimes (Li, 2022, Mandrolko et al., 2023). At the nanoscale, the Cox–Voinov model matches MD results for droplets where rolling motion is weak (Ca $<0.1$), with deviations observed only at high Capillary numbers (Mandrolko et al., 2023).

• Power-Law and Empirical Correlations: For liquid bridges and capillary tubes, power-law relations such as

$|\cos\theta_{d}-\cos\theta_s| = K\,\mathrm{Ca}^n$

hold, with system-dependent $K$, $n$ (commonly $n\in[0.16, 0.7]$) depending on line geometry, dynamics, and dissipation mechanisms (Shi et al., 2017, Bao et al., 2018). For example, in liquid bridges, $n_\text{rec}=0.16$, $n_\text{adv}=0.17$, showing a much weaker rate-dependence than classical 1D moving contact lines (Shi et al., 2017).

• Frictional Law & Energy Balance: The dynamic angle for a contact line with friction $\xi$ satisfies

$\cos\theta_d = \cos\theta_0 - (\xi/\gamma) U$

highlighting the linear relation between angular deviation and contact-line speed set by the molecular-scale friction (Liu et al., 2019). Multiphase mechanical energy balance yields a general dynamic law (Harvie, 2022):

$$\cos\theta_d \approx \cos\theta_{\rm stat} + \frac{\mu U}{\sigma_{12}\ln\left(\tfrac{r}{h}\right)}$$

where $\sigma_{12}$ is interfacial tension; $r$, $h$ are capillary and roughness lengthscales.

• Rate-Independent Hysteresis (Pinning Interval): Geometric/variational minimization approaches, as in the one-phase Bernoulli free-boundary problem, define admissible macroscopic angles within a pinning interval, with interface motion only at threshold (Feldman et al., 2023).

3. Microscopic Origins and Surface Effects

• Role of Surface Heterogeneity and Roughness: Contact angle hysteresis is fundamentally set by local energy barriers at the contact line resulting from surface roughness or chemical inhomogeneity (Harvie, 2022, Kumar et al., 2023). On atomically “smooth” walls (no lateral force corrugation), both static and dynamic hysteresis vanish (Liu et al., 2019); even minimal atomic roughness is sufficient to induce pinning, convex nanobending, and dynamic angle deviation.
• Topography-Driven Dissipation: For periodic micro-roughness (e.g., pillar arrays), the dissipated energy per unit projected area during a contact-line jump varies as $\phi\ln\phi$ ($\phi$: pillar area fraction), and the angle hysteresis obeys

$$\cos\theta_\text{adv} = r\cos\theta_e - \overline{D} ,\quad \cos\theta_\text{rec} = r\cos\theta_e + \overline{D}$$

with $r$ the roughness ratio and $\overline{D}$ the normalized dissipation (Kumar et al., 2023). As $\phi\to0$, $\overline{D}\to0$ and hysteresis disappears.

• Dynamic vs. Static Hysteresis: The presence of rate-dependent terms and the scaling of frozen-in volume disorder with $\Delta\theta^3$ in granular media, foams, and emulsions confirms that dynamic hysteresis is not purely a generalization of static pinning but encodes the rate-dependent competition between interfacial “memory” and local energy dissipation (Mani et al., 2015).
• Chemically Patterned and Electrowetting Substrates: For surfaces patterned with spatial variations of equilibrium angle $\theta_Y(x)$, both the width and skewness of hysteresis are set by the local pattern and its statistical distribution. The macroscopic advancing and receding angles reflect coarse-grained, harmonic averages over the substrate pattern (Zhang et al., 2021, Xu et al., 2020). In electrowetting, static hysteresis can be tuned dynamically by a voltage-dependent energy barrier scaling as $U_0^2\sqrt{p/R}$ (electrode pitch, droplet radius), enabling programmable hysteresis (Wang et al., 2020).

4. Experimental and Computational Observations

Dynamic contact angle hysteresis has been extensively probed and modeled both experimentally and numerically.

• Granular and Wetting Systems: In wet granular beds, the variance of bridge volumes at late times saturates to a finite, $\propto \Delta\theta^3$ value, due to persistent pinning—even as local pressures equilibrate with a $t^{-3/2}$ power-law (without hysteresis) (Mani et al., 2015).
• Liquid Bridges and Drop Deposition: Cyclic stretching/compression of liquid bridges shows pronounced, rate-dependent hysteretic force loops, with the area scaling with loading rate and triple-line geometry. A key observation is the reduced yield strength and history-dependent cohesion when $\Delta\theta > 0$ (Shi et al., 2017, Akbari et al., 2015).
• Micro- and Nanoscale Droplets: MD of nanodroplets subject to body force demonstrates that macroscale laws (Cox–Voinov) remain accurate at small scales provided $Ca < 0.1$ and rolling “circulation” is weak. At higher $Ca$, rolling dominates, leading to failure of hydrostatic models and non-linear, enhanced hysteresis (Mandrolko et al., 2023).
• Surface Evolver and Numerical Geometry: Mesoscopic simulations using lattice Boltzmann, SPH, and Surface Evolver capture dynamic hysteresis as emergent from the competition between interface relaxation, contact-line friction, and wetting energetics (Colosqui et al., 2012, Kumar et al., 2023, Wen et al., 2017).
• Oscillatory Motion: Under high-frequency oscillatory forcing, the dynamic hysteresis loop (phase shift between angle and line-speed) is governed not by viscous dissipation but by the elastic/memory properties of molecular surface layers. Phase lag $\tau$ and hysteresis area correlate strongly with static hysteresis (layer flexibility), but not with contact-line friction (Shen et al., 2024).

5. Modeling Frameworks and Predictive Laws

A variety of first-principles and semi-empirical models account for dynamic hysteresis:

• Variational Principles and Energy Minimization: Onsager's principle and minimizing movement schemes provide transparent ODE/PDE-based laws for dynamic hysteresis on both homogeneous and heterogeneous substrates. These systematically unite pinning, depinning, and the evolution of the macroscopic angle as a function of velocity and substrate pattern statistics (Xu et al., 2020, Feldman et al., 2023, Zhang et al., 2021).
• Mesoscopic and Kinetic Models: Lattice Boltzmann–based, SPH, and molecular-kinetic models predict both the static and dynamic branches of the $\theta$–$Ca$ relation and connect mesoscopic structure (disjoining pressure, slip length, interfacial forces) to macroscopic observable hysteresis, with no ad hoc boundary conditions required (Colosqui et al., 2012, Bao et al., 2018, Wen et al., 2017).
• Mechanical Energy-Balance: Small control volume MMEB formulations yield predictive, physically grounded laws separating the effects of energy dissipation, inertia, and kinetic flow at the contact line. These models extend directly to include new physics such as substrate compliance or surfactant activity (Harvie, 2022).
• Boundary Condition Formulations: In continuum simulations, dynamic hysteresis can be encoded via generalized Hocking laws, complex mobility (frequency-dependent) extensions, and Navier-type slip or GNBC boundary conditions with empirically fitted parameters capturing phase lag and loop area in oscillatory regimes (Shen et al., 2024).

6. Applications, Implications, and Outlook

Dynamic contact angle hysteresis is essential to the prediction and control of multiphase flows, adhesion, and wetting-driven actuation across diverse fields:

• Granular Mechanics and Wet Soil: Persistent heterogeneity in capillary bridge volumes enables the “memory” and history-dependent rheology on which the macroscopic strength of wet granular materials rests (Mani et al., 2015).
• Manufacturing and Printing: In liquid transfer, inkjet, or coating processes, the magnitude and velocity-dependence of dynamic hysteresis determine drop detachment, patterning fidelity, and operational speed limits (Akbari et al., 2015, Shi et al., 2017).
• Micro- and Nanofluidics: Precise understanding of dynamic hysteresis enables quantitative modeling of droplet motion, shape, and stability on patterned, rough, or externally controlled (e.g., electrowetting) surfaces (Wang et al., 2020, Li, 2022).
• Programmable Surfaces: Dynamic (voltage-tunable) hysteresis opens access to real-time control of drop manipulation, adhesion, and “sticky-on-demand” behavior in microfluidics and soft robotics (Wang et al., 2020).
• Experimental and Simulation Protocols: Accurate measurement and implementation of dynamic contact angle boundary conditions in numerical models require explicit accounting for hysteresis via the appropriate geometric, physical, and energy-based frameworks (Wen et al., 2017, Shen et al., 2024).

7. Summary Table: Key Hysteresis Regimes and Scaling Laws

System/Mechanism	Scaling law / dynamic CAH behavior	Reference
Wet granular matter	$\sigma_v^2(\infty) \sim (\Delta\theta)^3$	(Mani et al., 2015)
Liquid bridge, line-radius	$|\cos\theta_d - \cos\theta_s| = K\,Ca^n$	(Shi et al., 2017)
Patterned/heterogeneous	Effective angle: harmonic average pinned/depinned	(Zhang et al., 2021, Xu et al., 2020)
Nanoscale MD	Cox–Voinov holds for $Ca<0.1$	(Mandrolko et al., 2023)
Rough pillar arrays	Dissipation $\propto \phi\ln\phi$, CAH $= (\theta_\text{adv}-\theta_\text{rec})$ from $D(\phi)$	(Kumar et al., 2023)
Oscillatory motion, elastic layer	Hysteresis loop area $\propto$ static CAH (not friction)	(Shen et al., 2024)

Dynamic contact angle hysteresis thus embodies a unified manifestation of molecular-scale details, energetic pinning, and hydrodynamic flow—predictable via first-principles or empirically calibrated continuum models, and essential for rational design and quantitative control of wetting-driven phenomena.