Substructure in the phase-space of the Galaxy with Gaia

Abstract

The Gaia mission and its extensive catalogue of stars has marked the beginning of the Golden Age of Galactic dynamics. We have now access to the kinematics of more than a billion stars with which we can start deciphering the history of our Galaxy. Before Gaia, we knew that the phase-space of the Milky Way (MW) contained a wealth of substructure created by the different dynamical processes and components of the Galaxy, but we did not have enough kinematical data with which to characterise it with precision outside the Solar neighbourhood (SN).

In this thesis we set out to detect and describe the kinematic substructure of the MW, not only near the Sun but to the very limit of the Gaia sphere. Our goal is to characterise the phase-space features of our Galaxy, introducing novel empirical measurements that can help us understand the structure of the MW and its evolution, together with the dynamics of its different components (bar, spiral arms, halo, satellites), and the interactions among them over time. We will focus on the disc and the structures of the halo, in particular, the Sagittarius (Sgr) dwarf spheroidal, a galaxy undergoing full tidal disruption and most likely a significant perturber in the recent history of our Galaxy.

Our methodology is centred on the analysis of the Gaia DR2 data mostly using the wavelet transformation (WT) that is a robust statistical tool to reveal the substructure present in the velocity maps. We apply it on the plane of Galactocentric radial velocity (vR) against rotational velocity (vφ), venturing into regions beyond the SN, but also, and for the first time, to the proper motion histograms of distant stars all around the sky.

As a result, we find that the velocity field of the MW disc is highly structured not only at the SN but, at least, up to distances of ∼3kpc. We have characterised the shape of each structure and described their changes with position in the Galaxy. We report that the vφ of each structure decreases with Galactocentric radius in a different way, which allows us to associate tentatively each structure with a different dynamical mechanism, namely, the phase-mixing process most likely initiated by the last pericentre passage of Sgr, and the resonances with the non-axisymmetric components of the MW (bar and/or spiral arms).

Additionally, we detect two arch-like structures towards the anticentre corresponding to Monoceros and Anticentre stream (ACS). Thanks to our method, we are able to observe their morphology sharper than ever and produce a kinematic selection of member stars. From this, we obtain a precise characterisation of Monoceros and ACS that can be now used to compare quantitatively the models proposed to explain their origin. Based on our preliminary analysis, we favour the scenario in which they were caused by the repeated perturbations induced by Sgr, as opposed to being disrupted satellites as some authors propose.

With respect to Sgr, we have detected the kinematic signature of its stream across the whole sky, and measured the proper motion of this system along a large portion of the tidal tails, which was not possible before Gaia. We have shown also that the predictions of the models do not match our observations, and thus require revision. The resulting sample of candidate stars represents a vast catalogue of stars of different stellar types that can be used to model the orbit of Sgr, as well as their distribution along the stream and its link to the star formation history.

Finally, we complement this sample with the largest list of RR Lyrae in Sgr. In doing so, we have obtained a precise characterisation of the 3D distribution of its tidal debris as well as its tangential velocities, mostly unexplored until now as these require distances and proper motions simultaneously.

All together, we have provided a detailed characterisation of the kinematic substructure present throughout the MW that we can now use to study and model the evolution of our Galaxy and its components in a global, transverse manner which takes into account the coupling between the Sgr stream and the disc.

Des del moment en que la Via Làctia (VL) va néixer, diferents mecanismes dinàmics han anat interactuant i modificant la distribució d’estels a l’espai de fases. Aquesta tesi està enfocada a la detecció i caracterització de la subestructura cinemàtica present a les dades del satèl·lit Gaia, en diferents projeccions de l’espai de fases. De l’estudi del pla de velocitats Vr-Vphi hem obtingut la llista més completa de grups mòbils (grups d’estrelles amb cinemàtica coherent), amb la que hem estudiat com la seva localització dins del pla varia en funció d’on estiguem de la Galaxia. Els gradients en funció del radi i azimuth Galactocèntrics que hem mesurat ens han permès, per primer cop, fer una associació temptativa de cada estructura amb un mecanisme dinàmic diferent: ressonàncies amb la barra/braços espirals o phase-mixing, probablement induït per la interacció amb la galaxia nana de Sagitari (Sgr). En estudiar les parts llunyanes del disc, hem detectat les estructures de Monoceros i ACS a partir solament de la seva cinemàtica, a diferència dels estudis anteriors. Aquestes estructures, que gràcies als nostres resultats s’han pogut veure de forma més clara que mai, estan probablement causades també per Sgr i, amb les dades que aportem, ara serà possible testejar de forma acurada aquesta hipòtesi (contrastant-la amb d’altres que existeixen a la literatura). Finalment, aquesta mateixa metodologia ens ha permès també identificar la signatura cinemàtica del corrent estelar de Sgr de forma empírica i homogènia en tot el cel. La mostra resultant és la més gran existent actualment i aporta les dues components de la velocitat que fins ara havien faltat: moviments propis al llarg de tota la estructura. D’altra banda, aplicant un metodologia diferent, hem compilat la mostra més gran de RR Lyrae pertanyents al corrent de Sgr. Per aquestes estrelles es poden obtenir distàncies molt precises i, combinades a l’altre mostra de Sgr, permetran modelar de forma detallada l'òrbita passada d’aquesta galàxia nana. La combinació de tots els nostres resultats connecta la cinemàtica del disc amb la de l’halo i obra la porta a caracteritzar, per primer cop, la història de la VL de forma híbrida.

