Using state-of-the-art high-resolution fully GPU N-body simulations, we demonstrate for the first time that the infall of a dark matter rich satellite naturally explains a present black hole offset by sub-parsecs in M31. Observational data of the tidal features provide stringent constraints on the initial conditions of our simulations. The heating of the central region of M31 by the satellite via dynamical friction entails a significant black hole offset after the first pericentric passage. After having reached its maximum offset, the massive black hole sinks towards the M31 centre due to dynamical friction and it is determined to be offset by sub-parsecs as derived by observations.

Using fully GPU $N$-body simulations, we demonstrate for the first time that subhalos sink and transfer energy via dynamical friction into the centres of dwarf galaxies. This dynamical heating kicks any central massive black hole (MBH) out to tens of parsecs, especially at early epochs ($z$=1.5-3). This mechanism helps explain the observed off-center BHs in dwarf galaxies and also predicts that off-center BHs are more common in higher mass dwarf galaxies since dynamical friction becomes significantly weaker and BHs take more time to sink back towards the centres of their host galaxies. One consequence of off-center BHs during early epochs of dwarf galaxies is to quench any BH feedback.

We run state-of-the-art high resolution fully GPU N-body simulations to show that M31 harbours a dark matter core. Observational data in the inner halo of M31 provide stringent constraints on the initial conditions of our simulations. We demonstrate that an infalling satellite on a highly eccentric orbit is at the origin of the giant stellar stream and the shell-like features of M31. The infalling satellite heats up the central parts of its host and triggers a cusp-core transition in the dark matter halo of M31, generating a universal core independently of the initial conditions. Our results imply that cores in massive galaxies are a common and natural feature of cold dark matter haloes that have been initially cuspy but have accreted subhaloes on highly eccentric orbits

We performed a series of high-resolution $N$-body simulations to examine whether dark matter candidates in the form of primordial black holes (PBHs) can solve the cusp-core problem in low-mass dwarf galaxies. If some fraction of the dark matter in low-mass dwarf galaxies consists of PBHs and the rest is cold dark matter, dynamical heating of the cold dark matter by the PBHs induces a cusp-to-core transition in the total dark matter profile. The mechanism works for PBHs in the 25-100 M$_{\sun}$ mass window, consistent with the LIGO detections, but requires a lower limit on the PBH mass fraction of 1$\%$ of the total dwarf galaxy dark matter content. The cusp-to-core transition time-scale is between 1 and 8 Gyr. This time-scale is also a constant multiple of the relaxation time between cold dark matter particles and PBHs, which depends on the mass, the mass fraction and the scale radius of the initial density profile of PBHs. We conclude that dark matter cores occur naturally in halos comprised of cold dark matter and PBHs, without the need to invoke baryonic processes.

We use a fully GPU $N$-body code to demonstrate that dark matter minihalos, as a new component of globular clusters, resolve both the timing and cusp-core problems in Fornax if the five (or six) globular clusters were recently accreted ($\leq$ 3 Gyr ago) by Fornax. Under these assumptions, infall of these globular clusters does not occur and no star clusters form in the centre of Fornax in accordance with observations. We find that crossings of globular clusters that have DM minihalos near the Fornax centre induce a cusp-to-core transition of the dark matter halo and hence resolve the cusp-core problem in this dwarf galaxy. The dark matter core size depends on the frequency of globular cluster crossings. Our simulations clearly demonstrate also that between the passages, the dark matter halo can regenerate its cusp. Moreover, our models are in good agreement with constraints on the dark matter masses of globular clusters as our clusters lose a large fraction of their initial dark matter minihalos. These results provide circumstantial evidence for the universal existence of dark matter halos in globular clusters.

We re-investigate the Fornax cusp-core problem using observational results on the spatial and mass distributions of globular clusters (GCs) in order to put constraints on the dark matter profile. We model Fornax using high-resolution N-body simulations with entirely live systems, i.e. self-gravitating systems composed of stars and dark matter, which account correctly for dynamical friction and tidal effects between Fornax and the globular clusters. We test two alternative hypotheses, which are a cored and a cuspy halo for Fornax by exploring a reasonable range of initial conditions on globular clusters. For Fornax cored dark matter halo, we derive a lower limit on the core size of $r_{c}\gtrsim$ 0.5 kpc. Contrary to many previous works, we show also that for different initial conditions, a cuspy halo is not ruled out in our simulations based on observations of Fornax globular clusters.

Fornax is the most massive of the Milky Way dwarf spheroidal galaxies and has five globular clusters orbiting in a dense background of dark matter. Observational analyses suggest that globular clusters were initially much more massive and lost most of their stars to the Fornax field. We re-investigate the Fornax cusp-core problem, to clarify tensions between simulations and observations concerning the dark matter halo density profile. N-body simulations predict a centrally steep power-law density profile, while observations of the globular clusters seem to prefer that the dark matter halo density is constant at the center. For the first time, we ran pure N-body simulations with entirely live systems, i.e. self-gravitating systems composed of particles (including stars and dark matter). Only this numerical approach accounts correctly for dynamical friction and tidal effects between Fornax and the globular clusters. We show that a weak cusp ($r_{s}=$ 1.5 kpc) or a large core ($r_{c}=$ 848 pc) are not compatible with the current observed positions and masses of Fornax clusters. In contrast, a small dark matter core ($r_{c}=$ 282 pc) for Fornax naturally reproduces the cluster spatial and mass distributions over a wide range of initial globular cluster masses. We derive an upper limit of $r_{c}\lesssim$ 282 pc. This core size range favors only warm dark datter (WDM). It is also possible to obtain a compatible core size range from cold dark matter (CDM) theories, if the initial halo's central cusp is heated by gas to form a small core.