Source author record

M. Béthermin

M. Béthermin 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

43works
3topics
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

43 published item(s)

preprint2024arXiv

CONCERTO: Extracting the power spectrum of the [C II ] emission line

CONCERTO is the first experiment to perform a [CII] line intensity mapping survey to target $z>5.2$. Measuring the [CII] power spectrum allows us to study the role of dusty star-forming galaxies in the star formation history during the Reionization and post-Reionization. The main obstacle to this measurement is the contamination by bright foregrounds. We evaluate our ability to retrieve the [CII] signal in mock observations using the Simulated Infrared Dusty Extragalactic Sky. We compared two methods for dealing with the dust continuum emission from galaxies: the standard PCA and the arPLS method. For line interlopers, the strategy relies on masking low-redshift galaxies using external catalogues. As we do not have observations of CO or classical CO proxies ,we relied on the COSMOS stellar mass catalogue. To measure the power spectrum of masked data, we adapted the P of K EstimatoR and discuss its use on LIM data. The arPLS method achieves a reduction of the continuum background to a sub-dominant level of the [CII] at z=7 by a factor of>70. When using PCA, this factor is only 0.7. The masking lowers the power amplitude of line contamination down to $2 \times 10^2 Jy^2/sr$ This residual level is dominated by faint undetected sources. For our [CII] model, this results in a detection at z = 5.2 with a power ratio [CII]/(residual interlopers) = $62 \pm 32$ for a 22 % area survey loss. However, at z = 7, [C II ] / (residual interlopers)$=2.0 \pm 1.4$. Thanks to the large area covered by SIDES-Uchuu, we show that the power amplitude of line residuals varies by 12-15% for z=5.2-7. We present an end-to-end simulation of the extragalactic foreground removal that we ran to detect the [CII] at high redshift via its power spectrum. We show that dust continuum emission are not a limiting foreground for [CII] LIM. Residual CO and [CI] limits our ability to measure the [CII] power spectrum at z>7.

preprint2022arXiv

ALPINE: A Large Survey to Understand Teenage Galaxies

A multiwavelength study of galaxies is important to understand their formation and evolution. Only in the recent past, thanks to the Atacama Large (Sub) Millimeter Array (ALMA), were we able to study the far-infrared (IR) properties of galaxies at high redshifts. In this article, we summarize recent research highlights and their significance to our understanding of early galaxy evolution from the ALPINE survey, a large program with ALMA to observe the dust continuum and 158um C+ emission of normal star-forming galaxies at z = 4-6. Combined with ancillary data at UV through near-IR wavelengths, ALPINE provides the currently largest multiwavelength sample of post-reionization galaxies and has advanced our understanding of (i) the demographics of C+ emission; (ii) the relation of star formation and C+ emission; (iii) the gas content; (iv) outflows and enrichment of the intergalactic medium; and (v) the kinematics, emergence of disks, and merger rates in galaxies at z > 4. ALPINE builds the basis for more detailed measurements with the next generation of telescopes, and places itself as an important post-reionization baseline sample to allow a continuous study of galaxies over 13 billion years of cosmic time.

preprint2022arXiv

GOODS-ALMA 2.0: Starbursts in the main sequence reveal compact star formation regulating galaxy evolution prequenching

Compact star formation appears to be generally common in dusty star-forming galaxies (SFGs). However, its role in the framework set by the scaling relations in galaxy evolution remains to be understood. In this work we follow up on the galaxy sample from the GOODS-ALMA 2.0 survey, an ALMA blind survey at 1.1mm covering a continuous area of 72.42arcmin$^2$ using two array configurations. We derived physical properties, such as star formation rates, gas fractions, depletion timescales, and dust temperatures for the galaxy sample built from the survey. There exists a subset of galaxies that exhibit starburst-like short depletion timescales, but they are located within the scatter of the so-called main sequence of SFGs. These are dubbed starbursts in the main sequence and display the most compact star formation and they are characterized by the shortest depletion timescales, lowest gas fractions, and highest dust temperatures of the galaxy sample, compared to typical SFGs at the same stellar mass and redshift. They are also very massive, accounting for $\sim 60\%$ of the most massive galaxies in the sample ($\log (M_{\rm{*}}/M_{\odot}) > 11.0$). We find trends between the areas of the ongoing star formation regions and the derived physical properties for the sample, unveiling the role of compact star formation as a physical driver of these properties. Starbursts in the main sequence appear to be the extreme cases of these trends. We discuss possible scenarios of galaxy evolution to explain the results drawn from our galaxy sample. Our findings suggest that the star formation rate is sustained in SFGs by gas and star formation compression, keeping them within the main sequence even when their gas fractions are low and they are presumably on the way to quiescence.

preprint2022arXiv

High resolution spectral imaging of CO(7-6), [CI](2-1) and continuum of three high-z lensed dusty star-forming galaxies using ALMA

High-redshift dusty star-forming galaxies with very high star formation rates (500 -- 3000 M$_{\odot}$ yr$^{-1}$) are key to understanding the formation of the most extreme galaxies in the early Universe. Characterising the gas reservoir of these systems can reveal the driving factor behind the high star formation. Using molecular gas tracers like high-J CO lines, neutral carbon lines and the dust continuum, we can estimate the gas density and radiation field intensity in their interstellar media. In this paper, we present high resolution ($\sim$0.4$^{\prime\prime}$) observations of CO(7-6), [CI](2-1) and dust continuum of 3 lensed galaxies from the SPT-SMG sample at $z\sim$ 3 with the Atacama Large Millimeter/submillimeter Array. Our sources have high intrinsic star-formation rates ($>$850 M$_{\odot}$yr$^{-1}$) and rather short depletion timescales ($<$100 Myr). Based on the L$_{[\rm CI](2-1)}$/L$_{\rm CO(7-6)}$ and L$_{[\rm CI](2-1)}$/L$_{\rm IR}$ ratios, our galaxy sample has similar radiation field intensities and gas densities compared to other submillimetre galaxies. We perform visibility-based lens modelling on these objects to reconstruct the kinematics in the source plane. We find that the cold gas masses of the sources are compatible with simple dynamical mass estimates using ULIRG-like values of the CO-H$_2$ conversion factor $α_{\rm CO}$ but not Milky Way-like values. We find diverse source kinematics in our sample: SPT0103-45 and SPT2147-50 are likely rotating disks while SPT2357-51 is possibly a major merger. The analysis presented in the paper could be extended to a larger sample to determine better statistics of morphologies and interstellar medium properties of high-$z$ dusty star-forming galaxies.

