arxiv: v1 [astro-ph.ga] 14 Jan 2014

Transcription

1 Mon. Not. R. Astron. Soc., 1 18 (213) Printed 2 June 214 (MN LATEX style file v2.2) The ATLAS 3D project XXVI. H I discs in real and simulated fast and slow rotators arxiv: v1 [astro-ph.ga] 14 Jan 214 Paolo Serra, 1,2 Ludwig Oser, 3,4 Davor Krajnović, 5 Thorsten Naab, 3 Tom Oosterloo, 1,6 Raffaella Morganti, 1,6 Michele Cappellari, 7 Eric Emsellem, 8,9 Lisa M. Young, 1 Leo Blitz, 11 Timothy A. Davis, 8 Pierre-Alain Duc, 12 Michaela Hirschmann, 13 Anne-Marie Weijmans, 14,15 Katherine Alatalo, 11,16 Estelle Bayet, 7 Maxime Bois, 8,9 Frédéric Bournaud, 12 Martin Bureau, 7 Roger L. Davies, 7 P. T. de Zeeuw, 8,17 Sadegh Khochfar, 18,19 Harald Kuntschner, 8 Pierre-Yves Lablanche, 8,9 Richard M. McDermid, 2 Marc Sarzi, 21 and Nicholas Scott 22,23 1 Netherlands Institute for Radio Astronomy (ASTRON), Postbus 2, 799 AA Dwingeloo, The Netherlands 2 CSIRO Astronomy and Space Science, Australia Telescope National Facility, PO Box 76, Epping, NSW 171, Australia 3 Max-Planck-Institut für Astrophysik, Karl-Schwarzschild-Str. 1, Garching, Germany 4 Department of Astronomy, Columbia University, New York, NY 127, USA 5 Leibniz-Institut für Astrophysik Potsdam (AIP), An der Sternwarte 16, D Potsdam, Germany 6 Kapteyn Astronomical Institute, University of Groningen, Postbus 8, 97 AV Groningen, The Netherlands 7 Sub-Dept. of Astrophysics, Dept. of Physics, University of Oxford, Denys Wilkinson Building, Keble Road, Oxford, OX1 3RH, UK 8 European Southern Observatory, Karl-Schwarzschild-Str. 2, Garching, Germany 9 Université Lyon 1, Observatoire de Lyon, Centre de Recherche Astrophysique de Lyon and Ecole Normale Supérieure de Lyon, 9 avenue Charles André, F-6923 Saint-Genis Laval, France 1 Physics Department, New Mexico Institute of Mining and Technology, Socorro, NM 8781, USA 11 Department of Astronomy, Campbell Hall, University of California, Berkeley, CA 9472, USA 12 Laboratoire AIM Paris-Saclay, CEA/IRFU/SAp CNRS Université Paris Diderot, Gif-sur-Yvette Cedex, France 13 INAF - Astronomical Observatory of Trieste, via G.B. Tiepolo 11, I Trieste, Italy 14 Dunlap Institute for Astronomy & Astrophysics, University of Toronto, 5 St. George Street, Toronto, ON M5S 3H4, Canada 15 School of Physics and Astronomy, University of St Andrews, North Haugh, St Andrews KY16 9SS, UK 16 Infrared Processing and Analysis Center, California Institute of Technology, Pasadena, California 91125, USA 17 Sterrewacht Leiden, Leiden University, Postbus 9513, 23 RA Leiden, the Netherlands 18 Max-Planck Institut für extraterrestrische Physik, PO Box 1312, D Garching, Germany 19 Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh EH9 3HJ, UK 2 Gemini Observatory, Northern Operations Centre, 67 N. A ohoku Place, Hilo, HI 9672, USA 21 Centre for Astrophysics Research, University of Hertfordshire, Hatfield, Herts AL1 9AB, UK 22 Sydney Institute for Astronomy (SIfA), School of Physics, The University of Sydney, NSW 26, Australia 23 ARC Centre of Excellence for All-sky Astrophysics (CAASTRO), 44 Rosehill Street, Redfern, NSW 216, Australia Accepted... Received...; in original form... c 213 RAS

