Source author record

D. Burgarella

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

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

86 published item(s)

preprint2022arXiv

Comparison of the star formation of X-ray selected active galactic nuclei in eFEDS with star-forming galaxies

We use $\sim 1800$ X-ray Active Galactic Nuclei (AGN) in the eROSITA Final Equatorial-Depth Survey (eFEDS), that span over two orders of magnitude in X-ray luminosity, $\rm L_{X,2-10keV} \approx 10^{43-45}\,ergs^{-1}$, and compare their star-formation rate (SFR) relative to that of non-AGN star-forming systems, at $\rm 0.5<z<1.5$. For that purpose, we compile a reference galaxy catalogue with $\sim 17000$ sources. Both samples have the same photometric coverage, from optical to far-infrared. We construct the spectral energy distributions (SEDs) of all sources and fit them using the CIGALE code, utilizing the same templates and parametric grid for both samples. We account for the mass incompleteness and exclude quiescent systems from both datasets. These allow us to compare the SFR of the two populations in a uniform manner, minimising systematic effects. Based on our analysis, AGN at low and moderate L$_X$ ($\rm L_{X,2-10keV} < 10^{44}\,ergs^{-1}$), have SFR that is lower, or at most, equal to that of star-forming galaxies, in agreement with previous studies. The large number of luminous X-ray AGN available in the eFEDS dataset, enable us to expand our investigations at higher L$_X$ to test previous, tentative results. At $\rm L_{X,2-10keV} > 10^{44.2}\,ergs^{-1}$, the SFR of AGN appears enhanced, by $\sim 30\%$, compared to that of star-forming sources, for systems with stellar mass, $\rm 10.5 < log\,[M_*(M_\odot)] < 11.5$, confirming indications found in previous studies. The most massive sources $\rm log\,[M_*(M_\odot)] > 11.5$, present a flat SFR$_{norm}$-L$_X$ relation up to $\rm L_{X,2-10keV} \sim 10^{44.5}\,ergs^{-1}$, with SFR similar to that of star-forming galaxies. However, at higher L$_X$ ($\rm L_{X,2-10keV} \sim 10^{45}\,ergs^{-1}$), we find indications that the SFR of these massive AGN hosts may be enhanced compared to that of non-AGN systems.

preprint2022arXiv

IR characteristic emission and dust properties of star-forming galaxies at 4.5 $<$ z $<$ 6.2

The luminosity functions at z < 4 - 5 suggest that most galaxies have a relatively low stellar mass (logM_star = 10) and a low dust attenuation (A_FUV = 1.0). The physical properties of these objects are quite homogeneous. We used an approach where we combined their rest-frame far-infrared and submillimeter emissions and utilized the universe and the redshift as a spectrograph to increase the amount of information in a collective way. From a subsample of 27 ALMA-detected galaxies at z > 4.5, we built an infrared spectral energy distribution composite template. It was used to fit, with CIGALE, the 105 galaxies (detections and upper limits) in the sample from the FUV to the FIR. The derived physical parameters provide information to decipher the nature of the dust cycle and of the stellar populations in these galaxies. The derived IR composite template is consistent with the galaxies in the studied sample. A delayed star formation history with tau_main = 500 Myrs is slightly favored by the statistical analysis as compared to a delayed with a final burst or a continuous star formation history. The position of the sample in the star formation rate (SFR)- M_star diagram is consistent with previous papers. The redshift evolution of the log M_star versus A_FUV relation is in agreement with evolution in the redshift of this relation. This evolution is necessary to explain the cosmic evolution of the average dust attenuation of galaxies. Evolution is also observed in the L_dust/ L_FUV (IRX) versus UV slope beta_FUV diagram: younger galaxies have bluer beta_FUV. We modeled the shift of galaxies in the IRX versus the beta_FUV diagram with the mass-weighted age as a free parameter, and we provide an equation to make predictions.

preprint2020arXiv

A hyper luminous starburst at z=4.72 magnified by a lensing galaxy pair at z=1.48

[Abridged] We discovered in the Herschel Reference Survey an extremely bright IR source with $S_{500}$~120mJy (Red Virgo 4 - RV4). Based on IRAM/EMIR and IRAM/NOEMA detections of the CO(5-4), CO(4-3), and [CI] lines, RV4 is located at z=4.724, yielding a total observed L$_{IR}$ of 1.1+/-0.6x0$^{14}$L$_{\odot}$. At the position of the Herschel emission, three blobs are detected with the VLA at 10cm. The CO(5-4) line detection of each blob confirms that they are at the same redshift with the same line width, indicating that they are multiple images of the same source. In Spitzer and deep optical observations, two sources, High-z Lens 1 (HL1) West and HL1 East, are detected at the center of the three VLA/NOEMA blobs. These two sources are placed at z=1.48 with XSHOOTER spectra, suggesting that they could be merging and gravitationally lensing the emission of RV4. HL1 is the second most distant lens known to date in strong lensing systems. The Einstein radius of the lensing system is 2.2"+/-0.2 (20kpc). The high redshift of HL1 and the large Einstein radius are highly unusual for a strong lensing system. We present the ISM properties of the background source RV4. Different estimates of the gas depletion time yield low values suggesting that RV4 is a SB galaxy. Among all high-z SMGs, this source exhibits one of the lowest L$_{[CI]}$ to L$_{IR}$ ratios, 3.2+/-0.9x10$^{-6}$, suggesting an extremely short gas tdepl of only 14+/-5Myr. It also shows a relatively high L$_{[CI]}$ to L$_{CO(4-3)}$ ratio (0.7+/-0.2) and low L$_{CO(5-4)}$ to L$_{IR}$ ratio (only ~50% of the value expected for normal galaxies) hinting a low density of gas. Finally, we discuss that the short tdepl of RV4 can be explained by either a very high SFE, which is difficult to reconcile with major mergers simulations of high-z galaxies, or a rapid decrease of SF, which would bias the estimate of tdepl toward low value.

preprint2020arXiv

A Virgo Environmental Survey Tracing Ionised Gas Emission. VESTIGE VIII. Bridging the cluster-ICM-galaxy evolution at small scales

We measure FIR emission from tails of stripped dust following the ionised and atomic gas components in galaxies undergoing ram pressure stripping. We study the dust-to-gas relative distribution and mass ratio in the stripped interstellar medium and relate them to those of the intra-cluster medium, thus linking the cluster-ICM-galaxy evolution at small-scales. The galaxy sample consists of three Scd Virgo galaxies with stellar masses in the range $10^9\lesssim \mathrm{M_{*}} \lesssim 10^{10}\, \mathrm{M_{\odot}}$, and within 1 Mpc from the cluster centre, namely NGC 4330, NGC 4522, and NGC 4654. Through the analysis of VESTIGE H$α$, $Herschel$ SPIRE far-infrared, and VIVA HI data, we trace the spatial distribution of the tails and infer the dust and gas masses from the measured far-infrared 250 $μ$m and HI flux densities. Dust-to-gas mass ratios (DGRs) in the tails are analysed as a function of the galaxy mass, metallicity, and dust temperature. Along the stripped component, the dust distribution closely follows the HI and H$α$ emitting gas, all extending beyond the optical disc. In these regions, the DGRs are $2.0\pm0.6\times10^{-3}$, $0.7\pm0.1\times10^{-3}$, and $0.4\pm0.03\times10^{-3}$, for NGC 4330, NGC 4522, and NGC 4654, respectively, i.e. up to a factor of 15 less than the values measured in the main body of nearby galaxies. We also find a negative trend in the DGR as a function of the metallicity that can be explained in terms of a dust component more centrally concentrated in more metal-rich systems. Together with the finding that the stripped dust is cold, $T_{d} \lesssim 25\, K$, our results support an outside-in stripping scenario of the galaxy interstellar medium. This study shows that ram pressure stripping is a key mechanism in the building up of the Virgo intra-cluster component injecting dust grains into the ICM, thus contributing to its metal enrichment.

