## Further evidence for a population of dark-matter-deficient dwarf galaxies

Qi Guo ∗ 1 2 3 , , ; Huijie Hu 1 2 , ; Zheng Zheng 1 4 , ; Shihong Liao 1 3 , ; Wei Du 1 5 , ; Shude Mao 6 1 , ; Linhua Jiang ; Jing Wang ; Yingjie Peng ; Liang Gao 7 7 7 1 2 3 , , ; Jie Wang 1 2 3 , , ; Hong Wu 1 2 5 , ,

1 National Astronomical Observatories, Chinese Academy of Science, 20A Datun Road, Chaoyang

District, Beijing, 100101, China

2 University of Chinese Academy of Sciences, No. 19 A Yuquan Road, Beijing, 100049, China

3 Key Laboratory for Computational Astrophysics, National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100101, China

4 Key Laboratory of The Five-hundred-meter Aperture Spherical radio Telescope, National Astronomical Observatories, Chinese Academy of Sciences

5 Key Laboratory of Optical Astronomy, National Astronomical Observatories, Beijing, 100101, China

6 Physics Department and Tsinghua Center for Astrophysics, Tsinghua University, Beijing 100084, China

7 Kavli Institute for Astronomy and Astrophysics, Peking University, Beijing 100871, China

In the standard cosmological model, dark matter drives the structure formation and constructs potential wells within which galaxies may form. The baryon fraction in dark halos can reach the universal value (15.7%) in massive clusters and decreases rapidly as the mass of the system decreases 1,2 . The formation of dwarf galaxies is sensitive both to baryonic processes and the properties of dark matter owing to the shallow potential wells in which they form. In dwarf galaxies in the Local Group, dark matter dominates the mass content even within their optical-light half-radii (r e ∼ 1 kpc) 3,4 . However, recently it has been argued that not all dwarf galaxies are dominated by dark matter 5-7 . Here we report 19 dwarf galaxies

that could consist mainly of baryons up to radii well beyond r , at which point they are exe pected to be dominated by dark matter. Of these, 14 are isolated dwarf galaxies, free from the influence of nearby bright galaxies and high dense environments. This result provides observational evidence that could challenge the formation theory of low-mass galaxies within the framework of standard cosmology. Further observations, in particular deep imaging and spatially-resolved kinematics, are needed to constrain the baryon fraction better in such galaxies.

In the standard cosmology 8 , dark matter contributes around 26.6% of the Universe, while baryons only occupy 4.9%. The baryonic fraction is about the universal value, 15.7%, in galaxy clusters, and decreases rapidly towards lower masses. In low mass systems, the potential wells are shallow so that various baryonic processes, including background UV ionization and stellar feedback, could either prevent gas infall or expel gas out of the systems, leading to a low baryonic fraction. Observationally, it infers a typical mass-to-light ratio of 10-1000, using stellar dynamics of the dwarf satellite galaxies in the Local Group 9,10 . Such values are supported by modern simulations 11 . However, recently two dwarf galaxies were found to have baryonic mass comparable to the dark matter mass up to several kiloparsecs using the dynamics of their HI gas 7,12 . Similar results were reported for two ultra-diffuse dwarf galaxies in the NGC1052 group 5,6 using the dynamics of their globular clusters, of which the NGC1052-DF4 dwarfs may have almost no dark matter at all. Such relation between dark matter and baryons is also revealed in the population of tidal dwarf galaxies 13,14 . In all cases, these galaxies reside in (relatively) high density regions, and the environmental stripping/interaction might play a role.

Here we search for baryon dominated dwarf galaxies (r-band absolute magnitude, M r &gt; -18 ) in the field, using the catalogue of Arecibo Legacy Fast ALFA (ALFALFA) α .40 15 , which matches the HI selected galaxies to the Seventh Data Release of Sloan Digital Sky Survey (SDSS DR7) 16 PhotoPrimary optical catalogue. Signal-to-noise cut at HI spectrum, SNR HI &gt; 10 , is im-

posed. We further remove sources with suspicious optical images and HI spectra, and sources which have multiple Galaxy Evolution Explorer(GALEX) 17 images. The photometric and geometric properties are re-calculated for low surface brightness galaxies for SDSS DR7 pipeline systematically over-estimate the sky background for such systems 18,19 . In the lack of HI image, we use the r-band b/a ratio to correct for the inclination. In order to minimize the mis-alignment effect between the HI content and the r-band image, we select galaxies with r-band b/a ratio between 0.3 and 0.6, leading to a final parent sample of 324 galaxies. The stellar mass of each galaxy is derived from the optical photometry from SDSS and the distance measured by the ALFALFA 20 . Dynamical mass is estimated using the circular velocity V HI which is derived from 20% of the HI line width (w 20 ), and the HI radius (r HI , defined at HI surface density = 1M glyph[circledot] pc -2 ) which is converted from the HI mass based on the tight relation between HI radius and HI mass 21 . We verify the accuracy of such defined r HI and V HI using dwarf galaxies in the Local Irregulars That Trace Luminosity Extremes The HI Nearby Galaxy Survey (LITTLE THINGS), where the HI images have been well resolved and HI dynamics been detailed analyzed. Details are given in the Method.

Figure 1 shows images and spectra of two example dwarf galaxies: AGC 213086 and AGC 8915. The r-band luminosity of AGC 213086 is 3 07 . × 10 42 erg s -1 , corresponding to a stellar mass of 5 51 . +4 02 . -2 32 . × 10 8 M . glyph[circledot] The HI mass is 2 45 . +0 11 . -0 09 . × 10 9 M , corresponding to a radius of glyph[circledot] r HI = 14 37 . ± 1 023 . kpc. This radius is much larger than the half light radius in the optical wavelength, r e ∼ 3 39 . kpc, extending to the dark matter dominated region if AGC 213086 were a typical dwarf galaxy. The total baryonic mass (M bary ) within r HI is 3 82 . +0 34 . -0 32 . × 10 9 M . glyph[circledot] The circular velocity is 43 47 . ± 2 40 . km/s and the dynamical mass (M dyn ) enclosed within the r HI is 6 31 . +0 89 . -0 77 . × 10 9 M . Error bars account for various uncertainties in the data and in the method glyph[circledot] for the mass estimation. Assuming a typical NFW profile 22 the virial mass (M halo , the total mass within the virial radius, R vir , within which the mean density is 200 times the critical density) is 1 41 . × 10 10 M , and is glyph[circledot] 1 01 . × 10 10 M if assuming a Burkert profile glyph[circledot] 23 . The ratio between M dyn ( &lt; r HI )

