Satellite systems around galaxies in hydrodynamic simulations

We investigate the properties of satellite galaxies formed in N -body/SPH simulations of galaxy formation in the ΛCDM cosmology. The simulations include the main physical effects thought to be important in galaxy formation and, in several cases, produce realistic spiral discs. In total, a sample of nine galaxies of luminosity comparable to the Milky Way was obtained. At magnitudes brighter than the resolution limit,   M_V =−12  , the luminosity function of the satellite galaxies in the simulations is in excellent agreement with data for the Local Group. The radial number density profile of the model satellites, as well as their gas fractions also match observations very well. In agreement with previous N -body studies, we find that the satellites tend to be distributed in highly flattened configurations whose major axis is aligned with the major axis of the (generally triaxial) dark halo. In two out of three systems with sufficiently large satellite populations, the satellite system is nearly perpendicular to the plane of the galactic disc, a configuration analogous to that observed in the Milk Way. The discs themselves are perpendicular to the minor axis of their host haloes in the inner parts, and the correlation between the orientation of the galaxy and the shape of the halo persists even out to the virial radius. However, in one case the disc's minor axis ends up, at the virial radius, perpendicular to the minor axis of the halo. The angular momenta of the galaxies and their host halo tend to be well aligned.     

Dwarf galaxies: Important clues to galaxy formation

The smallest dwarf galaxies are the most straight forward objects in which to study star formation processes on a galactic scale. They are typically single cell star forming entities, and as small potentials in orbit around a much larger one they are unlikely to accrete much (if any) extraneous matter during their lifetime (either intergalactic gas, or galaxies) because they will typically lose the competition with the much larger galaxy. We can utilise observations of stars of a range of ages to measure star formation and enrichment histories back to the earliest epochs. The most ancient objects we have ever observed in the Universe are stars found in and around our Galaxy. Their proximity allows us to extract from their properties detailed information about the time in the early Universe into which they were born. A currently fashionable conjecture is that the earliest star formation in the Universe occurred in the smallest dwarf galaxy sized objects. Here I will review some recent observational highlights in the study of dwarf galaxies in the Local Group and the implications for understanding galaxy formation and evolution. This revised version was published online in August 2006 with corrections to the Cover Date.     

The mass of the Andromeda galaxy

This paper argues that the Milky Way galaxy is probably the largest member of the Local Group. The evidence comes from estimates of the total mass of the Andromeda galaxy (M31) derived from the three-dimensional positions and radial velocities of its satellite galaxies, as well as the projected positions and radial velocities of its distant globular clusters and planetary nebulae. The available data set comprises 10 satellite galaxies, 17 distant globular clusters and nine halo planetary nebulae with radial velocities. We find that the halo of Andromeda has a mass of together with a scalelength of 90 kpc and a predominantly isotropic velocity distribution. For comparison, our earlier estimate for the Milky Way halo is Although the error bars are admittedly large, this suggests that the total mass of M31 is probably less than that of the Milky Way . We verify the robustness of our results to changes in the modelling assumptions and to errors caused by the small size and incompleteness of the data set.
Our surprising claim can be checked in several ways in the near future. The numbers of satellite galaxies, planetary nebulae and globular clusters with radial velocities can be increased by ground-based spectroscopy, while the proper motions of the companion galaxies and the unresolved cores of the globular clusters can be measured using the astrometric satellites Space Interferometry Mission ( SIM ) and Global Astrometric Interferometer for Astrophysics ( GAIA ). Using 100 globular clusters at projected radii 20 R 50 kpc with both radial velocities and proper motions, it will be possible to estimate the mass within 50 kpc to an accuracy of 20 per cent. Measuring the proper motions of the companion galaxies with SIM and GAIA will reduce the uncertainty in the total mass caused by the small size of the data set to 22 per cent.     

The Morphological Evolution of Galaxy Satellites