preprint2020arXiv

Observational and theoretical constraints on the formation and early evolution of the first dust grains in galaxies at 5 < z < 10

The first generation of stars were born a few hundred million years after the big bang. These stars synthesized elements heavier than H and He, that are later expelled into the interstellar medium, initiating the rise of metals. Within this enriched medium, the first dust grains formed. This event is cosmological crucial for molecule formation as dust plays a major role by cooling low-metallicity star-forming clouds which can fragment to create lower mass stars. Collecting information on these first dust grains is difficult because of the negative alliance of large distances and low dust masses. We combine the observational information from galaxies at redshifts 5 < z < 10 to constrain their dust emission and theoretically understand the first evolutionary phases of the dust cycle. Spectral energy distributions (SEDs) are fitted with CIGALE and the physical parameters and their evolution are modelled. From this SED fitting, we build a dust emission template for this population of galaxies in the epoch of reionization. Our new models explain why some early galaxies are observed and others are not. We follow in time the formation of the first grains by supernovae later destroyed by other supernova blasts and expelled in the circumgalactic and intergalactic media. We have found evidence for the first dust grains formed in the universe. But, above all, this letter underlines the need to collect more data and to develop new facilities to further constrain the dust cycle in galaxies in the epoch of reionization.

preprint2020arXiv

X-CIGALE: fitting AGN/galaxy SEDs from X-ray to infrared

CIGALE is a powerful multiwavelength spectral energy distribution (SED) fitting code for extragalactic studies. However, the current version of CIGALE is not able to fit X-ray data, which often provide unique insights into AGN intrinsic power. We develop a new X-ray module for CIGALE, allowing it to fit SEDs from the X-ray to infrared (IR). We also improve the AGN fitting of CIGALE from UV-to-IR wavelengths. We implement a modern clumpy two-phase torus model, SKIRTOR. To account for moderately extincted type 1 AGNs, we implement polar-dust extinction. We publicly release the source code (named X-CIGALE). We test X-CIGALE with X-ray detected AGNs in SDSS, COSMOS, and AKARI-NEP. The fitting quality (as indicated by reduced $χ^2$) is good in general, indicating that X-CIGALE is capable of modelling the observed SED from X-ray to IR. We discuss constrainability and degeneracy of model parameters in the fitting of AKARI-NEP, for which excellent mid-IR photometric coverage is available. We also test fitting a sample of AKARI-NEP galaxies for which only X-ray upper limits are available from Chandra observations, and find that the upper limit can effectively constrain the AGN SED contribution for some systems. Finally, using X-CIGALE, we assess the ability of Athena to constrain the AGN activity in future extragalactic studies.

preprint2020arXiv

X-ray flux in the SED modelling: An application of X-CIGALE in the XMM-XXL field

X-CIGALE, built upon the spectral energy distribution (SED) code of CIGALE, implements important new features: The code accounts for obscuring material in the polars of the AGN and has the ability to fit X-ray fluxes. In this work, we use ~2500 spectroscopic, X-ray AGN from the XMM-XXL-North field and examine the improvements the new features bring in the SED modelling analysis. Based on our results, X-CIGALE successfully connects the X-ray with the UV luminosity in the whole range spanned by our sample (log Lx(2-10 keV) = (42 - 46) erg/s). The addition of the new features globally improves the efficiency of X-CIGALE in the estimation and characterization of the AGN component. The classification into type 1 and type 2 based on their inclination angle is improved, especially at redshifts lower than 1. The statistical significance of the AGN fraction, fracAGN, measurements is increased, in particular for luminous X-ray sources (log LX > 45 erg/s). These conclusions hold under the condition that (mid-) IR photometry is available in the SED fitting process. The addition of polar dust increases the AGN fraction and the efficiency of the SED decomposition to detect AGN among X-ray selected sources. X-CIGALE estimates a strong AGN (fracAGN > 0.3) in more than 90% of the infrared selected AGN and 75% of X-ray detected AGN not selected by IR colour criteria. The latter drops to ~50% when polar dust is not included. The ability of X-CIGALE to include X-ray information in the SED fitting process can be instrumental in the optimal exploitation of the wealth of data that current (eROSITA) and future (ATHENA) missions will provide us.

preprint2016arXiv

The FLARE mission: Deep and Wide-field 1-5$μ$m Imaging and Spectroscopy for the Early Universe: a proposal for M5 Cosmic Vision call

FLARE (First Light And Reionization Explorer) is a space mission that will be submitted to ESA (M5 call). Its primary goal (about 80% of lifetime) is to identify and study galaxies that dwell in the universe before the end of the reionization up to about z = 15, a redshift that might not be statistically reachable for JWST, Euclid and WFIRST. A secondary objective (about 20% of lifetime) is to survey star formation in the Milky Way. The strategy selected for FLARE optimizes the science return: imaging and spectroscopic integral-field observations will be carried out simultaneously on two parallel focal planes and over very wide instantaneous fields of view. FLARE will feature an instantaneous field of view of about 0.2 deg$^2$ with 0.2-arcsec pixels and an instantaneous integral-field spectroscopic field of view of about 1 arcmin$^2$ with R = 500 - 1000 and an angular resolution of about 0.4 arcsec. To detect first-light galaxies, the imaging and spectroscopic survey (parallel observations over about 6 years) will reach m$_{AB}$ = 28 and f$_λ$ = 10$^{-18}$ erg/cm$^2$/S. FLARE will help addressing two of ESA's Cosmic Vision themes: a) "How did the universe originate and what is it made of?" and b) "What are the conditions for planet formation and the emergence of life?" and more specifically, "From gas and dust to stars and planets". FLARE will provide to the ESA community a leading position to statistically study the early universe after JWST's deep but pin-hole surveys. Moreover, the instrumental development of wide-field imaging and wide-field integral-field spectroscopy in space will be a major breakthrough after making them available on ground-based telescopes.

preprint2016arXiv

The quenching of the star formation activity in cluster galaxies

We study the star formation quenching mechanism in cluster galaxies by fitting the SED of the Herschel Reference Survey, a complete volume-limited K-band-selected sample of nearby galaxies including objects in different density regions, from the core of the Virgo cluster to the general field. The SED are fitted using the CIGALE SED modelling code. The truncated activity of cluster galaxies is parametrised using a specific SFH with 2 free parameters, the quenching age QA and the quenching factor QF. These 2 parameters are crucial for the identification of the quenching mechanism which acts on long timescales if starvation while rapid and efficient if ram pressure. To be sensitive to an abrupt and recent variation of the star formation activity, we combine in a new way 20 UV to FIR photometric bands with 3 age-sensitive Balmer line absorption indices extracted from available medium-resolution integrated spectroscopy and with Halpha narrow band imaging data. The use of a truncated SFH significantly increases the quality of the fit in those objects whose atomic gas content has been removed during the interaction with the hostile cluster environment. The typical QA of the perturbed late-type galaxies is QA < 300 Myr whenever the activity of star formation is reduced by 50% < QF <= 80% and QA < 500 Myr for QF > 80%, while that of the quiescent early-types is QA ~ 1-3 Gyr. The fraction of late-types with a star formation activity reduced by QF > 80% and with an HI-deficiency parameter HI-def > 0.4 drops by a factor of ~ 5 from the inner half virial radius of the Virgo cluster, where the hot diffuse X-ray emitting gas of the cluster is located, to the outer regions. The efficient quenching of the star formation activity observed in Virgo suggests that the dominant stripping process is ram pressure. We discuss the implication of this result in the cosmological context of galaxy evolution.

