Source author record

Douglas Scott

Douglas Scott 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

160works
9topics
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

160 published item(s)

preprint2024arXiv

Cosmic Microwave Background Mini-review

A compact overview of the status of CMB anisotropy results and their cosmological interpretation. Section headings: Introduction; CMB Spectrum; Description of CMB Anisotropies; Cosmological Parameters; Physics of Anisotropies; Current Anisotropy Data; CMB Polarization; Complications; Constraints on Cosmological Parameters; Particle Physics Constraints; Fundamental Lessons; and Future Directions. This 2023 replacement is a substantial update compared with earlier versions of this review posted to the arXiv.

preprint2023arXiv

Brightest Cluster Galaxy Formation in the z=4.3 Protocluster SPT2349-56: Discovery of a Radio-Loud AGN

We have observed the z=4.3 protocluster SPT2349-56 with ATCA with the aim of detecting radio-loud active galactic nuclei (AGN) amongst the ~30 submillimeter galaxies identified in the structure. We detect the central complex of SMGs at 2.2\,GHz with a luminosity of L_2.2=(4.42pm0.56)x10^{25} W/Hz. The ASKAP also detects the source at 888 MHz, constraining the radio spectral index to alpha=-1.6pm0.3, consistent with ATCA non-detections at 5.5 and 9GHz, and implying L_1.4(rest)=(2.4pm0.3)x10^{26}W/Hz. This radio luminosity is about 100 times higher than expected from star formation, assuming the usual FIR-radio correlation, which is a clear indication of an AGN driven by a forming brightest cluster galaxy (BCG). None of the SMGs in SPT2349-56 show signs of AGN in any other diagnostics available to us (notably 12CO out to J=16, OH163um, CII/IR, and optical spectra), highlighting the radio continuum as a powerful probe of obscured AGN in high-z protoclusters. No other significant radio detections are found amongst the cluster members, consistent with the FIR-radio correlation. We compare these results to field samples of radio sources and SMGs, along with the 22 SPT-SMG gravitational lenses also observed in the ATCA program, as well as powerful radio galaxies at high redshifts. Our results allow us to better understand the effects of this gas-rich, overdense environment on early supermassive black hole (SMBH) growth and cluster feedback. We estimate that (3.3pm0.7)x10^{38} W of power are injected into the growing ICM by the radio-loud AGN, whose energy over 100Myr is comparable to the binding energy of the gas mass of the central halo. The AGN power is also comparable to the instantaneous energy injection from supernova feedback from the 23 catalogued SMGs in the core region of 120kpc projected radius. The SPT2349-56 radio-loud AGN may be providing strong feedback on a nascent ICM.

preprint2022arXiv

Candidate high-redshift protoclusters and lensed galaxies in the Planck list of high-z sources overlapping with Herschel-SPIRE imaging

The Planck list of high-redshift source candidates (the PHz catalogue) contains 2151 peaks in the cosmic infrared background, unresolved by Planck's 5 arcmin beam. Follow-up spectroscopic observations have revealed that some of these objects are $z \approx 2$ protoclusters and strong gravitational lenses, but an unbiased survey has not yet been carried out. To this end, we have used archival Herschel-SPIRE observations to study a uniformly-selected sample of 187 PHz sources. In contrast with follow-up studies that were biased towards bright, compact sources, we find that only one of our PHz sources is a bright gravitationally-lensed galaxy (peak flux $\gtrsim 300$ mJy), indicating that such objects are rarer in the PHz catalogue than previously believed ($< 1$ per cent). The majority of our PHz sources consist of many red, star-forming galaxies, demonstrating that typical PHz sources are candidate protoclusters. However, our new PHz sources are significantly less bright than found in previous studies and differ in colour, suggesting possible differences in redshift and star-formation rate. Nonetheless, 40 of our PHz sources contain $> 3σ$ galaxy overdensities, comparable to the fraction of $> 3σ$ overdensities found in earlier biased studies. We additionally use a machine-learning approach to identify less extreme (peak flux $\sim 100$ mJy) gravitationally-lensed galaxies among Herschel-SPIRE observations of PHz sources, finding a total of seven candidates in our unbiased sample, and 13 amongst previous biased samples. Our new uniformly-selected catalogues of $> 3σ$ candidate protoclusters and strong gravitational lenses provide interesting targets for follow up with higher-resolution facilities, such as ALMA and JWST.

preprint2022arXiv

Constraints on cosmic birefringence using $E$-mode polarisation

A birefringent universe could show itself through a rotation of the plane of polarisation of the cosmic microwave background photons. This is usually investigated using polarisation $B$ modes, which is degenerate with miscalibration of the orientation of the polarimeters. Here we point out an independent method for extracting the birefringence angle using only temperature and $E$-mode signals. We forecast that, with an ideal cosmic-variance-limited experiment, we could constrain a birefringence angle of $0.3^\circ$ with $3\,σ$ statistical significance, which is close to the current constraints using $B$ modes. We explore how this method is affected by the systematic errors introduced by the polarisation efficiency. In the future, this could provide an additional way of checking any claimed $B$-mode derived birefringence signature.

preprint2022arXiv

Rapid build-up of the stellar content in the protocluster core SPT2349$-$56 at $z\,{=}\,4.3$

The protocluster SPT2349$-$56 at $z\,{=}\,4.3$ contains one of the most actively star-forming cores known, yet constraints on the total stellar mass of this system are highly uncertain. We have therefore carried out deep optical and infrared observations of this system, probing rest-frame ultraviolet to infrared wavelengths. Using the positions of the spectroscopically-confirmed protocluster members, we identify counterparts and perform detailed source deblending, allowing us to fit spectral energy distributions in order to estimate stellar masses. We show that the galaxies in SPT2349$-$56 have stellar masses proportional to their high star-formation rates, consistent with other protocluster galaxies and field submillimetre galaxies (SMGs) around redshift 4. The galaxies in SPT2349$-$56 have on average lower molecular gas-to-stellar mass fractions and depletion timescales than field SMGs, although with considerable scatter. We construct the stellar-mass function for SPT2349$-$56 and compare it to the stellar-mass function of $z\,{=}\,1$ galaxy clusters, finding consistent shapes between the two. We measure rest-frame galaxy ultraviolet half-light radii from our HST-F160W imaging, finding that on average the galaxies in our sample are similar in size to typical star-forming galaxies at these redshifts. However, the brightest HST-detected galaxy in our sample, found near the luminosity-weighted centre of the protocluster core, remains unresolved at this wavelength. Hydrodynamical simulations predict that the core galaxies will quickly merge into a brightest cluster galaxy, thus our observations provide a direct view of the early formation mechanisms of this class of object.

preprint2020arXiv

Megaparsec-scale structure around the proto-cluster core SPT2349$-$56 at $z\,{=}\,4.3$

We present an extensive ALMA spectroscopic follow-up programme of the $z\,{=}\,4.3$ structure SPT2349$-$56, one of the most actively star-forming proto-cluster cores known, to identify additional members using their [C{\sc ii}] 158\,$μ$m and \mbox{CO(4--3)} lines. In addition to robustly detecting the 14 previously published galaxies in this structure, we identify a further 15 associated galaxies at $z\,{=}\,4.3$, resolving 55$\,{\pm}\,$5\,per cent of the 870-$μ$m flux density at 0.5\,arcsec resolution compared to 21\,arcsec single-dish data. These galaxies are distributed into a central core containing 23 galaxies extending out to 300\,kpc in diameter, and a northern extension, offset from the core by 400\,kpc, containing three galaxies. We discovered three additional galaxies in a red {\it Herschel\/}-SPIRE source 1.5\,Mpc from the main structure, suggesting the existence of many other sources at the same redshift as SPT2349$-$56 that are not yet detected in the limited coverage of our data. An analysis of the velocity distribution of the central galaxies indicates that this region may be virialized with a mass of (9$\pm$5)$\,{\times}\,$10$^{12}$\,M$_{\odot}$, while the two offset galaxy groups are about 30 and 60\,per cent less massive and show significant velocity offsets from the central group. We calculate the [C{\sc ii}] and far-infrared number counts, and find evidence for a break in the [C{\sc ii}] luminosity function. We estimate the average SFR density within the region of SPT2349$-$56 containing single-dish emission (a proper diametre of 720\,kpc), assuming spherical symmetry, to be roughly 4$\,{\times}\,10^4$\,M$_{\odot}$\,yr$^{-1}$\,Mpc$^{-3}$; this may be an order of magnitude greater than the most extreme examples seen in simulations.

preprint2020arXiv

Re-evaluating evidence for Hawking points in the CMB

We investigate recent claims for a detection of "Hawking points" (positions on the sky with unusually large temperature gradients between rings) in the cosmic microwave background (CMB) temperature maps at the 99.98% confidence level. We find that, after marginalization over the size of the rings, an excess is detected in Planck satellite maps at only an 87% confidence level (i.e., little more than 1$σ$). Therefore, we conclude that there is no statistically significant evidence for the presence of Hawking points in the CMB.

preprint2020arXiv

Rise of the Titans: Gas Excitation and Feedback in a Binary Hyper-Luminous Dusty Starburst Galaxy at z~6

We report new observations toward the hyper-luminous dusty starbursting major merger ADFS-27 (z=5.655), using ATCA and ALMA. We detect CO 2-1, 8-7, 9-8, 10-9 and H2O(321-221) emission, and a P-Cygni-shaped OH+(11-01) absorption/emission feature. We also tentatively detect H2O(321-312) and OH+(12-01) emission and CH+(1-0) absorption. We find a total cold molecular mass of M_gas = (2.1+/-0.2) x 10^11 (alpha_CO/1.0) Msun. We also find that the excitation of the star-forming gas is overall moderate for a z>5 dusty starburst, which is consistent with its moderate dust temperature. A high density, high kinetic temperature gas component embedded in the gas reservoir is required to fully explain the CO line ladder. This component is likely associated with the "maximum starburst" nuclei in the two merging galaxies, which are separated by only (140+/-13) km/s along the line of sight and 9.0 kpc in projection. The kinematic structure of both components is consistent with galaxy disks, but this interpretation remains limited by the spatial resolution of the current data. The OH+ features are only detected towards the northern component, which is also the one that is more enshrouded in dust and thus remains undetected up to 1.6 um even in our sensitive new HST/WFC3 imaging. The absorption component of the OH+ line is blueshifted and peaks near the CO and continuum emission peak while the emission is redshifted and peaks offset by 1.7 kpc from the CO and continuum emission peak, suggesting that the gas is associated with a massive molecular outflow from the intensely star-forming nucleus that supplies 125 Msun/yr of enriched gas to its halo.

preprint2020arXiv

SCUBA-2 Ultra Deep Imaging Eao Survey (Studies) III: Multi-wavelength properties, luminosity functions and preliminary source catalog of 450-$μ$m-selected galaxies

We construct a SCUBA-2 450-$μ$m map in the COSMOS field that covers an area of 300 arcmin$^{2}$ and reaches a 1$σ$ noise level of 0.65 mJy in the deepest region. We extract 256 sources detected at 450 $μ$m with signal-to-noise ratio $>$ 4.0 and analyze the physical properties of their multi-wavelength counterparts. We find that most of the sources are at $z\lesssim3$, with a median of $z = 1.79^{+0.03}_{-0.15}$. About $35^{+32}_{-25}$% of our sources are classified as starburst galaxies based on their total star-formation rates (SFRs) and stellar masses ($M_{\ast}$). By fitting the far-infrared spectral energy distributions, we find that our 450-$μ$m-selected sample has a wide range of dust temperatures (20 K $ \lesssim T_{\rm d} \lesssim$ 60 K), with a median of ${T}_{\rm d} = 38.3^{+0.4}_{-0.9}$ K. We do not find a redshift evolution in dust temperature for sources with $L_{\rm IR}$ > $10^{12}$ $\rm L_\odot$ at $z<3$. However, we find a moderate correlation where dust temperature increases with the deviation from the SFR-$M_{\ast}$ relation. The increase in dust temperature also correlates with optical morphology, which is consistent with merger-triggered starbursts in sub-millimeter galaxies. Our galaxies do not show the tight IRX-$β_{\rm UV}$ correlation that has been observed in the local Universe. We construct the infrared luminosity functions of our 450-$μ$m sources and measure their comoving SFR densities. The contribution of the $L_{\rm IR}$ > $10^{12}$ $\rm L_\odot$ population to the SFR density rises dramatically from $z$ = 0 to 2 ($\propto$ ($1+z$)$^{3.9\pm1.1}$) and dominates the total SFR density at $z \gtrsim 2$.

preprint2020arXiv

SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES) IV: Spatial clustering and halo masses of 450-$μ$m-selected sub-millimeter galaxies

We analyze an extremely deep 450-$μ$m image ($1σ=0.56$\,mJy\,beam$^{-1}$) of a $\simeq 300$\,arcmin$^{2}$ area in the CANDELS/COSMOS field as part of the SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES). We select a robust (signal-to-noise ratio $\geqslant 4$) and flux-limited ($\geqslant 4$\,mJy) sample of 164 sub-millimeter galaxies (SMGs) at 450-$μ$m that have $K$-band counterparts in the COSMOS2015 catalog identified from radio or mid-infrared imaging. Utilizing this SMG sample and the 4705 $K$-band-selected non-SMGs that reside within the noise level $\leqslant 1$\,mJy\,beam$^{-1}$ region of the 450-$μ$m image as a training set, we develop a machine-learning classifier using $K$-band magnitude and color-color pairs based on the thirteen-band photometry available in this field. We apply the trained machine-learning classifier to the wider COSMOS field (1.6\,deg$^{2}$) using the same COSMOS2015 catalog and identify a sample of 6182 450-$μ$m SMG candidates with similar colors. The number density, radio and/or mid-infrared detection rates, redshift and stellar mass distributions, and the stacked 450-$μ$m fluxes of these SMG candidates, from the S2COSMOS observations of the wide field, agree with the measurements made in the much smaller CANDELS field, supporting the effectiveness of the classifier. Using this 450-$μ$m SMG candidate sample, we measure the two-point autocorrelation functions from $z=3$ down to $z=0.5$. We find that the 450-$μ$m SMG candidates reside in halos with masses of $\simeq (2.0\pm0.5) \times10^{13}\,h^{-1}\,\rm M_{\odot}$ across this redshift range. We do not find evidence of downsizing that has been suggested by other recent observational studies.

preprint2020arXiv

The search for life and a new logic

Exploring the Universe is one of the great unifying themes of humanity. Part of this endeavour is the search for extraterrestrial life. But how likely is it that we will find life, or that if we do it will be similar to ourselves? And therefore how do we know where and how to look? We give examples of the sort of reasoning that has been used to narrow and focus this search and we argue that obvious extensions to that logical framework will result in greater success.

preprint2020arXiv

Understanding galaxy formation and evolution through an all-sky submillimetre spectroscopic survey

We illustrate the extraordinary discovery potential for extragalactic astrophysics of a far-IR/submm all-sky spectroscopic survey with a 3m-class space telescope. Spectroscopy provides both a 3D view of the Universe and allows us to take full advantage of the sensitivity of present-day instrumentation, overcoming the spatial confusion that affects broadband far-IR/submm surveys. Emission lines powered by star formation will be detected in galaxies out to $z \simeq 8$. It will provide measurements of spectroscopic redshifts, SFRs, dust masses, and metal content for millions of galaxies at the peak epoch of cosmic star formation and of hundreds of them at the epoch of reionization. Many of these galaxies will be strongly lensed; the brightness amplification and stretching of their sizes will make it possible to investigate (by means of follow-up with high-resolution instruments) their internal structure and dynamics on the scales of giant molecular clouds. This will provide direct information on the physics driving the evolution. Furthermore, the arc-min resolution of the telescope at submm wavelengths is ideal for detecting the cores of galaxy proto-clusters, out to the epoch of reionization. Tens of millions of these galaxy-clusters-in-formation will be detected at $z \simeq 2$-3, with a tail out to $z \simeq 7$, and thousands of detections at 6 < z < 7. Their study will allow us to track the growth of the most massive halos well beyond what is possible with classical cluster surveys (mostly limited to $z < 1.5$-2), tracing the history of star formation in dense environments and teaching us how star formation and galaxy-cluster formation are related across all epochs. Such a survey will overcome the current lack of spectroscopic redshifts of dusty star-forming galaxies and galaxy proto-clusters, representing a quantum leap in far-IR/submm extragalactic astrophysics.

preprint2019arXiv

Cosmology in front of the background: studying the growth of structure at CMB wavelengths

Canada has thriving communities in CMB (cosmic microwave background) studies, cosmology and submillimetre (submm) astronomy, with involvement in many facilities that featured prominently in previous Astronomy Long Range Plans. The standard cosmological model continues to be well fit using a small number of parameters. No one expects this model to be complete and so we need to continue to challenge it with data; moreover, it does not explain how galaxies and other structures form. So, how do we improve the precision of our understanding of structure formation within this model? Wavelengths from the microwave to the submm will be particularly fruitful for answering this question. That's because, in addition to the CMB anisotropies, there are other signals that can be extracted from large maps at these wavelengths - particularly the cosmic infrared and submm backgrounds, the thermal and kinetic Sunyaev-Zeldovich effects, and CMB lensing. Such signals carry a wealth of information about the cosmological model, as well as how dust, gas and star-formation evolve within dark-matter halos. Cross-correlations between these signals and those coming from the radio, optical and X-ray surveys, will provide even more information. Canadians are already members of teams for several related facilities and are working to be involved in others. In order for Canada to be fully engaged in exploiting the detailed information coming from these cosmological signatures, it is crucial that we find the resources to participate competitively in a combination of projects currently being planned. Examples include CMB-S4, CCAT-prime, AtLAST, a new camera for JCMT, balloon projects such as BFORE and a future ambitious CMB satellite.

preprint2019arXiv

Development Plans for the Atacama Large Millimeter/submillimeter Array (ALMA)

(abridged) The Atacama Large Millimeter/submillimeter Array (ALMA) was the top-ranked priority for a new ground-based facility in the 2000 Canadian Long Range Plan. Ten years later, at the time of LRP2010, ALMA construction was well underway, with first science observations anticipated for 2011. In the past 8 years, ALMA has proved itself to be a high-impact, high-demand observatory, with record numbers of proposals submitted to the annual calls and large numbers of highly cited scientific papers across fields from protoplanetary disks to high-redshift galaxies and quasars. The LRP2010 ALMA white paper laid out 8 specific metrics that could be used to judge the success of Canada's participation in ALMA. Among these metrics were publications (number; impact), collaborations (international; multi-wavelength), and student training. To call out one particular metric, Canadians are making excellent use of ALMA in training graduate students and postdocs: as of June 2018, 12 of 23 Canadian first-author papers were led by a graduate student, and a further 4 papers were led by postdocs. All 8 metrics argue for Canada's involvement in ALMA over the past decade to be judged a success. The successful achievement of these wide-ranging goals argues strongly for Canada's continuing participation in ALMA over the next decade and beyond. Looking forward, our community needs to: (1) maintain Canadian access to ALMA and our competitiveness in using ALMA; (2) preserve full Canadian funding for our share of ALMA operations; (3) identify components of ALMA development in which Canada can play a significant role, including stimulating expertise in submillimetre instrumentation to capitalize on future opportunities; and (4) keep Canadians fully trained and engaged in ALMA, as new capabilities become available, reaching the widest possible community of potential users.

preprint2019arXiv

The Canada-France Imaging Survey: Reconstructing the Milky Way Star Formation History from its White Dwarf Population

As the remnants of stars with initial masses $\lesssim$ 8 M$_{\odot}$, white dwarfs contain valuable information on the formation histories of stellar populations. In this paper, we use deep, high-quality, u-band photometry from the Canada France Imaging Survey (CFIS), griz photometry from Pan-STARRS 1 (PS1), as well as proper motions from Gaia DR2, to select 25,156 white dwarf candidates over $\sim$4500 deg$^2$ using a reduced proper motion diagram. We develop a new white dwarf population synthesis code that returns mock observations of the Galactic field white dwarf population for a given star formation history, while simultaneously taking into account the geometry of the Milky Way, survey parameters, and selection effects. We use this model to derive the star formation histories of the thin disk, thick disk, and stellar halo. Our results show that the Milky Way disk began forming stars (11.3 $\pm$ 0.5) Gyr ago, with a peak rate of (8.8 $\pm$ 1.4) M$_{\odot}$yr$^{-1}$ at (9.8 $\pm$ 0.4) Gyr, before a slow decline to a constant rate until the present day --- consistent with recent results suggesting a merging event with a satellite galaxy. Studying the residuals between the data and best-fit model shows evidence for a slight increase in star formation over the past 3 Gyr. We fit the local fraction of helium-atmosphere white dwarfs to be (21 $\pm$ 3) %. Incorporating this methodology with data from future wide-field surveys such as LSST, Euclid, CASTOR, and WFIRST should provide an unprecedented view into the formation of the Milky Way at its earliest epoch through its white dwarfs.

preprint2016arXiv

A new VLA/e-MERLIN limit on central images in the gravitational lens system CLASS B1030+074

We present new VLA 22-GHz and e-MERLIN 5-GHz observations of CLASS B1030+074, a two-image strong gravitational lens system whose background source is a compact flat-spectrum radio quasar. In such systems we expect a third image of the background source to form close to the centre of the lensing galaxy. The existence and brightness of such images is important for investigation of the central mass distributions of lensing galaxies, but only one secure detection has been made so far in a galaxy-scale lens system. The noise levels achieved in our new B1030+074 images reach 3 microJy/beam and represent an improvement in central image constraints of nearly an order of magnitude over previous work, with correspondingly better resulting limits on the shape of the central mass profile of the lensing galaxy. Simple models with an isothermal outer power law slope now require either the influence of a central supermassive black hole, or an inner power law slope very close to isothermal, in order to suppress the central image below our detection limit. Using the central mass profiles inferred from light distributions in Virgo galaxies, moved to z=0.5, and matching to the observed Einstein radius, we now find that 45% of such mass profiles should give observable central images, 10% should give central images with a flux density still below our limit, and the remaining systems have extreme demagnification produced by the central SMBH. Further observations of similar objects will therefore allow proper statistical constraints to be placed on the central properties of elliptical galaxies at high redshift.

preprint2016arXiv

ALMA observations of Lyman-alpha Blob 1: halo sub-structure illuminated from within

