Soluciones a ecuaciones geodésicas en el espacio AdS3

Estoy luchando con una tarea de mi profesor, me ha pedido que calcule numéricamente (mathematica) y analíticamente soluciones para la trayectoria de una partícula puntual en el espacio AdS3 con coordenadas globales dadas por métrica

d s 2 = R 2 ( aporrear 2 ( ρ ) d t 2 + d ρ 2 + pecado 2 ( ρ ) d ϕ 2 )
Calculé tres ecuaciones geodésicas de movimiento.

ρ ¨ + pecado ( ρ ) aporrear ( ρ ) ( t ˙ t ˙ ϕ ˙ ϕ ˙ ) = 0

ϕ ¨ + bata ( ρ ) ρ ˙ ϕ ˙ = 0

t ¨ + bronceado ( ρ ) ρ ˙ t ˙ = 0

Me doy cuenta de que están acoplados, pero honestamente no tengo ni idea de cómo voy a encontrar ρ ( τ ) , t ( τ ) y ϕ ( τ ) .

Primero, creo que es posible que falten algunos factores de 2 en el ϕ y t ecuaciones
En segundo lugar, un caso que se puede resolver fácilmente es cuando ρ está arreglado.
Hola Secavara, ¿podrías dar más detalles?
Oh, hice algunos garabatos para que pudiera estar equivocado, pero para ϕ y t Obtuve ϕ ¨ + 2 bata ( ρ ) ρ ˙ ϕ ˙ = 0 y t ¨ + 2 bronceado ( ρ ) ρ ˙ t ˙ = 0 . Y luego observe que si tomamos ρ ˙ = 0 , entonces podemos encontrar algunas geodésicas fácilmente.
En ese caso, t = a + b τ y ϕ = C ± b τ .

Respuestas (1)

Suponiendo que sus ecuaciones sean correctas (no las he verificado), pero las ODE acopladas como esta (y las ecuaciones geodésicas en general) deben resolverse numéricamente. El truco es emplear los solucionadores estándar de Runge-Kutta, es más conveniente tratar con sistemas de primer orden. En general, dadas las ecuaciones geodésicas:

X ¨ a + Γ b C a X ˙ b X ˙ C = 0 ,

podemos escribirlos como un sistema de primer orden:

X ˙ a = v a , v ˙ a = Γ b C a v b v C

Entonces, para su sistema anterior, defina X a = ρ ˙ a , y a = ϕ ˙ a , y z a = t ˙ a . Luego, junto con las ecuaciones restantes, tiene un sistema 6-D, de primer orden, y las soluciones dependerán de especificar 6 de tales condiciones iniciales.

Espero que esto ayude.

Difícilmente es un "truco" emplear solucionadores RK, ya que se encuentran entre los más populares y, a menudo, son los "a los que se recurre" antes de recurrir a métodos más elaborados.
@JamalS, creo que te tomaste mi oración demasiado literalmente. Los solucionadores de RK no funcionan de forma nativa en ODE de segundo orden, el "truco" es escribir ODE de segundo orden como ODE de primer orden y luego aplicar los métodos RK.
@Dr.IkjyotSinghKohli He logrado un progreso sustancial en mi problema, sin embargo, tengo problemas para especificar las condiciones iniciales que conducen a una solución estable de acuerdo con el solucionador que estoy usando. Gracias
@Cyphox32 hola. Sí. Esto es frecuentemente un problema en tales aplicaciones numéricas. Es posible que desee ver el uso de MATLAB/Octave/NumPY/SciPy en Python para este tipo de cosas. MATLAB tiene solucionadores numéricos muy potentes: ODe23s OdE45, etc. y es más fácil experimentar con condiciones iniciales.