Source author record

A. P. Thomson

A. P. Thomson 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

21works
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

21 published item(s)

preprint2022arXiv

Deep sub-arcsecond widefield imaging of the Lockman Hole field at 144 MHz

High quality low-frequency radio surveys have the promise of advancing our understanding of many important topics in astrophysics, including the life cycle of active galactic nuclei (AGN), particle acceleration processes in jets, the history of star formation, and exoplanet magnetospheres. Currently leading low-frequency surveys reach an angular resolution of a few arcseconds. However, this resolution is not yet sufficient to study the more compact and distant sources in detail. Sub-arcsecond resolution is therefore the next milestone in advancing these fields. The biggest challenge at low radio frequencies is the ionosphere. If not adequately corrected for, ionospheric seeing blurs the images to arcsecond or even arcminute scales. Additionally, the required image size to map the degree-scale field of view of low-frequency radio telescopes at this resolution is far greater than what typical soft- and hardware is currently capable of handling. Here we present for the first time (to the best of our knowledge) widefield sub-arcsecond imaging at low radio frequencies. We derive ionospheric corrections in a few dozen individual directions and apply those during imaging efficiently using a recently developed imaging algorithm (arXiv:1407.1943, arXiv:1909.07226). We demonstrate our method by applying it to an eight hour observation of the International LOw Frequency ARray (LOFAR) Telescope (ILT) (arXiv:1305.3550). Doing so we have made a sensitive $7.4\ \mathrm{deg}^2$ $144\ \mathrm{MHz}$ map at a resolution of $0.3''$ reaching $25\ μ\mathrm{Jy\ beam}^{-1}$ near the phase centre. The estimated $250,000$ core hours used to produce this image, fit comfortably in the budget of available computing facilities. This result will enable future mapping of the entire northern low-frequency sky at sub-arcsecond resolution.

preprint2022arXiv

Quasar Feedback Survey: Multi-phase outflows, turbulence and evidence for feedback caused by low power radio jets inclined into the galaxy disk

We present a study of a luminous, z=0.15, type-2 quasar (log [L([OIII])/(erg/s)]=42.8) from the Quasar Feedback Survey. It is classified as 'radio-quiet' (log [L(1.4 GHz)/(W/Hz)]=23.8); however, radio imaging reveals ~1 kpc low-power jets (log [Pjet/(erg/s)]=44) inclined into the plane of the galaxy disk. We combine MUSE and ALMA observations to map stellar kinematics and ionised and molecular gas properties. The jets are seen to drive galaxy-wide bi-conical turbulent outflows, reaching W80 = 1000-1300 km/s, in the ionised phase (traced via optical emission-lines), which also have increased electron densities compared to the quiescent gas. The turbulent gas is driven perpendicular to the jet axis and is escaping along the galaxy minor axis, reaching 7.5 kpc on both sides. Traced via CO(3-2) emission, the turbulent material in molecular gas phase is one-third as spatially extended and has 3 times lower velocity-dispersion as compared to ionised gas. The jets are seen to be strongly interacting with the interstellar medium (ISM) through enhanced ionised emission and disturbed/depleted molecular gas at the jet termini. We see further evidence for jet-induced feedback through significantly higher stellar velocity-dispersion aligned, and co-spatial with, the jet axis (<5 deg). We discuss possible negative and positive feedback scenarios arising due to the interaction of the low-power jets with the ISM in the context of recent jet-ISM interaction simulations, which qualitatively agree with our observations. We discuss how jet-induced feedback could be an important feedback mechanism even in bolometrically luminous 'radio-quiet' quasars.

preprint2020arXiv

An ALMA survey of the S2CLS UDS field: Optically invisible submillimetre galaxies

