Particle spectra in the integrated hydrokinetic model at RHIC Beam-Energy-Scan energies

Dit artikel onderzoekt de productie van lichte hadronen in Au+Au-botsingen bij RHIC-energieën met het uitgebreide geïntegreerde hydrokinetische model, waarbij wordt vastgesteld dat zowel een crossover als een eerste-orde faseovergang een vergelijkbare beschrijving van de spectra bieden, hoewel de grootste verschillen bij de laagste energie van 7,7 GeV optreden in de opbrengst van protonen en kaonen.

De kern

Stel je voor dat je twee enorme, zware vrachtwagens (de atoomkernen) tegen elkaar laat botsen, maar dan op een snelheid die bijna de lichtsnelheid bereikt. Dit gebeurt in een gigantische deeltjesversneller, de RHIC, in de Verenigde Staten. De wetenschappers in dit paper willen weten wat er precies gebeurt in die split-seconden na de botsing.

Ze kijken specifiek naar de "Beam Energy Scan" (BES), wat betekent dat ze de botsing niet altijd even hard laten gebeuren. Soms is het een zachte klap, soms een harde. Ze willen de "smakelijke" details van deeltjes die eruit vliegen begrijpen, vooral bij de lagere energieën.

Hier is de uitleg van hun onderzoek, vertaald naar alledaags Nederlands:

1. De Grote Kookpot: Quark-Gluon Plasma

Normaal gesproken zitten de bouwstenen van de materie (quarks en gluonen) gevangen in kleine cellen, net als kippen in een kooi. Deze cellen noemen we hadronen (zoals protonen).
Maar als je twee atoomkernen hard genoeg tegen elkaar botst, smelt de kooi. De kippen (quarks) lopen vrij rond in een soep. Dit noemen ze Quark-Gluon Plasma (QGP). Het is als een superheet, superdicht soepje van vrije deeltjes.

2. Het Probleem: De "Gordijnen" van de Tijd

Het probleem is dat dit soepje heel kort bestaat (minder dan een seconde, zelfs minder dan een miljardste van een seconde) en heel klein is. Het is alsof je probeert te begrijpen hoe een cake rijst, terwijl je alleen maar een flitsfoto maakt van het moment dat deeg in de oven gaat.

Om te weten wat er in die cake gebeurt, gebruiken de wetenschappers een computermodel. Ze noemen dit hun iHKMe-model. Dit is een soort "tijdmachine" die de hele reis van het deeg simuleert:

1. De botsing: De vrachtwagens raken elkaar.
2. Het wachten (Thermalisatie): Het deeg moet eerst even rusten en warm worden voordat het als soep begint te stromen.
3. Het stromen (Hydrodynamica): Het soepje zet uit en koelt af.
4. De afkoeling (Hadronisatie): Het soepje stolt weer tot vaste deeltjes (de kippen komen weer in kooitjes).
5. De uitloop: De deeltjes vliegen de detector in.

3. De Twee Soorten "Recepten" (Equations of State)

De wetenschappers twijfelen over één ding: hoe gedraagt dit soepje zich precies als het afkoelt? Ze proberen twee verschillende "recepten" (wiskundige regels) om te zien welke het beste werkt:

• Recept A (De Zachte Overgang): Het soepje verandert heel geleidelijk van vloeistof naar vast stof. Geen schok, gewoon een gladde overgang.
• Recept B (De Harde Schok): Het soepje ondergaat een plotselinge, harde verandering, alsof water ineens bevriest tot ijs. Dit zou kunnen betekenen dat er een "kritiek punt" is in het universum waar de regels veranderen.

4. Wat Vonden Ze?

Ze lieten hun computermodel draaien met deze twee recepten voor verschillende botsingsenergieën (van 7,7 tot 39 GeV).

• Bij hoge energieën (harde botsingen): Het maakt bijna niet uit welk recept je kiest. Het soepje zet zo snel uit dat de verschillen tussen de twee recepten verdwijnen. Het resultaat is hetzelfde.
• Bij lage energieën (zachte botsingen): Hier wordt het interessant! Bij de laagste energie (7,7 GeV) beginnen de twee recepten verschillende resultaten te geven. Vooral bij protonen en kaonen (soorten deeltjes) zie je verschillen. Het lijkt erop dat bij deze lagere energieën het "wachtgedeelte" (waar het deeg nog niet helemaal soep is) langer duurt en belangrijker wordt.

5. De Belangrijkste Leer: De "Wachttijd"

Een van de belangrijkste ontdekkingen is de tijd die het systeem nodig heeft om van chaos naar orde te gaan (de "thermalisatie").

• Ze ontdekten dat dit proces ongeveer 1 "fm/c" duurt. Dat is een ongelofelijk korte tijd (een femtoseconde), maar in de wereld van atomen is dat eeuwig.
• Ze vonden dat dit proces begint net voordat de twee atoomkernen volledig op elkaar liggen, en dat het ongeveer even lang duurt, ongeacht hoe hard je botst.

6. Waarom is dit belangrijk?

