Open Access
Issue
A&A
Volume 696, April 2025
Article Number A176
Number of page(s) 37
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/202453005
Published online 18 April 2025

© The Authors 2025

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1. Introduction

Hickson Compact Groups (HCGs) are systems of typically four to ten groups of galaxies in close proximity to each other. They were first catalogued by Paul Hickson in Hickson (1982). On large scales, they are located in low-density environments (Palumbo et al. 1995). Their galaxy members are characterised by low-velocity dispersion (∼200 km s−1, Hickson et al. 1992) but are separated by small distances similar to those in the centre of clusters, making HCGs an ideal testbed to study gas and star formation quenching processes. According to the evolutionary sequence proposed by Verdes-Montenegro et al. (2001), in phase 1, galaxies in HCGs have the majority of their neutral atomic hydrogen (H I) gas in the disc of the galaxies. In phase 2, gravitational interactions between member galaxies are thought to be the main driver in disrupting gas morphology and kinematics, displacing approximately 30% to 60% of the gas from the discs and forming structures such as tidal tails and bridges. Phase 3 HCGs are divided into two subcategories. In phase 3a, most or all the gas is found in the form of intra-group clouds or tidal tails. In phase 3b, the H I in the groups forms a large cloud with a single peaked global profile and a single velocity gradient. However, Jones et al. (2023) questioned the existence of this latter phase, as subsequent studies found that only HCG 49, among the 72 HCGs analysed by Verdes-Montenegro et al. (2001), fit this classification. Therefore, Jones et al. (2023) made a minor modification to the original evolutionary sequence of HCGs proposed by Verdes-Montenegro et al. (2001). The first modification concerns the threshold at which groups are classified as being phase 2 or phase 3. Jones et al. (2023) redefined phase 2 groups as those where 25% to 75% of the H I was associated with features not related to the disc of the galaxies. If more than 75% of the H I was associated with extended features, they classified the groups as phase 3. The second modification involved the elimination of the phase 3b classification and the introduction of phase 3c. They used this new designation scheme to specify HCGs that would have been classified as phase 1 but only one galaxy among the group members was detected in H I. In fact, phase 3c groups are thought to be systems at an advanced evolutionary stage that have recently acquired a new member galaxy. In this paper, we use the modified classification scheme by Jones et al. (2023).

The H I content of HCGs has been the subject of many studies in the literature. Williams & Rood (1987) observed 51 HCGs using the NRAO 91 m and NAIC 305 m telescopes. They found that HCGs contain, on average, half as much H I as loose groups of similar optical morphology and luminosity types. Their H I mass is lower than the expected H I mass from their optical morphology and luminosity types. Huchtmeier (1997) used the Effelsberg 100-m telescope to observe 54 HCGs at the 21 cm wavelength. They detected H I for 41 HCGs. The detection rate increased with the number of spiral galaxies in the groups. Huchtmeier (1997) expanded their study by incorporating the sample from Williams & Rood (1987), resulting in a total of 75 HCGs, including 61 detections and 14 upper limits. Huchtmeier (1997) observed many HI-deficient groups by comparing their integrated H I mass with their blue luminosity. Verdes-Montenegro et al. (2001) used single-dish data for 72 HCGs to examine their overall HI contents. In addition, they used high-resolution Very Large Array (VLA) mapping for 16 of these groups to investigate their H I spatial distributions and kinematics in more details. They found that, on average, HCGs contain only 40% of the expected H I mass based on the optical luminosity and morphological types of the member galaxies. Although the more deficient groups showed a higher X-ray detection rate than the less deficient ones, suggesting the importance of phase transformation for the missing H I, the results were not statistically significant; subsequent analysis showed no significant correlation between H I deficiency and X-ray emission (Rasmussen et al. 2008). Verdes-Montenegro et al. (2001) also found that the individual member galaxies have a higher deficiency than the group as a whole, indicating the importance of gas removal processes from galaxies in HCGs. High angular resolution VLA observations showed that 70% of the spiral galaxies had a disturbed morphology and kinematics. However, they did not find a one-to-one correlation between the amount of H I mass in tidal features and H I deficiency. This may indicate a rapid evolutionary sequence once the gas is dislodged from the member galaxies. Our data will be crucial to investigating this scenario.

Our paper is motivated by the Green Bank Telescope (GBT) observations of Borthakur et al. (2010) and the VLA-based analysis of Jones et al. (2023). Borthakur et al. (2010) observed the H I spectrum of 26 HCGs, and of these, only one was undetected. They detected diffuse H I emission, spread over 1000 km s−1, which was missed by previous observations. They quantified the newly observed H I emission in terms of excess mass, which was defined as the difference between the mass measured by the GBT and the mass measured by the VLA (Mexcess = MGBT − MVLA). The excess gas mass fraction (Mexcess = MGBT − MVLA/MGBT) was found to vary between 5% and 80%, and it was higher for more H I deficient groups. The result suggested that at least part of the missing H I was in the form of a faint, diffuse H I component, spread over a large velocity range. Hence, despite the detection of the diffuse component, there is still a significant discrepancy between the observed H I mass and the expected values, and the question of where the remaining H I in HCGs is located remains unanswered.

If the rest of the gas in HCGs is heated, we would expect to detect hot intra-group X-ray gas in H I deficient groups. However, previous observations have showed that either most of the detected X-ray gas was confined within the individual member galaxies (Hess et al. 2017) or the properties of the detected hot intra-group gas among the studied H I deficient sample were very diverse and not correlated with the H I properties (Rasmussen et al. 2008). In addition, being H I deficient does not necessarily imply the presence of hot intergalactic medium (IGM) gas. For example, HCG 30 is one of the most H I deficient group, yet no hot IGM has been observed (Rasmussen et al. 2008). Previous observations also ruled out thermal evaporation as the likely dominant mechanisms to remove H I gas in HCGs (Rasmussen et al. 2008). As suggested by Hess et al. (2017), the missing H I is either too diffuse to be detected or has been ionised and is in the form of a diffuse X-ray gas that is also undetected. The increased sensitivity of new radio telescopes is promising regarding the ability to detect the diffuse diffuse H I component that might have been missed by previous instruments such as the VLA or the Westerbork Synthesis Radio Telescope (WSRT). For example, using the Five-hundred-meter Aperture Spherical Telescope (FAST), Wang et al. (2023) detected faint, extended, diffuse H I emission around the tidally interacting NGC 4631 group that was missed by the WSRT. The column density of the diffuse H I was still above the critical column density for photo-ionisation and thermal evaporation. By combining data from the Karoo Array Telescope (KAT-7), WSRT, and the Arecibo Legacy Fast ALFA (ALFALFA) survey, Hess et al. (2017) also detected faint H I emission in HCG 44, far from the group centre. Thus, we expect that part of the missing H I in HCGs is in the form of extended, diffuse H I that was missed by previous observations.

In an effort to recover the missing H I in HCGs and understand more about the survival of H I in groups, we observed six HCGs at intermediate and advanced evolutionary stages with MeerKAT. The main scientific goals of the proposal were to

  1. understand the transition of HCGs from a complex of H I tidal structures to the most extreme phase, where galaxies seem to have entirely lost their H I,

  2. determine the distance from the core at which the H I can survive in groups and assess if magnetic fields aid in this survival within the harsh intra-group medium (IGrM),

  3. examine the effects of potential encounters with intruder galaxies located outside the field of view of the VLA,

  4. investigate the previously unexplored role of intra-group gas in the accelerated transition of galaxies from an active to quiescent state.

This paper describes the processing and general description of the data; specific goals as listed above will be addressed in future papers. This paper is organised as follows. We introduce the individual groups in Section 2. We describe the observations and data reduction in Section 3. We present the results in Section 4. We summarise the findings and give our conclusion in Section 5.

Reproducibility: This paper is fully reproducible. We used the Snakemake workflow management system Mölder et al. (2021) to organise our scripts in terms of ‘rules’ that contain the names of inputs and shell commands to produce different data products, figures, and tables. A single command was then used to execute the different rules sequentially and produce the final PDF paper. Snakemake is an advanced workflow management system that facilitates reproducibility and encourages transparency in scientific research by automating the execution of complex data analysis steps. This aligns with the FAIR Principles (i.e., findable, accessible, interoperable, and reusable) and Open Science policies.

2. Sample

To have a coherent view of the evolutionary sequence of H I in HCGs, we observed six HCGs in which three are in phase 2 and three are in phase 3. We omitted groups in phase 1 due to their similarities with galaxies in isolation. We selected phase 2 groups since the mechanisms responsible for gas removal, dispersion, or depletion are mostly active during this evolutionary stage. In addition, we included phase 3 groups, as they represent the most advanced stage of evolution, allowing us to investigate the ultimate fate of the H I gas. We list the properties of the groups in Table 1 and summarise previous findings below.

Table 1.

Group parameters.

2.1. HCG 16

The H I spectrum of HCG 16 was obtained by Borthakur et al. (2010) using the GBT. Besides, HCG 16 was previously mapped in H I with the VLA using the C and D array configurations (Verdes-Montenegro et al. 2001; Jones et al. 2019). Jones et al. (2019) found a total mass of log MHI/M = 10.53 ± 0.05. This mass is higher than the mass measured by the GBT since a significant fraction of the H I in HCG 16 was beyond the reach of the GBT field of view, which was centred at the core of the group. However, by weighting the VLA map with the GBT beam response, the measured H I mass obtained by the GBT is about 5% higher than that of the VLA (Borthakur et al. 2010). Verdes-Montenegro et al. (2001) found that most of the H I in HCG 16 was in the form of extended features such as tidal tails and bridges (see also Jones et al. 2019). The group members are connected by tidal features. The most notable one is the long tidal tail with a projected extent of 160 kpc, stretching southeast towards NGC 848 from the core of the group. Additionally, an H I tail was found extending from HCG 16c to a faint optical counterpart located north of HCG 16a (Jones et al. 2019; Román et al. 2021). HCG 16 was also observed in Hα and R-band by the Survey for Ionisation in Neutral Gas Galaxies (SINGG, Meurer et al. 2006). The observations revealed that all four members of HCG 16 had prominent high surface brightness (HSB) nuclear Hα emission, with at least three members having minor-axis outflow. The deep Chandra X-ray and VLA-GMRT 1.4 GHz radio continuum study of O’Sullivan et al. (2014) uncovered the presence of a ridge of diffuse X-ray emission linking the four member galaxies. Their map also showed similarities between the hot and cold gas structures, with the diffuse X-ray ridge generally overlapping with the H I filaments, especially the brighter ridge connecting NGC 838 and NGC 839. However, the X-ray ridge connecting NGC 838 and NGC 835 did not match well with the H I filaments. O’Sullivan et al. (2014) detected H I X-ray, and radio continuum emission from the same region, indicating the presence of a multi-phase IGM in HCG 16.

2.2. HCG 30

Previous VLA observations of HCG 30 failed to detect H I in the core of HCG 30 (Jones et al. 2023). However, the GBT spectrum of Borthakur et al. (2010) showed faint, broad H I emission, spread over the whole velocity range of the group, suggesting the presence of a diffuse H I component in the group. By comparing the VLA flux with that of the GBT, Borthakur et al. (2010) found an excess gas mass fraction of 23%, corresponding to an H I excess mass of 1.45 × 108 M. Despite being the most H I deficient group (Verdes-Montenegro et al. 2001), the Chandra X-ray observations of Rasmussen et al. (2008) found no clear evidence of the presence of diffuse IGM emission. Various tests, including an examination of emission across their S2 and S3 CCDs and a search for radial gradients in emission did not reveal significant variation or clear spatial variations in diffuse emission. This indicated a consistency with the local background. In addition, previous observations failed to detect radio continuum emission in the group, and little far infrared emission (FIR) was observed in the member galaxies (Verdes-Montenegro et al. 2007). The group members also have little Hα emission (Vílchez & Iglesias-Páramo 1998). Small amounts of Hα emission were seen in the central region of the barred early-type spiral galaxy HCG 30a and HCG 30c. However, the emission in HCG 30a was contaminated by a nearby bright star, leading to a large error in the measured Hα emission. The Hα emission in the barred late-type spiral galaxy HCG 30b was found in its bulge. No Hα emission was found in HCG 30d, a lenticular galaxy. Lastly, HCG 30 may have consumed all its molecular gas through star formation as evidenced by the lack of CO detection (Verdes-Montenegro et al. 1998).

2.3. HCG 31

