No Arabic abstract
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.
We present a framework for modelling the star-formation histories of galaxies as a stochastic process. We define this stochastic process through a power spectrum density with a functional form of a broken power-law. Star-formation histories are correlated on short timescales, the strength of this correlation described by a power-law slope, $alpha$, and they decorrelate to resemble white noise over a timescale that is proportional to the timescale of the break in the power spectrum density, $tau_{rm break}$. We use this framework to explore the properties of the stochastic process that, we assume, gives rise to the log-normal scatter about the relationship between star-formation rate and stellar mass, the so-called galaxy star-forming main sequence. Specifically, we show how the measurements of the normalisation and width ($sigma_{rm MS}$) of the main sequence, measured in several passbands that probe different timescales, give a constraint on the parameters of the underlying power spectrum density. We first derive these results analytically for a simplified case where we model observations by averaging over the recent star-formation history. We then run numerical simulations to find results for more realistic observational cases. As a proof of concept, we use observational estimates of the main sequence scatter at $zsim0$ and $M_{star}approx10^{10}~M_{odot}$ measured in H$alpha$, UV+IR and the u-band, and show that combination of these point to $tau_{rm break}=178^{+104}_{-66}$ Myr, when assuming $alpha=2$. This implies that star-formation histories of galaxies lose memory of their previous activity on a timescale of $sim200$ Myr, highlighting the importance of baryonic effects that act over the dynamical timescales of galaxies.
Using a sample of 11,478 spaxels in 34 galaxies with molecular gas, star formation and stellar maps taken from the ALMA-MaNGA QUEnching and STar formation (ALMaQUEST) survey, we investigate the parameters that correlate with variations in star formation rates on kpc scales. We use a combination of correlation statistics and an artificial neural network to quantify the parameters that drive both the absolute star formation rate surface density (Sigma_SFR), as well as its scatter around the resolved star forming main sequence (Delta Sigma_SFR). We find that Sigma_SFR is primarily regulated by molecular gas surface density (Sigma_H2) with a secondary dependence on stellar mass surface density (Sigma_*), as expected from an `extended Kennicutt-Schmidt relation. However, Delta Sigma_SFR is driven primarily by changes in star formation efficiency (SFE), with variations in gas fraction playing a secondary role. Taken together, our results demonstrate that whilst the absolute rate of star formation is primarily set by the amount of molecular gas, the variation of star formation rate above and below the resolved star forming main sequence (on kpc scales) is primarily due to changes in SFE.
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.
Using a compilation of 25 studies from the literature, we investigate the evolution of the star-forming galaxy (SFG) Main Sequence (MS) in stellar mass and star formation rate (SFR) out to $z sim 6$. After converting all observations to a common set of calibrations, we find a remarkable consensus among MS observations ($sim 0.1$ dex 1$sigma$ interpublication scatter). By fitting for time evolution of the MS in bins of constant mass, we deconvolve the observed scatter about the MS within each observed redshift bins. After accounting for observed scatter between different SFR indicators, we find the width of the MS distribution is $sim 0.2$ dex and remains constant over cosmic time. Our best fits indicate the slope of the MS is likely time-dependent, with our best fit $logtextrm{SFR}(M_*,t) = left(0.84 pm 0.02 - 0.026 pm 0.003 times tright) log M_* - left(6.51 pm 0.24 - 0.11 pm 0.03 times tright)$, with $t$ the age of the Universe in Gyr. We use our fits to create empirical evolutionary tracks in order to constrain MS galaxy star formation histories (SFHs), finding that (1) the most accurate representations of MS SFHs are given by delayed-$tau$ models, (2) the decline in fractional stellar mass growth for a typical MS galaxy today is approximately linear for most of its lifetime, and (3) scatter about the MS can be generated by galaxies evolving along identical evolutionary tracks assuming an initial $1sigma$ spread in formation times of $sim 1.4$ Gyr.
We compare various star formation rate (SFR) indicators for star-forming galaxies at $1.4<z<2.5$ in the COSMOS field. The main focus is on the SFRs from the far-IR (PACS-Herschel data) with those from the ultraviolet, for galaxies selected according to the BzK criterion. FIR-selected samples lead to a vastly different slope of the SFR-stellar mass ($M_*$) relation, compared to that of the dominant main sequence population as measured from the UV, since the FIR selection picks predominantly only a minority of outliers. However, there is overall agreement between the main sequences derived with the two SFR indicators, when stacking on the PACS maps the BzK-selected galaxies. The resulting logarithmic slope of the SFR-{$M_*$} relation is $sim0.8-0.9$, in agreement with that derived from the dust-corrected UV-luminosity. Exploiting deeper 24$mu$m-Spitzer data we have characterized a sub-sample of galaxies with reddening and SFRs poorly constrained, as they are very faint in the $B$ band. The combination of Herschel with Spitzer data have allowed us to largely break the age/reddening degeneracy for these intriguing sources, by distinguishing whether a galaxy is very red in B-z because of being heavily dust reddened, or whether because star formation has been (or is being) quenched. Finally, we have compared our SFR(UV) to the SFRs derived by stacking the radio data and to those derived from the H$alpha$ luminosity of a sample of star-forming galaxies at $1.4<z<1.7$. The two sets of SFRs are broadly consistent as they are with the SFRs derived from the UV and by stacking the corresponding PACS data in various mass bins.