We present new Atacama Large Millimeter/Submillimeter Array (ALMA) 850um continuum observations of the original Lyman-alpha Blob (LAB) in the SSA22 field at z=3.1 (SSA22-LAB01). The ALMA map resolves the previously identified submillimeter source into three components with total flux density S_850 = 1.68+/-0.06 mJy, corresponding to a star formation rate of ~150 M_sun/yr. The submillimeter sources are associated with several faint (m~27 mag) rest-frame ultraviolet sources identified in Hubble Space Telescope Imaging Spectrograph (STIS) clear filter imaging (~5850A). One of these companions is spectroscopically confirmed with Keck MOSFIRE to lie within 20 projected kpc and 250 km/s of one of the ALMA components. We postulate that some of these STIS sources represent a population of low-mass star-forming satellites surrounding the central submillimeter sources, potentially contributing to their growth and activity through accretion. Using a high resolution cosmological zoom simulation of a 10^13 M_sun halo at z=3, including stellar, dust and Ly-alpha radiative transfer, we can model the ALMA+STIS observations and demonstrate that Ly-alpha photons escaping from the central submillimeter sources are expected to resonantly scatter in neutral hydrogen, the majority of which is predicted to be associated with halo substructure. We show how this process gives rise to extended Ly-alpha emission with similar surface brightness and morphology to observed giant LABs.

preprint2016arXiv

Balloon-Borne Submillimeter Polarimetry of the Vela C Molecular Cloud: Systematic Dependence of Polarization Fraction on Column Density and Local Polarization-Angle Dispersion

We present results for Vela C obtained during the 2012 flight of the Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol). We mapped polarized intensity across almost the entire extent of this giant molecular cloud, in bands centered at 250, 350, and 500 μm. In this initial paper, we show our 500 μm data smoothed to a resolution of 2.5 arcminutes (approximately 0.5 pc). We show that the mean level of the fractional polarization p and most of its spatial variations can be accounted for using an empirical three-parameter power-law fit, p = p_0 N^(-0.4) S^(-0.6), where N is the hydrogen column density and S is the polarization-angle dispersion on 0.5 pc scales. The decrease of p with increasing S is expected because changes in the magnetic field direction within the cloud volume sampled by each measurement will lead to cancellation of polarization signals. The decrease of p with increasing N might be caused by the same effect, if magnetic field disorder increases for high column density sightlines. Alternatively, the intrinsic polarization efficiency of the dust grain population might be lower for material along higher density sightlines. We find no significant correlation between N and S. Comparison of observed submillimeter polarization maps with synthetic polarization maps derived from numerical simulations provides a promising method for testing star formation theories. Realistic simulations should allow for the possibility of variable intrinsic polarization efficiency. The measured levels of correlation among p, N, and S provide points of comparison between observations and simulations.

preprint2016arXiv

CMB-S4 Science Book, First Edition

This book lays out the scientific goals to be addressed by the next-generation ground-based cosmic microwave background experiment, CMB-S4, envisioned to consist of dedicated telescopes at the South Pole, the high Chilean Atacama plateau and possibly a northern hemisphere site, all equipped with new superconducting cameras. CMB-S4 will dramatically advance cosmological studies by crossing critical thresholds in the search for the B-mode polarization signature of primordial gravitational waves, in the determination of the number and masses of the neutrinos, in the search for evidence of new light relics, in constraining the nature of dark energy, and in testing general relativity on large scales.

preprint2016arXiv

Deep 3 GHz Observations of the Lockman Hole North with the Very Large Array - 1. Source extraction and uncertainty analysis

This is the first of two papers describing the observations and cataloguing of deep 3-GHz observations of the Lockman Hole North using the Karl G. Jansky Very Large Array. The aim of this paper is to investigate, through the use of simulated images, the uncertainties and accuracy of source-finding routines, as well as to quantify systematic effects due to resolution, such as source confusion and source size. While these effects are not new, this work is intended as a particular case study that can be scaled and translated to other surveys. We use the simulations to derive uncertainties in the fitted parameters, as well as bias corrections for the actual catalogue (presented in Paper 2). We compare two different source-finding routines, OBIT and AEGEAN, and two different effective resolutions, 8 and 2.75 arcsec. We find that the two routines perform comparably well, with OBIT being slightly better at de-blending sources, but slightly worse at fitting resolved sources. We show that 30 to 70 per cent of sources are missed or fit inaccurately once the source size becomes larger than the beam, possibly explaining source count errors in high-resolution surveys. We also investigate the effect of blending, finding that any sources with separations smaller than the beam size are fit as single sources. We show that the use of machine-learning techniques can correctly identify blended sources up to 90 per cent of the time, and prior-driven fitting can lead to a 70 per cent improvement in the number of de-blended sources.

preprint2016arXiv

Deep 3 GHz Observations of the Lockman Hole North with the Very Large Array -2. Catalogue and μJy source properties

This is the second of two papers describing the observations and source catalogues derived from sensitive 3-GHz images of the Lockman Hole North using the Karl G. Jansky Very Large Array. We describe the reduction and cataloguing process, which yielded an image with 8 arcsec resolution and instrumental noise of σ$_{\rm n}$=1.0 μJy beam$^{-1}$ rms (before primary beam corrections) and a catalogue of 558 sources detected above 5σ$_{\rm n}$. We include details of how we estimate source spectral indices across the 2-GHz VLA bandwidth, finding a median index of -0.76$\pm$0.04. Stacking of source spectra reveals a flattening of spectral index with decreasing flux density. We present a source count derived from the catalogue. We show a traditional count estimate compared with a completely independent estimate made via a P(D) confusion analysis, and find very good agreement. Cross-matches of the catalogue with X-ray, optical, infrared, radio, and redshift catalogues are also presented. The X-ray, optical and infrared data, as well as AGN selection criteria allow us to classify 10 per cent as radio-loud AGN, 28 percent as radio-quiet AGN, and 58 per cent as star-forming galaxies, with only 4 per cent unclassified.

preprint2016arXiv

HerMES: A search for high-redshift dusty galaxies in the HerMES Large Mode Survey - Catalogue, number counts and early results

Selecting sources with rising flux densities towards longer wavelengths from Herschel/SPIRE maps is an efficient way to produce a catalogue rich in high-redshift (z > 4) dusty star-forming galaxies. The effectiveness of this approach has already been confirmed by spectroscopic follow-up observations, but the previously available catalogues made this way are limited by small survey areas. Here we apply a map-based search method to 274 deg$^2$ of the HerMES Large Mode Survey (HeLMS) and create a catalogue of 477 objects with SPIRE flux densities $S_{500} > S_{350} >S_{250}$ and a $5 σ$ cut-off $S_{500} > $ 52 mJy. From this catalogue we determine that the total number of these "red" sources is at least an order of magnitude higher than predicted by galaxy evolution models. These results are in agreement with previous findings in smaller HerMES fields; however, due to our significantly larger sample size we are also able to investigate the shape of the red source counts for the first time. We have obtained spectroscopic redshift measurements for two of our sources using the Atacama Large Millimeter/submillimeter Array (ALMA). The redshifts z = 5.1 and z = 3.8 confirm that with our selection method we can indeed find high-redshift dusty star-forming galaxies.

preprint2016arXiv

How much cosmological information can be measured?

It has become common to call this the "era of precision cosmology," and hence one rarely hears about the finiteness of the amount of information that is available for constraining cosmological parameters. Under the assumption that the perturbations are purely Gaussian, the amount of extractable information (in terms of total signal-to-noise ratio for power spectrum measurements) is the same (up to a small numerical factor) as an accounting of the number of observable modes. For studies of the microwave sky, we are probably within a factor of a few of the amount of accessible information. To dramatically reduce the uncertainties on parameters will require three-dimensional probes, such as ambitious future redshifted 21-cm surveys. However, even there the available information is still finite, with the total effective signal-to-noise ratio on parameters probably not exceeding $10^7$. The amount of observable information will increase with time (but very slowly) into the extremely distant future.

preprint2016arXiv

Instrumental performance and results from testing of the BLAST-TNG receiver, submillimeter optics, and MKID arrays

Polarized thermal emission from interstellar dust grains can be used to map magnetic fields in star forming molecular clouds and the diffuse interstellar medium (ISM). The Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol) flew from Antarctica in 2010 and 2012 and produced degree-scale polarization maps of several nearby molecular clouds with arcminute resolution. The success of BLASTPol has motivated a next-generation instrument, BLAST-TNG, which will use more than 3000 linear polarization sensitive microwave kinetic inductance detectors (MKIDs) combined with a 2.5m diameter carbon fiber primary mirror to make diffraction-limited observations at 250, 350, and 500 $μ$m. With 16 times the mapping speed of BLASTPol, sub-arcminute resolution, and a longer flight time, BLAST-TNG will be able to examine nearby molecular clouds and the diffuse galactic dust polarization spectrum in unprecedented detail. The 250 $μ$m detector array has been integrated into the new cryogenic receiver, and is undergoing testing to establish the optical and polarization characteristics of the instrument. BLAST-TNG will demonstrate the effectiveness of kilo-pixel MKID arrays for applications in submillimeter astronomy. BLAST-TNG is scheduled to fly from Antarctica in December 2017 for 28 days and will be the first balloon-borne telescope to offer a quarter of the flight for "shared risk" observing by the community.

preprint2016arXiv

Multi-wavelength lens reconstruction of a $\textit{Planck}$ $\&$ $\textit{Herschel}$-detected star-bursting galaxy

We present a source-plane reconstruction of a ${\it Herschel}$ and ${\it Planck}$-detected gravitationally-lensed dusty star-forming galaxy (DSFG) at $z=1.68$ using {\it Hubble}, Sub-millimeter Array (SMA), and Keck observations. The background sub-millimeter galaxy (SMG) is strongly lensed by a foreground galaxy cluster at $z=0.997$ and appears as an arc of length $\sim 15^{\prime \prime}$ in the optical images. The continuum dust emission, as seen by SMA, is limited to a single knot within this arc. We present a lens model with source plane reconstructions at several wavelengths to show the difference in magnification between the stars and dust, and highlight the importance of a multi-wavelength lens models for studies involving lensed DSFGs. We estimate the physical properties of the galaxy by fitting the flux densities to model SEDs leading to a magnification-corrected star formation rate of $390 \pm 60$ M$_{\odot}$ yr$^{-1}$ and a stellar mass of $1.1 \pm 0.4\times 10^{11}$ M$_{\odot}$. These values are consistent with high-redshift massive galaxies that have formed most of their stars already. The estimated gas-to-baryon fraction, molecular gas surface density, and SFR surface density have values of $0.43 \pm 0.13$, $350 \pm 200$ M$_{\odot}$ pc$^{-2}$, and $\sim 12 \pm 7~$M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$, respectively. The ratio of star formation rate surface density to molecular gas surface density puts this among the most star-forming systems, similar to other measured SMGs and local ULIRGS.

preprint2016arXiv

Pi in the sky

Deviations of the observed cosmic microwave background (CMB) from the standard model, known as "anomalies", are obviously highly significant and deserve to be pursued more aggressively in order to discover the physical phenomena underlying them. Through intensive investigation we have discovered that there are equally surprising features in the digits of the number $π$, and moreover there is a remarkable correspondence between each type of peculiarity in the digits of $π$ and the anomalies in the CMB. Putting aside the unreasonable possibility that these are just the sort of flukes that appear when one looks hard enough, the only conceivable conclusion is that, however the CMB anomalies were created, a similar process imprinted patterns in the digits of $π$.

preprint2016arXiv

Star Formation in Luminous Quasars at 2<z<3

We investigate the relation between star formation rates ($\dot{M}_{s}$) and AGN properties in optically selected type 1 quasars at $2<z<3$ using data from Herschel and the SDSS. We find that $\dot{\rm{M}}_s$ remains approximately constant with redshift, at $300\pm100~\rm{M}_{\odot}$yr$^{-1}$. Conversely, $\dot{\rm{M}}_s$ increases with AGN luminosity, up to a maximum of $\sim600~\rm{M}_{\odot}$yr$^{-1}$, and with CIV FWHM. In context with previous results, this is consistent with a relation between $\dot{\rm{M}}_s$ and black hole accretion rate ($\dot{\rm{M}}_{bh}$) existing in only parts of the $z-\dot{\rm{M}}_{s}-\dot{\rm{M}}_{bh}$ plane, dependent on the free gas fraction, the trigger for activity, and the processes that may quench star formation. The relations between $\dot{\rm{M}}_s$ and both AGN luminosity and CIV FWHM are consistent with star formation rates in quasars scaling with black hole mass, though we cannot rule out a separate relation with black hole accretion rate. Star formation rates are observed to decline with increasing CIV equivalent width. This decline can be partially explained via the Baldwin effect, but may have an additional contribution from one or more of three factors; $M_i$ is not a linear tracer of L$_{2500}$, the Baldwin effect changes form at high AGN luminosities, and high CIV EW values signpost a change in the relation between $\dot{\rm{M}}_s$ and $\dot{\rm{M}}_{bh}$. Finally, there is no strong relation between $\dot{\rm{M}}_s$ and Eddington ratio, or the asymmetry of the CIV line. The former suggests that star formation rates do not scale with how efficiently the black hole is accreting, while the latter is consistent with CIV asymmetries arising from orientation effects.

preprint2016arXiv

Submillimeter Polarization Spectrum in the Vela C Molecular Cloud

Polarization maps of the Vela C molecular cloud were obtained at 250, 350, and 500um during the 2012 flight of the balloon-borne telescope BLASTPol. These measurements are used in conjunction with 850um data from Planck to study the submillimeter spectrum of the polarization fraction for this cloud. The spectrum is relatively flat and does not exhibit a pronounced minimum at λ~350um as suggested by previous measurements of other molecular clouds. The shape of the spectrum does not depend strongly on the radiative environment of the dust, as quantified by the column density or the dust temperature obtained from Herschel data. The polarization ratios observed in Vela C are consistent with a model of a porous clumpy molecular cloud being uniformly heated by the interstellar radiation field.

preprint2016arXiv

The information content of cosmic microwave background anisotropies

The cosmic microwave background (CMB) contains perturbations that are close to Gaussian and isotropic. This means that its information content, in the sense of the ability to constrain cosmological models, is closely related to the number of modes probed in CMB power spectra. Rather than making forecasts for specific experimental setups, here we take a more pedagogical approach and ask how much information we can extract from the CMB if we are only limited by sample variance. We show that, compared with temperature measurements, the addition of E-mode polarization doubles the number of modes available out to a fixed maximum multipole, provided that all of the TT, TE, and EE power spectra are measured. However, the situation in terms of constraints on particular parameters is more complicated, as we explain and illustrate graphically. We also discuss the enhancements in information that can come from adding B-mode polarization and gravitational lensing. We show how well one could ever determine the basic cosmological parameters from CMB data compared with what has been achieved with Planck, which has already probed a substantial fraction of the TT information. Lastly, we look at constraints on neutrino mass as a specific example of how lensing information improves future prospects beyond the current 6-parameter model.

preprint2016arXiv

The multiplicity of 250-$μ$m Herschel sources in the COSMOS field

We investigate the multiplicity of extragalactic sources detected by the Herschel Space Observatory in the COSMOS field. Using 3.6- and 24-$μ$m catalogues, in conjunction with 250-$μ$m data from Herschel, we seek to determine if a significant fraction of Herschel sources are composed of multiple components emitting at 250 $μ$m. We use the XID+ code, using Bayesian inference methods to produce probability distributions of the possible contributions to the observed 250-$μ$m flux for each potential component. The fraction of Herschel flux assigned to the brightest component is highest for sources with total 250-$μ$m fluxes < 45 mJy; however, the flux in the brightest component is still highest in the brightest Herschel sources. The faintest 250-$μ$m sources (30-45 mJy) have the majority of their flux assigned to a single bright component; the second brightest component is typically significantly weaker, and contains the remainder of the 250-$μ$m source flux. At the highest 250-$μ$m fluxes (45-110 mJy), the brightest and second brightest components are assigned roughly equal fluxes, and together are insufficient to reach 100 per cent of the 250-$μ$m source flux. This indicates that additional components are required, beyond the brightest two components, to reproduce the observed flux. 95 per cent of the sources in our sample have a second component that contains more than 10 per cent of the total source flux. Particularly for the brightest Herschel sources, assigning the total flux to a single source may overestimate the flux contributed by around 150 per cent.

preprint2016arXiv

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

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

preprint2016arXiv

The SCUBA-2 Cosmology Legacy Survey: The clustering of submillimetre galaxies in the UKIDSS UDS field

Submillimetre galaxies (SMGs) are among the most luminous dusty galaxies in the Universe, but their true nature remains unclear; are SMGs the progenitors of the massive elliptical galaxies we see in the local Universe, or are they just a short-lived phase among more typical star-forming galaxies? To explore this problem further, we investigate the clustering of SMGs identified in the SCUBA-2 Cosmology Legacy Survey. We use a catalogue of submillimetre ($850μ$m) source identifications derived using a combination of radio counterparts and colour/IR selection to analyse a sample of 914 SMGs in the UKIDSS Ultra Deep Survey (UDS), making this the largest high redshift sample of these galaxies to date. Using angular cross-correlation techniques, we estimate the halo masses for this large sample of SMGs and compare them with passive and star-forming galaxies selected in the same field. We find that SMGs, on average, occupy high-mass dark matter halos (M$_{\text{halo}} >10^{13}$M$_{\odot}$) at redshifts $z > 2.5$, consistent with being the progenitors of massive quiescent galaxies in present-day galaxy clusters. We also find evidence of downsizing, in which SMG activity shifts to lower mass halos at lower redshifts. In terms of their clustering and halo masses, SMGs appear to be consistent with other star-forming galaxies at a given redshift.

preprint2016arXiv

The SCUBA-2 Cosmology Legacy Survey: The EGS deep field I - Deep number counts and the redshift distribution of the recovered Cosmic Infrared Background at 450 and 850 um

We present deep observations at 450 um and 850 um in the Extended Groth Strip field taken with the SCUBA-2 camera mounted on the James Clerk Maxwell Telescope as part of the deep SCUBA-2 Cosmology Legacy Survey (S2CLS), achieving a central instrumental depth of $σ_{450}=1.2$ mJy/beam and $σ_{850}=0.2$ mJy/beam. We detect 57 sources at 450 um and 90 at 850 um with S/N > 3.5 over ~70 sq. arcmin. From these detections we derive the number counts at flux densities $S_{450}>4.0$ mJy and $S_{850}>0.9$ mJy, which represent the deepest number counts at these wavelengths derived using directly extracted sources from only blank-field observations with a single-dish telescope. Our measurements smoothly connect the gap between previous shallower blank-field single-dish observations and deep interferometric ALMA results. We estimate the contribution of our SCUBA-2 detected galaxies to the cosmic infrared background (CIB), as well as the contribution of 24 um-selected galaxies through a stacking technique, which add a total of $0.26\pm0.03$ and $0.07\pm0.01$ MJy/sr, at 450 um and 850 um, respectively. These surface brightnesses correspond to $60\pm20$ and $50\pm20$ per cent of the total CIB measurements, where the errors are dominated by those of the total CIB. Using the photometric redshifts of the 24 um-selected sample and the redshift distributions of the submillimetre galaxies, we find that the redshift distribution of the recovered CIB is different at each wavelength, with a peak at $z\sim1$ for 450 um and at $z\sim2$ for 850um, consistent with previous observations and theoretical models.

preprint2015arXiv

A Farewell to Falsifiability

Some of the most obviously correct physical theories - namely string theory and the multiverse - make no testable predictions, leading many to question whether we should accept something as scientific even if it makes no testable predictions and hence is not refutable. However, some far-thinking physicists have proposed instead that we should give up on the notion of Falsifiability itself. We endorse this suggestion but think it does not go nearly far enough. We believe that we should also dispense with other outdated ideas, such as Fidelity, Frugality, Factuality and other "F" words. And we quote a lot of famous people to support this view.

preprint2015arXiv

BFORE: The B-mode Foreground Experiment

The B-mode Foreground Experiment (BFORE) is a proposed NASA balloon project designed to make optimal use of the sub-orbital platform by concentrating on three dust foreground bands (270, 350, and 600 GHz) that complement ground-based cosmic microwave background (CMB) programs. BFORE will survey ~1/4 of the sky with 1.7 - 3.7 arcminute resolution, enabling precise characterization of the Galactic dust that now limits constraints on inflation from CMB B-mode polarization measurements. In addition, BFORE's combination of frequency coverage, large survey area, and angular resolution enables science far beyond the critical goal of measuring foregrounds. BFORE will constrain the velocities of thousands of galaxy clusters, provide a new window on the cosmic infrared background, and probe magnetic fields in the interstellar medium. We review the BFORE science case, timeline, and instrument design, which is based on a compact off-axis telescope coupled to >10,000 superconducting detectors.

preprint2015arXiv

Far-infrared observations of an unbiased sample of gamma-ray burst host galaxies

Gamma-ray bursts (GRBs) are the most energetic phenomena in the Universe; believed to result from the collapse and subsequent explosion of massive stars. Even though it has profound consequences for our understanding of their nature and selection biases, little is known about the dust properties of the galaxies hosting GRBs. We present analysis of the far-infrared properties of an unbiased sample of 20 \textit{BeppoSAX} and \textit{Swift} GRB host galaxies (at an average redshift of $z\,=\,3.1$) located in the {\it Herschel} Astrophysical Terahertz Large Area Survey, the {\it Herschel} Virgo Cluster Survey, the {\it Herschel} Fornax Cluster Survey, the {\it Herschel} Stripe 82 Survey and the {\it Herschel} Multi-tiered Extragalactic Survey, totalling $880$ deg$^2$, or $\sim 3$\% of the sky in total. Our sample selection is serendipitous, based only on whether the X-ray position of a GRB lies within a large-scale {\it Herschel} survey -- therefore our sample can be considered completely unbiased. Using deep data at wavelengths of 100\,--\,500$\,μ$m, we tentatively detected 1 out of 20 GRB hosts located in these fields. We constrain their dust masses and star formation rates (SFRs), and discuss these in the context of recent measurements of submillimetre galaxies and ultraluminous infrared galaxies. The average far-infrared flux of our sample gives an upper limit on SFR of $<114\,{\rm M}\odot\,\mbox{yr}^{-1}$. The detection rate of GRB hosts is consistent with that predicted assuming that GRBs trace the cosmic SFR density in an unbiased way, i.e. that the fraction of GRB hosts with $\mbox{SFR}>500\,{\rm M}\odot\,\mbox{yr}^{-1}$ is consistent with the contribution of such luminous galaxies to the cosmic star formation density.

preprint2015arXiv

GOODS-HERSCHEL: star formation, dust attenuation and the FIR-radio correlation on the Main Sequence of star-forming galaxies up to z~4

