Galactic and extragalactic probe of dark matter with LISA’s binary black holes near their galactic center

Abstract

The upcoming LISA mission will be able to detect gravitational waves from galactic and extragalactic compact binaries. Here, we report on LISA’s capability to probe dark matter around these binaries if the latter constitute black holes. By analyzing the variation in the chirp mass of the binary, we show that depending on the black hole masses, LISA should be able to probe their surrounding dark matter to a luminosity distance of $\approx 1$ Gpc if such binaries are observed within the inner $\approx 10$ pc of their galactic center for particle-like dark matter or near the galactic solitonic core for wave-like dark matter. However, for the latter, the density of dark matter near the galactic center must be higher than predicted from dark matter only simulations. Even if a null result is recorded during the course of observation of well-localized binaries, one can rule out certain parameter spaces of dark matter as being the dominant contributor to the matter budget of the Universe.

I Introduction

[lines=3]With the emergence of gravitational wave astronomy, the use of gravitational radiation as a probe of dark matter has gained prominence, with black holes being one of the prime targets [1, 2, 3, 4, 5, 6, 7]. Like all astrophysical bodies, black holes reside within a dark matter surrounding and would slowly accumulate this matter over time. However, in contrast to dark matter accumulation in celestial bodies, such as planets [8, 9, 10], stars [11, 12, 13, 14, 15, 16, 17, 18, 19], white dwarfs [20, 21, 22, 23], neutron stars [24, 25, 26, 27, 28, 29, 30, 31, 32] which have to rely on as of yet uncertain scattering processes, the presence of event horizon makes black holes unique sinks for their surrounding material. This renders the process of dark matter capture by black holes immune to the uncertain nature of dark matter’s interactions with standard-model particles.

A direct consequence of a steady dark matter accumulation is a continuous increase in the black hole’s mass. Thus, if the black hole is part of a binary system, this mass growth leads to a faster orbital inspiral driven by the emission of gravitational radiation compared to the case when the binary lives in a vacuum. Within the first order of post-Newtonian expansion, the gravitational waveform frequency evolution depends only on the binary’s chirp mass $\mathcal{M}$. To this end, by accurately measuring both the gravitational wave frequency and its rate of change — see Eq. (22) later — one can calculate the change in $\mathcal{M}$. However, detecting such subtle variations requires a prolonged period of observation.

Current ground-based gravitational wave detectors operating in the $\approx 10-10^{3}$ Hz frequency range observe stellar-mass binaries only for a few seconds before merger [33, 34]. In contrast, the upcoming space-based LISA mission will open a new window in the 0.1-100 mHz frequency regime [35], with planned observation times of 4-10 years, offering an excellent resolution into mass variation. Over this extended duration, LISA will be capable of observing binaries across a mass range spanning stellar to supermassive black holes (SMBHs). Moreover, for strain amplitudes with large signal-to-noise ratios (SNRs), many gravitational wave-emitting systems will be individually resolvable [36]. As a consequence, stellar-mass binary black holes (BBHs) present within the Milky Way will provide us with the opportunity to probe their surroundings. On the other hand, for extragalactic distances, intermediate-mass black holes (IMBHs) with masses $10^{2}-10^{4}M_{\odot}$ become particularly important.

To this end, in this work, we investigate the impact of a dark matter environment on the evolution of BBHs that radiate within the LISA frequency band and are individually resolvable (as a binary). We show that certain BBHs that reside near the center of their galaxy will strongly interact with their surrounding dark matter halo. Depending on the distance to the source, Depending on the distance to the source, for those with gravitational wave frequencies $f\gtrsim 1$ mHz, changes in their chirp mass can be inferred from the observed frequency evolution, allowing us to probe their environment. For the latter, we focus on two variants of dark matter, namely particle-like dark matter and (wave-like) scalar-field dark matter. We also discuss the possible impact of the presence of interstellar gas around the BBHs on our analysis.

It is anticipated that within its 4-10 yr of operation, LISA will detect $\mathcal{O}(10)-\mathcal{O}(100)$ BBHs distributed across the Milky Way [37]. Although it is difficult to make such estimates for the rate of IMBH encounters with their host SMBH, IMBHs will be detectable out to redshift two [38]. Depending on the masses of these black holes, we demonstrate that LISA can probe the surrounding dark matter of BBHs to a luminosity distance of $\approx 1$ Gpc provided these binaries are observed either within the inner $\approx 10$ pc of their galactic center (for particle-like dark matter) or near their galactic solitonic core (for wave-like dark matter). However, for the latter, the density of dark matter near the galactic center must be higher than predicted from dark matter only simulations. A null result during the course of observation of spatially well-localized binaries could still rule out certain parameter spaces of dark matter as being the dominant contributor to the matter budget of the Universe.

The remainder of this paper is organized as follows. In Section II, we discuss the formalism employed for investigating the capture of dark matter by black holes. In Section III, we study the dynamics of a BBH embedded in such dark matter environments and then discuss its observability in Section IV. In Section V, we present the result for galactic and extragalactic probes of dark matter, followed by a brief discussion in Section VI.

II Capture of dark matter by black holes

