Extraction of baryon number susceptibilities at finite density from heavy-ion collisions

Resumo Técnico: Extração Bayesiana de Suscetibilidades do Número de Bárons em Colisões de Íons Pesados

Enunciado do Problema
As flutuações de cargas conservadas, especificamente o número de bárons, servem como sondas sensíveis da estrutura de fase da Cromodinâmica Quântica (QCD), particularmente próximo a um hipotético ponto crítico (CP). Embora cálculos de QCD em rede (lattice QCD) baseados em primeiros princípios determinem com precisão as suscetibilidades do número de bárons ($\chi_B^n$) em potencial químico de bárons nulo ($\mu_B = 0$), seu comportamento em densidade finita depende de extrapolações ou abordagens teóricas alternativas. Experimentalmente, o programa Beam Energy Scan (BES) do Relativistic Heavy Ion Collider (RHIC) mede flutuações evento a evento de números de prótons líquidos em colisões Au-Au para sondar este regime. No entanto, a interpretação quantitativa desses dados é complicada por vários efeitos de não-equilíbrio e experimentais: os prótons medidos carregam apenas uma fração do número de bárons total, a aceitação experimental é limitada no espaço de momento, e a conservação exata do número de bárons global modifica não trivialmente as flutuações em relação às expectativas do grande cânonico. Estratégias anteriores tipicamente confrontaram os dados com cálculos de modelos específicos (por exemplo, Gás de Ressonância de Hádrons, hidrodinâmica com volume excluído ou modelos de transporte), em vez de inferir diretamente as suscetibilidades subjacentes.

Metodologia
Este trabalho apresenta, tanto quanto aos autores sabem, a primeira extração bayesiana das suscetibilidades do número de bárons ($\chi_B^2, \chi_B^3, \chi_B^4$) diretamente de dados de colisões de íons pesados. O arcabouço constrói um mapa direto (forward map) das suscetibilidades teóricas para os cumulantes de prótons medidos através de quatro etapas distintas:

1. Fundo Hidrodinâmico: Superfícies de hipersuperfícies hidrodinâmicas médias por evento são geradas usando o código MUSIC com condições iniciais 3D e particlização em uma densidade de energia constante ($\epsilon_{sw} = 0,26$ GeV/fm$^3$).
2. Freeze-out de Máxima Entropia (MaxEnt): As suscetibilidades locais de bárons na hipersuperfície são convertidas em cumulantes conjuntos de bárons-antibarons. O arcabouço parametriza a entrada de suscetibilidade local como razões relativas a uma linha de base de um Gás de Ressonância de Hádrons ideal (iHRG): $\gamma = (\chi_B^2/\bar{\chi}_B^2, \chi_B^3/\chi_B^1, \chi_B^4/\chi_B^2)$. A prescrição MaxEnt mapeia os campos hidrodinâmicos e suas flutuações em cumulantes irredutíveis de multiplicidade de hádrons.
3. Aceitação e Filtragem: A aceitação de Cooper-Frye e a filtragem de bárons-para-prótons mapeiam os cumulantes conjuntos para prótons e antiprótons aceitos via amostragem binomial independente, contabilizando a aceitação cinemática do detector STAR.
4. Conservação Exata: O Método de Aceitação de Subensemble (SAM-3.0) impõe a conservação exata do número de bárons global, convertendo cumulantes de grande cânonico em cumulantes canônicos adequados para comparação com dados experimentais.

A inferência bayesiana é realizada sobre dados do STAR para colisões Au-Au de 0–5% de centralidade em energias de colisão específicas dentro do intervalo BES-II de $\sqrt{s_{NN}} = 7,7–27$ GeV (especificamente 7,7, 9,2, 11,5, 14,6, 17,3, 19,6 e 27 GeV). Duas extrações complementares são conduzidas usando (i) razões de cumulantes fatoriais de prótons ($\hat{C}_n/\hat{C}_1$) e (ii) razões de cumulantes de prótons líquidos ($\kappa_n/\kappa_1$). Priores planos são aplicados às razões de suscetibilidade, e uma verossimilhança gaussiana é construída usando incertezas estatísticas e sistemáticas.

