Source author record

Matthieu Béthermin

Matthieu 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

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

21 published item(s)

preprint2022arXiv

CONCERTO: a breakthrough in wide field-of-view spectroscopy at millimeter wavelengths

CarbON CII line in post-rEionization and ReionizaTiOn (CONCERTO) is a low-resolution spectrometer with an instantaneous field-of-view of 18.6 arcmin, operating in the 130-310 GHz transparent atmospheric window. It is installed on the 12-meter Atacama Pathfinder Experiment (APEX) telescope at 5100 m above sea level. The Fourier transform spectrometer (FTS) contains two focal planes hosting a total of 4304 kinetic inductance detectors. The FTS interferometric pattern is recorded on the fly while continuously scanning the sky. One of the goals of CONCERTO is to characterize the large-scale structure of the Universe by observing the integrated emission from unresolved galaxies. This methodology is an innovative technique and is called line intensity mapping. In this paper, we describe the CONCERTO instrument, the effect of the vibration of the FTS beamsplitter, and the status of the CONCERTO main survey.

preprint2022arXiv

The ALPINE-ALMA [CII] survey: The infrared-radio correlation and AGN fraction of star-forming galaxies at z $\sim$ 4.4-5.9

We present the radio properties of 66 spectroscopically-confirmed normal star-forming galaxies (SFGs) at $4.4<z<5.9$ in the COSMOS field that were [C II] detected in the Atacama Large Millimeter Array (ALMA) Large Program to INvestigate [C II] at Early times (ALPINE). We separate these galaxies ("CII-detected-all") into lower redshift ("CII-detected-lz", $\langle z\rangle=4.5$) and higher redshift ("CII-detected-hz", $\langle z\rangle=5.6$) sub-samples and stack multi-wavelength imaging for each sub-sample from X-ray to radio bands. A radio signal is detected in the stacked 3 GHz image of CII-detected-all and -lz samples at $\gtrsim3σ$. We find that the infrared-radio correlation of our sample, quantified by $q_{\mathrm{TIR}}$, is lower than the local relation for normal SFGs at $\sim$3$σ$ significance level, and is instead broadly consistent with that of bright sub-mm galaxies at $2<z<5$. Neither of these samples show evidence of dominant AGN activity in their stacked Spectral Energy Distributions (SEDs), rest-frame UV spectra, or X-ray images. Although we cannot rule out the possible effect of the assumed spectral index and the applied infrared SED templates as at least partially causing these differences, the lower obscured fraction of star formation than at lower redshift can alleviate the tension between our stacked $q_{\mathrm{TIR}}$s and that of local normal SFGs. It is possible that the dust buildup, which primarily governs the IR emission in addition to older stellar populations, has not had enough time to occur fully in these galaxies, whereas the radio emission can respond on a more rapid timescale. Therefore, we might expect a lower $q_{\mathrm{TIR}}$ to be a general property of high-redshift SFGs.

preprint2021arXiv

Characterization of Two 2mm-detected Optically-Obscured Dusty Star-Forming Galaxies

The 2mm Mapping Obscuration to Reionization with ALMA (MORA) Survey was designed to detect high redshift ($z\gtrsim4$), massive, dusty star-forming galaxies (DSFGs). Here we present two, likely high redshift sources, identified in the survey whose physical characteristics are consistent with a class of optical/near-infrared (OIR) invisible DSFGs found elsewhere in the literature. We first perform a rigorous analysis of all available photometric data to fit spectral energy distributions and estimate redshifts before deriving physical properties based on our findings. Our results suggest the two galaxies, called MORA-5 and MORA-9, represent two extremes of the "OIR-dark" class of DSFGs. MORA-5 ($z_{\rm phot}=4.3^{+1.5}_{-1.3}$) is a significantly more active starburst with a star-formation rate of 830$^{+340}_{-190}$M$_\odot$yr$^{-1}$ compared to MORA-9 ($z_{\rm phot}=4.3^{+1.3}_{-1.0}$) whose star-formation rate is a modest 200$^{+250}_{-60}$M$_\odot$yr$^{-1}$. Based on the stellar masses (M$_{\star}\approx10^{10-11}$M$_\odot$), space density ($n\sim(5\pm2)\times10^{-6}$Mpc$^{-3}$, which incorporates two other spectroscopically confirmed OIR-dark DSFGs in the MORA sample at $z=4.6$ and $z=5.9$), and gas depletion timescales ($<1$Gyr) of these sources, we find evidence supporting the theory that OIR-dark DSFGs are the progenitors of recently discovered $3<z<4$ massive quiescent galaxies.

