К СТАТИСТИЧЕСКОЙ ТЕОРИИ ДИЭЛЕКТРИЧЕСКИХ СВОЙСТВ

advertisement
ДОКЛАДЫ АКАДЕМИИ НАУК РЕСПУБЛИКИ ТАДЖИКИСТАН
2013, том 56, №5
ФИЗИКА
УДК 537.266
Академик АН Республики Таджикистан С.Одинаев, Р.С.Махмадбегов*
К СТАТИСТИЧЕСКОЙ ТЕОРИИ ДИЭЛЕКТРИЧЕСКИХ СВОЙСТВ
РАСТВОРОВ ЭЛЕКТРОЛИТОВ
Академия наук Республики Таджикистан,
*
Таджикский национальный университет
На основе соотношения между токами смещения и проводимости, а также микроскопического определения плотности тока проводимости, найдено выражение между комплексными коэффициентами диэлектрической проницаемости и удельной электропроводностью. Разделив действительную и мнимую части последних, получены аналитические выражения для коэффициентов диэлектрической проницаемости и диэлектрических потерь растворов электролитов.
Ключевые слова: коэффициент удельной электропроводности – коэффициент диэлектрической
проницаемости – диэлектрические потери – потенциал межмолекулярного взаимодействия – радиальная функция распределения.
Электрические свойства диэлектрических веществ описываются посредством коэффициента
диэлектрической проницаемости. При слабых электрических полях последний не зависит от напряжѐнности поля и со временем остается постоянным, а в случае переменных полей коэффициент диэлектрической проницаемости зависит от частоты. В этом случае в диэлектрике появляется связь
между коэффициентами электропроводности, диэлектрической проницаемости и диэлектрических
потерь, которые зависят от частоты электрического поля. Определение этой связи для жидких диэлектрических веществ, в частности для водных растворов электролитов, на основе микроскопической теории является сложной задачей.
В работах [1-3] дана теория электропроводности водных растворов электролитов в случае
слабых стационарных электрических полей и исследовано влияние высокой частоты на электропроводность и коэффициент диэлектрической проницаемости растворов электролитов при наличии переменного внешнего электрического поля во времени. Подробно изложена теория диэлектрической
проницаемости 1 ( ) и диэлектрической потери  2 ( ) для однородных диэлектрических веществ,
находящихся под действием как статических, так и зависящих от времени внешних электрических
полей, не зависящих от пространственных координат. Выведены формулы Дебая для частотной зависимости комплексной диэлектрической проницаемости  () в различных приближениях и моделях.
В основу исследования теории Дебая положена гипотеза о том, что установление равновесия во времени происходит по экспоненциальному закону.
В работах [4-8] определены как статические  S , так и динамические  () коэффициенты диэлектрической проницаемости диэлектрических веществ, на основе молекулярно-кинетической тео
Адрес для корреспонденции: Одинаев Саидмухамад. 734025, Республика Таджикистан, г. Душанбе, пр. Рудаки, 33, Президиум АН РТ. E-mail: ods42@mail.ru.
381
Доклады Академии наук Республики Таджикистан
2013, том 56, №5
рии, методами временных корреляционных функций и коллективных переменных. В различных приближениях для 1 ( ) и  2 ( ) получены уравнения Коула-Коула, Девидсона-Коула, Гаврильяк и
Негами, Вильямса-Ваттса и Фуосса-Кирквуда. Уравнение Фуосса-Кирквуда является достаточно простым для объяснения экспериментальных результатов, а его параметры связаны с параметрами уравнений Коула-Коула и Девидсона-Коула. В определѐнных приближениях эти уравнения переходят в
уравнение Дебая.
При детальном анализе экспериментальных данных по диэлектрическим параметрам концентрированного раствора NaCl в [9] показано, что диэлектрические параметры описываются релаксационной моделью Коула-Коула. В [10] изучены статическая диэлектрическая проницаемость и диэлектрические спектры 1,3-пропандиола и его бинарных растворов с н-пропанолом в диапазоне частот от
9.4 до 72.4 ГГц, а также плотность в широком интервале температур. Обнаружено существование
дисперсии диэлектрической проницаемости 1 ( ) в диоле и в растворах в изученном интервале частот и температур, а также установлено существование спектра времѐн дипольной релаксации. Измерение диэлектрической проницаемости водных растворов акриламида в интервале концентраций 0-10
мол.% амида на трѐх частотах 9.4; 14 и 25 ГГц при температурах 288, 298 и 308 К показано в [11].
Построены диаграммы Коула-Коула и определены времена диэлектрической релаксации. Расчѐт эквивалентной электропроводности и исследование кинетики диффузии ионов, согласно данным измерений электропроводности водных растворов селенита натрия различной концентрации и температуры, приведены в [12, 13]. На основе закона действия масс в среднесферическом приближении в [14]
разработан подход для описания влияния ионной ассоциации на растворитель и диэлектрическую
проницаемость в растворах электролитов.
Следует отметить, что в основе всех экспериментальных и теоретических работ физики диэлектриков лежат формулы Дебая для исследования частотной зависимости динамических коэффициентов диэлектрической проницаемости 1 ( ) и диэлектрических потерь  2 ( ) . Эти коэффициенты в подынтегральных выражениях содержат функцию y( ) , которая описывает распределение времѐн релаксации и является неизвестной. Для подробного определения частотной дисперсии этих коэффициентов необходимо знание явного вида y( ) и еѐ связи со структурой диэлектрических веществ, то есть с радиальной функцией распределения [15].
Целью настоящей работы является исследование частотной дисперсии динамических коэффициентов диэлектрической проницаемости 1 ( ) и диэлектрических потерь  2 ( ) растворов электролитов и их связи со структурой раствора.
В качестве исходного воспользуемся системой уравнений Максвелла, которые связывают
средние плотности зарядов и токов (всех, как внешних, так и индуцированных) со средними значениями напряжѐнностей электрического и магнитного полей. Внешние заряды и токи можно задать,
но индуцированные или наведѐнные заряды и токи сами зависят от значений напряжѐнности полей.
Следует отметить что уравнения Максвелла выражают систему основных законов электромагнитных
процессов в произвольной среде и форма их записи является произвольной. При описании медленных
382
Физика
С.Одинаев, Р.С.Махмадбегов
процессов свойства среды характеризуются заданием коэффициентов удельной электропроводности,
диэлектрической и магнитной проницаемости, значения которых берутся из опыта.
В случае переменных полей, как и в статическом случае, для описания электрического поля
внутри диэлектрика требуются две полевые характеристики, такие как напряжѐнность электрическо-
 
 
го поля E(q, t ) и электрической индукции D(q, t ) :
 
 
D(q, t )   0 E (q, t ) ,
(1)
где  0 – электрическая постоянная,  – приведѐнный коэффициент диэлектрической проницаемости.

