MNRAS 000, 1–13 (0000) Preprint 11 January 2021 Compiled using MNRAS LATEX style file v3.0 Dark-matter-deficient dwarf galaxies form via tidal stripping of dark matter in interactions with massive companions R. A. Jackson,1★ S. Kaviraj,1 G. Martin,2,3, J. E. G. Devriendt,4 A. Slyz,4 J. Silk,4,5,6 Y. Dubois,5 S. K. Yi,7 C. Pichon,5,8 M. Volonteri,5 H. Choi,7 T. Kimm,7 K. Kraljic,9 and S. Peirani5,10 1Centre for Astrophysics Research, School of Physics, Astronomy and Mathematics, University of Hertfordshire, Hatfield, AL10 9AB, UK 2Steward Observatory, University of Arizona, 933 N. Cherry Ave, Tucson, AZ 85719, USA 3Korea Astronomy and Space Science Institute, 776 Daedeokdae-ro, Yuseong-gu, Daejeon 34055, Korea 4Dept of Physics, University of Oxford, Keble Road, Oxford OX1 3RH UK 5Institut d’Astrophysique de Paris, Sorbonne Universités, UMPC Univ Paris 06 et CNRS, UMP 7095, 98 bis bd Arago, 75014 Paris, France 6Department of Physics & Astronomy, The Johns Hopkins University, Baltimore, MD 21218, USA 7Department of Astronomy and Yonsei University Observatory, Yonsei University, Seoul 03722, Republic of Korea 8School of Physics, Korea Institute for Advanced Study (KIAS), 85 Hoegiro, Dongdaemun-gu, Seoul, 02455, Republic of Korea 9Institute for Astronomy, Royal Observatory, Edinburgh EH9 3HJ, UK 10Observatoire de la Coˆte d’Azur, CNRS, Laboratoire Lagrange, Bd de l’Observatoire, Université Coˆte d’Azur, CS 34229, 06304 Nice Cedex 4, France 11 January 2021 ABSTRACT In the standard ΛCDM paradigm, dwarf galaxies are expected to be dark-matter-rich, as baryonic feedback is thought to quickly drive gas out of their shallow potential wells and quench star formation at early epochs. Recent observations of local dwarfs with extremely low dark matter content appear to contradict this picture, potentially bringing the validity of the standard model into question. We use NewHorizon, a high-resolution cosmological simulation, to demonstrate that sustained stripping of dark matter, in tidal interactions between amassive galaxy and a dwarf satellite, naturally produces dwarfs that are dark-matter-deficient, even though their initial dark matter fractions are normal. The process of dark matter stripping is responsible for the large scatter in the halo-to-stellar mass relation in the dwarf regime. The degree of stripping is driven by the closeness of the orbit of the dwarf around its massive companion and, in extreme cases, produces dwarfs with halo-to-stellar mass ratios as low as unity, consistent with the findings of recent observational studies. ∼30 per cent of dwarfs show some deviation from normal dark matter fractions due to dark matter stripping, with 10 per cent showing high levels of dark matter deficiency (Mhalo/M★<10). Given their close orbits, a significant fraction of dark-matter-deficient dwarfs merge with their massive companions (e.g. ∼70 per cent merge over timescales of ∼3.5 Gyrs), with the dark-matter-deficient population being constantly replenished by new interactions between dwarfs and massive companions. The creation of these galaxies is, therefore, a natural by-product of galaxy evolution and their existence is not in tension with the standard paradigm. Key words: galaxies: evolution – galaxies: formation – galaxies: interactions – galaxies:dwarf – methods: numerical 1 INTRODUCTION In the standard ΛCDM paradigm, dwarf galaxies (M★ < 109.5 M ) are expected to be dark matter (DM) rich, because their shallow po- tential wells make it easier for processes like stellar and supernova feedback to drive gas out from their central regions at early epochs. This reduces their star formation rates and produces systems that ★ E-mail: r.jackson9@herts.ac.uk have a lower stellar content at a given DM halo mass (Dubois & Teyssier 2008; Di Cintio et al. 2017; Chan et al. 2018; Jackson et al. 2020). Dwarf galaxies are expected to form a natural exten- sion to the stellar-to-halo mass relation in massive galaxies (Moster et al. 2010; Read et al. 2017), with the Mhalo/M★ values expected to progressively increase towards lower stellar masses. However, some recent observational studies appear to challenge this picture. A growing number of studies suggest that some dwarf galaxies may deviate strongly from the expected stellar-to-halo mass rela- © 0000 The Authors ar X iv :2 01 0. 02 21 9v 2 [a str o- ph .G A] 7 Ja n 2 02 1 2 R. A. Jackson et al. tion (e.g. van Dokkum et al. 2018a, 2019; Guo et al. 2019; Hammer et al. 2020), with unexpectedly low DM fractions. For example, Guo et al. (2019) have found 19 nearby dwarf galaxies that are DM deficient. The majority of these galaxies (14/19) appear to be relatively isolated, without many nearby bright galaxies, implying that they might have been born DM deficient. van Dokkum et al. (2018a, 2019) have studied two dwarf galaxies in group environ- ments which could exhibit Mhalo/M★ values (within 7.6 kpc) close to unity, suggesting that they may have DM fractions that are around 400 times lower than that expected for galaxies of their stellar mass. There is still some controversy surrounding the validity of these findings (e.g. Oman et al. 2016; Martin et al. 2018b; Blakeslee & Cantiello 2018; Emsellem et al. 2019; Fensch et al. 2019a; Trujillo et al. 2019), mostly focusing on the reliability of obtaining accurate distance measurements for these systems. For example, in the case of the van Dokkum et al. (2018a) galaxy, a measured distance of 13 Mpc, as suggested by Trujillo et al. (2019), rather than the quoted ∼ 20Mpc in van Dokkum et al. (2018a), would makeMhalo/M★ > 20, making it a relatively normal dwarf galaxy (but see van Dokkum et al. (2018b) for an alternate view). The existence of galaxies that are deficient in DM could pose a serious challenge to the ΛCDM model, as it is difficult for galaxies that are rich in baryons to form in halos that are deficient in DM. It is, therefore, important to inves- tigate whether there are natural channels for the formation of such galaxies in the standard model. Some formation methods for DM deficient galaxies have been suggested in the recent literature. Awell-known channel for forming such systems are tidal dwarfs (Wetzstein et al. 2007; Bournaud et al. 2008a,b; Kaviraj et al. 2012; Kroupa 2012; Ploeckinger et al. 2015; Haslbauer et al. 2019). These dwarf galaxies are formed in the tidal tails that emerge as a result of gas-rich major mergers of massive galaxies, either through Jeans instabilities within the gas which lead to gravitational collapse and the formation of self-bound objects (Elmegreen et al. 1993), or a large fraction of the stellar material in the progenitor disk being ejected and providing a local potential well into which gas condenses and fuels star formation (e.g. Barnes & Hernquist 1992; Duc et al. 2004; Hancock et al. 2009). However, the contribution of tidal dwarfs to the galaxy population, particularly at the stellar masses of the DM deficient galaxies found by recent observational studies (∼109 M ), is extremely small (e.g. Kaviraj et al. 2012). It has been postulated that DM deficient galaxies could form via high velocity, gas-rich mergers of dwarf galaxies themselves (Silk 2019; Shin et al. 2020). In this scenario, DM deficient galax- ies form as these mergers separate DM from the warm disc gas, which is then compressed by tidal interactions and shocks to form stars. Another potential formation mechanism may involve anoma- lously weak stellar feedback, due to low star formation rate surface densities. As a result of this, significant amounts of gas are not ejected from the disc, resulting in a galaxy which has a high baryon fraction and therefore a relatively low DM fraction (e.g Mancera Piña et al. 2020). A further formation channel is DM stripping of satellite galax- ies, particularly in extreme environments such as clusters (Ogiya 2018; Jing et al. 2019; Niemiec et al. 2019). Indeed, both N-body simulations in a range of environments, from Milky-Way mass ha- los (Hayashi et al. 2003; Kravtsov et al. 2004; Diemand et al. 2007; Rhee et al. 2017; Buck et al. 2019) to clusters (Ghigna et al. 1998; Gao et al. 2004; Tormen et al. 2004; Nagai & Kravtsov 2005; van den Bosch et al. 2005; Giocoli et al. 2008; Xie & Gao 2015), and analytical models (e.g. Mamon 2000; Gan et al. 2010; Han et al. 2016; Hiroshima et al. 2018) have shown that DM stripping is ca- pable of removing parts of a galaxy’s halo in group and cluster environments. It has been suggested that this process could drive a large scatter in the stellar-to-halo mass relation for satellite galaxies in groups and clusters, moving them from their initial positions to- wards lower halo masses (Vale & Ostriker 2004; Smith et al. 2016; Bahé et al. 2017; Rhee et al. 2017; Niemiec et al. 2017, 2019), and therefore lower Mhalo/M★ values. A comprehensive analysis of whether DM deficient systems (dwarfs in particular) can form naturally as a by-product of the process of galaxy evolution demands a hydrodynamical simulation in a cosmological volume which has both high mass and spatial resolution. The hydrodynamics is required for spatially-resolved predictions for the DM and baryons. The cosmological volume enables us to probe galaxy populations in a statistical manner, tak- ing into account environmental effects (which idealised studies, for instance, cannot do), and is particularly important for makingmean- ingful comparisons to large observational surveys (e.g. forthcoming datasets like LSST, Robertson et al. 2019). In recent years, large hydrodynamical cosmological simula- tions, e.g. Horizon-AGN (Dubois et al. 2014b), Illustris (Vogels- berger et al. 2014), EAGLE (Schaye et al. 2015) and Simba (Davé et al. 2019), have been successful in reproducing many properties of (massive) galaxies over cosmic time (e.g. Kaviraj et al. 2017). For example, Saulder et al. (2020) have investigated the properties of massive (M★ > 109.5M ) DM deficient galaxies which are iso- lated in these large cosmological simulations, akin to those found in Guo et al. (2019). They find that these galaxies are probably regular objects that undergo un-physical processes at the boundary of the simulation box and are therefore artefacts. Jing et al. (2019) on the other hand have studied the formation of massive (109 M < M★ < 1010 M ) DM deficient galaxies in the EAGLE and Illustris sim- ulations, in denser environments. They find that a non-negligible fraction (2.6 per cent in EAGLE, and 1.5 per cent in Illustris) of satellite galaxies, in large groups and clusters (M200 > 1013 M ), are DM deficient in their central regions. These galaxies, which are not initially DM deficient, become so through the stripping of DM by tidal interactions with their host galaxy, an effect that has also been noted in Horizon-AGN (Volonteri et al. 2016)). However, given that these environments are relatively rare, and since themajority of dwarfs have lower stellarmasses than the galax- ies in these studies, it remains unclear from these studies (1) whether DM deficient galaxies can form in low-density environments which host the majority of objects, and (2) whether they can form in the dwarf regime where galaxies are expected to be significantly more DM dominated at early epochs, and where most observational stud- ies that have found DM deficient galaxies are focussed. It is worth noting that it is challenging to probe dwarf galaxy evolution using the large-scale cosmological simulations mentioned above due to their relatively low mass and spatial resolutions. For example, the DM mass resolution is ∼108 M in Horizon-AGN, EAGLE and Illustris and the spatial resolution of these simulations is around a kpc. Recall that the stellar scale height of the Milky Way, for exam- ple, is ∼300 pc (e.g. Kent et al. 1991; López-Corredoira et al. 2002; McMillan 2011), which implies that much better spatial resolution is needed to properly resolve dwarfs. An accurate exploration of the evolution of dwarf galaxies re- quires a cosmological simulation with significantly better mass and spatial resolution, in order to properly resolve the processes on the small spatial scales involved. In this study, we use the NewHorizon hydrodynamical cosmological simulation, to understand the for- mation of DM deficient galaxies in the stellar mass range M★ > 108 M . NewHorizon is a zoom-in of an average region within MNRAS 000, 1–13 (0000) Dark-matter-deficient dwarf galaxies 3 Horizon-AGN simulation, which has a volume of 142 comoving Mpc3). NewHorizon offers a maximum spatial resolution of 34 pc and mass resolutions of 104 M and 106 M in the stars and DM respectively, making it an ideal tool to study the evolution of dwarf galaxies (e.g. Jackson et al. 2020; Martin et al. 2020). At the lower limit of our stellar mass range, a typical low-redshift dwarf galaxy contains ∼10,000 stellar and ∼10,000 DM particles respectively. Our aims are to (1) study whether DM deficient dwarf galaxies form naturally in NewHorizon, (2) estimate the frequency with which they are created and (3) quantify the processes that produce these galaxies. This paper is structured as follows. In Section 2, we briefly describe the NewHorizon simulation and the selection of DM defi- cient galaxies. In Section 3, we study the processes that create these objects. We summarise our findings in Section 4. 2 SIMULATION The NewHorizon cosmological, hydrodynamical simulation (Dubois et al. 2020), is a high-resolution simulation, produced us- ing a zoom-in of a region of the Horizon-AGN simulation (Dubois et al. 2014b; Kaviraj et al. 2017, H-AGN hereafter). H-AGN em- ploys the adaptive mesh refinement code RAMSES (Teyssier 2002) and utilises a grid that simulates a 142 comoving Mpc3 volume (achieving a spatial resolution of 1 kpc), using 10243 uniformly- distributed cubic cells that have a constant mass resolution, using MPGrafic (Prunet et al. 2008). For NewHorizon, we resample this grid at higher resolution (using 40963 uniformly-distributed cubic cells), with the same cos- mology as that used in H-AGN (Ω𝑚 = 0.272, Ω𝑏 = 0.0455, ΩΛ = 0.728, H0 = 70.4 km s−1 Mpc−1 and 𝑛𝑠 = 0.967; Komatsu et al. (2011)). The high-resolution volume occupied by NewHorizon is a sphere which has a radius of 10 comoving Mpc, centred on a region of average density within H-AGN. NewHorizon has a DM mass resolution of 106 M (compared to 8×107 M in H-AGN), stellar mass resolution of 104 M (compared to 2×106 M in H- AGN) and a maximum spatial resolution of 34 pc (compared to 1 kpc in H-AGN)1. The simulation has been performed down to 𝑧 = 0.25. 2.1 Star formation and stellar feedback NewHorizon employs gas cooling via primordial Hydrogen and Helium, which is gradually enriched by metals produced by stellar evolution (Sutherland & Dopita 1993; Rosen & Bregman 1995). An ambient UV background is switched on after the Universe is re-ionized at 𝑧 = 10 (Haardt & Madau 1996). Star formation is assumed to take place in gas that has a hydrogen number density greater than 𝑛𝐻 >10 H cm−3 and a temperature lower than 2×104 K, following a Schmidt relation (Schmidt 1959; Kennicutt 1998). The star-formation efficiency is dependent on the local turbulent Mach number and the virial parameter 𝛼 = 2𝐸𝑘 /|𝐸𝑔 |, where E𝑘 is the kinetic energy of the gas and E𝑔 is the gas gravitational binding energy (Kimm et al. 2017). The probability of forming a star particle of mass M∗,𝑟𝑒𝑠=104 M is drawn at each time step using a Poissonian sampling method, as described in Rasera & Teyssier (2006). Each star particle represents a coeval stellar population,with 31 1 The gravitational force softening is equal to the local grid size. per cent of the stellar mass of this population (corresponding to stars moremassive than 6M ) assumed to explode as Type II supernovae, 5 Myr after its birth. This fraction is calculated using a Chabrier initial mass function, with upper and lower mass limits of 150 M and 0.1 M (Chabrier 2005). Supernova feedback takes the form of both energy and momentum, with the final radial momentum capturing the snowplough phase of the expansion (Kimm & Cen 2014). The initial energy of each supernova is 1051 erg and the supernova has a progenitor mass of 10 M . Pre-heating of the ambient gas by ultraviolet radiation from young O and B stars is included by augmenting the final radialmomentum from supernovae following Geen et al. (2015). 2.2 Supermassive black holes and black-hole feedback Supermassive black holes (SMBHs) are modelled as sink particles. These accrete gas and impart feedback to their ambient medium via a fraction of the rest-mass energy of the accreted material. SMBHs are allowed to form in regions that have gas densities larger than the threshold of star formation, with a seed mass of 104 M . New SMBHs do not form at distances less than 50 kpc from other existing black holes. A dynamical gas drag force is applied to the SMBHs (Ostriker 1999) and two SMBHs are allowed tomerge if the distance between them is smaller than 4 times the cell size, and if the kinetic energy of the binary is less than its binding energy. The Bondi-Hoyle-Lyttleton accretion rate determines the ac- cretion rate on to SMBHs, with a value capped at Eddington (Hoyle & Lyttleton 1939; Bondi & Hoyle 1944). The SMBHs release en- ergy back into the ambient gas via a thermal quasar mode and a jet radio mode, when accretion rates are above and below 1 per cent of the Eddington rate respectively (Dubois et al. 2012). The spins of these SMBHs are evolved self-consistently through gas accretion in the quasar mode and coalescence of black hole binaries (Dubois et al. 2014a), which modifies the radiative efficiencies of the accre- tion flow (following themodels of thin Shakura&Sunyaev accretion discs) and the corresponding Eddington accretion rate, mass-energy conversion, and bolometric luminosity of the quasar mode (Shakura & Sunyaev 1973). The quasar mode imparts 15 per cent of the bolo- metric luminosity as thermal energy into the surrounding gas. The radio mode employs a spin-dependent variable efficiency and spin up and spin down rates that follow the simulations of magnetically- choked accretion discs (see e.g. McKinney et al. 2012). 2.3 Identification of galaxies and merger trees We use the AdaptaHOP algorithm to identify DM halos (Aubert et al. 2004; Tweed et al. 2009). AdaptaHOP efficiently removes subhalos from principal structures and keeps track of the fractional number of low-resolution DM particles within the virial radius of the halo in question. We identify galaxies in a similar fashion, using the HOP structure finder applied directly on star particles (Eisen- stein & Hut 1998). The difference between HOP and AdaptaHOP lies in the fact that HOP does not remove substructures from the main structure, since this would result in star-forming clumps be- ing removed from galaxies. We then produce merger trees for each galaxy in the final snapshot at 𝑧 = 0.25, with an average timestep of ∼ 15 Myr, enabling us to track the main progenitor of every galaxy with high temporal resolution. Given that NewHorizon is a high resolution zoom of Horizon- AGN, it is worth considering the DM purity of galaxies, since higher-mass DM particles may enter the high resolution region of MNRAS 000, 1–13 (0000) 4 R. A. Jackson et al. Figure 1. The stellar mass (M★) vs the ratio of the DM halo (Mhalo) and stellar mass for galaxies in NewHorizon at 𝑧 = 0.25. In this study, we focus on galaxies with M★>108 M , which are well-resolved in the simulation to early epochs (galaxies that are less massive that this threshold are shown simply for completeness). We split this population into objects that are cen- trals (open black circles) and satellites (coloured circles). The satellites are further divided into ‘DM deficient’ galaxies, defined as objects that exhibit Mhalo/M★<10 (blue circles), satellites that coincide with the relatively tight locus defined by the centrals (orange circles) and ‘DM poor’ galaxies (green circles) which fall in between these two populations. The satellites that coin- cide with the centrals (orange circles) are essentially unstripped in their DM. We use these unstripped satellites as a control sample, since these galaxies show normal levels of DM. The grey dotted line denotes the mass threshold adopted for defining dwarf galaxies. While there are some galaxies in our sample above this threshold, the overwhelming majority of galaxies in our study are dwarfs. NewHorizon from the surrounding lower-resolution regions. Given the large mass difference, these DM particles may interact with galaxies they are passing through in unusual ways. The vast majority of galaxies affected by low DM purity exist at the outer edge of the NewHorizon sphere. TheDMdeficient galaxies studied in this paper all have DM halos with a purity of 100 per cent. 2.4 Selection of galaxies that are deficient in DM Figure 1 shows the stellar mass (M★) vs the ratio of the DM halo (Mhalo) and stellar mass, for galaxies in NewHorizon at 𝑧 = 0.25. We study galaxies which have stellar masses above 108 M , which remain well-resolved in the simulation to early epochs.We note that, while there are somegalaxies in our sample above themass threshold that we use to define dwarfs (M★ < 109.5 M ), the overwhelming majority of galaxies in our study are dwarfs. We split this galaxy population into objects that are centrals (open black circles) and satellites (coloured circles). Satellite galaxies are defined in the usual way, as systems whose DM halos have been accreted by a larger halo. The satellites are further divided into ‘DM deficient’ galaxies (blue circles), defined as objects that exhibitMhalo/M★<10, satellites that coincide with the relatively tight locus defined by the centrals (orange circles) and ‘DM poor’ galaxies (green circles) Figure 2.The evolution of themedian properties of the central,DMdeficient, DM poor and control populations for M★>108 M . The top, middle and bottom panels show the evolution in the stellar mass, DM halo mass and the Mhalo/M★ ratio respectively. While the stellar mass evolution is similar across all populations, the DM halo mass evolves differently, with the DM poor and DM deficient galaxies exhibiting sharp declines in their halo mass at late epochs (i.e. look-back times less than 8 Gyrs). This, in turn, drives the scatter in the Mhalo vs M★ relation seen in Figure 1. which fall in between these populations. The satellites that occupy the same region in the M★-Mhalo/M★ plane as the centrals (orange circles) are essentially unstripped in their DM. For the analysis below, where we study the reasons for the DM stripping in the DM deficient and DM poor populations, we use these unstripped satellites as a control sample, since these galaxies show ‘normal’ levels of DM (i.e. a similar DM content to that of centrals). Note that in this studywe defineMhalo to be the totalmass of the DM halo associated with the galaxy. This is in contrast to previous studies (e.g. Jing et al. 2019) that have mostly looked at the central regions of the DM halo (e.g. two times the half-mass radius in the stars). The fractions of DM deficient, DM poor, control and central galaxies in NewHorizon are 12, 19, 18 and 51 per cent respectively. Recall that the centrals and the controls (i.e. the unstripped satellites) show normal levels of DM, so that the fraction of galaxies that show some sort of a deficiency in their DM content is around 30 per cent. In Figure 2 we show the evolution of the median properties of each population. The top, middle and bottom panels show the evolution in the stellar mass, DMhalomass and theMhalo/M★ ratios respectively. While the stellar mass evolution is similar across all populations, the DM halo mass evolves differently, with the DM MNRAS 000, 1–13 (0000) Dark-matter-deficient dwarf galaxies 5 Figure 3. The locations of our different galaxy populations (DM deficient, DM poor and control) in NewHorizon at 𝑧 = 0.25, in three mutually orthogonal projections (xy, yz and xz) of the simulation volume. The centrals are shown using small grey circles. The DM deficient and DM poor populations form in all parts of the Universe, with the DM deficient galaxies generally showing a preference for the nodes of the cosmic web. Figure 4. The local 3D number density for the different populations vs look- back time for galaxies with M★>108 M . The local 3D number density is calculated using an adaptive kernel density estimation method (see text in Section 3 for details). Recall that NewHorizon does not contain rich clusters and therefore the highest density regions in the simulation correspond to large groups, with halo masses of∼1013 M . Both the DMdeficient and DM poor populations show elevated local densities, compared to the controls, indicating that local density may play a role in their formation. deficient and DM poor galaxies exhibiting declines in DM content at late epochs. The scatter in the M★ vs Mhalo/M★ relation seen in Figure 1 is therefore driven by the removal of DM in these galaxies, which produces systems that are poorer in DM for their stellar mass. While the gradual change from the locus of centrals and con- trols seen in Figure 1 indicates that the DM deficient and DM poor populations are not simply artefacts and are likely created by a per- sistent process that operates on these galaxies, we perform several checks to ensure that this is indeed the case. The galaxies in these populations are found in all parts of the simulation sphere (as dis- cussed further in Section 3 and Figure 3) and not just at the edges of this volume where boundary effects can produce spurious effects. They also exhibit 100 per cent purity in DM, so that their evolution is not driven by interactions with massive low-resolution DM parti- cles passing through them. Finally, as shown in the analysis below (Figure 5) the reduction of DM is very gradual, rather than abrupt, which indicates that this phenomenon is not a result of misallocation of particles between DM halos. 3 FORMATION OF DM DEFICIENT GALAXIES THROUGH STRIPPING OF DM IN TIDAL INTERACTIONS We begin our analysis by exploring the local environment of our dif- ferent satellite populations (DM deficient, DM poor, and control). Figure 3 shows the positions of these populations in the cosmic web. The DM deficient and DM poor populations form in all parts of the Universe, with the DM deficient galaxies generally showing a preference for the nodes of the cosmic web. To explore local galaxy density more quantitatively, we followMartin et al. (2018a) and de- fine a 3-D local number density of objects around each galaxy. The local density is calculated using an adaptive kernel density estima- tion method2, where the width of the kernel used for multivariate density estimation is responsive to the local density of the region, such that the error between the density estimate and the true density is minimised (Breiman et al. 1977; Ferdosi et al. 2011; Martin et al. 2018a). The density estimate takes into account all galaxies with stellar masses above 108 M . In Figure 4 we show the evolution of the local density occupied by our different populations. Recall that NewHorizon does not contain rich clusters and therefore the highest density regions in the simulation correspond to large groups with halo masses of ∼1013 M . Figures 3 and 4 indicate that the DM deficient and DM poor populations form in all regions of the Universe (note again that NewHorizon does not contain clusters or large voids so these envi- ronments are not probed here). The DM deficient galaxies generally show a preference for the nodes of the cosmic web. Both the DM deficient and DM poor populations show elevated local densities, compared to the controls, indicating that local density may play a role in their formation. Not unexpectedly, the dwarf centrals reside 2 Code available at github.com/garrethmartin/MBEtree MNRAS 000, 1–13 (0000) 6 R. A. Jackson et al. Figure 5. The rows show (from top to bottom) the evolution of the stellar mass, DM halo mass, the halo level (central [level 1] or satellite [level > 1]), perturbation index (PI) and the Mhalo/M★ ratio of the galaxy in question. The PI is split into contributions from ‘similar-mass’ galaxies (those that have 0.3 Mgal < M★ < 3 Mgal, where Mgal is the stellar mass of the galaxy in question), ‘smaller galaxies’ (those that have M★ < 0.3 Mgal) and ‘larger galaxies’ (those that have M★ > 3 Mgal). The ‘median dwarf PI’ indicates the median PI experienced by galaxies with stellar masses within 0.5 dex of the galaxy in question. In the bottom panel, the dashed black line indicates the median halo-to-stellar mass ratio of galaxies with stellar masses within 0.5 dex of the galaxy in question and the shaded region is its 1𝜎 dispersion. The left-hand column shows a DM deficient galaxy, while the right-hand column shows a control galaxy. The dotted red lines indicate times when galaxy mergers, with mass ratios greater than 10:1, take place. in lower density environments, where the probability of interacting with, and becoming the satellite of, a more massive galaxy is lower. The evolution of galaxies is influenced by a combination of in- ternal processes, like supernova (e.g. Kimm& Cen 2014) and AGN (e.g. Croton et al. 2006) feedback, and external processes such as tidal perturbations (e.g. Martin et al. 2019), mergers (e.g. Kaviraj 2014a,b; Kaviraj et al. 2019) and ram pressure (e.g. Hester 2006). Baryonic feedback and ram pressure are not capable of removing significant amounts of DM in galaxies, except in the very central regions through rapid expansion of gas (Navarro et al. 1996; Gov- ernato et al. 2012; Teyssier et al. 2013). Furthermore, the steady stripping of DM, which drives the creation of the DM poor and DM deficient populations, requires a dynamical process, suggesting that tidal perturbations (possibly including mergers) are likely to be im- portant in giving rise to these systems. In particular, interactions between a massive galaxy and a lower mass companion typically results in the stripping of material from the smaller companion. This is because, while the same tidal force acts on both objects, the smaller companion is not as strongly bound and stripping typically takes place from the outskirts of the system where the depth of the gravitational potential well is shallowest (Smith et al. 2016). To quantify the effect of tidal perturbations, we follow Martin et al. (2019) and Jackson et al. (2020) to define a dimensionless ‘perturbation index’ (PI) that quantifies the strength of the ambient tidal field around a galaxy: PI = ∑︁ i ( Mi Mhalo ) ( Rvir Di )3 , (1) where Mhalo is the halo mass of the galaxy in question, Rvir is MNRAS 000, 1–13 (0000) Dark-matter-deficient dwarf galaxies 7 Figure 6.Evolution of the DMdensity profile, for various lookback times, in the DM deficient galaxy described in the left-hand column of Figure 5. The evolution is shown from the point at which the galaxy becomes a satellite, which is also the point at which the DM stripping starts (‘max mass’), through to the end of the simulation. The DM stripping takes place in the outskirts of the dwarf where the gravitational potential well is shallowest and DM particles are more loosely bound. the virial radius of its DM halo, Mi is the DM halo mass of the 𝑖th perturbing galaxy and Di is the distance to the 𝑖th perturbing galaxy. We consider all perturbing galaxies within 3 Mpc of the object in question. We use the PI to explore the trends in the tidal perturbations in our different populations and study how these perturbations may be driving the DM stripping in the DM poor and DM deficient populations. In Figure 5, we visually illustrate the evolution of two typical objects (selected randomly) in the DM deficient and control populations. The left-hand column shows the DM deficient galaxy, while the right-hand column shows the control galaxy. Recall that the control galaxies are satellites in which DM is not stripped and that the DM deficient galaxies are systems which have the lowest values of Mhalo/M★ (Figure 1). These two populations therefore bracket the dwarf population as a whole. The rows show (from top to bottom) the evolution of the stellar mass, DM halo mass, halo level (central [level 1] or satellite [level > 1]), PI and the Mhalo/M★ ratio of the galaxy in question. The PI is split into contributions from ‘similar-mass’ galaxies (those that have 0.3 Mgal < M★ < 3 Mgal, where Mgal is the stellar mass of the galaxy in question), ‘smaller galaxies’ (those that have M★ < 0.3 Mgal) and ‘larger galaxies’ (those that have M★ > 3Mgal). While the DM deficient galaxy shows a steady increase in stellar mass (row 1), it exhibits sustained DM stripping after a look-back time of around 6 Gyrs (row 2). The stripping coincides with the galaxy transitioning from being a central to a satellite (row 3) and an increase in the PI frommassive galaxies, even as that from smaller galaxies decreases (row 4). The DM stripping drives the Mhalo/M★ ratio down, until it is below the region traced by the median ratio for galaxies which have stellar masses within ±0.5 dex (row 5). The principal difference between the DMdeficient galaxy and its control Figure 7. The total perturbation index (PI) from all companions vs that just from the massive companion, integrated over the period from when the dwarf becomes a satellite until the end of the simulation. The differ- ent satellite populations are shown using different colours. The histograms are normalised and the dashed lines indicate median values. DM deficient galaxies show both larger values for the total PI and the PI from massive companions, compared to their control counterparts (with the DM poor galaxies falling in between these two populations). Note also that the PI from the massive companion dominates the total PI in the DM deficient and DM poor populations. counterpart is that, even after the control object becomes a satellite, the PI due to massive galaxies does not increase and the control galaxy does not exhibit significant DM stripping. In Figure 6, we show the evolution of the DM density profile (calculated using all DM particles associated with the galaxy) of the DM deficient galaxy described in the left-hand column of Figure 5. The evolution is shown from the point at which the galaxy becomes a satellite (which is also the point at which the DM stripping starts), through to the end of the simulation. Not unexpectedly, the DM stripping typically takes place in the outskirts of the dwarf, where the gravitational potential well is shallower and DM particles are less well bound. The trends are identical in all DM deficient and DM poor galaxies. In Table 1, we present the four largest PI contributions from individual companions, for the DM deficient and control galaxy in Figure 5, at the point where each galaxy becomes a satellite. The left hand table, which corresponds to the DM deficient galaxy, shows that ∼80 per cent of the total PI comes from one larger companion (which is the central that hosts this dwarf), with the next largest PI contribution being an order of magnitude smaller. For the control galaxy, on the other hand, the largest PI contribution, while also derived from a larger companion (again the host central), is a factor of 5 smaller than that in the DM deficient galaxy and contributes only 25 per cent of the total PI. Indeed, many smaller companions contribute PI values which are within a factor of 2 of that from the larger companion. Note that, in all cases, the larger companion which dominates the PI in the DM deficient and DM poor galaxies is the central that hosts the dwarf satellite in question. While Figure MNRAS 000, 1–13 (0000) 8 R. A. Jackson et al. Figure 8. The DM halo (left) and stellar (centre) mass of the massive centrals that host our dwarf satellites and the halo mass ratio (right) between the dwarf satellite and its massive central. The dashed lines indicate the median values of the histograms. The distribution of massive central masses and halo mass ratios are similar for all dwarf populations and are therefore not the principal reason for the differences in perturbation index (PI) in these populations. Figure 9. The orbits of four galaxies in each of the DM deficient (shades of blue), DM poor (shades of green) and control (shades of orange) populations around their centrals. Individual galaxy orbits are shown using different shades of the colour in question (see legend). We show orbits from the point at which the dwarf in question becomes a satellite through to the end of the simulation. The orbital distances are normalised by the virial radius of the central, which is always shown at the origin. PI PI fraction Cumulative PI Mass ratio 0.48 0.78 0.78 5.75 0.05 0.08 0.86 0.01 0.01 0.02 0.88 0.88 0.01 0.01 0.89 0.01 PI PI fraction Cumulative PI Mass ratio 0.09 0.25 0.25 10.98 0.04 0.12 0.37 0.67 0.04 0.12 0.49 0.02 0.02 0.04 0.53 0.01 Table 1. The four largest perturbation index (PI) contributions from individual companions for the DM deficient (left hand table) and control (right hand table) galaxy shown in Figure 5, at the point where the galaxy becomes a satellite. Columns in each table show (from left to right) the PI contributed by the individual companion, the fraction this represents of the total PI, the cumulative fraction of PI contributed by this companion and those in previous rows and the mass ratio of the companion and the dwarf galaxy in question. The left hand table, which corresponds to the DM deficient galaxy, shows that ∼80 per cent of the total PI comes from one larger companion with the next largest PI contribution being an order of magnitude smaller. For the control galaxy the largest PI contribution, while also derived from a larger companion, is a factor of 5 smaller than that in the DM deficient galaxy and contributes only 25 per cent of the total PI, with many smaller companions contributing PI values which are within a factor of 2 of that from the larger companion. MNRAS 000, 1–13 (0000) Dark-matter-deficient dwarf galaxies 9 Figure 10. Left: Minimum distance between dwarfs and their massive centrals for our different satellite populations. The dashed lines indicate the median value of each distribution. Right: Orbital distances (normalised by the virial radius of the central) between the dwarf and their massive central at all timesteps in the period between the dwarf becoming a satellite and the end of the simulation. The DM deficient and DM poor populations exhibit closer orbits than their control counterparts which drives the higher perturbation index (PI) values seen in these populations. Figure 11. 𝑟 -band mock images of two example systems that bracket the types of interactions that lead to DM stripping. Each image is 200 proper kpc across. In all images, the DM deficient galaxy in question is at the centre of the image and indicated by the red arrow. The blue circle shows the virial radius of the DM deficient galaxy. The sequence (where time moves forward from left to right) shows ∼3 Gyrs of evolution. The top row shows a two body interaction, while the bottom row shows a more complex interaction where the dwarf interacts with a larger companion which dominates a group. In both cases, the dwarf undergoes significant DM stripping and ends up as a DM deficient galaxy. MNRAS 000, 1–13 (0000) 10 R. A. Jackson et al. 5 illustrates the trends seen in typical galaxies in the DM deficient and control populations, Figure 7 summarises the differences in the PI from massive galaxies seen in our different satellite populations, after the galaxies become satellites. We show both the total PI from all companions vs that just from the massive companion, integrated over the period fromwhen the dwarf becomes a satellite until the end of the simulation. DM deficient galaxies show both larger values for the total PI and the PI from massive companions, compared to their control counterparts (with the DM poor galaxies falling in between these two populations). Note also that the PI from the massive companion dominates the total PI in the DM deficient and DM poor populations. This indicates that the DM stripping in the DM deficient and DM poor galaxies is indeed driven by stronger interactions with a massive (central) companion. We proceed by exploring why the DM deficient dwarfs exhibit higher values of PI due to their massive companions. Figure 8 indicates that both the halo (left panel) and stellar (centre panel) masses of the central that host the dwarfs are similar. In addition, the mass ratios between the dwarfs and their host centrals (right panel), are comparable across all the different populations. Therefore, these quantities are unlikely to be driving the higher values of PI in the DM deficient and DM poor populations. However, given the strong dependence of the PI on the distance between galaxies (Equation 1), it is likely that the higher PI in DM deficient and DM poor galaxies is driven by orbits that bring them closer to their corresponding centrals. In Figure 9, we first illustrate this graphically by plotting the orbits of four galaxies (selected to span the mass range in our study) in each of the DM deficient, DM poor and control populations around their centrals. The orbital distances are normalised by the virial radius of the central and orbits are shown from the point at which the dwarf in question becomes a satellite to the end of the simulation. It is clear that DM deficient galaxies have orbits that bring them significantly closer to their massive central than in their control counterparts. Figure 10 presents this more quantitatively, by comparing both the orbital distances (normalised by the virial radius of the cen- tral) in the period after individual galaxies become satellites, and the minimum distances between the satellites and centrals during these orbits. DM deficient dwarfs spend most of their orbits at sig- nificantly smaller distances from their centrals compared to their control counterparts, with the DM poor galaxies lying in between these two populations. The same patterns are apparent in the mini- mum orbital distances. The larger PI values from larger galaxies in these populations are therefore driven by orbits which bring these galaxies closer to their massive central companions. The degree of DM stripping is correlated with these orbital distances, with members of the DM deficient population (which exhibits the largest amount of stripping) spending their orbits at much smaller distances than their DM poor and control counterparts. Given that the DM deficient galaxies are in tight orbits around nearby massive companions, it is worth considering how long these dwarfs may survive. We explore this by considering the evolution of the DM deficient population (selected in an identical way) at 𝑧 = 0.7 i.e. around 3.5 Gyrs before the epoch (𝑧 = 0.25) at which the analysis above is performed. The stellar mass distribution of DM deficient galaxies at both epochs is similar. We find that only 30 per cent of DM deficient galaxies that exist at 𝑧 = 0.7 still survive at 𝑧 = 0.25. This indicates that the creation of DM deficient galaxies is a constant process over cosmic time. In other words, satellites that are in close orbits are not only stripped of their DM, but many also do not survive after the stripping starts as they are accreted by their larger companions. Thus, the DM deficient population at a given redshift are largely systems that have formed recently enough that they still exist in the simulation at that epoch. We also consider whether dwarfs that may be candidates for being DM deficient systems could potentially be identified using quantities that are readily available in imaging survey data. As de- scribed above, the process that drives the creation of these objects are tidal interactions with more massive companions on very close orbits. In the appendix we show a version of the right-hand panel of Figure 10 without normalising the orbital distances of the dwarfs by the virial radius of the massive companion (Figure A1), because virial radii are not measurable quantities in imaging data. Combin- ing Figures 8, 10, 11 and A1 indicates that dwarfs that are found close to massive companions with M★ > 1010 M (Figure 8) at distances less ∼150 kpc (Figure A1) and show stellar tidal features (Figure 11) are likely to be good candidates for being DM deficient systems. As is typical of interactions between massive galaxies and dwarf companions (e.g. Kaviraj 2010, 2014b), the tidal features (e.g. the ones visible in Figure 11) are faint with surface bright- nesses typically fainter than 29 mag arcsec−2. This suggests that finding large samples of DM deficient candidates would ideally re- quire deep wide-area surveys, such as LSST (Robertson et al. 2019; Kaviraj 2020), which have limiting surface-brightnesses that are fainter than such values. 3.1 Comparison to observational studies We complete our study by considering our results in the con- text of recent observational work. DM deficient galaxies form in NewHorizon via tidal interactions in close orbits with a more mas- sive (central) companion. It is worth noting that the theoretical study of Macciò et al. (2020) comes to similar conclusions, albeit using a small sample of three zoom-in simulations of individual galax- ies. The formation mechanism described in our study is strongly supported by the recent empirical work of Montes et al. (2020), who have used deep optical imaging to demonstrate that the DM deficient galaxy NGC 1052-DF4 exhibits tidal tails, similar to those predicted in Figure 11. They show that these tidal features are driven by an interaction with its neighboring galaxy NGC 1035, and argue that the DM-deficient nature of this dwarf galaxy is caused by the DM stripping produced by this interaction. Our proposed mechanism is also consistent with the recent studies by van Dokkum et al. (2018a) and van Dokkum et al. (2019), who find multiple DM deficient dwarf galaxies in group environ- ments. These galaxies are likely to have been subjected to strong tidal forces that would have stripped their DM halo (see e.g. the lower panel of Figure 11 which shows a possible analog of this in NewHorizon). The formation of DMdeficient dwarfs via tidal inter- actions also appears consistent with the excess number of globular clusters (GCs) reported in these systems (e.g. Fensch et al. 2019a; Müller et al. 2020). Since tidal interactions between dwarfs and larger companions can trigger enhanced GC formation (e.g. Fensch et al. 2019b; Carleton et al. 2020; Somalwar et al. 2020), systems like DM deficient dwarfs, that undergo strong tidal interactions, could be expected to show an excess number of GCs. Note, how- ever, that the resolution of NewHorizon is not sufficient for us to directly study the formation of GCs. The formation channel outlined in this paper appears less well- aligned with the findings of Guo et al. (2019), who suggest that most of their DM deficient galaxies lie at distances greater than three times the virial radius of the nearest group or cluster. However, it is worth noting that the nearest groups and clusters in Guo et al. are significantly more massive than those in our study. It is possible that MNRAS 000, 1–13 (0000) Dark-matter-deficient dwarf galaxies 11 the tidal forces required to strip dwarfs of their DM can be produced at larger distances aroundmuchmore massive groups. Furthermore, without deep images it is difficult to identify tidal features which are the tell-tale signatures of the tidal stripping process that creates DMdeficient galaxies. Nevertheless, the relatively isolated nature of these galaxies, coupled with other potential mechanisms for produc- ing such objects, e.g. weak baryonic feedback (Mancera Piña et al. 2020) suggests that other channels maybe required to satisfactorily explain the entire DM deficient population. 4 SUMMARY In the standard ΛCDM paradigm, dwarf galaxies are expected to be DM-rich, because their shallow gravitational potential wells make it easier for processes like stellar and supernova feedback to deplete their gas reservoirs. This results in a rapid depletion of gas and reduction of star formation at early epochs, leaving these objects with relatively high DM fractions. However, recent observational work suggests that some local dwarfs exhibit DM fractions as low as unity, around 400 times lower than what is expected for systems of their stellar mass. The existence of such DM deficient galaxies appears to contradict our classical expectations of theDMproperties of dwarf galaxies, potentially bringing the validity of the standard paradigm into serious question. Understanding the origins of these galaxies, using a high-resolution cosmological simulationwhich can make realistic statistical predictions of dwarf galaxies, is, therefore, a key exercise. Here, we have used the NewHorizon cosmological hydrody- namical simulation to explain the formation of DM deficient galax- ies. We are able to perform this exercise, for the first time, in a statis- tical fashion, as the cosmological volume of NewHorizon allows us to study large numbers of dwarf galaxies, while its high spatial res- olution enables us to resolve these systems with the requisite detail. We have shown that interactions between massive central galaxies and dwarf satellites can drive sustained and significant stripping of DM from the dwarfs, which reduces their DM content. The level of stripping is determined by the details of the orbit, with dwarfs that are heavily stripped typically spending significant fractions of their time in tight orbits around their correspondingmassive central, after they turn into satellites. DM stripping is responsible for the large dispersion in the stellar-to-halo mass relation in the dwarf regime (see Figure 1), with the DM deficient and DM poor populations, which comprise ∼30 per cent of dwarfs, scattering off the tight main locus defined by the dwarf centrals and the dwarf satellites that remain unstripped. In extreme cases, this DM stripping produces dwarfs which exhibit Mhalo/M★ ratios as low as unity, consistent with the findings of recent observational studies. Given their close orbits, a significant fraction of DM deficient dwarfs will merge with their massive com- panions and disappear from the galaxy population (e.g.∼70 per cent of such dwarfs will merge over timescales of ∼3.5 Gyrs). But the DM deficient population is replenished by new interactions between dwarfs and massive companions. It is worth noting that our results are robust with respect to the details of the sub-grid recipes imple- mented in NewHorizon, as this formation mechanism is principally driven by gravitational forces. The formation of DM deficient galaxies through tidal stripping of DM, as hypothesised here, is strongly supported by recent ob- servational studies (e.g. Montes et al. 2020), which are located in similar environments to the DM deficient galaxies in NewHorizon. Observations of an excess of GCs in these galaxies, although not directly testable in NewHorizon, also appear to support this forma- tion mechanism (although, as noted above, contributions from other pathways may be required to explain the DM deficient population as a whole). The results of our study offer a route to identifying dwarf galaxies that may have low DM content, purely from data that is typically available in imaging surveys. Dwarfs that are found close to massive companions with M★ > 1010 M , at distances less than ∼150 kpc, and show stellar tidal features, are likely to be good candidates for being DM deficient systems. The surface brightness of the tidal features produced by such interactions (typically fainter than 29 mag arcsec−2) suggests that new and future deep wide surveys, like the Hyper Suprime-Cam SSP or LSST, could be used to identify large samples of such DM deficient dwarf candidates. Our study demonstrates that stripping of DM via a tidal inter- action, a process that takes place in all environments, can routinely create dwarfs that have DM fractions that deviate significantly from their initial values in the early Universe and, in extreme cases, pro- duces systems that are DM deficient. The existence of such galaxies is therefore an integral feature of galaxy evolution in the standard ΛCDM paradigm and their existence is not in tension with the pre- dictions of this model. ACKNOWLEDGEMENTS We are grateful to the anonymous referee for many constructive comments that helped us to improve the quality of the original manuscript. RAJ and SK acknowledge support from the STFC [ST/R504786/1, ST/S00615X/1]. SK acknowledges a Senior Re- search Fellowship from Worcester College Oxford. The research of AS and JD is supported by the Beecroft Trust and STFC. TK was supported in part by the National Research Foundation of Ko- rea (NRF-2017R1A5A1070354 and NRF-2020R1C1C100707911) and in part by the Yonsei University Future-leading Research Ini- tiative (RMS2-2019-22-0216). Some of the numerical work made use of the DiRAC Data Intensive service at Leicester, operated by the University of Leicester IT Services, which forms part of the STFC DiRAC HPC Facility (www.dirac.ac.uk). The equip- ment was funded by BEIS capital funding via STFC capital grants ST/K000373/1 and ST/R002363/1 and STFC DiRAC Operations grant ST/R001014/1.DiRAC is part of theNational e-Infrastructure. This research has used the DiRAC facility, jointly funded by the STFCand the Large Facilities Capital Fund ofBIS, and has been par- tially supported by grant Segal ANR- 19-CE31-0017 of the French ANR. This work was granted access to the HPC resources of CINES under the allocations 2013047012, 2014047012 and 2015047012 made by GENCI and was granted access to the high-performance computing resources of CINES under the allocations c2016047637 and A0020407637 from GENCI, and KISTI (KSC-2017-G2-0003). Large data transfer was supported by KREONET, which is managed and operated by KISTI. This work has made use of the Horizon cluster on which the simulation was post-processed, hosted by the Institut d’Astrophysique de Paris. We thank Stephane Rouberol for running it smoothly for us. DATA AVAILABILITY The data that underpins the analysis in this paper were provided by theHorizon-AGNandNewHorizon collaborations. The datawill be shared on request to the corresponding author, with the permission MNRAS 000, 1–13 (0000) 12 R. A. Jackson et al. of theHorizon-AGNandNewHorizon collaborations. Local density estimation code is available at github.com/garrethmartin/MBEtree. REFERENCES Aubert D., Pichon C., Colombi S., 2004, MNRAS, 352, 376 Bahé Y. M., et al., 2017, MNRAS, 470, 4186 Barnes J. E., Hernquist L., 1992, Nature, 360, 715 Blakeslee J. P., Cantiello M., 2018, Research Notes of the American Astro- nomical Society, 2, 146 Bondi H., Hoyle F., 1944, MNRAS, 104, 273 Bournaud F., Bois M., Emsellem E., Duc P. A., 2008a, Astronomische Nachrichten, 329, 1025 Bournaud F., Duc P. A., Emsellem E., 2008b, MNRAS, 389, L8 Breiman L., Meisel W., Purcell E., 1977, Technometrics, 19, 135 Buck T., Macciò A. V., Dutton A. A., Obreja A., Frings J., 2019, MNRAS, 483, 1314 Carleton T., Guo Y., Munshi F., Tremmel M., Wright A., 2020, arXiv e- prints, p. arXiv:2008.11205 Chabrier G., 2005, The Initial Mass Function: From Salpeter 1955 to 2005. p. 41, doi:10.1007/978-1-4020-3407-7_5 Chan T. K., Kereš D., Wetzel A., Hopkins P. F., Faucher-Giguère C. A., El-Badry K., Garrison-Kimmel S., Boylan-Kolchin M., 2018, MNRAS, p. 1094 Croton D. J., et al., 2006, MNRAS, 365, 11 Davé R., Anglés-Alcázar D., Narayanan D., Li Q., Rafieferantsoa M. H., Appleby S., 2019, MNRAS, 486, 2827 Di Cintio A., Brook C. B., Dutton A. A., Macciò A. V., Obreja A., Dekel A., 2017, MNRAS, 466, L1 Diemand J., Kuhlen M., Madau P., 2007, ApJ, 667, 859 Dubois Y., Teyssier R., 2008, A&A, 482, L13 Dubois Y., Devriendt J., Slyz A., Teyssier R., 2012, MNRAS, 420, 2662 Dubois Y., Volonteri M., Silk J., 2014a, MNRAS, 440, 1590 Dubois Y., et al., 2014b, MNRAS, 444, 1453 Dubois Y., et al., 2020, arXiv e-prints, p. arXiv:2009.10578 Duc P. A., Bournaud F., Masset F., 2004, A&A, 427, 803 Eisenstein D. J., Hut P., 1998, ApJ, 498, 137 Elmegreen B. G., Kaufman M., Thomasson M., 1993, ApJ, 412, 90 Emsellem E., et al., 2019, A&A, 625, A76 Fensch J., et al., 2019a, A&A, 625, A77 Fensch J., et al., 2019b, A&A, 628, A60 Ferdosi B. J., Buddelmeijer H., Trager S. C., Wilkinson M. H. F., Roerdink J. B. T. M., 2011, A&A, 531, A114 Gan J., Kang X., van den Bosch F. C., Hou J., 2010, MNRAS, 408, 2201 Gao L., White S. D. M., Jenkins A., Stoehr F., Springel V., 2004, MNRAS, 355, 819 Geen S., Rosdahl J., Blaizot J., Devriendt J., Slyz A., 2015, MNRAS, 448, 3248 Ghigna S., Moore B., Governato F., Lake G., Quinn T., Stadel J., 1998, MNRAS, 300, 146 Giocoli C., Tormen G., van den Bosch F. C., 2008, MNRAS, 386, 2135 Governato F., et al., 2012, MNRAS, 422, 1231 Guo Q., et al., 2019, Nature Astronomy, 4, 246 Haardt F., Madau P., 1996, ApJ, 461, 20 Hammer F., Yang Y., Arenou F., Wang J., Li H., Bonifacio P., Babusiaux C., 2020, ApJ, 892, 3 Han J., Cole S., Frenk C. S., Jing Y., 2016, MNRAS, 457, 1208 Hancock M., Smith B. J., Struck C., Giroux M. L., Hurlock S., 2009, AJ, 137, 4643 Haslbauer M., Dabringhausen J., Kroupa P., Javanmardi B., Banik I., 2019, A&A, 626, A47 Hayashi E., Navarro J. F., Taylor J. E., Stadel J., Quinn T., 2003, ApJ, 584, 541 Hester J. A., 2006, ApJ, 647, 910 Hiroshima N., Ando S., Ishiyama T., 2018, Phys. Rev. D, 97, 123002 Hoyle F., Lyttleton R. A., 1939,Mathematical Proceedings of the Cambridge Philosophical Society, 35, 592–609 Jackson R. A., et al., 2020, arXiv e-prints, p. arXiv:2007.06581 Jing Y., Wang C., Li R., Liao S., Wang J., Guo Q., Gao L., 2019, MNRAS, 488, 3298 Kaviraj S., 2010, MNRAS, 406, 382 Kaviraj S., 2014a, MNRAS, 437, L41 Kaviraj S., 2014b, MNRAS, 440, 2944 Kaviraj S., 2020, arXiv e-prints, p. arXiv:2001.01728 Kaviraj S., Darg D., Lintott C., Schawinski K., Silk J., 2012, MNRAS, 419, 70 Kaviraj S., et al., 2017, MNRAS, 467, 4739 Kaviraj S., Martin G., Silk J., 2019, MNRAS, 489, L12 Kennicutt Jr. R. C., 1998, ApJ, 498, 541 Kent S. M., Dame T. M., Fazio G., 1991, ApJ, 378, 131 Kimm T., Cen R., 2014, ApJ, 788, 121 Kimm T., Katz H., Haehnelt M., Rosdahl J., Devriendt J., Slyz A., 2017, MNRAS, 466, 4826 Komatsu E., et al., 2011, ApJS, 192, 18 Kravtsov A. V., Gnedin O. Y., Klypin A. A., 2004, ApJ, 609, 482 Kroupa P., 2012, Publ. Astron. Soc. Australia, 29, 395 López-Corredoira M., Cabrera-Lavers A., Garzón F., Hammersley P. L., 2002, A&A, 394, 883 Macciò A. V., Huterer Prats D., Dixon K. L., Buck T., Waterval S., Arora N., Courteau S., Kang X., 2020, arXiv e-prints, p. arXiv:2010.02245 Mamon G. A., 2000, in Combes F., Mamon G. A., Charmandaris V., eds, Astronomical Society of the Pacific Conference Series Vol. 197, Dy- namics of Galaxies: from the Early Universe to the Present. p. 377 (arXiv:astro-ph/9911333) Mancera Piña P. E., et al., 2020, MNRAS, 495, 3636 Martin G., Kaviraj S., Devriendt J. E. G., Dubois Y., Pichon C., Laigle C., 2018a, MNRAS, 474, 3140 Martin N. F., Collins M. L. M., Longeard N., Tollerud E., 2018b, ApJ, 859, L5 Martin G., et al., 2019, MNRAS, 485, 796 Martin G., et al., 2020, arXiv e-prints, p. arXiv:2007.07913 McKinney J. C., Tchekhovskoy A., Blandford R. D., 2012, Monthly Notices of the Royal Astronomical Society, 423, 3083 McMillan P. J., 2011, MNRAS, 414, 2446 Montes M., Infante-Sainz R., Madrigal-Aguado A., Román J., Monelli M., Borlaff A. S., Trujillo I., 2020, ApJ, 904, 114 Moster B. P., Somerville R. S., Maulbetsch C., van den Bosch F. C., Macciò A. V., Naab T., Oser L., 2010, ApJ, 710, 903 Müller O., et al., 2020, A&A, 640, A106 Nagai D., Kravtsov A. V., 2005, ApJ, 618, 557 Navarro J. F., Eke V. R., Frenk C. S., 1996, MNRAS, 283, L72 Niemiec A., et al., 2017, MNRAS, 471, 1153 Niemiec A., Jullo E., Giocoli C., Limousin M., Jauzac M., 2019, MNRAS, 487, 653 Ogiya G., 2018, MNRAS, 480, L106 Oman K. A., Navarro J. F., Sales L. V., Fattahi A., Frenk C. S., Sawala T., Schaller M., White S. D. M., 2016, MNRAS, 460, 3610 Ostriker E. C., 1999, ApJ, 513, 252 Ploeckinger S., Recchi S., Hensler G., Kroupa P., 2015, MNRAS, 447, 2512 Prunet S., Pichon C., Aubert D., Pogosyan D., Teyssier R., Gottloeber S., 2008, ApJS, 178, 179 Rasera Y., Teyssier R., 2006, A&A, 445, 1 Read J. I., Iorio G., Agertz O., Fraternali F., 2017, MNRAS, 467, 2019 Rhee J., Smith R., Choi H., Yi S. K., Jaffé Y., Candlish G., Sánchez-Jánssen R., 2017, ApJ, 843, 128 Robertson B. E., et al., 2019, Nature Reviews Physics, 1, 450 Rosen A., Bregman J. N., 1995, ApJ, 440, 634 Saulder C., Snaith O., Park C., Laigle C., 2020, MNRAS, 491, 1278 Schaye J., et al., 2015, MNRAS, 446, 521 Schmidt M., 1959, The Astrophysical Journal, 129, 243 Shakura N. I., Sunyaev R. A., 1973, A&A, 500, 33 Shin E.-j., Jung M., Kwon G., Kim J.-h., Lee J., Jo Y., Oh B. K., 2020, ApJ, 899, 25 Silk J., 2019, MNRAS, 488, L24 MNRAS 000, 1–13 (0000) Dark-matter-deficient dwarf galaxies 13 Figure A1. Orbital distances of our different satellite populations, in the period after individual dwarfs become satellites through to the end of the simulation. The dashed lines indicate the median value of each distribution. Smith R., Choi H., Lee J., Rhee J., Sanchez-Janssen R., Yi S. K., 2016, ApJ, 833, 109 Somalwar J. J., Greene J. E., Greco J. P., Huang S., Beaton R. L., Goulding A. D., Lancaster L., 2020, arXiv e-prints, p. arXiv:2008.02806 Sutherland R. S., Dopita M. A., 1993, ApJS, 88, 253 Teyssier R., 2002, A&A, 385, 337 Teyssier R., Pontzen A., Dubois Y., Read J. I., 2013, MNRAS, 429, 3068 Tormen G., Moscardini L., Yoshida N., 2004, MNRAS, 350, 1397 Trujillo I., et al., 2019, MNRAS, 486, 1192 Tweed D., Devriendt J., Blaizot J., Colombi S., Slyz A., 2009, A&A, 506, 647 Vale A., Ostriker J. P., 2004, MNRAS, 353, 189 Vogelsberger M., et al., 2014, Nature, 509, 177 Volonteri M., Dubois Y., Pichon C., Devriendt J., 2016, MNRAS, 460, 2979 Wetzstein M., Naab T., Burkert A., 2007, MNRAS, 375, 805 Xie L., Gao L., 2015, MNRAS, 454, 1697 van Dokkum P., et al., 2018a, Nature, 555, 629 van Dokkum P., Danieli S., Cohen Y., Romanowsky A. J., Conroy C., 2018b, ApJ, 864, L18 van Dokkum P., Danieli S., Abraham R., Conroy C., Romanowsky A. J., 2019, ApJ, 874, L5 van den Bosch F. C., Tormen G., Giocoli C., 2005, MNRAS, 359, 1029 APPENDIX A: ORBITAL DISTANCES OF DWARFS AROUND THEIR MASSIVE CENTRALS Figure A1 shows the orbital distances of dwarfs around their mas- sive centrals in our different satellite populations. This figure is a version of the right-hand panel of Figure 10, without the distance being normalised by the virial radius of the halo of the massive cen- tral. These distances are directly comparable to measured distances between massive galaxies and nearby dwarfs in deep-wide surveys in order to identify DM deficient galaxy candidates. This paper has been typeset from a TEX/LATEX file prepared by the author. MNRAS 000, 1–13 (0000)