Surface tension ($\gamma$)

Surface tension $\gamma$ is the tensile force per unit length across a line drawn in a liquid interface. For a clean interface with isotropy, it has the same value in every tangential direction: each side pulls the other along the surface, perpendicular to the line. A straight segment of length $\ell$ therefore transmits a force of magnitude $\gamma \ell$ from one interface. Here we assume constant $\gamma >0$, with fixed temperature and composition. The same coefficient measures the reversible work required to create surface area:The equality of units follows from $1\mathrm{J}=1\mathrm{N}\mathrm{m}$; the equality of the mechanical coefficients follows from the moving-rim argument below. ^{[ref] [ref]}

For a concrete example, take an ideal flat soap film spanning a rectangular wire frame with a straight sliding rim of length $\ell$. Pull the rim outwards in the plane of the film, perpendicular to the rim. Each of the film's two liquid–gas interfaces pulls the rim back towards the film with force $\gamma \ell$, so the outward holding force is $F=2\gamma \ell$. The two faces are two interfaces, each with its own surface tension. With the assumed toy value $\gamma =0.072\mathrm{N}{\mathrm{m}}^{-1}$ and $\ell =20\mathrm{mm}=0.020\mathrm{m}$, the required force is $2.88\mathrm{mN}$; one interface would contribute $1.44\mathrm{mN}$. This value of $\gamma$ is a stipulated model parameter, rather than a measured value for a particular soap solution.

Why does an interface resist being enlarged? Its molecules have a different environment from those in the bulk liquid; changing that environment has an energy cost. In the constant-surface tension model, reducing surface area lowers surface energy. This connects an easily measured rim force to the shape of a droplet, the Young–Laplace equation for its pressure, and the restoring mechanism of capillary waves. It is useful because one coefficient links forces, energy and geometry. ^[ref]

The moving-rim argument also explains the factor of two without memorizing a special formula. An outward displacement $dx$ creates surface area $\ell dx$ on each face, hence $d{A}_{\mathrm{total}}=2\ell dx$. Slow pulling without viscous dissipation requires work $Fdx=\gamma d{A}_{\mathrm{total}}$, giving $F=2\gamma \ell$. For one interface, the same argument gives $F=\gamma \ell$. A surface tension gradient instead produces Marangoni stress; the constant-$\gamma$ assumption excludes that additional effect.

Elastocapillarity is the coupling between surface tension and elastic deformation. Changing an interface can reduce surface energy while bending a compliant solid increases bending energy. Thin structures can deform strongly when the available capillary energy competes with their elastic resistance.

For an isotropic thin elastic plate, its bending stiffness iswhere $E$ is Young's modulus, $h$ the thickness and $\nu$ Poisson's ratio. A sheet bent over a size $L$ has curvature of order $1/L$. Its bending energy scales as $B(1/L^2)L^2\sim B$, while the capillary energy scale is $\gamma L^2$, where $\gamma$ is surface tension. Balancing these scales gives the elastocapillary lengthHere $B$ has units of energy and $\gamma$ energy per area, so the ratio has units of length squared. At fixed material properties, ${\ell }_{\mathrm{ec}}\propto h^{3/2}$; doubling thickness multiplies this length by $2\sqrt{2}$. Geometry and contact angle enter the more detailed wrapping criterion.^[ref]

The 2007 capillary origami experiments of Charlotte Py and collaborators provide a direct demonstration. They cut polydimethylsiloxane membranes, 40–80 micrometres thick, into planar shapes, placed them on a water-repellent support, added a water droplet and recorded the deformation during evaporation. Suitable triangular sheets folded around the shrinking drop into a tetrahedral enclosure. Smaller or stiffer sheets failed to close completely. The observations connected a controlled material thickness and initial shape to the competition between surface energy and bending energy.^[ref]

Capillary origami uses surface tension to fold an initially flat elastic plate around a liquid droplet. The fold trades a reduction in surface energy for an increase in bending energy; the initial outline helps select the resulting enclosure.^[ref]

In the 2007 experiments of Charlotte Py and collaborators, 40–80 micrometre polydimethylsiloxane sheets rested on a water-repellent support. A water droplet reached the sheet's corners; time-lapse images followed its evaporation. Squares changed from four-corner bending to a tubular shape; triangles formed a tetrahedron. Smaller or thicker sheets could bend and then reopen without enclosing the droplet.^[ref]

The elastocapillary length makes the thickness dependence quantitative. For fixed Young's modulus, Poisson's ratio and surface tension, the bending stiffness obeys $B\propto h^3$ and ${\ell }_{\mathrm{ec}}\propto h^{3/2}$. Halving thickness therefore reduces ${\ell }_{\mathrm{ec}}$ by $2^{3/2}$, about $2.83$. At the same lateral size $L$, the sheet elastocapillary energy ratio $\eta =\gamma L^2/B$ increases eightfold. This calculation explains why making the sheet thinner can enable wrapping even without changing its outline.

The sheet elastocapillary energy ratio compares available energy scales; it does not determine a unique final shape. The contact angle, initial liquid volume and outline also matter. A cut chosen to permit several panels to meet is a mechanical design choice, while the surface tension supplies the driving force. Predicting the complete folding path requires more than the scalar energy comparison. The same group reported cross-shaped templates forming a cube and flower-shaped templates forming nearly spherical enclosures.^[ref]

