Source author record

Darach Watson

Darach Watson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

32works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

32 published item(s)

preprint2026arXiv

On the origins of oxygen: ALMA and JWST characterise the multi-phase, metal-enriched, star-bursting medium within a 'normal' $z > 11$ galaxy

The unexpectedly high abundance of galaxies at $z > 11$ revealed by JWST has sparked a debate on the nature of early galaxies and the physical mechanisms regulating their formation. The Atacama Large Millimeter/submillimeter Array (ALMA) has begun to provide vital insights on their gas and dust content, but so far only for extreme 'blue monsters'. Here we present new, deep ALMA observations of JADES-GS-z11-0, a more typical (sub-$L^*$) $z > 11$ galaxy that bridges the discovery space of JWST and the Hubble Space Telescope. These data confirm the presence of the [O III] 88 $μ$m line at $4.5σ$ significance, precisely at the redshift of several faint emission lines previously seen with JWST/NIRSpec, while the underlying dust continuum remains undetected ($F_ν< 9.0 \, \mathrm{μJy}$), implying an obscured star formation rate (SFR) of $\text{SFR}_\text{IR} \lesssim 6 \, \mathrm{M_\odot \, yr^{-1}}$ and dust mass of $M_\text{dust} \lesssim 1.0 \times 10^{6} \, \mathrm{M_\odot}$ (all $3σ$). The accurate ALMA redshift of $z_\text{[O III]} = 11.1221 \pm 0.0006$ ($\gtrsim \! 5\times$ refined over NIRSpec) helps confirm that redshifts measured purely from the Lyman-$α$ break, even spectroscopically, should properly take into account the effects of potential damped Lyman-$α$ absorption (DLA) systems to avoid systematic overestimates of up to $Δz \approx 0.5$. The [O III] 88 $μ$m luminosity of $L_\text{[O III]} = (1.1 \pm 0.3) \times 10^{8} \, \mathrm{L_\odot}$, meanwhile, agrees well with the scaling relation for local metal-poor dwarfs given the SFR measured by NIRCam, NIRSpec, and MIRI. The spatially resolved MIRI and ALMA emission also underscores that JADES-GS-z11-0 is likely to consist of two low-mass components that are undergoing strong bursts of star formation yet are already pre-enriched in oxygen (~20-30% solar), only 400 Myr after the Big Bang.

preprint2022arXiv

ALMA reveals extended cool gas and hot ionized outflows in a typical star-forming galaxy at $z=7.13$

We present spatially-resolved morphological properties of [CII] 158 $μ$m, [OIII] 88 $μ$m, dust, and rest-frame ultraviolet (UV) continuum emission for A1689-zD1, a strongly lensed, sub-L* galaxy at $z=7.13$, by utilizing deep Atacama Large Millimeter/submillimeter Array (ALMA) and Hubble Space Telescope (HST) observations. While the [OIII] line and UV continuum are compact, the [CII] line is extended up to a radius of $r \sim 12$ kpc. Using multi-band rest-frame far-infrared (FIR) continuum data ranging from 52-400 $μ$m, we find an average dust temperature and emissivity index of $T_{\rm dust} = 41^{+17}_{-14}$ K and $β= 1.7^{+1.1}_{-0.7}$, respectively, across the galaxy. We find slight differences in the dust continuum profiles at different wavelengths, which may indicate that the dust temperature decreases with distance. We map the star-formation rate (SFR) via IR and UV luminosities and determine a total SFR of $37\pm 1~M_\odot~{\rm yr}^{-1}$ with an obscured fraction of $87\%$. While the [OIII] line is a good tracer of the SFR, the [CII] line shows deviation from the local $L_{\rm [CII]}$-SFR relations in the outskirts of the galaxy. Finally, we observe a clear difference in the line profile between [CII] and [OIII], with significant residuals ($\sim 5σ$) in the [OIII] line spectrum after subtracting a single Gaussian model. This suggests a possible origin of the extended [CII] structure from the cooling of hot ionized outflows. The extended [CII] and high-velocity [OIII] emission may both contribute in part to the high $L_{\rm [OIII]}$/$L_{\rm [CII]}$ ratios recently reported in $z>6$ galaxies.

preprint2022arXiv

The archival discovery of a strong Lyman-$α$ and [CII] emitter at z = 7.677

We report the archival discovery of Lyman-$α$ emission from the bright ultraviolet galaxy Y002 at $z=7.677$, spectroscopically confirmed by its ionized carbon [CII] 158$μ$m emission line. The Ly$α$ line is spatially associated with the rest-frame UV stellar emission ($M_{\rm UV}$~-22, 2x brighter than $M^\star_{\rm UV}$) and it appears offset from the peak of the extended [CII] emission at the current ~1" spatial resolution. We derive an estimate of the unobscured SFR(UV)=$(22\pm1)\,M_\odot$ yr$^{-1}$ and set an upper limit of SFR(IR)$<15\,M_\odot$ yr$^{-1}$ from the far-infrared wavelength range, which globally place Y002 on the SFR(UV+IR)-L([CII]) correlation observed at lower redshifts. In terms of velocity, the peak of the Ly$α$ emission is redshifted by $Δv$(Ly$α$)~500 km s$^{-1}$ from the systemic redshift set by [CII] and a high-velocity tail extends to up to ~1000 km s$^{-1}$. The velocity offset is up to ~3.5x higher than the average estimate for similarly UV-bright emitters at z~6-7, which might suggest that we are witnessing the merging of two clumps. A combination of strong outflows and the possible presence of an extended ionized bubble surrounding Y002 would likely facilitate the escape of copious Ly$α$ light, as indicated by the large equivalent width EW(Ly$α$)=$24^{+5}_{-6}$ Å. Assuming that [CII] traces the neutral hydrogen, we estimate a HI gas fraction of $M({\rm HI})/M_\star\gtrsim8$ for Y002 as a system and speculate that patches of high HI column densities could contribute to explain the observed spatial offsets between Ly$α$ and [CII] emitting regions. The low dust content, implied by the non-detection of the far-infrared continuum emission at rest-frame ~160 $μ$m, would be sufficient to absorb any potential Ly$α$ photons produced within the [CII] clump as a result of large HI column densities.

