Search

Recent high-resolution interferometric images of submillimetre galaxies (SMGs) reveal fascinatingly complex morphologies. This raises a number of questions: how does the relative orientation of a galaxy affect its observed submillimetre emission, and does this result in an 'orientation bias' in the selection and analysis of such galaxies in flux-limited cosmological surveys? We investigated these questions using the simba cosmological simulation paired with the dust radiative transfer code powderday. We selected eight simulated SMGs (S 850 ≳ 2 mJy) at z  = 2, and measured the variance of their 'observed' emission over 50 random orientations. Each galaxy exhibits significant scatter in its emission close to the peak of the thermal dust emission, with variation in flux density of up to a factor of 2.7. This results in an appreciable dispersion in the inferred dust temperatures and infrared luminosities (16th–84th percentile ranges of 5 K and 0.1 dex, respectively) and therefore a fundamental uncertainty in derived parameters such as dust mass and star formation rate (∼30 per cent for the latter using simple calibrations). Using a Monte Carlo simulation we also assessed the impact of orientation on flux-limited surveys, finding a bias in the selection of SMGs towards those with face-on orientations, as well as those at lower redshifts. We predict that the orientation bias will affect flux-limited single-dish surveys, most significantly at THz frequencies, and this bias should be taken into account when placing the results of targeted follow-up studies in a statistical context. [ABSTRACT FROM AUTHOR]

We investigate the Baryonic Tully–Fisher relation (BTFR) in the |$(100\, h^{-1}{\rm Mpc})^3$| simba hydrodynamical galaxy formation simulation together with a higher resolution |$(25\, h^{-1}{\rm Mpc})^3$| simba run, for over 10 000 disc-dominated, H  i -rich galaxies. We generate simulated galaxy rotation curves from the mass distribution, which we show yields similar results to using the gas rotational velocities. From this, we measure the galaxy rotation velocity V circ using four metrics: |$V_{\rm max}, V_{\rm flat}, V_{2R_e},$| and V polyex. We compare the predicted BTFR to the SPARC observational sample and find broad agreement. In detail, however, simba is biased towards higher V circ by up to 0.1 dex. We find evidence for the flattening of the BTFR in V circ > 300  km s−1 galaxies, in agreement with recent observational findings. simba 's rotation curves are more peaked for lower mass galaxies, in contrast with observations, suggesting overly bulge-dominated dwarf galaxies in our sample. We investigate for residuals around the BTFR versus H  i mass, stellar mass, gas fraction, and specific star formation rate, which provide testable predictions for upcoming BTFR surveys. simba 's BTFR shows sub-optimal resolution mass convergence, with the higher resolution run lowering V in better agreement with data. [ABSTRACT FROM AUTHOR]

We present a Bayesian full-spectral-fitting analysis of 75 massive (⁠|$M_* \gt 10^{10.3} \, \mathrm{M_\odot }$|⁠) UVJ-selected galaxies at redshifts of 1.0 < z < 1.3, combining extremely deep rest-frame ultraviolet spectroscopy from VANDELS with multiwavelength photometry. By the use of a sophisticated physical plus systematic uncertainties model, constructed within the bagpipes code, we place strong constraints on the star-formation histories (SFHs) of individual objects. We first constrain the stellar mass versus stellar age relationship, finding a steep trend towards earlier average formation time with increasing stellar mass (downsizing) of |$1.48^{+0.34}_{-0.39}$| Gyr per decade in mass, although this shows signs of flattening at |$M_* \gt 10^{11} \, \mathrm{M_\odot }$|⁠. We show that this is consistent with other spectroscopic studies from 0 < z < 2. This relationship places strong constraints on the AGN-feedback models used in cosmological simulations. We demonstrate that, although the relationships predicted by simba and illustristng agree well with observations at z = 0.1, they are too shallow at z = 1, predicting an evolution of ≲0.5 Gyr per decade in mass. Secondly, we consider the connections between green-valley, post-starburst, and quiescent galaxies, using our inferred SFH shapes and the distributions of galaxy physical properties on the UVJ diagram. The majority of our lowest-mass galaxies (⁠|$M_* \sim 10^{10.5} \, \mathrm{M_\odot }$|⁠) are consistent with formation in recent (z < 2), intense starburst events, with time-scales of ≲500 Myr. A second class of objects experience extended star-formation epochs before rapidly quenching, passing through both green-valley and post-starburst phases. The most massive galaxies in our sample are extreme systems: already old by z = 1, they formed at z ∼ 5 and quenched by z = 3. However, we find evidence for their continued evolution through both AGN and rejuvenated star-formation activity. [ABSTRACT FROM AUTHOR]