preprint2022arXiv

The ALPINE-ALMA [CII] Survey: Investigation of 10 Galaxies at $z\sim4.5$ with [OII] and [CII] Line Emission $-$ ISM Properties and [OII]$-$SFR Relation

We present $10$ main-sequence ALPINE galaxies (log($M/M_{\odot}$) = 9.2-11.1 and ${\rm SFR}=23-190\,{\rm M_{\odot}\,yr^{-1}}$) at $z\sim4.5$ with optical [OII] measurements from Keck/MOSFIRE spectroscopy and Subaru/MOIRCS narrow-band imaging. This is the largest such multi-wavelength sample at these redshifts, combining various measurements in the ultra-violet, optical, and far-infrared including [CII]$_{158{\rm μm}}$ line emission and dust continuum from ALMA and H$α$ emission from Spitzer photometry. For the first time, this unique sample allows us to analyze the relation between [OII] and total star-formation rate (SFR) and the interstellar medium (ISM) properties via [OII]/[CII] and [OII]/\halpha luminosity ratios at $z\sim4.5$. The [OII]$-$SFR relation at $z\sim4.5$ cannot be described using standard local descriptions, but is consistent with a metal-dependent relation assuming metallicities around $50\%$ solar. To explain the measured dust-corrected luminosity ratios of $L[OII]/L[CII] \sim 0.98^{+0.21}_{-0.22}$ and $L[OII]/LHa \sim -0.22^{+0.13}_{-0.15}$ for our sample, ionization parameters $\log(U)< -2$ and electron densities $\log(\rm n_e / {\rm [cm^{-3}]}) \sim 2.5-3$ are required. The former is consistent with galaxies at $z\sim2-3$, however lower than at $z>6$. The latter may be slightly higher than expected given the galaxies' specific SFR. The analysis of this pilot sample suggests that typical log($ M/M_{\odot})$ > 9 galaxies at $z\sim4.5$ to have broadly similar ISM properties as their descendants at $z\sim2$ and suggest a strong evolution of ISM properties since the Epoch of Reionization at $z>6$.

preprint2022arXiv

The ALPINE-ALMA [CII] Survey: the population of [CII]-undetected galaxies and their role in the $\mathrm{L_{[CII]}}$-SFR relation

The [CII] 158$~μ$m emission line represents so far one of the most profitable tools for the investigation of the high-redshift galaxies in the early Universe. Being one of the brightest cooling lines in the rest-frame far-infrared regime of star-forming galaxies, it has been successfully exploited as a tracer of star-formation rate (SFR) in local sources. The picture is more complex at higher redshifts, where its usability in this context is still under investigation. Recent results from the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) survey suggest that there is no (or weak) evolution of the L$\mathrm{_{[CII]}}$-SFR relation up to $z\sim6$ but their reliability is hampered by the presence of a large population of [CII] non-detected galaxies. In this work, we characterize the population of [CII] non-detections in ALPINE. By stacking their ALMA spectra, we obtain a signal detected at $\sim5.1σ$, resulting in a [CII] luminosity of $\mathrm{log(L_\mathrm{[CII]}}/\mathrm{L_{\odot}})$ $\sim7.8$. When combining this value with those from the [CII] detections, we find a $\mathrm{L_{[CII]}}$-SFR relation with a slope $b=1.14\pm0.11$, in agreement within the uncertainties both with the linear relation found in the local Universe, and with the previous findings from ALPINE at $z\sim5$. This suggests that the [CII] line can be considered a good tracer of star formation up to the distant Universe. Finally, we show that the galaxies of our sample that most deviate from the observed L$_\mathrm{[CII]}$-SFR relation could suffer from a less precise redshift estimation, perhaps artificially reducing their [CII] luminosity. In this respect, we claim that there is no evidence in favour of a deficit of [CII] content in high-z galaxies, in contrast with earlier studies.

preprint2021arXiv

CONCERTO at APEX: Installation and first phase of on-sky commissioning

CONCERTO (CarbON CII line in post-rEionisation and ReionisaTiOn) is a large field-of-view (FoV) spectro-imager that has been installed on the Cassegrain Cabin of Atacama Pathfinder EXperiment (APEX) telescope in April 2021. CONCERTO hosts 2 focal planes and a total number of 4000 Kinetic Inductance Detectors (KID), with an instantaneous FoV of 18.6 arcminutes in the range of 130-310 GHz. The spectral resolution can be easily tuned down to 1 GHz depending on the scientific target. The scientific program of CONCERTO has many objectives, with two main programs focused on mapping the fluctuations of the [CII] line intensity in the reionisation and post-reionisation epoch (4.5<z<8.5), and on studying galaxy clusters via the thermal and kinetic Sunyaev-Zel'dovich (SZ) effect. CONCERTO will also measure the dust and molecular gas contents of local and intermediate-redshift galaxies, it will study the Galactic star-forming clouds and finally it will observe the CO intensity fluctuations arising from 0.3<z<2 galaxies. The design of the instrument, installation at APEX and current status of the commissioning phase and science verification will be presented. Also we describe the deployment and first on-sky tests performed between April and June 2021.

preprint2021arXiv

GOODS-ALMA 2.0: Source catalog, number counts, and prevailing compact sizes in 1.1 mm galaxies

