2025-11-19T17:58:14.286609

Charge Susceptibility and Kubo Response in Hatsugai-Kohmoto-related Models

Ma, Zhao, Huang et al.
We study in depth the charge susceptibility for the band Hatsugai-Kohmoto (HK) and orbital (OHK) models. As either of these models describes a Mott insulator, the charge susceptibility takes on the form of a modified density response function with lower and upper Hubbard bands, thereby giving rise to a multi-pole structure. The particle-hole continuum consists of hot spots along the $ω$ vs $q$ axis arising from inter-band transitions. Such transitions, which are strongly suppressed in non-interacting systems, obtain here because of the non-rigidity of the Hubbard bands. This modified density response function gives rise to a plasmon dispersion that is inversely dependent on the momentum, resulting in an additional contribution to the conventional f-sum rule. This extra contribution originates from a long-range diamagnetic contribution to the current. This results in a non-commutativity of the long-wavelength ($q\rightarrow 0$) and thermodynamic ($L\rightarrow\infty$) limits. When the correct limits are taken, we find that the Kubo response computed with either open or periodic boundary conditions yields identical results that are consistent with the continuity equation contrary to recent claims. We also show that the long wavelength pathology of the current noted previously also plagues the Anderson impurity model interpretation of dynamical mean-field theory (DMFT).
academic

Зарядовая восприимчивость и отклик Кубо в моделях, связанных с Хацугаи-Кохмото

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

  • ID статьи: 2409.07522
  • Название: Charge Susceptibility and Kubo Response in Hatsugai-Kohmoto-related Models
  • Авторы: Yuhao Ma, Jinchao Zhao, Edwin W. Huang, Dhruv Kush, Barry Bradlyn, Philip W. Phillips
  • Классификация: cond-mat.str-el (сильнокоррелированные электронные системы)
  • Дата публикации: сентябрь 2024 г. (препринт arXiv)
  • Ссылка на статью: https://arxiv.org/abs/2409.07522

Аннотация

В данной работе проводится глубокое исследование зарядовой восприимчивости полосовой модели Хацугаи-Кохмото (ХК) и орбитальной модели ХК (ОХК). Поскольку обе модели описывают изолятор Мотта, их зарядовая восприимчивость проявляется в виде модифицированной функции плотностного отклика с нижней и верхней полосами Хаббарда, что приводит к мультипольной структуре. Спектр частица-дырка на плоскости ω против q содержит горячие точки, возникающие из межзонных переходов. Эти переходы сильно подавлены в невзаимодействующей системе, но здесь реализуются благодаря нежёсткости полос Хаббарда. Эта модифицированная функция плотностного отклика приводит к дисперсии плазмонов, обратно пропорциональной импульсу, что вызывает дополнительный вклад в традиционное f-сумм-правило. Этот дополнительный вклад возникает из-за дальнодействующего диамагнитного вклада тока, приводящего к некоммутативности пределов малых длин волн (q→0) и термодинамического предела (L→∞).

Научный контекст и мотивация

Проблемный фон

  1. Вызовы в сильнокоррелированных электронных системах: В сильнокоррелированных электронных системах (таких как высокотемпературные сверхпроводники на основе меди) традиционная теория жидкости Ферми не работает, требуется новая теоретическая база для понимания свойств зарядового отклика.
  2. Особенности изолятора Мотта: Изолятор Мотта вследствие сильных эффектов электронной корреляции демонстрирует расщепление зонной структуры на верхнюю и нижнюю полосы Хаббарда, что резко контрастирует с предсказаниями традиционной теории зон.
  3. Аномальные экспериментальные наблюдения: Недавние эксперименты по спектроскопии потерь энергии электронов с разрешением по импульсу (MEELS) в медных оксидах выявили:
    • Плоский спектр, независимый от температуры и импульса на средних энергетических масштабах
    • χ''∝q² при низких энергиях и затухание как q²/ω² при высоких энергиях
    • Плазмонные моды, независимые от плотности