preprint2020arXiv

An Unambiguous Separation of Gamma-Ray Bursts into Two Classes from Prompt Emission Alone

The duration of a gamma-ray burst (GRB) is a key indicator of its physics origin, with long bursts perhaps associated with the collapse of massive stars and short bursts with mergers of neutron stars.However, there is substantial overlap in the properties of both short and long GRBs and neither duration nor any other parameter so far considered completely separates the two groups. Here we unambiguously classify every GRB using a machine-learning, dimensionality-reduction algorithm, t-distributed stochastic neighborhood embedding (t-SNE), providing a catalog separating all Swift GRBs into two groups. Although the classification takes place only using prompt emission light curves,every burst with an associated supernova is found in the longer group and bursts with kilonovae in the short, suggesting along with the duration distributions that these two groups are truly long and short GRBs. Two bursts with a clear absence of a supernova belong to the longer class, indicating that these might have been direct-collapse black holes, a proposed phenomenon that may occur in the deaths of more massive stars.

preprint2020arXiv

Direct measurement of the [CI] luminosity to molecular gas mass conversion factor in high-redshift star-forming galaxies

The amount of cold, molecular gas in high-redshift galaxies is typically inferred from proxies of molecular hydrogen (H$_2$), such as carbon monoxide (CO) or neutral atomic carbon ([CI]) and molecular gas mass conversion factors. The use of these proxies, however, relies on modeling and observations that have not been directly measured outside the local universe. Here, we use recent samples of high-redshift gamma-ray burst (GRB) and quasar molecular gas absorbers to determine this conversion factor $α_{\rm [CI]}$ from the column density of H$_2$, which gives us the mass per unit column, and the [CI]($J=1$) column density, which provides the luminosity per unit column. This technique allows us to make direct measurements of the relative abundances in high-redshift absorption-selected galaxies. Our sample spans redshifts of z=1.9-3.4 and covers two orders of magnitude in gas-phase metallicity. We find that $α_{\rm [CI]}$ scales linearly with the metallicity: $\log α_{\rm [CI]} = -1.13\times \log(Z/Z_{\odot}) + 1.33$, with a scatter of $σ_{α_{\rm [CI]}} = 0.2$ dex. Using a sample of emission-selected galaxies at z~0-5, with both [CI] and CO line detections, we apply the $α_{\rm [CI]}$ conversion to derive independent estimates of the molecular gas mass and $α_{\rm CO}$. We find a remarkable agreement between the molecular gas masses inferred from the absorption-derived $α_{\rm [CI]}$ compared to typical $α_{\rm CO}$-based estimates, and an inferred metallicity evolution of $α_{\rm CO}$ that is consistent with $α_{\rm [CI]}$ and previous estimates from the literature. These results thus support the use of the absorption-derived $α_{\rm [CI]}$ conversion factor for emission-selected star-forming galaxies and demonstrate that both methods probe the same universal properties of molecular gas in the local and high-redshift universe.

preprint2017arXiv

A merger in the dusty, $z=7.5$ galaxy A1689-zD1?

The gravitationally-lensed galaxy A1689-zD1 is one of the most distant spectroscopically confirmed sources ($z=7.5$). It is the earliest known galaxy where the interstellar medium (ISM) has been detected; dust emission was detected with the Atacama Large Millimetre Array (ALMA). A1689-zD1 is also unusual among high-redshift dust emitters as it is a sub-L* galaxy and is therefore a good prospect for the detection of gaseous ISM in a more typical galaxy at this redshift. We observed A1689-zD1 with ALMA in bands 6 and 7 and with the Green Bank Telescope (GBT) in band $Q$. To study the structure of A1689-zD1, we map the mm thermal dust emission and find two spatial components with sizes about $0.4-1.7$\,kpc (lensing-corrected). The rough spatial morphology is similar to what is observed in the near-infrared with {\it HST} and points to a perturbed dynamical state, perhaps indicative of a major merger or a disc in early formation. The ALMA photometry is used to constrain the far-infrared spectral energy distribution, yielding a dust temperature ($T_{\rm dust} \sim 35$--$45$\,K for $β= 1.5-2$). We do not detect the CO(3-2) line in the GBT data with a 95\% upper limit of 0.3\,mJy observed. We find a slight excess emission in ALMA band~6 at 220.9\,GHz. If this excess is real, it is likely due to emission from the [CII] 158.8\,$μ$m line at $z_{\rm [CII]} = 7.603$. The stringent upper limits on the [CII]/$L_{\rm FIR}$ luminosity ratio suggest a [CII] deficit similar to several bright quasars and massive starbursts.

preprint2015arXiv

A dusty, normal galaxy in the epoch of reionization

