We investigate the connection between the morphology and internal kinematics of the stellar component of central galaxies with mass $M_\star > {10}^{9.5} {\rm M}_\odot$ in the EAGLE simulations. We compare several kinematic diagnostics commonly used to describe simulated galaxies, and find good consistency between them. We model the structure of galaxies as ellipsoids and quantify their morphology via the ratios of their principal axes. We show that the differentiation of blue star-forming and red quiescent galaxies using morphological diagnostics can be achieved with similar efficacy to the use of kinematical diagnostics, but only if one is able to measure both the flattening and the triaxiality of the galaxy. Flattened oblate galaxies exhibit greater rotational support than their spheroidal counterparts, but there is significant scatter in the relationship between morphological and kinematical diagnostics, such that kinematically-similar galaxies can exhibit a broad range of morphologies. The scatter in the relationship between the flattening and the ratio of the rotation and dispersion velocities ($v/\sigma$) correlates strongly with the anisotropy of the stellar velocity dispersion: at fixed $v/\sigma$, flatter galaxies exhibit greater dispersion in the plane defined by the intermediate and major axes than along the minor axis, indicating that the morphology of simulated galaxies is influenced significantly by the structure of their velocity dispersion. The simulations reveal that this anisotropy correlates with the intrinsic morphology of the galaxy's inner dark matter halo, i.e. the halo's morphology that emerges in the absence of dissipative baryonic physics. This implies the existence of a causal relationship between the morphologies of galaxies and that of their host dark matter haloes.

GalICS 2.0 is a new semianalytic code to model the formation and evolution of galaxies in a cosmological context. N-body simulations based on a Planck cosmology are used to construct halo merger trees, track subhaloes, compute spins and measure concentrations. The accretion of gas onto galaxies and the morphological evolution of galaxies are modelled with prescriptions derived from hydrodynamic simulations. Star formation and stellar feedback are described with phenomenological models (as in other semianalytic codes). GalICS 2.0 computes rotation speeds from the gravitational potential of the dark matter, the disc and the central bulge. As the rotation speed depends not only on the virial velocity but also on the ratio of baryons to dark matter within a galaxy, our calculation predicts a different Tully-Fisher relation from models in which the rotation speed is proportional to the virial velocity. This is why GalICS 2.0 is able to reproduce the galaxy stellar mass function and the Tully-Fisher relation simultaneously. Our results are also in agreement with halo masses from weak lensing and satellite kinematics, gas fractions, the relation between star formation rate (SFR) and stellar mass, the evolution of the cosmic SFR density, bulge-to-disc ratios, disc sizes and the Faber-Jackson relation.

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.