## Neutron Star Collapse From Accretion: a Probe of Massive Dark Matter Particles

Ning Liu 1, ∗ and Arvind Kumar Mishra

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00594v1-figure-4e3b9fe9-1c24-40a8-8b46-32ea3ff512f4.jpg)

1 Department of Physics and Institute of Theoretical Physics, Nanjing Normal University, Nanjing, 230021, China (Dated: August 2, 2024)

We explore the multi-scatter capturing of the massive dark matter (DM) particle inside the neutron star via a momentum-dependent dark matter-nucleon scattering cross-section. We find that the capturing enhanced for the positive velocity and momentum transfer dependent DM-nucleon scattering in comparison with the constant cross-section case. Further, a large capture of the DM particles can be thermalized and lead to black hole formation and, therefore, destroy the neutron star. Using the observation of the old neutron star in the DM-dominated region, we obtain strong constraints on massive DM parameters.

∗ liuning@njnu.edu.cn

† arvindm215@gmail.com

## I. INTRODUCTION

The various observations suggest that dark matter (DM) exists from its gravitational interaction with the baryons [1]). The ongoing DM searches experiment assumes some form of the DM-baryon interactions; however, irrespective of current technological progress, the micro-physical properties of the DM particles remain an enigma [2-4]. The direct DM-detection experiments are sensitive for O (1 -10 ) GeV DM masses. 4 However, if the DM particles are massive, i.e., m χ &gt; 10 6 GeV, or lighter, m χ &lt; GeV, then during DM-nucleon collision, the energy transfer to the nucleus becomes small; therefore the detection experiment becomes insensitive to these masses.

In this work, we aim to explore the constraints on the massive dark matter, m &gt; χ 10 6 GeV, from the neutron star observations. The neutron star provides a lucrative place to study the DM-nucleon interaction, thanks to their strong gravity and density, see a recent review [5]. NS situated in the DM-dominated regions can be used to investigate the constraints in two ways: First, after capturing, DM particles become thermalized in the NS (via DM-nucleon scattering or DM-annihilation to the visible particles), which leads to heating of the NS and, therefore, observation of isolated of NS with high temperature can constraint the DM parameters [6-16]. Second, after thermalization, the DM particles can form a self-gravitating core, which can collapse, form a black hole, and ultimately destroy the NS; therefore, observation of the old neutron star in a DM-dominated region (central region of galaxies) can used to probe the DM microphysics [17-31]. However, most of the previous work has been done using a constant cross-section, without form factor, and the rest frame of the DM halo. Furthermore, incorporating the above effects into account, we obtained a strong constraint on the low-mass DM particles from the survival of NS [32].

The constraint on the massive DM has been explored from heating for both a single and multiscatter capture formalism [33, 34]. In previous works, [21-23], the single scatter capture formalism has been applied to put constraints on the massive DM particles from the survival of NS, specifically using constant DM-nucleon cross-section without form factor. However, for actual estimation, one needs to consider the multiscatter capture formalism for massive DM particles [35-37]. For other methods to explore massive DM particles, please see Refs [38-42].

In this work, we explore the constraint on the massive DM particles from the neutron star survival. To explore the capture of massive DM particles on the neutron star, we adopt from the multiscatter capture approach. Here, we focus on the velocity-dependent and momentum-transfer dependent DM-nucleon cross-section, which, to our knowledge, was not explored for massive DM particles in the context of black hole formation (see, for low mass DM constraint in our previous work [32]). In addition, we have also considered the general relativistic effect, the finite size effect of the nucleon, and the relative motion of the neutron star and galaxy frame in the capture rate.

Further, considering the power law parametrization of momentum and velocity dependent DM-nucleon cross-section, we found that the capture rate is enhanced with a positive power of velocity and momentum. We obtained that pulsar data provide the strongest constraint on DM parameter space (both for bosonic and fermionic DM case) for the model σ ∝ v 2 rel , where v rel is the relative velocity of DM and nucleon. Using this model, pulsar observation strongly constrain bosonic DM mass range, 10 6 GeV ≤ m χ ≤ 6 × 10 8 GeV, and fermionic DM masses, 10 6 GeV ≤ m χ ≤ 10 10 GeV. We also report that although the DM accretion is higher for σ ∝ v 4 rel model, a large thermalization time puts weak constraints on the DM parameters. Interestingly, a cross-section model with a large momentum dependency, i.e., σ ∝ q 4 ( q represent the momentum transfer), provides no constraint on the DM parameters.

The arrangement of our work is as follows: In section II, we present a brief discussion on the DM particle capture rate using a multiscatter capture process. In Section III, we discuss the black hole formation from DM accretion. Further, we present DM-Nucleon cross-section and estimate the captured particles IV. In Section V, we present and discuss our main results. Finally, in Section VI, we conclude our work.

## II. MULTISCATTER CAPTURE OF DM PARTICLES ON NEUTRON STAR

We assume that the neutron star is situated in a dark matter-rich environment and explore the capture of the massive DM particles to the neutron star. In the case of massive DM particles, an energy transfer via a single interaction between the DM particles and nucleons will be insufficient to capture the DM particles. Thus, a multiple scattering will be required [35-37]. Inside the NS, the number of times by which DM scatter with the nucleons can be obtained as the ratio of the typical diameter of the NS, 2 R NS with the mean free path of the DM particles, λ MFP , given by N ∼ 2 R NS /λ MFP , where R NS is neutron star radius. Assuming λ MFP ∼ 1 / n ( N σ ), one can obtain, N ∼ 2 n N σR NS , where σ is DM-nucleon scattering cross-section, and n N is the nucleon density [35].

Further, the probability of DM particles, p N ( τ ) participating in N collisions with an optical depth, τ can be defined via Poisson distribution ( τ, N) and given by [35]

$$p _ { N } ( \tau ) = 2 \int _ { 0 } ^ { 1 } y e ^ { - y \tau } \frac { ( y \tau ) ^ { N } } { N! } d y$$

where y is related to the incidence angle subtended by DM particles while entering the star. The optical depth is defined as τ = 3 σ 2 σ sat , where σ is the DM-nucleon cross-section and σ sat is the saturation cross-section allowed during the scattering between the DM-nucleon interaction. In the context of NS, a saturation limit is dictated via geometrical cross-section, given by σ geo = πR 2 NS /N B where, R NS is the radius of NS, and N B is total number of the baryons, respectively. For NS, considering, R NS = 10 6 km and . N B = 1 7 . × 10 57 , we obtain σ sat = 2 × 10 -45 cm , which we 2 assumed to remains a constant with DM mass 1 .

In the multiscatter capture process, the capture rate for N scatter, C N is given by [35] 2

$$C _ { N } = \pi R _ { N S } ^ { 2 } p _ { N } ( \tau ) \int _ { v _ { e s c } } ^ { \infty } \frac { f ( u ) } { u ^ { 2 } } \ w ^ { 3 } g _ { N } ( w ) \ d w \.$$

where g N ( w ) corresponds to the probability of DM particles reducing their velocity lower than the escape velocity, v esc after N collisions, and f ( u ) is the DM velocity distribution function. In Eq. (2), u is the DM velocity at very large distance, and w = √ u 2 + v 2 esc . Therefore, the total capture rate after the N collision is obtained as

$$C _ { \text{total} } = \sum _ { N = 1 } ^ { \infty } C _ { N } \.$$

However, summation can be cut off at some value N max corresponding to condition, p N max ( τ ) ≈ 0. In our calculation, we approximate, g N ( w ) = Θ [ v esc (1 -⟨ z i ⟩ β + 2 ) -N/ 2 -w ] , where, β + = 4 m m N χ ( m N + m χ ) 2 in which, m χ and m N , represents the mass of the DM particles and nucleons, respectively [35]. Here we assume, ⟨ z i ⟩ = 1 2, [35] however for more / general estimation of g N ( w ), see Ref. [44]. Furthermore, the DM particle velocity distribution function is considered to be Maxwellian, given in the neutron star frame as [45]