For a thin elastic plate interacting with a liquid interface, the elastocapillary length isHere $B$ is the plate's bending stiffness, $\gamma$ the surface tension, $E$ Young's modulus, $h$ the thickness and $\nu$ Poisson's ratio. In this plate convention $B$ has units $\mathrm{N}\mathrm{m}$, while $\gamma$ has units $N/m$; their ratio has units of length squared. The quantity is different from the gravity-dependent capillary length.

To see the balance, bend a portion of sheet of lateral size $L$ through an angle of order one. Its curvature is of order $1/L$, so its bending energy scales as $B(1/L{)}^2{L}^2\sim B$. The competing surface energy scale is $\gamma L^2$. Their dimensionless quantity is the sheet elastocapillary energy ratioLarge $\eta$ means that the capillary energy scale is large relative to the bending energy scale. This dimensional analysis does not make $L={\ell }_{\mathrm{ec}}$ a universal threshold for complete closure in capillary origami.^[ref]

For chosen illustrative values $B=10^{-9}\mathrm{N}\mathrm{m}$ and $\gamma =0.05N/m$, the elastocapillary length is $\sqrt{2\times 10^{-8}}\mathrm{m}\simeq 141\mu \mathrm{m}$. A lateral size of $0.1\mathrm{mm}$ gives $\eta =0.5$, whereas $1\mathrm{mm}$ gives $\eta =50$. These are calculations with assumed parameters, not measurements of a particular material.

Two ways of changing size give opposite trends. Keep Young's modulus, Poisson's ratio and surface tension fixed, and let $s=L/L_0$ be a lateral scale factor. At fixed thickness, $B$ stays fixed and $\eta /{\eta }_0=s^2$. If all sheet dimensions shrink in the same proportion, then $h/h_0=s$, $B/B_0=s^3$ and $\eta /{\eta }_0=s^{-1}$. Halving lateral size alone decreases the sheet elastocapillary energy ratio fourfold; halving lateral size and thickness together doubles it. The phrase "smaller sheets" therefore needs a specified thickness rule.

The length can also be calibrated directly: the 2007 study fitted an elastic-loop model to a long polydimethylsiloxane strip held by an ethanol meniscus. For their complete-wetting configuration, loop width was $0.89{\ell }_{\mathrm{ec}}$. The coefficient belongs to that model and geometry.^[ref]

The sheet elastocapillary energy ratio compares capillary and bending energy scales for an elastic sheet:Here $L$ is a lateral size, $\gamma$ is surface tension, $B$ is plate bending stiffness and ${\ell }_{\mathrm{ec}}$ is the elastocapillary length. The ratio is a dimensionless quantity: $\gamma L^2$ and $B$ both have units of energy.

The estimate assumes bending through an angle of order one over size $L$. Then curvature scales as $1/L$ and integrating the bending energy density over area $L^2$ gives a scale $B$. A comparable change in interfacial area gives a surface energy scale $\gamma L^2$. The same balance underlies the capillary origami study's definition of elastocapillary length.^[ref]

For an original calculation, take $B=2\times 10^{-8}\mathrm{N}\mathrm{m}$, $\gamma =0.05N/m$ and $L=1\mathrm{mm}$. Then $\eta =2.5$. Doubling $L$ alone raises it to 10. Halving thickness at fixed $L$, Young's modulus and Poisson's ratio divides $B$ by eight and raises $\eta$ to 20.

Changing all dimensions gives another trend. Write $s=L/L_0$ and keep material properties and surface tension fixed. At fixed thickness, $\eta /{\eta }_0=s^2$; at fixed aspect ratio $h/L$, $B/B_0=s^3$ and $\eta /{\eta }_0=s^{-1}$. Thus halving $L$ alone gives one quarter of the initial ratio, while halving both $L$ and $h$ doubles it.

A larger ratio makes capillary energy more competitive with bending, but does not establish complete closure. Actual wetting, contact angle, droplet volume, contact-line pinning and accessible deformation paths also matter. The ratio compares scales rather than computing the full change in energy along a folding path.

A curved free interface of a liquid near a solid wall. Its shape reflects surface tension, contact angle and hydrostatic pressure.

A static liquid meniscus rising on a perfectly wetting vertical wall has dimensional height $z_0=\sqrt{2}{\ell }_c$, where ${\ell }_c=\sqrt{\sigma /(\rho g)}$ is the capillary length. In capillary-length coordinates, write the interface $x=h(z)$. The Young–Laplace equation gives $h_{zz}/(1+h_z^2{)}^{3/2}=z$ and the first integral $h_z/\sqrt{1+h_z^2}=z^2/2-1$. The integration constant follows from the horizontal far bath, where $z\to 0$ and $h_z\to -\infty$. Zero contact angle means $h=h_z=0$ at the wall, so $z_0=\sqrt{2}$. Expanding and integrating there gives $h\sim (z-z_0{)}^2/\sqrt{2}$, the parabolic shape used to match an entrained film.

Drying through a fluid's thermodynamic critical point to avoid traversing an ordinary liquid-vapor coexistence interface. Removing the receding meniscus can protect compliant microstructures from capillary force during release. Solvent exchange and compatibility with the structure remain practical constraints.

