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 star formation picture provided by a total of 10 497 individual far-infrared detections is supplemented by the stacking analysis of a mass complete sample of 62 361 starforming galaxies from the Hubble Space Telescope (HST) H band-selected catalogs of the CANDELS survey and from two deep ground-based K s band-selected catalogs in the GOODS-North and the COSMOS-wide field 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", which may be easily generalized to other experiments. The combination of direct UV and far-infrared UV-reprocessed light provides a complete census on the star formation rates (SFRs), allowing us to demonstrate that galaxies at z = 4 to 0 of all stellar masses (M * ) follow a universal scaling law, the so-called main sequence of star-forming galaxies. We find a universal close-to-linear slope of the log 10 (SFR)-log 10 (M * ) relation, with evidence for a flattening of the main sequence at high masses (log 10 (M * /M ) > 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, which mostly happens over the same redshift range. Within the main sequence, we measure a nonvarying SFR dispersion of 0.3 dex: at a fixed redshift and stellar mass, about 68% of star-forming galaxies form stars at a universal rate within a factor 2. 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 on the origin of present-day stars: more than two-thirds 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 billion years.

We present the deepest 100 to 500 μm far-infrared observations obtained with the Herschel Space Observatory as part of the GOODS-Herschel key program, and examine the infrared (IR) 3-500 μm spectral energy distributions (SEDs) of galaxies at 0 < z < 2.5, supplemented by a local reference sample from IRAS, ISO, Spitzer, and AKARI data. We determine the projected star formation densities of local galaxies from their radio and mid-IR continuum sizes. We find that the ratio of total IR luminosity to rest-frame 8 μm luminosity, IR8 (≡L tot IR /L 8 ), follows a Gaussian distribution centered on IR8 = 4 (σ = 1.6) and defines an IR main sequence for star-forming galaxies independent of redshift and luminosity. Outliers from this main sequence produce a tail skewed toward higher values of IR8. This minority population (<20%) is shown to consist of starbursts with compact projected star formation densities. IR8 can be used to separate galaxies with normal and extended modes of star formation from compact starbursts with high-IR8, high projected IR surface brightness (Σ IR > 3 × 10 10 L kpc −2 ) and a high specific star formation rate (i.e., starbursts). The rest-frame, UV-2700 Å size of these distant starbursts is typically half that of main sequence galaxies, supporting the correlation between star formation density and starburst activity that is measured for the local sample. Locally, luminous and ultraluminous IR galaxies, (U)LIRGs (L tot IR ≥ 10 11 L ), are systematically in the starburst mode, whereas most distant (U)LIRGs form stars in the "normal" main sequence mode. This confusion between two modes of star formation is the cause of the so-called "mid-IR excess" population of galaxies found at z > 1.5 by previous studies. Main sequence galaxies have strong polycyclic aromatic hydrocarbon (PAH) emission line features, a broad far-IR bump resulting from a combination of dust temperatures (T dust ∼ 15-50 K), and an effective T dust ∼ 31 K, as derived from the peak wavelength of their infrared SED. Galaxies in the starburst regime instead exhibit weak PAH equivalent widths and a sharper far-IR bump with an effective T dust ∼ 40 K. Finally, we present evidence that the mid-to-far IR emission of X-ray active galactic nuclei (AGN) is predominantly produced by star formation and that candidate dusty AGNs with a power-law emission in the mid-IR systematically occur in compact, dusty starbursts. After correcting for the effect of starbursts on IR8, we identify new candidates for extremely obscured AGNs.

