Neutron tagging of heavy lepton pair production at RHIC and LHC energies

Resumo Técnico: Marcação de Nêutrons na Produção de Pares de Léptons Pesados em Energias de RHIC e LHC

Problema e Motivação
A produção de pares de lépton-antilépton via fusão fóton-fóton ($\gamma\gamma \to \ell^+\ell^-$) em colisões ultraperiféricas (UPCs) de íons pesados serve como um teste rigoroso da Eletrodinâmica Quântica (QED) em campos eletromagnéticos fortes. Embora a produção de pares elétron-pósitron tenha sido extensamente estudada, descrições teóricas para pares de léptons mais pesados (múons e taus) em conjunto com topologias específicas de emissão de nêutrons permanecem menos desenvolvidas na literatura. Em UPCs, onde os núcleos passam um pelo outro em parâmetros de impacto maiores que a soma de seus raios, os intensos campos Coulomb Lorentz-boosted atuam como fontes de fótons virtuais. Esses fótons podem induzir a excitação eletromagnética mútua dos núcleos, levando à emissão de nêutrons via Ressonância de Dipolo Gigante (GDR) e outros mecanismos. Colaborações experimentais (STAR, ALICE, ATLAS, CMS) categorizam eventos baseando-se na detecção de nêutrons em Calorímetros de Grau Zero (ZDCs) em classes como $0n0n$ (sem nêutrons), $0nXn$ (nêutrons em um ZDC) e $XnXn$ (nêutrons em ambos os ZDCs). Um arcabouço teórico consistente é necessário para calcular as seções de choque para essas classes específicas marcadas por nêutrons para permitir a comparação direta com dados experimentais em diferentes energias de centro de massa ($\sqrt{s_{NN}} = 200$ GeV, 2.76 TeV e 5.02 TeV).

Metodologia
Os autores realizam cálculos de QED de ordem mais baixa (LO QED) para a produção de $\gamma\gamma \to \ell^+\ell^-$ ($\ell = e, \mu, \tau$) em colisões Au+Au e Pb+Pb. A metodologia envolve três componentes primários:

1. Probabilidade de Produção de Pares Livres: O processo é tratado considerando as trajetórias dos íons como fontes eletromagnéticas externas clássicas. A seção de choque é derivada de elementos da matriz S envolvendo diagramas de Feynman diretos e cruzados. Para lidar com as integrais multidimensionais envolvendo funções de Bessel, que exibem oscilações rápidas, os autores empregam uma técnica de Monte Carlo. Eles calculam uma função $F(q)$ (relacionada à transferência de momento) e a ajustam a uma forma exponencial simples, $F(q) = F(0)e^{-aq}$, usando ajuste de mínimos quadrados não lineares. Isso permite a integração analítica da dependência do parâmetro de impacto ($b$) para determinar a probabilidade de produção de pares, $P_{\ell^+\ell^-}(b)$.

2. Emissão de Nêutrons e Excitação Nuclear: A probabilidade de excitação nuclear é modelada integrando o fluxo de fótons multiplicado pela seção de choque de fotoabsorção total ($\sigma_{tot}$). A seção de choque de fotoabsorção é construída combinando um ajuste Lorentziano para a região da GDR (usando dados experimentais recentes para $^{197}$Au e $^{208}$Pb) com o modelo de Kossov para contribuições não-GDR (quasi-deuteron, ressonância $\Delta$, ressonância de Roper e Regge-Pomeron). O número médio de fótons absorvidos é usado dentro de um arcabouço de estatística de Poisson para determinar as probabilidades de emitir zero ($0n$), um ($1n$) ou pelo menos um ($Xn$) nêutron. Isso permite a categorização de eventos nas classes $0n0n$, $0n1n$, $1n1n$, $0nXn$e$XnXn$.

3. Probabilidade de Não-Interação Hadrônica: Para garantir que as colisões permaneçam ultraperiféricas, a probabilidade de nenhuma interação hadrônica, $P_{nohad}(b)$, é calculada usando o modelo de Glauber com distribuições de densidade nuclear de Woods-Saxon. Esta probabilidade é multiplicada pelas probabilidades de produção de pares e excitação para obter as seções de choque marcadas finais.

