Charmonium radiative transitions to dileptons from lattice QCD: The case of $h_c \to \eta_c \ell^+\ell^-$ and $\chi_{c1} \to J/\psi\,\ell^+\ell^-$

Dit artikel presenteert de eerste volledig dynamische rooster-QCD-berekeningen van de vervaltarieven en differentieel verdelingen voor de charmonium-overgangen $h_c \to \eta_c \ell^+\ell^-$ en $\chi_{c1} \to J/\psi\,\ell^+\ell^-$, waarbij de resultaten voor de $\chi_{c1}$-verval goed overeenkomen met experimentele data maar de voorspelling voor de $h_c \to \eta_c e^+ e^-$-verval ongeveer 3$\sigma$ afwijkt van de BESIII-meting.

De kern

Stel je voor dat deeltjesfysica een gigantisch, ingewikkeld bordspel is, waarbij de stukjes (de deeltjes) niet van hout of plastic zijn, maar van pure energie en krachten. In dit spel spelen zware deeltjes, genaamd charmonium, een belangrijke rol. Deze deeltjes zijn als een dansend koppel: een zware quark en zijn tegenhanger, een antiquark, die om elkaar heen draaien.

Deze wetenschappers hebben een nieuwe manier gevonden om te kijken hoe deze dansende koppel veranderen van vorm en energie, en ze hebben dit gedaan met een digitale simulatie die zo krachtig is dat het bijna net echt is.

Hier is het verhaal van hun ontdekking, vertaald in alledaags taal:

1. Het Probleem: De "Magische" Verdwijning

In de natuurkunde hebben we twee soorten "dansjes" die deze deeltjes kunnen doen:

• Het echte flitsje: Soms schiet een deeltje een echt foton (lichtdeeltje) uit, net als een cameraflits. Dit hebben we al vaak bestudeerd.
• Het virtuele flitsje (het nieuwe geheim): Soms schiet het deeltje geen echt licht uit, maar een "virtueel" foton. Dit is als een flits die nooit echt de camera bereikt, maar direct verandert in een paar nieuwe deeltjes: een elektron en een positron (of een muon en een antimuon).

Dit tweede soort dansje is lastiger te begrijpen. Het is alsof je probeert te voorspellen hoe een danser beweegt terwijl hij een onzichtbaar touw vasthoudt dat direct verandert in een nieuw danspaar. Tot nu toe waren onze theorieën hierover vaak gebaseerd op gissingen of benaderingen.

2. De Oplossing: De Supercomputer als "Tijdmachine"

De auteurs van dit paper hebben geen nieuwe deeltjesversneller gebouwd, maar ze hebben een supercomputer gebruikt om het heelal na te bootsen. Ze noemen dit Lattice QCD (Kwantumchromodynamica op een rooster).

• Het Rooster: Stel je het heelal voor als een gigantisch driedimensionaal raster (zoals een heel fijn honingraatpatroon).
• De Simulatie: Ze hebben de regels van de natuur (de krachten tussen de deeltjes) in dit raster gezet en laten de computer zien wat er gebeurt als deze dansende koppel veranderen.
• De "Twist": Ze gebruikten een speciale techniek (Twisted Mass) die ervoor zorgt dat de berekeningen niet "vervormen" door de grofheid van het raster, net zoals een goede lens een scherp beeld geeft zonder vervorming.

3. Wat hebben ze ontdekt? (De Dansstappen)

Ze hebben twee specifieke dansjes onderzocht:

1. De $\chi_{c1}$ naar $J/\psi$ dans: Dit is een dans waarbij het koppel van vorm verandert en een paar lichte deeltjes (elektronen of muonen) creëert.

   • Het resultaat: Hun berekeningen kwamen perfect overeen met wat de experimentatoren in China (BESIII) al hadden gemeten. Het was alsof hun voorspelling een perfecte kopie was van de realiteit. Dit geeft ons veel vertrouwen in de theorie.

2. De $h_c$ naar $\eta_c$ dans: Dit is een iets andere dans.

   • Het resultaat: Hier kregen ze een verrassing! Hun berekening gaf een resultaat dat ongeveer 3 keer zo groot was als wat de experimentatoren tot nu toe hadden gemeten.
   • Wat betekent dit? Dit is als een weerman die zegt: "Het regent 100 liter per uur," terwijl de regenmeter maar 30 liter aangeeft. Het kan zijn dat de meting nog niet precies genoeg is, of dat er iets raars gebeurt in de natuurkunde dat we nog niet begrijpen. Het is een spannende aanwijzing voor toekomstig onderzoek!

4. Waarom is dit belangrijk? (De "X-Factor")

Deze studie is uniek omdat ze niet alleen keken naar het totaal aantal deeltjes dat eruit komt, maar ook naar de snelheid en hoek waarmee ze eruit vliegen.

