Радиотехника и электроника, 2019, T. 64, № 8, стр. 819-825

Особенности низкотемпературной теплоемкости монокристаллов гранатов Er3 – xTmxAl5O12

Е. В. Чарная 1*, Е. В. Шевченко 1, Е. Н. Хазанов 2**, А. В. Таранов 2, А. М. Уляшев 3

1 Санкт-Петербургский государственный университет
198504 С.-Петербург, Петергоф, ул. Ульяновская, 3, Российская Федерация

2 Институт радиотехники и электроники им. В.А. Котельникова РАН
125009 Москва, ул. Моховая, 11, стр. 7, Российская Федерация

3 Сыктывкарский государственный университет им. Питирима Сорокина
167001 Сыктывкар, Октябрьский просп., 55, Российская Федерация

* E-mail: charnaya@mail.ru
** E-mail: khazanov@cplire.ru

Поступила в редакцию 16.05.2018
После доработки 16.05.2018
Принята к публикации 14.06.2018

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

Аннотация

Представлены результаты исследований теплоемкости монокристаллов гранатов Er3– xTmxAl5O12 (x = 0, 1, 2, 3) в магнитных полях до 9 T и в диапазоне температур 1.9…220 K. Температурные зависимости теплоемкости аппроксимированы суммой вкладов аномалий Шоттки, связанных с парамагнитными ионами Er3+ и Tm3+, и решеточных вкладов Дебая и Эйнштейна. По данным для теплоемкости рассчитаны энтропия и магнитная энтропия. Показано, что с ростом магнитного поля энтропия уменьшается, что указывает на возможность применения исследуемых гранатов в методе адиабатического размагничивания.

ВВЕДЕНИЕ

Алюмо-редкоземельные гранаты Re3Al5O12 (Re – различные редкоземельные элементы) являются широко известными лазерными и магнитооптическими кристаллами [1, 2]. В последнее время, однако, были выявлены новые свойства алюмо-редкоземельных гранатов, перспективные в плане других возможных применений. В работе [3] был обнаружен метамагнетизм в смешанном алюмо-эрбий-гольмиевом гранате, аналогичный метамагнитным явлениям, ранее наблюдавшимся в монокристалле галлий-тербиевого граната [4]. В работах [57] на основании изучения магнитокалорических эффектов в гранатах Re3Al5O12 было предложено использовать эти кристаллы для охлаждения методом адиабатического размагничивания. Монокристаллы гранатов с парамагнитными ионами обладают структурной стабильностью и магнитной обратимостью, что ставит их в один ряд с различными соединениями и металлоорганическими материалами, которые рассматриваются в настоящее время для магнитного охлаждения [8, 9]. Тем не менее, в литературе имеется ограниченная информация о термодинамических свойствах алюмо-редкоземельных гранатов и о влиянии магнитного поля на их теплоемкость. Теплоемкость кристалла Y3Al5O12 с 30% ионов Y, замещенных на ионы Er, была измерена в работе [5] в температурном интервале от 93 мК до 8 К в нулевом магнитном поле и в полях до 8 Tл. Максимум теплоемкости соответствовал температуре ≈266 мК в нулевом магнитном поле. Теплоемкость чистого алюмо-эрбиевого граната и твердых растворов Er3 –xYxAl5O12 в области температур от 1.9 до 220 K во внешних магнитных полях до 9 Tл была измерена в работе [6]. Полученная температурная зависимость теплоемкости интерпретирована с учетом вкладов решеточных колебаний и аномалий Шоттки. Кроме того, было рассмотрено расщепление нижнего крамерсовского дублета магнитными полями, создаваемыми соседними ионами. Насколько нам известно, теплоемкость чистого алюмо-тулиевого граната и смешанных алюмо-эрбий-тулиевых гранатов ранее не исследовали, несмотря на то что тулиевые гранаты широко применяются на практике [10]. Следует также отметить, что данные по низкотемпературной теплоемкости являются существенными при анализе кинетики тепловых фононов [11].