$$f _ { 0 } ( u ) = \left ( \frac { 3 } { 2 } \right ) ^ { \frac { 3 } { 2 } } \frac { 4 } { \sqrt { \pi } } \frac { \rho _ { \chi } } { m _ { \chi } } \frac { u ^ { 2 } } { u _ { 0 } ^ { 3 } } e ^ { - \frac { 3 } { 2 } \frac { u _ { 0 } ^ { 2 } } { u _ { 0 } ^ { 2 } } }$$

where ρ χ , and u 0 represent energy density, mass, and average velocity of the DM particles, respectively.

Furthermore, as on the surface of the neutron star, the DM velocity is high, so during the interaction, DM imparts large energy to the nucleons, see Ref.[46]. Therefore, in this analysis, we consider the finite-size effect of the nucleon during the DM-Nucleon interaction. Here, we adopt a conservative approach and assume the dipole form factor for DM-nucleon interaction (considered in the next subsection), given by [32, 47]

$$F ( q ^ { 2 } ) = \frac { \Lambda ^ { 4 } } { ( q ^ { 2 } + \Lambda ^ { 2 } ) ^ { 2 } }$$

where q ≃ √ 2 ( m m χ N m χ + m N ) v esc is momentum transfer during the DM-nucleon collision and m N is the nucleon mass. Further, we assume Λ ≃ 0 25 GeV throughout our analysis [32]. Furthermore, in our computation, we use the general . relativistic (GR) correction, which modifies the capture rate and escapes velocity. So, using the GR correction, the capture rate and escape velocity is modifies via [17, 35]

$$C _ { N } \to \frac { C _ { N } } { ( 1 - v _ { \text{esc} } ^ { 2 } ) } \quad \text{and} \quad v _ { \text{esc} } \to \sqrt { 2 \left [ 1 - ( 1 - v _ { \text{esc} } ^ { 2 } ) ^ { \frac { 1 } { 2 } } \right ] }.$$

## III. BLACK HOLE FORMATION FROM DM ACCRETION

After the accretion, the DM particles start interacting with the neutron star material and may get thermalized. In the context of the multi-capture process, the interaction of the DM particles with the neutron star is complex because, at each step, the DM particle transfers less momentum, and hence, the estimation of the thermalization time is complex. Now, we consider the simplified approach and assume the form of the thermalization time scale as [23]

$$t _ { t h } \simeq 1 0 7 0 0 \frac { \mu } { ( 1 + \mu ) ^ { 2 } } \left ( \frac { 1 0 ^ { 5 } K } { T _ { N S } } \right ) ^ { 2 } \left ( \frac { 1 0 ^ { - 4 5 } \text{cm} ^ { 2 } } { \sigma } \right ) \quad \text{years.}$$

1 Although we have assumed the saturation limit to be a constant quantity throughout the DM mass range, in general, it depends on the DM mass, see Refs. [5, 16, 43]. Therefore our estimation are somewhat conservative.

2 The derivation has been done is the assumption of multi-capture with v esc ( r ) ≃ v esc ( R NS ) and number density of the baryons are constant inside the neutron star, n B ( r ) = n B see Ref. [35]

where, µ = m /m χ N , and σ is DM-nucleon scattering cross-section. When the thermalization time scale is smaller than the lifetime of the NS ( t th &lt; t NS ), DM particles will be thermalized inside the NS and lie inside the thermal radius given by

$$r _ { \text{th} } = \left ( \frac { 9 T _ { \text{NS} } } { 8 \pi G _ { \text{N} } \rho _ { \text{B} } m _ { \chi } } \right ) ^ { \frac { 1 } { 2 } } \.$$

where G N is the Newton's gravitational constant, and ρ B is the baryon energy density. For a given density and temperature of the NS, r th ∝ √ 1 m χ , which implies that the thermal radius is smaller for massive DM particles and vice versa. Further, DM particles start self-gravitating whenever

$$\rho _ { \chi } > \rho _ { B } \text{ \ for \ } r \leq r _ { \text{th} }.$$

Therefore, with the help of the above condition, we define the DM particles dictated via N self , necessary for the self-gravitation of the DM particles inside the NS as

$$N _ { \text{self} } \simeq 4. 8 \times 1 0 ^ { 4 1 } \left ( \frac { 1. 4 \times 1 0 ^ { 1 4 } \text{gm/cm} ^ { 3 } } { \rho _ { B } } \right ) ^ { \frac { 1 } { 2 } } \left ( \frac { 1 0 0 \text{GeV} } { m _ { \chi } } \right ) ^ { \frac { 5 } { 2 } } \left ( \frac { T _ { \text{NS} } } { 1 0 ^ { 5 } \text{K} } \right ) ^ { \frac { 3 } { 2 } } \.$$

So the condition for DM self-gravitation suggest, N χ ≥ N Self .

Furthermore, when the number of DM particles becomes higher than those defined by the Chandrashekhar limit, then self-gravitating DM particles collapse and form a black hole. In the case of fermionic dark matter (FDM) particles, gravitational collapse requires gravity to surpass the fermi pressure. For non-interacting fermionic DM particles, the Chandrashekhar limit is obtained as [21]

$$N _ { \text{Chandra} } ^ { \text{Permissions} } \simeq 1. 8 \times 1 0 ^ { 5 1 } \left ( \frac { 1 0 0 \text{GeV} } { m _ { \chi } } \right ) ^ { 3 } \,.$$

From the above equation, it is important to emphasize that the massive DM particles form low-mass BH and vice versa. In this case, the BH formation occurs when DM particles start self-gravitating, i.e., N χ ≥ N Self , and captured DM particles are higher than the Chandrashekhar limit,i.e., N χ &gt; N Fermions Chandra . The mass of formed BH is obtained as M BH ∼ m N χ Fermions Chandra , and given by

$$M _ { \text{BH} } \simeq 1. 8 \times 1 0 ^ { 5 3 } \left ( \frac { 1 0 0 \text{GeV} } { m _ { \chi } } \right ) ^ { 2 } \text{GeV}.$$

From the above equation, it is important to emphasize that the massive DM particles form low-mass BH and vice versa. Furthermore, to form a BH in the bosonic dark matter (BDM) particle case, gravity must overcome pressure originating from the uncertainty relation. For non-interacting bosons, the Chandrashekhar limit is given by [21]

$$N _ { \text{Chandra} } ^ { \text{Bosons} } \simeq 1. 5 \times 1 0 ^ { 3 4 } \left ( \frac { 1 0 0 \text{GeV} } { m _ { \chi } } \right ) ^ { 2 } \,.$$

In bosonic DM case, BH forms as soon as DM starts self-gravitating,

$$N _ { \chi } \geq N _ { \text{Self} }, \text{ for } m _ { \chi } \leq 1 0 ^ { 2 1 } \left ( \frac { \rho _ { B } } { 1. 4 \times 1 0 ^ { 1 4 } \text{gm} / \text{cm} ^ { 3 } } \right ) \left ( \frac { 1 0 ^ { 5 } \text{K} } { T _ { \text{NS} } } \right ) ^ { 3 } \text{ GeV } \.$$

In this case, the mass of formed BH is obtained as M BH ∼ m N χ Self , and explicitly given by

$$M _ { \text{BH} } \simeq 4. 8 \times 1 0 ^ { 4 3 } \left ( \frac { 1. 4 \times 1 0 ^ { 1 4 } \text{gm} / \text{cm} ^ { 3 } } { \rho _ { B } } \right ) ^ { \frac { 1 } { 2 } } \left ( \frac { 1 0 0 \text{GeV} } { m _ { \text{X} } } \right ) ^ { \frac { 1 } { 2 } } \left ( \frac { T _ { \text{NS} } } { 1 0 ^ { 5 } \text{K} } \right ) ^ { \frac { 3 } { 2 } } \text{ GeV}.$$

From the above equation, we find that the massive DM particles and low neutron star temperature form low-mass BH and vice versa.

After its formation, BH starts growing via accreting the surrounding material and also reducing its mass via Hawking evaporation. The BH mass evolves as [21]

$$\frac { d M _ { \text{BH} } } { d t } \simeq 4 \pi \lambda _ { s } \left ( \frac { G _ { N } M _ { \text{BH} } } { v _ { s } ^ { 2 } } \right ) ^ { 2 } \rho _ { \text{B} } v _ { s } + \left ( \frac { d M _ { \text{BH} } } { d t } \right ) _ { \text{DM} } - \frac { 1 } { 1 5 6 0 \pi G _ { N } ^ { 2 } M _ { \text{BH} } ^ { 2 } } \$$

