Известия РАН. Серия физическая, 2021, T. 85, № 8, стр. 1074-1081

Различия в кинетических характеристиках зарядки сегнетоэлектриков и диэлектриков при ионном облучении

К. Е. Озерова 1, А. А. Татаринцев 1, Э. И. Рау 1*, К. Ф. Миннебаев 1, С. В. Зайцев 1

1 Федеральное государственное бюджетное образовательное учреждение высшего образования “Московский государственный университет имени М.В. Ломоносова”
Москва, Россия

* E-mail: rau@phys.msu.ru

Поступила в редакцию 12.03.2021
После доработки 05.04.2021
Принята к публикации 28.04.2021

Полный текст (PDF)

Аннотация

Проведены сравнительные исследования фундаментальных характеристик зарядки классических диэлектриков (ПMMA) и сегнетоэлектриков (LiNbO3) при ионном облучении. Эксперименты проводились при облучении образцов ионами Ar+ с энергией 5–10 кэВ. Показаны как общие закономерности кинетики зарядки, так и типичные различия в характеристиках тока эмитированных ионов, тока смещения (накопления зарядов) и потенциала поверхности. Первой характерной чертой зарядки сегнетоэлектриков является появление дополнительного тока в детектирующих устройствах – тока переполяризации. Второй особенностью является наличие двойного слоя зарядов на поверхности сегнетоэлектриков, существенно влияющего на все кинетические характеристики зарядки.

ВВЕДЕНИЕ

Процессы зарядки диэлектрических материалов при их облучении заряженными частицами вызывают большой интерес в радиационной физике, в космических исследованиях, во многих практических приложениях. Если проблемы зарядки диэлектриков при электронном облучении достаточно широко и полно освещены в литературе, то по ионной зарядке имеется только ограниченное число публикаций [17]. Аналогично выглядит ситуация по радиационной зарядке сегнетоэлектриков. Значительное число работ посвящено электронно-зондовой зарядке сегнетоэлектриков, в частности явлению переполяризации вектора спонтанной поляризации в задачах доменной инженерии, т.е. генерации доменов, зависимостям размеров доменов от доз облучения и энергии электронов [812]. И лишь единичные исследования касаются непосредственно ионно-лучевой записи одно- и двумерных доменов в полярных кристаллах, в частности LiNbO3 [1315].

Но в указанных работах не фигурировал такой важный аспект, как кинетика самой зарядки сегнетоэлектрика, проявленная в таких характеристиках, как генерируемых при ионном облучении поверхностных потенциалах, а также токи эмитированных ионов, токи смещения (зарядки) и токи переполяризации. Эти моменты как раз и являются предметом исследования в настоящей работе. При этом будет учитываться такой важный фактор, как изначальное наличие двойного слоя зарядов на поверхности поляризованного сегнетоэлектрика: зарядов поляризации в подповерхностном объеме и компенсирующих зарядов в адсорбирующем слое на поверхности кристалла [1618].

Известно, что монодоменные сегнетоэлектрики для удовлетворения условия минимума энергии системы должны обладать экранировкой спонтанной поляризации [16]. При этом экранирование для большинства сегнетоэлектриков происходит за счет двух механизмов:

1. “Внутреннее” экранирование за счет носителей заряда – электронов и дырок. Данный механизм экранирования имеет место как для сегнетоэлектриков-полупроводников, так и для чистых диэлектриков.

2. В отсутствие свободных зарядов внутри кристалла (идеальный диэлектрик) экранирование спонтанной поляризации может происходить за счет внешней среды (адсорбция ионов или других заряженных частиц).

