Open Access
Issue
A&A
Volume 690, October 2024
Article Number A69
Number of page(s) 19
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/202449441
Published online 30 September 2024

© The Authors 2024

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

In the hierarchical model of galaxy evolution, dwarf galaxies (M ≲ 109 M) are the most abundant in the Universe and form the building blocks of more massive galaxies, such as our own Milky Way (MW). Dwarf galaxies typically have dynamical masses more than four times larger than their baryonic mass (e.g. Das et al. 2020), which makes them ideal laboratories in which to test cosmological evolution and model dark matter (DM), which accounts for up to ∼27% of the total energy density of the Universe (e.g. Planck Collaboration VI 2020). Nevertheless, little is known about the formation and evolution of these low-mass galactic systems. Dwarf galaxies can be strongly and rapidly affected by their environments. Tidal interactions, for example, are known to generate a specific class of dwarf galaxies (i.e. tidal dwarfs) that in their formation are depleted of DM (e.g. Duc & Mirabel 1998; Read et al. 2017; Venhola et al. 2017; Buck et al. 2019; Gray et al. 2023). Conversely, hydro-dynamical interactions with the intergalactic medium (IGM) can strip and deplete the gaseous disks of these galaxies (e.g. McConnachie et al. 2007; Westmeier et al. 2011; Józsa et al. 2021; Yang et al. 2022). Dwarfs, then, help us build a complete picture of how the environment affects the availability of their gas reservoir, ultimately regulating their star formation history. They are fundamental to understanding the dynamical processes driving the assembly of the stellar and gaseous disks of galaxies.

The new generation of wide-field optical imaging surveys, such as the Dark Energy Survey (DES, e.g. Dey et al. 2019; Martínez-Delgado et al. 2023) or the VST Early-type Galaxy Survey (VEGAS, e.g. Capaccioli et al. 2015; Iodice et al. 2020) have enabled the discovery and study of the stellar component of hundreds of dwarf galaxies down to very low surface brightnesses, that is ∼27 − 29 mag/arcsec2 in the g-band (Forbes et al. 2020; Iodice et al. 2021). However, the analysis of the gaseous component, which is the best tracer for the kinematics and consequently the DM halo and environment, has been much more difficult. A galaxy at the low-mass end of the population (M ≲ 107 M), with a gas fraction of 50%, would require a sensitivity to neutral atomic gas masses of Mgas ≲ 5 × 106 M. Hence, a search for the H I in these sources has been limited so far only to the nearby Universe (≲10 Mpc) and by the sensitivities of H I surveys such as, for example, the Westerbork observations of Neutral Hydrogen in Irregular and SPiral galaxies (WHISP, Swaters et al. 2002), the Very Large Array survey of Advanced Camera for Surveys Nearby Galaxy Survey Treasury galaxies (VLA-ANGST, Ott et al. 2012), the ‘(almost) dark galaxies in the ALFALFA survey’ (Leisman et al. 2017), The Local Volume H I Survey (LVHIS; Koribalski et al. 2018), Apertif (Adams et al. 2022), and the Wide-field ASKAP L-band Legacy All-sky Blind Survey (WALLABY; Westmeier et al. 2022). Only in a handful of dwarf galaxies have resolved studies of the bulk of their H I kinematics been possible, such as, for example, in the sample of The H I Nearby Galaxy Survey of the Local Irregulars That Trace Luminosity Extremes (LITTLE THINGS; Hunter et al. 2012; Oh et al. 2015; Iorio et al. 2017), the dwarf galaxies of LVHIS and VLA-ANGST (Mancera Piña et al. 2021), or in, for example, Leo T (Ryan-Weber et al. 2008; Adams & Oosterloo 2018) and Leo-P (Bernstein-Cooper et al. 2014).

Recently, the combination of deep 21-cm H I observations (MH I ≲ 1019 cm−2) with deep optical photometry (μr ≳ 27 mag/arcsec2) allowed us to discover several new gas-rich and low-mass low-surface-brightness objects beyond the Local Group (LG). The MeerKAT Fornax Survey (Serra et al. 2023), for example, has investigated the H I counterparts of all dwarf galaxies in the Fornax cluster (D ∼ 20 Mpc) down to an H I mass limit of MH I ∼ 5 × 105 M (Kleiner et al. 2023), which were discovered and catalogued by the optical Fornax Deep Survey down to a surface brightness limit of ∼26 mag/arcsec2 in the r-band (Venhola et al. 2017, 2022).

In clusters, the environment strongly affects the H I content of dwarf galaxies. In Virgo, for example, the H I disks of dwarfs are stripped by ram pressure (RP) by the IGM even in the low-density outskirts of the cluster (Boselli et al. 2008a,b). In the less dense medium of Fornax, tidal and hydrodynamical interactions efficiently deplete the H I in dwarf galaxies in between a few tens of millions of years and a few hundred million years (Kleiner et al. 2023). The study of low-mass objects in groups and in the field is so far more incomplete, the main limitation being the difficulty in reaching very deep optical and H I surface brightness limits and in covering, at the same time, a wide-enough field of view to be able to discover these sources, which are often found hundreds of kilo-parsecs away from their host galaxies. The MeerKAT H I Observations of Nearby Galactic Objects: Observing Southern Emitters survey (MHONGOOSE, de Blok et al. 2016, 2020, 2024) is now providing us with the unique opportunity to discover low H I mass sources (MH I ∼ 106 M) within 1 × 1 deg2 of 30 nearby star forming galaxies (median distance D = 10.3 Mpc). Because of the sensitivity to diffuse H I emission over the full field of view guaranteed by MeerKAT, MHONGOOSE investigates the gaseous environment of these galaxies out to almost their virial radius. The coincidence of the H I gas detected in the environment with the optical counterpart detected by sensitive optical observations (μr ≳ 27 mag/arcsec2) allows us to construct a complete census of the low-mass gas-rich galaxy population around the target galaxies. The MHONGOOSE galaxies are mostly isolated or live in nearby sparse groups and span an H I mass range between 106 and 1011 M. The first data release of the shallow observations of the first 10% of the MHONGOOSE data already mapped tens of companion satellite galaxies (de Blok et al. 2024).

In this paper, we present the discovery of the lowest H I mass system (MH I = 1.68 × 106 M) with an optical counterpart detected in the Dorado group (D = 17.69 Mpc, e.g. Rampazzo et al. 2020) by the MHONGOOSE deep observations of the star-forming spiral galaxy NGC 1566. This group was observed in the neutral hydrogen by previous surveys (e.g. Kilborn et al. 2005; Elagali et al. 2019), but never with the depth and resolution provided by MHONGOOSE, which led to the discovery of this and other H I counterparts of low-mass galaxies in the group. The newly discovered H I source (MKT 042326-551621) is the cold gas counterpart of the low surface brightness galaxy (LSBG) catalogued, for example, in Tanoglidis et al. (2021) with source number 19696. The authors associated this source with the density peak of their distribution of LSBGs in the surroundings of the Fornax cluster, at 19 Mpc. For the very first time, because of the HI systemic velocity of this source, we can determine that LSB-D lies in the outskirts of the Dorado group. Given its location in Dorado, and its low-surface-brightness nature, we name this galaxy LSBG-Dorado-A, and for brevity, we hereafter refer to it as LSB-D.

Deep optical observations from VEGAS of the Dorado group allow us to study the optical counterpart of the H I detection. This galaxy is the least-massive low surface brightness galaxy (M ∼ 2.45 × 106 M) with highest gas fraction (fgas = Mgas/(Mgas + M) = 0.48) detected beyond the Local Group. Thanks to the high sensitivity and resolution of both MeerKAT and VST observations we are able to resolve and study in detail the H I gas of such a small object in Dorado and relate its morphology and kinematics to its star formation history and evolution within the group’s environment.

This paper is structured as follows: in Sect. 2, we present the MeerKAT H I and optical VST observations and show the morphology and kinematics of the H I in LSB-D and its optical and star formation properties. In Sect. 3 we discuss the location of this dwarf galaxy within the Dorado group, its baryon content and we present possible scenarios on the evolution of its H I disk and SF. Section 4 summarizes the results and examines the possible studies of low-mass sources that are being opened by these novel deep H I and optical observations.

2. Observations and results

2.1. MHONGOOSE neutral hydrogen observations

The 21-cm MeerKAT (Camilo et al. 2018; Mauch et al. 2020) observations presented in this paper were taken as part of the MHONGOOSE survey targeting the massive lopsided spiral galaxy NGC 1566 (NH I = 1.1 × 1010 M, DNGC1566 = 17.6 Mpc Elagali et al. 2019; Anand et al. 2021). NGC 1566 was observed for 55 hours in ten separate 5.5 hour tracks reaching H I column density sensitivities of NH I(3σ) ∼ 6 × 1018 cm−2 at 30″ resolution, over 16 km s−1.

A complete description of the data reduction of the MHONGOOSE observations is given in the survey paper (de Blok et al. 2024). Here, we summarise the main steps taken to produce the H I datacubes1 that led to the detection of the gas rich low surface brightness galaxy LSB-D. We reduce the data using CARACal (Józsa et al. 2020), a containerised pipeline written in Python3. CARACal is set up in a modular fashion allowing us to use tasks from different radio astronomical packages for each step of the data reduction. The data-reduction steps are performed on a single node of the MeerGas cluster (128 cores, 1 Tb RAM). All ten tracks are reduced following the same strategy and using the same set of parameters, except for a few minor exceptions.

The data reduction procedure is divided in three parts. The first part is performed on the single 5.5 hours tracks and consists of cross-calibration, RFI-flagging, self-calibration and continuum subtraction. The second part consists of combining all ten observations together in the uv-space, followed by the joint deconvolution and imaging to produce H I datacubes at multiple resolutions with a 1.4 km s−1 channel width, without further continuum subtraction. In the third part, we create the moment maps (i.e. H I flux density, velocity field and H I line-width map) of the imaged field of view (1.5 × 1.5 deg2) using the automated source finder SoFiA-2 (Westmeier et al. 2021).

We produce the final datacubes with six weighting and tapering combinations, that span resolutions from 90″ to 7″, with 3σ sensitivity limits ranging from 5 × 1017 cm−2 to 5 × 1020 cm−2 over 16 km s−1, respectively (further details are given in Appendix A). We use SoFiA-2 to search for H I detections with a smooth & clip algorithm, with Gaussian kernels 1, 2 and 3 times the spatial resolution of the cubes. The detection threshold is set to 3.5 times the noise (evaluated in 3D over a spectral window of 7 channels). All H I sources detected in the field of view have a detection reliability of 85%, based on the distribution of positive and negative noise peaks in the datacubes (see Serra et al. 2015; Westmeier et al. 2021, for details on the reliability estimate). We verify that SoFiA-2 identified all genuine H I emission by also visually inspecting the datacubes in virtual reality (Jarrett et al. 2021, iDaVIE).

The H I flux density distribution of NGC 1566 and of the galaxies detected by MeerKAT at resolutions between 90″ to 7″ is shown in Fig. 1 over a combination of the g, r, z optical observations of the Dorado group from the Dark Energy Camera Legacy Survey (dr10-DECaLS). The doublet of early-type galaxies NGC 1549 and NGC 1553 in the centre of the group is shown in the south-west, beyond the imaged MeerKAT field of view (dashed circle). New H I detections are marked in green. The first contour of each colour (from yellow-green to red) shows the S/N = 3 H I column density limit of datacubes at resolutions between 90″ and 12″. MeerKAT has a primary beam diameter of 1 degree (at this diameter the sensitivity of the array drops to 50%). This enables us to study not only NGC 1566 and its interactions with its nearby satellites, but also to investigate in detail its environment out to ∼360 projected-kpc, which corresponds approximately to 75% of its radius at 500 times the critical density of the Universe at the current epoch (r500 = 0.47 Mpc, Kilborn et al. 2005).

thumbnail Fig. 1.

