Энтропия модели клеточной жидкости с взаимодействием Кюри-Вейсса

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

• ID статьи: 2510.22700
• Название: Entropy of the cell fluid model with Curie-Weiss interaction
• Авторы: Р. В. Романик, О. А. Добуш, М. П. Козловский, И. В. Пилюк, М. А. Шпот
• Учреждение: Институт физики конденсированного состояния им. Юхновского, Национальная академия наук Украины
• Классификация: cond-mat.stat-mech (физика конденсированного состояния - статистическая механика)
• Дата подачи: 26 октября 2025 г.
• Ссылка на статью: https://arxiv.org/abs/2510.22700v1

Аннотация

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

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

Исследуемая проблема

Работа посвящена изучению свойств энтропии многозаселённой решёточной газовой модели, в частности модели клеточной жидкости с взаимодействием Кюри-Вейсса. Данная модель позволяет каждой ячейке содержать произвольное число частиц, в отличие от классического решёточного газа, где каждый узел может содержать не более одной частицы.

Значимость проблемы

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

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

• Исследования многозаселённых решёточных газов значительно менее распространены, чем исследования классического решёточного газа
• Отсутствуют систематические аналитические исследования зависимостей энтропии (в частности, соотношений энтропия-плотность и энтропия-химический потенциал) для многозаселённых моделей
• Существующая литература сосредоточена главным образом на квантовых системах (например, модель Хаббарда), тогда как исследования энтропии в классических многозаселённых системах практически отсутствуют

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

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

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

1. Аналитическое выражение энтропии: Впервые выведено точное аналитическое выражение энтропии для многозаселённой модели клеточной жидкости с взаимодействием Кюри-Вейсса, включая одночастичную энтропию S* и энтропию на одну ячейку S*_v
2. Параметрическое представление: Установлены параметрические уравнения энтропии как функции плотности, облегчающие анализ соотношения энтропия-плотность
3. Анализ характеристик фазовых переходов: Систематически исследовано поведение скачков энтропии при бесконечной последовательности фазовых переходов первого рода
4. Обнаружение минимумов при целочисленной плотности: Выявлено, что энтропия проявляет значительные минимумы вблизи целочисленных значений плотности, и обосновано, что это является универсальной характеристикой многозаселённых моделей
5. Сравнение между моделями: Путём сравнения с квантовыми системами, такими как модель Хаббарда и модель Бозе-Хаббарда, выявлены общие характеристики многозаселённых моделей

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

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

Исследование энтропийных свойств точечных частиц в открытой системе, где объём системы V разделён на N_v неперекрывающихся кубических ячеек объёмом v каждая. Цель состоит в получении:

• Одночастичной энтропии: S* = S/(k_B⟨N⟩)
• Энтропии на одну ячейку: S_v = S/(k_B N_v) как функций приведённой температуры T, приведённого химического потенциала μ* и приведённой плотности ρ*

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

Определение взаимодействия

Энергия двухчастичного взаимодействия определяется как: $\Phi_{N_v}(r_i, r_j) = -\frac{J_1}{N_v} + J_2 \sum_{l=1}^{N_v} I_{\Delta_l}(r_i)I_{\Delta_l}(r_j)$

• Первый член: глобальное притяжение Кюри-Вейсса (среднее поле) с параметром интенсивности J_1 > 0
• Второй член: локальное отталкивание между частицами в одной ячейке с параметром интенсивности J_2 > 0
• I_Δl(r): индикаторная функция ячейки Δ_l

Приведённые переменные

Используются стандартные безразмерные переменные:

• Приведённая температура: T* = k_B T/J_1
• Приведённый химический потенциал: μ* = μ/J_1
• Приведённая плотность: ρ* = ⟨N⟩/(N_v)
• Приведённое давление: P* = Pv/J_1
• Параметр: a = J_2/J_1, v* = v/λ³ (λ связана с тепловой длиной волны)

Основные формулы

Одночастичная энтропия

$S^*(T^*, \mu^*) = \left(\frac{3}{2} - \frac{\mu^*}{T^*}\right) - \frac{1}{T^*}\frac{K_1}{K_0} + \frac{K_0\ln K_0}{K_1} + \frac{a}{2T^*}\frac{K_2}{K_1}$

где специальные функции K_j определяются как: $K_j(T^*, \mu^*; \bar{y}_{max}) = \sum_{n=0}^{\infty} \frac{n^j(v^*T^{*3/2})^n}{n!} \exp\left[\left(\bar{y}_{max} + \frac{\mu^*}{T^*}\right)n - \frac{a}{2T^*}n^2\right]$

