Star Formation Rate-stellar Mass Relation Research Articles

The Art of Measuring Physical Parameters in Galaxies: A Critical Assessment of Spectral Energy Distribution Fitting Techniques

The study of galaxy evolution hinges on our ability to interpret multiwavelength galaxy observations in terms of their physical properties. To do this, we rely on spectral energy distribution (SED) models, which allow us to infer physical parameters from spectrophotometric data. In recent years, thanks to wide and deep multiwave band galaxy surveys, the volume of high-quality data have significantly increased. Alongside the increased data, algorithms performing SED fitting have improved, including better modeling prescriptions, newer templates, and more extensive sampling in wavelength space. We present a comprehensive analysis of different SED-fitting codes including their methods and output with the aim of measuring the uncertainties caused by the modeling assumptions. We apply 14 of the most commonly used SED-fitting codes on samples from the CANDELS photometric catalogs at z ∼ 1 and z ∼ 3. We find agreement on the stellar mass, while we observe some discrepancies in the star formation rate (SFR) and dust-attenuation results. To explore the differences and biases among the codes, we explore the impact of the various modeling assumptions as they are set in the codes (e.g., star formation histories, nebular, dust and active galactic nucleus models) on the derived stellar masses, SFRs, and A V values. We then assess the difference among the codes on the SFR–stellar mass relation and we measure the contribution to the uncertainties by the modeling choices (i.e., the modeling uncertainties) in stellar mass (∼0.1 dex), SFR (∼0.3 dex), and dust attenuation (∼0.3 mag). Finally, we present some resources summarizing best practices in SED fitting.

The Roles of Morphology and Environment on the Star Formation Rate–Stellar Mass Relation in COSMOS from 0 < z < 3.5

We investigate the relationship between environment, morphology, and the star formation rate (SFR)–stellar mass relation derived from a sample of star-forming (SF) galaxies (commonly referred to as the “star formation main sequence”, SFMS) in the COSMOS field from 0 < z < 3.5. We constructed and fit the far-UV–far-IR spectral energy distributions of our stellar-mass-selected sample of 111,537 galaxies with stellar and dust emission models using the public packages MAGPHYS and SED3FIT. From the best-fit parameter estimates, we construct the SFR–stellar mass relation as a function of redshift, local environment, NUVrJ color diagnostics, and morphology. We find that the shape of the main sequence derived from our color–color and specific-star-formation-rate-selected SF galaxy population, including the turnover at high stellar mass, does not exhibit an environmental dependence at any redshift from 0 < z < 3.5. We investigate the role of morphology in the high-mass end of the SFMS to determine whether bulge growth is driving the high-mass turnover. We find that SF galaxies experience this turnover independent of bulge-to-total ratio, strengthening the case that the turnover is due to the disk component’s specific SFR evolving with stellar mass rather than bulge growth.

Constraining galactic baryon cycle using the galaxy stellar-to-halo mass relations

ABSTRACT Galaxies display several well-behaved scaling relations between their properties, such as the star formation rate–stellar mass relation (the main sequence, MS) and the stellar mass–halo mass relation (SHMR). In principle, these scaling relations could imply different star formation histories (SFHs) of galaxies and different constraints on galaxy formation physics. In this paper, we derive the SFHs of galaxies by assuming that they always follow the SHMRs at different redshifts and use an empirical model to constrain key processes in their baryon cycle. It is found that, besides cold accretion due to halo growth, outflow of gas produced by stellar feedback has to be recycled to sustain the derived SFHs of galaxies. The recycled fraction is strongly affected by the baryon fraction in accreted low-mass haloes and the mass loading factor that quantifies the ratio between the galactic outflow rate and star formation rate. Our fiducial model predicts that around 20–60 per cent of outflow is recycled in $\sim 0.5\!-\!4\, \mathrm{Gyr}$, while simulations predict a slightly higher recycle fraction and a lower recycle time. We argue that strong constraints on the baryon cycle process can be obtained from future observation of the circum-galactic medium (CGM) of galaxies, such as the gas cooling rate of CGM. We also find that the implied SFHs from the SHMRs indicate that galaxies stay on the MS only for part of their lifetimes. Our model reproduces the evolution of the mass–metallicity relation as well.

The spatially resolved view of star formation in galaxy clusters

AbstractIntegral field spectroscopic studies of galaxies in dense environments, such as clusters and groups of galaxies, have provided new insights for understanding how star formation proceeds, and quenches. I present the spatially resolved view of the star formation activity and its link with the multiphase gas in cluster galaxies based on MUSE and multi-wavelength data of the GASP survey. I discuss the link among the different scales (i.e. the link between the spatially resolved and the global star formation rate-stellar mass relation), the spatially resolved signatures and the quenching histories of jellyfish (progenitors) and post-starburst (descendants) galaxies in clusters. Finally, I discuss the multi-wavelength view of star-forming clumps both in galaxy disks and in the tails of stripped gas.

The H α star formation main sequence in cluster and field galaxies at z ∼ 1.6