Primary beam corrected flux-density H I emission detected by MeerKAT within a 1.5 × 1.5 deg2 field centred on NGC 1566 (the imaged f.o.v. is marked by the dashed grey circle), overlaid on the DECaLS optical image in the g, r, z filters. The H I emission is a composite of flux-density maps derived from the multi-resolution datacubes. The resolutions shown by coloured contours are 12″ × 10″ (red), 25″×18″ (green), 32″×23″ (blue), 65″×64″ (purple), 94″×92″ (yellow-green). For resolutions between ∼90″ and 25″ two contours are shown, i.e., NH I = 3σ × 2n with n = 0, 2. The 12″×10″ resolution has further increasing contours (n = 0, 2, 4...). The PSF of the multi-resolution maps are shown in the bottom-right corner with the same colour coding. Details on the H I properties of the sources in the field are given in Table A.2.

The H I MHONGOOSE observations of NGC 1566 are the deepest and highest resolution 21-cm observations available so far of this galaxy (a summary of the sensitivities is shown in Table A.1). They improve our characterization of the low-column densities features of its H I disk of at least an order of magnitude (Elagali et al. 2019) and allow us to detect H I associated with several satellites, as well as diffuse H I tails tracing past interactions (as in NGC 1581, in the east of the field of view, Fig. 1). A full description of the cold gas distribution and kinematics of NGC 1566 and the other detections will be given in a following paper, a summary of their H I properties is shown in Table A.2.

Among the new sources, to the south-east of NGC 1566, at 186 kpc in projection at coordinates RA = 04h23m26s, Dec = −55d16m21s, lies the lowest H I detection associated with a galaxy made by MHONGOOSE so far, the aforementioned LSB-D (circled in red in Fig. 1). Deep optical observations allow us to identify that the H I gas is co-located with a low surface brightness galaxy. The H I systemic velocity of LSB-D (vsys = 1214 km s−1) is consistent with the average recession velocity of the Dorado group (vhel = 1269 km s−1) and the source seems to be located at its outskirts. Dorado is a loose group with centre coincident with the S0 galaxy NGC 1553 (visible in the south-west in Fig. 1), which forms a pair with the brightest group galaxy (BGG) NGC 1549 (Firth et al. 2006; Kourkchi & Tully 2017). LSB-D is located to the north-east of Dorado’s centre (68′ away, see Fig. 1), whereas NGC 1566 and its companions are at a distance of 35′. A summary of the main properties of Dorado is shown in Table 1.

Table 1.

Properties of the Dorado group.

The H I disk of LSB-D is detected at resolutions between 10 and 40 arcseconds with peak signal-to-noise Speak/N ≳ 5, while the gas becomes too beam-diluted to be detected at the lowest resolutions and below the detection limit at the highest resolution of our datacubes (7″). This highlights the compact but diffuse nature of the H I disk, see Sect. 3.2 for further details. In what follows, we analyse the results only from the 25″ × 18″ resolution datacube, with 1.4 km s−1 channel-width generated with no tapering and a robustness parameter of 1.0, (Briggs et al. 1999, see Table A.1). At this resolution, we fully recover the total extent of the H I, spatially resolving it with 3 beam elements, while at higher resolutions we resolve out part of the disk. The noise in the primary beam corrected cube close to LSB-D is 0.265 mJy beam−1. In Appendix B we show the channel maps of the H I emission of LSB-D.

2.2. H I properties of LSB-D

The integrated H I spectrum of LSB-D is shown in Fig. 2. The H I line spans over 17 channels (1.4 km s−1 wide) with a full-width at zero intensity (FWZI) of wFWZI = 24 km s−1. The profile is double-peaked suggesting the gas may be rotating within the galaxy and an extended tail is present at blueshifted velocities. The total flux density of the H I line is S1.4 GHzΔv = 22.9 mJy km s−1. The H I properties of LSB-D are summarized in Table 2.

thumbnail Fig. 2.

Integrated spectrum of LSB-D, with its noise per channel (grey shaded region). The vertical dashed line shows the systemic velocity of the source ( v sys radio = 1214 $ v^{\mathrm{radio}}_{\mathrm{sys}} = 1214 \, $km s−1), while the grey vertical lines show its full velocity range (wFWZI = 24 km s−1).

Table 2.

H I properties of LSB-D.

In a group, the projected line of sight velocity we measure from the H I line is the composition of the systemic velocity of the source and its relative projected motion within the group. The exact trajectory of LSB-D is unknown, but since its systemic velocity (vsys,  LSB-D = 1214 km s−1) is similar to that of the group centre (vsys,  NGC1553 = 1201 km s−1), and of the close-by BGG NGC 1549 (vsys,  NGC1549 = 1202 km s−1), it is reasonable to approximate its distance by that of the group centre (D = 17.69 Mpc)2, which has been estimated by averaging the radial velocities of their group members (Kourkchi & Tully 2017; Rampazzo et al. 2020). At this distance LSB-D lies in projection 351 kpc away from NGC 1553 and 181 kpc and 116 kpc from NGC 1566 and NGC 1581, respectively (vsys,  NGC1566 = 1496 km s−1, vsys,  NGC1581 = 1609 km s−1). Further information on the location of LSB-D within the group is given in Sect. 3.1.

Following standard equations (e.g. Meyer et al. 2017), the total flux of LSB-D corresponds to an H I mass of MH I = 1.68 × 106 M. So far, this is the lowest H I mass that has been associated with a low surface brightness galaxy beyond the Local Group. For comparison, the typical dwarf population in the LG has associated H I masses two orders of magnitudes higher than in LSB-D, with only a few sources with MH I ∼ 1.68 × 106 M (e.g., Bouchard et al. 2006; Grcevich & Putman 2009; Hunter et al. 2012; Spekkens et al. 2014; Iorio et al. 2017).

The left panel of Fig. 3 shows the H I and optical body of LSB-D in detail. The first contour shows the S/N = 3 flux-limit (before primary beam correction) at 25″ × 18″, which corresponds to a column density of NH I = 4.6 × 1018 cm−2. Further details on the spatial and velocity extent of LSB-D are visible in the channel maps (Fig. B.2). The H I gas is distributed in an irregular morphology that at first look resembles a disk with projected major axis a = 36″ (DH I = 2.8 ± 0.8 kpc), and minor axis b = 22″ = 1.8 ± 0.8 kpc. The major axis is misaligned with respect to the optical body (PAH I = 141° and PA = 190°). Along the minor axis, in the north-east, the H I elongates in a low-column density cloud (NH I = 3 − 9 × 1018 cm−2). This cloud contributes to the blue-shifted tail in the H I spectral profile (Fig. 2), and is visible in the channel maps at velocities 1209 km s−1 and 1212 km s−1. We highlight that also a low-mass H I cloud without associated stars may be present (at the detection limit of our observations) 8′ further north of LSB-D (41-projected kpc) at coordinates RA = 04h22m31s, Dec = −55d12m02s, The integrated spectrum is shown in Fig. B.1 and has full-width at zero intensity of 18 km s−1. The total estimated H I mass is 3.65 × 105 M, which is slightly higher than the 3σ detection limit over 16 km s−1 of the observations.

thumbnail Fig. 3.

Left Panel: DeCaLS optical image in the g, r, z filters overlaid with the H I contours from the 25″ × 18″ datacube used for the analysis in this paper. Levels increase as 4.6 × 1018 × 2n cm−2 (n = 0, 1, 2, 3), where the first contour marks the mean S/N = 3 detection limit. The PSF is shown in the bottom-right corner. Note the H I overlay in color with also the 32 × 23 ″ map. Right Panel: velocity field of the H I gas in LSB-D. The systemic velocity (vsys = 1214 km s−1) is marked by the thick black iso-velocity contour, the other contour levels are shown in the bottom-left corner. The dashed and dashed-dotted lines mark the directions of the major (141°) and minor axes (51°) of the H I disk, respectively.

The right panel of Fig. 3 shows the velocity field of LSB-D. The H I gas is characterised by a gradient of velocities along the major axis. One interpretation is that regular rotating motions are ongoing in the gas. Nevertheless, some asymmetries are present in the velocity field with a deep V-shape penetrating the minor axis in the south-west. Independent of this, the morphology of the H I shown in the left panel of Fig. 3 is suggestive of a disk-like distribution. Under this assumption we assume we can then infer an inclination from the axis ratio (b/a), following, for example Sakai et al. (2000):

i = cos 1 ( b / a ) 2 q 0 2 1 q 0 2 $$ \begin{aligned} i = \cos ^{-1}\sqrt{\frac{(b/a)^2-q_0^2}{1-q_0^2}} \end{aligned} $$(1)

where q0 = 0.13 − 0.27 is the intrinsic minor-to-major axis ratio of edge-on spirals. Given the irregularity of the H I morphology, we determine the range of possible inclinations assuming the extreme values of the q0 range and propagating the errors on the axis, which we estimate to be half the size of the beam minor axis, 9″ (0.8 kpc). We obtain an inclination of i = 50° ±15°.

To further investigate the kinematics of the H I in LSB-D, we draw the position-velocity (PV) diagram across the major (PA = 141°) and minor axis (PA = 51°) of the H I disk (marked in the velocity field of Fig. 3 with a dashed and dashed-dotted line, respectively). The PV-diagram across the major axis (Fig. 4, top panel) shows a steep gradient of velocities which disappears in the PV-diagram along the minor axis (bottom panel of the Figure), this may indicate that the H I is supported by rotation. Along the velocity gradient, the H I is detected only over three consecutive channels (of 1.4 km s−1 each) indicating that the velocity dispersion is low ( σ = w 50 / ( 2 ln 2 ) 3 ± 1.4 $ \sigma=w_{50}/(2\sqrt{\ln{2}}) \lesssim 3\pm1.4 $ km s−1) as typically found for the cold H I in rotating disks (T < 1000 K; Young et al. 2003; Warren et al. 2012). Under the assumption that rotation is (mostly) supporting the disk, a proxy for its rotational velocity can be inferred from the inclination-corrected velocity width (Sakai et al. 2000; Begum et al. 2008): vrot = 1/2 w20sin(i)−1(1 + z)−1 km s−1. From the smooth gradient of velocities across the major axis we measure w20 = 14 km s−1 and vrot = 9.0 ± 3.5 km s−1. Because of the difficulty in estimating the rotational velocity we estimate an error of ∼40%.

thumbnail Fig. 4.

Top Panel: position-velocity diagram taken along the major axis of the H I disk (PA = 141°) with a slit of 20″. Contour levels are at 5 × 1017 cm−2 and 1 × 1018 cm−2. Bottom Panel: position-velocity diagram taken along the minor axis of the H I disk (PA = 51°) with an analogous slit, equal to the resolution of the observations. Contours are as in the top panel.

Figure 4 shows that in the centre, a second peak is present at v = 1218 km s−1 beyond the range of rotational velocities traced by the smooth gradient. Without considering the beam-smearing due to the turnover velocities in the central resolution beam, the second peak causes the velocity dispersion to increase to σ = 5.5 km s−1. This central velocity dispersion is typical of low-mass systems (e.g. Iorio et al. 2017; Mancera Piña et al. 2021) and implies vrot/σv ∼ 2. This suggests that the central region of the H I disk has more turbulent kinematics, as also indicated by the shift of its centre with respect to the optical body (Fig. 3, left panel).

2.3. Deep optical data and stellar properties of LSB-D

Combining the deep H I observations with optical photometric observations we can state that LSB-D is without any reasonable doubt a gas rich low surface brightness galaxy. To properly characterize its stellar properties we rely on the deepest available photometric observations of the Dorado group from VEGAS (Spavone et al. 2017; Iodice et al. 2021). These observations reach sensitivities of the limiting surface brightness of μg = 28.9 mag/arcsec2, μr = 27.3 mag/arcsec2, in the g and r filters respectively, estimated for a point source at 5σ over a circular area with FWHM ∼ 1″. These limits are approximately two magnitudes deeper than the DECaLS surface brightness limit (Martínez-Delgado et al. 2023). The VEGAS observations were acquired using the VST/OmegaCAM instrument (run ID: 0102.A-0669(A)), which has a total field of view (FoV) of 1° ×1° and a spatial resolution of 0.214 arcsec/pixel (Kuijken 2011). The images were reduced using the dedicated AstroWISE pipeline (for details see McFarland et al. 2013; Venhola et al. 2017) specifically developed for the reduction of VST/OmegaCAM observations. The g and r bands were analysed by following the same procedure described in, for example, Ragusa et al. (2021, 2022, 2023), to which we refer for more details. In particular, the magnitudes are in the AB system, and are corrected for Galactic extinction using the coefficients from Schlafly & Finkbeiner (2011). The VEGAS g, r image of LSB-D is shown in Fig. 5.