Stel je voor dat je een detective bent die een misdaad probeert op te lossen. De "misdaad" is het zoeken naar het kritieke punt in het universum. Als je weet hoe het soepje zich gedraagt (welk recept klopt), kun je zeggen: "Ah, hier is het kritieke punt!"

Deze paper zegt: "We hebben een betere manier gevonden om te kijken naar de deeltjes die eruit vliegen. We hebben de 'knoppen' van onze computermodel (zoals de wachttijd en de dikte van de soep) afgesteld zodat ze precies overeenkomen met de echte metingen."

Conclusie in één zin:
De wetenschappers hebben een betere "tijdmachine" (computermodel) gebouwd om te simuleren hoe atoomkernen botsen; ze ontdekten dat bij lagere energieën de manier waarop het materiaal "opwarmt" cruciaal is, en dat twee verschillende theorieën over hoe materie zich gedraagt, bij lage energieën verschillende voorspellingen doen, wat hen dichter brengt bij het vinden van het mysterieuze kritieke punt in de natuurkunde.

Technische samenvatting

Hier is een gedetailleerde technische samenvatting van het artikel "Particle spectra in the integrated hydrokinetic model at RHIC Beam-Energy-Scan energies", geschreven in het Nederlands.

Titel: Deeltjesspectra in het geïntegreerd hydrokinetisch model bij energieën van de RHIC Beam-Energy-Scan

Auteurs: Narendra Rathod, Yuri Sinyukov, Musfer Adzhymambetov, en Hanna Zbroszczyk.

---

1. Probleemstelling en Context

Relativistische zware-ionenbotsingen creëren quark-gluonplasma (QGP), een toestand van materie die bij extreme temperaturen en dichtheden voorkomt. Hoewel het standaardmodel voor botsingen bij hoge energieën (zoals bij de LHC en de top-energieën van RHIC) goed gevestigd is, wordt het beeld complexer bij lagere energieën (het Beam Energy Scan of BES-programma van RHIC, $\sqrt{s_{NN}} = 7.7 - 39$ GeV).

De belangrijkste uitdagingen bij deze lagere energieën zijn:

• Langere overlaptijd: De tijd waarin de twee kernen elkaar overlappen is aanzienlijk langer (exceedt 1,5 fm/c bij 14,5 GeV) vergeleken met TeV-energieën ($\sim 10^{-3}$ fm/c). Dit leidt tot sterke inhomogeniteiten in zowel transversale als longitudinale richtingen, waardoor de aanname van boost-invariantie (Bjorken-beeld) faalt.
• Thermalisatie: De thermalisatie is trager en het hydrodynamische stadium kan korter zijn of zelfs ontbreken. Het systeem bevindt zich mogelijk langere tijd in een niet-thermisch evenwicht.
• Faseovergang: Het doel van het BES-programma is het vinden van het kritieke punt in het QCD-fasediagram en het identificeren van een eerste-orde faseovergang. Dit vereist een nauwkeurige modellering van de evolutie van het systeem, inclusief de overgang van niet-thermisch evenwicht naar hydrodynamica en terug naar hadronen.

Bestaande hybride modellen die hydrodynamica en microscopische transportmodellen scheiden in distincte ruimtetijd-fasen, zijn mogelijk niet adequaat voor deze energieën. Er is behoefte aan een model dat deze fasen continu integreert.

2. Methodologie: Het uitgebreide Integrated HydroKinetic Model (iHKMe)

De auteurs gebruiken een uitgebreide versie van het Integrated HydroKinetic Model (iHKMe), ontwikkeld om de volledige ruimtetijd-evolutie van het systeem te beschrijven, van de initiële kernoverlap tot de uiteindelijke decoupling van vrij stromende hadronen.

Kerncomponenten van het model:

• Pre-equilibrium en Thermalisatie: In plaats van een scherpe scheiding, integreert iHKMe microscopische (UrQMD cascade) en macroscopische (relativistische hydrodynamica) beschrijvingen gelijktijdig.
  • Het systeem begint ver van evenwicht (beschreven door UrQMD).
  • Tijdens het thermalisatiestadium ($\tau_0$ tot $\tau_{th}$) coëxisteren een niet-equilibriumpart (UrQMD) en een bijna-equilibriumpart (hydrodynamica).
  • Een gewichtsfunctie $P(\tau)$, gebaseerd op de Boltzmann-vergelijking in relaxatietijdsbenadering, regelt de overgang van UrQMD naar hydrodynamica.
• Vrije parameters: Het model bevat vrije parameters die gekalibreerd moeten worden, waaronder:
  • $\tau_0$: Start van thermalisatie.
  • $\tau_{th}$: Einde van thermalisatie.
  • $\tau_{rel}$: Relaxatietijd (snelheid van thermalisatie).
  • $\eta/s$: Verhouding van schuifviscositeit tot entropiedichtheid.
  • $\epsilon_{sw}$: Energie dichtheid waarbij overgang naar de hadronische "afterburner" (UrQMD) plaatsvindt.
  • $R$: Gaussische uitvlekparameter voor initiële condities.
