BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Michelle Furlong

Michelle Furlong contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.GAastro-ph.CO

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
23works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2016arXiv

A chronicle of galaxy mass assembly in the EAGLE simulation

We analyse the mass assembly of central galaxies in the EAGLE hydrodynamical simulations. We build merger trees to connect galaxies to their progenitors at different redshifts and characterize their assembly histories by focusing on the time when half of the galaxy stellar mass was assembled into the main progenitor. We show that galaxies with stellar mass $M_*<10^{10.5}M_{\odot}$ assemble most of their stellar mass through star formation in the main progenitor (`in-situ&#39; star formation). This can be understood as a consequence of the steep rise in star formation efficiency with halo mass for these galaxies. For more massive galaxies, however, an increasing fraction of their stellar mass is formed outside the main progenitor and subsequently accreted. Consequently, while for low-mass galaxies the assembly time is close to the stellar formation time, the stars in high-mass galaxies typically formed long before half of the present-day stellar mass was assembled into a single object, giving rise to the observed anti-hierarchical downsizing trend. In a typical present-day $M_*\geq10^{11}M_{\odot}$ galaxy, around $20\%$ of the stellar mass has an external origin. This fraction decreases with increasing redshift. Bearing in mind that mergers only make an important contribution to the stellar mass growth of massive galaxies, we find that the dominant contribution comes from mergers with galaxies of mass greater than one tenth of the main progenitor&#39;s mass. The galaxy merger fraction derived from our simulations agrees with recent observational estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Galaxies in the EAGLE hydrodynamical simulation and in the Durham and Munich semi-analytical models

We compare global predictions from the EAGLE hydrodynamical simulation, and two semi-analytic (SA) models of galaxy formation, L-GALAXIES and GALFORM. All three models include the key physical processes for the formation and evolution of galaxies and their parameters are calibrated against a small number of observables at $z\approx 0$. The two SA models have been applied to merger trees constructed from the EAGLE dark matter only simulation. We find that at $z\leq 2$, both the galaxy stellar mass functions for stellar masses $M_*<10^{10.5} M_{\odot}$ and the median specific star formation rates (sSFRs) in the three models agree to better than $0.4$~dex. The evolution of the sSFR predicted by the three models closely follows the mass assembly history of dark matter haloes. In both EAGLE and L-GALAXIES there are more central passive galaxies with $M_*<10^{9.5} M_{\odot}$ than in L-GALAXIES. This difference is related to galaxies that have entered and then left a larger halo and which are treated as satellites in GALFORM. In the range $0<z<1$, the slope of the evolution of the star formation rate density in EAGLE is a factor of $\approx 1.5$ steeper than for the two SA models. The median sizes for galaxies with $M_*>10^{9.5} M_{\odot}$ differ in some instances by an order of magnitude, while the stellar mass-size relation in EAGLE is a factor of $\approx 2$ tighter than for the two SA models. Our results suggest the need for a revision of how SA models treat the effect of baryonic self-gravity on the underlying dark matter. The treatment of gas flows in the models needs to be revised based on detailed comparison with observations to understand in particular the evolution of the stellar mass-metallicity relation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Subhalo abundance matching and assembly bias in the EAGLE simulation

Subhalo abundance matching (SHAM) is a widely-used method to connect galaxies with dark matter structures in numerical simulations. SHAM predictions agree remarkably well with observations, yet they still lack strong theoretical support. We examine the performance, implementation, and assumptions of SHAM using the EAGLE project simulations. We find that $V_{\rm relax}$, the highest value of the circular velocity attained by a subhalo while it satisfies a relaxation criterion, is the subhalo property that correlates most strongly with galaxy stellar mass ($M_{\rm star}$). Using this parameter in SHAM, we retrieve the real-space clustering of EAGLE to within our statistical uncertainties on scales greater than $2$ Mpc for galaxies with $8.77<\log_{10}(M_{\rm star}[M_\odot])<10.77$. Conversely, clustering is overestimated by $30\%$ on scales below $2$ Mpc for galaxies with $8.77<\log_{10}(M_{\rm star}[M_\odot])<9.77$ because SHAM slightly overpredicts the fraction of satellites in massive haloes compared to EAGLE. The agreement is even better in redshift-space, where the clustering is recovered to within our statistical uncertainties for all masses and separations. Additionally, we analyse the dependence of galaxy clustering on properties other than halo mass, i.e. the assembly bias. We demonstrate assembly bias alters the clustering in EAGLE by $20\%$ and that $V_{\rm relax}$ captures its effect to within $15\%$. We trace small differences in the clustering to the failure of SHAM as typically implemented, i.e. the $M_{\rm star}$ assigned to a subhalo does not depend on i) its host halo mass, ii) whether it is a central or a satellite. In EAGLE we find that these assumptions are not completely satisfied.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The EAGLE simulations of galaxy formation: public release of halo and galaxy catalogues

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)

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The fundamental plane of star formation in galaxies revealed by the EAGLE hydrodynamical simulations

We investigate correlations between different physical properties of star-forming galaxies in the &#34;Evolution and Assembly of GaLaxies and their Environments&#34; (EAGLE) cosmological hydrodynamical simulation suite over the redshift range $0\le z\le 4.5$. A principal component analysis reveals that neutral gas fraction ($f_{\rm gas, neutral}$), stellar mass ($M_{\rm stellar}$) and star formation rate (SFR) account for most of the variance seen in the population, with galaxies tracing a two-dimensional, nearly flat, surface in the three-dimensional space of $f_{\rm gas, neutral}-M_{\rm stellar}-\rm SFR$ with little scatter. The location of this plane varies little with redshift, whereas galaxies themselves move along the plane as their $f_{\rm gas, neutral}$ and SFR drop with redshift. The positions of galaxies along the plane are highly correlated with gas metallicity. The metallicity can therefore be robustly predicted from $f_{\rm gas, neutral}$, or from the $M_{\rm stellar}$ and SFR. We argue that the appearance of this &#34;fundamental plane of star formation&#34; is a consequence of self-regulation, with the plane&#39;s curvature set by the dependence of the SFR on gas density and metallicity. We analyse a large compilation of observations spanning the redshift range $0\lesssim \rm z\lesssim 2.5$, and find that such a plane is also present in the data. The properties of the observed fundamental plane of star formation are in good agreement with EAGLE&#39;s predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Baryon effects on the internal structure of LCDM halos in the EAGLE simulations

We investigate the internal structure and density profiles of halos of mass $10^{10}-10^{14}~M_\odot$ in the Evolution and Assembly of Galaxies and their Environment (EAGLE) simulations. These follow the formation of galaxies in a $Λ$CDM Universe and include a treatment of the baryon physics thought to be relevant. The EAGLE simulations reproduce the observed present-day galaxy stellar mass function, as well as many other properties of the galaxy population as a function of time. We find significant differences between the masses of halos in the EAGLE simulations and in simulations that follow only the dark matter component. Nevertheless, halos are well described by the Navarro-Frenk-White (NFW) density profile at radii larger than ~5% of the virial radius but, closer to the centre, the presence of stars can produce cuspier profiles. Central enhancements in the total mass profile are most important in halos of mass $10^{12}-10^{13}M_\odot$, where the stellar fraction peaks. Over the radial range where they are well resolved, the resulting galaxy rotation curves are in very good agreement with observational data for galaxies with stellar mass $M_*<5\times10^{10}M_\odot$. We present an empirical fitting function that describes the total mass profiles and show that its parameters are strongly correlated with halo mass.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Colours and luminosities of z=0.1 simulated galaxies in the EAGLE simulations

We calculate the colours and luminosities of redshift z = 0.1 galaxies from the EAGLE simulation suite using the GALAXEV population synthesis models. We take into account obscuration by dust in birth clouds and diffuse ISM using a two-component screen model, following the prescription of Charlot and Fall. We compare models in which the dust optical depth is constant to models where it depends on gas metallicity, gas fraction and orientation. The colours of EAGLE galaxies for the more sophisticated models are in broad agreement with those of observed galaxies. In particular, EAGLE produces a red sequence of passive galaxies and a blue cloud of star forming galaxies, with approximately the correct fraction of galaxies in each population and with g-r colours within 0.1 magnitudes of those observed. Luminosity functions from UV to NIR wavelengths differ from observations at a level comparable to systematic shifts resulting from a choice between Petrosian and Kron photometric apertures. Despite the generally good agreement there are clear discrepancies with observations. The blue cloud of EAGLE galaxies extends to somewhat higher luminosities than in the data, consistent with the modest underestimate of the passive fraction in massive EAGLE galaxies. There is also a moderate excess of bright blue galaxies compared to observations. The overall level of agreement with the observed colour distribution suggests that EAGLE galaxies at z = 0.1 have ages, metallicities and levels of obscuration that are comparable to those of observed galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Intrinsic alignments of galaxies in the EAGLE and cosmo-OWLS simulations

We report results for the alignments of galaxies in the EAGLE and cosmo-OWLS simulations as a function of galaxy separation and halo mass. The combination of these hydro-cosmological simulations enables us to span four orders of magnitude in halo mass ($10.7<log_{10}(M_{200}/[h^{-1}M_\odot])<15$) and a large range of separations ($-1<log_{10}(r/[h^{-1}Mpc])< 2$). We focus on two classes of alignments: the orientations of galaxies with respect to either the directions to, or the orientations of, surrounding galaxies. We find that the strength of the alignment is a strongly decreasing function of the distance between galaxies. The orientation-direction alignment can remain significant up to ~100 Mpc, for galaxies hosted by the most massive haloes in our simulations. Galaxies hosted by more massive subhaloes show stronger alignment. At a fixed halo mass, more aspherical or prolate galaxies exhibit stronger alignments. The spatial distribution of satellites is anisotropic and significantly aligned with the major axis of the main host halo. The major axis of satellite galaxies, when all stars are considered, are preferentially aligned towards the centre of the main host halo. The predicted projected direction-orientation alignment, $ε_{g+}(r_{p})$, is in broad agreement with recent observations when only stars within the typical observable extent of a galaxy are used to define galaxy orientations. We find that the orientation-orientation alignment is weaker than the orientation-direction alignment on all scales. Overall, the strength of galaxy alignments depends strongly on the subset of stars that are used to measure the orientations of galaxies and it is always weaker than the alignment of the dark matter haloes. Thus, alignment models that use halo orientation as a direct proxy for galaxy orientation will overestimate the impact of intrinsic alignments on weak lensing analyses.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Molecular hydrogen abundances of galaxies in the EAGLE simulations

We investigate the abundance of galactic molecular hydrogen (H$_2$) in the &#34;Evolution and Assembly of GaLaxies and their Environments&#34; (EAGLE) cosmological hydrodynamic simulations. We assign H$_2$ masses to gas particles in the simulations in post-processing using two different prescriptions that depend on the local dust-to-gas ratio and the interstellar radiation field. Both result in H$_2$ galaxy mass functions that agree well with observations in the local and high-redshift Universe. The simulations reproduce the observed scaling relations between the mass of H$_2$ and the stellar mass, star formation rate and stellar surface density. Towards high edshifts, galaxies in the simulations display larger H$_2$ mass fractions, and correspondingly lower H$_2$ depletion timescales, also in good agreement with observations. The comoving mass density of H$_2$ in units of the critical density, $Ω_{\rm H_2}$, peaks at $z\approx 1.2-1.5$, later than the predicted peak of the cosmic star formation rate activity, at $z\approx 2$. This difference stems from the decrease in gas metallicity and increase in interstellar radiation field with redshift, both of which hamper H$_2$ formation. We find that the cosmic H$_2$ budget is dominated by galaxies with $M_{\rm H_2}>10^9\,\rm M_{\odot}$, star formation rates $>10\,\rm M_{\odot}\,\rm yr^{-1}$ and stellar masses $M_{\rm stellar}>10^{10}\,\rm M_{\odot}$, which are readily observable in the optical and near-IR. The match between the H$_2$ properties of galaxies that emerge in the simulations and observations is remarkable, particularly since H$_2$ observations were not used to adjust parameters in EAGLE.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Recycled stellar ejecta as fuel for star formation and implications for the origin of the galaxy mass-metallicity relation

We use cosmological, hydrodynamical simulations from the EAGLE and OWLS projects to assess the significance of recycled stellar ejecta as fuel for star formation. The fractional contributions of stellar mass loss to the cosmic star formation rate (SFR) and stellar mass densities increase with time, reaching $35 \%$ and $19 \%$, respectively, at $z=0$. The importance of recycling increases steeply with galaxy stellar mass for $M_{\ast} < 10^{10.5}$ M$_{\odot}$, and decreases mildly at higher mass. This trend arises from the mass dependence of feedback associated with star formation and AGN, which preferentially suppresses star formation fuelled by recycling. Recycling is more important for satellites than centrals and its contribution decreases with galactocentric radius. The relative contribution of AGB stars increases with time and towards galaxy centers. This is a consequence of the more gradual release of AGB ejecta compared to that of massive stars, and the preferential removal of the latter by star formation-driven outflows and by lock up in stellar remnants. Recycling-fuelled star formation exhibits a tight, positive correlation with galaxy metallicity, with a secondary dependence on the relative abundance of alpha elements (which are predominantly synthesized in massive stars), that is insensitive to the subgrid models for feedback. Hence, our conclusions are directly relevant for the origin of the mass-metallicity relation and metallicity gradients. Applying the relation between recycling and metallicity to the observed mass-metallicity relation yields our best estimate of the mass-dependent contribution of recycling. For centrals with a mass similar to that of the Milky Way, we infer the contributions of recycled stellar ejecta to the SFR and stellar mass to be $35 \%$ and $20 \%$, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The alignment and shape of dark matter, stellar, and hot gas distributions in the EAGLE and cosmo-OWLS simulations

We report the alignment and shape of dark matter, stellar, and hot gas distributions in the EAGLE and cosmo-OWLS simulations. The combination of these state-of-the-art hydro-cosmological simulations enables us to span four orders of magnitude in halo mass ($11 < log_{10}(M_{200}/ [h^{-1}M_\odot]) < 15$), a wide radial range ($-2.3 < log_{10}(r/[h^{-1}Mpc ]) < 1.3$) and redshifts $0 < z < 1$. The shape parameters of the dark matter, stellar and hot gas distributions follow qualitatively similar trends: they become more aspherical (and triaxial) with increasing halo mass, radius and redshift. We measure the misalignment of the baryonic components (hot gas and stars) of galaxies with their host halo as a function of halo mass, radius, redshift, and galaxy type (centrals vs satellites and early- vs late-type). Overall, galaxies align well with the local distribution of the total (mostly dark) matter. However, the stellar distributions on galactic scales exhibit a median misalignment of about 45-50 degrees with respect to their host haloes. This misalignment is reduced to 25-30 degrees in the most massive haloes ($13 < log_{10}(M_{200}/ [h^{-1}M_\odot ]) < 15$). Half of the disc galaxies in the EAGLE simulations have a misalignment angle with respect to their host haloes larger than 40 degrees. We present fitting functions and tabulated values for the probability distribution of galaxy-halo misalignment to enable a straightforward inclusion of our results into models of galaxy formations based on purely collisionless N-body simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The APOSTLE simulations: solutions to the Local Group&#39;s cosmic puzzles

The Local Group of galaxies offer some of the most discriminating tests of models of cosmic structure formation. For example, observations of the Milky Way (MW) and Andromeda satellite populations appear to be in disagreement with N-body simulations of the &#34;Lambda Cold Dark Matter&#34; (ΛCDM) model: there are far fewer satellite galaxies than substructures in cold dark matter halos (the &#34;missing satellites&#34; problem); dwarf galaxies seem to avoid the most massive substructures (the &#34;too-big-to-fail&#34; problem); and the brightest satellites appear to orbit their host galaxies on a thin plane (the &#34;planes of satellites&#34; problem). Here we present results from APOSTLE (A Project Of Simulating The Local Environment), a suite of cosmological hydrodynamic simulations of twelve volumes selected to match the kinematics of the Local Group (LG) members. Applying the Eagle code to the LG environment, we find that our simulations match the observed abundance of LG galaxies, including the satellite galaxies of the MW and Andromeda. Due to changes to the structure of halos and the evolution in the LG environment, the simulations reproduce the observed relation between stellar mass and velocity dispersion of individual dwarf spheroidal galaxies without necessitating the formation of cores in their dark matter profiles. Satellite systems form with a range of spatial anisotropies, including one similar to that of the MW, confirming that such a configuration is not unexpected in ΛCDM. Finally, based on the observed velocity dispersion, size, and stellar mass, we provide new estimates of the maximum circular velocity for the halos of nine MW dwarf spheroidals.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The distribution of atomic hydrogen in EAGLE galaxies: morphologies, profiles, and HI holes

We compare the mass and internal distribution of atomic hydrogen (HI) in 2200 present-day central galaxies with M_star > 10^10 M_Sun from the 100 Mpc EAGLE Reference simulation to observational data. Atomic hydrogen fractions are corrected for self-shielding using a fitting formula from radiative transfer simulations and for the presence of molecular hydrogen using an empirical or a theoretical prescription from the literature. The resulting neutral hydrogen fractions, M_(HI+H2)/M_star, agree with observations to better than 0.1 dex for galaxies with M_star between 10^10 and 10^11 M_Sun. Our fiducial, empirical H2 model based on gas pressure results in galactic HI mass fractions, M_HI/M_star, that agree with observations from the GASS survey to better than 0.3 dex, but the alternative theoretical H2 formula leads to a negative offset in M_HI/M_star of up to 0.5 dex. Visual inspection reveals that most HI disks in simulated HI-rich galaxies are vertically disturbed, plausibly due to recent accretion events. Many galaxies (up to 80 per cent) contain spuriously large HI holes, which are likely formed as a consequence of the feedback implementation in EAGLE. The HI mass-size relation of all simulated galaxies is close to (but 16 per cent steeper than) observed, and when only galaxies without large holes in the HI disc are considered, the agreement becomes excellent (better than 0.1 dex). The presence of large HI holes also makes the radial HI surface density profiles somewhat too low in the centre, at Σ_HI > 1 M_Sun pc^-2 (by a factor of <~ 2 compared to data from the Bluedisk survey). In the outer region (Σ_HI < 1 M_Sun pc^-2), the simulated profiles agree quantitatively with observations. Scaled by HI size, the simulated profiles of HI-rich (M_HI > 10^9.8 M_Sun) and control galaxies (10^9.1 M_Sun > M_HI > 10^9.8 M_Sun) follow each other closely, as observed. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The distribution of neutral hydrogen around high-redshift galaxies and quasars in the EAGLE simulation

The observed high covering fractions of neutral hydrogen (HI) with column densities above $\sim 10^{17} \rm{cm}^{-2}$ around Lyman-Break Galaxies (LBGs) and bright quasars at redshifts z ~ 2-3 has been identified as a challenge for simulations of galaxy formation. We use the EAGLE cosmological, hydrodynamical simulation, which has been shown to reproduce a wide range of galaxy properties and for which the subgrid feedback was calibrated without considering gas properties, to study the distribution of HI around high-redshift galaxies. We predict the covering fractions of strong HI absorbers ($N_{\rm{HI}} \gtrsim 10^{17} \rm{cm}^{-2}$) inside haloes to increase rapidly with redshift but to depend only weakly on halo mass. For massive ($M_{200} \gtrsim 10^{12} {\rm M_{\odot}}$) halos the covering fraction profiles are nearly scale-invariant and we provide fitting functions that reproduce the simulation results. While efficient feedback is required to increase the HI covering fractions to the high observed values, the distribution of strong absorbers in and around halos of a fixed mass is insensitive to factor of two variations in the strength of the stellar feedback. In contrast, at fixed stellar mass the predicted HI distribution is highly sensitive to the feedback efficiency. The fiducial EAGLE simulation reproduces both the observed global column density distribution function of HI and the observed radial covering fraction profiles of strong HI absorbers around LBGs and bright quasars.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The EAGLE simulations of galaxy formation: calibration of subgrid physics and model variations

We present results from thirteen cosmological simulations that explore the parameter space of the &#34;Evolution and Assembly of GaLaxies and their Environments&#34; (EAGLE) simulation project. Four of the simulations follow the evolution of a periodic cube L = 50 cMpc on a side, and each employs a different subgrid model of the energetic feedback associated with star formation. The relevant parameters were adjusted so that the simulations each reproduce the observed galaxy stellar mass function at z = 0.1. Three of the simulations fail to form disc galaxies as extended as observed, and we show analytically that this is a consequence of numerical radiative losses that reduce the efficiency of stellar feedback in high-density gas. Such losses are greatly reduced in the fourth simulation - the EAGLE reference model - by injecting more energy in higher density gas. This model produces galaxies with the observed size distribution, and also reproduces many galaxy scaling relations. In the remaining nine simulations, a single parameter or process of the reference model was varied at a time. We find that the properties of galaxies with stellar mass <~ M* (the &#34;knee&#34; of the galaxy stellar mass function) are largely governed by feedback associated with star formation, while those of more massive galaxies are also controlled by feedback from accretion onto their central black holes. Both processes must be efficient in order to reproduce the observed galaxy population. In general, simulations that have been calibrated to reproduce the low-redshift galaxy stellar mass function will still not form realistic galaxies, but the additional requirement that galaxy sizes be acceptable leads to agreement with a large range of observables.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The EAGLE simulations of galaxy formation: the importance of the hydrodynamics scheme

We present results from a subset of simulations from the &#34;Evolution and Assembly of GaLaxies and their Environments&#34; (EAGLE) suite in which the formulation of the hydrodynamics scheme is varied. We compare simulations that use the same subgrid models without re-calibration of the parameters but employing the standard GADGET flavour of smoothed particle hydrodynamics (SPH) instead of the more recent state-of-the-art ANARCHY formulation of SPH that was used in the fiducial EAGLE runs. We find that the properties of most galaxies, including their masses and sizes, are not significantly affected by the details of the hydrodynamics solver. However, the star formation rates of the most massive objects are affected by the lack of phase mixing due to spurious surface tension in the simulation using standard SPH. This affects the efficiency with which AGN activity can quench star formation in these galaxies and it also leads to differences in the intragroup medium that affect the X-ray emission from these objects. The differences that can be attributed to the hydrodynamics solver are, however, likely to be less important at lower resolution. We also find that the use of a time step limiter is important for achieving the feedback efficiency required to match observations of the low-mass end of the galaxy stellar mass function.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The effect of baryons on the inner density profiles of rich clusters

We use the &#34;Evolution and assembly of galaxies and their environments&#34; (EAGLE) cosmological simulation to investigate the effect of baryons on the density profiles of rich galaxy clusters. We focus on EAGLE clusters with $M_{200}>10^{14}~M_\odot$ of which we have six examples. The central brightest cluster galaxies (BCGs) in the simulation have steep stellar density profiles, $ρ_*(r) \propto r^{-3}$. Stars dominate the mass density for $r<10~\rm{kpc}$, and, as a result, the total mass density profiles are steeper than the Navarro-Frenk-White (NFW) profile, in remarkable agreement with observations. The dark matter halo itself closely follows the NFW form at all resolved radii ($r\gtrsim3.0~\rm{kpc}$). The EAGLE BCGs have similar surface brightness and line-of-sight velocity dispersion profiles as the BCGs in the sample of Newman et al., which have the most detailed measurements currently available. After subtracting the contribution of the stars to the central density, Newman et al. infer significantly shallower slopes than the NFW value, in contradiction with the EAGLE results. We discuss possible reasons for this discrepancy, and conclude that an inconsistency between the kinematical model adopted by Newman et al. for their BCGs, which assumes isotropic stellar orbits, and the kinematical structure of the EAGLE BCGs, in which the orbital stellar anisotropy varies with radius and tends to be radially biased, could explain at least part of the discrepancy.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The link between the assembly of the inner dark matter halo and the angular momentum evolution of galaxies in the EAGLE simulation

We explore the co-evolution of the specific angular momentum of dark matter haloes and the cold baryons that comprise the galaxies within. We study over two thousand central galaxies within the reference cosmological hydrodynamical simulation of the &#34;Evolution and Assembly of GaLaxies and their Environments&#34; (EAGLE) project. We employ a methodology within which the evolutionary history of a system is specified by the time-evolving properties of the Lagrangian particles that define it at z=0. We find a strong correlation between the evolution of the specific angular momentum of today&#39;s stars (cold gas) and that of the inner (whole) dark matter halo they are associated with. This link is particularly strong for the stars formed before the epoch of maximum expansion and subsequent collapse of the central dark matter halo (turnaround). Spheroids are typically assembled primarily from stars formed prior to turnaround, and are therefore destined to suffer a net loss of angular momentum associated with the strong merging activity during the assembly of the inner dark matter halo. Stellar discs retain their specific angular momentum since they are comprised of stars formed mainly after turnaround, from gas that mostly preserves the high specific angular momentum it acquired by tidal torques during the linear growth of the halo. Since the specific angular momentum loss of the stars is tied to the galaxy&#39;s morphology today, it may be possible to use our results to predict, statistically, the assembly history of a halo given the morphology of the galaxy it hosts.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The unexpected diversity of dwarf galaxy rotation curves

We examine the circular velocity profiles of galaxies in ΛCDM cosmological hydrodynamical simulations from the EAGLE and LOCAL GROUPS projects and compare them with a compilation of observed rotation curves of galaxies spanning a wide range in mass. The shape of the circular velocity profiles of simulated galaxies varies systematically as a function of galaxy mass, but shows remarkably little variation at fixed maximum circular velocity. This is especially true for low-mass dark matter-dominated systems, reflecting the expected similarity of the underlying cold dark matter haloes. This is at odds with observed dwarf galaxies, which show a large diversity of rotation curve shapes, even at fixed maximum rotation speed. Some dwarfs have rotation curves that agree well with simulations, others do not. The latter are systems where the inferred mass enclosed in the inner regions is much lower than expected for cold dark matter haloes and include many galaxies where previous work claims the presence of a constant density &#34;core&#34;. The &#34;cusp vs core&#34; issue is thus better characterized as an &#34;inner mass deficit&#34; problem than as a density slope mismatch. For several galaxies the magnitude of this inner mass deficit is well in excess of that reported in recent simulations where cores result from baryon-induced fluctuations in the gravitational potential. We conclude that one or more of the following statements must be true: (i) the dark matter is more complex than envisaged by any current model; (ii) current simulations fail to reproduce the effects of baryons on the inner regions of dwarf galaxies; and/or (iii) the mass profiles of &#34;inner mass deficit&#34; galaxies inferred from kinematic data are incorrect.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Bent by baryons: the low mass galaxy-halo relation

The relation between galaxies and dark matter halos is of vital importance for evaluating theoretical predictions of structure formation and galaxy formation physics. We show that the widely used method of abundance matching based on dark matter only simulations fails at the low mass end because two of its underlying assumptions are broken: only a small fraction of low mass (below 10^9.5 solar masses) halos host a visible galaxy, and halos grow at a lower rate due to the effect of baryons. In this regime, reliance on dark matter only simulations for abundance matching is neither accurate nor self-consistent. We find that the reported discrepancy between observational estimates of the halo masses of dwarf galaxies and the values predicted by abundance matching does not point to a failure of LCDM, but simply to a failure to account for baryonic effects. Our results also imply that the Local Group contains only a few hundred observable galaxies in contrast with the thousands of faint dwarfs that abundance matching would suggest. We show how relations derived from abundance matching can be corrected, so that they can be used self-consistently to calibrate models of galaxy formation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Local Group galaxies emerge from the dark

The &#34;Lambda Cold Dark Matter&#34; (LCDM) model of cosmic structure formation is eminently falsifiable: once its parameters are fixed on large scales, it becomes testable in the nearby Universe. Observations within our Local Group of galaxies, including the satellite populations of the Milky Way and Andromeda, appear to contradict LCDM predictions: there are far fewer satellite galaxies than dark matter halos (the &#34;missing satellites&#34; problem), galaxies seem to avoid the largest substructures (the &#34;too big to fail&#34; problem), and the brightest satellites appear to orbit their host galaxies on a thin plane (the &#34;planes of satellites&#34; problem). We present results from the first hydrodynamic simulations of the Local Group that match the observed abundance of galaxies. We find that when baryonic and dark matter are followed simultaneously in the context of a realistic galaxy formation model, all three &#34;problems&#34; are resolved within the LCDM paradigm.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The chosen few: the low mass halos that host faint galaxies

Since reionization prevents star formation in most halos below 3 x 10^9 solar masses, dwarf galaxies only populate a fraction of existing dark matter halos. We use hydrodynamic cosmological simulations of the Local Group to study the discriminating factors for galaxy formation in the early Universe and connect them to the present-day properties of galaxies and halos. A combination of selection effects related to reionization, and the subsequent evolution of halos in different environments, introduces strong biases between the population of halos that host dwarf galaxies, and the total halo population. Halos that host galaxies formed earlier and are more concentrated. In addition, halos more affected by tidal stripping are more likely to host a galaxy for a given mass or maximum circular velocity, vmax, today. Consequently, satellite halos are populated more frequently than field halos, and satellite halos of 10^8 - 10^9 solar masses or vmax of 12 - 20 km/s, similar to the Local Group dwarf spheroidals, have experienced a greater than average reduction in both mass and vmax after infall. They are on closer, more radial orbits with higher infall velocities and earlier infall times. Together, these effects make dwarf galaxies highly biased tracers of the underlying dark matter distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The EAGLE project: Simulating the evolution and assembly of galaxies and their environments

We introduce the Virgo Consortium&#39;s EAGLE project, a suite of hydrodynamical simulations that follow the formation of galaxies and black holes in representative volumes. We discuss the limitations of such simulations in light of their finite resolution and poorly constrained subgrid physics, and how these affect their predictive power. One major improvement is our treatment of feedback from massive stars and AGN in which thermal energy is injected into the gas without the need to turn off cooling or hydrodynamical forces, allowing winds to develop without predetermined speed or mass loading factors. Because the feedback efficiencies cannot be predicted from first principles, we calibrate them to the z~0 galaxy stellar mass function and the amplitude of the galaxy-central black hole mass relation, also taking galaxy sizes into account. The observed galaxy mass function is reproduced to $\lesssim 0.2$ dex over the full mass range, $10^8 < M_*/M_\odot \lesssim 10^{11}$, a level of agreement close to that attained by semi-analytic models, and unprecedented for hydrodynamical simulations. We compare our results to a representative set of low-redshift observables not considered in the calibration, and find good agreement with the observed galaxy specific star formation rates, passive fractions, Tully-Fisher relation, total stellar luminosities of galaxy clusters, and column density distributions of intergalactic CIV and OVI. While the mass-metallicity relations for gas and stars are consistent with observations for $M_* \gtrsim 10^9 M_\odot$, they are insufficiently steep at lower masses. The gas fractions and temperatures are too high for clusters of galaxies, but for groups these discrepancies can be resolved by adopting a higher heating temperature in the subgrid prescription for AGN feedback. EAGLE constitutes a valuable new resource for studies of galaxy formation.

0 citations0 reviewsNo code linkedOpen access