Source author record

S. Dye

S. Dye 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

86works
8topics
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

86 published item(s)

preprint2022arXiv

A Call to Arms Control: Synergies between Nonproliferation Applications of Neutrino Detectors and Large-Scale Fundamental Neutrino Physics Experiments

The High Energy Physics community can benefit from a natural synergy in research activities into next-generation large-scale water and scintillator neutrino detectors, now being studied for remote reactor monitoring, discovery and exclusion applications in cooperative nonproliferation contexts. Since approximately 2010, US nonproliferation researchers, supported by the National Nuclear Security Administration (NNSA), have been studying a range of possible applications of relatively large (100 ton) to very large (hundreds of kiloton) water and scintillator neutrino detectors. In parallel, the fundamental physics community has been developing detectors at similar scales and with similar design features for a range of high-priority physics topics, primarily in fundamental neutrino physics. These topics include neutrino oscillation studies at beams and reactors, solar, and geological neutrino measurements, supernova studies, and others. Examples of ongoing synergistic work at U.S. national laboratories and universities include prototype gadolinium-doped water and water-based and opaque scintillator test-beds and demonstrators, extensive testing and industry partnerships related to large area fast position-sensitive photomultiplier tubes, and the development of concepts for a possible underground kiloton-scale water-based detector for reactor monitoring and technology demonstrations. Some opportunities for engagement between the two communities include bi-annual Applied Antineutrino Physics conferences, collaboration with U.S. National Laboratories engaging in this research, and occasional NNSA funding opportunities supporting a blend of nonproliferation and basic science R&D, directed at the U.S. academic community.

preprint2022arXiv

The Bright Extragalactic ALMA Redshift Survey (BEARS) I: redshifts of bright gravitationally-lensed galaxies from the Herschel ATLAS

We present spectroscopic measurements for 71 galaxies associated with 62 of the brightest high-redshift submillimeter sources from the Southern fields of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), while targeting 85 sources which resolved into 142. We have obtained robust redshift measurements for all sources using the 12-m Array and an efficient tuning of ALMA to optimise its use as a redshift hunter, with 73 per cent of the sources having a robust redshift identification. Nine of these redshift identifications also rely on observations from the Atacama Compact Array. The spectroscopic redshifts span a range $1.41<z<4.53$ with a mean value of 2.75, and the CO emission line full-width at half-maxima range between $\rm 110\,km\,s^{-1} < FWHM < 1290\,km\,s^{-1}$ with a mean value of $\sim$ 500kms$^{-1}$, in line with other high-$z$ samples. The derived CO(1-0) luminosity is significantly elevated relative to line-width to CO(1-0) luminosity scaling relation, which is suggestive of lensing magnification across our sources. In fact, the distribution of magnification factors inferred from the CO equivalent widths is consistent with expectations from galaxy-galaxy lensing models, though there is a hint of an excess at large magnifications that may be attributable to the additional lensing optical depth from galaxy groups or clusters.

preprint2020arXiv

[NII] fine-structure emission at 122 and 205um in a galaxy at z=2.6: a globally dense star-forming interstellar medium

We present new observations with the Atacama Large Millimeter/sub-millimeter Array of the 122um and 205um fine-structure line emission of singly-ionised nitrogen in a strongly lensed starburst galaxy at z=2.6. The 122/205um [NII] line ratio is sensitive to electron density, n_e, in the ionised interstellar medium, and we use this to measure n_e~300cm^-3 averaged across the galaxy. This is over an order of magnitude higher than the Milky Way average, but comparable to localised Galactic star-forming regions. Combined with observations of the atomic carbon (CI(1-0)) and carbon monoxide (CO(4-3)) in the same system, we reveal the conditions in this intensely star-forming system. The majority of the molecular interstellar medium has been driven to high density, and the resultant conflagration of star formation produces a correspondingly dense ionised phase, presumably co-located with myriad HII regions that litter the gas-rich disk.

preprint2020arXiv

Close-up view of a luminous star-forming galaxy at z=2.95

(Abridged) Exploiting the sensitivity and broad band width of NOEMA, we have studied the molecular gas and dust in the galaxy HerBS-89a, at z=2.95. High angular resolution images reveal a partial 1.0" diameter Einstein ring in the dust continuum emission and the molecular emission lines of 12CO(9-8) and H2O(2_02-1_11). We report the detection of the three fundamental transitions of the molecular ion OH+, seen in absorption; the molecular ion CH+(1-0) seen in absorption (and tentatively in emission); two transitions of amidogen (NH2), seen in emission; and HCN(11-10) and/or NH(1_2-0_1) seen in absorption. The NOEMA data are complemented with VLA data tracing the 12CO(1-0) emission line, which provides a measurement of the total mass of molecular gas and an anchor for a CO excitation analysis. In addition, we present HST imaging that reveals the foreground lensing galaxy in the near-infrared. Together with data from the GTC, we derive a photometric redshift of z(phot)~0.9 for the foreground lensing galaxy. Modelling the lensing of HerBS-89a, we reconstruct the dust continuum and molecular emission lines (magnified by a factor ~4-5) in the source plane. The 12CO(9-8) and H2O emission lines have comparable spatial and kinematic distributions; the source-plane reconstructions do not clearly distinguish between a one-component and a two-component scenario, but the latter accounts for the observed broad line widths. HerBS-89a is a powerful star forming galaxy with a dust-to-gas ratio delta(GDR)~80, a SFR = 614 +/- 59 Msun/yr and a depletion timescale tau(depl) = (3.4 +/- 1.0) 1e8 years. The OH+ and CH+ absorption lines, all have their main velocity component red-shifted by Δ(V)~100 km/s relative to the global CO reservoir. We argue that these absorption lines trace a rare example of gas inflow towards the center of the galaxy.

preprint2016arXiv

GAMA/H-ATLAS: Common star-formation rate indicators and their dependence on galaxy physical parameters

We compare common star-formation rate (SFR) indicators in the local Universe in the GAMA equatorial fields (around 160 sq. deg.), using ultraviolet (UV) photometry from GALEX, far-infrared (FIR) and sub-millimetre (sub-mm) photometry from H-ATLAS, and Halpha spectroscopy from the GAMA survey. With a high-quality sample of 745 galaxies (median redshift 0.08), we consider three SFR tracers: UV luminosity corrected for dust attenuation using the UV spectral slope beta (SFRUV,corr), Halpha line luminosity corrected for dust using the Balmer decrement (BD) (SFRHalpha,corr), and the combination of UV and IR emission (SFRUV+IR). We demonstrate that SFRUV,corr can be reconciled with the other two tracers after applying attenuation corrections by calibrating IRX (i.e. the IR to UV luminosity ratio) and attenuation in the Halpha (derived from BD) against beta. However, beta on its own is very unlikely to be a reliable attenuation indicator. We find that attenuation correction factors depend on parameters such as stellar mass, z and dust temperature (Tdust), but not on Halpha equivalent width (EW) or Sersic index. Due to the large scatter in the IRX vs beta correlation, when compared to SFRUV+IR, the beta-corrected SFRUV,corr exhibits systematic deviations as a function of IRX, BD and Tdust.

preprint2016arXiv

H-ATLAS: The Far-Infrared properties of galaxies in and around the Coma Cluster

We describe a far infrared survey of the Coma cluster and the galaxy filament it resides within. Our survey covers an area of $\sim$150 deg$^2$ observed by $Herschel$ H-ATLAS in five bands at 100, 160, 250, 350 and 500$μ$m. The SDSS spectroscopic survey ($m_{r} \le 17.8)$ is used to define an area (within the Virial radius) and redshift selected ($4268 < v < 9700$ km s$^{-1}$) sample of 744 Coma cluster galaxies - the Coma Cluster Catalogue (CCC). For comparison we also define a sample of 951 galaxies in the connecting filament - the Coma Filament Catalogue (CFC). The optical positions and parameters are used to define appropriate apertures to measure each galaxy's far infrared emission. We have detected 99 of 744 (13\%) and 422 of 951 (44\%) of the cluster and filament galaxies in the SPIRE 250$μ$m band. We consider the relative detection rates of galaxies of different morphological types finding that it is only the S0/Sa population that shows clear differences between the cluster and filament. We find no differences between the dust masses and temperatures of cluster and filament galaxies with the exception of early type galaxy dust temperatures, which are significantly hotter in the cluster than in the filament (X-ray heating?). From a chemical evolution model we find no evidence for different evolutionary processes (gas loss or infall) between galaxies in the cluster and filament.

preprint2016arXiv

Submillimeter H2O and H2O+ emission in lensed ultra- and hyper-luminous infrared galaxies at z ~ 2-4

(abridged) We report rest-frame submillimeter H2O emission line observations of 11 HyLIRGs/ULIRGs at z~2-4 selected among the brightest lensed galaxies discovered in the Herschel-ATLAS. Using the IRAM NOEMA, we have detected 14 new H2O emission lines. The apparent luminosities of the H2O emission lines are $μL_{\rm{H_2O}} \sim 6-21 \times 10^8 L_\odot$, with velocity-integrated line fluxes ranging from 4-15 Jy km s$^{-1}$. We have also observed CO emission lines using EMIR on the IRAM 30m telescope in seven sources. The velocity widths for CO and H2O lines are found to be similar. With almost comparable integrated flux densities to those of the high-J CO line, H2O is found to be among the strongest molecular emitters in high-z Hy/ULIRGs. We also confirm our previously found correlation between luminosity of H2O ($L_{\rm{H_2O}}$) and infrared ($L_{\rm{IR}}$) that $L_{\rm{H_2O}} \sim L_{\rm{IR}}^{1.1-1.2}$, with our new detections. This correlation could be explained by a dominant role of far-infrared (FIR) pumping in the H2O excitation. Modelling reveals the FIR radiation fields have warm dust temperature $T_\rm{warm}$~45-75 K, H2O column density per unit velocity interval $N_{\rm{H_2O}}/ΔV \gtrsim 0.3 \times 10^{15}$ cm$^{-2}$ km$^{-1}$ s and 100 $μ$m continuum opacity $τ_{100} > 1$ (optically thick), indicating that H2O is likely to trace highly obscured warm dense gas. However, further observations of $J\geq4$ H2O lines are needed to better constrain the continuum optical depth and other physical conditions of the molecular gas and dust. We have also detected H2O+ emission in three sources. A tight correlation between $L_{\rm{H_2O}}$ and $L_{\rm{H_2O^+}}$ has been found in galaxies from low to high redshift. The velocity-integrated flux density ratio between H2O+ and H2O suggests that cosmic rays generated by strong star formation are possibly driving the H2O+ formation.

preprint2016arXiv

The faint end of the 250 micron luminosity function at z < 0.5

Aims. We aim to study the 250 micron luminosity function (LF) down to much fainter luminosities than achieved by previous efforts. Methods. We developed a modified stacking method to reconstruct the 250 micron LF using optically selected galaxies from the SDSS survey and Herschel maps of the GAMA equatorial fields and Stripe 82. Our stacking method not only recovers the mean 250 micron luminosities of galaxies that are too faint to be individually detected, but also their underlying distribution functions. Results. We find very good agreement with previous measurements in the overlapping luminosity range. More importantly, we are able to derive the LF down to much fainter luminosities (around 25 times fainter) than achieved by previous studies. We find strong positive luminosity evolution \propto (1 + z)^4.89\pm1.07 and moderate negative density evolution \propto (1 + z)^-1.02\pm0.54 over the redshift range z=[0.02, 0.5].

preprint2016arXiv

The Herschel-ATLAS Data Release 1 Paper I: Maps, Catalogues and Number Counts

We present the first major data release of the largest single key-project in area carried out in open time with the Herschel Space Observatory. The Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) is a survey of 600 deg^2 in five photometric bands - 100, 160, 250, 350 and 500 um - with the PACS and SPIRE cameras. In this paper and a companion paper (Bourne et al. 2016) we present the survey of three fields on the celestial equator, covering a total area of 161.6 deg^2 and previously observed in the Galaxy and Mass Assembly (GAMA) spectroscopic survey. This paper describes the Herschel images and catalogues of the sources detected on the SPIRE 250 um images. The 1-sigma noise for source detection, including both confusion and instrumental noise, is 7.4, 9.4 and 10.2 mJy at 250, 350 and 500 um. Our catalogue includes 120230 sources in total, with 113995, 46209 and 11011 sources detected at >4-sigma at 250, 350 and 500 um. The catalogue contains detections at >3-sigma at 100 and 160 um for 4650 and 5685 sources, and the typical noise at these wavelengths is 44 and 49 mJy. We include estimates of the completeness of the survey and of the effects of flux bias and also describe a novel method for determining the true source counts. The H-ATLAS source counts are very similar to the source counts from the deeper HerMES survey at 250 and 350 um, with a small difference at 500 um. Appendix A provides a quick start in using the released datasets, including instructions and cautions on how to use them.

preprint2016arXiv

The Herschel-ATLAS Data Release 1 Paper II: Multi-wavelength counterparts to submillimetre sources