thumbnail Fig. 5.

VEGAS optical image in the g, r filters overlaid with the H I column density map. Orange contour levels are as the left panel of Fig. 3. In cyan and magenta contours we show the S/N = 3 and 5 contours from the NUV and FUV GALEX observations, respectively.

We infer the optical properties of LSB-D with an isophote fitting procedure in the g and r bands. The surface brightness profiles are shown in the top panel of Fig. 6. The VEGAS observations confirm that the peak of the optical distribution of LSB-D is shifted with respect to the H I peak by ∼0.4 kpc. We determine the stellar mass (M) from g and r magnitudes (Mr = −13.47, Mg = −13.04) following the procedure illustrated in Into & Portinari (2013). Briefly, based on a mass-to-light (M/L) - color relationship, it is possible to derive an estimate of the average galaxy’s M/L by using its g − r color. This, in turn, leads to an estimate, passing through the luminosity, of a stellar mass M (r) = 2.2 × 106 M and M (g) = 2.7 × 106 M, respectively. The stellar properties, along with the average g − r colour, are listed in Table 3.

thumbnail Fig. 6.

Top panel: LSB-D’s g and r radial surface brightness profile (cyan and red, respectively) inferred from the VEGAS observations. Bottom panel: radial g − r colour profile measured out to 8″ from the centre (∼670 pc).

Table 3.

Stellar properties of LSB-D.

LSB-D has blue colours, similar to other nearby dwarf galaxies (e.g. Honey et al. 2018), and small effective radius, reff = 0.5 kpc. Its radial colour profile (bottom panel of Fig. 6) shows a moderate reddening from the centre to the outskirts. Bluer colour in the centre suggest the presence of a younger stellar population or on-going star formation.

To further investigate this possibility, we examine if ultraviolet (UV) emission is associated with LSB-D in the Galaxy Evolution Explorer Nearby Galaxy Survey (Gil de Paz et al. 2007; Bianchi & Garcia 2014). Both in the far-UV (λ ≈ 154 nm, FUV) and near-UV (λ ≈ 231 nm, NUV) observations emission is detected in the H I disk of LSB-D, see the magenta and cyan contours in Fig. 5. Given its low-signal-to-noise, the UV observations on their own would unlikely have been considered a detection. Since the peak is present both in the FUV and in the NUV and is coincident with the centre and the stellar body of the galaxy, we consider the detection significant.

The higher velocity dispersion of the central H I gas may be due to the on-going or recent star formation suggested by the UV emission. The star formation rate (SFR) can be reliably estimated from the FUV emission (e.g. Leroy et al. 2019). We measure the UV-flux by centering a five-by-five pixel box on the peak, obtaining a total flux FFUV = 0.031 counts/sec. Correcting for extinction assuming a coefficient of AFUV = 0.13 from Wall et al. (2019) and given that the GALEX observation has an average background emission (as measured in a few areas nearby the source) of Fbackgroud = 0.00035 counts/sec, the corrected flux is FFUV = 0.025 counts/sec. We determine the SFR from the correlation with the FUV, as inferred from the cross-matched sample of GALEX-SDSS galaxies (Salim et al. 2007; Leroy et al. 2008):

SFR FUV unobscured = 0.68 × 10 28 L ν ( F U V ) M yr 1 $$ \begin{aligned} \mathrm{SFR}^\mathrm{unobscured}_{\rm FUV} = 0.68 \times 10^{-28} L_\nu (\mathrm {FUV}) \ \mathrm{M}_\odot \,\,\mathrm {yr}^{-1} \end{aligned} $$(2)

where L ν ( FUV ) $ L_\nu(\rm FUV) $ is in  erg s−1 Hz−1. For LSB-D, L ν ( FUV ) = 9.6 × 10 24 $ L_\nu(\rm FUV) = 9.6\times10^{24} $ erg s−1 Hz−1, which gives log(SFR)LSB − D = −3.2 [M yr−1]. LSB-D’s SFR value is a factor of only a few below the SFR of the lowest mass MHONGOOSE target galaxies (see Table 1 of de Blok et al. 2024) and it is similar to the SFR in the low-mass end (M ≲ 107 M) of dwarf star-forming galaxies in the Local Universe (McQuinn et al. 2015a,b; Marasco et al. 2023).

3. Discussion

In this section we discuss the possible evolution of LSB-D and its H I disk as it can be inferred from its location within the Dorado group (Sect. 3.1), by its gaseous and stellar content (Sect. 3.2) and by its baryon-to-dark-matter content (Sect. 3.3). In Sect. 3.4, we present possible scenarios on the formation and evolution of its H I disk and on the enhancement of SF.

3.1. LSB-D in the Dorado group

LSB-D lies 390 projected-kpc to the north-east of the centre of Dorado (formed by the galaxy pair NGC 1549/NGC 1553), while the closest (in projection) massive galaxies NGC 1566 and NGC 1581 are 186 kpc and 190 kpc to the north-west and north-east, respectively. Since LSB-D is at the edge of their projected virial radii it has likely not yet tidally interacted with them. Figure 1 shows that LSB-D is too far away from other more massive group members to be considered their satellite. The closest galaxy to LSB-D is MKT 042200-545440 (also associated for the first time to Dorado thanks to our MeerKAT observations) which is a late-type galaxy with similar H I mass (see Table A.2).

The phase-space diagram of Dorado can provide further insights related to the location of LSB-D and on its infall into the group. This diagram (see, for example, Jaffé et al. 2015) takes into account the distance and systemic velocity of the group members as well as the virial mass and velocity dispersion of Dorado (Mvir = 3.5 × 1010 M, Rvir = 528 kpc, σDorado = 282 km s−1) to understand where the sources are located in phase space with respect to the BGG, if they are already within its escape velocity, if they recently assembled into the group or if they have long been virialised. Figure 7 shows the phase-space diagram of all Dorado members as catalogued in Rampazzo et al. (2020, 2021, 2022) to which we add the sources newly associated with the group thanks to the MeerKAT observations, i.e. LSB-D, PGC 128826, PGC 075137, MKT 042200-545440, MKT 042230-550757, MKT 042326-551621. The Figure shows the caustic curves defined by the escape velocity of the group (which is derived from its virial mass and radius) and the virialized area of the group (dashed-dotted lines). As illustrated in Rhee et al. (2017), within this latter region it is more likely to find old members of the group between which interactions have already occurred, compared to recent infallers which are located between the virialized region and the caustic lines.

thumbnail Fig. 7.

Phase-space diagram of the Dorado group. LSB-D is shown in red, while previously known H I rich galaxies are shown in blue and new H I detections in green. The group’s centre (NGC 1553) is at the zero coordinates, the dashed-dotted lines mark the virialised region of the Dorado group. The caustic curves (dashed lines) are defined by the escape velocity of the group.

In the phase-space diagram, LSB-D is within the escape velocity of the group in the region of recent infallers. The recent infallers are expected to be more gas-rich than similar systems in the virialized region, since the environmental effects of the group (i.e. tidal interactions and/or ram pressure stripping) have not yet been strong enough or did not have enough time to deplete the sources of their gas.

To further investigate if LSB-D is a recent infaller of Dorado or if it had multiple orbits within the group, already, we estimate the probability distributions for orbital parameters and deprojected coordinates in group and clusters based on their phase space coordinates. In particular we follow the methodology illustrated in Oman et al. (2021, Sect. 3.2) (see also Oman et al. 2013). In brief, from an N-body simulation catalogue we select within an overdensity similar to Dorado all orbits of satellites with projected coordinates, host mass and satellite mass compatible with LSB-D. Further details on the selection constraints are given in Appendix C. From the distribution of deprojected coordinates of the selection, we derive the probability for the deprojected coordinates of LSB-D and for its orbit within the group. The median deprojected coordinates of LSB-D (with 16 − 84 percentile range) are r 3 D = 506 129 + 481 $ r_{\mathrm{3D}} = 506^{+481}_{-129} $ kpc, v 3 D = 509 242 + 208 $ v_{\mathrm{3D}} = 509^{+208}_{-242} $ km s−1. These are shown in Fig. C.1 along with the distribution of the deprojected coordinates of all selected orbits. The probability for LSB-D to be an interloper (r3D > 2.5Rvir) is ≲1%. The distribution of the pericentre of the selected orbits is bimodal with modes corresponding to either first or to second and later orbits. The probability for LSB-D to be past apocentre is ∼50%, but if it was within the first orbit the probability for LSB-D to be in the first infall before passing through the first pericentre is ∼27%, while the probability to be on the first orbit but past pericentre is ∼23%. In conclusion, the position in the phase space diagram and the predictions on its deprojected coordinates suggest that LSB-D is a recent infaller of Dorado. The simulations suggest that in an overdensity similar to Dorado, galaxies with properties similar to LSB-D with those projected coordinates are typically found on their first or at the beginning of the second orbit within the group.

3.2. Gas and stellar content of LSB-D and its star formation rate

LSB-D belongs to the low-mass end of the H I-detected dwarf population. Here we compare its H I properties with its stellar content and with those dwarf galaxies in the Local Universe for which we have reliable estimates. In particular, we select sources from LITTLE THINGS (Oh et al. 2015; Iorio et al. 2017), SPARC (Spitzer Photometry and Accurate Rotation Curves, Lelli et al. 2016) and The H I Nearby Galaxy Survey (THINGS, Walter et al. 2008) samples, local low surface brightness galaxies (from McGaugh 2011; McGaugh et al. 2017) and ultra-diffuse galaxies (UDGs, from Mancera Piña et al. 2019, 2022b), for which we have precise estimates of the gas and stellar masses, as well as rotational velocities. The top panel of Fig. 8 shows the stellar mass against the gas mass for the mentioned samples and for LSB-D, which nicely fits on the one-to-one relation (shown by the dashed line). This Figure also highlights the low star and gas mass of LSB-D of which only a few analogues in the Local Universe have been studied in detail.

thumbnail Fig. 8.

Top panel: gas mass versus stellar mass for a sample of galaxies in the Local Universe. LITTLE-THINGS, SPARC and THINGS galaxies are marked by triangles (in orange, grey, and blue, respectively). Local low surface brightness galaxies (LSBGs) from McGaugh (2011, 2017, 2021) are in green, and UDGs in cyan. LSB-D is marked by the red circle. The dashed line marks the 1 : 1 linear relation. Bottom panel: the baryonic Tully-Fisher relation of Local Group dwarf and spiral galaxies (dashed line; McGaugh et al. 2021). Symbols are as in the top panel. For the Local LSBGs sample, pressure supported dwarfs (plotted against their velocity dispersion, σv) are marked by light green squares, while green triangles mark rotationally supported galaxies, plotted against their rotational velocity, vrot.

The baryonic Tully-Fisher (BTFR) relation links the baryonic mass of a galaxy with its rotational velocity (e.g. McGaugh et al. 2000; Ponomareva et al. 2018; Lelli et al. 2019). Whereas the relation is well established in the high mass regime (e.g. McGaugh 2012; McGaugh et al. 2021; Lelli 2022), its nature in the low-mass end is still a matter of heated debate as the scatter appears to increase towards lower masses (Geha et al. 2006; Begum et al. 2008; Mancera Piña et al. 2019). We investigate where LSB-D is located on the BTFR to understand if dynamically its overall properties are typical of the ones of the low-mass low surface brightness galaxies.

In the bottom panel of Fig. 8 we show the BTFRs from McGaugh et al. (2021):

v 0 = 0.379 ( M bar M ) 1 / 4 km s 1 $$ \begin{aligned} v_0= 0.379\,\Bigg (\frac{M_{\rm bar}}{\mathrm{M}_\odot }\Bigg )^{1/4} \,\mathrm{km \,s}^{-1} \end{aligned} $$(3)

