A numerikus relativitáselmélet az általános relativitáselmélet egyik területe , amely numerikus módszereket és algoritmusokat fejleszt és használ fizikai folyamatok számítógépes szimulációjára erős gravitációs mezőkben, amikor az Einstein-egyenleteket numerikusan kell megoldani . A főbb fizikai rendszerek , amelyek leírásához numerikus relativitáselmélet szükséges, a relativisztikus asztrofizikához kapcsolódnak, ide tartoznak a gravitációs összeomlás , a neutroncsillagok , a fekete lyukak , a gravitációs hullámok és más objektumok és jelenségek, amelyek megfelelő leírásához a teljes általános elméletre kell hivatkozni. relativitáselmélet konvencionális közelítések nélkül gyenge mezők és alacsony sebességek (mint a Newtoni utáni terjeszkedésben és a perturbációelméletben az Einstein -egyenletek egzakt megoldásainak hátterében ) [1] .
Ezen a területen a modellezés speciális numerikus módszereket igényel az Einstein-egyenletek összetettsége és nemlinearitása miatt (például az egyenletek ábrázolásától függ az időbeli alakulásuk Cauchy-probléma megfogalmazásának hiperbolicitása és helyessége, valamint a kezdeti és peremfeltételek [2] ), valamint a legtöbb háromdimenziós feladathoz is - a nagy számítási teljesítmény csak a modern szuperszámítógépek számára elérhető . Jelenleg a numerikus relativitáselméletben a relativisztikus közeli kettős csillagok és a hozzájuk kapcsolódó gravitációs hullámok modellezése, valamint számos más matematikai és asztrofizikai probléma [1] területén relevánsak a kutatások .
A numerikus relativitáselmélet fő célja a gravitációs mezők tanulmányozása , amelyek pontos analitikai formája ismeretlen. A gravitációs mezők, amelyek formáját számítással keresik, lehetnek teljesen dinamikusak , stacionerek vagy statikusak, és tartalmazhatnak anyagmezőket [~ 1] vagy vákuumokat is. Stacionárius és statikus megoldások esetén numerikus módszerekkel vizsgálható ezen konfigurációk stabilitása. A dinamikus gravitációs mezők esetében viszont a probléma két részre osztható, amelyek különböző megoldási módszereket igényelnek: a kezdeti értékek problémájára és az evolúció problémájára [3] .
A numerikus relativitáselméletet a kozmológiai modellek , a gravitációs összeomlás kritikus jelenségei , valamint a fekete lyukakat és a neutroncsillagokat érintő folyamatok , különösen azok egyesülésének és perturbációjának tanulmányozására használják . Mindegyik esetben nyomon kell követni a téridő alakulását, amelyre az Einstein-egyenletek többféleképpen is ábrázolhatók. A legnépszerűbbek a Cauchy-probléma módszerei , de alkalmazzák a karakterisztikás módszert [4] és a Regge-számításon [5] alapuló módszereket is . A fenti módszerek mindegyike a gravitációs mező "pillanatfelvételével" kezdődik valamilyen hiperfelületen , azaz a kezdeti adatokból, majd az időben előrehaladva nyomon követi annak alakulását a következő közeli hiperfelületekre [6] .
Mint a numerikus elemzés minden problémájában, a numerikus relativitáselméletben is nagy figyelmet fordítanak a numerikus megoldások stabilitására és konvergenciájára, a megengedett kezdeti és peremfeltételekre. A numerikus relativitáselmélet sajátosságai a szelvény- és koordinátafeltételek jelenléte okozta bonyodalmak , valamint az Einstein-egyenletek különféle reprezentációi és azok befolyása a pontos numerikus megoldások megszerzésére.
A klasszikus térelméletben használt numerikus technikák közül sok nem alkalmazható az általános relativitáselméletben, ezért az ezen a területen végzett munka különbözik a numerikus relativitáselmélet területén végzett kutatásoktól. A nagy léptékű problémákban azonban a numerikus relativitáselmélet számos vonatkozásban osztozik más számítástechnikai tudományokkal, mint például a számítási folyadékdinamika , az elektrodinamika és a merev test mechanikája . A numerikus relativitáselméletben részt vevő tudósok gyakran dolgoznak együtt alkalmazott matematikusokkal, és kapcsolatba kerülnek a matematika olyan területeivel, mint a numerikus elemzés , a párhuzamos számítástechnika , a parciális differenciálegyenletek és a geometria [7] .
Albert Einstein 1915-ben publikálta az általános relativitáselmélet végső változatát [8] . Ez az elmélet, akárcsak az őt megelőző speciális relativitáselmélet , a teret és az időt egyetlen tárgyként írja le - téridőként , amelynek fejlődése engedelmeskedik Einstein egyenleteinek . Összekapcsolt nemlineáris parciális differenciálegyenletek rendszerét alkotják . Az egyenletek levezetése óta eltelt évszázadban ezeknek az egyenleteknek csak viszonylag kevés fizikailag releváns , egzakt analitikai megoldása vált ismertté , és ezek többsége a nagy szimmetria feltételezésével származtatott, ami leegyszerűsíti az egyenletek megoldását. , mint például a Friedmann-féle megoldások egy homogén és izotróp Univerzumra [9] .
A numerikus relativitáselmélet területe abból a vágyból alakult ki, hogy az Einstein-egyenletek általánosabb és fizikailag alkalmazható megoldásait numerikus közelítő megoldással tanulmányozzák. Egy ilyen megoldás szükséges feltétele volt, hogy egyetlen négydimenziós téridőt osztott háromdimenziós térre és egydimenziós időre, az úgynevezett 3 + 1 hasítást hajtsák végre . Sőt, sokféleképpen végrehajtható, ami jelentősen bonyolíthatja vagy leegyszerűsítheti a kapott egyenletek integrálásának problémáját. Az első meglehetősen sikeres szétválási kísérletet Richard Arnowitt, Stanley Deser és Charles Misner tette az 1950-es évek végén a hamiltoni formalizmusban a Dirac által jelzett úton . Ez az úgynevezett ADM formalizmust, az Arnowitt-Deser-Mizner formalizmust [10] alkotó egyenletek megszerzésében tetőzött . Bár technikai okok miatt ezek az egyenletek nem bizonyultak túl kényelmesnek a numerikus integrációhoz - csak gyengén hiperbolikusak , ezért ritkán használják valós számításokban - a numerikus relativitáselmélet gyakorlati megközelítéseinek túlnyomó többsége a 3 + 1-es felosztást alkalmazza, közel az ADM formalizmusban használt . Egy ilyen felosztás az Einstein-egyenletek újrafogalmazásához vezet egy Cauchy-probléma formájában, a kezdeti értékek korlátozásával, ami már számítógépeken is numerikus megoldásra alkalmas [11] .
A téridő koordinátái nem határozhatók meg egyértelműen, ezért a kezdeti hiperfelületen koordináták rögzítésekor is, amikor egy szomszédos hiperfelületre mozogunk, az idő- és térkoordináták különböző pontokon eltérő módon „tolhatók” (már a speciális elméletben). relativitáselmélet, az idő áramlásának iránya és sebessége nem esik egybe a különböző tehetetlenségi vonatkoztatási rendszerekben), ami a numerikus relativitáselmélet sajátossága. Ez a mérőszabadság - amely nem befolyásolja a fizikai folyamatokat, hanem csak megváltoztatja azok leírását a koordinátákban és ennek megfelelően a megoldandó egyenletekben - a mozgás- és eltolásfüggvények megválasztásának önkényességében nyilvánul meg , fix pontokkal "lökve" a pontokat. térbeli koordináták a kezdetitől a szomszédos hiperfelületig időben előre — , illetve oldalirányban térben — , ill. Ezen függvények megválasztásának lehetősége potenciális előnyt jelent az egyenletek numerikus megoldása szempontjából, de számos „természetes” választása ezeknek a koordináta- vagy szelvényfeltételeknek a megoldásokban numerikus instabilitásokat okoz, ami szimulációs szünetekhez vezet [12] .
Az ADM formalizmussal foglalkozó eredeti közlemények megjelenésekor a számítástechnika fejlődése nem tette lehetővé az egyenletek felhasználásával történő számításokat bármilyen ésszerű méretű problémára. Történelmileg az Einstein-egyenletek numerikus megoldására az első kísérletet Hahn és Lindqvist tette 1964-ben [13] , majd az 1970-es években Smarr [14] [15] és Eppley [16] . Ezek a korai próbálkozások Misner kiindulási adatainak tengelyirányban szimmetrikus terekben (más néven "2+1 dimenziók") való fejlődéséhez kapcsolódtak. Körülbelül ugyanebben az időben Zvi Piran megírta az első kódot, amely egy hengerszimmetrikus, gravitációs sugárzást kibocsátó rendszer fejlődését követte nyomon [17] . Ebben a fejlesztésben Piran számos, a numerikus relativitáselméletben manapság használt fogalmat kezdeményezett, mint például a szabad evolúció és a kényszerű evolúció, olyan módszereket, amelyek különböző módon közelítik meg a kezdeti adatkorlátok időbeni alakulásának problémáját [18] [19] . A szimmetria alkalmazása csökkentette a memória és a számítási teljesítmény szükségességét, így a tudósok az akkori szuperszámítógépeket használhatták a probléma megoldására [17] .
A valódi asztrofizikai problémára, a forgó összeomlásra az első reális számításokat az 1980-as évek elején végezték Richard Stark és Zvi Piran [20] , amelyben először számították ki a formálódó forgó fekete lyuk által kibocsátott gravitációs hullámokat. A publikáció óta eltelt közel két évtizedben a numerikus relativitáselméletben csak néhány új eredményt hoztak nyilvánosságra, valószínűleg azért, mert hiányoztak a problémák megoldásához elég erős számítógépek. Az 1990-es években a Binary Black Hole Grand Challenge Alliance sikeresen szimulált két fekete lyuk frontális ütközését a probléma tengelyirányú szimmetriájából adódó egyszerűsítések segítségével . Az utófeldolgozási szakaszban a csoport ki tudta számítani az eredményül kapott megoldás eseményhorizontját [21] .
Az Einstein-egyenletek teljes 3D-s térgeometriában való numerikus megoldására tett első ismert kísérletek egy nem forgó Schwarzschild-fekete lyukra összpontosítottak , amely az Einstein-egyenletek statikus és gömbszimmetrikus megoldása. Kiváló teszt a numerikus relativitáselmélet vizsgálatára, mivel egyrészt a megoldás egzakt analitikai formában ismert, amellyel a numerikus eredmények összehasonlíthatók, másrészt statikus, és minden nem forgó fekete lyuknak idővel konvergálnia kell hozzá. , harmadszor pedig a numerikus modellezés egyik legnehezebb objektumát tartalmazza - egy fizikai gravitációs szingularitást a közepén. Az egyik első kísérletet ennek a megoldásnak a számszerű megszerzésére Anninos és munkatársai tettek 1995-ben [22] . Ebben a munkában megjegyezték:
A 3D numerikus relativitáselmélet fejlődését részben hátráltatja az elegendő memóriával és számítási teljesítménnyel rendelkező számítógépek hiánya, amelyek jó felbontású 3D-s számításokat végeznének.
Eredeti szöveg (angol)[ showelrejt] A háromdimenziós numerikus relativitáselmélet fejlődését részben hátráltatta a megfelelő memóriával és számítási teljesítménnyel rendelkező számítógépek hiánya a 3D-s téridő-számítások elvégzéséhez.Az évek során a számítógépek erősebbé válása mellett a különböző kutatócsoportok alternatív technikákat fejlesztettek ki a számítástechnika hatékonyságának növelésére. Először is, a Lazarus -csoport olyan módszereket dolgozott ki, amelyek a fekete lyukak összeolvadására vonatkozó nemlineáris ADM-egyenleteket megoldó rövid szimulációk korai eredményeit használták fel, hogy kiindulási adatokat szolgáltassanak egy robusztusabb kódhoz, amely az egyetlen fekete lyuk perturbáció elméletének linearizált egyenletein alapul [23]. . Ezután a fekete lyuk modellezéssel kapcsolatban két technikát fejlesztettek ki, hogy elkerüljék a fizikai szingularitás meglétével kapcsolatos problémákat az egyenletek megoldásában: (1) az eliminációt és (2) a „prick” módszert [24] . Ezeknek a módszereknek a kombinációja a talált megfelelő koordináta-viszonyokkal 2005-ben áttörést hozott a bináris fekete lyukak modellezésében, amely Pretorius munkájával kezdődött [25] . Néhány évvel később az új módszerek numerikus stabilitása lehetővé tette a bináris fekete lyukak szinte tetszőleges konfigurációinak szimulálását, amelyek egymás körül tíz és száz fordulatot írnak le az egyesülés előtt. Emellett a numerikus relativitáselméletben elkezdték alkalmazni a számítási rács adaptív finomításának módszereit, amelyeket korábban a számítási folyadékdinamikában alkalmaztak [26] .
Project LazarusA Lazarus projektet (1998–2005) a Nagy Kihívás után fejlesztették ki, hogy asztrofizikailag releváns eredményeket nyerjenek ki a bináris fekete lyukak egyesülésének akkoriban elérhető rövid numerikus szimulációiból. Abban az időben az Einstein-féle bináris fekete lyukak tér-idő-egyenleteinek szuperszámítógépeken történő integrálására tett kísérletek – különféle instabilitások miatt – nem tudtak előrehaladni a rendszer egy teljes körforgásának befejezéséig. A projekt keretein belül a kutatók közelítő módszereket kombináltak egy lyukpár eggyé alakítása előtt ( poszt-newtoni trajektóriák ) és utána (egyes fekete lyukak zavarásai) magának a folyamatnak a teljes numerikus megoldásaival [23] .
Abban az időben a Lazarus-projekt megközelítése volt a legjobb megközelítés a bináris fekete lyukak problémájára, és számos, asztrofizikai alkalmazásokhoz kellően pontos eredményt adott, például a gravitációs hullámok által elszállított energia és szögimpulzus értékeit [27] ] [28] , valamint a különböző tömegű fekete lyukak egyesülésének lendülete [29] , valamint a kialakuló fekete lyuk végső tömegének, impulzusának és szögimpulzusának értékei [30] . A projekt módszerei lehetővé tették az egyesülés során kibocsátott gravitációs hullámok részletes formáinak kiszámítását is - ami fontos volt a gravitációs teleszkópok számára , és azt jósolták, hogy a fekete lyukak ütközését az univerzum legerősebb energiakitörései kísérhetik, amikor több energiát kell elérni. a másodperc töredéke alatt szabadul fel gravitációs sugárzás formájában, mint amennyit a galaxis összes csillaga kisugároz létezése során – a gravitációs sugárzás a rendszer kezdeti csökkentett tömegének több százalékát is elviszi [31] .
Kizárási módszerA kivágási technikában , amelyet először az 1990-es évek végén javasoltak [32] , a téridőnek az eseményhorizonton belüli , a fekete lyuk szingularitást körülvevő részét egyszerűen kizárják az evolúcióból . Ez elméletileg nem befolyásolhatja az eseményhorizonton kívüli döntést az ok-okozati összefüggés elve és a horizont tulajdonságai miatt - hiszen a horizont alatti fizikai kölcsönhatások nem lehetnek hatással a fizikára azon kívül. Így, ha egyszerűen nem oldja meg az egyenleteket a fekete lyukon belül, akkor is megkaphatja a pontos valós megoldást azon kívül. A belső dinamika „kiküszöbölése” lehetséges úgy, hogy a horizonton belüli határra kényszerítjük a szingularitást, a kimenő hullámok hiányának peremfeltételeit [33] .
Bár az eliminációs technika alkalmazása nagyon sikeres volt, két apró probléma van vele. Az első az, hogy körültekintően kell kiválasztani és használni a koordináta-feltételeket. Míg a fizikai hatások nem terjedhetnek ki a horizonton belülről kifelé, a koordinált hatások igen. Például, ha elliptikus koordináta-feltételeket írnak elő, a fekete lyukon belüli rácsváltozások azonnal kifelé terjedhetnek a horizonton [34] . Ez azt jelenti, hogy az eliminációs módszer alkalmazásához hiperbolikus típusú koordinátafeltételek alkalmazása szükséges, amelyekben a koordinátahatások jellemző terjedési sebessége kisebb vagy egyenlő, mint a fénysebesség (például harmonikus koordináta-feltételek alkalmazása esetén). [35] . A második probléma az, hogy mivel a fekete lyuk mozog, a kizárási tartománynak folyamatosan együtt kell mozognia vele [33] .
Az eliminációs módszert több éve fejlesztették, miközben új kalibrációs feltételeket találtak, amelyek növelik a megoldási eljárás stabilitását, és igazolták a kizárt régiók mozgási képességét a számítási rács mentén [36] [37] [38] [ 39] [40] [35] . Az első stabilan hosszú számítás két fekete lyuk pályájáról és egyesüléséről ezzel a technikával 2005-ben jelent meg [25] .
Az injekciós módszerA punkciós módszerben a megoldás egy analitikus részre [41] van felosztva , amely tartalmazza a fekete lyuk szingularitást - a punkciót, és egy numerikusan felépített részre, amely nem tartalmazza a szingularitást. Ez a módszer a Brill–Lindquist algoritmus [42] általánosítása nyugalmi fekete lyukakat tartalmazó kiindulási adatokra, és tovább általánosítható a Bowen–York algoritmusra [43] forgó és mozgó fekete lyukakat tartalmazó kezdeti adatok esetében. 2005-ig a "prick" módszer használatára vonatkozó összes publikált példa megkövetelte, hogy az összes szúrás koordinátáit rögzítsék a szimuláció teljes időtartama alatt. Természetesen az egymáshoz közel lévő fekete lyukak gravitációs erők hatására elmozdulnak, így a csapok rögzített koordinátái azt jelentik, hogy a koordinátarendszerek "megnyúlnak" vagy "eltorzulnak", ami bizonyos szakaszokban numerikus instabilitáshoz vezet. a szimuláció. Hasonló hatásokat okoz egy másik módszer alkalmazása - a szingularitások elkerülése, amikor a szimulációban az anyag összeomlásával fekete lyukak keletkeznek, és a koordinátafeltételeket úgy választják meg, hogy az időben fejlődő háromdimenziós hiperfelület ne érje el. a szingularitást a számítások végéig, körülötte megnyúlt "szarvat" képezve [44] .
A kutatók 2005-ben mutatták be először annak lehetőségét, hogy a tűszúrások egy koordinátarendszer mentén mozogjanak, ezzel megoldva a módszer néhány korai problémáját, ami lehetővé tette a fekete lyukak hosszú távú fejlődésének pontos nyomon követését [25] [45 ] ] [46] . Megfelelő koordinátafeltételek kiválasztásával és a szingularitás közelében lévő fizikai mezők durva analitikai közelítésével (mivel a fekete lyukból fizikai hatások nem szökhetnek ki, a közelítés érdessége nem lényeges), numerikus megoldásokat kaphatunk két fekete lyuk problémájára. egymás körül forognak, és a gravitációs sugárzásukat is pontosan kiszámítják [47] .
Adaptív háló finomításAz adaptív hálófinomítást mint numerikus módszert már jóval a numerikus relativitáselmélet megjelenése előtt használták a fizikában. Ebben használták először az 1980-as években Choptwick munkáiban, amikor a skaláris mező összeomlása során bekövetkező kritikus jelenségeket tanulmányozták , amikor a mező konfigurációi a fekete lyuk végső kialakulása és a tér végső tágulása között vannak. [48] [49] . Az eredeti munka egydimenziós volt, mivel gömbszimmetriát alkalmaztak, de aztán a módszert kétdimenziósra általánosították [50] . Kétdimenziós finomítási módszereket is alkalmaztak az inhomogén kozmológiák [51] [52] és a Schwarzschild fekete lyukak [53] vizsgálatára . Az adaptív finomítási módszerek mára a numerikus relativitáselmélet szabványos eszközévé váltak, ésaz ilyen események által keltett gravitációs hullámok terjedésének tanulmányozása mellett a fekete lyukak egyesülésének és más kompakt objektumok tanulmányozásában is használatosak [54] [55] .
A numerikus relativitáselméletről a mai napig több tucat és több száz dolgozat születik, amelyek az általános relativitáselmélet matematika, a gravitációs hullámok és az asztrofizika területén mutatják be az egymás körül forgó fekete lyukak problémájának megoldásával kapott eredmények széles skáláját. Az alkalmazott módszereket asztrofizikai kettős rendszerek, köztük neutroncsillagok, fekete lyukak [56] és fekete lyukhalmazok [57] vizsgálatára általánosították . Ezek a papírok többek között azt jósolják, hogy ha két forgó fekete lyuk egyesül, a létrejövő lyuk akár 4000, sőt akár 10 000 km/s sebességet is elérhet , ami lehetővé teszi, hogy bármely ismert galaxison túllépjen [58] [59] . A szimulációk azt is jósolják, hogy az összeolvadás során hatalmas energiafelszabadulás következik be, amely a nyugalmi össztömeg akár 8%-át is elérheti, valamint egy fekete lyuk forgástengelyének éles változásának lehetőségét , ami megmagyarázhatja a sugárirányok változásait. rádiógalaxisokban megfigyelhető [ 60] . Fontos kutatási irány az összeolvadó fekete lyukak gravitációs sugárzási formáinak katalógusának elkészítése is, amely nélkül ezeknek a jeleknek a keresése az olyan detektorok adataiban, mint a LIGO és a VIRGO , sokkal kevésbé érzékeny [61] .
A numerikus relativitáselmélet modern módszereinek pontossága közvetlenül a gravitációs hullámok felfedezése után a gyakorlatban is ellenőrizhetővé vált . A GW150914 jel 4%-os hibán belül megegyezik a numerikus relativitáselmélet előrejelzéseivel [62] .