where λ s and v s are the accretion eigenvalue and sound speed. In the r.h.s, the first and second term corresponds to the accretion rate obtained from baryon (Bondi-Hoyle accretion rate) and DM particles, and the last term represents Hawking evaporation. Assuming the DM accretion contribution is small, Eq. (16) suggests an initial critical BH mass given as M crit BH = 1 2 . × 10 37 GeV. Whenever formed BH with initial mass, M ini BH is larger than the critical mass, i.e., M ini BH &gt; M crit BH , BH will grow via accretion and destroy the host NS. However, in the opposite case, i.e., M ini BH &lt; M crit BH , black hole will evaporate and decrease its mass with the time.

Furthermore, the value of the critical mass of the formed BH can be applied to estimate an upper limit on the mass of DM particles, which can be probed from the BH formation in the core of NS. The initial formed BH mass formed by fermionic DM particles, M BH ∼ m N χ Fermions Chandra , and bosonic DM particles via M BH ∼ m N χ Self . Therefore assuming that BH grows via accreting the NS materials, we obtain the m &lt; . χ 2 5 × 10 8 GeV (for T NS = 10 7 Kelvin) for BDM particles, and m χ &lt; 10 10 GeV for the FDM particles. We emphasize that these constraints are obtained using the assumption that formed BH is growing via accreting the NS material and not considered DM accretion. Therefore, when DM accretion is taken into account (or NS temperature is changed for the BDM case), the constraint will be modified.

## IV. DM ACCRETION FROM NON-CONSTANT DM-NEUCLEON ELASTIC SCATTERING

In this work, we explore the DM-Nucleon interaction, which is dependent on relative velocity and momentum transfer during the collision [48-56], which are investigated in contest of the DM capturing on sun [51, 52, 57-60], white dwarfs [36] and neutron stars [9, 12-15]. In this work, we parameterize the relative velocity-dependent crosssection (VDCS) as a power law form, given by [32, 51]

$$\sigma _ { \chi N } ^ { \vdash } ( v _ { \text{rel} } ) = \sigma _ { \chi N } \left ( \frac { v _ { \text{rel} } } { u _ { 0 } } \right ) ^ { 2 \alpha } \simeq \sigma _ { \chi N } \left ( \frac { w } { u _ { 0 } } \right ) ^ { 2 \alpha },$$

where, v rel is the DM-nucleon relative velocity, and u 0 represent the normalized velocity. Here α depends on the DM models, and we specifically focus on the positive velocity dependence cases, i.e., α = 1 and 2. This kind of interaction , is called p-wave ( α = 1), and d-wave ( α = 2) can be originated for the cases when initial state particles have relative angular momentum of 1 or 2 units, see Refs. [50, 51],[56, 61, 62]. Here we consider the typical value, u 0 = 220 km sec / [32, 51].

Further, we parameterize the momentum-dependent cross-section (MDCS) in the power law form as[32, 51]

$$\sigma _ { \chi N } ^ { \text{md} } ( q ) = \sigma _ { \chi N } \left ( \frac { q } { q _ { 0 } } \right ) ^ { 2 \beta },$$

where q is the momentum transfer and q 0 is a normalized momentum. The momentum dependent cross-section can originate from various mechanisms such as dipole interaction between the DM and baryons, parity violating DMbaryon interaction, see Refs. [48-51, 63, 64]. Further, for calculation purposes, we assume, q 0 = 40 MeV [32, 51]).

Having been equipped with the form of the DM-baryon interaction, we will now calculate the probability and DM capture rate using the multicapture process. This can be obtained by using Eq.(1) and Eq.(2). We have estimated the multiplescatter capture rate using constant, MDCS and VDCS in Appendices A,B, and C, respectively. Below, we will briefly summarize the results:

Capture rate in constant cross-section case: For constant cross-section case, the expression of probability is obtained after using Eq.(1) as

$$P _ { N } ^ { \text{cont} } ( \tau _ { 0 } ) = \frac { 2 } { \tau _ { 0 } ^ { 2 } N! } \left [ \Gamma ( 2 + N ) - \Gamma ( 2 + N, \tau _ { 0 } ) \right ]$$

where, τ 0 = 3 σ χ N F q ( 2 ) 2 σ sat , and Γ( a, b ) represent the incomplete Gamma function. Therefore using the above probability, C N can be obtained as

$$C _ { N } ^ { \text{cont} } ( \tau _ { 0 } ) = \frac { R ^ { 2 } } { 1 - v _ { \text{esc} } ^ { 2 } } \frac { n _ { \chi } } { 6 u _ { 0 } } P _ { N } ^ { \text{cont} } ( \tau _ { 0 } ) f ( v _ { N }, \eta ) \Theta ( v _ { N } - v _ { \text{esc} } )$$

where, v N = v esc ( 1 -β + 2 ) -N 2 is the DM velocity after N scattering [35]. Further, the total capture rate after N collision is given by C cont tot ( τ 0 ) = ∑ ∞ N=1 C cont N ( τ 0 ). In this case, an analytic expression of the total capture rate is

Figure 1. The capturing probability as a function of the number of the DM-Nucleon scattering, N (Left Panel), and crosssection (Right panel) for constant DM-Nucleon scattering cross-section.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00594v1-figure-4f648920-e678-4d59-a02b-fbf120aa724d.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00594v1-figure-299b3122-3eb9-4d6b-b7e9-65433a3cb4a4.jpg)

obtained using various analytic approximations, see Refs. [43, 65, 66]. Further, to calculate the total collected DM particles, we considered DM particles to be asymmetric (neglecting DM annihilation). In this assumption, the number of collected DM particles by neutron star throughout its lifetime is given by N χ ( t NS ) = C cont tot ( τ 0 ) t NS . Here, we will not apply the analytical approximations and calculate N χ ( t NS ) numerically.

Capture rate in MDCS case: For MDCS, the optical depth is given by τ md = τ 0 ( q/q 0 ) 2 β . So, the probability after N collision is obtained as

$$P _ { N } ^ { \text{MDCS} } ( \tau _ { 0 }, q ) = \frac { 2 } { \tau _ { 0 } ^ { 2 } N! } \left ( \frac { q } { q _ { 0 } } \right ) ^ { - 4 \beta } \left [ \Gamma ( 2 + N ) - \Gamma \left ( 2 + N, \tau _ { 0 } ( q / q _ { 0 } ) ^ { 2 \beta } \right ) \right ] \,.$$

The form of C N in the MDCS case is given by

$$C _ { N } ^ { M D C S } ( \tau _ { 0 }, q ) = \frac { R ^ { 2 } } { 1 - v _ { e s c } ^ { 2 } } \frac { n _ { \chi } } { 6 u _ { 0 } } P _ { N } ^ { M D C S } ( \tau _ { 0 }, q ) f ( v _ { N }, \eta ) \Theta ( v _ { N } - v _ { e s c } )$$

Further, the total capture rate after N collision is given by C MDCS tot ( τ 0 , q ) = ∑ ∞ N=1 C MDCS N ( τ 0 , q ). Further, assuming the DM particles to be asymmetric, the total captured DM particles in the MDCS case is given by N χ ( t NS ) = C MDCS tot ( τ 0 , q ) t NS . In this case, we will calculate N χ ( t NS ) numerically.

Capture rate in VDCS case: For VDCS, the optical depth is given by τ vd = τ 0 ( w/u 0 ) 2 α . So, the probability after N collision is obtained as

$$P _ { N } ^ { V D C S } ( \tau _ { 0 }, w ) = \frac { 2 } { \tau _ { 0 } ^ { 2 } N! } \left ( \frac { w } { u _ { 0 } } \right ) ^ { - 4 \alpha } \left [ \Gamma ( 2 + N ) - \Gamma \left ( 2 + N, \tau _ { 0 } ( w / u _ { 0 } ) ^ { 2 \alpha } \right ) \right ] \,.$$

Using the above probability after simplification (see appendix C), we obtain an analytic expression for the capture rate as