2 2 Paolo Serra et al. ABSTRACT One quarter of all early-type galaxies (ETGs) outside Virgo host a disc/ring of H I with size from a few to tens of kpc and mass up to 1 9 M. Here we investigate whether this H I is related to the presence of a stellar disc within the host making use of the classification of ETGs in fast and slow rotators (FR/SR). We find a large diversity of H I masses and morphologies within both families. Surprisingly, SRs are detected as often, host as much H I and have a similar rate of H I discs/rings as FRs. Accretion of H I is therefore not always linked to the growth of an inner stellar disc. The weak relation between H I and stellar disc is confirmed by their frequent kinematical misalignment in FRs, including cases of polar and counterrotating gas. In SRs the H I is usually polar. This complex picture highlights a diversity of ETG formation histories which may be lost in the relative simplicity of their inner structure and emerges when studying their outer regions. We find that ΛCDM hydrodynamical simulations have difficulties reproducing the H I properties of ETGs. The gas discs formed in simulations are either too massive or too small depending on the star formation feedback implementation. Kinematical misalignments match the observations only qualitatively. The main point of conflict is that nearly all simulated FRs and a large fraction of all simulated SRs host corotating H I. This establishes the H I properties of ETGs as a novel challenge to simulations. Key words: ISM: kinematics and dynamics galaxies: elliptical and lenticular, cd galaxies: evolution galaxies: formation galaxies: kinematics and dynamics. 1 INTRODUCTION The observation of both neutral hydrogen (H I) and molecular gas (H 2, traced by CO) has shown that early-type galaxies (E/Ss, hereafter ETGs), while on average gas-poorer than spirals, can in some cases host a significant mass of cold gas (for early H I studies see, e.g., Gallagher, Faber & Balick 1975; Knapp et al. 1977; Bieging & Biermann 1977; Knapp, Turner & Cunniffe 1985; Wardle & Knapp 1986; for CO see Wiklind & Henkel 1989; Lees et al. 1991). What has remained unclear until recently is whether gas-rich ETGs are peculiar, rare objects or represent a significant fraction of the overall population (e.g., van Gorkom & Schiminovich 1997). This question has now been answered by large H I and H 2 surveys of morphologically-selected ETGs carried out as part of the ATLAS 3D project (Cappellari et al. 211a, hereafter Paper I). In Serra et al. (212, hereafter Paper XIII) we show that 4 per cent of all ETGs outside the Virgo cluster and 1 per cent inside Virgo host H I down to a mass of 1 7 M and a column density of cm 2, placing recent, previous results (Morganti et al. 26; di Serego Alighieri et al. 27; Grossi et al. 29; Oosterloo et al. 21) on a strong statistical basis. We find that in 2/3 of all detections (and 1/4 of all ETGs outside Virgo) the H I is distributed in a low-column-density disc or ring with size from a few to tens of kpc, and mass from 1 7 to more than 1 9 M (see also van Driel & van Woerden 1991; Oosterloo et al. 27). The remaining detected ETGs host H I on an unsettled configuration which usually reveals recent dynamical interaction within a group. Concerning molecular gas, Young et al. (211, hereafter Paper IV) detect CO in 22 per cent of all ETGs independent of galaxy luminosity or environment, with M(H 2 ) between 1 7 M (the typical sensitivity of their data) and 1 9 M (for recent, previous surveys of smaller samples see Combes, Young & Bureau 27; Sage, Welch & Young 27). Follow-up interferometric observations show that this gas is mostly distributed on discs and rings paolo.serra@csiro.au Dunlap Fellow with typical radius 1 kpc, just a fraction of the galaxy effective radius (Alatalo et al. 213; Davis et al. 213). Taken together, these studies establish that about half of all ETGs host between 1 7 and 1 9 M of cold gas (atomic and/or molecular) usually settled within the galaxy potential. Very sensitive observations are required to reach such a detection rate (for example, this is significantly higher than in the recent surveys of massive galaxies by Catinella et al. 21 and Saintonge et al. 211) but when these are performed we find that ETGs with a significant mass of cold gas are not peculiar, rare systems. The obvious implication is that in order to understand the assembly of ETGs as a family we need to understand how they get or retain their gas and what impact this has on the host, on its star formation history and on its structure. Here we focus on the latter aspect. In most previous studies of this kind ETGs are simply divided into ellipticals and lenticulars. When this is done, molecular gas is almost exclusively found in lenticulars (Lees et al. 1991, Welch, Sage & Young 21, Paper IV), consistent with the fact that CO is invariably associated with dust lanes. Furthermore, H I seems to be more abundant in lenticular than in elliptical objects at a typical M(H I) sensitivity of 1 9 M (Roberts 1975; Wardle & Knapp 1986; Bregman, Hogg & Roberts 1992; Sadler 21) although more sensitive observations of smaller samples do not reveal any significant difference (Grossi et al. 29; Oosterloo et al. 21). One complication is that the distinction between ellipticals and lenticulars can be ambiguous and is affected by projection effects (Jorgensen & Franx 1994). A more robust and physically meaningful approach is to classify ETGs on the basis of their internal kinematics (e.g., Davies et al. 1983; Cappellari et al. 27; Emsellem et al. 27). In this context a key result of the ATLAS 3D project is that the vast majority of all ETGs ( 85 per cent) are structurally similar to spirals. They rotate like discs (Krajnović et al. 211, hereafter Paper II), can be described as a family of oblate rotators on the λ R -ɛ diagram (Emsellem et al. 211, hereafter Paper III), and cover the same range of intrinsic flattening and bulge-to-disc ratio of Sa to Sc spirals (Cappellari et al. 211b, 213b, hereafter Paper c 213 RAS, MNRAS, 1 18

3 H I discs in fast and slow rotators 3 VII and Paper XX, respectively; Krajnović et al. 213a; Weijmans et al., submitted; see also Laurikainen et al. 211; Kormendy & Bender 212). We call these ETGs fast rotators (FRs). In practice, FRs include nearly all classical lenticulars plus a large number of galaxies morphologically misclassified as ellipticals because of unfavourable viewing angle (Paper III). The evolutionary history of most FRs must include the growth and survival of a stellar disc following accretion of cold gas from the surrounding medium as well as gas-rich minor and major mergers (Khochfar et al. 211, hereafter Paper VIII; Naab et al. 213, hereafter Paper XXV). The remaining 15 per cent of all ETGs do not exhibit the disc-like stellar rotation observed in FRs and occupy a region of the λ R -ɛ diagram corresponding to intrinsically lower specific angular momentum and flattening (Paper III). We call these galaxies slow rotators (SRs). We find that SRs form a heterogeneous family including systems with no rotation at all, non-disc-like rotation and kinematically distinct cores (KDCs). Massive SRs are thought to grow via (a small number of) major mergers at z > 1, which lower the system angular momentum and may create a KDC, followed by many minor mergers with low gas fraction, which preserve the KDC and make the stellar body rounder (Naab, Johansson & Ostriker 29; Bois et al. 211, hereafter Paper VI; Paper VIII; Paper XXV). This sequence of events is most likely to occur in a dense environment, in agreement with the increased fraction of SRs in the centre of the Virgo cluster (Paper VII). Within this picture one may naively think that SRs should lack a significant cold-gas phase, unlike FRs. Indeed, Paper IV shows that the CO detection rate is much lower in SRs than in FRs. However, the situation is not necessarily so simple. For example, the detection of kinematically-misaligned ionized gas (Davis et al. 211, hereafter Paper X), cold dust (Smith et al. 212) and stellar nuclear cusps (Krajnović et al. 213b, hereafter Paper XXIII) in SRs suggests that they too can accrete some gas during their life; and simulations show that at least some SRs may have formed in gas-rich mergers (Paper VI; Paper XXV). Furthermore, the mere fact that FRs are found at basically all environment densities (Paper VII) implies that there must be a large variety of evolutionary histories leading to their formation, and therefore potentially a large variety of resulting cold-gas properties. All that is required to mark the visual difference between a FR and a spiral at z = is that the former should have lower starformation-rate surface density on the stellar disc. In other words, they should be quenched owing to a lower M gas /M and/or starformation efficiency within the stellar body. Observations indicate that such quenching may be related to the growth of the stellar bulge: the cold gas content of the stellar body decreases as the bulge becomes more dominant (Catinella et al. 21; Saintonge et al. 211; Paper XX) and, additionally, the bulge may stabilize the remaining gas against star formation (Kawata, Cen & Ho 27; Martig et al. 29; Saintonge et al. 212; Martig et al. 213). Yet it is still possible for a FR to host a large mass of cold gas well outside the stellar body and at low surface density (e.g., Morganti et al. 26; Oosterloo et al. 27; Paper XIII). Such gas distributions are indicative that there are many ways of making a FR and that this information, lost in the apparent homogeneity of their shape and stellar kinematics, may emerge when observing their cold-gas phase. Furthermore, these tenuous gas systems are only weakly bound to the host and hence are easily perturbed. Their actual morphology is therefore a sensitive probe of the past dynamical history of the host, and indeed in Paper XIII we show that the H I morphology varies systematically with the density of the environment. In this spirit we study here the H I properties of ETGs as a function of their structure within the context of the ATLAS 3D project. The few previous attempts in this direction reveal no clear relation between H I and ETG structure but are based on significantly smaller samples (Morganti et al. 26; Oosterloo et al. 21). Here we aim at clarifying the situation by studying the large ATLAS 3D H I sample. We describe the sample and the H I data in Sec. 2, discuss H I detection rate, mass and morphology as a function of galaxy structure in Sec. 3, analyse the kinematical misalignment between H I and stars in Sec. 4, perform for the first time a comparison to the predictions of ΛCDM hydrodynamical simulations in Sec. 5, and summarize the results in Sec SAMPLE AND H I SURVEY The ATLAS 3D H I sample and the radio observations which this study is based on are described in detail in Paper XIII. Here we provide a short summary. We study a subset of the full ATLAS 3D sample of ETGs (Paper I; distance < 42 Mpc, M K < 21.5). This subset includes all galaxies with declination δ 1 except the four objects closest to Virgo A, resulting in a volume-limited sample of 166 ETGs (127 outside the Virgo cluster). We obtain Westerbork Synthesis Radio Telescope H I data for most of them either as part of the ATLAS 3D project or from previous projects (Morganti et al. 26; Józsa et al. 29; Oosterloo et al. 21). We make use of Very Large Array data for 1 galaxy with H I (Chung et al. 29) and of Arecibo spectra for 2 galaxies (all undetected; Giovanelli et al. 25). As mentioned in Sec. 1, our Westerbork observations allow us to detect H I down to a column density of a few times 1 19 cm 2 and a mass of 1 7 M. We detect H I associated with the target ETG in 53 galaxies (49 outside Virgo). In Paper XIII we classify detected ETGs in the following classes on the basis of the H I morphology: D = large discs (14 objects); R = large rings 1 (1); d = small discs (1); u = unsettled distributions (14); c = gas clouds (5). Gas distributions characterized by ordered rotation (D, R and d classes) are therefore the most common ones among H I-rich ETGs, amounting to 2/3 of all detections and 1/4 of all ETGs outside Virgo. As discussed in Paper XIII, the above H I classes are a way to simplify the observed continuum of gas morphologies, with a number of objects being intermediate (or in transition) between classes. We make use of the H I classes D, R, d, u and c in the rest of this paper. In Paper XIII we show H I images of all detections. Each image is obtained after creating a mask including all emission in the corresponding H I cube. The images shown in Paper XIII are the zero-th moment of the masked cubes. In the present paper we show and analyse for the first time the H I velocity fields obtained as the first moment of the same masked cubes. 3 RELATION BETWEEN H I AND ETG INNER STRUCTURE In this section we investigate whether the H I detection rate, mass and morphology of ETGs are related to their stellar kinematics as 1 In fact, in Paper XIII we group Ds and Rs in a single class. In table 1 of that paper we indicate galaxies where the H I is distributed on a ring. Here we use that table as a basis for the division into D and R objects. c 213 RAS, MNRAS, 1 18

