Simulación de colisión de dos galaxias para estudiar la formación de las galaxias enanas

esferoidales satélites de la Vía Láctea

Omar Alfonso Bohórquez Pacheco

Universidad Nacional de Colombia Facultad de Ciencias, Departamento de Física Bogotá, Colombia 2016

Simulación de colisión de dos galaxias para estudiar la formación de las galaxias enanas esferoidales satélites de la Vía Láctea

Omar Alfonso Bohórquez Pacheco

Trabajo de grado presentado como requisito parcial para optar al título de: Magister en Ciencias - Física

Director: Dr. rer. nat. Rigoberto Casas Miranda

Grupo de Investigación: Astrofísica

Universidad Nacional de Colombia Facultad de Ciencias, Departamento de Física Bogotá, Colombia 2016

A mis padres y en especial a mi abuela Esther Vega viuda de Pacheco, sin su ayuda no hubiera sido posible.

Posible es que me equivoque y tome por oro y diamantes lo que solo es cobre y vidrio.

René Descartes

Agradecimientos

Un gran agradecimiento al CECAD - Centro de Computación de Alto Desempeño de la Universidad Distrital “Francisco José de Caldas” y a mis amigos José Benavides, Diana Judith Cubillos Jara y Yeimy Camargo, sin su ayuda no habría sido posible.

Resumen y Abstract IX

Resumen

En la actualidad dentro del área de la astrofísica se presentan un sinnúmero de problemas sin resolver, entre ellos el problema del origen de las galaxias satélite de la VL. Estas galaxias se caracterizan por ser de tipo enana esferoidal. La mayoría de estas se encuentran distribuidas en una estructura tipo disco que se encuentra dispuesta de casi forma perpendicular al plano de la galaxia, esta estructura es conocida con el nombre de disco de satélites (DoS) o Vast Polar Structure Of Satellite Galaxies (VPOS).

Hasta el momento no se ha podido encontrar un modelo que dé cuenta de la cantidad y de la distribución espacial de estas galaxias. Sin embargo se han presentado varias propuestas de solución para el caso de la VL, una de las cuales propone que estas tuvieron origen en la colisión de dos galaxias de disco hace miles de millones de años. En este trabajo se llevaron a cabo simulaciones numéricas de N-cuerpos con el software Gadget2 para colisiones entre dos galaxias de disco que pudieron dar origen al disco de satélites de la VL.

Se puede decir que bajo las condiciones iniciales planteadas para estas simulaciones no es posible generar las características dinámicas y espaciales de lo que se denomina el disco de galaxias satélite la VL (DoS), sin embargo, esto no descarta por completo el modelo propuesto para este estudio, ya que si bien no refleja los resultados esperados, se pueden proponer trabajos futuros para realizar simulaciones y modificar los parámetros de impacto y las relaciones de masa y así poder enmarcar las condiciones físicas bajo las cuales sería posible este evento. Los resultados de este trabajo permiten acotar aún más el espacio de condiciones iniciales orbitales y de masa de eventuales progenitores del DoS de la VL.

Palabras clave: Grupo Local, Galaxias Enanas, Vía Láctea

X Colisión de dos Galaxias de Disco

Abstract

At present within the area of astrophysics there are a number of unresolved problems, including the problem of the origin of the satellite galaxies of the . These galaxies are characterized by dwarf spheroidal galaxies type. Most of them are distributed in a disk-like structure which is arranged almost perpendicular to the plane of the galaxy, this structure is known as satellite disk (DoS) or Polar Structure Of Vast Satellite Galaxies (VPOS).

So far there is not a model that fully reproduces the amount and spatial distribution of these galaxies. However there have been several proposed for the solutions, one of which suggests that these originated in the collision of two disk galaxies billions of years ago. In this work we have performed N-bodies numerical simulations of with Gadget2 of the collision between two disk galaxies that could give rise to disc of VL satellites.

It can be said that under the initial conditions set for these simulations is not possible to generate dynamic and spatial characteristics of what is called the satellite the disc of Milky Way (DoS) galaxies, however, this does not preclude in any way the model proposed for this study because although it does not reflect the expected results can be proposed for future work simulations and modify the parameters of impact and the mass ratios to frame the physical conditions under which this event possible. The results of this work allow limiting further orbital space initial conditions and eventual mass progenitors of DoS VL.

Keywords: , Dwarf Galaxies, Milky Way Contenido XI

Contenido

Pág.

Resumen ...... IX

Lista de figuras ...... XIV

Lista de tablas ...... XXIV

Introducción ...... 1

1. Galaxias, Grupo Local y Galaxias Satélites ...... 5 1.1 Galaxias ...... 5 1.1.1 Secuencia de Hubble ...... 6 1.1.2 Galaxias elípticas ...... 6 1.1.3 Galaxias espirales ...... 8 1.1.4 Galaxias irregulares ...... 9 1.1.5 Galaxias enanas ...... 9 1.1.6 Vía Láctea (VL) ...... 11 1.2 Grupo Local ...... ¡Error! Marcador no definido. 1.2.1 Distribución espacial del Grupo Local ...... 16 1.2.2 Propiedades del Grupo Local ...... 19 1.2.3 Galaxias enanas del Grupo Local ...... 22 1.3 Galaxias enanas satélite de la Vía Láctea ...... 25 1.3.1 Disco de satélites de la Vía Láctea (DoS) ...... 28 1.3.2 Hipótesis de formación del disco de galaxias (DoS) ...... 30

2. Marco teórico ...... 35 2.1 Componentes de una galaxia de disco ...... 35 2.1.1 Bulbo central ...... 36 2.1.2 Disco ...... 36 2.1.3 Halo estelar ...... 36 2.1.4 Halo de materia oscura ...... 36 2.2 Modelos teóricos para el potencial gravitacional de las estructuras (Bulbo, Disco y Halo de materia oscura) ...... 37 2.2.1 Potencial de Hernquist (Bulbo central) ...... 37 XII Colisión de dos Galaxias de Disco

2.2.2 Perfil de densidad (Disco) ...... 38 2.2.3 Potencial de Hernquist (Halo de materia oscura) ...... 39 2.3 Orbita pre-colisión ...... 40 2.3.1 Parámetro de impacto ...... 41 2.4 Colisión de galaxias ...... 42 2.4.1 Sistemas no colisiónales ...... 42 2.4.2 Sistemas estelares interactuantes (definiciones) ...... 43 2.4.3 Encuentro de galaxias en rotación ...... 44 2.4.4 Encuentros de discos planos (estabilidad) ...... 46 2.4.5 Fricción dinámica ...... 48 2.4.6 Algunas colisiones reales y simulación ...... 48

3. Condiciones iniciales ...... 51 3.1 Software Zeno ...... 51 3.2 Parámetros de la Vía Láctea...... 52 3.2.1 Condiciones iniciales baja resolución Vía Láctea ...... 53 3.2.2 Condiciones iniciales alta resolución Vía Láctea ...... 54 3.2.3 Masa del halo de materia oscura ...... 55 3.3 Galaxias satélites ...... 56 3.3.1 Redimensionamiento del Halo de materia oscura ...... 56 3.3.2 Redimensionamiento del disco ...... 56 3.3.3 Redimensionamiento del bulbo ...... 57 3.3.4 Corroborando el modelo ...... 58 3.3.5 Parámetros de las galaxias satélites ...... 59 3.3.5.1 Condiciones iniciales de las galaxias satélites de baja resolución ...... 61 3.3.5.2 Condiciones iniciales de las galaxias satélites de alta resolución ...... 65 3.4 Parámetros orbitales ...... 69 3.4.1 Simulaciones de baja resolución...... 69 3.4.1.1 Condiciones orbitales iniciales para simulaciones de baja resolución ..... 70 3.4.1.2 Condiciones orbitales iniciales de simulaciones de alta resolución ...... 73

4. Simulaciones y Análisis de Resultados ...... 77 4.1 Gadget2 ...... 77 4.2 Resultados simulaciones de baja resolución ...... 77 4.2.1 Simulaciones de baja resolución galaxia satélite 4% ...... 78 4.2.2 Simulaciones de baja resolución galaxia satélite 8% ...... 80 4.2.3 Simulaciones de baja resolución galaxia satélite 12% ...... 82 4.2.4 Simulaciones de baja resolución galaxia satélite 16% ...... 84 4.3 Análisis simulaciones baja resolución ...... 86 4.4 Resultados simulaciones de alta resolución ...... 86 4.4.1 Simulaciones de alta resolución galaxia satélite 4% ...... 86 4.4.2 Simulaciones de alta resolución galaxia satélite 8% ...... 88 Contenido XIII

4.4.3 Simulaciones de alta resolución galaxia satélite 12% ...... 90 4.4.4 Simulaciones de alta resolución galaxia satélite 16% ...... 92 4.5 Análisis de las simulaciones de alta resolución ...... 94 4.5.1 Análisis de grupos ...... 95 4.5.2 Grupos galaxia satélite 16 % ...... 95 4.5.3 Perfil de densidad Vía Láctea 16% ...... 97 4.5.4 Dispersión de velocidades en z para la Vía Láctea 16% ...... 98 4.5.5 Plano de la Vía Láctea 16% ...... 98 4.5.6 Dispersión de velocidades grupos del 16% ...... 99 4.5.7 Plano grupos 16% ...... 101 4.5.8 Planos y ángulo entre planos 16% ...... 101 4.5.9 Grupos para la galaxia satélite 12% ...... 103 4.5.10 Perfil de densidad Vía Láctea 12% ...... 105 4.5.11 Dispersión de velocidades en z para la Vía Láctea 12% ...... 106 4.5.12 Plano de la Vía Láctea 12% ...... 106 4.5.13 Dispersión de velocidades en grupos del 12% ...... 107 4.5.14 Plano grupos 12% ...... 108 4.5.15 Planos y ángulo entre planos 12% ...... 108

5. Conclusiones ...... 111

A. Apéndice: Redimensionamiento de estructuras (Cálculos) ...... 115

B. Anexo: Resultados de simulación ...... 121

Bibliografía ...... 143

Contenido XIV

Lista de figuras

Pág. Figura 1.1: Secuencia de Hubble. Las galaxias elípticas (E) vienen acompañadas por un número que hace referencia a la excentricidad de la galaxia, en las galaxias espirales (S) las letras a, b y c hacen referencia a la separación de los brazos, la letra a designa los brazos más cerrados y la letra c hace referencia a los brazos más separados, las galaxias espirales barradas (SB) las letras a, b y c hacen referencia a la longitud de la barra, la letra a designa barras cortas y la letra c barras más largas. Tomado de http://es.wikipedia.org/wiki/Secuencia_de_Hubble. 1 de abril de 2015...... 6 Figura 1.2 Galaxia elíptica M87 Tomado de https://www.astro.virginia.edu/ ...... 7 Figura 1.3: Galaxia espiral M31 Tomado de www.bbc.co.uk ...... 8 Figura 1.4: Galaxia irregular la Gran nube de Magallanes Tomado de http://www.surastronomico.com/ ...... 9 Figura 1.5: Galaxia enana de Fornax tomado de www.abc.es...... 10 Figura 1.6: VL visión artística. Tomada de http://apod.nasa.gov/ ...... 12 Figura 1.7: Esquema representativo de la disposición espacial de algunos de los objetos del Grupo Local. Tomada de http://www.australianscience.com.au ...... 14 Figura 1.8: La figura muestra que la mayoría de los miembros del LG se concentran en dos subgrupos que se centran en la VL y M31. El círculo de radio de 1180 kpc corresponde a la superficie de la velocidad cero del LG. El círculo de trazos con el radio 450 kpc muestra el radio de la esfera que contiene la mitad de todas las galaxias del Grupo Local (van den Bergh, 1999) ...... 16 Figura 1.9: Distribución de los miembros del Grupo Local, en esta sugiere la existencia de un orden espacial de tipo anillo en el Grupo Local (LG) (Sawa y Fujimoto 2005)...... 17 Figura 1.10: Distribución espacial tridimensional de los miembros del Grupo Local (Sawa and Fujimoto, 2005) vista desde la dirección ( l,b ) = (296° ,−11°). Esta es una vista de canto del plano orbital de la galaxia VL y de la M31. Los círculos rellenos oscuros denotan los miembros del grupo de la VL, los grises son los miembros de M31 y los círculos sin relleno no pertenecen a ninguno de los dos grupos. Se puede ver como los Contenido XV miembros del Grupo Local están distribuidos en una estructura tipo disco con un espesor finito de 50 – 100 kpc. El tamaño de cada círculo representa su brillo cualitativamente. . 18 Figura 1.11: Esta es una toma vista de frente en el plano del LG, o del plano orbital de la galaxia M31. (Sawa y Fujimoto 2005)...... 19 Figura 1.12: Distribución Luminosidad para (LG). La figura muestra que las galaxias irregulares del LG tienen una función de luminosidad más baja que la de las galaxias dSph del Grupo Local. (Van den Bergh, 1999)...... 20

Figura 1.13: Relación entre luminosidad M V y la metalicidad [Fe/H] para las galaxias que pertenecen al Grupo Local (LG). Se muestra la línea de correlación que viene dada

por la expresión M V = 20 − 5[Fe/H] (van den Bergh, 1999)...... 20 Figura 1.14: Ilustración de la luminosidad versus la longitud de escala del disco para las galaxias del Grupo Local. Las galaxias irregulares (cuadrados rellenos) están por encima de la línea de tendencia, las espirales (signo +) y las enanas esferoidales (círculos). Tenga en cuenta que parece que no hay diferencias sistemáticas entre galaxias dSph en el subgrupo M 31 (círculos negros) y los del subgrupo (MW) (círculos abiertos). La línea de

tendencia en la figura es la relación M V = -15-7 logaritmo de longitud de escala (van den Bergh, 1999)...... 21 Figura 1.15 : Escala de masa frecuente para algunas galaxias enanas de la VL en función de la luminosidad (Strigari et al. 2008)...... 24 Figura 1.16 : Metalicidad estelar de estrellas población II versus la luminosidad bariónica para diferentes clases de galaxias enanas como se indica en la leyenda. Galaxias con deficiencia de gas (símbolos rellenos) y ricas en gas (símbolos sin relleno). Galaxias con la misma luminosidad, las antiguas poblaciones de dSphs son más ricas en metales que las de dIrrs. Así, en contraste con dIrrs, las dSphs deben haber experimentado un rápido enriquecimiento temprano. Tenga en cuenta la ubicación de los llamados dSph /dirr galaxias de tipo transición, que combinan propiedades de población de dSphs con la formación estelar en curso y un contenido de gas medible en el diagrama, estas propiedades hacen de las dSphs unas plausibles progenitoras (Grebel, Gallagher, & Harbeck, 2003)...... 24 Figura 1.17 : Ubicación espacial de algunas de las galaxias satélite de la VL que se encuentran a una distancia menor de 250 kpc. Tomada de http://www.mso.anu.edu.au/ ...... 25 Figura 1.18 : Distribución radial de los satélites de la VL, en la parte inferior se muestra la frecuencia a acumulada, se aprecia la acumulación de satélites en los primeros 100 Kpc (McConnachie and Irwin, 2006)...... 26 XVI Colisión de dos Galaxias de Disco

Figura 1.19 : Posición de los 11 satélites más cercanos a la VL, vistos desde un punto localizado en el infinito y en l = 167°,91. El disco de la VL se representa con la línea horizontal con −25 ≤ X/pc ≤ 25 y el centro de coordenadas es el centro galáctico. La línea punteada marca el plano ajustado visto de perfil para N=11 (Kroupa et al., 2005)...... 28 Figura 1.20 : Imagen ilustrativa de la posible forma del plano que se ajusta a la distribución de disco de las galaxias satélite de la VL, Tomada de http://www.newscientist.com/article/ ...... 29 Figura 1.21 : La distribución de todas las galaxias conocidas hoy en día en el Grupo Local que se define por la esfera velocidad cero. Las galaxias fuera de esta esfera de radio alrededor de 1,5 Mpc se alejan del Grupo Local, mientras que dentro de la esfera caen hacia nosotros. En esta se presentan todas las galaxias satélite dentro de unos 250 kpc alrededor de la VL, el disco galáctico de la VL se ve casi de canto la dirección del polo norte galáctico es hacia arriba. Los satélites de la VL con movimientos propios conocidos y velocidades radiales se muestran como círculos de colores. Los satélites rojos se están alejando del observador, los azules se están acercando al observador. (Kroupa et al, 2014)...... 30 Figura 1.22 : Isotropía materia oscura fría (Izquierda), Anisotropía Kroupa et al (Derecha)...... 31 Figura 2.1: Elementos constitutivos de una galaxia de disco. Tomada de http: //atenea.pntic.mec.es/Antares/modulo7/m7_u103_0.html ...... 35 Figura 2.2: Curva de rotación observada para la galaxia M33. Tomada de http://www.learner.org/ ...... 37 Figura 2.3: Orbita pre-colisión seguida por la “Galaxia fantasma” (esquina superior izquierda) que se aproxima a la VL...... 40 Figura 2.4: En esta gráfica se muestra un parámetro de impacto de aproximadamente 35 kpc...... 42 Figura 2.5: Imagen donde se muestra en colisiones reales de galaxias, se aprecia las colas de marea generadas por la interacción entre dos galaxias. Tomada de ttp://www.extragalactic.info/ ...... 43 Figura 2.6: Disco con spin anti-paralelo al momentum angular de la órbita (Toomre 1972)...... 45 Figura 2.7: Disco con spin paralelo al momentum angular de la órbita (Toomre 1972). 46 Figura 2.8: Este par de galaxias, NGC 4676, también conocido como "Los Ratones" por sus colas de estrellas y gas, han colisionado y eventualmente se fusionarán en una sola galaxia. Corrientes de material han sido tirados fuera de las galaxias por la fuerza de Contenido XVII gravedad, lo que provocó nuevo nacimiento de estrellas. Tomada de http://hubblesite.org...... 49 Figura 2.9: En esta figura se aprecia que las simulaciones numéricas de N-cuerpos pueden dar cuenta de las interacciones que se aparecían en las observaciones...... 49 Figura 3.1: VL simulada y virializada 2Gyr con 75000 partículas, verde (bulbo) y rojo (disco)...... 53 Figura 3.2: Energía de la VL; energía cinética (verde), energía potencial (rojo) y energía total (azul)...... 53 Figura 3.3: Vía Láctea simulada y virializada 2Gyr con 150000 partículas, verde (disco) y rojo (bulbo)...... 54 Figura 3.4: Energía de la VL; energía cinética (verde), energía potencial (rojo) y energía total (azul)...... 54 Figura 3.5: En esta figura se aprecian los diferentes valores de masa estimados para el halo de materia oscura de la VL, las técnicas utilizadas para estos estimativos son diversas, entre estas se encuentran; fracción de bariones cósmica, coincidencia de abundancia y la dinámica de satélites distantes entre otras. Se resaltan los valores 12 encontrados por los modelos de materia oscura fría (M=1,2 x 10 M⊙) y materia oscura 12 caliente (2.6 x 10 M⊙) (Frenk and White 2014)...... 55 Figura 3.6: Características galaxia satélite (4%) virializada 2Gyr con 75000 partículas, verde (disco) y rojo (bulbo)...... 61 Figura 3.7: Energía de la galaxia satélite (4%), energía cinética (verde), energía potencial (rojo) y energía total (azul)...... 61 Figura 3.8: Características galaxia satélite (8%) virializada 2Gyr con 75000 partículas, verde (disco) y rojo (bulbo)...... 62 Figura 3.9: Energía de la galaxia satélite (8%), energía cinética (verde), energía potencial (rojo) y energía total (azul)...... 62 Figura 3.10: Características galaxia satélite (12%) virializada 2Gyr con 75000 partículas, verde (disco) y rojo (bulbo)...... 63 Figura 3.11: Energía de la galaxia satélite (12%), energía cinética (verde), energía potencial (rojo) y energía total (azul)...... 63 Figura 3.12: Características galaxia satélite (16%) virializada 2Gyr con 75000 partículas, verde (disco) y rojo (bulbo)...... 64 Figura 3.13: Energía de la galaxia satélite (16%), energía cinética (verde), energía potencial (rojo) y energía total (azul)...... 64 Figura 3.14: Características galaxia satélite (4%) virializada 2Gyr con 350300 partículas, verde (disco) y rojo (bulbo)...... 65 XVIII Colisión de dos Galaxias de Disco

Figura 3.15: Energía de la galaxia satélite (4%), energía cinética (verde), energía potencial (rojo) y energía total (azul)...... 66 Figura 3.16: Características galaxia satélite (8%) virializada 2Gyr con 350300 partículas, verde (disco) y rojo (bulbo)...... 66 Figura 3.17: Energía de la galaxia satélite (8%), energía cinética (verde), energía potencial (rojo) y energía total (azul)...... 67 Figura 3.18: Características galaxia satélite (12%) virializada 2Gyr con 350300 partículas, (verde) disco y rojo (bulbo)...... 67 Figura 3.19: Energía de la galaxia satélite (12%), energía cinética (verde), energía potencial (rojo) y energía total (azul)...... 68 Figura 3.20: Características galaxia satélite (16%) virializada 2Gyr con 350300 partículas, verde (disco) y rojo (bulbo)...... 68 Figura 3.21: Energía de la galaxia satélite (16%), energía cinética (verde), energía potencial (rojo) y energía total (azul)...... 69 Figura 3.22: Condiciones orbitales iniciales galaxia satélite 4%, orbita parabólica y parámetro de impacto de 55 kpc...... 71 Figura 3.23: Condiciones orbitales iniciales galaxia satélite 8%, orbita parabólica y parámetro de impacto de 50 kpc...... 71 Figura 3.24: Condiciones orbitales iniciales galaxia satélite 12%, orbita parabólica y parámetro de impacto de 45 kpc...... 72 Figura 3.25: Condiciones orbitales iniciales galaxia satélite 16%, orbita parabólica y parámetro de impacto de 40 kpc...... 72 Figura 3.26: Condiciones orbitales iniciales galaxia satélite 16%, orbita parabólica y parámetro de impacto de 35 kpc...... 73 Figura 3.27: Condiciones orbitales iniciales para todas las galaxias satélite en las simulaciones de alta resolución, orbita parabólica y parámetro de impacto de 45 kpc...... 74 Figura 4.1: Posición final en el plano X-Z de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 4% y las azules (disco) y morado (bulbo) de la VL, Se debe tener presente que por la distribución espacial de las partículas, estas no formarían un grupo...... 78 Figura 4.2: Posición final en el plano X-Y de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 4% y las azules (disco) y moradas (bulbo) de la VL...... 79 Contenido XIX

Figura 4.3: Posición final en el plano Y-Z de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 4% y las azules (disco) y moradas (bulbo) de la VL...... 79 Figura 4.4: Posición final en el plano X-Z de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 8% y las azules (disco) y morado (bulbo) de la VL, Se debe tener presente que por la distribución espacial de las partículas, estas no formarían un grupo...... 80 Figura 4.5: Posición final en el plano X-Y de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 8% y las azules (disco) y moradas (bulbo) de la VL...... 81 Figura 4.6: Posición final en el plano Y-Z de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 8% y las azules (disco) y moradas (bulbo) de la VL...... 81 Figura 4.7: Vía Láctea de color azul disco y morado bulbo, partículas galaxia satélite rojo disco y verde bulbo. En esta imagen se aprecia una cola de marea que crece hasta los 500 kpc aproximadamente...... 82 Figura 4.8: Posición final en el plano X-Y de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 12% y las azules (disco) y moradas (bulbo) de la VL ...... 83 Figura 4.9: Posición final en el plano Y-Z de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 12% y las azules (disco) y moradas (bulbo) de la VL...... 83 Figura 4.10: Vía Láctea de color azul (disco) y morado (bulbo), partículas galaxia satélite rojo disco y verde bulbo. En esta imagen se aprecia una cola de marea con gran densidad de partículas que crece hasta los 200 kpc aproximadamente que comparada con la de 500 kpc de la simulación de 12% es pequeña...... 84 Figura 4.11: Posición final en el plano X-Y de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 12% y las azules (disco) y moradas (bulbo) de la VL...... 85 XX Colisión de dos Galaxias de Disco

