## L etter to the E ditor

## The mass-size plane of EAGLE galaxies

M. S. Rosito , P. B. Tissera 1 2 3 4 ; ; , S. E. Pedrosa 1 5 ; , C. D. P. Lagos 6 7 ;

- 1 Instituto de Astronomía y Física del Espacio, CONICET-UBA, Casilla de Correos 67, Suc. 28, 1428, Buenos Aires, Argentina.
- 2 Departamento de Ciencias Físicas, 700 Fernandez Concha, Universidad Andrés Bello, Santiago, Chile.
- 3 Corresponding Investigator, IATE-CONICET, Laprida 927, Córdoba, Argentina.
- 4 Millenium Institute of Astronomy, 700 Fernandez Concha, Universidad Andrés Bello, Santiago, Chile.
- 5 Departamento de Física Teórica, Facultad de Ciencias, Universidad Autónoma de Madrid.
- 6 International Centre for Radio Astronomy Research (ICRAR), M468, University of Western Australia, 35 Stirling Hwy, Crawley, WA6009, Australia.
- 7 ARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D).

August 2, 2019

## ABSTRACT

Context. Current observational results show that both late-and-early-type galaxies follow tight mass-size planes, on which physical properties such as age, velocity dispersion and metallicities correlate with the scatter on the plane.

Aims.

We study the mass-size plane of galaxies in cosmological hydrodynamical simulations, as a function of velocity dispersion, age, chemical abundances, ellipticity and spin parameters with the aim at assessing to what extent the current cosmological paradigm

can reproduce these observations and provide a physical interpretation of them.

Methods. We select a sample of well-resolved galaxies from the (100 Mpc) 3 simulation of the EAGLE Project. This sample is composed by 508 spheroid-dominated galaxies and 1213 disc-dominated galaxies. The distributions of velocity dispersion, age, metallicity indicators and gradients and spin parameters across the mass-size plane are analysed. Furthermore, we study the relation between shape and kinematic parameters. The results are compared with observations.

