Source author record

M. J. Page

M. J. Page 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

145works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

145 published item(s)

preprint2022arXiv

The star-formation rates of QSOs

We examine the far-IR properties of a sample of 5391 optically selected QSOs in the 0.5<z<2.65 redshift range down to log[nuLnu,2500 (erg/s)]>44.7, using SPIRE data from Herschel-ATLAS. We split the sample in a grid of 74 luminosity-redshift bins and compute the average optical-infrared spectral energy distribution (SED) in each bin. By normalising an intrinsic AGN template to the AGN optical power (at 5100A) we decompose the total infrared emission (L_IR; 8-1000um) into an AGN (L_IR,AGN) and star-forming component (L_IR,SF). We find that the AGN contribution to L_IR increases as a function of AGN power which manifests as a reduction of the `far-IR bump' in the average QSO SEDs. We note that L_IR,SF does not correlate with AGN power; the mean star formation rates (SFRs) of AGN host galaxies are a function of redshift only and they range from ~6 Msun/yr at z~0 to a plateau of <200 Msun/yr at z~2.6. Our results indicate that the accuracy of far-IR emission as a proxy for SFR decreases with increasing AGN luminosity. We show that, at any given redshift, observed trends between infrared luminosity (whether monochromatic or total) and AGN power (in the optical or X-rays) can be explained by a simple model which is the sum of two components: (A) the infrared emission from star-formation, uncorrelated with AGN power and (B) the infrared emission from AGN, directly proportional to AGN power in the optical or X-rays.

preprint2020arXiv

Swift Multiwavelength Follow-up of LVC S200224ca and the Implications for Binary Black Hole Mergers

On 2020 February 24, during their third observing run ("O3"), the Laser Interferometer Gravitational-wave Observatory and Virgo Collaboration (LVC) detected S200224ca: a candidate gravitational wave (GW) event produced by a binary black hole (BBH) merger. This event was one of the best-localized compact binary coalescences detected in O3 (with 50%/90% error regions of 13/72 deg$^2$), and so the Neil Gehrels Swift Observatory performed rapid near-UV/X-ray follow-up observations. Swift-XRT and UVOT covered approximately 79.2% and 62.4% (respectively) of the GW error region, making S200224ca the BBH event most thoroughly followed-up in near-UV (u-band) and X-ray to date. No likely EM counterparts to the GW event were found by the Swift BAT, XRT, or UVOT, nor by other observatories. Here we report on the results of our searches for an EM counterpart, both in the BAT data near the time of the merger, and in follow-up UVOT/XRT observations. We also discuss the upper limits we can place on EM radiation from S200224ca, and the implications these limits have on the physics of BBH mergers. Namely, we place a shallow upper limit on the dimensionless BH charge, $\hat{q} < 1.4 \times10^{-4}$, and an upper limit on the isotropic-equivalent energy of a blast wave $E < 4.1\times10^{51}$ erg (assuming typical GRB parameters).

preprint2020arXiv

The XMM-Newton serendipitous survey IX. The fourth XMM-Newton serendipitous source catalogue

Sky surveys produce enormous quantities of data on extensive regions of the sky. The easiest way to access this information is through catalogues of standardised data products. {\em XMM-Newton} has been surveying the sky in the X-ray, ultra-violet, and optical bands for 20 years. The {\em XMM-Newton} Survey Science Centre has been producing standardised data products and catalogues to facilitate access to the serendipitous X-ray sky. Using improved calibration and enhanced software, we re-reduced all of the 14041 {\em XMM-Newton} X-ray observations, of which 11204 observations contained data with at least one detection and with these we created a new, high quality version of the {\em XMM-Newton} serendipitous source catalogue, 4XMM-DR9. 4XMM-DR9 contains 810795 detections down to a detection significance of 3 $σ$, of which 550124 are unique sources, which cover 1152 degrees$^{2}$ (2.85\%) of the sky. Filtering 4XMM-DR9 to retain only the cleanest sources with at least a 5 $σ$ detection significance leaves 433612 detections. Of these detections, 99.6\% have no pileup. Furthermore, 336 columns of information on each detection are provided, along with images. The quality of the source detection is shown to have improved significantly with respect to previous versions of the catalogues. Spectra and lightcurves are also made available for more than 288000 of the brightest sources (36\% of all detections).

preprint2020arXiv

The XMM-Newton serendipitous survey. X: The second source catalogue from overlapping XMM-Newton observations and its long-term variable content

The XMM-Newton Survey Science Centre Consortium (SSC) develops software in close collaboration with the Science Operations Centre to perform a pipeline analysis of all XMM-Newton observations. In celebration of the 20th launch anniversary, the SSC has compiled the 4th generation of serendipitous source catalogues, 4XMM. The catalogue described here, 4XMM-DR9s, explores sky areas that were observed more than once by XMM-Newton. It was constructed from simultaneous source detection on the overlapping observations, which were bundled in groups ("stacks"). Stacking leads to a higher sensitivity, resulting in newly discovered sources and better constrained source parameters, and unveils long-term brightness variations. As a novel feature, positional rectification was applied beforehand. Observations with all filters and suitable camera settings were included. Exposures with a high background were discarded, which was determined through a statistical analysis of all exposures in each instrument configuration. The X-ray background maps used in source detection were modelled via adaptive smoothing with newly determined parameters. Source fluxes were derived for all contributing observations, irrespective of whether the source would be detectable in an individual observation. From 1,329 stacks with 6,604 contributing observations over repeatedly covered 300 square degrees in the sky, 4XMM-DR9s lists 288,191 sources. 218,283 of them were observed several times. Most stacks are composed of two observations, the largest one comprises 352. The number of observations of a source ranges from 1 to 40. Auxiliary products like X-ray images, long-term light curves, and optical finding charts are published as well. 4XMM-DR9s is considered a prime resource to explore long-term variability of X-ray sources discovered by XMM-Newton. Regular incremental releases including new public observations are planned.

preprint2019arXiv

Multi-wavelength Campaign on NGC 7469 VI. Photoionisation Modelling of the Emission Line Regions and the Warm Absorber

Aims. We aim to investigate and characterise the photoionised X-ray emission line regions within NGC 7469. Methods. We apply the photoionisation model, PION, within the spectral fitting code SPEX to analyse the 640 ks RGS spectrum of NGC 7469 gathered during an XMM-Newton observing campaign in 2015. Results. We find the emission line region in NGC 7469 to be multiphased, consisting of two narrow components with ionisation parameters of $\log ξ= 0.4$ and 1.6. A third, broad emission component, with a broadening velocity of $v_b \sim 1400$ km s\textsuperscript{-1} and an outflow velocity of $v_{out} \sim -4500$ km s\textsuperscript{-1}, is required to fit the residuals in the O VII triplet, at around 22 Å. Assuming a volume filling factor of 0.1, the lower distance limits of the narrow emission line region components are estimated for the first time at 2.6 and 2.5 pc from the central black hole, whereas the broad component has an estimated lower bound distance between 0.004 to 0.03 pc, depending on the assumed plasma parameters. The collisionally ionised plasma from the star burst region in NGC 7469 has a plasma temperature of 0.32 keV and outflow velocity of $-280$ km s\textsuperscript{-1}, consistent with previous results in this campaign. In addition, we model the photoionised plasma of the warm absorber (WA) in NGC 7469, and find that it consists of three photoionised phases, with different values of $ξ$, $N_H$ and $v_{out}$. The upper bound distances of these WA components are 1.9, 0.3 and 0.6 pc, respectively, consistent with archival results. Conclusions. The environment of NGC 7469 is a complex mix of plasma winds absorbing and emitting X-rays. We find the picture painted by our results can be attributed to line emitting plasma located at distances ranging from near the black hole to the torus and beyond the ionised outflows.

preprint2016arXiv

AGN are cooler than you think: the intrinsic far-IR emission from QSOs

We present an intrinsic AGN SED extending from the optical to the submm, derived with a sample of unobscured, optically luminous (vLv(5100)>10^43.5 erg/s) QSOs at z<0.18 from the Palomar Green survey. The intrinsic AGN SED was computed by removing the contribution from stars using the 11.3um polycyclic aromatic hydrocarbon (PAH) feature in the QSOs' mid-IR spectra; the 1sigma uncertainty on the SED ranges between 12 and 45 per cent as a function of wavelength and is a combination of PAH flux measurement errors and the uncertainties related to the conversion between PAH luminosity and star-forming luminosity. Longwards of 20um the shape of the intrinsic AGN SED is independent of the AGN power indicating that our template should be applicable to all systems hosting luminous AGN (vLv(5100) or L_X(2-10keV) > 10^43.5 erg/s). We note that for our sample of luminous QSOs, the average AGN emission is at least as high as, and mostly higher than, the total stellar powered emission at all wavelengths from the optical to the submm. This implies that in many galaxies hosting powerful AGN, there is no `safe' broadband photometric observation (at lambda<1000um) which can be used in calculating star-formation rates without subtracting the AGN contribution. Roughly, the AGN contribution may be ignored only if the intrinsic AGN luminosity at 5100 Ang is at least a factor of 4 smaller than the total infrared luminosity (L_IR; 8-1000um) of the galaxy. Finally, we examine the implication of our work in statistical studies of star-formation in AGN host galaxies.

preprint2016arXiv

The 80 Ms follow-up of the X-ray afterglow of GRB 130427A challenges the standard forward shock model

GRB 130427A was the brightest gamma-ray burst detected in the last 30 years. With an equivalent isotropic energy output of $8.5\times10^{53}$ erg and redshift $z=0.34$, it uniquely combined very high energetics with a relative proximity to Earth. As a consequence, its X-ray afterglow has been detected by sensitive X-ray observatories such as XMM-Newton and Chandra for a record-breaking baseline longer than 80 million seconds. We present the X-ray light-curve of this event over such an interval. The light-curve shows a simple power-law decay with a slope $α= 1.309 \pm 0.007$ over more than three decades in time (47 ks - 83 Ms). We discuss the consequences of this result for a few models proposed so far to interpret GRB 130427A, and more in general the significance of this outcome in the context of the standard forward shock model. We find that this model has difficulty in explaining our data, in both cases of constant density and stellar wind circumburst media, and requires far-fetched values for the physical parameters involved.

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

X-ray narrow emission lines from the nuclear region of NGC 1365

NGC 1365 is a Seyfert 2 galaxy with a starburst ring in its nuclear region. In this work we look at the XMM Reflection Grating Spectrometer (RGS) data from four 2012-13, three 2007 and two 2004 observations of NGC 1365. We characterise the narrow-line emitting gas visible by XMM RGS and make comparisons between the 2012-13 spectra and those from 2004-07, already published. This source is usually absorbed within the soft X-ray band, with a typical neutral column density of >1.5 x 10$^{23}$ cm$^{-2}$, and only 1 observation of the 9 we investigate shows low enough absorption for the continuum to emerge in the soft X-rays. We stack all observations from 2004-07, and separately three of the four observations from 2012-13, analysing the less absorbed observation separately. We first model the spectra using gaussian profiles representing the narrow line emission. We fit physically motivated models to the 2012-13 stacked spectra, with collisionally ionised components representing the starburst emission and photoionised line emission models representing the AGN line emission. The collisional and photoionised emission line models are fitted together on top of a physical continuum and absorption model. The X-ray narrow emission line spectrum of NGC 1365 is well represented by a combination of two collisionally ionised and three photoionised phases of emitting gas, all with higher than solar nitrogen abundances. This physical model was fitted to the 2012-13 stacked spectrum, and yet also fits well to the 2004-07 stacked spectrum, without changing any characteristics of the emitting gas phases. Our 2004-07 results are consistent with previous emission line work using these data, with 5 additional emission lines detected in both this and the 2012-13 stacked spectra. We also estimate the distance of the X-ray line-emitting photoionised gas from the central source to be <300 pc.

preprint2015arXiv

Exploring the canonical behaviour of long gamma-ray bursts using an intrinsic multi-wavelength afterglow correlation

In this paper we further investigate the relationship, reported by Oates et al., 2012, between the optical/UV afterglow luminosity (measured at restframe 200s) and average afterglow decay rate (measured from restframe 200s onwards) of long duration Gamma-ray Bursts (GRBs). We extend the analysis by examining the X-ray light curves, finding a consistent correlation. We therefore explore how the parameters of these correlations relate to the prompt emission phase and, using a Monte Carlo simulation, explore whether these correlations are consistent with predictions of the standard afterglow model. We find significant correlations between: $\rm log\;L_{O,200\rm{s}}$ and $\rm log\;L_{X,200\rm{s}}$; $α_{O,>200\rm{s}}$ and $α_{X,>200\rm{s}}$, consistent with simulations. The model also predicts relationships between $\rm log\;E_{iso}$ and $\rm log\;L_{200\rm{s}}$, however, while we find such relationships in the observed sample, the slope of the linear regression is shallower than that simulated and inconsistent at $\gtrsim 3σ$. Simulations also do not agree with correlations observed between $\rm log\;L_{200\rm{s}}$ and $α_{>200\rm{s}}$, or $\rm log\;E_{iso}$ and $α_{>200\rm{s}}$. Overall, these observed correlations are consistent with a common underlying physical mechanism producing GRBs and their afterglows regardless of their detailed temporal behaviour. However, a basic afterglow model has difficulty explaining all the observed correlations. This leads us to briefly discuss alternative more complex models.

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

Submm-bright X-ray absorbed QSOs at z~2: insights into the co-evolution of AGN and star-formation

We have assembled a sample of 5 X-ray-absorbed and submm-luminous type 1 QSOs at $z \sim 2$ which are simultaneously growing their central black holes through accretion and forming stars copiously. We present here the analysis of their rest-frame UV to submm Spectral Energy Distributions (SEDs), including new Herschel data. Both AGN (direct and reprocessed) and Star Formation (SF) emission are needed to model their SEDs. From the SEDs and their UV-optical spectra we have estimated the masses of their black holes $M_{BH}\sim 10^{9}-10^{10}\,M_{\odot}$, their intrinsic AGN bolometric luminosities $L_{BOL}\sim(0.8 - 20)\times 10^{13} L_{\odot}$, Eddington ratios $L_{BOL}/L_{Edd}\sim 0.1 - 1.1$ and bolometric corrections $L_{BOL}/L_{X,2-10}\sim 30 - 500$. These values are common among optically and X-ray-selected type 1 QSOs (except for RX~J1249), except for the bolometric corrections, which are higher. These objects show very high far-infrared luminosities $L_{FIR}\sim$ (2 - 8)$\times10^{12}\,M_{\odot}$ and Star Formation Rates SFR$\sim 1000 M_{\odot}/$y. From their $L_{FIR}$ and the shape of their FIR-submm emission we have estimated star-forming dust masses of $M_{DUST}\sim 10^9\,M_\odot$. We have found evidence of a tentative correlation between the gas column densities of the ionized absorbers detected in X-ray (N$_{H_{ion}}$) and $SFR$. Our computed black hole masses are amongst the most massive known.

preprint2015arXiv

The central engine of GRB 130831A and the energy breakdown of a relativistic explosion

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

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 Swift-UVOT ultraviolet and visible grism calibration

We present the calibration of the Swift UVOT grisms, of which there are two, providing low-resolution field spectroscopy in the ultraviolet and optical bands respectively. The UV grism covers the range 1700-5000 Angstrom with a spectral resolution of 75 at 2600 Angstrom for source magnitudes of u=10-16 mag, while the visible grism covers the range 2850-6600 Angstrom with a spectral resolution of 100 at 4000 Angstrom for source magnitudes of b=12-17 mag. This calibration extends over all detector positions, for all modes used during operations. The wavelength accuracy (1-sigma) is 9 Angstrom in the UV grism clocked mode, 17 Angstrom in the UV grism nominal mode and 22 Angstrom in the visible grism. The range below 2740 Angstrom in the UV grism and 5200 Angstrom in the visible grism never suffers from overlapping by higher spectral orders. The flux calibration of the grisms includes a correction we developed for coincidence loss in the detector. The error in the coincidence loss correction is less than 20%. The position of the spectrum on the detector only affects the effective area (sensitivity) by a few percent in the nominal modes, but varies substantially in the clocked modes. The error in the effective area is from 9% in the UV grism clocked mode to 15% in the visible grism clocked mode .

