Lindblad-driven quarkonium production in heavy-ion collisions

Resumo Técnico: Produção de quarkônio impulsionada por Lindblad em colisões de íons pesados

Declaração do Problema
Estados de quarkônio pesado servem como sondas de precisão do Plasma de Quarks e Glúons (QGP) formado em colisões de íons pesados ultrarelativísticos. Sua produção é governada pela interplay entre o screening de Debye do potencial de quark pesado e a decoerência no meio, levando a uma supressão hierárquica baseada na energia de ligação. Embora o potencial no meio seja conhecido por ser de valor complexo (com uma parte real descrevendo o screening e uma parte imaginária codificando o amortecimento de Landau), faltava um quadro teórico unificado que descrevesse simultaneamente tanto a supressão de estados pré-formados quanto a recombinação (coalescência) de quarks pesados termalizados a partir de primeiros princípios. Abordagens anteriores frequentemente tratavam a supressão e a regeneração como mecanismos fenomenológicos independentes.

Metodologia
Os autores empregam um quadro de sistema quântico aberto baseado na equação de Lindblad para descrever a evolução da matriz de densidade quark-antiquark pesado ($Q\bar{Q}$) acoplada ao banho do QGP. A metodologia procede através das seguintes etapas:

1. Potencial no Meio: O estudo utiliza o modelo da lei de Gauss (Lafferty e Rothkopf) para parametrizar o potencial complexo no meio $V(r, T) = \text{Re } V + i \text{Im } V$. Este modelo combina o potencial de Cornell no vácuo com a teoria de perturbação de Loop Térmico Rígido (HTL), governado por um único parâmetro dependente da temperatura, a massa de Debye $m_D(T)$. A parte imaginária é rigorosamente regularizada para remover divergências não físicas no termo da corda.
2. Funções Espectrais e Propriedades de Dissociação: Em vez de resolver a equação de Schrödinger dependente do tempo com um potencial imaginário dependente da posição (o que é numericamente desafiador), os autores resolvem a equação de Schrödinger no espaço de frequências. Eles extraem funções espectrais no meio calculando a parte imaginária da função de Green retardada traçada.
   • Energias de Ligação ($E_n$): Determinadas a partir das posições dos picos das funções espectrais após subtrair as massas dos quarks constituintes.
   • Larguras de Decaimento Térmico ($\Gamma_n$): Extraídas da largura total à meia altura (FWHM) dos picos, utilizando um perfil de Breit-Wigner assimétrico para contabilizar a assimetria próxima ao limiar do contínuo.
   • Temperaturas de Dissociação ($T_d$): Identificadas como a temperatura na qual o pico espectral se funde ao contínuo (energia de ligação desaparece).
3. Probabilidade de Sobrevivência: Para um sistema submetido à expansão de Bjorken, a probabilidade de sobrevivência de um estado pré-formado é calculada integrando a largura de decaimento dependente da temperatura $\Gamma(T(t))$ ao longo do tempo.
4. Modelo de Recombinação: O quadro é estendido para incluir recombinação projetando os operadores de salto estocástico da equação de Lindblad no subespaço de estados ligados. Sob uma aproximação adiabática (onde o Hamiltoniano efetivo evolui lentamente em comparação com as escalas de tempo de salto quântico), uma taxa de coalescência é derivada. Esta taxa é proporcional ao quadrado do número de quarks pesados livres ($N^2$) e ao volume do meio ($V$), modulada por um fator de relaxação para contabilizar o tempo finito necessário para os quarks pesados se termalizarem.

Principais Contribuições

• Quadro Unificado: O artigo deriva tanto os mecanismos de supressão quanto de recombinação a partir de uma única equação de movimento (a equação de Lindblad) para a matriz de densidade $Q\bar{Q}$, garantindo consistência interna entre a parte imaginária do potencial (que impulsiona o decaimento) e os operadores de salto (que impulsionam a recombinação).
• Derivação a Partir de Primeiros Princípios: A abordagem evita combinações ad hoc de modelos independentes. O potencial complexo fixa as larguras de decaimento térmico, que determinam diretamente tanto a probabilidade de sobrevivência quanto a normalização da taxa de recombinação via o teorema flutuação-dissipação.
• Coalescência Adiabática: Os autores derivam um modelo específico de coalescência a partir da equação de Lindblad, incorporando um fator de relaxação para corrigir a termalização não instantânea de quarks pesados no meio.

Resultados
O estudo aplica este quadro aos estados de quarkônio de charm ($J/\psi, \psi(2S)$) e bottomônio ($\Upsilon(nS)$) em colisões Pb–Pb a $\sqrt{s_{NN}} = 5.02$ TeV.

• Temperaturas de Dissociação:
  • Quarkônio de Charm: $T_{J/\psi} \approx 372$ MeV; $T_{\psi(2S)} \lesssim 155$ MeV (dissocia-se perto ou abaixo da temperatura crítica).
  • Quarkônio de Bottom: $T_{\Upsilon(1S)} \approx 564$ MeV; $T_{\Upsilon(2S)} \approx 223$ MeV; $T_{\Upsilon(3S)} \approx 164$ MeV; $T_{\Upsilon(4S)} < 155$ MeV.
• Fator de Modificação Nuclear ($R_{AA}$):
  • $J/\psi$: O modelo prevê uma contribuição significativa de recombinação em colisões centrais ($N_{part} \gtrsim 200$), impulsionada pela alta multiplicidade de quarks charm nas energias do LHC. Esta regeneração compensa parcialmente a supressão, resultando em um $R_{AA}$ total próximo à unidade, consistente com os dados do ALICE e CMS.
  • $\Upsilon(1S)$: Devido à multiplicidade de quarks bottom ser aproximadamente duas ordens de magnitude menor que a de charm, o termo de recombinação é negligenciável. O $R_{AA}$ para $\Upsilon(1S)$ é dominado quase inteiramente pela supressão, com os resultados do modelo concordando razoavelmente bem com os dados do CMS sem parâmetros de ajuste adicionais.

Significado e Afirmações
Os autores afirmam que este trabalho fornece uma "descrição unificada, inspirada em primeiros princípios" da produção de quarkônio. Ao derivar tanto a supressão quanto a recombinação a partir da mesma dinâmica microscópica (o potencial complexo e os operadores de salto de Lindblad), o quadro elimina a necessidade de ajuste fenomenológico independente desses efeitos concorrentes. O artigo afirma que, uma vez que o potencial complexo é fixado (via o modelo da lei de Gauss), todos os efeitos do meio sobre o rendimento do quarkônio seguem sem parâmetros livres adicionais. Isso representa um passo significativo em direção a uma descrição verdadeiramente a partir de primeiros princípios do quarkônio em colisões de íons pesados, oferecendo uma alternativa mais fundamentada teoricamente aos modelos de transporte tradicionais onde a supressão e a regeneração são frequentemente combinadas manualmente.