The double dashed line in the Figure takes into account that there is a factor of two offset between the velocity estimated from rotation or pressure supported systems (v0 and σ, details are given in McGaugh et al. 2021). In Fig. 8 we added the samples from LITTLE-THINGS, THINGS, and SPARC keeping their selection parameters and estimates of the rotational velocity and baryonic masses (there is a negligible discrepancy with the BTFR estimated with different methods for the v0, Lelli et al. 2019). To this we also added the UDGs which are off the BTFR (Mancera Piña et al. 2019). LSB-D (shown in red) seems to be slightly offset with respect to the relation, even though because of the large uncertainties (in the inclination and rotational velocity of the H I disk) it can be consistent with the relation. We estimated LSB-D’s baryonic mass and rotational velocity as follows. Since we do not spatially resolve the source with enough resolution elements to properly characterize its rotation curve (i.e. ∼5 beams across the minor axis, see Bosma 1978; Sancisi et al. 2008), we estimate a rough dynamical mass enclosed in the LSB-D’s H I disk as M dyn , rot ( r H I ) r H I v rot 2 G 2.6 ± 1.0 × 10 7 $ M_{\mathrm{dyn, rot}}(\leq r_{{\mathrm{H}\,{\tt I}}}) \sim \frac{r_{{\mathrm{H}\,{\tt I}}}\,v_{\mathrm{rot}}^2}{G}\approx 2.6\pm1.0\times10^7 $ M. Accounting for helium (Mgas = 1.33 ⋅ MHI = 2.23 × 106 M), the baryonic mass is Mbar = Mgas + M = 4.68 × 106 M and the total gas fraction is fgas = Mgas/(Mgas + M) = 0.48. LSB-D is dark-matter dominated with its baryonic mass counting 5 − 10% of the dynamical mass. This is compatible, even though slightly higher, with other low surface brightness H I rich galaxies in the Local Universe (Giovanelli et al. 2013; Iorio et al. 2017). Considering that LSB-D may not be fully supported by rotation (since in the centre vrot/σv ∼ 2), if we infer the dynamical mass from the maximum velocity dispersion through the virial theorem (as in other low-mass pressure supported systems, e.g. McGaugh et al. 2021, also shown in the Figure), we obtain a similar result because of the low velocities of the gas along the disk.

The very low mass of LSB-D, the normality of its gaseous and stellar properties and compatibility with the baryonic Tully-Fisher relation shows the universality of the dynamical properties of galaxies which defines their evolution throughout the entire spectrum of masses, down to the smallest low surface brightness sources (albeit with larger scatter in this regime).

Overall, LSB-D shows normal gas, stellar and dynamical properties compared to Local Group galaxies. This suggests that LSB-D has undergone a ‘typical’ evolution throughout cosmic time, and that past recent interactions or the environment of Dorado have not influenced its halo or mass content significantly. Nevertheless, the optical and UV properties of LSB-D suggest, instead, that recently some event must have enhanced its star formation. Though gas rich by the standards of comparable low stellar mass galaxies (i.e. dwarf spheroidals in Fig. 8), LSB-D has slightly lower gas masses than the Local LSBGs of not too much greater stellar mass. However, the top panel of Fig. 9 shows that LSB-D is ≈0.8 dex above the M-SFR relation defined by the Local LSBGs and THINGS galaxies. While it lies among the typically dwarf irregular LITTLE-THINGS galaxies3, some of the offset between LSB-D LITTLE-THINGS and the LSBGs and THINGS sample is due to the different tracers used to measure SF (FUV and Hα, respectively). For low-mass galaxies with log(SFR)≲ − 2.5 [M yr−1] this discrepancy can vary between a factor 2 − 10 (Lee et al. 2009). However, in LITTLE-THINGS galaxies this discrepancy is at most a factor three, which is possibly due to a lack in sensitivity of the Hα measurements (Hunter et al. 2010). Part of the offset between LITTLE-THINGS galaxies and LSBGs may indicate that these dwarf irregulars have higher SFR than the Local LSBGs of comparable mass. Within this scenario, LSB-D appears to be more similar to the SF dwarf irregulars than to the more quiescent LSBGs, which also tend to live in less dense environments.

thumbnail Fig. 9.