• De "Onzichtbare" Kracht: Bij het echte flitsje (licht) zie je alleen bepaalde bewegingen. Maar bij het virtuele flitsje (de deeltjesparen) kun je ook een "lengte-richting" meten die normaal gesproken onzichtbaar blijft.
• De Meting: Ze hebben nu een "handleiding" (formules) gemaakt die precies beschrijft hoe deze deeltjes zich gedragen. Dit is cruciaal voor twee redenen:
  1. Nieuwe deeltjes zoeken: Als er in de toekomst een nieuw, onbekend deeltje (zoals een "donker foton") zou bestaan, zou het deze dansjes verstoren. Nu we weten hoe het moet zijn volgens de standaardtheorie, kunnen we elke afwijking direct zien als een teken van nieuwe natuurkunde.
  2. Toekomstige experimenten: De meetapparatuur in de toekomst kan zich nu richten op de specifieke hoeken en snelheden die deze auteurs voorspellen.

Samenvatting in één zin

Deze wetenschappers hebben met een digitale supercomputer de "dansstappen" van zware deeltjes na-gebootst, wat hen in staat stelde om met hoge precisie te voorspellen hoe deze deeltjes veranderen in elektronenparen; hun voorspelling klopte perfect met het ene experiment, maar wees op een mogelijke mysterieuze afwijking bij het andere, wat de deur openzet voor nieuwe ontdekkingen in de deeltjesfysica.

Technische samenvatting

Probleemstelling en Context

Elektromagnetische overgangen in zware quarkonia (gebonden toestanden van een zware quark en zijn antiquark) bieden een gevoelige test voor de interne structuur en dynamica van deze systemen. Traditionele benaderingen, zoals effectieve veldtheorieën (NRQCD, pNRQCD), zijn afhankelijk van een expansie in de snelheid van de zware quark en vereisen vaak truncaties die systematische onzekerheden introduceren die moeilijk te kwantificeren zijn.

Hoewel radiatieve vervalprocessen met een reëel foton (bijv. $h_c \to \eta_c \gamma$) reeds bestudeerd zijn op het rooster, blijven processen waarbij een virtueel foton wordt uitgezonden dat vervolgens vervalt in een lepton-antilepton paar (Dalitz-verval, $h_c \to \eta_c \ell^+\ell^-$ en $\chi_{c1} \to J/\psi \ell^+\ell^-$) theoretisch minder goed begrepen.

• Het specifieke probleem: Bij virtuele fotonen is de impuls-overdracht $q^2$ niet nul. Dit opent de mogelijkheid voor longitudinale polarisatie van het virtuele foton, wat bij reële fotonen ($q^2=0$) afwezig is. Deze longitudinale componenten zijn afhankelijk van aanvullende vormfactoren die niet uit radiatieve vervaldata kunnen worden afgeleid.
• Experimentele situatie: Experimenten (voornamelijk door BESIII) hebben metingen gedaan van deze vervalbreedtes, maar er zijn discrepanties en onzekerheden. Vooral voor de $h_c \to \eta_c e^+e^-$ overgang is er een spanning tussen eerdere experimentele schattingen en theoretische verwachtingen. Er is behoefte aan "first-principles" berekeningen vanuit QCD zonder modelaannames om deze metingen te valideren en nieuwe fysica (zoals donkere sector deeltjes) te kunnen uitsluiten of detecteren.

Methodologie

De auteurs voeren een grootrooster QCD (Lattice QCD) studie uit om de relevante hadronische matrixelementen direct vanuit de QCD-Lagrangiaan te berekenen.

1. Configuraties en Fermionen:

   • Gebruik van $N_f = 2+1+1$ dynamische Wilson-Clover "twisted-mass" fermionen gegenereerd door de Extended Twisted Mass Collaboration (ETMC).
   • Simulaties op vier verschillende roosterafstanden ($a$) variërend van ongeveer 0,09 fm tot 0,057 fm.
   • De lichte quarkmassa's ($u, d, s$) en de charm-quarkmassa zijn ingesteld op hun fysische waarden (behalve op het ruwste rooster waar de pionmassa iets zwaarder is).
   • Dit maakt een gecontroleerde extrapolatie naar de continue limiet ($a \to 0$) mogelijk.

2. Berekening van Matrixelementen:

   • De studie focust op de overgangen $1^{+-} \to 0^{-+}$ ($h_c \to \eta_c$) en $1^{++} \to 1^{--}$ ($\chi_{c1} \to J/\psi$).
   • Er worden drie-punts correlatiefuncties berekend in het ruststelsel van de vervallende meson.
   • Vormfactoren: De auteurs extraheren de volledige $q^2$-afhankelijkheid van de vormfactoren:
     • Voor $\chi_{c1} \to J/\psi$: De multipoolvormfactoren $E_1(q^2)$ en $M_2(q^2)$, plus de nieuwe longitudinale vormfactor $C_1(q^2)$.
     • Voor $h_c \to \eta_c$: De vormfactoren $F_1(q^2)$ en $\tilde{F}_2(q^2)$.
   • Technische verbeteringen: Om statistische ruis te minimaliseren, vooral bij lage impulsen waar de signalen klein zijn, gebruiken ze geavanceerde technieken zoals "zero-momentum subtraction" en correlaties tussen verschillende impulsen. Ze gebruiken ook "partially-quenched twisted boundary conditions" om de impuls continu af te stemmen.

