On the stellar halo metallicity profile of Milky Way-like galaxies in the Auriga simulations

A recent observational study of haloes of nearby Milky Way-like galaxies shows that only half (four out of eight) of the current sample exhibits strong negative metallicity ([Fe/H]) gradients. This is at odds with predictions from hydrodynamical simulations where such gradients are ubiquitous. In this Letter, we use high-resolution cosmological hydrodynamical simulations to study the [Fe/H] distribution of galactic haloes. We find that kinematically selected stellar haloes, including both in situ and accreted particles, have an oblate [Fe/H] distribution. Spherical [Fe/H] radial profiles show strong negative gradients within 100 kpc, in agreement with previous numerical results. However, the projected median [Fe/H] profiles along the galactic disc minor axis, typically obtained in observations, are significantly flatter. The median [Fe/H] values at a given radius are larger for the spherical profiles than for the minor axis profiles by as much as 0.4 dex within the inner 50 kpc. Similar results are obtained if only the accreted stellar component is considered indicating that the differences between spherical and minor axis profiles are not purely driven by heated disc star particles formed in situ. Our study highlights the importance of performing careful comparisons between models and observations of halo [Fe/H] distributions

The chemo-dynamical groups of Galactic globular clusters

Galaxie

On the relevance of chaos for halo stars in the solar neighbourhood II

In a previous paper based on dark matter only simulations we show that, in the approximation of an analytic and static potential describing the strongly triaxial and cuspy shape of Milky Way-sized haloes, diffusion due to chaotic mixing in the neighbourhood of the Sun does not efficiently erase phase space signatures of past accretion events. In this second paper we further explore the effect of chaotic mixing using multicomponent Galactic potential models and solar neighbourhood-like volumes extracted from fully cosmological hydrodynamic simulations, thus naturally accounting for the gravitational potential associated with baryonic components, such as the bulge and disc. Despite the strong change in the global Galactic potentials with respect to those obtained in dark matter only simulations, our results confirm that a large fraction of halo particles evolving on chaotic orbits exhibit their chaotic behaviour after periods of time significantly larger than a Hubble time. In addition, significant diffusion in phase space is not observed on those particles that do exhibit chaotic behaviour within a Hubble time

Aurigaia: mock Gaia DR2 stellar catalogues from the Auriga cosmological simulations

We present and analyse mock stellar catalogues that match the selection criteria and observables (including uncertainties) of the Gaia satellite data release 2 (DR2). The source are six cosmological high-resolution magneto-hydrodynamic ΛCDM zoom simulations of the formation of Milky Way analogues from the AURIGA project. Mock data are provided for stars with V 20 deg. The mock catalogues are made using two different methods: the public SNAPDRAGONS code, and a method based on that of Lowing et al. (2015) that preserves the phase-space distribution of the model stars. These publicly available catalogues contain five-parameter astrometry, radial velocities, multiband photometry, stellar parameters, dust extinction values, and uncertainties in all these quantities. In addition, we provide the gravitational potential and information on the origin of each star. By way of demonstration, we apply the mock catalogues to analyses of the young stellar disc and the stellar halo. We show that (i) the young outer stellar disc exhibits a flared distribution that is detectable in the height and vertical velocity distribution of A - and B -dwarf stars up to radii of ∼15 kpc, and (ii) the spin of the stellar halo out to 100 kpc can be accurately measured with Gaia DR2 RR Lyrae stars. These catalogues are well suited for comparisons with observations and should help to (i) develop and test analysis methods for the Gaia DR2 data, (ii) gauge the limitations and biases of the data, and (iii) interpret the data in the light of theoretical predictions from realistic ab initio simulations of galaxy formation in the ΛCDM cosmological model

Stellar dynamics around transient co-rotating spiral arms

Spiral density wave theory attempts to describe the spiral pattern in spiral galaxies in terms of a long-lived wave structure with a constant pattern speed in order to avoid the winding dilemma. The pattern is consequently a rigidly rotating, long-lived feature. We run an N-body/SPH simulation of a Milky Way-sized barred disk, and find that the spiral arms are transient features whose pattern speeds decrease with radius, in such a way that the pattern speed is almost equal to the rotation curve of the galaxy. We trace particle motion around the spiral arms. We show that particles from behind and in front of the spiral arm are drawn towards and join the arm. Particles move along the arm in the radial direction and we find a clear trend that they migrate toward the outer (inner) radii on the trailing (leading) side of the arm. Our simulations demonstrate that tat all radii where there is a co-rotating spiral arm the particles continue to be accelerated (decelerated) by the spiral arm for long periods, which leads to strong migration