Candidates for the modest galaxies that formed most of the stars in the early universe, at redshifts $z > 7$, have been found in large numbers with extremely deep restframe-UV imaging. But it has proved difficult for existing spectrographs to characterise them in the UV. The detailed properties of these galaxies could be measured from dust and cool gas emission at far-infrared wavelengths if the galaxies have become sufficiently enriched in dust and metals. So far, however, the most distant UV-selected galaxy detected in dust emission is only at $z = 3.25$, and recent results have cast doubt on whether dust and molecules can be found in typical galaxies at this early epoch. Here we report thermal dust emission from an archetypal early universe star-forming galaxy, A1689-zD1. We detect its stellar continuum in spectroscopy and determine its redshift to be $z = 7.5\pm0.2$ from a spectroscopic detection of the Lyα break. A1689-zD1 is representative of the star-forming population during reionisation, with a total star-formation rate of about 12M$_\odot$ yr$^{-1}$. The galaxy is highly evolved: it has a large stellar mass, and is heavily enriched in dust, with a dust-to-gas ratio close to that of the Milky Way. Dusty, evolved galaxies are thus present among the fainter star-forming population at $z > 7$, in spite of the very short time since they first appeared.

preprint2015arXiv

Extinction curve template for intrinsically reddened quasars

We analyze the near-infrared to UV data of 16 quasars with redshifts ranging from 0.71 $<$ $z$ $<$ 2.13 to investigate dust extinction properties. The sample presented in this work is obtained from the High $A_V$ Quasar (HAQ) survey. The quasar candidates were selected from the Sloan Digital Sky Survey (SDSS) and the UKIRT Infrared Deep Sky Survey (UKIDSS), and follow-up spectroscopy was carried out at the Nordic Optical Telescope (NOT) and the New Technology Telescope (NTT). To study dust extinction curves intrinsic to the quasars, from the HAQ survey we selected 16 cases where the Small Magellanic Cloud (SMC) law could not provide a good solution to the spectral energy distributions (SEDs). We derived the extinction curves using Fitzpatrick & Massa 1986 (FM) law by comparing the observed SEDs to the combined quasar template from Vanden Berk et al. 2001 and Glikman et al. 2006. The derived extinction, $A_V$, ranges from 0.2-1.0 mag. All the individual extinction curves of our quasars are steeper ($R_V=2.2$-2.7) than that of the SMC, with a weighted mean value of $R_V=2.4$. We derive an `average quasar extinction curve' for our sample by fitting SEDs simultaneously by using the weighted mean values of the FM law parameters and a varying $R_V$. The entire sample is well fit with a single best-fit value of $R_V=2.2\pm0.2$. The `average quasar extinction curve' deviates from the steepest Milky Way and SMC extinction curves at a confidence level $\gtrsim95\%$. Such steep extinction curves suggest a significant population of silicates to produce small dust grains. Moreover, another possibility could be that the large dust grains may have been destroyed by the activity of the nearby active galactic nuclei (AGN), resulting in steep extinction curves.

preprint2015arXiv

Simulations of the OzDES AGN Reverberation Mapping Project

As part of the OzDES spectroscopic survey we are carrying out a large scale reverberation mapping study of $\sim$500 quasars over five years in the 30 deg$^2$ area of the Dark Energy Survey (DES) supernova fields. These quasars have redshifts ranging up to 4 and have apparent AB magnitudes between $16.8<r<22.5$ mag. The aim of the survey is to measure time lags between fluctuations in the quasar continuum and broad emission line fluxes of individual objects in order to measure black hole masses for a broad range of AGN and constrain the radius-luminosity ($R-L$) relationship. Here we investigate the expected efficiency of the OzDES reverberation mapping campaign and its possible extensions. We expect to recover lags for $\sim$35-45\% of the quasars. AGN with shorter lags and greater variability are more likely to yield a lag, and objects with lags $\lesssim$6 months or $\sim$1 year are expected be recovered the most accurately. The baseline OzDES reverberation mapping campaign is predicted to produce an unbiased measurement of the $R-L$ relationship parameters for H$β$, Mg II $λ$2798, and C IV $λ$1549. However, extending the baseline survey by either increasing the spectroscopic cadence, extending the survey season, or improving the emission line flux measurement accuracy will significantly improve the $R-L$ parameter constraints for all broad emission lines.

preprint2015arXiv

Using machine learning to classify the diffuse interstellar bands

Using over a million and a half extragalactic spectra we study the correlations of the Diffuse Interstellar Bands (DIBs) in the Milky Way. We measure the correlation between DIB strength and dust extinction for 142 DIBs using 24 stacked spectra in the reddening range E(B-V) < 0.2, many more lines than ever studied before. Most of the DIBs do not correlate with dust extinction. However, we find 10 weak and barely studied DIBs with correlations that are higher than 0.7 with dust extinction and confirm the high correlation of additional 5 strong DIBs. Furthermore, we find a pair of DIBs, 5925.9A and 5927.5A which exhibits significant negative correlation with dust extinction, indicating that their carrier may be depleted on dust. We use Machine Learning algorithms to divide the DIBs to spectroscopic families based on 250 stacked spectra. By removing the dust dependency we study how DIBs follow their local environment. We thus obtain 6 groups of weak DIBs, 4 of which are tightly associated with C2 or CN absorption lines.

preprint2014arXiv

A dust-parallax distance of 19 megaparsecs to the supermassive black hole in NGC 4151

The active galaxy NGC 4151 has a crucial role as one of only two active galactic nuclei for which black hole mass measurements based on emission line reverberation mapping can be calibrated against other dynamical methods. Unfortunately, effective calibration requires an accurate distance to NGC 4151, which is currently not available. Recently reported distances range from 4 to 29 megaparsecs (Mpc). Strong peculiar motions make a redshift-based distance very uncertain, and the geometry of the galaxy and its nucleus prohibit accurate measurements using other techniques. Here we report a dust-parallax distance to NGC 4151 of $D_A = 19.0^{+2.4}_{-2.6}$ Mpc. The measurement is based on an adaptation of a geometric method proposed previously using the emission line regions of active galaxies. Since this region is too small for current imaging capabilities, we use instead the ratio of the physical-to-angular sizes of the more extended hot dust emission as determined from time-delays and infrared interferometry. This new distance leads to an approximately 1.4-fold increase in the dynamical black hole mass, implying a corresponding correction to emission line reverberation masses of black holes if they are calibrated against the two objects with additional dynamical masses.

