No Arabic abstract
We study the gas distribution in the Milky Way and Andromeda using a constrained cosmological simulation of the Local Group (LG) within the context of the CLUES (Constrained Local UniversE Simulations) project. We analyse the properties of gas in the simulated galaxies at $z=0$ for three different phases: `cold, `hot and HI, and compare our results with observations. The amount of material in the hot halo ($M_{hot}approx 4-5times10^{10},$M$_{odot}$), and the cold ($M_{cold}(rlesssim10,$kpc$)approx10^{8},$M$_{odot}$) and HI ($M_{HI}(rlesssim50,$kpc$)approx 3-4times10^8,$M$_{odot}$) components display a reasonable agreement with observations. We also compute the accretion/ejection rates together with the HI (radial and all-sky) covering fractions. The integrated HI accretion rate within $r=50,$kpc gives $sim$$0.2-0.3,$M$_{odot},$yr$^{-1}$, i.e. close to that obtained from high-velocity clouds in the Milky Way. We find that the global accretion rate is dominated by hot material, although ionized gas with $Tlesssim10^5,$K can contribute significantly too. The $net$ accretion rates of $all$ material at the virial radii are $6-8,$M$_{odot},$yr$^{-1}$. At $z=0$, we find a significant gas excess between the two galaxies, as compared to any other direction, resulting from the overlap of their gaseous haloes. In our simulation, the gas excess first occurs at $zsim1$, as a consequence of the kinematical evolution of the LG.
We apply a semi-analytic galaxy formation model to two high resolution cosmological N-body simulations to investigate analogues of the Milky Way system. We select these according to observed properties of the Milky Way rather than by halo mass as in most previous work. For disk-dominated central galaxies with stellar mass (5--7) x 10d10Msun, the median host halo mass is 1.4 x 10d12Msun, with 1 sigma dispersion in the range [0.86, 3.1] x 10d12Msun, consistent with dynamical measurements of the Milky Way halo mass. For any given halo mass, the probability of hosting a Milky Way system is low, with a maximum of ~20% in haloes of mass ~10d12Msun. The model reproduces the V-band luminosity function and radial profile of the bright (MV < -9) Milky Way satellites. Galaxy formation in low mass haloes is found to be highly stochastic, resulting in an extremely large scatter in the relation between MV (or stellar mass) for satellites and the depth of the subhalo potential well in which they live, as measured by the maximum of the rotation curve, Vmax. We conclude that the too big to fail problem is an artifact of selecting satellites in N-body simulations according to subhalo properties: in 10% of cases we find that three or fewer of the brightest (or most massive) satellites have Vmax > 30 km/s. Our model predicts that around half of the dark matter subhaloes with Vmax > 20 km/s host satellites fainter than MV = -9 and so may be missing from existing surveys.
We study the two main constituent galaxies of a constrained simulation of the Local Group as candidates for the Milky Way (MW) and Andromeda (M31). We focus on the formation of the stellar discs and its relation to the formation of the group as a rich system with two massive galaxies, and investigate the effects of mergers and accretion as drivers of morphological transformations. We use a state-of-the-art hydrodynamical code which includes star formation, feedback and chemical enrichment to carry out our study. We run two simulations, where we include or neglect the effects of radiation pressure from stars, to investigate the impact of this process on the morphologies and star formation rates of the simulated galaxies. We find that the simulated M31 and MW have different formation histories, even though both inhabit, at z=0, the same environment. These differences directly translate into and explain variations in their star formation rates, in-situ fractions and final morphologies. The M31 candidate has an active merger history, as a result of which its stellar disc is unable to survive unaffected until the present time. In contrast, the MW candidate has a smoother history with no major mergers at late times, and forms a disc that grows steadily; at z=0 the simulated MW has an extended, rotationally-supported disc which is dominant over the bulge. Our two feedback implementations predict similar evolution of the galaxies and their discs, although some variations are detected, the most important of which is the formation time of the discs: in the model with weaker/stronger feedback the discs form earlier/later. In summary, by comparing the formation histories of the two galaxies, we conclude that the particular merger/accretion history of a galaxy rather than its environment at the LG-scales is the main driver of the formation and subsequent growth or destruction of galaxy discs.
With Gaia Data Release 2, the astronomical community is entering a new era of multidimensional surveys of the Milky Way. This new phase-space view of our Galaxy demands new tools for comparing observations to simulations of Milky-Way-mass galaxies in a cosmological context, to test the physics of both dark matter and galaxy formation. We present ananke, a framework for generating synthetic phase-space surveys from high-resolution baryonic simulations, and use it to generate a suite of synthetic surveys resembling Gaia DR2 in data structure, magnitude limits, and observational errors. We use three cosmological simulations of Milky-Way-mass galaxies from the Latte suite of the Feedback In Realistic Environments (FIRE) project, which feature self-consistent clustering of star formation in dense molecular clouds and thin stellar/gaseous disks in live cosmological halos with satellite dwarf galaxies and stellar halos. We select three solar viewpoints from each simulation to generate nine synthetic Gaia-like surveys. We sample synthetic stars by assuming each star particle (of mass 7070 $M_{odot}$) represents a single stellar population. At each viewpoint, we compute dust extinction from the simulated gas metallicity distribution and apply a simple error model to produce a synthetic Gaia-like survey that includes both observational properties and a pointer to the generating star particle. We provide the complete simulation snapshot at $z = 0$ for each simulated galaxy. We describe data access points, the data model, and plans for future upgrades. These synthetic surveys provide a tool for the scientific community to test analysis methods and interpret Gaia data.
We examine the prevalence, longevity, and causes of planes of satellite dwarf galaxies, as observed in the Local Group. We use 14 Milky Way/Andromeda-(MW/M31) mass host galaxies from the FIRE-2 simulations. We select the 14 most massive satellites by stellar mass within 300 kpc of each host and correct for incompleteness from the foreground galactic disc when comparing to the MW. We find that MW-like planes as spatially thin and/or kinematically coherent as observed are uncommon, but they do exist in our simulations. Spatially thin planes occur in 1-2 per cent of snapshots during $z=0-0.2$, and kinematically coherent planes occur in 5 per cent of snapshots. These planes are generally transient, surviving for less than 500 Myr. However, if we select hosts with an LMC-like satellite near first pericentre, the fraction of snapshots with MW-like planes increases dramatically to 7-16 per cent, with lifetimes of 0.7-1 Gyr, likely because of group accretion of satellites. We find that M31s satellite distribution is much more common: M31s satellites lie within about 1 sigma of the simulation median for every plane metric we consider. We find no significant difference in average satellite planarity for isolated hosts versus hosts in LG-like pairs. Baryonic and dark matter-only simulations exhibit similar levels of planarity, even though baryonic subhaloes are less centrally concentrated within their host haloes. We conclude that planes of satellites are not a strong challenge to LCDM cosmology.
We present a direct comparison of the Pan-Andromeda Archaeological Survey (PAndAS) observations of the stellar halo of M31 with the stellar halos of 6 galaxies from the Auriga simulations. We process the simulated halos through the Auriga2PAndAS pipeline and create PAndAS-like mocks that fold in all observational limitations of the survey data (foreground contamination from the Milky Way stars, incompleteness of the stellar catalogues, photometric uncertainties, etc). This allows us to study the survey data and the mocks in the same way and generate directly comparable density maps and radial density profiles. We show that the simulations are overall compatible with the observations. Nevertheless, some systematic differences exist, such as a preponderance for metal-rich stars in the mocks. While these differences could suggest that M31 had a different accretion history or has a different mass compared to the simulated systems, it is more likely a consequence of an under-quenching of the star formation history of galaxies, related to the resolution of the Auriga simulations. The direct comparison enabled by our approach offers avenues to improve our understanding of galaxy formation as they can help pinpoint the observable differences between observations and simulations. Ideally, this approach will be further developed through an application to other stellar halo simulations. To facilitate this step, we release the pipeline to generate the mocks, along with the six mocks presented and used in this contribution.