Access the full text.
Sign up today, get DeepDyve free for 14 days.
Loading next page...
/lp/arxiv-cornell-university/constraining-light-fermionic-dark-matter-with-binary-pulsars-3vGIpZ8nS7
References
References for this paper are not available at this time. We will be adding them shortly, thank you for your patience.
- ISSN
- 2212-6864
- eISSN
- ARCH-3330
- DOI
- 10.1016/j.dark.2019.100343
- Publisher site
- See Article on Publisher Site
Abstract
A binary system embedded in a Dark Matter (DM) background may experience a change in its orbital period due to dynamical friction as the binary moves through a wind of DM particles. Since the orbital motion variation depends on the DM environment, such a phenomenon provides an intriguing way of constraining the properties of DM from timing pulsar observations. We compute such a perturbative effect on the binary evolution considering that DM is constituted of degenerate gas of free fermions. The analysis point out that the secular change of the orbital period is more sensitive, and likely measurable, to degenerate fermions with masses & 50 eV, depending slightly, but still being distinguishable, on the binary star configuration (e.g. NS-NS, NS-WD and WD-WD). Interestingly, we find that NS-NS binary systems with large orbital periods, P & 100 days, experience larger orbital period decays. We also show that this effect is clearly increased, under the former conditions, in binaries orbiting small DM halos, which correspond to extragalactic pulsars. This situation represents the best astrophysical scenario to test such effects of light fermionic DM. We use some available measurements of the orbital period time-derivative for long-period binaries in the Milky-Way to quantify more realistically this effect. For instance, measurements of the J1713+0747 pulsar set an upper bound on the fermion mass of m . 1 keV. This bound can be considerably improved by using pulsar timing observations of extragalactic pulsars. Under this perspective, high precision of timing pulsar observations will reveal whether DM dynamical friction effect may be tested with the upcoming generation of surveys leading to the possibility of constraining more strongly the properties of light fermionic DM. Keywords: Dark matter: degenerate gas of free fermions. Binary pulsars: Orbital period decay 1. Introduction [17], such as the core-cusp [18, 19] and the too big to fail [20, 21] problems. Despite great efforts in both direct [1–3], and indirect Non-interacting Fermionic DM models have been advo- [4, 5] detection searches and colliders [6, 7], the nature 1 cated for the formation of cores at the center of galaxies of DM still remains unrevealed, thus demanding more so- due to the degeneracy pressure which prevents gravita- phisticated and sensitive experiments. A less conservative tional collapse. This unequivocal consequence of the Pauli position on the DM problem is to assume that the ob- exclusion principle has been explored, with great astro- served non-baryonic gravitational effects might rather be physical interest, in order to describe adequately the ob- a manifestation of some sort of modified theory. See e.g served cores in dwarf galaxies of the Milky-Way [22–28]. Ref. [8] among the most attractive alternatives that may As discussed, fermionic DM may play an important account for observations. 2 role in the description of cored DM halo . We focus then Nevertheless, the observational evidences coming from on self-gravitating degenerate gas of free fermions, also cosmological and astrophysical scales, such as Ref. [9–11] known as light (sub-keV) fermions, as a simple but well and Ref. [12, 13] respectively, support convincingly the motivated realization of non-interacting fermionic DM to idea of a non-baryonic component in the Universe. Ac- describe the DM distribution in galaxies. For this end, we cordingly, there have been proposed some DM models with promising theoretical predictions; being among the most popular candidates, apart from the weakly interactive mas- Such a core corresponds to a polytropic core with index n = 3/2. See for instance Refs. [22, 23]. sive particles (WIMPs)[14]: ultra-light bosons (see e.g. It is important to state that there is still a debate in the literature [15]) and fermionic DM (e.g. sterile neutrinos [16]). These about the inner structure of galaxies and why some low-mass spiral latter candidates have been postulated to explain, with galaxies are well described by cored profiles while others not [29]. Trying to reconcile such a discrepancy is a challenge for DM-only great success, the main problems of the CDM paradigm models which require inevitability the help of baryonic physics or self-interactions. For a complete discussion of the small scale crisis and unified solutions see Ref. [30, 31]. Email address: gabrielphysics@gmail.com (L. Gabriel G´omez) Preprint submitted to Physics of the Dark Universe July 12, 2019 arXiv:1906.10316v2 [astro-ph.GA] 10 Jul 2019 mention briefly the main advances in the development of been found the conditions under which such a perturba- this topic. A comprehensive attempt to describe the kine- tive effect in compact-star binaries, for different DM mod- matical data of dwarf galaxies of the Milky-Way under this els, may become larger than the one due to gravitational approach was performed in Ref. [27]. In there, analysis of wave emission [41]. Furthermore, the DM dynamical fric- the best-fit of the observed velocity dispersion set a nar- tion effects of binary pulsars has been used as probes of a row range for the fermion mass of 100−200 eV, describing Galactic DM disc scenario. [42]. Importantly, these results thereby the structural properties of the classical satellite exhibit different theoretical predictions of orbital period galaxies. A step forward was done by Ref. [28] where a decays that make DM dynamical friction a promising tool more robust analysis of the kinematical data was carried from future observations of timing pulsars to constrain or out , along with the inclusion of a thermal envelope, yield- rule out DM models. ing a fermion mass lying in the range of 70− 400 eV. We As mentioned before, degenerate gas of free fermions conclude this review by mentioning that analysis of rota- have recently garnered attention into the astrophysical tion curves of the Milky-Way has also provided constraints context because they yield the possibility of describing suc- on the allowed range of the fermion mass of 75 − 104 eV cessfully the structural properties of dwarf galaxies of the [33]. Milky-Way [27]. In view of this, the astrophysical implica- Despite the compelling indirect evidences for DM through tions of degenerate fermionic DM in galactic halos deserve its gravitational effect, new astrophysical scenarios have further examinations. We study a novel astrophysical ap- brought the attention as a feasible alternative to study plication of such a DM model by assessing the dynamical the properties of DM. Binary pulsars have been used, for friction effect it produces on the orbital evolution of binary instance, to infer some information about the properties pulsars. We analyze the effect of varying orbital param- of the gravitational potential and the interstellar matter eters and the fermion mass on the secular change of the content of the Galaxy [34]. Interestingly, these objects orbital period such that we can devise the best astrophys- may also have a measurable change in the orbital period, ical scenario to constrain the properties of light fermionic as the two stars spiral inward toward each other, which DM. The secular change in the orbital period due to DM can be attributed to energy loss due to the emission of dynamical friction of some binary stars with measured or- gravity waves [35]. These are some of the few objects bital period time-derivatives is also evaluated. Thus, we with extraordinary accuracy in measurements which make show that it is possible to constrain the fermion mass by them excellent tools for studying gravity (in the strong using available timing pulsar observations. For the sake of field regime) and matter around them [34]. Furthermore, example, we use the J1713+0747 binary to set the upper such compact-stars binaries offer an invaluable opportu- limit m . 1 keV. Finally, we comment some prospects nity to test deviations of General Relativity though with of timing pulsar observations of extragalactic pulsars to no reported observational signals yet [36, 37]. improve such a bound in the near future that makes our Apart from this, binary pulsars may also offer inter- results testable predictions. esting possibilities to know about the medium they are We outline now the general structure of this paper. embedded like the intergalactic medium. For instance, dy- Secs. (2)-(4) are devoted to review the underlying approach namical friction, a drag force that is presented in the evo- of the presented astrophysical application. We start by de- lution of many astrophysical systems , may perturb the scribing in Sec. (2) the generalities of the fully degenerate orbital motion of binary pulsar as they move through a fermions in DM halos. We proceed in Sec. (3) with the def- wind of collisionless field of DM. This drag induces a wake inition of the Chandrasekhar’s dynamical friction formula. of medium particles on the object. As a result of this With this at hand, we present in Sec. (4) a general illus- perturbative effect, due to the interaction between each tration of how dynamical friction may originate a secular binary component with its respective wake, there may be change of the orbital period in binary stars. Later in the a secular change of the orbital period associated to this same section we consider such a degenerate gas of fermions phenomenon [39, 40]. to account for the DM background. Finally, we compute This interesting astrophysical scenario has been gain- and show in Sec. (5) changes in the orbital period in binary ing attention as an alternative way to constrain the prop- pulsars under the aforementioned setup. The main con- erties of DM from a purely gravitational interaction. The clusions of this study and the observational perspectives dynamical friction effect in binary pulsars has been pos- of this work are discussed at the end of this manuscript. tulated and studied in binary systems with a DM back- ground modeled by the NFW profile [40]. Recently, it has 2. General description of degenerate Fermionic DM See also Ref. [32] for further considerations in fitting the kine- We are interested in the strong degeneracy limit of non- matical data of dwarf spheroidal galaxies with degenerate fermions. interacting fermions, i.e, at zero temperature (T → 0), For a more complete discussion of how such ranges have been ob- however, we start by describing a system of self-gravitating tained, we refer the interested reader to the introduction of Ref. [28]. 5 fermions, keeping finite temperature T and chemical po- See e.g. [38] for the role of dynamical friction in the inspiral of dwarf galaxies within DM halos and the orbital evolution of black tential μ, with the Fermi-Dirac distribution function in the hole binaries in a stellar medium. 2 non-relativistic regime Making some analytical treatments, they present in par- ticular the mass-radius relation as f(p) = , (1) 1+β(E−μ) 5K ′ 5 −3 M = −4π θ (ξ )ξ R , (6) 8πG where β is the inverse temperature defined as β = 1/K T and E = p /2m is the single particle kinetic energy. K f B where ξ = 3.65 and θ (ξ ) = −0.203 are numerical con- 1 1 is the Boltzmann constant and m is the fermion mass. stants. It can be easily read off that the total mass scales Now, if we assume that DM is composed of fully degen- −8 −3 with the fermion mass and with the radius as M ∝ m R . erate free fermions, then they will occupy all the quantum From this latter relation is set that large self-gravitating states with momentum lower that the Fermi pressure p . systems can be reproduced only considering small values The Fermi pressure corresponds to the the highest level of m and vice versa. This model is then characterized by occupied at T = 0. In the non-relativistic regime, it is the two free parameters m and ρ we shall have in mind f 0 possible to write p = mv. Likewise, one may write the in all our computations. occupation number in terms of the momentum as n = R On the other hand, if one assumes a gravitationally 4π 2 g p dp, where g denotes the particle spin degen- h 0 bound DM object composed of a degenerate fermion gas, eracy and h is the Planck constant. Taking g = 2 for there exists a densest packing of the DM phase space dis- 3 3 our case, then n = 8πp /3h and, along with the mass tribution. Hence, the phase-space density of the DM must density ρ = mn, we get the the Fermi momentum p = 8 not exceed that of the degenerate Fermi gas . In our case, 3 1/3 (3h ρ/8πm ) . With this in mind, the pressure can be it will correspond exactly to the maximum phase space. If computed from the definition one demands that the Fermi velocity v does not exceed the escape velocity v of the self-gravitating degenerate 2/3 p 4 2 8π p h 3 5/3 gas of fermions of mass M and volume V = 4/3πR , p = dp = ρ . (2) 3 8/3 3h 2 2 8π 0 5m namely v ≤ v , it leads to p + m f ∞ 1/3 r Hence, the equation of state of a degenerate Fermi gas in 9Mh 2GM ≤ , (7) the non-relativistic regime follows the polytropic equation 2 3 4 32π R m R 2/3 γ 3 h P = Kρ , with K = and polytropic index 8/3 8π 5m which sets an upper bound for the Fermi velocity that any γ = 5/3. single particle should satisfy inside the DM bound object The density profile has to satisfy both the hydrostatic [44]. We shall keep in mind this relation as a condition equation that any single particle, with a given velocity distribution, dP GM(r) = − ρ(r), (3) must satisfy inside galactic halos. dr r and the continuity equation 3. DM dynamical friction in binary stars M(r) = 4πr ρ(r)dr, (4) We describe briefly the general scenario where an ob- ject experiences a dynamical friction as it moves through to describe a static self-gravitating system in thermody- a collisionless medium of field particles. There are numer- namic equilibrium. This set of equations can be written in ous applications of this (purely gravitational) effect in the the form of the Lane-Emden equation with index n = 3/2 astrophysical context since dynamical friction is a natu- ral force which manifests, for instance, in a change in the 1 d dθ 2 3/2 orbital evolution of most of the known objects (see foot- ξ = −θ(ξ) . (5) ξ ξ dξ note 4). This drag force induces a wake, as a result of the gravitational interaction of medium particles on the ob- We have followed closely the scaled numerical solution of ject, with a characteristic overdensity proportional to its this equation given by Ref. [27]. They introduced the di- 6 7 3/2 mass [38]. mensionless parameter θ and α so that ρ(ξ) = ρ θ(ξ) The drag force experienced by a test body of mass and r = ξα, where ρ is the central density and α = −1/3 m ≫ m, where m is the DM particle mass, and with 1/2 (5Kρ /8πG) with density vanishing at R = ξ α. orbital velocity v moving through the DM background characterized by the velocity distribution function f(u), is Only in this section we have used β and α to denote the inverse temperature and the dimensionless parameter respectively. We use however the same letters in the remainder of this paper as angles. This argument has been considered, since the pioneering work of Gunn & Tremaine [43], to put constraints on the fermion mass for a In our computations will be helpful to use, as a well approxima- tion, the function ρ(ξ) = ρ cos [ ξ] found in Ref. [32]. given DM dominated object. 3 given by the Chandrasekhar’s dynamical friction formula the dynamical friction force, making them more attractive [38, 45] targets to test this effect from observations. Thus, we can apply this approach under the above conditions to binary v˜ max systems such as NS-NS/NS-WD [46] and WD-WD [47]. 2 2 3 2 2 f = −4πG m m d uf(u) ln (v˜ − u ) fr,i i i Gm Some important assumptions and considerations are sum- marized to ensure the viability of this approach as follows: v∞ u + v˜ v˜ v˜ i i i + d uf(u) ln − 2 , (8) u − v˜ u v˜ • the medium is treated as a collisionless gas as long as v˜ i the mean free path of the DM background is much where the integral in the first term accounts for the frac- larger than the size the objects [39] tion of particles moving slower than the object, while the integral in the second term, referring to fast particles, is • the Chandrasekhar’s DF formula neglects the inter- limited by the escape velocity v . The above equation ∞ action of one component with its companions wake. for the dynamical friction force considers the orbital ve- This condition is guarantee provided that the orbital locity of each object with respect to the DM wind relative velocity is smaller that the velocity dispersion of the to the center of mass of the binary system v ˜ = v + v , i i w DM background such that the wake of one object dis- with v = v (cos α sin β, sin α sin β, cos β) and β and α perses before the arrival of its companion [40]. This w w being the angles between the wind velocity vector and the also implies that we must consider only binaries with perpendicular axis of the binary orbital plane and the pro- large orbital periods P > O(1) days jection of the wind velocity vector with an axes lying in • linearity of the Chandrasekhar’s DF formula requires the orbital plane, respectively. We follow the discussion L/a ≪ 1, where L is the size of the components wake presented in Ref. [41] in which they present at least two and a is the orbital separation. Here L is determined different contributions of wind velocities: bound and un- by Eqn. (11) (which is identified as R ) and is de- bound binaries to the galaxy potential depending on the fined as the radius of the sphere of gravitational in- magnitude of the kick velocity in binaries with NS com- fluence of the test body ponents. Thus, for the sake of generality of our results, we assume v = v + v so that we take values for it w rot T • the latter condition leads that binary stars can be −1 −1 ranging from 10 km s all the way to 1000 km s . Here considered as point masses v is the transversal velocity of the system. A crucial ingredient in the Chandrasekhar’s dynamical 4. Orbital decay of binary pulsars by dynamical friction formula (Eqn. (8)) is the Coulomb logarithm which friction accounts for the finite size of the system and is usually defined as the ratio of the maximum and minimum impact We review the general approach underlying the DM parameters for encounters, respectively b and b , i.e. max min dynamical friction effect on the orbital period evolution of binary systems. We follow the derivation of the perturbed max log Λ ≡ log . (9) Keplerian orbits, from the osculating formalism [48] as was min presented in Ref. [40], to compute the orbital period de- b is assumed typically of the order of the size of the max cay due to the aforementioned effect. We start by writing system (for our case the orbital separation a), and b is min the equation of motion for the two body system consider- defined as the impact parameter for a 90 deflection [38] ing the effects of the dynamical friction force of the DM background b ≈ a, b = max(r , R ), (10) max min h A Gm m 1 2 where r is the half-mass radius of the subject system and h m ¨r = ± r + f , (12) i i fr,i R is the radius at which a particle of the surrounding medium is affected by the sphere of gravitational influence where r = r − r . Accordingly the relative acceleration 2 1 of the test body, namely: between two bodies is written as Gm i GM R = . (11) A,i v˙ = − r + f, (13) v˜ For comparison purposes, we shall take values of b and where v = r˙ and f = a ηv + a v is the perturbing max 1 2 w b so that log Λ ≈ 20. dynamical friction force measured on the center of mass. min As discussed in Ref. [40], we deal with binary systems Note that in the absence of this force, the orbital velocity with orbital velocities smaller than the velocity dispersion obeys certainly a Keplerian motion v = Ω r , with Ω and 0 0 0 of the DM medium so that each binary component does r being the angular velocity and orbital separation, re- not interact with its respective companion’s wake. We are spectively. We have also introduced compact definitions: limited then to consider only binary systems with large or- η = μ/M, μ = m m /M and M = m + m , to describe p c p c bital periods which are, as we shall see, more sensitive to the binary system. Here, the subscripts p and c refer to 4 the primary component and its companion, respectively. Note that in the above definitions we have considered for With this in mind, the perturbed orbital elements, needed simplicity that the Coulomb logarithm in the dynamical to compute the orbital motion due to the drag force, can friction force is equal for the pulsar and its companion −1 be then computed. In particular, the expression for the which is a good approximation for v > 100 km s . change of the orbital separation reads We realize also that different values of the initial phase α r does not introduce significant changes in our results, then setting α = 0 is a suitable choice. a˙ = 2 S(t), (14) Based on the above formulation, we compute in the GM next section the secular change of P due to DM dynamical where a is the semiaxis major. It is important to high- friction solely. Thus, for the sake of generality of our con- light that the source term S(t), in the perturbed orbital clusions, we shall not include additional effects that may elements have been defined in terms of the dynamical fric- contribute to the secular change of the orbital period such tion force as well as the wind velocity vector according to as mass loss of the star components or accretion of DM par- Ref. [40] ticles onto the binary components. The justification relies on the fact that we are interested in long binaries periods, S(t) = a ηv − a v sin β sin(Ω t − α), (15) 1 2 w 0 low-mass and compact star binaries in whose case the for- In short, the drag force produces a change of the orbital mation of accretion disk is unlikely either by Roche-lobe separation with time leading to an intrinsic change of the overflow or stellar winds. Hence accretion of matter from orbital period P = 2π/Ω given by [49] one component into the other could occur only via Roche b 0 lobe overflow for extremely short binary periods near the P 3 a˙ merging process. We are also ignoring the effects of DM = . (16) P 2 r b 0 accretion because of the unknown cross section between DM and baryonic matter inside the stars. For interest- This simple relation along with the corresponding one for ing discussion on this latter issue in binary systems see the change of the orbital separation Eqn. (14) yields the Ref. [50]. Furthermore, we are not interested in binaries expression for the orbital period time-derivative located in globular clusters where stellar interactions with ˙ encounters are more likely neither in regions inside galax- P (t) = 3P [a η − a Γ sin β sin (Ω t − α)]. (17) b b 1 2 0 ies where gas and dust may perturb the binary system and Importantly, the measured secular change in the orbital induce an orbital decay [51]. period is obtained after averaging over one period P , namely 5. Numerical results: secular change in the orbital ˙ ˙ hP i = P (t)dt. (18) period b b Let us now compute the orbital period time-derivative What concerns to the overall formulation, we have intro- according to the underlying description presented in Secs. (2)- duced the same definitions as in Ref. [40] for both an eas- (4). To do so, firstly, we fix two of the three free parame- ier comparison of the results and simple implementation ters: v , P and m (as it shall be specified in each plot) w b f of the DM model in consideration Γ = v /v, Δ = Δ± 1, w ± and then perform the numerical computation varying the Δ = 1 − 4η. The coefficients a can be written in terms −1 9 remaining one in the pertinent range. This is: 10 km s of the integral velocity contribution function −1 . v . 1000 km s , 0.1 days . P . 100 days and 50 eV w b . m . 400 eV. The latter range is thus set because it b = , (19) v˜ describes properly galactic halos [27] and the structural properties of dwarf galaxies [27, 28, 32] being then of as- as trophysical interest for the purpose of this work. Thus, for m = 90, 200 eV, their associated central densities shall a = −A(b + b ), a = (b Δ + b Δ ), (20) 1 1 2 2 1 + 2 − −3 be ρ = 0.355, 5.6 GeV cm . The former value is within the suitable range to fit the Milky-Way Galactic observa- with A = 4πG Mρ(r) and ρ(r) is the mass density of a tions [33] while the latter one corresponds roughly to the degenerate gas of free fermions found in Sec.2. The depen- one inferred to describe a typical dwarf galaxy [27]. We dence of the dynamical friction force on the fermion mass consider also for the sake of example the following binary is intrinsically inhered through its relation with the den- system configurations: NS-WD with m = 1.3 M and p ⊙ sity profile. This incorporation is justified from the fact m = 0.2 M , NS-NS with m = m = 1.3 M and WD- c ⊙ p c ⊙ that we are treating with a gas of self-gravitating fermions WD with m = 0.5 M and m = 0.25 M . p ⊙ c ⊙ which take into account microphysical properties of a sys- tem such as the particle mass to describe the mean density. 5.1. The effect of the wind velocity Having established the range of orbital parameters and The quantity I is defined by the term in parenthesis in Eqn. (8). fermion mass, we start our analysis by investigating the 5 -15 changes occur as was investigated in Ref. [41] for different DM models. Nevertheless, for the sake of generality of our -16 conclusions, we shall take henceforth β = π/2 in all our numerical estimations. -17 Taking advantage on this preliminary result, let us do a -18 short commentary beforehand, on the effect of the fermion mass on P . We can infer, from Fig. 1, that it acts only as a -19 scaling factor over all the solutions. Hence, the displayed NS-WD solutions in both panels of Fig. 1 preserve the same be- NS-NS m = 90 eV -20 f 10 havior for any binary configuration as can be simply seen. WD-WD However, for the solution m = 200 eV (bottom panel), 10 50 100 500 1000 which reproduces a typical dwarf galaxy according to the -1 mass-radius relation Eqn. (6), the DM dynamical friction v [Km s ] effect is almost one order of magnitude larger that the one -14 10 found for the Milky-Way, i.e. for the solution m = 90 eV. This result yields us some insights about the most promis- -15 ing scenario to search for such an effect in host galaxies. -16 -17 -16 NS-WD NS-NS NS-WD -18 WD-WD NS-NS -17 m = 200 eV 10 WD-WD -19 m = 90 eV -18 f 10 50 100 500 1000 -1 v [Km s ] -19 Figure 1: Secular change of the orbital period as a function of the -20 DM velocity wind. For this computation we have adopted the values 10 P = 100 days and β = π/2. The curves refer to numerical solutions for different binary systems as indicated below: 1.3–0.2 M NS- 0.1 1 10 100 1000 WD (black solid line), 1.3–1.3 M NS-NS (blue dashed line) and 0.25–0.50 M WD-WD (red dotted line). Top Panel: This solution ⊙ P [days] −3 corresponds to m = 90 eV and ρ = 0.355 GeV cm . Bottom 0 -14 panel: This solution corresponds to m = 200 eV and ρ = 5.6 GeV f 0 NS-WD −3 cm . NS-NS -15 WD-WD effects of the wind velocity on the secular change of P b -16 m = 200 eV for different binary configurations. The results are shown in Fig. 1. The two plots correspond to a fermion mass of -17 m = 90 eV with central density ρ = 0.355 GeV cm f 0 (top panel) and m = 200 eV with ρ = 5.6 GeV cm f 0 -18 (bottom panel), respectively. We can see that P does not exhibit large changes in magnitude (in logarithmic scale) -19 −1 from v ≈ 40 km s for the NS-NS binary and from 0.1 1 10 100 1000 −1 v ≈ 100 km s for both NS-WD and WD-WD binaries. P [days] It is also appreciable multiple dips for a given binary star configuration. The physical meaning of those dips in all Figure 2: Secular change in the orbital period as a function of the −1 the plots is due to the change of sign of the orbital period orbital period. We have here adopted the values v = 100 km s and β = π/2. The curves refer to numerical solutions for different time derivative leading to an orbital shrinking to widening binary systems as indicated below: 1.3–0.2 M NS-WD (black solid or vice-versa, as appropriate. Such changes of sign occur line), 1.3–1.3 M NS-NS (blue dashed line) and 0.25–0.50 M WD- ⊙ ⊙ for certain values of the binary parameters and the angle β. WD (red dotted line). For these solutions we have initially fixed the −3 As a general trend, we realize that the trivial choice β = 0 values of the fermion mass: m = 90 eV with ρ = 0.355 GeV cm −3 (top panel) and m = 200 eV with ρ = 5.6 GeV cm (bottom does not provide such changes of sign. On the contrary, f 0 panel). β = π/2 leads to multiple changes of sign as was also identified in Refs. [40, 41]. We are not interested however in exploring the ranges of the parameters under which such P P b b P P b b 5.2. The effect of the orbital period -13 NS-WD Before proceeding with our analysis, we remind that NS-NS -14 the two chosen values of fermion masses, adopted here for WD-WD a self-gravitating degenerate gas of free fermions, lead to -15 a suitable description of both large (m ∼ 90 eV) and small (m ∼ 200 eV) galaxies. Moreover, we should take -16 in mind these values only as a reference ones since these might change when, for instance, the effect of a thermal -17 envelope is incorporated in the hydrostatic solutions [28] or when additional observational data is included to find the -18 best-fit model [33]. Hence, whichever the fermion mass is, it must describe self-consistently the galaxy size that hosts 50 100 200 the binary pulsar. m [eV] Turning now to our main discussion, we seek for the dependence of the orbital period time-derivative on the or- Figure 3: Secular change of the orbital period as a function of the fermion mass. As before, we have here adopted the values v = bital period. We illustrate this in Fig. 2 for the two partic- −1 100 km s and β = π/2. The curves refer to numerical solutions for ular choices of fermion masses m = 90 eV (top panel) and different binary systems as indicated below: 1.3–1.3 M NS-NS (top m = 200 eV (bottom panel) as indicated for different bi- panel), 1.3–0.2 M NS-WD (middle panel) and 0.25–0.50 M WD- ⊙ ⊙ nary star configurations. As expected for the DM dynam- WD (bottom panel). For this case, we have taken P = 100 days. halos. ical friction effect, the larger the orbital period the larger orbital period time-derivative in agreement with former studies [40, 41]. Notice also that NS-NS binary systems solutions exhibit a simple scaling dependence for a given experience from P & 0.4 days larger dynamical friction star configuration. We can see readily that these solutions and therefore larger orbital decay while for NS-WD and are in agreement with the ones displayed in Fig. 2 when WD-WD binaries, their orbital period time-derivatives are P = 100 days. Interestingly, NS-NS binaries lead to a essentially equal from P & 1 days and smaller than the slightly larger (but appreciable) secular changes of P as ones given by NS-NS binary systems for a fixed fermion can be also appreciated in all the plots. This makes then mass. These features can be observed in Fig. 2. this binary system an ideal target to test our theoretical Another crucial point, aiming to find an ideal scenario predictions. For any fixed orbital period, the larger the to test the DM dynamical friction effect, is the fact that fermion mass the larger the orbital period time-derivative such an effect could be observed, with major astrophysical which implies that this kinematic effect may be larger in expectations, in binaries with large orbital periods which binary systems orbiting inside small DM halos as the satel- is also a challenge of outstanding precision of pulsar ob- lite galaxies of the Milky-Way (e.g. m & 150 eV). Ac- servations. cordingly, this effect is almost reduced by one order of magnitude in Milky-Way like galaxies (e.g. m ∼ 90 eV). 5.3. The effect of the DM fermion mass Finally, we compute the secular change of P with re- 5.4. Upper limit on the fermion mass from Pulsar timimg spect to the mass of the degenerate fermion. To do so, we observations firstly have to set for the fermion masses their associated It is instructive at this point to use some available mea- central densities ρ whose associated configurations may surements of the orbital period time-derivative for long- describe the structure of galactic halos. As mentioned, period binaries in the Milky-Way to quantify more real- once the value of (m , ρ ) is set, it defines the size of the f 0 istically such an effect, as well as to put constraints on DM halo that harbor the binary pulsar. the fermion mass. We evaluate then this effect and show A first systematic approach could be to establish it it in Table (1) for different binaries-type with measured from the mass-radius relation Eqn. (6) but we find this orbital period time-derivatives . For such a computation is a weak correlation despite it describes successfully the we have considered a benchmark m = 1 keV in accor- galactic DM sizes. Instead, we use a much reliable strat- dance with the bounds derived below. As we have dis- egy as was outlined in [32] by using further astrophysical cussed, NS-NS binaries (which correspond to J1740-3052 constraints such as the halo radius and the (nearly con- and J1903+0327 in our case) show slightly larger orbital stant) central surface DM density of galaxies [52]. This period time-derivatives due to DM dynamical friction (col- −2 is Σ = ρ R(m , ρ ) with Σ ≈ 100 M pc . From 0 0 f 0 0 ⊙ umn 7) and can be used to put stringent bounds on the this scaling relation, we find the (m , ρ ) parameter space f 0 consistent with the desire configurations such that we can estimate more reliably the orbital decay of binary stars These binary stars with long orbital periods have depending on the type of their host galaxies. been identified from the ATNF Pulsar Catalogue The result is plotted in Fig. 3 for different binary con- http://www.atnf.csiro.au/research/pulsar/psrcat/ [58]. figurations as before. As we previously advertised, the P b obs ˙ ˙ Name Type m [M ] m [M ] P [days] |P | |P | Ref. p ⊙ c ⊙ b b −9 −13 J1740-3052 NS-NS 1.4 2 231.03 3 × 10 1.08 × 10 [53] −13 −15 J1713+0747 NS-WD 1.31 0.286 67.82 2.0 × 10 7.08 × 10 [54] −12 −16 J0437-4715 NS-WD 1.76 0.254 5.741 3.73 × 10 5.67 × 10 [55] −11 −15 J2019+2425 NS-WD 1.33 0.35 76.51 3 × 10 9.41 × 10 [56] −12 −14 J1903+0327 NS-NS 1.667 1.029 95.17 53 × 10 2.77 × 10 [57] Table 1: Limits in the measured orbital period time-derivatives for several binary systems in the Galaxy (column 6) [58] and the ones predicted by DM dynamical friction of degenerate light fermions (column 7). We have set m = 1 keV and assumed that the central density ρ at the f 0 binary position can be approximated to the DM density near the sun position, namely ρ = 0.355 GeV cm . Here, it is also shown the values of mass binaries (pulsars and their companions) and orbital periods. fermion mass . Moreover, it also depends on which con- capability of observing binary pulsars in extreme environ- figuration of the binary system observed has a compara- ments, with improved timing precision, and detecting pul- ble orbital period time-derivative with the one predicted sars in nearby galaxies [60, 61]. Specifically SKA expects by DM dynamical friction. In order to set a bound for to achieve a precision in the pulse arrival time by a fac- the fermion mass from timing pulsar measurements, we tor of 100 better than current measurements. SKA is also demand that the orbital period time derivative due to dy- expected to detect very long period pulsars in binary sys- namical friction does not excess the value limit of the ob- tems and to discover new systems up to O(10 ) in high served one for each binary system in consideration. Thus, DM densities, including binary systems in nearby galaxies we look for a maximum value of the fermion mass such within 5−10 Mpc, depending on the strength of the pulses that the mentioned condition is satisfied. For instance, by [62, 63]. Thus, the current bound can be improved up to using J1903+0327 NS-NS measurements we can estimate two orders of magnitude. Hence, we expect to reach the the upper bound for the fermions mass technological improvement in the future that permits us to measure with outstanding precision orbital period de- m . 2 keV, (21) cays, under the conditions found, to test our theoretical predictions. while the orbital period time-derivative measured in J1713+0747 In summary, we have studied in this work the effect NS-WD, which is roughly two order of magnitude smaller of DM dynamical friction in the orbital evolution of bi- than the J1903+0327 binary, yields nary pulsars in which DM halos are constituted of self- gravitating degenerate gas of free fermions. Thereby, we m . 1 keV, (22) have accessed to quantify the secular change in the orbital period for this DM model and therefore to devise the best which are, however, above the ones found in Refs. [27, 28, astrophysical scenario to constrain the properties of light 32] by means of velocity dispersion measurements. Nev- fermionic DM by using pulsar timing observations. ertheless, DM dynamical friction can be considerably in- creased, as we showed, in binary pulsars orbiting small host galaxies (and/or) with very long orbital periods 6. Conclusions whereby it would permit us to put more stringent bounds on the fermions mass. It is important however to men- We have assessed the conditions under which DM dy- tion that there has been discovered, up to the best of our namical friction effect of a degenerate and non-interacting knowledge, around 30 extragalactic pulsars inside satellite Fermi gas can be potentially tested with the help of new galaxies of the Milky-Way, but most of them with unmea- generation of pulsar surveys. The main conclusions of this sured orbital periods and orbital period time-derivatives . work are: Despite the current limitations of measuring the or- 1. DM dynamical friction can be considered as a (an bital period decay of extragalactic pulsars and binary pul- additional) kinematic effect in the observed orbital sars in the Galaxy with very long orbital periods, there decay of binary pulsars when high sensitivity of pul- are encouraging prospects for high-precision pulsar timing sar observations is remarkably attained. with the new generation of radio telescopes. In partic- 2. We have showed that this effect may become impor- ular, the Square Kilometre Array (SKA) will reach the tant in future measurements of the orbital period de- cay in binaries with characteristic long periods, i.e. Though a direct comparison with the other NS-WD binaries can for P & 100 days (see Fig. 2). not be done directly since they differ in orbital periods and orbital 3. We have quantified the effects of the DM dynam- period time-derivatives. It remains the possibility of considering self-interactions between ical friction in different binary star configurations. degenerate DM fermions that may increase the local density and lead More specifically, we found that NS-NS binary sys- to a larger orbital period time-derivative. This is however beyond tems experience slightly larger orbital period decays, the scope of this work. −15 i.e. P & 10 for P & 100 days and m = 200 eV. For more details see Ref. [59] and references therein. b b f 8 4. We have computed the DM dynamical friction effect References for some binary stars in the Milky-Way with mea- [1] D. S. Akerib, S. Alsum, H. M. Arau´jo, X. Bai, A. J. sured orbital period time-derivatives (see Table 1). Bailey, J. Balajthy, P. Beltrame, E. P. Bernard, A. Bern- From the J1713+0747 binary, the constraint on the stein, T. P. Biesiadzinski, Results from a Search for Dark Matter in the Complete LUX Exposure, Phys. fermion mass m . 1 keV have been put. Moreover, Rev. Lett.118 (2) (2017) 021303. arXiv:1608.07648, this bound can be considerably improved by using doi:10.1103/PhysRevLett.118.021303. timing pulsar measurements of very long-orbital pe- [2] PandaX-II Collaboration, :, X. Cui, A. Abdukerim, W. Chen, riod binary pulsars. X. Chen, Y. Chen, B. Dong, D. Fang, C. Fu, Dark Matter Results From 54-Ton-Day Exposure of PandaX-II Experiment, 5. We found another promising situation in which bi- arXiv e-prints (2017) arXiv:1708.06917arXiv:1708.06917. nary pulsars orbiting small galaxies, which corre- [3] E. Aprile, J. Aalbers, F. Agostini, M. Alfonsi, L. Althueser, sponds to degenerate fermionic DM halos with m ∼ f F. D. Amaro, M. Anthony, F. Arneodo, L. Baudis, B. Bauer- meister, Dark Matter Search Results from a One Ton-Year Ex- 200 eV, experience larger orbital period decay by dy- posure of XENON1T, Phys. Rev. Lett.121 (11) (2018) 111302. namical friction because of the enhancement of this arXiv:1805.12562, doi:10.1103/PhysRevLett.121.111302. effect with the fermion mass (see Figs. 1-3). Thus, [4] M. Aguilar, et al., Antiproton Flux, Antiproton-to-Proton Flux the fermion mass establishes the galaxy size that har- Ratio, and Properties of Elementary Particle Fluxes in Primary Cosmic Rays Measured with the Alpha Magnetic Spectrometer bor the binary pulsar. Hence, the established bound on the International Space Station, Phys. Rev. Lett. 117 (9) in Eqn. (22) for the fermion mass can also be im- (2016) 091103. doi:10.1103/PhysRevLett.117.091103. proved in such small host galaxies. Interestingly, the [5] A. Albert, B. Anderson, K. Bechtol, A. Drlica-Wagner, SKA’s sensitivity will be sufficient to detect pulsars M. Meyer, M. S´anchez-Conde, L. Strigari, M. Wood, T. M. C. Abbott, F. B. Abdalla, Searching for Dark Matter Anni- in nearby galaxies which makes of our results testable hilation in Recently Discovered Milky Way Satellites with predictions in the near future [60, 61]. Fermi-Lat, ApJ834 (2) (2017) 110. arXiv:1611.03184, 6. It has been shown that there are distinctive theoret- doi:10.3847/1538-4357/834/2/110. [6] M. Aaboud, G. Aad, B. Abbott, J. Abdallah, O. Abdi- ical predictions of the orbital period time-derivative nov, B. Abeloos, R. Aben, O. S. Abouzeid, N. L. Abraham, for certain DM model as the ones found in this work. H. Abramowicz, Search for dark matter in association with a Hence, we are reaching a phenomenological situation √ Higgs boson decaying to b-quarks in pp collisions at { s} = 13 in which we may constrain or even rule out DM mod- TeV with the ATLAS detector, Physics Letters B 765 (2017) 11– els by using pulsar timing measurements. 31. arXiv:1609.04572, doi:10.1016/j.physletb.2016.11.035. [7] CMS Collaboration, Search for associated production of dark matter with a Higgs boson decaying to b b-bar or DM dynamical friction is an appealing effect on the gamma-gamma at sqrt(s) = 13 TeV, arXiv e-prints (2017) binary evolution because it would permit us to put con- arXiv:1703.05236arXiv:1703.05236. straints on the local DM environment the binaries are em- [8] J. D. Bekenstein, Relativistic gravitation theory for the modified bedded due to the high-precision measurements which is a Newtonian dynamics paradigm, Phys. Rev. D70 (2004) 083509. arXiv:astro-ph/0403694, doi:10.1103/PhysRevD.70.083509. characteristic property in such systems. We found as the [9] D. Clowe, M. Bradaˇc, A. H. Gonzalez, M. Markevitch, S. W. main conclusion of this work that NS-NS binary systems Randall, C. Jones, D. Zaritsky, A Direct Empirical Proof with large orbital periods P & 100 days orbiting small DM of the Existence of Dark Matter, ApJ648 (2006) L109–L113. halos, composed of degenerate fermions of m ∼ 200 eV, arXiv:astro-ph/0608407, doi:10.1086/508162. [10] M. Tegmark, M. R. Blanton, M. A. Strauss, F. Hoyle, (which correspond to extragalactic pulsars), are the best D. Schlegel, R. Scoccimarro, M. S. Vogeley, D. H. Wein- astrophysical scenario to test the effect of dynamical fric- berg, I. Zehavi, A. Berlind, The Three-Dimensional Power tion of light fermionic DM once observational data of tim- Spectrum of Galaxies from the Sloan Digital Sky Sur- ing pulsar is available. Interestingly, there are promising vey, ApJ606 (2) (2004) 702–740. arXiv:astro-ph/0310725, doi:10.1086/382125. pulsar surveys that can reach an astonishing sensitivity to [11] Planck Collaboration, N. Aghanim, Y. Akrami, M. Ash- test the theoretical predictions based on DM models. down, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, Planck 2018 re- sults. VI. Cosmological parameters, arXiv e-prints (2018) Acknowledgements arXiv:1807.06209arXiv:1807.06209. [12] V. C. Rubin, W. K. Ford, Jr., Rotation of the Andromeda Neb- We thank anonymous referee for all the insightful ques- ula from a Spectroscopic Survey of Emission Regions, ApJ159 (1970) 379. doi:10.1086/150317. tions and suggestions that help to improve this manuscript. [13] Y. Sofue, Rotation Curve and Mass Distribution in the Galac- We would like also to offer our gratitude to Alexander tic Center - From Black Hole to Entire Galaxy, PASJ65. Gallego for careful and meticulous reading of this paper. arXiv:1307.8241, doi:10.1093/pasj/65.6.118. Special thanks go to Laura Becerra, Nicol´as Bernal and [14] G. Steigman, M. S. Turner, Cosmological con- straints on the properties of weakly interacting mas- Cl´ement Stahl for insightful comments and discussions on sive particles, Nuclear Physics B 253 (1985) 375–386. different subjects of this work. This project was supported doi:10.1016/0550-3213(85)90537-1. by the Postdoctoral Fellowship Program N 2018000101 of [15] L. Hui, J. P. Ostriker, S. Tremaine, E. Witten, Ultralight scalars as cosmological dark matter, Phys. Rev. D95 (2017) 043541. the Vicerrector´ıa de Investigaci´on y Extensio´n-UIS. arXiv:1610.08297, doi:10.1103/PhysRevD.95.043541. [16] S. Dodelson, L. M. Widrow, Sterile neutrinos as dark 9 matter, Physical Review Letters 72 (1994) 17–20. [37] J. Antoniadis, P. C. C. Freire, N. Wex, T. M. Tauris, R. S. arXiv:hep-ph/9303287, doi:10.1103/PhysRevLett.72.17. Lynch, M. H. van Kerkwijk, M. Kramer, C. Bassa, V. S. Dhillon, [17] A. Del Popolo, M. Le Delliou, Small Scale Problems of T. Driebe, J. W. T. Hessels, V. M. Kaspi, V. I. Kondratiev, the ΛCDM Model: A Short Review, Galaxies 5 (2017) 17. N. Langer, T. R. Marsh, M. A. McLaughlin, T. T. Pennucci, arXiv:1606.07790, doi:10.3390/galaxies5010017. S. M. Ransom, I. H. Stairs, J. van Leeuwen, J. P. W. Ver- [18] B. Moore, Evidence against dissipation-less dark matter from biest, D. G. Whelan, A Massive Pulsar in a Compact Rel- observations of galaxy haloes, Nature370 (1994) 629–631. ativistic Binary, Science 340 (2013) 448. arXiv:1304.6875, doi:10.1038/370629a0. doi:10.1126/science.1233232. [19] W. J. G. de Blok, The Core-Cusp Problem, Ad- [38] J. Binney, S. Tremaine, Galactic Dynamics: Second Edition, vances in Astronomy 2010 (2010) 5. arXiv:0910.3538, Princeton University Press, 2008. doi:10.1155/2010/789293. [39] J. D. Bekenstein, R. Zamir, Dynamical friction in binary sys- [20] M. Boylan-Kolchin, J. S. Bullock, M. Kaplinghat, Too big tems, ApJ359 (1990) 427–437. doi:10.1086/169075. to fail? The puzzling darkness of massive Milky Way [40] P. Pani, Binary pulsars as dark-matter probes, Phys. subhaloes, MNRAS415 (2011) L40–L44. arXiv:1103.0007, Rev. D92 (12) (2015) 123530. arXiv:1512.01236, doi:10.1111/j.1745-3933.2011.01074.x. doi:10.1103/PhysRevD.92.123530. [21] S. Garrison-Kimmel, M. Boylan-Kolchin, J. S. Bullock, E. N. [41] L. G. G´omez, J. A. Rueda, Dark matter dynamical friction ver- Kirby, Too big to fail in the Local Group, MNRAS444 (2014) sus gravitational wave emission in the evolution of compact-star 222–236. arXiv:1404.5313, doi:10.1093/mnras/stu1477. binaries, Phys. Rev. D96 (2017) 063001. arXiv:1706.06801, [22] P.-H. Chavanis, M. Lemou, F. M´ehats, Models of dark mat- doi:10.1103/PhysRevD.96.063001. ter halos based on statistical mechanics: The classical King [42] A. Caputo, J. Zavala, D. Blas, Binary pulsars as model, Phys. Rev. D91 (6) (2015) 063531. arXiv:1409.7838, probes of a Galactic dark matter disk, Physics of the doi:10.1103/PhysRevD.91.063531. Dark Universe 19 (2018) 1–11. arXiv:1709.03991, [23] C. R. Argu¨elles, J. A. Rueda, R. Ruffini, Theoretical evidence of doi:10.1016/j.dark.2017.10.005. 50 keV fermionic dark matter from galactic observables, ArXiv [43] S. Tremaine, J. E. Gunn, Dynamical role of light neutral lep- e-printsarXiv:1606.07040. tons in cosmology, Physical Review Letters 42 (1979) 407–410. [24] C. Destri, H. J. de Vega, N. G. Sanchez, Quantum WDM doi:10.1103/PhysRevLett.42.407. fermions and gravitation determine the observed galaxy struc- [44] A. Boyarsky, O. Ruchayskiy, D. Iakubovskyi, A lower tures, Astroparticle Physics 46 (2013) 14–22. arXiv:1301.1864, bound on the mass of dark matter particles, J. doi:10.1016/j.astropartphys.2013.04.004. Cosm. Astropart. Phys.3 (2009) 5. arXiv:0808.3902, [25] H. J. de Vega, P. Salucci, N. G. Sanchez, Observational doi:10.1088/1475-7516/2009/03/005. rotation curves and density profiles versus the Thomas- [45] S. Chandrasekhar, Dynamical Friction. I. General Considera- Fermi galaxy structure theory, MNRAS442 (2014) 2717–2727. tions: the Coefficient of Dynamical Friction., ApJ97 (1943) 255. arXiv:1309.2290, doi:10.1093/mnras/stu972. doi:10.1086/144517. [26] R. Ruffini, C. R. Argu¨elles, J. A. Rueda, On the core-halo distri- [46] J. Antoniadis, Gravitational Radiation from Compact Bi- bution of dark matter in galaxies, MNRAS451 (2015) 622–628. nary Pulsars, in: C. F. Sopuerta (Ed.), Gravitational arXiv:1409.7365, doi:10.1093/mnras/stv1016. Wave Astrophysics, Vol. 40 of Astrophysics and Space [27] V. Domcke, A. Urbano, Dwarf spheroidal galaxies as degenerate Science Proceedings, 2015, p. 1. arXiv:1407.3404, gas of free fermions, J. Cosm. Astropart. Phys.1 (2015) 002. doi:10.1007/978-3-319-10488-1_1. arXiv:1409.3167, doi:10.1088/1475-7516/2015/01/002. [47] J. J. Hermes, M. Kilic, W. R. Brown, D. E. Winget, C. Allende [28] L. Randall, J. Scholtz, J. Unwin, Cores in Dwarf Galax- Prieto, A. Gianninas, A. S. Mukadam, A. Cabrera-Lavers, S. J. ies from Fermi Repulsion, MNRAS467 (2017) 1515–1525. Kenyon, Rapid Orbital Decay in the 12.75-minute Binary White arXiv:1611.04590, doi:10.1093/mnras/stx161. Dwarf J0651+2844, ApJ757 (2012) L21. arXiv:1208.5051, [29] S.-H. Oh, C. Brook, F. Governato, E. Brinks, L. Mayer, doi:10.1088/2041-8205/757/2/L21. W. J. G. de Blok, A. Brooks, F. Walter, The Central Slope [48] E. Poisson, C. M. Will, Gravity, 2014. of Dark Matter Cores in Dwarf Galaxies: Simulations ver- [49] L. D. Landau, E. M. Lifshitz, Statistical physics. Pt.1, Pt.2, sus THINGS, Astron. J.142 (1) (2011) 24. arXiv:1011.2777, 1980. doi:10.1088/0004-6256/142/1/24. [50] C. F. B. Macedo, P. Pani, V. Cardoso, L. C. B. [30] S. Tulin, H.-B. Yu, Dark Matter Self-interactions and Crispino, Into the Lair: Gravitational-wave Signatures Small Scale Structure, Phys. Rept. 730 (2018) 1–57. of Dark Matter, ApJ774 (2013) 48. arXiv:1302.2646, arXiv:1705.02358, doi:10.1016/j.physrep.2017.11.004. doi:10.1088/0004-637X/774/1/48. [31] A. Del Popolo, M. Le Delliou, Small scale problems of the [51] B. J. Prager, S. M. Ransom, P. C. C. Freire, J. W. T. Hes- ΛCDM model: a short review, Galaxies 5 (1) (2017) 17. sels, I. H. Stairs, P. Arras, M. Cadelano, Using long-term mil- arXiv:1606.07790, doi:10.3390/galaxies5010017. lisecond pulsar timing to obtain physical characteristics of the [32] C. Di Paolo, F. Nesti, F. L. Villante, Phase-space mass bulge globular cluster Terzan 5, Astrophys. J. 845 (2) (2017) bound for fermionic dark matter from dwarf spheroidal 148. arXiv:1612.04395, doi:10.3847/1538-4357/aa7ed7. galaxies, MNRAS475 (2018) 5385–5397. arXiv:1704.06644, [52] F. Donato, G. Gentile, P. Salucci, C. Frigerio Martins, doi:10.1093/mnras/sty091. M. I. Wilkinson, G. Gilmore, E. K. Grebel, A. Koch, [33] J. Barranco, A. Bernal, D. Delepine, Constraining ultra light R. Wyse, A constant dark matter halo surface density in fermionic dark matter with galactic observations, arXiv e- galaxies, MNRAS397 (2009) 1169–1176. arXiv:0904.4054, printsarXiv:1811.11125. doi:10.1111/j.1365-2966.2009.15004.x. [34] D. R. Lorimer, Binary and Millisecond Pulsars, Living [53] E. C. Madsen, I. H. Stairs, M. Kramer, F. Camilo, G. B. Hobbs, Reviews in Relativity 11 (2008) 8. arXiv:0811.0762, G. H. Janssen, A. G. Lyne, R. N. Manchester, A. Possenti, doi:10.12942/lrr-2008-8. B. W. Stappers, Timing the main-sequence-star binary pulsar [35] R. A. Hulse, J. H. Taylor, Discovery of a pulsar in a binary J1740-3052, MNRAS425 (2012) 2378–2385. arXiv:1207.2202, system, ApJ195 (1975) L51–L53. doi:10.1086/181708. doi:10.1111/j.1365-2966.2012.21691.x. [36] M. Kramer, I. H. Stairs, R. N. Manchester, M. A. McLaugh- [54] W. W. Zhu, I. H. Stairs, P. B. Demorest, D. J. Nice, J. A. lin, A. G. Lyne, R. D. Ferdman, M. Burgay, D. R. Lorimer, Ellis, S. M. Ransom, Z. Arzoumanian, K. Crowter, T. Dolch, A. Possenti, N. D’Amico, J. M. Sarkissian, G. B. Hobbs, J. E. R. D. Ferdman, E. Fonseca, M. E. Gonzalez, G. Jones, M. L. Reynolds, P. C. C. Freire, F. Camilo, Tests of General Relativ- Jones, M. T. Lam, L. Levin, M. A. McLaughlin, T. Pennucci, ity from Timing the Double Pulsar, Science 314 (2006) 97–102. K. Stovall, J. Swiggum, Testing Theories of Gravitation Using arXiv:astro-ph/0609417, doi:10.1126/science.1132305. 21-Year Timing of Pulsar Binary J1713+0747, ApJ809 (2015) 10 41. arXiv:1504.00662, doi:10.1088/0004-637X/809/1/41. [55] J. P. W. Verbiest, M. Bailes, W. van Straten, G. B. Hobbs, R. T. Edwards, R. N. Manchester, N. D. R. Bhat, J. M. Sarkissian, B. A. Jacoby, S. R. Kulkarni, Precision Timing of PSR J0437-4715: An Accurate Pulsar Distance, a High Pul- sar Mass, and a Limit on the Variation of Newton’s Gravi- tational Constant, ApJ679 (2008) 675–680. arXiv:0801.2589, doi:10.1086/529576. [56] D. J. Nice, E. M. Splaver, I. H. Stairs, On the Mass and Incli- nation of the PSR J2019+2425 Binary System, ApJ549 (2001) 516–521. arXiv:astro-ph/0010489, doi:10.1086/319079. [57] P. C. C. Freire, C. G. Bassa, N. Wex, I. H. Stairs, D. J. Cham- pion, S. M. Ransom, P. Lazarus, V. M. Kaspi, J. W. T. Hes- sels, M. Kramer, On the nature and evolution of the unique binary pulsar J1903+0327, MNRAS412 (4) (2011) 2763–2780. arXiv:1011.5809, doi:10.1111/j.1365-2966.2010.18109.x. [58] R. N. Manchester, G. B. Hobbs, A. Teoh, M. Hobbs, The Australia Telescope National Facility Pulsar Catalogue, Astron. J.129 (2005) 1993–2006. arXiv:astro-ph/0412641, doi:10.1086/428488. [59] R. N. Manchester, G. Fan, A. G. Lyne, V. M. Kaspi, F. Craw- ford, Discovery of 14 Radio Pulsars in a Survey of the Magellanic Clouds, ApJ649 (2006) 235–242. arXiv:astro-ph/0604421, doi:10.1086/505461. [60] E. Keane, B. Bhattacharyya, M. Kramer, B. Stappers, E. F. Keane, B. Bhattacharyya, M. Kramer, B. W. Stappers, S. D. Bates, M. Burgay, S. Chatterjee, D. J. Champion, R. P. Eatough, J. W. T. Hessels, G. Janssen, K. J. Lee, J. van Leeuwen, J. Margueron, M. Oertel, A. Possenti, S. Ransom, G. Theureau, P. Torne, A Cosmic Census of Radio Pulsars with the SKA, in: Advancing Astrophysics with the Square Kilome- tre Array (AASKA14), 2015, p. 40. arXiv:1501.00056. [61] B. W. Stappers, E. F. Keane, M. Kramer, A. Possenti, I. H. Stairs, The prospects of pulsar timing with new-generation ra- dio telescopes and the Square Kilometre Array, Philosophical Transactions of the Royal Society of London Series A 376 (2018) 20170293. doi:10.1098/rsta.2017.0293. [62] R. Smits, M. Kramer, B. Stappers, D. R. Lorimer, J. Cordes, A. Faulkner, Pulsar searches and timing with the square kilometre array, Astron. Astrophys. 493 (2009) 1161–1170. arXiv:0811.0211, doi:10.1051/0004-6361:200810383. [63] K. Liu, J. P. W. Verbiest, M. Kramer, B. W. Stappers, W. van Straten, J. M. Cordes, Prospects for high-precision pulsar timing, MNRAS417 (4) (2011) 2916–2926. arXiv:1107.3086, doi:10.1111/j.1365-2966.2011.19452.x.
Journal
Astrophysics – arXiv (Cornell University)
Published: Jun 25, 2019