Abstract

The standard Lambda Cold Dark Matter model is considered to be a triumph of theoretical astrophysics but observations of the Milky Way and its system of satellite galaxies irresistibly signal that theory is incomplete on galactic and subgalactic scales. The Stromlo Missing Satellites (SMS) Survey is a critical endeavor to investigate at what level predictions of CDM cosmology are consistent with the observed matter distribution in the Milky Way halo. It will be the deepest, most extended search for optically elusive satellite galaxies to date, covering 20 000 square degrees of sky. The international SMS Survey collaboration will exploit 150 TB of CCD images in six filters acquired by the new SkyMapper telescope of the Australian National University over the next five years, expecting on completion photometric limits 0.5–1.0 mag fainter than the Sloan Digital Sky Survey. The primary objective of the program is to characterise the baryonic and dark matter components of a complete sample of MW satellites in the Southern hemisphere to provide stringent observational constraints for improving our understanding of how the Milky Way formed and what physical processes governed galaxy formation and evolution in general.

1. Satellite Galaxies

According to cosmological theory, overdensities of cold dark matter gravitationally collapsed and formed the first structures in the Universe one billion years after the Big Bang. The gravitational pull of those dark matter clumps drew in primordial baryons in the form of hydrogen gas, providing the seeds for galaxy formation. The observational Universe today is populated with galaxies, the prime repositories of stars, the shining baryonic matter. For obvious reasons, most of the detected and catalogued galaxies are intrinsically the largest and the brightest, those that can be seen from the greatest distance and are most easily studied against the night sky. Ironically, a major limitation on our ability to develop a physically consistent model that describes how galaxies evolved out of the dark matter and baryonic ingredients comes from our incomplete picture of the vicinity of the Milky Way, in particular from the lack of a good understanding of the phenomenon “satellite galaxies." This term is broadly used for dwarf galaxy companions of the Milky Way (Figure 1), some of which contain only a few thousand stars. They cover a stellar mass range of and consist of up to 99.9% of dark matter (e.g., [17]). At the low end of the mass scale, satellites exhibit extreme low star densities so they are completely resolved, that is, transparent, and hard to detect against the screen of MW foreground stars. A few examples of such ultrafaint satellites have been discovered around the Milky Way in the recent years. Cold Dark Matter theory tells us that those are the survivors of the accretion/merger process that assembled larger galaxies. Hence, if we want to shed light on the nature of dark matter and understand how galaxies formed we have to spend a much greater amount of effort on finding and physically characterising these faintest, most elusive galaxies that exist in the Universe.

2. Cold Dark Matter Theory to the Test

While the cold dark matter paradigm successfully explains the large-scale distribution of bright baryons (i.e., luminous galaxies) in the Universe (e.g., [8]), one crucial test of the model lies in its ability to predict the properties of dark matter on subgalactic scales, most fundamentally the number and distribution of substructures around the Milky Way and their optical manifestation, the satellite galaxies. A generic result of high-resolution numerical simulations of large galaxies [912] is that the dark matter halo of the Milky Way must be filled with hundreds or even thousands of smaller dark matter condensations [13]. It is thought that the majority of these subhalos will gravitationally collect sufficient primordial hydrogen gas, turn it into stars, and form MW satellite galaxies observable today. However, the number of bound dark matter subhaloes exceeds the entire population of currently known MW satellites by at least one order of magnitude. This discrepancy is discussed in the literature as the “missing satellites" or “substructure" problem [9, 14]. There has been no shortage of solutions proposed, ranging from a drastic modification of the dark matter power spectrum [15] or self-interacting dark matter [16] to the reduction of the amount of small-scale features through warm dark matter-like models (e.g., [1720]) to careful estimates of the effects of the epoch of photoionization on the ability of subhalos to retain baryons and evolve into luminous satellite galaxies (e.g., [21]).

