Понимание термодинамических параметров сорбции имеет ключевое значение для оценки самопроизвольности процесса, выбора оптимальных температурных режимов и прогнозирования эффективности очистки в реальных условиях эксплуатации. Листовой опад, рассматриваемый как перспективный природный сорбент для извлечения ионов тяжелых металлов из водных сред, требует детального изучения механизмов взаимодействия с загрязнителями. Термодинамический подход позволяет количественно описать эти взаимодействия, определить вклад физической и химической составляющих, а также установить связь между составом сорбента и его сорбционными характеристиками [1, 2]. Целью настоящей работы является проведение термодинамического анализа сорбции ионов тяжелых металлов (Fe²⁺, Fe³⁺, Cu²⁺, Pb²⁺) листовым опадом различных древесных пород, расчет основных термодинамических параметров и на их основе определение природы и механизма сорбционного процесса.

Термодинамический анализ сорбционных процессов базируется на расчете изменений энергии Гиббса (ΔG), энтальпии (ΔH) и энтропии (ΔS). Изменение энергии Гиббса является критерием самопроизвольности процесса: если ΔG < 0, сорбция протекает самопроизвольно; если ΔG > 0, процесс требует подвода энергии извне. Величина ΔG рассчитывается по уравнению:

ΔG = –RT ln K

где R – универсальная газовая постоянная (8,314 Дж/(моль·К)), T – абсолютная температура (К), K – константа равновесия сорбции [3, 4].

Изменение энтальпии ΔH характеризует тепловой эффект процесса. Отрицательные значения ΔH указывают на экзотермический характер сорбции (выделение тепла), положительные – на эндотермический (поглощение тепла). Величина ΔH позволяет судить о механизме взаимодействия: при |ΔH| < 20 кДж/моль преобладает физическая адсорбция, при 20 < |ΔH| < 80 кДж/моль – физико-химическая сорбция (ионный обмен, водородные связи), при |ΔH| > 80 кДж/моль – химическая сорбция с образованием ковалентных связей [5, 6].

Изменение энтропии ΔS отражает степень упорядоченности системы в процессе сорбции. Положительные значения ΔS указывают на увеличение беспорядка (например, при десорбции молекул растворителя с поверхности сорбента), отрицательные – на упорядочивание системы (фиксация ионов на активных центрах).

Расчет ΔH и ΔS выполняется с использованием линейной формы уравнения Вант-Гоффа:

ln K = –ΔH / (RT) + ΔS

Построение графика зависимости ln K от 1/T позволяет определить ΔH по углу наклона прямой, а ΔS – по отсекаемому отрезку на оси ординат [3, 4].

Объектами исследования служили образцы листового опада пяти древесных пород, широко распространенных в озеленении Санкт-Петербурга и Ленинградской области: березы повислой (Betula pendula), осины обыкновенной (Populus tremula), клена остролистного (Acer platanoides), дуба черешчатого (Quercus robur) и липы мелколистной (Tilia cordata). Сбор опада производился в период массового листопада, образцы высушивались до воздушно-сухого состояния и измельчались. В качестве сорбатов использовались водные растворы солей железа (II), железа (III), меди (II) и свинца (II) с исходной концентрацией 50 мг/л. Выбор этих металлов обусловлен их распространенностью в промышленных сточных водах и различной способностью к комплексообразованию.

Сорбционные эксперименты проводились в статических условиях при температурах 293, 303, 313, 323 и 333 К. Навеску сорбента массой 1,0 г помещали в 100 мл раствора металла с заданным pH (оптимальное значение 5,5, установленное в предварительных экспериментах) и выдерживали в течение времени, достаточного для достижения равновесия (120 мин). После завершения процесса раствор фильтровали, концентрацию металлов определяли методом атомно-абсорбционной спектроскопии. Сорбционную емкость рассчитывали как разность между начальной и равновесной концентрациями, отнесенную к массе сорбента.

Для каждой системы рассчитывали константу распределения Kd и константу равновесия K. Затем строили графики зависимости ln K от 1/T, определяли ΔH и ΔS, после чего вычисляли ΔG для каждой температуры.

Экспериментальные данные показали, что сорбционная емкость всех исследованных образцов листового опада закономерно изменяется с повышением температуры. Для опада дуба черешчатого, который характеризуется максимальным содержанием лигнина (60,2%) и танинов (12,46%), отмечено увеличение сорбционной емкости по ионам Fe³⁺ с 9,8 мг/г при 293 К до 12,4 мг/г при 333 К. Аналогичная тенденция, но с меньшими абсолютными значениями, наблюдалась для березы и клена, тогда как опад осины и липы демонстрировал наименьший рост сорбционной емкости с температурой, что связано с их меньшим содержанием лигнина [2, 7].

Расчет термодинамических параметров позволил количественно охарактеризовать процесс сорбции. В таблице 1 представлены значения ΔG, ΔH и ΔS для сорбции ионов Fe³⁺ исследуемыми образцами при 293 К.

Таблица 1

Термодинамические параметры сорбции ионов Fe³⁺ листовым опадом при 293 К

Отрицательные значения ΔG для всех образцов подтверждают самопроизвольный характер сорбции ионов железа листовым опадом. Наиболее отрицательное значение ΔG отмечено для дуба черешчатого (–7,6 кДж/моль), что согласуется с его высокой сорбционной активностью, обусловленной богатым химическим составом. Для осины и липы значения ΔG близки к нулю (–3,9…–4,1 кДж/моль), что указывает на меньшую термодинамическую предрасположенность к сорбции.