preprint2014arXiv

A panchromatic analysis of starburst galaxy M82: Probing the dust properties

(Abridged) We combine NUV, optical and IR imaging of the nearby starburst galaxy M82 to explore the properties of the dust both in the interstellar medium of the galaxy and the dust entrained in the superwind. The three NUV filters of Swift/UVOT enable us to probe in detail the properties of the extinction curve in the region around the 2175A bump. The NUV colour-colour diagram strongly rules out a Calzetti-type law, which can either reflect intrinsic changes in the dust properties or in the star formation history compared to starbursts well represented by such an attenuation law. We emphasize that it is mainly in the NUV region where a standard Milky-Way-type law is preferred over a Calzetti law. The age and dust distribution of the stellar populations is consistent with the scenario of an encounter with M81 in the recent 400 Myr. The radial gradients of the NUV and optical colours in the superwind region support the hypothesis that the emission in the wind cone is driven by scattering from dust grains entrained in the ejecta. The observed wavelength dependence reveals either a grain size distribution $n(a)\propto a^{-2.5}$, where $a$ is the size of the grain, or a flatter distribution with a maximum size cutoff, suggesting that only small grains are entrained in the supernovae-driven wind.

preprint2014arXiv

Herschel/PACS Observations of the Host Galaxy of GRB 031203

We present Herschel/PACS observations of the nearby (z=0.1055) dwarf galaxy that has hosted the long gamma ray burst (LGRB) 031203. Using the PACS data we have been able to place constraints on the dust temperature, dust mass, total infrared luminosity and infrared-derived star-formation rate (SFR) for this object. We find that the GRB host galaxy (GRBH) 031203 has a total infrared luminosity of 3x10^10 L_sun placing it in the regime of the IR-luminous galaxy population. Its dust temperature and specific SFR are comparable to that of many high-redshift (z=0.3-2.5) infrared (IR)-detected GRB hosts (T_dust>40K ; sSFR>10 Gyr^-1), however its dust-to-stellar mass ratio is lower than what is commonly seen in IR-luminous galaxies. Our results suggest that GRBH 031203 is undergoing a strong starburst episode and its dust properties are different to those of local dwarf galaxies within the same metallicity and stellar mass range. Furthermore, our measurements place it in a distinct class to the well studied nearby host of GRB 980425 (z=0.0085), confirming the notion that GRB host galaxies can span a large range in properties even at similar cosmological epochs, making LGRBs an ideal tool in selecting samples of star-forming galaxies up to high redshift.

preprint2014arXiv

Linking the X-ray and infrared properties of star-forming galaxies at z<1.5

We present the most complete study to date of the X-ray emission from star-formation in high redshift (median z=0.7; z<1.5), IR-luminous (L_IR=10^10-10^13 L_sun) galaxies detected by Herschel's PACS and SPIRE instruments. For our purpose we take advantage of the deepest X-ray data to date, the Chandra deep fields (North and South). Sources which host AGN are removed from our analysis by means of multiple AGN indicators. We find an AGN fraction of 18+/-2 per cent amongst our sample and note that AGN entirely dominate at values of log[L_X/L_IR]>-3 in both hard and soft X-ray bands. From the sources which are star-formation dominated, only a small fraction are individually X-ray detected and for the bulk of the sample we calculate average X-ray luminosities through stacking. We find an average soft X-ray to infrared ratio of log[L_SX/L_IR]=-4.3 and an average hard X-ray to infrared ratio of log[L_HX/L_IR]=-3.8. We report that the X-ray/IR correlation is approximately linear through the entire range of L_IR and z probed and, although broadly consistent with the local (z<0.1) one, it does display some discrepancies. We suggest that these discrepancies are unlikely to be physical, i.e. due to an intrinsic change in the X-ray properties of star-forming galaxies with cosmic time, as there is no significant evidence for evolution of the L_X/L_IR ratio with redshift. Instead they are possibly due to selection effects and remaining AGN contamination. We also examine whether dust obscuration in the galaxy plays a role in attenuating X-rays from star-formation, by investigating changes in the L_X/L_IR ratio as a function of the average dust temperature. We conclude that X-rays do not suffer any measurable attenuation in the host galaxy.

preprint2014arXiv

Traces of co-evolution in high z X-ray selected and submm-luminous QSOs

We present a detailed study of a X -ray selected sample of 5 submillimeter bright QSOs at $z\sim2$, where the highest rates of star formation (SF) and further growth of black holes (BH) occur. Therefore, this sample is a great laboratory to investigate the co-evolution of star formation and AGN. We present here the analysis of the spectral energy distributions (SED) of the 5 QSOS, including new data from Herschel PACS and SPIRE. Both AGN components (direct and reprocessed) and like Star Formation (SF) are needed to model its SED. From the SED and their UV-optical spectra we have estimated the mass of the black hole ($M_{BH} = 10^9 - 10^{10} M_{SUN}$) and bolometric luminosities of AGN ($L_{BOL} = (0.8-20) \times 10^{13} L_{SUN}$). These objects show very high luminosities in the far infrared range (at the H/ULIRG levels) and very high rates of SF (SFR = 400-1400 $M_{SUN}$/y). Known their current SFR and their BH masses, we deduce that their host galaxies must be already very massive, or would not have time to get to the local relation between BH mass and bulge. Finally, we found evidence of a possible correlation between the column density of ionized gas detected in X-rays ($NH_{ion}$) and SFR, which would provide a link between AGN and SF processes.

preprint2014arXiv

X-rays from the redshift 7.1 quasar ULAS J1120+0641

We present X-ray imaging and spectroscopy of the redshift z=7.084 radio-quiet quasar ULAS J112001.48+064124.3 obtained with Chandra and XMM-Newton. The quasar is detected as a point source with both observatories. The Chandra observation provides a precise position, confirming the association of the X-ray source and the quasar, while a sufficient number of photons is detected in the XMM-Newton observation to yield a meaningful X-ray spectrum. In the XMM-Newton observation the quasar has a 2-10 keV luminosity of 4.7+-0.9 times 10^44 erg/s and a spectral slope alpha = 1.6+0.4/-0.3 (where f_nu is proportional to nu^-alpha). The quasar appears to have dimmed in the 15 months between the two observations, with a 2-10 keV luminosity of $1.8+1.0/-0.7 times 10^45 erg/s during the Chandra observation. We derive optical to X-ray spectral slopes alpha_OX of 1.76+-0.07 and 1.54+0.09/-0.08 at the times of the XMM-Newton and Chandra observations respectively, consistent with the range of alpha_OX found in other quasars of comparable ultraviolet luminosity. The very soft X-ray spectrum suggests that the quasar is accreting above the Eddington rate, L/L_Edd = 5+15/-4, compared to L/L_Edd = 1.2+0.6/-0.5 derived from the rest-frame ultraviolet. Super-Eddington accretion would help to reduce the discrepancy between the age of the quasar implied by the small size of the ionized near zone in which it sits (<10^7 years), and the characteristic e-folding time (2.5 times 10^7 years if L/L_Edd=2). Such super-Eddington accretion would also alleviate the challenging constraints on the seed black hole mass provided that the quasar has been rapidly accreting throughout its history. The remnant of an individual population III star is a plausible progenitor if an average L/L_Edd>1.46 has been maintained over the quasar's lifetime.

preprint2013arXiv

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

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

preprint2013arXiv

A multiwavelength study of the Magellanic-type galaxy NGC 4449 - I. Modelling the spectral energy distribution, the ionization structure and the star formation history

[Abridged] We present an integrated photometric spectral energy distribution (SED) of the Magellanic-type galaxy NGC 4449 from the far-ultraviolet (UV) to the submillimetre, including new observations acquired by the Herschel Space Observatory. We include integrated UV photometry from the Swift Ultraviolet and Optical Telescope using a measurement technique which is appropriate for extended sources with coincidence loss. In this paper, we examine the available multiwavelength data to infer a range of ages, metallicities and star formation rates for the underlying stellar populations, as well as the composition and the total mass of dust in NGC 4449. We present an iterative scheme, which allows us to build an in-depth and multicomponent representation of NGC 4449 `bottom-up', taking advantage of the broad capabilities of the photoionization and radiative transfer code MOCASSIN (MOnte CArlo SimulationS of Ionized Nebulae). We fit the observed SED, the global ionization structure and the emission line intensities, and infer a recent SFR of 0.4 Msolar/yr and a total stellar mass of approximately 1e9 Msolar emitting with a bolometric luminosity of 5.7e9 Lsolar. Our fits yield a total dust mass of 2.9e6 Msolar including 2 per cent attributed to polycyclic aromatic hydrocarbons. We deduce a dust to gas mass ratio of 1/190 within the modelled region. While we do not consider possible additional contributions from even colder dust, we note that including the extended HI envelope and the molecular gas is likely to bring the ratio down to as low as ~ 1/800.

preprint2013arXiv

AGN in dusty hosts: implications for galaxy evolution

We present strong empirical evidence for a physical connection between the occurrence of a starburst (SB) and a luminous AGN phase. Drawing infrared (IR), X-ray, and optically selected samples from COSMOS, we find that the locus of type-2 AGN hosts in the optical colour-magnitude (U-V/V) and colour-colour (U-V/V-J) space significantly overlaps with that of IR-luminous (L_IR > 10^10 L_sun) galaxies. Based on our observations, we propose that, when simultaneously building their black hole and stellar masses, type-2 AGN hosts are located in the same part of colour-colour space as dusty star-forming galaxies. In fact, our results show that IR-luminous galaxies at z<1.5 are on average 3 times more likely to host a type-2 AGN (L_X > 10^42 erg/s) than would be expected serendipitously, if AGN and star-formation events were unrelated. In addition, the optical and infrared properties of the AGN/SB hybrid systems tentatively suggest that the AGN phase might be coeval with a particularly active phase in a galaxy's star-formation history. Interestingly, we also find a significant fraction of type-2 AGN hosts offset from the dusty galaxy sequence in colour-colour space, possibly representing a transitional or post-starburst phase in galaxy evolution. Our findings are consistent with a scenario whereby AGN play a role in the termination of star-formation in massive galaxies.

preprint2013arXiv

Chronos: A NIR spectroscopic galaxy survey. From the formation of galaxies to the peak of activity

Chronos is our response to ESA's call for white papers to define the science for the future L2, L3 missions. Chronos targets the formation and evolution of galaxies, by collecting the deepest NIR spectroscopic data, from the formation of the first galaxies at z~10 to the peak of formation activity at z~1-3. The strong emission from the atmospheric background makes this type of survey impossible from a ground-based observatory. The spectra of galaxies represent the equivalent of a DNA fingerprint, containing information about the past history of star formation and chemical enrichment. The proposed survey will allow us to dissect the formation process of galaxies including the timescales of quenching triggered by star formation or AGN activity, the effect of environment, the role of infall/outflow processes, or the connection between the galaxies and their underlying dark matter haloes. To provide these data, the mission requires a 2.5m space telescope optimised for a campaign of very deep NIR spectroscopy. A combination of a high multiplex and very long integration times will result in the deepest, largest, high-quality spectroscopic dataset of galaxies from z=1 to 12, spanning the history of the Universe, from 400 million to 6 billion years after the big bang, i.e. covering the most active half of cosmic history.

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: Cosmic Infrared Background Anisotropies and the Clustering of Dusty Star-Forming Galaxies

We present measurements of the auto- and cross-frequency power spectra of the cosmic infrared background (CIB) at 250, 350, and 500um (1200, 860, and 600 GHz) from observations totaling ~ 70 deg^2 made with the SPIRE instrument aboard the Herschel Space Observatory. We measure a fractional anisotropy dI / I = 14 +- 4%, detecting signatures arising from the clustering of dusty star-forming galaxies in both the linear (2-halo) and non-linear (1-halo) regimes; and that the transition from the 2- to 1-halo terms, below which power originates predominantly from multiple galaxies within dark matter halos, occurs at k_theta ~ 0.1 - 0.12 arcmin^-1 (l ~ 2160 - 2380), from 250 to 500um. New to this paper is clear evidence of a dependence of the Poisson and 1-halo power on the flux-cut level of masked sources --- suggesting that some fraction of the more luminous sources occupy more massive halos as satellites, or are possibly close pairs. We measure the cross-correlation power spectra between bands, finding that bands which are farthest apart are the least correlated, as well as hints of a reduction in the correlation between bands when resolved sources are more aggressively masked. In the second part of the paper we attempt to interpret the measurements in the framework of the halo model. With the aim of fitting simultaneously with one model the power spectra, number counts, and absolute CIB level in all bands, we find that this is achievable by invoking a luminosity-mass relationship, such that the luminosity-to-mass ratio peaks at a particular halo mass scale and declines towards lower and higher mass halos. Our best-fit model finds that the halo mass which is most efficient at hosting star formation in the redshift range of peak star-forming activity, z ~ 1-3, is log(M_peak/M_sun) ~ 12.1 +- 0.5, and that the minimum halo mass to host infrared galaxies is log(M_min/M_sun) ~ 10.1 +- 0.6.

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

Herschel reveals the obscured star formation in HiZELS Hα emitters at z=1.47

We describe the far-infrared (FIR; rest-frame 8--1000μm) properties of a sample of 443 Hα-selected star-forming galaxies in the COSMOS and UDS fields detected by the HiZELS imaging survey. Sources are identified using narrow-band filters in combination with broad-band photometry to uniformly select Hα (and [OII] if available) emitters in a narrow redshift slice at z = 1.47+/-0.02. We use a stacking approach in Spitzer, Herschel (from PEP and HerMES surveys) and AzTEC images to describe their typical FIR properties. We find that HiZELS galaxies with observed Hα luminosities of ~ 10^{8.1-9.1} Lo have bolometric FIR luminosities of typical LIRGs, L_FIR ~ 10^{11.48+/-0.05} Lo. Combining the Hα and FIR luminosities, we derive median SFR = 32+/-5 Mo/yr and Hα extinctions of A(Hα) = 1.0+/-0.2 mag. Perhaps surprisingly, little difference is seen in typical HiZELS extinction levels compared to local star-forming galaxies. We confirm previous empirical stellar mass (M*) to A(Hα) relations and the little or no evolution up to z = 1.47. For HiZELS galaxies, we provide an empirical parametrisation of the SFR as a function of (u-z)_rest colours and 3.6μm photometry. We find that the observed Hα luminosity is a dominant SFR tracer when (u-z)_rest ~< 0.9 mag or when 3.6μm photometry > 22 mag (Vega) or when M* < 10^9.7 Mo. We do not find any correlation between the [OII]/Hα and FIR luminosity, suggesting that this emission line ratio does not trace the extinction of the most obscured star-forming regions. The luminosity-limited HiZELS sample tends to lie above of the so-called `main sequence' for star-forming galaxies, especially at low M*. This work suggests that obscured star formation is linked to the assembly of M*, with deeper potential wells in massive galaxies providing dense, heavily obscured environments in which stars can form rapidly.

preprint2013arXiv

On the redshift evolution (0 < z < 4) of dust attenuation and of the total (UV+IR) star formation rate density

Using new homogeneous LFs in the FUV and in the FIR Herschel/PEP and Herschel/HerMES, we study the evolution of the dust attenuation with redshift. With this information in hand, we are able to estimate the redshift evolution of the total (FUV + FIR) star formation rate density SFRD_TOT. By integrating SFRD_TOT, we follow the mass building and analyze the redshift evolution of the stellar mass density (SMD). This letter aims at providing a complete view of star formation from the local universe to z = 4 and, using assumptions on earlier star formation history, compares this evolution to what was known before in an attempt to draw a homogeneous picture of the global evolution of star formation in galaxies. The main conclusions of this letter are: 1) the dust attenuation A_FUV is found to increase from z = 0 to z \sim 1.2 and then starts to decrease up to our last data point at z = 3.6; 2) the estimated SFRD confirms published results up to z = 2. At z > 2, we observe either a plateau or a small increase up to z = 3 and then a likely decrease up to z = 3.6; 3) the peak of A_FUV is delayed with respect to the plateau of SFRD_TOT and a likely origin might be found in the evolution of the bright ends of the FUV and FIR LFs; 4) using assumptions (namely exponential rise and linear rise with time) for the evolution of the star formation density from z = 3.6 to z_form = 10, we integrate SFRD_TOT and find a good agreement with the published SMDs.

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