$$C _ { N } ^ { \text{VDCS} } ( \tau _ { 0 }, w _ { 0 } ) = \frac { R ^ { 2 } } { 1 - v _ { \text{esc} } ^ { 2 } } \frac { n _ { \text{X} } } { 6 u _ { 0 } } P _ { N } ^ { \text{VDCS} } ( \tau _ { 0 }, w _ { 0 } ) f ( v _ { N }, \eta ) \Theta ( v _ { N } - v _ { \text{esc} } )$$

In this case, the total capture rate after N collision is obtained as C VDCS tot ( τ 0 , w 0 ) = ∑ ∞ N=1 C VDCS N ( τ 0 , w 0 ). Further, in the assumption that the DM particles are asymmetric, the total captured DM particles in this case is given by N χ ( t NS ) = C VDCS tot ( τ 0 , w 0 ) t NS . In this case, we will calculate N χ ( t NS ) numerically.

Furthermore, the estimation of total accreted DM particles using the general expression of capture rate and probability is tedious. Therefore, we try to restrict our calculation up to a certain accuracy as follows: We calculate the total sum up to N = N max , where N max ≤ τ [65]. Further, we use a simplified expression of probability in some extreme cases of optical depth. For a single scatter ( τ ≪ 1) and multiscatter ( τ ≫ 1), the probability is approximated

Figure 2. The probability of multicapture in the momentum-dependent (Left panel) and velocity-dependent (Right panel) DM-Nucleon scattering cross-section.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00594v1-figure-216bfe90-6397-41e6-9d35-d493447350c0.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00594v1-figure-0cca76ac-10ea-46b3-a174-2448c37e5fd0.jpg)

as [5]

$$P _ { N } ( \tau ) \approx \begin{cases} \frac { 2 \tau ^ { N } } { N! ( N + 2 ) } + O ( \tau ^ { N + 1 } ), & i f \ \tau \ll 1 & ( 2 5 a ) \\ \frac { 2 } { \tau ^ { 2 } } \left ( N + 1 \right ) \Theta \left ( \tau - N \right ), & i f \ \tau \gg 1 & ( 2 5 b ) \\. &. &. &. &. &. &. &. &. \end{cases}$$

Therefore, using the capture rate expression estimated above, we calculate the accreted DM particles for constant, MDCS, and VDCS cases numerically.

## V. RESULTS AND DISCUSSIONS

In this Section, we will present our results. For calculation purposes, we assume the NS parameters as M NS = 1 44 . M ⊙ , R NS = 10 6 km, . N B = 1 7 . × 10 57 , and v esc = 1 8 . × 10 5 km sec [21], otherwise specified explicitly. /

## A. Effect of DM-Nucleon interactions on the multiscatter probability

In Fig. 1, we plot the probability of the DM multiscatter capture as a function of a number of scattering, N , (left panel) and cross-section (right panel). It can be seen that the probability decreases as the number of scattering N increases and becomes zero for some value of N , say N max , i.e., p N max ( τ ) ≈ 0. Further, fixing the DM parameters, m χ = 10 6 GeV, ρ χ = 0 3GeV cm . / 3 and N max takes values 1 and 2, and 3 corresponding to σ χ N = 10 -47 cm , 2 σ χ N = 10 -46 cm , and 2 σ χ N = 10 -45 cm , respectively. 2 Also, increments in the cross-section lead to an enhanced probability and, therefore, enlarge N max . Further, the right panel of Fig 1 suggests that depending on the N , the probability becomes max at some σ crit and vanishes on low DM-Nucleon scattering cross-section, σ min such that σ min &lt; σ crit . Nevertheless, increasing the number of scattering will shift the peak probability position and σ min on a comparably higher cross-section. Furthermore, we also find that increasing the DM masses decreases the probability.

In Fig. 2, we plot the probability of the DM capture using MDCS, σ ∝ q 2 β (Left panel) and VDCS, σ ∝ v 2 α rel (Right panel) cases. From the Left panel of Fig. 2, it can be seen that the probability decreases as the number of scattering N increases and becomes zero for some value N max as discussed previously. However, β = 2 shows an interesting feature that probability increases first and decreases sharply after a large scattering. For our parametrization and fixing the other DM parameters, m ,ρ χ χ and σ χ N , N max takes values 6 , 30 and 1000 corresponding to β = 0 1 , , and 2, respectively. Further, in the VDCS case (Left panel of Fig. 2), DM capturing inside the NS requires more scattering, N max . For example, keeping DM parameters fixed, m χ = 10 6 GeV, ρ χ = 0 3GeV cm . / 3 and σ χ N = 10 -45 cm , DM capturing in 2 α = 2 case required ∼ 10 9 scattering. Further, for constant cross-section, the probability decreases promptly, leading to a smaller N max value in comparison with the MDCS and VDCS cases. It is important to emphasize that the value of N max depends on the DM parameters; therefore, changing it will modify the N max value.

Figure 3. The accreted number of the DM particles via multiscatter process. The Left and Right panels correspond to the MDCS and VDCS cases, respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00594v1-figure-533d352e-6c01-4422-9985-c6af56d3b087.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00594v1-figure-8987829d-ef87-4211-a87b-ece31f45f226.jpg)

Figure 4. The thermalization time scale as a function of DM mass.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00594v1-figure-b909290b-5758-449e-9f9d-3968f240e30d.jpg)

## B. Accreted DM particles via multi-capture process

Fig. 3 shows the accreted DM particles via multi-scatter capture for MDCS (Left panel), and VDCS (Right panel) cases. We found that for both cases, the number of the accreted DM particles, N χ , enhanced with the DM scattering, but decreased with the DM mass. For the MDCS case, N χ increases with the positive power of cross-section, β &gt; 0 and number of scattering, N . Here we see that for β = 1, N χ is same for N = 1 and N = 10. This is because for such a small DM-neutron scattering and also parameters used in the calculation (see Left panel of Fig. 3), DM capturing happens through single scattering, and therefore N χ remains constant for N &gt; 1. For the β = 2 case, DM capturing becomes possible via multiscattering, and therefore, N χ increases with N until it gets saturated at N max = 1000. Furthermore, in the VDCS case (Right panel of Fig. 3) and for α = 1, the number of captured DM particles is equal for N = 1 and N = 1000 (same reason as previously). For the strong VDCS case, α = 2, DM particles get captured via multiscattering, and therefore, N χ increases till its saturation at N max = 5000. From Fig. 3, we also point out that for positive MDCS and VDCS cases, the number of captured DM particles is higher than the constant cross-section case.

Figure 5. Constraints on the DM-nucleon scattering cross-section for constant cross-section obtained using the pulsar data J2124-3858. Here, orange, red, and sky-blue color regions are excluded from the BH formation via bosonic DM, fermionic DM, and non-thermalization of the captured DM particles. Further, olive and brown color regions represent the BH evaporation formed from BDM and FDM, respectively. The black dashed line is the saturation limit of the DM-Nucleon scattering crosssection.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00594v1-figure-7c522869-3672-443f-809c-6933a863a078.jpg)

## C. Thermalization time scale

To estimate a DM thermalization time for non-constant DM-nucleon interaction, we assume as follows: In VDCS and MDCS case, we replace σ via Eq.(17), and Eq.(18), and the relative velocity is assumed to be the thermal velocity inside the NS, i.e., v rel ≃ v th = √ 3 T NS /m χ . With this simplification, we plot the thermalization time scale for our various parameterizations in Fig. 4. We find that for constant and σ ∝ v 2 rel , σ ∝ q 2 , t th is remains constant but increases linearly with mass, i.e., t th ∝ m χ . This behavior is evident from Eq.(17), and Eq.(18). It can be seen that for a few cases, i.e. σ ∝ v 2 rel , σ ∝ v 4 rel and σ ∝ q 2 , DM quickly thermalized inside the NS, and therefore, these cases are interesting in the context of the BH formation. However, for the σ ∝ q 4 case, the thermalization scale is quite large; therefore, DM does not thermalize in the typical lifetime of the NS.

## D. Constraining massive DM from NS collapse

