A second galaxy missing dark matter in the NGC1052 group

TL;DR

This paper reports the discovery of a second galaxy in the NGC1052 group with little or no dark matter, similar to NGC1052-DF2, challenging existing galaxy formation models.

Contribution

It identifies a new galaxy with low dark matter content in the same group, providing evidence for a class of such galaxies and raising questions about their origins.

Findings

The galaxy has a velocity dispersion consistent with stars alone.

It has a similar distance and morphology to NGC1052-DF2.

The galaxy's velocity dispersion is lower than expected from standard dark matter models.

Abstract

The ultra-diffuse galaxy NGC1052-DF2 has a very low velocity dispersion, indicating that it has little or no dark matter. Here we report the discovery of a second galaxy in this class, residing in the same group. NGC1052-DF4 closely resembles NGC1052-DF2 in terms of its size, surface brightness, and morphology; has a similar distance of $D=19.9\pm 2.8$ Mpc; and also has a population of luminous globular clusters extending out to 7 kpc from the center of the galaxy. Accurate radial velocities of the diffuse galaxy light and seven of the globular clusters were obtained with the Low Resolution Imaging Spectrograph on the Keck I telescope. The velocity of the diffuse light is identical to the median velocity of the clusters, $v_{\rm sys}=\langle v_{\rm gc} \rangle=1445$ km/s, and close to the central velocity of the NGC1052 group. The rms spread of the observed velocities is very small at $\sigma_{\rm obs}=5.8$ km/s. Taking observational uncertainties into account we determine an intrinsic velocity dispersion of $\sigma_{\rm intr}=4.2^{+4.4}_{-2.2}$ km/s, consistent with the expected value from the stars alone ($\sigma_{\rm stars}\approx 7$ km/s) and lower than expected from a standard NFW halo ($\sigma_{\rm halo}\sim 30$ km/s). We conclude that NGC1052-DF2 is not an isolated case but that a class of such objects exists. The origin of these large, faint galaxies with an excess of luminous globular clusters and an apparent lack of dark matter is, at present, not understood.

Tables

Table 1. Table 1 NGC1052-DF4 Globular Clusters

Id	RA	Dec	$R$aaDistance from the center of the galaxy.	$M_{V,606}$	$V-I$	$v$	$r_h$	$ϵ$
	(J2000)	(J2000)	[kpc]bbFor an assumed distance of $D=20$ Mpc.	[mag]bbFor an assumed distance of $D=20$ Mpc.	[mag]cc$V_{606}-I_{814}$ from HST/ACS, in the AB system.	[km s-1]	[pc]
NGC1052-DF4	2${}_{}{}^{h}39_{}^m{15.11}^s$	$-8°6\prime 58\stackrel{\prime \prime }{.}6$	$\cdots$	$-15.0$	$0.32$	${1444.6}_{-7.7}^{+7.8}$	$1600$	0.11
GC-2726	2${}_{}{}^{h}39_{}^m{16.75}^s$	$-8°6\prime 16\stackrel{\prime \prime }{.}7$	4.69	$-9.2$	$0.38$	${1441.2}_{-4.8}^{+4.9}$	$4.9\pm 0.6$	$0.14\pm 0.04$
GC-2537	2${}_{}{}^{h}39_{}^m{12.53}^s$	$-8°6\prime 41\stackrel{\prime \prime }{.}4$	4.08	$-9.2$	$0.36$	${1451.0}_{-3.3}^{+3.6}$	$4.1\pm 0.5$	$0.08\pm 0.06$
GC-2239	2${}_{}{}^{h}39_{}^m{15.23}^s$	$-8°6\prime 53\stackrel{\prime \prime }{.}0$	0.57	$-8.6$	$0.34$	${1457.1}_{-5.5}^{+4.6}$	$5.4\pm 0.5$	$0.16\pm 0.04$
GC-1968	2${}_{}{}^{h}39_{}^m{15.25}^s$	$-8°6\prime 58\stackrel{\prime \prime }{.}8$	0.20	$-9.8$	$0.27$	${1445.4}_{-2.3}^{+2.6}$	$3.4\pm 0.3$	$0.36\pm 0.05$
GC-1790	2${}_{}{}^{h}39_{}^m{17.24}^s$	$-8°7\prime 06\stackrel{\prime \prime }{.}7$	3.17	$-9.0$	$0.31$	${1438.4}_{-4.6}^{+4.8}$	$3.2\pm 0.4$	$0.19\pm 0.08$
GC-1452	2${}_{}{}^{h}39_{}^m{18.23}^s$	$-8°7\prime 24\stackrel{\prime \prime }{.}1$	5.13	$-9.1$	$0.32$	${1445.5}_{-4.1}^{+4.0}$	$5.1\pm 0.6$	$0.04\pm 0.03$
GC-943	2${}_{}{}^{h}39_{}^m{16.98}^s$	$-8°8\prime 5\stackrel{\prime \prime }{.}3$	7.01	$-8.6$	$0.34$	${1445.1}_{-5.2}^{+5.0}$	$3.3\pm 1.1$	$0.41\pm 0.17$
BG-2844	2${}_{}{}^{h}39_{}^m{20.37}^s$	$-8°6\prime 18\stackrel{\prime \prime }{.}6$	$\cdots$	$\cdots$	$\cdots$	$z=0.2298$	$\cdots$	$\cdots$
BG-254	2${}_{}{}^{h}39_{}^m{11.55}^s$	$-8°9\prime 03\stackrel{\prime \prime }{.}0$	$\cdots$	$\cdots$	$\cdots$	$z=0.2557$	$\cdots$	$\cdots$