HCG 31 is one of the most studied compact groups in the literature due to the peculiar characteristics of its members, such as the presence of tidal tails, a merger, and prominent starbursts. All members of HCG 31 were previously detected in H I. Their H I content range from 10% to 80% of the expected values (Verdes-Montenegro et al. 2005). However, the group as a whole is not deficient in H I. The H I observations by Verdes-Montenegro et al. (2005) revealed many tidal features, containing 60% of the total HI mass in HCG 31. In addition, all member galaxies showed signs of interactions. The H I analysis by Jones et al. (2023) using VLA data indicated that most of the H I in HCG 31 was found in the IGrM, although no reliable separation of galaxies and tidal features was possible. Several members of HCG 31, namely E, F, H, and R, were considered to be TDGs, tidal dwarf candidates, or tidal debris (Richer et al. 2003; López-Sánchez et al. 2004; Mendes de Oliveira et al. 2006; Gómez-Espinoza et al. 2023). Despite the abundance of strong dynamical processes in HCG 31, only weak diffuse X-ray emission was detected in the group (Desjardins et al. 2013). The CO emission in HCG 31 is also faint, with enhanced CO emission found in the overlap region between HCG 31a and HCG 31c (Yun et al. 1997). This deficiency in CO was mainly attributed to tidal disruption in the group, rather than a consequence of low metal abundance.

2.4. HCG 90

HCG 90 was previously imaged in H I by the ATCA and VLA, with only NGC 7172 detected (Oosterloo & Iovino 1997; Jones et al. 2023), while weak H I emission in NGC 7176 was indicated by the GBT spectrum (Borthakur et al. 2010). However, all four members were detected in CO (Heckman et al. 1989; Huchtmeier & Tammann 1992; Verdes-Montenegro et al. 1997). CCD images taken with the 3.5 m New Technology Telescope at the European Southern Observatory (ESO) by White et al. (2003) showed that the three core members of HCG 90 are embedded in large diffuse optical light with surprisingly a very small amount of associated hot intra-cluster gas. The presence of the diffuse optical light was attributed to tidally stripped stars due to galaxy interactions. The amount of diffuse light in HCG 90 is unprecedented in comparison to either other observations or theoretical expectations involving multiple interacting systems. It has a narrow range of colour, consistent with an old stellar population. Miah et al. (2015) suggested that the diffuse light existed before the onset of the merger event in HCG 90. The Chandra map by Desjardins et al. (2013) showed bridges of diffuse X-ray emission connecting the three core members, as well as small common envelopes encompassing them. However, no diffuse X-ray emission associated with the IGM was found.

2.5. HCG 91

HCG 91 is heavily discussed in the literature due to the anomalous morphologies and kinematics of its members, including the presence of extended tidal tails, bridges, clumps, a double gaseous component, disturbed velocity field, and asymmetric rotation curves (Barnes & Webster 2001; Amram et al. 2003; Mendes de Oliveira et al. 2003; Vogt et al. 2015, 2016). The most pronounced tidal feature is the extended tidal tail of the Seyfert 1 galaxy, HCG 91a, pointing toward HCG 91c, and is visible both at optical (Eigenthaler et al. 2015) and radio wavelength (Amram et al. 2003; Vogt et al. 2016; Jones et al. 2023). Apart from HCG 91d, all other core members of HCG 91 have been detected previously in H I. The VLA maps by Jones et al. (2023) showed prominent H I tail extending to the east and curving north toward HCG 91c. In addition, HCG 91b and HCG 91c are connected by a weak extended H I emission. Overall, the VLA data showed that HCG 91 was moderately deficient in H I, containing 63% of the expected amount. At the time of writing, we are not aware of any diffuse X-ray emission measurement reported in the literature for this group.

2.6. HCG 97

Out of the five core members of HCG 97, only the approaching side of HCG 97b (IC5 359) was previously detected in H I in the VLA map of Jones et al. (2023), with the receding side having too low a signal-to-noise (S/N) to be included as real emission. Jones et al. (2023) suggested that HCG 97b has a disturbed morphology. This is in agreement with the LOFAR 144 MHz and VLA (1.4 GHz and 4.86 GHz) maps of Hu et al. (2023). They found a radio tail and an extended tail at the western side of HCG 97b. The radio contours are compressed at the southeastern side of the galaxy, suggesting ram-pressure stripping effects. Hu et al. (2023) used the 12m Atacama Large Millimeter Array (ALMA) and the 7m Atacama Compact Array (ACA) to trace molecular gas in HCG 97b. They detected CO emission in the disc of HCG 97 b, but not in the tails. The CO emission shows asymmetric morphology, with the southeast part showing an ‘upturn’, consistent with the optical morphology of the disc. Diffuse X-ray emission was previously detected in HCG 97 but without any associated diffuse optical light, optical tidal features or any other (optical) signs of interactions (Pildis et al. 1995). However, the outer stellar population of HCG 97a were found to be bluer than what was expected from an elliptical galaxy. In addition, its colour was consistent with the outer envelope of HCG 97d, suggesting a previous exchange of stars between the two galaxies.

3. Observations and data reduction

The data was taken between July 31, 2021, and January 02, 2022, using the full MeerKAT array and a minimum of 61 antennas (proposal id: MKT-20101). The L-band receiver (856–1711.974 MHz) was used in its 32k correlator mode. The observations were centred at 1389.1322 MHz and a channel width of 26.123 KHz (or 5.5 km s−1 at 1420.4 MHz) was used to divide the 856 MHz total bandwidth, resulting 32768 channels. The data was recorded with the four linear polarisation products of the MeerKAT antenna feeds (XX, XY, YX, and YY). For each observing run, the primary calibrator was observed for ∼8 minutes at the beginning. This was followed by a cycle of alternating observations between the gain calibrator (∼2 minutes) and the target (∼38 minutes), with the primary calibrator revisited every two hours. Polarisation calibrators were observed twice for ∼6 minutes. The total time spent per group was ∼6.25 hours (about 5.13 hours spent on source) or a total of 37.5 hours, including calibration overheads, at 8 seconds integration period for the six HCGs. Table 2 shows the calibrators used and the observing time spent on each target.

Table 2.

Observational properties.

The raw data were processed using the CARACal pipeline (Józsa et al. 2020). CARACal integrates various radio astronomical software that can be executed through a Stimela script (Makhathini 2018), enabling a user control over each step of the data processing. To produce a science-ready H I data cube, the raw data were passed through three main steps, called workers, in CARACal:

  1. flagging and cross-calibration,

  2. continuum imaging and self-calibration,

  3. continuum subtraction and line imaging.

CARACal produces various diagnostic plots that can be visually inspected to assess the quality of the data products. Our data was processed on a virtual machine of the Spanish Prototype of the Square Kilometre Array (SKA) Regional Center (espSRC, Garrido et al. 2022). The virtual machine has 24 CPUs and 186 GB RAM. The espSRC is hosted by the Institute of Astrophysics of Andalusia (IAA-CSIC). It serves as a test bed for future scientific activities with the SKA while promoting Open Science and FAIR Principles. We adopt a similar data reduction strategy as described by Serra et al. (2023). However, our observations consist of a single pointing rather than a mosaic. We describe our data reduction process below.

3.1. Flagging and cross-calibration

Briefly, flagging and cross-calibration were first applied to the calibrators before performing on-the-fly calibration and flagging on the target. The detailed steps are as follows. First, auto-correlations and shadowed antennas were flagged using the Common Astronomy Software Applications (CASA) task flagdata incorporated in CARACal. RFI was then flagged using AOFlagger (Offringa et al. 2012), which is the default option in CARACal. Next, cross-calibration was performed. For the primary calibrator, time-dependent delay calibration (K) was first carried out, followed by gain calibration (G) while applying the derived K term on the fly. After that, bandpass calibration (B) was run while applying the previous K and G terms on the fly. This sequence was repeated twice (order: KGBKGB in CARACal) before applying the primary calibrator delay and bandpass to the secondary calibrator (apply: KB in CARACal) and solving for gain amplitude and phase. Next, flagging was performed to remove spurious calibrated visibilities before solving again for the gain amplitude and phase, and additionally bootstrapping the flux scale from the primary calibrator (order: GAF and calmode: [ap, null, ap] in CARACal). For the target, on-the-fly calibration, applying the primary calibrator delays and bandpass, as well as the secondary calibrator gains, was carried out. Auto-correlations, shadowed antennas, and RFI were also flagged.

3.2. Continuum imaging and self-calibration

This step is part of the selfcal worker in CARACal to get continuum models and perform self-calibration. High signal-to-noise (S/N) continuum images and polarisation data will be described in a future paper. Here, we only describe steps relevant for H I line imaging. We use WSClean (Offringa et al. 2014) for imaging and Cubical (Kenyon et al. 2018) for self-calibration. Radio continuum emission was iteratively imaged and self-calibrated. At each step, WSClean was allowed to clean blindly down to a masking threshold set by us. At each new iteration, the previous clean mask was used to allow for a deeper cleaning. Multi-scale cleaning was used, with scales automatically selected by the algorithm. WSclean was run on a 2 deg × 2 deg image with 1″ pixel size using Briggs weighting of −1 with 0″ tapering. Three calibration iteration was used, requiring four WSClean imaging with different clean mask threshold (auto-mask option for WSClean): 30, 20, 10, and 5 times the rms noise level. The same cleaning threshold, 0.5 (auto-threshold parameter for WSClean), was used for each imaging option.

3.3. Continuum subtraction and line imaging

The continuum model generated in the previous step was subtracted from the target visibilities. Additionally, doppler-tracking corrections and a first-order polynomial fitting to the real and imaginary parts of the line-free channels were carried out with the CASA task mstransform in CARACal. Next, any erroneous visibilities produced by mstransform were flagged with AOflagger (Serra et al. 2023). After that, H I data cubes were produced with WSClean using different imaging parameters. Multi-scale cleaning was used using a major cycle clean gain of 0.2 and scales were automatically selected by WSClean. For each group, we made available cubes at 60″, as well as higher resolution cubes down to 20 kpc linear resolution, roughly corresponding to the scale length of the optical disc. Cubes at intermediate resolutions are also available. The 60″ cubes allow the study of extended emission; whereas the higher resolution cubes aid at separating galaxies from tidal tails or intra-group gas as previously done by Jones et al. (2023). To get the final data cubes for each group, we run CARACal to obtain a cube at 98″. Then, we run the Source Finding Application (SoFiA, Serra et al. 2015; Westmeier et al. 2021) outside of CARACal to get a better clean mask. After that, we run CARACal again using that mask to get our final lowest-resolution data cube. We again use SoFiA outside of CARACal to get another clean mask and use it to clean the next higher-resolution cube. We proceed in this manner until the highest resolution cube. To correct for primary beam attenuation, we used the model of Mauch et al. (2020). To convert the cubes from frequency to velocity, we used the MIRIAD task VELSW. Throughout this paper, we use the optical velocity definition. We also use the 60″ cube unless stated otherwise.

4. Results

The basic data products consist of data cubes, moment maps, and integrated spectra for the group as a whole, as well as source catalogues, cubelettes, moment maps, position velocity cuts, integrated spectra for each source detected by SoFiA, and 3D visualisation of each group. We used ViSL3D (Labadie-García et al. 2024), a new 3D visualisation software built upon X3D pathway (Vogt et al. 2016) to construct the 3D view of each cube, which can be used to facilitate the separation of different H I structures based on spatial and kinematical information. They are available at https://amiga.iaa.csic.es/x3d-menu/. For each of the groups, we provide two versions of the data cube: one showing the large-scale structures and one highlighting the central part of the group. The user can toggle between different iso-surface levels to make it easier to identify low-column density features. The technical details underlying their construction will be described in Labadie et al. (in prep.) A detailed 3D analysis of the different H I features in HCGs is beyond the scope of this paper. Figure 1 shows examples of data products for NGC 1622 at 43.8″ × 47.5″ resolution, a spiral galaxy previously detected in H I located at about 381 kpc southeast from the group centre. Below we describe the data for the individual groups. Overall, we observed even more extended tidal features in phase 2 groups than previously reported. For phase 3 groups, we did not detect any diffuse HI component, though we did identify some new high surface brightness features. Numerous galaxies were detected around the group centres, most notably in HCG 97 and HCG 91, suggesting that these groups may be embedded within larger structures. Additionally, some of the detected galaxies lack optical counterparts.

thumbnail Fig. 1.

Example SoFiA data products for NGC 1622; a spiral galaxy previously detected in H I located at about 381 kpc southeast from the group centre. From top-left to bottom-right panel: Signal-to-noise ratio of each individual pixels; column density map and R-band DeCaLS DR10 (Dey et al. 2019) optical image – the contour levels are (5.70 × 1018, 1.140 × 1019, 2.28 × 1019, 4.56 × 1019, 9.12 × 1019, 1.82 × 1020, 3.65 × 1020, 7.30 × 1020) cm−2; first moment map (velocity field); second moment map; major axis position velocity diagram; minor axis position velocity diagram; an example channel map; and integrated spectrum from a data cube where the noise has been masked. The ellipse at the bottom left corner of each plot shows the beam, 43.8″ × 47.5″.

4.1. GBT versus MeerKAT flux