Now, we will constrain the DM-nucleon parameter space represented via σ χ N -m χ using the survival of the neutron star situated in the DM-dominated region. We consider the cold and old neutron stars whose distance from the galaxy center, central temperature, and age are estimated accurately. For this purpose, we consider pulsar J21243858 situated at a distance 270 pc with surface temperature, T e &lt; 4 6 . × 10 K (central temperature, 5 T c = 2 5 . × 10 K 7 3 ), t age ∼ 7 81 . × 10 9 years. As this pulsar is close to the sun, we may approximate u s ≃ u ⊙ = 230 km sec [57], and / ρ χ = 0 3GeV cm . -3 , where u ⊙ represent the sun's speed relative to galactic halo.

Fig 5 shows a constraint on the DM-nucleon scattering cross-section for a constant cross-section using the pulsar data J2124-3858. Here, orange and red colors are excluded from the BH formation via bosonic and fermionic DM particles. The olive and brown color regions represent the BH evaporation regions for bosonic and fermionic DM particles, and in these regions, the constraint is relaxed. Further, sky-blue regions are excluded via non-thermalization of the captured DM particles, i.e., t th &gt; t NS . The black dashed line is the saturation limit of the DM-Nucleon scattering cross-section, and above this line, the constraints are not valid.

From Fig 5, we see that for BDM, the number of the accreted DM particles is high enough to form the BH across all the DM mass ranges. The higher DM particles form low mass BH, so for sufficiently large DM masses,

3 The central temperature can be estimated using the surface data from the analytic formula discussed in Ref. [67] .

Figure 6. Constraints on the DM-Nucleon scattering cross-section for momentum dependent cross-section ( σ ∝ q 2 (Left panel) and σ ∝ q 4 (Right panel) ) using the pulsar data J2124-3858. The representative colors carry the same meaning as earlier.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00594v1-figure-e2a133ce-5604-47b9-a501-540345774d68.jpg)

formed BH masses become smaller than the critical mass, i.e., M ini BH &gt; M crit BH , so BH start evaporating. Due to evaporation, constraints are relaxed; hence, the collapse of the neutron star may not constrain quite massive DM particles, i.e., m χ &gt; 6 × 10 8 GeV for bosonic DM. Therefore, most of the DM parameter space constraints are relaxed from the evaporation constraint. The BH formation from BDM constrains very limited parameter spaces 10 -45 cm 2 ≤ σ χ N ≤ 4 × 10 -46 cm 2 of BDM masses 5 × 10 7 GeV ≤ m χ ≤ 5 × 10 8 GeV. Further, we find that for FDM particles, BH formation happens on high DM masses, i.e. m &gt; χ 10 9 GeV. This is expected because overcoming the degeneracy pressure of the FDM requires more mass. Furthermore, BH formation in FDM does not constrain the DM microphysics as all the parameter space is relaxed from the saturation limit and evaporation constraints.

Fig. 6 shows the constraint on the DM-nucleon scattering cross-section for the MDCS, case from the pulsar observation. The Left and Right panel corresponds to σ ∝ q 2 ( N max = 10) and σ ∝ q 4 ( N max = 1000) cases, respectively. In this case, the DM accretion is enhanced (in comparison to the constant cross-section); therefore, the pulsar data constrain the massive DM particles. For σ ∝ q 2 , the pulsar data constraint on cross-section 10 -48 cm 2 ≤ σ χ N ≤ 10 -45 cm 2 for BDM masses 10 6 GeV ≤ m χ ≤ 6 × 10 8 GeV and 2 × 10 -47 cm 2 ≤ σ χ N ≤ 10 -45 cm 2 for FDM masses 10 8 GeV ≤ m χ ≤ 10 10 GeV. Although the pulsar data also constrains more massive DM particles, those parameter spaces are relaxed by Hawking evaporation. Furthermore, for σ ∝ q 4 case, DM accretion is even higher than σ ∝ q 2 case. However, in this model, the thermalization time is higher than the lifetime of the pulsar, which makes the BH formation highly unlikely, and therefore, this model does not constrain the DM parameters.

Furthermore, we constrain the velocity dependent DM-Nucleon scattering cross-section from the pulsar observation in Fig. 7. The Left and Right panel corresponds to σ ∝ v 2 rel ( N max = 1000) and σ ∝ v 4 rel ( N max = 5000) cases, respectively. For σ ∝ v 2 rel , the captured DM particles are higher than the constant case; therefore, it constrains the DM parameter space more strongly. This model constrain 10 -52 cm 2 ≤ σ χ N ≤ 10 -50 cm 2 for BDM masses, 10 6 GeV ≤ m χ ≤ 6 × 10 8 GeV, and 10 -50 cm 2 ≤ σ χ N ≤ 2 × 10 -47 cm 2 for FDM masses, 10 6 GeV ≤ m χ ≤ 10 10 GeV. Furthermore, due to large positive power, the number of accreted DM particles is highest for model σ ∝ v 4 rel (among all models discussed previously). Therefore we expect that this model will put the tightest constraint. However, as seen in Fig. 7, the thermalization constraint strongly restricts the DM parameter space and limits the cross-section, 10 -50 cm 2 ≤ σ χ N ≤ 10 -47 cm 2 for the BDM masses, 10 6 GeV ≤ m χ ≤ 6 × 10 8 GeV, and 10 -50 cm 2 ≤ σ χ N ≤ 3 × 10 -46 cm 2 for the FDM masses, 10 6 GeV ≤ m χ ≤ 10 10 GeV.

From the above discussion, we find that the pulsar data provide the strongest constraint on DM parameter space (both for bosonic and fermionic DM case) for the model σ ∝ v 2 rel .

Figure 7. Constraints on the DM-Nucleon scattering cross-section for velocity-dependent cross-section ( σ ∝ v 2 rel (Left panel) and σ ∝ v 4 rel (Right panel) ) using the pulsar data J2124-3858. The representative colors carry the same meaning as earlier.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00594v1-figure-901224fb-6cf2-4af9-9a51-e87d7f638c96.jpg)

## VI. CONCLUSION

A neutron star situated in a DM-dominated environment can efficiently capture the dark matter particles, and for large capturing, dark matter forms a black hole and destroys the host neutron star. In this work, we constrain the massive dark matter particles from the survival of the neutron star.

To estimate the DM capture rate of massive dark matter particles m χ &gt; 10 6 GeV, we adapted the multiscatter capture formalism to a non-constant DM-nucleon scattering cross-section. Adapting a model-independent approach, we assumed a power law parametrization of the velocity dependent and momentum dependent DM-Nucleon scattering cross-section, σ vd χ N ( v rel ) = σ χ N ( v rel /u 0 ) 2 α and σ md χ N ( q ) = σ χ N ( q/q 0 ) 2 β . Further, using the form factor and incorporating the general relativistic correction, we obtain analytic expressions for the capture rate in the neutron frame. Our estimation shows that the capture rate enhanced for the positive velocity and momentum dependent cross-section, and a number of the DM-Neucleon scattering, which leads to a black hole formation in a typical lifetime of the neutron star.

Further, applying the observation of pulsar data J2124-3858, we explored the constraint on the dark matter parameter space, i.e., σ χ N -m χ for non-constant cross-section. We obtained that pulsar data provide the strongest constraint on DM parameter space (both for bosonic and fermionic DM case) for the model σ ∝ v 2 rel . This model constrain 10 -52 cm 2 ≤ σ χ N ≤ 10 -50 cm 2 for bosonic DM masses, 10 6 GeV ≤ m χ ≤ 6 × 10 8 GeV, and 10 -50 cm 2 ≤ σ χ N ≤ 2 × 10 -47 cm 2 for fermionic DM masses, 10 6 GeV ≤ m χ ≤ 10 10 GeV. Interestingly, we also found that although for model σ ∝ v 4 rel , the DM accretion is higher, but due to a large thermalization time, it put weak constraints on the DM parameters. Furthermore, a large momentum dependent cross-section model, i.e., σ ∝ q 4 provides no constraint on the DM parameters.

