Численная теория относительности — одна из ветвей общей теории относительности , которая использует численные методы и алгоритмы для решения и анализа проблем. С этой целью суперкомпьютеры часто используются для изучения черных дыр , гравитационных волн , нейтронных звезд и многих других явлений, описанных общей теорией относительности Альберта Эйнштейна . В настоящее время активной областью исследований в области численной теории относительности является моделирование релятивистских двойных звезд и связанных с ними гравитационных волн.
Основной целью численной теории относительности является изучение пространств-времен , точная форма которых неизвестна. Пространства-времена, найденные таким образом вычислительным путем, могут быть полностью динамическими , стационарными или статическими и могут содержать поля материи или вакуум. В случае стационарных и статических решений численные методы также могут использоваться для изучения устойчивости равновесных пространств-времен. В случае динамических пространств-времен проблема может быть разделена на задачу начального значения и эволюцию, каждая из которых требует различных методов.
Численная теория относительности применяется во многих областях, таких как космологические модели , критические явления , возмущенные черные дыры и нейтронные звезды , а также слияние черных дыр и нейтронных звезд, например. В любом из этих случаев уравнения Эйнштейна можно сформулировать несколькими способами, которые позволяют нам развивать динамику. Хотя методы Коши получили большую часть внимания, также использовались характеристические методы и методы, основанные на исчислении Редже . Все эти методы начинаются со снимка гравитационных полей на некоторой гиперповерхности , начальных данных, и развивают эти данные на соседних гиперповерхностях. [1]
Как и во всех проблемах численного анализа, особое внимание уделяется стабильности и сходимости численных решений. В этой линии большое внимание уделяется калибровочным условиям , координатам и различным формулировкам уравнений Эйнштейна и их влиянию на способность производить точные численные решения.
Численные исследования относительности отличаются от работы над классическими теориями поля , поскольку многие методы, реализованные в этих областях, неприменимы в относительности. Однако многие аспекты являются общими с крупномасштабными проблемами в других вычислительных науках, таких как вычислительная гидродинамика , электромагнетизм и механика твердого тела. Численные релятивисты часто работают с прикладными математиками и черпают знания из численного анализа , научных вычислений , уравнений в частных производных и геометрии среди других математических областей специализации.
Альберт Эйнштейн опубликовал свою общую теорию относительности в 1915 году. [2] Она, как и его более ранняя теория специальной теории относительности , описывала пространство и время как единое пространство-время, подчиненное тому, что сейчас известно как уравнения поля Эйнштейна . Они образуют набор связанных нелинейных уравнений в частных производных (УЧП). Спустя более 100 лет с момента первой публикации теории известно относительно немного решений в замкнутой форме для уравнений поля, и из них большинство являются космологическими решениями, которые предполагают особую симметрию для уменьшения сложности уравнений.
Область численной теории относительности возникла из желания построить и изучить более общие решения уравнений поля путем приближенного решения уравнений Эйнштейна численно. Необходимым предшественником таких попыток было разложение пространства-времени обратно на разделенные пространство и время. Это было впервые опубликовано Ричардом Арновиттом , Стэнли Дезером и Чарльзом В. Мизнером в конце 1950-х годов в том, что стало известно как формализм ADM . [3] Хотя по техническим причинам точные уравнения, сформулированные в оригинальной статье ADM, редко используются в численном моделировании, большинство практических подходов к численной теории относительности используют «разложение 3+1» пространства-времени на трехмерное пространство и одномерное время, которое тесно связано с формулировкой ADM, поскольку процедура ADM переформулирует уравнения поля Эйнштейна в ограниченную задачу начального значения , которая может быть решена с использованием вычислительных методологий .
В то время, когда ADM опубликовали свою оригинальную статью, компьютерные технологии не поддерживали численное решение их уравнений для любой задачи сколько-нибудь существенного размера. Первая задокументированная попытка численно решить уравнения поля Эйнштейна, по-видимому, была предпринята SG Hahn и RW Lindquist в 1964 году [4] , за которыми вскоре последовали Larry Smarr [5] [6] и KR Eppley. [7] Эти ранние попытки были сосредоточены на развитии данных Мизнера в осесимметрии (также известной как «2+1 измерения»). Примерно в то же время Цви Пиран написал первый код, который развивал систему с гравитационным излучением, используя цилиндрическую симметрию. [8] В этом расчете Пиран заложил основу для многих концепций, используемых сегодня при развитии уравнений ADM, таких как «свободная эволюция» против «ограниченной эволюции», [ необходимо разъяснение ], которые имеют дело с фундаментальной проблемой обработки уравнений ограничений, возникающих в формализме ADM. Применение симметрии снизило требования к вычислительным ресурсам и памяти, связанные с этой задачей, что позволило исследователям получать результаты на суперкомпьютерах, доступных в то время.
Первые реалистичные расчеты вращающегося коллапса были выполнены в начале восьмидесятых Ричардом Старком и Цви Пираном [9], в которых впервые были рассчитаны формы гравитационных волн, возникающие в результате образования вращающейся черной дыры. В течение почти 20 лет после первых результатов было довольно мало других опубликованных результатов в численной теории относительности, вероятно, из-за отсутствия достаточно мощных компьютеров для решения этой проблемы. В конце 1990-х годов Альянс по Большому Вызову Двойной Черной Дыры успешно смоделировал лобовое столкновение двойной черной дыры . В качестве этапа постобработки группа вычислила горизонт событий для пространства-времени. Этот результат все еще требовал наложения и использования осесимметрии в расчетах. [10]
Некоторые из первых задокументированных попыток решить уравнения Эйнштейна в трех измерениях были сосредоточены на одной черной дыре Шварцшильда , которая описывается статическим и сферически симметричным решением уравнений поля Эйнштейна. Это обеспечивает превосходный тестовый случай в численной теории относительности, поскольку у нее есть решение в замкнутой форме, так что численные результаты можно сравнить с точным решением, поскольку оно статично и поскольку оно содержит одну из самых численно сложных особенностей теории относительности — физическую сингулярность . Одной из первых групп, попытавшихся смоделировать это решение, были Питер Аннинос и др. в 1995 году. [11] В своей статье они указывают, что
В последующие годы не только компьютеры стали мощнее, но и различные исследовательские группы разработали альтернативные методы для повышения эффективности вычислений. Что касается конкретно моделирования черных дыр, были разработаны два метода, чтобы избежать проблем, связанных с существованием физических сингулярностей в решениях уравнений: (1) вырезание и (2) метод «прокола». Кроме того, группа Лазаруса разработала методы для использования ранних результатов кратковременного моделирования, решающего нелинейные уравнения ADM, чтобы предоставить начальные данные для более стабильного кода на основе линеаризованных уравнений, полученных из теории возмущений . В более общем плане, методы адаптивного уточнения сетки , уже используемые в вычислительной гидродинамике , были введены в область численной теории относительности.
В технике вырезания, которая была впервые предложена в конце 1990-х годов [12], часть пространства-времени внутри горизонта событий, окружающего сингулярность черной дыры, просто не эволюционирует. Теоретически это не должно влиять на решение уравнений за пределами горизонта событий из-за принципа причинности и свойств горизонта событий (т. е. ничто физическое внутри черной дыры не может влиять на физику за пределами горизонта). Таким образом, если просто не решать уравнения внутри горизонта, то все равно можно будет получить действительные решения снаружи. Можно «вырезать» внутреннюю часть, налагая входящие граничные условия на границу, окружающую сингулярность, но внутри горизонта. Хотя реализация вырезания была очень успешной, у техники есть две небольшие проблемы. Первая заключается в том, что нужно быть осторожным с координатными условиями. Хотя физические эффекты не могут распространяться изнутри наружу, координатные эффекты могут. Например, если координатные условия были эллиптическими, изменения координат внутри могли бы мгновенно распространяться через горизонт. Это означает, что для распространения эффектов координат необходимы гиперболические координатные условия с характерными скоростями, меньшими, чем у света (например, с использованием гармонических координатных условий). Вторая проблема заключается в том, что по мере движения черных дыр необходимо постоянно корректировать местоположение области вырезания, чтобы она двигалась вместе с черной дырой.
Методика вырезания разрабатывалась в течение нескольких лет, включая разработку новых калибровочных условий, которые увеличили стабильность, и работу, которая продемонстрировала способность областей вырезания перемещаться по вычислительной сетке. [13] [14] [15] [16] [17] [18] Первая стабильная, долгосрочная эволюция орбиты и слияние двух черных дыр с использованием этой техники были опубликованы в 2005 году. [19]
В методе проколов решение разлагается на аналитическую часть [20] , которая содержит сингулярность черной дыры, и численно построенную часть, которая затем свободна от сингулярности. Это обобщение предписания Брилла-Линдквиста [21] для начальных данных черных дыр в состоянии покоя и может быть обобщено до предписания Боуэна-Йорка [22] для начальных данных вращающихся и движущихся черных дыр. До 2005 года все опубликованные случаи использования метода проколов требовали, чтобы координатное положение всех проколов оставалось фиксированным в ходе моделирования. Конечно, черные дыры, находящиеся поблизости друг от друга, будут стремиться двигаться под действием силы тяжести, поэтому тот факт, что координатное положение прокола оставалось фиксированным, означал, что сами системы координат становились «растянутыми» или «скрученными», и это обычно приводило к численным нестабильностям на каком-то этапе моделирования.
В 2005 году группа исследователей впервые продемонстрировала возможность перемещения проколов через систему координат, тем самым устранив некоторые из более ранних проблем с методом. Это позволило получить точные долгосрочные эволюции черных дыр. [19] [23] [24] Выбрав соответствующие координатные условия и сделав грубые аналитические предположения о полях вблизи сингулярности (поскольку никакие физические эффекты не могут распространяться из черной дыры, грубость приближений не имеет значения), можно было получить численные решения для задачи о двух черных дырах, вращающихся друг вокруг друга, а также точно вычислить гравитационное излучение (рябь в пространстве-времени), испускаемое ими. 2005 год был переименован в « annus mirabilis » численной теории относительности, через 100 лет после статей annus mirabilis по специальной теории относительности (1905).
Проект Lazarus (1998–2005) был разработан как пост-Grand Challenge метод извлечения астрофизических результатов из кратковременных полных численных симуляций двойных черных дыр. Он объединил методы аппроксимации до (пост-ньютоновские траектории) и после (возмущения одиночных черных дыр) с полными численными симуляциями, пытающимися решить уравнения поля Эйнштейна. [25] Все предыдущие попытки численно интегрировать в суперкомпьютерах уравнения Гильберта-Эйнштейна, описывающие гравитационное поле вокруг двойных черных дыр, приводили к сбою программного обеспечения до завершения одной орбиты.
Подход проекта Lazarus, тем временем, дал лучшее понимание проблемы бинарных черных дыр и выдал многочисленные и относительно точные результаты, такие как излучаемая энергия и угловой момент, испускаемые в последнем состоянии слияния, [26] [27] линейный импульс, излучаемый дырами неравной массы, [28] и конечная масса и спин оставшейся черной дыры. [29] Метод также вычислил подробные гравитационные волны, испускаемые в процессе слияния, и предсказал, что столкновение черных дыр является самым энергичным единичным событием во Вселенной, высвобождающим больше энергии за долю секунды в форме гравитационного излучения, чем целая галактика за свою жизнь.
Адаптивное уточнение сетки (AMR) как численный метод имеет корни, которые выходят далеко за рамки его первого применения в области численной теории относительности. Уточнение сетки впервые появляется в литературе по численной теории относительности в 1980-х годах в работе Чоптуика в его исследованиях критического коллапса скалярных полей . [30] [31] Первоначальная работа была в одном измерении, но впоследствии она была расширена до двух измерений. [32] В двух измерениях AMR также применялся для изучения неоднородных космологий , [33] [34] и для изучения черных дыр Шварцшильда . [35] Этот метод теперь стал стандартным инструментом в численной теории относительности и использовался для изучения слияния черных дыр и других компактных объектов в дополнение к распространению гравитационного излучения, генерируемого такими астрономическими событиями. [36] [37]
За последние несколько лет [ когда? ] были опубликованы сотни исследовательских работ, приведших к широкому спектру результатов математической теории относительности, гравитационных волн и астрофизики для проблемы орбитальной черной дыры. Этот метод был распространен на астрофизические двойные системы, включающие нейтронные звезды и черные дыры, [38] и множественные черные дыры. [39] Одним из самых удивительных предсказаний является то, что слияние двух черных дыр может придать остаточной дыре скорость до 4000 км/с, что может позволить ей вырваться из любой известной галактики. [40] [41] Моделирование также предсказывает огромное высвобождение гравитационной энергии в этом процессе слияния, составляющее до 8% от ее общей массы покоя. [42]