Энтропия на одну ячейку

$S_v^*(T^*, \mu^*) = \left(\frac{3}{2} - \frac{\mu^*}{T^*}\right)\frac{K_1}{K_0} - \frac{1}{T^*}\frac{K_1^2}{K_0^2} + \ln K_0 + \frac{a}{2T^*}\frac{K_2}{K_0}$

Параметрическое представление

Введена новая функция: $\bar{z}_{max}(T^*, \mu^*) = \bar{y}_{max} + \frac{\mu^*}{T^*}$

Определены: $\tilde{K}_j(T^*; \bar{z}_{max}) = \sum_{n=0}^{\infty} \frac{n^j(v^*T^{*3/2})^n}{n!} \exp\left[\bar{z}_{max}n - \frac{a}{2T^*}n^2\right]$

Параметрическое представление энтропии: $S^*(T^*; \bar{z}_{max}) = \frac{3}{2} - \bar{z}_{max} + \frac{\tilde{K}_0\ln\tilde{K}_0}{\tilde{K}_1} + \frac{a}{2T^*}\frac{\tilde{K}_2}{\tilde{K}_1}$

Плотность: $\rho^*(T^*; \bar{z}_{max}) = \frac{\tilde{K}_1(T^*; \bar{z}_{max})}{\tilde{K}_0(T^*; \bar{z}_{max})}$

Химический потенциал: $\mu^* = T^*\bar{z}_{max} - \frac{\tilde{K}_1(T^*; \bar{z}_{max})}{\tilde{K}_0(T^*; \bar{z}_{max})}$

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

1. Точная разрешимость: Использование среднепольевого характера взаимодействия Кюри-Вейсса для получения точного решения в термодинамическом пределе
2. Параметрическая техника: Введение z̄_max в качестве параметра элегантно преобразует сложные неявные функциональные соотношения в параметрические уравнения, облегчая численные расчёты и анализ
3. Условия экстремума: Определение ȳ_max через условия экстремума функций E_1 и E_2 обеспечивает термодинамическую стабильность
4. Разделение вклада идеального газа: Явное выделение вклада идеального газа (3/2 - μ*/T*) в выражение энтропии облегчает понимание эффектов взаимодействия

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

Выбор параметров

• Параметр взаимодействия: a = J_2/J_1 = 1.2 (отношение интенсивностей отталкивания и притяжения)
• Параметр объёма: v* = 5.0
• Диапазон температур:
  • T* = 0.40 (значительно выше критической температуры)
  • T* = 0.30 (немного выше критической температуры)
  • T* = 0.25 (немного ниже критической температуры)
  • T* = 0.20 (значительно ниже критической температуры)

Данные критических точек

Статья содержит подробные численные значения первых четырёх критических точек (см. Таблица 1):

• 1-я критическая точка: T_c = 0.254567, ρ_c = 0.513896
• 2-я критическая точка: T_c = 0.261881, ρ_c = 1.50557
• 3-я критическая точка: T_c = 0.265254, ρ_c = 2.50303
• 4-я критическая точка: T_c = 0.267242, ρ_c = 3.50191

Данные линий сосуществования

При температурах T* = 0.25 и T* = 0.20 предоставлены данные линий сосуществования первых четырёх фазовых переходов (см. Таблица 2), включая плотность, энтропию, химический потенциал и другие физические величины.

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

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

1. Соотношение энтропия-химический потенциал (Рисунок 1a)

• Надкритическая область (красная кривая T*=0.40 и зелёная кривая T*=0.30): одночастичная энтропия S* монотонно убывает с химическим потенциалом μ*
• Подкритическая область (чёрная кривая T*=0.25 и синяя кривая T*=0.20):
  • Энтропия фазы II ниже энтропии фазы I, энтропия фазы III ниже энтропии фазы II и т.д.
  • При фазовых переходах энтропия скачком переходит от более высокого значения к более низкому
  • За исключением фазы I, зависимость энтропии от химического потенциала в отдельных фазах может быть немонотонной

2. Соотношение энтропия-плотность (Рисунок 1b)

• Надкритическая температура: одночастичная энтропия S* является убывающей функцией плотности ρ*
• Подкритическая температура:
  • В фазе I энтропия остаётся убывающей
  • Ключевое открытие: в последующих фазах энтропия проявляет значительные минимумы вблизи целочисленных значений плотности (ρ* ≈ 1, 2, 3...)
  • При фазовых переходах первого рода энтропия и плотность претерпевают разрывы
  • В фазах высокой плотности (фаза V и выше) энтропия становится отрицательной (артефакт классического подхода)

