• ### The Astropy Project: Building an inclusive, open-science project and status of the v2.0 core package(1801.02634)

Jan. 16, 2018 astro-ph.IM
The Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project.
• ### Early Spectra of the Gravitational Wave Source GW170817: Evolution of a Neutron Star Merger(1710.05432)

Oct. 16, 2017 astro-ph.HE
On 2017 August 17, Swope Supernova Survey 2017a (SSS17a) was discovered as the optical counterpart of the binary neutron star gravitational wave event GW170817. We report time-series spectroscopy of SSS17a from 11.75 hours until 8.5 days after merger. Over the first hour of observations the ejecta rapidly expanded and cooled. Applying blackbody fits to the spectra, we measure the photosphere cooling from $11,000^{+3400}_{-900}$ K to $9300^{+300}_{-300}$ K, and determine a photospheric velocity of roughly 30% of the speed of light. The spectra of SSS17a begin displaying broad features after 1.46 days, and evolve qualitatively over each subsequent day, with distinct blue (early-time) and red (late-time) components. The late-time component is consistent with theoretical models of r-process-enriched neutron star ejecta, whereas the blue component requires high velocity, lanthanide-free material.
• ### The Old Host-Galaxy Environment of SSS17a, the First Electromagnetic Counterpart to a Gravitational Wave Source(1710.05439)

