Представьте, что вы осторожно достаете из морозильника запечатанную пластиковую бутылку дистиллированной воды: термометр показывает \(-8^\circ\text{C}\) (\(265{,}15\,\text{K}\)), однако вода внутри остается совершенно прозрачной и жидкой. Но стоит слегка ударить дном бутылки о стол — или вылить струю на крошечный осколок льда, — как ослепительно белая сеть игольчатых кристаллов пронзает жидкость со скоростью нескольких сантиметров в секунду, за две секунды превращая всё содержимое в непрозрачную снежную кашу (шугу).

Со школы нам известно правило: при нормальном атмосферном давлении (\(1\,\text{атм} = 101{,}325\,\text{кПа}\)) вода замерзает при \(0^\circ\text{C}\) (\(273{,}15\,\text{K}\)). Однако с точки зрения статистической термодинамики и физики конденсированного состояния \(0^\circ\text{C}\) — это строго равновесная температура плавления (\(T_m\)), при которой макроскопический гексагональный лёд (\(\text{I}_h\)) и жидкая вода имеют одинаковый химический потенциал (\(\mu_{\text{L}}(P, T_m) = \mu_{\text{S}}(P, T_m)\)), но вовсе не температура автоматической кристаллизации воды.

В отсутствие пыли и царапин на стенках сосуда сверхчистая вода остается жидкой при \(-15^\circ\text{C}\), \(-25^\circ\text{C}\) и вплоть до предела гомогенной нуклеации (\(T_H \approx -38^\circ\text{C}\dots{-41^\circ\text{C}}\), или \(232\text{--}235\,\text{K}\)), что доказано экспериментами с использованием фемтосекундных рентгеновских лазеров на свободных электронах (XFEL). Какой же микроскопический барьер не дает воде замёрзнуть сразу ниже \(0^\circ\text{C}\)?

Глава 1: Термодинамика метастабильного состояния

При постоянных температуре \(T\) и давлении \(P\) направление фазового перехода определяется энергией Гиббса \(G = H - TS\) (или молярным химическим потенциалом \(\mu\)). Поскольку молярная энтропия жидкой воды выше, чем у кристаллического льда (\(s_{\text{L}} > s_{\text{S}} > 0\)), из соотношения \(\left(\frac{\partial \mu}{\partial T}\right)_P = -s\) следует, что при всех \(T < T_m\) химический потенциал льда ниже (\(\mu_{\text{S}}(T) < \mu_{\text{L}}(T)\)), то есть лёд отвечает глобальному минимуму свободной энергии.

Тем не менее в интервале от \(0^\circ\text{C}\) до \(-38^\circ\text{C}\) переохлажденная вода (Supercooled Water) оказывается запертой в метастабильном состоянии (локальном минимуме) ландшафта свободной энергии \(G(\eta)\):

Поскольку вторая производная по параметру порядка \(\eta\) строго положительна, любые бесконечно малые тепловые флуктуации лишь повышают свободную энергию и подавляются термодинамическими возвращающими силами. Чтобы перейти в глобальный минимум (лёд), системе необходима локализованная флуктуация конечной амплитуды, способная преодолеть энергетический барьер активации \(\Delta G^* > 0\) — процесс нуклеации (зародышеобразования). Напротив, при плавлении льда (\(T > T_m\)) на его поверхности уже существует нанометровый квазижидкий слой предплавления (Premelting Layer), благодаря чему энергетический барьер отсутствует и лёд начинает таять строго при \(0^\circ\text{C}\).

Глава 2: Классическая теория нуклеации (CNT) — борьба объемного выигрыша и поверхностного штрафа

2.1 Точный вывод \(\Delta G(r)\), критического радиуса \(r^*\) и барьера \(\Delta G^*\)

Разработанная Фольмером, Вебером, Беккером и Дёрингом и обобщенная на переходы жидкость–кристалл Дэвидом Тёрнбуллом и Джоном Фишером (1949), классическая теория нуклеации (CNT) записывает изменение энергии Гиббса \(\Delta G(r)\) при образовании сферического зародыша льда радиуса \(r\) в воде с переохлаждением \(\Delta T \equiv T_m - T > 0\) как:

где \(\Delta g_v = \Delta H_f - T\Delta S_f \approx \frac{\Delta H_f \cdot \Delta T}{T_m}\) — объемная движущая сила кристаллизации (\(\Delta H_f \approx 3{,}06 \times 10^8\,\text{Дж/м}^3\), удельная теплота плавления \(L_f = 333{,}5\,\text{Дж/г}\)), а \(\gamma_{\text{SL}} \approx 28\text{--}32\,\text{мДж/м}^2\) — поверхностное натяжение на границе раздела фаз лёд–вода.

При малых радиусах (\(r \to 0\)) квадратичный поверхностный штраф (\(\propto r^2\)) доминирует над кубическим объемным выигрышем (\(\propto -r^3\)). Приравнивая производную к нулю \(\frac{\partial \Delta G(r)}{\partial r} = -4\pi r^2 \Delta g_v + 8\pi r \gamma_{\text{SL}} = 0\), находим критический радиус зародыша (\(r^*\)) и энергетический барьер нуклеации (\(\Delta G^*\)):