This paper is the second in a pair of articles presenting data release 1 (DR1) of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), the largest single open-time key project carried out with the Herschel Space Observatory. The H-ATLAS is a wide-area imaging survey carried out in five photometric bands at 100, 160, 250, 350 and 500$μ$m covering a total area of 600deg$^2$. In this paper we describe the identification of optical counterparts to submillimetre sources in DR1, comprising an area of 161 deg$^2$ over three equatorial fields of roughly 12$^\circ$x4.5$^\circ$ centred at 9$^h$, 12$^h$ and 14.5$^h$ respectively. Of all the H-ATLAS fields, the equatorial regions benefit from the greatest overlap with current multi-wavelength surveys spanning ultraviolet (UV) to mid-infrared regimes, as well as extensive spectroscopic coverage. We use a likelihood-ratio technique to identify SDSS counterparts at r<22.4 for 250-$μ$m-selected sources detected at $\geq$ 4$σ$ ($\approx$28mJy). We find `reliable' counterparts (reliability R$\geq$0.8) for 44,835 sources (39 per cent), with an estimated completeness of 73.0 per cent and contamination rate of 4.7 per cent. Using redshifts and multi-wavelength photometry from GAMA and other public catalogues, we show that H-ATLAS-selected galaxies at $z<0.5$ span a wide range of optical colours, total infrared (IR) luminosities, and IR/UV ratios, with no strong disposition towards mid-IR-classified AGN in comparison with optical selection. The data described herein, together with all maps and catalogues described in the companion paper (Valiante et al. 2016), are available from the H-ATLAS website at www.h-atlas.org.

preprint2016arXiv

The Herschel-ATLAS: a sample of 500μm-selected lensed galaxies over 600 square degrees

We present a sample of 80 candidate strongly lensed galaxies with flux density above 100mJy at 500μm extracted from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), over an area of 600 square degrees. Available imaging and spectroscopic data allow us to confirm the strong lensing in 20 cases and to reject it in one case. For other 8 objects the lensing scenario is strongly supported by the presence of two sources along the same line of sight with distinct photometric redshifts. The remaining objects await more follow-up observations to confirm their nature. The lenses and the background sources have median redshifts z_L = 0.6 and z_S = 2.5, respectively, and are observed out to z_L = 1.2 and z_S = 4.2. We measure the number counts of candidate lensed galaxies at 500μm and compare them with theoretical predictions, finding a good agreement for a maximum magnification of the background sources in the range 10-20. These values are consistent with the magnification factors derived from the lens modelling of individual systems. The catalogue presented here provides sub- mm bright targets for follow-up observations aimed at exploiting gravitational lensing to study with un-precedented details the morphological and dynamical properties of dusty star forming regions in z >~ 1.5 galaxies.

preprint2015arXiv

A multi-wavelength exploration of the [CII]/IR ratio in H-ATLAS/GAMA galaxies out to z=0.2

We explore the behaviour of [CII]-157.74um forbidden fine-structure line observed in a sample of 28 galaxies selected from ~50deg^2 of the H-ATLAS survey. The sample is restricted to galaxies with flux densities higher than S_160um>150mJy and optical spectra from the GAMA survey at 0.02<z<0.2. Far-IR spectra centred on this redshifted line were taken with the PACS instrument on-board the Herschel Space Observatory. The galaxies span 10<log(L_IR/Lo)<12 (where L_IR=L_IR[8-1000um]) and 7.3<log(L_[CII]/Lo)<9.3, covering a variety of optical galaxy morphologies. The sample exhibits the so-called [CII] deficit at high IR luminosities, i.e. L_[CII]/L_IR (hereafter [CII]/IR) decreases at high L_IR. We find significant differences between those galaxies presenting [CII]/IR>2.5x10^-3 with respect to those showing lower ratios. In particular, those with high ratios tend to have: (1) L_IR<10^11Lo; (2) cold dust temperatures, T_d<30K; (3) disk-like morphologies in r-band images; (4) a WISE colour 0.5<S_12um/S_22um<1.0; (5) low surface brightness Sigma_IR~10^8-9 Lo kpc^-2, (6) and specific star-formation rates of sSFR~0.05-3 Gyr^-1. We suggest that the strength of the far-UV radiation fields (<G_O>) is main parameter responsible for controlling the [CII]/IR ratio. It is possible that relatively high <G_O> creates a positively charged dust grain distribution, impeding an efficient photo-electric extraction of electrons from these grains to then collisionally excite carbon atoms. Within the brighter IR population, 11<log(L_IR/Lo)<12, the low [CII]/IR ratio is unlikely to be modified by [CII] self absorption or controlled by the presence of a moderately luminous AGN (identified via the BPT diagram).

preprint2015arXiv

Cross-correlation between the CMB lensing potential measured by Planck and high-z sub-mm galaxies detected by the Herschel-ATLAS survey

We present the first measurement of the correlation between the map of the cosmic microwave background (CMB) lensing potential derived from the \emph{Planck} nominal mission data and $z\gtrsim 1.5$ galaxies detected by the \emph{Herschel}-ATLAS (H-ATLAS) survey covering about $600\,\hbox{deg}^2$, i.e. about 1.4\% of the sky. We reject the hypothesis that there is no correlation between CMB lensing and galaxy detection at a $20\,σ$ significance, checking the result by performing a number of null tests. The significance of the detection of the theoretically expected cross-correlation signal is found to be $10\,σ$. The galaxy bias parameter, $b$, derived from a joint analysis of the cross-power spectrum and of the auto-power spectrum of the galaxy density contrast is found to be $b=2.80^{+0.12}_{-0.11}$, consistent with earlier estimates for H-ATLAS galaxies at similar redshifts. On the other hand, the amplitude of the cross-correlation is found to be a factor $1.62 \pm 0.16$ higher than expected from the standard model and also found by cross-correlation analyses with other tracers of the large-scale structure. The enhancement due to lensing magnification can account for only a fraction of the excess cross-correlation signal. We suggest that part of it may be due to an incomplete removal of the contamination of the CIB, that includes the H-ATLAS sources we are cross-correlating with. In any case, the highly significant detection reported here using a catalog covering only 1.4\% of the sky demonstrates the potential of CMB lensing correlations with submillimeter surveys.

preprint2015arXiv

H-ATLAS/GAMA and HeViCS - Dusty Early-Type Galaxies in Different Environments

The Herschel Space Observatory has had a tremendous impact on the study of extragalactic dust. Specifically, early-type galaxies (ETG) have been the focus of several studies. In this paper we combine results from two Herschel studies - a Virgo cluster study HeViCS and a broader, low-redshift H-ATLAS/GAMA study - and contrast the dust and associated properties for similar mass galaxies. This comparison is motivated by differences in results exhibited between multiple Herschel studies of early-type galaxies. A comparison between consistent modified blackbody derived dust mass is carried out, revealing strong differences between the two samples in both dust mass and dust-to-stellar mass ratio. In particular, the HeViCS sample lacks massive ETG with as high a specific dust content as found in H-ATLAS. This is most likely connected with the difference in environment for the two samples. We calculate nearest neighbour environment densities in a consistent way, showing that H-ATLAS ETG occupy sparser regions of the local Universe, whereas HeViCS ETG occupy dense regions. This is also true for ETG that are not Herschel-detected but are in the Virgo and GAMA parent samples. Spectral energy distributions are fit to the panchromatic data. From these we find that in H-ATLAS the specific star formation rate anticorrelates with stellar mass and reaches values as high as in our Galaxy. On the other hand HeViCS ETG appear to have little star formation. Based on the trends found here, H-ATLAS ETG are thought to have more extended star formation histories and a younger stellar population than HeViCS ETG.

preprint2015arXiv

Revealing the complex nature of the strong gravitationally lensed system H-ATLAS J090311.6+003906 using ALMA

We have modelled Atacama Large Millimeter/sub-millimeter Array (ALMA) long baseline imaging of the strong gravitational lens system H-ATLAS J090311.6+003906 (SDP.81). We have reconstructed the distribution of band 6 and 7 continuum emission in the z=3.042 source and we have determined its kinematic properties by reconstructing CO(5-4) and CO(8-7) line emission in bands 4 and 6. The continuum imaging reveals a highly non-uniform distribution of dust with clumps on scales of 200 pc. In contrast, the CO line emission shows a relatively smooth, disk-like velocity field which is well fit by a rotating disk model with an inclination angle of 40+/-5 degrees and an asymptotic rotation velocity of 320 km/s. The inferred dynamical mass within 1.5kpc is (3.5+/-0.5)x10^{10} M_sol which is comparable to the total molecular gas masses of (2.7+/-0.5)x10^{10} M_sol and (3.5+/-0.6)x10^{10} M_sol from the dust continuum emission and CO emission respectively. Our new reconstruction of the lensed HST near-infrared emission shows two objects which appear to be interacting, with the rotating disk of gas and dust revealed by ALMA distinctly offset from the near-infrared emission. The clumpy nature of the dust and a low value of the Toomre parameter of Q=0.3 suggest that the disk is in a state of collapse. We estimate a star formation rate in the disk of 470+/-80 M_sol/yr with an efficiency 65 times greater than typical low-redshift galaxies. Our findings add to the growing body of evidence that the most infra-red luminous, dust obscured galaxies in the high redshift Universe represent a population of merger induced starbursts.

preprint2014arXiv

Colour matters: the effects of lensing on the positional offsets between optical and submillimetre galaxies in Herschel-ATLAS

We report an unexpected variation in the positional offset distributions between Herschel-ATLAS sub-millimetre (submm) sources and their optical associations, depending on both 250-μm signal-to-noise ratio and 250/350-μm colour. We show that redder and brighter submm sources have optical associations with a broader distribution of positional offsets than would be expected if these offsets were due to random positional errors in the source extraction. The observation can be explained by two possible effects: either red submm sources trace a more clustered population than blue ones, and their positional errors are increased by confusion; or red submm sources are generally at high redshifts and are frequently associated with low-redshift lensing structures which are identified as false counterparts. We perform various analyses of the data, including the multiplicity of optical associations, the redshift and magnitude distributions in H-ATLAS in comparison to HerMES, and simulations of weak lensing, and we conclude that the effects are most likely to be explained by widespread weak lensing of Herschel-SPIRE sources by foreground structures. This has important consequences for counterpart identification and derived redshift distributions and luminosity functions of submm surveys.

preprint2014arXiv

Herschel-ATLAS: Modelling the first strong gravitational lenses

We have determined the mass-density radial profiles of the first five strong gravitational lens systems discovered by the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). We present an enhancement of the semi-linear lens inversion method of Warren & Dye which allows simultaneous reconstruction of several different wavebands and apply this to dual-band imaging of the lenses acquired with the Hubble Space Telescope. The five systems analysed here have lens redshifts which span a range, 0.22<z<0.94. Our findings are consistent with other studies by concluding that: 1) the logarithmic slope of the total mass density profile steepens with decreasing redshift; 2) the slope is positively correlated with the average total projected mass density of the lens contained within half the effective radius and negatively correlated with the effective radius; 3) the fraction of dark matter contained within half the effective radius increases with increasing effective radius and increases with redshift.

preprint2014arXiv

Herschel-ATLAS: Properties of dusty massive galaxies at low and high redshifts

We present a comparison of the physical properties of a rest-frame $250μ$m selected sample of massive, dusty galaxies from $0<z<5.3$. Our sample comprises 29 high-redshift submillimetre galaxies (SMGs) from the literature, and 843 dusty galaxies at $z<0.5$ from the Herschel-ATLAS, selected to have a similar stellar mass to the SMGs. The $z>1$ SMGs have an average SFR of $390^{+80}_{-70}\,$M$_\odot$yr$^{-1}$ which is 120 times that of the low-redshift sample matched in stellar mass to the SMGs (SFR$=3.3\pm{0.2}$ M$_\odot$yr$^{-1}$). The SMGs harbour a substantial mass of dust ($1.2^{+0.3}_{-0.2}\times{10}^9\,$M$_\odot$), compared to $(1.6\pm0.1)\times{10}^8\,$M$_\odot$ for low-redshift dusty galaxies. At low redshifts the dust luminosity is dominated by the diffuse ISM, whereas a large fraction of the dust luminosity in SMGs originates from star-forming regions. At the same dust mass SMGs are offset towards a higher SFR compared to the low-redshift H-ATLAS galaxies. This is not only due to the higher gas fraction in SMGs but also because they are undergoing a more efficient mode of star formation, which is consistent with their bursty star-formation histories. The offset in SFR between SMGs and low-redshift galaxies is similar to that found in CO studies, suggesting that dust mass is as good a tracer of molecular gas as CO.

preprint2014arXiv

Herschel-ATLAS/GAMA: SDSS cross-correlation induced by weak lensing

We report a highly significant ($>10σ$) spatial correlation between galaxies with $S_{350μ\rm m}\ge 30\,$mJy detected in the equatorial fields of the \textsl{Herschel} Astrophysical Terahertz Large Area Survey (H-ATLAS) with estimated redshifts $\gtrsim 1.5$, and SDSS or GAMA galaxies at $0.2\le z\le 0.6$. The significance of the cross-correlation is much higher than those reported so far for samples with non-overlapping redshift distributions selected in other wavebands. Extensive, realistic simulations of clustered sub-mm galaxies amplified by foreground structures confirm that the cross-correlation is explained by weak gravitational lensing ($μ<2$). The simulations also show that the measured amplitude and range of angular scales of the signal are larger than can be accounted for by galaxy-galaxy weak lensing. However, for scales $\lesssim 2\,$arcmin, the signal can be reproduced if SDSS/GAMA galaxies act as signposts of galaxy groups/clusters with halo masses in the range $10^{13.2}$--$10^{14.5} M_{\odot}$. The signal detected on larger scales appears to reflect the clustering of such halos.

preprint2014arXiv

Lens Models of Herschel-Selected Galaxies From High-Resolution Near-IR Observations

We present Keck-Adaptive Optics and Hubble Space Telescope high resolution near-infrared (IR) imaging for 500 um-bright candidate lensing systems identified by the Herschel Multi-tiered Extra-galactic Survey (HerMES) and Herschel Astrophysical Terahertz Survey (H-ATLAS). Out of 87 candidates with near-IR imaging, 15 (~17%) display clear near-IR lensing morphologies. We present near-IR lens models to reconstruct and recover basic rest-frame optical morphological properties of the background galaxies from 12 new systems. Sources with the largest near-IR magnification factors also tend to be the most compact, consistent with the size bias predicted from simulations and pre- vious lensing models for sub-millimeter galaxies. For four new sources that also have high-resolution sub-mm maps, we test for differential lensing between the stellar and dust components and find that the 880 um magnification factor (u_880) is ~1.5 times higher than the near-IR magnification factor (u_NIR), on average. We also find that the stellar emission is ~2 times more extended in size than dust. The rest-frame optical properties of our sample of Herschel-selected lensed SMGs are consistent with those of unlensed SMGs, which suggests that the two populations are similar.

preprint2014arXiv

Planning the Future of U.S. Particle Physics (Snowmass 2013): Chapter 8: Instrumentation Frontier

These reports present the results of the 2013 Community Summer Study of the APS Division of Particles and Fields ("Snowmass 2013") on the future program of particle physics in the U.S. Chapter 8, on the Instrumentation Frontier, discusses the instrumentation needs of future experiments in the Energy, Intensity, and Cosmic Frontiers, promising new technologies for particle physics research, and issues of gathering resources for long-term research in this area.

preprint2014arXiv

The dust budget crisis in high-redshift submillimetre galaxies

We apply a chemical evolution model to investigate the sources and evolution of dust in a sample of 26 high-redshift ($z>1$) submillimetre galaxies (SMGs) from the literature, with complete photometry from ultraviolet to the submillimetre. We show that dust produced only by low-intermediate mass stars falls a factor 240 short of the observed dust masses of SMGs, the well-known `dust-budget crisis'. Adding an extra source of dust from supernovae can account for the dust mass in 19 per cent of the SMG sample. Even after accounting for dust produced by supernovae the remaining deficit in the dust mass budget provides support for higher supernova yields, substantial grain growth in the interstellar medium or a top-heavy IMF. Including efficient destruction of dust by supernova shocks increases the tension between our model and observed SMG dust masses. The models which best reproduce the physical properties of SMGs have a rapid build-up of dust from both stellar and interstellar sources and minimal dust destruction. Alternatively, invoking a top-heavy IMF or significant changes in the dust grain properties can solve the dust budget crisis only if dust is produced by both low mass stars and supernovae and is not efficiently destroyed by supernova shocks.