We present an analysis of the host-galaxy environment of Swope Supernova Survey 2017a (SSS17a), the discovery of an electromagnetic counterpart to a gravitational wave source, GW170817. SSS17a occurred 1.9 kpc (in projection; 10.2") from the nucleus of NGC 4993, an S0 galaxy at a distance of 40 Mpc. We present a Hubble Space Telescope (HST) pre-trigger image of NGC 4993, Magellan optical spectroscopy of the nucleus of NGC 4993 and the location of SSS17a, and broad-band UV through IR photometry of NGC 4993. The spectrum and broad-band spectral-energy distribution indicate that NGC 4993 has a stellar mass of log (M/M_solar) = 10.49^{+0.08}_{-0.20} and star formation rate of 0.003 M_solar/yr, and the progenitor system of SSS17a likely had an age of >2.8 Gyr. There is no counterpart at the position of SSS17a in the HST pre-trigger image, indicating that the progenitor system had an absolute magnitude M_V > -5.8 mag. We detect dust lanes extending out to almost the position of SSS17a and >100 likely globular clusters associated with NGC 4993. The offset of SSS17a is similar to many short gamma-ray burst offsets, and its progenitor system was likely bound to NGC 4993. The environment of SSS17a is consistent with an old progenitor system such as a binary neutron star system.
• ### Light Curves of the Neutron Star Merger GW170817/SSS17a: Implications for R-Process Nucleosynthesis(1710.05443)

On 2017 August 17, gravitational waves were detected from a binary neutron star merger, GW170817, along with a coincident short gamma-ray burst, GRB170817A. An optical transient source, Swope Supernova Survey 17a (SSS17a), was subsequently identified as the counterpart of this event. We present ultraviolet, optical and infrared light curves of SSS17a extending from 10.9 hours to 18 days post-merger. We constrain the radioactively-powered transient resulting from the ejection of neutron-rich material. The fast rise of the light curves, subsequent decay, and rapid color evolution are consistent with multiple ejecta components of differing lanthanide abundance. The late-time light curve indicates that SSS17a produced at least ~0.05 solar masses of heavy elements, demonstrating that neutron star mergers play a role in r-process nucleosynthesis in the Universe.
• ### Swope Supernova Survey 2017a (SSS17a), the Optical Counterpart to a Gravitational Wave Source(1710.05452)

Oct. 16, 2017 astro-ph.HE
On 2017 August 17, the Laser Interferometer Gravitational-wave Observatory (LIGO) and the Virgo interferometer detected gravitational waves emanating from a binary neutron star merger, GW170817. Nearly simultaneously, the Fermi and INTEGRAL telescopes detected a gamma-ray transient, GRB 170817A. 10.9 hours after the gravitational wave trigger, we discovered a transient and fading optical source, Swope Supernova Survey 2017a (SSS17a), coincident with GW170817. SSS17a is located in NGC 4993, an S0 galaxy at a distance of 40 megaparsecs. The precise location of GW170817 provides an opportunity to probe the nature of these cataclysmic events by combining electromagnetic and gravitational-wave observations.
• The detection of GW170817 in both gravitational waves and electromagnetic waves heralds the age of gravitational-wave multi-messenger astronomy. On 17 August 2017 the Advanced LIGO and Virgo detectors observed GW170817, a strong signal from the merger of a binary neutron-star system. Less than 2 seconds after the merger, a gamma-ray burst (GRB 170817A) was detected within a region of the sky consistent with the LIGO-Virgo-derived location of the gravitational-wave source. This sky region was subsequently observed by optical astronomy facilities, resulting in the identification of an optical transient signal within $\sim 10$ arcsec of the galaxy NGC 4993. These multi-messenger observations allow us to use GW170817 as a standard siren, the gravitational-wave analog of an astronomical standard candle, to measure the Hubble constant. This quantity, which represents the local expansion rate of the Universe, sets the overall scale of the Universe and is of fundamental importance to cosmology. Our measurement combines the distance to the source inferred purely from the gravitational-wave signal with the recession velocity inferred from measurements of the redshift using electromagnetic data. This approach does not require any form of cosmic "distance ladder;" the gravitational wave analysis can be used to estimate the luminosity distance out to cosmological scales directly, without the use of intermediate astronomical distance measurements. We determine the Hubble constant to be $70.0^{+12.0}_{-8.0} \, \mathrm{km} \, \mathrm{s}^{-1} \, \mathrm{Mpc}^{-1}$ (maximum a posteriori and 68% credible interval). This is consistent with existing measurements, while being completely independent of them. Additional standard-siren measurements from future gravitational-wave sources will provide precision constraints of this important cosmological parameter.
• ### Gas inflow and outflow in an interacting high-redshift galaxy: The remarkable host environment of GRB 080810 at $z=3.35$(1705.01543)

July 17, 2017 astro-ph.GA
We reveal multiple components of an interacting galaxy system at $z\approx3.35$ through a detailed analysis of the exquisite high-resolution Keck/HIRES spectrum of the afterglow of a gamma-ray burst (GRB). Through Voigt-profile fitting of absorption lines from the Lyman-series, we constrain the neutral hydrogen column density to $N_{\mathrm{HI}} \leq 10^{18.35}$ cm$^{-2}$ for the densest of four distinct systems at the host redshift of GRB~080810, among the lowest $N_{\mathrm{HI}}$ ever observed in a GRB host, despite the line of sight passing within a projected 5 kpc of the galaxy centres. By detailed analysis of the corresponding metal absorption lines, we derive chemical, ionic and kinematic properties of the individual absorbing systems, and thus build a picture of the host as a whole. Striking differences between the systems imply that the line of sight passes through several phases of gas: the star-forming regions of the GRB host; enriched material in the form of a galactic outflow; the hot and ionised halo of a second, interacting galaxy falling towards the host at a line-of-sight velocity of 700 km s$^{-1}$; and a cool, metal-poor cloud which may represent one of the best candidates yet for the inflow of metal-poor gas from the intergalactic medium.
• ### The survival of gas clouds in the Circumgalactic Medium of Milky Way-like galaxies(1608.05416)

May 19, 2017 astro-ph.GA
Observational evidence shows that low-redshift galaxies are surrounded by extended haloes of multiphase gas, the so-called 'circumgalactic medium' (CGM). To study the survival of relatively cool gas (T < 10^5 K) in the CGM, we performed a set of hydrodynamical simulations of cold (T = 10^4 K) neutral gas clouds travelling through a hot (T = 2x10^6 K) and low-density (n = 10^-4 cm^-3) coronal medium, typical of Milky Way-like galaxies at large galactocentric distances (~ 50-150 kpc). We explored the effects of different initial values of relative velocity and radius of the clouds. Our simulations were performed on a two-dimensional grid with constant mesh size (2 pc) and they include radiative cooling, photoionization heating and thermal conduction. We found that for large clouds (radii larger than 250 pc) the cool gas survives for very long time (larger than 250 Myr): despite that they are partially destroyed and fragmented into smaller cloudlets during their trajectory, the total mass of cool gas decreases at very low rates. We found that thermal conduction plays a significant role: its effect is to hinder formation of hydrodynamical instabilities at the cloud-corona interface, keeping the cloud compact and therefore more difficult to destroy. The distribution of column densities extracted from our simulations are compatible with those observed for low-temperature ions (e.g. SiII and SiIII) and for high-temperature ions (OVI) once we take into account that OVI covers much more extended regions than the cool gas and, therefore, it is more likely to be detected along a generic line of sight.
• ### Photometric Observations of Supernova 2013cq Associated with GRB 130427A(1702.04762)

Feb. 15, 2017 astro-ph.HE
We observed the afterglow of GRB 130427A with the RATIR instrument on the 1.5-m Harold L. Johnson telescope of the Observatorio Astron\'omico Nacional on Sierra San Pedro M\'artir. Our homogenous $griZYJH$ photometry extends from the night of burst to three years later. We fit a model for the afterglow. There is a significant positive residual which matches the behavior of SN 1998bw in the $griZ$ filters; we suggest that this is a photometric signature of the supernova SN 2013cq associated with the GRB. The peak absolute magnitude of the supernova is $M_r=-18.43\pm0.11$.
• ### The Circumgalactic Medium of Submillimeter Galaxies. I. First Results from a Radio-Identified Sample(1607.00016)

Aug. 19, 2016 astro-ph.GA
We present the first results from an on-going survey to characterize the circumgalactic medium (CGM) of the massive high-redshift galaxies detected as submillimeter galaxies (SMGs). We constructed a parent sample of 163 SMG-QSO pairs with separations less than $\sim$36 arcsec by cross-matching far-infrared-selected galaxies from Herschel with spectroscopically confirmed QSOs. The Herschel sources were selected to match the properties of SMGs. We determined the sub-arcsecond positions of six Herschel sources with the Very Large Array and obtained secure redshift identification for three of those with near-infrared spectroscopy. The QSO sightlines probe transverse proper distances of 112, 157, and 198 kpc at foreground redshifts of 2.043, 2.515, and 2.184, respectively, which are comparable to the virial radius of the $\sim10^{13}$ Msun halos expected to host SMGs. High-quality absorption-line spectroscopy of the QSOs reveals systematically strong HI Lyman-alpha absorption around all three SMGs, with rest-frame equivalent widths of $\sim2-3$ \AA. However, none of the three absorbers exhibits compelling evidence for optically thick HI gas or metal absorption, in contrast to the dominance of strong neutral absorbers in the CGM of luminous $z \sim 2$ QSOs. The low covering factor of optically thick HI gas around SMGs tentatively indicates that SMGs may not have as prominent cool gas reservoirs in their halos as the co-eval QSOs and that they may inhabit less massive halos than previously thought.
• ### XQ-100: A legacy survey of one hundred 3.5 < z < 4.5 quasars observed with VLT/XSHOOTER(1607.08776)

Aug. 12, 2016 astro-ph.GA
We describe the execution and data reduction of the European Southern Observatory Large Programme "Quasars and their absorption lines: a legacy survey of the high-redshift universe with VLT/XSHOOTER" (hereafter `XQ-100'). XQ-100 has produced and made publicly available an homogeneous and high-quality sample of echelle spectra of 100 QSOs at redshifts z~3.5-4.5 observed with full spectral coverage from 315 to 2500 nm at a resolving power ranging from R~4000 to 7000, depending on wavelength. The median signal-to-noise ratios are 33, 25 and 43, as measured at rest-frame wavelengths 1700, 3000 and 3600 Angstrom, respectively. This paper provides future users of XQ-100 data with the basic statistics of the survey, along with details of target selection, data acquisition and data reduction. The paper accompanies the public release of all data products, including 100 reduced spectra. XQ-100 is the largest spectroscopic survey to date of high-redshift QSOs with simultaneous rest-frame UV/optical coverage, and as such enables a wide range of extragalactic research, from cosmology and galaxy evolution to AGN astrophysics.
• ### Upper limits from five years of blazar observations with the VERITAS Cherenkov telescopes(1603.02410)

March 8, 2016 astro-ph.HE
Between the beginning of its full-scale scientific operations in 2007 and 2012, the VERITAS Cherenkov telescope array observed more than 130 blazars; of these, 26 were detected as very-high-energy (VHE; E>100 GeV) {\gamma}-ray sources. In this work, we present the analysis results of a sample of 114 undetected objects. The observations constitute a total live-time of ~570 hours. The sample includes several unidentified Fermi-Large Area Telescope (LAT) sources (located at high Galactic latitude) as well as all the sources from the second Fermi-LAT catalog which are contained within the field of view of the VERITAS observations. We have also performed optical spectroscopy measurements in order to estimate the redshift of some of these blazars that do not have a spectroscopic distance estimate. We present new optical spectra from the Kast instrument on the Shane telescope at the Lick observatory for 18 blazars included in this work, which allowed for the successful measurement or constraint on the redshift of four of them. For each of the blazars included in our sample we provide the flux upper limit in the VERITAS energy band. We also study the properties of the significance distributions and we present the result of a stacked analysis of the data-set, which shows a 4 {\sigma} excess.
• ### Discovery of a Damped Lya Absorber at z = 3.3 along a galaxy sight-line in the SSA22 field(1512.04625)

Jan. 31, 2016 astro-ph.GA
Using galaxies as background light sources to map the Lya absorption lines is a novel approach to study Damped Lya Absorbers (DLAs). We report the discovery of an intervening z = 3.335 +- 0.007 DLA along a galaxy sight-line identified among 80 Lyman Break Galaxy (LBG) spectra obtained with our Very Large Telescope/Visible Multi-Object Spectrograph survey in the SSA22 field. The measured DLA neutral hydrogen (HI) column density is log (NHI/cm^{-2}) = 21.68 +- 0.17. The DLA covering fraction over the extended background LBG is > 70% (2 sigma), yielding a conservative constraint on the DLA area of > 1 kpc^2. Our search for a counterpart galaxy hosting this DLA concludes that there is no counterpart galaxy with star formation rate (SFR) larger than a few Msun yr^{-1}, ruling out an unobscured violent star formation in the DLA gas cloud. We also rule out the possibility that the host galaxy of the DLA is a passive galaxy with Mstar > 5 x 10^{10} Msun or a heavily dust-obscured galaxy with E(B-V) > 2. The DLA may coincide with a large-scale overdensity of the spectroscopic LBGs. The occurrence rate of the DLA is compatible with that of DLAs found in QSO sight-lines.
• ### Optical and near-infrared observations of SN 2013dx associated with GRB 130702A(1508.00575)

Jan. 6, 2016 astro-ph.HE
We present optical and near-infrared light curves and optical spectra of SN 2013dx, associated with the nearby (redshift 0.145) gamma-ray burst GRB 130702A. The prompt isotropic gamma-ray energy released from GRB 130702A is measured to be $E_{\gamma,iso}=6.4_{-1.0}^{+1.3}\times10^{50}$erg (1keV-10MeV in the rest frame), placing it intermediate between low-luminosity GRBs like GRB 980425/SN 1998bw and the broader cosmological population. We compare the observed $g'r'i'z'$ light curves of SN 2013dx to a SN 1998bw template, finding that SN 2013dx evolves ~20% faster (steeper rise time), with a comparable peak luminosity. Spectroscopically, SN 2013dx resembles other broad-lined Type Ic supernovae, both associated with (SN 2006aj and SN 1998bw) and lacking (SN 1997ef, SN 2007I, and SN 2010ah) gamma-ray emission, with photospheric velocities around peak of ~21,000km s$^{-1}$. We construct a quasi-bolometric ($g'r'i'z'yJ$) light curve for SN 2013dx, only the fifth GRB-associated SN with extensive NIR coverage and the third with a bolometric light curve extending beyond $\Delta t>40$d. Together with the measured photospheric velocity, we derive basic explosion parameters using simple analytic models. We infer a $^{56}$Ni mass of $M_{\mathrm{Ni}}=0.37\pm0.01$M$_{\odot}$, an ejecta mass of $M_{\mathrm{ej}}=3.1\pm0.1$M$_{\odot}$, and a kinetic energy of $E_{\mathrm{K}}=(8.2\pm0.43)\times10^{51}$ erg (statistical uncertainties only), consistent with previous GRB-associated SNe. When considering the ensemble population of GRB-associated SNe, we find no correlation between the mass of synthesized $^{56}$Ni and high-energy properties, despite clear predictions from numerical simulations that $M_{\mathrm{Ni}}$ should correlate with the degree of asymmetry. On the other hand, $M_{\mathrm{Ni}}$ clearly correlates with the kinetic energy of the supernova ejecta across a wide range of core-collapse events.
• ### The evolution of neutral gas in damped Lyman $\alpha$ systems from the XQ-100 survey(1511.05003)

Nov. 16, 2015 astro-ph.CO, astro-ph.GA
We present a sample of 38 intervening Damped Lyman $\alpha$ (DLA) systems identified towards 100 $z>3.5$ quasars, observed during the XQ-100 survey. The XQ-100 DLA sample is combined with major DLA surveys in the literature. The final combined sample consists of 742 DLAs over a redshift range approximately $1.6 < z_{\rm abs} < 5.0$. We develop a novel technique for computing $\Omega_{\rm HI}^{\rm DLA}$ as a continuous function of redshift, and we thoroughly assess and quantify the sources of error therein, including fitting errors and incomplete sampling of the high column density end of the column density distribution function. There is a statistically significant redshift evolution in $\Omega_{\rm HI}^{\rm DLA}$ ($\geq 3 \sigma$) from $z \sim 2$ to $z \sim$ 5. In order to make a complete assessment of the redshift evolution of $\Omega_{\rm HI}$, we combine our high redshift DLA sample with absorption surveys at intermediate redshift and 21cm emission line surveys of the local universe. Although $\Omega_{\rm HI}^{\rm DLA}$, and hence its redshift evolution, remains uncertain in the intermediate redshift regime ($0.1 < z_{\rm abs} < 1.6$), we find that the combination of high redshift data with 21cm surveys of the local universe all yield a statistically significant evolution in $\Omega_{\rm HI}$ from $z \sim 0$ to $z \sim 5$ ($\geq 3 \sigma$). Despite its statistical significance, the magnitude of the evolution is small: a linear regression fit between $\Omega_{\rm HI}$ and $z$ yields a typical slope of $\sim$0.17$\times 10^{-3}$, corresponding to a factor of $\sim$ 4 decrease in $\Omega_{\rm HI}$ between $z=5$ and $z=0$.
• ### Happy Birthday Swift: Ultra-long GRB141121A and its broad-band Afterglow(1510.00996)

Oct. 4, 2015 astro-ph.HE
We present our extensive observational campaign on the Swift-discovered GRB141121A, al- most ten years after its launch. Our observations covers radio through X-rays, and extends for more than 30 days after discovery. The prompt phase of GRB 141121A lasted 1410 s and, at the derived redshift of z = 1.469, the isotropic energy is E{\gamma},iso = 8.0x10^52 erg. Due to the long prompt duration, GRB141121A falls into the recently discovered class of UL-GRBs. Peculiar features of this burst are a flat early-time optical light curve and a radio-to-X-ray rebrightening around 3 days after the burst. The latter is followed by a steep optical-to-X-ray decay and a much shallower radio fading. We analyze GRB 141121A in the context of the standard forward-reverse shock (FS,RS) scenario and we disentangle the FS and RS contributions. Finally, we comment on the puzzling early-time (t ~3 d) behavior of GRB 141121A, and suggest that its interpretation may require a two-component jet model. Overall, our analysis confirms that the class of UL-GRBs represents our best opportunity to firmly establish the prominent emission mechanisms in action during powerful GRB explosions, and future missions (like SVOM, XTiDE, or ISS-Lobster) will provide many more of such objects.
• ### The central engine of GRB 130831A and the energy breakdown of a relativistic explosion(1509.09234)

Sept. 30, 2015 astro-ph.HE
Gamma-ray bursts (GRBs) are the most luminous explosions in the universe, yet the nature and physical properties of their energy sources are far from understood. Very important clues, however, can be inferred by studying the afterglows of these events. We present optical and X-ray observations of GRB 130831A obtained by Swift, Chandra, Skynet, RATIR, Maidanak, ISON, NOT, LT and GTC. This burst shows a steep drop in the X-ray light-curve at $\simeq 10^5$ s after the trigger, with a power-law decay index of $\alpha \sim 6$. Such a rare behaviour cannot be explained by the standard forward shock (FS) model and indicates that the emission, up to the fast decay at $10^5$ s, must be of "internal origin", produced by a dissipation process within an ultrarelativistic outflow. We propose that the source of such an outflow, which must produce the X-ray flux for $\simeq 1$ day in the cosmological rest frame, is a newly born magnetar or black hole. After the drop, the faint X-ray afterglow continues with a much shallower decay. The optical emission, on the other hand, shows no break across the X-ray steep decrease, and the late-time decays of both the X-ray and optical are consistent. Using both the X-ray and optical data, we show that the emission after $\simeq 10^5$ s can be explained well by the FS model. We model our data to derive the kinetic energy of the ejecta and thus measure the efficiency of the central engine of a GRB with emission of internal origin visible for a long time. Furthermore, we break down the energy budget of this GRB into the prompt emission, the late internal dissipation, the kinetic energy of the relativistic ejecta, and compare it with the energy of the associated supernova, SN 2013fu.
• ### The first ultraviolet quasar stacked spectrum at z=2.4 from WFC3(1503.02075)

April 7, 2015 astro-ph.GA
The ionising continuum from active galactic nuclei (AGN) is fundamental for interpreting their broad emission lines and understanding their impact on the surrounding gas. Furthermore, it provides hints on how matter accretes onto supermassive black holes. Using HST's Wide Field Camera 3 we have constructed the first stacked ultraviolet (rest-frame wavelengths 600-2500\AA) spectrum of 53 luminous quasars at z=2.4, with a state-of-the-art correction for the intervening Lyman forest and Lyman continuum absorption. The continuum slope ($f_\nu \propto \nu^{\alpha_\nu}$) of the full sample shows a break at ~912\AA\ with spectral index $\alpha_\nu=-0.61\pm0.01$ at $\lambda>912$\AA\ and a softening at shorter wavelengths ($\alpha_\nu=-1.70 \pm 0.61$ at $\lambda\leq 912$\AA). Our analysis proves that a proper intergalactic medium absorption correction is required to establish the intrinsic continuum emission of quasars. We interpret our average ultraviolet spectrum in the context of photoionisation, accretion disk models, and quasar contribution to the ultraviolet background. We find that observed broad line ratios are consistent with those predicted assuming an ionising slope of $\alpha_\mathrm{ion}=$-2.0, similar to the observed ionising spectrum in the same wavelength range. The continuum break and softening are consistent with accretion disk plus X-ray corona models when black hole spin is taken into account. Our spectral energy distribution yields a 30% increase to previous estimates of the specific quasar emissivity, such that quasars may contribute significantly to the total specific Lyman limit emissivity estimated from the Ly$\alpha$ forest at z<3.2.
• ### A detailed study of the optical attenuation of gamma-ray bursts in the Swift era(1412.6530)

March 3, 2015 astro-ph.HE
We present optical and near-infrared (NIR) photometry of 28 gamma-ray bursts (GRBs) detected by the \textit{Swift} satellite and rapidly observed by the Reionization and Transients Infrared/Optical (RATIR) camera. We compare the optical flux at fiducial times of 5.5 and 11 hours after the high-energy trigger to that in the X-ray regime to quantify optical darkness. 46$\pm$9 per cent (13/28) of all bursts in our sample and 55$\pm$10 per cent (13/26) of long GRBs are optically dark, which is statistically consistently with previous studies. Fitting RATIR optical and NIR spectral energy distributions (SEDs) of 19 GRBs, most (6/7) optically dark GRBs either occur at high-redshift ($z>4.5$) or have a high dust content in their host galaxies ($A_{\rm V} > 0.3$). Performing K-S tests, we compare the RATIR sample to those previously presented in the literature, finding our distributions of redshift, optical darkness, host dust extinction and X-ray derived column density to be consistent. The one reported discrepancy is with host galaxy dust content in the BAT6 sample, which appears inconsistent with our sample and other previous literature. Comparing X-ray derived host galaxy hydrogen column densities to host galaxy dust extinction, we find that GRBs tend to occur in host galaxies with a higher metal-to-dust ratio than our own Galaxy, more akin to the Large and Small Magellanic Clouds. Finally, to mitigate time evolution of optical darkness, we measure $\beta_{\rm OX,rest}$ at a fixed rest frame time, $t_{\rm rest}=1.5$ hours and fixed rest frame energies in the X-ray and optical regimes. Choosing to evaluate optical flux at $\lambda_{\rm rest}=0.25~\mu$m, we remove high-redshift as a source of optical darkness, demonstrating that optical darkness must result from either high-redshift, dust content in the host galaxy along the GRB sight line, or a combination of the two.
• ### Unveiling the Secrets of Metallicity and Massive Star Formation Using DLAs along Gamma-ray Bursts(1408.3578)

We present the largest, publicly available, sample of Damped Lyman-$\alpha$ systems (DLAs) along Gamma-ray Bursts (GRB) line of sights in order to investigate the environmental properties of long GRBs in the $z=1.8-6$ redshift range. Compared with the most recent quasar DLAs sample (QSO-DLA), our analysis shows that GRB-DLAs probe a more metal enriched environment at $z\gtrsim3$, up to $[X/H]\sim-0.5$. In the $z=2-3$ redshift range, despite the large number of lower limits, there are hints that the two populations may be more similar (only at 90\% significance level). Also at \hiz, the GRB-DLA average metallicity seems to decline at a shallower rate than the QSO-DLAs: GRB-DLA hosts may be polluted with metals at least as far as $\sim 2$kpc from the GRB explosion site, probably due to previous star-formation episodes and/or supernovae explosions. This shallow metallicity trend, extended now up to $z\sim5$, confirms previous results that GRB hosts are star-forming and have, on average, higher metallicity than the general QSO-DLA population. Finally, our metallicity measurements are broadly consistent with the hypothesis of two channels of GRB progenitors, one of which is mildly affected by a metallicity bias. The metallicity evolution of modeled GRB hosts agrees reasonably well with our data up to intermediate redshift, while more data are needed to constrain the models at $z\gtrsim 4$.
• ### Identifying high-redshift GRBs with RATIR(1312.3967)

April 21, 2014 astro-ph.CO, astro-ph.HE
We present a template fitting algorithm for determining photometric redshifts, $z_{\rm phot}$, of candidate high-redshift gamma-ray bursts (GRBs). Using afterglow photometry, obtained by the Reionization And Transients InfraRed (RATIR) camera, this algorithm accounts for the intrinsic GRB afterglow spectral energy distribution (SED), host dust extinction and the effect of neutral hydrogen (local and cosmological) along the line of sight. We present the results obtained by this algorithm and RATIR photometry of GRB 130606A, finding a range of best fit solutions $5.6 < z_{\rm phot} < 6.0$ for models of several host dust extinction laws (none, MW, LMC and SMC), consistent with spectroscopic measurements of the redshift of this GRB. Using simulated RATIR photometry, we find our algorithm provides precise measures of $z_{\rm phot}$ in the ranges $4 < z_{\rm phot} \lesssim 8$ and $9 < z_{\rm phot} < 10$ and can robustly determine when $z_{\rm phot}>4$. Further testing highlights the required caution in cases of highly dust extincted host galaxies. These tests also show that our algorithm does not erroneously find $z_{\rm phot} < 4$ when $z_{\rm sim}>4$, thereby minimizing false negatives and allowing us to rapidly identify all potential high-redshift events.
• ### A Population of Massive, Luminous Galaxies Hosting Heavily Dust-Obscured Gamma-Ray Bursts: Implications for the Use of GRBs as Tracers of Cosmic Star Formation(1301.5903)

Oct. 30, 2013 astro-ph.CO, astro-ph.HE
We present observations and analysis of the host galaxies of 23 heavily dust-obscured gamma-ray bursts (GRBs) observed by the Swift satellite during the years 2005-2009, representing all GRBs with an unambiguous host-frame extinction of A_V>1 mag from this period. Deep observations with Keck, Gemini, VLT, HST, and Spitzer successfully detect the host galaxies and establish spectroscopic or photometric redshifts for all 23 events, enabling us to provide measurements of the intrinsic host star-formation rates, stellar masses, and mean extinctions. Compared to the hosts of unobscured GRBs at similar redshifts, we find that the hosts of dust-obscured GRBs are (on average) more massive by about an order of magnitude and also more rapidly star-forming and dust-obscured. While this demonstrates that GRBs populate all types of star-forming galaxies including the most massive, luminous systems at z~2, at redshifts below 1.5 the overall GRB population continues to show a highly significant aversion away from massive galaxies and a preference for low-mass systems relative to what would be expected given a purely SFR-selected galaxy sample. This supports the notion that the GRB rate is strongly dependent on metallicity, and may suggest that the most massive galaxies in the Universe underwent a transition in their chemical properties ~9 Gyr ago. We also conclude that, based on the absence of unobscured GRBs in massive galaxies and the absence of obscured GRBs in low-mass galaxies, the dust distributions of the lowest-mass and the highest-mass galaxies are relatively homogeneous, while intermediate-mass galaxies (~10^9 M_sun) have diverse internal properties.
• ### Gemini Spectroscopy of the Short GRB 130603B Afterglow and Host(1306.2028)

July 29, 2013 astro-ph.CO
We present early optical photometry and spectroscopy of the afterglow and host galaxy of the bright short-duration gamma-ray burst GRB 130603B discovered by the Swift satellite. Using the Gemini South telescope, our prompt optical spectra reveal a strong trace from the afterglow superimposed on continuum and emission lines from the $z = 0.3568 \pm 0.0005$ host galaxy. The combination of a relatively bright optical afterglow (r' = 21.52 at $\Delta_t$= 8.4hr), together with an observed offset of 0\farcs9 from the host nucleus (4.8kpc projected distance at z=0.3568), allow us to extract a relatively clean spectrum dominated by afterglow light . The spatially resolved spectrum allows us to constrain the properties of the explosion site directly, and compare these with the host galaxy nucleus, as well as other short-duration GRB host galaxies. We find that while the host is a relatively luminous ($L \approx 0.8 L^{*}_{B}$), star-forming galaxy with solar metallicity, the spectrum of the afterglow exhibits weak CaII absorption features but negligible emission features. The explosion site therefore lacks evidence of recent star formation, consistent with the relatively long delay time distribution expected in a compact binary merger scenario. The star formation rate (both in an absolute sense and normalized to the luminosity) and metallicity of the host are both consistent with the known sample of short-duration GRB hosts and with recent results which suggest GRB130603B emission to be the product of the decay of radioactive species produced during the merging process of a NS-NS binary ("kilonova"). Ultimately, the discovery of more events similar to GRB130603B and their rapid follow-up from 8-m class telescopes will open new opportunities for our understanding of the final stages of compact-objects binary systems.
• ### The Firm Redshift Lower Limit of the Most Distant TeV-Detected Blazar PKS 1424+240(1304.4859)

April 23, 2013 astro-ph.HE
We present the redshift lower limit of z>0.6035 for the very-high-energy (VHE; E>100 GeV) emitting blazar PKS 1424+240 (PG 1424+240). This limit is inferred from Lyman beta and gamma absorption observed in the far-ultraviolet spectra from the Hubble Space Telescope/Cosmic Origins Spectrograph. No VHE-detected blazar has shown solid spectroscopic evidence of being more distant. At this distance, VHE observations by VERITAS are shown to sample historically large gamma-ray opacity values at 500 GeV, extending beyond tau=4 for low-level models of the extragalactic background light (EBL) and beyond tau=5 for high-levels. The majority of the z=0.6035 absorption-corrected VHE spectrum appears to exhibit a lower flux than an extrapolation of the contemporaneous LAT power-law fit beyond 100 GeV. However, the highest energy VERITAS point is the only point showing agreement with this extrapolation, possibly implying the overestimation of the gamma-ray opacity or the onset of an unexpected VHE spectral feature. A curved log parabola is favored when fitting the full range of gamma-ray data (0.5 to 500 GeV). While fitting the absorption-corrected VHE data alone results in a harder differential power law than that from the full range, the indices derived using three EBL models are consistent with the physically motivated limit set by Fermi acceleration processes.
• ### An Independent Measurement of the Incidence of MgII Absorbers along Gamma-Ray Burst Sightlines: the End of the Mystery?(1211.6528)

Nov. 28, 2012 astro-ph.CO
In 2006, Prochter et al. reported a statistically significant enhancement of very strong Mg II absorption systems intervening the sightlines to gamma-ray bursts (GRBs) relative to the in- cidence of such absorption along quasar sightlines. This counterintuitive result, has inspired a diverse set of astrophysical explanations (e.g. dust, gravitational lensing) but none of these has obviously resolved the puzzle. Using the largest set of GRB afterglow spectra available, we reexamine the purported enhancement. In an independent sample of GRB spectra with a survey path 3 times larger than Prochter et al., we measure the incidence per unit redshift of $\geq 1$\AA rest-frame equivalent width Mg II absorbers at $z \approx 1$ to be l(z)= 0.18 $\pm$ 0.06. This is fully consistent with current estimates for the incidence of such absorbers along quasar sightlines. Therefore, we do not confirm the original enhancement and suggest those results suffered from a statistical fluke. Signatures of the original result do remain in our full sample (l(z) shows an $\approx 1.5$ enhancement over l(z)QSO), but the statistical significance now lies at $\approx 90%$ c.l. Restricting our analysis to the subset of high-resolution spectra of GRB afterglows (which overlaps substantially with Prochter et al.), we still reproduce a statistically significant enhancement of Mg II absorption. The reason for this excess, if real, is still unclear since there is no connection between the rapid afterglow follow-up process with echelle (or echellette) spectrographs and the detectability of strong Mg II doublets. Only a larger sample of such high-resolution data will shed some light on this matter.