Celupkan cincin plastik ke dalam air sabun lalu tiup perlahan: selaput tipis cairan akan merentang, terlepas, dan dalam hitungan milidetik membentuk bola sempurna berwarna-warni pelangi yang melayang di udara. Atau amati busa sabun di dalam mangkuk menggunakan kaca pembesar: apa yang sekilas tampak seperti kumpulan gelembung acak ternyata tunduk pada keteraturan geometri yang mutlak. Di seluruh busa, dinding-dinding cairan tidak pernah bertemu empat atau lima sekaligus—mereka selalu bertemu bertiga di sepanjang saluran cairan (batas Plateau) dengan sudut tepat $120^\circ$ ($2\pi/3$), dan saluran-saluran tersebut selalu bertemu berempat di titik simpul dengan sudut tetrahedral tepat $\arccos(-1/3) \approx 109,4712^\circ$.

Mengapa gelembung sabun bebas selalu berbentuk bola? Mengapa ketika gelembung kecil dan gelembung besar dihubungkan dengan pipa, justru gelembung kecil yang menyusut dan memompa udaranya ke dalam gelembung besar—dan mengapa ketidakstabilan yang sama tidak membuat 300 juta alveolus di paru-paru kita kolaps setiap kali kita mengembuskan napas? Mengapa gelembung air murni atau alkohol murni pecah dalam milidetik, sedangkan sedikit surfaktan memberi lapisan film kemampuan menyembuhkan diri sendiri? Dan bagaimana konjektur Lord Kelvin yang bertahan selama 107 tahun tentang pembagian ruang tiga dimensi dengan luas permukaan minimum akhirnya dipatahkan pada tahun 1994 oleh struktur Weaire-Phelan?

Artikel ini memadukan termodinamika, geometri diferensial, kalkulus variasi, teori ukuran geometris, dan mekanika fluida antarmuka untuk mengungkap fisika di balik gelembung sabun dan busa.

1. Termodinamika Tegangan Permukaan dan Asimetri Molekuler

Di bagian dalam (bulk) air cair, setiap molekul $\text{H}_2\text{O}$ dikelilingi secara isotropis oleh molekul tetangga melalui jaringan ikatan hidrogen dan gaya van der Waals. Gaya total rata-rata terhadap waktu bernilai nol dan molekul berada pada keadaan energi potensial kohesif yang rendah (bilangan koordinasi $z_b \approx 12$, energi $-u_0$).

Pada antarmuka cair-uap, simetri ini rusak secara drastis: kerapatan molekul di fase gas sekitar seribu kali lebih kecil daripada di fase cair. Molekul di permukaan kehilangan hampir separuh tetangga penariknya ($z_s < z_b$) sehingga memiliki kelebihan energi potensial. Memindahkan molekul dari bagian dalam ke permukaan untuk memperluas luas antarmuka $A$ sebesar $dA$ membutuhkan kerja reversibel positif melawan gaya kohesi.

Diferensial total energi bebas Gibbs $G(T, P, A, \{N_i\})$ adalah:

Dari sini diperoleh definisi termodinamika eksak untuk tegangan permukaan $\gamma$:

Besaran ini memiliki dua makna fisis yang setara:

1. Energi bebas berlebih per satuan luas ($\text{J/m}^2$): Pada volume $V$ tetap, sistem meminimalkan energi bebas $\gamma A$ dengan meminimalkan luas permukaannya $A$. Berdasarkan ketaksamaan isoperimetrik tiga dimensi $A^3 \ge 36\pi V^2$, bola adalah satu-satunya bangun ruang yang memiliki luas permukaan minimum untuk volume tertentu.
2. Gaya tarik tangensial per satuan panjang ($\text{N/m}$): Menurut rumus Kirkwood-Buff (1949) dalam mekanika statistik, $\gamma$ adalah integral anisotropi antara tekanan normal $P_N$ dan tekanan tangensial $P_T(z)$ melintasi lapisan antarmuka: $$\gamma = \int_{-\infty}^{+\infty} \bigl[ P_N - P_T(z) \bigr]\,dz$$

2. Geometri Diferensial Permukaan Lengkung dan Persamaan Young-Laplace