preprint2013arXiv

COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses XIII: Time delays and 9-yr optical monitoring of the lensed quasar RX J1131-1231

We present the results from nine years of optically monitoring the gravitationally lensed z=0.658 quasar RX J1131-1231. The R-band light curves of the four individual images of the quasar were obtained using deconvolution photometry for a total of 707 epochs. Several sharp quasar variability features strongly constrain the time delays between the quasar images. Using three different numerical techniques, we measure these delays for all possible pairs of quasar images while always processing the four light curves simultaneously. For all three methods, the delays between the three close images A, B, and C are compatible with being 0, while we measure the delay of image D to be 91 days, with a fractional uncertainty of 1.5% (1 sigma), including systematic errors. Our analysis of random and systematic errors accounts in a realistic way for the observed quasar variability, fluctuating microlensing magnification over a broad range of temporal scales, noise properties, and seasonal gaps. Finally, we find that our time-delay measurement methods yield compatible results when applied to subsets of the data.

preprint2013arXiv

Far-infrared spectroscopy of a lensed starburst: a blind redshift from Herschel

We report the redshift of HATLAS J132427.0+284452 (hereafter HATLAS J132427), a gravitationally lensed starburst galaxy, the first determined 'blind' by the Herschel Space Observatory. This is achieved via the detection of [C II] consistent with z = 1.68 in a far-infrared spectrum taken with the SPIRE Fourier Transform Spectrometer. We demonstrate that the [C II] redshift is secure via detections of CO J = 2 - 1 and 3 - 2 using the Combined Array for Research in Millimeter-wave Astronomy and the Institut de Radioastronomie Millimetrique's Plateau de Bure Interferometer. The intrinsic properties appear typical of high-redshift starbursts despite the high lensing-amplified fluxes, proving the ability of the FTS to probe this population with the aid of lensing. The blind detection of [C II] demonstrates the potential of the SAFARI imaging spectrometer, proposed for the much more sensitive SPICA mission, to determine redshifts of multiple dusty galaxies simultaneously without the benefit of lensing.

preprint2013arXiv

Gravitational Lens Models Based on Submillimeter Array Imaging of Herschel-selected Strongly Lensed Sub-millimeter Galaxies at z>1.5

Strong gravitational lenses are now being routinely discovered in wide-field surveys at (sub)millimeter wavelengths. We present Submillimeter Array (SMA) high-spatial resolution imaging and Gemini-South and Multiple Mirror Telescope optical spectroscopy of strong lens candidates discovered in the two widest extragalactic surveys conducted by the Herschel Space Observatory: the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) and the Herschel Multi-tiered Extragalactic Survey (HerMES). From a sample of 30 Herschel sources with S_500>100 mJy, 21 are strongly lensed (multiply imaged), 4 are moderately lensed (singly imaged), and the remainder require additional data to determine their lensing status. We apply a visibility-plane lens modeling technique to the SMA data to recover information about the masses of the lenses as well as the intrinsic (i.e., unlensed) sizes (r_half) and far-infrared luminosities (L_FIR) of the lensed submillimeter galaxies (SMGs). The sample of lenses comprises primarily isolated massive galaxies, but includes some groups and clusters as well. Several of the lenses are located at z_lens>0.7, a redshift regime that is inaccessible to lens searches based on Sloan Digital Sky Survey spectroscopy. The lensed SMGs are amplified by factors that are significantly below statistical model predictions given the 500um flux densities of our sample. We speculate that this may reflect a deficiency in our understanding of the intrinsic sizes and luminosities of the brightest SMGs. The lensed SMGs span nearly one decade in L_FIR (median L_FIR=7.9x10^12 L_sun) and two decades in FIR luminosity surface density (median Sigma_FIR=6.0x10^11 L_sun kpc^-2). The strong lenses in this sample and others identified via (sub-)mm surveys will provide a wealth of information regarding the astrophysics of galaxy formation and evolution over a wide range in redshift.

preprint2013arXiv

H-ATLAS: estimating redshifts of Herschel sources from sub-mm fluxes

Upon its completion the Herschel ATLAS (H-ATLAS) will be the largest submillimetre survey to date, detecting close to half-a-million sources. It will only be possible to measure spectroscopic redshifts for a small fraction of these sources. However, if the rest-frame spectral energy distribution (SED) of a typical H-ATLAS source is known, this SED and the observed Herschel fluxes can be used to estimate the redshifts of the H-ATLAS sources without spectroscopic redshifts. In this paper, we use a subset of 40 H-ATLAS sources with previously measured redshifts in the range 0.5<z<4.2 to derive a suitable average template for high redshift H-ATLAS sources. We find that a template with two dust components T_c = 23.9 K, T_h = 46.9 K and ratio of mass of cold dust to mass of warm dust of 30.1) provides a good fit to the rest-frame fluxes of the sources in our calibration sample. We use a jackknife technique to estimate the accuracy of the redshifts estimated with this template, finding a root mean square of Delta z/(1+z) = 0.26. For sources for which there is prior information that they lie at z > 1 we estimate that the rms of Delta z/(1+z) = 0.12. We have used this template to estimate the redshift distribution for the sources detected in the H-ATLAS equatorial fields, finding a bimodal distribution with a mean redshift of 1.2, 1.9 and 2.5 for 250, 350 and 500 um selected sources respectively. \end{abstract}

preprint2013arXiv

H2O emission in high-z ultra-luminous infrared galaxies

Using IRAM PdBI we report the detection of H2O in six new lensed ultra-luminous starburst galaxies at high redshift, discovered in the Herschel H-ATLAS survey. The sources are detected either in the 2_{02}-1_{11} or 2_{11}-2_{02} H_2O emission lines with integrated line fluxes ranging from 1.8 to 14 Jy.km/s. The corresponding apparent luminosities are mu x L_H2O ~ 3-12 x 10^8 Lo, where mu is the lensing magnification factor (3 < mu < 12). These results confirm that H2O lines are among the strongest molecular lines in such galaxies, with intensities almost comparable to those of the high-J CO lines, and same profiles and line widths (200-900 km/s) as the latter. With the current sensitivity of PdBI, H2O can therefore easily be detected in high-z lensed galaxies (with F(500um) > 100 mJy) discovered in the Herschel surveys. Correcting the luminosities for lensing amplification, L_H2O is found to have a strong dependence on the IR luminosity, varying as ~L_IR^{1.2}. This relation which needs to be confirmed with better statistics, may indicate a role of radiative (IR) excitation of the H2O lines, and implies that high-z galaxies with L_IR >~ 10^13 Lo tend to be very strong emitters in H2O, that have no equivalent in the local universe.

preprint2013arXiv

Herschel-ATLAS: Correlations between Dust and Gas in Local Submm-Selected Galaxies

We present an analysis of CO molecular gas tracers in a sample of 500μm-selected Herschel-ATLAS galaxies at z<0.05 (cz<14990km/s). Using 22-500μm photometry from WISE, IRAS and Herschel, with HI data from the literature, we investigate correlations between warm and cold dust, and tracers of the gas in different phases. The correlation between global CO(3-2) line fluxes and FIR-submillimetre fluxes weakens with increasing IR wavelength (λ>60μm), as a result of colder dust being less strongly associated with dense gas. Conversely, CO(2-1) and HI line fluxes both appear to be better correlated with longer wavelengths, suggesting that cold dust is more strongly associated with diffuse atomic and molecular gas phases, consistent with it being at least partially heated by radiation from old stellar populations. The increased scatter at long wavelengths implies that submillimetre fluxes are a poorer tracer of SFR. Fluxes at 22 and 60μm are also better correlated with diffuse gas tracers than dense CO(3-2), probably due to very-small-grain emission in the diffuse interstellar medium, which is not correlated with SFR. The FIR/CO luminosity ratio and the dust mass/CO luminosity ratio both decrease with increasing luminosity, as a result of either correlations between mass and metallicity (changing CO/H2) or between CO luminosity and excitation [changing CO(3-2)/CO(1-0)].

preprint2013arXiv

Herschel-ATLAS: deep HST/WFC3 imaging of strongly lensed submillimeter galaxies

We report on deep near-infrared observations obtained with the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope (HST) of the first five confirmed gravitational lensing events discovered by the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). We succeed in disentangling the background galaxy from the lens to gain separate photometry of the two components. The HST data allow us to significantly improve on previous constraints of the mass in stars of the lensed galaxy and to perform accurate lens modelling of these systems, as described in the accompanying paper by Dye et al. We fit the spectral energy distributions of the background sources from near-IR to millimetre wavelengths and use the magnification factors estimated by Dye et al. to derive the intrinsic properties of the lensed galaxies. We find these galaxies to have star-formation rates of approximately 400 to 2000 M_sol/yr, with approximately (6-25)x10^10 M_sol of their baryonic mass already turned into stars. At these rates of star formation, all remaining molecular gas will be exhausted in less than 100 Myr, reaching a final mass in stars of a few 10^11 M_sol. These galaxies are thus proto-ellipticals caught during their major episode of star formation, and observed at the peak epoch z=1.5-3 of the cosmic star formation history of the Universe.

preprint2013arXiv

Neutrino mass hierarchy determination and other physics potential of medium-baseline reactor neutrino oscillation experiments

Medium-baseline reactor neutrino oscillation experiments (MBRO) have been proposed to determine the neutrino mass hierarchy (MH) and to make precise measurements of the neutrino oscillation parameters. With sufficient statistics, better than ~3%/\sqrt{E} energy resolution and well understood energy non-linearity, MH can be determined by analyzing oscillation signals driven by the atmospheric mass-squared difference in the survival spectrum of reactor antineutrinos. With such high performance MBRO detectors, oscillation parameters, such as \sin^22θ_{12}, Δm^2_{21}, and Δm^2_{32}, can be measured to sub-percent level, which enables a future test of the PMNS matrix unitarity to ~1% level and helps the forthcoming neutrinoless double beta decay experiments to constrain the allowed <m_{ββ}> values. Combined with results from the next generation long-baseline beam neutrino and atmospheric neutrino oscillation experiments, the MH determination sensitivity can reach higher levels. In addition to the neutrino oscillation physics, MBRO detectors can also be utilized to study geoneutrinos, astrophysical neutrinos and proton decay. We propose to start a U.S. R&D program to identify, quantify and fulfill the key challenges essential for the success of MBRO experiments.

preprint2013arXiv

The suppression of star formation by powerful active galactic nuclei

The old, red stars which constitute the bulges of galaxies, and the massive black holes at their centres, are the relics of a period in cosmic history when galaxies formed stars at remarkable rates and active galactic nuclei (AGN) shone brightly from accretion onto black holes. It is widely suspected, but unproven, that the tight correlation in mass of the black hole and stellar components results from the AGN quenching the surrounding star formation as it approaches its peak luminosity. X-rays trace emission from AGN unambiguously, while powerful star-forming galaxies are usually dust-obscured and are brightest at infrared to submillimetre wavelengths. Here we report observations in the submillimetre and X-ray which show that rapid star formation was common in the host galaxies of AGN when the Universe was 2-6 Gyrs old, but that the most vigorous star formation is not observed around black holes above an X-ray luminosity of 10^44 erg/s. This suppression of star formation in the host galaxies of powerful AGN is a key prediction of models in which the AGN drives a powerful outflow, expelling the interstellar medium of its host galaxy and transforming the galaxy's properties in a brief period of cosmic time.

preprint2012arXiv

A Comprehensive View of a Strongly Lensed Planck-Associated Submillimeter Galaxy

We present high-resolution maps of stars, dust, and molecular gas in a strongly lensed submillimeter galaxy (SMG) at z = 3.259. HATLAS12--00 is selected from the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) as a strong lens candidate mainly based on its unusually high 500um flux density (~300 mJy). It is the only high-redshift Planck detection in the 130 deg^2 H-ATLAS Phase 1 area. Keck Adaptive Optics images reveal a quadruply imaged galaxy in the K-band while the Submillimeter Array and the Extended Very Large Array show doubly imaged 880um and CO(1-0) sources, indicating differentiated distributions of the various components in the galaxy. In the source plane, the stars reside in three major kpc-scale clumps extended over ~1.6 kpc, the dust in a compact (~1 kpc) region ~3 kpc north of the stars, and the cold molecular gas in an extended (~7 kpc) disk ~5 kpc northeast of the stars. The emission from the stars, dust, and gas are magnified by ~17, 8, and 7 times, respectively, by four lensing galaxies at z ~ 1. Intrinsically, the lensed galaxy is a warm (T_dust ~ 40-65 K), hyper-luminous (L_IR ~ 1.7e13 Lsun; SFR ~ 2000 Msun/yr), gas-rich (M_gas/M_baryon ~ 70%), young (M_stellar/SFR ~ 20 Myr), and short-lived (M_gas/SFR ~ 40 Myr) starburst, without a significant active galactic nucleus. With physical properties similar to unlensed z > 2 SMGs, HATLAS12--00 offers a detailed view of a typical SMG through a powerful cosmic microscope.

preprint2012arXiv

A Detailed Gravitational Lens Model Based on Submillimeter Array and Keck Adaptive Optics Imaging of a Herschel-ATLAS Sub-millimeter Galaxy at z=4.243