Научная мотивация

  1. Теоретический пробел: Несмотря на то, что модель ХК широко изучается как точно решаемая платформа для физики Мотта, её функция плотностного отклика никогда ранее не подвергалась глубокому анализу.
  2. Фундаментальные теоретические вопросы: Необходимо понять в сильнокоррелированных системах:
    • Как изменяется спектр частица-дырка при наличии нежёстких полос Хаббарда
    • Остаётся ли в силе компенсация импульсного множителя, приводящая к плазменной частоте
    • Требуется ли модификация f-сумм-правила

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

  1. Первый полный расчёт зарядовой восприимчивости модели ХК: Предоставлены точные аналитические выражения для функции плотностного отклика полосовой модели ХК и орбитальной модели ХК.
  2. Обнаружение аномальной дисперсии плазмонов: Доказано, что плазменная частота в модели ХК ωₚ∼1/q, что принципиально отличается от поведения в традиционных металлах.
  3. Раскрытие модифицированного f-сумм-правила: Обнаружено, что благодаря дальнодействующему диамагнитному вкладу тока f-сумм-правило получает дополнительный член, нарушающий традиционный закон масштабирования q².
  4. Разрешение спора об определении оператора тока: Строгий вывод калибровочно-инвариантного оператора тока из принципа минимального взаимодействия, доказательство справедливости уравнения непрерывности.
  5. Установление связи с DMFT: Доказано, что однопримесная модель Андерсона в рамках DMFT также демонстрирует аналогичные патологии в длинноволновом пределе.

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

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

Исследование зарядовой восприимчивости χ(q,ω) в моделях, связанных с ХК, где:

  • Входные данные: импульс q и частота ω
  • Выходные данные: функция корреляции плотность-плотность
  • Ограничения: удовлетворение законам сохранения заряда и уравнению непрерывности

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

1. Полосовая модель ХК (bHK)

Гамильтониан имеет вид:

H_bHK = Σ_kσ ξ_k n_kσ + Σ_k U n_k↑ n_k↓

где ξ_k = ε_k - μ — кинетический член, U — независимое от импульса отталкивание плотность-плотность.

2. Орбитальная модель ХК (OHK)

H_OHK = Σ_k,α,α',σ g_αα'(k) c†_kασ c_kα'σ - μΣ_k,α n_kα,σ + Σ_k,α,α' U_αα' n_kα↑ n_kα'↓

3. Расчёт функции плотностного отклика

Использование свойства факторизации по импульсу в модели ХК:

χ⁰_HK(q,ω) = (1/N) Σ_kσ Σ_i,j=L,U [w^i_k+q w^j_k (f(ξ^j_k) - f(ξ^i_k+q))/(ω + i0⁺ - ξ^i_k+q + ξ^j_k)]

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

1. Точная обработка многозонной структуры

  • Отображение двухполюсной структуры ХК на эффективную многозонную систему
  • Различие между прямыми членами (L,L) и (U,U) и перекрёстными членами (L,U) и (U,L)
  • Перекрёстные члены приводят к вкладу порядка q⁰, характерному для систем Мотта

2. Тщательный анализ длинноволнового предела

Обнаружена критическая некоммутативность:

lim_q→0 lim_L→∞ ≠ lim_L→∞ lim_q→0

Это требует сначала взять предел q→0 в конечной системе, а затем термодинамический предел.

3. Построение калибровочно-инвариантного оператора тока

Через принцип минимального взаимодействия:

U_αβ(R₁+R₃-R₂-R₄) → U^A_αβ(R₁+R₃-R₂-R₄) = U_αβ(R₁+R₃-R₂-R₄)e^(-i∫A·dx)

Вывод оператора плотности тока, удовлетворяющего уравнению непрерывности.

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

Теоретическая вычислительная база

  1. Условия приближения:
    • Длинноволновый высокочастотный предел: q·P_i/(m(ω+U)) << 1
    • Предел низких температур: ∂f/∂ξ_k ≈ -δ(ξ_k)
    • Параболическая дисперсия: ξ_k = |k|²/2m - μ
  2. Обработка в RPA:
    • Точное решение взаимодействия ХК
    • Кулоновское взаимодействие в приближении RPA
    • Диэлектрическая функция: ε_RPA(q,ω) = 1 - V_q χ⁰_HK(q,ω)