В данной работе приведены результаты измерения теплоемкости монокристаллов Er3 –xTmxAl5O12 (x = 0, 1, 2, 3) в диапазоне температур от 1.9 до 220 K в магнитных полях до 9 Tл.

1. ХАРАКТЕРИСТИКИ ОБРАЗЦОВ И МЕТОД ИССЛЕДОВАНИЯ

Чистые и смешанные монокристаллы гранатов Er3– xTmxAl5O12 (x = 0, 1, 2, 3) были выращены методом направленной горизонтальной кристаллизации в молибденовом контейнере. Образцы вырезали из внутренних частей буль. Рентгеновские исследования показали, что образцы были однофазными и свободными от двойникования. Рассчитанные и измеренные плотности гранатов представлены в табл. 1. Незначительное различие между этими плотностями свидетельствует о высоком качестве выращенных кристаллов. Плотности однокомпонентных монокристаллов Er3Al5O12 и Tm3Al5O12 согласуются с данными, опубликованными в [12]. Образцы для измерения теплоемкости представляли собой пластинки, вырезанные перпендикулярно кубическим осям, толщиной около 0.3 мм и площадью поверхности около 0.2 см2.

Таблица 1.  

Рассчитанные и измеренные плотности гранатов Er3– xTmxAl5O12

x ρрас, г/см3 ρизм, г/см3
0 6.396 6.410
1 6.424 6.414
2 6.451 6.393
3 6.479 6.482

Алюмо-тулиевые и алюмо-эрбиевые гранаты образуют твердые растворы при всех концентрациях благодаря близости радиусов трехвалентных ионов эрбия Er3+ и тулия Tm3+ [2]. Редкоземельные ионы Er3+ и Tm3+ занимают в кристаллической решетке гранатов додекаэдрические c-позиции. Элементарная ячейка содержит 160 ионов (восемь формульных единиц). Гранаты имеют пространственную симметрию $Ia\bar {3}d.$

Теплоемкость измеряли с помощью комплекса для измерения физических свойств PPMS-9 + EverCool-II фирмы Quantum Design. Измерения проводили в температурном интервале 1.9…220 K в постоянных магнитных полях от 0 до 9 Tл.

2. ОБСУЖДЕНИЕ РЕЗУЛЬТАТОВ ИЗМЕРЕНИЙ

Данные для теплоемкости монокристаллов гранатов Er3– xTmxAl5O12 в нулевом магнитном поле представлены на рис. 1. На вставках к рисунку для наглядности показана теплоемкость в узких температурных диапазонах. Для чистого (х = 0) кристалла Er3Al5O12 результаты были опубликованы ранее [6] и приведены на рис. 1 для полноты. Теплоемкость гранатов Er3– xTmxAl5O12 при низкой температуре определяется преимущественно аномалиями Шоттки [5, 13], обусловленными низколежащими уровнями парамагнитных ионов Er и Tm. При более высоких температурах преобладает вклад решеточных колебаний [13, 14].

Рис. 1.

Температурные зависимости теплоемкости гранатов Er3– xTmxAl5O12 в нулевом магнитном поле. Стрелки показывают возрастание x. На вставках показана теплоемкость в узких интервалах температуры.

Приложение магнитного поля существенно влияет на низкотемпературную теплоемкость за счет изменений аномалий Шоттки. Температурные зависимости теплоемкостей кристаллов Er2TmAl5O12, ErTm2Al5O12 и Tm3Al5O12 ниже 40 K в полях до 9 Tл показаны на рис. 2а–2в соответственно. Видно, что максимумы теплоемкости, соответствующие вкладам Шоттки от ионов Er3+ и Tm3+, с ростом магнитного поля смещаются к высоким температурам и становятся более пологими.

Рис. 2.

Температурные зависимости теплоемкости гранатов Er2TmAl5O12 (а), ErTm2Al5O12 (б) и Tm3Al5O12 (в) во внешнем магнитном поле до 9 Тл в области температур ниже 40 K. Стрелки показывают возрастание поля. На вставках представлено поведение теплоемкостей в узком температурном интервале.