To compare our flux measurements with the single dish GBT ones by Borthakur et al. (2010), we followed the procedures outlined in that paper. They used the AIPS task PATGN to model the GBT beam response as a two-dimensional Gaussian with a FWHM of 9 . 1 $ 9{{\overset{\prime}{.}}}1 $ defined by:

f ( R ) = C 3 + ( C 4 C 3 ) exp [ ( ( 0.707 / C 5 ) 2 R 2 ) ] , $$ \begin{aligned} f(R) = C3 + (C4 - C3)\ \exp [-((0.707/C5)^{2}\ R^{2})], \end{aligned} $$(1)

where R is the radius, C3 is the minimum response of the beam, which is equal to zero, C4 is the maximum response of the beam, which is equal to 1, and C5 is the width of the distribution in arc second, which in our case, is equal to 9.1 × 60 / 2.355. The output of PATGN is a two-dimensional map containing the GBT beam correction factor. For the sake of reproducibility, we wrote a python function to generate the GBT beam response using Equation (1). We multiply it with every channel of the MeerKAT cube. Then, we derived the integrated spectrum from the full area of the resulting cube. Our procedure ensures that emission outside the primary beam of the GBT is tapered to allow a reasonable comparison with the MeerKAT spectrum. We used the original, unsmoothed spectrum of Borthakur et al. (2010) and divided the observed brightness temperature by the antenna gain correction factor, 1.65 K/Jy, measured by the authors to convert it to Jy. We compare the GBT and MeerKAT integrated spectra in Figure 2. For a better visualisation, we have smoothed the spectra from both telescopes to a common velocity resolution of 20 km s−1 using a boxcar kernel. For phase 2 groups, the agreement is good, the MeerKAT and GBT spectra agree well with each other in terms of shape and recovered flux. For phase 3 groups, the agreement is less obvious. For HCG 30, MeerKAT recovers less emission around the velocities of HCG30b, HCG30c, and HCG30d. For HCG 90, the GBT and MeerKAT spectra seem to agree with each other except around 3500 km s−1 where the GBT spectrum shows positive emission. For HCG 97, the GBT spectrum shows emission spreading across a large velocity range. However, that of MeerKAT shows narrower emission. In addition, at virtually each velocity of HCG 97, GBT recovers more flux than MeerKAT. Figure 3 and Table 3 compares the flux measured by MeerKAT and GBT within the GBT beam area. For phase 2 groups, there is a good agreement between the flux recovered by the GBT and MeerKAT within this area, with the largest flux difference being 22%. For the groups in phase 3, GBT recovers over 70% more flux than MeerKAT for HCG 30 and HCG 97, while for HCG 90, MeerKAT recovers 32% more flux than GBT. This is partly due to the presence of negative fluxes in the GBT spectrum of HCG 90. Additionally, the MeerKAT spectrum around the velocity of HCG 90a peaks well above that of the GBT. The total flux recovered by MeerKAT is significantly higher than that of the GBT due to the much larger field of view of MeerKAT. We summarise the observational parameters and derived H I mass for each group in Table 4 and compare them with the corresponding VLA data from Jones et al. (2023). To calculate the noise level listed in the table, we use the noise cubes generated by SoFiA. First, we take the median value of each velocity slice across the noise cube. Then we take the median of the values from the slices. The noise has been estimated using the non-primary beam corrected data. However, the tabulated mass is calculated from the global profile of the primary beam corrected cubes using the following formula:

thumbnail Fig. 2.

GBT versus MeerKAT integrated spectrum smoothed at 20 Jy km s−1. The pink solid lines indicate the MeerKAT integrated spectra. The solid black lines indicate the GBT spectra of Borthakur et al. (2010). The vertical dotted lines indicate the velocities of the galaxies in the core of each group. The horizontal blue lines indicate zero intensity values to guide the eyes.

M HI [ M ] = 235600 × D 2 × S i Δ v , $$ \begin{aligned} M_{HI}[M_{\odot }] = 235600 \times D^2 \times \sum S_{i}\Delta v, \end{aligned} $$(2)

thumbnail Fig. 3.

GBT versus MeerKAT integrated flux measured within the GBT beam area. The solid line is a line of equality.

Table 3.

GBT versus MeerKAT recovered H I flux.

Table 4.

Observational parameters and derived H I mass.

where ∑SiΔv is the sum of the flux over all channels and within the field of view of MeerKAT in units of Jy km s−1, and D is the distance to the group. As described by Serra et al. (2023), excluding channels containing known H I line emission while performing continuum subtraction using the UVLIN option of the CASA task mstransform, implemented in CARACaL, can introduce artefacts in the form of vertical stripes in a RA-velocity slice from the data cube. To check the possible presence of RFI or continuum residual emission in channels excluded from the visibility fit, we show a plot of the RA-velocity slice at the central declination of each HCGs cube. We describe the data products for each group below.

4.2. HCG 16

We show the RA-velocity slice of HCG 16 in the first panel of Figure 4. We found no obvious RFI or continuum residuals. We plot the median noise level of the data cube as a function of velocity in the second panel of the figure. The noise does not change much across the spectral channel. We plot the integrated spectrum of the group in the right panel of Figure 4. The spectrum was derived from the primary beam corrected cube where the noise was blanked using a 3D mask generated by SoFiA. The most significant discrepancies between the VLA and MeerKAT spectra are observed at the velocities corresponding to galaxies HCG 16c and HCG 16d. HCG 16c is the group’s most active star-forming galaxy and is strongly interacting with HCG 16d (Jones et al. 2023). The difference between the VLA and MeerKAT spectra can be attributed to MeerKAT’s capability to detect more tidal features than the VLA. Although MeerKAT detected more H I features, as is explained in the next section, the total mass recovered is remarkably similar to that observed by the VLA.

thumbnail Fig. 4.

Left panel: Velocity versus right ascension of HCG 16. Middle panel: Median noise values of each RA-DEC slice of the non-primary beam corrected 60″data cube of HCG 16 as a function of velocity. The horizontal dashed line indicates the median of all the noise values from each slice. Right panel: Integrated spectra; the blue solid line indicates the MeerKAT integrated spectrum of HCG 16; the red solid line indicates the VLA integrated spectrum of the group derived by Jones et al. (2023). The vertical dotted lines indicate the velocities of the galaxies in the core of the group. The spectra have been extracted from area containing only genuine H I emission.

We show example channel maps of the central part of HCG 16 in Figure 5. The rest is available online. As also mentioned in Jones et al. (2019), the NW tail, apparent in the moment maps, is not visible in the channel maps. In contrast, the NE tail is clearly visible and appears to be dislodged gas from HCG 16c. As can be seen from the channel map at ∼3966 km s−1, the SE tail is made of gas coming from both NGC 848 and HCG 16d, extending across about seven channels. The hook feature becomes more distinct at ∼4034 km s−1. The channel maps clearly show the double sided tails of HCG 848, that is, the hook-like feature and part of the gas that makes up the SE tail.

thumbnail Fig. 5.

Example channel maps of the primary beam corrected cube of HCG 16 overlaid on DECaLS DR10 R-band optical images. Contour levels are (1.5, 2, 2.5, 3, 4, 8, 16, 32) times the median noise level in the cube (0.58 mJy beam−1). The blue colours show contour levels below 3σ; the red colours represent contour levels at 3σ, or higher. More channel maps can be found here.

In Figure 6, we show maps of the H I sources detected by SoFiA within the MeerKAT field of view. The left panels show all sources detected by SoFiA. The right panels highlight the central part of the group. In total, SoFiA detected nine surrounding members; their properties are discussed in an accompanying paper. In the top panels of the figure, we show the column density maps of HCG 16 overplotted on top of DeCaLS R-band (Dey et al. 2019) optical images. The moment one map of the group is shown in the bottom panel of Figure 6, which is plotted in a way that each individual members are scaled by different colour-scaling factor to highlight any rotational patterns. The surrounding members show differential rotation with an apparently undisturbed H I morphology. However, the integrated iso-velocity contours of the centre of the group as a whole are disturbed. The individual members may still have their rotation. The most remarkable finding here is the presence of numerous tidal features and clumps, as well as a long, continuous tidal tail connecting the core member of the group with NGC 848. In addition, a hook-like structure extends southeast of NGC 848 before curving to the northwest in direction parallel to the main group. As shown in Figure 7, all core members are linked by tails or high column density bridges. There is a continuity in velocity between NGC 848S and the curved tail. The hook appears to be a stretched gas from NGC 848 due to its encounter with the main group. The presence of numerous tails and clumps are typical of compact groups at the intermediate stage, indicating substantial gas loss into the IGM. To better visualise the tails, we show segmented position-velocity diagrams, taken from the paths shown in Figure 7, in Figure 8 as previously done by Jones et al. (2019) using VLA data. Most of the tails were previously identified by the VLA map of Jones et al. (2019) and we adopt the author’s nomenclature to name them. However, the tails are more pronounced in the MeerKAT map. In addition, the curved tail (hook) was not visible in the VLA map.

thumbnail Fig. 6.

H I moment maps of HCG 16. The left panels show all sources detected by SoFiA. The right panels show sources within the rectangular box shown on the left to better show the central part of the group. The top panels show the column density maps with contour levels of (3.1  ×  1018, 6.2  ×  1018, 1.2  ×  1018, 2.5  ×  1019, 5.0  ×  1019, 9.9  ×  1019, 2.0  ×  1020) cm−2. The contours are overlaid on DECaLS DR10 R-band optical images. The bottom panels show the moment one map. Each individual source has its own colour scaling and contour levels to highlight any rotational component.

thumbnail Fig. 7.

Left panel: MeerKAT H I surface density map of HCG 16 showing the paths (black lines) from which the segmented position-velocity diagrams shown in Figure 8 are taken. The black circles show the positions of the nodes that make up the different slices. Right panel: Same moment zero map overplotted on DECaLS optical images to highlight the core members.

thumbnail Fig. 8.

Segmented position-velocity diagrams of HCG 16 taken from the paths shown in Figure 7. Left panel: MeerKAT data from this paper. Right panel: VLA data from Jones et al. (2019). Blue contours show emission at 3× the rms noise. Dashed lines show negative contours. The vertical black lines indicate the positions of the nodes that make up the slices from which the position velocity diagrams where taken. The yellow ellipses at the bottom left corner of each panels show the half-power beam width BMAJ BMIN $ \mathrm{\sqrt{BMAJ * BMIN}} $ and 20 km s−1 velocity width.

Screenshots of the 3D visualisation of the group is shown in Figure 9. The left panel highlights high column density HI, whereas the right panel shows superposition of iso-surface corresponding to both low and high column density gas. This kind of visualisation is very crucial to identify features not associated with the disc of the member galaxies. The connection between all member galaxies is evident. A high surface brightness bridge connects the core of the group and NGC 848. The hook is composed of lower column density gas, stretching from NGC 848 before bending upward.

thumbnail Fig. 9.

Screenshots of the three-dimensional visualisation of HCG 16. The left panel shows iso-surface level highlighting the high-column density gas. The right panel showcases the low-column density H I gas. The 2D greyscale image is a DeCaLS R-band optical image of the group. The blue circles indicate the position of the member galaxies. The 2D greyscale image is a DeCaLS R-band optical image of the group. Interactive visualisations of the group can be accessed online and at https://amiga.iaa.csic.es/x3d-menu/

4.3. HCG 31

The RA-velocity slice shown in the left panel of Figure 10 shows no apparent RFI or continuum residuals. In addition, the median noise values as a function of velocity, in the middle panel of the figure, do not change much. We compare the global spectrum of HCG 31 from the VLA with that from MeerkAT in the right panel of Figure 10. It’s clear that MeerKAT recovers more H I emission than VLA at virtually all the velocity range of the group. The largest difference in terms of recovered flux corresponds to the A+C complex. Therefore, most of the new tidal features we detected is expected to be coming from the A+C complex. This is not surprising since HCG 31A and HCG 31C are known to be in the process of merging and their interactions are expected to produce many tidal features into the IGM. Separating features in HCG 31 is extremely challenging (Verdes-Montenegro et al. 2005; Jones et al. 2023) and may require complex kinematic analysis, which we will investigate in a future paper.

thumbnail Fig. 10.

Left panel: Velocity versus right ascension of HCG 31. Middle panel: Median noise values of each RA-Dec slice of the non-primary beam corrected 60″ data cube of HCG 31 as a function of velocity. The horizontal dashed line indicates the median of all the noise values from each slice. Right panel: Integrated spectra; the blue solid line indicates the MeerKAT integrated spectrum of HCG 31; the red solid line indicates VLA integrated spectrum of the group derived by Jones et al. (2023). The vertical dotted lines indicate the velocities of the galaxies in the core of the group. The spectra have been extracted from area containing only genuine H I emission.

