Al sumergir un aro en agua jabonosa y soplar suavemente, una delgada lámina líquida se estira, se estrangula y se relaja en milisegundos hasta formar una esfera iridiscente perfecta flotando en el aire. Si observamos con una lupa la espuma de jabón en un recipiente, descubriremos que lo que parece un caos desordenado obedece en realidad a una geometría inflexible: las paredes líquidas nunca se reúnen de cuatro en cuatro ni de cinco en cinco, sino exclusivamente de tres en tres a lo largo de canales líquidos formando ángulos exactos de $120^\circ$ ($2\pi/3$), y dichas aristas se cruzan exclusivamente de cuatro en cuatro en los vértices bajo el ángulo tetraédrico exacto de $\arccos(-1/3) \approx 109,4712^\circ$.

¿Por qué una pompa de jabón libre es siempre esférica? ¿Por qué al conectar una pompa pequeña con una pompa grande mediante un tubo es la pompa pequeña la que se encoge e infla a la grande, y por qué esta misma inestabilidad no hace colapsar los 300 millones de alvéolos de nuestros pulmones en cada espiración? ¿Por qué las burbujas de agua pura o alcohol puro estallan en milisegundos mientras que una traza de surfactante otorga a la película una elasticidad autorreparadora? ¿Y cómo fue derrotada en 1994 la conjetura centenaria de Lord Kelvin sobre la partición óptima del espacio tridimensional gracias a la estructura de Weaire-Phelan?

En este artículo integramos la termodinámica, la geometría diferencial, el cálculo de variaciones, la teoría geométrica de la medida y la mecánica de fluidos interfacial para desvelar toda la ciencia de las pompas y las espumas.

1. Termodinámica de la tensión superficial y asimetría molecular

En el interior (seno o bulk) del agua líquida, cada molécula de $\text{H}_2\text{O}$ está rodeada isotrópicamente por moléculas vecinas unidas mediante puentes de hidrógeno y fuerzas de van der Waals. La fuerza neta promediada en el tiempo es nula y cada molécula se halla en un estado de baja energía potencial (número de coordinación $z_b \approx 12$, energía cohesiva $-u_0$).

En la interfaz líquido-vapor, esta simetría se rompe drásticamente: la densidad numérica de moléculas en la fase gaseosa es unas mil veces menor que en el líquido. Una molécula superficial pierde cerca de la mitad de sus vecinos atractivos ($z_s < z_b$) y posee un exceso de energía potencial. Llevar moléculas desde el interior hasta la superficie para incrementar el área interfacial en $dA$ exige realizar un trabajo reversible positivo contra las fuerzas cohesivas.

El diferencial total de la energía libre de Gibbs $G(T, P, A, \{N_i\})$ es:

de donde se obtiene la definición termodinámica exacta de la tensión superficial $\gamma$:

Esta magnitud posee dos interpretaciones físicamente equivalentes:

1. Exceso de energía libre por unidad de área ($\text{J/m}^2$): A volumen $V$ constante, el sistema minimiza su energía libre $\gamma A$ minimizando su superficie $A$. Por la desigualdad isoperimétrica tridimensional $A^3 \ge 36\pi V^2$, la esfera es el único cuerpo que minimiza el área superficial para un volumen encerrado dado.
2. Fuerza de tracción tangencial por unidad de longitud ($\text{N/m}$): Según la fórmula de Kirkwood-Buff (1949) de la mecánica estadística, $\gamma$ es la integral de la anisotropía entre la presión normal $P_N$ y la presión tangencial $P_T(z)$ a través de la zona interfacial: $$\gamma = \int_{-\infty}^{+\infty} \bigl[ P_N - P_T(z) \bigr]\,dz$$

2. Geometría diferencial de superficies curvas y ecuación de Young-Laplace