For the following investigation, we assume that the dark matter medium is at rest w.r.t. galactic rest frame but has a non-zero dispersion velocity. We consider two variants for dark matter, namely those that have masses $\gtrsim 1$ eV and can be treated as point particle-like and those that have masses $\lesssim 1$ eV resulting in wave-like properties that can be modeled via scalar fields. While the former has been successful on large scales within the $\Lambda$CDM paradigm, the latter is also naturally (i.e., independent of baryonic physics) able to account for the small-scale discrepancies [39] typically associated with particle-like dark matter [40, 41].

II.1 Particle-like dark matter

For a Schwarzschild black hole with mass $m$ that is moving at a supersonic velocity $v_{\infty}$ w.r.t. the asymptotic/unperturbed background medium with density $\rho_{\infty}$, matter within the impact parameter $r_{a}={2Gm}/{v_{\infty}^{2}}$ will be captured by the black hole. The resulting accretion cross-section can be estimated as $\sigma_{\rm HL}=\pi r_{a}^{2}$ which results in the Hoyle-Lyttleton mass accretion rate [42]

$$\dot{m}_{\rm HL}=\sigma_{\rm HL}\rho_{\infty}v_{\infty}=\frac{4\pi G^{2}m^{2}\rho_{\infty}}{v_{\infty}^{3}}\,.$$ (1)

Meanwhile, following Bondi [43], for a black hole that has $v_{\infty}=0$, the mass accretion rate reads

$$\dot{m}_{\rm B}=\frac{4\pi\lambda G^{2}m^{2}\rho_{\infty}}{c_{s}^{3}}\,;\quad\lambda=\frac{1}{4}\left(\frac{2}{5-3\gamma}\right)^{\frac{5-3\gamma}{2(\gamma-1)}}\,,$$ (2)

where $c_{s}=c\sqrt{\gamma\Theta_{\infty}}$ is the local sound speed and

$$\Theta_{\infty}=\frac{k_{B}T}{m_{\chi}c^{2}}\,,$$ (3)

is the dark matter dimensionless temperature. Above $T,m_{\chi},k_{B},c$ are the dark matter (effective) temperature, particle mass, Boltzmann constant, and the speed of light in vacuum, respectively.

In deriving Eq. (2), it is assumed that the background medium can be treated as a polytrope with the equation of state $P\propto\rho^{\gamma}$, where $\gamma\leq 5/3$ for a non-relativistic fluid, and that dark matter can be treated as a non-relativistic medium. A relativistic extension of the Bondi formalism was later conducted by Michel in Ref. [44] (for a review, see Appendix of Ref. [45]). However, for a medium with $\Theta_{\infty}\lesssim 10^{-5}$, the mass accretion rate can be well modelled by Eq. (2) (cf., [46]). Since the current upper bound on the dark matter’s dimensionless temperature is expected to be $\Theta_{\infty}\approx\mathcal{O}(10^{-8})$ [47], thus one can consider Eq. (2) to hold for dark matter accretion irrespective of it being cold or warm.

Since $\dot{m}_{\rm HL}\propto v_{\infty}^{-3}$, to generate a larger accretion rate, for the current work, we are interested in the scenario where the black hole either moves at a subsonic or marginally supersonic velocity. To this end, using Eq. (1) and (2), a simple interpolation gives us the Bondi-Hoyle [48] mass accretion rate,

$$\dot{m}_{\rm BH}=\frac{4\pi\lambda G^{2}m^{2}\rho_{\infty}}{\left(c_{s}^{2}+v^{2}_{\infty}\right)^{3/2}}\,,$$ (4)

which holds for both the subsonic and supersonic scenarios. The above expression for the accretion rate would be weakly sensitive to the spin of the black holes since the accretion cross-section is set at a large distance from the black hole. Thus, the performed calculations would not be significantly impacted by the assumption of Schwarzschild geometry.

To determine $\rho_{\infty}$, we use an NFW profile for the dark matter halo; see Eq. (36). In the inner region, the halo can experience adiabatic compression due to the presence of the SMBH, resulting in a spike in dark matter density [49, 50]. To give conservative estimates of $\dot{m}$, we disregard the presence of such a density spike in this work.

II.2 Scalar field dark matter

If dark matter is instead modeled as a scale field, far away from the black hole, the incoming wave can be treated as a plane wave. The treatment of such long wavelength ($\lambda_{\rm DM}\gg m$) scalar field absorption by a Schwarzschild black hole is provided in [51]. The resulting mass accretion rate reads

$$\dot{m}_{\rm S}=\sigma_{S}(m,m_{\chi},v_{\infty})\rho_{\infty}v_{\infty}\,,$$ (5)

where (with $\hbar=G=c=1$)

$$\sigma_{\rm S}=\frac{16\pi G^{2}m^{2}}{v_{\infty}}\frac{\xi}{1-e^{-\xi}};\quad\xi=2\pi Gmm_{\chi}\frac{1+v^{2}_{\infty}}{v_{\infty}\sqrt{1-v^{2}_{\infty}}}\,.$$ (6)

The capture rate of such long-wavelength scalar fields is expected to be relatively weak compared to point-particle dark matter and likely to be important when the black hole is present in dense dark matter environments, such as near the center of our galaxy. A generic prediction for such ultralight fields is that they will produce a solitonic core near the galactic center (e.g., [39, 52]). In the absence of self-interaction, the soliton is purely supported via quantum pressure, and its density can be approximated as [39]

$$\rho=\rho_{c}\left(1+\tilde{\lambda}\left(\frac{r}{r_{c}}\right)^{2}\right)^{-8}\,,$$ (7)