preprint2015arXiv

Dust attenuation up to z~2 in the AKARI North Ecliptic Pole Deep Field

(Abridged) We aim to study the evolution of dust attenuation in galaxies selected in the IR in the redshift range in which they are known to dominate the star formation activity in the universe. The comparison with other measurements of dust attenuation in samples selected using different criteria will give us a global picture of the attenuation at work in star-forming galaxies and its evolution with redshift. Using multiple filters of IRC instrument, we selected more than 4000 galaxies from their rest-frame emission at 8 microns, from z~0.2 to 2$. We built SEDs from the rest-frame UV to the far-IR by adding data in the optical-NIR and from GALEX and Herschel surveys. We fit SEDs with the physically-motivated code CIGALE. We test different templates for AGNs and recipes for dust attenuation and estimate stellar masses, SFRs, amount of dust attenuation, and AGN contribution to the total IR luminosity. The AGN contribution to the total IR luminosity is found to be on average approximately 10% with a slight increase with redshift. Dust attenuation in galaxies dominating the IR luminosity function is found to increase from z=0 to z=1 and to remain almost constant from z=1 to z=1.5. Conversely, when galaxies are selected at a fixed IR luminosity, their dust attenuation slightly decreases as redshift increases but with a large dispersion. The attenuation in our mid-IR selected sample is found ~ 2 mag higher than that found globally in the universe or in UV and Halpha line selections in the same redshift range. This difference is well explained by an increase of dust attenuation with the stellar mass, in global agreement with other recent studies. Starbursting galaxies do not systematically exhibit a high attenuation

preprint2015arXiv

Dust properties of Lyman break galaxies at $z\sim3$

We explore from a statistical point of view the far-infrared (far-IR) and sub-millimeter (sub-mm) properties of a large sample of LBGs (22,000) at z~3 in the COSMOS field. The large number of galaxies allows us to split it in several bins as a function of UV luminosity, UV slope, and stellar mass to better sample their variety. We perform stacking analysis in PACS (100 and 160 um), SPIRE (250, 350 and 500 um) and AzTEC (1.1 mm) images. Our stacking procedure corrects the biases induced by galaxy clustering and incompleteness of our input catalogue in dense regions. We obtain the full IR spectral energy distributions (SED) of subsamples of LBGs and derive the mean IR luminosity as a function of UV luminosity, UV slope, and stellar mass. The average IRX is roughly constant over the UV luminosity range, with a mean of 7.9 (1.8 mag). However, it is correlated with UV slope, and stellar mass. We investigate using a statistically-controlled stacking analysis as a function of (stellar mass, UV slope) the dispersion of the IRX-UVslope and IRX-M* plane. Our results enable us to study the average relation between star-formation rate (SFR) and stellar mass, and we show that our LBG sample lies on the main sequence of star formation at z~3.

preprint2015arXiv

Environment of the submillimeter-bright massive starburst HFLS3 at $z\sim$6.34

We describe the search for Lyman-break galaxies (LBGs) near the sub-millimeter bright starburst galaxy HFLS3 at $z$$=$6.34 and a study on the environment of this massive galaxy during the end of reionization.We performed two independent selections of LBGs on images obtained with the \textit{Gran Telescopio Canarias} (GTC) and the \textit{Hubble Space Telescope} (HST) by combining non-detections in bands blueward of the Lyman-break and color selection. A total of 10 objects fulfilling the LBG selection criteria at $z$$>$5.5 were selected over the 4.54 and 55.5 arcmin$^2$ covered by our HST and GTC images, respectively. The photometric redshift, UV luminosity, and the star-formation rate of these sources were estimated with models of their spectral energy distribution. These $z$$\sim$6 candidates have physical properties and number densities in agreement with previous results. The UV luminosity function at $z$$\sim$6 and a Voronoi tessellation analysis of this field shows no strong evidence for an overdensity of relatively bright objects (m$_{F105W}$$<$25.9) associated with \textit{HFLS3}. However, the over-density parameter deduced from this field and the surface density of objects can not excluded definitively the LBG over-density hypothesis. Moreover we identified three faint objects at less than three arcseconds from \textit{HFLS3} with color consistent with those expected for $z$$\sim$6 galaxies. Deeper data are needed to confirm their redshifts and to study their association with \textit{HFLS3} and the galaxy merger that may be responsible for the massive starburst.

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

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 spectral energy distribution of galaxies at z > 2.5: Implications from the Herschel/SPIRE color-color diagram

We use the Herschel SPIRE color-color diagram to study the spectral energy distribution (SED) and the redshift estimation of high-z galaxies. We compiled a sample of 57 galaxies with spectroscopically confirmed redshifts and SPIRE detections in all three bands at $z=2.5-6.4$, and compared their average SPIRE colors with SED templates from local and high-z libraries. We find that local SEDs are inconsistent with high-z observations. The local calibrations of the parameters need to be adjusted to describe the average colors of high-z galaxies. For high-z libraries, the templates with an evolution from z=0 to 3 can well describe the average colors of the observations at high redshift. Using these templates, we defined color cuts to divide the SPIRE color-color diagram into different regions with different mean redshifts. We tested this method and two other color cut methods using a large sample of 783 Herschel-selected galaxies, and find that although these methods can separate the sample into populations with different mean redshifts, the dispersion of redshifts in each population is considerably large. Additional information is needed for better sampling.

preprint2014arXiv

HerMES: Spectral Energy Distributions of Submillimeter Galaxies at z >4

We present a study of the infrared properties for a sample of seven spectroscopically confirmed submillimeter galaxies at $z>$4.0. By combining ground-based near-infrared, Spitzer IRAC and MIPS, Herschel SPIRE, and ground-based submillimeter/millimeter photometry, we construct their Spectral Energy Distributions (SED) and a composite model to fit the SEDs. The model includes a stellar emission component at $λ_{\rm rest} <$ 3.5$ μ$m; a hot dust component peaking at $λ_{rest} \sim$ 5$\,μ$m; and cold dust component which becomes significant for $λ_{\rm rest} >$ 50$\,μ$m. Six objects in the sample are detected at 250 and 350$ μ$m. The dust temperatures for the sources in this sample are in the range of 40$-$80 K, and their $L_{\rm FIR}$ $\sim$ 10$^{13}$ L$_{\odot}$ qualifies them as Hyper$-$Luminous Infrared Galaxies (HyperLIRGs). The mean FIR-radio index for this sample is around $< q > = 2.2$ indicating no radio excess in their radio emission. Most sources in the sample have 24$ μ$m detections corresponding to a rest-frame 4.5$ μ$m luminosity of Log$_{10}$(L$_{4.5}$ / L$_{\odot}$) = 11 $\sim$ 11.5. Their L$_{\rm 4.5}$/$L_{\rm FIR}$ ratios are very similar to those of starburst dominated submillimeter galaxies at $z \sim$ 2. The $L_{\rm CO}-L_{\rm FIR}$ relation for this sample is consistent with that determined for local ULIRGs and SMGs at $z \sim$ 2. We conclude that submillimeter galaxies at $z >$ 4 are hotter and more luminous in the FIR, but otherwise very similar to those at $z \sim$ 2. None of these sources show any sign of the strong QSO phase being triggered.

preprint2014arXiv

HerMES: The Rest-Frame UV Emission and A Lensing Model for the z=6.34 Luminous Dusty Starburst Galaxy HFLS3

