---
abstract: 'We present a new semi-analytical model of the population of satellite galaxies of the Milky Way, aimed at estimating the effect of the spatial structure of reionisation at galaxy scale on the properties of the satellites. In this model reionisation can be either: (A) externally-driven and uniform, or (B) internally-driven, by the most massive progenitor of the Milky-Way. In the latter scenario the propagation of the ionisation front and photon dilution introduce a delay in the photo-evaporation of the outer satellites’ gas with respect to the inner satellites. As a consequence, outer satellites experience a longer period of star formation than those in the inner halo. We use simple models to account for star formation, the propagation of the ionisation front, photo-evaporation and observational biases. Both scenarios yield a model satellite population at z=0 that matches the observed luminosity function and mass-to-light ratios. However, the predicted population for scenario (B) is significantly more extended spatially than for scenario (A), by about 0.3 dex in distance, resulting in a much better match to the observations. The survival of this structural signature imprinted by the local UV field during reionisation on the radial distribution of satellites makes it a promising tool for studying the reionisation epoch at galaxy scale in the Milky Way and nearby galaxies resolved in stars with forthcoming large surveys. However, more work is needed to determine how the effect reported here can be disentangled from that of cosmic variance between different realisations of Milky Way haloes.'
author:
- |
    Pierre Ocvirk$^{{1}}$ &Dominique Aubert$^{{1}}$\
    $^1$Observatoire Astronomique de Strasbourg, 11 rue de l’Université, 67000 Strasbourg, France.\
bibliography:
- 'mybib.bib'
date: 'Typeset ; Received 15 April 2011 / Accepted 1 August 2011'
title: 'A signature of the internal reionisation of the Milky Way?'
---

\[firstpage\]

Galaxy: formation - galaxies: dwarf, luminosity function - Local Group - cosmology: reionization - diffuse radiation

Introduction
============

Reionisation took place as the first stars of the universe brought the dark ages to an end. It seems now well established that reionisation also had a profound impact on galaxy formation. The metagalactic UV radiation field may be responsible for evaporating the gas of low-mass galaxies [@gnedin2000; @hoeft2006], affecting their star formation and therefore also the buildup of the galactic halo [@bekkichiba2005]. This seems to provide an acceptable solution for the “missing satellites problem” [@klypin1999; @moore1999], by inhibiting star formation in low mass galaxies at early times [@bullock2000].

Several recent semi-analytical models (SAMs) have shown that the observed luminosity function (hereinafter LF) of the satellites of the Milky Way (hereinafter MW) can be well reproduced under a number of simple assumptions for the baryonic physics [@koposov2009; @munoz2009; @busha2010; @maccio2010; @font2011]. They also reproduce the apparent common mass scale of $10^7 {\ensuremath{\mathrm{M}_{\odot}}}$ [@strigari2008], although the latter is hotly debated [@aden2009; @walker2010]. They suggest that the ultra-faint dwarf galaxies (hereinafter UFDs) discovered by the SDSS [@martin2004; @willman2005short; @zucker2006; @belokurov2007short; @irwin2007short; @walsh2007] are reionisation fossils, living in sub-haloes of about $10^{6-9} {\ensuremath{\mathrm{M}_{\odot}}}$. The success of these models invites us to reconsider UFDs, not as a problem, but as probes of the reionisation epoch. Most past studies implemented a uniform, instantaneous reionisation by simply suppressing star formation in the sub-haloes below ${{\rm z_{reion}}}$ [@koposov2009; @busha2010; @maccio2010]. Others tracked the evolution of the physical properties of the gas (temperature, neutral fraction) subject to the UV background from external and local sources, and SuperNova (SN) feedback, in order to determine self-consistently the satellite’s reionisation redshift [@munoz2009; @font2011]. [@lunnan2011] studied the effect of a complex, patchy reionisation on the satellites properties. In this letter, we investigate the effect of varying the structure of the UV field at reionisation on the properties of the satellites of the MW. We choose 2 deliberately simple configurations: a uniform background, and a unique, galacto-centric source. We show that the latter produces a significantly more extended satellite population.

