Ao mergulhar um aro em água com sabão e soprar suavemente, uma fina película líquida estica-se, destaca-se e relaxa em milissegundos numa esfera iridescente perfeita flutuando no ar. Se observarmos a espuma de sabão com uma lupa, veremos que o aparente caos obedece a uma ordem geométrica rigorosa: as películas líquidas nunca se encontram de quatro ou cinco, mas exclusivamente de três em três ao longo de canais líquidos formando ângulos exatos de $120^\circ$ ($2\pi/3$), e essas arestas cruzam-se exclusivamente de quatro em quatro nos vértices segundo o ângulo tetraédrico exato de $\arccos(-1/3) \approx 109,4712^\circ$.

Por que uma bolha de sabão livre é sempre esférica? Por que, ao conectar uma bolha pequena a uma bolha grande por um tubo, é a bolha pequena que encolhe e infla a grande — e como os 300 milhões de alvéolos dos nossos pulmões escapam desse colapso a cada expiração? Por que bolhas de água pura ou álcool puro estouram em milissegundos, enquanto um traço de surfactante confere à película uma elasticidade autorregenerativa? E como a conjectura centenária de Lord Kelvin sobre a partição do espaço tridimensional em células de mesmo volume e área mínima foi superada em 1994 pela estrutura de Weaire-Phelan?

Neste artigo, unimos termodinâmica, geometria diferencial, cálculo variacional, teoria geométrica da medida e mecânica dos fluidos interfacial para desvendar a física das bolhas e espumas.

1. Termodinâmica da tensão superficial e assimetria molecular

No interior (bulk) da água líquida, cada molécula de $\text{H}_2\text{O}$ é cercada isotropicamente por moléculas vizinhas ligadas por pontes de hidrogênio e forças de van der Waals. A força resultante média é nula e a molécula encontra-se num estado de baixa energia potencial (número de coordenação $z_b \approx 12$, energia coesiva $-u_0$).

Na interface líquido-vapor, essa simetria é quebrada: a densidade molecular na fase gasosa é cerca de mil vezes menor que no líquido. Uma molécula na superfície perde quase metade das suas vizinhas atrativas ($z_s < z_b$) e possui um excesso de energia potencial. Trazer moléculas do interior para a superfície para aumentar a área interfacial em $dA$ exige realizar trabalho reversível positivo contra as forças de coesão.

A diferencial total da energia livre de Gibbs $G(T, P, A, \{N_i\})$ é:

de onde decorre a definição termodinâmica exata da tensão superficial $\gamma$:

Ela possui duas interpretações fisicamente equivalentes:

1. Excesso de energia livre por unidade de área ($\text{J/m}^2$): A volume $V$ constante, um sistema minimiza sua energia livre $\gamma A$ minimizando sua área superficial $A$. Pela desigualdade isoperimétrica tridimensional $A^3 \ge 36\pi V^2$, a esfera é o único sólido que minimiza a área para um dado volume.
2. Força de tração tangencial por unidade de comprimento ($\text{N/m}$): Pela fórmula de Kirkwood-Buff (1949) da mecânica estatística, $\gamma$ é a integral da anisotropia entre a pressão normal $P_N$ e a pressão tangencial $P_T(z)$ através da camada interfacial: $$\gamma = \int_{-\infty}^{+\infty} \bigl[ P_N - P_T(z) \bigr]\,dz$$

2. Geometria diferencial de superfícies curvas e equação de Young-Laplace

Quando uma interface é curva, a tensão tangencial $\gamma$ gera uma força normal líquida apontando para o lado côncavo, exigindo um salto de pressão $\Delta P = P_{\text{int}} - P_{\text{ext}} > 0$. Sejam $R_1, R_2$ os dois raios de curvatura principais ($\kappa_1 = 1/R_1, \kappa_2 = 1/R_2$) de uma superfície suave $\Sigma \subset \mathbb{R}^3$. A curvatura média $H$ e a curvatura gaussiana $K$ são definidas por:

Um deslocamento normal virtual $\delta\xi$ de um elemento de área $dA = R_1 R_2\,d\theta_1\,d\theta_2$ produz uma variação de área $\delta(dA) = dA\left(\frac{1}{R_1} + \frac{1}{R_2}\right)\delta\xi$ e um volume varrido $\delta(dV) = dA\,\delta\xi$. O equilíbrio termodinâmico $\delta F = \gamma\,\delta(dA) - \Delta P\,\delta(dV) = 0$ fornece a equação de Young-Laplace (1805–1806):

Enquanto uma gota de água possui apenas uma interface ($\Delta P_{\text{gota}} = 2\gamma/R$), uma bolha de sabão no ar é uma fina casca líquida de espessura $h \ll R$ dotada de duas interfaces gás-líquido (interna de raio $R$ e externa de raio $R+h$), resultando no fator $4$:

3. Superfícies mínimas ($H = 0$), catástrofe do catenoide e problema de Plateau