Figura 4.12: Posición final en el plano Y-Z de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 12% y las azules (disco) y moradas (bulbo) de la VL...... 85 Figura 4.13: Posición final en el plano X-Z de los componentes del sistema después de una simulación de 8,8 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 4% y las azules (disco) y morado (bulbo) de la VL, Se debe tener presente que por la distribución espacial de las partículas, estas no formarían un grupo...... 87 Figura 4.14: Posición final en el plano X-Y de los componentes del sistema después de una simulación de 8,8 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 4% y las azules (disco) y moradas (bulbo) de la VL...... 87 Figura 4.15: Posición final en el plano Y-Z de los componentes del sistema después de una simulación de 8.8 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 4% y las azules (disco) y moradas (bulbo) de la VL...... 88 Figura 4.16: Posición final en el plano X-Z de los componentes del sistema después de una simulación de 9.5 Gyr, las partículas en rojo (Disco) y verde (Bulbo) son las correspondientes a la galaxia satélite de 8% y las azules (Disco) y morado (Bulbo) de la VL, Se debe tener presente que por la distribución espacial de las partículas, estas no formarían un grupo...... 89 Figura 4.17: Posición final en el plano X-Y de los componentes del sistema después de una simulación de 9,5 Gyr, las partículas en rojo (Disco) y verde (Bulbo) son las correspondientes a la galaxia satélite de 8% y las azules (Disco) y moradas (Bulbo) de la VL...... 89 Figura 4.18: Posición final en el plano Y-Z de los componentes del sistema después de una simulación de 9.5 Gyr, las partículas en rojo (Disco) y verde (Bulbo) son las correspondientes a la galaxia satélite de 8% y las azules (Disco) y moradas (Bulbo) de la VL...... 90 Figura 4.19: Posición final en el plano X-Z de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 12% y las azules (disco) y morado (bulbo) de la VL, por la distribución espacial de las partículas, estas podrían formar grupos (galaxias enanas satélite)...... 91 Contenido XXI

Figura 4.20: Posición final en el plano X-Y de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (Disco) y verde (Bulbo) son las correspondientes a la galaxia satélite de 12% y las azules (Disco) y moradas (Bulbo) de la VL...... 91 Figura 4.21: Posición final en el plano Y-Z de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 12% y las azules (disco) y moradas (bulbo) de la VL...... 92 Figura 4.22: Posición final en el plano X-Z de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 16% y las azules (disco) y morado (bulbo) de la VL, Se debe tener presente que por la distribución espacial de las partículas, estas podrían formar grupos (galaxias enanas satélite)...... 93 Figura 4.23: Posición final en el plano X-Y de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 16% y las azules (disco) y moradas (bulbo) de la VL...... 93 Figura 4.24: Posición final en el plano Y-Z de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 16% y las azules (disco) y moradas (bulbo) de la VL...... 94 Figura 4.25: Grupos generados después de la colisión en el plano X-Y para la galaxia satélite del 16% de masa de la VL...... 96 Figura 4.26: Grupos generados después de la colisión en el plano X-Z para la galaxia satélite del 16% de masa de la VL...... 96 Figura 4.27: Grupos generados después de la colisión en el plano Y-Z para la galaxia satélite del 16% de masa de la VL...... 97 Figura 4.28: Perfil de densidad del disco de la VL en azul antes de la colisión y rojo después, como se aprecia el perfil de densidad en ambas situaciones se comporta de manera similar en la región de interés entre 0 kpc y 16 kpc...... 97 Figura 4.29: Dispersión de velocidades en z (s) para la VL respecto a la coordenada radial del disco (R) en morado antes de la colisión y verde después, se observa un aumento en la dispersión de velocidades en z antes y después de la colisión...... 98 Figura 4.30: Plano de la VL obtenido con las posiciones de las partículas del disco de la Vía Láctea usando el método de mínimos cuadrados para dos variables...... 99 XXII Colisión de dos Galaxias de Disco

Figura 4.31: En esta gráfica se aprecia la dispersión de velocidades (S) de los 35 grupos encontrados en los remanentes de la galaxia satélite dejados después de la colisión, la media de dispersión para los grupos es de 234 km/s...... 100 Figura 4.32: En esta imagen se aprecian las dispersiones de velocidades de algunas galaxias enanas satélites de la VL respecto a R en (pc), se aprecian dispersiones de velocidad del orden de 10 km/s (Walker et al. 2008)...... 100 Figura 4.33: Plano de los grupos del 16% obtenido con las posiciones de los grupos encontrados usando el método de mínimos cuadrados para dos variables...... 101 Figura 4.34: Partículas rojas Vía Láctea, puntos azules grupos, plano verde Vía Láctea y plano morado grupos, longitudes en kpc...... 101 Figura 4.35: Partículas rojas Vía Láctea, puntos azules grupos, plano verde Vía Láctea y plano morado grupos, longitudes en kpc...... 102 Figura 4.36: Partículas rojas Vía Láctea, puntos azules grupos, plano verde Vía Láctea y plano morado grupos, se aprecia un ángulo entre los dos planos inferior a los 88°, longitudes en kpc...... 102 Figura 4.37: Grupos generados después de la colisión en el plano X-Y para la galaxia satélite del 12% de masa de la VL...... 104 Figura 4.38: Grupos generados después de la colisión en el plano X-Z para la galaxia satélite del 12% de masa de la VL...... 104 Figura 4.39: Grupos generados después de la colisión en el plano Y-Z para la galaxia satélite del 16% de masa de la VL...... 105 Figura 4.40: Perfil de densidad del disco de la VL en azul antes de la colisión y rojo después, como se aprecia el perfil de densidad en ambas situaciones se comporta de manera similar en la región de interés entre 0 kpc y 16 kpc...... 105 Figura 4.41: Dispersión de velocidades en z (s) para la VL respecto a la coordenada radial del disco (R) en morado antes de la colisión y verde después, se observa un aumento en la dispersión de velocidades en z antes y después de la colisión...... 106 Figura 4.42: Plano de la VL obtenido con las posiciones de las partículas del disco de la Vía Láctea usando el método de mínimos cuadrados para dos variables...... 106 Figura 4.43: En esta gráfica se aprecia la dispersión de velocidades media de los 32 grupos encontrados en los remanentes de la galaxia satélite dejados después de la colisión, el promedio de estas dispersiones es de 242 km/s...... 107 Figura 4.44: Plano de los grupos del 12% obtenido con las posiciones de los grupos encontrados usando el método de mínimos cuadrados para dos variables...... 108 Figura 4.45: Partículas rojas Vía Láctea, puntos azules grupos, plano verde Vía Láctea y plano morado grupos, longitudes en kpc...... 108 Contenido XXIII

Figura 4.46: Partículas rojas Vía Láctea, puntos azules grupos, plano verde Vía Láctea y plano morado grupos, longitudes en kpc...... 109 Figura 4.47: Partículas rojas Vía Láctea, puntos azules grupos, plano verde Vía Láctea y plano morado grupos. Se aprecia un ángulo entre los dos planos inferior a los 88°...... 109 Figura 4.48: Izquierda (Kroupa 2014), derecha simulación 12%, si bien se aprecia en la simulación una estructura tipo disco, se ven muy amontonadas las galaxias satélites sobre el plano de la VL, en las observaciones se aprecian distancias de hasta 200 kpc respecto al disco de la VL, En la simulación esta distancia llega ser máximo de 20 kpc respecto al disco...... 110

Contenido XXIV

Lista de tablas

Pág. Tabla 1.1 : Propiedades globales de la Vía Láctea (Binney and Tremaine 2008)...... 12 Tabla 1.2: Galaxias constituyentes del Grupo Local. 1. Nombre; 3. y 4. Ascensión y declinación 5. y 6. Longitud y latitud galáctica; 7. Tipo de galaxia, (8) Subgrupo al que pertenece dentro del Grupo Local (Mateo, 1998)...... 15 Tabla 1.3: Propiedades del Grupo Local (van den Berg 1999)...... 22 Tabla 1.4: Propiedades de galaxias enanas del Grupo Local – (1) Nombres; (2) Influencia gravitacional, (3) Magnitud absoluta en el azul; (4) Distancia desde la rama horizontal; (5) Tipo morfológico; (6) Índice de marea (Weisz et al., 2011)...... 23 Tabla 1.5: Galaxias satélite de la VL que se encuentran a una distancia menor de 250 kpc, The stromlo Milky Way satellites survey program leader: Helmut Jerjen (anu) Tomada de http://www.mso.anu.edu.au/ 6 de abril de 2015...... 27 Tabla 3.1 : Parámetros galácticos para la Vía Láctea (Frenk and White 2014), (Pawlowski 2011) y (Binney and Tremaine 2008)...... 52 Tabla 3.2 : Los parámetros de los modelos de galaxias objetivo y proyectil de disco (Pawlowski, 2011)...... 58 Tabla 3.3 : Tabla comparativa entre los datos usados por Pawlowsky y los datos obtenidos por el modelo usado por este estudio, se nota que la diferencia porcentual media entre los datos es de 4.74% lo cual es un buen indicativo para la utilización de este modelo...... 59 Tabla 3.4 : Se muestran los parámetros de las diferentes galaxias satélites utilizadas en este estudio, calculados con la metodología planteada en los puntos 3.3.1, 3.3.2 y 3.3.3. .. 60 Tabla 3.5 : Parámetros de las diferentes galaxias satélite utilizadas en este estudio para las orbitas de tipo parabólico en simulaciones de baja resolución...... 70 Tabla 3.6 : Tabla en la que se muestran los parámetros de las diferentes galaxias satélite utilizadas en este estudio para las orbitas de tipo parabólico planteadas...... 74 Tabla 4.1 : Número de grupos respecto a la distancia de enlace elegida...... 95 Tabla 4.2 : Número de grupos respecto a la distancia de enlace elegida...... 103

Introducción

Desde la antigüedad la humanidad ha sentido una fuerte curiosidad por el cielo y los astros, una muestra de esto se refleja en los monumentos celebres como Stonhenge en Gran Bretaña, la pirámide de Keops en Egipto o el parque arqueológico de San Agustín en Colombia. En la antigüedad se destacaron grandes hombres como Aristarco de Samos con su teoría heliocéntrica, Eratóstenes que midió el perímetro de la tierra y además construyo uno de los primeros astrolabios, entre otros grandes pensadores, estos fueron los inicios de la rama del conocimiento que hoy llamamos astronomía (Longair 2006).

Ya entrada la edad media los eruditos de la época como Nicolás Copérnico, Johannes Kepler entre otros, empezaron a desarrollar modelos geométricos a partir de las observaciones para tratar de predecir con mayor precisión el movimiento de los cuerpos celestes, ya teniendo este conocimiento en su mano el hombre tuvo que dar el siguiente gran paso y fue allí cuando apareció el genio de Isaac Newton quien en su libro “Principios matemáticos de la filosofía natural” postuló la ley de gravitación universal y las leyes del movimiento, fundamentos de lo que se conoce como mecánica clásica, además Newton descubrió los principios básicos de la óptica especialmente a lo que concierne al espectro de la luz blanca, de esta forma se forjaron los principios básicos en los cuales se cimienta la Astrofísica. Ya entrado el siglo XIX J. von Frauhofer analizando el espectro de la luz solar, descubrió que ciertas propiedades de los astros pueden ser estudiadas a través de las características del espectro electromagnético, en este momento surge la Astrofísica moderna. (Longair 2006)

Es claro que la Astrofísica evoluciona de forma paralela a como lo hacen las teorías físicas y la tecnología en cada una de las épocas de la historia, en la actualidad con el gran desarrollo de los instrumentos y técnicas de observación se ha llegado a entender el universo como un recinto donde se encuentran las más variadas condiciones de campo, energía y materia. Condiciones que en la cotidianidad científica terrestre no son de ninguna forma alcanzables. Dadas estas condiciones se puede decir que la astrofísica es la rama de la ciencia que estudia el comportamiento y evolución de los astros a la luz de las teorías físicas, entre las se encuentran la relatividad general, la mecánica cuántica, la física de partículas, entre otras.

2 Introducción

Sin embargo en la astrofísica se presentan un gran número de problemas sin resolver, entre estos se hallan la alta temperatura de la corona solar, la existencia de planetas alrededor de otros soles, el nacimiento de estrellas masivas, el origen de los rayos ultra energéticos, los estallidos de rayos gamma, agujeros negros masivos galácticos, la materia oscura, la energía oscura, el problema del origen de las galaxias satélite de la VL, entre otros. En el caso de las galaxias satélites estas se caracterizan por ser galaxias de tipo enana, mejor conocidas como galaxias esferoidales. La mayoría de estas se encuentran distribuidas en un disco perpendicular al plano de la galaxia; (Kroupa et al., 2005) este disco es conocido con el nombre de disco de satélites (DoS, del inglés Disc of Satellites). Encontrar un modelo que dé cuenta de la cantidad y de la distribución espacial de estas galaxias no ha sido posible, sin embargo, se han presentado varias propuestas de solución, una de las cuales propone que estas tuvieron origen en la colisión de la VL con otra galaxia de disco hace ya miles de millones de años (Okazaki and Taniguchi (2000)).

En la Cosmología, las galaxias son los constituyentes más básicos, ya que ayudan a la construcción de estructuras a gran escala en el universo, estas se acumulan en grupos o cúmulos de diferentes formas y constituyentes. En este caso la VL no es la excepción ya que esta forma parte de un grupo denominado el Grupo Local (LG), que al mismo tiempo forma parte del supercúmulo de Virgo. Las galaxias enanas son estructuras estelares pequeñas, presentan poco brillo superficial y poca cantidad de estrellas con respecto a su masa total. Su contenido metálico es bajo debido a que se componen básicamente de estrellas viejas. Se clasifican morfológicamente en dos grupos: las irregulares dIrrs (Dwarf Irregular), donde aún se forman estrellas, y las enanas esferoidales o dSphs (Dwarf Spheroidals) conocidas por ser la población más dominante debido a su cantidad.

La gran importancia de las galaxias enanas radica en que son la población más abundante en el universo local, razón por la cual se convierten en herramienta fundamental para entender la formación y evolución de galaxias de mayor tamaño. Estas galaxias enanas poseen un número de estrellas bajo con respecto a su masa total y esa falta de materia luminosa se ha interpretado hasta ahora como un halo de materia oscura envolvente. (Lokas 2008, 2011). Las galaxias satélites de la VL se encuentran en promedio orientadas sobre un plano que presenta una inclinación de 88° con respecto al disco de la galaxia (Metz et al., 2007). De la distribución de galaxias enanas satélite cercanas a la VL en la actualidad no se conoce con claridad su origen, su formación y distribución.

El estudio de las galaxias enanas también se hace relevante por su alto contenido de estrellas viejas que tienen una metalicidad baja y que pueden dar información de propiedades del universo en sus inicios, es por esta razón que el estudio de las dSph satélites de la VL es de gran valor. Además, su importancia también radica en que la distribución de galaxias satélite de la VL contradicen las predicciones hechas por el modelo cosmológico de la materia oscura fría. Esta teoría predice que la distribución de Introducción 3 subestructuras dentro de halos como el de la VL debería ser isotrópica, lo cual no coincide con las observaciones, (Kroupa et al. 2005) mostraron que la distribución de estos objetos es completamente anisotropica.

En el estudio del origen de las galaxias satélite de la VL el grupo de Astrofísica de la Universidad Nacional de Colombia la estudiante Diana Judith Cubillos Jara realiza un estudio titulado Estudio de la Formación de Galaxias Enanas Esferoidales Satélites de la VL Mediante la Simulación de la Colisión de dos Galaxias. En este trabajo se realizó una excelente revisión bibliográfica respecto al tema de investigación referente importante para el presente estudio, así como también una simulación de la colisión de dos galaxias de disco para una relación de masa entre la VL y la Galaxia Fantasma de 25%, además se tuvieron en cuenta parámetros de impacto de 6, 6.5 y 7 kpc, en este trabajo se encontró que la relación de masa del 25% genera un calentamiento bastante elevado en el disco de la VL, para finalizar esta disertación recomienda realizar simulaciones con relaciones de masa y parámetros de impacto diferentes.

El presente trabajo tiene como objeto realizar simulaciones numéricas de N-cuerpos de colisiones entre la VL (Vía Láctea) y la hipotética galaxia de disco (Galaxia fantasma) que pudieron haber dado origen al disco de satélites de la VL. Teniendo en cuenta lo expuesto anteriormente se realizará el estudio del origen de las galaxias enanas satélite de la VL mediante la simulación de la colisión de dos galaxias de disco, este trabajo se desarrolla primero definiendo las características morfológicas, de masa y de orientación espacial correspondientes a las enanas del Grupo Local, a continuación se mencionaran los referentes teóricos que permiten la construcción de los objetos de estudio, así como comprender la dinámica de su interacción, a continuación se proceden a determinar las condiciones iniciales para llevar a cabo las simulaciones de N-cuerpos. Después de hacer el análisis de las simulaciones realizadas, se revisa la distribución espacial y de número de partículas encontrados en los remanentes de la colisión y finalmente se concluye sobre qué características debe tener la (Galaxia fantasma) para que la colisión de dos galaxias de disco sea la posible causa que dio origen al disco de galaxias satélite de la VL.

1. Galaxias, Grupo Local y Galaxias Satélites

En este capítulo se exponen las características más generales de las galaxias, algunas de las características más relevantes de la (VL), la ubicación de la misma dentro de un grupo denominado como Grupo Local de galaxias, así mismo se pondrá en consideración las características de sus galaxias enanas satélites, su disposición espacial y su posible origen, por último se menciona el modelo definido para el presente estudio.

1.1 Galaxias

Una galaxia es un grupo de gas, polvo, cúmulos de estrellas, nebulosas y planetas unidos a través de su fuerza gravitatoria (Sparke & Gallagher 2000), etimológicamente hablando la palabra galaxia deriva del griego ( γαλαξίας ) que significa (lechoso). Las galaxias según su número de estrellas se pueden clasificar en galaxias enanas las cuales tienen un máximo de 10 8 estrellas y las galaxias gigantes que tienen hasta 10 14 estrellas, estas siempre orbitando respecto al centro de masa del sistema ( Uson, Boughn and Kuhn, 1990) . En la Cosmología las galaxias son los constituyentes más básicos del cosmos, ya que permiten la construcción de estructuras a gran escala en el universo, estas se acumulan en grupos o cúmulos de diferentes formas y constituyentes (Binney and Tremaine. 1994).

El conteo de galaxias asciende a un orden de 1.7 x10 11 galaxias en el universo conocido (Gott et al. 2005), en general las galaxias se encuentran organizadas en el espacio debido a los efectos gravitatorios producidos por los objetos circundantes, estas organizaciones generalmente se les llama asociaciones, que pueden ser grupos, cúmulos o supercúmulos (Kravtsov 2012). En este caso la VL no es la excepción ya que esta forma parte de un grupo denominado el Grupo Local LG, y este a su vez forma parte del supercúmulo de Virgo (Tully et al 2014 ). Las galaxias dependiendo de su morfología se clasifican según la secuencia de Hubble.

6 Colisión de dos galaxias de Disco

1.1.1 Secuencia de Hubble

La secuencia de Hubble fue inventada por él en el año de 1936 (Hubble 1936), con el objetivo de clasificar las galaxias según su morfología como se aprecia en la figura 1.1, dentro de esta secuencia se clasifican las galaxias en tres clases: elípticas (E), lenticulares (S0) y las espirales (S), en la actualidad se emplea una cuarta clase que hace referencia a las galaxias de tipo irregular (Irr).

Figura 1.1: Secuencia de Hubble. Las galaxias elípticas (E) vienen acompañadas por un número que hace referencia a la excentricidad de la galaxia, en las galaxias espirales (S) las letras a, b y c hacen referencia a la separación de los brazos, la letra a designa los brazos más cerrados y la letra c hace referencia a los brazos más separados, las galaxias espirales barradas (SB) las letras a, b y c hacen referencia a la longitud de la barra, la letra a designa barras cortas y la letra c barras más largas. Tomado de http://es.wikipedia.org/wiki/Secuencia_de_Hubble. 1 de abril de 2015.

Se debe tener presente que en el caso de la VL esta se presenta como una galaxia espiral barrada (SB), más adelante se hablara de las características fundamentales de la galaxia. A continuación se describe de manera cualitativa algunos de los rasgos de los diferentes tipos de galaxias.

1.1.2 Galaxias elípticas

Las galaxias de tipo elíptico presentan una forma esferoidal o una forma elipsoidal y no presentan de ninguna forma estructura de tipo espiral. En estas se encuentra poco gas, poco polvo y pocas estrellas jóvenes. En la secuencia de Hubble se clasifican desde E0 hasta E7 esto dependiendo del grado de aplastamiento de la estructura (Rodríguez 2001) (Binney and Merrifield 1998). Galaxias, Grupo Local y galaxias satélites 7

Estas se pueden denotar con la notación E p donde p es un entero que varía entre 1 y 7 el cual satisface la Ecuación (1.1)

 = 10 1 −  (1.1) ͤ ʠ ʡ

Siendo a y b los ejes mayor y menor respectivamente, como ejemplo, en el caso de una galaxia completamente esférica se tiene que a=b=1, si se aplica la ecuación (1.1) se tiene que p=0 por lo tanto este objeto se clasifica como una galaxia tipo E0. Ya en el caso más extremo si la estructura es demasiado aplastada se tendría que b/a=0,3, esto indicaría que es una galaxia tipo E7. (Rodríguez 2001, Binney and Merrifield 1998, Karttunen et al 1996).

Figura 1.2 Galaxia elíptica M87 Tomado de https://www.astro.virginia.edu/

El tamaño de estas estructuras varía considerablemente, las más grandes se les denomina galaxias elípticas supergigantes, esto es debido a que sus diámetros pueden llegar a tener miles de kpc de longitud, mientras que las galaxias más pequeñas son llamadas galaxias elípticas enanas, cuyos diámetros solo ascienden a las decenas de kpc de longitud. Las 12 galaxias supergigantes pueden llegar a tener masas hasta de 10 M⊙ y las enanas llegan 6 a tener masas del orden de 10 M⊙. Estas galaxias enanas son el tipo más común de galaxias. (Rodríguez 2001, Binney and Tremaine. 1994 y Binney and Merrifield 1998). 8 Colisión de dos galaxias de Disco

1.1.3 Galaxias espirales

Las galaxias espirales reciben su nombre debido a su estructura, esta presenta unos brazos de tipo espiral que se envuelven alrededor de un bulbo central. Su distribución geométrica es de tipo plano, la zona más brillante se encuentra en el centro de la estructura que se compone de un bulbo que tiene forma elipsoidal. Como ejemplo típico de estas galaxias tenemos la VL (MW) y la galaxia de Andrómeda (M31). En la zona del bulbo se encuentran estrellas que se mueven de forma aleatoria alrededor del centro de masa, esta región presenta una forma aproximadamente esférica debido al equilibrio que se presenta debido al movimiento aleatorio y la fuerza de atracción gravitatoria. El disco está compuesto por gas, polvo y estrellas, estas también se mueven respecto al centro de masa, en orbitas aproximadamente circulares, no obstante su movimiento se encuentra restringido al plano del disco. (Rodríguez 2001), (Binney and Tremaine. 1994).

Figura 1.3: Galaxia espiral M31 Tomado de www.bbc.co.uk

Entre las galaxias espirales se encuentran algunas diferencias en el tamaño de su núcleo o en la disposición de sus brazos de tipo espiral, estas diferencias son tenidas en cuenta al momento de su clasificación en la secuencia de Hubble. Casi en la mitad de las galaxias espirales se encuentra una barra que atraviesa el núcleo, estas galaxias son las denominadas espirales barradas (Binney and Tremaine. 1994).

En cuanto a las dimensiones se refiere las galaxias espirales tienen diámetros que van 10 desde los 7 kpc hasta los 35 kpc aproximadamente, y sus masas varían entre 10 M⊙ y 12 10 M⊙ , la VL y Andrómeda están entre las más grandes y masivas (Rodríguez 2001).

También se encuentran algunas galaxias que cuentan con bulbo y disco pero no muestran presencia de brazos espirales, en la secuencia de Hubble son las denominadas galaxias lenticulares . En estas estructuras el bulbo es demasiado grande y se cree que estos objetos Galaxias, Grupo Local y galaxias satélites 9

presentan una transición entre galaxia elíptica a galaxia espiral. (Rodríguez 2001, Graham and Worley 2008 y Karttunen et al 1996).

1.1.4 Galaxias irregulares