preprint2021arXiv

Mapping Obscuration to Reionization with ALMA (MORA): 2mm Efficiently Selects the Highest-Redshift Obscured Galaxies

We present the characteristics of 2mm-selected sources from the largest Atacama Large Millimeter and submillimeter Array (ALMA) blank-field contiguous survey conducted to-date, the Mapping Obscuration to Reionization with ALMA (MORA) survey covering 184arcmin$^2$ at 2mm. Twelve of the thirteen detections above 5$σ$ are attributed to emission from galaxies, eleven of which are dominated by cold dust emission. These sources have a median redshift of $\langle z_{\rm 2mm}\rangle=3.6^{+0.4}_{-0.3}$ primarily based on optical/near-infrared (OIR) photometric redshifts with some spectroscopic redshifts, with 77$\pm$11% of sources at $z>3$ and 38$\pm$12% of sources at $z>4$. This implies that 2mm selection is an efficient method for identifying the highest redshift dusty star-forming galaxies (DSFGs). Lower redshift DSFGs ($z<3$) are far more numerous than those at $z>3$ yet likely to drop out at 2mm. MORA shows that DSFGs with star-formation rates in excess of 300M$_\odot$ yr$^{-1}$ and relative rarity of $\sim$10$^{-5}$ Mpc$^{-3}$ contribute $\sim$30% to the integrated star-formation rate density between $3<z<6$. The volume density of 2mm-selected DSFGs is consistent with predictions from some cosmological simulations and is similar to the volume density of their hypothesized descendants: massive, quiescent galaxies at $z>2$. Analysis of MORA sources' spectral energy distributions hint at steeper empirically-measured dust emissivity indices than typical literature studies, with $\langleβ\rangle=2.2^{+0.5}_{-0.4}$. The MORA survey represents an important step in taking census of obscured star-formation in the Universe's first few billion years, but larger area 2mm surveys are needed to more fully characterize this rare population and push to the detection of the Universe's first dusty galaxies.

preprint2020arXiv

Deceptively cold dust in the massive starburst galaxy GN20 at $z\sim4$

We present new observations, carried out with IRAM NOEMA, of the atomic neutral carbon transitions [CI](1-0) at 492 GHz and [CI](2-1) at 809 GHz of GN20, a well-studied star-bursting galaxy at $z=4.05$. The high luminosity line ratio [CI](2-1)/[CI](1-0) implies an excitation temperature of $48^{+14}_{-9}$ K, which is significantly higher than the apparent dust temperature of $T_{\rm d}=33\pm2$ K ($β=1.9$) derived under the common assumption of an optically thin far-infrared dust emission, but fully consistent with $T_{\rm d}=52\pm5$ K of a general opacity model where the optical depth ($τ$) reaches unity at a wavelength of $λ_0=170\pm23$ $μ$m. Moreover, the general opacity solution returns a factor of $\sim 2\times$ lower dust mass and, hence, a lower molecular gas mass for a fixed gas-to-dust ratio, than with the optically thin dust model. The derived properties of GN20 thus provide an appealing solution to the puzzling discovery of starbursts appearing colder than main-sequence galaxies above $z>2.5$, in addition to a lower dust-to-stellar mass ratio that approaches the physical value predicted for starburst galaxies.

preprint2020arXiv

The ALPINE-ALMA [C II] survey: the luminosity function of serendipitous [C II] line emitters at $z\sim 5$