and M bary ( &lt; r HI ) is less than 2, suggesting baryon dominates over dark matter within r HI . The baryonic fraction within R vir is 0.271 (and 0.376 if assuming a Burkert profile), well above the universal value. In comparison, we show at the bottom a 'typical' dwarf galaxies, AGC 8915. It has M bary ( &lt; r HI ) = M gas ( &lt; r HI ) + M star ( &lt; r HI )= 3 14 . +0 18 . -0 17 . × 10 9 M glyph[circledot] glyph[lessmuch] M dyn ( &lt; r HI )= 2 647 . +0 20 . -0 18 . × 10 10 M glyph[circledot] and M halo = 1 83 . × 10 11 M glyph[circledot] (NFW profile ) / 1 04 . × 10 11 M glyph[circledot] (Burkert profile). Its M dyn -to-M bary ratio within r HI is ∼ 10, and the baryonic fraction within R vir is 0.0171 (NFW profile) / 0.0301 (Burkert profile), well below the universal value.

The main statistical results are presented in Figure 2. We find 19 dwarf galaxies out of 324 parent dwarf galaxies, a non-negligible fraction, with M dyn ( &lt; r HI ) &lt; 2 × M bary ( &lt; r HI ), i.e. M DM ( &lt; r HI ) &lt; M bary ( &lt; r HI ). These 19 dwarf galaxies have more extended r e distribution and r HI distribution compared to the parent sample. For each of them, r HI is much larger than r , suggesting e its dynamics is estimated at rather large radius where dark matter is expected to be dominant.

The median M dyn ( &lt; r HI )/M bary ( &lt; r HI ) (here after M dyn /M bary ) of the parent sample is 7.36 with the oneσ dispersion of [3.85, 11.03], well below the the typical values found in the dwarf galaxies in the Local Group (right panel of Figure 2). There is a clear excess of the baryonic dominated dwarf galaxies (here after BDDGs) compared to a Gaussian distribution in the log M dyn /M bary , suggesting BDDGs to be a distinct galaxy population rather than statistical outliers.

Figure 3 shows that as the parent sample, for most BDDGs M dyn ( &lt; r HI ) is dominated by gas masses, so that the potential mis-matching between the optical and HI sources could not affect the results substantially. The parent sample generally follows the same Tully-Fisher relation as spiral galaxies, yet BDDGs tends to have much lower circular velocities for given luminosities. It is the low circular velocities rather than systemic uncertainties (e.g. geometric structures, distance estimates and etc.) that leads to their small M dyn /M bary values. Assuming the NFW profile, we find the BDDGs reside in halos as small as 10 9 ∼ 10 11 M glyph[circledot] (Figure 4), consistent with the mass

range if assuming the Burkert profile. Their baryon fractions within R vir are mostly well above the universal value. This is unexpected for these systems have rather shallow potential well and are thus harder to maintain their baryons against various feedback. Hydro-dynamical simulations of galaxy formation show that interactions of galaxies at high density region could be responsible for the formation of baryonic dominated tidal dwarf galaxies, and the environmental stripping could play a role in forming the BDDGs around clusters 26,27 . However, only 5 of the 19 BDDGs are within/close to groups/clusters. The rest of the 14 BDDGs reside at distances farther than 3 times the virial radii of any surrounding groups/clusters (right panel in Figure 4). This makes the environmental stripping and origin from interaction of galaxies very unlikely. Our results suggest that a population of dwarf galaxies could form in a particular way such that much less dark matter is required than for those in the Local Group and those found in simulations.

One possible channel is that, under extreme conditions, strong supernova feedback could flatten the potential well by expelling gas. Along with this procedure dark matter could be removed to the outer region, in response to the gravity of the expelled gas. Further fueling cold gas might reach the potential well along filaments(cold flow), while dark matter grow in a more smooth and steady way. Under such special circumstance, it is possible to form BDDGs. Further studies with high resolution simulations are required to explore such possibility. A warm dark matter model could also likely lead to a relatively low density in the inner region of a given halo and thus could probably be more in line with the existence of the BDDGs.

Alternatively, the systematic errors in inclination may lead to an under-estimate of the circular velocity and thus an under-estimate of the dynamical mass, for there could be mis-alignment between the optical inclination and the inclination of the HI kinematic field. Applying the same mis-alignment distribution between the optical images 28 and the HI dynamics from the LITTLE THINGS and those from the most up-to-data simulations 29 to our sample, we find similar distribu-

tions of the M dyn /M bary , suggesting the existence of the BDDGs. In addition, most of our BDDGs are field dwarf galaxies for which less mis-alignment effect is expected. There could also exist systematic errors in optical inclinations, i.e. galaxies with high axis ratios could be bars in face-on disks. Recent work 30 shows a fraction of 10% of such cases for dwarf galaxies, thought the fraction could be lower in low dense environments where most of our BDDGs are found. Detailed HI interferometric mapping of these galaxies in the future are required to rule out such possibilities.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00046v2-figure-ce53cd8b-d4b9-4c91-9ba6-a1775651c163.jpg)

Figure 1: Image and Spectra of two galaxy examples. The square panels show the SDSS r -band image of two selected examples. Each panel spans 3' × 3'. The spectrum from ALFALFA of each object is shown next to the corresponding SDSS image, centered at the HI lines. The velocities are in unit of km/s. The spectral resolution of ALFALFA is 10 km/s. The signal-to-noise ratio is illustrated on the upper right corner of each spectrum. The upper one shows a dwarf galaxy with M bary ∼ M dyn , while the lower one shows a typical dwarf sample with M bary glyph[lessmuch] M dyn .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00046v2-figure-7edb69da-cc94-499f-94d8-f496e846ad58.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00046v2-figure-82059872-2d56-4c35-a50a-5c370acc7ce8.jpg)