Fueron llamadas irregulares por Hubble, ya que no presentaban ningún tipo de característica morfológica (esferoidal o circular) o dinámica (rotacional) similar a las galaxias ya estudiadas, de hecho estas galaxias presentan una estructura completamente desorganizada, al parecer muchas de estas galaxias están pasando por un ciclo de formación estelar intenso, para este tipo de galaxias Hubble definió dos clases, Irr I e Irr II , en la primera clase las galaxias presentan perfiles asimétricos y no presentan ningún abultamiento central o estructura espiral definida, presentan una gran cantidad de grupos individuales de estrellas jóvenes; en la segunda clase se presentan un poco más suavizadas, aunque no se definen con claridad estrellas individuales o algún tipo de cumulo estelar (Rodríguez 2001, Longair, M. S. 1998 y Karttunen et al 1996).

Figura 1.4: Galaxia irregular la Gran Nube de Magallanes Tomado de http://www.surastronomico.com/

1.1.5 Galaxias enanas

Las galaxias enanas como su nombre lo indica son estructuras pequeñas comparadas con las galaxias regulares, estas son galaxias de baja luminosidad, se consideran la población de galaxias más dominantes del universo actual y seguramente eran mucho más numerosas en épocas pasadas, además si estas estructuras se encuentran dominadas por 10 Colisión de dos galaxias de Disco

materia oscura (DM) como se piensa que ocurre con las del Grupo Local, tal vez estas galaxias podrían ser responsables de una gran parte de la masa de los cúmulos de galaxias y quizá de todo el universo (Mateo, 1998).

Las galaxias enanas de baja luminosidad tienden a ser pobres en metales, siendo estas las estructuras galácticas más simples conocidas. No obstante “simple” es un término relativo ya que la historia de formación estelar y enriquecimiento químico de estas galaxias es compleja, variada, y en la mayoría de los casos provocada y sostenida por mecanismos que aún no son muy claros (Mateo, 1998). Estas galaxias se clasifican básicamente en dos tipos: enanas esferoidales (dSph) y enanas irregulares (dIrr), usando los parámetros establecidos por la secuencia de Hubble (Tolstoy, 2003).

Figura 1.5: Galaxia enana de Fornax tomado de www.abc.es.

Las galaxias enanas pueden estar entre las galaxias individuales más oscuras conocidas. Ellas juegan un papel importante en la solución del problema de la materia oscura (DM), que ya ha puesto restricciones interesantes sobre la naturaleza y la distribución de la (DM), e incluso si el paradigma DM es válido para estos sistemas (Mateo, 1998). Sin embargo para resolver esta cuestión se dice que las galaxias enanas se caracterizan por tener un número de estrellas muy bajo con respecto a su masa total, lo cual justifica su baja luminosidad, y siendo esto interpretado como un halo de materia oscura que las envuelve. También partiendo de sus relaciones de masa-luminosidad que superan en un orden de magnitud a la de los cúmulos globulares, esto podría indicar que la materia oscura es la que domina estas estructuras (Da Costa, 1999).

Las galaxias enanas son una gran fuente de información para la comprensión de la formación y evolución del Grupo Local de galaxias, ya que la mayoría de galaxias satélite de las galaxias nodrizas de este grupo son de este tipo. Galaxias, Grupo Local y galaxias satélites 11

1.1.6 Vía Láctea (VL)

Se debe tener presente que el concepto de galaxia o de universos islas es relativamente nuevo, si bien este fue planteado por el filósofo Immanuel Kant en el siglo XVIII, fue solo hasta el siglo XX que esta idea se aceptó como una realidad. Esta discusión que llevaba ya varias décadas y en la cual habían realizado aportes grandes pensadores como: William Hershel, quien trabajó un tiempo en la idea, que después abandono; William Parsons quien retoma estas opiniones y observa la galaxia espiral M51 conocida como El Remolino, que se logró apreciar usando un telescopio de grandes proporciones llamado el Leviatán de Parsons, sin embargo no era claro si estos objetos pertenecían a la VL o estaban más allá de ella; ya en 1912 aparece en escena Henrieta Leavitt que descubre la relación periodo- luminosidad de las llamadas estrellas variables cefeidas que permitirían determinar la distancia a la cual se encuentran algunos objetos celestes. Comprendiendo esta relación el astrónomo Harlow Shapley (Shapley and Curtis 1921) estimó las dimensiones de la VL, y de un valor estimado de 10 kpc de diámetro este lo llevo a un valor de aproximadamente 100 kpc valor que era exageradamente alto para las dimensiones reales de la galaxia, esto llevo a una gran discusión que fue resuelta por el descubrimiento de una cefeida en la llamada nebulosa de Andrómeda por Edwin Hubble, este determinó que la distancia respecto a nosotros es del orden de 300 kpc, confirmando así que Andrómeda era una galaxia isla y no pertenecía a la VL. (Uribe 2001).

Nuestra galaxia recibe su nombre ya que en las noches bajo un cielo despejado se parece a una mancha lechosa que se extiende en ambas direcciones de la bóveda celeste. La VL es una galaxia que presenta una morfología de galaxia espiral barrada gigante (Binney and Tremaine. 1994). Esta galaxia se encuentra formada por cuatro componentes; el halo de materia oscura, el disco, el bulbo central y el halo estelar. El halo de materia oscura es una estructura de simetría aproximadamente esférica que domina la masa del sistema, que no es visible (Karttunen, Kroger, Poutanen and Donner 2003). El disco es una estructura plana que rota respecto al centro de masa, contiene la mayoría de estrellas jóvenes de la galaxia y en general la masa de polvo estelar y de gas de la galaxia (Karttunen, Kroger, Poutanen and Donner 2003). 12 Colisión de dos galaxias de Disco

Figura 1.6: VL visión artística. Tomada de http://apod.nasa.gov/

En cuanto al bulbo que se encuentra en la zona central de la VL es la zona de mayor densidad de estrellas de la galaxia. Esta región no se puede apreciar por la gran concentración de polvo y las imágenes se presentan de forma difusa. (Grebel 1998), el brillo superficial y la dispersión de velocidades es consecuente con la presencia de un 6 agujero negro en el centro que se estima tiene una masa de (3,9±0,3)x10 M⊙ (Binney and Tremaine 2008).

Se estima que la VL tiene de entre 10 9 y 10 12 estrellas, también cuenta con alrededor de 10 11 planetas (Cassan et al. 2012).

Propiedades globales (VL)

Longitud de escala radial del disco R d (2,5±0,5) kpc 10 Luminosidad del disco (2,5±1)x10 L⊙ 9 Luminosidad del bulbo (5±2)x10 L⊙ 9 Luminosidad total (3,0±1)x10 L⊙ 10 Masa del disco (4,5±0,5)x10 M⊙ 9 Masa del bulbo (4,5±1,5)x10 M⊙ 12 Masa del halo de materia oscura (Entre 2 -1,8 y 2+3,0)x10 M⊙ Halo oscuro radio de masa media (Entre 100 -80 y 100+100) kpc

Razón masa luminosidad disco ΥR (1,8±0,7) Υ⊙

Razón masa luminosidad total ΥR (Entre 70 -63 y 70+100) Υ⊙ 6 Masa agujero negro central (3,9±0,3)x10 M⊙ Clasificación de Hubble Sbc

Tabla 1.1: Propiedades globales de la Vía Láctea (Binney and Tremaine 2008). Galaxias, Grupo Local y galaxias satélites 13

Las estrellas y gases que constituyen el disco de la VL giran alrededor del centro galáctico a una velocidad aproximada de 220 km/s, esta velocidad de rotación circular constante va en contradicción de la dinámica de rotación Kepleriana, lo cual parece indicar la existencia de una materia no visible denominada materia oscura. (Imamura 2006 y Zwicky 1937).

Además, la VL se encuentra rodeada por unas galaxias pequeñas llamadas galaxias satélites (Metz, Kroupa, Jerjen. 2006), estas son parte central de este trabajo, se estudiará más adelante el tipo, disposición espacial y posible origen de estas estructuras. La VL forma parte del llamado Grupo Local de galaxias que a su vez es subestructura del supercúmulo de virgo, y este a su vez es subcomponente de supercúmulo de Laniakea (Brent Tully 2014).

1.2 Grupo Local

El termino Grupo Local fue acuñado por Edwin Hubble en el año de 1936, en el libro que lleva por nombre The Realm of the Nebulae (Hubble 1936, pp. 124–151) en este documento Hubble identificó algunos objetos que eran posibles miembros de lo que el denominaba Grupo Local, entre los objetos más importantes se encuentran la Vía Láctea, Andrómeda (M31), Triangulum (M33), la Gran Nube de Magallanes entre otros. En la actualidad la cifra de miembros del grupo ha ido en aumento y de los doce encontrados por Hubble, para el año 2003 ya había aumentado a alrededor de treinta y seis (Van den Bergh 2003).

El Grupo Local es una reunión de galaxias compuesta por aproximadamente 36 galaxias, la mayoría de las cuales son galaxias de tipo enana. Las más importantes son dos galaxias espirales gigantes llamadas Vía Láctea y Andrómeda (M31), estas tienen grupos más pequeños de galaxias llamadas sistema de galaxias satélites (Karachentsev and Kashibadze 2006). En estudios de los últimos años se han estimado los valores de las masas de esta

agrupación donde se reporta un radio de masa medio de R LG = 350 kpc y una masa total 12 de M LG = (2,3±0,6) x 10 M⊙ (Van den Bergh 1999). Este grupo forma parte del supercúmulo de Virgo.

14 Colisión de dos galaxias de Disco

Figura 1.7: Esquema representativo de la disposición espacial de algunos de los objetos del Grupo Local. Tomada de http://www.australianscience.com.au

En la tabla 1.2 se muestran un registro de las galaxias que de manera potencial podrían pertenecer al Grupo Local, en esta se encuentran las galaxias caracterizadas por su morfología, incluyendo las que se encuentran en estado de transición (Mateo 1998). Esta tabla se emplea por que representa uno de los primeros y más completos estudios respecto a las galaxias que constituyen el Grupo Local. Se debe resaltar que en la actualidad el número de galaxias que pertenecen a este grupo ha aumentado considerablemente, debido al reciente descubrimiento de las galaxias ultra tenues (UFDG) (Theis, Metz, Kroupa, Jerjen. 2009). Estas galaxias se caracterizan por tener poco luminosidad en comparación a sus compañeras de mayor brillo.

Galaxias, Grupo Local y galaxias satélites 15

Galaxia Otro nombre α2000 δ2000 l b Tipo Subgrupo (1) (2) (3) (4) (5) (6) (7) (8) WLM DDO 221 00 01 58 -15 27.8 75.9 -73.6 IrrIV -V LGC NGC 55 00 15 08 -39 13.2 332.7 -75.7 IrrIV LGC IC 10 UGC 192 00 20 25 +59 17.5 119.0 -3.3 dIrr M31 NGC 147 DDO 3 00 33 12 +48 30.5 119.8 -14.3 dSph/dE5 M31 And III 00 35 17 +36 30.5 119.3 -26.2 dSph M31 NGC 185 UGC 396 00 38 58 +48 20.2 120.8 -14.5 dSph/dE3p M31 NGC 205 M 110 00 40 22 +41 41.4 120.7 -21.1 E5p/dSph -N M31 M32 NGC 221 00 42 42 +40 51.9 121.2 -22.0 E2 M31 M31 NGC 224 00 42 44 +41 16.1 121.2 -21.6 SbI -II M31 And I 00 45 43 +38 00.4 121.7 -24.9 dSph M31 SMC NGC 292 00 52 44 -72 49.7 302.8 -44.3 IrrIV-V MW Sculptor 01 00 09 -33 42.5 287.5 -83.2 dSph MW LGS 3 Pisces 01 03 53 +21 53.1 126.8 -40.9 dIrr/dSph M31 IC 1613 DDO 8 01 04 54 +02 08.0 129.8 -60.6 IrrV M31/LGC And II 01 16 27 +33 25.7 128.9 -29.2 dSph M31 M 33 NGC 598 01 33 51 +30 39.6 133.6 -31.3 ScII -III M31 Phoenix 01 51 06 -44 26.7 272.2 -68.9 dIrr/dSph MW/LGC Fornax 02 39 59 -34 27.0 237.1 -65.7 dSph MW EGB 0427+63 UGCA 92 04 32 01 +63 36.4 144.7 +10.5 dIrr M31 LMC 05 23 34 -69 45.4 280.5 -32.9 IrrIII -IV MW Carina 06 41 37 -50 58.0 260.1 -22.2 dSph MW A DDO 69 09 59 24 +30 44.7 196.9 +52.4 dIrr MW/N3109 Sextans B DDO 70 10 00 00 +05 19.7 233.2 +43.8 dIrr N3109 NGC 3109 DDO 236 10 03 07 -26 09.5 262.1 +23.1 IrrIV -V N3109 Antlia 10 04 04 -27 19.8 263.1 +22.3 dIrr/dSph N3109 DDO 74 10 08 27 +12 18.5 226.0 +49.1 dSph MW Sextans A DDO 75 10 11 06 -04 42.5 246.2 +39.9 dIrr N3109 Sextans 10 13 03 -01 36.9 243.5 +42.3 dSph MW Leo II DDO 93 11 13 29 +22 09.2 220.2 +67.2 dSph MW GR 8 DDO 155 12 58 40 +14 13.0 310.7 +77.0 dIrr GR8 Ursa Minor DDO 199 15 09 11 +67 12.9 105.0 +44.8 dSph MW Draco DDO 208 17 20 19 +57 54.8 86.4 +34.7 dSph MW Milky Way 17 45 40 -29 00.5 0.0 0.0 Sbc MW Sagittarius 18 55 03 -30 28.7 5.6 -14.1 dSph -N MW SagDIG UKS 1927 -177 19 29 59 -17 40.7 21.1 -16.3 dIrr LGC NGC 6822 DDO 209 19 44 56 -14 48.1 25.3 -18.4 IrrIV -V LGC DDO 210 Aquarius 20 46 46 -12 51.0 34.0 -31.3 dIrr/dSph LGC IC 5152 22 02 42 -51 17.7 343.9 -50.2 dIrr LGC Tuscana 22 41 50 -64 25.2 322.9 -47.4 dSph LGC UKS2323 -326 UGCA 438 23 26 27 -32 23.3 11.9 -70.9 dIrr LGC Pegasus DDO 216 23 28 34 +14 44.8 94.8 -43.5 dIrr/dSph LGC

Tabla 1.2: Galaxias constituyentes del Grupo Local. 1. Nombre; 3. y 4. Ascensión y declinación 5. y 6. Longitud y latitud galáctica; 7. Tipo de galaxia, (8) Subgrupo al que pertenece dentro del Grupo Local (Mateo, 1998). 16 Colisión de dos galaxias de Disco

1.2.1 Distribución espacial del Grupo Local

Respecto a la distribución espacial del Grupo Local (LG) se han realizado varios trabajos, entre los que destacan el realizado por van den Bergh en el año 1999, este afirma que el Grupo Local es un clúster de tipo binario, ya que de manera general la mayoría de objetos se ven dominados por las dos influencias gravitatorias dominantes del grupo, estas son las correspondientes al subgrupos de la VL y de Andrómeda (M31).

Figura 1.8: La figura muestra que la mayoría de los miembros del LG se concentran en dos subgrupos que se centran en la VL y M31. El círculo de radio de 1180 kpc corresponde a la superficie de la velocidad cero del LG. El círculo de trazos con el radio 450 kpc muestra el radio de la esfera que contiene la mitad de todas las galaxias del Grupo Local (van den Bergh, 1999)

Otro trabajo relevante es el realizado por Sawa y Fujimoto en el año 2005, en el cual se presenta la distribución de los objetos en el cielo, revelando que los objetos que constituyen el Grupo Local no están ordenados aleatoriamente, sino al contrario presentan una distribución tipo anillo de un diámetro que asciende a los 180° de arco.

En la Figura 1.9 se aprecia la distribución bidimensional de los miembros del LG, la Corriente de Magallanes, y las colas de marea de Sagitario en la proyección del martillo en el cielo. Los datos de los miembros LG se deben a Mateo (1998), y no se dan aquí las Galaxias, Grupo Local y galaxias satélites 17

líneas de las longitudes y latitudes galácticas con el fin de enfatizar que estos objetos no se distribuyen de forma aleatoria, sino de alguna manera sistemática, como en un anillo de gran diámetro, excepto cuatro enanas en la dirección del centro galáctico (Sawa y Fujimoto 2005).

Figura 1.9: Distribución de los miembros del Grupo Local, en esta sugiere la existencia de un orden espacial de tipo anillo en el Grupo Local (LG) (Sawa y Fujimoto 2005).

En este mismo trabajo Sawa y Fujimoto (2005) realizan un análisis tridimensional de las ubicaciones de las galaxias del Grupo Local (LG). Aquí se resalta que las galaxias satélite de la VL (MW) y de Andrómeda (M31) se encuentran distribuidas en estructuras tipo disco y estas se ven cercanamente paralelas aunque forman un pequeño ángulo entre las dos. Los miembros de la VL y los de Andrómeda (M31) se ven alineados respecto a sus galaxias nodrizas, los de la VL un poco más que los de (M31) sin embargo la alineación entre los objetos de cada grupo no es perfecta ver Figura 1.10.

18 Colisión de dos galaxias de Disco

Figura 1.10: Distribución espacial tridimensional de los miembros del Grupo Local (Sawa and Fujimoto, 2005) vista desde la dirección ( l,b ) = (296° ,−11°). Esta es una vista de canto del plano orbital de la galaxia VL y de la M31. Los círculos rellenos oscuros denotan los miembros del grupo de la VL, los grises son los miembros de M31 y los círculos sin relleno no pertenecen a ninguno de los dos grupos. Se puede ver como los miembros del Grupo Local están distribuidos en una estructura tipo disco con un espesor finito de 50 – 100 kpc. El tamaño de cada círculo representa su brillo cualitativamente.

Lo que se aprecia en la Figura 1.10 también se ve en la figura 1.11, sino que visto desde la dirección (l, b) = (206°, -11°). Esta es una toma vista de frente en el plano del LG, o del plano orbital de la galaxia M31. El grupo de la VL y los miembros del grupo (M31), y los miembros fuera de estos grupos parecen formar un anillo fragmentario cuyo diámetro es ~ 800kpc y más (Sawa y Fujimoto 2005). Galaxias, Grupo Local y galaxias satélites 19

Figura 1.11: Esta es una toma vista de frente en el plano del LG, o del plano orbital de la galaxia M31. (Sawa y Fujimoto 2005).

1.2.2 Propiedades del Grupo Local

En estudios realizados por van den Bergh en 1999 se encuentran tres características fundamentales de las galaxias del Grupo Local, primero el comportamiento de la distribución de su luminosidad, segundo la relación existente entre la luminosidad M V y la metalicidad (FE/H) y por último la relación entre la longitud de escala del disco y la luminosidad.

Este encontró que el número de objetos del (LG) crece a medida que la luminosidad de los objetos viaja hacia magnitudes más bajas, donde el descenso del número de objetos en la

última parte de la Figura 1.12 para cuando las magnitudes son más débiles que M V = -10 se debe a la falta de datos.

20 Colisión de dos galaxias de Disco

Figura 1.12: Distribución Luminosidad para (LG). La figura muestra que las galaxias irregulares del LG tienen una función de luminosidad más baja que la de las galaxias dSph del Grupo Local. (Van den Bergh, 1999).

En segunda instancia se encuentra una relación entre la luminosidad y la metalicidad para las galaxias que constituyen el grupo como se aprecia en la Figura 1.12.

Figura 1.13: Relación entre luminosidad M V y la metalicidad [Fe/H] para las galaxias que pertenecen al Grupo Local (LG). Se muestra la línea de correlación que viene dada

por la expresión M V = 20 − 5[Fe/H] (van den Bergh, 1999).

En la Figura 1.12 se puede apreciar que la dependencia de la metalicidad y la luminosidad puede depender de que galaxias con masas muy pequeñas no puedan retener los metales formados en las explosiones de supernovas, lo que aparentemente se justifica con que pocos cúmulos de tipo globular muestran enriquecimiento interno, ya que la 5 mayoría de los metales se disipan en sistemas con masas inferiores a 10 M⊙ (van den Galaxias, Grupo Local y galaxias satélites 21

Bergh, 1999). En esta Figura se aprecia una anomalía en la galaxia M32 o le Gentil, ya que su metalicidad es baja y no coincide con la tendencia, esto se puede deber a que en el pasado una parte de su masa pudo ser desprendida por interacciones con Andrómeda (M31).

Por último y no por esto menos importante, van der Bergh encuentra la relación entre la longitud de escala del disco y su luminosidad, en este punto se encuentra que la mayoría de galaxias con gran longitud de escala presentan mayor luminosidad, esto al parecer por consecuencia del constante cambio de la luminosidad (van den Bergh, 1999) Ver Figura 1.14.

Figura 1.14: Ilustración de la luminosidad versus la longitud de escala del disco para las galaxias del Grupo Local. Las galaxias irregulares (cuadrados rellenos) están por encima de la línea de tendencia, las espirales (signo +) y las enanas esferoidales (círculos). Tenga en cuenta que parece que no hay diferencias sistemáticas entre galaxias dSph en el subgrupo M 31 (círculos negros) y los del subgrupo (MW) (círculos abiertos). La línea de

tendencia en la figura es la relación M V = -15-7 logaritmo de longitud de escala (van den Bergh, 1999).

A continuación se presenta un pequeño resumen de las propiedades más importantes del Grupo Local (LG) ver tabla 1.3.

22 Colisión de dos galaxias de Disco

Propiedades globales Grupo Local LG

Radio superficie velocidad cero Ro=1,18±0,16 Mpc

Radio de masa media Rh=350 kpc 10 Luminosidad total LG LV=4,2x10 L⊙

Dispersión de velocidades radial σr=61±8 km/s 12 Masa total LG M=(2,3±0,6)x10 M⊙

Razón masa luminosidad LG M/L V=44±14 12 Masa subgrupo de Andrómeda M(A)=1,15 -1,5x10 M⊙ 10 Luminosidad subgrupo Andrómeda LV=3,0x10 L⊙

Subgrupo Andrómeda M/L V M/L V=38 -50 12 Masa del subgrupo VL M(G)=0,46 -1,25x10 M⊙ 10 Luminosidad subgrupo VL LV=1,1x10 L⊙

Subgrupo VL M/L V M/L V=42 -144 Tabla 1.3: Propiedades del Grupo Local (van den Berg 1999).

1.2.3 Galaxias enanas del Grupo Local

Como se mencionó con anterioridad las galaxias enanas son las estructuras que se presentan en mayor número dentro del Grupo Local de galaxias (van den Bergh, 1999) ver tabla 1.4.

Las galaxias enanas del Grupo Local cumplen con las características de las galaxias mencionadas en el numeral 1.1.5, llevando esto a la conclusión que las galaxias enanas del Grupo Local en especial las (dSph) tienen una alta composición de estrellas antiguas, con 7 algún contenido de hidrogeno y helio, además estas presentan masas del orden de 10 M⊙ ver figura 1.15.

En el modelo de evolución jerárquico las galaxias enanas tienen un papel fundamental, ya que son estas las que permiten la formación de estructuras mucho más grandes ya sea por mezclas o por acreción de las mismas.

Como se presenta una gran dificultad para detectar nuevos integrantes del Grupo Local, en los últimos años se han incluido técnicas de tipo analítico que tienen como fin identificar las propiedades de las galaxias enanas que potencialmente podrían ser integrantes del Grupo Local. Entre las propiedades más relevantes para la categorización utilizando estas técnicas se tienen, la morfología, distribución espacial, luminosidad y contenido de gas e historia de evolución química (Angus, Kroupa. 2003), (Rocha, Annika, Peter, Bullock. 2011).

Galaxias, Grupo Local y galaxias satélites 23

Galaxia Influencia MB D Type Θ Nombre Gravitacional (Mpc) (1) (2) (3) (4) (5) (6) UrsaMin MW -7.13 0.08 dSph 3.3 LGS3 M31 -7.96 0.61 dTrans 1.7 And V M31 -8.41 0.78 dSph 2.8 Draco MW -8.74 0.09 dSph 3.0 Carina MW -8.97 0.10 dSph 2.7 Leo II MW -9.23 0.20 dSph 1.7 And III M31 -9.30 0.72 dSph 3.5 And II M31 -9.33 0.65 dSph 2.4 Antlia M31 -9.38 1.30 dTrans -0.1 Sculptor MW -9.77 0.08 dSph 2.8 Cetus M31 -10.18 0.77 dSph 0.5 Phoenix MW -10.22 0.41 dTrans 0.8 And VI M31 -10.80 0.83 dSph 1.7 And I M31 -10.87 0.76 dSph 3.7 Leo I MW -10.97 0.25 dSph 1.5 DDO210 M31 -11.09 0.94 dTrans 1.6 Pegasus M31 -11.47 0.95 dI 1.2 SagDIG MW -11.49 1.11 dI -0.3 Fornax MW -11.50 0.14 dSph 2.3 AndVII M31 -11.67 0.94 dSph 2.0 Leo A MW -11.70 0.79 dI 0.1 Sagittarius MW -12.80 0.03 dSph 4.0 Tucana MW -12.94 0.86 dSph -0.1 Sex A MW -13.71 1.30 dI -0.6 Sex B MW -13.88 1.40 dI -0.7 WLM M31 -13.95 0.93 dI 0.3 N185 M31 -14.76 0.61 dSph/dE 3.5 N147 M31 -14.79 0.72 dSph/dE 3.0 IC1613 M31 -15.57 0.74 dI 0.9