We present high-spatial resolution imaging obtained with the Submillimeter Array (SMA) at 880um and the Keck Adaptive Optics (AO) system at Ks-band of a gravitationally lensed sub-millimeter galaxy (SMG) at z=4.243 discovered in the Herschel-Astrophysical Terahertz Large Area Survey. The SMA data (angular resolution ~0.6") resolve the dust emission into multiple lensed images, while the Keck AO Ks-band data (angular resolution ~0.1") resolve the lens into a pair of galaxies separated by 0.3". We present an optical spectrum of the foreground lens obtained with the Gemini-South telescope that provides a lens redshift of z_lens = 0.595 +/- 0.005. We develop and apply a new lens modeling technique in the visibility plane that shows that the SMG is magnified by a factor of mu = 4.1 +/- 0.2 and has an intrinsic infrared (IR) luminosity of L_IR = (2.1 +/- 0.2) x 10^13 Lsun. We measure a half-light radius of the background source of r_s = 4.4 +/- 0.5 kpc which implies an IR luminosity surface density of Sigma_IR = (3.4 +/- 0.9) x 10^11 Lsun kpc^-2, a value that is typical of z > 2 SMGs but significantly lower than IR luminous galaxies at z~0. The two lens galaxies are compact (r_lens ~ 0.9 kpc) early-types with Einstein radii of theta_E1 = 0.57 +/- 0.01 and theta_E2 = 0.40 +/- 0.01 that imply masses of M_lens1 = (7.4 +/- 0.5) x 10^10 Msun and M_lens2 = (3.7 +/- 0.3) x 10^10 Msun. The two lensing galaxies are likely about to undergo a dissipationless merger, and the mass and size of the resultant system should be similar to other early-type galaxies at z~0.6. This work highlights the importance of high spatial resolution imaging in developing models of strongly lensed galaxies discovered by Herschel.

preprint2012arXiv

Blind detections of CO J = 1--0 in 11 H-ATLAS galaxies at z = 2.1--3.5 with the GBT/Zpectrometer

We report measurements of the carbon monoxide ground state rotational transition (12C16O J = 1--0) with the Zpectrometer ultra-wideband spectrometer on the 100-m diameter Green Bank Telescope. The sample comprises 11 galaxies with redshifts between z = 2.1 and 3.5 from a total sample of 24 targets identified by Herschel-ATLAS photometric colors from the SPIRE instrument. Nine of the CO measurements are new redshift determinations, substantially adding to the number of detections of galaxies with rest-frame peak submillimeter emission near 100um. The CO detections confirm the existence of massive gas reservoirs within these luminous dusty star-forming galaxies (DSFGs). The CO redshift distribution of the 350um-selected galaxies is strikingly similar to the optical redshifts of 850um-selected submillimeter galaxies (SMGs) in 2.1 < z < 3.5. Spectroscopic redshifts break a temperature-redshift degeneracy; optically thin dust models fit to the far-infrared photometry indicate characteristic dust temperatures near 34 K for most of the galaxies we detect in CO. Detections of two warmer galaxies and statistically significant nondetections hint at warmer or molecule-poor DSFGs with redshifts difficult determine from from Herschel-SPIRE photometric colors alone. Many of the galaxies identified by H-ATLAS photometry are expected to be amplified by foreground gravitational lenses. Analysis of CO linewidths and luminosities provides a method for finding approximate gravitational lens magnifications mu from spectroscopic data alone, yielding mu ~ 3--20. Corrected for magnification, most galaxy luminosities are consistent with an ultra-luminous infrared galaxy (ULIRG) classification, but three are candidate hyper-LIRGs with luminosities greater than 10^13 L_sun.

preprint2012arXiv

HerMES: Candidate Gravitationally Lensed Galaxies and Lensing Statistics at Submillimeter Wavelengths

We present a list of 13 candidate gravitationally lensed submillimeter galaxies (SMGs) from 95 square degrees of the Herschel Multi-tiered Extragalactic Survey, a surface density of 0.14\pm0.04deg^{-2}. The selected sources have 500um flux densities (S_500) greater than 100mJy. Gravitational lensing is confirmed by follow-up observations in 9 of the 13 systems (70%), and the lensing status of the four remaining sources is undetermined. We also present a supplementary sample of 29 (0.31\pm0.06deg^{-2}) gravitationally lensed SMG candidates with S_500=80--100mJy, which are expected to contain a higher fraction of interlopers than the primary candidates. The number counts of the candidate lensed galaxies are consistent with a simple statistical model of the lensing rate, which uses a foreground matter distribution, the intrinsic SMG number counts, and an assumed SMG redshift distribution. The model predicts that 32--74% of our S_500>100mJy candidates are strongly gravitationally lensed (mu>2), with the brightest sources being the most robust; this is consistent with the observational data. Our statistical model also predicts that, on average, lensed galaxies with S_500=100mJy are magnified by factors of ~9, with apparently brighter galaxies having progressively higher average magnification, due to the shape of the intrinsic number counts. 65% of the sources are expected to have intrinsic 500micron flux densities less than 30mJy. Thus, samples of strongly gravitationally lensed SMGs, such as those presented here, probe below the nominal Herschel detection limit at 500 micron. They are good targets for the detailed study of the physical conditions in distant dusty, star-forming galaxies, due to the lensing magnification, which can lead to spatial resolutions of ~0.01" in the source plane.

preprint2012arXiv

Herschel-ATLAS: Multi-wavelength SEDs and physical properties of 250 micron-selected galaxies at z < 0.5

We present a pan-chromatic analysis of an unprecedented sample of 1402 250 micron-selected galaxies at z < 0.5 (mean z = 0.24) from the Herschel-ATLAS survey. We complement our Herschel 100-500 micron data with UV-K-band photometry from the Galaxy And Mass Assembly (GAMA) survey and apply the MAGPHYS energy-balance technique to produce pan-chromatic SEDs for a representative sample of 250 micron selected galaxies spanning the most recent 5 Gyr of cosmic history. We derive estimates of physical parameters, including star formation rates, stellar masses, dust masses and infrared luminosities. The typical H-ATLAS galaxy at z < 0.5 has a far-infrared luminosity in the range 10^10 - 10^12 Lsolar (SFR: 1-50 Msolar/yr) thus is broadly representative of normal star forming galaxies over this redshift range. We show that 250 micron-selected galaxies contain a larger mass of dust at a given infra-red luminosity or star formation rate than previous samples selected at 60 micron from IRAS. We derive typical SEDs for H-ATLAS galaxies, and show that the emergent SED shape is most sensitive to specific star formation rate. The optical-UV SEDs also become more reddened due to dust at higher redshifts. Our template SEDs are significantly cooler than existing infra-red templates. They may therefore be most appropriate for inferring total IR luminosities from moderate redshift submillimetre selected samples and for inclusion in models of the lower redshift submillimetre galaxy populations.

preprint2012arXiv

Herschel-ATLAS: the far-infrared properties and star-formation rates of broad absorption line quasi-stellar objects

We have used data from the Herschel-ATLAS at 250, 350 and 500 μm to determine the far-infrared (FIR) properties of 50 Broad Absorption Line Quasars (BAL QSOs). Our sample contains 49 high-ionization BAL QSOs (HiBALs) and 1 low-ionization BAL QSO (LoBAL) which are compared against a sample of 329 non-BAL QSOs. These samples are matched over the redshift range 1.5 \leq z < 2.3 and in absolute i-band magnitude over the range -28 \leq M_{i} \leq -24. Of these, 3 BAL QSOs (HiBALs) and 27 non-BAL QSOs are detected at the > 5 sigma level. We calculate star-formation rates (SFR) for our individually detected HiBAL QSOs and the non-detected LoBAL QSO as well as average SFRs for the BAL and non-BAL QSO samples based on stacking the Herschel data. We find no difference between the HiBAL and non-BAL QSO samples in the FIR, even when separated based on differing BAL QSO classifications. Using Mrk 231 as a template, the weighted mean SFR is estimated to be \approx240\pm21 M_{\odot} yr^{-1} for the full sample, although this figure should be treated as an upper limit if AGN-heated dust makes a contribution to the FIR emission. Despite tentative claims in the literature, we do not find a dependence of {\sc C\,iv} equivalent width on FIR emission, suggesting that the strength of any outflow in these objects is not linked to their FIR output. These results strongly suggest that BAL QSOs (more specifically HiBALs) can be accommodated within a simple AGN unified scheme in which our line-of-sight to the nucleus intersects outflowing material. Models in which HiBALs are caught towards the end of a period of enhanced spheroid and black-hole growth, during which a wind terminates the star-formation activity, are not supported by the observed FIR properties.

preprint2012arXiv

Herschel-ATLAS: towards a sample of ~1000 strongly-lensed galaxies

While the selection of strongly lensed galaxies with 500μm flux density S_500>100 mJy has proven to be rather straightforward (Negrello et al. 2010), for many applications it is important to analyze samples larger than the ones obtained when confining ourselves to such a bright limit. Moreover, only by probing to fainter flux densities is possible to exploit strong lensing to investigate the bulk of the high-z star-forming galaxy population. We describe HALOS (the Herschel-ATLAS Lensed Objects Selection), a method for efficiently selecting fainter candidate strongly lensed galaxies, reaching a surface density of ~1.5-2 deg^-2, i.e. a factor of about 4 to 6 higher than that at the 100 mJy flux limit. HALOS will allow the selection of up to ~1000 candidate strongly lensed galaxies (with amplifications μ>2) over the full H-ATLAS survey area. Applying HALOS to the H-ATLAS Science Demonstration Phase field (~14.4 deg^2) we find 31 candidate strongly lensed galaxies, whose candidate lenses are identified in the VIKING near-infrared catalog. Using the available information on candidate sources and candidate lenses we tentatively estimate a ~72% purity of the sample. The redshift distribution of the candidate lensed sources is close to that reported for most previous surveys for lensed galaxies, while that of candidate lenses extends to substantially higher redshifts than found in the other surveys. The counts of candidate strongly lensed galaxies are also in good agreement with model predictions (Lapi et al. 2011). Even though a key ingredient of the method is the deep near-infrared VIKING photometry, we show that H-ATLAS data alone allow the selection of a similarly deep sample of candidate strongly lensed galaxies with an efficiency close to 50%; a slightly lower surface density (~1.45 deg^-2) can be reached with a ~70% efficiency.

preprint2012arXiv

Herschel-ATLAS: VISTA VIKING near-IR counterparts in the Phase 1 GAMA 9h data

We identify near-infrared Ks band counterparts to Herschel-ATLAS sub-mm sources, using a preliminary object catalogue from the VISTA VIKING survey. The sub-mm sources are selected from the H-ATLAS Phase 1 catalogue of the GAMA 9h field, which includes all objects detected at 250, 350 or 500 um with the SPIRE instrument. We apply and discuss a likelihood ratio (LR) method for VIKING candidates within a search radius of 10" of the 22,000 SPIRE sources with a 5 sigma detection at 250 um. We find that 11,294(51%) of the SPIRE sources have a best VIKING counterpart with a reliability $R\ge 0.8$, and the false identification rate of these is estimated to be 4.2%. We expect to miss ~5% of true VIKING counterparts. There is evidence from Z-J and J-Ks colours that the reliable counterparts to SPIRE galaxies are marginally redder than the field population. We obtain photometric redshifts for ~68% of all (non-stellar) VIKING candidates with a median redshift of 0.405. Comparing to the results of the optical identifications supplied with the Phase I catalogue, we find that the use of medium-deep near-infrared data improves the identification rate of reliable counterparts from 36% to 51%.

preprint2012arXiv

Herschel-ATLAS/GAMA: a census of dust in optically selected galaxies from stacking at submillimetre wavelengths

We use the Herschel-ATLAS survey to conduct the first large-scale statistical study of the submm properties of optically selected galaxies. Using ~80,000 r-band selected galaxies from 126 deg^2 of the GAMA survey, we stack into submm imaging at 250, 350 and 500μm to gain unprecedented statistics on the dust emission from galaxies at z < 0.35. We find that low redshift galaxies account for 5% of the cosmic 250μm background (4% at 350μm; 3% at 500μm), of which approximately 60% comes from 'blue' and 20% from 'red' galaxies (rest-frame g - r). We compare the dust properties of different galaxy populations by dividing the sample into bins of optical luminosity, stellar mass, colour and redshift. In blue galaxies we find that dust temperature and luminosity correlate strongly with stellar mass at a fixed redshift, but red galaxies do not follow these correlations and overall have lower luminosities and temperatures. We make reasonable assumptions to account for the contaminating flux from lensing by red sequence galaxies and conclude that galaxies with different optical colours have fundamentally different dust emission properties. Results indicate that while blue galaxies are more luminous than red galaxies due to higher temperatures, the dust masses of the two samples are relatively similar. Dust mass is shown to correlate with stellar mass, although the dust/stellar mass ratio is much higher for low stellar mass galaxies, consistent with the lowest mass galaxies having the highest specific star formation rates. We stack the 250μm/NUV luminosity ratio, finding results consistent with greater obscuration of star formation at lower stellar mass and higher redshift. Submm luminosities and dust masses of all galaxies are shown to evolve strongly with redshift, indicating a fall in the amount of obscured star formation in ordinary galaxies over the last four billion years.

preprint2012arXiv

Herschel-ATLAS/GAMA: a difference between star-formation rates in strong-line and weak-line radio galaxies

We have constructed a sample of radio-loud objects with optical spectroscopy from the Galaxy and Mass Assembly (GAMA) project over the Herschel-ATLAS Phase 1 fields. Classifying the radio sources in terms of their optical spectra, we find that strong-emission-line sources (`high-excitation radio galaxies') have, on average, a factor ~4 higher 250-micron Herschel luminosity than weak-line (`low-excitation') radio galaxies and are also more luminous than magnitude-matched radio-quiet galaxies at the same redshift. Using all five H-ATLAS bands, we show that this difference in luminosity between the emission-line classes arises mostly from a difference in the average dust temperature; strong-emission-line sources tend to have comparable dust masses to, but higher dust temperatures than, radio galaxies with weak emission lines. We interpret this as showing that radio galaxies with strong nuclear emission lines are much more likely to be associated with star formation in their host galaxy, although there is certainly not a one-to-one relationship between star formation and strong-line AGN activity. The strong-line sources are estimated to have star-formation rates at least a factor 3-4 higher than those in the weak-line objects. Our conclusion is consistent with earlier work, generally carried out using much smaller samples, and reinforces the general picture of high-excitation radio galaxies as being located in lower-mass, less evolved host galaxies than their low-excitation counterparts.

preprint2012arXiv

Measurements of CO redshifts with Z-Spec for lensed submillimeter galaxies discovered in the H-ATLAS survey

We present new observations from Z-Spec, a broadband 185-305 GHz spectrometer, of five sub-millimeter bright lensed sources selected from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) science demonstration phase (SDP) catalog. We construct a redshift finding algorithm using combinations of the signal-to-noise of all the lines falling in the Z-Spec bandpass to determine redshifts with high confidence, even in cases where the signal-to-noise in individual lines is low. We measure the dust continuum in all sources and secure CO redshifts for four out of five (z~1.5-3). In one source, SDP.17, we tentatively identify two independent redshifts and a water line, confirmed at z=2.308. Our sources have properties characteristic of dusty starburst galaxies, with magnification-corrected star formation rates of 10^(2-3) M_sun/yr. Lower limits for the dust masses (~a few 10^8 M_sun) and spatial extents (~1 kpc equivalent radius) are derived from the continuum spectral energy distributions, corresponding to dust temperatures between 54 and 69K. In the LTE approximation, we derive relatively low CO excitation temperatures (< 100 K) and optical depths (tau<1). Performing a non-LTE excitation analysis using RADEX, we find that the CO lines measured by Z-Spec (from J=4->3 to 10->9, depending on the galaxy) localize the best solutions to either a high-temperature / low-density region, or a low-temperature / high-density region near the LTE solution, with the optical depth varying accordingly. Observations of additional CO lines, CO(1-0) in particular, are needed to constrain the non-LTE models.