Submm-bright QSOs at z~2: signposts of co-evolution at high z

We have assembled a sample of 5 X-ray and submm-luminous z~2 QSOs which are therefore both growing their central black holes through accretion and forming stars copiously at a critical epoch. Hence, they are good laboratories to investigate the co-evolution of star formation and AGN. We have performed a preliminary analysis of the AGN and SF contributions to their UV-to-FIR SEDs, fitting them with simple direct (disk), reprocessed (torus) and star formation components. All three are required by the data and hence we confirm that these objects are undergoing strong star formation in their host galaxies at rates 500-2000 Msun/y. Estimates of their covering factors are between about 30 and 90%. In the future, we will assess the dependence of these results on the particular models used for the components and relate their observed properties to the intrinsice of the central engine and the SF material, as well as their relevance for AGN-galaxy coevolution.

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 Herschel PEP/HerMES Luminosity Function. I: Probing the Evolution of PACS selected Galaxies to z~4

We exploit the deep and extended far infrared data sets (at 70, 100 and 160 um) of the Herschel GTO PACS Evolutionary Probe (PEP) Survey, in combination with the HERschel Multi tiered Extragalactic Survey (HerMES) data at 250, 350 and 500 um, to derive the evolution of the restframe 35 um, 60 um, 90 um, and total infrared (IR) luminosity functions (LFs) up to z~4. We detect very strong luminosity evolution for the total IR LF combined with a density evolution. In agreement with previous findings, the IR luminosity density increases steeply to z~1, then flattens between z~1 and z~3 to decrease at z greater than 3. Galaxies with different SEDs, masses and sSFRs evolve in very different ways and this large and deep statistical sample is the first one allowing us to separately study the different evolutionary behaviours of the individual IR populations contributing to the IR luminosity density. Galaxies occupying the well established SFR/stellar mass main sequence (MS) are found to dominate both the total IR LF and luminosity density at all redshifts, with the contribution from off MS sources (0.6 dex above MS) being nearly constant (~20% of the total IR luminosity density) and showing no significant signs of increase with increasing z over the whole 0.8<z<2.2 range. Sources with mass in the 10< log(M/Msun) <11 range are found to dominate the total IR LF, with more massive galaxies prevailing at the bright end of the high-z LF. A two-fold evolutionary scheme for IR galaxies is envisaged: on the one hand, a starburst-dominated phase in which the SMBH grows and is obscured by dust, is followed by an AGN dominated phase, then evolving toward a local elliptical. On the other hand, moderately starforming galaxies containing a low-luminosity AGN have various properties suggesting they are good candidates for systems in a transition phase preceding the formation of steady spiral galaxies.

preprint2013arXiv

The rapid assembly of an elliptical galaxy of 400 billion solar masses at a redshift of 2.3

Stellar archeology shows that massive elliptical galaxies today formed rapidly about ten billion years ago with star formation rates above several hundreds solar masses per year (M_sun/yr). Their progenitors are likely the sub-millimeter-bright galaxies (SMGs) at redshifts (z) greater than 2. While SMGs' mean molecular gas mass of 5x10^10 M_sun can explain the formation of typical elliptical galaxies, it is inadequate to form ellipticals that already have stellar masses above 2x10^11 M_sun at z ~ 2. Here we report multi-wavelength high-resolution observations of a rare merger of two massive SMGs at z = 2.3. The system is currently forming stars at a tremendous rate of 2,000 M_sun/yr. With a star formation efficiency an order-of-magnitude greater than that of normal galaxies, it will quench the star formation by exhausting the gas reservoir in only ~200 million years. At a projected separation of 19 kiloparsecs, the two massive starbursts are about to merge and form a passive elliptical galaxy with a stellar mass of ~4x10^11 M_sun. Our observations show that gas-rich major galaxy mergers, concurrent with intense star formation, can form the most massive elliptical galaxies by z ~ 1.5.

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.

preprint2013arXiv

The use and calibration of read-out streaks to increase the dynamic range of the Swift Ultraviolet/Optical Telescope

The dynamic range of photon counting micro-channel-plate (MCP) intensified charged-coupled device (CCD) instruments such as the Swift Ultraviolet/Optical Telescope (UVOT) and the XMM-Newton Optical Monitor (XMM-OM) is limited at the bright end by coincidence loss, the superposition of multiple photons in the individual frames recorded by the CCD. Photons which arrive during the brief period in which the image frame is transferred for read out of the CCD are displaced in the transfer direction in the recorded images. For sufficiently bright sources, these displaced counts form read-out streaks. Using UVOT observations of Tycho-2 stars, we investigate the use of these read-out streaks to obtain photometry for sources which are too bright (and hence have too much coincidence loss) for normal aperture photometry to be reliable. For read-out-streak photometry, the bright-source limiting factor is coincidence loss within the MCPs rather than the CCD. We find that photometric measurements can be obtained for stars up to 2.4 magnitudes brighter than the usual full-frame coincidence-loss limit by using the read-out streaks. The resulting bright-limit Vega magnitudes in the UVOT passbands are UVW2=8.80, UVM2=8.27, UVW1=8.86, u=9.76, b=10.53, v=9.31 and White=11.71; these limits are independent of the windowing mode of the camera. We find that a photometric precision of 0.1 mag can be achieved through read-out streak measurements. A suitable method for the measurement of read-out streaks is described and all necessary calibration factors are given.

preprint2012arXiv

A correlation between intrinsic brightness and average decay rate of Swift UVOT GRB optical/UV light curves

We examine a sample of 48 Swift/UVOT long Gamma-ray Burst light curves and find a correlation between the logarithmic luminosity at 200s and average decay rate determined from 200s onwards, with a Spearman rank coefficient of -0.58 at a significance of 99.998% (4.2 sigma). We discuss the causes of the log L_200s - alpha_>200s correlation, finding it to be an intrinsic property of long GRBs, and not resulting from the selection criteria. We find two ways to produce the correlation. One possibility is that there is some property of the central engine, outflow or external medium that affects the rate of energy release so that the bright afterglows release their energy more quickly and decay faster than the fainter afterglows. Alternatively, the correlation may be produced by variation of the observers viewing angle, with observers at large viewing angles observing fainter and slower decaying light curves.

preprint2012arXiv

A Population of z>2 Far-Infrared Herschel-SPIRE selected Starbursts

We present spectroscopic observations for a sample of 36 Herschel-SPIRE 250-500um selected galaxies (HSGs) at 2<z<5 from the Herschel Multi-tiered Extragalactic Survey (HerMES). Redshifts are confirmed as part of a large redshift survey of Herschel-SPIRE-selected sources covering ~0.93deg^2 in six extragalactic legacy fields. Observations were taken with the Keck I Low Resolution Imaging Spectrometer (LRIS) and the Keck II DEep Imaging Multi-Object Spectrograph (DEIMOS). Precise astrometry, needed for spectroscopic follow-up, is determined by identification of counterparts at 24um or 1.4GHz using a cross-identification likelihood matching method. Individual source luminosities range from log(L_IR/Lsun)=12.5-13.6 (corresponding to star formation rates 500-9000Msun/yr, assuming a Salpeter IMF), constituting some of the most intrinsically luminous, distant infrared galaxies yet discovered. We present both individual and composite rest-frame ultraviolet spectra and infrared spectral energy distributions (SEDs). The selection of these HSGs is reproducible and well characterized across large areas of sky in contrast to most z>2 HyLIRGs in the literature which are detected serendipitously or via tailored surveys searching only for high-z HyLIRGs; therefore, we can place lower limits on the contribution of HSGs to the cosmic star formation rate density at (7+-2)x10^(-3)Msun/yr h^3Mpc^(-3) at z~2.5, which is >10% of the estimated total star formation rate density (SFRD) of the Universe from optical surveys. The contribution at z~4 has a lower limit of 3x10^(-3)Msun/yr h^3 Mpc^(-3), ~>20% of the estimated total SFRD. This highlights the importance of extremely infrared-luminous galaxies with high star formation rates to the build-up of stellar mass, even at the earliest epochs.

preprint2012arXiv

A Redshift Survey of Herschel Far-Infrared Selected Starbursts and Implications for Obscured Star Formation

We present Keck spectroscopic observations and redshifts for a sample of 767 Herschel-SPIRE selected galaxies (HSGs) at 250, 350, and 500um, taken with the Keck I Low Resolution Imaging Spectrometer (LRIS) and the Keck II DEep Imaging Multi-Object Spectrograph (DEIMOS). The redshift distribution of these SPIRE sources from the Herschel Multitiered Extragalactic Survey (HerMES) peaks at z=0.85, with 731 sources at z<2 and a tail of sources out to z~5. We measure more significant disagreement between photometric and spectroscopic redshifts (<delta_z>/(1+z)>=0.29) than is seen in non-infrared selected samples, likely due to enhanced star formation rates and dust obscuration in infrared-selected galaxies. We estimate that the vast majority (72-83%) of z<2 Herschel-selected galaxies would drop out of traditional submillimeter surveys at 0.85-1mm. We estimate the luminosity function and implied star-formation rate density contribution of HSGs at z<1.6 and find overall agreement with work based on 24um extrapolations of the LIRG, ULIRG and total infrared contributions. This work significantly increased the number of spectroscopically confirmed infrared-luminous galaxies at z>>0 and demonstrates the growing importance of dusty starbursts for galaxy evolution studies and the build-up of stellar mass throughout cosmic time. [abridged]

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

HerMES: A Statistical Measurement of the Redshift Distribution of Herschel-SPIRE Sources Using the Cross-correlation Technique

The wide-area imaging surveys with the {\it Herschel} Space Observatory at sub-mm wavelengths have now resulted in catalogs of order one hundred thousand dusty, star-burst galaxies. We make a statistical estimate of $N(z)$ using a clustering analysis of sub-mm galaxies detected at each of 250, 350 and 500 $μ$m from the Herschel Multi-tiered Extragalactic Survey (HerMES) centered on the Boötes field. We cross-correlate {\it Herschel} galaxies against galaxy samples at optical and near-IR wavelengths from the Sloan Digital Sky Survey (SDSS), the NOAO Deep Wide Field Survey (NDWFS) and the Spitzer Deep Wide Field Survey (SDWFS). We create optical and near-IR galaxy samples based on their photometric or spectroscopic redshift distributions and test the accuracy of those redshift distributions with similar galaxy samples defined with catalogs of the Cosmological Evolution Survey (COSMOS), as the COSMOS field has superior spectroscopy coverage. We model-fit the clustering auto and cross-correlations of {\it Herschel} and optical/IR galaxy samples to estimate $N(z)$ and clustering bias factors. The $S_{350} > 20$ mJy galaxies have a bias factor varying with redshift as $b(z)=1.0^{+1.0}_{-0.5}(1+z)^{1.2^{+0.3}_{-0.7}}$. This bias and the redshift dependence is broadly in agreement with galaxies that occupy dark matter halos of mass in the range of 10$^{12}$ to 10$^{13}$ M$_{\sun}$. We find that the redshift distribution peaks around $z \sim 0.5$ to 1 for galaxies selected at 250 $μ$m with an average redshift of $< z > = 1.8 \pm 0.2$. For 350 and 500 $μ$m-selected SPIRE samples the peak shifts to higher redshift, but the average redshift remains the same with a value of $1.9 \pm 0.2$.

preprint2012arXiv

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

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

preprint2012arXiv

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

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

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

preprint2012arXiv

Panchromatic observations of the textbook GRB 110205A: constraining physical mechanisms of prompt emission and afterglow

We present a comprehensive analysis of a bright, long duration (T90 ~ 257 s) GRB 110205A at redshift z= 2.22. The optical prompt emission was detected by Swift/UVOT, ROTSE-IIIb and BOOTES telescopes when the GRB was still radiating in the gamma-ray band. Nearly 200 s of observations were obtained simultaneously from optical, X-ray to gamma-ray, which makes it one of the exceptional cases to study the broadband spectral energy distribution across 6 orders of magnitude in energy during the prompt emission phase. By fitting the time resolved prompt spectra, we clearly identify, for the first time, an interesting two-break energy spectrum, roughly consistent with the standard GRB synchrotron emission model in the fast cooling regime. Although the prompt optical emission is brighter than the extrapolation of the best fit X/gamma-ray spectra, it traces the gamma-ray light curve shape, suggesting a relation to the prompt high energy emission. The synchrotron + SSC scenario is disfavored by the data, but the models invoking a pair of internal shocks or having two emission regions can interpret the data well. Shortly after prompt emission (~ 1100 s), a bright (R = 14.0) optical emission hump with very steep rise (alpha ~ 5.5) was observed which we interpret as the emission from the reverse shock. It is the first time that the rising phase of a reverse shock component has been closely observed. The full optical and X-ray afterglow lightcurves can be interpreted within the standard reverse shock (RS) + forward shock (FS) model. In general, the high quality prompt emission and afterglow data allow us to apply the standard fireball shock model to extract valuable information about the GRB including the radiation mechanism, radius of prompt emission R, initial Lorentz factor of the outflow, the composition of the ejecta, as well as the collimation angle and the total energy budget.

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 Dust & Gas Properties of M83

We examine the dust and gas properties of the nearby, barred galaxy M83, which is part of the Very Nearby Galaxy Survey. Using images from the PACS and SPIRE instruments of Herschel, we examine the dust temperature and dust mass surface density distribution. We find that the nuclear, bar and spiral arm regions exhibit higher dust temperatures and masses compared to interarm regions. However, the distribution of dust temperature and mass are not spatially coincident. Assuming a trailing spiral structure, the dust temperature peaks in the spiral arms lie ahead of the dust surface density peaks. The dust mass surface density correlates well with the distribution of molecular gas as traced by CO (J=3-2) images (JCMT) and the star formation rate as traced by H?2 with a correction for obscured star formation using 24 micron emission. Using HI images from THINGS to trace the atomic gas component, we make total gas mass surface density maps and calculate the gas-to-dust ratio. We find a mean gas-to-dust ratio of 84 \pm 4 with higher values in the inner region assuming a constant CO-to-H2 conversion factor. We also examine the gas-to-dust ratio using CO-to-H2 conversion factor that varies with metallicity.

preprint2012arXiv

The evolutionary connection between QSOs and SMGs: molecular gas in far-infrared luminous QSOs at z ~ 2.5

We present IRAM Plateau de Bure Interferometer observations of the 12CO(3-2) emission from two far-infrared luminous QSOs at z ~ 2.5 selected from the Herschel-ATLAS survey. These far-infrared bright QSOs were selected to have supermassive black holes (SMBHs) with masses similar to those thought to reside in sub-millimetre galaxies (SMGs) at z ~ 2.5; making them ideal candidates as systems in transition from an ultraluminous infrared galaxy phase to a sub-mm faint, unobscured, QSO. We detect 12CO(3-2) emission from both QSOs and we compare their baryonic, dynamical and SMBH masses to those of SMGs at the same epoch. We find that these far-infrared bright QSOs have similar dynamical but lower gas masses than SMGs. In particular we find that far-infrared bright QSOs have ~50+-23% less warm/dense gas than SMGs, which combined with previous results showing the QSOs lack the extended, cool reservoir of gas seen in SMGs, suggests that they are at a different evolutionary stage. This is consistent with the hypothesis that far-infrared bright QSOs represent a short (~1Myr) but ubiquitous phase in the transformation of dust obscured, gas-rich, starburst-dominated SMGs into unobscured, gas-poor, QSOs.

preprint2012arXiv

The gas-to-dust mass ratio of Centaurus A as seen by Herschel

We present photometry of the nearby galaxy NGC 5128 (Centaurus A) observed with the PACS and SPIRE instruments on board the Herschel Space Observatory, at 70, 160, 250, 350 and 500 μm, as well as new CO J = 3-2 observations taken with the HARP-B instrument on the JCMT. Using a single component modified blackbody, we model the dust spectral energy distribution within the disk of the galaxy using all five Herschel wavebands, and find dust temperatures of ~30 K towards the centre of the disk and a smoothly decreasing trend to ~20 K with increasing radius. We find a total dust mass of (1.59 \pm 0.05) \times 10^7 M\odot, and a total gas mass of (2.7 \pm 0.2) \times 10^9 M\odot. The average gas-to-dust mass ratio is 103 \pm 8 but we find an interesting increase in this ratio to approximately 275 toward the centre of Cen A. We discuss several possible physical processes that may be causing this effect, including dust sputtering, jet entrainment and systematic variables such as the XCO factor. Dust sputtering by X-rays originating in the AGN or the removal of dust by the jets are our most favoured explanations.