Tabla 1.4: Propiedades de galaxias enanas del Grupo Local – (1) Nombres; (2) Influencia gravitacional, (3) Magnitud absoluta en el azul; (4) Distancia desde la rama horizontal; (5) Tipo morfológico; (6) Índice de marea (Weisz et al., 2011).

Empezando con la morfología en general resaltan las galaxias de tipo esferoidal (dSph) y las irregulares (dIrr) (Mateo, 1998), existe una buena cantidad de observaciones que indican que algunas galaxias enanas están en proceso de cambio de galaxias del tipo irregular a galaxias del tipo esferoidal. (Grebel, 1998). Lo cual permite concluir que las galaxias más jóvenes presentan una forma irregular y las galaxias más arcaicas exhiben una morfología irregular (Huctchmeier, Karachentsev, Makarov, 2003).

24 Colisión de dos galaxias de Disco

Figura 1.15 : Escala de masa frecuente para algunas galaxias enanas de la VL en función de la luminosidad (Strigari et al. 2008).

En lo referente a la distribución espacial se ha encontrado que las galaxias enanas del Grupo Local tienen la tendencia de acumularse alrededor de galaxias nodriza y generar subgrupos dentro del Grupo Local, la evidencia de este evento se ve reflejada en la tabla 1.4 en la que las galaxias nodrizas del grupo son la VL y Andrómeda (M31) y las galaxias enanas se disponen en dos grupos de influencia gravitacional.

Figura 1.16 : Metalicidad estelar de estrellas población II versus la luminosidad bariónica para diferentes clases de galaxias enanas como se indica en la leyenda. Galaxias con deficiencia de gas (símbolos rellenos) y ricas en gas (símbolos sin relleno). Galaxias con la misma luminosidad, las antiguas poblaciones de dSphs son más ricas en metales que las de dIrrs. Así, en contraste con dIrrs, las dSphs deben haber experimentado un rápido enriquecimiento temprano. Tenga en cuenta la ubicación de los llamados dSph /dirr galaxias de tipo transición, que combinan propiedades de población de dSphs con la formación estelar en curso y un contenido de gas medible en el diagrama, estas propiedades hacen de las dSphs unas plausibles progenitoras (Grebel, Gallagher, & Harbeck, 2003). Galaxias, Grupo Local y galaxias satélites 25

Respecto a la luminosidad y la historia química se puede apreciar en la figura 1.16 que son factores relevantes para categorizar las galaxias enanas y permitir definir su pertenencia al Grupo Local de galaxias.

El estudio de las galaxias enanas del Grupo Local es de relevancia para poder escudriñar en los orígenes del universo, además estas se presentan como una oportunidad única de poder estudiar este tipo de estructuras ya que como se mencionó anteriormente por su baja luminosidad se hace muy difícil su apreciación en otros grupos de galaxias que se encuentran fuera del Grupo Local (Theis, Metz, Kroupa, Jerjen. 2009).

Asimismo el estudio de las galaxias enanas esferoidales satélites de la VL es de vital importancia ya que se presentan discrepancias en la cantidad de galaxias observadas y las predicciones de los modelos cosmológicos de materia oscura fría (CDM), los cuales predicen una cantidad mucho mayor de galaxias Kroupa et al. (2005).

1.3 Galaxias enanas satélite de la Vía Láctea

Teniendo presente que el Grupo Local tiene una dinámica dominada por la presencia de la VL (MW) y Andrómeda (M31), se ha podido determinar de manera observacional que estas estructuras presentan un conjunto de galaxias menores denominadas galaxias enanas satélite.

Figura 1.17 : Ubicación espacial de algunas de las galaxias satélite de la VL que se encuentran a una distancia menor de 250 kpc. Tomada de http://www.mso.anu.edu.au/ 26 Colisión de dos galaxias de Disco

La mayoría de satélites de la VL se encuentran en la franja de los 10 a los 250 kpc, mostrando una gran acumulación en la región de 10 kpc a 100 kpc ver figura 1.18.

Figura 1.18 : Distribución radial de los satélites de la VL, en la parte inferior se muestra la frecuencia a acumulada, se aprecia la acumulación de satélites en los primeros 100 Kpc (McConnachie and Irwin, 2006).

Nombre Distancia Documento de Descubrimiento (kpc) Canes Major 7.2 Martin et al. 2004, A remnant in Canis Major: the fossil of an in-plane accretion on to the Milky Way

Segue 3 17 Belokurov et al. 2010, Big Fish, Little Fish: Two New Ultra- Faint Satellites of the Milky Way Segue 1 23 Belokurov et al. 2007, Cats and Dogs, Hair and A Hero: A Quintet of New Milky Way Companions Sagittarius 24 Ibata, Gilmore & Irwin, 1994, A dwarf in Sagittarius 1995, Sagittarius: the nearest dwarf galaxy

Segue 2 34.7 Belokurov et al. 2009, The discovery of : a prototype of the population of satellitesof satellites Bootes II 43 Walsh, Jerjen & Willman, 2007, A Pair of Bootes: A New Milky Way Satellite Coma 44 Belokurov et al. 2007, Cats and Dogs, Hair and A Hero: A Quintet of New Milky Way Companions 45 Willman et al. 2005, A New Milky Way Companion: Unusual (SDSSJ1049+5103) or Extreme Dwarf Satellite? Galaxias, Grupo Local y galaxias satélites 27

Bootes III 46 Grillmair 2009, Four New Stellar Debris Streams in the LMC 50.8 - SMC 59.7 - Bootes 60 Belokurov et al. 2006, A Faint New Milky Way Satellite in Bootes Ursa Minor 66 A.G. Wilson of the Lowell Observatory in 1955,Sculptor-Type Systems in the Local Group of Galaxies

Sculptor (Scl) 79 discovered in 1938 by Harlow Shapley, A Stellar System of a New Type Draco 82 A.G. Wilson of the Lowell Observatory in 1955,Sculptor-Type Systems in the Local Group of Galaxies

Sextans 89 Mike Irwin, M.T. Bridgeland, P.S. Bunclark and R.G. McMahon, 1990 A new satellite galaxy of the Milky Way in the of Sextans Ursa Major (UMa) 100 Willman et al. 2005, A New Milky Way Dwarf Galaxy in Ursa Major Carina 103 Cannon, R. D., Hawarden, T. G., & Tritton, S. B., 1977, A new Sculptor-type dwarf elliptical galaxy in Carina

Hercules 140 Belokurov et al. 2007, Cats and Dogs, Hair and A Hero: A Quintet of New Milky Way Companions Fornax 140 discovered in 1938 by Harlow Shapley, described in "Two Stellar Systems of a New Kind", Nature, Vol. 142, p. 715 Canes Venatici II 150 Sakamoto & Hasegawa 2006, Discovery of a Faint Old Stellar System at 150 kpc Leo IV 160 Belokurov et al. 2007, Cats and Dogs, Hair and A Hero: A Quintet of New Milky Way Companions Pisces II 182 Belokurov et al. 2010, Big Fish, Little Fish: Two New Ultra- Faint Satellites of The Milky Way Leo II (Leo B) 208 Robert G. Harrington and Albert George Wilson, 1950, Two New Stellar Systems in Leo Canes Venatici 220 Zucker et al. 2006 A New Milky Way Dwarf Satellite in Canes Venatici Leo I 254 Robert G. Harrington and Albert George Wilson, 1950, Two New Stellar Systems in Leo

Tabla 1.5: Galaxias satélite de la VL que se encuentran a una distancia menor de 250 kpc, The stromlo Milky Way satellites survey program leader: Helmut Jerjen (anu) Tomada de http://www.mso.anu.edu.au/ 6 de abril de 2015.

28 Colisión de dos galaxias de Disco

Estas galaxias se encuentran en una disposición espacial que corresponde a una estructura de disco denotada como (DoS) que se encuentra orientada de forma casi perpendicular al disco de la VL, este disco se encuentra en su gran mayoría compuesto por galaxias enanas esferoidales (dSph), irregulares (dIrr) y ultra tenues (UFDG) (Metz et al. 2009b).

1.3.1 Disco de satélites de la Vía Láctea (DoS)

La forma en que se encuentran dispuestas las galaxias satélites de la VL se parece a una estructura de tipo disco que se encuentra inclinada aproximadamente 88° respecto al plano del disco de la galaxia. En estudios realizados por (Kroupa et al. 2005), sobre 16 galaxias satélite y partiendo de los datos reportados en la tabla 1.5, estos se comparan con una distribución de densidad isotrópica de ley de potencia, ͯ+ con = ͤ ͦ ͦ ͦ . Utilizando la prueba de Kolmogorov-Smirnov (KS), se puede =  +  +  ʚ͌ʛ  ͌ demostrar͌ ǭ͒ que͓ la͔ muestra enana acumulativa es consistente con una distribución de la densidad radial casi isotérmica, con un margen de confianza del 5 por ciento, siendo 1,8

Figura 1.19 : Posición de los 11 satélites más cercanos a la VL, vistos desde un punto localizado en el infinito y en l = 167°,91. El disco de la VL se representa con la línea horizontal con −25 ≤ X/pc ≤ 25 y el centro de coordenadas es el centro galáctico. La línea punteada marca el plano ajustado visto de perfil para N=11 (Kroupa et al., 2005). Galaxias, Grupo Local y galaxias satélites 29

Figura 1.20 : Imagen ilustrativa de la posible forma del plano que se ajusta a la distribución de disco de las galaxias satélite de la VL, Tomada de http://www.newscientist.com/article/

El trabajo realizado por Metz y Kroupa ha sido bastante extenso, trabajo en el cual se incluyen artículos de (Metz et al. 2007, 2009), trabajos en los cuales se mejora el cálculo del disco de satélites de la VL y se incluyen las nuevas galaxias descubiertas, las denominadas ultra tenues, entre los más recientes se encuentra un artículo del año 2014, donde hace un repaso general de los estudios realizados en los últimos años y concluye recalcando la necesidad de partículas de materia oscura exóticas, de este artículo se resalta la figura 1.21 (Ver página siguiente).

En la figura 1.21 se aprecia que las galaxias ultra tenues (UFDG) se ajustan también a una estructura de disco, que coincide con la estructura de disco de las galaxias satélites más luminosas, esto genera mayor relevancia a los estudios de Metz y Kroupa.

Queda una pregunta en el aire ¿Cuál es el origen de las galaxias satélite de la VL?, ¿Por qué presenta la VL un disco de galaxias satélite?

30 Colisión de dos galaxias de Disco

Figura 1.21 : La distribución de todas las galaxias conocidas hoy en día en el Grupo Local que se define por la esfera velocidad cero. Las galaxias fuera de esta esfera de radio alrededor de 1,5 Mpc se alejan del Grupo Local, mientras que dentro de la esfera caen hacia nosotros. En esta se presentan todas las galaxias satélite dentro de unos 250 kpc alrededor de la VL, el disco galáctico de la VL se ve casi de canto la dirección del polo norte galáctico es hacia arriba. Los satélites de la VL con movimientos propios conocidos y velocidades radiales se muestran como círculos de colores. Los satélites rojos se están alejando del observador, los azules se están acercando al observador. (Kroupa et al, 2014).

1.3.2 Hipótesis de formación del disco de galaxias (DoS)

Materia oscura fría

La materia oscura fría o CDM es uno de los tipos de materia oscura propuestos a principios de los años 80 para explicar la formación de estructura cósmica en el modelo del Big Bang.

La teoría de formación de estructura dentro de la cosmóloga de la materia oscura fría hace predicciones sobre sistemas como el de la VL con relación al número y forma de la distribución de satélites. Según esto la VL debería contener alrededor de 500 sub-halos con 8 masas del orden de M ≈ 10 M⊙ dentro de un radio de 500 kpc, situación que no concuerda con lo observado hasta el momento.

Con la teoría CDM la distribución de subestructuras dentro de halos como el de la Vía Láctea debería ser isotrópica, lo cual Kroupa et al. (2005) mostraron que no coincide con las observaciones. Galaxias, Grupo Local y galaxias satélites 31

Figura 1.22 : Isotropía materia oscura fría (Izquierda), Anisotropía Kroupa et al (Derecha).

En trabajos realizados por (Casallas, Cubillos y Casas 2009) y (Cubillos 2013) se hace una excelente recopilación de los diferentes modelos planteados por la comunidad científica para dar explicación a este fenómeno. A continuación se resaltarán los modelos más relevantes presentados en este documento.

Formación por acreción

En la actualidad se han formulado varias hipótesis de formación para el disco de galaxias satélite de la VL entre ellas, la caída o acreción de asociaciones de galaxias más pequeñas al halo de la VL, (Tully et al. 2006) o la atracción de galaxias enanas que incluían las nubes de Magallanes que formaban parte de un gran sistema magallánico (Nichols et al. 2011).

Mateo 1996 analiza la posibilidad de que el halo galáctico se haya formado por la acreción de galaxias enanas esferoidales basándose en la comparación de poblaciones estelares, cúmulos globulares y contenido de materia oscura de las galaxias enanas y del halo, mencionando algunas características fundamentales de las enanas esferoidales de la VL, su población estelar, variabilidad estelar y contenido de materia oscura.

Pre-reionización cosmológica

En el 2005, (Ricotti and Gnedin 2005) consideran las historias de formación de las enanas en el Grupo Local y comparan sus propiedades con las de galaxias simuladas que formaron todas sus estrellas antes de la reionización cosmológica donde proponen que las galaxias enanas del Grupo Local y todas las otras enanas del universo se formaron en tres diferentes trayectorias de evolución. “Fósiles verdaderos” que formaron la mayor parte de sus estrellas en la época de pre-reionización y tuvieron poca formación estelar, “Fósiles contaminados” que comienzan como fósiles verdaderos, pero tienen un episodio importante de formación de estrellas en el que continúan acretando masa y “sobrevivientes” que comenzaron a formar estrellas después de la reionización. 32 Colisión de dos galaxias de Disco

Gas isotérmico

Mashchenko et al. (2005) proponen un modelo de formación para las galaxias esferoidales enanas en el cual se supone que las estrellas se formaron de un gas isotérmico y en equilibrio hidrostático dentro de halos de materia oscura extendidos, encontrando que dicho modelo describe adecuadamente las propiedades de tres galaxias esferoidales satélites (Draco, Sculptor, y Carina).

Remoción por resonancia

D’Onghia et al. (2009) hacen simulaciones de encuentros entre galaxias enanas de disco y galaxias un poco más grandes, donde descubrieron que estos encuentros generan un proceso llamado “remoción por resonancia” que puede transformar dichas galaxias en enanas esferoidales.

Encuentros cercanos - colisión galáctica

Inicialmente Zwicky (1956) menciona la posibilidad de que perturbaciones violentas observadas en galaxias durante encuentros cercanos pueden formar filamentos de material estelar y como a partir de estos escombros de marea se podrían formar sistemas de galaxias enanas, también una posible relación entre la VL y las nubes de Magallanes mediante la existencia de un puente intergaláctico entre ellas.

En Toomre and Toomre (1972) muestran mediante varias simulaciones de encuentros de galaxias de disco como se producen puentes y colas de material estelar debido a las fuerzas de marea que intervienen en la colisión. Concluyen que las orbitas aproximadamente parabólicas son las que generan colas y puentes de marea más definidas, que pueden ser el origen de sistemas galácticos más pequeños, teniendo en cuenta que dichos encuentros sean cercanos y con velocidades relativamente bajas. También se muestra que estos puentes y colas de marea se forman más claramente si las masas de las galaxias en colisión son diferentes.

Okazaki and Taniguchi (2000), adoptando un escenario de interacción galáctica muestra que si solo unas pocas galaxias enanas se forman en cada colisión galáctica sería posible explicar las relaciones entre la morfología y la densidad tanto de las galaxias enanas como de las gigantes, y a partir de eso es posible analizar la relación entre la formación de galaxias enanas con la transformación de galaxias ricas en gas a galaxias de formación temprana como las elípticas, también analizan la formación de galaxias enanas de marea suponiendo que la interacción galáctica ocurre durante el curso de la formación de estructura jerárquica en el universo, para luego compararla con los datos observacionales. Galaxias, Grupo Local y galaxias satélites 33

Pawlowski et al. (2011) reporta formación de material contraorbitante que emerge de la colisión de dos galaxias de disco, utilizando para la escala de dichas galaxias, las características que tenía la VL hace aproximadamente 10 Gyr. Se muestra como a partir de la colisión de dos galaxias usando modelos de mezcla y fly-by en razones de masa 1-1 y 4-1 para la galaxia objetivo y en caída respectivamente, se producen poblaciones de partículas progradas y retrogradas que ocupan cierta región del espacio similar al ocupado por las galaxias enanas satélites de la VL. También encuentra similitudes en la dirección del momento angular orbital, excentricidades, distancias apocentricas, número de partículas pro y retrogradas, apoyando de esta forma el escenario en el que las enanas satélites de nuestra galaxia puedan ser el resultado de una colisión entre la joven VL y otra galaxia, en la cual durante la colisión se formaron escombros de marea en forma de gas que luego se convirtieron en las pequeñas galaxias satélites que hoy posee.

Por las razones expuestas a lo largo de este capítulo, el objetivo que tiene este estudio es realizar una simulación con el software GADGET2 (Springel et al. 2005) con el fin de determinar si la distribución espacial y de número de galaxias satélite de la VL tuvieron su origen a través de una colisión entre dos galaxias de disco. Este trabajo se fundamenta en el modelo en el cual la VL hace miles de millones de años interactuó con una galaxia de disco “Galaxia fantasma” 1 cuya masa se considera mucho menor que la de la VL, esta interacción genero unas colas de marea que se consideran como las posibles causantes de las galaxias satélite de la VL. Después de la colisión se analizan los escombros de marea en busca de grupos de partículas que se comporten de manera similar a las galaxias satélites de la VL tanto en su distribución espacial como en su número.

1 Se decidió darle este nombre a la galaxia que interactúa con la Vía Láctea, en conversación con mi compañero de grupo de investigación y amigo José Benavides, recordando hace unos años cuando trabaje como físico forense, cuando un vehículo atropella a una víctima solo queda evidencia física en la vía de dicho evento, pero nadie sabe en realidad que vehículo fue el que perpetro el hecho se le denomina (carro fantasma), igual como pasa en el caso de la (Galaxia Fantasma) nadie sabe quién es pero pudo dejar evidencia de su existencia.

2. Marco teórico

En el presente capítulo se pretende realizar una revisión a los diferentes modelos y teorías que enmarcan la realización de la simulación, entre estos modelos se encuentran los que dan cuenta de las características de potencial gravitatorio que presentan las diferentes estructuras que constituyen una galaxia espiral, en segunda instancia se presentan los modelos que dan cuenta de la órbita de las dos galaxias interactuantes y en última instancia se hace revisión de los modelos y teorías que hacen referencia a la interacción entre las galaxias.

2.1 Componentes de una galaxia de disco

Las galaxias de disco son estructuras que se encuentran constituidas por tres componentes estelares básicos; bulbo central, disco, halo de estrellas y un halo de materia oscura.

Figura 2.1: Elementos constitutivos de una galaxia de disco. Tomada de http: //atenea.pntic.mec.es/Antares/modulo7/m7_u103_0.html

36 Colis ión de dos galaxias de Disco

2.1.1 Bulbo central

El bulbo central es una estructura de simetría aproximadamente esférica, que se encuentra en el núcleo de la galaxia. Esta estructura se encuentra formada por estrellas jóvenes de (población I) y estrellas viejas (población II) (Binney and Tremaine 2008).

2.1.2 Disco

La estructura de disco es común en las galaxias de tipo espirales, en la VL esta estructura tiene un diámetro aproximado de 32 kpc y un grosor de alrededor de 0.6 kpc, el sol se encuentra a unos 8 kpc del centro de la galaxia. En la VL esta estructura se encuentra compuesta por estrellas (población I) ricas en metales, gas y polvo. Las órbitas de las estrellas que constituyen el disco son cuasi-circulares (Binney and Tremaine 2008).

2.1.3 Halo estelar

El halo es un sistema de estructura de tipo esferoidal, en la VL este halo tiene un diámetro que puede ser por lo menos de 50 kpc, se encuentra formado por estrellas viejas que son pobres en metales (Población II) las cuales se encuentran en dos grupos, el primer grupo corresponde a las que se encuentran en cúmulos globulares y las segundas se encuentran aisladas, estas se denominan estrellas de alta velocidad. (Binney and Tremaine 2008)

2.1.4 Halo de materia oscura

La materia oscura fue propuesta por Fritz Zwicky en 1933, este planteo la posibilidad de su existencia debido a la extraña forma presentada en la curva de rotación galáctica, ya que se esperaba que dicha curva obedeciera a la curva de rotación que proyectarían las leyes de Kepler para los planetas, es decir que a mayor distancia del centro menor velocidad de rotación y a menor distancia del centro mayor velocidad, sin embargo lo que se aprecia es que la velocidad de rotación de la galaxia es en promedio la misma para cualquier distancia radial respecto al centro ver figura 2.2. (Rubin, 1970). Se piensa que existe un halo de materia oscura de simetría aproximadamente esférica en el cual se encuentra embebida la materia bariónica, sus dimensiones son mucho mayores que la de las estructuras de tipo bariónico, en el caso de la VL este puede tener un diámetro aproximado de hasta 500 kpc. La naturaleza de este tipo de materia es desconocida aun sin embargo en la actualidad hay teorías que indican la existencia de pequeños cuerpos Marco Teórico 37

conocidos como WIMPS (Weakly Interacting Massive Particles) y MACHOS (Massive Compact Halo Objects), sin embargo las observaciones en busca de lente gravitacional indican que el número de MACHOS no es suficiente para explicar la masa que se demanda. (Diemand, Moore y Stadel 2005).

Figura 2.2: Curva de rotación observada para la galaxia M33. Tomada de http://www.learner.org/

2.2 Modelos teóricos para el potencial gravitacional de las estructuras (Bulbo, Disco y Halo de materia oscura)

Para cada una de las estructuras descritas anteriormente se requiere un modelo teórico, a continuación se presentan dichos modelos para el bulbo central, el disco y el halo de materia oscura.

2.2.1 Potencial de Hernquist (Bulbo central)

Para el bulbo central de la galaxia de disco se toma un potencial de tipo esférico conocido como el potencial de Hernquist (Hernquist, 1990 y1993).

El potencial de tipo esférico asociado al bulbo central de la galaxia de disco viene dado por: 38 Colis ión de dos galaxias de Disco

∅# = − , (2.1) ͇́+

ͦ ͕ dónde: M: Representa la masa del esferoide. G: La constante de gravitación universal. r: Punto donde se mide el potencial. a: Radio del núcleo asociado al esferoide (Longitud de escala radial).

A partir de este potencial se puede, resolviendo la ecuación de Poisson, encontrar el perfil de densidad del bulbo central así como su distribución de masa acumulada.

= (2.2) # 2 ( + )ͧ ͇͕  ͦ ͕ ͦ

=  2.3 +ͦ ͦ ͇ʚͦʛ ͇ ʚ ʛ ʚͦ ͕ʛ 2.2.2 Perfil de densidad (Disco)

Para la estructura de disco de las galaxias espirales se toma un perfil de densidad volumétrico que en la dirección del radio del disco R en coordenadas cilíndricas es conocido como tipo exponencial ∑ para n=1 para un perfil de Sérsic dado por; (Springel and White, 1999). ʚ͌ʛ

ͯ ∑ = Ɵė (2.4) 2 ͇ ͦ ʚ͌ʛ ͙ ͌ Expresión en la cual se tiene:

∑ͤ = ͦ (2.5) 2͇ ͌ y por su estructura vertical Z se adopta la decisión común de hoja isotérmica de (Spitzer 1942).

