Relativistic viscous hydrodynamics in the Density Frame for heavy-ion collisions

In het hart van de meest gewelddadige botsingen in het universum, waar atoomkernen met bijna de snelheid van het licht op elkaar botsen, gedraagt materie zich op manieren die onder normale omstandigheden onmogelijk lijken. Wanneer protonen of loodkernen botsen in faciliteiten zoals de Large Hadron Collider, creëren ze een vluchtige, superhete soep van subatomaire deeltjes die bekend staat als het quark-gluonplasma. Gedurende een fractie van een seconde stroomt dit plasma als een vloeistof, waarbij het uitzet en afkoelt voordat het bevriest in de deeltjes die detectoren uiteindelijk registreren. Om dit proces te begrijpen, vertrouwen natuurkundigen op een reeks wiskundige regels genaamd hydrodynamica, die beschrijven hoe vloeistoffen bewegen en veranderen. Het toepassen van deze regels op de extreme, relativistische snelheden van deze botsingen is echter lang een uitdaging geweest. De standaardmethoden die worden gebruikt om deze vloeistof te modelleren, vereisen vaak het toevoegen van extra, onzichtbare variabelen om de wiskunde stabiel te houden en te voorkomen dat de vergelijkingen instorten. Deze extra variabelen fungeren als een vangnet, maar ze introduceren een laag van complexiteit waardoor het moeilijk wordt om precies te weten welke eigenschappen van de vloeistof worden gemeten en welke slechts artefacten zijn van het wiskundige model.

Een team van onderzoekers heeft nu een nieuwe manier ontwikkeld om deze simulaties uit te voeren, waarbij ze de onnodige veiligheidsnetten weghalen om een schoner, directer model van het gedrag van de vloeistof te creëren. In een recente studie hebben Jean-François Paquet, Mayank Singh en Derek Teaney een methode geïntroduceerd die zij de "Density Frame" noemen. In plaats van de stroming van de vloeistof en de interne spanningen als afzonderlijke, evoluerende zaken bij te houden, behandelt deze nieuwe aanpak de energie en het momentum van de vloeistof als de enige zaken die in de loop van de tijd veranderen. Alles wat volgt, inclusief de interne wrijving die de stroming tegenwerkt, wordt direct berekend op basis van hoe de vloeistof op dat exacte moment beweegt. Dit elimineert de noodzaak voor extra, niet-fysische variabelen die eerdere modellen hebben vertroebeld. De onderzoekers testten deze nieuwe motor door deze te vergelijken met bekende oplossingen voor hoe deeltjes zich in een vacuüm gedragen en door deze naast de meest gebruikte simulatiesoftware in het vakgebied te draaien, bekend als MUSIC. Ze ontdekten dat de nieuwe methode kleine, scherpe rimpelingen in de vloeistof sneller gladstrijkt dan de oude methode, maar dat de uiteindelijke resultaten voor de deeltjes die uit de botsing voortkomen, opmerkelijk vergelijkbaar zijn. Dit suggereert dat de complexe extra stappen in de oude modellen mogelijk niet nodig zijn om het grote plaatje van de stroming van het quark-gluonplasma te begrijpen.

De kern van dit werk ligt in de manier waarop de onderzoekers het concept van entropie, of wanorde, binnen de vloeistof behandelen. In elke echte vloeistof creëert wrijving warmte en verhoogt het de wanorde, een proces dat continu moet plaatsvinden. De nieuwe methode is gebouwd rond een wiskundig instrument genaamd een "dissipatieve functie", die fungeert als een grootboek voor deze wanorde. Door hun computercode zo te ontwerpen dat elke stap in de simulatie deze wanorde automatisch op een consistente manier verhoogt, hebben de onderzoekers ervoor gezorgd dat hun model stabiel blijft, zelfs wanneer de vloeistof zeer dun of turbulent wordt. Deze stabiliteit maakte het mogelijk om de botsing te simuleren zonder de kunstmatige "regulatoren" te gebruiken die andere codes vereisen om te voorkomen dat de wiskunde crasht. In hun tests begonnen ze met een specifiek type initiële conditie genaamd "IP-Glasma", die de chaotische staat van de kernen vlak voordat ze botsen nabootst. Ze lieten de simulatie vooruit in de tijd lopen en observeerden hoe de temperatuur en de stroomsnelheid evolueerden.

