Ультрабыстрая динамика магнитных моментов в объемных ферромагнетиках

Основная информация

• ID статьи: 2502.07375
• Название: Ultrafast dynamics of moments in bulk ferromagnets
• Авторы: Mouad Fattouhi, Pascal Thibaudeau, Liliana D. Buda-Prejbeanu
• Классификация: cond-mat.mes-hall cond-mat.mtrl-sci physics.app-ph
• Дата публикации: 22 мая 2025 г. (препринт arXiv)
• Ссылка на статью: https://arxiv.org/abs/2502.07375

Аннотация

В данной работе разработана надежная и эффективная модель для исследования ультрабыстрой динамики магнитных материалов, возбуждаемых лазерными импульсами. Модель интегрирует динамическое уравнение Ландау-Лифшица-Блоха, квантовую тепловую ванну и двухтемпературную модель. Модель успешно применена к трем типичным материалам из литературы: никелю, кобальту и железу. Кроме того, анализ тензора ультрабыстрой динамической восприимчивости показывает, что недиагональные элементы демонстрируют специфические характеристики в зависимости от наличия или отсутствия постоянного внешнего магнитного поля.

Исследовательский контекст и мотивация

Решаемые фундаментальные проблемы

1. Точный контроль тепловых флуктуаций: Точное управление тепловыми флуктуациями соответствующих физических величин в миниатюрных магнитных устройствах для обеспечения надежной работы остается серьезной проблемой
2. Моделирование ультрабыстрой магнитной динамики: Из-за сложных процессов, включающих различные пространственные и временные масштабы, в настоящее время отсутствует единая модель ультрабыстрой магнитной динамики
3. Вычислительная эффективность: Существующие стохастические уравнения Ландау-Лифшица-Гильберта (sLLG) требуют множественных стохастических реализаций и усреднения, что приводит к низкой вычислительной эффективности

Научная значимость

1. Развитие ультрабыстрой спинтроники: Схемы возбуждения фемтосекундными лазерами способны быстро генерировать спиновые токи, открывая область ультрабыстрой спинтроники
2. ТГц спинтронные устройства: Прокладывает путь для разработки спинтронных устройств в диапазоне ТГц, что является первым практическим применением ультрабыстрой спиновой физики
3. Острая необходимость в теоретическом понимании: Оптимизация и надежный контроль ультрабыстрых магнитных устройств возможны только благодаря теоретическому пониманию

Ограничения существующих методов

1. Классическое уравнение LLB: Требует предварительного знания точных термодинамических функций, включая температурную зависимость обменной жесткости, спонтанной намагниченности и параметра затухания
2. Вычислительные затраты стохастических методов: Уравнение sLLG требует большого количества стохастических реализаций для получения надежных средних значений
3. Проблема теплового связывания: Способ связи температуры электронной системы со стохастической или детерминированной формой динамики атомных спинов остается нерешенной проблемой

Основные вклады

1. Динамическая модель LLB (dLLB): Предложена система уравнений dLLB, основанная на процедуре статистического усреднения, способная описывать динамику первого и второго моментов спиновых переменных
2. Интеграция квантовой тепловой ванны: Объединение квантовой теоремы флуктуации-диссипации с моделью dLLB с использованием статистики Бозе-Эйнштейна вместо классической статистики Максвелла-Больцмана
3. Высокоэффективный детерминированный метод: Вычислительная скорость метода dLLB более чем в 100 раз выше, чем у уравнения sLLG
4. Анализ ультрабыстрой динамической восприимчивости: Первый систематический анализ поведения диагональных и недиагональных элементов тензора ультрабыстрой динамической восприимчивости

Методологические подробности

Определение задачи

Исследование ультрабыстрого размагничивания и переммагничивания ферромагнитных материалов под воздействием лазерных импульсов, включая:

• Входные данные: Параметры лазерного импульса (плотность мощности, длительность импульса, длина волны), параметры материала
• Выходные данные: Эволюция намагниченности во времени, динамика тензора восприимчивости
• Ограничения: Теорема флуктуации-диссипации, сохранение энергии

Архитектура модели

1. Система уравнений dLLB

Для каждого атомного узла i модель dLLB сводится к следующей системе уравнений:

Уравнение динамики первого момента:

