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

Hoog boven de aarde regenen er constant onzichtbare deeltjes neer, ontstaan wanneer kosmische stralen de atmosfeer raken. Onder deze deeltjes bevinden zich neutrino's, spookachtige deeltjes die zelden met iets interageren, maar die wel geheimen in zich dragen over de meest gewelddadige gebeurtenissen in het universum. Om te begrijpen waar deze neutrino's vandaan komen, moeten wetenschappers hun oorsprong herleiden naar het moment waarop ze werden geboren: het verval van zware deeltjes zoals charm- en beauty-hadrons. Deze zware deeltjes worden gevormd wanneer hogenergetische bundels van protonen of pionen tegen een doelwit botsen, een proces dat niet alleen in de bovenste lagen van de atmosfeer plaatsvindt, maar ook in gespecialiseerde laboratoria waar onderzoekers deeltjesbundels in dichte blokken materie afvuren. Het voorspellen van hoeveel van deze zware deeltjes worden gemaakt, en exact hoe ze bewegen na hun creatie, is essentieel voor het berekenen van de flux van neutrino's die onze detectoren bereiken. Als de voorspellingen niet kloppen, wordt de achtergrondruis in experimenten die zoeken naar nieuwe, verborgen deeltjes onmogelijk te onderscheiden van een werkelijke ontdekking.

Decennialang hebben natuurkundigen vertrouwd op computersimulaties om deze botsingen te modelleren, gebruikmakend van complexe software die fungeert als een virtueel laboratorium. Eén dergelijke methode, Pythia, is al decennia lang de standaardtool, maar het was oorspronkelijk afgestemd op gegevens van botsingen met hoge energie bij de Large Hadron Collider, waar deeltjes op elkaar botsen bij energieën die veel hoger liggen dan wat fixed-target experimenten kunnen bereiken. Wanneer wetenschappers probeerden deze standaardinstelling te gebruiken om te voorspellen wat er gebeurt in fixed-target experimenten—waar bundels tegen stationaire doelwitten botsen bij lagere energieën—waren de resultaten vaak onjuist. De software had de neiging om te voorspellen dat de zware deeltjes te snel zouden bewegen en in grotere aantallen zouden worden geproduceerd dan wat daadwerkelijk in de data werd waargenomen. Deze mismatch creëerde een blinde vlek voor experimenten zoals SHiP, die zijn ontworpen om zeldzame, zwak interagerende deeltjes te zoeken door de resten van een massieve beam dump te analyseren. Zonder een nauwkeurigere kaart van hoe zware deeltjes worden geboren en zich gedragen in dit specifieke energiebereik, wordt de zoektocht naar nieuwe fysica gehinderd door onzekerheid.

Om dit op te lossen, heeft een team van onderzoekers van de Universiteit Gent een nieuwe set instructies ontwikkeld voor de Pythia-software, specifiek afgestemd op de omstandigheden die worden aangetroffen bij fixed-target botsingen. Ze noemen deze nieuwe configuratie FTFT. Het team heeft niet simpelweg geraden naar de juiste instellingen; in plaats daarvan hebben ze een enorme collectie echte metingen gebruikt uit experimenten die over de afgelopen veertig jaar zijn uitgevoerd. Deze metingen kwamen van bundels van zowel protonen als pionen die verschillende metalen doelwitten raakten, waarbij de snelheid en richting van de resulterende charm- en beauty-deeltjes werden geregistreerd. Door deze historische data in hun optimalisatietools te voeden, hebben de onderzoekers de interne knoppen van de simulatie bijgesteld totdat de virtuele botsingen overeenkwamen met de werkelijke botsingen. Ze concentreerden zich op twee hoofdzaken: de vorm van de deeltjesverdelingen, zoals hoe ver naar voren de deeltjes vliegen en hoeveel zijwaartse impuls ze meenemen, en het totale aantal geproduceerde deeltjes.

De resultaten van deze afstemming waren opmerkelijk. De standaard softwareinstellingen, die goed werken voor botsingen met de hoogste energieën, slaagden er niet in om de fixed-target data te beschrijven, omdat ze een spectrum van deeltjessnelheden voorspelden dat te agressief was. De nieuwe FTFT-tune corrigeerde dit door te wijzigen hoe de simulatie de afbraak van zware quarks in detecteerbare deeltjes en hoe meerdere interacties binnen een enkele botsing plaatsvinden, afhandelt. De onderzoekers ontdekten dat de zware deeltjes in fixed-target botsingen anders fragmenteren dan in de LHC, wat een zachter, meer geleidelijk afbraakproces in de simulatie vereist. Ze moesten ook aanpassen hoe de simulatie schaalt met energie, waarbij het referentiepunt van de enorme energieën van de LHC werd verlaagd naar het bereik dat relevant is voor fixed-target bundels. Zodra deze vormparameters waren vastgelegd, paste het team een definitieve schalfactor toe om het totale aantal geproduceerde deeltjes te laten overeenkomen met de experimentele metingen. Voor protonenbundels moest de simulatie met een factor van ongeveer 2,5 worden versterkt, terwijl voor pionenbundels de versterking ongeveer 2,0 was. Deze factoren houden rekening met de complexe fysica die het standaardmodel bij deze energieën mist.

