Numerická relativita je obor obecné teorie relativity , který vyvíjí a využívá numerické metody a algoritmy pro počítačovou simulaci fyzikálních procesů v silných gravitačních polích , kdy je potřeba řešit Einsteinovy rovnice numericky . Hlavní fyzikální systémy , jejichž popis vyžaduje numerickou relativitu, se vztahují k relativistické astrofyzice a zahrnují gravitační kolaps , neutronové hvězdy , černé díry , gravitační vlny a další objekty a jevy, pro jejichž adekvátní popis je nutné odkázat na úplnou obecnou teorii relativity bez konvenčních aproximací slabá pole a nízké rychlosti (jako v postnewtonských expanzích a poruchové teorii na pozadí exaktních řešení Einsteinových rovnic ) [1] .
Modelování v této oblasti vyžaduje speciální numerické metody kvůli složitosti a nelinearitě Einsteinových rovnic (např. hyperbolicita a správnost formulace Cauchyho problému jejich časového vývoje závisí na reprezentaci rovnic, stejně jako na počátečním a okrajové podmínky [2] ), a také pro většinu trojrozměrných úloh - vysoký výpočetní výkon dostupný pouze moderním superpočítačům . V současné době je v numerické relativitě relevantní výzkum v oblasti modelování relativistických blízkých dvojhvězd a souvisejících gravitačních vln, jakož i mnoha dalších matematických a astrofyzikálních problémů [1] .
Hlavním účelem numerické relativity je studium gravitačních polí , jejichž přesná analytická forma není známa. Gravitační pole, jejichž podobu hledáme výpočty, mohou být buď plně dynamická , nebo stacionární nebo statická a mohou obsahovat i materiálová pole [~ 1] nebo být vakuová. V případě stacionárních a statických řešení lze pro studium stability těchto konfigurací použít numerické metody. V případě dynamických gravitačních polí lze problém rozdělit na dvě části, které vyžadují různé metody řešení: problém počátečních hodnot a problém evoluce [3] .
Numerická teorie relativity se používá při studiu kosmologických modelů , kritických jevů při gravitačním kolapsu , stejně jako procesů zahrnujících černé díry a neutronové hvězdy , zejména jejich slučování a perturbace . V každém z těchto případů je nutné sledovat vývoj časoprostoru, pro který lze Einsteinovy rovnice znázornit několika způsoby. Nejoblíbenější jsou metody Cauchyho problému , ale používá se také metoda charakteristik [4] a metody založené na Reggeho kalkulu [5] . Všechny výše uvedené metody začínají „snímkem“ gravitačního pole na nějakém hyperpovrchu , tj. z počátečních dat, a poté sledují jeho vývoj k dalším blízkým hyperpovrchům, pohybujícím se v čase [6] .
Jako ve všech problémech numerické analýzy je i v numerické relativitě věnována velká pozornost stabilitě a konvergenci numerických řešení, přípustným počátečním a okrajovým podmínkám. Specifikem numerické relativity jsou komplikace způsobené přítomností kalibračních a souřadnicových podmínek , stejně jako různé reprezentace Einsteinových rovnic a jejich vliv na schopnost získat přesná numerická řešení.
Mnohé z numerických technik používaných v klasické teorii pole nejsou použitelné v obecné teorii relativity, čímž se práce v této oblasti liší od výzkumu v oblasti numerické relativity. Ve velkých problémech však numerická relativita sdílí mnoho aspektů s jinými výpočetními vědami, jako je výpočetní dynamika tekutin , elektrodynamika a mechanika tuhého těla . Vědci zabývající se numerickou relativitou často pracují po boku aplikovaných matematiků a přicházejí do kontaktu s oblastmi matematiky, jako je numerická analýza , paralelní výpočty , parciální diferenciální rovnice a geometrie [7] .
Albert Einstein publikoval konečnou verzi obecné teorie relativity v roce 1915 [8] . Tato teorie, stejně jako speciální teorie relativity , která jí předcházela , popisuje prostor a čas jako jediný objekt – časoprostor , jehož vývoj se řídí Einsteinovými rovnicemi . Tvoří systém sdružených nelineárních parciálních diferenciálních rovnic . Ve století, které uplynulo od odvození těchto rovnic, je známo pouze relativně malé množství fyzikálně relevantních exaktních analytických řešení těchto rovnic a většina z nich je odvozena za předpokladu vysoké symetrie, což zjednodušuje řešení rovnic. , jako jsou Friedmannova řešení pro homogenní a izotropní vesmír [9] .
Obor numerické relativity vznikl z touhy studovat obecnější a fyzikálně použitelná řešení Einsteinových rovnic tím, že je řešíme numericky přibližně. Nezbytnou podmínkou takového řešení bylo provést zpětné rozdělení jednoho čtyřrozměrného časoprostoru na rozdělený trojrozměrný prostor a jednorozměrný čas, tzv. rozdělení 3 + 1 . Navíc to může být provedeno mnoha různými způsoby, což může výrazně zkomplikovat nebo zjednodušit problém integrace výsledných rovnic. První docela úspěšný pokus o rozdělení provedli Richard Arnowitt, Stanley Deser a Charles Misner koncem 50. let 20. století v hamiltonovském formalismu podél cesty naznačené Diracem . Vyvrcholilo to získáním rovnic, které tvoří tzv. ADM formalismus, Arnowitt-Deser-Miznerův formalismus [10] . Ačkoli se z technických důvodů ukázalo, že tyto konkrétní rovnice nejsou pro numerickou integraci příliš vhodné - jsou pouze slabě hyperbolické , a proto se ve skutečných výpočtech používají zřídka - naprostá většina praktických přístupů k numerické relativitě používá dělení 3 + 1 blízké který se používá ve formalismu ADM . Takové rozdělení vede k přeformulování Einsteinových rovnic ve formě Cauchyho problému s omezeními na počáteční hodnoty, který je již přístupný numerickému řešení na počítačích [11] .
Souřadnice v časoprostoru nelze jednoznačně určit, proto i při fixaci souřadnic na počáteční hyperpovrchu, při pohybu na sousední hyperpovrchu, mohou být časové a prostorové souřadnice „tlačeny“ různými způsoby v různých bodech (již ve speciální teorii relativity, směr a rychlost toku času se v různých inerciálních vztažných soustavách neshodují), což je specifikum numerické relativity. Tato měřická volnost - která neovlivňuje fyzikální procesy, ale pouze mění jejich popis v souřadnicích a v souladu s tím i řešené rovnice - se projevuje libovolnou volbou funkce pohybu a posuvu , "tlačením" bodů s pevnými body. prostorové souřadnice od počátečního k sousednímu hyperpovrchu vpřed v čase — , respektive do stran v prostoru — . Možnost zvolit si tyto funkce je potenciální výhodou pro numerické řešení rovnic, ale bylo zjištěno, že mnoho "přirozených" voleb těchto souřadnicových nebo kalibračních podmínek způsobuje numerické nestability v řešeních, což vede k přerušení simulace [12] .
V době zveřejnění původních prací o formalismu ADM vývoj výpočetní techniky neumožňoval výpočty pomocí jejich rovnic pro jakýkoli problém jakékoli rozumné velikosti. Historicky první pokus o numerické řešení Einsteinových rovnic provedli Hahn a Lindqvist v roce 1964 [13] a poté v 70. letech Smarr [14] [15] a Eppley [16] . Tyto rané pokusy souvisely s vývojem Misnerových počátečních dat v osově symetrických prostorech (také známých jako „dimenze 2+1“). Přibližně ve stejné době napsal Zvi Piran první kód, který sledoval vývoj válcově symetrického systému vyzařujícího gravitační záření [17] . V tomto vývoji Piran inicioval mnoho konceptů, které se nyní používají v numerické relativitě, jako je volná evoluce a omezená evoluce, metody, které přistupují k problému evoluce počátečních datových omezení v čase různými způsoby [18] [19] . Použití symetrie snížilo požadované požadavky na paměť a výpočetní výkon, což vědcům umožnilo k řešení tohoto problému použít tehdejší superpočítače [17] .
První realistické výpočty skutečného astrofyzikálního problému, spinningového kolapsu, provedli počátkem 80. let 20. století Richard Stark a Zvi Piran [20] , ve kterých byly poprvé vypočteny gravitační vlny vyzařované formující se rotující černou dírou. Za téměř dvě desetiletí od této publikace bylo zveřejněno pouze několik nových výsledků v numerické relativitě, pravděpodobně kvůli nedostatku počítačů dostatečně výkonných k vyřešení těchto problémů. V 90. letech 20. století Binary Black Hole Grand Challenge Alliance úspěšně simulovala čelní srážku dvou černých děr pomocí zjednodušení vyplývajících z osové symetrie problému . Ve fázi následného zpracování byla skupina schopna vypočítat horizont událostí pro výsledné řešení [21] .
Některé z prvních známých pokusů o numerické řešení Einsteinových rovnic v plné 3D prostorové geometrii se zaměřily na nerotující Schwarzschildovu černou díru , která je statickým a sféricky symetrickým řešením Einsteinových rovnic. Je to vynikající test pro metody numerické relativity, protože za prvé je řešení známé v přesné analytické formě, se kterou lze porovnávat numerické výsledky, za druhé je statické a každá nerotující černá díra by se k němu měla časem sblížit. a za třetí obsahuje jeden z nejobtížnějších objektů pro numerické modelování - fyzickou gravitační singularitu uprostřed. Jeden z prvních pokusů získat toto řešení numericky provedl Anninos et al., v roce 1995 [22] . V této práci poznamenali:
Pokrok v 3D numerické relativitě je částečně brzděn nedostatkem počítačů s dostatkem paměti a výpočetního výkonu pro provádění výpočtů ve 3D s dobrým rozlišením.
Původní text (anglicky)[ zobrazitskrýt] Pokroku v trojrozměrné numerické relativitě částečně brání nedostatek počítačů s dostatečnou pamětí a výpočetním výkonem k provádění dobře vyřešených výpočtů 3D prostoročasů.V průběhu let, kromě toho, že počítače byly stále výkonnější, byly různými výzkumnými skupinami vyvinuty alternativní techniky ke zvýšení efektivity práce na počítači. Za prvé, skupina Lazarus vyvinula metody, které používaly rané výsledky krátkých simulací řešících nelineární rovnice ADM pro slučování černých děr, aby poskytla počáteční data pro robustnější kód založený na linearizovaných rovnicích teorie poruch jediné černé díry [23]. . Poté, s ohledem na modelování černých děr, byly vyvinuty dvě techniky, jak se vyhnout problémům spojeným s existencí fyzikální singularity při řešení rovnic: (1) eliminace a (2) „prick“ metoda [24] . Kombinace těchto metod s nalezenými vhodnými souřadnicovými podmínkami umožnila v roce 2005 učinit průlom v modelování binárních černých děr, který začal prací Pretoriuse [25] . O několik let později umožnila numerická stabilita nových metod simulovat téměř libovolné konfigurace binárních černých děr, které před sloučením popisují desítky a stovky otáček kolem sebe. V numerické relativitě se navíc začaly používat metody adaptivního zpřesňování výpočetní sítě, které se dříve používaly ve výpočetní dynamice tekutin [26] .
Projekt LazarusProjekt Lazarus (1998–2005) byl vyvinut po Velké výzvě jako technika pro extrakci astrofyzikálně relevantních výsledků z tehdy dostupných krátkých numerických simulací sloučení binárních černých děr. V té době všechny známé pokusy o integraci Einsteinových rovnic pro časoprostor binárních černých děr na superpočítačích kvůli různým druhům nestabilit nebyly schopny pokročit ani do dokončení jedné úplné rotace systému. V rámci projektu výzkumníci kombinovali přibližné metody před ( postnewtonské trajektorie ) a po transformaci dvojice děr v jednu (poruchy jednotlivých černých děr) s kompletním numerickým řešením samotného procesu [23] .
Přístup projektu Lazarus byl v té době nejlepším přístupem k problému binárních černých děr a poskytl velké množství výsledků dostatečně přesných pro astrofyzikální aplikace, jako jsou hodnoty energie a momentu hybnosti unášené gravitačními vlnami [27 ] [28] , dále hybnost při slučování černých děr různých hmotností [29] a hodnoty konečné hmotnosti, hybnosti a momentu hybnosti vznikající černé díry [30] . Metody projektu také umožnily vypočítat podrobné formy gravitačních vln emitovaných během sloučení – což bylo důležité pro gravitační dalekohledy , a předpovídaly, že srážky černých děr by měly být doprovázeny nejsilnějšími výbuchy energie ve vesmíru, když více energie se uvolní ve formě gravitačního záření ve zlomku sekundy, než vyzáří všechny hvězdy galaxie během své existence - gravitační záření odnese několik procent původní redukované hmoty systému [31] .
Metoda vyloučeníV technice excize , která byla poprvé navržena na konci 90. let 20. století [32] , je část časoprostoru uvnitř horizontu událostí obklopujícího singularitu černé díry jednoduše vyloučena z evoluce . Teoreticky by to nemělo mít vliv na rozhodování mimo horizont událostí kvůli principu kauzality a vlastnostem horizontu - protože žádné fyzikální interakce pod horizontem nemohou mít žádný vliv na fyziku mimo něj. Pokud tedy jednoduše nevyřešíte rovnice uvnitř černé díry, stále můžete získat přesné skutečné řešení mimo ni. Vnitřní dynamiku je možné „eliminovat“ tím, že na hranici uvnitř horizontu, obejmeme singularitu, vložíme okrajové podmínky nepřítomnosti vystupujících vln [33] .
Přestože bylo použití vylučovací techniky velmi úspěšné, má dva malé problémy. První je, že je třeba pečlivě vybírat a používat souřadnicové podmínky. Zatímco fyzikální efekty se nemohou šířit z horizontu ven, souřadnicové efekty ano. Pokud jsou například zavedeny podmínky eliptických souřadnic, změny mřížky uvnitř černé díry se mohou okamžitě šířit směrem ven přes horizont [34] . To znamená, že pro aplikaci eliminační metody je nutné použít souřadnicové podmínky hyperbolického typu, ve kterých jsou charakteristické rychlosti šíření souřadnicových efektů menší nebo rovné rychlosti světla (například při použití harmonických souřadnicových podmínek) [35] . Druhým problémem je, že jelikož se černá díra pohybuje, oblast vyloučení se musí neustále pohybovat v souladu s ní [33] .
Eliminační metoda byla vyvíjena několik let, přičemž byly nalezeny nové kalibrační podmínky zvyšující stabilitu postupu řešení a byla prokázána schopnost vyloučených oblastí pohybovat se po výpočetní síti [36] [37] [38] [ 39] [40] [35] . První stabilní dlouhý výpočet oběžné dráhy a sloučení dvou černých děr pomocí této techniky byl publikován v roce 2005 [25] .
Metoda vstřikováníU metody punkce je řešení rozděleno na analytickou část [41] , která obsahuje singularitu černé díry — punkce, a numericky konstruovanou část, která singularitu neobsahuje. Tato metoda je zobecněním Brill-Lindquistova algoritmu [42] pro počáteční data s černými dírami v klidu a může být dále zobecněna na Bowen-Yorkův algoritmus [43] pro počáteční data s rotujícími a pohyblivými černými dírami. Do roku 2005 všechny publikované příklady použití metody „píchnutí“ vyžadovaly, aby byly souřadnice všech vpichů po celou dobu trvání simulace pevně dané. Černé díry ve vzájemné těsné blízkosti se samozřejmě budou pohybovat vlivem gravitačních sil, takže pevné souřadnice kolíků znamenají, že se souřadnicové systémy "natahují" nebo "deformují", což vede k numerické nestabilitě v některých fázích simulace. Podobné efekty vyvolává i použití jiné metody - vyhýbání se singularitám, kdy černé díry vznikají v simulaci kolapsem hmoty a souřadnicové podmínky jsou voleny tak, aby trojrozměrná hyperplocha vyvíjející se v čase nedosáhla singularita až do konce výpočtů, tvořící kolem ní protáhlý "roh" [44] .
V roce 2005 výzkumníci poprvé prokázali možnost pohybu vpichů podél souřadnicového systému, čímž vyřešili některé z raných problémů metody, což umožnilo přesně sledovat dlouhodobý vývoj černých děr [25] [45 ] [46] . Volbou vhodných souřadnicových podmínek a provedením hrubých analytických aproximací fyzikálních polí v blízkosti singularity (protože z černé díry nemohou uniknout žádné fyzikální efekty, drsnost aproximace není důležitá), lze získat numerická řešení problému dvou černých děr. rotující kolem sebe, a také přesně vypočítat jejich gravitační záření [47] .
Adaptivní zjemnění sítěAdaptivní zjemňování sítě jako numerická metoda byla používána ve fyzice dlouho před vzestupem numerické relativity. V něm byl poprvé použit v 80. letech 20. století v dílech Choptwicka při studiu kritických jevů během kolapsu skalárního pole , kdy konfigurace pole jsou na samé hranici mezi konečným vytvořením černé díry a konečnou expanzí ve vesmíru. [48] [49] . Původní práce byla jednorozměrná, protože používali sférickou symetrii, ale pak byla metoda zobecněna na dvě dimenze [50] . Metody dvourozměrného zpřesňování byly také aplikovány na studium nehomogenních kosmologií [51] [52] a Schwarzschildových černých děr [53] . Metody adaptivního zpřesňování se nyní staly standardním nástrojem numerické relativity a kromě studia šíření gravitačních vln generovaných takovými událostmise používají při studiu slučování černých děr a dalších kompaktních objektů [54] [55] .
K dnešnímu dni jsou napsány desítky a stovky prací ročně o numerické relativitě, které prezentují širokou škálu výsledků v oblastech matematiky obecné teorie relativity, gravitačních vln a astrofyziky, získaných řešením problému černých děr rotujících kolem sebe. Použité metody byly zobecněny pro studium astrofyzikálních binárních systémů, včetně neutronových hvězd, černých děr [56] a souborů černých děr [57] . Tyto práce mimo jiné předpovídají, že když se dvě rotující černé díry spojí, výsledná díra může dosáhnout rychlosti až 4000 a dokonce až 10 000 km/s , což jí umožňuje překonat jakoukoli známou galaxii [58] [59] . Simulace také předpovídají obrovské uvolnění energie během slučování, které může být až 8 % celkové hmoty v klidu, a možnost prudké změny osy rotace černé díry , což může vysvětlit změny směrů výtrysků. pozorované v rádiových galaxiích [ 60] . Důležitou linií výzkumu je také vytvoření katalogu forem gravitačního záření ze splývajících černých děr, bez kterého je hledání těchto signálů v datech z detektorů jako LIGO a VIRGO mnohem méně citlivé [61] .
Přesnost moderních metod numerické relativity bylo možné v praxi ověřit ihned po objevu gravitačních vln . Bylo zjištěno, že signál GW150914 je ve shodě s předpovědí numerické relativity v rámci 4% chyby [62] .