We first present our model and its ingredients (Sec. \[s:model\]), and then our results in terms of LF, masses, and radial distribution of satellites (Sec. \[s:results\]). Our conclusions are presented in Sec. \[s:conclusion\].

A semi-analytical model of the satellites of the Milky Way {#s:model}
==========================================================

Our SAM is partly inspired by [@busha2010] and [@koposov2009] (hereinafter B10 and K09), but describes the reionisation epoch in more detail. We allow for a non-steady, non-uniform UV background and compute individual evaporation times for the gas hosted by the dark matter sub-haloes.

Ingredients {#s:ingredients}
-----------

The life of a model dark matter sub-halo in the range $10^{7-9}{\ensuremath{\mathrm{M}_{\odot}}}$ (typical of a MW satellite at z=0) can be summarized as follows: after its formation, it grows until it reaches a mass sufficient for forming stars. During reionisation, ionisation fronts (hereinafter I-front) of galactic or external origin propagate, and eventually overtake the sub-halo at a redshift ${{\rm z_{I}}}$. The latter starts to lose its gas via photo-evaporation, which rapidly suppresses star formation, between 0.4-2 Gyr after the Big Bang. Then it undergoes a purely passive evolution for the next $\approx 6$ Gyr, until it is accreted by the MW, around $z\approx 1$. Depending on its impact parameters and concentration it will be tidally and ram-pressure stripped. Finally, it may get disrupted, losing all of its stars to the halo, or conserve a fraction of its pre-accretion mass and stellar content, and survive as a satellite until today. The purpose of the SAM described here is to model this chain of physical processes, with a number of deliberate simplifications, in order to produce a population of satellites and its observables. The major ingredients are listed below, in the order of the modelling steps.

1.  [Sub-halo tracks: as in B10 we use the Via Lactea II (hereinafter VLII) data [@diemand2007]. The simulation represents a MW-size halo (about $2.10^{12} {\ensuremath{\mathrm{M}_{\odot}}}$) that has not had any recent major merger, making it a suitable host for a MW-like disk galaxy. The dark matter particle mass is $4.10^3 {\ensuremath{\mathrm{M}_{\odot}}}$, which allows about 20000 sub-haloes to be resolved, down to a few times $10^{5}{\ensuremath{\mathrm{M}_{\odot}}}$.]{}

2.  [Threshold mass for star formation ${\rm M_{t}}$: a sub-halo is allowed to form stars when it reaches ${{\rm T_{vir}}}=8.10^3$K, which translates to ${\rm M_{t}}\approx 10^{7-8} {\ensuremath{\mathrm{M}_{\odot}}}$. This reflects the typical halo mass where atomic hydrogen cooling becomes efficient, as in B10.]{}

3.  [Following B10, once the threshold mass is reached, the halo forms stars at a rate ${\rm SFR} = \epsilon ({{\rm{f_{cold}}}{{\rm M}}_7})^{\alpha}$, with $\epsilon = 10^{-5}$, ${\rm{f_{cold}}}=\Omega_{\rm b}/\Omega_{\rm m}$, $\alpha =2$, and ${{\rm M}}_7={{\rm M}}/10^7 {\ensuremath{\mathrm{M}_{\odot}}}$, where ${{\rm M}}$ is the dark matter mass of the sub-halo. ]{}

4.  [Reionisation redshifts: we define the reionisation redshift ${{\rm z_{I}}}$ as the time where the I-front [*reaches*]{} the sub-halo. The details for both external and internal scenarios are summarized in Sec. \[s:scenarios\]. At variance with most of the literature, star formation does [*not*]{} end at ${{\rm z_{I}}}$, as photo-evaporation is not considered instantaneous here.]{}

5.  [Photo-evaporation: we compute evaporation times for the gas in the haloes using Eq. (5) of [@iliev2005], as a function of the normalized flux ${\rm F_0}$, the dark matter mass ${{\rm M}}$ of the halo, and its collapse redshift ${{\rm z_{coll}}}$: $${\rm t_{ev}}= A  {{\rm M}}_{7}^{\rm B} {\rm F_0}^{\rm C+D\,\log_{10}({\rm F_0})} \left[ E + F \left( \frac{1+{{\rm z_{coll}}}}{10} \right) \right] {\rm Myr}  \, ,
    \label{eq:tevap}$$ where (A,B,C,D,E,F) are taken from [@iliev2005] for the case of a $5.10^{4}$ K black body. Expectedly, ${\rm t_{ev}}$ increases with increasing mass and decreasing UV flux. For simplicity we assume a constant ${{\rm z_{coll}}}=15$ for all haloes. Photo-evaporation starts at ${{\rm z_{I}}}$, and ${\rm t_{ev}}$ is converted to a $\Delta {\rm z_{ev}}$, so that evaporation finishes at a redshift ${\rm z_{ev}}={{\rm z_{I}}}+ \Delta {\rm z_{ev}}$. The SFR of the sub-halo is set to 0 for $z<{\rm z_{ev}}$. The values taken by ${\rm F_0}$ in our 2 reionisation scenarios are detailed in Sec. \[s:scenarios\].]{}

6.  The accretion redshift ${\rm z_{acc}}$ is computed for each sub-halo using the VLII data. We consider that accretion on the MW halo results in instantaneous gas removal (if any left) from the satellite through ram-pressure stripping. We set ${\rm SFR}=0$ for $z<{\rm z_{acc}}$. We can thus write the SFR$({\ensuremath{\mathrm{M}_{\odot}}}/{\rm yr})$ prescription more synthetically as:

    $$\displaystyle
      {\rm SFR}= 
      \quad
      \left\{\begin{array}{l}
       \epsilon ({\rm{f_{cold}}}{{\rm M}}_7 )^{\alpha} \mbox{ if}  
       \left\{\begin{array}{l}
       {{\rm T_{vir}}}\geq  8.10^3 {\rm K} \mbox{ and}\\
       z  >  {\rm max}({\rm z_{ev}},{\rm z_{acc}}) 
       \end{array} \right. \\
       0 \mbox{ otherwise} \, .
      \end{array}\right.
    \label{eq:SFR}$$

    This formalism gives rise to typical periods of star formation of several 100 Myr between $z=15-6$, with SFR=$10^{-6}-10^{-2} {\ensuremath{\mathrm{M}_{\odot}}}/{\rm yr}$ for the observable satellites.

7.  [Tidal stripping and disruption: A halo is considered disrupted if it loses more than 99% of its pre-accretion mass. However, for a surviving halo we assume that the stellar content remains untouched. This is similar to the approach of K09 and [@maccio2010].]{}

8.  [Stellar population models: integrating eq. (\[eq:SFR\]) over time yields the stellar mass and the individual star formation history (SFH) of the sub-haloes. We use the latter along with the [@BC03] models to assign luminosities to our satellite population, assuming a constant low metallicity Z=0.0004.]{}

Reionisation scenarios {#s:scenarios}
----------------------

We consider two scenarios, where reionisation of the IGM and the subsequent evaporation of the satellites’ gas is driven either by an external or a galacto-centric source. Note that here the 2 scenarios are mutually exclusive.

### [**Scenario (A): external reionisation**]{} {#s:external .unnumbered}

In this scenario, the UV background responsible for photo-evaporating the gas of the satellites is of extragalactic origin. We assume it is uniform, as in a majority of previous studies, and appears at ${{\rm z_{I}}}=15$ for all sub-haloes, corresponding to the beginning of the reionisation of the IGM [@kogut2003]. We tune the value of the background flux to a constant ${\rm F_0}= 0.03$, in the units of [@iliev2005], so as to match the observed LF of the satellites. This translates into ${\rm J_{21}}=0.03$, which is compatible with the fluxes predicted by [@aubert2010] for the considered redshift range, although in the upper envelope.

### [**Scenario (B): internal reionisation by a galacto-centric source**]{} {#scenario-b-internal-reionisation-by-a-galacto-centric-source .unnumbered}

While reionisation appears patchy in large scale computations [@zahn2007], the morphology of I-fronts at the scale of the MW progenitor (5 comoving Mpc) is rather unclear, and depends strongly on the Lyman continuum escape fraction ${\rm f_{esc}}$ [@iliev2011], which is currently very uncertain. Here we assume that star formation within the most massive progenitor of the MW provides the photons responsible for evaporating the satellites’ gas and the intervening IGM. This model results in a radially decreasing reionisation redshift profile in the IGM around the emitter, similar to those of Fig. 1 of [@lunnan2011]. However this model can not account for any stochasticity of the reionisation history of the IGM.

1.  [ we model the I-front ${{\rm R_{I}}}(z)$ as a cosmological Strömgren sphere [@barkana2001] with a clumping factor C=1. The SFR of the galacto-centric source is an extrapolation from Fig. 4 of [@salvadori2010]: ${\rm SFR(z)=SFR_{25} 10^{(25-z)/10}}$, with ${\rm SFR_{25}=SFR(z=25)=0.02}\, {\ensuremath{\mathrm{M}_{\odot}}}$/yr. We assume that a rate SFR=1[$\mathrm{M}_{\odot}$]{}/yr produces $1.6 \, 10^{53}$ ionising photons/s as in [@shapiro2004]. We set the start of star formation at $z=25$, as suggested by theoretical studies [@bromm2009]. Numerical simulations from [@razoumov2010] show that at z$>10$, ${\rm f_{esc}}$ is close to 1. Therefore we use a constant ${\rm f_{esc}}=1$ throughout the model.]{}

2.  [ for a satellite at a distance ${{\rm R_{GC}}}$ from the Galactic center, photo-evaporation starts at ${{\rm z_{I}}}$, when the expanding I-front reaches the satellite, i.e. ${{\rm R_{I}}}({{\rm z_{I}}})={{\rm R_{GC}}}({{\rm z_{I}}})$. We convert our SFR to ${\rm F_0}$ using ${\rm F_0}=16\,  \,SFR(z) \, ({10 \, {\rm kpc}}/{{{\rm R_{GC}}}})^2$, following the conversions of [@shapiro2004].]{}

Observational sample and biases
-------------------------------

We consider the same observed satellite population as in Tab. 2 of K09, consisting of the 11 classical and the 10 recently discovered SDSS satellites. The galacto-centric distances are taken from [@kroupa2010]. The observed luminosity function is corrected for the sky coverage of DR5 and the detection efficiency given by [@koposov2008]. Detected model satellites are required to be within 260 kpc of the galactic center (the maximum distance where the tip of the red giant branch is detected by SDSS). Moreover they need to be within the completeness radius at a given magnitude as given by eq.(8) of K09 at z=0. Finally, their detection efficiency is computed using eq.(8) of [@koposov2008].

Results {#s:results}
=======

The LF is degenerated with respect to the parameters $(\epsilon,{\rm F_0})$: a decrease in star formation efficiency $\epsilon$ can be compensated by a longer star formation period, allowed by a smaller UV background flux. This is similar to the $\epsilon-{\rm z_{re}}$ degeneracy reported in B10. In the present paper, ${\rm F_0}$ was set so as to achieve evaporation of the satellites at high redshift [ on average]{}, i.e. $<{\rm z_{ev}}>\approx 6$.

Luminosity function
-------------------

Model and observed cumulative LFs are compared in Fig. \[f:LF\]. Both models with reionisation agree equally well the observations. At he faint end, they differ slightly, but are both equally likely given the available data. At the bright end we underestimate the luminosity of the satellites, because of missing physics for the high-mass objects. For instance, MW’s most massive satellites, such as the large and small Magellanic clouds, are still forming stars at z=0, which is forbidden in our model due to the shutdown of star formation subsequent to accretion on the MW halo and ram-pressure stripping. Our internal and external reionisation scenarios are indistinguishable when considering the existing data on the LF alone. To illustrate the overall effect of reionisation, we also show the results of a model identical to model (A) but without reionisation. It overestimates the abundance of satellites by a factor 10 for all but the 2 brightest objects and is therefore not favored in our framework. We refer the reader to [@li2010] for an in-depth exploration of models without reionisation.

[![Observed and model cumulative luminosity functions. The shaded areas show the observational error bars corresponding to Poisson noise, corrected for biases for the SDSS satellites. Both models with reionisation provide a good fit to the observations.[]{data-label="f:LF"}](LFs_zfirststar=15_with_noreion_highepsilon.ps "fig:"){width="0.95\linewidth"}]{}

Masses
------

[![Masses of the faint satellites. [*Top:*]{} central masses (${\rm M_{300}}$), models and as derived by Strigari et al. (2008). [*Bottom:*]{} total masses of simulated satellites. []{data-label="f:m300"}](Masses_300_total_zfirststar=15.ps "fig:"){width="1.11\linewidth"}]{}

The central masses ($\rm M_{300}$, within the 300 central pc) of the model satellites are compared with observations in the top panel of Fig. \[f:m300\]. They are taken directly from the VLII data at z=0. Both reionisation scenarios reproduce the range of values of ${\rm M_{300}}$ and $\rm{M_V}$ derived by [@strigari2008]. The distribution of ${\rm M_{300}}$ is almost flat, but there is a slight ${\rm M_{300}} - \rm{M_V}$ correlation which is not clearly seen in the data. This is typical of this class of models and is seen in all published works so far ([@munoz2009; @maccio2010; @font2011];K09;B10), as noted by [@kroupa2010]. Eventhough the physical meaning of this relation is currently unclear [@walker2010], we include it as a comparison with previously published models. The bottom panel of Fig. \[f:m300\] shows the total dark matter masses of the satellites: despite the apparent common mass scale of ${\rm M_{300}} \approx 10^7 {\ensuremath{\mathrm{M}_{\odot}}}$, the satellites do indeed span a wide range in [*total*]{} mass, with a large scatter. As for the LF, the masses of the satellites do not allow discrimination between external and internal reionisation.

Radial distribution of satellites
---------------------------------

In order to quantify the impact of the limited sky coverage of the SDSS (about 1/5 in area), we produced mock surveys using our model satellite population: we include the 11 brightest to simulate the “classical” satellites and randomly pick 10 of the fainter ones within the detectable subsample (i.e. the ones that make up the LFs of Fig. \[f:LF\]). We show the mean (solid lines) and dispersion (shaded areas) of the cumulative radial distributions at z=0 obtained for 1000 of these mock surveys in Fig. \[f:crd\]. There is a strong difference in the predicted profiles of our two reionisation scenarios: the distribution for the internal reionisation model is shifted outwards by about 0.3 dex with respect to the external reionisation model. This shift is 2 times larger than the dispersion of the predicted profiles, and is therefore significant. This is the signature of internal reionisation: in the galacto-centric reionisation scenario the satellites of the inner halo see a more intense UV flux than their outer halo cousins, and therefore evaporate faster. Thus the outer halo satellites experience a longer star formation activity period and end up brighter. In this model, all the detected satellites beyond R$>$150 kpc see their photo-evaporation achieved later than z=6, the epoch of end of reionisation of the IGM derived by [@Fan2006].

[![Cumulative radial distribution profiles, observed (black), and predicted by the internal and external reionisation models at z=0. The shaded areas show the dispersion of profiles obtained by 1000 realisations of a mock DR5 survey. For better readability, the blue area is plotted on top of the red one, thereby partially masking it at the top and bottom of the figure. []{data-label="f:crd"}](CRDs_with_errorbars_zfirststar=15.ps "fig:"){width="0.99\linewidth"}]{}

Conclusions and discussion {#s:conclusion}
==========================

We developed a SAM of the satellite population of the MW based on the VLII simulation, simple recipes for star formation and radiative effects, and taking into account the observational biases. Our treatment of the propagation of the I-front within the early MW progenitor and photo-evaporation allows us to explore the effect of the structure of the UV field during reionisation on the properties of the satellite population. In the galacto-centric reionisation model, the satellites of the inner halo are exposed earlier to a nearer UV source, and therefore evaporate faster than the outer halo satellites. As a consequence star formation in the inner halo is suppressed sooner, whereas in the outer halo, satellites continue forming stars until up to several 100 Myr later. The resulting radial distribution of satellites is shifted outwards by 0.3 dex in radius with respect to the external reionisation model, and provides a very good match to the observations. The sensitivity of the radial distribution of satellites to the structure of the UV field during reionisation makes it an interesting probe for the study of reionisation at galaxy scale in the MW with Pan-STARRS and upcoming LSST. Provided that the VLII halo is representative of the MW, our result suggests that star formation at the center of our Galaxy is responsible for the photo-evaporation of its own satellites and the intervening gas. However this interpretation is complicated by cosmic variance in the assembly history of the MW, which can also shift the distribution, as shows Fig. 14 of K09. In order to assess the significance of our result for the reionisation history of the MW, we will need to couple our model to a larger set of high resolution simulations of the formation of the MW similar to the VLII, such as the AQUARIUS dataset [@aquarius], in the spirit of [@font2011; @lunnan2011].

Moreover, a number of simplifications were made in this work, which need to be improved on. Firstly, the [@iliev2005] photo-evaporation simulations were designed to be mostly reliable around $10^7 {\ensuremath{\mathrm{M}_{\odot}}}$. Thus, we need to check via further numerical simulations that the evaporation times we extrapolated up to the $>10^{8} {\ensuremath{\mathrm{M}_{\odot}}}$ regime remain realistic. We will also investigate the effect of endogenous reionisation (i.e. induced by stars formed within the dwarf galaxy itself), which has been neglected in this work. Secondly, we need to generalize the present framework to account for multiple sources in the internal scenario: some fraction of the less massive progenitor haloes may also contribute to the UV field. Even though they are clustered around the most massive progenitor, they could lead to a more patchy UV field around the proto-MW, as seen at larger scale in [@zahn2007]. The transient nature of the radiation field resulting from the short lifetime of massive stars is another untested aspect of the model. Accurate evaluation and modelling of these effects is crucial for the interpretation of Pan-STARRS and LSST surveys of the satellite galaxies of the MW. Since the latter are major progenitors of MW stellar halo and streams, as illustrated in [@gomez2010], accurate interpretation of Gaia data will render this detailed understanding of the reionisation of dwarf galaxies mandatory.

Acknowledgements {#acknowledgements .unnumbered}
================

PO is grateful to B. Famaey, C. Boily, R. Ibata, M. Mouhcine, J. Chardin for stimulating discussions in the course of this study. This work was carried out in the context of the LIDAU (LIght in the Dark Ages of the Universe) project, which is funded by the french ANR (Agence Nationale pour la Recherche, project code ANR-09-BLAN-0030). Additional thanks go to all LIDAU members and workshops participants for providing useful criticisms of the model and a motivating work atmosphere, to M. Busha for communicating details of his models, and J. Diemand for making the Via Lactea II data public. PO also thanks D. Munro for freely distributing his Yorick programming language (available at `http://yorick.sourceforge.net/`), which was used throughout these investigations.

\[lastpage\]
