We investigate the origin of the Hubble sequence by analysing the evolution of the kinematic morphologies of central galaxies in the EAGLE cosmological simulation. By separating each galaxy into disk and spheroidal stellar components and tracing their evolution along the merger tree, we find that the morphology of galaxies follows a common evolutionary trend. We distinguish three phases of galaxy formation. These phases are determined primarily by mass, rather than redshift. For M_star < 10^{9.5} M_sun galaxies grow in a disorganised way, resulting in a morphology that is dominated by random stellar motions. This phase is dominated by in-situ star formation, partly triggered by mergers. In the mass range 10^{9.5} M_sun < M_star < 10^{10.5} M_sun galaxies evolve towards a disk-dominated morphology, driven by in-situ star formation. The central spheroid (i.e. the bulge) at z = 0 consists mostly of stars that formed in-situ, yet the formation of the bulge is to a large degree associated with mergers. Finally, at M_star > 10^{10.5} M_sun growth through in-situ star formation slows down considerably and galaxies transform towards a more spheroidal morphology. This transformation is driven more by the buildup of spheroids than by the destruction of disks. Spheroid formation in these galaxies happens mostly by accretion at large radii of stars formed ex-situ (i.e. the halo rather than the bulge).

Galaxy evolution can be studied observationally by linking progenitor and descendant galaxies through an evolving cumulative number density selection. This procedure can reproduce the expected evolution of the median stellar mass from abundance matching. However, models predict an increasing scatter in main progenitor masses at higher redshifts, which makes galaxy selection at the median mass unrepresentative. Consequently, there is no guarantee that the evolution of other galaxy properties deduced from this selection are reliable. Despite this concern, we show that this procedure approximately reproduces the evolution of the average stellar density profile of main progenitors of M = 10^11.5 Msun galaxies, when applied to the EAGLE hydrodynamical simulation. At z > 3.5 the aperture masses disagree by about a factor two, but this discrepancy disappears when we include the expected scatter in cumulative number densities. The evolution of the average density profile in EAGLE broadly agrees with observations from UltraVISTA and CANDELS, suggesting an inside-out growth history for these massive galaxies over 0 < z < 5. However, for z < 2 the inside-out growth trend is stronger in EAGLE. We conclude that cumulative number density matching gives reasonably accurate results when applied to the evolution of the mean density profile of massive galaxies.

We investigate the relation between kinematic morphology, intrinsic colour and stellar mass of galaxies in the EAGLE cosmological hydrodynamical simulation. We calculate the intrinsic u-r colours and measure the fraction of kinetic energy invested in ordered corotation of 3562 galaxies at z=0 with stellar masses larger than $10^{10}M_{\odot}$. We perform a visual inspection of gri-composite images and find that our kinematic morphology correlates strongly with visual morphology. EAGLE produces a galaxy population for which morphology is tightly correlated with the location in the colour- mass diagram, with the red sequence mostly populated by elliptical galaxies and the blue cloud by disc galaxies. Satellite galaxies are more likely to be on the red sequence than centrals, and for satellites the red sequence is morphologically more diverse. These results show that the connection between mass, intrinsic colour and morphology arises from galaxy formation models that reproduce the observed galaxy mass function and sizes.

In this paper, we use stacking analysis to trace the mass-growth, colour evolution, and structural evolution of present-day massive galaxies ($\log(M_{*}/M_{\odot})=11.5$) out to $z=5$. We utilize the exceptional depth and area of the latest UltraVISTA data release, combined with the depth and unparalleled seeing of CANDELS to gather a large, mass-selected sample of galaxies in the NIR (rest-frame optical to UV). Progenitors of present-day massive galaxies are identified via an evolving cumulative number density selection, which accounts for the effects of merging to correct for the systematic biases introduced using a fixed cumulative number density selection, and find progenitors grow in stellar mass by $\approx1.5~\mathrm{dex}$ since $z=5$. Using stacking, we analyze the structural parameters of the progenitors and find that most of the stellar mass content in the central regions was in place by $z\sim2$, and while galaxies continue to assemble mass at all radii, the outskirts experience the largest fractional increase in stellar mass. However, we find evidence of significant stellar mass build up at $r<3~\mathrm{kpc}$ beyond $z>4$ probing an era of significant mass assembly in the interiors of present day massive galaxies. We also compare mass assembly from progenitors in this study to the EAGLE simulation and find qualitatively similar assembly with $z$ at $r<3~\mathrm{kpc}$. We identify $z\sim1.5$ as a distinct epoch in the evolution of massive galaxies where progenitors transitioned from growing in mass and size primarily through in-situ star formation in disks to a period of efficient growth in $r_{e}$ consistent with the minor merger scenario.

