The isolated dark matter-poor galaxy that ran away An example from IllustrisTNG

,


Introduction
The discovery of galaxies lacking dark matter (van Dokkum et al. 2018a(van Dokkum et al. , 2019) ) has sparked a debate on recent challenges to the standard cosmological model (ΛCDM).Yang et al. (2020) argued that self-interacting dark matter is a favorable scenario for explaining the formation and existence of such galaxies.Questioning the uncertainties of the data, modified Newtonian dynamics (MOND) has also been proposed (e.g., Famaey et al. 2018;Haghi et al. 2019).Also, Martin et al. (2018) questioned the validity of the data, arguing that the discovered galaxy is a typical dwarf, but this claim was rejected (e.g., van Dokkum et al. 2018b;Wasserman et al. 2018;Keim et al. 2022).Furthermore, Saulder et al. (2020) warned that isolated dark matter-deficient galaxies in simulations could represent a numerical artifact (arising from problems with the algorithm identifying bound structures and the periodic boundary conditions of the simulation boxes).
Numerous studies have since shown that these exotic galaxies are not in tension with the ΛCDM model and that they may naturally occur due to stronger tidal interactions (e.g., Ogiya 2018;Montes et al. 2020;Shin et al. 2020;Jackson et al. 2021;Macciò et al. 2021;Ogiya et al. 2022;Trujillo-Gomez et al. 2022).In particular, Moreno et al. (2022) estimated that ∼30% massive central galaxies have at least one dark matter-deficient satellite with stellar masses in the 10 8 −10 9 M range.Exploring the evolution of satellite galaxies in IllustrisTNG simulation boxes, Montero-Dorta et al. (2023) also considered dark matter-deficient systems and reinforced the notion that the tidal stripping is responsible for the lack of dark matter.
In this context, the recent discovery (Comerón et al. 2023) of a massive NGC 1277 galaxy lacking dark matter, while also lacking a more massive companion, was puzzling.These authors considered multiple possible scenarios to explain this phenomenon: either the galaxy lost its dark matter early or it had been lacking ab initio.To test the viability of the former assumption, we set our aims on the search for candidates in IllustrisTNG.From the theoretical point of view, the existence of backsplash galaxies (Borrow et al. 2023, and references within), namely, ones that had resided inside a cluster and have since migrated back outside, would favor the possibility that galaxies may, in fact, "run away" from dense environments.Moreover, Chilingarian & Zolotukhin (2015) argued that the observed isolated compact elliptical and isolated quiescent dwarf galaxies are tidally stripped systems that had run away from their host clusters.
Here, we present one such example from IllustrisTNG300: an isolated dark matter-poor galaxy that lost its extended dark matter envelope early.Since then, it has dwelled in the cluster environment, interacting with the cluster members without merging, while accelerating on its orbit; eventually, it left the cluster, spending the last ∼2 Gyr in isolation.This Letter is organized as follows.We describe our candidate selection criteria in Sect. 2. In Sect.3, we present the evolutionary path of the galaxy separated into stages.We discuss our results, along with the limitations and caveats of the simulation in Sect. 4. We summarize our work and give our concluding remarks in Sect. 5.

Candidate selection
We utilized the IllustrisTNG cosmological simulation suite1 (Weinberger et al. 2017;Pillepich et al. 2018a;Nelson et al. 2019).Choosing an appropriate simulation box was not a straightforward task.The largest simulation box, TNG300 (Marinacci et al. 2018;Naiman et al. 2018;Nelson et al. 2018;Pillepich et al. 2018b;Springel et al. 2018), has the poorest particle mass resolution.However, due to its size, it allows us to explore galaxies experiencing vastly different environments.Knowing that dark, matter-poor galaxies are not always low-mass objects (e.g., Comerón et al. 2023), the resolution limit should not pose a challenge.Moreover, the TNG300 box offers the largest sample of galaxies, which would help identify extremely rare present-day isolated dark, matter-poor galaxies.
We describe our selection criteria in the following.First, we considered only the subhalos of cosmological origin (that have formed due to the process of structure formation and collapse, and have SubhaloFlag = 12 ).Then, at z = 0, we extracted the ones that are isolated (i.e., they are the only subhalo within their parent halo) and that have the dark-matter-tobaryon-mass ratio of M DM /M B < 2. While the filter yielded 448 subhalos, the majority of them are poorly resolved.We chose the best-resolved, most massive subhalo, ID32549123 (herein, "the galaxy").The galaxy has, at z = 0, a total mass of M 1.56 × 10 10 M , with the mass of stellar component of M 6.8×10 9 M , no gas, and M DM /M B 1.31.The respective half-mass radii are R 0.5,DM = 2.95 kpc and R 0.5, = 2.45 kpc4 .

