Event-by-event fluctuations of elliptic flow in ultrarelativistic O+O collisions

Resumo Técnico: Flutuações de evento a evento do fluxo elíptico em colisões ultra-relativísticas de O+O

Definição do Problema
O artigo investiga a natureza do fluxo anisotrópico, especificamente o fluxo elíptico ($v_2$), em colisões de oxigênio-oxigênio (O+O) em $\sqrt{s_{NN}} = 5,36$ TeV. Embora a visão predominante para núcleos grandes (ex: Pb+Pb) seja que o fluxo elíptico surge primariamente da assimetria geométrica da região de sobreposição determinada pelo parâmetro de impacto, os autores postulam que as colisões O+O podem representar um regime distinto. Em sistemas menores, o número limitado de núcleons sugere que as flutuações evento a evento da deposição de energia podem dominar a geometria inicial, tornando as colisões O+O potencialmente mais análogas a sistemas pequenos (p+p) do que a colisões de íons pesados no que diz respeito à origem do fluxo. O estudo visa quantificar essas flutuações e determinar se um conjunto reduzido de modos de flutuação pode reproduzir com precisão as propriedades estatísticas do fluxo elíptico.

Metodologia
Os autores empregam um arcabouço numérico totalmente tridimensional combinando McDipper (estado inicial), MUSIC (evolução hidrodinâmica) e iSS (particulização).

• Estado Inicial: As configurações de núcleons para o núcleo de $^{16}$O são geradas usando Teoria de Campo Eficaz de Rede Nuclear (NLEFT). Essas configurações são rotacionadas aleatoriamente e colididas com um parâmetro de impacto aleatório. O código McDipper, baseado na teoria de campo de Condensado de Vidro Colorido (modo GBW), gera perfis de densidade de energia 3D $\Phi(x, y, \eta_s)$, incorporando hot spots subnucleônicos e flutuações de espessura.
• Hidrodinâmica: Os perfis servem como condições iniciais para o MUSIC, evoluindo o sistema com uma viscosidade de cisalhamento específica constante ($\eta/s = 0,16$) e viscosidade volumétrica dependente da temperatura até uma temperatura de congelamento de 155 MeV.
• Estrutura de Decomposição: Uma inovação metodológica chave é a decomposição do perfil de densidade de energia 3D em um perfil médio $\bar{\Psi}$ e um conjunto de modos de flutuação não correlacionados $\{\Psi_l\}$ com coeficientes $c_l$. Isso é realizado via uma decomposição de autovalores da matriz de autocorrelação das flutuações de densidade de energia. Esta abordagem é aplicada aqui pela primeira vez a estados iniciais 3D não invariantes por boost.
• Resposta Linear: Os autores utilizam um arcabouço de resposta linear onde observáveis (excentricidade $\varepsilon_2$ e fluxo $v_2$) são aproximados como funções lineares dos modos de flutuação. Eles comparam simulações "evento a evento" (5000 eventos por classe de centralidade) contra observáveis "linearizados" gerados através da amostragem de coeficientes distribuídos gaussianamente para um conjunto reduzido de modos.

Contribuições Principais e Resultados

1. Dominância das Flutuações: O estudo demonstra que, em colisões O+O, a excentricidade elíptica inicial ($\varepsilon_2$) é impulsionada primariamente pelas flutuações evento a evento da deposição de energia, em vez do parâmetro de impacto finito. Consequentemente, o estado inicial médio $\bar{\Psi}$ é quase de formato arredondado mesmo em colisões não centrais (centralidade de 20–30%), apresentando uma excentricidade média muito pequena ($\varepsilon_2 \approx 0,087$). Isso contrasta fortemente com as colisões de íons pesados, onde o estado médio é significativamente em forma de amêndoa.
2. Pesos de Flutuação: Descobriu-se que o peso relativo do estado médio ($\bar{w}$) em O+O é de apenas $\sim 2\%$, significativamente menor que os $\sim 13\%$ observados em colisões Pb+Pb. Isso confirma quantitativamente que as flutuações são a contribuição dominante para a geometria inicial em O+O.
3. Precisão da Decomposição de Modos: Os autores mostram que a distribuição de probabilidade conjunta de $\varepsilon_2$ e $v_2$ pode ser reproduzida com precisão razoável usando um pequeno conjunto de modos de flutuação (especificamente 10 modos). O "modelo linear" e a "aproximação de flutuação Gaussiana analítica" (G.f.a.) produzem resultados em estreito acordo com as simulações completas evento a evento para $\varepsilon_2 \lesssim 0,6$.
4. Limitações Não Lineares: Em excentricidades elevadas ($\varepsilon_2 > 0,6$), a aproximação linear desvia-se dos dados evento a evento devido à resposta hidrodinâmica não linear (especificamente termos cúbicos na relação $v_2 \approx \kappa_{2,2}\varepsilon_2 + \kappa_{2,222}\varepsilon_2^3$). Os modelos lineares falham em capturar a forma de "banana" da distribuição e a cauda mais larga de $v_2$ causada por essas não linearidades.
5. Descorrelação Longitudinal: O artigo observa que, embora o fluxo integrado possa ser capturado por poucos modos, observáveis sensíveis à estrutura longitudinal, como a descorrelação longitudinal da excentricidade ($r_2$), requerem um número muito maior de modos (centenas) para reproduzir os resultados evento a evento. Isso é atribuído à importância dos modos de flutuação ímpares em rapidez, que aparecem em ordens superiores no espectro de modos.

Significância e Alegações
O artigo alega que as colisões O+O não são meramente uma versão "em escala reduzida" das colisões de íons pesados. Em vez disso, a origem do fluxo elíptico em O+O é fundamentalmente diferente, sendo dominada por flutuações do estado inicial, em vez da geometria global do parâmetro de impacto. Isso sugere que as colisões O+O fornecem um laboratório único para estudar o surgimento de fenômenos coletivos onde as flutuações são o principal motor.

Os autores alegam modestamente que seu arcabouço de decomposição modo a modo é uma alternativa viável e de menor custo às simulações hidrodinâmicas 3D completas evento a evento para descrever o comportamento estatístico global do fluxo elíptico. Eles demonstram que um conjunto limitado de modos de flutuação é suficiente para reproduzir os valores médios e as distribuições de probabilidade de $v_2$ e $\varepsilon_2$, desde que se considere a dominância das flutuações sobre a geometria média. No entanto, eles reconhecem que a aproximação linear possui limitações em relação às respostas não lineares em altas excentricidades e que observáveis mais complexos (como a descorrelação longitudinal) requerem uma base de modos significativamente maior.