We investigate correlations between different physical properties of star-forming galaxies in the 'Evolution and Assembly of GaLaxies and their Environments' (EAGLE) cosmological hydrodynamical simulation suite over the redshift range 0 ≤ z ≤ 4.5. A principal component analysis reveals that neutral gas fraction (fgas,neutral), stellar mass (Mstellar) 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 fgas, neutral–Mstellar–SFR with little scatter. The location of this plane varies little with redshift, whereas galaxies themselves move along the plane as their fgas, 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 fgas, neutral, or from the Mstellar and SFR. We argue that the appearance of this 'Fundamental Plane of star formation' is a consequence of self-regulation, with the plane'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 ≲ z ≲ 3, 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's predictions.
We explore scaling relations between the physical properties of spatially resolved regions within the galaxies that emerge in the Evolution and Assembly of GaLaxies and their Environments (EAGLE) hydrodynamical, cosmological simulations. Using 1 kpc-scale spaxels, we compute the relationships between the star formation rate and stellar mass surface densities, i.e. the spatially resolved star-forming main sequence (rSFMS), and between the gas metallicity and the stellar mass surface density, i.e. the spatially resolved mass–metallicity relation (rMZR). We compare these with observed relations derived from integral field unit surveys and galaxy imaging. EAGLE reproduces the slope of the local (z ≈ 0.1) rSFMS well, but with a ≈− 0.15 dex offset, close to that found for the galaxy-integrated relation. The shape of the rMZR agrees reasonably well with observations, replicating the characteristic turnover at high surface density, which we show is due to AGN feedback. The residuals of the rSFMS and rMZR are negatively (positively) correlated at low (high) surface density. The rSFMS becomes shallower as the simulation evolves from z = 2 to 0.1, a manifestation of inside-out galaxy formation. The shape of the rMZR also exhibits dramatic evolution, from a convex profile at z = 2 to the observed concave profile at z = 0.1, such that the gas in regions of high stellar density is more enriched at higher redshift. The redshift independence of the relationship between the galaxy-wide gas fraction and metallicity in EAGLE galaxies is not preserved on 1 kpc scales, implying that chemical evolution is non-local due to the transport of gas and metals within galaxies.
Sonification is the technique of representing data with sound, with potential applications in astronomy research for aiding discovery and accessibility. Several astronomy-focused sonification tools have been developed; however, efficacy testing is extremely limited. We performed testing of astronify, a prototype tool for sonification functionality within the Barbara A. Mikulski Archive for Space Telescopes (MAST). We created synthetic light curves containing zero, one, or two transit-like signals with a range of signal-to-noise ratios (SNRs=3-100) and applied the default mapping of brightness to pitch. We performed remote testing, asking participants to count signals when presented with light curves as a sonification, visual plot, or combination of both. We obtained 192 responses, of which 118 self-classified as experts in astronomy and data analysis. For high SNRs (=30 and 100), experts and non-experts performed well with sonified data (85-100% successful signal counting). At low SNRs (=3 and 5) both groups were consistent with guessing with sonifications. At medium SNRs (=7 and 10), experts performed no better than non-experts with sonifications but significantly better (factor of ~2-3) with visuals. We infer that sonification training, like that experienced by experts for visual data inspection, will be important if this sonification method is to be useful for moderate SNR signal detection within astronomical archives and broader research. Nonetheless, we show that even a very simple, and non-optimised, sonification approach allows users to identify high SNR signals. A more optimised approach, for which we present ideas, would likely yield higher success for lower SNR signals.
We study present-day galaxy clustering in the eagle cosmological hydrodynamical simulation. eagle's galaxy formation parameters were calibrated to reproduce the redshift z = 0.1 galaxy stellar mass function, and the simulation also reproduces galaxy colours well. The simulation volume is too small to correctly sample large-scale fluctuations and we therefore concentrate on scales smaller than a few mega parsecs. We find very good agreement with observed clustering measurements from the Galaxy And Mass Assembly (GAMA) survey, when galaxies are binned by stellar mass, colour or luminosity. However, low-mass red galaxies are clustered too strongly, which is at least partly due to limited numerical resolution. Apart from this limitation, we conclude that eagle galaxies inhabit similar dark matter haloes as observed GAMA galaxies, and that the radial distribution of satellite galaxies, as a function of stellar mass and colour, is similar to that observed as well.
The cosmic spectral energy distribution (CSED) is the total emissivity as a function of wavelength of galaxies in a given cosmic volume. We compare the observed CSED from the UV to the submm to that computed from the EAGLE cosmological hydrodynamical simulation, post-processed with stellar population synthesis models and including dust radiative transfer using the SKIRT code. The agreement with the data is better than 0.15 dex over the entire wavelength range at redshift $z=0$, except at UV wavelengths where the EAGLE model overestimates the observed CSED by up to a factor 2. Global properties of the CSED as inferred from CIGALE fits, such as the stellar mass density, mean star formation density, and mean dust-to-stellar-mass ratio, agree to within better than 20 per cent. At higher redshift, EAGLE increasingly underestimates the CSED at optical-NIR wavelengths with the FIR/submm emissivity underestimated by more than a factor of 5 by redshift $z=1$. We believe that these differences are due to a combination of incompleteness of the EAGLE-SKIRT database, the small simulation volume and the consequent lack of luminous galaxies, and our lack of knowledge on the evolution of the characteristics of the interstellar dust in galaxies. The impressive agreement between the simulated and observed CSED at lower $z$ confirms that the combination of EAGLE and SKIRT dust processing yields a fairly realistic representation of the local Universe.
We have created a show about the Solar System, freely available for both planetariums and home viewing, where objects in space are represented with sound as well as with visuals. For example, the audience listens to the stars appear above the European Southern Observatory's Very Large Telescope and they hear the planets orbit around their heads. The aim of this show is that it can be enjoyed and understood, irrespective of level of vision. Here we describe how we used our new computer code, STRAUSS, to convert data into sound for the show. We also discuss the lessons learnt during the design of the show, including how it was imperative to obtain a range of diverse perspectives from scientists, a composer and representatives of the blind and vision impaired community.
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 iMaNGA project uses a forward-modelling approach to compare the predictions of cosmological simulations with observations from SDSS-IV/MaNGA. We investigate the dependency of age and metallicity radial gradients on galaxy morphology, stellar mass, stellar surface mass density ($\Sigma_*$), and environment. The key of our analysis is that observational biases affecting the interpretation of MaNGA data are emulated in the theoretical iMaNGA sample. The simulations reproduce the observed global stellar population scaling relations with positive correlations between galaxy mass and age/metallicity quite well and also produce younger stellar populations in late-type in agreement with observations. We do find interesting discrepancies, though, that can inform the physics and further development of the simulations. Ages of spiral galaxies and low-mass ellipticals are overestimated by about 2-4 Gyr. Radial metallicity gradients are steeper in iMaNGA than in MaNGA, a discrepancy most prominent in spiral and lenticular galaxies. Also, the observed steepening of metallicity gradients with increasing galaxy mass is not well matched by the simulations. We find that the theoretical radial profiles of surface mass density $\Sigma_*$ are steeper than in observations except for the most massive galaxies. In both MaNGA and iMaNGA [Z/H] correlates with $\Sigma_*$, however, the simulations systematically predict lower [Z/H] by almost a factor of 2 at any $\Sigma_*$. Most interestingly, for galaxies with stellar mass $\log M_*\leq 10.80 M_\odot$ the MaNGA data reveal a positive correlation between galaxy radius and [Z/H] at fixed $\Sigma_*$, which is not recovered in iMaNGA. Finally, the dependence on environmental density is negligible in both the theoretical iMaNGA and the observed MaNGA data.
We examine the evolution of intrinsic u-r colours of galaxies in the EAGLE cosmological hydrodynamical simulations, which has been shown to reproduce the observed redshift z=0.1 colour-magnitude distribution well. The median u-r of star-forming ('blue cloud') galaxies reddens by 1 mag from z=2 to 0 at fixed stellar mass, as their specific star formation rates decrease with time. A red sequence starts to build-up around z=1, due to the quenching of low-mass satellite galaxies at the faint end, and due to the quenching of more massive central galaxies by their active galactic nuclei (AGN) at the bright end. This leaves a dearth of intermediate-mass red sequence galaxies at z=1, which is mostly filled in by z=0. We quantify the time-scales of colour transition due to satellite and AGN quenching, finding that most galaxies spend less than 2 Gyr in the 'green valley'. On examining the trajectories of galaxies in a colour-stellar mass diagram, we identify three characteristic tracks that galaxies follow (quiescently star-forming, quenching and rejuvenating galaxies) and quantify the fraction of galaxies that follow each track.
We present results from a survey of galaxies in the fields of six z>3 Damped Lyman alpha systems (DLAs) using the Multi Unit Spectroscopic Explorer (MUSE) at the Very Large Telescope (VLT). We report a high detection rate of up to ~80% of galaxies within 1000 km/s from DLAs and with impact parameters between 25 and 280 kpc. In particular, we discovered 5 high-confidence Lyman alpha emitters associated with three DLAs, plus up to 9 additional detections across five of the six fields. The majority of the detections are at relatively large impact parameters (>50 kpc) with two detections being plausible host galaxies. Among our detections, we report four galaxies associated with the most metal-poor DLA in our sample (Z/Z_sun = -2.33), which trace an overdense structure resembling a filament. By comparing our detections with predictions from the Evolution and Assembly of GaLaxies and their Environments (EAGLE) cosmological simulations and a semi-analytic model designed to reproduce the observed bias of DLAs at z>2, we conclude that our observations are consistent with a scenario in which a significant fraction of DLAs trace the neutral regions within halos with a characteristic mass of 10^11-10^12 M_sun, in agreement with the inference made from the large-scale clustering of DLAs. We finally show how larger surveys targeting ~25 absorbers have the potential of constraining the characteristic masses of halos hosting high-redshift DLAs with sufficient accuracy to discriminate between different models.
'%3e%3cpath fill='%23012169' d='M0 0v30h60V0z'/%3e%3cpath stroke='%23fff' stroke-width='6' d='m0 0 60 30m0-30L0 30'/%3e%3cpath stroke='%23C8102E' stroke-width='4' d='m0 0 60 30m0-30L0 30' clip-path='url(%23b)'/%3e%3cpath stroke='%23fff' stroke-width='10' d='M30 0v30M0 15h60'/%3e%3cpath stroke='%23C8102E' stroke-width='6' d='M30 0v30M0 15h60'/%3e%3c/g%3e%3c/svg%3e)
'/%3e%3c/defs%3e%3cpath fill='%23012169' d='M0 0h10080v5040H0z'/%3e%3cpath stroke='%23fff' stroke-width='.6' d='m0 0 6 3m0-3L0 3' clip-path='url(%23b)' transform='scale(840)'/%3e%3cpath stroke='%23e4002b' stroke-width='.4' d='m0 0 6 3m0-3L0 3' clip-path='url(%23c)' transform='scale(840)'/%3e%3cpath stroke='%23fff' stroke-width='840' d='M2520 0v2520M0 1260h5040'/%3e%3cpath stroke='%23e4002b' stroke-width='504' d='M2520 0v2520M0 1260h5040'/%3e%3cg fill='%23fff'%3e%3cuse xlink:href='%23d' x='2520' y='3780'/%3e%3cuse xlink:href='%23a' x='7560' y='4200'/%3e%3cuse xlink:href='%23a' x='6300' y='2205'/%3e%3cuse xlink:href='%23a' x='7560' y='840'/%3e%3cuse xlink:href='%23a' x='8680' y='1869'/%3e%3cuse xlink:href='%23e' x='8064' y='2730'/%3e%3c/g%3e%3c/svg%3e)