preprint2014arXiv

Dusting off the diffuse interstellar bands: DIBs and dust in extragalactic SDSS spectra

Using over a million and a half extragalactic spectra we study the properties of the mysterious Diffuse Interstellar Bands (DIBs) in the Milky Way. These data provide us with an unprecedented sampling of the skies at high Galactic-latitude and low dust-column-density. We present our method, study the correlation of the equivalent width of 8 DIBs with dust extinction and with a few atomic species, and the distribution of four DIBs - 5780.6A, 5797.1A, 6204.3A, and 6613.6A - over nearly 15000 squared degrees. As previously found, DIBs strengths correlate with extinction and therefore inevitably with each other. However, we show that DIBs can exist even in dust free areas. Furthermore, we find that the DIBs correlation with dust varies significantly over the sky. DIB under- or over-densities, relative to the expectation from dust, are often spread over hundreds of square degrees. These patches are different for the four DIBs, showing that they are unlikely to originate from the same carrier, as previously suggested.

preprint2014arXiv

Investigating CXOU J163802.6-471358: a new pulsar wind nebula in the Norma region?

We present the first analysis of the extended source CXOU J163802.6--471358, which was discovered serendipitously during the {\em Chandra} X-ray survey of the Norma region of the Galactic spiral arms. The X-ray source exhibits a cometary appearance with a point source and an extended tail region. The complete source spectrum is fitted well with an absorbed power law model and jointly fitting the {\em Chandra} spectrum of the full source with one obtained from an archived {\em XMM-Newton} observation results in best fit parameters $N_{\rm H}$ $=1.5^{+0.7}_{-0.5}\times10^{23} \text{cm}{^{-2}}$ and $Γ=1.1^{+0.7}_{-0.6}$ (90$%$ confidence uncertainties). The unabsorbed luminosity of the full source is then $L_X\sim 4.8\times10^{33}d_{10}^2$ergs s$^{-1}$ with $d_{10}=d/10$kpc, where a distance of 10 kpc is a lower bound inferred from the large column density. The radio counterpart found for the source using data from the Molonglo Galactic Plane Survey epoch-2 (MGPS-2) shows an elongated tail offset from the X-ray emission. No infrared counterpart was found. The results are consistent with the source being a previously unknown pulsar driving a bow shock through the ambient medium.

preprint2014arXiv

Rapid formation of large dust grains in the luminous supernova SN 2010jl

The origin of dust in galaxies is still a mystery. The majority of the refractory elements are produced in supernova explosions but it is unclear how and where dust grains condense and grow, and how they avoid destruction in the harsh environments of star-forming galaxies. The recent detection of 0.1-0.5 solar masses of dust in nearby supernova remnants suggests in situ dust formation, while other observations reveal very little dust in supernovae the first few years after explosion. Observations of the bright SN 2010jl have been interpreted as pre-existing dust, dust formation or no dust at all. Here we report the rapid (40-240 days) formation of dust in its dense circumstellar medium. The wavelength dependent extinction of this dust reveals the presence of very large (> 1 micron) grains, which are resistant to destructive processes. At later times (500-900 days), the near-IR thermal emission shows an accelerated growth in dust mass, marking the transition of the supernova from a circumstellar- to an ejecta-dominated source of dust. This provides the link between the early and late dust mass evolution in supernovae with dense circumstellar media.

preprint2014arXiv

Ultraviolet Emission Lines in Young Low Mass Galaxies at z~2: Physical Properties and Implications for Studies at z>7

We present deep spectroscopy of 17 very low mass (M* ~ 2.0x10^6 Msun to 1.4x10^9 Msun) and low luminosity (M_UV ~ -13.7 to -19.9) gravitationally lensed galaxies in the redshift range z~1.5-3.0. Deep rest-frame ultraviolet spectra reveal large equivalent width emission from numerous lines (NIV], OIII], CIV, Si III], CIII]) which are rarely seen in individual spectra of more massive star forming galaxies. CIII] is detected in 16 of 17 low mass star forming systems with rest-frame equivalent widths as large as 13.5 Angstroms. Nebular CIV emission is present in the most extreme CIII] emitters, requiring an ionizing source capable of producing a substantial component of photons with energies in excess of 47.9 eV. Photoionization models support a picture whereby the large equivalent widths are driven by the increased electron temperature and enhanced ionizing output arising from metal poor gas and stars, young stellar populations, and large ionization parameters. The young ages implied by the emission lines and continuum SEDs indicate that the extreme line emitters in our sample are in the midst of a significant upturn in their star formation activity. The low stellar masses, blue UV colors, and large sSFRs of our sample are similar to those of typical z>6 galaxies. Given the strong attenuation of Ly-alpha in z>6 galaxies we suggest that CIII] is likely to provide our best probe of early star forming galaxies with ground-based spectrographs and one of the most efficient means of confirming z>10 galaxies with the James Webb Space Telescope.

preprint2013arXiv

On inferring extinction laws in z~6 quasars as signatures of supernova dust