We study the evolution of galaxy satellites with high resolutionN-body simulations. Satellites are modelled as replicas of typical low and high surface brightness galaxies (LSBs and HSBs). Encounters on high eccentricity orbits (as typical in hierarchical models of galaxy formation) strip LSBs of most of their stars and tend to decrease their surface brightness. In contrast, bar instability in HSBs leads to substantial loss of angular momentum of the stellar component and to an increase in central surface brightness. In both cases the remnant resembles a spheroidal galaxy with an exponential surface brightness profile. A simple modelling of colour evolution and interaction-driven star formation gives M/L ratios for the remnants that are roughly consistent with observations. These results suggest an evolutionary scenario for the dwarf galaxies in our Local Group, faint dSphs being the descendants of LSBs and brighter dSphs/dEs being the final state of HSB satellites. This revised version was published online in July 2006 with corrections to the Cover Date.     

The growth and assembly of a massive galaxy at z∼ 2

We study the stellar mass assembly of the Spiderweb galaxy  (MRC 1138−262)  , a massive   z = 2.2  radio galaxy in a protocluster and the probable progenitor of a brightest cluster galaxy. Nearby protocluster galaxies are identified and their properties are determined by fitting stellar population models to their rest-frame ultraviolet to optical spectral energy distributions. We find that within 150 kpc of the radio galaxy the stellar mass is centrally concentrated in the radio galaxy, yet most of the dust-uncorrected, instantaneous star formation occurs in the surrounding low-mass satellite galaxies. We predict that most of the galaxies within 150 kpc of the radio galaxy will merge with the central radio galaxy by   z = 0  , increasing its stellar mass by up to a factor of ≃2. However, it will take several hundred Myr for the first mergers to occur, by which time the large star formation rates are likely to have exhausted the gas reservoirs in the satellite galaxies. The tidal radii of the satellite galaxies are small, suggesting that stars and gas are being stripped and deposited at distances of tens of kpc from the central radio galaxy. These stripped stars may become intracluster stars or form an extended stellar halo around the radio galaxy, such as those observed around cD galaxies in cluster cores.     

Discs of satellites: the new dwarf spheroidals

The spatial distributions of the most recently discovered ultra-faint dwarf satellites around the Milky Way and the Andromeda galaxy are compared to the previously reported discs-of-satellites (DoS) of their host galaxies. In our investigation, we pay special attention to the selection bias introduced due to the limited sky coverage of Sloan Digital Sky Survey (SDSS). We find that the new Milky Way satellite galaxies follow closely the DoS defined by the more luminous dwarfs, thereby further emphasizing the statistical significance of this feature in the Galactic halo. We also note a deficit of satellite galaxies with Galactocentric distances larger than  100 kpc  that are away from the DoS of the Milky Way. In the case of Andromeda, we obtain similar results, naturally complementing our previous finding and strengthening the notion that the DoS are optical manifestations of a phase-space correlation of satellite galaxies.     

The growth of central and satellite galaxies in cosmological smoothed particle hydrodynamics simulations

We examine the accretion and merger histories of central and satellite galaxies in a smoothed particle hydrodynamics (SPH) cosmological simulation that resolves galaxies down to  7 × 10⁹ M_⊙  . Most friends-of-friends haloes in the simulation have a distinct central galaxy, typically 2–5 times more massive than the most massive satellite. As expected, satellites have systematically higher assembly redshifts than central galaxies of the same baryonic mass, and satellites in more massive haloes form earlier. However, contrary to the simplest expectations, satellite galaxies continue to accrete gas and convert it to stars; the gas accretion declines steadily over a period of 0.5–1 Gyr after the satellite halo merges with a larger parent halo. Satellites in a cluster mass halo eventually begin to lose baryonic mass. Typically, satellites in our simulation are 0.1–0.2 mag bluer than in models that assume no gas accretion on to satellites after a halo merger. Since   z = 1  , 27 per cent of central galaxies (above  3 × 10¹⁰ M_⊙  ) and 22 per cent of present-day satellite galaxies have merged with a smaller system above a 1:4 mass ratio; about half of the satellite mergers occurred after the galaxy became a satellite and half before. In effect, satellite galaxies can remain 'central' objects of halo substructures, with continuing accretion and mergers, making the transition in assembly histories and physical properties a gradual one. Implementing such a gradual transformation in semi-analytic models would improve their agreement with observed colour distributions of satellite galaxies in groups and with the observed colour dependence of galaxy clustering.     

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.     

