5. Грузинов, Б.Ф. Физика и техника полупроводников / [Б.Ф. Грузинов и др.]. - 1976. - № 10. - Вып. 3. - С. 497 -503.
6. Коржуев, М.А. Физика и техника полупроводников / М.А. Коржуев, Л.В. Аракчеева. - 1978. - № 12. - Вып. 11. - С. 2192 - 2196.
7. Коржуев, М.А. Физика и техника полупроводников / М.А. Коржуев, Л.Е. Шелимова, Н.Х. Абрикосов. - 1977. -№ 11. - Вып. 2. - С. 296 - 302.
8. Новикова, С.И. Физика твердого тела / [С.И. Новикова и др.]. - 1973. - № 15. - Вып. 11. - С. 3407 - 3411.
9. Katayama, S. Solid State Commun. / S. Katayama. -1976. - № 4; 19. - Р. 381 - 383.
10. Katayama, S. Phys. Rev. / S. Katayama, D.L. Mills. -1980. - B22. - № 1. - Р. 336 - 352.
11. Kristoffeland, N. Phys.stat. sol. / N. Kristoffeland, P. Konsin. - 1967. - 21. - K39; 1968. - 28, 731.
12. Kristoffel, N. Ferroelectrics / N. Kristoffeland, P. Konsin. - 1973. - 6, 3.
13. Konsin, P. Phys.stat.sol. (b) / P. Konsin. - 1976. - 76,
487.
14. Kobayashi, K.L.I. Solid State Commun. / [K.L.I. Ko-bayashi and etc.]. - 1975. - 17. - Р. 875 - 878.
15. Munemura, G. Solid State Commun. / G. Munemura, A. Morita. - 1978. - 28. - № 3. - Р. 273 - 275.
16. Takaoka, S. Phys. Rev. / S. Takaoka, K. Murase. -1979. - B20. - № 7. - Р. 2823 - 2833.
УДК 538.911
И.С. Касьяненко
Научный руководитель: доктор физико-математических наук, профессор В.И. Алтухов
МОДЕЛИРОВАНИЕ И РАСЧЕТ КОЭФФИЦИЕНТОВ ТЕПЛОПРОВОДНОСТИ И ТЕПЛОВОГО РАСШИРЕНИЯ СЕГНЕТОКЕРАМИКИ НА ОСНОВЕ ТВЕРДЫХ РАСТВОРОВ ЦИРКОНАТА-ТИТАНАТА СВИНЦА
На основе представления о мягкой моде предложена микроскопическая модель аномального температурного поведения теплофизических свойств сегнетокерамики на основе твердых растворов цирконата-титаната свинца (ЦТС). Результаты расчетов температурной зависимости коэффициентов теплопроводности и теплового расширения для сегнетокерамики ПКР-8 согласуются с данными опытов.
Коэффициент теплопроводности, коэффициент теплового расширения, аномальное поведение, мягкая мода, фазовый переход, сегнетокерамика.
Based on the concept of soft mode the paper presents а microscopic model of the anomalous temperature behavior of the thermal properties of ferroelectric on the basis of solid solutions of lead zirconate titanate (PZT). The results of calculations of the temperature dependence of thermal conductivity and thermal expansion for ferroelectric ceramics RCC-8 are consistent with the experimental data.
Thermal conductivity, thermal expansion coefficient, abnormal behavior, soft mode, phase transition, ferroelectric ceramics.
1. Введение. В работе, исходя из представления о роли мягкой моды юМ (Т) в сегнетоэлектрическом
фазовом переходе, проведены расчеты и дана интерпретация аномального поведения коэффициентов теплопроводности и теплового расширения около температуры фазового перехода ТС сегнетокерамики на основе ЦТС. Сегнетокерамики на основе твердых растворов РЪ(Т1,2г)03 (ЦТС) со структурой типа пе-ровскита обладают уникальными физическими свойствами и благодаря возможности их варьирования за счет изменения химического состава находят широкое применение в технике. К тому же рассматриваемая пьезокерамика относится к важному классу сег-неэлекстрических систем с разупорядоченными структурами, в которых могут реализоваться размытые фазовые переходы. Механизмы таких переходов в неоднородных и многокомпонентных системах являются сложными и недостаточно выясненными, что приводит к значительным трудностям при их описании. Недавно в [5], [6] были получены новые экспериментальные данные по теплофизическим
свойствам пьезокерамики на основе твердых растворов ЦТС (ПКР-8 и ПРК-7М) в области температур 300 - 800 °К и было обнаружено аномальное поведение коэффициентов теплопроводности X и теплового расширения а около сегнетоэлектрического фазового перехода при температуре ТС равной 600 и 460 °К для систем ПКР-8 и ПРК-7М соответственно.
В данном сообщении с учетом представления о мягкой моде юМ (Т) [1], [2], [13], предложена микроскопическая модель температурного поведения коэффициентов теплопроводности и теплового расширения в широкой области температур от 0 до 900 °К. Для сегнетокерамики ПКР-8 расчеты согласуются с данными опытов. Несмотря на ряд использованных приближений, данные расчетов качественно и количественно согласуются с экспериментом и с единых позиций мягкой моды объясняют аномальное поведение Х(Т) и а(Т) около ТС = 600 К для ПКР-8.
2. Коэффициент теплопроводности к(Т) ЦТС. Как видно из рис. 1 и 2, данные опытов показывают, что изменение теплопроводности с температурой имеет характер, свойственный для неупорядоченных и стеклообразных веществ, т.е. Х(Т) растет с увеличением температуры. При этом в точке Т = ТС на кривой Х(Т) наблюдается резкий излом, хотя в этой точке часто наблюдается особенность типа капса. Согласно [12], [14], [15] аномальное поведение Х(Т) около ТС обусловлено взаимодействием продольных акустических фононов с мягкой модой < (Т) = ®2 (Т)-82 ответственной за фазовый переход при Т = ТС и способной привести к заметному уменьшению теплопроводности около ТС, как при температурах меньших, так и при Т > ТС [3], [4], [16], [17]. На поведение Х(Т) в ЦТС влияет и скачок удельной теплоемкости АС (Т), связанный с тепловыми флуктуациями около ТС . Тогда коэффициент теплопроводности Х(Т) для сегнетокерамики на основе ЦТС в модели Дебая можно представить в виде [2], [9]:
Х(Т)= 3Т I Су (ю,Т)х(ю)-р(ю2)^®2, (1)
Су (ю,Т) = С0 (ю,Т) + АС (Т),
(2)
Со (<ю Т )= к\т^
± \2 Ню ( Ню Л 2 Йю ] тг
кТ
екТ -1 V У
АС (Т ) =
2т
а Т
2^/р2 + 4а-у- |Т - ТС
(3)
ю,
(Т) = ю2 (Т)-82, ю2 (Т) = гао \Тс - Т\,
1, Т > Тс,
2, Т < Тс.
(4)
Здесь С0 (ю,Т) - удельная теплоемкость для фононов с частотой ю ; АС (Т) - скачкообразное изменение теплоемкости за счет тепловых флуктуаций около температуры фазового перехода второго рода, близкого к трикритической точке ТС [9]; а , р , у -
коэффициенты соответствующего разложения термодинамического потенциала в теории Ландау; V -
средняя групповая скорость фононов; У0 - молярный объем; к - постоянная Больцмана. В модели
Дебая р(ю2 ) = 3ю/2ю3в (юБ - частота Дебая) и с переменной х = ю/юБ формулы (1) - (4) для тепло -проводности принимают вид:
ЧТ )=^Т
3Уо
2 1
9*[7] '1т(х)'
ех©/Т х4
(вх@1Т -1)
-йх +
4
Уо
1 + а0 г
Тс - Т
Б
(5)
Полное время релаксации фононов х-1 (х) = у( х) имеет вид [2]:
у( Х) = У0 (х) + Ут (х) ,
у 0 (х) = (А + В (х - у)4) -108 (сек-1), (у = юв -10-13),
(6)
где ут (х) - обратное время релаксации фононов за
счет рассеяния фононов на двухуровневых системах в неупорядоченных стеклообразных веществах Ут (х); А и В - коэффициенты в обратном времени
релаксации фононов у0 (х) в «затравочном - идеальном» поликристалите; ^х0 (х- усредненная по
частоте с удельной теплоемкостью время релаксации фононов в керамике; р(х) = 3х / 2 ; С1 - концентрация двухуровневых систем; х0 =ю0 / юБ - резонансная частота этих систем; © - температура Дебая. При учете скачкообразного изменения теплоемкости в (3), использовано усредненное время релаксации (х0 (х)) [7] с учетом монотонной зависимости
х0 (х) в неупорядоченных стеклообразных системах
с размытым фазовым переходом. Расчеты по формулам (5) - (6) дают хорошее качественное и количественное согласие с данными опытов (рис. 1). При этом были использованы следующие значение параметров теории: Х0 =(2,5/1,8)-1013 (моль/м - с2),
У0 = 1014 (с-1), * = 0,6-10-2 (Дж/с - м - К), й = 5 К, а0 = 1,32 (мэВ2/К), с, = 3, = 1,8-103, х0 = 0,42, А = 106, В = 0,2, © = 1000К (уменьшение или увеличение © на 20 - 30 % ведет лишь к незначительному изменению полученных результатов).
2
г =
3.5
3.5
Рис. 1. Теплопроводность Х(Т) сегнетокерамики на основе твердых растворов ЦТС (цирконата-титаната свинца РЬ(Т1, 7г) 03): 1 - расчет для кристаллов без вклада от процессов переброса (хи = 0 )
и аномальной теплоемкости при фазовом переходе и X0 = (2,5/1,8)-1013 (моль/м • с2); 2 - для поликристалла с вкладом в Х(Т) от аномальной теплоемкости АС(Т) около Тс = 600 °К при Х0 = (2,0/1,8) • 1013 (моль/м • с2) и ут (х) = 0 ; 3 - для сегнетокерамики ЦТС с резонансным рассеянием фононов на двухуровневых кластерах ут (х) (поликристаллах с локальной плотностью состояний Рт (х) = х21тg(х) для двухуровневых систем); 4 - точки - эксперимент [5], [6]
3. Модели для описания поведения коэффициента теплового расширения. В области размытого фазового перехода в ЦТС наблюдается резкое изменение коэффициента теплового расширения (КТР). В области около ТС КТР становится отрицательным с наименьшим значением при Т = ТС для ПКТ-8 порядка -1,1 -10 6 / К. С термодинамической точки зрения отрицательное значение КТР можно объяснить отрицательными значениями эффективного коэффициента Грюнайзена [8]. В целом поведение КТР а(Т) около ТС явно соответствует поведению мягкой моды юМ (Т) = ю02 (Т)-82[1], [2], [13] (82 - ангармоническая добавка для необращающейся при Т = ТС в нуль мягкой моды [2]). Отрицательные
значения КТР около ТС при этом можно связать с
неполным смягчением мягкой моды (4) за счет 82 и
мнимостью ее значений в узкой области температур
АТ = Т2 - Т1 = 17 К , где ю2 (Т )<8°. Для описания
поведения КТР твердых тел существует несколько разных теоретических подходов - моделей [8], [10], [11].
Поведение КТР сильно отличается для различных систем кристаллов, и для объяснения его температурного поведения требуются весьма различные
подходы и модели. Кроме того, имеется целая группа веществ, КТР которых в некоторой области температур имеет отрицательное значение. Такое поведение КТР кристаллов типа алмаза и сфалерита объясняется особенностями фононного спектра и знаком производных от частоты фононов ю по давлению
(dю / с1р ). В некоторой области низких частот возможны значения dю / dp < 0, и КТР будет иметь отрицательные значения. Необходимо также учитывать различие параметров Грюнайзена для поперечных уТ и продольных у1 акустических колебаний решетки. В [8] также отмечается, что знак у ; (I = Т, I) зависит от знака производной от упругих постоянных с у по половине периода решетки г (dc у / dr).
4. Расчет КТР в сегнетокерамике ЦТС. С учетом наличия в ЦТС сегнетоэлектрического фазового перехода с мягкой модой юМ (Т) для описания поведения КТР в ЦТС получаем:
а(Т) = аС (Т) + аг (Т), (Т) АЙюм (Т) АН 1-г—Т\ , 8 /п\
аС(Т) = —кт— = ~кт^г а° Т - Т +8o, (7)
11, T > Тс
[2, T < TC;
Т.. =
71, T > TC
T2.
T < Tc ;
av (Т ) =
Y r Xt
Vn
© f 1 x4ex©/Tdx
9R TU
i x©IT 1 \2
Г -
(8)
(9)
где * - постоянная уравнения состояния системы; 50 - отрицательное значение КТР в области мнимой частоты юМ (Т). Расчеты по формулам (7) - (9) дают хорошее качественное и количественное согласие с данными опытов (рис. 2). В расчете были использованы следующие значение параметров теории:
В> = У г Хг/У =(2,5/1,8 )-10-7 (1/Дж), v = 104 (м/с), У = 6 -10-6 (м3/моль), 80 = 3-10-6 (1/К), Т1 = 586 К , Т2 = 603 К.
На рис. 2 представлено сопоставление результатов вычислений Х(Т) и а(Т) с данными эксперимента [5], [6], при этом кривая 1 для Х(Т) получена в предположении асимметричного поведения скачка теплоемкости АС(Т) около Тс , т.е. при Т > Тс в соответствии с данными опытов были использованы значения АС (Т) = С0Т , где С0 = 5-10-3 Дж/с - м - К. Заметное (до 30 - 40 %) превышение при Т < 400 К рассчитанных значений Х(Т) над данными опытов, видимо, связано с необходимостью более точного определения плотности состояний рт (х) при больших значениях концентрации двухуровневых систем.
Как видно на рис. 2, результаты расчетов по предложенным выше моделям (5), (7) для Х(Т) и для а(Т) с учетом мягкой моды объясняют температурное поведение коэффициентов теплопроводности и
теплового расширения ЦТС и хорошо согласуются с данными опытов.
5. Обсуждение. Около температуры TC сегнето-электрического фазового перехода, как правило, на кривой теплопроводности Х(Т) наблюдаются типичные особенности различного вида [3], [4]:
(а) излом или особенность типа капса (скобки) направленные в «отрицательную» сторону, т.е. в сторону уменьшения значений Х(Т);
(б) явно выраженный в широкой области температур отрицательный минимум и относительно узкий (в виде провала), направленный вниз «пик»;
(в) аномалия типа «ступеньки» или шаг - точка разрыва первого рода;
(г) положительный, почти симметричный, мало выраженный пик на фоне широкой области уменьшения значений теплопроводности Х(Т);
(д) разрыв или точка разрыва второго рода при Т ^ Тс.
Такие особенности обусловлены критическим рассеянием акустических фононов на колебаниях мягкой моды и на «центральном пике». Это рассеяние и ведет к аномальному поведению Х(Т) при Т ^ Тс [3], [4], [12], [14] - [17]. На кривой ЦТ) в ЦТС в точке TC , однако, наблюдается нетипичный (положительный) излом, направленный в сторону увеличения значений Х(Т). По-видимому, это характерно для сильно неупорядоченных, стеклообразных и сегнетокерамических систем с размытыми фазовыми переходами. При таком поведении кривой Х(Т) естественно предположить, что в ЦТС важную роль играет скачок удельной теплоемкости AC (Т), обусловленный тепловыми флуктуациями около Тс [13]. В такой ситуации времена релаксации фо-
нонов
(т0 (ю), гс (ю)) осредняются [7],
r =
и
Рис. 2. Температурная зависимость коэффициентов теплопроводности и теплового расширения сегнетокерамики ЦТС (ПКР-8): кривая 1 - расчет Х(Т) с учетом скачка теплоемкости при Т = ТС (А = 1,22 -106, рт (х) = х 1т g (х) ); кривая 2 - расчет аномального поведения а(Т), обусловленного мягкой модой юМ (Т); точки и ромбики - данные экспериментов [5] и [6]
ведущим механизмом теплового сопротивления остается характерное для стекол резонансное рассеяние фононов на двухуровневых системах уm (х). При этом помимо аномалии на кривой ЦТ) около TC в расчете может появиться максимум (рис. 2) при температуре порядка 250 °К. Поэтому представляли бы значительный интерес экспериментальные данные по теплопроводности ЦТ), теплоемкости C (Т) и
КТР а(Т) ЦТС в более широкой области температур (50 - 1000 °К).
Согласно [8] отрицательная область значений коэффициента теплового расширения (КТР) в твердых телах может иметь весьма разную природу и разную форму, например, как в ПКР-7М [1], [2]. Это обстоятельство интерпретируется исходя из самых разных
подходов [8]. В нашей модели при ю0 (Т) < 82 имеем мнимые значения мягкой моды -8° < юМ (Т) < 0 (юМ (Т/8 ), что соответствует согласно [8] отрицательным значениям КТР с 80 =-3 -10-6 (1/К) в области Т < Т < Т2.
Остается неясным наблюдаемое слишком резкое возрастание а(Т) справа от Тс (Д = (1,602/1,38) • 10-7,
при Т <Тс и Д =(1,602/1,38)-10-4 (1/мэВ), при Т > Тс, рис. 2). Отметим, что еще более значительное возрастание а(Т) при Т > Тс наблюдается для ряда систем, в том числе в кристаллах триглицин-сульфата (ТГС) [9]. Несмотря на это, предложенные нами в работе модели для ЦТ) и для а(Т) с единых позиций мягкой моды хорошо качественно объясняют температурное поведение коэффициентов теплопроводности и теплового расширения ЦТС и результаты расчетов количественно согласуются с данными опытов.
Литература
1. Аксенов, В.Л. Рассеяние нейтронов сегнетоэлектри-ками / В.Л. Аксенов, Н.М. Плакида, С. Стаменкович. - М., 1984.
2. Алтухов, В.И. Основы теории кинетических свойств кристаллов с дефектами и фазовыми переходами: диэлектрики и сегнетоэлектрики / В.И. Алтухов. - Ставрополь, 2003.
3. Алтухов, В.И. Модели особенностей теплового сопротивления кристаллов с фазовыми переходами и дефектами / В.И. Алтухов, Б.А. Казаров, Н.В. Баландина // Изв. Самар. НЦ РАН. - 2007. - № 3. - С. 640 - 646.
4. Алтухов, В.И. Рассеяние фононов на точесных дефектах структуры, комплексах-наночастицах и типичные особенности теплового сопротивления реальных кристаллов и сегнетоэлектриков / В.И. Алтухов, А.Т. Ростова, Б.А. Казаров // Нано- и микросистемная техника. - 2006. - № 3. - С. 19 - 25 (I ч.); № 4. - С. 14 - 20 (II ч.).
5. Катаев, С.Н. Особенности теплопроводности и теплового расширения сегнетокерамики на основе ЦТС / [С.Н. Каллаев и др.] // Материалы международной конференции «Fizika-2005». - Баку, 2005. - С. 333 - 335.
6. Катаев, С.Н. Теплофизические свойства сегнетокерамики на основе ЦТС / [С.Н. Каллаев и др.] // Физика твердого тела. - 2006. - Т. 48. - Вып. 6. - С. 1099 - 1100.
7. Могилевский, Б.М. Теплопроводность полупроводников / Б.М. Могилевский, А.В. Чудновский. - М., 1972.
8. Новикова, С.И. Тепловое расширение твердых тел / С.И. Новикова. - М., 1974.
9. Струков, Б.А. Физические основы сегнетоэлектриче-ских явлений в кристаллах / Б.А. Струков, А.П. Леванюк. -М., 1995.
10. Ферми, Э. Молекулы и кристаллы / Э. Ферми. - М., 1947.
11. Френкель, Я.И. Введение в теорию металлов / Я.И. Френкель. - М. ; Л., 1950.
12. Jnoue M., J. Phys. Soc. / M., J. Jnoue. - 1968. - Jpn.
25.
13. Levanyuk, A.P. Defects and structural phase transitions / A.P. Levanyuk, A.S. Sigov // OPA. - Amsterdam, 1988.
14. Mante A.J., J. Volger Phys. Lett. / A.J., J. Mante. -
1967. - A24.
15. Steigmeier, E.F. Phys.Rev.168 / E.F. Steigmeier. -
1968.
16. Strukov, B.A. Heat transport properties of ferroelectrics and related materials / B.A. Strukov, A.A. Belov // Phase Transitions. - 1994. - V. 51. - P. 175 - 197.
17. Strukov, B.A. Study phonon scattering in displacive fer-roelectrics by means of heat conductivity measurements / B.A. Strukov, A.A. Belov, V.I. Altukhov // Ferroelectrics. - 1994. -V. 159. - P. 25 - 30.
УДК 662.641
П.В. Малыгин, В.К. Любое ИССЛЕДОВАНИЕ СТРУКТУРЫ, СОСТАВА И СВОЙСТВ ТОРФА
В статье рассмотрен вопрос использования торфа как перспективного вида топлива для выработки тепловой и электрической энергии на территории Северо-Западного региона. Приведены данные по запасам торфа в месторождениях Архангельской области, а также результаты исследования теплотехнических характеристик для некоторых из них. Проанализированы результаты элементного анализа торфа и торфяных брикетов, дифрактометрических исследований, а также данные микрофотографии, полученные на электронном микроскопе.
Торф, топливный брикет, элементный состав, теплотворная способность, кристалличность, электронная микроскопия.
The article considers the use of peat as a promising type of fuel for heat and power in the North-West region. The data on stocks of peat deposits of the Arkhangelsk region as well as the results of studying the thermal characteristics of some of them are presented. The results of elemental analysis of peat and peat briquettes, diffractometric research as well as microphotograph data obtained by an electron microscope have been analyzed.