Теплоемкость диэлектрических кристаллов с парамагнитными ионами может быть обусловлена несколькими независимыми вкладами: решеточным (фононным) вкладом, аномалиями Шоттки и аномалиями, вызванными магнитными и структурными фазовыми переходами [13]. Известно, что алюмогранаты не претерпевают структурных фазовых переходов вплоть до плавления [2]. Кроме того, кристаллическая симметрия гранатов сильно влияет на магнитное упорядочение. Магнитные моменты редкоземельных ионов в додекаэдрических позициях образуют гиперкагоме структуру, состоящую из двух взаимопроникающих треугольных подрешеток [3]. В такой структуре магнитные моменты не могут быть попарно антиферромагнитно ориентированы. Это приводит к магнитной фрустрации и нетривиальному отклику на приложенное магнитное поле, что наблюдалось, например, для Tb3Ga5O12 [3] и Er2HoAl5O12 [15, 16]. Гиперкагоме структура гранатов приводит также к понижению температур магнитного упорядочения по сравнению с оценками, полученными из закона Кюри–Вейсса. Антиферромагнитный фазовый переход в гранате Er3Al5O12 наблюдается при 0.155 K [17] и понижается с ростом x в твердых растворах Er3– xYxAl5O12 [3, 17]. Согласно теоретическим оценкам в чистом гранате Tm2Al5O12 электронное магнитное упорядочение не происходит из-за слабого взаимодействия между ионами Tm [18]. Таким образом, температура Нееля в смешанных эрбий-тулиевых гранатах не может превышать 0.155 K. Действительно, мы не наблюдали аномалии теплоемкости λ-типа, связанные с фазовыми переходами, вплоть до 1.9 K (см. рис. 1). Следовательно, для количественной интерпретации теплоемкости в исследуемых гранатах необходимо рассмотреть только вклады колебаний решетки и аномалий Шоттки.

Аномалии Шоттки в теплоемкости возникают за счет заселенности штарковских подуровней парамагнитных ионов в кристаллическом поле. Молярный вклад Шоттки при температуре T может быть записан в виде [13]

(1)
$\begin{gathered} {{C}_{S}} = \frac{{nR}}{{{{T}^{2}}}}\left[ {\frac{{\sum\limits_{i = 1} {\Delta _{i}^{2}{{g}_{i}}\exp \left( {{{ - {{\Delta }_{i}}} \mathord{\left/ {\vphantom {{ - {{\Delta }_{i}}} T}} \right. \kern-0em} T}} \right)} }}{{{{g}_{0}} + \sum\limits_{i = 1} {{{g}_{i}}\exp \left( {{{ - {{\Delta }_{i}}} \mathord{\left/ {\vphantom {{ - {{\Delta }_{i}}} T}} \right. \kern-0em} T}} \right)} }} - } \right. \\ - \,\,\left. {{{{\left( {\frac{{\sum\limits_{i = 1} {{{\Delta }_{i}}{{g}_{i}}\exp \left( {{{ - {{\Delta }_{i}}} \mathord{\left/ {\vphantom {{ - {{\Delta }_{i}}} T}} \right. \kern-0em} T}} \right)} }}{{{{g}_{0}} + \sum\limits_{i = 1} {{{g}_{i}}\exp \left( {{{ - {{\Delta }_{i}}} \mathord{\left/ {\vphantom {{ - {{\Delta }_{i}}} T}} \right. \kern-0em} T}} \right)} }}} \right)}}^{2}}} \right], \\ \end{gathered} $
где n – число парамагнитных ионов в базисе кристаллической решетки, R – газовая постоянная, ${{\Delta }_{i}}$ – разность энергий в кельвинах между i-м возбужденным и основным уровнями, g0 и gi – кратности вырождения основного и i-го возбужденного уровней.

Система уровней энергии трехвалентного иона эрбия в алюмо-иттриевом гранате приведена в работе [19]. Основной мультиплет 4I15/2 расщепляется во внутрикристаллическом электрическом поле на восемь крамерсовских дублетов. Первый возбужденный дублет отстоит от основного на 32 K (22 см–1). Энергии следующих двух дублетов равны 86 и 115 K. Остальные дублеты лежат выше 600 K.