We discuss the restframe UV emission from the starbursting galaxy HFLS3 at z=6.34, discovered in Herschel/SPIRE data due to its red color in the submm wavelengths from 250-500 um. The apparent inst. SFR of HFLS3 inferred from the total FIR luminosity measured with over 15 photometric data points between 100 to 1000 um is 2900 Msun/yr. Keck/NIRC2 Ks band adaptive optics imaging data showed two potential NIR counterparts near HFLS3. Previously, the northern galaxy was taken to be in the foreground at z=2.1 while the southern galaxy was assumed to HFLS3's NIR counterpart. New HST/WFC3 and ACS imaging data show both optically bright galaxies are in the foreground at z<6. A new lensing model based on HST data and mm-wave continuum emission yields a magnification of 2.2+/-0.3. The lack of multiple imaging constrains the lensing magnification to be lower than either 2.7 or 3.5 at the 95% confidence level for the two scenarios, which attribute one or two components to HFLS3 in the source plane. Correcting for gravitational lensing, the inst. SFR is 1320 Msun/yr with the 95% confidence lower limit around 830 Msun/yr. Using models for the restframe UV to FIR SED, the ave. SFR over the last 100 Myr is around 660 Msun/yr. The dust and stellar masses of HFLS3 from the same SED models are 3x10^8 Msun and ~5x10^10 Msun, respectively, with large systematic uncertainties on assumptions related to the SED model. With HST/WFC3 images we also find diffuse NIR emission about 0.5" (~3 kpc) SW of HFLS3 that remains undetected in the ACS data. The emission has a photometric redshift consistent with either z~6 or a dusty galaxy template at z~2. If at the same redshift as HFLS3 the detected diffuse emission could be part of the complex merger system that could be triggering the starburst. Alternatively, it could be part of the foreground structure at z~2.1 that is responsible for lensing of HFLS3.

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

Dusty Universe viewed by AKARI far infrared detector

We present the results of the analysis of multiwavelength Spectral Energy Distributions (SEDs) of far-infrared galaxies detected in the AKARI Deep Field-South (ADF--S) Survey. The analysis uses a carefully selected sample of 186 sources detected at the 90 $μ$m AKARI band, identified as galaxies with cross-identification in public catalogues. For sources without known spectroscopic redshifts, we estimate photometric redshifts after a test of two independent methods: one based on using mainly the optical -- mid infrared range, and one based on the whole range of ultraviolet -- far infrared data. We observe a vast improvement in the estimation of photometric redshifts when far infrared data are included, compared with an approach based mainly on the optical -- mid infrared range. We discuss the physical properties of our far-infrared-selected sample. We conclude that this sample consists mostly of rich in dust and young stars nearby galaxies, and, furthermore, that almost 25% of these sources are (Ultra)Luminous Infrared Galaxies. Average SEDs normalized at 90 $μ$m for normal galaxies (138 sources), LIRGs (30 sources), and ULIRGs (18 galaxies) a the significant shift in the peak wavelength of the dust emission, and an increasing ratio between their bolometric and dust luminosities which varies from 0.39 to 0.73.

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: dust attenuation and star formation activity in UV-selected samples from z~4 to z~1.5

We study the link between observed ultraviolet luminosity, stellar mass, and dust attenuation within rest-frame UV-selected samples at z~ 4, 3, and 1.5. We measure by stacking at 250, 350, and 500 um in the Herschel/SPIRE images from the HerMES program the average infrared luminosity as a function of stellar mass and UV luminosity. We find that dust attenuation is mostly correlated with stellar mass. There is also a secondary dependence with UV luminosity: at a given UV luminosity, dust attenuation increases with stellar mass, while at a given stellar mass it decreases with UV luminosity. We provide new empirical recipes to correct for dust attenuation given the observed UV luminosity and the stellar mass. Our results also enable us to put new constraints on the average relation between star formation rate and stellar mass at z~ 4, 3, and 1.5. The star formation rate-stellar mass relations are well described by power laws (SFR~ M^0.7), with the amplitudes being similar at z~4 and z~3, and decreasing by a factor of 4 at z~1.5 at a given stellar mass. We further investigate the evolution with redshift of the specific star formation rate. Our results are in the upper range of previous measurements, in particular at z~3, and are consistent with a plateau at 3<z<4. Current model predictions (either analytic, semi-analytic or hydrodynamic) are inconsistent with these values, as they yield lower predictions than the observations in the redshift range we explore. We use these results to discuss the star formation histories of galaxies in the framework of the Main Sequence of star-forming galaxies. Our results suggest that galaxies at high redshift (2.5<z<4) stay around 1 Gyr on the Main Sequence. With decreasing redshift, this time increases such that z=1 Main Sequence galaxies with 10^8<M_*/Msun<10^10 stay on the Main Sequence until z=0.

preprint2013arXiv

Mid- to far infrared properties of star-forming galaxies and active galactic nuclei

We study the mid- to far-IR properties of a 24um-selected flux-limited sample (S24 > 5mJy) of 154 intermediate redshift (<z>~0.15), infrared luminous galaxies, drawn from the 5MUSES survey. By combining existing mid-IR spectroscopy and new Herschel SPIRE submm photometry from the HerMES program, we derived robust total infrared luminosity (LIR) and dust mass (Md) estimates and infered the relative contribution of the AGN to the infrared energy budget of the sources. We found that the total infrared emission of galaxies with weak 6.2um PAH emission (EW<0.2um) is dominated by AGN activity, while for galaxies with EW>0.2um more than 50% of the LIR arises from star formation. We also found that for galaxies detected in the 250-500um Herschel bands an AGN has a statistically insignificant effect on the temperature of the cold dust and the far-IR colours of the host galaxy, which are primarily shaped by star formation activity. For star-forming galaxies we reveal an anti-correlation between the LIR-to-rest-frame 8um luminosity ratio, IR8 = LIR\L8, and the strength of PAH features. We found that this anti-correlation is primarily driven by variations in the PAHs emission, and not by variations in the 5-15um mid-IR continuum emission. Using the [NeIII]/[NeII] line ratio as a tracer of the hardness of the radiation field, we confirm that galaxies with harder radiation fields tend to exhibit weaker PAH features, and found that they have higher IR8 values and higher dust-mass-weighted luminosities (LIR/Md), the latter being a proxy for the dust temperature (Td). We argue that these trends originate either from variations in the environment of the star-forming regions or are caused by variations in the age of the starburst. Finally, we provide scaling relations that will allow estimating LIR, based on single-band observations with the mid-infrared instrument, on board the upcoming JWST.

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

Properties of star forming galaxies in AKARI Deep Field-South

The main aim of this work is the characterization of physical properties of galaxies detected in the far infrared (FIR) in the AKARI Deep Field-South (ADF-S) survey. Starting from a catalog of the 1 000 brightest ADF-S sources in the WIDE-S (90$μ$m) AKARI band, we constructed a subsample of galaxies with spectral coverage from the ultraviolet to the far infrared. We then analyzed the multiwavelength properties of this 90$μ$m selected sample of galaxies. For galaxies without known spectroscopic redshifts we computed photometric redshifts using the codes Photometric Analysis for Redshift Estimate (Le PHARE) and Code Investigating GALaxy Emission (CIGALE), tested these photometric redshifts using spectroscopic redshifts, and compared the performances of both codes. To test the reliability of parameters obtained by fitting spectral energy distributions, a mock cataloge was generated. We built a large multiwavelength catalog of more than 500 ADF-S galaxies. We successfully fitted Spectral Energy Distributions of 186 galaxies with $\rm{χ^2_{min}<4}$, and analyzed the output parameters of the fits. We conclude that our sample consists mostly of nearby actively star-forming galaxies, and all our galaxies have a relatively high metallicity. We estimated photometric redshifts for 113 galaxies from the whole ADF-S sample. Comparing the performance of Le PHARE and CIGALE, we found that CIGALE gives more reliable redshift estimates for our galaxies, which implies that including the IR photometry allows for substantial improvement of photometric redshift estimation.

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

