Search

We have developed a new tool to analyse galaxies in the EAGLE simulations as close as possible to observations.We investigated the evolution of their kinematic properties by means of the angular momentum proxy parameter,$ \lambda_{Re} $for galaxies with $M_{*} \ge 5 \times 10^{9} M_{\odot}$ in the RefL0100N1504 simulation up to redshift two (z = 2). Galaxies in the simulation show a wide variety of kinematic features, similiar to those found in integral-field spectroscopic studies. At z=0 the distribution of galaxies in the ${\lambda}_{Re}-{\epsilon}$ plane is also in good agreement with results from observations. Scaling relations at z = 0 indicate that there is critical mass, $M_{crit} = 10^{10.3} M_{\odot}$, that divides two different regimes when we include the ${\lambda}_{Re}$ parameter. The simulation shows that the distribution of galaxies in the ${\lambda}_{Re}-{\epsilon}$ plane evolves with time until z = 2 when galaxies are equally distributed both in ${\lambda}_{Re}$ and ${\epsilon}$. We studied the evolution of ${\lambda}_{Re}$ with time and found that there is no connection between the angular momentum at z = 2 and z = 0. All systems reach their maximum ${\lambda}_{Re}$ at z = 1 and then steadily lose angular momentum regardless of their merger history, except for the high star-forming systems that sustain that maximum value over time. The evolution of the Re in galaxies that have not experienced any merger in the last 10 Gyr can be explained by their level of gas accretion., Comment: 16 aojes, 9 Figures, accepted MNRAS

Received 12th May 2021. Accepted 17th May 2021. Published online 28th May 2021. Issue published 1st August 2021.

We test cosmological hydrodynamical simulations of galaxy formation regarding the properties of the Blue Cloud (BC), Green Valley (GV) and Red Sequence (RS), as measured on the 4000$\small{ \mathring {\mathrm A}}$ break strength vs stellar mass plane at $z=0.1$. We analyse the RefL0100N1504 run of EAGLE and the TNG100 run of IllustrisTNG project, by comparing them with the Sloan Digital Sky Survey, while taking into account selection bias. Our analysis focuses on the GV, within stellar mass $\log\,\mathrm{M_\star/M_{\odot}} \simeq 10-11$, selected from the bimodal distribution of galaxies on the D$_n$(4000) vs stellar mass plane, following Angthopo et al. methodology. Both simulations match the fraction of AGN in the green-valley. However, they over-produce quiescent GV galaxies with respect to observations, with IllustrisTNG yielding a higher fraction of quiescent GV galaxies than EAGLE. In both, GV galaxies have older luminosity-weighted ages with respect to the SDSS, while a better match is found for mass-weighted ages. We find EAGLE GV galaxies quench their star formation early, but undergo later episodes of star formation, matching observations. In contrast, IllustrisTNG GV galaxies have a more extended SFH, and quench more effectively at later cosmic times, producing the excess of quenched galaxies in GV compared with SDSS, based on the 4000$\small{ \mathring {\mathrm A}}$ break strength. These results suggest the AGN feedback subgrid physics, more specifically, the threshold halo mass for black hole input and the black hole seed mass, could be the primary cause of the over-production of quiescent galaxies found with respect to the observational constraints., 20 pages, 8 + 1 figures. Published in MNRAS. Slight amendments from original version (Fig. 3, 4, 7 and 8)

