A numerical framework for Newtonian-noise estimation at the Einstein Telescope: 2-D simulations beyond the plane-wave approximation

Dit artikel introduceert een numeriek framework op basis van spectrale-elementen-simulaties voor de schatting van Newtoniaanse ruis bij de Einstein Telescope, waarmee de beperkingen van de traditionele vlakke-gbenbenadering worden overwonnen en aantoont dat in een homogene omgeving het P-golf-aandeel lager is dan vaak wordt aangenomen, wat betere mogelijkheden voor ruisreductie suggereert.

De kern

De Einstein Telescoop en de "Onzichtbare Ruis": Een Simpele Uitleg

Stel je voor dat je probeert een fluisterend gesprek te horen in een drukke treinwagon. Dat is wat de Einstein Telescoop (ET) doet, maar dan met het heelal. Het is een gigantisch ondergronds apparaat dat bedoeld is om de allerzwakste trillingen van het heelal te horen: zwaartekrachtgolven. Deze golven komen van enorme gebeurtenissen, zoals botsende zwarte gaten.

Het probleem? De aarde zelf is niet stil. De grond trilt voortdurend door wind, verkeer, en zelfs door de zee. Deze trillingen maken een ruis die het signaal van het heelal kan verdoezelen. Dit noemen wetenschappers Newtoniaanse ruis.

In dit artikel leggen Patrick Schillings en zijn team uit hoe ze een nieuwe manier hebben bedacht om deze ruis beter te begrijpen en misschien zelfs te filteren.

1. Het Probleem: De Aarde is geen Perfecte IJsblok

Vroeger dachten wetenschappers dat ze de ruis konden voorspellen door te doen alsof de aarde een perfect, egaal stuk ijs is. Ze gebruikten simpele formules die aannamen dat de trillingen als perfecte, rechte golven door de grond gaan (zoals een laserstraal).

Maar in het echt is de aarde niet als ijs. Het is meer als een grote, onregelmatige cake met verschillende lagen, stukken steen en holtes. Wanneer trillingen door zo'n cake gaan, gebeuren er gekke dingen: ze botsen, ze breken, en ze veranderen van vorm. De oude simpele formules houden hier geen rekening mee.

2. De Oplossing: Een Digitale Zandbak

In plaats van alleen met formules te rekenen, hebben de onderzoekers een digitale zandbak (een computersimulatie) gebouwd.