A force caused by an interface through its contact-line surface tension and the capillary pressure it produces. Both contributions must be counted with consistent geometry. For a narrow, completely wetting liquid bridge between flat plates of wetted radius $a$ and gap $d\ll a$, the leading attractive force is $2\pi a^2\gamma /d$.

For two liquid half-spaces, integrate the magnitude of a dispersion attraction per area $A_H/(6\pi h^3)$ from a molecular contact cutoff $h_0$ to infinity. The work per area is $A_H/(12\pi h_0^2)$. Because cleavage creates two interfaces, each has the displayed surface tension in this approximation. The cutoff is essential: a continuum attraction cannot describe atomic overlap and short-range repulsion. This estimate neglects interfacial entropy, detailed density profiles and nonadditive interactions, and it cannot replace the metallic bonding or hydrogen bonds that dominate some liquids' cohesion.

For a film described by $dU=TdS+\gamma dA$, with material content fixed, its Helmholtz free energy obeys $dF=-SdT+\gamma dA$. Equality of mixed partial derivatives gives the Maxwell relation $(\partial S/\partial A{)}_T=-(\partial \gamma /\partial T{)}_A$. If $\gamma$ depends only on temperature, write ${\gamma }^{\prime }=d\gamma /dT$; thenHere $C_A$ is the heat capacity at constant area. For $C_A>0$ and ${\gamma }^{\prime }<0$, reversible adiabatic stretching cools the film, since $dT/dA=T{\gamma }^{\prime }/C_A<0$, although positive mechanical input can increase its total internal energy. If $\gamma ={\gamma }_0-\beta T$ with constant positive $\beta$ and constant $C_A$, the isentropic process satisfies $C_A\ln (T/T_{\ast })+\beta A=\mathrm{constant}$, where $T_{\ast }$ is an arbitrary fixed reference temperature. The work convention uses the total area conjugate to $\gamma$; a two-surface film described by projected area requires the corresponding factor of two consistently throughout.

For a liquid interface with constant surface tension $\gamma$, the mechanical surface energy is the energy associated with its surface area $A$:The additive constant $C$ does not affect the work needed to change surface area. For several interfaces with the same surface tension, $A$ means their total surface area; the two faces of a thin soap film must both be counted. ^[ref]

With assumed $\gamma =0.072\mathrm{N}{\mathrm{m}}^{-1}$, creating $1\mathrm{c}{\mathrm{m}}^2=10^{-4}{\mathrm{m}}^2$ of one interface requires $7.2\mu \mathrm{J}$ of reversible work. In the sliding-rim example illustrated under surface tension, moving a $20\mathrm{mm}$ rim outwards by $1\mathrm{mm}$ creates $4.0\times 10^{-5}{\mathrm{m}}^2$ across the two faces and requires $2.88\mu \mathrm{J}$. Independently, multiplying its $2.88\mathrm{mN}$ holding force by the displacement gives the same work.

This makes the preference for a round isolated droplet precise. At fixed volume $V$, the three-dimensional isoperimetric inequality giveswith equality for a sphere. Thus a sphere minimizes surface energy among sufficiently regular closed shapes enclosing the same volume, when surface tension is uniform and other energy contributions are negligible. This is an application of the isoperimetric inequality, rather than a proof of it. Walls, wetting, and gravity can change the shape by adding constraints or other energy contributions. ^[ref]

The differential rule $d{E}_{\mathrm{s}}=\gamma dA$ integrates directly to $E_{\mathrm{s}}=\gamma A+C$ because surface tension is constant. At fixed liquid volume, decreasing surface area can drive viscous fluid flow; the released surface energy can supply viscous dissipation. This mechanical bookkeeping is useful for predicting which deformations can release energy. If surface tension depends on temperature or composition, distinguish this mechanical area-work relation from the complete interfacial thermodynamic free energy and internal energy.

Wetting describes how a liquid maintains contact with or spreads over a solid in the presence of another phase, usually a gas. The competition among solid-liquid, solid-gas and liquid-gas interfacial free energies determines the equilibrium contact angle. For an ideal smooth homogeneous surface, Young's contact-angle equation giveswith $\theta$ measured through the liquid.

A small water drop that spreads over clean glass but beads on a water-repelling surface provides a simple comparison. A low contact angle indicates stronger spreading within this partial-wetting description; complete spreading requires considering the associated interfacial free-energy change. Surface roughness, contamination and contact-angle hysteresis can make the observed angle depend on preparation and whether the contact line is advancing or receding, so one photograph need not determine an equilibrium material constant.

A fully penetrated rough-surface wetting model with $\cos {\theta }^{\ast }=r\cos \theta$, where $r$ is actual wetted solid area divided by projected area. Roughness amplifies the planar surface tendency: it lowers contact angle when $\theta <90^{\circ }$.

The area ratio $r=A_{\mathrm{actual}}/A_{\mathrm{projected}}\ge 1$, not the projected fraction covered by particles. A hemispherical array of footprint coverage $\varphi$ on an exposed flat substrate has $r=1+\varphi$.