where $\tilde{\lambda}=2^{1/8}-1$. Additionally, the core’s central density can be written as [52]

$$\rho_{c}=\rho_{0}\left(\frac{Gm_{\chi}^{2}}{c^{2}\hbar^{2}}\right)^{3}M_{\rm sol}^{4}\,,`$$ (8)

where $M_{\rm sol}$ is the mass of the soliton and $\rho_{0}=0.00440$ [52]. The core radius $r_{c}$ is defined as $\rho(r_{c})=\rho_{c}/2$ and takes the form [39]

$$r_{c}=\left(\frac{\rho_{c}}{0.019{M}_{\odot}\ \text{pc}^{-3}}\right)^{-1/4}\left(\frac{m_{\chi}}{10^{-22}\text{eV}}\right)^{-1/2}\text{kpc }\,.$$ (9)

The virial velocity of the solitonic core takes the form

$$v_{\infty}\equiv v_{\mathrm{vir}}=\frac{GM_{\rm sol}m_{\chi}}{\hbar}w_{0}^{1/2}$$ (10)

with $w_{0}=0.10851$ [52]. Our target black holes will be going in a circular orbit around their galactic center with velocity $v_{c}$, with $v_{c}\ll v_{\mathrm{vir}}$ ($v_{c}$ calculated in Appendix A). Thus, we make the assumption that $v_{\infty}\equiv v_{\mathrm{vir}}$, which would otherwise be determined by $v_{c}$.

The masses of black holes considered here, along with the expression for $v_{\infty}$ in Eq. (10), suggests $\xi\ll 1$, implying that Eq. (5) can be reduced to

$$\frac{\mathrm{d}m_{\rm S}}{\mathrm{d}t}=\frac{16\pi(Gm)^{2}\rho_{\infty}}{c^{3}}\,.$$ (11)

Substituting the density of the soliton in Eq. (11) thus results in the accretion rate

$$\frac{\mathrm{d}m_{\rm S}}{\mathrm{d}t}=\frac{2.5\,M_{\odot}}{10^{3}\,\mathrm{yr}} \left(\frac{m}{10^{9}M_{\odot}}\right)^{2}\left(\frac{m_{\chi}}{10^{-22}\mathrm{eV}}\right)^{6}\left(\frac{M_{\rm sol}}{10^{10}M_{\odot}}\right)^{4}$$ (12)

$$\times\left(1+\lambda\left(\frac{r}{r_{c}}\right)^{2}\right)^{-8}\,.$$

Based on numerical simulations, $M_{\rm sol}$ can be calculated using the soliton - halo mass relation as [53]

$$M_{\rm sol}=1.25\times 10^{9}\left(\frac{M_{{\mathrm{halo}}}}{10^{12}{~{}M}_{\odot}}\right)^{1/3}\left(\frac{m_{\chi}}{10^{-22}\mathrm{~{}eV}}\right)^{-1}{~{}M}_{\odot}\,,$$ (13)

where $M_{\rm halo}$ is the mass of the dark matter halo. This implies that $M_{\rm sol}\propto m_{\chi}^{-1}$ and, as we later discuss, for Milky Way like galaxy with $M_{\rm halo}\approx 10^{12}M_{\odot}$ results in a negligible impact on $\mathcal{M}$ of the binary. To this end, for the sake of demonstration of the formalism discussed here, we also consider an alternative scenario where we fix the Milky Way’s soliton mass to $M_{\rm sol}=10^{9}M_{\odot}$ (the latter value is based on Ref. [54])).

III Dynamics of a BBH immersed in a dark matter surrounding

III.1 Mass evolution

Since LISA’s binaries will have an orbital period of $10-10^{4}$ s [35], this implies $r_{a}\gg a$, $a$ being the binary’s semi-major axis. Thus far away from the binary, the infalling point-particle dark matter will perceive the binary as a single object with mass $M=m_{1}+m_{2}$ [55, 56, 57]. This gives the total mass accretion rate on the BBH as

$$\dot{M}_{\rm BH}=\dot{m}_{1,\rm BH}+\dot{m}_{2,\rm BH}+2\delta\frac{4\pi\lambda G^{2}m_{1}m_{2}\rho_{\infty}}{\left(c_{s}^{2}+v^{2}_{\infty}\right)^{3/2}}\,,$$ (14)

where $0\leq\delta\leq 1$.

As the matter falls into the total potential of mass $M$, far away from the BBH, $\delta=1$. Then, to maintain a steady-state mass accretion rate, this captured mass would eventually have to be accreted. The mass contained in the last term in Eq. (14) would be distributed among the component masses. As the current treatment only serves as a test of LISA’s capability to detect dark matter, in the following, for simplicity, we present our results assuming $\delta=0$, which provides a floor for the expected evolution of the BBH mass.

A similar effect would manifest during the descent of wave-like dark matter on the binary’s collective gravitational potential, cf. Eq. (11). Nevertheless, as previously, we treat the infall and absorption of dark matter on each black hole separately for the latter case as well.

III.2 Orbital decay

III.2.1 Due to emission of gravitational radiation

The rate of orbital decay of a BBH is governed by the rate at which it emits gravitational radiation. For the most part of the evolution, the orbital decay of BBHs considered here can be modeled in the weak field regime. If averaged over the orbital period, the decay takes the form [58]

$$\left.\frac{da}{dt}\right|_{\rm GW}=-\frac{64}{5}\frac{G^{3}\mu M^{2}}{c^{5}a^{3}\left(1-e^{2}\right)^{7/2}}\left(1+\frac{73}{24}e^{2}+\frac{37}{96}e^{4}\right)+\mathcal{O}(\dot{m}),$$ (15)

where $\mu=m_{1}m_{2}/(m_{1}+m_{2})$ is the reduced mass of the binary and $e$ is the eccentricity. We note that due to the capture of dark matter, above the masses are time-varying and only the leading-order effect has been considered. Similarly, the eccentricity of the binary also changes as

$$\left.\frac{de}{dt}\right|_{\rm GW}=-\frac{304}{15}e\frac{G^{3}\mu M^{2}}{c^{5}a^{4}\left(1-e^{2}\right)^{5/2}}\left(1+\frac{121}{304}e^{2}\right)+\mathcal{O}(\dot{m})\,,$$ (16)

where the subscript “GW” indicates that the decay is mediated by gravitational wave emission.

III.2.2 Due to impact on orbital angular momentum

We are interested in BBHs whose orbital frequency lies in the range [$10^{-4},0.1$] Hz. At such frequencies, the orbital velocity of the binary component $v_{\rm orb}\gg v_{\infty}$. Regardless of the direction of motion of the BBH w.r.t. the background dark matter medium, averaged over one full orbit, this implies that dark matter accretion would contribute negligibly to the orbital angular momentum $L$ of the binary. Thus, we assume that outside of the loss in gravitational wave emission,

$$L=\sqrt{G\mu^{2}Mr\left(1-e^{2}\right)}$$ (17)

is a constant of the motion. For a mass-changing binary, this yields

$$\left.\frac{da}{dt}\right|_{L}=-a\left[\frac{d}{dt}\ln{(\mu^{2}M})-\frac{2e\dot{e}}{1-e^{2}}\right]\,,$$ (18)

where, if averaged over one period, following equation (25) in Ref. [59], we find that

$$\left.\frac{de}{dt}\right|_{L}=-\frac{e\dot{M}}{M}\,.$$ (19)

Thus, the net orbit decay rate can be written as $\dot{a}=\dot{a}|_{\rm GW}+\dot{a}|_{L}$ and the net eccentricity decay rate as $\dot{e}=\dot{e}|_{\rm GW}+\dot{e}|_{L}$.

III.3 Gravitational wave frequency evolution

For BBHs orbiting in a Keplerian orbit, the orbital frequency $f_{\rm orb}$ satisfies

$$f_{\rm orb}^{2}=\frac{GM}{4\pi^{2}a^{3}}\,.$$ (20)

Thus, if the binary accretes dark matter, the evolution of $f_{\rm orb}$ can be calculated as

$$\frac{df_{\rm orb}}{dt}=\frac{G^{1/2}}{4\pi}a^{-3/2}M^{-1/2}\left[\dot{M}-\frac{3M}{a}\dot{a}\right]\,,$$ (21)

where in the scenario when the BBH lives in a vacuum, we can set $\dot{M}=0$. From Eq. (21), one can then calculate the gravitational wave frequency evolution, noting that the gravitational wave frequency in the $n$th harmonic should be given by $f=nf_{\rm orb}$.

IV Observability

IV.1 An evolving chirp mass

To restate, we are interested in determining the changes in the evolution of the BBH properties, owing to the capture of surrounding dark matter. The leading order orbital evolution of a gravitational wave emitting binary is determined by its chirp mass $\mathcal{M}=\mu^{3/5}M^{2/5}$ and in the case where $\dot{M}=0$, $\mathcal{M}$ remains invariant over time. Observationally, for certain binaries, LISA may be able to measure their gravitational wave frequency $f$ as well as the frequency drift $\dot{f}$ yielding

$$\mathcal{M}=\frac{c^{3}}{G}\left(\frac{5}{96}\pi^{-8/3}f^{-11/3}\dot{f}\right)^{3/5}\,,$$ (22)

which LISA can measure with the resolution [60]

$$\sigma_{\mathcal{M}}\approx\frac{3}{5}\frac{\mathcal{M}}{\dot{f}}\sigma_{\dot{f}}\,.$$ (23)

Moreover, for circular binaries, $\dot{f}$ can be measured with the resolution

$$\sigma_{\dot{f}}\approx 0.43\left(\frac{\bar{\rho}_{2}}{10}\right)^{-1}T_{\rm obs}^{-2}\,,$$ (24)

where $\bar{\rho}_{2}$ is the SNR of the gravitational wave of the binary in the second orbital harmonic and is defined later in Eq. (29). This yields

$$\sigma_{\mathcal{M}}\approx 0.258\left(\frac{\bar{\rho}_{2}}{10}\right)^{-1}\frac{\mathcal{M}}{\dot{f}}T_{\rm obs}^{-2}\,,$$ (25)

implying that a longer observation time will lead to a smaller error in the measured value of $\mathcal{M}$. Assuming that we measure $\mathcal{M}$ for two different observation times $t_{1},t_{2}$ (with $t_{2}>t_{1}$), we can then calculate the change in the binary’s chirp mass over the observation time as $\Delta\mathcal{M}=\mathcal{M}(t_{2})-\mathcal{M}(t_{1})$, where the change in $\mathcal{M}$ is driven by the capture of dark matter. For a successful detection of this change, we require that the net measurement uncertainty remain smaller than the measured value, i.e.,

$$\Delta\mathcal{M}(t_{2})>\sigma_{\mathcal{M}(t_{1})}+\sigma_{\mathcal{M}(t_{2})}\,.$$ (26)

where $\sigma_{\mathcal{M}}(t_{1})>\sigma_{\mathcal{M}}(t_{2})$. For calculational simplicity, in the current work, we adopt a more stringent criterion, i.e.,

$$\Delta\mathcal{M}(t>t_{1})\geq 2\sigma_{\mathcal{M}(t_{1})}$$ (27)

for a successful detection.

IV.2 Strain calculation

For a given value of $f,\mathcal{M},e$ and source luminosity distance $D_{L}$, the dimensionless gravitational wave strain amplitude of the binary in the $n=2$ orbital harmonic (i.e., the most dominant harmonic at low eccentricities) is [61, 36]

$$h_{2}=\frac{8G^{5/3}\mathcal{M}^{5/3}\pi^{2/3}f^{2/3}}{5^{1/2}D_{L}c^{4}}\left[1-\frac{5}{2}e^{2}+\frac{35}{24}e^{4}+O\left(e^{6}\right)\right]\,.$$ (28)

In calculating $h_{2}$ for a given binary, we will treat $D_{L}$ to remain constant over the course of observation. Moreover, since we are not interested in a particular target, here we consider the gravitational wave SNR to be averaged over inclination, sky location, and gravitational wave polarization over an observing time $T_{\rm obs}$. This yields [36]

$$\bar{\rho_{2}}=\frac{h_{2}}{S_{n}(f)^{1/2}}\sqrt{T_{\rm obs}}\,.$$ (29)

where $S_{n}(f)$ is LISA’s noise power spectral density and is defined in Appendix B. For higher harmonics, $h_{n}\propto\sqrt{g(n,e)}/n$ and become increasingly important for $e\gtrsim 0.2$ (see figure 3 in Ref. [62], the expression for $g(n,e)$ can be also be found therein). For simplicity below, we present our calculation only for the case $e=0$ for which the binary radiates at $n=2$. As the intensity of the emitted radiation is directly proportional to $e$, so our results can be treated as the floor of LISA’s sensitivity.

V Result

Figure 1: Variation in the chirp mass $\mathcal{M}$ of a BBH with $m_{1},m_{2}=25M_{\odot}$ and $e=0$ that is located at a distance of $r_{\rm center}=1$ pc, 10 pc from the center of Milky Way. Colored lines represent $\Theta_{\infty}$ - the dimensionless temperature of particle-like dark matter. The shaded region above the dotted curves represents the region where variations in $\mathcal{M}$ of a BBH with initial gravitational wave frequency $f_{\rm in}$ (i.e., $f$ at the beginning of the observation run) are detectable by LISA.

Figure 2: Same as Fig. 1 but now for wave-like dark matter. The colors of the curves represent the mass of dark matter, and the BBH is assumed to be much closer to the galactic center. For reference, the Schwarzschild radius of Sgr A^∗ is $\approx 4\times 10^{-7}$ pc.

V.1 Galactic sources

Given the inverse dependence of $h_{2}$ on $D_{L}$, binaries within the galaxy will be the optimal target. As mentioned earlier, we assume such binaries to be in a circular orbit around the galactic center with velocity $v_{c}$, which will be the source of $v_{\infty}$ in Eq. (4). In addition to a large $M$, Eq. (4) and (11) suggests that a larger change in $\mathcal{M}$ requires $\rho_{\infty}$ to be large and $v_{c}$ to be small. Thus, the ideal location to probe dark matter using LISA would be near the galactic center.

The calculation of $v_{c}$ for black holes residing in the Milky Way’s disk is presented in Appendix A. Since the role of the galactic bulge remains uncertain in the inner region, we assume that it is the NFW halo that determines the value of $v_{c}$ near the galactic center. Nevertheless, we compare $v_{c}$ resulting from the presence of an NFW halo and the bulge independently in Fig. 5. As discussed in Section II, in contrast to particle-like dark matter, absorption of ultralight scalar fields will only be important when the BBH resides within the galactic soliton. For such a case, $v_{\rm vir}\gg v_{c}$, thus the effect on a non-zero $v_{c}$ can be ignored.

Fig. 1 shows the variation in $\mathcal{M}$ of a BBH with $m_{1},m_{2}=25M_{\odot}$ due to the presence of particle-like dark matter. The shaded area in the figure shows the region of detectability by LISA and is dependent on the emitted gravitational wave frequency. The change in $\mathcal{M}$ depends on the magnitude of $\Theta_{\infty}$ and $r_{\rm center}$ (the distance of the BBH from the galactic center) with smaller values resulting in more favorable outcomes. We note that the kink in the curve with initial frequency $f_{\rm in}=2$ mHz is an artifact and results from the interpolation carried over the LISA sensitivity data presented in Table 1 of Ref. [36].

Similar to Fig. 1, Fig. 2 shows the variation in $\mathcal{M}$ for the case where the BBH is immersed in a wave-like dark matter surrounding of $M_{\rm sol}=10^{9}M_{\odot}$ with the corresponding dark matter masses shown in the figure legend. As the mass of dark matter increases, the density of the solitonic core increases. However, this also decreases the value of $r_{c}$ in Eq. 9. Thus, we place the BBH closer to the galactic center for a larger value of $m_{\chi}$. For both Fig. 1 and 2, we find that within 4-10 yr of operation, LISA will be able to probe the surroundings of the BBHs if such systems are detected during the observational run. However, $M_{\rm sol}$ is expected to be a dynamical quantity that changes as one varies $m_{\chi}$ [53]. For the sake of demonstration, let us pick $m_{\chi}=10^{-19}$ eV as this mass results in a detectable effect in Fig. 2. Then, based on Eq. (13), the corresponding $M_{\rm sol}=1.25\times 10^{6}M_{\odot}$. This reduces the mass accretion rate in Eq. (12) by a factor of $\approx 10^{12}$ compared to the case where $M_{\rm sol}=10^{9}M_{\odot}$ and it gets worse as $m_{\chi}$ increases. As the value of $M_{\rm sol}$ has a direct impact on $\rho_{c}$ in Eq. (8), adopting a larger value of $M_{\rm sol}$ results in a larger $\rho$. This effect is akin to the impact of adiabatic compression of the soliton on dark matter density, where the compression is driven by the SMBH at the galactic center [50]. Thus, we conclude that wave-like dark matter is unlikely to cause any noticeable variation in $\mathcal{M}$ unless the soliton gets adiabatically compressed to result in larger values of $\rho$ near the galactic center.

Until now, we assumed the BBH mass ratio $q\equiv m_{2}/m_{1}=1$. To investigate the impact of a lower value of $q$, we reproduce Fig. 1 and 2 in Fig. 6 in Appendix C, assuming $q=0.1$. As anticipated, a lower mass ratio reduces the prospects of detectability, meaning for smaller $q$ values, a higher value of $M$ would be desirable to probe the surrounding dark matter.

Figure 3: Top two rows: Shaded region represents the parameter space where variations in $\mathcal{M}$ of a BBH (with total mass $M$ and $q=1$) are detectable by LISA. Colors represent $f_{\rm in}$ of the BBH at the onset of observation. The two linestyles represent the length of observation. Various regions are clipped from above as these BBHs merge before the end of observation. The BBHs are assumed to be located at $r_{\rm center}=1$ pc from their galactic center. The various values of $D_{L}$ are also shown near the lower left of each plot. Bottom two rows: Same as top two rows but representing wave-like dark matter with $r_{\rm center}=10^{-4}$ pc.

V.2 Extragalactic sources

LISA will be sensitive to IMBHs out to redshift two [38]. Thus, in this section, we apply the above-discussed methodology to larger values of $D_{L}$. For simplicity, we take the Milky Way’s dark matter profile as a prototype for all galaxies at extragalactic distances.

The resulting detectability prospects are shown in Fig. 3 and 7 (in Appendix C) for $q=1,0.1$ respectively, where we again set $M_{\rm sol}=10^{9}M_{\odot}$. For demonstration purposes, here we have chosen a subset of $f_{\rm in}$ values and assumed a fixed observation time of $T_{\rm obs}=4,10$ yr. At the end of $T_{\rm obs}$, $\Delta\mathcal{M}$ is calculated as $\Delta\mathcal{M}(T_{\rm obs})=\mathcal{M}(T_{\rm obs})-\mathcal{M}(T_{\rm obs}-0.1)$. Many BBHs merge within this timespan and are thus removed from the figure, as is evident from the clipped region at the top of the shaded area.

In principle, one could fill other regions of Fig. 3 and 7 by calculating $\Delta\mathcal{M}$ as one continuously varies $T_{\rm obs}$. These figures show that observing $\mathcal{M}$ variations allow to probe $\Theta_{\infty}\lesssim 10^{-8}$ for the case of particle dark matter and $m_{\chi}\in[10^{-20},5\times 10^{-19}]$ eV for wave-like dark matter. Such variations can be probed to a distance of $\approx 1$ Gpc. The latter depends on the value of $r_{\rm center}$, and one can imagine that the detection prospects would be more favorable for smaller values of $r_{\rm center}$.

VI Discussion

VI.1 Impact of a baryonic surrounding

Solving Eq. (4) yields the mass evolution of the black hole as

$$m(t)=\eta m_{i}\,;\quad\eta(t)=-\frac{1}{m_{i}\alpha t-m_{i}\alpha t_{i}-1};\quad\eta(t)\geq 1\,,$$ (30)

where $m_{i}$ is the mass of the black hole at some initial time $t_{i}$ and $\alpha=\frac{4\pi\lambda G^{2}\rho_{\infty}}{(c_{s}^{2}+v_{\infty}^{2})^{3/2}}$. The accretion of the interstellar medium (ISM) will be most efficient when the black hole comoves w.r.t. to its surrounding medium. Thus, for the case of baryonic matter accretion, we set $v_{\infty}=0$ in Eq. (30). Moreover, the value of $\rho_{\infty}$ for the baryonic environment will depend on the nature of ISM in the vicinity of the black hole. A favorable scenario occurs when the ISM is cold and hence contains neutral particles, typically hydrogen. Such a medium is referred to as a cold neutral medium (CNM) and contains an average particle number density $n_{H}\approx(20-50)$ cm^-3 and a temperature $T\approx 100$ K [63].

It is likely that most of the surrounding CNM would have already been accreted by the time the BBH enters the LISA frequency band. Nevertheless, it is worthwhile to investigate how the black hole mass evolution compares to the scenario when the mass growth is supported by dark matter and baryonic matter accretion, respectively. To this end, Fig. 4 shows the mass evolution of black holes with $m_{i}=100M_{\odot},\,1000M_{\odot}$ immersed in a baryonic and particle-like dark matter medium with variable properties.

Figure 4: The fractional change in a black hole’s mass with initial mass $m_{i}$ located at $r_{\rm center}=$ 1 pc, 10 pc from the galactic center. The LHS (RHS) figure shows $m_{i}=100M_{\odot}$ ($m_{i}=1000M_{\odot}$). In the legend, “DM” and “BM” stand for dark matter and baryonic matter, respectively. $\Theta_{\infty}$ represents dark matter dimensionless temperature and $n_{H}$ represents baryonic particle number density. Although we have considered only two values for $m_{i}$, a similar trend holds for other values as well.

Here, the ISM is considered to be cold with $T=100$ K. However, there is a possibility that baryonic matter may form an accretion disk around the binary (by borrowing angular momentum from the BBH’s orbit) and thus heat up the surrounding gas. This will cause the temperature to rise, and for $n_{H}\approx 1$, the baryonic matter will have a subdominant impact on the black hole once $T\gtrsim 200$ K.

To be quantitative, let us assume that only a fraction $\nu$ of the accreted mass is assimilated into the accretion disk (with the rest getting accreted directly). Then, one finds that the mass assimilation rate in the disk w.r.t. the Eddington rate should evolve as [64]

$$f=\epsilon\frac{\dot{m}}{\dot{m}_{\rm Edd}}=3.9\times 10^{-3}\epsilon\nu\left(\frac{n_{H}}{{\rm cm}^{3}}\right)\left(\frac{m}{M_{\odot}}\right)\left(\frac{100}{T}\right)^{3/2}\,.$$ (31)

For illustrative purposes, assuming $\nu=0.1$, $T=100$ K, $m\in[10M_{\odot},1000M_{\odot}]$ and the accretion efficiency $\epsilon=1/16$, yields $f\in(0.00024n_{H},0.024n_{H})$. Noting that for $f\gtrsim 0.001$, one excepts the accretion disk to radiate as a Shakura & Sunyaev accretion flow [65, 66], thus, for black holes with $m\gtrsim 100M_{\odot}$, the effect of feedback could be very strong (even for $n_{H}=1$ cm^-3), likely shutting off accretion of baryonic ISM due to a rise in its temperature.

VI.2 Implication of a null result and the issue of localization

Regardless of the presence of a baryonic environment near the black holes, LISA not detecting any $\mathcal{M}$ evolution during the course of observation would imply that the underlying dark-matter model used for performing the calculation here is not the dominant form of dark matter around the BBHs. While it can still be a subdominant component, it has to be below the threshold at which $\mathcal{M}$ evolution can be observed. Depending on the value of $M,q$, and $D_{L}$, the method proposed here can be used to perform such a check for a large parameter space of dark matter.

However, LISA alone will be unable to localize most BBHs to a subgalactic scale if observed at larger distances. Following equation 15 in Ref. [60] suggests that of the BBH mass range considered here, only the massive IMBHs within $D_{L}\lesssim 1$Mpc will have sufficient angular resolution to be localized to an $\approx 10$ pc scale. Thus, a null result will only be useful for such values of $D_{L}$. This is because, as we showed, BBHs at larger distances from their galactic center would not interact strongly with their surrounding dark matter. Thus, for them, one may not expect to see any appreciable $\mathcal{M}$ variation anyway.

There remains a possibility that if the BBH is located close to its galactic center (i.e., within $\lesssim$ 100 Schwarzschild radii of the SMBH), the SMBH-BBH system might radiate gravitational waves such that the inspiral of the BBH into the SMBH also becomes resolvable by LISA. This will result in two distinct gravitational wave signals, both coming from the same location in space but from two different sources. This can be a strong clue that the BBH signal may arise from the center of the host galaxy hosting the SMBH. In such a scenario, the null result mentioned above becomes a useful tool for much larger values of $D_{L}$.

VI.3 Conclusion

In this work, we have attempted to estimate the floor of the possible $\mathcal{M}$ variation caused by dark matter accretion onto LISA’s BBHs and whether such a variation is detectable. In particular, we set $\delta=0$ - Eq. (14), assumed circular orbits, i.e., $e=0$, and considered no adiabatic compression of the dark matter halo [49, 50], the latter leading to larger dark matter density near the galactic center (although see Section V.1).

Our analysis demonstrates that the BBH surrounding environment can be probed by measuring variation in the binary’s chirp mass over the course of LISA’s observation. As the impact of baryonic matter accretion may be subdominant for IMBHs (Section VI.1), any significant variation in their $\mathcal{M}$ value can be attributed to the presence of dark matter. These variations depend on the intrinsic properties of dark matter, allowing us to probe them up to a luminosity distance of approximately 1 Gpc, depending on the BBH’s location relative to the galactic center and the particle physics property of dark matter. For the case where the BBH is very close to the galactic SMBH, the presence of an AGN disk can contaminate the resulting gravitational wave signal. It remains to be seen if meaningful information about the surrounding dark matter can still be recovered under such a scenario.

Appendix A Galactic circular velocity

We calculate the galactic mid-plane ($z=0$) circular velocity as a function of distance $r$ from the galactic center as

$$v_{\mathrm{c}}^{2}(r)=\left.r\frac{\partial\Phi}{\partial r}\right|_{z=0}\,,$$ (32)

where the gravitational potential $\Phi$ is the sum of contributions from different components, denoted as $\Phi_{i}$. To account for $\Phi_{i}$, we include contributions from a spherical dark matter halo and a spherical Galactic bulge. When a component is described by its density $\rho_{i}$ - e.g., Eq. (36), we obtain the corresponding $\Phi_{i}$ through the Poisson equation

$$\nabla^{2}\Phi_{i}=4\pi G\rho_{i}\,.$$ (33)

In the inner region where $v_{c}$ is small, $\Phi$ could be dominated by the contribution from the bulge. The potential for the latter can be approximated as a Plummer potential and takes the form

$$\Phi_{\text{Plummer}}(r)=-\frac{GM_{\text{bulge }}}{\sqrt{r^{2}+r_{b}^{2}}}\,,$$ (34)

where $M_{\rm bulge}=1.067\times 10^{10}M_{\odot}$ and $r_{b}=0.3$ kpc is the cut-off radius [67]. The resulting velocity takes the form

$$v_{c}^{2}(r)=\frac{GM_{\rm bulge}r^{2}}{(r^{2}+r_{b}^{2})^{3/2}}$$ (35)

and is shown in Fig. 5.

Although $v_{c}$ for $r\in[5,25]$ kpc has been well measured [68], uncertainties remain regarding its value in the inner region of the galaxy. Moreover, there is no unanimous agreement on the nature of the galactic bulge (for small $r$), with some studies suggesting that it might be subdominant [69]. In such a case, it is the dark matter halo that dictates the value of $v_{c}$ in the inner region of the galaxy (which we assume to be true in the present work). To calculate the resulting value of $v_{c}$, we adopt an NFW profile [70] for the Milky Way’s dark matter halo as

$$\rho_{\mathrm{NFW}}(r)=\rho_{0}\frac{r_{s}}{r}\left(1+\frac{r}{r_{s}}\right)^{-2}\,,$$ (36)

where $\rho_{0}=0.052\,M_{\odot}$pc^-3, $r_{s}=8.1$ kpc are the normalisation constant and scale radius, respectively [71]. Eq. (36), thus yields

$$v_{c}^{2}(r)=4\pi G\rho_{0}\frac{r_{s}^{3}}{r}\left[\frac{r_{s}}{r_{s}+r}-\ln{\left(\frac{r_{s}}{r_{s}+r}\right)-1}\right]$$ (37)

and has been plotted in Fig. 5.

Figure 5: Galactic circular velocity of Milky Way under the scenario when $\Phi$ is dominated by the galactic bulge (blue) and the NFW halo (orange), respectively. The inset figure shows the region near the center in units of parsec.

Appendix B LISA sensitivity

To determine the sensitivity of the LISA towards the detection of a given binary, we follow the fitting relations provided in Ref. [36]. The LISA noise power spectral density $S_{n}$ has a contribution from two sources

$$S_{n}(f)=S_{n}^{\rm{ins}}(f)+S_{n}^{\rm{WDB}}(f)\,,$$ (38)

where the former term results from the instrument noise while the latter results due to the confusion noise from the population of unresolved galactic white dwarf binaries. Explicitly,

$$S_{n}^{\rm{ins}}(f)= A_{1}\left(P_{\rm{OMS}}+2\left[1+\cos^{2}\left(f/f_{\star}\right)\right]\frac{P_{\rm{acc}}}{(2\pi f)^{4}}\right)$$ (39)

$$\times\left(1+\frac{6}{10}\frac{f^{2}}{f_{\star}^{2}}\right)\,,$$

where

$$A_{1}=\frac{10}{3L^{2}},\,L=2.5{\rm Gm},\,f_{\star}=19.09\rm{mHz}\,,$$ (40)

$$P_{\rm{OMS}}=(1.5\times 10^{-11}{\rm~{}m})^{2}\left[1+\left(\frac{2\rm{mHz}}{f}\right)^{4}\right]\rm{Hz}^{-1}\,,$$

$$P_{\rm acc}=\left(3\times 10^{-15}\rm{~{}ms}^{-2}\right)^{2}\left[1+\left(\frac{0.4\rm{mHz}}{f}\right)^{2}\right]$$

$$\quad\quad\times{\left[1+\left(\frac{f}{8\rm{mHz}}\right)^{4}\right]\rm{Hz}^{-1}}\,.$$

Additionally,

$$S_{n}^{\rm{WDB}}=A_{2}f^{-7/3}e^{-f^{\alpha}+\beta f\sin(\kappa f)}$$ (41)

$$\quad\quad\quad\quad\quad\times\left[1+\tanh\left(\gamma\left(f_{k}-f\right)\right)\right]\,\,\rm{Hz}^{-1}\,,$$

where depending on magnitude of $T_{\rm obs}$ the parameters $A_{2},\alpha,\beta,\kappa,\gamma,f_{k}$ have been interpolated from Table 1 of Ref. [36].

Appendix C Additioanl figuers

Includes Fig. 6 and 7.

Figure 6: Same as Fig. 1 and 2 but with $q=0.1$ ($M=50M_{\odot}$).

Figure 7: Same as Fig. 3 but with $q=0.1$.