Results. The mass-weighted ages of the EAGLE galaxies are found to vary along lines of constant velocity dispersion on the mass-size plane, except for galaxies with velocity dispersion larger than GLYPH&lt;24&gt; 150km s GLYPH&lt;0&gt; 1 . Negative age gradients tend to be found in extended disc galaxies in agreement with observations. However, the age distributions of early-type galaxies show a larger fraction with inverted radial profiles. The distribution of metallicity gradients does not show any clear dependence on this plane. Galaxies with similar spin parameters ( GLYPH&lt;21&gt; ) display larger sizes as their dynamical masses increase. Stellar-weighted ages are found to be good proxies for GLYPH&lt;21&gt; in galaxies with low ellipticity ( " ). A bimodal distribution of GLYPH&lt;21&gt; is found so that the highGLYPH&lt;21&gt; peak is dominated by discs with young SPs while the second peak is mainly populated by slow rotators ( GLYPH&lt;21&gt; &lt; 0 2) with old stars. Our findings suggest that the physical : processes which regulate the star formation histories in galaxies might also a GLYPH&lt;11&gt; ect the angular moment budgets of gas and stars and as a consequence, their morphology.

Key words. galaxies: formation - galaxies: elliptical and lenticular, cD - galaxies: abundances - galaxies: kinematics and dynamics

## 1. Introduction

The formation and evolution of galaxies is a complex interplay of di GLYPH&lt;11&gt; erent physical processes such as rapid gas collapse, disc mergers, minor mergers and secular evolution, among others. These processes contribute to regulate the star formation activity and also to mix the chemical elements, polluting the interstellar medium (ISM) of galaxies. The formation of spiral galaxies can be understood within the scenario of global conservation of specific angular momentum (Fall &amp; Romanowsky 2013). However, in a hierarchical universe, mergers and interactions as well as other environmental e GLYPH&lt;11&gt; ects may redistribute the angular momentum of gas and stars, modifying the morphologies and overall properties of galaxies (e.g. Pedrosa &amp; Tissera 2015; Lagos et al. 2017, 2018b). Ellipticals are expected to have a complex formation history where the mentioned processes might also take place, albeit acting with di GLYPH&lt;11&gt; erent e GLYPH&lt;14&gt; ciencies and characteristics (Clauwens et al. 2018). In particular, the mass-size plane, which combines two fundamental properties of galaxies, could hold relevant information for the understanding of galaxy formation (Li et al. 2018, hereafter Li2018).

The galaxy assembly process is expected to leave imprints on chemical patterns on the stellar populations (SPs). As baryons are transformed into stars, stellar evolution takes place and new chemical elements are synthesised and ejected into the ISM, at di GLYPH&lt;11&gt; erent stages of evolution of the stellar progenitors. The newborn SPs a GLYPH&lt;11&gt; ect the chemical abundances of the ISM, building up the chemical history of galaxies. The growth of stellar mass and galaxy size via star formation, gas accretion or and galaxy / mergers are therefore expected to be related to the properties of the SPs such as chemical abundances, kinematics and ages.

With the recent help of integral-field spectroscopy (IFS) techniques, galaxy surveys such as MaNGA (Bundy et al. 2015), CALIFA (Sánchez et al. 2012) and SAMI (Bryant et al. 2015) have gathered information which allows the construction of the mass-size plane populated with enough galaxies to study the interdependencies on those properties with higher statistical signals than it was possible before (e.g. Ryden et al. 2001). Li2018 analysed the mass-size plane of early-type galaxies (ETGs) and spiral galaxies (LTGs) in the MaNGA survey, finding that both types of galaxies are on tight mass-size planes, which are consistent with the predictions of the virial theorem. Metallicities and

ages are found to vary systematically on the mass-size plane, along the direction of constant velocity dispersion. Galaxies of di GLYPH&lt;11&gt; erent morphologies follow mass-size planes with distinct dependences on other galaxy properties. Overall, Li2018 found that the velocity dispersion may be used as a proxy for the bulge-tototal ( B T = ) ratio (e.g. Cappellari et al. 2013). In their sample, ETGs of higher velocity dispersion are typically associated with older and higher metallicities SPs. Additionally, van de Sande et al. (2018) analysed the link between the rotation-to-dispersion velocity ratio of a galaxy and its ellipticity, finding that stellar age follows ellipticity ( " ) very well in oblate rotating spheroids. Graham et al. (2018) reported a bimodal distribution in the ( GLYPH&lt;21&gt;; " ) plane of galaxies in the MaNGA survey, with the secondary peak being dominated by slow rotator galaxies.

Cosmological simulations that include chemical evolution are fundamental tools to study the relation between dynamical, structural and chemical properties (e.g. Mosconi et al. 2001; Lia et al. 2002; Tissera et al. 2012; Taylor &amp; Kobayashi 2017). These simulations show that the formation of LTGs could be explained if overall global angular momentum conservation takes place. Hence, they mostly grow inside-out, leading to negative age and metallicity gradients (Tissera et al. 2016, 2019). These trends may be disturbed if processes such as galaxy mergers, bars and migration take place. The formation of ETGs could involve the action of a variety of mechanisms at di GLYPH&lt;11&gt; erent stages of their evolution (Naab 2013), such as minor and major mergers, as well as early filamentary, cold gas accretion. Indeed, galaxy mergers are considered the most e GLYPH&lt;14&gt; cient mechanisms to redistribute angular momentum and modify the morphologies of galaxies, among other galaxy properties (e.g. Bois et al. 2011; Perez et al. 2013; Naab et al. 2014; Lagos et al. 2018b). Lagos et al. (2018a) found that 70 per cent of the slow rotators have experienced at least one important merger event in the EAGLE simulations. They also found that wet mergers tend to increase the GLYPH&lt;21&gt; while dry mergers have the opposite e GLYPH&lt;11&gt; ect. Recently, van de Sande et al. (2019) analysed the structural, dynamical and stellar age populations of galaxies in several large-volume simulations, including the EAGLE simulations, and confront them with IFS observations. These authors reported a good match of the mass-size plane between observations and simulated galaxies but did not find a correlation between " and age that resembled observations.

In this paper, we analyse the mass-size relation of simulated galaxies extracted from the EAGLE simulations. The EAGLE simulations is able to reproduce relatively well the observed diversity of galaxy morphologies (Frenk et al. 2018, and references therein). We classify the analysed galaxies in spheroidaldominated (hereafter, E-SDGs) and disc-dominated (E-DDGs) systems according to the fraction of their stellar mass that is rotationally supported. In a previous work, Rosito et al. (2018) studied the structural relations and the mass growth history (MGHs) of dispersion-dominated galaxies selected from the EAGLEproject, finding that they are able to reproduce structural relations such as the mass-size relation and the fundamental plane. The MGHs analysis of Rosito et al. (2018) shows a trend for a coeval formation of the SPs in massive E-SDGs and for a weak outside-in formation in low-mass E-SDGs, driven by a rejuvenating star formation activity that took place in the central regions. On the other hand, the E-DDG galaxies tend to form inside-out as expected. We use 3D (real) and 2D (projected) distributions in order to analyse the intrinsic distributions as well as to compare with observations.

This paper is organised as follows. Section 2 describes the main characteristics of the simulations. Section 3 explores the mass-size plane as a function of galaxy properties. Finally, Section4 summarises our main findings.

## 2. Simulated galaxies

The EAGLE project is a suite of cosmological simulations 1 run with a modified version of GADGET-3 code that includes radiative cooling, star formation, chemical evolution and stellar and AGN feedback (Schaye et al. 2015). Details on the subgrid physics adopted to simulate these processes can be found in Schaye et al. (2015), Crain et al. (2015) and Rosas-Guevara et al. (2016). The simulations are produced assuming a GLYPH&lt;3&gt; Cold Dark Matter scenario with the Planck Cosmology (Planck Collaboration et al. 2014a,b). For this work, we use the 100 Mpc box-sized reference simulation that allows us to explore a large statistical sample of galaxies with di GLYPH&lt;11&gt; erent morphologies.

The classification of EAGLE galaxies into E-DDGs and ESDG has been done based on a dynamical criteria to define the rotational-and-dispersion-dominated stellar components (i.e. disc and bulge, respectively) as explained in Tissera et al. (2019). Hence, all galaxies have been rotated so that the plane of rotation is perpendicular to the total angular momentum vector. Galaxies are classified according to their B T = ratio, which directly quantifies the stellar mass fraction that is dominated by velocity dispersion. For consistency with Rosito et al. (2018), we adopt B T = = 0 5 to separate disc (E-DDG) from bulge (E-: SDG) dominated galaxies. The analysed sample comprises 508 E-SDGs and 1213 E-DDGs resolved with more than 10,000 star particles within the optical radius 2 . They are all central galaxies, which corresponds to the main galaxy in the potential well of a halo. Satellite galaxies are not analysed in this work.

We define the 3D half-mass radius R hm as the one that encloses 50 per cent of the stellar mass of the simulated galaxies. For each simulated galaxy, the median value of chemical abundances (O H and [O Fe]) are calculated within / / R hm. We also estimate the dynamical mass ( M dyn ), assuming virialization for both the observations and the simulations, as shown in Rosito et al. (2018).

The chemical abundance and age profiles are estimated within three times the gravitational softening (0.7 kpc) and R hm and normalised by R hm. Linear regression fits (in logarithm space) are performed to the normalized profiles, following the procedure given by Li2018 (within the same radial interval mentioned above).

We also estimated the projected half-mass radius, R 2D hm . The GLYPH&lt;21&gt; and " parameters are taken from Lagos et al. (2018a). These authors measured the r-band luminosity-weighted line-of-sight velocity stellar GLYPH&lt;21&gt; and " at the projected half-r-band luminosity radius, using random inclinations to best mimic observational procedures.

## 3. Analysis

We explore the distribution of velocity dispersion, age, metallicity and GLYPH&lt;21&gt; of the EAGLE galaxies on the mass-size plane. It is important to bear in mind that observed and simulated samples have been constructed using di GLYPH&lt;11&gt; erent selection criteria. Li2018 use the Sérsic index to distinguish ETGs and LTGs so that galaxies with n Sersic &gt; 2 5 are classified as early types. They also use : the dynamical mass and the major semi-axis of the half-light

1 We use the publicly available database by McAlpine et al. (2016).

2 The optical radius is defined as the one that encloses 80 per cent of the baryonic mass of the galaxy.

isophote. For the selected EAGLE sample, we use the dynamical B T = ratio to distinguish morphologies and the half-mass radius as the characteristic size, consistently with Rosito et al. (2018). We also use the mass-size relations for passive and active starforming (SF) galaxies reported by van der Wel et al. (2014). Passive SF galaxies are compared to E-SDGs while active SF are confronted with E-DDGs.

Figure 1 shows the stellar mass-size plane as a function of GLYPH&lt;27&gt; e for E-SDGs, E-DDGs and all galaxies together. In this case, we consider R 2D hm in order to achieve a better comparison with the observations. We use the two-dimensional Locally Weighted Regression method to obtain smoothed distributions 3 (LOESS Cleveland &amp; Devlin 1988). In Appendix B, figures with the real distributions are included for comparison. The stellar mass-size plane shows the expected trend for increasing sizes with stellar mass within the observational range reported by van der Wel et al. (2014). We find the best agreement for intermediate stellar mass galaxies GLYPH&lt;24&gt; 10 10 3 : M for both E-SDGs and E-DDGs (see GLYPH&lt;12&gt; also van de Sande et al. 2019).

Globally, it can be seen that larger size galaxies have lower GLYPH&lt;27&gt; e at a given stellar mass. E-DDGs increase systematically their stellar mass and size at fixed values of GLYPH&lt;27&gt; e . However, GLYPH&lt;27&gt; e is only a proxy of stellar mass in the case of E-SDGs. At a similar GLYPH&lt;27&gt; e , EDDGs are more extended than E-SDGs, possibly due to di GLYPH&lt;11&gt; erent assembly histories.

In Fig. 2, the dynamical mass-size plane is shown as a function of the stellar mass-weighted ages of the SPs following Li2018. The range of dynamical masses covered by the selected samples is smaller than the one provided by the observations of Li2018. The lack of extended massive discs as well as the narrow range of dynamical masses could be partially due to the restricted simulated volume, which is GLYPH&lt;25&gt; 80 smaller than the volume covered by MaNGA.

As can be seen, galaxies with similar median stellar ages move approximately along lines of constant GLYPH&lt;27&gt; e , in agreement with observations. However, the trend has a slightly flatter slope for galaxies with GLYPH&lt;27&gt; e &gt; 150km s GLYPH&lt;0&gt; 1 . Young SPs dominate in galaxies with low stellar masses and GLYPH&lt;27&gt; e . For a given dynamical mass, galaxies with overall old SPs are more compact, regardless of their morphology. These trends are in qualitative agreement with observational results shown in Li2018. We estimate the youngest galaxies to have median stellar ages of GLYPH&lt;24&gt; 3Gyr and GLYPH&lt;24&gt; 3 7 Gyr : for E-DDGs and E-SDGs, and the oldest ones of GLYPH&lt;24&gt; 12 5 for both : type of galaxies. The median ages of the both samples are GLYPH&lt;24&gt; 6 8 : Gyr and GLYPH&lt;24&gt; 8 7 Gyr, respectively. :

Fig. 3 shows the distribution of age and stellar metallicity gradients across the dynamical mass-size plane. As can be seen for the upper panels, massive E-SDGs tend to have shallow positive age gradients which are in global qualitative agree-/ ment with Li2018. Low-mass galaxies show a larger variety of age gradients that are not reported by Li2018, probably as a result of their later assembly histories (see Rosito et al. 2018). For E-DDGs, there is a clear trend for more extended galaxies to be more massive at a given age gradient. As the statistics are dominated by E-DDGs, when we analyse all the galaxies together, we find the trends to be similar to those of the E-DDGs. We note that the behaviour for E-DDGs is in good agreement with expectations of the inside-out formation model for disc galaxies. However, we acknowledge a lack of extended, massive galaxies in the EAGLE sample (van de Sande et al. 2019).

3 The colour-bars are defined by using the first and third quartile as limits.

The metallicity gradients are weakly correlated with dynamical mass (Tissera et al. 2019) and this is reflected on the mass-size plane shown in the lower panels in Fig 3. For the E-DDGs, metallicity gradients are globally negative. There is a weak change in the slopes at about M dyn GLYPH&lt;24&gt; 10 10 5 : M . This trend GLYPH&lt;12&gt; is determined mainly by a fraction of the EAGLE galaxies that have weak positive metallicity gradients around / GLYPH&lt;27&gt; e GLYPH&lt;24&gt; 150 km s GLYPH&lt;0&gt; 1 .

We also analyse the mass-size plane as a function of the spin parameter GLYPH&lt;21&gt; and the relation of the latter with " and the median stellar ages of galaxies. In Fig. 4, we show the distribution of GLYPH&lt;21&gt; on the dynamical mass-size plane (upper panels) in the same samples of Fig. 2 and Fig. 3 for comparison. As it can be seen from this figure, there is a clear dependence on GLYPH&lt;21&gt; . Overall, E-SDGs have lower GLYPH&lt;21&gt; than E-SDGs. For both morphological types, at fixed M dyn , GLYPH&lt;21&gt; varies smoothly with R hm (note that di GLYPH&lt;11&gt; erent colour-bars are used in each panel in order to highlight the trends). This figure shows that galaxies with similar GLYPH&lt;21&gt; are located on tracks with positive slopes in the mass-size relation, so that more massive galaxies are more extended than less massive systems at fixed GLYPH&lt;21&gt; . When all galaxies are analysed together, this trend remains. For the range of masses we analyse, E-SDGs tend to be more compact than E-SDGs at a given dynamical mass, with the larger di GLYPH&lt;11&gt; erences for the most extreme values of GLYPH&lt;21&gt; (Graham et al. 2018).

The lower panels of Fig. 4 displays GLYPH&lt;21&gt; GLYPH&lt;0&gt; " as a function of the stellar mass-weighted ages. As can be seen the tracks of constant age are flat for " &lt; GLYPH&lt;24&gt; 0 3, indicating that for more oblate sys-: tem, stellar age is found to be a good proxy of GLYPH&lt;21&gt; and vice-versa. For both morphological types, older SPs are found in galaxies with lower GLYPH&lt;21&gt; . This is consistent with these systems having experienced a larger number of (dry) mergers while those with larger GLYPH&lt;21&gt; values have progressively younger SPs, consistently having with more extended star formation histories. We note that galaxies with " &gt; GLYPH&lt;24&gt; 0 3 seem to be younger compared to those with : smaller ellipticity at fixed GLYPH&lt;21&gt; . However, no clean proxy for GLYPH&lt;21&gt; is found in these galaxies.

Following previous works, if we take GLYPH&lt;21&gt; GLYPH&lt;24&gt; 0 2 as a suitable : value to separate fast and slow rotators (based on the GLYPH&lt;21&gt; GLYPH&lt;0&gt; GLYPH&lt;15&gt; distribution of EAGLE galaxies presented; see Fig. 3 in Lagos et al. 2018a), we found that 55 per cent of the E-SDGs are slow rotators populated by old stars on average, while only 8 per cent of the E-DDGs are classified as slow rotators. If instead we use the definition adopted by Graham et al. (2018), the fractions are GLYPH&lt;24&gt; 28 per cent and GLYPH&lt;24&gt; 2 per cent, respectively. The fraction of fast rotator in E-SDGs is GLYPH&lt;24&gt; 71 per cent, in agreement with the results reported by Lagos et al. (2018a). However this fraction is lower than that found by Emsellem et al. (2011) for ATLAS 3D ( GLYPH&lt;24&gt; 86 per cent) though this could be caused by ATLAS 3D selecting (by construction) ETGs only. Slow rotators in the EAGLE simulations are detected mostly in dispersion-dominated systems with old SPs. They define the second peak in the GLYPH&lt;21&gt; histogram shown in the lower panel of Fig. 4. This histogram shows a bimodal distribution with the first peak populated by disc galaxies (with young SPs) while the second one, by slow-rotator galaxies (with old SPs). The EAGLE galaxies are found to reproduce well the MaNGA distribution reported by Graham et al. (2018).

## 4. Conclusions

Using a galaxy sample of simulated galaxies selected from the large (100 Mpc) 3 simulation of the EAGLE Project, we performed a statistical analysis of the dependence of the mass-size

Fig. 1. Velocity dispersion ( GLYPH&lt;27&gt; e ) LOESS-smoothed distributions on the stellar mass-size plane for E-SDGs (left panel), E-DDGs (middle panel) and all galaxies (right panel) in the EAGLE simulation at z = 0. The median relations for the EAGLE galaxies (pink rhombus) are also shown. For comparison the median relations for passive (left panel), active (middle panel) and all galaxies together (right panel) reported by van der Wel et al. (2014) are also included in black rhombus. Error bars correspond to the 16 and 84 percentiles for both simulated and observed data.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00416v1-figure-c4846beb-f57e-48d7-8781-9c69b70dbf17.jpg)

