Термоактивационная токовая спектроскопия электрически активных центров в эпитаксиальных монокристаллических пленках ферритов-гранатов (TmBi)3(FeGa)5O12:Ca2+
Аннотация
Настоящая работа посвящена изучению спектров термостимулированного тока проводимости образцов феррит-гранатовых пленок (TmBi)3(FeGa)5O12:Ca2+ с целью выявления в них электрически активных центров. Объекты исследования были получены методом ЖФЭ из раствора в расплаве на подложках Gd3Ga5O12 кристаллографической ориентации <111> при использовании бессвинцового Bi-содержащего флюса. Увеличение концентрации ионов Ca2+ в пленках достигали путем увеличения в растворе-расплаве концентрации CaO. Исследуемые пленки содержали ионы Ca2+ от 0,03 форм. ед. до 0,1 форм. ед. На универсальной установке термоактивационной токовой спектроскопии получены спектры исходных образцов и образцов, возбужденных воздействием ультрафиолетового облучения или обработкой в отрицательном коронном разряде. Проведен расчет параметров электрически активных дефектов и других электрофизических параметров.
Ключевые слова: термостимулированный ток проводимости, пленки ферритов-гранатов, электрически активные центры, глубокие уровни захвата, механизмы зарядовой компенсации, легирование кальцием, ТСТП-спектр, температурное положение
Введение
Для управления свойствами эпитаксиальных монокристаллических пленок ферритов-гранатов (ЭМПФГ) (в частности, для осуществления зарядовой компенсации ионов Me4+) широко используется легирование ионами Ca2+ [1]. При использовании легирования Ca2+ - Me4+ в кристаллической решетке пленки сложно достичь полного равенства концентрации ионов Ca2+ и Me4+, что приводит к реализации различных механизмов зарядовой компенсации и, в результате, к понижению эксплуатационных параметров пленок. Несмотря на наличие в научной литературе огромного количества публикаций по изучению поведения ионов Ca2+ в магнитных гранатах [2], на сегодняшний день нет единого взгляда на природу зарядовой компенсации кальция в номинально чистых кристаллах Y3Fe5O12 и сложнозамещенных оксидах на его основе.
Настоящая работа посвящена изучению спектров термостимулированного тока проводимости (ТСТП) образцов феррит-гранатовых пленок (TmBi)3(FeGa)5O12: Ca2+ с целью выяснения в данных материалах природы электрически активных центров и уточнения реально действующих механизмов зарядовой компенсации.
Объекты исследований и методы эксперимента
В качестве объектов исследования в работе использовались ЭМПФГ системы (TmBi)3(FeGa)5O12: Ca2+.
ЭМПФГ были получены методом жидкофазной эпитаксии (ЖФЭ) из раствора в расплаве на подложках Gd3Ga5O12 кристаллографической ориентации <111> по стандартной технологии при использовании бессвинцового Bi-содержащего флюса. Увеличения концентрации ионов Ca2+ в пленках достигали путем увеличения в растворе-расплаве концентрации CaO.
Использовалась установка жидкофазной эпитаксии фирмы LAS (Франция). Кристаллизация происходила из переохлажденного раствора гранатообразующих компонент в расплаве растворителя Bi2O3.
В процессе выращивания объектов исследования молярные параметры изменяли в следующих пределах: R1: 0,137-0,140; R2: 0,180-0,220; R3: 0,0622-0,0625; R4: 0,149-0,160; R5: 0-0,169. В результате проведенной работы получены ЭМПФГ (TmBi)3(FeGa)5O12 с содержанием ионов Ca2+ от 0,03 формульных единиц (ф.е.) до 0,1 ф.е. Толщина пленок составляла 1,5-33,4 микрон.
В таблице № 1 представлены основные технологические характеристики полученных в работе пленок (TmBi)3(FeGa)5O12: Ca2+.
Таблица № 1
Основные технологические характеристики роста полученных в работе пленок (TmBi)3(FeGa)5O12/Gd3Ga5O12<111>
№ п/п |
Пленка |
ТРОСТ, |
h, |
fРОСТ, мкм/мин |
Значения молярных коэффициентов |
||||
R1 |
R2 |
R3 |
R4 |
R5 |
|||||
1 |
23 |
850 |
1,5 |
0,25 |
0,140 |
0,180 |
0,0625 |
0,149 |
0 |
2 |
26 |
839 |
30,35 |
0,51 |
0,140 |
0,180 |
0,0625 |
0,149 |
0 |
3 |
28 |
836,4 |
2,69 |
0,53 |
0,139 |
0,180 |
0,0625 |
0,151 |
0,028 |
4 |
29 |
839 |
27,3 |
0,40 |
0,139 |
0,180 |
0,0625 |
0,151 |
0,028 |
5 |
30 |
839,5 |
1,62 |
0,324 |
0,139 |
0,180 |
0,0624 |
0,154 |
0,056 |
6 |
31 |
840 |
23,9 |
0,27 |
0,138 |
0,180 |
0,0624 |
0,153 |
0,056 |
7 |
33 |
836,7 |
22,81 |
0,228 |
0,138 |
0,220 |
0,0624 |
0,153 |
0,084 |
8 |
34 |
837 |
5,68 |
0,258 |
0,138 |
0,220 |
0,0624 |
0,153 |
0,084 |
9 |
35 |
835 |
33,4 |
0,24 |
0,138 |
0,220 |
0,0624 |
0,154 |
0,096 |
10 |
36 |
833 |
6,18 |
0,206 |
0,138 |
0,220 |
0,0623 |
0,155 |
0,112 |
11 |
37 |
831 |
23,7 |
0,16 |
0,137 |
0,220 |
0,0623 |
0,155 |
0,112 |
12 |
38 |
829 |
4,56 |
0,228 |
0,137 |
0,220 |
0,0623 |
0,158 |
0,141 |
13 |
39 |
827 |
31,7 |
0,22 |
0,137 |
0,220 |
0,0623 |
0,158 |
0,141 |
14 |
40 |
828,7 |
5,3 |
0,176 |
0,137 |
0,220 |
0,0622 |
0,160 |
0,169 |
Примечание: Ниже приведены молярные коэффициенты, характеризующие соотношение компонент шихты и их нумерация:
(1)
(2)
(3)
(4)
(5)
Уточнение состава ЭМПФГ проводилось на микроанализаторе «JEOL-JXA» путем рентгеноспектрального анализа пленок. Для предотвращения интенсивного накопления образцами заряда во время анализа, на поверхность последних напыляли тонкую пленку золота толщиной ~500 Å.
В таблице № 2 представлены результаты рентгеноструктурного анализа и рентгеноспектрального микроанализа полученных в настоящей работе пленок (TmBi)3(FeGa)5O12: Ca2+.
Таблица № 2
Результаты рентгеноструктурных и рентгеноспектральных исследований полученных в работе пленок (TmBi)3(FeGa)5O12: Ca2+
№ п/п |
Пленка |
Δа, |
Са2+, ф.е. |
Tm3+, ф.е. |
Bi3+, ф.е. |
Fe3+, ф.е. |
Ga3+, ф.е. |
1 |
23 |
- 27 |
0 |
2,250 |
0,712 |
3,534 |
1,346 |
2 |
26 |
- 1,8 |
0 |
2,250 |
0,712 |
3,534 |
1,346 |
3 |
28 |
+12,3 |
0,034 |
2,222 |
0,748 |
3,534 |
1,348 |
4 |
29 |
+ 4,9 |
0,034 |
2,222 |
0,748 |
3,534 |
1,348 |
5 |
30 |
+ 0,3 |
0,052 |
2,184 |
0,700 |
3,482 |
1,314 |
6 |
31 |
+ 1,2 |
0,052 |
2,184 |
0,700 |
3,482 |
1,314 |
7 |
33 |
- 1,2 |
0,068 |
2,182 |
0,722 |
3,326 |
1,500 |
8 |
34 |
- 1,4 |
0,068 |
2,182 |
0,722 |
3,326 |
1,500 |
9 |
35 |
+ 0,5 |
0,07 |
2,108 |
0,744 |
3,320 |
1,496 |
10 |
36 |
- 3,4 |
0,084 |
2,108 |
0,744 |
3,320 |
1,496 |
11 |
37 |
+ 6,4 |
0,084 |
2,128 |
0,742 |
3,350 |
1,498 |
12 |
38 |
+ 6,8 |
0,088 |
2,128 |
0,742 |
3,350 |
1,498 |
13 |
39 |
+ 9,2 |
0,088 |
2,066 |
0,774 |
3,316 |
1,460 |
14 |
40 |
+ 4,0 |
0,092 |
2,066 |
0,774 |
3,316 |
1,460 |
Анализ представленных в таблице № 2 данных позволил заключить, что для всех используемых в работе ЭМПФГ (за исключением образцов № 23 и № 26) характерно вхождение в пленку ионов Ca2+, причем концентрация последних растет с увеличением значения коэффициента R5. Так, при увеличении R5 от значения 0,028 до 0,169 концентрация ионов кальция растет со значения 0,034 ф.е. до 0,092 ф.е. соответственно.
Известно, что электрически активные дефекты создают в кристаллической решетке кристалла глубокие уровни захвата [3]. Эффективными методами исследования глубоких уровней захвата в кристаллах являются методы термоактивационной токовой спектроскопии (ТАТС), куда входит и метод ТСТП [3].
Регистрация спектров ТСТП проводилась на разработанной универсальной автоматизированной установке ТАТС в диапазоне 300 – 673 К при скорости линейного нагрева 0,05 – 0,15 К/с [4]. В качестве измерителя малых токов использовался прибор ИМТ-05. Для регистрации сигнала с исследуемого образца применялись прижимные стальные никелированные электроды с наконечниками из серебра. На поверхность исследуемого образца наносились электроды на основе (In,Ga)-эвтектики размерами 5 × 1 мм, расстояние между электродами составляло 1 мм. Наличие омического контакта проверялось перед каждой регистрацией спектра. Возбуждение образцов проводилось в униполярном коронном разряде и в ультрафиолете.
Анализ представленных в таблице № 2 данных позволил заключить, что для всех используемых в работе ЭМПФГ (за исключением образцов № 23 и № 26) характерно вхождение в пленку ионов Ca2+, причем концентрация последних растет с увеличением значения коэффициента R5. Так, при увеличении R5 от значения 0,028 до 0,169 концентрация ионов кальция растет со значения 0,034 ф.е. до 0,092 ф.е. соответственно.
Известно, что электрически активные дефекты создают в кристаллической решетке кристалла глубокие уровни захвата [3]. Эффективными методами исследования глубоких уровней захвата в кристаллах являются методы термоактивационной токовой спектроскопии (ТАТС), куда входит и метод ТСТП [3].
Регистрация спектров ТСТП проводилась на разработанной универсальной автоматизированной установке ТАТС в диапазоне 300 – 673 К при скорости линейного нагрева 0,05 – 0,15 К/с [4]. В качестве измерителя малых токов использовался прибор ИМТ-05. Для регистрации сигнала с исследуемого образца применялись прижимные стальные никелированные электроды с наконечниками из серебра. На поверхность исследуемого образца наносились электроды на основе (In,Ga)-эвтектики размерами 5 × 1 мм, расстояние между электродами составляло 1 мм. Наличие омического контакта проверялось перед каждой регистрацией спектра. Возбуждение образцов проводилось в униполярном коронном разряде и в ультрафиолете.
Результаты экспериментальных исследований и их обсуждение
Для исследования энергетического спектра структурных дефектов ЭМПФГ вначале регистрировались ТСТП-спектры исследуемых материалов в исходном состоянии, а затем регистрировались спектры образцов, возбужденных воздействием отрицательного коронного разряда или УФ-облучения в течение времени 5 - 20 минут. По пикам ТСТП проводился расчет энергии активации и других электрофизических параметров.
При измерениях термостимулированных токов для каждого значения времени обработки в отрицательном коронном разряде или после УФ-облучения были получены зависимости тока, протекающего в образце, от температуры.
Спектры регистрировались в режиме ТСТП и ТСТП КЗ (короткого замыкания), термостимулированную деполяризацию зарегистрировать не удалось.
Следует отметить, что в режиме ТСТП удавалось получить спектры только при незначительном напряжении измерения (0,5-5,0 мВ). При превышении измерительным напряжением указанного значения спектр ТСТП представлял собою экспоненциальную зависимость тока от температуры, что обусловлено интенсивным нарастанием «темнового» тока, существенно превышающего по значению ток, обусловленный высвобождением носителей с уровней захвата.
Способом определения энергии активации центров захвата по кривой «термовысвечивания» явился способ температурного положения пика в соответствии с методиками Урбаха и Бьюба [3].
С целью получения тонкой структуры спектра, разделения пиков с близкими значениями температурного положения, а также с целью избавления от «темнового» тока, «забивающего» ток глубоких уровней захвата, был применен метод ТСТП в режиме короткого замыкания (ТСТП КЗ).
Следует отметить, что на сегодняшний день в научной литературе имеется всего несколько работ по изучению магнитных гранатов (в том числе и ЭМПФГ) методами ТАТС [5-6].
На рис. 1 представлены характерные спектры ТСТП и ТСТП КЗ для пленки 38. Всего для данного состава ЭМПФГ нами было обнаружено 36 локальных электрически активных центров.
а б
в г
д е
Рис. 1. – Характерные спектры ТСТП и ТСТП КЗ для эпитаксиальных пленок ферритов-гранатов (TmBi)3(FeGa)5O12: Ca2+
Все данные о температурном положении пиков для спектров ТСТП всех пленок представлены в таблице № 3.
Таблица № 3
Температурное положение пиков* ТСТП КЗ эпитаксиальных пленок ферритов-гранатов (TmBi)3(FeGa)5O12: Ca2+
* - Возбуждение образцов: отрицательный коронный разряд, ультрафиолет
Температура, °C |
40-70 |
70-100 |
100-130 |
130-160 |
160-190 |
190-220 |
220-250 |
250-280 |
280-310 |
310-340 |
340-370 |
370-400 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
12 |
13 |
Время обработки, вид обработки |
|
|
|
|
|
|
|
|
|
|
|
|
29 |
||||||||||||
исходный |
|
82,1 |
|
|
|
|
|
|
|
324,8 |
|
|
5 мин – корона, |
|
|
|
140,4 |
|
|
226,3 |
|
|
|
|
|
5 мин – УФ, |
|
|
|
|
|
|
|
256,6 |
|
|
|
|
37 |
||||||||||||
5 мин – УФ, |
|
|
|
|
|
|
|
|
|
317,1 |
|
|
10 мин – корона, |
|
90,2 |
|
|
|
|
|
|
|
|
340,9 |
|
5 мин – УФ |
|
|
|
139,7 |
|
|
|
|
|
|
|
|
38 |
||||||||||||
исходный |
|
|
|
|
|
|
|
|
301 |
|
|
|
10 мин – корона, |
|
76,3 |
|
|
|
|
247,3 |
|
|
|
|
|
5 мин – корона, |
|
|
|
|
|
209,3 |
|
|
|
|
|
|
20 мин - корона |
|
|
|
|
|
216 |
|
|
|
|
|
|
10 мин - корона |
|
|
|
|
|
195,4 |
|
|
302,9 |
|
|
|
10 мин – УФ, |
|
|
|
|
|
|
|
|
|
335,5 |
|
|
39 |
||||||||||||
исходный |
|
|
115,3 |
|
|
|
|
|
|
313,9 |
|
|
40 |
||||||||||||
исходный, |
|
|
|
|
|
|
|
276,9 |
|
314,5 |
|
|
5 мин – УФ, |
|
|
|
|
189,3 |
|
|
|
|
|
342,2 |
|
5 мин – корона, |
|
|
|
150 |
|
|
|
|
281,1 |
|
|
|
10 мин – корона, |
|
|
|
|
|
|
234,4 |
|
|
|
363,1 |
|
По результатам Урбаха энергия активации центра захвата может быть рассчитана по следующей формуле (6):
(6)
где Et – энергия активации захвата; Тm – температура, соответствующая пику max интенсивности.
По данным Бьюба величина Еt может быть представлена как формула (7):
(7)
где k – константа Больцмана, k=8,617×10−5 эВ·К−1.
Результаты расчетов энергии активации высвобождения заряда по приведенным выше формулам представлены в таблице № 4.
Концентрация активных центров для каждого типа центра определялась по величине выделенного центром заряда. Величина выделенного центром заряда, в свою очередь, определялась по площади соответствующего центру пика ТСТП. Расчёт заряда для каждого пика ТСТП определялся по формуле (8):
(8)
где: S – площадь кривой под пиком; MX и MY – масштаб по осям X и Y соответственно; е – заряд электрона.
Концентрация электрически активных центров (для каждого типа центров) рассчитывалась по формуле (9):
(9)
где n – концентрация электрически активных центров, см-3; Q – рассчитанный по площади под кривой выделенный данным типом центров заряд, Кл; e – заряд электрона, Кл; SЭ – площадь электродов, см2; h – толщина пленки, см.
Данные расчета заряда и концентрации электрически активных центров представлены в таблице № 4.
Таблица № 4
Энергия активации ловушек и концентрация электрически активных центров в пленках (TmBi)3(FeGa)5O12: Ca2+ (расчет по пикам ТСТП). Возбуждение: отрицательный коронный разряд; ультрафиолет
Время обработки, мин |
Номер пика |
Температурное положение пика Тm |
Средняя энергия активации |
Заряд* 10-8, Кл |
Суммарный заряд* 10-8, Кл |
Концентрация центров 1016, см-3 |
|
°С |
К |
||||||
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
29 |
|||||||
исходный |
1 |
82,1 |
355,1 |
0,70 |
25,4 |
34,92 |
1,16 |
2 |
324,8 |
597,8 |
1,18 |
9,52 |
0,44 |
||
5 - корона |
1 |
140,4 |
413,4 |
0,81 |
12,1 |
20,14 |
0,55 |
2 |
226,3 |
499,3 |
0,98 |
8,04 |
0,37 |
||
5 - УФ |
1 |
256,6 |
529,6 |
1,04 |
18 |
18 |
0,82 |
37 |
|||||||
5 – УФ |
1 |
317,1 |
590,1 |
1,16 |
106 |
106 |
5,61 |
5 – УФ |
1 |
139,7 |
412,7 |
0,81 |
69,1 |
69,1 |
3,64 |
10 - корона |
1 |
90,2 |
363,2 |
0,72 |
6,72 |
31,42 |
0,04 |
2 |
340,9 |
613,9 |
1,21 |
24,7 |
1,31 |
||
38 |
|||||||
исходный |
1 |
301 |
574 |
1,13 |
45,7 |
45,7 |
12,5 |
20 - корона |
1 |
216 |
489 |
0,96 |
28,5 |
28,5 |
7,82 |
10 - УФ |
1 |
335,5 |
608,5 |
1,20 |
9,05 |
9,05 |
2,48 |
10 – корона (U=1 мВ) |
1 |
195,4 |
468,4 |
0,92 |
12 |
125 |
3,28 |
2 |
302,9 |
575,9 |
1,13 |
113 |
30,9 |
||
10 – корона (U=0 мВ) |
1 |
76,3 |
349,3 |
0,69 |
9,66 |
40,46 |
2,65 |
2 |
247,3 |
520,3 |
1,03 |
30,8 |
8,45 |
||
5 - корона |
1 |
209,3 |
482,3 |
0,95 |
9,02 |
9,02 |
2,47 |
39 |
|||||||
исходный |
1 |
115,3 |
388,3 |
0,77 |
6,51 |
8,97 |
0,26 |
2 |
313,9 |
586,9 |
1,16 |
2,46 |
0,10 |
||
40 |
|||||||
исходный |
1 |
276,9 |
549,9 |
1,08 |
49,1 |
173,1 |
11,6 |
2 |
314,5 |
587,5 |
1,16 |
124 |
29,2 |
||
5 - УФ |
1 |
189,3 |
462,3 |
0,91 |
2,87 |
104,87 |
0,68 |
2 |
342,2 |
615,2 |
1,21 |
102 |
24,1 |
||
5 - корона |
1 |
150 |
423 |
0,83 |
6,46 |
43,86 |
0,15 |
2 |
281,1 |
554,1 |
1,09 |
37,4 |
8,81 |
||
10 - корона |
1 |
234,4 |
507,4 |
1,00 |
25,5 |
56,7 |
6,01 |
2 |
363,1 |
636,1 |
1,25 |
31,2 |
7,36 |
* - количество выделенного данным типом центра заряда
Связывать обнаруженные пики ТСТП непосредственно с ионами Ca2+ нет оснований, поскольку последние, по литературным данным, создают в решетке ЖИГ акцепторный уровень, который «высвечивается» при температурах ниже комнатной. Поэтому, обнаруженные пики следует связывать с индуцируемыми ионами кальция дефектами. Такие дефекты создают в запрещенной зоне объектов исследования электрически активные центры, проявляющиеся на спектрах ТСТП в виде пиков тока.
Ионы Са2+ являются электроотрицательными по отношению к катионной подрешетке, вследствие чего для соблюдения электронейтральности должен работать один (или несколько) из следующих механизмов зарядовой компенсации:
1) Переход части ионов железа в состояние Fe4+;
2) Появление дырочных центров О-;
3) Появление кислородных вакансий V(O2-)2+; V(O2-)+ (F+-центр) и V(O2-)0 (F-центр).
Используя данные оптических исследований [1], можно сделать вывод, что ионы Fe4+ должны создавать в магнитных гранатах уровень захвата с максимумом в области 0,2 эВ, что соответствует температуре порядка 100 К.
Обнаруженные же нами центры имеют минимальное значение температурного положения 349-355 К. Поэтому, обнаруженные центры следует связывать с дефектами по кислороду и реализацией 2-го и 3-го из вышеуказанных механизмов зарядовой компенсации.
Отметим, что наличие большого количества активных центров – весьма характерно для кристаллов со структурой граната и ранее было обнаружено для немагнитных гранатов [7]. Такое многообразие обнаруженных центров следует связывать со сложной структурой граната. Тот факт, что, скажем, центру V(O2-)+(F+-центр) соответствует в решетке граната (TmBi)3(FeGa)5O12 значительное количество глубоких уровней захвата, объясняется большим количеством возможных вариантов набора «комплектов» ближайших соседей в первой и второй координационных сферах.
Заключение
Методами ТСТП и ТСТП КЗ впервые изучен энергетический спектр электрически активных центров в ЭМПФГ (TmBi)3(FeGa)5O12: Ca2+.
Установлено, что при увеличении концентрации ионов Ca2+ в пленках (TmBi)3(FeGa)5O12, энергетический спектр электрически активных центров сдвигается в зону более высоких температур.
Впервые экспериментально установлено, что спектр термостимулированных токов объектов исследования формируется 36-ю центрами, создающими в запрещенной зоне (TmBi)3(FeGa)5O12: Ca2+ электрически активные центры с энергетическим положением от 0,69 эВ до 1,21 эВ.
Предполагается, что обнаруженные локальные центры соответствуют однозарядным V(O2-)+(F+-центр) и нейтральным V(O2-)0(F-центр) кислородным вакансиям, находящимся в разном окружении в первой координационной сфере.
Работа выполнена в рамках проекта «Разработка нового класса наноразмерных материалов на основе пленочных магнитных электретов и мультиферроиков для сверхплотной магнитной и магнитооптической записи информации» (ГОСУДАРСТВЕННЫЙ КОНТРАКТ № 11.519.11.4026)
Литература
1. Рандошкин В.В., Червоненкис А.Я. Прикладная магнитооптика [Текст]: Монография / В.В. Рандошкин. – М.: Энергоатомиздат, 1990, – 320 с.
2. Костишин В.Г. Радиационно-стимулированные и короноэлектретные изменения структуры и свойств феррогранатовых гетерокомпоиций [Текст]: дис. докт. физ.-мат. наук: 01.04.10: – М.: МИСиС, 2009. – 264 с.
3. Гороховатский Ю.А., Бордовский Г.А. Термоактивационная токовая спектроскопия высокоомных полупроводников и диэлектриков [Текст]: Монография / Ю.А. Гороховатский. – М.: Наука: Главная редакция физико–математической литературы, 1991. – 248 с.
4. Савицкий Г.В., Костишин В.Г., Читанов Д.Н. и др. Установка термоактивационной токовой спектроскопии для исследования электрически активных дефектов в феррогранатовых гетерокомпозициях [Текст] // Тез. докл. IV–й Междунар. конф. По физике кристаллов «Кристаллофизика ХХI века», посвященной памяти М.П. Шаскольской. М.: НИТУ «МИСиС», 2010. – т. 2. – С. 361–362.
5. Находин Н.Г., Вознюк В.В. Термодеполяризация номинально чистого железоиттриевого граната [Текст] // УФЖ, 1986. – т. 31, – № 11, – С. 1750-1752.
6. Костюк В.Х., Костишин В.Г., Летюк Л.М., Подтяшкин Е.Я., Шипко М.Н. Исследование активных центров в феррит-гранатовых пленках с цилиндрическими магнитными доменами методом термостимулированных токов [Текст] // УФЖ, 1988. – т. 33, – № 2, – С. 1750-1752.
7. Шевчук В.Н., Костишин В.Г., Бугакова О.Е. Явления электропереноса в кристаллах Gd2,6Ca0,4Mg0,25Zr0,65Ga4,1O12 [Текст] // Неорганические материалы, 2000. – т. 36, – № 11. – С. 1369-1373.