blue straggler stars ( bsss ) were first identified in the galactic globular cluster ( ggc ) m3 by @xcite , as a group of stars bluer and brighter than the cluster turn - off stars , thus being more massive than the stars currently evolving at the turn - off ( to ) . since then
, bsss have been identified in a variety of stellar systems in very different environmental conditions including : the field @xcite , open clusters @xcite , globular clusters ( e.g. , * ? ? ?
* ) , and dwarf galaxies @xcite .
their puzzling position in the color - magnitude diagram ( cmd ) suggested that they do not fit in the traditional scheme of stellar evolution for single stars .
in fact , in order to populate that region of the cmd , these exotic objects must have experienced a physical process able to allow them to still stay in the core h - burning stage despite of their mass and the cluster age ( for a detailed historical review , see * ? ? ? * ) .
presently , two different physical mechanisms are favored for explaining the formation of bsss : _
i ) _ coalescence through direct stellar collisions ( col - bs , * ? ? ?
* ) , and _
ii ) _ evolution of primordial binary systems , in which mass transfer between the two components allows the rejuvenation of the secondary component ( mt - bss , * ? ? ?
it is also widely accepted that these two scenarios are not in competition , but they could both occur simultaneously in a stellar system .
in fact , the coexistence of bsss formed through both the collisional and mass transfer channels is commonly invoked to interpret the radial distribution of bsss in many ggcs .
the higher central concentration is explained by an efficient rate of stellar collisions in the densest regions , while the rise at large distance would be due to mt - bsss @xcite surviving in the outskirts .    while the scenario for bsss in galactic clusters is well settled , at least from the observational point of view , in contrast ,
this is not true for dwarf galaxies .
in fact , even though the first indication for the presence of bss candidates in the sextans dwarf dates back to @xcite , and , since then , several observational findings suggest the presence of bss candidates in several dsphs ( see , e.g. , * ? ? ? * and references therein ) , we are still faced with the problem of understanding if these candidates are _ genuine _ bsss or _
normal _ core h - burning stars belonging to an , albeit sparse , intermediate - age stellar population .
recently , a few studies have focused on the properties of bsss in nearby dsph galaxies in the local group ( lg ) .
@xcite presented wide field data for the galaxies draco and ursa minor , concluding that the population of bsss is , in both systems , compatible with the mt - bss scenario .
@xcite discussed the properties of bsss in sculptor and fornax .
sculptor presents similarities with both draco and ursa minor , as its bss population is also compatible with a mt - bss population . on the contrary
, the relatively strong central radial concentration of the bss candidate stars in fornax favors the presence of an intermediate - age stellar population ( in agreement with previous studies of the star formation history ( sfh ) of fornax ) .
unraveling the nature of the bss candidates has important implications for our understanding of the sfhs of these galaxies .
for example , whenever it is possible to discriminate between a sequence of bsss and a young ms population , it is possible to set a constraint on the age of the last star formation event in a galaxy . therefore , the bss populations in draco , ursa minor , and sculptor suggest that no star formation occurred in the last few gyr in these galaxies .    intriguingly , @xcite found a statistically significant anti - correlation between the frequency of bsss and the absolute visual magnitude ( m@xmath0 ) of 8 nearby dsphs .
the anti - correlation is valid over a large range of luminosities , and presents a different slope than that of ggcs .
in addition , @xcite ruled out the possibility that collisional binaries could have contributed to the observed bss populations and concluded that the bsss in dwarf spheroidal galaxies are mainly formed via mass exchange in primordial binaries .    in view of the important implications of the presence of a sizable sample of bsss in a dwarf galaxy , it is worthwhile to apply the approach outlined by @xcite to other dwarfs .
the hope is to eventually produce independent constraints useful for addressing the true nature of the bss candidates .
l|cc @xmath1 & @xmath2 & @xmath3 + @xmath4 & @xmath5 & @xmath6 + n@xmath7 & @xmath8 & @xmath9 + n@xmath10 & @xmath11 & @xmath12 + f@xmath13 & @xmath14 & @xmath15 + n@xmath16 & @xmath17 & @xmath18 +    in this paper , we investigate these issues further by analyzing the properties of candidate bsss in the isolated dsphs cetus and tucana .
the data sets are part of the lcid project ( _ local cosmology from isolated dwarfs _ ) , aimed at recovering the full sfhs of six isolated lg galaxies : cetus and tucana ( dsph ) , lgs  3 and phoenix ( dirr / dsph ) , ic  1613 and leo  a ( dirr ) .
the plan of the paper is as follows : in the next section we briefly describe the photometric data set and outline the criteria for selecting the candidate bsss in the observed cmds . in
3 we investigate the properties of the bss populations .
 4 presents the results of dynamical simulations .
the discussion and final remarks close the paper .
the observations used in this paper have already been presented in ( * ? ? ? * cetus ) and ( * ? ? ? * tucana ) . here
we summarize the main points relevant to the following discussion .
the images were collected with the acs camera aboard the hst , under the project _ the onset of star formation in the universe : constraints from nearby isolated dwarf galaxies _ ( pid 10505 , pi c. gallart ) .
a total of 25 and 32 orbits have been devoted to cetus and tucana , respectively , using the @xmath19 and @xmath20 filter pass bands .
the photometric reduction has been performed using the daophot / allframe package @xcite .
the final catalog has been calibrated to the vegamag system using the transformations presented in @xcite .
the adopted distance modulus and reddening are summarized in table [ tab : tab01 ] , together with other quantities calculated in this work .
figure [ fig : samples ] shows the cmds of the two galaxies , with the regions adopted to select the various samples studied superimposed . note
that we adopted the same boxes for both galaxies , applying a shift to take into account the differences in distance modulus and reddening .
the cmds of both galaxies are @xmath21 mag deeper than the to , reaching @xmath22 mag .
the blue plume of candidate bsss emerges clearly in the cmd , for @xmath23 mag , @xmath24 .
we defined a box in this region to include the bulk of these stars .
note that there are a number of objects between the bsss region and the red giant branch ( rgb ) .
we decided not to include these stars for two reasons : _
i ) _ it is possible that blended objects from the most populated to region are polluting this part of the cmd ; _ ii ) _ in this color range , we expect to find stars evolved off the main sequence
. also , we did not include the brightest objects ( @xmath25 ) because it is possible that they are extreme horizontal branch ( hb ) stars .
also note that we estimated a negligible contamination from both foreground galactic stars ( due both to the small area covered and the high galactic latitude of the two galaxies ) and background galaxies ( thanks to the careful cleaning of the catalogs , see @xcite ) .
figure [ fig : samples ] also shows the regions adopted to select stars in the hb , the rgb , and the progeny of bsss during the central helium burning phase .
this last box has been placed above the red hb , paying attention not to include the agb clump , which clearly appears in the cmds of both galaxies ( @xmath26 mag ) .     0.8 cm
in the box covering the blue plume , we have identified 940 and 1214 objects in cetus and tucana , respectively , which we consider as bss candidates . note that the present data cover a significant fraction of the body of tucana ( tidal radius @xmath27 ) , which was centered in the acs field ( field - of - view @xmath28 ) .
however , cetus occupies a larger area in the sky ( @xmath29 , @xcite ) , and the sampled field covers a smaller fraction .
moreover , the observed field was offset from the center so that the innermost regions of cetus were not observed @xcite .
if we take this into account , following the reasoning of @xcite to estimate the total number of rr lyrae stars , the total number of bss candidates in cetus is estimated to be @xmath30 5,000 stars .
the most important question is to assess the nature of these stars .
bsss have been found in very different environments , from the field , to low - density stellar systems such as open clusters and dwarf spheroidals , and in most globular clusters .
they are therefore common products , and it is reasonable to expect to find them in tucana and cetus as well .
therefore , the basic question we are trying to address is : are these candidates genuine bsss , or is there also a significant component of truly young stars ? in the following sections we will try to address this using different observables .    0.4 cm , but we present a zoom of the cmd of cetus ( left panel ) and tucana ( right panel ) , corresponding to the core he - burning stage . in each panel
the box represents the region used to select the candidate progeny of the blue plume stars .
four tracks are overplotted , for different mass and metallicity assumptions .
the filled stars mark the location of the anomalous cepheid variables detected in these galaxies @xcite .
, title="fig:",width=321 ]      the main objective of the lcid project was to derive accurate sfhs for a sample of six isolated galaxies in the local group , in order to understand their evolution in a cosmological context . to do this
, we adopted well - established techniques based on the comparison of the observed cmd with a synthetic one .
we devoted significant effort to investigating the error budget and also possible systematics affecting our solution , due to the use of different photometric packages , stellar evolutions libraries , sfh codes . the main conclusion
( see * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) of this analysis was that we can rely , within the error bars , on both the uniqueness and stability of the derived solutions .
therefore , the most important features recovered , such as the epoch and the duration of the main star formation events , including those at the oldest epoch , and the age - metallicity relations are solid results .
the sfhs presented in @xcite and @xcite , for cetus and tucana respectively , show that both galaxies are made of old and metal - poor populations .
in particular , both galaxies formed 90% of their stars at epochs older than 9  gyr ago , and the bulk of the stars have metallicities in the range @xmath31 .
in addition , both galaxies present a well defined age - metallicity relation in the sense of increasing metallicity with decreasing age .
the sfh reconstructions indicate a small population ( @xmath32 in mass ) of stars significantly younger than the majority of stars ( @xmath33 gyr ) but with metallicities consistent with the oldest stars ( @xmath34 ) .
thus , these stars do not follow the general age - metallicity relations in these two galaxies .
the analysis of the best - fit cmds show that they nicely populate the blue plume region of the cmd , brighter and bluer than the oldest ms turn - off .
the possibility that the metallicity of these stars has been understimated by the code is excluded with high confidence level .
in fact , higher metallicity and an age of @xmath35 gyr , would produce also a redder rgb sequence that is not observed . on the basis of the retrieved sfh
, @xcite made the hypothesis that the sequence of blue stars , brighter than the old ms turn - off , are consistent with a population of bsss , characterized by the same metal content as the metal - poor dominant stellar component .
an alternative interpretation would be that these are truly young metal - poor ms stars , formed as a consequence of gas accretion .
however , it seems very unlikely that two galaxies , that spent most of their life in isolation , accreted gas and formed stars of the same metallicity and at the same epoch , in one event that created roughly the same amount of mass ( 3% of the total stellar mass ) .
figure  [ fig : track ] shows the comparison of the tucana cmd with theoretical tracks from the basti database of stellar evolution models @xcite .
here we note that the analysis of the sfh for the lcid galaxies has been performed by using the same evolutionary framework ; so the present analysis is fully consistent with those concerning the sfhs .
the selected tracks span the range of metallicities suggested by the sfhs , z = 0.0001 , 0.0003 , 0.0006 , and 0.001 as an upper limit .
the figure discloses some interesting clues .
first , the black lines in the left panel show the zero age main sequence ( zams ) for two metallicities : z=0.0001 ( solid line ) and z = 0.001 ( dashed ) .
the comparison with the observed cmds shows that both zamss border the blue envelope of the bss sequence , independent of the assumed metal content .
moreover , the same panel shows that the maximum mass fitting the brightest objects is equal to about 1.5 @xmath36 ( red lines ) , again with a negligible dependence on the assumed metallicity ( z=0.0001 , 0.001 for the red solid and dashed lines , respectively ) .
the same figure also shows the location of two theoretical tracks corresponding to @xmath37 stellar models ( green lines ) , taken as representative of the stellar structures currently evolving at the to and sub - giant branch of both galaxies .
the filled symbols show the first point corresponding to the central helium burning phase .
it is evident - as also fully supported by the sfh analysis - that for the to and the sgb , a range of metallicities in the adopted evolutionary framework is required to properly interpret the observed morphology of the cmd .
this comparison reveals that the maximum mass of the stars in the blue plume is at most @xmath38 2 times larger than the mass of the stars at the to .
this occurrence is in good agreement with the properties of bsss in other stellar systems such as ggcs ( see , i.e. , * ? ? ?
* and references therein ) .
the central panel shows the comparison with stellar tracks for different masses ( @xmath39 , green and red lines , respectively ) , and different metallicities ( @xmath40-solid , @xmath41-dashed , @xmath42-dot - dashed ) .
the shaded area encloses the region corresponding to the central hydrogen burning phase , whose end is marked by the asterisks .
the bulk of bss candidates clearly occupies this region of the cmd .
the plot clearly shows that a spread in mass is needed to fit the observed distribution of blue plume stars .
a spread in metallicity is also possible but , as shown in the right panel of the same figure , values larger than @xmath43 seem unlikely .
the 1.0 @xmath44 track is too faint and too red to represent the bulk of the blue plume objects .
the 1.2 @xmath44 still gives a good representation , and can be considered as an upper limit both in term in of mass and metallicity .
in fact , despite the 1.5 @xmath44 still limits the brightest blue plume stars , the age of this model at the turn - off ( @xmath381.3gyr ) is in contrast with the finding of the sfh . on the other hand ,
higher metallicities would shift the tracks to redder color and fainter magnitude , significantly worsening the agreement in particular in the rgb phase .
thus , relatively metal - poor and relatively massive ( @xmath45 ) stellar structures are preferred to better represent the peak of the distribution .
therefore , the metallicity derived for the bss candidates is in good agreement with that estimated for the bulk of the older populations , as discussed in
[ sec : sfh ] , and not in agreement with the higher metallicities found at the end of the initial episodes of star formation in cetus and tucana .
this result can be considered a circumstantial evidence supporting the idea that - at least a large fraction of - the stars located along the blue plume could be genuine bsss ; in fact if they were true young stars formed in a late star formation burst , they would have a metallicity lower than that one could expect of the basis of a standard age - metallicity relation .
additional information can be provided by the study of the bss progeny in the core he - burning stage .
evolutionary models for ms stars predict that stars more massive than about @xmath46 are located in the red side of the hb at a magnitude level depending on their he core mass @xcite .
in particular , the right - hand panel of figure  [ fig : track ] suggests that the evolved stars in this mass range all clump in the same region of the cmd .
this implies that the progeny of bsss of similar mass , independent of the formation mechanism ( see , e.g. , * ? ? ?
* ) , should occupy a similar region in the cmd , redder and brighter than their low - mass counterparts .
this means that , on the one hand , the distribution of bss progeny in the cmd provides very little information about the mass distribution of the bss population above @xmath46 . on the other hand ,
the ratio of the number of stars can be compared with theoretical expectations based on the evolutionary lifetimes in the core h- and he - burning stages .
figure  [ fig : trackevol ] shows an expanded view of the hb region of the cetus ( left panel ) and tucana ( right panel ) cmds . superimposed to both plots are the evolutionary tracks , starting with the onset of the central helium burning , for two masses ( 1.0 -blue line- and 1.3 -red-@xmath36 ) and two metallicities ( @xmath40 -solid- and @xmath42 -dashed- ) .
the boxes corresponding to our selections are shown .
both panels show the presence of a sparse sequence above the red hb .
while a few of them , which are located close to theoretical prediction for the most metal - poor stars , fall outside our box , the majority of stars , inside the box , look more compatible with slightly more metal - rich tracks . note that this is in agreement with the sfh results , which gives for the two galaxies a mean metallicity of the order of @xmath47 ( cetus ) and @xmath42 ( tucana ) .    to further investigate the interpretation that the stars located inside the box showed in fig .
[ fig : trackevol ] can be really considered the progeny of the stars located along the ms blue plume , we compared the star counts ratio between the objects located inside the two boxes with theoretical predictions obtained taking into the account the information from the sfh .
* their fig 18 ) have already shown a synthetic cmd built from the sfh solution of cetus , where stars with colors and magnitudes typical of a population of bsss clearly appear .
we therefore counted the stars using the same boxes defined here , and estimated the same ratio using a best - fit cmd for both cetus and tucana .
the derived ratios are 0.08@xmath480.01 and 0.08@xmath480.01 for cetus and tucana respectively .
analogous estimates can be made from the observed cmd .
the number of bsss in cetus and tucana are 940 and 1214 , respectively , while the number of objects in the evolved phase are 90 and 123 .
this gives ratios of 0.10@xmath480.01 for cetus and 0.10@xmath480.01 for tucana .
this shows that there is good agreement with the empirical values .
this suggests that the selected stars are the evolved counterpart of blue plume objects .
moreover , this test strongly suggests that , if we are dealing with a pure bsss population , the analysis performed using models calculated for normal ms stars is only marginally affecting the results .
it is also worth noting that the progeny of bss - in particular the more massive ones - can cross the instability strip for radial pulsation during the core he - burning stage ; in this case they pulsate as anomalous cepheids @xcite .
when studying the populations of variable stars in the lcid galaxies , @xcite discovered three anomalous cepheids in cetus , and six in tucana .
the locations of these variable stars in the cmds of both galaxies are shown in fig . [
fig : trackevol ] .
they appear bluer and brighter than the expected bulk of the bss progeny .
this suggests that anomalous cepheids might have slightly larger masses than the typical bss .
0.4 cm    [ fig : cumul ]      the luminosity function of blue plume stars is a powerful diagnostic to both infer information on their nature and , in case of a bss population , to distinguish between the two main formation mechanisms .
in fact , if these stars belong to an intermediate - age population , we expect to find them more centrally concentrated ( see next section ) and also that the brightest ( youngest ) are preferentially located to the innermost regions . on the other hand
, it has been found that in ggcs the bss luminosity function changes as a function of radial distance @xcite , in the sense that the brightest bss stars are preferentially located in the central regions . since there is no doubt that in gc we only have bsss , with no contamination from intermediate - age stars ,
this occurrence is explained with a combination of two different formation mechanisms .
collisions , which are efficient in the inner regions of ggc , tend to create bluer and brighter bsss @xcite . in the outer parts of the ggcs ,
the bsss are predominantly created from primordial binaries . in dsph galaxies
the formation of col - bss is highly disfavored due to the very low stellar density , so that we do not expect to observe a similar behavior @xcite .
therefore , in case a central concentration of the brightest blue plume stars were observed , this would point to the presence of a genuine intermediate - age population .
figure  [ fig : lumfun_rad ] shows the bss candidates luminosity function for cetus ( left panel ) and tucana ( right panel ) .
the plot shows the luminosity function of stars inside ( dotted line ) and outside ( dashed line ) a distance of 1.5 times the core radius , @xmath49 . due to the different number of stars and area covered by the two samples ,
the curves were normalized to their area .
the shape of the luminosity function does not change significantly as a function of radius , and similarly we do not find any evidence for an obvious shift in the magnitude of the peak .
we find that the fraction of the brightest bsss ( @xmath50 , i.e. , within 0.5 mag of the brightest ) , is also constant with radius .
for the inner and outer regions we found @xmath51 and @xmath52 in the case of cetus , and @xmath53 and @xmath54 for tucana . note that the errors are solely poisson statistical errors .
thus , we find no evidence for segregation in either galaxy , supporting the absence of centrally concentrated bright stars , and supporting that the population of bsss in cetus and tucana is formed as a result of the evolution of binary systems , even in the innermost regions .
this last result can be further verified studying the radial distribution of candidate bsss compared to that of other stellar tracers .
it has been found in many ggcs that the radial profile of bsss has a central cusp ( c.f .
* ) , which falls to a minimum ( the zone of avoidance * ? ? ?
* ) , and then increases with radius .
this phenomenon has again been interpreted with the different formation mechanisms at work : the col - bsss are more centrally concentrated than mt - bss @xcite , which instead dominate in the outskirts of clusters . if col - bsss are not present in dsphs , we expect no central concentration , as verified by @xcite and @xcite for draco , ursa minor and sculptor .
figure [ fig : grad_rad ] presents the ratios @xmath55 , @xmath56 , and @xmath57 for cetus ( left ) and tucana ( right ) , as a function of the galactocentric distance expressed in units of core radius .
we divided the complete catalog into five ( cetus ) and six ( tucana ) radial bins , containing the same total number of stars , according to the position of stars in elliptical annuli .
the points in figure [ fig : grad_rad ] show the ratio in each elliptical radius , plotted as a function of the mean galactocentric distance of each annulus .
figure [ fig : grad_rad ] clearly shows the absence of strong central concentrations of bsss in both galaxies .
we tested the hypothesis that the ratios of the various samples are consistent with a flat distribution by comparing the ratios to the mean of the observed distribution . in the case of cetus
, we estimated non - reduced @xmath58 of 6.6 , 7.3 , 9.3 for @xmath59 , @xmath60 , and @xmath61 , respectively .
the null hypothesis therefore has corresponding probabilities of 0.15 , 0.12 , and 0.05 .
this suggests that @xmath59 and @xmath60 are consistent with a flat distribution , while @xmath61 is only marginally consistent . in the case of tucana
, we derive @xmath58 values of 14.8 , 10.6 , 12.8 ( @xmath59 , @xmath60 , @xmath61 ) , corresponding to probabilities of 0.01 , 0.06 , and 0.03 .
this means that @xmath59 and @xmath61 are not consistent with a flat distribution , while @xmath60 is only marginally consistent .
although @xmath61 has larger uncertainties than the other ratios , it is interesting to note that , qualitatively , @xmath61 shows larger fluctuations than the others .
moreover , in the range @xmath62 the trend for the two galaxies is opposite .
cetus presents an initial increase followed by a drop , while in tucana an initial decrease within the core radius is followed by the increase in the ratio of bsss relative to hb stars . at large distance , there is mild evidence of an increase of @xmath61 , especially for tucana where the field covers a large area in terms of @xmath49 .
this is in agreement with the predictions by @xcite .
another way to represent the radial distribution of different samples is to use the cumulative distributions , shown in figure [ fig : cumul ] for the bss , rgb , sgb , and hb stars .
a kolmogorov - smirnov test gives probabilities of 0.861 ( 0.004 ) , 0.48 ( 0.02 ) , and 0.58 ( 0.28 ) for cetus ( tucana ) that the bsss and the rgb / hb / sgb samples are drawn from the same parent population .
however , it is interesting to note that , starting from 3 @xmath49 , the bss cumulative distribution is systematically lower , in both galaxies , than the others .
this indicates higher bss frequencies at larger distances .
@xcite found that a linear correlation exists between the logarithmic frequency of bsss ( f@xmath63 , normalized to the number of hb stars ) and the absolute visual magnitude of nearby lg dwarfs .
interestingly enough , they found that this correlation is verified over a large range of * magnitudes * , and that its slope is significantly different , i.e. , less steep , than a similar relation found by @xcite for ggcs . presently , an underlying physical mechanism for this trend has not been identified
based on the selections presented in figure [ fig : samples ] , we estimated that f@xmath63 is equal to 0.05@xmath480.05 and 0.02@xmath480.04 for cetus and tucana , respectively .    in figure
[ fig : fm ] , these values are plotted together with data from other galaxies given in figure  2 of @xcite .
we adopted the absolute @xmath64 magnitude from @xcite for cetus ( @xmath65 mag ) and from @xcite for tucana ( @xmath66 mag ) .
[ fig : fm ] shows that tucana and cetus are consistent with the relationship derived from other lg dwarfs . using the analytic relation for @xmath67 from @xcite
, we obtain @xmath680.08 and 0.03@xmath480.08 for cetus and tucana , respectively , again , showing very good agreement with the observed values .
if we accept that the @xmath69 relation is significant , then this is an independent indication that the we are dealing with a genuine bss population .
if the blue plume stars were due to fluctuations in the sfh producing younger stars , one would not expect any relationship between @xmath67 and @xmath64 . on the other hand ,
if blue plume stars are bsss , they essentialy belong to the same population generating the hb stars . as already noted by @xcite ,
the trend is similar in the case of globular clusters , but with significantly shallower slope .
the existence of an anti - correlation with the luminosity ( * the higher the luminosity , the lower the bss frequency * ) , and hence with the mass of the cluster , might suggest that in more massive globular clusters the mechanisms affecting the binary systems are more efficient then in less massive ones , and this cause a deficiency of bss with respect of hb stars ( see * ? ? ? * ) . in dsph
, the shallower slope might indicate that the same mechanisms are still at work , but are less efficient due to the lower density of dwarf galaxies .    on the basis of present results
, it appears evident that if the reliability of @xmath70 relation was confirmed by increasing the number of lg dwarfs in the sample , it could represent a valuable tool for discriminating stellar systems hosting truly young ms stars from those with a significant population of old , genuine bsss .
for example , a preliminary analysis of a deep cmd of fornax dsph ( del pino et al .
2011 , in preparation , also shown in figure  [ fig : fm ] as a crossed - circle symbol ) suggests that the f@xmath63 estimate for this galaxy ( f@xmath71 ) is significantly higher than expected from its absolute magnitude ( f@xmath72 ) .
this is in agreement with the conclusion presented in @xcite , and with the fact that fornax hosts a population as young as 100  myr ( see also * ? ? ?
further clues about the origin of bsss can be inferred from their dynamical evolution ( see , e.g. , * ? ? ?
we thus perform a wide sample of dynamical simulations , following the evolution of mt - bss and col - bss populations in the potential of cetus and tucana .
we use the code bev @xcite , in the upgraded version described by @xcite .
the same code has been used to model bsss in draco , ursa minor @xcite , sculptor and fornax @xcite .
bev integrates the dynamical evolution of bsss under the influence of the host - galaxy potential , dynamical friction , distant encounters with other stars , and three - body encounters ( which are negligible for dsphs ) .
the potential of the host galaxy is assumed to be time - independent and is represented by a multi - mass king model .
the classes of mass are the same as in @xcite , and the to mass is assumed to be 0.8 @xmath73 ( see figure [ fig : track ] ) for both cetus and tucana . to calculate the potential , we input the number density of stars in the core ( @xmath74 and @xmath75 stars pc@xmath76 for cetus and tucana , respectively ) the velocity dispersion in the core ( we assume @xmath77 km s@xmath78 for both cetus and tucana ) , and we modify the central adimensional potential @xmath79 ( defined in * ? ? ? * ) , equating the simulated concentration and radial density profile with the observed values .
[ fig : simprofile ] shows the resulting models for cetus ( obtained for @xmath80 and @xmath81 ) and for tucana ( obtained for @xmath82 and @xmath83 ) . in the case of tucana ,
the simulated profile and the concentration of the model match the observations reported in @xcite . for the concentration of cetus
, we adopt @xmath84 , which best matches our data .
we note that the model for cetus has a relatively high value of @xmath79 , with respect to other dsphs ( generally , @xmath85 for globular clusters and @xmath86 for dsphs ) .    in our model , bsss
are generated with a given position , velocity , mass and lifetime .
the fundamental differences between col - bsss and mt - bsss are the following :    \(i ) the distribution of initial positions : the initial positions of mt - bsss are randomly chosen according to a probability distribution homogeneous in radial distance from the center .
this means that mt - bsss are initially distributed according to an isothermal sphere , since mt - bsss are expected to follow the same distribution as the primordial binaries ( see * ? ? ?
* ; * ? ? ?
( the spherical nature of the dsphs is not conclusively established , but is chosen as a reasonable starting point . ) the minimum and the maximum value of the distribution of initial radial distances of mt - bsss , @xmath87 and @xmath88 ( both defined as three - dimensional quantities ) , have been tuned in order to find the best - fitting simulation . on the other hand , col - bsss
are generated only inside @xmath49 ( i.e. , @xmath89 and @xmath90 in all the simulations ) , the only region were collisions might occur , although collisions are highly unlikely in dsphs ( the probability of a three - body interaction in the core of cetus and tucana is @xmath91 times lower than in a moderately dense globular cluster ) .
the initial distribution of col - bsss is not necessarily homogeneous in the radial distance , as they are not expected to follow the same distribution as binaries ( in our simulations , we assume a constant probability distribution between the center of the cluster and @xmath49 ) .
\(ii ) the distribution of initial velocities : the initial velocities of mt - bsss are drawn from the equilibrium distributions of dsph stars with the same masses as the mt - bsss .
col - bsss are assumed to be born with a natal kick ( see * ? ? ?
we adopt a kick velocity , @xmath92 , @xmath93 . for kick velocities
higher than @xmath94 most of col - bsss are ejected from the dsph , whereas for kick velocities lower than @xmath95 , col - bsss are almost indistinguishable from mt - bsss from a dynamical point of view .    in most of the runs ,
the masses of the bsss are assumed to be @xmath96 = 1.3 m@xmath97 .
we also calculated simulations with @xmath96 = 0.9 , 1.1 m@xmath97 , but the results are substantially unchanged .
each bss is evolved for a time @xmath98 , randomly selected from a homogeneous distribution between @xmath99 and @xmath100 .
the parameter @xmath101 is the lifetime of bsss ( see * ? ? ?
we calculated simulations with @xmath102 2 and 4 gyr ( 4 gyr is the age inferred from the derived sfh ) .
tables [ tab : tab02 ] and [ tab : tab03 ] show the parameters of the simulations and report also the fraction @xmath103 of col - bsss present in each run ( @xmath104 , where @xmath105 and @xmath106 are the number of col - bsss and of mt - bsss , respectively ) . the results of the @xmath58 analysis are also included . for each run ( we made 38 runs for cetus and 37 for tucana ) 10000 bsss have been simulated .
we show three values of the non - reduced @xmath58 , referring to @xmath107 ( @xmath108 , column 8 of tables [ tab : tab02 ] and [ tab : tab03 ] ) , @xmath109 ( @xmath110 , column 9 ) and @xmath111 ( @xmath112 , column 10 ) .
[ fig : dyn ] shows the radial distributions of @xmath113 , @xmath109 , and @xmath111 obtained for the best - matching simulations ( run c3 and run t7 for cetus and tucana , respectively ) , compared with the observed distributions ( the same as in fig .
[ fig : grad_rad ] ) .
we stress that only the radial distribution of @xmath114 comes from the simulations , and it is compared with the observed radial distribution of @xmath115 , @xmath116 and @xmath117 .
in fact , the adopted simulation method does not allow to evolve these specific stellar types , but only stars in a given class of mass ( e.g. , mapelli et al .
we also remind that the underlying potential of the dsph is analytic and time - independent .
the best - matching simulations reproduce well the observed radial distributions of @xmath107 , @xmath118 , and @xmath111 .
in fact , the values of the minimum non - reduced @xmath119 are @xmath120 , @xmath121 and @xmath122 for cetus ( run c3 ) and @xmath123 , @xmath124 and @xmath125 for tucana ( runs t5 , t7 , t14 , t15 , and t16 ) .
fig  [ fig : chisq ] shows the behavior of @xmath108 ( that of @xmath110 and of @xmath112 is very similar , see tables [ tab : tab02 ] and [ tab : tab03 ] ) as a function of the fraction of col - bsss ( @xmath103 ) . as it was reasonable to expect , adding a population of col - bsss significantly reduces the agreement between data and model : for both cetus and tucana , only simulations with @xmath126 give an acceptable @xmath119 .
the dependence of the results on the bsss mass ( in the range allowed by the observations ) is negligible .
this is in agreement with the findings of previous papers @xcite for bsss in draco , ursa minor , and sculptor .
the dependence on the lifetime of bsss is negligible for cetus ( where the runs with @xmath127 gyr are almost equivalent to the runs with @xmath128 gyr ) and slightly more important for tucana ( where the runs with @xmath127 gyr have generally lower values of @xmath119 with respect to the runs with @xmath128 gyr , but the difference is not statistically significant ) . given the range of masses expected for the bss , suggested by the comparison with stellar tracks , we assume 4 gyr as fiducial value .    in the case of cetus , for the best - matching runs @xmath129 and @xmath130 , indicating that primordial binaries were initially distributed ranging from the center of the dsph out to ( at least ) @xmath131 .
the tidal radius of cetus ( derived from our data ) is 24@xmath132 , i.e. , @xmath133 , much larger than @xmath134 .
however , this difference likely does not mean that bsss do not form beyond @xmath131 , but is a consequence of the fact that we do not have data beyond @xmath135 ( furthermore , the coverage of cetus is not complete even inside @xmath136 ) .
in the case of tucana , for which data are available almost up to the tidal radius , the best - matching value of @xmath137 is equal to the tidal radius .
interestingly , the best - matching values of @xmath138 for tucana are all larger than one core radius .    in summary ,
the results of dynamical simulations show that it is possible to explain the observed distribution of bss candidates , in both cetus and tucana , with a mt - bss population .
any population of col - bsss more numerous than the 10% of the total number of bsss reduces the agreement between data and simulations .
this result is not surprising , as basic dynamical calculations tend to exclude the existence of col - bsss in dsphs , but it is important to point out that the simulations can almost rule out the unphysical scenario where @xmath139% of bsss are col - bsss in both cetus and tucana .
cccccccccc c1 & 0.0 & 5.0 & 0.9 & 4 &  & 0 & 7.7 & 6.5 & 4.7 + c2 & 0.0 & 5.0 & 1.1 & 4 &  & 0 & 9.5 & 8.1 & 5.8 + * c3 * & * 0.0 * & * 5.0 * & * 1.3 * & * 4 * & *  * & * 0 * & * 5.8 * & * 4.9 * & * 3.5 * + c4 & 0.1 & 5.0 & 1.3 & 4 &  & 0 & 8.3 & 7.0 & 5.0 + c5 & 0.2 & 5.0 & 1.3 & 4 &  & 0 & 7.3 & 6.2 & 4.5 + c6 & 0.3 & 5.0 & 1.3 & 4 &  & 0 & 10 & 8.7 & 6.2 + c7 & 0.4 & 5.0 & 1.3 & 4 &  & 0 & 14 & 12 & 8.5 + c8 & 0.5 & 5.0 & 1.3 & 4 &  & 0 & 19 & 16 & 12 + c9 & 1.0 & 5.0 & 1.3 & 4 &  & 0 & 51 & 42 & 31 + c10 & 0.0 & 5.0 & 0.9 & 2 &  & 0 & 6.7 & 5.6 & 4.1 + c11 & 0.0 & 5.0 & 1.1 & 2 &  & 0 & 7.2 & 6.1 & 4.4 + c12 & 0.0 & 5.0 & 1.3 & 2 &  & 0 & 9.8 & 8.3 & 5.9 + c13 & 0.1 & 5.0 & 1.3 & 2 &  & 0 & 9.0 & 7.6 & 5.4 + c14 & 0.2 & 5.0 & 1.3 & 2 &  & 0 & 7.0 & 5.9 & 4.2 + c15 & 0.3 & 5.0 & 1.3 & 2 &  & 0 & 7.5 & 6.3 & 4.6 + c16 & 0.4 & 5.0 & 1.3 & 2 &  & 0 & 12 & 10 & 7.5 + c17 & 0.5 & 5.0 & 1.3 & 2 &  & 0 & 20 & 17 & 12 + c18 & 1.0 & 5.0 & 1.3 & 2 &  & 0 & 50 & 41 & 30 + c19 & 0.0 & 4.0 & 1.3 & 4 &  & 0 & 10 & 8.7 & 6.2 + c20 & 0.2 & 4.0 & 1.3 & 4 &  & 0 & 9.4 & 7.8 & 5.5 + c21 & 0.5 & 4.0 & 1.3 & 4 &  & 0 & 7.9 & 6.7 & 4.8 + c22 & 1.0 & 4.0 & 1.3 & 4 &  & 0 & 41 & 34 & 25 + c23 & 0.0 & 6.0 & 1.3 & 4 &  & 0 & 13 & 11 & 7.7 + c24 & 0.2 & 6.0 & 1.3 & 4 &  & 0 & 16 & 13 & 9.1 + c25 & 0.5 & 6.0 & 1.3 & 4 &  & 0 & 26 & 22 & 16 + c26 & 1.0 & 6.0 & 1.3 & 4 &  & 0 & 77 & 64 & 46 + c27 & 0.0 & 5.0 & 1.3 & 4 & 1 & 0.1 & 7.2 & 6.0 & 4.3 + c28 & 0.0 & 5.0 & 1.3 & 4 & 1 & 0.2 & 13 & 11 & 7.8 + c29 & 0.0 & 5.0 & 1.3 & 4 & 1 & 0.3 & 20 & 17 & 12 + c30 & 0.0 & 5.0 & 1.3 & 4 & 1 & 0.4 & 28 & 23 & 16 + c31 & 0.0 & 5.0 & 1.3 & 4 & 1 & 0.5 & 43 & 36 & 26 + c32 & 0.0 & 5.0 & 1.3 & 4 & 1 & 0.6 & 56 & 47 & 33 + c33 & 0.0 & 5.0 & 1.3 & 4 & 1 & 0.7 & 76 & 63 & 45 + c34 & 0.0 & 5.0 & 1.3 & 4 & 1 & 0.8 & 94 & 78 & 56 + c35 & 0.0 & 5.0 & 1.3 & 4 & 1 & 0.9 & 110 & 91 & 65 + c36 &  &  & 1.3 & 4 & 1 & 1.0 & 122 & 101 & 72 + c37 & 0.0 & 5.0 & 1.3 & 4 & 2 & 0.1 & 9.4 & 7.9 & 5.7 + c38 & 0.0 & 5.0 & 1.3 & 4 & 2 & 0.2 & 14 & 12 & 8.3 +    cccccccccc t1 & 0.0 & 5.4 & 1.3 & 4 & 0 &  & 101 & 78 & 60 + t2 & 0.5 & 5.4 & 1.3 & 4 & 0 &  & 51 & 40 & 30 + t3 & 1.0 & 5.4 & 1.3 & 4 & 0 &  & 21 & 17 & 12 + t4 & 1.5 & 5.4 & 1.3 & 4 & 0 &  & 5.7 & 4.6 & 3.4 + t5 & 1.7 & 5.4 & 0.9 & 4 & 0 &  & 2.6 & 2.1 & 1.6 + t6 & 1.7 & 5.4 & 1.1 & 4 & 0 &  & 3.2 & 2.6 & 1.9 + * t7 * & * 1.7 * & * 5.4 * & * 1.3 * & * 4 * & * 0 * & *  * & * 2.4 * & * 1.9 * & * 1.5 * + t8 & 2.0 & 5.4 & 1.3 & 4 & 0 &  & 7.7 & 6.1 & 4.8 + t9 & 2.5 & 5.4 & 1.3 & 4 & 0 &  & 18 & 14 & 11 + t10 & 0.0 & 5.4 & 1.3 & 2 & 0 &  & 106 & 83 & 63 + t11 & 0.5 & 5.4 & 1.3 & 2 & 0 &  & 55 & 43 & 33 + t12 & 1.0 & 5.4 & 1.3 & 2 & 0 &  & 23 & 18 & 14 + t13 & 1.5 & 5.4 & 1.3 & 2 & 0 &  & 2.5 & 2.1 & 1.5 + t14 & 1.7 & 5.4 & 0.9 & 2 & 0 &  & 2.3 & 1.8 & 1.4 + t15 & 1.7 & 5.4 & 1.1 & 2 & 0 &  & 1.9 & 1.5 & 1.1 + t16 & 1.7 & 5.4 & 1.3 & 2 & 0 &  & 1.6 & 1.3 & 0.95 + t17 & 2.0 & 5.4 & 1.3 & 2 & 0 &  & 4.8 & 3.7 & 2.9 + t18 & 2.5 & 5.4 & 1.3 & 2 & 0 &  & 16 & 13 & 9.8 + t19 & 0.0 & 4.0 & 1.3 & 4 & 0 &  & 139 & 109 & 83 + t20 & 0.5 & 4.0 & 1.3 & 4 & 0 &  & 83 & 65 & 49 + t21 & 1.0 & 4.0 & 1.3 & 4 & 0 &  & 41 & 32 & 24 + t22 & 1.5 & 4.0 & 1.3 & 4 & 0 &  & 15 & 12 & 8.6 + t23 & 1.7 & 4.0 & 1.3 & 4 & 0 &  & 5.5 & 4.5 & 3.3 + t24 & 2.0 & 4.0 & 1.3 & 4 & 0 &  & 6.3 & 4.9 & 3.8 + t25 & 2.5 & 4.0 & 1.3 & 4 & 0 &  & 11 & 8.7 & 6.7 + t26 & 1.7 & 5.4 & 1.3 & 4 & 1 & 0.1 & 7.8 & 6.0 & 4.5 + t27 & 1.7 & 5.4 & 1.3 & 4 & 1 & 0.2 & 24 & 18 & 14 + t28 & 1.7 & 5.4 & 1.3 & 4 & 1 & 0.3 & 45 & 34 & 26 + t29 & 1.7 & 5.4 & 1.3 & 4 & 1 & 0.4 & 98 & 76 & 58 + t30 & 1.7 & 5.4 & 1.3 & 4 & 1 & 0.5 & 140 & 109 & 84 + t31 & 1.7 & 5.4 & 1.3 & 4 & 1 & 0.6 & 199 & 156 & 120 + t32 & 1.7 & 5.4 & 1.3 & 4 & 1 & 0.7 & 247 & 195 & 149 + t33 & 1.7 & 5.4 & 1.3 & 4 & 1 & 0.8 & 297 & 236 & 180 + t34 & 1.7 & 5.4 & 1.3 & 4 & 1 & 0.9 & 330 & 262 & 200 + t35 &  &  & 1.3 & 4 & 1 & 1.0 & 377 & 301 & 229 + t36 & 1.7 & 5.4 & 1.3 & 4 & 2 & 0.1 & 5.8 & 4.5 & 3.3 + t37 & 1.7 & 5.4 & 1.3 & 4 & 2 & 0.2 & 24 & 18 & 15 +
the results discussed in the previous sections provide evidence that the two isolated dwarf galaxies cetus and tucana host a significant population of bsss , whose properties strongly suggest that they are the product of mass exchange in binary systems
. we did not find any evidence supporting the hypothesis that part of the objects populating the blue plume are intermediate - age or young ms stars . on the basis of these findings
, it appears that both cetus and tucana are similar to ursa minor , draco @xcite , and sculptor @xcite . in the following , we first focus on the comparison of tucana and sculptor , two dwarfs which show close similarities in terms of their sfh and present - day stellar content .
sculptor is a well established case of a dsph galaxy hosting two different old populations . @xcite
first noted evidence of two rgb bumps in sculptor .
@xcite presented a photometric and spectroscopic ( 235 rgb stars ) investigation of this galaxy , showing the presence of strong radial gradients and two stellar populations with different kinematic and chemical properties .
moreover , the radial distributions of the red and blue hb components are correlated with the properties of the rgb stars , suggesting that red stars are more centrally concentrated , more metal - rich , and have a higher velocity dispersion , while the blue component has a broader distribution , lower metallicity , and a smaller velocity dispersion .
@xcite showed the first evidence of stellar radial gradients in tucana ; after dividing both the hb and the rgb into red and blue components , the number of stars belonging to the red components decreases faster than the corresponding blue sequence with increasing galaxy radius .
moreover , the analysis of the properties of rr lyrae variables stars @xcite , the sfh @xcite , and the double rgb bump @xcite strongly support the coexistence of two old ( @xmath140 10 gyr ) populations in tucana , with slightly different ages , metallicities , and spatial distributions .
0.4 cm    figure [ fig : scltuc ] compares the ratio of the number of bsss versus red ( asterisks ) and blue ( squares ) hb stars as a function of galactocentric radius for the two galaxies .
the values for sculptor are taken from table 2 of @xcite .
interestingly , in both cases the bsss are more centrally concentrated than the blue component of the hb , and less concentrated than the red component .
we note also that , for each galaxy , the two curves cross at a similar galactocentric distance , @xmath141 @xmath49 . in the case of tucana , this trend is confirmed also for the bsss and the hb stars when compared to red and blue rgb stars ( figure [ fig : grad_rad2 ] , top left and top right , respectively ) .
therefore , if we compare the global n@xmath142/n@xmath143 ratio in figure [ fig : grad_rad ] with figure [ fig : scltuc ] and figure [ fig : grad_rad2 ] ( top panels ) , it can be inferred that the radial trends of the bsss are indeed driven by the fast decrease of the red hb stars in the central regions , associated with the younger , more metal - rich component , while the increase at larger radii is driven by the more homogeneously distributed bluest stars .
this means that the peculiar radial trend of the candidate bsss when they are normalized to the hb stars is not due to a peculiar distribution of these objects inside the galaxy but only to the existence of a significant metallicity gradient in tucana , a gradient that is well traced by both the red and blue hb stellar populations .
the similarities between these two galaxies are particularly intriguing because of the remarkably different environmental conditions : tucana spent most of its life in isolation , while sculptor is a close satellite of the milky way .
if all the bsss are descendants of primordial binaries , this suggests that the different environments of the two galaxies did not alter their binary distributions .
in contrast , figure [ fig : grad_rad2 ] ( bottom panels ) shows that the properties of different stellar samples in cetus appear more homogeneous .
this would be in agreement with other findings , as the properties of the rr lyrae variable stars and of the rgb bump , which indicate a single dominant population characterized by a more uniform distribution of ages and metallicities .
interesting features appear when comparing the distribution of the bsss with those of the red and blue stars of both the rgb and hb .
we defined these samples as follows .
the region selecting the rgb stars was simply split in two , similarly to what has been done for the analysis of the rgb bump in @xcite .
the red and blue hb are defined considering only the stars redder than @xmath144 and bluer than @xmath145 in the hb region , therefore mostly avoiding the rr lyrae instability strip .
figure [ fig : cumul2 ] compares the distributions of these four samples with the bsss for both tucana and cetus . in the case of cetus , we find that the bsss are significantly less concentrated than the red hb and more concentrated than the blue hb , while they more closely follow the distributions of both of the rgb components . in the case of tucana ,
the bsss closely follow the distributions of both blue components ( possibly more concentrated in the central regions ) , while they are significantly less concentrated than the red sequences .
we have demonstrated that the populations of bss candidates in tucana and cetus are consistent with all predictions of creation from mt - bsss .
in fact , it would be unprecedented to have a truncation in star formation without an associated population of mt - bsss .
unfortunately , our limited knowledge of the evolution of binary stellar systems prevents us from making accurate predictions of the expected numbers of mt - bsss . in sum
, we expect to see mt - bsss in dsph galaxies , and we detect them with consistent properties .
the remaining question is whether some fraction of the bss candidates could still be due to newly formed ms stars .
the inconsistency with the age - metallicity relationship could be due to infall of a gas cloud with a low metallicity coincident with that of the main population which then produces subsequent star formation .
however , the truncation at luminosities corresponding to twice the to mass implies a second unlikely coincidence if some of the bss candidates were true ms stars . while very difficult to rule out such a scenario completely , it is clear that newly formed ms stars are not necessary to explain all of the observations .
in this paper , we have presented the analysis of candidate bsss in two isolated dsphs in the periphery of the local group , cetus and tucana .
deep hst / acs data allowed us to identify 940 and 1214 bss candidates in cetus and tucana , respectively .
different indicators suggest that these are true bsss ; and , thus , there are no young main sequence stars in either galaxy .
the first evidence comes from the sfh we derived for both objects @xcite . not allowing for the presence of bsss , then a small population of relatively young ( 3 to 5 gyr old ) stars is detected , with metallicities similar to the oldest most metal - poor stars , and thus , much lower than expected from the age - metallicity relation of the dominant population .
the comparison with stellar tracks presented in this paper suggests that the bsss have typical masses @xmath146 , and smaller than @xmath147 , in agreement with the fact that the maximum expected mass of a bss is twice the mass presently evolving at the to ( @xmath148 ) .
moreover , the limited spread in metallicity found for the bss populations , perfectly in agreement with that of turn - off and subgiant stars , suggests that the use of stellar evolution models for single stars does not introduce strong biases in the analysis .
the analysis of the luminosity functions and radial distributions support the bss interpretation . on one hand
, we showed that the luminosity function does not change as a function of radius , and in particular the brightest ( youngest ) stars do not appear more centrally concentrated , as one would expect from a young main sequence . on the other hand
, the radial profiles do not show any central peak , as observed in globular clusters due to col - bss stars or as it would be in the case of a young component .
thus , the present analysis also suggests that the bsss in cetus and tucana formed from primordial binaries and , as in the case of other dsphs , col - bsss are unlikely to form .
this conclusion is strongly supported by the dynamical simulations we performed , taking into account the properties of the dwarfs under scrutiny .
these simulations indicate that the observed bss distributions are well reproduced when only mt - bsss are taken into account , while the fit worsens when even a small number of col - bsss ( @xmath149% ) is included in the dynamical simulation .
we can therefore safely conclude that it is highly likely that cetus and tucana host a genuine population of bss stars , and we confirm the conclusions presented in previous studies that both galaxies did not experience any star formation episodes in the last 8 gyr .
it is particularly interesting that the positions of cetus and tucana in the f@xmath63 vs.  @xmath64 plane fit very well in the relation presented by @xcite .
this , in turn , could be a powerful instrument to identify young main sequence stars in resolved spheroidal galaxies , discriminating from systems with only old populations from those also hosting younger stars .
the authors thank the anonymous referee for the useful comments , yazan momany for sharing the data used in figure  [ fig : fm ] , and giuliana fiorentino for useful discussions on the acs .
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 , the iac ( grant 310394 ) , the education and science ministry of spain ( grants aya2004 - 06343 , aya2007 - 3e3507 , and aya2010 - 16717 ) .
s.c . acknowledges the financial support of inaf through the prin inaf 2009 ( p.i . :
r. gratton ) .
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 .