preprint2012arXiv

The Herschel Exploitation of Local Galaxy Andromeda (HELGA) II: Dust and Gas in Andromeda

We present an analysis of the dust and gas in Andromeda, using Herschel images sampling the entire far-infrared peak. We fit a modified-blackbody model to ~4000 quasi-independent pixels with spatial resolution of ~140pc and find that a variable dust-emissivity index (beta) is required to fit the data. We find no significant long-wavelength excess above this model suggesting there is no cold dust component. We show that the gas-to-dust ratio varies radially, increasing from ~20 in the center to ~70 in the star-forming ring at 10kpc, consistent with the metallicity gradient. In the 10kpc ring the average beta is ~1.9, in good agreement with values determined for the Milky Way (MW). However, in contrast to the MW, we find significant radial variations in beta, which increases from 1.9 at 10kpc to ~2.5 at a radius of 3.1kpc and then decreases to 1.7 in the center. The dust temperature is fairly constant in the 10kpc ring (ranging from 17-20K), but increases strongly in the bulge to ~30K. Within 3.1kpc we find the dust temperature is highly correlated with the 3.6 micron flux, suggesting the general stellar population in the bulge is the dominant source of dust heating there. At larger radii, there is a weak correlation between the star formation rate and dust temperature. We find no evidence for 'dark gas' in M31 in contrast to recent results for the MW. Finally, we obtained an estimate of the CO X-factor by minimising the dispersion in the gas-to-dust ratio, obtaining a value of (1.9+/-0.4)x10^20 cm^-2 [K kms^-1]^-1.

preprint2012arXiv

The Herschel Exploitation of Local Galaxy Andromeda (HELGA). I: Global far-infrared and sub-mm morphology

We have obtained Herschel images at five wavelengths from 100 to 500 micron of a ~5.5x2.5 degree area centred on the local galaxy M31 (Andromeda), our nearest neighbour spiral galaxy, as part of the Herschel guaranteed time project "HELGA". The main goals of HELGA are to study the characteristics of the extended dust emission, focusing on larger scales than studied in previous observations of Andromeda at an increased spatial resolution, and the obscured star formation. In this paper we present data reduction and Herschel maps, and provide a description of the far-infrared morphology, comparing it with features seen at other wavelengths. We use high--resolution maps of the atomic hydrogen, fully covering our fields, to identify dust emission features that can be associated to M31 with confidence, distinguishing them from emission coming from the foreground Galactic cirrus. Thanks to the very large extension of our maps we detect, for the first time at far-infrared wavelengths, three arc-like structures extending out to ~21, ~26 and ~31 kpc respectively, in the south-western part of M31. The presence of these features, hosting ~2.2e6 Msol of dust, is safely confirmed by their detection in HI maps. Overall, we estimate a total dust mass of ~5.8e7 Msol, about 78% of which is contained in the two main ring-like structures at 10 and 15 kpc, at an average temperature of 16.5 K. We find that the gas-to-dust ratio declines exponentially as a function of the galacto-centric distance, in agreement with the known metallicity gradient, with values ranging from 66 in the nucleus to ~275 in the outermost region. [Abridged]

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

The SCUBA-2 Cosmology Legacy Survey: blank-field number counts of 450um-selected galaxies and their contribution to the cosmic infrared background

The first deep blank-field 450um map (1-sigma~1.3mJy) from the SCUBA-2 Cosmology Legacy Survey (S2CLS), conducted with the James Clerk Maxwell Telescope (JCMT) is presented. Our map covers 140 arcmin^2 of the COSMOS field, in the footprint of the HST CANDELS area. Using 60 submillimetre galaxies (SMGs) detected at >3.75-sigma, we evaluate the number counts of 450um-selected galaxies with flux densities S_450>5mJy. The 8-arcsec JCMT beam and high sensitivity of SCUBA-2 now make it possible to directly resolve a larger fraction of the cosmic infrared background (CIB, peaking at ~200um) into the individual galaxies responsible for its emission than has previously been possible at this wavelength. At S_450>5mJy we resolve (7.4[+/-]0.7)x10^-2 MJy/sr of the CIB at 450um (equivalent to 16[+/-]7% of the absolute brightness measured by COBE at this wavelength) into point sources. A further ~40% of the CIB can be recovered through a statistical stack of 24um emitters in this field, indicating that the majority (~60%) of the CIB at 450um is emitted by galaxies with S_450>2mJy. The average redshift of 450um emitters identified with an optical/near-infrared counterpart is estimated to be <z>=1.3, implying that the galaxies in the sample are in the ultraluminous class (L_IR~1.1x10^12 L_sun). If the galaxies contributing to the statistical stack lie at similar redshifts, then the majority of the CIB at 450um is emitted by galaxies in the LIRG class with L_IR>3.6x10^11 L_sun.

preprint2012arXiv

The XMM-Newton serendipitous ultraviolet source survey catalogue

The XMM-Newton Serendipitous Ultraviolet Source Survey (XMM-SUSS) is a catalogue of ultraviolet (UV) sources detected serendipitously by the Optical Monitor (XMM-OM) on-board the XMM-Newton observatory. The catalogue contains ultraviolet-detected sources collected from 2,417 XMM-OM observations in 1-6 broad band UV and optical filters, made between 24 February 2000 and 29 March 2007. The primary contents of the catalogue are source positions, magnitudes and fluxes in 1 to 6 passbands, and these are accompanied by profile diagnostics and variability statistics. The XMM-SUSS is populated by 753,578 UV source detections above a 3 sigma signal-to-noise threshold limit which relate to 624,049 unique objects. Taking account of substantial overlaps between observations, the net sky area covered is 29-54 square degrees, depending on UV filter. The magnitude distributions peak at 20.2, 20.9 and 21.2 in UVW2, UVM2 and UVW1 respectively. More than 10 per cent of sources have been visited more than once using the same filter during XMM-Newton operation, and > 20 per cent of sources are observed more than once per filter during an individual visit. Consequently, the scope for science based on temporal source variability on timescales of hours to years is broad. By comparison with other astrophysical catalogues we test the accuracy of the source measurements and define the nature of the serendipitous UV XMM-OM source sample. The distributions of source colours in the UV and optical filters are shown together with the expected loci of stars and galaxies, and indicate that sources which are detected in multiple UV bands are predominantly star-forming galaxies and stars of type G or earlier.

preprint2011arXiv

A Multi-wavelength MOCASSIN Model of the Magellanic-type Galaxy NGC 4449

We use the photoionisation and dust radiative transfer code MOCASSIN to create a model of the dwarf irregular galaxy NGC 4449. The best-matching model reproduces the global optical emission line fluxes and the observed spectral energy distribution (SED) spanning wavelengths from the UV to sub-mm, and requires the bolometric luminosity of 6.25e9 Lsolar for the underlying stellar component, M_d/M_g of 1/680 and M_d of 2.2e6 Msolar.

preprint2011arXiv

A strongly starforming group: three massive galaxies associated with a QSO

We present here photometric redshift confirmation of the presence of large scale structure around the z=1.82 QSO RXJ0941, which shows an overdensity of submm sources. Radio imaging confirms the presence of the submm sources and pinpoints their likely optical-NIR counterparts. Four of the five submm sources present in this field (including the QSO) have counterparts with redshifts compatible with z=1.82. We show that our photometric redshifts are robust against the use of different spectral templates. We have measured the galaxy stellar mass of the submm galaxies from their rest-frame K-band luminosity obtaining log(M*/Msun)~11.5+-0.2, slightly larger than the Schechter mass of present day galaxies, and hence indicating that most of the stellar mass is already formed. We present optical-to-radio spectral energy distributions (SEDs) of the five SCUBA sources. The emission of RXJ0941 is dominated by reprocessed AGN emission in the observed MIR range, while the starburst contribution completely dominates in the submm range. The SEDs of the other three counterparts are compatible with a dominant starburst contribution above ~24um, with star formation rates SFR~2000Msun/yr, central dust masses log(Mdust/Msun)~9+-0.5 and hence central gas masses log(Mgas/Msun)~10.7. There is very little room for an AGN contribution. From X-ray upper limits and the observed 24um flux, we derive a maximum 2-10keV X-ray luminosity of 1e44 erg/s for any putative AGN, even if they are heavily obscured. This in turn points to relatively small black holes with log(MBH/Msun)<~8 and hence stellar-to-black hole mass ratios about one order of magnitude higher than those observed in the present Universe: most of their central black hole masses are still to be accreted. Local stellar-to-black hole masses ratios can be reached if ~1.3% of the available nuclear gas mass is accreted.

preprint2011arXiv

An Updated Ultraviolet Calibration for the Swift/UVOT

We present an updated calibration of the Swift/UVOT broadband ultraviolet (uvw1, uvm2, and uvw2) filters. The new calibration accounts for the ~1% per year decline in the UVOT sensitivity observed in all filters, and makes use of additional calibration sources with a wider range of colours and with HST spectrophotometry. In this paper we present the new effective area curves and instrumental photometric zeropoints and compare with the previous calibration.

preprint2011arXiv

Averaging the AGN X-ray spectra from deep Chandra fields

The X-ray spectra of Active Galactic Nuclei (AGN) carry the signatures of the emission from the central region, close to the Super Massive Black Hole (SMBH). For this reason, the study of deep X-ray spectra is a powerful instrument to investigate the origin of their emission. The emission line most often observed in the X-ray spectra of AGN is Fe K. It is known that it can be broadened and deformed by relativistic effects if emitted close enough to the central SMBH. In recent statistical studies of the X-ray spectra of AGN samples, it is found that a narrow Fe line is ubiquitous, while whether the broad features are as common is still uncertain. We present here the results of an investigation on the characteristics of the Fe line in the average X-ray spectra of AGN in deep Chandra fields. The average spectrum of the AGN is computed using Chandra spectra with more than 200 net counts from the AEGIS, Chandra Deep Field North and Chandra Deep Field South surveys. The sample spans a broader range of X-ray luminosities than other samples studied with stacking methods up to z=3.5. We analyze the average spectra of this sample using our own averaging method, checking the results against extensive simulations. Subsamples defined in terms of column density of the local absorber, luminosity and z are also investigated. We found a very significant Fe line with a narrow profile in all our samples and in almost all the subsamples that we constructed. The equivalent width (EW) of the narrow line estimated in the average spectrum of the full sample is 74 eV. The broad line component is significantly detected in the subsample of AGN with L<1.43 1E44 cgs and z<0.76, with EW=108 eV. We concluded that the narrow Fe line is an ubiquitous feature of the X-ray spectra of the AGN up to z=3.5.The broad line component is significant in the X-ray spectra of the AGN with low luminosity and low z.

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

Dynamical Structure of the Molecular Interstellar Medium in an Extremely Bright, Multiply Lensed z~3 Submillimeter Galaxy Discovered with Herschel

We report the detection of CO 5-4, 3-2, and 1-0 emission in the strongly lensed, Herschel/SPIRE-selected submillimeter galaxy (SMG) HLSW-01 at z=2.9574+/-0.0001, using the Plateau de Bure Interferometer, the Combined Array for Research in Millimeter-wave Astronomy, and the Green Bank Telescope. The observations spatially resolve the molecular gas into four lensed images with a maximum separation of ~9", and reveal the internal gas dynamics in this system. We derive lensing-corrected CO line luminosities of L'(CO 1-0) = (4.17+/-0.41), L'(CO 3-2) = (3.96+/-0.20) and L'(CO 5-4) = (3.45+/-0.20) x 10^10 (mu_L/10.9)^-1 Kkm/s pc^2, corresponding to luminosity ratios of r_31 = 0.95+/-0.10, r_53 = 0.87+/-0.06, and r_51 = 0.83+/-0.09. This suggests a total molecular gas mass of Mgas = 3.3 x 10^10 (alpha_CO/0.8) (mu_L/10.9)^-1 Msun. The gas mass, gas mass fraction, gas depletion timescale, star formation efficiency, and specific star formation rate are typical for an SMG. The velocity structure of the gas reservoir suggests that the brightest two lensed images are dynamically resolved projections of the same dust-obscured region in the galaxy that are kinematically offset from the unresolved fainter images. The resolved kinematics appear consistent with the complex velocity structure observed in major, `wet' (i.e., gas-rich) mergers. Major mergers are commonly observed in SMGs, and are likely to be responsible for fueling their intense starbursts at high gas consumption rates. This study demonstrates the level of detail to which galaxies in the early universe can be studied by utilizing the increase in effective spatial resolution and sensitivity provided by gravitational lensing.

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: Lyman Break Galaxies individually detected at 0.7 < z < 2.0 in GOODS-N with Herschel/SPIRE

As part of the Herschel Multi-tiered Extragalactic Survey we have investigated the rest-frame far-infrared (FIR) properties of a sample of more than 4800 Lyman Break Galaxies (LBGs) in the Great Observatories Origins Deep Survey North field. Most LBGs are not detected individually, but we do detect a sub-sample of 12 objects at 0.7 < z < 1.6 and one object at z ~ 2.0. The ones detected by Herschel SPIRE have redder observed NUV-U and U-R colors than the others, while the undetected ones have colors consistent with average LBGs at z > 2.5. The UV-to-FIR spectral energy distributions of the objects detected in the rest-frame FIR are investigated using the code CIGALE to estimate physical parameters. We find that LBGs detected by SPIRE are high mass, luminous infrared galaxies. It appears that LBGs are located in a triangle-shaped region in the A_FUV vs. Log L_FUV diagram limited by A_FUV=0 at the bottom and by a diagonal following the temporal evolution of the most massive galaxies from the bottom-right to the top-left of the diagram. This upper envelop can be used as upper limits for the UV dust attenuation as a function of L_FUV}. The limits of this region are well explained using a closed-box model, where the chemical evolution of galaxies produces metals, which in turn lead to higher dust attenuation when the galaxies age.

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-ATLAS: The link between accretion luminosity and star formation in quasar host galaxies

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

preprint2011arXiv

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

Investigations of dust heating in M81, M83, and NGC 2403 with the Herschel Space Observatory

We use Spitzer Space Telescope and Herschel Space Observatory far-infrared data along with ground-based optical and near-infrared data to understand how dust heating in the nearby face-on spiral galaxies M81, M83, and NGC 2403 is affected by the starlight from all stars and by the radiation from star forming regions. We find that 70/160 micron surface brightness ratios tend to be more strongly influenced by star forming regions. However, the 250/350 micron and 350/500 micron surface brightness ratios are more strongly affected by the light from the total stellar populations, suggesting that the dust emission at >250 microns originates predominantly from a component that is colder than the dust seen at <160 microns and that is relatively unaffected by star formation activity. We conclude by discussing the implications of this for modelling the spectral energy distributions of both nearby and more distant galaxies and for using far-infrared dust emission to trace star formation.

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

Redshift Determination and CO Line Excitation Modeling for the Multiply-Lensed Galaxy HLSW-01

We report on the redshift measurement and CO line excitation of HERMES J105751.1+573027 (HLSW-01), a strongly lensed submillimeter galaxy discovered in Herschel/SPIRE observations as part of the Herschel Multi-tiered Extragalactic Survey (HerMES). HLSW-01 is an ultra-luminous galaxy with an intrinsic far-infrared luminosity of 1.4x10^(13) solar luminosities, and is lensed by a massive group of galaxies into at least four images with a total magnification of 10.9+/-0.7. With the 100 GHz instantaneous bandwidth of the Z-Spec instrument on the Caltech Submillimeter Observatory, we robustly identify a redshift of z=2.958+/-0.007 for this source, using the simultaneous detection of four CO emission lines (J = 7-6, J = 8-7, J = 9-8, and J = 10-9). Combining the measured line fluxes for these high-J transitions with the J = 1-0, J = 3-2 and J = 5-4 line fluxes measured with the Green Bank Telescope, the Combined Array for Research in Millimeter Astronomy, and the Plateau de Bure Interferometer, respectively, we model the physical properties of the molecular gas in this galaxy. We find that the full CO spectral line energy distribution is well described by warm, moderate-density gas with Tkin = 86-235 K and n(H2) = (1.1-3.5)x10^3 cm^(-3). However, it is possible that the highest-J transitions are tracing a small fraction of very dense gas in molecular cloud cores, and two-component models that include a warm/dense molecular gas phase with Tkin ~ 200 K, n(H2) ~ 10^5 cm^(-3) are also consistent with these data. Higher signal-to-noise measurements of the J(upper) > 7 transitions with high spectral resolution, combined with high spatial resolution CO maps, are needed to improve our understanding of the gas excitation, morphology, and dynamics of this interesting high-redshift galaxy.

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