Results
While the galaxy still contains a non-negligible amount of dark matter at present, the total mass of the dark matter component is comparable to its stellar counterpart, indicating it is a dark matter-poor galaxy.Additionally, about half of the total dark matter mass is confined within the stellar half-mass radius, suggesting a history of substantial stripping.Due to the outside-in nature of the tidal stripping (e.g., Diemand et al. 2007;Choi et al. 2009), the particles residing in the core could not be stripped without a significant disruption of the stellar component and its mass-loss (e.g., Smith et al. 2016;Łokas 2020).In the following, we present and discuss the galaxy's cosmological evolution according to its subsequent stages.

Mergers and early history
The galaxy formed at z = 10, with M 2.6 × 10 9 M and M DM /M B 5.88.Until z = 7.01, just before its merger phase, it acquired a substantial amount of dark matter through accretion, resulting in M 9.7 × 10 9 M and M DM /M B 10.81.
According to the SubLink (Rodriguez-Gomez et al. 2015) data, this galaxy merged with 16 subhalos.The merger phase spanned over 17 early snapshots and lasted for 1.84 Gyr.The galaxy had only one major merger, with a mass ratio of q 0.599 at z = 3.71.However, given that subsequent mergers have mass ratios significantly below q = 0.1, these may have been misclassified.They are likely to be accreted non-cosmological subhalos (i.e., clumps).Checking the merger history catalog (Rodriguez-Gomez et al. 2017;Eisert et al. 2023) confirms this assumption.The galaxy has had two mergers in total: one major and one minor, which were concluded in nearly subsequent snapshots.Considering this and typical merger timescales, these are not likely to be two separate events (i.e., binary mergers) but, rather, a single multiple merger.
After this phase, at z = 2.44 (just before the first infall), the galaxy reached M 1.17 ×10 11 M , with M DM /M B 6.9 and with the total masses of sub-components (dark matter, gas, and stars) of: M DM 1.03 × 10 11 M , M G 1.43 × 10 10 M , and M 5.94×10 8 M .The respective half-mass radii are: R 0.5,DM = 19.71kpc, R 0.5,G = 20.82kpc and R 0.5, = 3.36 kpc.The galaxy has also acquired a black hole particle, with M • 1.82 × 10 7 M .

First infall and intense phase
The galaxy starts its first infall into a small group with M 1.92 × 10 12 M , containing a little over 20 subhalos (herein, "members," when referring to subhalos within groups and clusters) at z = 2.32.As the infall starts, the galaxy immediately interacts with a group of already interacting galaxies, one of which is the most massive galaxy of the group.In this process, the galaxy abruptly gets stripped of about half its total mass from the outskirts, both dark matter and gas.
The evolution of the relevant parameters are shown in Fig. 1.Since the galaxy is almost constantly tidally disturbed during this intense phase, as a measure of its center, we used the position of a particle with the lowest potential energy, instead of a calculated center of mass.This provides a more accurate estimate, as the difference between the two can become significant due to tidal structures (e.g., Mitrašinović 2022).Aside from being tidally disturbed, during this period, the galaxy occasionally and briefly becomes a satellite of more massive cluster members as it interacts with them.This can be roughly estimated from the irregularities in its orbit and relative velocity curve.
During its first pericentric passage, the galaxy accretes some dark matter and gaseous particles, and the star formation rate (SFR) increases, producing stars whose total mass consequentially increases.At z = 1.6 (snapshot 38 and lookback time of 9.77 Gyr), the group merges with a nearby cluster of M 1.61 × 10 13 M containing a little over 100 members.Upon this event, the galaxy loses the remaining majority of its dark matter (at z = 1.53), while accreting a considerable amount of gas on its way to the first pericenter (z = 1.3, snapshot 43, and lookback time of 8.99 Gyr).This triggers a burst of star formation, sharply increasing the stellar mass to essentially presentday value.During this period, the dark matter to baryon mass ratio M DM /M B drops as low as 0.52.While a fraction of accreted gas is transformed into stars, its majority gets stripped between the second and third pericentric passage, and the galaxy becomes gas-free.As a consequence, M DM /M B fluctuates around 1 until the present (i.e., z = 0).
It is worth mentioning that the galaxy loses its acquired black hole particle near the first pericenter in the cluster.However, this is not a physical process but stems from the black hole re-centering scheme employed in the simulation.The problem was highlighted by Borrow et al. (2023) and it is known to be quite common among backsplash galaxies.