ͯ Ɵė ͦ , = ͦ ℎ , 2.6 4 ͇ ͮͤ  ʚ͌ ͮʛ ͤ ͙ ͙ͧ͗ ƴ Ƹ ʚ ʛ ͌ ͮ ͮ dónde: ρ(R,z): Densidad volumétrica de masa del disco

Md: Masa del disco.

Rd: Longitud escala del disco. Marco Teórico 39

z0: Longitud de escala vertical.

2.2.3 Perfil de Hernquist (Halo de materia oscura)

El halo de materia oscura de la galaxia espiral se modela como una estructura esférica que tiene un perfil de densidad del tipo Hernquist (Ecuación 2.2):

= , (2.9) # 2 ( + )ͧ ͇͕ dónde:  ͦ ͕ ͦ ρh = Densidad de masa del halo de materia oscura. M= Masa halo de materia oscura. a = Longitud de escala radial del halo

En esta ecuación es un parámetro que se define para evitar indeterminaciones de tipo matemático, lo cual se puede explicar físicamente como una estructura de núcleo, este ͕ parámetro también se puede asociar con la longitud de escala y el parámetro de concentración de un halo con perfil NFW (Navarro, Frenk and White 1995).

͗ = . 2 ln 1 + − /(1 + ) (2.10) ͕ ͦ ǭ ʞ ʚ ͗ʛ ͗ ͗ ʟ Siendo : = ͦͤͤ (2.11) ͗ ͦ . ͗ ͦ En esta ecuación el radio virial ͦͤͤ viene dado por la expresión dada en (Springel and White, 1999). ͦ  ͦͤͤ = , (2.12) 10 ͪ ( ) ͦ ͂ ͮ Siendo vc la velocidad circular y H(z) es el parámetro de Hubble dependiendo del corrimiento al rojo z.

Para el halo de materia oscura se tiene que la masa virial ͦͤͤ en términos de la velocidad circular es: ͧ  ͇ = (2.13) ͦͤͤ 10 ( ) ͪ Para las ecuaciones (2.10) y (2.11) la͇ velocidad circular local viene dada por la expresión ́͂ ͮ (Springel and White, 1999):

ͦ ͦͤͤ ͦͤͤ = (2.14) ͦͤͤ ͇́ ͪ ͦ 40 Colis ión de dos galaxias de Disco

2.3 Orbita pre-colisión

A continuación se presentan las características fundamentales que competen a la geometría de la colisión en sus condiciones iniciales. Inicialmente se considera una galaxia “Galaxia fantasma” con su estructura de disco dispuesta de forma perpendicular al disco de la VL y que se acerca a la VL a través de una órbita parabólica como se aprecia en la Figura 2.6.

Figura 2.3: Orbita pre-colisión seguida por la “Galaxia fantasma” (esquina superior izquierda) que se aproxima a la VL.

Se escogió esta orbita de tipo parabólico ya que en la literatura se encuentran referencias indicando que las orbitas parabólicas generan colas de marea bastante prominentes (Toomre & Toomre. 1972). La órbita se describe en términos de lo establecido por (Van de Kamp, 1964). Que parte de la expresión (Ecuación 2.15) de tipo polar que define la sección cónica

= , (2.15) 1 + ( ) ͤ ͦ dónde: ͙͗ͣͧ  r es la distancia de la “Galaxia fantasma” a la VL (Tomada la VL como foco de la órbita). p es la distancia de aproximación más cercana (Parámetro de impacto). e es la excentricidad de la órbita. θ es el ángulo del radio vector respecto al eje x (Ver figura 2.6).

Partiendo de la ecuación (2.15) se pueden calcular las componentes de velocidad paralela

- y perpendicular al radio vector W de la siguiente manera:

ͪ ͪ Marco Teórico 41

- = (2.16) ͦ͘ ͪ W = ͨ͘ (2.17) ͘ ͪ ͦ Ahora utilizando la velocidad areal A planteada porͨ͘ (Van de Kamp, 1964) ecuación (2.18) y del parámetro gravitacional µ ecuación (2.19) se tienen - y W

ͦ ͪ ͪ = , (2.18) 2 ͦ ͘ ̻ = ͨ͘+ , (2.19) con:  ́ʚ͇ ͡ʛ µ: Parámetro gravitacional, M: Masa de la VL, m: Masa Galaxia fantasma.

Partiendo de las expresiones (2.18) y (2.19) se tiene que:

- = (2.20)  ͪ ǰ ͙͙ͧ͢ ʚʛ ͤ

W = 1 + ( ) (2.21)  ͪ ǰ ʞ ͙͗ͣͧ  ʟ Estas ecuaciones son las que describen la órbitaͤ pre-colisión entre la “Galaxia fantasma” y la Vía Láctea.

2.3.1 Parámetro de impacto

Se denomina parámetro de impacto p en la colisión de los dos discos galácticos a la distancia que hay entre sus centros cuando se encuentran unidos por una línea vertical al eje y de la colisión. 42 Colis ión de dos galaxias de Disco

p

Figura 2.4: En esta gráfica se muestra un parámetro de impacto de aproximadamente 35 kpc.

2.4 Colisión de galaxias

2.4.1 Sistemas no colisiónales

Para estudiar sistemas en colisión lo primero que se debe suponer es que estos sistemas estelares son no colisiónales. La definición de sistema colisional o no colisional depende del tiempo de relajación, para definir el tiempo de relajación se debe considerar una estrella de radio R con velocidad atravesando una galaxia, a esta estrella le toma un tiempo para hacerlo, mientras cruza el campo de estrellas constituyentes de la -*.. = / ͪ ͨgalaxia ͌ estasͪ modifican lentamente la velocidad de la estrella debido a su interacción gravitacional, para que la velocidad cambie la estrella debe haber cruzado con estrellas del campo, el es el denominado logaritmo de Coulomb - '3 = /8 ͪ (Binney͢ ͈and ͢͠TremaineÛ 2008), este es el tiempo͢͠ Û necesario para que la estrella cambie su velocidad debido a la interacción con estrellas del campo, de lo anterior se ∆ = - '3 tiene que: ͪ ͪ ͢

- '3 = - '3 -*.. (2.22)

En conclusión se puede decir que todasͨ las estrella͢ sͨ en una galaxia se mueven como si el sistema fuera un campo gravitacional suave y no como una distribución de puntos, esto Marco Teórico 43

garantiza que el número de colisiones entre las estrellas que constituyen los sistemas galácticos interactuantes es demasiado bajo (Binney and Tremaine 2008).

2.4.2 Sistemas estelares interactuantes (definiciones)

En los sistemas galácticos cuando dos estructuras se aproximan la una a la otra, su dinámica cambia, esta situación es la denominada colisión de galaxias, los resultados de este evento dependen directamente de las masas de los objetos involucrados así como de su velocidad relativa de interacción.

Mientras dos galaxias están interactuando, la galaxia objetivo es la galaxia sobre la cual se centran los efectos de la colisión, las otras galaxias que se encuentran cayendo a la galaxia objetivo se llaman galaxias satélites. En el caso en que la masa de la galaxia objetivo es mucho mayor que el da la galaxia satélite a esta se le denomina galaxia anfitriona.

Figura 2.5: Imagen donde se muestra en colisiones reales de galaxias, se aprecia las colas de marea generadas por la interacción entre dos galaxias. Tomada de ttp://www.extragalactic.info/

Cuando se presenta la interacción entre dos galaxias, la galaxia anfitriona puede atraer las partículas más cercanas pertenecientes a la galaxia satélite, esto causado principalmente por fuerte interacción gravitatoria entre los objetos. Los fuertes cambios en la fuerza de interacción gravitatoria de los objetos producen colas de marea de estrellas y gas entre las galaxias como se observa en la Figura 2.5. En estas interacciones las partículas que formaban parte de la galaxia satélite se convierten en partículas de la galaxia anfitriona, 44 Colis ión de dos galaxias de Disco

este suceso tiene una gran dependencia con la distancia a la cual se encuentran las partículas de la galaxia satélite respecto al centro de masa de la galaxia anfitriona, esta distancia es la que se define como parámetro de impacto p. El radio rt de marea define la región para la cual ocurre la transferencia de masa entre las dos galaxias, se debe tener presente que este tipo de interacciones son las que determinan la evolución y el estado final de las galaxias.

Existen diferentes métodos para estimar el radio de marea, de forma numérica se dice que el radio de marea se encuentra en el momento en que una de las partículas de la galaxia satélite se desprende y se convierte en partícula de la galaxia anfitriona. También se puede calcular el radio de marea partiendo de las observaciones usando perfiles de densidad y velocidad de las estrellas. Se puede calcular realizando la gráfica del perfil de densidad y mirando el radio finito para el cual la densidad se hace cero (Binney and Tremaine 2008).

2.4.3 Encuentro de galaxias en rotación

La rotación de las galaxias al momento de su interacción o colisión es un factor muy relevante, como muestran (Toomre 1964 y Toomre 1972), se deben considerar dos casos extremos, primero cuando el espín de rotación de la galaxia es anti-paralelo al momentum angular de la órbita, esta situación es la que se denomina encuentro retrogrado y cuando el spin es paralelo al momentum angular se denomina encuentro progrado. Entre los 60s y los 70s Toomre y Toomre realizaron experimentos numéricos sencillos de tipo pionero como se aprecia en las Figuras 2.6 y 2.7. Marco Teórico 45

Figura 2.6: Disco con spin anti-paralelo al momentum angular de la órbita (Toomre 1972).

En la Figura 2.6 se puede apreciar una partícula masiva central que se ve orbitada por cinco anillos de partículas atrapados en su campo gravitatorio, este sistema se mueve a través de una órbita parabólica de forma retrograda en interacción con otro sistema en las mismas condiciones. Como se puede apreciar la interacción no es muy fuerte para los anillos interiores del sistema, sin embargo en la Figura 2.7 donde los sistemas presentan spines paralelos al momentum angular de sus orbitas, se hace evidente la destrucción de la mayoría de anillos que constituyen el sistema.

Esto nos muestra que los encuentros de tipo progrado son más violentos que los encuentros de tipo retrogrado. 46 Colis ión de dos galaxias de Disco

Figura 2.7: Disco con spin paralelo al momentum angular de la órbita (Toomre 1972).

Esto ocurre debido a que las fuerzas de interacción gravitacional de uno de los discos hacen que se transfiera el momento angular a los anillos del otro durante su encuentro. La transferencia de momento angular es menor en los encuentros de tipo retrogrado debido a la inercia rotacional de los anillos (Toomre 1972). En los encuentros progrados la transferencia de energía orbital es máxima sobre las partículas de los anillos del otro sistema.

2.4.4 Encuentros de discos planos (estabilidad)

Para poder hablar del criterio de estabilidad dinámica de un disco estelar definido por (Toomre 1972) se deben definir conceptos como; frecuencia epicíclica y dispersión de velocidades σ. 

En astrofísica, en particular en el estudio de los discos de acreción, la frecuencia epicicloidal es la frecuencia a la que una pequeña masa que es desplazada radialmente Marco Teórico 47

oscilará. Al considerar un disco astrofísico con rotación diferencial f, la frecuencia epicicloidal está dada por:  2Ω ͦ = ͦ (2.23) ͘  ʚ͌ ʛ Donde R es la coordenada radial, esta cantidad͌ ͌͘ puede usarse para examinar los "límites" de un disco de acreción, cuando ͦ es negativa entonces una pequeña perturbación en la órbita (supuesta circular) de la  masa se volverá inestable. Se tiene que para un disco kepleriano = Ω (Pringle and King 2007).  En cuanto a la dispersión de velocidades σ se define como la dispersión estadística de las velocidades sobre la velocidad media de un grupo de objetos, estos pueden ser una galaxia, un cúmulo de galaxias, o un supercúmulo entre otros. Midiendo las velocidades radiales de sus miembros, la dispersión de la velocidad de un clúster puede ser estimada y utilizada para obtener la masa del cúmulo a partir del teorema virial (Collins Dictionary of Astronomy 2000).

El criterio de estabilidad establecido por Toomre (1972) se lleva a cabo realizando una

comparación entre dispersión de velocidades radial σR y la frecuencia epicíclica. Las perturbaciones radiales sobre el disco producen fuerzas gravitatorias extras de perturbación interna y externa las cuales se denominan correspondientemente compresión y presión. Estas perturbaciones producen un aumento en la dispersión de velocidades esto generando que la aproximación epicíclica se rompa, lo cual genera que el disco pierda su estabilidad y pueda ser destruido. Este criterio se expresa analíticamente de la siguiente manera  = , (2.24) 3.36 ∑   ͋ ́ expresión en la que se tiene que ∑ es la densidad superficial del disco. Ahora, en la colisión de discos galácticos si la dispersión de velocidades se incrementa se dice que el disco se calentó, en un disco caliente si Q > 1 el disco será estable y resistirá las perturbaciones radiales, si esto no es así el disco será inestable y perderá su soporte rotacional. Los eventos que ocurren dentro de la interacción de dos galaxias de disco sobre cada una de las estructuras dependen de sus velocidades, parámetro de impacto, masas y su estabilidad dinámica al momento del encuentro (Quiroga 2013). 48 Colis ión de dos galaxias de Disco

2.4.5 Fricción dinámica

En astrofísica, la fricción dinámica, a veces llamada arrastre gravitacional, es la pérdida de momento y de energía cinética de los cuerpos en movimiento a través de las interacciones gravitacionales con la materia circundante en el espacio (Shandrasekhar 1943). La fricción dinámica aparece en todos los procesos de interacción de tipo galáctico.

Un caso especial es comúnmente utilizado donde hay una densidad uniforme en el campo de materia, con partículas de materia significativamente más ligeras ( ma << M) que la principal partícula bajo consideración y con una distribución de Maxwell para la velocidad de las partículas de materia. Bajo estas condiciones las partículas se ven atraídas por una fuerza de tipo gravitatoria grande, esto genera una sobre densidad de partículas en la trayectoria del objeto, así se genera una fuerza contraria a la trayectoria de desplazamiento del objeto, de esta forma el objeto pierde energía cinética y su velocidad decrece. En este caso, la fuerza de fricción dinámica es la siguiente (Shandrasekhar 1943).

dV Τ‘ V ΁ ͦ ͦ ͦ ΁ = −16π G MmΏlnΛ vΏ f(vΏ)dv Ώ ͧ (2.25) dt ͤ v΁ ʨǹ ʩ

Donde es el logaritmo de Coulomb y ( ·) es la función de distribución de velocidades de las partículas. ͢͠Û ¼ Ì Como se aprecia en la ecuación (2.25) la fricción dinámica tiene directa dependencia con la masa de los objetos circundantes, los cuales derivan de una distribución o perfil de densidad, lo cual indica que la fricción dinámica podría ser un indicador del perfil de densidad de un halo de materia oscura de una galaxia anfitriona. Suponiendo que el perfil de densidad del halo de materia oscura de una galaxia anfitriona fuera un perfil tipo Hernquist se puede demostrar que el tiempo de fricción dinámica es:

1 1 3 ͧ ͧ ͦ ͦ !-$ = − − + − , (2.26) 4 3 $ 2͕ $ ͨ Ƭ Ƴͦ ͦ Ʒ Ƴͦ ͦ Ʒư ͇͕́͢͠ Û donde a es el factor de escala del halo, r i es la posición del objeto en el tiempo t i y si r=0 se tiene que t= tfricc .

2.4.6 Algunas colisiones reales y simulación

A continuación se mostrará que las masas y formas de los sistemas galácticos que interactúan se pueden reproducir de una forma bastante aproximada por los experimentos Marco Teórico 49

de tipo numérico denominados como simulaciones, estas simulaciones permiten comprender toda la dinámica de formación y evolución de las galaxias.

Figura 2.8: Este par de galaxias, NGC 4676, también conocido como "Los Ratones" por sus colas de estrellas y gas, han colisionado y eventualmente se fusionarán en una sola galaxia. Corrientes de material han sido tirados fuera de las galaxias por la fuerza de gravedad, lo que provocó nuevo nacimiento de estrellas. Tomada de http://hubblesite.org.

Figura 2.9: En esta figura se aprecia que las simulaciones numéricas de N-cuerpos pueden dar cuenta de las interacciones que se encuentran en las observaciones.

Como se aprecia en la figura 2.12 es posible recrear con las leyes físicas conocidas en la tierra eventos de tipo galáctico, por lo menos en primera instancia en lo que compete a su morfología.

3. Condiciones iniciales de simulación

Teniendo presente que la Vía Láctea y las galaxias satélite simuladas en este estudio se encuentran constituidas por una componente estelar de bulbo y disco embebidas en un halo de materia oscura, se describirán las características principales del software usado para construir los objetos simulados, además, se mostrarán los parámetros empleados para construir las diferentes componentes de las galaxias y por último se mostrarán las condiciones iniciales de las orbitas elegidas para llevar cabo la simulación.

3.1 Software Zeno

Este software fue desarrollado por Joshua Barnes (Barnes 2011), se usa para generar condiciones iniciales de simulación, que permite generar estructuras de tipo astrofísico con N-cuerpos y diferentes tipos de potencial gravitacional, entre estos potenciales se encuentran los mencionados en el capítulo anterior: perfil de Hernquist (halo y bulbo) y perfil exponencial (disco).

El paquete de software ZENO integra códigos de N-cuerpos y simulación SPH (Smooth particle Hydrodynamics) con una gran variedad de programas para generar las condiciones iniciales y analizar las simulaciones numéricas. Escrito en C, el sistema ZENO es portable entre plataformas Mac, Linux y Unix. Es de uso activo en el Instituto de Astronomía (IFA), NRAO, STScI, y en otros lugares. El código fuente está disponible bajo la licencia GPL de GNU.

Desde el punto de vista del usuario, ZENO consiste en una colección de programas que se invocan y controlados directamente desde la línea de comandos de UNIX. Sin embargo, el sistema utiliza una variedad de técnicas para proporcionar un entorno de trabajo flexible y potente (Barnes 2011).

Zeno puede realizar una amplia gama de tareas de simulación y análisis, como ejemplo: Rutinas de generación de instantáneas crean distribuciones de partículas con diversas propiedades. Los sistemas con perfiles de densidad especificados por el usuario se pueden poner a interactuar en sistemas no colisiónales o en sistemas gaseosos; múltiples componentes esféricas y de disco se pueden establecer en equilibrio mutuo. Los códigos de simulación incluyen programas combinados N-cuerpo / SPH N-cuerpo, los códigos SPH

52 Colisión de dos galaxias de Disco

ofrecen una amplia gama de opciones para la física de gas, incluyendo modelos isotérmico, adiabático, y radiantes (Barnes 2011).

3.2 Parámetros de la Vía Láctea

Las condiciones iniciales de la VL se realizaron de forma independiente para sus tres componentes Halo, Disco y Bulbo, estas condiciones iniciales son tomadas de la literatura y hacen referencia a las observaciones. La VL fue creada para simulaciones de baja resolución (pocas partículas) y simulaciones de alta resolución (muchas partículas). Las simulaciones con pocas partículas permiten determinar los mejores parámetros iniciales para las simulaciones de alta resolución.

Parámetro Descripción Vía Láctea Unidades

12 Mh Masa halo de materia oscura 1,55 x 10 M⊙

Rvir o R 200 Radio virial del halo 244, 8 Kpc

aH Longitud de escala del halo 37 Kpc

10 Md Masa del disco 5,0 x 10 M⊙

RS Longitud de escala radial del disco 4 Kpc

RZ Longitud de escala vertical del disco 0, 3 Kpc

10 Mb Masa del bulbo 1, 8 x 10 M⊙

RSB Longitud de escala bulbo 0, 6 Kpc C Parámetro de concentración 12 -

Tabla 3.1: Parámetros galácticos para la Vía Láctea usados en la simulación (Frenk and White 2014), (Pawlowski 2011) y (Binney and Tremaine 2008).

Condiciones Iniciales 53

3.2.1 Condiciones iniciales baja resolución Vía Láctea

Figura 3.1: VL simulada y virializada 2Gyr con 75000 partículas, verde (bulbo) y rojo (disco).

Figura 3.2: Energía de la VL; energía cinética (verde), energía potencial (rojo) y energía total (azul).

54 Colisión de dos galaxias de Disco

3.2.2 Condiciones iniciales alta resolución Vía Láctea

Figura 3.3: Vía Láctea simulada y virializada 2Gyr con 150000 partículas, verde (disco) y rojo (bulbo).

Figura 3.4: Energía de la VL; energía cinética (verde), energía potencial (rojo) y energía total (azul).

Condiciones Iniciales 55

3.2.3 Masa del halo de materia oscura

En la literatura científica se encuentran muchos valores para los diferentes parámetros de la VL, pero quizás los que son más difíciles de encontrar por la falta de consenso son los parámetros que competen a las dimensiones del halo de materia oscura en el que se encuentran embebidos el disco y el bulbo de la galaxia, en especial con relación a su masa. Por las razones expuestas anteriormente es que en este estudio se realizó un promedio con los valores de masa del halo de materia oscura encontrados por las distintas técnicas reportadas en (Frenk and White 2014) ver Figura 3.5 con lo cual se determinó un valor 12 de M=1,55 x 10 M⊙ valor que se utilizó para realizar las distintas simulaciones.

Figura 3.5: En esta figura se aprecian los diferentes valores de masa estimados para el halo de materia oscura de la VL, las técnicas utilizadas para estos estimativos son diversas, entre estas se encuentran; fracción de bariones cósmica, coincidencia de abundancia y la dinámica de satélites distantes entre otras. Se resaltan los valores 12 encontrados por los modelos de materia oscura fría (M=1,2 x 10 M⊙) y materia oscura 12 caliente (2.6 x 10 M⊙) (Frenk and White 2014).

56 Colisión de dos galaxias de Disco

3.3 Galaxias satélites

En este trabajo se elaboró un modelo para construcción de galaxias satélite, este consiste en crear galaxias de menores o mayores dimensiones que una galaxia dada. Para esto se suponen los nuevos objetos como múltiplos de la galaxia conocida. Esta metodología se podría denominar como un redimensionamiento de la VL. A continuación se explica en que consiste esta técnica.

3.3.1 Redimensionamiento del Halo de materia oscura

En primera instancia se realiza el redimensionamiento del halo de materia oscura, en el caso del halo solo se redimensiona su masa en términos de su velocidad circular v C, para determinar las características del halo de la galaxia satélite se utilizan las siguientes expresiones (Springel and White, 1999).

ͧ  = (3.1) ͦͤͤ 10 ( ) ͪ ͇ Se procede a calcular una velocidad circular pará͂ el haloͮ de la galaxia enana que cumpla con la condición:

= (3.2) ͇ ͢ ͇ Donde n es el porcentaje de relación de masas entre las galaxias, M L la masa de la galaxia

anfitriona (VL) y M S la masa de la galaxia satélite (Galaxia fantasma).

3.3.2 Redimensionamiento del disco

En segunda instancia se redimensiona el disco galáctico, para este caso se deben redimensionar dos longitudes fundamentales para la estructura, primero la longitud de escala radial R S y segundo la longitud de escala vertical RZ utilizando las siguientes ecuaciones (Springel and White, 1999). ∑ͤ = ͦ (3.3) 2 ͇ Teniendo presente que la relación de masas k ͌ entre la masa del disco de la galaxia anfitriona M dL y la masa del disco de la galaxia satélite M dS es:

= (3.4)  ͇ ͟ ͇ Ahora asumiendo que ∑ͤ es igual para ambas galaxias se tiene que: Condiciones Iniciales 57

 =  (3.5) ͌ ͌ Donde R SL es la longitud de escala radial√͟ de la galaxia anfitriona y R SS es la longitud de escala radial de la galaxia satélite.

A continuación se presenta la forma en que se calculó la longitud de escala vertical para el disco de la galaxia satélite, para su determinación se toma como referencia el cociente

entre la longitud de escala radial de la galaxia anfitriona RSL y su longitud de escala

vertical z oL .

 = (3.6) ͌ͤ  ͮ Con el valor de τ y conociendo ya la longitud de escala radial del disco de la galaxia satélite R SS se procede a calcular z 0 para la misma como sigue:

 = ͤ (3.7) ͌ ͔  Donde z 0S es la longitud de escala vertical de la galaxia satélite.

3.3.3 Redimensionamiento del bulbo

En última instancia se redimensiona el bulbo, para este caso se utilizan las siguientes expresiones de densidad de masa volumétrica # del bulbo (Hernquist, 1990).