A liquid spreads macroscopically over a solid with zero macroscopic equilibrium contact angle. In the simple interfacial-energy model the threshold is ${\gamma }_{SV}-{\gamma }_{SL}={\gamma }_{LV}$, or zero spreading coefficient; positive coefficient favors a spreading film. Extrapolating a concentration-dependent liquid-vapour surface tension to this threshold assumes the other interfacial energies remain fixed and no intervening adsorption or micellization behavior changes the trend.

The excess interfacial free energy released per newly wetted area when a liquid film replaces a solid-vapour interface is the displayed coefficient. Positive $S$ favors macroscopic spreading in the ideal equilibrium surface-energy model, whereas negative $S$ permits a finite contact angle obeying Young's equation. Surface preparation, chemical composition and adsorption can change all three interfacial energies.

A Cassie–Baxter state is a composite wetting state in which a drop contacts solid asperities over area fraction $f_s$ and trapped gas over the remainder. An ideal area-weighted contact-angle model gives the displayed apparent angle from the smooth-solid Young's equation angle ${\theta }_Y$. Reducing solid contact can strongly enhance apparent hydrophobicity. Such a state can be metastable and can collapse when liquid enters the texture; roughness alone does not guarantee it.

The line where the interfaces of three phases meet, for example where a liquid film meets a solid and the surrounding gas. Continuum assumptions at this line can differ from the bulk lubrication approximation.

A contact line is pinned when its position remains fixed while the surrounding drop shape or contact angle changes. Surface heterogeneity or deposited material can sustain such a state over a range of angles. Evaporation from a pinned drop then requires compensating liquid motion and can generate outward solute transport associated with the coffee-ring effect. A model with a freely receding line and fixed contact angle has a different boundary constraint and need not predict the same flow direction.

The Bond number compares hydrostatic stress over a length $\ell$ with capillary stress. For a gas bubble of radius $a$ and negligible density, $\mathrm{Bo}=\rho g{a}^2/\gamma$. Small Bond number is the usual sufficient criterion for capillarity to resist gravity-scale deformations. For a clean spherical bubble with negligible internal dynamic viscosity settling in creeping flow, the capillary number is $\mathrm{Ca}=\mu U/\gamma =\mathrm{Bo}/3$. This last relation uses the clean-bubble terminal velocity; it does not apply at arbitrary Reynolds number or to an interface immobilized by surfactants.

For a broad, thin oil lens below a flat ice-water interface, define the effective interfacial energy $\Gamma ={\gamma }_{io}+{\gamma }_{ow}-{\gamma }_{iw}$ from the relevant surface tensions. At fixed oil volume $\mathcal{V}=Ah$, the energy is $E=\Gamma \mathcal{V}/h+\frac{1}{2}({\rho }_w-{\rho }_o)g\mathcal{V}h$. Setting its derivative to zero gives $h^2=2\Gamma /[({\rho }_w-{\rho }_o)g]$ for $\Gamma >0$, and the positive second derivative proves a minimum. A contact angle measured through the oil gives $\Gamma ={\gamma }_{ow}(1-\cos \vartheta )$. Oil-ice tension alone does not determine this effective quantity. Complete wetting or rough underside pockets invalidate this finite, flat-lens approximation. With densities $1025$ and $850\mathrm{kg}{\mathrm{m}}^{-3}$, a lens below $1\mathrm{cm}$ requires $0<\Gamma <0.08584\mathrm{N}{\mathrm{m}}^{-1}$.

The contact angle is measured through a liquid between its interface and a solid substrate at their contact line. In a gently sloping two-dimensional film on a horizontal substrate, its small-angle value is the limiting magnitude of the height slope. A zero contact angle means a tangential approach to the substrate; it does not imply bounded curvature, as the limiting zero-flux thermocapillary film profiles illustrate.

The planar interfacial balance ${\gamma }_{SV}-{\gamma }_{SL}={\gamma }_{LV}\cos \theta$ for the equilibrium contact angle through the liquid, with interfacial free energies defined per actual area. Roughness and retained air require additional area accounting.

A contact line can remain pinned while the observed contact angle changes between receding and advancing limits. Surface roughness and chemical heterogeneity can produce this hysteresis. Slowly adding or removing liquid from a sessile drop measures advancing and receding angles separately. An observed pinned angle therefore need not be the single equilibrium angle predicted by Young's equation for a perfectly homogeneous surface.

For a smooth homogeneous solid in thermodynamic equilibrium with a liquid and vapor, Young's equation relates their three interfacial tensions to the contact angle measured through the liquid. A virtual increase $dA$ of wetted area changes the interfacial free energy by $({\gamma }_{SL}-{\gamma }_{SV}+{\gamma }_{LV}\cos \theta )dA$. Setting this to zero proves the equation. Roughness, pinning and chemical heterogeneity can make an observed apparent angle differ from this ideal equilibrium value. The equation does not determine each solid–liquid surface energy from a contact angle alone, but two liquids on the same unchanged dry substrate allow the solid–vapor term to be eliminated by subtraction.

Subtracting Young's equations for two liquids on the same solid eliminates the unknown solid–vapour interfacial energy. The result assumes equilibrium angles and an unchanged solid surface. Swelling, contamination or contact-angle hysteresis can invalidate this cancellation as a comparison of equilibrium material properties.