Equations

$$\mathcal{L}=\prod _{i=1}^{i\le 7}\frac{1}{\sqrt{2\pi }{\sigma }_{\mathrm{eff}}}\exp \left[-0.5{\left(\frac{v_i-\mu }{{\sigma }_{\mathrm{eff}}}\right)}^2\right],$$

$$\delta (S(\hat{D}),S(D))<e$$

$$M_{\mathrm{TME}}=\frac{C}{G}⟨(\Delta v^{\prime }{)}^2{r}^{\alpha }⟩r_{\mathrm{out}}^{1-\alpha }.$$

$$C=\frac{4\Gamma \left(\frac{\alpha }{2}+\frac{5}{2}\right)}{\sqrt{\pi }\Gamma \left(\frac{\alpha }{2}+1\right)}\frac{\alpha +\gamma -2\beta }{\alpha +3-\beta (\alpha +2)},$$

Full text

††slugcomment: Accepted for publication in ApJ Letters

1. Introduction

Sensitive surveys using state of the art telescopes have identified large numbers of intrinsically-large galaxies with very low surface brightness (e.g., van Dokkum et al. 2015; Koda et al. 2015; van der Burg, Muzzin, & Hoekstra 2016). These “ultra diffuse galaxies” (UDGs), with sizes $R_{\rm e}>1.5$ kpc and central surface brightness $\mu_{g}>24$ mag arcsec*-2*, have been found in many different environments (including the Local Group; Martin et al. 2016; Torrealba et al. 2019), and have a wide range of properties (see, e.g., Merritt et al. 2016).

One of the most intriguing UDGs that have been studied so far is NGC1052-DF2 in the NGC 1052 group. Using a combination of Hubble Space Telescope HST Advanced Camera for Surveys (ACS) imaging and Keck spectroscopy, we determined that this galaxy has an unusual population of luminous globular cluster-like objects (van Dokkum et al. 2018c). Furthermore, from the radial velocities of ten of these globular clusters we determined that the galaxy appears to have little or no dark matter ($\lesssim 10^{8}$ M*⊙; van Dokkum et al. 2018b; Wasserman et al. 2018). Both aspects are surprising: the globular cluster luminosity function was thought to be universal (Rejkuba 2012), and a galaxy with a stellar mass of $\sim 2\times 10^{8}$ M⊙* should have a dark matter mass of $\sim 6\times 10^{10}$ M*⊙* (Behroozi, Wechsler, & Conroy 2013b).