We present the first [CII] 158 $μ$m luminosity function (LF) at $z\sim 5$ from a sample of serendipitous lines detected in the ALMA Large Program to INvestigate [CII] at Early times (ALPINE). A search performed over the 118 ALPINE pointings revealed several serendipitous lines. Based on their fidelity, we selected 14 lines for the final catalog. According to the redshift of their counterparts, we identified 8 out of 14 detections as [CII] lines at $z\sim 5$, and two as CO transitions at lower redshifts. The remaining 4 lines have an elusive identification in the available catalogs and we considered them as [CII] candidates. We used the 8 confirmed [CII] and the 4 [CII] candidates to build one of the first [CII] LFs at $z\sim 5$. We found that 11 out of these 12 sources have a redshift very similar to that of the ALPINE target in the same pointing, suggesting the presence of overdensities around the targets. Therefore, we split the sample in two (a "clustered" and "field" sub-sample) according to their redshift separation and built two separate LFs. Our estimates suggest that there could be an evolution of the [CII] LF between $z \sim 5$ and $z \sim 0$. By converting the [CII] luminosity to star formation rate we evaluated the cosmic star formation rate density (SFRD) at $z\sim 5$. The clustered sample results in a SFRD $\sim 10$ times higher than previous measurements from UV-selected galaxies. On the other hand, from the field sample (likely representing the average galaxy population) we derived a SFRD $\sim 1.6$ higher compared to current estimates from UV surveys but compatible within the errors. Because of the large uncertainties, observations of larger samples are necessary to better constrain the SFRD at $z\sim 5$. This study represents one of the first efforts aimed at characterizing the demography of [CII] emitters at $z\sim 5$ using a mm-selection of galaxies.

preprint2016arXiv

Metal enrichment in a semi-analytical model, fundamental scaling relations, and the case of Milky Way galaxies

Gas flows play a fundamental role in galaxy formation and evolution, providing the fuel for the star formation process. These mechanisms leave an imprint in the amount of heavy elements. Thus, the analysis of this metallicity signature provides additional constraint on the galaxy formation scenario. We aim to discriminate between four different galaxy formation models based on two accretion scenarios and two different star formation recipes. We address the impact of a bimodal accretion scenario and a strongly regulated star formation recipe. We present a new extension of the eGalICS model, which allows us to track the metal enrichment process. Our new chemodynamical model is applicable for situations ranging from metal-free primordial accretion to very enriched interstellar gas contents. We use this new tool to predict the metallicity evolution of both the stellar populations and gas phase. We also address the evolution of the gas metallicity with the star formation rate (SFR). We then focus on a sub-sample of Milky Way-like galaxies. We compare both the cosmic stellar mass assembly and the metal enrichment process of such galaxies with observations and detailed chemical evolution models. Our models, based on a strong star formation regulation, allow us to reproduce well the stellar mass to gas-phase metallicity relation observed in the local universe. However, we observe a systematic shift towards high masses. Our $Mstar-Zg-SFR relation is in good agreement with recent measurements: our best model predicts a clear dependence with the SFR. Both SFR and metal enrichment histories of our Milky Way-like galaxies are consistent with observational measurements and detailed chemical evolution models. We finally show that Milky Way progenitors start their evolution below the observed main sequence and progressively reach this observed relation at z = 0.

preprint2015arXiv

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

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

preprint2015arXiv

The influence of wavelength, flux, and lensing selection effects on the redshift distribution of dusty, star-forming galaxies

We interpret the large variety of redshift distributions of galaxies found by far-infrared and (sub-)millimeter deep surveys depending on their depth and wavelength using the Béthermin et al. (2012) phenomenological model of galaxy evolution. This model reproduces without any new parameter tuning the observed redshift distributions from 100 $μ$m to 1.4 mm, and especially the increase of the median redshift with survey wavelength. This median redshift varies also significantly with the depth of the surveys, and deeper surveys do necessarily not probe higher redshifts. Paradoxically, at fixed wavelength and flux limit, the lensed sources are not always at higher redshift. We found that the higher redshift of 1.4 mm-selected south pole telescope (SPT) sources compared to other SMG surveys is not only caused by the lensing selection, but also by the longer wavelength. This SPT sample is expected to be dominated by a population of lensed main-sequence galaxies and a minor contribution ($\sim$10\%) of unlensed extreme starbursts.