Example channel maps of HCG 31 are shown in Figure A.1, and the rest is available online. Visual inspections indicate that the southeastern and the northwestern elongated H I is most likely associated with the A+C complex. Part of the northern extension could be coming from HCG 31Q, though. The elongation of the southeastern tail is likely a result of gas being stripped away during the interaction between HCG 31A and HCG 31C, which then subsequently extended by a flyby encounter between member G and the A+C complex. There are no known optical counterparts associated with the extended H I emission. Thus, they are not the results of interactions with members outside the central galaxies. The highest H I peak is associated with HCG 31C, which also corresponds to sites of active star forming regions.

We show the moment maps of HCG 31 in Figure A.2. The central part of HCG 31 is detected as one source although it is composed of five late-type galaxies and four tidal dwarf candidates. The surrounding members seem to have differential rotation and present no obvious signs of interactions. However, the integrated iso-velocity contours of the central part of the group as a whole is clearly disturbed, although the H I disc of the individual core members may still have rotation. The H I moment map is elongated towards the southeast and the northwest, which was not visible in previous VLA maps. The DeCaLS image shows no optical counterparts within the elongated region. We show the column density map of the central part of HCG 31 from a 15.47″× 11.86″datacube in Figure A3. We also show the previously identified tidal fragments and tidal dwarf candidates mentioned in Mendes de Oliveira et al. (2006). The A+C complex, HCG 31B, and the tidal fragments/dwarf candidates are embedded in high column density features, with the highest contour corresponding to the overlap area between A and C and also to F1. All the tidal fragments and the tidal dwarf candidates are embedded in a tail connecting HCG 31 G and the other core members.

We show screenshots of the 3D visualisation of HCG 31 in Figure A4. The left panel emphasises regions of high H I column density, while the right panel displays a superposition of iso-surfaces representing both low and high column density gas. Similar to HCG 16, significant tidal features are observed in the group as a result of the ongoing interactions and mergers among the member galaxies.

4.4. HCG 91

The RA-velocity plot of HCG 91 is shown in the first panel of Figure B1. There are subtle vertical stripes at velocities corresponding to channels excluded during UVLIN fit. The effects can be seen in the noise-velocity plot in the middle panel of Figure B1 where the noise increases at the channels where known emission is excluded while performing UVLIN fit. Unfortunately, we cannot correct for such effects. This is inherent to UVLIN as described in Serra et al. (2023). We show the global profiles of HCG 91 from the VLA and MeerKAT in the right panel of Figure B1. MeerKAT detects much more H I emission than the VLA due to its larger field of view of MeerKAT and better sensitivity. As shown in the global profiles, there are many sources beyond the velocity range of the central part of the group and some or many of them might still belong to the group. A more careful study of them will be done in an accompanying paper (Sorgho et al. in prep.).

We show example channel maps of HCG 91 in Figure B2, which contain the SE tail of HCG 91a. Additional channel maps are available online. Some of the 2σ contours of HCG 91a bend toward the north, and appear to be aligned with its optical tails. In addition, the curved H I emission that apparently connects HCG 91c and HCG 91d in projection appears as broken contours of low-column density gas coming from both HCG 91c and HCG 91b. Lastly, the channel maps show no H I emission at the velocity and location of HCG 91d. This galaxy has never been detected in H I emission before (Jones et al. 2023).

We show the moment maps of HCG 91 in Figure B3. We have detected many sources in HCG 91 but only those within 1000 km s−1 from the systemic velocity are shown here. On projection, the core members and a far-off members at the west side, LEDA 749936, are embedded in one common H I envelop. Two H I bridges connect HCG 91b and HCG 91c. The first was identified previously by Jones et al. (2023) which they separated as intra-group gas. The second one extends to the west from HCG 91c before curving toward the north to join another extended feature of HCG 91b. In addition, it seems to elongate towards LEDA 749936. Another higher S/N emission connects HCG 91c and HCG 91a. All the core members of HCG 91, including the detection slightly far west are detected as one source by SoFiA, highlighting the complex interactions between the members. We also show the velocity fields of the galaxies in the core of HCG 91 in Figure B4 from a data cube at 24.0″ × 20.2″, as well as their DeCaLS optical image along with their H I iso-velocity contours. Their optical major axis position angles are well aligned with their H I major axis position angles. However, the outer disc of HCG 91b is warped.

HCG 91a has asymmetric double-horned profiles, typical of spiral galaxies. Its main disc still has a clear rotational pattern; however, its halo is disturbed. HCG 91c also has asymmetric double-horned profiles. In addition, its main disc has retained rotation despite the interaction with other members. The blue-shifted velocity field of HCG 91b is more or less regular, except the warp mentioned earlier. However, the red-shifted part is clearly disturbed due to the interaction with HCG 91c. The H I emission in LEDA 749936 has very low S/N and its velocity field is disturbed. It has a single peaked global profile, typical of dwarf galaxies. We detect many sources around the central part of HCG 91 that seem to have regular rotation. However, a group of four strongly interacting galaxies can be seen southeast of the core members. Their H I velocity fields and morphologies are clearly disturbed. In addition, two interacting galaxies are found further east. Their systemic velocities differ by about 1500 km s−1 compared to those of the core groups.

We show screenshots of the 3D visualisation of HCG 91 in Figure B5. Again, the left panel showcases regions of high H I column density, while the right panel illustrates a composite view, overlaying both low and high column density gas distributions. The H I distribution of the group appears to be less chaotic and exhibits fewer low surface brightness features than that of HCG 16 and HCG 91. HCG 91b and HCG 91c are connected by a hook-like high surface brightness features and a low column density bridge. However, there is no H I connection between HCG 91c and HCG 91a. The overall morphology of the member galaxies are clearly disturbed, especially that of HCG 91a. It is also clear from the 3D view that no H I has been detected in HCG 91d.

4.5. HCG 30

The RA-velocity plot of HCG 30 is presented in the left panel of Figure C1, showing no apparent signs of continuum residuals or cleaning artefacts. The median noise of the data cube as a function of velocity is shown in the middle panel of the figure. The comparison of the MeerKAT and VLA global profile is shown in the right panel of Figure C1. The larger field of view of MeerKAT and the detection of new features and a member in the central of part of HCG 30 makes the global profile of MeerKAT much brighter than that of the VLA.

We present the example channel maps of HCG 30 in Figure C2. These only show the channels corresponding to the receding side of HCG 30b, the rest is presented online. Overall, the emission is faint and is mostly below the 3σ detection threshold.

HCG 30 is one of the most H I deficient groups in our sample and previous observations failed to detect any H I emission in the galaxies making up its core. However, we detect two clump-like H I emission slightly offset from the optical centres of HCG30a and HCG30c as shown in Figure C3. One clump is located at 0.86′ (15 kpc) from HCG30a, while the other one is at 0.33′ (6 kpc) from HCG30c. They could be the remains of stripped H I from HCG30a and HCG30c. The H I detection in HCG 30a remains spurious and unresolved despite being a large spiral galaxy. We show overview plots of HCG 30c in C4. These plots were generated using the SoFiA Image Pipeline (SIP) by Hess K.M.1 (Hess et al. 2024). As shown in the figure, the H I emission in HCG30c is faint and is barely resolved.

The H I emission corresponding to HCG 30b was split as three sources by SoFiA. We combined the three individual cubelettes from SoFiA using the MIRIAD task IMCOMB and rederived the moment maps, which we show in Figure C5. The moment one map presents apparent signs of rotation but the iso-velocity contours do not show the spider patterns expected from Sa galaxies. Instead, they show systematically convex curvatures. We note that the emission along the minor axis is barely resolved. The global H I profile exhibits an enhanced intensity at either side of a more subdued central peak. However, the overall S/N of the profile is low. The systemic velocity we got from SIP, 4519 km s−1, corresponds well to the one quoted by Jones et al. (2023), 4508 km s−1, from its optical redshift.

We show a screenshot of the 3D visualisation of HCG 30 in Figure C6. The 3D plot does not show a hint of detection at the location of the core members. The tentative detection we observed in the moment maps needs to be confirmed by future more sensitive observations. In addition, the surrounding galaxies do not show any signs of morphological disturbance.

4.6. HCG 90

4.6.1. Data cube

The RA-velocity plot of HCG 90 is shown in the left panel of Figure D1, indicating no obvious continuum residuals. The median noise values of the non-primary beam corrected cube as a function of velocity is presented in the middle panel of the figure. We compare the VLA and the MeerKAT global profile of HCG 90 in the right panel of Figure D1. Since the VLA observations only detected H I emission in HCG 90a, the MeerKAT global profile appears much brighter than that of the VLA as MeerKAT detects emission both at the central part and in the vicinity of HCG 90.

We show example channel maps of the central part of HCG 90 in Figure D2. More channel maps, covering a wider velocity range, can be downloaded here. The tail described previously appears as broken contours of low-column density H I, with the brightest contour found at 2271 to 2293 km s−1, at the southern optical tail. No 3-σ contour is found at the location of HCG 90b, HCG 90c, and HCG 90d.

The moment maps of HCG 90 are shown in Figure D3. The maps indicate a previously undetected H I tail in the central region of HCG 90. When plotted against a deep DECaLS image, the H I tail seems to be aligned with the optical tail. It has a total mass of 2.81 × 108 M, and is about 128 kpc long. On projection, it is difficult to assess which members the tails are associated with. However, part of the tail coincides with the location of HCG90b and HCG90d. Previous Hα kinematics by Plana et al. (1998) shows evidence of ongoing interaction between HCG 90b and HCG 90d, with HCG 90d acting as a gas provider. Spatially, the northern tail coincides more with the location of HCG 90d than that of HCG 90b. We also show in Figure D3 a position-velocity map of the tail taken from a thick slice shown at the top right panel of the figure. A continuity in velocity is clearly visible, suggesting that this is not a projection of multiple features but rather a single, coherent structure. A more compact H I structure is found further north of HCG 90c, at the southwestern side of HCG 90a. We therefore hypothesise that the observed tail is part of a more extended structure connecting HCG 90c with the other members, which may have escaped our detection, or have already been dispersed by other processes such as ionisation or star formation. HCG 90a has been detected before by the VLA, and we also clearly detect it with MeerKAT.

Figure E1 presents the 3D view of HCG 90. The tail we see in the moment map is also visible in the 3D plot and appears to be a tidal remnant of HCG 90c. This group might be a good candidate to look for diffuse H I components as the detection of this tidal fragment suggests that it is at a less advanced evolutionary stage than HCG 30 and HCG 97.

4.7. HCG 97

The RA-velocity plot of HCG 97 is shown in the left panel of Figure E2, presenting some continuum residual emission that appears as positive and negative vertical stripes. These artefacts also manifest as increased noise levels starting around 7500 km s−1 and above as shown in the middle panel of the figure. The comparison between the VLA global profile and that of MeerKAT is presented in the right panel of Figure E2. Another H I emission peak beyond the velocity covered by the VLA is detected by MeerKAT. The moment maps are shown in Figure E3. Many galaxies are detected beyond the central part of HCG 97. However, no H I emission is detected in HCG 97a, HCG 97c, HCG 97d, and HCG 97e. Only HCG 97b is detected, whose approaching side was also detected by the VLA (Jones et al. 2023). Jones et al. (2023) suggested that HCG 97b might be disturbed since they only detected one side of the galaxy. However, we do not see any H I extension in our map despite the fact that MeerKAT detected the two sides of the galaxies. The receding part of this galaxy is indeed fainter as previously mentioned by Jones et al. (2023). HCG 97b is an edge-on spiral galaxy, but its velocity field is typical of a dwarf with solid body rotation curve. This is further corroborated by the position-velocity diagram shown in Figure E3, which shows that HCG 97b has a very steep inner rotation curve, which most spirals with flat rotation curve also have (Bosma 1978). If this galaxy has a flat rotation curve, then we are only probing its inner disc, corresponding to the rising part of the rotation curve. Its outer disc might be too faint to be detected. Thus, we do not exclude the possibility that this galaxy is interacting.

Figure E5 presents the 3D view of HCG 97. No signs of detection are observed at the location of the core members. Only HCG 97b is detected. The surrounding members are not morphologically disturbed. Our observations confirm that this group is at an advanced evolutionary stage.

5. Summary and conclusion

We have presented MeerKAT H I observations of six HCGs (3 phase 2 and 3 in phase 3), providing essential data to understand the transition between the two most evolved phases in the evolutionary sequence of these compact galaxy aggregations. We aimed to detect diffuse H I gas that was apparent in previous GBT observations but missed by the VLA telescope. Our observations have revealed significantly more extended tidal features in phase 2 groups compared to those detected by the VLA. In addition, we have detected new high surface brightness features in phase 3 groups. The presence of tidal features, such as H I tails, bridges, and clumps, suggests substantial gas loss into the IGrM due to gravitational interactions. When derived within the field of view of GBT, our measured H I flux is similar to that of the GBT. However, we have not detected the diffuse H I component apparent in previous GBT spectra for phase 3 groups despite MeerKAT’s superb sensitivity. This indicates that part of the missing H I could be too diffuse to be detected by MeerKAT and the rest might be ionised. Numerous surrounding galaxies have been detected for both phase 2 and phase 3 groups, most of which are normal disc galaxies. This suggests that these groups might be embedded in larger structures. MeerKAT detected significantly more extended H I gas in Phase 2 groups than the VLA, thanks to its superior short baseline coverage. Thus, a single-dish telescope such as FAST is expected to detect H I at even larger angular scales beyond MeerKAT’s coverage. Apart from data cubes, source catalogues, and moment maps, we have released 3D visualisation of our data. This offers a more advanced approach to visually inspect different H I features, which can help in separating the complex substructures present in phase 2 groups.