Конкретные численные значения (T* = 0.20, фаза II):

• При ρ* = 1.12438, S* = 0.996714
• При ρ* = 1.88058, S* = 0.601031
• Минимум существует вблизи ρ* ≈ 1.5

3. Соотношение энтропия на одну ячейку - химический потенциал (Рисунок 1c)

• Высокая температура (T*=0.40): S*_v сначала возрастает, затем убывает, существует максимум
• Вблизи критической температуры (T*=0.30): поведение сложное, вблизи критической точки наблюдается явный изгиб
• Низкая температура:
  • В фазе I S_v возрастает с увеличением μ
  • Фазы II и III имеют наибольшую энтропию на одну ячейку
  • Последующие фазы при фазовых переходах скачком переходят к более низким значениям

4. Соотношение энтропия на одну ячейку - плотность (Рисунок 1d)

• Высокая температура: S*_v сначала возрастает, затем убывает
• Низкая температура:
  • В фазе I возрастает
  • Каждая последующая фаза проявляет немонотонное поведение
  • Существуют минимумы вблизи целочисленных значений плотности (ρ* ≈ 1, 2, 3...)

Важные открытия

Универсальность минимумов при целочисленной плотности

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

1. Трёхмерная модель Хаббарда с фермионами16 (максимальное число заселения = 2):
   • Одноузловая энтропия имеет явное углубление при ρ*=1.0
   • Качественно сходна с результатом данной работы при T*=0.3
2. Одномерная модель Хаббарда с отталкиванием17,18:
   • При высокой температуре одноузловая энтропия имеет единственный максимум
   • При низкой температуре наблюдается значительный минимум при ρ*=1.0
3. Модель Бозе-Хаббарда19 (результаты до числа заполнения n=3.0):
   • Наблюдаются два углубления при n=1.0 и n=2.0
   • Соответствуют минимумам при ρ*≈1.0 и ρ*≈2.0 в данной работе

Единая схема

• Модель Хаббарда (ρ_max=2.0) → минимум при ρ=1.0
• Модель Бозе-Хаббарда (результаты до ρ*≤3.0) → минимумы при ρ*=1.0 и 2.0
• Модель Кюри-Вейсса данной работы (без ограничений на заселение) → бесконечная последовательность минимумов при ρ*≈1, 2, 3...

Явление отрицательной энтропии

При высокой плотности (фазы V и выше) или малых значениях v* энтропия становится отрицательной. Это артефакт классического подхода, аналогичный наблюдению Стиллинджера в модели с гауссовым ядром. Увеличение v* смещает это поведение в область более высоких плотностей, но качественная картина остаётся неизменной.

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

Многозаселённые решёточные газовые модели

1. Теория Ли-Янга1: основы статистической теории классического решёточного газа
2. Жидкости с мягким взаимодействием3-8:

• Системы с гауссовым ядром3
• Проникающие сферы4
• Двузаселённый решёточный газ7,8

3. Теория среднего поля9,10: модель Кюри-Вейсса с взаимодействием, на которой основана данная работа
4. Модели многослойной адсорбции11: решёточный газ, допускающий многозаселение
5. Модели кластерных кристаллов12,13: обобщённые модели степени 4

Квантовые многозаселённые модели

1. Модель Хаббарда16-18:

• Исследование трёхмерной модели Хаббарда с фермионами методом динамической теории среднего поля
• Исследование одномерной модели Хаббарда с отталкиванием методом функционала плотности

2. Модель Бозе-Хаббарда19: исследование методом квантового Монте-Карло

Сравнение с жидкостью Ван-дер-Ваальса

В Приложении A подробно проанализирована энтропия жидкости Ван-дер-Ваальса, обнаружено:

• Только переход из фазы I в фазу II может быть качественно сравним с жидкостью Ван-дер-Ваальса
• Это указывает на то, что данный фазовый переход первого рода относится к типу газ-жидкость

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

1. Аналитическая точность: предоставляет полностью аналитическое выражение энтропии, не зависящее от численных приближений
2. Системная полнота: одновременно исследует одночастичную энтропию и энтропию на одну ячейку, анализирует на нескольких термодинамических плоскостях
3. Обоснование универсальности: путём сравнения между моделями выявляет общие свойства многозаселённых систем
4. Бесконечная последовательность фаз: модель обладает бесконечной последовательностью фазовых переходов первого рода, что богаче, чем модели с конечным заселением

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

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

