# Термодинамика стратификации: от океанских озер CO2 до спинодального распада

Источник: https://www.youtube.com/watch?v=2Q5ZFIms9Bo
Канал: MIT OpenCourseWare
Опубликовано: 05.12.2024

---

В лекции профессора Массачусетского технологического института (MIT) Джана Паоло Беретты подробно рассматриваются физико-химические механизмы стратификации в газовых и жидких смесях под воздействием различных силовых полей. Автор анализирует, как гравитационные, центробежные и электрические силы влияют на пространственное распределение компонентов, и переходит к термодинамическому описанию фазовых переходов на основе модели Ван-дер-Ваальса. Центральной темой материала становится концепция метастабильных состояний и спинодального распада, объясняющая как масштабные природные явления, так и принципы работы фундаментальных научных приборов.

## 🧪 Разделение смесей при постоянных условиях: термодинамический парадокс?
[[JUMP:0:14]]

Процессы смешения и разделения компонентов лежат в основе множества химических и физических технологий. Один из ключевых вопросов термодинамики заключается в том, как осуществить разделение смеси при постоянных температуре и давлении. Известно, что энтропия смешения идеальных газов всегда положительна, а значит, обратный процесс — полное разделение смеси в условиях взаимного равновесия с окружающей средой (термостатом) — требует совершения минимальной работы.

Эта минимальная работа численно равна разности свободных энергий Гиббса конечного (разделенного) и исходного (смешанного) состояний смеси. Графически этот баланс можно представить через изменение энтропии: поскольку температура и внутренняя энергия идеальных газов при смешении не меняются, разделенное состояние обладает меньшей энтропией, чем исходная смесь. Чтобы осуществить разделение, необходимо отвести избыточную энтропию. Для этого требуется внешняя система (резервуар), которая поглотит эту энтропию, но резервуар не может принять энтропию без получения соответствующего количества энергии, равного произведению температуры резервуара на величину отводимой энтропии.

Теоретически такой процесс можно реализовать обратимым путем через предварительное охлаждение смеси. На практике для этого потребовалось бы сложное устройство:

* Сначала смесь обратимо охлаждается до температуры ниже температуры резервуара, при этом машина потребляет работу, равную доступной энергии данного состояния.
* Затем проводится адиабатический процесс без теплообмена, где с помощью полупроницаемых мембран (каждая из которых проницаема только для одного типа частиц) веса опускаются и совершают работу по парциальному сжатию компонентов смеси.

Поскольку использование идеальных полупроницаемых мембран в реальности крайне затруднено, в инженерной практике чаще применяют методы химической технологии. Проводятся реакции с селективными растворителями, которые избирательно связывают один из компонентов смеси. После разделения химическую реакцию обращают, чтобы регенерировать растворитель и вернуть его в исходное состояние. Тем не менее, даже в таком циклическом химическом процессе вычисленная термодинамическая минимальная работа остается главным ориентиром для расчета эффективности разделения по второму закону термодинамики.

## 🌍 Гравитационная стратификация: почему на Эвересте тяжело дышать
[[JUMP:9:34]]

Переходя к явлению стратификации, профессор Беретта предлагает рассмотреть распределение газов в атмосфере Земли (кислорода, азота и углекислого газа). В рамках термодинамической модели полную энергию элементарного объема флюида представляют как сумму его внутренней энергии в состоянии устойчивого равновесия и потенциальной энергии в однородном гравитационном поле, которая выражается классической формулой $mgz$. Локальная масса флюида изменяется с высотой из-за изменения химического состава смеси.

В этих условиях традиционный химический потенциал заменяется полным (гравихимическим) потенциалом, который учитывает как чисто химическую составляющую, так и вклад силы тяжести. Условием взаимного устойчивого равновесия между слоями флюида на разных высотах является строгое равенство их полных потенциалов. Принимая приближение, что локальные изменения температуры не влияют на химическое равновесие, равенство потенциалов записывается для высоты $z$ и $z + dz$. Для химической части используется модель идеальной смеси. Математические манипуляции приводят к дифференциальному уравнению, связывающему изменение мольной доли компонента, градиент давления и молярную массу.

