Global Bayesian Analysis of $\mathrm{J}/ψ$ Photoproduction on Proton and Lead Targets

Dit artikel presenteert een globale Bayesiaanse analyse van diffractieve $\mathrm{J}/\psi$-fotoproductie op proton- en looddoelwitten met behulp van een color glass condensate-raamwerk, waarbij wordt onthuld dat hoewel een gelijktijdige beschrijving van HERA- en LHC-data uitdagend is, het introduceren van een algemene $K$-factor het vermogen van het model om beide datasets te fitten aanzienlijk verbetert.

De kern

Stel je voor dat je het perfecte gebak probeert te bakken, maar je hebt twee heel verschillende recepten om te volgen: één voor een klein, delicaat cupcake'tje (dat een enkel proton voorstelt) en één voor een enorme, dichte bundt cake (die een zware loi-kern voorstelt).

In de wereld van de hogeregetische fysica gebruiken wetenschappers een theoretisch "receptenboek" genaamd de Color Glass Condensate (CGC) om te voorspellen hoe deze gebakjes zich gedragen wanneer ze worden geraakt door een lichtstraal (fotonen). Dit licht wordt gebruikt om een specif kind soort deeltje te creëren, een J/ψ (uitgesproken als "J-psi"), wat als een klein, zwaar kersentje bovenop de taart dient.

Het Probleem: Het Recept Past Niet Bij Beide Gebakjes

Lange tijd merkten natuurkundigen een frustrerend probleem op. Wanneer ze het CGC-recept gebruikten om de resultaten voor het cupcake'tje (proton) te voorspellen, werkte het perfect. De voorspellingen kwamen overeen met de gegevens van deeltjesversnellers zoals HERA en de LHC.

Echter, wanneer ze dat exacte zelfde recept gebruikten om de resultaten voor de bundt cake (loi-kern) te voorspellen, ging het mis. Het recept voorspelde dat de bundt cake veel te veel J/ψ-deeltjes zou produceren, vooral wanneer de botsingsenergie hoog was. Het was alsof het recept zei: "Voeg een kop suiker toe voor het cupcake'tje," en zonder de hoeveelheid te veranderen, ook zei: "Voeg een kop suiker toe voor de bundt cake," wat resulteerde in een gebak dat veel te zoet was.

De wetenschappers wilden weten: Is er een enkele set ingrediënten (parameters) die zowel het kleine cupcake'tje als de gigantische bundt cake tegelijkertijd kan verklaren?

Het Onderzoek: Een Bayesiaanse "Proeverij"

Om dit op te lossen, voerden de auteurs een Global Bayesian Analysis uit. Denk hierbij aan een superintelligente, geautomatiseerde proeverij.

1. De Ingrediënten (Parameters): Ze hadden een lijst met variabelen die ze konden aanpassen, zoals de "grootte" van het proton, hoe "luchtig" de binnenkant is, en hoe de ingrediënten mengen bij hoge snelheden.
2. De Simulator (De Emulator): Omdat het bakken van deze theoretische gebakjes een enorme hoeveelheid computerkracht vereist, bouwden ze een "slimme gokker" (een Gaussian Process emulator). Dit hulpmiddel leerde de uitkomst van het bakproces te voorspellen zonder dat de volledige, trage simulatie telkens opnieuw gedraaid hoefde te worden.
3. De Test: Ze draalden duizenden simulaties waarbij ze de ingrediënten aanpasten om te zien welke combinatie zowel het cupcake'tje als de bundt cake goed van smaak kon maken (om de experimentele data te matchen) op hetzelfde moment.

De Bevindingen: De "Magische Schalingsfactor"

Dit is wat zij ontdekten:

• Het Standaard Recept Faalde: Wanneer ze probeerden beide datasets te fitten met het standaard recept (zonder enige extra trucjes), lukte dat niet. De instellingen die het cupcake'tje perfect maakten, maakten de bundt cake te zoet (te veel deeltjes). De instellingen die de bundt cake perfect maakten, maakten het cupcake'tje te droog (te weinig deeltjes). De twee datasets leken verschillende "evolutiesnelheden" voor de energie te willen.
• De "K-factor" Oplossing: De doorbraak kwam toen ze een K-factor introduceerden. Stel je dit voor als een universele "volumeknop" of een "schalingsdraaiwijze" die je voor het hele recept omhoog of omlaag kunt draaien.
  • Toen ze deze draaiknop omlaag draaiden naar ongeveer 0,3 (wat betekent dat ze de voorspelde output met 70% verminderden), gebeurde er iets magisch.
  • Door de totale output te verlagen, werd het model gedwongen om de interne ingrediënten aan te passen (specifiek, de dichtheid van de "lijm" die de deeltjes bij elkaar houdt, te verhogen).
  • Deze hogere dichtheid creëerde een sterkere "nucleaire onderdrukking" (zoals een dichtere cake die weerstand biedt tegen het uit elkaar vallen), wat de groei van deeltjes in de loi-kern op natuurlijke wijze vertraagde.
  • Resultaat: Plotseling kon hetzelfde recept zowel het kleine proton als de grote loi-kern perfect beschrijven.