preprint2015arXiv

The nature of the [CII] emission in dusty star-forming galaxies from the SPT survey

We present [CII] observations of 20 strongly lensed dusty star forming galaxies at 2.1 < z < 5.7 using APEX and Herschel. The sources were selected on their 1.4 mm flux (S_1.4mm > 20 mJy) from the South Pole Telescope survey, with far-infrared (FIR) luminosities determined from extensive photometric data. The [CII] line is robustly detected in 17 sources, all but one being spectrally resolved. Eleven out of 20 sources observed in [CII] also have low-J CO detections from ATCA. A comparison with mid- and high-J CO lines from ALMA reveals consistent [CII] and CO velocity profiles, suggesting that there is little differential lensing between these species. The [CII], low-J CO and FIR data allow us to constrain the properties of the interstellar medium. We find [CII] to CO(1-0) luminosity ratios in the SPT sample of 5200 +- 1800, with significantly less scatter than in other samples. This line ratio can be best described by a medium of [CII] and CO emitting gas with a higher [CII] than CO excitation temperature, high CO optical depth tau_CO >> 1, and low to moderate [CII] optical depth tau_CII ~< 1. The geometric structure of photodissociation regions allows for such conditions.

preprint2014arXiv

Dust and gas in luminous proto-cluster galaxies at z=4.05: the case for different cosmic dust evolution in normal and starburst galaxies

We measure the dust and gas content of the three sub-millimeter galaxies (SMGs) in the GN20 proto-cluster at $z=4.05$ using new IRAM Plateau de Bure interferometer (PdBI) CO(4-3) and 1.2-3.3 mm continuum observations. All these three SMGs are heavily dust obscured, with UV-based star formation rate (SFR) estimates significantly smaller than the ones derived from the bolometric infrared (IR), consistent with the spatial offsets revealed by HST and CO imaging. Based also on evaluations of the specific SFR, CO-to-H$_2$ conversion factor and gas depletion timescale, we classify all the three galaxies as starbursts (SBs), although with a lower confidence for GN20.2b that might be a later stage merging event. We place our measurements in the context of the evolutionary properties of main sequence (MS) and SB galaxies. ULIRGs have 3-5 times larger $L'_{\rm CO}/M_{\rm dust}$ and $M_{\rm dust}/M_\star$ ratios than $z=0$ MS galaxies, but by $z\sim2$ the difference appears to be blurred, probably due to differential metallicity evolution. SB galaxies appear to slowly evolve in their $L'_{\rm CO}/M_{\rm dust}$ and $M_{\rm dust}/M_\star$ ratios all the way to $z>6$ (consistent with rapid enrichment of SB events), while MS galaxies rapidly increase in $M_{\rm dust}/M_\star$ from $z=0$ to 2 (due to gas fraction increase, compensated by a decrease of metallicities). While no IR/submm continuum detection is available for indisputably normal massive galaxies at $z>2.5$, we show that if metallicity indeed decrease rapidly for these systems at $z>3$ as claimed in the literature, we should expect a strong decrease of their $M_{\rm dust}/M_\star$, consistent with recent PdBI and ALMA upper limits. We conclude that the $M_{\rm dust}/M_\star$ ratio could be a powerful tool for distinguishing starbursts from normal galaxies at $z>4$.

preprint2014arXiv

Evolution of the dust emission of massive galaxies up to z=4 and constraints on their dominant mode of star formation