Для чистой однокомпонентной субстанции это уравнение переходит в известное гидростатическое соотношение Бернулли. Для идеального газа, у которого удельный объем зависит от давления, интегрирование дает классический закон экспоненциального убывания давления с высотой. В случае газовой смеси аналогичный закон выполняется для парциального давления каждого отдельного компонента. При этом показатель экспоненты напрямую зависит от молярной массы газа:

* Более тяжелые молекулы опускаются вниз быстрее, и их парциальное давление с высотой падает стремительнее.
* Молекулярный кислород имеет молярную массу 32 г/моль, а молекулярный азот — 28 г/моль. 

Из-за этой разницы парциальное давление кислорода уменьшается с высотой быстрее, чем давление азота, и мольная доля кислорода падает с 21% на уровне моря до более низких значений на больших высотах.

Общее давление смеси в любой точке определяется по закону Дальтона как сумма парциальных давлений. На вершине Эвереста (около 8848 метров) атмосферное давление составляет всего около 32% от давления на уровне моря. Это означает, что при каждом вдохе человек получает лишь треть от привычного количества воздуха. Данный эффект усугубляется падением мольной доли самого кислорода в составе разреженного воздуха.

## 🏔️ Климатические хроники: от Мауна-Лоа до жидких озер CO2 на дне океана
[[JUMP:27:03]]

Данные многолетних наблюдений наглядно иллюстрируют законы гравитационной стратификации. Профессор Беретта демонстрирует расчеты состава воздуха, верифицированные по сведениям высокогорной лаборатории Мауна-Лоа (3400 метров над уровнем моря). Сегодня концентрация углекислого газа (CO2) там составляет около 425 ppm (частей на миллион). Согласно термодинамическому расчету, на уровне моря из-за высокой молярной массы CO2 (44 г/моль) его концентрация должна превышать 500 ppm, тогда как на вершине Эвереста она значительно ниже.

Анализ исторических графиков с научного ресурса Climate 4 You подтверждает устойчивый рост концентрации CO2 в атмосфере с 1955 года по настоящее время. Одновременно фиксируется рост температуры верхнего 100-метрового слоя Мирового океана — зоны интенсивного турбулентного перемешивания. Между атмосферой и этим слоем непрерывно поддерживается стремление к термодинамическому равновесию, вследствие чего CO2 растворяется в воде. Расчет усложняется тем, что углекислый газ в воде ионизируется, превращаясь в бикарбонаты и карбонаты.

Однако на глубине около 450 метров физические условия меняются настолько, что CO2 переходит в жидкое состояние. С этого момента газовая стратификация сменяется стратификацией жидких смесей. Поскольку жидкий диоксид углерода плотнее морской воды, он опускается на самый качественный «склад» планеты — на океанское дно, где формируются настоящие подводные озера из жидкого и сверхкритического CO2.

В научной среде, как отмечает лектор, ведутся дискуссии вокруг временного лага между пиками температур океана и пиками концентрации атмосферного CO2. Сжигание ископаемого топлива поставляет огромные объемы газа, из которых около 40% поглощается океаном через механизмы межфазного равновесия. Помимо CO2, исследователей беспокоят другие парниковые газы. Газета MIT Tech Talk недавно опубликовала статью об элегазе — гексафториде серы (SF6), чей парниковый потенциал в 24 000 раза выше, чем у CO2. Будучи чрезвычайно тяжелой, эта молекула под действием гравитации должна интенсивно концентрироваться в приземных слоях, однако точные данные о ее высотном распределении на специализированных ресурсах вроде Climate 4 You пока не обобщены. 

Дополнительным фактором глобального климата выступает циклическая солнечная активность (число солнечных пятен). По мнению ряда исследователей климата, изменения солнечной активности влияют на интенсивность космических лучей, которые, в свою очередь, модулируют процесс образования облаков. Облака определяют альбедо Земли, отражая солнечный свет и вызывая охлаждение планеты.