• Equaties van Staat (EoS): Er worden twee verschillende EoS's getest om de gevoeligheid voor de aard van de faseovergang te onderzoeken:
  1. CO (Crossover): Een zachte overgang tussen QGP en hadron-resonantie-gas (HRG).
  2. PT (First-order Phase Transition): Een EoS met een eerste-orde faseovergang en een "zachte" toestand (softening).
• Kalibratie: De auteurs voeren 750 simulaties uit voor Au+Au botsingen bij $\sqrt{s_{NN}} = 14.5$ GeV (20-30% centraliteit) met Latijnse hyperkubus-sampling om de correlaties tussen parameters en observabelen te analyseren. Vervolgens worden de parameters geoptimaliseerd om de experimentele data van de STAR-collaboratie te beschrijven over het volledige BES-energiebereik.

3. Belangrijkste Bijdragen en Resultaten

A. Gevoeligheidsanalyse van parameters

• De initiële gladmakingsparameter ($R$) heeft een sterke invloed op hydrodynamische simulaties, wat de extractie van materiaaleigenschappen beperkt.
• De relaxatietijd ($\tau_{rel}$) gedraagt zich vergelijkbaar met de starttijd ($\tau_0$), maar heeft een minder uitgesproken effect.
• Antiprotonen zijn bijzonder gevoelig voor de schakelenergie-dichtheid ($\epsilon_{sw}$), omdat transportmodellen zoals UrQMD de productie van (anti)baryonen in de late fase onderdrukken.

B. Thermalisatietijd en Evolutie

• De thermalisatie begint kort voor de volledige overlap van de kernen en duurt ongeveer 1 fm/c over het hele BES-energiebereik.
• De starttijd $\tau_0$ schaal met de overlaptijd van de kernen ($\tau_{overlap}$), waarbij $\tau_0 \approx 0.75 \times \tau_{overlap}$.
• Bij lagere energieën (bijv. 7.7 GeV) is het hydrodynamische stadium langer, maar het systeem bereikt toch grotendeels thermisch evenwicht, behalve in de perifere regio's. Bij zeer lage energieën (onder 4.5 GeV) wordt verwacht dat het niet-thermische component dominant wordt.

C. Vergelijking van Equaties van Staat (CO vs. PT)

• Hogere energieën: Bij hogere BES-energieën (boven 14.5 GeV) zijn de verschillen in de transversale impuls-spectra ($p_T$) tussen de crossover- en de eerste-orde faseovergang-scenario's verwaarloosbaar. Beide EoS's beschrijven de data even goed met lichte aanpassingen van de parameters.
• Lagere energieën (7.7 GeV): De grootste verschillen treden op bij de laagste energie ($\sqrt{s_{NN}} = 7.7$ GeV). Hier tonen protonen en kaonen een duidelijke gevoeligheid voor de gekozen EoS en de freeze-out parameters. De verschillen worden toegeschreven aan de verschillende trajecten die het systeem volgt in het $T-\mu_B$ fasediagram.

D. Beschrijving van de Spectra

• Het model beschrijft de spectra van pionen en kaonen goed.
• Er is een systematische afwijking voor protonen en antiprotonen: het model overschat de protonopbrengst en onderschat de antiprotonopbrengst (tot wel 40% bij lage $p_T$). Dit wordt toegeschreven aan een onbalans in baryoncreatie- en annihilatieprocessen tijdens de hadronische afterburner-fase of onnauwkeurigheden in de baryon-dichtheidsverdeling op het particlisatie-hypervlak.

4. Significatie en Conclusie

De studie demonstreert dat het uitgebreide iHKMe-model een robuust kader biedt voor het analyseren van zware-ionenbotsingen in het BES-energiebereik.

• Unificatie: Het model slaagt erin om de overgang van niet-equilibrium naar hydrodynamica en terug naar transport naadloos te modelleren, wat essentieel is voor lagere energieën waar de traditionele scheiding van fasen faalt.
• EoS-onafhankelijkheid: Een cruciale bevinding is dat de transversale impuls-spectra op zich niet voldoende zijn om tussen een crossover en een eerste-orde faseovergang te onderscheiden bij de meeste BES-energieën, tenzij bij de laagste energie (7.7 GeV).
• Toekomstperspectief: Om de aard van de faseovergang en het kritieke punt definitief te bepalen, zijn aanvullende bulk-observabelen nodig, zoals elliptische stroom ($v_2$) en femtoscopie-analyse. De auteurs plannen om deze observabelen in toekomstig werk binnen hetzelfde iHKMe-raamwerk te onderzoeken.

Samenvattend biedt dit werk een gedetailleerde kalibratie van het iHKMe-model en benadrukt het dat, hoewel de huidige spectra goed worden beschreven, de gevoeligheid voor de onderliggende QCD-fysica (zoals de aard van de faseovergang) voornamelijk zichtbaar is bij de laagste BES-energieën en in specifieke deeltjesyields (protonen/kaonen).