To the best of our knowledge, this is the first constraint on the massive DM particles from the neutron star survival using the velocity and momentum dependent cross-section (in the multiscatter capture formalism). Nevertheless, we also point out that our calculation of capture rate has been done under certain assumptions, such as non-relativistic scattering of DM-nucleon, saturation cross-section is constant with DM mass (however, varies for massive DM particle), nucleon density (also escape velocity) is constant. After inclusion, these effects will change the constraint, and we leave these interesting aspects for the future.

## VII. ACKNOWLEDGEMENT

This work is supported by the National Natural Science Foundation of China (NNSFC) under grants No. 12275134, No. 12335005, and No. 12147228. AKM would like to thank Prof. Lei Wu and Prof. Chih-Ting Lu for the useful discussions.

## Appendix A: Calculation of multiscatter capture of DM particles via constant cross-section

After having equipped with the probability expression (see Section), we now estimate the multiscatter capture rate using Eq.(2). The derivation of Eq.(2) has been done is the assumption of multi-capture with v esc ( r ) ≃ v esc ( R NS ) and number density of the baryons are constant inside the neutron star, n B ( r ) = n B see Ref. [35]. To estimate the capture rate, we apply the Maxwell distribution function in the neutron star frame, Eq. (4).

In a constant cross-section case, the probability, P cont N ( τ 0 ) is independent of the DM velocity. So, applying Eq.(2), Eq.(4), and Eq.(19) the analytic expression of C N is obtained as

$$C _ { N } ^ { \text{cont} } ( \tau ) = \frac { \pi R ^ { 2 } } { 1 - v _ { \text{esc} } ^ { 2 } } \frac { 1 } { 6 \pi u _ { 0 } } \frac { \rho _ { \chi } } { m _ { \chi } } P _ { N } ^ { \text{cont} } ( \tau _ { 0 } ) f ( v _ { N }, \eta ) \Theta ( v _ { N } - v _ { \text{esc} } )$$

where

$$f$$

