a long gamma - ray burst ( grb ) is one of the most energetic explosions in the universe , which is observed via soft gamma - ray emission followed by afterglow in lower energy range .
it is now broadly agreed that at least some of long grbs originate in core - collapse of massive stars together with core - collapse supernovae ( cc sne ) . however , cc sne do not always accompany long grbs , and the criteria discriminating between long grbs and general cc sne is one of the most important questions regarding long grbs .
some theoretical studies on the origin of long grbs using stellar evolution models suggest that low - metallicity may be a necessary condition for a long grb to occur ( @xmath3 , e.g. @xcite ) .
metallicity distribution of long grb host galaxies provides us with an important clue to study the relation between metallicity and a long grb occurrence .
although metallicity of a long grb host galaxy is not necessarily identical to that of the progenitor ( e.g. , @xcite ) , galaxies with low - metallicity would have higher grb production efficiency [ i.e. , long grb rate to star formation rate ( sfr ) ratio ] than high - metal galaxies if low - metallicity is a necessary condition for a long grb to occur .
with a sample of 5 host galaxies ( 3 long grbs and 2 x - ray flashes , xrfs ) , ( , hereafter s06 ) showed that the metallicity distribution of the host galaxies at redshifts @xmath4 is significantly biased towards low - metallicities compared with star forming galaxies at similar redshifts , suggesting long grbs really occur in low - metallicity environment .
however their sample is too small to determine the relation between metallicity and long grb occurrence rate .
the sample number of spectroscopically studied long grb host galaxies has been dramatically increased during the 10 years after s06 ( e.g. , @xcite ) .
however , the relation between metallicity and long grb occurrence rate is not understood quantitatively .
one difficulty is that we can determine redshifts and/or identify host galaxies of only a small fraction of long grbs ( @xmath1 30% in recent observations ) , which is likely biased with respect to the distributions of redshifts and host galaxy properties ( the @xmath5-determination / host - identification bias ) .
furthermore , even when the redshift and the host galaxy of a long grb are known , the host galaxy is not always studied in detail .
the targets of the spectroscopic studies have been selected non - uniformly by numerous independent researchers making the effect of the sampling bias unevaluable ( the reporting bias , e.g. , @xcite ) .
there have been some previous long grb host observation campaigns invented to overcome the sampling biases ( so called unbiased surveys , @xcite ) . in the unbiased surveys ,
the samples are selected without using informations about redshifts and host galaxies to avoid the @xmath5-determination / host - identification bias , and hence the samples span a very wide range of redshifts ( @xmath6 0-6 ) . in such a wide range of redshift , properties of general galaxies largely evolve making it difficult to reach enough statistics to study the environmental requirements for a long grb occurrence by comparing the properties of grb host galaxies to those of general galaxies .
some recent studies pick out grb host galaxies in specific ranges of redshifts from the unbiased surveys and compared their properties to those of general galaxies at similar redshifts @xcite .
however , the properties of general galaxies , to which we compare long grb host galaxies to investigate the environmental requirements for a long grb occurrence , are still not well understood at high - redshifts .
metallicity measurement of high - redshift galaxies is an especially challenging issue ( at @xmath7 , e.g. , @xcite ) .    on the other hand ,
the sample number becomes very small at low - redshifts in the unbiased surveys ( @xmath8 grbs at @xmath9 ) because the cosmic long grb rate density is low and the comoving volume element is small , although long grb host galaxies at low - redshifts are of special importance to understand the relation between metallicity and long grb occurrence rate as we describe below .
hence , although the unbiased surveys are important to understand what bias we would suffer by selecting long grbs whose redshifts and host galaxies are known , it is difficult to unveil the environmental requirements for a long grb occurrence only with the unbiased surveys .
metallicity measurements of the host galaxies of all long grbs known at low - redshifts ( @xmath10 ) can dramatically improve our understanding of the relation between metallicity and long grb occurrence rate .
there are advantages of studying grb host galaxies at low - redshift in some aspects .    1 .   the success rate of grb redshift - determinations / host - identifications would be higher at lower redshifts ( i.e.
, the @xmath5-determination / host - identification bias is weaker ) , because faint galaxies can be detected without very deep observations and spectroscopies of host galaxies can provide redshift informations even when grb afterglows are not bright enough for spectroscopies .
the reporting bias can be eliminated by finding host galaxies of all grbs known in the redshift range and measuring metallicities of all host galaxies , which is difficult to do at higher redshifts ( @xmath11 ) .
3 .   a wealth of spectroscopic studies of general galaxies at low - redshifts , such as the sloan digital sky survey ( sdss ) , provides us with control sample of galaxy properties , to which we can compare long grb host properties to investigate the relation between long grb rate and metallicity .
the advantages and the disadvantages of a redshift selected study of grb host galaxies at low - redshifts can be summarized as in table  [ tab : procon ] , in comparison with the other strategies of grb host studies .    in this study , we present the overall metallicity distribution of the host galaxies of grbs known at @xmath12 that occurred before the end of march 2014 .
for a few of this low - redshift host galaxies , sufficient spectroscopic data to constrain their metallicity were not previously available in the literature , or signal - to - noise ratio ( s / n ) of important metallicity indicating emission lines was low .
we perform spectroscopic metallicity measurements of the host galaxies of three low - redshift grbs , 060614 , 090417b , and 130427a , to obtain better constraints on the metallicity distribution .
we compare the metallicity distribution of this low - redshift sample to the predictions of empirical models of galaxies , and constrain the relation between metallicity and grb occurrence rate .
this paper is organized as follows .
we describe our sample selection , observations , and analysis of the obtained spectra in section  [ sec : sample ] , [ sec : obs ] , and [ sec : linefit ] , respectively . in section
[ sec : prop ] , we discuss the properties of the low - redshift long grb host galaxies including metallicities . in section  [ sec : moedel ] , we describe models of galaxy metallicity distributions , to which we compare the low - redshift sample of host galaxies to constrain the relation between metallicity and grb occurrence rate . the results of the comparison are presented in section  [ sec : result ] . in section  [ sec : discussion ]
, we discuss uncertainties of our galaxy model , possible subpopulations in the low - redshift grbs , and sampling biases we may suffer by collecting long grbs with known redshifts .
we summarize our conclusions in section  [ sec : conclusion ] .    throughout this paper , we assume the fiducial cosmology with @xmath13 , @xmath14 , and @xmath15 70 km s@xmath16 mpc@xmath16 .
the magnitudes are given in the ab system .
[ tab : procon ]    @xmath17 this work .
we collect low - redshift grbs from online databases : swift gamma - ray burst table , gamma - ray burst online index ( grbox ) , j. greiner s grb table , and grb host studies ( ghosts ) .
although there are various classes of gamma - ray burst like events : e.g. , long grbs , short grbs , ultra - long grbs , and xrfs , we focus on host galaxies of long grbs in this study .
the classification of burst events are still controversial ( e.g. , @xcite ) , and important burst properties for the classification such as burst duration @xmath18 , and spectral peak energy @xmath19 are uncertain in some cases . however , for simplicity , we select long grbs simply by reported @xmath18 and @xmath19 as follows .
the prompt emission properties are taken from @xcite , @xcite , @xcite , and @xcite , in addition to the online databases mentioned above .
first , we collect bursts with 2 @xmath20 sec as long grbs , considering those with @xmath21 sec as short grbs and @xmath22 sec as ultra - long grbs
. we also exclude short grbs with extended emissions , grb 050709 , 050724 , 061210 , and 071227 from our sample . from the bursts which meet the duration criteria , we remove bursts with very soft spectra ( xrfs ) which may be a different population from long grbs .
xrfs are often defined according to their hardness ratio which depends on energy bands of x - ray detectors . to take into account burst events discovered by various instruments ,
we define bursts with @xmath23 kev as xrfs in this study . for some _ swift _ detected bursts ,
@xmath19 estimation is not available . in that case , we distinguish grbs and xrfs using photon spectral index in the energy range of _ swift_-bat ( @xmath24 ) .
we consider bursts with @xmath25 which is typical of long grbs with @xmath26 kev @xcite as long grbs , and those with @xmath27 as xrfs .    within the criteria described above , 15 long grbs occurred before the end of march 2014 at spectroscopically measured redshifts @xmath12 ( table  [ tab : grbs ] )
however , we exclude grb 050219 from further discussion in this study , leaving 14 low - redshift long grbs in our sample , because it is unclear whether this burst really has occurred at a low - redshift .
grb 050219 occurred close to ( but @xmath1 offset from ) an early - galaxy for which the redshift @xmath28 is measured @xcite .
given that long grbs result from death of massive stars , the association of grb 050219 with an early - type galaxy is surprising , although it is not strictly prohibited because early - type galaxies can harbour small amount of star formation ( e.g. , @xcite ) .
the significant offset of grb 050219 from its claimed host galaxy also suggests that the galaxy might be aligned with the grb by chance without any physical association with the burst .
it is possible that the remaining sample is also contaminated by galaxies that coincide with grbs by chance .
we discuss this issue in section  [ sec : contamination ] .    among the 14 grbs ,
grb 060614 , 090417b , 111225a , and 130427a were without sufficient spectroscopic observations to significantly constrain their host metallicity .
we perform spectroscopic observations of the host galaxies of grb 060614 , 090417b , and 130427a in this study .
there are also some recently published spectroscopic observations of the host galaxies of grb 060614 and 130427a .
we also discuss those observations in section  [ sec : obs ] .
the host galaxy of grb 111225a missed our spectroscopic observation because its redshift had not been measured as of the end of march 2014 , and no metallicity information is available for this host galaxy .
the redshift of this burst was reported @xmath1 2.3 years after the burst @xcite .
we take into account this host galaxy as a part of the uncertainty in the metallicity distribution of the grb host galaxies .
we discuss the implications for the sampling bias in the redshift selected sample from the late redshift report of this burst in section  [ sec : fainthost ] .
hereafter , grb means long grb unless otherwise stated .
[ tab : grbs ]    long grbs whose redshifts are spectroscopically determined to be @xmath12 .
short grbs , xrfs , and ultra - long grbs are not included .
+ @xmath17 if ` host galaxy ' , the redshift is determined by emission lines of the host galaxy .
if ` afterglow ' , the redshift is determined by absorption lines in the afterglow .
the numbers in the square brackets are the references as listed below .
+ references : 1 .
@xcite , 2 .
@xcite , 3 .
@xcite , 4 .
@xcite , 5 .
@xcite , 6 .
@xcite , 7 .
@xcite , 8 .
@xcite , 9 .
@xcite , 10 .
@xcite , 11 .
@xcite , 12 .
@xcite , 13 .
@xcite , 14 .
@xcite , 15 .
@xcite , 16 .
@xcite , 17 .
@xcite , 18 .
@xcite , 19 .
@xcite , 20 .
@xcite + @xmath29 bursts with confirmed sn associations @xcite + @xmath30
we do not include grb 050219 in the analyses in this paper , because the association of this bust with its host galaxy for which the redshift is measured is highly uncertain .
we perform optical spectroscopy of the host galaxies of grb 060614 , 090417b , and 130427a with the gemini multi - object spectrographs ( gmos , @xcite ) , and use flux ratios of the emission lines : [ o]@xmath313727 , h@xmath32 , [ o]@xmath314959 , [ o]@xmath315007 , h@xmath33 , and [ n]@xmath316584 , as metallicity indicators ( hereafter , [ o ] , [ o ] , and [ n ] mean [ o]@xmath313727 , [ o]@xmath315007 , and [ n]@xmath316584 , respectively , unless otherwise stated ) .
the apparent spatial extent of the three host galaxies are @xmath1 , and major part of their flux is collected in the slit .
the spectra obtained with gmos were calibrated using spectroscopic standard stars observed on different nights from the observations of the host galaxies , and hence the absolute flux scale is not accurate .
emission line fluxes of the grb 130427a host galaxy obtained during a target - of - opportunity program ( too , pi : n. kawai ) using subaru / focas @xcite are presented together .
we also performed imaging observations of the host galaxy of grb 111225a .
we describe the observations of each target in the following subsections .
see table  [ tab : obs ] for the summary of our observations .
the data were reduced in a standard manner using pyraf/ iraf , together with the gemini iraf package and the focasred package .
grb 060614 is known as a grb unassociated with a bright sn at @xmath34 @xcite .
although some spectroscopic observations of the grb 060614 host galaxy were performed soon after the burst @xcite , the metallicity measurements were uncertain because the fluxes of [ n ] and [ o ] lines were not well constrained .
recently , @xcite also reported metallicity measurement of the host galaxy of grb 060614 .
the metallicity they derived is broadly consistent with ours which we discuss in section  [ sec : metal ] .
we performed spectroscopy of the grb 060614 host galaxy with gmos - south using the slit and two different settings of grisms and order - cut filters .
one is the b600 grating ( no - filter ) which covers 35006000   , and the other is the r400 grating + the og515 order - cut filter covering 53009000   .
the signals were @xmath35 pixels binned .
the b600 spectroscopy was performed on 2015 april 25 ( ut ) .
the integration time was 1000 sec @xmath36 , and the slit position angle was set to the parallactic angle to minimize the effect of differential atmospheric refraction .
the r400 spectroscopy was performed on 2015 april 28 , and may 14 with the integration time of 1200 sec @xmath37 and the position angle set to the parallactic angle .
grb 090417b is an optically dark grb at @xmath38 with a very long duration of @xmath39 sec @xcite .
the spectroscopy was performed with gmos - north on 2015 march 25 .
we use the slit , the r831 grating , and the rg610 order - cut filter covering 73009200   . at this redshift ,
the h@xmath33 line coincides with a strong night - sky emission line at 8827   . to accurately subtract sky emission lines
, we utilize the nod and shuffle sky subtraction .
the signals were 2 pixels binned along the spectral direction .
the total integration time was 2640 sec ( 60 sec @xmath40 2 positions @xmath40 11 cycles @xmath40 2 sequences ) .
the slit position angle was set to the parallactic angle .
grb 130427a is a bright burst at @xmath41 associated with a broad - lined type ic sn 2013cq @xcite .
although detection of the host galaxy emission lines over the afterglow and the sn are reported by @xcite and  ( , hereafter k15 ) , the detection of the [ n ] line is marginal in either case .
the spectroscopic follow up observation of grb 130427a with subaru / focas was performed on 2013 may 17 ( i.e. , @xmath1 20 days after the burst ) , under a weather condition with cirrus clouds ( pi : n. kawai ) .
the slit , the 300b grating , and the y47 filter covering 47009000    were used . the slit position angle was set to to put both of the center of the host galaxy and the grb position into the slit .
the signals were 2 pixels binned along the spatial direction .
the integration time was 9600 sec .
unfortunately the atmospheric diffraction corrector ( adc ) could not be used at the time of the observation , however the spectroscopy was performed with elevation angles @xmath42 , and atmospheric refraction would be @xmath43 throughout the wavelength range .
although the s / n of the focas spectrum is higher than those previously reported , the significance of the detection of the [ n ] line was still low @xmath44 ( @xmath45 in k15 )
. furthermore , the strong night - sky emission line at 8827     overlaps with the [ n ] line of the host galaxy at this redshift . hence accurate subtraction of the sky emission is essential for the spectroscopy of this host galaxy .
we performed further spectroscopy of the grb 130427a host galaxy with gmos - north on 2015 march 27 to measure the [ n ] line flux more securely with higher significance of the line detection .
the total integration time was 2640 sec .
the slit position angle was set to the parallactic angle .
we used the nod and shuffle sky subtraction to subtract the sky emission line accurately .
the signals were 2 pixels binned along the spectral direction .
@xcite reported the redshift of grb 111225a , @xmath46 , in april 2014 ( @xmath1 2.3 years after the burst ) by analyzing the archival data with an updated software .
thus grb 111225a was not included in the target list of our spectroscopic observations which was compiled in march 2014 .
we perform imaging observations of the host galaxy of grb 111225a in @xmath47 , and @xmath48band filters with subaru / focas on 2015
september 22 , under a photometric condition .
the integration time is 1800 sec for @xmath49band , and 600 sec for each of @xmath50 , and @xmath48band .
signals were @xmath35 pixels binned .
photometric calibration is performed with respect to @xcite standard stars in the pg0231 + 051 field observed in the same night .
[ tab : obs ]    the observations of the grb host galaxies we present in this study .
+ @xmath17 use of nod & shuffle sky subtraction .
+ @xmath29 parallactic angle + @xmath30 the integration time of 1 nod & shuffle sequence consists of 60 sec @xmath40 2 positions @xmath40 11 cycles .
+ @xmath51 too , poor weather condition
to minimize the effect of stellar absorption features on the line flux measurements , we first subtract stellar spectral energy distribution ( sed ) models of the host galaxies from the observed spectra .
we perform sed fittings to optical / near - infrared ( nir ) broad - band photometries in the literature ( figure  [ fig : sedfit ] ) using the _ sedfit _ software package @xcite which utilizes population synthesis models by @xcite .
the extinction law by @xcite and the initial mass function ( imf ) by @xcite are assumed .
we examine six cases of stellar metallicity : @xmath52 , and five cases of star formation history : simple stellar population ( ssp ) , constant star formation , and exponentially decaying star formation with @xmath53 and 5 gyr .    as mentioned in section  [ sec : obs ] , absolute scale of our spectra is not accurate .
furthermore , it is possible that a fraction light is lost at the slit . in the following discussion ,
we rescale the observed spectra so that the continuum flux agrees with the best fit sed models .
the broad - band photometries , the best fit sed models , and the rescaled spectra are shown in figure  [ fig : sedfit ] .
the parameters of the best fit stellar sed models are listed in table  [ tab : pop ] .
we derive metallicities of the grb host galaxies using relative flux ratios between different emission lines as indicators ( section  [ sec : metal ] ) , and hence the derived metallicities are not significantly affected by the error of absolute flux scale .        [ tab : pop ]    the parameters of the best fitting stellar sed models @xcite that reproduce the broad - band photometries in the literature ( the upper left panel of figure  [ fig : sedfit ] ) .
the extinction law by @xcite and the imf by @xcite are assumed .
there are some archival spectra of the grb 060614 host galaxy as mentioned in section  [ sec : obs:060614 ] , and those by @xcite and @xcite are obtained with gmos - south using the slit and the r400 grating as well as our r400 spectrum .
although the gmos archival data alone are not deep enough to constrain the metallicity of the host galaxy , we stack them with our r400 spectra to maximize the s / n .
the gmos - south detectors were replaced between the archival observations and ours , the pixel scale remained @xmath54  / pix with the r400 grating .
the integration time of the archival observations is 1800 sec @xmath40 4 and 1200 sec @xmath40 4 , respectively .
we note that the detector replacement improved the efficiency at the wavelength of h@xmath33 and [ n ] at @xmath34 by @xmath1 35% .
we subtract the best fit sed models from the rescaled spectra , and measure emission line fluxes by fitting the lines with gaussian plus flat continuum models .
the residual spectra after the model subtraction and the gaussian fits are shown in figure  [ fig : fitha ] and figure  [ fig : fitoxhb ] . although [ o ] line is a doublet , we fit it with a single gaussian as well as the other lines , because the doublet is not resolved in our spectrum .
the obtained emission line fluxes are corrected for the foreground extinction in the milky way ( mw ) using the extinction map by @xcite and the extinction law by @xcite , and shown in table  [ tab : lines ] .
the continuum component in the spectrum obtained by the focas too observation of grb 130427a is contaminated with the afterglow and the sn ( figure  [ fig : too130427a ] ) .
hence we can not calibrate the flux scale comparing the continuum spectrum with the best fit stellar sed model .
instead , we firstly measure h@xmath33 flux without the absorption correction , and rescale the focas spectrum so that the absorption uncorrected h@xmath33 flux agrees with that obtained from the gmos spectrum .
the emission line fluxes of the other low - redshift grb host galaxies in the literature are also shown in table  [ tab : lines ] .
the host galaxies of grb 980425 , 060505 , and 120422a , are spatially resolved in the spectroscopic observations @xcite , and emission line fluxes are measured at multiple positions within the galaxies . in that case , we use the integrated flux over the whole galaxy or the flux measured at the peak of the light profile of the host galaxy , for consistency with the spectroscopic data of the other grb host galaxies which are not spatially resolved .        , but for the [ o ] , h@xmath32 , [ o ] lines of the grb 060614 host galaxy .
the spectrum is stacked with the the archival data @xcite in the wavelength ranges of h@xmath32 and [ o ] .
, width=302 ]
the emission line fluxes of the host galaxies in the low - redshift sample ( see section  [ sec : sample ] ) are listed in table  [ tab : lines ] . with these emission line measurements ,
we derive metallicity of ionized gas , 12+log@xmath55(o / h ) , which would be close to the metallicity of young stars in the galaxies .
we utilize both n2 index (=  @xmath56  [ n]@xmath316584/h@xmath33 ) and @xmath57 index [ =  ( [ o]@xmath313727+[o]@xmath314959+[o]@xmath315007)/h@xmath32 ] as indicators of metallicity .
various calibrations are proposed for the relation between the indiators and metallicity , and different calibrations are often inconsistent with each other @xcite . in this study , we utilize the method of  ( , kk04 ) , in which both n2 and @xmath57 indices are calibrated using the same photoionization model , to obtain consistent results from the two indicators .
the solar - metallicity @xmath59 corresponds to 12+log@xmath55(o / h ) = 8.69 @xcite .    in the kk04 calibration , the n2 and @xmath57 indices
are related to metallicity taking the effect of ionization parameter ( @xmath60 ) into account .
the @xmath60 parameter can be estimated from @xmath61 index [ =  ( [ o]@xmath314959+[o]@xmath315007)/[o]@xmath313727 ] .
the metallicity calibration by kk04 suffers from the model uncertainties of @xmath1 0.1 dex , and hence we consider error of a metallicity measurement to be 0.1 dex when the error propagation from the emission line fluxes is smaller than this value or flux error informations are not available .
the fluxes in table  [ tab : lines ] are corrected for the foreground extinction in the mw , but not for the extinction in the host galaxies themselves .
to derive metallicity , we correct the line fluxes for the host extinctions , assuming the intrinsic value of the balmer decrement h@xmath33/h@xmath32 = 2.86 and the extinction law by @xcite .
the positions of the host galaxies on the baldwin - phillips - terlevich ( bpt ) diagram @xcite with the extinction correction are shown in figure  [ fig : bpt ] . in the bpt diagram ,
the grb host galaxies show emission line ratios that are consistent to result purely from star forming activities .
however , @xcite showed that the spectrum of the host galaxy of grb 031203 is possibly affected by an active galactic nucleus ( agn ) investigating more numerous emission line fluxes than examined in the bpt diagram .
hence the derived properties of the grb 031203 host galaxy might be systematically affected by an agn .
we describe the metallicity determination of each host galaxy in the appendix .
the derived metallicities and the extinctions are listed in table  [ tab : metal ]  &  [ tab : msfr ] .
when both of the n2 and @xmath57 indices are valid for a galaxy , we use the @xmath57 index , although the two indices agree with each other except for the host galaxies of grb 031203 and 011121 .
the emission line ratio of the grb 031203 host galaxy might be affected by an agn as mentioned above , and the grb 011121 host galaxy is a remarkable outlier from the mass - metallicity ( @xmath62-@xmath63 ) relation of galaxies ( see section  [ sec : result ] ) .
the metallicity measurement of these objects should be considered with care .
[ tab : lines ]    fluxes are in units of @xmath64 erg s@xmath16@xmath65 , and corrected for the foreground extinction in the mw .
when the estimated flux is consistent with zero within the 2@xmath66 error , we show the 95% upper - limits instead .
+ @xmath17 the flux errors are @xmath1 10% .
+ @xmath29 flux error informations are not available .
+ @xmath30 the h@xmath33 and [ n ] fluxes are measured from the gmos spectrum , and the other line fluxes are measured from the focas spectrum .
indices do not agree well with each other for these host galaxies .
the vertical dotted line represents the [ n]/h@xmath33 ratio of the host galaxy of grb 090417b for which the [ o]/h@xmath32 ratio is not known . ,
width=302 ]    [ tab : metal ]    @xmath67 kk04 calibration method is utilized .
+ @xmath68 the estimation errors are @xmath69 0.1 unless otherwise stated .
+ @xmath70 the best estimates as described in the appendix .
+ @xmath71 possibly contaminated by an agn .
+ @xmath72 derived from the emission line fluxes by k15 .
+ @xmath73 assumed value ( see the appendix ) +      stellar mass ( @xmath62 ) and sfr of the low - redshift grb host galaxies are listed in table  [ tab : msfr ] .
the stellar masses of the host galaxies of grb 060614 , 090417b , and 130427a are derived by the sed fitting assuming the chabrier imf as described in section  [ sec : linefit ] , and are broadly consistent with previously reported results @xcite .
the grb 111225a host galaxy is detected only in 4 bands ( @xmath74 , and @xmath75 , see table  [ tab:111225a ] ) , and it is difficult to derive the properties of this host galaxy from its sed .
however , we can naively estimate its stellar mass as @xmath76 assuming typical mass - to - luminosity ratio of galaxies at similar redshifts @xcite .
we collect @xmath62 measurements of the other host galaxies from the literature . when @xmath62 in the literature is derived assuming an imf other than the chabrier imf , we rescale the mass assuming : @xmath77 @xcite , for consistency .
we derive sfr of the host galaxies from their extinction corrected h@xmath33 fluxes assuming the global kennicutt - schmidt law @xcite and the same scaling factor between the different imfs as mentioned above .
we note that major part of the errors of sfr is propageted from the error of @xmath78 in most of the cases .
the host galaxy of grb 060505 is spatially extended over a few @xmath40 10 arcsec , and only a part of the emission line fluxes is collected in the slit , and hence we adopt a sfr measured from ultraviolet photometries without dust correction instead @xcite .
although the dust - uncorrected sfr is much larger than that suggested from the h@xmath33 flux with the large slit - loss , the actual sfr with dust correction might be further high .
the emission lines of the grb 120714b host galaxy were detected over the optical afterglow ( @xcite ; k15 ) , and the fraction of the host galaxy light in the slit is not known .
hence we consider the obtained sfr for the host galaxies of grb 060505 and 120714b as lower - limits .
[ tab:111225a ]    magnitudes are in the ab system , and not corrected for the foreground extinction in the mw [ @xmath79 , @xcite ] .
+ @xmath17 2@xmath66-limit .
[ tab : msfr ]    @xmath67 estimated from the balmer decrement + @xmath68 the chabrier imf is assumed .
the numbers in the square brackets are references as listed below .
+ @xmath70 sfr is derived from the h@xmath33 flux in table  [ tab : lines ] and @xmath80 in the second column unless otherwise stated .
+ @xmath71 possibly contaminated by an agn .
+ @xmath72 the extinction correction is based on the result of the sed fitting : @xmath81 .
+ references : 1 .
@xcite , 2 .
@xcite , 3 .
@xcite , 4 .
@xcite , 5 . this work , 6 .
@xcite , 7 .
@xcite , 8 .
to investigate the relation between metallicity and grb occurrence , we compare the metallicity distribution of the low - redshift sample presented in section  [ sec : prop ] with the metallicity distribution of star forming galaxies at similar redshifts .
following @xcite and @xcite , we compute the metallicity distribution of star forming galaxies using the empirical formulation of stellar mass function [ @xmath82 , @xmath62-sfr relation ( so called `` galaxy main - sequence '' ) , and @xmath62-@xmath63 relation of galaxies . in this study , we assume @xmath83 by  ( , blue population ) , the @xmath62-sfr relation by  ( , see their equation  11 , scatter @xmath84 dex ) , and the @xmath62-@xmath63 relation computed with the kk04 calibration by @xcite .
we only consider the mass range log@xmath85 8.0 , below which the properties of galaxies are not well constrained .
this mass limit is also comparable to the lowest - mass grb host galaxies in the low - redshift sample . @xmath62 and sfr in the formulations
are rescaled to be consistent with the chabrier imf .
the scatter of the @xmath62-@xmath63 relation is largely dependent on @xmath62 in the sense that the scatter is larger at smaller @xmath62 .
here we assume a simple formulation : @xmath86 which is broadly consistent with the scatter of the @xmath62-@xmath63 relation discussed in @xcite .
it should be noted that @xcite derived the @xmath62-@xmath63 relation using galaxy spectra obtained by the sdss whose fiber spectrograph covers only the central of each galaxy .
hence it is possible that the systematic loss of light from the outskirts of galaxies affects the derived metallicities , because galaxies often have lower metallicity in their outskirts than at their centers ( so called metallicity gradient , e.g. , @xcite ) .
some of the low - redshift grb host galaxies might also suffer from similar problems .
however , @xcite examined correlation between measured metallicities and fraction of light covered by the spectrograph using a sample of galaxies with similar @xmath62 and sfr in the sdss , and found that the dependence of derived metallicities on the fiber coviring fraction is small ( @xmath87 dex , see figure  5 of @xcite ) .
thus we consider that the systematic loss of light from the outskirts of galaxies does not significantly affect our metallicity measurements    although it is suggested that there is a correlation between sfr and metallicity of galaxies with similar stellar masses @xcite , the observed correlation between sfr and @xmath63 correlation is not quantitatively understood and possibly affected by sample selections and metallicity calibration methods ( e.g. @xcite ) .
hence we do not consider the correlation in our baseline model , and discuss the effects of the sfr-@xmath63 correlation on our results in section  [ sec : msfrz ] .    in this study ,
we consider normalized probability distribution function ( pdf ) of grb host galaxy metallicities when we compare the model predictions to the observations .
we note that the absolute scale of the stellar mass function ( @xmath88 ) , and also the absolute scale and the scatter of the @xmath62-sfr relation do not affect the predicted metallicity pdf . @xcite
have shown that the choice of the formulations of the stellar mass function and the @xmath62-sfr relation does not significantly affect the predicted metallicity pdf , as far as the stellar mass function and the @xmath62-sfr relation that represent field star forming galaxies are selected .
observations of some nearby galaxies ( mw , the magellanic clouds , and m31 ) show that inter - stellar medium ( ism ) in a galaxy is not chemically homogeneous .
hence the metallicity of a grb progenitor star might be different from the metallicity of the host galaxy which we can measure by follow up spectroscopies .    to examine metallicity variation of star forming regions within a galaxy , we consider observed metallicity distributions of h regions in nearby galaxies .
although recent integral field unit ( ifu ) spectroscopies provides us with spatially resolved map of metallicity in some star forming galaxies , the spatial resolution of such observations are still limited to @xmath89 100 pc in most of the cases ( e.g. , @xcite ) .
hence the spectroscopy of individual h regions ( except for giant ones ) can be performed only in a few nearby galaxies .    in the upper panel of figure  [ fig : zint ]
, we show the metallicity distributions of h regions in the mw ( solar - neighborhood , @xcite ) , the large / small magellanic clouds ( lmc and smc , @xcite ) , and the andromeda galaxy ( m31 , @xcite ) .
it is possible that the metallicity estimates of the faint h regions in the m31 sample by @xcite have systematic errors @xcite , and hence we consider only h regions with h@xmath33 luminosity @xmath90 [ erg s@xmath16 ] in the m31 sample .
the h@xmath33 luminosities of the h regions are taken from @xcite .
see @xcite for the detail of the catalog matching between the samples of @xcite and @xcite .
the metallicities of the h regions can be represented by log - normal distributions but with different median value and dispersion in different galaxies .
the median and the scatter of the h region metallicity distributions in the nearby galaxies [ 12+log@xmath55(o / h)@xmath92 and @xmath93 are shown in the lower panel of figure  [ fig : zint ] .
although the sample number is small , it is seen that the galaxies with the higher median metallicity have larger scatter of the distribution . as a baseline model , we assume a broken linear relation between 12+log@xmath55(o / h)@xmath92 and @xmath94 as shown in the lower panel of figure  [ fig : zint ] : @xmath95 & ( -3.7 \leq { \rm log_{10}(\frac{o}{h})_{gal } } <
-3.4 ) \\      0.3 & -3.4 ( \leq { \rm log_{10}(\frac{o}{h})_{gal } } ) \\
\end{array}\right . , \end{aligned}\ ] ] considering 12+log@xmath55(o / h)@xmath92 equals to the representative metallicity of a galaxy which we obtain by spectroscopically observing the galaxy without spatially resolving it .
we also discuss the results with different @xmath94 formulations in section  [ sec : sigma_zint ] .    in figure
[ fig : zint ] , we plot the metallicity distribution of h regions in the grb 980425 host galaxy together with those in the local galaxies .
the grb 980425 host galaxy is the nearest grb host galaxy known and the only grb host galaxy for which metallicity is measured with a spatial resolution @xmath96 kpc , although the resolution @xmath97 pc is much larger than that achieved for the local galaxies ( e.g. , @xmath98 pc for m31 , @xcite ) .
the internal metallicity variation observed within the grb 980425 host galaxy is naively consistent with those seen in the magellanic clouds .
metallicity distribution of young stars which are formed in a galaxy might be significantly different from the metallicity distribution of h regions in the galaxy , if metallicity and sfr of h regions are correlated with each other .
@xcite examined the correlation between h@xmath33 luminosity and metallicity of h regions in m31 using the narrow band photometric data by @xcite and the spectroscopic data by @xcite .
although it is suggested that metallicity measurements of faint h regions with h@xmath33 luminosity @xmath99 [ erg s@xmath16 ] have luminosity dependent systematic error depending on the metallicity calibrator , no correlation was seen for h@xmath33 luminosity @xmath90 [ erg s@xmath16 ] .
hence , we consider that the correlation between metallicity and sfr of h regions would not be strong even if it exists , although the existence of the correlation is not robustly ruled out .
the metallicities of the nearby h regions discussed in this section are measured using different metallicity calibration methods than the kk04 calibration which we utilize in this study , and it is known that different metallicity calibrations are often inconsistent with each other @xcite .
we consider this possible inconsistency as a part of the uncertainty of @xmath94 which we discuss in section  [ sec : sigma_zint ] .
it is broadly agreed that the metallicity decreases as the galactocentric radius increases in many galaxies ( e.g. , @xcite ) .
however , we note that a scatter around the gradient at each radius is also significant .
the scatter of metallicity at each galactocentric radius can be comparable to ( or evan larger than ) the radial variation ( e.g. , @xcite ) , and lowest - metallicity h regions in a galaxy may reside close to the center of the galaxy .
it should be noted that the scatter of the metallicities in each galaxy may in part result from uncertainties of the metallicity measurements .
however , similar variation of metallicity is also found in the mw using stars with ages of @xmath100 gyr as a tracer of metallicity @xcite .
furthermore , in the m31 h region sample of @xcite , the weak auroral line [ o]@xmath314363 is detected for four h regions .
with the electron temperature ( @xmath101 ) method of metallicity measurement using the auroral line ( e.g. , @xcite ) , which is not affected by uncertainties of the photoionization models , they showed that the four h regions have 12+log@xmath55(o / h ) @xmath1 8.38.4 .
it is generally difficult to detect an auroral line of a high - metallicity h region , and it is not surprising that the four h regions with the auroral line detections are biased towards low - metallicities .
however , the @xmath101 method measurements provides a confirmation that some h regions in m31 have low - metallicity in reality .
[ erg s@xmath16 ] ) , and the host galaxy of grb 980425 which resides at @xmath102 .
the best fit log - normal distributions are plotted together .
_ lower panel _ : the median [ 12+log@xmath55(o / h)@xmath92 ] and the scatter ( @xmath94 ) of the h region metallicity distributions in the nearby galaxies .
the pentagonal symbol represents the grb 980425 host galaxy .
the correlations between 12+log@xmath55(o / h)@xmath92 and @xmath94 which we examine in this study : a broken linear model ( solid line ) , a single linear model ( dashed line ) , and constant models ( dotted lines ) are plotted together .
, width=302 ]      we parameterize the relation between grb production efficiency ( @xmath104/sfr ) and metallicity of a population of young stars . as discussed in section  [ sec : zint ] , the metallicity of the stellar population may be different from that of the host galaxy .
here we consider that @xmath105 is suppressed above a threshold metallicity ( o / h)@xmath106 by a factor of @xmath107 which is not necessarily zero ( step function model ) : @xmath108 we note that absolute scale of the efficiency ( i.e. , @xmath109 ) is marginalized when we consider normalized metallicity pdf of galaxies .
we also consider the cases with the power - law @xmath105 model : @xmath110^{\alpha } & { \rm o / h > ( o / h)_{cut } }     \end{array}\right .
, \label{eq : eff_pow}\end{aligned}\ ] ] and with the exponential @xmath105 model : @xmath111^{\beta } ) .
\label{eq : eff_exp}\ ] ] the schematic picture of the @xmath105 models is shown in figure  [ fig : eff ]    grb rate of a galaxy with a given set of sfr and 12+log@xmath55(o / h)@xmath92 ( median metallicity of star forming regions in the galaxy ) can be obtained by performing the following integration : @xmath112 where @xmath113 is the metallicity distribution of star forming regions in the galaxy which is represented by a log - normal distribution with a median value that corresponds to 12+log@xmath55(o / h)@xmath92 and a log scale scatter @xmath94 ( see section  [ sec : zint ] ) .    ) , power - law ( dashed , eq .
[ eq : eff_pow ] ) , and exponential ( dotted , eq .
[ eq : eff_exp ] ) models of @xmath105 as functions of metallicity .
, width=302 ]
assuming @xmath83 , the @xmath62-sfr relation , and the @xmath62-@xmath63 relation , we compute number density of galaxies in the parameter space of @xmath62 , sfr , @xmath63 : @xmath114 [ mpc@xmath115dex@xmath115 ] .
the expected metallicity distribution of grb host galaxies can be derived integrating @xmath114 over @xmath62 and sfr weighted with @xmath117 ( equation  [ eq : rgrb ] ) .    each @xmath105 model discussed in section  [ sec : eff ] has two free parameters : [ ( o / h)@xmath106 , @xmath107 ] , [ ( o / h)@xmath106 , @xmath33 ] , or [ ( o / h)@xmath106 , @xmath32 ] . the parameter ( o / h)@xmath106 determines the cutoff metallicity above which the grb production efficiency is suppressed , and @xmath107 , @xmath33 , or @xmath32 determines the sharpness of the cutoff .
we search for parameter sets that reproduces the observed metallicity distribution of the low - redshift sample , through the parameter ranges @xmath120 , @xmath121 , @xmath122 , and @xmath123 with intervals of @xmath124 , and @xmath125 .
the goodness of fit is determined using the kolmogorov - smirnov ( ks ) test against the `` median '' distribution shown in the upper left panel of figure  [ fig : result ] .
although the metallicity of the grb 111225a host galaxy is not measured and those of the host galaxies of grb 061021 and 090417b are not precisely constrained ( see section  [ sec : metal ] ) , the uncertainties of these metallicities do not significantly affect the overall metallicity distribution of grb host galaxies .
the resulting parameters from the fitting are not significantly changed within the range of this uncertainty .
the bestfit model distributions are plotted in the upper left panel of figure  [ fig : result ] together with the metallicity distribution of the low - redshift sample .
the best fitting metallicity distributions with the three different @xmath105 models are very similar to each other .
the sfr weighted metallicity distribution predicted from the same @xmath83 , the @xmath62-sfr relation , and the @xmath62-@xmath63 relation , which the grb host galaxies would follow if @xmath105 is not dependent on metallicity , is also shown together .
the acceptable range of the parameters are shown as contour maps of the k - s test probability @xmath126 in figure  [ fig : result ] .
the parameter ranges that the @xmath127 are :    step function : : :    12+log@xmath55(o / h)@xmath128 ,    @xmath129 power - law : : :    12+log@xmath55(o / h)@xmath130 ,    @xmath131 exponential : : :    12+log@xmath55(o / h)@xmath132 ,    @xmath133    the degeneracy between the two parameters is seen in the cases of the power - law and exponential models .
the best fit 12+log@xmath55(o / h)@xmath106 is @xmath1 8.3 , and the sharper cutoff reproduces the observations better in any of the models .
this means that majority of low - redshift grb host galaxies have higher metallicity than 12+log@xmath55(o / h)@xmath106 above which grbs can not be produced .
in other words , most low - redshift grbs take place in star forming regions whose local metallicity is much lower than the representative value of their host galaxies .
the model distribution with the step function model of @xmath105 does not significantly depend on @xmath107 when log@xmath134 , and the distribution with @xmath135 is quite similar to the best fit distribution shown in the upper left panel of figure  [ fig : result ] .
we also note that , in the limit of @xmath136 ( @xmath137 ) , the power - law ( exponential ) model of @xmath105 is identical to the step function model with @xmath135 . in the following ,
we only consider the step function model of @xmath105 with @xmath135 to examine the predicted properties of low - redshift grb host galaxies with our models .    in figure
[ fig : effgal ] , we compare the step function model of @xmath105 [ 12+log@xmath55(o / h)@xmath106 = 8.3 , @xmath135 ] with that convolved over galaxy scale ( @xmath139 ) as defined in equation  [ eq : rgrb ] .
@xmath139 steeply declines with increasing metallicity at around 12+log@xmath55(o / h)@xmath106 .
however , galaxies with 12+log@xmath55(o / h)@xmath140 12+log@xmath55(o / h)@xmath106
have @xmath141 due to the internal variation of metallicity , and their contribution to the cosmic grb rate density is significant because they play a dominant role in the cosmic star formation at low - redshifts .
the steep decline of @xmath139 at 12+log@xmath55(o / h)@xmath142 is also reported by @xcite based on a grb host galaxy sample collected from a wider range of redshifts where metallicity measurements of grb host galaxies are incomplete .
we note that @xmath143 likely evolve with redshift because it depends on the properties of galaxies ( e.g. , @xmath94 ) , unlike @xmath144 which would be determined by stellar physics .    in some cases
, grb positions are spatially resolved from the center ( or brightest part ) of the host galaxies ( e.g. , @xcite ) .
@xcite has shown that the metallicities at the grb positions are systematically lower than that of their host galaxies by @xmath145 dex , but not necessarily as low as suggested by our results .
@xcite pointed out that the actual metallicity of a grb explosion site can not be obtained with a spatial resolution @xmath146 kpc which is typical of the observations of grb positions at @xmath147 , because the length scale of the ism metallicity variation is @xmath96 kpc in nearby galaxies ( e.g. , @xcite ) .
our results suggest that the actual metallicity of the grb explosion sites ( or the progenitor star itself ) which is buried within the spatial resolution is significantly lower than the kpc scale metallicity currently observed , by up to @xmath1 0.5 dex .        ) of a stellar population with a single metallicity [ solid line , the step function model with 12+log@xmath55(o / h)@xmath106 = 8.3 and @xmath135 ] , and the efficiency convolved over galaxy scale taking the internal variation of metallicity within each galaxy into account ( dashed line ,
[ eq : rgrb ] ) . , width=302 ]    the @xmath117 weighted distribution of galaxies on the @xmath62-@xmath63 and @xmath62-sfr parameter planes , predicted using the step function model of @xmath105 with 12+log@xmath55(o / h)@xmath106 = 8.3 and @xmath135 , is plotted as contours in figure  [ fig:2ddist ] .
the sfr weighted model distribution , the assumed @xmath62-@xmath63 relation @xcite , the @xmath62-sfr relation @xcite , and the observed properties of the low - redshift sample ( table  [ tab : metal ]  &  [ tab : msfr ] ) are plotted together .
the observed grb host galaxies construct a sequence which is offset from the @xmath62-@xmath63 relation by 0.10.2 dex towards lower - metallicities , as previously noted by @xcite with a smaller incomplete sample , except for one remarkable outlier : the grb 011121 host galaxy , whose n2 and @xmath57 metallicities do not agree well with each other .
the @xmath117 weighted model distribution also shows similar offset on the @xmath62-@xmath63 plane due to the higher @xmath139 in galaxies with lower 12+log@xmath55(o / h)@xmath92 .
although the offset of the observed grb host galaxies may be slightly larger than that of the model distribution , this possible contradiction can be resolved by taking the sfr-@xmath63 correlation into account as we discuss in section  [ sec : msfrz ] .    in the @xmath62-sfr plane ,
the predicted galaxy distributions with the sfr and @xmath117 weightings are offset from the @xmath62-sfr relation by @xmath1 0.5 dex ( @xmath148 , as generally expected when the scatter of a relation is represented by a log - normal distribution ) .
the observed grb host galaxies are broadly consistent with the predicted offset , although the sfr scatter of the host galaxies with log@xmath149 is large .
the host galaxy of grb 031203 , which might have an agn , agrees well with the sequence of other grb host galaxies in the @xmath62-@xmath63 plane , while it is largely offset towards higher - sfr in the @xmath62-sfr plane .
the outliner in the @xmath62-@xmath63 plane , the grb 011121 host galaxy , does not show any peculiarity in the @xmath62-sfr plane .
the host galaxy with the lowest sfr is that of grb 060614 which is not associated with a sn .
the grb 060614 host galaxy does not show any peculiarity in the @xmath62-@xmath63 plane .     vs. @xmath63 ( _ upper panel _ ) , and @xmath62 vs. sfr ( _ lower panel _ ) .
the solid and dotted contours represent the predicted distributions weighted with the @xmath117 and sfr , respectively .
the contours are drawn at a logarithmic interval of 0.25 dex from the peaks of the distributions .
the vertical dashed lines indicate naively estimated @xmath62 of the host galaxy of grb 111225a for which 12+log@xmath55(o / h)@xmath92 and sfr is not measured .
the horizontal dashed lines represent 12+log@xmath55(o / h)@xmath92 and the lower - limit sfr of the host galaxy of grb 120714b whose @xmath62 is not known .
the dashed lines with gray shaded regions are the @xmath62-@xmath63 relation @xcite and the @xmath62-sfr relation @xcite .
the shaded regions represent 1@xmath66-scatter of the relations .
the following host galaxies are marked with specific symbols : grb 011121 ( square , @xmath62-@xmath63 outlier ) , 031203 ( solid circle , possible agn , @xmath62-sfr outlier ) , 060614 ( dashed circle , @xmath62-sfr outlier ) .
the vertical dotted line indicates the lower - limit of the @xmath62 range which we consider when we predict the metallicity distribution of grb host galaxies .
, width=302 ]
the model predictions of the metallicity distribution depend not only on the @xmath105 model , but also on the underlying assumptions of the properties of galaxies . in this section ,
we examine how our results are affected by some galaxy properties which is not well understood , namely the sfr-@xmath63 correlation and the variation of metallicity within a galaxy .
to investigate the effect of the sfr-@xmath63 correlation , we consider the @xmath62-sfr-@xmath63 relation by @xcite , which is an extension of the relation by @xcite towards lower-@xmath62 .
however , the metallicity calibration method used in @xcite is different from that used in this study ( kk04 ) . to examine the effect of the sfr-@xmath63 correlation in a consistent metallicity scale
, we compute the @xmath150 gradient of the @xmath62-sfr-@xmath63 relation around the galaxy main - sequence at each @xmath62 . and assume the gradient without changing the @xmath62-@xmath63 relation @xcite , which has been used in the previous sections .
it should be noted that the scatter of the @xmath62-@xmath63 relation at each @xmath62 may partly result from the sfr-@xmath63 correlation .
@xcite showed that the metallicity scatter around the @xmath62-sfr-@xmath63 relation is smaller by 50% than that of the @xmath62-@xmath63 relation .
thus we assume the scatter of the @xmath62-sfr-@xmath63 relation is 50% of that of the @xmath62-@xmath63 relation defined in equation  [ eq : sigmamz ] .
the resulting @xmath151 contour map and the bestfit metallicity distribution using the step function @xmath105 model are shown in figure  [ fig : gmodel_pks ] ( top panel ) and figure  [ fig : gmodel_dist ] , respectively .
both the best fit parameters and the goodness of the fit are not significantly different from the case without the sfr-@xmath63 correlation .
this situation is similar with the power - law and exponential models of @xmath105 .
the small effect of the sfr-@xmath63 correlation on the metallicity distribution is not surprising .
the sfr-@xmath63 gradient around the galaxy main - sequence is @xmath152 in the low-@xmath62 range ( log@xmath153 ) and shallower in higher-@xmath62 ranges .
given that the grb host galaxies have systematically higher - sfr than the main - sequence by @xmath1 0.5 dex ( section  [ sec : result ] ) , the expected effect of the sfr-@xmath63 correlation is @xmath154 , which is smaller than the uncertainties of the metallicity measurements .
the sfr and @xmath117 weighted model distributions of galaxies on the @xmath62-@xmath63 and @xmath62-sfr planes , predicted assuming the sfr-@xmath63 gradient are shown in figure  [ fig:2ddist_fmr ] . the @xmath105 model and the parameters are the same as in figure  [ fig:2ddist ] .
the sfr weighted galaxy distribution ( i.e. , no metallicity effect ) on the @xmath62-@xmath63 plane is also shifted towards lower - metallicities due to the sfr-@xmath63 correlation , but the shift is smaller than that of the @xmath117 weighted distribution .
furthermore , the grb host galaxies have systematically lower-@xmath62 than the typical @xmath62 of the sfr weighted population ( log@xmath155 ) , which can not be explained solely by the effect of sfr-@xmath63 correlation ( e.g. , @xcite ) .
although the sfr-@xmath63 correlation does not significantly affect the overall metallicity distribution of the grb host galaxies as mentioned above , the predicted number of grb host galaxies above the @xmath62-@xmath63 relation is reduced ( figure  [ fig:2ddist_fmr ] , top panel ) .
this improves the consistency between the predicted and observed distributions of grb host galaxies on the @xmath62-@xmath63 plane compared with the case without the sfr-@xmath63 correlation ( figure  [ fig:2ddist ] ) .    , but with different assumptions of the properties of galaxies .
the step function model of @xmath105 is utilized . _ ( a ) _ : the correlation between sfr and log@xmath55(o / h ) is incorporated . _
( b ) _ : @xmath94 is linearly correlated with 12+log@xmath55(o / h)@xmath92 . _
( c ) and ( d ) _ : @xmath94 is a constant value ( 0.1 and 0.3 dex , respectively ) .
, width=302 ]    , but with different assumptions of the properties of galaxies . the step function model of @xmath105 is utilized .
the plotted models are the same as those investigated in figure  [ fig : gmodel_pks ] .
the black histogram is the same as that in the upper left panel of figure  [ fig : result ] , while the gray histogram represents the metallicity distribution without the host galaxies of grb 011121 , 031203 , 060505 , and 060614 whose burst classifications or metallicity measurements are uncertain .
, width=302 ]    , but with the sfr-@xmath63 correlation .
the anomalies in the contour lines at log@xmath156 result from a discontinuity of the @xmath62-sfr-@xmath63 relation by @xcite .
, width=302 ]      as mentioned in section  [ sec : zint ] , the detailed investigation of the internal ism metallicity variation is performed only for a small number of nearby galaxies .
therefore , @xmath94 of general star forming galaxies and its correlation with other galaxy properties are highly uncertain . given that the nearby high - metallicity galaxies , namely the mw and m31 , have higher @xmath62 than the grb host galaxies in our sample , it is possible that ism in the grb host galaxies have different properties to those in the nearby galaxies .    to investigate how our results depend on the assumed @xmath94 , we perform the parameter fittings assuming other @xmath94 formulations than the baseline model . here
we consider a single linear relation between 12+log@xmath55(o / h)@xmath92 and @xmath94 , and also constant @xmath94 of 0.1 , 0.2 , and 0.3 dex ( figure  [ fig : zint ] , bottom panel ) .
however , it should be noted that the constant @xmath94 is disfavored by observations .
the @xmath151 contour maps and the predicted metallicity distributions are shown in figure  [ fig : gmodel_pks ] and figure  [ fig : gmodel_dist ] for the cases with the single linear relation , and the constant @xmath94 of 0.1 and 0.3 dex .
the single linear relation between 12+log@xmath55(o / h)@xmath92 and @xmath94 overpredicts the contribution of high - metallicity galaxies to the grb production even with the best fit parameters , due to the large @xmath94 of galaxies with 12+log@xmath55(o / h)@xmath157 .
however , the statistical significance of the overprediction is low ( @xmath158 ) .
the best fit parameters are similar to the case with the baseline model .
when @xmath159 dex independently of 12+log@xmath55(o / h)@xmath92 , the preferred cutoff metallicity is 12+log@xmath55(o / h)@xmath160 .
this result is consistent with the previous studies without the internal metallicity variation being explicitly treated in which the cutoff of the grb efficiency at 0.51@xmath59 is suggested ( e.g. , @xcite ) . even in this case ,
high @xmath105 in high - metallicity stellar population is disfavored ( @xmath161 ) .    when @xmath162 dex , the cutoff metallicity is 12+log@xmath55(o / h)@xmath163 , which is close to the results with the baseline model .
the results with the the constant @xmath164 0.2 dex are similar with those with 0.1 , and 0.3 dex , but with 12+log@xmath55(o / h)@xmath165 .      some of the low - redshift grbs are less energetic than typical grbs that are observed up to higher - redshifts ( hereafter classical grbs ) , and possibly are different kind of phenomena .
these bursts , so called low - luminosity grbs , have isotropic equivalent gamma - ray energy @xmath166 [ erg ] ( e.g. , @xcite ) .
@xcite have shown that @xmath167 of a grb might be dependent on the metallicity of its host galaxy based on observations of the host galaxies of 5 grbs / xrfs at redshifts @xmath4 , while later studies @xcite have shown that there is no significant correlation between them with a larger sample from wider range of redshifts ( up to @xmath168 ) .
because low - luminosity grbs can be detected only at low - redshifts @xmath169 , the comparison of the host galaxies and low - luminosity grbs and with those of classical grbs at @xmath170 might be affected by evolution of galaxy populations over the redshits .
thus it is interesting to revisit the correlation between @xmath167 and host metallicity with our sample limited to low - redshifts but much larger than that of @xcite .
@xmath167 of the low - redshift grbs and 12+log@xmath55(o / h)@xmath92 are shown in the left panel of figure  [ fig : eiso_z ] .
although the grbs in the pre-_swift _ era ( until 2004 ) shows a trend that grbs with larger @xmath167 occurr in lower - metallicity galaxies as reported in @xcite , we do not find significant correlation between @xmath167 and 12+log@xmath55(o / h)@xmath92 when grbs in the _ swift _ era are included ( correlation coefficient is -0.36 ) .
grb 060505 and 060614 are widely known as grbs which are not associated with sne @xcite , and it is possible that they are actually short grbs that do not result from core - collapse of massive stars despite their long duration @xmath171 sec ( e.g. , @xcite ) .
these bursts are marked with circles in figure  [ fig : eiso_z ] .
the host metallicities and @xmath167 of grb 060505 and 060614 are typical of low - redshift grbs in our sample .
the metallicity distributions of the grb host galaxies excluding grb hosts with possible systematic error in their metallicity measurements or with uncertain burst classifications ( grb 011121 , 031203 , 060505 , and 060614 ) is shown in figure  [ fig : gmodel_dist ] .
the metallicity distribution without these peculiar objects does not significantly differ from the distribution of the whole sample .          in this study
, we have presented the overall metallicity distribution of the host galaxies of all grbs known at @xmath12 . however , the fraction of grbs with spectroscopic redshifts are @xmath1 30% even in recent observations ( e.g. , @xcite ) , and hence we may suffer from the redshift determination bias by collecting grbs with known redshifts .      one possible source of the redshift determination bias is `` dark '' grbs ( e.g. , @xcite ) .
a `` dark '' grb is a grb with a fainter optical afterglow than expected from its x - ray spectrum .
roughly 1/3 of grbs are `` dark '' , and most of them ( @xmath172% ) result from dust extinction , while some may result from absorption by neutral hydrogen in inter - galactic medium at high - redshifts @xcite .
redshift determination of grbs largely relies on optical spectroscopy of afterglows , and it is difficult to determine redshifts of `` dark '' grbs . thus a sample of grb host galaxy with known redshifts may be biased against dusty galaxies which may tend to have higher - metallicity .
however , at low - redshifts , detection of a grb host galaxy is not quite challenging in most of the cases . assuming the stellar - mass to luminosity ratio derived by @xcite , a galaxy with @xmath173 , which is typical of the low - redshift grb host galaxies ,
would be observed with @xmath174band magnitude @xmath175 at @xmath176 .
once the host galaxy of a grb is detected , x - ray afterglow localization by _
swift_/xrt with typical precision of a few arcseconds , which is available for @xmath172% of the _ swift_/bat - detected grbs , can associate the grb to the host galaxy with a certain level of confidence @xcite . and
the redshift can be determined using emission lines of the host galaxy .    in reality , many redshifts of the low - redshift grbs are determined by the emission lines of the host galaxies ( table  [ tab : grbs ] ) , including that of a `` dark '' grb 090417b ( marked with a square in figure  [ fig : eiso_z ] ) , while redshift determination by host emission lines is rarely available at higher redshifts .
the fraction of the grbs whose redshifts are determined only by the host emission lines in the low - redshift sample ( 8/14 ) is naively comparable to the fraction of grbs without known redshifts in the _ swift_/bat - detected grbs .
thus the bias against grbs in dusty environment is not necessarily effective at low - redshifts .
it is possible that grbs whose redshifts are determined with the emission lines of their host galaxies are biased against grbs that occurred in galaxies with weak emission line fluxes .
emission line luminosities of a galaxy is largely dependent on sfr in the galaxy , which correlates with @xmath62 and 12+log@xmath55(o / h)@xmath92 , in the sense that galaxies with higher @xmath62 and 12+log@xmath55(o / h)@xmath92 have more luminous emission lines .
hence grbs that occurred in very low - metallicity host galaxies might be systematically missed in the low - redshift sample , especially at higher - end of the redshift range where detection of weak emission lines are difficult . however
, such redshift dependent metallicity bias is not seen in our sample ( the right panel of figure  [ fig : eiso_z ] , correlation coefficient is 0.007 ) .
it is also interesting to compare the expected rate of low - redshift grbs to the actual finding rate to evaluate the success rate of identifying low - redshift grbs .
based on the largest unbiased grb survey to date ( shoals ) , @xcite estimated the event rate density of bright grbs with @xmath177 [ erg ] at @xmath9 to be @xmath178 [ yr@xmath16gpc@xmath115 ] .
considering the field of view of _ swift_/bat ( 1.4 steradian ) and the operation time until the end of march 2014 ( 9.3 years since december 2004 ) , the expected number of bright grbs detected by _
swift_/bat at @xmath12 in this time period is @xmath179 . on the other hand ,
the number of _ swift_/bat - detected grbs with @xmath177 [ erg ] in our sample is 4 ( figure  [ fig : eiso_z ] ) , which suggests that the success rate of redshift determination is higher at low - redshifts compared with that of the overall population of grbs , although the estimation of the grb rate density is still uncertain .
the faint end of the grb luminosity function is still controversial affected by sample selections and model assumptions ( e.g. , @xcite ) .
hence it is difficult to estimate the actual rate and evaluate the redshift success rate of low - luminosity grbs .    on the other hand , it is also true that redshift determination of a grb has failed ( partly ) due to the weak emission line of the host galaxy .
as mentioned in section  [ sec : obs:111225a ] , the redshift of grb 111225a has not been determined until @xcite reports detections of weak emission lines over the afterglow together with a marginal absorption line using an updated analysis tool . the failure of the redshift determination in the immediate follow up observations after the burst suggests that grb host galaxies with @xmath180 might be missed in redshift selected samples , given that the stellar mass of the grb 111225a host galaxy is naively @xmath181 ( section  [ sec : host_msfr ] ) .      another interesting case is the identification of the host galaxy of grb 020819b .
grb 020819b has been considered to be a low - redshift burst at @xmath182 since radio afterglow observations by @xcite localized the burst on a `` blob '' which resides @xmath1 from a spiral galaxy at @xmath182 , and further observations of the host galaxy by @xcite suggested that the `` blob '' also resides at @xmath182 .
however , recent extensive observations of the host galaxy by @xcite showed that the `` blob '' is actually a dusty starforming galaxy at @xmath183 , and the low - redshift spiral galaxy is possibly a foreground object unrelated to the grb .
hence it is possible that a few of the low - redshift grbs in our sample whose redshifts are obtained only via their host galaxies actually reside in high - redshift galaxies which are too faint to be observed , and currently claimed host galaxies at low - redshifts are foreground objects .
although most of the low - redshift grbs are closely located to the center of their host galaxies within @xmath184 unlike grb 020819b , it is difficult to robustly confirm that a `` low - redshift '' grb is actually at a low - redshift when the grb does not have an afterglow measured redshift or a confirmed sn association ( grb 060505 , 080517 , 060614 , 050826 , and 090417b , see table  [ tab : grbs ] ) . furthermore , among the low - redshift grbs without robust redshift confirmations , grb 060505 and 080517 have offsets @xmath185 from the center of their host , although the large offsets are not surprising given their very low - redshifts @xmath186 .
it should be noted that the originally claimed host galaxy of grb 020819b at @xmath182 was a massive - spiral galaxy with the highest metallicity in the known grb host galaxies .
when a grb host galaxy sample is contaminated by unrelated forground galaxies , it is possible that the contaminants tend to be brighter galaxies because they can be found more easily than the actual host galaxies which are likely faint .
missing redshift informations for grbs in faint galaxies and contamination by bright foreground galaxies may bias the metallicity distribution of the overall population towards higher - metallicities , in contrast to what is expected in the case of the bias against `` dark '' grbs .
we have presented the overall metallicity distribution of the host galaxies of all grbs known at @xmath12 that occurred before the end of march 2014 , including the newly obtained emission line fluxes of the host galaxies of grb 060614 , 090417b , and 130427a .
the low - redshift sample of grb host galaxies with complete metallicity measurements is essential to study the nature of grb progenitor in some aspects : 1 ) the high success rate of redshift determinations and host identifications , 2 ) the reporting bias can be eliminated by completing the metallicity measurements , and 3 ) a wealth of comparison sample is available .
we have compared the metallicity distribution of the low - redshift sample with the predictions from the empirical formulations of the properties of low - redshift galaxies .
to predict the metallicity distribution , we adopted models of grb efficiency ( @xmath187 ) as functions of metallicity .
the metallicity of the progenitor stars considered here is not necessarily identical to that of its host galaxy . instead , we considered the metallicity variation within each galaxy motivated by the observations of h regions in nearby galaxies .
the three different formulations of @xmath105 , namely the step function , power - law , and exponential models are examined .
each of the three @xmath105 models has two free parameters that can be used to fit the metallicity distribution of the host galaxies . in either case , @xmath105 function with a sharp cutoff around 12+log@xmath55(o / h)@xmath188 8.3 ( @xmath189 )
reproduces the metallicity distribution of the low - redshift sample best .
this cutoff metallicity is naively consistent with those predicted by the stellar evolution models ( @xmath1 0.1@xmath190 , @xcite ) .
@xmath105 models with moderate or weak metallicity dependence [ @xmath191 are disfavored , although the statistical significance of the current constraints is still low .
this result is in contrast to the results of some previous studies which have suggested the cutoff metallicity @xmath1 0.5@xmath192 ( e.g. , @xcite ) .
this is because we take the internal variation of metallicity within each galaxy into account .
our galaxy models also indicate the cutoff metallicity around 0.51@xmath59 , when the internal metallicity variation is assumed to be smaller ( 0.20.1 dex ) in any galaxy
. however , the current model of the internal metallicity variation relies on the observed metallicity distributions of h regions in a small number of local galaxies , and the actual internal metallicity variation in grb host galaxies is hardly known .
this uncertainty possibly affect our results significantly .    the effect of the correlation beween sfr and metallicity of galaxies on the predicted metallicity distribution is smaller than the uncertainty of metallicity measurement .
however , it may improve the agreement of the predicted and observed distribution of the grb host galaxies on the parameter plane of @xmath62 versus @xmath63 .
the relation between the nature of grb progenitor stars and the observable properties of their host galaxy depends on the internal and global properties of galaxies .
better understanding on the general population of galaxies at each redshift is essential to obtain more robust constraints on the properties of grb progenitors , in addition to gathering larger less biased sample of grb host galaxies .    because we collected grbs with known redshifts ( @xmath0 ) in this study
, it is possible that the sample of the host galaxies is biased for a kind of galaxy population in which grb redshift determination is easy ( e.g. , galaxies with smaller dust content ) , although the grb redshift determination is generally easier at lower - redshifts .
however , the number of low - redshift grbs in our sample is close to the number expected from the latest unbiased surveys , suggesting a high success rate of redshift determinations at low - redshifts , although the statistical error is still large .
furthermore , the obtained metallicities do not show significant systematic trend with respect to burst energy ( @xmath167 ) and redshift .
thus there is no sign that the low - redshift sample of the grb host galaxies is biased compared to the actual population of grb host galaxies , although the tests are not robust yet .    on the other hand , the failure of the redshift determination of grb 111225a in the immediate follow up observations after the burst suggests that faint ( and thus low - metallicity ) host galaxies might be systematically missed in the redshift selected sample .
it is also possible that a few low - redshift galaxies which are aligned in the foreground of higher - redshift grbs contaminate the low - redshift sample as in the case of grb 020819b .
the more precise estimation of the grb rate density and their properties , including those of high - redshift bursts which might contaminate low - redshift samples , will enable us to evaluate sampling biases in the redshift selected sample robustly .
we thank gemini and subaru observatory staffs for their kind support for the observations .
yn is supported by the research fellowship for young scientists from the japan society for the promotion of science ( jsps ) .
this research has made use of the ghosts database ( www.grbhosts.org ) , which is partly funded by spitzer / nasa grant rsa agreement no .
thanks are also due to nathaniel r. butler , jochen greiner , and daniel a. perley , for producing the useful online databases .
we would also like to thank the anonymous referee for helpful comments .
* * grb 980425 host galaxy * : the n2 method indicates @xmath193 and log@xmath194 = 7.3 , while the @xmath57 method does not have any valid solution .
the [ o]@xmath314959 flux is not measured in @xcite , and we assume theoretically expected line ratio of [ o]@xmath314959/[o]@xmath315007 = 1/3 ( e.g. , @xcite ) .
the lack of @xmath57 solution generally suggests @xmath195 8.18.7 , which is consistent with the n2 solution in this case .
we adopt the n2 solution as the best estimate of the metallicity of this host galaxy . * * grb 060505 host galaxy * : the n2 method indicates @xmath196 and log@xmath194 = 7.2 ( the best estimate ) .
the @xmath57 method does not have valid solution . * * grb 080517 host galaxy * : the emission line strength is presented as equivalent widths in @xcite .
we derive the emission line fluxes assuming the continuum flux density to be 0.8 , 1.7 , and 1.5 @xmath197 erg s@xmath16@xmath65@xmath16 at the wavelength of [ o ] , h@xmath32[o ] , and h@xmath33[n ] , respectively .
the n2 , @xmath57 upper - branch , and @xmath57 lower - branch solutions are @xmath198 , 8.68@xmath199 , and 8.38@xmath199 , with log@xmath194 = 7.2 , 7.2 , and 7.1 , respectively .
we adopt the @xmath57 upper - branch solution as the best estimate . * * grb 031203 host galaxy * : the n2 method indicates @xmath200 ( log@xmath194 = 8.6 ) , and the @xmath57 method provides two solutions of @xmath201 and 8.32@xmath199 ( so called upper-/lower - branches , with log@xmath194 = 8.3 and 8.2 , respectively ) .
+ kk04 calibrated the metallicity indicators in a range of ionization parameter @xmath202 , and we need to extrapolate the relation when log@xmath194 is outside of this range . because the n2 index is sensitive to the ionization parameter , we do not consider the n2 solution for this host galaxy reliable , although the indicated super - solar metallicity suggests the upper - branch solution of the @xmath57 method is more likely than the lower - branch solution .
thus we adopt the @xmath57 upper - branch solution as the best estimate .
the @xmath57 upper - branch is less sensitive to the ionization parameter than the lower - branch and the n2 .
+ although the position of this galaxy on the bpt diagram is consistent with being a star forming galaxy ( figure  [ fig : bpt ] ) , @xcite showed that the host galaxy might have an agn .
the high @xmath60 parameter value also supports the existence of an agn in the galaxy . * * grb 060614 host galaxy * : due to the non - detection of [ n ] line , we can put only upper - limit on the n2 index ( n2 @xmath203 ) .
the [ o]@xmath314959 line coincides with an strong night - sky emission line , and the flux can not be measured from the observed spectra .
hence we assume the theoretically predicted line ratio of [ o]@xmath314959/[o]@xmath315007 = 1/3 .
the upper- and lower - branch solutions of the @xmath57 method are @xmath204 and 8.35@xmath199 with log@xmath194 = 7.4 and 7.3 , respectively .
assuming log@xmath194 = 7.37.4 , the n2 upper - limit indicates @xmath205 .
thus we adopt the @xmath57 lower - branch solution as the best estimate . *
* grb 030329 host galaxy * : the n2 , @xmath57 upper - branch , and @xmath57 lower - branch solutions are @xmath206 , 8.75@xmath199 , and 8.14@xmath199 , with log@xmath194 = 7.8 , 8.2 , and 7.8 , respectively .
we adopt the @xmath57 lower - branch solution . *
* grb 120422a host galaxy * : the n2 method indicates @xmath207 and log@xmath194 = 7.3 ( the best estimate ) .
the @xmath57 method does not have valid solution . * * grb 050826 host galaxy * : the n2 , @xmath57 upper - branch , and @xmath57 lower - branch solutions are @xmath208 , 8.86@xmath199 , and 8.14@xmath199 , with log@xmath194 = 7.5 , 7.5 , and 7.3 , respectively .
we adopt the @xmath57 upper - branch solution . * * grb 130427a host galaxy * : the n2 , @xmath57 upper - branch , and @xmath57 lower - branch solutions are @xmath209 , 8.67@xmath199 , and 8.33@xmath199 , with log@xmath194 = 7.4 , 7.4 , and 7.3 , respectively .
we adopt the @xmath57 upper - branch solution as the best estimate . however , we note that the focas too spectrum from which we measure the emission line fluxes of [ o]@xmath313727 , h@xmath32 , and [ o]@xmath2104959,5007 , is obtained without adc .
hence it is possible that the line fluxes suffer from differential atmospheric refraction . + to evaluate the effect of this systematic uncertainty , we also examine the independently measured emission line fluxes reported in k15 ( shown in table  [ tab : lines ] ) .
the k15 line fluxes are systematically weaker than those in the gmos+focas spectra possibly due to larger loss of light at the slit , however we can compare relative flux ratios between the lines .
+ with the k15 line fluxes , the n2 , @xmath57 upper - branch , and @xmath57 lower - branch solutions are @xmath211 , 8.70@xmath212 , and 8.34@xmath213 .
we note that the detection of the [ n ] line is marginal in k15 ( s / n = 2.5 ) .
although it is difficult to determine which of the @xmath57 solutions is true only with the k15 spectrum , if we adopt the upper - branch solution as suggested by our gmos spectrum , it agrees well with the estimate from the gmos+focas spectra . * * grb 090417b host galaxy * : the [ o ] , h@xmath32 , and [ o ] line fluxes are not known .
assuming @xmath214 log@xmath215 ( the @xmath60 parameter range covered by the other grb host galaxies in the sample excluding that of grb 031203 ) , the n2 method indicates @xmath216 = 8.68.9 . * * grb 061021 host galaxy * : the [ n ] line is not detected providing n2 @xmath217 .
the @xmath57 upper - branch and lower - branch solutions are @xmath218 and @xmath219 , respectively ( log@xmath194 = 7.5 in either case ) .
the n2 upper - limit is not deep enough to reject either of the @xmath57 solutions .
combining the @xmath57 upper - branch and lower - branch solutions , we estimate the metallicity of this host galaxy to be 8.288.66 . * * grb 011121 host galaxy * : the n2 , @xmath57 upper - branch , and @xmath57 lower - branch solutions are @xmath220 , 8.69@xmath199 , and 8.33@xmath199 , with log@xmath194 = 7.2 , 7.4 , and 7.2 , respectively . we adopt the @xmath57 lower - branch solution . * * grb 120714b host galaxy * : the n2 method indicates @xmath221 , although the detection of the [ n ] line is marginal . the @xmath57 method , which is not affected by the low - s / n of the [ n ] line , provides the upper - branch and the lower - branch solutions : @xmath222 and 8.43@xmath199 , respectively ( log@xmath194 = 7.5 in either case ) .
we adopt the @xmath57 upper - branch solution , which agrees with the n2 solution .