Articles

DELAYED STAR FORMATION IN ISOLATED DWARF GALAXIES: HUBBLE SPACE TELESCOPE STAR FORMATION HISTORY OF THE AQUARIUS DWARF IRREGULAR*

, , , , , , and

Published 2014 October 13 © 2014. The American Astronomical Society. All rights reserved.
, , Citation Andrew A. Cole et al 2014 ApJ 795 54 DOI 10.1088/0004-637X/795/1/54

0004-637X/795/1/54

ABSTRACT

We have obtained deep images of the highly isolated (d = 1 Mpc) Aquarius dwarf irregular galaxy (DDO 210) with the Hubble Space Telescope Advanced Camera for Surveys. The resulting color–magnitude diagram (CMD) reaches more than a magnitude below the oldest main-sequence turnoff, allowing us to derive the star formation history (SFH) over the entire lifetime of the galaxy with a timing precision of ≈10% of the lookback time. Using a maximum likelihood fit to the CMD we find that only ≈10% of all star formation in Aquarius took place more than 10 Gyr ago (lookback time equivalent to redshift z ≈ 2). The star formation rate increased dramatically ≈6–8 Gyr ago (z ≈ 0.7–1.1) and then declined until the present time. The only known galaxy with a more extreme confirmed delay in star formation is Leo A, a galaxy of similar MH i/M, dynamical mass, mean metallicity, and degree of isolation. The delayed stellar mass growth in these galaxies does not track the mean dark matter accretion rate from CDM simulations. The similarities between Leo A and Aquarius suggest that if gas is not removed from dwarf galaxies by interactions or feedback, it can linger for several gigayears without cooling in sufficient quantity to form stars efficiently. We discuss possible causes for the delay in star formation including suppression by reionization and late-time mergers. We find reasonable agreement between our measured SFHs and select cosmological simulations of isolated dwarfs. Because star formation and merger processes are both stochastic in nature, delayed star formation in various degrees is predicted to be a characteristic (but not a universal) feature of isolated small galaxies.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

1.1. Environment and the Evolution of Dwarf Galaxies

The question of the early evolution of low-mass galaxies has been a vexing one for galaxy formation theory in a wide range of hierarchical merger models for more than three decades (e.g., Dekel & Silk 1986; Kauffmann et al. 1993; Somerville 2002; Benson et al. 2003; Sawala et al. 2011; Kuhlen et al. 2012). During this timespan a range of problems has been identified and addressed in successive generations of models; these problems tend to relate to the efficiency with which gas is accreted, retained, cooled, and astrated in the lowest-mass dark matter halos.

Recurring issues surrounding the over- or under-production of stars in simulations include overcooling, reionization, the role of feedback, and the effects of tidal perturbations; a recent, detailed discussion of these issues in the context of Local Group dwarf galaxies may be found in Skillman et al. (2014). A variety of recently proposed solutions invoke a range of external heating and internal feedback processes, many of which are strongly mass and/or environment dependent (e.g., Hopkins et al. 2013, and references therein); the influence of numerical resolution on the simulations is also likely to play a role (e.g., Brooks & Zolotov 2014). In addition to supernova feedback, reionization of the Universe at redshifts ≳6 has been considered as a source of heating that might inhibit star formation in the smallest galaxies (e.g., Babul & Rees 1992). In the simplest formulation, reionization may quench star formation in galaxies with a total dark halo mass ≲ 1010M by warming gas to the point where it may be driven out of the shallow potential wells of these systems (e.g., Thoul & Weinberg 1996; Gnedin 2000; Ricotti et al. 2002); a review of the observational evidence for and against such a scenario has been given by Weisz et al. (2014b), who find that there is some evidence for negative feedback on star formation, but with large galaxy-to-galaxy variability and with some inconsistency between the quenching timescale and the epoch of reionization.

While some process (or set of processes) clearly inhibits the formation of dwarf galaxies compared to their proliferation in cold dark matter-only simulations, the nature of these processes is still an open question. Continued observational study of isolated dwarfs is necessary for at least two reasons; one is to probe the ancient star formation history (SFH) independent of local external perturbations, and another is to investigate the possible quenching or delay in star formation near the transition mass below which star formation is strongly suppressed. The most isolated galaxy previously observed with sufficient depth and resolution to distinguish the oldest main-sequence turnoff (MSTO) population is Leo A (DDO 69; Cole et al. 2007), a galaxy on the opposite side of the Local Group from Aquarius, but with similar dynamical mass. Leo A was found to be the dwarf galaxy with the youngest mean age and the smallest fraction of stars older than 10 Gyr of any Local Group member (Hidalgo et al. 2013; Weisz et al. 2014a, and references therein).

In this context, it is of great interest to know if the unusual SFH of Leo A presented in Cole et al. (2007) is an anomaly or a general feature of small, isolated galaxies. In this paper we show that the Aquarius dwarf has many similarities with Leo A, and explore some of the implications of this consistency. We present our observational program and data reductions in Section 2, and the stellar content of Aquarius is shown in Section 3. We derive the SFH of Aquarius in Section 4 and compare its initial low level of star formation and sudden increase some 6–8 Gyr ago to other Local Group dwarf galaxies in Section 5. We put these results into context with recent galaxy formation models in Section 6.

1.2. The Aquarius Dwarf Irregular, DDO 210

The Aquarius dwarf galaxy was discovered on Palomar sky survey prints by van den Bergh (1959) and identified by him as a low surface brightness, very blue galaxy at the limit of resolution into stars. Based on the radial velocity of neutral hydrogen, Fisher & Tully (1975) identified Aquarius as a possible Local Group member. Aquarius is extremely isolated; McConnachie (2012) lists it as one of just three Local Group galaxies with a free-fall time to the Local Group barycenter larger than a Hubble time given its radial velocity.11 Currently only two other galaxies are known to lie within ≈500 kpc of Aquarius; these are SagDIG at a distance of ≈340 kpc, and NGC 6822, ≈440 kpc away. Numerical action models of Local Group dynamics indicate that Aquarius is on its first infall into the Local Group from quite a large distance, and is exceedingly unlikely to have tidally interacted with any other known galaxy (e.g., Shaya & Tully 2013).

Early CCD photometric studies found Aquarius to have an extended history of star formation continuing to very recent times, albeit at a low level (Greggio et al. 1993). Despite a high fraction of cold H i (Young et al. 2003), virtually no current star formation was detected by Hα imaging surveys (van Zee et al. 1997). The consequent lack of bright, evolved stars (e.g., Elias & Frogel 1985; Marconi et al. 1990) confirms Aquarius to be one of the lowest-luminosity gas-rich galaxies in the Local Group (and also resulted in longstanding ambiguity about its true distance). Aquarius is classified as a transition type galaxy, with properties intermediate between the dwarf irregulars and spheroidals (e.g., Mateo 1998). van den Bergh (1994) drew attention to a possible link between its low stellar mass, high gas fraction, and isolated location in the context of the Local Group density–morphology relation (Einasto et al. 1974). The prospective connection to environment was further explored by Grebel et al. (2003) in the light of proposed gas-stripping mechanisms that might produce spheroidals from low-luminosity irregulars.