Ultraviolet to infrared emission of z>1 galaxies: Can we derive reliable star formation rates and stellar masses?

We seek to derive star formation rates (SFR) and stellar masses (M_star) in distant galaxies and to quantify the main uncertainties affecting their measurement. We explore the impact of the assumptions made in their derivation with standard calibrations or through a fitting process, as well as the impact of the available data, focusing on the role of IR emission originating from dust. We build a sample of galaxies with z>1, all observed from the UV to the IR (rest frame). The data are fitted with the code CIGALE, which is also used to build and analyse a catalogue of mock galaxies. Models with different SFHs are introduced. We define different set of data, with or without a good sampling of the UV range, NIR, and thermal IR data. The impact of these different cases on the determination of M_star and SFR are analysed. Exponentially decreasing models with a redshift formation of the stellar population z ~8 cannot fit the data correctly. The other models fit the data correctly at the price of unrealistically young ages when the age of the single stellar population is taken to be a free parameter. The best fits are obtained with two stellar populations. As long as one measurement of the dust emission continuum is available, SFR are robustly estimated whatever the chosen model is, including standard recipes. M_star measurement is more subject to uncertainty, depending on the chosen model and the presence of NIR data, with an impact on the SFR-M_star scatter plot. Conversely, when thermal IR data from dust emission are missing, the uncertainty on SFR measurements largely exceeds that of stellar mass. Among all physical properties investigated here, the stellar ages are found to be the most difficult to constrain and this uncertainty acts as a second parameter in SFR measurements and as the most important parameter for M_star measurements.

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

An investigation of star formation and dust attenuation in major mergers using ultraviolet and infrared data

Merger processes play an important role in galaxy formation and evolution. To study the influence of merger processes on the evolution of dust properties and cosmic star formation rate, we investigate a local sample of major merger galaxies and a control sample of isolated galaxies using GALEX ultraviolet (UV) and Spitzer infrared (IR) images. Through a statistical study, we find that dust attenuation in merger galaxies is enhanced with respect to isolated galaxies. We find this enhancement is contributed mainly by spiral galaxies in spiral-spiral (S-S) pairs, and increases with the increasing stellar mass of a galaxy. Combining the IR and UV parts of star formation rates (SFRs), we then calculated the total SFRs and specific star formation rates (SSFRs). We find the SSFRs to be enhanced in merger galaxies. This enhancement depends on galaxy stellar mass and the companion's morphology, but depends little on whether the galaxy is a primary or secondary component or on the separation between two components. These results are consistent with a previous study based only on IR images. In addition, we investigate the nuclear contributions to SFRs. SFRs in paired galaxies are more concentrated in the central part of the galaxies than in isolate galaxies. Our studies of dust attenuation show that the nuclear parts of pairs most resemble ULIRGs. Including UV data in the present work not only provides reliable information on dust attenuation, but also refines analyses of SFRs.

preprint2012arXiv

GOODS-Herschel:dust attenuation properties of UV selected high redshift galaxies

We study dust attenuation at UV wavelengths at high redshift, where the UV is redshifted to the observed visible. In particular, we search for a UV bump and related implications for dust attenuation determinations. We use data in the CDFS, obtained in intermediate and broad band filters by the MUSYC project, to sample the UV rest-frame of 751 galaxies with 0.95<z<2.2. When available, Herschel/PACS data (GOODS-Herschel project), and Spitzer/MIPS measurements, are used to estimate the dust emission. The SED of each source is fit using the CIGALE code. The amount of dust attenuation and the dust attenuation curve are obtained as outputs of the SED fitting process, together with other parameters linked to the SFH. The global amount of dust attenuation at UV wavelengths is found to increase with stellar mass and to decrease as UV luminosity increases. A UV bump at 2175A is securely detected in 20% of the galaxies, and the mean amplitude of the bump for the sample is similar to that observed in the LMC supershell region. This amplitude is found to be lower in galaxies with very high SSFRs, and 90% of the galaxies exhibiting a secure bump are at z<1.5. The attenuation curve is confirmed to be steeper than that of local starburst galaxies for 20$% of the galaxies. The large dispersion found for these two parameters describing the attenuation law is likely to reflect a wide diversity of attenuation laws among galaxies. The relations between dust attenuation, IR-to-UV flux ratio, and the slope of the UV continuum are derived for the mean attenuation curve found for our sample. Deviations from the average trends are found to correlate with the age of the young stellar population and the shape of the attenuation curve.(abriged)

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

Star Formation and Dust Extinction Properties of Local Galaxies as seen from AKARI and GALEX

An accurate estimation of the star formation-related properties of galaxies is crucial for understanding the evolution of galaxies. In galaxies, ultraviolet (UV) light emitted by recently formed massive stars is attenuated by dust, which is also produced by star formation (SF) activity, and is reemitted at mid- and far- infrared (IR) wavelengths. In this study, we investigate the star formation rate (SFR) and dust extinction using UV and IR data. We selected local galaxies which are detected at AKARI FIS 90 um and matched the IRAS IIFSCz 60 um select catalog. We measured FUV and NUV flux densities from GALEX images. We examined the SF and extinction of Local galaxies using four bands of AKARI. Then, we calculated FUV and total IR luminosities, and obtained the SF luminosity, L_{SF}, the total luminosity related to star formation activity, and the SFR. We find that in most galaxies, L_{SF} is dominated by L_{dust}. We also find that galaxies with higher SF activity have a higher fraction of their SF hidden by dust. In fact, the SF of galaxies with SFRs >20 M_{sun}/yr is almost completely hidden by dust. Our results boast a significantly higher precision with respect to previously published works, due to the use of much larger object samples from the AKARI and GALEX all sky surveys.

preprint2012arXiv

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

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

preprint2012arXiv

The Herschel Multi-tiered Extragalactic Survey: HerMES

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

preprint2011arXiv

AKARI/IRC Broadband Mid-infrared data as an indicator of Star Formation Rate

AKARI/Infrared Camera (IRC) Point Source Catalog provides a large amount of flux data at {\it S9W} ($9\ {\rm μm}$) and {\it L18W} ($18\ {\rm μm}$) bands. With the goal of constructing Star-Formation Rate(SFR) calculations using IRC data, we analyzed an IR selected GALEX-SDSS-2MASS-AKARI(IRC/Far-Infrared Surveyor) sample of 153 nearby galaxies. The far-infrared fluxes were obtained from AKARI diffuse maps to correct the underestimation for extended sources raised by the point-spread function photometry. SFRs of these galaxies were derived by the spectral energy distribution fitting program CIGALE. In spite of complicated features contained in these bands, both the {\it S9W} and {\it L18W} emission correlate with the SFR of galaxies. The SFR calibrations using {\it S9W} and {\it L18W} are presented for the first time. These calibrations agree well with previous works based on Spitzer data within the scatters, and should be applicable to dust-rich galaxies.

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

GOODS-Herschel A population of 24micron dropout sources at z < 2

