2025-11-13T11:31:10.165583

Localized fermions on the triangular lattice with Ising-like interactions

Regeciová, Kapcia
The model of localized fermions on the triangular lattice is analyzed in means of the Monte Carlo simulations in the grand canonical ensemble. The Hamiltonian of the system has a form of the extended Hubbard model (at the atomic limit) with nearest-neighbor Ising-like magnetic $J$ interactions and onsite Coulomb $U$ interactions. The model is investigated for both signs of $J$, arbitrary $U$ interaction and arbitrary chemical potential $μ$ (or, equivalently, arbitrary particle concentration $n$). Based on the specific heat capacity and sublattice magnetization analyses, the phase diagrams of the model are determined. For ferromagnetic case ($J<0$), the transition from the ordered phase (which is a standard ferromagnet and can be stable up to $k_{B}T/|J| \approx 0.61$) is found to be second-order (for sufficiently large temperatures $k_{B}T/|J| \gtrsim 0.2$) or first-order (for $-1<U/|J|<-0.65$ at the half-filling, i.e., $n=1$). In the case of $J>0$, the ordered phase occurs in a range of $-1/2<U/|J|<0$ (for $n=1$), while for larger $U$ the state with short-range order is also found (also for $n \neq 1$). The ordered phase is characterized by an antiferromagnetic arrangement of magnetic moments in two sublattices forming the hexagonal lattice. The transition from this ordered phase, which is found also for $μ\neq 0$ ($n \neq 1$) and $U/|J|>-1/2$ is always second-order for any model parameters. The ordered phase for $J>0$ can be stable up to $k_{B}T/|J| \approx 0.06$.
academic

Локализованные фермионы на треугольной решетке с изинговыми взаимодействиями

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

  • ID статьи: 2412.10808
  • Название: Localized fermions on the triangular lattice with Ising-like interactions
  • Авторы: Lubomíra Regeciová (Словацкая академия наук), Konrad Jerzy Kapcia (Университет Адама Мицкевича)
  • Классификация: cond-mat.stat-mech cond-mat.other cond-mat.quant-gas cond-mat.str-el physics.comp-ph
  • Дата публикации: 17 октября 2025 г.
  • Ссылка на статью: https://arxiv.org/abs/2412.10808

Аннотация

В данном исследовании методом Монте-Карло в канонической системе проанализирована модель локализованных фермионов на треугольной решетке. Гамильтониан системы имеет форму расширенной модели Хаббарда (атомный предел) и содержит ближайшие соседние изинговые магнитные взаимодействия J и кулоновское взаимодействие на узле U. Исследование рассматривает оба знака J, произвольное взаимодействие U и произвольный химический потенциал μ (или эквивалентно произвольную концентрацию частиц n). На основе анализа теплоемкости и намагниченности подрешеток определена фазовая диаграмма модели. Для ферромагнитного случая (J<0) переход из упорядоченной фазы может быть фазовым переходом второго рода (при достаточно высоких температурах kBT/|J|≳0.2) или первого рода (при половинном заполнении -1<U/|J|<-0.65). Для случая J>0 упорядоченная фаза появляется в диапазоне -1/2<U/|J|<0 (при n=1), а при больших значениях U обнаружено также состояние с короткодействующим порядком.

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

  1. Основная проблема: Исследование поведения магнитных фазовых переходов локализованных фермионов в геометрически фрустрированных системах, в частности фазовая диаграмма расширенной модели Хаббарда в атомном пределе на треугольной решетке.
  2. Значимость:
    • Треугольная решетка обладает внутренней геометрической фрустрацией, приводящей к невозможности одновременного удовлетворения всех антиферромагнитных взаимодействий
    • Данная модель представляет нетривиальное обобщение классической модели Изинга на фермионные системы
    • Имеет важное значение для понимания сильно коррелированных электронных систем и магнитных материалов
  3. Ограничения существующих методов:
    • Исследования гиперкубических решеток достаточно полны, но точные аналитические результаты для треугольной решетки отсутствуют
    • Приближение среднего поля переоценивает температуру фазового перехода
    • Отсутствуют детальные численные исследования эффектов геометрической фрустрации
  4. Исследовательская мотивация:
    • Заполнить пробел в численных исследованиях данной модели на треугольной решетке
    • Исследовать влияние геометрической фрустрации на поведение фазовых переходов
    • Обеспечить теоретическую основу для экспериментальных материалов (таких как Li2MnTeO6, KCeO2 и др.)

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

  1. Первое систематическое исследование методом Монте-Карло: Полное численное моделирование расширенной модели Хаббарда с изинговыми магнитными взаимодействиями на треугольной решетке
  2. Построение полной фазовой диаграммы: Определена трехмерная фазовая диаграмма в координатах температура-химический потенциал-сила взаимодействия для ферромагнитного (J<0) и антиферромагнитного (J>0) случаев
  3. Анализ природы фазовых переходов: Точное определение границ фазовых переходов первого и второго рода, обнаружено, что все переходы в антиферромагнитном случае являются переходами второго рода
  4. Выявление эффектов фрустрации: Обнаружено, что геометрическая фрустрация значительно влияет на структуру фазовой диаграммы, температура стабильности упорядоченной фазы в антиферромагнитном случае на порядок ниже, чем в ферромагнитном
  5. Идентификация состояния с короткодействующим порядком: Обнаружено и охарактеризовано наличие состояния с короткодействующим порядком в антиферромагнитном случае