We use deep panchromatic datasets in the GOODS-N field, from GALEX to the deepest Herschel far-infrared and VLA radio continuum imaging, to explore, using mass-complete samples, the evolution of the star formation activity and dust attenuation of star-forming galaxies to z~4. Our main results can be summarized as follows: i) the slope of the SFR-M correlation is consistent with being constant, and equal to ~0.8 at least up to z~1.5, while its normalization keeps increasing with redshift; ii) for the first time here we are able to explore the FIR-radio correlation for a mass-selected sample of star-forming galaxies: the correlation does not evolve up to z~4; iii) we confirm that galaxy stellar mass is a robust proxy for UV dust attenuation in star-forming galaxies, with more massive galaxies being more dust attenuated, strikingly we find that this attenuation relation evolves very weakly with redshift, the amount of dust attenuation increasing by less than 0.3 magnitudes over the redshift range [0.5-4] for a fixed stellar mass, as opposed to a tenfold increase of star formation rate; iv) the correlation between dust attenuation and the UV spectral slope evolves in redshift, with the median UV spectral slope of star-forming galaxies becoming bluer with redshift. By z~3, typical UV slopes are inconsistent, given the measured dust attenuation, with the predictions of commonly used empirical laws. Finally, building on existing results, we show that gas reddening is marginally larger (by a factor of around 1.3) than stellar reddening at all redshifts probed, and also that the amount of dust attenuation at a fixed ISM metallicity increases with redshift. We speculate that our results support evolving ISM conditions of typical star-forming galaxies such that at z~1.5 Main Sequence galaxies have ISM conditions getting closer to those of local starbursts.

preprint2015arXiv

H-ATLAS/GAMA: Quantifying the Morphological Evolution of the Galaxy Population Using Cosmic Calorimetry

Using results from the Herschel Astrophysical Terrahertz Large-Area Survey and the Galaxy and Mass Assembly project, we show that, for galaxy masses above approximately 1.0e8 solar masses, 51% of the stellar mass-density in the local Universe is in early-type galaxies (ETGs: Sersic n > 2.5) while 89% of the rate of production of stellar mass-density is occurring in late-type galaxies (LTGs: Sersic n < 2.5). From this zero-redshift benchmark, we have used a calorimetric technique to quantify the importance of the morphological transformation of galaxies over the history of the Universe. The extragalactic background radiation contains all the energy generated by nuclear fusion in stars since the Big Bang. By resolving this background radiation into individual galaxies using the deepest far-infrared survey with the Herschel Space Observatory and a deep near-infrared/optical survey with the Hubble Space Telescope (HST), and using measurements of the Sersic index of these galaxies derived from the HST images, we estimate that approximately 83% of the stellar mass-density formed over the history of the Universe occurred in LTGs. The difference between this and the fraction of the stellar mass-density that is in LTGs today implies there must have been a major transformation of LTGs into ETGs after the formation of most of the stars.

preprint2015arXiv

HerMES: ALMA Imaging of Herschel-selected Dusty Star-forming Galaxies

The Herschel Multi-tiered Extragalactic Survey (HerMES) has identified large numbers of dusty star-forming galaxies (DSFGs) over a wide range in redshift. A detailed understanding of these DSFGs is hampered by the limited spatial resolution of Herschel. We present 870um 0.45" resolution imaging from the Atacama Large Millimeter/submillimeter Array (ALMA) of 29 HerMES DSFGs with far-infrared (FIR) flux densities in between the brightest of sources found by Herschel and fainter DSFGs found in ground-based sub-millimeter (sub-mm) surveys. We identify 62 sources down to the 5-sigma point-source sensitivity limit in our ALMA sample (sigma~0.2mJy), of which 6 are strongly lensed (showing multiple images) and 36 experience significant amplification (mu>1.1). To characterize the properties of the ALMA sources, we introduce and make use of uvmcmcfit, a publicly available Markov chain Monte Carlo analysis tool for interferometric observations of lensed galaxies. Our lens models tentatively favor intrinsic number counts for DSFGs with a steep fall off above 8mJy at 880um. Nearly 70% of the Herschel sources comprise multiple ALMA counterparts, consistent with previous research indicating that the multiplicity rate is high in bright sub-mm sources. Our ALMA sources are located significantly closer to each other than expected based on results from theoretical models as well as fainter DSFGs identified in the LABOCA ECDFS Submillimeter Survey. The high multiplicity rate and low projected separations argue in favor of interactions and mergers driving the prodigious emission from the brightest DSFGs as well as the sharp downturn above S_880=8mJy.

preprint2015arXiv

HerMES: Current Cosmic Infrared Background Estimates Can be Explained by Known Galaxies and their Faint Companions at z < 4

We report contributions to cosmic infrared background (CIB) intensities originating from known galaxies and their faint companions at submillimeter wavelengths. Using the publicly-available UltraVISTA catalog, and maps at 250, 350, and 500 μm from the \emph{Herschel} Multi-tiered Extragalactic Survey (HerMES), we perform a novel measurement that exploits the fact that uncatalogued sources may bias stacked flux densities --- particularly if the resolution of the image is poor --- and intentionally smooth the images before stacking and summing intensities. By smoothing the maps we are capturing the contribution of faint (undetected in K_S ~ 23.4) sources that are physically associated, or correlated, with the detected sources. We find that the cumulative CIB increases with increased smoothing, reaching 9.82 +- 0.78, 5.77 +- 0.43, and 2.32 +- 0.19$\, \rm nW m^{-2} sr^{-1}$ at 250, 350, and 500 μm at 300 arcsec FWHM. This corresponds to a fraction of the fiducial CIB of 0.94 +- 0.23, 1.07 +- 0.31, and 0.97 +- 0.26 at 250, 350, and 500 μm, where the uncertainties are dominated by those of the absolute CIB. We then propose, with a simple model combining parametric descriptions for stacked flux densities and stellar mass functions, that emission from galaxies with log(M/Msun) > 8.5 can account for the most of the measured total intensities, and argue against contributions from extended, diffuse emission. Finally, we discuss prospects for future survey instruments to improve the estimates of the absolute CIB levels, and observe any potentially remaining emission at z > 4.

preprint2015arXiv

HerMES: Disentangling active galactic nuclei and star formation in the radio source population

We separate the extragalactic radio source population above ~50 uJy into active galactic nuclei (AGN) and star-forming sources. The primary method of our approach is to fit the infrared spectral energy distributions (SEDs), constructed using Spitzer/IRAC and MIPS and Herschel/SPIRE photometry, of 380 radio sources in the Extended Chandra Deep Field-South. From the fitted SEDs, we determine the relative AGN and star-forming contributions to their infrared emission. With the inclusion of other AGN diagnostics such as X-ray luminosity, Spitzer/IRAC colours, radio spectral index and the ratio of star-forming total infrared flux to k-corrected 1.4 GHz flux density, qIR, we determine whether the radio emission in these sources is powered by star formation or by an AGN. The majority of these radio sources (60 per cent) show the signature of an AGN at some wavelength. Of the sources with AGN signatures, 58 per cent are hybrid systems for which the radio emission is being powered by star formation. This implies that radio sources which have likely been selected on their star formation have a high AGN fraction. Below a 1.4 GHz flux density of 1 mJy, along with finding a strong contribution to the source counts from pure star-forming sources, we find that hybrid sources constitute 20-65 per cent of the sources. This result suggests that hybrid sources have a significant contribution, along with sources that do not host a detectable AGN, to the observed flattening of the source counts at ~1mJy for the extragalactic radio source population.

preprint2015arXiv

Probing the diffuse baryon distribution with the lensing-tSZ cross-correlation

Approximately half of the Universe's baryons are in a form that has been hard to detect directly. However, the missing component can be traced through the cross-correlation of the thermal Sunyaev-Zeldovich (tSZ) effect with weak gravitational lensing. We build a model for this correlation and use it to constrain the extended baryon component, employing data from the Canada France Hawaii Lensing Survey and the {\it Planck\/} satellite. The measured correlation function is consistent with an isothermal $β$-model for the halo gas pressure profile, and the 1- and 2-halo terms are both detected at the 4$σ$ level. In addition, we measure the hydrostatic mass bias $(1-b)=0.79^{+0.07}_{-0.10}$, which is consistent with numerical simulation results and the constraints from X-ray observations. The effective temperature of the gas is found to be in the range ($7\times10^{5}$--$3 \times10^{8}$)\,K, with approximately $50\%$ of the baryons appearing to lie beyond the virial radius of the halos, consistent with current expectations for the warm-hot intergalactic medium.

preprint2015arXiv

SEDeblend: A new method for deblending spectral energy distributions in confused imaging

For high-redshift submillimetre or millimetre sources detected with single dish telescopes, interferometric follow-up has shown that many are multiple submm galaxies blended together. Confusion-limited Herschel observations of such targets are also available, and these sample the peak of their spectral energy distribution in the far-infrared. Many methods for analysing these data have been adopted, but most follow the traditional approach of extracting fluxes before model spectral energy distributions are fit, which has the potential to erase important information on degeneracies among fitting parameters and glosses over the intricacies of confusion noise. Here, we adapt the forward-modelling method that we originally developed to disentangle a high-redshift strongly-lensed galaxy group, in order to tackle this problem in a more statistically rigorous way, by combining source deblending and SED fitting into the same procedure. We call this method "SEDeblend." As an application, we derive constraints on far-infrared luminosities and dust temperatures for sources within the ALMA follow-up of the LABOCA Extended Chandra Deep Field South Submillimetre Survey. We find an average dust temperature for an 870 micron-selected sample of (33.9+-2.4) K for the full survey. When selection effects of the sample are considered, we find no evidence that the average dust temperature evolves with redshift.

preprint2015arXiv

The HerMES sub-millimetre local and low-redshift luminosity functions

We used wide area surveys over 39 deg$^2$ by the HerMES collaboration, performed with the Herschel Observatory SPIRE multi-wavelength camera, to estimate the low-redshift, $0.02<z<0.5$, monochromatic luminosity functions (LFs) of galaxies at 250, 350 and 500$\,μ$m. SPIRE flux densities were also combined with Spitzer photometry and multi-wavelength archival data to perform a complete SED fitting analysis of SPIRE detected sources to calculate precise k-corrections, as well as the bolometric infrared (8-1000$\,μ$m) luminosity functions and their low-$z$ evolution from a combination of statistical estimators. Integration of the latter prompted us to also compute the local luminosity density (LLD) and the comoving star formation rate density (SFRD) for our sources, and to compare them with theoretical predictions of galaxy formation models. The luminosity functions show significant and rapid luminosity evolution already at low redshifts, $0.02<z<0.2$, with L$_{IR}^* \propto (1+z)^{6.0\pm0.4}$ and $Φ_{IR}^* \propto (1+z)^{-2.1\pm0.4}$, L$_{250}^* \propto (1+z)^{5.3\pm0.2}$ and $Φ_{250}^* \propto (1+z)^{-0.6\pm0.4}$ estimated using the IR bolometric and the 250$\,μ$m LFs respectively. Converting our IR LD estimate into an SFRD assuming a standard Salpeter IMF and including the unobscured contribution based on the UV dust-uncorrected emission from local galaxies, we estimate a SFRD scaling of SFRD$_0+0.08 z$, where SFRD$_0\simeq (1.9\pm 0.03)\times 10^{-2} [\mathrm{M}_\odot\,\mathrm{Mpc}^{-3}]$ is our total SFRD estimate at $z\sim0.02$.

preprint2015arXiv

The Herschel view of the dominant mode of galaxy growth from z=4 to the present day

We present an analysis of the deepest Herschel images in four major extragalactic fields GOODS-North, GOODS-South, UDS and COSMOS obtained within the GOODS-Herschel and CANDELS-Herschel key programs. The picture provided by 10497 individual far-infrared detections is supplemented by the stacking analysis of a mass-complete sample of 62361 star-forming galaxies from the CANDELS-HST H band-selected catalogs and from two deep ground-based Ks band-selected catalogs in the GOODS-North and the COSMOS-wide fields, in order to obtain one of the most accurate and unbiased understanding to date of the stellar mass growth over the cosmic history. We show, for the first time, that stacking also provides a powerful tool to determine the dispersion of a physical correlation and describe our method called "scatter stacking" that may be easily generalized to other experiments. We demonstrate that galaxies of all masses from z=4 to 0 follow a universal scaling law, the so-called main sequence of star-forming galaxies. We find a universal close-to-linear slope of the logSFR-logM* relation with evidence for a flattening of the main sequence at high masses (log(M*/Msun) > 10.5) that becomes less prominent with increasing redshift and almost vanishes by z~2. This flattening may be due to the parallel stellar growth of quiescent bulges in star-forming galaxies. Within the main sequence, we measure a non varying SFR dispersion of 0.3 dex. The specific SFR (sSFR=SFR/M*) of star-forming galaxies is found to continuously increase from z=0 to 4. Finally we discuss the implications of our findings on the cosmic SFR history and show that more than 2/3 of present-day stars must have formed in a regime dominated by the main sequence mode. As a consequence we conclude that, although omnipresent in the distant Universe, galaxy mergers had little impact in shaping the global star formation history over the last 12.5 Gyr.

preprint2015arXiv

The SCUBA-2 Cosmology Legacy Survey: ALMA resolves the bright-end of the sub-millimeter number counts

We present high-resolution 870-um ALMA continuum maps of 30 bright sub-millimeter sources in the UKIDSS UDS field. These sources are selected from deep, 1-square degrees 850-um maps from the SCUBA--2 Cosmology Legacy Survey, and are representative of the brightest sources in the field (median SCUBA2 flux S_850=8.7+/-0.4 mJy). We detect 52 sub-millimeter galaxies (SMGs) at >4-sigma significance in our 30 ALMA maps. In 61+/-17% of the ALMA maps the single-dish source comprises a blend of >=2 SMGs, where the secondary SMGs are Ultra--Luminous Infrared Galaxies (ULIRGs) with L_IR>10^12 Lo. The brightest SMG contributes on average 80+/-4% of the single-dish flux density, and in the ALMA maps containing >=2 SMGs the secondary SMG contributes 25+/-3% of the integrated ALMA flux. We construct source counts and show that multiplicity boosts the apparent single-dish cumulative counts by 20% at S_870>7.5mJy, and by 60% at S_870>12mJy. We combine our sample with previous ALMA studies of fainter SMGs and show that the counts are well-described by a double power-law with a break at 8.5+/-0.6mJy. The break corresponds to a luminosity of ~6x10^12Lsol or a star-formation rate of ~1000Mo/yr. For the typical sizes of these SMGs, which are resolved in our ALMA data with r=1.2+/-0.1kpc, this yields a limiting SFR density of ~100Msol/yr/kpc2. Finally, the number density of S_870>2mJy SMGs is 80+/-30 times higher than that derived from blank-field counts. An over-abundance of faint SMGs is inconsistent with line-of-sight projections dominating multiplicity in the brightest SMGs, and indicates that a significant proportion of these high-redshift ULIRGs must be physically associated.

preprint2014arXiv

Black hole accretion preferentially occurs in gas rich galaxies

We have investigated the gas content of a sample of several hundred AGN host galaxies at z$<$1 and compared it with a sample of inactive galaxies, matched in bins of stellar mass and redshift. Gas masses have been inferred from the dust masses, obtained by stacked Herschel far-IR and sub-mm data in the GOODS and COSMOS fields, under reasonable assumptions and metallicity scaling relations for the dust-to-gas ratio. We find that AGNs are on average hosted in galaxies much more gas rich than inactive galaxies. In the vast majority of stellar mass bins, the average gas content of AGN hosts is higher than in inactive galaxies. The difference is up to a factor of ten higher in low stellar mass galaxies, with a significance of 6.5$σ$. In almost half of the AGN sample the gas content is three times higher than in the control sample of inactive galaxies. Our result strongly suggests that the probability of having an AGN activated is simply driven by the amount of gas in the host galaxy; this can be explained in simple terms of statistical probability of having a gas cloud falling into the gravitational potential of the black hole. The increased probability of an AGN being hosted by a star-forming galaxy, identified by previous works, may be a consequence of the relationship between gas content and AGN activity, found in this paper, combined with the Schmidt-Kennicutt law for star formation.

preprint2014arXiv

Deep 3 GHz Number Counts from a P(D) Fluctuation Analysis

Radio source counts constrain galaxy populations and evolution, as well as the global star formation history. However, there is considerable disagreement among the published 1.4-GHz source counts below 100 microJy. Here we present a statistical method for estimating the microJy and even sub-microJy source count using new deep wide-band 3-GHz data in the Lockman Hole from the Karl G. Jansky Very Large Array (VLA). We analyzed the confusion amplitude distribution P(D), which provides a fresh approach in the form of a more robust model, with a comprehensive error analysis. We tested this method on a large-scale simulation, incorporating clustering and finite source sizes. We discuss in detail our statistical methods for fitting using Monte Carlo Markov chains, handling correlations, and systematic errors from the use of wide-band radio interferometric data. We demonstrated that the source count can be constrained down to 50 nJy, a factor of 20 below the rms confusion. We found the differential source count near 10 microJy to have a slope of -1.7, decreasing to about -1.4 at fainter flux densities. At 3GHz the rms confusion in an 8arcsec FWHM beam is ~ 1.2 microJy/beam, and a radio background temperature ~ 14mK. Our counts are broadly consistent with published evolutionary models. With these results we were also able to constrain the peak of the Euclidean normalized differential source count of any possible new radio populations that would contribute to the cosmic radio background down to 50 nJy.

preprint2014arXiv

Detailed modelling of a large sample of Herschel sources in the Lockman Hole: identification of cold dust and of lensing candidates through their anomalous SEDs

We have studied in detail a sample of 967 SPIRE sources with 5-sigma detections at 350 and 500 micron and associations with Spitzer-SWIRE 24 micron galaxies in the HerMES-Lockman survey area, fitting their mid- and far-infrared, and submillimetre, SEDs in an automatic search with a set of six infrared templates. For almost 300 galaxies we have modelled their SEDs individually to ensure the physicality of the fits. We confirm the need for the new cool and cold cirrus templates, and also of the young starburst template, introduced in earlier work. We also identify 109 lensing candidates via their anomalous SEDs and provide a set of colour-redshift constraints which allow lensing candidates to be identified from combined Herschel and Spitzer data. The picture that emerges of the submillimetre galaxy population is complex, comprising ultraluminous and hyperluminous starbursts, lower luminosity galaxies dominated by interstellar dust emission, lensed galaxies and galaxies with surprisingly cold (10-13K) dust. 11 % of 500$μ$m selected sources are lensing candidates. 70 % of the unlensed sources are ultraluminous infrared galaxies and 26 % are hyperluminous. 34 % are dominated by optically thin interstellar dust ('cirrus') emission, but most of these are due to cooler dust than is characteristic of our Galaxy. At the highest infrared luminosities we see SEDs dominated by M82, Arp220 and young starburst types, in roughly equal proportions.

preprint2014arXiv

How does pressure gravitate? Cosmological constant problem confronts observational cosmology

An important and long-standing puzzle in the history of modern physics is the gross inconsistency between theoretical expectations and cosmological observations of the vacuum energy density, by at least 60 orders of magnitude, otherwise known as the \textit{cosmological constant problem}. A characteristic feature of vacuum energy is that it has a pressure with the same amplitude, but opposite sign to its energy density, while all the precision tests of General Relativity are either in vacuum, or for media with negligible pressure. Therefore, one may wonder whether an anomalous coupling to pressure might be responsible for decoupling vacuum from gravity. We test this possibility in the context of the \textit{Gravitational Aether} proposal, using current cosmological observations, which probe the gravity of relativistic pressure in the radiation era. Interestingly, we find that the best fit for anomalous pressure coupling is about half-way between General Relativity (GR), and Gravitational Aether (GA), if we include \textit{Planck} together with \textit{WMAP} and BICEP2 polarization cosmic microwave background (CMB) observations. Taken at face value, this data combination excludes both GR and GA at around the 3-sigma level. However, including higher resolution CMB observations ("highL") or baryonic acoustic oscillations (BAO) pushes the best fit closer to GR, excluding the Gravitational Aether solution to the cosmological constant problem at the 4--5-sigma level. This constraint effectively places a limit on the anomalous coupling to pressure in the parametrized post-Newtonian (PPN) expansion, $ζ_4 = 0.105 \pm 0.049$ (+highL CMB), or $ζ_4 = 0.066 \pm 0.039$ (+BAO). These represent the most precise measurement of this parameter to date, indicating a mild tension with GR (for $Λ$CDM including tensors, with $ζ_4=0$), and also among different data sets.

preprint2014arXiv

Imaging the environment of a z = 6.3 submillimeter galaxy with SCUBA-2

We describe a search for submillimeter emission in the vicinity of one of the most distant, luminous galaxies known, HerMES FLS3 at z=6.34, exploiting it as a signpost to a potentially biased region of the early Universe, as might be expected in hierarchical structure formation models. Imaging to the confusion limit with the innovative, wide-field submillimeter bolometer camera, SCUBA-2, we are sensitive to colder and/or less luminous galaxies in the surroundings of HFLS3. We use the Millennium Simulation to illustrate that HFLS3 may be expected to have companions if it is as massive as claimed, but find no significant evidence from the surface density of SCUBA-2 galaxies in its vicinity, or their colors, that HFLS3 marks an over-density of dusty, star-forming galaxies. We cannot rule out the presence of dusty neighbours with confidence, but deeper 450-um imaging has the potential to more tightly constrain the redshifts of nearby galaxies, at least one of which likely lies at z>~5. If associations with HFLS3 can be ruled out, this could be taken as evidence that HFLS3 is less biased than a simple extrapolation of the Millennium Simulation may imply. This could suggest either that it represents a rare short-lived, but highly luminous, phase in the evolution of an otherwise typical galaxy, or that this system has suffered amplification due to a foreground gravitational lens and so is not as intrinsically luminous as claimed.

preprint2014arXiv

Map-making for large-format detector arrays on CCAT

CCAT is a large submillimetre telescope to be built near the ALMA site in northern Chile. A large-format KID camera, with up to 48,000 detectors at a single waveband sampled at about 1 kHz, will have a data rate about 50 times larger than SCUBA-2, the largest existing submillimetre camera. Creating a map from this volume of data will be a challenge, both in terms of memory and processing time required. We investigate how to extend SMURF, the iterative map-maker used for reducing SCUBA-2 observations, to a distributed-node parallel system, and estimate how the processing time scales with the number of nodes in the system.

preprint2014arXiv

Map-making for the Next Generation of Ground-based Submillimeter Instruments

Current ground-based submillimeter instruments (e.g. SCUBA-2, SHARC-2 and LABOCA) have hundreds to thousands of detectors, sampled at tens to hundreds of hertz, generating up to hundreds of gigabytes per night. Since noise is correlated between detectors and in time, due to atmospheric signals and temperature oscillations, naive map-making is not applicable. In addition, the size of the data sets makes direct likelihood based inversion techniques intractable. As a result, the data reduction approach for most current submm cameras is to adopt iterative methods in order to separate noise from sky signal, and hence effectively produce astronomical images. We investigate how today's map-makers scale to the next generation of instruments, which will have tens of thousands of detectors sampled at thousands of hertz, leading to data sets of challenging size. We propose strategies for reducing such large data sets.

preprint2014arXiv

The Deep Diffuse Extragalactic Radio Sky at 1.75 GHz