We present the results of a study utilizing ultradeep, rest-frame UV, spectroscopy to quantify the relationship between stellar mass and stellar metallicity for 681 star-forming galaxies at 2.5 < z < 5.0 (〈 z 〉 = 3.5 ± 0.6) drawn from the VANDELS survey. Via a comparison with high-resolution stellar population synthesis models, we determine stellar metallicities (Z ∗, here a proxy for the iron abundance) for a set of high signal-to-noise ratio composite spectra formed from subsamples selected by mass and redshift. Across the stellar mass range |$8.5\,\, \lt\,\, \mathrm{log}(\langle M_{\ast } \rangle /\rm {M}_{\odot })\,\, \lt\,\, 10.2$|⁠ , we find a strong correlation between stellar metallicity (Z ∗/Z⊙) and stellar mass, with stellar metallicity monotonically increasing from Z ∗/Z⊙ < 0.09 at |$\langle M_{\ast } \rangle = 3.2 \times 10^{8}\, \rm {M}_{\odot }$| to Z ∗/Z⊙ = 0.27 at |$\langle M_{\ast } \rangle = 1.7 \times 10^{10}\, \rm {M}_{\odot }$|⁠. In contrast, at a given stellar mass, we find no evidence for significant metallicity evolution across the redshift range of our sample. However, comparing our results to the z  = 0 stellar mass–metallicity relation for star-forming galaxies, we find that the 〈 z 〉 = 3.5 relation is consistent with being shifted to lower metallicities by ≃0.6 dex at all stellar masses. Contrasting our derived stellar metallicities with estimates of the gas-phase metallicities of galaxies at similar redshifts and stellar masses, we find evidence for enhanced |$\rm {O}/\rm {Fe}$| ratios in z  ≳ 2.5 star-forming galaxies of the order (O/Fe) ≳ 1.8 × (O/Fe)⊙. Finally, by comparing our results to the predictions of three cosmological simulations, we find that the 〈 z 〉 = 3.5 stellar mass–metallicity relation is consistent with current predictions for how outflow strength scales with galaxy stellar mass. This conclusion is supported by an analysis of one-zone analytic chemical evolution models, and suggests that the mass-loading parameter (⁠|$\eta =\dot{M}_{\mathrm{outflow}}/M_{\ast }$|⁠) scales as |$\eta \propto M_{\ast }^{\beta }$| with β ≃ −0.4. [ABSTRACT FROM AUTHOR]

Observed galaxies with high stellar masses or in dense environments have low specific star formation rates, i.e. they are quenched. Based on cosmological hydrodynamic simulations that include a prescription where quenching occurs in regions dominated by hot (>105.4 K) gas, we argue that this hot gas quenching in haloes >1012 M⊙ drives both mass quenching (i.e. central quenching) and environment quenching (i.e. satellite quenching). These simulations reproduce a broad range of locally observed trends among quenching, halo mass, stellar mass, environment, and distance to halo centre. Mass quenching is independent of environment because ~1012-1013 M⊙ 'mass quenching haloes' inhabit a large range of environments. On the other hand, environment quenching is independent of stellar mass because galaxies of all stellar masses may live in dense environments as satellites of groups and clusters. As in observations, the quenched fraction of satellites increases with halo mass and decreases with distance to the centre of the group or cluster. We investigate pre-processing in group haloes, ejected former satellites, and hot gas that extends beyond the virial radius. The agreement of our model with key observational trends suggests that hot gas in massive haloes plays a leading role in quenching low-redshift galaxies. [ABSTRACT FROM AUTHOR]

We measure new estimates for the galaxy stellar mass function and star formation rates for samples of galaxies at z ∼ 4, 5, 6 and 7 using data in the CANDELS GOODS South field. The deep near-infrared observations allow us to construct the stellar mass function at z ≥ 6 directly for the first time. We estimate stellar masses for our sample by fitting the observed spectral energy distributions with synthetic stellar populations, including nebular line and continuum emission. The observed UV luminosity functions for the samples are consistent with previous observations; however, we find that the observed MUV-M* relation has a shallow slope more consistent with a constant mass-to-light ratio and a normalization which evolves with redshift. Our stellar mass functions have steep low-mass slopes (α ≈ −1.9), steeper than previously observed at these redshifts and closer to that of the UV luminosity function. Integrating our new mass functions, we find the observed stellar mass density evolves from $\log _{10} \rho _{*} = 6.64^{+0.58}_{-0.89}$ at z ∼ 7 to 7.36 ± 0.06 M⊙ Mpc− 3 at z ∼ 4. Finally, combining the measured UV continuum slopes (β) with their rest-frame UV luminosities, we calculate dust-corrected star formation rates (SFR) for our sample. We find the specific SFR for a fixed stellar mass increases with redshift whilst the global SFR density falls rapidly over this period. Our new SFR density estimates are higher than previously observed at this redshift. [ABSTRACT FROM AUTHOR]