Fig. 2. Stellar mass-weighted average age (lower panels) LOESS-smoothed distributions on the dynamical mass-size plane for E-SDGs (left panel), E-DDGs (middle panel) and all galaxies (right panel) in the EAGLE simulation at z = 0. The dashed lines show the predicted distributions for systems with constant GLYPH&lt;27&gt; e at 100, 125, 150, 200, 225, 250 and 300 km s GLYPH&lt;0&gt; 1 (from left to right).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00416v1-figure-22660093-f24a-4d59-9787-24147cb44c99.jpg)

plane on a variety of galaxy properties. Our main results can be summarised as follow:

- -Early and late-type galaxies in the EAGLE simulation, so called E-SDGs and E-DDGs, are found to follow stellar mass-size planes in global agreement with observations, principally for galaxies with stellar masses of GLYPH&lt;24&gt; 10 10 3 : M . GLYPH&lt;12&gt; Smaller and larger mass galaxies show discrepancies but are within observed ranges reported by van der Wel et al. (2014). Larger-size galaxies have lower GLYPH&lt;27&gt; e values at fixed mass for simulated late-type galaxies. In the case of the simulated early-type systems, stellar mass is found to be a good proxy of GLYPH&lt;27&gt; e .
- -Galaxies dominated by old stars are more compact at fixed dynamical mass. Li2018 report GLYPH&lt;27&gt; e to be a good proxy for the stellar age. We find this to be on average valid for EAGLE galaxies. However, mean stellar ages show a slightly weaker trend than GLYPH&lt;27&gt; e for massive disc galaxies which is not observed in Li2018.
- -More extended discs have steeper negative age gradients while early-type galaxies show weak or more positive age