We present a study of diffuse extragalactic radio emission at $1.75\,$GHz from part of the ELAIS-S1 field using the Australia Telescope Compact Array. The resulting mosaic is $2.46\,$deg$^2$, with a roughly constant noise region of $0.61\,$deg$^2$ used for analysis. The image has a beam size of $150 \times60\,$arcsec and instrumental $\langleσ_{\rm n}\rangle= (52\pm5)\, μ$Jy beam$^{-1}$. Using point-source models from the ATLAS survey, we subtract the discrete emission in this field for $S \ge 150\, μ$Jy beam$^{-1}$. Comparison of the source-subtracted probability distribution, or \pd, with the predicted distribution from unsubtracted discrete emission and noise, yields an excess of $(76 \pm 23) \, μ$Jy beam$^{-1}$. Taking this as an upper limit on any extended emission we constrain several models of extended source counts, assuming $Ω_{\rm source} \le 2\,$arcmin. The best-fitting models yield temperatures of the radio background from extended emission of $T_{\rm b}=(10\pm7) \,$mK, giving an upper limit on the total temperature at $1.75\,$GHz of $(73\pm10)\,$mK. Further modelling shows that our data are inconsistent with the reported excess temperature of ARCADE2 to a source-count limit of $1\, μ$Jy. Our new data close a loop-hole in the previous constraints, because of the possibility of extended emission being resolved out at higher resolution. Additionally, we look at a model of cluster halo emission and two WIMP dark matter annihilation source-count models, and discuss general constraints on any predicted counts from such sources. Finally, we report the derived integral count at $1.4\,$GHz using the deepest discrete count plus our new extended-emission limits, providing numbers that can be used for planning future ultra-deep surveys.

preprint2014arXiv

The Next Generation BLAST Experiment

The Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol) was a suborbital experiment designed to map magnetic fields in order to study their role in star formation processes. BLASTPol made detailed polarization maps of a number of molecular clouds during its successful flights from Antarctica in 2010 and 2012. We present the next-generation BLASTPol instrument (BLAST-TNG) that will build off the success of the previous experiment and continue its role as a unique instrument and a test bed for new technologies. With a 16-fold increase in mapping speed, BLAST-TNG will make larger and deeper maps. Major improvements include a 2.5 m carbon fiber mirror that is 40% wider than the BLASTPol mirror and ~3000 polarization sensitive detectors. BLAST-TNG will observe in three bands at 250, 350, and 500 microns. The telescope will serve as a pathfinder project for microwave kinetic inductance detector (MKID) technology, as applied to feedhorn coupled submillimeter detector arrays. The liquid helium cooled cryostat will have a 28-day hold time and will utilize a closed-cycle $^3$He refrigerator to cool the detector arrays to 270 mK. This will enable a detailed mapping of more targets with higher polarization resolution than any other submillimeter experiment to date. BLAST-TNG will also be the first balloon-borne telescope to offer shared risk observing time to the community. This paper outlines the motivation for the project and the instrumental design.

preprint2014arXiv

The next-generation BLASTPol experiment

The Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol) is a suborbital mapping experiment designed to study the role magnetic fields play in star formation. BLASTPol has had two science flights from McMurdo Station, Antarctica in 2010 and 2012. These flights have produced thousands of polarization vectors at 250, 350 and 500 microns in several molecular cloud targets. We present the design, specifications, and progress towards the next-generation BLASTPol experiment (BLAST-TNG). BLAST-TNG will fly a 40% larger diameter primary mirror, with almost 8 times the number of polarization-sensitive detectors resulting in a factor of 16 increase in mapping speed. With a spatial resolution of 22 arcseconds and four times the field of view of BLASTPol, BLAST-TNG will bridge the angular scales between Planck's low resolution all-sky maps and ALMA's ultra-high resolution narrow fields. The new receiver has a larger cryogenics volume, allowing for a 28 day hold time. BLAST-TNG employs three arrays of Microwave Kinetic Inductance Detectors (MKIDs) with 30% fractional bandwidth at 250, 350 and 500 microns. In this paper, we will present the new BLAST-TNG instrument and science objectives.

preprint2014arXiv

Velocities hasten to tell us about the Universe

The peculiar velocities of galaxies are driven by gravity, and hence hold the promise of probing details of how gravity forms structures. In particular it is possible to constrain cosmological parameters and to test extensions to the standard model, such as modifications to the theory of gravity or the existence of primordial density perturbations which are non-Gaussian. This constraining power has been frustrated by systematic effects, but we appear to be entering an era when velocity measurements may finally be living up to their promise.

preprint2014arXiv

WMAP Observations of Planck ESZ Clusters

We examine the Sunyaev--Zeldovich (SZ) effect in the 7-yr Wilkinson Microwave Anisotropy Probe (WMAP) data by cross-correlating it with the Planck Early-release Sunyaev--Zeldovich catalog. Our analysis proceeds in two parts. We first perform a stacking analysis in which the filtered WMAP data are averaged at the locations of the 175 Planck clusters. We then perform a regression analysis to compare the mean amplitude of the SZ signal, Y_{500}, in the WMAP data to the corresponding amplitude in the Planck data. The aggregate Planck clusters are detected in the 7-yr WMAP data with a signal-to-noise ratio of 16.3. In the regression analysis we find that the SZ amplitude measurements agree to better than 25%: $a=1.23 \pm 0.18$ for the fit Y^{wmap}_{500} = aY^{planck}_{500}$.

preprint2013arXiv

A Roadmap for Canadian Submillimetre Astronomy

We survey the present landscape in submillimetre astronomy for Canada and describe a plan for continued engagement in observational facilities to ~2020. Building on Canada's decadal Long Range Plan process, we emphasize that continued involvement in a large, single-dish facility is crucial given Canada's substantial investment in ALMA and numerous PI-led submillimetre experiments. In particular, we recommend: i) an extension of Canadian participation in the JCMT until at least the unique JCMT Legacy Survey program is able to realize the full scientific potential provided by the world-leading SCUBA-2 instrument; and ii) involvement of the entire Canadian community in CCAT, with a large enough share in the partnership for Canadian astronomers to participate at all levels of the facility. We further recommend continued participation in ALMA development, involvement in many focused PI-led submillimetre experiments, and partnership in SPICA.

preprint2013arXiv

Cold Dust in Hot Regions

We mapped five massive star forming regions with the SCUBA-2 camera on the James Clerk Maxwell Telescope (JCMT). Temperature and column density maps are obtained from the SCUBA-2 450 and 850 $μ$m images. Most of the dense clumps we find have central temperatures below 20 K with some as cold as 8 K, suggesting that they have no internal heating due to the presence of embedded protostars. This is surprising, because at the high densities inferred from these images and at these low temperatures such clumps should be unstable, collapsing to form stars and generating internal heating. The column densities at the clump centres exceed 10$^{23}$ cm$^{-2}$, and the derived peak visual extinction values are from 25-500 mag for $β$ = 1.5-2.5, indicating highly opaque centres. The observed cloud gas masses range from $\sim$ 10 to 10$^{3}$ M$_{\odot}$. The outer regions of the clumps follow an $r^{-2.36\pm0.35}$ density distribution and this power-law structure is observed outside of typically 10$^{4}$ AU. All these findings suggest that these clumps are high-mass starless clumps and most likely contain high-mass starless cores.

preprint2013arXiv

Cosmic Mnemonics

Our current description of the large-scale Universe is now known with a precision undreamt of a generation ago. Within the simple standard cosmological model only six basic parameters are required. The usual parameter set includes quantities most directly probed by the cosmic microwave background, but the nature of these quantities is somewhat esoteric. However, many more numbers can be derived that quantify various aspects of our Universe. Using constraints from the Planck satellite, in combination with other data sets, we explore several such quantities, highlighting some specific examples.

preprint2013arXiv

Disentangling a group of lensed submm galaxies at z~2.9

MS$\,$0451.6$-$0305 is a rich galaxy cluster whose strong lensing is particularly prominent at submm wavelengths. We combine new SCUBA-2 data with imaging from Herschel SPIRE and PACS and HST in order to try to understand the nature of the sources being lensed. In the region of the "giant submm arc," we uncover seven multiply imaged galaxies (up from the previously known three), of which six are found to be at a redshift of $z\sim2.9$, and possibly constitute an interacting system. Using a novel forward-modelling approach, we are able to simultaneously deblend and fit SEDs to the individual galaxies that contribute to the giant submm arc, constraining their dust temperatures, far infrared luminosities and star formation rates. The submm arc first identified by SCUBA can now be seen to be composed of at least five distinct sources, four of these within the galaxy group at $z\sim2.9$. The total unlensed luminosity for this galaxy group is $(3.1\pm0.3) \times 10^{12}\,\mathrm{L}_\odot$, which gives an unlensed star formation rate of $(450\pm50)$ M$_\odot$ yr$^{-1}$. From the properties of this system, we see no evidence of evolution towards lower temperatures in the dust temperature versus far-infrared luminosity relation for high redshift galaxies.

preprint2013arXiv

Empirical modelling of the BLASTPol achromatic half-wave plate for precision submillimetre polarimetry

A cryogenic achromatic half-wave plate (HWP) for submillimetre astronomical polarimetry has been designed, manufactured, tested, and deployed in the Balloon-borne Large-Aperture Submillimeter Telescope for Polarimetry (BLASTPol). The design is based on the five-slab Pancharatnam recipe and it works in the wavelength range 200-600 micron, making it the broadest-band HWP built to date at (sub)millimetre wavelengths. The frequency behaviour of the HWP has been fully characterised at room and cryogenic temperatures with incoherent radiation from a polarising Fourier transform spectrometer. We develop a novel empirical model, complementary to the physical and analytical ones available in the literature, that allows us to recover the HWP Mueller matrix and phase shift as a function of frequency and extrapolated to 4K. We show that most of the HWP non-idealities can be modelled by quantifying one wavelength-dependent parameter, the position of the HWP equivalent axes, which is then readily implemented in a map-making algorithm. We derive this parameter for a range of spectral signatures of input astronomical sources relevant to BLASTPol, and provide a benchmark example of how our method can yield improved accuracy on measurements of the polarisation angle on the sky at submillimetre wavelengths.

preprint2013arXiv

Far-IR Emission From Dust-Obscured Galaxies

Dust-obscured galaxies (DOGs) are a UV-faint, IR-bright galaxy population that reside at z~2 and are believed to be in a phase of dusty star-forming and AGN activity. We present far-IR observations of a complete sample of DOGs in the 2 deg2 of COSMOS. The 3077 DOGs have <z>=1.9+/-0.3 and are selected from 24um and r+ observations using a color cut of r+ - [24]>=7.5 (AB mag) and S24>=100uJy. Based on the near-IR SEDs, 47% are star-formation dominated and 10% are AGN-dominated. We use SPIRE far-IR photometry from HerMES to calculate the IR luminosity and characteristic dust temperature for the 1572 (51%) DOGs that are detected at 250um (>=3σ). For the remaining 1505 (49%) that are undetected, we perform a median stacking analysis to probe fainter luminosities. Detected and undetected DOGs have average IR luminosities of (2.8+/-0.4) x 1012 LSun and (0.77+-0.08)x10^12LSun, and dust temperatures of 34+/-7 K and 37+/-3 K, respectively. The IR luminosity function of DOGs with S24>=100uJy is calculated, using far-IR observations and stacking. DOGs contribute 10-30% to the total star formation rate density of the Universe at z=1.5-2.5, dominated by 250um detected and bump DOGs. For comparison, DOGs contribute 30% to the star-formation rate density of the Universe for all galaxies at z=1.5-2.5 with S24>=100uJy. DOGs have a large scatter about the star-formation main sequence and their specific star-formation rates show that the observed phase of star-formation could be responsible for their observed stellar mass at z~2.

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

HerMES: A Deficit in the Surface Brightness of the Cosmic Infrared Background Due to Galaxy Cluster Gravitational Lensing

We have observed four massive galaxy clusters with the SPIRE instrument on the Herschel Space Observatory and measure a deficit of surface brightness within their central region after subtracting sources. We simulate the effects of instrumental sensitivity and resolution, the source population, and the lensing effect of the clusters to estimate the shape and amplitude of the deficit. The amplitude of the central deficit is a strong function of the surface density and flux distribution of the background sources. We find that for the current best fitting faint end number counts, and excellent lensing models, the most likely amplitude of the central deficit is the full intensity of the cosmic infrared background (CIB). Our measurement leads to a lower limit to the integrated total intensity of the CIB of I(250 microns) > 0.69_(-0.03)^(+0.03) (stat.)_(-0.06)^(+0.11) (sys.) MJy/sr, with more CIB possible from both low-redshift sources and from sources within the target clusters. It should be possible to observe this effect in existing high angular resolution data at other wavelengths where the CIB is bright, which would allow tests of models of the faint source component of the CIB.

preprint2013arXiv

HerMES: Candidate High-Redshift Galaxies Discovered with Herschel/SPIRE

We present a method for selecting $z>4$ dusty, star forming galaxies (DSFGs) using Herschel/SPIRE 250/350/500 $μm$ flux densities to search for red sources. We apply this method to 21 deg$^2$ of data from the HerMES survey to produce a catalog of 38 high-$z$ candidates. Follow-up of the first 5 of these sources confirms that this method is efficient at selecting high-$z$ DSFGs, with 4/5 at $z=4.3$ to $6.3$ (and the remaining source at $z=3.4$), and that they are some of the most luminous dusty sources known. Comparison with previous DSFG samples, mostly selected at longer wavelengths (e.g., 850 $μm$) and in single-band surveys, shows that our method is much more efficient at selecting high-$z$ DSFGs, in the sense that a much larger fraction are at $z>3$. Correcting for the selection completeness and purity, we find that the number of bright ($S_{500\,μm} \ge 30$ mJy), red Herschel sources is $3.3 \pm 0.8$ deg$^{-2}$. This is much higher than the number predicted by current models, suggesting that the DSFG population extends to higher redshifts than previously believed. If the shape of the luminosity function for high-$z$ DSFGs is similar to that at $z\sim2$, rest-frame UV based studies may be missing a significant component of the star formation density at $z=4$ to $6$, even after correction for extinction.

preprint2013arXiv

HerMES: deep number counts at 250, 350, and 500 microns in the COSMOS and GOODS-N fields and the build-up of the cosmic infrared background

ABRIGED Herschel/SPIRE has provided confusion limited maps of deep fields at 250, 350, and 500um, as part of the HerMES survey. Due to confusion, only a small fraction of the Cosmic Infrared Background can be resolved into individually-detected sources. Our goal is to produce deep galaxy number counts and redshift distributions below the confusion limit, which we then use to place strong constraints on the origins of the cosmic infrared background and on models of galaxy evolution. We individually extracted the bright SPIRE with a method using the positions, the flux densities, and the redshifts of the 24um sources as a prior, and derived the number counts and redshift distributions of the bright SPIRE sources. For fainter SPIRE sources, we reconstructed the number counts and the redshift distribution below the confusion limit using the deep 24um catalogs associated with photometric redshift and information provided by the stacking of these sources into the deep SPIRE maps. Finally, by integrating all these counts, we studied the contribution of the galaxies to the CIB as a function of their flux density and redshift. Through stacking, we managed to reconstruct the source counts per redshift slice down to ~2 mJy in the three SPIRE bands, which lies about a factor 10 below the 5sigma confusion limit. None of the pre-existing population models are able to reproduce our results at better than 3sigma. Finally, we extrapolate our counts to zero flux density in order to derive an estimate of the total contribution of galaxies to the CIB, finding 10.1, 6.5, and 2.8 nW/m2/sr at 250, 350, and 500um, respectively. These values agree well with FIRAS absolute measurements, suggesting our number counts and their extrapolation are sufficient to explain the CIB. Finally, combining our results with other works, we estimate the energy budget contained in the CIB between 8 and 1000um: 26 nW/m2/sr.

preprint2013arXiv

Independent constraints on local non-Gaussianity from the peculiar velocity and density fields

Primordial, non-Gaussian perturbations can generate scale-dependent bias in the galaxy distribution. This in turn will modify correlations between galaxy positions and peculiar velocities at late times, since peculiar velocities reflect the underlying matter distribution, whereas galaxies are a biased tracer of the same. We study this effect, and show that non-Gaussianity can be constrained by comparing the observed peculiar velocity field to a model velocity field reconstructed from the galaxy density field assuming linear bias. The amplitude of the spatial correlations in the residual map obtained after subtracting one velocity field from the other is directly proportional to the strength of the primordial non-Gaussianity. We construct the corresponding likelihood function use it to constrain the amplitude of the linear flow $β$ and the amplitude of local non-Gaussianity $f^{\rm NL}_{\rm local}$. Applying our method to two observational data sets, the Type-Ia supernovae (A1SN) and Spiral Field \textit{I}-band (SFI++) catalogues, we obtain constraints on the linear flow parameter consistent with the values derived previously assuming Gaussianity. The marginalised 1-D distribution of $|f^{\rm NL}_{\rm local}|$ does not show strong evidence for non-zero $f^{\rm NL}_{\rm local}$, and we set 95% upper limits $|f^{\rm NL}_{\rm local}|<51.4$ from A1SN and $|f^{\rm NL}_{\rm local}|<92.6$ from SFI++. These limits on $f^{\rm NL}_{\rm local}$ are as tight as any set by previous large-scale structure measurements. Our method can be applied to any survey with radial velocities and density field data, and provides an independent check of recent CMB constraints on $f^{\rm NL}_{\rm local}$, extending these to smaller spatial scales.

preprint2013arXiv

Lupus I Observations from the 2010 Flight of the Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry

The Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol) was created by adding polarimetric capability to the BLAST experiment that was flown in 2003, 2005, and 2006. BLASTPol inherited BLAST's 1.8 m primary and its Herschel/SPIRE heritage focal plane that allows simultaneous observation at 250, 350, and 500 μm. We flew BLASTPol in 2010 and again in 2012. Both were long duration Antarctic flights. Here we present polarimetry of the nearby filamentary dark cloud Lupus I obtained during the 2010 flight. Despite limitations imposed by the effects of a damaged optical component, we were able to clearly detect submillimeter polarization on degree scales. We compare the resulting BLASTPol magnetic field map with a similar map made via optical polarimetry (The optical data were published in 1998 by J. Rizzo and collaborators.). The two maps partially overlap and are reasonably consistent with one another. We compare these magnetic field maps to the orientations of filaments in Lupus I, and we find that the dominant filament in the cloud is approximately perpendicular to the large-scale field, while secondary filaments appear to run parallel to the magnetic fields in their vicinities. This is similar to what is observed in Serpens South via near-IR polarimetry, and consistent with what is seen in MHD simulations by F. Nakamura and Z. Li.

preprint2013arXiv

Minimal parameterizations for modified gravity

The increasing precision of cosmological data provides us with an opportunity to test general relativity (GR) on the largest accessible scales. Parameterizing modified gravity models facilitates the systematic testing of the predictions of GR, and gives a framework for detecting possible deviations from it. Several different parameterizations have already been suggested, some linked to classifications of theories, and others more empirically motivated. Here we describe a particular new approach which casts modifications to gravity through two free functions of time and scale, which are directly linked to the field equations, but also easy to confront with observational data. We compare our approach with other existing methods of parameterizing modied gravity, specifically the parameterized post-Friedmann approach and the older method using the parameter set $\{μ,γ\}$. We explain the connection between our parameters and the physics that is most important for generating cosmic microwave background anisotropies. Some qualitative features of this new parameterization, and therefore modifications to the gravitational equations of motion, are illustrated in a toy model, where the two functions are simply assumed to be constant parameters.

preprint2013arXiv

Panchromatic spectral energy distributions of Herschel sources

(abridged) Far-infrared Herschel photometry from the PEP and HerMES programs is combined with ancillary datasets in the GOODS-N, GOODS-S, and COSMOS fields. Based on this rich dataset, we reproduce the restframe UV to FIR ten-colors distribution of galaxies using a superposition of multi-variate Gaussian modes. The median SED of each mode is then fitted with a modified version of the MAGPHYS code that combines stellar light, emission from dust heated by stars and a possible warm dust contribution heated by an AGN. The defined Gaussian grouping is also used to identify rare sources. The zoology of outliers includes Herschel-detected ellipticals, very blue z~1 Ly-break galaxies, quiescent spirals, and torus-dominated AGN with star formation. Out of these groups and outliers, a new template library is assembled, consisting of 32 SEDs describing the intrinsic scatter in the restframe UV-to-submm colors of infrared galaxies. This library is tested against L(IR) estimates with and without Herschel data included, and compared to eight other popular methods often adopted in the literature. When implementing Herschel photometry, these approaches produce L(IR) values consistent with each other within a median absolute deviation of 10-20%, the scatter being dominated more by fine tuning of the codes, rather than by the choice of SED templates. Finally, the library is used to classify 24 micron detected sources in PEP GOODS fields. AGN appear to be distributed in the stellar mass (M*) vs. star formation rate (SFR) space along with all other galaxies, regardless of the amount of infrared luminosity they are powering, with the tendency to lie on the high SFR side of the "main sequence". The incidence of warmer star-forming sources grows for objects with higher specific star formation rates (sSFR), and they tend to populate the "off-sequence" region of the M*-SFR-z space.

preprint2013arXiv

SCUBA-2: iterative map-making with the Sub-Millimetre User Reduction Facility

The Submillimetre Common User Bolometer Array 2 (SCUBA-2) is an instrument operating on the 15-m James Clerk Maxwell Telescope, nominally consisting of 5120 bolometers in each of two simultaneous imaging bands centred over 450 and 850 um. The camera is operated by scanning across the sky and recording data at a rate of 200 Hz. As the largest of a new generation of multiplexed kilopixel bolometer cameras operating in the (sub)millimetre, SCUBA-2 data analysis represents a significant challenge. We describe the production of maps using the Sub-Millimetre User Reduction Facility (SMURF) in which we have adopted a fast, iterative approach to map-making that enables data reduction on single, modern, high-end desktop computers, with execution times that are typically shorter than the observing times. SMURF is used in an automated setting, both at the telescope for real-time feedback to observers, as well as for the production of science products for the JCMT Science Archive at the Canadian Astronomy Data Centre. Three detailed case studies are used to: (i) explore convergence properties of the map-maker using simple prior constraints (Uranus -- a point source); (ii) achieve the white-noise limit for faint point-source studies (extragalactic blank-field survey of the Lockman Hole); and (iii) demonstrate that our strategy is capable of recovering angular scales comparable to the size of the array footprint (approximately 5 arcmin) for bright extended sources (star-forming region M17).

preprint2013arXiv

The Herschel census of infrared SEDs through cosmic time

