Tiempo adecuado de vuelo La métrica de Schwarzschild es finita

Con el estudio de las geodésicas métricas de Schwarzschild, uno puede llegar fácilmente a la siguiente ecuación diferencial

d r d τ = C 2 ( 1 2 GRAMO METRO r )
que relaciona la coordenada radial y el tiempo propio fuera del horizonte de sucesos r H = 2 GRAMO METRO (Estoy usando, por supuesto, C = 1 ).
Imagine que un observador está cayendo desde la posición inicial r 0 = 2 GRAMO METRO + h (con h > 0 para estar fuera del horizonte de eventos) hasta el horizonte de eventos r H = 2 GRAMO METRO . El tiempo propio de vuelo es, gracias a la ecuación diferencial dada anteriormente, por
Δ τ = τ 0 τ H d τ = r 0 = 2 GRAMO METRO + h r H = 2 GRAMO METRO d r C 2 ( 1 2 GRAMO METRO r )
Sin embargo, la coordenada del tiempo de vuelo está dada por la relación d t = C ( 1 2 GRAMO METRO / r ) d τ (dónde C > 0 , ver por ejemplo http://gfm.cii.fc.ul.pt/events/lecture_series/general_relativity/gfm-general_relativity-lecture4.pdf ), resultando
Δ t = t 0 t H d t = r 0 = 2 GRAMO METRO + h r H = 2 GRAMO METRO C 1 2 GRAMO METRO r d r C 2 ( 1 2 GRAMO METRO r )
A partir de esto (ver por ejemplo la referencia dada arriba) se dice que " Es fácil ver, por evaluación directa de las integrales, que el tiempo propio de vuelo es finito, mientras que el tiempo de coordenadas de vuelo es infinito, la órbita El observador nunca verá al observador que cae alcanzar el horizonte de eventos, excepto asintótico ''.

He intentado calcular la primera integral usando Mathematica y el resultado que obtengo es absurdo, dando un tiempo propio complejo. Aquí está el código que usé:
Integrate[1/Sqrt[C^2 - (1 - 2*G*M/r)], {r, 2*G*M, 2*G*M + h}]
Y el resultado es:
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]

¿Qué me estoy perdiendo aquí? Gracias

La mejor manera de abordar las integrales es realizar algunas sustituciones agradables para simplificar los integrandos tanto como sea posible (como se sugiere en la segunda respuesta), dejar tantas variables tiende a estropear el CAS porque tiene que hacer suposiciones sobre esas cantidades. No sabe que M es una masa y debe ser positivo, por ejemplo, como puede ver en la expresión condicional. Una forma de ayudar con esto es establecer suposiciones usando Refinar como Refinar[ expresión, Suposiciones->{...}], donde las llaves contienen todas las suposiciones necesarias para limpiar la expresión.
El resultado analítico se proporciona aquí: physics.stackexchange.com/questions/426143/427025#427025 - El segundo gráfico también muestra que t (en azul) diverge en el horizonte mientras τ (en verde) no.
@Triatticus ¡Gracias! Probé con tu respuesta pero el resultado dado, aunque es real y tiene sentido, no concuerda con el valor esperado.
@safesphere ¡Gracias! Los resultados que das son correctos, esos son los resultados que estaba tratando de verificar con Mathematica. Sin embargo, todavía no sé por qué el comando no funciona.
¿Hay una pregunta de física aquí?

Respuestas (2)

Tampoco sé muy bien cómo esto es "fácil de ver" a partir de eso.
Sin embargo, puedes resolver la integral para el tiempo propio si primero observas que para C 2 = constante tenemos (con r s = 2 GRAMO METRO )

C 2 = 1 r s r 0

que se deriva de enchufar ( 4.2 ) en ( 4.1 ) y considerando una partícula que cae y que está en reposo en r = r 0 , significado r ˙ | r = r 0 = 0 .
Esto lo puedes usar para escribir tu primera ecuación como

r ˙ = r s 1 / 2 ( r 0 r r 0 r ) 1 / 2

A partir de esto, la integral de tiempo adecuada resultará ser

Δ τ = r s 1 / 2 r 0 r s ( r r 0 r 0 r ) 1 / 2 d r

que puedes resolver usando Mathematica o también introduciendo la parametrización

r ( η ) = r 0 2 ( 1 + porque η )

con η [ 0 , π ] .

¡Gracias! Sí, ya sabes, los artículos científicos están llenos de palabras como "trivial", "obvio", "muy claro", ... en momentos en los que no está ni cerca de ser trivial o fácil, jaja. Gracias por su respuesta, pero poniendo el comando en Mathematica con su expresión para Δ τ el programa no da resultado, no se que pasa.
¿Has probado también con la parametrización? Con eso, en realidad es bastante sencillo y se puede hacer con lápiz y papel. No tengo acceso a Mathematica y, por lo tanto, lamentablemente no puedo ayudarlo con la parte del software.
Sí, también lo hice con la parametrización y no calcula un resultado. ¡No te preocupes! Gracias por tu respuesta. Sé que se puede hacer analíticamente, pero quería saber qué está pasando con el software.

el resultado que obtengo es absurdo, dando un tiempo propio complejo.

Su resultado tiene Ien él, que supongo que es la notación de Mathematica para i = 1 , pero eso no significa que el resultado sea complejo. Es común con ciertas técnicas de integración obtener resultados que parecen complejos pero que en realidad resultan ser reales cuando los evalúas.

Está complicando innecesariamente sus resultados al tomar GRAMO y METRO como parámetros. Al hacer este tipo de integral, lo que hay que hacer es convertir tantas variables como sea posible en cantidades sin unidades, lo que generalmente tiene el efecto de eliminar constantes innecesarias como estas. Aquí quieres cambiar a la variable X = r / 2 GRAMO METRO . No lo introduzca simplemente en el Sistema de álgebra computacional (CAS) sin hacer esta preparación. Los CAS son estúpidos y producen resultados complicados, así que dale una oportunidad haciendo la misma configuración que harías normalmente si trabajas sin un CAS.

Uso los máximos CAS de código abierto, en lugar de mathematica. No tiene ningún problema en producir un resultado sin explícito i está en eso.

$ 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

Es posible que desee intentar configurarlo de la misma manera en Mathematica. Sin duda, será más simple y más fácil de interpretar después del cambio de variables.

En realidad, no hay garantía de que este resultado sea real solo porque no hay i está en eso. Tiene una raíz cuadrada que parece que podría ser imaginaria. La forma en que sabemos que debería ser real es que es una integral definida de un integrando real. Si luego desea verificar si realmente es real, intente usar algunos números aleatorios para C y h . Si el resultado resulta ser real, entonces eso no es una coincidencia. Si resulta complejo, entonces uno de los siguientes puede ser válido: (1) lo ha codificado incorrectamente, (2) hay un error en el CAS o (3) hay un problema con los cortes de rama para la tangente inversa.

Intenté usar su sustitución en unidades adimensionales, pero el resultado es el mismo que se indicó anteriormente y aún tiene el terrible explícito I = i = 1 . No entiendo que pasa, tal vez trato de usar Maxima u otro software y verifico ahí. ¡Gracias!