We aim to measure the average dust and molecular gas content of massive star-forming galaxies ($\rm > 3 \times 10^{10}\,M_\odot$) up to z=4 in the COSMOS field to determine if the intense star formation observed at high redshift is induced by major mergers or caused by large gas reservoirs. Firstly, we measured the evolution of the average spectral energy distributions as a function of redshift using a stacking analysis of Spitzer, Herschel, LABOCA, and AzTEC data for two samples of galaxies: normal star-forming objects and strong starbursts, as defined by their distance to the main sequence. We found that the mean intensity of the radiation field $< U >$ heating the dust (strongly correlated with dust temperature) increases with increasing redshift up to z$\sim$4 in main-sequence galaxies. We can reproduce this evolution with simple models that account for the decrease of the gas metallicity with redshift. No evolution of $< U >$ with redshift is found in strong starbursts. We then deduced the evolution of the molecular gas fraction (defined here as $\rm M_{\rm mol}/(M_{\rm mol}+M_\star)$) with redshift and found a similar, steeply increasing trend for both samples. At z$\sim$4, this fraction reaches $\sim$60%. The average position of the main-sequence galaxies is on the locus of the local, normal star-forming disks in the integrated Schmidt-Kennicutt diagram (star formation rate versus mass of molecular gas), suggesting that the bulk of the star formation up to z=4 is dominated by secular processes.

preprint2014arXiv

PRISM (Polarized Radiation Imaging and Spectroscopy Mission): An Extended White Paper

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to ESA in May 2013 as a large-class mission for investigating within the framework of the ESA Cosmic Vision program a set of important scientific questions that require high resolution, high sensitivity, full-sky observations of the sky emission at wavelengths ranging from millimeter-wave to the far-infrared. PRISM's main objective is to explore the distant universe, probing cosmic history from very early times until now as well as the structures, distribution of matter, and velocity flows throughout our Hubble volume. PRISM will survey the full sky in a large number of frequency bands in both intensity and polarization and will measure the absolute spectrum of sky emission more than three orders of magnitude better than COBE FIRAS. The aim of this Extended White Paper is to provide a more detailed overview of the highlights of the new science that will be made possible by PRISM

preprint2013arXiv

The redshift-evolution of the distribution of star formation among dark matter halos as seen in the infrared

[Abridged] Recent studies revealed a strong correlation between the star formation rate (SFR) and stellar mass of star-forming galaxies, the so-called star-forming main sequence. An empirical modeling approach (2-SFM) which distinguishes between the main sequence and rarer starburst galaxies is capable of reproducing most statistical properties of infrared galaxies. In this paper, we extend this approach by establishing a connection between stellar mass and halo mass with the technique of abundance matching. Based on a few, simple assumptions and a physically motivated formalism, our model successfully predicts the (cross-)power spectra of the cosmic infrared background (CIB), the cross-correlation between CIB and cosmic microwave background (CMB) lensing, and the correlation functions of bright, resolved infrared galaxies measured by Herschel, Planck, ACT and SPT. We use this model to infer the redshift distribution these observables, as well as the level of correlation between CIB-anisotropies at different wavelengths. We also predict that more than 90% of cosmic star formation activity occurs in halos with masses between 10^11.5 and 10^13.5 Msun. Taking into account subsequent mass growth of halos, this implies that the majority of stars were initially (at z>3) formed in the progenitors of clusters, then in groups at 0.5<z<3 and finally in Milky-Way-like halos at z<0.5. At all redshifts, the dominant contribution to the star formation rate density stems from halos of mass ~10^12 Msun, in which the instantaneous star formation efficiency is maximal (~70%). The strong redshift-evolution of SFR in the galaxies that dominate the CIB is thus plausibly driven by increased accretion from the cosmic web onto halos of this characteristic mass scale. Material available at http://irfu.cea.fr/Sap/Phocea/Page/index.php?id=537.

preprint2012arXiv

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

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

preprint2012arXiv

The early early type: discovery of a passive galaxy at z=3

