skip to main content


Title: Dynamics of intermediate-mass black holes wandering in the milky way galaxy using the illustris TNG50 simulation
ABSTRACT

The detection of Intermediate-Mass Black Holes (IMBHs) in dwarf galaxies is crucial to closing the gap in the wide mass distribution of black holes ($\sim 3 \, {\rm M_\odot }$ to $\sim 5 \times 10^{10} \, {\rm M_\odot }$). IMBHs originally located at the centre of dwarfs that later collide with the Milky Way (MW) could be wandering, undetected, in our Galaxy. We used TNG50, the highest resolution run of the IllustrisTNG project, to study the kinematics and dynamics of star clusters, in the appropriate mass range, acting as IMBH proxies in an MW analogue galaxy. We showed that $\sim 87{{\ \rm per\ cent}}$ of our studied IMBHs drift inward. The radial velocity of these sinking IMBHs has a median magnitude of $\sim 0.44 \, \mathrm{ckpc \, h^{-1} \, Gyr^{-1}}$ and no dependence on the black hole mass. The central $1 \, \rm ckpc \, h^{-1}$ has the highest number density of IMBHs in the galaxy. A physical toy model with linear drag forces was developed to explain the orbital circularization with time. These findings constrain the spatial distribution of IMBHs, suggesting that future searches should focus on the central regions of the Galaxy. Additionally, we found that the 3D velocity distribution of IMBHs with respect to the galactic centre has a mean of $\sim 180 \, \mathrm{km \, s^{-1}}$ and larger variance with decreasing radius. Remarkably, the velocity distribution relative to the local gas shows significantly lower values, with a mean of $\sim 88 \, \mathrm{km \, s^{-1}}$. These results are instrumental for predicting the accretion and radiation properties of IMBHs, facilitating their detection with future surveys.

 
more » « less
NSF-PAR ID:
10362664
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
511
Issue:
2
ISSN:
0035-8711
Page Range / eLocation ID:
p. 2229-2238
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    Intermediate-mass black holes (IMBHs, $10^{3\!-\!6} \, {\rm M_\odot }$), are typically found at the centre of dwarf galaxies and might be wandering, thus far undetected, in the Milky Way (MW). We use model spectra for advection-dominated accretion flows to compute the typical fluxes, in a range of frequencies spanning from radio to X-rays, emitted by a putative population of $10^5 \, {\rm M_\odot }$ IMBHs wandering in five realistic volume-weighted MW environments. We predict that $\sim 27{{\ \rm per\ cent}}$ of the wandering IMBHs can be detected in the X-ray with Chandra, $\sim 37{{\ \rm per\ cent}}$ in the near-infrared with the Roman Space Telescope, $\sim 49{{\ \rm per\ cent}}$ in the sub-mm with CMB-S4, and $\sim 57{{\ \rm per\ cent}}$ in the radio with ngVLA. We find that the brightest fluxes are emitted by IMBHs passing through molecular clouds or cold neutral medium, where they are always detectable. We propose criteria to facilitate the selection of candidates in multiwavelength surveys. Specifically, we compute the X-ray to optical ratio (αox) and the optical to sub-mm ratio, as a function of the accretion rate of the IMBH. We show that at low rates the sub-mm emission of IMBHs is significantly higher than the optical, UV, and X-ray emission. Finally, we place upper limits on the number N• of these objects in the MW: N• < 2000 and N• < 100, based on our detectability expectations and current lack of detections in molecular clouds and cold neutral medium, respectively. These predictions will guide future searches of IMBHs in the MW, which will be instrumental to understanding their demographics and evolution.

     
    more » « less
  2. ABSTRACT A possibility of DM being multicomponent has a strong implication on resolving decades-long known cosmological problems on small scale. In addition to elastic scattering, the model allows for inelastic interactions, which can be characterized by a ‘velocity kick’ parameter. The simplest 2cDM model with cross-section $0.01\lesssim \sigma /m\lt 1\, \textrm {cm}^{2}{ \rm g}^{-1}$ and the kick velocity $V_{\mathrm{ k}}\simeq 100\, \rm {km\, s}^{-1}$ have been shown to robustly resolve the missing satellites, core-cusp, and too-big-to-fail problems in N-body cosmological simulations tested on Milky Way (MW)-like haloes of a virial mass ${\sim}5 \times 10^{11}\, {\rm M_{\odot }}$ (Papers I & II). With the aim of further constraining the parameter space available for the 2cDM model, we extend our analysis to dwarf and galaxy cluster haloes with their virial mass of ∼107−108 and ${\sim}10^{13} - 10^{14}\, {\rm M_{\odot }}$, respectively. We find that σ0/m ≳ 0.1 cm2g−1 is preferentially disfavoured for both dwarfs and galaxy cluster haloes in comparison with observations, while σ0/m = 0.001 cm2g−1 causes little perceptible difference from that of the CDM counterpart for most of the cross-section’s velocity dependence studied in this work. Our main result is that within the reasonable set of parameters, the 2cDM model can successfully explain the observational trends seen in dwarf galaxy and galaxy cluster haloes, and the model leaves us an open window for other possible alternative DM models. 
    more » « less
  3. ABSTRACT

    Galactic nuclei are potential hosts for intermediate-mass black holes (IMBHs), whose gravitational field can affect the motion of stars and compact objects. The absence of observable perturbations in our own Galactic Centre has resulted in a few constraints on the mass and orbit of a putative IMBH. Here, we show that the Laser Interferometer Space Antenna (LISA) can further constrain these parameters if the IMBH forms a binary with a compact remnant (a white dwarf, a neutron star, or a stellar-mass black hole), as the gravitational-wave signal from the binary will exhibit Doppler-shift variations as it orbits around Sgr A*. We argue that this method is the most effective for IMBHs with masses $10^3\, \mathrm{ M}_\odot \lesssim M_{\rm IMBH}\lesssim 10^5\, \mathrm{ M}_\odot$ and distances of 0.1–2 mpc with respect to the supermassive black hole, a region of the parameter space partially unconstrained by other methods. We show that in this region the Doppler shift is most likely measurable whenever the binary is detected in the LISA band, and it can help constrain the mass and orbit of a putative IMBH in the centre of our Galaxy. We also discuss possible ways for an IMBH to form a binary in the Galactic Centre, showing that gravitational-wave captures of stellar-mass black holes and neutron stars are the most efficient channel.

     
    more » « less
  4. ABSTRACT We present and study a large suite of high-resolution cosmological zoom-in simulations, using the FIRE-2 treatment of mechanical and radiative feedback from massive stars, together with explicit treatment of magnetic fields, anisotropic conduction and viscosity (accounting for saturation and limitation by plasma instabilities at high β), and cosmic rays (CRs) injected in supernovae shocks (including anisotropic diffusion, streaming, adiabatic, hadronic and Coulomb losses). We survey systems from ultrafaint dwarf ($M_{\ast }\sim 10^{4}\, \mathrm{M}_{\odot }$, $M_{\rm halo}\sim 10^{9}\, \mathrm{M}_{\odot }$) through Milky Way/Local Group (MW/LG) masses, systematically vary uncertain CR parameters (e.g. the diffusion coefficient κ and streaming velocity), and study a broad ensemble of galaxy properties [masses, star formation (SF) histories, mass profiles, phase structure, morphologies, etc.]. We confirm previous conclusions that magnetic fields, conduction, and viscosity on resolved ($\gtrsim 1\,$ pc) scales have only small effects on bulk galaxy properties. CRs have relatively weak effects on all galaxy properties studied in dwarfs ($M_{\ast } \ll 10^{10}\, \mathrm{M}_{\odot }$, $M_{\rm halo} \lesssim 10^{11}\, \mathrm{M}_{\odot }$), or at high redshifts (z ≳ 1–2), for any physically reasonable parameters. However, at higher masses ($M_{\rm halo} \gtrsim 10^{11}\, \mathrm{M}_{\odot }$) and z ≲ 1–2, CRs can suppress SF and stellar masses by factors ∼2–4, given reasonable injection efficiencies and relatively high effective diffusion coefficients $\kappa \gtrsim 3\times 10^{29}\, {\rm cm^{2}\, s^{-1}}$. At lower κ, CRs take too long to escape dense star-forming gas and lose their energy to collisional hadronic losses, producing negligible effects on galaxies and violating empirical constraints from spallation and γ-ray emission. At much higher κ CRs escape too efficiently to have appreciable effects even in the CGM. But around $\kappa \sim 3\times 10^{29}\, {\rm cm^{2}\, s^{-1}}$, CRs escape the galaxy and build up a CR-pressure-dominated halo which maintains approximate virial equilibrium and supports relatively dense, cool (T ≪ 106 K) gas that would otherwise rain on to the galaxy. CR ‘heating’ (from collisional and streaming losses) is never dominant. 
    more » « less
  5. ABSTRACT

    Previous studies of fueling black holes in galactic nuclei have argued (on scales ${\sim}0.01{-}1000\,$pc) accretion is dynamical with inflow rates $\dot{M}\sim \eta \, M_{\rm gas}/t_{\rm dyn}$ in terms of gas mass Mgas, dynamical time tdyn, and some η. But these models generally neglected expulsion of gas by stellar feedback, or considered extremely high densities where expulsion is inefficient. Studies of star formation, however, have shown on sub-kpc scales the expulsion efficiency fwind = Mejected/Mtotal scales with the gravitational acceleration as $(1-f_{\rm wind})/f_{\rm wind}\sim \bar{a}_{\rm grav}/\langle \dot{p}/m_{\ast }\rangle \sim \Sigma _{\rm eff}/\Sigma _{\rm crit}$ where $\bar{a}_{\rm grav}\equiv G\, M_{\rm tot}(\lt r)/r^{2}$ and $\langle \dot{p}/m_{\ast }\rangle$ is the momentum injection rate from young stars. Adopting this as the simplest correction for stellar feedback, $\eta \rightarrow \eta \, (1-f_{\rm wind})$, we show this provides a more accurate description of simulations with stellar feedback at low densities. This has immediate consequences, predicting the slope and normalization of the MBH − σ and MBH − Mbulge relation, LAGN −SFR relations, and explanations for outliers in compact Es. Most strikingly, because star formation simulations show expulsion is efficient (fwind ∼ 1) below total-mass surface density $M_{\rm tot}/\pi \, r^{2}\lt \Sigma _{\rm crit}\sim 3\times 10^{9}\, \mathrm{M}_{\odot }\, {\rm kpc^{-2}}$ (where $\Sigma _{\rm crit}=\langle \dot{p}/m_{\ast }\rangle /(\pi \, G)$), BH mass is predicted to specifically trace host galaxy properties above a critical surface brightness Σcrit (B-band $\mu _{\rm B}^{\rm crit}\sim 19\, {\rm mag\, arcsec^{-2}}$). This naturally explains why BH masses preferentially reflect bulge properties or central surface densities (e.g. $\Sigma _{1\, {\rm kpc}}$), not ‘total’ galaxy properties.

     
    more » « less