4 4 Paolo Serra et al. Table 1. Number of galaxies in the various H I morphological classes (Paper XIII) as a function of stellar kinematics class (Paper III) and SR/FR classification (Paper II) for the full ATLAS 3D H I sample. Full Stellar kinematics Sample All a b c d e SR FR All Undet Det H I D R d u c Table 2. Number of galaxies in the various H I morphological classes (Paper XIII) as a function of stellar kinematics class (Paper III) and SR/FR classification (Paper II) for non-virgo ETGs in the ATLAS 3D H I sample. Outside Stellar kinematics Virgo All a b c d e SR FR All Undet Det H I D R d u c well as position on the λ R -ɛ diagram 2 and mass-size plane. In particular, we make use of the kinematical classification presented in Paper II. In that paper ETGs are divided in the following classes based on a quantitative analysis of their stellar velocity field 3 : a = no rotation; b = some rotation but with a velocity field not consistent with that of a rotating disc; c = KDC; d = two counterrotating stellar discs; e = velocity field consistent with that of a rotating disc. These classes are important to understand the classification of ETGs into FRs and SRs discussed in Sec. 1. In Paper III we define FRs and SRs as galaxies above and below the line λ R =.31 ɛ on the λ R -ɛ diagram, respectively. Fast rotators and e ETGs are essentially the same family and represent 85 per cent of all objects. On the other hand, the SR family is very heterogeneous and includes galaxies in all other groups (a to d ). 3.1 Fast versus slow rotators Figure 1 shows all galaxies on the λ R -ɛ diagram. In this figure we indicate galaxies with different stellar kinematics and H I morphology using different colours and markers, respectively. To first order Fig. 1 confirms the results of Morganti et al. (26) and Oosterloo et al. (21): the relation between H I properties and galaxy structure is not a simple one. We detect H I in galaxies at all locations on the λ R -ɛ diagram where there are ETGs, both FRs and SRs. Many of the H I discs and rings (D, R and d) are hosted by flat, nearly edge-on FRs (topright of the diagram) and their lower-inclination counterparts (at lower values of λ R and ɛ). This is reasonable as these galaxies have a pronounced stellar disc component and therefore must have experienced a considerable level of dissipation during their assembly. Their current H I discs/rings may be a remnant of that process. The more detailed analysis of Sec. 4 shows that this is indeed the case for about half of all gas-rich FRs, but also that in the remaining half the gas accretion episode traced by H I does not seem to be related to the formation of the current stellar disc. What might be more surprising is to find H I also among SRs (a, b and c objects only) and that in some of them the detected H I is distributed on a large disc or ring. This indicates that SRs too can experience gas accretion and dissipative processes during their 2 We adopt the λ R and ɛ measurements obtained within a 1 R e aperture in Paper III. Note that not all galaxies have their stellar kinematics measured out to 1 R e, as discussed in that paper. 3 Note the addition of a subscript to the kinematical classes name to avoid confusion with the H I classes defined in Sec. 2. formation, and gas-poor evolutionary paths are not the only way to make a galaxy with this internal structure. A similar conclusion is reached by Paper X based on the misalignment between ionized gas and stellar kinematics in these galaxies within 1 R e, by Smith et al. (212) based on the frequent detection of cold dust in SRs, by Paper XXIII based on the existence of SRs with a cuspy nuclear light profile, and by both idealized and ΛCDM hydrodynamical simulations (Paper VI; Paper XXV). Additional observations such as very deep optical imaging may provide further clues about the formation of SRs in presence of a cold gas component (Duc et al. 211). These qualitative results are illustrated in more quantitative terms in Table 1 for the full sample and in Table 2 for non-virgo galaxies only. Within the limited statistics available for SRs (and for classes a to d ), their H I detection rate is the same as that of FRs (and e objects). Namely, we detect H I in 41 ± 14 per cent of all SRs and 31 ± 5 per cent of all FRs, while outside Virgo the detection rate becomes 5 ± 18 per cent for SRs and 37 ± 6 per cent for FRs (error bars assume a binomial distribution). The fraction of galaxies hosting a rotating H I distribution (H I classes D, R and d) is 18 ± 9 per cent for SRs and 21 ± 4 per cent for FRs over the full sample, while outside Virgo it becomes 25 ± 13 per cent for SRs and 25 ± 5 per cent for FRs. The above paragraph describes the H I properties (detection rate and morphology) as a function of stellar kinematics. A complementary approach is to examine the stellar kinematics as a function of H I properties, but also in this case Tables 1 and 2 show that there is no clear trend. Within each H I class (including that formed by all undetected ETGs) the fraction of SRs and FRs is consistent with the global values of 15 and 85 per cent, respectively. The same is true considering all galaxies with a settled H I distribution together (D, R and d). We also compare the λ R /(.31 ɛ) distribution of D + R + d ETGs to that of undetected + c objects 4 and find that they are indistinguishable according to a two-sample KS test the probability that the two samples are not drawn from the same parent distribution is 78 per cent, insufficient to claim a significant difference. The only possible exception to this lack of trends is that small H I discs (d) include no SRs. This is consistent with the fact that H I in these galaxies seems to be strongly linked to the presence of molecular gas and star formation within a stellar disc (Paper XIII). 4 In this context c ETGs can be grouped together with non-detections since M(H I) is very low and the detected gas does not appear obviously associated with the host (Paper XIII). c 213 RAS, MNRAS, 1 18