## 🌀 Центрифуги и изотопы: стратификация во вращающихся системах
[[JUMP:35:14]]

Математический аппарат, разработанный для гравитационного поля, полностью применим к анализу стратификации во вращающихся системах, где на флюид действуют центробежные силы. Центробежная сила является инерционной (фиктивной) и может быть описана через отрицательную потенциальную энергию, пропорциональную квадрату угловой скорости и квадрату радиуса: $-\frac{1}{2} m \omega^2 r^2$. Полный потенциал компонента во вращающемся флюиде включает этот отрицательный энергетический вклад.

При установлении взаимного равновесия между слоями на разных расстояниях от оси вращения интегрирование уравнений состояния показывает, что давление в центробежном поле растет по экспоненциальному закону для газа и по параболическому — для жидкости. В газовых смесях парциальные давления компонентов перераспределяются аналогично гравитационному разделению.

Этот физический принцип критически важен для мировой энергетики, поскольку лежит в основе технологии обогащения урана. Природный уран содержит всего 0,71% fissile-изотопа U-235, тогда как для работы ядерных реакторов требуется повысить его концентрацию до 3–3,5%. Для разделения изотопы переводят в газообразную форму — гексафторид урана (UF6). Газ подается в высокоскоростные центрифуги. Ничтожная разница в молярной массе между UF6 с ураном-235 и ураном-238 оказывается достаточной, чтобы более тяжелые молекулы под действием центробежной силы концентрировались у периферии ротора.

Поскольку за один проход через центрифугу степень обогащения увеличивается незначительно, на практике выстраивают каскады из тысяч последовательных ступеней. На выходе из каскада получают целевой обогащенный уран (3%) и обедненный уран с содержанием U-235 около 0.2%. Хранение огромных объемов обедненного урана сопряжено со значительными рисками и затратами: стальные контейнеры со временем ржавеют, что требует постоянных инвестиций в безопасность.

Несмотря на эти издержки, ядерная энергетика остается экономически высококонкурентной. Профессор Беретта приводит личный пример: в Италии по итогам референдума 1989 года были остановлены все атомные электростанции. На одной из совершенно новых, но так и не запущенных АЭС спустя 35 лет продолжают работать около 200 человек, чья единственная задача — поддерживать безопасное законсервированное состояние выведенного из эксплуатации объекта.

## ⚡ Электрохимический потенциал: поведение заряженных частиц в поле
[[JUMP:43:34]]

Третьим примером поля, вызывающего пространственную стратификацию компонентов, является электростатическое поле. Если в смеси присутствуют заряженные частицы (ионы в электролите), в термодинамику вводится понятие электрохимического потенциала. Полный заряд компонента определяется валентностью ионов $z$ и постоянной Фарадея $F$, которая задает заряд одного моля электронов. Например, для иона натрия $Na^+$ величина $z = +1$, а для электрона $z = -1$.

Согласно современным стандартам Международной системы единиц (СИ), фундаментальные константы зафиксированы с абсолютной точностью:

* Постоянная Авогадро и постоянная Больцмана имеют строго определенные численные значения без округлений.
* Произведение числа Авогадро на заряд электрона дает постоянную Фарадея (около 96 500 Кл/моль).
* Произведение числа Авогадро на постоянную Больцмана дает универсальную газовую постоянную $R$. Профессор Беретта подчеркивает, что историческое название «газовая постоянная» вводит в заблуждение, поскольку эта величина играет фундаментальную роль в термодинамике жидких и твердых тел в той же мере, что и в газах.

Энергия, добавляемая к термодинамической системе за счет электрического поля, равна произведению заряда смеси на электростатический потенциал $\phi$. Градиент этого потенциала, взятый с обратным знаком, представляет собой вектор электрического поля. При расчете условий равновесия раствора электролита в электрическом поле равенство полных электрохимических потенциалов приводит к выражениям, математически идентичным формулам гравитационной стратификации.