Unusual extinction curves of high-redshift QSOs have been taken as evidence that dust is primarily produced by supernovae at high redshift. In particular, the 3000 A Todini-Ferrara-Maiolino kink in the extinction curve of the z = 6.20 SDSS J1048+4637 has been attributed to supernova dust. Here we discuss the challenges in inferring robust extinction curves of high-redshift QSOs and critically assess previous claims of detection of supernova dust. In particular, we address the sensitivity to the choice of intrinsic QSO spectrum, the need for a long wavelength baseline, and the drawbacks in fitting theoretical extinction curves. In a sample of 21 QSOs at z ~ 6 we detect significant ultraviolet extinction using existing broad-band optical, near-infrared, and Spitzer photometry. The median extinction curve is consistent with a Small Magellanic Cloud curve with A_1450 ~ 0.7 mag and does not exhibit any conspicuous (restframe) 2175 A or 3000 A features. For two QSOs, SDSS J1044-0125 at z = 5.78 and SDSS J1030+0524 at z = 6.31, we further present X-shooter spectra covering the wavelength range 0.9-2.5 um. The resulting non-parametric extinction curves do not exhibit the 3000 A kink. Finally, in a re-analysis of literature spectra of SDSS J1048+4637, we do not find evidence for a conspicuous kink. We conclude that the existing evidence for a 3000 A feature is weak and that the overall dust properties at high and low redshift show no significant differences. This, however, does not preclude supernovae from dominating the dust budget at high redshift.

preprint2013arXiv

The metals-to-dust ratio to very low metallicities using GRB and QSO absorbers; extremely rapid dust formation

Among the key parameters defining the ISM of galaxies is the fraction of the metals that are locked up in dust: the metals-to-dust ratio. This ratio bears not only on the ISM and its evolution, but particularly on the origin of cosmic dust. We combine extinction and abundance data from GRB afterglows, from QSO absorbers, as well as from galaxy-lensed QSOs, to determine the metals-to-dust ratios for lines-of-sight through a wide diversity of galaxies from blue, dwarf starbursts to massive ellipticals, across a vast range in redshift z=0.1-6.3, and nearly three orders of magnitude in column density and metal abundance. We thus determine the metals-to-dust ratio in a unique way, providing direct determinations of in situ gas and dust columns without recourse to assumptions with large uncertainties. We find that the metals-to-dust ratios in these systems are surprisingly close to the value for the local group (10^{21.3} cm-2 A_V mag-1), with a mean value of 10^{21.2} cm-2 A_V mag-1 and a standard deviation of 0.3 dex. There is no evidence of deviation from this mean ratio as a function of metallicity, even down to our lowest metallicity of 0.01 Z/Z_sun. The lack of any obvious dependence of the metals-to-dust ratio on either column density, galaxy type or age, redshift, or metallicity indicates a close correspondence between the formation of the metals and the formation of dust. Any delay between the formation of metals and dust must be shorter than the typical metal-enrichment times of these galaxies. Formation of the bulk of the dust in low mass stars is therefore ruled out by these data at any cosmic epoch. Furthermore, dust destruction must not dominate over formation/growth in virtually any galaxy environment. The correlation between metals and dust is a natural consequence of the formation of the bulk of dust in SNe [Abridged].

preprint2013arXiv

Thermal emission in the early X-ray afterglows of GRBs: following the prompt phase to the late times

Thermal radiation, peaking in soft X-rays, has now been detected in a handful of GRB afterglows and has to date been interpreted as shock break-out of the GRB's progenitor star. We present a search for thermal emission in the early X-ray afterglows of a sample of Swift bursts selected by their brightness in X-rays at early times. We identify a clear thermal component in eight GRBs and track the evolution. We show that at least some of the emission must come from highly relativistic material since two show an apparent super-luminal expansion of the thermal component. Furthermore we determine very large luminosities and high temperatures for many of the components-too high to originate in a SN shock break-out. Instead we suggest that the component may be modelled as late photospheric emission from the jet, linking it to the apparently thermal component observed in the prompt emission of some GRBs at gamma-ray and hard X-ray energies. By comparing the parameters from the prompt emission and the early afterglow emission we find that the results are compatible with the interpretation that we are observing the prompt quasi-thermal emission component in soft X-rays at a later point in its evolution.

preprint2012arXiv

Dust extinction bias in the column density distribution of gamma-ray bursts; high column density, low redshift GRBs are more heavily obscured

The afterglows of gamma-ray bursts (GRBs) have more soft X-ray absorption than expected from the foreground gas column in the Galaxy. While the redshift of the absorption can in general not be constrained from current X-ray observations, it has been assumed that the absorption is due to metals in the host galaxy of the GRB. The large sample of X-ray afterglows and redshifts now available allows the construction of statistically meaningful distributions of the metal column densities. We construct such a sample and show, as found in previous studies, that the typical absorbing column density (N_HX) increases substantially with redshift, with few high column density objects found at low to moderate redshifts. We show, however, that when highly extinguished bursts are included in the sample, using redshifts from their host galaxies, high column density sources are also found at low to moderate redshift. We infer from individual objects in the sample and from observations of blazars, that the increase in column density with redshift is unlikely to be related to metals in the intergalactic medium or intervening absorbers. Instead we show that the origin of the apparent increase with redshift is primarily due to dust extinction bias: GRBs with high X-ray absorption column densities found at $z\lesssim4$ typically have very high dust extinction column densities, while those found at the highest redshifts do not. It is unclear how such a strongly evolving N_HX/A_V ratio would arise, and based on current data, remains a puzzle.

preprint2012arXiv

Helium in natal HII regions: the origin of the X-ray absorption in gamma-ray burst afterglows