We present the discovery of a massive, quiescent galaxy at z=2.99. We have obtained a HST/WFC3 spectrum of this object and measured its redshift from the detection of a deep 4000A break consistent with an old population and a high metallicity. By stellar population modeling of both its grism spectrum and broad-band photometry, we derive an age of ~0.7 Gyr, implying a formation redshift of z>4, and a mass >10^11 Msun. Although this passive galaxy is the most distant confirmed so far, we find that it is slightly less compact than other z>2 early-types of similar mass, being overall more analogous to those z~1.6 field early-type galaxies. The discovery of this object shows that early-type galaxies are detectable to at least z=3 and suggests that the diversity of structural properties found in z=1.4-2 ellipticals to earlier epochs could have its origin in a variety of formation histories among their progenitors.

preprint2012arXiv

Where stars form and live at high redshift: clues from the infrared

The relation between dark matter halos and the loci of star formation at high redshift is a pressing question in contemporary cosmology. Matching the abundance of halos to the abundance of infrared (IR) galaxies, we explicit the link between dark matter halo mass (Mh), stellar mass (M*) and star-formation rate (SFR) up to a redshift of 2. Our findings are five-fold. First, we find a strong evolution of the relation between M* and SFR as a function of redshift with an increase of sSFR = SFR/M* by a factor ~30 between z=0 and z= 2.3. Second, we observe a decrease of sSFR with stellar mass. These results reproduce observed trends at redshift z>0.3. Third, we find that the star formation is most efficient in dark matter halos with Mh~5x10^11 Msun, with hints of an increase of this mass with redshift. Fourth, we find that SFR/Mh increases by a factor ~15 between z = 0 and z = 2.3. Finally we find that the SFR density is dominated by halo masses close to ~7x10^11 Msun at all redshift, with a rapid decrease at lower and higher halo masses. Despite its simplicity, our novel use of IR observations unveils some characteristic mass-scales governing star formation at high redshift.

preprint2011arXiv

Measurements of the spectral energy distribution of the cosmic infrared background

The extragalactic background light (EBL) is the relic emission of all processes of structure formation in the Universe. About half of this background, called the Cosmic Infrared Background (CIB) is emitted in the 8-1000 microns range, and peaks around 150 microns. It is due to the dust reemission from star formation processes and AGN emission. The CIB spectral energy distribution (SED) constraints the models of star formation in the Universe. It is also useful to compute the opacity of the Universe to the TeV photons. We present the different types of measurements of the CIB and discuss their strengths and weaknesses. 1. The absolute SED was measured by COBE, and by other experiments. These measurements are limited by the accuracy of the component separation, i.e. the foreground subtraction. 2. Robust lower limits are determined from the extragalactic number counts of infrared galaxies. These lower limits are very stringent up to 100 microns. At larger wavelengths, the rather low angular resolution of the instruments limits strongly the depth of the number counts. The "stacking" method determines the flux emitted at a given wavelength by a population detected at another wavelength, and provides stringent lower limits in the sub-mm range. It is complementary with other methods based on the statistical analysis of the map properties like the P(D) analysis. 3. Finally, upper limits can be derived from the high energy spectra of extragalactic sources. These upper limits give currently good constraints in the near- and mid-IR. Progress have been amazing since the CIB discovery about 15 years ago: the SED is much better known, and most of it can be accounted for by galaxies (directly or indirectly). Prospects are also exciting, with fluctuation analysis with Planck&Herschel, and forthcoming missions.

preprint2011arXiv

Modeling the evolution of infrared galaxies : clustering of galaxies in the Cosmic Infrared Background