The basic properties of DDO 210 are summarized in Table 1. The first tip of the red giant branch (TRGB) distance to Aquarius, 950 ± 50 kpc, was published by Lee et al. (1999), who confirmed both the Local Group membership and the unusually low luminosity of the galaxy. This distance was confirmed with Hubble Space Telescope (HST)/WFPC2 photometry by Karachentsev et al. (2002), who derived a distance of 940 ± 40 kpc; later reanalyses found distances of 968 ± 63 (Jacobs et al. 2009) and 977 ± 63 (Tully et al. 2009). McConnachie et al. (2005) found two possible locations for the TRGB and favored the more distant value of 1071 ± 39 kpc, arguing that smaller estimates could be the result of contamination of the brightest part of the red giant branch (RGB) by asymptotic giant branch (AGB) stars. We adopted the more conservative distance in planning the present observations, but in the course of color–magnitude diagram (CMD) modeling the best-fit solutions to the SFH were found to converge on a distance modulus of 24.95 ± 0.10; this corresponds to a distance of 977 ± 45 kpc and is in good agreement with our preliminary analysis of the TRGB magnitude in our data. A more complete reanalysis of the distance to Aquarius will be given in our analysis of the TRGB and variable star populations (E. D. Skillman et al., in preparation).

Table 1. Basic Properties of Aquarius

Parameter Value References
Galactic coordinates (ℓ, b) (deg) 34fdg1, −31fdg4 (1)
Distance (kpc) 977 ± 45 (2)
Line of sight reddening E(BV) (mag) 0.045 (3)
Absolute magnitude MV −10.6 (4)
Half-light radius (arcmin) 1.47 ± 0.04 (4)
Mass within r1/2 (M) 1.8$^{+0.9}_{-0.7}\times 10^{7}$ (5, 2)
Stellar mass (M) 1–2 × 106 (6, 4, 2)
MH i (M) 2.7 × 106 (7)
Current SFR (M yr−1) 1.6 × 10−4 (8)
Stellar metallicity, range (dex) −1.44, ± 0.35 (9)

References. Key to references: (1) Marconi et al. 1990; (2) this paper; (3) Schlafly & Finkbeiner 2011; (4) McConnachie et al. 2006; (5) Kirby et al. 2014; (6) Lee et al. 1999; (7) Begum & Chengalur 2004; (8) Lee et al. 2009; (9) Kirby et al. 2013.

Download table as:  ASCIITypeset image

The first stellar metallicities were obtained by Kirby et al. (2013), who measured spectra in the region of the NIR calcium triplet for 24 red giants, finding a mean metallicity 3.6% of solar (Z/Z = 0.036, or expressed logarithmically and referred to the iron abundance, [Fe/H] = −1.44). This value is consistent with the expectations from the dwarf galaxy mass–metallicity relation of Lee et al. (2006) and typical for a dwarf galaxy of the luminosity of Aquarius. Kirby et al. (2013) found a range of metallicities, with a 1σ spread about the mean of ±0.35 in the log, also typical for dwarfs with a prolonged history of star formation.

Kirby et al. (2014) measured the velocity dispersion of the Aquarius red giants to be σv = 7.9$^{+1.9}_{-1.6}$ km s−1 and from this derived a dynamical mass at the half-light radius of 2.0$^{+1.0}_{-0.8}\times 10^{7}$M and a mass to light ratio ϒ = 32$^{+16}_{-14}$. This result is in good agreement with results based on the H i rotation curve within 550 pc of the center of Aquarius, by Begum & Chengalur (2004). Under the assumption that the neutral hydrogen is distributed in a circularly symmetric disk with inclination iH i = 27°, Begum & Chengalur (2004) found Vcirc = 16 km s−1, leading to an enclosed mass within 570 pc of Mdyn = 3.4 × 107M.

The two dynamical mass estimates agree well within the errors, recognizing that the half-light radius used by Kirby et al. (2014) is 342 pc, about 62% of the distance to the farthest point on the H i rotation curve. Begum & Chengalur (2004) go further, and match the velocity profile to a standard Navarro, Frenk, and White (NFW) dark matter profile (Navarro et al. 1997), finding a best fit maximum velocity at the virial radius V200 ≈ 38 km s−1. This is a major extrapolation, given the virial radius of a few tens of kiloparsecs is many times farther out than the last measured data point. If taken at face value, this profile fit suggests a total mass within the virial radius on the order of 1010M, which allows at least rough comparisons to cosmological simulations.

2. OBSERVATIONS AND DATA REDUCTION