Combining reversible work of adhesion per area, $W={\gamma }_{SV}+{\gamma }_{LV}-{\gamma }_{SL}$, with Young's equation gives the displayed relation. It assumes the same smooth equilibrium interfaces; contact-angle hysteresis and dissipation are not included.

The capillary number compares viscous stresses with surface tension. Its velocity and viscosity conventions must be specified; a porous-medium displacement can use the Darcy velocity and the viscosity of the displaced fluid.

Capillary pressure is the pressure jump produced by a curved interface. For constant surface tension it is proportional to the interface curvature according to the Young–Laplace equation.

The threshold pressure difference required for a nonwetting fluid to invade a connected pore throat. In a simple cylindrical throat of radius $r$, the magnitude is approximately $2\gamma |\cos \theta |/r$, using surface tension $\gamma$ and contact angle $\theta$. A sealing caprock can fail to retain a buoyant phase if its pressure exceeds this threshold or if connected fractures bypass the fine pore network.

A surfactant is an interfacially active substance that changes surface tension. A concentration gradient therefore generates tangential Marangoni stress.

A cationic surfactant with a sixteen-carbon hydrocarbon tail and a positively charged trimethylammonium head group, accompanied by a ${\mathrm{Br}}^{-}$ counterion. It is an amphiphile that can form aqueous micelles. Its head-group charge has the opposite sign to that of sodium dodecyl sulfate.

An anionic surfactant with a twelve-carbon hydrocarbon tail and a negatively charged sulfate head group, accompanied by a sodium counterion. This separation of hydrophobic and hydrophilic regions makes it an amphiphile. It forms aqueous micelles under suitable concentration and temperature conditions and is used in sodium dodecyl sulfate polyacrylamide gel electrophoresis.

The material derivative of surface concentration combines tangential and normal area change, surface diffusion, and relaxation towards a reservoir concentration $C_0$. The rate $k$ represents exchange; setting $k=0$ recovers transport without exchange. The concentration is measured per interfacial area, so stretching dilutes it.

Consider a nearly spherical bubble of radius $a$ translating with velocity $\mathbf{U}$ through a liquid of dynamic viscosity $\mu$. Neglect internal gas viscosity and use exterior Stokes flow. Let the surfactant concentration be $C=C_0+C^{\prime }$ and linearize surface tension as $\gamma ={\gamma }_0-{\gamma }_1{C}^{\prime }$ with ${\gamma }_1>0$. Surface transport has surface diffusion coefficient $D_s$ and exchange-relaxation rate $k$.

In the bubble's translating frame, write its tangential surface velocity aswhere $\mathbf{n}$ points out of the bubble. The clean-bubble value is $A=-1/2$; a completely immobile interface has $A=0$. The dipolar surfactant distribution is $C^{\prime }=B\mathbf{U}\cdot \mathbf{n}$, withCombining that transport solution with the tangential Marangoni stress balance givesIncreasing $\lambda$ suppresses interfacial motion, whereas stronger diffusion or exchange reduces the concentration contrast and restores mobility. The linearization separately requires $|C^{\prime }|/C_0\ll 1$, or $2|A||\mathbf{U}|a/(k{a}^2+2D_s)\ll 1$. Small Reynolds number and capillary number are also needed for Stokes flow and a nearly spherical shape; a large immobilization parameter alone does not establish those assumptions.

The dipolar normal traction is balanced by buoyant hydrostatic pressure and the concentration-dependent normal capillary stress, in addition to the constant gas pressure and equilibrium Laplace pressure. A degree-one shape perturbation translates a sphere rather than deforming it; the balance fixes the terminal speed. A uniform gas pressure alone cannot balance this mode.

Linear steady transport on a translating sphere gives $B=2A{C}_{0}a/(k{a}^2+2D_s)$. The components of its unit normal vector are degree-one spherical harmonics. The condition $2|A|Ua/(k{a}^2+2D_s)\ll 1$ ensures the concentration perturbation is small relative to $C_0$.

Outward interfacial flow from a thinning region dilutes surfactant there and raises its surface tension. The resulting Marangoni stress pulls fluid back toward the thinning region, opposes further stretching, and slows rupture.

Consider a long-wave film-rupture model whose strong-Marangoni stress asymptotic growth law isHere $h_0$ is the base film thickness, $\mu$ is dynamic viscosity, $\gamma >0$ is the base surface tension, $V>0$ is the model's destabilizing interaction coefficient with dimensions of energy, and $K=k{h}_0$ is dimensionless wavenumber. The precise relation of $V$ to a microscopic interaction depends on the chosen film model.

Set $y=K^2$. Differentiating $y(1-\Gamma y)$ gives $1-2\Gamma y$, so the unique unstable-band maximum isThe second derivative is negative. Since $\Gamma \gg 1$, the selected $K_m\ll 1$ is consistent with the long-wave assumption. The growth timescale is not itself a rupture time or bubble lifetime: a perturbation starting at amplitude $\delta h_0$ must grow to a specified nonlinear scale, giving an additional logarithmic dependence on the initial disturbance, and the base thickness may evolve before that happens.

An insoluble surfactant remains on an interface. Without surface diffusion, its material concentration changes only through the local surface-area dilation rate.