gradients. The latter trend seems to be in agreement with the slight outside-in assembly histories of low-mass, early-type galaxies reported by Rosito et al. (2018).

- -Chemical abundances are found to have a shallower dependence on stellar mass or GLYPH&lt;27&gt; e than reported by observations (see Appendix A). Reproducing the observed trends with metallicity seem to be a major challenge for cosmological simulations.
- -Overall, we find that age can be used as a proxy for GLYPH&lt;21&gt; for more oblate systems as reported by van de Sande et al. (2019). Because the adopted morphological classification is based on the level of rotation in galaxies, this implies that morphology is a good proxy for age. However, " does not show a clear trend with stellar age in the EAGLE galaxies.
- -The distributions of EAGLE galaxies in the ( GLYPH&lt;21&gt;; " ) plane shows a bimodal distribution of GLYPH&lt;21&gt; where the extended, disc systems with young stars populate the first peak while the second one is dominated by slow rotator galaxies with old SPs.

Fig. 3. Age (upper panels) and [O H] (lower panels) gradients LOESS-smoothed distributions on the mass-size plane for E-SDGs (left panels), / E-DDGs (middle panels) and all galaxies (right panels) in EAGLE simulation at z = 0. Dashed lines are the same as given in Fig. 2.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00416v1-figure-ee212421-3b8e-46ad-a1e2-a8fced405bf5.jpg)