Submillimeter/millimeter observations of dusty star-forming galaxies with the Atacama Large Millimeter/submillimeter Array (ALMA) have shown that dust continuum emission generally occurs in compact regions smaller than the stellar distribution. However, it remains to be understood how systematic these findings are. Studies often lack homogeneity in the sample selection, target discontinuous areas with inhomogeneous sensitivities, and suffer from modest $uv$ coverage coming from single array configurations. GOODS-ALMA is a 1.1mm galaxy survey over a continuous area of 72.42arcmin$^2$ at a homogeneous sensitivity. In this version 2.0, we present a new low resolution dataset and its combination with the previous high resolution dataset from the survey, improving the $uv$ coverage and sensitivity reaching an average of $σ= 68.4μ$Jy beam$^{-1}$. A total of 88 galaxies are detected in a blind search (compared to 35 in the high resolution dataset alone), 50% at $S/N_{peak} \geq 5$ and 50% at $3.5 \leq S/N_{peak} \leq 5$ aided by priors. Among them, 13 out of the 88 are optically dark or faint sources ($H$- or $K$-band dropouts). The sample dust continuum sizes at 1.1mm are generally compact, with a median effective radius of $R_{e} = 0"10 \pm 0"05$ (a physical size of $R_{e} = 0.73 \pm 0.29$kpc at the redshift of each source). Dust continuum sizes evolve with redshift and stellar mass resembling the trends of the stellar sizes measured at optical wavelengths, albeit a lower normalization compared to those of late-type galaxies. We conclude that for sources with flux densities $S_{1.1mm} > 1$mJy, compact dust continuum emission at 1.1mm prevails, and sizes as extended as typical star-forming stellar disks are rare. The $S_{1.1mm} < 1$mJy sources appear slightly more extended at 1.1mm, although they are still generally compact below the sizes of typical star-forming stellar disks.

preprint2021arXiv

Optical and near-infrared observations of the SPT2349-56 proto-cluster core at z = 4.3

We present Gemini-S and {\it Spitzer}-IRAC optical-through-near-IR observations in the field of the SPT2349-56 proto-cluster at $z=4.3$. We detect optical/IR counterparts for only nine of the 14 submillimetre galaxies (SMGs) previously identified by ALMA in the core of SPT2349-56. In addition, we detect four $z\sim4$ Lyman-break galaxies (LBGs) in the 30 arcsec diameter region surrounding this proto-cluster core. Three of the four LBGs are new systems, while one appears to be a counterpart of one of the nine observed SMGs. We identify a candidate brightest cluster galaxy (BCG) with a stellar mass of $(3.2^{+2.5}_{-1.4})\times10^{11}\,{\rm M}_{\odot}$. The stellar masses of the eight other SMGs place them on, above, and below the main sequence of star formation at $z\approx4.5$. The cumulative stellar mass for the SPT2349-56 core is at least $(11.5\pm2.9)\times10^{11}\,{\rm M}_{\odot}$, a sizeable fraction of the stellar mass in local BCGs, and close to the universal baryon fraction (0.16) relative to the virial mass of the core ($10^{13}\,{\rm M}_{\odot}$). As all 14 of these SMGs are destined to quickly merge, we conclude that the proto-cluster core has already developed a significant stellar mass at this early stage, comparable to $z=1$ BCGs. Importantly, we also find that the SPT2349-56 core structure would be difficult to uncover in optical surveys, with none of the ALMA sources being easily identifiable or constrained through $g,r,$ and $i$ colour-selection in deep optical surveys and only a modest overdensity of LBGs over the extended core structure. SPT2349-56 therefore represents a truly dust-obscured phase of a massive cluster core under formation.

preprint2021arXiv

Searching for high-z DSFGs with NIKA2 and NOEMA

As the possible progenitors of passive galaxies at z=2-3, dusty star-forming galaxies (DSFGs) at z>4 provide a unique perspective to study the formation, assembly, and early quenching of massive galaxies in the early Universe. The extreme obscuration in optical-IR makes (sub)mm spectral scans the most universal and unbiased way to confirm/exclude the high-z nature of candidate dusty star-forming galaxies. We present here the status of the NIKA2 Cosmological Legacy Survey (N2CLS), which is the deepest wide-area single-dish survey in the millimeter searching for high-z DSFGs. We also introduce a joint-analysis method to efficiently search for the spectroscopic redshift of high-z DSFGs with noisy spectra and photometric data and present its success in identifying the redshift of DSFGs found in NIKA2 science verification data.

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

The [CII]/[NII] ratio in 3 < z < 6 sub-millimetre galaxies from the South Pole Telescope survey

We present Atacama Compact Array and Atacama Pathfinder Experiment observations of the [N II] 205 $μ$m fine-structure line in 40 sub-millimetre galaxies lying at redshifts z = 3 to 6, drawn from the 2500 deg$^2$ South Pole Telescope survey. This represents the largest uniformly selected sample of high-redshift [N II] 205 $μ$m measurements to date. 29 sources also have [C II] 158 $μ$m line observations allowing a characterization of the distribution of the [C II] to [N II] luminosity ratio for the first time at high-redshift. The sample exhibits a median L$_{[C II]}$ /L$_{[N II]}$ $\approx$ 11 and interquartile range of 5.0 to 24.7. These ratios are similar to those observed in local (U)LIRGs, possibly indicating similarities in their interstellar medium. At the extremes, we find individual sub-millimetre galaxies with L$_{[C II]}$ /L$_{[N II]}$ low enough to suggest a smaller contribution from neutral gas than ionized gas to the [C II] flux and high enough to suggest strongly photon or X-ray region dominated flux. These results highlight a large range in this line luminosity ratio for sub-millimetre galaxies, which may be caused by variations in gas density, the relative abundances of carbon and nitrogen, ionization parameter, metallicity, and a variation in the fractional abundance of ionized and neutral interstellar medium.

preprint2019arXiv

Probing the cold magnetized Universe with SPICA-POL (B-BOP)

