Non-Gaussian timing residual skewness in pulsar timing arrays

Diep in het stille ritme van het universum vindt een nieuwe vorm van luisteren plaats. Decennialang hebben astronomen telescopen gebruikt om licht te vangen, maar een ander instrument is opgekomen om naar het weefsel van de ruimte zelf te luisteren. Dit instrument is de pulsar timing array, een kosmisch netwerk dat de hele melkweg als detector behandelt. Het steunt op pulsars, de in elkaar geklapte, draaiende kernen van dode sterren. Deze stellaire resten zijn de meest precieze klokken van de natuur; ze draaien met een zodanige regelmaat dat hun pulsen de aarde bereiken met de consistentie van een metronoom. Wanneer een zwaartekrachtgolf — een rimpeling in de ruimtetijd veroorzaakt door massieve kosmische gebeurtenissen — tussen de aarde en een pulsar passeert, rekt en krimpt zij de ruimte daartussen uit. Deze minuscule vervorming verandert de tijd die een puls nodig heeft om aan te komen, wat een lichte vertraging of "timing residual" creëert. Door een collectie van deze pulsars te monitoren, kunnen wetenschappers het zwakke, collectieve gezoem van een stochastische zwaartekrachtgolfachtergrond detecteren, een zee van overlappende rimpelingen van talloze bronnen doorheen de kosmos.

De vraag die veel van dit onderzoek drijft, is de oorsprong van dat gezoem. Is het het koor van supermassieve zwarte gaten die om elkaar heen draaien in verre sterrenstelsels, of is het een gefluister van het allereerste begin van de tijd, een overblijfsel van de snelle expansie van het universum? Hoewel huidige metingen het bestaan van deze achtergrond hebben bevestigd, hebben ze nog niet onderscheid kunnen maken tussen deze mogelijkheden. De standaardmethode van detectie kijkt naar hoe de timing van verschillende pulsars met elkaar correleert, een patroon dat bekend staat als de Hellings-Downs curve. Deze twee-punts correlatie is echter als het luisteren naar het gemiddelde volume van een menigte; het vertelt je dat de menigte er is, maar niet wie er spreekt of wat ze zeggen. Om de bron te vinden, moeten wetenschappers zoeken naar complexere patronen, specifiek zoekend naar "niet-Gaussianiteit". In eenvoudige bewoordingen is een Gaussianistisch signaal een signaal waarbij de fluctuaties perfect willekeurig en symmetrisch zijn, zoals de ruis op een radio. Als het signaal afwijkt van deze perfecte willekeur en een specifieke soort asymmetrie of "skew" vertoont, zou dit de unieke vingerafdruk van de oorsprong kunnen onthullen, of dat nu de chaotische dans van zwarte gaten is of de fysica van het vroege universum.

Een team van onderzoekers begon onlangs met het onderzoek of een dergelijke afwijking gedetecteerd kon worden met behulp van deze pulsarnetwerken. Ze richtten zich op een specif kind van asymmetrie genaamd skewness (scheefheid), die meet of de timing residuals meer naar de late kant of juist naar de vroege kant neigen. Om dit te doen, stelden ze een wiskundig model voor waarbij de zwaartekrachtgolven grotendeels willekeurig zijn, maar een kleine, kwadratische correctie bevatten. Deze correctie vertegenwoordigt een zwakke interactie tussen verschillende delen van de golf, een scenario dat wordt gemotiveerd door theorieën over het vroege universum, maar ook van toepassing is op astrofysische bronnen. De onderzoekers probeerden dit model eerst op te bouwen met behulp van de fundamentele eigenschappen van de golven zelf, specifiek hun "heliciteit", die de richting van hun spin beschrijft. Ze berekenden hoe de timing residuals eruit zouden zien als dit model waar zou zijn. Echter, ze ontdekten een cruciaal gebrek: het resultaat hing af van hoe de wetenschappers ervoor kozen de oriëntatie van de spin van de golven te definiëren. Aangezien de fysieke realiteit van het universum niet verandert op basis van de willekeurige keuze van een waarnemer voor een coördinatensysteem, werd deze initiële aanpak als inconsistent beschouwd. Het model produceerde een signaal dat geen echte fysieke observeerbare was, omdat het verschoof wanneer het referentiekader werd geroteerd.

Onvermoeibaar hervormde het team hun aanpak. In plaats van het model op te bouwen vanuit de abstracte eigenschappen van de golven, construeerden ze het direct vanuit het signaal dat de detectoren daadwerkelijk meten. Ze definieerden een "detector-geprojecteerde" amplitude, die de eigenschappen van de golf combineert met de specifieke geometrie van de gezichtslijn van de pulsar. Deze nieuwe formulering zorgde ervoor dat de fasen van de golven correct wegvielen, waardoor het resultaat onafhankelijk werd van een willekeurige keuze van oriëntatie. Met dit gecorrigeerde model berekenden ze de verwachte skewness voor een enkele pulsar. Ze vereenvoudigden de berekening door zich te concentreren op het deel van het signaal dat afkomstig is van de positie van de aarde, waarbij ze de complexere en meer afstandelijke bijdrage van de pulsar zelf negeerden, en namen aan dat de achtergrond van golven in alle richtingen uniform was. Het resultaat van hun numerieke benchmarks was een zeer kleine waarde. Wanneer ze deze skewness normaliseerden tegen de typische grootte van de willekeurige fluctuaties, bleek het signaal extreem klein te zijn.

De studie concludeert dat hoewel een niet-Gaussianistisch signaal van dit specifieke type theoretisch mogelijk en wiskundig consistent is wanneer het gedefinieerd wordt vanuit het perspectief van de detector, de resulterende asymmetrie in de timing residuals waarschijnlijk te klein is om gemakkelijk van de achtergrondruis te onderscheiden met de huidige methoden. De onderzoekers hebben niet gevonden hoe ze het bestaan van dergelijke signalen zouden kunnen uitsluiten, maar hun werk laat zien dat de specifieke soort niet-Gaussianiteit die zij modelleerden, een skewness produceert die in de praktijk verwaarloosbaar is. Deze bevinding dient als een cruciale benchmark voor toekomstige zoektochten. Het vertelt de wetenschappelijke gemeenschap dat als zij op zoek gaan naar dit specifieke soort asymmetrie om tussen kosmologische en astrofysische bronnen te onderscheiden, zij voorbereid moeten zijn op een signaal dat ongelooflijk zwak is. Het werk verheldert het theoretische landschap door aan te tonen dat hoewel het idee om skewness te meten solide is, het specifieke model dat hier getest is een resultaat oplevert dat in de praktijk effectief onzichtbaar is tegen de standaard willekeurige fluctuaties van de zwaartekrachtgolfachtergrond.