Although these unexpected results were initially greeted with some skepticism (e.g., Martin et al. 2018; Laporte, Agnello, & Navarro 2018; Nusser 2018; Hayashi & Inoue 2018; Ogiya 2018; Trujillo et al. 2019), recent studies have corroborated the unusual nature of NGC1052-DF2: the distance to the galaxy was placed on surer footing ($D=19-20$ Mpc; van Dokkum et al. 2018a; Blakeslee & Cantiello 2018) and, crucially, the low mass of NGC1052-DF2 has been confirmed by measurements of the stellar velocity dispersion (Emsellem et al. 2019; Danieli et al. 2019). Specifically, Danieli et al. (2019) find $\sigma=8.5^{+2.3}_{-3.1}$ km s*-1* from a high resolution spectrum taken with the Keck Cosmic Web Imager.

At this point, the central unanswered question is whether NGC1052-DF2 is an isolated case or representative of a population of similar galaxies. This is important for judging the likelihood of interpretations that require unusual orbits or viewing angles (see, e.g., Ogiya 2018) and, most importantly, for judging the relevance of NGC1052-DF2 for our ideas about galaxy formation and the relation between dark matter and normal matter. With the important exception of tidal dwarfs (Bournaud et al. 2007; Gentile et al. 2007; Lelli et al. 2015), it is often thought that a gravitationally-dominant dark matter halo is the sine qua non for the formation of a galaxy (White & Rees 1978). If galaxies such as NGC1052-DF2 are fairly common we may have to revise our concept of what a galaxy is, and come up with alternative pathways for creating such large and relatively massive stellar systems.

Here we report the discovery of a galaxy that shares essentially all of NGC1052-DF2’s unusual properties, to a remarkable degree. It is in the same group, has a similar size, luminosity, and color, the same morphology, the same population of luminous globular clusters, and the same extremely low velocity dispersion.

2. NGC1052-DF4

NGC1052-DF4 is a low surface brightness galaxy in the field of the elliptical galaxy NGC 1052. It is part of a sample of 23 objects that we identified in images taken with the Dragonfly Telephoto Array (Abraham & van Dokkum 2014) and followed up with the ACS on HST (see Cohen et al. 2018). Unlike NGC1052-DF2, which had been described earlier by Fosbury et al. (1978) and Karachentsev et al. (2000), NGC1052-DF4 was discovered with Dragonfly.111This is a somewhat academic point as the galaxy is not particularly faint; it is clearly visible in Plate 1 of Fosbury et al. (1978) and in many other imaging datasets. The HST image is shown in Fig. 1.

Basic parameters, as derived from the HST imaging, are given in Cohen et al. (2018). Its surface brightness fluctuation distance is $D_{\rm sbf}=19.9\pm 2.8$ Mpc, and as highlighted in Fig. 5 of van Dokkum et al. (2018a) its color-magnitude diagram is very similar to that of NGC1052-DF2. We infer that the galaxy is part of the NGC 1052 group, and use $D=20$ Mpc for distance-dependent quantities. The galaxy is well-fitted by a 2D Sersic (1968) profile with a Sérsic index of $n=0.79$, an axis ratio of $b/a=0.89$, central surface brightness $\mu(V_{606},0)=23.7$, and major axis half-light radius $R_{\rm e}=1.6$ kpc. These properties place the galaxy just inside the UDG selection box (see Fig. 10 of Cohen et al. 2018). The total absolute magnitude is $M_{V,606}=-15.0\pm 0.1$, corresponding to $L_{V,606}=(7.7\pm 0.8)\times 10^{7}$ L*⊙. Assuming $(M/L)_{{\rm stars},V}=(2.0\pm 0.5)\,{\rm M}_{\odot}/{\rm L}_{\odot}$ (see van Dokkum et al. 2018b) the implied total stellar mass is $M_{\rm stars}=(1.5\pm 0.4)\times 10^{8}$ M⊙*.

At first glance NGC1052-DF4 seemed to lack the spectacular population of bright globular cluster-like objects that initially drew our attention to NGC1052-DF2 (see van Dokkum et al. 2018b, 2018c). However, careful inspection of the HST imaging data shows that the NGC1052-DF4 field actually does contain objects with similar sizes, colors, and apparent magnitudes as the globular clusters in NGC1052-DF2 – but that they are even more spread out relative to the body of the galaxy. This unexpected finding inspired us to obtain spectra of these candidate globular clusters, to determine whether they are actually in a single structure and, if so, to use them to constrain the mass of NGC1052-DF4.

