Star formation in galaxies is inefficient, and understanding how star formation is regulated in galaxies is one of the most fundamental challenges of contemporary astrophysics. Radiative cooling, feedback from supernovae and active galactic nuclei, large-scale dynamics and dissipation of turbulent energy act over various time and spatial scales, and all regulate star formation in a complex gas cycle. This paper presents the physics implemented in a new semi-analytical model of galaxy formation and evolution: G.A.S.. The fundamental underpinning of our new model is the development of a multi-phase interstellar medium in which energy produced by supernovae and active galactic nuclei maintains an equilibrium between the diffuse, hot, stable gas and a cooler, clumpy, low-volume filling factor gas. The hot gas is susceptible to thermal and dynamical instabilities. We include a description of how turbulence leads to the formation of giant molecular clouds through an inertial turbulent energy cascade. We explicitly model the evolution of the velocity dispersion at different scales of the cascade and account for thermal instabilities in the hot halo gas. We show that rapid and multiple exchanges between diffuse and unstable gas phases strongly regulates star-formation rates in galaxies because only a small fraction of the unstable gas is forming stars. For high mass halos, cooling is naturally regulated by the growth of thermal instabilities, so we do not need to implement strong AGN feedback. The characteristic timescales describing the gas cycle are in good agreement with observations. Results are in good agreement with the observed stellar mass function from z=6.0 to z=0.5. Thermal instabilities and the cascade of turbulent energy in the dense gas phase introduce a delay between gas accretion and star formation, which keeps galaxy growth inefficient in the early Universe.

Dust is a crucial component of the interstellar medium of galaxies. The presence of dust strongly affects the light produced by stars within a galaxy. As these photons are our main information vector to explore the stellar mass assembly and therefore understand a galaxy's evolution, modeling the luminous properties of galaxies and taking into account the impact of the dust is a fundamental challenge for semi-analytical models.We present the complete prescription of dust attenuation implemented in the new semi-analytical model: G.A.S. This model is based on a two-phase medium originating from a physically motivated turbulent model of gas structuring (G.A.S. I paper). Dust impact is treated by taking into account three dust components: Polycyclic Aromatic Hydrocarbons, Very Small Grains, and Big Grains. All three components evolve in both a diffuse and a fragmented/dense gas phase. Each phase has its own stars, dust content and geometry. Dust content evolves according to the metallicity of it associated phase.The G.A.S. model is used to predict both the UV and the IR luminosity functions from $z=9.0$ to $z=0.1$. Our two-phase ISM prescription catches very well the evolution of UV and IR luminosity functions. We note a small overproduction of the IR luminosity at low redshift ($z<0.5$). We also focus on the Infrared-Excess (IRX) and explore its dependency with the stellar mass, UV slope, stellar age, metallicity and slope of the attenuation curves. Our model predicts large scatters for relations based on IRX, especially for the IRX-$\beta$ relation. Our analysis reveals that the slope of the attenuation curve is more driven by absolute attenuation in the FUV band than by disk inclination. We confirm that the age of the stellar population and the slope of the attenuation curve can both shift galaxies below the fiducial star-birth relation in the IRX-$\beta$ diagram.

The [CII] fine structure transition at 158 microns is the dominant cooling line of cool interstellar gas, and is the brightest of emission lines from star forming galaxies from FIR through meter wavelengths. With the advent of ALMA and NOEMA, capable of detecting [CII]-line emission in high-redshift galaxies, there has been a growing interest in using the [CII] line as a probe of the physical conditions of the gas in galaxies, and as a SFR indicator at z>4. In this paper, we use a semi-analytical model of galaxy evolution (G.A.S.) combined with the code CLOUDY to predict the [CII] luminosity of a large number of galaxies at 4< z<8. At such high redshift, the CMB represents a strong background and we discuss its effects on the luminosity of the [CII] line. We study the LCII-SFR and LCII-Zg relations and show that they do not strongly evolve with redshift from z=4 and to z=8. Galaxies with higher [CII] luminosities tend to have higher metallicities and higher star formation rates but the correlations are very broad, with a scatter of about 0.5 dex for LCII-SFR. Our model reproduces the LCII-SFR relations observed in high-redshift star-forming galaxies, with [CII] luminosities lower than expected from local LCII-SFR relations. Accordingly, the local observed LCII-SFR relation does not apply at high-z. Our model naturally produces the [CII] deficit, which appears to be strongly correlated with the intensity of the radiation field in our simulated galaxies. We then predict the [CII] luminosity function, and show that it has a power law form in the range of LCII probed by the model with a slope alpha=1. The slope is not evolving from z=4 to z=8 but the number density of [CII]-emitters decreases by a factor of 20x. We discuss our predictions in the context of current observational estimates on both the differential and cumulative luminosity functions.