1. Успешное выведение аналитических выражений: получены точные аналитические выражения энтропии для многозаселённой модели клеточной жидкости с взаимодействием Кюри-Вейсса, включая одночастичную энтропию и энтропию на одну ячейку как функции температуры и химического потенциала
2. Энтропийные характеристики фазовых переходов: при температурах ниже критической энтропия проявляет разрывы при фазовых переходах первого рода, полностью соответствующие бесконечной каскадной области сосуществования
3. Минимумы при целочисленной плотности: энтропия проявляет значительные минимумы вблизи целочисленных значений плотности, непосредственно связанные с дискретным характером заселения ячеек, являясь универсальной термодинамической характеристикой многозаселённых решёточных систем
4. Универсальность между моделями: путём качественного сравнения с решёточными моделями различного микроскопического происхождения (Хаббард, Бозе-Хаббард) подтверждено, что все системы, допускающие более чем одночастичное заселение, развивают минимумы энтропии при целочисленной плотности
5. Идентификация типа фазового перехода: качественное сходство перехода из фазы I в фазу II с жидкостью Ван-дер-Ваальса указывает на то, что этот переход относится к типу газ-жидкость

Ограничения

1. Ограничения классического подхода:

• При высокой плотности или малых v* энтропия становится отрицательной
• Это внутренний артефакт классического статистического подхода
• Требуется квантовая коррекция для полной точности

2. Приближение среднего поля:

• Взаимодействие Кюри-Вейсса имеет среднепольевой характер
• Игнорируются короткодействующие флуктуации и пространственные корреляции
• Критические показатели имеют среднепольевые значения

3. Параметрическая зависимость:

• Результаты чувствительны к параметрам v* и a
• При малых v* число фаз с положительной энтропией ограничено

4. Отсутствие экспериментальной проверки:

• Как теоретическое исследование, отсутствует прямое сравнение с экспериментальными данными
• Экспериментальная реализация многозаселённых классических систем редка

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

1. Проверка на других многозаселённых моделях:

• Двузаселённый решёточный газ7,8 должен иметь углубление энтропии при ρ*=1.0
• Энтропийные характеристики модели многослойной адсорбции11
• Минимумы энтропии при целочисленной плотности в модели кластерных кристаллов12,13

2. Выход за пределы среднего поля:

• Рассмотрение коррекций короткодействующих флуктуаций
• Исследование эффектов конечного размера
• Анализ методом ренормализационной группы

3. Квантовое расширение:

• Разработка квантовой версии модели клеточной жидкости
• Количественное сравнение с моделью Бозе-Хаббарда

4. Практические приложения:

• Мягкие материалы (коллоиды, полимеры)
• Системы с кластерообразующим взаимодействием
• Теоретическое руководство для экспериментов с холодными атомами в оптических решётках

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

Достоинства

1. Методологическая инновативность

• Полностью аналитическое решение: получение точного аналитического выражения энтропии в многозаселённых моделях чрезвычайно редко и представляет важный теоретический эталон для этой области
• Параметрическая техника: элегантное разрешение сложных неявных функциональных соотношений путём введения параметра z̄_max
• Системность: одновременное исследование одночастичной энтропии и энтропии на одну ячейку на нескольких термодинамических плоскостях обеспечивает полную картину

2. Научные открытия

• Универсальность минимумов при целочисленной плотности: это наиболее важное открытие работы, выявляющее фундаментальную характеристику многозаселённых моделей
• Межмодельные связи: связь классических статистико-механических моделей с квантовыми многотельными системами демонстрирует глубокие физические инсайты
• Бесконечная последовательность фаз: систематическое исследование модели с бесконечной последовательностью фазовых переходов первого рода теоретически весьма сложно

3. Качество представления

• Ясные графики: четыре панели на Рисунке 1 систематически представляют различные аспекты энтропии с чёткой цветовой кодировкой
• Подробные численные данные: Таблицы 1 и 2 содержат точные значения критических точек и линий сосуществования, удобные для цитирования в последующих исследованиях
• Логическая строгость: от определения модели, аналитического вывода до численных результатов и физического обсуждения структура ясна

4. Теоретическая ценность

• Предоставляет теоретические инструменты для физики мягких материалов
• Предлагает новую перспективу для понимания термодинамики систем с кластерообразованием
• Устанавливает единую рамку понимания классических и квантовых многозаселённых систем

Недостатки

1. Недостаточная глубина физической интерпретации

• Микроскопический механизм минимумов энтропии при целочисленной плотности объяснён недостаточно
• Почему многозаселение приводит к этому универсальному поведению? Отсутствует более фундаментальное статистико-механическое обоснование
• Физический образ различных фаз (фаза I, II, III и т.д.) описан недостаточно