Eclipsing Binaries in Local Group Galaxies

A review is presented of the progress that has been made in the last 3 years towards quantifying the properties of high-mass detached and semi-detached eclipsing binaries in Local Group galaxies. Comparisons between these observational results on masses, radii, temperatures and luminosities for stars in detached binaries and evolution models for single stars at the appropriate metallicity are found to be very good. New evolution models for interacting binaries passing through case A mass exchange are being calculated, and indicate a requirement for some mass loss to find agreement with the observational data. The observational data on such semi-detached systems show similar properties to those in the Milky Way galaxy. The directly-determined distances to all these eclipsing binaries are proving to be most valuable for strengthening the distance scale amongst the Local Group 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.     

Star clusters around recoiled black holes in the Milky Way halo

Gravitational wave emission by coalescing black holes (BHs) kicks the remnant BH with a typical velocity of hundreds of  km s⁻¹  . This velocity is sufficiently large to remove the remnant BH from a low-mass galaxy but is below the escape velocity from the Milky Way (MW) galaxy. If central BHs were common in the galactic building blocks that merged to make the MW, then numerous BHs that were kicked out of low-mass galaxies should be freely floating in the MW halo today. We use a large statistical sample of possible merger tree histories for the MW to estimate the expected number of recoiled BH remnants present in the MW halo today. We find that hundreds of BHs should remain bound to the MW halo after leaving their parent low-mass galaxies. Each BH carries a compact cluster of old stars that populated the core of its original host galaxy. Using the time-dependent Fokker–Planck equation, we find that the present-day clusters are  ≲1 pc  in size, and their central bright regions should be unresolved in most existing sky surveys. These compact systems are distinguishable from globular clusters by their internal (Keplerian) velocity dispersion greater than 100 km s⁻¹ and their high mass-to-light ratio owing to the central BH. An observational discovery of this relic population of star clusters in the MW halo would constrain the formation history of the MW and the dynamics of BH mergers in the early Universe. A similar population should exist around other galaxies and may potentially be detectable in M31 and M33.     

Quantifying dust and the ultraviolet radiation density in the local Universe

A sample of local galaxies for which far-infrared and ultraviolet fluxes are available is used to estimate the characteristic dust extinction in galaxies and to test whether standard dust properties are plausible. Assuming galaxies can be characterized by a single dust optical depth (certainly not valid for galaxies with a dominant starburst component), the infrared excess and ultraviolet colours of local galaxies are found to be consistent with normal Milky Way dust, with a mean value for   E ( B − V )  of 0.16. A significant fraction of the dust heating is caused by older, lower-mass stars, and this fraction increases towards earlier galaxy types.
Analysis of  ( F _FIR/ F _UV)  versus ultraviolet colour diagrams for starburst galaxies in terms of a simple screen dust model does not support a Calzetti (1997) rather than a Milky Way extinction law, though the absence of the expected strong 2200- Å feature in several galaxies with IUE spectra does show that more detailed radiative transfer models are needed – probably with nonspherical geometry.
A simple treatment in which the  100/60-μm  flux ratio is used to subtract the optically thick starburst contribution to the far-infrared radiation results in lower extinction estimates for the optically thin cirrus component, with a mean   E ( B − V )  of 0.10.
The ultraviolet luminosity density, corrected for dust extinction, is derived and a value for the local mean star formation rate inferred. This is consistent with previous estimates from ultraviolet surveys and from  Hα  surveys.     

Environmental effects on satellite galaxies: the link between concentration, size and colour profile

