We present the evolution of galaxy sizes, from redshift 2 to 0, for actively star forming and passive galaxies in the cosmological hydrodynamical 100 super( 3) cMpc3 simulation of the EAGLE project. ...We find that the sizes increase with stellar mass, but that the relation weakens with increasing redshift. Separating galaxies by their star formation activity, we find that passive galaxies are typically smaller than active galaxies at a fixed stellar mass. These trends are consistent with those found in observations and the level of agreement between the predicted and observed size-mass relations is of the order of 0.1 dex for z < 1 and 0.2-0.3 dex from redshift 1 to 2. We use the simulation to compare the evolution of individual galaxies with that of the population as a whole. While the evolution of the size-stellar mass relation for active galaxies provides a good proxy for the evolution of individual galaxies, the evolution of individual passive galaxies is not well represented by the observed size-mass relation due to the evolving number density of passive galaxies. Observations of z ~ 2 galaxies have revealed an abundance of massive red compact galaxies, which depletes below z ~ 1. We find that a similar population forms naturally in the simulation. Comparing these galaxies with their z = 0 descendants, we find that all compact galaxies grow in size due to the high-redshift stars migrating outwards. Approximately 60 per cent of the compact galaxies increase in size further due to renewed star formation and/or mergers.
We investigate the evolution of galaxy masses and star formation rates in the Evolution and Assembly of Galaxies and their Environment (eagle) simulations. These comprise a suite of hydrodynamical ...simulations in a Λ cold dark matter cosmogony with subgrid models for radiative cooling, star formation, stellar mass-loss and feedback from stars and accreting black holes. The subgrid feedback was calibrated to reproduce the observed present-day galaxy stellar mass function and galaxy sizes. Here, we demonstrate that the simulations reproduce the observed growth of the stellar mass density to within 20 per cent. The simulations also track the observed evolution of the galaxy stellar mass function out to redshift z = 7, with differences comparable to the plausible uncertainties in the interpretation of the data. Just as with observed galaxies, the specific star formation rates of simulated galaxies are bimodal, with distinct star forming and passive sequences. The specific star formation rates of star-forming galaxies are typically 0.2 to 0.5 dex lower than observed, but the evolution of the rates track the observations closely. The unprecedented level of agreement between simulation and data across cosmic time makes eagle a powerful resource to understand the physical processes that govern galaxy formation.
Feedback from energy liberated by gas accretion on to black holes (BHs) is an attractive mechanism to explain the exponential cut-off at the massive end of the galaxy stellar mass function. Most ...previous implementations of BH accretion in hydrodynamical simulations of galaxy formation have assumed that BHs grow at an accretion rate that is proportion to the Bondi rate. A major concern is that the Bondi accretion rate is inappropriate when the accreting material has significant angular momentum. We present an improved accretion model that takes into account the circularization and subsequent viscous transport of infalling material, and implemented as a ‘subgrid’ model in hydrodynamic simulations. The resulting accretion rates are generally low in low mass (≲ 1011.5 M⊙) haloes, but show outbursts of Eddington-limited accretion during galaxy mergers. During outbursts these objects strongly resemble quasars. In higher mass haloes, gas accretion peaks at ∼10 per cent of the Eddington rate, which is thought to be conducive to the formation of radio jets. The resulting accretion rate depends strongly on the effective pressure of the gas surrounding the BH, which in turn depends strongly on halo mass. This induces a sharp transition in the importance of BH feedback. In small haloes, the growth of galaxies is regulated by star formation and supernova feedback, but above a halo mass of 1011.5 M⊙, rapid BH growth leads to the suppression of star formation and reduced growth of stellar mass with increasing halo mass.
Abstract
We present sensitive 870 μm continuum measurements from our ALMA programmes of 114 X-ray selected active galactic nuclei (AGN) in the Chandra Deep Field-South and Cosmic Evolution Survey ...fields. We use these observations in combination with data from Spitzer and Herschel to construct a sample of 86 X-ray selected AGN, 63 with ALMA constraints at z = 1.5–3.2 with stellar mass >2 × 1010 M⊙. We constructed broad-band spectral energy distributions in the infrared band (8–1000 μm) and constrain star-formation rates (SFRs) uncontaminated by the AGN. Using a hierarchical Bayesian method that takes into account the information from upper limits, we fit SFR and specific SFR (sSFR) distributions. We explore these distributions as a function of both X-ray luminosity and stellar mass. We compare our measurements to two versions of the Evolution and Assembly of GaLaxies and their Environments (EAGLE) hydrodynamical simulations: the reference model with AGN feedback and the model without AGN. We find good agreement between the observations and that predicted by the EAGLE reference model for the modes and widths of the sSFR distributions as a function of both X-ray luminosity and stellar mass; however, we found that the EAGLE model without AGN feedback predicts a significantly narrower width when compared to the data. Overall, from the combination of the observations with the model predictions, we conclude that (1) even with AGN feedback, we expect no strong relationship between the sSFR distribution parameters and instantaneous AGN luminosity and (2) a signature of AGN feedback is a broad distribution of sSFRs for all galaxies (not just those hosting an AGN) with stellar masses above ≈1010 M⊙.
ABSTRACT
We present the goals, design, and first results of the MUSE Ultra Deep Field (MUDF) survey, a large programme using the Multi Unit Spectroscopic Explorer (MUSE) instrument at the ESO Very ...Large Telescope. The MUDF survey is collecting ≈150 h on-source of integral field optical spectroscopy in a 1.5 × 1.2 arcmin2 region which hosts several astrophysical structures along the line of sight, including two bright z ≈ 3.2 quasars with close separation (≈500 kpc). Following the description of the data reduction procedures, we present the analysis of the galaxy environment and gaseous properties of seven groups detected at redshifts 0.5 < z < 1.5, spanning a large dynamic range in halo mass, $\log (M_h/\rm {M_\odot }) \approx 11 - 13.5$. For four of the groups, we find associated Mg ii absorbers tracing cool gas in high-resolution spectroscopy of the two quasars, including one case of correlated absorption in both sightlines at distance ≈480 kpc. The absorption strength associated with the groups is higher than what has been reported for more isolated galaxies of comparable mass and impact parameters. We do not find evidence for widespread cool gas giving rise to strong absorption within these groups. Combining these results with the distribution of neutral and ionized gas seen in emission in lower redshift groups, we conclude that gravitational interactions in the group environment strip gas from the galaxy haloes into the intragroup medium, boosting the cross-section of cool gas and leading to the high fraction of strong Mg ii absorbers that we detect.
We use the Galaxies-Intergalactic Medium Interaction Calculation (GIMIC) suite of cosmological hydrodynamical simulations to study the global structure and kinematics of stellar spheroids of Milky ...Way mass disc galaxies. Font et al. have recently demonstrated that these simulations are able to successfully reproduce the satellite luminosity functions and the metallicity and surface brightness profiles of the spheroids of the Milky Way and M31. A key to the success of the simulations is a significant contribution to the spheroid from stars that formed in situ. While the outer halo is dominated by accreted stars, stars formed in the main progenitor of the galaxy dominate at r≲ 30 kpc. In the present study, we show that this component was primarily formed in a protodisc at high redshift and was subsequently liberated from the disc by dynamical heating associated with mass accretion. As a consequence of its origin, the in situ component of the spheroid has different kinematics (namely net prograde rotation with respect to the disc) than that of the spheroid component built from the disruption of satellites. In addition, the in situ component has a flattened distribution, which is due in part to its rotation. We make comparisons with measurements of the shape and kinematics of local galaxies, including the Milky Way and M31, and stacked observations of more distant galaxies. We find that the simulated disc galaxies have spheroids of the correct shape (oblate with a median axial ratio of ∼0.6 at radii of ≲30 kpc, but note there is significant system-to-system scatter in this quantity) and that the kinematics show evidence for two components (due to in situ versus accreted), as observed. Our findings therefore add considerable weight to the importance of dissipative processes in the formation of stellar haloes and to the notion of a 'dual stellar halo'.
ABSTRACT
There is a severe tension between the observed star formation rate (SFR)–stellar mass (M⋆) relations reported by different authors at z = 1–4. In addition, the observations have not been ...successfully reproduced by state-of-the-art cosmological simulations that tend to predict a factor of 2–4 smaller SFRs at a fixed M⋆. We examine the evolution of the SFR–M⋆ relation of z = 1–4 galaxies using the skirt simulated spectral energy distributions of galaxies sampled from the Evolution and Assembly of GaLaxies and their Environments simulations. We derive SFRs and stellar masses by mimicking different observational techniques. We find that the tension between observed and simulated SFR–M⋆ relations is largely alleviated if similar methods are used to infer the galaxy properties. We find that relations relying on infrared wavelengths (e.g. 24 ${\rm \, \mu m}$, MIPS – 24, 70, and 160 ${\rm \, \mu m}$ or SPIRE – 250, 350, and 500 ${\rm \, \mu m}$) have SFRs that exceed the intrinsic relation by 0.5 dex. Relations that rely on the spectral energy distribution fitting technique underpredict the SFRs at a fixed stellar mass by −0.5 dex at z ∼ 4 but overpredict the measurements by 0.3 dex at z ∼ 1. Relations relying on dust-corrected rest-frame ultraviolet luminosities, are flatter since they overpredict/underpredict SFRs for low/high star-forming objects and yield deviations from the intrinsic relation from 0.10 to −0.13 dex at z ∼ 4. We suggest that the severe tension between different observational studies can be broadly explained by the fact that different groups employ different techniques to infer their SFRs.
We compare the results of various cosmological gas-dynamical codes used to simulate the formation of a galaxy in the Λ cold dark matter structure formation paradigm. The various runs (13 in total) ...differ in their numerical hydrodynamical treatment smoothed particle hydrodynamics (SPH), moving mesh and adaptive mesh refinement but share the same initial conditions and adopt in each case their latest published model of gas cooling, star formation and feedback. Despite the common halo assembly history, we find large code-to-code variations in the stellar mass, size, morphology and gas content of the galaxy at z= 0, due mainly to the different implementations of star formation and feedback. Compared with observation, most codes tend to produce an overly massive galaxy, smaller and less gas rich than typical spirals, with a massive bulge and a declining rotation curve. A stellar disc is discernible in most simulations, although its prominence varies widely from code to code. There is a well-defined trend between the effects of feedback and the severity of the disagreement with observed spirals. In general, models that are more effective at limiting the baryonic mass of the galaxy come closer to matching observed galaxy scaling laws, but often to the detriment of the disc component. Although numerical convergence is not particularly good for any of the codes, our conclusions hold at two different numerical resolutions. Some differences can also be traced to the different numerical techniques; for example, more gas seems able to cool and become available for star formation in grid-based codes than in SPH. However, this effect is small compared to the variations induced by different feedback prescriptions. We conclude that state-of-the-art simulations cannot yet uniquely predict the properties of the baryonic component of a galaxy, even when the assembly history of its host halo is fully specified. Developing feedback algorithms that can effectively regulate the mass of a galaxy without hindering the formation of high angular momentum stellar discs remains a challenge.
Abstract
We present a MUSE (Multi-Unit Spectroscopic Explorer) and KMOS (K-band Multi-Object Spectrograph) dynamical study 405 star-forming galaxies
at redshift z = 0.28–1.65 (median redshift
...$\bar{z}$
= 0.84). Our sample is representative of the star-forming ‘main sequence’, with star formation rates of SFR = 0.1–30 M⊙ yr−1 and stellar masses M
⋆ = 108–1011 M⊙. For 49 ± 4 per cent of our sample, the dynamics suggest rotational support, 24 ± 3 per cent are unresolved systems and 5 ± 2 per cent appear to be early-stage major mergers with components on 8–30 kpc scales. The remaining 22 ± 5 per cent appear to be dynamically complex, irregular (or face-on systems). For galaxies whose dynamics suggest rotational support, we derive inclination-corrected rotational velocities and show that these systems lie on a similar scaling between stellar mass and specific angular momentum as local spirals with j
⋆ = J/
$M_\star \propto M_\star ^{2/3}$
but with a redshift evolution that scales as j
⋆ ∝ M
$_\star ^{2/3}(1+z)^{-1}$
. We also identify a correlation between specific angular momentum and disc stability such that galaxies with the highest specific angular momentum (log(j
⋆/M
$_\star ^{2/3}$
) > 2.5) are the most stable, with Toomre Q = 1.10 ± 0.18, compared to Q = 0.53 ± 0.22 for galaxies with log(j
⋆/M
$_\star ^{2/3}$
) < 2.5. At a fixed mass, the Hubble Space Telescope morphologies of galaxies with the highest specific angular momentum resemble spiral galaxies, whilst those with low specific angular momentum are morphologically complex and dominated by several bright star-forming regions. This suggests that angular momentum plays a major role in defining the stability of gas discs: at z ∼ 1, massive galaxies that have discs with low specific angular momentum are globally unstable, clumpy and turbulent systems. In contrast, galaxies with high specific angular momentum have evolved into stable discs with spiral structure where star formation is a local (rather than global) process.
ABSTRACT
We present a study of the environment of 27 z = 3–4.5 bright quasars from the MUSE Analysis of Gas around Galaxies (MAGG) survey. With medium-depth Multi Unit Spectroscopic Explorer (MUSE) ...observations (4 h on target per field), we characterize the effects of quasars on their surroundings by studying simultaneously the properties of extended gas nebulae and Ly α emitters (LAEs) in the quasar host haloes. We detect extended (up to ≈100 kpc) Ly α emission around all MAGG quasars, finding a very weak redshift evolution between z = 3 and z = 6. By stacking the MUSE datacubes, we confidently detect extended emission of C iv and only marginally detect extended He ii up to ≈40 kpc, implying that the gas is metal enriched. Moreover, our observations show a significant overdensity of LAEs within 300 $\rm km~s^{-1}$ from the quasar systemic redshifts estimated from the nebular emission. The luminosity functions and equivalent width distributions of these LAEs show similar shapes with respect to LAEs away from quasars suggesting that the Ly α emission of the majority of these sources is not significantly boosted by the quasar radiation or other processes related to the quasar environment. Within this framework, the observed LAE overdensities and our kinematic measurements imply that bright quasars at z = 3–4.5 are hosted by haloes in the mass range $\approx 10^{12.0}\small{--}10^{12.5}~\rm M_\odot$.