Inclusive $J/ψ$ productions in pp collisions at $\sqrt{s}=$ 5.02, 7, and 13 TeV with the PACIAE model

Este estudo utiliza o modelo PACIAE 4.0, incorporando contribuições de estados de cor-singlete e cor-octeto do NRQCD, decaimentos fracos e recombinação de clusters, para investigar a produção inclusiva de $J/\psi$ em colisões pp a 5,02, 7 e 13 TeV, demonstrando um bom acordo com dados experimentais e quantificando as contribuições relativas dos diversos mecanismos de produção e efeitos de espalhamento.

O essencial

Imagine que o universo é uma grande cozinha de alta velocidade, onde partículas subatômicas são os ingredientes. Neste artigo, os cientistas estão tentando entender como uma "sobremesa" muito específica e rara, chamada J/ψ (pronuncia-se "J-psi"), é feita quando duas partículas chamadas prótons colidem em velocidades incríveis.

Aqui está a explicação do que eles fizeram, usando analogias do dia a dia:

1. O Grande Experimento: Colisões de "Fóruns"

Os cientistas usaram um modelo de computador chamado PACIAE 4.0. Pense nele como um simulador de voo super avançado, mas em vez de aviões, ele simula partículas.

• O Cenário: Eles simularam colisões de prótons em três níveis de energia diferentes (5,02, 7 e 13 TeV), que são como acelerar carros em velocidades cada vez mais altas.
• O Objetivo: Descobrir exatamente como a "sobremesa" J/ψ é criada e quantas delas aparecem em diferentes situações.

2. A Receita: Três Maneiras de Fazer a Sobremesa

Antes deste estudo, os cientistas olhavam apenas para uma receita básica. Agora, eles olharam para três receitas diferentes que podem criar o J/ψ:

1. A Receita Oficial (NRQCD): É a forma "padrão" e mais comum. Imagine que dois chefs (partículas de luz chamadas glúons) jogam ingredientes (quarks) juntos e eles se encaixam perfeitamente para formar o bolo. A maioria dos bolos vem daqui.
2. O Acidente Feliz (Cluster Collapse): Às vezes, os ingredientes estão tão perto um do outro no caos da cozinha que eles se grudam sozinhos, sem seguir uma receita estrita. É como se dois pedaços de massa colassem acidentalmente e virassem um bolo.
3. O Presente de Outro (Decaimento de B-hádrons): Alguns bolos não são feitos diretamente na colisão. Eles são "presentes" deixados por outros bolos maiores e mais pesados que se desmancham (decaindo) e deixam o J/ψ sobrando. É como se um bolo gigante explodisse e deixasse pedaços de J/ψ no chão.

O que eles descobriram?

• A Receita Oficial é a campeã, produzindo a maioria dos bolos.
• Porém, conforme a energia da colisão aumenta (o forno fica mais quente), os "acidentes felizes" e os "presentes de outros" tornam-se um pouco mais comuns, enquanto a receita oficial perde um pouquinho de sua fatia do bolo total.

3. O Mapa da Cozinha: Centro vs. Bordas

Os cientistas olharam para a produção de bolos em duas áreas da cozinha:

• No Centro (Rapidez Média): Onde a colisão é mais direta.
• Nas Bordas (Rapidez Frontal): Onde os ingredientes são lançados para os lados.

A Descoberta: Nas bordas, a "Receita Oficial" e os "Acidentes Felizes" são ainda mais comuns. Por quê? Porque nas bordas, há uma densidade enorme de ingredientes pequenos e rápidos que se juntam facilmente. Já os "Presentes de Outros" (que exigem ingredientes pesados e caros) são mais raros nas bordas, porque é difícil encontrar esses ingredientes pesados ali.

4. O Caos Pós-Colisão: A "Briga" na Cozinha

Depois que o bolo é feito, ele não fica parado. Ele pode bater em outros ingredientes que voam pela cozinha.

• Colisões de Partículas (antes do bolo): Os cientistas viram que, antes do bolo ser formado, as colisões entre os ingredientes não mudam muito o resultado final. É como se os chefs se esbarrassem antes de começar a cozinhar, mas o bolo final saísse igual.
• Colisões de Hádrons (depois do bolo): Aqui está a surpresa! Depois que o bolo (J/ψ) é feito, ele pode bater em outros ingredientes (como pães ou frutas) que estão voando. Essas batidas quebram alguns bolos.
  • Resultado: Cerca de 8% dos bolos são destruídos por essas batidas posteriores. É como se, ao sair da cozinha, o bolo fosse esmagado por uma porta que fecha. Isso é crucial para entender por que vemos menos bolos do que o esperado em alguns experimentos.