Вырождение основного мультиплета 3H6 трехвалентного иона Tm3+ в решетке граната полностью снимается, что дает совокупность 13 синглетов [20, 21]. Первый возбужденный синглет Tm3+ в алюмо-иттриевом гранате отстоит от основного уровня на 39 K (27 см–1), следующие четыре синглета имеют энергии 311, 347, 355 и 363 K, остальные синглеты имеют энергию выше 800 K.

При аппроксимации вкладов Шоттки энергии самых низких возбужденных уровней для ионов эрбия и тулия слегка изменяли для лучшего количественного согласия с экспериментом. Энергии остальных состояний принимали равными значениям, найденным в [19, 20]. При расчетах учитывали только уровни, отстоящие от основного уровня не более чем на 600 K.

Наблюдаемый рост теплоемкости в эрбиевом и эрбий-тулиевом гранатах при понижении температуры ниже 5 K не может быть связан возбужденными уровнями, поскольку они расположены выше 30 К (см. выше). Аналогичное возрастание теплоемкости ниже 5 K наблюдали для смешанных алюмо-иттрий-эрбиевых кристаллов [6]. Оно было интерпретировано в работах [5, 6] как результат расщепления на величину ${{\Delta }_{{{\text{д о п }}}}}$ нижнего крамерсовского дублета при снятии вырождения в магнитных полях, создаваемых соседними магнитными ионами. Это же предположение было использовано и в данной работе.

Решеточный вклад в теплоемкость обычно описывают комбинацией вкладов Дебая и Эйнштейна [15]. Модель Дебая дает для молярной теплоемкости выражение

(2)
${{C}_{{\text{D}}}}(T) = 3{{r}_{{\text{D}}}}R{{\left( {\frac{T}{\theta }} \right)}^{3}}\int\limits_0^{{\theta \mathord{\left/ {\vphantom {\theta T}} \right. \kern-0em} T}} {\frac{{{{x}^{4}}\exp (x)}}{{{{{\left( {\exp (x) - 1} \right)}}^{2}}}}dx} ,$
где θ – температура Дебая и rD – число колебательных мод, рассматриваемых в рамках модели Дебая. Вклад Эйнштейна может быть записан в виде
(3)
${{C}_{{\text{E}}}}(T) = {{r}_{{\text{E}}}}R{{\left( {\frac{{{{\theta }_{{\text{E}}}}}}{T}} \right)}^{2}}\frac{{\exp \left( {{{{{\theta }_{{\text{E}}}}} \mathord{\left/ {\vphantom {{{{\theta }_{{\text{E}}}}} T}} \right. \kern-0em} T}} \right)}}{{{{{\left( {\exp \left( {{{{{\theta }_{{\text{E}}}}} \mathord{\left/ {\vphantom {{{{\theta }_{{\text{E}}}}} T}} \right. \kern-0em} T}} \right) - 1} \right)}}^{2}}}},$
где ${{\theta }_{{\text{E}}}}$ – температура Эйнштейна и rE – число колебательных мод, рассматриваемых в рамках теории Эйнштейна.

Количественную интерпретацию температурной зависимости теплоемкости в нулевом магнитном поле проводили с использованием комбинации вкладов (1)–(3). Множитель rD в выражении (2) рассматривали как подгоночный параметр. Температуру Дебая θ полагали одинаковой для всех гранатов, поскольку в работе [22] было показано, что для алюмогранатов θ слабо зависит от состава. Для немагнитного граната Y3Al5O12 температура Дебая θ = 800 K была определена в работе [23]. Эта же величина θ была принята при аппроксимации низкотемпературной теплоемкости для смешанных эрбий-иттриевых и эрбий-гольмиевых гранатов как в работе [6], так и в данной работе.

Аномалии Шоттки рассчитывали для ионов Er3+ и Tm3+. Численную интерпретацию теплоемкости проводили для температур ниже 45 K, для которых использование только двух вкладов Эйнштейна давало хороший результат. Две температуры Эйнштейна ${{\theta }_{{{\text{E}}1}}}$ и ${{\theta }_{{{\text{E}}2}}},$ а также коэффициенты ${{r}_{{{\text{E}}1}}}$ и ${{r}_{{{\text{E}}2}}}$ играли роли подгоночных параметров.

