Comparisons of triple-differential cross sections for quasielastic-like $\nu_\mu$-hydrocarbon interactions using $\langle E_\nu\rangle \sim$ 3~GeV versus $\sim$ 6~GeV beams in MINERvA

Deze MINERvA-studie vergelijkt triple-differentiële dwarsdoorsneden voor quasi-elastische $\nu_\mu$-koolstof interacties in 3 GeV en 6 GeV bundels om neutrino-interactiemodellen te testen, waarbij discrepanties worden onthuld die wijzen op een overschatting van de eindtoestand-interacties voor protonen en geladen pionen in huidige simulaties.

De kern

Stel je voor dat je probeert uit te rekenen hoe snel een auto reed vlak voordat deze tegen een muur botste. Je kunt de auto niet meer zien, maar je kunt wel de snelheid van het puin meten dat van de muur wegvliegt en de hoek waaronder het de muur raakte. In de wereld van de deeltjesfysica doen wetenschappers iets soortgelijks met neutrino's—kleine, spookachtige deeltjes die bijna alles passeren.

Dit artikel gaat over een team van wetenschappers (de MINERvA-collaboratie) die een enorme detector hebben gebouwd om deze neutrino's te vangen en te bestuderen wat er gebeurt als ze tegen atomen botsen. Ze kijken specifiek naar een bepaald type botsing, genaamd "quasi-elastisch-achtig", waarbij een neutrino een kern raakt en deeltjes (zoals protonen) naar buiten slaat, waarbij de rest van de kern intact blijft maar wel geschokt wordt.

Hier is het verhaal van hun onderzoek, eenvoudig uitgelegd:

De twee verschillende "hamers"

Om hun theorieën te testen, gebruikten de wetenschappers niet slechts één enkele neutrino-bundel. Ze gebruikten twee verschillende "hamers" om het doelwit te raken:

1. De lage-energiebundel: Deze bundel is als een zachte tik. De neutrino's in deze bundel hebben een gemiddelde energie van ongeveer 3 GeV.
2. De medium-energiebundel: Deze bundel is een zware zwaai. De neutrino's hier zijn ongeveer twee keer zo energierijk, met een gemiddelde van 6 GeV.

De wetenschappers wilden zien of hun "instructiehandleiding" (de computermodellen die ze gebruiken om te voorspellen wat er gebeurt) op dezelfde manier werkte voor zowel de zachte tik als de zware zwaai.

Het mysterie van de "ontbrekende energie"

Wanneer een neutrino een atoom raakt, is het de bedoeling dat het specifieke deeltjes naar buiten slaat. Als je de snelheid en richting van de uitgaande deeltjes meet, zou je precies kunnen berekenen hoeveel energie de inkomende neutrino had. Het is als een perfect biljartspel waarbij je de snelheid van de keugel weet door te kijken naar waar de andere ballen naartoe gaan.

Echter, atomen zijn rommelig. Binnen de kern zijn deeltjes aan elkaar gebonden, en wanneer een botsing plaatsvindt, wordt het ingewikkeld:

• Sommige energie kan worden opgeslokt door de kern zelf.
• Sommige deeltjes kunnen vast komen te zitten of geabsorbeerd worden voordat ze kunnen ontsnappen.
• Soms komt een deeltje dat een proton zou moeten zijn, naar buiten als een neutron (wat onzichtbaar is voor hun detectoren).