We investigate the effect of a metallicity-dependent stellar initial mass function (IMF), as deduced observationally by Martin-Navarro et al. (2015c), on the inferred stellar masses and star formation rates (SFRs) of a representative sample of 186,886 SDSS galaxies. Relative to a Chabrier IMF, for which we show the implied masses to be close to minimal, the inferred masses increase in both the low- and high-metallicity regimes due to the addition of stellar remnants and dwarf stars, respectively. The resulting galaxy stellar mass function (GSMF) shifts toward higher masses by 0.5 dex, without affecting the high-mass slope (and thus the need for effective quenching). The implied low-redshift SFR density increases by an order of magnitude. However, these results depend strongly on the assumed IMF parametrisation, which is not directly constrained by the observations. Varying the low-end IMF slope instead of the high-end IMF slope, while maintaining the same dwarf-to-giant ratio, results in a much more modest GSMF shift of 0.2 dex and a 10 per cent increase in the SFR density relative to the Chabrier IMF. A bottom-heavy IMF during the late, metal-rich evolutionary stage of a galaxy would help explain the rapid quenching and the bimodality in the galaxy population by on the one hand making galaxies less quenched (due to the continued formation of dwarf stars) and on the other hand reducing the gas consumption timescale. We conclude that the implications of the observational evidence for a variable IMF could vary from absolutely dramatic to mild but significant.

Cumulative number density matching of galaxies is a method to observationally connect descendent galaxies to their typical main progenitors at higher redshifts and thereby to assess the evolution of galaxy properties. The accuracy of this method is limited due to galaxy merging and scatter in the stellar mass growth history of individual galaxies. Behroozi et al. (2013) have introduced a refinement of the method, based on abundance matching of observed galaxies to the Bolshoi dark-matter-only simulation. The EAGLE cosmological hydro-simulation is well suited to test this method, because it reproduces the observed evolution of the galaxy stellar mass function and the passive fraction. We find agreement with the Behroozi et al. (2013) method for the complete sample of main progenitors of z = 0 galaxies, but we also find a strong dependence on the current star formation rate. Passive galaxies with a stellar mass up to 10^10.75 Msun have a completely different median mass history than active galaxies of the same mass. This difference persists if we only select central galaxies. This means that the cumulative number density method should be applied separately to active and passive galaxies. Even then, the typical main progenitor of a z = 0 galaxy already spans two orders of magnitude in stellar mass at z = 2.

The ATLAS3D Survey has reported evidence for a non-universal stellar initial mass function (IMF) for early type galaxies (ETGs) (Cappellari et al. 2012, 2013b,a). The IMF was constrained by comparing stellar mass measurements from kinematic data with those from spectral energy distribution (SED) fitting. Here we investigate possible effects of scatter in the reported stellar mass measurements and their potential impact on the IMF determination. We find that a trend of the IMF mismatch parameter with the kinematic mass to light ratio, comparable to the trend observed by Cappellari et al. (2012), could arise if the Gaussian errors of the kinematic mass determination are typically 30%. Without additional data, it is hard to separate between the option that the IMF has a true large intrinsic variation or the option that the errors in the determination are larger than anticipated. A correlation of the IMF with other properties would help to make this distinction, but no strong correlation has been found yet. The strongest correlation is with velocity dispersion. However, it has a large scatter and the correlation depends on sample selection and distance measurements. The correlation with velocity dispersion could be partly caused by the colour-dependent calibration of the surface brightness fluctuation distances of Tonry et al. (2001). We find that the K-band luminosity limited ATLAS3D Survey is incomplete for the highest M/L galaxies below 10^10.3 M_sun. There is a significant IMF - velocity dispersion trend for galaxies with SED masses above this limit, but no trend for galaxies with kinematic masses above this limit. We also find an IMF trend with distance, but no correlation between nearest neighbour ETGs, which excludes a large environmental dependence. Our findings do not rule out the reported IMF variations, but they suggest that further study is needed.

Context. Analytical and numerical analysis of the SimpleX radiative transfer algorithm, which features transport on a Delaunay triangulation. Aims. Verify whether the SimpleX radiative transfer algorithm conforms to mathematical expectations, to develop error analysis and present improvements upon earlier versions of the code. Methods. Voronoi-Delaunay tessellation, classical Markov theory. Results. Quantitative description of the error properties of the SimpleX method. Numerical validation of the method and verification of the analytical results. Improvements in accuracy and speed of the method. Conclusions. It is possible to transport particles such as photons in a physically correct manner with the SimpleX algorithm. This requires the use of weighting schemes or the modification of the point process underlying the transport grid. We have explored and applied several possibilities.