Cylindrical cosmological simulations with StePS

Technische Samenvatting: Cilindrische Kosmologische Simulaties met StePS

Probleemstelling
Standaard kosmologische $N$-body simulaties evolueren materie doorgaans binnen een periodieke kubusvormige box, wat een drie-torus ($T^3 = S^1 \times S^1 \times S^1$) topologie oplegt. Hoewel deze topologie numeriek handig is en statistische homogeniteit behoudt, verbreekt het de statistische isotropie omdat de kubische geometrie alleen invariant is onder de octaëdrische groep ($O_h$), en niet onder de continue rotatiegroep $SO(3)$. Bij een gevolg is het impulsmoment in periodieke boxen over het algemeen niet behouden, en kan de kubische symmetrie subtiele anisotropieën in het zwaartekrachtsveld introduceren, wat de grootschalige clusteringstatistieken potentieel kan beïnvloeden, vooral wanneer het simulatievolume niet aanzienlijk groter is dan de schalen van belang. Bovendien zijn standaard periodieke grenzen minder geschikt voor het simuleren van systemen met inherente cilindrische symmetrieën, zoals kosmische filamenten of specifieke anisotrope kosmologische modellen.

Methodologie
De auteurs introduceren een gecompacteerd simulatieframework gebaseerd op de $S^1 \times \mathbb{R}^2$ (slab) topologische manifold. In deze geometrie is het simulatiedomein periodiek langs één as (de $z$-as, met periode $L_z$) maar open in de radiale richting ($\mathbb{R}^2$).

• Geometrische Compactificatie: De auteurs maken gebruik van een inverse stereografische projectie langs de radiale richting om het niet-compacte $\mathbb{R}^2$-vlak te compactificeren. Dit brengt de oneindige radiale coördinaat in kaart naar een eindig domein, waardoor een "zoom-in" strategie mogelijk is waarbij de ruimtelijke en massale resolutie hoog is nabij de centrale as en geleidelijk afneemt naar de buitenste grens toe. Dit onderdrukt randartefacten bij de vrije buitenrand.
• Bewegingsvergelijkingen: De Newtoniaanse bewegingsvergelijkingen worden opgelost in comoving coördinaten. De berekening van de zwaartekracht wordt aangepast aan de $S^1 \times \mathbb{R}^2$ topologie door de sommatie over periodieke beelden langs de $z$-as te verrichten. Om de trage convergentie van directe sommatie te beheersen, implementeren de auteurs een Ewald-sommatie methode die is aangepast voor cilindrische geometrie, waarbij de kracht wordt opgesplitst in real-space en reciproke-space componenten.
• Krachtberekeningsalgoritmen: De implementatie binnen het StePS-framework (STEreographically Projected cosmological Simulations) ondersteunt twee methoden voor krachtberekening:
  1. Directe $O(N^2)$: Brengt efficiënt naar GPU's, geschikt voor kleinere deeltjesaantallen of kleinere cilindrische volumes.
  2. Octree $O(N \log N)$: Geïmplementeerd voor CPU-clusters. Om systematische fouten als gevolg van de kubische aard van het octree-grid in een cilindrisch domein te beperken, past de code bij elke tijdstap willekeurige verschuivingen en rotaties van het domein toe. Een cruciale kalibratiestap corrigeert voor radiale krachtfouten om spontane radiale ineenstorting of expansie te voorkomen, wat de stabiliteit van de achtergronddichtheid waarborgt.
• Begincondities: De begincondities worden gegenereerd met behulp van een "glass" deeltjesverdeling, gecreëerd via repulsieve zwaartekracht in de $S^1 \times \mathbb{R}^2$ topologie. Perturbaties worden geïmprinted met behulp van de Zel'dovich-benadering of tweede-orde Lagrangiaanse perturbatietheorie (2LPT), aangepast aan de cilindrische geometrie en de niet-uniforme bemonstering van de zoom-in setup. De vermogensspectra worden geschat met behulp van de Feldman–Kaiser–Peacock (FKP) estimator om rekening te houden met de niet-triviale windowfunctie.

Belangrijkste Bijdragen