Результаты подгонки для гранатов Tm3Al5O12, ErTm2Al5O12 и Er2TmAl5O12 представлены на рис. 3a–3в, где показаны отдельные вклады и их сумма. Значения подгоночных параметров приведены в табл. 2. Как видно из рис. 3, вклады Шоттки доминируют при температурах ниже 10 K. В решеточной теплоемкости доминирует вклад Дебая, что следует из сопоставления коэффициентов ${{r}_{{\text{D}}}}$ и ${{r}_{{\text{E}}}}$ в табл. 2. Однако для согласования с экспериментальными данными вклады Эйнштейна необходимо учитывать даже при температурах меньших 45 K.

Рис. 3.

Аппроксимация теплоемкости в гранатах Tm3Al5O12 (a), ErTm2Al5O12 (б) и Er2TmAl5O12 (в). Кружки – экспериментальные данные. Сплошные кривые показывают полную рассчитанную теплоемкость (1), вклад Дебая (2), вклады Эйнштейна (3, 4) и аномалии Шоттки, индуцированные ионами Tm (5) и Er (6). Вставки демонстрируют доминирующую роль аномалий Шоттки при низких температурах.

Таблица 2.  

Коэффициенты ${{r}_{{\text{D}}}}$,${{r}_{{{\text{E}}1}}}$,${{r}_{{{\text{E}}2}}}$; температуры Эйнштейна ${{\theta }_{{{\text{E}}1}}}$ и ${{\theta }_{{{\text{E}}2}}}$; расщепление основного уровня энергии ${{\Delta }_{{{\text{д о п }}}}}$ и подгоночные значения энергий первого возбужденного уровня ионов эрбия ${{E}_{{{\text{1Er}}}}}$ и тулия${{E}_{{{\text{1}}\,{\text{Tm}}}}}$ для гранатов Er3– xTmxAl5O12

x ${{r}_{{\text{D}}}}$ ${{\theta }_{{{\text{E}}1}}},{\text{K}}$ ${{r}_{{{\text{E}}1}}}$ ${{\theta }_{{{\text{E}}2}}},{\text{K}}$ ${{r}_{{{\text{E}}2}}}$ ${{\Delta }_{{{\text{д о п }}}}},{\text{K}}$ ${{E}_{{{\text{1Er}}}}},{\text{K}}$ ${{E}_{{{\text{1}}\,{\text{Tm}}}}},{\text{K}}$
0 54.0 125 1.2 163 4.8 1.1 39
1 53.4 80 0.6 168 6.0 1.0 40 51
2 53.4 87 0.6 164 6.0 0.7 40 51
3 53.0 149 6.0 215 1.0 51

Приложение внешнего магнитного поля приводит к расщеплению вырожденных штарковских уровней и сдвигу синглетов. Для ионов эрбия магнитное поле расщепляет крамерсовские дублеты пропорционально полю. За счет индуцированного внешним полем расщепления увеличивается подъем теплоемкости гранатов, содержащих эрбий, при понижении температуры (см. рис. 2а, 2б). При 2 Tл проявляется пик теплоемкости около 3 K, который согласуется с Шоттки-подобной аномалией, наблюдаемой в смешанных алюмо-иттрий-эрбиевых гранатах [5, 6]. Пик теплоемкости сдвигается к более высоким температурам с ростом поля и значительно размывается. Чувствительность к магнитному полю синглетных уровней ионов Tm3+ обусловлена тем, что поле смешивает волновые функции 3H6 мультиплета [24]. В результате слабовыраженные аномалии Шоттки в тулиевом гранате также сдвигаются к высоким температурам при возрастании магнитного поля (рис. 4).

Рис. 4.

