The collision between the Milky Way and Andromeda

We use an N -body/hydrodynamic simulation to forecast the future encounter between the Milky Way and the Andromeda galaxies, given present observational constraints on their relative distance, relative velocity, and masses. Allowing for a comparable amount of diffuse mass to fill the volume of the Local Group, we find that the two galaxies are likely to collide in a few billion years – within the Sun's lifetime. During the interaction, there is a chance that the Sun will be pulled away from its present orbital radius and reside in an extended tidal tail. The likelihood for this outcome increases as the merger progresses, and there is a remote possibility that our Sun will be more tightly bound to Andromeda than to the Milky Way before the final merger. Eventually, after the merger has completed, the Sun is most likely to be scattered to the outer halo and reside at much larger radii (>30 kpc). The density profiles of the stars, gas and dark matter in the merger product resemble those of elliptical galaxies. Our Local Group model therefore provides a prototype progenitor of late-forming elliptical galaxies.     

Tidal evolution of discy dwarf galaxies in the Milky Way potential: the formation of dwarf spheroidals

We conduct high-resolution collisionless N -body simulations to investigate the tidal evolution of dwarf galaxies on an eccentric orbit in the Milky Way (MW) potential. The dwarfs originally consist of a low surface brightness stellar disc embedded in a cosmologically motivated dark matter halo. During 10 Gyr of dynamical evolution and after five pericentre passages, the dwarfs suffer substantial mass loss and their stellar component undergoes a major morphological transformation from a disc to a bar and finally to a spheroid. The bar is preserved for most of the time as the angular momentum is transferred outside the galaxy. A dwarf spheroidal (dSph) galaxy is formed via gradual shortening of the bar. This work thus provides a comprehensive quantitative explanation of a potentially crucial morphological transformation mechanism for dwarf galaxies that operates in groups as well as in clusters. We compare three cases with different initial inclinations of the disc and find that the evolution is fastest when the disc is coplanar with the orbit. Despite the strong tidal perturbations and mass loss, the dwarfs remain dark matter dominated. For most of the time, the one-dimensional stellar velocity dispersion, σ, follows the maximum circular velocity, V _max, and they are both good tracers of the bound mass. Specifically, we find that   M _bound∝ V ^3.5_max  and     in agreement with earlier studies based on pure dark matter simulations. The latter relation is based on directly measuring the stellar kinematics of the simulated dwarf, and may thus be reliably used to map the observed stellar velocity dispersions of dSphs to halo circular velocities when addressing the missing satellites problem.     

Contrasting the chemical evolution of the Milky Way and Andromeda

The chemical evolution history of a galaxy hides clues about how it formed and has been changing through time. We have studied the chemical evolution history of the Milky Way (MW) and Andromeda (M31) to find which are common features in the chemical evolution of disc galaxies as well as which are galaxy-dependent. We use a semi-analytic multizone chemical evolution model. Such models have succeeded in explaining the mean trends of the observed chemical properties in these two Local Group spiral galaxies with similar mass and morphology. Our results suggest that while the evolution of the MW and M31 shares general similarities, differences in the formation history are required to explain the observations in detail. In particular, we found that the observed higher metallicity in the M31 halo can be explained by either (i) a higher halo star formation efficiency (SFE), or (ii) a larger reservoir of infalling halo gas with a longer halo formation phase. These two different pictures would lead to (i) a higher [O/Fe] at low metallicities, or (ii) younger stellar populations in the M31 halo, respectively. Both pictures result in a more massive stellar halo in M31, which suggests a possible correlation between the halo metallicity and its stellar mass.     

The abundance and radial distribution of satellite galaxies

Using detailed mock galaxy redshift surveys (MGRSs) we investigate the abundance and radial distribution of satellite galaxies. The mock surveys are constructed using large numerical simulations and the conditional luminosity function (CLF), and are compared against data from the Two Degree Field Galaxy Redshift Survey (2dFGRS). We use Monte Carlo Markov chains to explore the full posterior distribution of the CLF parameter space, and show that the average relation between light and mass is tightly constrained and in excellent agreement with our previous models and with that of Vale & Ostriker. The radial number density distribution of satellite galaxies in the 2dFGRS reveals a pronounced absence of satellites at small projected separations from their host galaxies. This is (at least partly) owing to the overlap and merging of galaxy images in the 2dFGRS parent catalogue. Owing to the resulting close-pair incompleteness we are unfortunately unable to put meaningful constraints on the radial distribution of satellite galaxies; the data are consistent with a radial number density distribution that follows that of the dark matter particles, but we cannot rule out alternatives with a constant number density core. Marginalizing over the full CLF parameter space, we show that in a ΛCDM concordance cosmology the observed abundances of host and satellite galaxies in the 2dFGRS indicate a power spectrum normalization of  σ₈≃ 0.7  . The same cosmology but with  σ₈= 0.9  is unable to match simultaneously the abundances of host and satellite galaxies. This confirms our previous conclusions based on the pairwise peculiar velocity dispersions and the group multiplicity function.     