Using the Sloan Digital Sky Survey Data Release 4 group catalogue of Yang et al., we investigate sizes, concentrations, colour gradients and surface brightness profiles of central and satellite galaxies. We compare central and satellite galaxies at fixed stellar mass, in order to disentangle environmental from stellar mass dependencies. Early- and late-type galaxies are defined according to concentration. We find that at fixed stellar mass, late-type satellite galaxies have smaller radii and larger concentrations than late-type central galaxies. No such differences are found for early-type galaxies. We have also constructed surface brightness and colour profiles for the central and satellite galaxies in our sample. We find that late-type satellite galaxies have a lower surface brightness and redder colours than late-type central galaxies. We show that all observed differences between satellite and central galaxies can be explained by a simple fading model, in which the star formation in the disc decreases over time-scales of 2–3 Gyr after a galaxy becomes a satellite. Processes that induce strong morphological changes (e.g. harassment) and processes that strip the galaxy of its entire interstellar medium need not to be invoked in order to explain the environmental dependencies we find.     

The sizes of minivoids in the local Universe: an argument in favour of a warm dark matter model?

Using high-resolution simulations within the cold dark matter (CDM) and warm dark matter (WDM) models, we study the evolution of small-scale structure in the local volume, a sphere of 8-Mpc radius around the Local Group. We compare the observed spectrum of minivoids in the local volume with the spectrum of minivoids determined from the simulations. We show that the ΛWDM model can easily explain both the observed spectrum of minivoids and the presence of low-mass galaxies observed in the local volume, provided that all haloes with circular velocities greater than 20 km s⁻¹ host galaxies. On the contrary, within the ΛCDM model the distribution of the simulated minivoids reflects the observed one if haloes with maximal circular velocities larger than  35 km s⁻¹  host galaxies. This assumption is in contradiction with observations of galaxies with circular velocities as low as 20 km s⁻¹ in our local Universe. A potential problem of the ΛWDM model could be the late formation of the haloes in which the gas can be efficiently photoevaporated. Thus, star formation is suppressed and low-mass haloes might not host any galaxy at all.     

On the motion of the Local Group and its substructures

We consider the problem of the relative motion both of the substructures of the Local Group of galaxies (revealed via the S-tree method), and of the velocity of the Local Group itself. The existence of statistically significant bulk flow of the Milky Way subsystem is shown via a 3D reconstruction procedure, which uses information on the radial velocities of the galaxies but does not take account of their distances. Once the bulk motion of the substructures is estimated we also consider, in combination with the observed cosmic microwave background (CMB) dipole, the mean velocity of the Local Group itself. Assigning to the Local Group the mean motion of its main substructures, we evaluate its peculiar velocity in Milky Way frame V _LG→MW = (−7 ± 303, −15 ± 155, +177 ± 144) or 178 km s⁻¹ toward galactic coordinates l  = 245 and b  = +85. Combined with the CMB dipole V _MW→CMB, we obtain a Local Group velocity in CMB frame: V _LG→CMB = (−41 ± 303, −497 ± 155, 445 ± 144) or 668 km s⁻¹ towards l  = 265 and b  = 42. This estimation is in good agreement, within the 1 σ level, with the estimation of Yahil et al.     

Inhomogeneous chemical evolution of dwarf spheroidal galaxies

Stellar abundance pattern of n-capture elements such as barium is used as a powerful tool to infer how the star formation proceeded in dwarf spheroidal (dSph) galaxies. It is found that the abundance correlation of barium with iron in stars belonging to dSph galaxies orbiting the Milky Way, i.e., Draco, Sextans, and Ursa Minor have a feature similar to that in Galactic metal-poor stars. The common feature of these two correlations can be realized by our in homogeneous chemical evolution model based on the supernova-driven star formation scenario if dSph stars formed from gas with a velocity dispersion of ∼ 26 km s^-1. This velocity dispersion together with the stellar luminosities strongly suggest that dark matter dominated dSph galaxies. The tidal force of the Milky Way links this velocity dispersion with the currently observed value ≲ 10 km s^-1 by stripping the dark matter in dSph galaxies. As a result, the total mass of each dSph galaxy is found to have been originally ∼ 25 times larger than at present. In this model, supernovae immediately after the end of the star formation can expel the remaining gas over the gravitational potential of the dSph galaxy. This revised version was published online in August 2006 with corrections to the Cover Date.     