3. Spectroscopy

We observed compact objects in the NGC1052-DF4 field with the dual arm Low Resolution Imaging Spectrograph (LRIS) on the Keck I telescope. The sample selection is modeled upon the properties of the confirmed clusters in NGC1052-DF2. We used SExtractor (Bertin & Arnouts 1996) to measure total magnitudes, colors, and FWHM sizes of objects in the HST images, following the methodology outlined in van Dokkum et al. (2018c). Priority was given to objects with $I_{814}<23$, $0.20<V_{606}-I_{814}<0.43$, and $0\farcs 12<{\rm FWHM}<0\farcs 30$, as all 11 clusters from van Dokkum et al. (2018c) satisfy these criteria. In NGC1052-DF4, seven objects fall within these limits. They have a mean total magnitude of $\langle I_{814}\rangle=22.10$ with a $1\sigma$ rms spread of 0.39 mag. Objects just outside these selection limits were given lower priority.222There are three objects with $23.0<I_{814}<23.5$ that satisfy the color and size criteria.

All seven high priority objects could be fitted in a single multi-slit mask, along with four lower priority targets. This mask was observed on November 6 2018 for a total of 19,800 s, split over eleven 1,800 s exposures. Conditions were excellent. On the blue side the low resolution 300 lines mm*-1* grism blazed at 5000 Å was used, and on the red side the 1200 lines mm*-1* grating blazed at 9000 Å. Here only the high resolution red-side observations are discussed. The slit width of the mask is $0\farcs 8$, providing a resolution of $\sigma_{\rm instr}\approx 25$ km s*-1* near the Calcium triplet.

The data reduction follows the same procedures as those outlined in van Dokkum et al. (2018b). Briefly, the slit edges and the sky lines are used to correct the spectra for distortions. The sky line modeling and subtraction is done using the method of Kelson (2003), which minimizes interpolation-related residuals. The sky lines are also used for wavelength calibration. Each individual 1,800 s exposure is analyzed independently to limit the effects of flexure on the distortion modeling and wavelength calibration. The exposures are combined using optimal weighting, and 1D spectra are extracted by weighting each line in the 2D spectrum by the S/N ratio.

An inspection of the spectra shows that all seven bright globular cluster candidates indeed have strong absorption features at the approximate redshift of the NGC 1052 group. Two of the four lower priority objects turn out to be compact background galaxies and the remaining two are too faint for a redshift measurement. We show the spectra of the seven high priority targets in Fig. 2, focusing on the regions near the redshifted $\lambda\lambda 8542.09,8662.14$ Å lines of the Ca triplet.333Note that the spectrum of GC-2726 is blended with a star forming background galaxy. The weaker $\lambda 8498.02$ Å line is mostly masked, as it coincides with a strong sky line at the redshift of NGC1052-DF4. Despite a shorter total integration time the spectra are of higher quality than those of most of the clusters in NGC1052-DF2 (van Dokkum et al. 2018b), due to better seeing and photometric conditions during the observations. The median S/N ratio near the Ca triplet is 17 Å*-1*.

We also extract the spectrum of the diffuse galaxy light, in the following way. The central two clusters (GC-1968 and GC-2239) were observed through a relatively long slit that covers most of the extent of NGC1052-DF4. We create an average sky spectrum from the other nine slits and subtract this from the background spectrum of the central slit, after scaling. The resulting spectrum is shown in the top panel of Fig. 2. The Ca triplet lines are clearly detected.

4. Kinematics

4.1. Radial Velocities

