more than seven decades ago , @xcite discovered that globular clusters ( gcs ) in the milky way are divided into two distinct groups according to the mean period ( @xmath1 ) of type ab rr lyrae variables .
the oosterhoff group i gcs are relatively metal - rich , and have @xmath1 of @xmath00.55 day , while more metal - poor group ii gcs are observed to have longer @xmath1 of @xmath00.65 day ( see figure  [ fig1 ] ) . despite the importance of this phenomenon in our understanding of the milky way formation and population ii distance scale ( * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* and references therein ) , this still remains one of the long - standing problems in modern astronomy .
numerous investigations were carried out to find the origin of the oosterhoff dichotomy since its discovery . among them , several deserve more detailed description .
@xcite suggested  hysteresis mechanism " , which explains the dichotomy as a difference in mean temperature between the type ab variables in the two groups .
sandage ( 1981 ) , however , discovered the period - shift , the difference in period at given temperature , between the rr lyrae stars in the two groups , which led him to conclude that the dichotomy is caused by the luminosity difference . by employing the synthetic horizontal - branch ( hb ) models , @xcite found that rr lyrae variables highly evolved from the zero - age horizontal - branch ( zahb ) can reproduce the period - shift for the group ii gcs if hb morphology is sufficiently blue . while this evolution effect can play an important role in the oosterhoff dichotomy ( see * ? ? ?
* ) , it is not clear whether this effect alone could reproduce the observed dichotomy , because hb types of some group ii gcs are too red to have enough evolution effect .
all of these investigations were performed on an assumption that gcs host only single population ( see also * ? ? ?
recent observations and population models , however , have convincingly demonstrated the presence of multiple populations in most gcs ( * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* and references therein ) .
this has motivated @xcite to suggest that the origin of the oosterhoff dichotomy is intimately related to the presence of multiple populations in gcs . according to this new paradigm ,
the instability strip ( is ) in the metal - poor group ii gcs is populated by second - generation stars ( g2 ) with enhanced helium and cno abundances , while the rr lyrae stars in the relatively metal - rich group i gcs are produced mostly by first - generation stars ( g1 ) without these enhancements .
this  population - shift " within the is with increasing metallicity can create the observed period - shift between the two groups .
synthetic hb models presented by @xcite were focused on the two outer halo gcs representing two oosterhoff groups , m15 and m3 .
it was , therefore , not clear whether this population - shift effect applies to all gcs in addition to these two gcs , and also whether there is any difference between the inner and outer halo gcs .
the purpose of this paper is to extend our multiple population models to all metallicity regimes to see whether the population - shift effect can reproduce the observed period dichotomy , and to suggest , based on @xmath1 versus [ fe / h ] correlations , new constraints on the star formation time scales for the inner and outer halo gcs .
table  [ tab1 ] lists milky way gcs containing 3 or more type ab rr lyrae variables .
metallicity ( [ fe / h ] ) and galactocentric distance ( @xmath2 , in kpc ) are from 2010 edition of @xcite catalogue .
data for the mean periods ( in days ) and numbers of rr lyrae stars ( @xmath1 , @xmath3 , @xmath4 , and @xmath5 ) are from the sources listed in the table , where double mode rr lyrae variables , if any , are treated as c type variables ( see * ? ? ?
for the gcs without the source listed in the table , they are from the compilations by @xcite and @xcite .
updated data for the hb type @xcite are mostly from h. gim & y .- w . lee ( 2009 , private communication ) and from our own compilation .
figure  [ fig1 ] shows the oosterhoff - arp diagram based on the data in table  [ tab1 ] .
as is clear from this figure , oosterhoff group ii gcs are relatively metal - poor , while group i gcs are more metal - rich .
however , two most metal - rich gcs ( ngc 6388 and ngc 6441 ) have longest @xmath1 and are classified as group iii @xcite .
the red open circles are 2mass - gc02 and terzan 10 , two bulge gcs recently reported by @xcite .
note that there is an abrupt increase in @xmath1 at [ fe / h ] @xmath6 -1.5 dex , and that the inner halo gcs are in the mean more metal - rich than those in the outer halo .    in order to reproduce the observed period dichotomy ,
we have constructed a series of synthetic hb models with rr lyrae variables , following the technics developed by @xcite , and as updated by @xcite and @xcite to include the effects of multiple populations .
our models are based on the most updated yonsei - yale ( @xmath7 ) isochrones and hb evolutionary tracks with enhanced helium and cno abundances @xcite , all constructed under the assumption that [ @xmath8/fe ] = 0.3 dex .
these isochrones and evolutionary tracks are available from the corresponding author upon request .
@xcite mass - loss parameter @xmath9 was adopted to be 0.50 , which best reproduced observed hb morphology versus [ fe / h ] relation of the inner halo gcs when the mean age was assumed to be 13 gyr ( see * ? ? ? * for a similar result ) .
depends on the adopted absolute age for the gcs .
if we had adopted @xmath9 = 0.42 @xcite , the mean age for the inner ( old ) halo gcs required to fit the observed hb morphologies would be @xmath0 14 gyr . ] following @xcite , mass dispersion on the hb was adopted to be @xmath10 = 0.010 @xmath11 for each subpopulation .
for the rr lyrae variables , the first overtone blue edge of the instability strip was taken from @xcite with corrections for convection @xcite , which agrees , to within the uncertainty in the reddening ( @xmath12(b - v ) @xmath6 0.02 ) , with more recent calculations by @xcite .
as for the transition temperature between ab and c type variables , @xcite showed that the transition occurs close to the first overtone red edge in m15 ( group ii ) , whereas it occurs near the fundamental blue edge in m3 ( group i ) , which suggest that the hysteresis mechanism @xcite may play an important role in the  either - or " zone of the is . even if this theoretical scenario is able to reproduce many observational behaviors @xcite , some residual uncertainties on the first overtone red edge due to the significant dependence on the convection treatment may remain .
an alternative approximated approach to this problem , adopted in our calculation , is to use the @xcite s semi - empirical equation ( his eq .
17 ) for the blue edge of type ab variables as a transition temperature between the two types of rr lyrae stars .
it appears that some hysteresis effect is empirically reflected in his equation , as we obtain cooler temperatures ( @xmath13 0.008 ) for more metal - poor group ii gcs from his equation .
the width of the instability strip was taken to be @xmath14  =   0.060 @xcite , which best reproduces the observed value of @xmath1 in our models , and is acceptable to within the uncertainty of the observed width in m3 and m15 ( see figures 2 & 3 of * ? ? ?
the absolute value of @xmath1 would be increased by @xmath0 0.02 day if we had adopted @xmath14
=  0.070 , but this has only negligible effect on the relative difference in @xmath1 between the two oosterhoff groups .
the periods of fundamental mode rr lyrae stars @xmath15 were calculated from @xcite and those of first overtone mode variables @xmath16 were obtained from the equation by @xcite . in order to minimize the stochastic effect stemming from a finite number of stars
, we obtained mean values of 1000 simulations , each with the sample size of 10000 hb stars , at given parameter combination .
the major assumption made in @xcite was the presence of g2 , which is mildly enhanced in helium and cno abundances .
in terms of helium enhancement , this subpopulation is in between the g1 ( with normal helium abundance ) and super - he - rich third generation stars ( g3 , if any ) . both theories and observations suggest that these g2 and g3 would be enhanced in both helium and cno abundances from the chemical pollution and enrichment by intermediate - mass asymptotic - giant - branch stars and/or fast - rotating massive - stars @xcite . in general , helium
enhanced subpopulations , like g3 in our models , are placed on the bluer hb @xcite . however , when the helium enhancement is relatively small ( @xmath12@xmath17 @xmath6 0.01 ) , as is required for the rr lyrae variables in m15 @xcite , g2 can be placed on the relatively redder hb , because the age and cno effects can overwhelm the helium effect .
this is illustrated in figure  [ fig2 ] , where the effects of age , [ cno / fe ] , and helium abundance on the hb morphology are presented using our synthetic hb models .
first , the hb morphology is getting redder as age decreases because of the higher mean mass of hb stars at younger age .
similarly , the morphology of hb gets redder as cno abundance increases , since the efficiency of hydrogen shell burning and opacity increase with increasing [ cno / fe ] .
however , as helium abundance increases , the hb morphology gets bluer , mainly because helium - rich stars evolve more rapidly , and therefore , have lower masses at the tip of the red giant - branch ( rgb ) and hb for a given age .
the bottom panel shows the net effect of these three parameters operating simultaneously .
this confirms that the hb morphology of g2 can be redder than that of g1 when the helium enhancement is modest ( see also figure  [ fig7 ] below ) .
the effects of cno and helium abundances on the periods of rr lyrae stars are presented in figure 3 of @xcite .
in this section , we extend our models to all metallicity regimes , covering the inner and outer halo gcs , to see whether the population - shift effect can reproduce the observed period dichotomy .
our previous modeling @xcite shows that the relative placement of subpopulations ( g1 , g2 and g3 ) on the hb is qualitatively similar in most gcs , but parameter combinations ( age , @xmath17 , and [ cno / fe ] ) defining these subpopulations can be somewhat different from one gc to another .
therefore , here , it is our goal to find the representative parameter combinations for the inner and outer halo gcs respectively , which can best reproduce the observed @xmath1 versus [ fe / h ] correlations in those regimes , while maintaining reasonable agreements in the observed hb type versus [ fe / h ] relations .
first , we constructed , with a parameter combination found by @xcite from an outer halo gc m15 , a series of synthetic hb models by varying [ fe / h ] from -2.3 to -0.8 dex .
preliminary comparison of this model sequence with the outer halo gcs on the @xmath1 versus [ fe / h ] plane ( see figure  [ fig3 ] ) indicated that the presence of two distinct @xmath1 groups can be reproduced qualitatively , but some fine tuning of input parameters was required for the best - fit with the data .
a solid line in figure  [ fig3 ] is our best - fit models for the outer halo , and table  [ tab2 ] lists the best - fit input parameters obtained from this endeavour ( see figure  [ fig6 ] below for the corresponding models in the hb type versus [ fe / h ] plane ) . in our modeling , we provisionally adopted [ @xmath8/fe ] = 0.0 dex , because otherwise the sudden change in @xmath1 occurred at [ fe / h ] that is somewhat metal - poor compared to the observed trend , and the model locus was not well matched with the data at the metal - rich regime . of course
, this is not consistent with the observed [ @xmath8/fe ] = 0.2 - 0.3 dex @xcite , and therefore is equivalent to the shift of [ fe / h ] by 0.14 - 0.2 dex , which , however , is within the current uncertainty in the metallicity scale .
this shift is required because our @xmath8 , cno , and he enhanced hb tracks have excessively extended blueward loops at given [ fe / h ] .
we suspect that this is due to the implicit assumption of fully efficient semi - convection adopted in our hb stellar models , which was introduced in 70s to match the observed color spread on the hb and the primordial helium abundance ( see , e.g. , * ? ? ?
* ) , well before the discovery of he enhanced multiple populations in gcs . ]
note that @xmath12@xmath17 and @xmath18 between g1 and g2 mostly determine the values of @xmath1 for the two oosterhoff groups ( see figure 3 of * ? ? ?
* ) , while @xmath12t(g1 - g2 ) mainly controls the overall feature of the distribution in @xmath1 versus [ fe / h ] diagram .
our best - fit models were constructed under the assumption that the age difference between g1 and g2 , @xmath12t(g1 - g2 ) , is 1.4 gyr .
we found that this age difference provides the best match with the observed gcs in the outer halo , both for the groups i and ii . when the age difference is different from this value , for example , @xmath12t(g1 - g2 ) = 0.5 gyr as shown by the dashed line , the models fail to reproduce the observed distribution , which illustrates the sensitivity of our models to @xmath12t(g1 - g2 ) .
similarly , if the separation between g1 and g2 on the hb was obtained by increasing @xmath18 instead of @xmath12t(g1 - g2 ) , we would obtain @xmath12@xmath1 that is too large ( by @xmath0 0.06 day ) compared to the observations .
our best fit models ( the solid line ) show a  wavy " feature in the group ii regime ( -2.3 @xmath19 [ fe / h ] @xmath19 -1.8 dex ) , producing gcs having @xmath1 of @xmath00.65 day along this locus .
inspection of synthetic hb models indicates that this feature is dictated by the placement of g2 within the is and the contribution of evolved stars from either side of the is .
we have then applied this parameter combination determined from the outer halo to the inner halo gcs to see whether this can be similarly employed to all gcs .
the hb types for the inner halo gcs are systematically bluer than those for the outer halo gcs ( see figure  [ fig6 ] below ) , and the general consensus is that this global second parameter phenomenon stems from a small age difference between the inner and outer halo gcs @xcite .
therefore , in order to reproduce the hb types of the inner halo gcs , first , we have adopted @xmath12t(g1 - g2 ) to be 1.4 gyr as was found from the outer halo gcs , but constructed models under the condition that the inner halo gcs are @xmath01 gyr older than those in the outer halo .
other parameters , @xmath12@xmath17 , @xmath12@xmath20 , and mass - loss parameter @xmath9 , are held identical to the values adopted for the outer halo . in this case , however , we fail to reproduce the observed distribution of @xmath1 among the inner halo gcs ( see the dashed line in figure  [ fig4 ] ) . compared to the observations , the jump in @xmath1 occurs at more metal - rich regime ( [ fe / h ] @xmath6 -1 dex ) , and the  plateau " of @xmath1 distribution is formed at too low value of @xmath1 .
therefore , by constructing models under different assumptions on the value of @xmath12t(g1 - g2 ) , we found that models can best match the observations when @xmath12t(g1 - g2 ) = 0.5 gyr ( see the solid line in figure  [ fig4 ] ) .
the observed plateau of @xmath1 in -1.8 @xmath19 [ fe / h ] @xmath19 -1.6 dex is also naturally reproduced in this case , which is explained in detail in figures  [ fig9 ] and [ fig10 ] below .
our best - fit input parameters for the inner halo gcs are listed and compared with those for the outer halo in table  [ tab2 ] , and figure  [ fig5 ] compares our best - fit models with the observations both for the inner and outer halo gcs .
our models were constructed to reproduce not only the observed mean periods of rr lyrae stars , but also the hb types of halo gcs , as presented in figure  [ fig6 ] and table  [ tab3 ] .
for example , our best - fit models in figure  [ fig5 ] were constructed along the solid lines ( g1+g2 models ) in figure  [ fig6 ] . unlike @xmath1 ,
the hb type is sensitively affected by the presence of g3 in all metallicity regimes , and therefore , models including g3 are also presented in figure  [ fig6 ] ( see the dashed lines ) . as is clear from the figure , most of the observed gcs are reasonably well reproduced by our  g1+g2 " or  g1+g2+g3 " models .
note that some of the outer halo gcs show inner halo (  old halo " ) characteristics in terms of hb morphology . in order to illustrate how g1 , g2 and g3 would affect the final morphology of hb , we present in figure  [ fig7 ] the synthetic hb models at three different metallicities along the model lines in figure  [ fig6 ] . in the outer halo , where the separation between g1 and g2 is relatively large ,
gcs undergo the population - shift more clearly as a subpopulation occupying the is is progressively changing with metallicity ( see also figure 1 of * ? ? ?
* ) . in the inner halo gcs ,
however , we can see that g1 and g2 are more overlapped in the is because the separation between g1 and g2 is small . for the gcs with g3 ,
hb stars belonging to this additional subpopulation are represented as filled gray circles in figure  [ fig7 ] by employing the parameters listed in table  [ tab2 ] .
note , however , that helium abundances of g3 in real gcs would be somewhat different from one gc to another .    as described above , and shown in figure  [ fig8 ] ,
the most obvious features in our models for the inner halo gcs are ( 1 ) the sudden jump of @xmath1 at [ fe / h ] @xmath6 -1.4 dex , ( 2 ) the plateau of @xmath1 distribution in the metallicity range of -1.8 @xmath19 [ fe / h ] @xmath19 -1.6 dex , and ( 3 ) the second jump of @xmath1 at [ fe / h ] @xmath6 -1.8 dex .
we believe that this  two - stage jump " in @xmath1 is the most important effect we found for the complete understanding of the oosterhoff dichotomy in the inner halo gcs .
as compared in figure  [ fig8 ] , the first jump in @xmath1 at [ fe / h ] @xmath6 -1.4 dex ( as first discovered by * ? ? ?
* ; * ? ? ?
* ) is well reproduced by both
g1 only " and  g1+g2 " models , while the plateau in oosterhoff ii regime is only reproduced in the multiple population paradigm ( g1+g2 model ) . in order to explain this effect in more detail , we illustrate , in figure  [ fig9 ] , four synthetic hb models with evolutionary tracks in the relevant metallicity range of -2.0 @xmath19 [ fe / h ] @xmath19 -1.4 dex . at [ fe / h ]
@xmath21 -1.4 dex , type ab rr lyrae zone is mostly occupied by hb stars on or near zahb from both g1 and g2 , producing low values of mean luminosity and @xmath1 . as metallicity decreases ,
hb morphology is progressively getting bluer , and the zahb portion of g1 is crossing the fundamental blue edge ( fbe ) , leaving the zone mostly occupied by he and cno enhanced g2 , which increases @xmath1 to @xmath00.65 day ( see the second panel ) .
this state continues with little change in @xmath1 until [ fe / h ] @xmath6 -1.8 dex .
with a further decrease in metallicity , zahb portion of g2 is also getting out of the zone , and therefore only highly evolved stars from the blue side of the fbe penetrate back into the zone , increasing again mean luminosity and @xmath1 .
figure  [ fig10 ] shows a schematic diagram that explains these sequential departures of g1 and g2 zahbs toward higher @xmath22 with decreasing metallicity , which lead to the two - stage jump in @xmath1 .
note that this effect is another manifestation of the population - shift effect @xcite in the inner halo , where , unlike the outer halo , the separation between g1 and g2 is small on the hb . note further that our models for the inner halo gcs correctly predict that the population shift , and therefore the jump in @xmath1 , occurs at bluer hb type compared to the outer halo gcs ( see also figures  [ fig6 ] and [ fig7 ] ) .
as suggested by previous investigators @xcite , for the case of more metal - rich oosterhoff group iii gcs such as ngc 6441 and 6388 , the hb stars can not be populated within the is unless helium abundance is significantly enhanced .
therefore , we present our models including more helium - rich g3 in figure  [ fig11 ] , and the best - fit parameters adopted for g3 are also listed in table  [ tab2 ] .
our models with g3 are qualitatively similar to those in @xcite , but here we include the spread in helium abundance in order to reproduce the observations more realistically .
as metallicity increases , g1 and g2 cross the red edge of the is , mostly settling down in the red hb , while helium - rich g3 penetrate into the is from blue hb , producing rr lyrae stars with longer periods . as shown in figure
[ fig11 ] , the presence of g3 starts to have influences on @xmath1 from [ fe / h ] @xmath6 -1.5 dex for the outer halo , while they have effects from [ fe / h ] @xmath6 -1 dex for the inner halo . with a further increase in metallicity , more and more helium enhanced g3 are progressively penetrating into the is , continuously increasing @xmath1 .
as is evident from figure  [ fig11 ] , group iii gcs and some group i gcs in the outer halo are better reproduced by our models with g3 .
note , however , that our models with g3 were calibrated with only three gcs ( ngc 6441 , 6388 , and terzan 10 ) for which type ab variables are predicted to be fully produced by g3 , and therefore , parameters for g3 listed in table  [ tab2 ] are more uncertain compared to those for g1 and g2 .
we have shown that a possible explanation of the observed oosterhoff dichotomy can be naturally obtained when the multiple population models of @xcite are extended to all metallicity regimes covering the inner and outer halo gcs .
we found that the population - shift responsible for the oosterhoff dichotomy can be traced back to the two or three discrete episodes of star formation in gcs .
interestingly , in order to reproduce the observed period dichotomies for the inner and outer halo gcs , respectively , different star formation timescales are required .
the best - fit models suggest that the star formation in the inner halo gcs commenced and completed earlier within the timescale of @xmath00.5 gyr , while that for the outer halo gcs was delayed by @xmath00.8 gyr with more extended timescale ( @xmath01.4 gyr ) between g1 and g2 . despite the difference in detail ,
our models show that the oosterhoff period groups observed in the inner and outer halo gcs are all manifestations of the population - shift effect within the instability strip .    further supporting evidence for our models
is provided by the observed difference in mean periods ( @xmath1 - @xmath3 ) between type ab and type c variables ( see figure  [ fig12 ] ) .
interestingly , the oosterhoff group ii gcs are divided into two subgroups in this diagram , which is not obvious in @xmath1 ( see figure  [ fig11 ] ) .
following @xcite , we call these subgroups  oo iia " for gcs with [ fe / h ] @xmath21 -1.9 dex and  oo iib " for those with [ fe / h ] @xmath19 -1.9 dex
. our models can explain this rather clearly .
the c type variables in oo iia subgroup are populated by both g1 and g2 , and thus have relatively low @xmath3 , while those in oo iib are mostly composed of g2 which leads to relatively higher value of @xmath3 and smaller difference in @xmath1 - @xmath3 .
in addition to @xmath1 , the fraction of c type variables ( @xmath23 ) is , in general , also different between the two oosterhoff groups .
although our models predict @xmath24 0.4 - 0.5 for the group ii gcs and @xmath25 0.2 for more metal - rich group i gcs , which apparently agree with the observations , this result is subject to the uncertainties in the adopted blue edges , transition temperature between the two types of variables , and the width of the is .
inspection of our models indicates that the difference in @xmath23 between the two groups is mostly due to the uneven distribution in @xmath22 of rr lyrae stars within the is , which is qualitatively consistent with the conclusion from single population models by @xcite .
our models predict that there is an age difference of @xmath01 gyr on average between the inner and outer halo gcs , excluding some outer halo gcs with predominantly blue hb type which would be as old as the inner halo gcs .
this is consistent with the main - sequence age dating based on the recent hst data @xcite , which suggests an age difference of @xmath00.8 gyr between the corresponding samples of gcs .
note , however , that ages obtained from such studies do not yet explicitly account for the effects of multiple populations in gcs , and therefore , variations in he and cno abundances among subpopulations should be taken into account for better age determinations of gcs .
direct age dating of g1 and g2 in a given gc would be even more challenging , because , as illustrated in figure 3 of @xcite , the effect of younger age for g2 would be cancelled out by the effect of enhanced cno on the isochrones .
considering this , perhaps the most important result predicted by our models from the oosterhoff dichotomy is that the star formation timescale between g1 and g2 is more extended ( @xmath01.4 gyr ) in the outer halo gcs compared to that of the inner halo ( @xmath00.5 gyr ) .
this result for the outer halo gcs is qualitatively consistent with a more prolonged star formation history observed in the local group dwarf galaxies ( see , e.g. , * ? ? ?
this is most likely because of the suggested difference in the formation history , in which the inner halo was formed mostly by dissipational merger and collapse , while dissipationless accretion of stars from primordial building blocks and/or satellite galaxies were more important in the outer halo @xcite . in this regard ,
it is important to note that @xcite found a planar alignment of oo iib gcs , which suggests a captured origin from a satellite galaxy .
@xcite went further from this and suggested that most of the outer halo gcs would belong to the vast polar structure of satellite galaxies .    according to our results ,
the type ab variables in the oosterhoff group ii gcs are produced mostly by g2 , while those in group i gcs are more populated by g1 .
the fact that the rr lyrae variables with oosterhoff i and ii characteristics are similarly present in the halo field @xcite would therefore suggest that building blocks similar to gcs ( in terms of stellar populations ) provided the halo field with rr lyrae stars originated from both g1 and g2 .
unlike the super - helium - rich g3 that require very extreme formation conditions in the central region of a proto - gc ( see , e.g. , * ? ? ?
* ) , the g2 in our models experienced just a modest helium enhancement ( @xmath12@xmath17 @xmath6 0.01 ) , and therefore , they could have been formed ubiquitously in a proto - gc .
thus , it appears that both g1 and g2 were effectively provided to the halo field when merging and disruption of building blocks were much more active in the early phase of the milky way formation .
there are several remaining issues to be clarified related to this work .
for example , rr lyrae variables in the local group dwarf galaxies often show oosterhoff intermediate characteristics with @xmath1 of @xmath00.60 day ( see , e.g. , * ? ? ? * ) .
although figures  [ fig3 ] and [ fig4 ] indicate that some of our models ( see the dashed lines ) , with different assumptions on the ages of g1 and g2 , can be placed within the oosterhoff gap , further investigations are required to fully understand this phenomenon . as to the observational test , as emphasized by @xcite , spectroscopic observations for the na , he , and cno abundances can provide important tests on our placements of g1 , g2 , and g3 on the hb .
although the available observations @xcite are not in conflict with our models ( see * ? ? ?
* ) , especially because now our models ( see figure  [ fig7 ] ) show more flexibility on the placements of g1 and g2 depending on @xmath12t(g1 - g2 ) , spectroscopy for a large sample of rr lyrae and hb stars would certainly help to better understand the population - shift scenario suggested in our models .
if confirmed by further observations and models , our understanding for this long standing problem , the oosterhoff dichotomy , is now sufficiently well - advanced that the rr lyrae stars can be used as more reliable standard candle and galaxy formation tracers .
we thank the referee for a number of helpful suggestions .
we also thank hansung gim for his early contribution in hb morphology classification , and seok - joo joo for his fortran code , based on which the present code was developed to include multiple populations and the details of rr lyrae stars .
support for this work was provided by the national research foundation of korea to the center for galaxy evolution research .
carretta , e. , bragaglia , a. , gratton , r.  g. , et al .
2009 , a&a , 505 , 117 catelan , m. 2004 , in asp conf .
310 , variable stars in the local group , ed .
d. kurtz & k. pollard ( san francisco : asp ) , 113 catelan , m. 2009 , ap&ss , 320 , 261 clement , c.  m. , muzzin , a. , dufton , q. et al .
2001 , aj , 122 , 2587 cohen , j.  g. , briley , m.  m. , & stetson , p.  b. 2005 , aj , 130 , 1177                                                                                piotto , g. , bedin , l.  r. , anderson , j. , et al .
2007 , apjl , 661 , l53 pritzl , b.  j. , smith , h.  a. , stetson , p.  b. , et al .
2003 , aj , 126 , 1381 pritzl , b.  j. , venn , k.  a. , irwin , m. , et al .
2005 , aj , 130 , 2140                                                            ngc 362 & & -1.26 & 0.585 & 24 & 0.365 & 3 & -0.82 & 9.4 & 1 + ngc 1261 & & -1.27 & 0.555 & 13 & 0.328 & 5 & -0.71 & 18.1 & + ngc 1851 & & -1.18 & 0.571 & 21 & 0.311 & 8 & -0.31 & 16.6 & + ngc 2419 & & -2.15 & 0.662 & 38 & 0.367 & 37 & 0.76 & 89.9 & 2 + ngc 2808 & & -1.14 & 0.563 & 11 & 0.300 & 5 & 0.09 & 11.1 & 3 + ngc 3201 & & -1.59 & 0.554 & 72 & 0.312 & 5 & 0.08 & 8.8 & + ngc 4147 & & -1.80 & 0.525 & 5 & 0.339 & 10 & 0.55 & 21.4 & 4 + ngc 4590 & m68 & -2.23 & 0.622 & 14 & 0.376 & 32 & 0.44 & 10.2 & 5 + ngc 4833 & & -1.85 & 0.708 & 7 & 0.380 & 7 & 0.95 & 7.0 + ngc 5024 & m53 & -2.10 & 0.649 & 29 & 0.348 & 29 & 0.81 & 18.4 + ngc 5053 & & -2.27 & 0.672 & 5 & 0.354 & 4 & 0.54 & 17.8 + ngc 5139 & @xmath26 cen & -1.53 & 0.651 & 76 & 0.362 & 85 & 0.90 & 6.4 + ngc 5272 & m3 & -1.50 & 0.561 & 158 & 0.336 & 48 & 0.08 & 12.0 & 6 + ngc 5286 & & -1.69 & 0.656 & 30 & 0.333 & 22 & 0.74 & 8.9 & 7 + ngc 5466 & & -1.98 & 0.646 & 13 & 0.348 & 7 & 0.59 & 16.3 + ngc 5634 & & -1.88 & 0.660 & 12 & 0.339 & 7 & 0.94 & 21.2 & 8 + ngc 5824 & & -1.91 & 0.624 & 7 & & 0 & 0.91 & 25.9 + ngc 5904 & m5 & -1.29 & 0.551 & 91 & 0.318 & 35 & 0.39 & 6.2 + ngc 5986 & & -1.59 & 0.652 & 7 & 0.328 & 1 & 0.94 & 4.8 + ngc 6093 & m80 & -1.75 & 0.680 & 7 & 0.368 & 7 & 0.93 & 3.8 & 9 + ngc 6121 & m4 & -1.16 & 0.553 & 31 & 0.294 & 13 & -0.07 & 5.9 & 10 + ngc 6139 & & -1.65 & 0.662 & 3 & 0.417 & 1 & 0.97 & 3.6 + ngc 6171 & m107 & -1.02 & 0.538 & 15 & 0.287 & 7 & -0.76 & 3.3 + ngc 6229 & & -1.47 & 0.553 & 30 & 0.322 & 8 & 0.36 & 29.8 + ngc 6266 & m62 & -1.18 & 0.547 & 133 & 0.302 & 76 & 0.60 & 1.7 & 11 + ngc 6284 & & -1.26 & 0.588 & 6 & & 0 & 0.90 & 7.5 + ngc 6333 & m9 & -1.77 & 0.639 & 8 & 0.336 & 10 & 0.85 & 1.7 & 12 + ngc 6341 & m92 & -2.31 & 0.630 & 11 & 0.356 & 6 & 0.87 & 9.6 + ngc 6362 & & -0.99 & 0.547 & 18 & 0.293 & 17 & -0.45 & 5.1 + ngc 6388 & & -0.55 & 0.676 & 9 & 0.387 & 11 & -0.75 & 3.1 & 13 +
ngc 6402 & m14 & -1.28 & 0.564 & 39 & 0.330 & 15 & 0.82 & 4.0 + ngc 6426 & & -2.15 & 0.704 & 9 & 0.332 & 4 & 0.53 & 14.4 + ngc 6441 & & -0.46 & 0.759 & 42 & 0.404 & 25 & -0.81 & 3.9 & 13 + ngc 6558 & & -1.32 & 0.556 & 6 & 0.345 & 3 & 0.84 & 1.0 + ngc 6584 & & -1.50 & 0.560 & 34 & 0.304 & 8 & -0.09 & 7.0 + ngc 6626 & m28 & -1.32 & 0.577 & 8 & 0.312 & 2 & 0.96 & 2.7 + ngc 6642 & & -1.26 & 0.544 & 10 & 0.322 & 6 & 0.43 & 1.7 + ngc 6656 & m22 & -1.70 & 0.660 & 10 & 0.330 & 16 & 0.97 & 4.9 & 14 + ngc 6712 & & -1.02 & 0.557 & 7 & 0.338 & 2 & -0.66 & 3.5 + ngc 6715 & m54 & -1.49 & 0.584 & 95 & 0.351 & 33 & 0.75 & 18.9 & 15 + ngc 6723 & & -1.10 & 0.541 & 35 & 0.292 & 7 & 0.00 & 2.6 & 16 + ngc 6809 & m55 & -1.94 & 0.662 & 4 & 0.340 & 9 & 0.88 & 3.9 + ngc 6864 & m75 & -1.29 & 0.587 & 25 & 0.305 & 13 & -0.14 & 14.7 & 17 + ngc 6934 & & -1.47 & 0.574 & 68 & 0.302 & 10 & 0.21 & 12.8 + ngc 6981 & m72 & -1.42 & 0.568 & 36 & 0.320 & 7 & 0.07 & 12.9 & 18 + ngc 7006 & & -1.52 & 0.569 & 53 & 0.329 & 9 & -0.11 & 38.5 + ngc 7078 & m15 & -2.37 & 0.644 & 68 & 0.370 & 95 & 0.74 & 10.4 & 19 + ngc 7089 & m2 & -1.65 & 0.674 & 19 & 0.327 & 13 & 0.92 & 10.4 & 20 + ngc 7099 & m30 & -2.27 & 0.693 & 4 & & 0 & 0.94 & 7.1 & 21 + arp2 & & -1.75 & 0.584 & 8 & 0.292 & 1 & 0.53 & 21.4 & 8 + ic 4499 & & -1.53 & 0.580 & 63 & 0.350 & 34 & 0.12 & 15.7 + rup 106 & & -1.68 & 0.617 & 13 & & 0 & -0.82 & 18.5 + terzan 10 & & -1.00&0.660&10 & & & & 2.3&22 + 2mass - gc2 & & -1.08&0.590&3&&&&3.2&22 +        -2.25 & 0.677 & 0.72 & & 0.677 & 0.84 & & & 1.00 & & & 1.00 + -2.20 & 0.653 & 0.69 & & 0.653 & 0.82 & & 0.849 & 0.99 & & & 1.00 + -2.15 & 0.638 & 0.65 & & 0.639 & 0.80 & & 0.833 & 0.99 & & 0.833 & 0.99 + -2.10 & 0.630 & 0.62 & & 0.632 & 0.78 & & 0.819 & 0.99 & & 0.821 & 0.99 + -2.05 & 0.631 & 0.60 & & 0.633 & 0.77 & & 0.806 & 0.99 & & 0.808 & 0.99 + -2.00 & 0.641 & 0.56 & & 0.642 & 0.75 & & 0.794 & 0.98 & & 0.798 & 0.99 + -1.95 & 0.654 & 0.51 & & 0.656 & 0.72 & & 0.775 & 0.97 & & 0.772 & 0.98 + -1.90 & 0.667 & 0.43 & & 0.669 & 0.67 & & 0.749 & 0.96 & & 0.748 & 0.98 + -1.85 & 0.669 & 0.34 & & 0.674 & 0.62 & & 0.704 & 0.95 & & 0.705 & 0.97 + -1.80 & 0.664 & 0.29 & & 0.669 & 0.59 & & 0.672 & 0.95 & & 0.672 & 0.97 + -1.75 & 0.652 & 0.26 & & 0.656 & 0.57 & & 0.662 & 0.95 & & 0.664 & 0.97 + -1.70 & 0.639 & 0.23 & & 0.642 & 0.55 & & 0.659 & 0.95 & & 0.658 & 0.97 + -1.65 & 0.625 & 0.20 & & 0.627 & 0.53 & & 0.656 & 0.95 & & 0.658 & 0.97 + -1.60 & 0.602 & 0.13 & & 0.605 & 0.49 & & 0.648 & 0.94 & & 0.649 & 0.96 + -1.55 & 0.571 & 0.03 & & 0.574 & 0.42 & & 0.632 & 0.92 & & 0.632 & 0.95 + -1.50 & 0.553 & -0.09 & & 0.560 & 0.33 & & 0.609 & 0.89 & & 0.609 & 0.93 + -1.45 & 0.547 & -0.21 & & 0.558 & 0.24 & & 0.580 & 0.83 & & 0.581 & 0.90 + -1.40 & 0.549 & -0.35 & & 0.564 & 0.14 & & 0.560 & 0.77 & & 0.560 & 0.87 + -1.35 & 0.556 & -0.50 & & 0.572 & 0.02 & & 0.555 & 0.69 & & 0.555 & 0.82 + -1.30 & 0.559 & -0.64 & & 0.581 & -0.09 & & 0.553 & 0.58 & & 0.552 & 0.76 + -1.25 & 0.561 & -0.70 & & 0.586 & -0.14 & & 0.548 & 0.50 & & 0.549 & 0.71 + -1.20 & 0.563 & -0.75 & & 0.592 & -0.20 & & 0.545 & 0.40 & & 0.545 & 0.65 + -1.15 & 0.568 & -0.82 & & 0.601 & -0.26 & & 0.542 & 0.26 & & 0.543 & 0.57 + -1.10 & 0.574 & -0.88 & & 0.614 & -0.33 & & 0.543 & 0.06 & & 0.544 & 0.46 + -1.05 & 0.584 & -0.94 & & 0.631 & -0.42 & & 0.549 & -0.15 & & 0.550 & 0.33 + -1.00 & & & & 0.648 & -0.51 & & 0.562 & -0.38 & & 0.562 & 0.18 + -0.95 & & & & 0.667 & -0.59 & & 0.578 & -0.63 & & 0.577 & 0.02 + -0.90 & & & & 0.691 & -0.68 & & 0.595 & -0.86 & & 0.591 & -0.14 + -0.85 & & & & 0.713 & -0.76 & & & & & 0.600 & -0.25 + -0.80 & & & & 0.723 & -0.83 & & & & & 0.614 & -0.32 + -0.75 & & & & 0.747 & -0.89 & & & & & 0.627 & -0.39 + -0.70 & & & & 0.769 & -0.93 & & & & & 0.640 & -0.45 + -0.65 & & & & & & & & & & 0.649 & -0.52 + -0.60 & & & & & & & & & & 0.656 & -0.59 + -0.55 & & & & & & & & & & 0.669 & -0.66 + -0.50 & & & & & & & & & & 0.688 & -0.74 + -0.45 & & & & & & & & & & 0.708 & -0.81 + -0.40 & & & & & & & & & & 0.734 & -0.88 + -0.35 & & & & & & & & & & 0.761 & -0.94 +