Подробное описание методов

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

Исследование поведения фазовых переходов локализованной фермионной системы на треугольной решетке, входные параметры модели (J, U, μ, T), выходные данные - фазовая диаграмма системы и природа фазовых переходов.

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

Гамильтониан: H^=U∑in^i↑n^i↓+2Jz∑⟨i,j⟩s^is^j−(μ+U2)∑in^i\hat{H} = U\sum_i \hat{n}_{i\uparrow}\hat{n}_{i\downarrow} + \frac{2J}{z}\sum_{\langle i,j\rangle} \hat{s}_i\hat{s}_j - \left(\mu + \frac{U}{2}\right)\sum_i \hat{n}_i

где:

  • s^i=(n^i↑−n^i↓)/2\hat{s}_i = (\hat{n}_{i\uparrow} - \hat{n}_{i\downarrow})/2: оператор локального спина
  • n^i=n^i↑+n^i↓\hat{n}_i = \hat{n}_{i\uparrow} + \hat{n}_{i\downarrow}: оператор числа частиц
  • z=6z = 6: координационное число треугольной решетки

Состояния узла: Каждый узел может находиться в четырех состояниях:

  • ∣0⟩|0\rangle: вакансия
  • ∣↑⟩|\uparrow\rangle: одиночное заполнение (спин вверх)
  • ∣↓⟩|\downarrow\rangle: одиночное заполнение (спин вниз)
  • ∣↑↓⟩|\uparrow\downarrow\rangle: двойное заполнение

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

  1. Моделирование в канонической системе: Использование канонической системы с флуктуациями числа частиц для более точного описания системы, контролируемой химическим потенциалом
  2. Анализ конечного размера: Систематическое исследование эффектов конечного размера для систем L = 12×12 до 36×36
  3. Методы определения фазовых переходов:
    • Использование метода функции распределения энергии P(E) по Challa и др. для различения переходов первого и второго рода
    • Комбинированный анализ масштабного поведения теплоемкости и намагниченности подрешеток
  4. Разложение по подрешеткам: Разложение треугольной решетки на три эквивалентные подрешетки для точного анализа магнитной упорядоченности

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

Набор данных

  • Размер системы: Треугольные кластеры, L = 12×12, 24×24, 36×36 узлов
  • Диапазон параметров:
    • Температура: kBT/|J| = 0.005 - 1.0
    • Взаимодействие: U/|J| = -1.0 - 5.0
    • Химический потенциал: μ/|J| = -0.6 - 0 (с использованием симметрии частица-дырка)

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

  1. Теплоемкость: C=⟨ε2⟩−⟨ε⟩2/T2C = \langle\varepsilon^2\rangle - \langle\varepsilon\rangle^2/T^2
  2. Намагниченность подрешеток: mα=(6/L)∑i∈α⟨s^i⟩m_\alpha = (6/L)\sum_{i\in\alpha}\langle\hat{s}_i\rangle
  3. Концентрация частиц: n=⟨n^i⟩n = \langle\hat{n}_i\rangle
  4. Функция распределения энергии: P(E)P(E) для определения природы фазовых переходов