(1 + α²) ds_i/dt = ω_i × s_i - α(Σ_i - Tr(Σ_i)1)ω_i - 2D/(1 + α²) s_i

Уравнение динамики второго момента:

(1 + α²) dΣ_i/dt = 2ΩΣ_i - 2α([M^A_i, Σ_i] + Tr(M_i)(Σ_i - 2Γ_i) - (Tr(Σ_i) - 2Tr(Γ_i))M^S_i) - 2D/(1 + α²)(3Σ_i - 1Tr(Σ_i))

где:

• s_i ≡ ⟨S_i⟩: направление среднего магнитного момента
• Σ_i ≡ ⟨S_i ⊗ S_i⟩: тензор второго момента
• α: параметр поперечного затухания Гильберта
• D: амплитуда шума

2. Модель квантовой тепловой ванны

Классическое соотношение флуктуации-диссипации (CFDR):

D = (α/ℏ)k_B T

Квантовое соотношение флуктуации-диссипации (QFDR):

D = (α/ℏ) ∫₀^∞ [ℏω/(exp(ℏω/k_B T) - 1)] g_m(ω, T) dω

где g_m(ω, T) — плотность состояний магнонов.

3. Двухтемпературная модель (2TM)

C_e dT_e/dt = -G_{e-ph}(T_{ph} - T_e) + P(t)
C_{ph} dT_{ph}/dt = G_{e-ph}(T_{ph} - T_e)

Член лазерного возбуждения:

P(t) = (μF)/(ℓ_p τ) exp(-(t-t₀)²/(τ²/4log(2)))

Технические инновации

1. Метод статистического усреднения: Строгая процедура статистического усреднения, основанная на функциональном исчислении, избегающая необходимости предварительного предоставления термодинамических функций
2. Квантовые поправки: Использование четвертого порядка дисперсионного соотношения магнонов и плотности состояний магнонов с параметрами особенности Ван Хова
3. Самосогласованное замыкание: Достижение самосогласованного замыкания иерархии уравнений через уравнение динамики второго момента
4. Детерминированное решение: Избежание проблемы множественных реализаций, требуемых стохастическими методами

Экспериментальная установка

Параметры материалов

Исследованы три типичных 3d ферромагнитных металла:

Параметры моделирования

• Ni и Co: Элементарная ячейка 5×5×5, N=500 спинов атомов
• Fe: Элементарная ячейка 7×7×7, N=686 спинов атомов
• Граничные условия: Периодические граничные условия
• Численный метод: Схема Рунге-Кутты четвертого порядка

Критерии оценки

1. Точность предсказания температуры Кюри: Сравнение с экспериментальными значениями
2. Временные масштабы размагничивания: Характерное время восстановления намагниченности после лазерного импульса
3. Вычислительная эффективность: Ускорение относительно метода sLLG

Экспериментальные результаты

Основные результаты

1. Предсказание температуры Кюри

• Метод CFDR: Систематическое завышение температуры Кюри
• Метод QFDR: Отличное согласие с экспериментальными данными
• Механизм улучшения: Квантовая тепловая ванна захватывает физику магнонов при низких температурах

2. Ультрабыстрое размагничивание

Тестирование трех материалов при различных плотностях мощности лазера:

• Ni: Хорошее согласие с экспериментальными данными в диапазоне 22-50 Дж/м²
• Fe: Точное воспроизведение процессов размагничивания и переммагничивания в диапазоне 15-60 Дж/м²
• Co: Захват основных характеристик динамики в диапазоне 22-50 Дж/м²

3. Повышение вычислительной эффективности

• Ускорение: Более чем в 100 раз быстрее, чем метод sLLG
• Причина: Детерминированный метод избегает необходимости множественных стохастических реализаций

Анализ динамической восприимчивости

Без внешнего магнитного поля

• Диагональные элементы: χ_, χ_, χ_ стремятся к 1/3 в парамагнитной фазе (изотропия)
• Недиагональные элементы: χ_, χ_, χ_ остаются равными нулю

С внешним магнитным полем (B = 1.5T вдоль z)

• Диагональные элементы: Поведение аналогично случаю без поля
• Недиагональные элементы: Появляются пики и релаксационное поведение, отражающие корреляции между пространственными ориентациями

Зависимость от длительности импульса