preprint2012arXiv

Mining the Herschel-ATLAS: submillimeter-selected blazars in equatorial fields

The Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) provides an unprecedented opportunity to search for blazars at sub-mm wavelengths. We cross-matched the FIRST radio source catalogue with the 11655 sources brighter than 35 mJy at 500μm in the \sim 135 square degrees of the sky covered by the H-ATLAS equatorial fields at 9 h and 15 h, plus half of the field at 12 h. We found that 379 of the H-ATLAS sources have a FIRST counterpart within 10 arcsec, including 8 catalogued blazars (plus one known blazar that was found at the edge of one the H-ATLAS maps). To search for additional blazar candidates we have devised new diagnostic diagrams and found that known blazars occupy a region of the log(S500μm/S350μm) vs. log(S500μm/S1.4GHz) plane separated from that of the other sub-mm sources with radio counterparts. Using this diagnostic we have selected 12 further candidates that turn out to be scattered in the (r-z) vs. (u-r) plane or in the WISE colour-colour diagram proposed by Massaro et al. (2012), where known blazars are concentrated in well defined strips. This suggests that the majority of them either are not blazars or have spectral energy distributions contaminated by their host galaxies. A significant fraction of true blazars are found to be hosted by star-forming galaxies. This finding, supported by an analysis of blazars detected in Planck 545 and 857 GHz bands, is at odds with the notion that blazar hosts are passive ellipticals and indicates that the sub-mm selection is providing a novel prospect on blazar properties. Based on an inspection of the available photometric data, including the WISE all-sky survey, the unpublished VIKING survey and new radio observations, we tentatively estimate that there are 11 blazars with synchrotron flux density S500μm > 35mJy over the considered area. This result already allows us to constrain blazar evolution models.

preprint2012arXiv

Spitzer IRAC identification of Herschel-ATLAS SPIRE sources

We use spitzer-IRAC data to identify near-infrared counterparts to submillimeter galaxies detected with Herschel-SPIRE at 250um in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). Using a likelihood ratio analysis we identify 146 reliable IRAC counterparts to 123 SPIRE sources out of the 159. We find that, compared to the field population, the SPIRE counterparts occupy a distinct region of 3.6 and 4.5um color-magnitude space, and we use this property to identify a further 23 counterparts to 13 SPIRE sources. The IRAC identification rate of 86% is significantly higher than those that have been demonstrated with wide-field ground-based optical and near-IR imaging of Herschel fields. We estimate a false identification rate of 3.6%, corresponding to 4 to 5 sources. Among the 73 counterparts that are undetected in SDSS, 57 have both 3.6 and 4.5um coverage. Of these 43 have [3.6] - [4.5]> 0 indicating that they are likely to be at z > 1.4. Thus, ~ 40% of identified SPIRE galaxies are likely to be high redshift (z > 1.4) sources. We discuss the statistical properties of the IRAC-identified SPIRE galaxy sample including far-IR luminosities, dust temperatures, star-formation rates, and stellar masses. The majority of our detected galaxies have 10^10 to 10^11 L_sun total IR luminosities and are not intense starbursting galaxies as those found at z ~ 2, but they have a factor of 2 to 3 above average specific star-formation rates compared to near-IR selected galaxy samples.

preprint2012arXiv

The Herschel Multi-tiered Extragalactic Survey: HerMES

The Herschel Multi-tiered Extragalactic Survey, HerMES, is a legacy program designed to map a set of nested fields totalling ~380 deg^2. Fields range in size from 0.01 to ~20 deg^2, using Herschel-SPIRE (at 250, 350 and 500 μm), and Herschel-PACS (at 100 and 160 μm), with an additional wider component of 270 deg^2 with SPIRE alone. These bands cover the peak of the redshifted thermal spectral energy distribution from interstellar dust and thus capture the re-processed optical and ultra-violet radiation from star formation that has been absorbed by dust, and are critical for forming a complete multi-wavelength understanding of galaxy formation and evolution. The survey will detect of order 100,000 galaxies at 5σin some of the best studied fields in the sky. Additionally, HerMES is closely coordinated with the PACS Evolutionary Probe survey. Making maximum use of the full spectrum of ancillary data, from radio to X-ray wavelengths, it is designed to: facilitate redshift determination; rapidly identify unusual objects; and understand the relationships between thermal emission from dust and other processes. Scientific questions HerMES will be used to answer include: the total infrared emission of galaxies; the evolution of the luminosity function; the clustering properties of dusty galaxies; and the properties of populations of galaxies which lie below the confusion limit through lensing and statistical techniques. This paper defines the survey observations and data products, outlines the primary scientific goals of the HerMES team, and reviews some of the early results.

preprint2011arXiv

COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses IX. Time delays, lens dynamics and baryonic fraction in HE 0435-1223

We present accurate time delays for the quadruply imaged quasar HE 0435-1223. The delays were measured from 575 independent photometric points obtained in the R-band between January 2004 and March 2010. With seven years of data, we clearly show that quasar image A is affected by strong microlensing variations and that the time delays are best expressed relative to quasar image B. We measured Delta_t(BC) = 7.8+/-0.8 days, Delta_t(BD) = -6.5+/-0.7 days and Delta_t_CD = -14.3+/-0.8 days. We spacially deconvolved HST NICMOS2 F160W images to derive accurate astrometry of the quasar images and to infer the light profile of the lensing galaxy. We combined these images with a stellar population fitting of a deep VLT spectrum of the lensing galaxy to estimate the baryonic fraction, $f_b$, in the Einstein radius. We measured f_b = 0.65+0.13-0.10 if the lensing galaxy has a Salpeter IMF and f_b = 0.45+0.04-0.07 if it has a Kroupa IMF. The spectrum also allowed us to estimate the velocity dispersion of the lensing galaxy, sigma_ap = 222+/-34 km/s. We used f_b and sigma_ap to constrain an analytical model of the lensing galaxy composed of an Hernquist plus generalized NFW profile. We solve the Jeans equations numerically for the model and explored the parameter space under the additional requirement that the model must predict the correct astrometry for the quasar images. Given the current error bars on f_b and sigma_ap, we did not constrain H0 yet with high accuracy, i.e., we found a broad range of models with chi^2 < 1. However, narrowing this range is possible, provided a better velocity dispersion measurement becomes available. In addition, increasing the depth of the current HST imaging data of HE 0435-1223 will allow us to combine our constraints with lens reconstruction techniques that make use of the full Einstein ring that is visible in this object.

preprint2011arXiv

GAMA/H-ATLAS: The ultraviolet spectral slope and obscuration in galaxies

We use multiwavelength data from the Galaxy And Mass Assembly (GAMA) and Herschel ATLAS (H-ATLAS) surveys to compare the relationship between various dust obscuration measures in galaxies. We explore the connections between the ultraviolet (UV) spectral slope, $β$, the Balmer decrement, and the far infrared (IR) to $150\,$nm far ultraviolet (FUV) luminosity ratio. We explore trends with galaxy mass, star formation rate (SFR) and redshift in order to identify possible systematics in these various measures. We reiterate the finding of other authors that there is a large scatter between the Balmer decrement and the $β$ parameter, and that $β$ may be poorly constrained when derived from only two broad passbands in the UV. We also emphasise that FUV derived SFRs, corrected for dust obscuration using $β$, will be overestimated unless a modified relation between $β$ and the attenuation factor is used. Even in the optimum case, the resulting SFRs have a significant scatter, well over an order of magnitude. While there is a stronger correlation between the IR to FUV luminosity ratio and $β$ parameter than with the Balmer decrement, neither of these correlations are particularly tight, and dust corrections based on $β$ for high redshift galaxy SFRs must be treated with caution. We conclude with a description of the extent to which the different obscuration measures are consistent with each other as well as the effects of including other galactic properties on these correlations.

preprint2011arXiv

Gas and dust in a submillimeter galaxy at z = 4.24 from the Herschel ATLAS

We report ground-based follow-up observations of the exceptional source, ID141, one the brightest sources detected so far in the H-ATLAS cosmological survey. ID141 was observed using the IRAM 30-meter telescope and Plateau de Bure interferometer (PdBI), the Submillimeter Array (SMA) and the Atacama Pathfinder Experiment (APEX) submillimeter telescope to measure the dust continuum and emission lines of the main isotope of carbon monoxide and carbon ([C I] and [C II]). The detection of strong CO emission lines with the PdBI confirms that ID141 is at high redshift (z=4.243 +/- 0.001). The strength of the continuum and emission lines suggests that ID141 is gravitationally lensed. The width (Delta V (FWHM) ~ 800 km/s}) and asymmetric profiles of the CO and carbon lines indicate orbital motion in a disc or a merger. The properties derived for ID141 are compatible with a ultraluminous (L_FIR ~ 8.5 +/- 0.3 x 10^13/mu_L Lsun, where mu_L is the amplification factor, dense (n ~ 10^4 cm^-3) and warm (T_kin ~ 40K) starburst galaxy, with an estimated star-formation rate of (0.7 to 1.7) x 10^4/mu_L Msun/yr. The carbon emission lines indicate a dense (n ~ 10^4 cm^-3) Photo-Dominated Region, illuminated by a far-UV radiation field a few thousand times more intense than that in our Galaxy. In conclusion, the physical properties of the high-z galaxy, ID141, are remarkably similar to those of local ultraluminous infrared galaxies.

preprint2011arXiv

H-ATLAS/GAMA: Dusty early-type galaxies and passive spirals

We present the dust properties and star-formation histories of local submillimetre-selected galaxies in Herschel-ATLAS, classified by optical morphology. The early-type galaxies (ETGs) that are detected contain as much dust as typical spirals, and form a unique sample that has been blindly selected at submillimetre wavelengths. Comparing H-ATLAS galaxies to a control sample of optically selected galaxies, we find 5.5% of luminous ETGs are detected in H-ATLAS. The H-ATLAS ETGs contain a significant mass of cold dust: the mean dust mass is 5.5x10^7 Msun, with individual galaxies ranging from 9x10^5-4x10^8 Msun. This is comparable to that of spirals in our sample, and is an order of magnitude more dust than that found for the control ETGs, which have a median dust mass inferred from stacking of (0.8-4.0)x10^6 Msun. The ETGs detected in H-ATLAS have bluer NUV-r colours, higher specific star-formation rates and younger stellar populations than ETGs which are optically selected, and may be transitioning from the blue cloud to the red sequence. We also find that H-ATLAS and control ETGs inhabit similar low-density environments. We conclude that the dust in H-ATLAS and control ETGs cannot be solely from stellar sources, and a large contribution from dust formed in the ISM or external sources is required. Alternatively, dust destruction may not be as efficient as predicted. We also explore the properties of the most passive spiral galaxies in our sample with SSFR<10^-11/yr. We find these passive spirals have lower dust-to-stellar mass ratios, higher stellar masses and older stellar population ages than normal spirals. The passive spirals inhabit low density environments similar to those of the normal spiral galaxies in our sample. This shows that the processes which turn spirals passive do not occur solely in the intermediate density environments of group and cluster outskirts. (Abridged)

preprint2011arXiv

Herschel-ATLAS Galaxy Counts and High Redshift Luminosity Functions: The Formation of Massive Early Type Galaxies

Exploiting the Herschel-ATLAS Science Demonstration Phase (SDP) survey data, we have determined the luminosity functions (LFs) at rest-frame wavelengths of 100 and 250 micron and at several redshifts z>1, for bright sub-mm galaxies with star formation rates (SFR) >100 M_sun/yr. We find that the evolution of the comoving LF is strong up to z~2.5, and slows down at higher redshifts. From the LFs and the information on halo masses inferred from clustering analysis, we derived an average relation between SFR and halo mass (and its scatter). We also infer that the timescale of the main episode of dust-enshrouded star formation in massive halos (M_H>3*10^12 M_sun) amounts to ~7*10^8 yr. Given the SFRs, which are in the range 10^2-10^3 M_sun/yr, this timescale implies final stellar masses of order of 10^11-10^12 M_sun. The corresponding stellar mass function matches the observed mass function of passively evolving galaxies at z>1. The comparison of the statistics for sub-mm and UV selected galaxies suggests that the dust-free, UV bright phase, is >10^2 times shorter than the sub-mm bright phase, implying that the dust must form soon after the onset of star formation. Using a single reference Spectral Energy Distribution (SED; the one of the z~2.3 galaxy SMM J2135-0102), our simple physical model is able to reproduce not only the LFs at different redshifts > 1 but also the counts at wavelengths ranging from 250 micron to ~1 mm. Owing to the steepness of the counts and their relatively broad frequency range, this result suggests that the dispersion of sub-mm SEDs of z>1 galaxies around the reference one is rather small.

preprint2011arXiv

Herschel-ATLAS: detection of a far-infrared population around galaxy clusters

We report the detection of a significant excess in the surface density of far-infrared sources from the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) within ~1 Mpc of the centres of 66 optically-selected clusters of galaxies in the SDSS with <z>~0.25. From the analysis of the multiwavelength properties of their counterparts we conclude that the far-infrared emission is associated with dust-obscured star formation and/or active galactic nuclei within galaxies in the clusters themselves. The excess reaches a maximum at a radius of ~0.8 Mpc, where we find 1.0\pm0.3 S_250um>34 mJy sources on average per cluster above what would be expected for random field locations. If the far-infrared emission is dominated by star formation (as opposed to AGN) then this corresponds to an average star formation rate of ~7 M_sun/yr per cluster in sources with L_IR>5d10 L_sun. Although lensed sources make a negligible contribution to the excess signal, a fraction of the sources around the clusters could be gravitationally lensed, and we have identified a sample of potential cases of cluster-lensed Herschel sources that could be targeted in follow-up studies.

