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

Bien haut au-dessus de la Terre, des particules invisibles pleuvent constamment, nées lorsque les rayons cosmiques frappent l'atmosphère. Parmi elles se trouvent les neutrinos, ces particules fantomatiques qui interagissent rarement avec quoi que ce soit, mais qui portent les secrets des événements les plus violents de l'univers. Pour comprendre d'où viennent ces neutrinos, les scientifiques doivent remonter à l'origine de leur création : la désintégration d'hadrons lourds comme le charme et la beauté. Ces particules lourdes sont créées lorsque des faisceaux de protons ou de pions à haute énergie percutent une cible, un processus qui se produit non seulement dans la haute atmosphère, mais aussi dans des laboratoires spécialisés où les chercheurs projettent des faisceaux de particules dans des blocs de matière denses. Prédire combien de ces particules lourdes sont produites, et exactement comment elles se déplacent après leur création, est essentiel pour calculer le flux de neutrinos qui atteignent nos détecteurs. Si les prédictions sont erronées, le bruit de fond des expériences cherchant de nouvelles particules cachées devient impossible à distinguer d'une véritable découverte.

Pendant des décennies, les physiciens se sont appuyés sur des simulations informatiques pour modéliser ces collisions, utilisant des logiciels complexes qui agissent comme un laboratoire virtuel. L'un de ces programmes, Pythia, est l'outil standard depuis des décennies, mais il a été initialement ajusté pour correspondre aux données de collisions à haute énergie au Grand Collisionneur de Hadrons (LHC), où les particules s'entrechoquent à des énergies bien au-delà de ce que les expériences à cible fixe peuvent atteindre. Lorsque les scientifiques ont tenté d'utiliser ce réglage standard pour prédire ce qui se passe dans les expériences à cible fixe — où des faisceaux frappent des cibles stationnaires à des énergies plus basses — les résultats étaient souvent erronés. Le logiciel avait tendance à prédire que les particules lourdes se déplaceraient trop vite et seraient produites en plus grand nombre que ce qui était réellement observé dans les données. Ce décalage a créé un angle mort pour des expériences comme SHiP, qui sont conçues pour traquer des particules rares et faiblement interagissantes en analysant les débris d'un énorme dépôt de faisceau (beam dump). Sans une carte plus précise de la naissance et du comportement des particules lourdes dans cette gamme d'énergie spécifique, la recherche d'une nouvelle physique est entravée par l'incertitude.

Pour remédier à cela, une équipe de chercheurs de l'Université de Gand a développé un nouvel ensemble d'instructions pour le logiciel Pythia, spécifiquement adapté aux conditions rencontrées dans les collisions à cible fixe. Ils appellent cette nouvelle configuration FTFT. L'équipe n'a pas simplement deviné les bons réglages ; elle a plutôt utilisé une vaste collection de mesures réelles issues d'expériences menées au cours des quarante dernières années. Ces mesures provenaient de faisceaux de protons et de pions frappant diverses cibles métalliques, enregistrant la vitesse et la direction des particules de charme et de beauté résultantes. En injectant ces données historiques dans leurs outils d'optimisation, les chercheurs ont ajusté les curseలు internes de la simulation jusqu'à ce que les collisions virtuelles correspondent aux collisions réelles. Ils se sont concentrés sur deux aspects principaux : la forme des distributions de particules, comme la distance vers l'avant parcourue par les particules et la quantité de moment transverse qu'elles transportent, ainsi que le nombre total de particules produites.

Les résultats de cet ajustement ont été frappants. Les réglages standards du logiciel, qui fonctionnent bien pour les collisions de très haute énergie, ne parvenaient pas à décrire les données de cible fixe, prédisant un spectre de vitesses de particules trop agressif. Le nouvel ajustement FTFT a corrigé cela en modifiant la manière dont la simulation gère la fragmentation des quarks lourds en particules détectables et la façon dont les interactions multiples se produisent au sein d'une même collision. Les chercheurs ont découvert que les particules lourdes dans les collisions à cible fixe se fragmentent différemment de celles de l'LHC, nécessissant un processus de fragmentation plus doux et plus délicat dans la simulation. Ils ont également dû ajuster la manière dont la simulation évolue avec l'énergie, déplaçant le point de référence des énergies massives de l'LHC vers la plage pertinente pour les faisceaux à cible fixe. Une fois ces paramètres de forme verrouillés, l'équipe a appliqué un facteur d'échelle final pour faire correspondre le nombre total de particules produites aux mesures expérimentales. Pour les faisceaux de protons, la simulation devait être augmentée d'un facteur d'environ 2,5, tandis que pour les faisceaux de pions, le facteur était d'environ 2,0. Ces facteurs tiennent compte de la physique complexe que le modèle standard omet à ces énergies.