В простейшем случае плоского конденсатора (между анодом и катодом) электростатический потенциал $\phi$ изменяется линейно, а напряженность электрического поля постоянна. Интегрирование уравнений дает экспоненциальный профиль распределения концентраций: плотность положительных ионов экспоненциально убывает по мере удаления от отрицательного электрода (катода), а концентрация отрицательных ионов симметрично падает при удалении от положительного анода.

## 📉 Модель Ван-дер-Ваальса и границы устойчивости вещества
[[JUMP:55:00]]

Для описания реального поведения чистых веществ и смесей при фазовых переходах (кипении, конденсации) современные исследователи активно используют компьютерное моделирование методом молекулярной динамики или Монте-Карло, базирующееся на потенциале Леннард-Джонса. В физической литературе возникла концепция моделирования псевдосубстанции, получившей название «леннардджонсиум» (Lennardjonesium), свойства которой рассчитываются в безразмерных координатах плотности и температуры.

Критическая точка леннардджонсиума характеризуется универсальными безразмерными параметрами. Зная реальные критические давление и температуру вещества, можно составить систему уравнений и вычислить индивидуальные параметры потенциала (эффективный диаметр молекулы $\sigma$ и глубину потенциальной ямы $\epsilon$). Однако двухпараметрическая модель Леннард-Джонса не является абсолютно универсальной: она дает погрешность при определении критического объема. Безразмерный фактор сжимаемости $Z = \frac{p_c v_c}{R T_c}$ для леннардджонсиума равен 0,31. Для реального гелия-4 это значение совпадает идеально, но для углекислого газа оно отличается, а для воды составляет 0,23. Модель не учитывает многочастичные (трехтельные) взаимодействия, проявляющиеся в плотной жидкой фазе.

Исторической предтечей этих расчетов стала докторская диссертация Яна Дидерика Ван-дер-Ваальса (Нобелевская премия 1910 года), предложившего феноменологическое уравнение состояния с двумя параметрами:

* Параметр $b$ (кообъем) учитывает собственный объем молекул, ограничивающий свободное пространство для их движения.
* Параметр $a$ характеризует силы межмолекулярного притяжения.

Уравнение Ван-дер-Ваальса позволяет рассчитывать изотермическую сжимаемость, коэффициент изобарного расширения и коэффициент Джоуля-Томсона, определяющий, будет ли реальный газ охлаждаться или нагреваться при дросселировании. 

Если построить изотермы Ван-дер-Ваальса на диаграмме «давление — объем» ($p-v$) при температурах ниже критической, кривые приобретают волнообразный вид с максимумом и минимумом. Согласно принципу Ле Шателье — Брауна, условием механической устойчивости вещества является положительная сжимаемость, то есть производная давления по объему вдоль изотермы должна быть строго отрицательной: $\left(\frac{\partial p}{\partial v}\right)_T < 0$. На волнообразном участке изотермы Ван-дер-Ваальса присутствует область, где эта производная становится положительной. Данная область абсолютно неустойчива: вещество в этих координатах физически не может существовать как однородная фаза и распадается на плотную жидкость и разреженный пар.

## 🧼 Метастабильность и нобелевские камеры: как зафиксировать элементарные частицы
[[JUMP:1:14:42]]

Для строгого описания равновесия сосуществующих фаз Ван-дер-Ваальс использовал равенство химических потенциалов, давлений и температур жидкой и паровой фаз. Интегрирование уравнения $\left(\frac{\partial \mu}{\partial p}\right)_T = v$ позволяет рассчитать химический потенциал модели Ван-дер-Ваальса как функцию температуры и объема. Совместное решение уравнений равенства давлений и равенства химических потенциалов определяет строгие объемы насыщенного пара $v_g$ и насыщенной жидкости $v_f$ при заданной температуре.