5 H I discs in fast and slow rotators U6176 N3945 N2824 N3998 N3489 N5582 N3941 N415 N728 N436 N3626 U9519 N2685 N3838 N2764 N4111 N5422 N513 N426 N6798 N471 e regular rotation d 2σ λre N3457 N68 N2859 N332 N3499 N4262 N3182 N423 N3619 U548 N4278 N3193 N5173 N373 N5631 N5198 N5557 U396 N3414 N446 N123 N7465 N3522 N2594 N2768 N5866 N4694 H i class undet ɛ Re D R d u c c KDC b complex rotation a no rotation Figure 1. H I on the λ R -ɛ diagram. The colour scheme represents the kinematical classification of ETGs presented in Paper II. The markers shape represent the H I morphology as in Paper XIII with the small modification described in Sec. 2. The magenta line represents the edge-on view of anisotropic, oblate models with β =.65 ɛ (Cappellari et al. 27), and the black dashed lines represent the models with ɛ =.75 and.85 viewed under all possible inclinations from edge-on (top right) to face-on (bottom left). The green line represents the empirical separation between FRs and SRs λ R =.31 ɛ proposed in Paper III. We indicate the name of all H I detections except that of c galaxies. However, the number of ds is low and only 1 SR would be expected among them based on the global SR fraction. Better statistics would therefore be needed to confirm this result. Another important quantity to analyse is the H I mass. Figure 2 shows the M(H I)/M ratio plotted against the λ R /(.31 ɛ) ratio. Here and in the rest of this paper we adopt for M the value derived by Cappellari et al. (213a, hereafter Paper XV) based on dynamical modelling M JAM in that paper 5. One galaxy, PGC 71531, has no reliable M measurement and is excluded from this analysis. Figure 2 shows that H I-detected SRs are found over a wide 5 The results discussed in this section do not change if we normalize M(H I) to L K. This follows from the fact that the distribution of M /L K ratio is very narrow for ETGs, peaking at 1.2 M /L and with an rms of.4 M /L. This scatter is insufficient to significantly alter a plot like the one showed in Fig. 2, where galaxies spread over 3 orders of magnitude in M(H I)/M. range of M(H I)/M. They seem not to reach the large M(H I)/M values allowed for FRs but the overall low number of SRs makes this difference statistically insignificant. For example, the fraction of galaxies with M(H I)/M.2 is 23 ± 1 per cent within the SR family and 24±4 per cent for FRs. These fractions remain comparable to one another if we consider ETGs outside Virgo only, and these (and the following) results do not change if we adopt a different M(H I)/M threshold. Conversely, we define two ETG subsamples with M(H I)/M above and below.2, including 39 and 126 galaxies respectively. We find that they contain a similar fraction of SRs consistent with the global value of 15 per cent, both considering all galaxies or non-virgo galaxies only. Finally, according to a two-sample KS test the probability that the λ R /(.31 ɛ) distributions of the two subsamples are not drawn from the same parent distribution is just 4 per cent. c 213 RAS, MNRAS, 1 18

6 6 Paolo Serra et al. N7465 U9519 e regular rotation log1 M(Hı)/M N3522 N5631 U396 N5198 N5557 N3414 N446 N2685 N5173 N5582 U548 N6798 N373 N2764 U6176 N123 N423 N68 N2594 N3626 N4694 N513 N3941 N3457 N4262 N4111 N3619 N332 N426 N3838 N4278 N728 N3998 N2859 N436 N3193 N2824 N5422 N3499 N3489 N2768 N3182 N415 N5866 N471 N3945 Hı class undet. D R d u c d 2σ c KDC b complex rotation a no rotation λ Re /(.31 ɛ Re ) Figure 2. M(H I)/M plotted against the ratio λ R /(.31 ɛ), such that SRs and FRs are found to the left and right of the green, vertical line, respectively. Colours and markers are as in Fig H I properties within the FR family Considering now FRs alone, we find no systematic variation of H I properties as a function of position on the λ R -ɛ diagram. Specifically, we compare intrinsically flatter FRs to rounder ones. We do this by considering e galaxies only and splitting them in two groups above and below the line of.75 intrinsic ellipticity, respectively (see Fig. 1; we comment below on alternative ways of defining flat and round galaxies). This results in 47 flat and 89 round e ETGs (39 and 67 outside Virgo, respectively). We detect H I in 32 ± 8 per cent of all flat objects and 3 ± 6 per cent of all round ones. The detection rate of rotating H I distributions (D, R and d objects) is 26 ± 7 and 2 ± 5 per cent, respectively, with both flat and round ETGs exhibiting a mix of D, R and d objects. No significant difference emerges if we consider non-virgo galaxies only. Conversely, the fraction of flat and round e ETGs does not vary significantly as a function of H I class (also when grouping together all D, R and d systems). Within the error bars, we always find 35 per cent of flat objects and 65 per cent of round galaxies. Furthermore, the λ R /λ R,ɛ=.75 distribution of galaxies with an H I disc or ring (D, R and d) is statistically undistinguishable from that of galaxies with no H I (undetected and c). The two-sample KS-test probability that the two distributions are not drawn from the same parent distribution is 38 per cent. We also investigate possible variations of M(H I)/M within the FR family. We find that 19 ± 6 per cent of all flat galaxies have M(H I)/M.2, comparable to the 26 ± 5 per cent found for round objects. Conversely, we define an H I-rich and an H I-poor subsample of e ETGs adopting the same M(H I)/M separation as above. The two samples include 32 and 13 galaxies, respectively. The fraction of flat objects in the two groups is comparable (28 ± 9 and 37 ± 6 per cent, respectively). Finally, a two-sample KS test gives a probability of just 1 per cent that the λ R /λ R,ɛ=.75 distributions of the two subsamples are not drawn from the same parent distribution. The definition of flat and round FRs adopted above relies on a simple de-projection of galaxies shape along model lines on the λ R -ɛ diagram. A more accurate estimate of the intrinsic ellipticity of ETGs is derived in Paper XX based on the modelling presented in Paper XV and on the ellipticity measurement at large radius rec 213 RAS, MNRAS, 1 18

7 H I discs in fast and slow rotators 7 log 1 R e maj (kpc) Hı class undet. D. R d u.2 c log 1 M (M ) Figure 3. H I properties on the mass-size plane. Different H I morphologies are indicated as in previous figures (see legend on the bottom right) while markers are colour-coded according to the M(H I)/M value. The velocity dispersion is constant along the blue dashed lines for galaxies following the virial relation (σ e = 5, 1, 2 and 3 km s 1 from top-left to bottomright, respectively; see Cappellari et al. 26). The red solid line indicates the Zone of Exclusion discussed in Paper XX. ported in Paper II. However, we have verified that using these intrinsic ellipticity values does not change any of the above conclusions. Flat and round FRs keep having equal H I detection rates as well as equal fractions of H I discs and of galaxies with M(H I)/M.2. Similarly, the distribution of FR intrinsic ellipticity does not change as a function of H I mass or morphology. 3.3 Mass-size plane Although the λ R -ɛ diagram is very useful to separate SRs from FRs, Paper XX shows that a number of key ETG properties vary in a much clearer way on the mass-size plane instead. In particular, properties linked to the star formation history of galaxies such as their mass-to-light ratio, optical colour, Hβ absorption and M(H 2 )/M ratio correlate with bulge-to-disc ratio (traced by velocity dispersion) on this plane. We therefore investigate whether the ETG H I properties too exhibit any systematic variation on the mass-size plane. Figure 3 shows all galaxies on the mass size plane of Paper XX. This figure does not show any clear trend other than the decrease of the typical M(H I) value and the prevalence of disturbed H I morphologies at high galaxy mass, already reported in Paper XIII. As with Fig. 2, the main observational result is a very large scatter of H I properties, reinforcing our conclusion that H I does not appear to be linked in any clear way to the inner galaxy structure. 3.4 Diversity of gas accretion histories The lack of a relation between ETG H I properties (detection rate, mass and morphology) and ETG structure demonstrates the importance of cold gas as an independent probe of the evolution of these objects. Neutral hydrogen observations give us clues about the diversity of ETG assembly histories which are otherwise lost in the log1 M(Hı)/M relative simplicity of their inner shape and kinematics. Our H I survey demonstrates that there are very diverse evolutionary paths for FRs, which include the formation of both gas-poor objects (undetected in H I or with very low M(H I)/M ) as well as galaxies with a considerable mass of H I spread over a large area (and therefore with low mean column density; see Paper XIII), but that this diversity is not reflected in the structure of ETGs within the FR family. Furthermore, we find that SRs too can accrete cold gas during their life and that they do so at a similar rate as FRs. It is important to realise that SRs have a fundamentally different internal structure than FRs despite having similar H I properties (in a statistical sense). This means that the accretion of cold gas has not altered their structure, which remains characterized by a relatively round shape and low specific angular momentum within 1 R e. In other words, gas accretion has not resulted in the growth of a significant stellar disc component in the central galaxy region. The same might be true for FRs: the H I discs and rings in FRs do not need to be causally related to the stellar discs revealed by integral-field spectroscopy inside 1 R e, although they might well be in at least some cases. To investigate this aspect in more detail we now turn to the study of the kinematical misalignment between H I and stars in these objects. 4 KINEMATICAL MISALIGNMENT BETWEEN H I AND STARS In order to study the kinematical misalignment between H I and stars we define a sample made of all galaxies for which we can reliably measure the position angle (PA) of the H I kinematical major axis, PA(H I). This includes all galaxies in the D, R and d classes plus two objects hosting H I on an overall unsettled configuration (u) but where some of the gas shows ordered rotation NGC 123 and NGC 426. The total sample studied in this section includes therefore 36 ETGs. 4.1 Polar, corotating, counterrotating and warped H I Figure 4 shows H I images (from Paper XIII), H I velocity fields (shown for the first time here) and stellar velocity fields (already shown in Paper II) of all D, R and u objects in this sample (26 objects in total; the 1 galaxies classified as d are analysed separately below). The figure reveals a diversity of H I-stars misalignments, including corotating, polar and counterrotating gas distributions. In a number of cases the strong variation of PA(H I) with radius indicates the warping of the gas disc. More quantitatively, we derive PA(H I) as a function of radius from the H I velocity field using the KINEMETRY software (Krajnović et al. 26). In order to maximize the radius out to which PA(H I) can be measured we lower the default minimum coverage of an ellipse with data points to 3 per cent. As we are only interested in the PA values this is an acceptable choice. We compare PA(H I) to the PA of the stellar kinematical major axis, PA(stars), which is given in Paper II. Before proceeding we highlight two aspects which should be kept in mind in order to interpret the results of this analysis correctly. First, while PA(H I) is a function of radius and can be measured sometimes out to tens of R e, PA(stars) is an average value measured over an aperture of typically 1 R e (see Paper II for a detailed discussion). Therefore, we are not comparing the H I kinematics to the stellar kinematics at the same radius, but the c 213 RAS, MNRAS, 1 18

