Milky Way-like galaxy simulations in the Gaia era: disk large scale structures and baryonic content

[cat] En aquesta tesi, fent servir simulacions d'alta resolució i mitjançant diferents codis i tècniques de generació de condicions inicials hem pogut detectar que existeixen dos tipus diferents d'estructura espiral segons el seu perfil de rotació. Primer, en el cas de galàxies no barrades, les estructures roten a la mateixa velocitat que el disc, i segon, en el cas de galàxies barrades, aquestes estructures roten com un sòlid rígid. El segon resultat que presentem en aquesta tesi s'ha obtingut analitzant la cinemàtica estel·lar en simulacions de galàxies espirals. Per dur a terme aquest segon estudi hem utilitzat aproximacions analítiques i simulacions d'N-cossos pures i de partícules test. A partir d'aquests anàlisis hem pogut constatar que la desviació del vèrtex és un bon traçador de la posició de les estructures de densitat i també que el canvi de signe que pateix en creuar els pics de densitat de les espirals i les regions inter-braç ens permeten conèixer la posició de les principals ressonàncies, és a dir, la corrotació i la ressonància externa de Lindblad. Finalment, desprès de l'estudi exhaustiu dels codis i els processos físics relacionats amb la física de la component gasosa de les galàxies hem aconseguit obtenir una simulació cosmològica d'una galàxia molt semblant a la Via Làctia. Tot just hem començar a analitzar-la però ja hem obtingut resultats molt interessants tals com la presència de gran quantitat de gas calent a la regió de l'halo de matèria fosca o l'aparició de dues barres desalineades 90 graus a la zona del disc galàctic, una de jove, procedent del disc i una de vella, que és el fòssil d'una fusió de dues galàxies el·líptiques a un desplaçament cap el vermell de 3.[eng] Simulations have shown to be one of the best tools to study properties of galactic large scale structures and their effects on to the local kinematics of stars. The aims of this thesis are: i) To obtain realistic N-body models that allow us to analyse kinematics, dynamics and internal structure of non-axisymmetric components in galaxies, ii) to learn how to control numerical effects and also how to distinguish them from the proper physical ones, iii) to find observable parameters from stellar kinematics that can give us information about formation, evolution and nature of the large scale structures in galaxies and, iv) to test which of such methods can be used to distinguish among spiral arm natures in real galaxies. In this thesis, using high resolution simulations obtained with different codes and initial condition techniques, we find that exist two different behaviours for the rotation frequency of transient spiral arms like structures. Whereas unbarred disks present spiral arms nearly corotating with disk particles, strong barred models (bulged or bulge-less) quickly develop a bar-spiral structure dominant in density, with a pattern speed almost constant in radius (Roca-Fàbrega et al. 2013). Preliminary results also indicate that particles in barred models move inside the spiral structures. A second result we present in the thesis has been obtained mapping the kinematics of stars in the simulated galaxy disks with spiral arms using the velocity ellipsoid vertex deviation (lv). For this study we have used both test particle and high resolution N-body simulations. What we have found is that for all barred models, spiral arms rotate closely to a rigid body manner and there the vertex deviation values correlate with the density peaks position bounded by overdense and underdense regions. However, the most interesting result is that In such cases, vertex deviation sign changes from negative to positive when crossing the spiral arms toward disk rotation, in regions where the spiral arms are in between corotation (CR) and the Outer Lindblad Resonance (OLR). By contrast, when the arm sections are inside the CR and outside the OLR, lv changes from negative to positive. We propose that measurements of the vertex deviations pattern can be used to trace the position of the main resonances of the spiral arms (Roca-Fàbrega et al. 2014). Finally we present a new cosmological Milky Way like galaxy simulation that includes both the collisionless N-body and also the gas components. This simulation has been obtained using the adaptive mesh refinement (AMR) N-body code ART (Kravtsov et al 1999) plus the hydrodynamics and physical processes presented by Kravtsov et al 2003. The MW like system has been evolved inside a 28 Mpc cosmological box with a spatial resolution of 109 pc. At z=0 the system has an Mvir = 7.33·10^11Msun. We have observed how a well defined disk is formed inside the dark matter halo and the overall amount of gas and stars is comparable with MW observations. Several non-axisymmetric structures arise out of the disk. We have also observed that a huge reservoir of hot gas is present at large distances from the disk. Gas column density, emission and dispersion measures have been computed from inside the simulated disk at a position of 8 kpc from the center and in several different directions. Our preliminary results reveal that the distribution of hot gas is non-isotropic according with observations. After a careful analysis we confirm that due to the anisotropy in the gas distribution more than 50 random observations of different sky regions are needed to recover the real distribution of hot gas in the galactic halo

