The Dark Side of Neutron Stars

We review severe constraints on asymmetric bosonic dark matter based on observations of old neutron stars. Under certain conditions, dark matter particles in the form of asymmetric bosonic WIMPs can be effectively trapped onto nearby neutron stars, where they can rapidly thermalize and concentrate in the core of the star. If some conditions are met, the WIMP population can collapse gravitationally and form a black hole that can eventually destroy the star. Based on the existence of old nearby neutron stars, we can exclude certain classes of dark matter candidates.


I. INTRODUCTION
Compact stars such as neutron stars and white dwarfs can lead in general to two types of constraints regarding dark matter candidates. The first one has to do with annihilating dark matter that changes the thermal evolution of the star. Annihilation of Weakly Interacting Massive Particles (WIMPs) that are trapped inside compact stars, can lead to the production of significant amount of heat that can change the temperature of old stars [1][2][3][4]. Such a phenomenon can be in principle contrasted to observations. The second type of constraints is related to asymmetric dark matter [5][6][7][8][9][10][11][12]. Asymmetric dark matter is an attractive alternative to thermally produced dark matter especially due to the intriguing possibility of relating its asymmetry to the baryonic one. For recent reviews see [13,14]. Due to the asymmetry, WIMP annihilation is not significant in this case. If a certain amount of WIMPs is trapped inside the star, the WIMPs can quite rapidly thermalize and concentrate within a tiny radius in the core of the star. If the WIMP population grows significantly, WIMPs might become self-gravitating and they might collapse forming a mini black hole. Under certain conditions, the black hole might consume the rest of the star, thus leading to the ultimate destruction of the star. However, very old (older than a few billion years) nearby neutron stars have been well observed and studied. The simple presence of such verified old stars leads to the conclusion that no black hole has consumed the star and as we shall argue, this can lead to very severe constraints on the properties of certain types of asymmetric dark matter. We should also mention that additional constraints on asymmetric dark matter can be imposed on different ways (e.g. from asteroseismology [15][16][17], from effects on the transport properties of the neutron stars [18] and/or hybrid dark matter rich compact stars [19,20]).
One can easily figure out that fermionic WIMPs due to the fact that they have to overcome Fermi pressure, require a huge number in order to collapse i.e. * Electronic address: kouvaris@cp3.sdu.dk 3 where M pl and m are the Planck mass and WIMP mass respectively. This number of WIMPs is very difficult to be accumulated within a few billion years and with dark matter densities similar to the ones of the earth. However, this required number for gravitational collapse is reduced significantly in the case of attractive Yukawa forces among the WIMPs [8].