Deze "ontbrekende" of "onzichtbare" energie maakt het moeilijk om precies te weten hoe snel de oorspronkelijke neutrino ging. Dit is een enorm probleem voor experimenten die neutrino-oscillaties (hoe neutrino's van smaak veranderen) bestuderen, want als je de beginenergie niet weet, kun je de verandering niet nauwkeurig meten.

Het onderzoek: De handleiding controleren

De wetenschappers maten het botsingspuin in zowel de lage-energie als de medium-energiebundels. Ze keken naar drie dingen bij elke botsing:

1. Hoe snel de muon (de "zus" van de neutrino) zijwaarts ging.
2. Hoe snel het naar voren ging.
3. De totale energie van alle zichtbare protonen die naar buiten vlogen.

Ze vergeleken hun gegevens uit de echte wereld met de voorspellingen van hun computermodellen (specifiek een programma genaamd GENIE).

De bevindingen: De modellen zaten ernaast

De resultaten toonden een duidelijke discrepantie tussen de echte wereld en de computermodellen:

• Het probleem van de "overschatting": De computermodellen voorspelden dat er meer hoogenergetisch puin zou zijn dan de wetenschappers daadwerkelijk zagen. Het is also kind dat het model dacht dat de botsing veel gewelddadiger was dan deze in werkelijkheid was.
• De "onzichtbare" boosdoener: De modellen leken te overschatten hoe vaak deeltjes worden geabsorbeerd of "opgeslokt" door de kern (Final State Interactions). Ze dachten dat protonen en pionen (een ander type deeltje) vaker rond stuiterden en vast kwamen te zitten dan ze in werkelijkheid deden.
• Het gaat niet alleen om snelheid: Interessant genoeg veranderde de fout niet veel simpelweg omdat de energie van de bundel veranderde van 3 GeV naar 6 GeV. De fout was consistent over beide bundels. Dit suggereert dat het probleem niet ligt in hoe de modellen de snelheid van de neutrino afhandelen, maar eerder in hoe ze de rommeligheid binnen de kern (de momentumoverdracht) afhandelen.

De "dubbele ratio"-truc

Om dit te bewijzen, gebruikten de wetenschappers een slimme truc. Ze namen de ratio van de lage-energiegegevens ten opzichte van de medium-energiegegevens, en deelden dat vervolgens door de ratio van de modellen voor diezelfde bundels. Deze "dubbele ratio" werkt als een vergrootglas.

Als de modellen perfect waren, zou deze ratio een vlakke lijn bij 1,0 zijn. In plaats daarvan zakte de lijn onder de 1,0 in specifieke gebieden. Dit bevestigde dat de modellen te veel gebeurtenissen voorspelden waarbij deeltjes werden geabsorbeerd, vooral wanneer het puin een hoge energie had.

De conclusie

De paper concludeert dat, hoewel de wetenschappers een goed begrip hebben van het algemene gedrag van neutrino's, de huidige computermodellen die door grote experimenten (zoals DUNE en NOvA) worden gebruikt, overschatten hoeveel energie er verloren gaat binnen de kern tijdens deze botsingen.

Ze kwamen tot de conclusie dat de modellen aangepast moeten worden om rekening te houden met het feit dat deeltjes niet zo vaak worden geabsorbeerd of "vast komen te zitten" als de huidige software nu denkt. Totdat deze modellen zijn gecorrigeerd, kunnen wetenschappers die de eigenschappen van neutrino's meten, er bij hun berekeningen net naast zitten, vergelijkbaar met het proberen te raden van de snelheid van een auto op basis van puin waarvan de computer denkt dat het verder is gevlogen dan het in werkelijkheid deed.

Kortom: De wetenschappers hebben een betere kaart gemaakt van het "verkeer" binnen de atoomkern. Ze ontdekten dat de huidige kaarten (modellen) te pessimistisch zijn over hoeveel verkeer er vast komt te staan, en dat ze moeten worden bijgewerkt om overeen te komen met de realiteit die wordt gezien in zowel lage- als hoge-energiebotsingen.

Technische samenvatting

Technische Samenvatting: Vergelijkingen van triple-differentiële dwarsdoorsneden voor quasielastische-achtige $\nu_\mu$-hydrocarbon interacties

Probleem en Motivatie
Neutrino-oscillatie-experimenten zijn sterk afhankelijk van geladen-stroom (CC) quasielastische-achtige (QE-achtige) verstrooiingsgebeurtenissen om de neutrino-energie te reconstrueren. Echter, deze interacties op kernen (nuclei) worden gecompliceerd door slecht begrepen nucleaire effecten, waaronder de beweging van gebonden nucleonen, kort-bereik multinucleon-correlaties (2p2h), en final-state interacties (FSI), waarbij geproduceerde hadronen ondergaan aan verstrooiing of absorptie binnen de kern. Deze effecten kunnen de relatie tussen de zichtbare final-state energie en de invallende neutrino-energie vertroebelen. Hoewel eerdere metingen discrepanties tussen data en interactiemodellen hebben geïdentificeerd, blijft de afhankelijkheid van deze modelleringsfouten van de neutrino-energie een kritieke onzekerheid. Dit artikel adresseert de noodzaak om interactiemodellen te testen over verschillende neutrino-energiespectra om te bepalen of mismodellering een functie is van de neutrino-energie of de momentumoverdracht.

Methodologie
De MINERvA-collaboratie heeft data geanalyseerd van twee verschillende breedbandige neutrino-bundels geproduceerd door 120 GeV protonen die grafietdoelen raken bij Fermilab:

1. Low Energy (LE) Beam: Piekt nabij 3 GeV, met weinig neutrino's boven de 6 GeV.
2. Medium Energy (ME) Beam: Piekt nabij 6 GeV, met weinig neutrino's boven de 10 GeV.

Beide datasets werden verzameld met dezelfde fijnmazige MINERvA-detector en de MINOS Near Detector voor muon-reconstructie. De analyse richt zich op CC QE-achtige gebeurtenissen ($\nu_\mu + A \to \mu^- + \text{nucleonen} + A'$) waarbij de final state alleen een muon en nucleonen bevat. De studie meet triple-differentiële dwarsdoorsneden in termen van:

• Muon transversale momentum ($p_t$)
• Muon longitudinale momentum ($p_{||}$)
• De som van de kinetische energieën van alle final-state protonen ($\Sigma T_p$), die de calorimetrisch zichtbare terugslagenergie (recoil energy) vertegenwoordigt.

De analyse maakt gebruik van een simultane gezamenlijke fit aan signaal- en achtergrondgedomineerde samples (gecategoriseerd door Michel-elektronen en geïsoleerde clusters om $\pi^+$ en $\pi^0$ te identificeren) om het niveau van de achtergrond te beperken. Achtergrond-gecorrigeerde event-rates worden ontbonden (unfolded) met behulp van een iteratieve techniek. Het primaire referentiemodel is een gemodificeerde versie van de GENIE event generator (v2.8.4 voor LE, getuned om equivalent te zijn aan v2.12.6 voor signaalprocessen), die specifieke MINERvA-tunes bevat voor pionproductie en multinucleon (2p2h) versterkingen.

Belangrijkste Bijdragen en Resultaten
Het artikel presenteert een directe vergelijking van triple-differentiële dwarsneden tussen de LE en ME bundels, een innovatieve aanpak die de energieafhankelijkheid van nucleaire effecten isoleert.

• Discrepanties in Dwarsdoorsneden: Significante discrepanties worden waargenomen tussen de data en de MINERVAs basismodel-voorspellingen. Specifiek overschat het model de dwarsdoorsnede-sterkte in de regio van hoge $\Sigma T_p$ (boven 200–300 MeV) en lage $p_t$. Deze overschatting is het meest uitgesproken in de LE-bundel voor $p_t < 400$ MeV/c en in de ME-bundel voor $p_t < 330$ MeV/c.
• FSI-modelleringsproblemen: De data wijst op een overschatting van de final-state interaction (FSI) sterktes in de simulatie. Dit manifesteert zich als een overvoorspelling van gebeurtenissen waarbij protonen neutronen uitstoten (wat resulteert in "onzichtbaar" energieverlies) en een overvoorspelling van pionproductie gevolgd door absorptie.
• Energieafhankelijkheid: Een cruciale bevinding is dat de discrepanties tussen data en voorspelling grotendeels consistent zijn over de twee verschillende bundelenergieën wanneer ze worden bekeken als een functie van $p_t$ en $\Sigma T_p$. De dubbele ratio van (Data/Model)$_{LE}$ tot (Data/Model)$_{ME}$ blijft consistent onder eenheid in regio's die gedomineerd worden door pionabsorptie en hoge $\Sigma T_p$. Dit suggereert dat de mismodellering primair een functie is van de momentumoverdracht en de specifieke nucleaire processen (zoals pionabsorptie), in plaats van een sterke functie van de invallende neutrino-energie.
• Modelvergelijkingen: De studie vergelijkt de data met diverse configuraties van GENIE (inclusclusief instellingen voor DUNE), NEUT, en NuWro. Geen van de alternatieve modellen slaagt erin de vorm van de dwarsdoorsnede bij lage $p_t$ of de specifieke energieafhankelijkheid waargenomen in de dubbele ratio's te reproduceren. Hoewel sommige modellen beter presteren bij intermediaire kinematica, falen zij om de kerndiscrepanties in de QE-achtige regio op te lossen.

Significantie
Het artikel concludeert dat het vermogen om de zichtbare energieverdeling van neutrino-interacties te modelleren momenteel wordt beperkt door de behandeling van final-state interacties en pionabsorptie, in plaats door het neutrino-energiespectrum zelf. Door aan te tonen dat de mismodellering persisteert over een factor-twee verschil in bundelenergie, bieden de resultaten een stringente benchmark voor huidige en toekomstige neutrino-event generators. De bevindingen impliceren dat verbeteringen aan oscillatie-analyses zich moeten richten op het verfijnen van de modellering van intranucleaire processen (specifiek FSI en pionabsorptie) om de neutrino-energie correct af te leiden uit quasielastische-achtige samples, in plaats van aan te nemen dat de fouten simpelweg schalen met de bundelenergie. Deze resultaten dienen als een directe test voor de interactiemodellen die worden gebruikt om neutrino-energieën af te leiden in oscillatie-experimenten.