Согласно [3], если электрическое поле E (t ) меняется периодически со временем, то электри-

ческая индукция D(t ) также меняется периодически со временем, но будет существовать некоторый


сдвиг фаз между E и D . В этом случае диэлектрическая проницаемость  () в (1) становится комплексной
 ( )  1 ( )  i 2 ( ) ,
(2)
где 1 ( ) – динамический коэффициент диэлектрической проницаемости,  2 ( ) – коэффициент
диэлектрической потери, пропорционально потерям энергий в диэлектрике: при   0 ,  2 ( )  0 ,
1 ( )   S , а при    ,  2 ( )  0 , 1 ( )    , то есть насыщенное значение коэффициента
диэлектрической проницаемости при наивысших частотах.
Для получения комплексной диэлектрической проницаемости в виде (2) в качестве исходного
воспользуемся дифференциальным уравнением из системы уравнений Максвелла (формула (1.17)
работы [16]) в виде


D
E 
  0 
 jпр ,
t
t

(3)

где jпр   Е – плотность тока проводимости, созданная свободными индуцированными зарядами.


D
Здесь
– сумма тока смещения, вызванного изменением полного электрического поля E и индуt


цированного тока jпр , то есть весь ток, за исключением заданного внешнего тока jв . По аналогии с
формулой (2.124 а,в,с) работы [8], в первый член правой части (3) внесѐн   , учитывающий диэлектрическую проницаемость растворителя.
Совершая в (3) Фурье-преобразование по времени, с учѐтом формулы (1) имеем
 ( )    
i
 0
383
 ( ) .