Swift UV/Optical Telescope Imaging of Star Forming Regions in M81 and Holmberg IX

We present Swift UV/Optical Telescope (UVOT) imaging of the galaxies M81 and Holmberg IX. We combine UVOT imaging in three near ultraviolet (NUV) filters (uvw2: 1928 Å, uvm2: 2246 Å, and uvw1: 2600 Å) with ground based optical imaging from the Sloan Digital Sky Survey to constrain the stellar populations of both galaxies. Our analysis consists of three different methods. First we use the NUV imaging to identify UV star forming knots and then perform SED modeling on the UV/optical photometry of these sources. Second, we measure surface brightness profiles of the disk of M81 in the NUV and optical. Last we use SED fitting of individual pixels to map the properties of the two galaxies. In agreement with earlier studies we find evidence for a burst in star formation in both galaxies starting ~200 Myr ago coincident with the suggested time of an M81-M82 interaction. In line with theories of its origin as a tidal dwarf we find that the luminosity weighted age of Holmberg IX is a few hundred million years. Both galaxies are best fit by a Milky Way dust extinction law with a prominent 2175 Å bump. In addition, we describe a stacked median filter technique for modeling the diffuse background light within a galaxy, and a Markov chain method for cleaning segment maps generated by SExtractor.

preprint2011arXiv

The Dust Extinction Curves of Gamma-Ray Burst Host Galaxies

The composition and amount of interstellar dust within gamma-ray burst (GRB) host galaxies is of key importance when addressing selection effects in the GRB redshift distribution, and when studying the properties of their host galaxies. As well as the implications for GRB research, probing the dust within the high-z hosts of GRBs also contributes to our understanding of the conditions of the interstellar medium and star-formation in the distant Universe. Nevertheless, the physical properties of dust within GRB host galaxies continues to be a highly contended issue. In this paper we explore the mean extinction properties of dust within the host galaxies of a sample of 17 GRBs with total host galaxy visual extinction Av<1 (<Av>=0.4), covering a redshift range z=0.7-3.1. We find the average host extinction curve to have an ultraviolet slope comparable to that of the LMC, but with little evidence of a 2175Angs dust extinction feature as observed along Milky Way and LMC sightlines. We cannot at present rule out the presence of a 2175Angs feature, and both the standard SMC and LMC extinction curves also provide good fits to our data. However, we can reject an extinction curve that has a UV slope as flat as the mean Milky Way extinction curve, whilst also having a 2175Angs feature as prominent as seen in the mean Milky Way extinction curve. This is in contrast to the clear detection of a 2175Angs bump and the flatter extinction curves of some more heavily extinguished GRBs (Av>1), which may be indicative of there being a dependence between dust abundance and the wavelength dependence of dust extinction, as has been previously speculated.

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.

preprint2011arXiv

The nature of X-ray absorbed QSOs

There exists a significant population of broad line, z~2 QSOs which have heavily absorbed X-ray spectra. Follow up observations in the submillimetre show that these QSOs are embedded in ultraluminous starburst galaxies, unlike most unabsorbed QSOs at the same redshifts and luminosities. Here we present X-ray spectra from XMM-Newton for a sample of 5 such X-ray absorbed QSOs that have been detected at submillimetre wavelengths. We also present spectra in the restframe ultraviolet from ground based telescopes. All 5 QSOs are found to exhibit strong C IV absorption lines in their ultraviolet spectra with equivalent width > 5 Angstroms. The X-ray spectra are inconsistent with the hypothesis that these objects show normal QSO continua absorbed by low-ionization gas. Instead, the spectra can be modelled successfully with ionized absorbers, or with cold absorbers if they posess unusually flat X-ray continuum shapes and unusual optical to X-ray spectral energy distributions. We show that the ionized absorber model provides the simplest, most self-consistent explanation for their observed properties. We estimate that the fraction of radiated power that is converted into kinetic luminosity of the outflowing winds is typically ~4 per cent, in agreement with recent estimates for the kinetic feedback from QSOs required to produce the M - sigma relation, and consistent with the hypothesis that the X-ray absorbed QSOs represent the transition phase between obscured accretion and the luminous QSO phase in the evolution of massive galaxies.

preprint2011arXiv

The X-ray spectral properties of the AGN population in the XMM-Newton bright serendipitous survey

We present here a detailed X-ray spectral analysis of the AGN belonging to the XMM-Newton bright survey (XBS) that comprises more than 300 AGN up to redshift ~ 2.4. We performed an X-ray analysis following two different approaches: by analyzing individually each AGN X-ray spectrum and by constructing average spectra for different AGN types. From the individual analysis, we find that there seems to be an anti correlation between the spectral index and the sources' hard X-ray luminosity, such that the average photon index for the higher luminosity sources (> 10E44 erg/s) is significantly flatter than the average for the lower luminosity sources. We also find that the intrinsic column density distribution agrees with AGN unified schemes, although a number of exceptions are found (3% of the whole sample), which are much more common among optically classified type 2 AGN. We also find that the so-called "soft-excess", apart from the intrinsic absorption, constitutes the principal deviation from a power-law shape in AGN X-ray spectra and it clearly displays different characteristics, and likely a different origin, for unabsorbed and absorbed AGN. Regarding the shape of the average spectra, we find that it is best reproduced by a combination of an unabsorbed (absorbed) power law, a narrow Fe Kalpha emission line and a small (large) amount of reflection for unabsorbed (absorbed) sources. We do not significantly detect any relativistic contribution to the line emission and we compute an upper limit for its equivalent width (EW) of 230 eV at the 3 sigma confidence level. Finally, by dividing the type 1 AGN sample into high- and low-luminosity sources, we marginally detect a decrease in the narrow Fe Kalpha line EW and in the amount of reflection as the luminosity increases, the "so-called" Iwasawa-Taniguchi effect.

preprint2010arXiv

A Statistical Comparison of the Optical/UV and X-ray Afterglows of Gamma-Ray Bursts using the Swift Ultra-violet Optical and X-ray Telescopes

We present the systematic analysis of the UVOT and XRT light curves for a sample of 26 Swift Gamma-Ray Bursts (GRBs). By comparing the optical/UV and X-ray light curves, we found that they are remarkably different during the first 500s after the BAT trigger, while they become more similar during the middle phase of the afterglow, i.e. between 2000s and 20000s. If we take literally the average properties of the sample, we find that the mean temporal indices observed in the optical/UV and X-rays after 500s are consistent with a forward-shock scenario, under the assumptions that electrons are in the slow cooling regime, the external medium is of constant density and the synchrotron cooling frequency is situated between the optical/UV and X-ray observing bands. While this scenario describes well the averaged observed properties, some individual GRB afterglows require different or additional assumptions, such as the presence of late energy injection. We show that a chromatic break (a break in the X-ray light curve that is not seen in the optical) is present in the afterglows of 3 GRBs and demonstrate evidence for chromatic breaks in a further 4 GRBs. The average properties of these breaks cannot be explained in terms of the passage of the synchrotron cooling frequency through the observed bands, nor a simple change in the external density. It is difficult to reconcile chromatic breaks in terms of a single component outflow and instead, more complex jet structure or additional emission components are required.

preprint2010arXiv

An excess of star-forming galaxies in the fields of high-redshift QSOs

We present submillimetre and mid-infrared imaging observations of five fields centred on quasi-stellar objects (QSOs) at 1.7<z<2.8. All 5 QSOs were detected previously at submillimetre wavelengths. At 850 (450) um we detect 17 (11) submillimetre galaxies (SMGs) in addition to the QSOs. The total area mapped at 850 um is ~28 arcmin^2 down to RMS noise levels of 1-2 mJy/beam, depending on the field. Integral number counts are computed from the 850 um data using the same analytical techniques adopted by `blank-field' submillimetre surveys. We find that the `QSO-field' counts show a clear excess over the blank-field counts at deboosted flux densities of 2-4 mJy; at higher flux densities the counts are consistent with the blank-field counts. Robust mid-infrared counterparts are identified for all four submillimetre detected QSOs and ~60 per cent of the SMGs. The mid-infrared colours of the QSOs are similar to those of the local ULIRG/AGN Mrk 231 if placed at 1<z<3 whilst most of the SMGs have colours very similar to those of the local ULIRG Arp 220 at 1<z<3. Mid-infrared diagnostics therefore find no strong evidence that the SMGs host buried AGN although we cannot rule out such a possibility. Taken together our results suggest that the QSOs sit in regions of the early universe which are undergoing an enhanced level of major star-formation activity, and should evolve to become similarly dense regions containing massive galaxies at the present epoch. Finally, we find evidence that the level of star-formation activity in individual galaxies appears to be lower around the QSOs than it is around more powerful radio-loud AGN at higher redshifts.

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 : SPIRE detection of high redshift massive compact galaxies in GOODS-N field

We have analysed the rest-frame far infrared (FIR) properties of a sample of massive (Mstar > 10^11Msun) galaxies at 2<z<3 in the GOODS (Great Observatories Origins Deep Survey) North field using the Spectral and Photometric Imaging Receiver (SPIRE) instrument aboard the Herschel Space Observatory. To conduct this analysis we take advantage of the data from the HerMES key program. The sample comprises 45 massive galaxies with structural parameters characterised with HST NICMOS-3. We study detections at submm Herschel bands, together with Spitzer 24μm data, as a function of the morphological type, mass and size. We find that 26/45 sources are detected at MIPS-24μm and 15/45 (all MIPS-24μm detections) are detected at SPIRE-250μm, with disk-like galaxies more easily detected. We derive star formation rates (SFR) and specific star formation rates (sSFR) by fitting the spectral energy distribution (SED) of our sources, taking into account non-detections for SPIRE and systematic effects for MIPS derived quantities. We find that the mean SFR for the spheroidal galaxies (50-100 Msun*yr^-1) is substantially (a factor ~ 3) lower than the mean value presented by disk-like galaxies (250-300 Msun*yr^-1).

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: Herschel-SPIRE observations of Lyman Break Galaxies