- [1] White, S. D. M.; Navarro, J. F.; Evrard, A. E. &amp; Frenk, C. S. The baryon content of galaxy clusters: a challenge to cosmological orthodoxy. Natur. 366 , 429-433 (1993).
- [2] Fukugita, M.; Hogan, C. J. &amp; Peebles, P. J. E. The cosmic baryon budget. Astrophys. J. 503 , 518-530 (1998).
- [3] Walker, M. G. in Planets Stars and Stellar Systems Vol. 5 (eds Oswalt, T. D. &amp; Gilmore, G.) 1039-1089 (2013).
- [4] Strigari, L. E. et al. A common mass scale for satellite galaxies of the Milky Way. Natur. 454 , 1096-1097 (2008).

Figure 2: Comparison of baryonic mass and total mass enclosed within r HI . Left panel: Baryonic mass vs. total mass enclosed within r HI . Grey dots are the parent sample with M r &gt; -18 , and black dots with error bars are the 19 baryonic dominated dwarf galaxies. Red dash line indicates M dyn = 2 × M bary . Error bars account for various uncertainties in the data and in the method for the mass estimation. Right panel: distribution of log(M dyn /M bary ) of baryonic dominated dwarf galaxies (black histogram) and the parent sample (grey histogram). Blue curve is the Gaussian fit to the parent sample, with the mean value of 0.91 and one σ dispersion of 0.18. A clear excess of BDDGs is presented.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00046v2-figure-dd6f8449-063d-4980-a363-0897bd328f45.jpg)

Figure 3: Galaxy properties. Left panel : Gas mass vs. stellar mass relation. Grey dots are for the parent sample, and black dots are for the baryonic dominated dwarf galaxies. Most of the parent galaxies have gas masses dominant over stellar masses. Their low M dyn /M bary compared to those found in the dwarf galaxies in the Local Group can thus be partly explained by the removal of gas at high density region. Right panel : r-band Tully-Fisher relation. V HI is the rotation velocity measured using w 20 which also takes into account the inclination and the velocity dispersion. The Tully-Fisher relation for more massive galaxies are shown with blue dots with error bars 24,25 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00046v2-figure-b25e19a3-30fc-48bf-bdb7-71e32adf602b.jpg)

Figure 4: Halo mass and distance to the nearest clusters/groups. Left panel : Dynamical mass enclosed within r HI vs. r HI for baryon dominated galaxies. Black symbols are our 19 BDDGs. Solid and dashed curves are the typical dark halo mass profiles estimated assuming the NFW profile and the Burkert profile, ranking from 10 11 , 10 10 , 10 9 , to 10 8 M from top to bottom, respectively. glyph[circledot] Right panel: Distance of BDDGs to the nearest clusters/groups. The distance of each BDDG to its clusters/groups are normalized by the corresponding virial radius of each cluster/group. Solid horizontal line indicate the 3 × R vir value, a limit where the effect of cluster/group could persists. We have 14 out of the 19 BDDGs beyond 3 × R vir of their surrounding clusters/groups. Among the remaining 5 BDDGs, 4 are close to the Virgo cluster.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00046v2-figure-03af968e-b1d5-4cf9-87b1-1f46f4687c78.jpg)

- [5] van Dokkum P. et al. A galaxy lacking dark matter. Natur. 555 , 629-632 (2018).
- [6] van Dokkum P. et al. A second galaxy missing dark matter in the NGC 1052 group. Astrophys. J. 874 , L5-L13 (2019).
- [7] Oman, A.-K. et al. Missing dark matter in dwarf galaxies? Mon. Not. R. Astron. Soc. 460 , 3610 (2016)
- [8] Planck Collaboration et al. Planck 2018 results. VI. Cosmological parameters. Preprint at http://arxiv.org/abs/1807.06209 (2018).
- [9] Simon, J. D. &amp; Geha, M. The kinematics of the ultra-faint Milky Way satellites: Solving the missing satellite problem. Astrophys. J. 670 , 313-331 (2007).
- [10] Martin, N. F. et al. A Keck/DEIMOS spectroscopic survey of faint Galactic satellites: searching for the least massive dwarf galaxies. Mon. Not. R. Astron. Soc. 380 , 281-300 (2007).
- [11] Sawala, T. et al. Bent by baryons: the low-mass galaxy-halo relation. Mon. Not. R. Astron. Soc. 448 , 2941-2947 (2015).
- [12] Oh, S.-H. et al. High-resolution mass models of dwarf galaxies from LITTLE THINGS. Astronomical. J. 149 , 180 (2015)
- [13] Bournaud, F. et al. Missing mass in collisional debris from galaxies. Science. 316 , 1166 (2007)
- [14] Fritz, T. K. et al. Gaia DR2 proper motions of dwarf galaxies within 420 kpc: orbits, Milky Way mass, tidal influences, planar alignments, and group infall. Preprint at http://arxiv.org/abs/1805.00908 (2018).
- [15] Haynes, M. P. et al. The Arecibo Legacy Fast ALFA survey: the α .40 HI source catalog, its characteristics and their Impact on the derivation of the HI mass function. Astronomical.J. 142 , 170 (2011).

- [16] Abazajian, K. N. et al. The seventh data release of the Sloan Digital Sky Survey. Astrophys. J. Suppl. Ser. 182 , 543-558 (2009).
- [17] Martin, D. C. et al. The Galaxy Evolution Explorer: A space ultraviolet survey mission. Astrophys. J. 619 , L1-L6 (2005).
- [18] Lisker, T.; Grebel, E. K.; Binggeli, B. &amp; Katharina, G. Virgo cluster early-type dwarf galaxies with the Sloan Digital Sky Survey. III. subpopulations: distributions, shapes, origins. Astrophys. J. 660 , 1186-1197 (2007).
- [19] Zheng, Z. et al. The structure and stellar content of the outer disks of galaxies: A new view from the Pan-STARRS1 medium deep survey. Astrophys. J. 800 , 120 (2015).
- [20] Bell E. F., McIntosh D. H., Katz N. &amp; Weinberg M. D. The optical and near-infrared properties of galaxies. I. luminosity and stellar mass functions. Astrophys. J. Suppl. Ser. 149 , 289-312 (2003).
- [21] Wang J. et al. New lessons from the HI size-mass relation of galaxies. Mon. Not. R. Astron. Soc. 460 , 2143-2151 (2016).
- [22] Navarro J. F.; Frenk C. S. &amp; White S. D. M. A universal density profile from hierarchical clustering. Astrophys. J. 490 , 493-508 (1997)
- [23] Burkert, A. The structure of dark matter halos in Dwarf Galaxies. Astrophys. J. 447 , L25-L28 (1995).
- [24] Blanton, M. R.; Geha, M. &amp; West, A. A. Testing cold dark matter with the low-mass TullyFisher relation. Astrophys. J. 682 , 861-873 (2008).
- [25] Springob, C. M.; Haynes, M. P.; Giovanelli, R. &amp; Kent, B. R. A digital archive of H I 21 centimeter line spectra of optically targeted galaxies. Astrophys. J. Suppl. Ser. 160 , 149-162 (2005).