Properties of HI discs in the Auriga cosmological simulations

We analyse the properties of the H i gas distribution in the Auriga project, a set of magnetohydrodynamic cosmological simulations performed with the moving-mesh code arepo and a physics model for galaxy formation that succeeds in forming realistic late-type galaxies in the 30 Milky Way-sized haloes simulated in this project. We use a simple approach to estimate the neutral hydrogen fraction in our simulation set, which treats low-density and star-forming gas separately, and we explore two different prescriptions to subtract the contribution of molecular hydrogen from the total H i content. The H i gas in the vast majority of the systems forms extended discs although more disturbed morphologies are present. Notwithstanding the general good agreement with observed H i properties – such as radial profiles and the mass–diameter relation – the Auriga galaxies are systematically larger and more gas-rich than typical nearby galaxies. Interestingly, the amount of H i gas outside the disc plane correlates with the star formation rate, consistent with a picture where most of this extra-planar H i gas originates from a fountain-like flow. Our findings are robust with respect to the different assumptions adopted for computing the molecular hydrogen fraction and do not vary significantly over a wide range of numerical resolution. The H i modelling introduced in this paper can be used in future work to build artificial interferometric H i data cubes, allowing an even closer comparison of the gas dynamics in simulated galaxies with observations

The orbital phase space of contracted dark matter haloes

We study the orbital phase space of dark matter (DM) haloes in the AURIGA suite of cosmological hydrodynamics simulations of Milky Way (MW) analogues. We characterize haloes by their spherical action distribution, F (Jr, L), a function of the specific angular momentum, L, and the radial action, Jr, of the DM particles. By comparing DM-only and hydrodynamical simulations of the same haloes, we investigate the contraction of DM haloes caused by the accumulation of baryons at the centre. We find a small systematic suppression of the radial action in the DM haloes of the hydrodynamical simulations, suggesting that the commonly used adiabatic contraction approximation can result in an underestimate of the density by 3c 8 per cent. We apply an iterative algorithm to contract the AURIGA DM haloes given a baryon density profile and halo mass, recovering the true contracted DM profiles with an accuracy of 3c 15 per cent, that reflects halo-to-halo variation. Using this algorithm, we infer the total mass profile of the MW\u2019s contracted DM halo. We derive updated values for the key astrophysical inputs to DM direct detection experiments: the DM density and velocity distribution in the Solar neighbourhood

Vertical disc heating in Milky Way-sized galaxies in a cosmological context

Vertically extended, high velocity dispersion stellar distributions appear to be a ubiquitous feature of disc galaxies, and both internal and external mechanisms have been proposed to be the major driver of their formation. However, it is unclear to what extent each mechanism can generate such a distribution, which is likely to depend on the assembly history of the galaxy. To this end, we perform 16 high-resolution cosmological-zoom simulations of Milky Way-sized galaxies using the state-of-the-art cosmological magnetohydrodynamical code AREPO, and analyse the evolution of the vertical kinematics of the stellar disc in connection with various heating mechanisms. We find that the bar is the dominant heating mechanism in most cases, whereas spiral arms, radial migration and adiabatic heating from mid-plane density growth are all subdominant. The strongest source, though less prevalent than bars, originates from external perturbations from satellites/subhaloes of masses log10(M/M⊙) ≳ 10. However, in many simulations the orbits of newborn star particles become cooler with time, such that they dominate the shape of the age–velocity dispersion relation and overall vertical disc structure unless a strong external perturbation takes place

The dark matter component of the Gaia radially anisotropic substructure

We study the properties of the dark matter component of the radially anisotropic stellar population recently identified in the Gaia data, using magneto-hydrodynamical simulations of Milky Way-like halos from the Auriga project. We identify 10 simulated galaxies that approximately match the rotation curve and stellar mass of the Milky Way. Four of these have an anisotropic stellar population reminiscent of the Gaia structure. We find an anti-correlation between the dark matter mass fraction of this population in the Solar neighbourhood and its orbital anisotropy. We estimate the local dark matter density and velocity distribution for halos with and without the anisotropic stellar population, and use them to simulate the signals expected in future xenon and germanium direct detection experiments. We find that a generalized Maxwellian distribution fits the dark matter halo integrals of the Milky Way-like halos containing the radially anisotropic stellar population. For dark matter particle masses below approximately 10 GeV, direct detection exclusion limits for the simulated halos with the anisotropic stellar population show a mild shift towards smaller masses compared to the commonly adopted Standard Halo Model