Cluster dwelling and the isolated phase
We go on to consider the next phase when the host cluster merges with a nearby massive cluster at z = 0.58 (snapshot L1, page 2 of 6 time spent in this cluster, the galaxy blazes through, reaching the first and only pericenter around z = 0.31 (snapshot 77 and lookback time 3.62 Gyr) and its relative velocity peaking at v S −v G 3016.95 km s −1 .Around the pericenter, the galaxy loses a small fraction of its stellar mass and almost all its dark matter particles at the outskirts.On its way to the apocenter, it accretes back some stellar and dark matter particles until it engages in an interaction with multiple cluster members near the cluster outskirts at d/R 200 0.84.Upon that event, the galaxy suffers mild to moderate tidal stripping, accelerates, diverts its orbit, and gets on its way out of the cluster.This interaction lasts for no more than four consecutive snapshots, making it almost impossible to examine its characteristics in detail.To isolate the interaction's contribution from global cluster effects, we approximate the cluster contribution with smooth NFW potential (Navarro et al. 1997), which corresponds to the calculated cluster density profile.We calculated the galaxy's orbit in this idealized smooth cluster potential and the results are shown in Fig. 3.The orbit is diverted during the interaction, and after it ends, the galaxy gets kicked outside the cluster, as the difference between the realistic and expected distances becomes considerably high.At the same time, the L1, page 3 of 6 galaxy accelerates, and its real relative velocity compared to the expected values is about twice as high at all times, peaking at 1371 km s −1 when the expected relative velocity is 640 km s −1 .This process remarkably resembles the gravitational slingshot mechanism.
While it is still gravitationally bound to the host, the cluster merges with the most massive other cluster in the simulation (M 2.64 × 10 15 M and a little over 12 000 members) at z = 0.2 (snapshot 84 and lookback time = 2.48 Gyr).The galaxy is located on the opposite side (from where the cluster merger occurs), at a distance of d/R 200 3.06 from the new host.Since the cluster merger is not an instantaneous process, in a few subsequent snapshots, we can see that the galaxy is still gravitationally bound to its previous host, until it gets detected as isolated at z = 0.15 (snapshot 87 and lookback time 2 Gyr).
Since for ∼2 Gyr, it is dwelling in isolation, just outside the most massive cluster in the simulation box, with essentially no changes in its mass and with its stellar half-mass radius slightly stretching.Chaotic variations in the dark matter halfmass radius, R 0.5,DM , are due to particle resolution.The number of dark matter particles fluctuates around 140, as the galaxy constantly captures and loses single particles at its outskirts.The present-day stellar photometric data for SDSS bands, in particular g and r, are taken from Nelson et al. (2018).With (g − r) = (−23.153+ 24.002) 0.85, the galaxy is quite red and almost perfectly spherical, according to the stellar mass tensor obtained from the supplementary catalog (Genel et al. 2015).An independent calculation of the stellar mass tensor, going back to the time when the stellar component becomes well-resolved for such task, confirms that the galaxy remains almost perfectly spherical at all times.The result is expected, as the stars form at the very core of a tidally stripped dark matter halo.

Discussion
The galaxy presented in this work has had quite an interesting evolution, which cannot be considered typical in any way.With M DM /M B 1.31, it certainly is dark matter-poor, but the dark matter still accounts for a little over 50% of the total mass (with a large fraction of it confined within the stellar half-mass radius).This discrepancy with the observed dark matter-deficient galaxies could be because simulated galaxies typically have higher central dark matter fractions (e.g., Lovell et al. 2018), which are not easily stripped.Particle resolution can also play a role and this mass ratio should thus be considered an upper limit.
While its present-day color, shape, and extent make it easy to argue that the galaxy should be classified as a spherical dwarf (dSph), its stellar component is overly massive for a dwarf galaxy (Revaz & Jablonka 2018).Hence, we may claim the galaxy is a massive compact galaxy (MCG), but this would not satisfy the compact observational criterion given by Barro et al. (2013) or a more liberal one given by Damjanov et al. (2015), nor is it massive enough overall (e.g., Lohmann et al. 2023).Given that almost the entirety of its stellar component is ∼9 Gyr old and it has largely remained unchanged since its formation, the galaxy can represent an example of a so-called relic galaxy (Ferré-Mateu et al. 2015).Observationally confirmed relics (also known as red nuggets) are typically more massive than our example galaxy (Lisiecki et al. 2023), but this assumption may be due to bias.
Interestingly, since the galaxy lost the majority of its dark matter very early, it does not represent a tidally stripped remnant of a massive galaxy.Instead, its present-day stellar component formed in a single, brief starburst episode at the core of an already stripped dark matter halo.Even though the galaxy does not satisfy the criteria for compact-size galaxies, we argue that the early loss of an extended dark matter and the gaseous envelope has made it compact enough to avoid mergers in higherdensity environments.Moreover, since the stripping happened very early on, massive galaxy clusters had not yet formed.This, in addition to its fairly compact size, has helped the galaxy avoid mergers and limited its interactions to members of proto-clusters and later clusters through non-merger events, as it has continued to accelerate on its orbit.After one such interaction that appears to have played out as essentially a gravitational slingshot, the galaxy left the most massive cluster in the simulation.This is the first example found in simulations of both an isolated dark matter-poor galaxy that lost its extended envelope early and a fairly compact stellar system that managed to run away.As the loss of the extended envelope is related to the tidal stripping and the galaxy never crossed the boundaries of the simulation box, the example is physically grounded and does not represent a numerical artifact (Saulder et al. 2020).

Summary and conclusions
After the discovery of dark matter-deficient galaxies and the numerous studies that followed, a consensus formed around the notion that these exotic objects could naturally occur in the ΛCDM model due to stronger tidal interactions through the stripping process.They are typically satellites, with stellar masses in the 10 8 −10 9 M range, of more massive galaxies.Thus, the recent discovery of a massive galaxy lacking dark matter, without a more massive neighbor, was puzzling.The authors L1, page 4 of 6 responsible for the discovery discussed two possible scenarios for this to happen, one of which was that the galaxy lost its dark matter early.To test the viability of this scenario, we searched for examples in IllustrisTNG300 cosmological simulation.In this work, we present the evolutionary path of one such example.The first dark matter-poor galaxy in an isolated environment found in cosmological simulations, which is physically grounded and does not represent a numerical artifact.
The galaxy has, at present, a stellar mass of M 6.8 × 10 9 M , no gas, M DM /M B 1.31, and stellar half-mass radius of R 0.5, = 2.45 kpc.While its characteristics do not match those of the NGC 1277 galaxy, its evolutionary path is a perfect example of how a galaxy that loses its extended dark matter and gaseous envelope early, becoming dark matter-poor, can manage to run away from cluster environment and massive neighbors.
The galaxy lost the majority of its dark matter early (starting with a multi-body non-merger encounter), between z = 2.32 and z = 1.53 (about 9.6 Gyr ago), even before it formed its presentday stellar component.On its way to the first pericentric passage in the first proto-cluster, around 8.99 Gyr ago, it accreted a considerable amount of gas, triggering a burst of star formation and sharply increasing the stellar mass to essentially present-day value.While a fraction of accreted gas was then transformed into stars, the majority was stripped between the second and third pericentric passage and the galaxy eventually became gas-free (∼6.67 Gyr ago).As a consequence, M DM /M B fluctuates around 1 until present (i.e., z = 0).
After the host protocluster merged with a nearby massive galaxy cluster ∼5.72 Gyr ago, the galaxy blazed through the new cluster, seemingly unnoticed, with a single pericentric passage.On its way to the apocenter, it had a non-merger interaction with multiple cluster members near the outskirts, then accelerating, diverting its orbit, and eventually leaving the cluster.Since then, for the last ∼2 Gyr, it has remained largely unchanged and isolated, residing just outside the most massive cluster in the simulation.
For this scenario to be viable, a galaxy has to lose its extended dark matter envelope very early, that is, before massive galaxy clusters formed.In this way, a galaxy has become compact enough early on (although not necessarily compact by standard definitions and constraints), which has allowed it to avoid mergers in dense environments and to interact with cluster members only through non-merger events.During these interactions, it has been able to accelerate and divert its orbit, eventually ending up leaving the cluster.
The uniqueness of the galaxy presented in this work and its peculiar evolution suggest that the whole evolutionary path, particularly in the later phase, represents a series of particular circumstances.However, circumstantial events, despite being rare, are physically grounded.The upcoming public release of the MillenniumTNG simulations5 will primarily expand our knowledge and understanding of large-scale clustering as well as the formation and evolution of galaxy clusters (e.g., Bose et al. 2023;Hernández-Aguayo et al. 2023;Pakmor et al. 2023).As an appealing consequence, it will also provide fertile ground for the study of intercluster dynamics and the evolution of peculiar galaxy classes and events, as well as of environmental effects in general, with much better statistics than have been previously available.

MitrašinovićFig. 1 .Fig. 2 .
Fig.1.Total mass of galactic components (dark matter, stars, gas) shown from top to bottom, with their half-mass radii, R 0.5 (different colors represent different components, as per the legend); star formation rate (SFR); relative distance from the center of the group or cluster (blue, left y-axis), and the galaxy's relative velocity (orange, right yaxis).The dashed vertical line represents the moment when the group merges with a nearby cluster.

Fig. 3 .
Fig. 3. Relative distance from the cluster's center (top) and the galaxy's relative velocity (bottom).Solid blue line (data) represents the real values, as seen in Fig. 2. Dashed black line (extrapolated) represents the data obtained by calculating galaxy's orbit in an idealized smooth cluster potential.