Fast partial-sky spherical harmonic transforms

Titel: Snelle partiële-hemel sferische harmonische transformaties

Auteurs: Julien Carron en Martin Reinecke

1. Het Probleem

Sferische harmonische transformaties (SHT) zijn fundamenteel voor de analyse van data die op een boloppervlak zijn verdeeld, zoals de Kosmische Microgolf-Achtergrondstraling (CMB). De twee kernoperaties zijn:

1. Synthese: Het reconstrueren van een kaart op het boloppervlak vanuit harmonische coëfficiënten ($g_{\ell m}, c_{\ell m}$).
2. Adjoint Synthese: Het berekenen van harmonische coëfficiënten vanuit waarnemingen op het boloppervlak.

De uitdaging waar dit artikel op inzoomt, is het efficiënt uitvoeren van deze transformaties voor kleine, willekeurig gepositioneerde gebieden op de hemel (partiële hemel), zonder de beperkingen van traditionele methoden.

• Bestaande methoden:
  • Isolatitude pixelisatie (bijv. HEALPix): Vereist dat data op specifieke breedtegraden ligt. Dit is inefficiënt als data willekeurig is of als men alleen een klein gebied analyseert, omdat de hele bol moet worden verwerkt.
  • Vlakke-hemelse benadering (Flat-sky): Werkt goed voor zeer kleine gebieden, maar faalt bij grotere schalen of wanneer nauwkeurige bol-geometrie nodig is (bijv. voor lensing-reconstructie).
• Motivatie: De analyse van diepe polarisatie-data van de South Pole Telescope (SPT-3G) en het BICEP3-programma. Deze projecten vereisen iteratieve methoden voor CMB-lensing en het zoeken naar oorspronkelijke gravitatiegolven, wat duizenden transformaties op kleine, niet-symmetrische hemelgebieden vereist.

2. Methodologie

De auteurs bouwen voort op de Fourier-sfeer-methode van Reinecke et al. (2023), maar introduceren een nieuwe strategie voor kleine gebieden die de effectieve bandbreedte van de transformatie verlaagt.

Kernconcepten:

1. Coördinaten-herdefinitie en "Stretching":
   • Het te analyseren gebied wordt zo gepositioneerd dat het centrum op de noordpool ligt.
   • Het kleine gebied wordt "uitgerekt" (geschaald) zodat het een groter deel van de Fourier-sfeer beslaat. Hierdoor worden de fysieke kenmerken in het gebied groter, wat effectief de benodigde bandbreedte ($\ell_{max}$) verlaagt.
2. Apodisatie en Vensterfuncties:
   • Omdat het uitrekken de periodiciteit en gladheid van de functie verstoort, wordt een gewichtsfunctie (venster) $w(\Theta)$ toegepast.
   • Deze functie is 1 binnen het gebied van belang en daalt glad naar 0 aan de randen van het uitgebreide gebied.
   • Vensterkeuze: De auteurs gebruiken een "exponentiële van een halve cirkel" (ES-kern) of Slepian-functies (prolate spheroidal wave functions). Deze functies zijn optimaal omdat ze compact zijn in de ruimtelijke domein en een zeer snelle afname hebben in het frequentiedomein, wat de spectrale lekkage minimaliseert.
3. Bandbreedte-verlaging ($L_{max}$):
   • Door het gebied te verkleinen en te schalen, is het niet meer nodig om $2\ell_{max} + 1$ breedtegraden te berekenen. In plaats daarvan volstaat een kleiner aantal rings ($L_{max} + 1$) binnen het nieuwe, schaalbare domein.
   • $L_{max}$ wordt bepaald door de oorspronkelijke resolutie, de grootte van het gebied en de gewenste tolerantie voor fouten.
4. Polaire Optimalisatie:
   • Omdat het centrum van het gebied op de pool ligt, worden de breedtegraden (rings) dichter bij de pool korter (schaalbaar met $\sin \theta$).
   • Dit maakt het mogelijk om hoogfrequente longitudinale golven ($m$-waarden) die groter zijn dan $\ell_{max} \sin \theta$ te negeren. Dit vermindert de rekentijd aanzienlijk.

Implementatie:
De methode vereist slechts minimale aanpassingen aan de bestaande, geoptimaliseerde ducc0.sht code. De kernveranderingen zijn de herschaling van coördinaten en het toepassen van de gewichtsfunctie op de Fourier-rings.

3. Belangrijkste Bijdragen

• Nieuw algoritme: Een efficiënte methode voor SHT op willekeurige, kleine hemelgebieden die geen isolatitude pixelisatie vereist.
• Foutcontrole: Een rigoureuze analyse van de fouten die ontstaan door het uitrekken en apodiseren. De auteurs tonen aan dat de fouten eenvoudig kunnen worden beheerst door de keuze van het venster en de apodisatielengte, en dat ze onder een door de gebruiker gespecificeerde tolerantie kunnen worden gehouden.
• Eliminatie van compromissen: De methode combineert de voordelen van exacte bol-geometrie (geen vlakke-hemelse benadering) met de snelheid van lokale berekeningen, zonder de nadelen van willekeurige pixelisatie.

4. Resultaten

De auteurs hebben de prestaties getest op verschillende configuraties, variërend van grote gebieden (SPT-3G) tot zeer kleine gebieden (EDFS).

• Snelheidswinst:
  • SPT-3G Main Field (~1500 deg²): Snelheidswinst van factor 3 tot 4.
  • BICEP3 Footprint (~980 deg²): Snelheidswinst van factor 6 tot 9.
  • EDFS (Euclid Deep Field South, ~57 deg²): Snelheidswinst van factor 21 tot 43 (tot wel 85x voor bepaalde adjoint transformaties).
• Vergelijking met bestaande methoden:
  • De methode is sneller dan het gebruik van traditionele HEALPix-rings (die op isolatitude liggen) voor dezelfde data, zelfs als de data oorspronkelijk op die rings is gegenereerd.
  • Voor zeer hoge resoluties op kleine vlakke gebieden is de methode vergelijkbaar met of sneller dan een 2D FFT, maar dan met behoud van de exacte bol-geometrie.
• Fouten: De fouten zijn voorspelbaar en kunnen worden verkleind door de tolerantie ($\epsilon$) te verlagen, wat leidt tot een kleine toename in rekentijd.

5. Betekenis en Toekomstperspectief

• CMB-wetenschap: De methode maakt complexe, iteratieve CMB-lensing-reconstructies en zoektochten naar B-modi (voor oorspronkelijke gravitatiegolven) haalbaar op standaard hardware, zonder dat er zware High-Performance Computing (HPC) resources nodig zijn.
• Algemene toepasbaarheid: Het algoritme is ideaal voor elke toepassing die frequente SHT's vereist op beperkte hemelgebieden, zoals catalogus-analyses van grote structuren of cluster-lensing.
• Software: De implementatie is beschikbaar via de lenspyx bibliotheek en is een lichte modificatie van de bestaande ducc code, wat de adoptie voor andere onderzoeksgroepen vergemakkelijkt.

Conclusie:
Dit artikel presenteert een doorbraak in de numerieke efficiëntie van sferische harmonische transformaties voor partiële hemelgebieden. Door slim gebruik te maken van coördinaten-schaling, apodisatie en polaire optimalisatie, elimineert het de noodzaak voor compromissen tussen nauwkeurigheid en snelheid, wat essentieel is voor de volgende generatie kosmologische experimenten.