=  (3.8) # 2 ( + )ͧ ͇͕ Teniendo presente que la relación de masa entre el bulbo de la galaxia anfitriona y el ͦ ͕ ͦ  de la galaxia satélite se define como:  ͇

͇ =  (3.9) ͇ ì ͇ Asumiendo igualdad en la densidad la densidad de masa del halo de la galaxia anfitriona es igual a la densidad de masa del halo de la galaxia satélite, se pueden calcular la longitud de escala del bulbo así:

w =  (3.10) ͕ ͕ ǭ 58 Colisión de dos galaxias de Disco

Donde a L es la longitud de escala del bulbo de la galaxia anfitriona y a S es la longitud de escala del bulbo de la galaxia satélite.

3.3.4 Corroborando el modelo

Para contrastar este modelo con otros estudios realizados, se comparó con el realizado por (Pawlowski, Kroupa y de Boer, 2011) y titulado “Making Counter-Orbiting Tidal Debris The Origin of the Milky Way Disc of Satellites” en este artículo se realizan simulaciones de N-cuerpos para la colisión de dos galaxias de disco que pudieron dar origen al disco de galaxias satélites de la VL. Este reporta los parámetros correspondientes a los sistemas interactuantes en la siguiente tabla:

Tabla 3.2: Los parámetros de los modelos de galaxias objetivo y proyectil de disco (Pawlowski, 2011).

Se construye la galaxia denominada por Pawlowski (target) con las mismas dimensiones y razones de masa establecidas en el artículo, después usando la metodología planteada en los puntos anteriores, se construye la galaxia denominada (infalling) y se comparan los valores obtenidos en el artículo para la misma y los obtenidos por la metodología usada en este trabajo. En la tabla 3.3 se muestran los parámetros que se compararon:

Condiciones Iniciales 59

Parámetro Descripción Infall Infall Dif % Pawlowsky Own 10 10 Mh Masa halo de materia oscura 4,0 x 10 M⊙ 3,9 x 10 M⊙ 0,95 aH Longitud de escala del halo 6,9 kpc 6,9 kpc 0,0 Radio máximo halo 100 kpc 114,6 kpc 14.6 RmáxH 9 9 Md Masa del disco 4,0 x 10 M⊙ 3,9 x 10 M⊙ 0,25

RS Longitud de escala radial del 1,15 kpc 1,12 kpc 2,61 disco

Rmáxd Radio máximo del disco 5,75 5, 62 2.26

RZ Longitud de escala vertical del 0,20 0,175 kpc 12,5 disco

Tabla 3.3: Tabla comparativa entre los datos usados por Pawlowsky y los datos obtenidos por el modelo usado por este estudio, se nota que la diferencia porcentual media entre los datos es de 4.74% lo cual es un buen indicativo para la utilización de este modelo.

3.3.5 Parámetros de las galaxias satélites

Para definir los valores de masa de las galaxias satélite que van a entrar en interacción con la VL se tuvo en cuenta el estudio realizado por (Quiroga 2013) y titulado “Evolution of Disk Galaxies in ΛCDM”, en este estudio se concluye que satélites con relaciones de masa menores o iguales al 16% al instante de una colisión permiten que la galaxia anfitriona sobreviva, ya que la transferencia de momentum lineal y momentum angular es baja en estas interacciones.

Conociendo estas restricciones, en este trabajo se utilizaron galaxias satélites con porcentajes de masa de 4%, 8%, 12% y 16 % de la galaxia anfitriona. Los valores de masa, longitud de escala, etc. Se pueden encontrar en la tabla 3.4.

60 Colisión de dos galaxias de Disco

Parámetro Descripción 4% 8% 12% 16% Unidades

10 11 11 11 Mh Masa halo de 6,20 x 10 1,24 x 10 1,94 x 10 2,48 x 10 M⊙ materia oscura

Rvir o R 200 Radio virial del 83,7 105,5 122,4 132,9 kpc halo

aH Longitud de 8,0 10,1 11,7 12,7 kpc escala del halo 9 9 9 9 Md Masa del disco 2 x 10 4 x 10 6,25 x 10 8 x 10 M⊙

RS Longitud de 0,8 1,13 1,41 1,60 kpc escala radial del disco

RZ Longitud de 0,06 0,085 0,106 0,120 kpc escala vertical del disco 8 9 9 9 Mb Masa del bulbo 7,2 x 10 1,44 x 10 2,25 x 10 2,88 x 10 M⊙

RSB Longitud de 0, 205 0,26 0,30 0,33 kpc escala del bulbo C Parámetro de 21, 8 21, 8 21, 8 21, 8 - concentración

Tabla 3.4: Se muestran los parámetros de las diferentes galaxias satélites utilizadas en este estudio, calculados con la metodología planteada en los puntos 3.3.1, 3.3.2 y 3.3.3.

Para este estudio se realizaron dos tipos de simulaciones, una de ellas denominada de baja resolución en la cual se presentan un número de partículas inferior a 1 x 10 5 partículas por cuerpo y otras las denominadas de alta resolución que presentan 5 x 10 5 partículas. A continuación se presentan las características morfológicas y energéticas de las galaxias utilizadas para las simulaciones de baja resolución.

Condiciones Iniciales 61

3.3.5.1 Condiciones iniciales de las galaxias satélites de baja resolución

Galaxia satélite 4%

Figura 3.6: Características galaxia satélite (4%) virializada 2Gyr con 75000 partículas, verde (bulbo) y rojo (disco)

Figura 3.7: Energía de la galaxia satélite (4%), energía cinética (verde), energía potencial (rojo) y energía total (azul).

62 Colisión de dos galaxias de Disco

En las Figuras 3.3 y 3.4 se aprecian las características morfológicas y dinámicas de la galaxia que tiene un 4% de la masa de la VL.

Galaxia satélite 8%

Figura 3.8: Características galaxia satélite (8%) virializada 2Gyr con 75000 partículas, verde (bulbo) y rojo (disco)

Figura 3.9: Energía de la galaxia satélite (8%), energía cinética (verde), energía potencial (rojo) y energía total (azul).

En las Figuras 3.5 y 3.6 se aprecian las características morfológicas y dinámicas de la galaxia que tiene un 8% de la masa de la VL. Condiciones Iniciales 63

Galaxia satélite 12%

Figura 3.10: Características galaxia satélite (12%) virializada 2Gyr con 75000 partículas, verde (bulbo) y rojo (disco)

Figura 3.11: Energía de la galaxia satélite (12%), energía cinética (verde), energía potencial (rojo) y energía total (azul).

En las Figuras 3.7 y 3.8 se aprecian las características morfológicas y dinámicas de la galaxia que tiene un 8% de la masa de la VL.

64 Colisión de dos galaxias de Disco

Galaxia satélite 16%

Figura 3.12: Características galaxia satélite (16%) virializada 2Gyr con 75000 partículas, verde (bulbo) y rojo (disco)

Figura 3.13: Energía de la galaxia satélite (16%), energía cinética (verde), energía potencial (rojo) y energía total (azul).

En las Figuras 3.9 y 3.10 se aprecian las características morfológicas y dinámicas de la galaxia que tiene un 16% de la masa de la VL.

Como se puede apreciar la energía mecánica total es constante para cada sistema, además el sistema alcanza el equilibrio virial rápidamente. Condiciones Iniciales 65

3.3.5.2 Condiciones iniciales de las galaxias satélites de alta resolución

Las características morfológicas y dinámicas de las galaxias de alta resolución son las mismas que las de las galaxias de baja resolución, sin embargo las imágenes de más resolución permiten apreciar mejor las características visuales de las galaxias, además para la simulación generan una aproximación mejor especialmente con lo que tiene que ver con el número de partículas y las distancias entre las mismas las cuales permitirán un mejor análisis en la búsqueda de grupos en los remanentes de la colisión.

Galaxia satélite 4%

Figura 3.14: Características galaxia satélite (4%) virializada 2Gyr con 350300 partículas, verde (bulbo) y rojo (disco)

66 Colisión de dos galaxias de Disco

Figura 3.15: Energía de la galaxia satélite (4%), energía cinética (verde), energía potencial (rojo) y energía total (azul).

En las Figuras 3.11 y 3.12 se aprecian las características morfológicas y dinámicas de la galaxia que tiene un 4% de la masa de la VL.

Galaxia satélite 8%

Figura 3.16: Características galaxia satélite (8%) virializada 2Gyr con 350300 partículas, verde (bulbo) y rojo (disco) Condiciones Iniciales 67

Figura 3.17: Energía de la galaxia satélite (8%), energía cinética (verde), energía potencial (rojo) y energía total (azul).

En las Figuras 3.13 y 3.14 se aprecian las características morfológicas y dinámicas de la galaxia que tiene un 8% de la masa de la VL.

Galaxia satélite 12%

Figura 3.18: Características galaxia satélite (12%) virializada 2Gyr con 350300 partículas, verde (bulbo) y rojo (disco). 68 Colisión de dos galaxias de Disco

Figura 3.19: Energía de la galaxia satélite (12%), energía cinética (verde), energía potencial (rojo) y energía total (azul).

En las Figuras 3.15 y 3.16 se aprecian las características morfológicas y dinámicas de la galaxia que tiene un 12% de la masa de la VL.

Galaxia satélite 16%

Figura 3.20: Características galaxia satélite (16%) virializada 2Gyr con 350300 partículas, verde (bulbo) y rojo (disco) Condiciones Iniciales 69

Figura 3.21: Energía de la galaxia satélite (16%), energía cinética (verde), energía potencial (rojo) y energía total (azul).

En las Figuras 3.17 y 3.18 se aprecian las características morfológicas y dinámicas de la galaxia que tiene un 16% de la masa de la VL.

3.4 Parámetros orbitales

A continuación se presentan los parámetros para la órbita parabólica elegida para las diferentes galaxias satélites que interactuaran con la VL. En este punto vale la pena mencionar que se realizan primero unas simulaciones con pocas partículas (150000 partículas entre los dos objetos) para determinar los parámetros orbitales óptimos y así generar las colas de marea más grandes, en las cuales se desprenda la mayor cantidad de materia posible de la galaxia satélite.

3.4.1 Simulaciones de baja resolución

Estas simulaciones tienen como objetivo determinar cuál es el valor óptimo de parámetro de impacto para la colisión entre las galaxias satélite y la VL, los valores usados para este objetivo se muestran en la Tabla 3.5.

70 Colisión de dos galaxias de Disco

Parámetro Descr ipción 4% 8% 12% 16% Unidades P Parámetro de 55 50 45 40 Kpc impacto µ Parámetro 2,23 x 10 32 2,32 x 10 32 2,41 x 10 32 2,49 x 10 32 kg/Nm 2 Gravitatorio Θ Angulo medido 130 130 130 130 ° respecto a x X Posición en x 0 0 0 0 Kpc Y Posición en y -98,97 -89,97 -80,98 -71,98 Kpc Z Posición en z 117,95 107,23 96,50 85,78 Kpc

vx Velocidad en x 0 0 0 0 km/s

vy Velocidad en y 277,95 297,07 319,60 344,23 km/s

vz Velocidad en z -129,61 -138,53 -149,03 -160,51 km/s

Tabla 3.5: Parámetros de las diferentes galaxias satélite utilizadas en este estudio para las orbitas de tipo parabólico en simulaciones de baja resolución.

3.4.1.1 Condiciones orbitales iniciales para simulaciones de baja resolución

A continuación se presentan figuras que muestran las condiciones orbitales iniciales que se pueden apreciar en la tabla 3.5. En estas imágenes se pueden apreciar las características geométricas de la órbita incluyendo el parámetro de impacto (p).

Condiciones Iniciales 71

Galaxia satélite 4%

Figura 3.22: Condiciones orbitales iniciales galaxia satélite 4%, orbita parabólica y parámetro de impacto de 55 kpc.

Galaxia satélite 8%

Figura 3.23: Condiciones orbitales iniciales galaxia satélite 8%, orbita parabólica y parámetro de impacto de 50 kpc.

72 Colisión de dos galaxias de Disco

Galaxia satélite 12%

Figura 3.24: Condiciones orbitales iniciales galaxia satélite 12%, orbita parabólica y parámetro de impacto de 45 kpc.

Galaxia satélite 16%

Figura 3.25: Condiciones orbitales iniciales galaxia satélite 16%, orbita parabólica y parámetro de impacto de 40 kpc.

Condiciones Iniciales 73

Galaxia satélite 16%

Figura 3.26: Condiciones orbitales iniciales galaxia satélite 16%, orbita parabólica y parámetro de impacto de 35 kpc.

Como se puede apreciar en la tabla 3.5 las galaxias satélites presentan un parámetro de impacto diferente, para las distintas masas, esto se llevó a cabo teniendo presente que para parámetros de impacto iguales con distintas masas los cambios porcentuales en el parámetro gravitacional y en las velocidades orbitales son del orden del 10% máximo para el parámetro gravitacional y del 5% máximo para las velocidades orbitales, lo cual permite generalizar los resultados para los diferentes parámetros de impacto, es decir que si por ejemplo se lanza la galaxia de masa 8% con un parámetro de impacto de 45 kpc se obtendrán resultados muy similares a los obtenidos si se lanzara la galaxia de masa 16% con el mismo parámetro de impacto.

3.4.1.2 Condiciones orbitales iniciales de simulaciones de alta resolución

La tabla 3.5 muestra los parámetros utilizados en las simulaciones de alta resolución para las cuales se usó el parámetro de impacto obtenido de las simulaciones de baja resolución que se considera como óptimo para generar las colas de marea más prominentes después de la colisión. Como se puede apreciar en esta tabla el parámetro de impacto al que se hace referencia es un parámetro de 45 kpc.

74 Colisión de dos galaxias de Disco

Parámetro Descripción 4% 8% 12% 16% Unidades P Parámetro de 45 45 45 45 kpc impacto µ Parámetro 2,23 x 10 32 2,32 x 10 32 2,41 x 10 32 2,49 x 10 32 kg/Nm 2 Gravitatorio Θ Angulo medido 130 130 130 130 ° respecto a x X Posición en x 0 0 0 0 kpc Y Posición en y -80, 98 -80, 98 -80,98 -80, 98 kpc Z Posición en z 96, 50 96, 50 96 ,50 96, 50 kpc

vx Velocidad en x 0 0 0 0 km/s

vy Velocidad en y 307, 28 313, 14 319,60 324,54 km/s

vz Velocidad en z -143, 29 -146, 02 -149,03 -151, 33 km/s Tabla 3.6: Tabla en la que se muestran los parámetros de las diferentes galaxias satélite utilizadas en este estudio para las orbitas de tipo parabólico planteadas.

Galaxias satélites 4%, 8%, 12% y 16%

Figura 3.27: Condiciones orbitales iniciales para todas las galaxias satélite en las simulaciones de alta resolución, orbita parabólica y parámetro de impacto de 45 kpc.

Condiciones Iniciales 75

A partir de los datos expuestos a lo largo de este capítulo se llevaron a cabo las simulaciones, que se presentan en el capítulo cuatro, las cuales se presentan como el objetivo primario de este estudio, realizadas estas simulaciones se llevarán a cabo los análisis correspondientes a los resultados que arrojen, haciendo énfasis en la distribución espacial y cantidad de grupos que se forman en los remanentes de la colisión.

4. Simulaciones y Análisis de Resultados

En el presente capítulo se muestran los resultados obtenidos en las simulaciones realizadas con las relaciones de masa de 4, 8, 12 y 16% de la masa de la VL y además se realizan los análisis correspondientes a la disposición espacial de los restos de la colisión, así como en el número de grupos que se formaron dentro de los remanentes, los cuales indican la cantidad de galaxias satélite que se pudieron haber formado producto de la colisión, estos resultados y análisis serán contrastados con las observaciones reportadas en la literatura.

4.1 Gadget2

Gadget2 es un software de simulación numérica de N-cuerpos, recibe su nombre por las siglas en ingles (GAlaxies with and Gas intEracT), desarrollado y publicado por (Springel, Yoshida & White 2001a), software que ha permitido desarrollar simulaciones cosmológicas hasta con 10 10 partículas. Con Gadget2 se pueden desarrollar simulaciones enfocadas en eventos físicos como colisión o interacción de galaxias, el colapso adiabático de esferas auto-gravitantes de gas, formación cosmológica de clústeres de galaxias y formación de estructura a gran escala entre otras.

Gadget2 permite realizar simulaciones en paralelo sobre los procesadores o diferentes computadores de un sistema de cómputo, el software fundamenta su funcionamiento en un código de árbol en el cual las interacciones de gravitatorias se calculan con una expansión multipolar jerárquica, la dinámica del gas se calcula usando hidrodinámica de partículas suavizadas (SPH). La fuerza gravitatoria a corto alcance se calcula usando una expansión multipolar jerárquica con algoritmo de árbol y a largo alcance se usan técnicas de Fourier (TreePM o particle-mesh).

4.2 Resultados simulaciones de baja resolución

Se realizó una comparación entre las diferentes relaciones de masa (4%, 8%, 12% y 16%) respecto a los parámetros correspondientes a la órbita y se encontró que para todas las galaxias el cambio porcentual máximo en el parámetro gravitacional es del 10%, lo cual

78 Colisión de dos galaxias de Disco produce un cambio porcentual entre las velocidades de los objetos que no supera el 5%, por lo tanto se procederá a realizar las simulaciones para todas las relaciones de masa con diferentes parámetros de impacto, asumiendo que todas en general se comportaran de forma similar para cada parámetro.

A continuación se presentan los resultados de las simulaciones de baja resolución (ver tabla 3.5) recordando que estas simulaciones permiten determinar cual es el parametro de impacto óptimo para realizar la colisión entre la galaxia satélite y la galaxia anfitriona en simulaciones de alta resolución.

4.2.1 Simulaciones de baja resolución galaxia satélite 4%

A continuación se presentan los resultados obtenidos para la simulación de colisión entre la masa del 4% de la VL y la galaxia anfitriona, los datos empleados en esta simulación son los mostrados en la tabla 3.5.

Como se puede apreciar en las figuras 4.1-4.3 El número de partículas que quedan orbitando la VL en este caso es alrededor de 49 que sumadas sus masas alcanzan a tener alrededor del 30% de masa de una galaxia enana satélite de la VL. Se debe tener presente que por la distribución espacial de las partículas, estas no formarían un grupo.

Figura 4.1: Posición final en el plano X-Z de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 4% y las azules (disco) y morado (bulbo) de la VL, Se debe tener presente que por la distribución espacial de las partículas, estas no formarían un grupo.

Simulaciones y análisis de Resultados 79

Figura 4.2: Posición final en el plano X-Y de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 4% y las azules (disco) y moradas (bulbo) de la VL.

Figura 4.3: Posición final en el plano Y-Z de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 4% y las azules (disco) y moradas (bulbo) de la VL. 80 Colisión de dos galaxias de Disco

4.2.2 Simulaciones de baja resolución galaxia satélite 8%

A continuación se presentan los resultados obtenidos para la simulación de colisión entre la masa del 8% de la VL y la galaxia anfitriona, los datos empleados en esta simulación son los mostrados en la tabla 3.5.

El número de partículas que quedan orbitando la VL en este caso son alrededor de 94 que sumadas sus masas alcanzan a tener alrededor de 72% de masa de una galaxia enana satélite de la VL. Se debe tener presente que por la distribución espacial de las partículas estas no formarían un grupo ver figuras 4.4-4.6.

Figura 4.4: Posición final en el plano X-Z de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 8% y las azules (disco) y morado (bulbo) de la VL, Se debe tener presente que por la distribución espacial de las partículas, estas no formarían un grupo.

Simulaciones y análisis de Resultados 81

Figura 4.5: Posición final en el plano X-Y de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 8% y las azules (disco) y moradas (bulbo) de la VL.

Figura 4.6: Posición final en el plano Y-Z de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 8% y las azules (disco) y moradas (bulbo) de la VL. 82 Colisión de dos galaxias de Disco

4.2.3 Simulaciones de baja resolución galaxia satélite 12%

A continuación se presentan los resultados obtenidos para la simulación de colisión entre la masa del 12% de la VL y la galaxia anfitriona, los datos empleados en esta simulación son los mostrados en la tabla 3.5.

En esta simulación se encuentra una cola de marea que crece hasta los 500 kpc aproximadamente y que se presenta como una estructura de disco perpendicular al disco de la VL, esto se puede apreciar de manera más clara en la figura 4.7 en la que se ve la posición final de los objetos en el plano X-Z, sin embargo se puede notar que la mayor cantidad de materia de la galaxia satélite queda concentrada muy cerca de la VL, situación que no se aprecia en las observaciones. Bajo estas condiciones es factible hacer una búsqueda de grupos de partículas que se puedan comportar como galaxias enanas satélite de la VL.

Figura 4.7: Vía Láctea de color azul disco y morado bulbo, partículas galaxia satélite rojo disco y verde bulbo. En esta imagen se aprecia una cola de marea que crece hasta los 500 kpc aproximadamente.

Simulaciones y análisis de Resultados 83

Figura 4.8: Posición final en el plano X-Y de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 12% y las azules (disco) y moradas (bulbo) de la VL

. Figura 4.9: Posición final en el plano Y-Z de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 12% y las azules (disco) y moradas (bulbo) de la VL. 84 Colisión de dos galaxias de Disco

4.2.4 Simulaciones de baja resolución galaxia satélite 16%

A continuación se presentan los resultados obtenidos para la simulación de colisión entre la masa del 16% de la VL y la galaxia anfitriona, los datos empleados en esta simulación son los mostrados en la tabla 3.5.

En esta imagen se aprecia una cola de marea que crece hasta los 200 kpc aproximadamente Y que se presenta como una estructura de disco perpendicular al disco de la VL, sin embargo se puede notar que la mayor cantidad de materia de la galaxia satélite se queda concentrada muy cerca a la VL, situación que no se aprecia en las observaciones.

En esta simulación se encuentra una cola de marea que crece hasta los 200 kpc que comparada con la de 500 kpc de la simulación de 12% es pequeña. Y que se presenta como una estructura de disco perpendicular al disco de la VL, esto se puede apreciar de manera más clara en la figura 4.10 en la que se aprecia la posición final de los objetos en el plano X-Z, sin embargo se puede notar que la mayor cantidad de materia de la galaxia satélite queda concentrada muy cerca de la VL.

Figura 4.10: Vía Láctea de color azul (disco) y morado (bulbo), partículas galaxia satélite rojo disco y verde bulbo. En esta imagen se aprecia una cola de marea con gran densidad de partículas que crece hasta los 200 kpc aproximadamente que comparada con la de 500 kpc de la simulación de 12% es pequeña.

Simulaciones y análisis de Resultados 85

Figura 4.11: Posición final en el plano X-Y de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 12% y las azules (disco) y moradas (bulbo) de la VL.

Figura 4.12: Posición final en el plano Y-Z de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 12% y las azules (disco) y moradas (bulbo) de la VL. 86 Colisión de dos galaxias de Disco

4.3 Análisis simulaciones baja resolución

Partiendo de los resultados obtenidos en las simulaciones de baja resolución se determinó teniendo en cuenta las características morfológicas y el número de partículas que constituyen las colas de marea resultantes que el parámetro de impacto óptimo para producir el mayor desprendimiento de materia de la galaxia satélite es un parámetro de 45 kpc. Teniendo presente que el objetivo es que los remanentes de la colisión se dispersen lo más posible y así se alejen de las vecindades de la VL, se debe tener presente que en estos remanentes es donde se encuentra la materia prima para la eventual constitución de las galaxias enanas satélite.

4.4 Resultados simulaciones de alta resolución

Ya definido el parámetro de impacto de 45 kpc para las cuatro diferentes galaxias satélite (4%, 8%, 12% y 16%) se procedió a realizar las simulaciones de alta resolución con los parámetros mostrados en la tabla 3.6.

4.4.1 Simulaciones de alta resolución galaxia satélite 4%

A continuación se presentan los resultados obtenidos para la simulación de colisión entre la masa del 4 % de la VL y la galaxia anfitriona, los datos empleados en esta simulación son los mostrados en la tabla 3.6.

En esta simulación se encuentra que el tiempo de evolución de la misma es hasta 8.8 Gyr, situación que difiere un poco con las simulaciones de baja resolución. Para el tiempo de 8.8 Gyr la galaxia satélite del 4% solo ha tenido una interacción con la galaxia anfitriona, situación que no ha permitido la desintegración de la galaxia satélite, haciendo una proyección se estima que el tiempo de desintegración de esta galaxia se presentaría en un tiempo superior a un tiempo de Hubble lo cual no sería compatible con nuestro conocimiento del universo. Sin embargo como se aprecia en la figura 4.13 se alcanza a generar una alineación de los remanentes de la galaxia satélite sobre la estructura de la VL, teniendo presente que estos no alcanzan a formar una galaxia enana.

