Impact of Spacecraft Orbit Uncertainties and Velocity Mismodeling on the LISA Gravitational-Wave Response

Technische Samenvatting: Impact van Onzekerheden in de Baan van Ruimtevaartuigen en Foutieve Snelheidsmodellering op de Gravitatiegolfrespons van LISA

Probleemstelling
De Laser Interferometer Space Antenna (LISA) is ontworpen om gravitatiegolven (GW's) te detecteren in het millihertz-frequentiebereik met behulp van een constellatie van drie ruimtevaartuigen. Nauwkeurige analyse van GW-gegevens vereist een precieze kennis van de responsfunctie van de detector, die de inkomende GW-rek (strain) koppelt aan het gemeten fractionele frequentiesignaal. Deze respons is kritisch afhankelijk van de posities, snelheden en armlengtes van de ruimtevaartuigen. Hoewel eerdere studies de instrumentele kalibratie en trajectsimulatie hebben behandeld, bleven twee specifieke vragen over de getrouwheid van het responsmodel onbeantwoord:

1. Hoe propageren onzekerheden in de kennis van de banen en armlengtes van de ruimtevaartuigen in de GW-respons?
2. Wat is de impact van het verwaarlozen van de snelheidstermen van de ruimtevaartuigen (Doppler-boosting) in het responsmodel, in het bijzonder voor Galactische binaire bronnen?

Methodologie
De auteurs kwantificeren deze fouten met behulp van de "mismatch"-metriek, die de fractionele verlies in signaal-ruisverhouding (SNR) meet wanneer een incorrecte template wordt gebruikt. De studie maakt gebruik van twee complementaire orbitale scenario's:

1. Statisch Speelgoedmodel: Een vereenvoudigde, tijdonafhankelijke driehoekige configuratie wordt gebruikt om geometrische effecten te isoleren. Drie klassen van perturbaties worden sequentieel toegepast:

   • Armlengtevariatie: Onafhankelijke Gaussische perturbaties ($\sigma_L = 1$ m) om hoog-precieze onboard interferometrische afstandsbepaling te simuleren.
   • Rigide-lichaam rotatie: Willekeurige rotaties die positionele onzekerheden induceren gelijk aan $\sigma_{rot} = 50$ km.
   • Rigide-lichaam translatie: Globale translaties met $\sigma_{trans} = 50$ km om absolute positie-onzekerheden te simuleren vanuit grondgebaseerde radiometrische tracking.

2. Realistisch Evoluerend Scenario: Een hoog-getrouwe simulatie met de ESA GODOT- en MIDAS-softwarepakketten. Dit omvat:

   • Het genereren van een "waarheidstraject" met realistische perturbaties (zelfgravitatie, niet-gravitatieve krachten, klokruis).
   • Het uitvoeren van een uitgebreide Square-Root Information Filter (SRIF) om het traject te schatten op basis van gesimuleerde observaties (afstand, Doppler, pseudodistances).
   • Het bemonsteren van de resulterende baan-covariantie-matrix om een ensemble van traject-schattingen te genereren.
   • Het vergelijken van de LISA-respons berekend met deze geperturbeerde trajecten tegen een nominale respons.

3. Respons Mismodellering: De studie vergelijkt de volledige Doppler-gebooste respons (inclusclusief snelheidstermen $\dot{r}/c$) met een nulde-orde benadering waarbij de snelheden van de ruimtevaartuigen worden verwaarloosd. Dit wordt geëvalueerd met zowel een enkelvoudige frequentie monochromatische bron als een snelle Galactische binaire waveform model (jaxgb).

4. Parameterestimatie: Een Markov Chain Monte Carlo (MCMC) analyse wordt uitgevoerd op een Galactisch binair signaal ($f = 10^{-4}$ Hz, SNR = 200) die gedurende één jaar wordt waargenomen om de parameter-biases te beoordelen die voortvloeien uit de snelheidsmismodellering.

Belangrijkste Resultaten

• Baanonzekerheden:

  • In het realistische evoluerende scenario produceren onzekerheden in de baanbepaling verwaarloosbare mismatches over het gehele LISA-band. De slechtste mismatch is lager dan $10^{-7}$.
  • Fouten worden gedomineerd door hoogfrequente effecten ($>10$ mHz) waar de respons het meest gevoelig is voor de absolute posities van de ruimtevaartuigen.
  • Vergeleken met het statische speelgoedmodel laten de realistische simulaties zien dat armlengtefouten 3–6 keer kleiner zijn, positiefouten 8–20 keer kleiner en hoekfouten ~300 keer kleiner zijn.

• Snelheidsmismodellering:

  • Het verwaarlozen van de snelheden van de ruimtevaartuigen leidt tot aanzienlijk grotere fouten, met name bij lage frequenties.
  • Bij $\sim 10^{-4}$ Hz zijn de mismatches van de orde $10^{-4}$, met relatieve amplitude-fouten rond $10^{-2}$ en fasefouten die een radian bereiken.
  • De fout schaalt ongeveer als $f^{-1}$ bij lage frequenties, en neemt af naar $\sim 10^{-9}$ bij 100 mHz.
  • Snelheidsmismodellering domineert de foutenbudget met ongeveer twee ordes van grootte bij 10 mHz en negen ordes van grootte bij 0,1 mHz vergeleken met baanonzekerheden.

• Parameterestimatie:

  • Voor een Galactische binaire met $f = 10^{-4}$ Hz en SNR = 200, leidt het verwaarlozen van de snelheidstermen in de respons tot parameter-biases van minder dan 1$\sigma$.
  • Het verschil in de gemiddelde log-likelihood tussen de volledige respons en het snelheid-verwaarlozende model is $\Delta \log \mathcal{L} \sim 10$.

Betekenis en Claims
Dit werk biedt de eerste karakterisering van hoe fouten in de LISA-gravitatiegolfrespons propageren van GW-rek door de detectoroutput naar geschatte parameters. De auteurs concluderen dat:

• Baaninzicht: De huidige mogelijkheden voor baanbepaling (met onzekerheden van $\sim 50$ km) zijn voldoende, aangezien de resulterende mismatches te klein zijn om de wetenschappelijke exploitatie van LISA te beïnvloeden.
• Snelheidstermen: Het verwaarlozen van de snelheden van de ruimtevaartuigen is de belangrijkste bron van modelleringsfout voor Galactische binaire signalen, met name bij lage frequenties. Zelfs in het slechtste geteste scenario (SNR=200) blijven de resulterende parameter-biases echter onder 1$\sigma$.
• Aanbeveling: Voor hoog-getrouwe inferentie van laagfrequente Galactische binaire systemen moet de volledige Doppler-gebooste respons worden gehanteerd, aangezien dit geen extra computationele kosten met zich meebrengt vergeleken met de nulde-orde benadering.
• Toekomstig Werk: De auteurs suggereren dat de frequentieafhankelijke nulpunten van de respons kunnen dienen als een kalibratiediagnostiek, en dat toekomstige studies de analyse moeten uitbreiden naar Extreme Mass Ratio Inspirals (EMRIs), stellaire zwarte gat-binairsystemen en stochastische achtergronden.