175

Fig. 4. Upper panel: Spin parameter ( GLYPH&lt;21&gt; ) distribution on the mass-size plane for E-SDGs (left panel), E-DDGs (middle panel) and all galaxies (right panel) in the EAGLE simulation. Lines are the same as given in Fig. 2. Lower panel: GLYPH&lt;21&gt; GLYPH&lt;0&gt; GLYPH&lt;15&gt; plane coloured by stellar-mass weighted ages. We depict the slow rotators region according to the definition in Graham et al. (2018) (black lines). Lower left panel: Histogram of GLYPH&lt;21&gt; for the complete sample.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00416v1-figure-cc946444-23c8-40e0-9bd3-ff32a13c8ece.jpg)

These results suggest that the processes regulating star formation activity (and hence the age distributions) might also affect the redistribution of angular momentum of stars and gas, and hence, galaxy morphology. Despite the many successes of the current generation of cosmological hydrodynamical simulations, some important shortcomings remain. In the area of stellar populations, we find that some important ones are related to the metallicity profiles and age-morphology relation in galaxies. Further work on the multi-phase nature of the ISM in simulations could naturally lead to more realistic internal kinematics and profiles of galaxies and hence is key to explore.

Acknowledgements. We thank J. Trayford and J. van de Sande for their valuable comments to this work. PBT acknowledges partial funding by Fondecyt Regular 2015 - 1150334 and Proyecto Interno UNAB 2019. This project has received funding from the European Union Horizon 2020 Research and Innovation Programme under the Marie Sklodowska-Curie grant agreement No 734374 and the GALNET network funded by Conicyt. CL is funded by the ARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), through project number CE170100013. This work used the DiRAC Data Centric system at Durham University, operated by the Institute for Computational Cosmology on behalf of the STFC DiRAC HPC Facility (www.dirac.ac.uk). This equipment was funded by BIS National E-infrastructure capital grant ST K00042X 1, / / STFC capital grants ST H008519 1 and ST K00087X 1, STFC DiRAC Opera-/ / / / tions grant ST K003267 1 and Durham University. DiRAC is part of the National / / E-Infrastructure. We acknowledge PRACE for awarding us access to the Curie machine based in France at TGCC, CEA, Bruyeres-le-Chatel. This work used RAGNAR cluster funded by Fondecyt Regular 2015 - 1150334 at Universidad Andres Bello.