Surface diffusion transports particles along an interface down a concentration gradient. On an isotropic smooth interface with constant diffusivity $D_s$, its flux is $-D_s{\nabla }_{s}C$ and its concentration source is $D_s{\Delta }_{s}C$, using the surface gradient and surface Laplacian.

For concentration $C$ per unit surface area, conservation on a moving material interface givesThe material derivative follows the surface particles. The two dilution terms are tangential expansion and curvature-induced area change from normal velocity $u_n$; the last term is surface diffusion.

For radius $a$ and tangential velocity ${\mathbf{u}}_s={\mathbf{I}}_s\mathbf{Ax}$ with symmetric traceless $\mathbf{A}$, the surface divergence is $-3\mathbf{n}\cdot \mathbf{An}$. At small surface Péclet number $a^2|\mathbf{A}|/D_s$, neglecting advection of the concentration perturbation givesThe result follows because this quadratic is a degree-two spherical harmonic with surface Laplacian $-6/a^2$ times itself and zero spherical mean.

Soap consists of sodium or potassium salts of long-chain fatty acids. Its hydrophilic carboxylate head and hydrophobic hydrocarbon tail make it an amphiphile; saponification of triglycerides produces these salts.

A detergent is a substance or formulation that assists cleaning, commonly through surfactants with water-compatible and oil-compatible molecular regions. These molecules collect at interfaces and can form micelles that help disperse oily material in water. Agitation, concentration and the composition of both the soil and washing solution affect the process.

A simple comparison is to shake an oil-water mixture with and without a suitable detergent and observe how readily droplets disperse and how long they persist. Lowering surface tension can improve wetting, while adsorption and aggregation change how oily material is carried away. Stable dispersion and removal are different steps: the surrounding fluid must still transport the dispersed material.

Detergents can also disrupt a lipid bilayer. In a protease protection assay, removing that membrane barrier tests whether it caused the apparent protection of an enclosed protein.

The Marangoni effect is interfacial flow driven by a surface-tension gradient. Tangential stress balances the surface gradient ${\nabla }_s\gamma$.

For a clean spherical gas bubble with negligible internal dynamic viscosity and thermal conductivity, a weak temperature gradient produces the displayed migration velocity, where ${\gamma }_T=d\gamma /dT$. The exterior Papkovich–Neuber representation can use only the scalar potential $\chi =-a^3\mathbf{U}\cdot \mathbf{r}/(2r^3)$. Its surface normal velocity equals $\mathbf{U}\cdot \mathbf{n}$, and its tangential traction balances minus the surface surface tension gradient from the temperature dipole around an insulating sphere. For the usual ${\gamma }_T<0$ the bubble moves toward warmer liquid. Small Reynolds number, Péclet number, capillary number, and fractional surface-tension variation justify the approximation; surfactants can oppose the temperature-driven Marangoni stress.

Across a fluid interface, the jump in normal traction is the surface tension times curvature, while the tangential traction jump is the surface gradient of surface tension. For one common normal convention,

For an almost spherical inviscid bubble in exterior Stokes flow $\mathbf{Ex}$, let $\gamma ={\gamma }_0-{\gamma }_1(C-C_0)$ and $K=C_0{a}^2/(2D_s)$. The quadrupolar surfactant distribution on a spherical interface makes the tangential velocity ${\mathbf{u}}_s=\alpha {\mathbf{I}}_s\mathbf{Ex}$ satisfyThe Marangoni stress approaches the stress needed to suppress surface motion as $M$ grows. The quadrupolar shape perturbation $r=a(1+\mathbf{n}\cdot \mathbf{Dn})$ isThese expressions require small deformation and small surface Péclet number; large $M$ alone does not justify a small concentration perturbation.

Chemophoresis is the motion of a particle or droplet driven by a gradient in chemical concentration. For a surfactant-coated bubble, adsorption converts the chemical gradient into a Marangoni stress.

The Rayleigh–Plateau instability makes a sufficiently long cylindrical fluid thread unstable to axisymmetric radius perturbations. Surface tension lowers area by amplifying wavelengths longer than the circumference-scale threshold.

A thin liquid lining inside a rigid cylindrical tube has pressure perturbation $p_1=-\sigma (a^{-2}{h}_1+h_{1,xx})$ at leading order in $h_0/a$. Lubrication theory gives flux $q=-h_0^3{p}_x/(3\mu )$, so mass conservation gives the displayed growth rate. Long waves with $0<k^2{a}^2<1$ grow, while short waves are stabilized by curvature. The maximum occurs at $k^2{a}^2=1/2$.

For insoluble surfactant with $\sigma ={\sigma }_0-A\Gamma$ and positive base tension, let $q=k^2{a}^2$, $\lambda =h_0^3({\sigma }_0-A{\Gamma }_0)/(3\mu a^4)$ and $\alpha =h_{0}A{\Gamma }_0/(\mu a^2)$. The coupled film-height and advected-surfactant amplitudes give the displayed growth equation. The unstable band remains $0<q<1$. Weak surfactant elasticity gives ${\beta }_{\max }=\lambda /4-3\alpha /8+O({\alpha }^2/\lambda )$; strong elasticity gives ${\beta }_{\max }\sim \lambda /16$. The Marangoni stress reduces tangential surface motion; strong elasticity reduces lubrication mobility to the value for an effectively immobile interface.