8 8 Paolo Serra et al. NGC 123 (u) NGC 2594 (R) NGC 2764 (D) NGC 2859 (R) NGC 3414 (R) NGC 3522 (R) cm cm cm cm cm (1.4) 119 cm (2.3) 119 cm 2 ±21 ±15 ±18 ±15 ±18 ±13 ±13 ±21 ±15 ±18 ±15 ±18 ±13 ±13 NGC 3619 (D) -V NGC 2685 (D) NGC 3626 (D) NGC 3838 (R) NGC 3941 (R) NGC 3945 (R) NGC 3998 (D) NGC 426 (u) cm cm cm cm cm cm cm 2 ±15 ±21 ±17 ±17 ±26 ±24 ±2 ±15 ±21 ±17 ±17 ±26 ±24 ±2 -V/2 +V/2 +V Figure 4. From top to bottom (three panels per galaxy), H I images (1 1 kpc2 ), H I velocity fields (5 5 kpc2 ) and stellar velocity fields (1 1 kpc2 ) of all D and R ETGs, plus the two u s showing some ordered rotation (NGC 123 and NGC 426). The grey scale-bar represents 5 kpc in all figures. The black ellipses overlaid on all velocity fields have semi-major axis equal to 3 Re and ellipticity taken from Paper II. The black arrow overlaid on the H I velocity fields indicates the PA of the stellar kinematical major axis (receding side) from Paper II and scales with the size of the ellipse. In a few cases the stellar kinematical PA could not be determined reliably and no arrow is shown (see the text). For all velocity fields the velocity range is indicated on the top-right. This is kept fixed for the stellar and H I velocity field of a same galaxy in order to enable a comparison of the relative line-of-sight rotation amplitude. The colour scheme of the velocity fields is represented by the horizontal colour bar at the bottom of the figure. H I kinematics (out to large radius) to the inner stellar kinematics. A comparison to the stellar kinematics at larger radius would require deeper optical spectroscopy over a larger field (e.g., Proctor et al. 29; Weijmans et al. 29; Arnold et al. 211, 213) or the study of planetary nebulae (e.g., Romanowsky et al. 23; Coccato et al. 29; Cortesi et al. 213) or globular cluster kinematics (e.g., Kartha et al. 213; Pota et al. 213). A comparison between gas and stellar kinematics at a common, small radius is presented in Paper X using ionized and molecular gas data. Secondly, while PA(stars) is a well defined quantity for FRs, where it is typically constant within the observed field, its in- terpretation may be ambiguous in SRs, where the stellar kinematics is more complex or characterized by very low-level rotation. The SRs for which we can measure PA(H I) are NGC 3414, NGC 3522, NGC 5631 (all KDCs) and UGC 396 (b? ). In the case of NGC 3414 PA(stars) reflects the kinematics of the region outside (and counterrotating relative to) the KDC. To obtain the PA of the KDC one would have to subtract 18 from the value given in Paper II. On the contrary, in NGC 5631 PA(stars) represents the inner (i.e., the KDC s) stellar rotation. In NGC 3522 the error on PA(stars) is very large because of the overall low-level rotation. In this case we calculate the misalignment relative to the photomet- c 213 RAS, MNRAS, 1 18

9 H I discs in fast and slow rotators NGC 436 (R) NGC 423 (D) NGC 4262 (R) NGC 4278 (D) NGC 513 (D) NGC 5173 (D) 9 NGC 5582 (R) cm cm cm (.9) 119cm cm cm cm 2 ±26 ±14 ±18 ±22 ±12 ±7 ±21 ±26 ±14 ±18 ±22 ±12 ±7 ±21 NGC 5631 (D) -V NGC 6798 (D) UGC 396 (D) UGC 6176 (D) UGC 9519 (D) cm cm cm cm cm 2 ±17 ±19 ±11 ±14 ±11 ±17 ±19 ±11 ±14 ±11 -V/2 +V/2 +V Figure 4. Continued ric major axis, which coincides with the (uncertain) kinematical one. Similarly, PA(stars) is unconstrained in UGC 396 and we adopt the photometric PA as a reference. Note that using the photometric rather than kinematical PA would make no difference for NGC 3414, while NGC 5631 is nearly round on the sky and kinematics is the only possible choice. In Fig. A1 we show PA(H I) PA(stars) as a function of radius for all 26 galaxies in this subsample. This figure and Fig. 4 show that 3 per cent of the systems host corotating H I (NGC 2859, NGC 3838, NGC 3945, NGC 426, NGC 436, NGC 423, NGC 5582, UGC 6176), while counterrotating and polar H I distributions account each for 1 per cent of the sample (NGC 3626, NGC 3941 and NGC 6798 in the first group, NGC 2594, NGC 3998 and NGC 4262 in the second). Note however that the observed misalignment of polar systems depends on viewing angle and it is likely that the H I is on a polar configuration in a larger fraction of all ETGs (e.g., NGC 3619, NGC 513 c 213 RAS, MNRAS, 1 18 and UGC 9519). In most of the remaining systems the misalignment angle varies significantly with radius, indicating the presence of a warp. In the two most extreme cases the H I changes from polar to corotating (NGC 2685) and from counter- to corotating (NGC 5173) with increasing radius. The situation is somewhat different for ds. The H I in these objects is typically very faint and extends over only a couple of Westerbork beams. For this reason instead of using KINEMETRY we assume that PA(H I) does not vary with radius and determine its value by direct analysis of the H I cubes. Position-velocity diagrams along the adopted H I kinematical major axes are published in Paper XIII (their Fig. 5). We find that most (8/1) of the ds are corotating. The exceptions are one counterrotator (NGC 332) and one polar system (NGC 3499), as discussed in that paper. Considering all 36 galaxies together we find therefore that in 45 per cent of all ETGs with an H I disc or ring the gas (which can spread out to tens of Re ) is corotating with the stars inside 1 Re.