Aquarius was observed by HST under Cycle 20 program GO-12925 for 24 orbits spanning the time period between 2013 June 27 and June 30. Because of the low ecliptic latitude of Aquarius (β = 4fdg9) the observations were made when the angle between Aquarius and the Sun was near the optimum value of ≈135°. The program was split into 12 visits of two orbits each in order to obtain enough separation between exposures to reliably identify short-period variable stars. We followed the Local Cosmology from Isolated Dwarfs (LCIDs) program (Monelli et al. 2010b) in using the F475W (Sloan g') and F814W (Broad I) filters to optimally balance throughput and temperature sensitivity between the oldest MSTO stars and the base of the RGB. We used a subpixel dither pattern between each exposure in order to facilitate cosmic ray rejection, but did not attempt to fill in the gap between the two Advanced Camera for Surveys (ACS)/Wide Field Camera (WFC) chips since the fraction of stars lost in the gap is quite small. The total exposure times in (F475W, F814W) were (22,920, 33,480) seconds.

A three-color composite image built from our co-added images is shown in Figure 1. This tends to highlight the brighter, bluer stars among the foreground stars and many background galaxies. No obvious nebulosities are present, consistent with the lack of current massive star formation in Aquarius. While the distribution of young stars is clumpy and irregular, the red giants are more smoothly distributed, and extend out beyond the boundaries of the ACS/WFC field (McConnachie et al. 2006). We obtained parallel images with the Wide Field Camera 3 (WFC3) and did not convincingly detect any Aquarius members at that radius (≳ 6'), in agreement with the density profile derived by McConnachie et al. (2006).

Figure 1.

Figure 1. Central 3farcm6 × 2farcm3 portion of DDO 210 as seen by ACS/WFC through the F475W (22,920 s) and F814W (33,480 s) filters. The 30'' scale bar corresponds to a physical length of 142 pc at the distance of DDO 210.

Standard image High-resolution image

The charge-transfer efficiency corrected images were processed through the standard HST pipeline and photometry was done using the most recent version of DOLPHOT (Dolphin 2000). Extended objects and residual hot pixels were rejected based on their brightness profiles (|sharpness| < 0.32), and aperture corrections were derived based on relatively isolated stars picked from around the image. Stars that were found to suffer from excessive crowding noise (crowding parameter > 1.0) due to partially resolved bright neighbors were rejected, leaving a sample of 51,239 well-measured stars with signal-to-noise ratio ⩾ 5. We performed an artificial star test analysis with 500,000 stars in order to characterize the measurement errors and incompleteness. Artificial stars with the appropriate point spread function were added to the images with randomly chosen positions and magnitudes between 18 ⩽ m475 ⩽ 31 and −0.75 ⩽ (m475m814) ⩽ 3.50 and the images were re-photometered in identical fashion to the originals. The completeness is above 50% down to limits of (m475, m814) = (28.2, 27.9), reaching the oldest MSTO with high reliability. At these magnitude levels the typical photometric error is ≲0.1 mag, allowing us to resolve the old stellar sequences and make comparisons to isochrone models with a high degree of precision.

The resulting CMD is shown as a Hess diagram in Figure 2. The CMD is dominated by a very narrow RGB that terminates at m814 = 21 and a main sequence indicating stars with a wide range of ages. We planned our observations using the conservative 1071 kpc distance estimate of McConnachie (2012), leading us to expect high completeness only down to the level of the oldest subgiants, but because Aquarius is in fact ≈0.2 mag closer than this, we were able to achieve 50% completeness at the level of the ≈10 Gyr MSTO.

Figure 2.

Figure 2. Our ACS/WFC Hess diagram for Aquarius. Individual stars are plotted where the density is less than 4 stars per decimag2. The colormap is logarithmic in order to show the overall density distribution as well as the structure in the well-populated stellar sequences. Representative error bars are shown, as are the 50% and 90% completeness limits.

Standard image High-resolution image

3. STELLAR POPULATIONS OF THE AQUARIUS DWARF

The stellar populations of Aquarius have been studied several times over the past two decades, with increasing depth. Marconi et al. (1990) made the first CCD-based study of the galaxy and commented on its extremely faint, low surface brightness nature. The most thorough study of the stellar populations of Aquarius up to now has been the study of McConnachie et al. (2006), who combined Subaru and ESO Very Large Telescope (VLT) data to study the red clump/horizontal branch for the first time, discovering that the CMD of the upper giant branch and red clump implied a young overall age for the galaxy, with most of the stars younger than 4$^{+2}_{-1}$ Gyr and fewer than 30% older than 10 Gyr.

Our CMD is the deepest ever presented for Aquarius by more than 3 mag, opening a unique window on the intermediate-age (1–10 Gyr) and old (⩾10 Gyr) stellar populations. Major features in the luminosity function of the RGB are the very well defined red clump centered at m814 = 24.6 and the RGB bump at m814 = 24. The well-defined RGB bump located above the red clump indicates that the dominant stellar population has a mean age quite a bit younger than the Milky Way globular clusters, but not a very high metallicity (Monelli et al. 2010a). There are very few bright AGB stars above the RGB tip, consistent with the very low number of carbon stars reported by Battinelli & Demers (2000) and Gullieuszik et al. (2007), and suggesting that the dominant population is older than 1–2 Gyr.

The sequence of core helium-burning stars begins above the upper main sequence at m814 ≈ 21 and continues downward and to the right through the red clump, reversing to form an extended horizontal branch. The continuity of these sequences shows immediately that while stars of all ages from ancient (blue horizontal branch (BHB)) to ≲108 yr (blue supergiants) are present, Aquarius is dominated by a population tending to the old side of intermediate age (red clump/red horizontal branch; Cole 1999).

Overall, the impression given by the CMD is of a sparsely populated, but in many ways typical, dwarf irregular that has formed stars more or less continuously over a Hubble time. The upper main sequence, brighter than m814 ≈ 24, is sparsely sampled due to the low star formation rate (SFR) and stochastic sampling of the initial mass function (IMF) for stars of the corresponding intermediate and high masses. Based on the low number of main sequence stars brighter than m814 ≈ 22, Aquarius is experiencing a lull in star formation at present despite retaining enough gas to continue forming stars at a significant rate for several more Gyr (e.g., Hunter & Gallagher 1985).

A significant discovery in the context of this project is the BHB, detected as an extension of the red horizontal branch to colors m814–475 ≲ 0.5, overlapping the main sequence at magnitude m814 ≈25.7. These stars are unique signatures of an ancient, metal-poor stellar population, comparable to those seen in the classic satellite dwarf spheroidals of the Milky Way, for example the Sculptor dSph (Salaris et al. 2013), as well as in virtually all dwarf irregulars imaged to sufficient depth (e.g., IC 1613; Cole et al. 1999). By contrast, Leo A was found to have virtually no BHB stars of its own (Cole et al. 2007), despite its similarity to Aquarius in both present-day stellar (McConnachie 2012) and dynamical (Begum & Chengalur 2004; Kirby et al. 2014) masses.

RR Lyrae type variables are identified in the HST observations, and these stars represent a valuable independent probe of the distance and the metallicity of the oldest stellar populations. An analysis of the variable star populations, including the RR Lyraes, will be presented in a forthcoming paper. Because the fraction of BHB stars does not trivially scale to relative SFR at ages ≈10 Gyr or older (e.g., Salaris et al. 2013), the degree to which the intermediate-age population is dominant, and the age and metallicity of the dominant populations, cannot be determined from counts of stars on the BHB alone.

4. STAR FORMATION HISTORY

4.1. Methodology

The SFH presented here has been calculated using the code developed by Cole, which has been used extensively on other Local Group dwarfs and tested against other methods with good consistency (e.g., Skillman et al. 2003; Cole et al. 2007; Monelli et al. 2010b; Cignoni et al. 2012; Lianou & Cole 2013; Skillman et al. 2014). Here, we use the Padova-Trieste (PARSEC; Bressan et al. 2012) isochrones. These are preferred because they have the most up to date input physics, include post-helium-core-flash phases of evolution, cover a wide range of age and metallicity at user-definable grid spacings, and are tuned to the revised solar oxygen abundance (see, e.g., Bahcall et al. 2005). This specific choice plays a role in determining the precise age and duration of episodes of high or low SFR and the overall fit quality of the synthetic CMD, but the broad conclusions are robust against the choice of isochrone library owing to the depth of the photometry (Weisz et al. 2011).

The IMF is that of Chabrier (2003), which is a Salpeter-like power law for M/M ⩾ 0.8 and log-normal for lower masses. Binary stars are included as per the statistics presented in Duquennoy & Mayor (1991) and Mazeh et al. (1992), whereby 35% of the "stars" are single objects, and the rest are binary; in turn, 25% of the binaries are defined as "close" binaries in which the secondary mass is drawn from a flat IMF instead of the Chabrier IMF. The choice of IMF and binary statistics may slightly alter the mass normalization of the solutions but do not significantly affect the shape of the derived SFH, as discussed extensively in Monelli et al. (2010a) and Cignoni et al. (2012), among others. The distance and reddening were initially set to the values derived by McConnachie et al. (2006), but were allowed to vary as demanded by the fitting procedure in order to best match the observed CMD. The stability of the solutions has been confirmed by comparison to fits made independently with different choice of IMF and binary fraction, using the MATCH software package (Dolphin 2002).

The isochrone library contains stars of ages between 6.60 ⩽ log(Age/yr) ⩽ 10.12, spaced by intervals of 0.01, and metallicities from 1/150th to 1/10 solar, evenly spaced at intervals of ≈0.15 in log(Z). Multiple metallicities are allowed at a given age and we do not impose an age–metallicity relation (we do exclude metallicities above one-tenth solar from consideration). The isochrones at each metallicity are grouped into age bins in order to speed the computation; the age bins are initially taken to be evenly spaced by 0.10 in log(age), but the bins may grow or shrink depending on the noise level in the solutions that may justify finer or coarser spacing.

We adopted color and magnitude bins 0.05 and 0.10 mag wide, respectively. Because the data quality is high, we model the region with −0.5 ⩽ (m814m475) ⩽ 2.75 and 19 ⩽ m814 ⩽ 29, but we have checked that our conclusions are not significantly altered if these limits are shifted to brighter magnitudes. The fitting procedure accepted the adopted reddening value of E(BV) = 0.05 and the Tully et al. (2009) distance modulus (m − M)0 = 24.95, but a good fit could not be found using the longer distance modulus from McConnachie (2012). As a check we performed a TRGB fit and found good consistency for distances ≈977 ± 45 kpc; the distance issue will be explored further in our paper presenting the variable star light curves (E. D. Skillman et al., in preparation).

The best-fit synthetic CMD is shown alongside the binned data in Figure 3. The overall fit quality is very high, as shown by the residual significance plot in the right-hand panel. The RGB color and breadth, the magnitude, color, and morphology of the red clump, the upper main sequence color, and the location and slope of the subgiant branches (SGBs) are all well-reproduced. The locations of the blue loop and red supergiant stars are indistinct due to the small number statistics and foreground contamination, but appears to be reasonable for metallicities ZZ/20. The independent fits made using MATCH gave very similar solutions within the uncertainties.

Figure 3.

Figure 3. Observed (left) and best-fit Hess diagrams of Aquarius showing the solution returned by synthetic CMD analysis using the PARSEC isochrone set (Bressan et al. 2012) and software written by Cole. Left: data. Center: model. Right: residual significance, [(data−model)/data]1/2.

Standard image High-resolution image

4.2. The Age of Aquarius

We present the lifetime SFH of Aquarius in Figure 4. The major episode of star formation in Aquarius took place ≈6–8 Gyr ago, although star formation has been continuous (on Gyr timescales) over its entire lifetime. Although star formation at the earliest times is low, it is inconsistent with zero. However, the drop in SFR between ages of 10–12 Gyr is robust and significant; the SFR in this age bin is only nonzero at the ≈2σ level. After the burst of activity at intermediate ages the SFR dropped, and has declined further over the past ≈2 Gyr.

Figure 4.

Figure 4. Maximum likelihood star formation history of Aquarius, derived using the PARSEC isochrone set (Bressan et al. 2012) using software written by Cole. The dominant intermediate-age burst is robust against changes to the model details including choice of isochrone set and distance; see the text for details. The error bars include both random and systematic errors. For clarity the SFH younger than 1 Gyr has been combined into two bins.

Standard image High-resolution image

Aquarius is one of the least luminous gas-rich galaxies in the Local Group, but it does extend beyond the footprint of the ACS/WFC field (202''). McConnachie et al. (2006) trace the structure out to at least a radius of 4', so our field only covers ≈50% of the area over which Aquarius stars are known. This is potentially problematic for interpreting the age of Aquarius based on our SFH, given the ubiquity of age gradients and extended structures in even relatively isolated irregular galaxies—for example, the case of Sextans A and B (Dohm-Palmer et al. 1997; Bellazzini et al. 2014). It has already been noted that the young stars in Aquarius are slightly off-center from the H i contours (McConnachie et al. 2006), which may be suggestive of a possible recent disturbance.

However, there is good reason to believe that our SFH is giving a nearly complete picture of the history of Aquarius. The surface density profiles in McConnachie et al. (2006) drop exponentially, becoming indistinguishable from the foreground population at ≈4' along the major axis; the equivalent distance along the minor axis is 2'. McConnachie et al. (2006) give an exponential fit to the elliptical surface density profile of Aquarius; by integrating their profile over the area of the ACS/WFC footprint, we find that 80% of the Aquarius members brighter than V = 26 are contained in our field. The integral of our SFR(t) over the Hubble time gives a total astrated mass ≈4 × 106M; based on the profiles from McConnachie et al. (2006) we may be missing ≈106M. Even if all of the missed stars were formed prior to 10 Gyr ago, then at the very most the derived full-galaxy SFH would show a time-averaged SFR from 10 to 13.7 Gyr of ≲40% of the measured SFR from 6 to 8 Gyr ago. Even in this unrealistic "oldest case" scenario for the stars in the outskirts of Aquarius, the galaxy as a whole would be dominated by intermediate age stars.

As a further check for possible population gradients, and to attempt to characterize the outer structure of Aquarius, we obtained parallel imaging with the HST WFC3 simultaneously with our ACS images of the central field. The WFC3 field is centered ≈5farcm8 (1.65 kpc) from the center of the ACS/WFC field, nearly along the major axis. Despite a high degree of completeness down to a magnitude of ≈28, we were unable to identify any sign of the stellar sequences in Aquarius in these images. At most a few Aquarius stars per square arcminute could be in the WFC3 field; unless the stellar distribution is highly asymmetric and we have been very unlucky in our parallel field placement, no significant population of Aquarius stars is present at radii ≳ 5' (1.4 kpc), down to very strict limits around the old MSTO.

4.3. Age–Metallicity Relation

The metallicity evolution of Aquarius is summarized in Figure 5, and is in excellent agreement with the spectroscopic metallicities reported by Kirby et al. (2014) and with the photometric estimates in McConnachie et al. (2006). The error bars in Figure 5 show the rms of the mass-weighted metallicity in each age bin; because the isochrone models are separated by discrete intervals of ≈0.15–0.20 in log metallicity, the error bars are typically about this size, except where a very wide range of metallicities are required to reproduce the CMD, as in the oldest age bins. The wide metallicity range for ages ≳8 Gyr likely reflects both the rapid early increase in metallicity and the small number of stars in the CMD at those ages. Because of the rather coarse age resolution at the earliest times, we are unable to unambiguously assign an age to the youngest major populations with Z < 0.01 Z. After the period of rapid early enrichment, ubiquitous in synthetic CMD analyses of dwarf galaxies, the metallicity remains nearly constant over time, slowly rising to a value of ≈ Z/20, fairly typical for dwarf galaxies of this luminosity (e.g., van Zee et al. 1997). The mean metallicity at intermediate ages is in good agreement with the Local Group dwarf galaxy stellar mass–metallicity relationships reported by Berg et al. (2012) and by Kirby et al. (2013).

Figure 5.

Figure 5. Maximum likelihood age–metallicity relation of Aquarius, derived using the PARSEC isochrone set (Bressan et al. 2012). This scale assumes a solar heavy element mass fraction Z = 0.015. A range of metallicities is present at most ages; the error bars show the rms metallicity range of stars at a given age.

Standard image High-resolution image

4.4. Recent SFH

Narrowband imaging and optical spectroscopy at Hα by van Zee et al. (1997) revealed no current massive star formation activity, even though the presence of bright blue main sequence and evolved stars proves that star formation has occurred within the past 108 yr. Aquarius is one of the faintest, lowest SFR galaxies in the 11HUGS sample of gas-rich, star-forming galaxies within 11 Mpc (Kennicutt et al. 2008), with a large discrepancy between the inferred Hα and far-ultraviolet SFR owing to the flickering nature of star formation in dwarf irregulars, combined with stochastic sampling of the IMF at low SFR (Lee et al. 2009; Fumagalli et al. 2011; Weisz et al. 2012a).

A higher resolution view of the SFH over the past 1 Gyr is shown in Figure 6. Aquarius is in a period of low SFR, with seemingly random fluctuations on timescales of ≈108 yr. The recent drop in activity after an upward fluctuation is in excellent agreement with the estimates of current SFR reported in Lee et al. (2009), whose UV (∼100 Myr sensitivity; dashed line) and Hα (∼10 Myr sensitivity; circled star) results are overplotted on Figure 6. The enhancement at ages of 250–300 Myr has been noted by previous authors as discussed above, and is identifiable here by the increased density of blue loop stars at m814 ≈ 22. It is natural to interpret this increase in SFR lasting several tens of Myr as a stochastic fluctuation in SFR of the type frequently observed in dwarf galaxies, leading to the wide variety of Hα, UV, and IR star formation indicators in dwarf galaxies where the gas densities and total mass participating in star formation are small (e.g., Lee et al. 2009). The event is of shorter duration than a typical episode that would be identified as a dwarf galaxy starburst (McQuinn et al. 2010).

Figure 6.

Figure 6. Detailed look at the SFH of Aquarius over the past Gyr. The SFR appears to be dominated by stochastic fluctuations with a very low mean level. The enhancement at ≈250–300 Myr is easily visible in the CMD. The dashed line indicates the SFR derived from GALEX observations, and the circled star gives the Hα SFR (Lee et al. 2009).

Standard image High-resolution image

5. AQUARIUS IN CONTEXT WITH OTHER DWARF GALAXIES

The isolation of Aquarius means that it has never been affected by tidal encounters with the Milky Way or M31, and so it is highly unlikely that it was of much higher mass in the distant past. This removes one complicating factor in the interpretation of the SFHs of close satellites, which may have been of much higher mass in the distant past than they are now (e.g., Peñarrubia et al. 2008). On the contrary, a field dwarf like Aquarius is likely to have only accreted about half of its present-day mass prior to ∼10 Gyr ago (e.g., Brooks & Zolotov 2014).

In making comparisons to other galaxies of the Local Group, it is crucial to remember that systematic errors in the SFH at old ages are unavoidable when shallower CMDs are used (Weisz et al. 2011), and the most meaningful comparisons are to be made comparing CMDs of similar absolute depth (e.g., Weisz et al. 2014a). This tends to limit the gold standard sample of isolated comparison galaxies to the LCID targets (Hidalgo et al. 2013) and Leo T (Weisz et al. 2012b) but still allows some patterns to emerge. Grouped by the overall history of star formation, the galaxy most similar to Aquarius is the late-blooming dIrr Leo A (Cole et al. 2007), which although similarly isolated and of similar dynamical mass was delayed for longer and formed stars much more efficiently after igniting at an intermediate age.

5.1. Leo A (DDO 69)

In Leo A, the next most-isolated galaxy to be studied in equivalent detail, star formation also began at the oldest ages, but continued at only a very low level until ≈6 Gyr ago (Cole et al. 2007). Like Aquarius, Leo A has almost certainly never interacted with another (known) galaxy (e.g., Shaya & Tully 2013). The two galaxies are similar in their degree of isolation, total magnitude MB, color, and gas content, so any similarities that can be identified in their SFH are of great interest.

Comparing the dynamical masses within the half-light radii, Kirby et al. (2014) gives M1/2 = 15$^{+6}_{-5}\times 10^{6}$M for Leo A and 20$^{+10}_{-8}\times 10^{6}$M for Aquarius; the latter value may be reduced by 9% to account for the revised distance estimate presented in this paper (see Table 1). With that correction the two galaxies are seen to have nearly the same total mass, although they have somewhat different mass to light ratios, in part because of the large number of young stars in Leo A. Leo A also supports a higher current SFR, has a higher gas fraction (McConnachie 2012), and a slightly lower present-day metallicity (Kirby et al. 2013).

In both galaxies, the total stellar mass reported by different authors varies based on the adopted distance and IMF, and details of the integrated light photometry. Most reported stellar mass estimates for Leo A ultimately derive either from the work of Lee et al. (2003) or Vansevičius et al. (2004); correcting these to the consistent distances derived in Dolphin et al. (2002) and Cole et al. (2007) we conclude that the best estimate of the stellar mass of Leo A is ≈3 ± 1 ×106M. Fewer measurements are available for Aquarius; using the method of Lee et al. (2003) to convert from integrated magnitude and B − V color to stellar mass indicates that Aquarius has ≈1/2 (using the photometry of Lee et al. 1999) to ≈1/3 (using the photometry of McConnachie et al. 2006) of the stellar mass of Leo A, ≈1–2 × 106M.

The CMDs of Aquarius and Leo A, corrected for distance and reddening, are shown side by side in Figure 7, highlighting their overall similarity and differences in detail. Apart from the quite different HB/red clump morphology, the biggest difference is in the SGBs, where Leo A is much more strongly concentrated around M8140 ≈ +1.5–2, with only a scattering of fainter stars.

Figure 7.

Figure 7. Our CMD of Aquarius, corrected for distance and reddening, beside the analogous CMD of Leo A (Cole et al. 2007). Differences in the red clump/HB morphology and the extent of the subgiant branch make plain that Leo A and Aquarius do not share a common quenching time. While both galaxies contain genuinely ancient stars, the relative fraction of old stars is over twice as large in Aquarius.

Standard image High-resolution image

In order to quantify the similarities and differences in the SFH of Aquarius and Leo A, we rederive the Leo A SFH using the data from Cole et al. (2007) with the identical stellar evolution library and time bins used in our solution for Aquarius. The results are entirely consistent with the SFH reported in Cole et al. (2007); the two galaxies are directly compared in Figure 8. It is seen that the drought in star formation was broken earlier in Aquarius, but that the galaxy then rapidly dropped back down to nearly its lifetime average SFR. By contrast, Leo A formed stars at a very slow rate over much of its lifetime and only managed to start converting gas into stars in earnest ≈2 Gyr after Aquarius. The median age of star formation (the time by which 50% of the mass of all stars ever formed had been created) for each galaxy is shown with the vertical dotted lines; the median age of Aquarius is 6.8 Gyr, but for Leo A it is 4.2 Gyr.

Figure 8.

Figure 8. Cumulative fractional star formation of Aquarius (blue) and Leo A (orange) compared to each other and to other Local Group galaxies. The redshift scale at the top is calculated assuming a concordance ΛCDM cosmology. Left: the median ages of the two galaxies are indicated by vertical dotted lines. The dashed line shows a constant SFR, while the convex-upward solid curve shows a reference dark matter assembly history for a 1010M simulated galaxy (Fakhouri et al. 2010). The shaded portion of the graph shows the era of cosmological reionization, from z ≈ 6–14. Right: the cumulative SFH of Aquarius compared to seven other galaxies with comparably precise SFH data (see the text for details). The estimated stellar mass (log M/M; McConnachie 2012) is given for each galaxy. Aquarius and Leo A are the only two galaxies in the sample of eight to show SFRs significantly below their lifetime average for the first few Gyr after reionization.

Standard image High-resolution image

The absolute SFRs, the integral of SFR over lifetime, and metallicity evolution of the two galaxies are compared in Table 2. The two galaxies are clearly of a very similar nature, with a nearly equal total astrated mass. It is interesting that, with the two solutions run completely independently and with no enforcement of a chemical evolution law, Aquarius is found to be more metal-rich than Leo A at almost every age. For the intermediate age stars that dominate the red giant sample of Kirby et al. (2013), the metallicity difference between the two galaxies is of the same magnitude and in the same sense as observed, with Leo A about 0.1–0.2 dex more metal-poor. While Leo A is ≈1 mag brighter than Aquarius at the present day, our modeling finds that its total lifetime astrated mass has been less than Aquarius; this difference is also in the same sense, and of the same order of magnitude as the dynamical mass difference reported by Kirby et al. (2014). Presumably the difference in metallicity may be related back to the more gas-rich nature of Leo A, which is a factor of two to three more gas-rich than Aquarius (McConnachie 2012) and therefore less chemically evolved. At the oldest ages, Aquarius did not just form stars at a higher relative rate than Leo A, but had a higher absolute SFR (a fact that is reflected in its much more prominent BHB).

Table 2. Comparative SFH of Aquarius and Leo A

Lookback Time Redshifta Aquarius Leo A
(Gya) (z) SFRb (10−4 M yr−1) Mc (106M) Z/Zd SFR (10−4 M yr−1) M (106M) Z/Z
0–0.5 0–0.04 0.76 $^{+0.06}_{-0.06}$ 4.03 0.065 1.82 $^{+0.10}_{-0.10}$ 3.54 0.051
0.5–1 0.04–0.08 1.04 $^{+0.23}_{-0.23}$ 3.99 0.065 1.73 $^{+0.14}_{-0.15}$ 3.45 0.051
1–2 0.08–0.2 1.26 $^{+0.44}_{-0.37}$ 3.94 0.056 3.19 $^{+0.23}_{-0.20}$ 3.36 0.041
2–4 0.2–0.4 2.77 $^{+0.27}_{-0.44}$ 3.81 0.053 5.70 $^{+0.53}_{-0.55}$ 3.05 0.039
4–6 0.4–0.7 3.28 $^{+0.88}_{-0.27}$ 3.26 0.052 6.31 $^{+0.51}_{-0.80}$ 1.91 0.039
6–8 0.7–1.1 8.26 $^{+0.39}_{-1.87}$ 2.60 0.034 1.33 $^{+0.14}_{-0.69}$ 0.64 0.022
8–10 1.1–1.8 2.38 $^{+1.51}_{-0.26}$ 0.95 0.032 0.99 $^{+0.30}_{-0.10}$ 0.38 0.026
10–11.5 1.8–3 0.50 $^{+0.16}_{-0.27}$ 0.48 0.020 0.12 $^{+0.26}_{-0.10}$ 0.18 0.013
⩾11.5 ⩾3 2.22 $^{+0.96}_{-0.97}$ 0.40 0.007 0.90 $^{+0.22}_{-0.64}$ 0.16 0.007

Notes. aAssuming Ω = 1, ΩΛ = 0.73, H0 = 71 km s−1 Mpc−1; Wright (2006). bAverage star formation rate over the given time period. c∫SFR(t)dt from formation until the end of the given time period (excludes losses due to finite stellar lifetimes). dAverage metallicity (mass fraction of elements heavier than helium) of stars formed during the given time period, compared to the solar metallicity Z = 0.0152.

Download table as:  ASCIITypeset image

5.2. Other Local Group Galaxies

All other galaxies with comparably deep CMDs (in absolute magnitude) are much closer to the Milky Way12 and therefore complicated by the possibility that past interactions have strongly influenced their SFH. McConnachie (2012) shows that the irregular and transition type galaxies IC 1613, Leo T, Phoenix, and LGS-3 are all at least a factor of two closer to a massive spiral than Aquarius. The isolated systems WLM and Sagittarius DIG have been surveyed with WFPC2 (see Weisz et al. 2014a), but not to a depth sufficient to resolve the oldest MSTO; the deepest ACS data for UGC 4879 reach the level of the red clump (Jacobs et al. 2011).

The non-satellite dwarf spheroidals Tucana and Cetus are somewhat different. They differ from the gas-rich galaxies in their SFHs; this may or may not be related to a proposed history of tidal interactions (Fraternali et al. 2009; Lewis et al. 2007, for Tucana and Cetus respectively). Nevertheless while interactions may be suspected for Tucana and Cetus (e.g., Teyssier et al. 2012), the timing and identity of the main tidal perturber cannot be established with confidence. They are included here for completeness.

Most of the Local Group dwarf galaxies show either nearly constant long-term average SFR, or are dominated by much older bursts. The latter pattern is more common in the nearby satellites of the Milky Way and M31 and among the lowest-luminosity systems. The Tucana and Cetus dwarf spheroidals are nearly pure ancient populations (albeit with ages differing by ≈1 Gyr; Monelli et al. 2010b, 2010a). The dwarf irregular galaxies, which statistically tend to be found at large distances from both M31 and the Milky Way, show more uniform SFRs over cosmic time. The prototypical example with comparable ACS data could be IC 1613, a far larger galaxy (MV = −15.2) which shows nearly constant SFR (Skillman et al. 2014). At the other end of the mass scale, Leo T (MV = −8) is observed to show a slowly declining SFR beginning at the earliest times and continuing until ≈25 Myr ago without dramatic variation (Weisz et al. 2012b; Clementini et al. 2012).

From Figure 8, we see that there is a wide range of observed SFHs at a given galaxy mass. These galaxies are not satellites of either the Milky Way or M31, so no tidal or ram pressure quenching of star formation is expected, and indeed most of the systems show extended star forming lifetimes. Among the galaxies with a freefall time onto either the Milky Way or M31 less than a Hubble time, no examples of significantly delayed peak SFR are found. Leo A and Aquarius, the two most isolated galaxies in the sample, on the other hand, simmer for at least 3–4 Gyr after reionization before dramatically increasing their SFR. Present-day isolation cannot be the sole factor determining the early SFH, however, because the similarly distant and similarly massive dwarf spheroidals, Cetus and Tucana, have the highest fraction of stars older than 10 Gyr among the galaxies in our sample.

The comparison to other Local Group galaxies serves to reinforce the point that there is no simple relationship between current galaxy mass, isolation, and age of the main stellar population; it is physically intuitive to suppose that such a relationship would be clarified if the complete history of interaction and feedback heating for each galaxy could be known. If there are any trends, it seems that major star formation events are more evenly distributed in redshift than in lookback time, and galaxies with low baryonic mass are more prone to bursty behavior than more massive and/or gas-rich systems.

6. IMPLICATIONS FOR GALAXY FORMATION MODELS

A significant amount of recent progress has been made in cosmological structure formation models describing the formation and early evolution of dwarf galaxies. The gains in numerical resolution and treatment of various radiative and hydrodynamical effects are reflected in simulated galaxies that plausibly resemble observed dwarfs.

In the most thoroughly studied galaxies, their proximity means that environmental effects arising from close passages to the Milky Way must be considered when interpreting the SFH. This has weakened inferences about the effects of the UV background or internal feedback on the SFH of dwarfs. While the hypothesized effects are expected to be strongest for the lowest-mass galaxies, most of the extremely low-luminosity ("ultrafaint") dwarfs are deep within the Milky Way potential well and therefore present an undesirable degree of complexity despite being the obvious place to test for feedback or photoionization heating. Aquarius and Leo A have certainly never had a close flyby of either M31 or the Milky Way, and have likely never been within the virial radius of either galaxy (Shaya & Tully 2013); they are the two best "clean" galaxies we have for comparison to model predictions.

Fakhouri et al. (2010) trace the assembly of dark matter halos beginning with 106M clumps at redshift z = 15 and find that the merger rate is approximately constant per unit redshift, which means that typical dwarfs assembled most of their mass at very high lookback time. Their examination of the Millennium simulation finds that the median redshift of the last major merger for 1010M halos is z = 3.4 (tlookback = 11.8 Gyr). The SFH derived assuming that SFR follows the accretion rate for a typical halo from Fakhouri et al. (2010) is shown as the smooth curve in Figure 8.

While the median dwarf is expected to form early, Fakhouri et al. (2010) also find a broad distribution of accretion rates at all redshifts for small halos. Their models predict that for halos with present-day virial mass 1010M, roughly 1 in 3 experienced their last major merger after z = 2, and about 1 in 6 since z = 1. The implication is that while late accretion is not the most likely evolutionary path for a dwarf galaxy the size of Aquarius, it is far from unlikely. The fact that both of the first two isolated dwarfs to have been studied in sufficient detail show just such delays leads to the impression that either long delays are more common in isolated galaxies than in the galaxy population as a whole, or that star formation does not track dark matter accretion. If late mergers are the explanation for the delayed SFH, then neither component of a major merger could have had a dominant ancient stellar population.

While a simple prescription for baryon (re)heating after reionization might predict the complete quenching of star formation in dwarf galaxies at redshifts z ≈ 6, there is little evidence in the SFHs of observed galaxies for such a quenching (Weisz et al. 2014b). Rather, the effects of reionization, if any, appear to be more subtle; note in Figure 8 that no isolated Local Group dwarf begins vigorously forming stars within 1–2 Gyr after the end of reionization. While star formation does not appear to be quenched, an argument could be made that it is at least dampened by the UV photoheating: every galaxy in Figure 8 has a lower rate of star formation immediately after the end of reionization than before.

Cast in the framework of hierarchical galaxy assembly in a ΛCDM cosmology, it appears likely that the accretion history of baryons in Aquarius either differs significantly from that of the dark matter, or baryonic physics has disrupted the one-to-one link between mass accretion and star formation. The observed complexity of these relationships has proven challenging for non-fully cosmological models to capture, and predictions for stellar and dark matter mass growth are not always in agreement with observations (e.g., Leitner 2012; Behroozi et al. 2013; Lu et al. 2014a, 2014b; Weisz et al. 2014b).

The most recent generation of cosmological structure formation models features high numerical resolution and advanced prescriptions for feedback, star formation, variable UV background, and radiative cooling through both H2 and metal lines (e.g., Christensen et al. 2012; Shen et al. 2014; Hopkins et al. 2013). The effects of reionization and feedback are typically more subtle and varied in recent simulations than in previous work, and it appears that progress is being made on the overcooling front.

Figure 9 compares the SFH of Aquarius and Leo A to the predictions from three recent sets of simulations that highlight the problems associated with dwarf galaxy formation and the recent progress that has been made. The left panel shows six isolated model dwarfs from the hydrodynamical simulations of Sawala et al. (2011). The simulated galaxies have a dark matter halo mass of ≈1010M, comparable to or a bit larger than Aquarius. In the simulations, star formation is suppressed by the combined effects of feedback and reionization, but unlike Aquarius and Leo A, the SFR never recovers. Despite an advanced cooling and star formation treatment, the simulations form the majority of their stars before z ≈ 2 because SFR tracks the accretion rate. The simulated galaxies also form a minimum of an order of magnitude too much stellar mass compared to the observed dwarfs (5–10×107M; for comparison see Table 2).

Figure 9.

Figure 9. Cumulative fractional star formation of Aquarius (blue) and Leo A (orange), compared to cosmological simulations including baryons. Left: simulated dwarfs from Sawala et al. (2011) for galaxies with similar halo mass to Aquarius and Leo A; the SFH tracks the dark matter assembly and forms too many stars too early. Right: simulated dwarfs from Shen et al. (2014; solid lines, color-coded following their Figure 5) and from Brooks & Zolotov (2014; dashed lines). In both panels, the dark matter-only simulated 1010M halo of Fakhouri et al. (2010) is shown as the smooth curve.

Standard image High-resolution image

Even more recent models are shown in the right hand panel. The dashed lines show the three galaxies in the simulations of Brooks & Zolotov (2014) that have produced on the order of 106–7M worth of stars, similar to Aquarius. These simulated isolated dwarfs are gas-rich (MH i/M > 1), dark matter-dominated, and have experienced a bursty SFH. The most gas-rich, lowest stellar-mass dwarf is an obvious comparison object to Leo A, with some initial star formation prior to the end of reionization, followed by a long fallow period and a late burst. The other two simulated dwarfs also show an initial epoch of star formation with a post-reionization drought, although star formation restarts after a much briefer delay. Brooks & Zolotov (2014) conclude that some isolated dwarfs are capable of retaining gas without forming stars for several Gyr. A complete analysis of the relationship between mass accretion, feedback, gas cooling, and SFR in these simulations will be presented in a forthcoming paper (Weisz et al., in preparation), but the general agreement is quite promising and may indicate that a critical threshold in resolution has been crossed.

The solid lines in the right-hand panel of Figure 9 show the results of Shen et al. (2014), whose high-resolution smooth-particle hydrodynamical model also shows that the combination of reionization and the mechanical energy input by supernova feedback is critical to the prevention of gas cooling in small halos. They make a strong prediction that halos with total mass ≲109M never manage to cool enough to form stars; this threshold is just below the expected virial masses of Aquarius and Leo A. Four of their simulated galaxies span the range of plausible masses for our isolated dwarfs. The nearly constant SFR dwarfs "Bashful" and "Doc" inhabit more massive halos of ≈1–3×1010M. In Figure 9, these galaxies have formed the bulk of their stars far earlier than Aquarius or Leo A, and have also formed over an order of magnitude more stellar mass that the isolated Local Group dwarfs.

The other two star-forming galaxies from Shen et al. (2014), "Dopey" and "Grumpy," are at the lower end of the applicable halo mass range, and present a stark contrast to the more massive systems: they are extremely gas-rich and have formed one to two orders of magnitude fewer stars than Aquarius or Leo A. The low-mass simulated dwarfs exhibit a highly bursty SFH, and are characterized by an extreme delay of star formation—neither galaxy experienced any star formation until very recently, making them exaggerated versions of Leo A. The variance in SFH between any two given galaxies of similar total mass indicates a large degree of stochasticity in all three models considered here, a phenomenon which is also expected in the models of, for example, Starkenburg et al. (2013) and others.

The fact that realistic small galaxies form at all in the simulations represents an encouraging sign of progress, but there are still significant discrepancies. A lingering sign of the overcooling problem is that the extrapolated Begum & Chengalur (2004) NFW profile mass for Aquarius is ≈1010, similar to Shen's "Bashful" and "Doc," but these produce ten times as much stellar mass as Aquarius and still have a very large amount of gas left over. Scaling the mass down by 0.5–1 order of magnitude for the other two simulated galaxies reduces the number of stars formed by two to three orders of magnitude, indicating that we are in a mass range where the feedback and star formation efficiency prescriptions are exceptionally sensitive to halo mass. By contrast, the most and least massive galaxies in Figure 7 (IC 1613 and Leo T, respectively) share broadly similar patterns of relative SFR, are both gas-rich at the present time, and are both older than the intermediate-mass galaxies Aquarius and Leo A.

7. SUMMARY AND CONCLUSIONS

We have measured the SFH of the Aquarius dwarf galaxy based on photometry extending to the oldest MSTO obtained with the HST ACS. We have fit the photometric data using the method of creating synthetic CMDs based on theoretical isochrones and a maximum-likelihood comparison. We have found that Aquarius is dominated by stars aged less than 10 Gyr, and that metallicity evolution over the past ≈10 Gyr has been mild. We find in particular the following major features.

  • 1.  
    The SFH is not consistent with a strong, early burst of star formation prior to the end of reionization. However, a significant population of stars older than 10 Gyr is present based on the SFH, the magnitude extent of the subgiant branch, and the detection of a BHB.
  • 2.  
    The first major episode of star formation began ≈8 Gyr ago and lasted for ≈2 Gyr. Star formation has lasted for the full age of the universe, with no evidence for long breaks between active periods. Apart from the increase at 8 Gyr, the SFH is almost consistent with a constant lifetime SFR. In this sense Aquarius resembles Leo A, which is similarly isolated and of comparable mass, but experienced an even stronger delay. The presence of an old stellar population amounting to 4×105M in place prior to redshift z = 3 was evidently not accompanied by sufficient feedback to quench subsequent star formation.

Aquarius is the most isolated galaxy for which the full SFH has been measured with high precision and accuracy based on photometry reaching the oldest MSTO; it therefore makes an ideal test case to benchmark the evolution of small galaxies in isolation. Either Aquarius experienced a "slow drip" accretion history, in which gas trickled in over several gigayears before reaching a threshold density for star formation, or feedback effects on the gas have decoupled the SFR from the gas accretion rate.

The differences between Aquarius and Leo A suggest that some combination of inefficient cooling in small halos, heating due to feedback, and/or a stochastic history of baryon accretion is necessary to understand the evolution of the smallest galaxies. Because of the high level of variance observed in the SFH of isolated dwarf galaxies to date, observations of additional galaxies to comparable photometric depth will be necessary in order to begin to map correlations in the space of total mass, baryon fraction, angular momentum, and environment.

We thank Sijing Shen and Till Sawala for sharing the results of their simulations with us. Support for D.R.W. is provided by NASA through Hubble Fellowship grant HST-HF-51331.01 awarded by the Space Telescope Science Institute. A.B., E.D.S., and D.R.W. were supported in part by National Science Foundation grant No. PHYS-1066293 and the hospitality of the Aspen Center for Physics. R.L. acknowledges financial support from the DAGAL network from the People Programme (Marie Curie Actions) of the European Union's Seventh Framework Programme FP7/2007–2013/under REA grant agreement number PITN-GA-2011-289313. Additional support for this work was provided by NASA through grant number HST GO-12925 from the Space Telescope Science Institute, which is operated by AURA, Inc., under NASA contract NAS5-26555. This research made extensive use of NASA's Astrophysics Data System bibliographic services.

Footnotes

  • Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. These observations were obtained under program GO-12925.

  • 11 

    The others being the Sagittarius dwarf irregular (SagDIG, ESO 594-4) and VV124 (UGC 4879).

  • 12 

    Or in a few cases, M31.

Please wait… references are loading.
10.1088/0004-637X/795/1/54