• De Simulatie: Ze hebben een stukje aarde in de computer nagebootst. Ze vulden dit met een homogene "deeg" (een simpele, egaal materiaal) en lieten er 30 verschillende "stootjes" op vallen, alsof er 30 mensen tegelijk op de grond springen.
• De Berekening: De computer rekent precies uit hoe de grond trilt, hoe de dichtheid van het gesteente verandert, en welke zwaartekrachtkracht dit uitoefent op de zware gewichten (de testmassa's) binnenin de telescoop.

Het is alsof je in plaats van te gokken hoe een golfje water zich gedraagt in een bad met speelgoed, je eerst een perfecte video maakt van het water en dan pas kijkt wat er gebeurt.

3. De Verrassende Bevinding: Minder P-golven dan gedacht

Bij trillingen in de grond zijn er twee hoofdsoorten golven:

1. P-golven (Longitudinaal): Dit zijn golven die de grond in en uit duwen, alsof je een accordeon in en uit trekt. Deze veranderen de dichtheid van het gesteente en veroorzaken de meeste ruis.
2. S-golven (Transversaal): Dit zijn golven die de grond van links naar rechts schudden, alsof je een laken schudt.

De oude theorieën gingen er vaak van uit dat er ongeveer 1/3 P-golven en 2/3 S-golven waren. Maar in hun nieuwe, gedetailleerde simulatie ontdekten ze iets verrassends:

De verhouding was veel lager! Ze vonden dat er maar ongeveer 14% P-golven waren.

De Analogie:
Stel je voor dat je een orkest hebt waarvan je dacht dat er 3 trompetten (P-golven) en 6 violen (S-golven) waren. Maar toen je naar de digitale opname keek, bleek dat er maar 1 trompet was en 6 violen. Omdat de trompetten (P-golven) het meest storende geluid maken voor de telescoop, betekent dit dat de "muziek" (de ruis) eigenlijk stiller is dan we dachten.

4. Waarom is dit belangrijk?

Dit is een groot nieuws voor de toekomst van de Einstein Telescoop:

• Beter filteren: Als je weet dat er minder "storende trompetten" zijn, kun je je filters beter instellen. Het is makkelijker om de ruis te verwijderen als je precies weet wat je tegenkomt.
• Meer kans op succes: Omdat de ruis lager lijkt te zijn dan eerder gedacht, heeft de telescoop een grotere kans om de fluisterende signalen uit het heelal echt te horen.
• De volgende stap: De huidige simulatie was nog vrij simpel (zoals een vlakke cake). De volgende stap is om de computer te laten rekenen met de echte, complexe geologie van de plek waar de telescoop komt te staan (met bergen, dalen en verschillende steensoorten).

Samenvatting

De onderzoekers hebben een nieuwe, slimme manier bedacht om de "grondtrillingen" te simuleren die de Einstein Telescoop kunnen storen. Ze hebben ontdekt dat de aardse trillingen waarschijnlijk minder storend zijn dan we dachten. Dit betekent dat de kans groter is dat we in de toekomst nog diepere geheimen van het heelal kunnen onthullen.

Het is alsof ze een betere kaart hebben getekend van een mistig landschap, en ze ontdekten dat de mist niet zo dik is als we dachten. Daardoor kunnen we sneller en scherper kijken naar de sterren.

Technische samenvatting

Titel: Een numeriek kader voor de schatting van Newtoniaanse ruis bij het Einstein Telescope: 2D-simulaties voorbij de benadering van vlakke golven

Auteurs: Patrick Schillings, Shi Yao, Johannes Erdmann, Andreas Rietbrock
Publicatie: arXiv:2603.15424v1 (16 maart 2026)

---

1. Het Probleem

Het Einstein Telescope (ET) is een geplande ondergrondse gravitatiegolfobservatorium van de derde generatie, ontworpen om de detectiegevoeligheid uit te breiden tot enkele Hertz. Een kritieke beperking voor de lage-frequentiegevoeligheid (vooral in het bereik van 1,7 tot 6 Hz) is Newtoniaanse ruis.

• Oorzaak: Newtoniaanse ruis wordt veroorzaakt door zwaartekrachtskrachten die werken op de testmassa's van de interferometer, als gevolg van dichtheidsfluctuaties ($\delta\rho$) in de rotsen rondom de detector. Deze fluctuaties worden veroorzaakt door seismische golfbewegingen (P-golven en S-golven).
• Huidige beperkingen: Bestaande schattingen vertrouwen vaak op analytische of semi-analytische modellen die aannemen dat het medium homogeen is of uit horizontale lagen bestaat. Ze verwaarlozen vaak:
  • Geologische heterogeniteit (verspreiding, mode-conversie).
  • Complexe golfinteracties.
  • Correlaties tussen verschillende golfsoorten.
  • De eindige geometrieën van de omgeving.
• Noodzaak: Er is een behoefte aan een numeriek kader dat volledig in staat is om deze complexe interacties op te lossen en nauwkeurigere, locatie-specifieke schattingen te leveren voor de mitigatie van Newtoniaanse ruis.

2. Methodologie

De auteurs presenteren een nieuw numeriek kader dat seismische golfvelden simuleert en deze koppelt aan de berekening van de resulterende zwaartekrachtskrachten.

• Numerieke Oplosser: Het seismische golfveld wordt opgelost met Salvus, een sterk geparallelleerde, GPU-versnelde spectrale-elementen-oplosser (spectral-element solver). Deze methode combineert de geometrische aanpasbaarheid van eindige-elementen met hoge-orde polynoombasisfuncties voor hoge nauwkeurigheid in heterogene media.
• Fysica van Newtoniaanse Ruis:
  • De dichtheidsfluctuaties worden afgeleid uit de continuïteitsvergelijking: $\delta\rho = -\nabla \cdot (\rho \vec{\xi})$.
  • Dit wordt opgesplitst in een bulk-term (dichtheidsfluctuaties in het volume) en een surface-term (dichtheidsverschillen aan grensvlakken, zoals de wanden van de cavernes).
  • De zwaartekrachtskracht op de testmassa wordt berekend door integratie over het volume (bulk) en de oppervlakken (cavernewanden), samengevat in de "dipoolvergelijking".
• Simulatieopzet (Proof of Concept):
  • Dimensie: 2D-simulaties (voor eenvoud, maar het kader is uitbreidbaar naar 3D).
  • Medium: Homogeen elastisch medium (20 km x 8 km) met parameters die de bovenste korst representeren ($c_P = 3000$ m/s, $c_S \approx 1765$ m/s, $\rho \approx 2224$ kg/m³).
  • Scenario 1 (Validatie): Een enkele puntbron op het oppervlak.
  • Scenario 2 (Realistische benadering): Een array van 30 stochastische, willekeurig verdeelde bronnen op het oppervlak om stationaire omgevingsseismische ruis te simuleren.

3. Belangrijkste Bijdragen

1. Validatie van het Kader: De eerste succesvolle koppeling van een spectrale-elementen-seismische simulatie met een volledige Newtoniaanse-ruisberekening, inclusief zowel bulk- als cavernewand-bijdragen.
2. Validatie tegen Analytische Modellen: Het numerieke model reproduceert nauwkeurig de analytische resultaten voor vlakke golven in een homogeen medium, wat de betrouwbaarheid van de methode bevestigt.
3. Schatting van het P-golf-aandeel: Voor het eerst wordt het relatieve aandeel van P-golven versus S-golven in een door bronnen gegenereerd golfveld kwantitatief bepaald via divergentie en rotatie (curl) van het verplaatsingsveld, in plaats van te vertrouwen op aannames.
4. Uitbreidbaarheid: Het kader is ontworpen om eenvoudig te worden uitgebreid naar 3D-simulaties en kan lokale geologische modellen en topografie integreren.

4. Resultaten

• Validatie (Enkele Bron):

  • De simulatie toont een uitstekende overeenkomst met analytische voorspellingen voor zowel de bulk-bijdrage (P-golven) als de cavernewand-bijdrage (P- en S-golven).
  • De "dipoolvergelijking" (die zowel bulk als oppervlakte omvat) voorspelt de totale Newtoniaanse ruis correct, waarbij de S-golven uitsluitend via de cavernewand-bijdrage bijdragen aan de ruis.

• Array van 30 Bronnen (Omgevingsruis):

  • Ruisniveau: De geschatte amplitude spectraaldichtheid (ASD) van de Newtoniaanse ruis ligt binnen de grenzen van eerdere analytische schattingen, maar vertoont fluctuaties veroorzaakt door correlaties tussen P- en S-golven die uit dezelfde bronnen stammen.
  • P-golf Aandeel ($p$): Dit is een cruciale bevinding. In de simulatie van het homogene medium met 30 bronnen bleek het P-golf-aandeel $p \approx 0,137 \pm 0,02$.
  • Vergelijking: Dit is significant lager dan de vaak aangenomen waarde van $p = 1/3$ in eerdere analytische studies.

5. Betekenis en Conclusie

• Verbeterde Mitigatiekansen: Een lager P-golf-aandeel is gunstig voor de mitigatie van Newtoniaanse ruis. Omdat P-golven (longitudinaal) dichtheidsfluctuaties in het bulk-materiaal veroorzaken en S-golven (transversaal) voornamelijk alleen aan grensvlakken, is het scheiden en compenseren van deze signalen complexer bij een hoger P-golf-aandeel. Een lager aandeel suggereert dat de ruis lager kan zijn dan eerder gedacht en dat mitigatiestrategieën effectiever kunnen zijn.
• Toekomstperspectief: Hoewel deze resultaten gebaseerd zijn op een homogeen medium, biedt het kader de basis voor veel realistischere studies. De volgende stappen omvatten het integreren van 3D-heterogene geologieën, reflecties en mode-conversies om de specifieke omstandigheden van de ET-locaties nauwkeurig te modelleren.
• Toepasbaarheid: Het framework kan worden gebruikt om de ontwerpvereisten voor seismische arrays voor ruispredictie te optimaliseren en om de uiteindelijke gevoeligheid van het Einstein Telescope nauwkeuriger te voorspellen.

Conclusie: Dit werk zet een belangrijke stap weg van vereenvoudigde analytische modellen naar een robuust, numeriek kader dat complexe seismische interacties kan oplossen, wat essentieel is voor het realiseren van de lage-frequentiedoelen van het Einstein Telescope.