Nonequilibrium Phenomenology of Identified Particle Spectra in Heavy-Ion Collisions at LHC Energies

Resumo Técnico: Fenomenologia de Não-Equilíbrio de Espectros de Partículas Identificadas em Colisões de Íons Pesados em Energias do LHC

Definição do Problema
O artigo aborda uma discrepância persistente na descrição teórica de espectros de partículas identificadas em colisões de íons pesados (HIC) ultra-relativísticas em energias do CERN LHC. Enquanto modelos estatísticos descrevem com sucesso as razões de rendimento de partículas via um cenário de congelamento químico, eles falham em reproduzir com precisão os espectros de momento transversal ($p_T$) de píons e kaons em baixos momentos. Os modelos de funções de distribuição térmica padrão e códigos hidrodinâmicos convencionais subestimam sistematicamente o rendimento de píons na região de baixo $p_T$. Explicações anteriores atribuíram esse aumento a decaimentos de ressonâncias e ao fluxo radial, mas dados recentes de alta precisão da Colaboração ALICE indicam que esses mecanismos são insuficientes para descrever a sobrepopulação observada de estados de baixa energia.

Metodologia
Os autores empregam o método de Zubarev do Operador Estatístico de Não-Equilíbrio (NSO) para desenvolver uma abordagem sistemática para descrever a produção de píons fora do equilíbrio. O núcleo da metodologia envolve:

1. Estrutura Teórica: Utilizando o formalismo NSO para derivar funções de distribuição onde flutuações de longa duração (especificamente números de píons e kaons, $N_\pi$ e $N_K$) são tratadas como observáveis relevantes. Isso leva à introdução de potenciais químicos efetivos de não-equilíbrio ($\mu_\pi, \mu_K$) na função de distribuição, distintos dos parâmetros de equilíbrio.
2. Abordagens de Modelagem: O estudo testa quatro modelos distintos (A–D) para descrever o congelamento de partículas:
   • Modelo A: Um modelo de onda de choque (blast-wave) padrão sem decaimentos de ressonância ou interações de estado final.
   • Modelo B: Um modelo de onda de choque incorporando potenciais químicos efetivos de não-equilíbrio para píons e kaons dentro da distribuição de Bose.
   • Modelos C & D: Uma abordagem híbrida acoplando um gerador de partículas térmicas baseado em blast-wave ao modelo de transporte de hádrons SMASH. O Modelo C assume o equilíbrio químico total e o Modelo D inclui um potencial químico de píons efetivo para contabilizar estados sobrepopulados.
3. Análise de Dados: Os modelos são ajustados aos espectros de momento transversal de alta precisão de partículas identificadas ($\pi^\pm, K^\pm, p/\bar{p}$) medidos pela ALICE em colisões Pb+Pb centrais em $\sqrt{s_{NN}} = 2,76$ TeV.
4. Inferência Estatística: Métodos de inferência bayesiana, utilizando emuladores de processo gaussiano e amostragem de Cadeia de Markov Monte Carlo (MCMC), são aplicados para extrair os parâmetros termodinâmicos mais prováveis (temperatura, velocidade de fluxo, potenciais químicos) e avaliar as probabilidades dos modelos via fatores de Bayes.

Resultos Principais

• Modelo A (Padrão): Replica achados anteriores, resultando em uma temperatura de congelamento de $\sim 97,7$ MeV e velocidade de superfície de $\sim 0,88c$. Ele desvia-se significativamente dos dados da ALICE em baixo $p_T$.
• Modelo B (Potencial Químico de Não-Equilíbrio): A introdução de potenciais químicos efetivos ($\mu_\pi \approx 126$ MeV, $\mu_K \approx 403$ MeV) permite que o modelo descreva o aumento de baixo $p_T$ sem decaimentos de ressonância. Este modelo resulta em uma temperatura mais alta ($\sim 118$ MeV) e fornece um ajuste significativamente melhor (melhoria de $\chi^2$ por um fator de 6,7) sobre uma ampla faixa de $p_T$, incluindo regiões fora dos intervalos de ajuste.
• Modelos C & D (Acoplados ao Transporte): Quando decaimentos de ressonância e interações de estado final são explicitamente simulados via SMASH, a temperatura resultante é de $\sim 156$ MeV, consistente com a QCD em rede e modelos de hadronização estatística.
  • O Modelo C (equilíbrio químico) descreve bem os dados.
  • O Modelo D (não-equilíbrio $\mu_\pi$) resulta em um $\mu_\pi$ compatível com zero e não oferece melhoria significativa na qualidade do ajuste sobre o Modelo C. O fator de Bayes ($B^D_C \approx 0,43$) indica que os dados não favorecem fortemente a extensão de não-equilíbrio quando os decaimentos de ressonância são explicitamente incluídos.

Principais Contribuições

1. Aplicação Sistemática do NSO: O artigo demonstra a aplicação do método NSO de Zubarev à fenomenologia de HIC, justificando formalmente a introdução de potenciais químicos efetivos como multiplicadores de Lagrange para números de partículas conservados em um cenário de não-equilíbrio.
2. Explicação Alternativa para o Aumento de Baixo $p_T$: O estudo estabelece que um potencial químico de píons de não-equilíbrio (Modelo B) oferece uma alternativa eficiente à explicação convencional do aumento de baixo $p_T$ via decaimentos de ressonância. Ele descreve com sucesso a sobrepopulação de estados de baixo momento usando um framework térmico estendido por parâmetros de não-equilíbrio.
3. Validação de Decaimentos de Ressonância: Ao acoplar o modelo de blast-wave ao código de transporte SMASH, os autores confirmam que o tratamento explícito de decaimentos de ressonância e termalização (Modelo C) explica naturalmente o aumento de baixo $p_T$, tornando o potencial químico de não-equilíbrio explícito (Modelo D) redundante nesse contexto específico.

Significância e Alegações
Os autores alegam que seu trabalho fornece uma abordagem mais fundamental para o problema do aumento de píons de baixo momento, indo além de correções fenomenológicas. Eles demonstram que a "sobrepopulação" de estados de baixo $p_T$ pode ser descrita por:

• Introduzir um potencial químico de não-equilíbrio em um modelo térmico simplificado (Modelo B), que serve como uma alternativa eficiente para rastrear processos microscópicos de termalização.
• Ou, por simular explicitamente decaimentos de ressonância e interações de estado final (Modelo C).

O artigo conclui que, embora o Modelo B seja uma ferramenta fenomenológica bem-sucedida para descrever espectros de congelamento cinético, ele é uma aproximação. Os autores observam que todos os modelos apresentados (particularmente C e D) permanecem aproximações porque ainda não incorporam totalmente efeitos de meio (como o efeito Mott ou correlações de contínuo), que são essenciais em altas densidades. O objetivo final declarado é estender a abordagem NSO para incluir observáveis de fluxo e dinâmica de formação de ressonâncias para fornecer um cenário consistente que justifique o sucesso e as limitações dos modelos estatísticos térmicos. O trabalho sugere que os parâmetros de fugacidade usados em modelos térmicos podem ser rigorosamente derivados como multiplicadores de Lagrange dentro do formalismo NSO.