Milky Way type galaxies in a ΛCDM cosmology

We analyse a sample of 52 000 Milky Way (MW) type galaxies drawn from the publicly available galaxy catalogue of the Millennium Simulation with the aim of studying statistically the differences and similarities of their properties in comparison to our Galaxy. Model galaxies are chosen to lie in haloes with maximum circular velocities in the range 200–250 km s⁻¹ and to have bulge-to-disc ratios similar to that of the MW. We find that model MW galaxies formed 'quietly' through the accretion of cold gas and small satellite systems. Only ≈12 per cent of our model galaxies experienced a major merger during their lifetime. Most of the stars formed ' in situ ', with only about 15 per cent of the final mass gathered through accretion. Supernovae (SNe) and active galactic nuclei (AGN) feedback play an important role in the evolution of these systems. At high redshifts, when the potential wells of the MW progenitors are shallower, winds driven by SNe explosions blow out a large fraction of the gas and metals. As the systems grow in mass, SNe feedback effects decrease and AGN feedback takes over, playing a more important role in the regulation of the star formation activity at lower redshifts. Although model MW galaxies have been selected to lie in a narrow range of maximum circular velocities, they nevertheless exhibit a significant dispersion in the final stellar masses and metallicities. Our analysis suggests that this dispersion results from the different accretion histories of the parent dark matter haloes. Statistically, we also find evidences to support the MW as a typical Sb/Sc galaxy in the same mass range, providing a suitable benchmark to constrain numerical models of galaxy formation.     

On the Andromeda to Milky Way mass ratio

We have explored the hypothesis that the total mass ratio of the two main galaxies of the Local Group, the Andromeda galaxy (M31) and the Milky Way (MW), can be constrained by measuring the tidal force induced by the surrounding mass distribution, M31 included, on the MW. We argue that the total mass ratio between the two groups can be approximated, at least qualitatively, by finding the tidal radius where the internal binding force of the MW balances the external tidal force acting on it. Since M31 is the massive tidal 'perturber' of the local environment, we have used a wide range of M31 to MW mass-ratio combinations to compute the corresponding tidal radii. Of these, only a few match the distance of the zero-tidal shell, i.e. the shell identified observationally by the outermost dwarf galaxies which do not show any sign of tidal effects. This is the key to constraining the best mass-ratio interval of the two galaxies. Our results favour a solution where the mass ratio ranges from 2 to 3, implying a massive predominance of M31.     

Masses for the Local Group and the Milky Way

We use the very large Millennium Simulation of the concordance Λ cold dark matter cosmogony to calibrate the bias and error distribution of Timing Argument estimators of the masses of the Local Group and of the Milky Way. From a large number of isolated spiral–spiral pairs similar to the Milky Way/Andromeda system, we find the interquartile range of the ratio of timing mass to true mass to be a factor of 1.8, while the 5 and 95 per cent points of the distribution of this ratio are separated by a factor of 5.7. Here, we define true mass as the sum of the 'virial' masses, M ₂₀₀, of the two dominant galaxies. For present best values of the distance and approach velocity of Andromeda, this leads to a median likelihood estimate of the true mass of the Local Group of  5.27 × 10¹² M_⊙  or  log  M _LG/M_⊙= 12.72  , with an interquartile range of [12.58, 12.83] and a 5–95 per cent range of [12.26, 13.01]. Thus, a 95 per cent lower confidence limit on the true mass of the Local Group is  1.81 × 10¹² M_⊙  . A timing estimate of the Milky Way's mass based on the large recession velocity observed for the distant satellite Leo I works equally well, although with larger systematic uncertainties. It gives an estimated virial mass for the Milky Way of  2.43 × 10¹² M_⊙  with a 95 per cent lower confidence limit of  0.80 × 10¹² M_⊙  .