- [26] Jing. Y. J. et al. The dark matter deficit galaxies in hydrodynamical simulations. Preprint at https://arxiv.org/abs/1811.09070 (2018).
- [27] Haslbauer, M. et al. Galaxies lacking dark matter in the Illustris simulation. Astron. Astrophys. 626 , A47 (2019).
- [28] Hunter, D. A. et al. Little things. Astronomical. J. 144 , 134 (2012)
- [29] Vogelsberger, M. et al. Introducing the Illustris Project: simulating the coevolution of dark and visible matter in the Universe Mon. Not. R. Astron. Soc. 444 , 1518-1547 (2014).
- [30] Diaz-Garcia, S.; Salo, H.; Laurikainen, E. &amp; Herrera-Endoqui, M. Astronomy &amp; Astrophysics 587 , A160-A220 (2016)

## METHODS

ALFALFA. The Arecibo Legacy Fast ALFA (Arecibo L-band Feed Array) survey(ALFALFA) 31-33 is a blind extragalactic HI survey, covering a wide area of 7000 deg 2 of the sky at high Galactic latitudes. The beam size is about 3.5 arcmin. It has HI line spectra of more than 30,000 extragalactic sources with radial velocities &lt; 18000 km s -1 . The spectral resolution 34 is ∼ 10 km s -1 , enabling reliable measurements of ∼ 20 km s -1 , the typical velocity of dwarf galaxies. Here we use α .40 15 , a catalog which covers 40 % of the full ALFALFA survey area, including the northern hemisphere spring region: 07 h 30 m &lt; ra &lt; 16 30 h m , +04 ◦ &lt; dec &lt; +16 , and +24 ◦ ◦ &lt; dec &lt; +28 , ◦ and the northern hemisphere autumn region: 22 00 h m &lt; ra &lt; 03 00 h m , +14 ◦ &lt; dec &lt; +16 , and ◦ +24 ◦ &lt; dec &lt; +32 . There are 15855 HI detected sources, among which 15041 are extragalactic ◦ objects and the rest 814 are more likely Galactic high velocity cloud (HVC).

The catalogue provides w 50 , 50% peak width of HI spectrum, and the corresponding errors. However, we find that in most cases, the edge of neither side of the spectrum has a cut-off feature. Instead, the flux at each side drops gradually. Calculations based on w 50 could lead to underestimates of the dynamical masses enclosed in r HI . w 20 , 20% peak width of the HI spectrum, is a better indicator of the gas velocity. In practice, we treat systems with one peak and two peaks in their HI spectra differently. For those with one peak, wavelengths at which the fluxes are 20% of the peak values are identified to determine the w 20 raw , , while for those with two peaks, the 20%-fluxes are calculated for both peak values and the corresponding characteristic wavelengths are identified at each rising side of the spectrum to determine the w 20 raw , . We perform Gaussian/Lorentzian fitting at each rising side. The results are consistent with those using polynomial fits. For each source, we generate 10000 realizations using the rms given by the ALFALFA catalogue, and adopt the corresponding 1σ range as the uncertainty in w 20 ( w 20 err , ). Then, we subtract the instrumental broadening as following 35

$$w _ { 2 0 } = & \sqrt { w _ { 2 0, r a w } ^ { 2 } - ( \delta s ) ^ { 2 } }$$

where δs is the final velocity resolution after smoothing and we adopted the value of 10 km/s.

SDSS DR7. The optical galaxy catalogue is based on SDSS DR7 16 , but with an independent set of improvement in sky-background subtraction. SDSS DR7 contains images and spectra for millions of galaxies, centered at z ∼ 0.1. The photometric data completes at 24 r-band apparent magnitude and the extinction-corrected spectroscopy data are complete at 17.72.

It is found that the SDSS photometric pipeline systematically overestimated the sky backgrounds of galaxies with extended low surface brightness by ∼ 0.5 mag 18,36,37 . The surface brightness of dwarf galaxies is usually very low and their luminosities could thus be very sensitive to the sky subtraction. Here we select a sub-sample of low surface brightness galaxy (central surface brightness in the B band µ 0 (B) &gt; 22 5 . mag arcsec -2 ) 38 and use the g-band and r-band SDSS mosaic images from the Twelfth Data Releases of the Sloan Digital Sky Survey(SDSS DR12) 39 to re-calculate their photometry 19 . The method is summarized as following: we first use SEXtractor 40 to extract and mask all sources above 2 σ local background variation using the r-band/g-band image and then unmask the target galaxy by visual inspection. We grow the masks using the IDL function DILATE with a square shape operator, whose size is proportional to the size of the mask. Stars on top of our target galaxies are removed and filled with the average value of pixels around the star. We fit each galaxy's ellipticity profile using the ELLIPSE task in IRAF and use the ellipticity of the largest fitted ellipse to calculate the r-band axis ratio b/a. The Petrosian flux is then calculated by summing up all fluxes within 2 times the elliptical Petrosian radius 41,42 . For those with smaller surface brightness, we simply adopt the Petrosian magnitudes, and use exponential axes ratio expAB r (or de Vaucouleurs axes ratio deVAB ) if the likelihood fraction of de Vaucouleurs r model fracDeV r &lt; 0.5 (fracDeV r &gt; 0.5) as the b/a value from the released SDSS DR7 catalogue.

There are 12423 galaxies with optical counterparts (OCs) in the HI and optical cross-matched catalogue 15 . Here we focus on dwarf galaxies with absolute magnitude M r &gt; -18 and high signalto-noise ratio in HI spectra: SNR HI &gt; 10 . To avoid significant bias caused by inclination, we further require the r-band b/a to be within 0.3 and 0.6, corresponding to the inclination angle range of [54.7, 76.8] degree. The final parent sample consists of 324 galaxies.

