Nauka #29 – Ogólna teoria względności i zakrzywiona czasoprzestrzeń: od zasady równoważności do czarnych dziur

Ogólna teoria względności (OTW) to teoria grawitacji Einsteina z 1915 roku. Zastępuje newtonowską siłę grawitacji ideą, że masywne obiekty zakrzywiają geometrię czasoprzestrzeni, a swobodnie spadające ciała (w tym światło) poruszają się po najprostszych możliwych trajektoriach (geodezyjnych) w tej zakrzywionej geometrii. OTW przeszła każdy eksperymentalny test z niezwykłą precyzją na przestrzeni ponad stu lat. Ten wpis prowadzi od zasady równoważności aż po termodynamikę czarnych dziur.

1. Zasada równoważności i geodezyjne

Einsteinowska Zasada Równoważności (EEP) głosi, że skutki grawitacji są lokalnie nieodróżnialne od skutków przyspieszenia. Ściślej: w wystarczająco małym obszarze czasoprzestrzeni obowiązują prawa szczególnej teorii względności w swobodnie spadającym (inercjalnym) układzie odniesienia. Implikuje to, że grawitacja musi zakrzywiać światło (jeśli poziomy promień światła zakrzywia się w przyspieszającej rakiecie, musi się też zakrzywiać w polu grawitacyjnym) i powodować grawitacyjną dylatację czasu.

Interwał czasoprzestrzenny (metryka ogólna):
  ds² = g_μν dxμ dxν
  Sygnatura: (−,+,+,+)  [ds²<0 podobne do czasu, >0 podobne do przestrzeni, =0 zerowe]

Czas własny:  dτ² = −ds²/c²

Równanie geodezyjne (swobodny spadek w zakrzywionej czasoprzestrzeni):
  d²xµ/dτ² + Γµ_(νρ) (dxν/dτ)(dxρ/dτ) = 0

Symbole Christoffela (koneksja Levi-Civity):
  Γµ_(νρ) = (1/2) gµσ (∂_ν g_σρ + ∂_ρ g_σν − ∂_σ g_νρ)
  Symetryczne w dolnych indeksach: Γµ_(νρ) = Γµ_(ρν)
  Nie jest tensorem (znika w lokalnie inercjalnym układzie)

Pochodna kowariantna:
  ∇_ν Vµ = ∂_ν Vµ + Γµ_(νρ) Vρ   (wektor kontrawariantny)
  ∇_ν ω_µ = ∂_ν ω_µ − Γρ_(νµ) ω_ρ   (wektor kowariantny, 1-forma)

Transport równoległy:
  DVµ/dλ = (dVµ/dλ) + Γµ_(νρ)(dxν/dλ)Vρ = 0
  Geodezyjna = krzywa samo-równoległa: DUµ/dτ = 0  gdzie Uµ = dxµ/dτ
          

2. Krzywizna: Riemann, Ricci i tensor Einsteina

Tensor krzywizny Riemanna mierzy niepowodzenie transportu równoległego wzdłuż zamkniętej pętli — równoważnie, komutator pochodnych kowariantnych. Jego zwężenia, tensor Ricciego i skalar Ricciego, pojawiają się bezpośrednio w równaniach pola Einsteina.

Tensor Riemanna (16&sup4; = 256 składowych; 20 niezależnych dzięki symetrii):
  Rσ_(μνρ) = ∂_νΓσ_(ρµ) − ∂_ρΓσ_(νµ) + Γσ_(νλ)Γλ_(ρµ) − Γσ_(ρλ)Γλ_(νµ)

Mierzy krzywiznę przez komutator:
  [∇_μ, ∇_ν]Vρ = Rρ_(σμν) Vσ

Symetrie:
  R_μνρσ = −R_νμρσ = −R_μνσρ = R_ρσμν
  Tożsamość Bianchiego (algebraiczna): R_[μνρσ] = 0
  Tożsamość Bianchiego (różniczkowa): ∇_[λR|μν|ρσ] = 0

Tensor Ricciego (symetryczny, 10 niezależnych):
  R_μν = Rρ_(µρν)  (zwężenie po 1. i 3. indeksie)

Skalar Ricciego:
  R = gμνR_µν

Tensor Einsteina (symetryczny, bezdywergencyjny):
  G_μν = R_μν − (1/2)g_μνR
  ∇ν Gμν = 0  (skontraktowana tożsamość Bianchiego → zachowanie Tμν)

Tensor Weyla C_μνρσ:
  Bezśladowa część tensora Riemanna
  Rządzi siłami pływowymi w próżni (R_μν=0 nie oznacza Rσ_µνρ=0)
          

3. Równania pola Einsteina

Równania pola Einsteina (EFE) opisują, jak rozkład energii i pędu jest źródłem krzywizny czasoprzestrzeni. To 10 sprzężonych, nieliniowych równań różniczkowych cząstkowych dla tensora metryki g_μν.