Soft X-ray absorption in excess of Galactic is observed in the afterglows of most gamma-ray bursts (GRBs), but the correct solution to its origin has not been arrived at after more than a decade of work, preventing its use as a powerful diagnostic tool. We resolve this long-standing problem and find that He in the GRB's host HII region is responsible for most of the absorption. We show that the X-ray absorbing column density (N_Hx) is correlated with both the neutral gas column density and with the optical afterglow extinction (Av). This correlation explains the connection between dark bursts and bursts with high N_Hx values. From these correlations we exclude an origin of the X-ray absorption which is not related to the host galaxy, i.e. the intergalactic medium or intervening absorbers are not responsible. We find that the correlation with the dust column has a strong redshift evolution, whereas the correlation with the neutral gas does not. From this we conclude that the column density of the X-ray absorption is correlated with the total gas column density in the host galaxy rather than the metal column density, in spite of the fact that X-ray absorption is typically dominated by metals. The strong redshift evolution of N_Hx/Av is thus a reflection of the cosmic metallicity evolution of star-forming galaxies. We conclude that the absorption of X-rays in GRB afterglows is caused by He in the HII region hosting the GRB. While dust is destroyed and metals are stripped of all of their electrons by the GRB to great distances, the abundance of He saturates the He-ionising UV continuum much closer to the GRB, allowing it to remain in the neutral or singly-ionised state. Helium X-ray absorption explains the correlation with total gas, the lack of strong evolution with redshift as well as the absence of dust, metal or hydrogen absorption features in the optical-UV spectra.

preprint2012arXiv

The optically unbiased GRB host (TOUGH) survey. I. Survey design and catalogs

Long-duration gamma-ray bursts (GRBs) are powerful tracers of star-forming galaxies. We have defined a homogeneous subsample of 69 Swift GRB-selected galaxies spanning a very wide redshift range. Special attention has been devoted to making the sample optically unbiased through simple and well-defined selection criteria based on the high-energy properties of the bursts and their positions on the sky. Thanks to our extensive follow-up observations, this sample has now achieved a comparatively high degree of redshift completeness, and thus provides a legacy sample, useful for statistical studies of GRBs and their host galaxies. In this paper we present the survey design and summarize the results of our observing program conducted at the ESO Very Large Telescope (VLT) aimed at obtaining the most basic properties of galaxies in this sample, including a catalog of R and Ks magnitudes and redshifts. We detect the host galaxies for 80 % of the GRBs in the sample, although only 42 % have Ks-band detections, which confirms that GRB-selected host galaxies are generally blue. The sample is not uniformly blue, however, with two extremely red objects detected. Moreover, galaxies hosting GRBs with no optical/NIR afterglows, whose identification therefore relies on X-ray localizations, are significantly brighter and redder than those with an optical/NIR afterglow. Our spectroscopic campaign has resulted in 77 % now having redshift measurements, with a median redshift of 2.14 +- 0.18. TOUGH alone includes 17 detected z > 2 Swift GRB host galaxies suitable for individual and statistical studies. Seven hosts have detections of the Ly-alpha emission line and we can exclude an early indication that Ly-alpha emission is ubiquitous among GRB hosts, but confirm that Ly-alpha is stronger in GRB-selected galaxies than in flux-limited samples of Lyman break galaxies.

preprint2012arXiv

The properties of the 2175AA extinction feature discovered in GRB afterglows

The unequivocal, spectroscopic detection of the 2175 bump in extinction curves outside the Local Group is rare. To date, the properties of the bump have been examined in only two GRB afterglows (GRB 070802 and GRB 080607). In this work we analyse in detail the detections of the 2175 extinction bump in the optical spectra of the two further GRB afterglows: GRB 080605 and 080805. We gather all available optical/NIR photometric, spectroscopic and X-ray data to construct multi-epoch SEDs for both GRB afterglows. We fit the SEDs with the Fitzpatrick & Massa (1990) model with a single or broken PL. We also fit a sample of 38 GRB afterglows, known to prefer a SMC-type extinction curve, with the same model. We find that the SEDs of GRB 080605 and GRB 080805 at two epochs are fit well with a single PL with a derived extinction of A_V = 0.52(+0.13 -0.16) and 0.50 (+0.13 -0.10), and 2.1(+0.7-0.6) and 1.5+/-0.2 respectively. While the slope of the extinction curve of GRB 080805 is not well-constrained, the extinction curve of GRB 080605 has an unusual very steep far-UV rise together with the 2175 bump. Such an extinction curve has previously been found in only a small handful of sightlines in the MW. One possible explanation of such an extinction curve may be dust arising from two different regions with two separate grain populations, however we cannot distinguish the origin of the curve. We finally compare the four 2175 bump sightlines to the larger GRB afterglow sample and to Local Group sightlines. We find that while the width and central positions of the bumps are consistent with what is observed in the Local Group, the relative strength of the detected bump (A_bump) for GRB afterglows is weaker for a given A_V than for almost any Local Group sightline. Such dilution of the bump strength may offer tentative support to a dual dust-population scenario.

preprint2012arXiv

The stellar masses and specific star-formation rates of submillimetre galaxies

Establishing the stellar masses (M*), and hence specific star-formation rates (sSFRs) of submillimetre galaxies (SMGs) is crucial for determining their role in the cosmic galaxy/star formation. However, there is as yet no consensus over the typical M* of SMGs. Specifically, even for the same set of SMGs, the reported average M* have ranged over an order of magnitude, from ~5x10^10 Mo to ~5x10^11 Mo. Here we study how different methods of analysis can lead to such widely varying results. We find that, contrary to recent claims in the literature, potential contamination of IRAC 3-8 um photometry from hot dust associated with an active nucleus is not the origin of the published discrepancies in derived M*. Instead, we expose in detail how inferred M* depends on assumptions made in the photometric fitting, and quantify the individual and cumulative effects of different choices of initial mass function, different brands of evolutionary synthesis models, and different forms of assumed star-formation history. We review current observational evidence for and against these alternatives as well as clues from the hydrodynamical simulations, and conclude that, for the most justifiable choices of these model inputs, the average M* of SMGs is ~2x10^11 Mo. We also confirm that this number is perfectly reasonable in the light of the latest measurements of their dynamical masses, and the evolving M* function of the overall galaxy population. M* of this order imply that the average sSFR of SMGs is comparable to that of other star-forming galaxies at z>2, at 2-3 Gyr^-1. This supports the view that, while rare outliers may be found at any M*, most SMGs simply form the top end of the main-sequence of star-forming galaxies at these redshifts. Conversely, this argues strongly against the viewpoint that SMGs are extreme pathological objects, of little relevance in the cosmic history of star-formation.

