ОБЗОР
ФИЗИКА КОНДЕНСИРОВАННОГО СОСТОЯНИЯ ВЕЩЕСТВА
Квантовые основные состояния нитратов меди
О. С. Волкова 1'2'a, А. Н. Васильев1^, В. В. Ховайло3,4
1 Московский государственный университет имени М.В. Ломоносова, физический факультет, кафедра физики низких температур и сверхпроводимости.
Россия, 119991, Москва, Ленинские горы, д. 1, стр. 2.
2 Уральский федеральный университет имени первого Президента России Б. Н. Ельцина, физико-технологический институт, кафедра теоретической физики и прикладной математики.
Россия, 620002, Екатеринбург, ул. Мира, д. 19.
3 Национальный исследовательский технологический университет «МИСиС».
Россия, 119049, Москва, Ленинский пр-т, д. 4.
4 Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики. Россия, 197101, Санкт-Петербург, Кронверкский пр., д. 49.
E-mail: a [email protected], b [email protected] Статья поступила 31.05.2014, подписана в печать 29.07.2014.
Ключ к пониманию фундаментальных свойств материи лежит в исследованиях материи при низких температурах. В условиях, когда тепловые колебания не скрывают тонкие взаимодействия, открывается поле квантовых кооперативных явлений, не имеющих аналогов в классической физике. Именно эти явления — сверхпроводимость, низкоразмерный экзотический магнетизм, волны спиновой и зарядовой плотности, бозе-эйнштейновская конденсация — наиболее актуальны в физике конденсированного состояния. Изучение квантовых основных состояний материи проводится на примере нитратов переходных металлов, в которых магнитоактивные ионы формируют каркасы пониженной размерности: димеры, цепочки, спиновые лестницы, двумерные плоскости.
Ключевые слова: низкоразмерный магнетизм, нитраты переходных металлов.
УДК: 537.6, 53.043, 538.94. PACS: 75.10.Jm, 75.10.Pq, 75.30.Et.
Введение
Работы по изучению физических свойств новых магнитных материалов направлены на изучение наиболее общих эффектов и взаимодействий, формирующих квантовые основные состояния материи. В рамках единой научной картины мира могут рассматриваться не только проблемы физики конденсированного состояния, например, проблема сверхпроводимости, но и космологии и элементарных частиц (в этом плане полезно сопоставить поле Гинзбурга-Ландау и поле Хиггса-Энгле-ра-Браута). В пределах заявленной проблемы рассматриваются и такие «частные» вопросы, как нахождение общности между квантовыми кооперативными явлениями в конденсированных средах. Квантовые кооперативные явления — магнетизм, сверхпроводимость, волны зарядовой плотности, бозе-эйнштейновская конденсация и т. д. — образуют особый, наиболее актуальный раздел физики конденсированного состояния. Именно квантовые аспекты в поведении материи («запутывание» волновых функций, спин-поляризованный транспорт, экзотическая сверхпроводимость) лежат в основе или предполагаются к использованию в наиболее передовых технологиях. Большие усилия прилагаются для установления фундаментальной общности и глубинных аналогий между низкоразмерным и фрустрированным магнетизмом, с одной стороны, и сверхпроводимостью, с другой стороны. Эти два явления — магнетизм и сверхпроводимость, считавшиеся долгое время анти-
подами, — обнаруживают общие черты в объектах, находившихся ранее вне поля зрения теоретиков и экспериментаторов. Само обнаружение высокотемпературной сверхпроводимости в сложных оксидах переходных металлов, которые исходно являются антиферромагнитными изоляторами, полностью изменило вектор развития физики твердого тела. Интерес стали привлекать «новые магнетики», т. е. вещества с пониженной размерностью магнитной подсистемы и фрустрацией обменного взаимодействия. Стало ясно, что в некоторых из этих систем основным состоянием материи является спиновая жидкость, а свойства этого состояния и его элементарных возбуждений близки к свойствам электронной жидкости в сверхпроводниках. Этим, собственно, и определяется научная значимость заявленной проблемы. В настоящей работе описаны базовые характеристики новых низкоразмерных магнетиков — нитратов меди Си(Ш3)2 • Н20 и (Ш)[Си(Ы03)3].
1. Низкоразмерный антиферромагнетик ^^3)2 • H2O
Изучение тригидрата нитрата меди Си(Ы03)2 • • 2.5 Н2 О сыграло выдающуюся роль в становлении области низкоразмерного магнетизма. В этом соединении, как и в других солях меди, ионы Си2+ имеют спин 5 = 1/2, для которого проведено наибольшее количество теоретических исследований и наиболее ярко выражены квантовые эффекты. Важным для количе-
ственного описания этого соединения явилось то обстоятельство, что малость магнитокристаллической анизотропии меди существенно упростила сопоставление теории с экспериментом. В результате именно на примере Си(М03)2 • 2.5Н20 было достигнуто хорошее соответствие между экспериментальными данными [1-3] и теоретическими представлениями [4, 5] и выработаны основные концепции поведения магнитных систем пониженной размерности [6].
В тригидрате нитрата меди находящиеся в квадратном окружении ионы Си2+ образуют альтернированные изолированные (5 = 1/2) цепочки, основным состоянием которых является спиновый синглет. Экспериментально формирование синглетного основного состояния в термодинамических свойствах Си(М03)2 • 2.5Н20 проявляется в характерных немонотонных температурных зависимостях намагниченности М [1] и теплоемкости С [7] с широкими максимумами при температурах, определяемых масштабом обменного взаимодействия в цепочках. Под воздействием внешнего магнитного поля синглетное состояние в Си(М03)2• 2.5Н20 разрушается [8, 9]. В области пересечения основного (5 = 0) и возбужденного триплетного состояния (5 = 1) в первом критическом поле НС1 происходит упорядочение перпендикулярной магнитному полю спиновой компоненты намагниченности. В интервале магнитных полей НС1 - НС2 реализуется скошенное антиферромагнитное состояние, а при НС2 система претерпевает спин — флип-переход [10, 11]. Формирование дальнего магнитного порядка при разрушении синглетного основного состояния внешним магнитным полем интерпретируется в настоящее время как бозе-эйнштейновская конденсация магнонов [12, 13].
Наряду с устойчивой фазой тригидрата нитрата меди Си(М03)2 • 2.5Н20 существует составленное из тех же компонентов крайне неустойчивое на воздухе соединение — моногидрат нитрата меди Си(М03)2 • Н20 [14]. До последнего времени информация о его физических свойствах практически отсутствовала.
Моногидрат нитрата меди Си(М03)2 • Н20 имеет моноклинную решетку с параметрами а = 6.377(1), Ь = 8.548(1), с = 9.769(1) А в = 100.41(1)°, 2 = 4, пространственная группа Р21/с. В этой структуре атомы меди (II) координированы четырьмя атомами кислорода на расстояниях от 1.928 до 1.986 А, пятый апикальный кислород находится на расстоянии 2.332 А. В структуре Си(М03)2 • Н20 квадратные фрагменты Си04 связаны между собой М03-группами, образуя однородные, с магнитной точки зрения, зигзагообразные цепочки 5 = 1/2, которые связываются в слои через длинные связи Си-0. Вдоль оси а магнитоактивные слои взаимодействуют через водородные связи.
Кристаллическая структура Си(М03)2 • Н20 представлена на рис. 1, а. Если учесть наличие апикального кислорода над плоскостью плакетки, то структура моногидрата нитрата меди может быть представлена как сетка димеров в плоскости Ьс , топология которой представлена на рис. 1, б. Устройство собственно диме-ра на ионах Си2+ показано на рис. 1, в. Вдоль оси а магнитоактивные слои взаимодействуют через водородные связи. Описанная структура моногидрата нитрата меди Си(М03)2 • Н20 принципиально отличается
Рис. 1. а — Структура магнитоактивного слоя Си(Ы03)2 • Н20 в плоскости Ьс в полиэдрическом представлении. Вдоль оси а магнитоактивные слои связаны водородными связями. Овалом выделен магнитный димер, а проекция элементарной ячейки на плоскость Ьс показана сплошной линией. Атомы меди обозначены светлыми сферами и находятся в квадратном кислородном окружении, атомы азота в треугольном кислородном окружении обозначены большими темными сферами. б — Структура магнитного димера в Си(Ы03)2 • Н20. в — Топология магнитной подсистемы Си(Ы03)2 • Н20 в плоскости Ьс. Параметр /1 отвечает обменному взаимодействию в димере, параметр /2 — междимерному взаимодействию
от моноклинной 112/с 1 -структуры тригидрата нитрата меди Си(М03)2• 2.5Н20, основным мотивом магнитной подсистемы которого являются изолированные альтернированные 5 = 1/2-цепочки.
Моногидрат нитрата меди получали из нитратокуп-рата нитрозония (М0)Си(М03)3. Для этого образец (М0)Си(М03)3 растворяли при нагревании в концентрированной азотной кислоте. Полученный темно-зеленый раствор запаивали в толстостенные стеклянные ампулы и помещали в дьюар с теплой водой (60°С). В ходе постепенного остывания воды в дьюаре наблюдалось образование и рост голубых пластинчатых кристаллов моногидрата нитрата меди. Монокристаллы моногидрата нитрата меди Си(М03)2 • Н20 представляли собой прозрачные тонкие пластины голубого цвета, нормаль к плоскости которых была ориентирована вдоль оси а. Поперечные размеры пластин составляли единицы миллиметров, их толщина была на порядок меньше.
Намагниченность Си(М03)2• Н20 измерялась в магнитном поле 1 Тл, направленном либо в плоскости Ьс, либо вдоль оси а . При измерениях в плоскости Ьс магнитная восприимчивость хЬс моногидрата нитрата меди в широком интервале температур, как показано на рис. 2, описывается суммой независящего от температуры слагаемого х0 = 3 • 10_5 см _3моль_1 и гиперболического слагаемого хcw = N^2ц/^в(Т — В) с положительной температурой Вейса В = 6 К и эффективным магнитным моментом /лец= 1.97^в. При низких
0.5 0.4 т 0.3
ив §
г
% 0.2 0.1 0.0
Л . ч® 80
1 „ 60 § у
\ 8
1 § 40
- \ ^ \ 20 & ¿Г
■ т » N 1, - * П
0 10 20 30 40 **»»<.»» Г>К
10
20 30
Г, К
40
50
Рис. 2. Температурная зависимость магнитной восприимчивости моногидрата нитрата меди Си(Ы03)2 • Н20 в поле 1 Тл. На вставке показана температурная зависимость обратной магнитной восприимчивости этого же образца
температурах зависимость х(Г) заметно отклоняется от закона Кюри-Вейса в сторону меньших значений, что указывает на «включение» антиферромагнитных взаимодействий.
Магнитное поле напряженностью 1 Тл, как будет показано ниже, существенно изменяет параметры магнитной подсистемы Си(Ы03)2 • Н20. По этой причине на зависимости х(Г), представленной на рис. 2, не наблюдается характерного для антиферромагнетиков пика восприимчивости при температуре Нееля. Для выявления собственного поведения моногидрата нитрата меди были проведены представленные на рис. з измерения его магнитной восприимчивости в слабых магнитных полях при низких температурах. При = 3.6 К на зависимостях хьс (Г) наблюдается максимум. Дополнительная аномалия видна на этих зависимостях при температуре Г * = 2.6 К. Эта аномалия практически полностью подавляется слабым магнитным полем ~ 5 • 10~2 Тл. На вставке к рис. 3 показана температурная зависимость магнитной восприимчивости Ха(Г),
0.45
0.40
0.35
0.30
0.25
0.20
-----
/
V
0.38
т* 0.37
г
г 0.36
0.35
ж
I \
' т„ \
'ооосоооооооо
т*
3.2 Г, К 4.0
2.0
2.5
3.0
г, к
3.5
4.0
Рис. 3. Температурные зависимости магнитной восприимчивости моногидрата нитрата меди Си(Ы03)2 • Н20 в плоскости Ьс в полях 0.01 Тл (красные), и 0.1 Тл (квадраты) при низких температурах. На вставке дана температурная зависимость магнитной восприимчивости вдоль оси а в поле 0.01 Тл (полые кружки)
при Г * также наблюдается слабовыраженная аномалия. В целом спад магнитной восприимчивости от до низких температур существенно более выражен при ориентации магнитного поля в плоскости ьс , нежели при ориентации его вдоль оси а.
Полевые зависимости намагниченности М моногидрата нитрата меди при ориентации магнитного поля В вдоль оси а ив плоскости Ьс приведены на рис. 4. Эти зависимости выходят на разные уровни насыщения при несколько различающихся значениях внешнего магнитного поля. На вставке показаны начальные участки зависимостей Ма(Н) и МЬс (Н), измеренных при Г = 1.8 К. Видно, что зависимость Ма(Н) практически линейна, тогда как зависимость МЬс (Н) демонстрирует в слабых полях сублинейное поведение.
0.6 0.8 1.0 В, Тл
Рис. 4. Полевые зависимости намагниченности моногидрата нитрата меди Си(Ы03) 2• Н2 0 в плоскости ВС (синие кружки) и вдоль оси а (полые кружки) при Г = 1.8 К. На вставке — начальные участки этих зависимостей
Особенности магнитного отклика моногидрата нитрата меди Си(Ы03)2 • Н20 на внешнее магнитное поле наиболее ярко проявляются на представленных на рис. 5 полевых зависимостях производной намагниченности dM/dH при низких температурах. На зависимостях (1МЬс/dH от Н видны аномалии как в слабых
Рис. 5. Производные полевых зависимостей намагниченности моногидрата нитрата меди Си(Ы03)2 • Н20 в плоскости Ьс при вариации температуры от 1.8 до 3.2 К. На вставке показана производная полевой зависимости намагниченности, измеренная вдоль оси а
магнитных полях HC1 ~ 0.05 T, так и в сильных магнитных полях HC2 ~ 1 T. В то же время на зависимости dMa/dH от H, показанной на вставке к рис. 5, аномалия в слабых полях практически не видна, а аномалия в сильных полях выражена столь же ярко.
При измерениях теплоемкости образцы Cu(NO3)2 • • H2O извлекались из стеклянных ампул и монтировались на измерительной платформе PPMS Quantum Design. В этом случае образцы находились на воздухе в течение около одной минуты, что приводило к их частичному превращению в устойчивый тригидрат нитрата меди. Полная гидратация описанных образцов происходила в течение примерно десяти минут, что наблюдалось путем мониторинга их массы.
Температурная зависимость теплоемкости моногидрата нитрата меди Cu(NO3)2 • H2O, представленная на рис. 6, демонстрирует Л -типа аномалию при TN ~ 3.6 К, а также избыточный вклад в теплоемкость типа аномалии Шоттки при низких температурах. При измерениях на разных образцах вклады в теплоемкость от этих двух аномалий менялись относительно друг друга, так что при уменьшении величины Л-пика увеличивалась амплитуда аномалии Шоттки. При низких температурах теплоемкость C исследованного образца может быть представлена суммой двух слагаемых C1 = k1T3 и C2 = k2(A/T)2 ехр(-Д/Т). Первое из этих слагаемых описывает сумму фононного и магнонного вкладов в теплоемкость Cu(NO3)2 • H2O, а второе — вклад в теплоемкость типа аномалии Шоттки, обязанной, по-видимому, примесной фазе Cu(NO3)2 • 2.5H2O. Аппроксимация примесного вклада в теплоемкость зависимостью C2(T) при Д' = 3.8 К показана линией на рис. 6. Измеренная в работе [7] теплоемкость примесной фазы Cu(NO3)2 • 2.5H2O, нормированная на ее предполагаемое содержание в исследованном образце, также показана на рис. 6. Обработка экспериментальных данных, полученных в работе [7], зависимостью типа C2(T) приводит к значению Д" = 4.6 К. Различие в значениях Д' и Д" может быть связано с тем, что
1
- 1 . I 1
Л 1 / - ; \
I /
/ N •••...... •• •■
---------
0 2 4 б 8 10
Г, К
Рис. 6. Температурная зависимость теплоемкости моногидрата нитрата меди Си(Ы03)2 • Н20 (черные кружки). Сплошной кривой показан вклад аномалии типа Шоттки в теплоемкость исследованного образца. Для сравнения приведена нормированная теплоемкость тригидрата нитрата меди Си(Ы03)2 • 2.5Н20 (полые кружки)
аномалия типа Шоттки в исследованном образце наблюдается вблизи фазового перехода в Си(Ы03)2• Н20, где зависимость С1(Т) уже заметно отклоняется от кубической.
Вся совокупность представленных экспериментальных данных свидетельствует о том, что удаление 1.5 Н2О из тригидрата нитрата меди Си(Ы03)2 • • 2.5Н2 0 приводит к формированию крайне неустойчивого к поглощению воды соединения Си(Ы03)2 • Н20 с принципиально отличными от устойчивого соединения магнитными свойствами. В тригидрате нитрата меди Си(Ы03)2 • 2.5Н20 ионы Си2+ координированы четырьмя ионами О2~, два из которых принадлежат нитратным группам Ы03, а два других — молекулам воды Н2О. При дегидратации тригидрата нитрата меди удаляются молекулы воды, расположенные между магнитоактивными слоями, а также по одной молекуле воды у каждого катиона меди, участвующей в организации основного обменного взаимодействия. В результате этого происходит перестройка структуры магнитоактивного слоя таким образом, что ионы Си2+ оказываются координированными по-прежнему четырьмя ионами О2~, но в моногидрате нитрата меди Си(Ы03)2 • Н20 три из них принадлежат нитратным группам Ы03, а один — молекуле воды Н2О.
Аппроксимация высокотемпературной магнитной восприимчивости Си(Ы03)2 • Н20 законом Кюри-Вей-са дает положительную температуру Вейса, что свидетельствует о доминировании в системе ферромагнитного обменного взаимодействия. При низких температурах, однако, зависимость х_1(Т) существенно отклоняется от линейной, указывая на проявление антиферромагнитного обменного взаимодействия. При Тм = 3.6 К в Си(Ы03)2 • Н20, как следует из измерений теплоемкости (см. рис. 6) и намагниченности (см. рис. 3), устанавливается дальний антиферромагнитный порядок. При измерениях в магнитном поле, направленном вдоль оси а, магнитная восприимчивость слабо изменяется с температурой. В то же время при измерениях в магнитном поле, ориентированном в плоскости Ьс, магнитная восприимчивость сильно спадает с понижением температуры. Это позволяет предположить, что магнитные моменты ионов Си2+ направлены в плоскости Ьс. Наличие дополнительной аномалии в магнитной восприимчивости хЬс(Т) при Тзк = 2.7 К может интерпретироваться как проявление спиновой переориентации в плоскости Ьс. Надо полагать, что магнитокристаллическая анизотропия в плоскости Ьс весьма мала, поскольку пик хЬс(Т) при Тзк подавляется слабым магнитным полем. Слабость этого магнитного перехода проявляется также в том, что на температурной зависимости теплоемкости при Тзк не наблюдается никаких особенностей. Из полевых зависимостей намагниченности Си(Ы03)2 • Н20 вытекает, что в магнитном поле, ориентированном в плоскости Ьс, магнитная подсистема моногидрата нитрата меди испытывает последовательные спин-флоп и спин-флип превращения. В то же время в магнитном поле, ориентированном вдоль оси а, проявляется только спин-флип переход. Фазовая диаграмма моногидрата нитрата меди Си(Ы03)2 • Н20, построенная на основа-
нии измерений температурных и полевых зависимостей намагниченности Мьс, представлена на рис. 7.
Г, К
Рис. 7. Магнитная фазовая диаграмма моногидрата нитрата меди Си(ЫОэ)2 • Н2О в магнитном поле в плоскости Ьс. Точки на графике, обозначенные незакрашенными кружками, получены из измерений кривых намагничивания при фиксированной температуре; точки, обозначенные закрашенными кружками, получены из измерений температурных зависимостей намагниченности в фиксированном магнитном поле. Области существования парамагнитной, антиферромагнитных и спин-флоп-фазы обозначены РМ, ЛИ 1, ЛИ 2, БР
При отсутствии прямых исследований магнитной структуры Си(МОз)2 • Н2О с помощью рассеяния нейтронов можно лишь предположить, что магнитные моменты ионов Си2+ ориентированы в плоскости Ьс, будучи ферромагнитно связанными между собой. Вдоль оси а магнитоактивные слои связаны, по-видимому, антиферромагнитным взаимодействием. Спин-флоп и спин-флип-превращения тем самым описывают эволюцию антиферромагнитной подсистемы Си(ЫО3)2 • Н2О, подрешетки которой принадлежат различным слоям. В этой модели спин-переориентационный переход обязан спонтанному изменению направления вектора намагниченности в магнитоактивном слое.
2. Модель «флаг конфедерата» и основное состояние в (NO)[Cu(NO3)3]
Концепция резонирующих валентных связей [15] является исключительно важной в физике конденсированного состояния вещества не только для фрустриро-ванного магнетизма, но также в контексте высокотемпературной сверхпроводимости соединений меди [16]. Интересной моделью в этом плане является фрустри-рованная 5 = 1/2 квадратная решетка. В зависимости от отношения между обменом /1 с ближайшим соседом и обменом /2 со следующим за ближайшим соседом эта модель имеет антиферромагнитное основное состояние при слабой фрустрации и страйп или коллинеарное антиферромагнитное состояние при сильной фрустрации. Имеется узкая область между двумя фазами в диапазоне 0.4 < /2//1 < 0.6, где предполагается отсутствие дальнего магнитного порядка. Вместо этого может быть реализовано состояние спиновой жидкости или кристалл на валентных связях. Поиск экспериментальной
реализации таких моделей интенсивно продолжается среди оксидов меди и ванадия [17]. Альтернатива предлагаемым /1 -/2 моделям недавно была предложена Нерсесяном и Цвеликом (и названа моделью «флаг конфедерата») [18]. Она отличается от модели /1 -/2 пространственной анизотропией обмена с ближайшими соседями /-/' вдоль горизонтального и вертикального направлений и также наличием фрустрирующего обмена /2 вдоль диагоналей, как показано на рис. 8, а. Модель особенно интересна для отношения /'//2 = 2, где основное состояние является резонирующей валентной связью в анизотропном пределе, / » /' = 2/2 [18]. Этот результат представляется альтернативой основного состояния типа кристалла на валентных связях [19-22]. В любом случае отношение /'//2 = 2 ведет к исчезновению эффективного встроенного поля и делает теорию среднего поля взаимодействующих цепочек неприменимой [23]. Хотя, по-видимому, это условие требует некоторой подстройки взаимодействий, мы покажем, что оно практически точно реализовалось в купрате нитрозония (ЫО)[Си(ЫО3)3].
б
Рис. 8. а — Схематичное представление анизотропной модели «флаг конфедерата». б — Схематичное представление кристаллической структуры (ЫО)[Си(ЫО3)3]. /' = 2/2 . Большими сферами показаны катионы Си2+. Гантелями показаны группы ЫО +. Анионы показаны наклоненными плоскими треугольниками
Монокристаллы (ЫО)[Си(ЫО3)3] были получены посредством метода кристаллизации из азотнокислых растворов над осушителем [24]. Однофазная природа полученных образцов была подтверждена идентичностью экспериментальных дифрак-тограмм и теоретических, рассчитанных из рентгеновских данных для монокристаллического образца. Монокристаллы (ЫО)[Си(ЫО3)3] размером (3-6) х (1.5-2.5) х (0.5-1) мм в форме продолговатых пластин нестабильны на воздухе и могли быть безопасно исследованы только в запаянных стеклянных ампулах.
Кристаллическая структура (ЫО)[Си(ЫО3)3] представлена слабо связанными слоями, структура которых показана на рис. 8, б. Предположительно самое сильное взаимодействие / между ионами Си2+ (5 = 1/2) обеспечивается через группы , формируя бесконечные горизонтальные цепочки вдоль оси Ь . Эти цепочки связаны через ионы и ЫО + в плоскости Ьс
таким образом, что вертикальное обменное взаимодействие вдоль оси с в два раза больше, чем обменное взаимодействие вдоль диагонали; /2: имеется два сим-
метричных суперобменных пути, дающих вклад в I', в то время как только один путь (и эквивалентный по симметрии) дает вклад в 12. Предполагается, что межплоскостная связь вдоль оси а слаба и не фруст-рирована. Два эквивалентных пути обменного взаимодействия вдоль этой оси проходят через ЫОз-группу, а расстояние через апикальный кислород Си-О 2.539 А. Путь одиночного обменного взаимодействия Си-О через ЫО3 группы составляет 1.985 А. Слабость межплоскостной связи следует из магнитной неактивности -орбитали, ориентированной вдоль оси а. Тем самым топология магнитной подсистемы является реализацией модели «флаг конфедерата» I > I' = 212 (ср. рис. 8, а).
Температурная зависимость магнитной восприимчивости х (НО)[Си(ЫО3)3], измеренная в магнитном поле 0.1 Тл, ориентированном в плоскости Ьс, показана на рис. 9. При понижении температуры магнитная восприимчивость сначала увеличивается, проходит через широкий максимум и затем снова увеличивается, демонстрируя кюри-вейсовское поведение. Широкий максимум х(Т) может рассматриваться как признак низко-размерности магнитной подсистемы (ЫО)[Си(ЫО3)3]. Особенность при низких температурах значительно подавляется магнитным полем 5 Тл. На вставке к рис. 9 показана полевая зависимость намагниченности, измеренная при 2 К. Зависимость М от Н суперлинейна и могла бы быть представлена суммой линейного вклада предположительно от матрицы (ЫО)[Си(ЫО3)3] и нелинейного вклада. Обработка кюри-вейсовского вклада функцией х = Х0 + С/Т, где температурно-неза-висящий вклад х0 = 2.71 • 10~3 см ~3моль_1, константа Кюри С = 0.008 К/(см3моль) определяет содержание примесных центров С = ng2S(S + 1)Лд/3£в как п = 0.02 (для g -фактора g = 2.06, по данным ЭПР). Тогда как обработка кривой намагничивания, а именно полученное значение момента насыщения М8 = 0.007^в/форм. ед., позволяет определить содержание примесных центров из уравнения М$ = ngS^B
как n = 0.007, что сильно отличается от значения, полученного из обработки температурной зависимости. К тому же источник этих дефектов не ясен, так как метод приготовления исключает присутствие в структуре любых других катионов, кроме Cu2+ и NO +, и любых других анионов, кроме NO^T, в то время как имеется высокое оптическое качество кристаллов. Тем самым сложно объяснить рост температурной зависимости магнитной восприимчивости наличием невзаимодействующих парамагнитных центров. Возможно, система (NO)[Cu(NO3)3] находится вблизи некоторого магнитного упорядочения. Тогда разностную кривую, полученную вычитанием линейного хода из общей намагниченности, можно отнести за счет поведения ферромагнитно скоррелированных магнитных моментов и обработать модифицированной функцией Бриллюэна M(B) = Ms th(NMsB/kBNAT), где N — число ферромаг-нитно скоррелированных спинов [25]. Из проведенной обработки, показанной сплошной линией на рис. 9, число скоррелированных спинов составляет порядка 100, и момент насыщения 0.015^в/форм. ед.
Для того чтобы изучить спиновую восприимчивость матрицы, были выполнены измерения электронного парамагнитного резонанса (ЭПР) монокристаллических образцов для двух ориентаций внешнего магнитного поля: ориентации в плоскости, параллельной плакеткам CuO4, и ориентации перпендикулярной плакеткам CuO4, У и ± соответственно. Для обоих направлений ЭПР спектры содержали одну линию с формой близкой к лоренцевской. Аппроксимация экспериментального сигнала линией Лоренца позволяет точно определить интенсивность сигнала ЭПР IEsr , как показано на рис. 10, ширину линии от пика до пика AHpp и резонансное поле Hres. Тензор g -фактора, рассчитанный из резонансного поля как g = hv/ßBHres, дает величины gll = 2.06 и g± = 2.36. Здесь h — постоянная Планка, /лв — магнетон Бора. Величины полученного g-фак-тора типичны для иона Cu2+ в планарной квадратной координации лигандов [26]. Нужно отметить, что Hres
0.000
100
200
г, к
300
Рис. 9. Температурные зависимости магнитной восприимчивости (ЫО)Си(ЫО3)3, измеренные в разных магнитных полях. На вставке представлена полевая зависимость намагниченности — закрашенные точки, линейный ход, показанный штриховой линией, и сублинейный ход, полученный вычитанием из общей намагниченности линейного хода. Аппроксимация сплошной линией выполнена с помощью модифицированной функции Бриллюэна
1.61.41.2-Я 1.00.8-£-0.6-
о
fS ^^
& ts
s" 0.40.20.0-
/1
ß T=108 К
2.10 2.08 с» 2.06 2.04 2.02 2.40 ^ 2.35 2.30
'» о
Я,-2.06
д±* 2.36
-i-50
-I-100
150 Г, К
О 50 100 150 200 250 300 Г, К
I-
200
250
300
Рис. 10. Температурные зависимости интенсивности ЭПР спектров /еби , нормированные при комнатной температуре на магнитное поле, приложенное параллельно и перпендикулярно плоскости СиО4-плаке-ток. На вставке: температурные зависимости продольной gu и поперечной компонент g-фактора [28]
и соответственно g-факторы практически не зависят от температуры (вставка к рис. 10), тем самым нет никаких указаний для развития внутренних локальных полей из-за конкуренции с квазистатическим ближним или дальним порядком во всем изучаемом температурном интервале. Это подкрепляет гипотезу низко-размерности цепочки гейзенберговских спинов 5 = 1/2 в (ЫО)Си(ЫО3)3 в исследованном температурном интервале.
Интегральная интенсивность ЭПР сигнала /еэн прямо пропорциональна статической восприимчивости спинов, участвующих в резонансе [27]. Этот анализ позволяет выделить вклад от решетки спинов Хэрт в магнитную восприимчивость (ЫО)Си(ЫО3)3. Температурная зависимость /еэн, приведенная к ее величине при 295 К, показана на рис. 10 для || и ± ориентаций магнитных полей. Для обоих направлений полей эти кривые близки и, если сравнивать со статическими магнитными измерениями, как изображено на рис. 11, показывают даже более ярко выраженное низкоразмерное поведение. На рис. 11 сплошная линия была получена суммированием приведенной магнитной восприимчивости, полученной из данных ЭПР, и кюри-вейсовской зависимости, связанной с примесями. Можно видеть, что она достаточно хорошо описывает экспериментальные точки. Однако при самых низких температурах есть некоторые расхождения.
Для получения информации об обменных магнитных взаимодействиях магнитная восприимчивость спинов Хврт была вычислена диагонализацией в модели Нерсесяна-Цвелика при вариации а = ]'/]. Для этого использовались кластеры до 24 спинов с различными геометриями, а именно квадрат 4х4, лестница 8 х 2 и полоска 6 х 4 с периодическими граничными условиями. Из-за влияния эффектов размерности решение для магнитной восприимчивости (сплошные и пунктирная линии на правой панели рис. 12) является
0.006
>4
и
0
1 0.003
0.000
г, к
300
Рис. 11. Температурные зависимости магнитной восприимчивости, измеренной в поле 0.1 Тл (кружки), нормализованной магнитной восприимчивости из данных ЭПР для (ЫО)Си(ЫО3)3 (квадраты). Штриховой линией обозначен кюри-вейсовский примесный вклад
точным только для температур Т > 2Ттах, где Ттах обозначает температуру максимума магнитной восприимчивости. Тем самым довольно сложным оказывается описание низкотемпературной части магнитной восприимчивости. Произведение хтахТтах меньше подвержено влиянию эффектов размерности и более полезно в свете того, что известно значение для этой величины для полностью изолированных цепочек (а = 0) ХтахТтах = 0.0941 А/д^Ав [28]. Когда а = 0, то произведение является функцией а, как показано на левой панели рис. 12. Произведение хтахТтах практически линейно зависит от а, и угол наклона может быть определен как - 0.0558 (для кластера 8 х 2), - 0.0614 (для кластера 4 х 4) и - 0.0607 (для кластера 6 х 4 ), и тем самым само зависит от размерности кластера. Можно заметить, что результат для а = 0 ближе всего к кластеру 8 х 2, который близок по
Рис. 12. На правой панели приведены теоретические расчеты температурных зависимостей магнитной восприимчивости, полученные методом диагонализации (а = /'// характеризует межцепочечное взаимодействие). N характеризует размер кластера, для которого был выполнен расчет (4 х 4, 8 х 2, 6 х 4). На левой панели приведена зависимость хтахТтах/(Ад£2^В/£В) от параметра а. Результат для полностью разделенных
цепочек ( а = 0) показан квадратами [28]
форме к цепочке. Комбинируя эти результаты, получим хтахТтах/{Л^2^В,/кв) = 0.0941 - 0.06а. Используя экспериментально определенное методом ЭПР значение g-фактора g2 = 4.68, получим ХтахТтах = 0.165 -— 0.105а К/(см3моль). Используя экспериментальное значение хтахТтах = 0.163 ± 0.007 К/(см3моль), можно получить оценку —0.05 < а < 0.09.
Проведенный анализ тем самым указывает, что система находится в режиме слабых взаимодействий, и межцепочечные взаимодействия могут быть либо ферро-, либо антиферромагнитными и очень слабыми.
Для случая полностью разделенных цепочек температура максимума определяет магнитный обмен как Ттах = 0.64081, что дает значение I = 170 К [29].
В работе [30] система (ЫО)[Си(ЫО3)3] исследовалась с помощью рамановской спектроскопии. Как показано на рис. 13, левая панель, рамановские спектры, снятые при 295 и 10 К, содержат пики в нескольких интервалах частот. Пики в интервале (70-200 см_1) отвечают внешним колебаниям (поступательным и вращательным модам) группы СиО6 и внутренним колебаниям группы СиО6 (200-600 смнитратной
Рис. 13. Рамановские спектры (ЫО)Си(ЫО3)3, полученные с (хх)-поляризацией при 10 и 295 К. Заштрихованная область отражает температурную эволюцию спектра непрерывного рассеяния. Кристаллическая структура плоскости аЬ (ЫО)Си(ЫО3)3 показана на вставке (левая панель). Рамановские спектры в области маленьких частот с (хи)-поляризацией при разных Т (правая панель) [31]
200 300 400 500 600 Рамановский сдвиг, см1
Рис. 14. Магнитное рамановское рассеяние с конечной энергией максимума Етах, определенной из обработки данных после извлечения всех фононных вкладов. Спектры были сняты при фиксированных температурах в интервале 10-295 К. На вставке представлены температурные зависимости энергии максимума и интегрированной интенсивности [31]
группы ЫО- (700-1500 см-1) и группы нитрозония ЫО + (~2280 см-1). Обработка спектров лоренцевыми функциями позволяет выделить 41 фононную моду в рассматриваемом интервале частот. Нужно заметить, что фононные моды присутствуют в спектре начиная с 70 см-1, что отвечает 100 К, или энергетике обменных магнитных взаимодействий в цепочке. Есть заметная температурная зависимость фононов и непрерывного спектра рассеяния.
На правой панели рис. 13 представлены спектры в области низких частот. Видно, что с температурой довольно сильно изменяются частота, интенсивность и полуширина линий. Так, для внешних и внутренних колебаний группы СиО6 было обнаружено появление новых мод и отклонение от ангармонического приближения ниже 100 К. В целом предполагается наличие динамических изменений дистанций между ионами в решетке (ЫО)Си(ЫО3)3 при низких температурах.
Помимо этого спектры рамановского рассеяния (ЫО)Си(ЫО3)3 при высоких температурах содержат линию с полушириной большей фононного рассеяния, которая показана на рис. 13 закрашенной областью и отдельно на рис. 14. Этот магнитный континуум демонстрирует максимум при энергии Етах ~ 200 см -1. Интенсивность этого континуума резко падает с понижением температуры, как показано на вставке к рис. 14. Однако магнитный континуум перекрывается со многими острыми фононными линиями. Тем самым трудно установить: уменьшилась ли его интенсивность до нуля при конечной температуре. Энергия максимума имеет некоторую сдвижку приблизительно на 10% ниже 100 К, как показано на вставке к рис. 14.
В одномерных цепочках 5 = 1/2 рамановское рассеяние на спиновых возбуждениях может наблюдаться в виде широкого максимума с энергией Етах = 21, где 1 — интеграл обменного взаимодействия. Тем самым полученное значение Етах позволяет определить интеграл обменного взаимодействия магнитного взаимодействия 1 = 150 К. Это хорошо согласуется с данными, полученными из температурной зависимости магнитной восприимчивости.
Сильное уменьшение интенсивности рассеяния для однородных цепочек 5 = 1/2 ниже 100 К не предполагается. Напротив, спиновое рассеяние должно слабо зависеть от температуры, как наблюдалось ранее в низкоразмерной решетке кагоме [31]. Влияние температуры могло бы быть связано с подавлением свободных спинов из-за открытия спиновой щели, как наблюдалось ранее в СиОеО3 [32].
Однако в экспериментах по теплоемкости не было получено никаких доказательств существования спиновой щели в системе (ЫО)Си(ЫО3)3. Зависимость С/Т от Т2, показанная на рис. 15, обнаруживает большой линейный терм 7Т, наложенный на кубический терм /ЗТ3. Хотя решеточный вклад вероятнее всего может быть приписан к фононам, спиновые возбуждения в пределах цепочки могут быть ответственны за первый вклад. Причем при самых низких температурах теплоемкость увеличивается при понижении температуры, будучи незначительно подавленной внешнем магнитным полем.
Рис. 15. Температурная зависимость приведенной теплоемкости (ЫО)Си(ЫО3)3
Рис. 16. Температурные зависимости теплоемкости (ЫО)Си(ЫО3)3, измеренные в разных магнитных полях. Для ясности зависимости, полученные в разных полях АБ = 1 Тл, сдвинуты друг относительно друга. Сплошной линией показана граница существования магнитоупорядоченной фазы. На вставке представлена зависимость С/Т -Т2 при Т > Тм для Б = 0 Тл
В измерениях температурных зависимостей теплоемкости вплоть до низких температур ~ 0.3 К, представленных на рис. 16, было установлено наличие фазового перехода в системе (ЫО)Си(ЫО3)3, который представлен Л — аномалией при = 0.58 ± 0.01 К. При приложении внешнего магнитного поля аномальный максимум смещается, что указывает на магнитную природу перехода. Причем в полях до 6-7 Тл температура аномалии увеличивается вплоть до 0.69 К и затем уменьшается 0.68 К в поле 9 Тл. Похоже,
что потребуются очень сильные магнитные поля для полного подавления магнитного порядка в системе. При отсутствии фононного вклада в прозрачном изоляторе (ЫО)Си(ЫО3)3 всю теплоемкость при Т < Тм можно приписать магнонам. Тем самым интегрирование зависимости С/Т(Т) позволяет определить маг-
ТЛ
нитный вклад в энтропию как Smagn = $ С йТ = 1.96 •
0
• 10~2 Дж/(моль• К), что составляет от теоретического значения Smagn = пК ln(2S + 1) = 5.76 Дж/(моль• К2) величину порядка 0.003К ln(2S + 1) и указывает на низкоразмерность системы, в которой выделение магнитной энтропии произошло в температурной области широкого максимума на зависимости х(Т).
Как показано на вставке рис. 16, при Т > Тм зависимость С/Т от Т2 линейна, что позволяет определить коэффициент 7 при линейном вкладе в теплоемкость, который принадлежит одномерным антиферромагнитным магнонам как 7 = 0.042 Дж/(моль• К). Полученное значение 7 позволяет определить интеграл антиферромагнитного обменного взаимодействия в цепочке I, 7 = 2&В/(3!). Полученное значение составляет 132 К, которое также достаточно хорошо согласуется со значением определенным из температурной зависимости магнитной восприимчивости.
Трудно оценить природу основного состояния (ЫО)[Си(ЫО3)3] с точки зрения настоящих термодинамических данных в этой системе. Нерсесян и Цве-лик аргументировали, что для малого межцепочечного взаимодействия основное состояние остается неупорядоченным и реализуется хиральный п-поток спиновой жидкости при нулевой температуре [18], который фактически не противоречит настоящим результатам. Таким образом, основное состояние (ЫО)[Си(ЫО3)3] является антиферромагнитно упорядоченным. Возможно, основное состояние модели «флаг конфедерата» реализуется в данной системе, однако оказывается нестабильным к дополнительным взаимодействиям, таким, как нефрустрированные внутриплоскостные возбуждения или взаимодействия Дзялошинского-Мория, которые приводят к трехмерному дальнему магнитному порядку при более низких температурах.
Заключение
В заключение подчеркнем, что нитраты переходных металлов как отдельный класс слабоизученных материалов представляют собой замечательный полигон для проверки многих базовых положений физики низкоразмерного магнетизма. Это связано с тем, что в этом семействе соединений реализуется множество ранее не исследованных магнитных топологий, которые допускают существование (и даже сосуществование) ближнего и дальнего магнитного порядка при низких температурах.
Авторы выражают благодарность М. Абдель-Хафизу за участие в проведении физических измерений.
Работа выполнена при финансовой поддержке Мин-обрнауки РФ по программе повышения конкуренто-
способности НИТУ «МИСиС» среди ведущих мировых научно-образовательных центров на 2013-2020 гг.
Список литературы
1. Berger L., Friedberg S.A., Schriempf J.T. // Phys. Rev. 1967. 132. P. 1057.
2. Myers B.E., Berger L., Friedberg S.A // J. Appl. Phys. 1969. 40. P. 1149.
3. Amaya K., Tokunaga Y., Yamada R. et al. // Phys. Lett.
A. 1969. 28. P. 732.
4. Tachiki M, Yamada T. // J. Phys. Soc. Jpn. 1970. 28. P. 1413.
5. Bonner J.C., Friedberg S.A., Kobayashi H. et al. // Phys. Rev. B. 1983. 27. P. 248.
6. Carlin R.L. // Magnetochemistry. Berlin: Springer-Verlag, 1986.
7. Friedberg S.A., Raquet C.A. // J. Appl. Phys. 1968. 39. P. 1132.
8. van TolM.W., Diederix K.M., Poulis N.J. //Physica. 1973. 64. P. 363.
9. Diederix K.M., Groen J.P., Henkens L.S.J.M. et al. // Physica. B. 1978. 94. P. 9.
10. van TolM.W., Henkens L.S.J.M, Poulis N.J. // Phys. Rev. Lett. 1971. 27. P. 739.
11. Diederix K.M., Groen J.P., Klaassen T.O., Poulis N.J. // Physica. B. 1979. 96. P. 41.
12. Matsubara T, Matsuda H. // Prog. Theor. Phys. 1956. 16. P. 569.
13. Sachdev S. // Science 2000. 288. P. 475.
14. Troyanov S.I., Morozov I.V., Znamenkov K.O., Kore-nev Yu.M. // Z. Anorg. Allg. Chem. 1995. 621. P. 1261.
15. Anderson P.W. // Science. 1987. 235. P. 1196.
16. Bednorz J.G., Muller K.A. // Z. Physik 1986. 64. P. 189.
17. Tsirlin A.A., Rosner H. // Phys. Rev. B. 2009. 79. P. 214417.
18. Nersesyan A.A., Tsvelik A.M. // Phys. Rev. B 2003. 67. P. 024422.
19. Sindzingre P. // Phys. Rev. B. 2004. 69. P. 094418.
20. Starykh O.A., Balents L. // Phys. Rev. Lett. 2004. 93. P. 127202.
21. Moukouri S., Alvarez J.V. // cond-mat: 0403372.
22. HungH.-H, Gong C.-D, Chen Y.-C, Yang M.-F. // Phys. Rev. B. 2006. 73. P. 224433.
23. Schulz H.J. // Phys. Rev. Lett. 1996. 77. P. 2790.
24. Znamenkov K.O., Morozov I.V., Troyanov S.I. // Rus. J. Inorg. Chem. 2004. 49. P. 172.
25. Klingeler R., Buchner B., Cheong S.-W., Hucker M. // Phys. Rev. B. 2005. 72. P. 104424.
26. Pilbrow J.R. Transition Ion Electron Paramagnetic Resonance. Oxford: Clarendon Press, 1990.
27. Abragam A., Bleaney B. Electron Paramagnetic Resonance of Transition Ions. L.: Oxford University Press, 1970.
28. Volkova O, Morozov I., Shutov V. et al. // Phys. Rev. B. 2010. 82. P. 054413.
29. Bonner J.C., Fisher M.E. // Phys. Rev. A. 1964. 135. P. 650.
30. Johnston D.C., Kremer R.K., Troyer M. et al. // Phys. Rev.
B. 2000. 61. P. 9558.
31. Gnezdilov V., Lemmens P., Pashkevich Yu.G. et al. // Phys. Rev. B 2012. 85. P. 214403.
32. Wulferding D., Lemmens P., Sheib P. et al. // Phys. Rev. B. 2010. 82. P. 144412.
33. Gros C, Lemmens P., Vojta M. et al. // Phys. Rev. B. 2003. 67. P. 174405.
Quantum ground states of copper nitrates O.S. Volkova12a, A.N. Vasiliev1Ab, V.V. Khovailo34
1 Department of Low-Temperature Physics and Superconductivity, Faculty of Physics, M. V. Lomonosov Moscow State University, Moscow 119991, Russia.
2 Department of Theoretical Physics and Applied Mathematics, Institute of Physics and Technology, Ural Federal University, Ekaterinburg 620002, Russia.
3 National University of Science Technology «MISiS», Moscow 119049, Russia.
4 Saint Petersburg National Research University of Information Technologies, Mechanics and Optics, St. Petersburg 197101, Russia.
E-mail: a [email protected], b [email protected].
Low-temperature study of matter provides a key to understanding its fundamental properties. Under conditions where thermal fluctuations do not mask fine interactions, quantum cooperative phenomena that have no analogy in classical physics come into play. Among them are superconductivity, low-temperature exotic magnetism, spin and charge density waves, and Bose-Einstein condensation, which are the most important issues for modern condensed-matter physics. Study of the quantum ground states is illustrated by the example of transition metal nitrates in which magnetically active ions form lower-dimensionality frameworks, including dimers, chains, spin ladders, and 2D planes.
Keywords: low-dimensional magnetism, transition metal nitrates. PACS: 75.10.Jm, 75.10.Pq, 75.30.Et. Received 31 May 2014.
English version: Moscow University Physics Bulletin 6(2014). Сведения об авторах
1. Волкова Ольга Сергеевна — канд. физ.-мат. наук, доцент; тел.: (495) 939-48-11, e-mail: [email protected].
2. Васильев Александр Николаевич — доктор физ.-мат. наук, профессор, зав. кафедрой; тел.: (495) 939-38-25, e-mail: [email protected].
3. Ховайло Владимир Васильевич — доктор физ.-мат. наук, профессор; тел.: (495) 638-46-20, e-mail: [email protected].