Data availability

All the necessary scripts and data inputs are released along with this paper and can be found at https://github.com/ianjarog/hcg-data-paper; the instructions to download and process the data are also given there.

The channel maps of the groups, covering a wider velocity range than is presented here, are available at https://zenodo.org/records/14856489.

The 3D movies associated to Figures 9, A.4, B.5, C.6, E.1, and E.5 are available at https://www.aanda.org


Acknowledgments

Authors RI, LVM, AS, IL, CC, TW, BN, JM, SSE, JG acknowledge financial support from the grant PID2021-123930OB-C21 funded by MICIU/AEI/ 10.13039/501100011033 and by ERDF/EU, and the grant CEX2021-001131-S funded by MICIU/AEI/ 10.13039/501100011033, and the grant TED2021-130231B-I00 funded by MICIU/AEI/ 10.13039/501100011033 and by the European Union NextGenerationEU/PRTR, and acknowledge the Spanish Prototype of an SRC (espSRC) service and support funded by the Ministerio de Ciencia, Innovación y Universidades (MICIU), by the Junta de Andalucía, by the European Regional Development Funds (ERDF) and by the European Union NextGenerationEU/PRTR. The espSRC acknowledges financial support from the Agencia Estatal de Investigación (AEI) through the “Center of Excellence Severo Ochoa” award to the Instituto de Astrofísica de Andalucía (IAA-CSIC) (SEV-2017-0709) and from the grant CEX2021-001131-S funded by MICIU/AEI/ 10.13039/501100011033. Part of BN’s work was supported by the grant PTA2023-023268-I funded by MICIU/AEI/ 10.13039/501100011033 and by ESF+. IL also acknowledges financial support from PRE2021-100660 funded by MICIU/AEI /10.13039/501100011033 and by ESF+. JMS acknowledge financial support from the Spanish state agency MCIN/AEI/10.13039/501100011033 and by ’ERDF A way of making Europe’ funds through research grant PID2022-140871NB-C22. MCIN/AEI/10.13039/501100011033 has also provided additional support through the Centre of Excellence María de Maeztu’s award for the Institut de Ciències del Cosmos at the Universitat de Barcelona under contract CEX2019–000918–M. JR acknowledges financial support from the Spanish Ministry of Science and Innovation through the project PID2022-138896NB-C55. MEC acknowledges the support of an Australian Research Council Future Fellowship (Project No. FT170100273) funded by the Australian Government. TW acknowledges financial support from the grant CEX2021-001131-S funded by MICIU/AEI/ 10.13039/501100011033, from the coordination of the participation in SKA-SPAIN, funded by the Ministry of Science, Innovation and Universities (MICIU). A. del Olmo and J. Perea acknowledge financial support from the Spanish MCIU through project PID2022-140871NB-C21 by ‘ERDF A way of making Europe’, and the Severo Ochoa grant CEX2021- 515001131-S funded by MCIN/AEI/10.13039/501100011033. JMS acknowledge financial support from the Spanish state agency MCIN/AEI/10.13039/501100011033 and by ’ERDF A way of making Europe’ funds through research grant PID2022-140871NB-C22. MCIN/AEI/10.13039/501100011033 has also provided additional support through the Centre of Excellence María de Maeztu’s award for the Institut de Ciències del Cosmos at the Universitat de Barcelona under contract CEX2019-000918-M. EA and AB gratefully acknowledge support from the Centre National d’Etudes Spatiales (CNES), France. RGB acknowledges financial support from the Severo Ochoa grant CEX2021-001131- S funded by MCIN/AEI/ 10.13039/501100011033 and PID2022-141755NB-I00. JM acknowledges financial support from grant PID2023-147883NB-C21, funded by MCIU/AEI/ 10.13039/501100011033 OMS’s research is supported by the South African Research Chairs Initiative of the Departlent of Science and Technology and National Research Foundation (grant No. 81737). The MeerKAT telescope is operated by the South African Radio Astronomy Observatory, which is a facility of the National Research Foundation, an agency of the Department of Science and Innovation.

References

  1. Alatalo, K., Appleton, P. N., Lisenfeld, U., et al. 2015, ApJ, 812, 117 [NASA ADS] [CrossRef] [Google Scholar]
  2. Amram, P., Plana, H., Mendes de Oliveira, C., et al. 2003, A&A, 402, 865 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Amram, P., Mendes de Oliveira, C., Plana, H., et al. 2004, ApJ, 612, L5 [NASA ADS] [CrossRef] [Google Scholar]
  4. Amram, P., Mendes de Oliveira, C., Plana, H., et al. 2007, A&A, 471, 753 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Barnes, D. G., & Webster, R. L. 2001, MNRAS, 324, 859 [Google Scholar]
  6. Bitsakis, T., Charmandaris, V., Appleton, P. N., et al. 2014, A&A, 565, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Borthakur, S., Yun, M. S., & Verdes-Montenegro, L. 2010, ApJ, 710, 385 [NASA ADS] [CrossRef] [Google Scholar]
  8. Bosma, A. 1978, Ph.D. Thesis, Groningen Univ., The Netherlands [Google Scholar]
  9. Cluver, M. E., Appleton, P. N., Ogle, P., et al. 2013, ApJ, 765, 93 [NASA ADS] [CrossRef] [Google Scholar]
  10. Da Rocha, C., Mieske, S., Georgiev, I. Y., et al. 2011, A&A, 525, A86 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Desjardins, T. D., Gallagher, S. C., Tzanavaris, P., et al. 2013, ApJ, 763, 121 [NASA ADS] [CrossRef] [Google Scholar]
  12. Dey, A., Schlegel, D. J., Lang, D., et al. 2019, AJ, 157, 168 [Google Scholar]
  13. Eigenthaler, P., Ploeckinger, S., Verdugo, M., et al. 2015, MNRAS, 451, 2793 [NASA ADS] [CrossRef] [Google Scholar]
  14. Garrido, J., Darriba, L., Sánchez-Expósito, S., et al. 2022, J. Astron. Telescopes Instrum. Syst., 8, 011004 [Google Scholar]
  15. Gómez-Espinoza, D. A., Torres-Flores, S., Firpo, V., et al. 2023, MNRAS, 522, 2655 [Google Scholar]
  16. Heckman, T. M., Blitz, L., Wilson, A. S., et al. 1989, ApJ, 342, 735 [Google Scholar]
  17. Hess, K. M., Cluver, M. E., Yahya, S., et al. 2017, MNRAS, 464, 957 [NASA ADS] [CrossRef] [Google Scholar]
  18. Hess, K. M., Paoloserra, Boschman, L., Shen, A., & Julia, H. 2024, https://doi.org/10.5281/zenodo.14016958 [Google Scholar]
  19. Hickson, P. 1982, ApJ, 255, 382 [NASA ADS] [CrossRef] [Google Scholar]
  20. Hickson, P., Kindl, E., & Auman, J. R. 1989, ApJS, 70, 687 [NASA ADS] [CrossRef] [Google Scholar]
  21. Hickson, P., Mendes de Oliveira, C., Huchra, J. P., et al. 1992, ApJ, 399, 353 [NASA ADS] [CrossRef] [Google Scholar]
  22. Hu, D., Zajaček, M., Werner, N., et al. 2023, MNRAS., 522, 2393 [Google Scholar]
  23. Huchtmeier, W. K. 1997, A&A, 325, 473 [NASA ADS] [Google Scholar]
  24. Huchtmeier, W. K., & Tammann, G. A. 1992, A&A, 257, 455 [NASA ADS] [Google Scholar]
  25. Iglesias-Paramo, J., & Vilchez, J. M. 1997, ApJ, 479, 190 [Google Scholar]
  26. Iglesias-Páramo, J., & Vílchez, J. M. 2001, ApJ, 550, 204 [Google Scholar]
  27. Johnson, K. E., Vacca, W. D., Leitherer, C., et al. 1999, AJ, 117, 1708 [Google Scholar]
  28. Jones, M. G., Verdes-Montenegro, L., Damas-Segovia, A., et al. 2019, A&A, 632, A78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Jones, M. G., Verdes-Montenegro, L., Moldon, J., et al. 2023, A&A, 670, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Józsa, G. I. G., White, S. V., Thorat, K., et al. 2020, ASP Conf. Ser., 527, 635 [Google Scholar]
  31. Kenyon, J. S., Smirnov, O. M., Grobler, T. L., et al. 2018, MNRAS, 478, 2399 [NASA ADS] [CrossRef] [Google Scholar]
  32. Labadie-García, I., & Garrido, J., & Verdes-Montenegro, L. 2024, https://doi.org/10.5281/zenodo.13254756 [Google Scholar]
  33. Lisenfeld, U., Appleton, P. N., Cluver, M. E., et al. 2014, A&A, 570, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. López-Sánchez, Á. R., Esteban, C., & Rodríguez, M. 2004, ApJS, 153, 243 [Google Scholar]
  35. Makhathini, S. 2018, Ph.D. Thesis, Rhodes University, Drosty Rd, Grahams town, 6139, Eastern Cape, South Africa [Google Scholar]
  36. Mauch, T., Cotton, W. D., Condon, J. J., et al. 2020, ApJ, 888, 61 [Google Scholar]
  37. Mazzarella, J. M., & Boroson, T. A. 1993, ApJS, 85, 27 [Google Scholar]
  38. Mendes de Oliveira, C., Amram, P., Plana, H., et al. 2003, AJ, 126, 2635 [Google Scholar]
  39. Mendes de Oliveira, C. L., Temporin, S., Cypriano, E. S., et al. 2006, AJ, 132, 570 [CrossRef] [Google Scholar]
  40. Meurer, G. R., Hanish, D. J., Ferguson, H. C., et al. 2006, ApJS, 165, 307 [NASA ADS] [CrossRef] [Google Scholar]
  41. Miah, J. A., Sharples, R. M., & Cho, J. 2015, MNRAS, 447, 3639 [Google Scholar]
  42. Mölder, F., Jablonski, K. P., Letcher, B., et al. 2021, Sustainable data analysis with Snakemake, 10, 33 [Google Scholar]
  43. Offringa, A. R., van de Gronde, J. J., & Roerdink, J. B. T. M. 2012, A&A, 539, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Offringa, A. R., McKinley, B., Hurley-Walker, N., et al. 2014, MNRAS, 444, 606 [Google Scholar]
  45. Oosterloo, T., & Iovino, A. 1997, The Nature of Elliptical Galaxies; 2nd Stromlo Symposium, 116, 358 [Google Scholar]
  46. Ordenes-Briceño, Y., Taylor, M. A., Puzia, T. H., et al. 2016, MNRAS, 463, 1284 [CrossRef] [Google Scholar]
  47. O’Sullivan, E., Zezas, A., Vrtilek, J. M., et al. 2014, ApJ, 793, 73 [Google Scholar]
  48. Palumbo, G. G. C., Saracco, P., Hickson, P., et al. 1995, AJ, 109, 1476 [Google Scholar]
  49. Pildis, R. A., Bregman, J. N., & Schombert, J. M. 1995, AJ, 110, 1498 [Google Scholar]
  50. Plana, H., Mendes de Oliveira, C., Amram, P., et al. 1998, AJ, 116, 2123 [Google Scholar]
  51. Rasmussen, J., Ponman, T. J., Verdes-Montenegro, L., et al. 2008, MNRAS, 388, 1245 [NASA ADS] [CrossRef] [Google Scholar]
  52. Richer, M. G., Georgiev, L., Rosado, M., et al. 2003, A&A, 397, 99 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Román, J., Jones, M. G., Montes, M., et al. 2021, A&A, 649, L14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Rubin, V. C., Hunter, D. A., & Ford, W. K. 1990, ApJ, 365, 86 [Google Scholar]
  55. Serra, P., Westmeier, T., Giese, N., et al. 2015, MNRAS, 448, 1922 [Google Scholar]
  56. Serra, P., Maccagni, F. M., Kleiner, D., et al. 2023, A&A, 673, A146 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Torres-Flores, S., Mendes de Oliveira, C., Amram, P., et al. 2015, ApJ, 798, L24 [Google Scholar]
  58. Verdes-Montenegro, L., Yun, M., Perea, J., et al. 1997, Dark and Visible Matter in Galaxies and Cosmological Implications, 117, 530 [Google Scholar]
  59. Verdes-Montenegro, L., Yun, M. S., Perea, J., et al. 1998, ApJ, 497, 89 [Google Scholar]
  60. Verdes-Montenegro, L., Yun, M. S., Williams, B. A., et al. 2001, A&A, 377, 812 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Verdes-Montenegro, L., Del Olmo, A., Yun, M. S., et al. 2005, A&A, 430, 443 [CrossRef] [EDP Sciences] [Google Scholar]
  62. Verdes-Montenegro, L., Yun, M. S., Borthakur, S., et al. 2007, New Astron. Rev., 51, 87 [Google Scholar]
  63. Verdes-Montenegro, L., Vogt, F., Aubery, C., et al. 2017, Formation and Evolution of Galaxy Outskirts, 321, 241 [Google Scholar]
  64. Véron-Cetty, M.-P., & Véron, P. 2010, A&A, 518, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. Vílchez, J. M., & Iglesias-Páramo, J. 1998, ApJS, 117, 1 [Google Scholar]
  66. Vogt, F. P. A., Dopita, M. A., Borthakur, S., et al. 2015, MNRAS, 450, 2593 [NASA ADS] [CrossRef] [Google Scholar]
  67. Vogt, F. P. A., Owen, C. I., Verdes-Montenegro, L., et al. 2016, ApJ, 818, 115 [Google Scholar]
  68. Wang, J., Yang, D., Oh, S.-H., et al. 2023, ApJ, 944, 102 [NASA ADS] [CrossRef] [Google Scholar]
  69. Westmeier, T., Kitaeff, S., Pallot, D., et al. 2021, MNRAS, 506, 3962 [NASA ADS] [CrossRef] [Google Scholar]
  70. White, P. M., Bothun, G., Guerrero, M. A., et al. 2003, ApJ, 585, 739 [NASA ADS] [CrossRef] [Google Scholar]
  71. Williams, B. A., & Rood, H. J. 1987, ApJS, 63, 265 [NASA ADS] [CrossRef] [Google Scholar]
  72. Yun, M. S., Verdes-Montenegro, L., del Olmo, A., et al. 1997, ApJ, 475, L21 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Additional figures of HCG 31