5. Por que isso importa?

Antes, os cientistas tentavam entender a produção de J/ψ apenas olhando para colisões entre núcleos pesados (como chumbo), o que é como tentar entender a receita de um bolo olhando para uma festa gigante onde tudo está bagunçado.

Este estudo é importante porque eles olharam apenas para colisões simples (próton com próton), como se estivessem cozinhandos em uma cozinha limpa e controlada. Ao entender exatamente como a receita funciona sem a bagunça de uma colisão pesada, eles criaram uma "régua" perfeita. Agora, quando olham para colisões pesadas (o caos da festa), eles podem dizer com certeza: "Isso aqui é porque a receita é essa, ou é porque o calor da festa (o plasma de quarks e glúons) derreteu o bolo?"

Resumo em uma frase:
Os cientistas usaram um supercomputador para mapear todas as formas possíveis de criar uma partícula rara (J/ψ) em colisões simples, descobrindo que a maioria vem de uma receita padrão, mas que cerca de 8% são destruídas por batidas posteriores, ajudando a desvendar os segredos da matéria mais densa do universo.

Resumo técnico

Título: Produção Inclusiva de $J/\psi$ em Colisões $pp$a$\sqrt{s} = 5.02, 7$e$13$ TeV com o Modelo PACIAE

1. Problema e Contexto

A produção de mésons $J/\psi$ (o estado vetorial de quarkônio mais leve, composto por um par $c\bar{c}$) em colisões próton-próton ($pp$) é fundamental para a física de íons pesados. A supressão de $J/\psi$ em colisões núcleo-núcleo foi originalmente proposta como uma assinatura do Plasma de Quarks e Glúons (QGP). No entanto, para isolar os efeitos do QGP ("matéria nuclear quente"), é essencial compreender precisamente os efeitos da "matéria nuclear fria" e a produção de base em colisões $pp$, onde não há núcleo no estado inicial.

Embora existam vários modelos teóricos (Modelo de Singlete de Cor, NRQCD, Evaporação de Cor), há uma necessidade de simulações mais completas que integrem mecanismos perturbativos e não perturbativos, incluindo tanto contribuições de singlete quanto de octeto de cor, além de efeitos de recombinação e decaimentos de estados mais pesados. Estudos anteriores dos autores focaram apenas em processos de singlete de cor; este trabalho visa expandir essa análise para uma descrição inclusiva e mais realista.

2. Metodologia

Os autores utilizaram o modelo PACIAE 4.0 (Parton And Hadron Cascade Evolution), que é uma extensão do gerador de eventos PYTHIA 8.3. O modelo opera em quatro estágios sequenciais:

1. Inicialização de Partões: Geração do estado inicial de partões, desativando a fragmentação de cordas inicial e quebrando diquarks.
2. Espalhamento de Partões (PRS - Partonic Rescattering): Ocorre antes da hadronização. Inclui espalhamentos $2 \to 2$ modelados por QCD perturbativo (pQCD) de ordem líder, com correções empíricas (fator $K$).
3. Hadronização: Ocorre via fragmentação de cordas de Lund (padrão no PYTHIA) ou colapso de clusters. O PACIAE 4.0 incorpora o mecanismo de colapso de clusters (cluster collapse), onde um par $c\bar{c}$ próximo no espaço de fases forma um estado ligado diretamente se a energia potencial da corda for insuficiente para criar um par $q\bar{q}$ leve.
4. Espalhamento Hadrônico (HRS - Hadronic Rescattering): Ocorre após a hadronização. Simula colisões elásticas e inelásticas entre hádrons. Para o $J/\psi$, canais inelásticos específicos (ex: $J/\psi + \pi \to D + \bar{D}^*$) são incluídos, permitindo a absorção ou transformação do $J/\psi$.

Abordagem Teórica da Produção:

• NRQCD (Cromodinâmica Quântica Não Relativística): Inclui contribuições de singlete de cor e octeto de cor para a produção direta e para a formação de estados excitados ($\psi(2S)$, $\chi_{c0}$, $\chi_{c1}$, $\chi_{c2}$).
• Decaimentos: Considera o "feed-down" (contribuição indireta) de decaimentos de estados de quarkônio mais pesados e de decaimentos fracos de hádrons com quark bottom ($b$-hadrons), que constituem a componente "não-prompt".
• Parâmetros: Foram gerados 400 milhões de eventos inelásticos não difrativos. Os fatores $K$ foram ajustados para 0,68, 0,71 e 0,75 para as energias de 5,02, 7 e 13 TeV, respectivamente, para alinhar com dados experimentais.

