Habiendo leído de un trabajo realizado para intentar diferenciar las velocidades de rotación de estrellas en el plano galáctico, tanto de la forma kepleriana habitual desarrollando el lagrangiano de una distribución de masa de estrellas, vi que se intentaba dar la explicación del aplanamiento de la curva de rotación, de un modo erróneo... de allí me surgió la idea de desarrollar el tema y ver si con los supuestos newtonianos y keplerianos, de distribuciones de masa bariónica y oscura se puede lograr de la función que grafique la famosa curva aplanada que es la evidencia más notable de la existencia de materia oscura, sea lo que sea la materia oscura.

Partimos del Lagrangiano de un objeto de masa m, en un campo gravitatorio de una galaxia que en vez de dos brazos suponemos un disco de densidad y espesor constante con el radio y además con las siguientes características

radio del bulbo galáctico

radio total de la galaxia

espesor del disco galáctico (de sección constante en función del radio) entre

masa total de la galaxia

masa bariónica de la galaxia (vía láctea )

masa bariónica del bulbo=10%

masa bariónica del disco=90%

masa de materia oscura de la galaxia

densidad de masa bariónica en el bulbo

densidad de masa bariónica en los brazos

= densidad de masa oscura en la galaxia (más adelante se verá el cálculo)


en general se puede definir el lagrangiano , para cada parte de la galaxia un vista general



pero habrá que ir analizando cómo cada parte de la galaxia bulbo, disco , halo, afectan a la rotación en función del radio la velocidad de rotación de la masa m respecto al centro galáctico se puede descomponen en velocidades radiales y tangenciales


de preferencia prefiero expresar las velocidades lineales con la letra v, las aceleraciones lineales con a y las respectivas angulares con y quedando


reemplazando 3 en 1


aplicando la ecuación de Euler-Lagrange con respecto a las variables del sistema


tenemos para el ángulo de rotación




y para la distancia radial






reemplazando 6,7,8y9 en 5 obtenemos las siguientes igualdades






entonces deducimos la conservación del momento angular





y que









reordenando y quitando la masa m que aparece en todo los términos


reemplazando la 10 en 11


sabiendo que si se multiplica la ecuación 12 por y se intenta dejar todos los términos expresados en derivadas temporales




que acomodando



resultando



por lo que


esa velocidad deducida es una constante del sistema que intentaremos deducir sabiendo que lo que hasta ahora llamamos y son masas asociadas al total de la galaxia cuando pero cuando queremos calcular velocidades en el interior de la galaxia lo que afecta gravitacionalmente al giro de la masa m es la masa de la porción de galaxia cuyo radio es menor a r y con suposiciones de simetría podemos determinar esa masa en función del radio.



La masa de materia oscura para una distribución esférica (aunque es sólo un aproximación) hasta el radio r la calculamos así: de las predicciones del modelo modelo de materia oscura fría (CDM por Cold Dark Matter), la distribución de la materia oscura en el halo y toda la galaxia sigue una distribución gaussiana del tipo




por las condiciones de contorno con , tenemos

de donde la masa para un determinado radio viene dada por


la densidad volumétrica queda




haciendo cambio de variable

Si expandimos la serie de la integral en como un polinomio de Taylor (fuente):










que indica que la masa de materia oscura crece muy proporcional al radio ya que los dos últimos términos son despreciables frente al primero dando



que para un radio aporta masa nula

usando la ecuación 13 pero en los radio menores al radio del bulbo galáctico podemos despejar la velocidad radial y reemplazarla en 3





se puede ver que tanto como tienden a 0 cuando r tiende a cero, por lo que si en el centro de la galaxia la velocidad radial es nula entonces la constante también lo es...

de modo entre 0 y a la curva se aproxima a una recta ascendente, ya que la densidad de materia bariónica crece con el cubo del radio y la de materia oscura crece solo proporcional al radio esta se hace despreciable en pequeños radios









Ahora haremos un desarrollo similar pero para el disco galáctico, trabajos mejor fundamentados que este toman la densidad de masa constante en función del radio , pero tiene en cuenta que la distribución de estrellas perpendicularmente al disco tiene una distribución gaussiana, pero aquí solo aportaría ruido a la matemática medianamente compleja sin errores de cálculo apreciables solo tomaremos que es constante en el espesor de un anillo entre y .

replanteando la ecuación 17 para la siguiente porción de curva





usando definiciones y que




desarrollando aun más





que es una ecuación de la forma


comparada con la versión kepleriana donde no se tiene en cuenta la materia oscura


vemos que el término aplana la caída de la curva de rotación elevando la velocidad por encima de la versión Kepleriana , tal cual demuestran las observaciones.

un análisis cuantitativo en función de la masa total de la galaxia y







donde corrige las unidades ,
La contribución de la materia oscura a la velocidad de rotación es un orden mayor a la contribución de la materia bariónica, por lo que las curvas resultan casi planas , a partir de radio similares al del bulbo galáctico.

Como el módulo de la velocidad de rotación es casi constante, la velocidad angular del centro de la galaxia debe ser mayor que en la periferia, creando ante una inhomogeneidad en la distribución de estrella, la característica espiral...

Espero haber investigado correctamente, no haber copiado demasiados errores ajenos y contribuir a esclarecer de donde sale el aplanamiento de las curvas de rotación de las galaxias espirales.

Otra conclusión es la contribución de la materia oscura a la rotación del bulbo galáctico, por los valores que he alcanzado, no debería despreciarse, alcanzando mayores velocidades de rotación que del modo Kepleriano.
Haz clic en la imagen para ampliar  Nombre:	vrot.png Vitas:	1 Tamaño:	11,3 KB ID:	340873








Cualquier contribución o corrección es bienvenida.

Fuentes entre otras

https://francis.naukas.com/2013/06/1...os-galacticos/

https://arxiv.org/pdf/1302.2728.pdf

https://es.wikipedia.org/wiki/Ecuaci...Euler-Lagrange

https://es.wikipedia.org/wiki/Funci%C3%B3n_error

http://www.astro.ugto.mx/~rcoziol/Cu...mponentsMW.pdf

http://www.das.uchile.cl/~jose/as200...se15_2017A.pdf