Równania pola Einsteina:
  G_μν + Λ g_μν = (8πG/c&sup4;) T_μν

  G_μν = krzywizna czasoprzestrzeni (geometria)
  T_μν = tensor energii-pędu (zawartość materii/energii)
  Λ ≈ 1.1 × 10−&sup5;² m−² = stała kosmologiczna (ciemna energia)

Tensor energii-pędu doskonałej cieczy:
  Tµν = (ρ + p/c²) UµUν + p gµν
  Pył: p=0 → Tµν = ρc² UµUν
  ∇_ν Tµν = 0  ↔  równania Eulera + zachowanie energii

Postać z odwróconym śladem (przydatna do linearyzacji):
  Rµν = (8πG/c&sup4;)(Tµν − (1/2)gµνT)
  Próżnia: Tµν = 0 → Rµν = 0 (Ricci-płaska ≠ płaska)

Grawitacja linearyzowana hµν = gµν − ηµν  (|h| ≪ 1):
  Perturbacja z odwróconym śladem: h̄µν = hµν − (1/2)ηµνh
  Warunek cechowania Lorenza: ∂νh̄µν = 0

  □h̄µν = −(16πG/c&sup4;) Tµν   (równanie falowe)
  □ = ηµν∂µ∂ν = −∂&sub2;t/c² + ∇²   (d'Alembertian)

Granica newtonowska: g_00 ≈ −(1 + 2Φ/c²), Φ = potencjał grawitacyjny
  ∇²Φ = 4πGρ  (odtworzone równanie Poissona)
          

4. Rozwiązanie Schwarzschilda

Karl Schwarzschild znalazł pierwsze dokładne rozwiązanie EFE w grudniu 1915 roku — zaledwie kilka tygodni po opublikowaniu przez Einsteina ostatecznych równań — będąc na froncie wschodnim I wojny światowej. Metryka Schwarzschilda opisuje geometrię czasoprzestrzeni na zewnątrz sferycznie symetrycznej, nierotującej, nienaładowanej masy.

Metryka Schwarzschilda (symetria sferyczna, próżnia, M, J=0, Q=0):
  ds² = −(1−r_s/r)c²dt² + (1−r_s/r)−¹dr² + r²dΩ²
  r_s = 2GM/c² = promień Schwarzschilda

  Przykłady:
    Słońce:  r_s = 2,95 km  (R_słońce = 696 000 km → daleko poza horyzontem)
    Ziemia:  r_s = 8,9 mm   (R_ziemia = 6371 km)
    M87*:    r_s = 2 × 10^10 km  (M ≈ 6,5 × 10^9 M_słońce)

Potencjał efektywny dla cząstek masywnych:
  V_eff(r) = (c²/2)[(1−r_s/r)(1 + L²/r²c²) − 1]
  gdzie L = specyficzny moment pędu
  Najbardziej wewnętrzna stabilna orbita kołowa (ISCO): r_ISCO = 3r_s = 6GM/c²

Precesja orbity Merkurego:
  Δφ_prec = 6πGM / [ac²(1−e²)]  na obieg
  Merkury: a = 0,387 j.a., e = 0,206 → Δφ = 43,0 sekund łuku/wiek ✓

Grawitacyjne przesunięcie ku czerwieni:
  ν_obs/ν_emit = √(g_00(r_emit) / g_00(r_obs))
  z = (1 − r_s/r)^(−1/2) − 1  ≈ GM/(rc²) dla r ≫ r_s
  Pound-Rebka (1959): zmierzono z = 2,46 × 10^−15 przy h=22,5m  ✓

Sfera fotonowa: r_ps = (3/2)r_s = 3GM/c²
  Niestabilna orbita kołowa dla fotonów
  Promień cienia na obrazach EHT: r_shadow = 3√3 GM/c²
  Obraz M87* EHT (2019): średnica cienia 40±3 µas  (przewidywana 39±1 µas) ✓
          

5. Fale grawitacyjne

Przyspieszająca masa-energia emituje fale grawitacyjne — zmarszczki geometrii czasoprzestrzeni rozchodzące się z prędkością światła. Cechowanie poprzeczno-bezśladowe (TT) wybiera fizyczne stopnie swobody: dwie polaryzacje, h₊ i h_×, które naprzemiennie ściskają i rozciągają kierunki prostopadłe podczas przechodzenia fali.

Rozwiązanie w cechowaniu TT (dalekie pole, r ≫ λ_GW):
  h_ij^TT(t,r) = (2G/rc&sup4;) Ï_ij^TT(t_ret)
  gdzie t_ret = t − r/c  (czas opóźniony)
  I_ij = ∫ ρ x_i x_j d³x  (tensor momentu kwadrupolowego masy)

Całkowita moc promieniowana (wzór kwadrupolowy):
  P = −G/(5c&sup5;) ⟨Q̈_ij Q̈^ij⟩   gdzie Q_ij = I_ij − (1/3)δ_ij I_kk

Układ podwójny (dwie równe masy m, separacja a, orbita kołowa f_orb):
  P = (32/5) G^4 m^5 / (c^5 a^5)  (wzór Petersa)
  a(t) = a_0 (1 − t/t_merge)^(1/4)
  t_merge = (12/19) c_0^4/(B F(e)) ≈ 12/85 (c²0 a^4_0)/(4G^3 m² M)

LIGO GW150914 (14 września 2015):
  Dwie czarne dziury: m_1 ≈ 36 M_słońce, m_2 ≈ 29 M_słońce
  Szczytowe odkształcenie: h ~ 10^−²¹ przy 150 Hz
  Długość ramienia LIGO: L = 4 km → ΔL = hL/2 ~ 10^−¹&sup8; m ≈ 1/1000 promienia protonu
  Końcowa CD: M_f ≈ 62 M_słońce (3 M_słońce → energia GW: ~ 5 × 10^47 J)

Sieci zegarów pulsarowych (PTA):
  Tło fal nanoherców wykryte w 2023 (NANOGrav, EPTA, PPTA, CPTA)
  f ~1–100 nHz, T_obs ~ 12,5 lat linii bazowej, ~67 milisekundowych pulsarów

LISA planowana (2034, ESA/NASA):
  3 statki kosmiczne, ramię 2,5 Mkm, 0,1 mHz–1 Hz
  Cele: fuzje SMBH z<20, galaktyczne białe karły podwójne, fale EMRI
          

6. Czarne dziury i promieniowanie Hawkinga

Czarna dziura to obszar czasoprzestrzeni, w którym prędkość ucieczki przekracza c: żaden sygnał przyczynowy nie może wydostać się spoza horyzontu zdarzeń. Rozwiązanie Kerra (1963) rozszerza metrykę Schwarzschilda na wirujące czarne dziury i przewiduje ergosferę — obszar, w którym czasoprzestrzeń jest tak gwałtownie "wleczona", że nic nie może pozostać nieruchome. Stephen Hawking pokazał w 1974 roku, że kwantowa teoria pola w pobliżu horyzontu prowadzi do promieniowania termicznego, nadając czarnym dziurom skończoną temperaturę i czas życia.

Metryka Kerra (współrzędne Boyera-Lindquista, M, J=Ma, Q=0):
  ds² = −(1−r_s r/Σ)c²dt² − (2r_s ra sin²θ/Σ)c dt dφ
           + (Σ/Δ)dr² + Σdθ² + (r²+a²+r_s ra²sin²θ/Σ)sin²θ dφ²
  Σ = r² + a²cos²θ,  Δ = r² − r_s r + a²
  a = J/(Mc) = specyficzny moment pędu  (0 ≤ a ≤ GM/c)

Zewnętrzny horyzont zdarzeń (największy pierwiastek Δ=0):
  r_+ = GM/c² + √((GM/c²)² − a²)
  Ekstremalny Kerr: a = GM/c² → r_+ = GM/c² (najmniejszy możliwy promień CD)

Ergosfera (gdzie g_tt = 0):
  r_ergo = GM/c² + √((GM/c²)² − a²cos²θ)
  Między ergosferą a horyzontem można wydobyć moment pędu (proces Penrose'a)

Proces Penrose'a:
  Cząstka rozpada się wewnątrz ergosfery: jedna spada z ujemną E (względem nieskończoności)
  Druga ucieka z E_esc > E_początkowe → CD traci rotację
  Maksymalna wydobywalna energia: ~21% masy spoczynkowej (dla ekstremalnego Kerra)
  Proces Blandforda-Znajka: wersja magnetyczna, zasila dżety AGN

Promieniowanie Hawkinga (1974 — efekt kwantowy):
  T_H = ℏc³ / (8πGM k_B)
  CD o masie Słońca: T_H ≈ 60 nanokelwinów  (pomijalne wobec CMB 2,7 K)
  CD o masie asteroidy (M ~ 10^15 g): T_H ~ 10^11 K → wyparowuje właśnie teraz

Entropia czarnej dziury (Bekensteina-Hawkinga):
  S_BH = k_B A / (4 l_Pl²)   gdzie l_Pl = √(ℏG/c³) = 1,6 × 10^−35 m
  A = 4πr_+²  (pole powierzchni horyzontu zdarzeń; rośnie w procesach klasycznych: 2. zasada)

Paradoks informacji:
  Promieniowanie Hawkinga jest dokładnie termiczne → żadna informacja nie ucieka?
  Narusza unitarność mechaniki kwantowej (debata o krzywej Page'a)
  Rozwiązania: fuzzballe (teoria strun), firewall (AMPS), formuła wyspy (holografia)
          

Wypróbuj te symulacje