График разности химического потенциала и его значения на линии насыщения ($\mu - \mu_{sat}$) ниже критической температуры имеет характерную петлю. Точка пересечения ветвей соответствует истинному термодинамическому равновесию — насыщенному пару и насыщенной жидкости с минимальной свободной энергией Гиббса. Однако участки петли Ван-дер-Ваальса продолжаются дальше точек насыщения, сохраняя отрицательную производную $\left(\frac{\partial p}{\partial v}\right)_T < 0$. Эти состояния термодинамически невыгодны по сравнению с расслоенной системой, но механически устойчивы к малым возмущениям. Они называются метастабильными.

Метастабильные состояния разделяются на два типа:

1.  **Перегретая жидкость** — жидкость, находящаяся при давлении ниже давления насыщения (или при температуре выше температуры кипения) для данной изотермы. Она «не знает», что ей выгоднее испариться, поскольку для зарождения первой газовой микропузырьковой фазы требуется преодолеть энергетический барьер — затратить энергию на создание межфазной поверхности (поверхностное натяжение).
2.  **Пересыщенный (переохлажденный) пар** — пар, находящийся при давлении выше давления насыщения, которому для конденсации необходимо сформировать первую жидкую каплю.

Если метастабильной системе сообщить начальный триггер (вибрация, излучение), избыточная свободная энергия мгновенно высвобождается, и запускается лавинообразный, взрывной процесс фазового перехода.

Этот фундаментальный термодинамический эффект принес ученым еще две Нобелевские премии за создание детекторов элементарных частиц. Чарльз Вильсон изобрел камеру Вильсона (cloud chamber), работающую на переохлажденном метастабильном паре. Дональд Глейзер создал пузырьковую камеру (bubble chamber), заполненную перегретой метастабильной жидкостью. Пролетающая через объем прибора высокоэнергетическая заряженная частица ионизирует среду на своем пути. Эта энергия ионизации служит триггером, снимающим метастабильность: вдоль трека частицы мгновенно вскипают микропузырьки или конденсируются капли, визуализируя траекторию ее движения. Аналогичный механизм реализуется в масштабах планеты, когда космические лучи, бомбардирующие атмосферу Земли, поставляют энергию для нуклеации капель воды, стимулируя глобальное облакообразование.

## 💥 Спинодальный распад: взрывное разделение фаз
[[JUMP:1:34:42]]

Что произойдет, если перевести вещество из стабильного состояния в абсолютно неустойчивую область изотермы Ван-дер-Ваальса настолько быстро, чтобы система не успела претерпеть обычную нуклеацию (рост зародышей) на метастабильном этапе? Экспериментально это достигается мгновенным сбросом давления (квенчем) или резким охлаждением тончайшего слоя флюида, зажатого между стеклами в лабораторных условиях.

Граница, отделяющая область метастабильных состояний от зоны абсолютной термодинамической неустойчивости (где $\left(\frac{\partial p}{\partial v}\right)_T = 0$), называется спинодалью. Внутри спинодальной кривой вещество полностью теряет фазовую однородность. Компьютерное моделирование демонстрирует различные морфологические сценарии спинодального распада в зависимости от исходного объема:

* Если система переводится в неустойчивое состояние строго через критический удельный объем, разделение фаз происходит по «бинепрерывному» (bicontinuous) типу. Вещество мгновенно формирует сложнейшую переплетенную трехмерную структуру, где все области плотной жидкой фазы соединены между собой, и точно так же взаимосвязаны все области паровой фазы.
* Если перевод осуществляется со смещением относительно критического объема, вещество сразу распадается на изолированные фракции — капли жидкости или пузыри пара.

Процесс спинодального распада протекает чрезвычайно быстро, поскольку свободная энергия системы расходуется напрямую на спонтанное разделение фаз по всему объему. Сразу после завершения начального распада включается механизм коалесценции: мелкие пузыри или капли начинают сливаться друг с другом, образуя более крупные конгломераты. Таким способом термодинамическая система уменьшает суммарную площадь межфазной поверхности, минимизируя поверхностную энергию и приходя к итоговому стабильному расслоению на две макрофазы.