skip to main content

Title: Amplified J-factors in the Galactic Centre for velocity-dependent dark matter annihilation in FIRE simulations

We use FIRE-2 zoom cosmological simulations of Milky Way size Galaxy haloes to calculate astrophysical J-factors for dark matter annihilation and indirect detection studies. In addition to velocity-independent (s-wave) annihilation cross-sections 〈σv〉, we also calculate effective J-factors for velocity-dependent models, where the annihilation cross-section is either p-wave (∝ v2/c2) or d-wave (∝ v4/c4). We use 12 pairs of simulations, each run with dark matter-only (DMO) physics and FIRE-2 physics. We observe FIRE runs produce central dark matter velocity dispersions that are systematically larger than in DMO runs by factors of ∼2.5–4. They also have a larger range of central (∼400 pc) dark matter densities than the DMO runs (ρFIRE/ρDMO ≃ 0.5–3) owing to the competing effects of baryonic contraction and feedback. At 3 deg from the Galactic Centre, FIRE J-factors are 3–60 (p-wave) and 10–500 (d-wave) times higher than in the DMO runs. The change in s-wave signal at 3 deg is more modest and can be higher or lower (∼0.3–7), though the shape of the emission profile is flatter (less peaked towards the Galactic Centre) and more circular on the sky in FIRE runs. Our results for s-wave are broadly consistent with the range of assumptions in most indirect detection studies. We observe more » p-wave J-factors that are significantly enhanced compared to most past estimates. We find that thermal models with p-wave annihilation may be within range of detection in the near future.

« less
 ;  ;  ;  ;  ;  ;  ;  ;  
