Resumo Técnico: Estudando a Matéria de QCD em Colisões de Íons Pesados usando o Mecanismo de Cálculo MUSES

Declaração do Problema

A equação de estado (EoS) da matéria de QCD quente e densa é fundamental para descrever colisões de íons pesados em todas as energias de colisão. No entanto, determinar a EoS apresenta dois desafios primários:

1. Limitações Teóricas: Cálculos de QCD em Rede (LQCD) de primeira ordem são restritos a potenciais químicos bariônicos reais ($\mu_B$) nulos ou pequenos devido ao problema do sinal de Fermi. Embora a LQCD forneça dados confiáveis em altas temperaturas e baixos $\mu_B$, modelos fenomenológicos são necessários para densidades mais altas. Garantir que esses modelos correspondam às restrições da LQCD de uma maneira termodinamicamente consistente permanece um obstáculo significativo.
2. Complexidade Dinâmica: Colisões de íons pesados são sistemas inerentemente dinâmicos e fora do equilíbrio. Comparar dados experimentais com a teoria requer simulações hidrodinâmicas que evoluem múltiplas cargas conservadas (número bariônico, estranheza e carga elétrica). Essas simulações requerem entradas de EoS que não sejam apenas termodinamicamente consistentes, mas também invertíveis (convertendo de $(T, \mu_X)$ para variáveis hidrodinâmicas naturais como densidade de entropia $s$ e densidades de carga $n_X$) e capazes de lidar com transições de fase, incluindo potenciais pontos críticos (CP) e transições de primeira ordem.

Além disso, conectar as restrições de colisões de íons pesados com estrelas de nêutrons requer um framework unificado capaz de lidar com o diagrama de fase quadridimensional $(T, \mu_B, \mu_S, \mu_Q)$, uma vez que os sistemas exploram diferentes regiões deste espaço, apesar de possuírem densidades bariônicas potencialmente semelhantes.

Metodologia

Os autores utilizam o Mecanismo de Cálculo MUSES (Modular Unified Solver of the Equation of State), especificamente o novo lançamento Calliope. Este framework fornece um sistema de gerenciamento de fluxo de trabalho composível que orquestra uma coleção de módulos de software para computar EoSs e observáveis.

Principais Módulos e Fluxos de Trabalho

O artigo foca em módulos projetados para fenomenologia de íons pesados, expandindo o conjunto anterior de 10 módulos para 17:

• Fase Hadrônica: O módulo HRG (Hadron Resonance Gas), baseado no Thermal-FIST, computa a termodinâmica em todo o espaço 4D $(T, \mu_B, \mu_S, \mu_Q)$. Ele suporta interações ideais, de volume excluído e de Van der Waals Quântico (QvdW).
• Expansões de Lattice QCD:
  • 4D-Taylor (BQS): Expande a pressão em potências de $\hat{\mu}_X = \mu_X/T$ usando suscetibilidades estimadas em continuidade.
  • 4D-TExS: Generaliza o esquema de expansão $T'$ para 4D, mapeando potenciais químicos cartesianos para coordenadas esféricas para estender o alcance de validade além das expansões de Taylor padrão.
• Modelos Fenomenológicos:
  • Ising 2D-TExS: Combina a expansão $T'$ de rede 2D com o escalonamento crítico de Ising 3D para modelar um ponto crítico e uma linha de primeira ordem no plano $(T, \mu_B)$.
  • Holografia (NumRelHolo): Utiliza um modelo de correspondência gauge/gravidade Einstein-Maxwell-Dilaton (EMD) para descrever matéria fortemente acoplada, calibrado para LQCD em $\mu_B=0$ e estendido para densidade finita. Inclui coeficientes de transporte.
• Módulo de Síntese: Um novo procedimento baseado em mistura estatística e critérios de estabilidade permite fundir duas EoSs diferentes (ex: Holografia e HRG) em uma única EoS termodinamicamente consistente, cobrindo um intervalo mais amplo de $(T, \mu_X)$. Este procedimento pode gerar transições de crossover, transições de primeira ordem ou pontos críticos.
• Inversor de EoS: Um módulo que inverte numericamente tabelas de EoS de coordenadas $(T, \mu_X)$ para variáveis hidrodinâmicas naturais $(s, n_X)$ ou $(\epsilon, n_X)$, essencial para solvers hidrodinâmicos Lagrangianos como o CCAKE.

Framework de Simulação

Os autores realizam simulações hidrodinâmicas viscosas relativísticas usando o framework NuclearConfectionary (especificamente AMPT para condições iniciais, CCAKE 2.0 para hidrodinâmica e SMASH para o afterburner hadrônico).

• Colisões: Colisões centrais de Au+Au em $\sqrt{s_{NN}} = 7.7, 19.6,$ e $39$ GeV.
• Variáveis: As simulações rastreiam partículas SPH (Smoothed Particle Hydrodynamics) individuais através do diagrama de fase, permitindo a análise de trajetórias dinâmicas, mapas de calor e distribuições de freeze-out.
• Coeficientes de Transporte: A viscosidade de cisalhamento é definida como uma constante $\eta T/(\epsilon+P) = 0.08$. A viscosidade volumétrica inclui um componente de escalonamento crítico mimetizado via suscetibilidade bariônica $\chi_B^2$.