preprint2011arXiv

Herschel-ATLAS: First data release of the Science Demonstration Phase source catalogues

The Herschel-ATLAS is a survey of 550 square degrees with the Herschel Space Observatory in five far--infrared and submillimetre bands. The first data for the survey, observations of a field 4x4 sq. degrees in size, were taken during the Science Demonstration Phase, and reach a 5 sigma noise level of 33 mJy/beam at 250 microns. This paper describes the source extraction methods used to create the corresponding Science Demonstration Phase catalogue, which contains 6876 sources, selected at 250 microns, within ~14 sq. degrees. SPIRE sources are extracted using a new method specifically developed for Herschel data; PACS counterparts of these sources are identified using circular apertures placed at the SPIRE positions. Aperture flux densities are measured for sources identified as extended after matching to optical wavelengths. The reliability of this catalogue is also discussed, using full simulated maps at the three SPIRE bands. These show that a significant number of sources at 350 and 500 microns have undergone flux density enhancements of up to a factor of ~2, due mainly to source confusion. Correction factors are determined for these effects. The SDP dataset and corresponding catalogue will be available from http://www.h-atlas.org/.

preprint2011arXiv

Herschel-ATLAS: Rapid evolution of dust in galaxies in the last 5 billion years

We present the first direct and unbiased measurement of the evolution of the dust mass function of galaxies over the past 5 billion years of cosmic history using data from the Science Demonstration Phase of the Herschel-ATLAS. The sample consists of galaxies selected at 250μm which have reliable counterparts from SDSS at z < 0.5, and contains 1867 sources. Dust masses are calculated using both a single temperature grey-body model for the spectral energy distribution and also using a model with multiple temperature components. The dust temperature for either model shows no trend with redshift. Splitting the sample into bins of redshift reveals a strong evolution in the dust properties of the most massive galaxies. At z = 0.4 - 0.5, massive galaxies had dust masses about five times larger than in the local Universe. At the same time, the dust-to-stellar mass ratio was about 3-4 times larger, and the optical depth derived from fitting the UV-sub-mm data with an energy balance model was also higher. This increase in the dust content of massive galaxies at high redshift is difficult to explain using standard dust evolution models and requires a rapid gas consumption timescale together with either a more top-heavy IMF, efficient mantle growth, less dust destruction or combinations of all three. This evolution in dust mass is likely to be associated with a change in overall ISM mass, and points to an enhanced supply of fuel for star formation at earlier cosmic epochs.

preprint2011arXiv

Herschel-ATLAS: The link between accretion luminosity and star formation in quasar host galaxies

We use the science demonstration field data of the Herschel-ATLAS to study how star formation, traced by the far-infrared Herschel data, is related to both the accretion luminosity and redshift of quasars selected from the Sloan Digital Sky Survey and the 2SLAQ survey. By developing a maximum likelihood estimator to investigate the presence of correlations between the far-infrared and optical luminosities we find evidence that the star-formation in quasar hosts is correlated with both redshift and quasar accretion luminosity. Assuming a relationship of the form L_IR \propto L_QSO^θ (1 + z)^ζ, we find θ = 0.22 +/- 0.08 and ζ = 1.6 +/- 0.4, although there is substantial additional uncertainty in ζ of order +/- 1, due to uncertainties in the host galaxy dust temperature. We find evidence for a large intrinsic dispersion in the redshift dependence, but no evidence for intrinsic dispersion in the correlation between L_QSO and L_IR, suggesting that the latter may be due to a direct physical connection between star formation and black hole accretion. This is consistent with the idea that both the quasar activity and star formation are dependent on the same reservoir of cold gas, so that they are both affected by the influx of cold gas during mergers or heating of gas via feedback processes.

preprint2011arXiv

Observation of H2O in a strongly lensed Herschel-ATLAS source at z=2.3

The Herschel survey, H-ATLAS, with its large areal coverage, has recently discovered a number of bright, strongly lensed high-z submillimeter galaxies. The strong magnification makes it possible to study molecular species other than CO, which are otherwise difficult to observe in high-z galaxies. Among the lensed galaxies already identified by H-ATLAS, the source J090302.9-014127B (SDP.17b) at z = 2.305 is remarkable due to its excitation conditions and a tentative detection of the H2O 202-111 emission line (Lupu et al. 2010). We report observations of this line in SDP.17b using the IRAM interferometer equipped with its new 277- 371GHz receivers. The H2O line is detected at a redshift of z = 2.3049+/-0.0006, with a flux of 7.8+/-0.5 Jy km s-1 and a FWHM of 250+/-60 km s-1. The new flux is 2.4 times weaker than the previous tentative detection, although both remain marginally consistent within 1.6-sigma. The intrinsic line luminosity and ratio of H2O(202-111)/CO8-7 seem comparable with those of the nearby starburst/enshrouded-AGN Mrk 231, suggesting that SDP.17b could also host a luminous AGN. The detection of a strong H2O 202-111 line in SDP.17b implies an efficient excitation mechanism of the water levels that must occur in very dense and warm interstellar gas.

preprint2011arXiv

Physical conditions of the interstellar medium of high-redshift, strongly lensed submillimetre galaxies from the Herschel-ATLAS

We present Herschel-SPIRE Fourier Transform Spectrometer (FTS) and radio follow-up observations of two Herschel-ATLAS (H-ATLAS) detected strongly lensed distant galaxies. In one of the targeted galaxies H-ATLAS J090311.6+003906 (SDP.81) we detect [OIII] 88\mum and [CII] 158\mum lines at a signal-to-noise ratio of ~5. We do not have any positive line identification in the other fainter target H-ATLAS J091305.0-005343 (SDP.130). Currently SDP.81 is the faintest sub-mm galaxy with positive line detections with the FTS, with continuum flux just below 200 mJy in the 200-600 \mum wavelength range. The derived redshift of SDP.81 from the two detections is z=3.043 +/-0.012, in agreement with ground-based CO measurements. This is the first detection by Herschel of the [OIII] 88\mum line in a galaxy at redshift higher than 0.05. Comparing the observed lines and line ratios with a grid of photo-dissociation region (PDR) models with different physical conditions, we derive the PDR cloud density n ~ 2000 cm^{-3} and the far-UV ionizing radiation field G_0 ~ 200 (in units of the Habing field -- the local Galactic interstellar radiation field of 1.6x10^{-6} W/m^2). Using the CO derived molecular mass and the PDR properties we estimate the effective radius of the emitting region to be 500-700 pc. These characteristics are typical for star-forming, high redshift galaxies. The radio observations indicate that SDP.81 deviates significantly from the local FIR/radio correlation, which hints that some fraction of the radio emission is coming from an AGN. The constraints on the source size from millimiter-wave observations put a very conservative upper limit of the possible AGN contribution to less than 33%. These indications, together with the high [OIII]/FIR ratio and the upper limit of [OI] 63\mum/[CII] 158\mum suggest that some fraction of the ionizing radiation is likely to originate from an AGN.

preprint2011arXiv

Spitzer Imaging of Herschel-ATLAS Gravitationally Lensed Submillimeter Sources

We present physical properties of two submillimeter selected gravitationally lensed sources, identified in the Herschel Astrophysical Terahertz Large Area Survey. These submillimeter galaxies (SMGs) have flux densities > 100 mJy at 500 um, but are not visible in existing optical imaging. We fit light profiles to each component of the lensing systems in Spitzer IRAC 3.6 and 4.5 um data and successfully disentangle the foreground lens from the background source in each case, providing important constraints on the spectral energy distributions (SEDs) of the background SMG at rest-frame optical-near-infrared wavelengths. The SED fits show that these two SMGs have high dust obscuration with Av ~4 to 5 and star formation rates of ~100 M_sun/yr. They have low gas fractions and low dynamical masses compared to 850 um selected galaxies.

preprint2011arXiv

The 2010 Interim Report of the Long-Baseline Neutrino Experiment Collaboration Physics Working Groups

In early 2010, the Long-Baseline Neutrino Experiment (LBNE) science collaboration initiated a study to investigate the physics potential of the experiment with a broad set of different beam, near- and far-detector configurations. Nine initial topics were identified as scientific areas that motivate construction of a long-baseline neutrino experiment with a very large far detector. We summarize the scientific justification for each topic and the estimated performance for a set of far detector reference configurations. We report also on a study of optimized beam parameters and the physics capability of proposed Near Detector configurations. This document was presented to the collaboration in fall 2010 and updated with minor modifications in early 2011.

preprint2011arXiv

The environment and characteristics of low redshift galaxies detected by the Herschel-ATLAS

We investigate the ultraviolet and optical properties and environment of low redshift galaxies detected in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) science demonstration data. We use the Sloan Digital Sky Survey seventh release and the Galaxy And Mass Assembly database to select galaxies with r_Petro < 19.0 mag in the redshift range 0.02 < z < 0.2 and look for their submillimeter counterparts in H-ATLAS. Our results show that at low redshift, H-ATLAS detects mainly blue/star-forming galaxies with a minor contribution from red systems which are highly obscured by dust. In addition we find that the colour of a galaxy rather than the local density of its environment determines whether it is detectable by H-ATLAS. The average dust temperature of galaxies that are simultaneously detected by both PACS and SPIRE is 25K \pm 4K, independent of environment. This analysis provides a glimpse of the potential of the H-ATLAS data to investigate the submillimeter properties of galaxies in the local universe.

preprint2010arXiv

A search for debris disks in the Herschel ATLAS

Aims: We aim to demonstrate that the Herschel ATLAS (H-ATLAS) is suitable for a blind and unbiased survey for debris disks by identifying candidate debris disks associated with main sequence stars in the initial science demonstration field of the survey. We show that H-ATLAS reveals a population of far-infrared/sub-mm sources that are associated with stars or star-like objects on the SDSS main-sequence locus. We validate our approach by comparing the properties of the most likely candidate disks to those of the known population. Methods: We use a photometric selection technique to identify main sequence stars in the SDSS DR7 catalogue and a Bayesian Likelihood Ratio method to identify H-ATLAS catalogue sources associated with these main sequence stars. Following this photometric selection we apply distance cuts to identify the most likely candidate debris disks and rule out the presence of contaminating galaxies using UKIDSS LAS K-band images. Results: We identify 78 H-ATLAS sources associated with SDSS point sources on the main-sequence locus, of which two are the most likely debris disk candidates: H-ATLAS J090315.8 and H-ATLAS J090240.2. We show that they are plausible candidates by comparing their properties to the known population of debris disks. Our initial results indicate that bright debris disks are rare, with only 2 candidates identified in a search sample of 851 stars. We also show that H-ATLAS can derive useful upper limits for debris disks associated with Hipparcos stars in the field and outline the future prospects for our debris disk search programme.

preprint2010arXiv

First results from HerMES on the evolution of the submillimetre luminosity function

We have carried out two extremely deep surveys with SPIRE, one of the two cameras on Herschel, at 250 microns, close to the peak of the far-infrared background. We have used the results to investigate the evolution of the rest-frame 250-micron luminosity function out to z=2. We find evidence for strong evolution out to a redshift of around 1 but evidence for at most weak evolution beyond this redshift. Our results suggest that a significant part of the stars and metals in the Universe today were formed at z<1.4 in spiral galaxies.

preprint2010arXiv

Galaxy and Mass Assembly (GAMA): survey diagnostics and core data release

The Galaxy And Mass Assembly (GAMA) survey has been operating since February 2008 on the 3.9-m Anglo-Australian Telescope using the AAOmega fibre-fed spectrograph facility to acquire spectra with a resolution of R~1300 for 120,862 SDSS selected galaxies. The target catalogue constitutes three contiguous equatorial regions centred at 9h (G09), 12h (G12) and 14.5h (G15) each of 12 x 4 sq.deg to limiting fluxes of r < 19.4, r < 19.8, and r < 19.4 mag respectively (and additional limits at other wavelengths). Spectra and reliable redshifts have been acquired for over 98 per cent of the galaxies within these limits. Here we present the survey footprint, progression, data reduction, redshifting, re-redshifting, an assessment of data quality after 3 years, additional image analysis products (including ugrizYJHK photometry, Sersic profiles and photometric redshifts), observing mask, and construction of our core survey catalogue (GamaCore). From this we create three science ready catalogues: GamaCoreDR1 for public release, which includes data acquired during year 1 of operations within specified magnitude limits (February 2008 to April 2008); GamaCoreMainSurvey containing all data above our survey limits for use by the GAMA team and collaborators; and GamaCoreAtlasSv containing year 1, 2 and 3 data matched to Herschel-ATLAS Science Demonstration data. These catalogues along with the associated spectra, stamps and profiles can be accessed via the GAMA website: http://www.gama-survey.org/

preprint2010arXiv

Galaxy And Mass Assembly (GAMA): The input catalogue and star-galaxy separation

We describe the spectroscopic target selection for the Galaxy And Mass Assembly (GAMA) survey. The input catalogue is drawn from the Sloan Digital Sky Survey (SDSS) and UKIRT Infrared Deep Sky Survey (UKIDSS). The aim is to measure redshifts for galaxies in three 4x12 degree regions at 9h, 12h and 14.5h, on the celestial equator, with magnitude selections r < 19.4, z < 18.2 and K(AB) < 17.6 over all three regions, and r < 19.8 in the 12-h region. The target density is 1080 per sq. deg. in the 12-h region and 720 per sq. deg. in the other regions. The average GAMA target density and area are compared with completed and ongoing galaxy redshift surveys. The GAMA survey implements a highly complete star-galaxy separation that jointly uses an intensity-profile separator (delta(sg) = r-band psf mag - model mag) as per the SDSS) and a colour separator. The colour separator is defined as delta(sg,jk) = J - K - f(g-i), where f(g-i) is a quadratic fit to the J-K colour of the stellar locus over the range 0.3 < g-i < 2.3. All galaxy populations investigated are well separated with delta(sg,jk) > 0.2. From two years out of a three-year AAOmega program on the Anglo-Australian Telescope, we have obtained 79599 unique galaxy redshifts. Previously known redshifts in the GAMA region bring the total up to 98497. The median galaxy redshift is 0.2 with 99% at z < 0.5. We present some of the global statistical properties of the survey, including K-band galaxy counts, colour-redshift relations and preliminary n(z).