(4)
Доклады Академии наук Республики Таджикистан
2013, том 56, №5
Формула (4) определяет связи между комплексными коэффициентами диэлектрической проницаемости () и удельной электропроводностью  ( ) . Они являются аналогами формул (2.13) и
(А1.10)-(А1.12) в работе [3], выражения (2.124 а,в, с ) в работе [8], формул (2.42) работы [16] и
(1.16;1.17) работы [17]. Согласно [16], из соотношения (4) следует, что в случае быстрых процессов
нет разницы в описании свойств диэлектрика и проводника. И в том и в другом случаях свойства
среды можно характеризовать комплексной диэлектрической проницаемостью () . Также вместо
диэлектрической проницаемости можно для описания свойства среды использовать комплексную
электропроводность  ( ) . Поэтому попытаемся методом молекулярно-кинетической теории определить аналитическое выражение комплексного коэффициента удельной электропроводности  ( )
в зависимости от состава, термодинамических параметров и структуры раствора, то есть радиальной
функции распределения g (r ) , а на основе (4) комплексного коэффициента диэлектрической проницаемости () . Ранее в [18,19] частично исследована частотная дисперсия коэффициента удельной
электропроводности для ионных жидкостей и растворов электролитов.
В качестве исходного принимаем кинетическое уравнение для одночастичной функции рас-

пределения f a ( xa , t ) , структурных единиц сорта a растворов электролитов с обобщѐнным потенциалом Власова [20]:
f a pa f a
f a
f a0
 
 

 ea E (q1 , t )   Fa (q1 , t )   I c ( f a ) ,
t ma q1
pa
pa
(5)
где



Fa (q1 , t )    Fab (l ) (t   )  e (t  )  Fab (l )  Fab (c)   e (t  ) Fab (el ) d ,
(6)
b

Ic ( fa )  a 
pa
 p a
f a 
 2Фab r  r 





f a  kT (q1 , t )   , Fab (l )  nb 
g ab (r )  u  (q2 , t )  u  (q1 , t )  dr ,

2
2
pa 
r
r
 ma
 2 Kab r  r 



Fab (c)  nb 
g ab (r )  u  (q2 , t )  u  (q1 , t )  dr ,
2
2
r
r

e
e   2Фab r  r 


Fab (el )  nb  b  a  
g ab (r )dr E  (q1 , t ) , K ab (r ) 
2
2
r
 b  a  r
 
q2  q1

Φab ( y)
Z ab ( y )dy ;
y
n ,

 
xa  (q a , pa ) ,

  
r  q2  q1 / d ab ;
  g (r ) 
 g (r ) 
Z ab ( y)  n  ab    T  ab 
  n T
 T n

1
,  
 P / T
n CV

fa0 ( pa )  na (2 mkT )3/2 exp( pa2 / 2ma kT ) – равновесная максвелловская функция распределения;
 

ma , qa , p a – масса, координата и импульс частиц сорта a , p a  pa  ma (q1 , t ) – относительный
384
Физика
С.Одинаев, Р.С.Махмадбегов
импульс, ea  za e , eb  zb e , e – элементарный заряд; z a , zb ,  a ,  b – валентность и коэффициенты

трения частиц сортов a и b , E (q1, t ) – компоненты напряжѐнности электрического поля, Фab ( r ) –
потенциальная энергия взаимодействия между структурными единицами раствора, gab (r ) – равновесная радиальная функция распределения, описывающая равновесную структуру раствора, которую,
согласно
[5,21],
считаем
 
ua (q1 , t )
известной,
–
вектор
смещения
частиц
сорта
a,
1
2
  ab   ab
 kT (a  b ) / dab
a b – феноменологическая частота структурной релаксации, k –
постоянная Больцмана, d ab  (d a  d b ) / 2 , d a и db – диаметры структурных единиц раствора сортов
a и b . Правая часть уравнения (5), то есть I c ( f a ) – столкновительный член Фоккера-Планка, получена в приближении парных взаимодействий. Он обеспечивает необратимость исходного уравнения
во времени, то есть возможность описания диссипативных процессов в растворах. Ради удобства по-






ложим, что координаты частиц qa  q1 , qb  q2 , qc  q3 и т.д.. Следует отметить, что частицы рас-

твора взаимодействуют посредством потенциала Фab ( r , ) , который состоит из суммы энергии
межионных (катионов и анионов) Фij , ионно-молекулярных Фis и Фjs , а также межмолекулярных
Фss взаимодействий [5]. Здесь  s  ( s ,  s ) -полярные углы, описывающие ориентацию диполя относительно оси, соединяющие центры масс взаимодействующих частиц.