Using Herschel data from the deepest SPIRE and PACS surveys (HerMES and PEP) in COSMOS and GOODS (N+S), we examine the dust properties of IR-luminous (L_IR>10^10 L_sun) galaxies at 0.1<z<2 and determine how these evolve with cosmic time. The unique angle of this work is the rigorous analysis of survey selection effects, making this the first study of the star-formation-dominated, IR-luminous population within a framework almost entirely free of selection biases. We find that IR-luminous galaxies have SEDs with broad far-IR peaks characterised by cool/extended dust emission and average dust temperatures in the 25-45K range. Hot (T>45K) SEDs and cold (T<25K), cirrus-dominated SEDs are rare, with most sources being within the range occupied by warm starbursts such as M82 and cool spirals such as M51. We observe a luminosity-temperature (L-T) relation, where the average dust temperature of log[L_IR/L_sun]=12.5 galaxies is about 10K higher than that of their log[L_IR/L_sun]=10.5 counterparts. However, although the increased dust heating in more luminous systems is the driving factor behind the L-T relation, the increase in dust mass and/or starburst size with luminosity plays a dominant role in shaping it. Our results show that the dust conditions in IR-luminous sources evolve with cosmic time: at high redshift, dust temperatures are on average up to 10K lower than what is measured locally. This is manifested as a flattening of the L-T relation, suggesting that (U)LIRGs in the early Universe are typically characterised by a more extended dust distribution and/or higher dust masses than local equivalent sources. Interestingly, the evolution in dust temperature is luminosity dependent, with the fraction of LIRGs with T<35K showing a 2-fold increase from z~0 to z~2, whereas that of ULIRGs with T<35K shows a 6-fold increase.

preprint2013arXiv

The SCUBA-2 Cosmology Legacy Survey: demographics of the 450 μm-population

We investigate the multi-wavelength properties of a sample of 450-μm selected sources from the SCUBA-2 Cosmology Legacy Survey (S2CLS). A total of 69 sources were identified above 4σ in deep SCUBA-2 450-μm observations overlapping the UDS and COSMOS fields and covering 210 sq. arcmin to a typical depth of σ450=1.5 mJy. Reliable cross identification are found for 58 sources (84 per cent) in Spitzer and Hubble Space Telescope WFC3/IR data. The photometric redshift distribution (dN/dz) of 450μm-selected sources is presented, showing a broad peak in the redshift range 1<z<3, and a median of z=1.4. Combining the SCUBA-2 photometry with Herschel SPIRE data from HerMES, the submm spectral energy distribution (SED) is examined via the use of modified blackbody fits, yielding aggregate values for the IR luminosity, dust temperature and emissivity of <LIR>=10^12 +/- 0.8 L_sol, <T_D>=42 +/- 11 K and <β_D>=1.6 +/- 0.5, respectively. The relationship between these SED parameters and the physical properties of galaxies is investigated, revealing correlations between T_D and LIR and between β_D and both stellar mass and effective radius. The connection between star formation rate and stellar mass is explored, with 24 per cent of 450 μm sources found to be ``star-bursts'', i.e. displaying anomalously high specific SFRs. However, both the number density and observed properties of these ``star-burst'' galaxies are found consistent with the population of normal star-forming galaxies.

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 comparison of the galaxy peculiar velocity field with the PSCz gravity field-- A Bayesian hyper-parameter method

We constructed a Bayesian hyper-parameter statistical method to quantify the difference between predicted velocities derived from the observed galaxy distribution in the \textit{IRAS}-PSC$z$ redshift survey and peculiar velocities measured using different distance indicators. In our analysis we find that the model--data comparison becomes unreliable beyond $70 \hmpc$ because of the inadequate sampling by \textit{IRAS} survey of prominent, distant superclusters, like the Shapley Concentration. On the other hand, the analysis of the velocity residuals show that the PSC$z$ gravity field provides an adequate model to the local, $\le 70 \hmpc$, peculiar velocity field. The hyper-parameter combination of ENEAR, SN, A1SN and SFI++ catalogues in the Bayesian framework constrains the amplitude of the linear flow to be $β=0.53 \pm 0.014$. For an rms density fluctuations in the PSC$z$ galaxy number density $σ_8^{\rm gal}=0.42\pm0.03$, we obtain an estimate of the growth rate of density fluctuations $fσ_{8}(z\sim0) = 0.42 \pm 0.033$, which is in excellent agreement with independent estimates based on different techniques.

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

A new era of wide-field submillimetre imaging: on-sky performance of SCUBA-2

SCUBA-2 is the largest submillimetre wide-field bolometric camera ever built. This 43 square arc-minute field-of-view instrument operates at two wavelengths (850 and 450 microns) and has been installed on the James Clerk Maxwell Telescope on Mauna Kea, Hawaii. SCUBA-2 has been successfully commissioned and operational for general science since October 2011. This paper presents an overview of the on-sky performance of the instrument during and since commissioning in mid-2011. The on-sky noise characteristics and NEPs of the 450 and 850 micron arrays, with average yields of approximately 3400 bolometers at each wavelength, will be shown. The observing modes of the instrument and the on-sky calibration techniques are described. The culmination of these efforts has resulted in a scientifically powerful mapping camera with sensitivities that allow a square degree of sky to be mapped to 10 mJy/beam rms at 850 micron in 2 hours and 60 mJy/beam rms at 450 micron in 5 hours in the best weather.

preprint2012arXiv

A Population of Dust-Rich Quasars at z ~ 1.5

We report Herschel SPIRE (250, 350, and 500 micron) detections of 32 quasars with redshifts 0.5 < z < 3.6 from the Herschel Multi-tiered Extragalactic Survey. These sources are from a MIPS 24 micron flux-limited sample of 326 quasars in the Lockman Hole Field. The extensive multi-wavelength data available in the field permit construction of the rest-frame Spectral Energy Distributions (SEDs)from ultraviolet to the mid-infrared for all sources, and to the far-infrared (FIR) for the 32 objects. Most quasars with Herschel FIR detections show dust temperatures in the range of 25K to 60K, with a mean of 34K. The FIR luminosities range from 10^{11.3} to 10^{13.5} Lsun, qualifying most of their hosts as ultra- or hyper-luminous infrared galaxies. These FIR-detected quasars may represent a dust-rich population, but with lower redshifts and fainter luminosities than quasars observed at ~ 1 mm. However, their FIR properties cannot be predicted from shorter wavelengths (0.3--20 micron, rest-frame), and the bolometric luminosities derived using the 5100 A index may be underestimated for these FIR-detected quasars. Regardless of redshift, we observed a decline in the relative strength of FIR luminosities for quasars with higher near-infrared luminosities.

preprint2012arXiv

A unified empirical model for infrared galaxy counts based on observed physical evolution of distant galaxies

We reproduce the mid-infrared to radio galaxy counts with a new empirical model based on our current understanding of the evolution of main-sequence (MS) and starburst (SB) galaxies. We rely on a simple Spectral Energy Distribution (SED) library based on Herschel observations: a single SED for the MS and another one for SB, getting warmer with redshift. Our model is able to reproduce recent measurements of galaxy counts performed with Herschel, including counts per redshift slice. This agreement demonstrates the power of our 2 Star-Formation Modes (2SFM) decomposition for describing the statistical properties of infrared sources and their evolution with cosmic time. We discuss the relative contribution of MS and SB galaxies to the number counts at various wavelengths and flux densities. We also show that MS galaxies are responsible for a bump in the 1.4 GHz radio counts around 50 μJy. Material of the model (predictions, SED library, mock catalogs...) is available online at http://irfu.cea.fr/Sap/Phocea/Page/index.php?id=537.

preprint2012arXiv

Clustering properties of BzK-selected galaxies in GOODS-N: environmental quenching and triggering of star formation at z ~ 2

Using a sample of BzK-selected galaxies at z~2 identified from the CFHT/WIRCAM near-infrared survey of GOODS-North, we discuss the relation between star formation rate (SFR), specific star formation rate (SSFR), and stellar mass (M_{*}), and the clustering of galaxies as a function of these parameters. For star-forming galaxies (sBzKs), the UV-based SFR, corrected for extinction, scales with the stellar mass as SFR ~ M_{*}^{alpha} with alpha = 0.74+/-0.20 down to M_{*} ~ 10^{9} M_{solar}, indicating a weak dependence on the stellar mass of the SSFR. We also measure the angular correlation function and hence infer the correlation length for sBzK galaxies as a function of M_{*}, SFR, and SSFR, as well as K-band apparent magnitude. We show that passive galaxies (pBzKs) are more strongly clustered than sBzK galaxies at a given stellar mass, mirroring the color-density relation seen at lower redshifts. We also find that the correlation length of sBzK galaxies ranges from 4 to 20 h^{-1}Mpc, being a strong function of M_{K}, M_{*}, and SFR. On the other hand, the clustering dependence on SSFR changes abruptly at 2x10^{-9} yr^{-1}, which is the typical value for "main sequence" star-forming galaxies at z~2. We show that the correlation length reaches a minimum at this characteristic value, and is larger for galaxies with both smaller and larger SSFRs; a dichotomy that is only marginally implied from the predictions of the semi-analytical models. Our results suggest that there are two types of environmental effects at work at z~2. Stronger clustering for relatively quiescent galaxies implies that the environment has started to play a role in quenching star formation. At the same time, stronger clustering for galaxies with elevated SSFRs ("starbursts") might be attributed to an increased efficiency for galaxy interactions and mergers in dense environments.

preprint2012arXiv

Cosmic bulk flows on 50 {h}^{-1}Mpc scales: A Bayesian hyper-parameter method and multishells likelihood analysis

It has been argued recently that the galaxy peculiar velocity field provides evidence of excessive power on scales of $50\hmpc$, which seems to be inconsistent with the standard $Λ$CDM cosmological model. We discuss several assumptions and conventions used in studies of the large-scale bulk flow to check whether this claim is robust under a variety of conditions. Rather than using a composite catalogue we select samples from the SN, ENEAR, SFI++ and A1SN catalogues, and correct for Malmquist bias in each according to the IRAS PSCz density field. We also use slightly different assumptions about the small-scale velocity dispersion and the parameterisation of the matter power spectrum when calculating the variance of the bulk flow. By combining the likelihood of individual catalogues using a Bayesian hyper-parameter method, we find that the joint likelihood of the amplitude parameter gives $σ_8=0.65^{+0.47}_{-0.35}(\pm 1 σ)$, which is entirely consistent with the $Λ$CDM model. In addition, the bulk flow magnitude ($v \sim 310 \kms$) and direction, $(l,b)\sim (280^{\circ} \pm 8^{\circ}, 5.1^{\circ} \pm 6^{\circ})$, found by each of the catalogues are all consistent with each other, and with the bulk flow results from most previous studies. Furthermore, the bulk flow velocities in different shells of the surveys constrain $σ_{8}$--$Ω_{\rm{m}}$ to be ($1.01^{+0.26}_{-0.20},0.31^{+0.28}_{-0.14}$, SFI++) and ($1.04^{+0.32}_{-0.24},0.28^{+0.30}_{-0.14}$, ENEAR), which is consistent with {\it WMAP} 7-year best-fit values. We finally discuss the differences between our conclusions and those of the studies claiming the largest bulk flows.

preprint2012arXiv

Dimensionless cosmology

Although it is well known that any consideration of the variations of fundamental constants should be restricted to their dimensionless combinations, the literature on variations of the gravitational constant $G$ is entirely dimensionful. To illustrate applications of this to cosmology, we explicitly give a dimensionless version of the parameters of the standard cosmological model, and describe the physics of Big Bang Neucleosynthesis and recombination in a dimensionless manner. The issue that appears to have been missed in many studies is that in cosmology the strength of gravity is bound up in the cosmological equations, and the epoch at which we live is a crucial part of the model. We argue that it is useful to consider the hypothetical situation of communicating with another civilization (with entirely different units), comparing only dimensionless constants, in order to decide if we live in a Universe governed by precisely the same physical laws. In this thought experiment, we would also have to compare epochs, which can be defined by giving the value of any {\it one} of the evolving cosmological parameters. By setting things up carefully in this way one can avoid inconsistent results when considering variable constants, caused by effectively fixing more than one parameter today. We show examples of this effect by considering microwave background anisotropies, being careful to maintain dimensionlessness throughout. We present Fisher matrix calculations to estimate how well the fine structure constants for electromagnetism and gravity can be determined with future microwave background experiments. We highlight how one can be misled by simply adding $G$ to the usual cosmological parameter set.

preprint2012arXiv

Evidence for a wide range of UV obscuration in z ~ 2 dusty galaxies from the GOODS-Herschel survey

Dusty galaxies at z ~ 2 span a wide range of relative brightness between rest-frame mid-infrared (8um) and ultraviolet wavelengths. We attempt to determine the physical mechanism responsible for this diversity. Dust-obscured galaxies (DOGs), which have rest-frame mid-IR to UV flux density ratios > 1000, might be abnormally bright in the mid-IR, perhaps due to prominent AGN and/or PAH emission, or abnormally faint in the UV. We use far-infrared data from the GOODS-Herschel survey to show that most DOGs with 10^12 L_Sun < L_IR < 10^13 L_Sun are not abnormally bright in the mid-IR when compared to other dusty galaxies with similar IR (8--1000um) luminosities. We observe a relation between the median IR to UV luminosity ratios and the median UV continuum power-law indices for these galaxies, and we find that only 24% have specific star formation rates which indicate the dominance of compact star-forming regions. This circumstantial evidence supports the idea that the UV- and IR-emitting regions in these galaxies are spatially coincident, which implies a connection between the abnormal UV faintness of DOGs and dust obscuration. We conclude that the range in rest-frame mid-IR to UV flux density ratios spanned by dusty galaxies at z ~ 2 is due to differing amounts of UV obscuration. Of galaxies with these IR luminosities, DOGs are the most obscured. We attribute differences in UV obscuration to either: 1) differences in the degree of alignment between the spatial distributions of dust and massive stars, or 2) differences in the total dust content.

preprint2012arXiv

Evidence for Environmental Changes in the Submillimeter Dust Opacity

The submillimeter opacity of dust in the diffuse Galactic interstellar medium (ISM) has been quantified using a pixel-by-pixel correlation of images of continuum emission with a proxy for column density. We used three BLAST bands at 250, 350, and 500 μm and one IRAS at 100 μm. The proxy is the near-infrared color excess, E(J-Ks), obtained from 2MASS. Based on observations of stars, we show how well this color excess is correlated with the total hydrogen column density for regions of moderate extinction. The ratio of emission to column density, the emissivity, is then known from the correlations, as a function of frequency. The spectral distribution of this emissivity can be fit by a modified blackbody, whence the characteristic dust temperature T and the desired opacity σ_e(1200) at 1200 GHz can be obtained. We have analyzed 14 regions near the Galactic plane toward the Vela molecular cloud, mostly selected to avoid regions of high column density (N_H > 10^{22} cm^-2) and small enough to ensure a uniform T. We find σ_e(1200) is typically 2 to 4 x 10^{-25} cm^2/H and thus about 2 to 4 times larger than the average value in the local high Galactic latitude diffuse atomic ISM. This is strong evidence for grain evolution. There is a range in total power per H nucleon absorbed (re-radiated) by the dust, reflecting changes in the interstellar radiation field and/or the dust absorption opacity. These changes affect the equilibrium T, which is typically 15 K, colder than at high latitudes. Our analysis extends, to higher opacity and lower T, the trend of increasing opacity with decreasing T that was found at high latitudes. The recognition of changes in the emission opacity raises a cautionary flag because all column densities deduced from dust emission maps, and the masses of compact structures within them, depend inversely on the value adopted.

preprint2012arXiv

FMOS near-IR spectroscopy of Herschel selected galaxies: star formation rates, metallicity and dust attenuation at z~1

We investigate the properties (e.g. star formation rate, dust attentuation, stellar mass and metallicity) of a sample of infrared luminous galaxies at z \sim 1 via near-IR spectroscopy with Subaru-FMOS. Our sample consists of Herschel SPIRE and Spitzer MIPS selected sources in the COSMOS field with photometric redshifts in the range 0.7 < z-phot < 1.8, which have been targeted in 2 pointings (0.5 sq. deg.) with FMOS. We find a modest success rate for emission line detections, with candidate Hα emission lines detected for 57 of 168 SPIRE sources (34 per cent). By stacking the near-IR spectra we directly measure the mean Balmer decrement for the Hα and Hβ lines, finding a value of <E(B-V)> = 0.51\pm0.27 for <LIR> = 10^12 Lsol sources at <z> = 1.36. By comparing star formation rates estimated from the IR and from the dust uncorrected Hα line we find a strong relationship between dust attenuation and star formation rate. This relation is broadly consistent with that previously seen in star-forming galaxies at z ~ 0.1. Finally, we investigate the metallicity via the N2 ratio, finding that z ~ 1 IR-selected sources are indistinguishable from the local mass-metallicity relation. We also find a strong correlation between dust attentuation and metallicity, with the most metal-rich IR-sources experiencing the largest levels of dust attenuation.

preprint2012arXiv

GOODS-Herschel: Impact of Active Galactic Nuclei and Star Formation Activity on Infrared Spectral Energy Distributions at High Redshift

We explore the effects of active galactic nuclei (AGN) and star formation activity on the infrared (0.3-1000 microns) spectral energy distributions of luminous infrared galaxies from z = 0.5 to 4.0. We have compiled a large sample of 151 galaxies selected at 24 microns (S24 > 100 uJy) in the GOODS-N and ECDFS fields for which we have deep Spitzer IRS spectroscopy, allowing us to decompose the mid-IR spectrum into contributions from star formation and AGN activity. A significant portion (~25%) of our sample is dominated by an AGN in the mid-IR. Based on the mid-IR classification, we divide our full sample into four sub-samples: z~1 star-forming (SF) sources; z~2 SF sources; AGN with clear 9.7 micron silicate absorption; and AGN with featureless mid-IR spectra. From our large spectroscopic sample and wealth of multi-wavelength data, including deep Herschel imaging at 100, 160, 250, 350, and 500 microns, we use 95 galaxies with complete spectral coverage to create a composite spectral energy distribution (SED) for each sub-sample. We then fit a two-temperature component modified blackbody to the SEDs. We find that the IR SEDs have similar cold dust temperatures, regardless of the mid-IR power source, but display a marked difference in the warmer dust temperatures. We calculate the average effective temperature of the dust in each sub-sample and find a significant (~20 K) difference between the SF and AGN systems. We compare our composite SEDs to local templates and find that local templates do not accurately reproduce the mid-IR features and dust temperatures of our high redshift systems. High redshift IR luminous galaxies contain significantly more cool dust than their local counterparts. We find that a full suite of photometry spanning the IR peak is necessary to accurately account for the dominant dust temperature components in high redshift IR luminous galaxies.

preprint2012arXiv

HerMES: Unveiling obscured star formation - the far infrared luminosity function of ultraviolet-selected galaxies at z~1.5

We study the far-infrared (IR) and sub-millimeter properties of a sample of ultraviolet (UV) selected galaxies at z\sim1.5. Using stacking at 250, 350 and 500 um from Herschel Space Observatory SPIRE imaging of the COSMOS field obtained within the HerMES key program, we derive the mean IR luminosity as a function of both UV luminosity and slope of the UV continuum beta. The IR to UV luminosity ratio is roughly constant over most of the UV luminosity range we explore. We also find that the IR to UV luminosity ratio is correlated with beta. We observe a correlation that underestimates the correlation derived from low-redshift starburst galaxies, but is in good agreement with the correlation derived from local normal star-forming galaxies. Using these results we reconstruct the IR luminosity function of our UV-selected sample. This luminosity function recovers the IR luminosity functions measured from IR selected samples at the faintest luminosities (Lir ~ 10^{11} L_sun), but might underestimate them at the bright-end (Lir > 5.10^{11} L_sun). For galaxies with 10^{11}<Lir/L_sun<10^{13}, the IR luminosity function of a UV selection recovers (given the differences in IR-based estimates) 52-65 to 89-112 per cent of the star-formation rate density derived from an IR selection. The cosmic star-formation rate density derived from this IR luminosity function is 61-76 to 100-133 per cent of the density derived from IR selections at the same epoch. Assuming the latest Herschel results and conservative stacking measurements, we use a toy model to fully reproduce the far IR luminosity function from our UV selection at z\sim 1.5. This suggests that a sample around 4 magnitudes deeper (i.e. reaching u \sim 30 mag) and a large dispersion of the IR to UV luminosity ratio are required.

preprint2012arXiv

Herschel observations of a z~2 stellar mass selected galaxy sample drawn from the GOODS NICMOS Survey

We present a study of the far-IR properties of a stellar mass selected sample of 1.5 < z < 3 galaxies with log(M_*/M_sun) > 9.5 drawn from the GOODS NICMOS Survey (GNS), the deepest H-band Hubble Space Telescope survey of its type prior to the installation of WFC3. We use far-IR and sub-mm data from the PACS and SPIRE instruments on-board Herschel, taken from the PACS Evolutionary Probe (PEP) and Herschel Multi-Tiered Extragalactic Survey (HerMES) key projects respectively. We find a total of 22 GNS galaxies, with median log(M_*/M_sun) = 10.8 and z = 2.0, associated with 250 um sources detected with SNR > 3. We derive mean total IR luminosity log L_IR (L_sun) = 12.36 +/- 0.05 and corresponding star formation rate SFR_(IR+UV) = (280 +/- 40) M_sun/yr for these objects, and find them to have mean dust temperature T_dust ~ 35 K. We find that the SFR derived from the far-IR photometry combined with UV-based estimates of unobscured SFR for these galaxies is on average more than a factor of 2 higher than the SFR derived from extinction corrected UV emission alone, although we note that the IR-based estimate is subject to substantial Malmquist bias. To mitigate the effect of this bias and extend our study to fainter fluxes, we perform a stacking analysis to measure the mean SFR in bins of stellar mass. We obtain detections at the 2-4 sigma level at SPIRE wavelengths for samples with log(M_*/M_sun) > 10. In contrast to the Herschel detected GNS galaxies, we find that estimates of SFR_(IR+UV) for the stacked samples are comparable to those derived from extinction corrected UV emission, although the uncertainties are large. We find evidence for an increasing fraction of dust obscured star formation with stellar mass, finding SFR_IR/SFR_UV \propto M_*^{0.7 +/- 0.2}, which is likely a consequence of the mass--metallicity relation.

preprint2012arXiv

Power-Law Template for IR Point Source Clustering

We perform a combined fit to angular power spectra of unresolved infrared (IR) point sources from the Planck satellite (at 217, 353, 545 and 857 GHz, over angular scales 100 < l < 2200), the Balloon-borne Large-Aperture Submillimeter Telescope (BLAST; 250, 350 and 500 um; 1000 < l < 9000), and from correlating BLAST and Atacama Cosmology Telescope (ACT; 148 and 218 GHz) maps. We find that the clustered power over the range of angular scales and frequencies considered is well fit by a simple power law of the form C_l \propto l^-n with n = 1.25 +/- 0.06. While the IR sources are understood to lie at a range of redshifts, with a variety of dust properties, we find that the frequency dependence of the clustering power can be described by the square of a modified blackbody, nu^beta B(nu,T_eff), with a single emissivity index beta = 2.20 +/- 0.07 and effective temperature T_eff = 9.7 K. Our predictions for the clustering amplitude are consistent with existing ACT and South Pole Telescope results at around 150 and 220 GHz, as is our prediction for the effective dust spectral index, which we find to be alpha_150-220 = 3.68 +/- 0.07 between 150 and 220 GHz. Our constraints on the clustering shape and frequency dependence can be used to model the IR clustering as a contaminant in Cosmic Microwave Background anisotropy measurements. The combined Planck and BLAST data also rule out a linear bias clustering model.