preprint2012arXiv

TOUGH: Observational aspects of gamma-ray burst host galaxies

GRB-selected galaxies are broadly known to be faint, blue, young, star-forming dwarf galaxies. This insight, however, is based in part on heterogeneous samples of optically selected, lower-redshift galaxies. To study the statistical properties of GRB-selected galaxies we here introduce The Optically Unbiased GRB Host (TOUGH) complete sample of 69 X-ray selected Swift GRB host galaxies spanning the redshift range 0.03-6.30 and summarise the first results of a large observational survey of these galaxies.

preprint2011arXiv

Dust grain growth in the interstellar medium of galaxies at redshifts 4<z<6.5

To discriminate between different dust formation processes is a key issue in order to understand its properties. We analysed six submillimeter galaxies at redshifts 4<z<5 and nine quasars at 5<z<6.4. We estimated their dust masses from their (sub)millimeter emission and their stellar masses from the spectral energy distribution modelling or from the dynamical and gas masses obtained from the CO line detections. We calculated the dust yields per AGB star and per SN required to explain these dust masses and concluded that AGB stars are not efficient enough to form dust in the majority of these galaxies. SN could be responsible for dust production, but only if dust destruction in the SN shocks is not taken into account. Otherwise even SNe are not efficient enough, which advocates for some other dust production mechanism. We present the hypothesis that grain growth in the interstellar medium is responsible for bulk of the dust mass accumulation in these galaxies.

preprint2011arXiv

The extinction curves of star-forming regions from z=0.1 to 6.7 using GRB afterglow spectroscopy

GRB afterglows are well suited to extinction studies due to their brightness, simple power-law spectra and the occurrence of GRBs in distant star forming galaxies. In this paper we present results from the SED analysis of a sample of 41 GRB afterglows, from X-ray to NIR wavelengths. This is the largest sample of extinction curves outside the Local Group and, to date, the only extragalactic sample of absolute extinction curves based on spectroscopy. Visual extinction correlation with HI column density as well as total and gas-phase metal column density are examined. Approximately half the sample require a cooling break between the optical and X-ray regimes. The broken power-law SEDs show an average change in the spectral index of delta_beta=0.51 with a standard deviation of 0.02. This is consistent with the expectation from a simple synchrotron model. Of the sample, 63% are well described by the SMC-type extinction curve and have moderate or low extinction, with AV<0.65. Almost a quarter of our sample is consistent with no significant extinction (typically AV<0.1). The 2175AA extinction bump is detected unequivocally in 7% of our sample (3 GRBs), which all have A_V>1.0. We find an anti-correlation between gas-to-dust ratio and metallicity consistent with the Local Group relation. Our metals-to-dust ratios derived from the soft X-ray absorption are always larger (3-30 times) than the Local Group value, which may mean that GRB hosts may be less efficient at turning their metals into dust. However, we find that gas, dust, and metal column densities are all likely to be influenced by photo-ionization and dust destruction effects from the GRB. [abridged]

preprint2011arXiv

The Galactic dust-to-metals ratio and metallicity using gamma-ray bursts

The metallicity and dust-to-metals ratio of the Galaxy are fundamental parameters in understanding the ISM, but there is still uncertainty surrounding these parameters. In this paper, the dust-to-metals ratio in the Galaxy is determined using the photoelectric absorption of the X-ray afterglows of a sample of several hundred gamma-ray bursts (GRBs) to determine the metal column density in combination with Galactic dust maps to determine the line-of-sight dust extinction through the Galaxy in the direction of the GRB. GRB afterglows often have large extragalactic soft X-ray absorptions and therefore the GRB sample's upper-bound will define the Galactic dust-to-metals relation. Using a two-dimensional two-sample KS test, we determine this upper-bound and so derive the dust-to-metals ratio of the Galaxy. We find N_H = 2.2^{+0.3}_{-0.4}e21 cm^-2 A_V assuming solar, Anders & Grevesse (1989), metallicity. This result is consistent with previous findings using bright X-ray sources in the Galaxy. Using the same technique but substituting the HI maps from the Leiden-Argentine-Bonn survey for the dust maps, allows us to place a limit on the metallicity in the Galaxy. We find a metallicity consistent with the Anders & Grevesse (1989) solar values often used in X-ray fitting. Based on this and previous studies, we suggest that the metallicity of a typical ISM sightline through the Galaxy is ~0.25 dex higher than the current best estimate of the solar metallicity. We further show that the dust-to-gas ratio seems to be correlated with the total gas column density, and that this may be due to the metallicity gradient observed toward the Galactic centre. Based on the non-constant nature of the dust-to-gas ratio, we propose that the dust column density, at N_H = 2.2e21 cm^-2 A_V, represents a better proxy for the soft X-ray absorption column density than HI maps.

preprint2010arXiv

Cosmic evolution of submillimeter galaxies and their contribution to stellar mass assembly

