spiral disks form as a result of complex interactions of star formation , accretion , and gas .
these processes are recorded in a galaxy s stellar component .
the stellar metallicity , in particular , encodes information on a galaxy s overall chemical enrichment resulting from its history of star formation , gas accretion , and gas outflows .
spiral galaxies are commonly thought to form `` inside out '' with the inner regions forming earliest ( e.g. * ? ? ? * ) .
this effect has been observed by measuring a decreasing age with increasing radius in some nearby galaxies ( e.g. * ? ? ?
* ; * ? ? ?
a metallicity gradient is a natural consequence of inside out growth .
observationally , present - day metallicity gradients can be measured using hii regions or the atmospheres of young stars .
studies of the present - day metallicity gradients in the milky way and other nearby spirals do show lower metallicities at larger radii ( e.g. * ? ? ?
* ; * ? ? ?
* ) .    to measure metallicity gradients not at the present - day , but in the past , one must analyze stellar metallicities as a function of age .
the simultaneous measurement of age and metallicity is non - trivial , but can be done using either the observations of resolved stars or using integrated spectra .
recently , the califa survey @xcite used integral field spectroscopic data to study both the present - day and the mean stellar metallicity gradients in a large sample of galaxies @xcite .
the mean mass - weighted stellar metallicity gradient found among the califa galaxies is relatively shallow , @xmath3  dex / r@xmath4 , where r@xmath4 is the disk effective radius @xcite .
the galaxies in this califa subsample are very similar to m31 in hubble type and luminosity .
similar results based on photometric modeling were also found by @xcite for spiral galaxies in the virgo cluster . in both studies , there is no clear dependence on hubble type for spiral galaxies . in the milky way ,
the best gradients utilize measurements of 1000s of individual stars from large surveys ( e.g. sdss ) and found a gradient of @xmath5 - 0.06  dex / kpc @xcite .
in addition to the smooth gradients in metallicity , distinct metallicity structures are visible within galaxies .
metallicity substructure in the halos of galaxies have been linked to galactic merger events ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
within the disks of galaxies , structures with distinct metallicity signatures can arise due to bars , bulges and other dynamical structures @xcite .
the ability to resolve individual stars all at the same distance makes andromeda a unique object with which to study massive galaxy disks .
distance is the main source of uncertainty within the milky way and large - scale morphology observations are largely restricted to the solar neighborhood .
in contrast , stars can be resolved with ease at a wide range of radii in andromeda .
the panchromatic hubble andromeda treasury ( phat ; * ? ? ? * ) survey recently observed 117 million stars over a third of andromeda , providing a unique database for studying andromeda s stellar populations @xcite .    in this paper , we present a detailed analysis of the metallicities of andromeda s red giant branch ( rgb ) stars using phat photometry .
previous work on the rgb metallicities in andromeda have focused on the halo with limited information in the disk of the galaxy @xcite .    for this study
, we assume the following parameters : the center of m31 is at ra=10.68458and dec=1.2692@xcite ; m31 position angle and inclination angle are respectively 38and 74@xcite ; the distance modulus is 24.45@xmath60.05 mag or 776@xmath618 kpc ( mean distance , see @xcite for discussion ) . due to the small uncertainty in m31 s distance , and the primary dependence of the rgb metallicity on color and not luminosity , the distance uncertainty should have little effect on the analysis we present in this paper .
the foreground extinction is @xmath7 0.17 mag @xcite ; m31 s effective radius is 8.9@xmath60.8 kpc @xcite .
the paper is structured as follows : in section  [ sec : data ] , we explain our data sources , including the rgb star photometry and dust extinction measurements . in section  [ sec : principles - of - metallicities ] , we discuss our method for estimating metallicities from isochrone models and the uncertainties caused by photometric crowding and dust extinction .
we present and discuss our results on andromeda s metallicity gradient and a metal rich structure associated with the bar in section  [ sec : results ] .
the phat survey observed approximately 117 million stars with the hubble space telescope ( hst ) in ultraviolet ( f275w and f336w ) , optical ( f475w and f814w ) , and infrared ( f110w and f160w ) filters using the wfc3/uvis , acs / wfc and wfc3/ir hst camera , respectively .
the camera pointing angles divided the survey into 23  bricks each consisting of 18  overlapping ir - fields defined from the wfc3/ir camera footprint ( for details see @xcite ) .
3.4@xmath8 stars .
the red clump is visible at @xmath9 and @xmath10  24.3  mag .
the black dashed line indicates the rgb selection box used for this paper ( see section  [ sec : completeness - limit ] ) .
color indicates the metallicity .
the three color lines plot padova parsec1.2s isochrones at @xmath11=$]-2.0 , 0.0 , and 0.6  dex ( in violet , orange , and yellow , respectively ) at a fiducial age of 4  gyr .
the black star marks an example from our data sample at a deprojected radius of 6.1  kpc whose interpolated metallicity estimate is 0.033  dex ( see section  [ sec : metallicities ] ) .
the black arrow indicates 1  mag of dust extinction in @xmath12  band .
extinction causes stars metallicities to be overestimated as discussed at length in section  [ sec : dust - impact ] . ]
two versions of the photometry are publicly available .
first , there is camera - by - camera photometry derived separately on the pairs of filters in each camera producing a uv , optical , and infrared catalog that are then merged together to form a single six - filter catalogcamera - by - camera photometry can be downloaded at http://archive.stsci.edu/prepds/phat/datalist.html ] @xcite .
second , the `` six - filter photometry '' is derived across all filters simultaneouslysix filter photometry can be downloaded at http://archive.stsci.edu/missions/hlsp/phat/ ] @xcite . in both versions
, each star is parameterized by a position and magnitude , magnitude error , signal - to - noise ratio , crowding , and sharpness for each filter .
sharpness parameterizes the shape of each source relative to the point spread function ; and , crowding parameterizes the brightness of nearby sources .    in this paper
, we use the optical camera - by - camera photometry to derive rgb metallicity estimates .
the primary reason for using the camera - by - camera photometry over the six filter photometry is the availability of extensive artificial star tests ( also publicly available from mast archive ) .
these artificial star tests allow us to understand the photometric completeness and bias ( see section  [ sec : artificial - star - tests ] ) .
the optical color of rgb stars also is more sensitive to the metallicity relative to the ir photometry , while the uv photometry suffers from significant incompleteness for even the brightest rgb stars .
the smaller effect of dust extinction on the ir data makes it logical to include this data in our analysis , however , artificial stars are only available for the two - camera photometry in the ir , and this data is significantly shallower and crowding limited than the optical data ( see figure  11 of * ? ? ? * ) .
the six filter photometry is significantly deeper in the ir due to the ability to use the optical data to determine the positions of stars , but it is not computationally feasible to run artificial star tests over the full phat footprint for the six band photometry @xcite .
we create an optical only ( @xmath13 and @xmath14 ) photometric catalog with uniform coverage from the individual acs / wfc field images .
this involves two primary steps : defining non - overlapping spatial regions , and filling in the chip gaps of each individual field .
the boundary of each field is first defined using a grid of non - overlapping regions that roughly follow the individual wfc3/ir field footprints .
we additionally crop bricks 9 and 12 , which overlap with neighboring bricks , because of orientation differences described in section 3.5 of @xcite .
the phat tiling scheme allows us to use neighboring acs / wfc fields to fill in the chip gap of each acs / wfc field .
thus , our final catalog has uniform spatial coverage over the full phat survey region but does not make use of multiple overlapping observations .
this same catalog creation process is also used to deal with chip gaps in the the six - filter photometry as described in @xcite .
we apply quality cuts on our photometry .
the cuts use signal - to - noise , sharpness , and crowding in the f475w and f814w filters for each star .
we use the same parameterization which defines `` gst '' ( good star ) cuts from @xcite by applying all of the following equations :    @xmath15 + the variables correspond to the names given in the data files . for example , @xmath16 , @xmath17 , and @xmath18 respectively correspond to the signal - to - noise , sharpness , and crowding parameters for the @xmath13 filter .
we correct the magnitudes in @xmath13 and @xmath14 for foreground extinction using @xmath70.17  mag @xcite .
this extinction in @xmath19 is converted to @xmath13 and @xmath14 extinctions by multiplying with @xmath20 equal to 1.19119 and 0.60593 , respectively taken from http://stev.oapd.inaf.it/cgi-bin/cmd ] .
the corrected photometry is indicated by a subscript zero ( _ e.g. _ @xmath21 ) .
in addition to the `` gst '' cuts , we require that @xmath22 to remove stars with higher photometric bias and lower completeness .
the foreground reddening correction is applied before this brightness cut . after applying the spatial and photometric cuts
our sample contains @xmath57.0@xmath23 stars .    in figure
[ fig : data - selection - cmd ] , we present the optical color - magnitude diagram after quality cuts and foreground extinction corrections .
this color - magnitude diagram ( cmd ) shows the stellar density prior to the magnitude cut in gray scale contours .
the red - clump and bump of early - agb stars are both visible as peaks in the density below our selection region .
the dashed area shows the selection region in which we determine rgb metallicities . in this paper
, we estimate metallicities for rgb stars from their positions in the cmd .
isochrones at a range of metallicities are shown in figure  [ fig : data - selection - cmd ] , showing that stars of different metallicities are clearly separated in color ; we use this information to infer metallicities after considering a number of subtleties that can affect this measurement ( section  [ sec : principles - of - metallicities ] ) .
dust extinction has a significant impact on the observed photometry and thus our metallicity estimates .
we discuss this in more detail in section  [ sec : dust - impact ] . in this section
, we describe the source of our dust extinction measurements .
we use extinction measurements from @xcite ( hereafter ) .
uses a novel method to directly measure the dust extinction from the ir photometry ( filters @xmath24 and @xmath25 ) of phat rgb stars . in the ir cmd , unreddened rgb stars form a tight locus that is insensitive to age and metallicity .
dust extinction causes rgb stars to redden off this locus and form a second , broad rgb .
models the unreddened and reddened stars using three parameters : median extinction ( ) , spread of extinctions ( ) , and the fraction of reddened stars ( ) .
this latter parameter takes into account the geometry of the dust relative to the rgb stars ; 1- stars are assumed to be in front of the dust , and thus unaffected by it .
maps of each parameter are created with 2  pixels across the phat survey region . note that the central region of the phat survey ( brick 1 ) is not included in the dust maps due to extreme crowding .    in figure
[ fig : dust - maps ] , we present maps of the and used in this paper .
we rebinned and down - sampled the data into spatial bins used in our final results ( 0.01squares ) .
the width is approximately 0.3 everywhere and not shown because it has little bearing on interpreting the effect of extinction .
in contrast , the plays a significant role . increases in cause
the reddening to be more pronounced at the same mean extinction values .
for instance , in the case where @xmath26 and high , the median metallicity remains unchanged as stars in front of the thin dust layer are unaffected while those behind are reddened beyond selection bounds .
one feature visible in the map of is the increase towards the lower ( northwest ) edge of the survey . in this region
a higher fraction of stars appear behind the dust because of the disk inclination .
increasing means more stars will be extincted .
the final dust parameter , , is roughly 0.3 everywhere and is not shown here because its effect on our results is negligible . ]
in this section , we explain our method for deriving metallicity estimates and the sources of systematic and random uncertainty in these estimates . this section is organized as follows .    in section  [ sec : isochrone - models - and - age ]
, we discuss our choice of isochrone models , which we then use to interpolate the metallicity of individual stars based on their optical color and magnitude in section  [ sec : metallicities ] .
we then spatially bin the individual star metallicities to create a map of median metallicity . because of the large number of stars , the random sampling uncertainties are negligible .
we address systematic uncertainties in the median metallicity maps due to photometric bias and completeness ( [ sec : artificial - star - tests ] ) , dust extinction ( [ sec : dust - impact ] ) , and the uncertainties in the absolute metallicity scale due to isochrone choices and age ( [ sec : absolute - metallicities ] ) .
we correct the median metallicity map for the systematic effects of photometric bias , completeness and dust extinction and present results from this map in section  [ sec : results ] .
the intrinsic color of an rgb star depends primarily on its age and metallicity .
thus , given the photometry of a star and an assumed age , we can use isochrones to estimate a star s metallicity from its color and magnitude .
when estimating individual star metallicities ( section  [ sec : metallicities ] ) , we break the age - metallicity degeneracy using assumptions about the age of m31 s disk motivated by literature estimates . in this section ,
we discuss our choice of isochrone models and the possible impact of age on our metallicity estimates ; this impact is further explored throughout the paper .    for isochrone models , we use padova parsec1.2s isochrones @xcite . at ages
@xmath271  gyr , these isochrone models are available at metallicities ( @xmath28 ) between 0.0001 and 0.06 using @xmath29=0.0152 ) .
we use the obc bolometric corrections described in @xcite to convert the isochrones to photometric brightnesses in our optical bands .
we test dependency on the isochrone models by comparing models from parsec to basti @xcite and find the later has a typical offset of approximately + 0.3  dex .
this offset drops to @xmath50 for metal poor ( @xmath30  dex ) and metal rich ( @xmath31  dex ) stars .
because the majority of our stars have the typical offset , the resulting median metallicity map is only affected by a change in the absolute scale of the metallicity . as we discuss in section  [ sec : absolute - metallicities ] , absolute metallicity changes do not affect our results .
the isochrone models for the bright rgb stars have a well known degeneracy between age and metallicity @xcite . to break this degeneracy we adopt a single age motivated by the literature .
we assume that this age is either uniform throughout the disk , or that it varies linearly with radius .
current constraints on the age of m31 s disk are relatively minimal .
spectroscopic measurements of the age by @xcite show a mean age of @xmath512  gyr throughout the bulge , but this drops to @xmath58  gyr at their maximum radius of @xmath52  kpc ( 500  ) where the disk light starts to contribute significantly @xcite .
only weak constraints on the ancient star formation history ( sfh ) can be obtained from the phat data in the inner disk , but measurements by williams  et  al .
(  _ in prep _ ) suggest a mean sfh age of 10  gyr throughout the outer disk .
we can also get some insight into the age and age gradient present in m31 s disk by examining similar galaxies . as noted in the introduction ,
the subsample of galaxies in the califa survey analyzed by @xcite includes a large number of galaxies similar to m31 . taking the 25 galaxies within one magnitude of m31 in @xmath32 band luminosity ( califa magnitudes taken from
* and use @xmath33 for m31 ) and with hubble type sab to sbc , we find that the minimum mass - weighted age in the disks of these galaxies is 3  gyr , with typical values of @xmath57  gyr .
these measurements were made over similar physical scales ( 0.5 - 2 @xmath34 ) as the region we probe with the phat data .
these old mass - weighted ages of spiral disks are also consistent with previous work by @xcite .
while these previous studies suggest that the disk of m31 is old , the average age of stars observed along the upper rgb is not equivalent to the mass - weighted average age of a stellar system due to varying rgb lifetimes and the difference in masses of rgb stars with age .
the average age of stars along the upper rgb is normally _ younger _ than the mass - weighted average stellar age . by simulating a cmd using the star formation history derived by williams  et  al .
( _  in  prep _ )
, we calculate the mean age of rgb stars to be @xmath54  gyr throughout the disk .    for the primary median metallicity maps we assume a flat fiducial age of 4  gyr for our rgb stars . in section  [ sec : fiducial - age - effect ] , we show the metallicity gradient is not dependent on the exact fiducial age because the primary effect is a scaling of the absolute metallicity which leaves the gradient unchanged . in section
[ sec : radial - age - effect ] , we show that the median metallicity gradient can be changed if we assume a significant age gradient in the disk .
we estimate individual star metallicities by interpolating isochrones at the distance of andromeda ( @xmath35 ) in the @xmath36 vs @xmath21 color - magnitude diagram .
for interpolation , we use the linearndinterpolator function within the python package scipy .
this function finds the closest three isochrone points that create a triangle containing a star s color and magnitude using delaunay triangulation .
then it calculates the star s metallicity by linearly interpolating the metallicity at those three nearest isochrone points .
the isochrones used in this procedure are sampled at @xmath37 over the full range of metallicities considered , so the interpolation in metallicity is very small . in figure
[ fig : data - selection - cmd ] , we overplot the star metallicity on top of the stellar density contours in color and magnitude . in figure
[ fig : mdf - full ] , we present the metallicity distribution function ( mdf ) of our rgb stars using an fiducial age of 4  gyr .
we additionally show the mdf generated from stars in low extinction regions . from this figure
, we see that the dust has a small effect on the overall mdf of our sample .     mag ; see section  [ sec : dust - sources ] ) .
the median metallicity for both histograms is @xmath5 - 0.11  dex ( shown as a vertical line ) .
we note this mdf is not corrected for the systematic biases in the metallicity discussed in section [ sec : artificial - star - tests ] and [ sec : dust - impact ] , however , after correction , the mdfs look qualitatively very similar . ]    our primary goal in this paper is to examine the spatial variation in the metallicity distribution across andromeda s disk .
we present a spatial map of the median metallicity of the rgb in the top panel of figure  [ fig : uncorrected - med - metal ] .
we derive the median metallicity in 0.01(135  pc , projected ) square bins .
each bin contains between 100 and 19000 stars .
the upper limit becomes 3300 after we exclude the inner regions ( @xmath38  kpc ) because of completeness ( discussed in section  [ sec : completeness - limit ] ) .
the typical bin contains @xmath5560  stars .
in addition to the true variations in median rgb metallicity , this map shows features due to photometric bias , completeness , and dust .
we attempt to correct these effects in subsequent sections .    ) . _
bottom _ : model for the expected changes in the median metallicity due to dust extinction .
we interpret the increases and decreases in median metallicity in section  [ sec : dust - impact ] . for low extinction regions
the dust causes an increase of @xmath50.03  dex in the absolute metallicity and no significant spatial variations in the median metallicity . in our results ( figure  [ fig : corrected - median - metallicity - map ] ) , we subtract the two models from the raw median metallicity map to create a corrected version ( presented in figure  [ fig : corrected - median - metallicity - map ] ) .
note , all color bars have same range so relative differences can be assessed by eye .
]      photometric bias and completeness affect a star s observed photometry and thus the inferred metallicity estimate .
artificial stars allow us to evaluate and ultimately correct for these two systematics .
@xcite describes the optical camera artificial star tests in detail .
to summarize , @xmath39 artificial stars were distributed uniformly over each phat acs / wcs field of view , resulting in a total of @xmath54@xmath8 stars over the full survey .
the input @xmath40 colors and @xmath14 magnitudes for the artificial stars were split into two halves : the first half sample the full color and magnitude range of the data , while the other half sample the local observed color - magnitude diagram with an extrapolation to the faintest observed magnitudes .
each artificial star was inserted into and recovered from the data images independently using the dolphot package ( an updated version of the software presented in * ? ? ?
the same photometric cuts described in section  [ sec : data ] were applied to the output photometry .
we use these artificial stars to measure the photometric bias in our measurement in section  [ sec : bias ] , and to model the completeness in section  [ sec : completeness ] . then in section  [ sec : phot - corrections ] , we describe how we apply corrections for bias and completeness to our photometric metallicity estimates .      unresolved and bright stars in crowded regions cause a bias in the measured photometry that becomes increasingly large in more crowded regions .
this bias generally causes a star to appear brighter and with a color closer to the mean color of the overall galaxy due to blending with fainter undetected stars .    in figure
[ fig : cmd - bias - completeness ] , we show the average color and magnitude bias in bins across the cmd for a region of m31 at radii between 4 - 5  kpc where the bias is significant .
the bias of each individual artificial star is measured from the recovered magnitude minus the input magnitude . because the star s recovered cmd position
is changed from its original position , this photometric bias results in a bias in our metallicity estimates as well .
examining figure  [ fig : cmd - bias - completeness ] shows stars within our rgb region are typically detected to be brighter and bluer than their intrinsic cmd position .
our metallicity estimates for these stars are thus biased lower than the true metallicity , with the effect being more severe for high metallicity stars which are closer to the @xmath13 detection limit due to their intrinsically redder colors .    in section  [ sec : phot - corrections ] , we use the measured bias to partially correct for this effect . for computational efficiency
, we create a binned data set of the average bias .
we bin using 4 dimensions : two spatial dimensions with step size of 0.02(270  pc ) in both projected major and minor axis of m31 ; and , two color - magnitude dimensions with step size of 0.45 mag in both dimensions ( sizes for spatial and cmd bins seen in figures  [ fig : min - completeness ] and [ fig : cmd - bias - completeness ] respectively ) .
the result is a median of 3.3@xmath41 stars input for every spatial pixel and 36 input stars for every 4d bin . in each bin with
@xmath42 recovered artificial stars , we calculate the mean photometric bias for each filter from all the artificial stars .
stars in 4d  bins with @xmath43 artificial stars receive no correction value .
2 - 4  kpc , and therefore highlights the region most affected by completeness and bias that we consider in our analysis ( section  [ sec : completeness - limit ] ) .
our rgb selection box is indicated as a dotted line .
each bin is 0.45 mag on each side .
the completeness is the number of detected artificial stars divided by the number of those inserted in each cmd bin ( section  [ sec : completeness ] ) .
the arrows represent the magnitude and direction of bias for each bin ( section  [ sec : bias ] ) .
the bias generally shifts stars brighter and closer to the median color of the galaxy .
we correct our median metallicity determinations for completeness and bias as discussed in section  [ sec : phot - corrections ] . ]
the completeness is the number of detected stars divided by the total number of stars present at a given magnitude and color .
our completeness correlates with magnitude ; at fainter magnitudes we detect a smaller fraction of stars and are thus less complete .
because we require stars to be detected in both filters , each filter contributes to the completeness ; for blue stars the f814w magnitude determines the completeness ( thus creating a horizontal completeness limit on our cmd ) , while for redder stars , the f475w magnitude limits the completeness , creating a diagonal completeness limit on our cmd .
these completeness limits are visible in figure  [ fig : cmd - bias - completeness ] where the darkest pixels represent @xmath44 completeness .
the faint f475w magnitudes of metal - rich giants causes these stars to be most affected by completeness , and thus we detect a smaller fraction of the most metal - rich stars .
we measure the completeness from the artificial star tests as follows .
we define an artificial star as detected if it passes the photometry cuts defined in section  [ sec : data ] and has an absolute magnitude change less than 0.75  mag .
the completeness is then simply the fraction of artificial stars detected in a particular spatial and cmd volume .
we measure the completeness using the same 4d bins described previously in section  [ sec : bias ] .
to correct for completeness , we divide the total number of artificial stars recovered and inserted at a given cmd position using the _ input _ magnitudes ( i.e. after the bias correction , section  [ sec : phot - corrections ] ) .
this is somewhat different from the standard practice where completeness is calculated using the observed magnitudes to determine the detections and the input magnitudes to determine the total number of stars inserted at a given cmd position .    due to the small number of stars in each individual 4d bin , the errors on the completeness in an individual bin are large . to reduce these uncertainties , we model the completeness at a given color as a function of @xmath21 magnitude using a sigmoid function ( e.g. @xcite ) :    @xmath45
our modeling process is a loop with five steps : first , we select artificial stars within a 2-dimensional spatial bin .
second , we select artificial stars within a particular color bin and in the neighboring color bin(s ) . using neighboring bins
increases the number of stars used for each fit . also as a result the fit parameters become covariant .
this is beneficial because the covariance effectively smooths along the color dimension and makes the measurement more robust .
third , we bin the color - selected stars into independent @xmath21 magnitude bins .
the completeness for each 2-dimensional color - magnitude bin is the fraction of detected to input artificial stars .
an error in this measurement is also estimated .
fourth , we fit the completeness vs @xmath21 to equation  [ eq : completeness - model ] using @xmath46 minimization and determine @xmath47 and @xmath48 for that color bin .
fifth , we use the model to assign a fitted completeness to each 4d bin .    at the end of this process
, we have robust measurements of our completeness binned in 4d .
we show a spatial map of our completeness for rgb stars in figure  [ fig : min - completeness ] ; this figure is discussed further in section  [ sec : completeness - limit ] .      using the results from [ sec : bias ] and [ sec : completeness ] we correct our final median metallicity measurements for the photometric bias and completeness systematics
. to apply this correction , we first correct the observed photometry for photometric bias . then , when aggregating into the median metallicity , we correct for completeness by weighting each star by the inverse of the recovered completeness from section  [ sec : completeness ] . in this section , we further describe these corrections and the creation of a photometrically corrected median metallicity map .    to correct the observed star photometry for the photometric bias in each filter , we first bin the observed stars using the same 4d bins described in section  [ sec : bias ] .
then , we select each spatial bin and consider only the mean bias values in the associated cmd bins .
we interpolate the bias for each star in both filters using the cmd bin centers and mean bias values ( as we did for the metallicity in section  [ sec : metallicities ] ) .
finally , the observed star photometry is corrected for bias by subtracting interpolated bias from the observed photometry .
we use the bias corrected photometry to estimate a photometric completeness for each star .
the process is similar to interpolating the bias values for each star , only we now use the 4d recovered completeness data from section  [ sec : completeness ] .
we apply a minimum completeness limit of 0.05 so when weighting by completeness no star can dominate the calculations .
this limit affects @xmath49 of our stars .    for every spatial bin in our final maps ,
we calculate a bias and completeness corrected median metallicity .
we first re - derive metallicities for the stars using the bias corrected photometry . to calculate the median metallicity
, we weigh each star by the inverse of its completeness and then take the 50th percentile of the stars mdf .
the middle panel of figure  [ fig : uncorrected - med - metal ] shows the effect of bias and completeness on the median metallicity map ; specifically , it shows the median metallicity determined from the uncorrected photometry minus the median metallicity after we apply our bias and completeness correction .
the bias and completeness correction is largest near the center where the density of stars is high and is negligible in the outer parts of the galaxy .
we discuss the limits to these corrections in section  [ sec : completeness - limit ] .
this correction model is not exactly accurate due to uncertainties in the artificial star tests caused by stochasticity and assumptions in the prior distribution of the stars .
nevertheless , the model provides a reasonable and robust correction to this known systematic .      where the stellar density becomes high , photometric bias and completeness effects
become large . when the completeness drops to near 0%
, our bias and completeness corrections are no longer valid as they require the intrinsic distribution of stars to be at least partially detected .
this limits how close to andromeda s galactic center we can measure metallicities and limits our ability to detect the most metal - rich stars .
we use a minimum completeness metric to quantify exactly how close to the center our corrections are applicable . at each spatial position
, we measure the completeness along the 90th percentile density contour of our rgb stars ( shown as a black contour line on figure  [ fig : data - selection - cmd ] ) limited our rgb @xmath14 magnitude limit .
we use the minimum completeness along the contour to create the map shown in figure  [ fig : min - completeness ] . for almost all spatial positions
the minimum completeness is measured at the most red rgb tip .
this means the completeness along our @xmath14 magnitude limit is greater than this minimum completeness metric at practically all spatial positions .    from the minimum completeness map ( figure  [ fig : min - completeness ] ) , we define the limiting central region at the 50% contour which is @xmath54  kpc . at this limit ,
the bias and completeness correction is @xmath50.1  dex .
we exclude data within from our final analysis . beyond the limit
, we find the same correction has a @xmath50 effect on our metallicity graident ( section  [ sec : metallicity - gradient ] ) and no bearing on our observed metal rich bar structure ( section  [ sec : metal - rich - bar ] ) . we conclude this limiting region to be a reasonable choice . in all subsequent figures ,
we include the results at all radii but shade the excluded central region .    ) .
this metric is calculated by assigning a completeness to all the stars then taking the minimum from the 90th percentile  in cmd density  of our data .
the solid lines indicate 50% , 75% , and 90% minimum completeness . in general
, the stars most affected by completeness are the reddest , highest metallicity portion of the rgb . at the center of andromeda
these stars are completely undetectable .
as the sample becomes more incomplete it becomes impossible to correct properly for the bias and completeness ; we therefore only analyze the area with @xmath51% completeness ( see section  [ sec : completeness - limit ] ) . ]
dust extinction causes stars to appear dimmer and redder . for rgb stars ,
this causes us to overestimate the metallicity of reddened stars ( see arrow in figure  [ fig : data - selection - cmd ] ) .
we use information on the dust distribution ( see section  [ sec : dust - sources ] ) to both identify high extinction regions and to model the effect of dust extinction on the median metallicity .
to identify regions where the extinction is likely to cause biased metallicity estimates , we multiply the mean extinction by the fraction of reddened stars ( ) from within each spatial bin .
we find  @xmath52  mag is a good limit for identifying regions where extinction significantly biases our results .
we measure the significance of dust extinction on our median metallicity measurements by forward modeling the effect of dust extinction on an unreddened fiducial cmd .
to forward model our data , we first need a realistic input cmd .
we take this input cmd from unreddened regions in the survey . for the unreddened fiducial cmd
, we use stars in the low extinction region ( @xmath53  mag ) at @xmath5  14  kpc along the major axis ( visible in figure  [ fig : dust - maps ] ) .
our final result should be insensitive to low level extinction present in our fiducial cmd because we measure the relative change in median metallicity of the sample .
our results are very similar using other low extinction regions within the survey ( e.g. along the ne edge ) .
the primary weakness of our approach is that we do not account for the stellar population differences ( i.e.  young stars ) we expect between the dust - free and dusty regions of the galaxy ; trying to account for this change would require full star formation history modeling and is beyond the scope of this paper .
we expect that any young contaminants not included in our fiducial cmd would be brighter bluer stars reddened into our cmd region , and thus would likely reduce the median metallicity at moderate extinctions .
we apply reddening to this cmd to measure the resulting change in median metallicity .
the extinction for each individual star in the cmd is drawn from the extinction probability function from using their maps of the best fit , and dust parameters ( see figure  [ fig : dust - maps ] ) .
the stars are binned into 0.01spatial bins which correspond to the bin size used in our median metallicity maps . in each bin , we take the associated dust parameters and determine a random extinction for each star . then we apply that extinction to create a new reddened data set and re - evaluate the metallicity as described in section  [ sec : metallicities ] .
finally , we measure the change in median metallicity between new reddened data and the original unreddened sample .    in figure  [
fig : dmedian - vs - dust ] , we present the change in median metallicity as a function of for each spatial bin colored by .
these two parameters dominate the predicted change in median metallicity estimates ; changing has a @xmath54  dex change in the median metallicity estimates . as previously noted , for a particular the change in median metallicity increases with increasing .
lines of constant ( and constant = 0.3 ) are also plotted on figure  [ fig : dmedian - vs - dust ] . for these constant lines ,
the change in median metallicity first increases then decreases with increasing .
the increases are easily understood as stars becoming reddened and appearing more metal rich .
the decrease is caused by stars becoming so reddened they are moved beyond our rgb selection box and the most metal rich isochrone ; the remaining unreddened stars ( which are in front of the dust distribution ) yield a very similar median metallicity to the input cmd .
the individual points in figure  [ fig : dmedian - vs - dust ] correspond to individual spatial bins .
we also present these data as a map of the change in median metallicity ( lower panel of figure  [ fig : uncorrected - med - metal ] ) .
this map shows sharp changes due to dust especially along the upper ( northeastern ) edge of the survey and in the 10  kpc ring .
our final results avoid these regions by imposing a cut to remove these high extinction regions ( @xmath52  mag ) .    in the low extinction regions ,
the dust still has a low level effect .
we apply a dust correction by subtracting the change in median metallicity due to dust .
the changes in metallicity are small ( typically @xmath55  dex ) and  more importantly  distributed randomly over the survey area so the result is only a systematic on the absolute metallicity scale and not on spatial variations .
systematics on the absolute metallicity scale do not heavily impact our results .
[ m / h ] ) for each spatial bin in the data .
the simulation uses a fiducial cmd taken from a low - dust region and applies reddening to individual stars using the dust parameters ( , , and ) from the dust maps for every spatial bin .
the solid lines show changes in median metallicity at constant values of 0.1 , 0.2 , 0.3 , and 0.4 .
along one of these solid lines ( and assuming = 0.3 ) .
the dust increases the median metallicity until @xmath560.5 because the reddened stars receive an overestimated metallicity .
the median metallicity decreases for high extinction when the reddened stars are moved beyond the edge of our rgb box and no longer affect the measurement .
the points below the dashed line are low extinction regions with @xmath57 .
these are the regions used in our metallicity gradient measurement .
their @xmath58[m / h ] is @xmath59  dex , with typical values of 0.03  dex .
details on the dust simulations are given in section  [ sec : dust - impact ] ; a map of @xmath58[m / h ] is shown in the bottom panel of figure  [ fig : uncorrected - med - metal ] . ]
the absolute metallicities we determine have large uncertainties which we discuss in this section .    _ age : _ older populations have redder rgbs than younger populations of the same metallicity .
therefore , using older fiducial ages results in lower metallicity estimates .
this effect is quite significant , with metallicity estimates varying by @xmath50.4 dex when changing the age from 3 to 13  gyr ( we discuss this further in section  [ sec : fiducial - age - effect ] ) .
overall , this is the dominant uncertainty in the absolute metallicity . because we are primarily interested in the metallicity gradient
, an age gradient can significantly affect our result ; this is discussed further in section  [ sec : radial - age - effect ] .    _
isochrones : _ as noted in section  [ sec : isochrone - models - and - age ] , use of the basti isochrones results in a median metallicity that is @xmath50.3 dex higher than the parsec  1.2s isochrones we adopt .
_ dust : _ from our simulations in the previous section , we find the low level dust extinction present throughout the survey results in a slight overestimate of the intrinsic metallicity ; however this bias is only @xmath50.03  dex .
_ photometric bias and completeness : _ near the center , the bias and completeness correction raise our median metallicity by more than 0.1  dex .
unlike the effects above , this is primarily a bias in the relative metallicity between the center and outskirts , as the bias and completeness have no affect on the median metallicity in the outer disk regions .    because of the large uncertainties in the absolute metallicities , we choose to focus on relative measurements for our results . particularly , the metallicity gradient of andromeda and metallicity substructure of the disk . after applying corrections for known systematics , we present the median metallicity map and results in the next section .
in this section , we analyze the radial and spatial variation of the median rgb metallicity . our results are based on the spatial map of corrected median metallicities presented in figure  [ fig : corrected - median - metallicity - map ] .
our results focus on two interesting features in the corrected median metallicity map .
the first is an overall gradient towards lower metallicities at larger radii ( section  [ sec : metallicity - gradient ] ) .
the second feature is a region of high metallicity at radii of 3 - 7  kpc lying just above ( southeast of ) the major axis , which we associate with an elongated bar ( section  [ sec : metal - rich - bar ] ) .    )
was divided into annular bins with @xmath60 ; within these bins the mean of the map is shown as the solid line , while the shaded region shows the standard deviation of the pixels in the bin .
high extinction regions ( @xmath61 ) were excluded .
the red line shows the gradient including corrections for dust extinction and photometric bias and completeness , while in cyan the metallicities have only been corrected for dust extinction .
the bias and completeness correction changes the gradient by @xmath51@xmath2 , while the dust extinction correction changes the gradient by less than a 0.001@xmath2 .
the black line and shaded region shows the best linear fit at deprojected radii @xmath62  kpc  ( @xmath63 ; @xcite ) to the corrected metallicities , the best - fit slope is -0.020@xmath60.004 dex / kpc . the black shaded region at @xmath64  kpc where our rgb data becomes incomplete at high metallicities . ]
we present the first measurement of the rgb metallicity gradient in the m31 disk from 4 to 20  kpc .
we binned the median metallicity ( shown in figure  [ fig : corrected - median - metallicity - map ] ) into elliptical apertures with a size of 1  kpc along the major axis , and assuming @xmath65 @xcite to account for the disk s inclination .
we exclude the high extinction regions with @xmath61 to limit the effect of dust on our measurement ; this limit is shown as the white contour in the bottom panel of figure  [ fig : corrected - median - metallicity - map ] . after excluding these regions , we take the average and standard deviation of the median metallicity within each elliptical aperture to create the radial metallicity profile shown in figure  [ fig : feh - vs - radius ] .
we fit a gradient to the measured metallicity profile at @xmath62  kpc  ( @xmath63 ; @xcite ) .
this radial cut excludes the regions we distrust due to low completeness ( shown as the shaded region in figures  [ fig : corrected - median - metallicity - map ] and [ fig : feh - vs - radius ] ) .
we fit the mean and standard deviation of the metallicity in each radial bin using a monte carlo simulation with a linear fit . for each iteration
, we calculate the best fit slope and intercept using least squares optimization .
we repeat this and obtain a distribution of fitted slopes and intercepts ; this procedure also allows us to examine the covariance of the slope and intercept .
we show the best fit parameters and 1@xmath2 contours in figure  [ fig : fiducial - age - params ] for various fiducial ages including our primary result with 4  gyr .
we measure andromeda s metallicity gradient to be @xmath0  dex / kpc .
this gradient uses a 4  gyr fiducial age ( see section  [ sec : isochrone - models - and - age ] , fig  [ fig : fiducial - age - params ] ) and our corrected metallicities . the error on the gradient
is calculated using the monte carlo technique described above ; bootstrapping errors are significantly smaller than this quoted value .
the photometric bias , completeness and dust corrections have little impact on producing the observed gradient : the dust corrections change the gradient by @xmath66 , and the photometric bias and completeness corrections steepen the relation somewhat , but still agree at about the 1@xmath2 level .
we also measure a consistent metallicity gradient ( within 0.25@xmath2 ) when selecting stars just along the major axis which indicates the assumed inclination angle has little impact on the gradient .
the intercept ( central metallicity ) of this fit is [ m / h]@xmath67 ; though caution should be taken in interpreting the absolute metallicity as explained in section  [ sec : absolute - metallicities ] .
to facilitate comparisons with the literature , we also present our gradient with different units . in terms of andromeda s half light radius ( 8.9@xmath60.8  kpc ; @xcite ) ,
the gradient is @xmath68  dex/@xmath69 . using the radius of the isophote at a bolometric b - band of 25  mag ( @xmath70  kpc ;
@xcite ) , the gradient is @xmath71  dex / r25 .      as discussed in section  [ sec : absolute - metallicities ] , the assumed age influences the estimated metallicity of our stars .
we evaluate the sensitivity of our metallicity gradient results to changes in the fiducial age by repeating our analysis using isochrone models with ages of 2 , 4 , 6 , 8 , 10 , and 12  gyr . for each age
, we calculate a distribution of metallicity gradients and intercepts ; in figure  [ fig : fiducial - age - params ] , we present the best fit and 1@xmath2 contour of the distribution .
the metallicity gradient across all fiducial ages is consistent within 0.25@xmath2 .
this is because the rgb of isochrones at different ages have similar shapes in the cmd , with the primary difference being that older ages have a lower metallicity at the same rgb color .
this results in a significant change in the normalization of the metallicity gradient and the fitted intercept in metallicity , but little change in the slope .
we find a change in the intercept of roughly @xmath72  dex / gyr .
confidence after performing a monte carlo simulation to evaluate the slope and intercept .
our primary result presented in section  [ sec : metallicity - gradient ] uses a fiducial age of 4  gyr shown here as a diamond .
the age gradient for all fiducial ages are within 0.25@xmath2 .
changing the fiducial age does have an effect on the absolute metallicity which we discuss in section  [ sec : fiducial - age - effect ] .
this figure also demonstrates a covariance between the fitted slope and intercept of the metallicity gradient . ]
a gradient in the age has a significant impact on our metallicity gradient .
we evaluate this impact using the following experiment ; we fix the age of stars at the center of andromeda to 4  gyr .
then , we calculate the age for each star using its spatial position for age gradients of -0.35 , -0.25 , -0.1 , 0.0 , and 0.1  gyr / kpc . we calculate the metallicities similar to in section  [ sec : metallicities ] with age in addition to color and magnitude . in this application
, the delaunay triangulation finds the nearest four color , magnitude , and age isochrone points which contain a particular star s parameters and performs the linear interpolation of the metallicity between the points . using these new metallicities ,
we recalculate the median metallicity map and derive the distribution of best fit metallicity gradients and intercepts as described in section  [ sec : metallicity - gradient ] .    in figure
[ fig : feh - age - vs - radius ] , we present the results of this age gradient experiment .
as mentioned , as the age gradient becomes negative , the metallicity gradient flattens , while a positive age gradient steepens the gradient .
for an age gradient of -0.35  gyr / kpc the metallicity gradient is consistent with being flat .
this figure can be combined with figure  [ fig : fiducial - age - params ] to determine the inferred age gradient and intercept for a range of central ages and age gradients .
for instance , if the mean age of stars in our rgb box is 8  gyr in the center changing by -0.1  gyr / kpc , then our measured metallicity gradient would be @xmath5 - 0.018  dex / kpc with a central metallicity of @xmath5 - 0.04  dex .
we expect any age gradient in m31 to be small .
this is based on results from the phat data ( williams et al . _ in prep _ ) . also , @xcite find flat mass weighted age - gradients ( @xmath73dex/@xmath34 ) for a majority of galaxies in the califa sample , with the small age gradients being both positive and negative .
this translates to an age gradient less than 0.1  gyr / kpc using values for m31 and our fiducial age .
however , resolved star formation histories of m33 and ngc 300 ( lower mass , later - type disk galaxies than m31 ) have shown steeper negative age gradients , with younger ages at larger radii @xcite .
if such gradients exist in m31 , the younger stars at larger radii will cause a flatter metallicity gradient than using a single fiducial age at all radii .
however , resolved star formation histories of m33 and ngc 300 ( lower mass , later - type disk galaxies than m31 ) have shown steeper negative age gradients , with younger ages at larger radii @xcite .
if such gradients exist in m31 , the younger stars at larger radii will cause a flatter metallicity gradient than using a single fiducial age at all radii .
if any new information becomes available , figure  [ fig : feh - age - vs - radius ] can be used to update our metallicity gradient result based on the determined age gradient .    , but with a variation of stellar age with radius .
the age gradients shown have slopes of -0.35 , -0.25 , -0.1 , 0.0 , and 0.1  gyr / kpc with a central age of 10  gyr .
the metallicity gradient flattens as the age gradient becomes more negative with it being consistently 0
dex / kpc when the age gradient is -0.35
gyr / kpc . based on observations of external galaxies similar to andromeda
, we expect the age gradient to be @xmath74  gyr / kpc ( see section  [ sec : comparison - metallicity - gradient ] ) . ]      in this section , we compare our measured metallicity gradient to other metallicity gradient measurements in m31 and to metallicity gradients measured in other galaxies .
we find some evidence that the metallicity gradient may become flatter at large radii , but does not appear to vary with age .
the metallicity gradient of the young population of andromeda has been measured by @xcite using hii region abundances over the same radial region covered by the phat data .
they find a gradient of @xmath75  dex / kpc .
this gradient is consistent with our rgb metallicity gradient , suggesting that the age gradient has been stable in m31 for some time .    at larger radii , @xcite used abundance measurements of planetary nebula to measure a gradient with deprojected radii of 18 to 43  kpc .
these planetary nebulae have typical ages of 1 - 2  gyr , thus probing slightly younger stars than the rgb stars probed here .
they find a gradient of @xmath76  dex / kpc , significantly flatter than our gradient .
this may be evidence that the metallicity gradient flattens at larger radii .
another , more indirect approach to measuring metallicities is the stellar population modeling of integrated photometry of m31 from the androids survey @xcite .
their preliminary results seem to show a somewhat shallower gradient of @xmath5 - 0.01  dex / kpc out to @xmath530  kpc , with a possible steepening at even larger radii .
these comparisons provide conflicting support for the expectations due to radial migration ( e.g. * ? ? ?
on one hand , the shallower metallicity gradient in old populations at large radii could be explained if an increasing fraction of stars at large radii come from smaller radii .
however , radial migration also predicts that any metallicity gradient should get shallower over time , thus the similarity of the metallicity gradients in the young and old populations in the inner disk suggests that either the metallicity gradient was steeper in the past , or significant migration has not occurred in the inner disk .
higher quality stellar population modeling can be done with spectra , but published results include only the central bulge region .
spectral index measurements of m31 s bulge by @xcite show the metallicity at radii between 100  pc and 2  kpc is flat and approximately solar . at smaller radii ,
the metallicity rises sharply , with [ m / h]@xmath50.5 in the center . at the outer most radii
, these measurements include a substantial disk component .
our central metallicity extrapolated inwards to 2  kpc is roughly consistent with this measurement though the absolute metallicity scale in our measurement is sensitive to many systematics ( figure  [ fig : fiducial - age - params ] and section  [ sec : absolute - metallicities ] ) .
the recent califa survey observed @xmath77 face - on spiral galaxies using wide - field integrated field units ( ifu ) .
andromeda s size , morphology , and central velocity dispersion are typical among this sample . from the califa sample
, @xcite selected 62 face on galaxies and measured their stellar metallicities and ages .
the mean mass - weighted metallicity gradient of these galaxies is @xmath78  dex/@xmath69 or @xmath79  dex / kpc when scaled using m31 s effective radius ( @xmath80  kpc ; @xcite ) .
our metallicity gradient is slightly steeper than this mean metallicity gradient , but is well within the scatter observed of their galaxies .    )
are highlighted in black .
our m31 metallicity gradient ( assuming no age gradient ) is shown as the blue star .
dashed lines give the median of the full sample ( gray ) , and the m31-like sample ( black ) along both axes . ]
the metallicity gradient in the milky way ( mw ) disk appears to be significantly steeper than the gradient in m31 .
@xcite measure the metallicity gradient as a function of height ( z ) from the disk for segue survey stars . for stars near the @xmath81 they measure @xmath82  dex / kpc which drops to @xmath83  dex / kpc for @xmath84  kpc
( see figure  8 of @xcite for metallicity gradient vs z height ) . given
that the disk mass is dominated by stars below this height , the mw gradient does appear to be steeper than the m31 gradient .
a similarly steep mid - plane gradient is found in other survey data as well including apogee data @xcite , rave data @xcite and the gaia - eso survey @xcite .
the shallower gradient in m31 relative to the milky way is also discussed by @xcite .
we observe a region of high median metallicity likely associated with the end of the bar . in figure
[ fig : corrected - median - metallicity - map ] , this non - axisymmetric region lies just above the major axis between radii of @xmath53 and 7  kpc .
this region contains the highest metallicities found anywhere in our maps .
examination of contours of rgb density in the middle panel of figure  [ fig : corrected - median - metallicity - map ] shows a similar axial asymmetry .
this asymmetry coincides with m31 s bar as measured by @xcite at 8above the major axis .
@xcite also noted an overdensity in agb star counts near the end of the bar .
in addition to the presence of a stellar overdensity , a peculiarly high velocity dispersion ( @xmath5150  km / s ) is observed at the outer edge of this region @xcite .
this high dispersion could be explained if this represents the end of a misaligned bar .
this high metallicity region is not caused by dust but represents a real change in the stellar population . as previously noted in section  [ sec : dust - impact ] , dust causes the median metallicity to appear more metal rich for moderate extinctions of @xmath85  mag .
we observe the high metallicity region prior to the dust correction ( figure  [ fig : uncorrected - med - metal ] ) and it becomes even clearer after the correction .
even more noteworthy , the high metallicity is seen even in the very low extinction regions ( dust contours shown in figure  [ fig : corrected - median - metallicity - map ] ) .
this higher median metallicity region may be expected if either : ( 1 ) gas collects at the bar end resulting in higher enrichment there @xcite ; or , ( 2 ) the bar scatters stars from the inner part of the galaxy which have higher metallicity .
an additional possibility is there could be a shift in the rgb age in the bar region ; however , to get a high metallicity region , the average age would need to drop in the bar relative to its surroundings .
regardless of the mechanism for metal enrichment , this metallicity adds further evidence for a very long bar in m31 .
we present an analysis of the spatial variations in the metallicity of rgb stars in andromeda s disk using the phat survey .
this study is the first large study of rgb metallicities within the andromeda disk and covers a radial region between 4 and 20  kpc .
metallicities are estimated for 7.0@xmath23 rgb stars using their optical photometry and isochrone models .
we create a map of median metallicity and correct this map for the effects of dust extinction and photometric bias and completeness .
our final results are based on the median metallicity map after correcting for these systematics ( figure  [ fig : corrected - median - metallicity - map ] ) . from the corrected median metallicity map
we conclude : +    1 .   over the inner disk phat survey ( @xmath54 - 20  kpc )
, we measure a clear gradient in the median metallicity .
the measured gradient is [ m / h]=@xmath0  dex / kpc assuming a constant fiducial age of 4  gyr ( section  [ sec : metallicity - gradient ] ) .
this gradient has a less than 1  sigma sensitivity to systematic effects of dust extinction , photometric bias and completeness , and fiducial age changes .
the derived metallicity gradient is sensitive to any age gradient that may be present in andromeda , with negative age gradients resulting in shallower metallicity gradients .
a steep age gradient of -0.35  gyr / kpc would result in no metallicity gradient ; however , spectroscopic observations of other similar galaxies typically show much shallower age gradients @xcite .
we measure an enhanced metallicity region offset from the major axis .
we show this region is coincident with an overdensity of rgb stars and a region of enhanced dispersion @xcite . in section  [ sec : metal - rich - bar ] , we discuss how this likely is due to andromeda s bar @xcite .
we thank the anonymous referee for suggestions that improved the paper .
this work was supported by the space telescope science institute through go-12055 .
this research made use of codes in the scipy stack , as well as astropy , a community - developed core python package for astronomy ( astropy collaboration , 2013 ) .                                  , g. , de vaucouleurs , a. , corwin , jr .
, h.  g. , et  al .
1991 , third reference catalogue of bright galaxies .
volume i : explanations and references .
volume ii : data for galaxies between 0@xmath86 and 12@xmath86 .
volume iii : data for galaxies between 12@xmath86 and 24@xmath86 .