10 1 Paolo Serra et al. These galaxies are all FRs where the gas accretion episode traced by the H I might be related to the formation of the stellar disc. They are in many respects similar to spiral galaxies as they host cold gas and stars corotating on the same disc the only difference being the absence of bright, star-forming spiral arms owing to the low H I column density (Paper XIII). The remaining 55 per cent of this sample is composed of 1 per cent of counterrotating H I systems, 1 per cent of polar objects, and 35 per cent of ETGs with a complex, warped H I distribution, some of which may in fact be polar. The link between H I disc and stellar disc in these objects is less obvious. These results show that the existence of multiple kinematical sub-components in ETGs, well documented in their inner region for the stellar and ionized gas phases, holds also on much larger scales. Following our results on the fraction of ETGs with a kpcscale KDC ( 7 per cent including counterrotating cases; Paper II) and misaligned ionized-gas kinematics ( 36 per cent; Paper X; see also Bureau & Chung 26) we are now finding misaligned H I distributions on scales from a few to tens of kpc in 12 per cent of all ETGs ( 15 per cent outside the Virgo cluster; for previous studies of smaller samples of ETGs with misaligned H I see, e.g., van Gorkom, Schechter & Kristian 1987; Bertola, Buson & Zeilinger 1992; Bettoni et al. 21). Whether or not associated with the presence of a stellar disc, the H I detected in ETGs provides important clues on the origin of the cold interstellar medium detected in the inner region of these objects in the form of ionized and molecular gas. This is clear in Fig. A1, where we show at R = the data point corresponding to the ionized gas kinematics (when detected), which was already found to be invariably aligned with the CO kinematics (Paper X). The figure confirms previous claims that H I, H 2 and ionized gas share the same kinematics in the region where they overlap (Morganti et al. 26; Paper X). The few exceptions are NGC 123, NGC 3414, NGC 4278, NGC 513 and NGC The first hosts an unsettled gas distribution with very complex gas kinematics, and it is unlikely that the central value of PA(H I) derived from the velocity field is representative of the actual inner gas motion. Furthermore, Morganti et al. (26) show that in both NGC 3414 and NGC 4278 the ionized gas velocity field twists in such a way that the outer ionized gas kinematics is aligned with the inner H I kinematics. The same occurs in NGC Therefore, these three galaxies conform to the general rule too. We do not have a good explanation for the misalignment between H I and ionized gas in NGC 513, which represents the only genuine exception to the rule. Not shown in Fig. A1, d galaxies contribute to strengthening this result. The H I is always corotating with the CO and ionized gas, including the two cases where the H I is misaligned relative to the stars (NGC 332 and NGC 3499). 4.2 Misalignment in fast and slow rotators Figure 5 summarizes our findings by showing all H I misalignments on the M(H I)/M versus λ R /(.31 ɛ) plane of Fig. 2. For each galaxy we adopt here the median misalignment value measured from the radial profile shown in Fig. A1. The median value is a fair representation of the overall misalignment for all galaxies except the heavily warped NGC 2685 and NGC 5173, which we highlight in the figure. On the one hand Fig. 5 shows a significant population of FRs where H I and stars are kinematically aligned (dark blue symbols, with the just mentioned exception of NGC 2685 and NGC 5173). In these galaxies gas and stars are just two different but coupled components of the same disc. On the other hand, as already emphasized above, the H I discs and rings detected in FRs are not always associated with the inner stellar disc, not even for the flattest and fastest rotators (see for example NGC 3626, NGC 513 and UGC 9519). In many FRs the large scale H I kinematics is obviously not a continuation of the inner stellar kinematics to larger radius. On the contrary, counterrotating and polar H I distributions account for exactly 1/4 of all FRs analysed in this section, while another 1/4 consists of warped discs. The detected H I is therefore aligned with the inner stellar disc only in half of the cases. Neutral hydrogen observations show therefore that FRs, while being a dynamically simple and homogeneous family inside 1 R e, are much more complex and diverse in their outer regions, in agreement with results obtained using other kinematical tracers at large radius (e.g., Coccato et al. 29; Pota et al. 213). The other clear result is that in none of the four SRs the H I kinematics is aligned with the stellar kinematics or morphology. In fact, all these systems exhibit a (close to) polar H I distribution, although PA(H I) typically varies with radius (Fig. A1). We also note that the H I morphology of these systems, although disc- or ring-like, appears disturbed (Fig. 4), with some cases intermediate between a D/R and a u classification. The lack of very regular, settled H I discs/rings in SRs and the prevalence of polar systems may be the only difference from the FR family (H I-wise), although even among FRs many gas discs and rings appear disturbed or not perfectly settled (e.g., NGC 2764, NGC 3998, NGC 423, NGC 4278, NGC 513). A larger sample of H I-rich SRs would be needed to confirm this indication. The above results, together with the established weak relation between M(H I) and ETG luminosity (e.g., Knapp, Turner & Cunniffe 1985; Wardle & Knapp 1986; Paper XIII), confirm the severe disconnect between ETG H I properties and their stellar component. They show, as argued in Sec. 3.4, that H I provides clues about the complex evolution of ETGs which would go unnoticed using only other observables. Confirming this complexity, we also find that the H I kinematical misalignment does not depend in any clear way on ETG stellar mass, environment density (quantified in Paper VII) or H I mass. For example among the very H I-rich galaxies in poor environments we find NGC 5582 (corotating R), NGC 6798 (counterrotating D) and UGC 9519 (warped D). Low-mass, isolated galaxies include NGC 2594 (polar R), UGC 6176 (aligned D) and again UGC Among the massive ETGs we find relatively small, aligned H I discs/rings as in NGC 2859 and NGC 436 but also (close to) polar H I as in NGC 3414 and NGC 5631, and the warped D NGC The only exception is that the lowest-m(h I) systems are all ds (Paper XIII) and, as illustrated above, mostly kinematically aligned to the stellar rotation. There is of course some dependence of H I properties on ETG mass and environment across the entire ATLAS 3D sample. This is discussed in detail in Paper XIII, where we show that unsettled H I distributions become more common, and the typical H I mass decreases, as environment density and galaxy mass increase. Here we find that there is no additional trend with environment and stellar mass when studying the kinematical misalignment of settled H I discs. Ever since the work of Knapp, Turner & Cunniffe (1985) this decoupling between stars and H I in ETGs has been attributed to the external origin of the H I. The main point of this interpretation is that most of the present-day gas must be accreted once the formation of the (inner region of the) host ETG is basically complete. The accreted gas could come from the intergalactic medium or be brought in by small gas-rich satellites which individually do not contribute significantly to the stellar mass of the ETG but may bring a significant addition to its interstellar medium. In fact, even c 213 RAS, MNRAS, 1 18