The nature of galaxies selected at submillimeter wavelengths (SMGs, S_850 > 3 mJy), some of the bolometrically most luminous objects at high redshifts, is still elusive. In particular their star formation histories and source of emission are not accurately constrained. In this paper we introduce a new approach to analyse the SMG data. Namely, we present the first self-consistent UV-to-radio spectral energy distribution fits of 76 SMGs with spectroscopic redshifts using all photometric datapoints from ultraviolet to radio simultaneously. We find that they are highly star-forming (median star formation rate 713 MSun yr^-1 for SMGs at z>0.5), moderately dust-obscured (median A_V~2 mag), hosting significant stellar populations (median stellar mass 3.7x10^11 MSun) of which only a minor part has been formed in the ongoing starburst episode. This implies that in the past, SMGs experienced either another starburst episode or merger with several galaxies. The properties of SMGs suggest that they are progenitors of present-day elliptical galaxies. We find that these bright SMGs contribute significantly to the cosmic star formation rate density (~20%) and stellar mass density (~30-50%) at redshifts 2-4. Using number counts at low fluxes we find that as much as 80% of the cosmic star formation at these redshifts took place in SMGs brighter than 0.1 mJy. We find evidence that a linear infrared-radio correlation holds for SMGs in an unchanged form up to redshift of 3.6, though its normalization is offset from the local relation by a factor of ~2.1 towards higher radio luminosities. We present a compilation of photometry data of SMGs and determinations of cosmic SFR and stellar mass densities.

preprint2010arXiv

Dust grain growth in the interstellar medium of 5<z<6.5 quasars

We investigate whether stellar dust sources i.e. asymptotic giant branch (AGB) stars and supernovae (SNe) can account for dust detected in 5<z<6.5 quasars (QSOs). We calculate the required dust yields per AGB star and per SN using the dust masses of QSOs inferred from their millimeter emission and stellar masses approximated as the difference between the dynamical and the H_2 gas masses of these objects. We find that AGB stars are not efficient enough to form dust in the majority of the z>5 QSOs, whereas SNe may be able to account for dust in some QSOs. However, they require very high dust yields even for a top-heavy initial mass function. This suggests additional non-stellar dust formation mechanism e.g. significant dust grain growth in the interstellar medium of at least three out of nine z>5 QSOs. SNe (but not AGB stars) may deliver enough heavy elements to fuel this growth.

preprint2010arXiv

Rapid dust production in submillimeter galaxies at z>4?

The existence of submillimeter-selected galaxies (SMGs) at redshifts z>4 has recently been confirmed. Using simultaneously all the available data from UV to radio we have modeled the spectral energy distributions of the six known spectroscopically confirmed SMGs at z>4. We find that their star formation rates (average ~2500 MSun yr^{-1}), stellar (~3.6x10^{11} MSun) and dust (~6.7x10^{8} MSun) masses, extinction (A_V~2.2 mag), and gas-to-dust ratios (~60) are within the ranges for 1.7<z<3.6 SMGs. Our analysis suggests that infrared-to-radio luminosity ratios of SMGs do not change up to redshift ~5 and are lower by a factor of ~2.1 than the value corresponding to the local IR-radio correlation. However, we also find dissimilarities between z>4 and lower-redshift SMGs. Those at z>4 tend to be among the most star-forming, least massive and hottest (~60 K) SMGs and exhibit the highest fraction of stellar mass formed in the ongoing starburst (~45%). This indicates that at z>4 we see earlier stages of evolution of submillimeter-bright galaxies. Using the derived properties for z>4 SMGs we investigate the origin of dust at epochs less than 1.5 Gyr after the big bang. This is significant to our understanding of the evolution of the early universe. For three z>4 SMGs, asymptotic giant branch stars could be the dominant dust producers. However, for the remaining three only supernovae (SNe) are efficient and fast enough to be responsible for dust production, though requiring a very high dust yield per SN (0.15-0.65 MSun). The required dust yields are lower if a top-heavy initial mass function or significant dust growth in the interstellar medium is assumed. We estimate lower limits of the contribution of SMGs to the cosmic star formation and stellar mass densities at z~4-5 to be ~4% and ~1%, respectively.

preprint2005arXiv

The optical afterglow of the short gamma-ray burst GRB 050709

It has long been known that there are two classes of gamma-ray bursts (GRBs), mainly distinguished by their durations. The breakthrough in our understanding of long-duration GRBs (those lasting more than ~2 s), which ultimately linked them with energetic Type Ic supernovae, came from the discovery of their long-lived X-ray and optical afterglows, when precise and rapid localizations of the sources could finally be obtained. X-ray localizations have recently become available for short (duration <2 s) GRBs, which have evaded optical detection for more than 30 years. Here we report the first discovery of transient optical emission (R-band magnitude ~23) associated with a short burst; GRB 050709. The optical afterglow was localized with subarcsecond accuracy, and lies in the outskirts of a blue dwarf galaxy. The optical and X-ray afterglow properties 34 h after the GRB are reminiscent of the afterglows of long GRBs, which are attributable to synchrotron emission from ultrarelativistic ejecta. We did not, however, detect a supernova, as found in most nearby long GRB afterglows, which suggests a different origin for the short GRBs.

preprint2003arXiv

A very energetic supernova associated with the gamma-ray burst of 29 March 2003

Over the past five years evidence has mounted that long-duration (> 2 s) gamma-ray bursts (GRBs)--the most brilliant of all astronomical explosions--signal the collapse of massive stars in our Universe. This evidence was originally based on the probable association of one unusual GRB with a supernova, but now includes the association of GRBs with regions of massive star formation in distant galaxies, the appearance of supernova-like 'bumps' in the optical afterglow light curves of several bursts and lines of freshly synthesized elements in the spectra of a few X-ray afterglows. These observations support, but do not yet conclusively demonstrate, the idea that long-duration GRBs are associated with the deaths of massive stars, presumably arising from core collapse. Here we report evidence that a very energetic supernova (a hypernova) was temporally and spatially coincident with a GRB at redshift z = 0.1685. The timing of the supernova indicates that it exploded within a few days of the GRB, strongly suggesting that core-collapse events can give rise to GRBs, thereby favouring the 'collapsar' model.