dwarf galaxies are the most common type of galaxy in the universe and are thought to be the progenitors of larger structures .
therefore , a clear understanding of dwarf galaxy evolution is the first step toward a larger understanding of general galaxy evolution .
fortunately , dwarf galaxies are also comparatively simple systems that can be effectively modeled and simulated . in particular , the satellite galaxies of the milky way , by virtue of being nearby , compact and situated in the sparser regions of the galactic halo , have proven to be exceptional laboratories for studying models of stellar and chemical evolution ( see , e.g. mateo 1998 ) .
a useful tool for studying the dwarf galaxies is the comparison of their chemical properties to those of the well - studied milky way globular clusters .
some globular clusters show evidence of multiple stellar populations ( omega centauri ( lee et al .
1999 ) ; ngc  1851 ( milone et al . 2008 ) , ngc  2808 ( piotto et al .
others , particularly those associated with dsph galaxies or tidal streams , show peculiar chemical enrichment histories and/or young ages ( see , e.g. , letarte et al .
however , most _ inner halo _ milky way globular clusters formed their stars a hubble time ago in a single burst ( marin - franch et al .
since the early universe was dominated by massive stars , it is likely that only core collapse supernovae ( type ii sn ) were able to enrich the primordial interstellar medium ( ism ) before the globular clusters formed .
these events would have enriched the ism with both fe and @xmath2-elements , resulting in relatively high @xmath2-element to iron ratios [ @xmath2/fe ] . by contrast , objects that formed later , after type ia sn had time to enrich the ism with more fe , would have lower relative @xmath2-element ratios .
most dwarf galaxies have compound populations , having undergone slow sporadic star formation over their lifetimes ( see review in mateo 1998 ) . while their earliest populations were formed in the aftermath of the the first snii events , their later populations benefited from ongoing snia .
the result is a steep decline in @xmath2-element abundance ratios with decreasing age / increasing metallicity .
the most metal - poor stars have @xmath2-element ratios similar to globular clusters , while the more metal - rich populations are deficient in @xmath2-elements ( shetrone et al .
1998 , 2001 , 2003 ; tolstoy et al . 2003 ; venn et al . 2004 ; smith et al . 2006 ; helmi et al . 2006 ) .
the chemistry of dwarf spheroidal ( dsph ) stars is of particular importance in untangling the hierarchical formation of the milky way .
cdm cosmological models predict that galaxies like the milky way have grown through the merging of low - mass systems with perhaps some low - level merging continuing at the present time . among the current retinue of dsph satellite galaxies ,
the sagittarius dsph shows clear signs of having been torn apart by the milky way potential ( majewski et al . 2003 ) and ursa minor , leo i , coma berenices and ursa major ii also hint at a violent dynamical past ( palma et al .
2003 ; sohn et al . 2007 ;
simon & geha 2007 ) .
however , while the most metal - poor stars of the dsph galaxies are chemically similar to the most metal - poor stars of the galactic halo , the enrichment of snia kicks in at lower fe abundances in the dsph galaxies , resulting in consistent deficits in the @xmath2-element ratios of dsph stars compared to halo stars of similar [ fe / h ] ( venn et al . 2004 ) . this would indicate that the galactic halo formed from objects that had more rapid and therefore more snii - dominated chemical evolution and that objects like the current dsph galaxies have contributed little to the halo .
this argument is , however , undermined by recent studies indicating that the stars the sagittarius dsph is contributing to the halo differ from those in the residual core ( bellazzini et al .
2006 ; chou et al . 2007 ; siegel et al . 2007 ) .
if the dsph galaxies are the remnants of much larger objects that have been disrupted over a hubble time , the _ present _ stellar populations of the dsph galaxies may not reflect the _ initial _ stellar populations they would have preferentially contributed to the galactic halo .
moreover , if the stellar populations of the dsph galaxies have not evolved in isolation , their chemical abundances may have been affected by continual dynamical harassment .
these arguments over dsph stellar chemistry could be clarified by the examination of more distant dwarf galaxies that have not been perturbed by the milky way potential . if an isolated dsph were shown to have the same chemical abundance patterns as the relatively nearby dsphs , such a discovery would ( a ) strengthen the argument that the lack of halo stars with  dsph - like " abundance patterns suggests that the vast majority of the halo formed very quickly and that the dsphs have contributed very few stars since this initial construction phase ; and ( b ) show that the @xmath2-element deficits have nothing to do with dynamical stirring and are entirely the product of dsph stellar and chemical evolution .    in an effort to address these issues
, we now expand the study of chemical abundances to the leo ii dsph galaxy .
leo ii is the second most distant dsph galaxy assumed to be orbiting the milky way ( 218 kpc , bellazzini et al .
2005 ; siegel et al .
2008 , hereafter s08 ) .
it is predominantly metal - poor ( hodge 1982 ; demers & harris 1983 ; aaronson et al .
1983 ; azzopardi et al . 1985 ; suntzeff et al . 1986 ; lee 1995 ; demers & irwin 1993 ; koch et al .
2007a , b ; bosler et al .
2007 , hereafter bss07 ) and dominated by intermediate age populations ( mighell & rich 1996 ; dolphin et al .
2005 ; bellazzini et al .
2005 ; gullieuszik et al .
it has a handful of carbon stars ( azzopardi et al .
1985 ) , a rich population of rr lyrae variable stars ( swope 1967 , 1968 ; siegel & majewski 2000 ) , a high central velocity dispersion ( vogt et al .
1995 ; koch et al . 2007b ; bss07 ; s08 ) and an extended stellar distribution ( coleman et al . 2007 ; s08 ) .
at present , however , the indications are that is had undergone little , if any , interaction with the galactic potential ( koch et .
al 2007b ; coleman et al . 2007 ;
walker et al .
2007 ; s08 ) and its modest radial velocity ( + 26.2 km s@xmath3 , s08 ) would be consistent with a circular orbit that does not bring it close to the milky way .
the bss07 study cited above analyzed medium resolution nir spectra of a large number of leo ii rgb stars in an effort to measure radial velocities and metallicities .
these spectra are dominated by the calcium ii infrared triplet ( cat ) features but also have numerous weak fe , ti , and mg lines that have yet to be exploited . in this paper , we show that these weak lines can be used to measure metallicity and abundance ratios for stars of sufficient s / n through the use of synthetic model spectra .
this method has allowed us to study the chemical evolution of leo ii , to constrain the relationship between the equivalent width of the cat lines and the [ fe / h ] abundances as measured from fei lines and to make a preliminary examination of the age - metallicity relationship ( amr ) of the brightest leo ii rgb stars .
 2 of this paper discusses the additional reduction and analysis we have applied to the bss07 data .  3
compares the abundances ratios measured in leo ii to those measured in the globular clusters to provide insight into the chemical evolution of leo ii .
 4 examines the inferred ages of our bright leo ii stars while  5 summarizes our findings .
a full description of the observations and data reduction can be found in bss07 but we give a brief summary of the observations here .
low - dispersion spectra of red giants in the leo ii dsph were taken using the keck i 10-m telescope and lris ( oke et al . 1995 ) in 2002 and 2003 .
the lris was configured with the 1200 l / mm grating blazed at 7500 which gives a dispersion of 0.62  per pixel and resolution of 1.55 (r @xmath4 10 000 ) .
twenty to thirty stars were observed on each slit mask creating a total sample of 74 stars with snr between 10 and 44 .
traditional equivalent width analysis of these spectra proved impossible due to the heavy blending of many lines and the difficulty in setting the continuum measure .
however , these problems can be overcome by comparison of the observed spectra to synthetic spectra created by the 2007 version of the lte line - analysis code moog ( sneden 1973 ) .
the initial line list was obtained from the r. l. kurucz cd - rom 23 . in order to further calibrate this line list
we compared the list with the nearby red giant star , arcturus ( hinkle et al .
we adopted atmospheric parameters for arcturus (= 4280 , = 1.55 , @xmath5 = 1.61 km s@xmath3 , [ fe / h ] = -0.50 , [ ti / fe ] = + 0.26 , [ mg / fe ] = + 0.39 and [ ca / fe ] = + 0.21 ) by combining the results of mcwilliam ( 1990 ) , fulbright et al .
( 2006,2007 ) and koch & mcwilliam ( 2008 ) . since we are performing a differential analysis , we chose to adjust the gf values in our synthetic spectrum only for lines that were significantly discrepant from the catalog arcturus spectrum , rather than adjust all 8345 lines .
we smoothed the model spectra with a gaussian kernel of fwhm equivalent to the bss07 spectral resolution of 1.55  .
initial metallicity estimates for globular cluster and leo ii stars were taken from bss07 .
effective temperatures were calculated from interpolation of the temperature - color - metallicity relationship of ramirez & melendez ( 2005 ) and the @xmath6 photometry of s08 , the latter allowing the leo ii and cluster data to be analyzed using the same photometric system .
surface gravities were calculated from the temperatures and bolometric - corrected magnitudes assuming a mass of 0.8 @xmath7 .
we used distance moduli of 15.07 , 13.76 , 15.59 , 15.19 and 21.61 for m3 , m107 , ngc  1094 , ngc  4590 and leo ii , respectively .
values by approximately 0.04 , an insignificant correction . ]
the microturbulence was estimated as @xmath8 , in accordance with recent analyzes of globular cluster giants near the tip of the giant branch ( e.g. sneden et al .
2004 , lee , carney & habgood 2005 , ivans et al .
2001 , cohen & melendez 2005 ) .
microturbulent velocities ranged from 1.4 km s@xmath3 to 2.0 km s@xmath3 in the globular cluster sample and from 1.6 to 2.0 km s@xmath3 in the leo ii sample .
these parameters were fed into moog to produce preliminary synthetic spectra .
we then adjusted the [ fe / h ] metallicity of the atmospheres until the fe lines matched the observed spectrum .
a weighted mean of all the iron lines was used to determine the best [ fe / h ] value for each star .
the model was then set to this metallicity and the abundance ratios were determined by holding all elements constant while varying the desired @xmath2-element .
the best fit was determined by the output difference between the moog synthetic spectra and observed spectra .
we found that we were able to effectively measure abundances for stars with snr @xmath0 per pixel .
the setting of the continuum was critical for a proper analysis of the weak lines in these medium resolution spectra . to refine the continuum levels we created a synthetic spectrum for each star , divided it into the observed spectrum and fit the residual with a high order spline .
this fitted spline was then divided into the observed spectrum to remove high order terms in the continuum placement .
if our modeled spectrum changed significantly , e.g. when a new overall metallicity change was required , we repeated this procedure for the new model .
we attempted to determine abundance ratios for all of the moderate lines in the spectra .
these lines include fe , ca , na , mg , and ti .
figure 1 compares one of our highest snr spectra to several synthetic spectra with varied fe and ti abundances .
there are a large number of lines in this region with varying strength that can be used to determine abundances . for a more typical snr spectrum ,
most of these lines are not strong enough for abundance analysis .
we identified sufficient ti and fe lines in all of our bright stars but the na lines and all but one of the mg lines were absent in most of the spectra .
thus we limit our analysis to the stronger fe , ca , mg and ti lines listed in table 1 .
the only ca lines in the observed spectral region are the extremely strong ca triplet lines .
the cat lines are difficult to model with synthetic spectra ( hence the more common use of globular cluster - calibrated equivalent width analysis ) . at a constant ca abundance ,
the strength of the cat lines decreases as the electron pressure increases due to changes in the continuous opacities . for cool stellar atmospheres ,
the main sources of electrons are mg , fe , si , ca , na and al , but which element contributes the most depends in the effective temperature of the star and which layer of the atmosphere is surveyed .
for example , for a = 4200 k star at 0.01 optical depth , the main sources of electrons are mg and fe ; while for a = 4000 k star at 0.01 optical depth , the main sources of electrons are ca , mg , na and al .
thus , as the @xmath2-abundance declines , the continuous opacity changes to counteract the ca abundance decline .
for example , a 0.5 dex change in the electron contributors will increase the cat line strength enough to mimic a ca abundance increase of 0.38 and 0.22 dex for stars of = 4000 , = 0.7 , [ m / h ] = -1.1 and = 4200 , = 0.7 , [ m / h ] = -1.9 , respectively .
this means that as [ m / h ] rises , any decline in the [ ca / fe ] ratio would be muted in the measured cat lines . as a further complication
, the cat lines become less sensitive to the ca abundance as the lines strengthen .    for the globular cluster sample , the plez model atmospheres with [ @xmath2/fe]=+0.4
enhanced abundance ratios were appropriate and there was no need to change the model .
however , our initial analysis showed that the ti and mg ratios of the leo ii stars decline with increasing metallicity (  3 ) .
this resulted in stars with low @xmath2-abundances having erroneously high ca abundances . to derive accurate ca abundance ,
we were therefore required to iterate .
we first determined the fe abundance for each star and then adjusted the model appropriately .
we next determined the ti and mg abundances , from which we calculated a preliminary @xmath2-abundance .
finally , we adjusted the model abundance based on the average @xmath2-abundance at each star s [ fe / h ] .
the result was that for the more metal - rich leo ii stars , we used model atmospheres that were more metal - poor than the fe abundance would imply , thus compensating for the lack of electrons owing to the lower @xmath2-abundances .
measurement errors were determined for each fitted feature by adjusting the abundances up and down until the residuals of the fit were larger than the surrounding continuum regions , ie .
larger than the snr .
the measurement errors were then propagated to yield a single measurement error which was added in quadrature to the abundance errors from the modeling uncertainties
k for , @xmath10 dex for , @xmath11 for v@xmath12 and @xmath10 dex for [ m / h ] . typically the measurement errors dominated over the modeling errors particularly for [ mg / fe ] which was derived from a single spectral feature .    tables 2 and 3 list the globular cluster and leo ii stars from bss07 for which we have been able to determine abundances while tables 4 and 5 show the abundances derived from the moog software .
we have recalculated the signal to noise ratio ( snr ) from the residual of the best fit synthetic spectra to the observed spectra and list that snr in tables 2 and 3 .
figure 2 shows the photometrically derived effective temperatures and surface gravities for the leo ii sample .
there is a clear rgb sequence , corresponding to the photometric sequence upon which the temperatures and gravities are based .
the two outlier stars  195 and 166  are metal - poor radial velocity members .
they are , however , according to s08 , well - removed from the dominant leo ii rgb locus in color - magnitude space although both have washington+ddo51 photometry consistent with giant stars .
it is possible that these are interloping giants in the field , possibly from the metal - poor debris stream of the sagittarius dsph galaxy passing in front of leo ii ( s08 ) .
alternatively , they could represent a young , metal - poor population in leoii (  4 ) .
given their outlier status , we have chosen to remove these two stars from the analysis of mean trends in the leo ii sample .
figure 3 compares our derived [ fe / h ] metallicities against those derived from the cat lines by bss07 . on average ,
the cat - based metallicities are 0.17 dex higher than the our model [ fe / h ] metallicities with an rms of 0.18 dex .
the abundances are correlated but there is a positive residual which shows increasing significance with decreasing metallicity .
the offset is likely due to several factors .
first , all of our fe lines are fei lines and it is possible that metal - poor giants overionize fe ( thevenin & idiart 1999 ; asplund & garcia - perez 2001 ) . in m3 , for example , sneden et al .
( 2004 ) find an underabundance factor of -0.13 in the fei lines and they adopt the feii abundances as the true metallicity . in m68 lee ,
carney & habgood ( 2005 ) find an underabundance factor of -0.37 in the fei lines and also adopt the feii abundances as the true metallicity .
however , not all authors find these overionization factors .
cohen & melendez ( 2005 ) find that the fei and feii abundances agree to within uncertainties .
kraft & ivans ( 2003 ) find small and variable overionization in metal - poor globular clusters but of much smaller magnitude than thevenin & idiart ( 1999 ) predict . in light of the disagreement ,
we have chosen not to make any correction for overionization but will discuss the potential impact in later sections .
a second potential source of the discrepancy between the [ m / h]@xmath13 and our derived [ fe / h ] comes from the calibration of the cat metallicity scale , which in bss07 is based on the carretta & gratton ( 1997 ) metallicity scale .
the choice of calibration system can significantly impact the zero point of the metallicity .
battaglia et al . ( 2008 ) , for example , find a 0.1 dex difference between the carretta & gratton ( 1997 ) metallicity scale and their high resolution analysis over the range of metallicity spanned by leo ii s rgb stars .
we can explore the issue of the metallicity zero point by comparing the fei metallicities we derive for the bss07 globular clusters to the high resolution analyzes listed in pritzl , venn & irwin ( 2005 ) .
the average abundances we derive for ngc  1904 , ngc  4590 and m3 are -0.14 ( @xmath14 ) dex more metal - poor than those listed in pritzl et al .. however , their ngc  1904 abundances come solely from the study of gratton & ortolani ( 1989 ) .
pritzl et al .
find that the gratton & ortolani abundances are typically 0.16 more metal - rich than more recent analyzes .
correcting the ngc  1904 abundance down by 0.16 dex , we find our calculated fe i abundances for ngc  1904 , ngc  4590 and m3 are -0.08 ( 0.14 ) dex more metal - poor than the pritzl et al
. values  a zero point offset within the error of the weighted mean abundances listed in table 4
.    figures 4 and 5 show the derived [ ti / fe ] and [ mg / fe ] abundance ratios , respectively , for both the leo ii and globular cluster sample .
the leo ii and globular cluster star [ ti / fe ] abundances overlap at the poorest metallicities but depart dramatically at higher metallicities .
the leo ii [ mg / fe ] ratios are offset from the globular cluster sample and show a pronounced decline with increasing metallicity .
the large error - bars on the [ mg / fe ] limit the usefulness of the individual measurements but confirm the same trend seen in the [ ti / fe ] abundance ratios .
we note that the [ ti / fe ] ratios in m107 are as high if not higher than those found in the more metal - poor globular clusters .
such large [ ti / fe ] abundance ratios are not unprecedented .
sneden et al .
( 1994 ) found that the relatively metal - rich globular cluster m71 had high ti abundance ratios ( [ ti / fe ] = + 0.5 ) , although ramirez & cohen ( 2002 ) found " normal  ti ratios m71 .
figure 6 shows the [ caii / fei ] abundance ratios for the globular cluster and leo ii samples . there is a slight decline in the [ caii / fei ] ratio with increasing metallicity among the globular cluster sample .
the leo ii ca abundance ratios are offset from the globular cluster sample and the slope of the best fit line is slightly steeper than that of the cluster sample .
this is very similar to the trends in mg abundances seen in figure 5 .
for m3 , ngc  4590 and ngc  1904 , we compare our abundance ratios to the compilation of pritzl , venn & irwin ( 2005 ) .
our [ ti / fe ] ratios are @xmath15 dex larger , our [ mg / fe ] ratios are @xmath16 dex larger , and our [ ca / fe ] ratios are @xmath17 dex larger .
the only significant zero point is in our [ caii / fei ] abundances .
this is not entirely unexpected .
most analyzes compare ionized species to ionized species and neutral species to neutral species .
we would obviously prefer to make a similar comparison , but the previously noted lack of feii lines prevents this .
additionally , as noted above , the literature is divided on the over - ionization correction need to converted from fei to feii abundances .
such corrections are likely to be metallicity dependent and could be the cause of the slope in the globular cluster sample . because we do not make a correction ,
the absolute value of our [ caii / fei ] abundances is uncertain .
however , because our analysis is concerned with _ relative _ differences between the clusters and leo ii , such a zero point correction will have little impact on our conclusions .
figure 7 show the [ ca / h ] ratios derived from our analysis compared to the [ ca / h ] ratios derived from the cat by bss07 .
the overlap of the cluster and leo ii sample in this figure suggests that the corrections we made for the more metal - rich low-@xmath2 stars are consistent with ca abundances derived with the cat methodology .
the slopes in these samples are less steep than those of figure 3  as would be expected since the cat lines should be sensitive to the overall ca abundance .
however , there is a residual slope between the two [ ca / h ] scales , which remains unexplained .
figure 8 shows the mean of the [ ca / fe ] , [ ti / fe ] and [ mg / fe ] ratios .
the fit to the data with the lowest formal reduced chi square is a second order legendre polynomial which exhibits a steep decline in the @xmath2-abundance ratios of the leo ii stars with metallicities above [ fe / h]@xmath4 - 2.0 .
this is consistent with what is seen in other dsph galaxies ( see references in  1 ) and is likely due to slower chemical evolution and concordant increased contribution of snia to leo ii .
however , the data would also be consistent with a step function near [ fe / h]@xmath4 - 1.8 in which the metal - poor stars have globular cluster like @xmath2-abundances while the more metal - rich stars are underabundant in @xmath2 elements .
the slow enrichment implied by our abundance analysis suggests that leo ii should have an age - metallicity relationship ( amr ) in which the more metal - rich stars are younger than the more metal - poor stars and have therefore been subject to more enrichment by snia .
an amr would not be immediately obvious in color - magnitude diagrams due the photometric degeneracy of age and metallicity .
however , our precise spectroscopic survey affords the ability to remove abundance from the equation .
we can now combine precise broadband photometry with modern isochrones to approximate individual stellar ages .
this method has had some success in untangling the stellar populations of the omega centauri globular cluster , for example ( hilker et al .
2004 ; sollima et al . 2005 ) and was recently applied to the leo ii dsph galaxy by koch et al .
( 2007a ) .    the koch et al .
study found that leo ii has undergone steady star formation with little enrichment , starting 9 gyr ago and continuing up until 2 gyr ago .
more recently - formed stars appear to be somewhat enriched .
our analysis uses similar data but with the advantage that we have measured @xmath2-abundances , rather than assumed abundances .    because our stars are near the tip of the rgb , rather than the more age - sensitive turnoff , any analysis will necessarily be imprecise ( koch et al .
estimate uncertainties of 40% using these techniques ) .
the ideal stars for this analysis are those near the age - sensitive main sequence turnoff ( msto ) . however , given that leo ii s msto is well beyond the spectroscopic range of even the largest ground - based telescopes , the rgb is our best option for probing leo ii s age distribution .
figure 9 demonstrates the analysis .
photometry is taken from s08 , which has a precision of 0.01 mag in all passbands at the tip of the red giant branch ( trgb ) .
@xmath18 colors were calculated by converting the washington measures in s08 using transformations in majewski et al .
isochrones are taken from the dartmouth stellar evolution program ( dotter et al .
2007 ) and shifted to a distance modulus of @xmath19 and reddening of @xmath20=0.02 , in accordance with s08 s trgb distance derived from the same data .    cursory examination of the individual panels , broken down by passband and metallicity bin , shows no obvious amr in leo ii . to make a more precise estimate , we interpolated the isochrones in 0.01 magnitude intervals and , with metallicity fixed at the spectroscopic value and [ @xmath21 fixed to the mean [ @xmath22 $ ] at each star s [ fe / h ] , we found the isochrone nearest to each star s photometric position .
we estimated age uncertainties by offsetting the [ fe / h ] of each star by its uncertainty ( which dominates over s08 s small photometric uncertainties ) .
the left panels of figure 10 show the resultant age - metallicity distribution .
although the uncertainties are significant ( typically 2 - 3 gyr ) , some interesting trends can be seen .
the mean age of the stars is @xmath23 gyr but has a spread significantly greater than the uncertainty , implying ongoing low - level star formation .
the younger stars appear to be slightly more metal - rich and indeed , an amr that shows slight enrichment with age  either gradually or step - wise  would be more consistent with the data than a flat trend .
however , we can not statistically rule out a flat amr , which would have a @xmath24 per degree of freedom of 0.9 .
in fact , our derived amr is very similar to of koch et al .
it should be noted that the exploration of the rgb stars has number of complications that could affect the inferred amr . in particular :    @xmath25 leo ii is known to have prominent asymptotic giant branch ( agb ; mighell & rich 1996 ; gullieuszik et al . 2008 ) .
stars on the agb will overlap young rgb isochrones and perhaps 10% of our stars could be agb stars
. the most likely candidates would be the two apparently young metal - poor stars that are clear outliers in the bottom panels of figure 9 and both panels of figure 10 ( note that these are different stars than the two outliers rejected in figures 2 - 8 , which were excluded _ a priori _ from the age analysis ) .
we attempted to measure ages for these stars using the padova isochrones ( marigo et al .
2008 ) , which include an agb . however , these stars are too blue for all but the youngest agb isochrones ( at these magnitudes , the rgb and agb are separated by only a few .01 mag ) . removing these outlier stars from the sample , however , increases the contrast between the younger and older stars of leo ii , making enrichment more tenable .
but a flat amr would still be within the uncertainties ( @xmath26 ) .
@xmath25 any age analysis will be sensitive to the distance modulus of leo ii . adjusting the distance of leo ii shifts the absolute ages , but not the relative ages .
the exception is at short distance moduli , where the stars are too faint to correspond to anything but the maximum age in our isochrone set ( 15 gyr ) . if we assume that leo ii s stars must be younger than 15 gyr , then leo ii must be at a distance of at least @xmath27
. longer distances result in younger absolute ages but have little effect on the inferred amr .
@xmath25 the overionization correction mentioned in  2 could alter the amr by shifting the sample to a more metal - rich distribution .
we applied a 0.2 dex correction ( right panels of figure 10 ) and find that the amr tightens somewhat , producing a sharper peak in the age distribution at around 6 - 7 gyr .
however , a flat amr would still be within the uncertainties ( @xmath26 ) , even if the two outlier stars are removed .
@xmath25 our snr - limited sample is unable to probe the fainter more metal - rich stars in leo ii .
improved spectra of these stars would extend the amr over a greater range of metallicity , providing better constraint of the enrichment and allowing a comparison to toy models of enrichment .
in particular , it would determine if the apparent rise in abundance for the younger stars in figure 10 is merely statistical scatter or is the indicator of a larger trend .    even with these caveats , the age - metallicity distribution of leo ii stars
, while giving tantalizing hints of enrichment , remain statistically consistent with no enrichment
. a flat enrichment would be unexpected for leo ii , since we _ do _ find that the more metal - rich stars have lower @xmath2-element abundances , suggesting gradual enrichment by snia over the course of several gyr .    in the end , while this exercise demonstrates that precise photometry and spectroscopy can be used to probe the amr of objects out to 200 kpc , the enrichment of leo ii may be too subtle to be detected in our data .
spectroscopic examination of the more numerous and age - sensitive turnoff stars may be the only way to tease out the amr of this distant and enigmatic galaxy .
our analysis of the weaker spectral features in the bss07 spectra provides critical information about not only the leo ii dsph galaxy but other spectroscopic studies as well .
@xmath25 we find the fe - line metallicities disagree with cat - based metallicities , although it is likely this is a result of the overionization of fe at low metallicity .
we are unable to analyze the fainter more metal - rich stars to confirm the high - metallicity tail identified by bss07 .
@xmath25 we confirm that leo ii , like the other low - luminosity dsph galaxies , has globular cluster - like @xmath2-abundances in its most metal - poor stars but shows a declining @xmath2-abundance at higher metallicities , likely due to the increasing influence of snia .
@xmath25 leo ii s amr gives tantalizing hints of slow enrichment .
however , a flat amr with no enrichment can not be ruled out with the present data , given the inherent uncertainties in analyzing only the brightest stars in the galaxy .
nevertheless , this demonstrates that precise spectroscopic and photometric data can be used to get a general picture of the enrichment history of distant dsph galaxies .
leo ii s particular influence on models of dwarf galaxy evolution lies in its difference from the other dsph galaxies .
as noted in  1 , leo ii appears to have had little interaction with the milky way .
this would indicate that the @xmath2-abundance patterns found in the dsph galaxies are a universal feature of dwarf galaxies and have no relation to any dynamical interaction with a larger parent galaxy .
it also strengthens the contention that while objects like the present retinue of dsphs have contributed somewhat to the construction of the milky way  given the similarities between their most metal - poor stars and some chemically peculiar stars within the halo  the stark chemical dissimilarities between the bulk of the dsph and field halo stars indicates that this contribution has been small .
further investigation is needed into the zero point issues intrinsic in our analysis so that we may more confidently translate measured fei abundances into intrinsic feii abundances .
this could be accomplished by looking at existing high snr cat spectra of globular clusters , especially those known to be deficient in @xmath2-elements or those known or suspected to have been stripped from larger dsph galaxies .
once the zero point issues are understood , the technique we have unveiled in this paper can be applied to the large cat spectral datasets in the literature to derive @xmath2-abundance measures for thousands of dsph and field halo stars .
better and easier chemical fingerprinting of distant rgb stars has the potential to finally unravel the hierarchical history of the milky way .
moreover , this technique expands the range of current and future investigations into dsph chemistry .
it shows that current telescopes can probe dsph chemistry out to at least 200 kpc and indicates that next - generation telescopes will open up investigation into even more isolated dwarfs such as the phoenix dirr / dsph ( tobolewski et al .
, in prep . )
this will expand the study of galactic chemical evolution into an even more diverse array of galactic environments . at this time , however , the existing models of dwarf galaxy chemical evolution appear to be valid out to at least 200 kpc from the galaxy .
aaronson , m. , olszewski , e. w. & hodge , p. w. 1983 , , 267 , 271 asplund , m. , & garca prez , a.  e.  2001 , , 372 , 601 azzopardi , m. , lequeux , j. & westerlund , b. e. 1985 , , 144 , 388 battaglia , g. , irwin , m. , tolstoy , e. , hill , v. , helmi , a. , letarte , b. , jablonka , p. 2008
, , 383 , 183 bellazzini , m. , gennari , n. & ferraro , f. r. , 2005 , , 360 , 185 bellazzini , m. , correnti , m. , ferraro , f. r. , monaco , l. , & montegriffo , p. 2006a
, , 446 , 1 bosler , t.l .
, smecker - hane , t.a . , stetson , p.b .
2007 , , 378 , 318 ( bss07 ) carretta , e. & gratton , r.g .
1997 , , 121 , 95 chou , m .- y . , et al .
2007 , , 670 , 346 cohen , j.  g. , & melndez , j.  2005 , , 129 , 303 coleman , m. , jordi , k. , rix , h .- w . , grebel , e. k. & koch , a. , 2007 , , 134 , 1938 demers , s. & harris , w. e. 1983 , , 88 , 329 demers , s. & irwin , m. j. 1993 , , 261 , 657 dolphin , a. e. , weisz , d. r. , skillman , e. d. & holtzman , j. a. , 2005 , _ astro - ph/0506430 _ fulbright , j.  p. , mcwilliam , a. , & rich , r.  m.  2006 , , 636 , 821 fulbright , j.  p. , mcwilliam , a. , & rich , r.  m.  2007 , , 661 , 1152 gullieuszik , m. , held , e. v. , rizzi , l. , girardi , k. , marigo , p. & momany , y. , 2008 , , _ in press , astro - ph/0805.0735 _ gratton , r.  g. , & ortolani , s.  1989 , , 211 , 41 helmi , a. , et al .
2006 , , 651 , l121 hilker , m. , kayser , a. , richtler , t. , & willemsen , p.  2004
, , 422 , l9 hinkle , k. , wallace , l. , valenti , j. , & harmer , dianne 2000 , visible and near infrared atlas of the arcturus spectrum 37279300  , ( san francisco : asp ) .
hodge , p. w. , 1982 , , 87 , 1668 ivans , i.  i. , kraft , r.  p. , sneden , c. , smith , g.  h. , rich , r.  m. , & shetrone , m.  2001 , , 122 , 1438 koch , a. , grebel , e.  k. , kleyna , j.  t. , wilkinson , m.  i. , harbeck , d.  r. , gilmore , g.  f. , wyse , r.  f.  g. , & evans , n.  w.  2007a , , 133 , 270 koch , a. , grebel , e. k. , kleyna , j. t. , wilkinson , m. i. , harbeck , d. r. , gilmore , g. f. , wyse , r. f. g. & evans , n. w. 2007b , , 133 , 270 koch , a. , & mcwilliam , a. 2008 , 135 , 1551 kraft , r.  p. , & ivans , i.  i.  2003 , , 115 , 143 lee , m. g. , 1995 , , 110 , 1155 lee , y .- w . , joo , j .-
m . , sohn , y .-
j . , rey , s .- c . ,
lee , h .- c . , & walker , a.  r.  1999 , , 402 , 55 lee , j .- w .
, carney , b.  w. , & habgood , m.  j.  2005 , , 129 , 251 letarte , b. , hill , v. , jablonka , p. , tolstoy , e. , franois , p. , & meylan , g.  2006 , , 453 , 547 majewski , s. r. , ostheimer , j. c. , kunkel , w. e. , & patterson , r. j. , 2000 , , 120 , 2550 majewski , s. r. , skrutski , m. f. , weinberg , m. d. & ostheimer , j. c. , 2003 , 599 , 1082 marigo , p. , girardi , l. , bressan , a. , groenewegen , m.  a.  t. , silva , l. , & granato , g.  l.  2008 , , 482 , 883 marin - franch , a. , aparicio , a. , piotto , g. , rosenberg , a. , chaboyer , b. , sarajedini , a. , siegel , m. , anderson , j. , bedin , l. r. , dotter , a. , hempel , m. , king , i. , majewski , s. , milone , a. , & reid , i. n. , , _ submitted _ mateo , m.  l.  1998 , , 36 , 435 mighell , k. j. & rich , m. r. 1996 , , 111 , 777 milone , a.  p. , et al .
2008 , , 673 , 241 oke , j.  b. , et al .
1995 , , 107 , 375 palma , c. , majewski , s. r. , siegel , m. h. , patterson , r. j. , ostheimer , j. c. & link , r. , 2003 , 125 , 1352 piotto , g. , et al .
2007 , , 661 , l53 pritzl , b.  j. , venn , k.  a. , & irwin , m.  2005 , , 130 , 2140 ramrez , s.  v. , & cohen , j.  g.  2002 , , 123 , 3277 ramirez , i. & melendez , j. 2005 , , 626 , 465 shetrone , m.  d. , bolte , m. , & stetson , p.  b.  1998 , , 115 , 1888 shetrone , m.  d. , ct , p. , & sargent , w.  l.  w.  2001 , , 548 , 592 shetrone , m. , venn , k.  a. , tolstoy , e. , primas , f. , hill , v. , & kaufer , a.  2003 , , 125 , 684 siegel , m. h. & majewski , s. r. , 2000 , , 120 , 284 siegel , m.  h. , et al .
2007 , , 667 , l57 siegel , m. h. , majewski , s. r. , sohn , s. t. , shetrone , m.d . , munoz , r. r. , & patterson , r. j. 2008 , , _ submitted _ [ s08 ] simon , j. d. & geha , m. , 2007 , , 670 , 313 smith , g.  h. , siegel , m.  h. , shetrone , m.  d. , & winnick , r.  2006 , , 118 , 1361 sneden , c. , kraft , r.p . , langer , g.e .
, prosser , c.f .
, & shetrone , m.d .
1994 , , 107 , 1773 .
sneden , c. , kraft , r.  p. , guhathakurta , p. , peterson , r.  c. , & fulbright , j.  p.  2004
, , 127 , 2162 sohn , s. t. , majewski , s. r. , munoz , r. m. , kunkel , w. e. , johnston , k. v. , ostheimer , j. c. , guhathakurta , p. , patterson , r. j. , siegel , m. h. & cooper , m. , 2007 , , 663 , 960 sollima , a. , pancino , e. , ferraro , f.  r. , bellazzini , m. , straniero , o. , & pasquini , l.  2005 , , 634 , 332 suntzeff , n. b. , aaronson , m. , olszewski , e. w. & cook , k.h .
1986 , , 91 , 1091 swope , h. h. 1967 , , 79 , 439 swope , h. h. 1968 , , 73 , s204 thvenin , f. , & idiart , t.  p.  1999 , , 521 , 753 tolstoy , e. , venn , k.  a. , shetrone , m. , primas , f. , hill , v. , kaufer , a. , & szeifert , t.   2003 , , 125 , 707 venn , k.  a. , irwin , m. , shetrone , m.  d. , tout , c.  a. , hill , v. , & tolstoy , e.  2004 , , 128 , 1177 vogt , s. s. , mateo , m , olszewski , e. w. & keane , m. j. , 1995 , 109 , 151 walker , m. g. , mateo , m. , olszewski , e. w. , gnedin , o. y. , wang , x. , sen , b. & woodroofe , m. , 2007 , , _ in press , astro - ph/0708.0010 _ mcwilliam , a. 1990 , 74 , 1075    lrrrrrrr tii & 8378.25 & 3.72 & -2.30 + tii & 8382.53 & 0.82 & -1.63
+ tii & 8435.65 & 0.84 & -1.30 + fei & 8220.38 & 4.32 & 0.20 + fei & 8327.06 & 2.20 & -1.64 + fei & 8387.77 & 2.18 & -1.60 + fei & 8468.41 & 2.22 & -2.17 + fei & 8514.07 & 2.20 & -2.25 + fei & 8611.80 & 2.85 & -2.15 + fei & 8674.75 & 2.83 & -2.00 + fei & 8688.63 & 2.18 & -1.21 + mgi & 8806.76 & 4.35 & -0.13 + ca ii & 8498.02 & 1.69 & -1.31 + ca ii & 8542.09 & 1.70 & -0.36 + ca ii & 8662.14 & 1.69 & -0.62 +    lrrrrrr + 241 & 13.61 & 1.11 & 37.6 & 4419 & 0.94 & 1.77 + 131 & 13.02 & 1.20 & 34.9 & 4298 & 0.85 & 1.80 + 223 & 13.19 & 1.16 & 34.5 & 4350 & 0.72 & 1.85 + 160 & 14.06 & 1.52 & 31.2 & 3969 & 0.31 & 2.02 + 294 & 14.22 & 1.00 & 22.3 & 4586 & 1.29 & 1.62 + 181 & 13.95 & 1.25 & 21.3 & 4236 & 0.94 & 1.76 +   + hi82 & 12.59 & 1.33 & 24.6 & 4188 & 0.34 & 1.98 + hi119&13.62 & 0.96 & 21.1 & 4633 & 1.09 & 1.70 + hi239&14.19 & 0.87 & 14.2 & 4788 & 0.89 & 1.79 +   + 265 & 13.26 & 1.30 & 93.0 & 4155 & 0.83 & 1.81 + 250 & 14.11 & 0.95 & 82.1 & 4627 & 1.50 & 1.54 + 640 & 13.27 & 1.26 & 77.1 & 4199 & 0.87 & 1.79 + 589 & 12.87 & 1.33 & 76.6 & 4124 & 0.64 & 1.89 + 1217 & 14.00 & 1.05 & 76.4 & 4470 & 1.36 & 1.59 + 885 & 13.47 & 1.07 & 75.5 & 4441 & 1.13 & 1.69 + 334 & 13.24 & 1.20 & 55.3 & 4269 & 0.91 & 1.78 + 238 & 12.69 & 1.57 & 53.5 & 3918 & 0.35 & 2.01 +   + sl & 14.04 & 1.47 & 56.1 & 4581 & 1.37 & 1.59 + sf & 13.39 & 1.70 & 49.0 & 4225 & 0.87 & 1.79 + sr & 14.66 & 1.29 & 42.2 & 4910 & 1.80 & 1.41 + sh & 13.84 & 1.61 & 38.9 & 4356 & 1.15 & 1.68 + ss & 14.79 & 1.40 & 37.9 & 4703 & 1.74 & 1.43 + su & 14.78 & 1.28 & 37.5 & 4930 & 1.86 & 1.39 + s62 & 13.97 & 1.62 & 34.1 & 4341 & 1.19 & 1.66 +    lrrrcrrr 180 & 18.998 & 1.5117 & 43.7 & -1.60 & 3981 & 0.33 & 2.02 + 255 & 19.291 & 1.3114 & 40.0 & -1.66 & 4151 & 0.62 & 1.90 + 271 & 19.375 & 1.3561 & 39.5 & -1.48 & 4100 & 0.60 & 1.90 + 258 & 19.311 & 1.3507 & 36.9 & -1.57 & 4110 & 0.59 & 1.91 + 336 & 19.531 & 1.1669 & 36.8 & -2.03 & 4299 & 0.83 & 1.81 + 195 & 19.026 & 1.0915 & 36.2 & -1.60 & 4406 & 0.71 & 1.86 + 296 & 19.385 & 1.0488 & 35.7 & -2.02 & 4438 & 0.87 & 1.79 + 293 & 19.304 & 1.2481 & 35.5 & -1.84 & 4216 & 0.68 & 1.87 + 166 & 18.827 & 1.1629 & 35.3 & -1.52 & 4318 & 0.56 & 1.92 + 285 & 19.337 & 1.3657 & 35.0 & -1.50 & 4091 & 0.58 & 1.91 + 304 & 19.389 & 1.2711 & 34.7 & -1.81 & 4194 & 0.69 & 1.87 + 254 & 19.259 & 1.3981 & 33.6 & -1.43 & 4055 & 0.51 & 1.87 + 236 & 19.212 & 1.3330 & 32.6 & -1.90 & 4144 & 0.58 & 1.91 + 351 & 19.552 & 1.2835 & 32.5 & -1.40 & 4176 & 0.74 & 1.85 + 333 & 19.495 & 1.2278 & 32.4 & -1.45 & 4240 & 0.77 & 1.83 + 379 & 19.618 & 1.3171 & 32.1 & -1.52 & 4141 & 0.74 & 1.85 + 209 & 19.033 & 1.4749 & 31.7 & -1.59 & 4006 & 0.37 & 2.00 + 256 & 19.332 & 1.2784 & 30.4 & -1.54 & 4182 & 0.66 & 1.88 + 341 & 19.541 & 1.2780 & 29.2 & -1.43 & 4182 & 0.74 & 1.85 + 281 & 19.342 & 1.1127 & 28.6 & -2.05 & 4358 & 0.80 & 1.82 + 282 & 19.329 & 1.4611 & 27.6 & -1.52 & 4008 & 0.49 & 1.95 + 420 & 19.717 & 1.0824 & 27.1 & -1.46 & 4430 & 1.00 & 1.74 + 377 & 19.587 & 1.1853 & 25.8 & -1.73 & 4283 & 0.84 & 1.80 + 260 & 19.268 & 1.2811 & 25.7 & -1.57 & 4179 & 0.63 & 1.89 + 248 & 19.354 & 1.2399 & 25.1 & -1.51 & 4225 & 0.70 & 1.86 + 234 & 19.250 & 1.2238 & 24.5 & -1.80 & 4240 & 0.67 & 1.87 + 230 & 19.194 & 1.4458 & 24.1 & -1.56 & 4025 & 0.45 & 1.96 +    lrrrr + 241 & -1.70 ( 0.18 ) & 0.11 ( 0.27 ) & 0.30 ( 0.38 ) & 0.70 ( 0.17 ) + 131 & -1.89 ( 0.21 ) & 0.70 ( 0.25 ) & 0.45 ( 0.37 ) & 0.70 ( 0.20 ) + 223 & -1.55 ( 0.18 ) & 0.57 ( 0.22 ) & 0.50 ( 0.39 ) & 0.70 ( 0.17 ) + 160 & -2.06 ( 0.26 ) & -0.03 ( 0.30 ) & 0.70 ( 0.59 ) & 0.35 ( 0.20 ) + 294 & -1.62 ( 0.22 ) & 0.65 ( 0.27 ) & 0.60 ( 0.42 ) & 0.45 ( 0.17 ) + 181 & -2.02 ( 0.22 ) & 0.46 ( 0.25 ) & 0.20 ( 0.43 ) & 0.85 ( 0.24 ) + average & -1.77 ( 0.08 ) & 0.44 ( 0.10 ) & 0.44 ( 0.17 ) & 0.61 ( 0.08 ) +   + hi82 & -2.44 ( 0.16 ) & 0.13 ( 0.20 ) & 0.50 ( 0.30 ) & 0.50 ( 0.22 ) + hi119 & -2.25 ( 0.14 ) & 0.22 ( 0.20 ) & 0.70 ( 0.24 ) & 0.35 ( 0.27 ) + hi239 & -2.03 ( 0.18 ) & 0.28 ( 0.37 ) & -0.20 ( 0.35 ) & 0.30 ( 0.22 ) + average & -2.26 ( 0.09 ) & 0.19 ( 0.13 ) & 0.44 ( 0.17 ) & 0.39 ( 0.13 ) +   + 265 & -1.62 ( 0.14 ) & -0.08 ( 0.16 ) & 0.20 ( 0.25 ) & 0.60 ( 0.17 ) + 238 & -2.10 ( 0.19 ) & 0.30 ( 0.18 ) & 0.10 ( 0.32 ) & 0.45 ( 0.19 ) + 250 & -1.40 ( 0.14 ) & 0.38 ( 0.16 ) & 0.50 ( 0.28 ) & 0.40 ( 0.15 ) + 334 & -1.71 ( 0.14 ) & 0.36 ( 0.16 ) & 0.30 ( 0.22 ) & 0.75 ( 0.15 ) + 589 & -1.77 ( 0.19 ) & 0.43 ( 0.16 ) & 0.30 ( 0.27 ) & 0.65 ( 0.15 ) + 640 & -1.51 ( 0.14 ) & -0.02 ( 0.16 ) & 0.30 ( 0.28 ) & 0.55 ( 0.15 ) + 885 & -1.61 ( 0.15 ) & 0.40 ( 0.17 ) & 0.05 ( 0.34 ) & 0.45 ( 0.16 ) + 1217 & -1.68 ( 0.16 ) & 0.60 ( 0.20 ) & 0.10 ( 0.23 ) & 0.50 ( 0.18 ) + average & -1.64 ( 0.05 ) & 0.28 ( 0.06 ) & 0.24 ( 0.09 ) & 0.55 ( 0.06 ) +   + s62 & -1.03 ( 0.19 ) & 0.55 ( 0.23 ) & 0.50 ( 0.27 ) & 0.50 ( 0.20 ) + sf & -0.94 ( 0.18 ) & 0.60 ( 0.23 ) & 0.15 ( 0.26 ) & 0.40 ( 0.20 ) + sh & -1.18 ( 0.16 ) & 0.72 ( 0.20 ) & 0.35 ( 0.24 ) & 0.50 ( 0.17 ) + sl & -1.17 ( 0.21 ) & 0.68 ( 0.26 ) & 0.35 ( 0.33 ) & 0.50 ( 0.22 ) + sr & -1.01 ( 0.20 ) & 0.55 ( 0.24 ) & 0.25 ( 0.27 ) & 0.35 ( 0.21 ) + ss & -1.05 ( 0.19 ) & 0.89 ( 0.27 ) & 0.10 ( 0.29 ) & 0.40 ( 0.20 ) + su & -0.92 ( 0.23 ) & 0.13 ( 0.40 ) & -0.10 ( 0.33 ) & 0.15 ( 0.24 ) + average & -1.05 ( 0.07 ) & 0.63 ( 0.09 ) & 0.25 ( 0.11 ) & 0.41 ( 0.08 ) +    lrrrr 180 & -1.94 ( 0.18 ) & -0.16 ( 0.16 ) & & 0.25 ( 0.15 ) + 209 & -1.97 ( 0.17 ) & -0.13 ( 0.16 ) & -0.20 ( 0.30 ) & 0.30 ( 0.15 ) + 271 & -1.77 ( 0.18 ) & -0.10 ( 0.18 ) & 0.20 ( 0.30 ) & 0.50 ( 0.17 ) + 351 & -1.73 ( 0.22 ) & -0.17 ( 0.21 ) & 0.25 ( 0.28 ) & 0.30 ( 0.18 ) + 166 & -1.85 ( 0.18 ) & 0.18 ( 0.19 ) & 0.10 ( 0.30 ) & 0.60 ( 0.16 ) + 293 & -1.93 ( 0.18 ) & -0.05 ( 0.17 ) & -0.20 ( 0.15 ) & 0.35 ( 0.15 ) + 336 & -1.98 ( 0.21 ) & & 0.10 ( 0.32 ) & 0.20 ( 0.17 ) + 236 & -2.26 ( 0.22 ) & 0.33 ( 0.21 ) & 0.50 ( 0.33 ) & 0.40 ( 0.20 ) + 195 & -1.70 ( 0.19 ) & 0.36 ( 0.23 ) & 0.70 ( 0.30 ) & 0.42 ( 0.18 ) + 258 & -1.64 ( 0.19 ) & -0.17 ( 0.19 ) & 0.05 ( 0.30 ) & 0.30 ( 0.17 ) + 254 & -1.80 ( 0.19 ) & -0.10 ( 0.18 ) & 0.40 ( 0.27 ) & 0.20 ( 0.17 ) + 255 & -2.05 ( 0.19 ) & 0.30 ( 0.17 ) & -0.35 ( 0.39 ) & 0.45 ( 0.16 ) + 296 & -2.13 ( 0.19 ) & 0.28 ( 0.22 ) & 0.30 ( 0.33 ) & 0.40 ( 0.17 ) + 304 & -2.09 ( 0.18 ) & -0.02 ( 0.18 ) & 0.45 ( 0.27 ) & 0.35 ( 0.16 ) + 281 & -2.07 ( 0.21 ) & 0.35 ( 0.26 ) & 0.15 ( 0.36 ) & 0.35 ( 0.20 ) + 282 & -1.68 ( 0.19 ) & -0.34 ( 0.19 ) & -0.30 ( 0.37 ) & 0.05 ( 0.17 ) + 333 & -1.95 ( 0.21 ) & 0.28 ( 0.21 ) & 0.35 ( 0.30 ) & 0.50 ( 0.19 ) + 285 & -1.73 ( 0.18 ) & -0.03 ( 0.20 ) & 0.05 ( 0.28 ) & 0.45 ( 0.16 ) + 379 & -1.71 ( 0.21 ) & -0.03 ( 0.21 ) & 0.15 ( 0.33 ) & 0.30 ( 0.20 ) + 260 & -1.75 ( 0.21 ) & -0.17 ( 0.20 ) & & 0.35 ( 0.18 ) + 341 & -1.42 ( 0.19 ) & -0.18 ( 0.19 ) & -0.00 ( 0.30 ) & 0.15 ( 0.17 ) + 256 & -1.73 ( 0.18 ) & 0.05 ( 0.18 ) & -0.20 ( 0.16 ) & 0.35 ( 0.16 ) + 248 & -1.72 ( 0.20 ) & 0.03 ( 0.21 ) & 0.05 ( 0.18 ) & 0.40 ( 0.19 ) + 420 & -1.80 ( 0.22 ) & 0.24 ( 0.25 ) & 0.45 ( 0.32 ) & 0.45 ( 0.21 ) + 230 & -1.88 ( 0.24 ) & 0.07 ( 0.25 ) & 0.75 ( 0.34 ) & 0.50 ( 0.23 ) + 377 & -1.98 ( 0.20 ) & 0.09 ( 0.24 ) & 0.15 ( 0.38 ) & 0.40 ( 0.19 ) + 234 & -1.76 ( 0.19 ) & -0.20 ( 0.21 ) & -0.10 ( 0.33 ) & 0.25 ( 0.18 ) +


gas - phase and gas - grain reactions are considered the two main paths by which molecules form in the interstellar medium ( ism ) @xcite . as of february 2016 , almost 200 molecular species
have been observed in the ism and circumstellar shells ; a number that will only increase as observational techniques improve .
it is well known that the observed abundance of molecular hydrogen ( h@xmath0 ) can not be accounted for by gas - phase reactions @xcite .
thus , gas - grain reactions are invoked to account for the amount of h@xmath0 observed in the ism .
the main regions of focus for h@xmath0 formation are molecular clouds ( mcs ) and the cold neutral medium ( cnm ) .
these regions have temperatures that range from 10 k to 100 k , and make up a bit under half of the ism by mass @xcite .
the interplay of dense mcs with the cnm allows for the formation of h@xmath0 and plays a crucial role in star formation .
there are four main steps in the typical gas - grain reaction .
the reactants must first adsorb to the interstellar dust surface ; they then must diffuse across the substrate .
when reactants encounter each other , they must react to form a molecular species .
finally , the newly formed molecule must desorb from the dust surface .    in dense clouds , these interstellar grains / substrates consist of silicate or carbonaceous cores covered by a thin layer of ice about 0.1 @xmath1 m thick @xcite .
the icy layers are believed to predominantly be composed of h@xmath0o ; although ices such as co , co@xmath0 , nh@xmath2 , and ch@xmath3 certainly occur as well @xcite .
the h@xmath0o ice is generally regarded to be amorphous in structure @xcite .    the formation of h@xmath0 on dust grain surfaces is thought to occur through three mechanisms .
the first , dubbed the eley - rideal ( er ) mechanism , involves an impinging hydrogen atom directly hitting an already adsorbed hydrogen atom on the icy surface @xcite . in the langmuir - hinshelwood ( lh ) mechanism ,
two hydrogen atoms already adsorbed to the grain surface diffuse some distance until they meet one other and react to form h@xmath0 @xcite .
the distinction between these two mechanisms lies in the thermalization of the atom with respect to the grain . in the er mechanism
, the adsorbed atom has thermalized to the grain , and is directly hit by an impinging atom that has not yet thermalized .
by contrast , in the lh mechanism , both atoms have already thermalized with the grain . the final mechanism , known as the hot - atom mechanism , involves a hydrogen atom with a high kinetic energy striking the substrate and thermally hopping along the substrate surface until it strikes another already adsorbed hydrogen atom and forms h@xmath0 @xcite .
the common feature of these three mechanisms is the adsorption of at least one of the hydrogen atoms to the icy substrate .
a key quantity characterizing the adsorption properties of a given gas species on a given substrate is the probability that the gas molecule sticks to the substrate when the two collide . in computational studies ,
the sticking probability can be deduced as the fraction of incident gas molecules that stick to the substrate , in the limit of a large number of trials .
the sticking probability for a general species is @xmath4 where @xmath5 is the number of sticking events and @xmath6 is the number of trials .
the sticking probability depends on the temperature of the dust or ice grain ( @xmath7 ) and the incident kinetic energy ( @xmath8 ) of the colliding species when the species and surface are very far apart .    in many experimental and observational situations
, one does nt have access to the incident kinetic energy of the gaseous species , but rather the temperature of the gas ( @xmath9 ) .
the quantity characterizing the sticking in this case is called the sticking coefficient , @xmath10 , a function of both gas and grain temperatures .
the sticking coefficient can be calculated by averaging the sticking probability weighted by a maxwellian velocity distribution of the desired gas temperature .
hollenbach & salpeter ( 1970 ) derived a theoretical form for this sticking coefficient : @xmath11 where @xmath12 , @xmath13 , and @xmath14 .
the variable @xmath15 is the well depth of the adsorption potential , @xmath16 is the energy transferred in a single collision , and @xmath17 is close to unity for lambert s law and equals two in the case of isotropic scattering @xcite . in astrophysical models , this equation is often approximated with @xmath18 where @xmath19 is dependent on the incident species and the type of substrate involved @xcite .
some experiments ( e.g. manico et al .
2001 ) suggest that eqn .
( 3 ) underestimates the sticking coefficient .
the sticking coefficient is of great importance to astrophysical models , as it is a crucial component for the rate equation of various grain - mediated processes ( including h@xmath0 formation ) in the ism .
the goal of this work is therefore to obtain accurate sticking coefficients of atomic hydrogen sticking on amorphous water - ice for astrophysical modeling purposes .
there have been many such studies , experimental and theoretical , with the same goal , though h@xmath0 is primarily formed on warm amorphous carbon and silicate grains @xcite . to our knowledge
, the first simulational study was performed by buch & zhang ( 1991 ) , who studied the sticking of h and d atoms on an amorphous cluster of 115 water molecules using classical molecular dynamics ( md ) .
they found sticking probabilities and coefficients on 10 k ice for both h and d atoms .
a group of studies investigated the sticking and mobility of h atoms on amorphous water - ice slabs consisting of 1000 water molecules using classical md @xcite .
they performed simulations on 10 k and 70 k ice and showed that the sticking probability can not be assumed to be equal to unity for the range of gas temperatures encountered in the ism .
however , due to the computational restraints of that era , not many simulations were performed for each incident kinetic energy of the hydrogen atom .
al - halabi et al . ( 2002 , 2007 ) also performed classical md simulations on both crystalline and amorphous water - ice substrates .
their simulations confirmed that sticking depends on the incident kinetic energy of the h - atom and the temperature of the ice surface by carrying out a large number of trials .
veeraghattam et al .
( 2014 ) performed large numbers of classical md simulations on both 10 k and 70 k amorphous water - ice slabs consisting of 1000 water molecules .
they computed sticking probabilities and coefficients for their simulations that are reproduced and expanded upon in the current work .
experimentally , manico et al .
( 2003 ) , hornekaer et al .
( 2003 ) , and watanabe et al .
( 2010 ) all observed h@xmath0 formation on lab - generated substrates and measured the recombination energy .    in this study
we perform an extensive analysis of atomic hydrogen sticking to amorphous water - ice surfaces for a wide range of grain temperatures ( 10 k - 70 k ) and for a wide range of atomic hydrogen kinetic energies ( 5 k - 400 k ) .
therefore a key goal is to fill in the large gaps in predicted sticking coefficients from previous studies to provide richer , more in - depth information for astrophysical modelers .
here we also pay particularly close attention to the definition of what constitutes a sticking event .
the previous simulational studies mentioned above do not agree on a objective definition for a sticking event , and we have found that different criteria can lead to different results .
the appendix provides additional detail on the subject .
using classical molecular dynamics ( md ) , we study the interaction of a hydrogen atom with an amorphous water - ice slab @xcite . it is reasonable to question whether classical dynamics is appropriate for describing hydrogen
. studies of hydrogen in other contexts have shown that quantum dynamics is preferable for hydrogen , but at a very high computational cost . in the meantime , classical trajectory calculations have been shown to prove useful in describing physical processes involving hydrogen on surfaces @xcite .
but , as andersson et al .
( 2006 ) mention , quantum effects can be tricky to pin down and could be important as the kinetic energy of the hydrogen atom decreases .
nevertheless , at this stage the vastly more complicated and expensive quantum calculations seem unwarranted .
the present simulations require two pairwise potentials : ( 1 ) h@xmath0o - h@xmath0o and ( 2 ) h - h@xmath0o .
these potentials are differentiated to define a force field , which is then used to numerically integrate the coupled , classical equations of motion for the many - particle system using the leapfrog algorithm at a fixed time - step of @xmath20 fs . in non - spherical systems containing rigid molecules , like h@xmath0o , the rotational equations of motion
must also be solved @xcite .
the h@xmath0o - h@xmath0o potential interaction is used to simulate the water - ice substrate while the h - h@xmath0o potential is used to model the hydrogen interaction with that substrate .
these potentials are derived from quantum chemical calculations for a single pair of particles .
our system is treated either as a canonical ensemble ( nve ) or as a isobaric isothermal ensemble ( nvt ) , depending on the stage of the simulation .
we consider a system of 1000 water molecules to model the amorphous water - ice surface as a slab of finite thickness .
the water molecules are treated as rigid and exist within a @xmath21 @xmath22 @xmath22 @xmath23 simulation box with periodic boundary conditions in the @xmath24 and @xmath25 directions .
the interaction potential for these water molecules is a tip4p / ice configuration @xcite .
this tip4p / ice model is a slight modification of the well - known transferable intermolecular potential model ( tip4p ) @xcite .
the basic form consists of a lennard - jones ( lj ; lennard - jones & devonshire 1936 ) term added to a coulombic term : @xmath26.\ ] ]    the lj term applies to a site that is on the oxygen atom but in the coulombic term , each h@xmath0o molecule is modeled by three charge sites : a positive charge on each hydrogen atom and a negative charge on the h - o - h molecular angle bisector , 0.1577 @xmath23 from the oxygen atom .
the variable @xmath27 is the distance between two oxygen atoms and @xmath28 is the distance between charge sites on two different h@xmath0o molecules , where @xmath29 and @xmath30 each run over the three charge sites of the two molecules .
the lj constants have values of @xmath31 @xmath23 and @xmath32 kcal mol@xmath33 .
the charge on each hydrogen site is @xmath34 , where @xmath35 is the electron charge ; the charge on the h - o - h molecular angle bisector is @xmath36 .
the h - o - h angle has a fixed value of 104.52@xmath37 , @xmath38 @xmath39 which is the fixed length of the o - h bond .
we use the potential above to generate the amorphous ice slab from a quenched md simulation . within our simulation box ,
we restrict the initial molecular position in the slab to a @xmath21 @xmath22 @xmath40 @xmath23 region .
the restriction on the @xmath25 coordinate ensures that we generate a slab rather than a cubic arrangement of water molecules .
we give each molecule a random initial position and velocity within this region and the system is then numerically solved .
a 10 @xmath23 cutoff was used for the coulombic interaction and ewald summation was used for long - range electrostatics .
quenching the slab occurs as follows : initially , the system evolves for 150 ps in nvt at a temperature of 300 k using velocity rescaling to set the temperature .
the system then evolves in nve for 50 ps , after which the system is cooled to 200 k and run again in nvt for 150 ps , followed by another 50 ps nve run .
the temperature is then lowered to 100 k. the process is repeated at 100 k before we cool the slab to 70 k. the reason it is desirable to perform a simulation in nvt is to ensure that the temperature of the amorphous ice slab , which is only a small component of an actual interstellar ice mantle , stays approximately constant . as a result of the nvt simulation , the substrate will settle into an energetically favorable configuration .
the system is then run in nve so that the substrate can be fixed in that energy configuration .    to generate slabs of other temperatures
we sequentially cool this 70 k slab in 10 k intervals , i.e. from 70 k to 60 k to ... to 10 k. the system runs in nve for 10 ps at each interval and velocity rescaling occurs every 5 ps .
a restart file is output at the end of each interval , containing the positions and velocities of the water molecules to be used as the initial conditions of the next cooling interval .
we tested this method by generating a 10 k slab from scratch and comparing calculated sticking probabilities on that slab to sticking probabilities on a 10 k slab made from the 70 k slab .
the sticking probabilities we calculated were in very good agreement , and we determined that this method is sufficient for slab generation at a variety of temperatures .
note that in this method of making the 10 k slab , we also make slabs of the other temperatures studied along the way .
these slabs are then used for hydrogen atom trajectory calculations .
a single hydrogen atom trajectory is initialized as follows .
first , we assign the hydrogen atom random initial @xmath41,@xmath42 coordinates and a @xmath25 coordinate fixed at 18  , about 8 @xmath39 above the ice surface . because the ice is amorphous the surface of the slab is ill defined , and therefore we use an average surface for referencing h - atom / surface distances .
random initial velocity components that correspond to a particular kinetic energy are assigned as well .
the initial velocity in the @xmath25-direction must be negative so that the hydrogen atom is launched towards the slab .
after the initialization , the trajectory is then integrated numerically using the leapfrog algorithm .
the interaction between the hydrogen atom and the water molecules of the substrate is described by the potential derived by zhang et al .
( 1991 ) , which has a form given by @xmath43y_{l , m}(\theta,\phi).\ ] ] this potential was developed using open - shell hartee fock and fourth - order mller - plesset perturbation theory .
the _ ab initio _ data were then fit to a lj form given by equation ( 5 ) .
the variable @xmath27 represents the distance between the impinging hydrogen atom and the oxygen atom .
the @xmath44 and @xmath45 coordinates are the spherical angles of the incident hydrogen atom with respect to the water molecule .
the origin is at the oxygen atom and the polar axis coincides with the molecular symmetry axis . for a water molecule ,
the spherical coordinates of the two hydrogen atoms are given by @xmath46  , @xmath47 , and @xmath48 .
the @xmath49 and @xmath50 parameters for this potential are given in table 1 .
the simulation time for a single run is 10 ps and the time - step is 1 fs . the simulation time was selected in an attempt to avoid ambiguity in determining the sticking or scattering status of the hydrogen atom .
periodic boundary conditions were used for the @xmath41 and @xmath42 coordinates ( parallel to the slab surface ) , but the simulation box was non - periodic in the @xmath25-coordinate ( normal to the surface ) .
thus the hydrogen atom leaves the simulation box  and the simulation terminates  if the @xmath25-coordinate of the atom ever exceeds 20  . during the simulation , the system runs in nve , but the velocity of the substrate molecules
are rescaled every 5 ps to ensure that the slab temperature remains ( approximately ) constant .    for nine incident kinetic energies of hydrogen atoms ranging from @xmath515 k - 400 k
, we performed 300 runs apiece for each of the 7 substrate temperatures .
thus , simulation data for a single substrate temperature comprises of 2700 runs , for a total of 18,900 runs over the 7 substrate temperatures .
we chose to perform 300 runs for each h - atom to obtain statistically significant values for the sticking probability .
for each of the 18,900 runs , the sticking or scattering status was determined .
we defined a sticking event to be a trajectory whose final @xmath25 coordinate ( @xmath52 ) was below 13 @xmath23 and whose final total energy ( @xmath53 ) was below -100 k. the physical motivation for this definition is discussed in the appendix . to get an estimate of error for each batch of runs we randomly selected 50 runs out of the 300 and calculated the sticking probability .
we repeated this random sampling 25 times and computed the standard deviation of the probabilities to represent the error of our computations .
mc and cnm regions of the ism contain dust grains covered with ice which have temperatures that range from 10 k to 100 k. we have performed simulations on 7 different ice substrate temperatures ( 10 k - 70 k ) for 9 different incident kinetic energies ( 5 k - 400 k ) of the h - atom .
this wide selection of slab temperatures provides a thorough coverage of the temperature ranges of interest .
all of our sticking probability and sticking coefficient results are shown in table 2 .
for comparison with previous studies ( e.g. , veeraghattam et al . 2014 and al - halabi et al .
2007 ) , we plot the sticking probability ( @xmath54 ) vs. the incident kinetic energies on 10 k and 70 k ice in figures 1 and 3 , respectively . in addition , in figures 2 and 4 , are plots of the sticking coefficient ( @xmath55 ) vs. the gas temperature for 10 k and 70 k ice , respectively .      because , to our knowledge , there have been no studies which calculated sticking probabilities and coefficients on substrate temperatures other than 10 k and 70 k , we focus the majority of our discussion on those two slab temperatures . figure 1 displays the sticking probability plotted vs. the incident kinetic energy of the hydrogen atom .
virtually all the atoms  stick " to the substrate at very low incident kinetic energies ( i.e. , 5 k and 10 k ) .
then , as the incident kinetic energy of the h atom increases , the sticking probability decreases
. the probability decreases slowly at first and then at around @xmath56 k , the probability begins to decrease more rapidly .
our results ( represented by dots ) show good agreement with those results of veeraghattam et al .
( 2014 ) ( represented by triangles ) , in spite of the fact that that study used a tip4p ( rather than a tip4p / ice ) potential for the water - water interaction and used a different sticking criterion ( see appendix ) . in addition , our results show strong agreement with those of al - halabi et al .
( 2002 ) , who also utilized a tip4p potential . with that being said , the results of al - halabi et al . ( 2002 ) were on a crystalline water - ice slab , rather than an amorphous one .
the results of buch & zhang ( 1991 ) and those of masuda et al .
( 1997 ) are not in good agreement with our study or those mentioned in the preceding paragraph .
buch & zhang ( 1991 ) obtain sticking probabilities that are much lower than those obtained here , whereas the sticking probabilities of masuda et al .
( 1997 ) are in general higher .
buch & zhang ( 1991 ) used a small cluster of ice comprising of 115 water molecules .
the cluster is highly irregular and as a result it is more difficult for the h atom to stick to the cluster @xcite .
masuda et al . ( 1997 ) performed simulations on a slab comprising of the same number of water molecules as our own ( 1000 ) but used a different potential to describe the water - water interaction , tip2s @xcite .
we cite the same discrepancy that al - halabi et al .
( 2002 ) reported to explain the difference between their results and those of masuda et al .
the study of masuda et al .
( 1997 ) applied an incorrect expression for @xmath57 that accounts for the greater chance of adsorption in their simulations @xcite .
figure 2 shows the sticking coefficient as a function of gas temperature .
the sticking coefficient is very close to unity at low gas temperatures ( 5 and 10 k ) and decreases as the temperature of the gas increases .
again , our results show strong agreement with those of veeraghattam et al .
( 2014 ) , despite the technical differences in our simulations .
the small differences between our calculated coefficients could possibly be attributed to the different number of runs that we performed , or the difference in time step .
regardless , all of the results of veeraghattam et al .
( 2014 ) lie within the uncertainties of our own results .
our results for the 70 k amorphous water - ice present a much more interesting case .
these results largely disagree with those of previous studies .
figure 3 , which displays the sticking probability as a function of incident kinetic energy , illustrates the discrepancies .
our calculated sticking probabilities lie outside of the error bars of all three studies of comparison .
overall the trends are similar : in all three studies , the sticking probability decreases as the incident kinetic energy of the hydrogen atom increases . in contrast with the 10 k slab , the sticking probability for atomic hydrogen does not go to unity at low incident kinetic energy . since the h@xmath0o molecules of the 70 k slab have greater average rotational and intermolecular vibrational energy than those of the 10 k slab , correspondingly larger energy transfers when the slab and h - atom collide yield a reduced likelihood of sticking .
the same reasoning is invoked to explain the differences between our results and those of masuda et al .
( 1997 ) as for the 10 k amorphous water - ice .
however , the discrepancies between our results and those of veeraghattam et al .
( 2014 ) are more subtle . the use of a new potential , tip4p / ice instead of tip4p , does not account for the discrepancy .
we performed simulations using the tip4p potential to test this and found that our results were still much higher than those of veeraghattam et al .
in fact , we found that the use of the new potential makes very little quantitative difference at all in these simulations .
the only remaining difference between our results and those of veeraghattam et al .
( 2014 ) is the criterion used for assessing sticking events .
we defined a sticking event to be a run in which the h - atom had a final @xmath25-coordinate that was less than 13 @xmath23 _ and _ had a final total energy less that -100 k. by contrast , veeraghattam et al .
( 2014 ) defined a sticking event as a run in which the h - atom had a final @xmath25-coordinate that was less than 11 @xmath23 and did not consider the final total energy of the h - atom .
these details make little difference on the 10 k ice , because at the end of the simulation atoms are very rarely still hopping around with high amplitudes . however , on 70 k ice it takes longer for atoms to give up their kinetic energy to the slab , because the water molecules in the slab have higher rotational and translational kinetic energies . as a result
, an atom may still be bouncing on the surface of the substrate at the end of the run .
we have selected a higher cutoff @xmath25 coordinate to include these bouncing atoms .
we have also opted for a longer simulation time : 10 ps , compared to the 5 ps of veeraghattam et al .
( 2014 ) , in order to allow for sufficient relaxation on the warmer substrate .
we also performed a few trial simulations of 20 ps to confirm that 10 ps was an adequate selection .
we found that the additional 10 ps had a negligible impact on the sticking or non - sticking status of the h - atom .
in addition , we have applied an energy cutoff of -100 k in the h - atom total energy as part of the sticking criterion .
this was done to align with the sticking criterion used by by al - halabi et al .
( 2002 ) . as a result of the contrasting definition with veeraghattam et al .
( 2014 ) , we observe large discrepancies between calculated probabilities and coefficients .
however , we do see that the sticking probability tends to be in closer agreement with the results of veeraghattam et al .
( 2014 ) as the incident kinetic energy of the h - atom increases .
this is because the  bouncing atom " scenario play less of a role at higher incident kinetic energies , where the h atom generally either sticks , or leaves the simulation box within the simulation time .
there are not many data points of al - halabi et al .
( 2002 ) to compare with and those that are present disagree with our predicted sticking probabilities .
that fact that al - halabi et al .
( 2002 ) performed their simulations on crystalline ice is the most likely explanation for the discrepancy , since all other features of the calculations ( potentials , sticking criteria , etc . )
are either similar or the same .    because the sticking coefficient is a function of the sticking probability , our calculated sticking coefficients are much larger than those of veeraghattam et al .
( 2014 ) , reflecting the discrepancy in our sticking probabilities .
this is shown in figure 4 , where the sticking coefficient is plotted against the gas temperature .
again there is better agreement at higher gas temperatures , where the particular criterion of a sticking event has less of an effect .
having discussed the results on the low - end and high - end of the slab temperatures studied , it is reasonable to assume that the sticking probabilities and coefficients for the in - between slab temperatures lie somewhere in - between the above results .
table 2 confirms that that is indeed the case .
we provide these results for astrophysical modeling purposes . as the slab temperature increases , in general , the sticking probability and coefficient decrease for each respective h - atom incident kinetic energy .
figure 5 shows the sticking coefficient for each gas temperature ( 5 k - 400 k ) as a function of substrate temperature .
this plot demonstrates some key features of these simulations . for all but the highest gas temperatures , the sticking coefficient decreases with increasing substrate temperature , as one would expect .
however , the overall drop in sticking coefficient from the value at substrates temperature of 10 k to 70 k decreases as gas temperature increases . by gas temperatures of 300 k and 400 k
, the sticking coefficient is very nearly constant in substrate temperature , at least over the range of temperatures studied .
this behavior can be understood by considering the average kinetic energy of the gas atoms . at high gas temperature
the average kinetic energy is so high compared to that of the substrate h@xmath0o molecules that energy transfers to and from the substrate are insensitive to substrate temperature @xcite .      as shown previously by @xcite , the sticking coefficients for 10 k ice obtained in our md simulations gives results significantly larger than those predicted by @xcite .
the latter is adopted in most astrochemical models involving ice - surface chemistry . since the formation of h@xmath0 on the ice surface is linearly proportional to the sticking coefficient , the adoption of the current values
would be expected to result in both an increased efficiency in the molecular hydrogen formed on the ice surface as well as faster removal of atomic hydrogen from the gas phase .
further , since the current sticking coefficients on 70 k ice can be as much as 30% larger than those calculated earlier in @xcite , we can expect the h@xmath0 efficiency to be relatively enhanced even as the surface temperature increases ( see fig .
5 ) until water begins to evaporate from the surface above @xmath58100  k. such an effect was found for the sticking of h@xmath0 and other molecular species on amorphous ice in the time - dependent , hot core warm - up modeling of @xcite .
we have reported the sticking probabilities and coefficients for nine gas temperatures of atomic hydrogen ( 10 - 400 k ) across seven different substrate temperatures ( 10 k - 70 k ) .
we present our results so that they may be used in future astrophysical modeling simulations .
the sticking of hydrogen atoms on amorphous water - ice is of significant importance in the ism , and its understanding is necessary before more complicated molecular species , which form through grain - mediated reactions , can be understood .
this study shows in detail how grain temperature affects sticking probabilities and coefficients , thereby showing how grain temperature affects the rate of h@xmath0 formation .
it is important however to point out some of the limitations of the study .
first , the nature of the simulation ( i.e. , launching a single h atom at a surface and determining its adsorption ) is different from the actual adsorption process in the ism , in which a dilute cloud of gas comprising of h atoms passes through these dusty regions .
in addition , the lack of a fully quantum treatment may be cause for some concern , because quantum effects tend to be more important for lighter atoms .
this study has shown that the criterion for assessing a sticking event can greatly effect sticking probabilities . in the appendix of the paper
we discuss the methods by which we have arrived at our sticking criterion .
we have explored many different criteria ( e.g. velocity , standard deviation of position , etc . ) for this process and have determined that our sticking criterion is the most appropriate choice for wide application .
further studies involving different species which have previous data ( e.g. h@xmath0 on h@xmath0o ice ) could confirm our criterion for sticking .
the authors would like to thank vijay veeraghattam for his useful comments and guidance throughout this entire project .
we also would like to thank the referee whose comments greatly improved the final version of the manuscript .
abascal , j. l. f. , sanz , e. , fernandex , r. g. , & vega , c. 2005 jchph , 122 , 23451    al - halabi , a. , kleyn , a. w. , van dishoeck , e.f .
, & kroes , g. j. 2002 , jpcb , 106 , 6515    al - halabi , a.,van dishoeck , e. f. , & kroes , g. j. 2004 , jchph , 120 , 3358    al - halabi , a. & van dishoeck , e. f. 2007 , mnras , 382 , 1648    allen , m. p. & tildesley , d.j .
1987 , _ computer simulation of liquids _
( oxford : claredon press )    andersson , s. , al - halabi , a. , kroes , g. , & van dischoeck , e. f. 2006 , jchph , 124 , 064715    eley , d. d. 1941 , rspsa , 178 , 452    hagan , w. , tielens , a. g. g. m. , & greenberg , j.m . , 1981 ,
cp , 54 , 367    harris , j. & kasimo , b. 1981 , sursc , 105 , l281 he , j. , acharyya , k. , & vidali , g. 2016 , , 823 , 56    hinshelwood , c. n. 1930 , ann . rep .
, 27 , 11    hollenbach , d. & salpeter , e. e. 1970 , jchph , 53 , 79    hollenbach , d. , kaufman , m. f. , bergin , e. a. , & melnick , g. j. 2009 , apj , 690 , 1497    hornekaer , l. , baurichter , l. , petrunin , v. v. , field , d. , & luntz , a. c. 2003 , sci , 302 , 1943 iqbal , w. , acharyya , k. , & herbst , e. 2012 , apj , 751 , 58    jorgenson , w. l. 1982 , jchph , 77 , 4156    langmuir , i. 1922 , trans .
faraday soc .
, 17 , 621 le bourlot , j. , le petit , f. , pinto , c. , roueff , e. , roy , f. 2012 , a&a , 541 , a76    lennard - jones , j. e. & devonshire , a. f. 1936 , natur , 137 , 1069    masuda , k. , takahashi , j. , & mukai , t. 1997 , a&a , 330 , 773    manico , g. , raguni , g. , pirronello , v. , roser , j. e. , & vidali , g. 2001 , apjl , 548 , l253    perets , h. b. , biham , o. , manico , g. , et al . 2005 , apj , 627 , 850    takahashi , j. , masuda , k. , & nagaoka , m. 1999 , mrnas , 306 , 22    tielens , a. g. g. m. 2005 , _ the physics and chemistry of the interstellar medium _ , cambrige university press    veeraghattam , v. , manrodt , k. , lewis , s. , & stancil , p. c. 2014 , apj , 790 , 4    watanabe , n. , kimura , y. , kouchi , a. , et al .
2010 , apjl , 714 , l233    zhang , q. , sabelli , n. , & buch , v. 1991 , jchph , 95 , 1080
the purpose of this section is to provide physical reasoning for our definition of a sticking event .
a sticking event occurs when a hydrogen atom has adsorbed to the water - ice substrate .
there are three main outcomes of these simulations : back scattering events , trapping - desorption events , and adsorption events @xcite .
a back scattering or elastic event is one in which the atom leaves the simulation box before the simulation ends .
this event is a clear non - sticking event .
an adsorption event is one in which the atom very clearly sticks to the surface within the timescale of the simulation .
this is determined by examining a ) the final position of the hydrogen atom and b ) the final total energy of the atom .
clear sticking events have a final position that is either embedded in the substrate or very close to the top of the substrate and low total energy , indicating that the atom has given up most of its kinetic energy and is bound to the surface .
trapping - desorption events , however , are difficult to define within classical md .
this is a case where an atom appears to be adsorbed to the substrate , but then gains kinetic energy from the substrate which causes the atom to desorb .
the uncertainty with regards to the desorption is that it may occur after the timescale of the md simulation .
trapping - desorption events have time scales on the scale of _ nanoseconds _ whereas our simulation has a timescale of 10 _ picoseconds _ @xcite . in the instance of trapping - desorption events
, an atom may be thermally hopping with large amplitude at the end of the run , yet the atom will still be confined to the simulation box . in classical md simulations these trapping desorption runs must be designated as a sticking or non - sticking event .
this designation is dependent on the particular definition of adsorption .
it is worth briefly considering the definitions of adsorption in previous studies of this nature .
buch & zhang ( 1991 ) , defines a sticking trajectory as one that has a final energy that is less than -100 k. they also state that no questionable runs occurred in their simulations .
either the atom stuck , or it was scattered .
k. masuda et al .
1997 studies the sticking probability of hydrogen atoms on the icy mantle of dust grains .
they define a sticking event as one where an  incident h atom resided on the surface of the amorphous water ice after 5000 fs of simulation time . " whereas a scattering / non - sticking event was one in which the atom returned to its initial height at the start of the simulation @xcite .
this definition neglects the possibility of trapping - desorption trajectories and gives no quantitative value for what was considered to be the cutoff for a sticking event .
al - halabi et al . ( 2004 ) who studied co adsorption on amorphous water ice , defines a sticking event as one that exhibits more than one turning point in the @xmath25 coordinate such that 1 ) the final @xmath25 coordinate of position @xmath52 @xmath59 @xmath39 and 2 ) the final energy of the adsorbed co molecule is below @xmath60 , where @xmath7 is the temperature of the slab , and @xmath61 is the boltzmann constant .
the values for @xmath52 are converted from the original number to allow comparison to our simulation box , in which the top of the substrate has an average height of 10  .
al - halabi et al . (
2002 ) employed a similar definition , but does without the @xmath52 cutoff and changed the final total energy cutoff to be -100 k. finally , veeraghattam et al .
( 2014 ) says  a sticking event occurs when the h atom stops bouncing on the surface and remains stuck for the rest of the simulation .
" however , we found that the h atom does not stop bouncing on the surface within the simulation time .
the h atom , even if embedded in the icy substrate continues to bounce with small amplitudes on or within the substrate . upon further investigation ,
veeraghattam et al . (
2014 ) applied a restrictive height cutoff value ( @xmath62  ) .
we initially took a sticking event to be any such run that leaves the simulation box before the timescale of the simulation . in doing so
, we noticed a discrepancy between our results and those of veeraghattam et al .
figures 6 and 7 , which illustrate the reason for this discrepancy , show the sticking probability as a function of cutoff value @xmath52 for 10 and 70 k substrates , respectively .
we let @xmath52 equal 100 values between 0.0 and 1.0 in reduced units , where 1.0 corresponds to 20 @xmath23 and 0.5 corresponds to 10  .
we then compute the sticking probability at each @xmath52 . from figure 6
, we can see that on a 10 k substrate the sticking probabilities plateau fairly quickly .
the approximate top of the substrate is at 0.5 ( 10  ) in reduced units .
probabilities plateau around 0.55 - 0.6 , which is about 11 - 12  .
in addition , in the region of 0.55 - 0.6 the probabilities are not varying significantly .
although uncertainties are not plotted here for aesthetic reasons , the probabilities in this region would be within uncertainties .
compare this to figure 7 , in which probabilities plateau around 0.65 ( 13  ) , and the probabilities are very different from those calculated at 0.55 ( 11  ) .
figure 7 further demonstrates the reason for our higher sticking probabilities when compared to those results of veeraghattam et al .
( 2014 ) on 70 k ice , while figure 6 demonstrates the reason that our results agree on 10 k ice . as a result of figure 7
, we decided to select @xmath63 as our cutoff value in all sticking calculations .
however , based on the studies of al - halabi et al .
( 2002 , 2004 , and 2007 ) , we decided apply a second criterion for sticking . after exploring other parameters such as
the final mean @xmath25 velocity of the h atom in the last 0.3 ps of a run and the standard deviation in @xmath25 position of the h atom in the final 0.3 ps of a run , we decided that final total energy ( @xmath53 ) was the appropriate choice .
figures 8 and 9 plot the final total energy value against the final @xmath25 coordinate value for each of the 300 runs of a 150 k atom on 10 and 70 k ice , respectively .
atoms in the dense region towards the upper right of each plot are clear non - sticking events .
these are atoms that have left the simulation box within the timescale of the simulation .
atoms whose location is near to the bottom left corner of the plots are definite sticking events .
these are atoms that are embedded in the slab with virtually zero kinetic energy .
comparing the middle regions of figure 8 to figure 9 , one is able to notice some significant differences .
for instance , in figure 8 , the atoms on the 10 k substrate that have not left the simulation box ( i.e. those atoms not in the upper right hand corner ) are relatively densely packed when compared to those atoms on the 70 k substrate .
this shows that the higher temperature of the substrate is affecting where the atoms will be and what energy the atoms will have at the end of the simulation .
note also that on the 10 k substrate the atoms seem to be grouped together in clumps , whereas there is no grouping behavior on the 70 k substrate .
the reasons for this clumping behavior are not completely clear at the present but they are likely due to regions of our slab configuration which are energetically favorable for the h - atoms to settle .
there appear to be  pockets " in our slab that the atoms have a preference to congregate towards .
by contrast , in the 70 k substrate , there are no energetically favorable sites for the atoms to settle into so all regions of z@xmath64 and final total energy are equally likely . to summarize , the clumping behavior is likely due to the differences in the particular amorphous ice slab configuration between the 10 k and 70 k slabs .
further analysis would be necessary to confirm this claim .
based on figures 8 and 9 we can draw a box around the atoms which have stuck to the substrate using @xmath52 and @xmath53 .
the atoms in the box have adsorbed to the substrate whereas those atoms which are outside the box have not adsorbed to the substrate .
figure 10 displays how the selection of @xmath53 with a fixed @xmath65 affects sticking probability on 70 k ice . note that we only show the 70 k substrate here because it is the 70 k substrate that we are using to choose our sticking event parameters .
we see that the sticking probability plateaus around -500 k , but in order to be conservative about our selection and to agree with the parameter of al - halabi et al .
( 2002 ) , we select @xmath66 k.        c c c c c c 10 & 5 & 0.993 & 0.014 & 0.990 & 0.017 + & 10 & 0.993 & 0.017 & 0.980 & 0.023 + & 50 & 0.957 & 0.040 & 0.860 & 0.045 + & 100 & 0.917 & 0.049 & 0.647 & 0.046 + & 150 & 0.770 & 0.055 & 0.490 & 0.040 + & 200 & 0.567 & 0.068 & 0.381 & 0.033 + & 250 & 0.457 & 0.060 & 0.304 & 0.027 + & 300 & 0.373 & 0.042 & 0.248 & 0.023 + & 400 & 0.220 & 0.024 & 0.173 & 0.016 + 20 & 5 & 1.000 & 0.000 & 0.976 & 0.017 + & 10 & 0.977 & 0.030 & 0.948 & 0.022 + & 50 & 0.917 & 0.027 & 0.800 & 0.046 + & 100 & 0.843 & 0.064 & 0.594 & 0.047 + & 150 & 0.640 & 0.059 & 0.455 & 0.039 + & 200 & 0.553 & 0.066 & 0.361 & 0.032 + & 250 & 0.410 & 0.056 & 0.295 & 0.027 + & 300 & 0.297 & 0.037 & 0.247 & 0.022 + & 400 & 0.223 & 0.021 & 0.180 & 0.017 + 30 & 5 & 0.963 & 0.031 & 0.973 & 0.027 + & 10 & 0.967 & 0.027 & 0.979 & 0.021 + & 50 & 0.923 & 0.041 & 0.812 & 0.048 + & 100 & 0.800 & 0.062 & 0.613 & 0.050 + & 150 & 0.770 & 0.051 & 0.467 & 0.043 + & 200 & 0.523 & 0.059 & 0.365 & 0.035 + & 250 & 0.440 & 0.054 & 0.292 & 0.029 + & 300 & 0.377 & 0.064 & 0.239 & 0.024 + & 400 & 0.220 & 0.033 & 0.167 & 0.017 + 40 & 5 & 0.937 & 0.039 & 0.940 & 0.036 + & 10 & 0.937 & 0.034 & 0.941 & 0.037 + & 50 & 0.900 & 0.046 & 0.766 & 0.051 + & 100 & 0.743 & 0.057 & 0.572 & 0.050 + & 150 & 0.647 & 0.063 & 0.437 & 0.044 + & 200 & 0.547 & 0.057 & 0.345 & 0.037 + & 250 & 0.420 & 0.052 & 0.278 & 0.031 + & 300 & 0.310 & 0.047 & 0.230 & 0.027 + & 400 & 0.210 & 0.032 & 0.164 & 0.020 + 50 & 5 & 0.917 & 0.033 & 0.916 & 0.043 + & 10 & 0.910 & 0.047 & 0.914 & 0.052 + & 50 & 0.870 & 0.064 & 0.733 & 0.058 + & 100 & 0.707 & 0.057 & 0.541 & 0.050 + & 150 & 0.610 & 0.054 & 0.408 & 0.040 + & 200 & 0.510 & 0.064 & 0.317 & 0.032 + & 250 & 0.387 & 0.044 & 0.252 & 0.026 + & 300 & 0.307 & 0.040 & 0.206 & 0.022 + & 400 & 0.180 & 0.023 & 0.144 & 0.015 + 60 & 5 & 0.853 & 0.055 & 0.861 & 0.061 + & 10 & 0.873 & 0.062 & 0.844 & 0.066 + & 50 & 0.767 & 0.069 & 0.662 & 0.060
+ & 100 & 0.647 & 0.058 & 0.501 & 0.048 + & 150 & 0.530 & 0.059 & 0.388 & 0.038 + & 200 & 0.457 & 0.047 & 0.306 & 0.031 + & 250 & 0.410 & 0.034 & 0.247 & 0.025 + & 300 & 0.340 & 0.033 & 0.203 & 0.021 + & 400 & 0.200 & 0.021 & 0.144 & 0.015 + 70 & 5 & 0.747 & 0.061 & 0.768 & 0.058 + & 10 & 0.780 & 0.063 & 0.771 & 0.052 + & 50 & 0.730 & 0.052 & 0.635 & 0.060 + & 100 & 0.633 & 0.078 & 0.481 & 0.053 + & 150 & 0.540 & 0.061 & 0.373 & 0.044 + & 200 & 0.433 & 0.062 & 0.297 & 0.036 + & 250 & 0.370 & 0.030 & 0.242 & 0.029 + & 300 & 0.290 & 0.046 & 0.201 & 0.024 + & 400 & 0.197 & 0.030 & 0.145 & 0.018