Voor de nog zwaardere beauty-deeltjes was de situatie iets anders. Omdat de productie van beauty zo zeldzaam en moeilijk in detail te meten is bij deze energieën, kon het team de vorm van de verdeling niet afstemmen. In plaats daarvan hebben ze simpelweg een schalfactor toegepast op het totale aantal geproduceerde beauty-paren. Verrassend genoeg was de standaard simulatie voor beauty al vrij dicht bij de data, waarbij slechts een kleine aanpassing van ongeveer 4% voor protonenbundels en 19% voor pionenbundels nodig was. Dit suggereert dat het fundamentele proces van het creëren van zware beauty-paren beter begrepen wordt dan het proces van het creëren van charm-deeltjes in dit regime. De nieuwe tune pakte ook een specifiek raadsel aan met betrekking tot de lading van de deeltjes. Bij botsingen met pionenbundels is er een bekende voorkeur voor het produceren van negatieve charm-deeltjes boven positieve, een fenomeen dat wordt gedreven door de interne structuur van de pion zelf. De standaard simulatie had moeite om deze asymmetrie te vatten, maar de nieuwe FTFT-tune bracht, door een specifieke set regels voor de structuur van pionen op te nemen, de voorspelling veel dichter bij de realiteit, hoewel een kleine discrepantie bleef bestaan.

De implicaties van dit werk zijn onmiddellijk merkbaar voor de volgende generatie natuurkundige experimenten. Het SHiP-experiment, dat momenteel wordt gebouwd om op zoek te gaan naar verborgen deeltjes, vertrouwt zwaar op nauwkeurige voorspellingen van neutrino-achtergronden. De oude simulaties onderschatten de productie van charm-deeltjes, wat betekende dat de schattingen van de achtergrond waarschijnlijk te laag waren. Met de FTFT-tune zijn de voorspellingen voor de neutrino-flux en de achtergrondpercentages nu geworteld in een model dat rigoureus is getest tegen decennia aan data. Dit geeft het SHiP-team een veel steviger fundament voor hun zoektocht, waardoor ze met meer vertrouwen het onderscheid kunnen maken tussen een achtergrondfluctuatie en een echt signaal van nieuwe fysica. Naast SHiP is deze tune ook relevant voor het begrijpen van het prompt-component van atmosferische neutrino's, die worden geproduceerd wanneer kosmische stralen de aardatmosfeer raken. Aangezien deze atmosferische stormen gepaard gaan met veel secundaire botsingen bij energieën die vergelijkbaar zijn met fixed-target experimenten, biedt de nieuwe tune een betere manier om de zware deeltjes te modelleren die bijdragen aan de neutrino-flux die door telescopen wordt gemeten.

De onderzoekers merken er voorzichtig bij op dat, hoewel deze tune een significante verbetering is, het geen definitief antwoord is voor alle mogelijke scenario's. De afstemming werd specifiek uitgevoerd voor het energiebereik van 20 tot 42 GeV, en de parameters werden aangepast om aan de beschikbare data te voldoen, die voornamelijk betrekking hebben op charm-mesonen. Het gedrag van charm-baryonen, de zwaardere neven van de mesonen, werd niet beperkt door de data en blijft een gebied van onzekerheid. Bovendien gaat de tune ervan uit dat de productie van deze deeltjes lineair schaalt met de grootte van de kern van het doelwit, een aanname die wordt ondersteund door bestaande data, maar die verder getest kan worden. Ondanks deze beperkingen vertegenwoordigt het werk een cruciale stap voorwaarts in het overbruggen van de kloof tussen hogenergetische deeltjesversnellerfysica en het lagere energiegebied van fixed-target experimenten. Door de simulatie te funderen in de realiteit van eerdere metingen, heeft het team een instrument geleverd dat een bron van onzekerheid verandert in een betrouwbare voorspeller, waardoor de zoektocht naar de verborgen deeltjes van het universum wordt geleid door de meest nauwkeurige kaart die beschikbaar is.