The effect of dwarf galaxy disruption in semi-analytic models

We present results for a galaxy formation model that includes a simple treatment for the disruption of dwarf galaxies by gravitational forces and galaxy encounters within galaxy clusters. This is implemented a posteriori in a semi-analytic model by considering the stability of cluster dark matter subhaloes at   z = 0  . We assume that a galaxy whose dark matter substructure has been disrupted will itself disperse, while its stars become part of the population of intracluster stars responsible for the observed intracluster light. Despite the simplicity of this assumption, our results show a substantial improvement over previous models and indicate that the inclusion of galaxy disruption is indeed a necessary ingredient of galaxy formation models. We find that galaxy disruption suppresses the number density of dwarf galaxies by about a factor of 2. This makes the slope of the faint end of the galaxy luminosity function shallower, in agreement with observations. In particular, the abundance of faint, red galaxies is strongly suppressed. As a result, the luminosity function of red galaxies and the distinction between the red and the blue galaxy populations in colour–magnitude relationships are correctly predicted. Finally, we estimate a fraction of intracluster light comparable to that found in clusters of galaxies.     

Disk galaxy evolution: from the Milky Way to high-redshift disks

We develop a detailed model of the Milky Way (a `prototypical' disk galaxy) and extend it to other disks with the help of some simple scaling relations, obtained in the framework of Cold Dark Matter models. This phenomenological (`hybrid') approach to the study of disk galaxy evolution allows us to reproduce successfully a large number of observed properties of disk galaxies in the local Universe and up to redshift z ∼ 1. The important conclusion is that, on average, massive disks have formed the bulk of their stars earlier than their lower mass counterparts: the `star formation hierarchy' has been apparently opposite to the `dark matter assembly' hierarchy. It is not yet clear whether `feedback' (as used in semi-analytical models of galaxy evolution) can explain that discrepancy. This revised version was published online in August 2006 with corrections to the Cover Date.     

The evolution of barium and europium in local dwarf spheroidal galaxies

By means of a detailed chemical evolution model, we follow the evolution of barium (Ba) and europium (Eu) in four Local Group Dwarf Spheroidal (dSph) galaxies, in order to set constraints on the nucleosynthesis of these elements and on the evolution of this type of galaxies compared with the Milky Way. The model, which is able to reproduce several observed abundance ratios and the present-day total mass and gas mass content of these galaxies, adopts up-to-date nucleosynthesis and takes into account the role played by supernovae (SNe) of different types (II, Ia) allowing us to follow in detail the evolution of several chemical elements (H, D, He, C, N, O, Mg, Si, S, Ca, Fe, Ba and Eu). By assuming that Ba is a neutron-capture element produced in low-mass asymptotic giant branch stars by s-process but also in massive stars (in the mass range 10–30 M_⊙) by r-process, during the explosive event of SNe of Type II, and that Eu is a pure r-process element synthesized in massive stars also in the range of masses 10–30 M_⊙, we are able to reproduce the observed [Ba/Fe] and [Eu/Fe] as functions of [Fe/H] in all four galaxies studied. We confirm also the important role played by the very low star formation (SF) efficiencies (ν= 0.005–0.5 Gyr⁻¹) and by the intense galactic winds (6–13 times the star formation rate) in the evolution of these galaxies. These low SF efficiencies (compared to the one for the Milky Way disc) adopted for the dSph galaxies are the main reason for the differences between the trends of [Ba/Fe] and [Eu/Fe] predicted and observed in these galaxies and in the metal-poor stars of our Galaxy. Finally, we provide predictions for Sagittarius galaxy for which data of only two stars are available.