3. Extrapolatie en Fit:

   • De berekende vormfactoren worden geëxtrapoleerd naar de continue limiet met een lineaire fit in $a^2$.
   • Vervolgens worden de vormfactoren over het volledige kinematische bereik ($q^2$) gefit met een lineaire parametrisatie om de totale vervalbreedtes en differentieel spectra te berekenen.

Belangrijkste Resultaten

De studie levert de eerste volledig dynamische rooster-QCD voorspellingen op voor deze Dalitz-vervalprocessen.

1. Totale Vervalbreedtes (in keV):

• $\chi_{c1}$ kanalen:

  • $\Gamma(\chi_{c1} \to J/\psi e^+e^-) = 2.869(90)$ keV
  • $\Gamma(\chi_{c1} \to J/\psi \mu^+\mu^-) = 0.1993(72)$ keV
  • Vergelijking: Deze resultaten vertonen uitstekende overeenstemming met experimentele data van BESIII en PDG (binnen 1 standaardafwijking).

• $h_c$ kanalen:

  • $\Gamma(h_c \to \eta_c e^+e^-) = 5.45(19)$ keV
  • $\Gamma(h_c \to \eta_c \mu^+\mu^-) = 0.635(22)$ keV
  • Vergelijking: Hier is een significante spanning. De voorspelling voor het elektronkanaal is ongeveer 3$\sigma$ groter dan de recente BESIII-meting en ongeveer 2.5$\sigma$ groter dan de waarde afgeleid uit de PDG-gemiddelde.

2. Differentiële Verdelingen en Hoekobservabelen:

• De auteurs presenteren voorspellingen voor de differentieel vervalbreedte als functie van de invariante massa van het dilepton ($q^2$).
• Ze introduceren hoekobservabelen die specifiek gevoelig zijn voor de longitudinale vormfactoren ($C_1$ en $F_L$). Deze observabelen zijn niet toegankelijk via radiatieve vervalprocessen met reële fotonen en bieden een unieke test van de QCD-dynamica.
• Bij vergelijking van de differentieel verdeling van $\chi_{c1} \to J/\psi e^+e^-$ met BESIII-data (na correctie voor efficiëntie) is er over het algemeen goede overeenstemming, hoewel er in enkele lage-$q^2$ bins (rond 17, 37 en 47 MeV) spanningen van 2-3$\sigma$ optreden.

Significantie en Conclusies

1. Eerste Principes Berekening: Dit werk markeert de eerste volledig dynamische rooster-QCD berekening voor dilepton-vervalbreedtes in charmonium-overgangen, inclusief de volledige $q^2$-afhankelijkheid en de longitudinale componenten.
2. Benchmark voor Experimenten: De resultaten dienen als een robuuste benchmark voor toekomstige experimenten bij charm-fabrieken (zoals BESIII en toekomstige faciliteiten). De hoge precisie (enige procenten) maakt het mogelijk om kleine afwijkingen te detecteren.
3. Spanning bij $h_c$: De discrepantie tussen de roosterberekening en de experimentele meting voor $h_c \to \eta_c e^+e^-$ is een belangrijk resultaat. Het suggereert dat ofwel de experimentele metingen (of de afgeleide totale breedte van $h_c$) herzien moeten worden, of dat er mogelijk nieuwe fysica of onbegrepen systematische fouten in de experimentele analyse zitten.
4. Nieuwe Fysica Zoeken: De nauwkeurige vormfactoren die in dit werk zijn bepaald, zijn essentieel voor het zoeken naar "nieuwe mediators" (zoals donkere fotonen) in charmonium-vervallen. Zonder deze niet-perturbatieve QCD-input is het moeilijk om afwijkingen van het Standaardmodel betrouwbaar te identificeren.
5. Toekomstperspectief: De auteurs wijzen op de noodzaak van verdere experimentele metingen, vooral voor de $h_c$-kanalen en de differentieel verdelingen bij lage $q^2$, om de geobserveerde spanningen op te helderen.

Kortom, dit artikel levert een fundamentele bijdrage aan het begrip van elektromagnetische overgangen in zware quarkonia en biedt een cruciale theoretische basis voor het interpreteren van huidige en toekomstige experimentele data in de zoektocht naar nieuwe fysica.