Представьте, что вы осторожно достаете из морозильника запечатанную пластиковую бутылку дистиллированной воды: термометр показывает $-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).