LITTLE THINGS. 'Local Irregulars That Trace Luminosity Extremes, The HI Nearby Galaxy Survey'(LITTLE THINGS) 28 is a high resolution, very large array HI survey with angular resolution of ∼ 6 arcsec and velocity resolution of &lt; 2.6 km s -1 , aiming to obtain deep HI-line maps of dwarf galaxies. The high resolution enable a reliable measurement of the rotation curve 12 . In practice, Oh et al. 12 fitted a 2D tilted-ring model to the 2D bulk velocity fields extracted from the HI data. This provides geometric and kinematic parameters for each galaxy, including central position, position angle, inclination, systemic velocity and rotation velocity. The asymmetric drift is then made to improve the rotation curve measurement by taking into account the dynamical support due to random motions. In combination with Spitzer 3.6 µ mand optical U, B, and V images, the rotation curve is then used to analyze the dynamical contributions by baryons and dark matter. They found the inner slope of the dark matter density profiles deviate significantly from the cusplike mass distribution predicted by dark matter only simulations, but in line with those predicted by Λ CDM + hydro-dynamical simulations, for which the baryon feedback process could change the inner profile.

Mass algorithm. Four kinds of masses are used in this work, including stellar mass, HI mass, dynamical mass, and dark halo virial mass. Here we describe briefly how each of the mass is derived.

Gas mass. Gas in galaxies mainly consists of HI and Helium. The HI mass is calculated following Roberts (1962) 43

$$M _ { \text{HI} } = 2. 3 6 \times 1 0 ^ { 5 } d ^ { 2 } \int F \text{dv} [ M _ { \odot } ]$$

where F is the HI flux in Jy km s -1 , and v is the receding velocity. d is the distance 15 of the HI detected galaxies estimated using d = cz/H 0 for galaxies with cz &gt; 6000 km/s, where H 0 is the Hubble constant, 70km s -1 , and c is the speed of light. For galaxies with cz &lt; 6000 km/s, the distance is derived adopting the local universe peculiar velocity model 44 , using a combination of primary distances form the literature and secondary distances from the Tully-Fisher relation. The total gas mass is derived by assuming the same Helium-HI ratio as that derived from the Big Bang nucleosynthesis: M gas = 1 33 . × M HI (Y BBN P = 0 247 . ) 8 , where Y BBN P is the primordial He abundance from Planck in Λ CDM.

Stellar mass. For a given galaxy, we calculate the absolute magnitude in g and r band, using the photometry as described above and the distance from α. 40. Following Bell's algorithm, we have

$$\log _ { 1 0 } ( M / L ) = a _ { \lambda } + b _ { \lambda } \times ( M _ { g } - M _ { r } ).$$

Here we adopt a Kroupa initial mass function(IMF) 45 with a r = -0 306 . -0 15 . and b r = 1 097 . . The final error bars accounts for both observational errors in photometry and distance, and the method error of 0.1 dex as suggested by Bell 20 .

Dynamical mass. With the velocity and radius of the HI component, we calculate the dynamical mass enclosed in r HI by assuming a spherical symmetry: M dyn ( &lt; r HI ) =V 2 HI HI r /G, where G is the gravitational constant. Previous work 21 found that there is a tight relation between the HI mass and the HI radius r HI , defined at a surface density of 1M glyph[circledot] pc -2 .

$$\log _ { 1 0 } \mathfrak { r } _ { \text{HI} } = ( 0. 5 0 6 \pm 0. 0 0 3 ) \log _ { 1 0 } M _ { \text{HI} } - ( 3. 2 9 3 \pm 0. 0 0 9 ) - \log _ { 1 0 } 2$$

This relation has been proven successful in describing the r HI and M HI relation for various galaxies types 21 . In particular, we apply this r HI vs. M HI relation to the dwarf galaxies from LITTLE THINGS (Left panel of Extended Data Figure 1) and find they agree very well with the direct measurement using the resolved HI images.

The velocity V HI is estimated using w 20 : V HI = w 20 /2/sin(i), where i is the inclination angel.

$$\sin ( i ) = & \sqrt { \frac { 1 - ( \frac { b } { a } ) ^ { 2 } } { 1 - q _ { 0 } ^ { 2 } } }$$

The parameter q 0 is the intrinsic axial ratio of a galaxy seen edge-on. The value of q 0 could depend on galaxy morphology 46,47 . Here we use 0.2 as usually adopted in the literature. w 20 is a good indicator of the circular velocity given sharp edges at both sides of the given spectrum. However, for dwarf galaxies, velocity dispersion could contribute a non-negligible fraction to the total kinematic energy of the systems. This is reflected by the shape of the spectrum that at both sides the spectrum decline gradually, lacking the feature of a sharp cut-off. w 20 could thus carry both information of the circular velocity and of the velocity dispersion. We test the choice of using w 20 as an indicator of the asymmetric drift-corrected circular velocity using dwarf galaxies in the LITTLE THINGS (Middle panel of Extended Data Figure 1). It shows the w 20 is a good indicator of the corrected circular velocity, while w 50 usually under-estimates the corrected circular velocity. We then test the full method by performing the same analysis as for our BDDGs on the LITTLE THINGS galaxies and find the derived M dyn are consistent with the estimates using the tilted-ring model (Right panel in Extended Data Figure 1).

Note for we do not have HI images, we use the inclination angels estimated in the optical r-band. There could be mis-alignment between r-band image and HI velocity field, leading to biased estimate of the circular velocity. We restrict our parent sample to have r-band b/a between [0.3, 0.6] in order to avoid severe deviation of the derived velocities from the true ones. Using dwarf galaxies in the LITTLE THINGS, we show in the left panel of Extended Data Figure 2 the distribution of the angle between r-band image and the HI velocity filed. It shows that in most cases the mis-alignment is below 20 degree, thought there is a tail towards large angle. We apply the same distribution of the mis-alignment to our parent sample and generate 10000 realization for each parent galaxy. The resulting M dyn /M bary distribution is shown in the right panel of Extended Data Figure 2. It shows there is still an apparent excess of BDDGs at the low M dyn /M bary end. The existence of BDDGs is also supported when applying the same distribution of the mis-alignment as those in the simulations 48 .

If the BDDGs were very compact dwarf galaxies, it would be possible to have baryon mass dominated in the compact region. We compare the r-band r e and r HI of the BDDGs and those of the parent dwarf galaxies and find that the BDDGs are slightly more extended than the parent dwarf galaxies, both in the optical and in the HI contents (Extended Data Figure 3). In addition, r HI is larger than r e for all the BDDGs. Both the baryonic masses and dynamical masses are measured at radii within which it is expected to be dominated by dark matter if the BDDGs were typical dwarf galaxies.