preprint2010arXiv

GBT Zpectrometer CO(1-0) Observations of the Strongly-Lensed Submillimeter Galaxies from the Herschel ATLAS

The Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) has uncovered a population of strongly-lensed submillimeter galaxies (SMGs). The Zpectrometer instrument on the Green Bank Telescope (GBT) was used to measure the redshifts and constrain the masses of the cold molecular gas reservoirs for two candidate high-redshift lensed sources. We derive CO(1-0) redshifts of z=3.042+/-0.001 and z=2.625+/-0.001, and measure molecular gas masses of (1--3)x10^{10}Msun, corrected for lens amplification and assuming a conversion factor of alhpa=0.8 Msun(K km/s pc^2)^{-1}. We find typical L(IR)/L'(CO) ratios of 120+/-40 and 140+/-50 Lsun (K km/s pc^2)^{-1}, which are consistent with those found for local ULIRGs and other high-redshift SMGs. From analysis of published data, we find no evidence for enhanced L(IR)/L'(CO(1-0)) ratios for the SMG population in comparison to local ULIRGs. The GBT results highlight the power of using the CO lines to derive blind redshifts, which is challenging for the SMG population at optical wavelengths given their high obscuration.

preprint2010arXiv

HerMES: Far-infrared properties of known AGN in the HerMES fields

Nuclear and starburst activity are known to often occur concomitantly. Herschel-SPIRE provides sampling of the FIR SEDs of type 1 and type 2 AGN, allowing for the separation between the hot dust (torus) and cold dust (starburst) emission. We study large samples of spectroscopically confirmed type 1 and type 2 AGN lying within the Herschel Multi-tiered Extragalactic Survey (HerMES) fields observed during the science demonstration phase, aiming to understand their FIR colour distributions and constrain their starburst contributions. We find that one third of the spectroscopically confirmed AGN in the HerMES fields have 5-sigma detections at 250um, in agreement with previous (sub)mm AGN studies. Their combined Spitzer-MIPS and Herschel-SPIRE colours - specifically S(250)/S(70) vs. S(70)/S(24) - quite clearly separate them from the non-AGN, star-forming galaxy population, as their 24-um flux is dominated by the hot torus emission. However, their SPIRE colours alone do not differ from those of non-AGN galaxies. SED fitting shows that all those AGN need a starburst component to fully account for their FIR emission. For objects at z > 2, we find a correlation between the infrared luminosity attributed to the starburst component, L(SB), and the AGN accretion luminosity, L(acc), with L(SB) propto L(acc)^0.35. Type 2 AGN detected at 250um show on average higher L(SB) than type 1 objects but their number is still too low to establish whether this trend indicates stronger star-formation activity.

preprint2010arXiv

Herschel ATLAS: The angular correlation function of submillimetre galaxies at high and low redshift

We present measurements of the angular correlation function of galaxies selected from the first field of the H-ATLAS survey. Careful removal of the background from galactic cirrus is essential, and currently dominates the uncertainty in our measurements. For our 250 micron-selected sample we detect no significant clustering, consistent with the expectation that the 250 micron-selected sources are mostly normal galaxies at z<~ 1. For our 350 micron and 500 micron-selected samples we detect relatively strong clustering with correlation amplitudes A of 0.2 and 1.2 at 1', but with relatively large uncertainties. For samples which preferentially select high redshift galaxies at z~2-3 we detect significant strong clustering, leading to an estimate of r_0 ~ 7-11 h^{-1} Mpc. The slope of our clustering measurements is very steep, delta~2. The measurements are consistent with the idea that sub-mm sources consist of a low redshift population of normal galaxies and a high redshift population of highly clustered star-bursting galaxies.

preprint2010arXiv

Herschel ATLAS: The cosmic star formation history of quasar host galaxies

We present a derivation of the star formation rate per comoving volume of quasar host galaxies, derived from stacking analyses of far-infrared to mm-wave photometry of quasars with redshifts 0<z<6 and absolute I-band magnitudes -22>I_AB>-32. We use the science demonstration observations of the first ~16 deg^2 from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) in which there are 240 quasars from the Sloan Digital Sky Survey (SDSS) and a further 171 from the 2dF-SDSS LRG and QSO (2SLAQ) survey. We supplement this data with a compilation of data from IRAS, ISO, Spitzer, SCUBA and MAMBO. H-ATLAS alone statistically detects the quasars in its survey area at >5sigma at 250, 350 and 500um. From the compilation as a whole we find striking evidence of downsizing in quasar host galaxy formation: low-luminosity quasars with absolute magnitudes in the range -22>I_AB>-24 have a comoving star formation rate (derived from 100um rest-frame luminosities) peaking between redshifts of 1 and 2, while high-luminosity quasars with I_AB<-26 have a maximum contribution to the star formation density at z~3. The volume-averaged star formation rate of -22>I_AB>-24 quasars evolves as (1+z)^{2.3 +/- 0.7} at z<2, but the evolution at higher luminosities is much faster reaching (1+z)^{10 +/- 1} at -26>I_AB>-28. We tentatively interpret this as a combination of a declining major merger rate with time and gas consumption reducing fuel for both black hole accretion and star formation.

preprint2010arXiv

Herschel-ATLAS: blazars in the SDP field

To investigate the poorly constrained sub-mm counts and spectral properties of blazars we searched for these in the Herschel-ATLAS (H-ATLAS) science demostration phase (SDP) survey catalog. We cross-matched 500$μ$m sources brighter than 50 mJy with the FIRST radio catalogue. We found two blazars, both previously known. Our study is among the first blind blazar searches at sub-mm wavelengths, i.e., in the spectral regime where little is still known about the blazar SEDs, but where the synchrotron peak of the most luminous blazars is expected to occur. Our early results are consistent with educated extrapolations of lower frequency counts and question indications of substantial spectral curvature downwards and of spectral upturns at mm wavelengths. One of the two blazars is identified with a Fermi/LAT $γ$-ray source and a WMAP source. The physical parameters of the two blazars are briefly discussed.These observations demonstrate that the H-ATLAS survey will provide key information about the physics of blazars and their contribution to sub-mm counts.

preprint2010arXiv

Herschel-ATLAS: Dust temperature and redshift distribution of SPIRE and PACS detected sources using submillimetre colours

We present colour-colour diagrams of detected sources in the Herschel-ATLAS Science Demonstration Field from 100 to 500 microns using both PACS and SPIRE. We fit isothermal modified black bodies to the spectral energy distribution (SED) to extract the dust temperature of sources with counterparts in Galaxy And Mass Assembly (GAMA) or SDSS surveys with either a spectroscopic or a photometric redshift. For a subsample of 330 sources detected in at least three FIR bands with a significance greater than 3 $σ$, we find an average dust temperature of $(28 \pm 8)$K. For sources with no known redshift, we populate the colour-colour diagram with a large number of SEDs generated with a broad range of dust temperatures and emissivity parameters, and compare to colours of observed sources to establish the redshift distribution of this sample. For another subsample of 1686 sources with fluxes above 35 mJy at 350 microns and detected at 250 and 500 microns with a significance greater than 3$σ$, we find an average redshift of $2.2 \pm 0.6$.

preprint2010arXiv

Herschel-ATLAS: Statistical Properties of Galactic Cirrus in the GAMA-9 Hour Science Demonstration Phase Field

We study the Spectral Energy Distribution (SED) and the power spectrum of Galactic cirrus emission observed in the 14 deg^2 Science Demonstration Phase field of the Herschel-ATLAS using Herschel and IRAS data from 100 to 500 um. We compare the SPIRE 250, 350 and 500um maps with IRAS 100um emission, binned in 6' pixels. We assume a modified black-body SED with dust emissivity parameter beta (F ~ lambda^(-beta)) and a single dust temperature T_d, and find that the dust temperature and emissivity index varies over the science demonstration field as 10< T_rm < 25 K and 1 < beta< 4. The latter values are somewhat higher than the range of beta often quoted in the literature (1< beta< 2). We estimate the mean values of these parameters to be T_d=19.0 +/- 2.4 K and beta = 1.4 +/- 0.4. In regions of bright cirrus emission, we find that the dust has similar temperatures with T_d = 18.0 +/- 2.5 K, and similar values of beta, ranging from 1.4 +- 0.5 to 1.9+/- 0.5. We show that T_d and beta associated with diffuse cirrus emission are anti-correlated and can be described by the relationship: beta(T_d) = NT_d^alpha with [N=116+/-38, alpha=-1.4+/1 0.1]. The strong correlation found in this analysis is not just limited to high density clumps of cirrus emission as seen in previous studies, but is also seen in diffuse cirrus in low density regions. To provide an independent measure of $T_{\rm d}$ and $β$, we obtain the angular power spectrum of the cirrus emission in the {\it IRAS} and SPIRE maps, which is consistent with a power spectrum of the form P(k)=P_0(k/k_0)^gamma where gamma = â H R2.6+/-m 0.2 for scales of 50-200' in the SPIRE maps. The cirrus rms fluctuation amplitude at angular scales of 100' is consistent with a modified blackbody SED with T_d = 20.1+/- 0.9 K and beta = 1.3+/- 0.2, in agreement with the values obtained above.

preprint2010arXiv

Herschel-ATLAS: The far-infrared--radio correlation at z<0.5

We use data from the Herschel-ATLAS to investigate the evolution of the far-infrared--radio correlation over the redshift range 0<z<0.5. Using the total far-infrared luminosity of all >5sigma sources in the Herschel-ATLAS Science Demonstration Field and cross-matching these data with radio data from the Faint Images of the Radio Sky at Twenty-Centimetres (FIRST) survey and the NRAO VLA Northern Sky Survey (NVSS), we obtain 104 radio counterparts to the Herschel sources. With these data we find no evidence for evolution in the far-infrared--radio correlation over the redshift range 0<z<0.5, where the median value for the ratio between far-infrared and radio luminosity, $q_{\rm IR}$, over this range is $q_{\rm IR} = 2.40\pm 0.12$ (and a mean of $q_{\rm IR}=2.52 \pm 0.03$ accounting for the lower limits), consistent with both the local value determined from {\em IRAS} and values derived from surveys targeting the high-redshift Universe. By comparing the radio fluxes of our sample measured from both FIRST and NVSS we show that previous results suggesting an increase in the value of $q_{\rm IR}$ from high to low redshift may be the result of resolving out extended emission of the low-redshift sources with relatively high-resolution interferometric data, although AGN contamination could still play a significant role. We also find tentative evidence that the longer wavelength, cooler dust is heated by an evolved stellar population which does not trace the star-formation rate as closely as the shorter wavelength $\ltsim 250~μ$m emission or the radio emission, supporting suggestions based on detailed models of individual galaxies.

preprint2010arXiv

Mid-infrared spectroscopy of candidate AGN-dominated submillimeter galaxies