We present first results of a study of the submillimetre (rest frame far-infrared) properties of z~3 Lyman Break Galaxies (LBGs) and their lower-redshift counterparts BX/BM galaxies, based on Herschel-SPIRE observations of the Northern field of the Great Observatories Origins Deep Survey (GOODS-N). We use stacking analysis to determine the properties of LBGs well below the current limit of the survey. Although LBGs are not detected individually, stacking the infrared luminous LBGs (those detected with Spitzer at 24 microns yields a statistically significant submm detection with mean flux <S_{250}>= 5.9+/-1.4 mJy confirming the power of SPIRE in detecting UV-selected high-redshift galaxies at submillimetre wavelengths. In comparison, the Spitzer 24 microns detected BX/BM galaxies appear fainter with a stacked value of <S_{250}> = 2.7 +/-0.8 mJy. By fitting the Spectral Energy Distributions (SEDs) we derive median infrared luminosities, L_{IR}, of 2.8x10^{12} Lsun and 1.5x10^{11} Lsun for z~3 LBGs and BX/BMs, respectively. We find that $L_{IR} estimates derived from present measurements are in good agreement with those based on UV data for z~2 BX/BM galaxies, unlike the case for z~3 infrared luminous LBGs where the UV underestimates the true $L_{IR}. Although sample selection effects may influence this result we suggest that differences in physical properties (such as morphologies, dust distribution and extent of star-forming regions) between z ~3 LBGs and z~2 BX/BMs may also play a significant role.

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 and SCUBA-2 imaging and spectroscopy of a bright, lensed submillimetre galaxy at z = 2.3

We present a detailed analysis of the far-IR properties of the bright, lensed, z = 2.3, SMG, SMM J2135-0102, using new observations with Herschel, SCUBA-2 and the VLA. These data allow us to constrain the galaxy's SED and show that it has an intrinsic rest-frame 8-1000um luminosity, L(bol), of (2.3 +/- 0.2) x 10^12 L(sun) and a likely SFR of ~400 M(sun)/yr. The galaxy sits on the far-IR/radio correlation for far-IR-selected galaxies. At ~>70um, the SED can be described adequately by dust components with T(d) ~ 30 and 60K. Using SPIRE's Fourier Transform Spectrometer we report a detection of the [CII] 158um cooling line. If the [CII], CO and far-IR continuum arise in photo-dissociation regions, we derive a characteristic gas density, n ~ 10^3 cm^-3, and a far-UV radiation field, G_0, 10^3x stronger than the Milky Way. L([CII])/L(bol) is significantly higher than in local ULIRGs but similar to the values found in local star-forming galaxies and starburst nuclei. This is consistent with SMM J2135-0102 being powered by starburst clumps distributed across ~2 kpc, evidence that SMGs are not simply scaled-up ULIRGs. Our results show that SPIRE's FTS has the ability to measure the redshifts of distant, obscured galaxies via the blind detection of atomic cooling lines, but it will not be competitive with ground-based CO-line searches. It will, however, allow detailed study of the integrated properties of high-redshift galaxies, as well as the chemistry of their ISM, once more suitably bright candidates have been found.

preprint2010arXiv

Herschel photometric observations of the low metallicity dwarf galaxy NGC 1705

We present Herschel SPIRE and PACS photometeric observations of the low metallicity (Z ~ 0.35 solar) nearby dwarf galaxy, NGC 1705, in six wavelength bands as part of the Dwarf Galaxy Survey guaranteed time Herschel Key Program. We confirm the presence of two dominant circumnuclear IR-bright regions surrounding the central super star cluster that had been previously noted at mid-IR wavelengths and in the sub-mm by LABOCA. On constructing a global spectral energy distribution using the SPIRE and PACS photometry, in conjunction with archival IR measurements, we note the presence of an excess at sub-mm wavelengths. This excess suggests the presence of a significant cold dust component within NGC 1705 and was modeled as an additional cold component in the SED. Although alternative explanations for the sub-mm excess beyond 350 microns, such as changes to the dust emissivity cannot be ruled out, the most likely explanation for the observed submillimetre excess is that of an additional cold dust component.

preprint2010arXiv

Herschel photometric observations of the nearby low metallicity irregular galaxy NGC 6822

We present the first Herschel PACS and SPIRE images of the low-metallicity galaxy NGC6822 observed from 70 to 500 mu and clearly resolve the HII regions with PACS and SPIRE. We find that the ratio 250/500 is dependent on the 24 mu surface brightness in NGC6822, which would locally link the heating processes of the coldest phases of dust in the ISM to the star formation activity. We model the SEDs of some regions HII regions and less active regions across the galaxy and find that the SEDs of HII regions show warmer ranges of dust temperatures. We derive very high dust masses when graphite is used in our model to describe carbon dust. Using amorphous carbon, instead, requires less dust mass to account for submm emission due to its lower emissivity properties. This indicates that SED models including Herschel constraints may require different dust properties than commonly used.

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 observations of the disturbed galaxy NGC4438

We present Herschel-SPIRE observations of the perturbed galaxy NGC4438 in the Virgo cluster. These images reveal the presence of extra-planar dust up to ~4-5 kpc away from the galaxy's disk. The dust closely follows the distribution of the stripped atomic and molecular hydrogen, supporting the idea that gas and dust are perturbed in a similar fashion by the cluster environment. Interestingly, the extra-planar dust lacks a warm temperature component when compared to the material still present in the disk, explaining why it was missed by previous far-infrared investigations. Our study provides evidence for dust stripping in clusters of galaxies and illustrates the potential of Herschel data for our understanding of environmental effects on galaxy evolution.

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

Mapping the interstellar medium in galaxies with Herschel/SPIRE

The standard method of mapping the interstellar medium in a galaxy, by observing the molecular gas in the CO 1-0 line and the atomic gas in the 21-cm line, is largely limited with current telescopes to galaxies in the nearby universe. In this letter, we use SPIRE observations of the galaxies M99 and M100 to explore the alternative approach of mapping the interstellar medium using the continuum emission from the dust. We have compared the methods by measuring the relationship between the star-formation rate and the surface density of gas in the galaxies. We find the two methods give relationships with a similar dispersion, confirming that observing the continuum emission from the dust is a promising method of mapping the interstellar medium in galaxies.

preprint2010arXiv

Mid-infrared spectroscopy of candidate AGN-dominated submillimeter galaxies

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

preprint2010arXiv

On the electron energy distribution index of Swift Gamma-Ray Burst Afterglows

The electron energy distribution index, p, is a fundamental parameter of the synchrotron emission from a range of astronomical sources. Here we examine one such source of synchrotron emission, Gamma-Ray Burst afterglows observed by the Swift satellite. Within the framework of the blast wave model, we examine the constraints placed on the distribution of p by the observed X-ray spectral indices and parametrise the distribution. We find that the observed distribution of spectral indices are inconsistent with an underlying distribution of p composed of a single discrete value but consistent with a Gaussian distribution centred at p = 2.36 and having a width of 0.59. Furthermore, accepting that the underlying distribution is a Gaussian, we find the majority (>94%) of GRB afterglows in our sample have cooling break frequencies less than the X-ray frequency.

preprint2010arXiv

On the origin of M81 group extended dust emission

Galactic cirrus emission at far-infrared wavelengths affects many extragalactic observations. Separating this emission from that associated with extragalactic objects is both important and difficult. In this paper we discuss a particular case, the M81 group, and the identification of diffuse structures prominent in the infrared, but also detected at optical wavelengths. The origin of these structures has previously been controversial, ranging from them being the result of a past interaction between M81 and M82 or due to more local Galactic emission. We show that over of order a few arcminute scales the far-infrared (Herschel 250 &μ&m) emission correlates spatially very well with a particular narrow velocity (2-3 km/s) component of the Galactic HI. We find no evidence that any of the far-infrared emission associated with these features actually originates in the M81 group. Thus we infer that the associated diffuse optical emission must be due to galactic light back scattered off dust in our galaxy. Ultra-violet observations pick out young stellar associations around M81, but no detectable far-infrared emission. We consider in detail one of the Galactic cirrus features, finding that the far-infrared HI relation breaks down below arc minute scales and that at smaller scales there can be quite large dust temperature variations.

preprint2010arXiv

Paper II: Calibration of the Swift ultraviolet/optical telescope

The Ultraviolet/Optical Telescope (UVOT) is one of three instruments onboard the Swift observatory. The photometric calibration has been published, and this paper follows up with details on other aspects of the calibration including a measurement of the point spread function with an assessment of the orbital variation and the effect on photometry. A correction for large scale variations in sensitivity over the field of view is described, as well as a model of the coincidence loss which is used to assess the coincidence correction in extended regions. We have provided a correction for the detector distortion and measured the resulting internal astrometric accuracy of the UVOT, also giving the absolute accuracy with respect to the International Celestial Reference System. We have compiled statistics on the background count rates, and discuss the sources of the background, including instrumental scattered light. In each case we describe any impact on UVOT measurements, whether any correction is applied in the standard pipeline data processing or whether further steps are recommended.

preprint2010arXiv

Photometric redshifts for GRB afterglows from GROND and Swift/UVOT

We present a framework to obtain photometric redshifts (photo-zs) for gamma-ray burst afterglows. Using multi-band photometry from GROND and Swift/UVOT, photo-zs are derived for five GRBs for which spectroscopic redshifts are not available. We use UV/optical/NIR data and synthetic photometry based on afterglow observations and theory to derive the photometric redshifts of GRBs and their accuracy. Taking into account the afterglow synchrotron emission properties, we investigate the application of photometry to derive redshifts in a theoretical range between z~1 and z~12. The measurement of photo-zs for GRB afterglows provides a quick, robust and reliable determination of the distance scale to the burst, particularly in those cases where spectroscopic observations in the optical/NIR range cannot be obtained. Given a sufficiently bright and mildly reddened afterglow, the relative photo-z accuracy is better than 10% between z=1.5 and z~7 and better than 5% between z=2 and z=6. We detail the approach on 5 sources without spectroscopic redshifts observed with UVOT on-board Swift and/or GROND. The distance scale to those same afterglows is measured to be $z=4.31^{+0.14}_{-0.15}$ for GRB 080825B, $z=2.13^{+0.14}_{-0.20}$ for GRB 080906, $z=3.44^{+0.15}_{-0.32}$ for GRB 081228, $z=2.03^{+0.16}_{-0.14}$ for GRB 081230 and $z=1.28^{+0.16}_{-0.15}$ for GRB 090530. Combining the response from UVOT with ground-based observatories and in particular GROND operating in the optical/NIR wavelength regime, reliable photo-zs can be obtained from z ~ 1.0 out to z ~ 10, and possibly even at higher redshifts in some favorable cases, provided that these GRBs exist, are localized quickly, have sufficiently bright afterglows and are not heavily obscured.

preprint2010arXiv

Probing the molecular interstellar medium of M82 with Herschel-SPIRE spectroscopy

We present the observations of the starburst galaxy M82 taken with the Herschel SPIRE Fourier Transform Spectrometer. The spectrum (194-671 μm) shows a prominent CO rotational ladder from J = 4-3 to 13-12 emitted by the central region of M82. The fundamental properties of the gas are well constrained by the high J lines observed for the first time. Radiative transfer modeling of these high-S/N 12CO and 13CO lines strongly indicates a very warm molecular gas component at ~500 K and pressure of ~3x10^6 K cm^-3, in good agreement with the H_2 rotational lines measurements from Spitzer and ISO. We suggest that this warm gas is heated by dissipation of turbulence in the interstellar medium (ISM) rather than X-rays or UV flux from the straburst. This paper illustrates the promise of the SPIRE FTS for the study of the ISM of nearby galaxies.

preprint2010arXiv

Radial distribution of gas and dust in the two spiral galaxies M99 and M100

By combining Herschel-SPIRE data with archival Spitzer, HI, and CO maps, we investigate the spatial distribution of gas and dust in the two famous grand-design spirals M99 and M100 in the Virgo cluster. Thanks to the unique resolution and sensitivity of the Herschel-SPIRE photometer, we are for the first time able to measure the distribution and extent of cool, submillimetre (submm)-emitting dust inside and beyond the optical radius. We compare this with the radial variation in both the gas mass and the metallicity. Although we adopt a model-independent, phenomenological approach, our analysis provides important insights. We find the dust extending to at least the optical radius of the galaxy and showing breaks in its radial profiles at similar positions as the stellar distribution. The colour indices f350/f500 and f250/f350 decrease radially consistent with the temperature decreasing with radius. We also find evidence of an increasing gas to dust ratio with radius in the outer regions of both galaxies.

preprint2010arXiv

Selection of ULIRGs in Infrared and Submm Surveys

We examine the selection characteristics of infrared and sub-mm surveys with IRAS, Spitzer, BLAST, Herschel and SCUBA and identify the range of dust temperatures these surveys are sensitive to, for galaxies in the ULIRG luminosity range (12<log(LIR)<13), between z=0 and z=4. We find that the extent of the redshift range over which surveys are unbiased is a function of the wavelength of selection, flux density limit and ULIRG luminosity. Short wavelength (<200μm) surveys with IRAS, Spitzer/MIPS and Herschel/PACS are sensitive to all SED types in a large temperature interval (17-87K), over a substantial fraction of their accessible redshift range. On the other hand, long wavelength (>200μm) surveys with BLAST, Herschel/ SPIRE and SCUBA are significantly more sensitive to cold ULIRGs, disfavouring warmer SEDs even at low redshifts. We evaluate observations in the context of survey selection effects, finding that the lack of cold ULIRGs in the local (z<0.1) Universe is not a consequence of selection and that the range of ULIRG temperatures seen locally is only a subset of a much larger range which exists at high redshift. We demonstrate that the local luminosity-temperature (L-T) relation, which indicates that more luminous sources are also hotter, is not applicable in the distant Universe when extrapolated to the ULIRG regime, because the scatter in observed temperatures is too large. Finally, we show that the difference between the ULIRG temperature distributions locally and at high redshift is not the result of galaxies becoming colder due to an L-T relation which evolves as a function of redshift. Instead, they are consistent with a picture where the evolution of the infrared luminosity function is temperature dependent, i.e. cold galaxies evolve at a faster rate than their warm counterparts.

preprint2010arXiv

SPIRE imaging of M82: cool dust in the wind and tidal streams

M82 is a unique representative of a whole class of galaxies, starbursts with superwinds, in the Very Nearby Galaxy Survey with Herschel. In addition, its interaction with the M81 group has stripped a significant portion of its interstellar medium from its disk. SPIRE maps now afford better characterization of the far-infrared emission from cool dust outside the disk, and sketch a far more complete picture of its mass distribution and energetics than previously possible. They show emission coincident in projection with the starburst wind and in a large halo, much more extended than the PAH band emission seen with Spitzer. Some complex substructures coincide with the brightest PAH filaments, and others with tidal streams seen in atomic hydrogen. We subtract the far-infrared emission of the starburst and underlying disk from the maps, and derive spatially-resolved far-infrared colors for the wind and halo. We interpret the results in terms of dust mass, dust temperature, and global physical conditions. In particular, we examine variations in the dust physical properties as a function of distance from the center and the wind polar axis, and conclude that more than two thirds of the extraplanar dust has been removed by tidal interaction, and not entrained by the starburst wind.

preprint2010arXiv

The central region of spiral galaxies as seen by Herschel. M81, M99 and M100

With appropriate spatial resolution, images of spiral galaxies in thermal infrared (~10 micron and beyond) often reveal a bright central component, distinct from the stellar bulge, superimposed on a disk with prominent spiral arms. ISO and Spitzer studies have shown that much of the scatter in the mid-infrared colors of spiral galaxies is related to changes in the relative importance of these two components, rather than to other modifications, such as the morphological type or star formation rate, that affect the properties of the galaxy as a whole. With the Herschel imaging capability from 70 to 500 micron, we revisit this two-component approach at longer wavelengths, to see if it still provides a working description of the brightness distribution of galaxies, and to determine its implications on the interpretation of global far-infrared properties of galaxies.

preprint2010arXiv

The dust morphology of the elliptical Galaxy M86 with SPIRE

We present Herschel-SPIRE observations at 250-500um of the giant elliptical galaxy M86 and examine the distribution of the resolved cold dust emission and its relation with other galactic tracers. The SPIRE images reveal three dust components: emission from the central region; a dust lane extending north-south; and a bright emission feature 10kpc to the south-east. We estimate that approximately 10^6 solar masses of dust is spatially coincident with atomic and ionized hydrogen, originating from stripped material from the nearby spiral NGC4438 due to recent tidal interactions with M86. The gas-to-dust ratio of the cold gas component ranges from ~20-80. We discuss the different heating mechanisms for the dust features.

preprint2010arXiv

The environments of z~1 Active Galactic Nuclei at 3.6um

We present an analysis of a large sample of AGN environments at z~1 using stacked Spitzer data at 3.6um. The sample contains type-1 and type-2 AGN in the form of quasars and radio galaxies, and spans a large range in both optical and radio luminosity. We find, on average, that 2 to 3 massive galaxies containing a substantial evolved stellar population lie within a 200-300 kpc radius of the AGN, constituting a >8-sigma excess relative to the field. Secondly, we find evidence for the environmental source density to increase with the radio luminosity of AGN, but not with black-hole mass. This is shown first by dividing the AGN into their classical AGN types, where we see more significant over-densities in the fields of the radio-loud AGN. If instead we dispense with the classical AGN definitions, we find that the source over-density as a function of radio luminosity for all our AGN exhibits a positive correlation. One interpretation of this result is that the Mpc-scale environment is in some way influencing the radio emission that we observe from AGN. This could be explained by the confinement of radio jets in dense environments leading to enhanced radio emission or, alternatively, may be linked to more rapid black-hole spin brought on by galaxy mergers.

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 Space Observatory view of dust in M81

We use Herschel Space Observatory data to place observational constraints on the peak and Rayleigh-Jeans slope of dust emission observed at 70-500 microns in the nearby spiral galaxy M81. We find that the ratios of wave bands between 160 and 500 microns are primarily dependent on radius but that the ratio of 70 to 160 micron emission shows no clear dependence on surface brightness or radius. These results along with analyses of the spectral energy distributions imply that the 160-500 micron emission traces 15-30 K dust heated by evolved stars in the bulge and disc whereas the 70 micron emission includes dust heated by the active galactic nucleus and young stars in star forming regions.

preprint2010arXiv

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

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

preprint2010arXiv

The Herschel-SPIRE instrument and its in-flight performance

The Spectral and Photometric Imaging Receiver (SPIRE), is the Herschel Space Observatory`s submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 microns, and an imaging Fourier Transform Spectrometer (FTS) which covers simultaneously its whole operating range of 194-671 microns (447-1550 GHz). The SPIRE detectors are arrays of feedhorn-coupled bolometers cooled to 0.3 K. The photometer has a field of view of 4' x 8', observed simultaneously in the three spectral bands. Its main operating mode is scan-mapping, whereby the field of view is scanned across the sky to achieve full spatial sampling and to cover large areas if desired. The spectrometer has an approximately circular field of view with a diameter of 2.6'. The spectral resolution can be adjusted between 1.2 and 25 GHz by changing the stroke length of the FTS scan mirror. Its main operating mode involves a fixed telescope pointing with multiple scans of the FTS mirror to acquire spectral data. For extended source measurements, multiple position offsets are implemented by means of an internal beam steering mirror to achieve the desired spatial sampling and by rastering of the telescope pointing to map areas larger than the field of view. The SPIRE instrument consists of a cold focal plane unit located inside the Herschel cryostat and warm electronics units, located on the spacecraft Service Module, for instrument control and data handling. Science data are transmitted to Earth with no on-board data compression, and processed by automatic pipelines to produce calibrated science products. The in-flight performance of the instrument matches or exceeds predictions based on pre-launch testing and modelling: the photometer sensitivity is comparable to or slightly better than estimated pre-launch, and the spectrometer sensitivity is also better by a factor of 1.5-2.

preprint2010arXiv

The warm absorber and X-ray variability of the Seyfert 1 galaxy NGC 3516 as seen by the XMM-Newton RGS

We present a new analysis of the soft and medium energy X-ray spectrum of the Seyfert 1 galaxy NGC 3516 taken with the Reflection Grating Spectrometer (RGS) and European Photon Imaging Camera (EPIC) on board the XMM-Newton observatory. We examine four observations made in October 2006. We investigate whether the observed variability is due to absorption by the warm absorber and/or is intrinsic to the source emission. We analyse in detail the EPIC-pn and RGS spectra of each observation separately. The warm absorber in NGC 3516 is found to consist of three phases of ionisation, two of which have outflow velocities of more than 1000 km/s. The third phase (the least ionised one) is much slower at 100 km/s. One of the high ionisation phases, with log xi of 2.4, is found to have a partial covering fraction of about 60%. It has previously been suggested that the passage of a cloud, part of a disk wind, in front of the source (producing a change in the covering fraction) was the cause of a significant dip in the lightcurve during one of the observations. From our modelling of the EPIC-pn and RGS spectra, we find that variation in the covering fraction cannot be solely responsible for this. We show that intrinsic change in the source continuum plays a much more significant role in explaining the observed flux and spectral variability than originally thought.

preprint2009arXiv

Dust and Metal Column Densities in Gamma-Ray Burst Host Galaxies

In this paper we present the results from the analysis of a sample of 28 gamma-ray burst (GRB) afterglow spectral energy distributions, spanning the X-ray through to near-infrared wavelengths. This is the largest sample of GRB afterglow spectral energy distributions thus far studied, providing a strong handle on the optical depth distribution of soft X-ray absorption and dust-extinction systems in GRB host galaxies. We detect an absorption system within the GRB host galaxy in 79% of the sample, and an extinction system in 71% of the sample, and find the Small Magellanic Cloud (SMC) extinction law to provide an acceptable fit to the host galaxy extinction profile for the majority of cases, consistent with previous findings. The range in the soft X-ray absorption to dust-extinction ratio, N_{H,X}/Av, in GRB host galaxies spans almost two orders of magnitude, and the typical ratios are significantly larger than those of the Magellanic Clouds or Milky Way. Although dust destruction could be a cause, at least in part, for the large N_{H,X}/Av ratios, the good fit provided by the SMC extinction law for the majority of our sample suggests that there is an abundance of small dust grains in the GRB environment, which we would expect to have been destroyed if dust destruction were responsible for the large N_{H,X}/Av ratios. Instead, our analysis suggests that the distribution of N_{H,X}/Av in GRB host galaxies may be mostly intrinsic to these galaxies, and this is further substantiated by evidence for a strong negative correlation between N_{H,X}/Av and metallicity for a subsample of GRB hosts with known metallicity. Furthermore, we find the N_{H,X}/Av ratio and metallicity for this subsample of GRBs to be comparable to the relation found in other more metal-rich galaxies.

preprint2009arXiv

Swift and Fermi observations of the early afterglow of the short Gamma-Ray Burst 090510

We present the observations of GRB090510 performed by the Fermi Gamma-Ray Space Telescope and the Swift observatory. This is a bright, short burst that shows an extended emission detected in the GeV range. Furthermore, its optical emission initially rises, a feature so far observed only in long bursts, while the X-ray flux shows an initial shallow decrease, followed by a steeper decay. This exceptional behavior enables us to investigate the physical properties of the GRB outflow, poorly known in short bursts. We discuss internal shock and external shock models for the broadband energy emission of this object.

preprint2009arXiv

The central energy source of 70micron-selected galaxies: Starburst or AGN?

We present the first AGN census in a sample of 61 galaxies selected at 70microns, a wavelength which should strongly favour the detection of star-forming systems. For the purpose of this study we take advantage of deep Chandra X-ray and Spitzer infrared (3.6-160micron) data, as well as optical spectroscopy and photometry from the Deep Extragalactic Evolutionary Probe 2 (DEEP2) survey for the Extended Groth Strip (EGS) field. We investigate spectral line diagnostics ([OIII]/Hbeta and [NeIII]/[OII] ratios, Hdelta Balmer absorption line equivalent widths and the strength of the 4000Ang break), X-ray luminosities and spectral energy distributions (SEDs). We find that the 70micron sources are undergoing starburst episodes and are therefore characterised by a predominance of young stars. In addition, 13 per cent of the sources show AGN signatures and hence potentially host an AGN. When the sample is split into starbursts (SBs, 10^10<L_IR<10^11 L_solar), Luminous InfraRed Galaxies (LIRGs, 10^11<L_IR<10^12 L_solar) and UltraLuminous InfraRed Galaxies (ULIRGs,10^12<L_IR<10^13 L_solar), the AGN fraction becomes 0, 11 and 23 per cent respectively, showing an increase with total infrared luminosity. However, by examining the sources' panchromatic SEDs, we conclude that although the AGN is energetically important in 1 out of 61 objects, all 70micron-selected galaxies are primarily powered by star-formation.

preprint2009arXiv

The Comoving Infrared Luminosity Density: Domination of Cold Galaxies across 0<z<1

In this paper we examine the contribution of galaxies with different infrared (IR) spectral energy distributions (SEDs) to the comoving infrared luminosity density, a proxy for the comoving star formation rate (SFR) density. We characterise galaxies as having either a cold or hot IR SED depending upon whether the rest-frame wavelength of their peak IR energy output is above or below 90um. Our work is based on a far-IR selected sample both in the local Universe and at high redshift, the former consisting of IRAS 60um-selected galaxies at z<0.07 and the latter of Spitzer 70um selected galaxies across 0.1<z<1. We find that the total IR luminosity densities for each redshift/luminosity bin agree well with results derived from other deep mid/far-IR surveys. At z<0.07 we observe the previously known results: that moderate luminosity galaxies (L_IR<10^11 Lsun) dominate the total luminosity density and that the fraction of cold galaxies decreases with increasing luminosity, becoming negligible at the highest luminosities. Conversely, above z=0.1 we find that luminous IR galaxies (L_IR>10^11 Lsun), the majority of which are cold, dominate the IR luminosity density. We therefore infer that cold galaxies dominate the IR luminosity density across the whole 0<z<1 range, hence appear to be the main driver behind the increase in SFR density up to z~1 whereas local luminous galaxies are not, on the whole, representative of the high redshift population.

preprint2009arXiv

The Link Between SCUBA and Spitzer: Cold Galaxies at z<1

We show that the far-IR properties of distant Luminous and Ultraluminous InfraRed Galaxies (LIRGs and ULIRGs) are on average divergent from analogous sources in the local Universe. Our analysis is based on Spitzer MIPS and IRAC data of L_IR>10^10 L_solar, 70um-selected objects in the 0.1<z<2 redshift range and supported by a comparison with the IRAS Bright Galaxy Sample. The majority of the objects in our sample are described by Spectral Energy Distributions (SEDs) which peak at longer wavelengths than local sources of equivalent total infrared luminosity. This shift in SED peak wavelength implies a noticeable change in the dust and/or star-forming properties from z~0 to the early Universe, tending towards lower dust temperatures, indicative of strong evolution in the cold dust, `cirrus', component. We show that these objects are potentially the missing link between the well-studied local IR-luminous galaxies, Spitzer IR populations and SCUBA sources -- the z<1 counterparts of the cold z>1 SubMillimetre Galaxies (SMGs) discovered in blank-field submillimetre surveys. The Herschel Space Observatory is well placed to fully characterise the nature of these objects, as its coverage extends over a major part of the far-IR/submm SED for a wide redshift range.

preprint2009arXiv

The XMM-Newton Wide Angle Survey (XWAS): the X-ray spectrum of type-1 AGN

We discuss the broad band X-ray properties of one of the largest samples of X-ray selected type-1 AGN to date (487 objects in total), drawn from the XMM-Newton Wide Angle Survey. The objects cover 2-10 keV luminosities from ~10^{42}-10^{45} erg s^{-1} and are detected up to redshift ~4. We constrain the overall properties of the broad band continuum, soft excess and X-ray absorption, along with their dependence on the X-ray luminosity and redshift and we discuss the implications for models of AGN emission. We constrained the mean spectral index of the broad band X-ray continuum to <Gamma>=1.96+-0.02 with intrinsic dispersion sigma=0.27_{-0.02}^{+0.01}. The continuum becomes harder at faint fluxes and at higher redshifts and luminosities. The dependence of Gamma with flux is likely due to undetected absorption rather than to spectral variation. We found a strong dependence of the detection efficiency of objects on the spectral shape which can have a strong impact on the measured mean continuum shapes of sources at different redshifts and luminosities. We detected excess absorption in ~3% of our objects, with column densities ~a few x10^{22} cm^{-2}. The apparent mismatch between the optical classification and X-ray properties of these objects is a challenge for the standard AGN unification model. We found that the fraction of objects with detected soft excess is ~36%. Using a thermal model, we constrained the soft excess mean temperature and intrinsic dispersion to <kT>~100 eV and sigma~34 eV. The origin of the soft excess as thermal emission from the accretion disk or Compton scattered disk emission is ruled out on the basis of the temperatures detected and the lack of correlation of the measured temperature with the X-ray luminosity (abridged).

preprint2008arXiv

An extremely luminous X-ray outburst at the birth of a supernova

Massive stars end their short lives in spectacular explosions, supernovae, that synthesize new elements and drive galaxy evolution. Throughout history supernovae were discovered chiefly through their delayed optical light, preventing observations in the first moments (hours to days) following the explosion. As a result, the progenitors of some supernovae and the events leading up to their violent demise remain intensely debated. Here we report the serendipitous discovery of a supernova at the time of explosion, marked by an extremely luminous X-ray outburst. We attribute the outburst to the break-out of the supernova shock-wave from the progenitor, and show that the inferred rate of such events agrees with that of all core-collapse supernovae. We forecast that future wide-field X-ray surveys will catch hundreds of supernovae each year in the act of explosion, and thereby enable crucial neutrino and gravitational wave detections that may ultimately unravel the explosion mechanism.

preprint2007arXiv

Dust and Gas in the Local Environments of Gamma-Ray Bursts

Using a sample of gamma-ray burst (GRB) afterglows detected by both the X-Ray and the UV/Optical Telescopes (XRT and UVOT) on Swift, we modelled the spectral energy distributions (SEDs) to determine gas column densities and dust extinction in the GRB local environment. In six out of seven cases we find an X-ray absorber associated with the GRB host galaxy with column density (assuming solar abundances) ranging from (0.8 - 7.7)x10^{21}cm^{-2}. We determine the rest-frame visual extinction A_V using the SMC, LMC and Galactic extinction curves to model the dust in the GRB host galaxy, and this ranges from A_V = 0.12\pm 0.04 to A_V = 0.65^{+0.08}_{-0.07}. The afterglow SEDs were typically best fit by a model with an SMC extinction curve. In only one case was the GRB afterglow better modelled by a Galactic extinction curve, which has a prominent absorption feature at 2175angstrom. We investigate the selection effects present in our sample and how these might distort the true distribution of A_V in GRB host galaxies. We estimate that GRBs with no afterglow detected blueward of 5500angstrom have average rest-frame visual extinctions almost eight times those observed in the optically bright population of GRBs. This may help account for the ~1/3 of GRBs observed by Swift that have no afterglow detected by UVOT.

preprint2007arXiv

Early afterglow detection in the Swift observations of GRB 050801

We present results of Swift optical, UV and X-ray observations of the afterglow of GRB 050801. The source is visible over the full optical, UV and X-ray energy range of the Swift UVOT and XRT instruments.Both optical and X-ray lightcurves exhibit a broad plateau (Δt/t ~ 1) during the first few hundred seconds after the gamma-ray event. We investigate the multiwavelength spectral and timing properties of the afterglow, and we suggest that the behaviour at early times is compatible with an energy injection by a newly born magnetar with a period of a few tenths of a millisecond, which keeps the forward shock refreshed over this short interval by irradiation. Reverse shock emission is not observed. Its suppression might be due to GRB ejecta being permeated by high magnetic fields, as expected for outflows powered by a magnetar.Finally, the multiwavelength study allows a determination of the burst redshift, z=1.56.

preprint2007arXiv

Extreme Properties Of GRB061007: A Highly Energetic Or A Highly Collimated Burst?

GRB061007 is the brightest gamma-ray burst (GRB) to be detected by Swift and is accompanied by an exceptionally luminous afterglow that had a V-band magnitude <11.1 at 80s after the prompt emission. From the start of the Swift observations the afterglow decayed as a power law with a slope of α_X=1.66+/-0.01 in the X-ray and α_{opt}=1.64+/-0.01 in the UV/optical, up to the point that it was no longer detected above background in the optical or X-ray bands. The brightness of this GRB and the similarity in the decay rate of the X-ray, optical and gamma-ray emission from 100s after the trigger distinguish this burst from others and present a challenge to the fireball model. The lack of a cooling or jet break in the afterglow up to \~10^5s constrains any model that can produce the large luminosity observed in GRB061007, which we found to require either an excessively large kinetic energy or highly collimated outflow. Analysis of the multi-wavelength spectral and high-resolution temporal data taken with Swift suggest an early time jet-break to be a more plausible scenario. This must have occurred within 80s of the prompt emission, which places an upper limit on the jet opening angle of θ_j=0.8deg. Such a highly collimated outflow resolves the energy budget problem presented in a spherical emission model, reducing the isotropic equivalent energy of this burst to E_γ^{corr}=10^{50} ergs; consistent with other GRBs.

preprint2007arXiv

GRB 061121: Broadband spectral evolution through the prompt and afterglow phases of a bright burst

Swift triggered on a precursor to the main burst of GRB 061121 (z=1.314), allowing observations to be made from the optical to gamma-ray bands. Many other telescopes, including Konus-Wind, XMM-Newton, ROTSE and the Faulkes Telescope North, also observed the burst. The gamma-ray, X-ray and UV/optical emission all showed a peak ~75s after the trigger, although the optical and X-ray afterglow components also appear early on - before, or during, the main peak. Spectral evolution was seen throughout the burst, with the prompt emission showing a clear positive correlation between brightness and hardness. The Spectral Energy Distribution (SED) of the prompt emission, stretching from 1eV up to 1MeV, is very flat, with a peak in the flux density at ~1keV. The optical-to-X-ray spectra at this time are better fitted by a broken, rather than single, power-law, similar to previous results for X-ray flares. The SED shows spectral hardening as the afterglow evolves with time. This behaviour might be a symptom of self-Comptonisation, although circumstellar densities similar to those found in the cores of molecular clouds would be required. The afterglow also decays too slowly to be accounted for by the standard models. Although the precursor and main emission show different spectral lags, both are consistent with the lag-luminosity correlation for long bursts. GRB 061121 is the instantaneously brightest long burst yet detected by Swift. Using a combination of Swift and Konus-Wind data, we estimate an isotropic energy of 2.8x10^53 erg over 1keV - 10MeV in the GRB rest frame. A probable jet break is detected at ~2x10^5s, leading to an estimate of ~10^51 erg for the beaming-corrected gamma-ray energy.

preprint2007arXiv

Photometric Calibration of the Swift Ultraviolet/Optical Telescope

We present the photometric calibration of the Swift UltraViolet/Optical Telescope (UVOT) which includes: optimum photometric and background apertures, effective area curves, colour transformations, conversion factors for count rates to flux, and the photometric zero points (which are accurate to better than 4 per cent) for each of the seven UVOT broadband filters. The calibration was performed with observations of standard stars and standard star fields that represent a wide range of spectral star types. The calibration results include the position dependent uniformity, and instrument response over the 1600-8000A operational range. Because the UVOT is a photon counting instrument, we also discuss the effect of coincidence loss on the calibration results. We provide practical guidelines for using the calibration in UVOT data analysis. The results presented here supersede previous calibration results.

preprint2007arXiv

Swift observations of GRB 070110: an extraordinary X-ray afterglow powered by the central engine

We present a detailed analysis of Swift multi-wavelength observations of GRB 070110 and its remarkable afterglow. The early X-ray light curve, interpreted as the tail of the prompt emission, displays a spectral evolution already seen in other gamma-ray bursts. The optical afterglow shows a shallow decay up to ~2 d after the burst, which is not consistent with standard afterglow models. The most intriguing feature is a very steep decay in the X-ray flux at ~20 ks after the burst, ending an apparent plateau. The abrupt drop of the X-ray light curve rules out an external shock as the origin of the plateau in this burst and implies long-lasting activity of the central engine. The temporal and spectral properties of the plateau phase point towards a continuous central engine emission rather than the episodic emission of X-ray flares. We suggest that the observed X-ray plateau is powered by a spinning down central engine, possibly a millisecond pulsar, which dissipates energy at an internal radius before depositing energy into the external shock.

preprint2007arXiv

The XMM large scale structure survey: optical vs. X-ray classifications of active galactic nuclei and the unified scheme

Our goal is to characterize AGN populations by comparing their X-ray and optical classifications. We present a sample of 99 spectroscopically identified X-ray point sources in the XMM-LSS survey which are significantly detected in the [2-10] keV band, and with more than 80 counts. We performed an X-ray spectral analysis for all of these 99 X-ray sources. Introducing the fourfold point correlation coefficient, we find only a mild correlation between the X-ray and the optical classifications, as up to 30% of the sources have differing X-ray and optical classifications: on one hand, 10% of the type 1 sources present broad emission lines in their optical spectra and strong absorption in the X-rays. These objects are highly luminous AGN lying at high redshift and thus dilution effects are totally ruled out, their discrepant nature being an intrinsic property. Their X-ray luminosities and redshifts distributions are consistent with those of the unabsorbed X-ray sources with broad emission lines. On the other hand, 25/32 are moderate luminosity AGN, which are both unabsorbed in the X-rays and only present narrow emission lines in their optical spectra. The majority of them have an optical spectrum which is representative of the host galaxy. We finally infer that dilution of the AGN by the host galaxy seems to account for their nature. 5/25 have been defined as Seyfert 2. In conclusion, most of these 32 discrepant cases can be accounted for by the standard AGN unified scheme, as its predictions are not met for only 12% of the 99 X-ray sources. ABRIDGED

preprint2006arXiv

Swift observations of GRB050712

We present the results of X-ray and optical observations of GRB050712 performed by Swift. The X-ray lightcurve of this burst exhibits episodes of flares in the first 1000s, the same epoch at which the UVOT detected an optical counterpart. A shallow X-ray decay, with a decay slope of -0.73, followed and lasted ~70ks. This behaviour can be explained in terms of activity of GRB 'inner engine', with the possibility of that the last flare is caused by the interaction of the ejecta with the surrounding medium. We also find interesting spectral parameters for the X-ray emission. In particular, data suggest the presence of an intrinsic absorption in the first 1000s, which can be explained if circumbust medium clouds lie along the line of sight.

preprint2006arXiv

The SCUBA HAlf Degree Extragalactic Survey (SHADES) -- II. Submillimetre maps, catalogue and number counts

We present the maps, source catalogue and number counts of the largest, most complete and unbiased extragalactic submillimetre survey ever undertaken: the 850-micron SCUBA HAlf Degree Extragalactic Survey (SHADES). Using the Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT), SHADES mapped two separate regions of sky: the Subaru/XMM-Newton Deep Field (SXDF) and the Lockman Hole East (LH). These SCUBA maps cover 720 square arcmin with an RMS noise level of about 2 mJy and have uncovered >100 submillimetre galaxies. In order to ensure the utmost robustness of the resulting source catalogue, data reduction was independently carried out by four sub-groups drawn within the SHADES team, providing an unprecedented degree of reliability with respect to other SCUBA catalogues available from the literature. Individual source lists from the four groups were combined to produce a robust 120-object SHADES catalogue; an invaluable resource for follow-up campaigns aiming to study the properties of a complete and consistent sample of submillimetre galaxies. We present differential and integral source number counts of submillimetre galaxies and find that the differential counts are better fit with a broken power-law or a Schechter function than with a single power-law; the SHADES data alone significantly show that a break is required at several mJy, although the precise position of the break is not well constrained. We also find that an 850-micron survey complete down to 2 mJy would resolve 20-30 per cent of the Far-IR background into point sources. [abridged]