A.1. Channel maps

thumbnail Fig. A.1.

Example channel maps of the primary beam corrected cube of HCG 31 overlaid on DECaLS DR10 R-band optical images. Contour levels are (1.5, 2, 2.5, 3, 6, 9, 16, 32) times the median noise level in the cube (0.71 mJy beam−1). The blue colours show contour levels below 3σ; the red colours represent contour levels at 3σ, or higher. More channel maps are available online.

Figure A.1 shows example channel maps from the primary beam corrected data cube of HCG 31, overlaid on DECaLS DR10 R-band optical images.

A.2. Moment maps

Figure A.2 shows the column density maps and moment one (velocity field) of HCG 31. The left panel highlights the large-scale structure of the group, whereas the right panel shows its central part. To better illustrate the previously identified tidal dwarf galaxies within HCG 31 (see section 4.3), we present in Figure A3 a higher-resolution view of the group’s central region, derived from a data cube at 15.47″× 11.86″.

A.3. Three-dimensional visualisation

We present in Figure A4 a screenshot of the 3D visualisation of HCG 31. The top panel displays iso-surfaces corresponding to high-column-density gas, emphasising the group’s densest structures, while the bottom panel highlights the diffuse low-column-density H I gas. Blue circles mark the positions of the member galaxies. The 2D greyscale background is a DeCaLS R-band optical image of the group. An interactive version of these 3D cubes is publicly accessible at https://amiga.iaa.csic.es/x3d-menu/.

thumbnail Fig. A.2.

H I moment maps of HCG 31. The left panels show all sources detected by SoFiA. The right panels show sources within the rectangular box shown on the left to better show the central part of the group. The top panels show the column density maps with contour levels of (4.0  ×  1018, 8.0  ×  1018, 1.6  ×  1019, 3.2  ×  1019, 6.4  ×  1019, 1.3  ×  1020, 2.6  ×  1020) cm−2. The contours are overlaid on DECaLS DR10 R-band optical images. The bottom panels show the moment one map. Each individual source has its own colour scaling and contour levels to highlight any rotational component.

thumbnail Fig. A.3.

Column density map of the central part of HCG 31 from a 15.47″× 11.86″datacube, overlaid on DECaLS filter optical images to highlight the core members. The contour levels are (5.96  ×  1019, 1.19  ×  1020, 2.38  ×  1020, 4.77  ×  1020, 9.53  ×  1020, 1.91  ×  1021, 3.20  ×  1021) cm−2.

thumbnail Fig. A.4.

Screenshots of the three-dimensional visualisation of HCG 31. The top panel shows iso-surface level highlighting the high-column density gas. The bottom panel showcases the low-column density H I gas. The blue circles indicate the position of the member galaxies. The 2D greyscale image is a DeCaLS R-band optical image of the group. Interactive visualisations of the group can be accessed online and at https://amiga.iaa.csic.es/x3d-menu/.

Appendix B: Additional figures of HCG 91

B.1. Noise properties and global profiles

The RA–velocity plot of HCG 91 is shown in the first panel of Figure B1. The middle panel shows the median noise levels in each RA–DEC slice of the non primary beam corrected data cube as a function of velocity. The horizontal dashed line indicates the global median noise. The right panel compares the MeerKAT H I integrated spectra with archival VLA measurements from (Jones et al. 2023), where the vertical dotted lines denote the systemic velocities of the galaxies in the group’s core. All spectra were extracted from regions containing only genuine H I emission.

B.2. Channel maps

Figure B2 presents the channel maps of HCG 91 overlaid on DECaLS DR10 I-band optical image. More channel maps, covering a larger velocity range, can be found here.

B.3. Moment maps

The column density maps of HCG 91 is shown in Figure B3. The left panels show all detected sources within the field of view of MeerKAT. The right panels present the central part of the group. The H I maps are overlaid on DECaLS DR10 I-band optical images.

B.4. Velocity fields of the core members of HCG 91

The velocity fields of the galaxies in the core of HCG 91 are shown in Figure B4. The greyscale images show optical I-band DeCALS DR10 data.

B.5. Three-dimensional visualisation

A screenshot of the 3D visualisation of HCG 91 is presented in Figure B5. The left panel displays iso-surfaces that highlight regions of high H I column density, while the right panel overlays areas of low and high column density. The blue circles mark the positions of the group’s member galaxies, and the 2D greyscale background is a DECaLS R-band image. An interactive version of these cubes can be found at https://amiga.iaa.csic.es/x3d-menu/.

thumbnail Fig. B.1.

Left panel: Velocity versus right ascension of HCG 91. Middle panel: Median noise values of each RA-DEC slice of the non-primary beam corrected 60″ data cube of HCG 91 as a function of velocity. The horizontal dashed line indicates the median of all the noise values from each slice. Right panel: Integrated spectra; the blue solid line indicates the MeerKAT integrated spectrum of HCG 91; the red solid line indicates VLA integrated spectrum of the group derived by (Jones et al. 2023). The vertical dotted lines indicate the velocities of the galaxies in the core of the group. The spectra have been extracted from areas containing only genuine H I emission.

thumbnail Fig. B.2.

Example channel maps of the primary beam corrected cube of HCG 91 overlaid on DECaLS DR10 I-band optical image. Contour levels are (1.5, 2, 2.5, 3, 6, 9, 16, 32) times the median noise level in the cube (0.69 mJy beam−1). The blue colours show contour levels below 3σ; the red colours represent contour levels at 3σ, or higher. Additional channel maps can be downloaded here.

thumbnail Fig. B.3.

H I moment maps of HCG 91. The left panels show all sources detected by SoFiA. The right panels show sources within the rectangular box shown on the left to better show the central part of the group. The top panels show the column density maps with contour levels of (4.6  ×  1018, 9.2  ×  1018, 1.8  ×  1019, 3.7  ×  1019, 7.4  ×  1019, 1.5  ×  1020, 3.0  ×  1020) cm−2. The contours are overlaid on DECaLS DR10 I-band optical images. The bottom panels show the moment one map. Each individual source has its own colour scaling and contour levels to highlight any rotational component.

thumbnail Fig. B.4.

H I velocity fields of the galaxies in the core of HCG 91. The greyscale images show DeCaLS DR10 I-band optical data. The crosses indicate the optical centre, whereas the green contours highlight the optical disc.

thumbnail Fig. B.5.

Screenshots of the three-dimensional visualisation of HCG 91. The left panel shows iso-surface level highlighting the high-column density gas. The right panel showcases the low-column density H I gas. The blue circles indicate the position of the member galaxies. The 2D greyscale image is a DeCaLS R-band optical image of the group. Interactive visualisations of the group can be accessed online and at https://amiga.iaa.csic.es/x3d-menu/.

Appendix C: Extra figures of HCG 30

C.1. Noise properties and global profiles

Figure C1 illustrates the noise properties in HCG 30 and compares the MeerKAT and VLA integrated spectra of the group. The left panel displays the RA-velocity plot. The middle panel shows the median noise levels across individual RA–DEC slices of the non primary beam corrected data cube, with a horizontal dashed line indicating the global median noise. The right panel compares the H I integrated spectrum derived from the MeerKAT data (blue solid line) with the VLA measurement from Jones et al. (2023) (red solid line). The vertical dotted lines mark the systemic velocities of the core galaxies. The spectra were extracted from regions containing only genuine H I emission.

C.2. Channel maps

Figure C2 shows example channel maps from the primary beam corrected data cube of HCG 30. Additional channel maps can be retrieved from https://zenodo.org/records/14856489.

C.3. Moment maps

Figure C3 presents the H I moment maps of HCG 30. The left panels show all H I sources detected by SoFiA within the MeerKA field of view, while the right panels display zoomed-in views of the group’s core. The top panels show column density maps overlaid on DECaLS DR10 R-band optical images. The bottom panels display the moment-one (velocity field) map for each detected source, highlighting the rotational signatures of the surrounding galaxies.

C.4. The core members of HCG 30

Our tentative detection of H I in the core members of HCG 30 is illustrated in Figure C4 and Figure C5. These plots were generated by the SoFiA Image Pipeline2(Hess et al. 2024), an automated tool that takes as an input a SoFiA-generated source catalogue and produces overview plots of the H I properties of the sources. Figure C4 and Figure C5 show H I contours overlaid on DeCALS false-colour images, H I moment maps, pixel-by-pixel signal-to-noise ratio maps, position–velocity diagrams, and global H I spectra.

C.5. Three-dimensional visualisation

Figure C6 presents a screenshot of the 3D visualisation of HCG 30. The member galaxies are indicated by the blue circles, and the greyscale background is a DECaLS R-band optical image. An interactive 3D version of these cubes can be found at https://amiga.iaa.csic.es/x3d-menu/.

thumbnail Fig. C.1.

Left panel: Velocity versus right ascension of HCG 30. Middle panel: Median noise values of each RA-DEC slice of the non-primary beam corrected 60″ data cube of HCG 30 as a function of velocity. The horizontal dashed line indicates the median of all the noise values from each slice. Right panel: Integrated spectra; the blue solid line indicates the MeerKAT integrated spectrum of HCG 30; the red solid line indicates VLA integrated spectrum of the group derived by (Jones et al. 2023). The vertical dotted lines indicate the velocities of the galaxies in the core of the group. The spectra have been extracted from areas containing only genuine H I emission.

thumbnail Fig. C.2.

Example channel maps of the primary beam corrected cube of HCG 30 overlaid on DECaLS DR10 R-band optical image. Contour levels are (1.5, 2, 3, 4, 5, 6) times the median noise level in the cube (0.59 mJy beam−1). The blue colours show contour levels below 3σ; the red colours represent contour levels at 3σ, or higher. Additional channel maps can be accessed online.

thumbnail Fig. C.3.

H I moment maps of HCG 30. The left panels show all sources detected by SoFiA. The right panels show sources within the rectangular box shown on the left to better show the central part of the group. The top panels show the column density maps with contour levels of (6.0  ×  1018, 1.2  ×  1019, 2.4  ×  1019, 4.8  ×  1019, 9.6  ×  1019, 1.9  ×  1020) cm−2. The contours are overlaid on DECaLS DR10 R-band optical images. The bottom panels show the moment one map. Each individual source has its own colour scaling and contour levels to highlight any rotational component.

thumbnail Fig. C.4.

Top left: H I column density map overlaid on DECaLS optical image of HCG 30c. The contour levels are (2.10, 4.20, 8.41) × 1018 cm−2. Top centre: Moment-1 map. The arrows indicate the slices from which the position-velocity diagrams shown at the bottom panels were derived. Top right: Signal-to-noise ratio map. Bottom left: Global H I profile. Bottom centre: Major axis position-velocity diagram. Bottom right: Minor axis position-velocity diagram.