Pour les particules de beauté, encore plus lourdes, la situation était légèrement différente. Comme la production de beauté est si rare et difficile à mesurer en détail à ces énergies, l'équipe n'a pas pu ajuster la forme de la distribution. À la place, elle a simplement appliqué un facteur d'échelle au nombre total de paires de beauté produites. Étonnamment, pour la beauté, la simulation standard était déjà assez proche des données, ne nécessitant qu'un ajustement infime d'environ 4 % pour les faisceaux de protons et de 19 % pour les faisceaux de pions. Cela suggère que le processus fondamental de création de paires de beauté lourdes est mieux compris que le processus de création de particules de charme dans ce régime. Le nouvel ajustement a également traité un puzzle spécifique concernant la charge des particules. Dans les collisions avec des faisceaux de pions, il existe une préférence connue pour la production de particules de charme négatives par rapport aux positives, un phénomène dicté par la structure interne du faisceau de pions lui-même. La simulation standard peinait à capturer cette asymétrie, mais le nouvel ajustement FTFT, en incorporant un ensemble spécifique de règles sur la structure des pions, a rapproché la prédiction de la réalité, bien qu'un petit écart subsiste.

Les implications de ce travail sont immédiates pour la prochaine génération d'expériences de physique. L'expérience SHiP, qui est actuellement en cours de construction pour rechercher des particules cachées, repose lourdement sur des prédictions précises des bruits de fond de neutrinos. L'ancienne simulation sous-estimait la production de particules de charme, ce qui signifiait que les estimations du bruit de fond étaient probablement trop basses. Avec l'ajustement FTFT, les prédictions pour le flux de neutrinos et les taux de bruit de fond sont désormais ancrées dans un modèle qui a été rigoureusement testé par des décennies de données. Cela donne à l'équipe de SHiP une base beaucoup plus solide pour sa recherche, lui permettant de distinguer une fluctuation du bruit de fond d'un véritable signal de nouvelle physique avec une plus grande confiance. Au-delà de SHiP, cet ajustement est également pertinent pour comprendre la composante prompt des neutrinos atmosphériques, qui sont produits lorsque les rayons cosmiques frappent l'atmosphère terrestre. Puisque ces gerbes atmosphériques impliquent de nombreuses collisions secondaires à des énergies similaires à celles des expériences à cible fixe, le nouvel ajustement offre un meilleur moyen de modéliser les particules lourdes qui contribuent au flux de neutrinos mesuré par les télescopes.

Les chercheurs précisent avec prudence que, bien que cet ajustement soit une améliation significative, il ne constitue pas une réponse finale pour tous les scénarios possibles. L'ajustement a été réalisé spécifiquement pour la gamme d'énergie de 20 à 42 GeV, et les paramètres ont été ajustés pour correspondre aux données disponibles, qui concernent principalement les mésons de charme. Le comportement des baryons de charme, des cousins plus lourds des mésons, n'a pas été contraint par les données et reste une zone d'incertitude. De plus, l'ajustement suppose que la production de ces particules suit une mise à l'échelle linéaire avec la taille du noyau cible, une hypothèse soutenue par les données existantes, mais qui pourrait être testée davantage. Malgré ces limites, ce travail représente une étape cruciale pour combler le fossé entre la physique des collisionneurs à haute énergie et le régime des cibles fixes à plus basse énergie. En ancrant la simulation dans la réalité des mesures passées, l'équipe a fourni un outil qui transforme une source d'incertitude en un prédicteur fiable, garantissant que la recherche des particules cachées de l'univers soit guidée par la carte la plus précise disponible.