ISSN 0868-5886
НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2020, том 30, № 4, с. 32-45 ФИЗИКА ПРИБОРОСТРОЕНИЯ =
УДК 541.182:537
© Б. П. Шарфарец, В. Е. Курочкин, 2020
ГИГАНТСКАЯ ДИСПЕРСИЯ ДИЭЛЕКТРИЧЕСКОЙ ПРОНИЦАЕМОСТИ ДИСПЕРСНОЙ СИСТЕМЫ В ПЕРЕМЕННОМ ЭЛЕКТРИЧЕСКОМ ПОЛЕ. ОБЗОР ПОДХОДОВ, УЧИТЫВАЮЩИХ НАЛИЧИЕ ДВОЙНОГО СЛОЯ
В обзоре проанализированы два подхода к описанию поляризации дисперсных систем при приложении, в общем случае, переменного электрического поля. Рассмотрены модели поляризации Трухана и Духина -Шилова. В результате проведенного обзора следует вывод о том, что наиболее предпочтительной моделью, позволяющей существенно нарастить скорость электрофореза дисперсной фазы, является модель Духина -Шилова, при реализации которой возможно появление гигантской дисперсии диэлектрической проницаемости гетерогенной дисперсной системы при выполнении условия наличия тонкого двойного слоя вокруг непроводящих дисперсных частиц в проводящей дисперсионной среде.
Кл. сл.: диэлектрическая проницаемость, гигантская дисперсия диэлектрической проницаемости, электрофорез, приближение тонкого двойного слоя, гетерогенная среда, поляризация Максвелла - Вагнера, дисперсная система
ВВЕДЕНИЕ
В работе [1] в интересах увеличения скорости электрофореза дисперсной фазы была рассмотрена возможность использования эффекта большой низкочастотной дисперсии относительной диэлектрической проницаемости в дисперсной системе, состоящей из дисперсной фазы в виде шарообразных частиц с параметрами е2,ст2, помещенных в однородную дисперсионную жидкость с параметрами е1,ст1, (е{ и а{ (/ = 1, 2) — относительные диэлектрические проницаемости и удельные проводимости соответствующих сред) в условиях поляризации Максвелла - Вагнера (ПМВ) [2, 3]. ПМВ носит название макроскопической, т.к. в ней игнорируется наличие двойного слоя на границах раздела фаз и не учитывается тем самым диффузный характер распределения заряда на границах раздела, что ведет к ограниченности применения этой теории областью тонкого двойного слоя ка ^ 1, где к = 1/ Л0, а — радиус включения, Л0
— радиус экранирования Дебая. В настоящей работе приводится обзор двух моделей теории поляризации дисперсных систем в условиях развитого диффузного распределения заряда двойного слоя: теории микроскопической поляризации Трухана [4-6] и Духина - Шилова [7].
Вначале выпишем наиболее полную систему уравнений, описывающих процесс поляризации в общем виде [7, § 3.1, § 3.2]:
су
р— = ^Р + ^Ау + РА, (1)
ОТ
V • у = 0, (2)
^±+0г = 0, (3)
ОТ
Ар = - ^^ (в СГСЭ). (4)
е
Система уравнений, описывающая рассматриваемый процесс, состоит из гидродинамического блока, включающего модернизированное уравнение Навье - Стокса для малых чисел Рейнольдса (1) (называемое в этом случае уравнением Стокса) и уравнение непрерывности для несжимаемой жидкости (2). Уравнение (3) является уравнением непрерывности потока ионов (см., например, [7, с. 76]). В качестве вектора рассматривается плотность потока этих ионов [7, с. 76]:
Г = с~ушр - Е>^ с~ + с~у =
= ±с V Е - D ^ • с± + с± у. (5)
Наконец уравнение Пуассона (4) описывает поведение электрического потенциала р (—Vp = E; описание вектора E см. ниже) при наличии распределения плотности электрического заряда ре.
В уравнениях (1)-(4), кроме отмеченных, приняты следующие обозначения: р = const — плотность несжимаемой жидкости; v — вектор поля скорости жидкости; p — поле давления в жидкости; т] — динамическая вязкость жидкости; вектор электрической напряженности внешнего поля E0 ориентирован, например, по оси Oz E0 = (0,0, Ez) и имеет гармонический характер E0 = kEze~iat, где k — орт вдоль оси Oz; s — относительная диэлектрическая проницаемость; D — коэффициент диффузии иона; c — коэффициент концентрации; / — подвижность иона. В работах [4-7] функция источника реЕ0 в уравнении (1) отсутствует, взамен этого ставится граничное условие о том, что значение электрического поля в системе Е на бесконечности стремится к значению Е0: Е0 = lim Е, x — вектор текущей
|x| ^да
точки в исследуемой области; именно значение поля Е фигурирует в (5). Обусловлено это тем, что в среде присутствует также поляризованное электрическое поле, обусловленное потенциалом р из (4). Согласно [6, с. 8], поле Е имеет вид:
Е ( x, t ) = Е (x ) e~" р( x, t ) = р( x ) e~ic Е ( x ) = -Vp( x ).
(6)
Далее переходим к конкретному рассмотрению заявленных видов поляризации.
ПОЛЯРИЗАЦИЯ ТРУХАНА
Раздел написан по материалам работ [4-8], в которых присутствует более обширная библиография по рассматриваемому вопросу. В наиболее концентрированном виде материал представлен в [6]. Необходимо отметить, что работе Э.М. Тру-хана [4] предшествовала изданная значительно раньше работа Я.И. Френкеля [8], содержавшая, однако, подмеченные в [4] неточности применительно к результатам, касающимся поляризации в переменном электрическом поле. Поэтому именно работу [4] можно назвать первой корректной работой, посвященной микроскопической поляризации дисперсных систем в переменном электрическом поле с учетом наличия двойного слоя.
Автором теории была исследована система сферических проводящих включений в непроводя-
щей среде. Задача рассматривалась в предположении, что период колебаний T = 2ж]со, среднее столкновительное время релаксации подвижных носителей т1 и среднее время рекомбинации зарядов противоположного знака в проводящей фазе т2 связаны соотношением т1 > T > т2; это ограничивает область рассмотрения частот диапазоном частот f = С(2я): 103-106 < / < 10п-1012 Гц [6, с. 8]. Для взвеси проводящих сфер в изоляторе
точное выражение для комплексной совокупной
^ *
диэлектрической проницаемости е получается на основе соотношения е = б / Е, где б и Е — средние комплексные амплитуды электрической индукции и напряженности поля соответственно. Взаимодействием системы сфер пренебрегается, т.е. принимается, что величина удельного объема а шаровых включений а = (4/3) яа3N (а — радиус шаровых включений; N — число включений в единице объема или их концентрация) мала: а< 1.
Изначально принимается условие малости чис-Va
ла Пекле Ре = — < 1 (здесь V — характерный
размер скорости; а — характерный размер дисперсной частицы), что позволяет опустить конвекционный член с±V в (5) и свести систему (1)-(5) к двум несвязанным системам: гидродинамической (1), (2) и электрохимической (3)-(5).
Здесь приведем решение задачи (3) в предположении, что б += б"= б и |и+| = |и"| = И Концентрация ищется в виде
с±( X, t ) = (с0 + у+-( X )) (7)
при условии слабой поляризации
Г" << со. (8)
Здесь с0 — равновесная концентрация ионов (дисперсных частиц); у± — амплитуда отклонения концентраций от равновесного значения, которая полагается много меньшей равновесной концентрации.
Предположение (7), (8) позволяет линеаризовать задачу и получить для плотности заряда уравнение
Аре (X) = Х2ре (X) , (9)
где
2 2 Ш 1 / т
X = к +—, к = 1/Av
Л р, -с
(10)
С учетом условия конечности плотности заряда в центре сферы решение (9) имеет вид
Ре(г)=
А ^6
(
ch ( г )-
sh ( г )
(11)
1 - ^ + р п е
2 (
1 + 2е1 - 2р
—е а Е
2 «1" ^0'
(12)
ст = астп
V е2 У
3" г,
г = -
соеп
(15)
где г = отсчитывается от центра сферы. После
подстановки (11) в уравнение Пуассона автор теории поляризации получает уравнение для определения потенциала р, которое существенно отличается от уравнения (9). Далее автор находит распределение потенциала р в окрестностях сферы. На основании этого им получается выражение для дипольного момента Се* (дипольный момент позволяет представить дипольную составляющую рр электрического потенциала р сферы во внешней области рп = Се* соъб/г2) [6, с. 11]:
и монотонно растет с частотой и достигает предельного значения ст
ст = 9аст,
1+2 е
Л
(16)
'2 У
Физический смысл расхождения данной теории с теорией поляризации ПМВ заключается в том, что в теории ПМВ пренебрегается наличием диффузного заряда на границе раздела фаз, т.е. заряд полагается чисто поверхностным. Это справедливо только при условии ка ^ 1. В общем случае считать заряд поверхностным (ПМВ) неверно. В работах Я.И. Френкеля [8] и Э.М. Трухана [4 и др.] произведено обобщение теории ПМВ на случай произвольного соотношения между величиной дебаевского радиуса экранирования и размерами включений.
где
Р = -
( 3 + (^а )2) Л (Ха )
- 3Ха
е2^
2+(ха )2^ ( Ха )- 2^а
(13)
Для системы сферических включений при их удельном объеме а ^ 1 получено следующее выражение для результирующей диэлектрической проницаемости гетерогенной системы:
е = е
[1 + 3а (1 - е1 / е2 + р)]
(14)
1 + 2е1 / е2 - 2р где Р определяется из (13). Характер частотной зависимости е' = Re (е ) похож на соответствующую зависимость теории поляризации ПМВ.
Диэлектрические свойства такой системы характеризуются удельным инкрементом 3
з=
=з'_ 3-,
а
который, как видно из формулы, равен изменению диэлектрической проницаемости дисперсионной среды за счет присутствия дисперсной фазы на единицу удельного объема и на единицу исходной диэлектрической проницаемости дисперсионной среды.
Удельная проводимость среды ст при одинаковых подвижностях /+ = / = / равна
ПОЛЯРИЗАЦИЯ ДУХИНА - ШИЛОВА
Авторами теории, в отличие от модели, рассмотренной выше, была исследована система различных непроводящих включений в проводящей среде. Важным свойством границы раздела включений и дисперсионной жидкости, влияющим на поляризацию дисперсных систем и не учтенным ни в теории ПМВ, ни в работах [4] и [8], является заряженность поверхности раздела. Заряд на поверхности раздела обуславливает наличие двойного слоя и связанной с ним дополнительной поверхностной проводимости, которая влияет на распределение потенциала, возникающего во внешнем электрическом поле [7, с. 7]. Исследование влияния диффузионного механизма на поляризацию в предельном случае ка ^ 1 привело авторов к выводу о низкочастотном отставании поля поляризации от внешнего электрического поля, что и является причиной наблюдаемых на опыте гигантских значений низкочастотной диэлектрической проницаемости дисперсных систем [7, с. 8]. В то же время область рассмотрения теории поляризации Трухана начинается с частоты f > 103Гц (см. выше).
В [7, с. 42] приводится факт того, что наличие удельной поверхностной проводимости Л в области двойного слоя так же влияет на поляризационные поля сферических частиц радиусом а, как и увеличение их объемной проводимости ст на величину 2 Л а, т.е. результирующая объемная проводимость ст при наличии поверхностной про-
г
г
2
е1к
2
е1а
водимости равна с? = и + 2//а. В переменном поле возникает дисперсия поверхностной проводимости из-за того, что больцмановское равновесие в двойном слое уже не успевает восстанавливаться, что и приводит к дисперсии поверхностной проводимости. Отмечается также [7, с. 43], что наличие у частицы дисперсной фазы большого постоянного дипольного момента также скажется на диэлектрических свойствах суспензий и растворов макромолекул.
Большие размеры дисперсных частиц по сравнению с молекулами полярного растворителя обуславливают дисперсию диэлектрической проницаемости на низких частотах [7, с. 46]:
= ен +(ео -ен)
1
1 + (тс)
3Г/У ^
где т = —--время релаксации; Ур — объем
частицы; ^ — число Авогадро. Там же приведено выражение для дисперсии диэлектрической проницаемости и для частиц несферической формы. В этом случае диэлектрическая проницаемость характеризуется набором времен релаксации [7, с. 46]. Далее в [7, с. 51] приводится очень важное замечание о том, что в дисперсных системах с проводящей дисперсионной средой возможен механизм усиления влияния поляризации на действительную часть диэлектрической проницаемостиЭтот механизм объясняется в [7] следующим образом. При подаче переменного поля к дисперсной среде возникает процесс поляризации, который отстает от исходного электрического поля, что приводит к отставанию тока проводимости, вызванного поляризационным электрическим полем. Этот ток в силу его отставания от внешнего поля интерпретируется как ток смещения. Наличие же значительных токов смещения в среде на низких частотах воспринимается экспериментально как высокая диэлектрическая проницаемость
среды (см., например, выражение (1.54) работы [7]
*
для диэлектрической проницаемости е неиде-
* 4пи
ального конденсатора: е = е+--, где е — ди-
ш
электрическая проницаемость диэлектрика в конденсаторе, и — его удельная проводимость).
Далее в [7, с. 52, 53] приводится важное условие появления высоких значений диэлектрической проницаемости суспензий, связанных с поляризацией двойного слоя частиц:
ка » 1.
Наряду с этим условием для возникновения гигантской низкочастотной дисперсии необходимы немалые значения величин поверхностного заряда. Далее авторы отмечают, что только при условии ка ^ 1 действует проанализированный эффект гигантского усиления диэлектрических явлений [7, с. 55].
Рассмотренный в [7] механизм гигантского усиления влияния поляризации включений на диэлектрическую проницаемость дисперсной системы, по мнению авторов работы, должен иметь место для более широкого класса систем, чем суспензии и электролиты, для которых построена соответствующая теория в [7]. Здесь, видимо, уместно повторить, что речь идет действительно о гигантских значениях диэлектрической проницаемости, которая, например, в водных суспензиях достигает размеров порядка е~ 104. Такая низкочастотная дисперсия диэлектрической проницаемости наблюдалась для проводящих суспензий и растворов полиэлектролитов [7, с. 55].
Характерные признаки дисперсии диэлектрической проницаемости
- Очень высокие значения низкочастотного предела, достигающие 105.
- Большая величина времени релаксации, квадратичная зависимость времени релаксации от линейных размеров частиц. Оценка времени релаксации может быть сделана по формуле
т ~ 12 / б ;
где I — характерный размер частицы (макроиона); б — коэффициент его диффузии в растворе электролита.
- Повышенная проводимость противоионов двойного слоя.
Согласно [7, с. 56], обнаружена также зависимость низкочастотной диэлектрической проницаемости от величины рН.
Учет перечисленных условий позволяет объяснить высокие значения диэлектрической проницаемости на основе существующих представлений о строении диффузной части двойного слоя.
При разработке рассматриваемой теории за основу также берется система уравнений (1)-(5). Распределение потенциала и и концентрации с± в дисперсной среде ищется в виде [7, с. 63]
Р = Ро +фо + и1,
С±= со + Го±+ С . (17)
1) Напомним, что в методе Френкеля - Трухана дисперсионная среда принята непроводящей.
Здесь Фо — равновесное распределение электрического потенциала в двойном слое частицы в отсутствие внешнего электрического поля; ро — распределение электрического потенциала под
£
воздействием внешнего поля; и1 — линейная по внешнему полю поляризационная добавка (принимается предположение о достаточно малой величине внешнего поля для справедливости достаточности только добавки и1); с± — поляризационное распределение концентрации; с± — концентрации ионов в свободном электролите; у\ — отклонение концентраций с± в окрестностях двойного слоя; с± — линейные по внешнему полю поляризационные добавки к концентрациям.
Стационарная поляризация
Далее авторами рассматривается двойной слой сферической частицы в постоянном электрическом поле [7, § 3.1] (в этом случае уравнение непрерывности (3) сводится к стационарному случаю У]± = о). При этом принимается допущение о приближении тонкого двойного слоя (16), что позволяет рассматривать его как локально плоский. Поляризационный потенциал и концентрации ищутся в виде
и = р +Ф1,
±±± С1± = с1± + тт.
В результате непростых выкладок в [7] получены следующие результаты.
- Распределение Больцмана сохраняется для концентрации ионов и для поляризованного тонкого двойного слоя с тем отличием, что к обычному распределению Больцмана распределения концентрации ионов добавляется зависящий от полярного угла в предэкспоненциальный множитель (в = о соответствует направлению постоянного электрического поля), т.к. концентрация ионов изменяется вдоль внешней границы двойного слоя [7, с. 7о], и, следовательно, структура поляризованного двойного слоя совпадает со структурой равновесного двойного слоя (в отсутствие внешнего электрического поля).
- Напряженность поля, возникающая внутри двойного слоя при его поляризации, в ка раз превышает напряженность поля у его внешней границы; в подобной ситуации выполнение условий непрерывности потоков ионов внутри двойного слоя возможно лишь за счет почти полной компенсации нормальных компонент электромиграционных и диффузионных потоков ионов, что и приводит к необходимости выполнения распределения Больцмана [7, с. 75].
- При слабой поляризации двойного слоя (т.е. при выполнении условий с± ^ у± + С1± из (17)) автоматически следует выполнение условия мало-
уа
сти критерия Пекле Ре = — < 1 (см. выше), что
позволяет свести систему (1)-(5) к двум несвязанным системам: гидродинамической (1), (2) и электрохимической (3)-(5) [7, с. 76].
Поляризация в переменном электрическом поле
Далее авторы обобщают полученную для стационарного режима теорию на случай переменных полей не слишком высокой частоты, сохраняя при этом упрощающие факторы, учтенные в теории стационарной поляризации. За основу вновь берется система уравнений (1)-(5), однако с уравнением непрерывности в неусеченном виде (3). При этом вновь принимается допущение относительно малого перепада концентрации, ведущего к малости числа Пекле и, следовательно, к независимости уравнений (1), (2) и (3), (4).
Далее авторы находят ограничение по частоте внешнего электрического поля (выражение (III.88) [7, с. 77]), при котором сохраняется локальное равновесие в двойном слое и при поляризации в переменном электрическом поле для концентраций ионов можно использовать в двойном слое распределение Больцмана [7, с. 78] (квазистационарное решение). Вследствие малых чисел Пекле рассматриваемой задачи авторы заменяют условие приближенной локальной электронейтральности за пределами двойного слоя условием полной электронейтральности. На этом основании вводится функция с = с+/г"= с"/г+ , для которой справедлив анзац [7, с. 77]
с = со + с1 (г ,0) е а для потенциала вне частиц: р = р(г,в)е-с.
(18)
(19)
Условие квазистационарности выполняется при выполнении условия [7, с. 78]:
2ж 1
Т = — >—-, или ю^2пук2 , (2о)
С УК
где у — кинематическая вязкость дисперсионной среды.
Решение для р и с ищется в виде
с = со + dccr (г)cosвe
(21)
Р = 1Ро + dc
б' - б+
г+б+ + г " б-
. (г)С08в +
г2
+-рС08в \е-с,
Ро =- Е
( 1 а ^
Г +--;
2 Г V ^ ' /
С0%в ,
(22) (23)
где Е — амплитуда гармонического электрического поля Е = |Е|; Сс, Ср и сг определены в [7, с. 80].
После подстановки Ср в формулу (11.39) работы [7, с. 52] с точностью до величины (ка) 2 получается следующее выражение для невязки между диэлектрическими проницаемостями поляризованной дисперсной смеси и дисперсионной жидкости (т.е. диэлектрического инкремента Ае = е* -е1) [7, с. 81]:
Ае =
рассматриваемой дисперсной среды (s" = Im (V )) определяется выражением [7, с. 83]
„9/42 s = — as1 (ка ) х
(А1а2 - A,a1) (A1W + A, W2 )2 A1 (A + A1W )2 + (A1W + AW2)
(27)
9 . ч 2 (Да2 -M )(l + W + W2) = — as1 (ка) ---
4 ( a + aw )2 +(AW+AW2)
(24)
Из (27) видно, что величина s" мала на очень низких частотах при W ^ 1 и на высоких частотах при W ^ 1, достигая максимального значения при W ~ 1, при этом, согласно зависимости (27), частота, на которой достигается максимум диэлектрических потерь, т.е. при s" ^ max будет обрат-
Величина статического диэлектрического ин- но пропорциональна квадрату радиуса частицы
кремента Asst = lim As равна ) [7, с. 81]
9 , \ 2 Аа - A а Asst = — as1 (ка) 12 2 2 1 .
st 4 1V ' A,2
(25)
Здесь величины А1, А2, а1 и а2 определены в [7, с. 80], W = а/г0, где г0 = — характеризу-
ет толщину диффузного слоя при частоте с; Бе = (г+ + г ")/(г+ + г-D") [7, с 77]. Далее понадобится зависимость
w=a=
Л
а 2ф
42DJ ф j2De'
(26)
На частотах, когда W ~ 1, поляризация частицы существенно отстает по фазе от приложенного поля и действительная часть дипольного (поляризационного) момента изменяется наиболее быстро. Таким образом, медленный процесс установления распределения концентраций за пределами двойного слоя приводит к существенному отставанию поляризации частицы от внешнего De и
поля на частотах с--— . На этих частотах возни-
а
кают высокие значения при малой толщине двойного слоя, когда ка ^ 1 [7, с. 81].
Для симметричного электролита с точностью до величины (ка) 2 величина мнимой части е" комплексной диэлектрической проницаемости е
На рисунке приведена зависимость величины удельного инкремента реальной части диэлектрической проницаемости Ае / а от величины Ж,
2) Величина Аег следует из (24) тривиально по причине того, что Ж (с = 0) = 0. См. зависимость (26) для Ж (с).
[7, рис. 11]: зависимость инкремента реальной части диэлектрической проницаемости Ае / а (1, 2) и удельной величины мнимой части диэлектрической проницаемости е"/а (1', 2') от логарифма величины Ж для водных суспензий сфер при ка = 50. Кривые 1 и 1' получены при значении ^-потенциала в полтора раза большем, чем в случае кривых 2 и 2'
х
2
а
а
r
0
а фактически от л/С в логарифмическом масштабе, по оси абсцисс (величина а — по-прежнему величина удельного объема шаровых включений в дисперсионной жидкости). Заимствован рисунок 11 работы [7, с. 81]. Одновременно на рисунке приведены зависимости величины е"/ а . Из рисунка виден характер поведения величин Ае / а и е ""/ а , а также прямая зависимость этих величин от величины ^-потенциала.
ВЫВОДЫ
В результате проведенного обзора следует вывод о том, что при выполнении условия наличия тонкого двойного слоя вокруг непроводящих дисперсных частиц в проводящей дисперсионной среде возникает возможность появления гигантской дисперсии действительной составляющей диэлектрической проницаемости, которая может существенно изменить скорость электрофореза таких частиц.
Работа выполнена в ИАП РАН в рамках Государственного задания 075-00780-20-00 по теме № 0074-20190013 Министерства науки и высшего образования.
СПИСОК ЛИТЕРАТУРЫ
1. Шарфарец Б.П., Курочкин В.Е. О возможности использования гигантской дисперсии диэлектрической проницаемости дисперсной системы для увеличения скорости электрофореза дисперсных частиц // Научное приборостроение. 2о2о. Т. 3о, № 3. С. 38-44.
2. Максвелл Дж.К. Трактат об электричестве и маг-
нетизме. Т. 1. М.: Наука, 1989. 416 с.
3. Wagner K.W. Electricity of the dielectric behavior on the basis on the Maxwell theory // Arch. Elektrotech. 1914. No. 9. P. 371-392.
4. Трухан Э.М. Дисперсия диэлектрической проницаемости гетерогенных систем // Физ. тв. тела. 1962. Т. 4, № 12. С. 3496-3511.
5. Трухан Э.М. Электронный транспорт в субклеточных структурах и их моделях. Разработка радиофизических методов исследования и квантово-электронный анализ. Дис. ... д-ра. физ-мат. наук. МГУ, 1980.
6. Трухан Э.М Электрофизические методы исследования. Кондуктометрия неоднородных материалов: Учебно-методическое пособие. М.: МФТИ. 26 с.
7. Духин С.С., Шилов В.Н. Диэлектрические явления и двойной слой в дисперсных системах и полиэлектролитах. Киев: Наукова думка, 1972. 206 с.
8. Френкель Я.И. Электрические колебания в дисперсных системах // Коллоидн. журн. 1948. Т. 10, № 2. С. 148.
Институт аналитического приборостроения РАН, Санкт-Петербург
Контакты: Шарфарец Борис Пинкусович sharb@mail.ru
Материал поступил в редакцию 28.08.2020
ISSN0868-5886 NAUCHNOE PRIBOROSTROENIE, 2020, Vol. 30, No. 4, pp. 32-45
GIANT DISPERSION OF DIELECTRIC PERMEABILITY OF THE DISPERSE SYSTEM IN AN ALTERNATING ELECTRIC FIELD.
OVERVIEW OF APPROACHES TAKING INTO ACCOUNT PRESENCE OF A DOUBLE LAYER
V. E. Kurochkin, B. P. Sharfarets
Institute for Analytical Instrumentation of RAS, Saint Petersburg, Russia
This review analyzes two approaches to describing the polarization of disperse systems when an alternating electric field is applied in the general case. Polarization models of Trukhan and Dukhin - Shilov are considered. As a result of this review, it follows that the most preferable model, which makes it possible to significantly increase the rate of electrophoresis of the dispersed phase, is the Dukhin - Shilov model, in the implementation of which a giant dispersion of the dielectric constant of a heterogeneous dispersed system may appear if the condition of a thin double layer around non-conducting dispersed particles in a conducting dispersive medium.
Keywords: the dielectric constant, giant dispersion of the dielectric constant, electrophoresis, thin double layer approximation, heterogeneous environment, Maxwell—Wagner polarization, dispersion system
INTRODUCTION
In [1], for the aim of increasing the rate of electro-phoresis of the dispersed phase, the possibility of using the effect of large low-frequency dispersion of the relative dielectric permittivity in a disperse system is considered. The disperse system consists of a dispersed phase in the form of spherical particles with parameters s2,a2, placed in a homogeneous dispersion liquid with parameters s1,a1, (st h at (i = 1,2) that are the relative dielectric permittivity and the specific conductivity of the corresponding medium) in terms of the Maxwell-Wagner polarization (MWP) [2, 3]. MWP is called macroscopic because it ignores the presence of a double layer at the interfaces and thus does not take into account the diffuse nature of the charge distribution at the interfaces, which leads to the limited application of this theory to the region of a
thin double layer tea » 1, where k = 1 / XD
is the
Ov
p— = -Vp + rAv + PeE0: ot
V-v = 0,
Vj± +°C± = o, ot
Ap = -(in CGS).
(2)
(3)
(4)
inclusion radius, XD is the Debye length. This paper provides an overview of two models of the theory of polarization of disperse systems under conditions of a developed diffuse distribution of the charge of a double layer: the theories of microscopic polarization of Trukhan [4-6] and Dukhin-Shilov [7].
First, we write out the most complete system of equations describing the polarization process in general form [7, § 3.1 h § 3.2].
(1)
The system of equations describing the process under consideration consists of a hydrodynamic unit, which includes the modernized Navier-Stokes equation for small Reynolds numbers (1) (called the Stokes equation in this case) and the continuity equation for an incompressible fluid (2). Equation (3) describes the continuity of the ion flux (see, for example, [7, p. 76]). The flux density of these ions is considered as a vector j± [7, p. 76]
i = c~ v . - D~ V- c~ + c~ v =
J migr
= ±c±^±E - D±V - c ±+ c± v. (5)
Finally, the Poisson equation (4) describes the pattern of the electric potential p (-Vp = E; for a description of the vector E see below) in the presence of a distribution of the electric charge density pe.
In equations (1)-(4), in addition to those noted, the following designations are adopted: p = const — density of an incompressible fluid; v — the vector of the fluid velocity field; p — pressure field in the liquid; j — dynamic viscosity of the liquid; the vector of the electric intensity of the external field E0 is oriented,
s
a
for example, along the axis Oz E0 = (0,0, Ez) and has a harmonicity E0 = kEze~—, where k — the unit vector along the axis Oz; s — relative dielectric constant; D — ion diffusion coefficient; c — concentration factor; ¡i — ion mobility. In [4-7], the source function peE0 in equation (1) is absent; instead, the boundary condition is set that the value of the electric field in the system E at infinity tends to the value E0 : E0 = lim E , x is the vector of the current point
|x| —><X)
in the investigated area; it is the field E value that appears in (5). This is due to the fact that there is also a polarized electric field in the medium due to the potential p from (4). According to [6, p. 8], the field E has the form
E ( x, t ) = E ( x ) e~ia <p( x, t ) = <p( x ) e~ia E ( x ) = -Vp( x ).
(6)
Next, we turn to a specific consideration of the declared types of polarization.
THE TRUKHAN POLARIZATION
The section was written on the basis of the works [4-8], which contain a more extensive bibliography on the issue under consideration. The most concentrated material is presented in [6]. It should be noted that the work of EM. Trukhan [4] was preceded by the work of Ya.I. Frenkel [8], which, however, contained the inaccuracies noted in [4] as applied to the results concerning polarization in an alternating electric field. Therefore, it is work [4] that can be called the first correct work devoted to the microscopic polarization of disperse systems in an alternating electric field with due regard for the presence of a double layer.
The author of the theory investigated a system of spherical conducting inclusions in a non-conducting medium. The problem was considered under the as-
2^
sumption that the period of oscillations T = —, the
a
average collisional time of mobile carriers relaxation t1 and the average time of recombination of charges of opposite sign in the conducting phase t2 are related as t1 > T > t2 ; this limits the area of consideration of
frequencies to the frequency range f = — :
2n
103-106 < f < J0jj-1012 Hz [6, c. 8]. For a suspension of conducting spheres in an insulator, the exact expression for the complex total dielectric constant s* is obtained on the basis of the rela-
tion s* = D / E, where D and E are the average complex amplitudes of the electric induction and field strength, respectively. The interaction of the system of spheres is neglected, i.e. it is assumed that the value of the specific volume a of spherical inclusions 4n 3
a =— aN (a — radius of spherical inclusions;
N — number of inclusions per unit volume or their concentration) is small: a ^ 1.
Initially, the condition for the smallness of the Pec-
Va
let number Pe = — <C 1 is accepted (here V is the
characteristic size of the velocity; a is the characteristic size of the dispersed particle), which makes it possible to omit the convection term c±v in (5) and reduce system (1)-(5) to two unconnected systems: hy-drodynamic (1), (2) and electrochemical (3)-(5). Here we present the solution to problem (3) under
the assumption that D+ = D~ = D and 1 = | = X .
The concentration is sought in the form
c±(xt) = (Co + r (x))
under the condition of weak polarization
+ ^
7~ < Co.
(7)
(8)
Here c0 is the equilibrium concentration of ions (dispersed particles); y± — the amplitude of the deviation of concentrations from the equilibrium value, which is assumed to be much less than the equilibrium concentration.
Assumption of (7)-(8) allows linearizing the problem and obtaining the equation for the charge density
APe (x) = X2Pe (x),
where
ia
X = * +—, k = 1/V
(9)
(10)
Taking into account the condition of finiteness of the charge density at the center of the sphere, solution (9) has the form
Pe ( r ) =
A cosd
(
ch ( r )-
sh ( r )
(11)
where r = |x| is measured from the center of the
sphere. After substituting (11) into the Poisson equation, the author of the polarization theory obtains an equation for determining the potential p, which differs significantly from equation (9). Further, the author finds the distribution of the potential p in the
r
r
vicinity of the sphere. Based on this, he obtains an expression for the dipole moment d* (the dipole moment allows one to represent the dipole component (pD of the electric potential of the sphere in the outer
region
(D = de* cos#/r2 ) [6, p. 11]
d* =
p.
1 - ^ + p ж s
2 ( 2s 1 + ^ - 2p
Y2
-sa E„
(12)
theory neglects the presence of a diffuse charge at the interface, i.e. the charge is assumed to be purely superficial. This is true only if Ka » 1. In the general
case, it is incorrect (MWP) to consider the charge surface. In the works of Ya.I. Frenkel [8] and E.M. Trukhan [4, et al.] generalization of the MWP theory for the case of an arbitrary relationship between the value of the Debye length and the size of the inclusions is carried out.
THE DUKHIN - SHILOV POLARIZATION
where
P =
(3 + (Xa )2 ) th ( xa )
- 3Xa
S2^
2 + (xa )2th (xa )- 2Xa
(13)
For a system of spherical inclusions with their specific volume a □ 1, the following expression was obtained for the resulting dielectric constant of the heterogeneous system
s = s
[1 + 3a (1 - sj s2 + p)] 1 + 2s1 / s2 - 2P
(14)
where ( is obtained from (13). The nature of the frequency dependence s' = Re (s*) is similar to the corresponding dependence of the MWP theory.
The dielectric properties of such a system are characterized by the specific increment 8
* * 5 = s_=s!
s 's1 -1 = 5'-/5" a
which, as seen from the formula, is equal to the change in the dielectric permittivity of the dispersion medium due to the presence of the dispersed phase per unit of specific volume and per unit of the initial dielectric permittivity of the dispersion medium.
The specific conductivity of the medium a at the same mobility /+ = ¡u- = /u is
f r \
7 = ao„
V s2 J
S"t.
t = -
asn
(15)
and monotonically increases with frequency and reaches the limiting value aw
< = 9a<.
( ~ V
( s > 1 + 2 s
s2 J
(16)
The physical meaning of the discrepancy between this theory and the MWP theory is that the MWP
The authors of the theory, unlike the model considered above, investigated a system of various nonconducting inclusions in a conducting medium. An important property of the interface between inclusions and a dispersion liquid, which affects the polarization of dispersed systems and is taken into account neither in the MWP theory nor in the works [4] and [8 ], is the charge of the interface. The charge on the interface determines the presence of a double layer and the associated additional surface conductivity, which affects the distribution of the potential arising in an external electric field [7, p. 7]. Investigation of the influence of the diffusion mechanism on polarization in the limiting case ta »1 led the authors to the conclusion
about the low-frequency lag of the polarization field from the external electric field, which is the reason for the experimentally observed giant values of the low-frequency dielectric constant of dispersed systems [7, p. 8]. At the same time, the field of consideration of the Trukhan polarization theory begins with the frequency f > 103 Hz (see above).
In [7, p. 42], the fact is given that the presence of specific surface conductivity X in the double layer region also affects the polarization fields of spherical particles with a radius a, as well as an increase in
2X
their bulk conductivity a by an amount —, i.e. the
a
resulting bulk conductivity a in the presence of surface conductivity is a = a + — . In an alternating
a
field, dispersion of surface conductivity arises due to the fact that the Boltzmann equilibrium in the double layer has no time to be restored, which leads to dispersion of surface conductivity. It is also noted [7, p. 43] that the presence of a large constant dipole moment of a particle of a dispersed phase will also affect the dielectric properties of suspensions and solutions of macromolecules.
Large sizes of dispersed particles in comparison with molecules of a polar solvent cause the dispersion of the dielectric constant at low frequencies [7, p. 46]
2
s1k
2
1
= sh + (s0 "sh)
1 + (ra)
where r =
wna
----the relaxation time; V — par-
RT p F
ticle volume; NA — the Avogadro's number. The expression for the dispersion of the dielectric constant and for particles of non-spherical shape is also given there. In this case, the dielectric constant is characterized by a set of relaxation times [7, p. 46]. Further in [7, p. 51], a very important remark is made that in disperse systems with a conductive dispersive medium, a mechanism for enhancing the effect of polarization on the real part of the dielectric constant is possi-ble:).This mechanism is explained in [7] as follows. When an alternating field is applied to a dispersed medium, a polarization process occurs, which lags behind the initial electric field, which leads to a lag in the conduction current caused by the polarization electric field. This current, due to its lagging behind the external field, is interpreted as a displacement current. The presence of significant displacement currents in the medium at low frequencies is perceived experimentally as a high dielectric constant of the medium (see, for example, expression (1.54) in [7] for the dielectric constant of a non-ideal capacitors*:
s* = s + ^n<J , where s is the dielectric constant of ia
the dielectric in the capacitor, < is its specific conductivity).
Further in [7, p. 52, 53], an important condition for the appearance of high dielectric permittivities of suspensions associated with polarization of the double layer of particles is given
KÚ > 1.
Along with condition, for the appearance of a giant low-frequency dispersion, considerable values of the surface charge are required. Further, the authors note that only under condition kú ^ 1 the analyzed effect
of giant amplification of dielectric phenomena acts [7, p. 55].
The mechanism of the giant enhancement of the influence of polarization of inclusions on the dielectric permittivity of a dispersed system considered in [7], according to the authors of this work, should be applied for a wider class of systems than suspensions and electrolytes, for which the corresponding theory was constructed in [7]. Here, apparently, it is appropriate to repeat that the object of consideration is in-fact gigantic values of the dielectric constant, which,
for example, in aqueous suspensions reaches values of the order of s ~ 104. Such a low-frequency dispersion of the dielectric constant was observed for conducting suspensions and solutions of polyelectrolytes [7, p. 55].
The characteristic features of the dispersion of the dielectric constant
• Very high values of the low-frequency limit, reaching 105;
• Large value of the relaxation time, quadratic dependence of the relaxation time on the linear dimensions of the particles. The relaxation time can be estimated by the formula
t ~ l2 / D ,
where l — is the characteristic size of a particle (macroion); D — the coefficient of its diffusion in the electrolyte solution.
• Increased conductivity of double layer counterions.
According to [7, p. 56], the dependence of the low-frequency dielectric constant on the pH value was also found.
Taking these conditions into account makes it possible to explain the high values of the dielectric constant on the basis of existing ideas about the structure of the diffuse part of the double layer.
When developing the theory under consideration, the system of equations (1)-(5) is also taken as a basis. The distribution of potential U and concentration c± in a dispersed medium is sought in the form [7, p. 63]
V = <Po + ° 0 + U1-.
c±= co"+ Cf .
(17)
1) One should remember that in the Frenkel - Trukhan method, the dispersion medium is assumed to be non-conductive.
Here 00 is the equilibrium distribution of the electric potential in the double layer of the particle in the absence of an external electric field; p0 — distribution of electric potential under the influence of an external field; U1 — a polarization additive that is linear in the external field (the assumption is made about a sufficiently small value of the external field for the validity of the sufficiency of only the additive U1);
c± — polarization concentration distribution; c^ — concentration of ions in free electrolyte; — deviation of concentrations c^ in the vicinity of the double layer; cf — polarization additions to concentrations, that are linear in the external field.
Stationary polarization
Further, the authors consider a double layer of a spherical particle in a constant electric field [7, Section 3.1] (in this case, the continuity equation (3) is
1
s
reduced to the stationary case Vj+ = 0). In this case, the assumption is made about the approximation of a thin double layer (16), which allows us to consider it as locally flat. The polarization potential and concentrations are defined as follows
U = p + Ф1:
Cf= cf
+ П
As a result of complicated calculations in [7], the following results were obtained:
• the Boltzmann distribution is preserved for the ion concentration and for a polarized thin double layer, with the difference that a pre-exponential factor depending on the polar angle 0 is added to the standard Boltzmann distribution of the ion concentration (0 = 0 corresponds to the direction of a constant electric field), since the concentration of ions changes along the outer boundary of the double layer [7, p. 70], and, therefore, the structure of the polarized double layer coincides with the structure of the equilibrium double layer (in the absence of an external electric field);
• the field strength arising inside the double layer during its polarization is Ka times higher than the field strength at its outer boundary; in such a situation, the fulfillment of the conditions for the continuity of ion fluxes inside the double layer is possible only due to the almost complete compensation of the normal components of electromigration and diffusion fluxes of ions, which leads to the necessity of fulfilling the Boltzmann distribution [7, p. 75];
• when a polarization of the double layer is weak (i.e., when the conditions c+ » y+ + C+ from (17) are satisfied), the smallness of the Peclet criterion
Va
Pe = — ^ 1 (see above) automatically follows,
which makes it possible to reduce system (1)-(5) to two uncoupled systems: hydrodynamic ( 1), (2) and electrochemical (3)-(5) [7, p. 76].
Polarization in an alternating electric field
Further, the authors generalize the theory obtained for the stationary regime for the case of alternating fields of not too high frequency, while retaining the simplifying factors taken into account in the theory of stationary polarization. The system of equations (1)-(5) is again taken as a basis, but with the continuity equation in non-truncated form (3). In this case, the assumption is again made of a relatively small concentration difference, leading to a small Peclet number and, consequently, to the independence of equations (1), (2) and (3), (4).
Further, the authors find a limitation on the frequency of the external electric field (expression (III.88) [7, p. 77]), under which local equilibrium is
maintained in the double layer and, when polarized in an alternating electric field, for the ion concentration in the double layer can be used the Boltzmann distribution [7, p. 78] (quasi-stationary solution). Due to the small Peclet numbers of the problem under consideration, the authors replace the condition of approximate local electroneutrality outside the double layer by the condition of complete electroneutrality.
c+
On this basis, a function c = — = — is introduced for
z- z+
which the ansatz is valid [7, p. 77]
c = c0 + c1 ( r,O) e
(18)
and for the potential outside the particles
p = p(r,e)e~ia>. (19)
The quasi-stationarity condition is satisfied if the condition [7, p. 78]
2n 1 t
T = — >—-, or rn^lnvK , (20) a vk
where v is the kinematic viscosity of the dispersion medium.
The solution for p and c is obtained as follows
■ = c0 + dccr (r)cosOe
(21)
P = 1P0 + dc
D~ - D+
z+D + + z " D~
. (r)cosO +
d
+—prCOsO fe r I
Po =-E
( 1 a ^
r +--7
2 r v ^ '
cosO.
(22) (23)
where E is the amplitude of the harmonic electric field E = ; dc, dv and cr are defined in [7, p. 80].
After substitution of dv into formula (II.39) as
noted in [7, p.52], the following expression for the discrepancy between the dielectric permittivities of the polarized dispersed mixture and the dispersion liquid (ie, the dielectric increment As = s* - £1) is
obtained with an accuracy to (Ka) 2 [7, p. 81]:
As =
9 (Aq - M )(1 + W + W2) (24)
— CCS1 (Ka) 2.
4 (A + AW )2 + (AW + A2W2)
The value of the static dielectric increment
2)
Aest = lim As is equal to ) [7, p. 81]
a^-0
9 , s 2 Aa2 - A2a, Asst = — as1 (Ka ) 21
st 4 1V ' A,2
(25)
Here the values of A1, A2, a1 and a2 are defined in
[7, p. 80], W = a where r0 = J2De / a characterizes r0
the thickness of the diffuse layer at frequency
a;De = D+D"(z+ + z")/(z+D+ + z-d) [7, p. 77].
Next, let's turn to the dependency
w=a=
Vi
a a
V2De/ a j2De'
(26)
At frequencies when W ~ 1 the polarization of the particle lags significantly behind the applied field in phase and the real part of the dipole (polarization) moment changes most rapidly. Thus, the slow process of establishing the concentration distribution outside the double layer leads to a significant lag of the particle polarization from the external field at frequencies
a ~ D. At these frequencies, high values appear at a
a
small thickness of the double layer, when Ka » 1 [7, p. 81].
For a symmetric electrolyte with an accuracy to (Ka) 2, the value of the imaginary part s" of the complex dielectric constant s*of the considered dispersed medium (s" = Im (s *)) is determined by the expression [7, p. 83]
" 9 f S2 s = — as1 (Ka) x
(A1a2 - A2a1 )(AW + A2W2)2 A1 ( A1 + AW )2 + ( A1W + A2 W2 )
(27)
2) The value Asst follows from (24) trivially due to the fact that W (a = 0) = 0. See dependence (26) for W (a).
It can be seen from (27) that the value s" is small at very low frequencies if W ^ 1 and at high frequencies a if W ^ 1, reaching a maximum value if W ~ 1, while, according to dependence (27), the frequency at which the maximum of dielectric losses is reached, i.e. s" — max will be inversely proportional to the square of the particle radius.
[7, Fig. 11] Dependence of the increment of the real part of the dielectric constant As / a and the specific value of the imaginary part of the dielectric constant s"/ a on the logarithm of the value W for aqueous suspensions of spheres for Ka = 50 . Curves 1 and 1' are obtained at a value t of the potential that is one and a half times more than in the case of curves 2 and 2'; As / a — not dashed curves, s"/ a — dashed curves
Fig. shows the dependence of the value of the specific increment of the real part of the dielectric constant As / a on the value W , and actually on -n— , on a logarithmic scale along the abscissa (the value a is still the value of the specific volume of spherical inclusions in a dispersion liquid). Fig. is taken from work [7, p. 81]. At the same time, the figure shows the dependences of the value s"/ a . The figure shows the pattern of the values As / a and s"/ a , as well as the direct dependence of these values on the value of t -potential.
CONCLUSIONS
As a result of this review, it follows that when the condition for the presence of a thin double layer around nonconducting dispersed particles in a conducting dispersion medium is satisfied, a possibility of giant dispersion of the real component of the dielectric constant appears, which can significantly change the electrophoresis rate of such particles.
The work was carried out at the Institute for Analytical Instrumentation of the Russian Academy of Sciences within the framework of State Assignment 075-00780-20-00 on topic No. 0074-2019-0013 of the Ministry of Science and Higher Education.
a
r
0
X
2
REFERENCES
1. Sharfarets B.P., Kurochkin V.E. [On the attempt to use the giant dispersion of the dielectric constant of a disperse system to increase the speed of electrophoresis of dispersed particles]. Nauchnoe Priborostroenie [Scientific Instrumentation], 2020, vol. 30, no. 3, pp. 38-44. URL: DOI: 10.18358/np-30-3-i3844 (In Russ.).
2. Maxwell J.C. A treatise on electricity and magnetism. Vol. 1. Dover Publications, 1954. (Russ. Ed.: Maks-vell Dzh.K. Traktat ob elektrichestve i magnetizme. T. 1. Moscow, Nauka, 1989. 416 p.).
3. Wagner K.W. Electricity of the dielectric behavior on the basis on the Maxwell theory. Arch. Elektrotech., 1914, no. 9, pp. 371-392.
4. Truhan E.M. [Variance of dielectric constant of heterogeneous systems]. Fizika tverdogo tela [Physics of the Solid State], 1962, vol. 4, no. 12, pp. 3496-3511. (In Russ.).
5. Truhan E.M. Elektronnyj transport v subkletochnyh struktu-rah i ih modelyah. Razrabotka radiofizicheskih metodov is-
Contacts: Sharfarets Boris Pinkusovich, sharb@mail.ru
sledovaniya i kvantovo-elektronnyj analiz. Dis. ... dokt. fiz-mat. nauk [Electronic transport in subcellular structures and their models. Development of radiophysical research methods and quantum-electron analysis. Dock. phys. math. sci. diss.]. Moscow, Moscow State University, 1980. 44 p. (In Russ.).
6. Truhan E.M. Elektrofizicheskie metody issledovaniya. Konduktometriya neodnorodnyh materialov. Uchebnome-todicheskoe posobie [Electrophysical research methods. Conductometry of heterogeneous materials. Teaching aid]. Moscow, MFTI. 26 p. (In Russ.).
7. Duhin S.S., Shilov V.N. Dielektricheskie yavleniya i dvojnoj sloj v dispersnyh sistemah i polielektrolitah [Dielectric phenomena and a double layer in disperse systems and polyelectrolytes]. Kiev, Naukova dumka, 1972. 206 p. (In Russ.).
8. Frenkel Ya.I. [Electrical oscillations in disperse systems]. Kolloidnyj zhurnal [Colloid Journal], 1948, vol. 10, no. 2, pp. 148. (In Russ.).
Article received by the editorial office on 28.08.2020