We analyse a robust sample of 30 near-infrared-faint (K>25.3, 5 sigma) submillimetre galaxies selected across a 0.96 deg^2 field, to investigate their properties and the cause of their lack of detectable optical/near-infrared emission. Our analysis exploits precise identifications based on ALMA 870um continuum imaging, combined with the very deep near-infrared imaging from the UKIDSS-UDS survey. We estimate that K>25.3 submillimetre galaxies represent 15+/-2 per cent of the total population brighter than S870=3.6mJy, with an expected surface density of ~450/deg^2 above S870>1mJy. As such they pose a source of contamination in surveys for both high-redshift "quiescent" galaxies and very-high-redshift Lyman-break galaxies. We show that these K-faint submillimetre galaxies are simply the tail of the broader submillimetre population, with comparable dust and stellar masses to K<25.3 mag submillimetre galaxies, but lying at significantly higher redshifts (z=3.44+/-0.06 versus z=2.36+/-0.11) and having higher dust attenuation (Av=5.2+/-0.3 versus Av=2.9+/-0.1). We investigate the origin of the strong dust attenuation and find indications that these K-faint galaxies have smaller dust continuum sizes than the K<25.3 galaxies, as measured by ALMA, which suggests their high attenuation is related to their compact sizes. We find a correlation of dust attenuation with star-formation rate surface density (Sigma_SFR), with the K-faint submillimetre galaxies representing the higher-Sigma_SFR and highest-Av galaxies. The concentrated, intense star-formation activity in these systems is likely to be associated with the formation of spheroids in compact galaxies at high redshifts, but as a result of their high obscuration these are completely missed in UV, optical and even near-infrared surveys.

preprint2020arXiv

An ALMA survey of the SCUBA-2 Cosmology Legacy Survey UKIDSS/UDS field: Dust attenuation in high-redshift Lyman break Galaxies

We analyse 870um Atacama Large Millimetre Array (ALMA) dust continuum detections of 41 canonically-selected z~3 Lyman-break galaxies (LBGs), as well as 209 ALMA-undetected LBGs, in follow-up of SCUBA-2 mapping of the UKIDSS Ultra Deep Survey (UDS) field. We find that our ALMA-bright LBGs lie significantly off the locally calibrated IRX-beta relation and tend to have relatively bluer rest-frame UV slopes (as parametrised by beta), given their high values of the 'infrared excess' (IRX=L_IR/L_UV), relative to the average 'local' IRX-beta relation. We attribute this finding in part to the young ages of the underlying stellar populations but we find that the main reason behind the unusually blue UV slopes are the relatively shallow slopes of the corresponding dust attenuation curves. We show that, when stellar masses are being established via SED fitting, it is absolutely crucial to allow the attenuation curves to vary (rather than fixing it on Calzetti-like law), where we find that the inappropriate curves may underestimate the resulting stellar masses by a factor of ~2-3x on average. In addition, we find these LBGs to have relatively high specific star-formation rates (sSFRs), dominated by the dust component, as quantified via the fraction of obscured star formation ( f_obs = SFR_IR/SFR_(UV+IR)). We conclude that the ALMA-bright LBGs are, by selection, massive galaxies undergoing a burst of a star formation (large sSFRs, driven, for example, by secular or merger processes), with a likely geometrical disconnection of the dust and stars, responsible for producing shallow dust attenuation curves.

preprint2020arXiv

Cosmic evolution of molecular gas mass density from an empirical relation between $\rm L_{1.4GHz}$ and $\rm L^{\prime}_{CO}$

Historically, GHz radio emission has been used extensively to characterize the star-formation activity in galaxies. In this work, we look for empirical relations amongst the radio luminosity, the infrared luminosity, and the CO-based molecular gas mass. We assemble a sample of 278 nearby galaxies with measurements of radio continuum and total infrared emission, and the $^{12}$CO (J = 1-0) emission line. We find a correlation between the radio continuum and the CO emission line (with a scatter of 0.36 dex), in a large sample of different kind of galaxies. Making use of this correlation, we explore the evolution of the molecular gas mass function and the cosmological molecular gas mass density in six redshift bins up to $z = 1.5$. These results agree with previous semi-analytic predictions and direct measurements: the cosmic molecular gas density increases up to $z=1.5$. In addition, we find a single plane across five orders of magnitude for the explored luminosities, with a scatter of 0.27 dex. These correlations are sufficiently robust to be used for samples where no CO measurements exist.

preprint2020arXiv

Extended H$α$ over compact far-infrared continuum in dusty submillimeter galaxies -- Insights into dust distributions and star-formation rates at $z\sim2$