preprint2012arXiv

The complex physics of dusty star-forming galaxies at high redshifts as revealed by Herschel and Spitzer

We combine far-infrared photometry from Herschel (PEP/HERMES) with deep mid-infrared spectroscopy from Spitzer to investigate the nature and the mass assembly history of a sample of 31 Luminous and Ultraluminous Infrared Galaxies at z~1 and 2 selected in GOODS-S with 24 $μ$m fluxes between 0.2 and 0.5 mJy. We model the data with a self-consistent physical model (GRASIL) which includes a state-of-the-art treatment of dust extinction and reprocessing. We find that all of our galaxies appear to require massive populations of old (>1 Gyr) stars and, at the same time, to host a moderate ongoing activity of SF (SFR < 100 M$_{\odot}$/yr). The bulk of the stars appear to have been formed a few Gyr before the observation in essentially all cases. Only five galaxies of the sample require a recent starburst superimposed on a quiescent star formation history (SFH). We also find discrepancies between our results and those based on optical-only SED fitting for the same objects; by fitting their observed Spectral Energy Distributions with our physical model we find higher extinctions (by $Δ$A_{V} ~ 0.81 and 1.14) and higher stellar masses (by $Δ$Log(M*) ~ 0.16 and 0.36 dex) for z~1 and z~2 (U)LIRGs, respectively. The stellar mass difference is larger for the most dust obscured objects. We also find lower SFRs than those computed from L_{IR} using the Kennicutt relation due to the significant contribution to the dust heating by intermediate-age stellar populations through 'cirrus' emission (~73% and ~66% of total L_{IR} for z~1 and z~2 (U)LIRGs, respectively).

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.

preprint2012arXiv

What do gas-rich galaxies actually tell us about modified Newtonian dynamics?

It has recently been claimed that measurements of the baryonic Tully-Fisher relation (BTFR), a power-law relationship between the observed baryonic masses and outer rotation velocities of galaxies, support the predictions of modified Newtonian dynamics for the slope and scatter in the relation, while challenging the cold dark matter (CDM) paradigm. We investigate these claims, and find that: 1) the scatter in the data used to determine the BTFR is in conflict with observational uncertainties on the data; 2) these data do not make strong distinctions regarding the best-fit BTFR parameters; 3) the literature contains a wide variety of measurements of the BTFR, many of which are discrepant with the recent results; and 4) the claimed CDM "prediction" for the BTFR is a gross oversimplification of the complex galaxy-scale physics involved. We conclude that the BTFR is currently untrustworthy as a test of CDM.

preprint2011arXiv

A Monte Carlo Approach to Evolution of the Far-Infrared Luminosity Function with BLAST

We constrain the evolution of the rest-frame far-infrared (FIR) luminosity function out to high redshift, by combining several pieces of complementary information provided by the deep Balloon-borne Large-Aperture Submillimeter Telescope surveys at 250, 350 and 500 micron, as well as other FIR and millimetre data. Unlike most other phenomenological models, we characterise the uncertainties in our fitted parameters using Monte Carlo Markov Chains. We use a bivariate local luminosity function that depends only on FIR luminosity and 60-to-100 micron colour, along with a single library of galaxy spectral energy distributions indexed by colour, and apply simple luminosity and density evolution. We use the surface density of sources, Cosmic Infrared Background (CIB) measurements and redshift distributions of bright sources, for which identifications have been made, to constrain this model. The precise evolution of the FIR luminosity function across this crucial range has eluded studies at longer wavelengths (e.g., using SCUBA and MAMBO) and at shorter wavelengths (e.g., Spitzer), and should provide a key piece of information required for the study of galaxy evolution. Our adoption of Monte Carlo methods enables us not only to find the best-fit evolution model, but also to explore correlations between the fitted parameters. Our model-fitting approach allows us to focus on sources of tension coming from the combination of data-sets. We specifically find that our choice of parameterisation has difficulty fitting the combination of CIB measurements and redshift distribution of sources near 1 mm. Existing and future data sets will be able to dramatically improve the fits, as well as break strong degeneracies among the models. [abridged]

preprint2011arXiv

Contribution to the Diffuse Radio Background from Extragalactic Radio Sources

We examine the brightness of the Cosmic Radio Background (CRB) by comparing the contribution from individual source counts to absolute measurements. We use a compilation of radio counts to estimate the contribution of detected sources to the CRB in several different frequency bands.We apply a Monte Carlo Markov Chain technique to estimate the brightness values and uncertainties, paying attention to various sources of systematic error. We compare our results to absolute measurements from the ARCADE 2 experiment. At v = 150 MHz, 325 MHz, 408 MHz, 610 MHz, 1.4 GHz, 4.8 GHz, and 8.4 GHz our calculated contributions to the background sky temperature are 18, 2.8, 1.6, 0.71, 0.11, 0.0032, 0.0059 K, respectively. If the ARCADE 2 measurements are correct and come from sources, then there must be an additional population of radio galaxies, fainter than where current data are probing. More specifically, the Euclidean-normalized counts at 1.4 GHz have to have an additional bump below about 10 μJy.

preprint2011arXiv

Correlations in the (Sub)millimeter background from ACTxBLAST

We present measurements of the auto- and cross-frequency correlation power spectra of the cosmic (sub)millimeter background at: 250, 350, and 500 um (1200, 860, and 600 GHz) from observations made with the Balloon-borne Large Aperture Submillimeter Telescope, BLAST; and at 1380 and 2030 um (218 and 148 GHz) from observations made with the Atacama Cosmology Telescope, ACT. The overlapping observations cover 8.6 deg^2 in an area relatively free of Galactic dust near the south ecliptic pole (SEP). The ACT bands are sensitive to radiation from the CMB, the Sunyaev-Zel'dovich (SZ) effect from galaxy clusters, and to emission by radio and dusty star-forming galaxies (DSFGs), while the dominant contribution to the BLAST bands is from DSFGs. We confirm and extend the BLAST analysis of clustering with an independent pipeline, and also detect correlations between the ACT and BLAST maps at over 25sigma significance, which we interpret as a detection of the DSFGs in the ACT maps. In addition to a Poisson component in the cross-frequency power spectra, we detect a clustered signal at >4sigma, and using a model for the DSFG evolution and number counts, we successfully fit all our spectra with a linear clustering model and a bias that depends only on redshift and not on scale. Finally, the data are compared to, and generally agree with, phenomenological models for the DSFG population. This study represents a first of its kind, and demonstrates the constraining power of the cross-frequency correlation technique to constrain models for the DSFGs. Similar analyses with more data will impose tight constraints on future models.

preprint2011arXiv

Discovery of a Multiply-Lensed Submillimeter Galaxy in Early HerMES Herschel/SPIRE Data

We report the discovery of a bright ($f(250\mum) > 400$ mJy), multiply-lensed submillimeter galaxy \obj\ in {\it Herschel}/SPIRE Science Demonstration Phase data from the HerMES project. Interferometric 880\mum\ Submillimeter Array observations resolve at least four images with a large separation of $\sim 9\arcsec$. A high-resolution adaptive optics $K_p$ image with Keck/NIRC2 clearly shows strong lensing arcs. Follow-up spectroscopy gives a redshift of $z=2.9575$, and the lensing model gives a total magnification of $μ\sim 11 \pm 1$. The large image separation allows us to study the multi-wavelength spectral energy distribution (SED) of the lensed source unobscured by the central lensing mass. The far-IR/millimeter-wave SED is well described by a modified blackbody fit with an unusually warm dust temperature, $88 \pm 3$ K. We derive a lensing-corrected total IR luminosity of $(1.43 \pm 0.09) \times 10^{13}\, \mathrm{L}_{\odot}$, implying a star formation rate of $\sim 2500\, \mathrm{M}_{\odot}\, \mathrm{yr}^{-1}$. However, models primarily developed from brighter galaxies selected at longer wavelengths are a poor fit to the full optical-to-millimeter SED. A number of other strongly lensed systems have already been discovered in early {\it Herschel} data, and many more are expected as additional data are collected.

preprint2011arXiv

HerMES: detection of cosmic magnification of sub-mm galaxies using angular cross-correlation

Cosmic magnification is due to the weak gravitational lensing of sources in the distant Universe by foreground large-scale structure leading to coherent changes in the observed number density of the background sources. Depending on the slope of the background source number counts, cosmic magnification causes a correlation between the background and foreground galaxies, which is unexpected in the absence of lensing if the two populations are spatially disjoint. Previous attempts using submillimetre (sub-mm) sources have been hampered by small number statistics. The large number of sources detected in the {\it Herschel} Multi-tiered Extra-galactic Survey (HerMES) Lockman-SWIRE field enables us to carry out the first robust study of the cross-correlation between sub-mm sources and sources at lower redshifts. Using ancillary data we compile two low-redshift samples from SDSS and SWIRE with <z> ~ 0.2 and 0.4, respectively, and cross-correlate with two sub-mm samples based on flux density and colour criteria, selecting galaxies preferentially at z ~ 2. We detect cross-correlation on angular scales between ~1 and 50 arcmin and find clear evidence that this is primarily due to cosmic magnification. A small, but non-negligible signal from intrinsic clustering is likely to be present due to the tails of the redshift distribution of the sub-mm sources overlapping with those of the foreground samples.

preprint2011arXiv

HerMES: point source catalogues from deep Herschel-SPIRE observations

We describe the generation of single-band point source catalogues from submillimetre Herschel-SPIRE observations taken as part of the Science Demonstration Phase of the Herschel Multi-tiered Extragalactic Survey (HerMES). Flux densities are found by means of peak-finding and the fitting of a Gaussian point-response function. With highly-confused images, careful checks must be made on the completeness and flux density accuracy of the detected sources. This is done by injecting artificial sources into the images and analysing the resulting catalogues. Measured flux densities at which 50 per cent of injected sources result in good detections at (250, 350, 500) μm range from (11.6, 13.2, 13.1) mJy to (25.7, 27.1, 35.8) mJy, depending on the depth of the observation (where a `good' detection is taken to be one with positional offset less than one full-width half-maximum of the point-response function, and with the measured flux density within a factor of 2 of the flux density of the injected source). This paper acts as a reference for the 2010 July HerMES public data release.

preprint2011arXiv

Herschel/HerMES: The X-ray - Infrared correlation for star-forming galaxies at z~1

For the first time, we investigate the X-ray/infrared (IR) correlation for star-forming galaxies at z~1, using SPIRE submm data from the recently-launched Herschel Space Observatory and deep X-ray data from the 2Ms Chandra deep field north (CDFN) survey. We examine the X-ray/IR correlation in the soft X-ray (SX, 0.5-2 keV) and hard X-ray (HX, 2-10 keV) bands by comparing our z~1 SPIRE-detected star-forming galaxies (SFGs) to equivalently IR-luminous (L_IR >10^10 L_sun) samples in the local/low redshift Universe. Our results suggest that the X-ray/IR properties of the SPIRE SFGs are on average similar to those of their local counterparts, as we find no evidence for evolution in the L_SX/L_IR and L_HX/L_IR ratios with redshift. We note however, that at all redshifts, both L_SX/L_IR and L_HX/L_IR are strongly dependent on IR luminosity, with luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs,L_IR >10^11 L_sun) having up to an order of magnitude lower values than normal infrared galaxies (L_IR <10^11 L_sun). We derive a L_SX-L_IR relation and confirm the applicability of an existing L_HX-L_IR relation for both local and distant LIRGs and ULIRGs, consistent with a scenario where X-ray luminosity is correlated with the star-formation rate (SFR).

preprint2011arXiv

Measuring star formation in high-z massive galaxies: A mid-infrared to submillimeter study of the GOODS NICMOS Survey sample

We present measurements of the mean mid-infrared-to-submillimeter flux densities of massive (M\ast \approx 2 \times 10^11 Msun) galaxies at redshifts 1.7 < z < 2.9, obtained by stacking positions of known objects taken from the GOODS NICMOS Survey (GNS) catalog on maps: at 24 μm (Spitzer/MIPS); 70, 100, and 160μm (Herschel/PACS); 250, 350, 500μm (BLAST); and 870μm (LABOCA). A modified blackbody spectrum fit to the stacked flux densities indicates a median [interquartile] star-formation rate of SFR = 63 [48, 81] Msun yr^-1 . We note that not properly accounting for correlations between bands when fitting stacked data can significantly bias the result. The galaxies are divided into two groups, disk-like and spheroid-like, according to their Sersic indices, n. We find evidence that most of the star formation is occurring in n \leq 2 (disk-like) galaxies, with median [interquartile] SFR = 122 [100,150] Msun yr^-1, while there are indications that the n > 2 (spheroid-like) population may be forming stars at a median [interquartile] SFR = 14 [9,20] Msun yr^-1, if at all. Finally, we show that star formation is a plausible mechanism for size evolution in this population as a whole, but find only marginal evidence that it is what drives the expansion of the spheroid-like galaxies.

preprint2011arXiv

Modeling of the HERMES J105751.1+573027 submillimeter source lensed by a dark matter dominated foreground group of galaxies

We present the results of a gravitational lensing analysis of the bright $\zs=2.957$ sub-millimeter galaxy (SMG), HERMES J105751.1+573027 found in {\it Herschel}/SPIRE Science Demonstration Phase data from the Herschel Multi-tiered Extragalactic Survey (HerMES) project. The high resolution imaging available in optical and Near-IR channels, along with CO emission obtained with the Plateau de Bure Interferometer, allow us to precisely estimate the intrinsic source extension and hence estimate the total lensing magnification to be $μ=10.9\pm 0.7$. We measure the half-light radius $R_{\rm eff}$ of the source in the rest-frame Near-UV and $V$ bands that characterize the unobscured light coming from stars and find $R_{\rm eff,*}= [2.0 \pm 0.1]$ kpc, in good agreement with recent studies on the Submillimeter Galaxy population. This lens model is also used to estimate the size of the gas distribution ($R_{\rm eff,gas}= [1.1\pm0.5]$) kpc by mapping back in the source plane the CO (J=5-4) transition line emission. The lens modeling yields a relatively large Einstein radius $R_{\rm Ein}= 4\farcs10 \pm 0\farcs02$, corresponding to a deflector velocity dispersion of [$483\pm 16] \,\kms$. This shows that HERMES J105751.1+573027 is lensed by a {\it galaxy group-size} dark matter halo at redshift $\zl\sim 0.6$. The projected dark matter contribution largely dominates the mass budget within the Einstein radius with $f_{\rm dm}(<R_{\rm Ein})\sim 80%$. This fraction reduces to $f_{\rm dm}(<R_{\rm eff,G1}\simeq 4.5\kpc)\sim 47%$ within the effective radius of the main deflecting galaxy of stellar mass $M_{\rm *,G1}=[8.5\pm 1.6] \times 10^{11}\msun$. At this smaller scale the dark matter fraction is consistent with results already found for massive lensing ellipticals at $z\sim0.2$ from the SLACS survey.

preprint2011arXiv

No evidence for anomalously low variance circles on the sky

In a recent paper, Gurzadyan & Penrose claim to have found directions on the sky centred on which are circles of anomalously low variance in the cosmic microwave background (CMB). These features are presented as evidence for a particular picture of the very early Universe. We attempted to repeat the analysis of these authors, and we can indeed confirm that such variations do exist in the temperature variance for annuli around points in the data. However, we find that this variation is entirely expected in a sky which contains the usual CMB anisotropies. In other words, properly simulated Gaussian CMB data contain just the sorts of variations claimed. Gurzadyan & Penrose have not found evidence for pre-Big Bang phenomena, but have simply re-discovered that the CMB contains structure.

preprint2011arXiv

Polarized Radio Sources: A Study of Luminosity, Redshift and Infrared Colors

The Dominion Radio Astrophysical Observatory Deep Field polarization study has been matched with the Spitzer Wide-Area Infrared Extragalactic survey of the European Large Area Infrared Space Observatory Survey North 1 field. We have used VLA observations with a total intensity rms of 87 microJy beam^-1 to match SWIRE counterparts to the radio sources. Infrared color analysis of our radio sample shows that the majority of polarized sources are elliptical galaxies with an embedded active galactic nucleus. Using available redshift catalogs, we found 429 radio sources of which 69 are polarized with redshifts in the range of 0.04 < z <3.2. We find no correlation between redshift and percentage polarization for our sample. However, for polarized radio sources, we find a weak correlation between increasing percentage polarization and decreasing luminosity.

preprint2011arXiv

Precision cosmology defeats void models for acceleration

The suggestion that we occupy a privileged position near the centre of a large, nonlinear, and nearly spherical void has recently attracted much attention as an alternative to dark energy. Putting aside the philosophical problems with this scenario, we perform the most complete and up-to-date comparison with cosmological data. We use supernovae and the full cosmic microwave background spectrum as the basis of our analysis. We also include constraints from radial baryonic acoustic oscillations, the local Hubble rate, age, big bang nucleosynthesis, the Compton y-distortion, and for the first time include the local amplitude of matter fluctuations, σ_8. These all paint a consistent picture in which voids are in severe tension with the data. In particular, void models predict a very low local Hubble rate, suffer from an "old age problem", and predict much less local structure than is observed.

preprint2011arXiv

Predicted Constraints on Cosmic String Tension from Planck and Future CMB Polarization Measurements

We perform a Fisher matrix calculation of the predicted uncertainties on estimates of the cosmic string tension Gmu from upcoming observational data (namely, cosmic microwave background power spectra from the Planck satellite and an idealized future polarization experiment). We employ simulations that are more general than others commonly used in the literature, leaving the mean velocity of strings, correlation length of the string network, and "wiggliness" (which parametrizes smaller-scale structure along the strings) as free parameters that can be observationally measured. In a new code, StringFast, we implement a method for efficient computation of the C_l spectra induced by a network of strings, which is fast enough to be used in Markov Chain Monte Carlo analyses of future data. Performing a calculation with the string parameters left free results in projected constraints on Gmu that are larger than those obtained by fixing their values a priori, typically by a factor of ~2-7. We also find that if Gmu is equal to the current observational maximum, Planck will be able to make a confident detection of strings. However, if Gmu is two orders of magnitude smaller, even a perfect, lensing-free measurement of polarization power spectra will not be able to detect a nonzero string tension at better than 2 sigma confidence.

preprint2011arXiv

Sub-millimetre galaxies reside in dark matter halos with masses greater than 3x10^11 solar masses

The extragalactic background light at far-infrared wavelengths originates from optically-faint, dusty, star-forming galaxies in the universe with star-formation rates at the level of a few hundred solar masses per year. Due to the relatively poor spatial resolution of far-infrared telescopes, the faint sub-millimetre galaxies are challenging to study individually. Instead, their average properties can be studied using statistics such as the angular power spectrum of the background intensity variations. A previous attempt at measuring this power spectrum resulted in the suggestion that the clustering amplitude is below the level computed with a simple ansatz based on a halo model. Here we report a clear detection of the excess clustering over the linear prediction at arcminute angular scales in the power spectrum of brightness fluctuations at 250, 350, and 500 microns. From this excess, we find that sub-millimetre galaxies are located in dark matter halos with a minimum mass of log[M_min/M_sun ]= 11.5^+0.7_-0.2 at 350 microns. This minimum dark matter halo mass corresponds to the most efficient mass scale for star formation in the universe, and is lower than that predicted by semi-analytical models for galaxy formation.

preprint2011arXiv

The BLAST View of the Star Forming Region in Aquila (ell=45deg,b=0deg)

We have carried out the first general submillimeter analysis of the field towards GRSMC 45.46+0.05, a massive star forming region in Aquila. The deconvolved 6 deg^2 (3\degree X 2\degree) maps provided by BLAST in 2005 at 250, 350, and 500 micron were used to perform a preliminary characterization of the clump population previously investigated in the infrared, radio, and molecular maps. Interferometric CORNISH data at 4.8 GHz have also been used to characterize the Ultracompact HII regions (UCHIIRs) within the main clumps. By means of the BLAST maps we have produced an initial census of the submillimeter structures that will be observed by Herschel, several of which are known Infrared Dark Clouds (IRDCs). Our spectral energy distributions of the main clumps in the field, located at ~7 kpc, reveal an active population with temperatures of T~35-40 K and masses of ~10^3 Msun for a dust emissivity index beta=1.5. The clump evolutionary stages range from evolved sources, with extended HII regions and prominent IR stellar population, to massive young stellar objects, prior to the formation of an UCHIIR.The CORNISH data have revealed the details of the stellar content and structure of the UCHIIRs. In most cases, the ionizing stars corresponding to the brightest radio detections are capable of accounting for the clump bolometric luminosity, in most cases powered by embedded OB stellar clusters.

preprint2011arXiv

The Herschel Multi-Tiered Extragalactic Survey: SPIRE-mm Photometric Redshifts

We investigate the potential of submm-mm and submm-mm-radio photometric redshifts using a sample of mm-selected sources as seen at 250, 350 and 500 μm by the SPIRE instrument on Herschel. From a sample of 63 previously identified mm-sources with reliable radio identifications in the GOODS-N and Lockman Hole North fields 46 (73 per cent) are found to have detections in at least one SPIRE band. We explore the observed submm/mm colour evolution with redshift, finding that the colours of mm-sources are adequately described by a modified blackbody with constant optical depth τ = (ν/ν0)^β where β = +1.8 and ν0 = c/100 μm. We find a tight correlation between dust temperature and IR luminosity. Using a single model of the dust temperature and IR luminosity relation we derive photometric redshift estimates for the 46 SPIRE detected mm-sources. Testing against the 22 sources with known spectroscopic, or good quality optical/near-IR photometric, redshifts we find submm/mm photometric redshifts offer a redshift accuracy of |z|/(1+z) = 0.16 (< |z| >= 0.51). Including constraints from the radio-far IR correlation the accuracy is improved to |z|/(1 + z) = 0.15 (< |z| >= 0.45). We estimate the redshift distribution of mm-selected sources finding a significant excess at z > 3 when compared to ~ 850 μm selected samples.

preprint2010arXiv

A joint analysis of BLAST 250--500um and LABOCA 870um observations in the Extended Chandra Deep Field South

We present a joint analysis of the overlapping BLAST 250, 350, 500um, and LABOCA 870um observations (from the LESS survey) of the ECDF-S. Out to z~3, the BLAST filters sample near the peak wavelength of far-infrared (FIR) emission from galaxies (rest-frame wavelengths ~60--200um), primarily produced by dust heated through absorption in star-forming clouds. However, identifying counterparts to individual BLAST peaks is very challenging, given the large beams (FWHM 36--60"). In contrast, the 870um observations have a significantly smaller 19" FWHM beam, and are sensitive to higher redshifts (z~1--5, and potentially beyond) due to the more favourable negative K-correction. We use the LESS data, as well as deep Spitzer and VLA imaging, to identify 118 individual sources that produce significant emission in the BLAST bands. We characterize the temperatures and FIR luminosities for a subset of 69 sources which have well-measured submm SEDs and redshift measurements out to z~3. For flux-limited sub-samples in each BLAST band, and a dust emissivity index β=2.0, we find a median temperature T=30K (all bands) as well as median redshifts: z=1.1 (interquartile range 0.2--1.9) for S250 > 40mJy; z=1.3 (interquartile range 0.6--2.1) for S350 > 30mJy; and z=1.6 (interquartile range 1.3--2.3) for S500 > 20mJy. Taking into account selection effects for our survey (a bias toward detecting lower-temperature galaxies), we find no evidence for evolution in the local FIR-temperature correlation out to z~2.5. Comparing with star-forming galaxy SED templates, about 8% of our sample appears to exhibit significant excesses in the radio and/or mid-IR, consistent with those sources harbouring an AGN. We describe the following techniques in two appendices: our `matched filter' for identifying sources in the presence of point-source confusion; and our approach for identifying counterparts using likelihood ratios.

preprint2010arXiv

A panchromatic study of BLAST counterparts: total star-formation rate, morphology, AGN fraction and stellar mass

We carry out a multi-wavelength study of individual galaxies detected by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) and identified at other wavelengths, using data spanning the radio to the ultraviolet (UV). We develop a Monte Carlo method to account for flux boosting, source blending, and correlations among bands, which we use to derive deboosted far-infrared (FIR) luminosities for our sample. We estimate total star-formation rates for BLAST counterparts with z < 0.9 by combining their FIR and UV luminosities. Star formation is heavily obscured at L_FIR > 10^11 L_sun, z > 0.5, but the contribution from unobscured starlight cannot be neglected at L_FIR < 10^11 L_sun, z < 0.25. We assess that about 20% of the galaxies in our sample show indication of a type-1 active galactic nucleus (AGN), but their submillimeter emission is mainly due to star formation in the host galaxy. We compute stellar masses for a subset of 92 BLAST counterparts; these are relatively massive objects, with a median mass of ~10^11 M_sun, which seem to link the 24um and SCUBA populations, in terms of both stellar mass and star-formation activity. The bulk of the BLAST counterparts at z<1 appear to be run-of-the-mill star-forming galaxies, typically spiral in shape, with intermediate stellar masses and practically constant specific star-formation rates. On the other hand, the high-z tail of the BLAST counterparts significantly overlaps with the SCUBA population, in terms of both star-formation rates and stellar masses, with observed trends of specific star-formation rate that support strong evolution and downsizing.