Principais Contribuições

1. Lançamento do MUSES Calliope: O artigo introduz um novo conjunto de módulos e fluxos de trabalho especificamente adaptados para colisões de íons pesados, complementando os módulos existentes para estrelas de nêutrons.
2. Cobertura Estendida de EoS: Os autores demonstram a construção de EoSs que cobrem o espaço de fase relevante para o Beam Energy Scan (BES-II) do RHIC, incluindo regiões de alto $\mu_B$ e baixo $T$, onde as expansões de LQCD tipicamente falham.
3. Fusão Termodinâmica: A implementação do módulo de Síntese permite a fusão contínua de diferentes modelos de EoS (ex: Holografia e HRG) com pontos críticos ajustáveis, garantindo consistência termodinâmica através das fronteiras de fase.
4. Inversão e Estabilidade: O módulo Inversor de EoS é testado por sua capacidade de lidar com o mapeamento não linear entre $(T, \mu_X)$ e $(s, n_X)$, incluindo estratégias para lidar com transições de fase de primeira ordem e regimes metaestáveis/instáveis.
5. Análise de Trajetória Dinâmica: Ao utilizar um solver hidrodinâmico Lagrangiano, os autores mapeiam as trajetórias específicas das células de fluido através do diagrama de fase de QCD, identificando quais regiões são sondadas em diferentes energias de feixe e rapididades.

Resultados

• Validade e Cobertura da EoS:
  • A EoS 4D-TExS fornece uma cobertura significativamente melhor no plano $(s, n_B)$ comparada à expansão 4D-Taylor padrão, particularmente em temperaturas mais baixas e densidades bariônicas mais altas.
  • A EoS fundida Holografia+QvdW-HRG oferece a cobertura mais ampla, abrangendo efetivamente as fases hadrônica e desconfinada.
  • Simulações em $\sqrt{s_{NN}} = 19.6$ GeV usando a EoS 4D-Taylor mostram que apenas ~50-70% das células de fluido na rapidez média permanecem dentro da tabela de EoS válida, necessitando do uso de EoSs conformais de backup. Em contraste, os modelos 4D-TExS e Holografia+HRG mantêm >90% de cobertura na rapidez média para $\sqrt{s_{NN}} \ge 19.6$ GeV.
• Sensibilidade ao Ponto Crítico:
  • A localização do ponto crítico influencia significativamente quais energias de feixe são sensíveis à sua presença.
  • Para um ponto crítico em $\mu_B^c = 598$ MeV (consistente com previsões teóricas recentes), apenas a menor energia ($\sqrt{s_{NN}} = 7.7$ GeV) mostra uma fração significativa de células de fluido ($\sim$40%) passando perto da região crítica. Energias mais altas sondam a região crítica principalmente em grandes rapididades forward/backward.
  • Nenhum "lenteamento crítico" (agrupamento de trajetórias isentrópicas) significativo foi observado nas simulações, sugerindo que a região crítica pode estar orientada horizontalmente ao longo do eixo $\mu_B$, em vez do eixo $T$, ou que efeitos fora do equilíbrio suavizam o sinal.
• Transições de Primeira Ordem:
  • Simulações incluindo regimes metaestáveis e instáveis através de uma transição de primeira ordem não causaram falhas numéricas no solver CCAKE, pois a região de velocidade do som ao quadrado negativa ($c_s^2 < 0$) era extremamente pequena e de curta duração ($\langle \Delta \tau \rangle \approx 1.4 \times 10^{-5}$ fm/c).
  • No entanto, os autores observam que o tratamento atual da linha de primeira ordem é uma exploração preliminar e não uma teoria dinâmica microscópica completa de conversão de fase.
• Observáveis:
  • A escolha da EoS (e a localização do ponto crítico) teve um efeito mínimo na multiplicidade total de píons ($dN/dy$), sugerindo que parâmetros de normalização dominam esses observáveis.
  • A largura da distribuição de rapidez ($dN/d\eta$) foi encontrada como fortemente dependente do comprimento de espalhamento do estado inicial no AMPT, em vez dos detalhes da EoS.

Significância e Perspectivas

O artigo estabelece a infraestrutura necessária para estudos sistemáticos da matéria densa de QCD dentro de um framework unificado e de código aberto. Ao fornecer EoSs modulares e termodinamicamente consistentes que podem ser fundidas e invertidas, o mecanismo MUSES Calliope permite comparações mais robustas entre dados de colisões de íons pesados e modelos teóricos.

Os autores enfatizam que seu trabalho é um "primeiro passo" e uma "exploração preliminar", particularmente em relação ao tratamento de transições de fase de primeira ordem e à inclusão de escalonamento crítico em coeficientes de transporte. Eles declaram explicitamente que o tratamento sistemático de transições de primeira ordem em simulações dinâmicas requer mais trabalho teórico e numérico. A principal significância reside na capacidade de rastrear quais regiões do diagrama de fase de QCD são amostradas durante a evolução do fireball e em quantificar como a validade limitada de EoSs baseadas em expansão afeta as previsões hidrodinâmicas. Este framework está posicionado para facilitar futuras análises conjuntas conectando restrições de colisões de íons pesados e estrelas de nêutrons.