Температурные зависимости энтропии SВ (сплошные кривые) и магнитной энтропии Sm (штриховые кривые) для гранатов Tm3Al5O12 (а), ErTm2Al5O12 (б) и Er2TmAl5O12 (в) во внешнем магнитном поле до 9 Тл.

Данные для теплоемкости можно использовать для расчета энтропии и выделения магнитного вклада. Энтропия $S{}_{B}(T)$ в магнитном поле B при температуре T записывается в виде

(4)
${{S}_{B}}(T) = \int\limits_0^T {\frac{{{{C}_{B}}(t)}}{t}} dt,$
где CB(t) – теплоемкость в поле B при температуре t. Магнитный вклад в энтропию рассчитывается по формуле $S{}_{m}(T) = {{S}_{B}}(T) - {{S}_{0}}(T).$ Экспериментальные результаты, показанные на рис. 2в, дают возможность непосредственно рассчитать изменение энтропии при возрастании температуры от минимального значения T0 ~ 1.9 K, достижимого в экспериментах, до T для граната Tm3Al5O12. Поскольку теплоемкость ниже T0 очень мала и следовательно вклад в энтропию ниже этой температуры также мал, то рассчитанное изменение энтропии дает хорошее приближение для полной энтропии в соотношении (4). Температурные зависимости энтропии при нескольких значениях внешнего поля и соответствующие магнитные вклады в энтропию для Tm3Al5O12 показаны на рис. 4а до 50 K.

Ситуация в гранатах Er3– xTmxAl5O12 (x = 1, 2) значительно отличается. Теплоемкость в нулевом и малых полях возрастает с понижением температуры при низких температурах (см. рис. 1, 2а, 2б) и нельзя пренебречь энтропией ниже T0. В этом случае можно рассчитать энтропию, используя сделанные выше теоретические оценки для вкладов в теплоемкость. Результаты расчета представлены на рис. 4б и 4в.

Обычная процедура охлаждения методом адиабатического размагничивания включает изотермическое снижение энтропии за счет приложения внешнего поля и адиабатическое охлаждение при понижении поля. Такая процедура предполагает уменьшение энтропии с ростом магнитного поля. Согласно рис. 4–6 , все исследованные в данной работе гранаты удовлетворяют этому требованию. Однако для чистого тулиевого граната магнитокалорический эффект становится очень слабым выше 10 K. Величину магнитной энтропии в кристаллах Er3– xTmxAl5O12 можно сравнить с известными данными для магнитокалорических материалов, наиболее часто применяемых в методе магнитного охлаждения [8, 9, 25].

Из сравнения полученных результатов видно, что магнитная энтропия в смешанных алюмо-эрбий-тулиевых и в чистом тулиевом гранатах имеет такой же порядок, как и в галлий-диспрозиевом гранате, сульфате хрома-калия, редкоземельных алюминатах или в интерметаллике YbPt2Sn [22], хотя и меньше в несколько раз, чем в фосфате гадолиния [26] и в галлий-гадолиниевом гранате [25].

ЗАКЛЮЧЕНИЕ

Получены температурные зависимости теплоемкости в монокристаллах гранатов Er3– xTmxAl5O12 (x = 0, 1, 2, 3). Эти зависимости были аппроксимированы суммой вклада Дебая, двух вкладов Эйнштейна и вкладами Шоттки, обусловленными трехвалентными ионами эрбия и тулия. При расчете аномалии Шоттки от ионов эрбия учтено расщепление нижнего крамерсовского дублета магнитными полями, создаваемыми соседними парамагнитными ионами. Показано, что аномалии Шоттки доминируют при низких температурах. Магнитные поля до 9 Tл сдвигали аномалии Шоттки и вызывали их размытие. Магнитная энтропия в исследованных гранатах имела такой же порядок, как и в известных магнитокалорических материалах. Таким образом, гранаты Er3 –xTmxAl5O12 могут быть использованы в методе магнитного охлаждения.

БЛАГОДАРНОСТИ

Авторы благодарны А.А. Каминскому за предоставленные образцы смешанных составов эрбий-тулиевых алюмогранатов.

ФИНАНСИРОВАНИЕ РАБОТЫ