The Young–Laplace equation relates the pressure jump across a stationary interface to its curvature and surface tension. Choose a unit normal vector $\mathbf{n}$ pointing from the inside region to the outside region, and define $\kappa ={\nabla }_s\cdot \mathbf{n}$ to be positive for a convex sphere. In mechanical equilibrium, with constant surface tension and no additional normal viscous stresses,Here $k_1,k_2$ are the signed principal curvatures in this convention. Their sum is twice the corresponding signed mean curvature; stating the convention prevents a factor-of-two ambiguity. The Young–Laplace equation is a local normal stress balance for one interface. ^[ref]

A spherical liquid droplet of radius $R$ in a gas has one interface, with $k_1=k_2=1/R$. An ideal thin soap bubble has two liquid–gas interfaces. Adding their pressure jumps givesFor the soap bubble, $\Delta p$ now means the inner gas pressure minus the outer gas pressure. The second formula assumes equal surface tensions on both faces and film thickness negligible relative to $R$, so their radii can be identified. Retaining distinct radii gives $2\gamma /R_{\mathrm{inner}}+2\gamma /R_{\mathrm{outer}}$. A gas bubble inside bulk liquid has only one interface and therefore uses $2\gamma /R$. ^[ref]

With the assumed constant toy surface tension $\gamma =0.072\mathrm{N}{\mathrm{m}}^{-1}$, a droplet with $R=1\mathrm{mm}$ has excess pressure $144\mathrm{Pa}$; at $R=0.5\mathrm{mm}$ it has $288\mathrm{Pa}$. A thin soap bubble assigned the same $\gamma$ and $R=1\mathrm{mm}$ has excess pressure $288\mathrm{Pa}$. These are differences from the surrounding pressure, rather than absolute pressures. Halving the radius doubles the capillary pressure; doubling the interface count also doubles it at fixed radius and surface tension.

The spherical case has a short work argument. A virtual increase in radius gives $dV=4\pi R^{2}dR$ and $dA=8\pi RdR$. At mechanical equilibrium, the net pressure work balances the increase in surface energy:For a thin soap bubble, both interfaces gain approximately the same surface area, doubling the right-hand side and giving $4\gamma /R$. This is a virtual work balance between neighboring spherical configurations, rather than an assumption that a fixed-volume droplet actually expands. The general Young–Laplace equation replaces the spherical area-to-volume ratio by local curvature, making geometry directly measurable through pressure.

For an axisymmetric graph $z=h(r)$ with the normal pointing out of the liquid, the sum of principal curvatures is the displayed expression. In the small-slope limit it becomes $\kappa =-(h^{\prime \prime }+h^{\prime }/r)$, giving the capillary pressure $p_{\mathrm{liquid}}-p_{\mathrm{air}}=\gamma \kappa$. Regularity at the axis gives $h^{\prime }(0)=0$ and $h^{\prime }/r\to h^{\prime \prime }(0)$: the azimuthal contribution does not disappear at the axis. For example, the paraboloid $h=H(1-r^2/R^2)$ has $\kappa \approx 4H/R^2$, twice the meridional contribution.

The Gibbs--Thomson relation shifts the equilibrium chemical potential or composition at a curved interface by an amount proportional to its mean curvature. Small droplets therefore have a larger equilibrium solubility than flat interfaces.

For the stipulated spherical solid–liquid capillarity model, $T_m(r)=T_m(\infty )[1-2\gamma /(r\Delta H_V)]$, where $\gamma$ is the interfacial free energy and $\Delta H_V$ the latent heat density. It follows by matching $g_l-g_s\simeq \Delta H_V[T_m(\infty )-T]/T_m(\infty )$ to $2\gamma /r$. Support contacts, free surfaces, anisotropy, premelting and atomistic structure can modify this relation. A continuum extrapolation to negative melting temperatures is not physical.

Curvature depresses the equilibrium melting temperature of solid convex into liquid. Here $\mathcal{K}$ is the sum of principal curvatures, equal to twice a common convention for mean curvature, and $T_m$ is absolute temperature. The coefficient uses solid-liquid surface energy, mass density and latent heat. This stabilizes fine corrugations of a solidification front.

The capillary length measures when gravity and surface tension are comparably important in shaping a liquid interface. For a liquid of uniform mass density $\rho$, beneath a gas whose mass density is negligible, with uniform gravitational acceleration $g>0$ and constant surface tension $\gamma >0$,A deformation of height $a$ and horizontal length scale $L$ has hydrostatic pressure variation of order $\rho ga$ and small-slope capillary pressure of order $\gamma a/L^2$. Equating these gives $L\sim {\ell }_c$: surface tension smooths small features, whereas gravity matters for broad ones. For two resting fluids, replace $\rho$ by their positive mass density difference $\Delta \rho$ in this static balance. This change alone does not give the two-fluid wave dynamics, because both fluids can have inertia. ^[ref]

With the stipulated illustrative parameters $\gamma =0.072\mathrm{N}{\mathrm{m}}^{-1}$, $\rho =1000\mathrm{kg}{\mathrm{m}}^{-3}$, and $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$, the capillary length is ${\ell }_c=2.71\mathrm{mm}$. These parameters define the ideal calculations below; they are not a measurement of a particular sample.

