Опустите колечко в мыльный раствор и слегка подуйте: тонкая жидкая пленка растянется, сомкнется и за считанные миллисекунды превратится в переливающуюся всеми цветами радуги идеальную сферу. А если взглянуть через лупу на густую мыльную пену (водную пену), то кажущийся хаос обнаружит безупречный геометрический порядок: жидкие перегородки никогда не сходятся по четыре или по пять — они встречаются исключительно по три вдоль жидких каналов (границ Плато) под точным углом \(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)\).