Spitzer spectroscopy has revealed that ~80% of submm galaxies (SMGs) are starburst (SB) dominated in the mid-infrared. Here we focus on the remaining ~20% that show signs of harboring powerful active galactic nuclei (AGN). We have obtained Spitzer-IRS spectroscopy of a sample of eight SMGs which are candidates for harboring powerful AGN on the basis of IRAC color-selection (S8/S4.5>2; i.e. likely power-law mid-infrared SEDs). SMGs with an AGN dominating (>50%) their mid-infrared emission could represent `missing link' sources in an evolutionary sequence involving a major merger. First of all, we detect PAH features in all of the SMGs, indicating redshifts from 2.5-3.4, demonstrating the power of the mid-infrared to determine redshifts for these optically faint dusty galaxies. Secondly, we see signs of both star-formation (from the PAH features) and AGN activity (from continuum emission) in our sample: 62% of the sample are AGN-dominated in the mid-infrared with a median AGN content of 56%, compared with <30% on average for typical SMGs, revealing that our IRAC color selection has successfully singled out sources with proportionately more AGN emission than typical SB-dominated SMGs. However, we find that only about 10% of these AGN dominate the bolometric emission of the SMG when the results are extrapolated to longer infrared wavelengths, implying that AGN are not a significant power source to the SMG population overall, even when there is evidence in the mid-infrared for substantial AGN activity. When existing samples of mid-infrared AGN-dominated SMGs are considered, we find that S8/S4.5>1.65 works well at selecting mid-infrared energetically dominant AGN in SMGs, implying a duty cycle of ~15% if all SMGs go through a subsequent mid-infrared AGN-dominated phase in the proposed evolutionary sequence.

preprint2010arXiv

The AGN fraction of submm-selected galaxies and contributions to the submm/mm-wave extragalactic background light

We present a comparison of the SCUBA Half Degree Extragalactic Survey (SHADES) at 450, 850 and 1100 microns with deep guaranteed time 15 microns AKARI FU-HYU survey data and Spitzer guaranteed time data at 3.6-24 microns in the Lockman Hole East. The AKARI data was analysed using bespoke software based in part on the drizzling and minimum-variance matched filtering developed for SHADES, and was cross calibrated against Infrared Space Observatory (ISO) fluxes. Our stacking analyses find AKARI 15um galaxies with >~200 microJy contribute >10% of the 450 micron background, but only <4% of the 1100 micron background, suggesting that different populations contribute at mm-wavelengths. We confirm our earlier result that the ultra-deep 450 micron SCUBA-2 Cosmology Survey will be dominated by populations already detected by AKARI and Spitzer mid-infrared surveys. The superb mid-infrared wavelength coverage afforded by combining Spitzer and AKARI photometry is an excellent diagnostic of AGN contributions, and we find that (23-52)% of submm-selected galaxies have AGN bolometric fractions f_AGN>0.3.

preprint2010arXiv

The Detection of a Population of Submillimeter-Bright, Strongly-Lensed Galaxies

Gravitational lensing is a powerful astrophysical and cosmological probe and is particularly valuable at submillimeter wavelengths for the study of the statistical and individual properties of dusty starforming galaxies. However the identification of gravitational lenses is often time-intensive, involving the sifting of large volumes of imaging or spectroscopic data to find few candidates. We used early data from the Herschel Astrophysical Terahertz Large Area Survey to demonstrate that wide-area submillimeter surveys can simply and easily detect strong gravitational lensing events, with close to 100% efficiency.

preprint2010arXiv

The first release of data from the Herschel ATLAS: the SPIRE images

We have reduced the data taken with the Spectral and Photometric Imaging Receiver (SPIRE) photometer on board the Herschel Space Observatory in the Science Demonstration Phase (SDP) of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). We describe the data reduction, which poses specific challenges, both because of the sheer size of the data, and because only two scans are made for each region. We implement effective solutions to process the bolometric timelines into maps, and show that correlations among detectors are negligible, and that the photometer is stable on time scales up to 250 s. This is longer than the time the telescope takes to cross the observed sky region, and it allows us to use naive binning methods for an optimal reconstruction of the sky emission. The maps have equal contribution of confusion and white instrumental noise, and the latter is estimated to 5.3, 6.4, and 6.7 mJy/beam (1-σ), at 250, 350, and 500 μ{m}, respectively. This pipeline is used to reduce other H-ATLAS observations, as they became available, and we discuss how it can be used with the optimal map maker implemented in the Herschel Interactive Processing Environment (HIPE), to improve computational efficiency and stability. The SDP dataset is available from http://www.h-atlas.org/.

preprint2010arXiv

The HerMES SPIRE submillimeter local luminosity function

Local luminosity functions are fundamental benchmarks for high-redshift galaxy formation and evolution studies as well as for models describing these processes. Determining the local luminosity function in the submillimeter range can help to better constrain in particular the bolometric luminosity density in the local Universe, and Herschel offers the first opportunity to do so in an unbiased way by imaging large sky areas at several submillimeter wavelengths. We present the first Herschel measurement of the submillimeter 0<z<0.2 local luminosity function and infrared bolometric (8-1000 $μ$m) local luminosity density based on SPIRE data from the HerMES Herschel Key Program over 14.7 deg^2. Flux measurements in the three SPIRE channels at 250, 350 and 500 \mum are combined with Spitzer photometry and archival data. We fit the observed optical-to-submillimeter spectral energy distribution of SPIRE sources and use the 1/V_{max} estimator to provide the first constraints on the monochromatic 250, 350 and 500 \mum as well as on the infrared bolometric (8-1000 \mum) local luminosity function based on Herschel data. We compare our results with modeling predictions and find a slightly more abundant local submillimeter population than predicted by a number of models. Our measurement of the infrared bolometric (8-1000 \mum) local luminosity function suggests a flat slope at low luminosity, and the inferred local luminosity density, 1.31_-0.21^+0.24 x 10^8 Lsun Mpc^-3, is consistent with the range of values reported in recent literature.

preprint2010arXiv

The Herschel Multi-Tiered Extragalactic Survey: Source Extraction and Cross-Identifications in Confusion-Dominated SPIRE Images

We present the cross-identification and source photometry techniques used to process Herschel SPIRE imaging taken as part of the Herschel Multi-Tiered Extragalactic Survey (HerMES). Cross-identifications are performed in map-space so as to minimise source blending effects. We make use of a combination of linear inversion and model selection techniques to produce reliable cross-identification catalogues based on Spitzer MIPS 24 micron source positions. Testing on simulations and real Herschel observations show that this approach gives robust results for even the faintest sources S250~10 mJy. We apply our new technique to HerMES SPIRE observations taken as part of the science demostration phase of Herschel. For our real SPIRE observations we show that, for bright unconfused sources, our flux density estimates are in good agreement with those produced via more traditional point source detection methods (SussExtractor; Savage & Oliver et al. 2006) by Smith et al. 2010. When compared to the measured number density of sources in the SPIRE bands, we show that our method allows the recovery of a larger fraction of faint sources than these traditional methods. However this completeness is heavily dependant on the relative depth of the existing 24 micron catalogues and SPIRE imaging. Using our deepest multi-wavelength dataset in GOODS-N, we estimate that the use of shallow 24 micron in our other fields introduces an incompleteness at faint levels of between 20-40 per cent at 250 micron.

preprint2010arXiv

The Herschel-ATLAS: Evolution of the 250um luminosity function out to z=0.5

We have determined the luminosity function of 250um-selected galaxies detected in the ~14 sq.deg science demonstration region of the Herschel-ATLAS project out to a redshift of z=0.5. Our findings very clearly show that the luminosity function evolves steadily out to this redshift. By selecting a sub-group of sources within a fixed luminosity interval where incompleteness effects are minimal, we have measured a smooth increase in the comoving 250um luminosity density out to z=0.2 where it is 3.6+1.4-0.9 times higher than the local value.

preprint2010arXiv

The Herschel-ATLAS: Extragalactic Number Counts from 250 to 500 Microns

Aims. The Herschel-ATLAS survey (H-ATLAS) will be the largest area survey to be undertaken by the Herschel Space Observatory. It will cover 550 sq. deg. of extragalactic sky at wavelengths of 100, 160, 250, 350 and 500 microns when completed, reaching flux limits (5 sigma) from 32 to 145mJy. We here present galaxy number counts obtained for SPIRE observations of the first ~14 sq. deg. observed at 250, 350 and 500 microns. Methods. Number counts are a fundamental tool in constraining models of galaxy evolution. We use source catalogs extracted from the H-ATLAS maps as the basis for such an analysis. Correction factors for completeness and flux boosting are derived by applying our extraction method to model catalogs and then applied to the raw observational counts. Results. We find a steep rise in the number counts at flux levels of 100-200mJy in all three SPIRE bands, consistent with results from BLAST. The counts are compared to a range of galaxy evolution models. None of the current models is an ideal fit to the data but all ascribe the steep rise to a population of luminous, rapidly evolving dusty galaxies at moderate to high redshift.

preprint2010arXiv

The Herschel-ATLAS: The dust energy balance in the edge-on spiral galaxy UGC 4754

We use Herschel PACS and SPIRE observations of the edge-on spiral galaxy UGC 4754, taken as part of the H-ATLAS SDP observations, to investigate the dust energy balance in this galaxy. We build detailed SKIRT radiative models based on SDSS and UKIDSS maps and use these models to predict the far-infrared emission. We find that our radiative transfer model underestimates the observed FIR emission by a factor of two to three. Similar discrepancies have been found for other edge-on spiral galaxies based on IRAS, ISO, and SCUBA data. Thanks to the good sampling of the SED at FIR wavelengths, we can rule out an underestimation of the FIR emissivity as the cause for this discrepancy. Instead we support highly obscured star formation that contributes little to the optical extinction as a more probable explanation.

preprint2010arXiv

Which halos host Herschel-ATLAS galaxies in the local Universe?

We measure the projected cross-correlation between low redshift (z < 0.5) far-IR selected galaxies in the SDP field of the Herschel-ATLAS (H-ATLAS) survey and optically selected galaxies from the Galaxy and Mass Assembly (GAMA) redshift survey. In order to obtain robust correlation functions, we restrict the analysis to a subset of 969 out of 6900 H-ATLAS galaxies, which have reliable optical counterparts with r<19.4 mag and well-determined spectroscopic redshifts. The overlap region between the two surveys is 12.6 sq. deg; the matched sample has a median redshift of z ~ 0.2. The cross-correlation of GAMA and H-ATLAS galaxies within this region can be fitted by a power law, with correlation length r_0 ~ 4.63 +/- 0.51 Mpc. Comparing with the corresponding auto-correlation function of GAMA galaxies within the SDP field yields a relative bias (averaged over 2-8 Mpc) of H-ATLAS and GAMA galaxies of b_H/b_G ~ 0.6. Combined with clustering measurements from previous optical studies, this indicates that most of the low redshift H-ATLAS sources are hosted by halos with masses comparable to that of the Milky Way. The correlation function appears to depend on the 250 um luminosity, L_250, with bright (median luminosity νL_250 ~ 1.6 x 10^10 L_sun) objects being somewhat more strongly clustered than faint (νL_250 ~ 4.0 x 10^9 L_sun) objects. This implies that galaxies with higher dust-obscured star formation rates are hosted by more massive halos.

preprint2009arXiv

AzTEC Half Square Degree Survey of the SHADES Fields -- I. Maps, Catalogues, and Source Counts

We present the first results from the largest deep extragalactic millimetre-wavelength survey undertaken to date. These results are derived from maps covering over 0.7 deg^2, made at 1.1mm, using the AzTEC continuum camera mounted on the James Clerk Maxwell Telescope. The maps were made in the two fields originally targeted at 0.85mm with SCUBA in the SHADES project, namely the Lockman Hole East (mapped to a depth of 0.9-1.3 mJy rms) and the Subaru XMM Deep Field (1.0-1.7 mJy rms). The wealth of existing and forthcoming deep multi-frequency data in these two fields will allow the bright mm source population revealed by these images to be explored in detail in subsequent papers. Here we present the maps themselves, a catalogue of 114 high-significance sub-millimetre galaxy detections, and a thorough statistical analysis leading to the most robust determination to date of the 1.1mm source number counts. Through careful comparison, we find that both the COSMOS and GOODS North fields, also imaged with AzTEC, contain an excess of mm sources over the new 1.1mm source-count baseline established here. In particular, our new AzTEC/SHADES results indicate that very luminous high-redshift dust enshrouded starbursts (S_{1.1} > 3 mJy) are 25-50% less common than would have been inferred from these smaller surveys, thus highlighting the potential roles of cosmic variance and clustering in such measurements. We compare number count predictions from recent models of the evolving mm/sub-mm source population to these SMG surveys, which provide important constraints for the ongoing refinement of semi-analytic and hydrodynamical models of galaxy formation, and find that all recent models over-predict the number of bright sub-millimetre galaxies found in this survey.

preprint2009arXiv

Discovery of a redshift 6.13 quasar in the UKIRT Infrared Deep Sky Survey

Optical and near-infrared (NIR) spectra are presented for ULAS J131911.29+095951.4 (hereafter ULAS J1319+0950), a new redshift z = 6.127 +/- 0.004 quasar discovered in the Third Data Release (DR3) of the UKIRT Infrared Deep Sky Survey (UKIDSS). The source has Y_Vega = 19.10 +/- 0.03, corresponding to M_1450,AB = -27.12, which is comparable to the absolute magnitudes of the z ~= 6 quasars discovered in the SDSS. ULAS J1319+0950 was, in fact, registered by SDSS as a faint source with z_AB = 20.1 +/- 0.1, just below the signal--to--noise ratio limit of the SDSS high-redshift quasar survey. The faint z-band magnitude is a consequence of the weak Ly alpha/N V emission line, which has a rest-frame equivalent width of only ~20 A and provides only a small boost to the z-band flux. Nevertheless, there is no evidence from this UKIDSS-based search for a significant population of high-redshift quasars with weak emission lines. The Ly alpha optical depth to ULAS J1319+0950 is consistent with that measured towards similarly distant SDSS quasars, implying that results from optical- and NIR-selected quasars may be combined. Also presented is a new NIR-spectrum of the previously discovered UKIDSS quasar ULAS J020332.38+001229.2, which reveals the object to be a broad absorption line quasar. The new spectrum shows that the emission line previously identified as Ly alpha is actually N V, leading to a revised redshift of z = 5.72, rather than z = 5.86 as previously estimated.

preprint2009arXiv

The Herschel ATLAS

The Herschel ATLAS is the largest open-time key project that will be carried out on the Herschel Space Observatory. It will survey 510 square degrees of the extragalactic sky, four times larger than all the other Herschel surveys combined, in five far-infrared and submillimetre bands. We describe the survey, the complementary multi-wavelength datasets that will be combined with the Herschel data, and the six major science programmes we are undertaking. Using new models based on a previous submillimetre survey of galaxies, we present predictions of the properties of the ATLAS sources in other wavebands.

preprint2007arXiv

The UKIRT Infrared Deep Sky Survey First Data Release

The First Data Release (DR1) of the UKIRT Infrared Deep Sky Survey (UKIDSS) took place on 2006 July 21. UKIDSS is a set of five large near-infrared surveys, covering a complementary range of areas, depths, and Galactic latitudes. DR1 is the first large release of survey-quality data from UKIDSS and includes 320 sq degs of multicolour data to (Vega) K=18, complete (depending on the survey) in three to five bands from the set ZYJHK, together with 4 sq degs of deep JK data to an average depth K=21. In addition the release includes a similar quantity of data with incomplete filter coverage. In JHK, in regions of low extinction, the photometric uniformity of the calibration is better than 0.02 mag. in each band. The accuracy of the calibration in ZY remains to be quantified, and the same is true of JHK in regions of high extinction. The median image FWHM across the dataset is 0.82 arcsec. We describe changes since the Early Data Release in the implementation, pipeline and calibration, quality control, and archive procedures. We provide maps of the areas surveyed, and summarise the contents of each of the five surveys in terms of filters, areas, and depths. DR1 marks completion of 7 per cent of the UKIDSS 7-year goals.

preprint2003arXiv

The evolution of faint AGN between z~1 and z~5 from the COMBO-17 survey

We present a determination of the optical/UV AGN luminosity function and its evolution, based on a large sample of faint (R<24) QSOs identified in the COMBO-17 survey. Using multi-band photometry in 17 filters within 350nm < lambda_obs < 930nm, we could simultaneously determine photometric redshifts with an accuracy of sigma_z<0.03 and obtain spectral energy distributions. The redshift range covered by the sample is 1.2<z<4.8, which implies that even at z~3, the sample reaches below luminosities corresponding to M_B = -23, conventionally employed to distinguish between Seyfert galaxies and quasars. We clearly detect a broad plateau-like maximum of quasar activity around z~2 and map out the smooth turnover between z~1 and z~4. The shape of the LF is characterised by some mild curvature, but no sharp `break' is present within the range of luminosities covered. Using only the COMBO-17 data, the evolving LF can be adequately described by either a pure density evolution (PDE) or a pure luminosity evolution (PLE) model. However, the absence of a strong L*-like feature in the shape of the LF inhibits a robust distinction between these modes. We present a robust estimate for the integrated UV luminosity generation by AGN as a function of redshift. We find that the LF continues to rise even at the lowest luminosities probed by our survey, but that the slope is sufficiently shallow that the contribution of low-luminosity AGN to the UV luminosity density is negligible. Although our sample reaches much fainter flux levels than previous data sets, our results on space densities and LF slopes are completely consistent with extrapolations from recent major surveys such as SDSS and 2QZ.