11 H I discs in fast and slow rotators U9519 log1 M(Hı)/M N3522 N5631 U396 N3414 N2685 N5173 N5582 U548 N6798 N2764 U6176 N123 N423 N2594 N3626 N513 N3941 N4262 N3619 N332 N426 N3838 N4278 N3998 N2859 N436 N2824 N5422 N3499 N3489 N Hı-star misalignment (deg) 4 N5866 N3182 N471 N λ Re /(.31 ɛ Re ) Hı class undet. D R d u c 45 Figure 5. H I kinematical misalignment on the M(H I)/M versus λ R /(.31 ɛ) plane. This figure is the same as Fig. 2 but here we highlight only the 36 galaxies for which we are able to study the H I kinematical misalignment (see Sec. 4 for details on this subsample; the remaining galaxies are shown with light grey markers). We colour-code them according to the median misalignment angle calculated over the full radial extent of the H I disc or ring. See Fig. A1 for the variation of the misalignment angle as a function of radius in individual galaxies. Note that NGC 2685 and NGC 5173, the two most heavily warped H I systems in the sample, have a small median misalignment angle which is not representative of their true complexity. For this reason we highlight them here with a red box. when gas is acquired at the same time as the formation of the stellar body, e.g. in gas-rich major mergers, there does not need to be a strong coupling between the stellar body of the remnant and the larger-scale H I distribution. Merging and accretion events are expected to occur in a ΛCDM Universe, where galaxies grow hierarchically, and our results may be seen as a confirmation of these theories. Qualitative statements of this kind have been made by previous authors and, at this stage, do not represent an advancement in our understanding of ETG formation. For this reason we attempt to take a step forward and perform a quantitative comparison to galaxy formation simulations in the next section. 5 COMPARISON TO COSMOLOGICAL SIMULATIONS We compare ETGs in the ATLAS 3D sample to galaxies formed in two sets of hydrodynamical simulations performed by Oser et al. (21) and Hirschmann et al. (213). These are high-resolution re-simulations of 5 dark-matter haloes extracted from a large, ΛCDM, dark-matter-only simulation. In each re-simulation baryonic processes such as hydrogen and helium cooling, star formation and supernova feedback are implemented in the presence of a uniform UV background. The two sets of simulations have different implementations of supernova feedback. In the first set the feedback is purely thermal (hereafter, NoW simulations because they do not include galactic stellar winds). These are the simulations analysed by Paper XXV to study the generic formation paths of SRs and FRs in a cosmological context. Simulations in the second set start from identical initial conditions but include an additional empirical model for galactic outflows driven by winds from starforming regions (Oppenheimer & Davé 26; we will refer to these as W simulations). All details on the galactic wind implementation and a comparison of galaxy properties between the two models are presented in Hirschmann et al. (213). Here we study for the first c 213 RAS, MNRAS, 1 18

12 12 Paolo Serra et al. time the distribution and kinematics of cold gas in the simulations. As we will show, these two sets give profoundly different predictions for galaxies H I properties. We refer to Oser et al. (21) and Hirschmann et al. (213) for a full description of the simulations. For the purpose of this work it is important to first investigate whether the simulated galaxies are a reasonable match to the ETGs we want to compare them to. On the one hand, they cover the same M range (from 1 1 to a few times 1 11 M ). Furthermore, the NoW simulations are in reasonable agreement with the observed mass-size relation of ETGs (Oser et al. 212) and exhibit a diversity of stellar kinematical properties similar to that of real ETGs (Paper XXV). On the other hand, not all simulated galaxies, in particular the W sample, would necessarily meet the morphological selection of the ATLAS 3D sample, which requires the absence of large-scale spiral arms in a galaxy optical image (Paper I). This is essentially a selection against strong star-formation on a disc while about half of all simulated galaxies have relatively large specific star formation rate (> 1 11 yr 1 ) at z =. Similar values are measured in only a handful of the observed ETGs (Davis et al., submitted). We note however that the uncertainty on the simulated star formation rate is large because of uncertainties in, e.g., star formation and feedback physics and their actual implementation. Additionally, neither of the simulated samples include AGNs, which could be an efficient way to suppress star formation. Therefore, we decide not to exclude star-forming galaxies from the comparison. Instead, we highlight this as a caveat of this work and stress that more sophisticated simulations would be needed to fully replicate the ATLAS 3D selection. The simulations do not include a full treatment of different gas phases and the correspondence between simulated gas particles and H I is not obvious. As a first step we define cold gas following the criterion of Hirschmann et al. (212), which is based on a selection of gas particles on the density-temperature plane. This is preferred to a simple temperature selection since dense star-forming gas (which we want to include in the cold-gas budget) can be relatively hot in the NoW simulations. Additionally, we assume that this cold gas is made of two phases, atomic and molecular (ignoring that some of it may in fact be ionized). We separate the two phases on the gas face-on view by applying equation 39 of Krumholz, McKee & Tumlinson (29). In their equation s is the H I+H 2 column density in M pc 2 (see Schruba et al. 211 for a comparison between this theoretical relation and observations). This results in a relatively sharp transition from an atomic- to a molecular dominated regime at a total gas column density of 1 M pc 2 (Wong & Blitz 22; Bigiel et al. 28). The resulting face-on H I images are used to calculate H I mass and column density distribution for all simulated galaxies (at the typical gas resolution of.5 kpc). In order to study the morphology and kinematics of the H I we construct seven different projections of each simulated galaxy relative to an hypothetical observer: face on (relative to the stellar kinematics); rotated by 3, 6 and 9, respectively, about the stellar morphological major axis in the face-on view; and rotated by 3, 6 and 9, respectively, about the stellar morphological minor axis in the face-on view. When constructing these projections we do not separate H I from H 2 as Krumholz s prescription is strictly valid only when the gas is seen face-on. This does not affect our conclusions since we are mainly concerned with the kinematics of the gas, which is the same for all cold-gas phases. All gas mass and surface density calculations as well the analysis of the H I and stellar kinematics presented below are performed after removal of star and gas particles associated with satellite darkmatter sub-haloes according to the SUBFIND algorithm of Springel, Yoshida & White (21). Furthermore, we exclude from the calcu- log1 M(H i)/m NoW simulations W simulations spiral M(H i)/m histogram H i undetected ETGs H i detected ETGs log 1 M (M ) Figure 6. M(H I)/M plotted against M for simulated galaxies (black circles and crosses) and real ETGs from Paper XIII. For the latter, the blue area shows the kernel density estimate derived using all 53 galaxies with detected H I while the red dots indicate H I non-detections and are therefore upper limits on M(H I)/M. The dotted histogram along the vertical axis represents the M(H I)/M distribution of spiral galaxies derived as explained in the text. lation of gas mass and surface density the discrete cold gas clouds which are a known artefact of smoothed-particle-hydrodynamics simulations and which are not part of the gas discs we wish to study (Sijacki et al. 212). These clouds, visible in some of the images discussed in Sec. 5.1, have a typical mass above 1 7 M and may therefore have a significant impact on the value of M(H I) if included in the calculation. Before showing the result of the comparison between simulated and real ETGs we note that the ATLAS 3D sample is volume limited and complete down to M K = 21.5 mag (Paper I). On the contrary, the sample of simulated galaxies is not volume limited (it is a set of individual simulations rather than a single, large simulated volume) and does not follow the observed galaxy luminosity function. Therefore, we cannot compare the statistical distribution of real and simulated galaxy properties. Instead, we investigate whether the diversity of ETG H I properties is reproduced by the simulations. 5.1 H I mass and morphology The first step of our comparison is presented in Fig. 6, where we plot M(H I)/M against M for galaxies in the two sets of simulations. Stellar mass values are taken from Paper XXV for the NoW simulations and are measured in an identical way for the W simulations. For comparison we show the distribution of real ETGs distinguishing between H I detections (blue density field) and non-detections (red dots). We also show the M(H I)/M histogram of spiral galaxies. This is obtained from the M(H I)/L K histogram shown in Paper XIII assuming that for M > 1 1 M the typical spiral has M /L K =.8 M /L (Bell et al. 23). The figure shows that the NoW simulations produce galaxies whose H I content is to first order comparable to that of real ETGs and lower than that of spirals. The main discrepancy between NoW simulated and real ETGs is that extremely gas-poor galaxies with M(H I)/M below a few times 1 4 are underreprec 213 RAS, MNRAS, 1 18

