Сразу после Большого взрыва во Вселенной наступила короткая фаза, когда кварки и глюоны свободно перемещались, ещё не объединившись в адроны, такие как протоны, нейтроны и мезоны. Это состояние, называемое кварк-глюонной плазмой, просуществовало недолго, пока температура не опустилась примерно до 20 триллионов Кельвин, после чего произошла адронизация.
Группа исследователей из Италии представила новые расчёты уравнения состояния плазмы, которые показывают, как велика была роль сильного взаимодействия до формирования адронов. Их работа опубликована в журнале Physical Review Letters.
Уравнение состояния квантовой хромодинамики (КХД)
Уравнение состояния КХД описывает коллективное поведение частиц, испытывающих сильное взаимодействие — газ сильно взаимодействующих частиц в равновесии, с неизменным количеством частиц и чистой энергией. Это аналогично хорошо известному простому уравнению состояния атомов в газе, $PV=nRT$, но не может быть так просто обобщено.
Однако, подобно классическому газу, совокупность частиц КХД в равновесии имеет температуру, давление, плотность энергии и энтропию, и может претерпевать фазовые переходы.
Первые важные фазовые переходы:
* Прямое КХД не объясняет своё уравнение состояния. Теория возмущений, основное оружие в арсенале физиков, где члены диаграммы Фейнмана вычисляются с помощью степеней константы связи, не работает для сильного взаимодействия, как для электромагнитного взаимодействия и квантовой электродинамики (КЭД), где константой связи является малая постоянная тонкой структуры, приблизительно равная 1/137.
* В КХД константа связи не мала. Теория неабелева (в отличие от фотона КЭД, который не несёт электромагнитного заряда, переносчики силы КХД, глюоны, действительно имеют цветовой заряд (на самом деле два — цвет и антицвет).
* Более того, константа связи КХД изменяется в зависимости от энергии взаимодействия — на малых расстояниях взаимодействие мало, но на больших расстояниях сила огромна, что проявляется в асимптотической свободе.
Поэтому физики обратились к решётчатой КХД для вычисления уравнения состояния.
Решётчатая КХД
В решётчатой КХД пространство-время делится на дискретные точки на четырёхмерном кубе, а свойства пространства-времени КХД вычисляются точечно непертурбативным образом, без использования диаграмм Фейнмана.
В конце концов, расстояние между точками пространства-времени становится всё меньше и меньше, но для выполнения большинства расчётов решёточной КХД всё равно требуются суперкомпьютеры.
Используя решёточную КХД, исследователи из Университета Милан-Бикокка и Национального института ядерной физики (INFN) в Италии попытались определить уравнение состояния КХД от температуры 3 ГэВ до электрослабого перехода.
Они сосредоточились на сильно взаимодействующей системе безмассовых частиц кварков, где большая часть их массы связана с окружающими их глюонными полями и составляет менее 500 МэВ/с$^2$ при таком температурном масштабе (так что примерно равна нулю относительно энергии плазмы).
Исследователи говорят, что «вычислительная стратегия полностью нова, и мы фокусируемся на теории с тремя ароматами безмассовых кварков», разработанной в 2022 году тремя соавторами. По сути, новая стратегия использует метод Монте-Карло для изучения решёточной КХД от низких до очень высоких температур на основе первых принципов, получая численные результаты путём случайного отбора проб.
После вычислений группа получила уравнение состояния для плотности энтропии кварк-глюонной плазмы от температур 3 ГэВ до 165 ГэВ для трёх ароматов кварков, вплоть до температуры электрослабого перехода, выразив его численно в виде полинома седьмой степени (суммы степеней) константы связи сильного взаимодействия, которая сама является функцией температуры.
Они взяли предел численно, чтобы фактически уменьшить расстояние решётки до нуля, так что их результаты применимы к реальному миру. «Решётчатые артефакты оказываются довольно мягкими», — заключают они. Это значительное улучшение по сравнению с предыдущими симуляциями кварк-глюонной плазмы, которые ограничивались температурами ниже 1 ГэВ.
Из плотности энтропии давление и плотность энергии были рассчитаны с помощью стандартных термодинамических уравнений. Они также определили, что давления, которые они рассчитали, не могут быть точно описаны моделью слабо взаимодействующих кварков и глюонов, что указывает на то, что сильное взаимодействие оказало влияние на раннюю Вселенную раньше, чем считалось ранее.
Чтобы пойти дальше, им нужны более быстрые компьютеры или больше компьютерного времени: «представленные здесь численные результаты действительно могут быть систематически улучшены в будущем за счёт вложения дополнительных вычислительных ресурсов».