For a sinusoid with wavenumber $k>0$, the relevant horizontal scale is $L=1/k$, while its wavelength is $\lambda =2\pi /k$. Thus the Bond number for this scale and the reciprocal restoring-term ratio areEquality is at $k_c=1/{\ell }_c$, hence ${\lambda }_c=2\pi {\ell }_c=17.02\mathrm{mm}$ for these parameters. Using $L=\lambda$ instead produces a different numerical Bond number; state the length convention. The capillary length itself is not the equality wavelength. This distinction helps choose whether a measured ripple should be modeled as a surface gravity wave, a nearly pure capillary wave, or a gravity-capillary wave. ^[ref]

A dynamic meniscus is a curved transition region joining a moving thin film to a bulk fluid reservoir.

In a thin wall-withdrawal model with large dimensionless surface-tension parameter $\hat{C}$, matching a parabolic static meniscus fixes inner longitudinal scale $ϵ={\hat{C}}^{-1/3}$ and thickness scale $\delta ={ϵ}^2$. The leading inner equation is $v_0{H}^{\prime }+(H^3{H}^{(3)}{)}^{\prime }/3=0$, matching $H\sim {\bar{z}}^2/\sqrt{2}$ below and a bounded film above. Scaling $H=v_0^{2/3}g$, $\bar{z}=v_0^{1/3}\zeta$ preserves the quadratic matching coefficient and gives the universal equation $3g^{\prime }+(g^3{g}^{(3)}{)}^{\prime }=0$. Its selected bounded solution tends to $0.948$, giving the displayed dimensionless law. It applies to positive slowly varying wall speed in the leading matched regime; the numerical constant comes from that universal inner problem, not a dimensional identity.

A capillary wave is a wave along a fluid interface whose restoring stress is predominantly supplied by surface tension. A curved free surface requires a pressure jump by the Young–Laplace equation; increasing its corrugation also costs surface energy. In the ideal limit, inertia carries the interface through its flat position and the resulting oscillation propagates. When gravity also contributes appreciably, use the combined gravity-capillary wave model. Within the capillary wave regime, shortening a ripple makes it travel faster. ^[ref]

The laser-diffraction measurement of capillary waves below turns this mechanism into a test: set the frequency, measure the wavelength independently with light, and ask whether the dispersion relation follows. The theory is useful both for predicting how ripple wave packets separate and for inferring surface tension from their motion.

A harmonic velocity potential decaying with depth, together with linearized kinematic and Bernoulli boundary conditions and capillary pressure, gives ${\omega }^2=gk+(T/\rho )k^3$. The corresponding phase speed combines long-wave gravity and short-wave capillarity contributions.

For a nearly flat interface, the extra surface energy in the small-slope approximation is $(\gamma /2)\int |\nabla h{|}^{2}dA$. Each nonzero Fourier mode therefore has thermal variance $⟨|h_{\mathbf{q}}{|}^2⟩=k_{B}T/(A\gamma q^2)$ in a convention with height reconstructed as the sum of its mode amplitudes. Summing over a two-dimensional annulus of wavevectors givesThe smallest wavevector is set by the lateral size or pinning, and the largest by the breakdown of the continuum description, commonly a molecular length. A fixed mean height removes the zero mode. Gravity supplies an additional restoring term and changes the long-wavelength cutoff when it is retained. Roughness here means a height variance; its square root is the root mean square amplitude, even when the mean height fluctuation is zero.

A gravity-capillary wave has restoring contributions from both gravity and surface tension. Consider an initially resting, horizontally unbounded liquid layer of constant depth $h$, with uniform mass density $\rho$, incompressible flow, inviscid flow and irrotational flow. The upper gas exerts constant pressure and has negligible inertia. Assume a clean interface with constant surface tension $\gamma$: no surfactant stresses, imposed current or wind. The wave amplitude $a$ must satisfy $ka\ll 1$ and $a\ll h$, so linearized boundary conditions apply. Deep water means $kh\gg 1$, rather than a particular depth in metres. ^[ref]

Let $z=0$ be the undisturbed free surface, with liquid below it. Potential flow gives a velocity potential $\varphi$ with $\mathbf{u}=\nabla \varphi$, and incompressible flow gives the Laplace equation. The linearized kinematic boundary condition for a free-surface graph and the Unsteady Bernoulli equation, including gravity and the Young–Laplace equation, giveHere $\eta (x,t)$ is the interface displacement; the signed small-slope pressure jump is $p_{\mathrm{liquid}}-p_{\mathrm{gas}}=-\gamma {\eta }_{xx}$. A wave crest has ${\eta }_{xx}<0$, so its capillary pressure is positive. Choose a right-travelling Fourier mode of positive wavenumber $k$ and angular frequency $\omega$:The factor $e^{kz}$ decays into the deep liquid, satisfies the Laplace equation, and makes ${\varphi }_z={\eta }_t$ at $z=0$. Substituting in the dynamic boundary condition leavesThis derivation displays the two restoring contributions with the same sign. Dimensional analysis also checks the result: $\gamma /\rho$ has units ${\mathrm{m}}^3{\mathrm{s}}^{-2}$, so both $gk$ and $(\gamma /\rho )k^3$ have units ${\mathrm{s}}^{-2}$. ^[ref]