Follow-up observations at high-angular resolution of submillimeter galaxies showed that the single-dish sources are comprised of a blend of several galaxies. Consequently, number counts derived from low and high angular resolution observations are in disagreement. This demonstrates the importance of resolution effects and the need to have realistic simulations to explore them. We built a new 2deg^2 simulation of the extragalactic sky from the far-infrared to the submillimeter. It is based on an updated version of the two star-formation mode galaxy evolution model. Using global galaxy properties, we use the abundance matching technique to populate a dark-matter lightcone and thus simulate the clustering. We produce maps and extract the sources, and show that the limited angular resolution of single-dish instruments have a strong impact on (sub)millimeter continuum observations. Taking into account these resolution effects, we are reproducing a large set of observables, including number counts, redshift distributions, and cosmic infrared background power spectra. Our simulation describes consistently the number counts from single-dish telescopes and interferometers. In particular, at 350 and 500 um, we find that number counts measured by Herschel between 5 and 50 mJy are biased towards high values by a factor 2, and that redshift distributions are biased towards low z. We also show that the clustering has an important impact on the Herschel pixel histogram used to derive number counts from P(D) analysis. Finally, we demonstrate that the large number density of red Herschel sources found in observations but not in models could be an observational artifact caused by the combination of noise, resolution effects, and the steepness of color and flux density distributions. Our simulation (SIDES) is available at http://cesam.lam.fr/sides

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.

We present an extended version of the semi-analytical model, GalICS. Like its predecessor, eGalICS applies a post-treatment of the baryonic physics on pre-computed dark-matter merger trees extracted from an N-body simulation. We review all the mechanisms that affect, at any given time, the formation and evolution of a galaxy in its host dark-matter halo. We mainly focus on the gas cycle from the smooth cosmological accretion to feedback processes. To follow this cycle with high accuracy we introduce some novel prescriptions: i) a smooth baryonic accretion with two phases: a cold mode and a hot mode built on the continuous dark-matter accretion. In parallel to this smooth accretion, we implement the standard photoionisation modelling to reduce the input gas flow on the smallest structures. ii) a complete monitoring of the hot gas phase. We compute the evolution of the core density, the mean temperature and the instantaneous escape fraction of the hot atmosphere by considering that the hot gas is in hydrostatic equilibrium in the dark-matter potential well, and by applying a principle of conservation of energy on the treatment of gas accretion, supernovae and super massive black hole feedback iii) a new treatment for disc instabilities based on the formation, the migration and the disruption of giant clumps. The migration of such clumps in gas-rich galaxies allows to form pseudo-bulges. The different processes in the gas cycle act on different time scales, and we thus build an adaptive time-step scheme to solve the evolution equations. The model presented here is compared in detail to the observations of stellar-mass functions, star formation rates, and luminosity functions, in a companion paper.

Semi-analytical models (SAMs) are currently the best way to understand the formation of galaxies within the cosmic dark-matter structures. While they fairly well reproduce the local stellar mass functions, correlation functions and luminosity functions, they fail to match observations at high redshift (z > 3) in most cases, particularly in the low-mass range. The inconsistency between models and observations indicates that the history of gas accretion in galaxies, within their host dark-matter halo, and the transformation of gas into stars, are not well followed. Hereafter, we briefly present a new version of the GalICS semi-analytical model. We explore the impacts of classical mechanisms, such as supernova feedback or photoionization, on the evolution of the stellar mass assembly. Even with a strong efficiency, these two processes cannot explain the observed stellar mass function and star formation rate distribution and some other relations. We thus introduce an ad-hoc modification of the standard paradigm, based on the presence of a \textit{no-star-forming} gas component, and a concentration of the star-forming gas in galaxy discs. The main idea behind the existence of the no-star-forming gas reservoir is that only a fraction of the total gas mass in a galaxy is available to form stars. The reservoir generates a delay between the accretion of the gas and the star formation process. This new model is in much better agreement with the observations of the stellar mass function in the low-mass range than the previous models, and agrees quite well with a large set of observations, including the redshift evolution of the specific star formation rate. However, it predicts a large fraction of no-star-forming baryonic gas, potentially larger than observed, even if its nature has still to be examined in the context of the missing baryon problem.

[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.

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