Aims. Photometric data of galaxies covering the rest-frame wavelength range from far-UV to far-IR make it possible to derive galaxy properties with a high reliability by fitting the attenuated stellar emission and the related dust emission at the same time.Methods. For this purpose we wrote the code CIGALE (Code Investigating GALaxy Emission) that uses model spectra composed of the Maraston (or PEGASE) stellar population models, synthetic attenuation functions based on a modified Calzetti law, spectral line templates, the Dale & Helou dust emission models, and optional spectral templates of obscured AGN. Depending on the input redshifts, filter fluxes were computed for the model set and compared to the galaxy photometry by carrying out a Bayesian-like analysis. CIGALE was tested by analysing 39 nearby galaxies selected from SINGS. The reliability of the different model parameters was evaluated by studying the resulting expectation values and their standard deviations in relation to the input model grid. Moreover, the influence of the filter set and the quality of photometric data on the code results was estimated. Results. For up to 17 filters with effective wavelengths between 0.15 and 160 μm, we find robust results for the mass, star formation rate, effective age of the stellar population at 4000 Å, bolometric luminosity, luminosity absorbed by dust, and attenuation in the far-UV. Details of the star formation history (excepting the burst fraction) and the shape of the attenuation curve are difficult to investigate with the available broad-band UV and optical photometric data. A study of the mutual relations between the reliable properties confirms the dependence of star formation activity on morphology in the local Universe and indicates a significant drop in this activity at about 10 11 M towards higher total stellar masses. The dustiest galaxies in the SINGS sample are present in the same mass range. On the other hand, the far-UV attenuation of our sample galaxies does not appear to show a significant dependence on star formation activity. Conclusions. The results for our SINGS test sample demonstrate that CIGALE can be a valuable tool for studying basic properties of galaxies in the near and distant Universe if UV-to-IR data are available.

Context. Measuring how the physical properties of galaxies change across cosmic times is essential to understand galaxy formation and evolution. With the advent of numerous ground-based and space-borne instruments launched over the past few decades we now have exquisite multi-wavelength observations of galaxies from the far-ultraviolet (FUV) to the radio domain. To tap into this mine of data and obtain new insight into the formation and evolution of galaxies, it is essential that we are able to extract information from their spectral energy distribution (SED). Aims. We present a completely new implementation of Code Investigating GALaxy Emission (CIGALE). Written in python, its main aims are to easily and efficiently model the FUV to radio spectrum of galaxies and estimate their physical properties such as star formation rate, attenuation, dust luminosity, stellar mass, and many other physical quantities. Methods. To compute the spectral models, CIGALE builds composite stellar populations from simple stellar populations combined with highly flexible star formation histories, calculates the emission from gas ionised by massive stars, and attenuates both the stars and the ionised gas with a highly flexible attenuation curve. Based on an energy balance principle, the absorbed energy is then reemitted by the dust in the mid-and far-infrared domains while thermal and non-thermal components are also included, extending the spectrum far into the radio range. A large grid of models is then fitted to the data and the physical properties are estimated through the analysis of the likelihood distribution. Results. CIGALE is a versatile and easy-to-use tool that makes full use of the architecture of multi-core computers, building grids of millions of models and analysing samples of thousands of galaxies, both at high speed. Beyond fitting the SEDs of galaxies and parameter estimations, it can also be used as a model-generation tool or serve as a library to build new applications.

The Herschel Multi‐tiered Extragalactic Survey (HerMES) is a legacy programme designed to map a set of nested fields totalling ∼380 deg2. Fields range in size from 0.01 to ∼20 deg2, using the Herschel‐Spectral and Photometric Imaging Receiver (SPIRE) (at 250, 350 and 500 μm) and the Herschel‐Photodetector Array Camera and Spectrometer (PACS) (at 100 and 160 μm), with an additional wider component of 270 deg2 with SPIRE alone. These bands cover the peak of the redshifted thermal spectral energy distribution from interstellar dust and thus capture the reprocessed optical and ultraviolet radiation from star formation that has been absorbed by dust, and are critical for forming a complete multiwavelength understanding of galaxy formation and evolution. The survey will detect of the order of 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.

We aim to measure the average dust and molecular gas content of massive star-forming galaxies (>3 × 10 10 M ) 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 is 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 = 4 in main-sequence galaxies. We can reproduce this evolution with simple models that account for the decrease in 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 M mol /(M mol + M )) with redshift and found a similar, steeply increasing trend for both samples. At z ∼ 4, this fraction reaches ∼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.

Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA.

We compare the dust attenuation properties of two samples of galaxies purely selected in the near-ultraviolet (NUV) band (1750-2750Å, λ mean = 2310Å) and in the far-infrared (FIR) at 60 µm . These samples are built using the GALEX