УДК 537.32
Влияние пироэлектрического эффекта на фазовый переход в
сегнетоэлектрике
Старков А.С., Пахомов О.В., Старков И.А.*, Иванов В.Л.
ifmoenergo@gmail .com Санкт-Петербургский национальный исследовательский университет информационных
технологий, механики и оптики Институт холода и биотехнологий
*Institute for Microelectronics, Technische Universität Wien, A-1040 Wien, Austria
В работе представлена математическая модель для описания температурного гистерезиса в сегнетоэлектрических материалах. Модель основана на учете пироэлектрического эффекта в модели Ландау. Получены уравнения, описывающие температурный гистерезис сегнетоэлектрика, зависимость поляризации от скорости изменения температуры, а также характер температурной зависимости для электрокалорического эффекта в сегнетоэлектрике. Учет пироэлектрического тока в феноменологической модели Ландау приводит не только к сглаживанию температурной зависимости спонтанной поляризации, но и к исчезновению фазового перехода второго рода в рассматриваемом материале.
Ключевые слова: сегнетоэлектрики, электрокалорический эффект, пироэлектрический эффект.
В работе рассматривается обобщенная модель описания температурного гистерезиса в сегнетоэлектрических материалах. Модель основана на учете пироэлектрического эффекта в модели Ландау. Получены уравнения, описывающие температурный гистерезис сегнетоэлектрика, зависимость поляризации от скорости изменения температуры, а также характер температурной зависимости электрокалорического эффекта.
Измерения, проведенные в последние годы, показывают наличие температурного гистерезиса для широкого класса сегнетоэлектриков. Поляризация, диэлектрическая проницаемость и теплоемкость [1-3] зависят не только от температуры, но и от того, увеличивается она или уменьшается. Более того, результаты зависят от скорости нагрева [4,5]. Подобное явление не может быть описано в рамках термодинамической модели Ландау. Цель нашего исследования - создание обобщенной модели, описывающей температурный гистерезис.
В случае собственного сегнетоэлектрического фазового перехода свободная энергия Гиббса имеет вид [6]
аР2 ЬР* сР6 С= — + — + — -ЕР. (1)
Здесь Р - поляризация, Е - электрическое поле, о. = а0(Т — Тс}5 а0, Ъ, с- не
зависящие от температуры Т постоянные, Тс - температура Кюри. Для фазового
перехода первого рода Ь<0, в трикритической точке Ь=0. Для фазового
перехода второго рода Ь>0 и слагаемым с Р6 пренебрегают. Уравнение
дО/ = о состояния ' д Р имеет вид
аР + ЬР2 +сР5 -Е = 0. (2)
При учете динамических процессов уравнение (2) заменяют либо на уравнение Ландау-Халатникова [7] дР дС
= — ^^ (3)
дь дР - 1
либо на уравнение демпфированного осциллятора [6]
д2Р дР дС
= (4>
С точки зрения термодинамики неравновесных процессов соотношение (3) является естественным обобщением (2) на нестационарный случай. Уравнение (4) используется для описания колебаний подрешетки сегнетоэлектрика типа смещения. Коэффициент демпфирования считается малой величиной, т.к. иначе колебания отсутствуют. В дальнейшем мы будем рассматривать только уравнение релаксации (3), поскольку для осцилляционного уравнения (4) формулы выводятся аналогичным образом.
Физическая интерпретация (3) такова. Изменение поляризации
/ = дР/
эквивалентно наличию тока ' '£?£. Этот ток приводит к возникновению
динамической напряженности поля /;'фг «/, где а - «внутреннее»
а /
сопротивление. Величина /\а\ определяет время релаксации. Наибольшей
величиной сс. обладают сегнетоэлектрики-релаксоры.
Будем рассматривать поляризацию Р как функцию двух независимых переменных Е и Т. По формуле полной производной дР дРдЕ дРдТ
— =--+--. (5)
дг дЕ дг дг дг у '
Первое слагаемое в правой части (5) есть обычный электрический ток, а второе - пироток, который ранее в модели Ландау не рассматривался. Влияние электрического тока в (5) на поляризацию изучено в [6], поэтому исследуем вклад пиротока. Будем считать, что напряженность электрического поля Е-постоянна, а изменяется только температура T(t) = Vt+To, где V - скорость изменения температуры, Т0 - начальная температура. В этом случае уравнение (3) принимает вид дР _ дС
^дТ ~ дР' ®
где в=aV. Имеющиеся экспериментальные данные [5] указывают на существенно более сложный характер зависимости коэффициента в от скорости
V. По-видимому, это связано с необходимостью учета механических напряжений, что приводит к появлению третьего слагаемого в правой части (5)
- пьезотока. В отличие от (2) уравнение (6) является дифференциальным и позволяет более точно описывать зависимости, полученные экспериментальным путем.
В большинстве случаев пироток является малым и можно выписать приближенное решение (6)
Р**Р0 + Рруг, (7)
Где Р0 - квазиизотермическая поляризация (решение уравнения (2)), а Рруг
- пироэлектрическая добавка к Р0.
? ЛТ (8)
Р -'руг
дТ
д2в дР2
Т.к. коэффициент /? для нагревания /?+ и охлаждения р имеет разный
Р+ и охлаждении Р~ имеет вид
знак, то разность поляризации при нагревании [8]
р+_р-= (Л+-Юа0Р0
(9)
{аШР2 + ЪсР*)2
Отметим, что в осцилляционной модели (4) величина температурного гистерезиса от коэффициента ^ не зависит.
Полученные формулы, описывающие зависимость поляризации от скорости нагрева (7), (8) и температурный гистерезис (9), становятся непригодными, если д2С
(10)
дР:
= 0.
Р=Рг,
Для фазового перехода второго рода равенство (10) выполняется при Т=Тс. Точное решение уравнения (6) для С=0, пригодное при всех значениях температуры, удается получить при Е=0. Спонтанная поляризация при в>0 может быть записана как [8]
■ю=
%
вхр
~ао(Т-ТС)2
т
(11)
где ег/^) - интеграл вероятности.
Для фазовых переходов первого рода условие (10) выполнено при температурах Т1 и Т2 , когда уравнение (2) имеет кратный корень (рис. 1). В случае нагрева, тогда для уравнения (2) поляризация при Т=Т2 меняется скачкообразно от значения Р2=Р(Т2) до 0. Решение (6) вблизи температуры Т2 будем искать в виде Р=Р2+дР. Тогда ¿Р должно удовлетворять уравнению
= K^spy + K2{T-r2)f
где Ki = &2G/dPз = + 20CP23 Д2 = '^дТ
где J4
dzG
Решение (12) имеет вид
дР =
(13)
где иф - функция Эйри.
Таким образом, для фазовых переходов первого рода ниже температуры Т2 поляризация отличается от квазиизотермической на величину Рруг (8); в окрестности температуры Т2 описывается через функции Эйри (13), и далее экспоненциально приближается к устойчивому значению.
В данной работе приводится обобщенная модель описания температурного гистерезиса в сегнетоэлектрических материалах. Учет пиротока в модели Ландау приводит к сглаживанию зависимости Р(Т), т.е. к исчезновению фазового перехода второго рода. Уравнения (7)-(9) описывают температурный гистерезис и зависимость поляризации от скорости изменения температуры. Разность значений физических величин при нагревании и охлаждении пропорциональна параметру в. В свою очередь, для фазовых переходов второго рода в точке Кюри поляризация не обращается в 0, а отлична от нее в малой
смещается от точки Кюри в сторону более высоких температур.
Данная работа выполнена при поддержке государственной аналитической программы «Развитие потенциала высшей школы 2009-2011». Регистрационный номер 2.1.2/5063.
Список литературы
1. T. M. Correia et al., Appl. Phys. Lett. , 95(18), 182904 (2009).
2. M.E. Guzhva, V.V. Lemanov and P.A. Markovin, Phys. Solid State, 43 (11), 2146 (2001)
3. E.D. Yakushkin, Phys. Solid State, 46 (2) , 335 2004.
4. F.Yan et al. Appl. Phys. Lett., 81 (24), 4580 (2002).
5. A.Dedyk et al., J. Vac. Sci. Technol. B, 29 (1), 01A501-1 (2010).
6. Гинзбург В. Л. Теория сегнетоэлектрических явлений // УФН. 1949. Т. 38. C. 490.
7. Ландау Л.Д., Халатников И.М. // ДАН СССР. 1954. Т. 96. С. 469
окрестности порядка
Максимум электрокалорического эффекта
8. А. С. Старков, О. В. Пахомов, И. А. Старков, Письма в ЖЭТФ, 91:10 (2010), 556-560
Influence of pyroelectric effect on phase transition in a ferroelectric
material
Starkov A.S., Pakhomov O.V., Starkov I.A.*, Ivanov V.L. ifmoenergo@gmail .com
Saint-Petersburg National Research University of Information Technologies, Mechanics and Optics
Institute of Refrigeration and Biotechnologies *Institute for Microelectronics, Technische Universität Wien, A-1040 Wien, Austria
In operation the mathematical model for the description of a temperature hysteresis in ferroelectric materials is provided. The model is based on the accounting of pyroelectric effect in Landau's model. The equations describing a temperature hysteresis of a ferroelectric material, dependence of polarization on the speed of change of temperature, and also nature of temperature dependence for electrocaloric effect in a ferroelectric material are received. The accounting of pyroelectric current in Landau's phenomenological model leads not only to smoothing of temperature dependence of spontaneous polarization, but also to second order phase change disappearance in a researched material. Keywords: ferroelectric material, electrocaloric effect, pyroelectric effect.