УДК 544.636, 541.135.5
ЭЛЕКТРОХИМИЧЕСКИЕ СВОЙСТВА АКТИВИРОВАННОГО ДРЕВЕСНОГО УГЛЯ
В ЩЕЛОЧНОМ ЭЛЕКТРОЛИТЕ
Б. П. Бахматюк, А. С. Курепа
Национальный университет «Львовская политехника», Украина
E-mail: bgbakhm@rambler.ru Поступила в редакцию 07.12.11 г.
В работе исследованы ёмкостные и кинетические свойства нанопористого углеродного материала с максимумом распределения по радиусам пор 1.37 им, полученного из древесины. Используя данные рассеяния рентгеновских лучей и электрохимической импедансной спектроскопии, исследовано влияние изменения пористой структуры и электронного строения активированного углеродного материала на механизм и кинетику процессов заряда-разряда в 7.6 М растворе КОН. Установлено два различных механизма заряда пористой структуры исследованного материала, зависящие от электродного потенциала и химических потенциалов ионов электролита. Первый механизм представляет собой процесс электростатической адсорбции гидратированных ионов электролита, а второй является процессом электросорбции ионов Н+ или ОН~. Показано, что за две секунды можно осуществить цикл заряд-разряд ёмкости 95 Ф/г, что составляет половину максимальной ёмкости исследованного материала.
Ключевые слова: суперконденсатор, импеданс, кинетика, нанопористый углеродный материал, электрическая эквивалентная схема, двойной электрический слой.
We have studied the capacitance and kinetic properties of nanoporous carbon material which was obtained from wood. The maximum radius pore distribution of the investigated material was 1.37 nm. The data from X-ray scattering and electrochemical impedance spectroscopy was used to investigate the influence of porous structure change and electronic structure of activated carbon material on the mechanism and kinetics of charge-discharge at 7.6m KOH solution. It was proved that depending on the electrode potential and chemical potentials of electrolyte ions there are two different mechanisms of the charge of the porous structure of investigated material. The first mechanism is a process of electrostatic adsorption of hydrated electrolyte ions, and the second is a process electrosorption of H+ or OH. It was shown the cycle of charge-discharge capacity of 95 F/g can be made for two seconds and this is half of the maximum capacity of the material under study.
Key words: ultracapacitor, impedance, kinetic, nanoporous carbon material, electrical equivalent circuit, electric double layer.
ВВЕДЕНИЕ
Повышенный интерес к электрохимическим свойствам активированных углеродных материалов обусловлен разработкой двухслойных конденсаторов (ДСК), основанных на электродах из активированного углеродного материала (АУМ). Первый патент на суперконденсатор был получен в 1957 г. [1]. В СССР первые ДСК были разработаны под руководством Н. С. Лидоренко и А. М. Иванова [2]. Они были названы «молекулярными накопителями электрической энергии» или «ионисторами». ДСК работают на хорошо обратимых процессах заряда-разряда электростатической и специфической ёмкости на границе раздела высокоразвитой поверхности АУМ с электролитом [3]. Одним из основных электролитов для таких конденсаторов является 7.6 М КОН, позволяющий получить максимальное рабочее напряжение 1.0 В [4]. В работе [5] было определено четыре области электродных потенциалов на границе раздела обеззоленного и обесгазированного активированного угля с 0.01 М КОН, что обусловлено процессами формирования двойного электрического слоя (ДЭС) гидратированными ионами, адсорбцией
водорода и гидроксильных групп, а также выделением газообразного кислорода. Эти адсорбционные механизмы формирования ДЭС были использованы в исследованиях [6] для интерпретации работы ДСК на основе промышленных активированных углей в 7.6 М КОН без поправки на изменение химических потенциалов ионов водорода и гидроксила. Кроме механизмов процессов формирования ДЭС, ёмкость системы зависит от доступности пористой структуры АУМ для ионов электролита [7, 8]. Распределение пористой структуры по размерам пор в значительной степени влияет на кинетику процессов заряд-разряд ДСК [9]. Поэтому в данной работе исследовано влияние изменения электронной и пористой структуры АУМ и концентрации электролита на механизм, ёмкостные характеристики и кинетику заряда ДЭС в 7.6 М КОН.
МЕТОДИКА ЭКСПЕРИМЕНТА
В работе исследован обеззоленный древесный уголь (зольность 0.17%), полученный высокотемпературной карбонизационной активацией водным паром при 950°С древесины. Процентное содержа-
© БАХМАТЮК Б. П., КУРЕПА А. С., 2011
ние вычислялось по массе золы, которая осталась в лодочке после полного сгорания угля. Для обез-золения активированного древесного угля (АДУ) он был подвергнут химической обработке в кислотах. Сначала его выдерживали в концентрированной плавиковой кислоте на протяжении 22-24 часов. После чего уголь фильтровали и промывали до достижения рН = 6-7. Затем такая же обработка была сделана при помощи азотной и хлороводородной кислот. После обработки уголь сушили при температуре 100-105°С до постоянной массы. Зольность полученного таким образом АДУ уменьшилась с 2.3 до 0.17%. Для десорбции активного материала и электродов была применена термовакуумная обработка, которую проводили в температурном диапазоне 200-250°С при остаточном давлении 10-1 10-2 атм.
Для изучения пористой структуры АДУ был использован метод малоуглового рассеяния рентгеновских лучей [10]. Эксперименты проводились на рентгеновском дифрактометре ДРОН-3 в Си-Ка-из-лучении (X = 1.5418 Е), монохроматизированном отражением от плоскости (111) монокристалла Ge в режиме прохождения пучка рентгеновского излучения. Толщина исследованных образцов составляла 1 мм. Для ограничения зоны паразитного рассеяния монокристалла-монохроматора и уменьшения фона рассеяния на воздухе перед образцом и детектором устанавливались щелевые устройства перпендикулярно к направлению рентгеновского пучка. Использование качественного монокристалла Ge и ограничительной системы первичного и рассеянного излучения позволяет проводить измерения спектров малоуглового рассеяния начиная со значений волнового вектора 0.01 Е-1. Перед детектором устанавливалась щель 0.1 мм, что соответствует пространственному разрешению А(2©) = 0.03°. Регистрация рассеянного излучения проводилась в режиме сканирования в угловом интервале 0.25^.00° с шагом 0.05°, время экспозиции — 100 с. Для получения дифракционных кривых в широком угловом интервале перед детектором устанавливали щель размером 1.00 мм.
Компоненты электродной массы — активированный уголь, ацетиленовая сажа (измельчённые до размеров частиц 60-80 мкм и 3-5 мкм соответственно) смешивались в пропорциях 85 мас.%: 10 мас.% с добавкой 5 мас.% связующего агента, в качестве которого применялся поливинилиденфто-рид. К полученной смеси добавляли ацетон, доводя её до пастообразной консистенции. Пастообразная электродная композиция прессовалась на металлическую основу (никель или нержавеющая сталь). Для электрохимических исследований использовались электроды с рабочей поверхностью 0.5 см2 и массой активного материала 4-10-3 г.
Исследование электрохимических свойств проводилось методом электрохимической импедансной спектроскопии (в диапазоне частот 10~2-8-103 Гц с амплитудой переменного тока 5 мВ) при потенци-остатических условиях с помощью измерительного комплекса «AUTOLAB» фирмы «ECO CHEMIE» (Голландия). Импедансные данные моделировались эквивалентной электрической схемой (ЭЭС), границы раздела электрод-электролит с помощью компьютерной программы ZView 2.
Измерения осуществлялись в закрытых откаченных насосом до остаточного давления трёхэлек-тродных электрохимических ячейках в 7.6 M водном растворе КОН. В качестве противоэлектрода использовали неполяризованный цинковый электрод с свежозачищенной рабочей поверхностью 50 см2. Рабочий электрод и противоэлектрод разделялись сепаратором из нетканого полипропилена. В качестве электрода сравнения был использован хлоридсереб-ряный электрод, приготовленный на насыщенном KCl, а все электродные потенциалы пересчитыва-лись относительно стандартного водородного электрода. Электрод сравнения был соединён с ячейкой при помощи агар-агарового солевого мостика. Ёмкость ДЭС определялась по формуле
С = -( 2-jt-/-ImZ)
-1
(1)
где / — частота, 1т2 — мнимая составляющая импеданса. Значения ёмкостей для построения вольт-фарадных зависимостей брались на частоте 10"2 Гц.
Химические потенциалы ионов в растворе были определены из уравнения:
= [1® +RTlna±,
(2)
где а+ — средняя ионная активность электролита [11]. Они представлены в шкале абсолютных энергий.
РЕЗУЛЬТАТЫ И ИХ ОБСУЖДЕНИЕ
На основе экспериментальных данных метода малоуглового рассеяния рентгеновских лучей изучена пористая структура (рис. 1).
Для расчёта пористой структуры (таблица) были использованы известные соотношения:
тк =
24(0)/%3'
ГПк
Rgk
(3)
(4)
рх ' Рт
Здесь тк — объёмная доля пор к-й фракции, 4(0) — интенсивность рассеяния при 5 = 0, Rg¡c —
электронный радиус инерции пор, Б к — удельная поверхность пор к-й фракции (сферическое приближение), рх = 1.98 г/см3 — плотность углерода, рт = 0.82 г/см3 — плотность исходного материала.
К,, А
Рис. 1. Функция распределения пор по размерам для исследованного материала. Пунктирные кривые — гаусовское распределение пор соответствующего размера
Существенным отличием пористой структуры АДУ является смещение максимума функции распределения по размерам пор в область мезопор (13.7 А) по сравнению с активированным углем, полученным из стиролбензольного сополимера (АУС) с удельной поверхностью = 1073 м2/г, исследованным в работе [12], для которого этот максимум находится в микропористом диапазоне.
Построенные экспериментальные диаграммы Найквиста и зависимости Боде при электродных потенциалах 0.1—0.9 В (рис. 2) имеют типичный вид для доминирующего ёмкостного процесса на границе раздела, что выражается прямой под углом, близким к 90°, на диаграмме Найквиста и фазовым углом около -80° на частоте 10~2 Гц, что соответствует тангенсу диэлектрических потерь 0.18 и тепловым потерям мощности 0.17. По данным работы [3] значение псевдоёмкости может достигать 5-10% от общей ёмкости. Хорошо известная способность к квазиобратному электрохимическому восстановлению кислородсодержащих поверхностных групп в щелочных растворах вызывает увеличение ёмкости
вследствие их псевдоёмкостного характера проявляющегося только при достаточно низкой частоте 10~3 Гц [13], а потому влияние этих поверхностных групп на механизм образования двойного электрического слоя и кинетику процессов при заряде-разряде не было учтено.
гч 12 Г
0 2 4 6 8 10 12
Яе2, Ом-см2
а
б
Рис. 2. Типичные диаграмма Найквиста (а) и зависимость Боде (б) для древесного угля в 30%-ном КОН при Е = -0.9 В
Измеренную вольт-фарадную зависимость в потенциостатических условиях на частоте 10~2 Гц (рис. 3, а) можно разделить на три участка с различными геометрическими наклонами, соответствующие различным механизмам формирования ДЭС. Из приведённых зависимости видно, что катодная поляризация до Е = -0.9 В приводит к более чем двукратному росту максимальной ёмкости до С = 189 Ф/г, а положительная поляризация при £ = 0.08В — к уменьшению ёмкости до С = 58 Ф/г. Об уменьшении ёмкости при положительной поляризации свидетельствуют и литературные источники
8
Пористая структура АДУ О
____ Е
№ Rgk, А 4(0), у. Е. ть, части целого м2/г
1 13.7 771 0.722 875
2 36.9 3688 0.177 79
3 85.2 25822 0.101 20
Всего 1 974
[6]. Катодную часть (-0.11—0.9 В) вольт-фарадной зависимости (см. рис. 3, а) можно разбить на две части с разным наклоном, соответствующие различным механизмам адсорбции. Измеренный электродный потенциал неполяризованного электрода -—0.11 В близок к потенциалу плоских зон -0.05 В АУМ при рН = 12 [5]. А. И. Фрумкин [14] отмечал способность адсорбированного на поверхности активированного угля кислорода при контакте с водой переходить в раствор в виде ионов ОН", оставляя обеднённую электронами поверхность, что можно представить известным уравнением [15]:
СхО + Н20 + 2е~ = Сх+ 20Н", (5)
где СхО — поверхность АУМ с хемосорбированным кислородом, е~ —электрон, Сх — положительно заряженная поверхность АУМ.
Е, В
а
б
Рис. 3. Вольт-фарадная зависимость, полученная в потенцио-статических условиях на чатоте 10~2 Гц (а) и зависимость обратного квадрата ёмкости (б) для электрода из древесного угля в 7.6 М КОН
Состояние поверхности АУМ при потенциале неполяризованного электрода является результатом химического взаимодействия адсорбированного кислорода с водой согласно реакции (5). Уменьшение ёмкости на анодной ветви (-0.11—0.01 В) вольт-фарадной зависимости (см. рис. 3, а) связано с формированием на поверхности АУМ при данных условиях обеднённой полупроводниковой области Мотта—Шоттки [16]. Доказательством этого может служить прямолинейная зависимость вольт-фарадной кривой в координатах С~2—Е (рис. 3, б). Экстраполяция этой прямой к С-2 = 0 даёт положение уровня Ферми АДУ цс = -4.14 эВ. Учитывая рассчитанные химические потенциалы ионов №+ = -3.61 эВ и цон- = "2.87 эВ в 7.6 М КОН, а также соотношение для условия разблокирования поверхности для катодного и анодного процесса электро сорбции водорода и гидроксила, получаем равновесные потенциалы адсорбции Н+ и ОН" из следующего соотношения:
фк=(^с-^н+)/е и Фа = (v-он - \>с)/е, (6)
где фк, фа — равновесные электродные потенциалы адсорбции Н+ и ОН" соответственно. Рассчитанное значение электродного потенциала для электро сорбции ионов водорода по реакции:
Сх + Н++ е~ = СХН равно фк = -0,49В. (7)
А высокий энергетический уровень иона гидроксила, например в 6 М КОН цон- = -2.96 эВ, даёт значение фа = 1.18 В для процесса электро сорбции гидроксида по реакции:
Сх(ОН) + е" = Сх + ОН", (8)
В работах [5, 15] показано, что угольный электрод в щелочной среде при положительной поляризации функционирует в качестве кислородного газового электрода. А поэтому, по всей вероятности, процесс электро сорбции ОН" протекает с образованием молекулы воды по известной реакции (5). Протеканием процесса (5) обусловлено изменение наклона вольт-фарадной кривой (рис. 4, а)
при -0.01 В, а промежуток на ней -0.1--0.01 В
является перенапряжением этого процесса. Как видно на вольт-фарадной зависимости (см. рис. 3, а), равновесный потенциал электро сорбции водорода -0.49 В находится в начале второй её части, что может служить подтверждением вклада этого механизма в заряд ДЭС, который на катодной ветви сформирован гидратированными ионами калия.
На рис. 4 представлены частотные зависимости удельной ёмкости, которые кроме максимальных ёмкостей дают важную информацию о кинетике
процесса через время заряда-разряда / = /_1. Частотная зависимость при электродном потенциале неполяризованного электрода Е = -0.11 В показывает способность к самопроизвольному заряду пористой структуры С = 82 Ф/г, что происходит по электростатическому механизму. На частотных зависимостях ёмкости (рис. 4, а) можно выделить частотный промежуток 3-10 Гц, соответствующим временам цикла 0.1-0.33 с, где значение ёмкости почти постоянно — 28, 50, 73 Ф/г, и последующий частотный отрезок 0.01-3 Гц, где заряд достигает максимальных ёмкостей. Очень быстрый заряд-разряд (г = 0.1-0.33 с) связан с крупнейшими микро-и мезопорами. Заряд микропор меньших радиусов проходит значительно медленнее до г = 100 с, что ещё соответствует режимам работы ДСК. Как можно увидеть на рис. 4, а половину максимальной ёмкости 95 Ф/г можно получить на частоте 0.5 Гц, что соответствует времени цикла конденсатора 2 с. Частотные зависимости ёмкости (рис. 4, б) можно разбить на две основные области, которые определяются доминирующим процессом на поверхности электрода. Убывание частоты от 103 до 10~2 Гц приводит к росту ёмкости (см. рис. 4, б), а этот диапазон соответствует ёмкостному процессу с эффективным распределённым сопротивлением Я = 1.8 Ом-см2: на частоте 10~2 Гц оно равно 0.48 Ом-см2 , на частоте 103 Гц — 1.32 Ом-см2. А при частотах 103104 Гц не происходит существенного заряда пористой структуры (см. рис. 4, б), а растёт действительная составляющая импеданса. При этих частотах доминирующим является омический процесс, что соответствует общему вкладу сопротивлений электрода, объёма электролита и электролита в порах.
Характеристика омического процесса есть действительная составляющая импеданса на частоте 103 Гц или эффективное последовательное сопротивление С^эпс), равное 0.48 Ом-см2 при Е = -0.9 В.
Для установления более детальной взаимосвязи влияния пористой структуры на ёмкостные и кинетические параметры заряда-разряда, импедансные данные моделировались ЭЭС (рис. 5, а), учитывающей процессы внутренней релаксации пористой структуры в процессе заряда-разряда [17]. Это даёт возможность учитывать изменение внутренних характеристик пористой структуры, таких как внутреннее сопротивление и ёмкость, при работе ДСК. Порам определённого диаметра соответствует определённая ЛС-цепочка. Каждая ЛС-цепочка характеризуется своей постоянной времени (т = КС). Как было показано в работе [18], поры с близкими значениями постоянной времени при работе пористой структуры проявляют себя одинаково. Хорошо известно, что исследования энергетики адсорбции заставило дифференцировать класс микропор, выделив в нём узкие поры радиусом близким к 0.5 нм («микропоры»), в которых энергия адсорбционного взаимодействия с молекулами выше, чем энергия адсорбции в широких порах. А поры радиусом 0.5-1.0 нм («супермик-ропоры») по энергии адсорбционного взаимодействия и по характеру заполнения адсорбированными молекулами обладают переходными свойствами от микропор к мезопорам, приближаясь к тем или иным в зависимости от размера [19, 20]. Общепризнанным является то, что электрохимически доступными для заряда ионами в 30%-ном растворе КОН являются все поры с диаметром большим 0.5 нм [13]. Учитывая эти известные обстоятельства и используя
200
в
о"
150
100
200
2 5
_1_
J
9 10
/, Гц
в
о"
150
100
Рис. 4. Типичные частотные зависимости ёмкости древесного угля в 7.6 М КОН в диапазоне до 10 Гц (а): Е = -0.9 В (кривая 1), Е = -0.11 В (кривая 2) и Е = 0.08 В (кривая 3); и во всём частотном диапазоне (б) для Е = -0.9 В
линейную трансмиссионную модель [21] для моделирования работы пористой структуры в щелочном растворе, достаточно воспользоваться электрической цепью, состоящей из трёх параллельно соединённых ЯпСп -цепочек. Всем порам с диаметром большим 2 нм (мезо- и макропоры) отвечает первая К\С\-цепочка, ко второй цепочке (^2^2) относятся супер-микропоры, размер которых находится в пределах от 1 до 2 нм, а к третьей (Яз Сз -цепочке) — микропоры диапазона от 0.5 до 1 нм. Естественно, что размерный эффект иона будет влиять на максимальную ёмкость ДЭС, общепризнанным значением которой для отрицательной поляризации является 20 мкФ/см2. Соответственно по данным работы [22] ион калия может иметь две гидратные оболочки с радиусами 2.66 и 3.6 А первой и второй оболочек соответственно. А в работе [23] показана способность терять определённое количество молекул воды гидратной оболочки в процессе абсорбции в порах меньшего диаметра. Поэтому в микропорах ион калия имеет одну гидратную оболочку.
Рис. 5. Электрическая эквивалентная схема для моделирования работы пористой структуры (а) и вольт-фарадная зависимость (б), полученная как результат моделирования работы пористой структуры
Вольт-фарадная зависимость (рис. 5, б), построенная по определённым моделированным параметрам ЭЭС, показывает несколько меньшее значение максимальной общей ёмкости С = 175 Ф/г
по сравнению с данными, полученными на частоте Ю-2 Гц (см. рис. 3, а). Это, по всей вероятности, связано с влиянием электросорбции водорода, которая не учитывается моделью. Как хорошо видно на рис. 6, а основной вклад (до 142 Ф/г) в максимальную общую ёмкость делают мезопоры (^¡Сь рис. 6, кривые 1) и большие микропоры (ЯгСг, рис. 6, кривые 2) с малой постоянной времени т < 1 с. И наименее эффективная работа, с точки зрения ДСК, у меньших микропор ((ЯзСз) рис. 6, кривые 3).
в 80
о
70 60 50 40 30 20
А Ъ
а
J_1_I_1_I_1_I_iA- - А-А |_I
10
-1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2
Е, В
101
10и
10"
~о- -о-
-о-о---
-1.0 -0.S
-0.6 -0.4
-0.2
0.0
0.2
Е, В
Рис. 6. Вольт-фарадная зависимость для каждого типа пор (а) и Я„С„-время для каждого типа пор (б), полученные как результат моделирования работы пористой структуры исследованного материала электрической эквивалентной схемой
Необходимо отметить, что исследованный в работе [12] АУС имеет меньшую максимальную ёмкость 130 Ф/г по сравнению с исследованным в данной работе АДУ при значительно больших временах релаксации и Я, что может быть объяснено увеличением вклада мезопор для древесного угля и различным электронным строением материалов.
ВЫВОДЫ
1. Показано, что использование раствора 7.6 М КОН приводит к изменению химических потенциалов ионов, к расширению области потенциалов заряда двойного электрического слоя по чисто электростатическому механизму (-0.01 В—0.49 В) и электросорбцией водорода (-0.49—0.9 В). Положительная поляризация в области потенциалов -0.01-0.1 В приводит к электро сорбции гидроксильных групп с образованием молекул воды.
2. Рассчитанный потенциал разблокирования поверхности катодного процесса срк = -0.49 В вследствие электросорбции ионов водорода совпадает с началом изменения наклона вольт-фарадной зависимости.
3. Пористая структура исследованного материала позволяет быстро заряжать ёмкость до 140 Ф/г с постоянной времени т < 1 с. Время цикла заряд-разряд для заряда максимальной ёмкости 189 Ф/г составляет 100 с.
4. Мезопористость древесного угля позволяет значительно повысить ёмкостные показатели от 130 до 189 Ф/г и кинетические параметры конденсаторов двойного электрического слоя по сравнению с микропористым углеродным материалом с удельной поверхностью около 1000 м2/г.
СПИСОК ЛИТЕРАТУРЫ
1. Пат. 2800616 US (to General Electric), CI. 317-230, 1957.
2. Лидоренко H. С. II Докл. АН СССР. 1974. Т. 216. С.1261-1268.
3. Conway В. Е. Electrochemical Supercapacitors. New York: Plenum Publishing, 1999.
4. Волъфкович Ю. M., Сердюк Т. M. И Электрохимия. 2002. Т. 38. С. 1043-1068.
5. Бурштейн P. X., Вилинская В. С., Загудаева Н. М., Тарасевич М. Р. // Электрохимия. 1974. Т. 7. С. 1094-1097.
6. Belyakov A. I., Brintsev А. М., Khodyrevskaya N. И Proc. 14th Intern. Seminar on Double Layer Capacitors and Hybrid Energy Storage Devices. Deerfield Beach (USA), 2004. P. 84-94.
7. Endo M., Takeda Т., Kim Y J., Koshiba K, Ishii К. И Carbon Science. 2001. Vol. 1. P. 117-128.
8. Shi H. // Electrochim. Acta. 1996. Vol. 41. P. 1633-1639.
9. Kotz R„ Carlen M. II Electrochim. Acta. 2000. Vol. 45. P. 2483-2498.
10. Кют P. H., Сморгонская Э. А., Гордеев С. К, Гречинская А. В., Данишевский А. М. И ФТТ. 1999. Т. 41. С. 1484-1491.
11. Справочник по электрохимии / под ред. А. М. Сухотина. Л.: Химия. Ленингр. отд-ние, 1981.
12. Bakhmatyuk В. P., Venhryn B.Ya., Grygorchak I. I., Micov M. M., Kulyk Yu. О. И Electrochim. Acta. 2007. Vol. 52. P. 6604-6610.
13. Frackowiak E„ Beguin F. II Carbon. 2001. Vol. 39. P. 937-950.
14. Фрумкин A. H. Потенциалы нулевого заряда. M.: Наука, 1979. С. 104-105.
15. Пономаренко Е. А., Фрумкин А. Н., Бурштейн P. X. И Изв. АН СССР. 1963. № 9. С.1550-1555.
16. Гуревич Ю. Я., Плесков Ю. В. Фотоэлектрохимия полупроводников. М.: Наука, 1983.
17. Размерзотов В. Ю. И HaymBi BicTi НТУУ «КП1». 2009. № 3. С. 107-112.
18. Maletin Y A., Strizhakova N. G., Izotov V. Y., Mironova A. A., Danilin V. V. И New Promising Electrochemical Systems for Rechargeable. Batteries, 1996. P. 363-372.
19. Кулъский JI. А. Обработка воды на водопроводах пылевидным активированным углем. Киев.: Наук, думка, 1965.
20. Когановский А. М., Клименко Н. А., Левченко Т. М., Рода И. Г. Адсорбция органических веществ из воды. Л.: Химия. Ленингр. отд-ние, 1990.
21. Hyun-Kon Song, Hee-Young Hwang, Kung-Hong Lee, H. Le Dao И Electrochim. Acta. 2000. Vol. 45. P. 2241-2257.
22. Родникова M. H., Засыпкин С. А., Маленков Г. Г. И Докл. АН. 1992. Т. 324, № 2. С. 368-371.
23. Carritto-Trippa М., San-RomanbM. L., Hernacdez-Cobosc J., Saint-Martinc I H. И Biophysical Chemistry. 2006. Vol. 124, iss. 3, 1 December. 2006. P. 243-250.