No Arabic abstract
We argue that the interplay between cosmic rays, the initial mass function, and star formation plays a crucial role in regulating the star-forming main sequence. To explore these phenomena we develop a toy model for galaxy evolution in which star formation is regulated by a combination of a temperature-dependent initial mass function and heating due to starlight, cosmic rays and, at very high redshift, the cosmic microwave background. This produces an attractor, near-equilibrium solution which is consistent with observations of the star-forming main sequence over a broad redshift range. Additional solutions to the same equations may correspond to other observed phases of galaxy evolution including quiescent galaxies. This model makes several falsifiable predictions, including higher metallicities and dust masses than anticipated at high redshift and isotopic abundances in the Milky Way. It also predicts that stellar mass-to-light ratios are lower than produced using a Milky Way-derived IMF, so that inferences of stellar masses and star formation rates for high redshift galaxies are overestimated. In some cases, this may also transform inferred dark matter profiles from core-like to cusp-like.
The origin of the star forming main sequence ( i.e., the relation between star formation rate and stellar mass, globally or on kpc-scales; hereafter SFMS) remains a hotly debated topic in galaxy evolution. Using the ALMA-MaNGA QUEnching and STar formation (ALMaQUEST) survey, we show that for star forming spaxels in the main sequence galaxies, the three local quantities, star-formation rate surface density (sigsfr), stellar mass surface density (sigsm), and the h2~mass surface density (sigh2), are strongly correlated with one another and form a 3D linear (in log) relation with dispersion. In addition to the two well known scaling relations, the resolved SFMS (sigsfr~ vs. sigsm) and the Schmidt-Kennicutt relation (sigsfr~ vs. sigh2; SK relation), there is a third scaling relation between sigh2~ and sigsm, which we refer to as the `molecular gas main sequence (MGMS). The latter indicates that either the local gas mass traces the gravitational potential set by the local stellar mass or both quantities follow the underlying total mass distributions. The scatter of the resolved SFMS ($sigma sim 0.25$ dex) is the largest compared to those of the SK and MGMS relations ($sigma sim$ 0.2 dex). A Pearson correlation test also indicates that the SK and MGMS relations are more strongly correlated than the resolved SFMS. Our result suggests a scenario in which the resolved SFMS is the least physically fundamental and is the consequence of the combination of the SK and the MGMS relations.
Using data from four deep fields (COSMOS, AEGIS, ECDFS, and CDFN), we study the correlation between the position of galaxies in the star formation rate (SFR) versus stellar mass plane and local environment at $z<1.1$. To accurately estimate the galaxy SFR, we use the deepest available Spitzer/MIPS 24 and Herschel/PACS datasets. We distinguish group environments ( $M_{halo}sim$10$^{12.5-14.2}$$M_{odot}$) based on the available deep X-ray data and lower halo mass environments based on the local galaxy density. We confirm that the Main Sequence (MS) of star forming galaxies is not a linear relation and there is a flattening towards higher stellar masses ( $M_*>10^{10.4-10.6}$ $M_{odot}$), across all environments. At high redshift ( $0.5<z<1.1$ ), the MS varies little with environment. At low redshift ( $0.15<z<0.5$ ), group galaxies tend to deviate from the mean MS towards the region of quiescence with respect to isolated galaxies and less-dense environments. We find that the flattening of the MS toward low SFR is due to an increased fraction of bulge dominated galaxies at high masses. Instead, the deviation of group galaxies from the MS at low redshift is caused by a large fraction of red disk dominated galaxies which are not present in the lower density environments. Our results suggest that above a mass threshold ( $sim10^{10.4}-10^{10.6}$$M_{odot}$ ) stellar mass, morphology and environment act together in driving the evolution of the SF activity towards lower level. The presence of a dominating bulge and the associated quenching processes are already in place beyond $zsim$1. The environmental effects appear, instead, at lower redshifts and have a long time-scale.
Deep far-infrared (FIR) cosmological surveys are known to be affected by source confusion, causing issues when examining the main sequence (MS) of star forming galaxies. This has typically been partially tackled by the use of stacking. However, stacking only provides the average properties of the objects in the stack. This work aims to trace the MS over $0.2leq z<6.0$ using the latest de-blended Herschel photometry, which reaches $approx10$ times deeper than the 5$sigma$ confusion limit in SPIRE. This provides more reliable star formation rates (SFRs), especially for the fainter galaxies, and hence a more reliable MS. We built a pipeline that uses the spectral energy distribution (SED) modelling and fitting tool CIGALE to generate flux density priors in the Herschel SPIRE bands. These priors were then fed into the de-blending tool XID+ to extract flux densities from the SPIRE maps. Multi-wavelength data were combined with the extracted SPIRE flux densities to constrain SEDs and provide stellar mass (M$_{star}$) and SFRs. These M$_{star}$ and SFRs were then used to populate the SFR-M$_{star}$ plane over $0.2leq z<6.0$. No significant evidence of a high-mass turn-over was found; the best fit is thus a simple two-parameter power law of the form log(SFR)$=alpha[$log(M$_{star})-10.5]+beta$. The normalisation of the power law increases with redshift, rapidly at $zlesssim1.8$, from $0.58pm0.09$ at $zapprox0.37$ to $1.31pm0.08$ at $zapprox1.8$. The slope is also found to increase with redshift, perhaps with an excess around $1.8leq z<2.9$. The increasing slope indicates that galaxies become more self-similar as redshift increases, implying that the specific SFR of high-mass galaxies increases over $z=0.2$ to $z=6.0$, becoming closer to that of low-mass galaxies. The excess in the slope at $1.8leq z<2.9$, if present, coincides with the peak of the cosmic star formation history.
The analytic equilibrium model for galaxy evolution using a mass balance equation is able to reproduce mean observed galaxy scaling relations between stellar mass, halo mass, star formation rate (SFR) and metallicity across the majority of cosmic time with a small number of parameters related to feedback. Here we aim to test this data-constrained model to quantify deviations from the mean relation between stellar mass and SFR, i.e. the star-forming galaxy main sequence (MS). We implement fluctuation in halo accretion rates parameterised from merger-based simulations, and quantify the intrinsic scatter introduced into the MS under the assumption that fluctuations in star formation follow baryonic inflow fluctuations. We predict the 1-sigma MS scatter to be ~ 0.2 - 0.25 dex over the stellar mass range 10^8 Mo to 10^11 Mo and a redshift range 0.5 < z < 3 for SFRs averaged over 100 Myr. The scatter increases modestly at z > 3, as well as by averaging over shorter timescales. The contribution from merger-induced star formation is generally small, around 5% today and 10 - 15% during the peak epoch of cosmic star formation. These results are generally consistent with available observations, suggesting that deviations from the MS primarily reflect stochasticity in the inflow rate owing to halo mergers.
By using a set of different SFR indicators, including WISE mid-infrared and Halpha emission, we study the slope of the Main Sequence (MS) of local star forming galaxies at stellar masses larger than 10^{10} M_{odot}. The slope of the relation strongly depends on the SFR indicator used. In all cases, the local MS shows a bending at high stellar masses with respect to the slope obtained in the low mass regime. While the distribution of galaxies in the upper envelope of the MS is consistent with a log-normal distribution, the lower envelope shows an excess of galaxies, which increases as a function of the stellar mass but varies as a function of the SFR indicator used. The scatter of the best log-normal distribution increases with stellar mass from ~0.3 dex at 10^{10} M_{odot} to ~0.45 at 10^{11} M_{odot}. The MS high-mass end is dominated by central galaxies of group sized halos with a red bulge and a disk redder than the lower mass counterparts. We argue that the MS bending in this region is due to two processes: i) the formation of a bulge component as a consequence of the increased merger activity in groups, and ii) the cold gas starvation induced by the hot halo environment, which cuts off the gas inflow onto the disk. Similarly, the increase of the MS scatter at high stellar masses would be explained by the larger spread of star formation histories of central group and cluster galaxies with respect to lower mass systems.