¿Existe una expresión analítica aproximada para la distribución de masa radial de la Vía Láctea?

Encontré la imagen a continuación en el artículo de Space.com Este mapa en color en 3D de 1.700 millones de estrellas en la Vía Láctea es el mejor jamás creado , aunque no es el mapa mencionado en el título .

El título de esta imagen dice:

Esta imagen de velocidad radial muestra el movimiento de 7 mil millones de estrellas. Los colores van del azul (estrellas que se mueven a 50 km/s hacia nosotros) al rojo (estrellas que se mueven a 50 km/s de nosotros). El color blanco muestra cuando, en promedio, las estrellas no se mueven en la línea de visión con respecto a nosotros. Las estrellas que se quedan atrás mientras orbitan el centro de la Vía Láctea parecen estar alejándose de nosotros, y las que están acelerando parecen estar viajando hacia nosotros. Crédito: ESA/Gaia/DPAC

Si imaginas una banda a lo largo del ecuador galáctico, la velocidad dominante muestra dos "picos" positivos y dos negativos, con un cruce por cero en la dirección del centro galáctico.

Solo por diversión, quería ver si podía reproducir este comportamiento con un cálculo simple basado en un cálculo 2D asumiendo un movimiento circular y una distribución de densidad radial. ρ ( r ) que luego podría usar para calcular una distribución de velocidad de rotación v ( r ) , pero rápidamente me di cuenta de que no tengo idea de cómo sería el perfil de densidad.

  1. A los efectos de este sencillo ejercicio, ¿cuál sería una expresión analítica que coincidiera aproximadamente con el perfil de densidad radial de la Vía Láctea, proyectado sobre su plano ecuatorial?

  2. Para distribuciones esféricamente simétricas, el teorema de la capa de Newton permite tratar toda la masa dentro de una esfera definida por el radio de una órbita como si estuviera en el centro e ignorar toda la masa en la capa fuera de ese radio. ¿Hay algo parecido a esto para una distribución radial dentro de un plano?

velocidad aparente de la Vía Láctea vista desde la Tierra (Gaia)

Funciona al revés. Obtienes la distribución de densidad mirando los datos de velocidad.
@RobJeffries No estoy tratando de ser un astrónomo, sino que solo intento hacer un ejercicio matemático para obtener una mejor comprensión, y así puedo trabajar en la dirección que quiera. Quiero trabajar un problema simple en la dirección indicada.
Es posible que desee consultar esta pregunta que no aborda específicamente su pregunta aquí, pero puede brindarle un buen punto de partida.
@zephyr es un artículo útil y también bastante interesante para leer, ¡gracias por la sugerencia!

Respuestas (2)

A los efectos de este sencillo ejercicio, ¿cuál sería una expresión analítica que coincidiera aproximadamente con el perfil de densidad radial de la Vía Láctea, proyectado sobre su plano ecuatorial?

El ejemplo más simple del trabajo de los astrónomos es el perfil de densidad radial de una esfera isotérmica única (SIS). Se llama así porque es esféricamente simétrico (y, por lo tanto, aplicable a un plano 2D para sus propósitos) y todos los objetos orbitan con la misma velocidad (y por lo tanto tienen la misma "temperatura", por lo tanto, isotérmicos). El perfil de densidad toma la forma:

ρ ( r ) = v 2 4 π GRAMO r 2

dónde v es la velocidad de rotación. Tenga en cuenta que puede ver otras formulaciones que usan σ v en vez de v . En este caso, están utilizando la dispersión de velocidad , que es ligeramente diferente a la velocidad de rotación.

Se han encontrado otros perfiles de densidad más realistas ejecutando simulaciones del Universo y haciendo coincidir ecuaciones funcionales con los perfiles de densidad de las galaxias resultantes. Tales resultados populares son el perfil NFW y el perfil Einasto .

El perfil NFW es una función de dos parámetros, dada por

ρ ( r ) = ρ 0 r R S ( 1 + r R S ) 2

dónde ρ 0 y R S y dos, parámetros dependientes del halo.

El perfil de Einasto es nuevamente un modelo de dos parámetros dado por

ρ ( r ) Exp ( A r α )

dónde A y α son parámetros configurables.

Para distribuciones esféricamente simétricas, el teorema de la capa de Newton permite tratar toda la masa dentro de una esfera definida por el radio de una órbita como si estuviera en el centro e ignorar toda la masa en la capa fuera de ese radio. ¿Hay algo parecido a esto para una distribución radial dentro de un plano?