Wanneer ze de resultaten van de nieuwe Density Frame vergeleken met de standaard MUSIC-code, merkten ze een duidelijk verschil op in de manier waarop de twee de interne wrijving van de vloeistof afhandelden. In de standaardcode heeft de weerstand van de vloeistof tegen stroming, bekend als schuifspanning (shear stress), een klein beetje tijd nodig om zich aan te passen aan veranderingen in de stroming, vergelijkbaar met een zwaar object dat een moment nodig heeft om in beweging te komen of te stoppen. In de nieuwe Density Frame past deze weerstand zich onmiddellijk aan. Dit verschil was nog prominenter in de bulkviscositeit, oftewel de weerstand van de vloeistof om samengedrukt of uitgerekt te worden. In het standaardmodel past deze weerstand zich langzaam aan over de tijd, vooral nabij het punt waar de vloeistof overgaat van een heet plasma naar een koeler gas. In het nieuwe model reageert het onmiddellijk. Ondanks deze verschillen in het midden van het proces, was de uiteindelijke uitkomst verrassend consistent. Wanneer de onderzoekers keken naar de deeltjes die aan het einde van de botsing werden geproduceerd, zoals het totale aantal gecreëerde deeltjes, hun gemiddelde snelheid en hoe ze in verschillende richtingen verspreidden, produceerden de twee methoden bijna identieke resultaten.

De studie onderzocht ook wat er gebeurt als de vloeistof geen weerstand heeft tegen het worden samengedrukt, een scenario waarbij de bulkviscositeit op nul wordt gezet. In dat geval kwamen de twee methoden bijna perfect overeen. Wanneer ze een temperatuurafhankelijke bulkviscositeit toevoegden, die het complexe gedrag van de vloeistof nabij zijn overgangspunt nabootst, verschenen er kleine verschillen, met name in de buitenste randen van de botsingszone. Deze verschillen waren echter over het algemeen kleiner dan het effect van de bulkviscositeit zelf. De onderzoekers concludeerden dat de complexe, niet-hydrodynamische modi waar de oude modellen op vertrouwen om stabiel te blijven, blijkbaar geen belangrijke rol spelen bij het bepalen van de uiteindelijke waarneembare eigenschappen van de botsing. Deze bevinding is significant omdat het suggereert dat de huidige onzekerheid bij het meten van de eigenschappen van de vloeistof mogelijk niet komt door de fysica van de vloeistof zelf, maar door de extra wiskundige bagage die de simulatietools met zich meedragen.

Door deze extra variabelen te verwijderen, biedt de Density Frame een directer pad naar het begrijpen van de fundamentele eigenschappen van het quark-gluonplasma. De onderzoekers hebben hun code openbaar gemaakt en gedocumenteerd, waardoor andere wetenschappers deze kunnen gebruiken als vervanging voor de hydrodynamische motoren in bestaande simulatiepakketten. Dit opent de deur naar een nieuw soort analyse waarbij de enige parameters die worden aangepast de werkelijke fysische eigenschappen van de vloeistof zijn, zoals de viscositeit en de toestandsvergelijking, in plaats van de willekeurige tijdsvertragingen die in de wiskundige modellen zijn ingebouwd. De auteurs suggereren dat deze schonere aanpak bijzonder nuttig zal zijn voor het bestuderen van kleinere botsingssystemen, waarbij de vloeistof minder ontwikkeld is en de invloed van deze extra wiskundige modi prominenter kan zijn. Uiteindelijk biedt dit werk een robuustere en transparantere tool voor natuurkundigen om de ware transporteigenschappen van het quark-gluonplasma te extraheren uit de gegevens verzameld bij hoogenergetische versnellers, wat ons dichter brengt bij een volledig begrip van de meest extreme vloeistof in het universum.