Measurement of the Z $\to$ $\mu^+\mu^-$ angular coefficients in pp collisions at $\sqrt{s}$ = 13 TeV as functions of transverse momentum and rapidity

Dit artikel presenteert een meting van de acht hoekige polarisatiecoëfficiënten ($A_0$ tot $A_7$) voor Drell-Yan $Z \to \mu^+\mu^-$-productie in 13 TeV proton-protonbotsingen met behulp van 140 fb$^{-1}$ aan CMS-data, waarbij dubbel-differentieel resultaten in transversale impuls en rapiditeit worden geleverd die worden vergeleken met voorspellingen van QCD op de volgende-na-volgende-leidende orde.

De kern

Het Grote Plaatje: Een Kosmische Flipspelmachine

Stel je de Large Hadron Collider (LHC) bij CERN voor als 's werelds krachtigste flipspelmachine. Wetenschappers slaan protonen met bijna de lichtsnelheid tegen elkaar. Meestal veroorzaken deze botsingen een chaotische brij van deeltjes. Soms echter ontstaat er bij de botsing een zwaar, onstabiel deeltje genaamd een Z-boson, dat direct uiteenvalt in twee muonen (zware neven van elektronen).

Dit artikel gaat over CMS, een van de gigantische detectoren die dit flipspel in de gaten houden. Het team telde niet alleen hoe vaak dit gebeurde; ze wilden begrijpen hoe de muons eruit vlogen. Schoten ze rechtuit? Draaiden ze? Hebben ze de voorkeur gegeven aan één richting boven een andere?

Het Doel: Het Meten van de "Spin" van de Botsing

De wetenschappers maten acht verschillende getallen (gemarkeerd als $A_0$ tot en met $A_7$). Denk aan deze getallen als een gedetailleerd rapport over de "houding" of polarisatie van het Z-boson voordat het explodeerde.

• De Analogie: Stel je een vuurwerkexplosie voor in de lucht. Als het perfect symmetrisch explodeert, vliegen de vonken uit in een perfecte bol. Als het scheef staat of draait, kunnen de vonken meer naar links, meer naar boven, of in een spiraal vliegen.
• De Meting: De acht coëfficiënten ($A_0$ tot $A_7$) vertellen ons precies hoe die "vorm" van de explosie eruitziet. Ze onthullen of het Z-boson draaide, wiebelde, of in een bepaalde richting "uitgerekt" was.

Hoe Ze Het Dedden: De "Dubbelcheck"

Het team keek naar 140 biljoen botsingen (140 fb⁻¹ aan data) die werden opgenomen tussen 2016 en 2018. Ze keken niet alleen naar de hele hoop data; ze sneed het op als een brood om te zien of de "spin" veranderde afhankelijk van hoe hard de protonen tegen elkaar botsten.

1. Snelheid (Transversale Impuls): Ze keken naar muons die langzaam zijwaarts bewogen versus diegene die zeer snel bewogen.
2. Hoek (Rapidity): Ze keken naar muons die rechtuit vlogen versus diegene die onder een scherpe hoek vlogen.

Door de hoeken van de muons in deze specifieke stukken te meten, konden ze de acht coëfficiënten met extreme precisie berekenen.

De Regels van het Spel: De "Lam-Tung"-Regel

Het artikel bespreekt een beroemde regel in de natuurkunde, de Lam-Tung-relatie.

• De Analogie: Denk aan een regel die zegt: "Als je een bal recht omhoog gooit, moet hij recht naar beneden komen." In de wereld van de deeltjesfysica zouden, op het eenvoudigste niveau van berekening, twee van deze coëfficiënten ($A_0$ en $A_2$) elkaar perfect moeten opheffen ($A_0 - A_2 = 0$).
• De Realiteit: Het artikel bevestigt dat deze regel goed standhoudt bij lage snelheden, maar naarmate de botsingen energiekker worden (hogere impuls), begint de regel te bezwijken. Dit is geen falen; het is een kenmerk! Het vertelt ons dat de "rommelige" delen van de botsing (zoals extra deeltjes die worden uitgestoten) beginnen te tellen.

De Resultaten: Data versus Theorie

De wetenschappers vergeleken hun metingen met de beste beschikbare computersimulaties (de "theoretische voorspellingen").