## References

Bois, M., Emsellem, E., Bournaud, F., et al. 2011, MNRAS, 416, 1654 Bryant, J. J., Owers, M. S., Robotham, A. S. G., et al. 2015, MNRAS, 447, 2857 Bundy, K., Bershady, M. A., Law, D. R., et al. 2015, ApJ, 798, 7 Cappellari, M., McDermid, R. M., Alatalo, K., et al. 2013, MNRAS, 432, 1862 Clauwens, B., Schaye, J., Franx, M., &amp; Bower, R. G. 2018, MNRAS, 478, 3994 Cleveland, W. S. &amp; Devlin, S. J. 1988, Journal of the American Statistical Association, 83, 596

Crain, R. A., Schaye, J., Bower, R. G., et al. 2015, MNRAS, 450, 1937 De Rossi, M. E., Bower, R. G., Font, A. S., Schaye, J., &amp; Theuns, T. 2017,

MNRAS, 472, 3354 Emsellem, E., Cappellari, M., Krajnovi´ c, D., et al. 2011, MNRAS, 414, 888 Fall, S. M. &amp; Romanowsky, A. J. 2013, ApJ, 769, L26 Frenk, C. S., Theuns, T., Trayford, J. W., Correa, C., &amp; Schaye, J. 2018, MNRAS, 483, 744

