🌊 Fale grawitacyjne — spirala podwójnego układu i detekcja LIGO
Symuluj spiralne zbliżanie się podwójnej czarnej dziury: obserwuj świergot (chirp) rosnący w częstotliwości i amplitudzie oraz to, jak interferometr LIGO wykrywa przechodzącą falę.
Podobne symulacje
Fizyka i równania
Częstotliwość fali grawitacyjnej podczas spiralowania: f(τ) ∝ (GM_chirp/c³)^(−5/8) · τ^(−3/8), gdzie τ to czas do zlania się obiektów. Masa chirpowa M_c = (m₁m₂)^(3/5)/(m₁+m₂)^(1/5). Odkształcenie h(t) ∝ (πf)^(2/3)·cos(2πft). Częstotliwość zlania się na orbicie ISCO: f_ISCO = c³/(6^(3/2)·π·GM_total).
Detekcja LIGO
LIGO mierzy odkształcenie h = ΔL/L — ułamkową zmianę długości ramienia interferometru wywołaną przez przechodzące fale grawitacyjne. Dla GW150914 h_peak ≈ 10⁻²¹, co oznacza, że ramiona o długości 4 km zmieniły się o ~10⁻¹⁸ m — mniej niż rozmiar protonu. LIGO musi odróżnić ten sygnał od szumu sejsmicznego, fluktuacji termicznych i kwantowego szumu śrutowego.
Sygnał ćwierkania (chirp)
Gdy dwie czarne dziury spiralnie zbliżają się do siebie, tracą energię na promieniowanie grawitacyjne, więc krążą coraz szybciej i coraz bliżej. Powstaje charakterystyczny „chirp": rosnąca częstotliwość i amplituda kulminujące w momencie zlania się. Masę chirpową można wyodrębnić z ewolucji częstotliwości z niezwykle wysoką precyzją — to najczulszy pomiar wykonywany przez LIGO.
O tej symulacji
Ten symulator renderuje ostatnią fazę spiralnego zbliżania się podwójnej czarnej dziury jako shader fragmentów WebGL działający w czasie rzeczywistym. Dwie masy punktowe m₁ i m₂ krążą wokół wspólnego środka w odległości, która maleje w czasie (od 0,85 do 0,04 jednostek shadera w miarę jak postęp zmienia się od 0 do 1), więc orbitalna częstość kątowa ω rośnie w przybliżeniu jak separacja−3/2 — stylizowany substytut rzeczywistej utraty energii na promieniowanie grawitacyjne. Częstotliwość emitowanej fali wynosi 2ω, a jej kwadrupolowy wzorzec promieniowania (zależność kątowa cos 2θ) w widoczny sposób odkształca tło w postaci siatki czasoprzestrzennej, kolorując rozciąganie na niebiesko, a ściskanie na czerwono. Masa chirpowa Mc = (m₁m₂)3/5/(m₁+m₂)1/5 oraz częstotliwość zlania się fISCO = c³/(63/2πGMtotal) są obliczane na podstawie rzeczywistych stałych fizycznych G, c i M☉ i wyświetlane na żywo obok animacji.
Co przedstawia
Dwa świecące punkty (czarne dziury) spiralnie zbliżają się do siebie na odkształconej siatce reprezentującej czasoprzestrzeń. W miarę zbliżania się, koncentryczne fale odkształcenia grawitacyjnego rozchodzą się na zewnątrz, rosnąc częstotliwością i amplitudą — „chirp" — aż jasny błysk zlania się oznacza najbardziej wewnętrzną stabilną orbitę kołową (ISCO) i krótkie wygaszanie (ringdown).
Jak z tego korzystać
Przeciągnij suwaki m₁ i m₂ (w masach Słońca), aby zmienić układ podwójny, lub kliknij ustawienie wstępne — GW150914, GW170817, Lekkie BH lub Masywne — aby przejść od razu do znanej konfiguracji. Użyj Pauzy/Restartu oraz suwaka prędkości, aby kontrolować odtwarzanie, i obserwuj panel statystyk pokazujący na żywo masę chirpową, częstotliwość orbitalną/fali grawitacyjnej, separację i fazę spiralowania.
Czy wiesz, że...
Pierwsza detekcja LIGO, GW150914 (2015), pochodziła od dwóch czarnych dziur o masach około 36 i 29 mas Słońca, które zlały się 1,3 miliarda lat świetlnych od nas, wytwarzając szczytowe odkształcenie rzędu 10⁻²¹ — zniekształcenie tysiące razy mniejsze niż proton na 4-kilometrowym ramieniu detektora. GW170817, zlanie się gwiazd neutronowych, zaobserwowano także w świetle, dając początek astronomii wielomesengerowej.
Najczęściej zadawane pytania
Czym jest masa chirpowa i dlaczego to kluczowa mierzalna wielkość?
Masa chirpowa Mc = (m₁m₂)3/5/(m₁+m₂)1/5 to jedyna kombinacja mas obu składników, która determinuje, jak szybko częstotliwość fali grawitacyjnej rośnie podczas spiralowania. Ponieważ ewolucję częstotliwości można mierzyć z niezwykłą precyzją czasową, LIGO potrafi wyznaczyć Mc ze znacznie większą dokładnością niż m₁ i m₂ z osobna — to pierwsza i najbardziej wiarygodna liczba wyodrębniana z każdej detekcji, a ten symulator przelicza ją na żywo na podstawie wartości suwaków.
Dlaczego zarówno częstotliwość, jak i amplituda fali rosną w miarę zbliżania się zlania (tzw. „chirp")?
Gdy obie masy tracą energię na promieniowanie grawitacyjne, ich orbita nieuchronnie się kurczy, a dynamika keplerowska wymusza wzrost częstotliwości orbitalnej (a więc i częstotliwości fali grawitacyjnej) wraz ze zmniejszaniem się separacji. Ponieważ ciaśniejsze, szybsze orbity promieniują też silniej, amplituda fali rośnie równocześnie z jej częstotliwością — połączony wzrost częstotliwości i amplitudy nadaje sygnałowi charakterystyczny rosnący dźwięk i wygląd „chirpu", dokładnie zgodnie z zakodowaną w tej symulacji zależnością omega ∝ separacja−3/2.
Czym jest ISCO i dlaczego zlanie się następuje właśnie tam?
Najbardziej wewnętrzna stabilna orbita kołowa (ISCO) to najmniejsza orbita wokół zwartego obiektu, na której nadal może istnieć stabilna droga kołowa; bliżej niż to ogólna teoria względności nie pozostawia żadnej stabilnej orbity, a oba ciała wpadają na siebie. Ta symulacja uruchamia błysk zlania się, gdy postęp spiralowania zbliża się do 1 (separacja ISCO), i oblicza odpowiadającą jej rzeczywistą częstotliwość jako fISCO = c³/(63/2πGMtotal), która dla czarnych dziur o masach gwiazdowych mieści się w najczulszym paśmie LIGO, rzędu kilkuset Hz.
Jak shader przekształca fizykę w niebiesko-czerwoną falującą siatkę?
Odległość i kąt każdego piksela od środka układu podwójnego zasilają opóźnioną w czasie fazę (uwzględniającą propagację fali na zewnątrz), pomnożoną przez kwadrupolowy czynnik kątowy cos 2θ, będący poprawnym wzorcem kątowym m=2 dla promieniowania grawitacyjnego pary orbitujących mas. Ta połączona wartość to lokalne odkształcenie, wykorzystywane do promieniowego przesunięcia („odkształcenia") tła siatki i jej kolorowania — na niebiesko tam, gdzie matematyka shadera wskazuje rozciąganie, na czerwono tam, gdzie wskazuje ściskanie — co daje bezpośredni wizualny odpowiednik wzorca polaryzacji „+", do wykrywania którego zbudowano prawdziwe interferometry.
Czym różnią się ustawienia wstępne GW150914, GW170817, Lekkie BH i Masywne?
GW150914 (36+29 mas Słońca) odtwarza pierwszą detekcję LIGO z 2015 roku; GW170817 (około 1,4+1,4 masy Słońca) reprezentuje zlanie się dwóch gwiazd neutronowych — lżejszych obiektów, które spiralują znacznie dłużej i zostały również zaobserwowane w całym spektrum elektromagnetycznym jako kilonowa. Lekkie BH (5+5) i Masywne (50+50) pozwalają porównać, jak całkowita masa wpływa na masę chirpową, czas do zlania się (cięższe układy podwójne zlewają się szybciej w tym modelu) oraz częstotliwość ISCO, która skaluje się odwrotnie proporcjonalnie do masy całkowitej.