SPICA, the cryogenic infrared space telescope recently pre-selected for a `Phase A' concept study as one of the three remaining candidates for ESA's fifth medium class (M5) mission, is foreseen to include a far-infrared polarimetric imager (SPICA-POL, now called B-BOP), which would offer a unique opportunity to resolve major issues in our understanding of the nearby, cold magnetized Universe. This paper presents an overview of the main science drivers for B-BOP, including high dynamic range polarimetric imaging of the cold interstellar medium (ISM) in both our Milky Way and nearby galaxies. Thanks to a cooled telescope, B-BOP will deliver wide-field 100-350 micron images of linearly polarized dust emission in Stokes Q and U with a resolution, signal-to-noise ratio, and both intensity and spatial dynamic ranges comparable to those achieved by Herschel images of the cold ISM in total intensity (Stokes I). The B-BOP 200 micron images will also have a factor ~30 higher resolution than Planck polarization data. This will make B-BOP a unique tool for characterizing the statistical properties of the magnetized interstellar medium and probing the role of magnetic fields in the formation and evolution of the interstellar web of dusty molecular filaments giving birth to most stars in our Galaxy. B-BOP will also be a powerful instrument for studying the magnetism of nearby galaxies and testing galactic dynamo models, constraining the physics of dust grain alignment, informing the problem of the interaction of cosmic rays with molecular clouds, tracing magnetic fields in the inner layers of protoplanetary disks, and monitoring accretion bursts in embedded protostars.

preprint2016arXiv

A giant Ly$α$ nebula in the core of an X-ray cluster at $z=1.99$: implications for early energy injection

We present the discovery of a giant $\gtrsim$100~kpc Ly$α$ nebula detected in the core of the X-ray emitting cluster CL~J1449+0856 at $z=1.99$ through Keck/LRIS narrow-band imaging. This detection extends the known relation between Ly$α$ nebulae and overdense regions of the Universe to the dense core of a $5-7\times10^{13}$ M$_{\odot}$ cluster. The most plausible candidates to power the nebula are two Chandra-detected AGN host cluster members, while cooling from the X-ray phase and cosmological cold flows are disfavored primarily because of the high Ly$α$ to X-ray luminosity ratio ($L_{\mathrm{Lyα}}/L_{\mathrm{X}} \approx0.3$, $\gtrsim10-1000\times$ higher than in local cool-core clusters) and by current modeling. Given the physical conditions of the Ly$α$-emitting gas and the possible interplay with the X-ray phase, we argue that the Ly$α$ nebula would be short-lived ($\lesssim10$ Myr) if not continuously replenished with cold gas at a rate of $\gtrsim1000$ M$_{\odot}$ yr$^{-1}$. We investigate the possibility that cluster galaxies supply the required gas through outflows and we show that their total mass outflow rate matches the replenishment necessary to sustain the nebula. This scenario directly implies the extraction of energy from galaxies and its deposition in the surrounding intracluster medium, as required to explain the thermodynamic properties of local clusters. We estimate an energy injection of the order of $\thickapprox2$ keV per particle in the intracluster medium over a $2$ Gyr interval. In our baseline calculation AGN provide up to $85$% of the injected energy and 2/3 of the mass, while the rest is supplied by supernovae-driven winds.

preprint2016arXiv

An ALMA view of the interstellar medium of the z=4.77 lensed starburst SPT-S J213242-5802.9

We present ALMA detections of the [NII] 205$μ$m and CO(12-11) emission lines, and the tentative detection of [CI] $^3$P$_1$ - $^3$P$_0$ for the strongly lensed ($μ$=5.7$\pm$0.5) dusty, star-forming galaxy SPT-S J213242-5802.9 (hereafter SPT2132-58) at z=4.77. The [NII] and CO(12-11) lines are detected at 11.5 and 8.5 $σ$ level, respectively, by our band-6 observations. The [CI] line is detected at 3.2 $σ$ after a re-analysis of existing band-3 data. The [CI] luminosity implies a gas mass of 3.8$\pm$1.2$\times$10$^{10}$ M$_\odot$, and consequently a very short depletion timescale of 34$\pm$13 Myr and a CO-luminosity-to-gas-mass conversion factor $α_{\rm CO}$ of 1.0$\pm$0.3 M$_\odot$ (K km s$^{-1}$ pc$^{2}$)$^{-1}$. SPT2132-58 is an extreme starburst with an intrinsic star formation rate of 1100$\pm$200 M$_\odot$/yr. We find a [CII]/[NII] ratio of 26$\pm$7, which is the highest reported at z$>$4. This suggests that SPT2132-58 hosts an evolved interstellar medium (0.5 Z$_\odot$<Z<1.5 Z$_\odot$), which may be dominated by photodissociation regions. The CO(2-1) and CO(5-4) transitions have lower CO-to-far-infrared ratios than local and high-redshift samples, while CO(12-11) is similar to these samples, suggesting the presence of an additional very excited component or an AGN.

preprint2016arXiv

The redshift distribution of dusty star forming galaxies from the SPT survey

We use the Atacama Large Millimeter/submillimeter Array (ALMA) in Cycle 1 to determine spectroscopic redshifts of high-redshift dusty star-forming galaxies (DSFGs) selected by their 1.4mm continuum emission in the South Pole Telescope (SPT) survey. We present ALMA 3mm spectral scans between 84-114GHz for 15 galaxies and targeted ALMA 1mm observations for an additional eight sources. Our observations yield 30 new line detections from CO, [CI] , [NII] , H_2O and NH_3. We further present APEX [CII] and CO mid-J observations for seven sources for which only a single line was detected in spectral-scan data from ALMA Cycle 0 or Cycle 1. We combine the new observations with previously published and new mm/submm line and photometric data of the SPT-selected DSFGs to study their redshift distribution. The combined data yield 39 spectroscopic redshifts from molecular lines, a success rate of >85%. Our sample represents the largest data set of its kind today and has the highest spectroscopic completeness among all redshift surveys of high-z DSFGs. The median of the redshift distribution is z=3.9+/-0.4, and the highest-redshift source in our sample is at z=5.8. We discuss how the selection of our sources affects the redshift distribution, focusing on source brightness, selection wavelength, and strong gravitational lensing. We correct for the effect of gravitational lensing and find the redshift distribution for 1.4mm-selected sources with a median redshift of z=3.1+/-0.3. Comparing to redshift distributions selected at shorter wavelengths from the literature, we show that selection wavelength affects the shape of the redshift distribution.

preprint2015arXiv

CO excitation of normal star forming galaxies out to z=1.5 as regulated by the properties of their interstellar medium

We investigate the CO excitation of normal star forming disk galaxies at z=1.5 using IRAM PdBI observations of the CO[2-1], CO[3-2] and CO[5-4] transitions for 4 galaxies, including VLA observations of CO[1-0] for 3 of them, with the aim of constraining the average state of H2 gas. Exploiting prior knowledge of the velocity range, spatial extent and size of the CO emission we measure reliable line fluxes with S/N>4-7 for individual transitions. While the average CO Spectral Line Energy Distribution (SLED) has a sub-thermal excitation similar to the Milky Way (MW) up to CO[3-2], we show that the average CO[5-4] emission is four times stronger than assuming MW excitation. This demonstrates the presence of an additional component of more excited, denser and possibly warmer molecular gas. The ratio of CO[5-4] to lower-J CO emission is lower than in local (U)LIRGs and high-redshift SMGs, correlating closely with the average intensity of the radiation field <U> and with the star formation surface density, but not with the SF efficiency (SFE). The CO[5-4] luminosity correlates linearly with LIR over 4 orders of magnitudes, with z=1.5 BzK galaxies following the same trend as local spirals and (U)LIRGs and high redshift star bursting SMGs. The CO[5-4] luminosity is thus empirically related to the dense gas, and might be a more convenient way to probe it than standard high--density tracers that are much fainter than CO. We see excitation variations among our sample galaxies, that can be linked to their evolutionary state and clumpiness in optical rest frame images. In one galaxy we see spatially resolved excitation variations, where the more highly excited part of the galaxy corresponds to the location of massive SF clumps. This supports to models that suggest that giant clumps are the main source of the high excitation CO emission in high redshift disk-like galaxies.

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

GOODS-$Herschel$: identification of the individual galaxies responsible for the 80-290$μ$m cosmic infrared background

We propose a new method of pushing $Herschel$ to its faintest detection limits using universal trends in the redshift evolution of the far infrared over 24$μ$m colours in the well-sampled GOODS-North field. An extension to other fields with less multi-wavelength information is presented. This method is applied here to raise the contribution of individually detected $Herschel$ sources to the cosmic infrared background (CIRB) by a factor 5 close to its peak at 250$μ$m and more than 3 in the 350$μ$m and 500$μ$m bands. We produce realistic mock $Herschel$ images of the deep PACS and SPIRE images of the GOODS-North field from the GOODS-$Herschel$ Key Program and use them to quantify the confusion noise at the position of individual sources, i.e., estimate a "local confusion noise". Two methods are used to identify sources with reliable photometric accuracy extracted using 24$μ$m prior positions. The clean index (CI), previously defined but validated here with simulations, which measures the presence of bright 24$μ$m neighbours and the photometric accuracy index (PAI) directly extracted from the mock $Herschel$ images. After correction for completeness, thanks to our mock $Herschel$ images, individually detected sources make up as much as 54% and 60% of the CIRB in the PACS bands down to 1.1 mJy at 100$μ$m and 2.2 mJy at 160$μ$m and 55, 33, and 13% of the CIRB in the SPIRE bands down to 2.5, 5, and 9 mJy at 250$μ$m, 350$μ$m, and 500$μ$m, respectively. The latter depths improve the detection limits of $Herschel$ by factors of 5 at 250$μ$m, and 3 at 350$μ$m and 500$μ$m as compared to the standard confusion limit. Interestingly, the dominant contributors to the CIRB in all $Herschel$ bands appear to be distant siblings of the Milky Way ($z$$\sim$0.96 for $λ$$<$300$μ$m) with a stellar mass of $M_{\star}$$\sim$9$\times$10$^{10}$M$_{\odot}$.

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

Probing star formation in the dense environments of z~1 lensing halos aligned with dusty star-forming galaxies detected with the South Pole Telescope

We probe star formation in the environments of massive $\sim10^{13}\,M_{\odot}$ dark matter halos at redshifts of $z$$\sim$$1$. This star formation is linked to a sub-millimetre clustering signal which we detect in maps of the Planck High Frequency Instrument that are stacked at the positions of a sample of high-redshift ($z$$>$$2$) strongly-lensed dusty star-forming galaxies (DSFGs) selected from the South Pole Telescope (SPT) 2500 deg$^2$ survey. The clustering signal has sub-millimetre colours which are consistent with the mean redshift of the foreground lensing halos ($z$$\sim$$1$). We report a mean excess of star formation rate (SFR) compared to the field, of $(2700\pm700)\,M_{\odot}\,{yr}^{-1}$ from all galaxies contributing to this clustering signal within a radius of 3.5' from the SPT DSFGs. The magnitude of the Planck excess is in broad agreement with predictions of a current model of the cosmic infrared background. The model predicts that 80$\%$ of the excess emission measured by Planck originates from galaxies lying in the neighbouring halos of the lensing halo. Using Herschel maps of the same fields, we find a clear excess, relative to the field, of individual sources which contribute to the Planck excess. The mean excess SFR compared to the field is measured to be ($370\pm40)$$\,M_{\odot}\,{yr}^{-1}$ per resolved, clustered source. Our findings suggest that the environments around these massive $z$$\sim$$1$ lensing halos host intense star formation out to about $2\,$Mpc. The flux enhancement due to clustering should also be considered when measuring flux densities of galaxies in Planck data.

preprint2015arXiv

Satellite content and quenching of star formation in galaxy groups at z~1.8

We study the properties of satellites in the environment of massive star-forming galaxies at z~1.8 in the COSMOS field, using a sample of 215 galaxies on the main sequence of star formation with an average mass of 10^11 Msun. At z>1.5, these galaxies typically trace halos of mass >10^13 Msun. We use optical-near-infrared photometry to estimate stellar masses and star formation rates (SFR) of centrals and satellites down to ~6*10^9 Msun. We stack data around 215 central galaxies to statistically detect their satellite halos, finding an average of ~3 galaxies in excess of the background density. We fit the radial profiles of satellites with simple beta-models, and compare their integrated properties to model predictions. We find that the total stellar mass of satellites amounts to 68% of the central galaxy, while SED modeling and far-infrared photometry consistently show their total SFR to be 25-35% of the central's rate. We also see significant variation in the specific SFR of satellites within the halo with, in particular, a sharp decrease at <100 kpc. After considering different potential explanations, we conclude that this is likely an environmental signature of the hot inner halo. This effect can be explained in the first order by a simple free-fall scenario, suggesting that these low-mass environments can shut down star formation in satellites on relatively short timescales of ~0.3 Gyr.

preprint2015arXiv

Stellar masses and star formation rates of lensed dusty star-forming galaxies from the SPT survey

To understand cosmic mass assembly in the Universe at early epochs, we primarily rely on measurements of stellar mass and star formation rate of distant galaxies. In this paper, we present stellar masses and star formation rates of six high-redshift ($2.8\leq z \leq 5.7$) dusty, star-forming galaxies (DSFGs) that are strongly gravitationally lensed by foreground galaxies. These sources were first discovered by the South Pole Telescope (SPT) at millimeter wavelengths and all have spectroscopic redshifts and robust lens models derived from ALMA observations. We have conducted follow-up observations, obtaining multi-wavelength imaging data, using {\it HST}, {\it Spitzer}, {\it Herschel} and the Atacama Pathfinder EXperiment (APEX). We use the high-resolution {\it HST}/WFC3 images to disentangle the background source from the foreground lens in {\it Spitzer}/IRAC data. The detections and upper limits provide important constraints on the spectral energy distributions (SEDs) for these DSFGs, yielding stellar masses, IR luminosities, and star formation rates (SFRs). The SED fits of six SPT sources show that the intrinsic stellar masses span a range more than one order of magnitude with a median value $\sim$ 5 $\times 10^{10}M_{\Sun}$. The intrinsic IR luminosities range from 4$\times 10^{12}L_{\Sun}$ to 4$\times 10^{13}L_{\Sun}$. They all have prodigious intrinsic star formation rates of 510 to 4800 $M_{\Sun} {\rm yr}^{-1}$. Compared to the star-forming main sequence (MS), these six DSFGs have specific SFRs that all lie above the MS, including two galaxies that are a factor of 10 higher than the MS. Our results suggest that we are witnessing the ongoing strong starburst events which may be driven by major mergers.

preprint2015arXiv

The imprint of rapid star formation quenching on the spectral energy distributions of galaxies

[Abridged] In high density environment, the gas content of galaxies is stripped, leading to a rapid quenching of their star formation activity. This dramatic environmental effect is generally not taken into account in the SFHs usually assumed to perform spectral energy distribution (SED) fitting of these galaxies, yielding to a poor fit of their stellar emission and, consequently, a biased estimate of the SFR. We aim at reproducing the SFH of galaxies that underwent a rapid star formation quenching using a truncated delayed SFH that we implemented in the SED fitting code CIGALE. We show that the ratio between the instantaneous SFR and the SFR just before the quenching ($r_{SFR}$) is well constrained as long as rest frame UV data are available. This SED modelling is applied to the Herschel Reference Survey (HRS) containing isolated galaxies and sources falling in the dense environment of the Virgo cluster. The latter are HI-deficient due to ram pressure stripping. We show that the truncated delayed SFH successfully reproduces their SED while typical SFH assumptions fail. A good correlation is found between $r_{SFR}$ and HI-def, the parameter quantifying the gas deficiency of cluster galaxies, meaning that SED fitting results can be used to provide a tentative estimate of the gas deficiency of galaxies for which HI observations are not available. The HRS galaxies are placed on the SFR-$M_*$ diagram showing that the HI-deficient sources lie in the quiescent region confirming previous studies. Using the $r_{SFR}$ parameter, we derive the SFR of these sources before quenching and show that they were previously on the main sequence relation. We show that the $r_{SFR}$ parameter is also well recovered for deeply obscured high redshift sources, as well as in absence of IR data. SED fitting is thus a powerful tool to identify galaxies that underwent a rapid star formation quenching.

preprint2014arXiv

Modeling the connection between ultraviolet and infrared galaxy populations across cosmic times

Using a phenomenological approach, we self-consistently model the redshift evolution of the ultraviolet (UV) and infrared (IR) luminosity functions across cosmic time, as well as a range of observed IR properties of UV-selected galaxy population. This model is an extension of the 2SFM (2 star-formation modes) formalism, which is based on the observed "main-sequence" of star-forming galaxies, i.e. a strong correlation between their stellar mass and their star formation rate (SFR), and a secondary population of starbursts with an excess of star formation. The balance between the UV light from young, massive stars and the dust-reprocessed IR emission is modeled following the empirical relation between the attenuation (IRX for IR excess hereafter) and the stellar mass, assuming a scatter of 0.4\,dex around this relation. We obtain a good overall agreement with the measurements of the IR luminosity function up to z~3 and the UV luminosity functions up to z~6, and show that a scatter on the IRX-M relation is mandatory to reproduce these observables. We also naturally reproduce the observed, flat relation between the mean IRX and the UV luminosity at L$_{\rm UV}>$10$^{9.5}$ L$_\odot$. Finally, we perform predictions of the UV properties and detectability of IR-selected samples and the vice versa, and discuss the results in the context of the UV-rest-frame and sub-millimeter surveys of the next decade.

preprint2014arXiv

Probing the galaxy-halo connection in UltraVISTA to $z\sim2$

We use percent-level precision photometric redshifts in the UltraVISTA-DR1 near-infrared survey to investigate the changing relationship be- tween galaxy stellar mass and the dark matter haloes hosting them to $z\sim2$. We achieve this by measuring the clustering properties and abundances of a series of volume-limited galaxy samples selected by stellar mass and star-formation activity. We interpret these results in the framework of a phenomenological halo model and numerical simulations. Our measurements span a uniquely large range in stellar mass and redshift and reach below the characteristic stellar mass to $z\sim2$. Our results are: 1. At fixed redshift and scale, clustering amplitude depends monotonically on sample stellar mass threshold; 2. At fixed angular scale, the projected clustering amplitude decreases with redshift but the co-moving correlation length remains constant; 3. Characteristic halo masses and galaxy bias increase with increasing median stellar mass of the sample; 4. The slope of these relationships is modified in lower mass haloes; 5. Concerning the passive galaxy population, characteristic halo masses are consistent with a simply less-abundant version of the full galaxy sample, but at lower redshifts the fraction of satellite galaxies in the passive population is very different from the full galaxy sample; 6. Finally we find that the ratio between the characteristic halo mass and median stellar mass at each redshift bin reaches a peak at ${\log(M_h/M_\odot)}\sim12.2$ and the position of this peak remains constant out to $z\sim2$. The behaviour of the full and passively evolving galaxy samples can be understood qualitatively by considering the slow evolution of the characteristic stellar mass in the redshift range probed by our survey.

preprint2014arXiv

Regularity underlying complexity: a redshift-independent description of the continuous variation of galaxy-scale molecular gas properties in the mass-star formation rate plane

Star-forming galaxies (SFGs) display a continuous distribution of specific star formation rates (sSFR) which can be approximated by the superposition of two log-normal distributions. The 1st of these encompasses the main sequence (MS) of SFGs, the 2nd one a rarer population of starbursts (SB). We show that the sSFR-distribution of SBs can be regarded as the result of a physical process (plausibly merging) taking the mathematical form of a log-normal boosting kernel and enhancing star formation activity. We explore the utility of splitting the star-forming population into MS and SB galaxies - an approach we term "2-Star Formation Mode" (2-SFM) framework - for understanding their molecular gas properties. Variations of star formation efficiency (SFE) and gas fractions among SFGs take a simple, redshift-independent form, once these quantities are normalized to the value of an average MS galaxy. The change in SFE for galaxies undergoing a starburst event scales supra-linearly with the SFR-increase, as expected for mergers. This implies a continuous distribution of galaxies in the Schmidt-Kennicutt plane that separates more clearly into loci for SBs and normal galaxies than observed in SFR vs. M* space. SBs with the largest deviations (>10-fold) from the MS, like many local ULIRGs, are not average SBs, but even rarer events with progenitors having larger gas fractions than typical MS galaxies. We statistically infer that the gas fractions of typical SBs decrease by a factor of 2 to 3 with respect to their direct MS progenitors, as expected to occur in short-term SFR-boosts during which internal gas reservoirs are drained more quickly than gas is accreted from the cosmic web. We predict variations of the CO-to-H2 conversion factor in the SFR-M* plane and provide evidence that the higher sSFR of distant galaxies is a direct consequence of larger gas fraction in these systems.

preprint2013arXiv

Extragalactic millimeter-wave point source catalog, number counts and statistics from 771 square degrees of the SPT-SZ Survey

We present a point source catalog from 771 square degrees of the South Pole Telescope Sunyaev Zel'dovich (SPT-SZ) survey at 95, 150, and 220 GHz. We detect 1545 sources above 4.5 sigma significance in at least one band. Based on their relative brightness between survey bands, we classify the sources into two populations, one dominated by synchrotron emission from active galactic nuclei, and one dominated by thermal emission from dust-enshrouded star-forming galaxies. We find 1238 synchrotron and 307 dusty sources. We cross-match all sources against external catalogs and find 189 unidentified synchrotron sources and 189 unidentified dusty sources. The dusty sources without counterparts are good candidates for high-redshift, strongly lensed submillimeter galaxies. We derive number counts for each population from 1 Jy down to roughly 9, 5, and 11 mJy at 95, 150, and 220 GHz. We compare these counts with galaxy population models and find that none of the models we consider for either population provide a good fit to the measured counts in all three bands. The disparities imply that these measurements will be an important input to the next generation of millimeter-wave extragalactic source population models.

preprint2013arXiv

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

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

preprint2013arXiv

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

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

preprint2013arXiv

The deepest Herschel-PACS far-infrared survey: number counts and infrared luminosity functions from combined PEP/GOODS-H observations

We present results from the deepest Herschel-PACS (Photodetector Array Camera and Spectrometer) far-infrared blank field extragalactic survey, obtained by combining observations of the GOODS (Great Observatories Origins Deep Survey) fields from the PACS Evolutionary Probe (PEP) and GOODS-Herschel key programmes. We describe data reduction and the construction of images and catalogues. In the deepest parts of the GOODS-S field, the catalogues reach 3-sigma depths of 0.9, 0.6 and 1.3 mJy at 70, 100 and 160 um, respectively, and resolve ~75% of the cosmic infrared background at 100um and 160um into individually detected sources. We use these data to estimate the PACS confusion noise, to derive the PACS number counts down to unprecedented depths and to determine the infrared luminosity function of galaxies down to LIR=10^11 Lsun at z~1 and LIR=10^12 Lsun at z~2, respectively. For the infrared luminosity function of galaxies, our deep Herschel far-infrared observations are fundamental because they provide more accurate infrared luminosity estimates than those previously obtained from mid-infrared observations. Maps and source catalogues (>3-sigma) are now publicly released. Combined with the large wealth of multi-wavelength data available for the GOODS fields, these data provide a powerful new tool for studying galaxy evolution over a broad range of redshifts.

preprint2013arXiv

The evolution of the dust temperatures of galaxies in the SFR$-M_{\ast}$ plane up to $z$$\,\thicksim\,$$2$

[Abridged] We study the evolution of the dust temperatures of galaxies in the SFR-M* plane up to z~2 using observations from the Herschel Space Observatory. Starting from a sample of galaxies with reliable star-formation rates (SFRs), stellar masses (M*) and redshift estimates, we grid the SFR-M* parameter space in several redshift ranges and estimate the mean Tdust of each SFR-M*-z bin. Dust temperatures are inferred using the stacked far-infrared flux densities of our SFR-M*-z bins. At all redshifts, Tdust increases with infrared luminosities (LIR), specific SFRs (SSFR; i.e., SFR/M*) and distances with respect to the main sequence (MS) of the SFR-M* plane (i.e., D_SSFR_MS=log[SSFR(galaxy)/SSFR_MS(M*,z)]). The Tdust-SSFR and Tdust-D_SSFR_MS correlations are statistically more significant than the Tdust-LIR one. While the slopes of these three correlations are redshift-independent, their normalizations evolve from z=0 and z~2. We convert these results into a recipe to derive Tdust from SFR, M* and z. The existence of a strong Tdust-D_SSFR_MS correlation provides us with information on the dust and gas content of galaxies. (i) The slope of the Tdust-D__SSFR_MS correlation can be explained by the increase of the star-formation efficiency (SFE; SFR/Mgas) with D_SSFR_MS as found locally by molecular gas studies. (ii) At fixed D_SSFR_MS, the constant Tdust observed in galaxies probing large ranges in SFR and M* can be explained by an increase or decrease of the number of star-forming regions with comparable SFE enclosed in them. (iii) At high redshift, the normalization towards hotter temperature of the Tdust-D_SSFR_MS correlation can be explained by the decrease of the metallicities of galaxies or by the increase of the SFE of MS galaxies. All these results support the hypothesis that the conditions prevailing in the star-forming regions of MS and far-above-MS galaxies are different.

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

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

The contribution of starbursts and normal galaxies to infrared luminosity functions at z < 2

We present a parameter-less approach to predict the shape of the infrared (IR) luminosity function (LF) at redshifts z < 2. It requires no tuning and relies on only three observables: (1) the redshift evolution of the stellar mass function for star-forming galaxies, (2) the evolution of the specific star formation rate (sSFR) of main-sequence galaxies, and (3) the double-Gaussian decomposition of the sSFR-distribution at fixed stellar mass into a contribution (assumed redshift- and mass-invariant) from main-sequence and starburst activity. This self-consistent and simple framework provides a powerful tool for predicting cosmological observables: observed IR LFs are successfully matched at all z < 2, suggesting a constant or only weakly redshift-dependent contribution (8-14%) of starbursts to the star formation rate density. We separate the contributions of main-sequence and starburst activity to the global IR LF at all redshifts. The luminosity threshold above which the starburst component dominates the IR LF rises from log(LIR/Lsun) = 11.4 to 12.8 over 0 < z < 2, reflecting our assumed (1+z)^2.8-evolution of sSFR in main-sequence galaxies.

preprint2012arXiv

The Herschel Multi-tiered Extragalactic Survey: HerMES

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

preprint2012arXiv

The hidden "AGN main sequence": Evidence for a universal black hole accretion to star formation rate ratio since z~2 producing a M_BH-M* relation

Using X-ray stacking analyses we estimate the average amounts of supermassive black hole (SMBH) growth taking place in star-forming galaxies (SFGs) at z~1 and z~2 as a function of galaxy stellar mass (M*). We find the average SMBH growth rate follows remarkably similar trends with M* and redshift as the average star-formation rates (SFRs) of their host galaxies (i.e., (i.e., dM_BH/dt ~ M*^(0.86+/-0.39) for the z~1 sample and dM_BH/dt ~ M*^(1.05+/-0.36) for the z~2 sample). It follows that the ratio of SMBH growth rate to SFR is (a) flat with respect to M* (b) not evolving with redshift and (c) close to the ratio required to maintain/establish a SMBH to M* ratio of ~10^{-3} as also inferred from today's M_BH-M_bulge relationship. We interpret this as evidence that SMBHs have, on average, grown in-step with their host galaxies since at least z~2, irrespective of host galaxy mass and AGN triggering mechanism. As such, we suggest that the same secular processes that drive the bulk of star formation are also responsible for the majority of SMBH growth. From this, we speculate that it is the availability of gas reservoirs that regulate both cosmological SMBH growth and star formation.

preprint2010arXiv

Submillimeter Number Counts at 250, 350 and 500 microns in BLAST data

BLAST (Balloon-borne Large-Aperture Submillimeter Telescope) performed the first deep and wide extragalactic survey at 250, 350 and 500 um. The extragalactic number counts at these wavelengths are important constraints for modeling the infrared galaxies evolution. [...] We use three methods to identify the submillimeter sources. 1) Blind extraction. [...] The photometry is computed with a new simple and quick PSF fitting routine (FASTPHOT). [...] 2) Extraction with prior. [...] 3) A stacking analysis. [...] With the blind extraction, we reach 97, 83 and 76 mJy at resp. 250, 350 and 500 um with a 95% completeness. With the prior extraction, we reach 76 mJy (resp. 63 and 49 mJy) at 250 um (resp. 350 and 500 um). With the stacking analysis, we reach 6.2 mJy (resp. 5.2 and 3.5 mJy) at 250 um (resp. 350 and 500 um). The differential submillimeter number counts are derived, and start showing a turnover at flux densities decreasing with increasing wavelength. There is a very good agreement with the P(D) analysis of Patanchon et al. (2009). At bright fluxes (>100 mJy), the Lagache et al. (2004) and Le Borgne et al. (2009) models slightly overestimate the observed counts, but there is a very good agreement near the peak of differential number counts. [...] Counts are available at: http://www.ias.u-psud.fr/irgalaxies/downloads.php

preprint2010arXiv

The Cosmic Far-Infrared Background Buildup Since Redshift 2 at 70 and 160 microns in the COSMOS and GOODS fields

The Cosmic Far-Infrared Background (CIB) at wavelengths around 160 μm corresponds to the peak intensity of the whole Extragalactic Background Light, which is being measured with increasing accuracy. However, the build up of the CIB emission as a function of redshift, is still not well known. Our goal is to measure the CIB history at 70 μm and 160 μm at different redshifts, and provide constraints for infrared galaxy evolution models. We use complete deep Spitzer 24 μm catalogs down to about 80 μJy, with spectroscopic and photometric redshifts identifications, from the GOODS and COSMOS deep infrared surveys covering 2 square degrees total. After cleaning the Spitzer/MIPS 70 μm and 160 μm maps from detected sources, we stacked the far-IR images at the positions of the 24 μm sources in different redshift bins. We measured the contribution of each stacked source to the total 70 and 160 μm light, and compare with model predictions and recent far-IR measurements made with Herschel/PACS on smaller fields. We have detected components of the 70 and 160 μm backgrounds in different redshift bins up to z ~ 2. The contribution to the CIB is maximum at 0.3 <= z <= 0.9 at 160μm (and z <= 0.5 at 70 μm). A total of 81% (74%) of the 70 (160) μm background was emitted at z < 1. We estimate that the AGN relative contribution to the far-IR CIB is less than about 10% at z < 1.5. We provide a comprehensive view of the CIB buildup at 24, 70, 100, 160 μm. IR galaxy models predicting a major contribution to the CIB at z < 1 are in agreement with our measurements, while our results discard other models that predict a peak of the background at higher redshifts. Our results are available online http://www.ias.u-psud.fr/irgalaxies/ .

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.