preprint2010arXiv

A Pilot Study for the SCUBA-2 'All-Sky' Survey

We have carried out a pilot study for the SCUBA-2 'All-Sky' Survey, SASSy, a wide and shallow mapping project at 850 microns, designed to find rare objects, both Galactic and extragalactic. Two distinct sets of exploratory observations were undertaken, and used to test the SASSy approach and data reduction pipeline. The first was a 0.5 by 0.5 degrees map around the nearby galaxy NGC 2559. The galaxy was easily detected at 156 mJy, but no other convincing sources are present in the map. Comparison with other galaxies with similar wavelength coverage indicates that NGC 2559 has relatively warm dust. The second observations cover 1 square degree around the W5-E HII region. As well as diffuse structure in the map, a filtering approach was able to extract 27 compact sources with signal-to-noise greater than 6. By matching with data at other wavelengths we can see that the SCUBA-2 data can be used to discriminate the colder cores. Together these observations show that the SASSy project will be able to meet its original goals of detecting new bright sources which will be ideal for follow-up observations with other facilities.

preprint2010arXiv

Deconvolution of Images from BLAST 2005: Insight into the K3-50 and IC 5146 Star-Forming Regions

We present an implementation of the iterative flux-conserving Lucy-Richardson (L-R) deconvolution method of image restoration for maps produced by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST). We have analyzed its performance and convergence extensively through simulations and cross-correlations of the deconvolved images with available highresolution maps. We present new science results from two BLAST surveys, in the Galactic regions K3-50 and IC 5146, further demonstrating the benefits of performing this deconvolution. We have resolved three clumps within a radius of 4.'5 inside the star-forming molecular cloud containing K3-50. Combining the well-resolved dust emission map with available multi-wavelength data, we have constrained the Spectral Energy Distributions (SEDs) of five clumps to obtain masses (M), bolometric luminosities (L), and dust temperatures (T). The L-M diagram has been used as a diagnostic tool to estimate the evolutionary stages of the clumps. There are close relationships between dust continuum emission and both 21-cm radio continuum and 12CO molecular line emission. The restored extended large scale structures in the Northern Streamer of IC 5146 have a strong spatial correlation with both SCUBA and high resolution extinction images. A dust temperature of 12 K has been obtained for the central filament. We report physical properties of ten compact sources, including six associated protostars, by fitting SEDs to multi-wavelength data. All of these compact sources are still quite cold (typical temperature below ~ 16 K) and are above the critical Bonner-Ebert mass. They have associated low-power Young Stellar Objects (YSOs). Further evidence for starless clumps has also been found in the IC 5146 region.

preprint2010arXiv

Evolution of Dust Temperature of Galaxies through Cosmic Time as seen by Herschel

We study the dust properties of galaxies in the redshift range 0.1<z<2.8 observed by the Herschel Space Observatory in the field of the Great Observatories Origins Deep Survey-North as part of PEP and HerMES key programmes. Infrared (IR) luminosity (L_IR) and dust temperature (T_dust) of galaxies are derived from the spectral energy distribution (SED) fit of the far-infrared (FIR) flux densities obtained with PACS and SPIRE instruments onboard Herschel. As a reference sample, we also obtain IR luminosities and dust temperatures of local galaxies at z<0.1 using AKARI and IRAS data in the field of the Sloan Digital Sky Survey. We compare the L_IR-T_dust relation between the two samples and find that: the median T_dust of Herschel-selected galaxies at z>0.5 with L_IR>5x10^{10} L_\odot, appears to be 2-5 K colder than that of AKARI-selected local galaxies with similar luminosities; and the dispersion in T_dust for high-z galaxies increases with L_IR due to the existence of cold galaxies that are not seen among local galaxies. We show that this large dispersion of the L_IR-T_dust relation can bridge the gap between local star-forming galaxies and high-z submillimeter galaxies (SMGs). We also find that three SMGs with very low T_dust (<20 K) covered in this study have close neighbouring sources with similar 24-\mum brightness, which could lead to an overestimation of FIR/(sub)millimeter fluxes of the SMGs.

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

HerMES: Deep Galaxy Number Counts from a P(D) Fluctuation Analysis of SPIRE Science Demonstration Phase Observations

Dusty, star forming galaxies contribute to a bright, currently unresolved cosmic far-infrared background. Deep Herschel-SPIRE images designed to detect and characterize the galaxies that comprise this background are highly confused, such that the bulk lies below the classical confusion limit. We analyze three fields from the HerMES programme in all three SPIRE bands (250, 350, and 500 microns); parameterized galaxy number count models are derived to a depth of ~2 mJy/beam, approximately 4 times the depth of previous analyses at these wavelengths, using a P(D) (probability of deflection) approach for comparison to theoretical number count models. Our fits account for 64, 60, and 43 per cent of the far-infrared background in the three bands. The number counts are consistent with those based on individually detected SPIRE sources, but generally inconsistent with most galaxy number counts models, which generically overpredict the number of bright galaxies and are not as steep as the P(D)-derived number counts. Clear evidence is found for a break in the slope of the differential number counts at low flux densities. Systematic effects in the P(D) analysis are explored. We find that the effects of clustering have a small impact on the data, and the largest identified systematic error arises from uncertainties in the SPIRE beam.

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

HerMES: Halo Occupation Number and Bias Properties of Dusty Galaxies from Angular Clustering Measurements

We measure the angular correlation function, w(theta), from 0.5 to 30 arcminutes of detected sources in two wide fields of the Herschel Multi-tiered Extragalactic Survey (HerMES). Our measurements are consistent with the expected clustering shape from a population of sources that trace the dark matter density field, including non-linear clustering at arcminute angular scales arising from multiple sources that occupy the same dark matter halos. By making use of the halo model to connect the spatial clustering of sources to the dark matter halo distribution, we estimate source bias and halo occupation number for dusty sub-mm galaxies at z ~ 2. We find that sub-mm galaxies with 250 micron flux densities above 30 mJy reside in dark matter halos with mass above (5\pm4) x 10^12 M_sun, while (14\pm8)% of such sources appear as satellites in more massive halos.

preprint2010arXiv

HerMES: SPIRE galaxy number counts at 250, 350 and 500 microns

Emission at far-infrared wavelengths makes up a significant fraction of the total light detected from galaxies over the age of Universe. Herschel provides an opportunity for studying galaxies at the peak wavelength of their emission. Our aim is to provide a benchmark for models of galaxy population evolution and to test pre-existing models of galaxies. With the Herschel Multi-tiered Extra-galactic survey, HerMES, we have observed a number of fields of different areas and sensitivity using the SPIRE instrument on Herschel. We have determined the number counts of galaxies down to ~20 mJy. Our constraints from directly counting galaxies are consistent with, though more precise than, estimates from the BLAST fluctuation analysis. We have found a steep rise in the Euclidean normalised counts at <100 mJy. We have directly resolved 15% of the infrared extra-galactic background at the wavelength near where it peaks.

preprint2010arXiv

HerMES: SPIRE Science Demonstration Phase Maps

We describe the production and verification of sky maps of the five SPIRE fields observed as part of the Herschel Multi-tiered Extragalactic Survey (HerMES) during the Science Demonstration Phase (SDP) of the Herschel mission. We have implemented an iterative map-making algorithm (SHIM; The SPIRE-HerMES Iterative Mapper) to produce high fidelity maps that preserve extended diffuse emission on the sky while exploiting the repeated observations of the same region of the sky with many detectors in multiple scan directions to minimize residual instrument noise. We specify here the SHIM algorithm and outline the various tests that were performed to determine and characterize the quality of the maps and verify that the astrometry, point source flux and power on all relevant angular scales meets the needs of the HerMES science goals. These include multiple jackknife tests, determination of the map transfer function and detailed examination of the power spectra of both sky and jackknife maps. The map transfer function is approximately unity on scales from one arcminute to one degree. Final maps (v1.0), including multiple jackknives, as well as the SHIM pipeline, have been used by the HerMES team for the production of SDP papers.

preprint2010arXiv

HerMES: The SPIRE confusion limit

We report on the sensitivity of SPIRE photometers on the Herschel Space Observatory. Specifically, we measure the confusion noise from observations taken during the Science Demonstration Phase of the Herschel Multi-tiered Extragalactic Survey. Confusion noise is defined to be the spatial variation of the sky intensity in the limit of infinite integration time, and is found to be consistent among the different fields in our survey at the level of 5.8, 6.3 and 6.8 mJy/beam at 250, 350 and 500 microns, respectively. These results, together with the measured instrument noise, may be used to estimate the integration time required for confusion-limited maps, and provide a noise estimate for maps obtained by SPIRE.

preprint2010arXiv

HerMES: The Submillimeter Spectral Energy Distributions of Herschel/SPIRE-Detected Galaxies

We present colours of sources detected with the Herschel/SPIRE instrument in deep extragalactic surveys of the Lockman Hole, Spitzer-FLS, and GOODS-N fields in three photometric bands at 250, 350 and 500 micrometers. We compare these with expectations from the literature and discuss associated uncertainties and biases in the SPIRE data. We identify a 500 micrometer flux limited selection of sources from the HerMES point source catalogue that appears free from neighbouring/blended sources in all three SPIRE bands. We compare the colours with redshift tracks of various contemporary models. Based on these spectral templates we show that regions corresponding to specific population types and redshifts can be identified better in colour-flux space. The redshift tracks as well as the colour-flux plots imply a majority of detected objects with redshifts at 1<z<3.5, somewhat depending on the group of model SEDs used. We also find that a population of S_250/S_350<0.8 at fluxes above 50 mJy as observed by SPIRE is not well represented by contemporary models and could consist of a mix of cold and lensed galaxies.

preprint2010arXiv

Herschel unveils a puzzling uniformity of distant dusty galaxies

The Herschel Space Observatory enables us to accurately measure the bolometric output of starburst galaxies and active galactic nuclei (AGN) by directly sampling the peak of their far-infrared (IR) emission. Here we examine whether the spectral energy distribution (SED) and dust temperature of galaxies have strongly evolved since z~2.5. We use Herschel deep extragalactic surveys from 100 to 500um to compute total IR luminosities in galaxies down to the faintest levels, using PACS and SPIRE in the GOODS-North field (PEP and HerMES key programs). We show that measurements in the SPIRE bands can be used below the statistical confusion limit if information at higher spatial resolution is used to identify isolated galaxies whose flux is not boosted by bright neighbors. Below z~1.5, mid-IR extrapolations are correct for star-forming galaxies with a dispersion of only 40% (0.15dex), therefore similar to z~0 galaxies. This narrow distribution is puzzling when considering the range of physical processes that could have affected the SED of these galaxies. Extrapolations from only one of the 160um, 250um or 350um bands alone tend to overestimate the total IR luminosity. This may be explained by the lack of far-IR constraints around and above ~150um (rest-frame) on local templates. We also note that the dust temperature of luminous IR galaxies around z~1 is mildly colder by 10-15% than their local analogs and up to 20% for ULIRGs at z~1.6. Above z=1.5, distant galaxies are found to exhibit a substantially larger mid- over far-IR ratio, which could either result from stronger broad emission lines or warm dust continuum heated by a hidden AGN. Two thirds of the AGNs identified in the field with a measured redshift exhibit the same behavior as purely star-forming galaxies. Hence a large fraction of AGNs harbor star formation at very high SFR and in conditions similar to purely star-forming galaxies.

preprint2010arXiv

Herschel-SPIRE, Far-Infrared Properties of Millimetre-Bright and -Faint Radio Galaxies

(abridged) We present the first study of the farIR properties of high redshift, radio-selected ULIRGs using deep observations obtained with SPIRE from the Herschel Multi-tiered Extragalactic Survey (HerMES). These galaxies span a large range of 850um fluxes from submillimetre-luminous ~10mJy "SCUBA galaxies -- SMGs" to ~1.5mJy from stacked SCUBA non-detections, thus likely representing a complete distribution of ULIRG spectral energy distributions. From Keck spectroscopic surveys in the Lockman-North field we identified a sample of 31 SMGs and 37 submillimetre-faint, optically-faint radio galaxies (OFRGs), all with radio-inferred IR luminosities >10^12 Lsun. These galaxies were cross-identified with SPIRE 250, 350 and 500um catalogs based on fluxes extracted at 24um positions in the SWIRE survey, yielding a sample of more than half of the galaxies well detected in at least two of the SPIRE bandpasses. By fitting greybody dust models to the SPIRE photometry together with SCUBA 850um measurements, we infer dust temperatures and far-infrared luminosities. The OFRGs detected by SPIRE have median <T_d>= 41+-5 K and the SMGs have <T_d>= 34+-5 K, both in reasonable agreement with previous (pre-Herschel) estimates, reaffirming that the local FIR/radio correlation holds (at least for this subset of high-z ULIRGs) at high redshift. Our observations firstly confirm that a substantial fraction of OFRGs exhibit large infrared luminosities corresponding to SFRs of ~400 Msun/yr. The SPIRE observations secondly confirm the higher dust temperatures for these OFRGs than similarly selected SMGs, consistent with early predictions of the submm-faint radio populations. Our observations also clearly confirm the large infrared luminosities of most SMGs selected with S850um>5 mJy and radio and strong 24um detections, corresponding to SFRs of ~700 Msun/yr.

preprint2010arXiv

Induced CMB quadrupole from pointing offsets

Recent claims in the literature have suggested that the {\it WMAP} quadrupole is not primordial in origin, and arises from an aliasing of the much larger dipole field because of incorrect satellite pointing. We attempt to reproduce this result and delineate the key physics leading to the effect. We find that, even if real, the induced quadrupole would be smaller than claimed. We discuss reasons why the {\it WMAP} data are unlikely to suffer from this particular systematic effect, including the implications for observations of point sources. Given this evidence against the reality of the effect, the similarity between the pointing-offset-induced signal and the actual quadrupole then appears to be quite puzzling. However, we find that the effect arises from a convolution between the gradient of the dipole field and anisotropic coverage of the scan direction at each pixel. There is something of a directional conspiracy here -- the dipole signal lies close to the Ecliptic Plane, and its direction, together with the {\it WMAP} scan strategy, results in a strong coupling to the $Y_{2,\,-1}$ component in Ecliptic co-ordinates. The dominant strength of this component in the measured quadrupole suggests that one should exercise increased caution in interpreting its estimated amplitude. The {\it Planck} satellite has a different scan strategy which does not so directly couple the dipole and quadrupole in this way and will soon provide an independent measurement.

preprint2010arXiv

Let's talk about varying G

It is possible that fundamental constants may not be constant at all. There is a generally accepted view that one can only talk about variations of dimensionless quantities, such as the fine structure constant $α_{\rm e}\equiv e^2/4πε_0\hbar c$. However, constraints on the strength of gravity tend to focus on G itself, which is problematic. We stress that G needs to be multiplied by the square of a mass, and hence, for example, one should be constraining $α_{\rm g}\equiv G m_{\rm p}^2/\hbar c$, where $m_{\rm p}$ is the proton mass. Failure to focus on such dimensionless quantities makes it difficult to interpret the physical dependence of constraints on the variation of G in many published studies. A thought experiment involving talking to observers in another universe about the values of physical constants may be useful for distinguishing what is genuinely measurable from what is merely part of our particular system of units.

preprint2010arXiv

Origins of the extragalactic background at 1mm from a combined analysis of the AzTEC and MAMBO data in GOODS-N

We present a study of the cosmic infrared background, which is a measure of the dust obscured activity in all galaxies in the Universe. We venture to isolate the galaxies responsible for the background at 1mm; with spectroscopic and photometric redshifts we constrain the redshift distribution of these galaxies. We create a deep 1.16mm map (sigma ~ 0.5mJy) by combining the AzTEC 1.1mm and MAMBO 1.2mm datasets in GOODS-N. This combined map contains 41 secure detections, 13 of which are new. By averaging the 1.16mm flux densities of individually undetected galaxies with 24um flux densities > 25uJy, we resolve 31--45 per cent of the 1.16mm background. Repeating our analysis on the SCUBA 850um map, we resolve a higher percentage (40--64 per cent) of the 850um background. A majority of the background resolved (attributed to individual galaxies) at both wavelengths comes from galaxies at z > 1.3. If the ratio of the resolved submillimeter to millimeter background is applied to a reasonable scenario for the origins of the unresolved submillimeter background, 60--88 per cent of the total 1.16mm background comes from galaxies at z > 1.3.

preprint2010arXiv

SCUBA-2 Data Processing

SCUBA-2 is the largest submillimetre array camera in the world and was commissioned on the James Clerk Maxwell Telescope (JCMT) with two arrays towards the end of 2009. A period of shared-risks observing was then completed and the full planned complement of 8 arrays, 4 at 850 microns and 4 at 450 microns, are now installed and ready to be commissioned. SCUBA-2 has 10,240 bolometers, corresponding to a data rate of 8 MB/s when sampled at the nominal rate of 200 Hz. The pipeline produces useful maps in near real time at the telescope and often publication quality maps in the JCMT Science Archive (JSA) hosted at the Canadian Astronomy Data Centre (CADC).

preprint2010arXiv

Spatial and temporal tuning in void models for acceleration

There has been considerable interest in recent years in cosmological models in which we inhabit a very large, underdense void as an alternative to dark energy. A longstanding objection to this proposal is that observations limit our position to be very close to the void centre. By selecting from a family of void profiles that fit supernova luminosity data, we carefully determine how far from the centre we could be. To do so, we use the observed dipole component of the cosmic microwave background, as well as an additional stochastic peculiar velocity arising from primordial perturbations. We find that we are constrained to live within 80 Mpc of the centre of a void--a somewhat weaker constraint than found in previous studies, but nevertheless a strong violation of the Copernican principle. By considering how such a Gpc-scale void would appear on the microwave sky, we also show that there can be a maximum of one of these voids within our Hubble radius. Hence, the constraint on our position corresponds to a fraction of the Hubble volume of order 10^{-8}. Finally, we use the fact that void models only look temporarily similar to a cosmological-constant-dominated universe to argue that these models are not free of temporal fine-tuning.

preprint2010arXiv

Spitzer MIPS 24 and 70 micron Imaging near the South Ecliptic Pole: Maps and Source Catalogs

We have imaged an 11.5 sq. deg. region of sky towards the South Ecliptic Pole (RA = 04h43m, Dec = -53d40m, J2000) at 24 and 70 microns with MIPS, the Multiband Imaging Photometer for Spitzer. This region is coincident with a field mapped at longer wavelengths by AKARI and the Balloon-borne Large Aperture Submillimeter Telescope. We discuss our data reduction and source extraction procedures. The median depths of the maps are 47 microJy/beam at 24 micron and 4.3 mJy/beam at 70 micron. At 24 micron, we identify 93098 point sources with signal-to-noise ratio (SNR) >5, and an additional 63 resolved galaxies; at 70 micron, we identify 891 point sources with SNR >6. From simulations, we determine a false detection rate of 1.8% (1.1%) for the 24 micron (70 micron) catalog. The 24 and 70 micron point-source catalogs are 80% complete at 230 microJy and 11 mJy, respectively. These mosaic images and source catalogs will be available to the public through the NASA/IPAC Infrared Science Archive.

preprint2010arXiv

Submillimetre observations of galaxy clusters with BLAST: the star-formation activity in Abell 3112

We present observations at 250, 350, and 500 um of the nearby galaxy cluster Abell 3112 (z=0.075) carried out with BLAST, the Balloon-borne Large Aperture Submillimeter Telescope. Five cluster members are individually detected as bright submillimetre sources. Their far-infrared SEDs and optical colours identify them as normal star-forming galaxies of high mass, with globally evolved stellar populations. They all have B-R colours of 1.38+/-0.08, transitional between the blue, active population and the red, evolved galaxies that dominate the cluster core. We stack to determine the mean submillimetre emission from all cluster members, which is determined to be 16.6+/-2.5, 6.1+/-1.9, and 1.5+/-1.3 mJy at 250, 350, and 500 um, respectively. Stacking analyses of the submillimetre emission of cluster members reveal trends in the mean far-infrared luminosity with respect to cluster-centric radius and Ks-band magnitude. We find that a large fraction of submillimetre emission comes from the boundary of the inner, virialized region of the cluster, at cluster-centric distances around R_500. Stacking also shows that the bulk of the submillimetre emission arises in intermediate-mass galaxies (L<L*), with Ks magnitude ~1 mag fainter than the giant ellipticals. The results and constraints obtained in this work will provide a useful reference for the forthcoming surveys to be conducted on galaxy clusters by Herschel.

