Моделирование токов поляризации сегнетокерамики, возникающих при одновременном приложении электрического поля и механических
напряжений
1 2 Л.В. Жога , В.В. Коренева
1 Волгоградский государственный технический университет 2Волгоградская государственная академия физической культуры
Аннотация: Исходя из модели движения доменных границ и реологической модели, получены уравнения зависимости поляризации от приложенной механической, электрической или одновременной электромеханической нагрузки. Показано, что нелинейность тока в сегнетокерамике обусловлена зависимостью коэффициентов модели от тока сквозной проводимости и тока, связанного с эмиссией зарядов из ловушек. Полученные коэффициенты позволяют переходить от поляризации (макроскопического параметра) к доменным границам (мезоскопическому масштабу). Параметры модели зависят от изменения доменной структуры, ее взаимодействия с дефектами, поэтому время релаксации тока в сегнетокерамике описаны изменением соответствующих времен для 180°- доменов и не 180°- доменов. Выбранная реологическая модель и подход с учетом механизмов движения доменных границ даст возможность описывать скорость изменения поляризации, используя поведение дефектов (дислокаций, доменных границ). Модель позволит с единой точки зрения описать экспериментальные закономерности поведения токов, при электромеханическом воздействии.
Ключевые слова: сегнетокерамика, компьютерное моделирование, реологическая модель, ток поляризации, плотность тока поляризации, доменная граница, 180°- домены, не 180°- домены, механическая нагрузка, электрическая нагрузка, время релаксации.
Введение. Сегнетоэлектрические материалы находят широкое применение во многих областях современной техники: гидроакустике, оптике, радио- и измерительной технике [1, 2]. Условия эксплуатации поликристаллических диэлектриков во внешних механических и электрических полях делают актуальным исследование свойств таких материалов и их устойчивости к экстремальным нагрузкам, которая зависит от процессов переполяризации и динамики доменной структуры.
Для исследования кинетики переключения поляризации в сегнетоэлектриках применяется метод анализа переходных токов, возникающих при приложении импульсов электрического поля и содержащих информацию о закономерностях эволюции доменной структуры
[3]. Для поликристаллических сегнетоэлектриков этот метод, по-видимому, является единственно возможным, поскольку возможности прямого наблюдения движения доменов крайне ограничены, но в связи с этим, актуальна задача моделирования свойств сегнетоэлектриков, вклада движения различных доменных границ (ДГ) в переходный ток при воздействии не столько одного электрического поля (ЭП) или механической нагрузки (МН), сколько при приложении этих нагрузок одновременно.
Методика моделирования. Уравнения зависимости поляризации от приложенной нагрузки можно получить исходя из модели движения доменных границ под действием электрического поля, описанной для сегнетоэлектриков [4]:
тд . + к. = 2Р5Е (1)
¥ 8§д (х2) + к2х2 = 2рЕ, (2)
где 7} и Р - эффективные коэффициенты вязкого и сухого трения,
соответственно, кЛ и к2 - эффективные коэффициенты квазиупругой силы,
х1, х2 и , х2 ~ средние смещения и скорости доменной границы,
соответственно, Р3 - спонтанная поляризация образца, Е - напряженность
электрического поля.
Уравнения определяют вклад в эффективную диэлектрическую проницаемость и эффективные диэлектрические потери релаксационного (1), осуществляемого за счет непрерывного изменения порядка, и гистерезисного (2) осуществляемого за счет скачкообразного изменения направления поляризации в «заметаемых» доменной границей элементарных ячейках сегнетоэлектрика, механизмов движения доменных границ. Однако в представленной модели не описывается влияние механической нагрузки.
и
Обозначив Ь, d - размеры доменной области и учитывая связь Q = рз = рМ и уравнение (2) полученное в [4]:
Получаем уравнение изменения поляризации:
Fd kd _л _
-Р i +-f—:—f—T Pi = 2PsE
2Ps + f jCK dt + f jrdt 1 2PS + f jCK dt + f jrdt
4PS2 + 2ps (f jCKdt + \ jrdt)
FPi + kPi =-^---E (3)
d
Задавая время релаксации, как zi =F и
k
4Ps 2 + 2Ps (f jк dt + f jrdt)
коэффициенту =-----, уравнение (3) приобретает вид:
kd
Ti pi + Pi = aiE
При приложении механической нагрузки к образцу происходит деформация зерен сегнетокерамики, при этом возникают дополнительные локальные механические напряжения в местах соприкосновения зерен. При прохождении таких мест 90° доменные границы испытывают дополнительное сопротивление, что приводит к локальным изменениям потенциального профиля [5].
Такие механизмы рассмотрены в рамках уравнения движения доменных границ [6]:
F sgn (х)+kx = anusaiJkl +«2i2psE (4)
где uS и Gij - спонтанная деформация и механическое нагружение соответственно, коэффициент а11 определяется геометрией механического напряжения и спонтанной деформации, коэффициент а21 зависит от взаимной ориентации векторов напряженности электрического поля и вектора спонтанной поляризации [6].
Учитывая связь заряда и поляризации из (4), получаем уравнение изменения поляризации при приложении механической нагрузки:
и
Рё кё _л _
Р2 + ~-г—.-г~ТР2 =аии3ат +а212Р5Е;
2Р5 +1 ]СК Я +1 ]ГЖ 2 2Р5 +1 ]ск Л +1 ]ГЖ
хии, (2р5 +1 ]скЛ +1jrd) 2 + 2Р (Т ]ск Л + Г ))
РР 2 + кР2 =-—---°г]ы +-^---Е
Задавая время релаксации, как т2 = —, коэффициент
к
«2! (4Р2 + 2Р (Г ]СК Ж + Г ]гл} апи8 (2Р? + Г ]ск ё + Г ]гж)
а2 =----- и коэффициент р2 =-----,
кё кё
уравнение (4) приобретает вид:
т 2Р2 + Р2 =«2Е + ¡32а
Вид тока в сегнетокерамике обусловлен зависимостью коэффициентов модели а1 ,а2, /32 от плотности тока сквозной проводимости и плотности тока, связанного с эмиссией зарядов из ловушек. Полученные коэффициенты позволяют переходить от поляризации (макроскопического параметра) к ДГ (мезоскопическому масштабу).
Плотность тока при изменении поляризации сегнетокерамики (СК) можно описать системой уравнений:
р! + Р! =а!Е (5)
Р 2 + Р2 =а2 Е + 02° (6)
] = Р1 + Р 2 (7)
где параметры тх, т2 зависят от изменения доменной структуры и ее взаимодействия с дефектами [4], поэтому тх будем обуславливать изменением соответствующих времен для 180°- доменов и т2 для не 180°-доменов; р1 - поляризация, обусловленная наличием 180° - доменов и р2 -поляризация, обусловленная наличием не 180° доменов; коэффициенты ах, а2 и р2 определяются геометрией взаимного расположения векторов напряженностей электрического поля Е , компонентов тензора механических напряжений и направлений векторов спонтанной поляризации р в 180°
и
доменах и не 180° доменах; у - плотность тока при изменении поляризации. В модели используются усредненные по всем кристаллитам величины. Компоненты вектора Е и тензора напряжений определяются прикладываемым воздействием к образцу.
Механические напряжения растяжения при нагружении образцов рассчитывали по формуле для изгиба круглых пластин с недеформируемой нейтральной плоскостью [5].
3 1 + у а = — ■
Ъ 1 -V Ъ2 -а2Л 1п — +---
У а 1 + у 2с2 у
Р
2 ттк1
где а - радиус нагружающего пуансона; Ь - радиус нижней опоры; с -радиус образца; V - коэффициент Пуассона; И - толщина образца, Р -приложенная сила.
Моделирование нелинейных зависимостей между электрическим полем и механическими напряжениями с одной стороны, и поляризацией и деформацией, с другой, относится к разделу реологических моделей.
В основе электромеханической аналогии для описания поляризационных эффектов в сегнетоэлектрических телах лежат соответствия: обобщенной координаты - электрическому заряду; обобщенной скорости - току; коэффициента упругой податливости -емкости; обобщенной силы - электродвижущей силе. При переходе к сплошным средам силы заменяются механическими напряжениями, перемещения - деформациями и т.д. В итоге можно прийти к соответствиям: а ^ Е, £ ^ р, где Е - напряженность электрического поля; р -поляризация. Упругому элементу сопоставляется конденсатор, элементу вязкого трения - сопротивление, элементу сухого трения - двуполярный стабилитрон [7].
Реологические формулы для двуполярного стабилитрона и конденсатора можно представить следующим образом [7]:
и
Ее = 1 Ре; Е0 = Г р0 к
где индексами «е» и «0» отмечены индуцированные и неупругие составляющие соответственно, г, к - некоторые константы, физический смысл которых следует определить. Соединяя эти элементы последовательно, согласно уравнениям Е = Ее + Е0 и р = Ре = р0 , получаем:
гкР + Р = кЕ
заменяя гк = т, а константу к = а получим окончательно:
тР + Р = аЕ
Полученное уравнение должно учитывать поляризацию, обусловленную движением 180° и не 180° доменов. С другой стороны, зависимость поляризации в сегнетокерамике от внешнего электрического поля или механической нагрузки (в общем случае при одновременном их действии) влияет на вид плотности тока [8], при этом, наряду с током, обусловленным переполяризацией сегнетоэлектрических доменов, возникает «поток деформации», вызванный передеформацией сегнетоэластических доменнов (упругих двойников) [8]. В связи с этим введено уравнение, учитывающее зависимость поляризации от внешней механической нагрузки. Причем в уравнении так же должно учитываться влияние внешнего электрического поля.
Физический смысл коэффициентов ах, а2 и р2, (ст = ст1 -механические напряжения (здесь <?1, <?2 - решеточная и доменная составляющие механического напряжения соответственно)) позволяет использовать известные их зависимости от электрического поля и внешних механических напряжений как в интервале действия линейного пьезоэффекта, так в области внешних нагрузок, когда линейность перестает существовать.
Размерность коэффициентов, используемых в уравнениях (5-7) модели:
[к ] =
Г г- 1 тт 2
И= "
_p_
.Ее_
Кл2 Ф
9
м ■ Н м
Р
2
Н ■ м ■ с м ■ с
Кл2 Ф
или [г ] =
Е0 р
3
кг ■ м
——- = Ом ■ м ;
А2 ■ с3
[*] = [кг] = Ф ■ ^ = с .
м Ф
Если время релаксации процесса т = кг = ££0р, то р = г - удельное
сопротивление материала. Коэффициент /¡2 можно определить через связь
1 2
„1 м _
к/а = /2а, где = — = —. Тогда /2 = к/1 = ££0gj, где g1, - пьезоэлектрический
Е Кл ^
коэффициент напряжения.
Анализ изменения поляризации с помощью уравнений (5-7) не учитывает токи проводимости. Но, как показывает обзор литературы, реализовать такие условия в эксперименте практически невозможно. В модели (5-7) две составляющие поляризации: упругая и пластическая (необратимая), изменения поляризации после внешнего воздействия ЭП можно объяснить движением доменных и межфазных границ [9, 10]. После приложения механической нагрузки начинается процесс разрушения. Доменные границы являются препятствием разрушению (росту трещины) до тех пор, пока они закреплены на стопорах (заряженных дефектах). Если механизм движения гистерезисный, то при приложении критической нагрузки происходит срыв доменных границ с дефектов и трещина растет скачком. Релаксационное движение доменных границ с постоянным трением приводит к росту трещин за счет создания механических напряжений.
Вывод. Выбранная реологическая модель [7] и подход с учетом механизмов движения доменных границ [3, 4, 6] даст возможность с помощью уравнений (5-7) описывать скорость изменения поляризации, используя поведение дефектов (дислокаций, доменных границ). Модель позволит с единой точки зрения описать экспериментальные закономерности
поведения токов при электромеханическом воздействии и даст возможность провести анализ вкладов 180° и не 180° доменных стенок в токи в сегнетокерамике.
Литература
1. Фиговский О.Л. Нанотехнологии для новых материалов // Инженерный вестник Дона, 2012, №3. URL: ivdon.ru/ru/magazine/archive/n3y2012/1048.
2. Панич А.А., Мараховский М.А., Мотин Д.В. Кристаллические и керамические пьезоэлектрики // Инженерный вестник Дона, 2011, №1. URL: ivdon.ru/ru/magazine/archive/n1y2011/325.
3. Жога Л.В., Нестеров В.Н., Терех В.В., Габриэлян А.В. Исследование поведения токов поляризации методом компьютерного моделирования движения доменных стенок // Известия Волгоградского государственного технического университета. Серия «Электроника, измерительная техника, радиотехника и связь», вып. 7: межвуз. сб. науч. ст., 2013. № 3(106). Вып.7. С. 41-45.
4. Нестеров В.Н. Динамика доменных и межфазных границ в сегнетоэлектрических твердых растворах на основе цирконата - титаната свинца (компьютерный анализ): дис. канд. физ. - мат. наук. ВАСА. Волгоград, 1997. 263 с.
5. Бондаренко Е.И., Тополов В.Ю., Турик Л.В. Внутренние механические напряжения и электрический пробой поликристаллического титаната бария // ЖТФ. 1992. Т. 62. № 12. С. 155-158.
6. Нестеров В.Н., Кочергин И.В., Жога Л.В. Абсорбционный ток в нагруженной сегнетокерамике // Физика твердого тела. 2009. Т. 51, Вып. 7. С. 1439-1441.
7. Белоконь А.В., Скалиух А.С. Математическое моделирование необратимых процессов поляризации. М.: ФИЗМАТЛЛИТ, 2010. 328 с.
8. Шур В.Я., Кожевников В.Л., Пелегов Д.В., Иванов Р.К. Формирование и движение фрактальных доменных стенок в сегнетоэлектриках // Вестник ВГТУ. Сер. Материаловедение. 2000. В. 1.8. С. 36-40.
9. Boddu V., Endres F., Steinmann P. Molecular dynamics study of ferroelectric domain nucleation and domain switching dynamics // Scientific Reports, 2017, 7:806 URL: nature.com/articles/s41598-017-01002-0.pdf.
10. Lee Dong Hyun, Lee Younghwan, Yang Kun, Park Ju Yong, Kim Se Hyun, Reddy Pothala Reddi Sekhar, Materano M., Mulaosmanovic H., Mikolajick T., Jones J. L., Schroeder U., Park Min Hyuk. Domains and domain dynamics in fluorite-structured ferroelectrics // Applied Physics Reviews, 2021, 8:021312 URL: aip.scitation.org/doi/10.1063/5.0047977.
References
1. Figovskij O.L. Inzhenernyj vestnik Dona, 2012, №3. URL: ivdon.ru/ru/magazine/archive/n3y2012/1048.
2. Panich A.A., Marahovskij M.A., Motin D.V. Inzhenernyj vestnik Dona, 2011, №1. URL: ivdon.ru/ru/magazine/archive/n1y2011/325.
3. Zhoga L.V., Nesterov V.N., Tereh V.V., Gabrijeljan A.V. Izvestija Volgogradskogo gosudarstvennogo tehnicheskogo universiteta. Serija «Jelektronika, izmeritel'naja tehnika, radiotehnika i svjaz'», vyp. 7: mezhvuz. sb. nauch. st., 2013. № 3(106). Vyp.7. pp. 41-45.
4. Nesterov V.N. Dinamika domennyh i mezhfaznyh granic v segnetojelektricheskih tverdyh rastvorah na osnove cirkonata - titanata svinca (komp'juternyj analiz) [Dynamics of domain and interphase boundaries in ferroelectric solid solutions based on lead zirconate-titanate (computer analysis)]: dis. kand. fiz. - mat. nauk. VASA. Volgograd, 1997. 263 p.
5. Bondarenko E.I., Topolov V.Ju., Turik L.V. ZhTF. 1992. T. 62. № 12. pp. 155-158.
М Инженерный вестник Дона, №3 (2022) ivdon.ru/ru/magazine/arcliive/n3y2022/7523
6. Nesterov V.N., Kochergin I.V., Zhoga L.V. Fizika tverdogo tela. 2009. T. 51, Vyp. 7. pp. 1439-1441.
7. Belokon' A.V., Skaliuh A.S. Matematicheskoe modelirovanie neobratimyh processov poljarizacii [Mathematical modeling of irreversible polarization processes]. M.: FIZMATLLIT, 2010. 328 p.
8. Shur V.Ja, Kozhevnikov V.L., Pelegov D.V., Ivanov R.K. Vestnik VGTU. Ser. Materialovedenie. 2000. V. 1.8. pp. 36-40.
9. Boddu V., Endres F., Steinmann P. Scientific Reports, 2017, 7:806. URL: nature.com/articles/s41598-017-01002-0.pdf.
10. Lee Dong Hyun, Lee Younghwan, Yang Kun, Park Ju Yong, Kim Se Hyun, Reddy Pothala Reddi Sekhar, Materano M., Mulaosmanovic H., Mikolajick T., Jones J. L., Schroeder U., Park Min Hyuk. Applied Physics Reviews, 2021, 8:021312 URL: aip.scitation.org/doi/10.1063/5.0047977.