• Het Goede Nieuws: Voor de meeste coëfficiënten stemde de data uit de echte wereld zeer goed overeen met de computermodellen. Dit betekent dat ons huidige begrip van hoe deze deeltjes met elkaar interageren stevig is.
• De Interessante Spanning: In het middengedeelte van de snelheden was de data voor één specifieke coëfficiënt ($A_0$) iets hoger dan de computer voorspelde (ongeveer 3 standaardafwijkingen afwijkend). Het is alsof een weersvoorspelling een kans van 50% op regen voorspelde, maar het daadwerkelijk 80% van de tijd regende. Het is geen ramp, maar het suggereert dat het computermodel misschien een klein detail mist.
• De "Geest"-Coëfficiënten: Drie van de coëfficiënten ($A_5, A_6, A_7$) zouden nul of zeer dicht bij nul moeten zijn. De data toonde aan dat ze inderdaad miniem waren, consistent met nul, hoewel één ervan ($A_6$) een klein, vaag hintje vertoonde van niet-nul te zijn. Dit is als een fluistering horen in een stille kamer; het is er, maar je hebt zeer gevoelige oren nodig om het te horen.

Waarom Dit Belangrijk Is

Dit artikel is in wezen een hoge precisie-kalibratiecheck voor de wetten van de natuurkunde.

1. Het Begrijpen van de "Lijm": Deze metingen helpen ons de "partonische dynamica" te begrijpen – hoe de kleine bouwstenen binnenin het proton (quarks en gluonen) zich gedragen wanneer ze botsen.
2. Het Testen van de Theorie: Door de "spin" van het Z-boson te vergelijken met complexe wiskunde (Kwantumchromodynamica), testen de wetenschappers ons begrip van het universum op zijn sterkte. Als de wiskunde niet overeenkomt met de data, betekent het dat we nieuwe natuurkunde moeten uitvinden.
3. De Referentiewaarde: Dit artikel biedt een nieuwe, ultra-precieze "liniaal" voor toekomstige experimenten. Elke nieuwe theorie moet deze acht getallen kunnen verklaren.

Samenvatting

Kortom, het CMS-team nam een enorme snapshot van deeltjesbotsingen, mat de exacte hoeken van de resulterende deeltjes en berekende acht getallen die de "spin" van de gebeurtenis beschrijven. Ze ontdekten dat hoewel onze huidige theorieën grotendeels correct zijn, er kleine, fascinerende afwijkingen zijn in het middensnelheidsbereik die fysici op hun tenen houden, waardoor het zoeken naar een dieper begrip van het universum doorgaat.

Technische samenvatting

1. Probleem en Motivatie

Het Drell–Yan-proces (DY) ($pp \to \gamma^*/Z \to \ell^+\ell^-$) is een fundamentele sonde voor het begrijpen van de onderliggende partonische dynamica van proton-protonbotsingen en de productie-mechanismen van elektroweak vectorbosonen. De hoekverdeling van de deeltjes in de eindtoestand in het ruststelsel van het boson bevat informatie over de polarisatietoestanden van het virtuele foton of het Z-boson.

Deze verdeling wordt geparametriseerd door acht hoekpolarisatiecoëfficiënten, $A_0$ tot en met $A_7$, gedefinieerd in het Collins–Soper (CS)-stelsel.

• Fysieke Doelen:
  • $A_0$ en $A_2$: Parametreren longitudinale polarisatie en interferentie tussen transversale toestanden. Hun verschil ($A_0 - A_2$) test de Lam–Tung-relatie ($A_0 - A_2 = 0$), die geldt op Leading Order (LO) in QCD maar naar verwachting wordt geschonden door hogere-orde QCD-straling en intrinsieke transversale impuls van partonen.
  • $A_4$: Gerelateerd aan de voor-achterwaartse asymmetrie ($A_{FB}$) en gebruikt om de effectieve zwakke mengingshoek ($\sin^2 \theta_{eff}^W$) te meten.
  • $A_1, A_3$: Bevatten interferentie tussen longitudinale en transversale toestanden en hangen af van axiale/vectoriële koppelingen.
  • $A_5, A_6, A_7$: Definieren T-odd asymmetrieën. Deze worden op LO en NLO verwacht nul te zijn, maar kunnen niet-nul zijn op Next-to-Next-to-Leading Order (NNLO) door een-lusprocessen.