Simulaciones y análisis de Resultados 87

Figura 4.13: Posición final en el plano X-Z de los componentes del sistema después de una simulación de 8,8 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 4% y las azules (disco) y morado (bulbo) de la VL, Se debe tener presente que por la distribución espacial de las partículas, estas no formarían un grupo.

Figura 4.14: Posición final en el plano X-Y de los componentes del sistema después de una simulación de 8,8 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 4% y las azules (disco) y moradas (bulbo) de la VL. 88 Colisión de dos galaxias de Disco

Figura 4.15: Posición final en el plano Y-Z de los componentes del sistema después de una simulación de 8.8 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 4% y las azules (disco) y moradas (bulbo) de la VL.

4.4.2 Simulaciones de alta resolución galaxia satélite 8%

A continuación se presentan los resultados obtenidos para la simulación de colisión entre la masa del 8% de la VL y la galaxia anfitriona, los datos empleados en esta simulación son los mostrados en la tabla 3.6.

En esta simulación ocurre una situación similar al presentado en la simulación del 4%, sin embargo en esta los eventos ocurren un poco más rápido, ya que para el momento de terminar la simulación ya la galaxia del 8% ha interactuando por segunda vez con la VL, sin embargo se puede inferir que para un tiempo de Hubble la galaxia satélite aún no se ha fragmentado por completo, esto indicando que no es físicamente posible este evento en el universo conocido. En este caso también se puede apreciar un alineamiento de las partículas de la galaxia satélite sobre la VL, no obstante estas partículas no formarían una galaxia enana Como se aprecia en la figura 4.16.

Simulaciones y análisis de Resultados 89

Figura 4.16: Posición final en el plano X-Z de los componentes del sistema después de una simulación de 9.5 Gyr, las partículas en rojo (Disco) y verde (Bulbo) son las correspondientes a la galaxia satélite de 8% y las azules (Disco) y morado (Bulbo) de la VL, Se debe tener presente que por la distribución espacial de las partículas, estas no formarían un grupo.

Figura 4.17: Posición final en el plano X-Y de los componentes del sistema después de una simulación de 9,5 Gyr, las partículas en rojo (Disco) y verde (Bulbo) son las correspondientes a la galaxia satélite de 8% y las azules (Disco) y moradas (Bulbo) de la VL. 90 Colisión de dos galaxias de Disco

Figura 4.18: Posición final en el plano Y-Z de los componentes del sistema después de una simulación de 9.5 Gyr, las partículas en rojo (Disco) y verde (Bulbo) son las correspondientes a la galaxia satélite de 8% y las azules (Disco) y moradas (Bulbo) de la VL.

4.4.3 Simulaciones de alta resolución galaxia satélite 12%

A continuación se presentan los resultados obtenidos para la simulación de colisión entre la masa del 12% de la VL y la galaxia anfitriona, los datos empleados en esta simulación son los mostrados en la tabla 3.6.

En este caso se presenta la desintegración total de la galaxia satélite y además este evento se presenta en un intervalo de tiempo que es compatible con el tiempo estimado para la edad del universo. Como se puede apreciar en la figura 4.19 la materia que constituye la galaxia satélite se encuentra dispersa alrededor de la VL, se generan colas de marea que tienen una tendencia perpendicular al disco de la VL, no obstante la mayor parte de las partículas de galaxia satélite se encuentran muy amontonadas cerca de la VL, a pesar de esto esta simulación permite un análisis en busca de grupos que se podrían denominar galaxias enanas satélite de la VL.

Simulaciones y análisis de Resultados 91

Figura 4.19: Posición final en el plano X-Z de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 12% y las azules (disco) y morado (bulbo) de la VL, por la distribución espacial de las partículas, estas podrían formar grupos (galaxias enanas satélite).

Figura 4.20: Posición final en el plano X-Y de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (Disco) y verde (Bulbo) son las correspondientes a la galaxia satélite de 12% y las azules (Disco) y moradas (Bulbo) de la VL. 92 Colisión de dos galaxias de Disco

Figura 4.21: Posición final en el plano Y-Z de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 12% y las azules (disco) y moradas (bulbo) de la VL.

4.4.4 Simulaciones de alta resolución galaxia satélite 16%

A continuación se presentan los resultados obtenidos para la simulación de colisión entre la masa del 16% de la VL y la galaxia anfitriona, los datos empleados en esta simulación son los mostrados en la tabla 3.6.

En primera instancia se aprecia que los escombros se encuentran a lo largo de una franja que presenta perpendicularidad respecto al plano de la VL, segundo y no menos importante se puede apreciar que la galaxia satélite se desgarro completamente, el tercer punto es que el tiempo de ocurrencia del evento es compatible con un tiempo de Hubble, pese a esto se aprecia un gran amontonamiento de materia de la galaxia satélite en vecindades de la VL lo cual no es compatible con las observaciones.

Simulaciones y análisis de Resultados 93

Figura 4.22: Posición final en el plano X-Z de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 16% y las azules (disco) y morado (bulbo) de la VL, Se debe tener presente que por la distribución espacial de las partículas, estas podrían formar grupos (galaxias enanas satélite).

Figura 4.23: Posición final en el plano X-Y de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 16% y las azules (disco) y moradas (bulbo) de la VL. 94 Colisión de dos galaxias de Disco

Figura 4.24: Posición final en el plano Y-Z de los componentes del sistema después de una simulación de 10 Gyr, las partículas en rojo (disco) y verde (bulbo) son las correspondientes a la galaxia satélite de 16% y las azules (disco) y moradas (bulbo) de la VL.

4.5 Análisis de las simulaciones de alta resolución

Para las simulaciones del 4% y 8% se puede decir que si bien es cierto que al cambiar el parámetro de impacto el cambio porcentual de la velocidad inicial y el del parámetro gravitacional no es superior al 5%, sin embargo el cambio en el porcentaje de masa si genera una disminución en las fuerzas gravitacionales entre las galaxia satélite y la galaxia anfitriona, esto causando una disminución en la aceleración adquirida por la galaxia satélite según la segunda ley de Newton, de esta forma el tiempo de ocurrencia de los eventos se dilata, estas colisiones se comportan de manera similar a las del 12% y 16 %, pero tienen un tiempo mucho más prolongado para su ocurrencia, un tiempo que incluso llega ser mayor a la edad del universo.

En el caso de las simulaciones del 12% y 16% se tiene una primera compatibilidad con la concepción de la cosmología moderna ya que el evento ocurre en un tiempo inferior a un tiempo de Hubble, es decir ocurre en un tiempo que es compatible con la edad del universo, sin embargo para poder determinar si en este evento se pueden generar Simulaciones y análisis de Resultados 95

galaxias enanas satélite se debe realizar un análisis de grupos que permita encontrar estas estructuras dentro de los remanentes de la colisión.

4.5.1 Análisis de grupos

A continuación se procede a buscar grupos de partículas dentro de los remanentes de la colisión los cuales se puedan definir como candidatos a galaxia satélite, para esta búsqueda se utiliza el programa amigos de amigos (FoF), para así poder determinar su dinámica y compararla con la reportada por las observaciones.

El código FoF es un código en C ++ desarrollado por la universidad de Washington, FoF es un buscador simple de grupos, que utiliza el método amigos-de-amigos para encontrar grupos. Una partícula pertenece a un grupo de amigos-de-amigos si está dentro de alguna longitud de enlace de cualquier otra partícula en el grupo. En el método de amigos-de- amigos (FoF), se especifica una longitud de enlace, y se identifican todos los pares de partículas con una separación igual o menor a la longitud de enlace. Tales pares se designan los amigos y los grupos se definen como conjuntos de partículas que están conectadas por una o más relaciones de amistad, es decir, los amigos de los amigos.

A continuación se muestran los resultados obtenidos con el método FoF para las simulaciones del 16% y 12% respectivamente.

4.5.2 Grupos galaxia satélite 16 %

Se realizaron análisis con el código FoF para diferentes longitudes de enlace haciendo un barrido desde una longitud de 0,1 kpc hasta una longitud de 1kpc, y se obtuvieron los resultados que se aprecian en la tabla 4.1.

Longitud de enlace Número de grupos (kpc)

0. 3 1 0. 4 7

0. 5 11 0. 6 14 0. 7 17 0. 8 25 0. 9 28 1. 0 35

Tabla 4.1: Número de grupos respecto a la distancia de enlace elegida. 96 Colisión de dos galaxias de Disco

La distancia de enlace que genera un número de grupos que se hace más compatible con las observaciones es la de 1,0 kpc, estos grupos tienen una masa media de 1.49x10 7⊙ y se pueden apreciar en las figuras 4.25, 4.26 y 4.27.

Figura 4.25: Grupos generados después de la colisión en el plano X-Y para la galaxia satélite del 16% de masa de la VL.

Figura 4.26: Grupos generados después de la colisión en el plano X-Z para la galaxia satélite del 16% de masa de la VL. Simulaciones y análisis de Resultados 97

Figura 4.27: Grupos generados después de la colisión en el plano Y-Z para la galaxia satélite del 16% de masa de la VL.

4.5.3 Perfil de densidad Vía Láctea 16%

Para determinar la validez de la colisión se procede a realizar la gráfica del perfil de densidad del disco de la VL para establecer si esta estructura se ha calentado después de la colisión, es decir si se ha modificado el perfil de tipo exponencial del disco, ver figura 4.28.

Figura 4.28: Perfil de densidad del disco de la VL en azul antes de la colisión y rojo después, como se aprecia el perfil de densidad en ambas situaciones se comporta de manera similar en la región de interés entre 0 kpc y 16 kpc. 98 Colisión de dos galaxias de Disco

Como se aprecia en la figura 4.28 la VL conserva su perfil de densidad exponencial en la región de interés lo que es favorable para este estudio sin embargo como se aprecia en las figuras 4.25 a 4.27 la VL presenta un calentamiento en su escala vertical, razón por la cual se procede a realizar el estudio de la dispersión de velocidades en z para el disco de la Vía Láctea.

4.5.4 Dispersión de velocidades en Z para la Vía Láctea 16%

Figura 4.29: Dispersión de velocidades en Z (s) para la VL respecto a la coordenada radial del disco (R) en morado antes de la colisión y verde después, se observa un aumento en la dispersión de velocidades en z antes y después de la colisión.

En este caso se encuentra un aumento medio de la dispersión de velocidades en z en promedio del 96%, lo cual indica un calentamiento en el disco después de la colisión.

4.5.5 Plano de la Vía Láctea 16%

A continuación se presenta el plano de la VL, plano que se calculó usando el método de mínimos cuadrados para dos variables.

Simulaciones y análisis de Resultados 99

Figura 4.30: Plano de la VL obtenido con las posiciones de las partículas del disco de la Vía Láctea usando el método de mínimos cuadrados para dos variables.

Usando este método se encontró que la ecuación del plano de la VL es:

Z = −0.3712x + 0.4746y − 317.6422 (4.1)

Con las ecuaciones de este plano y el plano de los grupos del 16% se calculará el ángulo entre los dos.

4.5.6 Dispersión de velocidades grupos del 16%

En esta sección se muestran las dispersiones de velocidades respecto a la visual de las diferentes galaxias satélites encontradas en este estudio y se contrastan estos resultados con los datos reportados por las observaciones (Walker et al. 2008). Para llevar a cabo el cálculo de estas dispersiones se realizó la rotación de la Vía Láctea y de las galaxias satélites o grupos empleando una matriz de rotación, los ángulos de rotación que se emplearon fueron los vistos para la Vía Láctea en los planos YZ y XZ, luego se procedió a centrar la VL y los grupos respecto al origen, se calculó la posición del centro de masa de cada uno de los grupos respecto al sol el cual se ubicó en la coordenada (0,-8,0), se convirtió este vector r de posición en un vector unitario y a continuación se realizó el producto punto entre y el vector de velocidad para = ( , , ) ͦ̂ = ( 3, 4, 5) determinar la proyección ͦ̂ Vproy ͬ deͭ ͮ sobre y por último se calculóͪľ ͐ la͐ desviación͐ estándar de estas proyecciones que esͪľ lo que se denͦ̂ omina dispersión de velocidades S.

+-*4 = 3 + 4 + 5 (4.2)

͐ ͐ ͬ ͐ ͭ ͐ ͮ 100 Colisión de dos galaxias de Disco

Figura 4.31: En esta gráfica se aprecia la dispersión de velocidades (S) de los 35 grupos encontrados en los remanentes de la galaxia satélite dejados después de la colisión, la media de dispersión para los grupos es de 234 km/s.

Figura 4.32: En esta imagen se aprecian las dispersiones de velocidades de algunas galaxias enanas satélites de la VL respecto a R en (pc), se aprecian dispersiones de velocidad del orden de 10 km/s (Walker et al. 2008).

La dispersión de velocidades de los grupos encontrados en la simulación presenta una diferencia de 224 km/s respecto a la dispersión reportada por las observaciones. Simulaciones y análisis de Resultados 101

4.5.7 Plano grupos 16%

A continuación se presenta el plano generado por los grupos del, este se calculó usando el método de mínimos cuadrados para dos variables.

Figura 4.33: Plano de los grupos del 16% obtenido con las posiciones de los grupos encontrados usando el método de mínimos cuadrados para dos variables.

Usando este método se encontró que la ecuación del plano de los grupos es:

Z = −0.1833x + 0.2400y − 217.2953 (4.3)

4.5.8 Planos y ángulo entre planos 16%

En este apartado se muestran el plano de la VL y el plano de los grupos en tres dimensiones.

x

Figura 4.34: Partículas rojas Vía Láctea, puntos azules grupos, plano verde Vía Láctea y plano morado grupos, longitudes en kpc. 102 Colisión de dos galaxias de Disco

y

Figura 4.35: Partículas rojas Vía Láctea, puntos azules grupos, plano verde Vía Láctea y plano morado grupos, longitudes en kpc.

Figura 4.36: Partículas rojas Vía Láctea, puntos azules grupos, plano verde Vía Láctea y plano morado grupos, se aprecia un ángulo entre los dos planos inferior a los 88°, longitudes en kpc.

Partiendo delas ecuaciones 4.1 y 4.3 se dedujeron los vectores normales y a los ͥ ͦ planos y se determinó el ángulo entre estos dos vectores que resulta ser el͢ʌʌ ʌʌmismoľ ͢ʌʌʌʌľ ángulo entre los planos:

ͥ ∙ ͦ = (4.4) ͥ ͦ ͢ʌʌʌʌľ ͢ʌʌʌʌľ El ángulo encontrado entre el plano̻ͦ͗̽ͣͧ de la ʦVL y el ʧplano de los grupos es de ͢ ͢ = 14.3°, ángulo que es muy pequeño comparado con el ángulo de 88° reportado por las observaciones (Kroupa 2014). Este ángulo se presenta por el amontonamiento que tienen los grupos sobre la VL.

Simulaciones y análisis de Resultados 103

Por la pérdida de estabilidad que presenta la VL en la colisión con la galaxia satélite del 16 % se descartan las condiciones iniciales de este evento como posibles causantes del origen del disco de satélites de la VL (DoS).

4.5.9 Grupos para la galaxia satélite 12%

De manera análoga a lo realizado con la simulación del 16% se procede a buscar grupos en la colisión de la VL con la galaxia satélite del 12% de su masa, usando el código FoF para distintas distancias de enlace se encontraron los resultados registrados en la tabla 4.2. Longitud de Número de grupos enlace (kpc)

0. 4 5 0. 5 14 0. 6 17 0. 7 19 0. 8 24

0. 9 30 1. 0 32

Tabla 4.2: Número de grupos respecto a la distancia de enlace elegida.

Con una distancia de enlace de 1,0 kpc se encuentran 32 grupos los cuales son candidatos a ser galaxias enanas satélite de la VL, este número de grupos se aproxima al número de galaxias enanas satélite observadas, estos grupos tienen una masa media de

1.48x10 7⊙ y pueden ser apreciados en las figuras 4.37, 4.38 y 4.39.

104 Colisión de dos galaxias de Disco

Figura 4.37: Grupos generados después de la colisión en el plano X-Y para la galaxia satélite del 12% de masa de la VL.

Figura 4.38: Grupos generados después de la colisión en el plano X-Z para la galaxia satélite del 12% de masa de la VL.

Simulaciones y análisis de Resultados 105

Figura 4.39: Grupos generados después de la colisión en el plano Y-Z para la galaxia satélite del 16% de masa de la VL.

4.5.10 Perfil de densidad Vía Láctea 12%

En este caso el perfil de densidad exponencial de la VL también se mantiene, sin embargo como se aprecia en las 4.38 y 4.39 el disco presenta un calentamiento en su componente vertical razón por la cual se hace necesario analizar la dispersión de velocidades de la galaxia en z.

Figura 4.40: Perfil de densidad del disco de la VL en azul antes de la colisión y rojo después, como se aprecia el perfil de densidad en ambas situaciones se comporta de manera similar en la región de interés entre 0 kpc y 16 kpc.

106 Colisión de dos galaxias de Disco

Como se aprecia en la figura 4.40 la VL conserva su perfil de densidad exponencial en la región de interés lo que es favorable para este estudio sin embargo como se aprecia en las figuras 4.38 y 4.39 la VL presenta un calentamiento en su escala vertical, razón por la cual se procede a realizar el estudio de la dispersión de velocidades en z para el disco de la Vía Láctea.

4.5.11 Dispersión de velocidades en Z para la Vía Láctea 12%

Figura 4.41: Dispersión de velocidades en Z (s) para la VL respecto a la coordenada radial del disco (R) en morado antes de la colisión y verde después, se observa un aumento en la dispersión de velocidades en z antes y después de la colisión.

En este caso se encuentra un aumento medio de la dispersión de velocidades en z es en promedio del 80%, lo cual indica un calentamiento en el disco después de la colisión.

4.5.12 Plano de la Vía Láctea 12%

Figura 4.42: Plano de la VL obtenido con las posiciones de las partículas del disco de la Vía Láctea usando el método de mínimos cuadrados para dos variables. Simulaciones y análisis de Resultados 107

Usando este método se encontró que la ecuación del plano de la VL es:

Z Ɣ Ǝ0.0106x Ǝ 0.1586y ƍ 3.0880 ʚ4.5ʛ

Con las ecuaciones de este plano y el plano de los grupos del 12% se calculará el ángulo entre los dos.

4.5.13 Dispersión de velocidades en grupos del 12%

En esta sección se muestran las dispersiones de velocidades de las diferentes galaxias satélites encontradas en este estudio y se contrastan estos resultados con los datos reportados por las observaciones (Walker et al. 2008).

Figura 4.43: En esta gráfica se aprecia la dispersión de velocidades media de los 32 grupos encontrados en los remanentes de la galaxia satélite dejados después de la colisión, el promedio de estas dispersiones es de 242 km/s.

Como se aprecia en las Figura 4.32 la dispersión de velocidades de las galaxias satélites de la VL es baja comparada con los valores de dispersión de velocidad encontrados en este estudio para las galaxias satélites, se debe tener presente que al presentarse una diferencia de aproximadamente 220 km/s en la dispersión nos indica que las partículas constituyentes de las galaxias satélites de la VL se encuentran ligadas de manera gravitacional, mientras que las partículas que constituyen las galaxias satélites encontradas en este estudio no se encuentran ligadas, siendo posible que se presente simplemente un solapamiento de estas partículas, pero que aún no sean sistemas estables.

108 Colisión de dos galaxias de Disco

4.5.14 Plano grupos 12%

Figura 4.44: Plano de los grupos del 12% obtenido con las posiciones de los grupos encontrados usando el método de mínimos cuadrados para dos variables.

Usando este método se encontró que la ecuación del plano de los grupos es:

= 0.1377 − 0.4586 + 94.8128 (4.6)

4.5.15 Planos yͮ ángulo ͬentre planosͭ 12%

En este apartado se presentan en tres dimensiones los planos correspondientes a la VL y a los grupos.

x y

Figura 4.45: Partículas rojas Vía Láctea, puntos azules grupos, plano verde Vía Láctea y plano morado grupos, longitudes en kpc. Simulaciones y análisis de Resultados 109

y

Figura 4.46: Partículas rojas Vía Láctea, puntos azules grupos, plano verde Vía Láctea y plano morado grupos, longitudes en kpc.

Figura 4.47: Partículas rojas Vía Láctea, puntos azules grupos, plano verde Vía Láctea y plano morado grupos. Se aprecia un ángulo entre los dos planos inferior a los 88°.

Partiendo de lo explicado en el apartado 4.5.8 y de la ecuación 4.4 se calculó el ángulo entre el plano de la VL y el plano de los grupos del 12%, el cual tiene un valor de = 17.4°, este ángulo es pequeño comparado con el de 88° reportado por las observaciones (Kroupa 2014). Este ángulo se presenta por el amontonamiento que tienen los grupos sobre la VL.

Si se compara la distribución espacial de las galaxias satélite encontradas en este estudio en la Fig. 4.38 con la distribución espacial observada (Kroupa 2014) se puede apreciar una gran diferencia respecto a lo que se refiere a la distancia a la que se encuentran las galaxias satélites respecto al plano de la VL.

110 Colisión de dos galaxias de Disco

Figura 4.48: Izquierda (Kroupa 2014), derecha simulación 12%, si bien se aprecia en la simulación una estructura tipo disco, se ven muy amontonadas las galaxias satélites sobre el plano de la VL, en las observaciones se aprecian distancias de hasta 200 kpc respecto al disco de la VL, En la simulación esta distancia llega ser máximo de 20 kpc respecto al disco.

5. Conclusiones

En este trabajo se realizaron simulaciones numéricas de N-cuerpos para la colisión de dos galaxias de disco con el fin determinar las condiciones iniciales de simulación que permitirían recrear la dinámica y distribución espacial de las galaxias satélite de la Vía Láctea. Se plantea un modelo en el que estas galaxias hayan tenido origen a través de la colisión de dos galaxias de disco.

En primera instancia se especifican las relaciones de masa entre las galaxias a interactuar (Vía Láctea y galaxia fantasma), para lo cual se determinan relaciones de 4%, 8%, 12% y 16%, estos valores se eligen con el fin de evitar inestabilidades después de la colisión en el disco de la Vía Láctea. Elegidas estas relaciones se disponen las galaxias de manera perpendicular una respecto a la otra, buscando que los escombros de marea queden dispuestos de manera perpendicular al plano de la Vía Láctea.

Con el fin de construir las denominadas galaxias fantasmas se desarrolla un modelo que permite construir estos objetos como múltiplos de una galaxia conocida, en términos de sus características geométricas (longitudes de escala). Los resultados del modelo son positivos respecto a simulaciones en la que los objetos se crean en términos del momentum angular, ya que la diferencia porcentual media entre los elementos característicos de la galaxia es del 4.74%.

Para el camino de la galaxia fantasma se elige un trayecto de tipo parabólico progrado que permita el mayor desprendimiento de materia posible de la galaxia fantasma que genere grandes colas de marea. Para definir esta trayectoria se requiere de un valor denominado parámetro de impacto el cual se determinó a partir de simulaciones de baja resolución (150000 partículas) las cuales permitieron enmarcar este en un valor de 45 kpc.

Determinadas estas condiciones iniciales se realizaron cuatro simulaciones una por cada relación de masa. Para estas simulaciones se encontraron dos comportamientos diferentes, para el grupo de simulaciones del 4% y 8% se encontró que bajo las condiciones iniciales planteadas anteriormente el tiempo de ocurrencia del evento es superior al tiempo de Hubble, esto debido a la relación de masas que genera fuerza de interacción entre las

112 Colisión de dos galaxias de Disco galaxias que son menores que las experimentadas por las relaciones del 12% y 16% respectivamente, esto como consecuencia genera aceleraciones menores sobre las galaxias fantasmas, lo cual por supuesto trae como consecuencia un aumento en el tiempo de ocurrencia de los eventos, situación que permite descartar estas dos condiciones iniciales como posibles causantes del origen del disco de galaxias satélite de la Vía Láctea