Distinguishing Between Photoionized and Collisionally Ionized Gas in the Circumgalactic Medium

Most studies of highly ionized plasmas have historically assumed ions are either in photoionization equilibrium, PIE, or collisional ionization equilibrium, CIE, sometimes including multiple phases with different relevant mechanisms. Simulation analysis packages, on the other hand, tend to use precomputed ion fraction tables which include both mechanisms, among others. Focusing on the low-density, high temperature phase space likely to be most relevant in the circumgalactic medium, in this work we show that most ions can be classified as 'PI' or 'CI' on an ion-by-ion basis. This means that for a cloud at a particular point in phase space, some ions will be created primarily by PI and others by CI, with other mechanisms playing only very minor roles. Specifically, we show that ions are generally CI if the thermal energy per particle is greater than $\sim6$\% of their ionization energy, and PI otherwise. We analyse the accuracy of this ansatz compared to usual PIE/CIE calculations, and show the surprisingly minor dependence of this conclusion on redshift and ionizing background.Comment: 12 pages, submitted to MNRA

Kinematics of symmetric Galactic longitudes to probe the spiral arms of the Milky Way with Gaia

Aims. We model the effects of the spiral arms of the Milky Way on the disk stellar kinematics in the Gaia observable space. We also estimate the Gaia capabilities of detecting the predicted signatures. Methods. We use both controlled orbital integrations in analytic potentials and self-consistent simulations. We introduce a new strategy to investigate the effects of spiral arms, which consists of comparing the stellar kinematics of symmetric Galactic longitudes (+l and −l), in particular the median transverse velocity as determined from parallaxes and proper motions. This approach does not require the assumption of an axisymmetric model because it involves an internal comparison of the data. Results. The typical differences between the transverse velocity in symmetric longitudes in the models are of the order of ~2 km s-1, but can be larger than 10 km s-1 for certain longitudes and distances. The longitudes close to the Galactic centre and to the anti-centre are those with larger and smaller differences, respectively. The differences between the kinematics for +l and −l show clear trends that depend strongly on the properties of spiral arms. Thus, this method can be used to quantify the importance of the effects of spiral arms on the orbits of stars in the different regions of the disk, and to constrain the location of the arms, main resonances and, thus, pattern speed. Moreover, the method allows us to test different origin scenarios of spiral arms and the dynamical nature of the spiral structure (e.g. grand design versus transient multiple arms). We estimate the number of stars of each spectral type that Gaia will observe in certain representative Galactic longitudes, their characteristic errors in distance and transverse velocity, and the error in computing the median velocity as a function of distance. We will be able to measure the median transverse velocity exclusively with Gaia data, with precision smaller than ~1 km s-1 up to distances of ~4-6 kpc for certain giant stars, and up to ~2-4 kpc and better kinematic precision (≲0.5 km s-1) for certain sub-giants and dwarfs. These are enough to measure the typical signatures seen in the models. Conclusions. The Gaia catalogue will allow us to use the presented approach successfully and improve significantly upon current studies of the dynamics of the spiral arms of our Galaxy. We also show that a similar strategy can be used with line-of-sight velocities, which could be applied to Gaia data and to upcoming spectroscopic surveys

Gaia DR2 reveals a star formation burst in the disc 2-3 Gyr ago

We use Gaia DR2 magnitudes, colours and parallaxes for stars with G1.53 Msun and α2≈1.3 for the mass range between 0,5 and 1,53 Msun. This is the first time that we consider a non-parametric SFH for the thin disc in the Besancon Galaxy Model. This new step, together with the capabilities of the Gaia DR2 parallaxes to break degeneracies between different stellar populations, allow us to better constrain the SFH and the IMF

OVI Traces Photoionized Streams With Collisionally Ionized Boundaries in Cosmological Simulations of z∼1z \sim 1 Massive Galaxies