Using extremely deep PACS 100- and 160um Herschel data from the GOODS-Herschel program, we identify 21 infrared bright galaxies previously missed in the deepest 24um surveys performed by MIPS. These MIPS dropouts are predominantly found in two redshift bins, centred at z ~0.4 and ~1.3. Their S_100/S_24 flux density ratios are similar to those of local LIRGs and ULIRGs, whose silicate absorption features at 18um (at z ~ 0.4) and 9.7um (at z ~ 1.3) are shifted into the 24um MIPS band at these redshifts. The high-z sub-sample consists of 11 infrared luminous sources, accounting for ~2% of the whole GOODS-Herschel sample and putting strong upper limits on the fraction of LIRGs/ULIRGs at 1.0<z<1.7 that are missed by the 24um surveys. We find that a S_100/S_24 > 43 colour cut selects galaxies with a redshift distribution similar to that of the MIPS dropouts and when combined with a second colour cut, S_16/S_8 > 4, isolates sources at 1.0 < z < 1.7. We show that these sources have elevated specific star formation rates (sSFR) compared to main sequence galaxies at these redshifts and are likely to be compact starbursts with moderate/strong 9.7um silicate absorption features in their mid-IR spectra. Herschel data reveal that their infrared luminosities extrapolated from the 24um flux density are underestimated, on average, by a factor of ~3. These silicate break galaxies account for 16% (8%) of the ULIRG (LIRG) population in the GOODS fields, indicating a lower limit in their space density of 2.0 \times 10^(-5) Mpc^(-3). Finally, we provide estimates of the fraction of z < 2 MIPS dropout sources as a function of the 24-, 100-, 160-, 250- and 350um sensitivity limits, and conclude that previous predictions of a population of silicate break galaxies missed by the major 24um extragalactic surveys have been overestimated.

preprint2011arXiv

GOODS-Herschel Measurements of the Dust Attenuation of Typical Star-Forming Galaxies at High Redshift: Observations of UV-Selected Galaxies at z~2

We take advantage of the sensitivity and resolution of Herschel at 100 and 160 micron to directly image the thermal dust emission and investigate the infrared luminosities, L(IR), and dust obscuration of typical star-forming (L*) galaxies at high redshift. Our sample consists of 146 UV-selected galaxies with spectroscopic redshifts 1.5<z<2.6 in the GOODS-North field. Supplemented with deep Very Large Array (VLA) and Spitzer imaging, we construct median stacks at the positions of these galaxies at 24, 100, and 160 micron, and 1.4 GHz. The comparison between these stacked fluxes and a variety of dust templates and calibrations implies that typical star-forming galaxies with UV luminosities L(UV)>1e10 Lsun at z~2 are luminous infrared galaxies (LIRGs) with a median L(IR)=(2.2+/-0.3)e11 Lsun. Typical galaxies at 1.5<z<2.6 have a median dust obscuration L(IR)/L(UV) = 7.1+/-1.1, which corresponds to a dust correction factor, required to recover the bolometric star formation rate (SFR) from the unobscured UV SFR, of 5.2+/-0.6. This result is similar to that inferred from previous investigations of the UV, H-alpha, 24 micron, radio, and X-ray properties of the same galaxies studied here. Stacking in bins of UV slope implies that L* galaxies with redder spectral slopes are also dustier, and that the correlation between UV slope and dustiness is similar to that found for local starburst galaxies. Hence, the rest-frame 30 and 50 micron fluxes validate on average the use of the local UV attenuation curve to recover the dust attenuation of typical star-forming galaxies at high redshift. In the simplest interpretation, the agreement between the local and high redshift UV attenuation curves suggests a similarity in the dust production and stellar and dust geometries of starburst galaxies over the last 10 billion years.

preprint2011arXiv

GOODS-Herschel: an infrared main sequence for star-forming galaxies

We present the deepest far-IR observations obtained with Herschel and examine the 3-500um SEDs of galaxies at 0<z<2.5, supplemented by a local reference sample from IRAS, ISO, Spitzer and AKARI data. We find that the ratio of total IR luminosity to rest-frame 8um luminosity, IR8 (=Lir/L8), follows a Gaussian distribution centered on IR8=4 and defines an IR main sequence (MS). A minority population (<20 %) of outliers producing a tail skewed toward higher values of IR8 consist of starbursts (SB) with compact projected star formation densities. IR8 can be used to separate galaxies with normal and extended modes of star formation from compact SBs with high-IR8, high projected IR surface brightness (>3x10^10 Lsun kpc^-2) and a high specific SFR (i.e., SBs). The rest-frame, UV-2700A size of these distant SBs is typically half that of MS galaxies, supporting the correlation between star formation density and SB activity that is measured for the local sample. Locally, (U)LIRGs are systematically in the SB mode, whereas most distant (U)LIRGs form stars in the "normal" MS mode. This confusion between two modes of star formation is the cause of the so-called "mid-IR excess" population of galaxies found at z>1.5 by previous studies. MS galaxies have strong PAH emission line features, a broad far-IR bump resulting from a combination of dust temperatures (Tdust~15-50 K), and an effective Tdust~31 K, as derived from the peak wavelength of their IR SED. Galaxies in the SB regime instead exhibit weak PAH EW and a sharper far-IR bump with an effective Tdust~40 K. Finally, we present evidence that the mid-to-far IR emission of X-ray AGNs is predominantly produced by star formation and that candidate dusty AGNs with a power-law emission in the mid-IR systematically occur in compact, dusty SBs. After correcting for the effect of SBs on IR8, we identify new candidates for extremely obscured AGNs.

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/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

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

Spectral energy distributions of an AKARI-SDSS-GALEX sample of galaxies

The nearby universe remains the best laboratory to understand physical properties of galaxies and is a reference for any comparison with high redshift observations. The all sky (or very large) surveys performed from the ultraviolet (UV) to the far-infrared (far-IR) gives us large datasets with a very large wavelength coverage to perform a reference study. We want to investigate dust attenuation characteristics as well as star formation rate (SFR) calibrations on a sample of nearby galaxies observed over 13 bands from 0.15 to 160 microns. A sample of 363 galaxies is built from the AKARI /FIS all sky survey cross-correlated with SDSS and GALEX surveys. Broad band spectral energy distributions are fitted with the CIGALE code optimized to analyze variations in dust attenuation curves and SFR measurements and based on an energetic budget between the stellar and dust emission. Our galaxy sample is primarily selected in far-IR and mostly constituted of massive, actively star forming galaxies. There is some evidence for a dust attenuation law slightly steeper than that used for starburst galaxies but we are unable to constrain the presence or not of a bump at 220 nm. We confirm that a time dependent dust attenuation is necessary to perform the best fits. Various calibrations of the dust attenuation in the UV as a function of UV-optical colours are discussed. A calibration of the current star formation rate combining UV and total IR emissions is proposed with an accurate estimation of dust heating by old stars: for the whole sample 17 % of the total dust luminosity is not related to the recent star formation

preprint2011arXiv

Spitzer Imaging of Herschel-ATLAS Gravitationally Lensed Submillimeter Sources

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

preprint2011arXiv

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

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

preprint2010arXiv

A search for debris disks in the Herschel ATLAS

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

preprint2010arXiv

Detection of the Cosmic Far-Infrared Background in the AKARI Deep Field South