thumbnail Fig. C.5.

Top left: H I column density map overlaid on DECaLS optical image of HCG 30b. The contour levels are (0.29, 0.57, 1.14, 2.28, 4.57) × 1019 cm−2. Top centre: Moment-1 map, the contour levels are (4384, 4444, 4504, 4564, 4624, 4684, 4744) km s−1. The arrows indicate the slices from which the position-velocity diagrams shown at the bottom panels were derived. Top right: Signal-to-noise ratio map. Bottom left: Global H I profile. Bottom centre: Major axis position-velocity diagram. Bottom right: Minor axis position-velocity diagram.

thumbnail Fig. C.6.

A screenshot of the three-dimensional visualisation of HCG 30. The blue circles indicate the position of the member galaxies. The 2D greyscale image is a DeCaLS R-band optical image of the group. Interactive visualisations of the group can be accessed online and at https://amiga.iaa.csic.es/x3d-menu/.

Appendix D: Supplementary figures of HCG 90

D.1. Noise properties and global profiles

Figure D1 shows the noise variations and global profiles of HCG 90. The left panel shows a RA-velocity slice. The middle panel plots the median noise levels per RA–DEC slice. The horizontal dashed line marks the overall median noise level. The right panel presents the MeerKAT integrated spectrum (blue solid line) compared with the VLA data (red solid line) originally presented by Jones et al. (2023). The vertical dotted lines indicate the systemic velocities of the galaxies in the core of the group.

thumbnail Fig. D.1.

Left panel: Velocity versus right ascension of HCG 90. Middle panel: Median noise values of each RA-DEC slice of the non-primary beam corrected 60″ data cube of HCG 90 as a function of velocity. The horizontal dashed line indicates the median of all the noise values from each slice. Right panel: Integrated spectra; the blue solid line indicates the MeerKAT integrated spectrum of HCG 90; the red solid line indicates VLA integrated spectrum of the group derived by (Jones et al. 2023). The vertical dotted lines indicate the velocities of the galaxies in the core of the group. The spectra have been extracted from areas containing only genuine H I emission.

D.2. Channel maps

Figure D2 displays representative channel maps from the primary beam corrected data cube of HCG 90, superimposed on enhanced DECaLS optical images. To enhance faint optical features, G- and R-band images were combined, and pixel scaling was adjusted. The rest of the channel maps can be found here.

thumbnail Fig. D.2.

Example channel maps of the primary beam corrected cube of HCG 90 overlaid on enhanced DECaLS optical images. To improve the signal-to-noise ratio of the optical data and highlight the optical tail, we have added G-band and R-band images and aligned the pixel scale to 8×8 (0.27×8 arcsec). Contour levels are (1.5, 2, 3, 4, 5, 6) times the median noise level in the cube (0.66 mJy beam−1). The yellow and green colours show contour levels below 3σ. The red colours represent contour levels at 3σ or higher. Additional figures can be downloaded here.

D.3. Moment maps

Figure D3 shows the column density and moment one (velocity field) maps of HCG 90. The left panels present all sources detected by SoFiA, whereas the right panels focus on the central region. The top panels overlay H I column density contours on DECaLS optical images. The bottom panels display the velocity field with individual colour scalings for each source to highlight the rotational patterns of the surrounding galaxies.

thumbnail Fig. D.3.

H I moment maps of HCG 90. The left panels show all sources detected by SoFiA. The right panels show sources within the rectangular box shown on the left to better show the central part of the group. The top panels show the column density maps with contour levels of (3.9  ×  1018, 7.7  ×  1018, 1.5  ×  1019, 3.1  ×  1019, 6.2  ×  1019, 1.2  ×  1020, 2.5  ×  1020) cm−2. The contours are overlaid on DECaLS optical images. The bottom panels show the moment one map. Each individual source has its own colour scaling and contour levels to highlight any rotational component.

D.4. H I tails in HCG 90

Figure D4 presents an overview plots of the H I tails in HCG 90. The top-left panel displays the HI column density map, with the blue region marking the slice used to derive the position-velocity diagram in the bottom-left panel of the figure. The top-centre panel shows the velocity field (moment-1 map), and the top-right panel presents the signal-to-noise ratio map. The bottom panels show a position-velocity diagram and a (smoothed) global profile of the tails.

thumbnail Fig. D.4.

Top left: H I column density map of the tails of HCG 90. The contour levels are (0.35, 0.69, 1.38, 2.77) × 1019 cm−2. The blue area indicates the slice from which the position-velocity diagram shown at the bottom panel was derived. Top centre: Moment-1 map, the contour levels are (2300, 2320, 2340, 2360, 2380, 2400, 2420, 2440, 2460, 2480, 2500, 2520) km s−1. Top right: Signal-to-noise ratio map. Pixels above 6-σ are all shown in red. The black contour represent the lowest column density contour plotted at the top left panel of the figure. Bottom left: Position-velocity diagram taken from the slice shown at the top right panel. Bottom right: Global H I profile. The black line is a smoothed version of the profile using at 20 km s−1 velocity resolution.

D.5. Three-dimensional visualisation

Figure E1 presents a screenshot of the 3D visualisation of HCG 90, where the blue circles mark the positions of the member galaxies. The background image is a DeCaLS R-band optical image of the group. An interactive version of these 3D cubes is available online at https://amiga.iaa.csic.es/x3d-menu/.

Appendix E: Additional figures of HCG 97

E.1. Noise properties and global profiles

Figure E2 shows the RA-velocity plot of HCG 97 (left), the median noise level across the velocity channels (middle), and the MeerKAT and VLA integrated spectra (right).

thumbnail Fig. E.1.

A screenshot of the three-dimensional visualisation of HCG 90. The blue circles indicate the position of the member galaxies. The 2D greyscale image is a DeCaLS R-band optical image of the group. Interactive visualisations of the group can be accessed online and at https://amiga.iaa.csic.es/x3d-menu/.

thumbnail Fig. E.2.

Left panel: Velocity versus right ascension of HCG 97. Middle panel: Median noise values of each RA-DEC slice of the non-primary beam corrected 60″ data cube of HCG 97 as a function of velocity. The horizontal dashed line indicates the median of all the noise values from each slice. Right panel: Integrated spectra; the blue solid line indicates the MeerKAT integrated spectrum of HCG 97; the red solid line indicates VLA integrated spectrum of the group derived by (Jones et al. 2023). The vertical dotted lines indicate the velocities of the galaxies in the core of the group. The spectra have been extracted from areas containing only genuine H I emission.

E.2. Moment maps

Figure E3 presents the column density maps (top panels) and moment one (bottom panels) maps of HCG 97. The left panels show all sources detected by SoFiA within the MeerKAT field of view, whereas the right panels present zoomed-in views of the group’s centre. The greyscale images are DECaLS R-band optical data.

thumbnail Fig. E.3.

H I moment maps of HCG 97. The left panels show all sources detected by SoFiA. The right panels show sources within the rectangular box shown on the left to better show the central part of the group. The top panels represent the column density map with contour levels of (3.5  ×  1018, 7.0  ×  1018, 1.4  ×  1019, 2.8  ×  1019, 5.6  ×  1019, 1.1  ×  1020, 2.2  ×  1020) cm−2. The contours are overlaid on DECaLS R-band optical images. The bottom panels show the moment one map. Each individual source has its own colour scaling and contour levels to highlight any rotational component.

E.3. H I in HCG 97b

Overview plots of the H I in HCG 97b are shown in Figure E4. These images were generated by SIP, showing moment maps, signal-to-noise ratio, position-velocity diagrams and an integrated spectrum.

E.4. Three-dimensional visualisation

Figure E5 shows a screenshot of the 3D visualisation of HCG 97. The blue circles mark the member galaxies and the 2D greyscale background image is a DECaLS R-band optical data. An interactive version of the data cubes is also accessible at https://amiga.iaa.csic.es/x3d-menu/.

thumbnail Fig. E.4.

Top left: H I column density map overlaid on DECaLS optical image of HCG 97b. The contour levels are (2.10, 4.20, 8.41) × 1018 cm−2. Top centre: Moment-1 map. The arrows indicate the slices from which the position-velocity diagrams shown at the bottom panels were derived. Top right: Signal-to-noise ratio map. Bottom left: Global H I profile. Bottom centre: Major axis position-velocity diagram. Bottom right: Minor axis position-velocity diagram.

thumbnail Fig. E.5.

A screenshot of the three-dimensional visualisation of HCG 97. The blue circles indicate the position of the member galaxies. The 2D greyscale image is a DeCaLS R-band optical image of the group. Interactive visualisations of the group can be accessed online and at https://amiga.iaa.csic.es/x3d-menu/.

All Tables

Table 1.

Group parameters.

Table 2.

Observational properties.

Table 3.

GBT versus MeerKAT recovered H I flux.

Table 4.

Observational parameters and derived H I mass.

All Figures

thumbnail Fig. 1.

Example SoFiA data products for NGC 1622; a spiral galaxy previously detected in H I located at about 381 kpc southeast from the group centre. From top-left to bottom-right panel: Signal-to-noise ratio of each individual pixels; column density map and R-band DeCaLS DR10 (Dey et al. 2019) optical image – the contour levels are (5.70 × 1018, 1.140 × 1019, 2.28 × 1019, 4.56 × 1019, 9.12 × 1019, 1.82 × 1020, 3.65 × 1020, 7.30 × 1020) cm−2; first moment map (velocity field); second moment map; major axis position velocity diagram; minor axis position velocity diagram; an example channel map; and integrated spectrum from a data cube where the noise has been masked. The ellipse at the bottom left corner of each plot shows the beam, 43.8″ × 47.5″.

In the text
thumbnail Fig. 2.

GBT versus MeerKAT integrated spectrum smoothed at 20 Jy km s−1. The pink solid lines indicate the MeerKAT integrated spectra. The solid black lines indicate the GBT spectra of Borthakur et al. (2010). The vertical dotted lines indicate the velocities of the galaxies in the core of each group. The horizontal blue lines indicate zero intensity values to guide the eyes.

In the text
thumbnail Fig. 3.

GBT versus MeerKAT integrated flux measured within the GBT beam area. The solid line is a line of equality.

In the text
thumbnail Fig. 4.

Left panel: Velocity versus right ascension of HCG 16. Middle panel: Median noise values of each RA-DEC slice of the non-primary beam corrected 60″data cube of HCG 16 as a function of velocity. The horizontal dashed line indicates the median of all the noise values from each slice. Right panel: Integrated spectra; the blue solid line indicates the MeerKAT integrated spectrum of HCG 16; the red solid line indicates the VLA integrated spectrum of the group derived by Jones et al. (2023). The vertical dotted lines indicate the velocities of the galaxies in the core of the group. The spectra have been extracted from area containing only genuine H I emission.

In the text
thumbnail Fig. 5.

Example channel maps of the primary beam corrected cube of HCG 16 overlaid on DECaLS DR10 R-band optical images. Contour levels are (1.5, 2, 2.5, 3, 4, 8, 16, 32) times the median noise level in the cube (0.58 mJy beam−1). The blue colours show contour levels below 3σ; the red colours represent contour levels at 3σ, or higher. More channel maps can be found here.

In the text
thumbnail Fig. 6.

H I moment maps of HCG 16. The left panels show all sources detected by SoFiA. The right panels show sources within the rectangular box shown on the left to better show the central part of the group. The top panels show the column density maps with contour levels of (3.1  ×  1018, 6.2  ×  1018, 1.2  ×  1018, 2.5  ×  1019, 5.0  ×  1019, 9.9  ×  1019, 2.0  ×  1020) cm−2. The contours are overlaid on DECaLS DR10 R-band optical images. The bottom panels show the moment one map. Each individual source has its own colour scaling and contour levels to highlight any rotational component.

In the text
thumbnail Fig. 7.

Left panel: MeerKAT H I surface density map of HCG 16 showing the paths (black lines) from which the segmented position-velocity diagrams shown in Figure 8 are taken. The black circles show the positions of the nodes that make up the different slices. Right panel: Same moment zero map overplotted on DECaLS optical images to highlight the core members.

In the text
thumbnail Fig. 8.

Segmented position-velocity diagrams of HCG 16 taken from the paths shown in Figure 7. Left panel: MeerKAT data from this paper. Right panel: VLA data from Jones et al. (2019). Blue contours show emission at 3× the rms noise. Dashed lines show negative contours. The vertical black lines indicate the positions of the nodes that make up the slices from which the position velocity diagrams where taken. The yellow ellipses at the bottom left corner of each panels show the half-power beam width BMAJ BMIN $ \mathrm{\sqrt{BMAJ * BMIN}} $ and 20 km s−1 velocity width.

In the text
thumbnail Fig. 9.