Resultos Principais

• Suscetibilidade de Segunda Ordem ($\chi_B^2$): A análise produz restrições rigorosas sobre a suscetibilidade de segunda ordem normalizada pelo valor iHRG ($\chi_B^2/\bar{\chi}_B^2$), com incertezas relativas de 5–10%. Os valores extraídos exibem uma dependência não-monotônica da energia de colisão através das energias medidas: um aumento ($\sim 1,06–1,07$) na menor energia ($\sqrt{s_{NN}} = 7,7$ GeV), um mínimo ($\sim 0,73–0,82$) em energias intermediárias ($11,5–19,6$ GeV) e uma supressão plana ($\sim 0,8–0,9$) em energias mais altas.
• Comparação com QCD em Rede (Lattice QCD): Quando mapeados para a linha de freeze-out químico, os valores extraídos de $\chi_B^2$ mostram concordância quantitativa com estimativas baseadas em lattice QCD (especificamente a equação de estado 4D-TExS) para $\mu_B \lesssim 300$ MeV. No entanto, em $\mu_B$ maiores (correspondente a $\sqrt{s_{NN}} = 7,7$ GeV), os valores extraídos desviam-se para cima da banda de lattice, mostrando um aumento de até $\sim 50\%$.
• Suscetibilidades de Ordem Superior: As suscetibilidades de terceira e quarta ordem ($\chi_B^3, \chi_B^4$) são apenas fracamente restringidas devido a perdas de eficiência no mapeamento de bárons para prótons e distorções da conservação de bárons. Enquanto a razão de quarta ordem mostra uma preferência local por valores negativos em $\sqrt{s_{NN}} = 19,6$ GeV ($\chi_B^4/\chi_B^2 = -33 \pm 14$), as grandes incertezas impedem conclusões definitivas sobre a estrutura de fase da QCD a partir desses momentos de ordem superior isoladamente.
• Dominância de Correlação: Um cenário mínimo assumindo apenas correlações irredutíveis de dois bárons (definindo $\Delta\chi_B^3 = \Delta\chi_B^4 = 0$) reproduz com sucesso o pico não-monotônico observado na razão de cumulantes fatoriais de prótons $\hat{C}_3/\hat{C}_1$ próximo a $\sqrt{s_{NN}} \approx 11$ GeV. Isso sugere que a estrutura observada surge da interação entre a dependência energética de $\chi_B^2$ e a conservação exata do número de bárons, em vez de exigir correlações irredutíveis de três ou quatro bárons.

Significância e Alegações
O artigo afirma fornecer a primeira determinação bayesiana de suscetibilidades do número de bárons baseada em dados ao longo da linha de freeze-out químico, evitando suposições sobre uma Equação de Estado (EoS) específica. Os autores afirmam que seu arcabouço isola com sucesso a suscetibilidade de segunda ordem, encontrando-a consistente com a QCD em rede em densidades de bárons menores, mas indicando um potencial aumento nas flutuações do número de bárons em densidades mais altas ($\mu_B \gtrsim 300$ MeV) em relação a extrapolações não-críticas.

O estudo conclui modestamente que, embora os dados possam ser amplamente descritos sem correlações irredutíveis de múltiplos bárons além da segunda ordem, isso não descarta sua existência; em vez disso, as medições atuais são primariamente sensíveis à contribuição de dois bárons. Os autores observam que o aumento observado em baixas energias de colisão merece escrutínio com dados futuros de alvo fixo, e que a extração atual é condicional ao setup hidrodinâmico e de freeze-out fixo, representando uma estimativa de limite inferior para as incertezas em uma análise global. O arcabouço é apresentado como uma ferramenta geral para restringir a EoS de densidade finita da QCD e pode ser estendido para incluir energias de alvo fixo, correções de flutuação de volume e outras cargas conservadas uma vez que as incertezas de modelagem estejam melhor controladas.