We introduce the Virgo Consortium's Evolution and Assembly of GaLaxies and their Environments (EAGLE) project, a suite of hydrodynamical simulations that follow the formation of galaxies and supermassive black holes in cosmologically representative volumes of a standard Λ cold dark matter universe. 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 active galactic nuclei (AGN) in which thermal energy is injected into the gas without the need to turn off cooling or decouple 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 present-day galaxy stellar mass function and the amplitude of the galaxy-central black hole mass relation, also taking galaxy sizes into account. The observed galaxy stellar mass function is reproduced to ≲ 0.2 dex over the full resolved mass range, 10⁸ < M_*/M_⊙ ≲ 10¹¹, 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 C iv and O vi. While the mass–metallicity relations for gas and stars are consistent with observations for M_* ≳ 10⁹ M_⊙ (M_* ≳ 10¹⁰ M_⊙ at intermediate resolution), they are insufficiently steep at lower masses. For the reference model, the gas fractions and temperatures are too high for clusters of galaxies, but for galaxy groups these discrepancies can be resolved by adopting a higher heating temperature in the subgrid prescription for AGN feedback. The EAGLE simulation suite, which also includes physics variations and higher resolution zoomed-in volumes described elsewhere, constitutes a valuable new resource for studies of galaxy formation.

Aims. The goal of this paper is to analyse the behaviour of the gas-to-dust mass ratio (G/D) of local Universe galaxies over a wide metallicity range. We especially focus on the low-metallicity part of the G/D vs metallicity relation and investigate several explanations for the observed relation and scatter.Methods. We assembled a total of 126 galaxies, covering a 2 dex metallicity range and with 30% of the sample with 12 + log(O/H)≤ 8.0. We homogeneously determined the dust masses with a semi-empirical dust model including submm constraints. The atomic and molecular gas masses have been compiled from the literature. We used two X CO scenarios to estimate the molecular gas mass: the Galactic conversion factor, X CO,MW, and a X CO that depends on the metallicity X CO,Z (∝Z -2). We modelled the observed trend of the G/D with metallicity using two simple power laws (slope of –1 and free) and a broken power law. Correlations with morphological type, stellar masses, star formation rates, and specific star formation rates are also discussed. We then compared the observed evolution of the G/D with predictions from several chemical evolution models and explored different physical explanations for the observed scatter in the G/D values.Results. We find that out of the five tested galactic parameters, metallicity is the main physical property of the galaxy driving the observed G/D. The G/D versus metallicity relation cannot be represented by a single power law with a slope of –1 over the whole metallicity range. The observed trend is steeper for metallicities lower than ~8.0. A large scatter is observed in the G/D values for a given metallicity: in metallicity bins of ~0.1 dex, the dispersion around the mean value is ~0.37 dex. On average, the broken power law reproduces the observed G/D best compared to the two power laws (slope of –1 or free) and provides estimates of the G/D that are accurate to a factor of 1.6. The good agreement of observed values of the G/D and its scatter with respect to metallicity with the predicted values of the three tested chemical evolution models allows us to infer that the scatter in the relation is intrinsic to galactic properties, reflecting the different star formation histories, dust destruction efficiencies, dust grain size distributions, and chemical compositions across the sample. Conclusions. Our results show that the chemical evolution of low-metallicity galaxies, traced by their G/D, strongly depends on their local internal conditions and individual histories. The large scatter in the observed G/D at a given metallicity reflects the impact of various processes occurring during the evolution of a galaxy. Despite the numerous degeneracies affecting them, disentangling these various processes is now the next step.

The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015–2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14 000 deg² of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance.

The Herschel ATLAS is the largest open-time key project that will be carried out on the Herschel Space Observatory. It will survey 570 deg(2) of the extragalactic sky, 4 times larger than all the other Herschel extragalactic surveys combined, in five far-infrared and submillimeter bands. We describe the survey, the complementary multiwavelength data sets that will be combined with the Herschel data, and the six major science programs we are undertaking. Using new models based on a previous submillimeter survey of galaxies, we present predictions of the properties of the ATLAS sources in other wave bands

The primary source of emission of active galactic nuclei (AGNs), the accretion disc, is surrounded by an optically and geometrically thick dusty structure (‘the so-called dusty torus’). The infrared radiation emitted by the dust is nothing but a reprocessed fraction of the accretion disc emission, so the ratio of the torus to the AGN luminosity (L_torus/L_AGN) should corresponds to the fraction of the sky obscured by dust, i.e. the covering factor. We undertook a critical investigation of the L_torus/L_AGN as the dust covering factor proxy. Using state-of-the-art 3D Monte Carlo radiative transfer code, we calculated a grid of spectral energy distributions (SEDs) emitted by the clumpy two-phase dusty structure. With this grid of SEDs, we studied the relation between L_torus/L_AGN and the dust covering factor for different parameters of the torus. We found that in the case of type 1 AGNs the torus anisotropy makes L_torus/L_AGN underestimate low covering factors and overestimate high covering factors. In type 2 AGNs L_torus/L_AGN always underestimates covering factors. Our results provide a novel easy-to-use method to account for anisotropy and obtain correct covering factors. Using two samples from the literature, we demonstrated the importance of our result for inferring the obscured AGN fraction. We found that after the anisotropy is properly accounted for, the dust covering factors show very weak dependence on L_AGN, with values in the range of ≈0.6–0.7. Our results also suggest a higher fraction of obscured AGNs at high luminosities than those found by X-ray surveys, in part owing to the presence of a Compton-thick AGN population predicted by population synthesis models.