Работа поддержана Министерством образования и науки РФ в рамках Государственного задания на 2019 г. и выполнена при частичной финансовой поддержке Российского фонда фундаментальных исследований (гранты № 18-07-00191, 16-07-00181).

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

  1. Koechner W. Solid State Laser Engineering. Berlin: Springer, 2006.

  2. Каминский А.А. Лазерные кристаллы. М.: Наука, 1975.

  3. Шевченко Е.В., Чарная Е.В., Хазанов Е.Н. и др. // ФТТ. 2017.Т. 59. № 4. С. 717.

  4. French J.D., Zhao J., Harmer M.P. et al. // J. Amer. Ceram. Soc. 1994. V. 77. № 11. P. 2857.

  5. Kushino A., Aoki Y., Yamasaki N.Y. et al. // J. Appl. Phys. 2001. V. 90. № 11. P. 5812.

  6. Shevchenko E.V., Charnaya E.V., Khazanov E.N. et al. // J. Alloys and Compounds. 2017. V. 717. P. 183.

  7. Li R., Numazawa T., Hashimoto T. et al. Advances in Cryogenic Engineering Materials / Eds. by Timmerhaus K.D., Fast R.W., Clark A.F., Reed R.P. N.Y.: Springer US, 1986. P. 287.

  8. Gschneidner Jr.K.A., Pecharsky V.K., Tsokol A.O. // Rep. Prog. Phys. 2005. V. 68. № 6. P. 1479.

  9. Tishin A.M., Bozkova L.P. // J. Appl. Phys. 1997. V. 81. № 2. P. 1000.

  10. Veissier L., Thiel C.W., Lutz T. et al. // Phys. Rev. B 2016. V. 94. № 20. P. 205133.

  11. Хазанов Е.Н., Таранов А.В., Шевченко Е.В., Чарная Е.В. // ЖЭТФ. 2015. Т. 148. № 1(7). С. 56.

  12. Портной К.И., Тимофеева Н.И. Кислородные соединения редкоземельных элементов М.: Металлургия, 1986.

  13. Tari A. The Specific Heat of Matter at Low Temperatures. L.: Imperial College Press, 2003.

  14. Gutowska M.U., Wieckowski J., Szewczyk A. et al. // J. Alloys and Compounds. 2016. V. 670. P. 175.

  15. Babin V., Chernenko K., Hanus M. et al. // Phys. Stat. Sol. B. 2017. V. 254. № 4. P. 1600570.

  16. Shevchenko E.V., Charnaya E.V., Lee M.K. et al. // Phys. Lett. A 2017, 381, 330.

  17. Багдасаров К.С., Додокин А.П., Сорокин А.А. // ФТТ. 1988. Т. 30. № 6. С. 1840.

  18. Perlov D.D., Sorokin A.A., Feodorov V.A., Dodokin A.P. // J. Mag. Mag. Mater. 1993. V. 123. № 1–2. P. 187.

  19. Burdick G.W., Gruber J.B., Nash K.L. et al. // Spectroscopy Lett. 2010. V. 43. № 5. P. 406.

  20. Gruber J.B., Hills M.E., Macfarlane R.M. et al. // Phys. Rev. B 1989. V. 40. № 14. P. 9464.

  21. Tiseanu C., Lupei A., Lupei V. // J. Phys.: Cond. Matt. 1995. V. 7. № 44. P. 8477.

  22. Beghi M.G., Bottani C.E., Russo V. // J. Appl. Phys. 2000. V. 87. № 4. P. 1769.

  23. Huang Z., Feng J., Pan W. // Solid State Sci. 2012. V. 14. № 9. P. 1327.

  24. Валиев У.В., Gruber J.B., Гапдукхаков И.Р. и др. // Оптика и спектроскопия. 2009. Т. 106. № 6. С. 937.

  25. Jang D., Gruner T., Steppke A. et al. // Nature Commun. 2015. V. 6. P. 8680.

  26. Palacios E., Rodrígues-Velamazán J.A., Evangelisti M. et al. // Phys. Rev. B. 2014. V. 90. № 21. P. 214423.

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