Детали реализации

  • Алгоритм: Алгоритм Метрополиса Монте-Карло
  • Шаги равновесия: 10610^6 шагов
  • Шаги статистики: 10610^6 шагов
  • Сканирование температуры: Постепенное охлаждение от высоких температур к низким
  • Граничные условия: Периодические граничные условия

Результаты экспериментов

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

Ферромагнитный случай (J < 0):

  • Максимальная температура фазового перехода: kBT/|J| ≈ 0.605 (при U → ∞, близко к 0.607 классической модели Изинга)
  • Положение тройной критической точки: kBT/|J| ≈ 0.18-0.20, практически независимо от U и μ
  • Природа фазового перехода: переход первого рода при низких температурах, второго рода при высоких температурах
  • Упорядоченная фаза: стандартная ферромагнитная фаза, все намагниченности подрешеток равны

Антиферромагнитный случай (J > 0):

  • Максимальная температура фазового перехода: kBT/|J| ≈ 0.06 (значительно ниже ферромагнитного случая)
  • Диапазон существования упорядоченной фазы: -1/2 < U/|J| < 0 (при половинном заполнении)
  • Природа фазового перехода: все переходы являются переходами второго рода
  • Структура упорядоченной фазы: две подрешетки образуют антиферромагнитный порядок, третья подрешетка имеет нулевую намагниченность

Абляционные исследования

Проверка поведения в термодинамическом пределе путем изменения размера системы:

  • Фазовый переход первого рода: функция распределения энергии имеет двухпиковую структуру
  • Фазовый переход второго рода: теплоемкость проявляет λ-образную сингулярность, пик растет с размером системы
  • Короткодействующий порядок: теплоемкость проявляет широкий пик, не масштабирующийся с размером системы

Анализ конкретных случаев

Типичное поведение фазовых переходов:

  1. U = -0.7, μ = -0.10: переход первого рода F-NO, сопровождаемый разрывным скачком намагниченности и числа частиц
  2. U = -0.6, μ = -0.12: переход второго рода F-NO, непрерывное изменение параметра порядка
  3. U = -0.2, μ = -0.15: переход второго рода AF-NO, антиферромагнитный порядок двух подрешеток

Экспериментальные находки

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

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

  1. Классическая модель Изинга: Точное решение на треугольной решетке (Houtappel, 1950) обеспечивает эталон для данного исследования
  2. Расширенная модель Хаббарда: Вариационные и численные исследования на гиперкубических решетках обеспечивают методологическую основу
  3. Фрустрированные магнитные системы: Теоретические и экспериментальные исследования антиферромагнитных треугольных решеток обеспечивают контекст для интерпретации результатов
  4. Эксперименты с холодными атомами: Реализация ультрахолодных атомных газов в оптических решетках обеспечивает возможность экспериментальной проверки модели

Выводы и обсуждение

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

  1. Геометрическая фрустрация треугольной решетки значительно влияет на поведение фазовых переходов локализованной фермионной системы
  2. Ферромагнитные и антиферромагнитные взаимодействия приводят к совершенно различным структурам фазовой диаграммы и природе фазовых переходов
  3. Упорядоченная фаза в антиферромагнитном случае обладает уникальной антиферромагнитной структурой двух подрешеток
  4. Все антиферромагнитные фазовые переходы являются непрерывными, в то время как в ферромагнитном случае наблюдается переход от первого рода ко второму роду

Ограничения

  1. Ограничение размера системы: Из-за ограничений вычислительных ресурсов максимальная система составляет 36×36
  2. Приближение атомного предела: Пренебрежение членом прыжка частиц может влиять на количественные предсказания для реальных материалов
  3. Приближение классического спина: Может быть недостаточно точным в области сильных квантовых флуктуаций

Будущие направления

  1. Исследования с включением более дальнодействующих взаимодействий
  2. Применение методов квантового Монте-Карло
  3. Количественное сравнение с реальными материалами
  4. Экспериментальная проверка в системах холодных атомов

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

Преимущества

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

Недостатки

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

Влияние

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

Области применения

  1. Теоретическое моделирование магнитных материалов на треугольной решетке
  2. Теоретическое руководство для экспериментов с холодными атомами в оптических решетках
  3. Исследование фазовых переходов в сильно коррелированных электронных системах
  4. Эталонное тестирование для фрустрированных квантовых многотельных систем

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

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