Radial velocities of the diffuse light and the seven clusters are determined by fitting the data with a synthetic 11 Gyr, [Fe/H]=$-1$ stellar population synthesis model (Conroy, Gunn, & White 2009; Choi et al. 2016), convolved to the instrumental resolution. The fit is performed in the redshifted Calcium triplet region 8520 Å$<\lambda<8740$ Å, using the emcee MCMC algorithm (Foreman-Mackey et al. 2013). The fit is regularized by dividing the data and template by a polynomial of order $100^{-1}\Delta\lambda$ [Å]. Two free parameters are fitted: the velocity and an additive continuum offset, to account for any template mismatch. With the exception of the blended spectrum of GC-2726 this parameter is typically close to zero. Errors are determined from simulations. In each simulation the residuals of the fits are randomly shuffled and the velocity fit is repeated. Residuals from sky lines are shuffled separately from the rest of the spectrum. The width of the distribution of the resulting velocities is taken as the uncertainty in the fit (see van Dokkum et al. 2018b). The results are not sensitive to the details of these procedures; owing to the high S/N ratio of the spectra the velocities and the associated uncertainties are very stable.444We note that we do not apply any corrections for possible slit alignment errors; inspection of alignment check images showed no evidence for offsets, but we cannot exclude systematic errors at the level of a few km/s.

The velocity of the diffuse light is $v_{\rm sys}=1444.6^{+7.8}_{-7.7}$ km s*-1*. This is very close to the average velocity of other galaxies in the NGC 1052 group. Including NGC1052-DF2 there are 22 galaxies in the NASA Extragalactic Database within a radius of two degrees centered on NGC 1052 in the velocity range $0<v<2500$ km s*-1*. All are in the range 1241 km s*-1* $<v<$1805 km s*-1*, with a biweight mean of $v_{\rm group}=1438\pm 25$ km s*-1* and width $\sigma_{\rm group}=128\pm 19$ km s*-1*. The velocities of the globular clusters are listed in Table 1 and displayed in the left panel of Fig. 3. The random uncertainties are small; the mean $\pm 1\sigma$ error is $\pm 4$ km s*-1*. The median (mean) velocity is $\langle v\rangle=1445$ km s*-1* (1446 km s*-1*), identical to the systemic velocity of the galaxy. This confirms that the clusters are associated with the galaxy, and suggests that the globular cluster system is at rest with respect to the stars.

4.2. Velocity Dispersion

The velocity range of the seven clusters is extremely small, echoing our earlier result for NGC1052-DF2. The observed rms, before correcting for observational uncertainties, is only 5.8 km s*-1* (compared to 10 km s*-1* for NGC1052-DF2). We use two methods to determine the intrinsic dispersion $\sigma_{\rm intr}$ and its associated uncertainties. Both methods use a generative model whose parameters are constrained by assessing the probability of measuring the observed velocity distribution. The model is a simple Gaussian with the center and width as free parameters. The classical method, used extensively for determining the kinematics of dwarf galaxies in the Local Group from the velocities of individual stars (e.g., Martin et al. 2007), is to construct the likelihood function:

[TABLE]

with $v_{i}$ the velocities of the individual tracers, $\mu$ the mean of the model, and $\sigma_{\rm eff}^{2}=\sigma_{\rm intr}^{2}+e_{i}^{2}$ with $\sigma_{\rm intr}$ the model dispersion and $e_{i}$ the uncertainty in velocity $v_{i}$ (calculated by averaging the positive and negative error bars). The likelihood, marginalized over $\mu$, is shown by the solid line in the right panel of Fig. 3. The likelihood analysis gives $\sigma_{\rm intr}=3.8^{+3.4}_{-2.6}$ km s*-1*, where the uncertainties contain 68 % of the probability. The 90 % (95 %) confidence level upper limit is 8.6 km s*-1* (10.4 km s*-1*).

The second method is approximate Bayesian computation (ABC; Tavare et al. 1997), which approximates the posterior distribution with simulations. For each set of model parameters ($\mu$, $\sigma_{\rm intr}$) a large number of simulated datasets $\hat{D}=(\hat{v}_{1},\hat{v}_{2},\ldots,\hat{v}_{7})$ are created by randomly drawing velocities from the model and linearly perturbing these values with errors that are themselves randomly drawn from Gaussians of width ($e_{1},e_{2},\ldots,e_{7}$). These $N$ datasets are compared to the actual data $D$ through a summary statistic $S$. Simulations that satisfy the criterion

[TABLE]

