Термодинамика сильных взаимодействий

Квантовая хромодинамика при конечной температуре и плотности


Сильные взаимодействия при конечной температуре и/или барионной плотности описываются квантовой хромодинамикой (КХД) в условиях термодинамического равновесия. Центральными характеристиками являются:

  • Температура T
  • Химический потенциал барионного числа μB
  • Объём системы V
  • Энергетическая плотность ε
  • Давление p
  • Плотность энтропии s

КХД-плазма, состоящая из кварков и глюонов, описывается как термодинамическая система, свойства которой можно исследовать методами статистической физики. Важной целью является определение уравнения состояния (УоС) и фазовых переходов, особенно перехода от обычной ядерной материи к кварк-глюонной плазме (КГП).


Канонический формализм и функция распределения

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

Z(T, μ, V) = Tr[e−( − μ)/T]

где — гамильтониан системы, — оператор числа частиц (или барионного числа), μ — химический потенциал. Все макроскопические термодинамические величины выводятся из логарифма Z:

$$ p = \frac{T}{V} \ln Z, \quad \varepsilon = -\frac{1}{V} \left( \frac{\partial \ln Z}{\partial \beta} \right) + \mu n, \quad s = \frac{\varepsilon + p - \mu n}{T} $$

где β = 1/T, а n — плотность барионного числа.


Решеточная КХД при конечной температуре

Вычисление функции распределения для КХД возможно только численно на решётке (решёточная КХД). Температура в решёточной теории определяется как:

$$ T = \frac{1}{N_t a} $$

где a — шаг решётки, Nt — число временных слоёв. Для расчётов используют периодические граничные условия по времени (для бозонов) и анти-периодические (для фермионов), что позволяет изучать систему при заданной температуре.

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

  • Трассировка энергии-импульса: Θμμ = ε − 3p
  • Энтропийная плотность: $s = \frac{\varepsilon + p}{T}$
  • Сжимаемость: $\chi_B = \frac{\partial n_B}{\partial \mu_B}$

Уравнение состояния и фазовый переход

Из расчётов следует, что при температуре около Tc ≈ 150 − 160 MeV происходит переход из конфинированной фазы (где активны адроны) в деконфинированную фазу кварков и глюонов. Характер этого перехода зависит от числа кварков и их масс:

  • При реальных массах — перекрестный переход (crossover)
  • При массе лёгких кварков mq → 0второго рода переход
  • При высоких μB — возможен фазовый переход первого рода и критическая точка

Температурная зависимость давления и энергии в КГП на решётке демонстрирует быстрое нарастание при T > Tc, с последующим насыщением значений, близких к ультралучевой идеальной газовой плазме:

$$ \varepsilon_{\text{ideal}} = \left( N_g + \frac{7}{8} N_f N_c \right) \frac{\pi^2}{30} T^4 $$

где Ng = 2(Nc2 − 1) — число степеней свободы глюонов, Nf — число ароматов, Nc = 3.


Химический потенциал и фазовая диаграмма

Введение ненулевого барионного химического потенциала μB важно для описания:

  • Ядерной материи в ядрах и нейтронных звёздах
  • Коллизий тяжёлых ионов на низких и средних энергиях

Однако прямая решёточная симуляция при μB > 0 затруднена из-за сигнатурной проблемы (комплексность детерминанта фермионного оператора). Используют приближения:

  • Разложение по малому μB/T
  • Метод аналитического продолжения
  • Имитация с помощью моделей (PNJL, HRG и др.)

Фазовая диаграмма в плоскости (T, μB) содержит:

  • Линию перекрестного перехода при μB ≈ 0
  • Потенциальную критическую точку на μB ≈ 200 − 500 MeV
  • Линию первого рода при больших μB

Адронный газ и модель резонансов

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

$$ p(T,\mu) = \sum_i \pm \frac{d_i}{6\pi^2} \int_0^\infty \frac{p^4 dp}{\sqrt{p^2 + m_i^2}} \left( \exp\left[ \frac{\sqrt{p^2 + m_i^2} \mp \mu_i}{T} \right] \pm 1 \right)^{-1} $$

где знак «+» для фермионов, «−» для бозонов, di — степень вырождения, μi — химпотенциал, mi — масса. Модель HRG успешно воспроизводит данные решёточной КХД при T < Tc, включая флуктуации зарядов и химические потенциальные производные давления.


Флуктуации и кумулянты

Особое значение в изучении термодинамики КХД имеют флуктуации консервативных зарядов: барионного, странного, электрического. Флуктуации характеризуются кумулянтами:

$$ \chi_n^X = \frac{\partial^n (p/T^4)}{\partial (\mu_X/T)^n} $$

где X = B, Q, S. Кумулянты 2-го порядка соответствуют дисперсиям, 4-го — эксцессам, и т.д. В экспериментах по тяжёлым ионам сравнение с теоретическими кумулянтами позволяет:

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

Вязкость и неравновесные эффекты

Помимо равновесной термодинамики, важна оценка транспортных коэффициентов, особенно отношения сдвиговой вязкости к энтропии η/s. Считается, что кварк-глюонная плазма близка к идеальной жидкости:

$$ \frac{\eta}{s} \gtrsim \frac{1}{4\pi} $$

— нижний предел, выведенный в теориях с голографической двойственностью. Современные результаты из феноменологии коллизий тяжёлых ионов (ALICE, RHIC) подтверждают малую вязкость КГП.


Связь с экспериментами и астрофизикой

Изучение термодинамики сильных взаимодействий находит подтверждение и применение в:

  • Анализе коллизий тяжёлых ионов (RHIC, LHC, FAIR, NICA)
  • Моделировании фазового перехода в ранней Вселенной
  • Описании вещества в нейтронных звёздах и гравитационно-волновых событиях

Данные по термодинамике, полученные на решётке, являются необходимыми входами в гидродинамическое моделирование, кинетические расчёты и астрофизические симуляции.


Таким образом, термодинамика КХД представляет собой фундаментальную область, где квантовая теория поля тесно переплетается с методами статистической физики и экспериментальной ядерной и астрофизики.