A Pythia8 Tune for Open Charm and Beauty Production in Fixed-Target Collisions

Muito acima da Terra, partículas invisíveis caem constantemente, nascidas quando raios cósmicos atingem a atmosfera. Entre elas estão os neutrinos, partículas fantasmagóricas que raramente interagem com qualquer coisa, mas que carregam segredos sobre os eventos mais violentos do universo. Para entender de onde vêm esses neutrinos, os cientistas devem rastrear suas origens até o momento em que nasceram: o decaimento de partículas pesadas como os hádrons charm e beauty. Essas partículas pesadas são criadas quando feixes de alta energia de prótons ou píons colidem com um alvo, um processo que acontece não apenas na atmosfera superior, mas também em laboratórios especializados onde pesquisadores disparam feixes de partículas contra blocos densos de matéria. Prever quantos desses hádrons pesados são produzidos, e exatamente como eles se movem após a criação, é essencial para calcular o fluxo de neutrinos que alcançam nossos detectores. Se as previsões estiverem erradas, o ruído de fundo em experimentos que buscam novas partículas ocultas torna-se impossível de distinguir de uma descoberta genuína.

Por décadas, físicos confiaram em simulações computacionais para modelar essas colisões, usando softwares complexos que atuam como um laboratório virtual. Um desses programas, o Pythia, tem sido a ferramenta padrão por décadas, mas foi originalmente ajustado para corresponder a dados de colisões de alta energia no Grande Colisor de Hádrons (LHC), onde as partículas colidem em energias muito além do que experimentos de alvo fixo podem alcançar. Quando os cientistas tentaram usar essa configuração padrão para prever o que acontece em experimentos de alvo fixo — onde feixes atingem alvos estacionários em energias mais baixas — os resultados eram frequentemente errados. O software tendia a prever que as partículas pesadas se moveriam rápido demais e seriam produzidas em maior número do que o realmente observado nos dados. Esse descompasso criou um ponto cego para experimentos como o SHiP, que são projetados para caçar partículas raras e fracamente interagentes através da análise dos detritos de um enorme despejo de feixe (beam dump). Sem um mapa mais preciso de como as partículas pesadas nascem e se comportam nesse intervalo de energia específico, a busca por nova física é prejudicada pela incerteza.

Para corrigir isso, uma equipe de pesquisadores da Universidade de Ghent desenvolveu um novo conjunto de instruções para o software Pythia, adaptado especificamente para as condições encontradas em colisões de alvo fixo. Eles chamam essa nova configuração de FTFT. A equipe não simplesmente adivinhou as configurações corretas; em vez disso, utilizaram uma vasta coleção de medições do mundo real de experimentos conduzidos ao longo dos últimos quarenta anos. Essas medições vieram de feixes tanto de prótons quanto de píons atingindo vários alvos metálicos, registrando a velocidade e a direção das partículas charm e beauty resultantes. Ao alimentar esses dados históricos em suas ferramentas de otimização, os pesquisadores ajustaram os controles internos da simulação até que as colisões virtuais coincidissem com as reais. Eles focaram em duas coisas principais: a forma das distribuições das partículas, como a distância que as partículas voam para a frente e quanto momento lateral elas carregam, e o número total de partículas produzidas.

Os resultados desse ajuste foram impressionantes. As configurações padrão do software, que funcionam bem para as colisões de altíssima energia, falharam em descrever os dados de alvo fixo, prevendo um espectro de velocidades de partículas excessivamente agressivo. O novo ajuste FTFT corrigiu isso ao alterar a forma como a simulação lida com a fragmentação de quarks pesados em partículas detectáveis e como ocorrem múltiplas interações dentro de uma única colisão. Os pesquisadores descobriram que as partículas pesadas em colisões de alvo fixo se fragmentam de forma diferente de como o fazem no LHC, exigindo um processo de fragmentação mais suave e gentil na simulação. Eles também tiveram que ajustar como a simulação escala com a energia, movendo o ponto de referência das energias massivas do LHC para a faixa relevante para feixes de alvo fixo. Uma vez que esses parâmetros de forma foram travados, a equipe aplicou um fator de escala final para igualar o número total de partículas produzidas às medições experimentais. Para feixes de prótons, a simulação precisou ser impulsionada por um fator de aproximadamente 2,5, enquanto para feixes de píons, o impulso foi de cerca de 2,0. Esses fatores explicam a física complexa que o modelo padrão ignora nessas energias.