Halo mass. Assuming a typical NFW profile 22 , r HI and M dyn ( &lt; r HI ) are related as

$$M _ { d y n } ( < r _ { \text{HI} } ) & = \int _ { 0 } ^ { r _ { \text{HI} } } 4 \pi r ^ { 2 } \rho ( r ) d r \\ & = 4 \pi \rho _ { 0 } R _ { s } ^ { 3 } \left [ \ln \left ( \frac { R _ { s } + r _ { \text{HI} } } { R _ { s } } \right ) - \frac { r _ { \text{HI} } } { R _ { s } + r _ { \text{HI} } } \right ],$$

where R s is a free parameter. Given r HI and M dyn ( &lt; r HI ) , for each galaxy, we obtain the corre-

## sponding R . The virial mass is expressed as s

$$M _ { \text{vir} } = \int _ { 0 } ^ { \text{R} _ { \text{vir} } } 4 \pi \text{r} ^ { 2 } \rho ( \text{r} ) \text{dr} = 4 \pi \rho _ { 0 } \text{R} _ { \text{s} } ^ { 3 } \left [ \ln \left ( 1 + c \right ) - \frac { c } { 1 + c } \right ],$$

$$\rho _ { 0 } = \frac { 2 0 0 \rho _ { c r } } { 3 } \frac { c ^ { 3 } } { \ln ( 1 + c ) - \frac { c } { 1 + c } },$$

where ρ 0 is a characteristic density, ρ cr is the critical density, and c is the concentration parameter which relates to R s and R vir via c = R vir R s . The concentration parameter follows a tight halo massconcentration relation 49 .

$$\log _ { 1 0 } c = 0. 9 0 5 - 0. 1 0 1 ( \log _ { 1 0 } M _ { \text{vir} } / 1 0 ^ { 1 2 } h ^ { - 1 } M _ { \odot } ).$$

Here h = H /(100 km s 0 -1 Mpc -1 )=0.7. Inserting Eq. 9 to Eq. 7, one could obtain the M vir .

Assuming a Burkert profile 23 , r HI and M dyn ( &lt; r HI ) are related as

$$M _ { \text{dyn} } ( < r _ { \text{HI} } ) & = \int _ { 0 } ^ { r _ { \text{HI} } } 4 \pi r ^ { 2 } \rho ( r ) \text{d} r \\ & = 2 \pi \rho _ { 0 } r _ { 0 } ^ { 3 } \left [ \ln \left ( 1 + \frac { r _ { \text{HI} } } { r _ { 0 } } \right ) + \frac { 1 } { 2 } \ln \left ( 1 + \frac { r _ { \text{HI} } ^ { 2 } } { r _ { 0 } ^ { 2 } } \right ) - a r c t a n ( \frac { r _ { \text{HI} } } { r _ { 0 } } ) \right ],$$

where r 0 and ρ 0 are free parameters. With the relation between r 0 and M vir 50

$$\log _ { 1 0 } ( r _ { 0 } / \text{kpc} ) = 0. 6 6 - 0. 5 8 ( \log _ { 1 0 } M _ { \text{vir} } / 1 0 ^ { 1 1 } M _ { \odot } ), & & ( 1 1 )$$

one could derive M vir for any given combination of r HI and M dyn ( &lt; r HI ) .

Galaxy clusters. The galaxies groups and clusters catalogues are created based on the SDSS DR12 39 and the 2MASS Redshift Survey (2MRS) 51 , using a friends-of-friends group finder algorithm 52 . Various observational biases have been taken into account, including Malmquist bias, the 'Fingers of God' and etc. The groups/clusters mass depends on the total luminosity, the luminosity distance, the velocity dispersion, the group radius and the number of detected group members. This method has been carefully calibrated and optimized on wide-angle mock catalogues from the Millennium simulation 53,54 . For any group/cluster included in these catalogues, where available, group parameters in the published literature are adopted primarily.

The properties of the 19 BDDGs are summarized in Extend Data Table 1.