Численные параметры

  • Дисперсия плотной связи на квадратной решётке
  • t = 0,1 эВ, t_p = 0,03 эВ, U = 0,14 эВ
  • Систематическое исследование размера системы и граничных условий

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

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

1. Аномальное поведение дисперсии плазмонов

Обнаружено, что плазменная частота подчиняется:

ω_p ∼ 1/q

Это резко контрастирует с постоянной ω_p в традиционных металлах.

2. Модифицированное f-сумм-правило

Традиционное f-сумм-правило:

∫ dω ω Im χ(q,ω) = -π n_e q²/m

В модели ХК получает дополнительный вклад:

∫ dω ω Im χ_HK(q,ω) = (4πU/N) Σ_k (⟨n_k↓n_k↑⟩ - ⟨n_k↓⟩⟨n_k↑⟩)

3. Структура горячих точек в спектре частица-дырка

  • Значительные горячие точки появляются при ω = U, соответствуя межзонным переходам
  • Контрастирует с подавленными межзонными переходами в невзаимодействующей системе
  • Интенсивность горячих точек пропорциональна объёму области одиночной занятости

Абляционные эксперименты

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

  • Прямые члены: χ_dir ∝ q²/ω², подобно свободным электронам
  • Перекрёстные члены: χ_cro ∝ 1/ω², обеспечивают доминирующий вклад

Проверка согласованности граничных условий

Доказано, что при открытых и периодических граничных условиях:

  • Форма оператора тока одинакова
  • Уравнение непрерывности выполняется везде
  • Поведение в длинноволновом пределе идентично

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

История развития модели ХК

  1. Исходная модель ХК (Hatsugai & Kohmoto, 1992): Предложена точно решаемая сильнокоррелированная модель
  2. Орбитальное обобщение (последние годы): Развита многоорбитальная версия, улучшающая приближение модели Хаббарда
  3. Исследование топологических свойств: Применение к топологическим изоляторам Мотта и т.д.

Исследование функций плотностного отклика

  1. Численные исследования DCA (Dong et al.): Численные расчёты на мнимых частотах, требующие аналитического продолжения
  2. Эксперименты MEELS: Недавние наблюдения в медных оксидах
  3. Традиционная теория: Ограничения теории RPA и жидкости Ферми

Спор об операторе тока

Недавно Guerci и др. подвергли сомнению определение оператора тока в модели ХК; данная работа даёт полное уточнение.

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

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

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

Ограничения

  1. Приближение RPA: Обработка кулоновского взаимодействия остаётся на уровне RPA
  2. Упрощение модели: Модель ХК, хотя точно решаема, всё ещё упрощена относительно реальных материалов
  3. Ограничение размерности: Строгие результаты справедливы при d=∞, конечные размерности требуют коррекции

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

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

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

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

  1. Теоретическая строгость: Предоставлена первая полная теоретическая база для зарядового отклика модели ХК
  2. Математическая точность: Все расчёты точны без неконтролируемых приближений
  3. Глубокие физические идеи: Раскрыты новые явления и механизмы в физике Мотта
  4. Разрешение споров: Решена важная полемика об определении оператора тока
  5. Ясное изложение: Логичная структура и чёткое выражение

Недостатки

  1. Связь с экспериментом: Отсутствует прямое сравнение с конкретными экспериментами
  2. Численная проверка: Некоторые аналитические результаты лишены независимой численной верификации
  3. Применимость: Применимость к реальным материалам требует дальнейшей проверки

Влияние

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

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

  1. Сильнокоррелированные электронные системы: Особенно изоляторы Мотта и связанные материалы
  2. Исследование теоретических моделей: Как важный инструмент для изучения эффектов сильной корреляции
  3. Верификация численных методов: Как эталонный тест для методов типа DMFT

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

Данная работа цитирует 48 важных источников, охватывающих модель ХК, сильнокоррелированные электронные системы, теорию DMFT и другие смежные области, обеспечивая прочную теоретическую базу для исследования.


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