1. Nieuwe Topologie Implementatie: Het artikel presenteert de eerste implementatie van kosmologische $N$-body simulaties in de $S^1 \times \mathbb{R}^2$ manifold, waarmee de StePS-code wordt uitgebreid voorbij de eerdere $T^3$ en sferische $\mathbb{R}^3$ mogelijkheden.
2. Conserveringswetten: De auteurs tonen aan dat in deze topologie de Lagrangiaan invariant is onder translatie langs de cilinderas en rotatie in het transversale vlak ($SO(2)$). Bijgevolg zijn het lineaire momentum langs de $z$-as en het impulsmoment langs dezelfde as behouden, wat de problemen met symmetriebreking van de $T^3$ topologie aanpakt.
3. Algoritmische Aanpassingen: Het werk beschrijft de noodzakelijke wijzigingen voor krachtberekening (cilindrische Ewald-sommatie), het genereren van begincondities (cilindrische glass) en halo-detectie (een aangepaste spherical-overdensity finder voor $S^1$ periodiciteit en radiale resolutievariatie).
4. Validatie: Een volledige $\Lambda$CDM-simulatie is uitgevoerd om de methode te valideren, waarbij de vorming van een kosmisch web wordt aangetoond dat consistent is met standaard kosmologische verwachtingen.

Resultaten

• Simulatie-executie: Er is een simulatie gedraaid met $2,4 \times 10^7$ deeltjes in een cilinder met hoogte $L_z = 1,0$ Gpc en straal $R_{sim} = 500$ Mpc, gebruikmakend van Planck 2018 $\Lambda$CDM-parameters.
• Structuurvorming: De simulatie reproduceert succesvol het niet-lineaire kosmische web (voids, muren, filamenten en halos). De zoom-in strategie maakt een hoge-getrouwheid resolutie van laag-massa halos nabij de centrale as mogelijk, terwijl grotere volumes meer massieve halos op grotere radii oplossen.
• Vermogensspectra: De geschatte materie-vermogensspectra komen goed overeen met lineaire CAMB-voorspellingen op grote schalen en de niet-lineaire Mead et al. (2021) "halofit" emulator op kleine schalen, wat de nauwkeurigheid van de methode bij het vastleggen van structuurgroei bevestigt.
• Halo-eigenschappen: De auteurs identificeerden $\sim 1,1 \times 10^5$ donkere materie halos bij $z=0$ en berekenden 68 parameters voor elk, waaronder massa, concentratie en impulsmoment.
• Anisotropietest: Gemotiveerd door recente claims over de anisotropie van de spin van sterrenstelsels, hebben de auteurs de verdeling van halo-spinassen ten opzichte van de periodieke $z$-as getest. Met behulp van Kolmogorov-Smirnov-testen vonden zij dat voor $L_z = 1,0$ Gpc de spinverdeling consistent is met statistische isotropie. Zij concluderen dat deze topologie op deze schaal de gerapporteerde asymmetrieën in de spins van sterrenstelsels niet van nature verklaart.

Betekenis en Claims
Het artikel beweert dat de $S^1 \times \mathbb{R}^2$ topologie een gespecialiseerd alternatief biedt voor standaard kubische simulaties, met name voor het bestuderen van systemen met cilindrische symmetrieën (bijv. filamenten) of modellen waarbij de conservering van het impulsmoment langs een specifieke as fysiek relevant is. De auteurs benadrukken dat hoewel de methode onderscheidende systematische effecten introduceert — zoals een radiaal variërende resolutie en selectie-effecten — het een zelfconsistent framework biedt dat vrij is van de kubische anisotropie van $T^3$ boxen.

De auteurs zijn bescheiden in hun claims en stellen dat zij niet verwachten dat dit algoritme algemene kosmologische $T^3$ simulaties zal vervangen. In plaats daarvan dient het als een complementair instrument voor specifieke wetenschappelijke vragen waarbij kubische symmetrie een beperkende factor is. Het werk benadrukt dat het anisotrope krachtveld van de topologie artefacten kan introduceren als de periodieke lengte te klein is ($L_z \lesssim 200$ Mpc), een onderwerp dat gereserveerd is voor toekomstige gedetailleerde karakterisering. De code wordt als open-source vrijgegeven (StePS v2) om verdere exploratie van deze geometrieën te faciliteren.