Gas depletion in Local Group dwarfs on 250 kpc scales: Ram pressure stripping assisted by internal heating at early times
A recent survey of the Galaxy and M31 reveals that more than of dwarf galaxies within kpc of their host galaxy are deficient in H gas. At such an extreme radius, the coronal halo gas is an order of magnitude too low to remove H gas through ram-pressure stripping for any reasonable orbit distribution. However, all dwarfs are known to have an ancient stellar population ( Gyr) from early epochs of vigorous star formation which, through heating of H, could allow the hot halo to remove this gas. Our model looks at the evolution of these dwarf galaxies analytically as the host-galaxy dark matter halo and coronal halo gas builds up over cosmic time. The dwarf galaxies—treated as spherically symmetric, smooth distributions of dark matter and gas—experience early star formation, which sufficiently heats the gas allowing it to be removed easily through tidal stripping by the host galaxy, or ram-pressure stripping by a tenuous hot halo (n cm at kpc). This model of evolution is able to explain the observed radial distribution of gas-deficient and gas-rich dwarfs around the Galaxy and M31 if the dwarfs fell in at high redshifts (–).
In recent years, the Sloan Digital Sky Survey (York et al., 2000) has more than doubled the number of known dwarf galaxies that orbit the Galaxy (Willman et al., 2005a, b; Zucker et al., 2006b; Belokurov et al., 2006; Zucker et al., 2006a; Belokurov et al., 2007; Irwin et al., 2007; Walsh et al., 2007; Belokurov et al., 2010; Willman et al., 2010). These newly discovered dwarf galaxies, being among the most dark-matter dominated structures in the universe, help address the perceived imbalance between the number of predicted low mass dark matter substructure and those observed today, the “Missing Satellites Problem” (Klypin et al., 1999; Moore et al., 1999). As low mass substructures are the building blocks upon which galaxies are formed under the CDM paradigm, the ability to probe dwarf galaxies provides invaluable knowledge into the history and formation of the Local Group (Tolstoy et al., 2009; Karlsson et al., 2011).
Although the dark matter component can be tracked through the use of high resolution N-body simulations (Diemand et al., 2007, 2008; Springel et al., 2008; Wetzel, 2010), a complete model of the evolution of these dwarfs is still elusive. All of the known dwarfs show signs of old stellar populations (Tolstoy et al., 2009) with more of the star formation occurring in discrete bursty periods of star formation—of order of total star formation (Lee et al., 2009)—compared to more massive galaxies which experience roughly continuous star formation. This bursty behaviour, with bursts separated by gigayears, does not necessarily require interactions to trigger star formation (Brosch et al., 2004) with the blow-out and subsequent infall of neutral gas sufficient to create bursts in isolated dwarfs (Valcke et al., 2008; Quillen & Bland-Hawthorn, 2008).
In a recent development, semi-analytic models of galaxy formation, combined with results from high-resolution dark matter simulations, have been used to model the physical properties of a Galactic dwarf galaxy system (Li et al., 2010). These models however do not yet account for the underabundance of H detected in these dwarf galaxies over the last years (Einasto et al., 1974; Grcevich & Putman, 2009, hereafter GP09). The dependence of the H deficiency on galactocentric radius is evidence of significant tidal and/or ram-pressure stripping of the dwarf galaxies which GP09 attribute towards close-in, potentially highly eccentric orbits, allowing the H to be removed by a hot halo surrounding the Galaxy. This removal for close pericentres ( kpc) has been explained by a combination of tidal and ram pressure forces (Mayer et al., 2006).
But we now show that unassisted stripping fails by an order of magnitude to explain the phenomenon of dwarf galaxy depletion, on scales of kpc as observed today, for any reasonable orbit families (§2). In §3, we include the effects of early star formation and find that feedback-assisted stripping is essential to explain the observed effect. Using models of dwarf infall, we then compare results from our model to the observed properties of the H abundance in dwarfs today in §4. Finally, a discussion of the limitations of this model and what impact that will have on the results is given in §5.
2 Dwarf Galaxy Model
A dwarf galaxy with an eccentric orbit will experience a wide range of extremes as it travels along its orbit, from a low velocity through sparse gas at its apogalacticon to a large velocity through dense gas with a strong external radiation field at its perigalacticon. All these environments must be modelled. We describe this wide range of environments by considering four scenarios: the orbital path of the dwarf around the host galaxy, discussed in §2.1; the heating and cooling of the gas present, discussed in §2.2; how tidal or ram pressure stripping removes gas from the dwarf in §2.3; and how this stripping can be influenced by stellar formation and feedback, detailed in §2.4. All four, to varying degrees, will be influenced by how much gas the dwarf galaxy has at the beginning of its orbit and the properties of the dark matter halo that envelopes the dwarf.
The profile of the dwarf galaxy dark-matter halo is not a simple choice, with no solid consensus on which profile most accurately describes the low-mass subhalos present in dwarf galaxies. Observationally, the dark-matter profile of dwarfs cannot be distinguished between a cuspy profile and that of a cored, constant-density profile (Walker et al., 2009). Both of these profiles also arise through simulations with the cuspy Einasto profile arising through large N-body, dark-matter only, simulations (Springel et al., 2008). With the addition of baryons, it has been argued that an even steeper central core—through the process of baryonic contraction—can be present than predicted by the dissipationless -body simulations (Blumenthal et al., 1986; Del Popolo, 2009; Napolitano et al., 2010). However, by coupling the baryons with the dark matter, and allowing dynamical and angular momentum transfer to occur, constant density cores with Burkert like profiles (Burkert, 1995) are able to develop in simulation (El-Zant et al., 2001; Del Popolo, 2009). Since our major focus is on dark matter-halos with low gas fraction, whereby both baryon–dark-matter coupling and baryonic contraction will have a smaller impact, we ignore these competing factors and use the Einasto profile.
The choice of profile for the host galaxy is less contentious, and is chosen to have an evolving Einasto profile. Although the host contains more baryons than the dwarf galaxies we are interested in, due to its much larger mass baryonic contraction processes have a lower net effect. Thus the evolution of the host galaxy is dominated by the infall of dwarfs onto the galaxy. The evolution of the host of this galaxy is assumed to be similar to that of the Milky Way, with a growth rate described by McBride et al. (2009) leading to a final virial mass of M (for an Einasto model with the same circular velocity at virial radius as the NFW model in Smith et al., 2007), to give a mass at redshift of
We assume that any cold gas is of a constant density, and gas in the warm and hot phases in the dwarf is in hydrostatic equilibrium with the density of the gas then given by
with the scaled radius, the halo scale radius, the characteristic scale velocity of the dark matter halo, the Einasto free parameter, and the lower incomplete gamma function.
The central density of the gas is calculated by assuming the calculated central density of Leo T in GP09 of cm is comprised of gas in pressure equilibrium with a mass ratio of cold gas to warm gas, matching the total ratio of masses within Leo T from which this value was calculated (Ryan-Weber et al., 2008). This gives a warm gas density of cm and cold gas density of cm. The hot phase central density is assumed to be in pressure equilibrium with the cold and warm phases.
2.1 Orbit Of The Dwarf
The orbit of the dwarf galaxy within the larger host galaxy is given by the shape of the potential. As the host galaxy evolves, the halo scale radius and scale velocity will change according to (Sternberg et al., 2002)
where is the characteristic overdensity of dark matter halos, is the mean baryonic density of the universe, is the concentration parameter and is the mass profile as given in Nichols & Bland-Hawthorn (2009).
The characteristic overdensity is given by (Bryan & Norman, 1998)
and the mean baryonic density by
For an Einasto profile, the concentration parameter, , also referred to as , is (Duffy et al., 2008)
In this growing potential there is no angular dependence and therefore the angular momentum is conserved; such that
where is the per unit mass angular momentum squared and and are the peri and apogalacticon radii that would occur if the potential was static.
We present our results in terms of both orbit circularity and eccentricity in recognition of the fact that both are in common use when discussing CDM model. We define the eccentricity of the orbit in terms of the ratio of the apogalacticon to the perigalacticon, , and the corresponding circularity to be . These parameters are not conserved throughout the orbit due to the growing potential, and when referred to we use the final state of these parameters.
We also add dynamical friction as a non-conservative force, with a magnitude of the instantaneous dynamical friction timescale multiplied by the velocity (Zhao, 2004)
giving the equations of motion as
with the additional term in the radial direction arising from the background contribution to the density of the dark matter, which becomes non-negligible beyond .
2.2 Heating And Cooling
Throughout the passage of the orbit, the gas present within a dwarf galaxy will be exposed to a variety of radiation fields, with the most powerful being the radiation field from the host galaxy, the extragalactic background light and the stellar radiation field produced from high mass, short lived stars formed within the dwarf galaxy. To examine if the changing orbit itself is enough to produce the observed galactocentric effects, we first examine the system without the internal stellar radiation field.
Without internal star formation, two contributions to the radiation field will dominate the heating and ionization of the dwarf galaxy, namely the extragalactic background field and the radiation field from the host galaxy which will dominate when the dwarf is within kpc of the host galaxy’s centre. For the extragalactic background field we use the redshift-dependent fields of Faucher-Giguère et al. (2009) and assume that all the radiation enters the dwarf galaxy radially. To a good approximation, the radiation field of the host galaxy is assumed to be a delta function at eV with the number of photons impacting the dwarf galaxy at a distance from the disk given by (Bland-Hawthorn & Maloney, 1999)
with the factor arising from a escape fraction vertically from the disk and isotropized over the halo. SFR is the star formation rate at redshift , given by (Just & Jahreiß, 2010)
where is the time since disk formation, which we assume formed Gyr before the present day.
These two radiation fields combine to heat and ionize the cold gas that is present in the dwarf galaxy. The gas is ionized according to the absorbed radiation field, for which a radial radiation field yields
with the radial size of the cold gas in the dwarf and is the combined radiation field of the Galaxy and the extragalactic UV field at a distance , frequency and redshift in ergs s cm Hz.
These radiation fields will also heat the gas, with the heating (considering only H and He) given by (Wolfire et al., 1995)
where is the heating per primary electron given by Wolfire et al. (1995).
As more particles are ionized than can be warmed, any excess ionizations (after recombination) are converted to heating with all eV going to heat the gas. The transfer from cold neutral gas to warm ionized gas, is then given by
where and are the temperatures of the cold and warm phase respectively, and is the central density and scaled radii of the edge of the warm phase respectively, and is the Einasto gas distribution function. Correspondingly .
Without any other source of heating (see §2.4), the cooling prevents any significant heating in the model, keeping most gas as cold gas with only a small amount of warm gas in the centre.
2.3 Tidal and Ram Pressure Stripping
Gas in the warm or cold phase may be stripped through a combination of tidal forces and that of ram pressure arising from the orbit through the hot halo of the host galaxy.
The tidal forces will not only potentially remove gas from the dwarf galaxy but will strip down the dark matter from the halo that surrounds the dwarf. The tidal extent of the dwarf’s halo, neglecting the contribution of baryonic mass, is given by (Hayashi et al., 2003)
where is the tidal radius of the dwarf, and is the orbital radius of the dwarf from the host galaxies centre. For an Einasto profile, this tends not to have an analytic solution, and for reasons of speed is only solved initially, before using the derivative to calculate the change in tidal radius over time. Any gas that extends above the tidal radius is considered to be no longer bound, and hence lost, in order to simplify calculation. A more extensive treatment involving re-accretion onto dwarfs requires hydrodynamical treatments that are now under way. Within our models, tidal stripping does not produce significant gas loss with Galactic and extragalactic ionizing fields keeping most gas deep within the tidal radius.
Working at the same time as the tidal stripping is ram pressure stripping. Here we use the McCarthy et al. (2008) formalism whereby the ram pressure removes shells of gas from the dwarf galaxy, assuming that the gas distribution behaves as an isothermal sphere beyond the cloud edge if
This formalism, originally used for a hot halo, can be carried over to the warm and cold phases by assuming spherical symmetry and a density drop off outside the cloud edge.
We note that when using the realistic approximation , this is times more difficult to strip than the Gunn & Gott (1972) equation used by GP09, i.e.
If the stripping does occur, then the rate of stripping will be dependent upon the speed of the shock induced in the gas (Mori & Burkert, 2000), such that
Here we assume that the halo stays spherically symmetric at all times, and hence the rate of the change in radius will be half the speed of the shock, i.e. .
The density of the host galaxy halo will be dependent upon the temperature of the halo. We use an isothermal halo with the temperature set at the virial temperature of the halo at any given redshift. We set an unchanging central density to give a density of cm at kpc at the present day. These values are consistent with isothermal halo parameters from Battaglia et al. (2005), the values used by Bland-Hawthorn et al. (2007) (which uses K, cm at kpc), and the high-entropy Galactic halo model of Kaufmann et al. (2009).
For this stripping, the phases are assumed to be independent, and the cold phase can be removed at the same time as the warm phase if the ram pressure is large enough. The rate of change of mass is then simply
With only Galactic and extragalactic ionizing fields, ram-pressure, like tidal stripping, does not produce large changes in gas mass even over cosmic time. Even with more favourable conditions to ram pressure stripping, such as the Gunn-Gott criterion and a lower temperature halo (resulting in a denser inner portion), only a small amount of gas is lost—less than M—corresponding to the cold gas heated by the Galactic and extragalactic UV fields to a warm phase and then lost to ram-pressure stripping.
These results are in conflict with much more substantial gas loss has been measured in hydrodynamical simulations with ram pressure alone able to strip the gas from close in ( kpc), low mass halos (Mayer et al., 2006). However, at larger perigalacticons, dwarfs in Mayer et al. (2006) are still able to retain much of their centrally concentrated gas, suggesting another form of heating is required to allow the removal of gas from dwarfs with larger perigalacticons (e.g. the perigalacticon of Sextans and Draco is expected to be kpc Lux et al., 2010).
2.4 Star Formation and Feedback
A dwarf consisting of cold, dense gas is resistant to tidal and ram-pressure stripping with only a thin skin, ionized by the Galactic of extragalactic UV fields being removed.
Early star formation however, will heat this cold gas, raising it in the potential well and making it more easily stripped. This warm ionized gas, at the same pressure as cold gas will occupy times the volume, and being much less dense will become even easier to strip than its height in the potential well would indicate.
The star formation in dwarfs is considered to consist of periods of low-level star formation, during which short bursts are induced which increase the star formation by a factor of , consistent with dwarfs observed by Lee et al. (2009). We assume these bursts are triggered by perigalacticon passages, induced by shocks created through tidal interactions (Pasetto et al., 2010). There is also evidence that bursts may be triggered by the re-accretion of heated and expanded gas (Valcke et al., 2008), however, much of this gas will be stripped while in the warm phase and far from the centre preventing the re-accretion and subsequent starburst from occurring.
The base star formation rate is taken to be similar to the dwarfs that surround M31, with a star formation rate (Kaisin & Karachentsev, 2006)
where is a constant calculated by assuming that the gas will be completely depleted at a time Gyr. For an initial gas mass of M, .
The UV–X-ray spectrum produced by the star formation in the dwarfs—which is responsible for the majority of heating—was calculated with Starburst99 (Leitherer et al., 1999); when a burst occured the change in the spectrum was calculated by summing many short bursts together to produce an approximately continuous change in the spectrum. Above the Lyman limit, this spectrum (in ergs s cm Hz sr) and its changes during a burst was well fitted by
where changes linearly over a Myr period at the beginning and end of each burst from (ergs s cm Hz sr) at the low level continuous star formation rate to (ergs s cm Hz sr) during a burst.
The addition of star formation will also result in supernova superheating the gas well above the temperature of the warm phase Any power from the supernova is assumed to go into heating gas into a hot ( K) phase. In reality, the energy from the supernova will also go into expanding the cold and warm phases, making them more susceptible to stripping, but due to the limitations of a smooth density profile we assume this energy is instead used to heat the gas.
We find the power output of a supernova for continuous, non-bursty star formation is given by
After a burst has commenced, the power output by the supernovae increases over Myr, with the power over this period of increasing output given by
where is the time in Myr since the burst commenced.
During the increased star formation associated with the bursty period (lasting yr Sharp & Bland-Hawthorn, 2010) the power output is
After the burst has ended and star formation drops back down to the continuous star formation rate, the power output subsides as massive stars die off. Over a period of Myr, the power output declines according to
where is the time since the starburst has ended in Myr.
In this simplified model, we assume that any power from the supernova is dissipated in the gas and goes into heating the cold and warm gas at efficiency and split between the cold and warm phase in proportion to their mass. The mass loss rate will then be the power dissipated into that phase divided between the energy difference between that phase and the hot phase.
Although the efficiency of supernova heating is typically much lower, there is only a small amount of material that is heated into the hot phase, with much more material being lost through the warm phase. In this limit, the snowplough phase of supernova expansion, where a large portion of energy goes, may have a noticeable contribution by raising gas outside of the potential well and we keep the efficiency at unity to partially compensate for this effect.
Under these conditions the internal star formation greatly exceeds the extragalactic UV background at early times, making it the dominant source of heating. The cooling of gas from the warm phase to the cold phase can approach the level of UV–X-ray heating, however, the amount of warm gas that is able to be held on by the halo is limited, preventing most gas from recooling. The relative rates of heating by the combined radiation fields (internal star formation, Galactic and extragalactic) and the cooling rate is shown in Fig. 1 for a model with perigalacticon of kpc and circularity of (eccentricity of ). Very little gas is moved to the hot phase from supernova bursts, much more gas is heated to the warm phase which quickly extends beyond the tidal radius, resulting in large amounts of gas loss from tidal stripping. Only in more massive dwarfs with correspondingly larger central densities will sufficient quantities of warm gas be able to cool back to a cold phase to allow for significant supernova heating.
2.5 Method and Initial Conditions
The gas and dark matter masses were tracked by solving the time dependent equations—equations (1)–(28)—with the odeint routine (Press et al., 1992) and encoded in python. Once a gas phase dropped below M, it was no longer tracked and assumed to no longer have any gas in that phase.
The dwarfs were initiated with M for the cold phase and M for the warm and hot phases. The dark matter halo of each dwarf galaxy was set with a virial mass of M for models beginning at and and a virial mass of M for . For all starting points, the predicted mass within pc is about M, consistent with the mass inferred for dwarf galaxies (Strigari et al., 2008).
The orbit of the dwarf galaxy was set so that it ended at apogalacticon, with the initial conditions found by integrating backwards without dynamical friction initially, and then with dynamical friction along the previous orbit until convergence of at all radii was found. The model was run for dwarf galaxies that at their present day orbits will have a pericentre of between and kpc, and a circularity of between and . We assume each point is similar to those around it, and each model was at the centre of a kpc, circularity box, for those with a pericentre kpc, a kpc, circularity box for pericentre kpc and a kpc, circularity box for kpc, for a total range of – kpc and – in circularity. A sample of these orbits (Fig. 2) illustrates the effect of varying eccentricity and pericentre in a growing host halo with dynamical friction.
Beginning at – allows most halos to retain a large portion of their cold gas, while beginning at higher redshift means most halos have been stripped of nearly all their cold gas. In all cases halos lose over of gas mass to warm phase stripping resulting from early heating by internal star formation. As could be expected, highly circular orbits at low pericentre show the strongest gas loss, with the final mass displayed in Fig. 3. This is large gas loss attributable to an on-average larger number of bursts and the increased halo density at low pericentres. We note that finishing at apogalacticon will always result in later bursts for the same pericentre, but the results remain if the models begin at apogalacticon in which case the bursts will always be earlier.
A large portion of the parameter space exists in which a halo will always be within kpc and retain over M of gas. Numerous other halos will spend a portion of their orbit within this kpc limit, but as the model finishes at apogalacticon, these subhalos may have experienced less pericentres than would be expected if they were within kpc today.
For those subhalos which begin at (Fig. 4) no halos are always within kpc and maintain over M of H; although some halos should be detectable with over M of gas remaining that will always be within kpc.
The low circularity, high pericentre drop in gas mass is a result of the model finishing at apogalacticon. This effect is much more pronounced for halos beginning at displayed in Fig. 5, here the initial velocities may exceed km s. The long time in orbit, and lower virial mass of the dwarf galaxies in this run means that nearly all halos that retain gas will not have an apogalacticon within kpc, that is, any observed halos with gas that began at will only be observed within kpc for a portion of their orbit.
In all models, at early redshifts the high rate of star formation produces large quantities of warm gas which is removed by tidal forces existing on the dwarf galaxy, and warm gas at the edge of the tidal radius is then quickly removed by stripping of the hot halo. This initial period of extremely high mass loss accounts for over of all mass loss as seen in Fig. 6. This initial period of gas loss is followed by a steady state, seen more easily in the mass loss rate versus time in Fig. 7, is reached where cold gas is not warmed rapidly, and a small amount of warm gas can survive stripping, with most of the gas being stripped from dwarfs which are approaching their perigalacticon and experiencing bursts of star formation and the increased impact of the Galactic UV field.
At low redshifts, the rate of stripping plummets as the initial conditions require that dwarf galaxies be at their apogalacticon at . Distributing the dwarf galaxies in phase along their orbits today has a minimal effect on the results presented in Figs. 3–5. The steep rate of gas loss for dwarf galaxies beginning at is also a byproduct of our requirement that dwarfs finish at apogalacticon today. This leads to the excluded region (marked ‘RAD’ in Fig. 5). The reason is that the dwarf is forced to orbit at unrealistic velocities to get back to its apogalacticon within an evolving potential well. This artifact has only a minimal effect on our results.
4 Comparison with observations and completeness
We now determine the fraction of dwarf galaxies retaining gas within kpc of the Milky Way and M31. In order to estimate this, the completeness of the parameter space searched needs to be calculated. We assume that the infalling halos follow a distribution identical to that of Wetzel (2010), which used a dissipationless N-body simulation with CDM cosmology (, , , and ) with a particle mass of M. We note that although dwarf galaxy sized halos cannot be represented in this model at we extend the distribution on the basis that for higher masses there exists no dependence on mass. For a Milky Way sized halo at these distributions are
Taking the circularity and pericentre to be independent, the halos in the parameter space searched will comprise of all halos with a pericentre kpc and of all halos, with the distribution within the parameter space shown in Fig. 8. We split the distribution into four main regions: Regions A and B consists of those halos which were within our parameter space and always/sometimes within kpc respectively. Regions C and D (halos with circularitys ) consists of those halos that were outside our parameter space and always/sometimes within kpc respectively. Regions C and D were excluded from the parameter space as the large number of pericentres and low tidal gradients would give an unrealistic amount of bursts as well as the large number of pericentres being computationally expensive. The shaded region of low circularity/high eccentricity consists of halos which were outside our parameter space and spend so little time within kpc that they were ignored. The shaded region of low pericentre is halos that would be stripped due to their very close passages, but the no disk assumption begins to break down as well as the increased computational cost means these halos were not modelled.
In our model all dwarfs end at apogalacticon, however, the results are substantially the same if you begin them at apogalaction (for low redshifts where the orbit does not change significantly over time), we hence assume for analysis that they are at a random point along their orbit, and that this orbit is elliptical, as opposed to the standard Rosetta orbit an Einasto profile will produce. Under these assumptions, a dwarf galaxy with a pericentre of and eccentricity will be within kpc with probability
Within each point in the model we then weight it by its contribution to all halos, with the weight then
where is half the size of the box in kpc, and is half the size of the box in circularity.
We assume that all halos with a pericentre kpc cannot retain gas, and that subhalos with a circularity —halos from regions C and D—will have the same amount of gas, as halos with the same pericentre and a circularity of . Halos with a circularity were assumed to contain gas, however, these halos have very eccentric orbits and hence will normally have a near zero weight.
We then calculate what fraction of observed dwarf galaxies contain gas, looking at all galaxies below a given galactocentric radius in kpc bins out to kpc and compare this to the fraction of galaxies—with confirmed detections and excluding the SMC and LMC discussed below—that contain gas from GP09. We use three cutoffs to distinguish between gas-deficient and gas-rich dwarfs. A lower cutoff of M is above the upper-boundary of a number of GP09 galaxies A mid-range cut-off of M is above the upper-boundary of even more galaxies in the GP09 sample. Finally an upper cut-off of M, this cut-off is below the H mass of any confirmed gas-rich galaxy in the GP09 sample. These cut offs do not accurately represent the physical constraints on observations, with the observed being a combination of low cut-offs for nearby Galactic satellites to higher-cutoffs for far-off or M31 satellites.
Under our model, the LMC and SMC being much more massive than the typical dwarf— M (Schommer et al., 1992), M (Harris & Zaritsky, 2006)—will be able to retain significant quantities of warm gas against ram-pressure stripping. This protection from stripping would allow a significant quantity of gas to survive until the present day, with only tidal interactions—potentially those between the LMC and SMC—removing material from either dwarf galaxy. Although tidal stripping is included in our model, it is too simplistic to model the multi-body system present in the Magellanic clouds, with hydrodynamical simulations required to simulate the LMC-SMC-Galaxy system (Lin et al., 1995; Bland-Hawthorn et al., 2007; Nidever et al., 2008; Besla et al., 2010).
For models beginning at (Fig. 9) we see a large over estimation for all gas cut-off masses. For models beginning at (Fig. 10) again the predicted fraction exceeds the observed fraction of dwarfs, but predicts a lower fraction than models beginning at . Discovery of new gas-rich dwarfs far from the Galaxy may raise the observed fraction, but this may not explain the over-abundance predicted close in, where gas-rich dwarfs would be more detectable than further gas-deficient dwarfs due to more ongoing star formation.
For dwarf galaxies that begin their orbits at (Fig. 11), the M and M cut-offs again overestimates the fraction of halos which will contain gas. The M cut-off closely follows the shape of the GP09 sample, and is always below this fraction. At any distance we would expect the Grcevich & Putman (2009) fraction to be slightly higher than reality, as the gas poor dwarfs being difficult to see may not be fully accounted for. That the toy model represents the gas fraction best for halos beginning at adds support to the idea that the Milky Way dwarfs were formed at or before this redshift (Lux et al., 2010).
5 Summary & Conclusion
This simple toy model reproduces the observed fraction of gas-rich dwarfs, however, there are several factors that were assumed to be negligible which could affect the amount of gas that survives. The accretion of gas onto dwarfs is unaccounted for in this model. The prospect of low redshift accretion (Ricotti, 2009) in particular would greatly increase the chance of a dwarf surviving with gas to the present day, and if included for all dwarfs would likely result in an overestimation of the gas fraction even for dwarfs beginning at late redshifts. The effects of dust through photoelectric heating and cooling is also ignored. Due to the low metallicity environments that dwarf galaxies typically have, this effect will be smaller than in larger galaxies, but may still be an important source of heating or cooling. We also assume that the early extragalactic UV field is uniform in space, a clumpy radiation field around the time of reionization, may have a large impact on the amount of gas that remains cold and protected from ram pressure stripping, this early gas loss could greatly impact the survival around the earliest pericentre passages. The use of a smooth medium for the gas—compared to a more realistic fractal medium—minimises the cooling of warm gas and allows it to extend beyond the tidal radius. Even with a large filling factor, the gas will cool much more quickly via metal line cooling than in the smooth medium minimising the large gas loss at the beginning. Potentially the biggest limitation is the assumption that all of the supernova energy goes into heating the gas to an extremely hot state, much of this energy likely goes into raising cold and warm gas out of the potential well of the dwarf, allowing it to be much more easily stripped. In particular this will predominantly impact the lowest gas masses, where the smooth gas assumption versus a fractal medium is more likely to have a large impact (Bland-Hawthorn et al., 2007).
Due to the large initial gas loss, our star formation rates do not represent that of the majority of dwarfs, with most stars forming early in the dwarfs life, as opposed to a roughly continuous star formation rate with some bursts (Weisz et al., 2011). This makes our model more suited to explaining dwarfs with these early periods of star formation where a majority of stars are formed, e.g. BK5N, KDG52 (Weisz et al., 2011).
Even with these limitations, we believe that the toy model provides strong support towards internal heating due to early star formation, allowing the gas to more easily stripped due to the larger scale heights and lower density of warm and hot gas. In particular we were able to reproduce the fraction of dwarf galaxies that retain H assuming that the dwarf galaxies infall at a redshift of consistent with Lux et al. (2010).
The covering fraction of warm gas stripped from the dwarfs is expected to be low, with the contrails containing approximately the same amount of warm gas as exists inside the cold streams that cover of the projected area of galaxies (Faucher-Giguere & Keres, 2010). This gas breaking away in small clumps of warm material will form a warm component of the hot halo which is stabilised against heat conduction from the halo by cooling (Vieser & Hensler, 2007b, a) before falling onto the Galaxy as a warm rain. Even with the majority of dwarf galaxies entering at early redshifts, the large timescale required to complete large pericentre orbits ( Gyr for kpc) means that gas from any late infalling dwarfs—such as those from or —will only be falling onto the disk of the Galaxy today. A search for this gas in H around local galaxies is now being undertaken using the Maryland-Magellan Tunable Filter (Veilleux et al., 2010) and the Grantecan Osiris Tunable Filter (Cepa et al., 2003) where the expected emission measure of the gas ( cm pc) may be just large enough to be detectable.
- Battaglia et al. (2005) Battaglia, G., et al. 2005, MNRAS, 364, 433
- Belokurov et al. (2006) Belokurov, V., et al. 2006, ApJ, 647, L111
- Belokurov et al. (2007) —. 2007, ApJ, 654, 897
- Belokurov et al. (2010) —. 2010, ApJ, 712, L103
- Besla et al. (2010) Besla, G., Kallivayalil, N., Hernquist, L., van der Marel, R. P., Cox, T. J., & Kereš, D. 2010, ApJ, 721, L97
- Bland-Hawthorn & Maloney (1999) Bland-Hawthorn, J., & Maloney, P. R. 1999, ApJ, 510, L33
- Bland-Hawthorn et al. (2007) Bland-Hawthorn, J., Sutherland, R., Agertz, O., & Moore, B. 2007, ApJ, 670, L109
- Bland-Hawthorn et al. (2011) Bland-Hawthorn, J., Sutherland, R., & Karlsson, T. 2011, ApJ, submitted
- Blumenthal et al. (1986) Blumenthal, G. R., Faber, S. M., Flores, R., & Primack, J. R. 1986, ApJ, 301, 27
- Brosch et al. (2004) Brosch, N., Almoznino, E., & Heller, A. B. 2004, MNRAS, 349, 357
- Bryan & Norman (1998) Bryan, G. L., & Norman, M. L. 1998, ApJ, 495, 80
- Burkert (1995) Burkert, A. 1995, ApJ, 447, L25+
- Cardone et al. (2005) Cardone, V. F., Piedipalumbo, E., & Tortora, C. 2005, MNRAS, 358, 1325
- Cepa et al. (2003) Cepa, J., et al. 2003, in Presented at the Society of Photo-Optical Instrumentation Engineers (SPIE) Conference, Vol. 4841, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, ed. M. Iye & A. F. M. Moorwood, 1739–1749
- Dalgarno & McCray (1972) Dalgarno, A., & McCray, R. A. 1972, ARA&A, 10, 375
- Del Popolo (2009) Del Popolo, A. 2009, ApJ, 698, 2093
- Diemand et al. (2007) Diemand, J., Kuhlen, M., & Madau, P. 2007, ApJ, 657, 262
- Diemand et al. (2008) Diemand, J., Kuhlen, M., Madau, P., Zemp, M., Moore, B., Potter, D., & Stadel, J. 2008, Nature, 454, 735
- Duffy et al. (2008) Duffy, A. R., Schaye, J., Kay, S. T., & Dalla Vecchia, C. 2008, MNRAS, 390, L64
- Einasto et al. (1974) Einasto, J., Saar, E., Kaasik, A., & Chernin, A. D. 1974, Nature, 252, 111
- El-Zant et al. (2001) El-Zant, A., Shlosman, I., & Hoffman, Y. 2001, ApJ, 560, 636
- Faucher-Giguere & Keres (2010) Faucher-Giguere, C., & Keres, D. 2010, ArXiv e-prints
- Faucher-Giguère et al. (2009) Faucher-Giguère, C., Lidz, A., Zaldarriaga, M., & Hernquist, L. 2009, ApJ, 703, 1416
- Grcevich & Putman (2009) Grcevich, J., & Putman, M. E. 2009, ApJ, 696, 385
- Gunn & Gott (1972) Gunn, J. E., & Gott, III, J. R. 1972, ApJ, 176, 1
- Harris & Zaritsky (2006) Harris, J., & Zaritsky, D. 2006, AJ, 131, 2514
- Hayashi et al. (2003) Hayashi, E., Navarro, J. F., Taylor, J. E., Stadel, J., & Quinn, T. 2003, ApJ, 584, 541
- Irwin et al. (2007) Irwin, M. J., et al. 2007, ApJ, 656, L13
- Just & Jahreiß (2010) Just, A., & Jahreiß, H. 2010, MNRAS, 402, 461
- Kaisin & Karachentsev (2006) Kaisin, S. S., & Karachentsev, I. D. 2006, Astrophysics, 49, 287
- Karlsson et al. (2011) Karlsson, T., Bromm, V., & Bland-Hawthorn, J. 2011, Rev. Mod. Phys., In prep.
- Kaufmann et al. (2009) Kaufmann, T., Bullock, J. S., Maller, A. H., Fang, T., & Wadsley, J. 2009, MNRAS, 396, 191
- Klypin et al. (1999) Klypin, A., Kravtsov, A. V., Valenzuela, O., & Prada, F. 1999, ApJ, 522, 82
- Lee et al. (2009) Lee, J. C., Kennicutt, R. C., José G. Funes, S. J., Sakai, S., & Akiyama, S. 2009, ApJ, 692, 1305
- Leitherer et al. (1999) Leitherer, C., et al. 1999, ApJS, 123, 3
- Li et al. (2010) Li, Y., De Lucia, G., & Helmi, A. 2010, MNRAS, 401, 2036
- Lin et al. (1995) Lin, D. N. C., Jones, B. F., & Klemola, A. R. 1995, ApJ, 439, 652
- Lux et al. (2010) Lux, H., Read, J. I., & Lake, G. 2010, MNRAS, 406, 2312
- Martin (1998) Martin, C. L. 1998, ApJ, 506, 222
- Mayer et al. (2006) Mayer, L., Mastropietro, C., Wadsley, J., Stadel, J., & Moore, B. 2006, MNRAS, 369, 1021
- McBride et al. (2009) McBride, J., Fakhouri, O., & Ma, C. 2009, MNRAS, 398, 1858
- McCarthy et al. (2008) McCarthy, I. G., Frenk, C. S., Font, A. S., Lacey, C. G., Bower, R. G., Mitchell, N. L., Balogh, M. L., & Theuns, T. 2008, MNRAS, 383, 593
- Moore et al. (1999) Moore, B., Ghigna, S., Governato, F., Lake, G., Quinn, T., Stadel, J., & Tozzi, P. 1999, ApJ, 524, L19
- Mori & Burkert (2000) Mori, M., & Burkert, A. 2000, ApJ, 538, 559
- Napolitano et al. (2010) Napolitano, N. R., Romanowsky, A. J., & Tortora, C. 2010, MNRAS, 405, 2351
- Nichols & Bland-Hawthorn (2009) Nichols, M., & Bland-Hawthorn, J. 2009, ApJ, 707, 1642
- Nidever et al. (2008) Nidever, D. L., Majewski, S. R., & Burton, W. B. 2008, ApJ, 679, 432
- Pasetto et al. (2010) Pasetto, S., Grebel, E. K., Berczik, P., Chiosi, C., & Spurzem, R. 2010, ArXiv e-prints
- Press et al. (1992) Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical recipes in FORTRAN. The art of scientific computing
- Quillen & Bland-Hawthorn (2008) Quillen, A. C., & Bland-Hawthorn, J. 2008, MNRAS, 386, 2227
- Ricotti (2009) Ricotti, M. 2009, MNRAS, 392, L45
- Ryan-Weber et al. (2008) Ryan-Weber, E. V., Begum, A., Oosterloo, T., Pal, S., Irwin, M. J., Belokurov, V., Evans, N. W., & Zucker, D. B. 2008, MNRAS, 384, 535
- Schommer et al. (1992) Schommer, R. A., Suntzeff, N. B., Olszewski, E. W., & Harris, H. C. 1992, AJ, 103, 447
- Schure et al. (2009) Schure, K. M., Kosenko, D., Kaastra, J. S., Keppens, R., & Vink, J. 2009, A&A, 508, 751
- Sharp & Bland-Hawthorn (2010) Sharp, R. G., & Bland-Hawthorn, J. 2010, ApJ, 711, 818
- Smith et al. (2007) Smith, M. C., et al. 2007, MNRAS, 379, 755
- Springel et al. (2008) Springel, V., et al. 2008, MNRAS, 391, 1685
- Sternberg et al. (2002) Sternberg, A., McKee, C. F., & Wolfire, M. G. 2002, ApJS, 143, 419
- Strigari et al. (2008) Strigari, L. E., Bullock, J. S., Kaplinghat, M., Simon, J. D., Geha, M., Willman, B., & Walker, M. G. 2008, Nature, 454, 1096
- Tolstoy et al. (2009) Tolstoy, E., Hill, V., & Tosi, M. 2009, ARA&A, 47, 371
- Valcke et al. (2008) Valcke, S., de Rijcke, S., & Dejonghe, H. 2008, MNRAS, 389, 1111
- Veilleux et al. (2010) Veilleux, S., et al. 2010, AJ, 139, 145
- Vieser & Hensler (2007a) Vieser, W., & Hensler, G. 2007a, A&A, 475, 251
- Vieser & Hensler (2007b) —. 2007b, A&A, 472, 141
- Walker et al. (2009) Walker, M. G., Mateo, M., Olszewski, E. W., Peñarrubia, J., Wyn Evans, N., & Gilmore, G. 2009, ApJ, 704, 1274
- Walsh et al. (2007) Walsh, S. M., Jerjen, H., & Willman, B. 2007, ApJ, 662, L83
- Weisz et al. (2011) Weisz, D. R., et al. 2011, ArXiv e-prints
- Wetzel (2010) Wetzel, A. R. 2010, ArXiv e-prints
- Willman et al. (2010) Willman, B., Geha, M., Strader, J., Strigari, L. E., Simon, J. D., Kirby, E., & Warres, A. 2010, ArXiv e-prints
- Willman et al. (2005a) Willman, B., et al. 2005a, AJ, 129, 2692
- Willman et al. (2005b) —. 2005b, ApJ, 626, L85
- Wolfire et al. (1995) Wolfire, M. G., Hollenbach, D., McKee, C. F., Tielens, A. G. G. M., & Bakes, E. L. O. 1995, ApJ, 443, 152
- York et al. (2000) York, D. G., et al. 2000, AJ, 120, 1579
- Zhao (2004) Zhao, H. 2004, MNRAS, 351, 891
- Zucker et al. (2006a) Zucker, D. B., et al. 2006a, ApJ, 650, L41
- Zucker et al. (2006b) —. 2006b, ApJ, 643, L103