Wat Niet Werkte

De wetenschappers probeerden ook andere chique aanpassingen aan het recept, zoals:

• Het vormgeven van het proton van een gladde bal naar iets meer grillig.
• Het toevoegen of verwijderen van "hot spots" (klontjes energie) binnenin het proton.
• Het wegfilteren van hoogfrequente ruis.

Ze ontdekten dat geen van deze chique aanpassingen zo goed hielp als simpelweg de volumeknop omlaag draaien (de K-factor). De data gaf een sterke voorkeur aan de eenvoudige schalingsoplossing boven deze complexe structurele veranderingen.

De Kern van het Verhaal

Het artikel concludeert dat hoewel het Color Glass Condensate-raamwerk krachtig is, het momenteel een "correctiefactor" (de K-factor) nodig heeft om zowel protonen als zware kernen tegelijkertijd te beschrijven.

Dit suggereert dat ons huidige begrip van de "niet-perturbatieve" delen van het recept (de rommelige, complexe delen van hoe deeltjes aan elkaar binden) of de hogere-orde effecten (de subtiele chemische reacties in de oven) nog niet volledig begrepen wordt. De K-factor fungeert als een tijdelijke oplossing voor deze ontbrekende stukjes, waardoor de theorie nu in staat is de data te matchen, maar het wijst er tegelijkertijd op dat de onderliggende theorie verdere verfijning nodig heeft om uit te leggen waarom die knop zo ver naar beneden gedraaid moet worden.

Kortom: Dezelfde natuurkundige regels gelden voor zowel kleine als grote doelwitten, maar ons huidige wiskundige "recept" heeft een globale volumeregeling nodig om de verhoudingen voor beiden juist te krijgen.

Technische samenvatting

Technische Samenvatting: Mondiale Bayesiaanse Analyse van $J/\psi$ Fotoproductie op Proton- en Looddoelwitten

Probleemstelling
Het artikel behandelt een hardnekkige spanning in de hogere-energiefysica van de kwantumchromodynamica (QCD) met betrekking tot de beschrijving van diffractieve $J/\psi$-fotoproductie. Hoewel Color Glass Condensate (CGC)-gebaseerde modellen succesvol zijn in het beschrijven van exclusieve $J/\psi$-productie in $\gamma + p$-botsingen bij HERA en de LHC, overschatten ze doorgaans de dwarsdoorsneden voor coherente en incoherente productie in $\gamma + \text{Pb}$-botsingen bij hoge centrum-van-massa-energieën ($W$). Eerdere berekeningen, die worden beperkt door protongegevens, slagen er niet in om de waargenomen mate van nucleaire onderdrukking in looddoelwitten te reproduceren; zij voorspellen een snellere toename met $W$ dan de experimentele data aangeeft. De centrale vraag is of een enkele set modelparameters simultaan zowel $\gamma + p$- als $\gamma + \text{Pb}$-data kan beschrijven binnen het standaard CGC-raamwerk, of dat de spanning onverzoenlijk is zonder de theoretische aannames te wijzigen.

Methodologie
De auteurs voeren een mondiale Bayesiaanse analyse uit om de waarschijnlijkheidsverdelingen van modelparameters ($\theta$) te beperken door theoretische voorspellingen te vergelijken met experimentele metingen ($y_{\text{exp}}$).

• Theoretisch Raamwerk: De berekening maakt gebruik van het CGC-raamwerk, specif specifiek het McLerran-Venugopalan (MV) model op een initiële schaal $x_P = 0,01$, geëvolueerd via de JIMWLK-vergelijking. De verstrooiingsamplitude wordt berekend met behulp van het Boosted Gaussian-model voor de niet-perturbatieve vector-meson-golffunctie en een dipool-doelwitamplitude geconstrueerd uit Wilson-lijnen.
• Modeluitbreidingen: Naast de "standaard" opstelling (vaste parameters voor protonvorm, hot spots en UV-demping) introduceren de auteurs verschillende uitbreidingen om hun flexibiliteit te testen:
  • Een flexibele transversale protonvorm-parameter ($\omega$).
  • Een variabel aantal sub-nucleaire hot spots ($N_q$).
  • Een UV-dempingsfactor ($v_{\text{UV}}$) om hoogfrequente modi te onderdrukken.
  • Een globale schalingfactor ($K$) toegepast op alle dwarsdoorsneden om onzekerheden in de golffunctie of ontbrekende hogere-orde effecten op te vangen.
• Bayesiaanse Inferentie: De analyse maakt gebruik van de pocoMC sampler en Gaussian Process (GP) emulatoren (specifiek de PCGP-emulator uit het surmise-pakket en Scikit-learn GP) om het computationeel kostbare voorwaartse model te benaderen. De trainingsstrategie omvat Latin Hypercube Designs (LHD) en actieve leerprocessen om de nauwkeurigheid van de emulator over de 11-dimensionale parameterruimte te waarborgen.
• Datasets: De analyse incorporeert coherente en incoherente differentiële en geïntegreerde dwarsdoorsnededata voor $\gamma + p$ (H1, ZEUS, ALICE, LHCb) en $\gamma + \text{Pb}$ (ALICE, CMS) over een breed bereik van $W$.

