Распределение ионов кобальта Cо2+ в монокристаллах шпинели Li0.5Ga2.5O4
- Authors: Шаповалов В.В.1, Шаповалов В.А.2, Дрокина Т.В.3, Воротынов А.М.3, Вальков В.И.2
-
Affiliations:
- Организация “Математика для Америки”
- ФГБНУ Донецкий физико-технический институт им. А.А. Галкина
- Институт физики им. Л.В. Киренского РАН
- Issue: Vol 125, No 1 (2024)
- Pages: 32-38
- Section: ЭЛЕКТРИЧЕСКИЕ И МАГНИТНЫЕ СВОЙСТВА
- URL: https://journals.rcsi.science/0015-3230/article/view/259710
- DOI: https://doi.org/10.31857/S0015323024010052
- EDN: https://elibrary.ru/ZREKJH
- ID: 259710
Cite item
Full Text
Abstract
Показано распределение ионов Cо2+ по подрешеткам и структурно неэквивалентным положениям в элементарной ячейке кристаллической решетки монокристалла литийгаллиевой шпинели Li0.5Ga2.5O4. Такое распределение определяет свойства как моно-, так и нанокристаллических веществ. Распределение обеспечивается специальной технологией и проявляется в спектрах электронного парамагнитного резонанса (ЭПР). Распределение ионов Cо2+ зависит от структурной и магнитной неэквивалентности. Структурная и магнитная неэквивалентность формирует многоминимумность потенциала кристаллического поля в элементарных ячейках монокристаллов в местах расположения ионов Cо2+. Ионы Cо2+ находятся в комплексах с тетраэдрическим и октаэдрическим окружением ионами кислорода. Обнаружены и исследованы три типа спектров ЭПР ионов Cо2+. Спектр Co2+тетр обязан иону Co2+, замещающему ион Ga3+, находящийся в тетраэдрическом кислородном окружении. Спектр Co2+окт, находящегося в кристаллическом поле аксиальной симметрии, принадлежит иону Co2+, замещающему ион Li+, находящийся в октаэдрическом кислородном окружении. Спектр Co2+окт, находящегося в кристаллическом поле низкой симметрии, принадлежат иону Co2+, замещающему ион Ga3+, находящийся в октаэдрическом кислородном окружении. Ближайшее катионное окружение иона создает ромбические искажения за счет разной валентности Li+ и Ga3+. Исследование угловых зависимостей спектров показало наличие в элементарных ячейках по четыре и по 12 магнитно-неэквивалентных положений.
Full Text
ВВЕДЕНИЕ
Активное изучение шпинелей в настоящее время обусловлено их широкими научными и технологическими применениями [1–4]. В материаловедении, как правило, используются шпинели с определенными примесями. Такие примеси располагаются в октаэдрических и тетраэдрических комплексах элементарных ячеек кристаллических решеток шпинелей. В комплексах магнитные ионы Co²⁺находятся в кислородном окружении. Распределение ионов кобальта Co²⁺ зависит от структурной и магнитной неэквивалентности ионов и от многоминимумности потенциала кристаллического поля в элементарных ячейках монокристаллов. От распределения ионов зависят свойства материалов, применяемых в приборах и лечебных препаратах. Однако технология анализа отдельного комплекса в материале до сих пор не разработана. Поэтому актуальны исследования магнитного иона в элементарной ячейке монокристалла, что возможно методом электронного парамагнитного резонанса в монокристалле путем исследований угловых и температурных зависимостей спектров ЭПР магнитных ионов. В работах [5–13] подробно изучены нанокристаллы Li0,5Ga2,5O4, допированные кобальтом.
В настоящее время наблюдается рост научного интереса к соединениям со структурой шпинели благодаря возможностям синтеза нанокристаллических образцов [7]. Но при исследовании нанокристаллических образцов появляются затруднения, которых нет при исследовании монокристаллов. Наиболее характерное свойство монокристалла – анизотропия, зависимость большинства физических свойств от направления. В нанокристаллических материалах нет такой зависимости, так как в реальном соединении анизотропия свойств усредняется. Интерпретация свойств нанокристаллических и монокристаллических материалов, имеющих одинаковый состав, является актуальной задачей [8–17]. Однако в этих работах не учитываются свойства отдельного комплекса с магнитным 3dn- ионом в используемых композитных материалах (монокристаллы, нанокристаллы, полимеры и др.). В нанокомплексе исследуемый металлический ион находится в лигандном окружении. Такой нанокомплекс представляет собой часть элементарной ячейки материала. В каждом нанокомплексе окружение иона имеет кристаллографическую симметрию, такую же, как и в монокристалле. Магнитные оси 3dn-ионов располагаются, как правило, вдоль кристаллографических осей комплекса. В работах [8–17] не учитывается существование структурной неэквивалентности вводимых магнитных ионов в элементарную ячейку кристаллической решетки. Поэтому объективная информация об исследуемых соединениях затруднительна. До сих пор имеется ряд проблем в получении достоверных характеристик – определении местоположения иона, его ближайшего окружения, параметров электрического поля в месте расположения иона и др. Анализ структуры комплекса с металлическим ионом в композитных материалах [8–18] делали без исследований параметров электрического кристаллического поля в месте нахождения металлического иона, что возможно методом ЭПР.
Для исследования комплексов в монокристалле Li0,5Ga2,5O4 с ионом кобальта мы применили один из наиболее информативных методов исследования электронной структуры соединений с точечными примесями – метод ЭПР [19]. В данной работе проведены экспериментальные исследования по обнаружению и изучению спектров ЭПР ионов Co²⁺ в монокристаллах литийгаллиевой шпинели Li0.5Ga2.5O4 в интервале температур от гелиевых до комнатных.
Рис. 1. Элементарная ячейка монокристалла литийгаллиевой шпинели Li0,5Ga2,5O4. Ближайшее окружение ионов кобальта Cо²⁺ (позиции А, В) – кислородные тетраэдры и октаэдры.
Элементарная ячейка монокристалла литийгаллиевой шпинели Li0,5Ga2,5O4 изображена на рис. 1. Элементарная ячейка содержит кислородные октаэдры и тетраэдры, внутри которых расположены магнитные зонды – ионы кобальта Cо²⁺, занимающие А (тетраэдрические позиции) и В (октаэдрические позиции).
МЕТОДИКА ЭкспериментА
Методика эксперимента и техника измерений определяется физической задачей, для успешного решения которой необходимо выполнить ряд требований, предъявленных к экспериментальной установке. Вышеизложенные факты дали основание для выбора рабочей частоты радиоспектрометра v = 37.5 ГГц (λ = 8 мм). Магнитные поля для выполнения условий g = 2 должны быть более 13.5 кгс.
Для выполнения поставленной задачи выбран супергетеродинный тип спектрометра с отражательным резонатором, на котором возможно получение предельно высокой чувствительности и который хорошо приспособлен для работы при низких температурах. Одним из основных элементов радиоспектрометра является низкотемпературный прямоугольный резонатор отражательного типа. Подстройка частоты резонатора по диапазону клистрона осуществляется бесконтактным поршнем. Для ориентации магнитных осей в резонаторе предусмотрено устройство для вращения образцов с выведенным лимбом для отсчета углов. Кристалл наклеивается на поворотный столик, укрепленный на узкой боковой стенке резонатора в пучности микроволнового магнитного поля Н1, имеющего вертикальное направление. Так как постоянное поле при вращении магнита имеет горизонтальное направление, то всегда выполняется условие H1 ^ H0.
Для решения поставленной физической задачи необходимо было в широких пределах менять рабочую температуру исследуемого кристалла. Поэтому резонатор помещали в четырехстенный гелиевый сосуд Дьюара, который заканчивается трехстенным отростком – пальцем, помещенным между наконечниками магнита. Жидкий гелий обладает диэлектрической проницаемостью и поэтому при попадании в резонатор не изменяет существенно частоту. Однако за счет теплоподвода жидкий гелий в резонаторе кипит, что приводит к значительным помехам при измерении спектров ЭПР. Поэтому резонатор изолировали от гелиевой ванны тонкостенным сосудом из нержавеющей стали. Для изменения температуры от 4.2° вплоть до азотных Т = 77 К применяли нагревательную печь из константановой проволоки R = 200 Ом, которую наматывали на тело резонатора. Питание печи производится от источника постоянного тока. Для уменьшения теплоотвода от резонатора применяется волноводный мост из нержавеющей стали, а резонатор вакуумируется. Аналогично повышали температуру от Т = 77 К до комнатной при заливке в сосуд Дьюара жидкого азота. Датчиками температуры в диапазонах температур 4.2–77 К и 55–300 К служили помещенные на внешней стенке резонатора в непосредственной близости от образца угольное и медное сопротивление. Градуировка датчиков температуры произведена по методике, описанной в литературе [20].
Результаты
Электронная конфигурация двухвалентного кобальта 3d7, L = 3 и S = 3/2. Основным термом Co²⁺, находящегося в свободном состоянии, согласно правилу Хунда является 4F9/2. Следовательно, основное спектроскопическое состояние 4F9/2 семикратно вырождено по орбитальному квантовому числу (2L+1) и четырехкратно вырождено по спину (2S+1). В чисто кубическом кристаллическом поле октаэдрической симметрии вырождение частично снимается, нижним уровнем является орбитальный триплет, каждый уровень которого четырехкратно вырожден по спину. При воздействии кристаллических полей с симметрией ниже кубической и спин-орбитальной связи триплет Г4 расщепляется на шесть крамерсовых дублетов, расстояние между которыми порядка сотен обратных сантиметров. Поэтому из-за сильной спин-орбитальной связи спектр ЭПР Co²⁺ в октаэдрическом окружении можно наблюдать только при низких температурах на нижнем дублете с эффективным спином S ’ = 1/2 (истинный спин S = 3/2). Спектр ЭПР обладает сильной анизотропией и существенным отклонением g-фактора от g -фактора для чисто спинового состояния.
Монокристаллы со структурой шпинели Li0.5Ga2.5O4 выращиваются методом спонтанной кристаллизации из раствора окислов в расплавленных солях. Технология получения монокристаллов Li0.5Ga2.5O4 аналогична технологии роста соответствующих ферритов [21]. При выращивании использовали следующий температурный режим: нагрев до 1150 °C и выдержка при этой температуре 10 ч. Охлаждение от 1150 °C до 1000 °C по программе 6 град/ч. При 1000 °C выдержка 30 ч, охлаждение до 600 °C по программе 2 град/ч, от 600 °C до 400 °C – 6 град/ч. При 400 °C печь выключали. Такой температурный режим обеспечивает получение кристаллов в упорядоченном состоянии. Выращенные монокристаллы литийгаллиевой шпинели имели форму октаэдров с размером грани до 4 мм. Содержание вводимых парамагнитных ионов меняли в пределах 0.01–0.5 вес.%. Распределение Cо²⁺ происходит в процессе роста монокристалла и обеспечивается заданием параметров роста монокристалла.
Спектр ЭПР ионов Co²⁺ исследовали на частоте 37 ГГц в интервале температур 4.2–77 К. Исследование спектра ЭПР Co²⁺ Li0.5Ga2.5O4 при Т=4.2 К в плоскостях {110} и {111} показало наличие трех типов спектров ЭПР в кристаллическом поле аксиальной и ромбической симметрии и Co²⁺тетр в поле аксиальной симметрии.
Измерения спектра ЭПР Co²⁺ромб проведены при Т=4.2 К. При Т=77 К спектр не наблюдается. Спектр ЭПР состоит из разрешенной сверхтонкой структуры (СТС) вдоль осей z и x и не разрешен для H0 // y. Исследование угловой зависимости спектра показало наличие двенадцати магнитно-неэквивалентных положений Co²⁺ромб в элементарной ячейке. Такая низкая симметрия кристаллического поля приводит к существованию трех значений g-факторов в трех главных ориентациях, а также трех значений констант Aij – тензора сверхтонкого взаимодействия. Каждое положение Co²⁺ описывается спин-гамильтонианом:
. (1)
с эффективным спином S’ = 1/2 и ядерным спином I = 7/2.
Ось z кристаллического поля отклонена на угол 14° от оси <111>, ось x совпадает с направлением <110>. Экспериментальные значения констант спин-гамильтониана в трех главных ориентациях:
(2)
Линия ЭПР расщепляется на 8 компонент СТС.
Рис. 2. Угловая зависимость положений линий спектра ЭПР для ионов Co²⁺ аксиальной симметрии в плоскости {110} при Т = 4.2К [22]. Буквами о и m рядом с цифрой на линии обозначены октаэдрические и тетраэдрические позиции Co²⁺ соответственно.
Измерения “аксиального” спектра Co²⁺окт в кислородном октаэдре проведены при Т = 4.2 К. При 77 К спектр ЭПР не наблюдается. Анализ угловой зависимости спектра (рис. 2) говорит о наличии четырех магнитно-неэквивалентных положений ионов Co²⁺окт в элементарной ячейке. Спектр ЭПР в параллельной ориентации состоит из хорошо разрешенной сверхтонкой структуры из восьми линий, которая обусловлена ядерным спином I = 7/2. В перпендикулярной ориентации СТС не разрешена. Спектр разрешается при отклонении магнитного поля от перпендикулярной ориентации в плоскости {110} на угол ±6°. Параллельная ориентация кристаллического поля совпадает с осью <111>. Для проверки аксиальности спектра снимали угловую зависимость в плоскости {111}, которая является гранью октаэдра. На рис. 2 показана угловая зависимость для ионов Co²⁺, находящихся в электрическом кристаллическом поле аксиальной симметрии в плоскости {110} при температуре Т = 4.2 К. Линии, обозначенные цифрами без штрихов, соответствуют первой линии СТС, со штрихами – восьмой. Кривые 1, 2 показывают угловую зависимость для ионов, параллельные ориентации которых совпадают с направлениями <111>, лежащими в исследуемой плоскости {110} (рис. 3). Кривые 3, 4 и 3’, 4’ представляют собой спектр ЭПР от двух магнитно неэквивалентных ионов, параллельные ориентации которых находятся в плоскости, которая перпендикулярна исследуемой плоскости {110}.
Спектр каждого иона Co²⁺ описывается спин-гамильтонианом (3):
(3)
с эффективным спином S’ = 1/2 и ядерным спином I = 7/2.
Получены следующие экспериментальные значения констант спин-гамильтониана (4):
(4)
Изучение угловых зависимостей положения линий спектра ЭПР при 4.2 К в плоскостях {110} и {111} показало наличие Co²⁺тетр в поле аксиальной симметрии.
Тетраэдрический Co²⁺. Основное спектроскопическое состояние 4F9/2. В кубическом поле тетраэдрической симметрии нижним уровнем является орбитальный синглет, который отстоит от ближайшего орбитального триплета на расстоянии порядка 104 см‑1 и четырехкратно вырожден по спину. Нижний орбитальный синглет расщепляется на два спиновых дублета, нижним из которых является уровень с Sz = ±½, спектр четко виден как при гелиевых, так и при азотных температурах, причем интенсивность спектра при гелиевых температурах гораздо выше.
Рис. 3. Плоскость {110} правильного октаэдра.
Спектр ЭПР при Т = 4.2 К исследован в плоскости {110} (рис. 3). Параллельная ориентация совпадает с направлением <111>, сверхтонкая структура ЭПР в этой ориентации разрешена. При отклонении от параллельной ориентации на угол ± 4 °CТС, состоящая из восьми компонент согласно I = 7/2, исчезает. При Т = 77 К СТС не разрешается. Изучение угловой зависимости положения линий спектра ЭПР показало наличие четырех магнитно неэквивалентных положений ионов в элементарной ячейке. Спектр каждого иона описывается гамильтонианом аксиальной симметрии. Эффективный спин S = 1/2, так как мы наблюдали переход только на нижнем уровне Sz = ± ½. Константы гамильтониана при Т = 4.2 К:
(5)
Обсуждение результатов
Обсудим экспериментальные результаты по исследованию ионов Co²⁺ в монокристалле Li0.5Ga2.5O4. Было обнаружено и изучено два октаэдрических центра Co²⁺, которым соответствуют спектры аксиальной и ромбической симметрии. Соотношение интенсивностей спектров октаэдрических Co²⁺акс: Co²⁺ромб = 2:1 при Т = 4.2 К. Разная интенсивность спектров позволила четко их разделить.
Третий тип спектров обязан иону Co²⁺, находящемуся в тетраэдрическом окружении. Основанием для этого предположения послужили следующие факты:
Обнаруженные и исследованные три типа спектров ЭПР позволяют сделать следующие заключения:
Объяснение распределения Co²⁺ в монокристаллах Li0.5Ga2.5O4 представлены на рис. 4. Показана зависимость потенциала кристаллического поля E от расстояния R в относительных единицах. Расстояние R отложено в относительных единицах, т. к. концентрация ионов Co²⁺мала и расстояния между ионами в элементарных ячейках различны.
Изучение угловых зависимостей спектров ЭПР позволило определить форму потенциала кристаллического поля в месте расположения ионов Co²⁺.
Рис. 4. Зависимость потенциала кристаллического поля E от расстояния R. Минимумы расположены вдоль осей типа [111]. Показаны тетраэдрические (t) и октаэдрические (o) узлы с ионами Co²⁺ в элементарной ячейке.
Вдоль кристаллографических направлений типа [111] форма многоминимумного потенциала имеет эквивалентные минимумы, разделенные потенциальными барьерами. На рис. 4 стрелками показаны направления <111>. Эти направления перпендикулярны плоскостям (111) октаэдров и тетраэдров. Приведены четыре минимума. Магнитное поле Н0 параллельно главной магнитной оси z ионa Co²⁺, направление показано стрелкой. В элементарной ячейке имеются магнитно неэквивалентные положения ионов Co²⁺, расположенные в тетраэдрическом и октаэдрическом кислородном окружении. Показаны такие комплексы в виде октаэдров и тетраэдров. Положения ионов Co²⁺ обозначены буквами о и t. Показаны тетраэдрические – Co(1t), Co(2t), Co(3t), Co(4t) и октаэдрические – Co(1o), Co(2o), Co(3o), Co(4o) неэквивалентные положения ионов кобальта.
Выводы
В элементарной ячейке монокристалла Li0,5Ga2,5O4 допированные ионы кобальта Co²⁺имеют различные предпочтительные возможности для распределения в четырех тетраузлах, замещая ионы Ga³⁺, четырех октаузлах аксиальной симметрии, замещая ионы Li+, и в двенадцати октаузлах низкой симметрии, замещая ионы Ga³⁺. Такое распределение ионов кобальта определяет свойства веществ, как монокристаллических, так и нанокристаллических. Таким образом, существует возможность получения как минимум двадцати типов материалов литийгаллиевой шпинели Li0.5Ga2.5O4, допированной ионами Co²⁺, с заданными свойствами. Достижение такого результата может быть обеспечено с помощью определенной технологии изготовления как моно-, так и нанокристаллических материалов. В случае смешанного распределения ионов Cо²⁺ возможно получение гораздо большего числа новых типов материалов.
Данная работа финансировалась за счет средств бюджета институтов и организации. Никаких дополнительных грантов на проведение или руководство данным конкретным исследованием получено не было.
Авторы данной работы заявляют, что у них нет конфликта интересов.
About the authors
В. В. Шаповалов
Организация “Математика для Америки”
Email: vashapovalov1@mail.ru
United States, Нью-Йорк
В. А. Шаповалов
ФГБНУ Донецкий физико-технический институт им. А.А. Галкина
Author for correspondence.
Email: vashapovalov1@mail.ru
Russian Federation, ул. Розы Люксембург, 72, Донецк, 283048
Т. В. Дрокина
Институт физики им. Л.В. Киренского РАН
Email: vashapovalov1@mail.ru
Russian Federation, Академгородок, 50/12, Красноярск, 660036
А. М. Воротынов
Институт физики им. Л.В. Киренского РАН
Email: vashapovalov1@mail.ru
Russian Federation, Академгородок, 50/12, Красноярск, 660036
В. И. Вальков
ФГБНУ Донецкий физико-технический институт им. А.А. Галкина
Email: vashapovalov1@mail.ru
Russian Federation, ул. Розы Люксембург, 72, Донецк, 283048
References
- Tsurkan V., Nidda H., Deisenhofer J., Lunkenheimer P., Loidl A. On the complexity of spinels: Magnetic, electronic, and polar ground states // Phys. Reports. 2021. V. 926. P. 1–86.
- Maigny L., Dupont M. Spinels: Occurrences, Physical Properties and Applications / Nova Sci. Publishers. Inc: New York. USA. 2013.
- Ganesh I.A. Review on Magnesium Aluminate (MgAl2O4) Spinel: Synthesis, Processing and Applications // Int. Mater. Rev. 2013. V. 58. P. 63–112.
- Zou Y., Gréaux S., Irifune T., Li B., Higo Y. Unusual Pressure Effect on the Shear Modulus in MgAl2O4 Spinel // J. Phys. Chem. C2013. V. 117. P. 24518–24526.
- Riesen H., Yildirim B. Persistent spectral hole-burning in diffuse reflection: application to nanocrystalline LiGa5O8: Co²⁺ // J. Phys. Chem. Lett. 2010. V. 1. P. 2380–2384.
- Yildirim B. Spectroscopic and Structural Investigations of Transition Metal Ion Doped Nanocrystalline LiGa5O8. http://handle.unsw.edu.au/1959.4/ 52631
- Sickafus K.E., Wills J.M. Spinel Coumpounds: Structure and Property Relations // J. Am. Ceram. Soc. 1999. V. 82. P. 3279–3292.
- Macfarlane R.M., Vial J.C. Photon-gated spectral hole burning in LiGa₅O₈: Co²⁺ // Phys. Rev. B. 1986. V. 34. P. 1–4.
- Malyarevich A.M., Yumashev K.V. Saturable absorbers based on tetrahedrally coordinated transition-metal ions in crystals // J. Appl. Spectr. 2009. V. 76. P. 1–43.
- Donegan J.F., Anderson F.G., Bergin F.J., Glynn T.J., Imbusch G.F. Optical and magnetic-circular-dichroism optically-detected-magnetic-resonance study of the Co²⁺ ion in LiGa₅O₈ // Phys. Rev. B. 1992. V. 45. P. 563–573.
- Ishihara T., Tanaka K., Hirao K., Soga N. Synthesis and optical properties of transparent LiGa₅O₈: Co²⁺ films by the sol-gel method // Bull. Chem. Soc. Jap. 1997. V. 70. P. 1347–1351.
- Lever A.B.P. Inorganic Electronic Spectroscopy. Elsevier Science Publisher: Amsterdam, 1984. P. 863.
- Riesen H., Badek K., Stevens-Kalceff M. Correlation between inhomogeneous width and crystallite size: ²E→ ⁴A₂ luminescence of Co(II) in nanocrystalline ZnAl₂O₄ // Chem. Phys. Lett. 2010. V. 501. P. 103–107.
- Burda C., Chen X.B., Narayanan R., El-Sayed M.A. Chemistry and properties of nanocrystals of different shapes // Chem. Rev. 2005. V. 105. P. 1025–1102.
- Talapin D.V., Lee J.S., Kovalenko M.V., Shevchenko E.V. Prospects of colloidal nanocrystals for electronic and optoelectronic applications // Chem. Rev. 2010. V. 110, P. 389–458.
- Träger F. Editorial: Optical properties of nanoparticles // Appl. Phys. B. 2001. V. 73. P. 291.
- Jaramillo T.F., Baeck S.H., Kleiman- Shwarsctein A., Choi K.S., Stucky G.D., McFarland E.W. Automated electrochemical synthesis and photoelectrochemical characterization of Zn₁₋ₓCoₓO thin films for solar hydrogen production // J. Combinat. Chem. 2005. V. 7. P. 264–271.
- Huppertz H., Hering S.A., Zvoriste C.E., Lauterbach S., Oeckler O., Riedel R., Kinski I. High-Pressure Synthesis, Electron Energy-Loss Spectroscopy Investigations, and Single Crystal Structure Determination of a Spinel-Type Gallium Oxonitride Ga₂.₇₉□₀.₂₁(O₃.₀₅N₀.₇₆□₀.₁₉)// Chem. Mater. 2009. V. 21. P. 2101–2107.
- Abragam A., Bleaney B. Electron Paramagnetic Resonance of Transition Ions. Clarendon press: Oxford, U.K. 1970. P. 928.
- Уайт Т.К. Экспериментальная техника в физике низких температур. М.: Физматгиз, 1961.
- Селезнев В.Н., Пухов И.К., Дрокин А.К, Шаповалов В.А. Магнитная кристаллографическая анизотропия монокристаллов литиевого и литий-цинкового ферритов с малыми добавками кобальта // ФТТ. 1970. Т. 12. Вып. 3. С. 885–891.
- Кожухарь А.Ю., Селезнев В.Н., Цинцадзе Г.А., Шаповалов В.А. Электронный парамагнитный резонанс ионов Co²⁺ в литий-галлиевой шпинели // ФТТ. 1973. Т. 15. Вып. 3. С. 706–709.
Supplementary files