Para as partículas ainda mais pesadas, as de beleza (beauty), a situação era ligeiramente diferente. Como a produção de beleza é tão rara e difícil de medir detalhadamente nessas energias, a equipe não pôde ajustar a forma da distribuição. Em vez disso, aplicaram apenas um fator de escala ao número total de pares de beleza produzidos. Surpreendentemente, para a beleza, a simulação padrão já estava bastante próxima dos dados, exigindo apenas um pequeno ajuste de cerca de 4% para feixes de prótons e 19% para feixes de píons. Isso sugere que o processo fundamental de criação de pares de beleza pesados é melhor compreendido do que o processo de criação de partículas charm neste regime. O novo ajuste também abordou um enigma específico em relação à carga das partículas. Em colisões com feixes de píons, há uma preferência conhecida pela produção de partículas charm negativas sobre as positivas, um fenômeno impulsionado pela própria estrutura interna do feixe de píons. A simulação padrão tinha dificuldade em capturar essa assimetria, mas o novo ajuste FTFT, ao incorporar um conjunto específico de regras sobre como os píons são estruturados, aproximou muito a previsão da realidade, embora uma pequena discrepância tenha permanecido.

As implicações deste trabalho são imediatas para a próxima geração de experimentos de física. O experimento SHiP, que está sendo construído para buscar partículas ocultas, depende fortemente de previsões precisas de fundos (backgrounds) de neutrinos. As simulações antigas subestimavam a produção de partículas charm, o que significava que as estimativas de fundo provavelmente estavam muito baixas. Com o ajuste FTFT, as previsões para o fluxo de neutrinos e as taxas de fundo estão agora fundamentadas em um modelo que foi rigorosamente testado contra décadas de dados. Isso dá à equipe do SHiP uma base muito mais sólida para sua busca, permitindo que eles distingam entre uma flutuação de fundo e um sinal verdadeiro de nova física com maior confiança. Além do SHiP, este ajuste também é relevante para entender o componente prompt dos neutrinos atmosféricos, que são produzidos quando raios cósmicos atingem a atmosfera terrestre. Como essas chuvas atmosféricas envolvem muitas colisões secundárias em energias semelhantes às de experimentos de alvo fixo, o novo ajuste oferece uma maneira melhor de modelar as partículas pesadas que contribuem para o fluxo de neutrinos medido por telescópios.

Os pesquisadores observam cautelosamente que, embora este ajuste seja uma melhoria significativa, não é uma resposta final para todos os cenários possíveis. O ajuste foi realizado especificamente para a faixa de energia de 20 a 42 GeV, e os parâmetros foram ajustados para se adequar aos dados disponíveis, que dizem respeito principalmente a mésons charm. O comportamento dos bárions charm, seus primos mais pesados, não foi restringido pelos dados e permanece uma área de incerteza. Além disso, o ajuste assume que a produção dessas partículas escala linearmente com o tamanho do núcleo do alvo, uma suposição apoiada por dados existentes, mas que pode ser testada futuramente. Apesar dessas limitações, o trabalho representa um passo crucial para preencher a lacuna entre a física de colisores de alta energia e o regime de alvo fixo de menor energia. Ao fundamentar a simulação na realidade das medições passadas, a equipe forneceu uma ferramenta que transforma uma fonte de incerteza em um preditor confiável, garantindo que a busca pelas partículas ocultas do universo seja guiada pelo mapa mais preciso disponível.