during the last decade , there has been a dramatic paradigm shift on the the definition of the gc systems . despite the formerly accepted idea of chemical homogeneity
, the variations in the lighter elemental abundances in several gcs in our galaxy had been known for several decades .
the lick - texas group is one of those who undertook the systematic study of lighter elemental abundances in several gcs in our galaxy since 1990 s ( e.g. * ? ? ?
* ; * ? ? ?
thanks to the advent of high performance multi - object high - resolution spectrographs mounted on the large aperture telescopes , it is now possible for us to look into the detailed substructure of elemental abundance distributions of the milky way gc systems ( e.g. * ? ? ?
* ) . the decades - long lighter elemental variation issue in gc stars is now considered to be a generic feature of normal gcs in our galaxy , most likely engraved during the multi - phase normal gc formation ( e.g. * ?
* ; * ? ? ?
* ; * ? ? ?
contrary to the normal gc system , one of the key features of the peculiar gcs , such as @xmath3 cen , is the spread or the distinctive substructure in the metallicity distributions ( e.g. * ? ? ?
* ; * ? ? ?
* ) , where the heavy elements must have been supplied by supernovae ( sne ) .
to retain ejecta from energetic sne explosions , such peculiar gc systems must have been much more massive in the past and they are generally thought to be the remaining core of a disrupted dwarf galaxy and accreted to our galaxy later in time , expected from the hierarchical merging paradigm in the @xmath4cmd cosmological model ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
the existence of these peculiar gcs have major implications in the context of near field cosmology ( e.g. * ? ? ?
* ) ; are they one of the original building blocks to our galaxy , mitigating so - called the `` missing satellite problem ? ''
@xcite ; how did they relate to the formation of the galactic halo and numerous streamers ?
( e.g. * ? ? ?
these are examples of the outstanding problems that we have to challenge in the next decade .    to measure metallicity of stars in gcs from high - resolution spectroscopy , the lte ( local thermodynamic equilibrium )
analysis is being widely used for the sake of convenience , where the final results critically depend on the stellar atmosphere models with a few appropriate input stellar parameters , such as effective temperature , surface gravity and turbulent velocity , and the oscillator strengths of the absorption lines .
although simple , the derivation of stellar elemental abundances is not a trivial task even for nearby bright stars .
the recent study of @xcite may highlight the current situation .
they showed that the interferometric effective temperatures for nearby k giant stars do not agree with those from spectroscopic observations , suggesting a missing source of opacities in stellar atmosphere models .
the situation would be even worse for fainter stars , such as those in gcs .
another line of difficulty is that changes in surface gravity can mimic the chemical compositions in the regime of rgb stars in gcs ( for example , see * ? ? ?
* ) , where h@xmath5  is the major source of the continuum opacity and the h@xmath5  population varies with the electron pressure , therefore , the surface gravity of rgb stars .    the spread in the metallicity distribution of m22 has been a controversial topic for many years .
the recent several studies of the cluster have found that m22 has a bimodal heavy elemental abundance distribution ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
the high - resolution spectroscopic elemental abundance measurements of rgb stars in the peculiar gc m22 by @xcite and @xcite showed a distinctive bimodal metallicity distribution .
their results were based on the spectroscopic @xmath6  and @xmath7 , which require the excitation and ionization equilibria , @xmath8 and @xmath9 where @xmath10  and @xmath11  are iron abundances from and lines and @xmath12 is the excitation potential .
it has been frequently suspected that the iron abundance from the line suffers from non lte ( nlte ) effects ( see , for example , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
since metal - poor stars have much weaker metal absorption in the ultraviolet , more non - local ultraviolet flux can penetrate from the deeper layers .
this flux is vital in determining the ionization equilibrium of the atoms , resulting in deviations from lte @xcite . in this regard
, the traditional spectroscopic surface gravity determination method may be in error , which led @xcite to propose to use the photometric gravity for elemental abundance study of rgb stars in gcs , where the bolometric correction could be the dominant source of uncertainty for the photometric gravity . in principle , the use of the [ fe / h]@xmath13  abundance as for the metallicity scale of rgb stars in gcs is most likely an appropriate approach since  is by far the dominant species and , therefore , the number of  atoms is unaffected by the nlte effect . in practice , however , the [ fe / h]@xmath13  abundances of rgb stars sensitively depends not only on the surface gravity and effective temperature but also on the metallicity of the input atmosphere model , which also affects the continuum opacity
. we will show later that an iterative procedure is useful to reduce the error raised by the incorrect metallicity of the input atmosphere model .    in their recent study , ( * ? ? ?
* mu15 hereafter ) re - analyzed m22 rgb stars of @xcite using three different approaches : ( 1 ) spectroscopic @xmath6  and @xmath7 ( method 1 ) ; ( 2 ) spectroscopic @xmath6  and photometric @xmath7 ( method 2 ) ; and ( 3 ) photometric @xmath6  and @xmath7 ( method 3 ) .
they confirmed a bimodal iron abundance distribution of m22 rgb stars by @xcite , when they relied on spectroscopic @xmath6  and @xmath7 ( method 1 ) .
oddly enough , when they used photometric @xmath7 ( method 2 and 3 by mu15 ) , the allegedly well - established bimodal iron distribution of m22 disappeared in [ fe / h]@xmath13 .
using the photometric gravity is most likely a correct approach but how can this be interpreted ?
we will show later that it is likely that mu15 used incorrect surface gravity and the metallicity of input atmosphere models in their analysis and the separation in [ fe / h]@xmath13  can be brought out more fully if different methods to compute these parameters are used .    in this paper , we revisit the internal metallicity distribution of m22 .
we developed new methods to estimate the surface gravity and we found that there exists substantial metallicity difference between the @xmath0-@xmath1  and @xmath0-@xmath2  groups-@xmath1  ( calcium - weak ) and the @xmath0-@xmath2  ( calcium - strong ) groups are defined to be rgb stars with smaller or larger @xmath14 index values , respectively , at a given @xmath15 magnitude .
they are equivalent to the @xmath2-process - poor ( @xmath0-@xmath1 ) and the @xmath2-process - rich ( @xmath0-@xmath2 ) groups classified by @xcite . ] in m22 @xcite . throughout this paper
, metallicity refers to [ fe / h]@xmath13 , unless specified .
before turning to the metallicity distribution of m22 , we would like to revisit the critical issues on the lte analysis ; surface gravity , effective temperature and the metallicity of the input atmosphere model .      as mentioned above , @xcite suggested to use [ fe / h]@xmath13  derived from photometric gravity as metallicity of rgb stars .
they pointed out that [ fe / h]@xmath13  is essentially independent of nlte effects , such as overionization by non - local uv flux since is the dominant species in gc rgb stars .    in figure
[ fig : ki03 ] , we show plots of @xmath16[fe / h ]  against [ fe / h]@xmath17 , [ fe / h]@xmath13 , @xmath6  and @xmath7 for the six group 1 clusters  and @xmath7 based on colors and absolute magnitudes . ] of @xcite . of particular concern is figure  [ fig : ki03 ] ( a ) and ( b ) , where each gc is showing its own correlation between [ fe / h]@xmath17  versus @xmath16[fe / h ]  ( or [ fe / h]@xmath13  versus @xmath16[fe / h ] ) .    in figure
[ fig : fit ] , we show plots of @xmath7 versus @xmath6 , @xmath6  versus @xmath18 , and @xmath7 versus @xmath18  for the six group 1 gcs . also shown are victoria - regina model isochrones for 12 gyr @xcite .
to estimate the @xmath19  level of the model isochrones , we used our previous relation @xcite , @xmath20 } + 1.5 ) + ( 0.52 \pm 0.13).\ ] ] note that this metallicity - luminosity relation for rr lyrae variables gives @xmath21 = 18.54 @xmath22 0.13 mag for lmc .
it also should be emphasized that the adopted age of the model isochrones does not affect our results presented in this work .    in plots of @xmath7 versus @xmath6  and @xmath6  versus @xmath18 ,
the loci of the model isochrones are in excellent agreement with the observations .
it should not be a surprise that , at a given temperature , the surface gravity of the metal - poor stars is lower and @xmath18  is smaller than those of the metal - rich stars .
this implies that applying a single @xmath6  versus @xmath7 relation for groups of stars with heterogeneous metallicity may result in incorrect surface gravity estimates and , as a consequence , incorrect elemental abundances .
this is especially true when one adopt [ fe / h]@xmath13  as metallicity of gc rgb stars .
the  line opacity does not vary with surface gravity since the almost all iron atoms are populated in the first ionized level . on the other hand ,
the h@xmath5  continuum opacity is sensitively dependent on the electron pressure and , therefore , surface gravity . as a consequence , in the regime of the linear part of the curve of growth , at fixed metallicity ,
the equivalent widths ( ews ) of  lines increases with decreasing surface gravity . in the figure ,
the surface gravity dependency on the metallicity of the metal - poor rgb stars at a fixed temperature is given by @xmath23}}{\partial \log g } \approx 1.8\ { \rm dex}.\ ] ] in other words , a change in the surface gravity by 0.1 in the cgs unit corresponds to a change in metallicity by about 0.2 dex in the figure .    on the other hand , a rather tight relation in @xmath7 versus @xmath18   can be found , suggesting that the @xmath18  magnitude can be a reddening- and distance - independent surface gravity indicator for rgb stars in gcs .
( see also , * ? ? ? * for the use of @xmath24 as a temperature indicator of heavily reddened metal - rich gc palomar  6 ) .
but caution is advised that the observed @xmath18  magnitude can be vulnerable to the foreground differential reddening effect .      as mentioned previously
, @xcite suggested to use the photometric surface gravity .
however , ironically , under the condition that lte is not valid , they concluded to use spectroscopic temperature from  lines .
we would like to discuss the utility of the spectroscopic effective temperature under the assumption that the conclusion made by @xcite is valid .    in figure
[ fig : log ] , we show plots of the @xmath10  against the excitation potential for ngc  6752-mg10 by @xcite .
as they noted , the quality of their spectra is superb ; a resolving power of @xmath25 = 110,000 and s / n @xmath26 150 per pixel near 5140  . to perform a lte abundance analysis , we used the 2014 version of the stellar line analysis program moog @xcite and we interpolated @xmath27-enhanced kurucz atmosphere models with new opacity distribution functions using a fortran program kindly provided by dr .
mcwilliam ( 2005 , private communication ) . using the weak lines , @xmath28 )  @xmath29 , we derived a spectroscopic temperature for this star by forcing the condition of the excitation equilibrium , i.e.  @xmath30@xmath10/@xmath31 = 0 , and we obtained the effective temperature of 4275 k. note that our effective temperature for this star is slightly cooler than those from the line - by - line differential analysis with respect to ngc  6752-mg9 and ngc  6752-mg6 by @xcite , 4291 k and 4295 k , respectively .
we derived the linear fits to the data , assuming @xmath10  @xmath32 slope@xmath33 , and we show them with thin dashed lines in the figure . in the parentheses in each panel
, we also show the @xmath7 value in the cgs unit , the slope in the excitation potential versus the iron abundance , [ fe / h]@xmath17  and [ fe / h]@xmath13 .
as can be seen , once the spectroscopic effective temperature is correctly determined , the assumption of the excitation equilibrium holds for rather wide range of the surface gravity , @xmath34 1.20 in this case .
therefore , it can be concluded that the slightly incorrect input surface gravity does not significantly affect the spectroscopic effective temperature in our results presented here .
it is well - known fact that the inferred metallicity of gc rgb stars from high resolution spectroscopy is critically dependent on the stellar parameters of the input atmosphere model . here
, we show how [ fe / h]@xmath17  and [ fe / h]@xmath13  behave against the changes in the input stellar parameters .
we also would like to demonstrate the importance of the iterative derivation of the metallicity .    in figure
[ fig : fehvspam ] , we show [ fe / h]@xmath17  and [ fe / h]@xmath13  of 9 rgb stars in ngc  6752 by @xcite . in each panel ,
the red crosses denote our spectroscopic [ fe / h]@xmath17  and [ fe / h]@xmath13 which were derived using the weak  and  lines measured by @xcite .
it should be noted that our [ fe / h ]  values are consistent with those of @xcite . using our spectroscopic @xmath6 , @xmath7 , and [ fe / h ]  as reference grids , we examine how the changes in the stellar parameters of the input atmosphere models affect the resultant metallicity .
we run moog using the kurucz atmosphere models , whose stellar parameters are different from the reference grids by @xmath16@xmath6  = @xmath22200 k , @xmath16[fe / h ]  = @xmath221.0 , and @xmath35 = @xmath220.3 .
we use the [ fe / h]@xmath13  abundances returned from moog as the reference metallicity and we calculate new atmosphere models , with which we run moog again .
we iterate this process five times and we show our results in figure  [ fig : fehvspam ] with blue solid lines .    the figure shows that the inferred [ fe / h]@xmath17  abundances with sufficient numbers of iterations are only affected by the changes in the effective temperature . as shown , [ fe / h]@xmath17  abundances of individual stars against the changes in the metallicity and the surface gravity of the input model atmosphere converge to their spectroscopic [ fe / h]@xmath17  values within 2 or 3 iterations
. however , the iterations with the effective temperature offsets of @xmath22200k fail to regain spectroscopic [ fe / h]@xmath17 , implying incorrect @xmath6  estimate results in irrecoverable deviations in the derived metallicity .
the [ fe / h]@xmath13  behaves differently against the changes in input stellar parameters .
the inferred [ fe / h]@xmath13  abundances with the effective temperature offsets shifted in the opposite direction of the changes in the [ fe / h]@xmath17  abundance
. a rather simple explanation of the temperature effect in cool stars can be found in @xcite , for example .
similar to [ fe / h]@xmath17 , the spectroscopic metallicity can be regained with the iterative derivation of the [ fe / h]@xmath13  abundances against the changes in the metallicity of the input atmosphere model but , however , the difference in [ fe / h]@xmath13  between the spectroscopic metallicity and that from the first iteration is much larger than that can be seen in [ fe / h]@xmath17 .
as we have already discussed earlier , the [ fe / h]@xmath13  abundance sensitively depends on the surface gravity of the input atmosphere model .
the figure shows that the spectroscopic metallicity can not be regained with incorrect surface gravity .
we also performed the same procedures using the [ fe / h]@xmath17  abundances of individual stars as reference metallicity to calculate the input atmosphere model used in the next iteration and we obtained the same results shown above .
figure  [ fig : delfeh ] summarize our exercise .
for both [ fe / h]@xmath17  and [ fe / h]@xmath13 , the inferred metallicity from the first iteration with offsets in the stellar parameters could be very different from those with correct stellar parameters . however , some discrepancies vanish after 2 or 3 iteration processes .
our exercise demonstrates that , with the iterative procedure , the [ fe / h]@xmath17  abundance depends only on the effective temperature , while the [ fe / h]@xmath13  abundance depends both on the effective temperature and the surface gravity .
it also shows the importance of having correct stellar parameters , especially for the lte analysis of [ fe / h]@xmath13  abundances , as we have discussed earlier .
the vivid evidence of the bimodal metallicity distribution of rgb stars in m22 from narrow - band photometry and high - resolution spectroscopy can be found in @xcite , @xcite and @xcite . in our earlier study of the cluster
, we extensively discussed that there are several observational lines of evidence that can not be easily explained without invoking a bimodal metallicity distribution between two groups of stars , namely the @xmath0-@xmath1  and the @xmath0-@xmath2  groups as shown in figure  [ fig : m22cmd ] ( e.g. , see figures  9 , 13  21 of * ? ? ? * and references therein ) ;    * the  h&k absorption strengths rgb stars at a given @xmath15 magnitude in m22 from both narrow - band photometry @xcite and low - resolution spectroscopy @xcite show a bimodal distribution . * the infrared  triplet by @xcite
also show a bimodal distribution among rgb stars in m22 . * the @xmath36 versus @xmath15 cmd as shown in figure  [ fig : m22cmd ]
( see also figure  19 of * ? ? ?
* ) also requires a bimodal metallicity distribution in m22 rgb stars .
the variation in the lighter elements only , such as cno , can not explain this distinct double @xmath36 rgb sequences of m22 .
the differential foreground reddening effect can not reproduce the observed multi - color cmds accordingly @xcite . *
the @xmath15 magnitude of the rgb bump , @xmath37 , of the @xmath0-@xmath2  group is significantly fainter than that of the @xmath0-@xmath1  group , which strongly suggests that the @xmath0-@xmath2  group is more metal - rich than the @xmath0-@xmath1  group is .
the difference in the @xmath37  between the two groups can not be explained by the differential foreground reddening effect @xcite . *
the slope of the @xmath0-@xmath2  rgb stars in the @xmath38 versus @xmath15 cmd is significantly larger than that of the @xmath0-@xmath1  rgb stars , indicative of the metal - rich nature of the @xmath0-@xmath2  group @xcite .
* the cn - ch positive correlation superposed on two separate cn - ch anticorrelations @xcite can be expected naturally if m22 is composed of two groups of stars with heterogeneous metallicities @xcite . *
finally , another piece of evidence can be found in the metallicity distribution of the blue horizontal branch ( bhb ) stars in m22 @xcite . in figure
[ fig : hb ] , we show the metallicity distributions for m22 bhb stars .
both the lte and the nlte treatments show the similar degree of metallicity spreads in [ fe / h]@xmath17  and [ fe / h]@xmath13(six stars for [ fe / h]@xmath17  and seven stars for [ fe / h]@xmath13 ) .
as shown in table  3 of @xcite , their derived [ fe / h]@xmath13  value is less sensitive to the effective temperature and to the surface gravity , @xmath16[fe / h]@xmath13= @xmath220.06 dex and @xmath220.01 dex for @xmath16@xmath6  = @xmath22170 k and @xmath39 = @xmath220.20 respectively , and we are likely seeing the real metallicity spread in m22 bhb stars .
as mentioned above , mu15 re - analyzed m22 rgb stars of @xcite using three different approaches .
we show 17 rgb and asymptotic giant branch ( agb ) stars studied by mu15 in figure  [ fig : m22cmd ] . among 12 rgb stars tagged by mu15 , six rgb stars belong to each of two rgb groups , according to our previous classification of rgb stars based on the @xmath14 index at a given @xmath15 magnitude .
figure  [ fig : pam ] shows a plot of the @xmath40 @xmath6  versus @xmath7 of spectroscopic target stars , where the filled symbols are for method 1 ( spectroscopic @xmath6  and @xmath7 ) and the open symbols are for method 2 ( spectroscopic @xmath6  and photometric @xmath7 ) of mu15 . also shown
are victoria - regina model isochrones for 12 gyr with [ fe / h ] = @xmath411.40 , @xmath411.70 and @xmath412.00 dex @xcite .
the most conspicuous feature of the figure is that the positions of stars from method 2 by mu15 are rather well aligned in a narrow strip between model isochrones with [ fe / h ] = @xmath411.40 and @xmath411.70 dex on the @xmath6
@xmath7  plane , which is most likely due to the adopted _ single _ @xmath6
@xmath7  relation by mu15 .
it should be worth to point out that , because the cmd of mu15 has a single rgb sequence , their photometric effective temperature and surface gravity will also define a single isochrone .
the broad band colors adopted by them are not sensitive to small [ fe / h ]  differences . on the other hand ,
positions of stars from method 1 occupy rather wide ranges in surface gravity at given effective temperature .
we performed a lte abundance analysis using ews of m22 rgb stars measured by mu15 to calculate the the mean iron abundance dependence on model atmospheres . in table
[ tab : dep ] , we show our result for eight rgb stars in m22 with 0.5 @xmath42 1.5 , where @xmath43 is the photometric surface gravity in the cgs unit . as shown in the table , the change in the surface gravity by @xmath44 0.2 results in @xmath16[fe / h]@xmath13  @xmath45 0.1 dex at the fixed @xmath6 , in the sense that the [ fe / h]@xmath13  value increases with the surface gravity
.. ] if we take this at face value , a single @xmath6  @xmath7  relation , which is valid only for mono - metallic stellar systems such as normal gcs in our galaxy , may be responsible for the narrow uni - modal [ fe / h]@xmath13  distribution of m22 rgb stars as claimed by mu15 as shown in their figure  2 .    in figure
[ fig : method2 ] , we show the metallicity distribution of m22 rgb stars using table  4 of mu15 ( method 2 ) .
note again that we use 8 rgb stars with 0.5 @xmath42 1.5 only ( corresponding to 11.5 @xmath46 12.75 mag in figure  [ fig : m22cmd ] ) to avoid the potential effect raised by very different surface gravity in the stellar atmosphere model calculations . according to our population classification scheme for m22 based on our @xmath14 index at a given @xmath15 magnitude ,
stars 61 , 71 , 200068 and 200076 belong to the @xmath0-@xmath1  group , while stars 51 , 88 , 200025 and 200101 belong to the @xmath0-@xmath2  group .
figure  [ fig : method2 ] ( a ) shows that each rgb group has different mean iron abundance both in [ fe / h]@xmath17  and [ fe / h]@xmath13 .
it is very interesting to note that the differences in the mean iron abundances between the two groups are larger than a 2.5 @xmath47 level both in [ fe / h]@xmath17  and [ fe / h]@xmath13 , although the separation in the mean [ fe / h]@xmath13  values is as small as @xmath45 0.05 dex . in table
[ tab : ks ] , we show our results ( mu15 m2 ) .    figure  [ fig : method2 ] ( b )
( e ) show the cumulative metallicity distributions and the generalized histograms of the metallicity distributions for each group .
as shown , the [ fe / h]@xmath17  distribution show two distinct peaks , while that for [ fe / h]@xmath13  shows a broad single peak in the overall metallicity distribution , similar to those obtained by mu15 .
but it should be emphasized that the [ fe / h]@xmath13  distribution by mu15 is composed of two separate mono - metallic distributions ; only the mean [ fe / h]@xmath13  values from the method 2 by mu15 for each group of stars happen to be similar .
we performed a student s @xmath48-test to see if the metallicity distributions of the two groups of stars are identical .
we obtained that the significance levels to reject the hypothesis that the mean [ fe / h ]  values of the @xmath0-@xmath1  and the @xmath0-@xmath2  groups are identical , are 1.93 % and 8.98 % for [ fe / h]@xmath17  and [ fe / h]@xmath13 , respectively .
we also performed a randomization test .
the significance levels to reject the hypothesis for being an identical metallicity distribution from bootstrap method are 0.00 % for both [ fe / h]@xmath17  and [ fe / h]@xmath13 , strongly suggest that the metallicity distributions for the @xmath0-@xmath1  and the @xmath0-@xmath2  groups by mu15 are not statistically identical .
it should not be a surprise because @xcite already discussed many aspects of heterogeneous nature between the two groups as summarized in
[ sec3 ] .
figure  [ fig : method2 ] ( f  g ) shows @xmath16[fe / h ]  (= [ fe / h]@xmath13  @xmath41 [ fe / h]@xmath17 ) against [ fe / h]@xmath17  and [ fe / h]@xmath13  and they are very intriguing .
it should be reminded that we chose stars with a narrow range of @xmath15 magnitude for both groups to avoid the potential effect raised by the very different surface gravity .
as shown in the figure , the discrepancies between [ fe / h]@xmath17  and [ fe / h]@xmath13  of the @xmath0-@xmath1  group rgb stars are preferentially much greater than those of the @xmath0-@xmath2  group , reaching as large as @xmath16[fe / h ]  @xmath45 0.4 dex . for comparison
, we show @xmath16[fe / h ]  of m22 along with six group 1 gcs by @xcite in figure  [ fig : ki03 ] , where the much greater @xmath16[fe / h ]   values of the @xmath0-@xmath1  rgb stars in m22 can be clearly seen .
if the two groups of stars have the same metallicity and the same surface gravity so that they suffer similar degree of nlte effect , one would expect to see the similar degree of @xmath16[fe / h ]  for both groups , in sharp contrast to the results of method 2 by mu15 .
we will show later that , with a mock peculiar gc , the incorrect surface gravity estimate by mu15 and the incorrect metallicity of input model atmospheres are responsible for the discrepancy in @xmath16[fe / h ] .
another aspect needs to consider is the initial metallicity used during the atmosphere model calculations . in figure  [ fig : mu15iter ] , we show plots for eight m22 rgb stars , similar to figure  [ fig : fehvspam ] for ngc  6752 . in parentheses of each panel ,
we show the method for the metallicity derivation by mu15 ( m1 or m2 ) and the reference metallicity for the iterations ( i for [ fe / h]@xmath17  and ii for [ fe / h]@xmath13 ) . in the figure ,
the crosses denote the mu15 s [ fe / h]@xmath17  and [ fe / h]@xmath13  values from m1 and m2 methods and blue and red colors denote the @xmath0-@xmath1  and @xmath0-@xmath2rgb stars , respectively . using [ fe / h]@xmath17  and [ fe / h]@xmath13  abundances as trial metallicities and @xmath6  and @xmath7 from m1 and m2 methods by mu15
, we performed the iterative derivations of [ fe / h]@xmath17  and [ fe / h]@xmath13  for individual stars .
as shown , [ fe / h]@xmath17  values do not vary significantly with the number of iterations with respect to the spectroscopic [ fe / h ] , meaning that the effective temperature adopted by mu15 is correct .
on the other hand , the discrepancy in [ fe / h]@xmath13  for @xmath0-@xmath2  rgb stars from method 2 are preferentially larger , indicating that the surface gravities of the @xmath0-@xmath2  rgb stars adopted by mu15 is most likely underestimated .
the figure also indicates that the separation in the [ fe / h]@xmath13  abundances between the @xmath0-@xmath1  and the @xmath0-@xmath2  groups becomes larger with the iteration processes .
finally , comparisons of ews between the two groups may also help to elucidate the underlying metallicity distributions of the cluster . in figure
[ fig : ew](a ) , we show the line - by - line ew difference between stars 88 ( @xmath0-@xmath2 ) and 200076 ( @xmath0-@xmath1 ) .
both stars have similar visual magnitudes and colors , ( @xmath15 , @xmath49 ) = ( 12.54 , 0.89 ) for the star 88 and ( 12,39 , 0.90 ) for the star 200076 .
therefore , if there is no significant foreground differential reddening and if the both stars have the same metallicity , they should have similar @xmath6 , @xmath7 and furthermore the similar ew strengths .
as shown in the figure , the ews of the @xmath0-@xmath2  group star 88 are systematically larger than those of the @xmath0-@xmath1  group stars 200076 . a comparison between the mean ews of the four @xmath0-@xmath2  group stars ( 51 , 88 , 200025 and 200101 ) and the four @xmath0-@xmath1  group stars ( 61 , 71 , 200068 and 200076 ) show the same trend that the mean ews of the @xmath0-@xmath2  group are larger than those of the @xmath0-@xmath1  group , 12.0 @xmath22 0.4 m  for lines and 1.8 @xmath22 0.4 m  for lines .
note that the eight stars above have similar visual magnitudes and colors and they should have similar @xmath6 , @xmath7 and [ fe / h ]  if they belong to a single stellar population .
a simple explanation why  lines are less sensitive to changes in metallicity as follows . since both the fraction of  atoms and the h@xmath5  continuum opacity of rgb stars depend on the electron pressure ( i.e. , metallicity or surface gravity ) and , furthermore , the two effects are expected to cancel out , the ews of  lines grow with metallicity at fixed effective temperature . on the other hand , the  atoms are the dominant species and only the electron pressure has an effect on the h@xmath5  continuum opacity , which has an opposite effect on the growth of ews with metallicity . therefore , the ews of  lines grows at a slower rate .
we suspect that the metallicity measurements by mu15 may be slightly incorrect and we devise new methods to derive metallicity of rgb stars in appropriate and consistent manners .
in our previous study , @xcite showed that a combination of two normal gcs , m55 and ngc  6752 , can reproduce many aspects of peculiar photometric characteristics of m22 . in figure
[ fig : m55n6752cmd ] , we show a composite cmds for m55 and ngc  6752 , which may highlight the importance of the choice of photometric passbands to distinguish multiple stellar populations in gcs . for ngc  6752 stars ,
we add offsets of 0.030 , @xmath410.010 , @xmath410.005 , and 0.700 mag in ( @xmath49 ) , @xmath36 , @xmath14 , , @xmath50 , and @xmath51 . ] and @xmath15 , respectively , in order to place ngc  6752 stars in the m55 colors and magnitude scale . in the figure
, we also show rgb stars studied from high - resolution spectroscopy of the clusters @xcite and the victoria - regina model isochrones for 12 gyr with [ fe / h ]  = @xmath411.84 and @xmath411.53 .    despite the difference in metallicity ,
the rgb sequence of ngc  6752 is in excellent agreement with that of m55 in @xmath52 versus @xmath15 cmds , and it is difficult to discern different populations . in such a case , one can be easily misled to adopt a single @xmath6
@xmath7  relation for heterogeneous stellar populations . on the other hand ,
the distinct double rgb sequences can be clearly seen in @xmath36 versus @xmath15 and @xmath14 versus @xmath15 cmds , where the necessity for double @xmath6
@xmath7  relations is obvious .
first , we derived the metallicity distributions of m55 and ngc  6752 by employing two separate @xmath6
@xmath7  relations , following the procedure recommended by @xcite . with the ew measurements by @xcite for and lines
, we performed a lte abundance analysis . applying the color - temperature relation and the equation for the bolometric correction given by @xcite , we calculated the photometric effective temperature and the surface gravity using our own strmgren  photometry . during our calculations
, we adopted [ fe / h ]  = @xmath411.95 and @xmath411.55 as the input metallicities for m55 and ngc  6752 , respectively @xcite , and we used the distance moduli and foreground reddening values by @xcite . during our analysis , we used weak lines only , @xmath28 )  @xmath29 , for both and in order to minimize the effect of the adopted micro - turbulent velocity on the derived metallicity .
we show our results in table  [ tab : mock ] and figure  [ fig : photorelation ] .
our [ fe / h]@xmath13  measurements for the clusters are in good agreement with those by @xcite , @xmath16[fe / h]@xmath13  = @xmath410.01 @xmath22 0.03 dex for m55 and @xmath16[fe / h]@xmath13  = 0.05 @xmath22 0.02 dex for ngc  6752 in the sense of current work minus those of @xcite .
these small differences are thought to be mainly due to slightly different final color - temperature relations .
@xcite used the @xmath53 color - temperature relations given by @xcite to derive the initial @xmath6  but they applied their own @xmath6  versus @xmath15 magnitude relation to derive their final adopted @xmath6 .    in figure
[ fig : photorelation ] ( b ) and ( c ) , we show @xmath16[fe / h ]  (= [ fe / h]@xmath13  @xmath41 [ fe / h]@xmath17 ) against [ fe / h]@xmath17  and [ fe / h]@xmath13 , where the extent of @xmath16[fe / h ]  for both clusters are in good agreement with those of six gcs studied by @xcite as shown in figure  [ fig : ki03 ] .
figure  [ fig : photorelation ] ( d ) shows a plot of @xmath6  versus @xmath7 along with the victoria - regina model isochrones for 12 gyr with [ fe / h ]  = @xmath411.84 and @xmath411.55 @xcite . at a given @xmath6 ,
the stars in the metal - poor cluster m55 have lower surface gravity , although the mean difference in @xmath7 between m55 and ngc  6752 is not as large as that can be inferred from model isochrones . using separate relations , the differences in the mean [ fe / h]@xmath17  and [ fe / h]@xmath13are 0.44 @xmath22 0.02 and 0.47 @xmath22 0.03 , respectively , and they are in excellent agreement .
following the same procedure described earlier , we performed iterative derivations of [ fe / h]@xmath17  and [ fe / h]@xmath13  for the clusters .
as shown in figure  [ fig : photorelation ] and table  [ tab : mock ] , [ fe / h]@xmath17  and [ fe / h]@xmath13  from iterative processes are in good agreement with those without the iterative process .
we also make use of the ew measurements for the  and  lines for ngc  6752 rgb stars by @xcite .
we show our results in figure  [ fig : m55n6752comp2 ] and table  [ tab : mock ] .
note that the oscillator strengths for the individual lines are slightly different between those adopted by @xcite and @xcite . for our results presented in figure  [ fig : m55n6752comp2 ]
, we used the oscillator strengths by @xcite because they provided more lines .
as summarized in table  [ tab : mock ] , [ fe / h]@xmath17  and [ fe / h]@xmath13  abundances using @xmath54-values by @xcite and @xcite are slightly different , but the differences in the mean values are no larger than 0.04 dex .
therefore , it can be said that the choice of the set of oscillator strengths does not affect our primary results presented here . the differences in [ fe / h]@xmath17   ( @xmath16[fe / h]@xmath17  = [ fe / h]@xmath17@xmath55
@xmath41 [ fe / h]@xmath17@xmath56 ) and [ fe / h]@xmath13(@xmath16[fe / h]@xmath13  = [ fe / h]@xmath13@xmath55 @xmath41 [ fe / h]@xmath13@xmath56 ) are 0.38 - 0.45 dex and 0.44 - 0.50 dex , respectively .      assuming that our mock gc ( i.e.  m55 + ngc  6752 ) is a mono - metallic gc with [ fe / h ]  @xmath57 dex ( i.e.  that of ngc  6752 ) , we calculated photometric effective temperatures and surface gravities of individual stars using the @xmath52 color - temperature relation and the equation for the bolometric correction given by @xcite .
using the weak and lines only , @xmath28 )  @xmath58 , we derived the metallicity of individual stars in both clusters and we show our results in table  [ tab : mock ] and figure  [ fig : photorelation ] .
the mean [ fe / h]@xmath17  abundance of m55 remains unchanged , while that of [ fe / h]@xmath13  increases almost 0.15 dex compared to the results from correct input stellar parameters for m55 presented in
[ sec5.1 ] .
the difference in the effective temperature between the two methods ( i.e.  two separate relations versus a single relation for m55 and ngc  6752 ) is negligibly small , @xmath59  = 12 k , in the sense that the mean effective temperature of m55 rgb stars is slightly warmer when the correct @xmath52 color - temperature relation is used . as we have discussed earlier ,
the [ fe / h]@xmath17  abundance is relatively insensitive to changes in the surface gravity and the metallicity of the input model atmosphere since the effects due to change in the number of  species and that in h@xmath5  continuum opacity are expected to cancel .
however , the surface gravity becomes larger , @xmath39 @xmath45 0.1 , and the metallicity of the input model atmospheres is higher ( from @xmath411.95 to @xmath411.55 dex ) when a single relation with respect to ngc  6752 is used .
both effects greatly enhance the h@xmath5  continuum opacity , while the fraction of  to the total number of iron atoms is unaffected since  is by far the dominant species . as a consequence , the mean [ fe / h]@xmath13  abundance of m55 appears to be enhanced at given ews .    in figure
[ fig : photorelation ] ( l ) , we show a plot of @xmath6  versus @xmath7 with a single relation , where rgb stars in both clusters are aligned well on a single locus .
figure  [ fig : photorelation ] ( j ) and ( k ) show @xmath16[fe / h ]  against [ fe / h]@xmath17  and [ fe / h]@xmath13  showing large discrepancies in m55 rgb stars , reminiscent of the m22 @xmath0-@xmath1  stars from method 2 of mu15 as shown in figure  [ fig : method2 ] .
the metallicity distributions of [ fe / h]@xmath17  and [ fe / h]@xmath13  from a single @xmath6
@xmath7  relation shown in figure  [ fig : m55n6752comp1 ] ( b ) and ( c ) are intriguing , since the two peaks in the [ fe / h]@xmath13  distribution become less conspicuous with a single relation with respect to ngc  6752 .
our result with a mock gc strongly suggests that applying a single photometric relation in order to derive the effective temperatures and surface gravities of individual stars in peculiar gcs with heterogeneous metallicities and perhaps ages , such as m22 , may result in slightly incorrect metallicity scales and distributions .
we also performed a traditional analysis using the spectroscopic effective temperatures and the surface gravities , which requires the excitation and ionization equilibria of iron abundances .
our results are shown in figures  [ fig : m55n6752comp1 ] ( i  l ) and [ fig : m55n6752comp2 ] ( i  l ) .
as shown in table  [ tab : mock ] , our mean spectroscopic [ fe / h ]  values are in good agreement with those of @xcite and @xcite . the difference in metallicity between m55 and ngc  6752 is @xmath16[fe / h]@xmath13  = 0.32  0.42 dex , depending on the data sets .
it should be mentioned that our spectroscopic [ fe / h ]  value of ngc  6752 from @xcite is about 0.1 dex higher than that from @xcite , which is consistent with the iron abundances of [ fe / h ]  = @xmath411.56 dex @xcite and @xmath411.65 dex @xcite for the cluster .
the origin of this discrepancy of the mean metallicity of ngc  6752 is beyond the scope of this study and we decline to discuss this matter further .
as shown in figure  [ fig : fit ] , the model isochrones can provide a useful means to derive the stellar parameters . in figure  [ fig:17gc ]
, we show the similar plots for 17 gcs from the homogeneous elemental abundance study by @xcite .
also shown are the victoria - regina isochrones for the age of 12 gyr and they appear to be in excellent agreement with observations .    we devise a new strategy to derive evolutionary surface gravities of rgb stars in gcs under the assumption that the excitation equilibrium of lines is applicable in such stellar atmospheres and , furthermore , excitation equilibrium holds for rather wide range of the surface gravity .
as shown in figure  [ fig : fit ] and [ fig:17gc ] , at a given effective temperature , the surface gravity increases with metallicity , and , as a consequence , the metallicity , especially [ fe / h]@xmath13 , without the proper estimates of the surface gravity may not be correct .    using the spectroscopic temperature and the [ fe / h]@xmath13  abundance from photometric method as initial input parameters , we interpolated the victoria - regina model isochrones to obtain the evolutionary surface gravity at the fixed effective temperature
. then we derive the updated metallicity by running moog using the model atmosphere with the spectroscopic effective temperature and the evolutionary surface gravity in an iterative manner until the derived metallicity converged to within the internal measurement error between consecutive measurements , which usually requires 2 to 3 iterations .
we show our new stellar parameters in figures  [ fig : m55n6752comp1 ] ( p ) and [ fig : m55n6752comp2 ] ( p ) and metallicity distributions in figures  [ fig : m55n6752comp1 ] ( n  o ) and [ fig : m55n6752comp2 ] ( n  o ) .
the difference in metallicity between m55 and ngc  6752 becomes @xmath16[fe / h]@xmath13  = 0.44  0.56 dex and our results are shown in table  [ tab : mock ] .
we also note that using the photometric temperatures of individual stars does not change our results presented here .
this approach should be reddening- and distance - independent and therefore it would be useful to derive surface gravity of gc stars with varying foreground reddening , such as rgb stars in m22 .      as shown in figure  [ fig : fit ] and [ fig:17gc ] , at a given metallicity , the @xmath15 magnitude differences from the hb , @xmath18 , of individual stars in gcs are well correlated with the surface gravities and the effective temperatures . similar to the previous approach , at a given @xmath18  and metallicity
we determine the evolutionary surface gravity and effective temperature simultaneously by interpolating the victoria - regina model isochrones . for this purpose ,
we use our own photometry of the clusters and [ fe / h]@xmath13  derived from the photometric stellar parameters as an initial guess as have done previously .
then we derive the updated metallicity by running moog using the model atmosphere with the evolutionary effective temperature and the surface gravity in an iterative manner until the derived metallicity converged to within the internal measurement error between consecutive measurements .
we show our results in figures  [ fig : m55n6752comp1 ] ( q  t ) and [ fig : m55n6752comp2 ] ( q  t ) .
this approach provides similar results as those from the spectroscopic method and the method relying on the evolutionary surface gravity .
the difference in metallicity between m55 and ngc  6752 becomes @xmath16[fe / h]@xmath13  = 0.44  0.52 dex as shown in table  [ tab : mock ] .
the merit of using @xmath18  is that it is also a reddening- and distance - independent parameter
. however it can be vulnerable to the differential foreground reddening effect of the individual stars .
it should be kept in mind that the main idea to deliver in our study is to demonstrate the importance of having appropriate stellar parameters for the lte abundance analysis in multiple stellar populations .
following the same procedures as for m55 and ngc  6752 , we derive the iron abundances of the two groups of stars in m22 in four different manners .
our results are consistent with the idea that the two groups of stars in m22 have different mean iron abundances as @xcite , @xcite and @xcite already showed .
we derive the metallicity of m22 rgb stars based on the photometric effective temperature and surface gravity from our strmgren  photometry of the cluster using the relations by @xcite . during our calculations
, we adopted the apparent visual distance modulus of 13.60 mag , @xmath60  = 0.34 and [ fe / h ]  = @xmath411.65 for m22 @xcite . as it was done before , we made use of the weak lines only , @xmath28 )  @xmath29 , for both and in order to minimize the effect of the adopted micro - turbulent velocity on the metallicity . in table  [ tab : ks ] and figures  [ fig : m22comb ] and [ fig : m22combdiff ] , we show our results . for non - differential analysis ,
the differences in the mean metallicity are @xmath16[fe / h]@xmath17  = 0.239 @xmath22 0.057 dex and @xmath16[fe / h]@xmath13  = 0.096 @xmath22 0.048 dex without iteration , and @xmath16[fe / h]@xmath17  = 0.233 @xmath22 0.048 dex and @xmath16[fe / h]@xmath13  = 0.108 @xmath22 0.052 dex after fifth iteration .
note that our results are consistent with those from the method 2 by mu15 . in panels ( b ) and ( c ) of figures  [ fig : m22comb ] and [ fig : m22combdiff ] , we show empirical distributions of the mean [ fe / h]@xmath17  and [ fe / h]@xmath13  for the two populations from the bootstrap method , strongly suggest that the metallicity distributions of the two groups of stars in m22 are different . as shown in table  [ tab : ks ] , the significance levels to reject the hypothesis that the mean [ fe / h]@xmath17  values of the @xmath0-@xmath1  and the @xmath0-@xmath2  groups are identical are lower than 1% . however , those for [ fe / h]@xmath13  are rather large , @xmath45 12 % , for the non - differential analysis . it should not be confused that this rather large significance levels do not indicate that two groups of stars in m22 belong to the same population , but the lte analysis of the heterogeneous groups of star with a single @xmath6
relation may be in error .
we also calculate the line - by - line differential iron abundances since the numbers of iron lines being measured by mu15 for individual stars are different .
we selected the star 51 to be the reference star since its @xmath6  and @xmath7 are close to the average for the sample . also the stellar parameters for this star
both from the photometric and spectroscopic methods by mu15 agree well as shown in figure  [ fig : pam ] . for our differential abundance measurements
, we did not adjust the stellar parameters of individual stars with respect to the reference star as have done by @xcite , for example , and we intended to calculate the proper metallicity offset differences among the sample stars with given stellar parameters .
the differences in the mean metallicity from the differential analysis are @xmath16[fe / h]@xmath17  = 0.203 @xmath22 0.039 dex and @xmath16[fe / h]@xmath13  = 0.101 @xmath22 0.045 dex without iteration , and @xmath16[fe / h]@xmath17  = 0.220 @xmath22 0.043 dex and @xmath16[fe / h]@xmath13  = 0.129 @xmath22 0.051 dex with fifth iterations , consistent with those from non - differential analysis . as shown in table  [ tab : ks ] , the separation of the mean [ fe / h]@xmath13  values between the @xmath0-@xmath1  and the @xmath0-@xmath2  groups is larger than 2.0  2.5 @xmath47 levels .
we calculated the significance levels to reject the hypothesis that the mean [ fe / h ]  values of the two groups of stars in m22 are identical .
we obtained the significance levels lower than 1% for [ fe / h]@xmath17  and 7.5% for [
fe / h]@xmath13 , strongly suggesting that they are different .
next , we derived the metallicity of individual stars based on the spectroscopic @xmath6  and @xmath7 .
the differences in the mean metallicity between the two groups of stars are @xmath16[fe / h]@xmath17  = 0.203 @xmath22 0.050 dex and @xmath16[fe / h]@xmath13  = 0.204 @xmath22 0.052 dex for non - differential analysis and @xmath16[fe / h]@xmath17  = 0.194 @xmath22 0.044 dex and @xmath16[fe / h]@xmath13  = 0.228 @xmath22 0.053 dex for differential analysis , making the separation of the mean [ fe / h]@xmath13  values between the two groups larger than 3.9 @xmath47 to 4.3 @xmath47 levels .
not surprisingly , our results are consistent with those obtained by @xcite , who relied on the traditional spectroscopic stellar parameters . as shown in table  [
tab : ks ] , the significance levels to reject the hypothesis that the mean [ fe / h ]  values of the two groups of stars in m22 are identical are very low , indicating that they are different .
the metallicity based on the evolutionary stellar parameters also suggest that the metallicity distributions of each group of stars are indeed different .
following the same procedure described in  [ chap : fixedt ] , we obtained the differences in the mean metallicity of @xmath16[fe / h]@xmath17  = 0.224 @xmath22 0.061 dex and @xmath16[fe / h]@xmath13  = 0.168 @xmath22 0.066 dex for non - differential analysis and @xmath16[fe / h]@xmath17  = 0.191 @xmath22 0.043 dex and @xmath16[fe / h]@xmath13  = 0.172 @xmath22 0.060 dex for differential analysis . the mean [ fe / h]@xmath13  values of each group are different more than 2.5 @xmath47 to 2.9 @xmath47 levels and the low significance levels of being identical distributions also confirm that they are different .
similar conclusion can be drawn when we use the evolutionary @xmath7 and @xmath6  at a given @xmath18 , where we used @xmath19  = 14.15 mag for m22 @xcite .
following the same procedure described in  [ chap : fixedvvhb ] , we obtained the differences in the mean metallicity of @xmath16[fe / h]@xmath17  = 0.216 @xmath22 0.047 dex and @xmath16[fe / h]@xmath13  = 0.128 @xmath22 0.062 dex for non - differential analysis and @xmath16[fe / h]@xmath17  = 0.181 @xmath22 0.034 dex and @xmath16[fe / h]@xmath13  = 0.132 @xmath22 0.056 dex for differential analysis . similar to the results shown above , the mean [ fe / h]@xmath13  values of each group are different more than 2.1 @xmath47 to 2.4 @xmath47 levels and the low significance levels of being identical distributions also confirm that they are different .
as shown in figures  [ fig : m22comb ] and [ fig : m22combdiff ] , it should be emphasized that the substructures not only in the [ fe / h]@xmath17  but also in the [ fe / h]@xmath13  distributions are notable , indicating that m22 contains multiple stellar populations with heterogeneous metallicities .
the precision elemental abundance measurement of individual stars in ggs is not a trivial task , especially for a peculiar gc with multiple stellar populations with heterogeneous metallicity distributions . in the context of lte analysis , our demonstrations with a mock peculiar gc composed of the two normal gcs ,
ngc  6752 and m55 , showed that the internal absolute and the relative metallicity scales are vulnerable to the incorrect treatment of the input stellar parameters among multiple stellar populations with different metallicity . in particular , photometric surface gravity without taking care of proper metallicity effect can result in the [ fe / h]@xmath13  measurement error of as large as 0.1  0.2 dex .
as we have discussed earlier , this is because the  line opacity does not vary with surface gravity since the almost all iron atoms are populated in the first ionized level , while the h@xmath5  continuum opacity is sensitively dependent on the electron pressure and , therefore , surface gravity . as a consequence , changes in surface gravity can mimic the [ fe / h]@xmath13  abundance of rgb stars in gcs . in this regard
, we developed methods independent of the traditional spectroscopic analysis approach , which is demanding excitation and ionization equilibria in and elements , to make use of the evolutionary surface gravity .
the metallicity scales from these new approaches are in good agreement with those of previous studies by others @xcite .    from our study of a mock peculiar gc
, it is worth to mention three comments concerning the metallicity scale of the multiple stellar populations in a gc .
first , the use of narrow band photometry , such as @xmath36 and @xmath14 which are sensitive to metallicity and less sensitive to interstellar reddening , is beneficial to discern small [ fe / h ]  differences .
second , our results show that our adaptive methods to estimate appropriate surface gravity would be essential in deriving the absolute and the relative [ fe / h]@xmath13  scale for the multiple stellar populations in peculiar gcs .
third , it is very interesting to note that our new methods appear to provide similar metallicity scale as that from the traditional spectroscopic analysis , suggesting that , the metallicity scale from the widely used traditional spectroscopic approach which makes use of excitation and ionization equilibria in and elements , is valid , at least in the relative sense .
contrary to the conclusion made by mu15 , our re - examination of m22 rgb stars showed that the peculiar gc m22 is composed of two groups of stars with heterogeneous metallicities , confirming our previous results and those of others @xcite , and the m22 saga will continue .
this work has been supported by the center for galaxy evolution research ( grant no .
2010 - 0027910 ) and the basic science research program ( grant no .
2016r1a2b4014741 ) through the national research foundation of korea ( nrf ) funded by the korea government ( msip ) .
lee thanks eugenio carretta for kindly providing ew measurements for ngc  6752 and m55 rgb stars and david yong and donghoh kim for their helpful discussions .
he also appreciates the anonymous referee for critical reading and for providing helpful suggestions that greatly improved the manuscript .
marino , a.  f. , milone , a.  p. , sneden , c. , bergemann , m. , kraft , r.  p. , wallerstein , g. , cassisi , s. , aparicio , a. , asplund , m. , bedin , r.  l. , hilker , m. , lind , k. , momany , y. , piotto , g. , roederer , i.  u. , stetson , p.  b. , & zoccali , m.  2012 , , 541 , 15    marino , a.  f. , sneden , c. , kraft , r.  p. , wallerstein , g. , norris , j.  e. , da costa , g. , milone , a.  p. , ivans , i.  i. , gonzalez , g. , fulbright , j.  p. , hilker , m. , piotto , g. , zoccali , m. , & stetson , p.  b.  2011 , , 532 , a8                      crrrrr [ fe / h]@xmath17 & 0.044 @xmath22 0.003 & @xmath410.044 @xmath22 0.003 & & @xmath410.007 @xmath22 0.002 & 0.000 @xmath22 0.001 + [ fe / h]@xmath13 & @xmath410.037 @xmath22 0.006 & 0.036 @xmath22 0.005 & & 0.090 @xmath22 0.004 & @xmath410.084 @xmath22 0.004 +    crrrrrrrrrr mu15 m2 & 0.232 @xmath22 0.066 & 3.52 & 1.93 & 0.00 & & 0.053 @xmath22 0.022 & 2.43 & 8.98 &
0.00 + mu15 m2 ( 5@xmath61 iter . )
& 0.217 @xmath22 0.076 & 2.84 & 3.02 & 3.14 & & 0.090 @xmath22 0.041 & 2.20 & 7.24 & 5.57 + @xmath6  + @xmath7 ( a99 , no iter . )
& 0.239 @xmath22 0.057 & 4.20 & 0.76 & 0.00 & & 0.096 @xmath22 0.048 & 2.01 & 9.11 & 11.27 + @xmath6  + @xmath7 ( a99 , 5@xmath61 iter . ) & 0.233 @xmath22 0.048 & 4.24 & 0.72 & 0.00 & & 0.108 @xmath22 0.052 & 1.80 & 12.32 & 11.42 + spectroscopic & 0.203 @xmath22
0.050 & 3.48 & 1.32 & 0.00 & & 0.204 @xmath22 0.052 & 3.38 & 1.49 & 3.02 + fixed @xmath6 & 0.224 @xmath22 0.061 & 3.67 & 1.04 & 0.00 & & 0.168 @xmath22
0.066 & 2.54 & 5.44 & 2.85 + fixed @xmath18 & 0.216 @xmath22 0.047 & 4.57 & 0.52 & 0.00 & & 0.128 @xmath22 0.062 & 2.08 & 8.56 & 8.46 + & & & & & & & & & +   + mu15 m2 & 0.213 @xmath22 0.064 & 3.30 & 1.92 & 0.00 & & 0.077 @xmath22
0.033 & 2.35 & 8.01 & 2.89 + mu15 m2 ( 5@xmath61 iter .
) & 0.205 @xmath22 0.069 & 2.96 & 2.60 & 2.71 & & 0.112 @xmath22 0.036 & 3.08 & 3.10 & 0.00 + @xmath6  + @xmath7 ( a99 , no iter . ) & 0.203 @xmath22 0.039 & 5.23 & 0.40 & 0.00 & & 0.101 @xmath22 0.045 & 2.25 & 6.66 & 5.55 + @xmath6  + @xmath7 ( a99 , 5@xmath61 iter . ) & 0.220 @xmath22 0.043 & 4.19 & 0.90 & 0.00 & & 0.129 @xmath22
0.051 & 2.19 & 7.33 & 5.74 + spectroscopic & 0.194 @xmath22 0.044 & 3.79 & 0.92 & 0.00 & & 0.228 @xmath22
0.053 & 3.74 & 1.09 & 0.00 + fixed @xmath6 & 0.191 @xmath22
0.043 & 4.43 & 0.45 & 0.00 & & 0.172 @xmath22 0.060 & 2.86 & 3.42 & 2.86 + fixed @xmath18 & 0.181 @xmath22 0.034 & 5.36 & 0.41 & 0.00 & & 0.132 @xmath22 0.056 & 2.36 & 5.67 & 5.61 +    cccccccccc phot . & c09 ( no ) & @xmath412.087 @xmath22 0.018 & @xmath411.936 @xmath22 0.020 & & @xmath411.645 @xmath22 0.014 & @xmath411.471 @xmath22 0.015 & & 0.443 @xmath22 0.023 & 0.465 @xmath22 0.025 + & c09 ( 5@xmath61 ) & @xmath412.090 @xmath22 0.016 & @xmath411.927 @xmath22 0.028 & & @xmath411.642 @xmath22 0.014 & @xmath411.432 @xmath22 0.024 & & 0.448 @xmath22 0.021 & 0.496 @xmath22 0.037 + & y13 ( no ) & @xmath412.067 @xmath22 0.016 & @xmath411.896 @xmath22 0.021 & & @xmath411.687 @xmath22 0.005 & @xmath411.463 @xmath22 0.006 & & 0.379 @xmath22 0.017 & 0.443 @xmath22 0.022 + & y13 ( 5@xmath61 ) & @xmath412.073 @xmath22 0.014 & @xmath411.869 @xmath22 0.029 & & @xmath411.683 @xmath22 0.005 & @xmath411.414 @xmath22 0.010 & & 0.390 @xmath22 0.015 & 0.454 @xmath22 0.031 + & & & & & & & & & + phot . &
c09 ( no ) & @xmath412.114 @xmath22 0.017 & @xmath411.789 @xmath22 0.022 & & @xmath411.645 @xmath22 0.014 & @xmath411.471 @xmath22 0.015 & & 0.469 @xmath22 0.022 & 0.318 @xmath22 0.027 + & c09 ( 5@xmath61 ) & @xmath412.110 @xmath22 0.016 & @xmath411.883 @xmath22 0.030 & & @xmath411.642 @xmath22 0.014 & @xmath411.432 @xmath22 0.024 & & 0.469 @xmath22 0.021 & 0.451 @xmath22 0.038 + & y13 ( no ) & @xmath412.093 @xmath22 0.016 & @xmath411.748 @xmath22 0.022 & & @xmath411.687 @xmath22 0.005 & @xmath411.463 @xmath22 0.006 & & 0.406 @xmath22 0.017 & 0.285 @xmath22 0.023 + & y13 ( 5@xmath61 ) & @xmath412.092 @xmath22 0.015 & @xmath411.825 @xmath22 0.031 & & @xmath411.683 @xmath22 0.005 & @xmath411.414 @xmath22 0.010 & & 0.409 @xmath22 0.016 & 0.411 @xmath22 0.032 + & & & & & & & & & + spec . &
c09 & @xmath412.003 @xmath22 0.020 & @xmath412.002 @xmath22 0.019 & & @xmath411.583 @xmath22 0.027 & @xmath411.584 @xmath22 0.025 & & 0.420 @xmath22 0.033 & 0.418 @xmath22 0.031 + & y13 & @xmath411.952 @xmath22 0.014 & @xmath411.956 @xmath22 0.017 & & @xmath411.657 @xmath22 0.003 & @xmath411.635 @xmath22 0.004 & & 0.295 @xmath22 0.014 & 0.321 @xmath22 0.017 + & & & & & & & & & + fixed @xmath6 & c09 & @xmath412.001 @xmath22 0.019 & @xmath412.026 @xmath22 0.030 & & @xmath411.579 @xmath22 0.025 & @xmath411.455 @xmath22 0.025 & & 0.422 @xmath22 0.031 & 0.571 @xmath22 0.039 + & y13 & @xmath411.966 @xmath22 0.015 & @xmath411.890 @xmath22 0.028 & & @xmath411.655 @xmath22 0.003 & @xmath411.463 @xmath22 0.010 & & 0.311 @xmath22 0.015 & 0.426 @xmath22 0.030 + & & & & & & & & & + fixed @xmath18 & c09 & @xmath411.999 @xmath22 0.009 & @xmath412.006 @xmath22 0.027 & & @xmath411.651 @xmath22 0.011 & @xmath411.467 @xmath22 0.021 & & 0.348 @xmath22 0.014 & 0.539 @xmath22 0.034 + & y13 & @xmath412.005 @xmath22 0.008 & @xmath411.943 @xmath22 0.029 & & @xmath411.690 @xmath22 0.004 & @xmath411.465 @xmath22 0.010 & & 0.315 @xmath22 0.009 & 0.479 @xmath22 0.031 +