В. В. Маслов, Т. А. Трюхан, С. В. Барышников
УШИРЕНИЕ ФАЗОВОГО ПЕРЕХОДА ДЛЯ ТРИГЛИЦИНСУЛЬФАТА В ПОРИСТЫХ МАТРИЦАХ
Исследованы диэлектрические свойства триглицинсульфата в различных пористых матрицах с размерами пор от 4 мкм до 50 А. Обнаружено, что степень размытия сегнетоэлектрического фазового перехода зависит от типа пористой матрицы. Предложено объяснение на основе феноменологической теории сегнетоэлектричества.
Ключевые слова: сегнетоэлектрики, диэлектрическая проницаемость, размерные эффекты.
V. Maslov, T. Tryuhan, S. Baryshnikov
PHASE TRANSITION BROADENING FOR TRIGLITSINSULFAT IN POROUS MATRICES
Dielectric properties of TGS in different porous matrixes with pore sizes from 4 mkm to 5,1 nm were examined and compared with the results obtained for their bulk counterparts. It was found out that the broadening offerroelectric phase transition depends on the type of porous matrix. The results are explained on the basis of Landau—Ginzburg—Devonshire theory.
Keywords: ferroelectrics, dielectric permittivity, size effects.
Многочисленные исследования, проведенные для малых частиц, выявили зависимость их физических свойств от размеров. Одним из способов получения малых частиц является внедрение исследуемого вещества в пористые матрицы, характерный размер пор которых лежит в микрометровом или нанометровом диапазоне. Физические свойства малых частиц в таких матрицах связаны с размерами и геометрией сетки пор. Кроме того, существенную роль играют степень заполнения пористой матрицы, взаимодействие частиц со стенками матрицы и связь частиц между собой. В совокупности эти факторы приводят к тому, что характеристики частиц в ограниченной геометрии могут значительно отличаться от характеристик как соответствующих объемных материалов, так и изолированных малых частиц.
Сегнетоэлектрические свойства частиц, введенных в нанопористые матрицы, изучены сравнительно слабо. Наибольшее количество публикаций посвящено исследованиям малых частиц нитрита натрия, в порах синтетических опалов, пористых стекол и молекулярных решеток MCM-41 (см. работы [6, c. 177; 3, c. 1430] и ссылки в них). Имеются несколько работ, посвященных сег-нетовой соли и триглицинсульфата (ТГС). В частности, в работе [3, c. 1430] было обнаружено, что для пористых пленок Al2O3 с включениями ТГС в зависи-
мости емкости от температуры наблюдаются два максимума. Первый из них примерно соответствует температурной области перестройки доменной структуры (39-41 оС), второй, основной максимум, смещен относительно точки Кюри объемного монокристалла ТГС на 5-6 оС в сторону высоких температур. Сведения по исследованию ТГС в пористых матрицах другого типа в литературе отсутствуют, поэтому целью нашей работы было исследование свойств ТГС в зависимости от типа матрицы.
Триглицинсульфат представляет собой классический сегнетоэлектрик второго рода и уже в течение долгого времени является объектом активных теоретических и экспериментальных исследований. Выше температуры Кюри (Тс = 49 °С) кристалл ТГС имеет моноклинную симметрию и принадлежит к центросимметричному классу 2/m. Ниже Тс зеркальная плоскость исчезает и кристалл принадлежит к полярной точечной группе 2 моноклинной системы. Полярная ось лежит вдоль моноклинной (2-го порядка) оси b. Параметры решетки при комнатной температуре равны: а = 9,15 Ä, b = 12,69 Ä, с = 5,73 Ä. На элементарную ячейку приходится более 100 атомов. Структура ТГС сложна и представляет собой сетку молекул глицина CH2NH2COOH и тетраэдров SO4, связанных между собой водородными связями типа — О-Н...О и N-Н.
В качестве пористых матриц в наших исследованиях использовались: фильтровальная бумага с размером пор порядка 4000 нм (4-103 Ä), окисные пленки Al2O3, полученные путем анодирования химически чистого алюминия в щавелевой кислоте со средним размером пор 100 нм (103 Ä), и силикатные матрицы SBA-15 с размером пор 5,1 нм (51 Ä). Внедрение триглицинсульфата в поры проводилось из насыщенного водного раствора. Для сравнения проводились также исследования соответствующих объемных образцов.
Для измерения электрических параметров образцов использовался цифровой измеритель импеданса E7-12 с рабочей частотой 106 Гц. В качестве электродов применялась In-Ga паста. Для измерения температуры применялся электронный термометр CENTER-304 с хромель-алюмелевой термопарой. Точность измерения температуры составляла 0,1 °С. Исследования проводились в температурном интервале от 20 °С до 100 °С. Для исключения влияния адсорбированной воды измерения проводились в вакууме.
Проведенные исследования показали, что максимум диэлектрической проницаемости для всех образцов не отличается более чем на 1-2 оС от температуры фазового перехода для объемного поликристаллического образца (рис. 1). Во всех матричных образцах наблюдается размытие максимума s'(T), которое определяется не размерами пор, а степенью упорядочения матрицы. Наибольшее размытие s'(T) имеют образцы из фильтровальной бумаги, где — наибольший разброс по размерам и направлениям пор. Меньшую ширину фазового перехода имеют образцы из окисной пленки Al2O3, и еще меньшая ширина s'(T) — для силикатных матриц SBA-15 где размеры пор калиброваны с точностью до 0,1 Ä.
Для образцов ТГС в силикатных матрицах SBA-15 кроме уширения фазового перехода наблюдается активационный рост s'(T) с ростом температуры (рис. 1), что можно объяснить ионно-миграционной поляризацией за счет увеличения числа ионов, принадлежащих поверхности.
1,2
--О-■1
Рис. 1. Температурная зависимость относительного изменения диэлектрической проницаемости для ТГС в различных матрицах:
1— фильтровальная бумага (~ 4 мкм), 2 — пористая пленка Al2O3 (~ 100 нм), 3 — силикатная матрица SBA-15 (—5,1 нм), 4 — прессованный объемный ТГС
20 30 40 50 60 70
1, "С
80
90 ЮС
Размытие фазового перехода является общим свойством всех сегнетоэлек-трических кристаллов при увеличении в них неоднородности поляризации, причиной которой могут являться: доменные границы, зерна в керамике, поверхностные слои, тепловые флуктуации, а также дефекты и примеси. Это размытие означает, что в сильно неоднородных кристаллах сингулярность в поведении материальных констант исчезает. По-видимому, это связано с тем, что в таких кристаллах фазовый переход при температуре То происходит не по всему объему образца, а появляется некоторое распределение локальных температур перехода, значения которых зависят от распределения неоднородностей и дефектов по различным локальным областям кристалла [4, с. 95].
Размытие фазового перехода можно объяснить и в рамках феноменологической теории Ландау—Гинзбурга [2, с. 490]. Вклад в плотность свободной энергии сегнетоэлектрика, обусловленный поляризацией Р, может быть представлен в виде
где а0, в, у — температурно независимые коэффициенты Ландау, То — температура Кюри. Отметим, что при в > 0, у > 0 выражение (1) описывает фазовый переход второго рода, а при в < 0, у> 0 — переход первого рода.
Разложение Ландау для плотности свободной энергии неоднородной системы содержит не только различные степени параметра порядка, но также и производные от параметра порядка по координатам. В случае идеализированной изотропной системы первые производные параметра порядка могут входить в разложение только в виде скалярной комбинации (§га^)2, и поэтому мы можем в самом общем виде записать [5, с. 197]:
Используя вариационный подход и минимизируя свободную энергию (2), можно найти уравнение состояния в электрическом поле:
Р =1аР2 + - вР4 +1 гР6 - ЕР, а = а (Т - Т),
2 4 6
(1)
Р(((г),Т) = 2а(Т)Р2(г) + -4вР» +13(Т)(§гааР(г))2 - ЕР. (2)
Е = ао (Т - То)Р(г) + вР3(г) + 23(Т) (ДР(г)).
(3)
Чтобы получить диэлектрическую восприимчивость х, нам необходимо найти производную по поляризации от соотношения (3):
Или, ограничиваясь первыми членами разложения и выражая диэлектрическую проницаемость, получим:
Из вышеизложенного следует, что для неоднородных систем в знаменателе закона Кюри—Вейсса появляется дополнительное слагаемое, приводящее к размытию кривой в'(Т). Степень размытия будет определяться неоднородностью композита матрица—сегнетоэлектрик, которую проблематично выразить в общем виде, так как она в каждом конкретном случае зависит от типа матрицы и вида сегнетоэлектрика.
Рассмотрим простейшую задачу, связанную с границей доменов: представим себе, что бесконечный кристалл состоит из двух доменов. Будем считать, что в одном полупространстве (г > 0) Р > 0, а в другом (г < 0) Р < 0, т. е. доменная стенка совпадает с плоскостью ху. В первом приближении спонтанная поляризация Р(г) между этими областями примет форму ступеньки, поскольку при таком распределении в каждой точке кристалла значение параметра порядка равно либо +|Ро|, либо -|Ро|, а только этим двум значениям Р и отвечает минимум термодинамического потенциала. Поэтому наличие областей с другими значениями Р означало бы энергетический проигрыш — увеличение термодинамического потенциала по сравнению с его значением при Р = Ро. В действительности, однако, создав ступеньку, мы отнюдь не избежали энергетического проигрыша; наличие ступеньки означает бесконечную величину (5Р/&), и, стало быть, вклад в плотность термодинамического потенциала члена 0,55(дР/&) бесконечно велик в области ступеньки. Чтобы его уменьшить, нужно «размазать» ступеньку, при этом появятся области, где |Р| ^ |Ро|. В результате то распределение Р, при котором увеличение термодинамического потенциала, обусловленное появлением доменной стенки, будет минимальным, имеет вид «размазанной» ступеньки. Явный вид этого распределения для полидоменного образца можно найти из условия минимальности термодинамического потенциала кристалла:
Отыскание функции Р = Р^), которая минимизирует интеграл (6), — типичная задача вариационного исчисления; известно, что такая функция должна удовлетворять уравнению Эйлера:
Г1 = ^ = ао (Т - То) + 3вР2(г) + 2ЩТ)-Р(АР(г)). аР аР
(4)
е'(Т) * 1 +---------------------------- ,
4 7 ао (Т - То) + /(Р(г, Т)) ’
/ (Р(г), Т) = 23(Т) -Р (АР(г, Т)).
(5)
где
(6)
Отсюда получим нелинейное дифференциальное уравнение второго порядка
aP + вРъ-8
д2 P дz2
= 0.
Рис. 2. Распределение поляризации на доменной границе
Регулярная процедура решения таких уравнений отсутствует. Пример решения такой задачи для плоской доменной стенки можно найти в работе [5, с. 197], где показано, что на границе P(z) удовлетворяет реше-
нию P( z) = + — Ш
радиус корреляции.
В реальных кристаллах зарождение доменов и ориентация Ро в зернах поликристаллов зависит от многих параметров, например: от дефектности кристалла, от соотношения размера зерна и радиуса корреляции, скорости охлаждения, от коорцитивного поля и т. д.
Для понимания изменения температуры фазового перехода при уменьшении размера пор обратимся к теории ГАК (Гинзбурга—Андерсона—Кохрэна), где природа сегнетоэлектричества связывается с существованием в кристалле низкочастотного температурно-зависимого поперечного оптического колебания «мягкой моды».
При сокращении размеров кристалла в одном, в двух или в трех измерениях (пленка, нить, сфера) из-за отсутствия трансляционной симметрии при анализе динамики решеточных колебаний нельзя уже использовать циклические условия Борна—Кармана. Последние должны быть заменены граничными условиями на поверхности. В результате, прежде всего, не реализуются длинные волны с Л/2>ё (д<п/ф, где ё — меньший размер кристалла, а q = 2п/Л — волновой вектор [1, с. 138]. В обычных кристаллах оптические ветви ЬО и ТО почти не имеют дисперсии, и на них такое отсечение почти не влияет. В сегне-тоэлектриках благодаря большой роли дальнодействующих кулоновских сил частота ветви ТО для длинных волн резко падает и размер ё может существенно влиять на спектр решеточных колебаний. Последнее приводит к тому, что при определенных размерах частиц длинноволновые колебания не возникают и сег-нетоэлектрические свойства уже не будут реализовываться.
Строго говоря, приведенные выше выводы получены: во-первых, для изолированных сегнетоэлектрических частиц; во-вторых, для веществ, где «мягкая мода» имеет малое затухание. Для сегнетоэлектриков типа порядок—беспорядок сегнетоэлектрическое колебание имеет релаксационный характер с большим коэффициентом затухания, и говорить о длине волны и граничных условиях не совсем корректно.
Из вышеизложенного можно предположить, что в сегнетоэлектриках типа ТГС механизм, связанный с ограничением длины волны «мягкой моды», не будет иметь места, хотя остаются механизмы, связанные с изменением соотношения количества атомов в объеме и на поверхности, с механическим зажатием нанокристаллов в порах и с взаимодействием отдельных частиц.
Разберем возможность взаимного влияния отдельных сегнетоэлектриче-ских частиц. А для этого рассмотрим массив наночастиц, находящихся в пористой матрице на некотором расстоянии друг от друга. В этом случае свободная энергия Гиббса системы из наночастиц может быть записана в виде суммы энергии частиц и энергии их связи. Для описания фазового перехода в малой изолированной частице обычно используется тот же подход, основанный на разложении свободной энергии Ландау [2, с. 490], где поляризация играет роль параметра порядка. Кроме этого необходимо учесть дополнительный вклад, обусловленный поверхностной энергией с соответствующими граничными условиями и деполяризующее поле.
Таким образом, свободную энергию системы сегнетоэлектрических частиц можно записать в виде суммы объемного и поверхностного интегралов:
где Р; — поляризация, которая является функцией температуры Т и координат, г — является экстрополяционной длиной, которая определяется градиентом поляризации на границе частицы в соответствии со следующим законом:
где п является единичной нормалью к поверхности частицы. В соотношении (7) деполяризующие поля не учитываются и предполагается, что в системе происходит фазовый переход II рода.
Чтобы учесть связь между частицами, мы предположим, что расстояние между сегнетоэлектрическими частицами в несколько раз превышает размеры частиц. С учетом этого энергию диполь-дипольного взаимодействия 1-й и ]-й частиц можно записать в виде [7, с. 75]:
где gi,gj — коэффициенты, учитывающие уменьшение диполь-дипольного взаимодействия за счет зарядов на поверхности частиц. Для упрощения дальнейших расчетов можно неоднородную поляризацию малых частиц заменить эффективной однородной поляризацией р;, полученной путем усреднения по объему. В результате такой замены термодинамический потенциал для сегнето-электрика с фазовым переходом II рода с учетом энергии взаимодействия частиц можно записать в виде
где V — объем 1-й частицы, коэффициент а01 и Уы линейно зависит от температуры а; = ао; [Т-То,(У,)].
Таким образом, свободная энергия системы взаимодействующих малых частиц будет определяться соотношением
(7)
дР __ Р
дп г
(8)
и
(9)
j L i* j
Для большого регулярного массива одинаковых малых частиц задача сводится к сумме равных энергий отдельных частиц и энергии сцепления [7, с. 75]:
Из соотношения (10) следует, что температура фазового перехода системы связанных частиц по сравнению с температурой фазового перехода в изолированных малых частицах будет определяться соотношением
Величина коэффициента К определяется формой, взаимным расположением частиц в матрице и направлением поляризации. В зависимости от этих характеристик знак К может быть отрицательным или положительным.
1. Бурсиан Э. В. Нелинейный кристалл. Титанат бария. М.: Наука, 1974. 295 с.
2. Гинзбург В. Л. Теория сегнетоэлектрических явлений // УФН. 1949. Т. 38. C. 490.
3. Рогазинская О. В., Миловидова С. Д., Сидоркин А. С., Чернышев В. В., Бабичева Н. Г. Свойства нанопористого оксида алюминия с включениями триглицинсульфата и сегнетовой соли // ФТТ. 2009. Т. 51. C. 1430.
4. Струков Б. А. Фазовые переходы в сегнетоэлектрических кристаллах с дефектами // Соросовский образовательный журнал. 1996. № 12. C. 95.
5. Струков Б. А., Леванюк А. П. Физические основы сегнетоэлектрических явлений в кристаллах. М.: Наука, 1995. 302 с.
6. Baryshnikov S. V., Cheng Tien, Charnaya E. V., Lee M. K., Michel D., Bohlmann W., Stukova E. V. Dielectric Properties of Mesoporous Sieves Filled with NaNO2 // Ferroelectrics
2008. 363. P. 177.
7. Charnaya E. V., Pirozerskii A. L., Cheng Tien, Lee M. K. Ferroelectricity in an Array of Electrically Coupled Confined Small Particles // Ferroelectrics. 2007. 350. P. 75.
1. Bursian Je. V. Nelinejnyj kristall. Titanat barija. M.: Nauka, 1974. 295 s.
2. Ginzburg V. L. Teorija segnetoelektricheskih javlenij // UFN. 1949. T. 38. C. 490.
3. Rogazinskaja O. V., Milovidova S. D., Sidorkin A. S., Chernyshev V. V., Babicheva N. G. Svojstva nanoporistogo oksida aluminija s vkluchenijami triglicinsul'fata i segnetovoj soli // FTT.
2009. T. 51. C. 1430.
4. Strukov B. A. Fazovye perehody v segnetoelektricheskih kristallah s defektami // Sorosovskij obrazovatel'nyj zhurnal. 1996. № 12. C. 95.
5. StrukovB. A., LevanukA. P. Fizicheskie osnovy segnetoelektricheskih javlenij v kristallah. M.: Nauka, 1995. 302 s.
6. Baryshnikov S. V., Cheng Tien, Charnaya E. V., Lee M. K., Michel D., Bohlmann W., Stukova E. V. Dielectric Properties of Mesoporous Sieves Filled with NaNO2 // Ferroelectrics. 2008. 363. P. 177.
7. Charnaya E. V., Pirozerskii A. L., Cheng Tien, Lee M. K. Ferroelectricity in an Array of Electrically Coupled Confined Small Particles // Ferroelectrics. 2007. 350. P. 75.
(10)
СПИСОК ЛИТЕРАТУРЫ
REFERENCES