The Cluster-EAGLE project: Velocity bias and the velocity dispersion-mass relation of cluster galaxies

We use the Cluster-EAGLE simulations to explore the velocity bias introduced when using galaxies, rather than dark matter particles, to estimate the velocity dispersion of a galaxy cluster, a property known to be tightly correlated with cluster mass. The simulations consist of 30 clusters spanning a mass range 14.0 ≤ log 10 (M 200 c /M ⊙ ) ≤ 15.4, with their sophisticated subgrid physics modelling and high numerical resolution (subkpc gravitational softening), making them ideal for this purpose. We find that selecting galaxies by their total mass results in a velocity dispersion that is 5-10 per cent higher than the dark matter particles. However, selecting galaxies by their stellar mass results in an almost unbiased ( < 5 per cent) estimator of the velocity dispersion. This result holds out to z = 1.5 and is relatively insensitive to the choice of cluster aperture, varying by less than 5 per cent between r 500 c and r 200m . We show that the velocity bias is a function of the time spent by a galaxy inside the cluster environment. Selecting galaxies by their total mass results in a larger bias because a larger fraction of objects have only recently entered the cluster and these have a velocity bias above unity. Galaxies that entered more than 4 Gyr ago become progressively colder with time, as expected from dynamical friction. We conclude that velocity bias should not be a major issue when estimating cluster masses from kinematic methods

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

We present results from a subset of simulations from the "Evolution and Assembly of GaLaxies and their Environments" (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

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

3D-electrical resistivity tomography monitoring of salt transport in homogeneous and layered soil samples

Monitoring transport of dissolved substances in soil deposits is particularly relevant where safety is concerned, as in the case of geo-environmental barriers. Geophysical methods are very appealing, since they cover a wide domain, localising possible preferential flow paths and providing reliable links between geophysical quantities and hydrological variables. This paper describes a 3D laboratory application of Electrical Resistivity Tomography (ERT) used to monitor solute transport processes. Dissolution and transport tests on both homogeneous and heterogeneous samples were conducted in an instrumented oedometer cell. ERT was used to create maps of electrical conductivity of the monitored domain at different time intervals and to estimate concentration variations within the interstitial fluid. Comparisons with finite element simulations of the transport processes were performed to check the consistency of the results. Tests confirmed that the technique can monitor salt transport, infer the hydro-chemical behaviour of heterogeneous geomaterials and evaluate the performances of clay barrier

The masses and density profiles of halos in a LCDM galaxy formation simulation

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 $\Lambda$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

The impact of angular momentum on black hole accretion rates in simulations of galaxy formation

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

Molecular hydrogen abundances of galaxies in the EAGLE simulations

We investigate the abundance of galactic molecular hydrogen (H$_2$) in the "Evolution and Assembly of GaLaxies and their Environments" (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, $\Omega_{\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

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 "Evolution and Assembly of GaLaxies and their Environments" (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 "knee" 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