Hoewel eerdere metingen bestonden bij $\sqrt{s}=8$ TeV en in het voorwaartse gebied (LHCb), was een uitgebreide, dubbel-differentiaal meting van de volledige set coëfficiënten ($A_0$–$A_7$) bij $\sqrt{s}=13$ TeV met hoge statistiek nodig om theoretische modellen te beperken en NNLO QCD-predicties met grotere precisie te testen.

2. Methodologie

Data en Selectie van Gebeurtenissen

• Dataset: Proton-proton botsingsdata verzameld door de CMS-detector bij $\sqrt{s} = 13$ TeV tijdens 2016–2018, overeenkomend met een geïntegreerde lichtopbrengst van 140 fb$^{-1}$.
• Signaal: Dilepton-gebeurtenissen ($Z \to \mu^+\mu^-$) met een invariant massa $81 < m_{\mu\mu} < 101$ GeV (Z-poolgebied).
• Kinematische Cuts:
  • Transversale impuls van het dimuon ($p_T^{\mu\mu}$): 0 tot 400 GeV.
  • Rapidity van het dimuon ($|y_{\mu\mu}|$): $< 2.4$.
  • Toepassing van muon-isolatie- en identificatiecriteria om hoge zuiverheid te waarborgen.
  • Trigger-eisen: Single-muon triggers met $p_T > 24$ of $27$ GeV, afhankelijk van de periode van dataverzameling.

Simulatie en Correcties

• Signaalsimulatie: Genereerd met POWHEG v2.0 + MiNNLOPS gekoppeld aan PYTHIA 8 (voor parton-showering/hadronisatie) en PHOTOS++ (voor QED-straling in de eindtoestand). De NNPDF3.1 PDF-set werd gebruikt.
• Schatting van Achtergrond: Monte Carlo-simulaties voor diboson (WW, WZ, ZZ), $t\bar{t}$ en $W$+jets-processen. Bijdragen van achtergronden zijn klein (0,03% bij lage $p_T$ tot ~2,9% bij hoge $p_T$).
• Correcties:
  • Schaalfactoren: Afgeleid via "tag-and-probe" om te corrigeren voor efficiënties van trigger, identificatie en isolatie.
  • Pileup Reweighting: Om de gesimuleerde pileup-verdeling te laten overeenkomen met de data.
  • Kinematische Reweighting: Gesimuleerde gebeurtenissen werden gewogen in $(p_T^{\mu\mu}, y_{\mu\mu})$-bakken om de kinematische verdelingen van de data te benaderen, waardoor modelafhankelijkheid wordt geminimaliseerd.
  • Prefiring Correctie: Adresseren van misalignement van bundelkruisingen in de Level-1 trigger.

Extractie van Coëfficiënten

De extractie maakt gebruik van een template-methode om rekening te houden met detectoracceptatie, resolutie en migratie tussen $p_T$-bakken.

1. Definitie: De differentiaalwerkzame doorsnede wordt uitgebreid als een som van trigonometrische polynomen $f_i(\theta^*, \phi^*)$ gewogen door coëfficiënten $P_i$ (waarbij $A_i = P_i/P_8$).
2. Templates: Negen 2D-templates ($T_i$) werden geconstrueerd uit gereconstrueerde MC-gebeurtenissen. Elke template komt overeen met een specifieke hoekterm $f_i$, gewogen om polarisatie-informatie uit de gebeurtenisgeneratie te verwijderen terwijl detector-effecten behouden blijven.
3. Fitten: Een twee-dimensionale gebinde maximum-likelihood fit (met MINUIT) werd uitgevoerd op de waargenomen hoekverdelingen ($\cos \theta^*, \phi^*$) in elk van de 16 kinematische bakken (8 $p_T$-bakken $\times$ 2 rapidity-bakken). De fit extrahert de coëfficiënten $P_i$ door het verschil tussen data en de lineaire combinatie van templates te minimaliseren.

3. Belangrijkste Bijdragen