preprint2004arXiv

Mutiwavelength Observations of Radio Galaxy 3C 120 with XMM-Newton

We present XMM-Newton observations of the radio galaxy 3C 120. The hard X-ray spectrum contains a marginally resolved Fe I K-alpha emission line with FWHM=9,000 km/s and an equivalent width of 57 eV. The line arises via fluorescence in a broad-line region with covering fraction of 0.4. There is no evidence of relativistically broad Fe K-alpha, contrary to some previous reports. The normal equivalent widths of the X-ray and optical emission lines exclude a strongly beamed synchrotron component to the hard X-ray and optical continua. There is an excess of 0.3-2 keV soft X-ray continuum over an extrapolation of the hard X-ray power-law, which may arise in a disk corona. Analysis of an archival Chandra image shows that extended emission from the jet and other sources contributes <3% of the total X-ray flux. A break in the X-ray spectrum below 0.6 keV indicates an excess neutral hydrogen column density of N_H=1.57 * 10^21 cm^{-2}. However, the neutral absorber must have an oxygen abundance of <1/50 of the solar value to explain the absence of an intrinsic or intervening O I edge. There is no ionized absorption in the soft X-ray spectrum, but there is a weak, narrow O VIII Ly-alpha emission line. We do not detect previously claimed O VIII absorption from the intervening intergalactic medium. Radio observations at 37 GHz show a fast, high frequency flare starting 8 days after the XMM-Newton observation. However, this has no obvious effect on the X-ray spectrum. The X-ray spectrum, including the soft excess, became harder as the X-ray flux decreased, with an estimated pivot energy of 40 keV. The UV and soft X-ray fluxes are strongly correlated over the 120 ks duration of the XMM-Newton observation. This is qualitatively consistent with Comptonization of UV photons by a hot corona. (Abridged)