are retained, with $\delta$ the absolute distance between the summary statistics and $e$ a small positive number. For a given choice of ($\mu$, $\sigma_{\rm intr}$) the posterior is proportional to the number of simulations that are retained. ABC does not assume a functional form of the likelihood, and summary statistics can be chosen that are best suited to particular situations (see, e.g., van Dokkum et al. 2018b). Another advantage is that it does not suffer from the “small sample bias” discussed in Laporte et al. (2018). We calculate the ABC posterior using $N=10^{4}$, $e=0.1$, and the rms as the summary statistic (Fig. 3). ABC gives a similar result as the likelihood: $\sigma_{\rm intr}=4.2^{+4.4}_{-2.2}$ km s*-1*. The small difference may reflect the likelihood’s sensitivity to small sample bias (see Laporte et al. 2018).

4.3. Implied Mass

Quantitative constraints on the halo mass are highly uncertain with seven tracers and require extensive modeling (see, e.g., Laporte et al. 2018; Wasserman et al. 2018). Here we simply test the hypothesis that there is no dark matter halo and all the mass is in the form of stars. Following Beasley et al. (2016) and van Dokkum et al. (2018b) we estimate the mass within the outermost globular cluster using the tracer mass estimator (TME) method of Watkins, Evans, & An (2010):

[TABLE]

Here $\Delta v^{\prime}=f^{-1}(v-1445.7)$ are the velocities of the individual tracers, with $f=\sigma_{\rm obs}/\sigma_{\rm intr}=1.4_{-0.7}^{+1.5}$. The parameter $\alpha$ is the slope of the potential and $C$ is a constant given by

[TABLE]

with $\gamma$ the power-law slope of the 3D density profile of the tracers and $\beta=1-\sigma_{\rm t}^{2}/\sigma^{2}_{\rm r}$ the Binney anisotropy parameter. For simplicity we assume that $\beta=0$ and that the globular clusters trace the potential, so that $\gamma=\alpha+2$. If all the mass is in stars the potential is similar to that of a point mass for most of the globular clusters; hence we use $\alpha=1$ and $\gamma=3$.

With these assumptions we find an enclosed mass within 7 kpc of $M_{\rm TME}=0.4^{+1.2}_{-0.3}\times 10^{8}$ M*⊙. The total stellar mass is $M_{\rm stars}=(1.5\pm 0.4)\times 10^{8}$ M⊙* (see § 2), and we conclude that we cannot reject the hypothesis that there is no dark matter in this system. Another way to phrase this is that the intrinsic dispersion of $\sigma_{\rm intr}=4.2^{+4.4}_{-2.2}$ km s*-1* is consistent with that expected from the stars alone ($\sigma_{\rm stars}\approx 7.3^{+0.9}_{-1.1}$ km s*-1*; Wolf et al. 2010). We note that this result is independent of the precise distance to NGC1052-DF4: as $\sigma_{\rm stars}\propto{}(M_{\rm stars}/R_{\rm e})^{0.5}$ it scales with distance as $\sigma_{\rm stars}\propto(D^{2}/D)^{0.5}\propto{}D^{0.5}$. For reference, the expected dispersion from an NFW halo of the expected mass is $\sigma_{\rm halo}\approx 30$ km s*-1* (Łokas & Mamon 2001; Behroozi et al. 2013a).

5. Discussion

In this Letter we have presented a doppelgänger of the dark matter-deficient galaxy NGC1052-DF2. NGC1052-DF4 is in the same group as NGC1052-DF2 and has a similar size, luminosity, morphology, globular cluster population, and velocity dispersion. The immediate implication is that NGC1052-DF2 is not an isolated case but that a class of such galaxies exists, and given how little we know about galaxies in the UDG parameter space it may well be that they are fairly common.

The discovery of NGC1052-DF4 does not bring us much closer to understanding how such galaxies are formed, although it does effectively rule out “tail of the distribution” explanations for NGC1052-DF2. Suggestions that the true velocity dispersion is in the upper 10 % of the posterior distribution (Martin et al. 2018), that the galaxy could be an asymmetric thin disk seen exactly face-on (see van Dokkum et al. 2018b), is on a precisely-tuned orbit (Ogiya 2018), or that it was formed in a carefully orchestrated sequence of events (e.g., Fensch et al. 2019), are far less likely now that there is a second system.