Screenshots of the three-dimensional visualisation of HCG 16. The left panel shows iso-surface level highlighting the high-column density gas. The right panel showcases the low-column density H I gas. The 2D greyscale image is a DeCaLS R-band optical image of the group. The blue circles indicate the position of the member galaxies. The 2D greyscale image is a DeCaLS R-band optical image of the group. Interactive visualisations of the group can be accessed online and at https://amiga.iaa.csic.es/x3d-menu/

In the text
thumbnail Fig. 10.

Left panel: Velocity versus right ascension of HCG 31. Middle panel: Median noise values of each RA-Dec slice of the non-primary beam corrected 60″ data cube of HCG 31 as a function of velocity. The horizontal dashed line indicates the median of all the noise values from each slice. Right panel: Integrated spectra; the blue solid line indicates the MeerKAT integrated spectrum of HCG 31; the red solid line indicates VLA integrated spectrum of the group derived by Jones et al. (2023). The vertical dotted lines indicate the velocities of the galaxies in the core of the group. The spectra have been extracted from area containing only genuine H I emission.

In the text
thumbnail Fig. A.1.

Example channel maps of the primary beam corrected cube of HCG 31 overlaid on DECaLS DR10 R-band optical images. Contour levels are (1.5, 2, 2.5, 3, 6, 9, 16, 32) times the median noise level in the cube (0.71 mJy beam−1). The blue colours show contour levels below 3σ; the red colours represent contour levels at 3σ, or higher. More channel maps are available online.

In the text
thumbnail Fig. A.2.

H I moment maps of HCG 31. The left panels show all sources detected by SoFiA. The right panels show sources within the rectangular box shown on the left to better show the central part of the group. The top panels show the column density maps with contour levels of (4.0  ×  1018, 8.0  ×  1018, 1.6  ×  1019, 3.2  ×  1019, 6.4  ×  1019, 1.3  ×  1020, 2.6  ×  1020) cm−2. The contours are overlaid on DECaLS DR10 R-band optical images. The bottom panels show the moment one map. Each individual source has its own colour scaling and contour levels to highlight any rotational component.

In the text
thumbnail Fig. A.3.

Column density map of the central part of HCG 31 from a 15.47″× 11.86″datacube, overlaid on DECaLS filter optical images to highlight the core members. The contour levels are (5.96  ×  1019, 1.19  ×  1020, 2.38  ×  1020, 4.77  ×  1020, 9.53  ×  1020, 1.91  ×  1021, 3.20  ×  1021) cm−2.

In the text
thumbnail Fig. A.4.

Screenshots of the three-dimensional visualisation of HCG 31. The top panel shows iso-surface level highlighting the high-column density gas. The bottom panel showcases the low-column density H I gas. The blue circles indicate the position of the member galaxies. The 2D greyscale image is a DeCaLS R-band optical image of the group. Interactive visualisations of the group can be accessed online and at https://amiga.iaa.csic.es/x3d-menu/.

In the text
thumbnail Fig. B.1.

Left panel: Velocity versus right ascension of HCG 91. Middle panel: Median noise values of each RA-DEC slice of the non-primary beam corrected 60″ data cube of HCG 91 as a function of velocity. The horizontal dashed line indicates the median of all the noise values from each slice. Right panel: Integrated spectra; the blue solid line indicates the MeerKAT integrated spectrum of HCG 91; the red solid line indicates VLA integrated spectrum of the group derived by (Jones et al. 2023). The vertical dotted lines indicate the velocities of the galaxies in the core of the group. The spectra have been extracted from areas containing only genuine H I emission.

In the text
thumbnail Fig. B.2.

Example channel maps of the primary beam corrected cube of HCG 91 overlaid on DECaLS DR10 I-band optical image. Contour levels are (1.5, 2, 2.5, 3, 6, 9, 16, 32) times the median noise level in the cube (0.69 mJy beam−1). The blue colours show contour levels below 3σ; the red colours represent contour levels at 3σ, or higher. Additional channel maps can be downloaded here.

In the text
thumbnail Fig. B.3.

H I moment maps of HCG 91. The left panels show all sources detected by SoFiA. The right panels show sources within the rectangular box shown on the left to better show the central part of the group. The top panels show the column density maps with contour levels of (4.6  ×  1018, 9.2  ×  1018, 1.8  ×  1019, 3.7  ×  1019, 7.4  ×  1019, 1.5  ×  1020, 3.0  ×  1020) cm−2. The contours are overlaid on DECaLS DR10 I-band optical images. The bottom panels show the moment one map. Each individual source has its own colour scaling and contour levels to highlight any rotational component.

In the text
thumbnail Fig. B.4.

H I velocity fields of the galaxies in the core of HCG 91. The greyscale images show DeCaLS DR10 I-band optical data. The crosses indicate the optical centre, whereas the green contours highlight the optical disc.

In the text
thumbnail Fig. B.5.

Screenshots of the three-dimensional visualisation of HCG 91. The left panel shows iso-surface level highlighting the high-column density gas. The right panel showcases the low-column density H I gas. The blue circles indicate the position of the member galaxies. The 2D greyscale image is a DeCaLS R-band optical image of the group. Interactive visualisations of the group can be accessed online and at https://amiga.iaa.csic.es/x3d-menu/.

In the text
thumbnail Fig. C.1.

Left panel: Velocity versus right ascension of HCG 30. Middle panel: Median noise values of each RA-DEC slice of the non-primary beam corrected 60″ data cube of HCG 30 as a function of velocity. The horizontal dashed line indicates the median of all the noise values from each slice. Right panel: Integrated spectra; the blue solid line indicates the MeerKAT integrated spectrum of HCG 30; the red solid line indicates VLA integrated spectrum of the group derived by (Jones et al. 2023). The vertical dotted lines indicate the velocities of the galaxies in the core of the group. The spectra have been extracted from areas containing only genuine H I emission.

In the text
thumbnail Fig. C.2.

Example channel maps of the primary beam corrected cube of HCG 30 overlaid on DECaLS DR10 R-band optical image. Contour levels are (1.5, 2, 3, 4, 5, 6) times the median noise level in the cube (0.59 mJy beam−1). The blue colours show contour levels below 3σ; the red colours represent contour levels at 3σ, or higher. Additional channel maps can be accessed online.

In the text
thumbnail Fig. C.3.

H I moment maps of HCG 30. The left panels show all sources detected by SoFiA. The right panels show sources within the rectangular box shown on the left to better show the central part of the group. The top panels show the column density maps with contour levels of (6.0  ×  1018, 1.2  ×  1019, 2.4  ×  1019, 4.8  ×  1019, 9.6  ×  1019, 1.9  ×  1020) cm−2. The contours are overlaid on DECaLS DR10 R-band optical images. The bottom panels show the moment one map. Each individual source has its own colour scaling and contour levels to highlight any rotational component.

In the text
thumbnail Fig. C.4.

Top left: H I column density map overlaid on DECaLS optical image of HCG 30c. The contour levels are (2.10, 4.20, 8.41) × 1018 cm−2. Top centre: Moment-1 map. The arrows indicate the slices from which the position-velocity diagrams shown at the bottom panels were derived. Top right: Signal-to-noise ratio map. Bottom left: Global H I profile. Bottom centre: Major axis position-velocity diagram. Bottom right: Minor axis position-velocity diagram.

In the text
thumbnail Fig. C.5.

Top left: H I column density map overlaid on DECaLS optical image of HCG 30b. The contour levels are (0.29, 0.57, 1.14, 2.28, 4.57) × 1019 cm−2. Top centre: Moment-1 map, the contour levels are (4384, 4444, 4504, 4564, 4624, 4684, 4744) km s−1. The arrows indicate the slices from which the position-velocity diagrams shown at the bottom panels were derived. Top right: Signal-to-noise ratio map. Bottom left: Global H I profile. Bottom centre: Major axis position-velocity diagram. Bottom right: Minor axis position-velocity diagram.

In the text
thumbnail Fig. C.6.

A screenshot of the three-dimensional visualisation of HCG 30. The blue circles indicate the position of the member galaxies. The 2D greyscale image is a DeCaLS R-band optical image of the group. Interactive visualisations of the group can be accessed online and at https://amiga.iaa.csic.es/x3d-menu/.

In the text
thumbnail Fig. D.1.

Left panel: Velocity versus right ascension of HCG 90. Middle panel: Median noise values of each RA-DEC slice of the non-primary beam corrected 60″ data cube of HCG 90 as a function of velocity. The horizontal dashed line indicates the median of all the noise values from each slice. Right panel: Integrated spectra; the blue solid line indicates the MeerKAT integrated spectrum of HCG 90; the red solid line indicates VLA integrated spectrum of the group derived by (Jones et al. 2023). The vertical dotted lines indicate the velocities of the galaxies in the core of the group. The spectra have been extracted from areas containing only genuine H I emission.

In the text
thumbnail Fig. D.2.

Example channel maps of the primary beam corrected cube of HCG 90 overlaid on enhanced DECaLS optical images. To improve the signal-to-noise ratio of the optical data and highlight the optical tail, we have added G-band and R-band images and aligned the pixel scale to 8×8 (0.27×8 arcsec). Contour levels are (1.5, 2, 3, 4, 5, 6) times the median noise level in the cube (0.66 mJy beam−1). The yellow and green colours show contour levels below 3σ. The red colours represent contour levels at 3σ or higher. Additional figures can be downloaded here.

In the text
thumbnail Fig. D.3.

H I moment maps of HCG 90. The left panels show all sources detected by SoFiA. The right panels show sources within the rectangular box shown on the left to better show the central part of the group. The top panels show the column density maps with contour levels of (3.9  ×  1018, 7.7  ×  1018, 1.5  ×  1019, 3.1  ×  1019, 6.2  ×  1019, 1.2  ×  1020, 2.5  ×  1020) cm−2. The contours are overlaid on DECaLS optical images. The bottom panels show the moment one map. Each individual source has its own colour scaling and contour levels to highlight any rotational component.

In the text
thumbnail Fig. D.4.

Top left: H I column density map of the tails of HCG 90. The contour levels are (0.35, 0.69, 1.38, 2.77) × 1019 cm−2. The blue area indicates the slice from which the position-velocity diagram shown at the bottom panel was derived. Top centre: Moment-1 map, the contour levels are (2300, 2320, 2340, 2360, 2380, 2400, 2420, 2440, 2460, 2480, 2500, 2520) km s−1. Top right: Signal-to-noise ratio map. Pixels above 6-σ are all shown in red. The black contour represent the lowest column density contour plotted at the top left panel of the figure. Bottom left: Position-velocity diagram taken from the slice shown at the top right panel. Bottom right: Global H I profile. The black line is a smoothed version of the profile using at 20 km s−1 velocity resolution.

In the text
thumbnail Fig. E.1.

A screenshot of the three-dimensional visualisation of HCG 90. The blue circles indicate the position of the member galaxies. The 2D greyscale image is a DeCaLS R-band optical image of the group. Interactive visualisations of the group can be accessed online and at https://amiga.iaa.csic.es/x3d-menu/.

In the text
thumbnail Fig. E.2.

Left panel: Velocity versus right ascension of HCG 97. Middle panel: Median noise values of each RA-DEC slice of the non-primary beam corrected 60″ data cube of HCG 97 as a function of velocity. The horizontal dashed line indicates the median of all the noise values from each slice. Right panel: Integrated spectra; the blue solid line indicates the MeerKAT integrated spectrum of HCG 97; the red solid line indicates VLA integrated spectrum of the group derived by (Jones et al. 2023). The vertical dotted lines indicate the velocities of the galaxies in the core of the group. The spectra have been extracted from areas containing only genuine H I emission.

In the text
thumbnail Fig. E.3.

H I moment maps of HCG 97. The left panels show all sources detected by SoFiA. The right panels show sources within the rectangular box shown on the left to better show the central part of the group. The top panels represent the column density map with contour levels of (3.5  ×  1018, 7.0  ×  1018, 1.4  ×  1019, 2.8  ×  1019, 5.6  ×  1019, 1.1  ×  1020, 2.2  ×  1020) cm−2. The contours are overlaid on DECaLS R-band optical images. The bottom panels show the moment one map. Each individual source has its own colour scaling and contour levels to highlight any rotational component.

In the text
thumbnail Fig. E.4.

Top left: H I column density map overlaid on DECaLS optical image of HCG 97b. The contour levels are (2.10, 4.20, 8.41) × 1018 cm−2. Top centre: Moment-1 map. The arrows indicate the slices from which the position-velocity diagrams shown at the bottom panels were derived. Top right: Signal-to-noise ratio map. Bottom left: Global H I profile. Bottom centre: Major axis position-velocity diagram. Bottom right: Minor axis position-velocity diagram.

In the text
thumbnail Fig. E.5.

A screenshot of the three-dimensional visualisation of HCG 97. The blue circles indicate the position of the member galaxies. The 2D greyscale image is a DeCaLS R-band optical image of the group. Interactive visualisations of the group can be accessed online and at https://amiga.iaa.csic.es/x3d-menu/.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.