Belangrijkste Bijdragen

1. Mondiaal Bayesiaans Raamwerk: Dit werk biedt de eerste simultane Bayesiaanse analyse van $\gamma + p$- en $\gamma + \text{Pb}$-diffractieve $J/\psi$-data binnen een verenigd CGC-raamwerk, waarbij modelonzekerheden systematisch worden voortgeplant.
2. Kwantificering van Spanning: De studie toont rigoureus aan dat het standaardmodel (zonder schalingfactor) niet beide datasets simultaan kan beschrijven. Afzonderlijke fits aan $\gamma + p$- en $\gamma + \text{Pb}$-data leveren tegenstrijdige beperkingen op voor de snelheid van de energie-evolutie en de verzadigingsschalen.
3. Modelvergelijking via Bayes-factoren: De auteurs gebruiken Bayes-factoren om modellen kwantitatief te vergelijken, waarbij ze de standaardopstelling vergelijken met diverse uitbreidingen. Deze statistische benadering bepaalt welke modeluitbreidingen daadwerkelijk door de data worden ondersteund versus welke slechts overgeparameteriseerd zijn.
4. Resolutie via Schalingfactor: De analyse identificeert dat de introductie van een globale schalingfactor $K$ de enige uitbreiding is die sterk wordt ondersteund door de data, terwijl andere structurele wijzigingen (protonvorm, aantal hot spots, UV-demping) niet significant verbeteren.

Resultaten

• Falen van het Standaardmodel: Wanneer $K$ vast staat op 1, kan het model beide systemen niet simultaan beschrijven. Een fit aan de $\gamma + p$-data overschat de $\gamma + \text{Pb}$-dwarsdoorsneden bij hoge $W$, terwijl een fit aan de $\gamma + \text{Pb}$-data de $\gamma + p$-dwarsdoorsneden onderschat. De posterieure verdelingen voor de globale fit vertonen een dubbelpiekstructuur in de JIMWLK-evolutieparameters, wat de spanning tussen de twee datasets aangeeft.
• Voorkeursmodel: Het model uitgebreid met een vrije schalingfactor $K$ levert een significant betere fit ($\chi^2/\text{dof} \approx 1,82$) vergeleken met het standaardmodel ($\chi^2/\text{dof} \approx 5,94$). De Bayes-factor bevoordeelt de inclusie van $K$ sterk ($\ln B_{A/B} \approx -18,56$).
• Parameterecten: De geprefereerde waarde voor de schalingfactor is $K \approx 0,3$. Deze reductie in de algehele normalisatie van de dwarsdoorsnede wordt gecompenseerd door een voorkeur voor een grotere kleurladendichtheid (kleinere $Q_s/(g^2\mu)$-ratio), wat leidt tot een grotere effectieve verzadigingsschaal.
• Fysische Interpretatie: De grotere verzadigingsschaal versterkt niet-lineaire effecten, wat resulteert in een sterkere nucleaire onderdrukking en een tragere energieafhankelijkheid voor het looddoelwit, wat overeenkomt met de experimentele observaties.
• Andere Uitbreidingen: Variaties in de protonvorm ($\omega$), het aantal hot spots ($N_q$) en de UV-demping ($v_{\text{UV}}$) bleken ofwel gedefavoriseerd of slechts zwak geprefereerd te zijn ten opzichte van het standaardmodel. Het toevoegen van deze parameters aan het met $K$ uitgebreide model verminderde de Bayesiaanse bewijskracht verder.

Betekenis en Claims
Het artikel concludeert dat binnen het huidige CGC-gebaseerde raamwerk een simultane beschrijving van $\gamma + p$- en $\gamma + \text{Pb}$-data uitdagend is en een significante herschaling van de dwarsdoorsneden vereist ($K \sim 0,3$). De auteurs interpreteren de noodzaak hiervan als een indicatie dat de onzekerheden in de niet-perturbatieve vector-meson-golffunctie en/of ontbrekende hogere-orde (NLO) correcties substantieel zijn en momenteel onveraccounted blijven in de leading-order formalisme.

De auteurs claimen bescheiden dat hun resultaten aangeven dat berekeningen die volledig in het lineaire regime worden uitgevoerd (of zonder deze specifieke niet-lineaire aanpassingen), niet beide datasets simultaan kunnen beschrijven. Zij suggereren dat NLO-correcties dit probleem wellicht kunnen oplossen, maar dat volledige NLO-berekeningen nog niet beschikbaar zijn om dit te bevestigen. De studie motiveert toekomstig werk om de theoretische beschrijving van vector-meson-productie te verbeteren, inclusief NLO-uitbreidingen en betere beschrijvingen van niet-perturbatieve fysica, en benadrukt de noodzaak van data van intermediaire kernsystemen (bijv. $\gamma + \text{O}$, $\gamma + \text{Xe}$) om verzadigingseffecten verder te verkennen.