Award ID(s):
1910346 1752913 2045928 2108318
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
p. 55-70
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT We explore the isothermal total density profiles of early-type galaxies (ETGs) in the IllustrisTNG simulation. For the selected 559 ETGs at z = 0 with stellar masses $10^{10.7}\, \mathrm{M}_{\odot } \leqslant M_{\ast } \leqslant 10^{11.9}\, \mathrm{M}_{\odot }$, the total power-law slope has a mean of 〈γ′〉 = 2.011 ± 0.007 and a scatter of $\sigma _{\gamma ^{\prime }} = 0.171$ over the radial range 0.4–4 times the stellar half-mass radius. Several correlations between γ′ and galactic properties including stellar mass, effective radius, stellar surface density, central velocity dispersion, central dark matter fraction, and in situ-formed stellar mass ratio are compared to observations and other simulations, revealing that IllustrisTNG reproduces many correlation trends, and in particular, γ′ is almost constant with redshift below z = 2. Through analysing IllustrisTNG model variations, we show that black hole kinetic winds are crucial to lowering γ′ and matching observed galaxy correlations. The effects of stellar winds on γ′ are subdominant compared to active galactic nucleus (AGN) feedback, and differ due to the presence of AGN feedback from previous works. The density profiles of the ETG dark matter haloes are well described by steeper than NFW profiles, and they are steeper in the full physics (FP) run than their counterpartsmore »in the dark matter-only (DMO) run. Their inner density slopes anticorrelate (remain constant) with the halo mass in the FP (DMO) run, and anticorrelate with the halo concentration parameter c200 in both the types of runs. The dark matter haloes of low-mass ETGs are contracted whereas high-mass ETGs are expanded, suggesting that variations in the total density profile occur through the different halo responses to baryons.« less
  2. ABSTRACT We analyse the cold dark matter density profiles of 54 galaxy haloes simulated with Feedback In Realistic Environments (FIRE)-2 galaxy formation physics, each resolved within $0.5{{\ \rm per\ cent}}$ of the halo virial radius. These haloes contain galaxies with masses that range from ultrafaint dwarfs ($M_\star \simeq 10^{4.5}\, \mathrm{M}_{\odot }$) to the largest spirals ($M_\star \simeq 10^{11}\, \mathrm{M}_{\odot }$) and have density profiles that are both cored and cuspy. We characterize our results using a new, analytic density profile that extends the standard two-parameter Einasto form to allow for a pronounced constant density core in the resolved innermost radius. With one additional core-radius parameter, rc, this three-parameter core-Einasto profile is able to characterize our feedback-impacted dark matter haloes more accurately than other three-parameter profiles proposed in the literature. To enable comparisons with observations, we provide fitting functions for rc and other profile parameters as a function of both M⋆ and M⋆/Mhalo. In agreement with past studies, we find that dark matter core formation is most efficient at the characteristic stellar-to-halo mass ratio M⋆/Mhalo ≃ 5 × 10−3, or $M_{\star } \sim 10^9 \, \mathrm{M}_{\odot }$, with cores that are roughly the size of the galaxy half-light radius, rc ≃ 1−5 kpc. Furthermore,more »we find no evidence for core formation at radii $\gtrsim 100\ \rm pc$ in galaxies with M⋆/Mhalo < 5 × 10−4 or $M_\star \lesssim 10^6 \, \mathrm{M}_{\odot }$. For Milky Way-size galaxies, baryonic contraction often makes haloes significantly more concentrated and dense at the stellar half-light radius than DMO runs. However, even at the Milky Way scale, FIRE-2 galaxy formation still produces small dark matter cores of ≃ 0.5−2 kpc in size. Recent evidence for a ∼2 kpc core in the Milky Way’s dark matter halo is consistent with this expectation.« less
  3. ABSTRACT We study effects of heating by dark matter (DM) annihilation on black hole gas accretion. We observe that, for reasonable assumptions about DM densities in spikes around supermassive black holes, as well as DM masses and annihilation cross-sections within the standard WIMP model, heating by DM annihilation may have an appreciable effect on the accretion on to Sgr A* in the Galactic Centre. Motivated by this observation we study the effects of such heating on Bondi accretion, i.e. spherically symmetric, steady-state Newtonian accretion on to a black hole. We consider different adiabatic indices for the gas, and different power-law exponents for the DM density profile. We find that typical transonic solutions with heating have a significantly reduced accretion rate. However, for many plausible parameters, transonic solutions do not exist, suggesting a breakdown of the underlying assumptions of steady-state Bondi accretion. Our findings indicate that heating by DM annihilation may play an important role in the accretion onto supermassive black holes at the centre of galaxies, and may help explain the low accretion rate observed for Sgr A*.

    As part of the mm-Wave Interferometric Survey of Dark Object Masses (WISDOM), we present a measurement of the mass of the supermassive black hole (SMBH) in the nearby early-type galaxy NGC 0383 (radio source 3C 031). This measurement is based on Atacama Large Millimeter/sub-millimeter Array (ALMA) cycle 4 and 5 observations of the 12CO(2–1) emission line with a spatial resolution of 58 × 32 pc2 (0.18 arcsec × 0.1 arcsec). This resolution, combined with a channel width of 10 km s−1, allows us to well resolve the radius of the black hole sphere of influence (measured as RSOI = 316 pc  =  0.98 arcsec), where we detect a clear Keplerian increase of the rotation velocities. NGC 0383 has a kinematically relaxed, smooth nuclear molecular gas disc with weak ring/spiral features. We forward model the ALMA data cube with the Kinematic Molecular Simulation (KinMS) tool and a Bayesian Markov Chain Monte Carlo method to measure an SMBH mass of (4.2 ± 0.7) × 109 M⊙, a F160W-band stellar mass-to-light ratio that varies from 2.8 ± 0.6 M⊙/L$_{\odot ,\, \mathrm{F160W}}$ in the centre to 2.4 ± 0.3 M⊙$/\rm L_{\odot ,\, \mathrm{F160W}}$ at the outer edge of the disc and a molecular gas velocity dispersion of 8.3 ± 2.1 km s−1(all 3σ uncertainties). We also detect unresolved continuum emission across the full bandwidth, consistent with synchrotron emission from an active galactic nucleus. Thismore »work demonstrates that low-J CO emission can resolve gas very close to the SMBH ($\approx 140\, 000$ Schwarzschild radii) and hence that the molecular gas method is highly complimentary to megamaser observations, as it can probe the same emitting material.

    « less
  5. Abstract

    We have observed the compact H ii region complex nearest to the dynamical center of the Galaxy, G−0.02−0.07, using ALMA in the H42α recombination line, CS J = 2–1, H13CO+J = 1–0, and SiO v = 0, J = 2–1 emission lines, and the 86 GHz continuum emission. The H ii regions HII-A to HII-C in the cluster are clearly resolved into a shell-like feature with a bright half and a dark half in the recombination line and continuum emission. The analysis of the absorption features in the molecular emission lines show that H ii-A, B, and C are located on the near side of the “Galactic center 50 km s−1 molecular cloud” (50MC), but HII-D is located on the far side of it. The electron temperatures and densities ranges are Te = 5150–5920 K and ne = 950–2340 cm−3, respectively. The electron temperatures in the bright half are slightly lower than those in the dark half, while the electron densities in the bright half are slightly higher than those in the dark half. The H ii regions are embedded in the ambient molecular gas. There are some molecular gas components compressed by a C-type shock wave around the H ii regions. From the line width of the H42αmore »recombination line, the expansion velocities of HII-A, HII-B, HII-C, and HII-D are estimated to be Vexp = 16.7, 11.6, 11.1, and 12.1 km s−1, respectively. The expansion timescales of HII-A, HII-B, HII-C, and HII-D are estimated to be tage ≃ 1.4 × 104, 1.7 × 104, 2.0 × 104, and 0.7 × 104 yr, respectively. The spectral types of the central stars from HII-A to HII-D are estimated to be O8V, O9.5V, O9V, and B0V, respectively. These derived spectral types are roughly consistent with the previous radio estimation. The positional relation among the H ii regions, the SiO molecule enhancement area, and Class-I maser spots suggest that a shock wave caused by a cloud–cloud collision propagated along the line from HII-C to HII-A in the 50MC. The shock wave would have triggered the massive star formation.

    « less