Using data from ALMA and near-infrared (NIR) integral field spectrographs including both SINFONI and KMOS on the VLT, we investigate the two-dimensional distributions of H$α$ and rest-frame far-infrared (FIR) continuum in six submillimeter galaxies at $z\sim2$. At a similar spatial resolution ($\sim$0.5" FWHM; $\sim$4.5 kpc at $z=2$), we find that the half-light radius of H$α$ is significantly larger than that of the FIR continuum in half of the sample, and on average H$α$ is a median factor of $2.0\pm0.4$ larger. Having explored various ways to correct for the attenuation, we find that the attenuation-corrected H$α$-based SFRs are systematically lower than the IR-based SFRs by at least a median factor of $3\pm1$, which cannot be explained by the difference in half-light radius alone. In addition, we find that in 40% of cases the total $V$-band attenuation ($A_V$) derived from energy balance modeling of the full ultraviolet(UV)-to-FIR spectral energy distributions (SEDs) is significantly higher than that derived from SED modeling using only the UV-to-NIR part of the SEDs, and the discrepancy appears to increase with increasing total infrared luminosity. Finally, considering all our findings along with the studies in the literature, we postulate that the dust distributions in SMGs, and possibly also in less IR luminous $z\sim2$ massive star-forming galaxies, can be decomposed into three main components; the diffuse dust heated by older stellar populations, the more obscured and extended young star-forming HII regions, and the heavily obscured central regions that have a low filling factor but dominate the infrared luminosity in which the majority of attenuation cannot be probed via UV-to-NIR emissions.

preprint2020arXiv

Searching for obscured AGN in z $\sim$ 2 submillimetre galaxies

Submillimetre-selected galaxies (SMGs) at high redshift ($z$ $\sim$ 2) are potential host galaxies of active galactic nuclei (AGN). If the local Universe is a good guide, $\sim$ 50$\%$ of the obscured AGN amongst the SMG population could be missed even in the deepest X-ray surveys. Radio observations are insensitive to obscuration; therefore, very long baseline interferometry (VLBI) can be used as a tool to identify AGN in obscured systems. A well-established upper limit to the brightness temperature of 10$^5$ K exists in star-forming systems, thus VLBI observations can distinguish AGN from star-forming systems via brightness temperature measurements. We present 1.6 GHz European VLBI Network (EVN) observations of four SMGs (with measured redshifts) to search for evidence of compact radio components associated with AGN cores. For two of the sources, e-MERLIN images are also presented. Out of the four SMGs observed, we detect one source, J123555.14, that has an integrated EVN flux density of 201 $\pm$ 15.2 $μ$Jy, corresponding to a brightness temperature of 5.2 $\pm$ 0.7 $\times$ 10$^5$ K. We therefore identify that the radio emission from J123555.14 is associated with an AGN. We do not detect compact radio emission from a possible AGN in the remaining sources (J123600.10, J131225.73, and J163650.43). In the case of J131225.73, this is particularly surprising, and the data suggest that this may be an extended, jet-dominated AGN that is resolved by VLBI. Since the morphology of the faint radio source population is still largely unknown at these scales, it is possible that with a $\sim$ 10 mas resolution, VLBI misses (or resolves) many radio AGN extended on kiloparsec scales.

preprint2020arXiv

SuperCLASS -- I. The Super CLuster Assisted Shear Survey: Project overview and Data Release 1

The SuperCLuster Assisted Shear Survey (SuperCLASS) is a legacy programme using the e-MERLIN interferometric array. The aim is to observe the sky at L-band (1.4 GHz) to a r.m.s. of 7 uJy per beam over an area of ~1 square degree centred on the Abell 981 supercluster. The main scientific objectives of the project are: (i) to detect the effects of weak lensing in the radio in preparation for similar measurements with the Square Kilometre Array (SKA); (ii) an extinction free census of star formation and AGN activity out to z~1. In this paper we give an overview of the project including the science goals and multi-wavelength coverage before presenting the first data release. We have analysed around 400 hours of e-MERLIN data allowing us to create a Data Release 1 (DR1) mosaic of ~0.26 square degrees to the full depth. These observations have been supplemented with complementary radio observations from the Karl G. Jansky Very Large Array (VLA) and optical/near infra-red observations taken with the Subaru, Canada-France-Hawaii and Spitzer Telescopes. The main data product is a catalogue of 887 sources detected by the VLA, of which 395 are detected by e-MERLIN and 197 of these are resolved. We have investigated the size, flux and spectral index properties of these sources finding them compatible with previous studies. Preliminary photometric redshifts, and an assessment of galaxy shapes measured in the radio data, combined with a radio-optical cross-correlation technique probing cosmic shear in a supercluster environment, are presented in companion papers.

preprint2020arXiv

SuperCLASS -- III. Weak lensing from radio and optical observations in Data Release 1

We describe the first results on weak gravitational lensing from the SuperCLASS survey: the first survey specifically designed to measure the weak lensing effect in radio-wavelength data, both alone and in cross-correlation with optical data. We analyse 1.53 square degrees of optical data from the Subaru telescope and 0.26 square degrees of radio data from the e-MERLIN and VLA telescopes (the DR1 data set). Using standard methodologies on the optical data only we make a significant (10 sigma) detection of the weak lensing signal (a shear power spectrum) due to the massive supercluster of galaxies in the targeted region. For the radio data we develop a new method to measure the shapes of galaxies from the interferometric data, and we construct a simulation pipeline to validate this method. We then apply this analysis to our radio observations, treating the e-MERLIN and VLA data independently. We achieve source densities of 0.5 per square arcmin in the VLA data and 0.06 per square arcmin in the e-MERLIN data, numbers which prove too small to allow a detection of a weak lensing signal in either the radio data alone or in cross-correlation with the optical data. Finally, we show preliminary results from a visibility-plane combination of the data from e-MERLIN and VLA which will be used for the forthcoming full SuperCLASS data release. This approach to data combination is expected to enhance both the number density of weak lensing sources available and the fidelity with which their shapes can be measured.

preprint2020arXiv

The e-MERLIN Galaxy Evolution Survey (e-MERGE): Overview and Survey Description

We present an overview and description of the eMERLIN Galaxy Evolution survey (eMERGE) Data Release 1 (DR1), a large program of high-resolution 1.5 GHz radio observations of the GOODS-N field comprising $\sim140$ hours of observations with eMERLIN and $\sim40$ hours with the Very Large Array (VLA). We combine the long baselines of eMERLIN (providing high angular resolution) with the relatively closely-packed antennas of the VLA (providing excellent surface brightness sensitivity) to produce a deep 1.5 GHz radio survey with the sensitivity ($\sim 1.5μ$Jy beam$^{-1}$), angular resolution ($0.2"$--$0.7"$) and field-of-view ($\sim15' \times 15'$) to detect and spatially resolve star-forming galaxies and AGN at $z\gtrsim 1$. The goal of eMERGE is to provide new constraints on the deep, sub-arcsecond radio sky which will be surveyed by SKA1-mid. In this initial publication, we discuss our data analysis techniques, including steps taken to model in-beam source variability over a $\sim20$ year baseline and the development of new point spread function/primary beam models to seamlessly merge eMERLIN and VLA data in the $uv$ plane. We present early science results, including measurements of the luminosities and/or linear sizes of $\sim500$ galaxes selected at 1.5 GHz. In combination with deep Hubble Space Telescope observations, we measure a mean radio-to-optical size ratio of $r_{\rm eMERGE}/r_{\rm HST}\sim1.02\pm0.03$, suggesting that in most high-redshift galaxies, the $\sim$GHz continuum emission traces the stellar light seen in optical imaging. This is the first in a series of papers which will explore the $\sim$kpc-scale radio properties of star-forming galaxies and AGN in the GOODS-N field observed by eMERGE DR1.

preprint2019arXiv

The East Asian Observatory SCUBA--2 survey of the COSMOS field: unveiling 1147 bright sub-millimeter sources across 2.6 square degrees

We present sensitive 850$μ$m imaging of the COSMOS field using 640hr of new and archival observations taken with SCUBA-2 at the East Asian Observatory's James Clerk Maxwell Telescope. The SCUBA-2 COSMOS survey (S2COSMOS) achieves a median noise level of $σ_{850μ{\mathrm{m}}}$=1.2mJy/beam over an area of 1.6 sq. degree (MAIN; HST/ACS footprint), and $σ_{850μ{\mathrm{m}}}$=1.7mJy/beam over an additional 1 sq. degree of supplementary (SUPP) coverage. We present a catalogue of 1020 and 127 sources detected at a significance level of >4$σ$ and >4.3$σ$ in the MAIN and SUPP regions, respectively, corresponding to a uniform 2% false-detection rate. We construct the single-dish 850$μ$m number counts at $S_{850}$>2mJy and show that these S2COSMOS counts are in agreement with previous single-dish surveys, demonstrating that degree-scale fields are sufficient to overcome the effects of cosmic variance in the $S_{850}$=2-10mJy population. To investigate the properties of the galaxies identified by S2COSMOS sources we measure the surface density of near-infrared-selected galaxies around their positions and identify an average excess of 2.0$\pm$0.2 galaxies within a 13$''$ radius (~100kpc at $z$~2). The bulk of these galaxies represent near-infrared-selected SMGs and/or spatially-correlated sources and lie at a median photometric redshift of $z$=2.0$\pm$0.1. Finally, we perform a stacking analysis at sub-millimeter and far-infrared wavelengths of stellar-mass-selected galaxies ($M_{\star}$=10$^{10}$-10$^{12}{\rm M_{\odot}}$) from $z$=0-4, obtaining high-significance detections at 850um in all subsets (SNR=4-30), and investigate the relation between far-infrared luminosity, stellar mass, and the peak wavelength of the dust SED. The publication of this survey adds a new deep, uniform sub-millimeter layer to the wavelength coverage of this well-studied COSMOS field.

preprint2016arXiv

The SCUBA-2 Cosmology Legacy Survey: 850um maps, catalogues and number counts

We present a catalogue of nearly 3,000 submillimetre sources detected at 850um over ~5 square degrees surveyed as part of the James Clerk Maxwell Telescope (JCMT) SCUBA-2 Cosmology Legacy Survey (S2CLS). This is the largest survey of its kind at 850um, probing a meaningful cosmic volume at the peak of star formation activity and increasing the sample size of submillimetre galaxies selected at 850um by an order of magnitude. We describe the wide 850um survey component of S2CLS, which covers the key extragalactic survey fields: UKIDSS-UDS, COSMOS, Akari-NEP, Extended Groth Strip, Lockman Hole North, SSA22 and GOODS-North. The average 1-sigma depth of S2CLS is 1.2 mJy/beam, approaching the SCUBA-2 850um confusion limit, which we determine to be ~0.8 mJy/beam. We measure the single dish 850um number counts to unprecedented accuracy, reducing the Poisson errors on the differential counts to approximately 4% at S_850~3mJy. With several independent fields, we investigate field-to-field variance, finding that the number counts on 0.5-1 degree scales are generally within 50% of the S2CLS mean for S_850>3mJy, with scatter consistent with the Poisson and estimated cosmic variance uncertainties, although there is a marginal (2-sigma) density enhancement in the GOODS-North field. The observed number counts are in reasonable agreement with recent phenomenological and semi-analytic models. Finally, the large solid angle of S2CLS allows us to measure the bright-end counts: at S_850>10mJy there are approximately ten sources per square degree, and we detect the distinctive up-turn in the number counts indicative of the detection of local sources of 850um emission and strongly lensed high-redshift galaxies. Here we describe the data collection and reduction procedures and present calibrated maps and a catalogue of sources; these are made publicly available.

preprint2015arXiv

Dusty starbursts and the formation of elliptical galaxies: A SCUBA-2 survey of a z=1.46 cluster

We report the results of a deep SCUBA-2 850- and 450-$μ$m survey for dust-obscured ultra-luminous infrared galaxies (U/LIRGs) in the field of the z=1.46 cluster XCS J2215.9-1738. We detect a striking overdensity of sub-millimeter sources coincident with the core of this cluster: $\sim 3-4 \times$ higher than expected in a blank field. We use the likely radio and mid-infrared counterparts to show that the bulk of these sub-millimeter sources have spectroscopic or photometric redshifts which place them in the cluster and that their multi-wavelength properties are consistent with this association. The average far-infrared luminosities of these galaxies are $(1.0\pm0.1) \times 10^{12} L_{\odot}$, placing them on the U/LIRG boundary. Using the total star formation occurring in the obscured U/LIRG population within the cluster we show that the resulting mass-normalized star-formation rate for this system supports previous claims of a rapid increase in star-formation activity in cluster cores out to $z\sim1.5$, which must be associated with the on-going formation of the early-type galaxies which reside in massive clusters today.

preprint2015arXiv

Storm in a "Teacup": a radio-quiet quasar with ~10kpc radio-emitting bubbles and extreme gas kinematics

We present multi-frequency (1-8 GHz) VLA data, combined with VIMOS IFU data and HST imaging, of a z=0.085 radio-quiet type 2 quasar (with L(1.4GHz)~5e23 W/Hz and L(AGN)~2e45 erg/s). Due to the morphology of its emission-line region, the target (J1430+1339) has been referred to as the Teacup AGN in the literature. We identify "bubbles" of radio emission that are extended ~10-12 kpc to both the east and west of the nucleus. The edge of the brighter eastern bubble is co-spatial with an arc of luminous ionized gas. We also show that the Teacup AGN hosts a compact radio structure, located ~0.8 kpc from the core position, at the base of the eastern bubble. This radio structure is co-spatial with an ionized outflow with an observed velocity of v=-740 km/s. This is likely to correspond to a jet, or possibly a quasar wind, interacting with the interstellar medium at this position. The large-scale radio bubbles appear to be inflated by the central AGN, which indicates that the AGN can also interact with the gas on >~10 kpc scales. Our study highlights that even when a quasar is formally "radio-quiet" the radio emission can be extremely effective for observing the effects of AGN feedback.

preprint2015arXiv

Tracing cool molecular gas and star formation on $\sim 100$pc scales within a $z=2.3$ galaxy

We present new, high-angular resolution interferometric observations with the Karl G. Jansky Very Large Array of $^{12}$CO $J=1-0$ line emission and 4-8 GHz continuum emission in the strongly lensed, $z=2.3$ submillimetre galaxy, SMM J21352-0102. Using these data, we identify and probe the conditions in $\sim 100$pc clumps within this galaxy, which we consider to be potential giant molecular cloud complexes, containing up to half of the total molecular gas in this system. In combination with far-infrared and submillimetre data, we investigate the far-infrared/radio correlation, measuring $q_{IR} = 2.39 \pm 0.17$ across SMM J21352. We search for variations in the properties of the interstellar medium throughout the galaxy by measuring the spatially-resolved $q_{IR}$ and radio spectral index, $α_{\rm radio}$, finding ranges $q_{IR} = [2.1, 2.6]$ and $α_{\rm radio} = [-1.5, -0.7]$. We argue that these ranges in $α_{\rm radio}$ and $q_{IR}$ may reflect variations in the age of the ISM material. Using multi-$J$ $^{12}$CO data, we quantitatively test a recent theoretical model relating the star-formation rate surface density to the excitation of $^{12}$CO, finding good agreement between the model and the data. Lastly, we study the Schmidt-Kennicutt relation, both integrated across the system and within the individual clumps. We find small offsets between SMM J21352 and its clumps relative to other star-forming galaxy populations on the Schmidt-Kennicutt plot - such offsets have previously been interpreted as evidence for a bi-modal star-formation law, but we argue that they can be equally-well explained as arising due to a combination of observational uncertainties and systematic biases in the choice of model used to interpret the data.

preprint2014arXiv

An ALMA survey of submillimetre galaxies in the Extended Chandra Deep Field South: radio properties and the far-infrared/radio correlation

We present a study of the radio properties of 870$μ$m-selected submillimetre galaxies (SMGs), observed at high resolution with ALMA in the Extended Chandra Deep Field South. From our initial sample of 76 ALMA SMGs, we detect 52 SMGs at $>3σ$ significance in VLA 1400MHz imaging, of which 35 are also detected at $>3σ$ in new 610MHz GMRT imaging. Within this sample of radio-detected SMGs, we measure a median radio spectral index $α_{610}^{1400} = -0.79 \pm 0.06$, (with inter-quartile range $α=[-1.16,-0.56]$) and investigate the far-infrared/radio correlation via the parameter $q_{\rm IR}$, the logarithmic ratio of the rest-frame 8-1000$μ$m flux and monochromatic radio flux. Our median $q_{\rm IR} = 2.56 \pm 0.05$ (inter-quartile range $q_{\rm IR}=[2.42,2.78]$) is higher than that typically seen in single-dish 870$μ$m-selected sources ($q_{\rm IR} \sim 2.4$), which may reflect the fact that our ALMA-based study is not biased to radio-bright counterparts, as previous samples were. Finally, we search for evidence that $q_{\rm IR}$ and $α$ evolve with age in a co-dependent manner, as predicted by starburst models: the data populate the predicted region of parameter space, with the stellar mass tending to increase along tracks of $q_{\rm IR}$ versus $α$ in the direction expected, providing the first observational evidence in support of these models.

preprint2014arXiv

The SCUBA-2 Cosmology Legacy Survey: ALMA resolves the rest-frame far-infrared emission of sub-millimeter galaxies

We present high-resolution (0.3'') ALMA 870um imaging of 52 sub-millimeter galaxies (SMGs) in the Ultra Deep Survey (UDS) field and investigate the size and morphology of the sub-millimeter (sub-mm) emission on 2-10kpc scales. We derive a median intrinsic angular size of FWHM=0.30$\pm$0.04'' for the 23 SMGs in the sample detected at a signal-to-noise ratio (SNR) >10. Using the photometric redshifts of the SMGs we show that this corresponds to a median physical half-light diameter of 2.4$\pm$0.2kpc. A stacking analysis of the SMGs detected at an SNR <10 shows they have sizes consistent with the 870um-bright SMGs in the sample. We compare our results to the sizes of SMGs derived from other multi-wavelength studies, and show that the rest-frame ~250um sizes of SMGs are consistent with studies of resolved 12CO (J=3-2 to 7-6) emission lines, but that sizes derived from 1.4GHz imaging appear to be approximately two times larger on average, which we attribute to cosmic ray diffusion. The rest-frame optical sizes of SMGs are around four times larger than the sub-millimeter sizes, indicating that the star formation in these galaxies is compact relative to the pre-existing stellar distribution. The size of the starburst region in SMGs is consistent with the majority of the star formation occurring in a central region, a few kpc in extent, with a median star formation rate surface density of 90$\pm$30Msol/yr/kpc$^2$, which may suggest that we are witnessing an intense period of bulge growth in these galaxies.

preprint2013arXiv

A Dust-Obscured Massive Maximum-Starburst Galaxy at a Redshift of 6.34

Massive present-day early-type (elliptical and lenticular) galaxies probably gained the bulk of their stellar mass and heavy elements through intense, dust-enshrouded starbursts - that is, increased rates of star formation - in the most massive dark matter halos at early epochs. However, it remains unknown how soon after the Big Bang such massive starburst progenitors exist. The measured redshift distribution of dusty, massive starbursts has long been suspected to be biased low in redshift owing to selection effects, as confirmed by recent findings of systems out to redshift z~5. Here we report the identification of a massive starburst galaxy at redshift 6.34 through a submillimeter color-selection technique. We unambiguously determined the redshift from a suite of molecular and atomic fine structure cooling lines. These measurements reveal a hundred billion solar masses of highly excited, chemically evolved interstellar medium in this galaxy, which constitutes at least 40% of the baryonic mass. A "maximum starburst" converts the gas into stars at a rate more than 2,000 times that of the Milky Way, a rate among the highest observed at any epoch. Despite the overall downturn of cosmic star formation towards the highest redshifts, it seems that environments mature enough to form the most massive, intense starbursts existed at least as early as 880 million years after the Big Bang.

preprint2013arXiv

The SCUBA-2 Cosmology Legacy Survey: Ultraluminous star-forming galaxies in a z=1.6 cluster

We analyse new SCUBA-2 submillimeter and archival SPIRE far-infrared imaging of a z=1.62 cluster, Cl0218.3-0510, which lies in the UKIDSS/UDS field of the SCUBA-2 Cosmology Legacy Survey. Combining these tracers of obscured star formation activity with the extensive photometric and spectroscopic information available for this field, we identify 31 far-infrared/submillimeter-detected probable cluster members with bolometric luminosities >1e12 Lo and show that by virtue of their dust content and activity, these represent some of the reddest and brightest galaxies in this structure. We exploit Cycle-1 ALMA submillimeter continuum imaging which covers one of these sources to confirm the identification of a SCUBA-2-detected ultraluminous star-forming galaxy in this structure. Integrating the total star-formation activity in the central region of the structure, we estimate that it is an order of magnitude higher (in a mass-normalised sense) than clusters at z~0.5-1. However, we also find that the most active cluster members do not reside in the densest regions of the structure, which instead host a population of passive and massive, red galaxies. We suggest that while the passive and active populations have comparable near-infrared luminosities at z=1.6, M(H)~-23, the subsequent stronger fading of the more active galaxies means that they will evolve into passive systems at the present-day which are less luminous than the descendants of those galaxies which were already passive at z~1.6 (M(H)~-20.5 and M(H)~-21.5 respectively at z~0). We conclude that the massive galaxy population in the dense cores of present-day clusters were already in place at z=1.6 and that in Cl0218.3-0510 we are seeing continuing infall of less extreme, but still ultraluminous, star-forming galaxies onto a pre-existing structure.

preprint2012arXiv

VLA imaging of 12CO J=1-0 and free-free emission in lensed submillimetre galaxies

We present a study using the Karl G. Jansky Very Large Array of 12CO J=1-0 emission in three strongly lensed submillimetre-selected galaxies (SMMJ16359, SMMJ14009 and SMMJ02399) at z=2.5-2.9. These galaxies span L(IR) = 10^11 - 10^13 Lsun, offering an opportunity to compare the interstellar medium of LIRGs and ULIRGs at high redshift. We estimate molecular gas masses in the range (2-40) x 10^9 Msun using a method that assumes canonical underlying brightness temperature ratios for star-forming and non-star-forming gas phases and a maximal star-formation efficiency. A more simplistic method using X(CO) = 0.8 yields gas masses twice as high. The observed CO(3-2)/CO(1-0) brightness temperature ratio for SMMJ14009, r(3-2)/(1-0) = (0.95 \pm 0.12), is indicative of warm star-forming gas, possibly influenced by the central AGN. We search for 12CO(1-0) emission in the Lyman-break galaxy, A2218 #384, located at z=2.517 in the same field as SMMJ16359, and assign a 3-sigma gas mass limit of <6 x 10^8 Msun. We use rest-frame 115-GHz free-free flux densities in SMMJ14009 and SMMJ02399 - measurements tied directly to the photionisation rate of massive stars and made possible by the VLA's bandwidth - to estimate star-formation rates of 400-600 Msun/yr and to estimate the fraction of L(IR) due to the AGN.

preprint2010arXiv

Tracing the molecular gas in distant submillimetre galaxies via CO(1-0) imaging with the EVLA

We report the results of a pilot study with the EVLA of 12CO J=1-0 emission from four SMGs at z=2.2-2.5, each with an existing detection of CO J=3-2. Using the EVLA's most compact configuration we detect strong, broad J=1-0 line emission from all of our targets. The median line width ratio, sigma(1-0)/sigma(3-2) = 1.15 +/- 0.06, suggests that the J=1-0 is more spatially extended than the J=3-2 emission, a situation confirmed by our maps which reveal velocity structure in several cases and typical sizes of ~16 kpc FWHM. The median Tb ratio is r(3-2/1-0) = 0.55 +/- 0.05, noting that our value may be biased high because of the J=3-2-based sample selection. Naively, this suggests gas masses ~2x higher than estimates made using higher-J transitions of CO, with the discrepency due to the difference in assumed Tb ratio. We also estimate masses using the 12CO J=1-0 line and the observed global Tb ratios, assuming standard underlying Tb ratios as well as a limiting SFE, i.e. without calling upon X(CO). Using this new method, we find a median molecular gas mass of (2.5 +/- 0.8) x 10^10 Msun, with a plausible range stretching 3x higher. Even larger masses cannot be ruled out, but are not favoured by dynamical constraints: the median dynamical mass for our sample is (2.3 +/- 1.4) x 10^11 Msun. We examine the Schmidt-Kennicutt relation for all the distant galaxy populations for which CO J=1-0 or J=2-1 data are available, finding small systematic differences. These have previously been interpreted as evidence for different modes of star formation, but we argue that these differences are to be expected, given the still considerable uncertainties. Finally, we discuss the morass of degeneracies surrounding molecular gas mass estimates, the possibilities for breaking them, and the future prospects for imaging and studying cold, quiescent molecular gas at high redshifts [abridged].