Out-of-equilibrium contributions to charm hadrons in a fluid-dynamic approach

Resumo Técnico: Contribuições fora do equilíbrio para hádrons de charme em uma abordagem de dinâmica de fluidos

Enunciado do Problema
Os quarks de charme servem como sondas únicas para o Plasma de Quarks e Glúons (QGP) devido à sua grande massa, o que garante que sejam produzidos primariamente em espalhamentos duros iniciais e não alcancem o equilíbrio térmico via produção térmica dentro da vida útil da bola de fogo. Estruturas de dinâmica de fluidos anteriores modelaram com sucesso a dinâmica dos quarks de charme assumindo a proximidade do equilíbrio térmico local. No entanto, esses modelos negligenciaram correções fora do equilíbrio em dois estágios críticos: o momento inicial da colisão (após a fase de livre fluxo ou free-streaming) e a hipersuperfície de congelamento (freeze-out). A omissão das correções de congelamento levou a dependências não físicas das multiplicidades de hádrons de charme no coeficiente de difusão espacial ($D_s$), impedindo uma determinação sistemática dos coeficientes de transporte de charme. Além disso, a distribuição inicial era frequentemente assumida como em equilíbrio térmico local, ignorando a natureza fora do equilíbrio da produção inicial de charme.

Metodologia
Os autores estendem uma estrutura de dinâmica de fluidos para quarks de charme desenvolvida anteriormente, computando explicitamente as correções fora do equilíbrio tanto no estágio inicial quanto na superfície de congelamento.

1. Inicialização do Estágio Inicial:

   • O modelo introduz uma fase de livre fluxo entre o tempo de produção de charme e o tempo de inicialização hidrodinâmica ($\tau_0$).
   • Utilizando o método das características, os autores resolvem analiticamente as equações cinéticas para partículas viajando em trajetórias de linha reta durante esta fase.
   • Esta abordagem modifica a densidade de charme inicial e gera uma corrente de difusão não nula ($\nu^\mu_q$) em $\tau_0$, derivada da dependência espacial da densidade de colisões binárias ($n_{coll}$) e da distribuição de momento dos espalhamentos duros iniciais.
   • Condições de correspondência de Landau são empregadas para fixar o integral da função de distribuição total de charme em termos da seção de choque diferencial de produção.

2. Correções de Congelamento (Abordagem de múltiplas espécies):

   • Os autores empregam uma abordagem de múltiplas espécies dentro do Gás de Ressonâncias de Hadrons de charme (HRGc), tratando o sistema como um fluido composto por todas as espécies de hádrons de charme interagindo com componentes leves do QGP.
   • A correção fora do equilíbrio para a função de distribuição de cada espécie de hadron $i$ ($\delta f_{i,p}$) é derivada usando o método de momentos e a aproximação de Navier-Stokes.
   • A correção é expandida em momentos irredutíveis, truncada no posto $l=1$ (corrente de difusão), assumindo que as contribuições de pressão volumétrica e tensão de cisalhamento dos hádrons de charme são negligenciáveis em comparação com os componentes leves.
   • Um resultado analítico chave (Eq. 23) expressa a correção como:
     $$\delta f_{i,p} = f^{(0)}_{i,p} \frac{1}{P_i} \bar{\kappa}_i \nu^\mu_q p_{\langle \mu \rangle}$$
     onde $\bar{\kappa}_i$ é um coeficiente de difusão normalizado dependente da carga e densidade da espécie, e $\nu^\mu_q$ é a corrente de difusão de charme total. Esta formulação permite a evolução apenas da corrente total em vez das correntes individuais de cada espécie.

3. Evolução e Espectros:

   • O sistema é evoluído usando uma estrutura de dinâmica de fluidos 1+1D (FluiduM) com $D_s T$ constante.
   • Os rendimentos invariantes e os espectros de momento transversal ($p_T$) são computados via procedimento de Cooper-Frye, incorporando as correções fora do equilíbrio derivadas.
   • O modelo é testado contra dados experimentais do ALICE (colisões Pb-Pb em $\sqrt{s_{NN}} = 5,02$ TeV, centralidade 0–10%) usando tanto as previsões pQCD de FONLL quanto as seções de choque de charme medidas pelo ALICE como entradas.

Principais Contribuições e Resultados

• Resolução de Dependências Não Físicas: A inclusão de correções fora do equilíbrio no congelamento resolve a dependência não física das multiplicidades de hádrons de charme no coeficiente de difusão espacial $D_s$. Sem essas correções, a multiplicidade variava significativamente com $D_s$; com elas, os rendimentos integrados permanecem consistentes com o caso não difusivo e com o número total de quarks de charme iniciais em todo o intervalo de $D_s$ testado.
• Corrente de Difusão Inicial: A inicialização por livre fluxo gera uma corrente de difusão em $\tau_0$ que é uma ordem de magnitude maior do que as correntes produzidas em tempos tardios. No entanto, o estudo constata que essa condição inicial tem um impacto negligenciável na evolução subsequente dos campos de charme e nos rendimentos invariantes finais em comparação com um cenário de corrente inicial zero.
• Comparação Modelo-Dados:
  • Usando a seção de choque FONLL, o modelo reproduz os rendimentos integrados de mésons de charme dentro das incertezas.
    im Observa-se uma tensão no rendimento do bárion $\Lambda_c^+$, que é subestimado. Seguindo a literatura anterior, os autores aplicam um "aumento de bárion de charme" ao triplicar os pesos estatísticos de bárions de charme excitados. Com este ajuste e utilizando a seção de choque medida pelo ALICE, o modelo alcança compatibilidade com os rendimentos experimentais e espectros de $p_T$ até $p_T \approx 3$ GeV/c.
• Limites de Aplicabilidade: O artigo identifica um limite fundamental para a abordagem de dinâmica de fluidos baseado na positividade da função de distribuição.
  • A correção fora do equilíbrio $\delta f_{i,p}$ pode tornar-se negativa (devido ao perfil radial negativo da corrente de difusão).
  • Quando $|\delta f_{i,p}| \approx f^{(0)}_{i,p}$, a função de distribuição total torna-se negativa, tornando a interpretação física inválida.
  • Este "cruzamento por zero" do espectro de $p_T$ ocorre em valores de $p_T$ mais baixos conforme $D_s T$ aumenta. Para $D_s T \gtrsim 0,2$, o modelo produz espectros negativos em $p_T \sim 2$ GeV/c.
  • Consequentemente, o arcabouço é confiável apenas para $D_s T \lesssim 0,1$ em todo o intervalo de $p_T$, ou para $D_s T$ mais altos apenas dentro de uma região restrita de baixo $p_T$ (por exemplo, $p_T \lesssim 2,5$ GeV/c para $D_s T \approx 0,24$).

Significância
O artigo estabelece um arcabouço consistente de dinâmica de fluidos para hádrons de charme que considera efeitos fora do equilíbrio tanto no estágio inicial quanto no estágio final da colisão. Ao corrigir a dependência não física dos rendimentos em relação aos coeficientes de transporte, o trabalho lança a base necessária para futuras extrações sistemáticas do coeficiente de difusão de charme (por exemplo, via análise Bayesiana). No entanto, os autores concluem modestamente que a atual abordagem linearizada é limitada pela quebra da suposição perturbativa ($|\delta f| \ll f^{(0)}$) em coeficientes de difusão e momentos transversais mais elevados. Eles sugerem que trabalhos futuros devem explorar métodos alternativos, como a abordagem de máxima entropia, para garantir uma função de distribuição positiva em todo o intervalo cinemático.