Cuando una interfaz es curva, las fuerzas tangenciales $\gamma$ generan una componente normal neta hacia el centro de curvatura, provocando un salto de presión $\Delta P = P_{\text{int}} - P_{\text{ext}} > 0$. Sean $R_1, R_2$ los dos radios de curvatura principales ($\kappa_1 = 1/R_1, \kappa_2 = 1/R_2$) de una superficie suave $\Sigma \subset \mathbb{R}^3$. La curvatura media $H$ y la curvatura gaussiana $K$ se definen como:

Al aplicar un desplazamiento normal virtual $\delta\xi$ a un parche infinitesimal $dA = R_1 R_2\,d\theta_1\,d\theta_2$, la variación de área de primer orden es $\delta(dA) = dA\left(\frac{1}{R_1} + \frac{1}{R_2}\right)\delta\xi$ y el volumen barrido es $\delta(dV) = dA\,\delta\xi$. La condición de equilibrio $\delta F = \gamma\,\delta(dA) - \Delta P\,\delta(dV) = 0$ proporciona la ecuación de Young-Laplace (1805–1806):

Mientras que una gota de agua en el aire tiene una sola interfaz ($\Delta P_{\text{gota}} = 2\gamma/R$), una pompa de jabón en el aire es una capa líquida delgada de espesor $h \ll R$ que posee dos interfaces gas-líquido (interna en $R$ y externa en $R+h$), lo que introduce un factor $4$:

3. Superficies mínimas ($H = 0$), catástrofe del catenoide y problema de Plateau

Una película de jabón abierta tendida sobre un bastidor de alambre sin encerrar aire tiene la misma presión atmosférica en ambas caras ($\Delta P = 0$). La ecuación de Young-Laplace exige entonces que su curvatura media sea idénticamente nula en todo punto:

Las superficies con $H \equiv 0$ se denominan superficies mínimas. Para una gráfica $z = u(x,y)$, minimizar el funcional de área $A[u] = \iint \sqrt{1+u_x^2+u_y^2}\,dx\,dy$ conduce a la ecuación en derivadas parciales de Euler-Lagrange deducida por Lagrange en 1760:

Todo punto no plano de una película abierta es localmente un punto de silla hiperbólico ($K < 0$). La superficie mínima de revolución tendida entre dos anillos coaxiales de radio $R$ separados por una distancia $2h$ es el catenoide (Euler, 1744):

Al separar los anillos más allá de la relación crítica ($\coth u_c = u_c \approx 1,19968$):

la solución real desaparece mediante una bifurcación silla-nodo y el cuello de jabón se estrangula, saltando a dos discos planos separados (solución discontinua de Goldschmidt). El problema de Plateau sobre la existencia general de superficies mínimas fue resuelto en 1930 por Jesse Douglas (primera Medalla Fields en 1936) y Tibor Radó.

4. La paradoja de las dos pompas, maduración de Ostwald y el surfactante pulmonar (DPPC)

Si conectamos mediante un tubo con válvula una pompa pequeña ($R_{\text{pequeña}}$) y una pompa grande ($R_{\text{grande}}$), al abrir la válvula la pompa pequeña se encoge rápidamente e infla aún más a la pompa grande. Según $\Delta P = 4\gamma/R$, cuanto menor es el radio, mayor es la presión interna:

  flowchart TD
    A["Dos pompas de jabón conectadas por un tubo con válvula:<br/>Pompa pequeña (R_pequeña) y Pompa grande (R_grande)"] --> B["Ecuación de Young-Laplace: ΔP = 4γ / R<br/>La sobrepresión es inversamente proporcional al radio R"]
    B --> C["Presión pompa pequeña: ΔP_pequeña = 4γ / R_pequeña<br/>Presión pompa grande: ΔP_grande = 4γ / R_grande"]
    C --> D["Desequilibrio: ΔP_pequeña > ΔP_grande<br/>(¡Mayor curvatura genera mayor presión interna!)"]
    D --> E["Válvula abierta: el gas fluye de la pompa pequeña<br/>hacia la pompa grande"]
    E --> F["Retroalimentación positiva:<br/>La pompa pequeña se reduce más, su presión sube y colapsa"]
    D --> G["Crisis fisiológica en los pulmones:<br/>300 millones de alvéolos de radios distintos (50–200 μm)"]
    G --> H["Los neumocitos tipo II secretan surfactante pulmonar (DPPC)"]
    H --> I["Al contraerse el alvéolo, la monocapa de DPPC se comprime,<br/>reduciendo γ drásticamente hasta 0,5–2 mN/m"]
    I --> J["Se igualan las presiones (2γ_pequeña / R_pequeña ≈ 2γ_grande / R_grande)<br/>¡impidiendo el colapso alveolar (atelectasia)!"]