$$& \text{where} \\ & f ( v _ { N }, \eta ) = \frac { 1 } { \eta } \left [ e ^ { - \frac { 3 ( c ) ^ { 2 } _ { 2 } - v _ { 2 } ^ { 2 } s _ { 2 } ) } { 2 a _ { 0 } ^ { 0 } } + \frac { 2 \sqrt { v _ { N } ^ { 2 } - v _ { 2 } ^ { 2 } s _ { 2 } } } { u _ { 0 } } + \eta ^ { 2 } \right ] \times \left [ - \sqrt { 6 } \pi u _ { 0 } \left ( - 1 + \frac { 2 \sqrt { v _ { N } ^ { 2 } } - v _ { 2 } ^ { 2 } s _ { 2 } } \eta } { u _ { 0 } } \right ) \sqrt { v _ { N } ^ { 2 } - v _ { 2 } ^ { 2 } s _ { 2 } } \\ & \left ( 1 + e ^ { \frac { \sqrt { ( v _ { N } ^ { 2 } - v _ { 2 } ^ { 2 } s _ { 2 } ) } } { u _ { 0 } } } - 2 e ^ { - \frac { 3 2 } { 2 a _ { 0 } ^ { 0 } } \left [ ( v _ { N } ^ { 2 } - v _ { 2 } ^ { 2 } s _ { 2 } ) + 2 \sqrt { ( v _ { N } ^ { 2 } - v _ { 2 } ^ { 2 } s _ { 2 } u _ { 0 } \eta } \right ] } \right ) u _ { 0 } \eta + \pi \left ( 3 v _ { 2 } ^ { 2 } + u _ { 0 } ^ { 2 } ( 1 + 3 \eta ^ { 2 } ) \right ) e ^ { \frac { 3 2 } { 2 a _ { 0 } ^ { 0 } } \left [ v _ { N } ^ { 2 } - v _ { 2 } ^ { 2 } + 2 \sqrt { v _ { N } ^ { 2 } - v _ { 2 } ^ { 2 } s _ { 2 } } \eta + u _ { 0 } ^ { 2 } \eta ^ { 2 } \right ] } \\ & \left \{ ( 2 E r f \sqrt { 3 } ^ { 2 } \frac { \eta + E r f } { e r f } \left ( \frac { \sqrt { \frac { 3 } { 2 } } ( \sqrt { v _ { N } ^ { 2 } - v _ { 2 } ^ { 2 } s _ { 2 } } - u _ { 0 } \eta ) } { u _ { 0 } } \right ) + E r f \left ( \frac { \sqrt { \frac { 3 } { 2 } } ( \sqrt { v _ { N } ^ { 2 } - v _ { 2 } ^ { 2 } } + u _ { 0 } \eta ) } { u _ { 0 } } \right ) \right \} \right ], \\ & \quad / \quad \ \ _ { \ } \ \cdot \ \, - N / 2 \quad \ \ _ { \ } \ \, - \quad \, - \quad \,.$$

where, v N = v esc ( 1 -β + 2 ) -N 2 / is the DM velocity after N scattering, and erf(z) = 2 √ 2 ∫ z 0 e -y 2 dy represent the error function. In the above equation, ρ χ is the DM energy density surrounding the NS, and N B is the total number of baryons inside the NS. We checked that in the rest frame of DM halo, i.e., η → 0, and point-like structure of nucleon (neglecting the form factor), our expression of C const χ N matches with the Ref. [68].

## Appendix B: Derivation of multiscatter capture of DM particles via MDCS

In this case, σ md χ N ∝ q 2 β , where q ≃ √ 2 ( m m χ N m χ + m N ) v esc . So, P MDCS N ( τ 0 , q ) is independent of the initial DM velocity. Using this, we obtain the analytic form of the capture rate after N scattering for the arbitrary momentum dependence cross section. The form of C N in the MDCS case is given by

$$C _ { N } ^ { \text{MDCS} } ( \tau _ { 0 }, q ) = \frac { \pi R ^ { 2 } } { 1 - v _ { \text{esc} } ^ { 2 } } \frac { 1 } { 6 \pi u _ { 0 } } \frac { \rho _ { \chi } } { m _ { \chi } } P _ { N } ^ { \text{MDCS} } ( \tau _ { 0 }, q ) f ( v _ { N }, \eta ) \Theta ( v _ { N } - v _ { \text{esc} } )$$

## Appendix C: Derivation of multiscatter capture of DM particles via VDCS

In this case, σ md χ N ∝ w 2 β , where w = √ u 2 + v 2 esc . Therefore the probability P VDCS N ( τ 0 , w ) depends on the DM velocity, see Eq. (23). Using this probability expression, no reliable analytic (or numerical) solution was obtained. Therefore, to calculate the capture rate, we follow a conservative approach and approximate w ≃ w 0 = √ u 2 0 + v 2 esc in the probability expression. After applying this approximation, we can easily integrate and obtain an analytic expression for the capture rate as

$$C _ { N } ^ { \text{VDCS} } ( \tau _ { 0 }, w _ { 0 } ) = \frac { \pi R ^ { 2 } } { 1 - v _ { \text{esc} } ^ { 2 } } \frac { 1 } { 6 \pi u _ { 0 } } \frac { \rho _ { \text{X} } } { m _ { \text{X} } } P _ { N } ^ { \text{VDCS} } ( \tau _ { 0 }, w _ { 0 } ) f ( v _ { N }, \eta ) \Theta ( v _ { N } - v _ { \text{esc} } )$$

[1] G. Bertone, D. Hooper, and J. Silk, 'Particle dark matter: Evidence, candidates and constraints,' Phys. Rept. 405 (2005) 279-390, arXiv:hep-ph/0404175 .

[2] M. Klasen, M. Pohl, and G. Sigl, 'Indirect and direct search for dark matter,' Prog. Part. Nucl. Phys. 85 (2015) 1-32, arXiv:1507.03800 [hep-ph] .

- [3] J. Liu, X. Chen, and X. Ji, 'Current status of direct dark matter detection experiments,' Nature Phys. 13 no. 3, (2017) 212-216, arXiv:1709.00688 [astro-ph.CO] .
- [4] C. P´ erez de los Heros, 'Status, Challenges and Directions in Indirect Dark Matter Searches,' Symmetry 12 no. 10, (2020) 1648, arXiv:2008.11561 [astro-ph.HE] .
- [5] J. Bramante and N. Raj, 'Dark matter in compact stars,' Phys. Rept. 1052 (2024) 1-48, arXiv:2307.14435 [hep-ph] .
- [6] A. Reisenegger, P. Jofre, R. Fernandez, and E. Kantor, 'Rotochemical Heating of Neutron Stars: Rigorous Formalism with Electrostatic Potential Perturbations,' Astrophys. J. 653 (2006) 568-572, arXiv:astro-ph/0606322 .
- [7] C. Kouvaris, 'WIMP Annihilation and Cooling of Neutron Stars,' Phys. Rev. D 77 (2008) 023006, arXiv:0708.2362 [astro-ph] .
- [8] N. Raj, P. Tanedo, and H.-B. Yu, 'Neutron stars at the dark matter direct detection frontier,' Phys. Rev. D 97 no. 4, (2018) 043006, arXiv:1707.09442 [hep-ph] .
- [9] N. F. Bell, G. Busoni, and S. Robles, 'Heating up Neutron Stars with Inelastic Dark Matter,' JCAP 09 (2018) 018, arXiv:1807.02840 [hep-ph] .
- [10] K. Yanagi, N. Nagata, and K. Hamaguchi, 'Cooling Theory Faced with Old Warm Neutron Stars: Role of Non-Equilibrium Processes with Proton and Neutron Gaps,' Mon. Not. Roy. Astron. Soc. 492 no. 4, (2020) 5508-5523, arXiv:1904.04667 [astro-ph.HE] .
- [11] K. Hamaguchi, N. Nagata, and K. Yanagi, 'Dark Matter Heating vs. Rotochemical Heating in Old Neutron Stars,' Phys. Lett. B 795 (2019) 484-489, arXiv:1905.02991 [hep-ph] .
- [12] N. F. Bell, G. Busoni, S. Robles, and M. Virgato, 'Improved Treatment of Dark Matter Capture in Neutron Stars II: Leptonic Targets,' JCAP 03 (2021) 086, arXiv:2010.13257 [hep-ph] .
- [13] R. Garani, A. Gupta, and N. Raj, 'Observing the thermalization of dark matter in neutron stars,' Phys. Rev. D 103 no. 4, (2021) 043019, arXiv:2009.10728 [hep-ph] .
- [14] A. Joglekar, N. Raj, P. Tanedo, and H.-B. Yu, 'Dark kinetic heating of neutron stars from contact interactions with relativistic targets,' Phys. Rev. D 102 no. 12, (2020) 123002, arXiv:2004.09539 [hep-ph] .
- [15] M. Fujiwara, K. Hamaguchi, N. Nagata, and J. Zheng, 'Capture of electroweak multiplet dark matter in neutron stars,' Phys. Rev. D 106 no. 5, (2022) 055031, arXiv:2204.02238 [hep-ph] .
- [16] M. Baryakhtar, J. Bramante, S. W. Li, T. Linden, and N. Raj, 'Dark Kinetic Heating of Neutron Stars and An Infrared Window On WIMPs, SIMPs, and Pure Higgsinos,' Phys. Rev. Lett. 119 no. 13, (2017) 131801, arXiv:1704.01577 [hep-ph] .
- [17] I. Goldman and S. Nussinov, 'Weakly Interacting Massive Particles and Neutron Stars,' Phys. Rev. D 40 (1989) 3221-3230.
- [18] A. Gould, B. T. Draine, R. W. Romani, and S. Nussinov, 'Neutron stars: Graveyard of charged dark matter,' Physics Letters B 238 no. 2-4, (Apr., 1990) 337-343.
- [19] C. Kouvaris and P. Tinyakov, 'Can Neutron stars constrain Dark Matter?,' Phys. Rev. D 82 (2010) 063531, arXiv:1004.0586 [astro-ph.GA] .
- [20] A. de Lavallaz and M. Fairbairn, 'Neutron Stars as Dark Matter Probes,' Phys. Rev. D 81 (2010) 123521, arXiv:1004.0629 [astro-ph.GA] .
- [21] S. D. McDermott, H.-B. Yu, and K. M. Zurek, 'Constraints on Scalar Asymmetric Dark Matter from Black Hole Formation in Neutron Stars,' Phys. Rev. D 85 (2012) 023519, arXiv:1103.5472 [hep-ph] .
- [22] J. Bramante and T. Linden, 'Detecting Dark Matter with Imploding Pulsars in the Galactic Center,' Phys. Rev. Lett. 113 no. 19, (2014) 191301, arXiv:1405.1031 [astro-ph.HE] .
- [23] R. Garani, Y. Genolini, and T. Hambye, 'New Analysis of Neutron Star Constraints on Asymmetric Dark Matter,' JCAP 05 (2019) 035, arXiv:1812.08773 [hep-ph] .
- [24] J. Bramante, K. Fukushima, J. Kumar, and E. Stopnitzky, 'Bounds on self-interacting fermion dark matter from observations of old neutron stars,' Phys. Rev. D 89 no. 1, (2014) 015010, arXiv:1310.3509 [hep-ph] .
- [25] J. Bramante, K. Fukushima, and J. Kumar, 'Constraints on bosonic dark matter from observation of old neutron stars,' Phys. Rev. D 87 no. 5, (2013) 055012, arXiv:1301.0036 [hep-ph] .
- [26] N. F. Bell, A. Melatos, and K. Petraki, 'Realistic neutron star constraints on bosonic asymmetric dark matter,' Phys. Rev. D 87 no. 12, (2013) 123507, arXiv:1301.6811 [hep-ph] .
- [27] T. G¨ uver, A. E. Erkoca, M. Hall Reno, and I. Sarcevic, 'On the capture of dark matter by neutron stars,' JCAP 05 (2014) 013, arXiv:1201.2400 [hep-ph] .
- [28] C. Kouvaris, P. Tinyakov, and M. H. G. Tytgat, 'NonPrimordial Solar Mass Black Holes,' Phys. Rev. Lett. 121 no. 22, (2018) 221102, arXiv:1804.06740 [astro-ph.HE] .
- [29] G.-L. Lin and Y.-H. Lin, 'Analysis on the black hole formations inside old neutron stars by isospin-violating dark matter with self-interaction,' JCAP 08 (2020) 022, arXiv:2004.05312 [hep-ph] .
- [30] R. Garani, D. Levkov, and P. Tinyakov, 'Solar mass black holes from neutron stars and bosonic dark matter,' Phys. Rev. D 105 no. 6, (2022) 063019, arXiv:2112.09716 [hep-ph] .
- [31] D. Singh, A. Gupta, E. Berti, S. Reddy, and B. S. Sathyaprakash, 'Constraining properties of asymmetric dark matter candidates from gravitational-wave observations,' Phys. Rev. D 107 no. 8, (2023) 083037, arXiv:2210.15739 [gr-qc] .
- [32] C.-T. Lu, A. K. Mishra, and L. Wu, 'Constraining Bosonic Dark Matter-Baryon Interactions from Neutron Star Collapse,' arXiv:2404.07187 [hep-ph] .
- [33] N. F. Bell, G. Busoni, S. Robles, and M. Virgato, 'Thermalization and Annihilation of Dark Matter in Neutron Stars,' arXiv:2312.11892 [hep-ph] .

- [34] F. Anzuini, N. F. Bell, G. Busoni, T. F. Motta, S. Robles, A. W. Thomas, and M. Virgato, 'Improved treatment of dark matter capture in neutron stars III: nucleon and exotic targets,' JCAP 11 no. 11, (2021) 056, arXiv:2108.02525 [hep-ph] .
- [35] J. Bramante, A. Delgado, and A. Martin, 'Multiscatter stellar capture of dark matter,' Phys. Rev. D 96 no. 6, (2017) 063002, arXiv:1703.04043 [hep-ph] .
- [36] H. Steigerwald, V. Marra, and S. Profumo, 'Revisiting constraints on asymmetric dark matter from collapse in white dwarf stars,' Phys. Rev. D 105 no. 8, (2022) 083507, arXiv:2203.09054 [astro-ph.CO] .
- [37] R. K. Leane and J. Smirnov, 'Dark matter capture in celestial objects: treatment across kinematic and interaction regimes,' JCAP 12 (2023) 040, arXiv:2309.00669 [hep-ph] .
- [38] LHAASO Collaboration, Z. Cao et al. , 'Constraints on Heavy Decaying Dark Matter from 570 Days of LHAASO Observations,' Phys. Rev. Lett. 129 no. 26, (2022) 261103, arXiv:2210.15989 [astro-ph.HE] .
- [39] D. Song, K. Murase, and A. Kheirandish, 'Constraining decaying very heavy dark matter from galaxy clusters with 14 year Fermi-LAT data,' JCAP 03 (2024) 024, arXiv:2308.00589 [astro-ph.HE] .
- [40] S. Das, J. A. Carpio, and K. Murase, 'Probing superheavy dark matter through lunar radio observations of ultrahigh-energy neutrinos,' arXiv:2405.06382 [hep-ph] .
- [41] LZ Collaboration, J. Aalbers et al. , 'New constraints on ultraheavy dark matter from the LZ experiment,' Phys. Rev. D 109 no. 11, (2024) 112010, arXiv:2402.08865 [hep-ex] .
- [42] S. Bhattacharya, A. L. Miller, and A. Ray, 'Continuous Gravitational Waves: A New Window to Look for Heavy Non-annihilating Dark Matter,' arXiv:2403.13886 [hep-ph] .
- [43] C. Ilie, C. Levy, and J. Diks, 'The effectiveness of exoplanets and Brown Dwarfs as sub-GeV Dark Matter detectors,' JCAP 04 (2024) 082, arXiv:2312.13979 [astro-ph.CO] .
- [44] B. Dasgupta, A. Gupta, and A. Ray, 'Dark matter capture in celestial objects: Improved treatment of multiple scattering and updated constraints from white dwarfs,' JCAP 08 (2019) 018, arXiv:1906.04204 [hep-ph] .
- [45] W. H. Press and D. N. Spergel, 'Capture by the sun of a galactic population of weakly interacting, massive particles,' Astrophys. J. 296 (Sept., 1985) 679-684.
- [46] N. F. Bell, G. Busoni, T. F. Motta, S. Robles, A. W. Thomas, and M. Virgato, 'Nucleon Structure and Strong Interactions in Dark Matter Capture in Neutron Stars,' Phys. Rev. Lett. 127 no. 11, (2021) 111803, arXiv:2012.08918 [hep-ph] .
- [47] Y. Ema, F. Sala, and R. Sato, 'Neutrino experiments probe hadrophilic light dark matter,' SciPost Phys. 10 no. 3, (2021) 072, arXiv:2011.01939 [hep-ph] .
- [48] S. Chang, A. Pierce, and N. Weiner, 'Momentum Dependent Dark Matter Scattering,' JCAP 01 (2010) 006, arXiv:0908.3192 [hep-ph] .
- [49] J. Fan, M. Reece, and L.-T. Wang, 'Non-relativistic effective theory of dark matter direct detection,' JCAP 11 (2010) 042, arXiv:1008.1591 [hep-ph] .
- [50] J. Kumar and D. Marfatia, 'Matrix element analyses of dark matter scattering and annihilation,' Phys. Rev. D 88 no. 1, (2013) 014035, arXiv:1305.1611 [hep-ph] .
- [51] A. C. Vincent, A. Serenelli, and P. Scott, 'Generalised form factor dark matter in the Sun,' JCAP 08 (2015) 040, arXiv:1504.04378 [hep-ph] .
- [52] A. C. Vincent, P. Scott, and A. Serenelli, 'Updated constraints on velocity and momentum-dependent asymmetric dark matter,' JCAP 11 (2016) 007, arXiv:1605.06502 [hep-ph] .
- [53] J. Ooba, H. Tashiro, and K. Kadota, 'Cosmological constraints on the velocity-dependent baryon-dark matter coupling,' JCAP 09 (2019) 020, arXiv:1902.00826 [astro-ph.CO] .
- [54] K. Maamari, V. Gluscevic, K. K. Boddy, E. O. Nadler, and R. H. Wechsler, 'Bounds on velocity-dependent dark matter-proton scattering from Milky Way satellite abundance,' Astrophys. J. Lett. 907 no. 2, (2021) L46, arXiv:2010.02936 [astro-ph.CO] .
- [55] M. A. Buen-Abad, R. Essig, D. McKeen, and Y.-M. Zhong, 'Cosmological constraints on dark matter interactions with ordinary matter,' Phys. Rept. 961 (2022) 1-35, arXiv:2107.12377 [astro-ph.CO] .
- [56] K. K. Boddy, G. Krnjaic, and S. Moltner, 'Investigation of CMB constraints for dark matter-helium scattering,' Phys. Rev. D 106 no. 4, (2022) 043510, arXiv:2204.04225 [astro-ph.CO] .
- [57] W.-L. Guo, Z.-L. Liang, and Y.-L. Wu, 'Direct detection and solar capture of dark matter with momentum and velocity dependent elastic scattering,' Nucl. Phys. B 878 (2014) 295-308, arXiv:1305.0912 [hep-ph] .
- [58] A. C. Vincent and P. Scott, 'Thermal conduction by dark matter with velocity and momentum-dependent cross-sections,' JCAP 04 (2014) 019, arXiv:1311.2074 [astro-ph.CO] .
- [59] A. C. Vincent, P. Scott, and A. Serenelli, 'Possible Indication of Momentum-Dependent Asymmetric Dark Matter in the Sun,' Phys. Rev. Lett. 114 no. 8, (2015) 081302, arXiv:1411.6626 [hep-ph] .

[60]

G. Busoni, A. De Simone, P. Scott, and A. C. Vincent, 'Evaporation and scattering of momentum- and velocity-dependent dark matter in the Sun,'

JCAP

10

(2017) 037, arXiv:1703.07784 [hep-ph]

.

- [61] K. Sigurdson, M. Doran, A. Kurylov, R. R. Caldwell, and M. Kamionkowski, 'Dark-matter electric and magnetic dipole moments,' Phys. Rev. D 70 (2004) 083501, arXiv:astro-ph/0406355 . [Erratum: Phys.Rev.D 73, 089903 (2006)].
- [62] C. Dvorkin, K. Blum, and M. Kamionkowski, 'Constraining Dark Matter-Baryon Scattering with Linear Cosmology,' Phys. Rev. D 89 no. 2, (2014) 023519, arXiv:1311.2937 [astro-ph.CO] .
- [63] B. Feldstein, A. L. Fitzpatrick, and E. Katz, 'Form Factor Dark Matter,' JCAP 01 (2010) 020, arXiv:0908.2991 [hep-ph] .

- [64] S. Chang, N. Weiner, and I. Yavin, 'Magnetic Inelastic Dark Matter,' Phys. Rev. D 82 (2010) 125011, arXiv:1007.4200 [hep-ph] .
- [65] C. Ilie, J. Pilawa, and S. Zhang, 'Comment on 'Multiscatter stellar capture of dark matter',' Phys. Rev. D 102 no. 4, (2020) 048301, arXiv:2005.05946 [astro-ph.CO] .
- [66] C. Ilie, 'Closed-form expressions for multiscatter Dark Matter capture rates,' arXiv:2402.07713 [astro-ph.CO] .
- [67] E. H. Gudmundsson, C. J. Pethick, and R. I. Epstein, 'Neutron star envelopes,' Astrophys. J. 259 (1982) L19-L23.
- [68] A. Gould, 'Resonant Enhancements in Weakly Interacting Massive Particle Capture by the Earth,' Astrophys. J. 321 (Oct., 1987) 571.
Ning Liu Arvind Kumar Mishra
Neutron Star Collapse From Accretion: a Probe of Massive Dark Matter Particles
We explore the multi-scatter capturing of the massive dark matter (DM)
particle inside the neutron star via a momentum-dependent dark matter-nucleon
scattering cross-section. We find that the capturing enhanced for the positive
velocity and momentum transfer dependent DM-nucleon scattering in comparison
with the constant cross-section case. Further, a large capture of the DM
particles can be thermalized and lead to black hole formation and, therefore,
destroy the neutron star. Using the observation of the old neutron star in the
DM-dominated region, we obtain strong constraints on massive DM parameters.