Собственное время полета Метрика Шварцшильда конечна

Изучая метрические геодезические Шварцшильда, можно легко прийти к следующему дифференциальному уравнению

г р г т "=" С 2 ( 1 2 г М р )
которая связывает радиальную координату и собственное время за пределами горизонта событий р ЧАС "=" 2 г М (использую, конечно, с "=" 1 ).
Представьте, что наблюдатель падает из исходного положения р 0 "=" 2 г М + час час > 0 для того, чтобы быть вне горизонта событий) вплоть до горизонта событий р ЧАС "=" 2 г М . Собственное время полета, благодаря приведенному выше дифференциальному уравнению, равно
Δ т "=" т 0 т ЧАС г т "=" р 0 "=" 2 г М + час р ЧАС "=" 2 г М г р С 2 ( 1 2 г М р )
Однако координатное время полета определяется соотношением г т "=" С ( 1 2 г М / р ) г т (где С > 0 , см., например, http://gfm.cii.fc.ul.pt/events/lecture_series/general_relativity/gfm-general_relativity-lecture4.pdf ), в результате
Δ т "=" т 0 т ЧАС г т "=" р 0 "=" 2 г М + час р ЧАС "=" 2 г М С 1 2 г М р г р С 2 ( 1 2 г М р )
Отсюда (см., например, приведенную выше ссылку) говорится, что « непосредственным вычислением интегралов легко видеть, что собственное время полета конечно, а координатное время полета бесконечно, орбита наблюдатель никогда не увидит, как падающий наблюдатель достигает горизонта событий, за исключением асимптотического ''.

Я попытался вычислить первый интеграл с помощью Mathematica, и результат, который я получил, абсурден, что дает сложное собственное время. Вот код, который я использовал:
Integrate[1/Sqrt[C^2 - (1 - 2*G*M/r)], {r, 2*G*M, 2*G*M + h}]
И вывод:
ConditionalExpression[(G M (2 Sqrt[C^2 - C^4] - I Log[2] - I Log[((-I + 2 I C^2 + 2 Sqrt[C^2 - C^4]) G M)/Sqrt[1 - C^2]]))/(1 - C^2)^(3/2) + (-Sqrt[1 - C^2] (h + 2 G M) Sqrt[C^2 - h/(h + 2 G M)] + I G M Log[2] + I G M Log[(-I h + I C^2 h - I G M + 2 I C^2 G M + Sqrt[1 - C^2] h Sqrt[C^2 - h/(h + 2 G M)] + 2 Sqrt[1 - C^2] G M Sqrt[C^2 - h/(h + 2 G M)])/Sqrt[1 - C^2]])/(1 - C^2)^(3/2), ((G M)/h != 0 && Re[(GM)/h] >= 0) ||Re[(G M)/h] < -(1/2) || (G M)/h \[NotElement] Reals]

Что мне здесь не хватает? Спасибо

Лучший способ приблизиться к интегралам - выполнить несколько приятных замен, чтобы максимально упростить подынтегральные выражения (как предложено во втором ответе), оставляя так много переменных, как правило, искажает CAS, потому что он должен делать предположения об этих величинах. Он не знает, что M является массой и должно быть положительным, например, как вы можете видеть в условном выражении. Один из способов помочь в этом — установить допущения, используя Refine as Refine[выражение, Допущения->{...}], где фигурные скобки содержат все допущения, необходимые для очистки выражения.
Аналитический результат приведен здесь: physics.stackexchange.com/questions/426143/427025#427025 . Второй график также показывает, что т (синим цветом) расходится на горизонте, а т (зеленым цветом) нет.
@Triatticus Спасибо! Я попытался с вашим ответом, но полученный результат, хотя он действительно реален и имеет смысл, не согласуется с ожидаемым значением.
@safesphere Спасибо! Результаты, которые вы даете, верны, это результаты, которые я пытался проверить с помощью Mathematica. Однако я до сих пор не знаю, почему команда не работает.
Здесь есть вопросы по физике?

Ответы (2)

Я тоже не совсем понимаю, как это "легко увидеть" из того.
Однако вы можете решить интеграл для собственного времени, сначала заметив, что для С 2 "=" константа у нас (с р с "=" 2 г М )

С 2 "=" 1 р с р 0

что следует из включения ( 4.2 ) в ( 4.1 ) и рассматривая падающую частицу, которая покоится в р "=" р 0 , значение р ˙ | р "=" р 0 "=" 0 .
Это вы можете использовать, чтобы написать свое первое уравнение как

р ˙ "=" р с 1 / 2 ( р 0 р р 0 р ) 1 / 2

Отсюда интеграл по собственному времени будет равен

Δ т "=" р с 1 / 2 р 0 р с ( р р 0 р 0 р ) 1 / 2 г р

которые вы можете решить с помощью Mathematica или также вводя параметризацию

р ( η ) "=" р 0 2 ( 1 + потому что η )

с η е [ 0 , π ] .

Спасибо! Да, вы знаете, научные статьи — это множество слов, таких как «тривиальный», «очевидный», «очень ясный», ... в моменты, когда это даже близко не может быть тривиальным или простым, ха-ха. Спасибо за ваш ответ, но поместите команду в Mathematica с вашим выражением для Δ т программа не дает результата, я не знаю, что происходит.
Вы тоже пробовали использовать параметризацию? При этом это на самом деле довольно просто и может быть сделано с ручкой и бумагой. У меня нет доступа к Mathematica, поэтому, к сожалению, я не могу помочь вам с программной частью.
Да, я сделал это и с параметризацией, и она не вычисляет результат. Не волнуйся! Спасибо за Ваш ответ. Я знаю, что это можно сделать аналитически, но я хотел знать, что происходит с программным обеспечением.

результат, который я получаю, абсурден, давая комплексу правильное время.

В вашем результате есть I, что я предполагаю, это обозначение Mathematica для я "=" 1 , но это не означает, что результат сложный. Для некоторых методов интеграции характерно получение результатов, которые выглядят сложными, но на самом деле оказываются реальными, когда вы их оцениваете.

Вы излишне усложняете свои результаты, принимая г и М как параметры. При выполнении такого рода интеграла нужно преобразовать как можно больше переменных в безразмерные величины, что обычно приводит к устранению ненужных констант, подобных этим. Здесь вы хотите перейти к переменной Икс "=" р / 2 г М . Не вводите его в систему компьютерной алгебры (CAS), не выполнив эту подготовку. CAS глупы и производят сложный вывод, поэтому дайте ему шанс побороться, выполнив ту же настройку, которую вы обычно делаете, работая без CAS.

Я использую CAS maxima с открытым исходным кодом, а не mathematica. У него нет проблем с получением результата без явного я в нем.

$ maxima -q --batch-string="assume(h>0 and (c^2-1)<0 and c>0); integrate(-(c^2-(1-1/x))^(-1/2),x,1+h,1);"

(%i1) assume(h > 0 and c^2-1 < 0 and c > 0)
                                     2
(%o1)                       [h > 0, c  < 1, c > 0]
(%i2) integrate(-(c^2-(1-1/r))^((-1)/2),r,1+h,1)
                                   2                     2           2
                 2       sqrt(1 - c ) sqrt(h + 1) sqrt((c  - 1) h + c )
(%o2) (sqrt(1 - c ) atan(----------------------------------------------)
                                        2           2
                                      (c  - 1) h + c  - 1
     2                         2           2     4      2
 + (c  - 1) sqrt(h + 1) sqrt((c  - 1) h + c ))/(c  - 2 c  + 1)
                                 2
             2       c sqrt(1 - c )     3
   sqrt(1 - c ) atan(--------------) + c  - c
                          2
                         c  - 1
 - ------------------------------------------
                  4      2
                 c  - 2 c  + 1

Возможно, вы захотите попробовать настроить его таким же образом в mathematica. Безусловно, после замены переменных интерпретировать будет проще и легче.

На самом деле нет никакой гарантии, что этот результат окажется реальным только потому, что нет я в нем. В нем действительно есть один квадратный корень, который выглядит так, как будто он может быть мнимым. То, как мы знаем, что оно должно быть реальным, состоит в том, что это определенный интеграл от действительного подынтегрального выражения. Если вы хотите проверить, действительно ли это правда, попробуйте просто использовать случайные числа для С и час . Если результат окажется реальным, то это не совпадение. Если получается сложно, то может иметь место одно из следующего: (1) вы закодировали его неправильно, (2) есть ошибка в CAS или (3) есть проблема с ответвлениями для арктангенса.

Я пытался использовать вашу замену в безразмерных единицах, но результат тот же, что и выше, и все еще имеет ужасное явное я "=" я "=" 1 . Я не понимаю, что происходит, может быть, я пытаюсь использовать Максиму или другой программный продукт и проверить это там. Спасибо!