13 H I discs in fast and slow rotators 13 M163 (NoW) M163 (W) M19 (NoW) M19 (W) M329 (NoW) M329 (W) velocity field stars velocity field gas gas overlay stars cm 2 ±28 ± cm 2 ±25 ± cm 2 ±21 ± cm 2 ±33 ± cm 2 ±31 ± cm 2 ±26 ±1 -V -V/2 +V/2 +V Figure 7. Three examples of simulated galaxies. For each galaxy we show the result of the NoW and the W simulation next to each other to highlight their differences. Each column represents one simulated object four rows per galaxy viewed at an inclination of 6. The first row corresponds to an extremely deep optical image ( 29 mag arcsec 2 in V band). The remaining rows are equivalent to those shown in Fig. 4 for real ETGs. Namely, the second row corresponds to a SDSS-like optical image with H I contours overlaid; the third row shows the intensity-weighted H I velocity field; the last row shows the intensity-weighted stellar velocity field. The linear scale of the various panels is the same as in Fig. 4, with the grey scale-bar representing 5 kpc. Galaxies are sorted according to decreasing dark-matter halo mass (left to right). The outer stellar isophote corresponds to a stellar mass surface density of 5 M pc 2. The colour scheme of the velocity fields is represented by the horizontal colour bar at the bottom of the figure. sented in the simulations. This results in an overall smaller range of M(H I)/M at any given M compared to the observations. Nevertheless, both the fact that such range is large in the simulations and that the typical M(H I)/M decreases with M (above a stellar mass of M ) are important points of agreement with the observed ETG sample. The situation is different for the W simulations. These galaxies are typically gas-richer than NoW ones (Hirschmann et al. 213). They populate the high-m(h I)/M end of the observed ETG distribution. At M > 1 11 M we find a number of objects whose H I content is too large by a factor of 1 to 1 compared to the data. We conclude that W systems are, as a family, too H I rich and that NoW galaxies are a better (albeit not perfect) match to real ETGs in terms of H I mass. Having established this we can now analyse the H I morphology of the simulated galaxies. We show three examples in Fig. 7 and another three in Fig. 8. For each galaxy we show the NoW and the W simulation next to each other. These images and velocity fields are immediately comparable to those in Fig. 4 and are drawn adopting the same linear scale and gas contour levels. The only difference between observed and simulated galaxy figures is that for the latter we also show the equivalent of a deep optical image (top panel, see figure caption). Figures 7 and 8 are meant to show a few cases representative of the diversity of galaxies produced by the simulations as well as the difference between the NoW and W simulation for a same halo. Indeed, it is immediately obvious that changing the star formation feedback implementation has a strong impact on the properties of a simulated galaxy. The same object can be a FR in the NoW simulation and a SR in the W one (M163) or vice versa (M19, M329). Furthermore, it can be surrounded by a large gas distribution in the W simulation but have the gas disc entirely confined inside the stellar body in the NoW simulation (e.g., M19, M721). In fact, the larger size of the H I distribution in the W simulations is a general result valid for all galaxies and a consequence of the mechanical nature of the W feedback, which pushes gas to large radius while keeping it cold. A visual comparison between Fig. 4 and Figs. 7 and 8 shows that these large W gas discs match reasonably well the size of the H I discs and rings in gas-rich ETGs. They too, like the observed systems, exhibit slight warps in the outer regions, possibly owing to their large size. The NoW gas distributions are on the contrary too small. The fact that the NoW simulations are a better match to the observed ETG M(H I) distribution (Fig. 6) but that W objects are a better match to the observed size of the H I discs may seem contradictory. In fact, this is easily understood noting that in all simuc 213 RAS, MNRAS, 1 18

14 14 Paolo Serra et al. M549 (NoW) M549 (W) M721 (NoW) M721 (W) M1196 (NoW) M1196 (W) velocity field stars velocity field gas gas overlay stars cm 2 ±23 ± cm 2 ±17 ± cm 2 ±24 ± cm 2 ±23 ± cm 2 ±24 ± cm 2 ±19 ±13 -V -V/2 +V/2 +V Figure 8. Three more examples of simulated galaxies. See caption of Fig. 7. lations the H I column density is on average higher than observed. For this reason H I discs with the right mass are too small (NoW) and H I discs with the right size are too massive (W). This result, which is clear from a comparison of the H I contours in Fig. 4 with those in Figs. 7 and 8, is shown more quantitatively in Fig. 9. In this figure the black lines represent the H I column density distribution in simulated galaxies while the red solid line shows the distribution in real galaxies (represented by the best-fitting Schechter function derived in Paper XIII). This mismatch is consistent with the fact that simulated galaxies appear to be forming stars at a rate higher than real ETGs at z = (Hirschmann et al. 213). Although not the focus of this study, it is worth noting in this context that all simulated galaxies host cold gas above a column density of 1 21 cm 2. This gas, which we do not include in our calculation of H I mass and column density, is more centrally concentrated than the H I and can be related to the molecular gas found in the central region of ETGs (Paper IV). Its mass is typically between and M with a handful of exceptions at lower mass and no significant difference between NoW and W simulations. Both sets of simulations miss therefore the large population of ETGs with M(H 2 ) < M reported in Paper IV while on the gas-rich side of the distribution the molecular gas masses are in reasonable agreement with the observations. We find that the molecular-to-total cold-gas fraction is larger in the NoW than in the W simulations. This is a consequence of their higher average coldgas column density coupled with our method for separating atomic from molecular gas. 5.2 H I discs and misalignment in simulated fast and slow rotators Figures 7 and 8 show that not only the H I mass and morphology but also the gas kinematical misalignment is affected by the feedback implementation. For example, M549 hosts a small counterrotating gas disc in the NoW simulation and a large, corotating one in the W simulation; M721 hosts a small corotating disc in the NoW simulation and a very large, counterrotating one in the W simulation. This kind of difference between the H I-misalignment of NoW and W galaxies is very common within the sample of simulations analysed here. It is caused by the cumulative effect of galactic stellar winds on a galaxy and the satellites which it accretes as the simulation progresses. Taken together, the NoW and the W simulations seem to reproduce qualitatively the large variety of gas morphologies and kinematics of real H I-rich ETGs. In this sense it is useful to associate some of the examples in Figs. 7 and 8 to ETGs in the ATLAS 3D sample. Bearing in mind that these matches are not exact and may not hold at a more quantitative level (as we illustrate below), M163-NoW is for example a FR with a small polar gas disc which may be related to NGC M163-W, a bright SR with a polar gas disc slightly larger than the stellar body, could be compared to NGC M549-NoW hosts a small, counterrotating gas disc similar to NGC 332. M1196-NoW and many W galaxies (e.g., M19, M1196) represent cases of large, corotating H I discs around FRs (similar for example to NGC 5582 and NGC 3838). M329-W exhibits an H I warp and tails comparable to those of NGC Counterrotating gas discs around FRs (e.g., NGC 6798) c 213 RAS, MNRAS, 1 18