El teorema de Shell para la gravedad no se extiende a un anillo 2D. Sin embargo, diré que cuando se habla de las órbitas de las estrellas en las galaxias, la masa de las estrellas fuera de la órbita de una estrella generalmente se considera insignificante. La razón principal de esto es que es la Materia Oscura la que comprende la mayor parte de la masa de una galaxia y la que más contribuye a definir la órbita de una estrella en una galaxia. A menudo se supone que el halo de materia oscura es esféricamente simétrico, en cuyo caso se aplica el teorema de la capa de Newton y la masa que le preocupa al determinar la órbita de una estrella es la masa del interior del halo de materia oscura en la órbita de la estrella.

Gracias, tengo curiosidad sobre el Teorema de Shell; ¿realmente se aplica a una distribución simétrica cilíndrica 2D en nuestro universo 3D en oposición a un universo 2D (teórico)? Todavía no he encontrado una fuente para esto.
-1 Prueba pendiente de que el teorema de la capa se aplica a un disco (o incluso a distribuciones elipsoidales no esféricas). Creo que el "teorema del anillo" de Newton se aplicaría en 2D si la gravedad disminuyera como r 1 .
@RobJeffries Ahora que lo miro, creo que tienes razón. Editaré mi respuesta.
Ok, pero ¿la mayor parte de la masa, digamos el interior de la órbita del Sol, es materia oscura? De lo contrario, entonces tiene una distribución de masa no esférica complicada y el teorema de la capa no se aplica. ¿No es por eso que las órbitas estelares no están cerradas?

@Rob Jeffries mencionó que "Obtienes la distribución de densidad al observar los datos de velocidad". También creo que esto es lo que está buscando, así que le daré algunos detalles de cálculo.

Suponiendo simetría esférica y movimiento circular, la gravedad es igual al movimiento circular como

α GRAMO METRO metro R 2 = metro v 2 R
dónde GRAMO es constante gravitacional, METRO es la masa encerrada a la distancia radial R , v es la velocidad tangencial (no radial), metro es una masa de prueba, y α es una constante para el potencial efectivo, que depende de la forma supuesta del potencial. Entonces, podemos expresar la masa METRO como el perfil de densidad ρ . Con la simetría esférica, el perfil ρ METRO R 3 . Por lo tanto,
α GRAMO ρ R 2 = v 2 .

Dado que por observación podemos construir la curva de rotación que es v = F ( R ) , el perfil de densidad es entonces una función que depende únicamente de R : ρ = gramo ( R ) , es decir, la distribución de masa radial.

Algunas notas incluyen i) la masa METRO incluye materia oscura; ii) v es la velocidad tangencial, no la velocidad radial como se presenta en la figura que mencionaste.

¡Rob Jeffries es genial! Pero esta es una respuesta a una pregunta que no hice, y definitivamente no es una respuesta a la pregunta formulada. Vea mi comentario y verifique la pregunta nuevamente.
"A los efectos de este sencillo ejercicio, ¿cuál sería una expresión analítica que coincidiera aproximadamente con el perfil de densidad radial de la Vía Láctea, proyectado en su plano ecuatorial?" R: "Más o menos", ρ = gramo ( R ) dado v = C o norte s t a norte t
Dado que los comentarios se consideran temporales en Stack Exchange, sería mejor si aborda la pregunta directamente en la respuesta, en lugar de comentarios adicionales. Además, como no tengo gramo ( R ) útil, lo que estoy pidiendo es una expresión analítica que coincida aproximadamente con la distribución de la Vía Láctea.
Sí. Tienes razón. Realmente me olvidé de eso. Estaba pensando en objetos distantes.
Y gramo ( R ) tiene una forma funcional como se define allí arriba. Simplemente no lo escribí.
¿Cómo es aproximadamente la distribución de masa radial de la Vía Láctea? ¿Es plano, trapezoidal, triangular, gaussiano, parabólico, tipo erf? ¡No tengo ni idea!
Dado v = C o norte s t a norte t , ρ R 2 .
¿Qué tan cerca está eso de la situación de la Vía Láctea? Estoy buscando una afirmación que describa nuestra galaxia. Si esto es lo que parece, incluya una declaración a tal efecto en su respuesta, porque esto es lo que estoy preguntando, y presentar esto como un "dado" en lugar de decirlo por adelantado significa que podría ser cierto o no. .
La curva de rotación es aproximadamente constante. Para acercarse más o menos a un determinado perfil de galaxia, debe ajustar esos parámetros. No conozco los parámetros exactos para MW. Así que dejaré mi respuesta así como un concepto general.
Buscando en Google, acabo de encontrar Modelado de distribución masiva de la galaxia de la Vía Láctea usando el mapa de mil millones de estrellas de Gaia, por ejemplo. Citar esto y posiblemente incluir una o ambas partes de la figura 1 ( i.stack.imgur.com/oKYNy.png y i.stack.imgur.com/d7dl9.png ) o algo similar y agregar las palabras "aproximadamente plano" sería hazlo. Sólo un pensamiento.