Para probar la validez de las simulaciones del 12% y 16% se procede a analizar el perfil de densidad de la Vía Láctea para garantizar su estabilidad, en este análisis se encontró que la Vía Láctea conserva su perfil de densidad exponencial en las dos simulaciones, sin embargo, la estructura de disco de la Vía Láctea presenta un calentamiento en la dirección de la longitud de escala vertical (z). Para determinar de qué orden es este calentamiento se hizo el estudio de la dispersión de velocidades en z para lo cual se encontró que para la simulación del 12% se presenta un aumento promedio en la dispersión de velocidades del 80% y para la simulación del 16% se encontró que es en promedio del 96%.

En las simulaciones del 12% y 16% se encontraron grupos que podrían ser denominados como galaxias satélites, para la simulación del 12 % se encontró un máximo de 32 grupos con una distancia de enlace de 1kpc y para la simulación de 16% se encontró un máximo de 35 grupos con una distancia de enlace de 1kpc, este número de grupos es compatible con el número de grupos reportado por las observaciones. La masa promedio de estos grupos para las simulaciones del 12% y 16% es de 1,485x10 7 M⊙ valor de masa que es conforme con las observaciones.

Se estudió la dispersión de velocidades para los grupos obtenidos para las dos simulaciones y se halló para los grupos de la simulación del 12% que la dispersión de velocidades media respecto a la visual es de 241 km/s y para la simulación del 16% se encontró una dispersión de velocidades media respecto a la visual de 234 km/s, esto no es coincidente con las observaciones que reportan dispersiones de velocidades medias del orden de 10 km/s. Este valor tan elevado en las dispersiones de velocidades encontradas en las simulaciones es un indicativo de que estas estructuras muy posiblemente no se encuentran ligadas gravitacionalmente.

Para establecer la disposición espacial de los grupos se realizó el cálculo del plano en el cual se encuentran dispuestos los grupos usando el método de mínimos cuadrados para dos variables, para la simulación del 12% se estableció que este plano forma un ángulo de 17.4° respecto al plano de la Vía Láctea y para la simulación del 16% se determinó que este ángulo es de 14.3° esta escenario no es concurrente con las observaciones que definen que el ángulo entre el plano de los grupos y la Vía Láctea es de 88°, esta situación se presenta por el amontonamiento de los grupos sobre la VL, ya que en las observaciones se encuentra que hay galaxias satélites que se encuentran a distancias respecto al centro de Conclusiones 113

la VL de hasta 200 kpc, mientras que en las simulaciones la mayor distancia respecto al centro de la VL es de 20 kpc. Para estas simulaciones se podría afirmar que no hay estructura de disco de galaxias satélites (DoS).

Se puede decir que bajo las condiciones iniciales planteadas para estas simulaciones no es posible generar las características dinámicas y espaciales de lo que se denomina el disco de galaxias satélite la Vía Láctea (DoS), sin embargo, esto no descarta por completo el modelo propuesto para este estudio ya que si bien no refleja los resultados esperados, se pueden proponer trabajos futuros para realizar simulaciones y modificar los parámetros de impacto y las relaciones de masa para poder enmarcar las condiciones físicas bajo las cuales es posible este evento. Los resultados de este trabajo permiten acotar aún más el espacio de condiciones iniciales orbitales y de masa de eventuales progenitores del DoS de la VL.

A. Apéndice: Redimensionamiento de estructuras (Cálculos)

Como se explicó en el capítulo 3 se realizó un redimensionamiento de las diferentes estructuras que constituyen la galaxia satélite (Galaxia fantasma), a continuación se presentan los cálculos que permiten llegar a las expresiones que definen las dimensiones del halo de materia oscura, el disco y el bulbo de la galaxia satélite dependiendo del porcentaje de relación de masas entre la VL y la galaxia satélite.

A.1. Redimensionamiento Halo de materia oscura

Para determinar las expresiones que definen el halo de materia oscura de la galaxia satélite se supone una región con una densidad de para la cual se produce el 200 - colapso del halo, donde se define como: -   3 ͦ( )

- = ( . 1) ͂8 ͮ  ̻ ́ Ahora definimos la masa del halo en términos de su densidad volumétrica y su - volumen asumiendo una simetría esférica, obteniendo que: 

4 ͧ 200 ͦͤͤ - ͦͤͤ = ( . 2) ͦ 3  Reemplazando se tiene que: ͇ ̻ -  4 ͧ 200 3 ͦ( ) ͦͤͤ ͦͤͤ = ( . 3) ͦ 3 ͂8 ͮ ͇ ʦ ʧ ̻ Realizando el proceso algebraico se tiene que: ́

100 ͧ ͦ( ) ͦͤͤ ͦͤͤ = ( . 4) ͦ ͂ ͮ ͇ ̻ ́ Ahora teniendo presente que la velocidad virial del halo es una velocidad circular en el radio virial , entonces se tiene: ͦͤͤ ͦ

116 Colisión de dos galaxias de Disco

ͦ ͦͤͤ ͦͤͤ = ( . 5) ͇́ͦͤͤ ͪ ̻ Por lo tanto: ͦ

ͦͤͤ ͦͤͤ = ( . 6) ͇́ͧ ͦ ͦͤͤ ̻ ͪ ͧ ͧ ͧ = ͦͤͤ ( . 7) ͦͤͤ ́ ͇ͪ ͦ ͦͤͤ ̻ ͪ Reemplazando esta expresión en la ecuación de M se tiene:

100 ͦ( ) ͧ ͧ ͦͤͤ = ͪ ( . 8) ͂ ͮ ́ ͇ͦͤͤ ͇ ʦ ʧ ̻ ́ ͪ ͪ ͦ = ͦͤͤ ( . 9) ͦͤͤ 100 ́ͪͧ ͦ( ) ͇ ̻ ́ ͂ ͮ Desarrollando se tiene que:

ͧ ͦͤͤ ͦͤͤ = ( . 10) 10 ͪ ( ) ͇ ̻ ́͂ ͮ Donde es la misma velocidad circular de la ecuación (3.1), ahora teniendo esta ͦͤͤ  expresión se define el radio virial como sigue: ͪ ͦͤͤ ͪ ͦ ͧ ͦͤͤ ͦͤͤ = ͧ ( . 11) ͦͤͤ́ 10 ͪ ( ) ͦ ʦ ʧ ̻ ͪ ́͂ ͮ Así es: ͦͤͤ ͦ ͦͤͤ ͦͤͤ = ( . 12) 10ͪ ( ) ͦ ̻ ͂ ͮ Las ecuaciones (A.10) y (A.12) son las que se emplearon para redimensionar el halo de materia oscura de la galaxia satélite (Galaxia fantasma) a partir del porcentaje de relación de masas entre las dos galaxias de disco interactuantes.

Apéndice A. Redimensionando Estructuras 117

A.2. Redimensionamiento del Disco

En este apartado se muestran los cálculos requeridos para redimensionar el disco de la galaxia satélite (Galaxia fantasma), Asumiendo que es igual para la galaxia satélite ∑ͤ como para la VL, definiendo como sigue: ∑ͤ

∑ͤ = ͦ ( . 13) 2͇ Ahora, como es el mismo para ambas galaxias se tiene que: ̻ ∑ͤ ͌

 ͦ = ͦ ( . 14) 2͇  2͇  ̻ ͌ ͌ Como se conoce que la relación porcentual de masas de los discos de la galaxia satélite y de la VL es k definida como sigue:  ͇ ͇ = ( . 15) ͇  ͟ ̻ ͇ Sustituyendo (A.15) en (A.14) se tiene:

  ͦ = ͦ ( . 16) 2͇͟  2͇  ̻ ͌ ͌ Realizando los correspondientes cálculos se obtiene (A.16)

 =  ( . 16) ͌ ͌ ̻ √͟ Esta expresión permite determinar la longitud de escala radial de la galaxia satélite en  términos de la longitud de escala de la VL .  ͌ ͌ Para determinar la escala de longitud vertical del disco de la galaxia satélite se supone ͤ un comportamiento análogo al de la galaxia anfitriona respecto a sus dimensiones͔ de escala vertical y escala radial , se asume que el cociente definido por la ecuación ͤ  (A.17) es el mismoͮ para las dos galaxias.͌ 

 = ( . 17) ͌ͤ  ̻ ͮ 118 Colisión de dos galaxias de Disco

Para la igualdad se tiene:

 =  ( . 18) ͌ͤ ͌ͤ ̻ ͮ ͔ Realizando el proceso algebraico se tiene que:

 = ͤ ( . 20) ͌ ͔ ̻  Las ecuaciones (A.16) y (A.20) son las expresiones que se usaron para simular la galaxia satélite Galaxia fantasma) a partir del porcentaje de relación de masas entre las dos galaxias de disco interactuantes k y la relación entre las dimensiones geométricas vertical y radial .  A.3. Redimensionamiento del Bulbo

Para determinar las dimensiones del bulbo de la galaxia satélite en términos de las dimensiones del bulbo de la VL usando el porcentaje de relación de masa ecuación (A.21) y el perfil de densidad de Hernquist ecuación (A.22). ì

=  ( . 21) ͇ ì ̻ ͇ = ( . 22)  2 ( + )ͧ ͇͕  ̻ Ahora teniendo presente que r se puede escribir:ͦ ͕ ͦ = ( . 23)

Donde a es la longitud de escala radial delͦ bulbo͕͢ y n es el cociente relación entre a y̻ r se tiene que:

= ( . 24)  2 ( + )ͧ ͇͕  ̻ Cancelando a y factorizando en el denominador͕͢ ͕ ͕͢ se tiene:

= ( . 25)  2 ( 1 + )ͧ ͇  ̻ Entonces: ͢ ͕ʚ ͢ʛ

= ( . 26)  2 ͧ(1 + )ͧ ͇  ̻ ͕͢ ͢ Apéndice A. Redimensionando Estructuras 119

Ahora asumiendo que la densidad de masa del bulbo de la galaxia satélite es igual a la  densidad del bulbo de la VL ; con aL y aS las correspondientes longitudes de escala   radiales. 

 =  ( . 27) 2 ͇ͧ(1 + )ͧ 2 ͇ͧ(1 + )ͧ  ̻ ͕͢ ͢ ͕͢ ͢ Y haciendo uso de se encuentra: ì  =  ( . 28) 2 ͇ͧ(1 + )ͧ 2 ͇ͧ(1 + )ͧ ̻ ͕͢ ͢ ͢ ͕ ͢ Por lo tanto:

1 1

ͧ = ͧ  En consecuencia; ͕ ͕

w =  ( . 29) ͕ ͕ ̻

Usando la expresión (A.29) se redimensionoǭ el bulbo de la galaxia satélite  a partir de la relación de masa respecto a la VL. ͕ Con las expresiones demostradas anteriormente se construyeron las galaxias satélites (Galaxia fantasma) empleadas para las diferentes simulaciones.

B. Anexo: Resultados de simulación

A continuación se presentan los resultados de simulación, se muestran imágenes con las distintas posiciones de la galaxia satélite y de la VL, además en estas imágenes se pueden apreciar los planos X-Y, X-Z y Y-Z, estos planos se pueden apreciar por columnas para las simulaciones del 4%, 8%, 12% y 16%.

122 Colisión de dos galaxias de Disco

B.1. Resultados simulación 4%

Anexo B. Resultados de Simulación 123

124 Colisión de dos galaxias de Disco

Anexo B. Resultados de Simulación 125

126 Colisión de dos galaxias de Disco

Anexo B. Resultados de Simulación 127

B.2. Resultados simulación 8%

128 Colisión de dos galaxias de Disco

Anexo B. Resultados de Simulación 129

130 Colisión de dos galaxias de Disco

Anexo B. Resultados de Simulación 131

132 Colisión de dos galaxias de Disco

B.3. Resultados simulación 12%

Anexo B. Resultados de Simulación 133

134 Colisión de dos galaxias de Disco

Anexo B. Resultados de Simulación 135

136 Colisión de dos galaxias de Disco

Anexo B. Resultados de Simulación 137

B.4. Resultados simulación 16%

138 Colisión de dos galaxias de Disco

Anexo B. Resultados de Simulación 139

140 Colisión de dos galaxias de Disco

Anexo B. Resultados de Simulación 141

142 Colisión de dos galaxias de Disco

Bibliografía

Angus, G. W., Diaferio, A., and Kroupa, P. (2011). Using dwarf satellite proper motions to determine their origin . MNRAS, 416:1401–1409.

Barnes, J. E. 2011, Astrophysics Source Code Library, record ascl:1102.027

Binney, J.; Merrifield, M. (1998). Galactic Astronomy . Princeton: Princeton University Press. ISBN: 978-0-691-02565-0. OCLC 39108765.

Binney. J, Tremaine. S. (2008). Galactic Dynamics , Princeton Univ.Press, U.K. Second Edition. ISBN: 987-0691130279 p. 904

Casallas A. (2012) Estudio de la formación del disco de satélites de la Vía Láctea (DoS) como un grupo de progenitores que entra al halo de la Galaxia. Facultad de Ciencias, Departamento de Física. Universidad Nacional de Colombia. Bogotá.

Casallas A., Cubillos D., Casas Miranda R. (2009). Galaxias Enanas del Grupo Local . Momento-Revista de Física, No 38.

Cassan, A. et al. (2012). “One or more bound planets per Milky Way star from microlensing observations” . Nature 481 (7380): 167–169. arXiv:1202.0903. Bibcode:2012Natur.481..167C.

Chandrasekhar, S. (1943), "Dynamical Friction. I. General Considerations: the Coefficient of Dynamical Friction" , ApJ vol. 97, 255–262.

Collins Dictionary of Astronomy, 2nd Ed.; Harper Collins Publishers; 2000; pp.444,449

Cubillos D. (2013). Estudio de la Formación de Galaxias Enanas Esferoidales Satélites de la Vía Láctea Mediante la Simulación de la Colisión de dos Galaxias. Facultad de Ciencias, Departamento de Física. Universidad Nacional de Colombia. Bogotá.

144 Colisión de dos galaxias de disco

D’Onghia, E., Besla, G., Cox, T. J., and Hernquist, L. (2009). Resonant stripping as the origin of dwarf spheroidal galaxies . Nature, 460:605–607.

Da Costa, G. S. (1999). The Dwarf Spheroidal Galaxies in the Galactic Halo . In Gibson, B. K., Axelrod, R. S., and Putman, M. E., editors, The Third Stromlo Symposium: The Galactic Halo, volume 165 of Astronomical Society of the Pacific Conference Series, page 153.

Diemand, J.; Moore, B. & Stadel, J. (2005). “ Earthmass dark-matter haloes as the first structures in the early ”. Nature 433. pp. 389-391.doi:10.1038/nature10684. PMID 22237108.

Graham, A.W., Worley, C.C. (2008). Inclination- and dust-corrected galaxy parameters: Bulge-to-disc ratios and size-luminosity relations MNRAS, 388, 1708.

Grebel, E. K., Gallagher, J. S., & Harbeck, D. (2003). The Progenitors of Dwarf Spheroidal Galaxies ApJ, vol. 125, 1926

Grebel. E.K, (1998). Star Formation Histories of Local Group Dwarf Galaxies , Arxiv:9806191v2 [astro-ph].

Hernquist, L. (1990). An analytical model for spherical galaxies and bulges . ApJ, 356:359– 364.

Hernquist, L., Spergel, D. N., and Heyl, J. S. (1993). Structure of Merger Remnants. III. Phase-Space Constraints . , 416:415.

Hubble, E. P. (1936). The Realm of the Nebulae . New Haven: Yale University Press. LCCN 36018182.

Huctchmeier.W, I. D. Karachentsev, V. E. Karachentseva, D.I. Makarov. (2006). A Catalog of Neighboring Galaxies . The Astronomical Journal, 132:729–748, August ApJ, 132, 729.

Imamura, Jim (2006). “Mass of the Milky Way Galaxy” . University of Oregon. Archived from the original on March 1, 2007. Retrieved May 10, 2007.

Karachentsev, I. D.; Kashibadze, O. G. (2006). “Masses of the local group and of the M81 group estimated from distortions in the local velocity field” . Astrophysics 49 (1): 3–18. Bibliografía 145

Karttunen. H, Kroger. P, Oja. H, Poutanen. M. And Donner. (2003) M. Fundamental Astronomy. Springer. 5th edition, New York.

Kravtsov, Andrey V.; Borgani, Stefano, Formation of Galaxy Clusters . Annual Review of Astronomy and Astrophysics, vol. 50, p.353-409.

Kroupa, C. Theis, and C. M. Boily (2005). The great disk of milky-way satellites and cosmological sub-structures. Astronomy and Astrophysics, 431(2):517–521.

Kroupa, P (2014). Lessons from the Local Group (and beyond) on dark matter . arXiv:1409.6302v1 [astro-ph.GA].

Łokas, E. L. (2008). The Milky Way’s Neighbors . Academia: Magazine of the Polish Academy of Sciences, 2:12–15.

Lokas, E. L. (2011). Dark matter in dwarf galaxies of the Local Group. ArXiv e-prints.

Longair, M. S. (1998). Galaxy Formation . New York: Springer. ISBN 3-540-63785-0.

Longair, Malcolm S. (2006), The Cosmic Century: A History of Astrophysics and Cosmology , Cambridge: Cambridge University Press, ISBN 9780521474368.

Mashchenko, S., Couchman, H. M. P., and Sills, A. (2005). Modeling Star Formation in Dwarf Spheroidal Galaxies: A Case for Extended Dark Matter Halos . ApJ, 624:726–741.

Mateo, M. (1996). Dwarf Spheroidal Galaxies and the Formation of the Galactic Halo . In Morrison, H. L. and Sarajedini, A., editors, Formation of the Galactic Halo...Inside and Out, volumen 92 of Astronomical Society of the Pacific Conference Series, page 434.

Mateo, M. L. (1998). Dwarf Galaxies of the Local Group . Annual Review of Astronomy & Astrophysics, 36:435–506.

McConnachie, A.W. and Irwin, M. J. (2006). The satellite distribution of M31 . MNRAS, 365:902–914.

Metz, M., Kroupa, P., and Jerjen, H. (2007). The spatial distribution of the Milky Way and Andromeda satellite galaxies . MNRAS, 374:1125–1145.

Metz, M., Kroupa, P., and Jerjen, H. (2007). The spatial distribution of the Milky Way and Andromeda satellite galaxies . MNRAS, 374:1125–1145. 146 Colisión de dos galaxias de disco

Metz, M., Kroupa, P., Theis, C., Hensler, G., and Jerjen, H. (2009b). Did the Milky Way Dwarf Satellites Enter The Halo as a Group? . ApJ, 697:269–274.

Metz. M, Kroupa. P.Did. (2006). The Milky Way Dwarf Satellites enter the Halo as a Group? . Arxiv:0610933v1 [astro-ph]

Navarro. F, Frenk. C, White. S. (1995). The Structure of Cold Dark Matter Halos . Arxiv:astro-ph 9508025

Nichols, M., Colless, J., Colless, M., and Bland-Hawthorn, J. (2011). Accretion of the Magellanic System onto the Galaxy. ApJ, 742:110.

Okazaki, T. and Taniguchi, Y. (2000). Dwarf Galaxy Formation Induced by Galaxy Interactions . ApJ, 543:149–152.

Pawlowski, M. S., Kroupa, P., and de Boer, K. S. (2011). Making counter-orbiting tidal debris. The origin of the Milky Way disc of satellites? Astronomy and Astrophysics, 532:A118.

Pringle and King (2007), Astrophysical Flows, p161. Cambridge University Press.

Quiroga L. (2013) Evolution of Disk Galaxies in ΛCDM. Facultad De Ciencias Exactas Y Naturales, Instituto De Física. Universidad de Antioquia. Medellín.

Cubillos D. (2013). Estudio de la Formación de Galaxias Enanas Esferoidales Satélites de la Vía Láctea Mediante la Simulación de la Colisión de dos Galaxias. Facultad de Ciencias, Departamento de Física. Universidad Nacional de Colombia. Bogotá.

R. Brent Tully; Helene Courtois; Yehuda Hoffman; Daniel Pomarède (2014). “The of galaxies” . Nature 513 (7516).

Ricotti, M. and Gnedin, N. Y. (2005). Formation Histories of Dwarf Galaxies in the Local Group. ApJ, 629:259–267

Rocha. M, Annika. H, Peter. G, and Bullock. J. (2011). Infall Times for Milky Way Satellites From Their Present-Day Kinematics , Astrophysics and Cosmology. Arxiv:1110;04664v1 [astroph.CO].

Rodriguez. A (2001), Astronomía para todos , Universidad nacional de Colombia, Unibiblos, Bogotá. ISBN 958-701-104-X. cap 19. Bibliografía 147

Rubin, Vera C.; Ford, W. Kent, Jr. (1970). “Rotation of the Andromeda Nebula from a Spectroscopic Survey of Emission Regions” . ApJ 159: 379–403.

Sawa, T. and Fujimoto, M. (2005). A Dynamical Model for the Orbit of the M31 and the Origin of the Local Group of Galaxies . PASJ, 57:429–446.

Shapley, H.; Curtis, H. D. (1921). “The Scale of the Universe” . Bulletin of the National Research Council 2: 171–217. Bibcode:1921BuNRC...2..171S.

Sparke, L. S.; Gallagher III, J. S. (2000). Galaxies in the Universe: An Introduction. Cambridge University Press. ISBN 0-521-59740-4.

Spitzer, Lyman, Jr. (1942). The Dynamics of the Interstellar Medium. III. Galactic Distribution. Astrophysical Journal, vol. 95, p.329.

Springel, V. and White, S. D. M. (1999). Tidal tails in cold dark matter cosmologies . MNRAS, 307:162–178.

Springel, V., Di Matteo, T., and Hernquist, L. (2005). Modelling feedback from stars and black holes in galaxy mergers . 361:776–794.

Springel, V., Yoshida, N., and White, S. D. M. (2001). GADGET: a code for collisionless and gasdynamical cosmological simulations . 6:79–117.

Springel, V. (2005). The cosmological simulation code GADGET-2. , 364:1105–1134.

Strigari, L. E., Bullock, J. S., Kaplinghat, M., Simon, J. D., Geha, M., Willman, B., & Walker, M. G. (2008a). , A common mass scale for satellite galaxies of the Milky Way . Nature, 454, 1096

Strigari, L. E., Koushiappas, S. M., Bullock, J. S., Kaplinghat, M., Simon, J. D., Geha, M., & Willman, B. (2008b), The most dark-matter-dominated galaxies: predicted gamma- ray signals from the faintest Milky Way dwarfs ApJ, 678, 614

Theis. C, Metz. Hensler. G, M, Kroupa. P, Jerjen. H. (2009). Discs of Satellites: The new dwarf spheroidals . Arxiv:0901;1658v1 [astro-ph.GA].

Tolstoy, E. (2003). Dwarf galaxies: Important clues to galaxy formation. Astrophysics & Space Science, 284:579–588. 148 Colisión de dos galaxias de disco

Toomre & Toomre. (1964). On the gravitational stability of a disk of stars . The Astrophysical Journal.

Toomre, A. and Toomre, J. (1972). Galactic Bridges and Tails . ApJ, 178:623–666.

Tully, R. B. (1982). “The Local Supercluster” . ApJ 257: 389–422.

Tully, R. B., Rizzi, L., Dolphin, A. E., Karachentsev, I. D., Karachentseva, V. E., Makarov, D. I.,Makarova, L., Sakai, S., and Shaya, E. J. (2006). Associations of Dwarf Galaxies . ApJ, 132:729–748.

Uson, J. M.; Boughn, S. P.; Kuhn, J. R. (1990). “The central galaxy in Abell 2029 – An old supergiant” . Science 250 (4980): 539–540. Bibcode:1990Sci...250..539U. doi:10.1126/science.250.4980.539. van de Kamp, P. (1964). Elements of Astromechanics . W. H. Freeman and Company. van den Bergh, S. (1999). The Local Group of galaxies . Journal of the Royal Astronomical Society of Canada, 93:175. van den Bergh, Sidney (2003). “History of the Local Group” . To be published in: “The Local Group as an Astrophysical Laboratory” (Cambridge University Press): 5042.

Weisz, D. R., Dolphin, A. E., Dalcanton, J. J., Skillman, E. D., Holtzman, J., Williams, B. F., Gilbert, K. M., Seth, A. C., Cole, A., Gogarten, S. M., Rosema, K., Karachentsev, I. D., McQuinn,K. B. W., and Zaritsky, D. (2011). How Typical Are the Local Group Dwarf Galaxies? ApJ,743:8.

Zwicky, F. (1956). Multiple Galaxies . Ergebnisse der exakten Naturwissenschaften, 29:344–385.

Zwicky, F. (1937). On the Masses of Nebulae and of Clusters of Nebulae, ApJ, 86, 217