ABSTRACT We calculate H α-based star formation rates and determine the star formation rate–stellar mass relation for members of three Spitzer Adaptation of the Red-Sequence Cluster Survey (SpARCS) clusters at z ∼ 1.6 and serendipitously identified field galaxies at similar redshifts to the clusters. We find similar star formation rates in cluster and field galaxies throughout our range of stellar masses. The results are comparable to those seen in other clusters at similar redshifts, and consistent with our previous photometric evidence for little quenching activity in clusters. One possible explanation for our results is that galaxies in our z ∼ 1.6 clusters have been accreted too recently to show signs of environmental quenching. It is also possible that the clusters are not yet dynamically mature enough to produce important environmental quenching effects shown to be important at low redshift, such as ram-pressure stripping or harassment.

The Roles of Mass and Environment in the Quenching of Galaxies. II.

Abstract We take advantage of an analytic model of galaxy formation coupled to the merger tree of an N-body simulation to study the roles of environment and stellar mass in the quenching of galaxies. The model has been originally set in order to provide the observed evolution of the stellar mass function as well as reasonable predictions of the star formation rate–stellar mass relation, from high redshift to the present time. We analyze the stellar mass and environmental quenching efficiencies and their dependence on stellar mass, halo mass (taken as a proxy for the environment), and redshift. Our analysis shows that the two quenching efficiencies are dependent on redshift and stellar and halo mass, and that the halo mass is also a good proxy for the environment. The environmental quenching increases with decreasing redshift and is inefficient below log M * ∼ 9.5, reaches the maximum value at log M * ∼ 10.5, and decreases again, becoming poorly efficient at very high stellar mass (log M * ≳ 11.5). Central and satellites galaxies are mass quenched differently: for the former, the quenching efficiency depends very weakly on redshift but strongly on stellar mass; for the latter, it strongly depends on both stellar mass and redshift in the range . According to the most recent observational results, we find that the two quenching efficiencies are not separable: intermediate-mass galaxies, as well as intermediate/massive galaxies in more massive halos, are environmentally quenched faster. At stellar masses lower than log M * ≲ 9.5, both quenching mechanisms become inefficient, independently of the redshift.

The origin of scatter in the star formation rate–stellar mass relation

Observations have revealed that the star formation rate (SFR) and stellar mass (M$_{\rm star}$) of star-forming galaxies follow a tight relation known as the galaxy main sequence. However, what physical information is encoded in this relation is under debate. Here, we use the EAGLE cosmological hydrodynamical simulation to study the mass dependence, evolution and origin of scatter in the SFR-M$_{\rm star}$ relation. At $z=0$, we find that the scatter decreases slightly with stellar mass from 0.35 dex at M$_{\rm star} \approx 10^9$ M$_{\odot}$ to 0.30 dex at M$_{\rm star} \gtrsim 10^{10.5}$ M$_{\odot}$. The scatter decreases from $z=0$ to $z=5$ by 0.05 dex at M$_{\rm star} \gtrsim 10^{10}$ M$_{\odot}$ and by 0.15 dex for lower masses. We show that the scatter at $z=0.1$ originates from a combination of fluctuations on short time-scales (ranging from 0.2-2 Gyr) that are presumably associated with self-regulation from cooling, star formation and outflows, but is dominated by long time-scale ($\sim 10$ Gyr) variations related to differences in halo formation times. Shorter time-scale fluctuations are relatively more important for lower-mass galaxies. At high masses, differences in black hole formation efficiency cause additional scatter, but also diminish the scatter caused by different halo formation times. While individual galaxies cross the main sequence multiple times during their evolution, they fluctuate around tracks associated with their halo properties, i.e. galaxies above/below the main sequence at $z = 0.1$ tend to have been above/below the main sequence for $\gg1$ Gyr.

SDSS IV MaNGA: Dependence of Global and Spatially Resolved SFR–M∗ Relations on Galaxy Properties

Abstract The galaxy integrated Hα star formation rate–stellar mass relation, or SFR(global)–M *(global) relation, is crucial for understanding star formation history and evolution of galaxies. However, many studies have dealt with SFR using unresolved measurements, which makes it difficult to separate out the contamination from other ionizing sources, such as active galactic nuclei and evolved stars. Using the integral field spectroscopic observations from SDSS-IV MaNGA, we spatially disentangle the contribution from different Hα powering sources for ∼1000 galaxies. We find that, when including regions dominated by all ionizing sources in galaxies, the spatially resolved relation between Hα surface density (ΣHα (all)) and stellar mass surface density (Σ*(all)) progressively turns over at the high Σ*(all) end for increasing M *(global) and/or bulge dominance (bulge-to-total light ratio, B/T). This in turn leads to the flattening of the integrated Hα(global)–M *(global) relation in the literature. By contrast, there is no noticeable flattening in both integrated Hα(H ii)–M *(H ii) and spatially resolved ΣHα (H ii)–Σ*(H ii) relations when only regions where star formation dominates the ionization are considered. In other words, the flattening can be attributed to the increasing regions powered by non-star-formation sources, which generally have lower ionizing ability than star formation. An analysis of the fractional contribution of non-star-formation sources to total Hα luminosity of a galaxy suggests a decreasing role of star formation as an ionizing source toward high-mass, high-B/T galaxies and bulge regions. This result indicates that the appearance of the galaxy integrated SFR–M * relation critically depends on their global properties (M *(global) and B/T) and relative abundances of various ionizing sources within the galaxies.