Graham, M. T., Cappellari, M., Li, H., et al. 2018, MNRAS, 477, 4711 Lagos, C. d. P., Schaye, J., Bahé, Y., et al. 2018a, MNRAS, 476, 4327 Lagos, C. d. P., Stevens, A. R. H., Bower, R. G., et al. 2018b, MNRAS, 473, 4956

Lagos, C. d. P., Theuns, T., Stevens, A. R. H., et al. 2017, MNRAS, 464, 3850 Li, H., Mao, S., Cappellari, M., et al. 2018, MNRAS, 476, 1765 Lia, C., Portinari, L., &amp; Carraro, G. 2002, MNRAS, 330, 821

- McAlpine, S., Helly, J. C., Schaller, M., et al. 2016, Astronomy and Computing, 15, 72
- Mosconi, M. B., Tissera, P. B., Lambas, D. G., &amp; Cora, S. A. 2001, MNRAS, 325, 34

Naab, T. 2013, in IAU Symposium, Vol. 295, The Intriguing Life of Massive Galaxies, ed. D. Thomas, A. Pasquali, &amp; I. Ferreras, 340-349

Naab, T., Oser, L., Emsellem, E., et al. 2014, MNRAS, 444, 3357

Pedrosa, S. E. &amp; Tissera, P. B. 2015, A&amp;A, 584, A43

- Perez, J., Valenzuela, O., Tissera, P. B., &amp; Michel-Dansac, L. 2013, MNRAS, 436, 259

Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2014a, A&amp;A, 571, A1 Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2014b, A&amp;A, 571, A16 Rosas-Guevara, Y., Bower, R. G., Schaye, J., et al. 2016, MNRAS, 462, 190 Rosito, M. S., Tissera, P. B., Pedrosa, S. E., &amp; Rosas-Guevara, Y. 2018, ArXiv e-prints [ arXiv:1811.11062 ]

Ryden, B. S., Forbes, D. A., &amp; Terlevich, A. I. 2001, MNRAS, 326, 1141 Sánchez, S. F., Kennicutt, R. C., Gil de Paz, A., et al. 2012, A&amp;A, 538, A8 Schaye, J., Crain, R. A., Bower, R. G., et al. 2015, MNRAS, 446, 521 Scott, N., Brough, S., Croom, S. M., et al. 2017, MNRAS, 472, 2833 Taylor, P. &amp; Kobayashi, C. 2017, MNRAS, 471, 3856

- Tissera, P. B., Machado, R. E. G., Sanchez-Blazquez, P., et al. 2016, A&amp;A, 592, A93

Tissera, P. B., Rosas-Guevara, Y., Bower, R. G., et al. 2019, MNRAS, 482, 2208 Tissera, P. B., White, S. D. M., &amp; Scannapieco, C. 2012, MNRAS, 420, 255 van de Sande, J., Lagos, C. D. P., Welker, C., et al. 2019, MNRAS, 484, 869 van de Sande, J., Scott, N., Bland-Hawthorn, J., et al. 2018, Nature Astronomy, 2, 483

van der Wel, A., Franx, M., van Dokkum, P. G., et al. 2014, ApJ, 788, 28

Article number, page 6 of 10

## Appendix A: Supplementary information on EAGLE galaxies

Appendix A.1: The mass-size plane as a function of metallicity

As a function of chemical abundances, we find very weak trends either using (O H) or [O Fe] as shown in Fig. A.1. This is con-/ / sistent with previous results (De Rossi et al. 2017; Tissera et al. 2019). There are also quite large dispersions in the values which make di GLYPH&lt;14&gt; cult disentangle clear trends (see also Tissera &amp; Solar in preparation for a discussion on the metallicity dispersion). The [O Fe] distribution shows a slight trend for massive E-SDGs / to have higher values. This is expected since massive E-SDGs would tend to form stars in short and strong starbursts (Rosito et al. 2018), and hence the chemical abundances are mainly determined by SNe II nucleosynthesis (Scott et al. 2017). It is also interesting to note that, at a given stellar mass, there is a trend for more massive galaxies to have lower [O Fe]. Nevertheless, / we point out for both type of galaxies no systematic dependencies of chemical abundances with mass or GLYPH&lt;27&gt; e are detected.

## Appendix B: Real distributions

In this section, we include the collection of figures shown in the main paper and in Appendix A, using the real distributions instead of LOESS-smoothed ones. This set of figures are given only for comparison.

Fig. A.1. Median (O H) / + 12 (upper panels) and [O Fe] (lower panels) LOESS-smoothed distributions on the mass-size plane for E-SDGs (left / panel), E-DDGs (middle panel) and all galaxies (right panel) in the EAGLE sample. Lines are the same as given in Fig. 2.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00416v1-figure-ab92a796-0a30-43b0-96e0-06b52a88e917.jpg)

