Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: for many years , the dominant method for discovering massive clusters of galaxies at @xmath11 , has been via x - ray emission from the hot gas in their dark matter potential wells . over the years , the _ rosat _ and _ xmm - newton _ telescopes , in particular , have yielded a handful of well - studied examples e.g. , rdcs j0910 + 5422 @xcite , rdcs j1252.9 - 2927 @xcite , rx j0848 + 4452 ( lynx e ; @xcite ) , xmmu j2235 - 25 @xcite , and xmmxcs j2215.9 - 1738 @xcite . in recent years , it has been realised that , by incorporating deep infrared ( ir ) observations , existing optical imaging techniques can be adapted to successfully detect clusters at redshifts competitive to the existing x - ray surveys . however , because ir observations covering only modest areas ( 120 arcmin@xmath1 - 7.25 deg@xmath1 ) have been available , the ( cluster or candidate ) systems detected to date have been less massive than those discovered from the x - ray surveys ( @xcite ; although see @xcite ) . massive clusters of galaxies are rare , and one requires as widefield a survey as possible to detect them . the largest area _ space telescope survey is the 50 square degree seven passband ( 3.6 , 4.5 , 5.8 , 8.0 , 24 , 70 , 160 @xmath12 ) spitzer wide - area infrared extragalactic survey ( swire ) legacy survey @xcite . we have obtained deep @xmath13band imaging @xcite , and combined this with the pre - existing infrared array camera ( irac ; @xcite ) observations from the swire survey , aiming to select clusters at @xmath14 using a two filter @xmath15 $ ] infrared adaptation of the well - proven optical cluster red sequence ( rs ) method @xcite . the sparcsgillianw / sparcs ] survey @xcite , is now complete and has an effective area ( defined as the usable overlap with swire after excluding chip gaps , regions near bright stars , etc ) , of 41.9 deg@xmath16 . the sparcs catalog contains several hundred cluster candidates at @xmath14 , by far the largest , homogeneously selected sample of its kind . in @xcite and @xcite we presented an overview of the sparcs survey , and reported spectroscopic confirmation of three clusters at @xmath17 and @xmath18 . in this paper , we report spectroscopic confirmation of three additional clusters at @xmath19 and @xmath20 .    the paper is organized as follows : in  [ observations ] , we describe the imaging and spectroscopic observations , and the data reduction ; in  [ analysis ] , we describe the spectroscopic catalog ; in  [ results ] , we present our results for the three individual clusters , and estimate the cluster velocity dispersions and dynamical masses ; and in  [ discussion ] , we present a discussion and the main conclusions based on our results . we assume a @xmath21-cdm cosmology with @xmath22 and @xmath23 , and h@xmath24 km s@xmath5 mpc@xmath5 . the sparcs / swire survey is comprised of six fields ( see table 1 of * ? ? ? * ) . the three clusters described in this paper were identified as high significance cluster candidates in the 7.9 deg@xmath1 elais - n1 field . a full description of the sparcs data reduction , cluster candidate detection algorithm and catalogs will appear in muzzin et al . 2010 , in preparation . we provide only a brief overview of the main details here . sparcs @xmath0 observations of the northern fields were observed by cfht / megacam for a total integration time of 6000s per pointing . photometry was performed on both the @xmath25 and irac mosaics using the sextractor photometry package @xcite . the total 3.6 @xmath12 magnitude of all sources in the field was computed using a large aperture , equal to the geometric mean radius of the sextractor isophotal aperture . further details of the photometric pipeline may be found in @xcite and @xcite . the depth of the @xmath0 data varies from pointing to pointing depending on the seeing and the sky background , however , the mean 5 @xmath26 depth for extended sources in the elais - n1 field was @xmath27 vega ( 24.2 ab ) . the @xmath28 color of all sources was computed using an aperture of diameter three irac pixels ( @xmath29 ) . no aperture corrections were applied to the @xmath25 photometry because the @xmath30 color aperture was much larger than the median seeing . however , the irac point spread function ( psf ) has broad wings compared to a typical ground - based seeing profile , and it was necessary to apply aperture corrections to the measured 3.6 @xmath12 magnitude before computing the color of each galaxy . galaxy clusters were identified in the sparcs survey , using an algorithm very similar to that described in detail in @xcite , as applied to the spitzer first look survey ( fls ; @xcite ) . for cluster detection in both the fls and sparcs surveys , we used an infrared adaptation of the cluster red - sequence ( crs ) technique @xcite . this algorithm maps the density of galaxies in a survey within narrow color slices , giving greater weight to brighter galaxies , and flagging the highest overdensities as candidate clusters . the one important difference here , compared to @xcite , is that the latter paper used an @xmath31 $ ] color to detect clusters at @xmath32 in the fls . the deeper irac data in swire ( @xmath33s frames ) , compared to the fls ( @xmath34s frames ) , combined with the @xmath15 $ ] color choice , allow sparcs to detect clusters to higher redshift than was possible with the fls dataset .    from analysis of the richnesses and colors , sparcs j161641 + 554513  and sparcs j161037 + 552417 were identified as rich cluster candidates at @xmath11 , and sparcs j161315 + 564930  as an unusually rich cluster candidate at @xmath35 ( table  [ tab_clusters ] ) . spectroscopic observations of sparcs j161315 + 564930 , sparcs j161641 + 554513  and sparcs j161037 + 552417  were obtained with the low resolution imaging spectrograph ( lris ; @xcite ) on the keck i telescope on the nights of april 3th and april 4th , 2008 . the seeing on both nights was about @xmath36 . lris has beam - splitters that separate the light between two sides : red @xcite and blue @xcite . on the red side , we used the 400 line mm@xmath5 grating which gives a coverage of @xmath37   centered at 8500  , with a dispersion of 1.86  / pixels . all the masks were designed with slits of 11 width , which gave a resolution of @xmath387 - 9   ( fwhm ) over the wavelength interval @xmath386600 - 9600  . for the redshift range of these particular clusters , all of the prominent spectral features e.g. , the 4000  - rest - frame break or [ o@xmath39 emission line , fall in the lris red side wavelength window . we observed one mask each in the case of clusters sparcs j161315 + 564930  and sparcs j161641 + 554513 , and two in the case of sparcs j161037 + 552417 . the total exposure time was 6000s ( @xmath40 ) for sparcs j161315 + 564930 , 8400s ( @xmath41 ) for sparcs j161641 + 554513 , and 8400s ( @xmath42 ) and 7000s ( @xmath43 ) for mask @xmath44 and mask @xmath45 of sparcs j161037 + 552417  ( table  1 ) . the masks contained about 30 slits each ( five of these were alignment stars and the remaining 25 were galaxies , selected according to a sliding scale of priority ) .    to reduce the number of slits placed on obvious foreground galaxies , we prioritized slits on galaxies with colors near each cluster s red - sequence . to avoid significant selection bias , we used a very broad cut around the red - sequence , intended to include both star forming and non - star forming systems . slits were placed on galaxies with priorities from 1 to 4 . priority 1 was galaxies with colors within 0.2 magnitude of the rs , and with @xmath46 < 17.0 $ ] ( vega ) . priority 2 was galaxies with colors between 0.2 and 0.6 magnitudes bluer than the rs , and with @xmath47 < 18.7 $ ] . priority 3 was galaxies with same colors as priority 1 , but with @xmath47 < 18.7 $ ] , and priority 4 was galaxies with same colors as priority 2 , but with @xmath47 < 18.7 $ ] . priorities 1 to 4 roughly correspond to bright red - sequence , blue cloud , faint red - sequence , and faint blue cloud galaxies respectively . approximately ten priority 1 galaxies could be accomodated per mask for sparcs j161315 + 564930 , and five per mask for sparcs j161641 + 554513  and sparcs j161037 + 552417 .    to reduce the lris data , we adapted custom software , based on that developed by @xcite to reduce vlt / fors2 data . overscan and bias corrections were applied to both the calibration ( bias , flats and lamp arcs ) and to the science ( clusters and standard stars ) frames using the iraf _ lrisbias _ and _ lccdproc _ tasks developed by g. d. wirth and c. fassnacht . after these corrections had been applied , the data were next corrected for geometric distortions across the dispersion axis , and then separated into individual slitlets and reduced using standard long - slit techniques . each individual slitlet was processed using an algorithm similar to that implemented in _ bogus _ ( developed by d. stern , a. j. bunker and s. a. stanford)stern / homepage / bogus.html ] . the flatfield slitlets were normalized and applied to the corresponding science slitlets . the individual 2d spectra were then background - subtracted and fringe - corrected before being stacked to produce a final co - added 2d spectrum . slits longer than 90 pixels ( @xmath48 ) were further corrected for  long - slit " distortions . the 1d spectrum of the source ( or sources ) in each slitlet was then extracted using standard iraf tasks . a set of near lamp exposures were used to wavelength - calibrate the observations . lris spectra can suffer significant flexure which may introduce wavelength offsets as large as @xmath49   in the wavelength solution , in the most extreme cases . the offsets in the wavelength calibration due to flexure were corrected by comparing to the wavelengths of known skylines . we estimate the residual uncertainties in wavelength calibration to be at most 1  , which correspond to a redshift uncertainty of @xmath50 at @xmath51  . finally , a relative flux calibration was performed , using a sensitivity function spanning the range 5600 - 9400  which was obtained from long - slit observations of the standard stars feige 34 and feige 67 @xcite . a redshift solution was sought for each object by cross - correlating @xcite each of the lris spectra with observed galaxy templates @xcite using the task _ xcsao _ in the iraf _ rvsao _ package @xcite . if a redshift could be obtained for a spectrum , it was assigned a quality flag , `` q '' , with one of two possible values : 0 or 1 . high - confidence redshifts , unambiguously determined by two or more clear features in the continuum or in absorption , or by obvious emission lines such as [ o@xmath39 ( @xmath52 ) were assigned @xmath53 . lower - confidence redshifts based on ambiguous identification of a few spectral features , often in low s / n spectra , were assigned @xmath54 . in addition to a quality flag , an emission line flag , `` e '' , was also assigned to each galaxy , based on the presence or absence of emission line features . spectra displaying one or more emission feature such as h@xmath55 , [ o@xmath56(@xmath574959,@xmath575007 ) or h@xmath58 in the case of low - redshift ( @xmath59 ) galaxies , or [ o@xmath39(@xmath573727 ) in the case of higher redshift galaxies , were assigned @xmath60 . spectra showing no sign of excess emission above the stellar continuum , were assigned @xmath61 . the total number of galaxies ( cluster members and foreground or background galaxies , excluding stars ) extracted from the fields of sparcs j161315 + 564930 , sparcs j161037 + 552417  and sparcs j161641 + 554513  were 22 , 19 and 34 ( see tables [ tab_en1_240 ] , [ tab_en1_349 ] and [ tab_en1_381 ] ) . occasionally , one or more serendipitous sources fell within a slit . the breakdown by quality flag was 20 , 13 and 14 @xmath53 galaxies , and 2 , 6 and 20 @xmath54 galaxies , respectively . example spectra for a subsample of cluster members ( see section  [ cl240 ] for the definition of a `` cluster member '' ) are shown in figure  [ cl_specs ] . the left column shows examples of sparcs j161315 + 564930  members , the center column shows examples of sparcs j161641 + 554513  members , and the right column shows exampels of sparcs j161037 + 552417  members . the i d of each object is indicated ( tables [ tab_en1_240 ] , [ tab_en1_349 ] and [ tab_en1_381 ] ) . fluxes are in relative units , and smoothed by 7 pixels ( 1 pix @xmath62  ) . prominent spectral features are indicated by vertical lines . [ results ]      in determining cluster membership , and calculating a velocity dispersion and mass , only high confidence , @xmath53 , galaxies were considered . definitive cluster membership was determined using the code of @xcite , which is based on the shifting - gap technique of @xcite . this procedure uses both galaxy angular position and radial velocity information to exclude near - field interlopers . the squares in the left panel of figure  [ z_histo_240 ] show the 14 @xmath53 galaxies identified as members of cluster sparcs j161315 + 564930  ( see table  [ tab_en1_240 ] ) by the shifting - gap technique . for sparcs j161315 + 564930 , these @xmath53 cluster members fall in the range @xmath63 , indicated by the vertical dashed lines shown in the right panel of figure  [ z_histo_240 ] , which shows the redshift histogram for the sparcs j161315 + 564930  field . the properties of the cluster members and non - members are summarized in table  [ tab_en1_240 ] . the total 3.6 @xmath12 magnitude ( column 4 in table  [ tab_en1_240 ] ) , and @xmath28 color ( column 5 in table  [ tab_en1_240 ] ) were calculated as described in section  [ photometry ] . in determining membership , we did not require that any @xmath54 galaxies satisfy the shifting - gap criteria ; if their redshifts fell in the range @xmath63 we included them as  cluster members " in table  [ tab_en1_240 ] , but re - emphasize that we do not utilize them in estimating the cluster mean redshift or velocity dispersion . the total number of cluster members is 16 ( 14 with @xmath53 and 2 with @xmath54 ) , and six foreground / background galaxies . all confirmed members of cluster sparcs j161315 + 564930  were passive ( @xmath61 ) galaxies .    in some cases , a redshift was obtained which did not correspond to any galaxy in our photometric catalog . in the case of faint galaxies , this was because a spectroscopic redshift was obtained for a strong emission line galaxy whose continuum fell below the detection threshold of the @xmath25 catalog . in the case of bright galaxies , this was because of blending issues with the irac @xmath64 @xmath65 psf . these galaxies were assigned both a @xmath46$]-band magnitude and a @xmath66 $ ] color of 99 in table  [ tab_en1_240 ] . figure  [ proj_dist_240 ] shows @xmath67 $ ] color composites of cluster sparcs j161315 + 564930 . the @xmath68 data were obtained from wfc on the isaac newton telescope , and are available with the swire public data release @xcite . the white squares ( green circles ) overlaid on the right panel show the 16 cluster members ( and foreground / background galaxies in the @xmath69 fov ) with spectroscopically - confirmed redshifts from keck / lris ( see table  [ tab_en1_240 ] ) . once cluster membership was established , both the redshift and rest - frame velocity dispersion , @xmath70 , of sparcs j003550 - 431224 were then calculated iteratively using the `` robust estimator '' , @xmath71 @xcite . the robust estimator has been shown to be less sensitive than the standard deviation to outliers which may persist even after rejecting interlopers using the shifting - gap technique . the actual estimator used depends on the number of cluster members and is either the biweight estimator for datasets with at least 15 members , or , as here in the case of sparcs j161315 + 564930  with 14 @xmath53 members , the gapper estimator . the gapper estimator is discussed more fully in @xcite , @xcite and @xcite . the line - of - sight rest - frame velocity dispersion , @xmath72 , was calculated directly from the vector of spectroscopic redshift measurements , , as    @xmath73    where @xmath74 is the speed of light , and @xmath75 is the estimated dispersion of the measured redshifts with respect to the center of the distribution , @xmath76 .    a mean redshift of @xmath77 and a velocity dispersion of @xmath78 km s@xmath5 were calculated for sparcs j161315 + 564930  ( table  [ tab_clusters ] ) . the uncertainty on the latter was determined using jackknife resampling of the data . for comparison , a gaussian with an rms of 1230 km s@xmath5 has been overlaid on the redshift histogram in the right panel of figure  [ z_histo_240 ] . the line - of - sight rest - frame velocity dispersion can be used to calculate a dynamical estimate of @xmath79 , the radius at which the mean interior density is 200 times the critical density , @xmath80 , and @xmath81 , the mass contained within @xmath79 . in the spherical collapse model , at redshift @xmath82 , @xmath83 can be calculated from    @xmath84    where @xmath85 , is the hubble parameter at redshift @xmath82 , and    @xmath86    based on its velocity dispersion of @xmath78 km s@xmath5 , we estimate an @xmath87 mpc and a dynamical mass of @xmath88 m@xmath89 for sparcs j161315 + 564930  ( table  [ tab_clusters ] ) . although this mass is preliminary ( see section  [ richness ] ) , it seems likely that sparcs j161315 + 564930  is an unusually massive cluster , perhaps the most massive cluster in the entire sparcs survey , and comparable in mass to that of cluster ms 1054 - 03 at @xmath90 @xcite . the @xmath66 $ ] vs. @xmath0 color - magnitude diagram for all galaxies ( gray circles ) within a radius of @xmath83 ( @xmath91 mpc ) of the center of sparcs j161315 + 564930  is shown in figure  [ colmag_240 ] . this radius is approximately equal to the virial radius . spectroscopically - confirmed galaxies are shown by colored symbols . in this figure ( and in figures  [ colmag_349 ] and [ colmag_381 ] ) , cluster members ( @xmath53 or @xmath54 ) are shown by red circles and foreground / background galaxies by blue squares . note that there are several cluster members or foreground / background galaxies shown in table  [ tab_en1_240 ] , which fall within a projected radius of @xmath83 of the cluster center , but for which a color could not be determined . these galaxies do not appear in figure  [ colmag_240 ] ( or in figures  [ colmag_349 ] or [ colmag_381 ] , in the cases of tables  [ tab_en1_349 ] and tables  [ tab_en1_381 ] ) . seven @xmath53 galaxies were determined to be cluster members of sparcs j161641 + 554513  by the shifting - gap technique . these galaxies are shown by squares in the left panel of figure  [ z_histo_349 ] . these galaxies lie in the redshift range @xmath92 , indicated by the vertical dashed lines in the right panel of figure  [ z_histo_349 ] . table  [ tab_en1_349 ] summarizes the  cluster members " , the ten @xmath53 and @xmath54 galaxies with redshifts in this range , and the nine foreground / background galaxies . the right panel of figure  [ z_histo_349 ] shows the redshift histogram for the field of sparcs j161641 + 554513 . a gaussian with an rms of 950 km s@xmath5 has been overlaid . figure  [ proj_dist_349 ] shows @xmath67 $ ] color composites of cluster sparcs j161641 + 554513 . the white squares ( green circles ) overlaid on the right panel show the cluster members ( foreground / background galaxies ) with spectroscopically - confirmed redshifts ( see table  [ tab_en1_349 ] ) which fall within the @xmath93 fov of the image .    of the seven @xmath53 cluster members , five were classified as passive ( @xmath61 ) , one was classified as emission line ( @xmath60 ) , and one was classifed as an agn ( @xmath94 ) . the upper left white square in figure  [ proj_dist_349 ] corresponds to the agn . there were also three passive @xmath54 cluster members ( table  [ tab_en1_349 ] ) . the spectrum of the confirmed agn ( object i d @xmath95 854892 in table  [ tab_en1_349 ] ) is shown in the lowest panel of the center column in figure  [ cl_specs ] . this source shows nev(@xmath573346,@xmath573426 ) and neiii(@xmath573869,@xmath573968 ) in emission , which are common in agn , as well as prominent [ oii](@xmath573727 ) emission and high - order balmer lines ( h@xmath96 , h@xmath97 and h6 ) also in emission . the cluster mean redshift and velocity dispersion , were estimated iteratively from the seven @xmath53 members , using the gapper estimator . the mean redshift of sparcs j161641 + 554513  was calculated to be @xmath98 . the velocity dispersion was calculated to be @xmath99 km s@xmath5 , which corresponds to @xmath100 mpc , and a dynamical mass of @xmath101 m@xmath89(table  [ tab_clusters ] ) . the @xmath66 $ ] vs. @xmath0 color - magnitude diagram for all galaxies ( gray circles ) within a radius of @xmath83 ( @xmath102 mpc ) of the center of sparcs j161641 + 554513  is shown in figure  [ colmag_349 ] . the red circles and blue squares indicate those cluster members and foreground / background galaxies which lie within a projected distance of @xmath83 of the cluster center , and for which a color could be determined . seven @xmath53 galaxies were also determined to be cluster members of sparcs j161037 + 552417  by the shifting - gap technique . these galaxies are indicated by squares in the left panel of figure  [ z_histo_381 ] . two galaxies , indicated by crosses in figure  [ z_histo_381 ] ( i d @xmath95 s 727869 and 734082 in table  [ tab_en1_381 ] ) , were identified as near - field interlopers . the seven @xmath53 members lie in the redshift range @xmath103 , shown by the vertical dashed lines in the right panel of figure  [ z_histo_381 ] . table [ tab_en1_381 ] summarizes the  cluster members " , the ten @xmath53 or @xmath54 galaxies with redshifts in this range , and the 24 foreground / background galaxies . the right panel of figure  [ z_histo_381 ] shows the redshift histogram for the field of sparcs j161037 + 552417 . a gaussian with an rms of 410 km s@xmath5 has been overlaid . figure  [ proj_dist_381 ] shows @xmath67 $ ] color composites of cluster sparcs j161037 + 552417 . the white squares ( green circles ) overlaid on the right panel show the cluster members ( foreground / background galaxies ) with spectroscopically - confirmed redshifts ( see table  [ tab_en1_381 ] ) which fall within the @xmath104 fov of the image . of the seven @xmath53 cluster members , all were classified as emission line ( @xmath60 ) . of the three @xmath54 cluster members , one was classified as passive ( @xmath61 ) , and two as emission line ( @xmath60 ) . the cluster mean redshift and velocity dispersion , were calculated iteratively from the seven @xmath53 members , using the gapper estimator . the mean redshift of sparcs j161037 + 552417  was estimated to be @xmath105 and the velocity dispersion , @xmath106 km s@xmath5 , which corresponds to @xmath107 mpc and a dynamical mass of @xmath108 m@xmath89 ( table  [ tab_clusters ] ) .    sparcs j161037 + 552417  is significantly less massive than sparcs j161315 + 564930  or sparcs j161641 + 554513 . the much smaller @xmath83 estimated for sparcs j161037 + 552417 , combined with the geometric limitations of lris with respect to the redshift number density yield measurable from a single mask , resulted in a yield of only four spectroscopically confirmed members within @xmath83 from two masks . moreover , because of blending issues with irac s @xmath64 @xmath65 psf , a reliable @xmath66 $ ] color could not be determined for one of these galaxies ( i d @xmath95 4000013 in table  [ tab_en1_381 ] ) . figure  [ colmag_381 ] shows the @xmath66 $ ] vs. @xmath0 color - magnitude diagram for all galaxies ( gray circles ) within a radius of @xmath109 ( @xmath110 kpc ) of the center of cluster sparcs j161037 + 552417 . the three red circles denote the cluster members with i d @xmath95 s 723814 , 722784 , and 722712 in table  [ tab_en1_381 ] . color - magnitude diagrams for sparcs j161315 + 564930 , sparcs j161037 + 552417  and sparcs j161037 + 552417  were presented in figures  [ colmag_240 ] , [ colmag_349 ] , and [ colmag_381 ] . column 2 of table  [ tab_properties ] shows the @xmath15 $ ] color of the red sequence for sparcs j161315 + 564930 , sparcs j161037 + 552417  and sparcs j161037 + 552417 . the color is calculated from the mean color of the spectroscopically confirmed red - sequence cluster members . also shown in column 2 of table  [ tab_properties ] is the color of the rs for three additional clusters , previously reported in @xcite and @xcite . these clusters are sparcs j163435 + 402151  at @xmath111 , sparcs j163852 + 403843  at @xmath112 , and sparcs j003550 - 431224  at @xmath113 ( column 6 ) . for all six clusters a systematic uncertainty in the rs color of 0.15 magnitude has been assumed ( this uncertainty reflects the fact that , at present , we are using the zeropoints provided by elixir for the @xmath0  observations . we expect , in the future , to be able to reduce these photometric uncertainties , using our own internal calibration ) .    columns 3 , 4 and 5 of table  [ tab_properties ] show the redshift that would be estimated for each cluster based on the measured rs color ( column 2 ) , assuming a solar metallicity single burst bc03 model and a formation redshift of either @xmath114 , 4 or 10 . as can be seen from table  [ tab_properties ] ( and the left panel of figure  [ zsandmasses ] ) , the photometric redshift inferred from the measured @xmath66 $ ] color has a slight dependence on one s choice of formation redshift , although differences in color between the models at @xmath115 are fairly small ( @xmath116 between @xmath117 and @xmath118 ) . utilizing the colors from the @xmath119 model , the three new clusters presented here were assigned preliminarily redshift estimates of @xmath120 , @xmath121 , and @xmath122 ( table  [ tab_properties ] ) . these photometrically estimated redshifts are very similar to , albeit slightly lower than , the spectroscopically determined values of @xmath123 , @xmath124 , and @xmath125 ( column 6 in table  [ tab_properties ] ) . the @xmath126 $ ] color vs. spectroscopic redshift for all six sparcs clusters in table  [ tab_properties ] is plotted in the left panel of figure  [ zsandmasses ] . the solid , dotted and and dashed lines show the bc03 model colors as a function of redshift for formation redshifts of @xmath114 , 4 and 10 . it is clear from figure  [ zsandmasses ] that the agreement between the model colors and the observations is very good .    with a larger sample of clusters , there may turn out to be small but real discrepancies between the models and the measured rs colors . in the left panel of figure  [ zsandmasses ] , the @xmath119 model can be seen to be slightly redder than the observations in the case of five clusters ( or equivalantly , the inferred photometric redshift can be seen to be slightly lower than the spectroscopic redshift ) , but slightly bluer in the case of one cluster ( sparcs j003550 - 431224 ) . these small offsets in color between the models and the observations , can also be seen directly from figures  [ colmag_240 ] , [ colmag_349 ] , and [ colmag_381 ] . the solid lines in these three figures show the rs color predicted by the bc03 @xmath119 model _ at the spectroscopic redshift of the cluster _ , and can be seen to be slightly redder than the observed color . a more detailed comparison between the model predictions and the observations will be made in a future paper employing a larger sample of sparcs clusters .    despite the aforementioned caveats , and the issue of degeneracies between the photometric redshift of the clusters and the formation redshifts of their galaxies , our overall conclusion is that the inferred one - color photometric redshifts and the spectroscopic redshifts are in excellent ( @xmath127 ) general agreement . in addition to estimating the masses for sparcs j161315 + 564930 , sparcs j161641 + 554513 , and sparcs j161037 + 552417  from the galaxy line - of - sight velocity dispersion , we also estimated the masses from the richness of the clusters , using the @xmath128 richness parameter . @xcite introduced @xmath128 , an adaptation of the @xmath129 richness parameter , intended to utilize two - band photometry to increase the contrast of the cluster with the background , and therefore provide a measurement of the richness that is less sensitive to foreground / background large scale structures . @xmath129 is the amplitude of the three - dimensional , cluster center - galaxy spatial correlation function , @xmath130 @xcite .    instead of counting galaxies in a single passband , @xmath128 is obtained by counting galaxies in a color slice centered on the location of each cluster s red - sequence in the @xmath66 $ ] color - magnitude diagram . in computing @xmath128 , we used a slice bounded in color by @xmath66 = \pm0.3 $ ] of the best - fit rs color returned by the cluster finding algorithm ( @xcite ) , and bounded in magnitude by @xmath131 , where @xmath132 is the bc03 @xmath119 model prediction of the characteristic magnitude of a galaxy _ at the photometric redshift _ corresponding to that rs color ( table  [ tab_properties ] ) . the background galaxy counts were determined from the color distribution in the entire 7.9 deg@xmath1 elais - n1 field , minus the regions known to contain galaxy clusters . the @xmath128 richnesses of the three clusters were computed to be @xmath133 mpc@xmath134 ( sparcs j161315 + 564930 ) , @xmath135 mpc@xmath134 ( sparcs j161641 + 554513 ) , and @xmath136 mpc@xmath134 ( sparcs j161037 + 552417 ) . based on the empirical calibration of @xmath137 vs. @xmath138 determined by @xcite in the k - band for 15 cnoc1 clusters at @xmath139 , these richnesses correspond to @xmath140 for sparcs j161315 + 564930 , @xmath141 for sparcs j161641 + 554513 , and @xmath142 for sparcs j161037 + 552417 .    for comparison , columns 7 and 8 of table  [ tab_properties ] show the dynamical mass , @xmath143 , and richness mass , @xmath144 , estimates for all six sparcs clusters . although the uncertainties associated with both of the mass estimators are large , they are consistent with each other at the 1@xmath145 level for five out of the six clusters and at the 2@xmath145 level for the sixth . the agreement between the two mass estimators can be seen in the right panel of figure  [ zsandmasses ] . based on all six sparcs clusters spectroscopically confirmed to date , our conclusion is that , in addition to there being excellent agreement between the photometric and spectroscopic redshifts , there is also reasonable agreement between the cluster dynamical and richness mass estimates . the dynamical masses estimates should be considered preliminary at this stage . uncertainties will reduce as more data becomes available from a large spectroscopic follow - up program of sparcs clusters currently being carried out at the gemini telescopes . collectively , the six sparcs clusters confirmed to date ( the three clusters presented in @xcite and @xcite , plus the three clusters presented here ) , demonstrate that , given the availability of infrared observations , the rs technique is an efficient and effective method of detecting _ bona fide _ massive galaxy clusters at @xmath146 . moreover , our studies of these six clusters are showing that it is possible to infer fundamental parameters such as cluster redshift and mass _ from the survey data itself _ ( see also @xcite ) .    at @xmath147 , both the optical red - sequence cluster surveys , rcs-1 @xcite and rcs-2 @xcite , and the sloan digital sky survey ( sdss ; @xcite ) have shown that it is feasible to measure cosmological parameters from the evolution of the cluster mass function @xcite . in order to do this efficiently , the survey data themselves are used to detect clusters , and also to estimate the redshift and the mass of those clusters . the redshift is estimated from the red sequence color @xcite , and the mass is estimated from the optical richness @xcite . the fact that sparcs is also now demonstrating the practicality of estimating redshifts and masses at @xmath146 from the survey data alone is heartening for the current generation of surveys aiming to utilize optical - infrared high redshift cluster observations to constrain cosmological parameters e.g. , sparcs , the ukirt infrared deep - sky survey deep extragalactic survey ( ukidss dxs ; @xcite),and the irac shallow cluster survey ( iscs ; @xcite ) . these optical - ir surveys will provide complementary samples to those selected using the sunyaev - zeldovich effect , e.g. , the south pole telescope survey ( spt ; @xcite ) , the atacama cosmology telescope ( act ; @xcite ) , and the atacama pathfinder experiment ( apex ; @xcite ) . the complete sparcs catalog contains several hundred cluster candidates at @xmath14 . with new large , homogeneous , reliable @xmath11 catalogs becoming available from sparcs and other surveys in the very near future , the prospects look bright for high redshift cluster and cluster galaxy evolution studies in the coming years . +    we thank the referee for useful comments . the authors wish to recognize and acknowledge the very significant cultural role and reverence that the summit of mauna kea has always had within the indigenous hawaiian community . we are most fortunate to have the opportunity to conduct observations from this mountain . this work is based in part on archival data obtained with the spitzer space telescope , which is operated by the jet propulsion laboratory , california institute of technology under a contract with nasa . support for this work was provided by an award issued by jpl / caltech . gw also gratefully acknowledges support from nsf grant ast-0909198 , and from the college of natural and agricultural sciences at ucr . andreon , s. , maughan , b. , trinchieri , g. , & kurk , j.  2008 , , 507 , 147 becker , m.  r. , et al . 2007 , , 669 , 905 beers , t.  c. , flynn , k. , & gebhardt , k.  1990 , , 100 , 32 bertin , e. , & arnouts , s.  1996 , , 117 , 393 blindert , k.  2006 , ph.d . thesis , university of toronto , canada bremer , m.  n. , et al . 2006 , , 371 , 1427 brodwin , m. , et al . 2006 , , 651 , 791 bruzual , g. , & charlot , s.  2003 , , 344 , 1000 carlstrom , j.  e. , et al . 2009 , arxiv:0907.4445 demarco , r. , et al . 2005 , , 432 , 381 demarco , r. , et al . 2007 , , 663 , 164 dobbs , m. , et al . 2006 , new astronomy review , 50 , 960 eisenhardt , p.  r.  m. , et al . 2008 , , 684 , 905 fadda , d. , girardi , m. , giuricin , g. , mardirossian , f. , & mezzetti , m.  1996 , , 473 , 670 fazio , g.  g. , et al . 2004 , , 154 , 10 gilbank , d.  g. , yee , h.  k.  c. , ellingson , e. , gladders , m.  d. , barrientos , l.  f. , & blindert , k.  2007 , , 134 , 282 gioia , i.  m. , braito , v. , branchesi , m. , della ceca , r. , maccacaro , t. , & tran , k .- v . 2004 , , 419 , 517 girardi , m. , biviano , a. , giuricin , g. , mardirossian , f. , & mezzetti , m.  1993 , , 404 , 38 girardi , m. , demarco , r. , rosati , p. , & borgani , s.  2005 , , 442 , 29 gladders , m.  d. , yee , h.  k.  c. , majumdar , s. , barrientos , l.  f. , hoekstra , h. , hall , p.  b. , & infante , l.  2007 , , 655 , 128 gladders , m.  d. , & yee , h.  k.  c.  2000 , , 120 , 2148 gladders , m.  d. , & yee , h.  k.  c.  2005 , , 157 , 1 goto , t. , et al . 2008 , , 60 , 531 hoekstra , h. , franx , m. , & kuijken , k.  2000 , , 532 , 88 jee , m.  j. , white , r.  l. , ford , h.  c. , blakeslee , j.  p. , illingworth , g.  d. , coe , d.  a. , & tran , k .- v .  h.  2005 , , 634 , 813 kinney , a.  l. , calzetti , d. , bohlin , r.  c. , mcquade , k. , storchi - bergmann , t. , & schmitt , h.  r.  1996 , , 467 , 38 koester , b.  p. , et al . 2007 , , 660 , 239 kosowsky , a.  2003 , new astronomy review , 47 , 939 koyama , y. , kodama , t. , tanaka , m. , shimasaku , k. , & okamura , s.  2007 , , 382 , 1719 krick , j.  e. , surace , j.  a. , thompson , d. , ashby , m.  l.  n. , hora , j.  l. , gorjian , v. , & yan , l.  2008 , , 686 , 918 kurk , j. , et al . 2008 , astronomical society of the pacific conference series , 399 , 332 kurtz , m.  j. , mink , d.  j. , wyatt , w.  f. , fabricant , d.  g. , torres , g. , kriss , g.  a. , & tonry , j.  l.  1992 , astronomical data analysis software and systems i , 25 , 432 lacy , m. , et al . 2005 , , 161 , 41 lawrence , a. , et al . 2007 , , 379 , 1599 lonsdale , c.  j. , et al . 2003 , , 115 , 897 mccarthy , j.  k. , et al . 1998 , , 3355 , 81 mccarthy , p.  j. , et al . 2007 , , 664 , l17 mullis , c.  r. , rosati , p. , lamer , g. , bhringer , h. , schwope , a. , schuecker , p. , & fassbender , r.  2005 , , 623 , l85 muzzin , a. , yee , h.  k.  c. , hall , p.  b. , & lin , h.  2007 , , 663 , 150 muzzin , a. , wilson , g. , lacy , m. , yee , h.  k.  c. , & stanford , s.  a.  2008 , , 686 , 966 muzzin , a. , et al . 2009 , , 698 , 1943 oke , j.  b.  1990 , , 99 , 1621 oke , j.  b. , et al . 1995 , , 107 , 375 rosati , p. , stanford , s.  a. , eisenhardt , p.  r. , elston , r. , spinrad , h. , stern , d. , & dey , a.  1999 , , 118 , 76 rosati , p. , et al .  2004 , , 127 , 230 rozo , e. , et al . 2009 , , 703 , 601 rozo , e. , et al . 2009 , arxiv:0902.3702 ruhl , j. , et al . 2004 , , 5498 , 11 shapley , a.  e. , steidel , c.  c. , pettini , m. , & adelberger , k.  l.  2003 , , 588 , 65 shupe , d.  l. , et al . 2008 , , 135 , 1050 stanford , s.  a. , elston , r. , eisenhardt , p.  r. , spinrad , h. , stern , d. , & dey , a.  1997 , , 114 , 2232 stanford , s.  a. , holden , b. , rosati , p. , eisenhardt , p.  r. , stern , d. , squires , g. , & spinrad , h.  2002 , , 123 , 619 stanford , s.  a. , et al . 2005 , , 634 , l129 stanford , s.  a. , et al . 2006 , , 646 , l13 , j. , shupe , d.  l. , fang , f. , lonsdale , c.  j. , & gonzalez - solares , e. 2005 , ssc website release tran , k .- v . h. , franx , m. , illingworth , g.  d. , van dokkum , p. , kelson , d.  d. , blakeslee , j.  p. , & postman , m.  2007 , , 661 , 750 tonry , j. , & davis , m.  1979 , , 84 , 1511 van breukelen , c. , et al . 2007 , , 382 , 971 wilson , g. , et al . 2009 , , 698 , 1934 yee , h.  k.  c. , & ellingson , e.  2003 , , 585 , 215 yee , h.  k.  c. , & lpez - cruz , o.  1999 , , 117 , 1985 yee , h.  k.  c. , gladders , m.  d. , gilbank , d.  g. , majumdar , s. , hoekstra , h. , ellingson , e. , & the rcs-2 collaboration 2007 , arxiv : astro - ph/0701839 york , d.  g. , et al . 2000 , , 120 , 1579 zatloukal , m. , rser , h .- j . , wolf , c. , hippelein , h. , & falter , s.  2007 , , 474 , l5      sparcs j161315 + 564930 & 16:13:14.6 & + 56:49:29.9 & 1 & 16 & 14 & @xmath148 & @xmath149 & @xmath150 & @xmath151 + sparcs j161641 + 554513 & 16:16:41.3 & + 55:45:12.5 & 1 & 10 & 7 & @xmath152 & @xmath153 & @xmath154 & @xmath155 + sparcs j161037 + 552417 & 16:10:36.5 & + 55:24:16.6 & 2 & 10 & 7 & @xmath156 & @xmath157 & @xmath158 & @xmath159 + [ tab_clusters ]      1070485 & 243.32851 & 56.821579 & @xmath160 & @xmath161 & @xmath162 & 0 & 0 + 1071085 & 243.31200 & 56.828140 & @xmath163 & @xmath164 & @xmath165 & 0 & 0 + 1071707 & 243.34280 & 56.835220 & @xmath166 & @xmath167 & @xmath168 & 0 & 0 + 1072162 & 243.36000 & 56.840450 & @xmath169 & @xmath170 & @xmath171 & 0 & 0 + 1068630 & 243.33870 & 56.802010 & @xmath172 & @xmath173 & @xmath174 & 0 & 0 + 1072641 & 243.30150 & 56.845970 & @xmath175 & @xmath176 & @xmath177 & 0 & 0 + 2000012 & 243.31107 & 56.826047 & @xmath178 & @xmath178 & @xmath179 & 0 & 0 + 1070775 & 243.31090 & 56.824970 & @xmath180 & @xmath181 & @xmath182 & 0 & 0 + 1066451 & 243.39220 & 56.778770 & @xmath183 & @xmath184 & @xmath185 & 0 & 1 + 2000007 & 243.34483 & 56.772531 & @xmath178 & @xmath178 & @xmath186 & 0 & 0 + 1066059 & 243.35040 & 56.775040 & @xmath187 & @xmath188 & @xmath189 & 0 & 0 + 1069805 & 243.34230 & 56.814430 & @xmath190 & @xmath191 & @xmath192 & 0 & 0 + 1074041 & 243.28720 & 56.862221 & @xmath193 & @xmath194 & @xmath195 & 0 & 1 + 1072364 & 243.29550 & 56.842720 & @xmath196 & @xmath197 & @xmath198 & 0 & 0 + 1069395 & 243.32159 & 56.810188 & @xmath199 & @xmath200 & @xmath201 & 0 & 0 + 1073458 & 243.30260 & 56.855709 & @xmath202 & @xmath173 & @xmath203 & 0 & 0 + 1068289 & 243.33870 & 56.797932 & @xmath204 & @xmath205 & @xmath206 & 1 & 0 + 1071481 & 243.32381 & 56.832401 & @xmath207 & @xmath167 & @xmath208 & 0 & 0 + 1070036 & 243.27380 & 56.817139 & @xmath209 & @xmath210 & @xmath211 & 0 & 0 + 1067689 & 243.27271 & 56.791660 & @xmath212 & @xmath213 & @xmath214 & 1 & 0 + 1065709 & 243.34480 & 56.771900 & @xmath215 & @xmath216 & @xmath217 & 1 & 0 + 2000006 & 243.34517 & 56.771469 & @xmath178 & @xmath178 & @xmath218 & 1 & 0 + [ tab_en1_240 ]      848097 & 244.17630 & 55.748390 & @xmath219 & @xmath220 & @xmath221 & 0 & 1 + 853440 & 244.18610 & 55.764381 & @xmath166 & @xmath222 & @xmath223 & 0 & 0 + 849838 & 244.17340 & 55.753590 & @xmath224 & @xmath225 & @xmath226 & 0 & 0 + 849110 & 244.19260 & 55.751389 & @xmath227 & @xmath228 & @xmath229 & 0 & 1 + 851083 & 244.17180 & 55.757141 & @xmath230 & @xmath231 & @xmath232 & 0 & 0 + 837239 & 244.26920 & 55.715519 & @xmath233 & @xmath234 & @xmath235 & 1 & 0 + 846319 & 244.25529 & 55.743149 & @xmath236 & @xmath237 & @xmath238 & 0 & 1 + 3000009 & 244.26935 & 55.715914 & @xmath178 & @xmath178 & @xmath239 & 0 & 0 + 839226 & 244.27161 & 55.721340 & @xmath240 & @xmath241 & @xmath242 & 0 & 0 + 854892 & 244.21670 & 55.768921 & @xmath243 & @xmath244 & @xmath245 & 2 & 0 + 833647 & 244.16330 & 55.704411 & @xmath246 & @xmath247 & @xmath248 & 1 & 0 + 828383 & 244.22971 & 55.688580 & @xmath249 & @xmath250 & @xmath251 & 1 & 0 + 840274 & 244.26221 & 55.724548 & @xmath252 & @xmath253 & @xmath254 & 1 & 0 + 856147 & 244.18530 & 55.772751 & @xmath255 & @xmath256 & @xmath257 & 1 & 0 + 844301 & 244.21809 & 55.736950 & @xmath258 & @xmath259 & @xmath260 & 1 & 0 + 842789 & 244.15030 & 55.732430 & @xmath261 & @xmath262 & @xmath263 & 0 & 1 + 857305 & 244.18820 & 55.776299 & @xmath264 & @xmath265 & @xmath266 & 0 & 0 + 851794 & 244.17329 & 55.759270 & @xmath267 & @xmath234 & @xmath268 & 0 & 1 + 860365 & 244.18150 & 55.786140 & @xmath209 & @xmath269 & @xmath270 & 0 & 1 + [ tab_en1_349 ]      4000013 & 242.65226 & 55.405353 & @xmath178 & @xmath178 & @xmath271 & 1 & 0 + 723814 & 242.65221 & 55.404598 & @xmath272 & @xmath273 & @xmath274 & 1 & 0 + 722784 & 242.63251 & 55.401890 & @xmath275 & @xmath276 & @xmath277 & 0 & 1 + 718577 & 242.61940 & 55.391060 & @xmath278 & @xmath279 & @xmath280 & 1 & 0 + 722712 & 242.66769 & 55.401711 & @xmath281 & @xmath282 & @xmath283 & 1 & 1 + 4000026 & 242.60676 & 55.475206 & @xmath178 & @xmath178 & @xmath284 & 1 & 1 + 710318 & 242.60490 & 55.370651 & @xmath285 & @xmath286 & @xmath287 & 1 & 0 + 4000010 & 242.67026 & 55.393053 & @xmath178 & @xmath178 & @xmath288 & 1 & 0 + 728525 & 242.69780 & 55.416550 & @xmath289 & @xmath290 & @xmath291 & 1 & 0 + 4000032 & 242.60455 & 55.371836 & @xmath178 & @xmath178 & @xmath292 & 1 & 0 + 744358 & 242.57840 & 55.457119 & @xmath293 & @xmath294 & @xmath295 & 1 & 0 + 4000012 & 242.65180 & 55.404111 & @xmath178 & @xmath178 & @xmath296 & 1 & 0 + 723736 & 242.64661 & 55.404339 & @xmath297 & @xmath298 & @xmath299 & 0 & 1 + 709346 & 242.61549 & 55.368011 & @xmath300 & @xmath301 & @xmath302 & 0 & 1 + 714998 & 242.59080 & 55.381821 & @xmath303 & @xmath304 & @xmath305 & 1 & 1 + 710721 & 242.60280 & 55.371490 & @xmath306 & @xmath307 & @xmath308 & 0 & 0 + 751852 & 242.60730 & 55.476471 & @xmath309 & @xmath310 & @xmath311 & 0 & 1 + 734249 & 242.66650 & 55.430779 & @xmath312 & @xmath313 & @xmath314 & 0 & 1 + 747661 & 242.62151 & 55.465469 & @xmath315 & @xmath316 & @xmath317 & 1 & 1 + 740510 & 242.69310 & 55.447041 & @xmath255 & @xmath318 & @xmath319 & 0 & 1 + 708779 & 242.61410 & 55.366550 & @xmath320 & @xmath321 & @xmath322 & 1 & 1 + 725337 & 242.69600 & 55.408329 & @xmath323 & @xmath324 & @xmath325 & 0 & 1 + 4000015 & 242.60952 & 55.453100 & @xmath178 & @xmath178 & @xmath326 & 0 & 1 + 730800 & 242.66510 & 55.422260 & @xmath327 & @xmath328 & @xmath329 & 0 & 0 + 720317 & 242.64140 & 55.395630 & @xmath330 & @xmath331 & @xmath332 & 0 & 1 + 745152 & 242.66029 & 55.459110 & @xmath333 & @xmath334 & @xmath335 & 0 & 1 + 753137 & 242.63040 & 55.479912 & @xmath336 & @xmath337 & @xmath338 & 0 & 1 + 738292 & 242.63560 & 55.441189 & @xmath339 & @xmath340 & @xmath341 & 0 & 0 + 727869 & 242.68291 & 55.414661 & @xmath333 & @xmath342 & @xmath343 & 1 & 0 + 734082 & 242.71609 & 55.430450 & @xmath344 & @xmath345 & @xmath343 & 1 & 0 + 712763 & 242.63080 & 55.376591 & @xmath346 & @xmath347 & @xmath348 & 1 & 1 + 742510 & 242.67841 & 55.452202 & @xmath349 & @xmath350 & @xmath351 & 1 & 1 + 737166 & 242.70931 & 55.438370 & @xmath352 & @xmath353 & @xmath354 & 1 & 1 + 735789 & 242.69299 & 55.434792 & @xmath355 & @xmath356 & @xmath357 & 1 & 1 + [ tab_en1_381 ]      sparcs j161315 + 564930 & @xmath358 & @xmath359 & @xmath360 & @xmath361 & @xmath148 & @xmath151 & @xmath362 + sparcs j161641 + 554513 & @xmath363 & @xmath364 & @xmath365 & @xmath366 & @xmath152 & @xmath155 & @xmath367 + sparcs j161037 + 552417 & @xmath368 & @xmath369 & @xmath370 & @xmath371 & @xmath156 & @xmath159 & @xmath372 + sparcs j163435 + 402151 & @xmath373 & @xmath374 & @xmath375 & @xmath376 & @xmath377 & @xmath378 & @xmath379 + sparcs j163852 + 403843 & @xmath380 & @xmath381 & @xmath382 & @xmath383 & @xmath384 & @xmath385 & @xmath386 + sparcs j003550 - 431224 & @xmath387 & @xmath388 & @xmath389 & @xmath390 & @xmath391 & @xmath392 & @xmath379 + [ tab_properties ]
Abstract:

Answer: the spitzer adaptation of the red - sequence cluster survey ( sparcs ) is a @xmath0-passband imaging survey of the 50 deg@xmath1 spitzer swire legacy fields , designed with the primary aim of creating the first large , homogeneously selected sample of massive clusters at @xmath2 . sparcs uses an infrared adaptation of the two - filter cluster red - sequence technique . in this paper we report keck / lris spectroscopic confirmation of two new exceptionally rich galaxy clusters , sparcs j161315 + 564930  at @xmath3 , with 14 high - confidence members and a rest - frame velocity dispersion of @xmath4 km s@xmath5 , and sparcs j161641 + 554513  at @xmath6 , with seven high - confidence members ( including one agn ) and a rest - frame velocity dispersion of @xmath7 km  s@xmath5 . we also report confirmation of a third new system , sparcs j161037 + 552417  at @xmath8 , with seven high - confidence members and a rest - frame velocity dispersion of @xmath9 km s@xmath5 . these three new spectroscopically confirmed clusters further demonstrate the efficiency and effectiveness of two - filter imaging for detecting _ bona fide _ galaxy clusters at high redshift . we conclude by demonstrating that prospects are good for the current generation of surveys aiming to estimate cluster redshifts and masses at @xmath10 _ directly _ from optical - infrared imaging .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: stellar mass loss is of critical importance to the later stages of stellar evolution . stars evolving up the red giant branch ( rgb ) lose mass via stellar winds . the remaining mass of a star leaving the rgb tip is the primary determinant of its position on the horizontal branch ( hb ; e.g.  @xcite ) . subsequent mass loss on the asymptotic giant branch ( agb ) ejects the star s entire hydrogen envelope , creating a post - agb star and ( perhaps ) a planetary nebula ( pn ) . the timescale and end - point of agb evolution may therefore not be determined by the nuclear burning rate , but by mass loss @xcite . winds from these stars are the most significant producers of interstellar dust and therefore provided the heavy elements present in population i stars , including the sun and solar system @xcite . debate exists around many finer points of the mass loss process , particularly its variation with fundamental parameters such as luminosity , temperature and metallicity . firstly , the amount of mass loss that occurs on the rgb remains largely unmeasured . rgb mass loss appears to determine the ` second parameter ' ( after metallicity ) that defines the morphology of hbs in globular clusters ( gcs ) @xcite . secondly , it is not proven where dust formation stars , and whether rgb stars can form significant amounts of dust in their ejecta . if many rgb stars do produce dust , it would add a large number of dust factories ejecting dust into a galaxy s interstellar medium ( ism ) ( cf . the missing dust source in many galaxies ; e.g.  @xcite ) . @xcite suggest that dust production occurs episodically over a large portion of the rgb , while @xcite have presented evidence that dust production is confined to the rgb tip , though gaseous mass loss is still present all along both the rgb and agb @xcite . dust can also help drive a wind through radiation pressure on grains ( e.g.  @xcite ) , though it is not certain whether winds which lack high - opacity amorphous carbon can be accelerated in this manner @xcite . dustless stars typically lose mass via magneto - acoustic processes in their chromospheres @xcite . finally , the composition of circumstellar dust and the chemical processes by which it forms are only partly defined . of particular interest is the rle that metallic iron grains play in dust formation , and how metallicity affects dust production and composition @xcite .    in this study , we aim to address these points . to do this , we explore dust production in the galactic globular cluster 47 tucanae ( ngc 104 ) . globular clusters are unique laboratories with which we can study late - stage evolution of low - mass stars : globular clusters contain stars that are among the oldest in the universe , and most clusters are comprised of a single ( or at least dominating ) stellar population with little internal variance in age or chemical composition . comparisons between clusters therefore probe differences in mass loss created by metallicity , while comparisons within clusters probe differences created by stellar temperature and luminosity . throughout this paper , we will refer to similar studies : @xcite , which covers the cluster @xmath0 centauri ( ngc 5139 ; hereafter paper i ) ; @xcite , which covers the cluster ngc 362 ( hereafter paper ii ) ; and a work which accompanies this paper ( submitted apjs ; hereafter paper iii ) . the cluster 47 tuc itself is among the galaxy s most massive and , co - incidentally , one of the closest to us . its distance has been estimated photometrically as 41904700 pc @xcite and kinematically as 4000 pc @xcite , with 4500 pc being the median distance determination in the literature . its mass is 69 @xmath1 m@xmath2 @xcite . its metallicity is determined to be roughly [ fe / h ] = 0.7 , though various studies have placed it anywhere from [ fe / h ] = 0.89 to 0.66 @xcite . the cluster suffers from little interstellar reddening ( @xmath3 mag  @xcite ) and is estimated to be over 11 gyr in age @xcite . the cluster hosts an ionised intra - cluster medium ( icm ; @xcite ) . this ionised medium is undetected in any other cluster ( m15 has a neutral medium with a dusty component ; @xcite ) . integrating the observed intracluster electron density over the cluster core yields @xmath40.1 m@xmath2 of h ii , representing @xmath55% of the intracluster hydrogen @xcite . the total stellar mass loss from 47 tuc s stars therefore gives us insight into the fate of gas and dust expelled by stars in this cluster . paper iii more - precisely examines which stars have mid - infrared excess , indicative of dust , finding again that we can not corroborate the presence of dust around the fainter giants . in this work , we use the data and results from paper iii to derive mass - loss rates for individual stars and ( by summation ) the cluster as a whole , and investigate the composition of the circumstellar dust produced by analysing archival mid - infrared ( mid - ir ) spectra . the data for this study come from the accompanying paper iii , which identifies those stars in 47 tuc which have infrared excess , indicative of circumstellar dust . this list of 22 objects is an exhaustive list of those stars in the cluster which we can be relatively certain are producing dust , given currently available data . these objects are listed in table [ mdottable ] , which also includes the parameters we determine in later sections . @l@c@c@c@  c@  c@  c@  c@  c@  c@  c@c@  c@  c@  c@  c@  cl@ star & @xmath6 & @xmath7 & @xmath8 & @xmath9 & @xmath10 & & @xmath11 & @xmath12 & @xmath13 & @xmath14 & & notes +   & ( mag ) & ( km / s ) & ( days ) & ( k ) & ( l@xmath15 & @xmath16 & [ 3.6 ] & [ 8 ] & ( % ) & 10@xmath17m@xmath2yr@xmath18 ) & ( kms@xmath18 & ( k ) & sil . & fe & al@xmath19o@xmath20 & feo &   + v1 & 4.03 & 20 & 221 & 3623 & 4824 & 0.936 & 1.173 & 2.213 & 3.6 & 21.0 & 4.0 & 900 & 12 & 88 & 0 & 0 & + v8 & 1.6 & 16 & 155 & 3578 & 3583 & 0.904 & 1.210 & 2.135 & 5.9 & 14.7 & 4.0 & 900 & 7 & 85 & 7 & 1 & + v2 & 2.78 & 23 & 203 & 3738 & 3031 & 1.103 & & & 4.3 & 12.0 & 3.8 & 900 & 5 & 95 & 0 & 0 & @xmath21 + v3 & 4.15 & 22 & 192 & 3153 & 2975 & 0.954 & 1.198 & & 3.2 & 9.4 & 3.2 & 1000 & 0 & 100 & 0 & 0 & @xmath22 + v4 & 1.5 & 18 & 165 & 3521 & 2603 & 0.932 & 1.167 & 1.902 & 5.6 & 11.7 & 3.6 & 900 & 2 & 95 & 2 & 1 & @xmath4rgb - tip + v26 & 0.3 & & & 3500 & 2541 & 1.073 & 0.807 & 0.995 & 0.9 & 5.9 & 3.8 & 1000 & & 100 & & & @xmath23 + lw10 & 1.2 & & 221 : & 3543 & 2324 & 1.073 & 0.876 & 1.017 & 0.6 & 4.2 & 3.2 & 1000 & 0 & 100 & 0 & 0 & + v21 & 1.0 & 7 & 76 : & 3575 & 2301 & 0.977 & 1.083 & 1.355 & 1.1 & 4.9 & 3.6 & 1000 & 0 & 93 & 5 & 2 & + lw9 & 1.1 & & 74 & 3374 & 2204 & 0.973 & 1.003 & 1.228 & 1.1 & 6.0 & 3.3 & 950 & & 100 & & & + v27 & 0.3 & & 69 & 3374 & 2140 & 1.023 & 0.929 & 1.218 & 1.4 & 3.3 & 3.6 & 1050 & & 100 & & & + lee1424 & & & 65 & 3565 & 2122 & 1.006 & 1.182 & 1.390 & @xmath244 & 4.4 & 4.0 & 1050 & & 100 & & & + a19 & 0.7 & & 60 & 3526 & 2096 & 1.034 & 0.992 & 1.240 & 0.8 & 3.6 & 4.2 & 1100 & & 100 & & & + lw12&0.6&&116 & 3713 & 2079 & 1.082 & 0.967 & 1.093 & @xmath251 & low + x03 & & & & 3816 & 1640 & & 0.958 & 1.156 & 3.2 & 5.1 & 1.8 & 750 & & 100 & & & @xmath26 + lw19 & 0.1 & & 40 & 3738 & 1638 & 1.078 & 0.991 & 1.080 & @xmath251 & low + v20 & 0.8 & & 232 : & 3602 & 1575 & 1.053 & 0.928 & 1.039 & @xmath251 & low&&&&&&&@xmath27 + v22 & 0.5 & & 62 & 3684 & 1528 & & 2.504 & 2.834 & @xmath251 & low & & & & & & & @xmath23 + v6 & 0.6 & 7 & 48 & 3763 & 1374 & 1.001 & 1.000 & 1.167 & @xmath251 & low & & & & & & & + v5 & 0.7 & 8 & 50 & 3741 & 1363 & 1.003 & 1.100 & & @xmath251 & low&&&&&&&@xmath23 + v18 & 0.3 & 5 & 83 : & 3692 & 1297 & 0.989 & 1.004 & & 2.4 & 5.3 & 2.4 & 800 & 33 & 67 & 0 & 0 & @xmath26 + v23 & 0.5 & & 52 & 3775 & 1259 & 1.023 & 0.985 & & @xmath251 & low&&&&&&&@xmath23 + v13 & & 12 & 40 & 3657 & 1029 & 0.990 & 1.106 & & 1.2 & 4.1 & 2.2 & 800 & 0 & 85 & 15 & 5 & @xmath26 +   +      accurate  or even approximate  mass - loss rates are notoriously difficult to measure . one must assume that the wind has certain properties . perhaps the most crucial of these is the relationship between wind velocity and distance from the star . we here assume a constant , spherically - symmetric , purely radiatively - driven wind , starting from a near - stationary , dust - forming layer . as discussed previously ( paper i ; @xcite ) , this may not be the most appropriate description , but suffices due to lack of contrary evidence : we return to this point in the discussion .    before calculating an estimated mass - loss rate , we can create a measure of the amount of circumstellar dust surrounding a star , irrespective of most modelling uncertainties . to do this , we define a quantity @xmath11 , the fraction of the star s luminosity that is reprocessed by circumstellar dust . for low optical depths , this should be directly proportional to the surface area of dust grains , therefore approximately proportional to the mass of circumstellar dust . the value of @xmath11 can be determined by taking the difference between the observed sed and model photospheric sed in the infrared ( yielding the infrared excess ) and dividing this by the total photospheric luminosity ( the value given in paper iii ) . we list this quantity in table [ mdottable ] . errors in the value of @xmath11 are introduced by difficulties in determining the integrated luminosity from photometry alone and from the uncertainty in the photospheric spectrum ( both temperature and luminosity ) . the former can be approximated as @xmath28 , or roughly 50% of @xmath11 , but is dominated by photometric uncertainty at short ( 2.210-@xmath29 m ) wavelengths . the latter gives a constant uncertainty of around @xmath30% . we reduce this latter error by fixing the photospheric luminosity to the @xmath31-band flux for all stars except v26 and x03 , which are fixed to the 3.6-@xmath29 m flux . the value @xmath11 could not reliably be measured for several stars with low mass - loss rates . the stars appear to fall into three main categories : those with large ( @xmath45% ) reprocessed flux ( v14 , v8 ) , those with low ( @xmath42% ) reprocessed flux concentrated toward long wavelengths ( v13 and v18 ) , and those stars with low ( @xmath322% ) reprocessed flux concentrated toward shorter wavelengths ( the remainder ) .      given the limitations of our understanding and of the model presented above , we stress that the mass - loss rates and velocities we derive are based only on one model for wind acceleration , but should be comparative both internally and with other works which use these same assumptions .    under these assumptions , we calculate mass - loss rates by modelling the sed using the radiative transfer package , dusty @xcite . as inputs , we use our interpolated marcs model of the stellar spectrum , and a standard mrn grain size distribution @xcite . within our dusty model , we vary the dust composition , temperature in the dust - forming region ( the inner edge of the dust envelope ) and optical depth of the wind . we change these parameters iteratively until a best fit is achieved . we relate the dusty mass - loss rate ( @xmath33 ) to the actual mass - loss rate ( @xmath34 ) as previously ( paper i ; paper ii ) , using the following formula @xcite : @xmath35 , \label{mdoteq2}\end{aligned}\ ] ] for luminosity @xmath10 , dust - to - gas ratio @xmath36 , grain density @xmath37 , and terminal wind velocity @xmath38 . we assume that @xmath39 g @xmath40 for silicates , 4 g @xmath40 for aluminium oxides , 5.7 g @xmath40 for feo , and 7 g @xmath40 for metallic iron . for these three species , we use the optical constants from @xcite ; @xcite ; and @xcite , respectively . the errors in the above assumptions are substantial , and are discussed in section [ mdoterrorsect ] . the dust - to - gas ratio , @xmath36 , can be estimated as : @xmath41 } ,    \label{psieqn}\ ] ] which is approximately 1/1000 for 47 tuc . this assumes that the fraction of condensed metals is identical for all metallicities . we discuss the accuracy of this assumption in section [ mdoterrorsect ] . table [ mdottable ] shows our calculated mass - loss rates and inferred wind parameters and dust compositions alongside literature pulsation properties . the striking fact about the dust composition is the fraction of iron grains present . the presence of metallic iron in these stars is debatable , as its featureless infrared emission can be mimicked by many other phenomena , including large silicate grains and emission from thick circumstellar molecular layers @xcite . we consider iron to be the most likely source of emission , though the fractions of metallic iron listed in table [ mdottable ] may not be particularly accurate , and we discuss this further in  [ dustcompsect ] . many of the mass - loss rates are too low to be sensibly measured : we can merely say that some circumstellar dust probably exists . those rates that can be measured have their seds shown in figure [ sedsfig ] . the photometry is taken from the sources cited in paper iii . mid - ir spectra are included from @xcite ( hereafter vlmo ) , and @xcite ( hereafter lph ) . these objects are summarised in the following luminosity - ordered list :    = 1em    v1 : this star is clearly identified as the cluster s dominant dust producer by vlmo , lph and @xcite ( hereafter itm ) . v1 has clear excess at all wavelengths @xmath42 m . this excess , which we attribute to metallic iron dust , dominates even over its strong silicate feature . v8 : no silicate feature was seen in this star by vlmo , but a clear silicate and aluminium oxide feature was seen by lph . it has previously been suggested ( lph ; paper i ) that v8 may have variable mass loss . however , the _ akari _ photometry at 3 , 4 , 7 , 11 , 18 and 23 @xmath29 m shows no variation in the ( dominant ) iron contribution . this means that if the mass loss is variable , strong variation only occurs in the silicate production . v8 has a complex ir spectrum showing peaks from silicates , aluminium oxide , and from crystalline minerals at 9.7 , 11.3 , 13.1 , 20 and 32 @xmath29 m . we discuss the carrier(s ) of these bands in  [ dustcompsect ] . note that the amplitude of the 10-@xmath29 m silicate peak is substantially over - estimated in our model compared to the 20-@xmath29 m silicate peak . this implies a lower temperature for the silicate grains than the metallic iron grains , which we also discuss in  [ dustcondsect ] . v2 : a silicate dust feature was observed by lph , though it was too weak to have been seen by vlmo . iron dominates the ir flux , however . the 10- and 20-@xmath29 m silicate peaks are not well modelled , as in v8 ( see also  [ dustcondsect ] ) . v3 : no silicate emission was seen by either lph or vlmo , but strong , featureless ( metallic iron ) emission is present , which is confirmed by _ akari _ photometry . v4 : the spectrum in lph shows complex dust mineralogy , akin to that in v8 . this star is roughly at the rgb - tip , as estimated from the luminosity function . v26 : blended with a star half its brightness , 2.2@xmath43 to the east , the near - ir , _ akari _ and 24-@xmath29 m photometry are affected . the amount of ir excess is difficult to determine here , but is likely to be quite small . no discernible features are present in the spectrum from vlmo . lw10 : also within the dense cluster core , it appears that the _ akari _ photometry ( at 9 & 11 @xmath29 m ) for this object is affected . no features are detectable in the spectrum from vlmo . v21 : observed by lph , the spectrum shows the same features as v8 . these features , and the silicate contribution , are considerably weaker in v21 . v11 : no obvious features were seen in the spectrum of vlmo , this star was only covered by _ spitzer _ at 3.6 and 5.8 @xmath29 m . 9-@xmath29 m photometry suggests it may have some excess there , but both the _ spitzer _ 8-@xmath29 m photometry and the vlmo spectrum ( which covers 9-@xmath29 m ) show no excess . the derived mass - loss rate can probably be treated as a rough upper limit . lw9 & v27 : an excess consistent with circumstellar dust is present in both the _ spitzer _ and _ akari _ photometry . the chemistry can not be well - constrained by this photometry , but it appears that the flux contribution from components other than iron is comparatively small . lee1424 : lacks any data redwards of 9 @xmath29 m , making constraining the mass - loss rate , dust temperature and dust chemistry very difficult . the values given should be treated as order - of - magnitude estimates only . a19 : as lw9 & v27 , however the _ akari _ 7-@xmath29 m data suggests little to no mid - ir excess in this star . v20 & v22 : not shown in figure [ sedsfig ] , they have no near - ir photometry attributed to them , meaning it is difficult to determine the amount of ir excess present . despite this , the dust temperature is relatively - well constrained by the presence of _ spitzer _ 24-@xmath29 m photometry . x03 : located in the cluster core , this star suffers from the same problems as v20 & v22 , but clearly displays a stronger excess , making an approximate determination of the mass - loss rate and wind characteristics possible . v18 : observed by both lph and vlmo , this star has a clear silicate feature . this feature was not claimed in vlmo due to poor signal - to - noise , but the detection appears to match that from lph well . it shows substantial mass loss compared to other stars of similar luminosity . the 10- and 20-@xmath29 m silicate peaks are not well modelled , as in v8 ( see also  [ dustcondsect ] ) . v13 : highlighted in both lph and itm , the spectrum of this star shows an excess similar to that of v8 . the contribution of iron dust is comparatively small , though would appear to be non - zero . this star is isolated from the cluster core , but has poor - quality photometry due to its comparative faintness .      the accuracy of the wind parameters we find is hard to determine . the major sources of error are listed below :    1 .   a 30% uncertainty inherent in dusty due to an unspecified gravitational correction . this uncertainty is exceeded in stars with @xmath44 km s@xmath18 , where the effect of stellar gravity will change the variation of dust particle density with radius from the star ( the stellar mass is unspecified in dusty ) . 2 .   uncertainty in the underlying photospheric spectrum . as the iron dust contribution is essentially a modified blackbody , one can fit a dust temperature to the 3.6- to 8-@xmath29 m data and extrapolate this fit backwards to @xmath16-band , where the dust contribution is essentially zero . the photospheric contribution can then be set from the @xmath45 photometry , but the accuracy of this fit depends on the stellar photometry . ( this would equally hold for other possible dust species , e.g.  amorphous carbon . ) uncertainty in the amount of ir excess . apart from identifying how much flux is generated by the photosphere , missing photometry and photometric errors cause uncertainty in the combined flux of the ( photosphere + dust ) contributions . 4 .   uncertainty in the dust composition . it is not clear for a number of these stars ( v26 , lw9 , v27 , lee1424 , a19 , x03 and any star marked `` low '' ) what the dust chemistry is . many of these stars have insufficient coverage in the 923 @xmath29 m region to determine whether they contain any dust components other than metallic iron . lw10 and v11 are not covered by lph : ground - based spectra by vlmo show no obvious spectral features , so we presume a pure metallic iron wind , but the spectra are of low signal - to - noise . the presence of silicates can significantly increase the implied mass - loss rate due to their low opacity . any carbonaceous dust present would also be hidden by the modified blackbody of metallic iron . carbonaceous dust typically has a higher opacity : if it is present , the implied mass - loss rate would be decreased . uncertainty in grain size distribution and grain shape ( assumed here to be spherical , with radius @xmath46 ) . the relationship between the grains absorption cross - section to radiation ( @xmath47 ) and their mass ( @xmath48 ) . a high grain porosity would also imply a decreased dust density and therefore a higher mass - loss rate . particularly , differences in grain properties between the iron and silicate grains could explain why we model such a high proportion of iron grains in the wind ( see  [ dustcompsect ] ) . uncertainty in dust temperature . the dust temperature can only be determined to within @xmath4100 k. as @xmath49 , a 15% error in temperature leads to an error of up to @xmath460% in the flux emitted by the circumstellar dust , depending on the implied error in the radius of the dust shell . this will be most significant in stars that are losing the least mass , as their dust temperatures are more difficult to determine . uncertainty in the wind driving mechanism . we presume here that the winds are sustained by radiation pressure alone . the very low outflow velocities present in at least three cases ( x03 , v18 and v13 ) would suggest that there is insufficient radiation pressure for it to be the sole driver of mass loss . additional energy input , for example from pulsation , would increase the wind velocity and ( by virtue of eq . ( [ mdoteq2 ] ) ) decrease the mass - loss rate . there is a related uncertainty in the dust - to - gas ratio , which is assumed to scale with metallicity . these uncertainties can not be quantified statistically . lllll parameter & @xmath50 & @xmath51 & parameter variation used & other parameters affected + grain size@xmath52 & @xmath5310% & @xmath5330% & low : @xmath54 , @xmath55 nm , @xmath56 nm & @xmath57 k , 90% iron , 10% silicates +   &   &   & assumed : @xmath54 , @xmath58 nm , @xmath59 nm & @xmath60 k , 88% iron , 12% silicates +   &   &   & high : @xmath61 , @xmath58 nm , @xmath59 nm & @xmath60 k , 80% iron , 20% silicates + grain density & 50% & 50% & low : 1.63 g @xmath40 & +   &   &   & assumed ( solid ) : 6.52 g @xmath40 & + dust formation & @xmath5310% & @xmath5320% & low : 750 k & @xmath62 , 92% iron , 8% silicates +  temperature &   &   & assumed : 850 k & @xmath63 , 88% iron , 12% silicates +   &   &   & high : 950 k & @xmath64 , 82% iron , 18% silicates + dust : gas ratio & @xmath65 & @xmath66 & low : 1 / 5000 & +   &   &   & assumed : 1 / 1000 & +   &   &   & high : 1 / 200 & + velocity distribution & @xmath4 + 500% & n / a & assumed : 0 @xmath67 v@xmath68 km s@xmath18 & @xmath60 k +   &   &   & high : 10 km s@xmath18 ( constant ) & @xmath69 k + photosphere & @xmath5312% & @xmath532% & assumed @xmath53 5% in luminosity & @xmath70 decreased from 0.45 to 0.42 +   & @xmath5315% & @xmath536% & assumed @xmath53 250 k in temperature & + calculation error@xmath71 & @xmath5330% &   &   &   + total & @xmath72 & @xmath73 &   &   +   +    the easiest way to see how these uncertainties affect a star s determined mass - loss rate is to propagate them through the dusty code . the results for v1 ( table [ mdoterrortable ] suggest that , in absolute terms , our mass - loss rates are little more than an educated guess . however , since a large number of other studies use the same assumptions , one can use these rates to trace differences between stars , providing these properties are largely invariant between stars . the parameters listed in table [ mdoterrortable ] have different effects . the grain size , for example , has little effect on the mass - loss rate but significant effects on other wind properties . conversely , the dust - to - gas ratio have significant effects on the mass - loss rate and velocity but little effect on the other wind parameters ( temperature , optical depth and minerology ) . clearly the largest effect on the mass - loss rate is whether the wind has an initial velocity driven by , for example , pulsation . determining the radial velocity profile of the wind would reduce the error substantially . while this has not been done accurately for such low - metallicity , dust - producing giants , inferences from chromospheric line modelling ( e.g.  @xcite ; @xcite ; @xcite ) appears to show that wind velocities do not decrease substantially from @xmath410 km s@xmath18 in the stages shortly prior to dust formation . we point out , however , that the highest mass - loss rates ( and by inference the highest wind velocities ) may be ruled out on evolutionary grounds (  [ dustcompsect ] ) . we emphasise that the errors derived above are not formal and are only estimates . these compare relatively well to the adopted error of @xmath45@xmath74 for silicate - dominated winds ( e.g.  @xcite ) , where the grain properties are somewhat better determined due to the long history of their study and constraints from their spectral features .    in relative terms , however , we can probably assume that the star - to - star variations are much smaller . in particular , the driving mechanism and hence the wind s initial velocity should not differ as markedly as the uncertainties we assume . likewise , the dust - to - gas ratio will only depend on the condensation fraction in the wind , which should not vary as much as a factor of 5@xmath74 between stars at the same metallicity . this reduces the overall error among stars in the same cluster to a factor of @xmath42@xmath74 and allows us to make direct comparisons between stars . the wind parameters listed in table [ mdottable ] pose several problems for stellar wind theory , which we discuss in this section . firstly , the projected terminal velocities we derive are very low , implying radiation pressure from the star is not effective in driving the wind : we must consider whether other processes may contribute to , or dominate , the wind acceleration . secondly , when one considers the evolutionary timescales of the stars in question , our derived mass - loss rates imply that the stars are losing much more material than is predicted from stellar theory . finally , we can use these problems to constrain some of the wind parameters listed in table [ mdoterrortable ] . mass - loss rate , pulsation amplitude and pulsation period are all known to correlate with stellar evolution . this is demonstrated in figures [ hrdvarsfig ] , [ plfig ] and [ mdotrelfig ] . these can be used to help solve the above problems .      a major problem when calculating mass - loss rates from metal - poor stars is that the terminal velocities predicted from radiative transfer codes are very low . this is a particular worry with dusty , as no correction is input for stellar gravity . in this situation , it is impossible not to have mass lost from a star when dust is present , whereas the real situation may be that radiation pressure on dust is insufficient for it to overcome the stellar gravity field . assuming the dust layers around these stars do have some outward velocity , a different method of wind acceleration may be required . dusty provides a minimum mass for a radiative transfer model . if a star is less massive than this , the error due to the unspecified correction for stellar gravity will be greater than the 30% inherent in the model . as this error increases to 100% , it is unclear whether radiation pressure is capable of driving the wind . this is a particular concern for the less - luminous stars with high mass - loss rates : x03 , v18 and v13 . while expansion velocities of 34 km s@xmath18 are not unreasonable when compared to halo carbon stars @xcite , the lower velocities we derive for those stars losing less mass become increasingly unreliable due to this error . @xcite predict that stars with metallicities below [ fe / h ] = 1 can not drive winds via radiation pressure , and must rely on acoustic ( pulsation ) driving to eject mass from the star ( see also @xcite ) . at [ fe / h ] = 0.7 , 47 tuc is near this hypothetical limit , meaning we must consider acceleration by pulsation and other sources of energy . the velocity amplitude of the photosphere in these variables ranges from 423 km s@xmath18 , with v1 , v2 and v3 showing line doubling in near - ir and/or line emission in the optical due to the propagation of shocks in the atmosphere ( table [ mdottable ] ; @xcite ; @xcite ) . while this velocity greatly exceeds our calculated wind velocity of @xmath44 km s@xmath18 , the bulk material does not necessarily attain this velocity . nevertheless , the presence of shocks in the extended atmosphere implies that energy is being transferred to larger radii , and implies that this energy is available for driving mass loss .    a comparison between the presence of circumstellar dust and pulsation reveals some interesting results . all the stars listed in table [ mdottable ] are variable : lee1424 = acvs002331 - 7222.6 , while x03 shows a variation of @xmath75 mag ( n.  matsunaga , private communication ) . conversely , 22 of 43 variable stars show at least probable evidence of ir excess . dust - producing stars therefore appear to be ( in 47 tuc ) exclusively variable , and variable stars tend to be dust producing . in both figures [ hrdvarsfig ] & [ plfig ] there appears no easy way of discerning rgb - pulsators from agb - pulsators : the outliers in figure [ hrdvarsfig ] are such due to crowding in the cluster core ( paper iii ) , while the pulsation periods in figure [ plfig ] do not show the bimodality that one might expect from having rgb and agb stars at different masses . however , a bimodality may exist in mass - loss rate with a transition occurring around 25003000 l@xmath2 ( figure [ mdotrelfig ] ; cf .  a similar possible bimodality in the gas - production rate in @xcite ) . the five stars with elevated mass - loss rates are also those with the highest pulsation velocity amplitudes and among those with the longest pulsation periods . as the presence of pulsation and its amplitude both appear to correlate with mass - loss rate , an obvious conclusion would be that pulsation energy is dissipated in the extended atmosphere and helps drive the wind , though we caution that the strongest variables are also the most - luminous stars with the highest capacity for driving a wind through radiation pressure .    further evidence that pulsation sets the mass - loss rate can be drawn from the apparent invariance of mass - loss over the narrow range of temperatures we probe ( figure [ mdotrelfig ] ) . the photospheric temperature sets the radius at which dust can condense and therefore ( if other parameters are invariant ) the density of the wind . higher - density winds lead to more - efficient dust formation , as the timescale for collision between dust - forming molecules is shorter . an invariance of mass - loss rate with temperature thus means that the formation of dust has little impact on the mass - loss rates present , providing further evidence that the mass - loss rate is ` hard - coded ' into the wind before the dust - formation radius . if this is true , the rle of radiation pressure on dust may simply be to affect the outflow velocity of the wind . if the wind is acoustically driven by pulsations , then the wind velocity and thus mass - loss rate is liable to be several times greater than the values listed in table [ mdottable ] , since these sources of energy are not included in our model . this poses potential problems , as mass - loss rates this high would lead to the dissipation of the star s atmosphere on much shorter timescales than is inferred from star counts . this problem may be solved if one invokes different parameters for dust in the wind , which could lead to mass - loss rates much smaller than the listed values ( see  [ dustcompsect ] ) .      in paper iii , we discussed the integrated rgb mass loss for a typical cluster star . we noted it was likely @xmath320.22 m@xmath2 and possibly , to first order , invariant with metallicity . this invariance would further suggest that the primary driver of rgb mass loss is not related to dust production . we had also estimated in paper iii , based on the cluster s luminosity function , that one star reaches the agb tip every 80 kyr .    these figures , however , do not correlate with our derived mass - loss rates . we find here that the cluster s @xmath422 dust - producing stars are losing a total of @xmath761.2 @xmath77 m@xmath2 yr@xmath18 of dust - containing wind . above the rgb tip ( @xmath762000 l@xmath2 ) , where dust production is ubiquitous , 13 stars are producing @xmath78 m@xmath2 yr@xmath18 . this implies that a staggering @xmath760.8 m@xmath2 is lost in the @xmath761 myr which a star spends on the agb above the rgb tip . clearly , this is impossible , as this is the entire mass of the star . this argues for a mass - loss rate of order 48 times lower than we measure . if additional wind driving mechanisms are important , this factor will increase further .    mass loss via dust - containing winds between 1000 and 2000 l@xmath2 can be similarly determined . here , nine out of the 56 stars in this luminosity range are producing @xmath42 @xmath79 m@xmath2 yr@xmath18 , implying 0.16 m@xmath2 is lost over a 4.5 myr timescale . if one assumes that mass loss from these stars is also lower than measured , and that some or all of these stars are agb stars , this implies that mass loss via non - dusty winds may be of similar or even greater importance in removing the estimated @xmath320.22 m@xmath2 from the star between the middle of the rgb and middle of the agb .    at the rates implied in table [ mdottable ] , the iron in the atmosphere of v1 would be entirely ejected within 1000040000 years , depending on the remaining envelope mass ( 0.080.3 m@xmath2 ) . v8 , v2 and v3 would similarly disperse the iron in their atmospheres within timescales of @xmath2590000 years . the padova and dartmouth isochrones used in paper iii predict that these stars ( especially v2 and v3 ) still have several hundred kyr left on the agb , implying that our iron mass - loss rates are ( again ) a factor of several too high . .maximum elemental abundances in dust [ cols= " < , > , > , > , > , > , > " , ]     in  [ windaccsect ] , we find radiation pressure may not be the dominant driver of winds from metal - poor stars , and is perhaps incapable of driving them in some cases . the simplest solution is to provide momentum from another source , such as stellar pulsations , in order to increase the wind velocity . unfortunately , this also has the effect of increasing the mass - loss rate , which is already higher than stellar evolutionary theory can accommodate (  [ dustevolsect ] ) . while we expect another momentum source to be present in the most - metal - poor stars ( simply because radiation pressure on dust is that inefficient ) , the limits implied by the mass - loss rates imply this is not the case in 47 tuc .    a solution to this dilemma should address both these problems by decreasing the mass - loss rate without significantly decreasing the wind velocity . examining table [ mdoterrortable ] , the simplest change would be to bring the dust - to - gas ratio closer to unity , which would decrease the mass - loss rate and increase the wind velocity . alternatively , we may consider measures that would increase the wind velocity without appreciable increasing the mass - loss rate , such as changing the size distribution of grains . the dust - to - gas ratio , @xmath36 , given in eq . ( [ psieqn ] ) , assumes that the dust fraction scales with metallicity ( cf . this will hold true if every dust species condenses with the same efficiency . navely , one may imagine that the fraction of each element condensed would decrease toward low metallicity simply because there is less opportunity for metal - metal reactions to take place in a metal - poor environment . from our conventional understanding , it is difficult to see how the dust - to - gas ratio in the wind could be substantially closer to parity . we can , however , calculate a maximum amount of dust per unit gas volume , by taking the major dust - forming elements ( o , c , fe , si , mg , ca , al ) plus some other potential dust - forming metals ( ni , ti , cr , mn ) and determining how much dust they can form . for 47 tuc , we assume the abundances of most elements are scaled from the solar values of @xcite by [ z / h ] = @xmath80 dex , following @xcite . we take the following exceptions , however : [ fe / h ] = 0.7 dex @xcite ; [ he / h ] @xmath76 0.02 dex @xcite ; [ o / fe ] = + 0.2 dex @xcite ; [ mg / fe ] @xmath76 [ si / fe ] @xmath76 [ ti / fe ] = + 0.23 dex and [ ca / fe ] = 0.0 dex @xcite . table [ psitable ] shows the stages of this calculation , by which we conclude that @xmath81 . here we make the simplistic assumptions that all dust - forming elements start out in the gas phase , that all carbon forms co , that all silicon forms enstatite ( mgsio@xmath20 ) , that iron and nickel condense entirely as pure metals , and that the leftover metals all form their respective oxides . chemical reactions must be calculated on the basis of number density , but these must be converted to mass densities in order to calculate a dust - to - gas ( mass ) ratio . for reference , the same calculation using solar metallicities arrives at @xmath82 . this approach neglects several points . carbon may be enhanced in these stars due to dredge - up 2 m@xmath2 ) , carbon - enhanced stars of @xmath40.8 m@xmath2 exist in globular clusters ( e.g.  @xcite ) . ] , but unless it is sufficiently abundant ( c / o @xmath5 0.75 ) to bind all the remaining free oxygen , this does not affect the end result . we also neglect water condensation , which could use up free oxygen , but this would serve to decrease @xmath36 . forsterite ( mg@xmath19sio@xmath83 ) may condense out in preference to enstatite , and calcium and aluminium may form into silicates before enstatite , but these would also decrease @xmath36 where they affect it at all . the formation of feo and iron - rich silicates would increase @xmath36 , but it seems that only a small percentage of iron condenses into non - metallic forms in these stars . our stated value of @xmath84 is therefore a relatively robust maximum value . table [ mdottable ] suggests that the number density of iron grains to other dust species in the winds is very high ( typically 10:1 , where it can be measured ) , but the analysis shown in table [ psitable ] would suggest it is of order 1:2 if all dust is condensed . this would suggest that the fraction of silicate dust that is condensed is quite low ( @xmath45% ) , however we have not yet accounted for two possibilities : that the temperatures of the grain populations may be limited , or that the size distributions ( specifically the surface - area - to - volume ratios ) of silicates and iron may be different . differing grain temperatures between dust species may exist when the wind is sparse enough that dust - gas collisions can not equalise the thermal energy . this would mean that more - opaque grains , like iron , are warmer ; more - transparent grains , like silicates are colder for a given radius , and could explain the difficulties we have in fitting both the metallic - iron continuum , and the 10- and 20-@xmath29 m silicate bands simultaneously (  [ mdotsect ] ) . more significantly , however , if thermal energy can not be equalised via collisions , this will lead to a significant velocity drift between the gas and dust components of the wind . this will substantially decrease the projected mass - loss rate , while the dust outflow velocity will also increase accordingly . if the gas outflow velocity is below the turbulent velocity in the wind , this may allow gas to preferentially fall back , making the star more metal poor . post - agb stars and planetary nebulae are known to be metal - poor for similar reasons ( e.g.  @xcite ) , but it is not clear that this process is happening while the stars are still on the agb where the outflow velocities are lower and the wind densities at the dust - formation radius are greater . conversely , if dust grain temperatures are equalised and the formation temperatures of silicates and iron are roughly similar , we may examine the effect of different grain sizes by taking two limiting cases . if we take the limiting case that both silicates and iron grains follow the same continuous distribution of ellipsoids , the silicate condensation fraction is @xmath45% , and we arrive at a value of @xmath85 to 1/3000 if all the iron is condensed this would increase our derived @xmath34 values by @xmath7660% while decreasing @xmath38 by @xmath7630% . this is the opposite effect to that needed to reconcile our observations with evolutionary theory , therefore we consider this scenario unlikely . in the converse limit , silicate and iron dust has condensed efficiently , but the iron grains have average radii 20@xmath74 smaller . this increases @xmath34 by a moderate ( 3060% ) amount , but increases @xmath38 by much more ( 100200% ) , depending on the exact size distribution . if one is prepared to accept that iron grains grow preferentially in one direction , e.g.  if magnetic effects are important , needle - like cylinders would be created ( cf . this would lead to a larger @xmath86 , but in this case without an increased @xmath34 . however , if identical fractions of silicates and iron are condensed , the length / radius ratio of the iron grains would need to be a staggering @xmath761200 , meaning that the grains must be smaller as well . for example , a population of grains 5@xmath74 smaller with a length / radius ratio of 10 allows a similar fraction of silicates and iron to be condensed with only a slight increase in @xmath34 . there is no obvious mechanism that would cause iron to be _ preferentially _ condensed at low metallicities , though magnetic effects and/or the basic problem of forming grains may be important ( an iron grain only needs two iron atoms to bond , whereas a silicate grain requires a series of chemical reactions to form mgsio@xmath20 ) . higher partial pressures of iron may also assist here , though it is not clear why this would happen in a metal - poor environment . note that in all cases above , we have assumed grains are homogeneous : our calculation of the wind parameters in  [ mdotsect ] assumes all grains are a uniform mix of the included species , while our discussion above takes the opposite limit : where grains of each species exist separately in pure form . in reality , grains aggregate and form inclusions inside each other , which complicates the above arguments . in general , however , the evolutionary constraints provided in  [ dustevolsect ] would suggest that dust condensation efficiency is high compared to solar - metallicity stars ( suggesting @xmath87 ) , and the iron grains are smaller and/or more - elongated than both the silicate grains and the assumed mrn distribution . a combination of these effects is likely required to reduce @xmath34 by the factor of 48 required without reducing the outflow velocities , @xmath86 , we derive ( table [ mdottable ] ) .      where we are able to determine it , the composition of circumstellar dust and its variation with luminosity is similar to that we see in other clusters in @xcite ( see also figure [ fehfig ] ) . the evolution of dust composition in 47 tuc broadly agrees with that seen in other clusters at this metallicity . metallic iron appears to be ubiquitous around dust - producing stars , while amorphous silicates are produced by most stars without any obvious dependance on luminosity . the published spectra of several stars ( v4 , v8 , v13 , and weakly in v21 ; see @xcite ) show sharp features around 9.7 , 11.3 , 13.1 , 20 and 32 @xmath29 m . @xcite attribute the 20-@xmath29 m feature to iron - rich wstite ( mg@xmath88fe@xmath89o ) , where @xmath90 ) . recent work by niyogi , speck & onaka ( submitted ) suggests that all these features , with the exception of the 13.1 @xmath29 m feature ( which is probably due to an aluminium oxide  @xcite ) , can be explained by iron - rich crystalline silicates . here , the 20-@xmath29 m feature is still produced by an fe - o bond , but it exists in a large crystalline lattice of fayalite ( fe@xmath19sio@xmath83 ) or ferrosilite ( fesio@xmath20 ) . iron - rich silicates are not expected to form in circumstellar environments , as their condensation temperatures are below that of both magnesium - rich silicates and metallic iron @xcite . if iron - bearing crystalline silicates are a significant component of some winds , they must exist alongside metallic iron and amorphous silicate species . one way this could occur is if cooling of the wind proceeds at a faster rate than dust condensation , which may occur in the tenuous winds of low - luminosity , metal - poor stars . figure [ fehfig ] shows that globular cluster stars which show a 20-@xmath29 m feature may be confined to luminosities below @xmath44000 l@xmath2 . it is difficult to say much about their metallicity dependance due to insufficient data , especially given silicate production appears to decline below [ fe / h ] @xmath4 1 . interestingly , the winds found in such stars ( x03 , v13 and v18 ) appear to be quite cool . in more - luminous stars , we find wind temperatures to be higher ( @xmath41000 k ) . the condensation of metallic iron before silicates remains a mystery : for pressures where iron condenses around @xmath41000 k , magnesium - rich silicates should have already formed closer to the star ( e.g.  @xcite , their figure 2 ) . significant departures from lte , or unusually - high partial pressures of iron ( caused , perhaps , by magnetic confinement or clumping in the wind ) may offer ways to circumvent this problem , but to invoke them as solutions would be speculative without further data . by whichever process metallic iron is created , if radiation pressure is effective in wind driving , metallic iron condensation should have a dramatic effect on the wind . its condensation should cause a rapid acceleration of the wind , dropping the density and preventing other dust species from condensing . this may be what we see in moderately - luminous stars , such as v3 , lw10 and v11 .    if radiation pressure is not an effective driver , the general absence of other dust species may be explained if silicates form using iron grains as seed nuclei . in this case , denser or more - metal - rich winds effectively condense silicates onto iron grains via collisions , while less - dense or metal - poor winds expand before much condensation occurs , leaving bare iron grains . this scenario may explain the iron grains found in the interplanetary dust particles known as gems ( glasses with embedded metal and sulphides ; @xcite ) .    in the most - evolved stars ( v1 , v8 , v2 and v4 ) , however , we do see silicate emission making an increasingly - strong presence . meanwhile , there is a hint that the temperature in the dust formation zone ( as traced by metallic iron condensation ) may be decreasing . as noted in  [ dustcondsect ] , different dust species may also have different temperatures . more detailed modelling of the winds of these stars , allowing for thermal decoupling of the grain species , is required if we are to speculate on factors which cause this temperature change . it could , however , be linked to stronger pulsation of the star and opacity shielding within the wind . in a cooler wind , however , silicates once again condense out in increasing preference to metallic iron , which contradicts our tentative finding that the silicate grains are cooler than the iron grains in these stars . in this work , we have used spectral energy distributions and infrared spectroscopy to establish mass - loss rates and investigate mineralogy around the dust - producing stars in 47 tuc . we summarise our conclusions as follows :        the total mass - loss rate is estimated to be @xmath41.2 @xmath74 10@xmath91 m@xmath2 yr@xmath18 , mainly in the form of metallic iron , but this is probably over - estimated by a factor of 48 due to unconstrained properties of the dust grains .    it remains unclear whether radiation pressure on dust is capable of driving a wind alone , due to the low projected terminal velocities of the stellar winds . the existence of circumstellar dust appears to correlate with the presence of stellar pulsation . we theorise that such pulsation may inject energy into the wind acoustically . variations in the wind chemistry broadly concur with those found in @xcite , with silicates found preferentially in high - luminosity , high - metallicity winds , while iron is found preferentially in low - luminosity , low - metallicity winds . iron may condense preferentially at low metallicity . alternatively , iron grains may be coated by silicates in dense , metal - rich winds , but this coating process may not be effective in winds around lower - luminosity and metallicity stars .    * acknowledgments : * we are grateful to thomas lebzelter for sharing his original data with us . this paper uses observations made using the _ spitzer space telescope _ ( operated by jpl , california institute of technology under nasa contract 1407 and supported by nasa through jpl ( contract number 1257184 ) ) ; observations using _ akari _ , a jaxa project with the participation of esa ; and data products from the two microns all sky survey , which is a joint project of the university of massachusetts and ipac / cit , funded by nasa and the nsf .                                                                        , s. 2007 , in astronomical society of the pacific conference series , vol . 378 , why galaxies care about agb stars : their importance as actors and probes , ed . f.  kerschbaum , c.  charbonnel , & r.  f. wing , 145
Abstract:

Answer: dust production among post - main - sequence stars is investigated in the galactic globular cluster 47 tucanae ( ngc 104 ) based on infrared photometry and spectroscopy . we identify metallic iron grains as the probable dominant opacity source in these winds . typical evolutionary timescales of agb stars suggest the mass - loss rates we report are too high . we suggest that this is because the iron grains are small or elongated and/or that iron condenses more efficiently than at solar metallicity . comparison to other works suggests metallic iron is observed to be more prevalent towards lower metallicities . the reasons for this are explored , but remain unclear . meanwhile , the luminosity at which dusty mass loss begins is largely invariant with metallicity , but its presence correlates strongly with long - period variability . this suggests that the winds of low - mass stars have a significant driver that is not radiation pressure , but may be acoustic driving by pulsations .