Abstract
We present the stellar population properties of 69 short gamma-ray burst (GRB) host galaxies, representing the largest uniformly modeled sample to date. Using the Prospector stellar population inference code, we jointly fit photometry and/or spectroscopy of each host galaxy. We find a population median redshift of (68% confidence), including nine photometric redshifts at z ≳ 1. We further find a median mass-weighted age of tm =
Gyr, stellar mass of log(M*/M⊙) =
, star formation rate of SFR =
M⊙ yr−1, stellar metallicity of log(Z*/Z⊙) =
, and dust attenuation of
mag (68% confidence). Overall, the majority of short GRB hosts are star-forming (≈84%), with small fractions that are either transitioning (≈6%) or quiescent (≈10%); however, we observe a much larger fraction (≈40%) of quiescent and transitioning hosts at z ≲ 0.25, commensurate with galaxy evolution. We find that short GRB hosts populate the star-forming main sequence of normal field galaxies, but do not include as many high-mass galaxies as the general galaxy population, implying that their binary neutron star (BNS) merger progenitors are dependent on a combination of host star formation and stellar mass. The distribution of ages and redshifts implies a broad delay-time distribution, with a fast-merging channel at z > 1 and a decreased neutron star binary formation efficiency from high to low redshifts. If short GRB hosts are representative of BNS merger hosts within the horizon of current gravitational wave detectors, these results can inform future searches for electromagnetic counterparts. All of the data and modeling products are available on the Broadband Repository for Investigating Gamma-ray burst Host Traits website.
Export citation and abstract BibTeX RIS
1. Introduction
Short-duration gamma-ray bursts (GRBs) are some of the most luminous explosions in the universe and are observed over a large range of cosmological distances (z ≈ 0.1–2.5; Selsing et al. 2018; Paterson et al. 2020). Host galaxy associations to short GRBs have been pivotal in uncovering the true nature of GRB progenitors, as they provide redshifts and information on their stellar populations and the types of stars that produced them. A few of the first short GRBs detected soon after the launch of the Neil Gehrels Swift Observatory (Swift; Gehrels et al. 2004) originated in older, quiescent host galaxies (GRBs 050509B and 050724; Berger et al. 2005; Gehrels et al. 2005; Bloom et al. 2006), implying that short GRBs are derived from older stellar progenitors than long-duration GRBs, with an assumed long formation timescale. However, within a few years, studies showed that a majority of short GRBs (≈75%) were associated with younger, star-forming galaxies, indicating that their progenitors have a spread of formation timescales (Leibler & Berger 2010; Fong & Berger 2013; Berger 2014). Coupled with their energy scales (Gehrels et al. 2008), inferred event rates (Fong et al. 2015), locations within their host galaxies (Fong & Berger 2013; Tunnicliffe et al. 2014), lack of associated supernovae (SNe; e.g., Fox et al. 2005; Hjorth et al. 2005; Soderberg et al. 2006), and claimed detections of the first r-process kilonovae (Berger et al. 2013a; Tanvir et al. 2013), there is a wealth of circumstantial evidence that short GRBs originate from compact object mergers, specifically either in binary neutron star (BNS) or neutron star black hole (NSBH) mergers.
The LIGO/Virgo discovery of gravitational wave (GW) BNS merger GW170817 and its associated electromagnetic (EM) counterpart, short GRB 170817A (Abbott et al. 2017a; Goldstein et al. 2017; Savchenko et al. 2017), provided a definitive connection between short GRBs and compact object binary mergers for the first time. The host of GW170817/GRB 170817A, NGC4993, was the oldest by several gigayears and one of the most quiescent galaxies compared to other short GRB hosts and the general field galaxy population (Blanchard et al. 2017; Levan et al. 2017; Palmese et al. 2017). Due to the older stellar population in this host, the "delay time" between the formation of the binary system and the BNS merger was likely several gigayears. The characterization of NGC4993, together with the cosmological population of short GRB host environments, motivates us to ask the following questions:
- 1.What is the relationship between short GRB progenitors and recent star formation and stellar mass in galaxies?
- 2.What are short GRB formation channels and timescales?
- 3.What are the differences between the environments of cosmological short GRBs and those of local universe BNS/NSBH mergers?
To fully understand the short GRB progenitor from formation to merger and the environmental factors that are at play, it is imperative to characterize the properties of short GRB host galaxies over cosmic time. For example, the stellar population ages and redshifts of the host galaxies serve as vital constraints on the delay-time distribution (DTD), while star formation rates (SFRs) and stellar masses describe the conditions under which BNS/NSBH mergers form. Previous studies determined a peak in redshift at z ≈ 0.5 (Berger 2014), with detection of just a few short GRBs existing at z > 1.0 (Paterson et al. 2020), where one may expect to find more short delay-time progenitors. More recently, O'Connor et al. (2022) suggested, based on photometric redshift estimates, that there is a larger sample of short GRBs at z > 1.0. Additional studies have deduced a median stellar population age of ∼1 Gyr, and SFRs and stellar masses that span a wide range (Berger 2009, 2014; Leibler & Berger 2010; Fong & Berger 2013; Nugent et al. 2020), indicating no clear trends with a single galaxy property. This is unlike other explosive transients, which originate from massive stars (e.g., long GRBs and superluminous supernovae, hereafter SLSNe) and have been shown to occur in low-metallicity, low-mass, young, and actively star-forming galaxies (Svensson et al. 2010; Perley et al. 2013, 2016; Wang & Dai 2014; Vergani et al. 2015; Niino et al. 2017; Hatsukade et al. 2018; Schulze et al. 2021). Type Ia SNe, which have older white dwarf progenitors, are found in both star-forming and quiescent galaxies and generally trace the properties of the field galaxy population (Oemler & Tinsley 1979; Mannucci et al. 2005; Sullivan et al. 2006; Pan et al. 2014; Wiseman et al. 2021).
By comparing short GRB host environments to the general field galaxy population, we can also determine if their progenitors have preferences for specific galaxy traits. For instance, the short GRB host population and its adherence to standard galaxy relations, including the star-forming main sequence (SFMS) and mass–metallicity relation (Gallazzi et al. 2005; Speagle et al. 2014; Whitaker et al. 2014; Leja et al. 2022), can offer insight on whether or not short GRB hosts are unusual. The frequency of short GRB host types is also informative, as we expect ≈80% of detectable field galaxies out to z = 0.3 to be star-forming (Dressler 1980; Martis et al. 2016; Leja et al. 2022); thus, adherence or deviation from this fraction tells us how dependent the progenitor is on recent star formation.
A uniform and holistic study of short GRB stellar populations, which takes advantage of nearly two decades of short GRB host identifications, has yet to be conducted in a systematic way. Given the steady flow of well-localized cosmological short GRBs, and the present era of GW discovery, it is timely to conduct a large sample study. Over the past few years, stellar population modeling has also rapidly progressed to include more sophisticated star formation histories (SFHs), informed statistics, and the ability to jointly fit photometric and spectroscopic data (e.g., Nugent et al. 2020; Johnson et al. 2021).
Here, we present the stellar population properties of 69 confidently associated short GRB host galaxies. This represents the second in a series of two papers. Paper I, Fong et al. (2022), focuses on the photometric and spectroscopic catalogs, host galaxy associations, spectroscopic redshifts, and galactocentric offsets. This paper, Paper II, focuses on spectral energy distribution (SED) modeling of these data, their inferred stellar population properties, and implications for the progenitors. In Section 2, we describe the host sample and how hosts were selected for this study. We outline our stellar population modeling methods in Section 3. We discuss the inferred stellar population properties and redshift distribution in Section 4, and compare short GRB host properties to several known galaxy relations, including the SFMS, mass-weighted Schechter function, and the mass–metallicity relation. We compare properties intrinsic to the GRB, including their offsets, and their host in Section 5. We examine short GRB hosts in the context of other transient hosts and the host of GW170817/GRB 170817 in Section 6. We list several possible sample biases and potential missing populations in Section 7. Finally, our conclusions are given in Section 8. We house all of the data and modeling products described in this work on the Broadband Repository for Investigating Gamma-ray burst Host Traits (BRIGHT) website. 15
Unless otherwise stated, all observations are reported in the AB magnitude system and have been corrected for Galactic extinction in the direction of the GRB (Cardelli et al. 1989; Schlafly & Finkbeiner 2011). We employ a standard WMAP9 cosmology of H0 = 69.6 km s−1 Mpc−1, ΩM = 0.286, and Ωvac = 0.714 (Hinshaw et al. 2013; Bennett et al. 2014).
2. Short GRB Host Sample
We use the short GRB host sample described in Fong et al. (2022). 16 This sample contains 84 short GRB hosts with broadband photometric detections and upper limits; 58% of the sample has spectroscopic redshifts and 50% of the sample has available spectra with spectral line detections. Host associations in this sample are defined through the probability of chance coincidence method (Pcc; Bloom et al. 2002), which uses an optical magnitude of the host compared to the distance from the host to the short GRB X-ray, optical, or radio afterglow position. A low Pcc translates to a higher likelihood that the short GRB is associated with the galaxy. For this study, we only consider hosts where Pcc ≲ 0.20, above which we cannot conclusively assign a host galaxy to a burst. Furthermore, as the goal of this study is to model the SEDs of the host galaxies, we only use the hosts that are detected in ≥3 photometric bands. Thus, our sample here comprises 69 out of the 84 host associations found in Fong et al. (2022).
From here on, we define the "Gold Sample" as that with the most robust host associations (Pcc < 0.02). Similarly, the "Silver Sample" (0.02 < Pcc < 0.09) and the "Bronze Sample" (0.09 < Pcc < 0.20) represent moderately robust associations (see Fong et al. 2022 for a more detailed explanation). We also classify hosts based on whether they have photometric (20 hosts) or spectroscopic redshifts (49 hosts). There are 41 hosts (31 with spectroscopic redshifts) in the Gold Sample, 19 hosts (11 with spectroscopic redshifts) in the Silver Sample, and nine (seven with spectroscopic redshifts) hosts in the Bronze Sample. We show this breakdown in Figure 1.
Figure 1. The fractions of star-forming (SF; blue), transitioning (T; yellow), and quiescent (Q; red) short GRB host galaxies based on the results in Table 2 for several subgroups of the sample. We see that in the full sample, a majority (∼84%) of short GRBs occur in star-forming galaxies, and ∼16% occur in transitioning or quiescent galaxies. These fractions also stay fairly consistent for the Gold, Silver, and Bronze samples. We find a lower star-forming fraction for the spectroscopic redshift sample than the photometric redshift sample, likely because the photometric redshift sample reaches higher redshifts, where star-forming galaxies are much more prevalent.
Download figure:
Standard image High-resolution image3. Stellar Population Modeling
To determine the stellar population properties of the host galaxies, we use the Python-based SED modeling code, Prospector (Leja et al. 2017; Johnson et al. 2021). We jointly fit photometric and spectroscopic data (when available) to determine properties including redshift (when not already determined from the spectrum), stellar mass, stellar population age, gas and stellar metallicities, and dust attenuation (see Table 1 for a review of fit properties). We apply the nested sampling fitting routine dynesty (Speagle 2020) within Prospector to all available observational data to produce posterior distributions of the sampled properties. Prospector applies the MIST (Paxton et al. 2018) models and MILES spectral library (Falcón-Barroso et al. 2011) through FSPS (Flexible Stellar Population Synthesis) and python-FSPS to produce model photometry and spectra (Conroy et al. 2009; Conroy & Gunn 2010).
Table 1. Prospector Parameters and Prior Distributions
Parameter | Definition | Prior |
---|---|---|
FITTING PARAMETERS | ||
z | redshift |
![]() |
tage (Gyr) | age of the galaxy at the time of observation |
![]() |
τ | e-folding time of delayed-τ SFH |
![]() |
![]() | total mass formed |
![]() |
![]() | stellar metallicity |
![]() |
τV,1 | optical depth of young stellar light |
![]() |
τV,2 | optical depth of old stellar light |
![]() |
N0 | spectral normalization factor |
![]() |
![]() | gas-phase metallicity |
![]() |
Ugas | gas ionization parameter |
![]() |
τAGN | mid-IR optical depth |
![]() |
fAGN | fraction of AGN luminosity in galaxy |
![]() |
Njitter | noise inflation factor for spectra |
![]() |
DERIVED PARAMETERS | ||
tm | mass-weighted age in gigayears |
tm
= ![]() |
![]() | stellar mass |
M* ≈ ![]() |
SFR | star formation rate in M⊙ yr−1 | SFR(t) = ![]() |
AV | total dust attenuation in mag | AV = 1.086 × (τV,1 + τV,2) |
Note. A list of all possible free parameters used in our Prospector fits, their definitions, and their prior distributions. represents a uniform distribution, while
represents a normal distribution. We also list the derived parameters used in this analysis. We note that z is only set free for hosts that do not have a known redshift, N0,
, and Ugas are only used in fits including spectroscopy, the active galactic nucleus (AGN) parameters were only included for one host (GRB 150101B) that has a known AGN, and Njitter was only applied to hosts with high signal-to-noise ratio (S/N) spectra.
Download table as: ASCIITypeset image
For all Prospector fits, we use a Chabrier initial mass function (Chabrier 2003), Milky Way extinction law (Cardelli et al. 1989), a parametric delayed-τ SFH (SFH ∝ t ∗ e−t/τ ) defined by the e-folding factor τ, and include the effects of nebular emission and line strengths through a nebular emission model (Byler et al. 2017). We choose a Milky Way extinction law, as it has been shown to model the majority of galaxy SED shapes well, and moreover the majority of the hosts in this sample. We also include the Gallazzi et al. (2005) mass–metallicity relation of galaxies in the fits, which ensures realistic mass values are being sampled for a given sampled metallicity. For star-forming hosts, we determine the dust attenuation from old (≥107 yr; Conroy & Gunn 2010) and young stars, respectively, through the commonly used 2:1 ratio (Calzetti et al. 2000; Price et al. 2014; Leja et al. 2019), as younger stars tend to attenuate twice the amount of dust as old stars. For known quiescent hosts, we assume there is no dust attenuation of young stellar light (parameter τV,1 = 0).
The observed host photometry is modeled in Prospector by integrating over the wavelength coverage of the filter transmission curves. For most host galaxies, we use the standard Sloan Digital Sky Survey (griz; Doi et al. 2010), Bessell (BVRI), Two Micron All Sky Survey (2MASS; JHK; Skrutskie et al. 2006), Spitzer, Wide-field Infrared Survey Explorer (WISE; Wright et al. 2010), and Wide Field Camera 3 WFC3/IR and WFC3/UVIS (Dressel 2012) transmission curves, available in the SedPy Python package (10.5281/zenodo.4582722). We apply telescope-specific filter transmission curves when photometry is from the Keck Observatory (LRIS, DEIMOS, MOSFIRE), the MMT Observatory (MMIRS), and UKIRT (WFCAM), as these differ from the standard filter sets. The instruments used for the photometric data of the hosts are reported in Table 1 in Fong et al. (2022).
For all hosts, the main fitted parameters are redshift (z) if it is not known, mass formed (MF ), age of the galaxy at the time of observation (tage), the SFH e-folding factor τ, optical depth due to dust of young (τV,1) and old (τV,2) stellar light, and stellar metallicity (Z*; see Table 1). We perform three main types of fits depending on the data available, each of which contains unique parameter specifications and sometimes added stellar population properties: (i) joint photometric and spectroscopic fits with known redshifts (see Figure 2, left panel); (ii) photometric fits with known redshifts; and (iii) photometric fits with no known redshifts (see Figure 2, right panel).
Figure 2. Representative fits from two of the 69 short GRBs in our sample. Prospector produced model spectra (purple lines) and photometry (purple squares) compared to the observed spectra (green line) and error (gray line) and photometry (green circles) for the joint spectroscopy and photometry fit of GRB 150728A (left panel) and the redshift-free, photometry-only fit of GRB 160411A (right panel). We highlight these fits to show the accuracy in measurement of the spectral continuum and line locations (GRB 150728A) and the photometric colors for redshift determination (GRB 160411A). In total, we have performed 40 joint spectroscopic and photometric fits and 29 photometry-only fits. These modeling products are all available on the BRIGHT website.
Download figure:
Standard image High-resolution imageFirst, we perform joint spectroscopic and photometric Prospector fits for all hosts with at least one spectral line detection of signal-to-noise ratio (S/N) ≥5 (we find diminishing returns for host galaxies with lower S/Ns). In our sample, we have 40 hosts with spectra that meet these conditions (for an example, see Figure 2, left panel). We determine redshifts of these hosts from spectral line detection (see Fong et al. 2022) and fix the z-parameter to the spectroscopic redshift and the maximum of tage to the lookback time at that redshift (tlookback(z)). We include additional free parameters that are dependent on spectral line detection in these fits to better estimate the observed spectral line strengths and continuum: gas-phase metallicity (Zgas), the dimensionless gas ionization parameter that measures the ratio of hydrogen ionizing photons density to hydrogen density (Ugas), and a parameter that normalizes the model to the observed continuum (N0). The N0 parameter is also used to marginalize over flux calibration uncertainties and slit losses. The model spectral continuum is built from an nth-order Chebyshev polynomial. 17 In addition, for some hosts we apply a spectral noise inflation model to properly weight the photometric observations against the high-S/N spectrum. 18 The host of GRB 150101B contains a known active galactic nucleus (AGN; Fong et al. 2016; Xie et al. 2016). Thus, we also include several AGN parameters to properly account for the optical depth in the mid-IR (tAGN) and fraction of the AGN luminosity in the galaxy (fAGN).
The second type of fit is for short GRB hosts with known spectroscopic redshifts, but no available or very low-S/N spectra for the Prospector fits. For a majority of these hosts, we found spectroscopic redshifts in previously published works, but could not find publicly available spectra. For both of the hosts of GRBs 140930B and 181123B (Paterson et al. 2020), their spectral continua have low S/N and are not meaningful in the Prospector fits. Thus, we perform a photometry-only fit for these hosts, where we set the z- and maximum tage parameters in the same manner as the joint spectroscopic and photometric fits. Furthermore, we fix Zgas = Z⊙ and Ugas = −2.0 and remove nebular emission lines in the model spectra, as there are no spectral lines to fit (Leja et al. 2017).
The third type of fit is for the 20 hosts with unknown redshifts; in these cases, spectroscopy is not feasible due to the apparent faintness of the hosts or there were no spectral lines detected in their observed spectrum (Figure 2, right panel). We therefore leave the redshift as a free parameter and allow it to range uniformly between 0.1 ≤ z ≤ 3.0. We choose z = 3.0 as the maximum possible redshift, as Swift's sensitivity to short GRBs steeply drops beyond z ≈ 1, likely due to detectability of the bursts' luminosity (Lien et al. 2016). Additionally, we set the maximum of tage to be the maximum age of the universe at the sampled redshift. For the hosts of GRBs 130515A, 160411A, and 180418A, we see large increases in flux between two photometric bands (r and i, z and J, and z and J, respectively), which is a clear indication of the 4000 Å break. Therefore, we tighten their redshift ranges to only allow for redshifts that give a 4000 Å break within those wavelength ranges. For GRB 180418A, we have further knowledge of its possible redshift (z > 1.0) from Rouco Escorial et al. (2021). For the host of GRB 210726A, we restrict the redshift range to z < 1.0, as solutions at z > 1.0 violate a deep U-band upper limit and its afterglow luminosity suggests a less likely z > 1.0 origin (G. Schroeder et al. 2022, in preparation). For the hosts of GRBs 170127B and 210726A , we use the upper limits in the J and K bands and the U and J bands, respectively, in the fits to better constrain photometric redshift estimate.
As a point of comparison, we also fit the photometric data of NGC4993 (the host of GRB 170817A). We fit the GALEX, PS1, 2MASS, and WISE photometric data (applying the respective filter transmission curves) and redshift of NGC4993 listed in Blanchard et al. (2017), in which a nonparametric SFH Prospector model fit was achieved. Parametric fits and results are preferred for this work, as they better establish uniformity among data sets with an inconsistent amount and quality of data. We also note that typically, nonparametric fits result in mass values that are 25%–100% larger and age values that are three to five times older (Leja et al. 2019). Although a Prospector fit with a nonparametric SFH was performed in Blanchard et al. (2017), for consistency, we perform a parametric SFH fit using the specifications for the photometry and known redshift fits. We do not include NGC4993 as a host in our catalog, as the GRB was discovered and associated in a different way than the rest of the host population.
For a direct comparison to commonly used physical properties of galaxies, we derive a number of properties from our fitting parameters (Table 1). For instance, we convert mass formed (MF ) to stellar mass (M*) using Equation (3) in Nugent et al. (2020) and referenced in Table 1. We also use the posteriors on tage and τ to derive the posterior on the mass-weighted age (tm ). Mass-weighted ages importantly do not overestimate the contribution from young stars, a major caveat of light-weighted or single stellar population ages (Conroy 2013). We derive the present-day SFR from the posteriors on tage, τ, and MF . Finally, we derive the total dust attenuation AV from the optical depths: τV,1 and τV,2.
We list all possible sampled properties, their prior distributions, and the derived parameters M*, tm , SFR, and AV in Table 1. For each host, we report the median and 68% credible interval of the posterior in several relevant stellar population properties in Table 2. We report upper limits when the 99.7% credible interval is consistent with the lower limit of the property range.
Table 2. Short GRB Host Galaxy Stellar Population Properties
GRB | z | Galaxy Type | tm (Gyr) | log(M*/M⊙) | SFR (M⊙ yr−1) | log(sSFR) (yr−1) | log(Z*/Z⊙) | AV (mag) | Fit Type |
---|---|---|---|---|---|---|---|---|---|
GOLD SAMPLE | |||||||||
050509B | 0.225 | Q |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
050709 | 0.161 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
050724 | 0.257 | Q |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
051221A | 0.546 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
060614 | 0.125 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
060801 | 1.13 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
061006 | 0.461 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
061210 | 0.41 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
070429B | 0.902 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
070714B | 0.923 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
070724 | 0.457 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
070809 | 0.473 | T |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
071227 | 0.381 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
090510 | 0.903 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
100117 | 0.914 | Q |
![]() |
![]() |
![]() | <−10.80 |
![]() |
![]() | S |
100206A | 0.407 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
101224A | 0.454 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
120804A |
![]() | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
121226A |
![]() | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
130603B | 0.357 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
140129B | 0.43 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
140903A | 0.353 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
141212A | 0.596 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
150101B | 0.134 | Q |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
150120A | 0.46 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
150728A | 0.461 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
160411A |
![]() | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
160525B |
![]() | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
170428A | 0.453 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
170728B | 1.272 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
170817 b | 0.01 | Q |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
180418A |
![]() | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
180618A |
![]() | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
180727A |
![]() | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
181123B | 1.754 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
200219A |
![]() | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
200522A | 0.554 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
200907B |
![]() | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
210323A | 0.733 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
210726A |
![]() | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
211023B | 0.862 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
211211A | 0.076 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
SILVER SAMPLE | |||||||||
051210 |
![]() | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
070729 |
![]() | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
090515 | 0.403 | Q |
![]() |
![]() |
![]() | <−11.64 |
![]() |
![]() | P |
100625A | 0.452 | Q |
![]() |
![]() |
![]() | <−11.64 |
![]() |
![]() | S |
101219A | 0.718 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
111117A | 2.211 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
120305A | 0.225 | T |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
130515A |
![]() | T |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
130822A | 0.154 | T |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
140930B | 1.465 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
150831A | 1.18 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
151229A |
![]() | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
160303A |
![]() | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
160624A | 0.484 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
160821B | 0.162 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
161001A |
![]() | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
170127B |
![]() | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
180805B | 0.661 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
191031D |
![]() | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
BRONZE SAMPLE | |||||||||
050813 | 0.719 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
080123 | 0.495 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
140622A | 0.959 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
160408A |
![]() | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
161104A | 0.793 | Q |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
170728A | 1.493 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
200411A |
![]() | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | P |
201221D | 1.055 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
210919A | 0.242 | SF |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() | S |
Notes. We present the medians and 68% credible region, or 99.7th percentile upper limits of the Prospector-derived stellar population properties for all 69 short GRB hosts in our sample, and include the host of GW/GRB 170817 for reference. All values, except redshift when it is known, are set free in the fits. Spectroscopic redshifts are detailed in Fong et al. (2022). Hosts with spectra used in the Prospector fits have fit type "S," whereas hosts that only have photometry available for fits have fit type "P." The Gold Sample hosts have the lowest Pcc , and Bronze Sample hosts have the highest Pcc .
a Photometric redshifts determined through Prospector. b Included as a point of comparison, but not included in the short GRB host catalog.A machine-readable version of the table is available.
4. Stellar Population Properties
4.1. Redshift Distributions
Here, we discuss the results of our Prospector fitting of 69 short GRB host galaxies (Tables 2–3). We show the cumulative and posterior distributions in redshift in Figure 3. We report 20 photometric redshifts, 10 of which are completely new results (we show the photometric redshift SEDs in Appendix A, Figure A1). We build the cumulative distribution function (CDF) for the entire host sample from 5000 random draws from individual host Prospector-derived posterior distributions. We include the redshifts of GRBs 090426A (z = 2.609), 160410A (z = 1.717), and 150423A (z = 1.394) in our sample as spectroscopic redshifts, as all bursts have known redshifts from their afterglows (Antonelli et al. 2009; Fong et al. 2013; Selsing et al. 2019; Agüí Fernández et al. 2021). To build the CDF and account for the individual measurement uncertainties, including both spectroscopic and photometric redshifts, we perform 5000 draws of their respective redshift values (for spectroscopic redshifts) or from the Prospector posteriors (for photometric redshifts). We choose 5000 draws as there is no change in the results if we increase the number draws beyond 5000. From the 5000 CDFs, we then determine the median CDF from the median value of the realizations in each redshift bin. We find that the spectroscopic redshifts range over 0.1 ≤ z ≤ 2.2, with a median of (68% confidence interval; Figure 3). The addition of the photometric redshift extends the redshift range to z = 2.8, and increases the median to z = 0.64. This is not surprising, as the photometric redshifts tend to capture z > 1.0 (with a median z ∼ 0.93). This is the so-called "redshift desert," a region for which it is difficult to obtain spectroscopic redshifts given that the primary identifiable spectral lines are shifted to infrared wavelengths, and these more distant sources are also apparently faint.
Figure 3. Left: the cumulative distribution function (CDF) in redshift for the Gold, Silver, and Bronze samples, as well as 5000 realizations of the CDF based on random draws from the individual posterior distributions. The thicker lines represent the median CDF while the thin lines show the individual realizations to demonstrate the sample spread. We find that there is an increasing redshift median from the Gold to Bronze Samples, as indicated by the arrows. Middle: the CDF and 5000 realizations in redshift of all short GRB hosts in our sample (purple), along with the photometric redshift (blue) and spectroscopic redshift sample (green). Right: redshift histograms of all hosts (purple), spectroscopic redshift hosts (green), and photometric redshift hosts (blue). The "All" hosts and photometric redshift distributions are built from the individual posterior distributions as described in Section 4.1. The photometric sample hosts most notably fill in the "redshift desert" (z > 1.0), and the photometric redshift sample has a median ∼0.5 greater than the spectroscopic redshift sample.
Download figure:
Standard image High-resolution imageTable 3. Stellar Population Properties for Short GRB Host Subsamples
Sample | z | tm (Gyr) | log(M*/M⊙) | SFR (M⊙ yr−1) | log(sSFR) (yr−1) | log(Z*/Z⊙) | AV (mag) |
---|---|---|---|---|---|---|---|
ALL |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
GOLD |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
SILVER |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
BRONZE |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
SPECTROSCOPIC REDSHIFT |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
PHOTOMETRIC REDSHIFT |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
STAR-FORMING |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
QUIESCENT/TRANSITIONING |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
Note. The median and 68% credible interval stellar population values for the full sample of short GRB hosts, as well as specific subsamples of hosts based on Pcc association (Gold, Silver, Bronze), redshift determination (spectroscopic or photometric), and host type (star-forming or transitioning and quiescent).
Download table as: ASCIITypeset image
Given that 10 of our photometric redshifts are at z > 1, we also want to ensure that the photometric redshifts are not simply sampling the redshift prior (see Table 1), which could skew them toward higher redshifts. We compare all photometric redshift posterior distributions to their respective prior distributions (which are uniform in linear space) and run an Anderson–Darling (AD) test, with the null hypothesis that both are derived from the same underlying distribution. The AD test best captures the overall shape of the distributions, including any tails. If the p-value of the null hypothesis PAD < 0.05, we can reject the null hypothesis. We find that for all short GRB hosts with photometric redshifts, the PAD ≲ 0.001, implying the posterior distributions are not derived from their uniform priors. Thus, the fact that the photometric redshifts are higher than the spectroscopic sample is a real effect and not simply a product of the more limited data. We also ran redshift-free fits for several to test the photometric redshift capabilities of Prospector. For the host of GRB 140930B (which had a very low S/N detection of [O ii]), we find a photometric redshift of , which is only 0.09 less than its spectroscopic redshift, although with large error. For the host of GRB 150831A, we find a photometric redshift of
, which is 0.09 lower than the true redshift, but within error. Thus, we infer that Prospector generally finds a photometric redshift of the host within its 68% credible region of the distribution, and is able to capture the low- or high-redshift nature of the host. Our analysis also highlights the importance of including both spectroscopic and photometric redshifts in the full redshift distribution of short GRB hosts. We also find an increasing median redshift between the Gold, Silver, and Bronze samples, as shown in Figures 3 and 4, although they are consistent within the 68% confidence intervals. This is likely caused by the fact that brighter, lower-redshift short GRBs are easier to associate to host galaxies than much fainter, higher-redshift galaxies, as the Pcc test, in part, depends on the brightness of the host (for more, see Fong et al. 2022).
Figure 4. From top to bottom: median (star) and 8% confidence interval values (line) for z, tm
, , and SFR for the Gold, Silver, and Bronze Samples. There is little difference in the median and 68% confidence interval for the properties in each sample; however, redshift does appear to increase from the Gold to Bronze Samples. Fong et al. (2022) also showed that there are weaker associations to hosts as redshift increases, due to faintness of the galaxies.
Download figure:
Standard image High-resolution imageIn Figure 5, we show the cumulative and posterior distributions for several stellar population properties. The CDFs are derived in the same manner as the redshift distributions. For the entire host sample, we find median stellar population property values of log(M*/M⊙) = ,
Gyr, SFR =
M⊙ yr−1,
, and
mag. Notably, we find that these properties do not significantly change among the Gold, Silver, and Bronze Samples (Figure 4), demonstrating that they are fairly impervious to the confidence in host association. We note that the Gold Sample spans a broader range in mass-weighted age, while the Bronze Sample has a slightly elevated SFR (Figure 4); these can naturally be explained as manifestations of the aforementioned redshift differences between the samples.
Figure 5. From left to right, top to bottom: CDFs of all (purple), star-forming (blue), and transitioning and quiescent (pink) short GRB hosts for log(M*/M⊙), tm in Gyr, log(SFR) in M⊙ yr−1, log(Z*/Z⊙), and AV in mag. Thick solid lines represent medians on the CDF, while 5000 realizations of the CDF in each distribution are plotted as thin lines to demonstrate the spread in each distribution. The probability distributions of all hosts are shown in the panel above the CDFs. Medians in each distribution are denoted by color-coded arrows from the bottom, and are also denoted in the legend. We find that the median short GRB hosts tend to be massive (∼109–1010 M⊙), moderately old (∼0.8 Gyr), albeit with large difference in ages between the star-forming and transitioning/quiescent hosts, have moderately low ongoing SFR (∼1.5 M⊙ yr−1), have a uniform distribution of stellar metallicities, and have low dust attenuation (∼0.4 mag).
Download figure:
Standard image High-resolution image4.2. Star Formation Classification
In order to put the short GRB host sample in context with a basic classification scheme, we systematically classify each host by the degree of star formation. We use the canonical galaxy classifications: "star-forming" comprises galaxies on the SFMS, "transitioning" comprises galaxies transitioning off the SFMS, and "quiescent," comprises galaxies off the SFMS. For this classification, we use the specific SFR (sSFR = SFR/M* in units of yr−1) and the redshift. We determine this classification using Equation (2) in Tacchella et al. (2022): , where tH
(z) is the Hubble time at redshift z. For hosts with photometric redshifts, we find
at every sampled z and sSFR. We use the following divisions for classification, following the methods of Tacchella et al. (2022):
describes star-forming galaxies,
defines the galaxies transitioning from star-forming to quiescent, and
represents quiescent galaxies. We define the hosts using the mode of the classification for their distribution of
values.
In Figure 1, we present the percentage of galaxies in each classification in the entire host sample, as well as the Gold, Silver, Bronze, spectroscopic, and photometric redshift samples. The majority of hosts in the entire sample are classified as star-forming galaxies (∼84%), while ∼6% are transitioning and ∼10% are quiescent. We draw errors on the fractions for the various subsamples to determine their statistical significance by randomly drawing "galaxies" from the full sample (choosing the number of galaxies based on the size of the group), calculate the star-forming fraction, and repeat this process 10,000 times. We find a slightly higher star-forming fraction for the Gold Sample of 87.8% (transitioning fraction of 2.4% and quiescent fraction 9.8%). However, we find a 68% confidence interval of ±5.7% from the random draws; thus, it is consistent with the full sample. This would imply that the most robust host associations are a fair representation of all included hosts. We find slightly different star-forming fractions compared to the full sample in the Silver Sample (∼73.6%) and more consistent fractions in the Bronze Sample (∼89%; Figure 1). To check whether the differences are statistically significant, we draw 19 galaxies (the number of Silver Sample hosts) from the full sample, and find a star-forming fraction of ≈84% ± 8% (68% confidence); the Silver star-forming fraction is just below the 68% confidence interval. This could be because Silver Sample hosts are more offset from their host galaxies (Fong et al. 2022), and larger physical offsets are noticed in non-star-forming hosts (see Section 5.2 for more details). When drawing nine galaxies (the number of Bronze Sample hosts) from the full sample, we find a median star-forming fraction of ≈88% ± 12% (68% confidence), consistent with the full, Gold, and Silver Samples. We find generally that the differences in the distributions are not statistically significant, demonstrating that the distribution in host galaxy type is not strongly dependent on the robustness of host association. However, the difference in the star-forming fraction of the Silver Sample may signify that quiescent and transitioning hosts have higher values of Pcc due to their faintness and/or the short GRB offset.
We also find that the photometric redshift sample contains a much higher percentage of star-forming galaxies (∼95%) than the full sample and the spectroscopic redshift sample. As there are 20 photometric redshifts, we do find that the star-forming fraction is out of the 68% confidence interval. As shown in Figure 3, the photometric sample covers a significantly higher redshift range, and the difference in the star-forming fraction is likely due to the prevalence of star-forming galaxies at z > 1. From the photometry in Fong et al. (2022), we find that many of these galaxies have little color variation, or difference in apparent magnitudes, between photometric bands, signaling no clear spectral breaks, typically a trait of actively star-forming galaxies. These faint galaxies, if truly at higher redshifts, are thus likely only visible due to the young, massive stars that dominate the SED light in the rest-frame UV (Conroy 2013).
Previous studies based on smaller samples of events used the "early-type" and "late-type" classifications, or morphological classifications such as elliptical and disk-dominated. Leibler & Berger (2010) found that ≈60% of their 19 short GRB hosts are late-type galaxies, and Fong et al. (2013) found that ≈61% of their 26 hosts are late-type (which includes the 19 bursts mentioned in the former work). We find that when we draw 19 and 26 "galaxies" from our population 10,000 times, the Leibler & Berger (2010) and Fong et al. (2013) star-forming fractions are both outside of the 68% confidence interval. The larger star-forming fractions found here are likely due to the much larger host sample, inclusion of a significant number of photometric redshift bursts at z > 1, and the more quantitative and uniform classification used in this work. Theoretical simulations of BNS mergers, such as that of Chu et al. (2022), which also uses the early- and late-type definitions, find similar star-forming host fractions of ≈81%–84%. Mandhai et al. (2022) predicted a lower star-forming fraction for BNS and NSBH hosts (≈60%–70%); however, they include simulated hosts that are too faint to be detected with ground-based imaging. It is still, however, likely that a majority of the nonhost associated GRBs described in Fong et al. (2022) are too faint because they exist at higher redshifts, where star-forming galaxies dominate. Overall, our work is consistent with previous studies that star-forming galaxies make up a majority of short GRB hosts.
4.3. Star-forming Main Sequence
We now determine how well short GRBs trace the normal field galaxy population by comparing their SFRs and stellar masses to the SFMS. The SFMS is a well-studied, redshift-dependent galaxy correlation that is observed to be followed by star-forming galaxies as they gain stellar mass, and is where galaxies spend a majority of their lifetime (Speagle et al. 2014; Whitaker et al. 2014; Leja et al. 2022). By comparing a host population to the SFMS, we can infer how the short GRB progenitor population traces a combination of stellar mass and star formation in galaxies. In the case of short GRBs, this can inform how galaxy-targeted searches following GW events could be performed (Section 6.2).
The SFMS is often parameterized by a single or broken power-law in sSFR–M* space (Speagle et al. 2014; Whitaker et al. 2014; Leja et al. 2022). Here, we use Equations (9) and 10 and the "ridge" values listed in Table 1 in Leja et al. (2022). This relation provides an excellent comparison to our data set as it is computed from the Prospector-derived SFRs and stellar masses of the COSMOS2015 (0.2 < z < 0.8; Laigle et al. 2016) and 3D–Hubble Space Telescope (HST; z > 0.5; Skelton et al. 2014) galaxy surveys. The COSMOS2015 and 3D-HST photometric surveys contain >67,000 galaxies, all with ≥17 bands of photometry, including the near-IR (NIR) wavelengths, which allows for more accurate determinations of these properties. The inferred stellar population properties from the galaxies in these surveys make an ideal comparison set as SFRs and stellar masses, and thus the SFMS, are determined across the range of redshifts relevant for our short GRB sample (0.2 < z < 3.0). Instead of using SFR versus M* to describe the SFMS, we use sSFR versus M*, as sSFR normalizes the amount of star formation per unit M*, which is useful when comparing galaxies over a wide range of stellar masses.
We show the comparison of the short GRB host population to the SFMS and COSMOS2015 and 3D-HST field galaxies in Figure 6, divided into five redshift bins (z < 0.5, 0.5 < z < 1.0, 1.0 < z < 1.5, 1.5 < z < 2.0, and 2.0 < z < 3.0). We find that across all redshifts, short GRB hosts tend to populate the entire SFMS, demonstrating that short GRBs are good tracers of star formation given their host stellar masses. We find that star-forming hosts are in the range of , quiescent hosts are in the range
, with transitioning hosts generally between the two populations (Table 3). We also note that Leja et al. (2022) used a nonparametric SFH to determine the SFMS. Were we to use the same SFH, the masses would likely increase by ∼0.2 dex and sSFR would decrease slightly, with more significant differences at higher redshifts (Leja et al. 2019). However, these differences would not be significant enough to change the trend noticed between short GRB hosts and the SFMS; however, they might explain why short GRB hosts at higher redshifts tend to lie above the SFMS.
Figure 6. 3D-HST and COSMOS2015 Prospector derived sSFRs and stellar masses (Brammer et al. 2012; Laigle et al. 2016; Momcheva et al. 2016; Leja et al. 2022) compared to the short GRB host Prospector-derived sSFRs and stellar masses over 0.0 < z < 3.0. The 3D-HST and COSMOS2015 data is plotted as a normalized histogram in gray. We use the definitions of star-forming (blue circles), transitioning (yellow squares), and quiescent (red triangles) galaxies from Tacchella et al. (2022) to plot the short GRB host data. Upper limits are shown as triangles, and hosts with photometric redshifts are shown as white-filled shapes. We also circle the host of GRB 170817 (NGC4993). We plot the broken-power-law SFMS equations derived in Leja et al. (2022) for each redshift bin. Across all redshifts, short GRB hosts track the SFMS well, implying they are good tracers of star formation given their stellar mass. The bottom-right plot shows the star-forming (SF) and quiescent/transitioning (Q/T) fractions of short GRBs in comparison to those of the 3D-HST and COSMOS2015 field galaxies. The short GRB host star-forming and non-star-forming fractions are comparable to the field galaxy's fraction at all redshifts, with a slight dominance of non-star-forming hosts at z < 0.25.
Download figure:
Standard image High-resolution imageFor comparison, we also include NGC4993, the quiescent host of GW170817/GRB 170817A (Abbott et al. 2017b; Blanchard et al. 2017; Levan et al. 2017; Pan et al. 2017; Palmese et al. 2017). While NGC4993 is not the most-massive host compared to the short GRB host sample at z < 0.5, it is by far the most quiescent. We find that compared to the short GRB host sample, GW170817 occurred in a unique environment, as the host has low star formation (sSFR ≈10−12 yr−1; ≈10−3 yr−1 less than the entire host population and ≈10−1 yr−1 less than the non-star-forming hosts) for its stellar mass (M* ≈ 4.1 × 1010 M⊙). The majority of galaxies at this stellar mass and redshift (≈75% in the 3D-HST sample) are expected to be star-forming. We also compare GRB the hosts of GRBs 060614 and 211211A to the host population, as they are low-redshift long GRBs with significant evidence for an NS merger origin (Della Valle et al. 2006; Gal-Yam et al. 2006; Rastinejad et al. 2022). Both hosts have some of the lowest masses and SFRs compared to the short GRB host population (bottom 12% and 19%, respectively). However, they are consistent with the short GRB host population. Their stark contrast to the stellar mass, sSFR, and optical luminosity (Fong et al. 2022) of NGC4993 highlights the range of environments of compact object mergers.
In the final panel of Figure 6, we show the fraction of quiescent and transitioning galaxies compared to the 3D-HST and COSMOS2015 data sets. We use the stellar mass and sSFR data in Leja et al. (2019) to determine the quiescent and transitioning fractions of the field galaxy population with the same method described in Section 4.2. We find that the short GRB host fractions are comparable with the field galaxy population at all redshifts, and begin to dominate the population at z < 0.25 (although numbers are small in this regime). As galaxies at z < 0.25 are easily observable, we suspect that this does represent the true rate of short GRBs in non-star-forming galaxies. Since quiescent hosts also have larger mass-weighted ages, this population could be indicative of a higher fraction of long-delay-time progenitors. The lack of quiescent and transitioning galaxies at high redshifts (z > 1.0) is likely due to observational limitations, as these galaxies are much fainter and thus much more difficult to detect with ground-based imaging.
Overall, we find short GRB hosts tend to have the expected stellar mass given their SFRs and thus tend to trace the SFMS. Furthermore, the fractions and stellar masses for star-forming, transitioning, and quiescent hosts are also within the observed range of field galaxies at their respective redshifts.
4.4. Stellar Mass Distribution
We next compare the short GRB host stellar masses to that of the normal field galaxy population to understand how they trace the stellar mass distribution in galaxies. We divide the host population into three redshift bins (z < 0.5, 0.5 < z < 1.0, and z > 1.0) and compute the stellar mass CDF from 5000 realizations of the Prospector-derived posterior distributions and the medians of those distributions. We use a mass-weighted Schechter function of the COSMOS2015 galaxy survey, detailed in Davidzon et al. (2017), to represent the stellar mass distribution of field galaxies within the three redshift bins. The Schechter function is an empirical parameterization for the number of galaxies at a given stellar mass, and the Davidzon et al. (2017) function ensures that the field galaxy stellar mass distribution is not biased by any observational limits of the COSMOS2015 survey (e.g., correcting for missing low-luminosity or high-redshift galaxies). The Schechter function in Davidzon et al. (2017) was built from stellar masses produced by a delayed-τ SFH from the SED modeling code GALAXEV (Simha et al. 2014). As choice of SFH has a greater effect on stellar population property values than other galaxy modeling assumptions, our Prospector-produced stellar masses are comparable to theirs due to common SFH. To test whether the occurrence rate is proportional to stellar mass, we convert the Schechter function into a total mass budget by multiplying it by a uniform distribution of stellar masses between 107 and 1013 M⊙ and calculate the CDF. We show the distributions in Figure 7. We find that the medians of the field galaxy stellar masses in our three redshift bins are higher by ≈0.3–0.5 dex than those of the host distributions.
Figure 7. CDFs of the short GRB host stellar masses (teal) compared to the mass-weighted Schechter function of the COSMOS survey's active galaxies (orange; Davidzon et al. 2017) across three redshift bins spanning 0 < z < 3.0. For the short GRB host mass distributions, we plot 5000 realizations on the CDF along with the median CDF (teal solid line). The Schechter function is used to represent the luminosity-corrected distribution of normal field galaxy stellar masses. Insets: we determine if the short GRB host and the mass-weighted Schechter functions are drawn from the same underlying distribution through an Anderson–Darling (AD) test, the distribution of PAD values for all realizations. The dashed lines in the insets show PAD = 0.05, with arrows in the direction of the point below which we reject the null hypothesis that they are drawn from the same underlying mass distribution. Across all redshift bins, we can reject the null hypothesis that short GRB hosts trace the stellar mass distribution in galaxies.
Download figure:
Standard image High-resolution imageFor a quantitative comparison to the short GRB host sample, we perform AD tests with the null hypothesis that both populations are derived from the same underlying stellar mass distribution. We perform an AD test for each of the 5000 realizations of the CDF and the Schechter function to build a distribution of probabilities (PAD), shown in the inset of Figure 7. For all redshift bins, all PAD < 0.05, and therefore, we can reject the null hypothesis that short GRB hosts trace the stellar mass distribution of field galaxies. Given that there are more high-mass galaxies in the normal field galaxy population than in the short GRB host sample, this also implies that high-mass galaxies are not necessarily more likely to host short GRB progenitors. Furthermore, given that short GRB hosts are tracing the SFMS, this indicates that short GRBs do not trace stellar mass alone; instead they trace a combination of star formation and stellar mass.
We note that if the true short GRB host population traces the stellar mass distribution in galaxies, and the discrepancy is purely an observational selection effect, this would require missing a fraction of high-mass hosts. As an exercise, we add a number of additional galaxies to the short GRB sample to determine how many more high-mass hosts are needed for the observed host population to trace the stellar mass distribution. For z < 0.5, we find that we would need to nearly double the number of high-mass galaxies (four to six more 1010.5
M⊙; >50% of PAD > 0.05). However, since high-mass galaxies at z < 0.5 are easily detectable given their higher luminosities, it is very unlikely that we are missing this many hosts in our sample. If anything, there is an observational bias against identifying low-mass galaxies as short GRB hosts, which would further drive the distributions apart and strengthen our conclusions.
4.5. Stellar Population Age
The stellar population age is an important property in understanding short GRB progenitors, as it can be used as a loose proxy of the formation and merger timescales, or the delay time, of the progenitor (e.g., Safarzadeh et al. 2019). Indeed, past work has shown that the wide range of stellar population ages is consistent with the range of compact object merger timescales (e.g., Shin & Berger 2007; Zheng & Ramirez-Ruiz 2007; O'Shaughnessy et al. 2008; Leibler & Berger 2010; Berger 2014). In Figure 5, we show the mass-weighted age (tm
) median CDF and 5000 realizations of the Prospector-derived posteriors for all, star-forming, and transitioning and quiescent hosts. As expected, we find a large difference in ages between the star-forming and non-star-forming populations; star-forming hosts have Gyr (median and 68% credible interval), whereas non-star-forming hosts have
Gyr. As there are more star-forming hosts in the sample, the entire population is weighted more heavily by the younger hosts and thus has a median age as a population of
Gyr. We find that the shape of the posterior distributions of the entire host population roughly follows that of a power-law decline, which is the expected shape of the delay-time distribution of BNS mergers (Peters 1964; Nakar et al. 2006; Lee & Jeong 2010; Hao & Yuan 2013). Using the age distribution as a proxy for progenitor delay times, this is consistent with previous findings that there exists a broad range of delay times for the expected progenitors of short GRBs, including BNS and NSBH mergers (Belczynski et al. 2002; Guetta & Piran 2005; Nakar et al. 2006; Hao & Yuan 2013; Wanderman & Piran 2015; Anand et al. 2018).
To understand age evolution with redshift, we show the stellar population age normalized by age of the universe at each redshift (tH
(z)) versus redshift in Figure 8. We find that at z > 1, the population of hosts is dominated by stellar populations that are only ≈20% the age of the universe. We also find that all but three hosts with photometric redshift determinations have ages <2 Gyr, are heavily clustered at ≲1 Gyr (Figure 9), and are more skewed toward higher redshifts. This is consistent with the association of short GRB progenitors at z > 1 to more recent star formation, as the youngest stellar population ages are just slightly above the stellar evolutionary timescale of NS progenitors Myr). At face value, higher-redshift host galaxies are closer in time to the epoch of peak star formation of the universe (Madau & Dickinson 2014), and thus, the dominance of younger stellar population ages at z > 1 is likely indicative of the fast-merging binary population comprising mostly BNS progenitors with tight orbital separations at formation and short delay times (Belczynski et al. 2002; Beniamini & Piran 2019; Safarzadeh et al. 2019). This also implies a larger population of fast-merging, z > 1 BNS systems than previously confirmed (e.g., Selsing et al. 2018; Paterson et al. 2020; O'Connor et al. 2022), as we have nearly 15 hosts at z > 1.0. Taking into account completeness, we believe 25%–44% of Swift short GRBs are from z > 1 (Fong et al. 2022). In contrast, at z < 1, the ages span a wide range compared to tH
(z). In particular, the growing number of older stellar populations with tm
/tH
(z) ≈ 0.8 indicates that these represent hosts of binary systems with longer delay times that formed during the epoch of peak star formation but did not merge until much more recently (Figure 8). We note that the longer delay times are most probably derived from the intrinsic properties of isolated binary systems, including orbital separation, rather than delayed mergers from dynamical formation in, e.g., globular clusters, as it is highly unlikely for these BNS systems to merge within a Hubble time (Bloom et al. 2002; Ye et al. 2020). We also find a much wider range of tm
/tH
values at z < 1 than at z > 1, consistent with the fact that binary production has decreased and we are no longer dominated by very fast-merging channels at z < 1. Taken together, we find this likely implies that BNS/NSBH production has decreased since z ≈ 1, as the population in the near universe is no longer dominated by fast-merging systems as at higher redshifts. We use these properties to quantify the DTD parameters in Zevin et al. (2022). A similar DTD is observed with the Galactic BNS population in which there are many more binaries expected to merge quickly (<1 Gyr; Beniamini & Piran 2019).
Figure 8. The mass-weighted ages (tm ) of short GRB hosts normalized by the age of the universe at their redshift (tH (z)) vs. redshifts, color-coded by their host type (blue is star-forming, yellow is transitioning, and red is quiescent). White-filled symbols are photometric redshifts and single color symbols are spectroscopic redshifts. We find that closer to the epoch of star formation (z = 2), there are more hosts that are much younger than tH (z). As redshift decreases, we see that the population of hosts begins to be dominated by those with ages closer to tH (z). This would imply that production of the short GRB progenitors is decreasing with redshift, as the short GRBs found at lower redshifts likely formed closer to the epoch of star formation and have a longer delay time.
Download figure:
Standard image High-resolution imageFigure 9. Stellar masses and mass-weighted ages of short GRB hosts, with the same color and shape schematic as in Figure 8. Hosts with photometric redshifts are represented with white-filled shapes. We mark NGC4993 and the hosts of GRBs 060614 and 211211A for reference. We find that quiescent host galaxies make up the oldest and most-massive hosts in the sample, and that the hosts of short GRBs span nearly the full range of galaxy properties.
Download figure:
Standard image High-resolution imageWe also explore the ages in the context of host stellar mass, which we show in Figure 9. As expected, we find that as the host age increases, the host mass also increases and that quiescent hosts comprise the oldest and most-massive galaxies in the sample. This is in alignment with standard galaxy evolution, as galaxies gain in mass as they evolve, through major and minor mergers (Cattaneo et al. 2011; Whitaker et al. 2014). It is unclear if galaxy mergers have an effect on binary production, although NGC4993 (GRB 170817A) did have a recent galaxy merger (Palmese et al. 2017; Ebrová et al. 2020; Kilpatrick et al. 2022). Thus, the wide range of short GRB host masses distribution reinforces the wide distribution of delay times possible for their progenitors. Overall, short GRBs originate in galaxies spanning the full range of properties. Although the parametric SFH assumption is known to produce younger ages than nonparametric SFH models by a factor of ≈3–5 (Papovich et al. 2001; Carnall et al. 2019), this will not affect the relative fraction of young and old stellar populations.
Overall, we find that the mass-weighted ages of short GRB hosts reinforce a wide range of progenitor delay times, with more support for fast-merging systems at z > 1, and a decrease in binary production toward lower redshifts (see Zevin et al. 2022 for more). Ground-based GW detectors will not be able to reach these cosmological distances for ≳15 yr. Therefore, they will not be able to fully sample the DTD of BNS/NSBH mergers, such as is possible with short GRBs.
4.6. Mass and Stellar Metallicity
Finally, we explore the stellar mass–metallicity relation of short GRB hosts. Previous studies on the effects of binary compact object merger formation and metallicity have found that lower-metallicity environments during star formation lead to more compact binaries (Belczynski et al. 2010; O'Shaughnessy et al. 2017). Stellar metallicity has been shown to trace both the star formation and age of galaxies, with more evolved stellar populations being more metal-rich, thus giving insight into the types of stars in an environment (Gallazzi et al. 2005; Panter et al. 2008). Studies have also shown that stellar metallicity increases with increasing stellar mass, as more-massive galaxies contain older stellar populations and thus more metals.
In Figure 10, we show the derived log(M*/M⊙) and log(Z*/Z⊙) in comparison to the Gallazzi et al. (2005) mass–metallicity (MZ) relation, as well as the COSMOS2015 and 3D-HST field galaxies color-coded by their star formation classification. We find that most short GRB hosts fall within the expected metallicity range given their stellar mass and star formation classification. Figure 10 also shows that BNS/NSBH mergers are able to form over a wide range of metallicity environments, even in the local universe; this is highlighted with the low-redshift hosts of GRBs 170817A (NGC4993) and 211211A. This suggests that short GRB formation efficiency across all redshifts is not dependent on the amount of metals in their environment. If short GRBs only occurred in low-metallicity environments, similar to long GRBs and SLSNe, we would witness the entire population very connected to recent star formation. Conversely, a dependence on higher-metallicity hosts would suggest that short GRBs overwhelmingly come from less star-forming populations. The spread in host metallicities therefore shows that short GRBs are derived from both fast-merging populations linked to recent star formations and slowly merging progenitors not linked to star formation. Furthermore, it shows that the short GRB progenitor formation efficiency is also independent of metallicity, which is consistent with predictions for BNS systems (Chruslinska et al. 2018; Giacobbo & Mapelli 2018; Neijssel et al. 2019). However, NSBH systems have a higher formation efficiency in lower-metallicity environments (Klencki et al. 2018; Broekgaarden et al. 2022); thus, the MZ relation of short GRB hosts does not follow expectations if all events are derived from NSBH mergers.
Figure 10. The Gallazzi et al. (2005) stellar mass–metallicity relation (MZ; solid black lines represent the median and dashed black lines represent the 68% credible region) overplotted with short GRB hosts. The colors of both the 3D-HST and COSMOS field galaxies (small circles) and the short GRB hosts (large, outlined circles) denote their star formation classification, with blue representing star-forming, yellow representing transitioning, and pink representing quiescent galaxies. The open circles represent photometric redshift hosts. Short GRB hosts follow the MZ relation, fall within the expected range for their star formation classification, and span a wide range of metallicities. We also highlight NGC4993 and the host of GRB 211211A for reference, which occur in very different places in the MZ-plane. Overall, GRB-producing NS mergers form independently of the metallicity of their stellar populations.
Download figure:
Standard image High-resolution image5. GRB and Locations versus Host Properties
We next explore how GRB properties, such as T90, prompt fluence (15–150 keV), presence of detectable extended emission, optical afterglow detection, afterglow luminosities, and galactocentric offsets, compare with their host properties to understand if host properties influence the GRB properties.
5.1. Duration, Fluence, and Extended Emission
We first compare host type (star formation classification) to T90 and fluence in Figure 11. We find values for T90 and fluences in the Swift BAT3 catalog 19 (Lien et al. 2016). While short GRBs are typically defined with T90 ≲ 2 s, nominally at the T90 ≈ 2 s boundary between short and long GRBs (which are derived from young, massive stellar progenitors and are known to only occur in star-forming galaxies; Leibler & Berger 2010; Svensson et al. 2010; Perley et al. 2013; Wang & Dai 2014; Vergani et al. 2015; Niino et al. 2017), we might expect some contamination in our sample. In Figure 11, we find that for very short durations, T90 ≲ 0.1 s, and low fluences, ≤3.1 × 10−8 erg cm−1, almost all hosts are quiescent, while almost none are star-forming in this region. Since γ-ray fluence has been shown to correlate with X-ray and optical afterglow luminosity (Gehrels et al. 2008), this could be a product of the lower ISM densities of quiescent hosts. However, including GRB 170817A and the extended emission short GRB 050724A, quiescent hosts span the full range of parameter space. The interval of bursts with T90 ≈ 0.1–2 s is fully occupied by star-forming and transitioning hosts, except for the quiescent hosts of GRBs 100625A and 100117, (Figure 11) but span 2 orders of magnitude in fluence. Thus, while we do not find any compelling evidence for significant contamination from massive star-originated events (see also Fong et al. 2022), we do note the prevalence of quiescent hosts with lower fluences and shorter durations.
Figure 11. The T90 and fluences of short GRBs, marked by the star formation classification of their hosts. We see that all types of hosts (star-forming, transitioning, and quiescent) span the full ranges, implying that host type does not influence these GRB properties. Short GRBs in quiescent hosts, however, do tend to have the lowest fluences. We mark GRBs 060614, 170817, and 211211A for reference.
Download figure:
Standard image High-resolution imageWe further investigate how the hosts of extended emission bursts compare to those of the normal short GRB host population. Several events in our sample have extended emission detectable in the γ-rays, which results in T90 ≫ 2 s (Norris & Bonnell 2006; Lien et al. 2016; Fong et al. 2022). 20 While the physical origin of extended emission is uncertain, it has been suggested that they may be derived from NSBH merger systems, on the basis of their potentially smaller galactocentric offsets, and larger γ-ray energies (Troja et al. 2008; Gompertz et al. 2020). By default, extended emission bursts share similar γ-ray properties to long GRBs given their long durations and larger fluences, but unlike long GRBs, they are not accompanied by Type Ic SNe. We find that the normal short GRB host population and distribution of extended emission bursts are very similar in stellar mass and age. For the stellar masses and mass-weighted ages, we find pAD > 0.25 between the two populations, and thus cannot negate that they are derived from the same underlying mass and age distributions. We also do not find any difference between the star formation classifications between each population: we randomly draw 12 hosts (the number of extended emission hosts) from the normal short GRB host population and find a star-forming fraction of 83.3% ± 10.5 (1σ). This is consistent within the 68% credible interval to the observed star-forming fraction of extended emission hosts of 91.6%. Thus, consistent with the conclusions of Fong et al. (2022), the host properties alone do not clearly support a distinct progenitor for extended emission bursts.
5.2. Galactocentric Offsets
We next explore the short GRB offsets compared to star formation classification, host stellar mass, and stellar age. Nugent et al. (2020) found that the short GRBs in galaxy clusters (GRBs 050509B, 090515, and 161104A), which host the oldest and most-massive galaxies in the universe, were significantly more offset from their hosts than the rest of the population, suggesting that the host type or age might affect the progenitor's offset (see, also, Berger et al. 2007; Shin & Berger 2007). We compare the projected physical (in kiloparsecs) and host-normalized offsets (in effective radii units re
) from Fong et al. (2022) of the star-forming and non-star-forming (quiescent and transitioning) hosted short GRBs in Figure 12. We examine both the observed offset CDF and offset distribution including uncertainties, which is built from 5000 realizations of a Rice distribution (Equation (2) in Blanchard et al. 2016; see Section 6.2 of Fong et al. 2022) using the offsets and their 1σ uncertainties in Fong et al. (2022). We find that the short GRBs in star-forming hosts have median offsets of kpc and
(sampled distributions give
kpc and
). Meanwhile, short GRBs in non-star-forming hosts have larger observed offsets of
kpc and
(sampled distributions give
kpc and
; Figure 12).
Figure 12. CDFs of observed physical (left) and host-normalized (right) offsets for star-forming (dark blue) and quiescent/transitioning (red) hosts. The CDFs produced from the Rice distributions are shown in the respective lighter colors for each host star formation classification. We show medians with the arrows and highlight the offset of GRB 170817 with a dashed black line. We find significant differences in physical offsets between the star-forming and non-star-forming hosts, but more similar distributions when looking at the host-normalized offsets. Taken together, this likely implies that short GRB progenitors in non-star-forming galaxies have migrated farther than those in star-forming hosts.
Download figure:
Standard image High-resolution imageOverall, the short GRBs in non-star-forming hosts are observed to be farther from their hosts than short GRBs in star-forming hosts (Figure 12; pAD = 0.004). However, we find that the differences are less clear for host-normalized offset distributions (pAD = 0.25). The observed offsets result from a combination of the systemic velocity of the binary progenitor and its delay time (Fryer et al. 1999). Given similar systemic velocities of binary progenitors in both star-forming and non-star-forming hosts, the observed physical offset differences can be naturally explained by the longer delay times in non-star-forming hosts (commensurate with their older ages; See, Figure 5). Galaxies can also grow significantly over approximately gigayear timescales, and thus, short GRBs from non-star-forming hosts may have gotten kicked when the galactic potential was much different than it is today. However, even if all short GRB progenitors migrated an approximately equal physical distance from their host centers, this would translate to less distance traversed within the quiescent or transitioning hosts, as their radii are larger. Rather, the differences in host-normalized offsets between star-forming and non-star-forming hosts are not expected to be as prominent as in physical space. Finally, we note that we only have host-normalized offsets for the most secure host associations, thus creating a slight selection effect toward bursts at smaller offsets.
5.3. Afterglow Properties
We also examine if host properties are correlated with short GRB optical afterglow detectability and luminosities. Naively, we expect that star-forming hosts will have higher ISM densities on average and, thus, brighter optical afterglows (see Granot & Sari 2002) than those in quiescent or transitioning hosts. However, we see no difference in distributions between optically detected and nondetected GRBs, finding PAD > 0.24 for all tested properties.
We further compare the optical afterglows of short GRBs to their host types. We collect optical afterglow luminosities or upper limits from several minutes to several months after the burst (Fong et al. 2015; Rastinejad et al. 2021 and references therein), and convert to luminosities in the host frame using the redshifts found in this work. Here, optical refers to afterglow detections taken in optical and NIR photometric filters, which de-redshifted lie in the optical wavelengths, approximately at the same rest-frame wavelength. In total, we have 49 optical afterglow limits and detections, 82.4% coming from short GRBs in star-forming galaxies, 7.8% in transitioning, and 9.8% in quiescent, which roughly matches the breakdown of our full sample. We compare all afterglows at a common rest-frame time of tRF = 5 hr, as the afterglow has not significantly faded., interpolating or extrapolating from their nearest detection or deepest upper limit when necessary. We build Gaussian distributions of the flux of the detection using their median and 3σ uncertainties. For the 3σ upper limits, we build a top-hat function from a lower limit of 1038 erg s−1 (representing ≳100 times lower than the faintest detected afterglow) to the value of the upper limit. We draw randomly from these distributions 1000 times and build CDFs of the luminosities based on whether the galaxy is star-forming or not. We show the median and 1σ confidence interval of these distributions in Figure 13. While it appears that short GRBs in quiescent and transitioning hosts have fainter afterglow luminosities than those in star-forming hosts, we cannot rule out the null hypothesis that they are drawn from the same underlying distribution (). Thus, we do not find statistical support for fainter optical afterglows in non-star-forming host galaxies.
Figure 13. The optical afterglow luminosities of short GRBs at a common rest-frame time of 5 hr in the host frame. While the optical afterglows appear to be fainter in quiescent and transitioning hosts (blue), we find no statistically significant difference when compared to those in star-forming (blue) hosts, implying that luminosity and detectability are not strongly dependent on global host properties.
Download figure:
Standard image High-resolution imageFinally, we explore the circumburst densities of short GRBs along with their host properties. We find circumburst densities for 29 short GRB hosts (≈80% occur in star-forming hosts) in Fong et al. (2015) and Schroeder et al. (2020) and compare them to the hosts' properties. We find no correlation between any host property and the circumburst density, and find that short GRBs in star-forming, transitioning, and quiescent galaxies span the full range of circumburst densities. We also do not see any trend between circumburst densities and the short GRB offsets. As these densities only probe the local environment around the GRB, we infer from this that short GRBs in star-forming and quiescent galaxies both can occur in low-density environments. It may also indicate that we do not have the sensitivity to detect any difference between short GRB densities in either host group.
6. Comparisons to Transient Environments and Implications for Gravitational Wave Follow-up
6.1. Supernovae and Long GRBs
To provide context for the short GRB host population and the conditions needed to form their progenitors, we now compare their host stellar population properties to those of other well-studied transients. We use the hosts of Type Ia SNe, Ca-rich transients, Type II SNe, SLSNe, and long GRBs, as these transients originate from a wide variety of progenitor models. We use the stellar population properties in Pan et al. (2014; Type Ia SNe at z < 0.09), Wiseman et al. (2021; Type Ia SNe at 0.2 < z < 0.6), Dong et al. (2022; Ca-rich transients), Svensson et al. (2010), Perley et al. (2013), Vergani et al. (2015), Wang & Dai (2014), Niino et al. (2017; long GRBs), and Schulze et al. (2021; SLSNe and Type II SNe). We caution that these studies use a variety of methods for determining stellar population properties, including applying different SED modeling codes, stellar mass and SFR relations based on colors, and elliptical and spiral galaxy stellar population templates, and thus, the very precise comparisons here are not possible. We choose stellar mass and sSFR as a basis for comparison, as they are the most cataloged host stellar population properties for the other transients, and stellar mass is generally fairly impervious to modeling differences. We choose to only compare the stellar masses and sSFRs at z < 0.5 for all transients, as this is where there exists the most overlap among short GRBs and the comparison samples. This is also the redshift range where a majority of short GRBs lie; thus, we have a large sample to compare to the other transient hosts.
We plot the CDFs of the other transient hosts from their inferred log(M*/M⊙) and sSFRs, and compare to short GRBs in Figure 14. We initially find that long GRBs and SLSNe, which both are derived from young massive star progenitors, occur in hosts with smaller stellar masses and higher rates of star formation than those of short GRBs. Type Ia SNe, which have older stellar progenitors, appear to share similar host properties to those of short GRBs. We also find that the stellar mass distributions of Ca-rich transients are very similar to those of short GRBs, although only a portion of Ca-rich transients may truly be derived from older stellar progenitors, while the rest are likely from massive stars. To further quantify this, we perform AD tests and can reject the null hypothesis that short GRB hosts are derived from the same stellar mass and sSFR distributions as Type II, long GRB, and SLSNe host galaxies. In particular, this reinforces the distinct nature of short GRB environments from those that bear massive star explosions. In contrast, we cannot reject the null hypothesis that short GRB hosts' sSFRs are derived from the same underlying distributions as those of Type Ia and Ca-rich transients, as well as the stellar masses between short GRBs and Ca-rich transients (Type Ia SNe hosts appear to have slightly higher stellar masses). Since we know that Type Ia SNe originate from evolved single- or double-degenerate progenitors and have an observed power-law DTD (p(t) ∝ t−1; Peters 1964; Belczynski et al. 2002; Totani et al. 2008; Maoz et al. 2012; Graur et al. 2014; Beniamini & Piran 2019), it comes as no surprise that their global host environments are similar to those of short GRBs. While the progenitors of Ca-rich transients are uncertain, with models ranging from white dwarf origins in globular clusters (Shen et al. 2019) to core collapse from a stripped-envelope low-mass massive star (De et al. 2021), it is clear that their progenitors do not select starkly different environments from short GRBs. We also note that the Ca-rich sample only extends to z ∼ 0.02, over which galaxies in general have lower inferred sSFRs and higher stellar masses; thus, our observations could also be a product of selection effects of the current Ca-rich host sample.
Figure 14. CDFs of log(M*/M⊙) (left) and sSFR (right) for short GRB hosts (purple), Type Ia SNe hosts (green), Type II SNe hosts (orange), Ca-rich transients hosts (red), SLSNe hosts (pink), and long GRB hosts (blue). For the short GRB host CDFs, we draw 5000 random samples from the respective Prospector posterior distribution, and plot the medians of the resulting CDFs as well as the 68% credible region to represent uncertainty (in light purple). We perform AD tests between all 5000 short GRB host distributions and the other transient hosts' distributions and show the percentage of tests that reject the null hypothesis that they are derived from the same underlying distribution in the legend. We find Type Ia, Ca-rich transient, and short GRB hosts' sSFRs to be less than those of other transients and find their stellar masses to be higher.
Download figure:
Standard image High-resolution image6.2. Comparison to GW170817 and Implications for Gravitational Wave Follow-Up
The joint detection of GW170817 with the short GRB 170817A provides direct evidence on the origins of some short GRBs. Thus, the host environments of short GRBs may provide information on the environmental properties of local, GW-detected mergers. In this work, we have compared the environment of GW170817, the only GW-detected BNS merger, with a host galaxy (NGC4993), and we summarize our findings here. With an age of ≈10.41 Gyr and an sSFR ≈−12.34 yr−1, NGC4993 is older and has a lower sSFR than all galaxies in the sample. It is fairly massive, at the ∼83.5th percentile for stellar mass although not the most-massive host. The fact that NGC4993 is unique galaxy to host NS mergers is also supported by simulations of BNS merger environments. Indeed, Mandhai et al. (2022) found NGC4993 to have one of the lowest SFRs compared to its sample of simulated BNS and NSBH mergers at z < 0.4, and instead found much higher probabilities that star-forming galaxies will host these mergers. While NGC4993 might be unusual compared to the hosts at higher redshifts, it may not be that extraordinary of a host given galaxy evolution, as the low-redshift galaxy population has a more significant fraction of quiescent galaxies (see Figure 6). Thus, as GW detectors are still constrained to local universe BNS and NSBH events (z < 0.04 for BNS mergers and <0.08 for NSBH mergers in LIGO/Virgo O4; Abbott et al. 2018) and given that short GRB hosts trace the SFMS, we might expect the majority of NS merger hosts will be star-forming, but with a more prominent fraction of quiescent hosts than at higher redshifts. If all local mergers are derived from NGC4993-like galaxies, this would indicate that the population may solely be tracing stellar mass, in which case we could infer that local mergers exclusively represent a long-delay-time population. On the other hand, GRB 211211A has a host with properties that contrast greatly from those of NGC4993. Thus, if this event is truly derived from an NS merger, this would indicate that mergers in the local universe also trace a combination of stellar mass and star formation, similar to the cosmological merger population.
One of the major challenges of identifying EM counterparts to GW events is the large localization regions that need to be searched. Many optical efforts have turned to galaxy-targeted searches, which search around potential host galaxies for counterparts. At present, such follow-up searches generally rank galaxies by either their B-band luminosities, which is correlated with star formation, or less frequently their K-band luminosity, which traces the stellar mass (White et al. 2011; Blanchard et al. 2017; Palmese et al. 2017; Dalya et al. 2018; Ebrová et al. 2020; Kilpatrick et al. 2022). The results from our catalog, if representative of the local universe population, suggest that ranked-based searches should not singularly depend on either the B- or K-band luminosities. Our results suggest that follow-up searches that invoke these observables should involve a combination of these luminosities, as the most likely host for a BNS event will not always be the most-massive nor the most star-forming galaxy. Our study also demonstrates that short GRBs span a diverse range of environments, and thus, ranked-based searches will not always be successful in identifying the most likely hosts.
7. Sample Biases
Finally, we discuss potential biases in the full sample and subsamples of short GRB host galaxies. Fong et al. (2022) discussed selection effects and the use of the probability of chance coincidence (Pcc) method to associate short GRBs to their hosts. There, we used inclusive criteria for host identification and found more high-redshift (e.g., z > 1) and low-luminosity hosts than previously reported. In this paper, we further find that there are almost no hosts with stellar masses of ≲108 M⊙ in our sample (except for GRB 060614 210726A's hosts), consistent with the lack of low-luminosity hosts despite our survey being sensitive enough (Fong et al. 2022). This may imply that short GRBs do not often occur in low-mass or dwarf galaxies because they produce less of the total stars and BNS systems in the universe, or that it is difficult to associate short GRBs to such hosts below the mass completeness limit. As a test, we create a Prospector model SED at z = 0.6 (the median redshift of host sample), MF = 108 M⊙, and the medians of the rest of the stellar population properties. We find that the host would have an i-band magnitude of ≈26 mag AB, at the limit of ground-based telescopes. If a progenitor experienced a significant kick from such a galaxy, and therefore exhibits even a moderate offset, the Pcc method would not generally favor such a host.
We also note that there are 21 short GRBs in the Fong et al. (2022) sample that either do not have confirmed hosts (six short GRBs) or do not have enough data for a Prospector fit (17 short GRBs), despite the fact that follow-up observations for these bursts were not limited by any significant observing constraint. Many of the latter hosts are simply too faint to be detected in multiple photometric bands with a single HST detection, and therefore have low luminosities and/or high redshifts (and likely low masses). If we assume an extreme case that all of these hosts exist at z ≥ 1.0, the redshift median of the short GRB sample would increase to z ≈ 0.9. Further, if all were low mass (≤108.5 M⊙), the mass median would decrease to log(M/M⊙) ≈ 9.4. Finally, if all of these hosts were quiescent galaxies, the star-forming fraction would decrease to ≈64%, which is more consistent with the results from Mandhai et al. (2022). However, given that we have evidence for disk morphologies for several of these hosts (see Fong & Berger 2013), it is unlikely that all of them are quiescent.
We next explore how results may change if we only analyze the Gold Sample hosts. As shown in Figure 4, there are no significant differences in stellar population properties between the Gold, Silver, and Bronze Samples, except potentially in redshift. Figure 8 in Fong et al. (2022) also shows that Pcc increases (poorer association) with increasing redshift, thus removing the Silver and Bronze Samples from any host analysis might bias our results slightly in favor of more low-redshift, brighter hosts. Since there are several redshift-dependent parameters explored in this paper, including stellar mass, SFR, and age, only examining the Gold Sample might decrease the significance of fast-merging short GRB progenitors occurring in young, star-forming galaxies.
A more significant bias in understanding short GRB hosts occurs if we only explore the spectroscopic redshift as opposed to the full sample. We find a significant difference in the redshift ranges (no significant changes between other stellar population properties) between the spectroscopic and photometric redshift samples, as the photometric redshift sample contains many z ≳ 1.0 high-redshift hosts. By not including the photometric redshift sample, we are undoubtedly excluding short GRB progenitors that formed and merged in the early universe, which has significant implications for their DTD. As short GRBs contribute to the r-process element nucleosynthesis, this may also affect our understanding of chemical enrichment of the universe. Therefore, in the absence of spectroscopic confirmation with, e.g., the James Webb Space Telescope, the inclusion of the photometric redshift sample will play a crucial role in understanding the environments and progenitors of short GRBs.
8. Conclusions
In this paper, we modeled 69 short GRB host galaxies with Prospector to determine their redshifts and stellar population properties, including stellar mass, mass-weighted age, (specific) SFR, metallicity, and dust attenuation. We used the host galaxies described in Fong et al. (2022), which are divided into three subsamples based on host association confidence: Gold (best associated), Silver, and Bronze Samples. For this work, we only model host galaxies when there are at least three host detections in different photometric filters. As the largest catalog of short GRB host properties to date, this sample has tremendous implications for their progenitor systems. We summarize our main conclusions below:
- 1.Short GRBs span a spectroscopically confirmed redshift range of 0.1 ≤ z ≤ 2.211, with a median of z ≈ 0.47. When including the short GRBs with photometric redshifts, which fills in the "redshift desert" at z ≳ 1.0, the median increases to z ≈ 0.64. The inclusion of the photometric redshift sample is important in quantifying the number of high-redshift, fast-merging progenitors. The Silver and Bronze Sample hosts span higher redshifts than the Gold Sample hosts.
- 2.The full short GRB host sample has
, SFR
M⊙ yr−1,
Gyr,
, and
mag. The Gold, Silver, and Bronze Samples all have stellar population properties that are consistent with the full sample.
- 3.Approximately 84% of short GRB hosts are classified as star-forming galaxies. Star-forming hosts dominate the host population at z > 0.25, and the star-forming fraction roughly follows that of the observed field galaxy population at these redshifts. At z < 0.25, we find a growing presence of quiescent and transitioning hosts. The frequency of short GRB hosts across redshift is fully consistent with the field galaxy population.
- 4.We find that across all redshifts, short GRBs have SFRs consistent with their stellar masses and populate the range of the SFMS. However, they do not trace the stellar mass distribution of field galaxies alone. This implies that short GRBs trace a combination of recent star formation and stellar mass, but their rate is not singularly dependent on either property.
- 5.Stellar population ages of short GRBs span a wide range (≈0.1–9 Gyr), implying a spread in progenitor delay times. The majority of ages are ≤1 Gyr. At z > 1, their ages are consistently only ≲20% of a Hubble time at their respective redshifts, indicating the dominance of fast-merging systems at these epochs. As redshift decreases, the population of hosts expands to include those with ages closer to the Hubble times at those redshifts. This indicates that systems at low redshifts likely formed closer to the epoch of peak star formation in the universe and have long delay times. This also suggests a decreasing NS binary formation efficiency from high to low redshifts.
- 6.Short GRB progenitors do not have a strong preference for stellar metallicity, but instead their hosts have a spread in metallicities that matches the field galaxy population. As the amount of metal content is inversely correlated to the amount of star formation, this shows that there are some short GRB progenitors directly linked to recent star formation, as well as progenitors very removed from it.
- 7.We find that most short GRBs in quiescent hosts appear to have shorter T90 values and lower fluences than those in star-forming hosts. We also see evidence for larger migrations overall from short GRBs in non-star-forming hosts than those in star-forming hosts. We find no correlation between host properties and short GRB optical afterglow detection, optical afterglow luminosity, and circumburst density. We also find that the population of bursts with extended emission is drawn from the same underlying distribution in terms of all tested properties, indicative of the same or similar progenitors for both populations.
- 8.Short GRB hosts have lower sSFRs and higher stellar masses than those of SLSNe, and long GRBs and Type II SNe, while they have similar such properties to the hosts of Type Ia SNe and Ca-rich transient hosts.
- 9.Compared to the full sample of short GRB hosts, NGC4993 (host of GW170817/GRB 170817A) is older than any other host and has a lower sSFR. It is fairly massive, at the 84% level for stellar mass. While quiescent, massive galaxies such as NGC4993 are more prevalent in the local horizon probed by GW detectors, the short GRB host population exhibits broader diversity even at low redshifts. Indeed, if the short GRB host population is representative of GW-detected merger hosts, then a combination of B- and K-band luminosities should be used in targeting potential host galaxies.
Overall, our work demonstrates that short GRB progenitors can form in a variety of environments, across a wide range of redshifts and SFHs. We find evidence for both a short delay-time population, which appears in young and actively star-forming stellar populations and are the majority of observed events, and a long-delay-time population, which becomes more prevalent at lower redshifts in older and quiescent galaxies. Our work paves the way for quantitative constraints on the true DTD by convolving the Prospector derived SFHs with the stellar population ages and redshifts (see Zevin et al. 2022).
Through modeling with Prospector, we have significantly augmented the known ranges of redshifts and stellar population properties probed by short GRBs. As NS mergers that lead to short GRBs are known to produce some (if not most or all) of the r-process elements in the universe, understanding their true redshift distribution allows us to better comprehend when in the universe's history these elements were created and the contribution from mergers. Furthermore, we can use this host stellar population catalog to better associate GW-detected BNS and NSBH mergers to their host galaxies, especially as the horizon of ground-based GW detectors increases. In light of the fact that short GRBs trace a combination of stellar mass and recent star formation, future galaxy-targeted searches for EM counterparts could be optimized to account for this information. At present, short GRBs are detected more frequently and out to farther cosmological distances than GW-detected mergers. Thus, this catalog provides a legacy sample for comparison that will not be matched by GW detections for many years to come. Furthermore, this catalog of host properties holds the promise to quantify specific NS binary traits, such as their natal kicks, orbital separations, and delay times.
We have provided a comprehensive catalog of short GRB host properties based on nearly two decades of discoveries. However, there remain bursts with no confirmed or very faint hosts that need characterization. In addition, deriving the true spectroscopic redshifts for the photometric sample will solidify the redshift distribution especially at z > 1. The James Webb Space Telescope will be pivotal in both endeavors, providing unprecedented high-quality spectra for inference on their stellar populations.
We thank the anonymous referee for valuable insight and suggestions. A.E.N. acknowledges support from the Henry Luce Foundation through a Graduate Fellowship in Physics and Astronomy. The Fong Group at Northwestern acknowledges support by the National Science Foundation under grant Nos. AST-1814782, AST-1909358 and CAREER grant No. AST-2047919. W.F. gratefully acknowledges support by the David and Lucile Packard Foundation, the Alfred P. Sloan Foundation, and the Research Corporation for Science Advancement through Cottrell Scholar Award 28284. Support for M.Z. was provided by NASA through the NASA Hubble Fellowship grant HST-HF2-51474.001-A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Incorporated, under NASA contract NAS5-26555. Y.D. is supported by the National Science Foundation Graduate Research Fellowship under grant No. DGE-1842165. A.L. is supported by the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation program (grant agreement No. 725246).
This work is based on observations taken by the 3D-HST Treasury Program (GO 12177 and 12328) with the NASA/ESA HST, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555.
This research has made use of the HST-COSMOS database, operated at CeSAM/LAM, Marseille, France.
This research was supported in part through the computational resources and staff contributions provided for the Quest high performance computing facility at Northwestern University, which is jointly supported by the Office of the Provost, the Office for Research, and Northwestern University Information Technology.
Facility: Swift (XRT) - .
Software: Prospector (Leja et al. 2017), Python-fsps (Conroy et al. 2009; Conroy & Gunn 2010), Dynesty (Speagle 2020).
Appendix A: Photometric Redshift Sample Prospector Fits
As discussed in Section 4.1, the photometric redshifts inferred from the Prospector fitting have a large impact on the full redshift distribution of short GRB host galaxies, as they are noticed to exist at higher redshifts than the spectroscopic redshift sample. We show the model SEDs from the stellar population fits in Figure A1. as further evidence that their photometric redshifts are plausible.
Figure A1. Here, we present the Prospector-produced model spectra (purple lines), photometry (purple squares), and observed photometry (green circles) for the 20 short GRBs with no spectroscopic redshift. Upper limits are denoted with green triangles. We find excellent agreement among the model and observed SED, suggesting that photometric redshifts are in agreement with the true redshifts of the hosts.
Download figure:
Standard image High-resolution imageAppendix B: Literature Comparison of Past GRB Host Properties
We compare the short GRB host properties found in this work to previously published properties to explore the differences in stellar population modeling techniques that result in different values for the properties of the hosts. In Nugent et al. (2020), we found Prospector-derived masses, ages, and redshifts for 19 hosts (GRBs 050509B, 050709, 050724, 050813, 051210, 051221A, 060801, 061006, 061210, 070429B, 070714B, 070724, 070809, 071227, 080123, 090510, 090515, 100117, and 161104A). The Prospector model differs slightly from the one used in this work, as it does not include the mass–metallicity relation, old to young dust attenuation (2:1) ratio, nebular emission-line fitting (the gas-phase metallicity and gas ionization parameters), or the spectroscopic noise inflation parameter (when needed). Furthermore, our Prospector fits include the spectra for 12 of the GRB hosts in Nugent et al. (2020; see Table 2), which better informs the metallicity, age, and SFR of the host galaxies, whereas Nugent et al. (2020) only included the spectrum for GRB 161104A. We also find stellar population properties for 100206A (SFR, age, and mass; Berger 2014; Chrimes et al. 2018), 100625A, 101219A, 111117A, 120804A (SFRs, ages, and masses; Berger 2014), 130603B (SFR and mass; Cucchiara et al. 2013), 141212A (SFR and mass; Chrimes et al. 2018; Pandey et al. 2019), 150101B (SFR, age, and mass; Fong et al. 2016), 150120A (SFR and mass; Chrimes et al. 2018), 160821B (SFR and mass; Troja et al. 2019), 181123B (SFR, mass, and age; Paterson et al. 2020), and 200522A (SFR, mass, and age; Fong et al. 2021). We note that Paterson et al. (2020) and Fong et al. (2021) also used Prospector modeling, but did not include the mass–metallicity relation, as in this work. The rest of the stellar population properties are derived from B-band luminosities, spectral line strength, and/or other SED models with different assumed SFHs, including single stellar populations (SSPs) and a declining SFH with bursts of star formation.
We compare the available ages, SFRs, and stellar masses, to the median values in this work. We find large differences in the ages and SFRs of the hosts, which is likely caused by a combination of including spectra in our SED fits, mass-weighted ages, and SED as opposed to emission-line SFRs. Age estimates from SSPs are known to systematically underestimate the true age of the stellar populations, which was the method used to determine the ages for a majority of the previously published values. Leja et al. (2019) also showed that choice of SFH model has the largest effect on the inferred mass-weighted age, thus likely accounting for these differences. Out of the 27 new age estimates, 16 were derived from both the photometry and spectra of the hosts, whereas only two of the previously published values used both. Similarly, we performed joint spectroscopic and photometric fits for eight out of the 12 hosts for new SFR estimates, which had only been done for one of the hosts previously. The addition of the available spectra in our fits better constrains the ages and SFRs (and metallicity) of the galaxies than traditional photometry-only-based SED fittings, as shown in Johnson et al. (2021) and Tacchella et al. (2022). Indeed, we find that the largest changes in both age and SFR are with the hosts that have a joint spectroscopic and photometric fit. The largest changes in age are with the hosts of GRBs 050509B (≈8 times larger) and 100206A (≈40 times larger), and the largest change in SFR is with the host of 160821B (≈6 times lower). The hosts of GRBs 050509B, 100206A, and 160821B all have spectra with strong spectral lines used in their fits, which likely affected the SFR measurement. We also find that a majority (≈55%) of the new age estimates are higher, and a majority (≈59%) of SFR measurements are lower than previously determined. In total, the changes in stellar mass, age, and SFR are well understood from comparisons of different stellar population modeling techniques, and the uncertainties from Prospector should properly account for these differences.
We find that the parameter that differs the most between our work and previously literature is the SFR. Thus, we further explore SFR estimates in previously published works when they vary significantly from our inferred SFRs. We note that while we do find differences in the literature of short GRB hosts, we purposely used a method derived in the same way as a comparison field galaxy population to avoid comparing different SFR measurements. We determine that that the majority of discrepancies in SFR are due to using the emission lines exclusively for an estimate, an SED fit with a different model assumption and exclusively using the host photometry, or an SFR inferred from radio data alone. With the exception of GRB 071227, we find that all SFR estimates determined from emission-line strengths are higher than the SED-inferred SFR, as is expected for galaxies (Leja et al. 2019). For example, the hosts of GRBs 050709 and 101219A have emission-line SFRs ≈10 and ≈4 times larger, respectively, than the SED-inferred SFRs presented here (Fox et al. 2005; Covino et al. 2006; Prochaska et al. 2006; Fong et al. 2013; Nicuesa Guelbenzu et al. 2021). The host of GRB 071227 has a lower reported emission-line SFR in D'Avanzo et al. (2007) than in this work. There was no available spectrum to use in the Prospector fit for this host; hence, the SFR estimate was more unconstrained, likely leading to this difference. Furthermore, the radio-derived SFR for this host is much larger, (Nicuesa Guelbenzu et al. 2014) implying there is obscured star formation. As our fit includes NIR data, it likely accounts for this dust, hence why the derived SFR is larger than the emission-line SFR. We find that SFR estimates also substantially change with SED fitting model assumptions, including SFH choice and the use of single versus composite stellar population models (SSP versus CSP). Perley et al. (2012) and Chrimes et al. (2018) attempted an exponentially declining SFH for the host of GRB 100206A, a simpler SFH than the delayed-tau model chosen here, which resulted in SFRs two to three times larger than in this work. The use of SSPs, which have been shown to less accurately model the true stellar population properties of galaxies and result in younger and less-massive galaxies than expected, in de Ugarte Postigo et al. (2014) and Pandey et al. (2019) for the host of GRB 130603B resulted in SFR estimates ≈10 times greater than those found here. We find that the usage of simpler SFHs, such as exponentially declining SFHs, and the sole usage of photometry in the SED fits leads to higher inferred SFR estimates. For the host of GRB 100206A, the SFR estimate is two to three times larger in Perley et al. (2012) and Chrimes et al. (2018), and the host of GRB 130603B is ≈10 times greater in de Ugarte Postigo et al. (2014) and Pandey et al. (2019). Finally, as expected for obscured SF probed by radio observations (e.g., Perley & Perley 2013; Perley et al. 2015; Gatkine et al. 2020), we find the radio-inferred SFRs in the literature always lead to higher values than those from SEDs. This includes the hosts of GRB 071227 (Nicuesa Guelbenzu et al. 2021), 100206A (Klose et al. 2019), and 120804A (Berger et al. 2013b; for the latter burst, the lower redshift presented here likely also contributes to a lower inferred SFR).
We find the least variation in stellar mass. As mass is best determined through galaxy photometry and not the spectra, we expect little variation in stellar mass calculations between different SED methods (Johnson et al. 2021). The host of GRB 150101B appears to have the largest change in stellar mass estimate, increasing 1.2 times from the value in Fong et al. (2016). The previous estimate used a grid-searched SSP model to derive the stellar mass, which has been shown to underestimate the true mass value (Conroy 2013). We note that an SSP model for the host of GW/GRB 170817 resulted in a higher inferred stellar mass estimate in Levan et al. (2017) than that found in our work. Furthermore, this host has a known AGN (Xie et al. 2016), which, if not accounted for in the SED fitting, may affect the amount of dust in the infrared, which will lead to errors in the mass inference. We also find larger differences in stellar mass, whether increasing or decreasing, with the galaxies with previously unknown redshifts or current photometric redshifts, where the true SED shape is overall more uncertain. Overall, it does appear that mass estimates from various techniques do result in consistent values.
We finally compare our photometric redshifts to those found in O'Connor et al. (2022). O'Connor et al. (2022) also used Prospector to determine the redshifts for GRBs 140129B (we reported a spectroscopic redshift in Fong et al. 2022), 170728B (we reported a spectroscopic redshift in Fong et al. 2022), 180618A, 151229A, and 200411A. We also identify a different host galaxy for GRB 191031D due to our deeper NIR imaging in Fong et al. (2022). They used a restrictive redshift prior range of 0 ≤ z ≤ 1.5 for their analysis. We find this range of sampled values too conservative as Swift has detected short GRBs out to much farther confirmed redshifts (e.g., GRB 181123B at z = 1.7 and GRB 111117A at z = 2.2; Selsing et al. 2018; Paterson et al. 2020), suggesting that short GRBs can in principle be detected out to z ≈ 3. Furthermore, many of our photometric redshift sample hosts have posterior distributions that extend beyond z = 1.5; thus, we would not be fully sampling their possible redshifts and other stellar population properties by restricting the maximum to z = 1.5.
We emphasize that although we show a majority of the values are consistent with previously published values, it is best practice to compare properties derived from the same technique. Furthermore, as our SED fits include all available photometry, spectroscopy, and redshifts for 69 short GRB hosts, as well as informed priors, our sample represents the largest, uniformly modeled catalog to date.
Footnotes
- 15
- 16
We note that this sample also includes the long-duration GRBs 060614 and 211211A as there is significant evidence that they are derived from neutron star (NS) mergers and not young, stellar progenitors.
- 17
We typically set n = 10, although we increase it to 12 for spectra that cover a wider observed wavelength range (≥4000 Å) and have larger fluctuations in the observed continua, and decrease it to n = 6 for spectra that have a smaller wavelength coverage (≤3000 Å) or very flat continua. We note that using a higher n than necessary overfits the continuum, which can affect true spectral line strength.
- 18
These include GRBs 050724, 051221A, 061006, 101224A, 120305A, 140903A, 150728A, 170428A, and 201221D.
- 19
- 20
These GRB are 050724, 061006, 060614, 061210, 070714B, 080123, 090510, 170728B, 180418A, 180618A, 180805B, and 200219A.