- [31] Giovanelli, R. et al. The Arecibo Legacy Fast ALFA Survey. I. science goals, survey design, and strategy. Astronomical.J. 130 , 2598 (2005).
- [32] Giovanelli, R. ALFALFA: An exploration of the z=0 HI universe. Il. Nuovo. Cimento. B 122 , 1097-1108 (2007).
- [33] Haynes, M. P. ALFALFA: The search for (almost) dark galaxies and their space distribution. Il. Nuovo. Cimento. B 122 , 1109-1114 (2007).
- [34] Giovanelli, R. et al. The Arecibo Legacy Fast ALFA Survey. II. results of precursor observations. Astronomical.J. 130 , 2613 (2005).
- [35] Catinella, B. et al. The GALEX Arecibo SDSS Survey. VI. Second data release and updated gas fraction scaling relations. A. &amp;A. 544 , A65 (2012).
- [36] He, Y. Q. et al. Photometric properties and luminosity function of nearby massive early-type galaxies. Astrophys. J. 773 , 37 (2013).
- [37] Liu, F. S.; Xia, X. Y.; Mao, S.; Wu, H. &amp; Deng, Z. G. Photometric properties and scaling relations of early-type brightest cluster galaxies. Mon. Not. R. Astron. Soc. 385 , 23-29 (2008).
- [38] Du, Wei et al. Low surface brightness galaxies selected from the 40% sky area of the ALFALFA H I survey. I. sample and statistical properties. Astronomical.J. 149 , 199 (2015).
- [39] Alam, S. et al. The eleventh and twelfth data releases of the Sloan Digital Sky Survey: final data from SDSS-III. Astrophys. J. Suppl. Ser. 219 , 12 (2015).
- [40] Bertin, E. &amp; Arnouts. SExtractor: Software for source extraction. A.&amp;A. S. 117 , 393-404 (1996).
- [41] Blanton, M. R. et al. The luminosity function of galaxies in SDSS commissioning data. Astronomical. J. 121 , 2358-2380 (2001).
- [42] Yasuda, N. et al. Galaxy number counts from the Sloan Digital Sky Survey commissioning data. Astronomical. J. 122 , 1104-1124 (2001).
- [43] Roberts, M. S. The neutral hydrogen content of late-type spiral galaxies. Astronomical. J. 67 , 437-446 (1962).
- [44] Masters, K. L. Galaxy flows in and around the Local Supercluster. Ph.D dissertation 0058 , 0606 (2005) AAT 3192074.
- [45] Kroupa, P. The initial mass function of stars: Evidence for uniformity in variable systems. Science. 295 , 82-91 (2002)
- [46] Giovanelli, R. et al. The I band Tully-Fisher relation for cluster galaxies: data presentation. Astronomical. J. 113 , 22-52 (1997)
- [47] Tully, B. R. et al. The extragalactic distance database. Astronomical. J. 138 , 323-331 (2009)
- [48] Starkenburg S. et al. On the origin of star-gas counterrotation in low-mass galaxies. Astrophys. J. 878 , 143-155 (2019).
- [49] Dutton A. A. &amp; Macci ` o A. V. Cold dark matter haloes in the Planck era: evolution of structural parameters for Einasto and NFW profiles. Mon. Not. R. Astron. Soc. 441 , 3359-3374 (2014).
- [50] Salucci P. et al. The universal rotation curve of spiral galaxies II. The dark matter distribution out to the virial radius. Mon. Not. R. Astron. Soc. 378 , 41-44 (2007).
- [51] Huchra, J. P. et al. The 2MASS Redshift Survey - description and data release. Astrophys. J. Suppl. Ser. 199 , 26 (2012)
- [52] Saulder, C.; van Kampen, E.; Chilingarian, I.; Mieske, S. &amp; Zeilinger, W. W. The matter distribution in the local Universe as derived from galaxy groups in SDSS DR12 and 2MRS. A.&amp;A. 596 , A14 (2016).
- [53] Guo, Qi et al. From dwarf spheroidals to cD galaxies: simulating the galaxy population in a ΛCDM cosmology. Mon. Not. R. Astron. Soc. 413 , 101-131 (2011).
- [54] Guo, Qi et al. Galaxy formation in WMAP1 and WMAP7 cosmologies. Mon. Not. R. Astron. Soc. 428 , 1351-1365 (2013).
- [55] Rines, K.; Geller, M.; Kurtz, M. &amp; Diaferio, A. CAIRNS: The Cluster and Infall Region Nearby Survey. I. Redshifts and Mass Profiles. Astronomical. J. 126 , 2152-2170 (2003).

## Acknowledgements

This work is supported by the National Key R&amp;D Program of China (Nos 2018YFA0404503, 2016YFA0400703 and 2016YFA0400702), the National Natural Science Foundation of China (No.11573033, 11622325, 11133003, 11425312, 11773001, 11721303, 11733006, 11703036) and the National Natural Science Foundation of

China (grant numbers 11573033, 11622325, 11133003, 11425312, 11773001, 11721303, 11733006, 11703036 and U1931110). Q.G. and L.G. acknowledge support from the Royal Society Newton Advanced Fellowships. Z.Z acknowledges the support by the Open Project Program of the Key Laboratory of FAST, Chinese Academy of Sciences.

## Author contributions

Q.G. led and played a part in all aspects of the analysis, and wrote the manuscript. H.H compiled the data and carried out most of the data reduction and analysis. Z.Z. measured the optical photometry of low-surface brightness galaxy. S.L carried out most of the analysis related to various mass determinations. All authors contributed to the analysis, and to the writing of the manuscript.

## Additional information

Supplementary information is available in the on-line version of the paper.

## Competing financial interests

The authors declare that they have no competing financial interests. Correspondence and requests for materials should be addressed to Q.G. (guoqi@nao.cas.cn).

Data availability. The ALFALFA data are available in the The Arecibo Legacy Fast ALFA Survey (http://egg.astro.cornell.e and the SDSS data are available at the Sloan Digital Sky Survey (https://www.sdss.org/) . The other data that support the results of this study are available from the corresponding author upon reasonable request.

Code availability. We use standard data reduction tools in the IDL and Python environments, and the publicly available codes SExtractor (https://www.astromatic.net/software/sextractor) for this study.

Extended Data Table 1: BDDGs data . Column 1: ALFALFA galaxy ID; Cols. 2and 3: equatorial coordinates; Col. 4: distance judged by ALFALFA group; Col. 5 heliocentric velocity; Col. 6: velocity width 20% of the peak flux; Col. 7: error of w 20 ; Cols. 8 and 9: g-band, r-band magnitude; Col. 10: effective radius in r-band; Col. 11: neutral hydrogen mass; Col. 12: stellar mass; Col.13: dynamical mass within r HI ; Col. 14: virial mass; Col. 15: r-band minor-to-major axis ratio; Col. 15: distance to the nearest group over virial radius

| dist gr /R vir,gr                      | 2.707 9.224 0.334 0.617 1.885 4.761 5.408 9.818 7.306 0.38 6.4 8.924 4.701 20.122 6.399 5.913 3.271 4.191 4.051                                                                                                                                                      |
|----------------------------------------|----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| b/a ratio                              | 0.524 0.55 0.491 0.4 0.6 0.444 0.57 0.334 0.469 0.486 0.551 0.576 0.518 0.326 0.592 0.515 0.365 0.563 0.464                                                                                                                                                          |
| log M vir [M glyph[circledot] ]        | 10.231 9.876 10.599 10.653 9.515 9.712 9.567 9.984 10.149 9.363 9.96 10.098 10.387 9.593 8.904 9.489 10.003 9.708 9.759                                                                                                                                              |
| log M dy [M glyph[circledot] ]         | 9.444 9.592 9.061 8.981 9.046 9.092 9.08 9.706 9.8 8.864 9.547 9.467 10.124 9.204 8.503 9.042 9.778 9.264 9.445                                                                                                                                                      |
| HI log M stellar [M glyph[circledot] ] | 8.951 8.731 8.749 9.046 7.702 8.667 8.273 8.66 8.741 6.915 9.02 9.133 8.76 8.105 6.196 8.276 7.066 7.863 8.504                                                                                                                                                       |
| log M [M glyph[circledot] ]            | 8.58 9.44 7.82 7.82 8.64 9.21 8.64 9.45 9.39 8.46 9.47 9.65 9.78 8.88 8.34 8.67 9.61 8.84 9.19                                                                                                                                                                       |
| re [kpc]                               | -17.75 0.59 -17.89 4.93 0.84 1.54 0.77 5.19 1.64 2.86 3.39 0.74 4.61 0.72 2.88 2.17 0.22 0.39 1.29 1.3 2.13                                                                                                                                                          |
| M r [mag]                              | -17.2 -17.71 -15.21 -17.56 -16.31 -17.33 -17.62 -14.07 -17.75 -18.0 -17.95 -16.73 -11.83 -16.47 -14.99 -15.72 -17.01                                                                                                                                                 |
| M g [mag]                              | -17.34 -17.74 -16.77 -17.2 -15.01 -17.34 -16.0 -17.03 -17.36 -14.18 -17.27 -17.52 -17.79 -16.72 -11.77 -16.2 -15.29 -15.57 -16.74                                                                                                                                    |
| w20er [km/s]                           | 2.03 1.54 4.45 2.6 0.83 0.31 1.21 4.41 4.33 0.74 1.98 1.18 8.53 4.41 2.9 0.99 4.91 3.06 5.49                                                                                                                                                                         |
| w20 [km/s]                             | 80.36 56.63 82.37 79.03 46.12 39.08 49.27 72.5 78.35 45.38 52.82 42.47 87.79 56.83 29.59 48.36 70.87 54.51 58.68                                                                                                                                                     |
| distance V helio [km/s]                | 1156 3342 2737 849 694 1690 1595 7125 6145 1005 5249 5260 10034 3863 2860 1280 7608 4374 6840                                                                                                                                                                        |
| DEC [deg ( ◦ )] [Mpc]                  | 9.98695 20.4 24.47278 50.3 13.36306 16.7 12.35028 16.7 3.84222 16.6 7.8625 27.6 27.56833 25.1 13.97361 106.7 8.78806 92.8 13.78139 17 14.24972 77.4 4.90972 76.9 6.28583 145.4 32.18834 56.6 14.26417 41.2 26.66472 23.0 27.26361 112.7 24.96306 64.8 25.16611 101.9 |
| RA [deg ( ◦ )]                         | 171.47292 179.77542 189.76666 191.18918 192.44583 221.33916 144.44835 163.98709 171.80333 189.99374 213.8425 212.19499 233.61292 341.93207 351.52747 224.41417 190.81416 179.17708 186.755                                                                           |
|                                        | 6980 ∗ 7817 7920 7983 9500 191707 205215 213086 † 220901 241266 † 242440 258421 321435 331776 733302 749244 † 749445 † 749457                                                                                                                                        |
| AGCNr                                  | AGC 6438 AGC AGC AGC AGC AGC AGC AGC AGC AGC AGC AGC AGC AGC AGC AGC AGC AGC AGC                                                                                                                                                                                     |

The projected distance between the AGC 6980 and Abell 1367 is about 8 Mpc, well

.

55

, ra = 176.15083, dec=19.77194 )

cz = 6459

In the literature AGC 6980 was defined as Abell 1367's member(Leo cluster:

*

of the Abell 1367.

vir beyond 3 R

The minimum projected distance of the BDDG to its surroundings clusters with relatively small difference in redshifts (cz BDDG - cz cluster &lt; 10 × R cluster vir , ) is &lt; 3 × R vir .

†

Extended Data Figure 1: Comparison of the galaxy properties between different methods for dwarf galaxies in LITTLE THINGS. Left panel: Comparison of the r HI using our method with that obtained using tilted-ring method by LITTLE THINGS team. Middle panel: Comparison of the V HI using our method with that obtained using tilted-ring method by LITTLE THINGS team. Right panel: Comparison of the dynamical mass estimated using our method with those derived using the tilted-ring method. Black and grey dots are measurements based on w 20 and w 50 , respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00046v2-figure-04982e7c-556d-4c2c-a347-0e3f538fd51a.jpg)

Extended Data Figure 2: Distribution of the M dyn -to-M bary ratio with different misalignments between optical image and HI velocity field. Left panel: The distribution of the misalignment between optical image and HI velocity field for dwarf galaxies in LITTLE THINGS(red)/Illustris(blue). Middle/Right panel: Distribution of the M dyn -to-M bary ratio for our parent sample after applying the same distribution of misalignment between optical image and HI velocity field from the LITTLE THINGS/Illustris in the left panel. Grey histogram is based on 10000 random realizations and blue curve is the Gaussian fit. Samples classified as BDDGs are high-lighted with black histogram.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00046v2-figure-59d78e05-786f-458f-bdb4-ab5ac2a7ed75.jpg)

Extended Data Figure 3: Comparison of r HI and r e between the parent sample and the BDDGs. Left panel: r HI and r e relation of the parent sample(grey dots) and the BDDGs(black dots). Upper right panel: the distribution of r HI , grey histogram and blue histogram represent parent sample and BDDGs respectively. Bottom right panel: the distribution of r , grey histogram e and blue histogram represent parent sample and BDDGs respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00046v2-figure-1f87381e-f62e-44f5-a896-58a9633e2475.jpg)
Qi Guo Huijie Hu Zheng Zheng Shihong Liao Wei Du Shude Mao Linhua Jiang Jing Wang Yingjie Peng Liang Gao Jie Wang Hong Wu
Further evidence for a population of dark-matter-deficient dwarf galaxies
In the standard cosmological model, dark matter drives the structure
formation and constructs potential wells within which galaxies may form. The
baryon fraction in dark halos can reach the universal value (15.7%) in massive
clusters and decreases rapidly as the mass of the system decreases. The
formation of dwarf galaxies is sensitive both to baryonic processes and the
properties of dark matter owing to the shallow potential wells in which they
form. In dwarf galaxies in the Local Group, dark matter dominates the mass
content even within their optical-light half-radii (r_e ~ 1 kpc). However,
recently it has been argued that not all dwarf galaxies are dominated by dark
matter. Here we report 19 dwarf galaxies that could consist mainly of baryons
up to radii well beyond r_e, at which point they are expected to be dominated
by dark matter. Of these, 14 are isolated dwarf galaxies, free from the
influence of nearby bright galaxies and high dense environments. This result
provides observational evidence that could challenge the formation theory of
low-mass galaxies within the framework of standard cosmology. Further
observations, in particular deep imaging and spatially-resolved kinematics, are
needed to constrain the baryon fraction better in such galaxies.