preprint2004arXiv

Relativistic O VIII Emission and Ionized Outflow in NGC 4051

We present XMM-Newton RGS observations of the soft X-ray spectrum of NGC 4051, and explore their implications for the inner accretion disk and ionized outflow in the active galactic nucleus. We fit the soft X-ray excess with a relativistically broadened O VIII recombination spectrum, including the entire line series and recombination continuum. This plus an underlying power law continuum provides a much better fit to the soft X-ray spectrum than a single temperature or disk blackbody plus power law. The emission line profiles, computed for a Kerr metric around a maximally rotating black hole, reveal a sharply peaked disk emissivity law and inner radius smaller than 1.7 R_G. The spectrum also includes narrow absorption and emission lines from C, N, O, Ne, and Fe in an ionized outflow. Outflow column densities are relatively low and do not create significant edges in the spectrum. The small amount of absorption bolsters confidence in the detection of relativistic emission line features. The narrow-line emitter has a large (76%) global covering fraction, leading to strong forbidden lines and filling in of the resonance absorption lines. We also find broad C VI Ly-alpha and very broad O VII emission from the broad-line region. The narrow and broad-line regions span large ranges in ionization parameter, and may arise in a disk outflow. The ionized absorber has a large ionization range which is inconsistent with pressure equilibrium in a multiphase medium. The mass outflow rate exceeds the accretion rate by a factor of one thousand.

preprint2003arXiv

Multiwavelength studies of the Seyfert 1 galaxy NGC7469. II - X-ray and UV observations with XMM-Newton

We present an XMM-Newton observation of NGC 7469, including studies of the X-ray and UV variability, 0.2-10 keV spectral continuum, Fe K-alpha emission line and the first-ever high-resolution X-ray spectrum of the soft X-ray warm absorber. We compare the properties of this X-ray warm absorber with the UV warm absorber as seen in a FUSE observation one year previously. The 0.2-10 keV spectral continuum is best fitted by a power-law plus two blackbody model. An Fe K-alpha emission line is visible which consists of a single narrow component and is well-modelled by a simple gaussian. Narrow absorption and emission lines in the soft X-ray RGS spectrum demonstrate the existence of a multi-phase warm absorber with a range in log xi of ~ 2 to ~ -2 where xi is in erg cm s^-1. The warm absorber is blueshifted by several hundred km s^-1. The highest-ionisation phase of the absorber is the best constrained and has an overall equivalent Hydrogen column of order 10^20 cm^-2; we find that its ionisation parameter is consistent with that of the warm emitter which generates the narrow emission lines. We identify this high ionisation absorber with the low-velocity phase of the UV absorber observed by FUSE.

preprint2002arXiv

Radio observations of the XMM-Newton/Chandra 13hr deep survey field

Our VLA observations of the XMM-Newton/Chandra 13hr deep survey field (see Page et al., this proceedings) result in one of the two deepest 1.4GHz radio maps ever made. Within the 15' radius field covered by the deep X-ray data (0.19 sq deg), a total of 556 radio sources are detected, down to a 4 sigma flux limit of 28uJy. Of the 214 Chandra sources, 55 have radio counterparts. The sub-arcsecond accuracy of the VLA and Chandra positions enable us to determine with high confidence the sources common to both surveys. Here we present the relationship between the X-ray and radio source populations at the faintest radio flux limits yet probed by such a study. We discuss how the X-ray/radio relationship differs as a function of optical morphology, ie between unresolved `stellar' objects and well resolved galaxies. We then discuss the origin of the X-ray and radio emission, ie AGN, starburst or a mixture of both, in these two classes of object.

preprint2002arXiv

Submillimeter Evidence for the Coeval Growth of Massive Black Holes and Galaxy Bulges

The correlation, found in nearby galaxies, between black hole mass and stellar bulge mass implies that the formation of these two components must be related. Here we report submillimeter photometry of eight x--ray absorbed active galactic nuclei which have luminosities and redshifts characteristic of the sources that produce the bulk of the accretion luminosity in the universe. The four sources with the highest redshifts are detected at 850 microns, with flux densities between 5.9 and 10.1 milliJanskies, and hence are ultraluminous infrared galaxies. Interpreting the submillimeter flux as emission from dust heated by starbursts, these results suggest that the majority of stars in spheroids were formed at the same time as their central black holes built up most of their mass by accretion, accounting for the observed demography of massive black holes in the local universe. The skewed rate of submillimeter detection with redshift is consistent with a high redshift epoch of star formation in radio quiet active galactic nuclei, similar to that seen in radio galaxies.

preprint2002arXiv

The XMM-Newton Survey Science Centre Medium Sensitivity Survey

We present preliminary results on the XMM-Newton Survey Science Centre medium sensitivity survey (XMS), with 0.5-4.5 keV flux limit 2 10^{-14} erg/cm2/s. At present, 19 fields have been examined with a total of 239 X-ray sources. Identifications for about 2/3 of these reveal that type 1 AGN dominate, but type 2 AGN, absorption-line galaxies and clusters as well as stars are also present. We also discuss on a few X-ray selected Broad-Absorption-Line QSOs found in this survey.

preprint2002arXiv

X-ray and optical properties of X-ray sources in the 13hr XMM-Newton/Chandra deep survey

The 13hr XMM-Newton/Chandra deep survey is the first of two extremely deep XMM-Newton fields observed by the XMM-OM consortium. A 120 ks Chandra mosaic, covering 0.2 square degrees, provides sensitive, confusion-free point source detection with sub-arcsecond positions, while the 200 ks XMM-Newton observation provides high quality X-ray spectroscopy over the same sky area. We have optical spectroscopic identifications for 70 X-ray sources. Of these, 42 are broad emission-line AGN with a wide range of redshifts. The optical counterparts of a further 23 sources are narrow emission line galaxies and absorption line galaxies. These 23 sources all lie at z < 1 and typically have lower X-ray luminosities than the broad-line AGN. About half of them show significant X-ray absorption and are almost certainly intrinsically absorbed AGN. However some of them have unabsorbed, AGN-like, power-law components in their X-ray spectra, but do not show broad emission lines in their optical spectra. These sources may be weak, unobscured AGN in bright galaxies and their existence at low redshifts could be a consequence of the strong cosmological evolution of AGN characteristic luminosities.

preprint2000arXiv

Soft X-ray emission lines from a relativistic accretion disk in MCG -6-30-15 and Mrk 766

XMM-Newton Reflection Grating Spectrometer (RGS) spectra of the Narrow Line Seyfert 1 galaxies MCG -6-30-15 and Mrk 766 are physically and spectroscopically inconsistent with standard models comprising a power-law continuum absorbed by either cold or ionized matter. We propose that the remarkably similar features detected in both objects in the 5 - 35 A band are H-like oxygen, nitrogen, and carbon emission lines, gravitationally redshifted and broadened by relativistic effects in the vicinity of a Kerr black hole. We discuss the implications of our interpretation, and demonstrate that the derived parameters can be physically self-consistent.

preprint1998arXiv

Clustering of X-ray selected Active Galactic Nuclei

A total of 235 Active Galactic Nuclei (AGN) from two different soft X-ray surveys (the ROSAT Deep Survey -DRS- and the ROSAT International X-ray Optical Survey -RIXOS-) with redshifts between 0 and 3.5 are used to study the clustering of X-ray selected AGN and its evolution. A 2 sigma significant detection of clustering of such objects is found on scales <40-80/h Mpc in the RIXOS sample, while no clustering is detected on any scales in the DRS sample. Assuming a single power law model for the spatial correlation function (SCF), quantitative limits on the AGN clustering have been obtained: a comoving correlation length 1.5<~ r_0 <~ 3.3/h Mpc is implied for comoving evolution, while 1.9 <~ r_0 <~ 4.8 for stable clustering and 2.2 <~ r_0 <~ 5.5 for linear evolution. These values are consistent with the correlation lengths and evolutions obtained for galaxy samples, but imply smaller amplitude or faster evolution than recent UV and optically selected AGN samples. We also constrain the ratio of bias parameters between X-ray selected AGN and IRAS galaxies to be <~1.7 on scales <~ 10/h Mpc, a significantly smaller value than is inferred from local large-scale dynamical studies.

preprint1997arXiv

Optical and X-ray properties of the RIXOS AGN: II - Emission lines

We present the optical and UV emission line properties of 160 X-ray selected AGN taken from the RIXOS survey (including Halpha, Hbeta, [OIII]5007, MgII2798 and CIII]1909). This sample is believed to contain a mixture of absorbed and unabsorbed objects, with column densities up to 4e21 cm-2. Although the distribution of the [OIII] EW for the RIXOS AGN is typical of optically selected samples, the Balmer line EWs are relatively low. This is consistent with the presence of a dust absorber between the broad and narrow line regions (eg. a molecular torus), and intrinsically weak optical line emission. We find Baldwin effects in CIII] and MgII, and a positive response of the MgII line to its ionizing continuum. There is a strong correlation between the EW and FWHM of MgII, which may be similar to that seen in other samples for Hbeta. We demonstrate that this is consistent with models which suggest two line-emitting zones, a `very broad line region' (VBLR) and an `intermediate line region' (ILR). The correlation between EW and FWHM in MgII may be a physical characteristic of the ILR or it may reflect a geometric dependence. We found no correlation between the Hbeta FWHM and the slope of the X-ray spectrum, however this may be due to the effects of dust absorption which suppresses the broad Hbeta component, masking any relationship. The Halpha FWHM does tend to be narrow when alpha_X is soft, and broadens as alpha_X hardens, although the formal probability for this correlation is low (91 per cent). If the distribution of alpha_X in the RIXOS sample reflects the level of intrinsic absorption in these AGN, the data suggest a possible link between the velocity of the Balmer line-emitting region and the amount of absorbing material beyond.

preprint1996arXiv

Optical and X-ray properties of the RIXOS AGN: I - The continua

We present measurements of the optical and X-ray continua of 108 AGN (Seyfert 1s and quasars) from the Rosat International X-ray/Optical Survey (RIXOS). The sample covers a wide range in redshift (0<z<3.3), in X-ray spectral slope (-1.5 <ax<2.6) and in optical-to-X-ray ratio (0.4$<aox<1.5). A correlation is found between ax and aox; similar correlations have recently been reported in other X-ray and optical samples. We also identify previously unreported relationships between the optical slope (aopt) and ax (particularly at high redshifts) and between aopt and aox. These trends show that the overall optical-to-X-ray conti nuum changes from convex to concave as ax hardens, demonstrating a strong beha vioural link between the optical/UV big blue bump (BBB) and the soft X-ray ex cess, which is consistent with them being part of the same spectral component. By constructing models of the optical-to-X-ray continuum, we demonstrate that the observed correlations are consistent with an intrinsic spectrum which is ab sorbed through different amounts of cold gas and dust. The intrinsic spectrum is the sum of an optical-to-soft X-ray `big bump' component and an ax=1 power law; the column density of the cold gas ranges from 0 to ~4E21 cm-2 and the dust-to-gas ratio is assumed to be Galactic. The `big bump' may be represented by a T~1E6 K thermal bremsstrahlung or an accretion disk with a surrounding hot corona. The scatter in the data can accommodate a wide range in big bump tempe rature (or black hole mass) and strength. A source for the absorbing gas may be the dusty, molecular torus which lies beyond the broad line-emitting regions, although with a much lower column density than observed in Seyfert 2 galaxies. Alternatively, it may be the bulge of a spiral host galaxy or an elliptical host galaxy.

preprint1996arXiv

The Luminosity Function Evolution of Soft X--ray selected AGN in the RIXOS survey

A sample of 198 soft X--ray selected active galactic nuclei (AGN) from the ROSAT International X--ray Optical Survey (RIXOS), is used to investigate the X--ray luminosity function and its evolution. RIXOS, with a flux limit of 3E-14 erg s-1 cm-2 (0.5 to 2.0 keV), samples a broad range in redshift over 20 deg^2 of sky, and is almost completely identified; it is used in combination with the Einstein Extended Medium Sensitivity Survey (EMSS), to give a total sample of over 600 AGN. We find the evolution of AGN with redshift to be consistent with pure luminosity evolution (PLE) models in which the rate of evolution slows markedly or stops at high redshifts z>1.8. We find that this result is not affected by the inclusion, or exclusion, of narrow emission line galaxies at low redshift in the RIXOS and EMSS samples, and is insensitive to uncertainties in the conversion between flux values measured with ROSAT and Einstein. We confirm, using a model independent Ve/Va test, that our survey is consistent with no evolution at high redshifts.