Ketika suatu antarmuka melengkung, gaya tarik tangensial $\gamma$ menghasilkan resultan gaya normal menuju pusat kelengkungan (sisi cekung), sehingga menciptakan loncatan tekanan $\Delta P = P_{\text{dalam}} - P_{\text{luar}} > 0$. Misalkan $R_1, R_2$ adalah dua jari-jari kelengkungan utama ($\kappa_1 = 1/R_1, \kappa_2 = 1/R_2$) dari permukaan mulus $\Sigma \subset \mathbb{R}^3$. Kelengkungan rata-rata $H$ dan kelengkungan Gauss $K$ didefinisikan sebagai:

Pergeseran normal virtual $\delta\xi$ pada elemen luas $dA = R_1 R_2\,d\theta_1\,d\theta_2$ menghasilkan variasi luas orde pertama $\delta(dA) = dA\left(\frac{1}{R_1} + \frac{1}{R_2}\right)\delta\xi$ dan volume tersapu $\delta(dV) = dA\,\delta\xi$. Syarat kesetimbangan $\delta F = \gamma\,\delta(dA) - \Delta P\,\delta(dV) = 0$ menghasilkan persamaan Young-Laplace (1805–1806):

Sementara tetesan air di udara hanya memiliki satu antarmuka ($\Delta P_{\text{tetes}} = 2\gamma/R$), gelembung sabun di udara adalah cangkang cairan tipis dengan ketebalan $h \ll R$ yang memiliki dua antarmuka gas-cair (bagian dalam berjari-jari $R$ dan bagian luar berjari-jari $R+h$), sehingga menghasilkan faktor $4$:

3. Permukaan Minimal ($H = 0$), Katastrofe Katenoid, dan Masalah Plateau

Selaput sabun terbuka yang direntangkan pada kerangka kawat tanpa mengurung udara memiliki tekanan atmosfer yang sama di kedua sisinya ($\Delta P = 0$). Persamaan Young-Laplace mengharuskan kelengkungan rata-rata bernilai nol di setiap titik:

Permukaan dengan $H \equiv 0$ disebut permukaan minimal. Dalam bentuk grafik $z = u(x,y)$, meminimalkan fungsional luas $A[u] = \iint \sqrt{1+u_x^2+u_y^2}\,dx\,dy$ menghasilkan persamaan diferensial parsial Euler-Lagrange (Lagrange, 1760):

Setiap titik tak-datar pada selaput sabun terbuka secara lokal berbentuk pelana hiperbolik ($K < 0$). Permukaan minimal putar di antara dua cincin koaksial berjari-jari $R$ yang terpisah sejauh $2h$ adalah katenoid (Euler, 1744):

Saat kedua cincin ditarik menjauh melewati rasio kritis ($\coth u_c = u_c \approx 1,19968$):

solusi real menghilang melalui bifurkasi pelana-simpul; leher selaput sabun putus dan melompat menjadi dua cakram datar terpisah (solusi diskontinu Goldschmidt). Masalah Plateau untuk kurva batas umum diselesaikan pada tahun 1930 oleh Jesse Douglas (peraih Medali Fields pertama tahun 1936) dan Tibor Radó.

4. Paradoks Dua Gelembung, Pematangan Ostwald, dan Surfaktan Paru (DPPC)

Jika gelembung sabun kecil ($R_{\text{kecil}}$) dan gelembung besar ($R_{\text{besar}}$) dihubungkan melalui pipa berkatup, saat katup dibuka gelembung kecil akan menyusut cepat dan meniupkan udaranya ke gelembung besar! Berdasarkan $\Delta P = 4\gamma/R$, semakin kecil jari-jari gelembung, semakin tinggi tekanan di dalamnya:

  flowchart TD
    A["Dua Gelembung Sabun Dihubungkan Pipa Berkatup:<br/>Gelembung Kecil (R_kecil) dan Gelembung Besar (R_besar)"] --> B["Persamaan Young-Laplace: ΔP = 4γ / R<br/>Tekanan berlebih berbanding terbalik dengan jari-jari R"]
    B --> C["Tekanan gelembung kecil: ΔP_kecil = 4γ / R_kecil<br/>Tekanan gelembung besar: ΔP_besar = 4γ / R_besar"]
    C --> D["Ketidakseimbangan: ΔP_kecil > ΔP_besar<br/>(Kelengkungan tajam menghasilkan tekanan internal lebih tinggi!)"]
    D --> E["Katup Dibuka: Gas mengalir dari gelembung kecil bertekanan tinggi<br/>ke gelembung besar bertekanan rendah"]
    E --> F["Umpan Balik Positif:<br/>Gelembung kecil makin menyusut, tekanannya melonjak lalu kolaps"]
    D --> G["Krisis Fisiologis di Paru-Paru Mamalia:<br/>300 juta alveolus dengan jari-jari berbeda (50–200 μm)"]
    G --> H["Pneumosit Tipe II mensekresikan Surfaktan Paru (DPPC)"]
    H --> I["Saat alveolus menyusut, monolayer DPPC terkompresi rapat,<br/>menurunkan γ secara drastis hingga 0,5–2 mN/m"]
    I --> J["Tekanan Laplace menjadi setara (2γ_kecil / R_kecil ≈ 2γ_besar / R_besar)<br/>mencegah kolaps paru (atelektasis)!"]