2. Ограничения численного диапазона

• Представлены результаты только для первых 4-5 фаз
• Поведение более высоких фаз (фаза VI и выше) обсуждено недостаточно
• Выбор v*=5.0 не имеет достаточного обоснования

3. Связь с экспериментом

• Полностью теоретическая работа без сравнения с экспериментом
• Не обсуждается, как наблюдать предсказанные энтропийные характеристики в реальных системах
• Соответствие параметров a=1.2, v*=5.0 реальным материалам неясно

4. Обсуждение ограничений модели

• Недостаточно обсуждены ограничения, связанные с дальнодействующим характером взаимодействия Кюри-Вейсса
• Почему рассматривается трёхмерный случай? Почему не одномерный или двумерный?
• Область применимости приближения точечных частиц не определена

5. Математическая строгость

• Условия сходимости специальных функций K_j не обсуждены подробно
• Математическая строгость термодинамического предела рассмотрена кратко
• Единственность и глобальность условия экстремума (11) не обоснованы достаточно

Оценка влияния

Вклад в область

• Теоретический эталон: предоставляет важный аналитический эталон для исследований многозаселённых моделей
• Принцип универсальности: минимумы энтропии при целочисленной плотности могут стать стандартной характеристикой для идентификации и классификации многозаселённых систем
• Междисциплинарный мост: связывает классическую статистическую механику, физику мягких материалов и теорию квантовых многотельных систем

Практическая ценность

• Средняя: как фундаментальное теоретическое исследование, прямое применение ограничено
• Может направлять:
  • Понимание фазового поведения коллоидных систем
  • Теоретические предсказания для экспериментов с холодными атомами
  • Проектирование новых мягких материалов

Воспроизводимость

• Высокая: предоставлены полные аналитические формулы и параметры
• Численные расчёты могут быть непосредственно повторены
• Рекомендуется авторам опубликовать вычислительный код для облегчения применения

Применимые сценарии

Теоретические исследования

1. Эталонное тестирование многозаселённых решёточных моделей
2. Учебный пример теории среднего поля
3. Объект исследования теории фазовых переходов

Потенциальные приложения

1. Мягкие материалы:

• Проникающие коллоидные частицы
• Растворы звездообразных полимеров
• Дендримерные макромолекулы

2. Системы холодных атомов:

• Бозе-газ в оптических решётках
• Многокомпонентные атомные смеси

3. Явления адсорбции:

• Многослойная газовая адсорбция
• Образование поверхностных кластеров

4. Развитие теоретических методов:

• Тестирование методов, выходящих за пределы среднего поля
• Разработка новых численных алгоритмов

Ключевые ссылки

1. Lee & Yang (1952): основы статистической теории решёточного газа и модели Изинга
2. Stanley (1971): классический учебник по фазовым переходам и критическим явлениям
3. Stillinger (1976): исследование фазовых переходов в системах с гауссовым ядром
4. Kozlovskii et al. (2020, 2022): предыдущие работы авторов по модели Кюри-Вейсса с клеточной жидкостью
5. De Leo et al. (2011): исследование термодинамики трёхмерной модели Хаббарда
6. Pollet et al. (2008): исследование модели Бозе-Хаббарда методом квантового Монте-Карло
7. Hill (1956): классический учебник статистической механики, справочник по энтропии идеального газа

Общая оценка

Это высококачественная работа по теоретической статистической механике, вносящая важный вклад в исследование энтропии многозаселённых решёточных моделей. Основная ценность работы заключается в: (1) предоставлении полностью аналитического выражения энтропии; (2) обнаружении и обосновании универсальности минимумов энтропии при целочисленной плотности; (3) установлении единой рамки понимания различных типов многозаселённых систем. Несмотря на некоторые ограничения (главным образом приближение среднего поля и отсутствие экспериментальной проверки), как фундаментальное теоретическое исследование работа демонстрирует глубокие физические инсайты и строгую математическую обработку. Данная работа станет важным справочником для исследований многозаселённых моделей и может стимулировать дальнейшие исследования в области физики мягких материалов и холодных атомов.

Рекомендуемый рейтинг: ⭐⭐⭐⭐☆ (4/5)

• Теоретическая инновативность: ⭐⭐⭐⭐⭐
• Экспериментальная релевантность: ⭐⭐⭐☆☆
• Методологическая строгость: ⭐⭐⭐⭐☆
• Ясность представления: ⭐⭐⭐⭐⭐
• Потенциал влияния: ⭐⭐⭐⭐☆