the present paper is part of the local cosmology from isolated dwarfs ( lcid ) project .
we have obtained deep _ hubble space telescope _ ( hst ) photometry , reaching the oldest main sequence turn - off points , of six isolated dwarf galaxies in the local group : ic  1613 , leo  a , cetus , tucana , , and phoenix .
five galaxies were observed with the advanced camera for surveys ( acs , * ? ? ? * ) , while phoenix was observed with the wide field and planetary camera-2 ( wfpc2 , * ? ? ? * ) .
the main goal of the lcid project is to derive the star formation histories ( sfhs ) , age - metallicity relations ( amrs ) , variable star populations , and stellar population gradients of this sample of galaxies .
our objective is to study their evolution at early epochs and to probe effects of cosmological processes , such as the cosmic uv background subsequent to the onset of star formation in the universe or physical processes such as the gas removal by supernovae ( sne ) feedback .
our sample consists of field dwarfs which were chosen in an effort to study systems as free as possible from environmental effects due to strong interactions with a host , massive galaxy .
the sfh is a powerful tool to derive fundamental properties of dwarf galaxies and their evolution @xcite , but to study the earliest epochs of star formation , deep cmds , reaching the oldest main sequence turn - offs , are required ( cf .
our impressions of the evolution of dwarf galaxies are biased by studies of the nearby , gas - poor dsph companions of the milky way .
because of their proximity , deep cmds are obtainable from ground - based observatories .
in contrast , the gas - rich , dirr galaxies , which are found at greater distances , have few studies with comparable photometric depth .
the magellanic clouds represent the one exception to this generalization , but , because of
they are substantially more massive than the typical dwarf galaxy and have potentially complex histories due to their current interactions ( see * ? ? ?
* and references therein ) , they present less than ideal targets for study .    because of the larger distances to the dirrs , until the lcid project , none have had resolved star studies which reach down to the oldest main sequence turnoff stars .
thus , our view of the sfhs of dwarf irregular galaxies has been shaped by relatively indirect measures .
the work by @xcite represents a seminal contribution to our understanding of the sfhs of dirrs .
based on galaxy mass estimates , blue luminosities , and h@xmath1 luminosities , they demonstrated that most irregular galaxies were consistent with nearly constant star formation over their lifetimes ( as opposed to the larger spiral galaxies which showed declining star formation rates , hereafter sfrs ) .
most importantly , they pointed out that `` the constant sfr history implies that the simple classical model in which star formation is proportional to gas density in a closed system can not be correct for irregular galaxies . ''
modern observations of nearby dirrs have been collected and summarized by @xcite , and indicate that , on average , sfrs are higher at early times , but that there are some dirrs for which constant sfrs are a good approximation .
recently , there have been theoretical papers emphasizing the difficulty in making models of galaxies which have nearly constant sfrs .
@xcite , @xcite , @xcite , and @xcite have all highlighted the difficulty of producing dwarf galaxies with properties comparable to those observed in the current universe .
the degree of failure is greatest in the amount of mass converted into stars , which is of order one magnitude too large in cosmological simulations .
this overproduction of stars in dwarfs is an extreme symptom of the `` over - cooling problem '' ( c.f .
* ) faced by all galaxy modeling .    without deep imaging of resolved stars in dirr galaxies ,
we lack sufficient time resolution to study sfhs at the earliest times .
for example , the earliest time bin used by @xcite is 4 gyr in duration . at this time resolution
, we can not distinguish between sfrs that are constant at all times , or sfrs that show considerable variation ( most importantly evidence for an early dominant episode of sf ) .
the initial time bin of 4 gyr covers the time range from before re - ionization up to a redshift of @xmath22 ( which is approaching the peak of in the history of universal star formation , e.g. , * ? ? ?
* ) , and represents this entire important range with a single average number .
the observations of ic  1613 presented here resolve this initial period and represent a small step into this relatively unexplored territory .      here
we present our analysis of a deep hst / acs observation of a field in ic  1613 . as we discuss in section [ repfield ] ,
the sfh that we derive for this field is likely a good representation for the entire galaxy .
ic  1613 is a low - luminosity , local group , dirr galaxy which is one of the nearest gas - rich irregular galaxies ( for a review of the properties of ic  1613 see * ? ? ?
because of its proximity , ic  1613 offers the opportunity to reconstruct a detailed sfh of a relatively isolated and non - interacting dwarf irregular galaxy .
ic  1613 also has very low foreground and internal reddening ( although it lies within 5 degrees of the ecliptic ) .
there have been several determinations of its distance based on cepheid and rr  lyr variable stars and the tip of the red giant branch .
using all three methods , @xcite derived a distance of 730 kpc . in a comparison of all of the literature values ( which included their own measurement using rr lyrae and cepheid data of 770 kpc ) , @xcite
determine a mean distance of 760 kpc .
the most recent determination is by @xcite who derived a distance of 740 kpc using cepheids . for consistency with the other lcid studies
, we will adopt the distance estimate from @xcite which corresponds to scales of 221 pc arcminute@xmath3 , 3.7 pc arcsecond@xmath3 , and @xmath2 0.2 pc pixel@xmath3 .
ic  1613 is not considered to be a satellite of either the milky way or m31 .
@xcite included ic  1613 in the diffuse `` local group cloud , '' and @xcite determined a distance of 517 kpc and a velocity of @xmath490 km s@xmath3 relative to the local group barycenter .
as such , it is located very close to the zero velocity surface for m  31 and well within the zero velocity surface for the local group @xcite . given its isolated position and velocity , ic  1613
has not had any recent interactions , although interactions with other galaxies long ago can not be ruled out . to date , there are no proper motion studies of ic  1613 , which would be very valuable in determining its potential interaction history .
the physical parameters of ic  1613 were summarized in @xcite .
these properties are normal for an i m  v galaxy with a moderate luminosity ( @xmath515.2 ) and a small value of the maximum amplitude of the rotation curve ( v@xmath6 @xmath7 25 km s@xmath3 ; * ? ? ?
its sfr of 0.003  @xmath8  yr@xmath3 @xcite is also normal for its type and luminosity .
the most recent ism oxygen abundance measurement was made by @xcite , who obtained a spectrum where the @xmath94363 auroral line of [ ] was detected resulting in a measurement of 12 @xmath10 log ( o / h ) @xmath7 7.62 @xmath11 0.05 .
this corresponds to 8.5% of the solar oxygen abundance ( as determined by * ? ? ?
* ) , which is slightly less than that of the smc , and normal for a galaxy of its luminosity ( e.g. , * ? ? ?
* ; * ? ? ?
this combination of proximity and normality makes ic  1613 one of the best opportunities to study the properties of a dwarf star - forming galaxy that is relatively isolated ( as is typical for i m  v galaxies ) .
the overall structure of ic  1613 is also typical of an i m  v galaxy .
ground - based studies have provided an overview of the stellar population distributions .
@xcite studied the distribution of luminous cool stars from j and k - band imaging and found agb stars covering a wide range in age in all of their inner galaxy fields . @xcite conducted a wide field survey of ic  1613 for c and m stars , and found the ( intermediate age ) c stars extended out to 15 arcminutes , well beyond the regions where star formation currently is active .
recently , @xcite have conducted a wide field optical survey of ic  1613 , and trace red giant branch ( rgb ) stars out to radii @xmath12 16.5 arcminutes ( @xmath2 3.6 kpc ) , showing the galaxy to be more extended than previously thought .
the resolved stellar populations of ic  1613 have been studied with the hst twice in the past , both times using the wfpc2 camera .
@xcite studied a central field and found ic  1613 to be a smoothly evolving galaxy with a relatively constant sfr over the last gyr .
horizontal branch ( hb ) stars were detected , indicating the presence of an old population .
@xcite obtained deep imaging for a field located 7.4 arcminutes southwest of the center . while that imaging was not quite deep enough to reach to the oldest main sequence turnoff stars , greatly limiting the time resolution at the oldest ages , the derived sfrs were constant within a factor of three over the entire lifetime of the galaxy .
in this paper , we present the sfh of ic  1613 obtained from observations with the acs on the hst .
the photometry reaches the oldest main sequence turn - offs of the galaxy , allowing us to obtain an accurate sfh even for the oldest stellar populations .
@xcite have already used these observations to conduct a study of the variable star content of ic  1613 .
the structure of the paper is as follows : in  [ secred ] the observations and data reduction are discussed and the cmd is presented .
the derived sfh of ic  1613 is presented in
[ secsfhresult ] and is compared with those of other lcid galaxies in  [ seclcid ] .
the implications of the sfh of ic  1613 for galaxy modeling , and , in particular , the over - cooling problem are discussed in  [ seccosmo ] .
the main conclusions of the work are summarized in  [ seccon ] . as with the previous lcid papers ,
cosmological parameters of @xmath13 , @xmath14 , and a flat universe with @xmath15 are assumed ( i.e. , * ? ? ?
* ) .    ) at 45 ( @xmath2 1.0 kpc ) and the half - light radius ( r@xmath16 ) at 65 ( @xmath2 1.4 kpc )
as can be seen from the figure , the hst acs field is located between the two .
also indicated are the positions of the two previous hst wfpc2 fields ( chevrons ) from @xcite ( inner ) and @xcite ( outer ) . [ f1],width=325 ]
the acs observations of ic  1613 were obtained between august 18 and 20 , 2006 .
the f475w and f814w bands were selected as the most efficient combination to trace age differences at old ages , since they provide the smallest relative error in age and metallicity in the main - sequence and sub - giant regions ( see c.  gallart et al .
2013 , in prep ) .
total integration times were 31,489 s in f475w and 27,119 s in f814w .
the observations were organized into 12 visits of 2 orbits each , and each orbit was split into one f475w and one f814w exposure ( in order to maximize sampling of variable star light curves ) .
the observing log is reported in @xcite .
dithers of a few pixels between exposures were introduced to minimize the impact of pixel - to - pixel sensitivity variations ( `` hot pixels '' ) in the ccds .
the observed field of ic  1613 is shown in figure [ f1 ] . at the adopted distance to ic  1613
, the footprint of the acs covers 0.56 kpc@xmath17 .
the optical scale length of ic  1613 is 29 @xcite and the stellar distribution can be traced out beyond 15 arcminutes , so the 34 @xmath18 34 format of the acs covers only a fraction of the area of ic  1613 ( @xmath2 9% ) .
we analyzed the images taken directly from the stsci pipeline ( bias , flat - field , and image distortion corrected ) .
two psf - fitting photometry packages , daophot / allframe @xcite and dolphot @xcite , were used independently to obtain the photometry of the resolved stars .
non - stellar objects and stars with discrepant and large uncertainties were rejected based on estimations of profile sharpness and goodness of fit .
see @xcite for more details about the photometry reduction procedures .
individual photometry catalogs were calibrated using the equations provided by @xcite .
the zero - point differences between the two sets of photometry are small ( @xmath19 0.04 mag ) and typical for obtaining hst photometry with different methods @xcite .
we direct the reader to extensive photometry reduction comparisons of lcid observations discussed in @xcite and @xcite . for simplicity ,
the rest of this paper is based on only the dolphot photometry dataset which contains 165,572 stars .
signal - to - noise limitations , detector defects , and stellar crowding can all impact the quality of the photometry of resolved stars with the resulting loss of stars , changes in measured stellar colors and magnitudes , and systematic uncertainties . to characterize these observational effects
, we injected @xmath20 artificial stars in the observed images and obtained their photometry in an identical manner as for the real stars . @xcite and @xcite
provide detailed descriptions of the procedures we adopt for the characterization and simulation of these observational effects .
, evenly spaced by factors of 2 .
the left axis shows observed magnitudes corrected for extinction .
the right axis shows absolute magnitudes .
a distance modulus of @xmath21 and an extinction @xmath22 , @xmath23 have been used .
the lines across the bottom of the right panel show the 0.25 , 0.50 , 0.75 and 0.90 completeness levels .
three isochrones from the basti stellar evolution library have been over plotted for comparison .
[ f2],width=325 ]    the cmd of ic  1613 is shown in figure [ f2 ] .
individual stars are plotted in the left panel and density levels are shown in the right panel .
the left axis shows magnitudes in the acs photometric system corrected for extinction .
absolute magnitudes are given on the right axis using the adopted values for the distance modulus ( @xmath24 ) and extinctions ( @xmath25 and @xmath26 mag ) from @xcite .
the completeness factor , @xmath27 , the rate of recovered artificial stars as a function of magnitude and color , is shown in figure [ f2 ] , right panel , at @xmath28 and 0.90 . finally , in order to highlight the main features of the cmd , three isochrones from the basti stellar evolution library @xcite are also shown in the right panel as labeled in the inset .    as shown in figure [ f2 ]
, our photometry reaches below the oldest main - sequence turn - off at the 50% completeness limit , allowing for very strong constraints on the oldest epochs of star formation .
these observations are @xmath29 mag fainter than the deepest cmd previously obtained for this galaxy @xcite .
since there are no early isolated episodes of elevated sfr in the sfh , there is no obvious indication of the time resolution in the early sfh of ic  1613 shown in figure [ f2 ] .
however , extensive modeling with similar lcid datasets has shown that observations of this quality yield a time resolution of @xmath2 1 gyr @xcite .
consistent with earlier observations , a comparison of the observations with the over - plotted isochrones and the presence of an extended horizontal branch indicate that a very old , very low metallicity stellar population is present in the galaxy . the gap produced in the hb by the rr - lyrae variables and a red clump ( rc ) at the red end of the hb
are clearly defined .
an ms with stars younger than 200 myr is clearly apparent .
some blue and red core helium - burning stars might also be present above the rc ( @xmath30 ; @xmath31 ) .
the rgb bump can be also observed at @xmath32 and @xmath33 mag ; this broad feature can be explained by an extended sfh @xcite .
as can be seen from figure [ f1 ] , the hst acs field of view covers only a small fraction of ic  1613 ( @xmath2 9% ) .
the position of the field was chosen mainly as a balance between optimizing the number of stars but minimizing the effects of crowding at the photometric limit .
however , the resulting position lies just inside of the half - light radius of ic  1613 , which is quite fortunate .
ideally , we would like to characterize the global sfh of ic  1613 from the observed field .
stellar population gradients are universal in dwarf galaxies and have been observed over a wide range of physical characteristics . in all cases ,
the gradients are in the sense that the mean age of the stellar population is younger toward the center of the galaxy ( see discussion in * ? ? ?
* for further details ) .
since the observed field is close to the half - light radius , we make a type of `` mean value '' assumption and assume that the sfh derived from the observed field is a relatively good representation of the global sfh for ic  1613 . the balance between outer fields ( with relatively more older stars ) and inner fields ( with relatively more younger stars ) should not be too different from the sfh of a field at the half - light radius . for the rest of the paper
the discussion will be based on this assumption and the sfh history of the observed field and the global sfh of ic  1613 will be used interchangeably .
this assumption is supported by the comparison of hst observations of ic  1613 from different fields presented in  [ comp3 ] .
note that our main result , the absence of a dominant star formation episode at early times , would only strengthen by the addition of more younger stars at smaller radii .
the opposite possibility ( a relative surplus of older stars in the inner parts of the galaxy ) is constrained by the direct observation of the number of rgb stars in the central field from the earlier hst / wfpc2 observations ( * ? ? ? * see discussion in section  [ comp3 ] ) .
following the other lcid studies , we first use the iac method , consisting of iac - star @xcite , iac - pop @xcite , and minniac @xcite in order to solve for the sfh and amr
. details of the methodology can be found in @xcite .
specifically , the solution has been obtained using the basti @xcite stellar evolutionary libraries and the bolometric corrections of @xcite , and the uncertainties are the statistical uncertainties and do not included certain systematics ( e.g. , the use of a different stellar library ) .    , of ic  1613
derived as described in the text .
the sfh as a function of age @xmath34 and metallicity @xmath35 are shown projected on the @xmath36 and @xmath37 planes , respectively .
dashed lines give the statistical error intervals .
the age - metallicity relationship is the projection onto the look - back time
metallicity plane .
[ f3],width=325 ]    figure [ f3 ] shows a plot of the sfh of ic  1613 , @xmath38 , as a function of both @xmath39 and @xmath40 .
@xmath34 and @xmath35 are also shown in the @xmath41 and @xmath42 planes , respectively .
the projection of @xmath38 on the age - metallicity plane shows the amr , including the metallicity dispersion .
figure [ f3 ] shows the presence of star formation at all times , and the gradual build up of metallicity with time .
note especially the relatively constant sfr over the first 6 gyr of the history of ic  1613 ( variations of @xmath2 30% or less in sfr ) .
while the iac - pop solution suggests a slight ( factor of 2 ) decrease in average sfr about 8 gyr ago ( @xmath43 @xmath44 1 ) , there is no evidence for a dominant episode of star formation before the end of reionization ( @xmath2 12.8 gyr , or @xmath43 @xmath44 6 , * ? ? ?
* ) or around the age of the peak sfr density of the universe ( @xmath2 10 gyr , or @xmath43 @xmath44 2 , * ? ? ?
we emphasize that there is no evidence of a dominant episode of star formation in very early times ( lookback time @xmath12 11 gyr ) , despite our ability to resolve such a feature .
( upper panel ) , metallicity ( middle panel ) , and cumulative stellar mass fraction ( lower panel ) of ic  1613 .
thin lines give the uncertainties .
vertical dotted - dashed lines indicate the times corresponding to the 10@xmath45 , 50@xmath45 and 95@xmath45 percentiles of @xmath34 , i.e. , the times for which the cumulative fraction of mass converted into stars was 0.1 , 0.5 and 0.95 of the current total .
the cumulative total mass of stars formed is @xmath46 @xmath47 m@xmath48 .
the mean metallicity , @xmath49[m / h]@xmath50 , is @xmath51 .
the average star formation rate , @xmath52 , is @xmath53 m@xmath48 yr@xmath3 .
a redshift scale is given in the upper axis .
[ f4],width=325 ]    in figure [ f4 ] we show @xmath34 , the amr ( with their associated errors ) , and the cumulative stellar mass fraction of ic  1613 , as a function of time .
three vertical dashed lines show the ages of the 10th , 50th , and 95th - percentile of the integral of @xmath34 .
the main features of the sfh of ic  1613 are nearly continuous star formation and an accompanying smooth increase in the stellar metallicity .
the sfr varies by at most a factor of two from its mean value over the entire lifetime of the galaxy .
most importantly , there is no strong initial epoch of star formation , as is seen in the lcid dsph galaxies and transition ( dsph / dirr ) galaxies .
this is reflected in the cumulative stellar mass fraction plot which is close to the diagonal line of a constant sfr .    .
the residuals are in units of poisson uncertainties.[f5],width=325 ]    figure [ f5 ] shows the observed cmd , the best - fit cmd , and the corresponding residuals .
residuals are given in units of poisson uncertainties , obtained as @xmath54 , where @xmath55 and @xmath56 are the number of stars in bin @xmath57 of an uniform grid defined on the observed and calculated cmds , respectively .
because the solution is heavily determined by the main sequence stars ( i.e. , the rgb and hb stars are not used in the solution , see discussion in * ? ? ?
* ) , the residuals along the ms are negligible .
there are significant residuals at the faintest levels ( where the incompleteness is large ) and in the horizontal branch ( which is typically not well fit by models ) .
, width=325 ]    as for the other galaxies in the lcid sample , we have obtained the sfh of ic  1613 using the match @xcite and cole @xcite methods in addition to the iac method .
these two methods use @xcite for the stellar evolution library for this comparison ( for a description of the main features of these methods and the particulars of their application to lcid galaxies , see * ? ? ?
in this exercise , we have allowed all parameters ( e.g. , distance , reddening , etc . ) to be solved for optimally independently .
that is , we are not trying to compare codes , but rather , trying to make an assessment of the systematic errors that may arise due to choices of different stellar libraries , codes , etc .
there are three main sources of uncertainties in deriving sfhs ( see , e.g. , discussion in * ? ? ?
* ) : the input observational uncertainties , the statistical uncertainties of the solution , and the systematic uncertainty of the limitations of our knowledge of stellar evolution . in a specific test of creating synthetic photometry from one stellar library and using a different stellar library to derive a sfh , @xcite showed that the uncertainties associated with our knowledge of stellar evolution dominated those of the numerical methodology or the observation uncertainties .
they recommend always using more than one stellar evolution library when analyzing a real population to test for these effects . @xcite and
@xcite have also shown that , for sufficiently deep photometry , the stellar evolution uncertainties , as approximated by differences between stellar evolution model and their resulting libraries , represent the dominant systematic uncertainty for derived sfhs .
( see also the discussion in * ? ? ? * concerning the effects of varying heavy element abundance patterns . )
thus , using different stellar libraries is integral to this comparison .
the sfh derived using match features uncertainties calculated following the prescriptions in @xcite . in this case
, random uncertainties were generated using a hybrid monte carlo ( hmc ) process @xcite , with implementation details as described by @xcite .
the result of this markov chain monte carlo routine is a sample of 10,000 sfhs whose density is proportional to the probability density .
( that is , the density of samples is highest near the maximum likelihood point . )
upper and lower random error bars for any given value ( e.g. , cumulative stellar mass fraction at a particular point in time ) are calculated by identifying the boundaries of the highest - density region containing 68% of the samples , with the value 68% adopted as it is the percentage of a normal distribution falling between the @xmath11 1 @xmath58 bounds .
figure [ f6 ] shows the cumulative stellar mass fraction of ic  1613 as obtained with the three different methods .
the agreement between the three methods is very good and all three are best described as consistent with nearly constant star formation .
the one small difference between the sfhs is the offset at intermediate ages between the solution based on the basti models and the two based on the padua models .
a comparison using the same code with two different stellar libraries shows that about half of this difference is due to the differences between the stellar libraries ( which have differences in physical inputs and differ in the treatment of core convective overshoot during the central h burning stage ) and the other half is likely attributable to how different sections of the cmd are binned and weighted .
this comparison strongly supports the conclusion that the star formation in ic  1613 has been relatively constant for the lifetime of the galaxy .
the similar behavior of the three solutions allows us to be confident that our conclusions are independent of the stellar evolution libraries and sfh solution method .
most importantly , an early dominant star formation event is clearly ruled out .
note also that ic  1613 has clearly formed the vast majority of its stars after the epoch of reionization which occurred in the first gyr of the history of the universe ( cf .
* ) .      , width=325 ]
, wfpc2-center corresponds to the 9 orbit observation presented in @xcite , wfpc2-outer corresponds to the 24 orbit observation presented in @xcite , and acs - outer corresponds to new acs observations presented here .
note that the deeper acs observations have allowed significantly smaller uncertainties at the earliest times .
all three sfhs are consistent with a constant sfr over the lifetime of the galaxy , but only the deeper acs observations clearly rule out a dominant episode of early star formation . [ f8],width=325 ]    relatively deep imaging of ic  1613 has been obtained with the hst three times : a central field obtained with the wfpc2 camera ( 9 orbits ) with a depth in i - band of @xmath2 26.5 mag .  presented by @xcite , a deep ( 24 orbits ) outer field ( radius @xmath44 1.4 kpc ) obtained with the wfpc2 camera with a depth in i - band of @xmath2 27.3 presented by @xcite , and our new imaging of an outer field ( radius @xmath44 1.2 kpc ) obtained with the acs camera with a depth in i - band of @xmath2 28.7 .
figure [ f7 ] presents a comparison of the photometry for these three observations . in making figure [ f7 ] ,
identical quality cuts were made on the data , so that they are directly comparable ( although the wfpc observations were obtained with an f555w filter and the present acs observations were obtained with an f475w filter ) .
figure [ f8 ] shows the cumulative stellar mass fraction of ic  1613 as obtained from the three different observations .
the sfhs for the wfpc2 fields are newly derived from the archival data @xcite .
all three sfhs were obtained using the match code @xcite and the padua @xcite stellar evolution library .
the errorbars were derived following the prescription outlined in @xcite .
this comparison to sfhs derived from shallower data needs to be approached carefully .
clearly , the larger uncertainties inherent in the shallower observations are reflected by larger errorbars on the individual points .
however , since the dominant component of the uncertainties are systematic , and since these can only be estimated and not calculated directly , these should be interpreted as `` best estimates '' .
predicting early sfhs from photometry which does not reach the oldest main sequence turn - offs is vulnerable to systematic biases .
there is no substitute for the deep observations with the resulting strong constraints on the old age sfh .
nonetheless , this comparison allows us to test our assumption of the representative nature of the acs observation for the global sfh of ic  1613 .
the main impression from figure [ f8 ] is that the sfhs for all three observations are best described as consistent with nearly constant star formation , but only the deeper acs observations clearly rule out a dominant episode of early star formation .
thus , figure [ f8 ] does provide strong support for the assumption that the sfh derived from the acs observations can be assumed to be representative of the global sfh for ic  1613 .
given the caveats from the previous paragraph we will refrain from any further discussion of this comparison .
in this section , we will compare the sfh of ic  1613 with that of leo  a ( the other dirr galaxy of the lcid sample , * ? ? ? * ) , phoenix and ( the two dsph / dirr transition galaxies , * ? ? ? * ; * ? ? ? * ) and with those of cetus and tucana ( the two dsphs , * ? ? ? * ; * ? ? ? * ) . note that the lcid sample covers a range in dynamical masses . from @xcite the dynamical masses of ic  1613 ,
leo  a , lgs-3 , cetus , and tucana are : 1.1 @xmath11 0.2 @xmath18 10@xmath59 @xmath8 , 1.5 @xmath60 @xmath18 10@xmath61 @xmath8 , 2.7 @xmath62 @xmath18 10@xmath61 @xmath8 , 4.0 @xmath63 @xmath18 10@xmath61 @xmath8 , and 7.1 @xmath11 1.2 @xmath18 10@xmath61 @xmath8 , respectively . for phoenix ,
we derive a value of @xmath22.7 @xmath18 10@xmath61 @xmath8 , assuming a velocity dispersion of @xmath28 km s@xmath3 @xcite and a half - light radius of 454 pc @xcite .
, width=325 ]    figure [ f9 ] compares the sfh and the amr of ic  1613 with the other lcid galaxies . in the top panel of figure [ f9 ] ,
the sfrs ( @xmath64 ) have been normalized such that the area under each curve is 1 ( i.e. , the integral of the sfr over the lifetime of the galaxy is 1 ) .
the differences between the sfhs as a function of morphological type is quite remarkable .
both dirr galaxies lack the dominant initial episode of star formation seen in both the dsph and transition type galaxies .
although the sfhs are the main topic of this study , the amrs are also solved for in deriving the sfhs and are included here for completeness .
both dirr galaxies also show relatively little metal enrichment during the first @xmath2 6 gyr , with gradual enrichment following .
note that @xcite highlighted the lack of early chemical enrichment in the transition galaxy lgs-3 and hypothesized that the differences between the transition galaxies and dsphs ( which show early chemical enrichment ) might be understood if the dsphs were initially more massive systems than the transition galaxies . under this hypothesis ,
the delayed metal enrichment in the transition galaxies is due to higher losses of metals during times of higher sfrs for the lower mass transition galaxies . in this light , it is interesting that ic  1613 ( and also leo  a ) show delayed metal enrichment , but with no initial high sfr . in this regard , further comparisons of the sfh and amrs of transition galaxies and dirr galaxies would be be great interest .    ,
width=325 ]    before comparing to theoretical models , we introduce one last comparative figure . in figure
[ f10 ] , we show a comparison between the sfhs of the six lcid galaxies as cumulative stellar mass fractions . in the upper panel ,
only the statistical uncertainties are shown . in the lower panel , the systematic uncertainties are included following the methodology of @xcite . because all of the lcid galaxies have been observed to comparable depth , systematics in the models should have similar impacts to all of the galaxies .
thus , it is likely appropriate to make comparisons using the upper panel .
nonetheless , the lower panel shows the larger uncertainties encountered when trying to account for systematics and , even with the larger uncertainties , the six galaxies are shown to each have distinctive features in their sfhs .
portraying the sfhs as cumulative stellar mass fractions ( as opposed to sfr as a function of time ) is the optimal way to compare observations to theoretical models for several reasons . variations in observed sfrs can be strongly affected by time binning and the changing time resolution as a function of lookback time . often , it is possible to have very different impressions of a single sfh simply by changing the time binning .
it is possible to match the observational time binning by reducing the resolution in the models , but using the cumulative stellar mass fraction as the diagnostic avoids this problem altogether .
it is also possible to compare galaxies at any arbitrary value of the cumulative stellar mass fraction , as opposed to choosing particular values to focus on . in the comparisons that follow
, we will use the cumulative stellar mass fraction as the sole diagnostic .
note that there is one obvious failing of the cumulative stellar mass fraction as the sole diagnostic , and that is the lack of information about the absolute masses of the systems .
in the following comparisons , we will provide information about the masses of both the observed and modeled systems .
recently , @xcite , @xcite , @xcite , and @xcite have all highlighted the difficulty of producing dwarf galaxies in simulations with properties comparable to those observed in the current universe .
together , the introductions to their papers give a comprehensive overview of the problems with modeling dwarf galaxy evolution .
to summarize , there are two major problems .
the first problem is the observed abundance of low - mass galaxies .
the observed slope of the low - mass galaxy luminosity function is shallow relative to the slope of the halo mass function and this difference seems to be a result of an extreme inefficiency of galaxy - scale star formation over cosmic times .
the well known problem of `` the missing satellites '' @xcite is one manifestation of this problem at very low luminosities / halo masses .
the second problem relates to _ when _ stars are formed in galaxies .
a natural assumption is that the time - scale for global star formation is related to the time - scale of baryonic accretion onto galaxies .
however , low - mass halos assemble almost all of their ( dark matter ) mass at high redshift ( e.g. , * ? ? ?
* ) , while essentially all field dwarf galaxies show star formation continuing to the present day ( e.g. , * ? ? ?
ic  1613 and leo  a are extreme examples of this , with essentially a constant sf rate across cosmic time for ic  1613 and delayed star formation in the case of leo  a.    in order to suppress the abundance of low - mass galaxies , most theoretical models impose strong feedback in small halos ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
two processes can dramatically affect the formation and evolution of dwarf - sized halos : heating from the ultraviolet radiation arising from cosmic reionization and feedback from internal supernovae .
both processes are , in principle , capable of completely halting the star formation in a dwarf halo and even fully removing all of the galaxy s gas . employing these feedback mechanisms
while tying star formation to the collection of baryons has the effect of predicting that essentially all star formation in low - mass halos happens at early times .
both semi - analytic models and hydrodynamical models fail to satisfactorily reproduce the evolution of low mass galaxies in the sense that stellar mass is over - produced per dark matter halo mass .
stars are produced too quickly at early times , resulting in stellar mass fractions that are too high by an order of magnitude .
model galaxies usually do not have the high gas mass fractions commonly observed in present day dwarfs ( e.g. , * ? ? ?
* ) , and , as a result , star formation falls off too fast and the colors of simulated dwarf galaxies are too red at @xmath43 @xmath7 0 .
this is a manifestation of the `` over - cooling '' problem which is a challenge for all modeling efforts , but which is exaggerated at low masses .
the missing satellites problem is largest at the lowest masses , and ic  1613 is massive enough that it could be expected to emerge from reionization with its gas intact .
thus , ic  1613 does not provide as strong a test of models for this problem as the other , less massive , galaxies in the lcid sample .
however , our deep hst observations have produced sfhs for ic  1613 and leo  a with small uncertainties even at earliest times and have completely ruled out the possibility of an early , dominant phase of star formation .
_ this represents a real challenge for any model where star formation follows mass assembly .
_       @xmath7 0 of @xmath2 10@xmath65 m@xmath66 , which is comparable to the estimated halo mass for ic  1613 .
the difference in the early sfhs between the models and ic  1613 is very clear .
, width=325 ]    @xcite have highlighted the shortcomings of models for dwarf galaxies where star formation follows directly from mass assembly .
they presented a set of high - resolution hydrodynamical simulations of the formation and evolution of isolated dwarf galaxies including the most relevant physical effects , namely metal - dependent cooling , star formation , feedback from type ii and ia sne , and uv background radiation .
their results are very useful for a direct comparison with observations .
they identify and study six halos with present day dark matter halo masses of @xmath44 10@xmath65 m@xmath48 with stellar masses ranging between @xmath67 m@xmath48 and @xmath68 m@xmath48 . for comparison ,
ic  1613 has an estimated stellar mass of @xmath2 10@xmath59 m@xmath48 @xcite , and the observed maximum rotation velocity of 25 km s@xmath3 @xcite . if we associate this peak rotation velocity with the maximum circular velocity of ic  1613 s host dark matter halo , this corresponds to a virial mass .
] of @xmath69 .
this is likely to be a lower limit , however , as the peak of the circular velocity curve for the dark matter halo may be attained at radii beyond the extent of the stellar or gaseous tracers .
abundance matching models find that galaxies with stellar masses of @xmath70 should be hosted by halos with virial masses of @xmath71 ( e.g. , * ? ? ?
* ) . given these estimates , the models of sawala et al .
provide an excellent sample to compare to the derived sfh of ic  1613 .    as figure [ f11 ] demonstrates clearly , the early evolution of the model galaxies all depart significantly from that of ic  1613 .
in essence , the models all have sfrs which depend directly on the gas content of the galaxies . as highlighted by @xcite , since , in the current paradigm , most of the mass of a low mass galaxy is in place well before @xmath43 @xmath7 1 ( e.g. , * ? ? ? * ) , then any prescription in which star formation follows gas content is going to build up most of the stellar mass before @xmath43 @xmath7 1 .
this discrepancy between models and observations , while demonstrated here for ic  1613 , is beginning to appear to be the norm .
although the numbers are very limited , so far , all gas - rich dwarfs for which there exist sufficiently deep hst observations ( i.e. , leo a , leo t , smc : * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) show no evidence for a dominant early episode of star formation .
there have been various attempts to improve models to better accommodate the inefficient star formation in dwarf galaxies . @xcite
used sph + n - body simulations and showed that supernova feedback could disrupt enough gas to temporarily quench star formation .
episodic star formation follows from the cycling of gas into a hot halo with subsequent cooling and infall .
@xcite account for the prolonged star formation observed in dwarfs by adding a stochastic variation in the density threshold of the star formation law .
essentially , this simply reduces the efficiency of star formation by hand , but the result is a significantly improved match to the observed sfhs of dwarfs .
recently , @xcite have demonstrated that the thermal feedback from early star formation can effectively decouple star formation from mass assembly , thus producing more realistic sfhs .
@xcite and @xcite have proposed a different physical approach by investigating the conversion of atomic to molecular gas and its affect on the efficiency of star formation .
this is well motivated by the observational work of @xcite and @xcite , which clearly demonstrate that star formation follows the molecular gas content and not the total gas content . in the specific case of dwarf galaxies , @xcite
show that models which incorporate a star formation prescription regulated by the local abundance of molecular hydrogen lead to the less efficient star formation that is desired .
unfortunately , their first attempts are not completely successful as they state , `` like most cosmological galaxy formation simulations to date , our simulated galaxies suffer from the so - called baryonic overcooling problem , resulting in unrealistically high central densities ( and hence strongly peaked circular velocity curves ) and stellar mass fractions in our high - mass halos that are too large compared to observations . ''
however , this appears to be a very promising avenue for future exploration ( see also * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
@xcite used a semi - analytic galaxy formation model to investigate the properties of the satellites of milky way - like galaxies and were able to match the star formation histories of several dwarf satellites .
the extended sfhs of these satellites were a result of a gas density threshold for star formation .
galaxies could have large reservoirs of gas lying just below the threshold , leading to inefficient star formation ( requiring a minor accretion or interaction event ) .
@xmath7 0 of @xmath2 4 @xmath18 10@xmath65 m@xmath66 and @xmath2 1 @xmath18 10@xmath65 m@xmath66 , respectively , comparable to the estimated halo mass for ic  1613 .
`` bashful '' and `` doc '' have nearly constant star formation and lack a dominant early episode of star formation as is observed in ic  1613 .
the two model galaxies with a delayed onset of star formation are `` dopey '' and `` grumpy '' with smaller halo masses at @xmath43 @xmath7 0 of @xmath2 3 @xmath18 10@xmath72 m@xmath66 and @xmath2 2 @xmath18 10@xmath72 m@xmath66 , respectively , ( which are considerably less massive than ic  1613 , but might be more directly comparable to leo  a ) .
, width=325 ]    @xcite have produced new simulations including supernova feedback , a star formation recipe based on a high gas density threshold , metal - dependent radiative cooling , turbulent diffusion of metals and thermal energy , and uv background radiation .
they reproduce the observed stellar mass and cold gas content , the sfhs , and the metallicities typical of field dwarfs in the local volume . in figure
[ f12 ] , we reproduce the cumulative stellar mass fraction for four of the model galaxies in @xcite .
two of the galaxies ( `` bashful '' and `` doc '' ) have virial masses which are comparable to ic  1613 .
these two models show nearly constant star formation with no dominant early episode of star formation , as observed in ic  1613 . in detail , the sfhs of these two models start later and rise faster during intermediate ages .
the slower start of the models as compared to that of the measured sfhs may be attributed to differences in how time is defined in the two methods .
that is , the simulations start at a fixed , given time , whereas the observationally derived ages come from stellar evolution libraries , and therefore depend on their underlying physics ( e.g. , nuclear reaction rates ) .
the two lower mass galaxies ( `` dopey '' and `` grumpy '' ) show no early star formation at all , and only start forming stars at intermediate ages .
these two models correspond to halo masses significantly less than that of ic  1613 , but might be more comparable to that of leo  a ( cf .
figure [ f10 ] ) .
note , however , that the stellar masses are even more discrepant .
that is , dopey and grumpy have stellar masses of 10@xmath73 and 5 @xmath18 10@xmath73 m@xmath66 compared to 10@xmath59 and 6 @xmath18 10@xmath74 m@xmath66 for ic  1613 and leo  a.    naively , the late onset of star formation in the dopey and grumpy models would appear to be in contradiction to observations in which all dwarf galaxies observed to date show evidence for an old ( @xmath12 10 gyr ) population , particularly as evidenced by the presence of rr lyrae stars .
for example , rr lyrae stars have been observed in leo  a @xcite .
however , @xcite hypothesize that such systems may be the very metal poor systems necessary for the production of extremely metal deficient galaxies such as i  zw  18 ( see , e.g. , discussion and references in * ? ? ?
one of the successes of recent modeling of dwarfs has been the ability of simulations to solve the cusp / core problem for dwarf galaxies ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
the dark matter halos of observed dwarf galaxies show nearly constant density cores , while steep central dark - matter profiles are expected from cdm models ( e.g. , * ? ? ?
the resolution to this conundrum lies in strong outflows from supernovae to remove low - angular - momentum gas from the centers of dwarf galaxies .
this solves two problems as bulgeless dwarf galaxies are formed with shallow central dark - matter profiles , and a large fraction of the original baryons are lost from the galaxy , resulting in the lower baryon fractions ( relative to the universal fraction ) which are observed ( e.g. , * ? ? ?
what is not clear is whether the current ideas concerning how to produce inefficient star formation also lead to cored dark matter halos .
clearly the strong outflows need to occur when the mass assembly is taking place , which is early on .
if the vast majority of the star formation is delayed until later times , after the mass assembly is virtually finished , then the central outflow solution to the cusp / core problem may not arise self - consistently .
a related point is the energetic requirements for creating bulgeless galaxies and forming cores in state - of - the - art hydrodynamic simulations .
simulations that are currently successful require ( 1 ) highly efficient coupling between the energy released in a supernova explosion and the ism of a galaxy  values typically adopted are 0.41 @xcite and ( 2 ) an artificial prevention of gas cooling for approximately @xmath75 years in a region surrounding a supernova explosion , meant to capture processes that may occur at sub - grid scales and can not be directly included in the simulation @xcite .
the success of simulations in producing bulgeless , cored dwarf galaxies with relatively low - level , extended star formation can be therefore seen as predictions for the effective coupling of supernova explosions with the ism and for the efficacy of additional feedback processes ( such as stellar winds ) at heating the ism on small scales .
these predictions of very efficient coupling may be confirmed or refuted with future generations of simulations that model the relevant processes directly .    finally , although we have emphasized the difficulty of producing a nearly constant sfr as seen in ic  1613 , observations of dwarf galaxies in the local volume point toward a large diversity of sfhs .
determining the critical parameters which drive that diversity is paramount .
how can it be that dwarfs with halos of similar mass , which presumably have similar accretion histories and early formation times , have vastly different sfhs ?
is environment and interaction history the main driver of the diversity ( see , e.g. , * ? ? ?
* ; * ? ? ?
* ) ? from an observational point of view , strong constraints will depend on accurate determinations of the early sfhs of a larger number of galaxies within and beyond the local group from both deep photometry as presented here and spectroscopic studies of individual stars ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
clearly the sample of one presented here , and the very small number of comparably observed galaxies presents a great limitation to our knowledge of the evolution of dwarf galaxies .
we have presented the sfh of the dirr galaxy ic  1613 , based on deep hst photometry obtained with the acs .    *
the sfh with relatively small uncertainties has been obtained for the entire lifetime of the galaxy .
the solution shows that the sfh of ic  1613 is consistent with a constant sfr over its entire lifetime . *
most or all the star formation was produced in ic  1613 after the reionization epoch , assumed to occur @xmath76 gyr ago .
there is no evidence of an early dominant episode of star formation in ic  1613 . *
a comparison of the derived sfh of ic  1613 with the models of @xcite reinforce their observation that models where star formation follows mass assembly form too many stars too early .
this well known aspect of the so called `` over - cooling problem '' appears to be universal , but now , with the deep hst photometry of ic  1613 , we see that the problem reaches back to the very earliest times in the evolution of galaxies .
* there are proposed solutions to the over - cooling problem for dwarf galaxies .
the solutions discussed in this paper rely on an efficient coupling of supernova feedback to the ism .
the predictions of very efficient coupling may be confirmed or refuted with future generations of simulations that model the relevant processes directly and through future observations of dwarf galaxy sfhs with similar quality to the lcid studies .
we would like to thank till sawala for sharing the numerical results of his modeling , and greg stinson , piero madau , alyson brooks , and the referee for helpful comments .
support for this work was provided by nasa through grant go-10515 from the space telescope science institute , which is operated by aura , inc .
, under nasa contract nas5 - 26555 .
support for drw is provided by nasa through hubble fellowship grant hst - hf-51331.01 awarded by the space telescope science institute .
the computer network at iac operated under the condor software license has been used .
authors sh , aa , cg , and mm are funded by the iac ( grant p3/94 ) and , with sc , ps , and eb , by the science and technology ministry of the kingdom of spain ( grant aya2007 - 3e3507 ) and economy and competitiveness ministry of the kingdom of spain ( grant aya2010 - 16717 ) .
this research has made use of nasa s astrophysics data system bibliographic services and the nasa / ipac extragalactic database ( ned ) , which is operated by the jet propulsion laboratory , california institute of technology , under contract with the national aeronautics and space administration .