II. ASYMMETRIC BOSONIC DARK MATTER
In the case of asymmetric bosonic WIMPs, the necessary WIMP number for collapse is much smaller because there is no Fermi pressure and only the uncertainty principle keeps particles from collapsing. The collapse takes place once the momentum becomes smaller than the selfgravitational potential energy.
where M = N m is the total mass of the WIMP cloud. A more accurate and generic estimate that includes the effect of self-interactions gives [21] M crit = 2 π Although self-interactions between WIMPs can be quite general in nature, without loss of generality, we can assume that the self-interaction can be approximated well by a λφ 4 interaction term. At the no interaction limit λ = 0 we trivially get the critical mass mentioned above (up to factors of order one). The accretion of WIMPs for a typical 1.4M 10 km neutron star taking into account relativistic effects has been calculated in [3]. The total mass of WIMPs accreted is where ρ dm is the local dark matter density, and the "efficiency" factor f = 1 if the WIMP-nucleon cross section satisfies σ > 10 −45 cm 2 , and f = σ/(10 −45 cm 2 ) if σ < 10 −45 cm 2 .
One can easily check that M acc can be larger than M crit practically for all masses larger than ∼ 100 keV. To form a black hole, satisfying the condition (2) is necessary but it is not sufficient. One should make sure that after the WIMPs have been captured, they slow down and thermalize with nuclear matter concentrating within a small thermal radius. Failing to satisfy this condition, even if the condition (2) is satisfied, does not necessarily lead to the formation of a black hole, since WIMPs would not be confined in a tiny region. The thermalization time scale has been estimated in [5] and [3] As one can observe, despite the Pauli blocked interactions between WIMPs and nucleons, unless they are very heavy, WIMPs thermalize in less than a year. Having thermalized with nuclear matter, WIMPs concentrate in the center of the star within a thermal radius that can be easily obtained by use of the virial theorem where k is the Boltzmann constant, T c is the temperature at the core of the star, and ρ c = 5 × 10 38 GeV/cm 3 is a typical value for the neutron star core density.
Once the WIMPs are thermalized and if sufficient number is accumulated in the star, there are two different events that take place, the time order of which depends on the WIMP mass. One is the self-gravitation of the WIMP sphere and the second is the formation of a Bose Einstein condensate (BEC). Self-gravitation takes place once the mass of the WIMP sphere inside the thermal radius becomes larger than the mass of the neutron star within the same radius. In other words, this happens once WIMPs start feeling strongly their own gravitational field. For this to happen the WIMP sphere should have a mass that satisfies On the other hand, BEC formation takes place once the WIMP number density is One can easily check that for WIMPs roughly lighter than 10 TeV, the accumulated WIMPs within r th meet first the condition for BEC formation. We are going to study these two cases (m < 10 TeV and m > 10 TeV) separately since events unfold with different order. For WIMPs lighter than 10 TeV, one can estimate that the total number of WIMPs needed to form a BEC is N BEC 2 × 10 36 . Any accumulated WIMPs on top of this number goes directly to the ground state of the BEC state. The radius of the BEC state is (8) As it can be seen, r BEC << r th and therefore WIMPs in the ground state can become self-gravitating much faster than what Eq. (6) predicts. In fact we can appreciate this if we substitute r th by r BEC in Eq. (6). This leads to the condition If Eqs. (2), (7), and (9) are satisfied, a black hole is going to be formed. Once the black hole is formed, its fate is determined by its initial mass M crit . One the one hand, the black hole is accreting dark matter and nuclear matter from the core of the star. This tends to increase the black hole mass. On the other hand, emission of photons and particles in general via Hawking radiation tends to reduce the mass of the black hole. The black hole mass evolution is determined by where c s is the sound speed at the core of the star, and f is a dimensionless number that in general depends on the number of particle species emitted and the rate of rotation of the black hole. We have used a spherically symmetric Bondi accretion of matter into the black hole. By inspection of Eq. (10) it is apparent that there is a critical value of the black hole mass M above which accretion always wins, while below, Hawking radiation reduces the mass of the black hole which in turn it increases even further the rate of Hawking radiation leading eventually to the evaporation of the black hole. This critical mass has been estimated if one considers only photons in [6] M > 5.7 × 10 36 GeV.
The mass becomes slightly larger [10] if one includes also other species that can be emitted (e.g gravitons, neutrinos, quarks, leptons etc). Comparison of Eq. (2) (with λ = 0) to Eq. (11) shows that WIMP masses larger than 16 GeV lead to black hole masses below the limit of Eq. (11). This means that for masses larger than 16 GeV, black holes evaporate and their effect is to heat up the star as they evaporate. However this does not lead to a dramatic effect like the destruction of the star. This 16 GeV mass limit becomes slightly smaller if more Hawking radiation modes are included. Finally there is one last constraint that should be satisfied. WIMP masses cannot be arbitrarily small because for small WIMP masses, after WIMPs have thermalized, those in the tail of the Maxwell-Boltzmann distribution have large enough velocities to escape from the star. This evaporation effect can be ignored for WIMP masses m > 2 keV [6]. If the accreted dark matter mass within a billion years M acc is larger than M crit of Eqs. (2), and (7), (9), and (11) are satisfied, the WIMPs form a black hole that can destroy the star. There are some subtle issues regarding how fast the black hole consumes the star that have been addressed to some extend in [6]. The constraints on asymmetric bosonic dark matter are depicted in Fig. 1. As it can be seen, depending on the WIMPnucleon cross section, WIMP candidates from 100 keV up to roughly 16 GeV are severely constrained by the existence of nearby old neutron stars. The constrained region is bound at 100 keV due to the fact that below that mass accretion is not sufficient to acquire M crit from Eq. (2). These constraints can be enlarged down to 2 keV (the limit from WIMP evaporation we mentioned before) as long as we consider old neutron stars in globular clusters with ρ dm 30 GeV/cm 3 . Now we can consider the case where the WIMP mass is larger than 10 TeV and therefore self-gravitation of the WIMP sphere happens before BEC formation. As we mentioned above, black holes of critical mass (2) with WIMP masses roughly larger than ∼ 16 GeV, do not survive due to Hawking radiation. Therefore one should expect that black holes of M crit (of Eq. (2)) formed out of 10 TeV WIMPs (or heavier) would evaporate quite fast. However, since self-gravitation takes place before BEC, and the self-gravitating mass of Eq. (6) for m > 10 TeV is much larger than the crucial mass for the survival of the black hole of Eq. (11), there were speculations in the literature [7,9,10] that constraints can be imposed also for m > 10 TeV. The claim was that instead of forming a black hole of M crit that is below the surviving threshold for Hawking radiation, a much larger black hole coming from the collapse of the self-gravitating WIMP sphere M sg forms, that due to its larger mass can grow and destroy the star, thus imposing constraints on this part of the parameter space of asymmetric bosonic dark matter. However we review here the argument that was put forward in [23] that demonstrates that the formation of smaller (non-surviving) black holes of mass M crit is unavoidable and therefore the M sg instead of collapsing to a single large black hole, it forms a series of black holes of M crit that evaporate one after the other, thus resulting to no constraint for WIMP masses with m > 10 TeV.
In order for the WIMP sphere to collapse, the whole mass should be confined within the Schwarzschild radius r s = 2GM of the black hole. The density of WIMPs just before forming the black hole would be n BH ∼ 3(32πG 3 M 2 sg m) −1 ∼ 10 74 cm −3 (GeV/m)(M sg /10 40 GeV) −2 . It is easy to see that this density is higher from the density required for BEC formation of Eq. (7). This means that unless the WIMP sphere collapses violently and rapidly, it should pass from a density where BEC is formed. As the selfgravitating WIMP sphere of mass M sg contracts, at some point it will reach the density where BEC is formed. Any further contraction of the WIMP sphere will not lead to an increase in the density of the sphere. The density remains that of BEC. The formation of BEC happens on time scales of order [22] t BEC ∼ /k B T ∼ 10 −16 s, i.e. practically instantaneously. Further shrinking of the WIMP sphere results in increasing the mass of the condensate rather than the density of non-condensed WIMPs. This process happens at a time scale which is determined by the cooling time of the WIMP sphere as discussed below. As we shall show, this cooling time is the relevant time scale for the BEC formation. As in the previous case, the ground state will start being populated with WIMPs which at some point will become self-gravitating themselves. This of course will happen not when Eq. (9) is satisfied. Eq. (9) was derived as the WIMP ground state becomes denser than the surrounding nuclear matter (since the dark matter that is not in the ground state of the BEC is less dense). Here, the condition is that the density of the ground state of the BEC should be larger than the density of the surrounding dark matter (that is already denser than the nuclear matter at this point). The condition reads Once the BEC ground state obtains this mass, the ground state starts collapsing within the collapsing WIMP sphere. Any contraction of the WIMP sphere does not change the density of the sphere but only the density of the ground state. M BEC, sg is smaller than M crit and therefore the BEC ground state cannot form a black hole yet. However as the ground state gets populated at some point it reaches the point where its mass is M crit and this leads to the formation of a black hole of mass M crit and not M sg . The evaporation time for such a black hole of M crit (composed of WIMPs of mass m) is The only way that such a black hole can be maintained in life is by adding fast new matter inside in a rate that is higher than the Hawking radiation. As we showed earlier accretion of nuclear matter from the core of the star is not sufficient for m larger than roughly 16 GeV. Accretion of dark matter from the rest of the WIMP sphere is also not sufficient. The accretion of non-interacting collisionless particles in the non-relativistic limit onto a black hole is given by [24] where ρ dm is mn BEC , and v ∞ is the average WIMP velocity far away from the black hole. Using the virial theorem, we take v ∞ = GM/r, M being the total mass of the WIMP sphere and r is the radius of the WIMP sphere. For the first black hole M = M sg and r = r 0 where r 0 is the radius of the WIMP sphere when it reaches the BEC density i.e.
For WIMP sphere masses ranging from Eq. (2) to Eq. (6) the Hawking radiation dominates overwhelmingly over the WIMP accretion despite the large WIMP density, and therefore the black hole evaporation time is given accurately by Eq. (13).
There are another two ways where the black hole can be maintained in life i.e. if the WIMP capture time scale and/or the cooling time scale are smaller than τ . The capture time scale is defined as the time it takes for the star to capture a number of WIMPs with a total mass of M crit . If this time is smaller than τ , a new black hole of M crit will form before the extinction of the previous one. Therefore it is possible that gradually the mass of the black hole could increase after successive black holes form before previous ones have evaporated. However, using Eqs. (2) and (3) we get the capture time scale which is larger than τ for m > 10 TeV even for neutron stars located in high dark matter density globular cluster cores. Capture rates of new WIMPs cannot replenish the black hole mass rate loss to Hawking radiation.
Finally one should estimate the cooling time scale for the WIMP sphere. Once self-gravitating, the WIMP sphere does not collapse to a black hole just because WIMPs still have a lot of energy. This energy can be lost via collisions with nucleons. This way the WIMPs lower their energy and the WIMP sphere contracts. Let us assume that a black hole of mass M crit has been formed and it is going to evaporate in τ . Once formed the black hole is doomed to evaporate unless the WIMP sphere contracts fast enough, so it replenishes the BEC ground state with another M crit mass that collapses and forms a second black hole, before the first one has evaporated. In that case, the mass of the black hole would be reinforced with new mass making possible for the black hole to grow. As we shall show the time scale it takes for the WIMP sphere to populate again the BEC ground state is larger than the evaporation time and therefore the black hole cannot survive. The energy that the WIMP sphere has to lose in order to admit M crit in the BEC ground state is where M is the total WIMP mass and r is the radius of the WIMP-sphere. For the first black hole M = M sg and r = r 0 . This energy must be lost by WIMP-nucleon collisions. The time it takes on average for a WIMPnucleon collision (taking into account the Pauli blocking effect for the degenerate nucleons) is where p F 1 GeV is the nucleon Fermi momentum at the core of the star, n the number density of nucleons, m N the nucleon mass, σ the WIMP-nucleon cross section, and is the kinetic energy of the WIMP. can be estimated by use of the virial theorem One should keep in mind that each WIMP-nucleon collision leads to a loss of energy δ = 2(m N /m) .Therefore in order for N = M/m WIMPs to lose the excess δE energy colliding in τ col time scales with nucleons losing each time on average energy δ , it takes This time is shorter by a factor of a few than the black hole evaporation time of Eq. (13). Note, however, the strong dependence of both quantities on the WIMP mass m. Already for masses m > 13 TeV the black hole evaporation time becomes shorter. The cooling time derived above refers to the formation of the first black hole. Subsequent black holes require progressively longer times. This is easily seen from Eq. (20) because t cool scales as 1/M 5/3 (recall that r 0 scales as M 1/3 from Eq. (15)). Thus, the more black holes have formed and evaporated, the smaller is the remaining WIMP mass M and the longer the time needed to form the next one.

III. CONCLUSIONS
We reviewed the current status of constraints on nonself-interacting asymmetric bosonic dark matter based on observations of old neutron stars. We showed that severe constraints can be imposed for these candidates from the keV to the GeV range, excluding all candidates with WIMP-nucleon cross sections like the ones favorited by the DAMA experiment within this mass range. We also reviewed the arguments why these constraints cannot be extended to candidates with masses in the TeV range or higher.