Esta es una pregunta interesante y, a menudo, las preguntas interesantes no pueden responderse fácilmente con ningún conocimiento actual, pero esta puede responderse hasta cierto punto. Voy a repasar los conceptos básicos de la teoría orbital y describir cómo pueden aplicarse a las galaxias y cómo difiere de los sistemas keplerianos. Debe tener una comprensión razonable de la física newtoniana (después de todo, las órbitas derivan precisamente de las leyes de Newton) y un fuerte conocimiento de las matemáticas. Si no tiene estas cosas, salte al final de cada sección donde trataré de resumir los puntos importantes detrás de las matemáticas.
Una nota rápida sobre notación matemática que usaré. Un punto sobre un símbolo indica una derivada del tiempo (por ejemplo, ) y los símbolos en negrita sin cursiva son cantidades vectoriales (por ejemplo, ). Vamos a ir al grano.a˙F
La ecuación orbital del movimiento
Considere una masa como una posición y moviéndose con un movimiento descrito por . Esta masa experimenta una fuerza que es solo una función de la distancia radial, , desde el centro del sistema de coordenadas. El objetivo aquí es determinar la ecuación de movimiento que puede describir la órbita de la masa debido a esta fuerza. Esta ecuación se puede usar para resolver . Según la ley de Newton, la ecuación de movimiento puede definirse inicialmente comomrr˙F(r)rr(θ)
F(r)=ma=m(r¨−rθ˙2)
Tenga en cuenta que en este caso es simplemente el componente radial de y es el ángulo acimutal del cuerpo en un sistema de coordenadas esféricas. Te lo dejo a ti para determinar cómo dividir la aceleración en los dos componentes anteriores, bajo el sistema de coordenadas apropiado. Intentemos eliminar nuestra dependencia para que solo tengamos una función de . Esto se puede lograr utilizando la conservación del momento angular. El momento angular por unidad de masa viene dado por para que . Esto darrθθrℓ=r2θ˙θ˙=ℓ/r2
F(r)=m(r¨−ℓ2/r3)
Esta es ahora una ecuación diferencial que nos permite resolver , pero queremos por lo que debemos hacer alguna conversión. Vuelva a parametrizar definiendo (la razón se aclarará en un momento) y determinando en términos de y .r(t)r(θ)u≡1/rr¨uθ
ddt(r)=ddt(1u)=1u2dudt=1u2dudθdθdt=−θ˙u2dudθ=−ℓdudθ
Tenga en cuenta la sustitución de . Ahora diferencie nuevamente para determinar .ℓ=r2θ˙=θ˙/u2r¨
d2dt2(r)=−ℓddt(dudθ)=−ℓdθdtddθ(dudθ)=ℓθ˙d2udθ2=−ℓ2u2d2udθ2
Poniendo esto en nuestra expresión para la ecuación de movimiento y haciendo la transformación que finalmente dar=1/u
F(1/u)=m(−ℓ2u2d2udθ2−ℓ2u3)
Escribiendo en una forma más conveniente finalmente llegamos a
d2udθ2+u=−F(1/u)mℓ2u2
Recuerda que es la masa del cuerpo, , es el momento angular por unidad de masa, es una fuerza puramente radial que actúa sobre el cuerpo, y y son las posiciones de coordenadas radiales y azimutales de la masa.mu(θ)≡1/r(θ)ℓFrθ
Punchline : El resultado final aquí es una ecuación general de movimiento para un cuerpo en órbita de acuerdo con alguna fuerza arbitraria. Esto podría ser gravedad, electromagnética, una fuerza de resorte o cualquier otra cosa que decidamos. Se deriva a propósito bajo suposiciones generales y no restrictivas y, con suerte, se puede ver que se puede usar para comprender el movimiento orbital de una estrella que orbita en una galaxia de disco. El objetivo con esta ecuación debe ser conectar su fuerza (lo que sea que sea) y resolver para . A partir de ahí, es fácil determinar .u(θ)r(θ)
Movimiento Kepleriano
Antes de comenzar a mirar el movimiento orbital en una galaxia, veamos el movimiento kepleriano estándar para que tengamos algo con lo que compararnos. El movimiento kepleriano se deriva de suponer que nuestra masa está orbitando una sola masa puntual y bajo la influencia de la gravedad simple. En ese caso, podemos escribir nuestra fuerza como y, por lo tanto, , donde es una constante, definida aquí por simplicidad matemática. Tenga en cuenta que es la constante gravitacional. La ecuación orbital general, bajo esta fuerza, ahora se convierte enmMF(r)=kr−2F(1/u)=ku2k≡GMmG
d2udθ2+u=−kmℓ2
Esta es una ecuación diferencial estándar no homogénea de segundo orden con una función de fuerza constante. Si conoce su ecualizador diferencial, debe conocer la solución casi de inmediato.
u(θ)=kmℓ2+Acos(θ−θ0)
En esta ecuación, es una constante desconocida, y representa el ángulo inicial de la órbita, que arbitrariamente podemos elegir que sea cero. Nuestro objetivo final es obtener así que hagamos eso. Sin embargo, voy a hacer algunos pasos en uno y dejar que trabajes en las matemáticas que intervienen. para y el momento angular como donde es la masa reducida de nuestro sistema . También declararé, sin pruebas, que donde es la excentricidad de la órbita.Aθ0r(θ)k=GMmL=ℓμμe=A(mℓ2/k)e
r(θ)=L2/GMμ21+ecos(θ)
Remate : hemos encontrado una ecuación final que representa el movimiento orbital de una masa bajo la influencia de una gravedad debido a un punto similar a la masa . Si conoce sus cosas, verá que esta ecuación describe con precisión las secciones cónicas, dependiendo del valor de . Si , obtienes movimiento circular (ya que convierte en una constante). Si , obtienes movimiento elíptico, es movimiento parabólico y es hiperbólico.Mee=0r(θ)0<e<1e=1e>1
Como era de esperar, el movimiento kepleriano (es decir, tener una fuerza central tal que ) dio como resultado cónicas, que es precisamente la primera ley de Kepler. Las leyes segunda y tercera de Kepler se derivan más o menos de los mismos supuestos. Es lógico entonces, que cualquier sistema donde no cumpla con ninguna de las leyes de Kepler. Las órbitas no son cónicas perfectas (por ejemplo, elipses, círculos, etc.), no barren áreas iguales en tiempos iguales, y el estándar ciertamente no se aplica.F∝r−2F∝̸r−2P2∝a3
Movimiento orbital en una galaxia
Su pregunta describe correctamente la situación de las estrellas (o cualquier cosa realmente) que orbitan en una galaxia. Las estrellas no están orbitando masas centrales, puntiagudas. Están incrustados dentro de la materia bariónica y oscura que comprende la galaxia y están orbitando a través de ella. Es un concepto bien conocido en física que las distribuciones de masa esféricamente simétricas no tienen una atracción gravitacional neta en los objetos interiores a esa distribución, lo que significa que para las estrellas en una galaxia, la masa que afecta a su órbita es la masa interior a su radio. Si ese radio cambia, ¡la masa cambia!
La fuerza central sobre nuestra estrella seguirá siendo la gravedad, pero la masa que actúe sobre ella será todo el interior de la masa hasta cierto radio, denotado por . Podemos ver que . Si queremos determinar la fuerza que actúa sobre nuestra estrella (y, por lo tanto, la órbita exacta, a través de la ecuación diferencial anterior), primero tenemos que averiguar cuál es la masa interior a cierto radio. Esto se puede lograr utilizando la ecuación de continuidad de masa .MrF(r)=GMr(r)m/r2
dMrdr=4πr2ρ(r)
Básicamente, puede calcular todo el interior de la masa en integrando toda la densidad de masa en función de . Aquí, por supuesto, necesita una buena ecuación para . Un perfil de densidad simple pero físicamente poco realista es la Esfera Isotérmica Singular (SIS), mientras que una ecuación más realista pero matemáticamente compleja podría ser el perfil NFW o el perfil de Einasto .rrρ(r)
Ahora he presentado todos los pasos que necesitas para descubrir el movimiento orbital en una galaxia, pero tengo que decir que no es bonito. Sin embargo, podemos ver parte del caso más simple, el del SIS.
Esfera isotérmica única
Para la esfera isotérmica individual, tienes que , donde es la velocidad de rotación de tu estrella. Este perfil se basa en un hecho crucial de las galaxias de disco su perfil de rotación es plano! Esto estaba bien establecido, por ejemplo, por Ruben et al. 1978 . He reproducido una figura de este documento a continuación que muestra la curva de rotación para varias galaxias. ¡El punto importante aquí es que esto muestra que es constante y no depende del radio! (Asumiendo que no estamos cerca del bulto o centro galáctico. Esa es una bestia completamente diferente).ρ(r)=v2/(4πGr2)v−v
Con esta información crucial, podemos resolver integrando (que te dejo a ti). El resultado es queMrρ(r)
Mr=v2rG
Esto significa que tu fuerza está dada por
F(r)=v2rmr2=v2mr⇒F(1/u)=v2mu∝ku
Puede ver aquí que, a diferencia del caso kepleriano, nuestra fuerza es proporcional a lugar de . Puede hacer este proceso con otros perfiles de densidad (como el NFW o Einasto I enumerados anteriormente), pero terminará con el mismo resultado.r−1r−2
Si está tan inclinado, podría optar por conectar esto a la ecuación de movimiento orbital anterior y resolverlo, pero ahora está trabajando con una ecuación diferencial no lineal y las cosas pueden volverse desordenadas rápidamente.
Punchline : No estoy seguro de si esto realmente responde a su pregunta o no. Te conduje parcialmente por la madriguera del conejo, pero espero que puedas apreciar lo complejo que se vuelve rápidamente. Todo el trabajo anterior utilizaba suposiciones y simplificaciones generales. Supongo que la respuesta breve a todo esto es que las estrellas orbitan galaxias en una órbita compleja pero cerrada que no se describe fácilmente con precisión (incluso para nuestra propia galaxia) a través de ecuaciones calculables. Podemos aproximarnos y hacer nuestro mejor esfuerzo para trabajar con las matemáticas, pero al final es una aproximación. Sin embargo, en las aproximaciones más crudas, también puede considerar que una órbita como nuestra estrella sea circular y termine con ella.