ABSTRACT We examine the cosmic growth of the red sequence in a cosmological hydrodynamic simulation that includes a heuristic prescription for quenching star formation that yields a realistic passive galaxy population today. In this prescription, haloes dominated by hot gas are continually heated to prevent their coronae from fuelling new star formation. Hot coronae primarily form in haloes above ∼1012 M⊙, so that galaxies with stellar masses ∼1010.5 M⊙ are the first to be quenched and move on to the red sequence at z > 2. The red sequence is concurrently populated at low masses by satellite galaxies in large haloes that are starved of new fuel, resulting in a dip in passive galaxy number densities around ∼1010 M⊙. Stellar mass growth continues for galaxies even after joining the red sequence, primarily through minor mergers with a typical mass ratio ∼1:5. For the most massive systems, the size growth implied by the distribution of merger mass ratios is typically approximately two times the corresponding mass growth, consistent with observations. This model reproduces mass-density and colour-density trends in the local Universe, with essentially no evolution to z = 1, with the hint that such relations may be washed out by z ∼ 2. Simulated galaxies are increasingly likely to be red at high masses or high local overdensities. In our model, the presence of surrounding hot gas drives the trends with both mass and environment. [ABSTRACT FROM AUTHOR]

Galaxies with stellar masses near M* contain the majority of stellar mass in the universe, and are therefore of special interest in the study of galaxy evolution. The Milky Way (MW) and Andromeda (M31) have present-day stellar masses near M*, at 5 × 1010M☼ (defined here to be MW-mass) and 1011M☼ (defined to be M31-mass). We study the typical progenitors of these galaxies using the FOURSTAR Galaxy Evolution Survey (ZFOURGE). ZFOURGE is a deep medium-band near-IR imaging survey, which is sensitive to the progenitors of these galaxies out to z ∼ 3. We use abundance-matching techniques to identify the main progenitors of these galaxies at higher redshifts. We measure the evolution in the stellar mass, rest-frame colors, morphologies, far-IR luminosities, and star formation rates, combining our deep multiwavelength imaging with near-IR Hubble Space Telescope imaging from Cosmic Near-IR Deep Extragalactic Legacy Survey (CANDELS), and Spitzer and Herschel far-IR imaging from Great Observatories Origins Deep Survey-Herschel and CANDELS-Herschel. The typical MW-mass and M31-mass progenitors passed through the same evolution stages, evolving from blue, star-forming disk galaxies at the earliest stages to redder dust-obscured IR-luminous galaxies in intermediate stages and to red, more quiescent galaxies at their latest stages. The progenitors of the MW-mass galaxies reached each evolutionary stage at later times (lower redshifts) and with stellar masses that are a factor of two to three lower than the progenitors of the M31-mass galaxies. The process driving this evolution, including the suppression of star formation in present-day M* galaxies, requires an evolving stellar-mass/halo-mass ratio and/or evolving halo-mass threshold for quiescent galaxies. The effective size and SFRs imply that the baryonic cold-gas fractions drop as galaxies evolve from high redshift to z ∼ 0 and are strongly anticorrelated with an increase in the Sérsic index. Therefore, the growth of galaxy bulges in M* galaxies corresponds to a rapid decline in the galaxy gas fractions and/or a decrease in the star formation efficiency. [ABSTRACT FROM AUTHOR]

We present measurements of the specific angular momentum content (j) of galaxies drawn from the Simba cosmological hydrodynamic simulations. For the stellar, H i and baryonic matter components we demonstrate the existence of extremely tight relations between j and the mass contained within the radius at which the H i mass surface density decreases to 1 M ⊙ pc−2. These relations are broadly consistent with a variety of empirical measurements. We confirm the observational result that the scatter in the stellar j – M relation is driven largely by H i content, and measure the dependence of its scatter on the deviations of galaxies from other important scaling relations. For a given stellar mass, H i -rich/poor galaxies have more/less-than-average stellar specific angular momentum. A similar, yet weaker, correlation exists for H i mass fraction. Overall, our results demonstrate the utility of the Simba simulations as a platform for understanding and contextualising the data and results from forthcoming large galaxy surveys. [ABSTRACT FROM AUTHOR]