We analyse the distribution and origin of OVI in the Circumgalactic Medium (CGM) of dark-matter haloes of $\sim 10^{12}$ M$_\odot$ at $z\sim1$ in the VELA cosmological zoom-in simulations. We find that the OVI in the inflowing cold streams is primarily photoionized, while in the bulk volume it is primarily collisionally ionized. The photoionized component dominates the observed column density at large impact parameters ($\gtrsim 0.3 R_{\rm vir}$), while the collisionally ionized component dominates closer in. We find that most of the collisional OVI, by mass, resides in the relatively thin boundaries of the photoionized streams. We discuss how the results are in agreement with analytic predictions of stream and boundary properties, and their compatibility with observations. This allows us to predict the profiles of OVI and other ions in future CGM observations and provides a toy model for interpreting them.Comment: 22 pages, 16 figures, submitted to MNRA

A novel method to bracket the corotation radius in galaxy discs:vertex deviation maps

We map the kinematics of stars in simulated galaxy disks with spiral arms using the velocity ellipsoid vertex deviation (l$_v$). We use test particle simulations, and for the first time, fully self-consistent high resolution N-body models. We compare our maps with the Tight Winding Approximation model analytical predictions. We see that for all barred models spiral arms rotate closely to a rigid body manner and the vertex deviation values correlate with the density peaks position bounded by overdense and underdense regions. In such cases, vertex deviation sign changes from negative to positive when crossing the spiral arms in the direction of disk rotation, in regions where the spiral arms are in between corotation (CR) and the Outer Lindblad Resonance (OLR). By contrast, when the arm sections are inside the CR and outside the OLR, l$_v$ changes from negative to positive.We propose that measurements of the vertex deviations pattern can be used to trace the position of the main resonances of the spiral arms. We propose that this technique might exploit future data from Gaia and APOGEE surveys. For unbarred N-body simulations with spiral arms corotating with disk material at all radii, our analysis suggests that no clear correlation exists between l$_v$ and density structures

GARROTXA Cosmological Simulations of Milky Way-sized Galaxies: General Properties, Hot-gas Distribution, and Missing Baryons

We introduce a new set of simulations of Milky Way (MW)-sized galaxies using the AMR code ART + hydrodynamics in a Λ cold dark matter cosmogony. The simulation series is called GARROTXA and it follows the formation of a halo/galaxy from z = 60 to z = 0. The final virial mass of the system is ¿7.4 × 1011 M ⊙. Our results are as follows. (a) Contrary to many previous studies, the circular velocity curve shows no central peak and overall agrees with recent MW observations. (b) Other quantities, such as M\_\ast (6 × 1010 M ⊙) and R d (2.56 kpc), fall well inside the observational MW range. (c) We measure the disk-to-total ratio kinematically and find that D/T = 0.42. (d) The cold-gas fraction and star formation rate at z = 0, on the other hand, fall short of the values estimated for the MW. As a first scientific exploitation of the simulation series, we study the spatial distribution of hot X-ray luminous gas. We have found that most of this X-ray emitting gas is in a halo-like distribution accounting for an important fraction but not all of the missing baryons. An important amount of hot gas is also present in filaments. In all our models there is not a massive disk-like hot-gas distribution dominating the column density. Our analysis of hot-gas mock observations reveals that the homogeneity assumption leads to an overestimation of the total mass by factors of 3-5 or to an underestimation by factors of 0.7-0.1, depending on the used observational method. Finally, we confirm a clear correlation between the total hot-gas mass and the dark matter halo mass of galactic systems

The intricate Galaxy disk: velocity asymmetries in Gaia-TGAS

We use Gaia-TGAS data to compare the transverse velocities in Galactic longitude (coming from proper motions and parallaxes) in the Milky Way disk for negative and positive longitudes as a function of distance. The transverse velocities are strongly asymmetric and deviate significantly from the expectations for an axisymmetric galaxy. The value and sign of the asymmetry changes at spatial scales of several tens of degrees in Galactic longitude and about 0.5 kpc in distance. The asymmetry is statistically significant at 95% confidence level for 57% of the region probed, which extends up to 1.2 kpc. A percentage of 24% of the region shows absolute differences at this confidence level larger than 5 km s-1 and 7% larger than 10 km s-1. The asymmetry pattern shows mild variations in the vertical direction and with stellar type. A first qualitative comparison with spiral arm models indicates that the arms are probably not the main source of the asymmetry. We briefly discuss alternative origins. This is the first time that global all-sky asymmetries are detected in the Milky Way kinematics beyond the local neighbourhood and with a purely astrometric sample