• Докритические кластеры (\(r < r^*\), «эмбрионы»): \(\frac{\partial \Delta G}{\partial r} > 0\). Вследствие эффекта Гиббса–Томсона их локальная температура плавления \(T_m(r) = T_m(1 - \frac{2\gamma_{\text{SL}}}{\Delta H_f r})\) ниже температуры воды \(T\), поэтому они самопроизвольно растворяются.
• Сверхкритические кластеры (\(r > r^*\), «ядра / зародыши»): \(\frac{\partial \Delta G}{\partial r} < 0\). Присоединение новых молекул снижает \(\Delta G\), запуская лавинообразный необратимый рост кристалла.

  flowchart TD
    A["Переохлажденная жидкая вода (T < Tm, ΔT = Tm - T > 0)<br/>Метастабильный локальный минимум"] --> B["Тепловые флуктуации формируют кластер льда<br/>ΔG(r) = -(4/3)πr³Δgv + 4πr²γSL"]
    B --> C{"Сравнение радиуса r с<br/>критическим радиусом r* = 2γSL / Δgv"}
    C -- "r < r* (Докритический эмбрион)" --> D["∂ΔG/∂r > 0 (Доминирует поверхностный штраф)<br/>Локальная точка плавления Tm(r) < T"]
    D --> E["Самопроизвольно тает и<br/>возвращается в жидкую фазу"]
    E --> A
    C -- "r > r* (Сверхкритический зародыш)<br/>Преодолевает барьер ΔG* = 16πγSL³ / 3(Δgv)²" --> F["∂ΔG/∂r < 0 (Доминирует объемный выигрыш)<br/>Локальная точка плавления Tm(r) > T"]
    F --> G["Необратимый макроскопический рост<br/>(Дендриты и выделение скрытой теплоты)"]

Глава 3: Гомогенная и гетерогенная нуклеация — фактор смачивания Тёрнбулла \(f(\theta)\)

Обычная водопроводная вода замерзает при \(-2^\circ\text{C}\dots{-5^\circ\text{C}}\) благодаря гетерогенной нуклеации на стенках или пылинках. В модели сферического сегмента Тёрнбулла зародыш с краевым углом смачивания \(\theta\) (\(\cos\theta = \frac{\gamma_{\text{sub,L}} - \gamma_{\text{sub,S}}}{\gamma_{\text{SL}}}\) по уравнению Юнга) имеет тот же радиус кривизны \(r^*\), но его барьер снижается на геометрический фактор смачивания \(f(\theta)\):

При \(\theta = 30^\circ\) величина \(f(30^\circ) \approx 0{,}013\) (барьер падает до \(1{,}3\%\)); при \(\theta = 10^\circ\) фактор \(f(10^\circ) \approx 1{,}7 \times 10^{-4}\)! Два выдающихся примера эпитаксии:

1. Иодид серебра (\(\beta\text{-AgI}\)) при засеве облаков (Бернард Воннегут, 1947): параметры гексагональной решетки \(\beta\text{-AgI}\) (\(a = 4{,}580\,\text{\AA}, c = 7{,}494\,\text{\AA}\)) отличаются от льда \(\text{I}_h\) (\(a = 4{,}518\,\text{\AA}, c = 7{,}356\,\text{\AA}\)) всего на \(1{,}4\%\), вызывая замерзание капель уже при \(-3^\circ\text{C}\).
2. Белки нуклеации льда (INPs, ген inaZ) бактерии Pseudomonas syringae: \(\beta\)-спиральные повторы на мембране выстраивают молекулы воды в структуру льда и запускают кристаллизацию при \(-1{,}5^\circ\text{C}\dots{-2{,}0^\circ\text{C}\) (применяются в снежных пушках).

Глава 4: Мгновенное замерзание при ударе, дендриты Маллинза–Секерки и закон 10% шуги

При ударе по бутылке с переохлажденной водой (\(-8^\circ\text{C}\)) акустическая волна вызывает акустическую кавитацию: схлопывание микропузырьков по уравнению Рэлея–Плессета порождает локальные скачки давления порядка гигапаскалей (\(\sim\text{ГПа}\)) и адиабатическое охлаждение, мгновенно преодолевая барьер \(\Delta G^*\).

Выделяющаяся при росте льда скрытая теплота плавления (\(L_f = 333{,}5\,\text{Дж/г}\)) нагревает межфазную границу до \(0^\circ\text{C}\). Из-за морфологической неустойчивости Маллинза–Секерки (1963) любой выступ, выдающийся в холодную воду (\(-8^\circ\text{C}\)), попадает в зону более крутого градиента температуры \(|\nabla T|\), быстрее сбрасывает тепло и стремительно растет в виде ветвящихся игл — дендритов — со скоростью \(5\text{--}10\,\text{см/с}\).

Поскольку процесс длится всего \(\sim 2\,\text{с}\), он является адиабатическим (\(\Delta H_{\text{total}} = 0\)): скрытая теплота, выделенная замёрзшей долей массы \(f_{\text{ice}}\), нагревает всю воду в бутылке ровно до \(T_m = 0^\circ\text{C}\), после чего кристаллизация прекращается:

В действительности замерзает лишь \(10{,}1\%\) воды! Этот трехмерный дендритный каркас удерживает остальные \(89{,}9\%\) жидкой воды при \(0^\circ\text{C}\) подобно губке, образуя мягкую снежную кашу (для \(100\%\)-го адиабатического замерзания потребовался бы предел гиперпереохлаждения \(\Delta T_{\text{hyp}} = L_f / c_{p,\text{L}} \approx 79{,}2\,\text{K}\), недостижимый при \(1\,\text{атм}\)).

Глава 5: Перегрев и взрывное вскипание, физика облаков и стеклование (\(T_g\))

1. Перегрев (Superheating) и пузырьковая камера Глейзера: В гладкой кружке в микроволновой печи вода нагревается выше \(100^\circ\text{C}\) без кипения, так как капиллярное давление Юнга–Лапласа \(\Delta P = \frac{2\gamma_{\text{LV}}}{r}\) раздавливает пузырьки пара меньше критического радиуса \(r_b^* = \frac{2\gamma_{\text{LV}}}{P_{\text{V}}(T) - P_{\text{L}}} \approx \frac{2\gamma_{\text{LV}} T_b}{\rho_{\text{V}} L_v \Delta T_{\text{sh}}}\). Добавление ложки сахара или кофе вносит центры парообразования и вызывает взрывное вскипание (bumping). В 1952 году Дональд Глейзер (Нобелевская премия 1960 г.) использовал перегретый жидкий водород в пузырьковой камере для регистрации треков элементарных частиц.
2. Процесс Вегенера–Бержерона–Финдайзена и антифризные белки (AFPs): Ниже \(0^\circ\text{C}\) давление насыщенного пара над льдом меньше, чем над переохлажденной водой (\(e_s^{\text{ice}}(T) < e_s^{\text{water}}(T)\)), поэтому кристаллы льда в облаках растут за счет испарения капель; в авиации при ударе крупных переохлажденных капель (SLD) о крыло мгновенно замерзает лишь доля \(f_{\text{ice}} = c_p\Delta T/L_f\), а остальная вода растекается, образуя опасный прозрачный гололёд (Glaze Ice). Полярные рыбы вырабатывают антифризные гликопротеины (AFGPs), которые закрепляются на гранях льда с шагом \(d\), навязывая кривизну \(r \approx d/2\) и понижая точку замерзания по закону Гиббса–Томсона (\(\Delta T_{\text{TH}} = \frac{4\gamma_{\text{SL}} T_m}{\Delta H_f d}\)).
3. Витрификация (стеклование) и «Ничейная земля» (No-Man’s Land): В формуле скорости нуклеации \(J(T) \propto \exp(-\Delta G_{\text{diff}}/k_B T)\exp(-\Delta G^*/k_B T)\) гигантский рост вязкости \(\eta(T)\) при низких температурах формирует С-образный «нос» на TTT-диаграмме (Время–Температура–Превращение).

  flowchart TD
    A["Стабильная жидкая вода (T > Tm = 273,15 K)"] --> B{"Скорость охлаждения |dT/dt| и<br/>гетерогенные примеси"}
    B -- "Обычное охлаждение + стенки/пыль" --> C["Гетерогенная нуклеация (-1 °C … -20 °C)<br/>Барьер снижен фактором f(θ)"]
    C --> D["Гексагональный кристаллический лёд Ih<br/>(Дендритный рост)"]
    B -- "Чистые микрокапли + обычное охлаждение" --> E["Глубокое переохлаждение до гомогенного предела<br/>TH ≈ -38 °C … -41 °C (232–235 K)"]
    E --> D
    B -- "Сверхбыстрая закалка (Hyperquenching)<br/>|dT/dt| > 10⁶ … 10⁷ K/s" --> F["Обходит нос TTT-диаграммы и пересекает<br/>Ничейную землю (150–235 K)!"]
    F --> G["Переходит температуру стеклования Tg ≈ 136 K (-137 °C)<br/>Вязкость η > 10¹² Па·с"]
    G --> H["Витрифицированный аморфный лёд (LDA / HGW)<br/>Стеклообразная вода без кристаллов (Основа Cryo-EM)"]

При сверхбыстром охлаждении со скоростью выше критической \(\dot{T}_c > 10^6\,\text{К/с}\) (метод погружения в жидкий этан, разработанный Жаком Дюбоше — Нобелевская премия по химии 2017 г. за криоэлектронную микроскопию, Cryo-EM) вода проскакивает «Ничейную землю» (\(150\text{--}235\,\text{K}\)) и ниже температуры стеклования \(T_g \approx 136\,\text{K}\) (\(-137^\circ\text{C}\)) затвердевает без кристаллизации в аморфный лёд низкой плотности (LDA).