Теория ДЛФО: устойчивость коллоидов и барьер

Теория ДЛФО (Дерягина-Ландау-Фервея-Овербека) - это фундаментальный инструмент коллоидной химии, позволяющий предсказать, будет ли суспензия устойчивой или слипнется в осадок. Она объединяет два противоположных взаимодействия между частицами: дальнодействующее электростатическое отталкивание двойного электрического слоя и притяжение Ван-дер-Ваальса. Пока отталкивание создаёт достаточно высокий энергетический барьер, частицы не сталкиваются - суспензия живёт. Добавьте соль или снизьте потенциал поверхности - барьер рухнет, и коагуляция пойдёт быстро. Покрутите ползунки калькулятора ниже, чтобы увидеть, как меняется профиль энергии при разных условиях, - затем разберём каждый член формулы строго.

Откуда берутся два слагаемых ДЛФО

Суммарная энергия взаимодействия двух коллоидных частиц (в приближении сфера-плоскость) в теории ДЛФО записывается как сумма двух вкладов:

$V_{\text{tot}}(h) = V_{\text{EDL}}(h) + V_{\text{vdW}}(h),$

где $h$ - расстояние между поверхностями частиц. Оба слагаемых по-разному зависят от $h$, и их конкуренция формирует характерный профиль с барьером и минимумами.

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

$V_{\text{vdW}}(h) = -\frac{A_H \, R}{12 \, h},$

где $A_H$ - константа Гамакера, характеризующая силу Ван-дер-Ваальса для данной пары материалов в данной среде. Для оксида кремния в воде $A_H \approx 1 \times 10^{-20}$ Дж, для золота в воде $A_H \approx 4 \times 10^{-19}$ Дж - на порядок больше, поэтому золотые наночастицы гораздо труднее стабилизировать.

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

$V_{\text{EDL}}(h) = 64\pi R\,\varepsilon_0\varepsilon_r\left(\frac{k_BT}{ze}\right)^2 \tanh^2\!\!\left(\frac{ze\psi_0}{4k_BT}\right)\exp(-\kappa h),$

где $\psi_0$ - поверхностный потенциал, $z$ - валентность иона, $\kappa$ - обратная дебаевская длина.

Рассчитать для своей задачи

Длина Дебая и роль соли

Ключевой параметр электростатического члена - дебаевская длина $\kappa^{-1}$, которая определяет, на каком расстоянии экранируется поверхностный заряд. Для симметричного 1:1 электролита в воде при 25°C удобная формула:

$\kappa^{-1}\,[\text{нм}] = \frac{0{,}304}{\sqrt{I\,[\text{моль/л}]}}.$

При $I = 1$ мМ длина Дебая $\approx 9{,}6$ нм - отталкивание чувствуется на расстоянии десятка нанометров. При $I = 100$ мМ она падает до $\approx 0{,}96$ нм - электростатика экранируется почти на контакте, барьер исчезает. Именно поэтому добавление соли разрушает коллоидные системы: соль повышает ионную силу, уменьшает дебаевскую длину и «прижимает» отталкивающую кривую к оси $h$. Между тем $V_{\text{vdW}}$ не зависит от ионной силы, поэтому при высокой солёности притяжение побеждает везде.

Рассчитать для своей задачи

Результат зависит ещё и от валентности ионов: правило Шульце-Гарди гласит, что критическая концентрация коагуляции (ККК) убывает примерно как $z^{-6}$. Это означает, что трёхвалентные ионы (Al$^{3+}$, Fe$^{3+}$) коагулируют коллоид в тысячи раз эффективнее, чем одновалентные (Na$^+$, K$^+$) - коллоидная химия промышленного очищения воды опирается именно на это.

Температура тоже влияет на длину Дебая через кТ в числителе формулы для $\kappa^{-1}$: при нагреве длина Дебая незначительно растёт, что немного увеличивает барьер. Однако более существенный температурный эффект связан с тем, что вязкость воды падает, ускоряя броуновское движение частиц. В результате при нагреве суспензия, балансирующая у порога устойчивости, может начать коагулировать быстрее - не из-за изменения профиля энергии, а из-за увеличения частоты столкновений.

