Опустите колечко в мыльный раствор и слегка подуйте: тонкая жидкая пленка растянется, сомкнется и за считанные миллисекунды превратится в переливающуюся всеми цветами радуги идеальную сферу. А если взглянуть через лупу на густую мыльную пену (водную пену), то кажущийся хаос обнаружит безупречный геометрический порядок: жидкие перегородки никогда не сходятся по четыре или по пять — они встречаются исключительно по три вдоль жидких каналов (границ Плато) под точным углом $120^\circ$ ($2\pi/3$), а сами эти каналы сходятся в узлах исключительно по четыре под точным тетраэдрическим углом $\arccos(-1/3) \approx 109{,}4712^\circ$.
Почему свободный мыльный пузырь всегда принимает форму сферы? Почему при соединении маленького и большого мыльных пузырей трубкой именно маленький пузырь сдувается и накачивает большой — и почему эта же самая неустойчивость не заставляет 300 миллионов альвеол в наших легких схлопываться при каждом выдохе? Почему пузыри из чистой воды или чистого спирта лопаются за миллисекунды, тогда как следы поверхностно-активного вещества (ПАВ) наделяют пленку способностью к самозалечиванию? И как столетняя гипотеза лорда Кельвина о разбиении трехмерного пространства на ячейки равного объема с минимальной площадью поверхности была опровергнута в 1994 году структурой Уэйра — Фелана?
В этой статье мы объединим термодинамику, дифференциальную геометрию, вариационное исчисление, геометрическую теорию меры и гидродинамику межфазных границ, чтобы раскрыть математическую физику мыльных пузырей и пен.
1. Термодинамика поверхностного натяжения и молекулярная асимметрия
Внутри объема (bulk) жидкой воды каждая молекула $\text{H}_2\text{O}$ изотропно окружена соседними молекулами и связана с ними сеткой водородных связей и силами Ван-дер-Ваальса. Средняя по времени равнодействующая сила в объеме равна нулю, и молекула находится в глубокой потенциальной яме (координационное число $z_b \approx 12$, энергия когезии $-u_0$).
На границе раздела жидкость — пар эта симметрия резко нарушается: плотность молекул в газовой фазе примерно на три порядка ниже, чем в жидкости. Поверхностная молекула теряет около половины притягивающих соседей ($z_s < z_b$) и обладает избыточной потенциальной энергией. Чтобы перенести молекулы из объема на поверхность и увеличить площадь границы раздела $A$ на величину $dA$, необходимо совершить положительную обратимую работу против сил когезии.
Полный дифференциал энергии Гиббса $G(T, P, A, \{N_i\})$ имеет вид:
$$dG = -S\,dT + V\,dP + \gamma\,dA + \sum_{i} \mu_i\,dN_i$$Отсюда следует строгое термодинамическое определение поверхностного натяжения $\gamma$:
$$\gamma \equiv \left( \frac{\partial G}{\partial A} \right)_{T, P, N_i} = \left( \frac{\partial F}{\partial A} \right)_{T, V, N_i}$$Эта величина имеет две физически эквивалентные трактовки:
- Избыточная свободная энергия на единицу площади ($\text{Дж/м}^2$): При фиксированном объеме $V$ система стремится к минимуму свободной энергии $\gamma A$, минимизируя площадь поверхности $A$. Согласно трехмерному изопериметрическому неравенству $A^3 \ge 36\pi V^2$, единственным телом с минимальной площадью поверхности при заданном объеме является сфера.
- Касательная сила натяжения на единицу длины ($\text{Н/м}$): Согласно формуле Кирквуда — Баффа (1949) в статистической механике, $\gamma$ выражается через интеграл анизотропии тензора давления — разности нормального давления $P_N$ и тангенциального давления $P_T(z)$ поперек межфазного слоя: $$\gamma = \int_{-\infty}^{+\infty} \bigl[ P_N - P_T(z) \bigr]\,dz$$
| Жидкость / межфазная система (при $20^\circ\text{C}$ с воздухом) | Поверхностное натяжение $\gamma$ ($\text{мН/м}$) | Доминирующее взаимодействие | Устойчивость тонких пленок и пузырей |
|---|---|---|---|
| Жидкая ртуть ($\text{Hg}$) | $486{,}5$ | Металлическая связь | Пленки не образуются; мгновенно собирается в капли |
| Чистая вода ($\text{H}_2\text{O}$) | $72{,}8$ | Сетка водородных связей + диполи | Лопается за миллисекунды (нет упругости Марангони) |
| Чистый этанол ($\text{C}_2\text{H}_5\text{OH}$) | $22{,}3$ | Слабые Н-связи + дисперсионные силы | Однокомпонентная пленка неустойчива («слезы вина») |
| Мыльный раствор / раствор ДСН (вблизи ККМ) | $25{,}0 \sim 35{,}0$ | Монослой амфифильного ПАВ | Стабильна от минут до часов (эффект Гиббса — Марангони) |
| Легочный сурфактант (сжатый монослой ДПФХ) | $0{,}5 \sim 2{,}0$ | Плотно упакованные насыщенные фосфолипиды | Сверхстабильна (предотвращает спадение альвеол) |
2. Дифференциальная геометрия кривых поверхностей и уравнение Юнга — Лапласа
Когда граница раздела искривлена, касательные силы натяжения $\gamma$ создают результирующую нормальную силу, направленную в сторону вогнутости, что вызывает скачок давления $\Delta P = P_{\text{вн}} - P_{\text{нар}} > 0$. Пусть $R_1, R_2$ — два ортогональных главных радиуса кривизны ($\kappa_1 = 1/R_1, \kappa_2 = 1/R_2$ — главные кривизны) гладкой поверхности $\Sigma \subset \mathbb{R}^3$. Тогда средняя кривизна $H$ и гауссова кривизна $K$ определяются как:
$$H \equiv \frac{\kappa_1 + \kappa_2}{2} = \frac{1}{2}\left(\frac{1}{R_1} + \frac{1}{R_2}\right), \qquad K \equiv \kappa_1 \kappa_2 = \frac{1}{R_1 R_2}$$При бесконечно малом виртуальном смещении $\delta\xi$ вдоль нормали элементарного участка $dA = R_1 R_2\,d\theta_1\,d\theta_2$ приращение площади в первом порядке равно $\delta(dA) = dA\left(\frac{1}{R_1} + \frac{1}{R_2}\right)\delta\xi$, а заметаемый объем равен $\delta(dV) = dA\,\delta\xi$. Из условия равновесия $\delta F = \gamma\,\delta(dA) - \Delta P\,\delta(dV) = 0$ получаем уравнение Юнга — Лапласа (1805–1806):
$$\Delta P = P_{\text{вн}} - P_{\text{нар}} = \gamma\left(\frac{1}{R_1} + \frac{1}{R_2}\right) = 2\gamma H$$В то время как капля воды имеет только одну границу раздела ($\Delta P_{\text{капля}} = 2\gamma/R$), мыльный пузырь в воздухе представляет собой тонкую жидкую оболочку толщиной $h \ll R$ и имеет две границы раздела газ — жидкость (внутреннюю радиуса $R$ и внешнюю радиуса $R+h$), откуда возникает множитель $4$:
$$\Delta P_{\text{пузырь}} = (P_{\text{вн}} - P_{\text{пленка}}) + (P_{\text{пленка}} - P_{\text{нар}}) = \frac{2\gamma}{R} + \frac{2\gamma}{R+h} \approx \frac{4\gamma}{R}$$3. Минимальные поверхности ($H = 0$), катастрофа катеноида и проблема Плато
Открытая мыльная пленка, натянутая на проволочный контур без замкнутого объема воздуха, испытывает одинаковое атмосферное давление с обеих сторон ($\Delta P = 0$). Уравнение Юнга — Лапласа требует, чтобы ее средняя кривизна тождественно обращалась в нуль в каждой точке:
$$H = \frac{1}{2}\left(\frac{1}{R_1} + \frac{1}{R_2}\right) = 0 \quad \Longleftrightarrow \quad \kappa_1 = -\kappa_2 \quad \Longrightarrow \quad K = -\kappa_1^2 \le 0$$Поверхности с $H \equiv 0$ называются минимальными поверхностями. Для графика функции $z = u(x,y)$ минимизация функционала площади $A[u] = \iint \sqrt{1+u_x^2+u_y^2}\,dx\,dy$ дает нелинейное уравнение Эйлера — Лагранжа (Лагранж, 1760):
$$\bigl(1 + u_y^2\bigr) u_{xx} - 2 u_x u_y u_{xy} + \bigl(1 + u_x^2\bigr) u_{yy} = 0$$Любая неплоская точка открытой мыльной пленки локально является гиперболическим седлом ($K < 0$). Классическая минимальная поверхность вращения, натянутая между двумя соосными кольцами радиуса $R$ на расстоянии $2h$, — это открытый Эйлером в 1744 году катеноид:
$$r(z) = c \cosh\!\left(\frac{z}{c}\right)$$При раздвижении колец свыше критического отношения ($\coth u_c = u_c \approx 1{,}19968$):
$$\left(\frac{2h}{R}\right)_{\text{крит}} = \frac{2u_c}{\cosh u_c} \approx 1{,}32548$$вещественное решение исчезает в результате седло-узловой бифуркации, шейка катеноида стремительно схлопывается и пленка распадается на два плоских диска (разрывное решение Гольдшмидта). Общая проблема Плато была решена в 1930 году Джесси Дугласом (первая Филдсовская премия 1936 года) и Тибором Радо.
4. Парадокс двух пузырей, созревание Оствальда и легочный сурфактант (ДПФХ)
Если соединить трубкой с краном маленький мыльный пузырь ($R_{\text{мал}}$) и большой пузырь ($R_{\text{бол}}$), то при открытии крана маленький пузырь быстро сдуется, перекачав свой воздух в большой пузырь! Ведь согласно $\Delta P = 4\gamma/R$ давление внутри пузыря обратно пропорционально его радиусу:
$$R_{\text{мал}} < R_{\text{бол}} \quad \Longrightarrow \quad \Delta P_{\text{мал}} = \frac{4\gamma}{R_{\text{мал}}} \;>\; \frac{4\gamma}{R_{\text{бол}}} = \Delta P_{\text{бол}}$$
flowchart TD
A["Два мыльных пузыря на концах трубки с краном:<br/>Малый пузырь (R_мал) и Большой пузырь (R_бол)"] --> B["Уравнение Юнга — Лапласа: ΔP = 4γ / R<br/>Избыточное давление обратно пропорционально радиусу R"]
B --> C["Давление в малом: ΔP_мал = 4γ / R_мал<br/>Давление в большом: ΔP_бол = 4γ / R_бол"]
C --> D["Дисбаланс давлений: ΔP_мал > ΔP_бол<br/>(Большая кривизна создает большее давление внутри!)"]
D --> E["Открытие крана: воздух перетекает из малого пузыря<br/>в большой пузырь"]
E --> F["Положительная обратная связь:<br/>Малый пузырь сжимается, его давление растет, и он схлопывается"]
D --> G["Физиологическая угроза в легких человека:<br/>300 миллионов альвеол разного радиуса (50–200 мкм)"]
G --> H["Пневмоциты II типа секретируют легочный сурфактант (ДПФХ)"]
H --> I["При сжатии альвеолы монослой ДПФХ уплотняется,<br/>снижая натяжение γ почти до нуля (0,5–2 мН/м)"]
I --> J["Давления Лапласа выравниваются (2γ_мал / R_мал ≈ 2γ_бол / R_бол),<br/>предотвращая спадение легких (ателектаз)!"]
В реальной пене диффузия газа сквозь жидкие перегородки из маленьких пузырьков в большие вызывает укрупнение пены — оствальдовское созревание. В 1952 году Джон фон Нейман с помощью теоремы Гаусса — Бонне доказал, что в двумерной сухой пене площадь $A_n$ ячейки с $n$ сторонами меняется по закону фон Неймана:
$$\frac{dA_n}{dt} = K_0 (n - 6)$$Ячейки с числом сторон $n < 6$ сжимаются и исчезают, шестиугольники ($n = 6$) стабильны, а ячейки с $n > 6$ растут (трехмерное обобщение получили Макферсон и Сроловиц в 2007 году).
В легких человека около 300 миллионов альвеол ($R \approx 50\text{–}200\,\text{мкм}$) соединены бронхиолами. Если бы поверхностное натяжение $\gamma$ было постоянным, мелкие альвеолы мгновенно опорожнялись бы в крупные (ателектаз). Пневмоциты II типа предотвращают это, выделяя легочный сурфактант на основе дипальмитоилфосфатидилхолина (ДПФХ / DPPC): при выдохе нерастворимый монослой ДПФХ сжимается так плотно, что $\gamma(R)$ падает до $0{,}5\text{–}2\,\text{мН/м}$, уравнивая давление во всех альвеолах.
5. Законы Жозефа Плато (1873), доказательство Джин Тейлор (1976) и теорема о двойном пузыре
Сформулированные в 1873 году ослепшим бельгийским физиком Жозефом Плато и строго доказанные в 1976 году Джин Тейлор (Jean E. Taylor) методами геометрической теории меры:
- 1-й закон Плато (3 пленки под углом $120^\circ$): Мыльные пленки сходятся исключительно по три вдоль границ Плато под равными углами $120^\circ$ ($2\pi/3$) ($\sum_{i=1}^3 \hat{\mathbf{t}}_i = \mathbf{0} \implies \cos\theta = -1/2$). Четырехкратное пересечение под углом $90^\circ$ ($L_{\text{крест}} = 4\sqrt{2} \approx 5{,}657$) неустойчиво и мгновенно расщепляется на два тройных узла по $120^\circ$ дерева Штейнера ($L_{\text{Штейнер}} = 2(1+\sqrt{3}) \approx 5{,}464$, T1-переход).
- 2-й закон Плато (4 ребра под углом $109{,}47^\circ$): Границы Плато сходятся в вершинах исключительно по четыре под тетраэдрическим углом Маральди: $$\theta_{\text{Maraldi}} = \arccos\!\left(-\frac{1}{3}\right) \approx 109{,}47122^\circ \quad \left(109^\circ 28' 16''\right)$$
- Теорема о двойном пузыре (Hutchings, Morgan, Ritoré, Ros 2000): Когда два пузыря радиусов $R_1 > R_2$ сливаются, общая перегородка (септа) выгибается внутрь большего пузыря с радиусом кривизны $R_{\text{септа}}$: $$\frac{1}{R_{\text{септа}}} = \frac{1}{R_2} - \frac{1}{R_1} \quad \Longleftrightarrow \quad R_{\text{септа}} = \frac{R_1 R_2}{R_1 - R_2}$$
6. Гипотеза Кельвина (1887) против структуры Уэйра — Фелана (1994)
В 1887 году лорд Кельвин поставил задачу о разбиении пространства $\mathbb{R}^3$ на ячейки равного объема ($V=1$) с минимальной площадью границ. Из-за геометрической фрустрации Кельвин предложил слегка искривленный усеченный октаэдр (14-гранник из 6 квадратов и 8 шестиугольников, ячейка Вороного ОЦК-решетки) с нормированной площадью $A_{\text{Kelvin}} \approx 5{,}30628$.
Спустя 107 лет, в 1994 году, физики Денис Уэйр и Роберт Фелан (Тринити-колледж, Дублин) опровергли гипотезу Кельвина, открыв структуру Уэйра — Фелана (вдохновленную клатратными гидратами I типа), кубическая элементарная ячейка которой состоит из 8 равновеликих многогранников двух типов: 2 неправильных пентагондодекаэдров ($5^{12}$) и 6 тетракайдекаэдров ($5^{12}6^2$ с 2 шестиугольниками и 12 пятиугольниками). Поскольку $88{,}9\%$ граней являются пятиугольниками (внутренний угол $108^\circ$ близок к $109{,}47^\circ$), площадь сокращается на $0{,}34\%$ до $A_{\text{WP}} \approx 5{,}28825$, что легло в основу архитектуры Пекинского национального плавательного комплекса («Водяной куб»).
| 3D-структура пены, заполняющая пространство | Ячеек в элем. ячейке $Z$ | Топология граней (среднее $\langle F \rangle$) | Законы Плато? | Нормированная площадь $A$ ($V=1$) | Разница с Кельвином |
|---|---|---|---|---|---|
| Простая кубическая решетка | $1$ | $6$ квадратов ($\langle F \rangle = 6{,}00$) | Нет (неустойчива) | $6{,}00000$ | $+13{,}07\%$ |
| Ромбододекаэдр (ГЦК) | $1$ | $12$ ромбов ($\langle F \rangle = 12{,}00$) | Нет (6-лучевые вершины распадаются) | $5{,}34539$ | $+0{,}74\%$ |
| Структура Кельвина (1887, ОЦК) | $1$ | $6$ квадратов + $8$ шестиугольников ($\langle F \rangle = 14{,}00$) | Да (изогнутые шестиугольники) | $5{,}30628$ | $0{,}00\%$ (База) |
| Структура Уэйра — Фелана (1994, A15) | $8$ ($2+6$) | $88{,}9\%$ пятиугольников + $11{,}1\%$ шестиугольников ($\langle F \rangle = 13{,}50$) | Да (минимальная кривизна) | $5{,}28825$ | $-0{,}34\%$ (Рекорд) |
7. Упругость Гиббса — Марангони, тонкопленочная интерференция и расклинивающее давление ДЛФО
Низкое поверхностное натяжение само по себе не стабилизирует пузырь (чистый этанол лопается мгновенно). Секрет живучести мыльной пленки — упругость Гиббса — Марангони:
flowchart TD
A["Внешнее возмущение (поток воздуха, дренаж, тепловая волна)<br/>локально растягивает и утончает мыльную пленку"] --> B["Локальная площадь ΔA растет;<br/>поверхностная концентрация ПАВ Γ локально падает"]
B --> C["Снижение Γ вызывает рост локального поверхностного<br/>натяжения γ_тонк выше окружающего γ_толст!"]
C --> D["Возникает касательный градиент ∇_s γ<br/>(напряжение сдвига Марангони τ_M = ∇_s γ направлено к тонкому участку)"]
D --> E["Конвекция Марангони:<br/>Зона высокого γ втягивает монослой ПАВ и увлекает за собой объемную воду"]
E --> F["Толщина пленки h самовосстанавливается за миллисекунды,<br/>предотвращая разрыв!"]
F --> G["Гравитационный дренаж утончает пленку до длин волн видимого света (100–700 нм)"]
G --> H["Тонкопленочная интерференция рождает радужные полосы"]
H --> I["При h < 30 нм расклинивающее давление ДЛФО Π(h) уравновешивает Ван-дер-Ваальса;<br/>сдвиг фазы π создает 'ньютоновскую черную пленку' (h ≈ 4 нм)"]
При локальном растяжении пленки поверхностная концентрация ПАВ $\Gamma$ падает, а локальное натяжение $\gamma_{\text{тонк}}$ возрастает. Возникающее напряжение Марангони $\boldsymbol{\tau}_M = \nabla_s \gamma$ и модуль упругости Гиббса
$$E_G \equiv 2A\left(\frac{\partial\gamma}{\partial A}\right)_{N_s} = -2\Gamma\left(\frac{\partial\gamma}{\partial\Gamma}\right)_T > 0$$немедленно подтягивают жидкость обратно в утонченную зону (тот же эффект объясняет «слезы вина» на стенках бокала). Когда толщина $h$ уменьшается до оптических длин волн, возникает тонкопленочная интерференция ($2n_{\text{пленка}}h\cos\theta_t + \lambda/2 = m\lambda$). При $h \ll \lambda/4 \approx 30\,\text{нм}$ потеря полуволны ($\pi$) при отражении от передней границы вызывает деструктивную интерференцию для всего видимого спектра: пленка становится черной (ньютоновская черная пленка, $h \approx 4\,\text{нм}$), удерживаемая от разрыва расклинивающим давлением Дерягина (ДЛФО / DLVO) $\Pi(h) = \Pi_{\text{vdW}}(h) + \Pi_{\text{el}}(h) + \Pi_{\text{sr}}(h)$.