Наличие поверхностных электронных уровней может внести существенный или даже определяющий вклад в экранирование, и таким образом для сегнетоэлектриков LiNbO3, LiTaO3 и других справедливы оба механизма экранирования [16]. Изначально на Z-срезе находятся экранирующие заряды: внутренние с эффективной поверхностной плотностью зарядов σi, а также внешние с плотностью σo. Такое распределение зарядов создает в экранирующем объеме толщиной ld поле: Ee = σ0/2εε0 + |Ps – σi|/2εε0, которое будет вызывать поляризацию этого экранирующего слоя с эффективным вектором поляризации Pe. Толщина области пространственного заряда определяется концентрацией свободных носителей заряда ld$n_{0}^{{{{ - 1} \mathord{\left/ {\vphantom {{ - 1} 2}} \right. \kern-0em} 2}}}.$ Согласно [19] дебаевская длина экранировки в системе может быть рассчитана по следующему выражению ${{l}_{D}} = \sqrt {\frac{{kT{{\varepsilon }_{k}}}}{{16\pi {{e}^{2}}{{n}_{s}}}}} \cdot \ln \left( {\frac{{{{n}_{s}}}}{{{{n}_{i}}}}} \right),$ где εk – диэлектрическая проницаемость в точке Кюри, ni – концентрация электронов в квазинейтральной области кристалла, ns – концентрация электронов в области изгиба зон (экранирующий слой). При этом ns может быть определено из формулы ns = = niexp[eΔχ/kT], где eΔχ – величина изгиба зоны. Для ниобата лития в работе [20] приводится значение ns = 1019–1021 см–3, что дает значение lD, порядка 100–400 нм при eΔχ = 0.4 эВ. Таким образом, основная зарядка при облучении ионами Ar+ с энергией до 10 кэВ будет происходить в слое, в котором нет исходного вектора спонтанной поляризации. В этом верхнем слое, согласно [20], будет образовываться экранирующее деполяризующее поле 104–105 В/см, что совпадает с оценками [19].

Наличие такого диэлектрического слоя на поверхности сегнетоэлектрика сказывается на процесс локальной переполяризации вектора спонтанной поляризации и образовании доменов. Так в работе [21] авторы указывают, что максимальное поле, находящееся на расстоянии от доменной стенки, равном толщине диэлектрического слоя, приводит к образованию нанодоменов перед движущейся стенкой [22, 23]. Таким образом, формирование периодической доменной структуры может быть обусловлено существованием эффективного диэлектрического слоя с толщиной, равной расстоянию между краем области пространственного заряда и поверхностью ниобата лития (около 1 мкм, что несколько больше наших оценок) [18]. При этом в [24] указывается, что доменная структура в виде полос образована со средним периодом, равным 110 нм. И именно такое значение может быть связано с толщиной внутреннего диэлектрического слоя в конгруэнтном ниобате лития.

Как отмечается в работе [25] процесс переключения вектора спонтанной поляризации и образования доменов облегчается за счет процессов зарождения (нуклеация) и прорастания множества доменов. В качестве примера можно привести BaTiO3, в котором экспериментальное значение коэрцитивного поля равняется 1 кВ/см, а расчетное 195 кВ/см. При этом стоит подчеркнуть, что процесс зародышеобразования доменов возникает в основном на поверхностных дефектах. Так в статье [26] приводятся данные о том, поле переключения образца ниобата лития, легированного магнием, составило 46 кВ/см, что примерно в пять раз меньше, чем коэрцитивное поле конгруэнтного кристалла ниобата лития и согласуется с результатами работы [27]. На влияние на процесс переполяризации поверхностных дефектов указывалось в исследовании [28], где исходный образец LiNbO3 покрыт тонким диэлектрическим материалом, обладающим очень высокой концентрацией ловушек, что позволяло локализовывать падающие электроны в очень ограниченном объеме слоя. В результате захвата заряда в диэлектрическом слое, находящемся на поверхности сегнетоэлектрического кристалла, возникало сильное электрическое поле, достаточное чтобы вызвать инверсию вектора спонтанной поляризации и образование домена в объеме LiNbO3. Эти домены прорастали на 350 мкм вглубь кристалла. В [28] подмечалось, что при облучении образца электронным пучком сплошной области, возникающий домен состоит из многочисленных нанодоменных областей. Практически мгновенное зарождение большого количества маленьких доменов с плотностью, достигающей 1000 мм–2 после приложения переключающего поля, было показано в более ранней работе [29]. При этом отмечалось, что величина тока переполяризации, регистрируемого в процессе эксперимента пропорциональна площади возникших доменов, и каждый элементарный акт перестройки доменной структуры – слияние доменов или быстрый рост ступеней – сопровождается генерацией элементарного импульса тока переключения. На значительное влияние площади возникших доменов и скорости движения доменной стенки на ток переполяризации указывается в теоретической работе [30].