Uma película de sabão aberta esticada num arame sem aprisionar ar está sujeita à mesma pressão atmosférica em ambos os lados ($\Delta P = 0$). A equação de Young-Laplace impõe que sua curvatura média seja identicamente nula em todos os pontos:

Superfícies com $H \equiv 0$ são chamadas de superfícies mínimas. Na forma $z = u(x,y)$, minimizar o funcional de área $A[u] = \iint \sqrt{1+u_x^2+u_y^2}\,dx\,dy$ leva à equação diferencial parcial de Euler-Lagrange obtida por Lagrange em 1760:

Todo ponto não plano de uma película aberta é localmente um ponto de sela hiperbólico ($K < 0$). A superfície mínima de revolução entre dois anéis coaxiais de raio $R$ separados por $2h$ é o catenoide (Euler, 1744):

Ao afastar os anéis além da razão crítica ($\coth u_c = u_c \approx 1,19968$):

a solução real desaparece numa bifurcação sela-nó e o pescoço da película rompe-se em dois discos planos (solução descontínua de Goldschmidt). O problema de Plateau foi resolvido em 1930 por Jesse Douglas (primeira Medalha Fields em 1936) e Tibor Radó.

4. O paradoxo das duas bolhas, amadurecimento de Ostwald e surfactante pulmonar (DPPC)

Conectando uma bolha pequena ($R_{\text{pequena}}$) e uma bolha grande ($R_{\text{grande}}$) por um tubo com válvula, ao abrir a válvula a bolha pequena encolhe rapidamente e infla ainda mais a bolha grande! Pela lei $\Delta P = 4\gamma/R$, quanto menor o raio, maior a pressão interna:

  flowchart TD
    A["Duas bolhas de sabão conectadas por um tubo com válvula:<br/>Bolha pequena (R_pequena) e Bolha grande (R_grande)"] --> B["Equação de Young-Laplace: ΔP = 4γ / R<br/>A sobrepressão é inversamente proporcional ao raio R"]
    B --> C["Pressão bolha pequena: ΔP_pequena = 4γ / R_pequena<br/>Pressão bolha grande: ΔP_grande = 4γ / R_grande"]
    C --> D["Desequilíbrio: ΔP_pequena > ΔP_grande<br/>(Maior curvatura gera maior pressão interna!)"]
    D --> E["Válvula aberta: o gás flui da bolha pequena de alta pressão<br/>para a bolha grande de baixa pressão"]
    E --> F["Retroalimentação positiva:<br/>A bolha pequena encolhe mais, sua pressão dispara e ela colapsa"]
    D --> G["Crise fisiológica nos pulmões:<br/>300 milhões de alvéolos de raios variados (50–200 μm)"]
    G --> H["Pneumócitos tipo II secretam surfactante pulmonar (DPPC)"]
    H --> I["Ao contrair o alvéolo, a monocamada de DPPC comprime-se,<br/>reduzindo γ drasticamente para 0,5–2 mN/m"]
    I --> J["As pressões de Laplace igualam-se (2γ_pequena / R_pequena ≈ 2γ_grande / R_grande)<br/>impedindo o colapso alveolar (atelectasia)!"]

Em espumas reais, a difusão de gás das bolhas pequenas para as grandes através das películas causa o amadurecimento de Ostwald (coarsening). Em 1952, John von Neumann provou via teorema de Gauss-Bonnet que, numa espuma seca 2D, a área $A_n$ de uma bolha de $n$ lados evolui segundo a lei de von Neumann:

Bolhas com $n < 6$ lados encolhem e somem, hexágonos ($n = 6$) permanecem estacionários e bolhas com $n > 6$ crescem (generalizado para 3D por MacPherson e Srolovitz em 2007).

Nos pulmões humanos, cerca de 300 milhões de alvéolos ($R \approx 50\text{–}200\,\mu\text{m}$) comunicam-se entre si. Se a tensão superficial $\gamma$ fosse constante, os alvéolos menores esvaziariam nos maiores e colapsariam (atelectasia). Os pneumócitos tipo II impedem isso secretando o surfactante pulmonar rico em dipalmitoilfosfatidilcolina (DPPC): na expiração, a compressão da monocamada insolúvel de DPPC derruba $\gamma(R)$ para $0,5\text{–}2\,\text{mN/m}$, estabilizando todos os alvéolos.

5. Leis de Joseph Plateau (1873), prova de Jean Taylor (1976) e teorema da bolha dupla

Formuladas em 1873 pelo físico belga Joseph Plateau e provadas rigorosamente em 1976 por Jean E. Taylor via teoria geométrica da medida:

1. 1.ª lei de Plateau (3 películas a $120^\circ$): As películas de sabão encontram-se exclusivamente de três em três ao longo das bordas de Plateau em ângulos iguais de $120^\circ$ ($2\pi/3$) ($\sum_{i=1}^3 \hat{\mathbf{t}}_i = \mathbf{0} \implies \cos\theta = -1/2$). Uma junção quádrupla a $90^\circ$ ($L_{\text{cruz}} = 4\sqrt{2} \approx 5,657$) é instável e divide-se instantaneamente em duas junções triplas a $120^\circ$ numa árvore de Steiner ($L_{\text{Steiner}} = 2(1+\sqrt{3}) \approx 5,464$, transição T1).
2. 2.ª lei de Plateau (4 bordas a $109,47^\circ$): As bordas de Plateau encontram-se exclusivamente de quatro em quatro nos vértices segundo o â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 da bolha dupla (Hutchings, Morgan, Ritoré, Ros 2000): Quando duas bolhas de raios $R_1 > R_2$ coalescem partilhando um septo, este curva-se para dentro da bolha maior com raio de curvatura $R_{\text{septo}}$: $$\frac{1}{R_{\text{septo}}} = \frac{1}{R_2} - \frac{1}{R_1} \quad \Longleftrightarrow \quad R_{\text{septo}} = \frac{R_1 R_2}{R_1 - R_2}$$

6. A conjectura de Kelvin (1887) contra a estrutura de Weaire-Phelan (1994)

Em 1887, Lord Kelvin perguntou como particionar o espaço $\mathbb{R}^3$ em células de igual volume ($V=1$) minimizando a área interfacial. Devido à frustração geométrica, Kelvin propôs o octaedro truncado (tetracaidecaedro de 6 quadrados e 8 hexágonos, célula de Voronoi BCC) levemente curvado, com área normalizada $A_{\text{Kelvin}} \approx 5,30628$.

Em 1994, Denis Weaire e Robert Phelan (Trinity College Dublin) derrubaram a conjectura de Kelvin após 107 anos com a estrutura de Weaire-Phelan (inspirada nos clatratos hidratos tipo I), cuja célula unitária cúbica reúne 8 poliedros de mesmo volume de dois tipos: 2 dodecaedros pentagonais irregulares ($5^{12}$) e 6 tetracaidecaedros ($5^{12}6^2$ com 2 hexágonos e 12 pentágonos). Com $88,9\%$ de faces pentagonais (ângulo interno de $108^\circ$ próximo a $109,47^\circ$), a área cai $0,34\%$ para $A_{\text{WP}} \approx 5,28825$, estrutura consagrada no Centro Aquático Nacional de Pequim («Cubo d’Água»).

7. Elasticidade de Gibbs-Marangoni, interferência em película fina e pressão de disjunção DLVO

Uma baixa tensão superficial isolada não estabiliza uma bolha (o etanol puro estoura instantaneamente). A estabilidade das películas de sabão decorre da elasticidade de Gibbs-Marangoni:

  flowchart TD
    A["Perturbação externa (corrente de ar, drenagem, onda térmica)<br/>estica e afina localmente a película de sabão"] --> B["A área local ΔA expande-se;<br/>a concentração superficial de surfactante Γ cai localmente"]
    B --> C["A menor densidade Γ faz a tensão superficial local<br/>γ_fina subir acima da γ_espessa vizinha!"]
    C --> D["Surge um gradiente tangencial ∇_s γ<br/>(tensão de cisalhamento de Marangoni τ_M = ∇_s γ rumo à região fina)"]
    D --> E["Convecção de Marangoni:<br/>A zona de alto γ puxa a monocamada e arrasta água subjacente por viscosidade"]
    E --> F["A espessura h autorregenera-se em milissegundos,<br/>impedindo a ruptura capilar!"]
    F --> G["A drenagem gravitacional afina a película até comprimentos de onda visíveis (100–700 nm)"]
    G --> H["A interferência óptica produz faixas iridescentes de arco-íris"]
    H --> I["Para h < 30 nm, a pressão de disjunção DLVO Π(h) equilibra van der Waals;<br/>a inversão de fase π forma o 'filme negro de Newton' (h ≈ 4 nm)"]

Quando uma película afina localmente, a concentração superficial de surfactante $\Gamma$ diminui e a tensão superficial local $\gamma_{\text{fina}}$ aumenta. A tensão de Marangoni $\boldsymbol{\tau}_M = \nabla_s \gamma$ e o módulo de elasticidade de Gibbs

puxam líquido de volta para a região afinada (o mesmo efeito produz as «lágrimas do vinho»). Quando a espessura $h$ atinge o espectro visível, ocorre interferência em película fina ($2n_{\text{película}}h\cos\theta_t + \lambda/2 = m\lambda$). Abaixo de $h \ll \lambda/4 \approx 30\,\text{nm}$, o salto de fase de $\pi$ na reflexão frontal gera interferência destrutiva em todo o espectro visível e a película torna-se negra (filme negro de Newton, $h \approx 4\,\text{nm}$), sustentada pela pressão de disjunção DLVO $\Pi(h) = \Pi_{\text{vdW}}(h) + \Pi_{\text{el}}(h) + \Pi_{\text{sr}}(h)$.