preprint2010arXiv

The Balloon-Borne Large Aperture Submillimeter Telescope (BLAST) 2005: A 10 deg^2 Survey of Star Formation in Cygnus X

We present Cygnus X in a new multi-wavelength perspective based on an unbiased BLAST survey at 250, 350, and 500 micron, combined with rich datasets for this well-studied region. Our primary goal is to investigate the early stages of high mass star formation. We have detected 184 compact sources in various stages of evolution across all three BLAST bands. From their well-constrained spectral energy distributions, we obtain the physical properties mass, surface density, bolometric luminosity, and dust temperature. Some of the bright sources reaching 40 K contain well-known compact H II regions. We relate these to other sources at earlier stages of evolution via the energetics as deduced from their position in the luminosity-mass (L-M) diagram. The BLAST spectral coverage, near the peak of the spectral energy distribution of the dust, reveals fainter sources too cool (~ 10 K) to be seen by earlier shorter-wavelength surveys like IRAS. We detect thermal emission from infrared dark clouds and investigate the phenomenon of cold ``starless cores" more generally. Spitzer images of these cold sources often show stellar nurseries, but these potential sites for massive star formation are ``starless" in the sense that to date there is no massive protostar in a vigorous accretion phase. We discuss evolution in the context of the L-M diagram. Theory raises some interesting possibilities: some cold massive compact sources might never form a cluster containing massive stars; and clusters with massive stars might not have an identifiable compact cold massive precursor.

preprint2010arXiv

The DRAO Planck Deep Fields: the polarization properties of radio galaxies at 1.4 GHz

We present results of deep polarization imaging at 1.4 GHz with the Dominion Radio Astrophysical Observatory as part of the DRAO Planck Deep Fields project. This deep extragalactic field covers 15.16 square degrees centered at RA = 16h 14m and DEC = 54d 56', has an angular resolution of 42" x 62" at the field center, and reaches a sensitivity of 55 microJy/beam in Stokes I and 45 microJy/beam in Stokes Q and U. We detect 958 radio sources in Stokes I of which 136 are detected in polarization. We present the Euclidean-normalized polarized differential source counts down to 400 microJy. These counts indicate that sources have a higher degree of fractional polarization at fainter Stokes I flux density levels than for brighter sources, confirming an earlier result. We find that the majority of our polarized sources are steep-spectrum objects with a mean spectral index of -0.77, and there is no correlation between fractional polarization and spectral index. We also matched deep field sources to counterparts in the Faint Images of the Radio Sky at Twenty Centimeters catalogue. Of the polarized sources, 77% show structure at the arc-second scale whereas only 38% of the sources with no detectable polarization show such structure. The median fractional polarization is for resolved sources is 6.8%, while it is 4.4% for compact objects. The polarized radio sources in our deep field are predominantly those sources which are resolved and show the highest degrees of fractional polarization, indicating that the lobe dominated structure may be the source of the highly polarized sources. These resolved radio galaxies dominate the polarized source counts at P_0 = sqrt(Q^2 + U^2) < 3 mJy.

preprint2010arXiv

The far-infrared/radio correlation as probed by Herschel

We set out to determine the ratio, q(IR), of rest-frame 8-1000um flux, S(IR), to monochromatic radio flux, S(1.4GHz), for galaxies selected at far-IR and radio wavelengths, to search for signs that the ratio evolves with redshift, luminosity or dust temperature, and to identify any far-IR-bright outliers - useful laboratories for exploring why the far-IR/radio correlation is generally so tight when the prevailing theory suggests variations are almost inevitable. We use flux-limited 250-um and 1.4-GHz samples, obtained in GOODS-N using Herschel (HerMES; PEP) and the VLA. We determine bolometric IR output using ten bands spanning 24-1250um, exploiting data from PACS and SPIRE, as well as Spitzer, SCUBA, AzTEC and MAMBO. We also explore the properties of an L(IR)-matched sample, designed to reveal evolution of q(IR) with z, spanning log L(IR) = 11-12 L(sun) and z=0-2, by stacking into the radio and far-IR images. For 1.4-GHz-selected galaxies, we see tentative evidence of a break in the flux ratio, q(IR), at L(1.4GHz) ~ 10^22.7 W/Hz, where AGN are starting to dominate the radio power density, and of weaker correlations with z and T(d). From our 250-um-selected sample we identify a small number of far-IR-bright outliers, and see trends of q(IR) with L(1.4GHz), L(IR), T(d) and z, noting that some of these are inter-related. For our L(IR)-matched sample, there is no evidence that q(IR) changes significantly as we move back into the epoch of galaxy formation: we find q(IR) goes as (1+z)^gamma, where gamma = -0.04 +/- 0.03 at z=0-2; however, discounting the least reliable data at z<0.5 we find gamma = -0.26 +/- 0.07, modest evolution which may be related to the radio background seen by ARCADE2, perhaps driven by <10uJy radio activity amongst ordinary star-forming galaxies at z>1.

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: 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-SPIRE Legacy Survey (HSLS): the scientific goals of a shallow and wide submillimeter imaging survey with SPIRE

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

preprint2010arXiv

The Hubble Space Telescope GOODS NICMOS Survey: Overview and the Evolution of Massive Galaxies at 1.5 < z < 3

We present the details and early results from a deep near-infrared survey utilising the NICMOS instrument on the Hubble Space Telescope centred around massive M_* > 10^11 M_0 galaxies at 1.7 < z < 2.9 found within the Great Observatories Origins Deep Survey (GOODS) fields. The GOODS NICMOS Survey (GNS) was designed to obtain deep F160W (H-band) imaging of 80 of these massive galaxies, as well as other colour selected objects such as Lyman-break drop-outs, BzK objects, Distant Red Galaxies, EROs, Spitzer Selected EROs, BX/BM galaxies, as well as sub-mm galaxies. We present in this paper details of the observations, our sample selection, as well as a description of features of the massive galaxies found within our survey fields. This includes: photometric redshifts, rest-frame colours, and stellar masses. We furthermore provide an analysis of the selection methods for finding massive galaxies at high redshifts, including colour selection, and how galaxy populations selected through different methods overlap. We find that a single colour selection method cannot locate all of the massive galaxies, with no one method finding more than 70 percent. We however find that the combination of these colour methods finds nearly all the massive galaxies, as selected by photometric redshifts with the exception of apparently rare blue massive galaxies. By investigating the rest-frame (U-B) vs. M_B diagram for these galaxies we furthermore show that there exists a bimodality in colour-magnitude space at z < 2, driven by stellar mass, such that the most massive galaxies are systematically red up to z~2.5, while lower mass galaxies tend to be blue. We also discuss the number densities for galaxies with stellar masses M_* > 10^11 M_0, whereby we find an increase of a factor of eight between z = 3 and z = 1.5, demonstrating that this is an epoch when massive galaxies establish most of their mass.

preprint2010arXiv

The Submillimetre Universe

Submillimetre continuum radiation allows us to probe cold objects, particularly the earliest, dusty phases of star formation, high-redshift galaxies and circumstellar disks. The submillimetre window gives a unique view of the physical and dynamical conditions in the neutral and molecular interstellar medium. In the next decade a combination of wide-field surveys with single-dish telescopes and targeted follow-up with ALMA and other facilities should enable rapid progress in answering questions about the origins of planetary systems, stars and galaxies.

preprint2010arXiv

Tilted Physics: A Cosmologically Dipole-Modulated Sky

Physical constants and cosmological parameters could vary with position. On the largest scales such variations would manifest themselves as gradients across our Hubble volume, leading to dipole-modulation of the cosmic microwave anisotropies. This generically leads to a correlation between adjacent multipoles in the spherical harmonics expansion of the sky, a distinctive signal which should be searched for in future data sets.

preprint2009arXiv

BLAST: the far-infrared/radio correlation in distant galaxies

We investigate the correlation between far-infrared (FIR) and radio luminosities in distant galaxies, a lynchpin of modern astronomy. We use data from the Balloon-borne Large Aperture Submillimetre Telescope (BLAST), Spitzer, the Large Apex BOlometer CamerA (LABOCA), the Very Large Array (VLA) and the Giant Metre-wave Radio Telescope (GMRT) in the Extended Chandra Deep Field South (ECDFS). For a catalogue of BLAST 250-micron-selected galaxies, we re-measure the 70--870-micron flux densities at the positions of their most likely 24-micron counterparts, which have a median [interquartile] redshift of 0.74 [0.25, 1.57]. From these, we determine the monochromatic flux density ratio, q_250 = log_10 (S_250micron / S_1400MHz), and the bolometric equivalent, q_IR. At z ~= 0.6, where our 250-micron filter probes rest-frame 160-micron emission, we find no evolution relative to q_160 for local galaxies. We also stack the FIR and submm images at the positions of 24-micron- and radio-selected galaxies. The difference between q_IR seen for 250-micron- and radio-selected galaxies suggests star formation provides most of the IR luminosity in ~< 100-uJy radio galaxies, but rather less for those in the mJy regime. For the 24-micron sample, the radio spectral index is constant across 0 < z < 3, but q_IR exhibits tentative evidence of a steady decline such that q_IR is proportional to (1+z)^(-0.15 +/- 0.03) - significant evolution, spanning the epoch of galaxy formation, with major implications for techniques that rely on the FIR/radio correlation. We compare with model predictions and speculate that we may be seeing the increase in radio activity that gives rise to the radio background.

preprint2009arXiv

BLAST05: Power Spectra of Bright Galactic Cirrus at Submillimeter Wavelengths

We report multi-wavelength power spectra of diffuse Galactic dust emission from BLAST observations at 250, 350, and 500 microns in Galactic Plane fields in Cygnus X and Aquila. These submillimeter power spectra statistically quantify the self-similar structure observable over a broad range of scales and can be used to assess the cirrus noise which limits the detection of faint point sources. The advent of submillimeter surveys with the Herschel Space Observatory makes the wavelength dependence a matter of interest. We show that the observed relative amplitudes of the power spectra can be related through a spectral energy distribution (SED). Fitting a simple modified black body to this SED, we find the dust temperature in Cygnus X to be 19.9 +/- 1.3 K and in the Aquila region 16.9 +/- 0.7 K. Our empirical estimates provide important new insight into the substantial cirrus noise that will be encountered in forthcoming observations.

preprint2009arXiv

Gravity heats the Universe

Structure in the Universe grew through gravitational instability from very smooth initial conditions. Energy conservation requires that the growing negative potential energy of these structures is balanced by an increase in kinetic energy. A fraction of this is converted into heat in the collisional gas of the intergalactic medium. Using a toy model of gravitational heating we attempt to link the growth of structure in the Universe and the average temperature of this gas. We find that the gas is rapidly heated from collapsing structures at around $z {\sim} 10$, reaching a temperature ${>} 10^6$K today, depending on some assumptions of our simplified model. Before that there was a cold era from $z\sim100$ to $\sim10$ in which the matter temperature is below that of the Cosmic Microwave Background.

preprint2009arXiv

Matter temperature after cosmological recombination

The temperature of the atomic matter in the Universe is held to that of the Cosmic Background radiation until decoupling at z~100. After this it cools faster than the radiation (\propto(1+z)^2 rather than (1+z)) and would have fallen to about 20mK today if astrophysical feedback processes had not heated up the interglactic medium. We show how the derivative of the Compton coupling equation helps numerically to follow the decoupling process.

preprint2009arXiv

On the X-ray properties of submm-selected galaxies

We present an analysis of the X-ray properties of 35 submm galaxies(SMGs) in the CDF-N. Using a sample of robust 850micron-selected galaxies, with accurate positions from Spitzer and/or radio counterparts, we find 16 objects (45+/-8%) with significant X-ray detections in the 2Ms Chandra data. 6 of these SMGs (~17+/-6%) have measured X-ray luminosities or upper limits consistent with those expected based on the far-infrared (FIR) or radio-derived star formation rate (SFR), and hence with the X-rays coming solely from star formation. In another 7 sources (20+/-7%) a dominant AGN contribution to the X-ray emission is required, while in 3 more it is unclear whether stellar process or accretion are responsible. Stacking of the X-ray undetected SMGs reveals a highly significant detection. If due to star formation, this corresponds to an average X-ray derived SFR of ~150 Msun/yr. We deduce that the AGN fraction in SMGs based on X-ray observations is 20-29 (+/-7) %, which is towards the lower limit of previous estimates. Spectral analysis shows that in general the SMGs are not heavily obscured in the X-ray but most of the SMGs classfied as AGN show absorption with N_H in excess of 10^22 cm^-2. Of the secure AGN, the bolometric luminosity appears to be dominated by the AGN in only 3 cases. In ~85% of the SMGs, the X-ray spectrum effectively rules out an AGN contribution that dominates the bolometric emission, even if the AGN is Compton thick. The evidence therefore suggests that intense star formation accounts for both the FIR and X-ray emission in most SMGs. We argue that, rather than having an especially high AGN fraction or duty cycle, SMGs have a high X-ray detection rate at very faint fluxes partly because of their high star formation rates and, in rarer cases, because the submm emission is from an AGN.(abridged)

preprint2009arXiv

Strong lensing of submillimetre galaxies: A tracer of foreground structure?

The steep source counts and negative K-corrections of bright submillimetre galaxies (SMGs) suggest that a significant fraction of those observed at high flux densities may be gravitationally lensed, and that the lensing objects may often lie at redshifts above 1, where clusters of galaxies are difficult to detect through other means. In this case follow-up of bright SMGs may be used to identify dense structures along the line of sight. Here we investigate the probability for SMGs to experience strong lensing, using the latest N-body simulations and observed source flux and redshift distributions. We find that almost all high redshift sources with a flux density above 100mJy will be lensed, if they are not relatively local galaxies. We also give estimates of the fraction of sources experiencing strong lensing as a function of observed flux density. This has implications for planning follow-up observations for bright SMGs discovered in future surveys with SCUBA-2 and other instruments. The largest uncertainty in these calculations is the maximum allowed lensing amplification, which is dominated by the presently unknown spatial extent of SMGs.

preprint2009arXiv

The BLAST Survey of the Vela Molecular Cloud: Physical Properties of the Dense Cores in Vela-D

The Balloon-borne Large-Aperture Submillimeter Telescope (BLAST) carried out a 250, 350 and 500 micron survey of the galactic plane encompassing the Vela Molecular Ridge, with the primary goal of identifying the coldest dense cores possibly associated with the earliest stages of star formation. Here we present the results from observations of the Vela-D region, covering about 4 square degrees, in which we find 141 BLAST cores. We exploit existing data taken with the Spitzer MIPS, IRAC and SEST-SIMBA instruments to constrain their (single-temperature) spectral energy distributions, assuming a dust emissivity index beta = 2.0. This combination of data allows us to determine the temperature, luminosity and mass of each BLAST core, and also enables us to separate starless from proto-stellar sources. We also analyze the effects that the uncertainties on the derived physical parameters of the individual sources have on the overall physical properties of starless and proto-stellar cores, and we find that there appear to be a smooth transition from the pre- to the proto-stellar phase. In particular, for proto-stellar cores we find a correlation between the MIPS24 flux, associated with the central protostar, and the temperature of the dust envelope. We also find that the core mass function of the Vela-D cores has a slope consistent with other similar (sub)millimeter surveys.

preprint2006arXiv

Cosmic Conspiracies

The now standard vanilla-flavoured LambdaCDM model has gained further confirmation with the release of the 3-year WMAP data combined with several other cosmological data-sets. As the parameters of this standard model become known with increasing precision, more of its bizarre features become apparent. Here we describe some of the strangest of these ostensible coincidences. In particular we appear to live (within 1sigma) at the precise epoch when the age of the Universe multiplied by the Hubble parameter H_0 t_0 = 1.

preprint2006arXiv

The Hubble Deep Field North SCUBA Super-map IV - Characterizing submillimetre galaxies using deep Spitzer imaging

We present SEDs, Spitzer colours, and IR luminosities for 850 micron selected galaxies in the GOODS-N field. Using the deep Spitzer Legacy images and new data and reductions of the VLA-HDF radio data, we find statistically secure counterparts for 60 per cent (21/35) of our submm sample, and identify tentative counterparts for another 12 objects. This is the largest sample of submm galaxies with statistically secure counterparts detected in the radio and with Spitzer. We find that in most cases the 850 micron emission is dominated by a single 24 micron source. A composite rest-frame SED shows that the submm sources peak at longer wavelengths than those of local ULIRGs of the same luminosity and therefore appear to be cooler. The SEDs of submm galaxies are also different from those of their high redshift neighbours, the near-IR selected BzK galaxies, whose mid-IR to radio SEDs are more like those of local ULIRGs. Using 24 micron. 850 micron and 1.4 GHz observations, we fit templates that span the mid-IR through radio to derived the integrated IR luminosity of the submm galaxies. By themselves, 24 micron and radio fluxes are able to predict LIR reasonably well because they are relatively insensitive to temperature. However, the submm flux by itself consistently overpredicts LIR when using spectral templates which obey the local ULIRG temperature-luminosity relation. The shorter Spitzer wavelengths sample the stellar bump at the redshifts of the submm sources, and we find that the Spitzer photometry alone provides a model independent estimate of the redshift. Using X-ray and mid-IR data, only 5 per cent of our secure counterparts show strong evidence for an active galactic nucleus dominating the IR luminosity.

preprint2006arXiv

The Infrared Properties of Submillimeter Galaxies: Clues From Ultra-Deep 70 Micron Imaging

We present 70 micron properties of submillimeter galaxies (SMGs) in the Great Observatories Origins Deep Survey (GOODS) North field. Out of thirty submillimeter galaxies (S_850 > 2 mJy) in the central GOODS-N region, we find two with secure 70 micron detections. These are the first 70 micron detections of SMGs. One of the matched SMGs is at z ~ 0.5 and has S_70/S_850 and S_70/S_24 ratios consistent with a cool galaxy. The second SMG (z = 1.2) has infrared-submm colors which indicate it is more actively forming stars. We examine the average 70 micron properties of the SMGs by performing a stacking analysis, which also allows us to estimate that S_850 > 2 mJy SMGs contribute 9 +- 3% of the 70 micron background light. The S_850/S_70 colors of the SMG population as a whole is best fit by cool galaxies, and because of the redshifting effects these constraints are mainly on the lower z sub-sample. We fit Spectral Energy Distributions (SEDs) to the far-infrared data points of the two detected SMGs and the average low redshift SMG (z_{median}= 1.4). We find that the average low-z SMG has a cooler dust temperature than local ultraluminous infrared galaxies (ULIRGs) of similar luminosity and an SED which is best fit by scaled up versions of normal spiral galaxies. The average low-z SMG is found to have a typical dust temperature T = 21 -- 33 K and infrared luminosity L_{8-1000 micron} = 8.0 \times 10^11 L_sun. We estimate the AGN contribution to the total infrared luminosity of low-z SMGs is less than 23%.

preprint2005arXiv

The Standard Cosmological Model

The Standard Model of Particle Physics (SMPP) is an enormously successful description of high energy physics, driving ever more precise measurements to find "physics beyond the standard model", as well as providing motivation for developing more fundamental ideas that might explain the values of its parameters. Simultaneously, a description of the entire 3-dimensional structure of the present-day Universe is being built up painstakingly. Most of the structure is stochastic in nature, being merely the result of the particular realisation of the "initial conditions" within our observable Universe patch. However, governing this structure is the Standard Model of Cosmology (SMC), which appears to require only about a dozen parameters. Cosmologists are now determining the values of these quantities with increasing precision in order to search for "physics beyond the standard model", as well as trying to develop an understanding of the more fundamental ideas which might explain the values of its parameters. Although it is natural to see analogies between the two Standard Models, some intrinsic differences also exist, which are discussed here. Nevertheless, a truly fundamental theory will have to explain both the SMPP and SMC, and this must include an appreciation of which elements are deterministic and which are accidental. Considering different levels of stochasticity within cosmology may make it easier to accept that physical parameters in general might have a non-deterministic aspect.

preprint2004arXiv

An 850-micron SCUBA map of the Groth Strip and reliable source extraction

We present an 850-micron map and list of candidate sources in a 70 arcmin^2 sub-area of the Groth Strip observed using SCUBA. We initially detect 7 candidate sources with signal-to-noise ratios (SNRs) between 3.0 and 3.5 and 4 candidate sources with SNR > 3.5. Simulations suggest that on average in a map this size one expects 1.6 false positive sources for SNR > 3.5 and 4.5 for 3 < SNR < 3.5. Flux boosting in maps is a well known effect and we have developed a simple Bayesian prescription for estimating the unboosted flux distribution and used this method to determine the best flux estimates of our sources. This method is easily adapted for any other modest signal-to-noise survey in which there is prior knowledge of the source counts. We performed follow-up photometry in an attempt to confirm or reject 5 of our source candidates. We confirm the reality of 2 of the SCUBA sources, although at lower levels than suggested in the map and we find that the photometry results are consistent with and confirm the de-boosted map fluxes. Our final candidate source list contains 3 sources, including the 2 confirmed detections and 1 further candidate source with SNR > 3.5 which has a reasonable chance of being real.

preprint1996arXiv

Cold dark matter models with high baryon content

Recent results have suggested that the density of baryons in the Universe, OmegaB, is much more uncertain than previously thought, and may be significantly higher. We demonstrate that a higher OmegaB increases the viability of critical-density cold dark matter (CDM) models. High baryon fraction offers the twin benefits of boosting the first peak in the microwave anisotropy power spectrum and of suppressing short-scale power in the matter power spectrum. These enable viable CDM models to have a larger Hubble constant than otherwise possible. We carry out a general exploration of high OmegaB CDM models, varying the Hubble constant h and the spectral index n. We confront a variety of observational constraints and discuss specific predictions. Although some observational evidence may favour baryon fractions as high as 20 per cent, we find that values around 10 to 15 per cent provide a reasonable fit to a wide range of data. We suggest that models with OmegaB in this range, with h about 0.5 and n about 0.8, are currently the best critical-density CDM models.

preprint1995arXiv

Cold Dark Matter Resuscitated?

The Cold Dark Matter (CDM) model has an elegant simplicitly which makes it very predictive, but when its parameters are fixed at their `canonical' values its predictions are in conflict with observational data. There is, however, much leeway in the initial conditions within the CDM framework. We advocate a re-examination of the CDM model, taking into account modest variation of parameters from their canonical values. We find that CDM models with $n=0.8$--0.9 and $h=0.45$--0.50 can fit the available data. Our ``best fit'' CDM model has $n=0.9$, $h=0.45$ and $C_2^{T}/C_2^{S}=0.7$. We discuss the current state of observations which could definitely rule out this model.