One pathway for creating dark matter-deficient galaxies is by forming them out of gas that was expelled from a disk with a high baryon fraction, through a tidal interaction (e.g., Duc & Mirabel 1998; Gentile et al. 2007). Although some form of tidal origin is perhaps the most plausible explanation for these objects, NGC1052-DF2 and NGC1052-DF4 lack two of the key identifying features of “classical” tidal dwarf galaxies. First, as the gas originated in a dense disk it should be pre-enriched. Therefore, tidal dwarfs should have a high metallicity for their mass (Duc & Mirabel 1998), and this is not the case for NGC1052-DF2 (van Dokkum et al. 2018c; Fensch et al. 2019). Second, although we identify several tidal features associated with NGC 1052 in the Dragonfly imaging, there is no evidence for debris in the vicinity of NGC1052-DF2 or NGC1052-DF4, although this has been reported around other old tidal dwarfs (see Duc et al. 2014).

More broadly, the environment of NGC1052-DF4 does not shed much light on its origins. The galaxy is highlighted in a wide field of view in Fig. 4. As noted in § 4.1 the systemic velocity of NGC1052-DF4 is almost identical to the average of the NGC 1052 group galaxies. It is at a projected distance of $28\farcm 5$ (165 kpc) from NGC 1052 itself, a factor of two further than NGC1052-DF2, and $26\arcmin$ (150 kpc) from the spiral galaxy NGC 1042 (which is almost certainly also a group member; see van Dokkum et al. 2019). It is close (23 kpc) in projection to NGC 1035, which has a radial velocity of $cz=1241$ km s*-1*. Given their velocity difference of 204 km s*-1* it is unlikely that NGC1052-DF4 is a satellite of this low mass disk galaxy (Truong et al. 2017), which means that their 3D distance in the group is probably much larger than their projected distance. Apart from the relatively large systemic velocity of NGC1052-DF2 ($1805$ km s*-1*; Danieli et al. 2019) there is nothing obviously “special” about the two galaxies in relation to other group members, or about the NGC 1052 group when compared to other structures.

Although it is hazardous to draw conclusions from such correlations, it seems likely that the low dark matter mass in these galaxies is somehow related to their unprecedented globular cluster systems (see van Dokkum et al. 2018c). The combined number of confirmed clusters with $M_{V}<-8.5$ in NGC1052-DF2 and NGC1052-DF4 is nearly the same as the number of globular clusters in the Milky Way to that limit (18 versus 23), despite the factor of 100 difference in stellar mass between them (see Fig. 5). The seven confirmed clusters in NGC1052-DF4 make up 3 % of its total luminosity, and the two most distant clusters by themselves make up $\approx 70$ % of NGC1052-DF4’s luminosity at $R>5$ kpc. As shown in the bottom panel of Fig. 5 the half-light radii of the clusters (measured in the same way as described in van Dokkum et al. 2018c) are somewhat smaller than those in NGC1052-DF2, and similar to those of luminous Milky Way clusters. The median size of the seven clusters is $\langle r_{\rm h}\rangle=4.1$ pc.

Looking ahead, we can determine the stellar velocity dispersion of NGC1052-DF4 (see Emsellem et al. 2019; Danieli et al. 2019), constrain its dark matter mass (e.g., Laporte et al. 2019; Wasserman et al. 2018), and assess the implications for alternative gravity (van Dokkum et al. 2018b; Famaey, McGaugh, & Milgrom 2018): taken together, NGC1052-DF2 and NGC1052-DF4 seem in tension with recent predictions from Modified Newtonian Dynamics (Müller, Famaey, & Zhao 2019). Following the adage “one is an exception but two is a population” this new object provides impetus for characterizing the properties of diffuse, dark matter-deficient galaxies as a class. We are performing wide field surveys with the Dragonfly Telephoto Array to identify other candidates.

We thank the anonymous referees for excellent comments, including the suggestion to look for diffuse galaxy light in the sky background. Support from STScI grants HST-GO-13682 and HST-GO-14644, as well as NSF grants AST-1312376, AST-1613582, and AST-1515084 is gratefully acknowledged. AJR is a Research Corporation for Science Advancement Cottrell Scholar.