Таким образом, поле, достаточное для локального переключения вектора спонтанной поляризации и образования доменов меньше, чем коэрцитивное поле материала сегнетоэлектрика. Количество зарождающихся доменов может быть достаточно большое [29], а основой вклад в регистрируемый ток переключения поляризации будет вносить увеличивающаяся эффективная площадь зарожденных доменов в основном за счет срастания доменов. В виду того, что в представленных в [28] экспериментах облучается исходно поверхностный экранирующий слой конгруэнтного сегнетоэлектрика, то влияние исходного среза монодоменного сегнетоэлектрика (–Z-среза и +Z-среза) не оказывает влияние на результирующие токовые и зарядовые характеристики. При этом, будет образовываться большое число неглубоких нанодоменов в области облучения (преимущественно из поверхностных дефектов), а при дальнейшем накоплении заряда они будут расти и объединяться в более крупные. Такое слияние зародышей исходных доменов будет давать основной вклад в полученные токовые характеристики с подложки образца, что указывается в работах [29, 30].

ЭКСПЕРИМЕНТАЛЬНЫЕ РЕЗУЛЬТАТЫ

Экспериментальное исследование процесса зарядки сегнетоэлектриков при ионном облучении проводилось с использованием комплексного метода, опубликованного в работах [6, 7]. В процессе эксперимента исследовались основные токовые характеристики – токи эмиссии Iσ(t), токи смещения, переполяризации и утечки IL + D(t), а также регистрировались энергетические спектры вторичных ионов и их сдвиги по энергетической шкале в процессе зарядки. При этом для предотвращения засева вторичными электронами области положительной зарядки в ионной пушке были использованы три диафрагмы. По сдвигу спектра вторичных ионов определялась зависимость потенциал поверхности от времени облучения VS(t). Методика одновременного измерения основных параметров зарядки диэлектриков реализована в вакуумной камере с источником ионов дугового типа с двойным уплотнением плазмы (дуоплазматрон Арденне). Ускоряющее напряжение ионного пучка Ar+ в процессе эксперимента могло варьироваться от E0 = 1 до 10 кэВ.

В качестве образца для исследования зарядовых характеристик классического диэлектрика был выбран образец полимерного резиста ПММА (полиметилметакрилат). Выбор продиктован тем, что полимерные пленки используется иногда в качестве покрытия при создании двухмерных решеток доменов на поляризованных сегнетоэлектриках [18]. Размеры массивного образца 10 × 10 × 0.5 мм, диэлектрическая константа εr ≈ 2.6, удельное сопротивление ρr = 1013 Ом ⋅ см, плотность ρm = 1.2 г ⋅ см3. Результаты экспериментов представлены на рис. 1. На нем приведены временные зависимости потенциала заряжающей поверхности VS, ток эмиттированных положительных ионов Iσ, интегральный ток с подложки IL +d. Необычное поведение коэффициента эмиссионной характеристики σ ≈ 1 может быть вызвано, как предположено в [7], двумя причинами. Во-первых, усилением коэффициента упругого отражения ионов при росте VS, т.е. при уменьшении энергии падающих ионов EL = E0qVS по мере зарядки поверхности. Второй возможной причиной может быть ионизация вторичных нейтральных распыленных атомов вещества в поле зарядов над поверхностью образца.

Рис. 1.

Кинетические характеристики зарядки полимерного диэлектрика ПММА под воздействием ионного облучения: потенциала поверхности Vs, токов утечки и смещения IL + d, тока эмиттированных ионов Iσ. Аккумулированный заряд Q(t) представлен двумя характеристиками.