Top panel: star formation rate as a function of stellar mass for LSB-D (marked by the red circle, and for a sample of galaxies in the Local Universe. LITTLE-THINGS (Oh et al. 2015; Iorio et al. 2017) and THINGS galaxies are marked by orange and blue triangles. Local low surface brightness galaxies (LSBGs) from McGaugh (2011, 2017) are in green. The dashed line shows a linear relation fitted betwee the LSBGs and THINGS samples. Central panel: specific star formation rate as a function of stellar mass for the same sample as in the top panel. The dashed-dotted line marks the inverse Hubble time, i.e. galaxies on the relation reach their current stellar mass with their average SFR since the Big Bang. Bottom panel: depletion time as a function of stellar mass for LSB-D and the samples of galaxies in the Local Universe. The dashed line marks the linear fit relation for the Local LSBGs sample.

When converting the SFR to specific star formation rate (SFR/M) we see that while Local Universe galaxies lie around or below the inverse Hubble time relation, they would have reached their current stellar mass with their average sSFR since the Big Bang. LSB-D lies above the relation suggesting that it would have doubled its stellar mass in a Hubble time with its current sSFR (Fig. 9, central panel). Hence, for its stellar mass the current LSB-D’s SFR is too high to have always been on-going. Moreover, the depletion time of its gas reservoir (tdep = Mgas/SFR) is much smaller than a Hubble time (Fig. 9, bottom panel) indicating that if LSB-D had formed stars at its current rate it would have already consumed all its gas.

3.3. The baryon-to-dark-matter content of LSB-D

The baryon-to-halo mass relation of galaxies (BHMR) is a powerful tool to study the efficiency with which galaxies turn their available baryons into stars and cold gas. In Fig. 10, we place LSB-D on the BHMR and we compare it with expectations based on abundance-matching (Calette et al. 2021) and against observations of a sample of dwarf and massive nearby disk galaxies spanning five orders of magnitude in baryonic mass from Mancera Piña et al. (2022a). For the latter galaxies the baryon fraction4 is estimated via accurate mass modelling from their rotation curve decomposition (see Mancera Piña et al. 2022a, for details).

thumbnail Fig. 10.

Mbar/M200 ratio normalized to the cosmological baryon fraction (fbar, cosmic). LSB-D is compared against a set of nearby disk galaxies (Mancera Piña et al. 2022a) and the semi-empirical relation from Calette et al. (2021), which uses abundance matching to link M with M200 and empirical correlations for Mstar/Mgas. Galaxies are colour coded according to their gas-to-baryon content (Mgas/Mbar).

For LSB-D we do not have a rotation curve, so we approximate its approximate location in Fig. 10 as follows. We estimate the circular speed of the galaxy as vcirc = vrot. This means that we ignore corrections for line broadening due to turbulence as well as asymmetric drift (which would decrease and increase our vcirc, respectively, e.g. Verheijen 2001; Iorio et al. 2017), which we cannot properly characterize because of the spatial resolution of our observations and considering that both corrections counteract each other. To estimate the uncertainty of vcirc we consider the error propagation in vrot (taking into account our error in the measurement in vrot and in the inclination uncertainty), and we add (in quadrature) an additional term of 3 km s−1 associated with the asymmetric drift correction based on the typical value of such correction in low-mass systems (e.g. Lelli et al. 2014; Iorio et al. 2017; Mancera Piña et al. 2021). We compare our value of vcirc against a set of circular speed profiles of Navarro-Frenk-White (NFW, Navarro et al. 1996) haloes of different halo masses and with concentration parameters following the concentration-mass relation from N-body cosmological simulations (Dutton & Macciò 2014). Given that this relation shows a tail towards low concentrations at low masses (Mancera Piña et al. 2022b; Kong et al. 2022), we also consider NFW profiles with concentration parameters 2σ below the nominal expectation. We find that our estimate of vcirc can be matched with the circular speed of a halo of ∼108M. Using this value, we add LSB-D to Fig. 10, though it is worth noticing that the uncertainties in vcirc allow for haloes in a broader mass range (7 ≲ log(M200/M)≲8.75). LSB-D has baryon-to-halo mass properties similar to dwarf galaxy DDO210 (e.g. LITTLE THINGS, Hunter et al. 2012).

Abundance matching techniques suggest that the BHMR relation steeply and monotonically increases as a function of M up to Mhalo ∼ 1012 M (e.g. Moster et al. 2013), but the behaviour at low-masses is only extrapolated from the relation in the 109–1011 M range. The properties of LSB-D, along with DDO210, suggest the relation may be flatter, and that low surface brightness galaxies may retain their gas within the dark-matter halo while assembling their stellar content. This agrees with the BHMR relation obtained by Calette et al. (2021) who used both abundance matching to link M with M200 and empirical correlations to link M to Mgas, or the relation by Posti et al. (2019) based on H I rotation curve measurements. Alternatively, LSB-D and DDO210 may represent the part of the low-mass population with relatively large baryonic masses for their halo mass, while most of the sources with low baryonic masses are undetected.

LSB-D and DDO210 are two of the few known sources of low-stellar mass for which the stellar-to-mass halo content can be estimated accurately, hence they provide an observational constraint to the low-mass end (≲5 × 107 M) of the stellar-to-mass halo relation. Overall, our analysis on LSB-D adds to the emerging picture of a large scatter in the baryon-content of low-mass galaxies due to feedback from star formation, which modern hydrodynamical simulations should aim to reproduce (McGaugh 2012).

Even though the gravitational potential and the baryon content of LSB-D and DDO210 are similar, their H I properties have some striking differences. First, the minimum H I column density detected in DDO210 is 6 × 1019 cm−2 (which corresponds to a S/N = 3 in the moment-zero map; Iorio et al. 2017) and its H I surface brightness profile peaks in the centre at 6 × 1020 cm−2. LSB-D has H I column densities one order of magnitude lower, peaking at 3.3 × 1019 cm−2, and it would have been at the detection limit of DDO210’s VLA observations (NH I (3σ, 4 km s−1) = 1.3 × 1019 cm−2; Hunter et al. 2012). Second, LSB-D’s H I disk is much more extended (rH I, LSB−D = 1.4 kpc), than the one of DDO210 (rH I, DD0210 = 0.3 kpc). LSB-D and DDO210 likely had a different formation history, which one disk being much more diffuse than the other. LSB-D appears closer in projection to its massive neighbours than DDO210. In the Local Group, DDO210 and other dwarf galaxies with MH I ≳ 106 M are located beyond 300 kpc from the Milky Way or Andromeda. Within 300 kpc dwarf satellites have low amounts of H I (MH I ≲ 106 M) below the detection limit of currently available observations (e.g. Grcevich & Putman 2009; Spekkens et al. 2014), or no H I at all. This also suggest that LSB-D’s interaction with the environment of Dorado, either began recently or has not been sufficiently disruptive to deplete the galaxy of H I.

Different scenarios may explain the absence of galaxies similar to LSB-D close to the Milky Way or Andromeda. Those low surface brightness galaxies may not yet have entered the group (such as DDO210), or they could have already lost their H I because of stronger interactions, or they could exist but because of their low-column density gas (≲1019 cm−2) they have not yet been detected. Dorado is not the only environment where we find several gas-rich low-mass satellites. Until wide-field H I observations of the LG with comparable sensitivities and resolutions to MHONGOOSE will become available, it will be unclear if the evolution of Milky Way’s satellites is representative of satellites in other groups, as suggested by the Satellites Around Galactic Analogs (SAGA) Survey (Geha et al. 2017), and if they provide a good representation of dwarfs in the local universe (Weisz et al. 2011). Even though the H I of the MW may interfere with the detection of close H I rich satellites which could have low radial velocities and not have yet accreted into the galaxy (Garrison-Kimmel et al. 2014).

3.4. On the enhancement of star formation in LSB-D

LSB-D is one of the few galaxies in the 106 M stellar and gas mass range known beyond the Local Group. Its overall baryon and halo properties are typical of the dwarf population, nevertheless it presents two main peculiarities: the H I disk has column densities on average one order of magnitude lower than the H I disks of Local LSBGs (∼1018 − 19 cm−2) and in its centre UV emission indicates recent star formation. Here, we discuss on the possible evolutionary scenarios that shape LSB-D as we see it now.

LSB-D is located in the outskirts of the Dorado group and likely only recently entered in the Dorado group (see Sect. 3.1). Possibly the interaction with the group’s environment is linked to the enhancement of the star formation. With the information available we cannot connect the timescales of the infall to the triggering of the SF, but we can make some considerations on the environmental processes that may have occurred. LSB-D could be the remnant of a tidal interaction that has occurred between NGC 1581 and NGC 1566, but tidal dwarfs usually lack dark matter and the distances and especially relative velocities between LSB-D and these galaxies may suggest this is not the most likely scenario. The tidal interaction nevertheless strongly unsettled the H I disks of the NGC 1581 and NGC 1566. Several H I streams and clouds are found around both galaxies (Fig. 1). In between NGC 1581 and LSB-D (at 45 projected-kpc from the latter) an H I cloud is visible at the detection limit of our observations (see Sect. 2.1 and Fig. B.1) with small relative velocity to LSB-D (Δv = 78 km s−1). The location of LSB-D on the BHMR (Sect. 3.3) suggests that the galaxy has a high fraction of baryons in its halo, ∼50% of which are in the gas phase. It is possible that while entering the group, LSB-D caught in its halo part of the H I clouds remnant of the interaction. This H I is now condensing onto the irregular disk we see now. This would explain why the disk may not fully rotating, the higher H I line-width in the centre and the extended tail in the NW may be traces of this condensation. While condensing, part of the H I in the centre may have triggered a new episode of SF. Nevertheless, very little is known about the evolution of these low-mass systems, numerical simulations have also shown that H I disks of this mass can be retained by the galaxy’s halo and re-settle in a rotating disk triggering a new episode of star-formation (Rey et al. 2024). In this picture, LSB-D could also have retained its H I within its disk over cosmic time and we are seeing a phase of re-settlement and condensation of its disk.

Even though the origin of LSB-D’s H I disk is unclear, different indications predict that it will be short-lived. Star formation feedback can disrupt an H I disk similar to LSB-D in a few hundred Myr (Rey et al. 2024), and the interaction with the intergalactic medium can also strip the gas off the galaxy. While most common in the dense environment of clusters, also in groups and the field ram pressure (RP) can be an efficient hydro-dynamical mechanism to strip galaxies (e.g. Poggianti et al. 2017; Vulcani et al. 2021). RP is likely removing H I from dwarf MW satellites and in the Sculptor group (Westmeier et al. 2011, 2015), possibly triggering SF in the dwarf galaxies Pegasus, Wolf-Lundmark-Melotte and Holmberg II in the outskirts (∼900 kpc) of the LG and M81 groups (e.g. McConnachie et al. 2007; Ianjamasimanana et al. 2020; Yang et al. 2022; Bureau & Carignan 2002). In the Fornax and Virgo clusters, RP could have depleted of their gas dwarf galaxies with similar dynamical masses to LSB-D (Boselli et al. 2008a; Venhola et al. 2022; Kleiner et al. 2023; Boselli et al. 2023). This would explain the abundance of dwarf gas poor ellipticals compared to other types of galaxies (e.g. Boselli et al. 2008b; Venhola et al. 2022). Interestingly, Dorado has tens of known low surface brightness galaxies within 150 projected-kpc from the two central galaxies (Carrasco et al. 2001). None are known to have gas, and the two that were in the MeerKAT field of view had no detectable H I. This may suggest that the hydro-dynamical interaction with the intra-group medium of Dorado has an effect on its low-mass members.

RP may already be shaping the H I disk of LSB-D. Ram pressure exerted by the IGM would explain the offset H I disk with respect to the optical body, its elongated tail in the direction opposite of infall, and possibly the recent enhancement of the SF. RP is efficient when the pressure (PRP) exerted by the IGM on a galaxy moving at a relative velocity, v, overcomes the pressure imposed on the gas disk by a galaxy’s own gravitational potential (Pgrav, e.g. Gunn et al. 1972):

P RP = ρ IGM v 2 2 π G Σ tot Σ gas P grav $$ \begin{aligned} P_{\rm RP} =\rho _{\rm IGM}v^2\gtrsim 2\pi G\Sigma _{\rm tot} \Sigma _{\rm gas} \simeq P_{\rm grav} \end{aligned} $$(4)

where ρIGM is the IGM density, given by the product of the number of atoms in the IGM and the mean mass of a particle in the ionised medium, μ = 0.75mp, Σtot = Mtot/(πRtot2) is the total surface density (stars plus gas) and Σgas is the ISM surface density of the disk. At the outermost edge of LSB-D’s disk (rH I = 2.1 kpc), the gas has column density of 2.8 × 1018 cm−2, which corresponds to ΣHI = 0.022 M pc−2. Correcting for helium abundance we obtain Σgas = 0.029 M pc−2. As illustrated in Sect. 3.2, the total baryonic mass of LSB-D is Mtot = Mbar = 4.43 × 106 M. Hence, ram pressure would be effective on LSB-D if it was greater than PRP ≳ 2.9 × 10−15 dyn cm−2.

In Sect. 3.1 we estimated the median deprojected velocity LSB-D’s velocity of LSB-D based on a set of N-body simulation of satellites within a Dorado-like group, ( v = 509 242 208 $ v = 509^{208}_{-242} $ km s−1). Assuming this velocity, for RP to take effect onto LSB-D the IGM of Dorado would have to have densities nIGM = ρIGM ≳ 3 × 10−7 cm−3. In these low density environments (similar to the outskirts of Virgo, Steyrleithner et al. 2020) RP has been observed to efficiently strip the gas off low H I mass galaxies (e.g. Westmeier et al. 2015; Ramatsoku et al. 2019; Boselli et al. 2023), hence is possible that also in Dorado ram-pressure from the IGM is shaping the HI disk of LSB-D.

Ram-pressure may enhance star-formation in the disks and in the stripped tails of galaxies (Poggianti et al. 2017; Vulcani et al. 2018), however beyond the Local Universe this mechanism has only been observed in galaxies of masses M ≳ 109 M (Boselli et al. 2008b; Fritz et al. 2017). In the outskirts of the Virgo cluster (8 × 10−5 ≲ nIGM ≲ 10−7 atoms cm−3), simulations have shown that SF can be caused by RP in galaxies with similar masses to LSB-D (Mbar = 106 M, MDM = 108 M) travelling through the medium with velocities as low as 200 km s−1 (Steyrleithner et al. 2020). Similarly, in the MW halo’s outskirts RP of the same order of magnitudes (≳10−15 dyn cm−2) may trigger SF in low-mass satellites (Samuel et al. 2023). Possibly, the same phenomenon is on-going in LSB-D. If that were the case, it would be the first time that this phenomenon is observed in a source of this low dynamical mass in the outskirts of a low-density group, like Dorado. Stripping must have, nevertheless, only just recently begun, otherwise in a few tens Myr the galaxy would have likely lost its H I, as is observed in the low-mass Fornax cluster (Kleiner et al. 2023). This also highlights the short-lived fate of LSB-D’s H I disk.

4. Conclusions and future prospects

In this paper we presented deep new MeerKAT 21-cm observations of the field of NGC 1566 from MHONGOOSE, which led to the discovery of the gas-rich LSBG LSB-D in the Dorado group. Combining the H I observations with deep optical images from VEGAS we determined the H I and stellar properties of one of the lowest-mass H I-rich galaxies thus far known outside the Local Group (Mbar = Mgas + M = 4.43 × 106 M, fgas = 0.48, Sect. 2.1). The H I disk has low column densities, the gas is compressed in the southern edge, and has an outer asymmetry in the opposite direction. Overall, the disk is likely supported by rotation. In the centre, where the H I line-width is higher (σ = 5.5 km s−1, Sect. 2.2), the stellar body has bluer colours and coincident UV emission indicates on-going star formation (Sect. 2.3). The SFR of LSB-D is similar to that of star-forming dwarf galaxies in the Local Group.

The gas and stellar properties of LSB-D allow us to draw a picture of its evolution within Dorado. In the past, LSB-D likely had a first episode of star formation that shaped the stellar body as we see it now. During this episode and throughout the subsequent evolution through cosmic time, LSB-D retained a high percentage of baryons and in particular H I within the halo or this gas was accreted onto its halo, as is indicated by its gas fraction and baryon-to-halo mass content (see Fig. 10). Possibly, the H I disk has not been consumed by further star formation because its column densities are too low for the gas to cool into dense molecular clouds and generate a new episode of star formation (Kanekar et al. 2011; Maccagni et al. 2017). LSB-D evolved mostly unperturbed until recent times, when it began its first infall into the Dorado group (Sect. 3.1). During its infall, LSB-D did not have a major interaction with other nearby galaxies but likely the environment caused the formation of the H I disk and the recent enhancement of SF. In Sect. 3.4, we suggested that on the one hand, the orbit within the group may have caused the H I in the halo to condense in the disk, and start a new epoch of star formation, or LSB-D may have caught some sparse H I clouds from Dorado’s IGM which are now forming a disk, and stars. On the other hand, RP from the low-density IGM may also be stripping the H I disk, while compressing it in the centre and possibly triggering the new episode of SF (Sect. 3.4). In all scenarios, LSB-D’s H I disk will likely be short-lived and survive only a few hundred million years.

In the future, Integral Field Unit observations in the optical band may be able to shed new light on the star formation history of LSB-D and provide information on its stellar kinematics. The high-resolution (≲1″) study of the stellar body of LSB-D (and of its kinematics with respect to the H I), as well as of its ionised gas and chemical composition, will enable us to obtain a precise measurement of its star formation rate, determine which phenomena affect its multi-phase ISM, and study its evolution within Dorado.

Given its very low mass, the discovery of LSB-D beyond the Local Group, and the detailed study of its H I and stellar properties, demonstrates that the combination of MHONGOOSE observations with deep optical images opens a new parameter space of studies of LSBGs. In the observations presented in this paper, four other sources with H I mass lower than 2 × 107 M have also been detected in the proximity of NGC 1566. MHONGOOSE is observing ≈24 deg2 with high sensitivity to low-density gas and low H I masses, and is commensally observing the same fields with the VST with high-surface-brightness sensitivities. This will enable the discovery and detailed analysis of several more low-mass gas-rich galaxies in the field environment. The first data release (DR1, with five hours of observing time per target, 10 $ \sqrt{10} $ less sensitive than the observations shown in this paper) already discovered five new sources (de Blok et al. 2024), amongst the 46 additional H I satellites detected around its targets. The full survey with ten times the observing time of MHONGOOSE-DR1 will certainly discover more objects like LSB-D.

The studies of ultra-low-mass H I rich galaxies have so far been limited to the Local Group, which alone cannot provide a representative picture of the evolution of dwarf galaxies. The study of LSB-D, and the many more that will be discovered in the future, will make possible Local Group cosmological studies in different environments in the nearby Universe. Extending the study done on LSB-D to all newly discovered sources will allow us to populate the low-mass end of the BHMR and study the scatter (and its origin) in the BTF relation. Moreover, these studies enable a direct comparison with the kinematics, gas content, and environmental effects of dwarf satellite galaxies present in large-volume hydrodynamical simulations (Rohr et al. 2023; Samuel et al. 2023; Jones et al. 2024), which will quantify the main differences between the observed and simulated galaxy assembly processes, thus enabling us to provide observational constraints for a more detailed description of these phenomena in the next generation of simulations.


1

The datacubes span approximately 720 channels (or a 1000 km s−1 velocity range) centred on the target galaxy.

2

At this distance, the image scale is 85 parsec/arcsec assuming a ΛCDM cosmology, with Hubble constant H0 = 70 km s−1 Mpc−1 and ΩΛ = 0.7 and ΩM = 0.3.

3

In this analysis we used the same samples of Fig. 9 for which measurements of the SFR and gas and stellar masses were available; SFR for THINGS are taken from Leroy et al. (2008).

4

fbar, the ratio between their baryonic (Mbar) to halo (M200) mass normalized by the cosmological baryon fraction fbar, cosmic = 0.16.

Acknowledgments

The authors thank the anonymous referee for the useful comments and suggestions. We thank Dr. Teymoor Saifollahi for giving the suggestion to write this paper. This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement no. 882793). PEMP acknowledges the support from the Dutch Research Council (NWO) through the Veni grant VI.Veni.222.364 KAO acknowledges support by the Royal Society through a Dorothy Hodgkin Fellowship (DHF\R1\231105) and by STFC through grant ST/T000244/1. This work used the DiRAC@Durham facility managed by the Institute for Computational Cosmology on behalf of the STFC DiRAC HPC Facility (https://www.dirac.ac.uk). The equipment was funded by BEIS capital funding via STFC capital grants ST/K00042X/1, ST/P002293/1, ST/R002371/1 and ST/S002502/1, Durham University and STFC operations grant ST/R000832/1. DiRAC is part of the National e-Infastructure. AB acknowledges support from the Centre National d’Etudes Spatial (CNES), France. MK acknowledges support from the National Research Foundation of Korea (NRF) grant funded by the Korea government (Ministry of Science and ICT: MSIT) (No. NRF-RS-2023-00243222, NRF-RS-2022-00197685). This paper makes use of MeerKAT data. 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. (Part of) the data published here have been reduced using the CARACal pipeline, partially supported by ERC Starting grant number 679627 “FORNAX”, MAECI Grant Number ZA18GR02, DST-NRF Grant Number 113121 as part of the ISARP Joint Research Scheme, and BMBF project 05A17PC2 for D-MeerKAT. Information about CARACal can be obtained online under the URL: https://caracal.readthedocs.io.

References

  1. Adams, E. A. K., & Oosterloo, T. A. 2018, A&A, 612, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Adams, E. A. K., Adebahr, B., de Blok, W. J. G., et al. 2022, A&A, 667, A38 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Anand, G. S., Lee, J. C., Van Dyk, S. D., et al. 2021, MNRAS, 501, 3621 [Google Scholar]
  4. Begum, A., Chengalur, J. N., Karachentsev, I. D., & Sharina, M. E. 2008, MNRAS, 386, 138 [CrossRef] [Google Scholar]
  5. Bernstein-Cooper, E. Z., Cannon, J. M., Elson, E. C., et al. 2014, AJ, 148, 35 [NASA ADS] [CrossRef] [Google Scholar]
  6. Bianchi, L., & Garcia, M. 2014, Adv. Space Res., 53, 973 [NASA ADS] [CrossRef] [Google Scholar]
  7. Boselli, A., Boissier, S., Cortese, L., & Gavazzi, G. 2008a, ApJ, 674, 742 [Google Scholar]
  8. Boselli, A., Boissier, S., Cortese, L., & Gavazzi, G. 2008b, A&A, 489, 1015 [CrossRef] [EDP Sciences] [Google Scholar]
  9. Boselli, A., Serra, P., de Gasperin, F., et al. 2023, A&A, 676, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Bosma, A. 1978, PhD Thesis, University of Groningen, The Netherlands [Google Scholar]
  11. Bouchard, A., Carignan, C., & Staveley-Smith, L. 2006, AJ, 131, 2913 [NASA ADS] [CrossRef] [Google Scholar]
  12. Briggs, D. S., Schwab, F. R., & Sramek, R. A. 1999, ASP Conf. Ser., 180, 127 [NASA ADS] [Google Scholar]
  13. Buck, T., Macciò, A. V., Dutton, A. A., Obreja, A., & Frings, J. 2019, MNRAS, 483, 1314 [NASA ADS] [CrossRef] [Google Scholar]
  14. Bureau, M., & Carignan, C. 2002, AJ, 123, 1316 [NASA ADS] [CrossRef] [Google Scholar]
  15. Calette, A. R., Rodríguez-Puebla, A., Avila-Reese, V., & Lagos, C. d. P. 2021, MNRAS, 506, 1507 [NASA ADS] [CrossRef] [Google Scholar]
  16. Camilo, F., Scholz, P., Serylak, M., et al. 2018, ApJ, 856, 180 [Google Scholar]
  17. Capaccioli, M., Spavone, M., Grado, A., et al. 2015, A&A, 581, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Carrasco, E. R., Mendes de Oliveira, C., Infante, L., & Bolte, M. 2001, AJ, 121, 148 [NASA ADS] [CrossRef] [Google Scholar]
  19. Das, M., McGaugh, S. S., Ianjamasimanana, R., Schombert, J., & Dwarakanath, K. S. 2020, ApJ, 889, 10 [Google Scholar]
  20. de Blok, W. J. G., Adams, E. A. K., Amram, P., et al. 2016, MeerKAT Science: On the Pathway to the SKA, 7 [Google Scholar]
  21. de Blok, W. J. G., Athanassoula, E., Bosma, A., et al. 2020, A&A, 643, A147 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. de Blok, W. J. G., Healy, J., Maccagni, F. M., et al. 2024, A&A, 688, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Dey, A., Schlegel, D. J., Lang, D., et al. 2019, AJ, 157, 168 [Google Scholar]
  24. Duc, P. A., & Mirabel, I. F. 1998, A&A, 333, 813 [Google Scholar]
  25. Dutton, A. A., & Macciò, A. V. 2014, MNRAS, 441, 3359 [Google Scholar]
  26. Elagali, A., Staveley-Smith, L., Rhee, J., et al. 2019, MNRAS, 487, 2797 [NASA ADS] [CrossRef] [Google Scholar]
  27. Firth, P., Evstigneeva, E. A., Jones, J. B., et al. 2006, MNRAS, 372, 1856 [NASA ADS] [CrossRef] [Google Scholar]
  28. Forbes, D. A., Dullo, B. T., Gannon, J., et al. 2020, MNRAS, 494, 5293 [Google Scholar]
  29. Fritz, J., Moretti, A., Gullieuszik, M., et al. 2017, ApJ, 848, 132 [NASA ADS] [CrossRef] [Google Scholar]
  30. Garrison-Kimmel, S., Boylan-Kolchin, M., Bullock, J. S., & Lee, K. 2014, MNRAS, 438, 2578 [NASA ADS] [CrossRef] [Google Scholar]
  31. Geha, M., Blanton, M. R., Masjedi, M., & West, A. A. 2006, ApJ, 653, 240 [NASA ADS] [CrossRef] [Google Scholar]
  32. Geha, M., Wechsler, R. H., Mao, Y.-Y., et al. 2017, ApJ, 847, 4 [NASA ADS] [CrossRef] [Google Scholar]
  33. Gil de Paz, A., Boissier, S., Madore, B. F., et al. 2007, ApJS, 173, 185 [Google Scholar]
  34. Giovanelli, R., Haynes, M. P., Adams, E. A. K., et al. 2013, AJ, 146, 15 [NASA ADS] [CrossRef] [Google Scholar]
  35. Gray, L. M., Rhode, K. L., Leisman, L., et al. 2023, AJ, 165, 197 [NASA ADS] [CrossRef] [Google Scholar]
  36. Grcevich, J., & Putman, M. E. 2009, ApJ, 696, 385 [NASA ADS] [CrossRef] [Google Scholar]
  37. Gunn, J. E., Gott, J., & Richard, I. 1972, ApJ, 176, 1 [Google Scholar]
  38. Honey, M., van Driel, W., Das, M., & Martin, J. M. 2018, MNRAS, 476, 4488 [NASA ADS] [CrossRef] [Google Scholar]
  39. Hunter, D. A., Elmegreen, B. G., & Ludka, B. C. 2010, AJ, 139, 447 [NASA ADS] [CrossRef] [Google Scholar]
  40. Hunter, D. A., Ficut-Vicas, D., Ashley, T., et al. 2012, AJ, 144, 134 [CrossRef] [Google Scholar]
  41. Ianjamasimanana, R., Namumba, B., Ramaila, A. J. T., et al. 2020, MNRAS, 497, 4795 [NASA ADS] [CrossRef] [Google Scholar]
  42. Into, T., & Portinari, L. 2013, MNRAS, 430, 2715 [Google Scholar]
  43. Iodice, E., Cantiello, M., Hilker, M., et al. 2020, A&A, 642, A48 [EDP Sciences] [Google Scholar]
  44. Iodice, E., Spavone, M., Raj, M. A., et al. 2021, ArXiv e-prints [arXiv:2102.04950] [Google Scholar]
  45. Iorio, G., Fraternali, F., Nipoti, C., et al. 2017, MNRAS, 466, 4159 [NASA ADS] [Google Scholar]
  46. Jaffé, Y. L., Smith, R., Candlish, G. N., et al. 2015, MNRAS, 448, 1715 [Google Scholar]
  47. Jarrett, T. H., Comrie, A., Marchetti, L., et al. 2021, Astron. Comput., 37, 100502 [NASA ADS] [CrossRef] [Google Scholar]
  48. Jones, M. G., Sand, D. J., Karunakaran, A., et al. 2024, ApJ, 966, 93 [NASA ADS] [CrossRef] [Google Scholar]
  49. Józsa, G. I. G., White, S. V., Thorat, K., et al. 2020, ASP Conf. Ser., 527, 635 [Google Scholar]
  50. Józsa, G. I. G., Thorat, K., Kamphuis, P., et al. 2021, MNRAS, 501, 2704 [CrossRef] [Google Scholar]
  51. Kanekar, N., Braun, R., & Roy, N. 2011, ApJ, 737, L33 [NASA ADS] [CrossRef] [Google Scholar]
  52. Kilborn, V. A., Koribalski, B. S., Forbes, D. A., Barnes, D. G., & Musgrave, R. C. 2005, MNRAS, 356, 77 [CrossRef] [Google Scholar]
  53. Kleiner, D., Serra, P., Maccagni, F. M., et al. 2023, A&A, 675, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Kong, D., Kaplinghat, M., Yu, H.-B., Fraternali, F., & Mancera Piña, P. E. 2022, ApJ, 936, 166 [NASA ADS] [CrossRef] [Google Scholar]
  55. Koribalski, B. S., Wang, J., Kamphuis, P., et al. 2018, MNRAS, 478, 1611 [NASA ADS] [CrossRef] [Google Scholar]
  56. Kourkchi, E., & Tully, R. B. 2017, ApJ, 843, 16 [NASA ADS] [CrossRef] [Google Scholar]
  57. Kuijken, K. 2011, The Messenger, 146, 8 [NASA ADS] [Google Scholar]
  58. Lee, J. C., de Paz, A. G., Tremonti, C., et al. 2009, ApJ, 706, 599 [NASA ADS] [CrossRef] [Google Scholar]
  59. Leisman, L., Haynes, M. P., Janowiecki, S., et al. 2017, ApJ, 842, 133 [NASA ADS] [CrossRef] [Google Scholar]
  60. Lelli, F. 2022, Nat. Astron., 6, 35 [NASA ADS] [CrossRef] [Google Scholar]
  61. Lelli, F., Verheijen, M., & Fraternali, F. 2014, A&A, 566, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Lelli, F., McGaugh, S. S., & Schombert, J. M. 2016, AJ, 152, 157 [Google Scholar]
  63. Lelli, F., McGaugh, S. S., Schombert, J. M., Desmond, H., & Katz, H. 2019, MNRAS, 484, 3267 [NASA ADS] [CrossRef] [Google Scholar]
  64. Leroy, A. K., Walter, F., Brinks, E., et al. 2008, AJ, 136, 2782 [Google Scholar]
  65. Leroy, A. K., Sandstrom, K. M., Lang, D., et al. 2019, ApJS, 244, 24 [Google Scholar]
  66. Maccagni, F. M., Morganti, R., Oosterloo, T. A., Geréb, K., & Maddox, N. 2017, A&A, 604, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  67. Maccagni, F. M., Ruffa, I., Loni, A., et al. 2023, A&A, 675, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  68. Mancera Piña, P. E., Fraternali, F., Adams, E. A. K., et al. 2019, ApJ, 883, L33 [CrossRef] [Google Scholar]
  69. Mancera Piña, P. E., Posti, L., Fraternali, F., Adams, E. A. K., & Oosterloo, T. 2021, A&A, 647, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Mancera Piña, P. E., Fraternali, F., Oosterloo, T., et al. 2022a, MNRAS, 514, 3329 [CrossRef] [Google Scholar]
  71. Mancera Piña, P. E., Fraternali, F., Oosterloo, T., et al. 2022b, MNRAS, 512, 3230 [CrossRef] [Google Scholar]
  72. Marasco, A., Belfiore, F., Cresci, G., et al. 2023, A&A, 670, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  73. Martínez-Delgado, D., Cooper, A. P., Román, J., et al. 2023, A&A, 671, A141 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  74. Mauch, T., Cotton, W. D., Condon, J. J., et al. 2020, ApJ, 888, 61 [Google Scholar]
  75. McConnachie, A. W., Venn, K. A., Irwin, M. J., Young, L. M., & Geehan, J. J. 2007, ApJ, 671, L33 [NASA ADS] [CrossRef] [Google Scholar]
  76. McFarland, J. P., Verdoes-Kleijn, G., Sikkema, G., et al. 2013, Exp. Astron., 35, 45 [Google Scholar]
  77. McGaugh, S. S. 2011, Phys. Rev. Lett., 106, 121303 [CrossRef] [Google Scholar]
  78. McGaugh, S. S. 2012, AJ, 143, 40 [Google Scholar]
  79. McGaugh, S. S., Schombert, J. M., Bothun, G. D., & de Blok, W. J. G. 2000, ApJ, 533, L99 [Google Scholar]
  80. McGaugh, S. S., Schombert, J. M., & Lelli, F. 2017, ApJ, 851, 22 [Google Scholar]
  81. McGaugh, S. S., Lelli, F., Schombert, J. M., et al. 2021, AJ, 162, 202 [NASA ADS] [CrossRef] [Google Scholar]
  82. McQuinn, K. B. W., Mitchell, N. P., & Skillman, E. D. 2015a, ApJS, 218, 29 [NASA ADS] [CrossRef] [Google Scholar]
  83. McQuinn, K. B. W., Skillman, E. D., Dolphin, A. E., & Mitchell, N. P. 2015b, ApJ, 808, 109 [NASA ADS] [CrossRef] [Google Scholar]
  84. Meyer, M., Robotham, A., Obreschkow, D., et al. 2017, PASA, 34, 52 [Google Scholar]
  85. Moster, B. P., Naab, T., & White, S. D. M. 2013, MNRAS, 428, 3121 [Google Scholar]
  86. Navarro, J. F., Frenk, C. S., & White, S. D. M. 1996, ApJ, 462, 563 [Google Scholar]
  87. Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 [Google Scholar]
  88. Oh, S. H., Kim, W.-T., & Lee, H. M. 2015, ApJ, 807, 73 [NASA ADS] [CrossRef] [Google Scholar]
  89. Oman, K. A., Hudson, M. J., & Behroozi, P. S. 2013, MNRAS, 431, 2307 [Google Scholar]
  90. Oman, K. A., Bahé, Y. M., Healy, J., et al. 2021, MNRAS, 501, 5073 [CrossRef] [Google Scholar]
  91. Ott, J., Stilp, A. M., Warren, S. R., et al. 2012, AJ, 144, 123 [NASA ADS] [CrossRef] [Google Scholar]
  92. Planck Collaboration VI. 2020, A&A, 641, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  93. Poggianti, B. M., Moretti, A., Gullieuszik, M., et al. 2017, ApJ, 844, 48 [Google Scholar]
  94. Ponomareva, A. A., Verheijen, M. A. W., Papastergis, E., Bosma, A., & Peletier, R. F. 2018, MNRAS, 474, 4366 [Google Scholar]
  95. Posti, L., Fraternali, F., & Marasco, A. 2019, A&A, 626, A56 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  96. Ragusa, R., Spavone, M., Iodice, E., et al. 2021, A&A, 651, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  97. Ragusa, R., Mirabile, M., Spavone, M., et al. 2022, Front. Astron. Space Sci., 9, 852810 [CrossRef] [Google Scholar]
  98. Ragusa, R., Iodice, E., Spavone, M., et al. 2023, A&A, 670, L20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  99. Ramatsoku, M., Serra, P., Poggianti, B. M., et al. 2019, MNRAS, 487, 4580 [Google Scholar]
  100. Rampazzo, R., Ciroi, S., Mazzei, P., et al. 2020, A&A, 643, A176 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  101. Rampazzo, R., Mazzei, P., Marino, A., et al. 2021, JAA, 42, 31 [NASA ADS] [Google Scholar]
  102. Rampazzo, R., Mazzei, P., Marino, A., et al. 2022, A&A, 664, A192 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  103. Read, J. I., Iorio, G., Agertz, O., & Fraternali, F. 2017, MNRAS, 467, 2019 [NASA ADS] [Google Scholar]
  104. Rey, M. P., Orkney, M. D. A., Read, J. I., et al. 2024, MNRAS, 529, 2379 [NASA ADS] [CrossRef] [Google Scholar]
  105. Rhee, J., Smith, R., Choi, H., et al. 2017, ApJ, 843, 128 [Google Scholar]
  106. Rohr, E., Pillepich, A., Nelson, D., et al. 2023, MNRAS, 524, 3502 [NASA ADS] [CrossRef] [Google Scholar]
  107. Ryan-Weber, E. V., Begum, A., Oosterloo, T., et al. 2008, MNRAS, 384, 535 [NASA ADS] [CrossRef] [Google Scholar]
  108. Sakai, S., Mould, J. R., Hughes, S. M. G., et al. 2000, ApJ, 529, 698 [NASA ADS] [CrossRef] [Google Scholar]
  109. Salim, S., Rich, R. M., Charlot, S., et al. 2007, ApJS, 173, 267 [NASA ADS] [CrossRef] [Google Scholar]
  110. Samuel, J., Pardasani, B., Wetzel, A., et al. 2023, MNRAS, 525, 3849 [NASA ADS] [CrossRef] [Google Scholar]
  111. Sancisi, R., Fraternali, F., Oosterloo, T., & van der Hulst, T. 2008, A&A Rev., 15, 189 [CrossRef] [Google Scholar]
  112. Schaye, J., Crain, R. A., Bower, R. G., et al. 2015, MNRAS, 446, 521 [Google Scholar]
  113. Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103 [Google Scholar]
  114. Serra, P., Westmeier, T., Giese, N., et al. 2015, MNRAS, 448, 1922 [Google Scholar]
  115. Serra, P., Maccagni, F. M., Kleiner, D., et al. 2019, A&A, 628, A122 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  116. Serra, P., Maccagni, F. M., Kleiner, D., et al. 2023, A&A, 673, A146 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  117. Spavone, M., Capaccioli, M., Napolitano, N. R., et al. 2017, A&A, 603, A38 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  118. Spekkens, K., Urbancic, N., Mason, B. S., Willman, B., & Aguirre, J. E. 2014, ApJ, 795, L5 [NASA ADS] [CrossRef] [Google Scholar]
  119. Steyrleithner, P., Hensler, G., & Boselli, A. 2020, MNRAS, 494, 1114 [Google Scholar]
  120. Swaters, R. A., van Albada, T. S., van der Hulst, J. M., & Sancisi, R. 2002, A&A, 390, 829 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  121. Tanoglidis, D., Drlica-Wagner, A., Wei, K., et al. 2021, ApJS, 252, 18 [Google Scholar]
  122. Venhola, A., Peletier, R., Laurikainen, E., et al. 2017, A&A, 608, A142 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  123. Venhola, A., Peletier, R. F., Salo, H., et al. 2022, A&A, 662, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  124. Verheijen, M. A. W. 2001, ApJ, 563, 694 [NASA ADS] [CrossRef] [Google Scholar]
  125. Vulcani, B., Poggianti, B. M., Gullieuszik, M., et al. 2018, ApJ, 866, L25 [Google Scholar]
  126. Vulcani, B., Poggianti, B. M., Moretti, A., et al. 2021, ApJ, 914, 27 [NASA ADS] [CrossRef] [Google Scholar]
  127. Wall, R. E., Kilic, M., Bergeron, P., et al. 2019, MNRAS, 489, 5046 [NASA ADS] [CrossRef] [Google Scholar]
  128. Walter, F., Brinks, E., de Blok, W. J. G., et al. 2008, AJ, 136, 2563 [Google Scholar]
  129. Warren, S. R., Skillman, E. D., Stilp, A. M., et al. 2012, ApJ, 757, 84 [NASA ADS] [CrossRef] [Google Scholar]
  130. Weisz, D. R., Dolphin, A. E., Dalcanton, J. J., et al. 2011, ApJ, 743, 8 [NASA ADS] [CrossRef] [Google Scholar]
  131. Westmeier, T., Braun, R., & Koribalski, B. S. 2011, MNRAS, 410, 2217 [Google Scholar]
  132. Westmeier, T., Staveley-Smith, L., Calabretta, M., et al. 2015, MNRAS, 453, 338 [NASA ADS] [CrossRef] [Google Scholar]
  133. Westmeier, T., Kitaeff, S., Pallot, D., et al. 2021, MNRAS, 506, 3962 [NASA ADS] [CrossRef] [Google Scholar]
  134. Westmeier, T., Deg, N., Spekkens, K., et al. 2022, PASA, 39, e058 [NASA ADS] [CrossRef] [Google Scholar]
  135. Xu, J.-L., Zhu, M., Yu, N., et al. 2023, ApJ, 944, L40 [NASA ADS] [CrossRef] [Google Scholar]
  136. Yang, Y., Ianjamasimanana, R., Hammer, F., et al. 2022, A&A, 660, L11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  137. Young, L. M., van Zee, L., Lo, K. Y., Dohm-Palmer, R. C., & Beierle, M. E. 2003, ApJ, 592, 111 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: MHONGOOSE datacubes sensitivities and H I detections in Dorado

As mentioned in Sect. 2, we have identified six optimal weighting and tapering combinations that span the resolution range of MeerKAT (7″–90″) at optimal sensitivities with which we generate the H I datacubes. The lowest resolution datacube (94″×92″) is generated with Briggs robust = 1.0 and taper to 90″, hence is named r10_t90. With robust = 0.5 and 60″ tapering we reach 65″×64″ resolution (r05_t60). Higher resolution datacubes are generated with different robust weights (robust = 1.5, 1.0, 0.5, 0.0) and no tapering reaching resolutions of 32″×23″(r15), 25″×23″ (r10), 12″×10″ (r05) and 7.6″×7.2″ (r00). The noise level matches the expected theoretical values, as illustrated in the survey paper (de Blok et al. 2024). The 3σ H I column density sensitivities in the centre of the cube are computed over 16 km s−1 and vary between 5 × 1017 cm−2 at 90″ resolution and 5 × 1019 cm−2 at 7″, for emission filling the beam uniformly. Figure A.1 shows the overlay of the H I flux density maps derived from the datacubes at all resolutions. New H I detections are marked in green.

Table A.1 summarizes the main properties of the datacubes. The datacube used to study LSB-D’s H I disk is r10, marked in bold in the Table. A summary of the H I properties of the sources detected by the MHONGOOSE observations is given in Table A.2.

Table A.1.

Properties of the MHONGOOSE H I datacubes of NGC 1566

Table A.2.

H I properties of the sources detected by MHONGOOSE in Dorado.

thumbnail Fig. A.1.

Primary beam corrected flux-density H I emission (in orange) detected by MeerKAT, overlaid with the 2 × 2 degrees DECaLS optical image in the g,r,z filters. The H I emission is an overlay of the flux-density maps extracted from the multi-resolution datacubes (see Sect. 2). The PSFs of the multi-resolution maps varying from 90″ to 7″ are shown in white in the bottom-right corner. The target of the MeerKAT observations, NGC 1566, is marked in orange, sources with previously known H I are marked in white, while new H I detections are marked in green. The doublet of early type galaxies in the centre of the Dorado group is marked in cyan. The newly discovered H I low surface brightness galaxy object of this paper is marked by the red circle.

Appendix B: H I clouds in Dorado

Figure B.1, shows the integrated H I profile of a ‘dark’ cloud that may be present at the detection limit of our observations at coordinates: 04h22m31s, −55d12m02s. At 25 × 18″, this emission line is unresolved, it is centred at vsys = 1292 km s−1 and has full width at zero intensity of wFWZI = 18 km s−1. It is best detected after hanning-smoothing the datacube over three channels. The H I integrated flux S1.4GHzΔv = 5.5 mJy km s−1. At the distance of Dorado, this corresponds to an H I mass of MH I = 3.6 × 105 M. This value is 20% higher than the S/N = 3 H I mass detection limit for a source of 16 km s−1 in those regions of the primary beam corrected datacube. If real this cloud may be the remnant of the past interaction between NGC 1566 and NGC 1581. The projected distance between this H I cloud and LSB-D is ∼8.86′ (∼45 kpc) and their difference in line of sight velocities is Δv = 78 km s−1. Even though rare, similar ‘dark’ low-column density H I clouds (NH I ≲3 × 1019 cm−2) have been detected as remnants of past interactions in groups (see, for example, Serra et al. 2019; Maccagni et al. 2023). Some ‘dark’ clouds also show indications of rotation (e.g. Xu et al. 2023), which could also occur in this newly detected cloud, as suggested by the double peaked profile.

In Fig. B.2 we show the channel maps of the H I emission of LSB-D detected at 25″ × 18″and 32″ × 23″resolution. The centre of the optical body of the galaxy is marked by the red cross. One every two channel is shown, the H I disk is visible between 1207 km s−1 and 1220 km s−1, while the cloud in the north-east is visible at channels 1209 km s−1 and 1212 km s−1.

thumbnail Fig. B.1.

Integrated spectrum the ‘dark’ cloud detected at coordinates 04h22m31s, −55d12m02s. The vertical dashed-dotted line shows the central velocity of the H I line ( v sys radio = 1292 $ v^{\mathrm{radio}}_{\mathrm{sys}} = 1292 $ km s−1), while the grey vertical lines marks its full velocity range (wFWZI = 18 km s−1).

thumbnail Fig. B.2.

Channel maps of the H I emission of LSB-D maps from the r10 and r15 datacubes (orange and blue contours, respectively). One every two channels is shown, contour levels increase as 2.5σ × 2n (n = 0, 1). The red cross marks the centre of the optical body of the galaxy. The PSFs of the datacubes are shown in the bottom left corner.

Appendix C: Deprojected coordinates of LSB-D derived from orbit catalogues

In Sect. 3 we provide the probability for LSB-D to be a recent infaller of the Dorado group based on the probability distribution of its orbital parameters and deprojected coordinates of galaxies in groups and clusters of an N-body simulation. Here, were describe more in detail the methodology which follows Oman et al. (2021). The input orbit catalogue is compiled from an N-body simulation of the same (100 Mpc)3 volume used in the EAGLE (Schaye et al. 2015) project (with a particle mass of 1.2 × 107 M) run to a final scale factor of two. Dorado’s virial mass and radius (see Table 1) correspond to an overdensity Δ vir = M vir / ( ρ crit 4 3 π r vir 3 ) = 220 $ \Delta_{\mathrm{vir}}=M_{\mathrm{vir}} / (\rho_{\mathrm{crit}}\frac{4}{3}\pi r_{\mathrm{vir}}^3) = 220 $. Assuming a Navarro et al. (1997) density profile and a concentration parameter c200 = 10 we convert to the virial definition used in Oman et al. (2021) and calculate the corresponding 3D velocity dispersion using their Eq. 1. The values used in this context are: Mvir = 4.1 × 1013 M; rvir = 873 kpc; σ3D = 571 km s−1. Assuming NGC 1553 as the group centre gives normalised offsets of r/rvir = 0.40 and (vsys − vGroup)/σ3D = 0.02. We select satellites from the orbit catalogue where: (i) the host mass is within 0.2 dex of Mvir; (ii) the normalised radial offset is within ±0.1 of that of LSB-D; (iii) the absolute value of the normalised LoS velocity offset is less than 0.12 (the absolute value reflects the degeneracy between the sign of the velocity offset and whether LSB-D is in the foreground or background of the group); (iv) the satellite halo peak mass (maximum bound mass at any earlier time) is within 0.5 dex of 108 M. The last element of the selection, (iv), is the most uncertain, but the H I kinematics and stellar mass of LSB-D suggest a total mass of about 108 M – and the probability distributions obtained from the orbit library are only weakly sensitive to this choice.

thumbnail Fig. C.1.

Deprojected coordinates of LSB-D derived from orbit catalogues following the methodology of Oman et al. (2021) as detailed in Sec. 3.1. The grayscale map shows a histogram of the deprojected coordinates for simulated satellites with projected coordinates, host mass and satellite mass compatible with those of LSB-D. The red lines mark the median and 16th-84th percentiles of the marginalised probability distributions for the deprojected radius r3D and current speed v3D.

All Tables

Table 1.

Properties of the Dorado group.

Table 2.

H I properties of LSB-D.

Table 3.

Stellar properties of LSB-D.

Table A.1.

Properties of the MHONGOOSE H I datacubes of NGC 1566

Table A.2.

H I properties of the sources detected by MHONGOOSE in Dorado.

All Figures

thumbnail Fig. 1.

Primary beam corrected flux-density H I emission detected by MeerKAT within a 1.5 × 1.5 deg2 field centred on NGC 1566 (the imaged f.o.v. is marked by the dashed grey circle), overlaid on the DECaLS optical image in the g, r, z filters. The H I emission is a composite of flux-density maps derived from the multi-resolution datacubes. The resolutions shown by coloured contours are 12″ × 10″ (red), 25″×18″ (green), 32″×23″ (blue), 65″×64″ (purple), 94″×92″ (yellow-green). For resolutions between ∼90″ and 25″ two contours are shown, i.e., NH I = 3σ × 2n with n = 0, 2. The 12″×10″ resolution has further increasing contours (n = 0, 2, 4...). The PSF of the multi-resolution maps are shown in the bottom-right corner with the same colour coding. Details on the H I properties of the sources in the field are given in Table A.2.

In the text
thumbnail Fig. 2.

Integrated spectrum of LSB-D, with its noise per channel (grey shaded region). The vertical dashed line shows the systemic velocity of the source ( v sys radio = 1214 $ v^{\mathrm{radio}}_{\mathrm{sys}} = 1214 \, $km s−1), while the grey vertical lines show its full velocity range (wFWZI = 24 km s−1).

In the text
thumbnail Fig. 3.

Left Panel: DeCaLS optical image in the g, r, z filters overlaid with the H I contours from the 25″ × 18″ datacube used for the analysis in this paper. Levels increase as 4.6 × 1018 × 2n cm−2 (n = 0, 1, 2, 3), where the first contour marks the mean S/N = 3 detection limit. The PSF is shown in the bottom-right corner. Note the H I overlay in color with also the 32 × 23 ″ map. Right Panel: velocity field of the H I gas in LSB-D. The systemic velocity (vsys = 1214 km s−1) is marked by the thick black iso-velocity contour, the other contour levels are shown in the bottom-left corner. The dashed and dashed-dotted lines mark the directions of the major (141°) and minor axes (51°) of the H I disk, respectively.

In the text
thumbnail Fig. 4.

Top Panel: position-velocity diagram taken along the major axis of the H I disk (PA = 141°) with a slit of 20″. Contour levels are at 5 × 1017 cm−2 and 1 × 1018 cm−2. Bottom Panel: position-velocity diagram taken along the minor axis of the H I disk (PA = 51°) with an analogous slit, equal to the resolution of the observations. Contours are as in the top panel.

In the text
thumbnail Fig. 5.

VEGAS optical image in the g, r filters overlaid with the H I column density map. Orange contour levels are as the left panel of Fig. 3. In cyan and magenta contours we show the S/N = 3 and 5 contours from the NUV and FUV GALEX observations, respectively.

In the text
thumbnail Fig. 6.

Top panel: LSB-D’s g and r radial surface brightness profile (cyan and red, respectively) inferred from the VEGAS observations. Bottom panel: radial g − r colour profile measured out to 8″ from the centre (∼670 pc).

In the text
thumbnail Fig. 7.

Phase-space diagram of the Dorado group. LSB-D is shown in red, while previously known H I rich galaxies are shown in blue and new H I detections in green. The group’s centre (NGC 1553) is at the zero coordinates, the dashed-dotted lines mark the virialised region of the Dorado group. The caustic curves (dashed lines) are defined by the escape velocity of the group.

In the text
thumbnail Fig. 8.

Top panel: gas mass versus stellar mass for a sample of galaxies in the Local Universe. LITTLE-THINGS, SPARC and THINGS galaxies are marked by triangles (in orange, grey, and blue, respectively). Local low surface brightness galaxies (LSBGs) from McGaugh (2011, 2017, 2021) are in green, and UDGs in cyan. LSB-D is marked by the red circle. The dashed line marks the 1 : 1 linear relation. Bottom panel: the baryonic Tully-Fisher relation of Local Group dwarf and spiral galaxies (dashed line; McGaugh et al. 2021). Symbols are as in the top panel. For the Local LSBGs sample, pressure supported dwarfs (plotted against their velocity dispersion, σv) are marked by light green squares, while green triangles mark rotationally supported galaxies, plotted against their rotational velocity, vrot.

In the text
thumbnail Fig. 9.

Top panel: star formation rate as a function of stellar mass for LSB-D (marked by the red circle, and for a sample of galaxies in the Local Universe. LITTLE-THINGS (Oh et al. 2015; Iorio et al. 2017) and THINGS galaxies are marked by orange and blue triangles. Local low surface brightness galaxies (LSBGs) from McGaugh (2011, 2017) are in green. The dashed line shows a linear relation fitted betwee the LSBGs and THINGS samples. Central panel: specific star formation rate as a function of stellar mass for the same sample as in the top panel. The dashed-dotted line marks the inverse Hubble time, i.e. galaxies on the relation reach their current stellar mass with their average SFR since the Big Bang. Bottom panel: depletion time as a function of stellar mass for LSB-D and the samples of galaxies in the Local Universe. The dashed line marks the linear fit relation for the Local LSBGs sample.

In the text
thumbnail Fig. 10.

Mbar/M200 ratio normalized to the cosmological baryon fraction (fbar, cosmic). LSB-D is compared against a set of nearby disk galaxies (Mancera Piña et al. 2022a) and the semi-empirical relation from Calette et al. (2021), which uses abundance matching to link M with M200 and empirical correlations for Mstar/Mgas. Galaxies are colour coded according to their gas-to-baryon content (Mgas/Mbar).

In the text
thumbnail Fig. A.1.

Primary beam corrected flux-density H I emission (in orange) detected by MeerKAT, overlaid with the 2 × 2 degrees DECaLS optical image in the g,r,z filters. The H I emission is an overlay of the flux-density maps extracted from the multi-resolution datacubes (see Sect. 2). The PSFs of the multi-resolution maps varying from 90″ to 7″ are shown in white in the bottom-right corner. The target of the MeerKAT observations, NGC 1566, is marked in orange, sources with previously known H I are marked in white, while new H I detections are marked in green. The doublet of early type galaxies in the centre of the Dorado group is marked in cyan. The newly discovered H I low surface brightness galaxy object of this paper is marked by the red circle.

In the text
thumbnail Fig. B.1.

Integrated spectrum the ‘dark’ cloud detected at coordinates 04h22m31s, −55d12m02s. The vertical dashed-dotted line shows the central velocity of the H I line ( v sys radio = 1292 $ v^{\mathrm{radio}}_{\mathrm{sys}} = 1292 $ km s−1), while the grey vertical lines marks its full velocity range (wFWZI = 18 km s−1).

In the text
thumbnail Fig. B.2.

Channel maps of the H I emission of LSB-D maps from the r10 and r15 datacubes (orange and blue contours, respectively). One every two channels is shown, contour levels increase as 2.5σ × 2n (n = 0, 1). The red cross marks the centre of the optical body of the galaxy. The PSFs of the datacubes are shown in the bottom left corner.

In the text
thumbnail Fig. C.1.

Deprojected coordinates of LSB-D derived from orbit catalogues following the methodology of Oman et al. (2021) as detailed in Sec. 3.1. The grayscale map shows a histogram of the deprojected coordinates for simulated satellites with projected coordinates, host mass and satellite mass compatible with those of LSB-D. The red lines mark the median and 16th-84th percentiles of the marginalised probability distributions for the deprojected radius r3D and current speed v3D.

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.