Os cálculos são realizados dentro das restrições de espaço de fase fiducial definidas pelos experimentos STAR, ALICE, ATLAS e CMS, incluindo cortes em pseudorapidez ($|\eta|$), momento transversal ($p_T$), massa invariante ($M_{\ell\ell}$) e momento transversal do par ($p_{T,\ell\ell}$).

Resultos Principais
O estudo apresenta seções de choque diferenciais e integradas para a produção de pares de elétrons, múons e taus, subdivididas por classes de emissão de nêutrons:

• Pares de Elétrons (STAR, $\sqrt{s_{NN}}=200$ GeV): A seção de choque inclusiva calculada é 5850 $\mu$b. A seção de choque marcada $XnXn$ é calculada como 274.5 $\mu$b, o que é aproximadamente 4-5% superior à medição do STAR (261 $\mu$b) e 24.5% superior à previsão do STARLight. As distribuições diferenciais para massa invariante e momento transversal mostram bom acordo com os dados do STAR, embora o modelo subestime levemente a região de alto $p_T$.
• Pares de Elétrons (ALICE, $\sqrt{s_{NN}}=2.76$ TeV): Para as regiões de baixa e alta massa, as seções de choque inclusivas calculadas (147 $\mu$b e 87.7 $\mu$b) estão 3-5% abaixo das medições do ALICE, mas 13-15% acima das previsões do STARLight. O canal $0n0n$ domina a seção de choque.
• Pares de Elétrons e Múons (ATLAS, $\sqrt{s_{NN}}=5.02$ TeV):
  • Para elétrons, a seção de choque inclusiva é 233 $\mu$b, sendo 8.3% acima da medição do ATLAS (215 $\mu$b).
  • Para múons, a seção de choque inclusiva é 36.5 $\mu$b, 7.0% acima da medição do ATLAS (34.1 $\mu$b).
  • Em ambos os casos, o canal $0n0n$ fornece a contribuição dominante. As distribuições diferenciais em massa invariante e rapidez geralmente concordam bem com os dados do ATLAS, com razões próximas da unidade, embora as desvios aumentem em massas invariantes mais altas.
• Pares de Tau (CMS, $\sqrt{s_{NN}}=5.02$ TeV): Os autores calculam uma seção de choque inclusiva de 1043 $\mu$b para a produção de $\tau^+\tau^-$. Este resultado está em estreito acordo com o cálculo LO de Shao e d'Enterria usando o fator de forma de carga nuclear (1060 $\mu$b), mas é significativamente superior aos resultados fiduciais extrapolados do CMS (580–850 $\mu$b). O canal $0n0n$ contribui com 802 $\mu$b.

Significância e Alegações
O artigo alega fornecer um arcabouço teórico consistente para o cálculo da produção de pares de léptons marcados por nêutrons através de uma ampla gama de massas de léptons e energias de colisão. Uma contribuição primária é a extensão da análise de marcação de nêutrons de pares elétron-pósitron para pares de múons e taus mais pesados, que receberam menos atenção na literatura. Os autores demonstram que sua abordagem de QED de ordem mais baixa, combinada com um modelo de emissão de nêutrons baseado em Poisson e exclusão hadrônica baseada em Glauber, produz resultados que são, em geral, consistentes com os dados experimentais de STAR, ALICE, ATLAS e CMS dentro das incertezas teóricas esperadas.

Os autores notam modestamente que não incluem correções de QED de ordem superior (como radiação de estado final), embora citem literatura sugerindo que esses efeitos são menores para as seções de choque totais. Eles também reconhecem que seu tratamento das frações de emissão de nêutrons como constantes através dos bins cinemáticos é uma simplificação, o que pode explicar as ligeiras discrepâncias na forma das distribuições diferenciais em massas invariantes mais altas. No entanto, o estudo valida a utilidade desta abordagem simplificada para descrever a maior parte da produção de pares de léptons em UPCs e fornece uma linha de base para futuros estudos de precisão.