The Neutral Hydrogen Properties of Galaxies in Gas-rich Groups
Robert D\v{z}ud\v{z}ar (1), Virginia Kilborn (1), Gerhardt Meurer (2),, Sarah M. Sweet (1, 3), Michael Drinkwater (4), Kenji Bekki (2), Fiona, Audcent-Ross (2), Baerbel Koribalski (5), Ji Hoon Kim (6), Mary Putman (7),, Emma Ryan-Weber (1, 3), Martin Zwaan (8)

TL;DR
This study analyzes neutral hydrogen properties in 27 galaxies within nine gas-rich, low-mass groups, revealing that central galaxies have typical HI content while satellites often show reduced gas likely due to tidal stripping, with internal properties dominating over environmental effects.
Contribution
It provides new insights into HI content and star formation in galaxies within early-stage, gas-rich groups, highlighting the role of internal dynamics over environmental influences.
Findings
Central galaxies have normal HI gas fractions.
Satellite galaxies often show reduced HI due to tidal stripping.
Large HI disks correlate with high angular momentum.
Abstract
We present an analysis of the integrated neutral hydrogen (HI) properties for 27 galaxies within nine low mass, gas-rich, late-type dominated groups which we denote "Choirs". We find that majority of the central Choir galaxies have average HI content: they have a normal gas-mass fraction with respect to isolated galaxies of the same stellar mass. In contrast, we find more satellite galaxies with a lower gas-mass fraction than isolated galaxies of the same stellar mass. A likely reason for the lower gas content in these galaxies is tidal stripping. Both the specific star formation rate and the star formation efficiency of the central group galaxies are similar to galaxies in isolation. The Choir satellite galaxies have similar specific star formation rate as galaxies in isolation, therefore satellites that exhibit a higher star formation efficiency simply owe it to their lower gas-mass…
Click any figure to enlarge with its caption.
Figure 10
Figure 11
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6
Figure 7
Figure 8
Figure 9
Figure 12
Figure 13
Figure 14
Figure 15| Group ID | 1.5/750/EW [h] | Phase Calibrator | fC [MHz] | [′′] | rms [mJy beam-1] |
| (1) | (2) | (3) | (4) | (5) | (6) |
| J0205-55 | 10.49/8.45/9.22 | PKS 0252-712 | 1391 | 45.12 58.94 | 1 |
| J0258-74 | 3.52/3.45/3.95 | PKS 0252-712 | 1398 | 37.49 45.57 | 1.5 |
| J0400-52 | 6.49/3.85/6.12 | PKS 0420-625 | 1371 | 44.96 56.76 | 1.2 |
| J1051-17 | 8.84/2.46/4.34 | PKS 1127-145 | 1393 | 29.77 108.53 | 1.3 |
| J1250-20 | 7.09/2.07/3.74 | PKS 1245-197 | 1384 | 30.79 102.8 | 2 |
| J2027-51 | 5.15/1.50/2.76 | PKS 1933-587 | 1393 | 37.52 50.53 | 1.65 |
| VLA:DnC | |||||
| J1059-09 | 3.18 | J1130-1449 | 1383 | 26.27 44.04 | 1.6 |
| J1026-19 | 2.98 | J1130-1449 | 1379 | 26.10 46.34 | 1.4 |
| J1408-21 | 2.98 | J1337-1257 | 1380 | 28.62 46.44 | 1.4 |
| HyperFit relation | 1 |
|---|---|
| log(M/M⋆) = (4.490.28) (0.490.03) log M⋆ | 0.14 |
| log sSFR = (8.260.39) (0.1950.039) log M⋆ | 0.19 |
| log SFE = (12.770.50) (0.290.05) log M⋆ | 0.25 |
| log M = (4.590.28) (0.250.01) MR | 0.143 |
| Property | Choir central | Choir satellite | Isolated galaxies | ||
|---|---|---|---|---|---|
| galaxies | galaxies | log M⋆ [M⊙] < 8.8 | 8.8 < log M⋆ [M⊙] < 10.5 | M⋆ [M⊙] > 10.5 | |
| [yr-1] | -9.50.3 | -9.680.49 | -10.270.39 | -9.930.18 | -9.520.16 |
| [dex] | -0.990.33 | -0.520.27 | 0.370.19 | -0.280.21 | -0.920.12 |
| [yr-1] | -10.470.17 | -10.210.34 | -9.890.23 | -10.140.12 | -10.470.12 |
| [Gyr] | 4.34.5 | 12.725.5 | 26.720.7 | 9.34.4 | 3.51.2 |
| Galaxy ID | R.A. | Dec. | D | MR | log M⋆ | log M | DEF | SFR | t | D | q | Morphology |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| HIPASS+ | [hh mm ss.s] | [dd mm ss.s] | [Mpc] | [mag] | [M⊙] | [M⊙] | [dex] | [M⊙ yr-1] | [Gyr] | [kpc] | ||
| (1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) | (11) | (12) | (13) |
| J0205-55:S1 | 02 05 05.48 | -55 06 42.54 | 93 | -23.010.32 | 11.440.16 | 10.430.09 | -0.07 | 9.91.9 | 2.7 | 98.54.9 | 0.09 | SAB(r)bc |
| J0205-55:S2 | 02 04 50.78 | -55 13 01.55 | 93 | -20.770.26 | 10.350.13 | 9.840.09 | -0.04 | 2.990.61 | 2.3 | 47.63.6 | 0.2 | SB(s)cd pec |
| J0258-74:S1 | 02 58 06.48 | -74 27 22.79 | 70 | -22.130.24 | 11.010.12 | 10.270.03 | -0.13 | 3.20.6 | 5.8 | 82.45.3 | 0.26 | SAB(rs)bc HII |
| J0258-74:S2 | 02 58 52.43 | -74 25 53.25 | 70 | -19.480.27 | 9.730.12 | 9.80.1 | -0.33 | 0.350.07 | 18.4 | 52.82.9 | S | |
| J0258-74:S3 | 02 58 42.76 | -74 26 03.55 | 70 | -18.510.33 | 9.260.16 | 8.50.2 | 0.72 | 0.130.03 | 2.5 | S | ||
| J0400-52:S1 | 04 00 40.82 | -52 44 02.71 | 151 | -22.310.25 | 11.100.12 | 10.660.09 | -0.48 | 0.50.4 | 91.9 | 77.33.8 | 0.26 | SA(rs)cd pec |
| J1250-20:S1 | 12 50 52.84 | -20 22 15.65 | 114 | -22.540.24 | 11.210.12 | 10.480.09 | -0.24 | 6.61.3 | 4.6 | 106.44.7 | 0.17 | SA(s)bc pec HII |
| J1250-20:S2 | 12 50 40.91 | -20 20 06.22 | 114 | -22.060.24 | 10.980.12 | 9.60.1 | 0.56 | 3.70.8 | 1.0 | 29.29.1 | S | |
| J2027-51:S1 | 20 28 06.39 | -51 41 29.83 | 87 | -21.730.25 | 10.820.13 | 10.20.1 | -0.14 | 5.41.1 | 2.7 | 56.44.2 | 0.24 | (R’)SA(s)bc |
| J2027-51:S2 | 20 27 31.97 | -51 39 20.81 | 87 | -21.910.24 | 10.910.12 | 10.00.1 | 0.07 | 4.30.9 | 2.4 | 42.43.8 | 0.21 | SAB(s)bc pec |
| J2027-51:Sm* | 20 28 08.4 | -51 28 31 | 87 | 9.490.13 | 25.93.8 | [Irr] | ||||||
| J1051-17:S1 | 10 51 37.45 | -17 07 29.23 | 83 | -21.780.25 | 10.840.12 | 10.430.02 | -0.38 | 1.50.3 | 17.7 | 99.66.3 | 0.67 | (R? + PR?) |
| J1051-17:S2 | 10 51 15.11 | -17 00 29.44 | 83 | -21.850.24 | 10.880.12 | 9.450.05 | 0.62 | 2.20.4 | 1.3 | SABb? | ||
| J1051-17:S3 | 10 51 35.94 | -16 59 16.80 | 83 | -18.160.26 | 9.080.13 | 8.730.07 | 0.41 | 0.1040.021 | 5.2 | dS | ||
| J1051-17:NED1* | 10 52 15.9 | -17 07 48.6 | 83 | 9.480.03 | SA(s)a | |||||||
| J1051-17:NED6* | 10 50 49.6 | -17 14 07 | 83 | 9.010.03 | [Irr?] | |||||||
| J1059-09:S1 | 10 59 16.25 | -09 47 38.16 | 122 | -23.240.24 | 11.550.12 | 10.160.10 | 0.26 | 8.91.8 | 1.6 | 96.36.3 | SAB(rs)b pec | |
| J1059-09:S2 | 10 59 06.77 | -09 45 04.38 | 122 | -20.200.24 | 10.070.13 | 9.60.1 | 0.02 | 1.30.3 | 3.3 | 50.43.8 | Irr | |
| J1059-09:S3 | 10 59 15.61 | -09 48 59.41 | 122 | -21.450.25 | 10.680.12 | 10.20.1 | -0.19 | 3.70.7 | 3.9 | Sab pec sp | ||
| J1059-09:S5 | 10 59 30.98 | -09 44 25.26 | 122 | -19.940.28 | 9.950.14 | 9.10.1 | 0.45 | 0.80.2 | 1.7 | S | ||
| J1059-09:S6 | 10 59 08.46 | -09 43 14.49 | 122 | -19.330.28 | 9.650.14 | 8.90.1 | 0.46 | 0.370.07 | 2.5 | Irr | ||
| J1026-19:S1 | 10 26 40.81 | -19 03 04.01 | 135 | -22.740.25 | 11.310.12 | 10.30.1 | -0.02 | 4.50.9 | 4.6 | 136.46.5 | SAB(s)bc: pec Sbrst | |
| J1026-19:S2 | 10 26 50.07 | -19 04 31.77 | 135 | -20.600.26 | 10.270.13 | 9.40.1 | 0.37 | 1.30.3 | 1.9 | Irr | ||
| J1026-19:S5a | 10 26 42.07 | -19 07 35.07 | 135 | -16.940.41 | 8.50.2 | 9.30.1 | -0.44 | 0.030.01 | dIrr | |||
| J1408-21:S1 | 14 08 42.04 | -21 35 49.82 | 128 | -23.150.24 | 11.510.12 | 10.250.22 | 0.14 | 7.11.4 | 2.5 | 134.73.7 | SB(rl)c | |
| J1408-21:S3 | 14 08 41.04 | -21 37 40.97 | 128 | -20.720.25 | 10.330.12 | 9.70.1 | 0.08 | 0.80.2 | 6.2 | S0 | ||
| J1408-21:Sm* | 14 08 46.6 | -21 27 09 | 128 | 9.50.2 | [Edge-on] |
| Galaxy ID | R.A. | Dec. | D | log M | MR | log M⋆ | DEF | SFR | t | Morphology |
|---|---|---|---|---|---|---|---|---|---|---|
| HIPASS+ | [hh mm ss.s] | [dd mm ss] | [Mpc] | [M⊙] | [mag] | [M⊙] | [dex] | [M⊙ yr-1] | [Gyr] | |
| (1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) | (11) |
| J0328-08a | 03 28 05.6 | -08 23 28 | 17.4 | 9.780.09 | -20.150.28 | 10.050.14 | -0.14 | 0.650.13 | 9.3 | SA(s)cd |
| J0412+02 | 04 12 47.3 | +02 22 44 | 70.6 | 10.20.1 | -22.270.24 | 11.080.12 | -0.03 | 3.2110.643 | 5.0 | SB(s)b |
| J0430-01b | 04 30 55.8 | -01 59 13 | 36.1 | 9.50.1 | -19.450.26 | 9.710.12 | -0.07 | 0.5010.101 | 6.8 | SAB(s)b pec |
| J0512-57 | 05 12 09.8 | -57 24 34 | 19.3 | 9.250.09 | -19.140.24 | 9.560.12 | 0.14 | 0.3570.07 | 5.0 | SB(s)m |
| J0515-41 | 05 15 01.9 | -41 21 33 | 14.5 | 9.030.09 | -17.280.24 | 8.660.12 | -0.11 | 0.0650.013 | 16.7 | IB(s)m |
| J1019-17 | 10 19 39.0 | -17 43 48 | 11 | 8.340.12 | -15.250.26 | 7.670.13 | 0.07 | 0.0030.01 | 64.5 | Sm |
| J1051-19 | 10 51 24.6 | -19 53 07 | 31 | 9.760.09 | -20.140.24 | 10.040.12 | -0.12 | 0.670.13 | 8.6 | SB(s)cd |
| J1120-21 | 11 20 12.3 | -21 28 32 | 18.3 | 9.220.09 | -17.950.27 | 8.980.13 | -0.13 | 0.0890.018 | 18.7 | (R’)SAc pec |
| J1123-08 | 11 23 28.1 | -08 38 52 | 56.5 | 10.220.09 | -22.590.24 | 11.240.12 | -0.03 | 7.7321.547 | 2.2 | SB(r)bc |
| J1558-10 | 15 58 22.3 | -10 31 48 | 14.8 | 8.780.11 | -16.410.25 | 8.230.12 | -0.08 | 0.0190.039 | 30.9 | BCD |
| J1621-02ab | 16 21 45.8 | -02 16 34 | 26.3 | 9.750.09 | -21.880.24 | 10.890.12 | 0.33 | 2.4760.495 | 2.3 | SA(s)cd |
| J1954-58 | 19 54 23.4 | -58 42 47 | 31.1 | 10.070.09 | -21.850.24 | 10.880.12 | 0.18 | 3.630.73 | 3.2 | SAB(rs)c |
| J2039-63c | 20 39 07.3 | -63 48 50 | 23.6 | 8.690.14 | -17.080.27 | 8.560.13 | -0.14 | 0.1220.024 | 4.0 | dS0p |
| J2127-60 | 21 27 23.2 | -60 01 22 | 24.9 | 9.940.09 | -20.790.25 | 10.360.12 | -0.03 | 1.1950.239 | 7.3 | SAB(rs)c |
| J2152-55 | 21 52 12.3 | -55 35 12 | 43.5 | 10.250.09 | -22.460.25 | 11.170.12 | 0.1 | 3.6210.724 | 4.9 | SAB(rs)bc |
| J2157-60 | 21 57 52.4 | -60 18 27 | 24.2 | 8.760.12 | -17.040.29 | 8.540.14 | -0.06 | 0.0330.007 | 17.7 | IB(s)m |
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
The Neutral Hydrogen Properties of Galaxies in Gas-rich Groups
Robert Džudžar1, Virginia Kilborn1, Gerhardt Meurer2, Sarah M. Sweet1,3, Michael Drinkwater4, Kenji Bekki2, Fiona Audcent-Ross2, Baerbel Koribalski5, Ji Hoon Kim6, Mary Putman7, Emma Ryan-Weber1,3, Martin Zwaan8, Joss Bland-Hawthorn9, Michael Dopita10, Marianne T. Doyle-Pegg11, Ed Elson12, Kenneth Freeman10, Dan Hanish13, Tim Heckman14, Robert Kennicutt14, Pat Knezek15, Martin Meyer2,3, Chris Smith16, Lister Staveley-Smith2, Rachel Webster17, Jessica Werk18
1 Centre for Astrophysics and Supercomputing, Swinburne University of Technology, PO Box 218, Hawthorn, VIC 3122, Australia
2 International Centre for Radio Astronomy Research, ICRAR M468, 35 Stirling Highway, Crawley, WA 6009, Australia
3 ARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D)
4 School of Mathematics and Physics, University of Queensland, Brisbane, QLD 4072, Australia
5 Australia Telescope National Facility, CSIRO, P.O. Box 76, Epping, NSW 1710, Australia
6 Subaru Telescope, National Astronomical Observatory of Japan, 650 North A’ohoku Place, Hilo, HI 96720, USA
7 Department of Astronomy, Columbia University, 550 West 120th Street, New York, 10027, USA
8 European Southern Observatory, Karl-Schwarzschild-Strasse 2, D-85748 Garching bei München, Germany
9 Sydney Institute for Astronomy, University of Sydney, Sydney, NSW 2006, Australia
10 Research School of Astronomy and Astrophysics (RSAA), Australian National University, Cotter Road, Weston Creek, ACT 2611, Australia
11 Department of Physics, University of Queensland, Brisbane, QLD 4072, Australia
12 Department of Physics & Astronomy, University of the Western Cape, Cape Town 7535, South Africa
13 Department of Physics and Astronomy, John Hopkins University, 3400 North Charles Street, Baltimore, MD 21218, USA
14 Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK
15 National Science Foundation, 4201 Wilson Boulevard, Arlington, Virginia 22230, USA
16 Cerro Tololo Inter-American Observatory (CTIO), Casilla 603, La Serena, Chile
17 School of Physics, University of Melbourne, VIC 3010, Australia
18 Astronomy Department, University of Washington, 3910 15th Ave. NE, Seattle, WA 98195-0002, USA E-mail: [email protected]
(Accepted XXX. Received YYY; in original form ZZZ)
Abstract
We present an analysis of the integrated neutral hydrogen (H i) properties for 27 galaxies within nine low mass, gas-rich, late-type dominated groups which we denote “Choirs". We find that majority of the central Choir galaxies have average H i content: they have a normal gas-mass fraction with respect to isolated galaxies of the same stellar mass. In contrast, we find more satellite galaxies with a lower gas-mass fraction than isolated galaxies of the same stellar mass. A likely reason for the lower gas content in these galaxies is tidal stripping. Both the specific star formation rate and the star formation efficiency of the central group galaxies are similar to galaxies in isolation. The Choir satellite galaxies have similar specific star formation rate as galaxies in isolation, therefore satellites that exhibit a higher star formation efficiency simply owe it to their lower gas-mass fractions. We find that the most H i massive galaxies have the largest H i discs and fall neatly onto the H i size-mass relation, while outliers are galaxies that are experiencing interactions. We find that high specific angular momentum could be a reason for galaxies to retain the large fraction of H i gas in their discs. This shows that for the Choir groups with no evidence of interactions, as well as those with traces of minor mergers, the internal galaxy properties dominate over the effects of residing in a group. The probed galaxy properties strengthen evidence that the Choir groups represent the early stages of group assembly.
keywords:
galaxies: general – galaxies: evolution – galaxies: groups – galaxies: ISM – galaxies: interactions – galaxies: groups: individual: HIPASSJ1250-20
††pubyear: 2018††pagerange: The Neutral Hydrogen Properties of Galaxies in Gas-rich Groups–LABEL:lastpage
1 Introduction
Today, we observe the hierarchical structure of the Universe (Press & Schechter, 1974; White & Rees, 1978) ranging from sparse voids to the densest of galaxy clusters. Small groups are the most common environment where galaxies reside (e.g. Tully 1987; Eke et al. 2004; Knobel et al. 2015; Saulder et al. 2016). Since environment has an impact on galaxy evolution, it is important to understand how galaxies evolve within such structures. Environmental impact can be traced via galaxy properties such as the observed morphology-density relation (the decline of the fraction of late-type galaxies with increasing environmental density) in clusters (Dressler, 1980; Goto et al., 2003) as well as in groups (Postman & Geller, 1984; Serra et al., 2012). Furthermore, studies have shown that galaxy’s color and morphology are connected with the galaxy environment (Whitmore et al., 1993; Kauffmann et al., 2004; Blanton et al., 2005; Brough et al., 2006; Melnyk et al., 2014), and that the star formation rate suppression becomes evident at group densities (Lewis et al., 2002; Gómez et al., 2003; Davies et al., 2015; Barsanti et al., 2018). Nevertheless, the relative contribution of nature vs. nurture (internal properties of a galaxy versus an environment in which galaxy resides) in small gas-rich groups remains an open question with arguments on both sides (e.g. Serra et al. 2012; Knobel et al. 2015; Odekon et al. 2016; Carollo et al. 2016; Brown et al. 2016; Janowiecki et al. 2017; Spindler et al. 2018).
Environmental effects can be probed by mapping the distribution of the neutral hydrogen (H i) content. The H i gas in spiral galaxies usually extends from 1.5–2 times the size of their observed stellar disc measured at the R-band 25 mag arcsec*-2* surface brightness level (e.g. Bosma 1981; Hoffman et al. 1996; Pisano et al. 2000; Bigiel & Blitz 2012). Due to the extent of the H i disc, it is more susceptible than the stellar disc to disturbances via external gravitational (tidal interactions, galaxy mergers) and hydrodynamical (ram-pressure stripping) interactions. These properties make the H i gas an excellent tracer of the physical processes that are affecting galaxies (e.g. Gunn & Gott 1972; Yun et al. 1994; Barnes et al. 2001; Verdes-Montenegro et al. 2001; Hess et al. 2017).
Tidal interactions and galaxy-galaxy mergers are often observed in galaxy groups (e.g. Yun et al. 1994; Rots et al. 1990). Such features are prominent in Hickson Compact Groups (HCG, compact configurations with four or more galaxy members within a magnitude difference of 3.0, obtained using red prints from the Palomar Observatory Sky Survey, Hickson 1982) and theoretically they are often connected to a transformation of late-type galaxies into early-type galaxies (Coziol & Plauchu-Frayn, 2007). Galaxy interactions are also observed in small (fewer members, Barnes et al. 2001) and loose (sparser, Oosterloo et al. 2018) groups. There are also indications that ram-pressure stripping can remove gas from galaxies in galaxy groups (Rasmussen et al., 2006; Westmeier et al., 2011; Wolter et al., 2015; Brown et al., 2017; Stevens & Brown, 2017) and not just in clusters (Gunn & Gott, 1972; Chung et al., 2009). Enormous effort is aimed towards understanding how such interactions relate to the observed H i content in galaxies and conversely how H i distribution and kinematics uncover the physics of interaction (e.g. Bekki 2008; Pisano et al. 2011; Mihos et al. 2012; Elagali et al. 2018; Bosma 2017, and references therein).
Spiral galaxies have a lower H i fraction in high density environments, for instance, near the centers of galaxy clusters (Giovanelli & Haynes 1985; Chung et al. 2009) or in HCGs (Verdes-Montenegro et al., 2001), than spiral galaxies that reside in the field. Verdes-Montenegro et al. (2001) proposed an evolutionary path for gas in compact galaxy groups. They proposed that the H i gas evolution in the compact groups occurs through a gradual removal of the gas from galaxies, increasing its content in the intragroup medium. This evolutionary picture should be expanded since gas poor galaxies are also found in intermediate density environments (Dénes et al., 2016; Hess & Wilcots, 2013) as well as in loose groups (Kilborn et al., 2005, 2009). With the observations of the gas content in group galaxies we are gathering evidence that gas pre-processing starts to be evident. However, it is still unclear at which density scale the environment starts affecting the gas content of galaxies.
Referring back to hierarchical structure formation theory, a galaxy embedded in the centre of a dark matter halo (central galaxy) grows as it accretes satellite galaxies, over time becoming the most massive and the brightest galaxy (White & Rees, 1978; Brough et al., 2006). In this scenario, the satellite galaxies are moving with respect to the centre of group potential and experiencing environmental effects, thus in theory they are going through a different evolutionary path from that of the central galaxy (White & Rees, 1978; Yang et al., 2007; Brown et al., 2017; White & Rees, 1978; Brough et al., 2006; Stevens & Brown, 2017). Observationally, and in simulations, it is often difficult to distinguish between the central and satellite galaxy, especially in a small groups (Berlind et al., 2006; Campbell et al., 2015).
Comparing central galaxies in groups with isolated galaxies, Janowiecki et al. (2017) found that in low mass groups, central galaxies have a higher gas content than similar galaxies in isolation which they attributed either to inflows of gas from the cosmic web, or, mergers with H i-rich galaxies; whilst in higher mass groups central galaxies have similar gas content as galaxies in isolation. This picture is consistent with the H i-rich galaxies from the Bluedisk sample (Wang et al., 2015) for which they found that H i-rich galaxies are embedded within H i-rich environments. Whether accretion from gas-rich mergers or accretion from the cosmic web is the main reservoir from where galaxies accrete the majority of their gas is still unclear. The contribution from gas-rich mergers is thought to be rather small (Sancisi et al., 2008; Di Teodoro & Fraternali, 2014), and the current observations are not sensitive enough (Popping et al., 2009) to probe the cold mode accretion from cosmic web filaments, which in theoretical models is the dominant reservoir from where galaxies obtain their gas (Kereš et al., 2005).
In this paper we investigate the extent to which the galaxy group environment influences evolution of the galaxies in small gas-rich groups, through a comparison of the group galaxies to a sample of isolated galaxies.
This paper is structured as follows. Section 2 describes our Choirs sample and shows Choir group HIPASSJ1250-20 in detail, describes the selection of isolated galaxies and summarises the observations of group galaxies and data reduction. Section 3 describes the main results of this paper: comparison of group galaxies with isolated galaxies in terms of gas-mass fraction, specific star formation rate, star formation efficiency, atomic depletion time and deficiency. This section also presents the size-mass relation and examines global stability parameter of the Choir galaxies. Section 4 discusses the implications of the obtained results, and Section 5 summarises the results and concludes.
Throughout this paper the assumed cosmology is H0 = 70 km s*-1* Mpc*−1*, = 0.3, and = 0.7.
2 Sample and data
In this work we make use of resolved interferometric H i data for 9 Choir groups, obtained with the Australia Telescope Compact Array (ATCA) and the Karl G. Jansky Very Large Array (VLA). A brief summary of the observations is provided in Table 1. We examine the properties of the individual galaxies within the Choir groups and compare them to the sample of isolated galaxies.
2.1 Choir group galaxies
Our groups are selected from the Survey of Ionization in Neutral Gas Galaxies (SINGG), the H narrow-band imaging follow-up program (Meurer et al., 2006) for the H i Parkes All Sky Survey (HIPASS; Barnes et al. 2001). Galaxies for the SINGG survey were selected based on their H i mass and thus the sample is biased towards more distant, massive galaxies (Meurer et al., 2006; Sweet et al., 2013). The HIPASS galaxies were separated into log(M) bins with width of 0.2 dex and then the nearest galaxies in each log(M) bin were selected. Additional cuts were made in the selection: galaxies were required to have a HIPASS flux with S/N larger than 3.8; galaxies with Galactic latitude lower than 30 deg were not selected to avoid Galactic plane; galaxies with systemic velocities smaller than 200 km s*-1* were not selected to avoid the Milky Way and high-velocity clouds. Detailed explanations can be found in Meurer et al. (2006); Audcent-Ross et al. (2018) and Meurer et al. (in prep.).
27 per cent of the HIPASS detections in the SINGG sample correlate with multiple H emitting galaxies111It should be noted that the number of the single and multiple detections relates only to the field of view of 14.7 (due to the large beamsize of the Parkes telescope at 21cm). We are investigating those cases where four or more H emitting galaxies were detected within a SINGG field-of-view (around seven per cent of the SINGG fields). 15 groups were denoted as SINGG Choirs in Sweet et al. (2013) and they are small groups of four to 10 members. Choir groups masses range between 1011 and 1012 M*⊙* (measured as the sum of H i mass and stellar mass). Choir group properties such as total group mass, specific star formation rates, surface brightness, H equivalent width and group H i deficiency were examined with respect to the entire SINGG sample in Sweet et al. (2013). They show that the galaxies in the Choir groups are very similar to the overall SINGG sample. This paper focuses on the integrated properties of nine Choir groups for which neutral hydrogen content was mapped.
2.2 Isolated SINGG galaxies
In order to examine whether the galaxies within the Choir groups are affected by the group environment, we compile a sub-sample of SINGG galaxies that are isolated and thus unlikely to have been affected by neighbouring galaxies within the past few Gyr. We compile such a sample because we found that outside the SINGG field of view (14.7), single detected galaxies are members of different environments ranging from being isolated to being in clusters. Galaxies can reside in the different environments due to the fact that nearby groups may break up into multiple H i detections and may in part be captured in SINGG as single galaxies (Sweet et al., 2013).
Multiple catalogues of isolated galaxies exist; these have been created by imposing different isolation criteria. The Catalogue of Isolated Galaxies (CIG) (Karachentseva, 1973) imposed the criteria that a galaxy is isolated if it is a single galaxy within a projected distance relative to neighbours diameter: , where D diameter of the primary galaxy, D is diameter of the neighbouring galaxies; thus, the projected distance should satisfy criteria: .
Applying similar criteria to the one from CIG sample, albeit using Ks-band (instead of Zwicky magnitudes) for determining galaxy angular diameters, Karachentseva et al. (2010) made the 2MIG catalogue of isolated galaxies. Verdes-Montenegro et al. (2005) created the AMIGA (Analysis of the interstellar Medium of Isolated GAlaxies) sample of isolated galaxies which have not been affected by other galaxies in the past 3 Gyr (assuming a field velocity of 190 km s*-1*, as adapted from Tonry et al. (2000)). The AMIGA sample was later revised by Argudo-Fernández et al. (2013) who imposed additional constraints to those from Karachentseva (1973) and Verdes-Montenegro et al. (2005) with the advantage of having the photometric SDSS-DR9 catalogue, thus creating a sample of 636 isolated galaxies. The SIGRID sample (the Small Isolated Gas-Rich Irregular Dwarfs), compiled by Nicholls et al. (2011), catalogued galaxies that have not had interactions within the last 5 Gyr. Using their own isolation criteria, Reda et al. (2004) created a sample of isolated early type galaxies with the constraints that there are no galaxies within 700 km s*-1* in recession velocity, within 0.67 Mpc in the plane of the sky, and within 2 B-band magnitudes of the isolated galaxy.
To compare the Choir group galaxies to those in isolation, we compiled a catalogue of isolated galaxies from the SINGG sample. Our isolation criteria closely follows those found in the literature, including a projected distance from the neighbouring galaxies, recessional velocity difference with neighbouring galaxies, as well as proximity to a galaxy cluster. For each SINGG galaxy we firstly determine galaxies which are single within 0.67 Mpc in the plane of the sky and within a velocity of 750 km s*-1*, as imposed in Reda et al. (2004). For this analysis we used NED (NASA/IPAC Extragalactic Database) to find all known galaxies with redshifts around each galaxy from the SINGG sample, thus the data varies in sensitivity. Applying restrictions on projected distance and velocity offset we obtained a sample of 25 galaxies. We can approximate the time since the last encounter for these galaxies, assuming that the galaxy perturber has an average diameter D 38 kpc and field velocity of 190 km s*-1*. Thus a perturber of the same size would need 4 Gyr to traverse the distance equal to 20Da. This time is similar to the one found for the AMIGA, SIGRID and 2MIG samples (we do find that some galaxies from our sample are also members of the SIGRID Melnyk et al. (2015) and 2MIG Karachentseva et al. (2010) samples). Secondly, we remove six galaxies which are within the projected radius of the Virgo cluster and within two times its velocity dispersion, adopting properties of the Virgo cluster given by Nicholls et al. (2011). By removing galaxies which are in the vicinity of the cluster we discarded possible interlopers which could be infalling cluster galaxies. We additionally flagged three galaxies, one with a noisy H i spectrum that could led to an error in the H i mass estimation, one due to large uncertainty in its stellar mass and one with three H clumps spreading outside of its stellar disc. With this procedure we created a sample of 16 isolated galaxies. The isolated galaxies and their properties are shown in Table 5. We acknowledge that it is possible, in spite of these isolation criteria, to have contamination in a sample of isolated galaxies due to low surface brightness galaxies. However within the 14.7 FOV of these galaxies in SINGG, emission line galaxies below R-band apparent magnitude of 20 were not detected (Meurer et al., 2006), thus contaminants should have negligible effects.
2.3 ATCA & VLA observations and data reduction
In this paper, we utilise observations obtained with the CFB 64M-32k (34 kHz resolution), CFB 1M-0.5k (0.5 kHz resolution) correlator configuration on the Compact Array Broadband Backend (CABB; Wilson et al. 2011) as well as a set of observations (HIPASSJ0205-55a) with the previous correlator (16 kHz resolution).
We also utilise observations obtained with the Karl G. Jansky Very Large Array (VLA, project 13A-207) in the hybrid DnC configuration, having a total bandwidth of 4 MHz, divided into 1024 channels with a width of 3.906 kHz.
The ATCA and VLA data were reduced using standard procedures with the MIRIAD (Sault et al., 1995) and CASA (McMullin et al., 2007) software packages, respectively. For the ATCA data we first removed radio frequency interference and then calibrated (bandpass, flux and phase) each data set separately. In the next step, we performed a continuum subtraction using a linear fit to line-free channels in each observational data set. We combined all baseline configurations for each data set in the Fourier transformation using task INVERT. This task was utilised excluding antenna six (CA06), using Brigg’s robust parameter of 0.5 and a channel width of 5 km s*-1*. The resulting cube for each data set was cleaned, restored and primary beam corrected. With the task MOMENT we created moment 0, 1 and 2 maps with a 3 clipping.
The VLA data were reduced in the same manner as the ATCA data using CASA. Firstly, we flagged data, calibrated the observations and then subtracted the continuum. Secondly, we cleaned our data cubes using a robust parameter of 0.5 and then we corrected each data cube for the primary beam. Finally, we created moment 0, 1 and 2 maps with a 3 clipping.
In this work we utilise the obtained integrated galaxy properties: H i fluxes, masses and diameters, and we show ATCA H i intensity distribution of Choir HIPASSJ1250-20. The full sample of the moments maps will be published and analysed in a subsequent paper.
3 Results
3.1 The HIPASSJ1250-20 group
We present the HIPASSJ1250-20 (hereafter J1250) Choir group and its complexity as an illustration for our discussion throughout this work. The optical observations of J1250 were carried out with the Dark Energy Camera (DECam) on the CTIO Blanco 4-m telescope in the g, r, i and z bands. The exposure times of the observations were 3000, 1920, 3000 and 6486 seconds respectively, these are the deepest optical images of this group (based solely on the exposure time calculator, these images reach 26 mag arcsec*-2* with S/N3). The optical images of similar depth for eight Choir groups will be included in the next paper. In Figure 1, we show the DECam g, r and i colour composite image of J1250 group and labels of the group members from S1 to S7. These members were presented in Sweet et al. (2013): two large spiral galaxies (S1 and S2; respectively ESO575-G006 and ESO575-G004), three dwarf galaxies (S3, S4 and S5) and two small and compact H emitters that they note as tidal dwarf candidates (S6 and S7).
We observed J1250 group with the Australia Telescope Compact Array and mapped its H i content. The H i observations are described in Section 2.3 and summarised in in Table 1. Figure 2 shows the H i column density contours overlaid on the optical colour composite image of J1250. The lowest shown H i column density, the white contour, is a 3 detection and the H i distribution reveals a possible interaction between S1 and S2 galaxy; we find similar interaction through H i distribution in the HIPASSJ1408-21 group. We must note that because of the large elongation of the synthesized beam (as the source is at 20 degrees in Declination), this could be a result of the beam smearing. We find no evidence of interaction between S1 and S2 galaxy in the optical or H imaging at the sensitivity with which they were imaged (while we do find interaction in optical and H in the HIPASSJ1059-09 and HIPASSJ1026-19 group). We find that the of the H i emission is off centre in S2 galaxy. The total H i mass of J1250 is log M [M*⊙] = 10.510.11 (Sweet et al., 2013), which is in agreement with the mass quoted in the HIPASS catalogue: log M [M⊙] = 10.390.11 (Meyer et al., 2004). The slight difference in obtained mass is due to the difference in adopted group distance, Meyer et al. (2004) uses 100 Mpc, while we are using 114 Mpc as in Sweet et al. (2013). With our ATCA observations we recover the H i flux detected with the single dish telescope and the sum of the masses of (S1 + S2) is log M [M⊙] = 10.530.2. Using the ATCA data we can separate the H i flux in S1 (H i emission between 7630 and 7860 km s-1*) and S2 (H i emission between 7685 and 7880 km s*-1*) galaxies and we find that the HIPASS source is mistakenly assigned to S2 galaxy in the HyperLeda database. We show that the majority of the H i gas appears to be associated with S1 galaxy, moreover the S1 is an H i-rich galaxy whilst S2 is an H i-poor galaxy. The dwarf galaxies S3, S6 and S7 are within the H i envelope associated with S1 galaxy. A small H i component associated with S2 galaxy coincides also with S4 dwarf galaxy (top right part in Figure 2). Galaxy marked S5 is not detected within the ATCA primary beam. We do not have spectra of S3, S4, S5 and S6 galaxies however they are assumed to be the group members as they are detected in SINGG H imaging. This imaging used filters with a bandwidth of 30 Å which corresponds to 3000 km s*-1* (Meurer et al., 2006; Sweet et al., 2013).
Both, S1 and S2 galaxies have comparable stellar masses, while the H i mass of S1 is 7.6 times larger than the H i mass of S2. This is the one of the groups where two galaxies have an equal probability of being the central galaxy (see next Section). The gas-mass fraction of S1 is -0.73 dex and it is within the 2 scatter from the fit to the isolated galaxies while the gas-mass fraction of S2 is -1.42 dex, and it is more than 2 lower from the fit to the isolated galaxies (see Figure 4).
Figure 3 shows the DECam g band image which we adaptively smoothed to enhance the faint stellar streams around the galaxies in J1250 (Illustrating the discussion in Section 3.6). The stellar streams around S1 galaxy (bottom left side in Figure 3) have a blue colour and they are visible in the H emission (Figure B10 in Sweet et al. (2013)), also S6 and S7 galaxies are detected in GALEX All Sky Survey far ultraviolet and near ultraviolet images. The stellar streams around S1 galaxy are possibly caused by a disrupted satellite galaxy. The streams around S2 galaxy (top right side in Figure 3) have a red colour, they do not show emission in the H imaging and they are possibly caused by a tidal interaction with S4 galaxy.
In this paper we present results from the analysis of similar quality H i data (listed in Table 1) for the other eight groups.
3.2 Galaxy properties
Analysing the interferometric observations from ATCA and VLA, we detected H i in 27 galaxies within nine Choir groups. For each galaxy we obtain their H i mass from their H i spectra using the standard formula:
[TABLE]
where D is the distance to the galaxy (Mpc) and F is the integrated flux density (Jy km s*-1*). The distance to the each galaxy is assumed to be equal to the galaxy’s group distance, as used in Sweet et al. (2013) and Meurer et al. (2006), based on the multipole attractor model (Mould et al., 2000). To determine the errors on the integrated flux density we follow the prescriptions by Koribalski et al. (2004), whilst we assume a 10 per cent error on the distance following Meurer et al. (2006). We use standard error propagation on the integrated flux density and the distance to determine the uncertainty of the H i mass. The Choir galaxy H i properties are summarised in Table 4.
We use dust-corrected R-band AB magnitudes (Meurer et al., 2006) to determine stellar masses (M*⋆) using the relation between M⋆*/L and absolute R-band magnitude, as derived for the SINGG sample, equation (4) in Wong et al. (2016). We note that using near-infrared bands for obtaining stellar masses would be more precise since they are less affected by dust extinction and recent star formation, however we do not have Ks fluxes for the entire sample of galaxies. In order to check the reliability of R-band stellar masses, we compute Ks stellar masses where possible, using images from the the 2 Micron All Sky Survey (2MASS; Skrutskie et al. 2006) and deriving stellar masses as prescribed by Wen et al. (2013). On average, there is good agreement between the stellar masses obtained with the R-band and Ks-band, with a mean difference of 0.3 dex.
We derive the star formation rates (SFR), as in Sweet et al. (2013), using the H luminosities which are corrected for Galactic extinction and [NII] contamination (see details of the flux measurement in Meurer et al. 2006). To check for any Active Galactic Nuclei (AGN) contamination we identify the dominant source of ionization in the Choir galaxies using the ‘BPT’ emission line diagnostic diagram of Baldwin et al. (1981). Using the emission line fluxes from Sweet et al. (2014), we examine the position of the Choir galaxies in BPT diagram using the line ratios: [NII]/H versus [OIII]/H (Kauffmann et al., 2003; Kewley et al., 2001) and identify four galaxies whose emission is dominated by AGN. We also use available classifications from Sweet et al. (2013) which is noted from the NED (NASA/IPAC Extragalactic Database) and find three more galaxies dominated by AGN radiation. We excluded these galaxies from the analysis of the specific star formation.
We show the relations between our sample of isolated galaxies and Choir galaxies in terms of gas-mass fraction (M/M*⋆), specific star formation rate (sSFR = SFR/M⋆*) and star formation efficiency (SFE = SFR/M). These relations are analysed by dividing Choir galaxies into group central and satellite galaxies. We denote the brightest group galaxy (based on R-band absolute magnitude) as the central galaxy. For groups with two galaxies of similar R-band magnitudes within uncertainties: HIPASSJ1051-17, HIPASSJ2027-51 and HIPASSJ1250-20, two galaxies have an equal probability to be central galaxy. In these three cases we denote both galaxies as centrals. All other galaxies within the same group are treated as satellites.
Figure 4 presents the H i fraction versus stellar mass relation for the sample of isolated galaxies. We fit the isolated galaxies in this plane with HyperFit222hyperfit.icrar.org package (Robotham & Obreschkow, 2015) and use the dispersion about this fit (vertical scatter) to parameterise how far off of this relationship the Choir galaxies deviate. In terms of gas-mass fraction, only one Choir central galaxy has a higher gas-mass fraction than the isolated galaxies. There are two centrals with smaller gas-mass fraction than the isolated galaxies (that is, which deviate more than 2 from the fitted relation) however, these two central galaxies are special cases. These two galaxies are part of groups (HIPASSJ1051-17 and HIPASSJ1250-20) where we have two central galaxies with similar stellar mass thus they have 50 per cent probability of being central galaxies in the first place. Satellites are more scattered around the relation; we find seven satellite galaxies (63 per cent) more than 2 from the relation. Five of these satellite galaxies have lower H i content, and two have higher H i content with respect to the similar mass isolated galaxies. These are indications that in groups with one central galaxy, the H i gas will be depleted firstly in galaxies with the lower stellar masses (such satellites are in HIPASSJ1059-09, HIPASSJ0258-74 and HIPASSJ1026-19 group). In HIPASSJ1026-19 there is an ongoing tidal interaction between S1 and S2 thus tidal stripping maybe the main gas removal mechanism that is acting on S2. In groups with two central galaxies the gas consumption is more complicated; central galaxies can have similar or vastly different gas content.
Figure 5 shows the specific star formation rate for the entire Choirs sample and isolated galaxies, where colours represent morphological classification of Choir galaxies. We show the HyperFit fit to the isolated galaxies and the associated 1 and 2 scatter. It has been shown that the suppression of star formation begins at group densities Lewis et al. (2002). We do find that the majority of the dwarf satellite galaxies (M 109.5 M*⊙*) have lower specific star formation rate than the isolated dwarf galaxies. Comparing the star formation of central and satellite galaxies, we find that central galaxies are less scattered, only one deviate more than 2 from the fit. We find the galaxies that deviate more than 2 from the relation and have a high specific star formation rate are three dwarf satellites and two peculiar spirals. However, we find larger number of galaxies with a lower specific star formation rate (compared to galaxies of the same stellar mass found in isolation). These galaxies are mostly dwarf satellite galaxies, as well as lenticular and spirals satellite galaxies.
The star formation efficiency for Choirs and isolated galaxies is shown in Figure 6. We fit the isolated galaxies using the HyperFit package. Two central galaxies (HIPASSJ1051-17:S1 and HIPASSJ1250-20:S2) deviate more than 2 from the fit. HIPASSJ1051-17:S1 exhibits lower star formation efficiency due to a higher gas-mass fraction, and HIPASSJ1250-20:S2 has higher specific star formation efficiency due to a lower gas-mass fraction. We illustrate H i fraction offset versus specific star formation rate offset in Figure 7. On the x-axis the galaxy’s offset from the fit to the specific star formation rate ( log sSFR) and on the y-axis the galaxy’s offset from the fit to the H i fraction ( log M/M*⋆*). The best-fit relations are placed in Table 2.
We summarise the mean properties of the Choir central, Choir satellite and isolated galaxies in the Table 3. The mean properties of isolated galaxies are separated into three bins of stellar mass which are matched to correspond to the stellar mass limits of Choir satellites (8.8 < log M*⋆* [M*⊙] < 10.5) and Choir central galaxies (log M⋆* [M*⊙*] > 10.5). The mean value of the gas-mass fraction for satellite galaxies and the ones in isolation differ, they are 0.52 and 0.28 dex respectively; however, due to the large scatter these values are in agreement within 1. Although there is no difference in the mean values, the scatter is considerably larger for the group galaxies, and actually 1/3 are somewhat H i deficient (as described in Section 3.5). Examining the H i fraction, specific star formation rate and star formation efficiency of Choir galaxies yielded a number of galaxy outliers which should be investigated on a group by group basis, which we will do in a future paper.
3.3 Size-mass relation
There is a tight correlation between the H i diameter of a galaxy defined at a surface density of 1 M*⊙* pc*-2* and H i mass (Broeils & Rhee, 1997; Verheijen & Sancisi, 2001). The derived size-mass relation by Broeils & Rhee (1997) for 108 galaxies; log D = 0.51 log(M) 3.32 was recently updated by Wang et al. (2016) using a sample of 562 galaxies who found roughly the same result: log D = 0.506 log(M) 3.293. It has been shown that the D-M relation holds for different galaxy morphologies and for H i-rich or relatively H i-poor galaxies (Broeils & Rhee, 1997; Wang et al., 2016).
We examine Choir galaxies and their position on the D-M relation. We measure the H i diameter from the H i column density maps following prescriptions by Wang et al. (2016). Firstly we de-project H i column density maps by multiplying them by b/a, where is the inclination angle of the galaxy’s disc and b/a is its axis ratio obtained from R-band imaging (Meurer et al., 2006; Sweet et al., 2013). We map the contour density of 1.2 1020 cm*-2* equal to 1 M*⊙* pc*-2*, and measure the diameter of the contour by fitting it with an ellipse using the KMPFIT of the Kapteyn python package (Terlouw & Vogelaar, 2015), see Appendix B for example maps and the obtained fit for HIPASSJ1250-20 group. We use standard error propagation of the uncertainties given by the fit to place errors on the H i diameter. We then correct the obtained H i diameters (D for beam smearing effects (the best suitable for circular beam shape), as described in Wang et al. (2013, 2016): , where the B and B are respectively synthesized beam sizes of the major and minor axes, D is uncorrected and D is corrected H i diameter. The correction for most of the Choir galaxies is small and it reduces the diameter up to 5 per cent. The only exception is the galaxy HIPASSJ1250-20:S2 for which the correction amounts to 32 per cent. This large correction is due to the fact that the major axis of the synthesized beam is larger than the measured H i diameter at a surface density of 1 M*⊙* pc*-2*. There were several galaxies for which the H i diameter is not measured, these are small satellite galaxies and one small central galaxy HIPASSJ1051-17:S2 - which is part of a group with two centrals. These galaxies were excluded from Figure 8 because their H i emission was not detected at a surface density of 1 M*⊙* pc*-2*.
Figure 8 shows the relation between the total H i mass and the diameter of the H i disc of the measured Choir galaxies. The H i size-mass relation prescribed by Wang et al. (2016) is shown as a solid line in Figure 8, and the dashed lines correspond to their 3 scatter which translates into 0.18 dex. We find three out of 15 measured Choir galaxies to be scattered outside 3 of the relation, and all three outliers are disturbed systems (black star symbols in Figure 8). The two systems with the largest diameters: HIPASSJ1408-21:S1+S3 and HIPASSJ1026-19:S1+S2 are merging galaxies and cannot be separated into individual galaxies at a surface density of 1 M*⊙* pc*-2*. Thus the measured H i diameter encompasses both of them. HIPASSJ0400-52:S1 is also disturbed, it shows an H i offset from the stellar center and has a large asymmetry of the H i distribution (Džudžar et al. in prep). Outliers from Wang et al. (2016) also show a broad range of abnormalities, such as highly H i deficient galaxies, merger events and galaxies with an asymmetric H i distribution within a disc. On the other hand, galaxies that have such abnormalities do not necessarily have to reside off the H i size-mass relation. HIPASSJ1250-20:S2 and HIPASSJ1059:S1+S3 are such examples. The H i distribution of HIPASSJ1250-10:S2 is offset from the stellar centre, and extends towards the larger companion (see Figure 2), and HIPASSJ1059:S1+S3 are tidally interacting galaxies.
The mean value of the D for the central Choir galaxies is 72 kpc, while for the satellites it is 51 kpc. The three largest central galaxies are HIPASSJ1250-20:S1, HIPASSJ1051-17:S1 and HIPASSJ0205-55:S1 with D sizes of 106, 100 and 98 kpc respectively; their size is comparable to the largest galaxies from the Bluedisk sample from Wang et al. (2013).
3.4 Atomic gas depletion time
Assuming continuous gas consumption with the current rate of the star formation determined as measured by H, we determined atomic gas depletion time (tdep) for Choir galaxies:
[TABLE]
Figure 9 shows that all except one central and two satellite Choir galaxies have atomic depletion times smaller than the Hubble time (13.5 Gyr). Overall, the mean gas depletion time of all Choir galaxies ( 8.5 Gyr) is smaller than that obtained for the sample of isolated galaxies ( 12.9 Gyr), two thirds of the isolated galaxies have depletion times smaller than the Hubble time. We did not find any starburst time scales (smaller than 1 Gyr) that would indicate fast gas exhaustion and thus relatively rapid transformation of a gas-rich galaxy into a gas-poor one. The two centrals that have the smallest depletion time are HIPASSJ1250-20:S2 and HIPASSJ1051-17:S2, having depletion times of 1.0 Gyr and 1.3 Gyr respectively. The short depletion times of HIPASSJ1250-20:S2 and HIPASSJ1051-17:S2 are attributed to their low gas-mass fraction (as seen in Figure 4: log(M/M*⋆*) 1.4 dex), whilst their specific star formation rates are comparable to those of isolated galaxies of the same stellar mass.
3.5 H i deficiency
We probe the H i content of the Choir galaxies to determine whether it differs from the H i content of the isolated galaxies. It has been shown that spiral galaxies have less H i in high density environments, for instance, near a galaxy cluster center (Chung et al., 2009) or in Hickson Compact Groups (HCG, Verdes-Montenegro et al. 2001), than spiral galaxies that reside in the field. However, gas poor galaxies are also found in intermediate density environments (Dénes et al., 2016; Hess & Wilcots, 2013) as well as in loose groups (Kilborn et al., 2005, 2009). The H i deficiency in high and intermediate density environments is related to the gas removal mechanisms either by gravitational interactions (tidal stripping e.g. Yun et al. 1994), or hydrodynamical interactions (ram pressure stripping e.g. Gunn & Gott 1972; Rasmussen et al. 2006).
The H i deficiency is defined as the difference between the derived H i mass from a scaling relation and the observed H i mass (Haynes & Giovanelli, 1984):
[TABLE]
where is the logarithm of the derived H i mass based on the optical properties of a galaxy and is the logarithm of the measured H i mass. A galaxy is considered to be H i excess or H i deficient with deficiency values below dex or above dex (Verdes-Montenegro et al., 2001; Kilborn et al., 2009). Lutz et al. (2017) uses a slightly higher value of dex whilst Dénes et al. (2016) and Wolfinger et al. (2016) use dex to define the most H i deficient and the most H i excess, respectively.
The global H i deficiency of Choir groups, based on the total H i group mass obtained from HIPASS, was studied by Sweet et al. (2013) and they found that these groups are not significantly H i deficient. This is in contrast to Hickson Compact Groups where Verdes-Montenegro et al. (2001) obtained average H i deficiency of 0.4 dex (although HCGs are optically selected groups and morphologically more evolved). Sweet et al. (2013) concluded that the Choir groups represent an early stage of a group assembly, similar to the M81 group.
In contrast, we derive the H i deficiency for the individual Choir galaxies, comparing H i masses obtained from the resolved ATCA and VLA imaging with the presumed H i masses based on the SINGG R-band magnitudes. We obtain the H i deficiency through the scaling relation from the HyperFit fit to the SINGG galaxies, see Figure 10. The obtained scaling relation is:
[TABLE]
We find that our relation is slightlyshallower than the scaling relation of Dénes et al. (2014) which was based on the optical photometry from HOPCAT (optical counterparts for H i Parkes All Sky Survey Catalogue from Doyle et al. 2005). The difference in slope could be due to several reasons: a) An intrinsic difference between HOPCAT R-band magnitude (from SuperCosmos plates Doyle et al. (2005)) and SINGG R-band magnitudes (from digital recordings Meurer et al. (2006)); b) We are using more strict isolation criteria than Dénes et al. (2014), thus it is less contaminated than the HOPCAT sample.
We compare the optical and H i properties of Choir galaxies to the entire SINGG galaxy sample and find that the Choir galaxies are scattered within its overall distribution (see Figure 10). We adopt the H i deficiency criteria from the literature, where galaxies with DEF above 0.3 dex are H i excess, and below 0.3 dex are H i deficient. We find that 52 per cent of our galaxies have normal H i content (between and dex). We find that one Choir central galaxy out of 11 is H i excess galaxy, (DEF 0.38 dex, as seen in Figure 10). Two centrals are H i deficient, having DEF of 0.56 and 0.62 dex. We find that two satellites are H i excess, four are within the scatter and five are H i deficient. HIPASSJ0258-74:S3 has more than four times less H i than expected from the scaling relation, which makes it the most H i deficient galaxy in our sample. Analysing separately central and satellite galaxies, we obtain average deficiency values for galaxies within Choir groups: DEF 0.06 dex and DEF 0.13 dex, respectively.
The deficiency of the entire SINGG sample is shifted towards being more gas-poor with the mean deficiency of dex. The obtained result for the SINGG sample is expected since galaxies in the sample are within different environments, while isolated galaxies are “nurture-free" and thus their gas content is not influenced by the environment. Choir satellite galaxies have a similar or lower H i mass than isolated galaxies of the same absolute magnitude, with two exceptions. Conversely, Choir central galaxies have similar H i mass than isolated galaxies of the same absolute magnitude, with three exceptions. The mean deficiency of the isolated galaxies is similar to that of the Choir central galaxies, while satellites are on average more H i deficient. This difference in H i deficiency is indicating that satellites are going to be depleted first in gas-rich groups, and likely reason for H i deficiency in Choir groups is tidal stripping.
3.6 Physical processes
What physical processes could be acting in Choir groups? We find groups which are showing signs of past or ongoing gravitational interactions as well as groups without evidence of interaction. We comment on a possible gas accretion from the gas-rich mergers and gas refuelling from the cosmic web. We also determine the global stability parameter of the Choir galaxies to determine whether the angular momenta (an internal galaxy process) is responsible for keeping gas in the disc of gas-rich galaxies.
3.6.1 Galaxy interactions and gas-accretion
We visually examined deep333Preliminary surface brightness limits, based on the DECam exposure time calculator is between 25-27 mag arcsec*-2*. optical DECam (Dark Energy Camera on the CTIO Blanco 4-m telescope) g band images and we found stellar streams as possible tracers of recent accretion of gas from minor gas-rich mergers (e.g. see Figure 3). In two groups faint stellar streams are visible: HIPASSJ0400-51 has faint features with g band magnitude of 22 mag and surface brightness 25 mag arcsec*-2*; HIPASSJ1250-20 has faint features with g band magnitude of 23 mag and surface brightness 26 mag arcsec*-2*, thus it is possible that they accreted gas through past gas-rich mergers. We presented HIPASSJ1250-20 group in Section 3.1 as an example of a possible accretion onto S1 galaxy. Gas accretion from minor mergers is shown not to be significant source of refuelling the gas content in galaxies (Sancisi et al., 2008; Di Teodoro & Fraternali, 2014) since it solely cannot be enough to sustain long star formation. Thus, the largest source of accretion in theory is the ‘cold mode’ accretion (Kereš et al., 2005), where gas is accreted from the cosmic web filaments. However, the current observations are not sensitive enough to probe this extremely low H i column density 1016 cm*-2* (Popping et al., 2009) at which H i gas within filaments is expected to be detected.
In four Choir groups (HIPASSJ1026-19, HIPASSJ1059-09, HIPASSJ1250-20, HIPASSJ1408-21), tidal interactions are evident either in optical images or in H i gas distribution or both, and yet these groups do not show extreme gas depletion. Groups with evident tidal interactions have galaxies with lower gas-mass fractions which could have experienced tidal stripping in the past. The best cases for ongoing tidal stripping as a main mechanism for gas removal is in group HIPASSJ1250-20 and HIPASSJ1026-19. Since the amount of the gas is not heavily depleted we can assume that these interactions are in an early phase. We will present in depth analysis of these groups in a subsequent paper.
3.6.2 Global stability parameter q
Rather than invoking accretion as gas-replenishment, since Choir groups are in an early stage of group formation, we can ask ourselves are the Choir galaxies instead retaining their gas content due to their intrinsic properties? Obreschkow et al. (2016) show that angular momentum regulates the H i content of a galaxy, and show a relation between the mass fraction and the global stability parameter for a sample of isolated disc galaxies. Using a sample of H i eXtreme (H ix) galaxies, Lutz et al. (2017) have shown that extremely gas-rich galaxies accumulate and retain gas (building up large H i discs) due to high baryonic specific angular momentum which makes them less efficient in forming stars.
We evaluate the global stability parameter for eight Choir galaxies: six centrals and two satellites using the approximation given by the Obreschkow et al. (2016) in order to determine whether gas-rich group galaxies are governed by this property. Obreschkow et al. (2016) define the global stability parameter q as: q j/(GM), where j is the baryonic specific angular momentum of the disc and its value can be approximated as rv; is a velocity dispersion and as in Lutz et al. (2018) we assume its value of 11 km s*-1*; G is the gravitational constant and M is the baryonic mass. The baryonic mass is defined as sum of the stellar mass, H i mass and the H2 mass: M = M*⋆* + 1.35 (M + M). The factor 1.35 accounts for the helium and the H2 fraction (M) is assumed to be 4 per cent (Obreschkow et al., 2016). In order to evaluate where galaxies lie on the mass fraction versus global stability parameter diagram we need the following galaxy properties: Ks-band stellar mass, rotational velocities and the galaxy H i radius. With the available data we can calculate stability parameter for eight Choir galaxies. We derive the Ks-band magnitudes from the 2 Micron All Sky Survey (2MASS) images (Skrutskie et al., 2006), and calculate stellar masses following Wen et al. (2013). We determine the rotational velocities using the 20% widths of H i emission line profiles following the prescription of Meyer et al. (2008). The measurement of the H i radius is described in the section 3.3.
In Figure 11 we show that all gas-rich Choir galaxies follow the analytical model derived by Obreschkow et al. (2016) for isolated disc galaxies, and are in agreement with the central, star-forming, disc-dominated galaxies from the semi-analytic model of galaxy formation (DARK SAGE; Stevens et al. 2018). The Choir galaxies with high H i mass fractions have similar q values to the H ix galaxies.
In the DARK SAGE semi-analytic model galaxies are seen to lie below the f-q relation due to merger-quenching (central galaxies) and ram-pressure stripping (satellite galaxies). These effects drastically reduce f in galaxies while q is only slightly affected (Stevens et al., 2018), which would cause galaxies to move downwards on the Figure 11. Observationally, the study of the f-q relation in compact groups and clusters have not been performed yet.
We conclude that, at least partially, the gas-rich Choir galaxies are keeping their gas in the disc due to higher specific angular momenta. This result shows that in small gas-rich groups, such as Choirs, the internal galaxy properties are dominating over a group impact on gas-rich galaxies since they follow the f-q relation.
4 Discussion
The main cause of disturbances in galaxy groups are tidal interactions, and it is still an open question whether ram-pressure stripping occurs in groups (Westmeier et al., 2011; Rasmussen et al., 2012; Brown et al., 2016; Stevens & Brown, 2017). Studying individual galaxies in large loose groups, Kilborn et al. (2009) found that the H i deficient galaxies are likely to be within a projected distance of 1 Mpc from the group centre. As galaxy groups grow in number of members they are more likely to have H i deficient galaxies (Hess & Wilcots, 2013). Based on the degree of H i disturbance and deficiency in 72 systems, Verdes-Montenegro et al. (2001) proposed an evolutionary scenario for Hickson compact groups (HCG) in phases: phase one is the least evolved and shows a relatively unperturbed H i distribution; in phase two, tidal features become apparent due to galaxy interactions; and in phase three, H i is stripped from some galaxies, ending up either in the intragroup medium or forming a large envelope around the group. Combining observations from the Very Large Array and Green Bank Telescope of 22 HCGs, Borthakur et al. 2010 found that the amount of detected diffuse gas was larger in the final phase of the group evolution with respect to phase one or two. We are starting to combine the puzzle of galaxy group evolutionary stages, from young gas-rich groups that are coming together for the first time (Stierwalt et al., 2017), all the way to the old fossil groups with a large central elliptical galaxy but no obvious companions (Ponman et al., 1994; Jones et al., 2003).
Our H i selected “Choirs" groups usually consist of one or two large spiral galaxies and a number of smaller galaxies. The central galaxies have on average normal gas content, unlike that typically found in more compact groups or clusters. Here we have been considering possible scenarios that can explain their H i gas content. H i-rich central galaxies could be a result of H i-rich mergers and/or gas accretion from the intrafilament medium (Janowiecki et al., 2017; Kleiner et al., 2017; Wang et al., 2015). However, gas accretion from minor mergers is thought not to be significant, Sancisi et al. (2008) and Di Teodoro & Fraternali (2014) estimate (for spiral galaxies from the WHISP sample) a maximum gas accretion rate to be 0.1–0.2 and 0.28 M*⊙* yr*-1* respectively, which is not enough to sustain star formation. We find a number of dwarf galaxies around the large Choir spiral galaxies, however we do not detect H i gas in all these dwarf galaxies. The H i flux from the HIPASS survey yielded the total flux of each group. By a comparison of the total flux from the HIPASS and the resolved H i fluxes of each detected galaxy in a group from the ATCA, we find that the majority of the H i content is within the central galaxies.
Sweet et al. (2013) stated that eight groups from the Choir sample can be considered as Local Group (LG) analogous in terms of their morphology and magnitude. They note that in terms of compactness Choirs are more compact than the LG and in terms of H i content Choirs are more gas-rich than LG. Our resolved study of H i content is in agreement with Sweet et al. (2013). We make a hypothesis of how LG would look like if it was placed at a mean Choir distance of 87 Mpc. The sum of the H i masses of Andromeda, M33, Milky Way and neglecting all other dwarf galaxies would be 13.4109 M*⊙* (Cram et al., 1980; Carignan et al., 2006; Tavakoli, 2012; Nakanishi & Sofue, 2016) which is above the 3 detection limit in HIPASS (Sweet et al., 2013). In order to detect both large spiral galaxies, Andromeda and Milky Way within the SINGG FOV, they would have to have a certain orientation in space. The separation between them is 784 kpc (Stanek & Garnavich, 1998) which is 34 arcmin (face-on), at the distance of 87 Mpc, thus the separation is larger than SINGG FOV (15 arcmin). We would be able to see them within single a SINGG FOV if we are viewing the Milky Way - Andromeda system under angle between 55 and 85 degrees, because in this case the projected separation between the galaxies would be between 2 and 13 arcmin. Only in the stated case, the LG would be considered as a group with two central galaxies, since Milky Way and Andromeda have similar R-band magnitudes (Sweet et al., 2013). We would be able to see one smaller satellite M33 and a dwarf satellite Large Magellanic Cloud (LMC) in SINGG because they are massive enough and are star-forming galaxies (Gratier et al., 2010; Garrison-Kimmel et al., 2018), although we would probably not resolve H i in LMC with our ATCA observations. In all other cases, we would observe Milky Way and Andromeda as single galaxies, (or they would be aligned one in front of the other) and their detection would be below the detection limit of 3 in HIPASS.
We compare our gas-rich galaxies to gas-rich galaxies from the H ix sample (Lutz et al., 2018) in terms of their global stability parameter q. We find that the Choir galaxies with the higher global stability parameter, comparable to that of H ix galaxies, have stabilized the gas against the gravitational collapse which led to a build up of their large H i sizes and masses. This led us to conclude that the evolution of these Choir galaxies is primarily governed by their internal properties and not by the group environment.
With upcoming surveys, such as WALLABY (Koribalski, 2015), the number of gas-rich galaxy groups will increase and we will be able to conduct an in depth statistical study in order to understand the most common environment in the Universe.
5 Summary and Conclusions
In this paper we present an analysis of integrated H i properties for 27 galaxies within nine low mass, gas-rich, late-type dominated groups denoted “Choirs". We compare the Choir group galaxies to a sample of isolated galaxies in terms of gas-mass fraction, specific star formation rate, star formation efficiency, depletion time and H i deficiency.
We find that central Choir galaxies have on average normal gas-mass fraction with respect to the isolated galaxies of the same stellar mass. The specific star formation rates of these central galaxies are similar to the galaxies in isolation. Using the H i size-mass relation we find that the most massive galaxies have the largest H i discs and fall neatly onto the Wang et al. (2016) relation. We find that all outliers are disturbed systems that are experiencing interactions. We find that satellite Choir galaxies have on average lower gas-mass fraction with respect to the isolated galaxies of the same stellar mass, with a number of exceptions. The specific star formation rate of Choir satellite galaxies is scattered with respect to the ones in isolation and we find that the majority of the dwarf satellite galaxies (M 109.5 M*⊙*) have lower specific star formation rate than the isolated dwarf galaxies.
We comment on the physical processes that are acting in Choir groups:
(a) We find possible traces of recent or ongoing minor mergers in HIPASSJ0400-51:S1 and HIPASSJ1250-20:S1 (see Figure 3) which could be responsible for an overall increase in their H i content. We find traces of past interactions in HIPASSJ1051-17 as the S1 galaxy is a polar ring galaxy (Kilborn et al. in prep.).
(b) In four galaxy groups: HIPASSJ1026-19, HIPASSJ1059-09, HIPASSJ1250-20 and HIPASSJ1408-21, tidal interactions are evident either in optical images, H images or H i gas distribution or in some groups all of them. Tidal interactions maybe the main gas removal mechanism in galaxies HIPASSJ1026-19:S2 and HIPASSJ1250-20:S2.
(c) We do not find evidence of current tidal inderactions in HIPASSJ0205-55, HIPASSJ0258-74 and HIPASSJ2027-51.
We determine the specific angular momentum for eight Choir galaxies (see Figure 11) and find that seven of them have comparable values to the galaxies from the H ix sample (Lutz et al., 2018). It is possible that a higher specific angular momentum could explain how galaxies retain the large fraction of H i gas in their discs. We find that Choir galaxies lie on the f-q relation thus it is possible that normal secular evolution is prevailing over the environmental influences on these galaxies.
We provided further evidence that the Choirs are in the early stages of a group assembly as stated in Sweet et al. (2013). This scenario is supported by the fact that we find galaxies with average gas content and specific star formation rate comparable to those galaxies found in isolation. In groups with two large galaxies of similar stellar mass, we find that either the H i gas content is evenly distributed amongst them or the H i gas content is vastly different: one galaxy being gas-rich, other galaxy being gas-poor. In groups with one large spiral galaxy and a number of smaller galaxies, we see that the H i-poor galaxies are those of a smaller stellar mass. These are indications that the H i gas will be depleted firstly in galaxies with the lower stellar mass. The overall H i content of Choir groups is higher than that of HCGs; nevertheless, we use the proposed evolutionary scenario for HCGs by Verdes-Montenegro et al. (2001), and we see that some Choir groups presented in this paper would be placed into phase 1: groups that contain average gas content (and gas-rich galaxies) and do not exhibit signs of interaction; and early phase 2: groups that contain galaxies that interact and contain H i deficient galaxies. We find no phase 3 Choir groups.
In a subsequent paper we will present an analysis of each group individually and we will discuss outliers, galaxy interactions and the distribution of atomic and molecular gas in Choir groups. The majority of galaxies in the local Universe are within environments such as Choirs thus it is important to study this group regime in order to have a complete understanding of a galaxy evolution.
6 Acknowledgements
We thank the anonymous referee whose very detailed and constructive comments and suggestions improved this paper.
RD is supported by a Swinburne University SUPRA postgraduate scholarship.
RD would like to thank Katinka Gereb for her help with the VLA data reduction using CASA; Katharina Lutz for providing values of the global stability parameter and mass fraction for the H ix galaxies; and Chandrashekar Murugeshan for the helpful discussions.
This work benefited from the trip to KIAA conference, RD would like to thank the Astronomical Society of Australia for providing a student grant for this trip.
The Australia Telescope Compact Array is part of the Australia Telescope National Facility which is funded by the Australian Government for operation as a National Facility managed by CSIRO.
This paper includes archived data obtained through the Australia Telescope Online Archive (http://atoa.atnf.csiro.au).
Parts of this research were conducted by the Australian Research Council Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), through project number CE170100013.
The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc.
Based on observations at Cerro Tololo Inter-American Observatory, National Optical Astronomy Observatory (NOAO Prop. AAT/13A/02; PI: Sarah M. Sweet), which is operated by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation.
This research has made use of the NASA/IPAC Ex- tragalactic Database (NED), which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration.
This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation.
This research has made use of the NASA/ IPAC Infrared Science Archive, which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration.
This research has made use of the VizieR catalogue access tool, CDS, Strasbourg, France. The original description of the VizieR service was published in A&AS 143, 23.
This research has made use of python https://www.python.org and python packages: astropy (Astropy Collaboration et al., 2013), matplotlib http://matplotlib.org/ (Hunter, 2007), APLpy https://aplpy.github.io/, NumPy http://www.numpy.org/ (Van Der Walt et al., 2011) and SciPy https://www.scipy.org/ (Jones et al., 2001).
Appendix A Additional tables: properties of the Choir and isolated galaxies.
Appendix B The deprojection and measurement of the H i diameters in HIPASSJ1250-20 group
The method of measuring H i diameters is described in Section 3.3, here we show the observed and deprojected intensity maps (Figure 12 (a) and (b) respectively) of HIPASSJ1250-20 group and the results from the measurement of the H i diameters (Figure 13: (a) HIPASSJ1250-20:S1; (b) HIPASSJ1250-20). In Figure 12 and Figure 13, images have pixel scale (X and Y), where one pixel corresponds to 1 kpc. Images used here are regridded using DECam optical image as a map, this was performed in MIRIAD with the task REGRID. The measured D without the regridding has larger uncertainty, also the smaller galaxies (such as S2) due to larger beam size, have large beam corrections. The measured D using normal and regridded images are in agreement: D[S1]=105.78.7; D[S1]=106.44.7; D[S2]=24.115.2; D[S2]=29.29.1.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Argudo-Fernández et al. (2013) Argudo-Fernández M., et al., 2013, A&A , 560, A 9 · doi ↗
- 2Astropy Collaboration et al. (2013) Astropy Collaboration et al., 2013, A&A , 558, A 33 · doi ↗
- 3Audcent-Ross et al. (2018) Audcent-Ross F. M., et al., 2018, MNRAS , 480, 119 · doi ↗
- 4Baldwin et al. (1981) Baldwin J. A., Phillips M. M., Terlevich R., 1981, PASP , 93, 5 · doi ↗
- 5Barnes et al. (2001) Barnes D. G., et al., 2001, MNRAS , 322, 486 · doi ↗
- 6Barsanti et al. (2018) Barsanti S., et al., 2018, Ap J , 857, 71 · doi ↗
- 7Bekki (2008) Bekki K., 2008, MNRAS: Letters , 390, L 24 · doi ↗
- 8Berlind et al. (2006) Berlind A. A., et al., 2006, Ap JS , 167, 1 · doi ↗