We report the discovery of BOSS-EUVLG1 at z = 2.469, by far the most luminous, almost un-obscured star-forming galaxy known at any redshift. First classified as a QSO within the Baryon Oscillation Spectroscopic Survey, follow-up observations with the Gran Telescopio Canarias reveal that its large luminosity, MUV ≃ -24.40 and log(LLyα/erg s-1) ≃ 44.0, is due to an intense burst of star formation, and not to an active galactic nucleus or gravitational lensing. BOSS-EUVLG1 is a compact (reff ≃ 1.2 kpc), young (4-5 Myr) starburst with a stellar mass log(Mz.ast;/M) = 10.0-0.1 and a prodigious star formation rate of ≃1000 M yr-1. However, it is metal- and dust-poor [12 + log(O/H) = 8.13 ± 0.19, E(B - V) ≃ 0.07, log(LIR/LUV) < -1.2], indicating that we are witnessing the very early phase of an intense starburst that has had no time to enrich the ISM. BOSS-EUVLG1 might represent a short-lived (, With funding from the Spanish government through the "María de Maeztu Unit of Excellence" accreditation (MDM-2017-0737)

We present the results of a study investigating the dust attenuation law at $z\simeq 5$, based on synthetic spectral energy distributions (SEDs) calculated for a sample of N=498 galaxies drawn from the First Billion Years (FiBY) simulation project. The simulated galaxies at $z\simeq 5$, which have M$_{1500} \leq -18.0$ and $7.5 \leq \rm{log(M/M}_{\odot}\rm{)} \leq 10.2$, display a mass-dependent $\alpha$-enhancement, with a median value of $[\alpha/\rm{Fe}]_{z=5}~\simeq~4~\times~[\alpha/\rm{Fe}]_{Z_{\odot}}$. The median Fe/H ratio of the simulated galaxies is $0.14\pm0.05$ which, even including the effects of nebular continuum, produces steep intrinsic UV continuum slopes; $\langle \beta_{i} \rangle = -2.4 \pm 0.05$. Using a set of simple dust attenuation models, in which the wavelength-dependent attenuation is assumed to be of the form $A(\lambda) \propto \lambda^{n}$, we explore the parameter values which best reproduce the observed $z=5$ luminosity function (LF) and colour-magnitude relation (CMR). We find that a simple model in which the absolute UV attenuation is a linearly increasing function of log stellar mass, and the dust attenuation slope ($n$) is within the range $-0.7 \leq n \leq-0.3$, can successfully reproduce the LF and CMR over a wide range of stellar population synthesis model (SPS) assumptions. This range of attenuation curves is consistent with a power-law fit to the Calzetti attenuation law in the UV ($n=-0.55$), and other similarly `grey' star-forming galaxy attenuation curves recently derived at $z\simeq2$. In contrast, attenuation curves as steep as the Small Magellanic Cloud (SMC) extinction curve ($n=-1.24$) are formally ruled out. Finally, we show that our models are consistent with recent 1.3mm ALMA observations of the Hubble Ultra Deep Field (HUDF), and predict the form of the $z\simeq5$ IRX$-\beta$ relation., Comment: 23 pages, 19 figures, accepted for publication in MNRAS

We examine the properties of barred disc galaxies in a LCDM cosmological hydrodynamical simulation from the EAGLE project. Our study follows the formation of 269 discs identified at z = 0 in the stellar mass range 10.6 < log Mstr /M < 11. These discs show a wide range of bar strengths, from unbarred discs to weak bars to strongly barred systems (= 20%). Bars in these systems develop after redshift = 1.3, on timescales that depend sen- sitively on the strength of the pattern. Strong bars develop relatively quickly (in a few Gyr, = 10 disc rotation periods) in systems that are disc dominated, gas poor, and have declining rotation curves. Weak bars develop more slowly in systems where the disc is less gravitation- ally important, and are still growing at z = 0. Unbarred galaxies are comparatively gas-rich discs whose rotation speeds do not exceed the maximum circular velocity of the halos they inhabit. Bar lengths compare favourably with observations, ranging from 0.2 to 0.8 times the radius containing 90% of the stars. Bars slow down remarkably quickly as they grow, causing the inner regions of the surrounding dark halo to expand. At z = 0 strong bars have corota- tion radii roughly ten times the bar length. Such slow bars are inconsistent with the few cases where pattern speeds have been measured or inferred observationally, a discrepancy that, if confirmed, might prove a challenge for disc galaxy formation in LCDM., 11 pages, 10 figures

Dwarfs are the largest population of galaxies in number in the nearby Universe. Deep spectroscopic data are still missing to obtain a better understanding of their formation and evolution processes. This study shows the results obtained from a spectroscopic campaign in the Perseus cluster. We have obtained 963 new galaxy spectra. We have measured the recessional velocity of the galaxies by using a cross-correlation technique. These data have been used to obtain the cluster membership, the dynamics of the galaxies, and the spectroscopic luminosity function (LF) of the cluster. The cluster membership was obtained by using the peak + gap technique, reporting a total of 403 galaxies as cluster members within 1.4r200. The mean velocity and velocity dispersion of the cluster galaxies are Vc = 5258 km s−1 and σc = 1040 km s−1, respectively. We obtained M200 = 1.2 × 1015 M⊙ and r200 = 2.2 Mpc for this cluster. The clusters members were classified blue and red according to their g − r stellar colour. The velocity dispersion of these two families of galaxies is different, indicating that the blue galaxies can be classified as recently accreted into the cluster. We present the spectroscopic galaxy LF of the cluster. This function turned to be flat: α = 0.99 ± 0.06. In addition, blue and red galaxies show similar densities in the faint end of the LF. This indicates that Perseus does not have a population of red dwarf galaxias as large as other nearby clusters. We have compared the LF of the Perseus cluster with other spectroscopic LFs of nearby clusters and those from cosmological simulations. This comparison shows that the spectroscopic LF of nearby galaxy cluster is far from universal.

We present the public data release of halo and galaxy catalogues extracted from the EAGLE suite of cosmological hydrodynamical simulations of galaxy formation. These simulations were performed with an enhanced version of the GADGET code that includes a modified hydrodynamics solver, time-step limiter and subgrid treatments of baryonic physics, such as stellar mass loss, element-by-element radiative cooling, star formation and feedback from star formation and black hole accretion. The simulation suite includes runs performed in volumes ranging from 25 to 100 comoving megaparsecs per side, with numerical resolution chosen to marginally resolve the Jeans mass of the gas at the star formation threshold. The free parameters of the subgrid models for feedback are calibrated to the redshift z=0 galaxy stellar mass function, galaxy sizes and black hole mass - stellar mass relation. The simulations have been shown to match a wide range of observations for present-day and higher-redshift galaxies. The raw particle data have been used to link galaxies across redshifts by creating merger trees. The indexing of the tree produces a simple way to connect a galaxy at one redshift to its progenitors at higher redshift and to identify its descendants at lower redshift. In this paper we present a relational database which we are making available for general use. A large number of properties of haloes and galaxies and their merger trees are stored in the database, including stellar masses, star formation rates, metallicities, photometric measurements and mock gri images. Complex queries can be created to explore the evolution of more than 10^5 galaxies, examples of which are provided in appendix. (abridged), Comment: 12 pages, 3 figures + appendices. The simulation data is available at http://www.eaglesim.org/database.php (comments welcome). Paper accepted to AstronomyComputing. Note there is a typo in the published version, Section 4.1, replace 3/2 with 2/3 for velocity dispersion

We present the results of a new study of dust attenuation at redshifts $3 < z < 4$ based on a sample of $236$ star-forming galaxies from the VANDELS spectroscopic survey. Motivated by results from the First Billion Years (FiBY) simulation project, we argue that the intrinsic spectral energy distributions (SEDs) of star-forming galaxies at these redshifts have a self-similar shape across the mass range $8.2 \leq$ log$(M_{\star}/M_{\odot}) \leq 10.6$ probed by our sample. Using FiBY data, we construct a set of intrinsic SED templates which incorporate both detailed star formation and chemical abundance histories, and a variety of stellar population synthesis (SPS) model assumptions. With this set of intrinsic SEDs, we present a novel approach for directly recovering the shape and normalization of the dust attenuation curve. We find, across all of the intrinsic templates considered, that the average attenuation curve for star-forming galaxies at $z\simeq3.5$ is similar in shape to the commonly-adopted Calzetti starburst law, with an average total-to-selective attenuation ratio of $R_{V}=4.18\pm0.29$. We show that the optical attenuation ($A_V$) versus stellar mass ($M_{\star}$) relation predicted using our method is consistent with recent ALMA observations of galaxies at $2, Comment: 16 pages, 12 figures, accepted for publication in MNRAS

We present a gravitational lensing and X-ray analysis of a massive galaxy cluster and its surroundings. The core of MACS\,J0717.5+3745 ($M(R, 18 pages, 6 figures, 8 tables, Accepted for publication in MNRAS

For the same stellar mass, physically smaller star-forming galaxies are also metal richer (Ellison et al. 2008). What causes the relation remains unclear. The central star-forming galaxies in the EAGLE cosmological numerical simulation reproduce the observed trend. We use them to explore the origin of the relation assuming that the physical mechanism responsible for the anti-correlation between size and gas-phase metallicity is the same in the simulated and the observed galaxies. We consider the three most likely causes: (1) metal-poor gas inflows feeding the star-formation process, (2) metal-rich gas outflows particularly efficient in shallow gravitational potentials, and (3) enhanced efficiency of the star-formation process in compact galaxies. Outflows (2) and enhanced star-formation efficiency (3) can be discarded. Metal-poor gas inflows (1) cause the correlation in the simulated galaxies. Galaxies grow in size with time, so those that receive gas later are both metal poorer and larger, giving rise to the observed anti-correlation. As expected within this explanation, larger galaxies have younger stellar populations. We explore the variation with redshift of the relation, which is maintained up to, at least, redshift 8., Comment: Accepted for publication in ApJ. 23 pages and 18 figures

Galaxies in clusters are strongly affected by their environment. They evolve according to several physical mechanisms that are active in clusters. Their efficiency can strongly depend on the orbital configuration of the galaxies. Our aim is to analyse the orbits of the galaxies in the cluster Abell 85, based on the study of the galaxy velocity anisotropy parameter. We have solved the Jeans equation under the assumption that the galaxies in A85 are collisionless objects, within the spherically symmetric gravitational potential of the virialized cluster. The mass of the cluster was estimated with X-ray and caustic analyses. We find that the anisotropy profile of the full galaxy population in A85 is an increasing monotonic function of the distance from the cluster centre: on average, galaxies in the central region (r/r200 < 0.3) are on isotropic orbits, while galaxies in the outer regions are on radial orbits. We also find that the orbital properties of the galaxies strongly depend on their stellar colour. In particular, blue galaxies are on less radial orbits than red galaxies. The different families of cluster galaxies considered here have the pseudo phase-space density profiles Q(r) and Qr(r) consistent with the profiles expected in virialized dark matter halos in $N$-body simulations. This result suggests that the galaxies in A85 have reached dynamical equilibrium within the cluster potential. Our results indicate that the origin of the blue and red colour of the different galaxy populations is the different orbital shape rather than the accretion time., 15 pages, 15 figures. Accepted for publication at MNRAS

Long gamma-ray burst's (LGRB's) association to the death of massive stars suggest they could be used to probe the cosmic star formation history (CSFH) with high accuracy, due to their high luminosities. We utilise cosmological simulations from the First Billion Years project to investigate the biases between the CSFH and the LGRB rate at z>5, assuming various different models and constraints on the progenitors of LGRBs. We populate LGRBs using a selection based on environmental properties and demonstrate that the LGRB rate should trace the CSFH to high redshifts. The measured LGRB rate suggests that LGRBs have opening angles of theta_jet=0.1 deg, although the degeneracy with the progenitor model cannot rule out an underlying bias. We demonstrate that proxies that relate the LGRB rate with global LGRB host properties do not reflect the underlying LGRB environment, and are in fact a result of the host galaxy's spatial properties, such that LGRBs can exist in galaxies of solar metallicity. However, we find a class of host galaxies that have low stellar mass and are metal-rich, and that their metallicity dispersions would not allow low-metallicity environments. Detection of hosts with this set of properties would directly reflect the progenitor's environment. We predict that 10% of LGRBs per year are associated with this set of galaxies that would have forbidden line emission that could be detected by instruments on the James Webb Space Telescope. Such a discovery would place strong constraints on the collapsar model and suggest other avenues to be investigated., Comment: 13 pages, 8 figures, 1 table, accepted for publication in MNRAS

We compare the predicted conditional mass function (CMF) of dark matter halos from two theoretical prescriptions against numerical N-body simulations, both in overdense and underdense regions and at different Eulerian scales ranging from $5$ to $30\,h^{-1}\,$Mpc. In particular, we consider in detail a locally-implemented rescaling of the unconditional mass function (UMF) already discussed in the literature, and also a generalization of the standard rescaling method described in the extended Press-Schechter formalism. First, we test the consistency of these two rescalings by verifying the normalization of the CMF at different scales, and showing that none of the proposed cases provides a normalized CMF. In order to satisfy the normalization condition, we include a modification in the rescaling procedure. After this modification, the resulting CMF generally provides a better description of numerical results. We finally present an analytical fit to the ratio between the CMF and the UMF (also known as the matter-to-halo bias function) in underdense regions, which could be of special interest to speed-up the computation of the halo abundance when studying void statistics. In this case, the CMF prescription based on the locally-implemented rescaling provides a slightly better description of the numerical results when compared to the standard rescaling., Comment: 20 pages, 13 figures. Accepted for publication in MNRAS

We compare the results of thirteen cosmological gasdynamical codes used to simulate the formation of a galaxy in the LCDM structure formation paradigm. The various runs differ in their hydrodynamical treatment (SPH, moving-mesh and AMR) but share the same initial conditions and adopt their latest published model of 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 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 disk is discernible in most simulations, though 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 observation. 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 disk component. Our conclusions hold at two different numerical resolutions. Some differences can also be traced to the numerical techniques: 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 disks remains a challenge.

We explore the dynamical signatures imprinted by baryons on dark matter haloes during the formation process using the OverWhelmingly Large Simulations (OWLS), a set of state-of-the-art high-resolution cosmological hydrodynamical simulations. We present a detailed study of the effects of the implemented feedback prescriptions on the orbits of dark matter particles, stellar particles and subhaloes, analysing runs with no feedback, with stellar feedback and with feedback from supermassive black holes. We focus on the central regions (0.25 r_{200}) of haloes with virial masses ~ 6 x 10^{13} (~ 7 x 10^{11}) Msun/h at z = 0(2). We also investigate how the orbital content (relative fractions of the different orbital types) of these haloes depends on several key parameters such as their mass, redshift and dynamical state. The results of spectral analyses of the orbital content of these simulations are compared, and the change in fraction of box, tube and irregular orbits is quantified. Box orbits are found to dominate the orbital structure of dark matter haloes in cosmological simulations. There is a strong anticorrelation between the fraction of box orbits and the central baryon fraction. While radiative cooling acts to reduce the fraction of box orbits, strong feedback implementations result in a similar orbital distribution to that of the dark matter only case. The orbital content described by the stellar particles is found to be remarkably similar to that drawn from the orbits of dark matter particles, suggesting that either they have forgotten their dynamical history, or that subhaloes bringing in stars are not biased significantly with respect to the main distribution. The orbital content of the subhaloes is in broad agreement with that seen in the outer regions of the particle distributions.

We present a new deep spectroscopic catalogue for Abell 85, within 3.0 $\times$ 2.6 Mpc$^2$ and down to $M_{r} \sim M_{r}^* +6$. Using the Visible Multi-Object Spectrograph at the Very Large Telescope (VIMOS@VLT) and the AutoFiber 2 at the William Herschel Telescope (AF2@WHT), we obtained almost 1,430 new redshifts for galaxies with $m_r \leq 21$ mag and $\langle \mu_{e,r} \rangle \leq 24$ mag arcsec$^{-2}$. These redshifts, together with SDSS-DR6 and NED spectroscopic information, result in 460 confirmed cluster members. This dataset allows the study of the luminosity function (LF) of the cluster galaxies covering three orders of magnitudes in luminosities. The total and radial LFs are best modelled by a double Schechter function. The normalized LFs show that their bright ($M_{r} \leq -21.5$) and faint ($M_{r}\geq -18.0$) ends are independent of clustercentric distance and similar to the field LFs unlike the intermediate luminosity range ($-21.5 \leq M_{r} \leq -18.0$). Similar results are found for the LFs of the dominant types of galaxies: red, passive, virialized and early-infall members. On the contrary, the LFs of blue, star forming, non-virialized and recent-infall galaxies are well described by a single Schechter function. These populations contribute to a small fraction of the galaxy density in the innermost cluster region. However, in the outskirts of the cluster, they have similar densities to red, passive, virialized and early-infall members at the LF faint end. These results confirm a clear dependence of the colour and star formation of Abell 85 members in the cluster centric distance., Comment: 15 pages, 12 figures, accepted for publication on MNRAS

The intergalactic magnetic field within filaments should be less polluted by magnetised outflows from active galaxies than magnetic fields in clusters. Therefore, filaments may be a better laboratory to study magnetic field amplification by structure formation than galaxy clusters which typically host many more active galaxies. We present highly resolved cosmological AMR simulations of magnetic fields in the cosmos and make predictions about the evolution and structure of magnetic fields in filaments. Comparing our results to observational evidence for magnetic fields in filaments suggests that amplification of seed fields by gravitational collapse is not sufficient to produce IGM fields. Finally, implications for cosmic ray transport are discussed., Comment: accepted by ApJL

Feedback from energy liberated by gas accretion onto black holes (BHs) is an attractive mechanism to explain the exponential cut-off at the massive end of the galaxy stellar mass function (SMF). Semi-analytic models of galaxy formation in which this form of feedback is assumed to suppress cooling in haloes where the gas cooling time is large compared to the dynamical time do indeed achieve a good match to the observed SMF. Furthermore, hydrodynamic simulations of individual halos in which gas is assumed to accrete onto the central BH at the Bondi rate have shown that a self-regulating regime is established in which the BH grows just enough to liberate an amount of energy comparable to the thermal energy of the halo. However, this process is efficient at suppressing the growth not only of massive galaxies but also of galaxies like the Milky Way, leading to disagreement with the observed SMF. The Bondi accretion rate, however, is inappropriate when the accreting material has angular momentum. We present an improved accretion model that takes into account the circularisation and subsequent viscous transport of infalling material and include it as a "subgrid" model in hydrodynamic simulations of the evolution of halos with a wide range of masses. The resulting accretion rates are generally low in low mass ($\lsim 10^{11.5} \msun$) halos, but show outbursts of Eddington-limited accretion during galaxy mergers. During outbursts these objects strongly resemble quasars. In higher mass haloes, gas accretion occurs continuously, typically at $~10$ % of the Eddington rate, which is conducive to the formation of radio jets. The resulting dependence of the accretion behaviour on halo mass induces a break in the relation between galaxy stellar mass and halo mass in these simulations that matches observations., 24 pages, 14 figures. Accepted to MNRAS for publication. Updated version with typos corrected

Local extremely metal-poor (XMP) galaxies are of particular astrophysical interest since they allow us to look into physical processes characteristic of the early Universe, from the assembly of galaxy disks to the formation of stars in conditions of low metallicity. Given the luminosity-metallicity relationship, all galaxies fainter than Mr < -13 are expected to be XMPs. Therefore, XMPs should be common in galaxy surveys. However, they are not, because several observational biases hamper their detection. This work compares the number of faint XMPs in the SDSS-DR7 spectroscopic survey with the expected number, given the known biases and the observed galaxy luminosity function. The faint end of the luminosity function is poorly constrained observationally, but it determines the expected number of XMPs. Surprisingly, the number of observed faint XMPs (around 10) is over-predicted by our calculation, unless the upturn in the faint end of the luminosity function is not present in the model. The lack of an upturn can be naturally understood if most XMPs are central galaxies in their low-mass dark matter halos, which are highly depleted in baryons due to interaction with the cosmic ultraviolet background and to other physical processes. Our result also suggests that the upturn towards low luminosity of the observed galaxy luminosity function is due to satellite galaxies., Accepted for publication in ApJ

We use numerical simulations to investigate how the statistical properties of dark matter (DM) haloes are affected by the baryonic processes associated with galaxy formation. We focus on how these processes influence the spin and shape of a large number of DM haloes covering a wide range of mass scales, from galaxies to clusters at redshifts zero and one, extending to dwarf galaxies at redshift two. The haloes are extracted from the OverWhelmingly Large Simulations, a suite of state-of-the-art high-resolution cosmological simulations run with a range of feedback prescriptions. We find that the median spin parameter in DM-only simulations is independent of mass, redshift and cosmology. At z=0 baryons increase the spin of the DM in the central region (, 15 pages, 8 figures; accepted for publication in MNRAS

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 $\Lambda$CDM 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 simulation also tracks 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.4 dex lower than observed, but the evolution of the rates track the observations closely. The unprecedented level of agreement between simulation and data makes EAGLE a powerful resource to understand the physical processes that govern galaxy formation., Comment: Accepted for publication by MNRAS; 18 pages, 7 figures; comments welcome

We compare the results of thirteen cosmological gasdynamical codes used to simulate the formation of a galaxy in the LCDM structure formation paradigm. The various runs differ in their hydrodynamical treatment (SPH, moving-mesh and AMR) but share the same initial conditions and adopt their latest published model of 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 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 disk is discernible in most simulations, though 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 observation. 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 disk component. Our conclusions hold at two different numerical resolutions. Some differences can also be traced to the numerical techniques: 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 disks remains a challenge., 28 pages, 16 figures, MNRAS accepted. High resolution figures and additional images can be found at http://www.aip.de/People/cscannapieco/aquila/