3. Contribuições Principais

• Expansão do Modelo: Primeira aplicação do PACIAE 4.0 a colisões $pp$ em múltiplas energias do LHC, incorporando simultaneamente mecanismos de singlete e octeto de cor dentro do formalismo NRQCD.
• Análise Inclusiva Completa: Separação quantitativa das contribuições de:
  • Produção direta (NRQCD).
  • Feed-down de estados excitados ($\psi(2S)$, $\chi_c$).
  • Colapso de clusters.
  • Componente não-prompt (decaimento de $b$-hadrons).
• Estudo de Dependência: Análise detalhada da dependência da produção em função da energia ($\sqrt{s}$) e da rapidez (central vs. frontal).
• Quantificação de Efeitos de Rescattering: Isolamento e medição do impacto específico do espalhamento de partões (PRS) e hadrônico (HRS) na yield final de $J/\psi$.

4. Resultados Chave

A. Concordância com Dados Experimentais:

• As seções de choque diferenciais em momento transversal ($p_T$) simuladas concordam bem com os dados experimentais (ALICE) nas rapididades central e frontal para todas as energias estudadas. A discrepância é geralmente inferior a 30%.

B. Contribuições Relativas (Energia e Rapidez):

• Dominância do NRQCD: O canal NRQCD é a principal fonte de $J/\psi$ (aprox. 74-86% da yield total).
• Tendência com Energia: À medida que a energia de colisão aumenta:
  • A fração relativa do NRQCD diminui.
  • As frações do colapso de clusters e do processo não-prompt (decaimento de $b$) aumentam. Isso ocorre porque a seção de choque para produção de pares $b\bar{b}$ cresce mais rapidamente com a energia do que a de $c\bar{c}$.
• Efeito de Rapidez:
  • Na rapidez frontal, a contribuição do NRQCD e do colapso de clusters é aumentada em comparação com a rapidez central.
  • A componente não-prompt é suprimida na rapidez frontal devido à escassez de partões com grande fração de momento ($x$), necessária para a produção de quarks pesados ($b$).

C. Decomposição do NRQCD:

• Baixo $p_T$: A produção direta e o decaimento de $\chi_{c2}$ dominam.
• Alto $p_T$: A produção direta combinada com o decaimento de $\chi_{c1}$ torna-se dominante.
• Outros Estados: O feed-down de $\psi(2S)$ e $\chi_{c0}$ contribui com frações menores (5-8% e 2-3%, respectivamente).
• Dependência de Energia nos Subcanais: A fração de produção direta aumenta modestamente com a energia, enquanto a contribuição de $\chi_{c1}$ e $\psi(2S)$ cresce, e a de $\chi_{c0}$ e $\chi_{c2}$ declina. Isso reflete a transição de mecanismos de espalhamento duro (singlete) para fragmentação (octeto) em energias mais altas.

D. Efeitos de Rescattering (PRS e HRS):

• Espalhamento de Partões (PRS): Tem um impacto negligenciável na yield de $J/\psi$ (variação < 1%). Isso ocorre porque o modelo não inclui espalhamento inelástico de $J/\psi$ com partões (ex: $J/\psi + g \to c + \bar{c}$).
• Espalhamento Hadrônico (HRS): Causa uma supressão de aproximadamente 8% na seção de choque inclusiva.
  • Suprime significativamente a produção via NRQCD e colapso de clusters devido a canais inelásticos com píons, nêutrons e prótons.
  • Não afeta a componente não-prompt (já que os hádrons $b$ não sofrem rescattering no modelo).

5. Significância

Este trabalho representa um avanço significativo na compreensão dos mecanismos de produção de quarkônio em colisões de alta energia. Ao integrar o modelo NRQCD completo (singlete + octeto) com dinâmicas de rescattering hadrônico e mecanismos de colapso de clusters, o estudo fornece uma base teórica robusta para:

1. Estabelecer uma linha de base precisa para colisões $pp$, essencial para interpretar a supressão de $J/\psi$ em colisões de íons pesados (efeitos de QGP).
2. Oferecer restrições quantitativas para os Elementos de Matriz de Longo Alcance (LDMEs) do NRQCD através de varreduras de energia e rapidez.
3. Demonstrar que o espalhamento hadrônico é um mecanismo de supressão relevante mesmo em colisões $pp$, o que deve ser considerado em análises de dados do LHC.

Os resultados são inteiramente novos e fornecem uma ferramenta quantitativa para desvendar a complexa interação entre processos perturbativos e não perturbativos na formação de quarkônios.