Анализ значений ΔH показывает, что все исследованные системы характеризуются отрицательной энтальпией, то есть процесс сорбции является экзотермическим и протекает с выделением тепла. Для дуба черешчатого величина ΔH составляет –32,5 кДж/моль, что попадает в диапазон 20–80 кДж/моль и свидетельствует о преобладании физико-химического механизма сорбции с участием ионного обмена и комплексообразования. Для березы, клена, осины и липы значения |ΔH| находятся в интервале 14–18 кДж/моль, что соответствует физической адсорбции с ван-дер-ваальсовыми взаимодействиями и водородными связями. Таким образом, вклад хемосорбционной составляющей максимален для дуба, что объясняется высоким содержанием фенольных гидроксилов лигнина и орто-дигидроксильных групп танинов, способных к образованию прочных комплексов с ионами металлов [5, 7].

Значения ΔS для всех образцов отрицательны, что указывает на уменьшение энтропии системы в процессе сорбции. Это объясняется тем, что при фиксации ионов металлов на активных центрах поверхности сорбента происходит упорядочивание системы: ионы переходят из объема раствора, где они находятся в подвижном состоянии, в связанное состояние на поверхности твердой фазы. Наиболее отрицательное значение ΔS отмечено для дуба (–84,7 Дж/(моль·К)), что связано с более жесткой фиксацией ионов на многочисленных функциональных группах.

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

Таблица 2

Термодинамические параметры сорбции различных ионов металлов листовым опадом дуба черешчатого при 293 К

Наиболее отрицательные значения ΔG характерны для ионов свинца и меди, что указывает на их более высокое сродство к функциональным группам листового опада по сравнению с железом. Это объясняется большей способностью Pb²⁺ и Cu²⁺ к комплексообразованию с фенольными и карбоксильными группами [6]. Значения ΔH для этих ионов также более отрицательны, что свидетельствует о более прочном связывании. Для трехвалентного железа отмечены более отрицательные значения ΔG и ΔH по сравнению с двухвалентным, что связано с более высоким зарядом и, как следствие, большей способностью к гидролизу и образованию гидроксокомплексов, которые эффективно взаимодействуют с поверхностью сорбента [2].

Полученные термодинамические параметры позволяют установить четкую связь между химическим составом листового опада и его сорбционными свойствами. Опад дуба черешчатого, содержащий максимальное количество лигнина (60,2%) и танинов (12,46%), демонстрирует наиболее отрицательные значения ΔG и наибольшие значения |ΔH|, что указывает на значительный вклад хемосорбционной составляющей в процесс извлечения ионов металлов. Фенольные гидроксилы лигнина и орто-дигидроксильные группы танинов выступают в роли эффективных хелатирующих агентов, образуя с катионами металлов устойчивые пяти- и шестичленные циклы [5, 7]. Опад березы и клена, содержащий умеренное количество этих компонентов, характеризуется промежуточными термодинамическими показателями, а опад осины и липы, с минимальным содержанием лигнина, проявляет наименьшее сродство к металлам и преобладающий физический механизм сорбции.

Влияние природы иона металла на термодинамику сорбции проявляется в закономерности: чем выше заряд иона и больше его способность к комплексообразованию, тем более отрицательны значения ΔG и ΔH. Это объясняется тем, что ионы с более высоким зарядом создают более сильное электростатическое поле, которое притягивает полярные функциональные группы сорбента, а также более активно участвуют в реакциях гидролиза с образованием гидроксокомплексов, способных к эффективному взаимодействию с поверхностью [2, 6].

Практическое значение полученных результатов заключается в возможности прогнозирования эффективности сорбционной очистки при различных температурах. Поскольку все исследованные системы характеризуются экзотермическим характером процесса (ΔH < 0), повышение температуры будет снижать равновесную сорбционную емкость. Это означает, что для достижения максимальной эффективности очистки целесообразно проводить процесс при умеренных температурах (20–30 °C). Исключение могут составлять случаи, когда требуется ускорить кинетику сорбции за счет повышения температуры, однако это будет сопровождаться некоторым снижением равновесной емкости. Для горячих сточных вод (например, после технологических процессов) необходимо учитывать снижение сорбционной способности и при необходимости увеличивать дозу сорбента или время контакта.

Проведенный термодинамический анализ сорбции ионов тяжелых металлов листовым опадом пяти древесных пород позволяет сделать следующие выводы. Для всех исследованных систем характерны отрицательные значения энергии Гиббса (ΔG от –3,9 до –9,1 кДж/моль), что подтверждает самопроизвольный характер сорбции, причем наиболее отрицательные значения ΔG отмечены для опада дуба черешчатого, что согласуется с его высоким содержанием лигнина и танинов. Процесс сорбции является экзотермическим (ΔH < 0), при этом значения энтальпии для дуба черешчатого составляют –28,3…–38,6 кДж/моль в зависимости от сорбируемого иона, что указывает на преобладание физико-химического механизма с участием ионного обмена и комплексообразования. Для остальных пород значения |ΔH| находятся в диапазоне 14–18 кДж/моль, что соответствует физической адсорбции. Отрицательные значения изменения энтропии (ΔS от –36,4 до –100,2 Дж/(моль·К)) свидетельствуют об упорядочивании системы при фиксации ионов металлов на поверхности сорбента. Установлено, что термодинамические параметры зависят как от вида древесной породы (в первую очередь от содержания лигнина и танинов), так и от природы сорбируемого иона (наибольшее сродство проявляют ионы Pb²⁺ и Cu²⁺, наименьшее – Fe²⁺). Полученные результаты имеют практическое значение для прогнозирования эффективности очистки в различных температурных условиях и обоснованного выбора сорбента для конкретных технологических задач.