Di dalam busa nyata, difusi gas menembus dinding cairan dari gelembung kecil ke gelembung besar menyebabkan pematangan Ostwald (coarsening). Pada tahun 1952, John von Neumann membuktikan lewat teorema Gauss-Bonnet bahwa dalam busa kering 2D, laju perubahan luas $A_n$ dari sel bersegi-$n$ mengikuti hukum von Neumann:

Sel dengan $n < 6$ sisi menyusut dan lenyap, segi enam ($n = 6$) tetap stabil stasioner, dan sel dengan $n > 6$ terus membesar (generalisasi 3D dibuktikan oleh MacPherson dan Srolovitz pada tahun 2007).

Di dalam paru-paru manusia terdapat sekitar 300 juta alveolus ($R \approx 50\text{–}200\,\mu\text{m}$). Jika tegangan permukaan $\gamma$ cairan pelapis alveolus bernilai konstan, alveolus kecil akan langsung mengempis ke dalam alveolus besar (atelektasis). Pneumosit tipe II mencegah hal ini dengan mensekresikan surfaktan paru yang kaya akan dipalmitoilfosfatidilkolin (DPPC): saat alveolus mengecil ketika mengembuskan napas, lapisan tunggal DPPC yang tak larut air terkompresi rapat sehingga menurunkan $\gamma(R)$ hingga mendekati nol ($0,5\text{–}2\,\text{mN/m}$), menyeimbangkan tekanan di seluruh alveolus.

5. Hukum Joseph Plateau (1873), Bukti Jean Taylor (1976), dan Teorema Gelembung Ganda

Dirumuskan secara eksperimental oleh fisikawan Belgia Joseph Plateau pada tahun 1873 dan dibuktikan secara matematis oleh Jean E. Taylor pada tahun 1976 melalui teori ukuran geometris:

1. Hukum Pertama Plateau (3 Film pada $120^\circ$): Selaput sabun selalu bertemu bertiga di sepanjang batas Plateau dengan sudut sama besar $120^\circ$ ($2\pi/3$) ($\sum_{i=1}^3 \hat{\mathbf{t}}_i = \mathbf{0} \implies \cos\theta = -1/2$). Persimpangan empat film pada $90^\circ$ ($L_{\text{silang}} = 4\sqrt{2} \approx 5,657$) tidak stabil dan langsung terbelah menjadi dua persimpangan tiga film $120^\circ$ membentuk pohon Steiner ($L_{\text{Steiner}} = 2(1+\sqrt{3}) \approx 5,464$, transisi T1).
2. Hukum Kedua Plateau (4 Tepi pada $109,47^\circ$): Batas-batas Plateau selalu bertemu berempat di titik sudut dengan sudut tetrahedral Maraldi: $$\theta_{\text{Maraldi}} = \arccos\!\left(-\frac{1}{3}\right) \approx 109,47122^\circ \quad \left(109^\circ 28' 16''\right)$$
3. Teorema Gelembung Ganda (Hutchings, Morgan, Ritoré, Ros 2000): Ketika dua gelembung berjari-jari $R_1 > R_2$ bergabung berbagi dinding sekat (septum), sekat tersebut melengkung ke dalam gelembung besar dengan jari-jari kelengkungan $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. Konjektur Kelvin (1887) vs. Struktur Weaire-Phelan (1994)

Pada tahun 1887, Lord Kelvin menanyakan cara membagi ruang $\mathbb{R}^3$ menjadi sel-sel bervolume sama ($V=1$) dengan luas permukaan minimum. Karena frustrasi geometris, Kelvin mengusulkan oktahedron terpotong (tetrakaidekahedron dengan 6 persegi dan 8 segi enam, sel Voronoi BCC) yang sedikit melengkung dengan luas ternormalisasi $A_{\text{Kelvin}} \approx 5,30628$.

Setelah bertahan 107 tahun, pada tahun 1994 Denis Weaire dan Robert Phelan (Trinity College Dublin) mematahkan konjektur Kelvin dengan struktur Weaire-Phelan (terinspirasi hidrat klatrat tipe I), yang sel satuan kubiknya terdiri atas 8 polihedron bervolume sama dari dua jenis: 2 dodekahedron pentagonal tak beraturan ($5^{12}$) dan 6 tetrakaidekahedron ($5^{12}6^2$ dengan 2 segi enam dan 12 segi lima). Karena $88,9\%$ sisinya adalah segi lima (sudut dalam $108^\circ$ sangat dekat dengan $109,47^\circ$), luas permukaannya turun $0,34\%$ menjadi $A_{\text{WP}} \approx 5,28825$, yang diabadikan pada arsitektur Pusat Akuatik Nasional Beijing (“Water Cube”).

7. Elastisitas Gibbs-Marangoni, Interferensi Film Tipis, dan Tekanan Disjoining DLVO

Tegangan permukaan yang rendah saja tidak mampu menstabilkan gelembung (etanol murni langsung pecah). Ketahanan selaput sabun berasal dari elastisitas Gibbs-Marangoni:

  flowchart TD
    A["Gangguan Eksternal (angin, drainase gravitasi, gelombang termal)<br/>merentangkan dan menipiskan selaput sabun secara lokal"] --> B["Luas antarmuka lokal ΔA membesar;<br/>konsentrasi permukaan surfaktan Γ turun secara lokal"]
    B --> C["Penurunan Γ membuat tegangan permukaan lokal<br/>γ_tipis naik melebihi γ_tebal di sekitarnya!"]
    C --> D["Muncul gradien tegangan permukaan tangensial ∇_s γ<br/>(tegangan geser Marangoni τ_M = ∇_s γ menuju titik tipis)"]
    D --> E["Konveksi Marangoni:<br/>Zona γ tinggi menarik monolayer surfaktan dan menyeret air di bawahnya"]
    E --> F["Ketebalan film h pulih sendiri dalam milidetik,<br/>mencegah pecahnya gelembung!"]
    F --> G["Drainase gravitasi menipiskan film ke panjang gelombang tampak (100–700 nm)"]
    G --> H["Interferensi film tipis menghasilkan pita warna pelangi"]
    H --> I["Pada h < 30 nm, tekanan disjoining DLVO Π(h) menahan gaya van der Waals;<br/>pembalikan fase π membentuk 'Film Hitam Newton' (h ≈ 4 nm)"]

Saat selaput sabun meregang dan menipis secara lokal, kerapatan surfaktan $\Gamma$ menurun sehingga tegangan permukaan lokal $\gamma_{\text{tipis}}$ meningkat. Tegangan Marangoni $\boldsymbol{\tau}_M = \nabla_s \gamma$ dan modulus elastisitas Gibbs

segera menarik cairan kembali ke bagian yang menipis (efek Marangoni yang sama menjelaskan fenomena “air mata anggur” pada dinding gelas). Saat ketebalan $h$ menyusut ke orde panjang gelombang cahaya tampak, pantulan dari permukaan depan dan belakang menghasilkan interferensi film tipis ($2n_{\text{film}}h\cos\theta_t + \lambda/2 = m\lambda$). Di bawah $h \ll \lambda/4 \approx 30\,\text{nm}$, pembalikan fase $\pi$ pada permukaan depan menyebabkan interferensi destruktif di seluruh spektrum tampak sehingga film tampak hitam pekat (film hitam Newton, $h \approx 4\,\text{nm}$), yang ditopang oleh tekanan disjoining DLVO $\Pi(h) = \Pi_{\text{vdW}}(h) + \Pi_{\text{el}}(h) + \Pi_{\text{sr}}(h)$.