More recently focus has shifted to the question whether the known Milky Way satellites can actually originate from a population of dark matter subhalos. Even if the suppression of star formation by cosmic reionization and other mechanisms was quite effective in the early Universe and thus the observed satellite galaxies today are comprising only a subset of the actually present CDM substructures, their number-density distribution should be still consistent with an isotropic or mildly triaxial parent distribution. However, Kroupa and collaborators [2224] showed that the spatial distribution of the luminous Milky Way satellites that was known to be anisotropic for more than 30 years [2529] is well described by a plane of finite thickness, statistically inconsistent with the distribution of a cosmological dark matter substructure population at the 99.5% confidence level (see Figure 1). First evidence for a similar disc-of-satellites in M31 has been found by [24] that is now awaiting further testing with the results from the Pan-Andromeda Archaeological Survey (PI McConnachie; http://www.nrc-cnrc.gc.ca/eng/projects/hia/pandas.html).

Although theoretical cosmologists argue that the distribution of the few galaxies in such a “Great Disk of Satellites" almost perpendicular to the MW plane is not improbable (e.g., [30]), debating this issue remains difficult due to the incompleteness of the census of Milky Way satellites, particularly in the crucial parts of the sky away from the proposed plane, regions that are coincidently not covered by the SDSS Survey. Currently, little is known about what physics is concealed behind the “missing satellites" and “Great Disk of Satellites" phenomena but researchers agree that resolving the mystery is going to have a major impact on our understanding of dark matter and its role in the formation process and evolution of galaxies over the lifetime of the Universe.

For example, Kroupa [31] showed that dSph satellite solutions are obtained from ancient tidal dwarf galaxies, implying that the disc-of-satellites can be naturally understood as the remnant of ancient tidal arms that formed during the early assembly of the Milky Way. The almost exact match of one of the author's models with the Hercules satellite [32] is noteworthy. If a significant number of MW dSphs are confirmed to be of tidal origin this would also imply the substructure crisis to become even more severe. Alternative explanations have been put forward to reconcile this anisotropic distribution within the CDM paradigm. So Li and Helmi [33] and D'Onghia and Lake [34] proposed a scenario where the majority of satellites entered the Milky Way halo in a group rather than individually to explain the spatial and dynamical peculiarities of its satellite distribution. However, no dwarf galaxy associations exist in the vicinity of the Local Group today that are sufficiently compact to produce a disc-of-satellites [35]. On the theoretical side, Libeskind et al. [36] found that one third of their N-body simulations of MW-type halos are in agreement with the observed flattened spatial arrangement and the coherence of angular momenta of the Milky Way satellites. Accurate proper motion estimates for more satellites are needed to properly test the consistency within CDM galaxy formation model.

3. Searching for the Invisibles

Because incompleteness in the census of Milky Way satellites prevents any serious scrutiny of the currently best cosmological model, a breakthrough can only be expected once observers established the number and spatial distribution of satellite galaxies around the Milky Way including accurate estimates of those that remain hidden due to detection limits [37, 38]. Indirect evidence for a large number of undetected Milky Way, satellites was provided by Willman et al. [39]. The comparison of the spatial distributions of satellites around M31, the Milky Way and a simulated galaxy showed that the radial distribution of M31 satellites and the 15 oldest subhaloes of the simulated galaxy are in agreement, while the Milky Way satellites fall below expectations within the limits of observational completeness. Moreover, if one looks at the distribution of satellites as a function of Galactic latitude, the number of Milky Way satellites is less than the simulation predicts at low latitudes.

Northern Hemisphere
Willman et al. [40, 41] conduced a first systematic modern search for Milky Way satellites based on the Sloan Digital Sky Survey (SDSS; [42]) covering 20 percent of the entire sky. The analyses of resolved stars in both the SDSS and the Two Micron All Sky Survey (2MASS) uncovered the Willman 1 and Ursa Major galaxies. Since then, the SDSS collaboration [4348] and other groups [49, 50] have found already more new Milky Way satellites than were detected in the previous 70 years while the exploitation of the SDSS data is still in process. In an attempt to reach to even lower luminosity systems, Walsh et al. [38, Figure  5] identified a total of 30 ultrafaint candidates in the DR6 footprint (see Figure 2). Establishing the true nature of these stellar overdensities will deeply probe into unknown territory of the galaxy luminosity function.

Southern Hemisphere
This part of the sky remains completely unexplored to faint magnitude levels with modern digital imaging data. The great potential of finding new Milky Way satellites in the Southern hemisphere was demonstrated when the analysis of stars from 2MASS led to the discovery of the Canis Major dwarf galaxy (see [5153]; but see [54] for an update on the controversy of the nature of the CMa phenomenon). Covering twice the area of the SDSS footprint, the Southern hemisphere is expected to conceal at least 30 satellites more luminous than Boo II, Willman 1, or Segue 2.

4. The Stromlo Missing Satellites Survey

The Stromlo Missing Satellites (SMS) Survey (see http://msowww.anu.edu.au/~jerjen/SMS_Survey.html) was initiated to carry out the deepest, most extended search for Milky Way satellites over the next five years, covering the entire 20 000 square degrees of the Southern hemisphere. The international SMS Survey collaboration will analyse the digital images obtained by the Australian National University (ANU) SkyMapper Telescope and will use specialized data mining tools on the catalogues produced from the images by the SkyMapper team. Extensive follow-up observations of newly detected satellite candidates are absolutely essential and require the most powerful telescopes to gain an accurate picture of their physical nature. Recent follow-up works on one of the most luminous MW satellites, the CMa dwarf [54, and references therein], and the least luminous Segue 1 [55, 56] highlight some of the many challenges.

4.1. The ANU SkyMapper Telescope

Large CCD arrays are finally becoming available to construct wide-field CCD cameras at optical telescopes. Given this new research window of wide-field astronomy it is not surprising that almost every major observatory worldwide is in the process of developing survey telescopes base on that technology: Pan-STARRS of the University of Hawaii, the near-infrared 4 m Survey Telescope VISTA of the European Southern Observatory, and on an even larger financial and time scale, the 8.4 m Large Synoptic Survey Telescope (http://www.lsst.org/lsst), a wide-field telescope facility that is planned to have first light in 2015.

In this latest development, the ANU 1.35 m SkyMapper telescope (http://www.mso.anu.edu.au/skymapper/) at Siding Spring Observatory represents an investment in Australian frontier technologies of A$13 million. It is among the first of this new breed of specialised telescopes capable of scanning the night sky more quickly and sensitively than ever before featuring a CCD mosaic camera with 0.5 arcsec pixels and a 5.7-sq degree field of view. The primary goal of the facility is to perform the Southern Sky Survey (hereafter S3) over the next five years, producing 150 Terabytes of CCD images and a reference catalogue for one billion stars and galaxies. This program is a six-colour photometry survey based on the SDSS filter set ( , and ) plus a Strömgren-like filter. It will cover all 20 000 square degrees south of the celestial equator in 3845 individual pointings with photometric limits 0.5–1.0 magnitudes fainter than SDSS. The absolute stellar photometry accuracy is expected to be three percent. The S3 is carried out in six full scans of the sky each pointing of 90 seconds duration and a total integration time of 540 seconds. On completion, the S3 is expected to be 70 times more sensitive than 2MASS and 1.5 times deeper than the northern hemisphere SDSS. A single S3 epoch requires approximately six months to complete. For a detailed description of the SkyMapper observing cadence, the reduction pipeline, photometric and astrometry calibration, and accuracy we refer to Keller et al. [57].

4.2. Needles in the Haystack

The S3 catalogue is the ideal data set to continue the success story of finding new dwarf galaxies in the Milky Way halo. Out to the virial radius of 250 kpc, Milky Way satellites are completely resolved into stars. Hence, a satellite initially reveals its presence only through an unusual concentration of faint stars in the sky behind a foreground screen of disk stars, embedded in halo and possibly tidal stream stars, and superimposed on a sea of mid-to-high redshift galaxies. This excess of stars will also populate a well-defined locus in the colour-magnitude diagram depending on the satellite’s star formation history (see Figure 3).

In practice, searching for new satellites in the multidimensional parameter space of galactocentric distance, satellite SFH, star density, and angular size is a computer intensive process; the S3 catalogue will contain coordinates and photometric parameters for approximately one billion stars but only about belong to a satellite galaxy. The SMS team has already a working data mining algorithm in place that was originally developed by Willman et al. [41] and subsequently extended and improved by Walsh et al. [38]. In essence, it is a cluster detection algorithm with capabilities to analyse colour-magnitude diagrams of any given patch of sky. The important preselection of stars in the colour-magnitude diagram employs a comprehensive library of isochrone filters based on Padova evolutionary tracks [58, 59]. This selection process is repeated to scan through various distance slices. The concept of this search strategy has been extensively tested [38] by taking advantage of the publicly available Sloan Digital Sky Survey data (DR5 and DR6). In single-blind experiments every Northern Hemisphere Milky Way satellite was recovered as well as the previously unknown satellite Boötes II [50, 60].

4.3. Satellites, Light Beacons, or Extreme Globular Clusters

Stochastic fluctuations in the distribution of MW halo stars are one of the major sources of confusion in the search for satellite candidates. Our algorithm has been carefully calibrated to separate these unwanted sources from true stellar overdensities. Detection thesholds were set above the 1 in 10 000 chance for a false positive. However, we cannot exclude a priori the possibility that some of the stellar overdensities turn out to be extremely low luminosity globular clusters, such as Koposov 1 and 2 [61], or stellar condensations in large-scale stellar streams [45, 62, 63]. Any of these alternatives would be equally important to confirm.

Ultrafaint satellites such as Willman 1, Boo II, and Segue 2 with absolute magnitudes [48, 66] are less luminous than the median Milky Way globular cluster. Hence, the expected numbers of RGB and HB stars are naturally small (see Figure 3) making it necessary to obtain accurate photometry fainter than the Main Sequence Turnoff (MSTO) to confirm/negate spurious detections to the highest possible degree. As a point of reference, the MSTO of a 10 Gyr old, stellar population at 150 kpc, is at magnitude. Ideally, one aims at 1.0 mag below the MSTO ( ) in both photometric bands to have a sufficiently large number of stars available for deriving a mean age and metallicity of the population from isochrone fitting and to measure the distance modulus. Deep photometry combined with the fact that these stellar systems have a typical half-light angular size of 10–20 arcmin (see Figure 2) demands a large aperture telescope and an CCD camera with a wide field of view. Imagers like IMACS at the 6.5 m Magellan telescope (  sq arcmin FOV) or Suprime-Cam at the 8m Subaru Telescope (  sq arcmin FOV) come immediately to mind.

4.4. Kinematics and Dark Matter Contents

The importance of Milky Way satellites extends far beyond the simple construction of a complete list with spatial coordinates. Stellar velocity measurements over the last 20 years continue to show that the luminous satellites all have central velocity dispersions of the order of 10 km s-1 (e.g., Sculptor [67, 68]; Fornax [69]; Carina [70]; Sextans [71]; Ursa Minor [7274]; Draco [7375]; Leo II [76]). These results support the picture that those stellar systems are completely governed by dark matter—up to 99.9% of the total mass (e.g., [17]). Consequently, each individual satellite represents a prime laboratory to study dark matter physics in great details.

Once a candidate has been confirmed with deep imaging, red-giant candidates in the satellite will be targeted to the faintest possible limit employing multiobject spectrographs such as AAOmega at the AAT and the Michigan/MIKE Fiber Spectrograph on 6.5 m Magellan telescope. Measuring routinely radial velocities of stars ( ) to 2-3 km s-1 precision over the last couple of years enabled not only to measure the systemic velocity of a satellite but also its internal velocity dispersion (dynamical mass) assuming spherical symmetry, dynamic equilibrium, velocity isotropy, and a constant stellar mass-to-light ratio [4, 5, 55, 77]. After the SMS Survey completion we hope to have measured these quantities for as many as 30 satellites providing a stringent test for the notion that dwarf satellites have a common total mass of within the central 300 parsec [7, 70, 78, 79].

In cases where sufficient stellar spectra can be measured across the surface of the dwarf satellite, its projected velocity dispersion profile and the underlying mass density profile can also be derived [39, 50]. These results will allow a statistically robust comparison with the dark matter halo profiles determined from high-resolution N-body simulations to address the core-cusp problem, that observed galactic dark matter haloes show a density profile with a flat core (e.g., [80, 81]) while CDM haloes have divergent density (a cusp) at the centre [82, 83].

As a valuable byproduct of radial velocity measurements the high resolution spectra of the more luminous stars in these dwarf satellites will be very useful for a chemical abundance analysis (e.g., [8487]). This will enable to test the “galactic building block” hypothesis in which the Milky Way galaxy grew through the merging of smaller dwarf galaxies. Current observations of the four dwarf satellites Sculptor, Sextans, Fornax, and Carina have shown a significantly different chemical composition from stars in the Milky Way halo (e.g., [88, 89]), leading to the suggestion that not the bright, but rather the newly discovered ultrafaint MW satellites are the predicted building blocks (e.g., [90]).

5. Possible Outcomes

A cosmological model is only as good as it is capable of describing what is observed. On galaxy scales these are the internal and external physical properties of dwarf satellite galaxies. From cold dark matter simulations we expect the virialized extent of the Milky Way's halo (radius   250 kpc), which has an almost spherically symmetric shape [91, 92] to contain hundreds of subhalos following a power-law radial distribution. After the accomplishment of a full analysis of the SDSS, SkyMapper, and Pan-STARRS data, a highly complete census of MW satellites over the entire sky will be available, based on data mining algorithms that reliably pick up dwarfs out to that galactocentric distance. If the number of satellites will be still relatively small and their spatial configuration isotropic, the conclusion would be that they comprise indeed only a subset of the existing low-mass DM halos orbiting the Milky Way [9397]. The ratio of observed satellites to dark matter subhalo numbers would tell us about the effectiveness of processes such as gas physics, photoionization, and stellar feedback to suppress or evacuate baryons from low mass dark matter halos. First attempts of detailed comparisons between the Milky Way satellite luminosity function down to its faint end and subhalo predictions from CDM simulations have been presented over the last year [98102]. On the one hand, the MW dwarf satellites would be superb targets to discriminate between competing dark matter particle models including the search for annihilation signals (e.g., [11, 103]) using the upcoming new generation of ground-based Imaging Air Cherenkov Telescopes like MAGIC II and CTA. On the other hand, the reconciliation of CDM theory and observations at that level would prompt the development of innovative techniques to search for the “missing" pure dark CDM substructures in the Milky Way halo (e.g., [104]).

However, if after all the observational effort the sky distribution of any number of genuine satellites is still consistent with a planar or highly anisotropic configuration, the result would be severely at odds with the apparent sphericity of the Milky Way dark matter halo, and even more at odds with an oblate halo having the same orientation as the MW disk. Such a finding would lend strong support to the hypothesis that the majority of MW dwarf satellites have a noncosmological origin and thus lack significant amounts of dark matter. They could have formed from a gas-rich parent satellite on an eccentric near-polar orbit that interacted with the proto-MW [22] or during the dense early stages of the Universe, when the proto-Milky Way and M31 were pulled apart by the general cosmological expansion [105].

In this debate it is interesting to note that the list of the most luminous Milky Way satellites is essentially complete, modulo the possibility of a bright satellite lurking behind the Galactic disk. Accordingly, these “classical" satellites will always show a disk-like arrangement independently of how many low luminous systems will be added to the population in the future. Hence there is a fair chance that the MW dwarf satellites may themselves be divided into two subclasses with different origins. While we look forward in great anticipation to the Stromlo Missing Satellites Survey and similar programs there is no doubt that MW dwarf satellites will play a central role in near-field cosmology for decades to come.

Acknowledgments

The author thanks both referees for their comments that helped to improve the presentation and the quality of this paper. He also acknowledges the financial support from the Go8-DAAD-Australia/Germany Joint Research Co-operation Scheme and from the Access to Major Research Facilities Programme which is a component of the International Science Linkages established under the Australian Government Innovation Statement, Backing Australia's Ability.