-
Dynamical friction from self-interacting dark matter
Authors:
Moritz S. Fischer,
Laura Sagunski
Abstract:
Context. Merging compact objects such as binary black holes provide a promising probe for the physics of dark matter (DM). The gravitational waves emitted during inspiral potentially allow to detect DM spikes around black holes. This is because the dynamical friction force experienced by the inspiraling black hole alters the orbital period and thus the gravitational wave signal. Aims. The dynamica…
▽ More
Context. Merging compact objects such as binary black holes provide a promising probe for the physics of dark matter (DM). The gravitational waves emitted during inspiral potentially allow to detect DM spikes around black holes. This is because the dynamical friction force experienced by the inspiraling black hole alters the orbital period and thus the gravitational wave signal. Aims. The dynamical friction arising from DM can potentially differ from the collisionless case when DM is subject to self-interactions. This paper aims to understand how self-interactions impact dynamical friction. Methods. To study the dynamical friction force, we use idealized N-body simulations, where we include self-interacting dark matter. Results. We find that the dynamical friction force for inspiraling black holes would be typically enhanced by DM self-interactions compared to a collisionless medium (ignoring differences in the DM density). At lower velocities below the sound speed, we find that the dynamical friction force can be reduced by the presence of self-interactions. Conclusions. DM self-interactions have a significant effect on the dynamical friction for black hole mergers. Assuming the Chandrasekhar formula may underpredict the deceleration due to dynamical friction.
△ Less
Submitted 29 May, 2024;
originally announced May 2024.
-
Dianoga SIDM: galaxy cluster self-interacting dark matter simulations
Authors:
Antonio Ragagnin,
Massimo Meneghetti,
Francesco Calura,
Giulia Despali,
Klaus Dolag,
Moritz S. Fischer,
Carlo Giocoli,
Lauro Moscardini
Abstract:
This work aims at assessing the impact of DM self-interactions on the properties of galaxy clusters. In particular, the goal is to study the angular dependence of the cross section by testing rare (large angle scattering) and frequent (small angle scattering) SIDM models with velocity-dependent cross sections. We re-simulate six galaxy cluster zoom-in initial conditions with a dark matter only run…
▽ More
This work aims at assessing the impact of DM self-interactions on the properties of galaxy clusters. In particular, the goal is to study the angular dependence of the cross section by testing rare (large angle scattering) and frequent (small angle scattering) SIDM models with velocity-dependent cross sections. We re-simulate six galaxy cluster zoom-in initial conditions with a dark matter only run and with a full-physics setup simulations that includes a self-consistent treatment of baryon physics. We test the dark matter only setup and the full physics setup with either collisionless cold dark matter, rare self-interacting dark matter, and frequent self-interacting dark matter models. We then study their matter density profiles as well as their subhalo population. Our dark matter only SIDM simlations agree with theoretical models, and when baryons are included in simulations, our SIDM models substantially increase the central density of galaxy cluster cores compared to full-physics simulations using collisionless dark matter. SIDM subhalo suppression in full-physics simulations is milder compared to the one found in dark matter only simulations, because of the cuspier baryionic potential that prevent subhalo disruption. Moreover SIDM with small-angle scattering significantly suppress a larger number of subhaloes compared to large angle scattering SIDM models. Additionally, SIDM models generate a broader range of subhalo concentration values, including a tail of more diffuse subhaloes in the outskirts of galaxy clusters and a population of more compact subhaloes in the cluster cores.
△ Less
Submitted 1 April, 2024;
originally announced April 2024.
-
Numerical challenges for energy conservation in N-body simulations of collapsing self-interacting dark matter haloes
Authors:
Moritz S. Fischer,
Klaus Dolag,
Hai-Bo Yu
Abstract:
Dark matter (DM) haloes can be subject to gravothermal collapse if DM is not collisionless but has strong self-interactions. When the scattering can efficiently transfer heat from the centre to the outskirts, the central region of the halo collapses and reaches densities much higher than those for collisionless DM. This phenomenon is potentially observable in studies of strong lensing. Current the…
▽ More
Dark matter (DM) haloes can be subject to gravothermal collapse if DM is not collisionless but has strong self-interactions. When the scattering can efficiently transfer heat from the centre to the outskirts, the central region of the halo collapses and reaches densities much higher than those for collisionless DM. This phenomenon is potentially observable in studies of strong lensing. Current theoretical efforts are motivated by observations of surprisingly dense substructures. A comparison with them requires accurate predictions. One method to obtain such predictions is to use N-body simulations. The collapsed haloes are extreme systems that pose severe challenges to state-of-the-art codes used to model self-interacting dark matter (SIDM). We investigate the root of such problems with a focus on energy non-conservation. We run N-body simulations with and without DM self-interactions of an isolated DM-only halo and change numerical parameters relevant to the accuracy of the simulation. We find that not only the numerical scheme for the DM self-interactions can lead to energy non-conservation but also the modelling of gravitational interaction and the time integration are problematic. The issues we found are: (a) particles changing their time step in a non-time-reversible manner; (b) the asymmetry in the tree-based gravitational force evaluation; (c) SIDM velocity kicks break the symplectic nature and time symmetry. Tuning the parameters of the simulation to achieve a high accuracy allows for conserving energy not only at early stages of the evolution but also at late ones. However, the cost of the simulations becomes prohibitively large. Some problems making the simulations of the gravothermal collapse phase inaccurate, can be overcome by choosing appropriate numerical schemes. However, others remain challenging. Our findings motivate further work on addressing these challenges.
△ Less
Submitted 1 March, 2024;
originally announced March 2024.
-
Simulations of galaxy cluster mergers with velocity-dependent, rare and frequent self-interactions
Authors:
V. M. Sabarish,
Marcus Brüggen,
Kai Schmidt-Hoberg,
Moritz S. Fischer,
Felix Kahlhoefer
Abstract:
Self-interacting dark matter (SIDM) has been proposed to solve small-scale problems in $Λ$CDM cosmology. In previous work, constraints on the self-interaction cross-section of dark matter have been derived assuming that the self-interaction cross-section is independent of velocity. However, a velocity-dependent cross-section is more natural in most theories of SIDM. Using idealized $N$-body simula…
▽ More
Self-interacting dark matter (SIDM) has been proposed to solve small-scale problems in $Λ$CDM cosmology. In previous work, constraints on the self-interaction cross-section of dark matter have been derived assuming that the self-interaction cross-section is independent of velocity. However, a velocity-dependent cross-section is more natural in most theories of SIDM. Using idealized $N$-body simulations without baryons, we study merging clusters with velocity-dependent SIDM. In addition to the usual rare scattering in the isotropic limit, we also simulate these systems with anisotropic, small-angle (frequent) scatterings. We find that the collision-less brightest cluster galaxy (BCG) has an offset from the DM peak that grows at later stages. Finally, we also extend the existing upper bounds on the velocity-independent, isotropic self-interaction cross-section to the parameter space of rare and frequent velocity-dependent self-interactions by studying the central densities of dark matter-only isolated haloes. For these upper-bound parameters, the DM-BCG offsets just after the first pericentre in the dark matter-only simulations are found to be $\leq$ 10 kpc. On the other hand, because of BCG oscillations, we speculate that the distribution of BCG offsets in a relaxed cluster is a statistically viable probe. Therefore, this motivates further studies of BCG off-centring in hydrodynamic cosmological simulations.
△ Less
Submitted 1 July, 2024; v1 submitted 11 October, 2023;
originally announced October 2023.
-
Cosmological and idealized simulations of dark matter haloes with velocity-dependent, rare and frequent self-interactions
Authors:
Moritz S. Fischer,
Lenard Kasselmann,
Marcus Brüggen,
Klaus Dolag,
Felix Kahlhoefer,
Antonio Ragagnin,
Andrew Robertson,
Kai Schmidt-Hoberg
Abstract:
Dark matter self-interactions may have the capability to solve or at least mitigate small-scale problems of the cosmological standard model, Lambda Cold Dark Matter. There are a variety of self-interacting dark matter models that lead to distinguishable astrophysical predictions and hence varying success in explaining observations. Studies of dark matter (DM) density cores on various mass scales s…
▽ More
Dark matter self-interactions may have the capability to solve or at least mitigate small-scale problems of the cosmological standard model, Lambda Cold Dark Matter. There are a variety of self-interacting dark matter models that lead to distinguishable astrophysical predictions and hence varying success in explaining observations. Studies of dark matter (DM) density cores on various mass scales suggest a velocity-dependent scattering cross-section. In this work, we investigate how a velocity dependence alters the evolution of the DM distribution for frequent DM scatterings and compare to the velocity-independent case. We demonstrate that these cases are qualitatively different using a test problem. Moreover, we study the evolution of the density profile of idealized DM haloes and find that a velocity dependence can lead to larger core sizes and different time-scales of core formation and core collapse. In cosmological simulations, we investigate the effect of velocity-dependent self-interaction on haloes and satellites in the mass range of $\approx 10^{11} - 10^{14}$ M$_\odot$. We study the abundance of satellites, density, and shape profiles and try to infer qualitative differences between velocity-dependent and velocity-independent scatterings as well as between frequent and rare self-interactions. We find that a strongly velocity-dependent cross-section can significantly amplify the diversity of rotation curves, independent of the angular dependence of the differential cross-section. We further find that the abundance of satellites in general depends on both the velocity dependence and the scattering angle, although the latter is less important for strongly velocity-dependent cross-sections.
△ Less
Submitted 18 March, 2024; v1 submitted 11 October, 2023;
originally announced October 2023.
-
A quantitative comparison between velocity dependent SIDM cross sections constrained by the gravothermal and isothermal models
Authors:
Shengqi Yang,
Fangzhou Jiang,
Andrew Benson,
Yi-Ming Zhong,
Charlie Mace,
Xiaolong Du,
Zhichao Carton Zeng,
Annika H. G. Peter,
Moritz S. Fischer
Abstract:
One necessary step for probing the nature of self-interacting dark matter (SIDM) particles with astrophysical observations is to pin down any possible velocity dependence in the SIDM cross section. Major challenges for achieving this goal include eliminating, or mitigating, the impact of the baryonic components and tidal effects within the dark matter halos of interest -- the effects of these proc…
▽ More
One necessary step for probing the nature of self-interacting dark matter (SIDM) particles with astrophysical observations is to pin down any possible velocity dependence in the SIDM cross section. Major challenges for achieving this goal include eliminating, or mitigating, the impact of the baryonic components and tidal effects within the dark matter halos of interest -- the effects of these processes can be highly degenerate with those of dark matter self-interactions at small scales. In this work we select 9 isolated galaxies and brightest cluster galaxies (BCGs) with baryonic components small enough such that the baryonic gravitational potentials do not significantly influence the halo gravothermal evolution processes. We then constrain the parameters of Rutherford and Moller scattering cross section models with the measured rotation curves and stellar kinematics through the gravothermal fluid formalism and isothermal method. Cross sections constrained by the two methods are consistent at $1σ$ confidence level, but the isothermal method prefers cross sections greater than the gravothermal approach constraints by a factor of $\sim3$.
△ Less
Submitted 26 April, 2024; v1 submitted 8 May, 2023;
originally announced May 2023.
-
The role of baryons in self-interacting dark matter mergers
Authors:
Moritz S. Fischer,
Nils-Henrik Durke,
Katharina Hollingshausen,
Claudius Hammer,
Marcus Brüggen,
Klaus Dolag
Abstract:
Mergers of galaxy clusters are promising probes of dark matter (DM) physics. For example, an offset between the DM component and the galaxy distribution can constrain DM self-interactions. We investigate the role of the intracluster medium (ICM) and its influence on DM-galaxy offsets in self-interacting dark matter models. To this end, we employ Smoothed Particle Hydrodynamics + N-body simulations…
▽ More
Mergers of galaxy clusters are promising probes of dark matter (DM) physics. For example, an offset between the DM component and the galaxy distribution can constrain DM self-interactions. We investigate the role of the intracluster medium (ICM) and its influence on DM-galaxy offsets in self-interacting dark matter models. To this end, we employ Smoothed Particle Hydrodynamics + N-body simulations to study idealized setups of equal- and unequal-mass mergers with head-on collisions. Our simulations show that the ICM hardly affects the offsets arising shortly after the first pericentre passage compared to DM-only simulations. But later on, e.g. at the first apocentre, the offsets can be amplified by the presence of the ICM. Furthermore, we find that cross-sections small enough not to be excluded by measurements of the core sizes of relaxed galaxy clusters have a chance to produce observable offsets. We found that different DM models affect the DM distribution and also the galaxy and ICM distribution, including its temperature. Potentially, the position of the shock fronts, combined with the brightest cluster galaxies, provides further clues to the properties of DM. Overall our results demonstrate that mergers of galaxy clusters at stages about the first apocentre passage could be more interesting in terms of DM physics than those shortly after the first pericentre passage. This may motivate further studies of mergers at later evolutionary stages.
△ Less
Submitted 4 July, 2023; v1 submitted 15 February, 2023;
originally announced February 2023.
-
Sensitivity of Halo Shape Measurements
Authors:
Moritz S. Fischer,
Lucas M. Valenzuela
Abstract:
Shape measurements of galaxies and galaxy clusters are widespread in the analysis of cosmological simulations. But the limitations of those measurements have been poorly investigated. In this paper, we explain why the quality of the shape measurement does not only depend on the numerical resolution, but also on the density gradient. In particular, this can limit the quality of measurements in the…
▽ More
Shape measurements of galaxies and galaxy clusters are widespread in the analysis of cosmological simulations. But the limitations of those measurements have been poorly investigated. In this paper, we explain why the quality of the shape measurement does not only depend on the numerical resolution, but also on the density gradient. In particular, this can limit the quality of measurements in the central regions of haloes. We propose a criterion to estimate the sensitivity of the measured shapes based on the density gradient of the halo and to apply it to cosmological simulations of collisionless and self-interacting dark matter. By this, we demonstrate where reliable measurements of the halo shape are possible and how cored density profiles limit their applicability.
△ Less
Submitted 14 February, 2023; v1 submitted 22 September, 2022;
originally announced September 2022.
-
Cosmological simulations with rare and frequent dark matter self-interactions
Authors:
Moritz S. Fischer,
Marcus Brüggen,
Kai Schmidt-Hoberg,
Klaus Dolag,
Felix Kahlhoefer,
Antonio Ragagnin,
Andrew Robertson
Abstract:
Dark matter (DM) with self-interactions is a promising solution for the small-scale problems of the standard cosmological model. Here we perform the first cosmological simulation of frequent DM self-interactions, corresponding to small-angle DM scatterings. The focus of our analysis lies in finding and understanding differences to the traditionally assumed rare DM (large-angle) self scatterings. F…
▽ More
Dark matter (DM) with self-interactions is a promising solution for the small-scale problems of the standard cosmological model. Here we perform the first cosmological simulation of frequent DM self-interactions, corresponding to small-angle DM scatterings. The focus of our analysis lies in finding and understanding differences to the traditionally assumed rare DM (large-angle) self scatterings. For this purpose, we compute the distribution of DM densities, the matter power spectrum, the two-point correlation function and the halo and subhalo mass functions. Furthermore, we investigate the density profiles of the DM haloes and their shapes. We find that overall large-angle and small-angle scatterings behave fairly similarly with a few exceptions. In particular, the number of satellites is considerably suppressed for frequent compared to rare self-interactions with the same cross-section. Overall we observe that while differences between the two cases may be difficult to establish using a single measure, the degeneracy may be broken through a combination of multiple ones. For instance, the combination of satellite counts with halo density or shape profiles could allow discriminating between rare and frequent self-interactions. As a by-product of our analysis, we provide - for the first time - upper limits on the cross-section for frequent self-interactions.
△ Less
Submitted 5 September, 2022; v1 submitted 4 May, 2022;
originally announced May 2022.
-
Unequal-mass mergers of dark matter haloes with rare and frequent self-interactions
Authors:
Moritz S. Fischer,
Marcus Brüggen,
Kai Schmidt-Hoberg,
Klaus Dolag,
Antonio Ragagnin,
Andrew Robertson
Abstract:
Dark matter (DM) self-interactions have been proposed to solve problems on small length scales within the standard cold DM cosmology. Here, we investigate the effects of DM self-interactions in merging systems of galaxies and galaxy clusters with equal and unequal mass ratios. We perform N-body DM-only simulations of idealized setups to study the effects of DM self-interactions that are elastic an…
▽ More
Dark matter (DM) self-interactions have been proposed to solve problems on small length scales within the standard cold DM cosmology. Here, we investigate the effects of DM self-interactions in merging systems of galaxies and galaxy clusters with equal and unequal mass ratios. We perform N-body DM-only simulations of idealized setups to study the effects of DM self-interactions that are elastic and velocity-independent. We go beyond the commonly adopted assumption of large-angle (rare) DM scatterings, paying attention to the impact of small-angle (frequent) scatterings on astrophysical observables and related quantities. Specifically, we focus on DM-galaxy offsets, galaxy--galaxy distances, halo shapes, morphology, and the phase--space distribution. Moreover, we compare two methods to identify peaks: one based on the gravitational potential and one based on isodensity contours. We find that the results are sensitive to the peak finding method, which poses a challenge for the analysis of merging systems in simulations and observations, especially for minor mergers. Large DM-galaxy offsets can occur in minor mergers, especially with frequent self-interactions. The subhalo tends to dissolve quickly for these cases. While clusters in late merger phases lead to potentially large differences between rare and frequent scatterings, we believe that these differences are non-trivial to extract from observations. We therefore study the galaxy/star populations which remain distinct even after the DM haloes have coalesced. We find that these collisionless tracers behave differently for rare and frequent scatterings, potentially giving a handle to learn about the micro-physics of DM.
△ Less
Submitted 18 January, 2022; v1 submitted 21 September, 2021;
originally announced September 2021.
-
N-body simulations of dark matter with frequent self-interactions
Authors:
Moritz S. Fischer,
Marcus Brüggen,
Kai Schmidt-Hoberg,
Klaus Dolag,
Felix Kahlhoefer,
Antonio Ragagnin,
Andrew Robertson
Abstract:
Self-interacting dark matter (SIDM) models have the potential to solve the small-scale problems that arise in the cold dark matter paradigm. Simulations are a powerful tool for studying SIDM in the context of astrophysics, but it is numerically challenging to study differential cross-sections that favour small-angle scattering, as in light-mediator models. Here, we present a novel approach to mode…
▽ More
Self-interacting dark matter (SIDM) models have the potential to solve the small-scale problems that arise in the cold dark matter paradigm. Simulations are a powerful tool for studying SIDM in the context of astrophysics, but it is numerically challenging to study differential cross-sections that favour small-angle scattering, as in light-mediator models. Here, we present a novel approach to model frequent scattering based on an effective drag force, which we have implemented into the N-body code gadget-3. In a range of test problems, we demonstrate that our implementation accurately models frequent scattering. Our implementation can be used to study differences between SIDM models that predict rare and frequent scattering. We simulate core formation in isolated dark matter haloes, as well as major mergers of galaxy clusters and find that SIDM models with rare and frequent interactions make different predictions. In particular, frequent interactions are able to produce larger offsets between the distribution of galaxies and dark matter in equal-mass mergers.
△ Less
Submitted 21 June, 2021; v1 submitted 18 December, 2020;
originally announced December 2020.
-
Time-implicit schemes in fluid dynamics? -- Their advantage in the regime of ultra-relativistic shock fronts
Authors:
Moritz S. Fischer,
Ahmad A. Hujeirat
Abstract:
Relativistic jets are intrinsic phenomena of active galactic nuclei (AGN) and quasars. They have been observed to also emanate from systems containing compact objects, such as white dwarfs, neutron stars and black hole candidates. The corresponding Lorentz factors, $Γ$, were found to correlate with the compactness of the central objects. In the case of quasars and AGNs, plasmas with $Γ$-factors la…
▽ More
Relativistic jets are intrinsic phenomena of active galactic nuclei (AGN) and quasars. They have been observed to also emanate from systems containing compact objects, such as white dwarfs, neutron stars and black hole candidates. The corresponding Lorentz factors, $Γ$, were found to correlate with the compactness of the central objects. In the case of quasars and AGNs, plasmas with $Γ$-factors larger than $8$ were detected. However, numerically consistent modelling of propagating shock-fronts with $Γ\geq 4$ is a difficult issue, as the non-linearities underlying the transport operators increase dramatically with $Γ$, thereby giving rise to a numerical stagnation of the time-advancement procedure or alternatively they may diverge completely. In this paper, we present a unified numerical solver for modelling the propagation of one-dimensional shock fronts with high Lorentz factors. The numerical scheme is based on the finite-volume formulation with adaptive mesh refinement (AMR) and domain decomposition for parallel computation. It unifies both time-explicit and time-implicit numerical schemes within the framework of the pre-conditioned defect-correction iteration solution procedure. We find that time-implicit solution procedures are remarkably superior over their time-explicit counterparts in the very high $Γ$-regime and therefore most suitable for consistent modelling of relativistic outflows in AGNs and micro-quasars.
△ Less
Submitted 26 April, 2020;
originally announced April 2020.