Пользуясь методом импульсных моментов одночастичной функции распределения f a ( xa , t ) и

определением вектора плотности тока ja , на основе уравнения (5) для ja (q1 , t ) получено дифференциальное уравнение. Совершая в нѐм Фурье-преобразование по времени и сравнивая полученный
результат с Фурье-образом дифференциального закона Ома, для комплексного коэффициента удельной электропроводности (в изотропном случае), согласно  ( ) в [19], имеем:
 ( )  
a

na
ea2 
 ab
Gelab (0)  ,
1  
1  i a  a 
b ea (1  i ab )

(7)
a b dab2
где  a 
,   
, a  2 a , ab  2 ab ,
 a ab
(a  b )kT
ma
 a0 
na ea2
a

1
ab

e
na ea2 a
e   2Фab (r , )

g ab (r , )dr .
; Gelab (0)  d ab nb  b  a  
2
r
ma
 b  a  
Разделяя действительную и мнимую части  ( ) в выражении (7), имеем:


 ( )    a0 ( ) 1   1  ( a )( ab )  Gelab ( )  
a


b

 
i   a0 ( ) a 1   1  ab
a
a
b 

385

 ab
 Gel ( )  ,


(8)
Доклады Академии наук Республики Таджикистан
2013, том 56, №5
где
 a0
 ab
Gelab (0)
ab
 ( ) 
; Gel ( ) 
.
1  ( a )2
1  ( ab )2 ea
0
a
(9)
А теперь, учитывая (8) в (4), для комплексного коэффициента диэлектрической проницаемости () имеем:


 ( )      01   a0 ( ) a 1   (1   ab /  a )Gelab ( )  


.




i 01 1   a0 ( ) 1   (1  ( a )( ab ))Gelab ( )  
b


a
a
b
(10)
Сравнивая выражения (10) и (2), для динамических коэффициентов диэлектрической проницаемости 1 ( ) и диэлектрической потери  2 ( ) получим следующие аналитические выражения:


1 ( )      01  a0 ( ) a 1   (1   ab /  a )Gelab ( )  ,

a
b



 2 ( )   01 1  a0 ( ) 1   (1  ( a )( ab ))Gelab ( )  .
a

b