Для дальнейшего изложения нам потребуется знание величины аккумулируемого заряда Q(t). Это значение находится из экспериментальных графиков следующим образом. Примем, что в равновесном состоянии насыщения процесса зарядки ток утечки IL = 0. Тогда накопление заряда в образце находится при интегрировании тока смещения Id: $Q = \int_0^{{{t}_{0}}} {{{I}_{d}}(t)dt} .$ По второму варианту Q(t) находится из экспериментального графика VS(t) по соотношению для равномерно заряженной плоскости в свободном пространстве: Vs = = Q0εra, где a – линейный размер стороны заряженного участка, ε0, εr – диэлектрические константы вакуума и материала диэлектрика соответственно. Оба вычисленных значения Q приведены на рис. 2а. Как видно, они довольно хорошо совпадают друг с другом, что говорит о правильности выбранного алгоритма определения зарядов Q(t).

Рис. 2.

Экспериментально измеренные зависимости потенциала Vs(t), токов с полусферического коллектора I1 и I2 с подложки образца-сегнетоэлектрика LiNbO3 +Z-среза (а). Рассчитанные значения аккумулированного заряда Q, токов смещения Id и сумма тока утечки и поляризации IL + p (б).

В качестве сравнительного образца классического сегнетоэлектрика нами выбран кристалл LiNbO3 +Z-среза. Размеры кристалла 10 × 10 × × 0.5 мм, εr = 45, ρm = 4.62 г ⋅ см3, ρr = 1013–1014 Ом ⋅ см, поляризации PS = 10–4–10–5 Кл ⋅ см–2, коэрцитивное поле Fc = 2 ⋅ 105 В ⋅ см–1 [20]. Облучение проводилось ионами Ar+ энергией E0 = = 6 кэВ, током I0 = 0.1 нА, площадь облучения 10–2 см2. Экспериментальные зависимости поверхностного потенциала VS(t), тока I1 с полусферического коллектора и тока I2 с подложки образца представлены на рис. 2а. Видны разительные отличия представленных характеристик от данных на рис. 1. Причиной различия является тот факт, что детектируемый экспериментально сигнал I1 не является для сегнетоэлектрика чистым током эмиттированных ионов Iσ, а ток I2 с подложки не является в чистом виде током смещения Id. В данном случае имеем суперпозицию токов I1 = Iσ + Ip, I2 = Id + IL + Ip, где Ip есть ток поляризации-деполяризации участка, облучаемого ионами. Чтобы выделить вклады токов Iσ, Id и IL + Ip по отдельности, поступим следующим образом. Из экспериментального графика VS(t) находим зависимость Q(t) = VS(t)C(t) = VS(ta2(t0εr/h, где h – толщина образца. Площадь облучаемого образца несколько увеличивается за счет дрейфа зарядов по поверхности, т.е. за счет роста параметра a(t). Затем определяем ток смещения (аккумуляции) зарядов Id = ∂Q(t)/∂t, и токи утечки и переполяризации IL + P = I2Id. Все рассчитанные таким образом характеристики представлены на рис. 2б. Таким образом, при трактовке экспериментов следует подходить с осторожностью при определении взаимосвязанных и взаимозависимых характеристиках зарядки сегнетоэлектриков, особенно при определении поглощенных доз облучения потоками заряженных частиц.

МОДЕЛЬ ЗАРЯДКИ ПРИ ИОННОМ ОБЛУЧЕНИИ И ОБЪЯСНЕНИЕ ПОЛУЧЕННЫХ РЕЗУЛЬТАТОВ

Для построения модели зарядки и переполяризации сегнетоэлектрика ниобата лития при ионном облучении рассмотрим область, в которой происходит зарядка. Так как глубина проникновения ионов Ar+ при энергиях ≤10 кэВ, рассчитанная с использованием SRIM, менее 10 нм, то зарядка осуществляется только в этом верхнем слое.

Исходно сегнетоэлектрик содержит три поверхностные плотности заряда: внутреннюю плотность зарядов σi экранирования, внешнюю с плотностью σo (за счет адсорбированных частиц), а также заряд, обусловленный спонтанной поляризации PS (рис. 3а). Эта распределенная система зарядов вызывает внутреннее поле величиной порядка 104–105 В/см [19, 20]. Для энергии первичных ионов Ar+E0 < 10 кэВ глубина проникновения будет меньше, чем глубина слоя экранировки R < lD. Таким образом, зарядка будет происходить в целом в области сегнетоэлектрика, где вектор спонтанной поляризации PS будет экранирован. В начальный момент не будет значительного ухода внешних экранирующих зарядов при увеличении потенциала зарядки (рис. 3б). При дальнейшем накоплении заряда положительная подрешетка сегнетоэлектрического кристалла будет смещаться дальше от области накопления заряда в объеме облучения, а отрицательная подрешетка будет наоборот, смещаться ближе к области облучения. При этом может зарождаться большая группа доменов, способная прорасти до противоположной стороны, как по отдельности, так и в результате их срастания. Хотя, ожидается, что здесь также преобладает срастание доменов, в результате которого будут генерироваться импульсы тока, объединяющиеся в общий ток переключения [29, 30] (см. рис 2). При этом вклад тока поляризации в процессе срастания нанодоменов в измеряемые токовые характеристики становиться нулевым, и соответствующие токовые характеристики принимают вид, близкий к классическому диэлектрику [6, 7]. Прорастание же домена вглубь не вносит существенного вклада в токовые характеристики главным образом из-за того, что скорость движения доменных стенок вглубь значительно меньше, чем скорость доменных стенок в процессе срастания доменов [29]. В дальнейшем, большое число нанодоменов может слиться в один домен (рис. 3в). При дальнейшем накоплении заряда такой домен может прорастать значительно глубже, чем область толщины слоя, экранирующего вектор спонтанной поляризации (см. рис. 3г). Здесь возникает доменная стенка “голова-к-голове”, способная перемещаться в сильном поле [31]. И таким образом конечная глубины прорастания домена будет определяться итоговым полем, созданным в результате положительной зарядки сегнетоэлектрика.

Рис. 3.

Схема зарядки +Z среза LiNbO3 и зарождения доменов при поляризации: исходный +Z-срез (а), начало облучения (б), зарождение домена (в), этап прорастания домена (г).

Что касается нахождения распределения поля в облучаемой зоне сегнетоэлектрика, то следует учитывать двухслойное распределение плотности зарядов ρ в приповерхностной области. Так для +Z-среза LiNbO3 у поверхности чистого образца расположен слой положительных зарядов, а над поверхностью адсорбируется экранирующий слой ионов с отрицательным зарядом. Для определения электростатических полей E в таких двухслойных (в том числе дипольных) системах обычно применяется уравнение Пуассона:

$\frac{{dE\left( y \right)}}{{dy}} = \frac{{\rho \left( y \right)}}{{{{\varepsilon }_{0}}{{\varepsilon }_{r}}}};~\,\,\,E\left( y \right) = {{\frac{1}{{{{\varepsilon }_{0}}\varepsilon }}}_{r}}\int {\rho \left( y \right)} dy,$
где ρ(y) – нормированное распределение плотности зарядов по глубине x в образце. В настоящей работе примем следующую аппроксимацию распределения ρ(y):
$~\frac{{\rho \left( y \right)}}{{{{\rho }_{0}}}} = \left( {A - B{{y}^{2}}} \right)\exp \left[ { - B{{{({{y - 1} \mathord{\left/ {\vphantom {{y - 1} B}} \right. \kern-0em} B})}}^{2}}} \right],$
где ρ0 – плотность зарядов на поверхности, а за x0 принимаем значение глубины, где происходит смена знаков зарядов, то есть y = x/x0.

Решение уравнения Пуассона приводит к следующему выражению для распределения поля, нормированного на его значение на поверхности:

$\begin{gathered} \frac{{~~E\left( y \right)}}{{{{E}_{0}}}} = \frac{{\left( {By + 1} \right)\exp \left[ { - {{{{{\left( {By - 1} \right)}}^{2}}} \mathord{\left/ {\vphantom {{{{{\left( {By - 1} \right)}}^{2}}} B}} \right. \kern-0em} B}} \right]}}{{2B}} - \\ - \,\,\frac{{~\sqrt \pi \left[ {\left( {2A - 1} \right)B - 2} \right]{\text{erf}}\left[ {{{\left( {1 - By} \right)} \mathord{\left/ {\vphantom {{\left( {1 - By} \right)} {\sqrt B }}} \right. \kern-0em} {\sqrt B }}} \right]}}{{4{{B}^{{3/2}}}}}, \\ \end{gathered} $
где A и B – подгоночные параметры, зависящие от материала образца и параметра x0. На рис. 4 приводятся в качестве примера распределения ρ/ρ0 и E/E0 при параметрах A = 0.5, B = 3.33. Они качественно довольно хорошо совпадают с аналогичными результатами, приведенными в [32, 33].

Рис. 4.

Распределение по глубине образца плотности зарядов ρ/ρ0 и электрического поля E/E0 в двухслойной модели заряженных сегнетоэлектриков при параметрах распределения A = 0.5, B = 3.33.

ЗАКЛЮЧЕНИЕ

На основе представленных сравнительных результатов комплексного исследования зарядки диэлектрика ПММА и +Z-среза сегнетоэлектрика LiNbO3 при их облучении пучком ионов Ar+ следуют следующие выводы.

1. Выявлено большое различие в кинетических характеристиках зарядки диэлектриков и сегнетоэлектриков. Оно обусловлено тем, что в процессе зарядки сегнетоэлектрика добавляется существенная компонента токов и зарядов поляризации-деполяризации, отсутствующих для классических диэлектриков.

2. По результатам экспериментальных характеристик зарядки сегнетоэлектрика предложено качественное описание процесса образование локальных доменов в +Z-срезе LiNbO3 под влиянием ионного облучения.

3. Показан алгоритм нахождения истинных характеристик тока эмиттированных ионов, аккумулируемого заряда и тока переполяризации на основе экспериментально измеренных значений поверхностного потенциала и токов с двух детекторов – полусферы и подложки.

4. Предложено аппроксимированное выражение для распределения плотности зарядов в двухслойной модели зарядов, на основе которого из уравнения Пуассона находится распределение электрического поля в поляризованном сегнетоэлектрике.

Работа выполнена в рамках проекта РФФИ № 18-02-00813а. Авторы благодарны профессору Волк Т.Р. за полезные консультации.

Список литературы

  1. Nagatomi T., Kuwayama T., Toshino K. et al. // J. Appl. Phys. 2009. V. 106. Art. No. 104912.

  2. Shi J., Fama M., Teolis B., Baragiola R. // Nucl. Instrum. Meth. Phys. Res. B. 2010. V. 268. P. 2888.

  3. Croccolo F., Ricardi C. // J. Microscopy. 2008. V. 229. P. 39.

  4. Zoory M.J., Mohammed I.A., Alwan B.J. // J. Chem. Biol. Phys. Sci. 2015. V. 5. P. 2822.

  5. Yogev S., Levin J., Molotskii M. et al. // J. Appl. Phys. 2008. V. 103. Art. No. 064107.

  6. Миннебаев К.Ф., Рау Э.И., Татаринцев А.А. // ФТТ. 2019. Т. 61. С. 1090; Minnebaev K.F., Rau E.I., Tatarintsev A.A. // Phys. Sol. State. 2019. V. 6. No. 6. P. 1013.

  7. Rau E.I., Tatarintsev A.A., Zykova E.Yu. et al. // Vacuum. 2020. V. 177. Art. No. 109373.

  8. Volk T.R., Kokhanchik L.S., Gainutdinov R.V. et al. // J. Adv. Dielectr. 2018. V. 8. Art. No. 183001.

  9. Kokhanchik L.S., Volk T.R. // Appl. Phys. Lett. 2013. V. 110. P. 367.

  10. Kokhanchik L.S., Borodin M.V., Burimov N.I. et al. // Ferroelectrics. 2011. V. 411. P. 71.

  11. Shur V.Ya., Chezganov D.S., Smirnov M.M. et al. // Appl. Phys. Lett. 2014. V. 105. Art. No. 052908.

  12. Tatarintsev A.A., Markovets K.E., Rau E.I. // J. Phys. D. 2019. V. 52. Art. No. 115104.

  13. Li X., Terabe K., Hatano H. et al. // Japan J. Appl. Phys. 2005. V. 44. Art. No. L-1550.

  14. Chezganov D.S., Shur V.Ya., Vlasov E.O. et al. // Appl. Phys. Lett. 2017. V. 110. Art. No. 082903.

  15. Chezganov D.S., Vlasov E.O., Pashnina E.A. et al. // Ferroelectrics. 2020. V. 559. P. 66.

  16. Фридкин В.М. Сегнетоэлектрики-полупроводники. Москва: Наука, 1976. 406 с.

  17. Kalinin S.V., Bonnell D. // Phys. Rev. B. 2001. V. 63. Art. No. 125411.

  18. Shur V.Ya., Chezganov D.S., Akhmatkhanov A.R., Kuznetsov D.K. // Appl. Phys. Lett. 2015. V. 106. No. 23. Art. No. 232902.

  19. Кузьминов Ю.С. Электрооптический и нелинейнооптический кристалл ниобата лития. М.: Наука, 1987. 126 с.

  20. Volk T., Wohlecke M. Lithium niobate. Defects, photorefraction and ferroelectric switching. Berlin: Springer, 2008.

  21. Vlasov E., Chezganov D., Chuvakova M., Shur V.Ya. // Hindawi Scanning. 2018. V. 2018. Art. No. 780982.

  22. Shur V.Y. Nucleation theory and applications. Wiley-VCH Verlag, 2005. P. 178.

  23. Shur V.Ya. // J. Mater. Sci. 2006. V. 41. No. 1. P. 199.

  24. Shur V.Ya., Shishkin E., Rumyantsev E. et al. // Ferroelectrics. 2004. V. 304. No. 1. P. 111.

  25. Барфут Дж., Тейлор Дж. Полярные диэлектрики и их применения. М.: Мир, 1981. С. 259.

  26. Chen Y., Yan W., Guo J. et al. // Appl. Phys. Lett. 2005. V. 87. Art. No. 212904.

  27. Mizuuchi K., Morikawa A., Sugita T., Yamamoto K. // J. Appl. Phys. 2004. V. 96. P. 6585.

  28. Glickman Y., Winebrand E., Arie A., Rosenman G. // Appl. Phys. Lett. 2006. V. 88. Art. No. 011103.

  29. Шур В.Я., Николаев Е.В., Шишкин Е.И. и др. // ФТТ. 2002. Т. 44. № 11. С. 2055; Shur V.Ya., Nikolaev E.V., Shishkin E.I. et al. // Phys. Sol. St. 2002. V. 44. No. 11. P. 2151.

  30. Масловская А.Г., Барабаш Т.К. // Поверхн. Рентг., синхротр., нейтрон. иссл. 2012. № 1. С. 42; Maslovskaya A.G., Barabash T.K. // J. Surf. Invest. X-Ray, Synchrotron. Neutron. Tech. 2012. V. 6. No. 1. P. 37.

  31. Eliseev E. A., Morozovska A. N., Svechnikov G.S. et al. // Phys. Rev. B. 2011. V. 83. № 23. Art. No. 235313.

  32. Meyza X., Goeuriot D., Guerret-Piécourt C. // J. Appl. Phys. 2003. V. 94. P. 5384.

  33. Cazaux J. // J. Appl. Phys. 2004. V. 95. P. 731.

Дополнительные материалы отсутствуют.