Star-forming galaxies are a highly biased tracer of the underlying dark matter density field. Their clustering can be studied through the cosmic infrared background anisotropies. These anisotropies have been measured from 100 \mum to 2 mm in the last few years. In this paper, we present a fully parametric model allowing a joint analysis of these recent observations. In order to develop a coherent model at various wavelengths, we rely on two building blocks. The first one is a parametric model that describes the redshift evolution of the luminosity function of star-forming galaxies. It compares favorably to measured differential number counts and luminosity functions. The second one is a halo model based description of the clustering of galaxies. Starting from a fiducial model, we investigate parameter degeneracies using a Fisher analysis. We then discuss how halo of different mass and redshift, how LIRGs and ULIRGs, contribute to the CIB angular power spectra. From the Fisher analysis, we conclude that we cannot constrain the parameters of the model of evolution of galaxies using clustering data only. The use of combined data of C\ell, counts and luminosity functions improves slightly the constraints but does not remove any degeneracies. On the contrary, the measurement of the anisotropies allows us to set interesting constraints on the halo model parameters, even if some strong degeneracies remain. Using our fiducial model, we establish that the 1-halo and 2-halo terms are not sensitive to the same mass regime. We also illustrate how the 1-halo term can be misinterpreted with the Poisson noise term. Conclusions. We present a new model of the clustering of infrared galaxies. However such a model has a few limitations, as the parameters of the halo occupation suffer from strong degeneracies.

preprint2011arXiv

Modeling the evolution of infrared galaxies: A Parametric backwards evolution model

We aim at modeling the infrared galaxy evolution in an as simple as possible way and reproduce statistical properties among which the number counts between 15 microns and 1.1 mm, the luminosity functions, and the redshift distributions. We then aim at using this model to interpret the recent observations (Spitzer, Akari, BLAST, LABOCA, AzTEC, SPT and Herschel), and make predictions for future experiments like CCAT or SPICA. This model uses an evolution in density and luminosity of the luminosity function with two breaks at redshift ~0.9 and 2 and contains the two populations of the Lagache et al. (2004) model: normal and starburst galaxies. We also take into account the effect of the strong lensing of high-redshift sub-millimeter galaxies. It has 13 free parameters and 8 additional calibration parameters. We fit the parameters to the IRAS, Spitzer, Herschel and AzTEC measurements with a Monte-Carlo Markov chain. The model ajusted on deep counts at key wavelengths reproduces the counts from the mid-infrared to the millimeter wavelengths, as well as the mid-infrared luminosity functions. We discuss the contribution to the cosmic infrared background (CIB) and to the infrared luminosity density of the different populations. We also estimate the effect of the lensing on the number counts, and discuss the recent discovery by the South Pole Telescope (SPT) of a very bright population lying at high-redshift. We predict confusion level for future missions using a P(D) formalism, and the Universe opacity to TeV photons due to the CIB.

preprint2010arXiv

Spitzer Deep and Wide Legacy Mid- and Far-Infrared Number Counts and Lower Limits of Cosmic Infrared Background

Aims: We want to place stronger lower limits on the Cosmic Infrared Background (CIB) brightness at 24, 70 and 160 um and measure the extragalactic number counts at these wavelengths in an homogeneous way from various surveys. Methods: Using Spitzer legacy data over 53.6 deg^2 of various depths, we build catalogs with the same extraction method at each wavelength. [...] Number count uncertainties are estimated with a counts-in-cells moment method, in order to take into account galaxy clustering. Furthermore, we use a stacking analysis to estimate number counts of sources not detected at 70 and 160 um. [...] Results: Number counts reach 35 muJy (resp. 3.5 mJy and 40 mJy) at 24 um (resp. 70 and 160um). We reach deeper flux densities of 0.38 mJy at 70 and 3.1 at 160 um using a stacking analysis. [...] Conclusions: These counts are integrated to estimate new CIB firm lower limits of 2.21-0.09+0.09 nW.m^-2.sr^-1 (resp. 5.2-0.4+0.4 and 9.0-1.1+1.1 nW.m^-2.sr^-1) at 24 um (resp. 70 and 160 um), and extrapolated to give new estimates of the CIB due to galaxies of 2.78-0.16+0.18 nW.m^-2.sr^-1 (resp. 6.4-0.6+0.7 and 14.6-2.9+7.1 nW.m^-2.sr^-1). Products (PSF, counts, CIB contributions, software) are publicly available for download at http://www.ias.u-psud.fr/irgalaxies/