Reconciling the Diversity and Uniformity of Galactic Rotation Curves with Self-Interacting Dark Matter
Total Page:16
File Type:pdf, Size:1020Kb
PHYSICAL REVIEW X 9, 031020 (2019) Reconciling the Diversity and Uniformity of Galactic Rotation Curves with Self-Interacting Dark Matter † Tao Ren,1 Anna Kwa,2 Manoj Kaplinghat,2,* and Hai-Bo Yu1, 1Department of Physics and Astronomy, University of California, Riverside, California 92521, USA 2Department of Physics and Astronomy, University of California, Irvine, California 92697, USA (Received 23 October 2018; revised manuscript received 25 April 2019; published 7 August 2019) Galactic rotation curves exhibit diverse behavior in the inner regions while obeying an organizing principle; i.e., they can be approximately described by a radial acceleration relation or the modified Newtonian dynamics phenomenology. We analyze the rotation curve data from the SPARC sample and explicitly demonstrate that both the diversity and uniformity are naturally reproduced in a hierarchical structure formation model with the addition of dark matter self-interactions. The required concentrations of the dark matter halos are fully consistent with the concentration-mass relation predicted by the Planck cosmological model. The inferred stellar mass-to-light (3.6 μm) ratios scatter around 0.5 M⊙=L⊙,as expected from population synthesis models, leading to a tight radial acceleration relation and a baryonic Tully-Fisher relation. The inferred stellar-halo mass relation is consistent with the expectations from abundance matching. These results provide compelling arguments in favor of the idea that the inner halos of galaxies are thermalized due to dark matter self-interactions. DOI: 10.1103/PhysRevX.9.031020 Subject Areas: Astrophysics, Cosmology ϒ ¼ 0 5 I. INTRODUCTION constant stellar mass-to-light ratio ⋆;disk . M⊙=L⊙ and ϒ⋆ ¼ 0.7 M⊙=L⊙ in the 3.6-μm band. The scatter Galactic rotation curves of spiral galaxies show a variety ;bulge of behaviors in the inner parts, even across systems with in this radial acceleration relation (RAR) is around similar halo and stellar masses, which lacks a self-con- 0.1 dex, and the tightness of this relation has been sistent explanation in the standard cold dark matter (CDM) interpreted as a signature of MOND [19]. model [1–13]. Along with this diversity, a long-standing It has long been argued that the acceleration scale observation is that many rotation curves can be understood (including the cH0 dependence) can emerge from hierar- in terms of modified Newtonian dynamics (MOND) chical structure formation predicted in CDM [20,21]. phenomenology [14,15] (see Ref. [16] for a review); Recent hydrodynamical simulations of galaxy formation i.e., there exists a characteristic gravitational acceleration with CDM have clearly shown that a RAR emerges −10 2 [22–24]. However, these simulated galaxies do not re- scale, g† ≈ 10 m=s ∼ cH0=7, with H0 being the present Hubble expansion rate, below which the observed present the full range of diversity in the SPARC data set, pffiffiffiffiffiffiffiffiffiffiffiffi and they cannot yet explain the rotation curves of low- and acceleration can be approximated as g†g , with g bar bar high-surface-brightness galaxies simultaneously. being the baryonic acceleration (a.k.a. Milgrom’slaw). In this paper, we show that self-interacting dark matter More recently, McGaugh et al. [17] analyzed the Spitzer (SIDM) provides a unified way to understand the diverse Photometry and Accurate Rotation Curves (SPARC) data rotation curves of spiral galaxies while reproducing set [18] and showed that there is a tight relation between the RAR with a small scatter. We analyze the SPARC the total gravitational acceleration at any radius and the data set based on the SIDM halo model proposed in acceleration contributed by the baryons, assuming a Refs. [25,26] and demonstrate three key observations leading to this result. σ ∼ 1 2 *Corresponding author. (i) For cross section per unit mass =m cm =g, [email protected] dark matter self-interactions thermalize the inner † Corresponding author. regions at distances less than about 10% of the virial [email protected] radius of galactic halos, while the outer regions remain unchanged. Thus, SIDM inherits essential Published by the American Physical Society under the terms of features of the ΛCDM hierarchical structure for- the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to mation model such as the halo concentration-mass the author(s) and the published article’s title, journal citation, relation, which sets the characteristic acceleration and DOI. scale of halos. 2160-3308=19=9(3)=031020(12) 031020-1 Published by the American Physical Society REN, KWA, KAPLINGHAT, and YU PHYS. REV. X 9, 031020 (2019) (ii) In the inner halo, thermalization ties dark matter and We match this isothermal profile to a Navarro-Frenk- baryon distributions together [25,27,28], and the White (NFW) form [33,34] at r1, where a dark matter SIDM halo can naturally accommodate the diverse particle has scattered Oð1Þ timesovertheageofthe range of “cored” and “cusped” central density pro- galaxy, assuming continuity in both the density and the files, depending on how the baryons are distributed. enclosed mass at r1.Inthisway,theisothermal ρ σ Combined with the scatter in the concentration-mass parameters ( 0, v0) directly map onto the NFW param- ρ relation, this provides the diversity required to ex- eters (rs, s)or(rmax, Vmax). This model provides an plain the rotation curves [26,29,30]. approximate way to calculate the SIDM distribution in a (iii) For the same σ=m that addresses the diversity halo if its CDM counterpart is known, and vice versa. It problem, the baryon content of the galaxies and correctly predicts the halo central density and its the mass model of their host halos also lead to the scalings with the outer halo properties, stellar profiles, RAR with a scatter as small as the one in Ref. [17].In andcrosssection,asconfirmedinbothisolatedand ϒ our SIDM fits, the inferred stellar ⋆;disk values for cosmological N-body simulations with and without individual galaxies have a distribution peaked toward baryons; see, e.g., Refs. [26,28,30,35,36].Seethe 0.5 M⊙=L⊙, as expected from stellar population Appendix and the Supplemental Material [32] for a synthesis models [31]. detailed description of the model and additional com- The rest of the paper is organized as follows. In parisons between model predictions and cosmological Sec. II, we present the SIDM fits to 135 galaxies from simulations. the SPARC sample, which exemplify the full range of We adopt two independent but complementary the diversity. In Sec. III, we show the radial acceleration approaches to perform the analysis. In the controlled relation and the distribution of the stellar mass-to-light sampling (CS) approach, we demand that the host halos ratios from our SIDM fits, compared to the MOND fits. follow the concentration-mass relation within a 2σ range In Sec. IV, we discuss the host halo properties and the predicted in cosmological simulations [37]. We model the origin of the acceleration scale. In Sec. V,weshowthe stellar distribution as an axisymmetric thin disk as in predicted stellar-halo mass relation and the baryonic Ref. [29], which directly enters into the calculation of Tully-Fisher relation (BTFR). We comment on future the density profile of SIDM through the gravitational directions and conclude in Sec. VI. In the Appendix, we potential ΦðR; zÞ. In the CS fits, we start with the outer provide detailed information about the model and the NFW halo and find the SIDM density profile that matches fitting procedure. In the Supplemental Material [32],we its mass and density at r1. In the second approach, we use present SIDM and MOND fits to 135 individual galaxies the Markov chain Monte Carlo (MCMC) sampling (MS) to from the SPARC sample and additional results that explore the full likelihood. To save computational time, we support the main text, including model fits to simulated assume spherical symmetry by spreading the mass within halos. the disk at radius R into a sphere of the same radius [25,28]. The rotation curves generated from two approaches II. THE DIVERSITY OF GALACTIC agree well, and the differences in the fits are small (see ROTATION CURVES Supplemental Material [32]). For our main results, we show inferences from both of the approaches. We select 135 out of 175 galaxies in the full SPARC In Fig. 1, we show the SIDM fits to the diverse rotation sample based on the criteria that they must have a curves from the controlled sampling with σ=m ¼ 3 cm2=g. recorded value for the flat part of the rotation curve, Vf . In each panel, galaxies are selected to have similar flat In our sample, 87, 42, and 6 galaxies have quality rotation velocities at their outermost data points. The rise up flags 1, 2, and 3, respectively. It spans a wide range to Vf within their central regions displays a wide variety of of galaxy masses and inner shapes of rotation curves slopes, and the SIDM halo model provides equally good 20 300 with Vf ranging from km=sto km=s. In fitting fits to the shallow and steeply rising rotation curves. The to the data, we utilize the analytical SIDM halo fits for the other galaxies in the sample are as good as those model [26,29], where we assume the dark matter in Fig. 1 (see Supplemental Material [32]). distribution in the inner halo follows the isothermal The success of the SIDM halo model stems from a density profile, combination of the following effects. First, SIDM thermal- ization ties the baryon and dark matter distributions together. ρ ð Þ¼ρ (½Φ ð0 0Þ − Φ ð Þ σ2 ) ð Þ iso R; z 0 exp tot ; tot R; z = v0 ; 1 For low-surface-brightness galaxies, thermalization leads to a shallow density core and a circular velocity profile that ρ σ – where 0 is the central dark matter density, v0 is the rises mildly with radius [38 44].