Fig. B.1. Velocity dispersion ( GLYPH&lt;27&gt; e ) LOESS-smoothed distributions on the stellar mass-size plane for E-SDGs (left panel), E-DDGs (middle panel) and all galaxies (right panel) in the EAGLE simulation at z = 0. The median relations for the EAGLE galaxies (pink rhombus) are also shown. For comparison the median relations for passive (left panel), active (middle panel) and all galaxies together (right panel) reported by van der Wel et al. (2014) are also included in black rhombus Error bars correspond to the 16 and 84 percentiles for both simulated and observed data. See Fig. 1 for the smoothed version.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00416v1-figure-3c8c6f03-7803-4393-83b1-7cb5da5db41f.jpg)

Fig. B.2. Stellar mass-weighted average age (lower panels) LOESS-smoothed distributions on the dynamical mass-size plane for E-SDGs (left panel), E-DDGs (middle panel) and all galaxies (right panel) in the EAGLE simulation at z = 0. The dashed lines show the predicted distributions for systems with constant GLYPH&lt;27&gt; e at 100, 125, 150, 200, 225, 250 and 300 km s GLYPH&lt;0&gt; 1 (from left to right). See Fig. 2 for the smoothed version.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00416v1-figure-9e9ab595-be11-48a6-ad18-79db2fc62760.jpg)

Fig. B.3. Age (upper panels) and [O H] (lower panels) gradients LOESS-smoothed distributions on the mass-size plane for E-SDGs (left panels), / E-DDGs (middle panels) and all galaxies (right panels) in EAGLE simulation at z = 0. Dashed lines are the same as given in Fig. 2. See Fig. 3 for the smoothed version.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00416v1-figure-859d9feb-1f97-440e-af4b-c6bee01327b7.jpg)

Fig. B.4. Upper panel: Spin parameter ( GLYPH&lt;21&gt; ) distribution on the mass-size plane for E-SDGs (left panel), E-DDGs (middle panel) and all galaxies (right panel) in the EAGLE simulation. Lines are the same as given in Fig. 2. Lower panel: GLYPH&lt;21&gt; GLYPH&lt;0&gt; GLYPH&lt;15&gt; plane coloured by stellar-mass weighted ages. We depict the slow rotators region according to the definition in Graham et al. (2018) (black lines). Lower left panel: Histogram of GLYPH&lt;21&gt; for the complete sample. See Fig. 4 for the smoothed version.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00416v1-figure-9756e9f5-ff88-46a0-af7e-81fa2282533c.jpg)

Fig. B.5. Median (O H) / + 12 (upper panels) and [O Fe] (lower panels) distributions on the mass-size plane for E-SDGs (left panel), E-DDGs / (middle panel) and all galaxies (right panel) in the EAGLE sample. Lines are as given in Fig. 1. See Fig. A.1 for the smoothed version.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00416v1-figure-caf295f2-80a2-4aac-900b-cc1d44e89f5a.jpg)
María Sol Rosito Patricia Tissera Susana Pedrosa Claudia Lagos
The mass-size plane of EAGLE galaxies
Current observational results show that both late-and-early-type galaxies
follow tight mass-size planes, on which physical properties such as age,
velocity dispersion and metallicities correlate with the scatter on the plane.
We study the mass-size plane of galaxies in cosmological hydrodynamical
simulations, as a function of velocity dispersion, age, chemical abundances,
ellipticity and spin parameters with the aim at assessing to what extent the
current cosmological paradigm can reproduce these observations and provide a
physical interpretation of them. We select a sample of well-resolved galaxies
from the (100 Mpc)^3 simulation of the EAGLE Project. This sample is composed
by 508 spheroid-dominated galaxies and 1213 disc-dominated galaxies. The
distributions of velocity dispersion, age, metallicity indicators and gradients
and spin parameters across the mass-size plane are analysed. Furthermore, we
study the relation between shape and kinematic parameters. The results are
compared with observations. The mass-weighted ages of the EAGLE galaxies are
found to vary along lines of constant velocity dispersion on the mass-size
plane, except for galaxies with velocity dispersion larger than aprox 150 km
s^(-1) . Negative age gradients tend to be found in extended disc galaxies in
agreement with observations. However, the age distributions of early-type
galaxies show a larger fraction with inverted radial profiles. The distribution
of metallicity gradients does not show any clear dependence on this plane.
Galaxies with similar spin parameters ({\lambda}) display larger sizes as their
dynamical masses increase. Stellar-weighted ages are found to be good proxies
for {\lambda} in galaxies with low ellipticity ({\epsilon}). Abridged