La astronomía Galáctica comenzó su época dorada el día 25 de abril de 2018, coincidiendo con los inicios de esta tesis, al hacerse públicos los datos de los primeros 22 meses de observación del satélite Gaia de la Agencia Espacial Europea. Dichos datos contienen posiciones y movimientos propios para más de 1.3 miles de millones de estrellas de nuestra Galaxia, la Vía Láctea (VL). Para una porción de estas, las más brillantes y cercanas, Gaia también proporciona distancias precisas y velocidades radiales, lo cual supone la muestra 6D (posiciones y velocidades) más grande hasta la fecha con alrededor de 7 millones de estrellas. El objetivo de esta tesis es usar todo el potencial de estos datos para estudiar en detalle la sub-estructura cinemática de las estrellas que forman la VL, tanto en el disco como en el halo y sus satélites. Se cree que dicha sub-estructura está inducida en su mayoría por elementos internos como la barra o los brazos espirales, o bien por interacciones con galaxias enanas que, al ser devoradas por nuestra Galaxia, perturban gravitacionalmente el disco. En este sentido, debemos fijarnos en la galaxia enana de Sagitario (Sgr), actualmente en pleno proceso de disrupción, claramente visible en forma de corriente estelar (producido por las estrellas que ha ido perdiendo con el tiempo) que da, al menos, una vuelta entera en el cielo. Estudios recientes demuestran que esta galaxia satélite fue lo suficientemente masiva en el pasado como para perturbar significativamente la VL.

Primero, usamos la muestra 6D para detectar la sub-estructura cinemática presente en los mapas de velocidad, tanto cerca del Sol, como en regiones inexploradas del disco de la VL. El resultado es la lista más completa de estructuras cinemáticas y su precisa caracterización morfológica. Al poder explorar regiones alejadas del Sol, hemos sido capaces de seguir las estructuras a medida que nos movíamos por la Galaxia y medir sus cambios de localización en el plano de las velocidades. Esto ha llevado al descubrimiento de que las estructuras cambian con el radio Galactocéntrico de formas distintas, algunas manteniendo un momento angular común (lo cual podríamos asociar con resonancias de la barra/brazos espirales) y otras compartiendo energías/frecuencias orbitales (posiblemente vinculadas con perturbaciones inducidas por la galaxia enana de Sgr).

Luego expandimos nuestra búsqueda a las partes más lejanas de la Galaxia usando un método novedoso que permite procesar de forma rápida toda la muestra de Gaia al completo. Para ello, buscamos sobre-densidades en los histogramas de movimientos propios (ya no usamos la velocidad radial) de todas las estrellas lejanas que encontramos en cada uno de los HEALpix de nivel 5 en que hemos dividido el cielo. Como resultado, hemos detectado dos estructuras en forma de arco en el hemisferio norte Galáctico, en la dirección del anticentro. Dichas estructuras corresponden a Monoceros y ACS, siendo esta la primera vez que han sido detectadas únicamente por su cinemática y sin limitarse a un tipo de estrellas concreto, como había sido siempre el caso hasta ahora. Tras caracterizar su morfología, cinemática y la distribución de sus estrellas en los diagramas de color-magnitud, concluimos que estas estructuras parecen estar asociadas al disco tal como algunos estudios apuntan, en lugar de al halo como otros trabajos defendían. Así pues, la

causa más probable para explicar sus propiedades observadas es que Sgr, al interactuar gravitatoriamente con la VL, impulsó las estrellas que habitan en el borde del disco separándolas así de este.

Aplicando la misma metodología usada para detectar Monoceros y ACS, hemos sido capaces de producir la muestra más grande existente de estrellas pertenecientes al corriente estelar de Sgr. Dicha muestra aporta movimientos propios a lo largo de la mayor parte de las colas de marea, lo cual hasta ahora no había sido posible al nivel de precisión actual, y que permite conocer con precisión los cambios en la cinemática que la corriente experimenta a lo largo de su órbita. Igual que antes, nuestra muestra no se restringe a una población estelar concreta, lo cual nos posibilitará estudiar con detalle la composición a lo largo de las colas de marea. Esto, unido a la precisa caracterización de su cinemática, permitirá modelar con gran precisión la historia de formación de Sgr.

Finalmente, expandimos la muestra de estrellas de Sgr con la lista más grande de RR Lyrae hasta la fecha. A pesar de ser una muestra pequeña en comparación a la anterior, debido básicamente a la baja proporción de este tipo de estrellas, el hecho de ser candelas estándar permite calcular su distancia con gran precisión. De esta manera, somos capaces de determinar la distribución espacial en 3D de las estrellas que forman la corriente y, al añadir los movimientos propios, obtener por primera vez las velocidades tangenciales a lo largo de las colas de Sgr. Con esta muestra y la anterior, hemos demostrado que la cinemática de modelos ampliamente usados hasta ahora no se ajusta a las observaciones, lo cual nos lleva a la necesaria revisión del potencial gravitatorio de la VL. Los nuevos datos aportados en esta tesi servirán para ajustar de forma más precisa que nunca la órbita de Sgr al poder contar de finalmente con los movimientos propios, a la vez que distancias precisas, para una gran cantidad de estrellas.

Los tres aspectos estudiados en esta tesis (sub-estructura cinemática en el disco, Monoceros-ACS, y la órbita de Sgr) están, muy probablemente, estrechamente ligados entre sí. El potencial gravitatorio de la VL da forma y moldea la cinemática de cada uno de ellos individualmente pero, al mismo tiempo, es Sgr a través de su órbita la responsable, probablemente, de la aparición de algunas de las estructuras observadas. Se crea así una conexión entre los tres elementos que, gracias a los datos aportados en esta tesis, permitirá realizar por primera vez un ajuste híbrido a los parámetros que definen las características de nuestra Galaxia así como su evolución en el tiempo.