(11)
(12)
Формулы (11) и (12) содержат вклады как трансляционной  a0 () , так и дискретных времѐн
структурной релаксации Gelab () . Выражения Gelab () , наряду с явной зависимостью с частотой процесса, состава, термодинамических параметров, ещѐ в подынтегральных выражениях, согласно (9),
содержат потенциал взаимодействия между структурными единицами раствора Фab (r , ) и равновесной радиальной функцей распределения g (r, ) . Правые части этих выражений являются аналогами формул (10)-(12) работы [3], с той лишь разницей, что здесь подынтегральные выражения являются известными функциями, когда в [3] функция y( ) является неизвестной. При определѐнном
выборе модели раствора, на основе [5,21], можно выбрать явный вид Фab (r , ) и g (r, ) , на основе
которого можно провести численный расчѐт 1 ( ) и  2 ( ) в широком интервале частот и термодинамических параметров состояния.
Поступило 05.04.2013г.
Л И Т Е РАТ У РА
1. Семенченко В.К. Физическая теория растворов. – М.-Л.: Госиздат технико-теоретической литературы, 1941, 382 с.
2. Харнед Г., Оуэн Б. Физическая химия растворов электролитов. – М.: Изд. Иностранной литературы, 1952, 628 с.
3. Фрѐлих Г. Теория диэлектриков. – М.: Изд. Иностранной литературы, 1960, 251 с.
386
Физика
С.Одинаев, Р.С.Махмадбегов
4. Браун В. Диэлектрики. – М.: Изд. Иностранной литературы, 1961, 326 с.
5. Юхновский И.Р., Головко М.Ф. Статистическая теория классических равновесных систем. – Киев:
Наукова думка, 1980, 372 с.
6. Богородицкий Н.П., Пасынков В.В, Тареев Б.М. Электрические материалы. Л.: Энергия, 1985,
304 с.
7. Потапов А.А., Мецик М.С. Диэлектрическая поляризация. – Иркутск: Изд. Иркут. унив-та, 1986,
264 с.
8. Barther J.M.G., Krienke H., Kunz W. Physical chemistry of electrolyte solution. – Modern aspects.
Steinkopff: Darmstadt; New York: Springer, 1998, 401p.
9. Шарков Е.А. – Журн. физ. химии, 1984, т. 58, №7, с. 1705-1710
10. Журавлев В.И., Дуров В.А., Усачѐва Т.М., Шахпаронов М.И. – Журн. физ. химии, 1985, т. 59, №7,
с. 1677-1680
11. Ястремский П.С., Родникова М.Н., Маркова В.Г., Калаева М.И. – Журн. физ. химии, 1998, т. 62,
№3, с. 779-781
12. Шапошник В.А. – Вестник ВГУ, серия: Химия, биология, фармация, 2003, №2, с. 81-82.
13. Влаев Л.Т., Гениева С.Д. – Журн. структур. химии, 2004, т. 45, №5, с. 870-871
14. Holovko M.F., Kapko V. – Acta Chim. Seor., 2009, v.56, pp. 203-208
15. Одинаев С., Махмадбегов Р.С. – Изв. АН РТ. Отд. физ.-мат, хим., геол. и техн. н., 2013, №1 (150),
с. 47-59
16. Климонтович Ю.Л. Статистическая теория неравновесных процессов в плазме. – М.: МГУ, 1964,
281 с.
17. Агранович В.М., Гинзбург В.Л. Кристаллооптика с учѐтом пространственной дисперсии и теория
экситонов. – М.: Наука 1965, 374 с.
18. Odinaev S., Odjimamadov I. Modern phys. Letters B, v. 15, nos. gandro, 2001,
19. Odinaev S., Odjimamadov I. Condensed Matter Physics, 2004, v.7, №4(40), pp.735-740.
20. Odinaev S., Odjimamadov I. – About one kinetic equation with generalized Vlasov s potential. // Jnt.
conf. PLMMP, Kyiv, 2001, Abstracts September 14-19, p. 35.
21. Каплан И.Г. Межмолекулярная взаимодействия. – М.: БИНОМ. Лаборатория знаний, 2012, 394 с.
С.Одинаев, Р.С.Мањмадбегов*
ОИДИ НАЗАРИЯИ СТАТИСТИКИИ ХОСИЯТЊОИ ДИЭЛЕКТРИКИИ
МАЊЛУЛЊОИ ЭЛЕКТРОЛИТЇ
Академияи илмњои Љумњурии Тољикистон,
*Донишгоњи
миллии Тољикистон
Дар асоси ифодаи байни љараёнњои кўчиш ва индуксионї, инчунин ифодаи муайян намудани зичии љараёни гузаронандагї, робитаи байни коэффисиентњои (зарибњои) мавњуми нуфузпазирии диэлектрикї ва љараёни гузаронандагии хос ёфта шудаанд. Бо ёрии људо намудани
ќисми њаќиќї ва мавњуми коэффисиентњои мазкур, ифодањои аналитикї барои зарибњои нуфузпазирии диэлектрикї ва талафшавии диэлектрикии мањлулњои электролитїњосил карда шудаанд.
387
Доклады Академии наук Республики Таджикистан
2013, том 56, №5
Калимањои калидї: зариби электрнаќлкунии хос – зариби нуфузпазирии диэлектрикї – талафшавии
диэлектрикї – потенсиали таъсири мутаќобила – функсияи таќсимоти радиалї.
S.Odinaev, R.S.Mahmadbegov*
TO STATISTICAL THEORIES OF THE DIELECTRIC PROPERTIES OF
SOLUTIONS OF ELECTROLYTES
Academy of Sciences of the Republic of Tajikistan,
*
Tajik National University
On the basis of a relation between dielectric displacement and conductivity currents, also microscopic definition of density of a current of conductivity, is obtained expression between complex coefficients the
dielectric of permeability and a direct-current conductivity. Parting the real and imaginary parts of the last
analytical expressions for coefficients the dielectric of permeability and the dielectric losses of solutions of
electrolytes are gained.
Key words: direct-current conductivity coefficient – the dielectric of permeability coefficient – the dielectric
of losses – intermolecular interaction potential – an allocation radial distribution function.
388
Download