We report the detection and measurement of the absolute brightness and spatial fluctuations of the cosmic infrared background (CIB) with the AKARI satellite. We have carried out observations at 65, 90, 140 and 160 um as a cosmological survey in AKARI Deep Field South (ADF-S), which is one of the lowest cirrus regions with contiguous area on the sky. After removing bright galaxies and subtracting zodiacal and Galactic foregrounds from the measured sky brightness, we have successfully measured the CIB brightness and its fluctuations across a wide range of angular scales from arcminutes to degrees. The measured CIB brightness is consistent with previous results reported from COBE data but significantly higher than the lower limits at 70 and 160 um obtained with the Spitzer satellite from the stacking analysis of 24-um selected sources. The discrepancy with the Spitzer result is possibly due to a new galaxy population at high redshift obscured by hot dust. From power spectrum analysis at 90 um, three components are identified: shot noise due to individual galaxies; Galactic cirrus emission dominating at the largest angular scales of a few degrees; and an additional component at an intermediate angular scale of 10-30 arcminutes, possibly due to galaxy clustering. The spectral shape of the clustering component at 90 um is very similar to that at longer wavelengths as observed by Spitzer and BLAST. Moreover, the color of the fluctuations indicates that the clustering component is as red as Ultra-luminous infrared galaxies (ULIRGs) at high redshift, These galaxies are not likely to be the majority of the CIB emission at 90 um, but responsible for the clustering component. Our results provide new constraints on the evolution and clustering properties of distant infrared galaxies.

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

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

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

preprint2010arXiv

H-ATLAS: PACS imaging for the Science Demonstration Phase

We describe the reduction of data taken with the PACS instrument on board the Herschel Space Observatory in the Science Demonstration Phase of the Herschel-ATLAS (H-ATLAS) survey, specifically data obtained for a 4x4-deg^2 region using Herschel's fast-scan (60 arcsec/s) parallel mode. We describe in detail a pipeline for data reduction using customised procedures within HIPE from data retrieval to the production of science-quality images. We found that the standard procedure for removing Cosmic-Ray glitches also removed parts of bright sources and so implemented an effective two-stage process to minimise these problems. The pronounced 1/f noise is removed from the timelines using 3.4- and 2.5-arcmin boxcar high-pass filters at 100 and 160-um. Empirical measurements of the point-spread function (PSF) are used to determine the encircled energy fraction as a function of aperture size. For the 100- and 160-um bands, the effective PSFs are ~9 and ~13 arcsec (FWHM), and the 90-per-cent encircled energy radii are 13 and 18 arcsec. Astrometric accuracy is good to ~<2 arcsec. The noise in the final maps is correlated between neighbouring pixels and rather higher than advertised prior to launch. For a pair of cross-scans, the 5-sigma point-source sensitivities are 125-165 mJy for 9-13-arcsec radius apertures at 100-um and 150-240 mJy for 13-18-arcsec radius apertures at 160-um.

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: 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: 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 ATLAS: The angular correlation function of submillimetre galaxies at high and low redshift

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

preprint2010arXiv

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

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

preprint2010arXiv

Herschel 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-ATLAS: blazars in the SDP field

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

preprint2010arXiv

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

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

preprint2010arXiv

Herschel-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

Lyman `bump' galaxies - II. A possible signature of massive extremely metal-poor or metal-free stars in z=3.1 Ly-alpha emitters

(Abridged) Deep NB359 imaging with Subaru by Iwata et al. have detected surprisingly strong Lyman continuum (LyC; ~900A in the rest-frame) from some LAEs at z=3.1. However, the redshifts might be misidentified due to a narrow wavelength coverage in previous spectroscopy. We here present new deep spectroscopy covering the observed 4,000-7,000A with VLT/VIMOS and Subaru/FOCAS of 8 LAEs detected in NB359. All the 8 objects have only one detectable emission line around 4,970A which is most likely to be Ly-A at z=3.1, and thus, the objects are certainly LAEs at the redshift. However, 5 of them show a ~0.''8 spatial offset between the Ly-A emission and the source detected in NB359. No indications of the redshifts of the NB359 sources are found although it is statistically difficult that all the 5 LAEs have a foreground object accounting for the NB359 flux. The rest 3 LAEs show no significant offset from the NB359 position. Therefore, they are truly LyC emitting LAEs at z=3.1. We also examine the stellar population which simultaneously accounts for the strength of the LyC and the spectral slope of non-ionizing ultraviolet of the LAEs. We consider the latest statistics of Lyman limit systems to estimate the LyC optical depth in the IGM and an additional contribution of the bound-free LyC from photo-ionized nebulae to the LyC emissivity. As a result, we find that stellar populations with metallicity Z>=1/50Z_sun can explain the observed LyC strength only with a very top-heavy initial mass function (IMF; <m>~50 M_sun). However, the critical metallicity for such an IMF is expected to be much lower. A very young (~1 Myr) and massive (~100 M_sun) extremely metal-poor (Z<=5e-4Z_sun) or metal-free (so-called Population III) stellar population can reproduce the observed LyC strength. The required mass fraction of such `primordial' stellar population is ~1--10% in total stellar mass of the LAEs.

preprint2010arXiv

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

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

preprint2010arXiv

The far-infrared/radio correlation as probed by Herschel

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

preprint2010arXiv

The HerMES SPIRE submillimeter local luminosity function

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

preprint2010arXiv

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

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

preprint2010arXiv

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

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

preprint2010arXiv

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

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

preprint2010arXiv

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

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

preprint2009arXiv

Analysis of galaxy SEDs from far-UV to far-IR with CIGALE: Studying a SINGS test sample

Photometric data of galaxies covering the rest-frame wavelength range from far-UV to far-IR make it possible to derive galaxy properties with a high reliability by fitting the attenuated stellar emission and the related dust emission at the same time. For this purpose we wrote the code CIGALE (Code Investigating GALaxy Emission) that uses model spectra composed of the Maraston (or PEGASE) stellar population models, synthetic attenuation functions based on a modified Calzetti law, spectral line templates, the Dale & Helou dust emission models, and optional spectral templates of obscured AGN. Depending on the input redshifts, filter fluxes are computed for the model set and compared to the galaxy photometry by carrying out a Bayesian-like analysis. CIGALE was tested by analysing 39 nearby galaxies selected from SINGS. The reliability of the different model parameters was evaluated by studying the resulting expectation values and their standard deviations in relation to the input model grid. Moreover, the influence of the filter set and the quality of photometric data on the code results was estimated. For up to 17 filters between 0.15 and 160 mum, we find robust results for the mass, star formation rate, effective age of the stellar population at 4000 A, bolometric luminosity, luminosity absorbed by dust, and attenuation in the far-UV. A study of the mutual relations between the reliable properties confirms the dependence of star formation activity on morphology in the local Universe and indicates a significant drop in this activity at about 10^11 M_sol towards higher total stellar masses. The dustiest sample galaxies are present in the same mass range. [abridged]

preprint2009arXiv

Detections of Lyman Continuum from Star-forming Galaxies at z~3 Through Subaru/Suprime-Cam Narrow-band Imaging

Knowing the amount of ionizing photons from young star-forming galaxies is of particular importance to understanding the reionization process. Here we report initial results of Subaru/Suprime-Cam deep imaging observation of the SSA22 proto-cluster region at z=3.09, using a special narrow-band filter to optimally trace Lyman continuum (LyC) from galaxies at z~3. The unique wide field-of-view of Suprime-Cam enabled us to search for ionizing photons from 198 galaxies (73 Lyman break galaxies (LBGs) and 125 Ly-alpha emitters (LAEs)) with spectroscopically measured redshifts z~3.1. We detected LyC from 7 LBGs, as well as from 10 LAE candidates. Some of the detected galaxies show significant spatial offsets of LyC from non-ionizing UV emission. For some LBGs the observed non-ionizing UV to LyC flux density ratios are smaller than values expected from population synthesis models with a standard Salpeter initial mass function (IMF) with moderate dust attenuation (which is suggested from the observed UV slopes), even if we assume very transparent IGM along the sightlines of these objects. This implies an intrinsically bluer spectral energy distribution, e.g, that produced by a top-heavy IMF, for these LBGs. The observed flux desity ratios of non-ionizing UV to LyC of 7 detected LBGs range from 2.4 to 23.8 and the median is 6.6. The observed flux density ratios of the detected LAEs are even smaller than LBGs, if they are truly at z~3.1. We find that the median value of the flux density ratio for the deteced LBGs suggest that their escape fractions is likely to be higher than 4%, if the LyC escape is isotropic. The results imply that some of the LBGs in the proto-cluster at z~3 have the escape fraction significantly higher than that of galaxies (in a general field) at z~1 studied previously.