При различных длительностях лазерного импульса (50 fs, 200 fs, 800 fs):

• Диагональные элементы: Нечувствительны к изменению длительности импульса
• Недиагональные элементы: Амплитуда пика уменьшается с увеличением длительности импульса, ширина увеличивается

Связанные работы

Классификация теоретических моделей ультрабыстрой магнетики

1. Прямое когерентное взаимодействие: Прямое взаимодействие лазерного импульса со степенями свободы спина
2. Локальная динамика спина: Локальная динамика спина, вызванная лазерным нагревом или возбуждением
3. Нелокальный перенос углового момента: Например, сверхдиффузионный спиновый транспорт

Существующие тепловые спиновые модели

• Стохастические методы: Прямое моделирование стохастической формы атомных спинов
• Детерминированные методы: Такие как тепловая макроспиновая модель и уравнение LLB
• Методы атомного масштаба: Атомная динамика спина, основанная на известном гамильтониане

Преимущества данной работы

1. Предсказательность: В отличие от классического уравнения LLB, не требует предварительных знаний
2. Эффективность: Значительное повышение вычислительной скорости по сравнению с методом sLLG
3. Точность: Квантовые поправки обеспечивают более точное описание низкотемпературного и критического поведения

Заключение и обсуждение

Основные выводы

1. Модель dLLB-QFDR: Метод dLLB с квантовой тепловой ванной успешно предсказывает тепловые эффекты в ферромагнитных материалах
2. Предсказание температуры Кюри: Квантовые поправки значительно улучшают предсказание температуры Кюри и кривых намагниченности при низких температурах
3. Ультрабыстрая динамика: Модель точно захватывает основные характеристики ультрабыстрого размагничивания 3d ферромагнитных элементов
4. Вычислительные преимущества: Детерминированный метод имеет значительные преимущества в вычислительной эффективности по сравнению со стохастическими методами

Ограничения

1. Приближение среднего поля: Игнорирование пространственных корреляций может влиять на точное описание критического поведения
2. Упрощенная 2TM: Предположение об одинаковых электронной и спиновой температурах может быть недостаточно точным
3. Диапазон материалов: В настоящее время проверено только на простых 3d ферромагнитных металлах

Направления будущих исследований

1. Расширение применения: Метод может быть распространен на ферримагнетики, антиферромагнетики и многослойные структуры
2. Улучшение 2TM: Учет более сложной температурной зависимости и нелинейных эффектов
3. Пространственные корреляции: Включение высших моментов флуктуаций для улучшения описания критического поведения

Глубокая оценка

Достоинства

1. Теоретическая инновация: Строгий и систематический вывод уравнений dLLB на основе статистического усреднения
2. Квантовые поправки: Введение QFDR значительно улучшает физическую точность модели
3. Вычислительная эффективность: Достигнутое ускорение детерминированного метода имеет важное практическое значение
4. Систематическая верификация: Достаточная проверка на различных материалах и условиях

Недостатки

1. Теоретические приближения: Приближение среднего поля может ограничить точность модели в критической области
2. Зависимость от параметров: Модель требует множество параметров материала, некоторые из которых сложно определить экспериментально
3. Область применимости: В настоящее время применима в основном к простым ферромагнитным металлам, применимость к сложным магнитным материалам требует проверки

Влияние

1. Научный вклад: Предоставляет новый теоретический инструмент и вычислительный метод для ультрабыстрой магнетики
2. Практическая ценность: Эффективный вычислительный метод способствует проектированию и оптимизации устройств
3. Воспроизводимость: Подробное описание методов и параметров облегчает воспроизведение и применение другими исследователями

Сценарии применения

1. Проектирование ультрабыстрых спинтронных устройств: Теоретическое проектирование магнитных устройств в диапазоне ТГц
2. Исследование лазерно-индуцированной магнитной динамики: Фундаментальные исследования взаимодействия фемтосекундных лазеров с магнитными материалами
3. Анализ тепловых эффектов магнитных материалов: Теоретическое предсказание свойств магнитных материалов при конечной температуре

Библиография

Статья цитирует 73 соответствующих источника, охватывающих важные работы в области ультрабыстрой магнетики, спиновой динамики, статистической физики и других областей, обеспечивая прочную теоретическую основу и экспериментальные данные для сравнения.