En una espuma real, la difusión de gas a través de las películas desde las burbujas pequeñas (alta presión) hacia las grandes (baja presión) produce el engrosamiento o maduración de Ostwald. En 1952, John von Neumann demostró mediante el teorema de Gauss-Bonnet que en una espuma seca 2D el área $A_n$ de una celda de $n$ lados evoluciona según la ley de von Neumann:

Las celdas con $n < 6$ lados se encogen y desaparecen, los hexágonos ($n = 6$) permanecen estacionarios y las celdas con $n > 6$ crecen (generalizado a 3D por MacPherson y Srolovitz en 2007).

En los pulmones humanos, unos 300 millones de alvéolos ($R \approx 50\text{–}200\,\mu\text{m}$) están interconectados. Si la tensión superficial $\gamma$ fuera constante, los alvéolos pequeños se vaciarían en los grandes y colapsarían (atelectasia). Los neumocitos tipo II evitan esto secretando surfactante pulmonar rico en dipalmitoilfosfatidilcolina (DPPC): al contraerse el alvéolo durante la espiración, la monocapa insoluble de DPPC se comprime tanto que $\gamma(R)$ cae hasta $0,5\text{–}2\,\text{mN/m}$, estabilizando todos los alvéolos.

5. Leyes de Joseph Plateau (1873), prueba de Jean Taylor (1976) y teorema de la doble burbuja

Formuladas en 1873 por el físico belga ciego Joseph Plateau y demostradas rigurosamente en 1976 por Jean E. Taylor mediante teoría geométrica de la medida:

1. 1.ª ley de Plateau (3 películas a $120^\circ$): Las películas de jabón se unen exclusivamente de tres en tres a lo largo de los bordes de Plateau formando ángulos iguales de $120^\circ$ ($2\pi/3$) ($\sum_{i=1}^3 \hat{\mathbf{t}}_i = \mathbf{0} \implies \cos\theta = -1/2$). Un cruce cuádruple a $90^\circ$ ($L_{\text{cruz}} = 4\sqrt{2} \approx 5,657$) es inestable y se divide instantáneamente en dos uniones triples a $120^\circ$ formando un árbol de Steiner ($L_{\text{Steiner}} = 2(1+\sqrt{3}) \approx 5,464$, transición T1).
2. 2.ª ley de Plateau (4 bordes a $109,47^\circ$): Los bordes de Plateau se cruzan exclusivamente de cuatro en cuatro en los vértices bajo el ángulo tetraédrico de Maraldi: $$\theta_{\text{Maraldi}} = \arccos\!\left(-\frac{1}{3}\right) \approx 109,47122^\circ \quad \left(109^\circ 28' 16''\right)$$
3. Teorema de la doble burbuja (Hutchings, Morgan, Ritoré, Ros 2000): Cuando dos pompas de radios $R_1 > R_2$ se unen compartiendo un tabique (septum), este se curva hacia el interior de la pompa mayor con radio de curvatura $R_{\text{septum}}$: $$\frac{1}{R_{\text{septum}}} = \frac{1}{R_2} - \frac{1}{R_1} \quad \Longleftrightarrow \quad R_{\text{septum}} = \frac{R_1 R_2}{R_1 - R_2}$$

6. La conjetura de Kelvin (1887) frente a la estructura de Weaire-Phelan (1994)

En 1887, Lord Kelvin planteó cómo dividir el espacio $\mathbb{R}^3$ en celdas de igual volumen ($V=1$) con área interfacial mínima. Debido a la frustración geométrica, Kelvin propuso el octaedro truncado (tetracaidecaedro de 6 cuadrados y 8 hexágonos, celda de Voronoi BCC) ligeramente curvado, con área normalizada $A_{\text{Kelvin}} \approx 5,30628$.

En 1994, Denis Weaire y Robert Phelan (Trinity College Dublin) refutaron la conjetura de Kelvin tras 107 años con la estructura de Weaire-Phelan (inspirada en los clatratos hidratos tipo I), cuya celda unidad cúbica contiene 8 poliedros de igual volumen de dos tipos: 2 dodecaedros pentagonales irregulares ($5^{12}$) y 6 tetracaidecaedros ($5^{12}6^2$ con 2 hexágonos y 12 pentágonos). Con un $88,9\%$ de caras pentagonales (cuyo ángulo interno de $108^\circ$ casi coincide con $109,47^\circ$), el área se reduce un $0,34\%$ hasta $A_{\text{WP}} \approx 5,28825$, estructura inmortalizada en el Centro Acuático Nacional de Pekín («Cubo de Agua»).

7. Elasticidad de Gibbs-Marangoni, interferencia de película delgada y presión de disyunción DLVO

Una baja tensión superficial por sí sola no estabiliza una burbuja (el etanol puro estalla al instante). La estabilidad de las películas de jabón reside en la elasticidad de Gibbs-Marangoni:

  flowchart TD
    A["Perturbación externa (aire, drenaje, onda térmica)<br/>estira y adelgaza localmente la película de jabón"] --> B["El área local ΔA se expande;<br/>la concentración superficial de surfactante Γ cae localmente"]
    B --> C["¡La menor densidad Γ hace que la tensión superficial local<br/>γ_delgada suba por encima de γ_gruesa circundante!"]
    C --> D["Aparece un gradiente tangencial ∇_s γ<br/>(esfuerzo cortante de Marangoni τ_M = ∇_s γ hacia la zona delgada)"]
    D --> E["Convección de Marangoni:<br/>La zona de alto γ tira de la monocapa y arrastra agua subyacente por viscosidad"]
    E --> F["¡El espesor h se autorrepara en milisegundos,<br/>evitando la ruptura capilar!"]
    F --> G["El drenaje gravitatorio adelgaza la película a longitudes de onda visibles (100–700 nm)"]
    G --> H["La interferencia óptica produce bandas iridiscentes de arcoíris"]
    H --> I["Para h < 30 nm, la presión de disyunción DLVO Π(h) frena a van der Waals;<br/>el desfase π forma la 'película negra de Newton' (h ≈ 4 nm)"]

Cuando una película se estira y adelgaza localmente, la concentración superficial de surfactante $\Gamma$ disminuye y la tensión superficial local $\gamma_{\text{delgada}}$ aumenta. El esfuerzo de Marangoni $\boldsymbol{\tau}_M = \nabla_s \gamma$ y el módulo de elasticidad de Gibbs

arrastran inmediatamente líquido hacia la zona adelgazada para repararla (el mismo efecto produce las «lágrimas del vino»). Al reducirse el espesor $h$ al rango visible, las reflexiones en ambas caras generan interferencia de película delgada ($2n_{\text{película}}h\cos\theta_t + \lambda/2 = m\lambda$). Por debajo de $h \ll \lambda/4 \approx 30\,\text{nm}$, el cambio de fase de $\pi$ en la cara frontal produce interferencia destructiva en todo el espectro visible y la película se vuelve negra (película negra de Newton, $h \approx 4\,\text{nm}$), sostenida por la presión de disyunción DLVO $\Pi(h) = \Pi_{\text{vdW}}(h) + \Pi_{\text{el}}(h) + \Pi_{\text{sr}}(h)$.