preprint2009arXiv

Star formation and dust extinction properties of local galaxies from AKARI-GALEX All-Sky Surveys: First results from most secure multiband sample from FUV to FIR

The AKARI All-Sky Survey provided the first bright point source catalog detected at 90um. Starting from this catalog, we selected galaxies by matching AKARI sources with those in the IRAS PSCz. Next, we have measured total GALEX FUV and NUV flux densities. Then, we have matched this sample with SDSS and 2MASS galaxies. By this procedure, we obtained the final sample which consists of 607 galaxies. If we sort the sample with respect to 90um, their average SED shows a coherent trend: the more luminous at 90um, the redder the global SED becomes. The M_r--NUV-r color-magnitude relation of our sample does not show bimodality, and the distribution is centered on the green valley between the blue cloud and red sequence seen in optical surveys. We have established formulae to convert FIR luminosity from AKARI bands to the total infrared (IR) luminosity L_TIR. With these formulae, we calculated the star formation directly visible with FUV and hidden by dust. The luminosity related to star formation activity (L_SF) is dominated by L_TIR even if we take into account the far-infrared (FIR) emission from dust heated by old stars. At high star formation rate (SFR) (> 20 Msun yr^-1), the fraction of directly visible SFR, SFR_FUV, decreases. We also estimated the FUV attenuation A_FUV from FUV-to-total IR (TIR) luminosity ratio. We also examined the L_TIR/L_FUV-UV slope (FUV- NUV) relation. The majority of the sample has L_TIR/L_FUV ratios 5 to 10 times lower than expected from the local starburst relation, while some LIRGs and all the ULIRGs of this sample have higher L_TIR/L_FUV ratios. We found that the attenuation indicator L_TIR/L_FUV is correlated to the stellar mass of galaxies, M*, but there is no correlation with specific SFR (SSFR), SFR/M*, and dust attenuation L_TIR/L_FUV. (abridged)

preprint2009arXiv

The Herschel ATLAS

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

preprint2009arXiv

The infrared emission of ultraviolet selected galaxies from z = 0 to z=1

We select galaxies in UV rest-frame at z=0, z~0.7 and z~1 together with a sample of LBGs at z~1, the samples are built in order to sample the same range of luminosity at any redshift. The evolution of the IR and UV luminosities with z is analysed for individual galaxies as well as in terms of luminosity functions. The L_IR/L_UV ratio is used to measure dust attenuation. This ratio does not exhibit a strong evolution with z for the bulk of our sample galaxies but some trends are found for galaxies with a strong dust attenuation and for UV luminous sources: galaxies with L_IR/L_UV>10 are more frequent at z>0 than at z=0 and the largest values of L_IR/L_UV are found for UV faint objects; conversely the most luminous galaxies of our samples (L_UV> 2 10^{10} L_sun$), detected at z=1, exhibit a lower dust attenuation than the fainter ones. L_IR/L_UV increases with the K rest-frame luminosity of the galaxies at all the redshifts considered and shows a residual anti-correlation with L_UV. The most massive and UV luminous galaxies exhibit quite large specific star formation rates. LBGs exhibit systematically lower dust attenuation than UV selected galaxies of same luminosity but similar specific star formation rates. The analysis of the UV+IR luminosity functions leads to the conclusion that up to z = 1 most of the star formation activity of UV selected galaxies is emitted in IR. Whereas we are able to retrieve all the star formation from our UV selection at z=0.7, at z = 1 we miss a large fraction of galaxies more luminous than ~ 10^{11} L_sun. The effect is found larger for Lyman Break Galaxies.

preprint2009arXiv

The Science Case for PILOT I: Summary and Overview

PILOT (the Pathfinder for an International Large Optical Telescope) is a proposed 2.5 m optical/infrared telescope to be located at Dome C on the Antarctic plateau. Conditions at Dome C are known to be exceptional for astronomy. The seeing (above ~30 m height), coherence time, and isoplanatic angle are all twice s good as at typical mid-latitude sites, while the water-vapour column, and the atmosphere and telescope thermal emission are all an order of magnitude better. These conditions enable a unique scientific capability for PILOT, which is addressed in this series of papers. The current paper presents an overview of the optical and instrumentation suite for PILO and its expected performance, a summary of the key science goals and observational approach for the facility, a discussion of the synergies between the science goals for PILOT and other telescopes, and a discussion of the future of Antarctic astronomy. Paper II and Paper III present details of the science projects divided, respectively, between the distant Universe (i.e., studies of first light, and the assembly and evolution of structure) and the nearby Universe (i.e., studies of Local Group galaxies, the Milky Way, and the Solar System).

preprint2009arXiv

The Science Case for PILOT II: the Distant Universe

PILOT (the Pathfinder for an International Large Optical Telescope) is a proposed 2.5 m optical/infrared telescope to be located at Dome C on the Antarctic plateau. The atmospheric conditions at Dome C deliver a high sensitivity, high photometric precision, wide-field, high spatial resolution, and high-cadence imaging capability to the PILOT telescope. These capabilities enable a unique scientific potential for PILOT, which is addressed in this series of papers. The current paper presents a series of projects dealing with the distant (redshift >) Universe, that have been identified as key science drivers for the PILOT facility. The potential for PILOT to detect the first populations of stars to form in the early Universe, via infrared projects searching for pair-instability supernovae and gamma-ray burst afterglows, is investigated. Two projects are proposed to examine the assembly and evolution of structure in the Universe: an infrared survey searching for the first evolved galaxies at high redshift, and an optical survey aimed at characterising moderate-redshift galaxy clusters. Finally, a large-area weak-lensing survey and a program to obtain supernovae infrared light-curves are proposed to examine the nature and evolution of dark energy and dark matter.

preprint2004arXiv

Dust attenuation in the nearby Universe: comparison between galaxies selected in the ultraviolet or in the infrared

We compare the dust attenuation properties of two samples of galaxies purely selected in the near-ultraviolet (NUV) band (1750-2750 A, lambda_m = 2310 A) and in the far-infrared (FIR) at 60micron. These samples are built using the GALEX and IRAS sky surveys over ~600 square degrees. The NUV selected sample contains 95 galaxies detected down to NUV=16mag (AB system). 83 galaxies in this sample are spirals or irregulars and only two of them are not detected at 60micron. The FIR selected sample is built from the IRAS PSCz catalog complete down to 0.6Jy. Among the 163 sources, we select 118 star forming galaxies well measured by IRAS, all but 1 are detected in NUV and 14 galaxies are not detected in the far-ultraviolet (FUV) band (1350-1750 A, lambda_m = 1530 A). The dust to ultraviolet (NUV and FUV) flux ratio is calibrated to estimate the dust attenuation at both wavelengths. The median value of the attenuation in NUV is found to be ~1 mag for the NUV selected sample versus ~2 mag for the FIR selected one. Within both samples, the dust attenuation is found to correlate with the luminosity of the galaxies. Almost all the NUV selected galaxies and 2/3 of the FIR selected sample exhibit a lower dust attenuation than that expected from the tight relation found previously for starburst galaxies between the dust attenuation and the slope of the ultraviolet continuum. The situation is inverse for one third of the FIR selected galaxies: their extinction is higher than that deduced from their FUV-NUV color and the relation valid for starbursts.