Энергетический барьер и устойчивость

Профиль $V_{\text{tot}}(h)$ обычно имеет характерную форму: глубокий первичный минимум при $h \to 0$ (притяжение на контакте), энергетический барьер $V_{\text{max}}$ при $h \sim$ нескольких нанометров и, возможно, неглубокий вторичный минимум при $h \sim 10$-$20$ нм.

Устойчивость системы определяется высотой барьера в единицах тепловой энергии $k_BT$:

• $V_{\text{max}} > 25\,k_BT$ - суспензия устойчива длительно (годы), броуновское движение не преодолевает барьер.
• $V_{\text{max}} = 10$-$25\,k_BT$ - медленная коагуляция (часы-дни), система метастабильна.
• $V_{\text{max}} < 10\,k_BT$ - быстрая коагуляция (минуты), система фактически нестабильна.
• Барьера нет - мгновенная коагуляция (диффузионно-контролируемая, $k \sim 10^{11}$ л/(моль·с)).

Вторичный минимум $V_{\text{sec}}$ интересен сам по себе: если $|V_{\text{sec}}| \sim 1$-$5\,k_BT$, частицы собираются в рыхлые агрегаты, удерживаемые этим минимумом. Такие агрегаты обратимы - они разрушаются встряхиванием или разбавлением. Это принципиально иное поведение по сравнению с коагуляцией в первичный минимум, из которого выйти без специальных условий невозможно.

Потенциал поверхности и дзета-потенциал

В экспериментах прямой доступ к $\psi_0$ затруднён - его не измеришь вольтметром. Вместо него измеряют дзета-потенциал $\zeta$ - потенциал на плоскости скольжения (граница между иммобилизованным и подвижным слоями жидкости). При умеренных ионных силах $\zeta \approx \psi_d$ (потенциал на внешней границе плотного слоя Штерна), и его подставляют в формулу ДЛФО вместо $\psi_0$.

Правило большого пальца: $|\zeta| > 30$ мВ - суспензия устойчива за счёт электростатики; $|\zeta| < 15$ мВ - неустойчива, нужны стерические стабилизаторы (полимеры, ПАВ). Дзета-потенциал измеряется методом электрофореза, и именно он - первая характеристика, которую снимают при разработке наночастиц для медицины, пищевых продуктов или красителей.

Рассчитать для своей задачи

Критическая концентрация коагуляции

Важным практическим результатом теории ДЛФО является предсказание критической концентрации коагуляции (ККК) - концентрации соли, при которой барьер $V_{\text{max}}$ обращается в нуль. Из условий $V_{\text{tot}} = 0$ и $dV_{\text{tot}}/dh = 0$ одновременно можно показать, что ККК убывает как:

$c_{\text{ккк}} \propto \frac{\varepsilon^3 (k_BT)^5}{A_H^2 \, z^6 \, e^6}.$

Главный вывод - зависимость от $z^{-6}$ (правило Шульце-Гарди): замена одновалентного электролита на двухвалентный снижает ККК в $2^6 = 64$ раза, на трёхвалентный - в $3^6 = 729$ раз. Это не просто теоретическая формула: именно поэтому квасцы KAl(SO$_4$)$_2$ со своим Al$^{3+}$ коагулируют взвеси в питьевой воде при концентрациях в сотни раз меньших, чем обычная поваренная соль.

Ограничения теории ДЛФО

Теория ДЛФО даёт правильную физическую картину, но имеет ряд ограничений, которые важно помнить при практическом применении:

• Модель пластинок: строгий вывод сделан для плоских поверхностей; для сфер используется приближение Дерягина, которое справедливо при $\kappa R \gg 1$ (радиус частицы много больше длины Дебая).
• Точечные ионы: теория Пуассона-Больцмана игнорирует конечный размер ионов. При высоких концентрациях солей это приводит к систематическим ошибкам - нужны поправки (модель МСА, теория Горн).
• Нет специфической адсорбции: в модели ДЛФО ионы взаимодействуют с поверхностью только электростатически. Реальные многозарядные ионы (Al$^{3+}$, цитрат) могут адсорбироваться специфически, меняя знак потенциала.
• Нет стерической составляющей: если поверхность покрыта полимером, добавляется стерический барьер, который ДЛФО не учитывает (это зона расширенной теории ХКDLVO).
• Динамические эффекты: теория статическая; гидродинамическое взаимодействие, вязкость и скорость сближения влияют на реальную скорость коагуляции, но не на профиль энергии.

Частые ошибки

• Не переводят ионную силу в моль/л перед подстановкой в формулу $\kappa^{-1} = 0{,}304/\sqrt{I}$. 10 мМ - это 0,01 моль/л, результат в нанометрах.
• Путают $V_{\text{tot}}$ и $V_{\text{EDL}}$: суммарная кривая - это сумма обоих членов, а не просто электростатика. При малых $h$ притяжение всегда побеждает, создавая первичный минимум.
• Игнорируют знак: $V_{\text{vdW}} < 0$ (притяжение), $V_{\text{EDL}} > 0$ (отталкивание). Если перепутать знак Ван-дер-Ваальса, барьер окажется вдвое выше реального.
• Подставляют $\zeta$ без коррекции на высокой ионной силе: при $I > 100$ мМ дзета-потенциал значительно отличается от $\psi_0$ из-за сжатия плотного слоя; подстановка $\zeta$ вместо $\psi_0$ даёт завышенный барьер.
• Путают первичный и вторичный минимум: коагуляция во вторичный минимум обратима; в первичный - необратима. Это определяет, поможет ли разбавление или механическое перемешивание.

FAQ

Что такое константа Гамакера и откуда её брать? Константа Гамакера $A_H$ характеризует суммарную силу дисперсионных взаимодействий для данной пары материалов через данную среду. Её значения табулированы: для SiO$_2$/вода/SiO$_2$ $\approx 1{,}0 \times 10^{-20}$ Дж, для TiO$_2$/вода/TiO$_2$ $\approx 6 \times 10^{-20}$ Дж, для Au/вода/Au $\approx 4 \times 10^{-19}$ Дж. Чем больше $A_H$, тем сильнее притяжение и тем труднее стабилизировать суспензию.

Почему при добавлении многовалентной соли коагуляция наступает при очень малых концентрациях? По правилу Шульце-Гарди критическая концентрация коагуляции убывает как $z^{-6}$. Физически это означает, что трёхвалентный ион, попадая в двойной слой, нейтрализует заряд поверхности в девять раз эффективнее, чем одновалентный, при том же числе ионов. Поэтому 0,1 мМ CaCl$_2$ может разрушить суспензию, устойчивую при 100 мМ NaCl.

Как теория ДЛФО используется на практике? В водоподготовке - для подбора дозы коагулянта (квасцы, сульфат железа) при осаждении взвешенных частиц. В фармацевтике - при создании стабильных суспензий наночастиц для инъекций (нужен $|\zeta| > 30$ мВ и минимальное содержание соли). В пищевой промышленности - для управления стабильностью эмульсий и коллоидов (майонез, молоко). В керамике - для предотвращения агрегации порошков при формовании.

Коротко

Теория ДЛФО описывает устойчивость коллоидных систем через конкуренцию двух сил: отталкивания двойного электрического слоя $V_{\text{EDL}} \propto \exp(-\kappa h)$ и притяжения Ван-дер-Ваальса $V_{\text{vdW}} \propto -1/h$. Дебаевская длина $\kappa^{-1} \approx 0{,}304/\sqrt{I}$ нм определяет, насколько далеко «видит» электростатика; добавление соли сжимает двойной слой и снижает барьер. Система устойчива, пока барьер $V_{\text{max}} > 25\,k_BT$; ниже $10\,k_BT$ - коагуляция неизбежна. Калькулятор выше позволяет отследить, как каждый параметр - потенциал поверхности, константа Гамакера, ионная сила, размер частицы - меняет форму профиля и высоту барьера.