• Eerste Volledige Set bij 13 TeV: Dit is de eerste meting van de volledige set hoekcoëfficiënten ($A_0$–$A_7$) in het centrale rapiditeitsgebied ($|y_{\mu\mu}| < 2.4$) bij $\sqrt{s}=13$ TeV.
• Dubbel-Differentiaal Precisie: Resultaten worden gegeven in 8 bakken van $p_T^{\mu\mu}$ en 2 bakken van $|y_{\mu\mu}|$, wat een multidimensionaal beeld van de hoekstructuur biedt.
• NNLO Vergelijking: De resultaten worden vergeleken met state-of-the-art theoretische voorspellingen op NNLO in QCD (met POWHEG+MiNNLOPS en MADGRAPH5_aMC@NLO).
• Beperkingen T-odd Asymmetrie: Biedt de strengste beperkingen tot nu toe voor de T-odd coëfficiënten ($A_5, A_6, A_7$) in het centrale gebied.

4. Resultaten

Algemene Overeenkomst

• De gemeten coëfficiënten komen over het algemeen overeen met NNLO QCD-predicties binnen de onzekerheden.
• De grote dataset (140 fb$^{-1}$) heeft de statistische onzekerheden aanzienlijk verminderd in vergelijking met eerdere Run 1-metingen.

Specifieke Observaties

• $A_0$ en $A_2$:
  • $A_0$ toont een afwijking van ongeveer 3 standaardafwijkingen in het intermediaire $p_T$-bereik (10–35 GeV) tussen data en de MADGRAPH5_aMC@NLO-predictie, terwijl POWHEG+MiNNLOPS goed overeenkomt.
  • $A_2$ wordt lichtjes overschat door MC-simulaties bij hoge $p_T$ ($>50$ GeV).
• Lam–Tung Relatie ($A_0 - A_2$):
  • De relatie wordt geschonden, zoals verwacht door hogere-orde QCD-effecten.
  • POWHEG+MiNNLOPS-predicties komen goed overeen met de gemeten $A_0 - A_2$-waarden in het bereik 10–35 GeV.
  • Beide MC-predicties onderschatten de meting voor $p_T > 35$ GeV.
• $A_1$:
  • Wordt significant overschat door POWHEG+MiNNLOPS in het gebied $1 < |y_{\mu\mu}| < 2.4$ voor $p_T < 35$ GeV.
  • MADGRAPH5_aMC@NLO overschat $A_1$ ook bij lage $p_T$ in beide rapiditeitsgebieden.
• T-odd Coëfficiënten ($A_5, A_6, A_7$):
  • Waarden zijn van de orde $10^{-3}$ en consistent met nul in de meeste bakken.
  • De meest significante afwijking is voor $A_6$ in de bak $55 < p_T < 80$ GeV en $1.2 < |y_{\mu\mu}| < 2.4$, met een afwijking van 2,8$\sigma$ van nul. Dit is consistent met NNLO-verwachtingen maar kleiner van grootte dan eerdere ATLAS-observaties bij 8 TeV.
• Systematische Onzekerheden:
  • Dominant voor $A_0$ en $A_1$ bij lage $p_T$ ($<55$ GeV) en voor $A_2, A_3, A_4$ tot $p_T < 200$ GeV.
  • Statistische onzekerheden domineren voor andere coëfficiënten en hogere $p_T$-gebieden.

5. Betekenis

• Theoretische Validatie: De resultaten bieden een rigoureuze benchmark voor NNLO QCD-berekeningen, wat bevestigt dat hogere-orde correcties essentieel zijn om de hoekstructuur van Drell–Yan-productie te beschrijven.
• Parton-dynamica: De metingen bieden nieuwe beperkingen voor Parton Distribution Functions (PDFs) en de modellering van intrinsieke transversale impuls van partonen.
• Toekomstige Fysica: De precisie en multidimensionale binningsindeling die hier zijn vastgesteld, dienen als een kritieke referentie voor toekomstige fenomenologische studies en voor het testen van zelfs hogere-orde berekeningen (bijv. N$^3$LO).
• Detectorprestaties: Demonstreert het vermogen van de CMS-detector om subtiele hoekcorrelaties met hoge precisie te meten, zelfs in aanwezigheid van complexe detector-effecten en pileup.

Samenvattend vertegenwoordigt dit artikel een mijlpaal in precisie elektroweak-fysica, waarbij de geldigheid van NNLO QCD-beschrijvingen van het productie-mechanisme van het Z-boson wordt bevestigd, terwijl specifieke kinematische gebieden worden belicht waar theoretische modellen mogelijk verdere verfijning vereisen.