The stellar mass function and star formation rate–stellar mass relation of galaxies at z ∼ 4–7

We investigate the evolution of the star formation rate-stellar mass relation (SFR-M*) and Galaxy Stellar Mass Function (GSMF) of z ~ 4-7 galaxies, using cosmological simulations run with the smoothed particle hydrodynamics code P-GADGET3(XXL). We explore the effects of different feedback prescriptions (supernova driven galactic winds and AGN feedback), initial stellar mass functions and metal cooling. We show that our fiducial model, with strong energy-driven winds and early AGN feedback, is able to reproduce the observed stellar mass function obtained from Lyman-break selected samples of star forming galaxies at redshift 6 < z < 7. At z ~ 4, observed estimates of the GSMF vary according to how the sample was selected. Our simulations are more consistent with recent results from K-selected samples, which provide a better proxy of stellar masses and are more complete at the high mass end of the distribution. We find that in some cases simulated and observed SFR-M* relations are in tension, and this can lead to numerical predictions for the GSMF in excess of the GSMF observed. By combining the simulated SFR(M*) relationship with the observed star formation rate function at a given redshift, we argue that this disagreement may be the result of the uncertainty in the SFR-M* (Luv-M*) conversion. Our simulations predict a population of faint galaxies not seen by current observations.

The star formation main sequence and stellar mass assembly of galaxies in the Illustris simulation

Understanding the physical processes that drive star formation is a key challenge for galaxy formation models. In this paper, we study the tight correlation between the star formation rate (SFR) and stellar mass of galaxies at a given redshift, how halo growth influences star formation, and star formation histories of individual galaxies. We study these topics using Illustris, a state-of-the-art cosmological hydrodynamical simulation of galaxy formation. Illustris reproduces the observed relation (the star formation main sequence, SFMS) between SFR and stellar mass at redshifts z = 0 and 4, but at intermediate redshifts of z ≃ 1–2, the simulated SFMS has a significantly lower normalization than reported by observations. The scatter in the relation is consistent with the observed scatter. However, the fraction of outliers above the SFR–stellar mass relation in Illustris is less than that observed. Galaxies with halo masses of ∼1012 M⊙ dominate the SFR density of the Universe, in agreement with the results of abundance matching. Furthermore, more-massive galaxies tend to form the bulk of their stars at high redshift, which indicates that ‘downsizing’ occurs in Illustris. We also studied the star formation histories of individual galaxies, including the use of a principal component analysis decomposition. We find that for fixed stellar mass, galaxies that form earlier have more-massive black holes at z = 0, indicating that star formation and black hole growth are tightly linked processes in Illustris. While many of the properties of normal star-forming galaxies are well reproduced in the Illustris simulation, forming a realistic population of starbursts will likely require higher resolution and probably a more sophisticated treatment of star formation and feedback from stars and black holes.

How to distinguish starbursts and quiescently star-forming galaxies: the ‘bimodal’ submillimetre galaxy population as a case study

In recent work (arXiv:1101.0002) we have suggested that the high-redshift (z ~ 2-4) bright submillimetre galaxy (SMG) population is heterogeneous, with major mergers contributing both at early stages, where quiescently star-forming discs are blended into one submm source (`galaxy-pair SMGs'), and late stages, where mutual tidal torques drive gas inflows and cause strong starbursts. Here we combine hydrodynamic simulations of major mergers with 3-D dust radiative transfer calculations to determine observational diagnostics that can distinguish between quiescently star-forming SMGs and starburst SMGs via integrated data alone. We fit the far-IR SEDs of the simulated galaxies with the optically thin single-temperature modified blackbody, the full form of the single-temperature modified blackbody, and a power-law temperature-distribution model. The effective dust temperature, T_dust, and power-law index of the dust emissivity in the far-IR, \beta, derived can significantly depend on the fitting form used, and the intrinsic \beta\ of the dust is not recovered. However, for all forms used here, there is a T_dust above which almost all simulated galaxies are starbursts, so a T_dust cut is very effective at selecting starbursts. Simulated merger-induced starbursts also have higher L_IR/M_gas and L_IR/L_FUV than quiescently star-forming galaxies and lie above the star formation rate-stellar mass relation. These diagnostics can be used to test our claim that the SMG population is heterogeneous and to observationally determine what star formation mode dominates a given galaxy population. We comment on applicability of these diagnostics to ULIRGs that would not be selected as SMGs. These `hot-dust ULIRGs' are typically starburst galaxies lower in mass than SMGs, but they can also simply be SMGs observed from a different viewing angle.